Articles | Volume 18, issue 9
Biogeosciences, 18, 2957–2979, 2021
Biogeosciences, 18, 2957–2979, 2021

Research article 17 May 2021

Research article | 17 May 2021

Climate change and elevated CO2 favor forest over savanna under different future scenarios in South Asia

Climate change and elevated CO2 favor forest over savanna under different future scenarios in South Asia
Dushyant Kumar, Mirjam Pfeiffer, Camille Gaillard, Liam Langan, and Simon Scheiter Dushyant Kumar et al.
  • Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Senckenberganlage 25, 60325 Frankfurt am Main, Germany

Correspondence: Dushyant Kumar (


South Asian vegetation provides essential ecosystem services to the 1.7 billion inhabitants living in the region. However, biodiversity and ecosystem services are threatened by climate and land-use change. Understanding and assessing how ecosystems respond to simultaneous increases in atmospheric CO2 and future climate change is of vital importance to avoid undesired ecosystem change. Failed reaction to increasing CO2 and climate change will likely have severe consequences for biodiversity and humankind. Here, we used the adaptive dynamic global vegetation model version 2 (aDGVM2) to simulate vegetation dynamics in South Asia under RCP4.5 and RCP8.5, and we explored how the presence or absence of CO2 fertilization influences vegetation responses to climate change. Simulated vegetation under both representative concentration pathways (RCPs) without CO2 fertilization effects showed a decrease in tree dominance and biomass, whereas simulations with CO2 fertilization showed an increase in biomass, canopy cover, and tree height and a decrease in biome-specific evapotranspiration by the end of the 21st century. The predicted changes in aboveground biomass and canopy cover triggered transition towards tree-dominated biomes. We found that savanna regions are at high risk of woody encroachment and transitioning into forest. We also found transitions of deciduous forest to evergreen forest in the mountain regions. Vegetation types using C3 photosynthetic pathway were not saturated at current CO2 concentrations, and the model simulated a strong CO2 fertilization effect with the rising CO2. Hence, vegetation in the region has the potential to remain a carbon sink. Projections showed that the bioclimatic envelopes of biomes need adjustments to account for shifts caused by climate change and elevated CO2. The results of our study help to understand the regional climate–vegetation interactions and can support the development of regional strategies to preserve ecosystem services and biodiversity under elevated CO2 and climate change.

1 Introduction

Global climate has been identified as the primary determinant of large-scale natural vegetation patterns (Overpeck et al.1990). Climate change has affected global vegetation pattern in the past and caused numerous shifts in plant species distribution over the last few decades (Chen et al.2011; Thuiller et al.2008). It is expected to have even more pronounced effects in the future and may lead to drastically increasing species extinction rates in various ecosystems (Brodie et al.2014). Natural ecosystems have been and continue to be exposed to increased climate variability and abrupt changes caused by increased intensity and frequency of extreme events such as heat waves, drought and flooding (Herring et al.2018). At the same time, they are under severe pressure due to anthropogenic disturbance and land conversion. Rising levels of atmospheric CO2 are a strong driver of climate-induced vegetation changes (Allen et al.2014). Anthropogenic CO2 emissions account for approximately 66 % of the total anthropogenic greenhouse forcing (Forster et al.2007) and are thus largely responsible for contemporary and future global climate change (Parry et al.2007). Rising CO2 is expected to alter distributions of plant species and ecosystems (Parry et al.2007) both indirectly through its influence on global temperatures and precipitation patterns (Cao et al.2010), two main drivers of vegetation dynamics, and directly via its physiological effects on plants (Nolan et al.2018). It is therefore of vital importance to understand how ecosystems respond to simultaneous increases in atmospheric CO2 and temperature, to changes in precipitation regime, and to altered ecosystem water balance in order to avoid critical ecosystem disruptions and the resulting consequences for biodiversity and humankind.

Increases in temperatures, decreases in precipitation and changes in precipitation seasonality can cause loss of vegetation biomass. Plants using C3 photosynthetic pathway are often not saturated at the current atmospheric CO2, whereas plants using the C4 photosynthetic pathways are already at their physical optimum at current atmospheric CO2 levels (Ehleringer and Cerling2002). The physiology of C3 plants implies that elevated atmospheric CO2 improves their ability for carbon uptake due to the CO2 fertilization (Woodrow and Berry1988) and enhances carbon sequestration (Leakey et al.2009; Norby and Zak2011) as well as plant water use efficiency (Soh et al.2019). This has also been observed in long-term free-air carbon dioxide enrichment (FACE) experiments (Norby and Zak2011). Thus, elevated CO2 influences photosynthesis and thereby affects other physiological processes such as respiration, decomposition (Doherty et al.2010), evapotranspiration (ET) and biomass accumulation (Frank et al.2015). Increasing CO2 concentration has been associated with woody cover increase in structurally open tropical biomes such as grasslands and savannas (Stevens et al.2017). This widespread proliferation of woody plants into arid and semiarid ecosystems has been attributed to increased water use efficiency in C3 plants that facilitates woody sapling establishment and growth due to higher drought tolerance (Kgope et al.2010; Stevens et al.2017). These CO2 effects on plant growth and competition can alter community structure (height distribution), ecosystem productivity, climatic niches of ecosystems and biome boundaries (Nolan et al.2018; Wingfield2013). Change in vegetation distribution and altered vegetation structure feed back on climate by altering fluxes of energy, moisture, and CO2 between land and atmosphere (Friedlingstein et al.2006). Feedback mechanisms also involve vegetation‐mediated changes in albedo, surface roughness, land‐-atmosphere fluxes and evapotranspiration (Field et al.2007; Richardson et al.2013).

Enhanced plant growth due rising CO2 implies rapid leaf area development and more total leaf area could translate into higher transpiration (Leakey et al.2009). However, elevated CO2 concentrations may decrease leaf stomatal conductance to water vapor, which could reduce transpiration. Evapotranspiration (ET) is a key ecophysiological process in the soil–vegetation–atmosphere continuum (Feng et al.2017). Annually, 64 % of the total global land-based precipitation is returned to the atmosphere through ET (Zhang et al.2016). Environmental change and concurrent vegetation changes alter ET and affect water availability (Mao et al.2015), especially in arid and semiarid regions. In these regions, ET affects surface and subsurface processes such as cloud development, land surface temperature and groundwater recharge (Fisher et al.2011).

South Asia is home to approximately 1.7 billion people and is one of the regions most vulnerable to climate change (Eckstein et al.2018). It hosts four of the world's biodiversity hotspots (Myers et al.2000) and harbors different biome types ranging from tropical in the south to temperate in the north at the fringe of the Himalayas. These hotspots are characterized by high levels of diversity and endemism, and they are threatened by climate change and anthropogenic land use (Deb et al.2017). For instance, woody encroachment due to rising CO2 threatens South Asian savannas (Kumar et al.2020), and sifting cultivation in the northeastern part of South Asia threatens biodiversity (Bera et al.2006).

Due to the absence of long-term field experiments such as FACE experiments, in the dominant biomes of the region, modeling studies are valuable tools to close existing knowledge gaps. Dynamic global vegetation models (DGVMs, Prentice et al.2007) are particularly well suited to address questions that focus on vegetation response to changing environmental drivers, e.g., climate and CO2. While most DGVM studies in South Asia focused on the vulnerability of forests to climate change (Chaturvedi et al.2011; Ravindranath et al.2006, 1997), they often overlooked the severely threatened savanna biome. These studies were further limited by the utilization of models with fixed ecophysiological parameters and traits, e.g., fixed carbon allocation values to assign carbon to plant biomass pools, fixed specific leaf area (SLA) and pre-defined bioclimatic limits that were derived from contemporary climatology in order to constrain the spatial distribution of plant functional types (PFTs). Moreover, many DGVMs used in these studies do not account for life history, eco-evolutionary processes and trait variability among individual plants (Kumar and Scheiter2019). While some global-scale studies have investigated the potential effect of increasing CO2 on natural vegetation, carbon sequestration and biome boundaries (e.g., Hickler et al.2006; Sato et al.2007; Smith et al.2014), detailed modeling studies focusing explicitly on different biomes in South Asia have not been conducted. The physiological effects of increased CO2 and climate change on South Asian vegetation are uncertain and need to be addressed in order to improve understanding of regional ecosystem functioning as well as implications for biodiversity conservation.

To address the knowledge gaps in existing studies, we used the aDGVM2 (adaptive dynamic global vegetation model version 2), an individual- and trait-based vegetation model that combines elements of traditional DGVMs (Prentice et al.2007) with newly implemented approaches for selection and trait filtering. In aDGVM2, environmental conditions select for the plants with trait value combinations that make them successful under these conditions. Therefore, plant communities that are adapted to site-specific environmental conditions dynamically assemble and emerge as a reaction to the environmental forcing (Langan et al.2017; Scheiter et al.2013). Originally, aDGVM2 had been tested for Amazonia (Langan et al.2017) and Africa (Gaillard et al.2018; Pfeiffer et al.2019). In order to adapt it to South Asian ecosystems and their diversity, we included C3 grasses, improved ecophysiological processes such as the leaf energy budget in order to estimate leaf temperature, implemented separate temperature sensitivities for C3 and C4 photosynthetic capacity (Vcmax), and included snow in the water balance model.

In this study we used the updated version of aDGVM2 and addressed the following questions:

  1. How do projected changes in climate and CO2 following two representative concentration pathways (RCP8.5 and RCP4.5, Meinshausen et al.2011) change the distribution, boundaries and climatic niches of biomes in South Asia?

  2. How does the relationship between projected biomass, ET, temperature and precipitation change in response to CO2 fertilization?

  3. What is the sensitivity of predicted changes in relation to presence and absence of CO2 fertilization?

Based on our results we analyzed climate–vegetation interactions to improve our understanding of how to manage and mitigate impacts on biomes under climate change and increasing CO2.

2 Methods

2.1 Description of the study region

Approximately 1.7 billion people populate South Asia, i.e., the Indian subcontinent, Afghanistan and Myanmar. South Asia incorporates a wide range of bioclimatic zones with distinctive biomes, ecosystem types and species (Rodgers and Panwar1988). Climatic conditions are controlled by interactions between the South Asian summer monsoon system and the region's complex topography. The climatic envelope ranges from tropical arid and semiarid regions in the west to humid tropical regions supporting rain forests in the northeast and temperate vegetation at the fringe of the Himalayas. Excluding the Himalayan regions, South Asia has a mean annual temperature of approximately 24 C with very low spatial variability. Mean annual precipitation (MAP) is 1190 mm, ranging from less than 500 mm in the warm desert zone in the west to more than 3500 mm in the northeast. The steep elevation gradients ranging from sea level to 8800 m result in a rich diversity of ecosystems that can alternate in areas of a few hundred square kilometers. Topography is recognized as a strong driver of ecological patterns, for example those related to forest structure and composition, floristic diversity and soil fertility (Gallardo-Cruz et al.2009; Jucker et al.2018; Sinha et al.2018). South Asia hosts four major global biodiversity hotspots, namely the Western Ghats, Himalayas, India and Myanmar, and Sri Lanka (Myers et al.2000). These hotspots include a wide diversity of ecosystems such as mixed wet evergreen, dry evergreen, deciduous and montane forests. Further vegetation types are alluvial grasslands and subtropical broadleaf forests along the foothills of the Himalayas, temperate broadleaf forests in the mid-hills, mixed conifer and conifer forests in the higher hills, savanna in the Deccan region and southern part of Malaysia, and alpine meadows above the tree line.

2.2 Model description

For this study we used aDGVM2 (Scheiter et al.2013; Langan et al.2017; Gaillard et al.2018), a DGVM with a dynamic trait approach. In the Supplement we summarize main features of aDGVM2 and explain how the physiological effects of changing CO2 concentration and rising temperature are simulated in a process-based way in the aDGVM2 by the implemented photosynthesis routine. To adapt the aDGVM2 to the requirements of the study region, we incorporated new subroutines into the model. We improved the representation of (a) the water balance by including snow, (b) the carboxylation rate, and (c) leaf temperature, and (d) we included C3 grasses (previous model versions only simulated C4 grasses).

  • a.

    Water balance.

    In aDGVM2, the soil water module is based on the tipping-bucket concept. As the model was originally developed with a strong focus on tropical and subtropical forest and savanna regions, the original model version only considered water input in the form of rain (see Langan et al.2017). In the updated model version, precipitation is assigned as snow when daily mean air temperature drops below 0 C. Snow accumulates on the soil surface or is added on top of an existing snowpack. The snowpack persists as long as air temperature remains below 0 C. Once temperature rises above 0 C, water from snowmelt is added to the soil water pool and becomes available to plants. This process may improve the water availability for plants at the beginning of spring, for example in the Himalayan region. Snowmelt (Smelt, mm/day) is calculated following Choudhury et al. (1998) as

    (1) S melt = 1.5 + K m P precip ( T a - T snow ) S pack ,

    where Km is the coefficient of snowmelt (0.007 mm/day/C), Spack is the depth of the snowpack (mm) and is equivalent to the accumulated solid portion of precipitation, Ta is the daily mean air temperature (C), Pprecip is precipitation (mm/day), and Tsnow is the maximum temperature where precipitation falls as snow (0 C). We do not consider insulation effects of the snowpack in the model.

  • b.

    Carboxylation rate.

    In earlier versions of aDGVM2, leaf-level photosynthesis was calculated at the population level; i.e., it was assumed that all plants of a simulated vegetation stand have the same leaf-level photosynthetic rate. Only C3- and C4-type photosynthesis were distinguished. We therefore implemented new routines to calculate photosynthesis at a daily time step for each individual plant. We further incorporated an empirical relation between specific leaf area (ASLA, mm2/mg) and leaf nitrogen content per unit area (Na, g/m2) following Sakschewski et al. (2015):

    (2) N a = 6.89 A SLA - 0.571 .

    The standard maximum carboxylation rate of rubisco per leaf area (Vcmax,25, µmol/m2/s) was derived from the TRY database (Kattge and Knorr2007) by Sakschewski et al. (2015) and is calculated as

    (3) V cmax , 25 = 31.62 N a 0.501 ,

    where Vcmax,25 is Vcmax at 25 C. In the model, ASLA is linked to the matric potential at 50 % loss of xylem conductance (P50; see Langan et al.2017). The trade-off between ASLA and Vcmax mediated by leaf traits (Na) introduces variability in the spectrum of tree growth strategies in aDGVM2. ASLA is linked to leaf longevity (LL) in aDGVM2, such that it affects the leaf turnover rates (represented by Eq. 72, in Appendix, Langan et al., 2017). Leaves with high ASLA have shorter LL and higher turnover rates than leaves with low ASLA (and vice versa). The correlation between ASLA, P50 and LL represents the trade-off between two opposing resource strategies, i.e., conservation vs. rapid acquisition of soil water and nutrients (Wright et al.2005). Trees that invest more carbon into their leaves (low ASLA) enhance their structural stability but have lower leaf turnover to mitigate the higher initial carbon investment.

    The effect of temperature on photosynthesis is well described (Kirschbaum2004), and temperature may influence photosynthesis both directly, via temperature dependency of enzyme-mediated metabolic rates of carboxylation and the Calvin cycle (Sharkey et al.2007), and indirectly, via its effect on transpiration and plant water uptake and transport (Urban et al.2017). The maximum carboxylation rate (Vcmax) increases with temperature until it reaches an optimum, and it decreases again at temperatures above the optimum (Kattge and Knorr2007) due to reductions in enzyme activity. Above 30 C the electron transport chain is gradually inhibited, and at temperatures above 40 C the denaturation of rubisco and associated proteins becomes relevant (Lloyd et al.2008). The temperature dependency of the carboxylation rate (Vcmax) is expressed as

    (4) V cmax = V cmax , 25 2 0.1 ( T leaf - 25 ) ( 1 + e 0.3 ( T low - T leaf ) ) ( 1 + e 0.3 ( T leaf - T upp ) ) ,

    where Tleaf is the leaf temperature in C (see next paragraph for calculation). The photosynthetic model of Collatz et al. (1992) and Collatz et al. (1992) assumes specific values of Tupp and Tlow for C3 and C4 plants, respectively (Tables S1 and S2 in the Supplement). These temperature ranges from −10 to 36 C and from 13 to 45 C for C3 and C4 photosynthetic pathways, respectively, allow plants to grow most efficiently in their plant-specific climatic niches.

  • c.

    Leaf temperature.

    We calculate leaf temperature following the leaf-level energy budget concept (Gates1968). Leaf-level photosynthesis, activity of leaf enzymes and transpiration depend on leaf temperature (Tleaf, C), calculated as

    (5) T leaf = T air + R n - λ E r gb ρ C P ,

    where Tair is air temperature (C), Rn is net radiation absorbed by the leaf (MJ/m2/day), λ is latent heat of vaporization (MJ/kg), E is evapotranspiration (m/day), rgb is the boundary layer resistance (m/s), ρ is the air density (kg/m3) derived from atmospheric pressure (101.325 kPa at sea level) that is scaled according to the elevation and Tair, and CP is the specific heat of dry air (MJ/kg/C). Leaf temperature is used to calculate the temperature dependence of Vcmax used in the photosynthesis model routines in Eq. (4). Absorbed net radiation (Rn), rgb and E are model state variables calculated from climate input used in aDGVM2 (Tair, long-wave and shortwave radiation, and ρ). The values of latent heat of vaporization (λ) and CP are 2.45 MJ/kg and 2.71 MJ/kg/C, respectively, and are assumed to be constant parameters in this model version.

  • d.

    C3 grasses.

    C3 grasses were not included in previous aDGVM2 versions (Gaillard et al.2018; Langan et al.2017; Pfeiffer et al.2019; Scheiter et al.2013). We therefore implemented C3 grasses, following the approach used for C4 grasses in previous model versions, but adjusted the photosynthetic pathway (see Appendix S2 in Langan et al.2017). C3 and C4 grasses use a different leaf-level photosynthesis model (Farquhar et al.1980) following the implementations of Collatz et al. (1991, 1992). The optimum temperature ranges for carboxylation for C3 and C4 grasses are also different (Table S1). As C3 grasses have higher cold tolerance than C4 grasses (Liu and Osborne2008), we implemented frost intolerance for C4 grasses but not for C3 grasses. Frost is assumed to damage the tissue of C4 grasses, and in aDGVM2 we kill 10 % of the living leaf biomass of C4 grasses per frost day independent of frost severity.

2.3 Model forcing data

2.3.1 Climate data

We used GFDL-ESM2M climate data for the period 1950 to 2099 from the Inter-Sectoral Impact Model Inter-comparison Project (ISIMIP2), as historical climate simulated by GFDL-ESM2M showed satisfactory performance for South Asia (McSweeney and Jones2016). The general circulation model (GCM) output was bias-corrected in the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) and downscaled to a spatial resolution of 0.5× 0.5 (Warszawski et al.2014). We used average, maximum and minimum air temperatures, precipitation, surface downwelling shortwave radiation and long-wave radiation, near-surface wind speed, and relative humidity at a daily temporal resolution. We used two representative concentration pathways, namely RCP4.5 and RCP8.5 (Meinshausen et al.2011). These scenarios assume increases in radiative forcing of 4.5 and 8.5 W/m2 by 2100 (Van Vuuren et al.2011) and increases in atmospheric CO2 concentrations to 560 and 970 ppm by 2100, respectively (Van Vuuren et al.2011).

