Land change is a cause and consequence of global environmental change1,2. Changes in land use and land cover considerably alter the Earth’s energy balance and biogeochemical cycles, which contributes to climate change and—in turn—affects land surface properties and the provision of ecosystem services1,2,3,4. However, quantification of global land change is lacking. Here we analyse 35 years’ worth of satellite data and provide a comprehensive record of global land-change dynamics during the period 1982–2016. We show that—contrary to the prevailing view that forest area has declined globally5—tree cover has increased by 2.24 million km2 (+7.1% relative to the 1982 level). This overall net gain is the result of a net loss in the tropics being outweighed by a net gain in the extratropics. Global bare ground cover has decreased by 1.16 million km2 (−3.1%), most notably in agricultural regions in Asia. Of all land changes, 60% are associated with direct human activities and 40% with indirect drivers such as climate change. Land-use change exhibits regional dominance, including tropical deforestation and agricultural expansion, temperate reforestation or afforestation, cropland intensification and urbanization. Consistently across all climate domains, montane systems have gained tree cover and many arid and semi-arid ecosystems have lost vegetation cover. The mapped land changes and the driver attributions reflect a human-dominated Earth system. The dataset we developed may be used to improve the modelling of land-use changes, biogeochemical cycles and vegetation–climate interactions to advance our understanding of global environmental change1,2,3,4,6.
This is a preview of subscription content
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
Turner, B. L. II, Lambin, E. F. & Reenberg, A. The emergence of land change science for global environmental change and sustainability. Proc. Natl Acad. Sci. USA 104, 20666–20671 (2007).
Foley, J. A. et al. Global consequences of land use. Science 309, 570–574 (2005).
Le Quéré, C. et al. Global carbon budget 2016. Earth Syst. Sci. Data 8, 605–649 (2016).
Alkama, R. & Cescatti, A. Biophysical climate impacts of recent changes in global forest cover. Science 351, 600–604 (2016).
FAO. Global Forest Resources Assessment 2015 (UN Food and Agriculture Organization, Rome, 2015).
Bonan, G. B. & Doney, S. C. Climate, ecosystems, and planetary futures: The challenge to predict life in Earth system models. Science 359, eaam8328 (2018).
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).
Feng, M. et al. Earth science data records of global forest cover and change: assessment of accuracy in 1990, 2000, and 2005 epochs. Remote Sens. Environ. 184, 73–85 (2016).
DeFries, R. S. et al. Mapping the land surface for global atmosphere–biosphere models: toward continuous distributions of vegetation’s functional properties. J. Geophys. Res. 100, 20867–20882 (1995).
Gibbs, H. K. et al. Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proc. Natl Acad. Sci. USA 107, 16732–16737 (2010).
Margono, B. A., Potapov, P. V., Turubanova, S., Stolle, F. & Hansen, M. C. Primary forest cover loss in Indonesia over 2000–2012. Nat. Clim. Change 4, 730–735 (2014).
Ordway, E. M., Asner, G. P. & Lambin, E. F. Deforestation risk due to commodity crop expansion in sub-Saharan Africa. Environ. Res. Lett. 12, 044015 (2017).
Hicke, J. A. et al. Postfire response of North American boreal forest net primary productivity analyzed with satellite observations. Glob. Change Biol. 9, 1145–1157 (2003).
van Mantgem, P. J. et al. Widespread increase of tree mortality rates in the western United States. Science 323, 521–524 (2009).
McManus, K. M. et al. Satellite-based evidence for shrub and graminoid tundra expansion in northern Quebec from 1986 to 2010. Glob. Change Biol. 18, 2313–2323 (2012).
Mitchard, E. T. & Flintrop, C. M. Woody encroachment and forest degradation in sub-Saharan Africa’s woodlands and savannas 1982–2006. Phil. Trans. R. Soc. Lond. B 368, 20120406 (2013).
Brandt, M. et al. Human population growth offsets climate-driven increase in woody vegetation in sub-Saharan Africa. Nat. Ecol. Evol. 1, 0081 (2017).
Harsch, M. A., Hulme, P. E., McGlone, M. S. & Duncan, R. P. Are treelines advancing? A global meta-analysis of treeline response to climate warming. Ecol. Lett. 12, 1040–1049 (2009).
Potapov, P. V. et al. Eastern Europe’s forest cover dynamics from 1985 to 2012 quantified from the full Landsat archive. Remote Sens. Environ. 159, 28–43 (2015).
Piao, S. et al. Detection and attribution of vegetation greening trend in China over the last 30 years. Glob. Change Biol. 21, 1601–1609 (2015).
Birdsey, R., Pregitzer, K. & Lucier, A. Forest carbon management in the United States: 1600–2100. J. Environ. Qual. 35, 1461–1469 (2006).
Donohue, R. J., McVicar, T. R. & Roderick, M. L. Climate-related trends in Australian vegetation cover as inferred from satellite observations, 1981–2006. Glob. Change Biol. 15, 1025–1039 (2009).
Liu, Y. Y. et al. Changing climate and overgrazing are decimating Mongolian steppes. PLoS ONE 8, e57599 (2013).
Herrick, J. E. et al. National ecosystem assessments supported by scientific and local knowledge. Front. Ecol. Environ. 8, 403–408 (2010).
Evenson, R. E. & Gollin, D. Assessing the impact of the green revolution, 1960 to 2000. Science 300, 758–762 (2003).
Ying, Q. et al. Global bare ground gain from 2000 to 2012 using Landsat imagery. Remote Sens. Environ. 194, 161–176 (2017).
Myneni, R. B., Keeling, C. D., Tucker, C. J., Asrar, G. & Nemani, R. R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 386, 698–702 (1997).
Forkel, M. et al. Codominant water control on global interannual variability and trends in land surface phenology and greenness. Glob. Change Biol. 21, 3414–3435 (2015).
Zhu, Z. et al. Greening of the Earth and its drivers. Nat. Clim. Chang e 6, 791–795 (2016).
FAO. Global Ecological Zoning for the Global Forest Resources Assessment 2000 http://www.fao.org/docrep/006/ad652e/ad652e00.htm (UN Food and Agriculture Organization, Rome, 2001).
Gessner, U., Machwitz, M., Conrad, C. & Dech, S. Estimating the fractional cover of growth forms and bare surface in savannas. A multi-resolution approach based on regression tree ensembles. Remote Sens. Environ. 129, 90–102 (2013).
This study was funded by the NASA Making Earth System Data Records for Use in Research Environments (MEaSUREs) Program (NNX13AJ35A), Gordon and Betty Moore Foundation (5131), Norwegian Climate and Forests Initiative through the World Resources Institute’s Global Forest Watch project, DOB Ecology through the World Resources Institute’s Global Restoration Initiative, the NASA Land-Cover and Land-Use Change (LCLUC) Program (NNX15AK65G), and the NASA Carbon Monitoring Systems Program (NNX13AP48G). We thank T. Loveland, B. Pengra and P. Olofsson for making their tree cover validation data available, C. Dimiceli for assistance with vegetation continuous field development and Z. Song for assistance with AVHRR calibration.
Nature thanks M. Forkel, L. Zhou and the other anonymous reviewer(s) for their contribution to the peer review of this work.
The authors declare no competing interests.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Satellite-derived, long-term (1982–2016) changes in land cover show strong coupling and symmetry in change detection.
a, Global map of co-located ∆TC and ∆SV. Pixels showing a statistically significant trend (n = 35 years, two-sided Mann–Kendall test, P < 0.05) in both TC and SV are depicted on the map. b, Global map of co-located ∆TC and ∆BG. c, Global map of co-located ∆SV and ∆BG. d, From left to right, intensity plot of change area for ΔTC versus ΔSV, ΔTC versus ΔBG and ΔSV versus ΔBG, corresponding to a, b and c, respectively. To create these intensity plots, paired per cent change layers (Fig. 1b) are used to construct a 2D histogram with bin size of 1% for both axes. Then, the total change area in each bin is calculated and plotted.
Extended Data Fig. 2 Long-term (1982–2016) gross land-change dynamics vary considerably between biomes.
a–p, Gross land-change dynamics per biome. Mountain systems (c, f, i, n) all exhibit larger area of TC gain than TC loss, larger area of SV loss than SV gain and larger area of BG loss than BG gain. q, Geographical distribution of all biomes, from a previous publication30, reproduced with permission. See Fig. 3 for other biomes and Extended Data Table 2 for change area estimates.
Extended Data Fig. 3 Attributing direct human impact versus indirect drivers to detected changes in land cover.
Indirect drivers include both natural drivers and human-induced climate change. a, Spatial distribution of the probability sample used for the attribution estimates (n = 1,500). b, Direct human impact (DHI) of each sample unit interpreted using a time-series of high-resolution images in Google Earth. c, Estimated DHI as a per cent of all change area at the global scale. Global average is calculated by weighting the human impact of each type by each respective global total area provided in Extended Data Table 1. The standard error (SE) for the estimated per cent of DHI is provided in the parentheses. d, e, Estimated DHI at the continental and biome scales. See Extended Data Fig. 4 for some representative sample examples.
Screenshots are taken from Google Earth. Each panel is 0.05° × 0.05° in size, corresponding to one AVHRR pixel. a, Deforestation for industrial agriculture expansion in Mato Grosso, Brazil (11.275° S, 52.125° W). b, Expanding shifting agriculture in northern Zambia (11.625° S, 28.625° E). c, Intensification of small-holder agriculture in Punjab, Pakistan (30.025° N, 71.675° E). d, Short vegetation gain in low-intensity agricultural lands in northern Nigeria (12.825° N, 7.825° E). e, Short vegetation increase due to effective fire suppression in pasture lands in Omaheke, Namibia31 (22.175° S, 18.925° E). f, Managed pasture lands in western Kazakhstan (49.475° N, 47.725° E). g, Forestry in southern Finland (61.075° N, 24.475° E). h, Urbanization in Shanghai, China (30.925° N, 121.175° E). i, Oil extraction in New Mexico, USA (32.875° N, 104.275° W). j, Herbaceous vegetation increase owing to glacial retreat in Chuy, Kyrgyzstan (42.575° N, 74.775° E). k, Bare ground cover variation along Mar Chiquita lake shore in Cordoba, Argentina (30.675° S, 63.025° W). l, Forest fires in Saskatchewan, Canada (55.225° N, 102.225° W). m, Tree cover increase in unpopulated savannahs in Western Equatoria, South Sudan16,17 (6.575° N, 27.725° E). n, Climate-change-driven woody encroachment in Quebec, Canada15 (59.475° N, 73.225° W). Examples a–i show various types of land use, whereas examples j–n do not show visible signs of human activity. Map data: Google, DigitalGlobe, CNES/Airbus, Landsat/Copernicus.
a, Trends in TC cover. b, Trends in SV cover. c, Trends in BG cover. The following steps were taken for each cover type using TC as the example. The TC gain layer (Fig. 1b) was overlaid on the annual TC% stack to compute annual global TC area within the gain mask (solid dark blue lines); the TC loss layer (Fig. 1b) was overlaid on the annual TC% stack to compute annual global TC area within the loss mask (solid dark red lines). Gross gain estimates from 1986 to 2016 are marked by blue arrows and dashed lines; gross loss estimates from 1986 to 2016 are marked by red arrows and dashed lines. See Extended Data Table 1 for exact gross change estimates.
a, Spatial distribution of annual mean root-mean-square-deviation (RMSD) of TC between 1982 and 2016. b, Spatial distribution of annual mean RMSD of BG between 1982 and 2016. c, Spatial distribution of ∆TC uncertainty. d, Spatial distribution of ∆BG uncertainty. e, Normalized frequency distribution of ∆TC uncertainty. f, Normalized frequency distribution of ∆BG uncertainty. TC, BG and associated RMSD values are outputs of regression tree models. Uncertainty is represented by the ratio of long-term TC (or BG) change estimates to respective RMSD estimates. Positive values of the ratio metric represent the uncertainties of gains and negative values represent the uncertainties of losses. A greater absolute value indicates lower uncertainty, and vice versa. Area under the frequency distribution equals 1. The frequency distributions suggest that tree cover gain exceeds tree cover loss and bare ground loss exceeds bare ground gain for any threshold level (for example, dashed lines), hence the observed trends (a net gain in tree cover and a net loss in bare ground cover over the study period) are valid.
This file contains the Supplementary Methods and Supplementary Figures 1-2.
Estimates of 1982 land cover area, 1982-2016 annual net land-cover change and 1982-2016 gross land-cover change at national scale. Consistent with Extended Data Tables 1 and 2, annual net land-cover change (slope) and 1982 land-cover area were estimated using Theil-Sen regression of time series of annual land-cover area per country. Lower and upper slopes represent the 90% confidence interval. Reported P value is for the two-sided Mann-Kendall test for trend, with P < 0.05 used to define statistically significant, and a sample size of n = 35 years. Gross land-cover change was estimated based on per-pixel non-parametric trend analysis. Per-pixel loss and gain were summed to derive gross loss and gain at the aggregated scale. Country boundary shapefile was downloaded from Global Administrative Areas (GADM) (http://www.gadm.org). Only countries with land area greater than 105 km2 (based on 0.05º grid) are included.
Annual phenological metrics that were used to generate the AVHRR vegetation continuous fields (VCF) land-cover product.
About this article
Cite this article
Song, XP., Hansen, M.C., Stehman, S.V. et al. Global land change from 1982 to 2016. Nature 560, 639–643 (2018). https://doi.org/10.1038/s41586-018-0411-9
- Land Change Map
- Climate Domain
- Bare Ground (BG)
- Indirect Drivers
- Vegetation Continuous Fields
Malaria Journal (2021)
First comprehensive quantification of annual land use/cover from 1990 to 2020 across mainland Vietnam
Scientific Reports (2021)
Nature Sustainability (2021)
Bat responses to changes in forest composition and prey abundance depend on landscape matrix and stand structure
Scientific Reports (2021)