Land-use models and integrated assessment models provide scenarios of land-use and land-cover (LULC) changes following pathways or storylines related to
different socioeconomic and environmental developments. The large diversity of available scenario projections leads to a recognizable variability
in impacts on land ecosystems and the levels of services provided. We evaluated 16 projections of future LULC until 2040 that reflected different
assumptions regarding socioeconomic demands and modeling protocols. By using these LULC projections in a state-of-the-art dynamic global vegetation model,
we simulated their effect on selected ecosystem service indicators related to ecosystem productivity and carbon sequestration potential,
agricultural production and the water cycle. We found that although a common trend for agricultural expansion exists across the scenarios, where and
how particular LULC changes are realized differs widely across models and scenarios. They are linked to model-specific considerations of some
demands over others and their respective translation into LULC changes and also reflect the simplified or missing representation of processes
related to land dynamics or other influencing factors (e.g., trade, climate change). As a result, some scenarios show questionable and possibly
unrealistic features in their LULC allocations, including highly regionalized LULC changes with rates of conversion that are contrary to or
exceed rates observed in the past. Across the diverging LULC projections, we identified positive global trends of net primary productivity
(
The recently presented Intergovernmental Panel on Climate Change (IPCC) Special Report on Climate Change and Land (IPCC, 2019) highlighted unprecedented rates of land and freshwater use and
biodiversity loss and underpinned existing socioeconomic, ecological and climatic challenges such as increasing per capita food consumption, land degradation and an accumulation of climate extreme events. The Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) global assessment report published earlier in 2019 (IPBES, 2019) also reported
deteriorating levels of most ecosystem services (ESs) and natural capital due to past and current human activities. The cumulative contribution of
land-use and land-cover (LULC) change to global
Future LULC and changes therein are modeled based on initial conditions of land use along with LULC history and different assumptions about possible socioeconomic and environmental developments regarding population growth, international cooperation, consumption preferences or technological developments. All of these are represented differently in land-use models (LUMs) or integrated assessment models (IAMs, e.g., DeFries et al., 2004; Meiyappan et al., 2014; van Vliet et al., 2016). However, these models play a central role in assessing possible climate change mitigation and adaptation or conservation strategies in terms of total land demand, investment and maintenance costs, and direct and indirect socioeconomic and ecological effects (e.g., Humpenöder et al., 2014; Popp et al., 2014; Reilly et al., 2012).
In total, the diversity of models, initial model conditions, socioeconomic pathways, climate mitigation targets, processes and process feedbacks considered in the LULC modeling procedure leads to a large number of diverging land-use projections. This reflects not only the fact that the future is unknown but also a large uncertainty introduced by the model structure itself (e.g., Alexander et al., 2017; Prestele et al., 2016; Schmitz et al., 2014; Stehfest et al., 2019; van Vliet et al., 2016). By evaluating a large set of LULC projections (75 and 43, respectively), Alexander et al. (2017) and Prestele et al. (2016) attributed a significant share of the uncertainty in global and regional LULC projections to the model initial conditions, resulting in part from different LULC definitions (especially for pastures; see also, e.g., Verburg et al., 2011), followed by the model structure, scenario storyline and other factors. Alexander et al. (2017) identified the differences in projected global LULC associated with the modeling approach to be at least as great as the differences due to scenario variations. In a regional-level analysis, Prestele et al. (2016) found the highest uncertainty in land-use projections generally at the boundaries of boreal and tropical forests. LULC projections have also been evaluated in a number of model intercomparison studies, in which models simulated the same scenario storylines based on harmonized drivers, in order to focus on the uncertainty in LULC changes resulting from structural differences between the models (e.g., Von Lampe et al., 2014; Popp et al., 2017; Schmitz et al., 2014; Stehfest et al., 2019).
The large uncertainties in LULC projections affect the confidence in projected changes in ecosystem functioning globally, which critically underpins the supply of future ESs available to human societies. In the same ways as the effects of climate change, the uncertainties arising from different LULC projections need to be identified and understood to adapt ecosystems in a sustainable way and possibly counteract critical regional trends. Studies have focused on the vulnerability of ecosystems and their services to changes in climate (e.g., Ahlström et al., 2012; Huntingford et al., 2011; Ostberg et al., 2013; Scholze et al., 2006), land use on a global or regional scale (e.g., Arora and Boer, 2010; Foley et al., 2005; Jantz et al., 2015; Krause et al., 2017; Lawler et al., 2014; Sterling et al., 2013), and a combination of climate and land-use effects (e.g., Dunford et al., 2015; Kim et al., 2018; Krause et al., 2019; Rabin et al., 2020). Moreover, uncertainties arising from different ES quantification methods were estimated (e.g., Schulp et al., 2014). These studies have already begun to document that diverging LULC projections are as important as diverging climate change scenarios with respect to the degree of impact on ecosystems.
Thus, we expand these previous studies here by bringing together a larger number of LULC scenarios and by critically examining the resulting variability in diverging LUM projections based on recent historical observations. We intend to also highlight how different LULC patterns impact ecosystems and related ES indicators. This supports the interpretation of conclusions derived from LUMs and IAMs towards policy decisions regarding factors such as intensification, conservation or climate change mitigation options. A broad range of future LULC projections from different LUMs and different socioeconomic assumptions is important, given the unknown future. Nevertheless, critically assessing the spatial pattern and rates of change can support their interpretation in terms of plausibility.
Our basis was 16 projections of future land use from five LUMs or IAMs with different modeling protocols and socioeconomic pathways. Their scenario storylines span a wide range of world views and policies, with some implemented to achieve a certain climate mitigation or conservation target while others focus only on basic demands for agricultural commodities, built-up area and so on. Models and scenarios were assessed based on their underlying demands, modeling protocols (e.g., assumptions involved, allocation strategies) and the projected spatially explicit land-use futures that they describe. We then used the 16 land-use projections as input for simulations with a state-of-the-art dynamic global vegetation model (DGVM) to analyze their effects on ecosystem functionality and six selected ES indicators linked to the productivity and carbon (C) sequestration potential of ecosystems (net primary productivity, vegetation C, soil and litter C), agricultural production (crop production) and the water cycle (evapotranspiration and annual runoff). We focused on changes until 2040, i.e., the near to medium future.
Overview of land-use models and scenarios used in this study. Main references are given for each model and the scenarios applied in this study; the reader is also referred to the references therein for further details.
shifting cultivation rates;
deforestation for agricultural land included in wood harvest statistics;
priority for land conversion (agricultural land taken from primary or secondary).
crop type and rotations;
shifting cultivation rates;
historical wood harvest statistics and future forest transitions;
biomass density and recovery rates;
priority for land conversion.
Continued.
FAOref: demands with respect to tons of crop production, size of livestock and population (as built-up area) as expected by the FAO until 2040. Demands projected to 24 world regions by the integrated assessment model IMAGE.
CStor: in addition to FAOref there is extra demand for C storage implemented with a no-carbon-loss approach. This results in the favorable consideration of land systems with high C storage capacity such as dense forests in the allocation procedure.
Bdiv: in addition to FAOref there is extra demand for protected areas. Protected areas are used that follow national biodiversity targets according to Aichi Target 11. Forest and grassland proportions were selected to account for protected areas in the CLUMondo land systems “dense forest”, “mosaic grassland and forest”, “mosaic grassland and bare” and “natural grassland”. Favorable consideration of land systems with high amount of protected areas in the allocation procedure. Land systems in land areas assigned to the International Union for Conservation of Nature (IUCN) categories I to IV at the beginning of the simulation procedure were maintained throughout the scenario.
Base: land-use change is driven by increased food demand and population growth according to SSP2.
ADAFF: demands are the same as in the Base scenario with an additional CDR target of 130
BECCS: demands are the same as in the Base scenario with an additional CDR target of 130
Base: land-use change is driven by increased food demand and population growth according to SSP2.
ADAFF: demands are the same as in the Base scenario with an additional CDR target of 130
BECCS: demands are the same as in the Base scenario with an additional CDR target of 130
26Be: scenario following RCP2.6 where the cultivation of bioenergy crops with CCS sequestering carbon significantly balances C so that global warming saturates below 2
45Aff: scenario following RCP4.5 assuming that global greenhouse gas emission prices support climate change mitigation by massive afforestation programs around the world with the aim of preserving large C stocks in forests. Vast expansion of forested area slightly decreases the area of arable land, but assumed technological improvements and more efficient global trading systems maintain agricultural production to satisfy demands. Simulated by the dynamic recursive economic model GCAM.
60Stab: scenario following RCP6.0 assuming stabilization with a medium population growth throughout the 21st century; climate change mitigation efforts only become effective in the last third of the century. Intensification and expansion of cropland area cover the increased food demand. Simulated by AIM and other models.
85Pop: scenario following RCP8.5 assuming rapid population growth to 12 billion people in 2100 and a global temperature rise of more than 4
SSP1-26: scenario combining SSP1 and RCP2.6 of an environmentally friendly world with low population growth, high urbanization, relatively low demand for animal products and high agricultural productivity. Land use and land-use change is relatively low, and policies targeting the reduction of atmospheric greenhouse gases (including BECCS and afforestation) are implemented so that an additional forcing is limited to 2.6
SSP3-70: scenario combining SSP3 and RCP7.0 that is focused on regional development. High population growth in developing countries as well as slow economic development and per capita food consumption lead to the high expansion of agricultural areas rather than agricultural land intensification under a relatively high level of climate change with a radiative forcing of 7.0
SSP5-85: scenario combining SSP5 and RCP8.5 of a world characterized by strong economic growth that is based on fossil fuels, with low population growth, high urbanization, and high food demand per capita with high agricultural productivity due to technological progress. Strong expansion of global cropland and the level of climate change is high with an assumed radiative forcing of 8.5
We used a total set of 16 land-use scenarios originating from five different LUMs or IAMs. The models differ in their underlying demands, modeling protocols and technical aspects (e.g., number of represented land-use classes, time horizons), which are summarized in this section and in Table 1. Although some of the models considered here are IAMs including a land-use component, we refer to all models in this study as global LUMs for the remainder of this paper because their projected LULC change is the target of this analysis. This includes the two versions of the Land-Use Harmonization (LUH) project prominently applied in many studies of the last and the upcoming IPCC reports, although these land-use products are based on the outputs of several LUMs and IAMs.
The CLUMondo model (van Asselen and Verburg, 2013) applies 30 land system types to model LULC changes. Land systems define typical combinations of shares of cropland, grassland, bare land and built-up land as well as a specific management intensity (e.g., extensive cropland with few livestock). Land systems are dynamically allocated based on local suitability, spatial restrictions and the competition between land systems to fulfill demands that were created exogenously by the IMAGE model on the level of world regions. Trade between world regions is excluded in CLUMondo. Eitelberg et al. (2016) designed three CLUMondo scenarios: a reference scenario following the development of basic demands as expected by the Food and Agriculture Organization (FAO) and two scenarios that additionally included a policy target of reducing deforestation and greenhouse gas emissions with a higher ecosystem carbon storage and international policy targets for the prevention of biodiversity loss.
The Integrated Model to Assess the Global Environment (IMAGE) is an IAM framework including sub-models representing the energy system, agricultural economy, land use and the climate system (Stehfest et al., 2014). LULC allocation is done following an assessment and ranking of land's suitability to fulfill demands. From IMAGE, a LULC baseline projection following increased food demand and population growth according to the Shared Socioeconomic Pathway (SSP) 2 (O'Neill et al., 2014) is available as well as two additional scenarios involving land-based climate change mitigation, either via the conservation and expansion of global forest area (ADAFF) or bioenergy crop cultivation and subsequent carbon capture and storage (BECCS) (Krause et al., 2017).
MAgPIE is a global land-use model of the agricultural sector (Lotze-Campen et al., 2008; Popp et al., 2014). It optimizes spatially explicit land-use patterns in a recursive dynamic way to satisfy given commodity demands at minimal production costs while meeting biophysical and socioeconomic constraints. Options to fulfill increasing demands are intensification (yield-increasing technologies), cropland and pasture expansion, and international trade. Future land-use projections of the MAgPIE model follow the same storylines as described for IMAGE (Krause et al., 2017).
The Hurtt et al. (2011) modeling approach (LUH1) combines a historic land-use reconstruction with national statistics of historical wood harvest and assumptions regarding shifting cultivation in some tropical regions, and it harmonizes these data with a set of four future LULC scenarios. Each scenario was produced by a different IAM with individual demands and strategies for allocating LULC in response to the demands. The four scenarios follow very different socioeconomic storylines that are combined with the emissions and climate change assumptions of the Representative Concentration Pathways (RCPs). LUH1 scenarios are not tied to the SSPs, as those were only introduced in 2014. These scenarios have frequently been used in the modeling community, especially for the work in the IPCC AR5 (e.g., O'Neill et al., 2014; van Vuuren et al., 2014). Here, we used the version of LUH1 which had the historical dataset extended until 2014 (Le Quéré et al., 2015) as well as future trajectories following IAM implementations of RCPs 2.6, 4.5, 6.0 and 8.5.
The Land-Use Harmonization v2 (LUH2; v2.1, Hurtt et al., 2020) has been developed for the Coupled Model Intercomparison Project Phase 6 (CMIP6; Eyring et al., 2016). It follows a similar methodology as in LUH1, although on a higher spatial resolution and over a longer time domain, using updated historical land-use reconstructions along with updated models of past and future land transitions and management (e.g., wood harvest, crop rotations and shifting cultivation) and extending the number of scenarios by combining RCPs with SSPs. Similar to LUH1, future LULC transitions in LUH2 are based on land-use projections from different IAMs that each follow their own strategy for allocating LULC in response to demands. Of the eight scenarios that have been harmonized in LUH2 with historical data, three scenarios were selected (SSP1-26, SSP3-70 and SSP5-85) to span the range from low to high radiative forcing as in the LUH1 scenarios in combination with diverging land-use trends according to the SSPs.
The process-based dynamic global vegetation model LPJ-GUESS simulates vegetation dynamics in response to climate, atmospheric
LPJ-GUESS was run at a 0.5
We modeled the ice-free land surface and included only those grid cells in our simulations for which all LUMs provided data. Table S1 in the Supplement provides the detailed simulation setup with all forcing data for the LPJ-GUESS simulations. The differences in the modeling protocol of CLUMondo/LUH1/LUH2 and IMAGE/MAgPIE simulations will affect the base level of ES indicators in 2000–2004 to some degree, although the impacts of slightly diverging historical model periods, spin-up and historical climate would have mostly disappeared by the beginning of the 21st century (baseline period). Larger effects arise from the intrinsic differences in the individual LUMs (see also Alexander et al., 2017). In principle, differences in the baseline land-cover maps could spill over to the simulated degree of change in the future scenarios. For instance, the presence or absence of natural vegetation in the baseline maps might translate into variable degrees of future (semi)natural vegetation regrowth. However, this would only be an important consideration when comparing similar scenarios (and their underpinning storylines related to factors such as sustainability). The alternative approach of harmonizing the different projections to the same starting point of land cover would artificially mask some of the simulated differences in ES indicators which would be contrary to our objectives. Therefore, LUM data were taken as they are, with each LUM scenario providing a seamless transition from historical to future, which is needed to simulate vegetation and carbon cycle responses.
Some of the variables assessed in LPJ-GUESS would also be computed in the models that deliver the LULC change scenarios – most notably crop yields and some carbon-cycle- or water-cycle-related variables. The spatial patterns of these would differ in the LUMs and LPJ-GUESS. However, this does not affect our analysis: here we take the LULC change projections in a unidirectional approach to assess impacts on ecosystem processes; we do not compare similar ecosystem output variables across different model types.
LULC fractions were taken as net annual transitions from the LUMs and aggregated to the three land-use types – cropland, pasture and natural land – used
by LPJ-GUESS (see Table S2 in the Supplement) and to the spatial resolution of 0.5
Changes in ecosystem function and services were assessed using a suite of regulating and provisioning ES indicators: net primary productivity (NPP,
foremost an indicator for ecosystem productivity and C sequestration related to global climate regulation, also used as indicator for ecosystem
health), C storage (natural capital that underpins and is closely related to C sequestration and global climate regulation), crop production
(contributing to food supply), annual water runoff (indicator for water availability but also related to flood regulation) and evapotranspiration
(indicator for regional climate regulation). C storage was investigated for vegetation, soil and litter C as well as total C, with the latter also including
carbon stored via carbon capture and storage (CCS) in BECCS scenarios. All plant and crop functional types contribute to an ecosystems' NPP. Therefore, NPP and crop production are
positively correlated. All variables are direct outputs of LPJ-GUESS simulations. Baseline crop yields of each run were scaled to FAO observed yields
in 1997–2003 as in Krause et al. (2017). Yields respond to changes in climate and
Absolute land area of croplands, pastures and natural areas between 1850 and 2100 for 16 scenarios of five land-use models, and the detailed relative changes in LU from 2000 to 2040 analyzed in this study.
Categories of dominant land-use change from 2000–2004 to 2036–2040 for each of 16 land-use scenarios. The legend is identical for all plots.
Total area of cropland, pasture and natural land as well as the change therein from 2000–2004 to 2036–2040 for 16 land-use scenarios. For each scenario, the left column gives the global total area for 2000–2004 (upper value) and 2036–2040 (lower value) and the right column gives the change from 2000–2004 to 2036–2040 in absolute terms (upper value) and as a percentage relative to the level in 2000–2004 (lower value). Gray shading indicates a positive or negative trend. Total area under change from 2000 to 2040 is given in absolute terms and as a percentage of the global ice-free land area considered in this study (see Sect. 2). Minor deviations in numbers may occur due to rounding.
Figures 1 and 2 summarize LULC changes in the 16 projections. The scenarios projected total changes in LULC of about 4.5 % to 11.4 % of the global ice-free land surface from 2000 to 2040 (Table 2), corresponding to different transitions between crop, pasture and natural land. While different socioeconomic assumptions realized as land-use change projections by the same model led only to small variations in the outcomes in terms of absolute rates (see Fig. 1, Table 2) and spatial patterns of LULC changes (see Fig. 2), the variation in LULC change between models was much more important, including for similar socioeconomic scenarios. This highlights the importance of the differences in modeling strategies. In this regard, we realize that the outcomes as described are indicative of the models' behaviors for the particular scenarios considered in this study, which are not necessarily but very likely representative of the models' general behaviors upon projecting future LULC patterns. LUH1 and 2 are exceptions with regard to the changes in between scenarios, because individual scenario data originate from different LUMs; therefore, their data differ substantially between all scenarios.
In comparison to the other LUMs, CLUMondo shows rather small-scale LULC changes spread across large parts of the world (Fig. 2). The three scenarios have the same demands for livestock, crop production, and so on, and the additional objectives in terms of C uptake and storage or biodiversity conservation did not introduce large variations. Therefore, differences between the CLUMondo scenarios are small (see Table 2). The biodiversity scenario leads to the most land area changes due to land system classes that diverge to either system intensification in some regions or extensification in others. Demands are fulfilled by almost linear trends until 2040 based on the assumed scenario storyline and demand estimates (Fig. 1).
Land-use changes in the three IMAGE scenarios also affect most of the productive land areas globally. Compared with CLUMondo, spatial patterns are different and the spread of percentage area changes across scenarios is higher (globally 5 %–10 % difference per LULC class for IMAGE scenarios; Table 2). Global trends are not linear; some scenarios even reverse their historic trend (e.g., IMAGE_ADAFF scenario for pasture) or accelerate it (e.g., IMAGE_BECCS scenario for cropland), possibly driven by the introduction of new land-use policies. In IMAGE, food production meeting the underlying societal demand has large priority. The IMAGE_Base scenario accordingly increases pasture area at the cost of natural land (presumably to satisfy demand for animal products in the underlying “SSP2” world) while cropland increases only slightly, presumably because yield increases satisfy the increasing food demand of a growing population. The afforestation and reforestation scenario IMAGE_ADAFF partially reverses IMAGE_Base by expanding natural land at the cost of pasture land, whereas IMAGE_BECCS, in addition to pasture expansion as in IMAGE_Base, also expands cropland areas at the cost of natural land. Spatially, the distribution of land-use classes among all scenarios differs little, indicating that demands from the IMAGE_Base scenario (e.g., population growth, diets, food demand, trade) outweigh specific scenario demands.
In MAgPIE, land changes only occur in specific regions or countries, although the changes in these regions are massive (southeastern Argentina and southern Brazil, some countries in eastern Africa, and parts of southern and eastern Asia), with the dominant change by far being cropland expansion. The three MAgPIE scenarios differ relatively little in time and space: only the afforestation and reforestation scenario again shows some very local natural area expansion. Trends over time are linear. Decisions regarding where land-use change takes place to meet food and feed demand strongly depend on minimizing the costs of land conversion. Here, some countries seem to provide substantially cheaper commodity prices than others, explaining the radical changes seen in the regions as listed above (compare also Fig. 2). It is noteworthy that MAgPIE and IMAGE derive potential crop yields and ecosystem C densities from the same DGVM (LPJmL, Bondeau et al., 2007), even though internal yield scaling and forest growth curves are implemented differently. However, their spatial patterns are quite different, emphasizing the role of individual strategies to translate demands under similar biophysical constraints into LULC patterns. Furthermore, the land demand to meet the same carbon dioxide removal (CDR) target was found to be larger in IMAGE than in MAgPIE (Krause et al., 2018).
In contrast, land changes in all LUH1 scenarios are large and occur in most of the productive land areas globally, reflecting both the highly diverging socioeconomic storylines as well as their implementation by different IAMs (see Table 2). Trends over time are nonlinear but involve multiple break points or gradual slopes. Interestingly, LUH1_26Be, which was developed by the IMAGE model (Hurtt et al., 2011), focuses on a broad expansion of croplands in tropical regions, whereas IMAGE_BECCS (although most likely implementing a different degree of bioenergy growth) includes a massive relocation of pastures and croplands to tropical and subtropical areas, respectively (Fig. 2). LUH1_45Aff and even more so LUH1_60Stab focus on the massive expansion of natural areas in all global regions where forests can be sustained, whereas LUH1_85Pop expands pastures and secondarily croplands in tropical and subtropical areas. The attribution of specific spatial LULC patterns to model allocation strategies vs. scenario storylines is impossible for LUH1, and in the same way also for LUH2, because underlying IAMs and storylines differ between each scenario.
LUH2 scenarios also differ substantially, corresponding to the very different SSP storylines and RCPs combined with their origin from different IAMs. LUH2's SSP1-26 (also implemented by IMAGE but with different socioeconomic assumptions than IMAGE_BECCS and LUH1_26Be), which includes options for both bioenergy crops and forest regrowth, shows expansion of natural areas mostly in temperate (and some boreal) regions of the northern latitudes and also in Australia, some cropland expansion and a reduction in pastures (Figs. 1 and 2). LUH2_SSP3-70, in contrast, results in a massive cropland expansion in some regions, in combination with a relocation of pastures, while LUH2_SSP5-85 (implemented by REMIND/MAgPIE) shows very large and concentrated regional dynamics, with cropland expansion similar to the MAgPIE scenarios as presented above.
In summary, most scenarios only agreed on a trend for cropland expansion at the cost of natural or pasture areas in terms of total area (see also
Fig. 1). Moreover, the scenarios showed very diverse patterns with respect to where and how these changes were realized. Thus, the deviation in LULC changes from
2000–2004 to 2036–2040 across the scenarios (Fig. S1 in the Supplement) showed major disagreement for cropland, pasture and natural
areas. Standard deviations of changes in land area
Change and uncertainty in ES indicators from 2000–2004 to 2036–2040 as a percentage relative to the base level in 2000–2004 across 16 land-use scenarios. Black bars give the average and standard deviation of the relative changes across all scenarios. See Table S3 for absolute levels and changes for each scenario.
The left column shows the categories of the average level of the provision of selected ES indicators in 2000–2004 and the relative change until 2036–2040 averaged over 16 land-use scenarios. Thresholds for categories for the average ES indicator level follow the 33rd and 67th percentiles for each ES indicator, whereas the change in indicators is given in 5 % steps for all indicators to allow for comparability. Note that for soil and litter C, the lowest category is negative. At high levels of ES indicator provision (highest 33 % of values), blue cells mark regions where little change is expected on average until 2040, and red cells mark regions where high changes are expected until 2040. The yellow category marks regions where base levels in ES provision are low; therefore, relative changes in these regions can be very high but are of minor importance. The percentage of global land area in each category is indicated. Cells where the average indicator level in 2000–2004 is zero are colored white and excluded from the statistical analysis. The right column gives the variability of the percent change in each ES indicator for each cell, which was calculated as the standard deviation of the changes in the ES indicator from 2000–2004 to 2036–2040 that was derived for each of the 16 land-use scenarios individually. Regions where the base level in 2000–2004 was below the 33rd percentile (yellow to green cells in the left column) were excluded from regional variability maps and colored gray in order to focus on cells with relevant ES indicator provision. Note that the legend scaling is different for vegetation C and crop yield production. Purple regions indicate a standard deviation in the predicted relative changes in this indicator higher than 10 % of the indicator level in 2000–2004 (30 % for vegetation C and 90 % for crop yield production). See Fig. S2 for NPP, total C storage and evapotranspiration.
The 16 land-use scenarios resulted in very diverse levels of ES indicators in 2000–2004 and changes therein until 2036–2040 simulated with LPJ-GUESS
(see Fig. 3 and Table S3 in the Supplement for all results given in the following). Figure 4 shows the spatial distribution of categories in ES
indicator levels and their changes until 2036–2040, averaged across the 16 scenarios. We decided to also investigate averages in order to explore some overall
emerging trends in ES indicators that result from the combined effects of climate and land-use change on ecosystem functionality. The average maps are
complemented by the regional variability in ES indicators (right column in Fig. 4, see also Sect. 4.2) as a measure of the large between-scenario
variability in ES indicators. Where regional variability is low, differences in LULC across scenarios are small and ES indicator changes can solely be
attributed to climatic changes and/or changing
The declining trend in natural areas (average decline of 0.9 %
Crop production was simulated to increase on average across all 16 scenarios by 31.2 %
Annual water runoff was simulated to increase on average by 9.3 %
Changes in evapotranspiration are closely linked to the calculations of runoff, although their effects are opposed, with higher evapotranspiration
rates contributing to reduced surface runoff (e.g., Piao et al., 2007) but also to biophysical cooling (e.g., Anderson et al., 2011)
Evapotranspiration rates increased in the CLUMondo, LUH1 and LUH2 scenarios on average by 2.6 %
The 16 scenarios project total changes in LULC of between 4.5 % and 11.4 % of the global ice-free land surface from 2000 to 2040. These rates are of a magnitude comparable to those observed for the past. For instance, four historical LULC reconstructions for the 1960–2000 period estimated changes in LULC of 7.5 %–12.8 % (Hurtt et al., 2011; Klein Goldewijk, 2016; Klein Goldewijk et al., 2011; Ramankutty et al., 2008, see Table S4 in the Supplement). A more recent global historical LULC change reconstruction by Winkler et al. (2021) estimated net LULC changes of 13.8 % for 1960–2015. This reconstruction uses a data-driven approach that is strongly based on remotely sensed information and provides higher spatial, temporal and thematic resolution than previous reconstructions. For a shorter time period, Liu et al. (2018) identified 3.4 % net LULC changes for 1992–2015 based on the ESA CCI Land-cover product. The picture becomes more complex when gross transitions, rather than net transitions, are considered, as gross changes (e.g., from forest to cropland in parts of a grid location and cropland to forest in another, over one time step) can be substantially larger than the net change. For instance, the change rates of Winkler et al. (2021) are 36.4 % when gross and multiple LULC changes are being considered individually. By contrast, LUM projections have a simplified representation of these more realistic LULC dynamics, omitting short-term, two-directional or small-scale transitions (such as under shifting cultivation, e.g., Heinimann et al., 2017). The improved representation of gross land-use changes in historical and future LULC reconstructions would be an important development to better account for LULC dynamics and their impacts on ecosystems (e.g., Bayer et al., 2017). However, such efforts are still hampered by a limited process understanding and data availability at the global scale.
Even given different initial states in LULC and different socioeconomic pathways of the 16 scenarios, we critically assess the spatial patterns, directions and rates of regional change based on past LULC changes. While it is not completely impossible, of course, we argue that a speed and magnitude which extremely oppose trends observed in the past seem at least questionable. Scenarios projecting future demands under reference or business-as-usual assumptions (CLUMondo_FAOref and IMAGE/MAgPIE_Base) are expected to continue recent historic LULC trends at least during the first part of the simulation period. Nevertheless, economic growth assumptions, demographic considerations and limitations with respect to land availability in the scenarios will likely cause some divergence from historical LULC trends. Historic trends are indeed continued in some of the broader regional projected trends, such as the expansion of natural areas in mid to high northern latitudes (CLUMondo_FAOref and IMAGE_Base; see Fig. 2) or the expansion of croplands in subtropical areas (CLUMondo_FAOref and to some degree also IMAGE/MAgPIE_Base) (e.g., Hansen, 2013). In contrast, some projected LULC changes have no historic precedent, such as the extensive pasture increase in tropical Africa (IMAGE_Base), the large cropland expansion in some selected African countries (MAgPIE_Base), or cropland expansion in the Mediterranean and Middle East (IMAGE_Base).
Scenarios including more drastic changes in the socioeconomic system (e.g., all LUH1 and LUH2 scenarios), including specific conservation measures or large-scale land-based climate change mitigation efforts (e.g., CLUMondo_CStor/Bdiv and IMAGE/MAgPIE_ADAFF/BECCS), affect future LULC patterns in very different ways, compared with reference or business-as-usual scenarios. Whether or not simulated abrupt LULC changes or a rapid reversion of regional historic LULC trends are realistic is difficult to judge, as analogous historical evidence is scarce or absent. Indeed, some rapid land-use changes have occurred in the past, caused by unexpected disruptions in markets or governance structures such as Brazil's soy moratorium combined with the enforcement of related policies (e.g., Nepstad et al., 2014; Gibbs et al., 2015) or the collapse of the Soviet Union (e.g. Hostert et al., 2011). However, capturing such unexpected LULC changes in global LUM projections is nearly impossible. Still, in response to most policy interventions, work has suggested that transitions in land use across regions tend to occur rather smoothly and with time lags of years to few decades (i.e., spanning a notable part of our simulation period) due to delayed policy uptake (e.g., Brown et al., 2019). In this context, large-area and relatively rapid regional change rates could be assessed critically, such as (1) forest regrowth on pasture and cropland areas with more than a 40 % area change from 2000–2004 to 2036–2040 in southeastern South America (LUH1_45Aff) or entire subtropical Africa (LUH1_60Stab); (2) massive cropland increases exceeding a 40 % total area change, e.g., in southeastern South America and eastern Africa (MAgPIE_ADAFF/BECCS and LUH2_SSP5-85); and (3) pasture expansion exceeding 20 % of the total area in tropical and subtropical Africa (e.g., IMAGE_BECCS, LUH2_SSP3-70). The first example reverses current deforestation trends (compare, e.g., Curtis et al., 2018; Hansen et al., 2013), whereas in the latter two examples, simulated rates of change are substantially larger than observed in these regions in recent decades (e.g., Hansen et al., 2013; Klein Goldewijk, 2016; compare Fig. S3 in the Supplement). LULC scenarios assuming significant amounts of forest regrowth should also be seen in light of the recent evidence provided by Holl and Brancalion (2020), pointing out manifold problems attached to tree planting and, therefore, calling for the prioritization of natural forest protection.
In general, CLUMondo and IMAGE were more capable of capturing small-scale changes within heterogeneous regions. Given the complexity in which land changes are being observed (e.g., Curtis et al., 2018; Hansen et al., 2013), the capacity to simulate changes at the small scale is likely more realistic. For instance, for a case study in Brazil, Stepanov et al. (2020) found that a spatially explicit regionally LUM creating small-scale changes simulated observed LULC patterns much better than the larger-scale changes in an economic model. By contrast, regional patterns in other LUMs tended to be fairly broad in extent or were limited to regions that were small in area. This is not consistent with what can be learned from global multi-temporal remote sensing.
When assessing the plausibility of regional LULC allocations, a good indicator is the expansion of agricultural areas (e.g., Salmon et al., 2015) based on existing yield potential (i.e., fertile areas currently not used or with a low share of croplands) or based on existing yield gaps (i.e., because of poor management or limits posed by the socioeconomic environment). We find that the LUMs use the regions with existing yield capacities in different ways for their allocation of LULC types. Cropland areas with currently relatively large yield gaps, such as in Brazil's Cerrado and west or east Africa (e.g., Mueller et al., 2012), were used by CLUMondo for further cropland expansion, which also seems plausible considering past LULC trends and continuing economic growth in response to an increasing population in these regions. In IMAGE simulations (especially IMAGE_Base/BECCS), these regions were typically converted into pastures and croplands were instead expanded in northern Africa and western Asia (especially Syria and Iraq). Besides current political turmoil in north Africa and west Asia, these allocations seem less plausible considering rather unfertile soils and existing yield gaps that are predominantly linked to low fertilizer inputs and water scarcity (e.g., Pala et al., 2011). In addition, the expansion of cropland area in IMAGE in northern Africa (especially Libya and Egypt; see Fig. 2) notably exceeds the current cropland extent (e.g., Fritz et al., 2015). It seems doubtful that this can indeed be achieved given that existing yield potential in these regions is low due to biophysical limitations during the crop-growing season which is not projected to change in the future. Other examples where LUMs diverge notably include central India, where all IMAGE scenarios used areas with large yield gaps as indicated by Mueller et al. (2012) for cropland expansion. In contrast, all CLUMondo scenarios used these croplands for forest regrowth, thereby ignoring their potential to contribute to fulfill increasing food demand. The massive (but very regional) cropland expansion of all MAgPIE scenarios and LUH2_SSP5-85 includes regions in eastern Africa, where yield gaps are high, but also regions in southeastern South America that are already under cropland usage to a high degree and attained yields are high (e.g., Fritz et al., 2015; Mueller et al., 2012). This suggests that in the MAgPIE model, economic considerations dominate the allocation of LULC classes rather than existing biophysical capacities. Only a few scenarios (e.g., LUH1_26Be, LUH2_SSP1-26 and LUH2_SSP5-85) expand cropland area in continental eastern Europe, especially along the “Chernozem Belt” into Russia, where soils are fertile and closing yield gaps would be expected to lead to large returns in terms of enhanced productivity (see Mueller et al., 2012). Other scenarios saw no potential for cropland expansion in these regions and simulated forest regrowth.
It has to be noted that in regions where attained yields are already currently relatively close to potential yields (e.g., northeastern North America, western Europe, and some parts of southern and southeastern Asia; see Mueller et al., 2012), yields may decline in the future due to climate change (e.g., Elliott et al., 2014; Funk and Brown, 2009; Lobell et al., 2009; Moore and Lobell, 2015; Pugh et al., 2016), and unless this can be counteracted by different management (e.g., increased irrigation) or different crop varieties, production may need to shift to other regions. However, the degree of climate change over our simulation period is too small to discern negative impacts on yields as well as associated climate-driven crop area changes.
Aside from biophysical considerations, LULC changes in response to changing economic conditions as projected in the scenarios require support from adequate regional policy, production and trading systems as well as appropriate technological capabilities (e.g., Lambin et al., 2003; Meyfroidt et al., 2019). Weak governance structures, for instance, would allow highly market-driven LULC changes, especially arising from changing demand elsewhere. LUM projections involving large-scale separation of LULC types across global regions (i.e., large regions of just one LULC class, e.g., IMAGE_Base/BECCS, all MAgPIE scenarios and LUH1_60Stab, LUH2_SSP5-85) could be interpreted to be detrimental to regional to national food production systems. This would include, for instance, subsistence farming systems over wide parts of Africa, which are an essential pillar of Africa's food supply and will continue to be so in the future (e.g., Sulser et al., 2015). In addition, technological capacities would have to be in place to support crop production at the projected locations, possibly including factors such as irrigation, the use of appropriate machinery and fertilization (see, e.g., Barrett and Toman, 2010; Lambin et al., 2014; Nilsson and Persson, 2012; Wang et al., 2016). At present, no global LUM is set up to consider governance aspects, such as land tenure rights or location-specific management, or transportation and trade capacities (apart from assumptions made in the economic core of LUMs that apply to large regions), which emphasizes an important need for further development.
To better understand regional LULC allocations and the related impacts on ecosystems and ES indicators, the availability of more spatial information from LUMs related to regional-scale assumptions on technological progress, flows of food import, export and local production would be useful. In this context it seems worthwhile for the land-use community to evaluate simulated future land-use changes against historic trends in spatial, temporal and thematic aspects. This may avoid some of the questionable, possibly unrealistic, land-use change effects seen in this study. However, current data products of historic land-use change are often themselves associated with high uncertainty in historic trends, due to data limitations. Improved historical products that merge multiple data sources could support the evaluation of future projected LULC changes. In addition, a clear declaration of scenarios showing potentially possible regional or global LULC pathways in comparison to those showing LULC changes going beyond historical exemplars would be desirable.
It is well known that different climate trajectories (e.g., for different RCPs) will greatly affect ecosystems. Even climate change projections for a single RCP when realized with different Earth system models will result in large variability in computed ecosystem outcomes (e.g., Ahlström et al., 2012). The fact that similarly large variability can be introduced by land-use change (within or between, for example, an SSP) is less established and, as such, an important outcome of this study.
The LULC patterns observed in the 16 scenarios suggested a general prioritization of food and feed demand affecting croplands and pastures (e.g., production of crops, meat) over those related to natural land dynamics (e.g., C storage, biodiversity). This is not surprising given that aspects that could impact “non-food” demands, such as C prices, are not considered in many of the scenarios' underlying storylines. In LUMs where both a baseline scenario and scenarios with additional demands were simulated (CLUMondo, IMAGE and MAgPIE), the overall trends in LULC and ES indicator changes were chiefly determined by the demands of the baseline scenarios and their model-specific implementation. Specific additional demands, such as for land-based mitigation or conservation, mostly resulted in only small deviations of LULC patterns and ES changes from the baseline scenarios. Deviations between scenarios of different models were much larger.
Scenarios that included specific climate change mitigation targets (CLUMondo_CStor, IMAGE_ADAFF/BECCS and MAgPIE_ADAFF/BECCS) resulted in larger total
C storage in LPJ-GUESS but lower crop production compared with the baseline scenarios. In LUMs, the mitigation scenarios include technology-driven
yield increases, which are higher than those assumed in LPJ-GUESS (see Sect. 2). From a C storage perspective, concern about the C storage potential
calculated by LUMs was raised by Krause et al. (2018), who could not reproduce the cumulative C uptake that was calculated in IMAGE and MAgPIE BECCS
and ADAFF when applying their LULC change to ecosystem models. On average only 62 % of C uptake was achieved with LPJ-GUESS, and only about 55 % was achieved
when three other DGVMs were used in addition. Likewise, Harper et al. (2018) also found the IMAGE C storage potential from BECCS to be achieved by
less than 25 % when two IMAGE mitigation scenarios were simulated with the JULES DGVM. These discrepancies likely arise from different assumptions in
LUMs and DGVMs related to growth rates and C uptake of regrowing forests and bioenergy crops, changes in soil C stocks upon LULC change, legacy
effects of previous land-use changes and some further processes such as disturbances (e.g., forest fires). By contrast, when using the CLUMondo LULC
patterns, all three scenarios led to an increase in total C storage in LPJ-GUESS, thereby even exceeding the no-net carbon loss target that was
implemented in CLUMondo_CStor. This might be explained by the joint effects of N deposition,
A number of studies have begun to identify in more detail how different assumptions in LUMs might affect LULC projections. Stehfest et al. (2019) recently provided a comprehensive sensitivity analysis of the socioeconomic drivers that were projected across five SSP-based storylines by six agro-economic models/IAMs and found very diverging sensitivities across models. The study highlights the existing variability in LULC modeling and emphasizes the need for more empirical research on crucial factors in the LULC modeling process such as long-term drivers of LULC change or the representation of land-use regulation and trade. The spread between LULC projections of different models could also be reduced by joint calibration and validation standards, but these do not currently exist (e.g., van Vliet et al., 2016). Reducing the spread between models, especially for similar scenario assumptions, would provide an important step forward in understanding LUMs' LULC patterns and in identifying possibly implausible allocations and would also support the assessment of calculated impacts on ecosystem functioning and ES indicators. A comprehensive comparison of the socioeconomic drivers of LUMs with the ES indicator levels simulated with DGVMs based on the LULC patterns from LUMs would be needed. This would provide deeper insights into the dependencies between drivers, modeling strategies and resulting ES provisions and would contribute to identifying the quality of the representation of interactions between socioeconomic and environmental systems including relevant feedback mechanisms.
Despite some very diverging features, the 16 land-use scenarios, on the level of global totals (see Fig. 3, Tables 2 and S3), resulted
in a positive change in NPP, vegetation C (all except one scenario), crop production and annual water runoff, and in a negative change in soil and
litter C stocks (all except one scenario) from 2000–2004 to 2036–2040 when simulated with LPJ-GUESS. Diverging trends were predicted for
evapotranspiration. Emerging trends are the result of joint climate change, increasing
Regionally, large variability in ES indicators (see right column in Figs. 4 and S2 in the Supplement) reflected the diverging LULC scenarios (effects
from scenario storylines and model-specific implementation) and their interactions with climatic changes. Only in areas colored white does the absence of
variability in changes in ES indicators indicate dominating climate change and
Tropical forest regions, or at least major regions therein, were identified as hotspots of variability across the LULC projections with large areas
showing variability of
In tropical savannas and temperate shrubland and grassland regions, the diverging LULC projections resulted in high variabilities in vegetation productivity, vegetation C storage, annual water runoff and crop production. Although 15 out of 16 land-use scenarios increase cropland area in these regions, they diverged widely with respect to the exact location of cropland expansion, resulting in these high variabilities in ES indicators. Savannas have been highlighted before as particularly vulnerable to future conversions of natural vegetation into cropland or pasture (e.g., Shin et al., 2019) because of large population growth in many savanna regions, their climatic suitability for agriculture and the relatively small effort needed for conversion considering the relatively low woody vegetation cover. Parts of these regions have already experienced intense cropland expansion in recent decades, such as South America's Cerrado and Chaco or African savannas (e.g., Aleman et al., 2016; Hansen et al., 2013; Noojipady et al., 2017). Processes affecting ecosystem functionality and services in these subtropical to semiarid regions are also particularly important in view of rapid population growth and associated demands (Alexandratos and Bruinsma, 2012) as well as their role in the global C cycle and climate dynamics (e.g., dominant role for the trend and interannual variability of the global land C sink, see Ahlström et al., 2015, with their semiarid class widely corresponding to our tropical savannas and temperate shrublands and grasslands classes). Therefore, high variability in the provision of ES indicators related to vegetation productivity, water availability and food supply may have severe consequences on ecological, economic and social systems in these regions.
In both temperate and boreal forest regions, the regional variability of vegetation productivity and vegetation C storage was high because the increased photosynthetic productivity and longer growing seasons under warmer climates in higher latitudes were either dampened or amplified by diverging LULC projections. While LULC changes dominate this combined positive effect in temperate regions, they are much smaller in boreal regions and the climatic effect dominates (see Table S5 in the Supplement, compare also high variability of changes in NPP, vegetation and soil/litter C in boreal regions for multiple DGVMs, GCMs and emission pathways in Nishina et al., 2015, under fixed LULC). Counteracting these positive trends, studies highlight new challenges for temperate and boreal forests emerging from climatic changes and anthropogenic disturbance with the capacity for severe ecosystem-level damage, such as droughts, insects, fire regimes and pathogens (e.g., de Groot et al., 2012; Millar and Stephenson, 2015; Park et al., 2014). As only fire is explicitly simulated as a disturbance process in LPJ-GUESS, whereas other forms of disturbance are subsumed in a stochastic background disturbance, this could further increase the regional variability in ES indicators (see, e.g., Pugh et al., 2019). Regional variability was also high for crop production due to diverging extensification and intensification trends across temperate and boreal regions. Land-use changes and associated regional variability in ES indicators in the cold (tundra) and warm desert regions are not significant in a global context because of the climatically constrained low productivity, although these effects are enhanced at least in tundra regions through warmer temperatures.
A direct correlation of the per grid cell changes in ES indicators with the corresponding changes in cropland, pasture and natural land fraction could
reveal the sensitivity of different ES indicators to changes in LULC. Across all scenarios and for the considered biomes, these relationships
suggested, for instance, an approximate 1.5 % increase in vegetation C per percent increase in natural land fraction and between a 12 % and 24 %
increase in crop production per percent increase in the cropland fraction (see Table S8 in the Supplement). Emergent responses of other ES indicators
to changes in LULC fractions were mostly low (slope of regression lines close to 0) across the biomes (see Table S8). None of the identified
relationships provided high reliability (highest
The IPBES report on plausible futures of nature identified some general trends across their scenarios (e.g., continued increase in managed land, increases in material and decreases in regulating and nonmaterial nature's contributions to people) and rated the knowledge base and confidence in effects of interactions of future LULC and climatic changes on biodiversity and ecosystem functioning as “established but incomplete” (see Shin et al., 2019). By evaluating 16 scenarios of five structurally different LUMs, we covered a large variability existing in currently available spatially explicit projections of LULC and, therefore, contributed to extending the existing knowledge base on the variability of future ES indicators. However, conclusions drawn here with regard to projected changes in LULC and ES indicators are inherently dependent on the selected set of LUMs and scenarios, the evaluation time period and the simulation setup. In addition, we did not consider climatic variability resulting from different GCM implementations of climate under RCP2.6 nor other possible emission pathways, which both add uncertainty to future ES indicator levels. Previous studies (e.g., Friend et al., 2014; Krause et al., 2019; Nishina et al., 2015; Pugh et al., 2018) have explored the variability of using different ecosystem models with the same or multiple LULC and climate pathways, whereas Ahlström et al. (2012) and Schaphoff et al. (2006) used climate forcing data from multiple GCMs following the same emission pathway instead of only one climate model to quantify ES indicators. Our sensitivity test using climate inputs from five GCMs and the four LUH1 scenarios indicated uncertainties of between 0.1 % and 8.6 % for the ES indicators considered here (Table S6). The fact that crop type distribution and N fertilization was taken directly from the LUMs for IMAGE and MAgPIE simulations was considered as an extension of the LULC input from the LUMs and, therefore, as a part of the LULC scenarios of these models. In this study, we considered the near to medium future until 2040. Estimates of ES indicator levels and their variability for the far future until 2100 are likely different.
We conclude that LUMs and IAMs have some fundamental limitations in capturing all relevant processes related to LULC changes which in some scenarios result in questionable and potentially unrealistic features in their regional LULC allocations and their global and regional trends. More spatial information from LUMs related to regional assumptions (e.g., on technological progress, flows of food import, export and local production) in addition to LULC patterns would be helpful to better understand their regional LULC allocations. Limitations also include technical aspects of LUMs, such as adequate spatial resolution, representation of net vs. gross LULC changes and the number of LULC classes. We can only deduce a generally high influence of model-specific singularities from our findings, although we could not attribute the discrepancy in LULC patterns to the individual factors in the LULC modeling process based on the available LUM data.
The variability across the 16 LULC scenarios entails a high variability in the trends of most ES indicators investigated in this study. On a regional level, this emphasizes the role of tropical forests, especially the borders of the currently forested tropics, as regions with the most uncertain future developments that have the potential to significantly alter regional ecosystem functionality and services in a way that is substantial in a global context. Among the investigated ES indicators, the identification of crop production as the indicator that was associated with the highest uncertainties by far in terms of global totals and regional variabilities, highlights the diverging production targets as well as the different regional strategies of the investigated models to fulfill future demands for crop production.
Our results stress that information from individual LUMs or IAMs, which is used for policy support towards sustainable transformations, should be complemented by further information such as the variability of ES indicators arising through different LULC projections. This would provide a wider context that is essential with respect to policy making. For regional decision-making, this is especially important in regions with highly diverging trends in land-use scenarios and, therefore, large variability in ES indicators. The issue may become even more relevant in the second half of the century, when LULC changes for climate change mitigation and adaptation are likely to intensify. Ultimately, we need to find improved ways to achieve a better integration between models targeting the different aspects in the cycle of socioeconomic developments and their direct, indirect and cumulative implications on natural systems.
LUH2 data are available at
The supplement related to this article is available online at:
ADB and AA conceived and designed the experiments. PHV and RF made CLUMondo data available and converted them for this study. PA prepared LUH2 input data for LPJ-GUESS and assisted with the LPJ-GUESS model simulations. RM, AK and ADB carried out the model simulations and led the data analysis with contributions from all authors. RF contributed to the analysis of the spatial–temporal patterns of land-use models. ADB led the writing of the paper with contributions from all authors.
The authors declare that they have no conflict of interest.
The authors thank the research groups for making climate model data (ISI-MIP) and land-use model data available to the community. We thank Karina Winkler for providing data on historical LULC changes.
This research has been supported by the European Commission's 7th Framework Programme (grant no. 308393, OPERAs, and grant no. 603542, LUC4C), and the Helmholtz Association and its ATMO research program. The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
This paper was edited by Somnath Baidya Roy and reviewed by Edouard Davin and one anonymous referee.