2.3.2 Projected changes in temperature and precipitation

Mean annual precipitation (MAP) from GFDL-ESM2M does not show a clear trend when averaged for South Asia under RCP4.5 and RCP8.5, due to high inter-annual variability of precipitation (Fig. S1 in the Supplement). Yet, there are region-specific differences in precipitation change. The Western Ghats which are located between 73–77 E and 8–21 N and the eastern Himalayan region are projected to become wetter under both RCP4.5 and RCP8.5, whereas the western part of the region is projected to become drier by the end of the century under both RCPs (Fig. S2). MAP is projected to increase by more than 600 mm in the eastern Himalayas and Western Ghats but predicted to decrease by 400–600 mm in the western and central area of the region (Fig. S2). By the end of the 21st century, the mean annual temperature (MAT) of South Asia is expected to increase between ca. 1 and 3.5 C under RCP4.5 and between 1 and 6 C under RCP8.5, relative to the average temperature in the baseline period of 2000–2009 (Figs. S1 and S2). The western parts of the region and the Himalayan mountains are projected to experience higher increases in temperature than the rest of the region (Fig. S2).

2.3.3 Soil and elevation data

Soil data were obtained from FAO (, last access: March 2016,  Nachtergaele et al.2009) and include information on soil properties and types. The soil properties include parameters required by aDGVM2: volumetric water-holding capacity, soil hydraulic conductivity, soil bulk density, soil depth, soil texture, soil carbon content, soil wilting point and field capacity (for details see Fig. S3b and Langan et al.2017). A digital elevation model (DEM) at 90 m spatial resolution was obtained from the Shuttle Radar Topography Mission (SRTM, Fig. S3a, last access: March 2016,  Jarvis et al.2008). It was resampled to a spatial resolution of 0.5× 0.5 to match the spatial resolution of climate data. In the model, elevation is used to calculate the surface pressure at a given altitude, which is used to scale up air density and partial pressure of oxygen. The partial pressure of oxygen is used to estimate the CO2 compensation point of photosynthesis (Eq. 2 of Appendix S2 in Langan et al.2017). We did not use slope and aspect in the model.

2.4 Model simulation protocol

To understand how climate and CO2 fertilization interact to influence the future vegetation state in South Asia, we simulated all combinations of two climate scenarios (RCP4.5 and RCP8.5) and two CO2 scenarios (CO2 fertilization enabled or disabled, four scenarios in total). We simulated potential natural vegetation between 1950 and 2099 using daily climate data for RCP4.5 and RCP8.5 (see Sect. 2.3.1). For both scenarios, simulations were run with a CO2 increase in line with RCP4.5 (hereafter RCP4.5+eCO2) and RCP8.5 (hereafter RCP8.5+eCO2) and with the same climate data but fixed CO2 after 2005 at 375 ppm for RCP4.5 (hereafter RCP4.5+fCO2) and RCP8.5 (hereafter RCP8.5+fCO2). Fixing the CO2 concentration after 2005 mimics a situation where CO2 fertilization would not occur and vegetation only responds to the climate signal. All simulations were conducted with natural fire as implemented in aDGVM2 and at 0.5× 0.5 spatial resolution. The aDGVM2 simulates 1 ha stands that are assumed to be representative for the vegetation at a larger scale; i.e., we assume that the stand-level vegetation homogeneously covers the grid cell. The “representative hectare approach” is a concession to computational limitation, as photosynthesis and physiological processes are simulated individually for all individual plants of a stand (up to 36 000 individuals). It balances adequate representation of trait diversity among individuals against computational constraints. Also due to computational constraints, we did not conduct replicate simulations.

To ensure that simulated vegetation had sufficient time to adapt to prevailing environmental conditions, we conducted simulations for 650 years, split into a 500-year spin-up phase and a 150-year transient phase. For the spin-up phase, we randomly sampled years of the first 30 years of daily climate data (1950 to 1979). For the transient phase, we used the sequence of daily climate data between 1950 and 2099. Trial simulations showed that a 500-year spin-up period is sufficient to ensure that vegetation is in a dynamic equilibrium state with environmental drivers.

2.5 Model benchmarking and evaluation

For benchmarking of aDGVM2 simulation results, we used five different remote sensing products: aboveground biomass (t/ha, Saatchi et al.2011), tree height (m, Simard et al.2011), tree cover (percent, Friedl et al.2010), MODIS evapotranspiration (mm/year, Zhang et al.2010) and natural vegetation type (Ramankutty et al.2010). We used a 10-year average of MODIS ET and compared it to a 10-year average of model-simulated ET (mm/year; 2000–2009). All remote sensing data sets were aggregated to a 0.5× 0.5 spatial resolution to match the spatial resolution of model simulations by calculating the mean of all values within each 0.5 grid cell or using nearest-neighbor aggregation in the case of vegetation type (“raster” package in R, Hijmans and van Etten2012). We first compared model results and observations assuming that the entire study region is covered by natural vegetation (Fig. 1). Then we repeated the comparisons only for areas with predominantly natural cover; i.e., we masked out areas with more than 50 % managed land (Fig. S4, land cover classes 7 “Cultivated and Managed Vegetation” and 9 “Urban and Built-up” in  Tuanmu and Jetz2014). We calculated the normalized mean squared error (NMSE) and coefficient of determination (R2) to quantify agreement between data and simulated variables.

Figure 1Comparison between aDGVM2 results and data for (a) simulated biomass and Saatchi et al. (2011) biomass, (b) simulated tree height and Simard et al. (2011) tree height, (c) simulated tree cover and Friedl et al. (2011) tree cover, and (d) simulated evapotranspiration and Zang et al. (2010) evapotranspiration. In the figure the first column shows the remote sensing products, the second column shows aDGVM2 results, the third column shows the difference between simulation and data, and the fourth column shows the scatter plot between simulated state variables against benchmarking data. NMSE and RMSE are the normalized mean square error and root mean square error, respectively. In the fourth column, each point represents one grid cell in the study region. For results with masked land-use cover see Supplement Fig. S4.

2.6 Biome classification

The aDGVM2 simulates state variables such as biomass and canopy cover of individual plants in simulated vegetation stands (1 ha, which is a representative of a grid cell). We used woody canopy area, abundance of shrubs and trees, and grass biomass to classify the simulated vegetation into biome types (Fig. S5). We used 10-year averages of state variables for the periods 2000–2009, 2050–2059 and 2090–2099 to represent the 2000s, 2050s and 2090s, respectively. We classified areas with woody canopy cover below 5 % as barren if grass biomass was below 100 kg/ha and as grassland if grass biomass exceeded 100 kg/ha. Grassland was classified as C3 grassland or C4 grassland based on predominance of C3 or C4 grass biomass. Simulated woody individuals were classified as trees if they had three or less stems and as shrubs if they had four or more stems (see Supplement). The canopy cover of woody plants and grass biomass was used to separate woodland and savanna biomes. Grid cells with tree canopy cover greater than shrub canopy cover, tree canopy cover between 5 % and 45 %, and grass biomass below 100 kg/ha were classified as woodland. Grid cells with the same woody cover characteristics but grass biomass higher than 100 kg/ha were classified as savanna. Savanna was further separated into C3 savanna and C4 savanna based on the predominance of C3 or C4 grass biomass. Areas with canopy cover greater than 45 % were classified as forest if tree cover exceeded shrub cover or shrubland if shrub cover exceeded tree cover, irrespective of grass biomass. Forests were subdivided into evergreen and deciduous forest based on the dominance of canopy area of both tree phenology types. In aDGVM2, whether a plant is deciduous or evergreen is decided by a trait. Biomes considered in this study were hence C3 grassland, C4 grassland, shrubland, woodland, deciduous forest, evergreen forest, C3 savanna and C4 savanna.

Biomes differ in the amount of precipitation they receive and their temperatures. Whittaker plots describe the boundaries of observed biomes with respect to temperature and precipitation. We used the “plotbiomes” R package ( by to create Whittaker plots based on Ricklefs (2008) and Whittaker (1975). We overlaid the simulated biomes on Whittaker plots to assess climatic niches of biomes under current climate to determine shifts in climatic niches by the end of this century as a result of climate change and elevated CO2 under both RCPs (see Sect. 3.6).

2.7 Calculation of biome-level evapotranspiration

For analyzing evapotranspiration change we calculated the amount of water transpired per unit leaf biomass. Simulated ET and leaf biomass for woody plants, C3 grass, and C4 grass were summed and scaled to the grid level, taking latitudinal variation of grid cell area into account. Absolute change in evapotranspiration quantity can either result from the change in biome area, from a change in total amount of leaf biomass over time or from changes in water use efficiency. In order to eliminate the effects caused by change in biome area and leaf biomass, we calculated biome-level evapotranspiration by normalizing evapotranspiration with biome-level leaf biomass (Eq. 6). Due to the normalization, differences in evapotranspiration at the biome level are comparable between different biomes and independent from biome attributes such as spatial extent and biome-level biomass. The biome-level evapotranspiration is calculated as the ratio of total annual ET over total leaf biomass for all respective biomes:

(6) E biome = i = 1 G ( E grid , i A grid , i ) i = 1 G ( B grid , i A grid , i ) ,

where Ebiome is biome-level ET (mm/kg/year), 1, 2, …, G represents the grid cells of the biome; Agrid,i is the area of grid cell i (m2); Egrid,i is the total evapotranspiration of grid cell i (mm/year); and Bgrid,i is the leaf biomass of grid cell i (kg/m2). Choosing to normalize evapotranspiration to leaf biomass integrates over both increased water use efficiency and soil water availability constraints. It is therefore suitable to characterize overall change in the water balance over time at the biome level, as it not only indicates water used to produce new biomass (as gross primary production over transpiration would express), but also includes water required to sustain existing biomass. We calculated the percentage change in Ebiome for respective scenarios between the 2010s and 2050s and between the 2010s and the 2090s.

Figure 2Comparison between simulated and observed biome patterns. (a) Simulated dominant biome type, (b) Sankey diagram showing overlap between simulated biomes and potential natural vegetation cover (ISLSCP II, Ramankutty et al., 2010), and (c) potential natural vegetation. The Sankey graph shows how aDGVM2 biomes and PNV classes overlap.

3 Results

3.1 Model performance and contemporary vegetation patterns

The aDGVM2 captured contemporary large-scale patterns of biomass, canopy cover, tree height and evapotranspiration. Model results agreed well with remote sensing products used for benchmarking (Figs. 1 and 2). R2 was 0.61, 0.45, 0.6 and 0.71, and NMSE was 0.48, 0.78, 0.4 and 1.07 for biomass (Saatchi et al.2011), tree height (Simard et al.2011), tree cover (Friedl et al.2010) and evapotranspiration (Zhang et al.2010), respectively (Fig. 1). Data–model agreement improved when masking out managed land (Tuanmu and Jetz2014). R2 increased to 0.66, 0.71, 0.67 and 0.80, while NMSE decreased to 0.43, 0.30, 0.61 and 1.03 for biomass, tree height, tree cover and evapotranspiration, respectively (Fig. S4). The model performed well in areas with higher fractional cover of natural vegetation, such as the Himalayas, Western Ghats and the northeast of the region, although the model overestimated biomass and canopy area in the Brahmaputra basin, which lies between 28–34 N and 90–96.5 E in the northeast of the study region (Fig. 1a, c,  Kumar et al.2020).

The model simulated evergreen forests along the Himalayan mountains, the southern part of the Western Ghats and Sri Lanka, whereas deciduous forest was simulated in the northern Western Ghats, central India and the southern parts of Myanmar (Fig. 2a). Savanna was simulated in the southern, northern, and western parts of India and some regions of central Myanmar. Shrublands were simulated in the arid regions of Pakistan, the western parts of India and Afghanistan. The aDGVM2 simulated woodland in the west of central India and grassland in the drier regions (Fig. 2a). A large proportion of simulated deciduous forest area is in good agreement with that in maps of potential natural vegetation (PNV, Fig. 2b, c). However a large proportion of simulated savanna area is represented as deciduous forest in the map of PNV (Fig. 2b).

3.2 Projected changes in biome distribution pattern

The aDGVM2 projected increasing trends for canopy cover and aboveground biomass in response to climate change and CO2 and hence changes in biome type, predominantly from savanna and grassland to deciduous forest (Fig. 3a, b). Simulations showed an increase in the area covered by evergreen and deciduous forests under both scenarios with eCO2, in contrast to simulations under both scenarios with fCO2, where CO2 was fixed after 2005 (Table 1). Under RCP4.5+eCO2, evergreen and deciduous forest cover increased by 3.1 % and 21.2 % until the 2050s and by 38.0 % and 59.1 % until the 2090s, respectively. Under RCP8.5+eCO2, evergreen and deciduous forest increased by 24.8 % and 45.4 % until the 2050s and by 46.5 % and 60.2 % until the 2090s, respectively. The model simulated a small increase in forest area for RCP4.5+fCO2, where the area increased by 7.9 % and 14.4 % for evergreen and deciduous forest until the 2090s, respectively. Evergreen forests were mainly simulated along the Himalayas, Western Ghats and eastern parts of the study region under current conditions (2000s, Fig. 3a) but expanded into the south of peninsular India in future periods (2050s and 2090s) under RCP4.5+eCO2. Deciduous forest cover also increased in future periods in central India and along the Himalayas (Figs. 3 and S6).

Table 1Biome cover (in %) for the 2000s, 2050s and 2090s, and percent (%) change in biome cover from the 2000s to the 2050s and the 2000s to the 2090s under RCP4.5 and RCP8.5 with fixed and elevated CO2. Δ indicates percent change in biome cover between time periods.

Download Print Version | Download XLSX

Figure 3Simulated biome distribution for the 2000s, 2050s and 2090s under (a) RCP4.5+eCO2 and (c) RCP4.5+fCO2. The Sankey diagrams show the fractional cover of biomes and transitions between biomes from the 2000s to the 2050s and the 2050s to the 2090s under (b) RCP4.5+eCO2 and (d) RCP4.5+fCO2. See Fig. S6 for the simulated biome distribution under RCP8.5.

The extent of C4 savanna showed a significant decrease under scenarios with eCO2, although in RCP4.5+eCO2 it showed an increase by 12.1 % between the 2010s and the 2050s (Table 1, Fig. 3). Simulated C4 savanna area decreased by 14.1 % relative to the 2000s until the 2090s under RCP4.5+eCO2. Under RCP8.5+eCO2 the model projected a decrease in C4 savanna area of 21.6 % and 32.2 % until the 2050s and the 2090s, respectively. The area covered by C4 savanna increased under both RCPs with fCO2 (Table 1). C4 savannas were mainly located in the northern plain and peninsular India in the baseline period. However, these areas were replaced by deciduous forests in the northern plain and central India and by evergreen forests in peninsular India and in the southeast of the region by the 2090s under eCO2 scenarios (Figs. 3a and S6a). The model simulated a decrease in area covered by woodland, shrubland, grasslands and C3 savanna by the 2090s under all scenarios (Table 1, Fig. 3). Simulations showed an increase in barren areas in the western part of the region under all scenarios (Figs. 3 and S6, Table 1).

3.3 Projected changes in biomass at the biome level

The aDGVM2 predicted an increase in mean biomass for evergreen and deciduous forest in the eCO2 scenarios for both RCPs (Table 2). Under RCP4.5+eCO2, mean aboveground biomass in evergreen and deciduous forest increased by 8.1 % and 14.4 % by the 2050s and 3.8 % and 15.7 % by the 2090s, relative to the baseline period. The increase is even higher under RCP8.5+eCO2 (Table 2). The mean biomass of woodland decreased under both RCPs except for the 2050s with eCO2 scenarios. The mean biomass of grassland increased under RCP4.5 but decreased for C4 grassland under RCP8.5 for both fCO2 and eCO2 scenarios. Shrublands in the western part of the study region showed an increase in mean biomass under eCO2 scenarios except for the 2050s under both RCPs and a decrease under fCO2 for both RCPs (Table 2). Our results showed that under RCP4.5 and RCP8.5 biomass decreased in the areas along the Himalayas, as well as in the central, northeastern and western parts of the study region by the end of the century. Modeled biomass decrease is higher under RCP8.5 in these regions (Figs. 4 and S7). Biomass in the central and southeastern part of the region is projected to increase under both RCPs with eCO2 until the 2050s and 2090s and to decrease in southern India and in parts of western South Asia (Figs. 4 and S7). We found increased biomass in Afghanistan, western Pakistan, Nepal, and the southern part of Myanmar and decreased biomass in the western arid part of the study region under both RCPs for both eCO2 and fCO2 (Fig. 5), though the magnitude of change is different (Figs. 4 and S7). There were few areas in the western part of the study region where the model predicted increased biomass only under fCO2 for both RCPs (Fig. 5). In large parts of the study region, biomass increased under eCO2 for both RCPs but decreased under fCO2; that is, CO2 fertilization compensates for climate-change-induced biomass diebacks in these regions (Fig. 5).

Table 2Mean biomass (in t/ha) within biomes for the 2000s, 2050s and 2090s, and percent (%) change in biomass from the 2000s to the 2050s and the 2000s to the 2090s under RCP4.5 and RCP8.5 with fixed and elevated CO2. Δ indicates percentual biomass changes between time periods.

Download Print Version | Download XLSX

Figure 4Projected change in biomass (t/ha), canopy cover (%) and evapotranspiration (ET, mm/year) between the 2000s and 2050s and between the 2000s and the 2090s under (a) RCP4.5+eCO2 and (b) RCP4.5+fCO2. See Fig. S7 for projected change of these variables under RCP8.5.

Figure 5Maps showing areas where CO2 fertilization compensates for biomass dieback caused by climate change between the 2000s and the 2090s under (a) RCP4.5 and (b) RCP8.5 and (c) aboveground biomass between 1950 and 2099 for South Asia in different scenarios.

3.4 Projected changes in evapotranspiration at the biome level

The response of simulated Ebiome varies in different biomes under both RCP4.5 and RCP8.5 (Table 3). Under the RCP4.5+fCO2 scenario the model predicted a decrease in ET in all biomes except for deciduous forest and shrubland where it increased by 1 % and 2.1 % until the 2050s and by 0.3 % and 11.9 % by the 2090s, respectively. Simulated Ebiome under RCP8.5+fCO2 for deciduous forest and shrubland increased by 4.2 % and 5.2 % until the 2050s and by 5.2 % and 16.4 % until the 2090s, respectively. The model also predicted increased Ebiome for C4 grassland, evergreen forest and C4 savanna until the 2090s under RCP8.5+fCO2 (Table 3). Comparisons of the RCP4.5+fCO2 and RCP8.5+fCO2 scenarios indicated that the former had a higher Ebiome than the latter scenario across all biomes because precipitation decrease is higher in the RCP8.5 scenario than in the RCP4.5 scenario. Under both RCPs with eCO2, the model predicted a decrease in Ebiome across all biomes, except for a marginal increase in shrubland under RCP4.5 and deciduous forest under RCP8.5 until the 2050s and the 2090s (Table 3). In general, scenarios with eCO2 showed lower biome-specific evapotranspiration across most of the biomes compared to simulations with fCO2.

Table 3Biome-level ET normalized to biomass (Ebiome, mm/kg/year) for the 2000s, 2050s and 2090s, and percent (%) change in Ebiome from the 2000s to the 2050s and the 2000s to the 2090s under RCP4.5 and RCP8.5 with fixed and elevated CO2. Δ indicates percentual ET changes between time periods.

Download Print Version | Download XLSX

3.5 Response of mean ET and mean aboveground biomass to climate change

The model predicted a larger increase in absolute annual mean ET (mm/year) under eCO2 than fCO2 for both RCP scenarios due to the corresponding increase in biomass (Figs. 4 and S7). We compared the spatially averaged annual values over all of South Asia of simulated absolute ET with MAP over the period from 1951 to 2099 and found a statically significant relation (p value < 0.005). We found that absolute ET was positively correlated with MAP under all four scenarios (Figs. 6a and S8a) but weakly correlated with MAT (Figs. 6b and S8b). For a given MAP, the spatially averaged annual value of aboveground biomass (AGBM) was lower in scenarios with fCO2 than scenarios with eCO2 under both RCPs (Figs. 6c and S8c). The spatially averaged annual value of AGBM decreased beyond a MAT of 23.5 C for both RCPs with fCO2, whereas it increased beyond 23.5 C under both RCP scenarios with eCO2 (Figs. 6d and S8d).

Figure 6Relationship between (a) evapotranspiration (ET) and mean annual precipitation (MAP), (b) ET and mean annual temperature (MAT), (c) mean aboveground biomass and MAP, and (d) mean aboveground biomass and MAT under RCP4.5. The lines (both solid and dotted) in all figures represent the best-fit regression line. Data points represent spatially averaged ET (a, b) and biomass (c, d) over all of South Asia for each year from 1950 to 2099. See Fig. S8 for results under RCP8.5.


3.6 Impact of climate change on climatic niches of biomes

The climate niches of simulated biomes broadly overlapped with the biome niches in the Whittaker scheme (Figs. 7 and S9, Ricklefs2008; based on Whittaker1975). Under RCP4.5+eCO2 and RCP8.5+eCO2, the aDGVM2 simulated shifts of climatic niches of biomes. Evergreen and deciduous forest biomes were predicted to invade the niche space of savannas under eCO2 scenarios (Figs. 7 and S9). Savannas in turn were predicted to expand their climatic niche to MAT > 30 C by 2099, a climatic space that was essentially not occupied by current biomes. Under RCP8.5+eCO2 in the 2090s, forests completely occupied the climate space, which is currently occupied by savanna (Fig. S9).

Figure 7Simulated climate niches of biomes for the (a) 2000s, (b) 2050s and (c) 2090s under RCP4.5+eCO2 and (d) 2000s, (e) 2050s and (f) 2090s under RCP4.5+fCO2. The simulated biomes are overlaid on the climate envelopes of Whittaker's biomes and are plotted following Ricklefs (2008) and Whittaker (1975). See Fig. S9 for projected change in climatic niches of biomes under RCP8.5.


In both scenarios with fCO2, savanna occupied the climate space delineated by MAT > 25 C and MAP between 500 and 1500 mm and did not experience major replacement by forest. The model predicted that savanna expansion in climate space was higher under RCP8.5+fCO2 than under RCP4.5+fCO2 (Figs. 7 and S9). Other biomes also experienced shifts in their climate space (Fig. 7). However, the results showed that for both the current and future period, grasslands and shrublands occupied the region with low MAP (<500 mm), and woodland occupied low MAP (<800 mm) regions, which correspond to the western arid and semiarid regions of the study region under the scenario with eCO2 (Fig. 7).

4 Discussion

4.1 Impact of climate change and elevated CO2 on biomes and biomass

Our simulations for RCP4.5+eCO2 and RCP8.5+eCO2 showed a strong positive response of vegetation growth, i.e., increases in biomass, canopy cover and canopy height. Mean biomass in most biomes was projected to increase, but the magnitude of increase differed considerably between different scenarios (Table 2). Projected change in canopy cover resulted in biome transitions. Under future conditions, the spatial distribution, extent and biomass of evergreen forests mostly remained at the current state, and evergreen forests were more resistant to climate change than deciduous forests. Expansion of deciduous forest into open biomes due to increasing woody cover resulted in significant loss of savanna area in the Deccan region under both RCPs with eCO2 by the end of the century. Transition from deciduous forests to evergreen forest was simulated for the mountain regions of South Asia (Scheiter et al.2020), i.e., the Himalayas and the Western Ghats, where precipitation was predicted to increase. The trade-offs between specific leaf area (ASLA) and leaf longevity (LL) result in the emergence of evergreen behavior in wet regions of South Asia. In the wet tropics, higher LL allows achieving a constant positive carbon balance from photosynthesis throughout the year and increases the residence time of nutrients and carbon in the plant and therefore enhances the photosynthetic gain per unit carbon and nutrient investment in leaves (Kikuzawa and Lechowicz2011). The deciduous behavior is advantageous in dry regions, as in the Deccan region, because trees that do not invest much carbon into their leaves per unit dry mass (higher ASLA and lower LL) lose less investment when shedding them during the dry season. Phenology change as a result of climate change has already been observed (Buitenwerf et al.2015; Cleland et al.2007). In Scheiter et al. (2020), we showed that climate change supports transitions to tall evergreen vegetation in tropical Asia and found increases in the abundance of evergreen plants and decreases in the abundance of deciduous plants in mainland Southeast Asia, central India and Pakistan. This relative advantage of evergreen plants over deciduous plants under elevated CO2 in aDGVM2 can be explained by the fact that increased intrinsic water use efficiency under eCO2 in evergreen plants is higher than in deciduous plants as demonstrated by Soh et al. (2019). Previous modeling studies also support aDGVM2 result showing transitions from deciduous to evergreen vegetation. With the BIOME4 model, Ravindranath et al. (2006) simulated the response of forest to Special Report on Emissions Scenarios (SRES) A2 and B2 and reported similar changes toward evergreen phenology. A study by Chaturvedi et al. (2011) using the IBIS model also predicted transitions toward evergreen forest.

Woody encroachment in many ecosystems is attributed to rising CO2, and this is supported by studies based on both field observations (e.g., FACE experiments) and satellite data (Brienen et al.2015; Fischlin et al.2007; Piao et al.2006; Schimel et al.2015; Archer et al.2017; Stevens et al.2017). The aDGVM2 also supports these findings, i.e., increasing canopy cover and woody biomass under the eCO2 condition, and agrees with the reported greening trend in South Asia during the last three decades (Wang et al.2017). Elevated CO2 affects plants by increasing their photosynthetic rate, growth rate and water use efficiency, leading to an increase in biomass (Leakey et al.2009; Norby and Zak2011). Increased photosynthetic rates under elevated CO2 are due to an increase in the rate of rubisco carboxylation for C3 plants, with a concurrent decrease in photorespiratory losses of carbon (Long et al.2004). Due to the improved carboxylation efficiency, C3 plants can respond by reducing stomatal conductance, thereby reducing transpirational losses, improving leaf water status and water use efficiency, and favoring leaf area growth (Long et al.2004; Norby and Zak2011). Evidence from both observation and modeling of forest dynamic suggests that combined effects of eCO2 and increased water use efficiency include increases in forest growth and canopy greening, as well as widespread increases in woody plant biomass and potential forest carbon sink. However, it is still unclear how the CO2 responses scale to the ecosystem level (Hickler et al.2015) and how nutrient limitation from the soil may influence ecosystem responses to eCO2. Körner (2015) argued that carbon from the atmosphere can only be converted into biomass if other factors such as nutrients, temperature and water are not limiting. In addition, the benefit of eCO2 can be downregulated by broad-scale forest die-off due to frequent drought and warmer temperature (Choat et al.2018; Mcdowell et al.2016), as well as tree mortality due to negative tree physiological responses, negative carbon balance and accelerated pest attacks. Rising background mortality rates combined with projected increases in intensity, frequency and duration of drought (Huang et al.2016) increase the uncertainty regarding positive effects of eCO2.

In the long run, whether ecosystems act as a carbon source or sink can be estimated using models that consider all factors that are relevant in the carbon cycle and its associated factors (Fatichi et al.2014; Körner2015). However, Terrer et al. (2019) showed that the global-scale response to eCO2 from experiments is similar to past changes in greenness (Piao et al.2019) and biomass (Sitch et al.2015) in response to eCO2. This suggests that CO2 will likely continue to stimulate plant biomass in the future despite the constraining effect of soil nutrients; however Terrer et al. (2019) also argued that the empirical relationships with soil nutrients can be powerful for explaining large-scale patterns of eCO2 responses, despite ecosystem-level uncertainties. According to our simulations we can conclude that natural vegetation of South Asia likely will remain a carbon sink in the future (Fig. 5).

4.2 Impact of climate change and fixed CO2 on biomes and biomass

Under both fCO2 scenarios, the spatial distribution of savanna areas remained in its contemporary state. Central India and the Deccan Plateau showed a transition of deciduous forest to savanna, because forest canopy opened up due to tree mortality caused by increasing temperature and reduced MAP. This indicates that plants experience temperature and drought stress under fixed CO2. These stresses were compensated for by CO2 fertilization in eCO2 scenarios where the aDGVM2 simulated increased biomass and woody encroachment in areas affected by climate-induced dieback in fCO2 simulations. This aDGVM2 behavior agrees with results presented by Lapola et al. (2009), who modeled biome shifts from forest to savanna in the absence of CO2 fertilization for the Amazon region. Changes in precipitation regimes are likely to have a strong influence particularly in arid and semi‐arid regions, such as grasslands (Verstraete et al.2009). The complex interactions of inter-annual precipitation variability and precipitation seasonality can result in rapid ecosystem transitions (e.g., between alternative stable states with high and low vegetation biomass; Holmgren and Scheffer2001). The decrease in simulated AGBM after MAT increases beyond 23.5 C under fCO2 scenarios can be explained by the longer exposure of vegetation to temperatures beyond the optimum temperature range of C3 photosynthesis during the main growing season. This effect was further enhanced by decreasing MAP and the absence of CO2 fertilization. This implies that the increase in MAT above 23.5 C together with weak CO2 fertilization would have negative consequences for carbon sequestration. The sensitivity of biomass to temperature and CO2 change has been investigated in many studies (Norby and Luo2004; Song et al.2019; Sperry et al.2019). A meta-analysis by Lin et al. (2010) showed that warming significantly increased biomass by 12.3 % (with a 95 % confidence interval of 8.4 %–16.3 %) across all the terrestrial plants included. This observation is consistent with our model results. Biomass showed a positive relation with MAT, which did not change with mean annual precipitation or experimental duration or CO2 enrichment (Lin et al.2010). These findings are also supported by previous studies by Rustad et al. (2001), Walker et al. (2006), and Dormann and Woodin (2002), which have revealed that warming generally increases terrestrial plant biomass, indicating enhanced terrestrial carbon uptake via plant growth. Previous modeling studies using Biome‐BGC (Running and Hunt Jr1993), Century (Parton et al.1993) and TEM (Tian et al.1999) have shown an increase in productivity when both climate change and CO2 effects were considered. However, the increase was smaller when only climate change effects were considered, and both Biome‐BGC and TEM suggest that without CO2 fertilization, average productivity would decline relative to the current annual average as shown by our result (Fig. 6d). This suggests complexity and challenges in seeking general patterns of terrestrial plant growth in a future warmer climate condition. It also implies that we need a better understanding of impacts of heat stress on vegetation and how it interacts with drought and CO2 fertilization. It is also unclear to what degree thermal acclimation may counteract some of the negative effects on plant growth caused by higher temperatures (Lombardozzi et al.2015).

4.3 Impact of climate change and CO2 change on climatic niches of biomes

Elevated CO2 has a major impact on the climatic niche space of biomes. Our simulations showed forest invasion into the niche space currently occupied by savanna by the end of the century. The expansion of forests to drier areas corresponds to a widening of their climate niche space under eCO2. This expansion is mainly driven by eCO2 and is corroborated by the fact that in the absence of CO2 fertilization the climatic niche of biomes is stable; i.e., biomes occupy the same niche space under current and future conditions. These findings imply that the bioclimatic boundaries used to define biome niche space are not static but are specific for given CO2 levels. Therefore, the thresholds of Whittaker's biomes need to be redefined for a high-CO2 world such that they encompass the altered climatic envelopes of biomes under elevated CO2 in the future (Fig. 7). The shift in niche space can be attributed to the shift in plant communities caused by the combined effect of climate change and elevated CO2, which increases plant water use efficiency, allowing them to grow under drier conditions. These community shifts can also lead to a change in the characteristics of biomes by altering community structure and ecosystem functions (Chapin et al.1997).

4.4 Effect of CO2 on ET and its interaction with climate change

Climate change has direct effects on hydrological processes (Liu et al.2008). ET and water deficit influence plant productivity and distribution (Stephenson1998). Higher biomass coincided with increased absolute amounts of ET for eCO2 scenarios in some parts of the study region under both RCPs by the end 21st century (Figs. 5 and S7). This change can be attributed to higher leaf biomass accumulated in plants enabled by increased photosynthetic efficiency under eCO2. The higher amount of leaf biomass offsets the water-saving effect caused by reduced stomatal conductance due to improved water use efficiency under eCO2 scenarios and resulted in reduced ET per unit leaf biomass (Warren et al.2011). Our results showed that the strength of the CO2 fertilization effect is relevant when attempting to determine Ebiome at the biome level during the 21st century. Biome-specific ET decrease was less pronounced under RCP4.5 due to a lower concentration of atmospheric CO2 compared to RCP8.5. Our simulated decrease in ET in response to climate change and increasing CO2 concentration agrees with Kergoat et al. (2002), who have reported decreased ET under elevated CO2 concentration in a chamber experiment. However, reduced ET under eCO2 can reduce regional-scale atmospheric humidity and thereby enhance the vapor pressure deficit (VPD) between leaves and atmosphere, a driving force for water loss, which may partially counteract CO2-induced reduction of ET due to decreased stomatal conductance. Due to stomatal closure, photosynthetic rates under soil water stress conditions decline in aDGVM2 when atmospheric VPD increases. The projected increase in air temperature also increases the saturated water vapor pressure. As a result VPD will increase, given that increase in actual vapor pressure is limited by soil water availability while the increase in saturated vapor pressure is not (Yuan et al.2019), and potential evapotranspiration will increase with temperature (Warren et al.2011). As future climate projections vary spatially and temporally, there was high model uncertainty on how ET will respond to changes in precipitation and temperature.

4.5 Implication of the projected change for conservation

Changes in biome types imply changes in biodiversity, ecosystem function and productivity. Each biome is characterized by a range of distinctive ecological processes and functions. For instance, the distribution of forest ecosystem in the mountains is largely regulated by the altitude and climatic factors (Saikia et al.2017). They have high species richness and needed to be protected from the ever-increasing anthropogenic pressure and climate change. Open biomes such as grassland and savanna support high biodiversity (Parr et al.2014). Pronounced increases in tree density in grasslands and savannas will alter vegetation structure and reduce grassland biodiversity. Such changes will negatively affect savanna-specific ecosystem services such as grazing potential, tourism and wildlife habitat availability (Parr et al.2012). The threat posed to the biodiversity of Asian savannas by climate change is aggravated by inadequate management policies that misinterpret them as degraded forest (Ratnam et al.2016). In this context, management policies aim to afforest open biomes, although paleoecological evidence indicates that these open biomes are not degraded forest but ancient ecosystems (Kumar et al.2020; Ratnam et al.2016). Moreover, increased woody cover can negatively affect water resources in the semiarid regions of the study area. Acharya et al. (2018) have shown that increased woody cover hinders the downward movement of water, causing increased water inception, which has negative effects on groundwater recharge. It is therefore necessary to control the abundance of woody plants in semiarid regions to control stream flow and enhance groundwater recharge (Bednarz et al.2001).

In South Asia, biodiversity hotspots have a very unique topography, where climate varies strongly over short distances. As global biodiversity hotspots, mountain forest ecosystems in the Western Ghats, the Himalayas and northeastern part of the study area (India and Myanmar) are particularly vulnerable to climate change (Myers et al.2000) and need targeted management action to mitigate adverse effects. Conservation of these hotspots requires consideration of many different attributes of plant communities, ecosystems, landscapes, and plant diversity; of how they will change; and of how their ecosystem services are valued.

Conservation methods and policies that can accommodate minimal losses of ecosystem services and provide robust strategies to mitigate climate change impacts should be developed and implemented. In this context, DGVMs facilitate the exploration of vegetation–climate interactions by providing detailed results for different management and climate scenarios. Such an exploration of different possible scenarios is necessary to develop optimized mitigation and conservation strategies for protected areas in biodiversity hotspots. The value of DGVM modeling results lies in their potential to provide insights into multiple future trajectories. Based on the most likely trajectories, the results can be used to tailor best-practice strategies for decision makers that need to manage conservation areas or protected areas (Boulangeat et al.2012).

4.6 Limitations of this modeling study

Our simulation results are constrained by the model formulation and the assumptions underlying aDGVM2. Disagreement between model results and data used for benchmarking can be attributed to the fact that the aDGVM2 simulates potential natural vegetation, whereas remote sensing products integrate land use. This implies that enhancing the model to simulate observed land cover patterns would require additional information on anthropogenic impacts. Anthropogenic activities such as deforestation, habitat conversion and urbanization can modify the interactions between climate, plant communities and biomes (Hansen et al.2001).

In addition data–model disagreement can be explained by model uncertainties and processes currently not considered in aDGVM2. For instance, aDGVM2 uses carbon allocation parameters that are not easily measurable in the field, which limits the evaluation of simulated mechanisms. The model currently does not account for carbon that plants invest into nutrient acquisition (e.g., mycorrhiza) or into defenses against predation and pathogens (Zemunik et al.2015). There is insufficient ecophysiological data from the study region, which are required for parameterization of trait ranges used to simulate regional plant communities (Kumar and Scheiter2019). The complexity of the interactions between global change and biomes as well as biodiversity is difficult to model in the absence of such data. While the model currently captures effects related to CO2 fertilization and temperature, associated mortality reasons such as pests attack and heat damage to plant tissues are insufficiently represented in the models. The low resolution of input data, both soil and climate data, also limits the model's capability to capture high-resolution regional heterogeneity in vegetation distribution. Further, the strength of CO2 fertilization modeled in aDGVM2 may be overestimated because the effect of nutrient limitation on productivity is not included in this version of aDGVM2 (Körner et al.2005; Terrer et al.2018). Despite these caveats, we are nonetheless confident in capturing general patterns of future global change and its consequences for biomes and their boundaries in South Asia.

5 Conclusions

The model reproduced the contemporary distribution of biomes, biomass, evapotranspiration and tree height. We investigated the impact of eCO2 and climate change on South Asian biomes and found that climate change and CO2 fertilization in combination are substantial drivers of biome change and that elevated CO2 concentrations altered the climatic envelope of biomes in addition to causing increases in biomass, tree height and canopy cover. Continued biomass increase indicates that South Asia's natural vegetation will likely remain a carbon sink in the 21st century. Our results also imply that woody encroachment poses a threat to open biomes and causes the transition of savanna biomes to deciduous forest in the future. The savanna biome is important in the context of biodiversity conservation. We showed that bioclimatic niches of biomes are not static but are specific for given CO2 concentrations. We therefore argue that Whittaker plots used to illustrate niches of biomes need to be adjusted for future climate conditions. We also found that the simulated decrease in biomass-specific ET is more pronounced in scenarios with eCO2 than in scenarios with fCO2, which indicates that water use efficiency will likely increase due to CO2 fertilization.

The biome transitions simulated under eCO2 and changing climate indicate the need to adjust ecosystem management, mitigation strategies and conservation policies for protected areas to allow targeted long-term management. To understand the significance of ecological responses to climate change, it is essential to improve and expand biological monitoring activities (Loreau et al.2001). To achieve this, the most vulnerable biomes that we identified in this study could be proposed as high-priority targets for programs that monitor vegetation–climate interactions, productivity and biodiversity (Proença et al.2017).

Code availability

The aDGVM2 code as well as scripts to conduct the model experiments and analyze the results are available upon request. Please contact any of the authors.


The supplement related to this article is available online at:

Author contributions

DK and SS conceived the study; DK included changes to aDGVM2 for the study region, conducted model simulations, analyzed results, created figures and led the writing. All authors contributed to the development of aDGVM2 and commented on the text.

Competing interests

The authors declare that they have no conflict of interest.


Dushyant Kumar and Simon Scheiter thank the Deutsche Forschungsgemeinschaft (DFG) for funding (Emmy Noether grant). Mirjam Pfeiffer acknowledges the German Federal Ministry of Education and Research (BMBF) for funding (SPACES II initiative, SALLnet project).

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. SCHE 1719/2-1) and the Bundesministerium für Bildung und Forschung (grant no. 01LL1802B).

Review statement

This paper was edited by Eyal Rotenberg and reviewed by Hisashi Sato and two anonymous referees.


Acharya, B. S., Kharel, G., Zou, C. B., Wilcox, B. P., and Halihan, T.: Woody plant encroachment impacts on groundwater recharge: A review, Water, 10, 1466,, 2018. a

Allen, M. R., Barros, V. R., Broome, J., Cramer, W., Christ, R., Church, J. A., Clarke, L., Dahe, Q., Dasgupta, P., and Dubash, N. K.: IPCC fifth assessment synthesis report-climate change 2014 synthesis report, Intergovernmental Panel on Climate Change, Geneva, Switzerland, 2014. a

Archer, S. R., Andersen, E. M., Predick, K. I., Schwinning, S., Steidl, R. J., and Woods, S. R.: Woody plant encroachment: causes and consequences, in: Rangeland systems, Springer, Cham, Switzerland, 25–84,, 2017. a

Bednarz, S. T., Dybala, T., Muttiah, R. S., Rosenthal, W., and Dugas, W. A.: Brush/water yield feasibility studies, Blackland Research Center, Temple, Texas, USA, 2001. a

Bera, S. K., Basumatary, S. K., Agarwal, A., and Ahmed, M.: Conversion of forest land in Garo Hills, Meghalaya for construction of roads: A threat to the environment and biodiversity, Curr. Sci. India, 91, 281–284, 2006. a

Boulangeat, I., Philippe, P., Abdulhak, S., Douzet, R., Garraud, L., Lavergne, S., Lavorel, S., Van Es, J., Vittoz, P., and Thuiller, W.: Improving plant functional groups for dynamic models of biodiversity: at the crossroads between functional and community ecology, Global Change Biol., 18, 3464–3475,, 2012. a

Brienen, R. J., Phillips, O. L., Feldpausch, T. R., Gloor, E., Baker, T. R., Lloyd, J., Lopez-Gonzalez, G., Monteagudo-Mendoza, A., Malhi, Y., and Lewis, S. L.: Long-term decline of the Amazon carbon sink, Nature, 519, 344–348, 2015. a

Brodie, J. F., Aslan, C. E., Rogers, H. S., Redford, K. H., Maron, J. L., Bronstein, J. L., and Groves, C. R.: Secondary extinctions of biodiversity, Trends Ecol. Evol., 29, 664–672,, 2014. a

Buitenwerf, R., Rose, L., and Higgins, S. I.: Three decades of multi-dimensional change in global leaf phenology, Nat. Clim. Change, 5, 364–368,, 2015. a

Cao, L., Bala, G., Caldeira, K., Nemani, R., and Ban-Weiss, G.: Importance of carbon dioxide physiological forcing to future climate change, P. Natl. Acad. Sci. USA, 107, 9513–9518,, 2010. a

Chapin, F. S., Walker, B. H., Hobbs, R. J., Hooper, D. U., Lawton, J. H., Sala, O. E., and Tilman, D.: Biotic control over the functioning of ecosystems, Science, 277, 500–504,, 1997. a

Chaturvedi, R. K., Gopalakrishnan, R., Jayaraman, M., Bala, G., Joshi, N. V., Sukumar, R., and Ravindranath, N. H.: Impact of climate change on Indian forests: a dynamic vegetation modeling approach, Mitig. Adapt. Strat. Gl., 16, 119–142,, 2011. a, b

Chen, I.-C., Hill, J. K., Ohlemüller, R., Roy, D. B., and Thomas, C. D.: Rapid range shifts of species associated with high levels of climate warming, Science, 333, 1024–1026,, 2011. a

Choat, B., Brodribb, T. J., Brodersen, C. R., Duursma, R. A., López, R., and Medlyn, B. E.: Triggers of tree mortality under drought, Nature, 558, 531–539,, 2018. a

Choudhury, B. J., DiGirolamo, N. E., Susskind, J., Darnell, W. L., Gupta, S. K., and Asrar, G.: A biophysical process-based estimate of global land surface evaporation using satellite and ancillary data II, Regional and global patterns of seasonal and annual variations, J. Hydrol., 205, 186–204,, 1998. a

Cleland, E. E., Chuine, I., Menzel, A., Mooney, H. A., and Schwartz, M. D.: Shifting plant phenology in response to global change, Trends Ecol. Evol., 22, 357–365,, 2007. a

Collatz, G. J., Ball, J. T., Grivet, C., and Berry, J. A.: Physiological and environmental regulation of stomatal conductance, photosynthesis and transpiration: a model that includes a laminar boundary layer, Agr. Forest Meteorol., 54, 107–136,, 1991. a

Collatz, G. J., Ribas-Carbo, M., and Berry, J. A.: Coupled photosynthesis-stomatal conductance model for leaves of C4 plants, Funct. Plant Biol., 19, 519–538,, 1992. a, b, c

Deb, J., Phinn, S. R., Butt, N., and McAlpine, C. A.: Summary of climate change impacts on tree species distribution, phenology, forest structure and composition for each of the 85 studies reviewed, The University of Queensland [Data Collection],, 2017. a

Doherty, R. M., Sitch, S., Smith, B., Lewis, S. L., and Thornton, P. K.: Implications of future climate and atmospheric CO2 content for regional biogeochemistry, biogeography and ecosystem services across East Africa, Global Change Biol., 16, 617–640,, 2010. a

Dormann, C. and Woodin, S. J.: Climate change in the Arctic: using plant functional types in a meta-analysis of field experiments, Funct. Ecol., 16, 4–17,, 2002. a

Eckstein, D., Hutfils, M., and Winges, M.: Global climate risk index 2019: Who suffers most from extreme weather events? Weather-related loss events in 2017 and 1998 to 2017, Germanwatch, Bonn, Germany, 2018. a

Ehleringer, J. R. and Cerling, T. E.: C3 and C4 photosynthesis, Encyclopedia of Global Environmental Change, 2, 186–190, 2002. a

Farquhar, G. D., von Caemmerer, S. V., and Berry, J. A.: A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species, Planta, 149, 78–90,, 1980. a

Fatichi, S., Leuzinger, S., and Körner, C.: Moving beyond photosynthesis: from carbon source to sink-driven vegetation modeling, New Phytol., 201, 1086–1095,, 2014. a

Feng, H., Zou, B., and Luo, J.: Coverage-dependent amplifiers of vegetation change on global water cycle dynamics, J. Hydrol., 550, 220–229,, 2017. a

Field, C. B., Lobell, D. B., Peters, H. A., and Chiariello, N. R.: Feedbacks of terrestrial ecosystems to climate change, Annu. Rev. Env. Resour., 32, 1–29,, 2007. a

Fischlin, A., Midgley, G. F., Price, J. T., Leemans, R., Gopal, B., Turley, C. M., Rounsevell, M. D. A., Dube, P., Tarazona, J., and Velichko, A.: Impacts adaptation and vulnerability, in: Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Parry, M. L., Canziani, O. F., Palutikof, J. P., van der Linden, P. J., and Hanson, C. E., Cambridge University Press, Cambridge, UK, 391–431, 2007. a

Fisher, J. B., Whittaker, R. J., and Malhi, Y.: ET come home: potential evapotranspiration in geographical ecology, Global Ecol. Biogeogr., 20, 1–18,, 2011. a

Forster, P., Ramaswamy, V., Artaxo, P., Berntsen, T., Betts, R., Fahey, D. W., Haywood, J., Lean, J., Lowe, D. C., Myhre, G., Nganga, J., Prinn, R., Raga, G., Schulz, M., and Van Dorland, R.: Changes in Atmospheric Constituents and in Radiative Forcing, in: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2007. a

Frank, D., Reichstein, M., Bahn, M., Thonicke, K., Frank, D., Mahecha, M. D., Smith, P., Van der Velde, M., Vicca, S., and Babst, F.: Effects of climate extremes on the terrestrial carbon cycle: concepts, processes and potential future impacts, Global Change Biol., 21, 2861–2880,, 2015. a

Friedl, M. A., Sulla-Menashe, D., Tan, B., Schneider, A., Ramankutty, N., Sibley, A., and Huang, X.: MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets, Remote Sens. Environ., 114, 168–182,, 2010. a, b

Friedlingstein, P., Cox, P., Betts, R., Bopp, L., von Bloh, W., Brovkin, V., Cadule, P., Doney, S., Eby, M., and Fung, I.: Climate–carbon cycle feedback analysis: results from the C4MIP model intercomparison, J. Climate, 19, 3337–3353,, 2006. a

Gaillard, C., Langan, L., Pfeiffer, M., Kumar, D., Martens, C., Higgins, S. I., and Scheiter, S.: African shrub distribution emerges via a trade‐off between height and sapwood conductivity, 45, 2815–2826,, 2018. a, b, c

Gallardo-Cruz, J. A., Pérez-García, E. A., and Meave, J. A.: β-Diversity and vegetation structure as influenced by slope aspect and altitude in a seasonally dry tropical landscape, Landscape Ecol., 24, 473–482,, 2009. a

Gates, D. M.: Transpiration and leaf temperature, Ann. Rev. Plant Physio., 19, 211–238,, 1968. a

Hansen, A. J., Neilson, R. P., Dale, V. H., Flather, C. H., Iverson, L. R., Currie, D. J., Shafer, S., Cook, R., and Bartlein, P. J.: Global change in forests: responses of species, communities, and biomes: interactions between climate change and land use are projected to cause large shifts in biodiversity, BioScience, 51, 765–779,[0765:GCIFRO]2.0.CO;2, 2001. a

Herring, S. C., Christidis, N., Hoell, A., Kossin, J. P., Schreck III, C. J., and Stott, P. A.: Explaining extreme events of 2016 from a climate perspective, B. Am. Meteorol. Soc., 99, 1–157,, 2018. a

Hickler, T., Prentice, I. C., Smith, B., Sykes, M. T., and Zaehle, S.: Implementing plant hydraulic architecture within the LPJ Dynamic Global Vegetation Model, Global Ecol. Biogeogr., 15, 567–577,, 2006. a

Hickler, T., Rammig, A., and Werner, C.: Modelling CO2 impacts on forest productivity, Current Forestry Reports, 1, 69–80,, 2015. a

Hijmans, R. J. and van Etten, J.: raster: Geographic analysis and modeling with raster data, R package version 2.0–12, 2012. a

Holmgren, M. and Scheffer, M.: El Niño as a window of opportunity for the restoration of degraded arid ecosystems, Ecosystems, 4, 151–159,, 2001. a

Huang, J., Yu, H., Guan, X., Wang, G., and Guo, R.: Accelerated dryland expansion under climate change, Nat. Clim. Change, 6, 166–171,, 2016. a

Jarvis, A., Reuter, H. I., Nelson, A., and Guevara, E.: Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90 m Database, 2008. a

Jucker, T., Bongalov, B., Burslem, D. F., Nilus, R., Dalponte, M., Lewis, S. L., Phillips, O. L., Qie, L., and Coomes, D. A.: Topography shapes the structure, composition and function of tropical forest landscapes, Ecol. Lett., 21, 989–1000,, 2018. a

Kattge, J. and Knorr, W.: Temperature acclimation in a biochemical model of photosynthesis: a reanalysis of data from 36 species, Plant Cell Environ., 30, 1176–1190,, 2007. a, b

Kergoat, L., Lafont, S., Douville, H., Berthelot, B., Dedieu, G., Planton, S., and Royer, J.-F.: Impact of doubled CO2 on global-scale leaf area index and evapotranspiration: Conflicting stomatal conductance and LAI responses, J. Geophys. Res.-Atmos., 107, 4808,, 2002. a

Kgope, B. S., Bond, W. J., and Midgley, G. F.: Growth responses of African savanna trees implicate atmospheric CO2 as a driver of past and current changes in savanna tree cover, Austral Ecol., 35, 451–463,, 2010. a

Kikuzawa, K. and Lechowicz, M. J.: Ecology of leaf longevity, in: Ecology of Leaf Longevity, 1–6, Springer, Tokyo, 2011. a

Kirschbaum, M. U. F.: Direct and indirect climate change effects on photosynthesis and transpiration, Plant Biology, 6, 242–253,, 2004. a

Körner, C.: Paradigm shift in plant growth control, Curr. Opin. Plant Biol., 25, 107–114,, 2015. a, b

Körner, C., Asshoff, R., Bignucolo, O., Hättenschwiler, S., Keel, S. G., Peláez-Riedl, S., Pepin, S., Siegwolf, R. T., and Zotz, G.: Carbon flux and growth in mature deciduous forest trees exposed to elevated CO2, Science, 309, 1360–1362,, 2005. a

Kumar, D. and Scheiter, S.: Biome diversity in South Asia – How can we improve vegetation models to understand global change impact at regional level?, Sci. Total Environ., 671, 1001–1016,, 2019. a, b

Kumar, D., Pfeiffer, M., Gaillard, C., Langan, L., Martens, C., and Scheiter, S.: Misinterpretation of Asian savannas as degraded forest can mislead management and conservation policy under climate change, Biol. Conserv., 241, 108293,, 2020. a, b, c

Langan, L., Higgins, S. I., and Scheiter, S.: Climate-biomes, pedo-biomes or pyro-biomes: which world view explains the tropical forest–savanna boundary in South America?, J. Biogeogr., 44, 2319–2330,, 2017. a, b, c, d, e, f, g, h, i

Lapola, D. M., Priess, J. A., and Bondeau, A.: Modeling the land requirements and potential productivity of sugarcane and jatropha in Brazil and India using the LPJmL dynamic global vegetation model, Biomass Bioenerg., 33, 1087–1095,, 2009. a

Leakey, A. D., Ainsworth, E. A., Bernacchi, C. J., Rogers, A., Long, S. P., and Ort, D. R.: Elevated CO2 effects on plant carbon, nitrogen, and water relations: six important lessons from FACE, J. Exp. Bot., 60, 2859–2876,, 2009. a, b, c

Lin, D., Xia, J., and Wan, S.: Climate warming and biomass accumulation of terrestrial plants: a meta-analysis, New Phytol., 188, 187–198,, 2010. a, b

Liu, M., Tian, H., Chen, G., Ren, W., Zhang, C., and Liu, J.: Effects of Land-Use and Land-Cover Change on Evapotranspiration and Water Yield in China During 1900–2000, J. Am. Water Resour. As., 44, 1193–1207,, 2008. a

Liu, M.-Z. and Osborne, C. P.: Leaf cold acclimation and freezing injury in C3 and C4 grasses of the Mongolian Plateau, J. Exp. Bot., 59, 4161–4170,, 2008. a

Lloyd, J., Bird, M. I., Vellen, L., Miranda, A. C., Veenendaal, E. M., Djagbletey, G., Miranda, H. S., Cook, G., and Farquhar, G. D.: Contributions of woody and herbaceous vegetation to tropical savanna ecosystem productivity: a quasi-global estimate, Tree Physiol., 28, 451–468,, 2008. a

Lombardozzi, D. L., Bonan, G. B., Smith, N. G., Dukes, J. S., and Fisher, R. A.: Temperature acclimation of photosynthesis and respiration: A key uncertainty in the carbon cycle-climate feedback, Geophys. Res. Lett., 42, 8624–8631,, 2015. a

Long, S. P., Ainsworth, E. A., Rogers, A., and Ort, D. R.: Rising atmospheric carbon dioxide: plants FACE the future, Annu. Rev. Plant Biol., 55, 591–628,, 2004. a, b

Loreau, M., Naeem, S., Inchausti, P., Bengtsson, J., Grime, J. P., Hector, A., Hooper, D. U., Huston, M. A., Raffaelli, D., and Schmid, B.: Biodiversity and ecosystem functioning: current knowledge and future challenges, Science, 294, 804–808,, 2001. a

Mao, J., Fu, W., Shi, X., Ricciuto, D. M., Fisher, J. B., Dickinson, R. E., Wei, Y., Shem, W., Piao, S., and Wang, K.: Disentangling climatic and anthropogenic controls on global terrestrial evapotranspiration trends, Environ. Res. Lett., 10, 094008,, 2015. a

Mcdowell, N. G., Williams, A., Xu, C., Pockman, W., Dickman, L., Sevanto, S., Pangle, R., Limousin, J., Plaut, J., Mackay, D. S., and Ogee, J.: Multi-scale predictions of massive conifer mortality due to chronic temperature rise, Nat. Clim. Change, 6, 295–300,, 2016. a

McSweeney, C. F. and Jones, R. G.: How representative is the spread of climate projections from the 5 CMIP5 GCMs used in ISI-MIP?, Climate Services, 1, 24–29,, 2016. a

Meinshausen, M., Smith, S. J., Calvin, K., Daniel, J. S., Kainuma, M. L. T., Lamarque, J.-F., Matsumoto, K., Montzka, S. A., Raper, S. C. B., and Riahi, K.: The RCP greenhouse gas concentrations and their extensions from 1765 to 2300, Climatic Change, 109, 213,, 2011. a, b

Myers, N., Mittermeier, R. A., Mittermeier, C. G., Da Fonseca, G. A., and Kent, J.: Biodiversity hotspots for conservation priorities, Nature, 403, 853–858,, 2000. a, b, c

Nachtergaele, F., van Velthuizen, H., Verelst, L., Batjes, N., Dijkshoorn, K., Van Engelen, V., Fischer, G., Jones, A., Montanarella, L., and Petri, M.: Harmonized world soil database (version 1.1), FAO, Rome, Italy, IIASA, Laxenburg, Austria, 2009. a

Nolan, C., Overpeck, J. T., Allen, J. R., Anderson, P. M., Betancourt, J. L., Binney, H. A., Brewer, S., Bush, M. B., Chase, B. M., and Cheddadi, R.: Past and future global transformation of terrestrial ecosystems under climate change, Science, 361, 920–923,, 2018. a, b

Norby, R. J. and Luo, Y.: Evaluating ecosystem responses to rising atmospheric CO2 and global warming in a multi-factor world, New Phytol., 162, 281–293,, 2004. a

Norby, R. J. and Zak, D. R.: Ecological lessons from free-air CO2 enrichment (FACE) experiments, Annu. Rev. Ecol. Evol. S., 42, 181–203,, 2011. a, b, c, d

Overpeck, J. T., Rind, D., and Goldberg, R.: Climate-induced changes in forest disturbance and vegetation, Nature, 343, 51–53,, 1990. a

Parr, C. L., Gray, E. F., and Bond, W. J.: Cascading biodiversity and functional consequences of a global change–induced biome switch, Divers. Distrib., 18, 493–503,, 2012. a

Parr, C. L., Lehmann, C. E., Bond, W. J., Hoffmann, W. A., and Andersen, A. N.: Tropical grassy biomes: misunderstood, neglected, and under threat, Trends Ecol. Evol., 29, 205–213,, 2014. a

Parry, M., Parry, M. L., Canziani, O., Palutikof, J., Van der Linden, P., and Hanson, C.: Climate change 2007 – impacts, adaptation and vulnerability: Working group II contribution to the fourth assessment report of the IPCC, Cambridge University Press, Cambridge, UK, 2007. a, b

Parton, W., Scurlock, J., Ojima, D., Gilmanov, T., Scholes, R., Schimel, D. S., Kirchner, T., Menaut, J.-C., Seastedt, T., Garcia Moya, E., and Kamnalrut, A.: Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide, Global Biogeochem. Cy., 7, 785–809,, 1993. a

Pfeiffer, M., Langan, L., Linstädter, A., Martens, C., Gaillard, C., Ruppert, J. C., Higgins, S. I., Mudongo, E. I., and Scheiter, S.: Grazing and aridity reduce perennial grass abundance in semi-arid rangelands – Insights from a trait-based dynamic vegetation model, Ecol. Model., 395, 11–22,, 2019. a, b

Piao, S., Friedlingstein, P., Ciais, P., Zhou, L., and Chen, A.: Effect of climate and CO2 changes on the greening of the Northern Hemisphere over the past two decades, Geophys. Res. Lett., 33, L23402,, 2006. a

Piao, S., Wang, X., Park, T., Chen, C., Lian, X., He, Y., Bjerke, J. W., Chen, A., Ciais, P., Tømmervik, H., and Nemani, R. R.: Characteristics, drivers and feedbacks of global greening, Nat. Rev. Earth Environ., 1, 14–27,, 2019. a

Prentice, I. C., Bondeau, A., Cramer, W., Harrison, S. P., Hickler, T., Lucht, W., Sitch, S., Smith, B., and Sykes, M. T.: Dynamic global vegetation modeling: quantifying terrestrial ecosystem responses to large-scale environmental change. In Terrestrial ecosystems in a changing world, 175–192, Springer, Berlin, Heidelberg, 2007. a, b

Proença, V., Martin, L. J., Pereira, H. M., Fernandez, M., McRae, L., Belnap, J., Böhm, M., Brummitt, N., García-Moreno, J., and Gregory, R. D.: Global biodiversity monitoring: from data sources to essential biodiversity variables, Biol. Conserv., 213, 256–263,, 2017. a

Ramankutty, N., Foley, J. A., Hall, F., Collatz, G., Meeson, B., Los, S., Brown De Colstoun, E., and Landis, D.: ISLSCP II Potential Natural Vegetation Cover, ORNL DAAC, Oak Ridge, Tennessee, USA,, 2010. a

Ratnam, J., Tomlinson, K. W., Rasquinha, D. N., and Sankaran, M.: Savannahs of Asia: antiquity, biogeography, and an uncertain future, Philos. T. Roy. Soc. B, 371, 20150305,, 2016. a, b

Ravindranath, N. H., Somashekhar, B. S., and Gadgil, M.: Carbon flow in Indian forests, Climatic Change, 35, 297–320,, 1997. a

Ravindranath, N. H., Murali, K. S., and Sudha, P.: Community forestry initiatives in Southeast Asia: a review of ecological impacts, International Journal of Environment and Sustainable Development (IJESD), 5, 1–11,, 2006. a, b

Richardson, A. D., Keenan, T. F., Migliavacca, M., Ryu, Y., Sonnentag, O., and Toomey, M.: Climate change, phenology, and phenological control of vegetation feedbacks to the climate system, Agr. Forest Meteorol., 169, 156–173,, 2013. a

Ricklefs, R. E.: The economy of nature, Macmillan, New York, W.H. Freeman, 2008. a, b, c

Rodgers, W. A. and Panwar, H. S.: Planning a wildlife protected area network in India, A report. Wildlife Institute of India, Dehradun, 1988. a

Running, S. W. and Hunt Jr, E. R.: Generalization of a forest ecosystem process model for other biomes, BIOME-BCG, and an application for global-scale models, in: Scaling Physiological Processes: Leaf to Globe, edited by: Ehleringer, J. R. and Field, C. B., Academic Press Inc., San Diego, CA, USA, 141–158, 1993. a

Rustad, L., Campbell, J., Marion, G., Norby, R., Mitchell, M., Hartley, A., Cornelissen, J., and Gurevitch, J.: A meta-analysis of the response of soil respiration, net nitrogen mineralization, and aboveground plant growth to experimental ecosystem warming, Oecologia, 126, 543–562,, 2001. a

Saatchi, S. S., Harris, N. L., Brown, S., Lefsky, M., Mitchard, E. T., Salas, W., Zutta, B. R., Buermann, W., Lewis, S. L., and Hagen, S.: Benchmark map of forest carbon stocks in tropical regions across three continents, P. Natl. Acad. Sci. USA, 108, 9899–9904,, 2011. a, b

Saikia, P., Deka, J., Bharali, S., Kumar, A., Tripathi, O. P., Singha, L. B., Dayanandan, S., and Khan, M. L.: Plant diversity patterns and conservation status of eastern Himalayan forests in Arunachal Pradesh, Northeast India, Forest Ecosystems, 4, 28,, 2017. a

Sakschewski, B., Bloh, W., Boit, A., Rammig, A., Kattge, J., Poorter, L., Peñuelas, J., and Thonicke, K.: Leaf and stem economics spectra drive diversity of functional plant traits in a dynamic global vegetation model, Global Change Biol., 21, 2711–2725,, 2015. a, b

Sato, H., Itoh, A., and Kohyama, T.: SEIB–DGVM: A new Dynamic Global Vegetation Model using a spatially explicit individual-based approach, Ecol. Model., 200, 279–307,, 2007. a

Scheiter, S., Langan, L., and Higgins, S. I.: Next-generation dynamic global vegetation models: learning from community ecology, New Phytol., 198, 957–969,, 2013. a, b, c

Scheiter, S., Kumar, D., Corlett, R. T., Gaillard, C., Langan, L., Lapuz, R. S., Martens, C., Pfeiffer, M., and Kyle, T. W.: Climate change promotes transitions to tall evergreen vegetation in tropical Asia, Global Change Biol., 26, 5106–5124,, 2020. a, b

Schimel, D., Stephens, B. B., and Fisher, J. B.: Effect of increasing CO2 on the terrestrial carbon cycle, P. Natl. Acad. Sci. USA, 112, 436–441,, 2015. a

Sharkey, T. D., Bernacchi, C. J., Farquhar, G. D., and Singsaas, E. L.: Fitting photosynthetic carbon dioxide response curves for C3 leaves, Plant Cell Environ., 30, 1035–1040,, 2007. a

Simard, M., Pinto, N., Fisher, J. B., and Baccini, A.: Mapping forest canopy height globally with spaceborne lidar, J. Geophys. Res.-Biogeo., 116, G04021,, 2011. a, b

Sinha, S., Badola, H. K., Chhetri, B., Gaira, K. S., Lepcha, J., and Dhyani, P. P.: Effect of altitude and climate in shaping the forest compositions of Singalila National Park in Khangchendzonga Landscape, Eastern Himalaya, India, J. Asia-Pac. Biodivers., 11, 267–275,, 2018. a

Sitch, S., Friedlingstein, P., Gruber, N., Jones, S. D., Murray-Tortarolo, G., Ahlström, A., Doney, S. C., Graven, H., Heinze, C., Huntingford, C., Levis, S., Levy, P. E., Lomas, M., Poulter, B., Viovy, N., Zaehle, S., Zeng, N., Arneth, A., Bonan, G., Bopp, L., Canadell, J. G., Chevallier, F., Ciais, P., Ellis, R., Gloor, M., Peylin, P., Piao, S. L., Le Quéré, C., Smith, B., Zhu, Z., and Myneni, R.: Recent trends and drivers of regional sources and sinks of carbon dioxide, Biogeosciences, 12, 653–679,, 2015. a

Smith, B., Wårlind, D., Arneth, A., Hickler, T., Leadley, P., Siltberg, J., and Zaehle, S.: Implications of incorporating N cycling and N limitations on primary production in an individual-based dynamic vegetation model, Biogeosciences, 11, 2027–2054,, 2014. a

Soh, W. K., Yiotis, C., Murray, M., Parnell, A., Wright, I. J., Spicer, R. A., Lawson, T., Caballero, R., and McElwain, J. C.: Rising CO2 drives divergence in water use efficiency of evergreen and deciduous plants, Science Advances, 5, eaax7906,, 2019. a, b

Song, J., Wan, S., Piao, S., Knapp, A. K., Classen, A. T., Vicca, S., Ciais, P., Hovenden, M. J., Leuzinger, S., Beier, C., and Kardol, P.: A meta-analysis of 1119 manipulative experiments on terrestrial carbon-cycling responses to global change, Nat. Ecol. Evol., 3, 1309–1320,, 2019. a

Sperry, J. S., Venturas, M. D., Todd, H. N., Trugman, A. T., Anderegg, W. R., Wang, Y., and Tai, X.: The impact of rising CO2 and acclimation on the response of US forests to global warming, P. Natl. Acad. Sci. USA, 116, 25734–25744,, 2019. a

Stephenson, N.: Actual evapotranspiration and deficit: biologically meaningful correlates of vegetation distribution across spatial scales, J. Biogeogr., 25, 855–870,, 1998. a

Stevens, N., Lehmann, C. E., Murphy, B. P., and Durigan, G.: Savanna woody encroachment is widespread across three continents, Global Change Biol., 23, 235–244,, 2017. a, b, c

Terrer, C., Vicca, S., Stocker, B. D., Hungate, B. A., Phillips, R. P., Reich, P. B., Finzi, A. C., and Prentice, I. C.: Ecosystem responses to elevated CO2 governed by plant–soil interactions and the cost of nitrogen acquisition, New Phytol., 217, 507–522,, 2018. a

Terrer, C., Jackson, R. B., Prentice, I. C., Keenan, T. F., Kaiser, C., Vicca, S., Fisher, J. B., Reich, P. B., Stocker, B. D., Hungate, B. A., and Penuelas, J.: Nitrogen and phosphorus constrain the CO2 fertilization of global plant biomass, Nat. Clim. Change, 9, 684–689,, 2019. a, b

Thuiller, W., Albert, C., Araújo, M. B., Berry, P. M., Cabeza, M., Guisan, A., Hickler, T., Midgley, G. F., Paterson, J., Schurr, F. M., and Sykes, M. T.: Predicting global change impacts on plant species' distributions: future challenges, Perspect. Plant Ecol., 9, 137–152,, 2008. a

Tian, H., Melillo, J., Kicklighter, D., McGuire, A., and Helfrich, J.: The sensitivity of terrestrial carbon storage to historical climate variability and atmospheric CO2 in the United States, Tellus B, 51, 414–452,, 1999. a

Tuanmu, M.-N. and Jetz, W.: A global 1 km consensus land-cover product for biodiversity and ecosystem modelling, Global Ecol. Biogeogr., 23, 1031–1045,, 2014. a, b

Urban, J., Ingwers, M., McGuire, M. A., and Teskey, R. O.: Stomatal conductance increases with rising temperature, Plant Signaling & Behavior 12, e1356534,, 2017. a

Van Vuuren, D. P., Edmonds, J., Kainuma, M., Riahi, K., Thomson, A., Hibbard, K., Hurtt, G. C., Kram, T., Krey, V., and Lamarque, J.-F.: The representative concentration pathways: an overview, Climatic Change, 109, 5,, 2011. a, b

Verstraete, M. M., Scholes, R. J., and Smith, M. S.: Climate and desertification: looking at an old problem through new lenses, Front. Ecol. Environ., 7, 421–428,, 2009. a

Walker, M. D., Wahren, C. H., Hollister, R. D., Henry, G. H., Ahlquist, L. E., Alatalo, J. M., Bret-Harte, M. S., Calef, M. P., Callaghan, T. V., Carroll, A. B., and Epstein, H. E.: Plant community responses to experimental warming across the tundra biome, P. Natl. Acad. Sci. USA, 103, 1342–1346, 2006. a

Wang, X., Wang, T., Liu, D., Guo, H., Huang, H., and Zhao, Y.: Moisture-induced greening of the South Asia over the past three decades, Global Change Biol., 23, 4995–5005,, 2017. a

Warren, J. M., Pötzelsberger, E., Wullschleger, S. D., Thornton, P. E., Hasenauer, H., and Norby, R. J.: Ecohydrologic impact of reduced stomatal conductance in forests exposed to elevated CO2, Ecohydrology, 4, 196–210, 2011. a, b

Warszawski, L., Frieler, K., Huber, V., Piontek, F., Serdeczny, O., and Schewe, J.: The inter-sectoral impact model intercomparison project (ISI–MIP): project framework, P. Natl. Acad. Sci. USA, 111, 3228–3232,, 2014. a

Whittaker, R. H.: Communities and Ecosystems, 2nd Edn., MacMillan Publishing Co., New York, 1975. a, b, c

Wingfield, J. C.: Ecological processes and the ecology of stress: the impacts of abiotic environmental factors, Funct. Ecol., 27, 37–44,, 2013. a

Woodrow, I. E. and Berry, J.: Enzymatic regulation of photosynthetic CO2, fixation in C3 plants, Annu. Rev. Plant Phys., 39, 533–594, 1988. a

Wright, I. J., Reich, P. B., Cornelissen, J. H., Falster, D. S., Groom, P. K., Hikosaka, K., Lee, W., Lusk, C. H., Niinemets, Ü., and Oleksyn, J.: Modulation of leaf economic traits and trait relationships by climate, Global Ecol. Biogeogr., 14, 411–421,, 2005. a

Yuan, W., Zheng, Y., Piao, S., Ciais, P., Lombardozzi, D., Wang, Y., Ryu, Y., Chen, G., Dong, W., Hu, Z., and Jain, A. K.: Increased atmospheric vapor pressure deficit reduces global vegetation growth, Science Advances, 5, eaax1396,, 2019.  a

Zemunik, G., Turner, B. L., Lambers, H., and Laliberté, E.: Diversity of plant nutrient-acquisition strategies increases during long-term ecosystem development, Nat. Plants, 1, 15050,, 2015. a

Zhang, K., Kimball, J. S., Nemani, R. R., and Running, S. W.: A continuous satellite-derived global record of land surface evapotranspiration from 1983 to 2006, Water Resour. Res., 46, W09522,, 2010. a, b

Zhang, Y., Peña-Arancibia, J. L., McVicar, T. R., Chiew, F. H., Vaze, J., Liu, C., Lu, X., Zheng, H., Wang, Y., and Liu, Y. Y.: Multi-decadal trends in global terrestrial evapotranspiration and its components, Sci. Rep.-UK, 6, 19124,, 2016. a

Short summary
In this paper, we investigated the impact of climate change and rising CO2 on biomes using a vegetation model in South Asia, an often neglected region in global modeling studies. Understanding these impacts guides ecosystem management and biodiversity conservation. Our results indicate that savanna regions are at high risk of woody encroachment and transitioning into the forest, and the bioclimatic envelopes of biomes need adjustments to account for shifts caused by climate change and CO2.
Final-revised paper