1. Introduction
About 41% of Earth’s land surface is composed of drylands (areas suffering from permanent or seasonal water deficiency; Lian et al. 2021). Drylands are some of the most sensitive areas to global warming and nurture over 38% of Earth’s population (Huang et al. 2016; Yao et al. 2020). About half of dryland residents live under the poverty threshold according to the United Nations (Reynolds et al. 2007; Lian et al. 2021). Therefore, the well-being of people, associated ecosystem services, and societal goods are particularly vulnerable to anthropogenic warming (Huang et al. 2017). Thus, past and future changes in drylands under global warming are a growing cause for concern (Dai 2013; Park et al. 2018; Pokhrel et al. 2021; Vicente-Serrano et al. 2020; Li et al. 2021; Zhang et al. 2017; Li et al. 2020; Zhang et al. 2019; Yao et al. 2020). Exploring the driving effects of anthropogenic climate change on changes in global dryland extent can greatly enhance our understanding in dryland evolution under future warming. Detection and attribution analysis separating the anthropogenic climate change signal in global dryland changes will deepen our overall knowledge of human activities on dryland extent and of how that may change with anthropogenic climate change.
Spatial and temporal changes in dryland extent and intensity (aridification) have been extensively examined. Previous studies have assessed changes in the spatial patterns of drylands between the current and future periods (Feng and Fu 2013; Berg et al. 2017; Cheng and Huang 2016; Cheng et al. 2015; Schlaepfer et al. 2017; Greve et al. 2014; Burrell et al. 2020; Gu et al. 2019c; Dai 2013; Pokhrel et al. 2021; Huang et al. 2016; Park et al. 2018). Feng and Fu (2013) expected a drier future over most land areas defined by atmospheric aridity. Dai (2013) projected soil moisture deficit of the top 10-cm layer in North America, South America, Europe, southern Africa, and Australia. Enhanced soil moisture drying occurred in humid transitional regions and then expanded to dry transitional regions (Cheng and Huang 2016). Pokhrel et al. (2021) indicated that many global land areas, especially in the Southern Hemisphere, are likely to witness a future reduction in terrestrial water storage. However, global land areas are also projected to experience wetting when using an ecohydrological index (based on surface soil moisture, transpiration, and leaf area index) as a proxy for drylands (Berg and McColl 2021). Based on the analysis from more than 300 combinations of various precipitation, evapotranspiration, and potential evaporation datasets, Greve et al. (2014) found that the “dry gets drier, wet gets wetter” pattern only occurs in 10.8% of global land areas, and the opposite pattern can be found in 9.5% of global land areas. Despite these studies, it is still unknown whether spatial changes in dryland areas can be attributed to human-induced climate change or natural climate variability.
The results of past and future changes in dryland extent differ in previous studies using different indices to define drylands from atmospheric, hydrological, agricultural, and ecological viewpoints (Huang et al. 2016; Lian et al. 2021; Dai 2013; Pokhrel et al. 2021; Berg and McColl 2021). Based on atmospheric aridity indices, dryland expansion is expected to accelerate, with a 23% increase under a high-emission scenario (Huang et al. 2016; Feng and Fu 2013). The self-calibrated Palmer drought severity index shows an 8% increase of global drought areas from 1950 to 2010 (Dai 2013). Land areas with decreases in terrestrial water storage and dry-season water availability are also expanding (Padrón et al. 2020; Pokhrel et al. 2021). In contrast, Berg and McColl (2021) argued that no expansions in global drylands are projected under greenhouse warming when using an ecohydrological index, which is consistent with greening and enhanced vegetation activity (He et al. 2019; Fensholt et al. 2012; Zhu et al. 2016; Piao et al. 2020). Recently, Lian et al. (2021) reviewed the characteristics of dryland aridity changes using atmospheric, hydrological, agricultural, and ecological indices, and confirmed that dryland areas are projected to increase in the future except when using the ecology-based index. Nevertheless, current research has not yet provided quantitative evidence to support human-induced dryland expansion or quantified the role of different external forcings.
Since both the spatial patterns and extent of dryland areas are projected to change under future warming, an equally important factor is the potentially affected dryland population. Many studies have estimated the exposure of population to droughts at global and regional scales (Chen et al. 2018; Liu et al. 2018; AghaKouchak et al. 2021; Kasprzyk et al. 2009), but few have focused on changes in dryland population (Koutroulis 2019; Stavi et al. 2021; Huang et al. 2016). For example, Huang et al. (2016) projected dryland population growth by 50% in the future under a high-emission scenario. Koutroulis (2019) estimated that 1.9 billion people could avoid living in drylands if the warming increase is reduced from 4° to 1.5°C. However, these studies did not separate the roles of the spatial evolution of drylands from alterations in dryland population (see Fig. 1). This separation can help us take a closer look at the potential effects of dryland changes on population and provide decision-makers with more specific options for adaptation and mitigation.
Schematic diagram of population changes arising from population density and spatial extent of drylands between a warm future (i.e., 20-yr) period and the reference period 1961–90, respectively. (left) The spatial extent of drylands identified during the reference period and a future warming period of +3.5°C under SSP585. The letters are as follows: A indicates the shared dryland extent between the two periods, A0 (Ai) is the extent in the 1961–90 reference period (in the future ith warm period), and P0 (Pi) is the population density in the reference period (in the ith warm period).
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
In this study, following the protocol that is widely used to define drylands from atmospheric aridity, we employ the aridity index [AI; the ratio of the annual total precipitation (PRCP) to potential evapotranspiration (PET)] to identify drylands as the areas with AI ≤ 0.65. Then, we address the following questions:
-
Can spatial patterns of changes in AI be attributed to anthropogenic climate change?
-
Is there quantitative evidence to support human-induced dryland expansion, and to what extent does anthropogenic climate change contribute to this expansion?
-
What are the separate roles of dryland expansion and population alteration in the changes of residents in drylands?
2. Data
a. Observation-based meteorological data
Observation-based datasets used to calculate the aridity index.
Due to the differences in spatial resolution and record length among the four datasets, we used bilinear interpolation (Gu et al. 2019a) to regrid these datasets as a common resolution of 0.5°, and extracted data from 1950–2014 as the common study period. The ensemble mean (OBS-EM) of these four datasets is calculated for the following detection and attribution analysis.
b. CMIP6 data
The recently released monthly model simulations from phase 6 of the Coupled Model Intercomparison Project (CMIP6) are obtained for the following analyses. The variables provided by CMIP6 are surface upward latent heat flux, surface upward sensible heat flux, wind speed, relative humidity, and minimum/maximum near-surface air temperature under difference forcings: natural (NAT) forcing, greenhouse gases (GHG) forcing, and historical (ALL) forcing (including NAT and anthropogenic forcings). Most of the models under these forcings have simulations ending in 2014. The simulations from the PiControl (CTL) experiments (where the GHG level is set to preindustrial conditions) are also obtained. Results based on the CTL experiments indicate the impact of natural climate variability because they do not include any forcing. Information on the models used in this study can be found in Table 2.
Model simulations used to calculate the aridity index. An asterisk (*) indicates that the corresponding model has the data for a given experiment; a long dash (—) indicates that it does not.
We also obtained projections for 2015–2100 from CMIP6 under four SSP-RCP (shared socioeconomic pathway–representative concentration pathway) scenarios, namely, SSP126, SSP245, SSP370, and SSP585. For example, SSP126 is the combination of SSP1 and RCP26. RCPs represent the radiative forcing paths, with the radiative forcing reaching 2.6 W m−2 (RCP26), 4.5 W m−2 (RCP45), 7.0 W m−2 (RCP70), or 8.5 W m−2 (RCP85) at the end of the twenty-first century, (Taylor et al. 2012). SSPs represent various global socioeconomic development scenarios, including sustainability (taking the green road; SSP1), middle of the road (SSP2), regional rivalry (a rocky road; SSP3), and fossil-fueled development (taking the highway; SSP5) (Riahi et al. 2017). Populations are projected to increase around 2050, 2060, and 2070 in the SSP1, SSP5, and SSP2, respectively, and then decline by 2100, while this increase in population will persist through 2100 in the SSP3. To ensure the same weight among these models, only the models with realization “r1i1p1fn” are used; for more information, refer to Table 3. The PM equation (Allen et al. 1998; Yang et al. 2019; Scheff and Frierson 2014; Fu and Feng 2014) is used to estimate the PET based on model simulations. Note that Rn* is calculated by sensible and latent heat in modeled data (Fu and Feng 2014; Yang et al. 2019; Scheff and Frierson 2014).
Model projections used to calculate the aridity index. An asterisk (*) indicates that the corresponding model has the data for a given experiment; a long dash (—) indicates that it does not.
c. Population data
The yearly population data (0.5° resolution) are collected from the InterSectoral Impact Model Intercomparison Project (ISIMIP). The ISIMIP provides historical (1861–2005) and future (2006–2100) population data under the SSP1, SSP2, SSP3, and SSP5, respectively. The population data are available at https://data.isimip.org/.
3. Methods
a. Correlation-based attribution approach
The attribution of changes in the spatial distribution of ΔAI to anthropogenic climate change is assessed by following a correlation-based attribution approach (Gudmundsson et al. 2017b; Qian and Zhang 2015; Padrón et al. 2020). We estimate the spatial Spearman correlation of ΔAI between observations and model simulations to quantify responses of spatial change to external forcings. The correlations between observations and CTL simulations depict the distribution of impacts only from natural climate variability. The null hypothesis is that there is no signal of anthropogenic forcing in the changes in spatial patterns of AI; that is, these changes are attributed to natural climate change only. The null hypothesis is rejected (an anthropogenic signal is detectable) if the spatial correlations between observations and simulations under ALL/GHG forcing exceed the 95th percentile of correlations between observations and CTL simulations. Furthermore, the presence of an anthropogenic signal is confirmed if the spatial correlations between observations and simulations under NAT forcing are below the 95th percentile.
b. The optimal fingerprint method
The spatial extent of global drylands (i.e., areas with AI ≤ 0.65) is estimated annually based on observations and simulations under different forcings. The Spearman correlations of dryland area extent (time series) between observations and CTL, ALL, GHG, and NAT simulations are estimated, respectively. Similarly, we first employ the correlation-based attribution approach to detect whether an anthropogenic signal is detectable in temporal changes of dryland area extent (Gudmundsson et al. 2017a; Padrón et al. 2020).
c. Population assessment in drylands
4. Results
a. Detection of anthropogenic climate change in spatial patterns of change
Figure 2 shows spatial patterns (excluding Antarctica and Greenland) of changes in AI, PRCP, and PET between the mean values during 1990–2014 and 1950–74 from multidata mean observations and multimodel mean simulations under historical forcings (i.e., ALL, GHG, and NAT). Consistent spatial patterns can be found when the periods are 1985–2009 and 1950–74 (Fig. S1 in the online supplemental material), 1985–2009 and 1955–79 (Fig. S2), and 1990–2014 and 1955–79 (Fig. S3). Observations indicate the global land surface has mostly witnessed decreases in AI (62.3% of land grids) during the recent period relative to during the past period (Fig. 2a). Areas with decreases in AI are mainly in the East Asia monsoon region, Africa, and eastern Australia, where a decreasing trend in AI during 1948–2005 was also found by Huang et al. (2016). The observed PRCP shows decreases in only 48% of land grid cells, which are also mainly located in the East Asia monsoon region, Africa, and eastern Australia (Fig. 2e), while increases in observed PET are widespread across the globe (84% of land grids; Fig. 2i). The prominent increases in PET are favorable to decreases in atmospheric aridity. However, the spatial patterns of changes in observed AI are more closely controlled by PRCP, as seen by their highly consistent spatial patterns (i.e., 83.2% of land grids have the same change direction; Figs. 2a,e). For observed AI and PET, only 34.5% of the total land grids have the same sign (Figs. 2a,i). Feng and Fu (2013) found that increasing PET dominates the areas with translation from humid to dry subhumid and from dry subhumid to semiarid while increasing precipitation plays a dominant role in the regions with reducing aridity.
(a)–(d) Spatial patterns of changes in the mean of aridity index (AI), (e)–(h) annual total precipitation (PRCP; mm), and (i)–(l) annual total potential evapotranspiration (PET; mm) between the two periods (i.e., 1990–2014 and 1950–74, respectively) from multidata mean observations (OBS-EM) and historical simulations (i.e., ALL, GHG, and NAT forcings, respectively). The blue curve on the right side of panels indicates zonal average change. The stippling indicates areas where 75% of the models agree on the sign of the change (observations indicate the agreement of all observed datasets).
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
There is an agreement in the sign of the changes in PRCP (PET) between observations and simulations under ALL forcing; 55.7% (75%) of land areas have the same direction of change. This agreement in PRCP (PET) is even higher between observations and GHG simulations, and the corresponding percentage is 56.7% (83.6%). However, simulated PRCP (PET) that excludes human-induced impacts under NAT forcing shows much weaker consistency with observations. The above results imply that anthropogenic climate change is likely to be the underlying reason for the observed pattern of change in PRCP and PET. Previous studies have shown that the spatial distribution of long-term changes in precipitation and evapotranspiration at the global scale is well captured by model simulations (Min et al. 2011; Madakumbura et al. 2021; J. Liu et al. 2021).
Nevertheless, the sign of spatial patterns in changes in observed AI is poorly captured by historical ALL simulations (Figs. 2a,b). Only 47.2% of land grids show consistent directions of change between observed and ALL simulated AI; opposite directions of change are found in North America, South America, northern Africa, and the East Asia monsoon region. The poor consistency of spatial change patterns between observations and ALL simulations was also found in previous studies (Huang et al. 2016; Feng and Fu 2013). The spatial patterns of change are highly similar between ALL- and GHG-simulated AI (Figs. 2b,c), implying that GHG forcing plays an important role in the changes under ALL forcing (which includes GHGs and other external forcings). We notice that the increase of observed AI is captured in simulations under NAT forcing, such as North America and South America (Figs. 2a,d), but not ALL forcing. The possible reasons for the poor ability of models to simulate the spatial patterns of changes in AI are model uncertainties and noises from natural climate variability (Monerie et al. 2020; Fyfe et al. 2021). As natural climate variability affects both observations and simulations, the agreement of spatial change patterns between them could be reduced by climate noise. On the other hand, the global land is dominated by increases in PRCP and PET under both ALL and GHG forcings, which could lead to mixed changes in P/PET (i.e., AI).
To formally detect an anthropogenic signal in the spatial change patterns, we assess whether the Spearman correlation of spatial patterns between observations and simulations under external forcings is greater than the correlations obtained from natural climate variability (Fig. 3). We find that the spatial correlation of PRCP (PET) between the ensemble mean of multiple observations (OBS-EM) and GHG simulations (red circle on the far right boxplot) exceeds the 90th (80th) percentile of correlations between observations and the 226 CTL simulation segments (Figs. 3b,c). This finding suggests that the GHG signal is detectable in the spatial patterns of changes in PRCP (PET) at the 90% (80%) confidence level (although both are lower than 95% confidence level). The spatial correlation of PRCP (PET) between OBS-EM and NAT simulations smaller than 80th percentile of correlations between OBS-EM and CTL simulations further confirms the impacts of human-induced climate change.
Correlation-based attribution analyses for changes in spatial patterns of (a) AI, (b) PRCP, and (c) PET (as illustrated in Fig. 2). The gray bars represent the distribution of spatial correlations between observations and all 65-yr nonoverlapping segments from CTL simulations. The orange, red, and blue dots indicate the spatial correlations between observations and the multimodel mean simulations from ALL, GHG, and NAT forcings, respectively.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
As shown in previous studies, anthropogenic climate change is detectable in the temporal increases in global precipitation and evapotranspiration (Wu et al. 2013; J. Liu et al. 2021). Based on the Clausius–Clapeyron relationship, anthropogenic warming is expected to increase atmospheric water holding capacity, and thus enhance precipitation totals (Donat et al. 2016; Wu et al. 2013). For global drylands, annual precipitation totals increase by ∼15% per 1°C of global warming (Donat et al. 2016). Besides the thermodynamic effects (i.e., changes in atmospheric moisture) of global warming, warming-driven dynamic effects also work on precipitation changes (Paik et al. 2020; Pfahl et al. 2017). For example, Pfahl et al. (2017) indicated that the thermodynamics cause a consistent increase in precipitation over most global areas, while dynamics amplify or weaken this increase depending on regions. Global hydrological cycle is also driven by anthropogenic-induced perturbations to the surface energy budget (Wu et al. 2010, 2013). GHGs can increase net surface radiation by intercepting more outgoing longwave radiation, and then enhance evapotranspiration (Dai et al. 2018; J. Liu et al. 2021). Rising atmospheric CO2 can reduce stomatal conductance in plants and lead to a lower transpiration rate, but simultaneously also increase transpiration by expanding plant leaf area (Cui et al. 2020; Ukkola et al. 2016; Yang et al. 2019). Our results indicate that it is likely (>80% probability) that GHG forcing is driving not only the observed increases but also the spatial changes of PRCP and PET.
Although a detectable GHG signal is found in spatial changes of PRCP and PET, this signal is not detected in the observed changes in AI (Fig. 3a). The spatial correlation between AI observations and GHG simulations is lower than the 80th percentile of correlations from CTL simulations. There is a very low spatial correlation between AI observations and simulations under ALL forcing, further confirming the poor ability of models to capture observed spatial change patterns (as mentioned above). This weak spatial correlation between observations and simulations under ALL forcing is also found in PRCP and PET. Possible reasons are that 1) model uncertainties may cause a disagreement between observations and simulations (Fyfe et al. 2021; Padrón et al. 2020), 2) the counteracting effect of anthropogenic aerosols and GHG may contribute to a more complex spatial change (Lau et al. 2017; Touma et al. 2021), and 3) the uncertainties among the four observation datasets may affect the response of model simulations to observed spatial patterns (Sun et al. 2018). There is an evident discrepancy in the response of different observational datasets to external forcings, for example, the GHG signal is detectable in the spatial changes of AI based on the UD and GLDAS (TerraClimate) dataset at 80% confidence level (Fig. 3a).
Although models have limited ability to simulate regional patterns of change, they do capture the temporal evolution of global land areas with drying atmospheric aridity (ΔAI < 0) during a running 20-yr period relative to the reference period of 1961–90 (Fig. 4). In recent decades, observations show a significant increase in the global proportion of land areas with drying atmospheric aridity. The simulations under ALL forcing barely capture this increasing land area; they underestimate the observed increasing rate, which the GHG simulations simulate well. When excluding human-induced forcings (NAT simulations), the increase in land area with drying atmospheric aridity is not detected. These results imply that an increasing proportion of land area is experiencing atmospheric drying, potentially caused by human-induced climate change. Land areas with drying atmospheric aridity are projected to continue expanding in the future under all four scenarios (although SSP126 may see a decrease in the second half of the twenty-first century).
Temporal evolution of the global proportion of land areas undergoing aridification [ΔAI < 0; see Eq. (1)] based on (a) observations and historical simulations and (b) projections under the four scenarios. Areas with ΔAI < 0 have drying atmospheric aridity during a running 20-yr window relative to the reference period 1961–90. The purple ribbon for observations indicates the extent of maximum and minimum values from the four observed datasets. The colored ribbons for the models indicate the 25th and 75th percentiles of the multimodel simulations.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
b. Detection and attribution of temporal changes in dryland areas
Following the increase in land areas with drying atmospheric aridity in recent decades (Fig. 4a), we then focus on changes in the extent of dryland areas (i.e., AI ≤ 0.65) in both observations and simulations (Fig. 5). As shown in previous studies (Huang et al. 2016; Dai 2013), the extent of drylands significantly increased at a rate of 5.65% century−1 in observations and 1.59% century−1 in simulations under ALL forcing during 1950–2014. This increase in dryland extent is more evident in GHG simulations (i.e., 4.47% century−1) than in ALL simulations. Both NAT and CTL simulations show a slight decrease in dryland area extent (Figs. 5d,e) and therefore do not explain the observed dryland expansion. Additionally, the observed trend and the GHG simulated trend are outside the range or at the high end (i.e., above the 99th percentile) of the trend expected from natural climate variability (see gray bar plots in Fig. 5f), which is quantified from CTL simulations. These results indicate that the GHG emission signal is strong enough to be detected in the observed increase in dryland extent.
Long-term changes in dryland areas in observations and model simulations. (a)–(e) The changes estimated from observations and model simulations under ALL, GHG, NAT, and CTL forcings, respectively. In (a), the gray shadow is the spread of maximum and minimum values from the four observed datasets. In (b)–(e), the gray shadow indicates the 25th and 75th percentiles of multimodel simulations. (f) The colored vertical lines show the trends in observations and model simulations under external forcings; blue text is the corresponding percentiles in the probability density distribution of trends in all 65-yr nonoverlapping segments from CTL simulations (see gray bars).
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
We first employ the correlation-based attribution approach (Gudmundsson et al. 2017a; Padrón et al. 2020) to conduct a quantitative assessment of anthropogenic impacts (Fig. 6). The Spearman correlations of drylands area time series between observations and ALL simulations are relatively high, that is, above the 85th percentile of correlations between observations and the CTL simulations (for the composite observations, Princeton, and UD and GLDAS datasets). This percentile is notably higher (i.e., at least above the 98th percentile) for the correlations between observations and GHG simulations (for the composite observations and all the individual datasets), suggesting that there is a high probability that GHG emissions are driving the observed increase in dryland extent. Furthermore, the NAT simulations are weakly correlated with the observed pattern (i.e., well within the range spanned by CTL simulations), implying that there is a low probability that the NAT forcing is driving the observed pattern.
Correlation-based attribution analyses for long-term changes in dryland extent based on (a) CRU, (b) TerraClimate, (c) Princeton, (d) UD, and (e) multimean observations. The bar plots are the Spearman correlations of dryland area time series between observations and CTL simulations during the period 1950–2014. The colored vertical lines show correlations between observations and model simulations under external forcings; blue texts are the corresponding percentiles in the probability density distribution of correlations between observations and CTL simulations.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
In addition to the correlation analysis, we employ the optimal fingerprinting method to further quantify the contributions of external forcings to the long-term changes in observed dryland extent (Fig. 7). We first regressed the observations on the individual external forcings (i.e., ALL, GHG, and NAT), respectively, and conducted one-signal detection and attribution (Figs. 7a,d). Both the OLS and TLS estimations show that the scaling factors of ALL and GHG forcings are significantly above zero at the 95% confidence level, while this is not the case for the NAT scaling factor. This is consistent with the results from the correlation-based detection, and strengthens the evidence that the signal of GHG emissions is detectable in dryland expansion. Both the scaling factors of ALL and GHG forcings are larger than 1 (i.e., around 2), indicating that model simulations underestimate the amplitude of observed dryland expansion.
Detection and attribution of long-term changes in dryland area by using the optimal fingerprint method. The scaling factors estimated from (a),(b) OLS and (d),(e) TLS methods, respectively. The scaling factors based on (a),(d) one-signal and (b),(e) two-signal analysis, and the error bars indicate the 5th and 95th percentile values of the corresponding scaling factor. (c),(f) The corresponding attribution results, except for the purple bar plot, which is the trend in multidata mean observations. The error bars in (c) and (f) indicate the 25%–75% interval of the bar plot values.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
The long-term changes in dryland extent show a very strong linear trend in both observations and simulations under ALL and GHG forcings (Figs. 5a–c). If this strong linear trend is distinguishable above the climate noise, this indicates that the one-signal analysis could be dominated by one forcing (Chen and Sun 2017; Gu et al. 2019b). Therefore, we conducted a two-signal analysis by regressing the observations on both ALL and NAT (GHG and NAT) forcings, simultaneously (Kong et al. 2020). As the observed patterns could be affected by several forcings, the two-signal analysis can also tell us whether the responses of ALL and GHG forcings can be separated from the NAT forcing (Figs. 7b,d). We find that the GHG scaling factor is significantly above zero in both OLS and TLS estimations, but the NAT factor is not. This result indicates that the response of GHG forcing is distinct from the NAT forcing and confirms that the result of one-signal analysis is robust.
We further estimate the contributions of individual external forcings to the observed change in dryland extent (Figs. 7c,f). Over the period 1950–2014, observed dryland areas have expanded by 3.67%. The ALL forcing (including both anthropogenic forcings and NAT forcing) has driven dryland expansion by 2.36% (1.81%) since 1950, according to the TLS (OLS) estimation. In contrast, the dryland expansion attributed to GHG forcing is 4.63% (4.52%) since 1950, which is comparable to the amplitude of change in observations.
c. Changes in dryland population
Since GHG emissions drive both global warming and dryland expansion, we assess the evolution of dryland extent and population with future warming levels (Fig. 8). All four future emission scenarios project increasing dryland extent with global warming (Figs. 8a,b). The dryland areas are projected to increase by 1%, 2%, and 4.8% relative to the reference period of 1961–90, with warming levels of 1.5°, 2°, and 3.5°C, relative to preindustrial conditions, respectively. This projection is consistent with previous studies (Huang et al. 2016; Koutroulis 2019). Huang et al. (2016) projected 23% and 11% increases in dryland areas by the end of the twenty-first century under RCP8.5 and RCP4.5 relative to the reference period of 1961–90. Koutroulis (2019) used the outputs of a high-resolution model (i.e., HadGEM3A with a spatial resolution of ∼60 km) and projected the areal coverage of drylands would increase by 7% by the end of twenty-first century under RCP8.5.
Projected changes in (a),(b) dryland area and (c),(d) population during future warming periods relative to the reference period of 1961–90 under the four future emission scenarios. The light-colored curves are the projections from individual models and the corresponding dark-colored curves are the projections from multimodel mean.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
Both studies (Huang et al. 2016; Koutroulis 2019) also projected increases in the population of drylands during this continued expansion. Koutroulis (2019) indicated that an additional 1.9 billion people would live in dryland areas with an increase in warming from 1.5° to 4°C. However, our projections show that the increases in dryland population persist among all warming levels only under the SSP370 scenario (Figs. 8c,d). The population is projected to increase to 12.6 billion by the end of the twenty-first century in SSP370, which is much higher than other SSPs (i.e., 7 billion for SSP1, 9.4 billion for SSP2, and 7 billion for SSP5; Samir and Lutz 2017; Riahi et al. 2017). Under SSP370, the persistently growing population and dryland expansion lead to continuous increases in relative change in dryland population. For SSP245 and SSP585, we observe a shift in dryland population from increases to decreases when the warming level reaches above ∼3.5°C. These changes in dryland population are also found in China where the population exposure to drought is projected to increase up to the 2040s and then decrease until the end of twenty-first century under the four future scenarios, with the exception of SSP370 (Chen et al. 2021).
We hence investigate whether population density and/or the spatial extent of drylands dominate this future shift in dryland population (Fig. 9). The growth in dryland extent is projected to contribute to the increases in dryland population under all four warming scenarios (Figs. 9b,d). Projected increases caused by spatial extent under SSP370, for instance, are 5.1%, 6.4%, and 14.1% at warming levels of 1.5°, 2°, and 3.5°C, respectively. In comparison with dryland expansion, changes in population density play an even greater role in dryland population growth (Figs. 9a,c). For example, the population increases induced by population density under SSP370 are 69%, 106%, and 160% at warming levels of 1.5°, 2°, and 3.5°C, respectively. However, this persistent increase in population induced by changes in population density is found only under SSP370 (Figs. 9a,c). For SSP245 and SSP585, the trajectory is projected to start decreasing beyond a warming level of ∼3.5°C, due to the reduction in dryland population density. These results suggest that warming above ∼3.5°C under SSP245 and SSP585 could make the dryland areas uninhabitable by humans, and thus lead to a reduction in dryland population.
Projected change in dryland population arising from changes in (a),(c) population density and (b),(d) spatial extent of drylands during a future warming period relative to the reference period 1961–90 under the four future emission scenarios. The light-colored curves are the projections from individual models and the corresponding dark-colored curves are the projections from the multimodel mean.
Citation: Journal of Climate 35, 20; 10.1175/JCLI-D-22-0103.1
5. Discussion and conclusions
In this study, we conduct a formal detection and attribution analysis of the observed spatiotemporal changes in global dryland extent by comparing multiple observed and simulated datasets using a correlation-related approach and optimal fingerprint method. Correlation-related attribution analysis shows it is unlikely that the spatial changes of atmospheric aridity (i.e., AI) between the recent period (1990–2014) and the past period (1950–74) are caused by GHG emissions. However, the growing fraction of land areas with increasing atmospheric aridity (i.e., ΔAI < 0) is captured by model simulations only when the models are forced with GHG emissions, implying that anthropogenic climate change may be driving dryland expansion. Formal detection and attribution analysis further support that a GHG signal is detectable (at 95% confidence level) in the observed expansion of global drylands over the past decades. Specifically, the observed extent of global drylands increased by 3.67% from 1950 to 2014, and a growth of ∼4.5% can be attributed to GHG emissions. Detection and attribution analysis quantitatively supports that anthropogenic climate change is the dominant contributor to global dryland expansion, which will increase population in drylands. In comparison with the increase driven by dryland expansion, changes in dryland population density play a more important role. A reduction in dryland population density is expected to drive a subsequent decrease in dryland population if future warming exceeds ∼3.5°C under the SSP245 and SSP585 scenarios.
As mentioned in the introduction, the spatial patterns of change and dryland expansion differ when using other indices to define drylands. Besides AI, Lian et al. (2021) also used vapor pressure deficit (VPD), soil moisture, runoff, and gross primary production (GPP) to define drylands and found that VPD- and soil moisture-based drylands show obvious expansion in the future while runoff-based (GPP-based) dryland extent remains stable (significantly decreases). The consistent obvious dryland expansion between AI and VPD and soil moisture implies that the anthropogenic climate change signal may be also detectable in dryland expansion identified based on VPD and soil moisture. Our previous study has detected the influence of human activities on the past decrease in global soil moisture (Gu et al. 2019b). However, no expansion is projected in runoff- and GPP-based global drylands, indicating that greenhouse warming signals may not be detectable in observed dryland expansion. This is because the effects of increasing atmospheric CO2 on runoff and vegetation are still ambiguous and debatable (Zhang et al. 2021; Lian et al. 2021; Berg and McColl 2021; Yang et al. 2019). Here, we focus on atmospheric aridity and use the aridity index to define drylands. Other indices based on agricultural, hydrological, and ecological attributes are beyond the scope of our study.
The robustness of our results is mainly affected by uncertainties in observations and model simulations. Spatial changes in the aridity index, precipitation, and potential evapotranspiration exhibit clear differences among the four observed datasets (Fig. S5). The Spearman correlation of spatial change between each pair of the four datasets is within 0.46–0.73 for the aridity index, 0.48–0.87 for precipitation, and 0.11–0.62 for potential evapotranspiration (Fig. S5). Observational uncertainties do affect the detection of a GHG emissions signal. For example, for the aridity index, a GHG signal is not detectable in the multimean observed changes, but is detected (>90% probability) in the UD and GLDAS dataset (Fig. 3a). In contrast, for precipitation, the GHG signal is detectable (>90% probability) in the multimean observed change, but not in the UD and GLDAS dataset (Fig. 3b). The inconsistencies among these observed precipitation datasets are described in Sun et al. (2018).
We also note the inconsistency of spatial changes in the aridity index between observations and model simulations (see Figs. 2a,b). The poor ability of models to capture spatial changes in the aridity index is also found in Huang et al. (2016). Alongside the changes in AI-based drylands, soil organic carbon storage and emission are also changed; however, these carbon cycle processes are not included in some climate models, which hinders models from accurately reproducing these spatial patterns of change (Huang et al. 2016). The models also underestimate the amplitude of the observed dryland expansion (see the scaling factor of ALL forcing, which exceeds 1 in Figs. 7a and 7d). We use a correlation-based approach and anomalies to reduce the potential impact of biases in the absolute magnitude of simulations on our detection and attribution results (Gudmundsson et al. 2017b; Qian and Zhang 2015), and relative changes in comparison with the reference period to reduce biases in the dryland population estimates.
Human-induced dryland expansion is accompanied by serious land degradation (Reynolds et al. 2007; Burrell et al. 2020). The change of above 5 million km2 of drylands toward desertification was driven by anthropogenic climate change during 1982–2015 (Burrell et al. 2020). Dryland expansion and degradation mean less water availability, and population growth increase water demand. The interaction between less water availability and growing water use would put great pressure on ecosystems and human settlements in drylands. Therefore, adaptation measures are needed to mitigate the adverse effects of water shortages, such as demand management, water markets, water-saving techniques, and limiting carbon emissions (AghaKouchak et al. 2021). Policymakers and societies should prioritize water demands and uses to minimize the impacts on ecosystems and human settlements. Researchers could focus on the integration of water supplies, needs, and management under climate change, which is important to sustain the development of humans and ecosystems in drylands (AghaKouchak et al. 2021).
Acknowledgments.
This study is supported by the National Key R&D Program of China (Grant 2018YFA0605603), the National Natural Science Foundation of China (Grant 41901041), the Fundamental Research Funds for the Central Universities (lzujbky-2021-kb12), and the fund for State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment (Grant SKLLQG2108). S. Feng is supported by the Key Laboratory of Environmental Change and Natural Disaster, Ministry of Education, Beijing Normal University (Grant 12800-312230014). X. Gu is supported by the China Scholarship Council. L.J. Slater is supported by U.K. Research and Innovation (MR/V022008/1). The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Data availability statement.
The Climatic Research Unit datasets are available at https://data.ceda.ac.uk//badc/cru/data/cru_ts/, the TerraClimate datasets at https://data.nkn.uidaho.edu/dataset/monthly-climate-and-climatic-water-balance-global-terrestrial-surfaces-1958-2015, the Princeton datasets at http://hydrology.princeton.edu, the University of Delaware datasets at http://climate.geog.udel.edu/∼climate/html_pages/download.html#P2017, and the Global Land Data Assimilation System datasets at https://doi.org/10.5067/9SQ1B3ZXP2C5. The simulated data from phase 6 of the Coupled Model Intercomparison Project are available at https://esgf-node.llnl.gov/search/cmip6/. Population datasets are available at https://data.isimip.org/.
REFERENCES
Abatzoglou, J. T., S. Z. Dobrowski, S. A. Parks, and K. C. Hegewisch, 2018: TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data, 5, 170191, https://doi.org/10.1038/sdata.2017.191.
AghaKouchak, A., and Coauthors, 2021: Anthropogenic drought: Definition, challenges, and opportunities. Rev. Geophys., 59, e2019RG000683, https://doi.org/10.1029/2019RG000683.
Allen, M. R., and S. F. B. Tett, 1999: Checking for model consistency in optimal fingerprinting. Climate Dyn., 15, 419–434, https://doi.org/10.1007/s003820050291.
Allen, M. R., and P. A. Stott, 2003: Estimating signal amplitudes in optimal fingerprinting, part I: Theory. Climate Dyn., 21, 477–491, https://doi.org/10.1007/s00382-003-0313-9.
Allen, R. G., L. S. Pereira, D. Raes, and M. Smith, 1998: Crop evapotranspiration—Guidelines for computing crop water requirements. FAO Irrigation and Drainage Paper 56, Food and Agriculture Organization of the United Nations, 300 pp., http://www.fao.org/docrep/X0490E/X0490E00.htm.
Beaudoing, H., and M. Rodell, 2019: GLDAS Noah Land Surface Model L4 monthly 0.25 × 0.25 degree V2.0, Goddard Earth Sciences Data and Information Services Center (GES DISC), accessed 2020, https://doi.org/10.5067/9SQ1B3ZXP2C5.
Berg, A., and K. A. McColl, 2021: No projected global drylands expansion under greenhouse warming. Nat. Climate Change, 11, 331–337, https://doi.org/10.1038/s41558-021-01007-8.
Berg, A., J. Sheffield, and P. C. D. Milly, 2017: Divergent surface and total soil moisture projections under global warming. Geophys. Res. Lett., 44, 236–244, https://doi.org/10.1002/2016GL071921.
Burrell, A. L., J. P. Evans, and M. G. De Kauwe, 2020: Anthropogenic climate change has driven over 5 million km2 of drylands towards desertification. Nat. Commun., 11, 3853, https://doi.org/10.1038/s41467-020-17710-7.
Chen, H., and J. Sun, 2017: Contribution of human influence to increased daily precipitation extremes over China. Geophys. Res. Lett., 44, 2436–2444, https://doi.org/10.1002/2016GL072439.
Chen, J., Y. Liu, T. Pan, Y. Liu, F. Sun, and Q. Ge, 2018: Population exposure to droughts in China under the 1.5°C global warming target. Earth Syst. Dyn., 9, 1097–1106, https://doi.org/10.5194/esd-9-1097-2018.
Chen, L., and Coauthors, 2021: Future drought in CMIP6 projections and the socioeconomic impacts in China. Int. J. Climatol., 41, 4151–4170, https://doi.org/10.1002/joc.7064.
Cheng, S., and J. Huang, 2016: Enhanced soil moisture drying in transitional regions under a warming climate. J. Geophys. Res. Atmos., 121, 2542–2555, https://doi.org/10.1002/2015JD024559.
Cheng, S., X. Guan, J. Huang, F. Ji, and R. Guo, 2015: Long-term trend and variability of soil moisture over East Asia. J. Geophys. Res. Atmos., 120, 8658–8670, https://doi.org/10.1002/2015JD023206.
Cook, B. I., J. E. Smerdon, R. Seager, and S. Coats, 2014: Global warming and 21st century drying. Climate Dyn., 43, 2607–2627, https://doi.org/10.1007/s00382-014-2075-y.
Cui, J., S. Piao, C. Huntingford, X. Wang, X. Lian, A. Chevuturi, A. G. Turner, and G. J. Kooperman, 2020: Vegetation forcing modulates global land monsoon and water resources in a CO2-enriched climate. Nat. Commun., 11, 5184, https://doi.org/10.1038/s41467-020-18992-7.
Dai, A., 2013: Increasing drought under global warming in observations and models. Nat. Climate Change, 3, 52–58, https://doi.org/10.1038/nclimate1633.
Dai, A., T. Zhao, and J. Chen, 2018: Climate change and drought: A precipitation and evaporation perspective. Curr. Climate Change Rep., 4, 301–312, https://doi.org/10.1007/s40641-018-0101-6.
Donat, M. G., A. L. Lowry, L. V. Alexander, P. A. O’Gorman, and N. Maher, 2016: More extreme precipitation in the world’s dry and wet regions. Nat. Climate Change, 6, 508–513, https://doi.org/10.1038/nclimate2941.
Ekström, M., P. D. Jones, H. J. Fowler, G. Lenderink, T. A. Buishand, and D. Conway, 2007: Regional climate model data used within the SWURVE project—1: Projected changes in seasonal patterns and estimation of PET. Hydrol. Earth Syst. Serv., 11, 1069–1083, https://doi.org/10.5194/hess-11-1069-2007.
Feng, S., and Q. Fu, 2013: Expansion of global drylands under a warming climate. Atmos. Chem. Phys., 13, 10 081–10 094, https://doi.org/10.5194/acp-13-10081-2013.
Fensholt, R., and Coauthors, 2012: Greenness in semi-arid areas across the globe 1981–2007—An Earth Observing Satellite based analysis of trends and drivers. Remote Sens. Environ., 121, 144–158, https://doi.org/10.1016/j.rse.2012.01.017.
Fu, Q., and S. Feng, 2014: Responses of terrestrial aridity to global warming. J. Geophys. Res. Atmos., 119, 7863–7875, https://doi.org/10.1002/2014JD021608.
Fyfe, J. C., V. V. Kharin, B. D. Santer, J. N. S. Cole, and N. P. Gillett, 2021: Significant impact of forcing uncertainty in a large ensemble of climate model simulations. Proc. Natl. Acad. Sci. USA, 118, e2016549118, https://doi.org/10.1073/pnas.2016549118.
Greve, P., B. Orlowsky, B. Mueller, J. Sheffield, M. Reichstein, and S. I. Seneviratne, 2014: Global assessment of trends in wetting and drying over land. Nat. Geosci., 7, 716–721, https://doi.org/10.1038/ngeo2247.
Gu, X., J. Li, Y. D. Chen, D. Kong, and J. Liu, 2019a: Consistency and discrepancy of global surface soil moisture changes from multiple model-based data sets against satellite observations. J. Geophys. Res. Atmos., 124, 1474–1495, https://doi.org/10.1029/2018JD029304.
Gu, X., Q. Zhang, J. Li, V. P. Singh, J. Liu, P. Sun, and C. Cheng, 2019b: Attribution of global soil moisture drying to human activities: A quantitative viewpoint. Geophys. Res. Lett., 46, 2573–2582, https://doi.org/10.1029/2018GL080768.
Gu, X., Q. Zhang, J. Li, V. P. Singh, J. Liu, P. Sun, C. He, and J. Wu, 2019c: Intensification and expansion of soil moisture drying in warm season over Eurasia under global warming. J. Geophys. Res. Atmos., 124, 3765–3782, https://doi.org/10.1029/2018JD029776.
Gudmundsson, L., S. I. Seneviratne, and X. Zhang, 2017a: Anthropogenic climate change detected in European renewable freshwater resources. Nat. Climate Change, 7, 813–816, https://doi.org/10.1038/nclimate3416.
Gudmundsson, L., S. I. Seneviratne, and X. Zhang, 2017b: Anthropogenic climate change detected in European renewable freshwater resources. Nat. Climate Change, 7, 813–816, https://doi.org/10.1038/nclimate3416.
Harris, I., T. J. Osborn, P. Jones, and D. Lister, 2020: Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Sci. Data, 7, 109, https://doi.org/10.1038/s41597-020-0453-3.
He, B., S. Wang, L. Guo, and X. Wu, 2019: Aridity change and its correlation with greening over drylands. Agric. For. Meteor., 278, 107663, https://doi.org/10.1016/j.agrformet.2019.107663.
Hu, C., J. Xia, D. She, L. Li, Z. Song, and S. Hong, 2021: A new framework for the identification of flash drought: Multivariable and probabilistic statistic perspectives. Int. J. Climatol., 41, 5862–5878, https://doi.org/10.1002/joc.7157.
Huang, J., H. Yu, X. Guan, G. Wang, and R. Guo, 2016: Accelerated dryland expansion under climate change. Nat. Climate Change, 6, 166–171, https://doi.org/10.1038/nclimate2837.
Huang, J., and Coauthors, 2017: Dryland climate change: Recent progress and challenges. Rev. Geophys., 55, 719–778, https://doi.org/10.1002/2016RG000550.
Kasprzyk, J. R., P. M. Reed, B. R. Kirsch, and G. W. Characklis, 2009: Managing population and drought risks using many-objective water portfolio planning under uncertainty. Water Resour. Res., 45, W12401, https://doi.org/10.1029/2009WR008121.
Kong, D., X. Gu, J. Li, G. Ren, and J. Liu, 2020: Contributions of global warming and urbanization to the intensification of human-perceived heatwaves over China. J. Geophys. Res. Atmos., 125, e2019JD032175, https://doi.org/10.1029/2019JD032175.
Koutroulis, A. G., 2019: Dryland changes under different levels of global warming. Sci. Total Environ., 655, 482–511, https://doi.org/10.1016/j.scitotenv.2018.11.215.
Lau, W. K. M., K.-M. Kim, and L. Ruby Leung, 2017: Changing circulation structure and precipitation characteristics in Asian monsoon regions: Greenhouse warming vs. aerosol effects. Geosci. Lett., 4, 28, https://doi.org/10.1186/s40562-017-0094-3.
Li, C., B. Fu, S. Wang, L. C. Stringer, Y. Wang, Z. Li, Y. Liu, and W. Zhou, 2021: Drivers and impacts of changes in China’s drylands. Nat. Rev. Earth Environ., 2, 858–873, https://doi.org/10.1038/s43017-021-00226-z.
Li, R., N. Chen, X. Zhang, L. Zeng, X. Wang, S. Tang, D. Li, and D. Niyogi, 2020: Quantitative analysis of agricultural drought propagation process in the Yangtze River Basin by using cross wavelet analysis and spatial autocorrelation. Agric. For. Meteor., 280, 107809, https://doi.org/10.1016/j.agrformet.2019.107809.
Lian, X., and Coauthors, 2021: Multifaceted characteristics of dryland aridity changes in a warming world. Nat. Rev. Earth Environ., 2, 232–250, https://doi.org/10.1038/s43017-021-00144-0.
Liu, J., J. Zhang, D. Kong, X. Feng, S. Feng, and M. Xiao, 2021: Contributions of anthropogenic forcings to evapotranspiration changes over 1980–2020 using GLEAM and CMIP6 simulations. J. Geophys. Res. Atmos., 126, e2021JD035367, https://doi.org/10.1029/2021JD035367.
Liu, M., X. Xu, B. R. Scanlon, A. Y. Sun, and K. Wang, 2021: A modified evaporation model indicates that the effects of air warming on global drying trends have been overestimated. J. Geophys. Res. Atmos., 126, e2021JD035153, https://doi.org/10.1029/2021JD035153.
Liu, W., F. Sun, W. H. Lim, J. Zhang, H. Wang, H. Shiogama, and Y. Zhang, 2018: Global drought and severe drought-affected populations in 1.5°C and 2°C warmer worlds. Earth Syst. Dyn., 9, 267–283, https://doi.org/10.5194/esd-9-267-2018.
Madakumbura, G. D., C. W. Thackeray, J. Norris, N. Goldenson, and A. Hall, 2021: Anthropogenic influence on extreme precipitation over global land areas seen in multiple observational datasets. Nat. Commun., 12, 3944, https://doi.org/10.1038/s41467-021-24262-x.
Marvel, K., B. I. Cook, C. J. W. Bonfils, P. J. Durack, J. E. Smerdon, and A. P. Williams, 2019: Twentieth-century hydroclimate changes consistent with human influence. Nature, 569, 59–65, https://doi.org/10.1038/s41586-019-1149-8.
Min, S.-K., X. Zhang, F. W. Zwiers, and G. C. Hegerl, 2011: Human contribution to more-intense precipitation extremes. Nature, 470, 378–381, https://doi.org/10.1038/nature09763.
Monerie, P.-A., C. M. Wainwright, M. Sidibe, and A. A. Akinsanola, 2020: Model uncertainties in climate change impacts on Sahel precipitation in ensembles of CMIP5 and CMIP6 simulations. Climate Dyn., 55, 1385–1401, https://doi.org/10.1007/s00382-020-05332-0.
Padrón, R. S., and Coauthors, 2020: Observed changes in dry-season water availability attributed to human-induced climate change. Nat. Geosci., 13, 477–481, https://doi.org/10.1038/s41561-020-0594-1.
Paik, S., S. K. Min, X. Zhang, M. G. Donat, A. D. King, and Q. Sun, 2020: Determining the anthropogenic greenhouse gas contribution to the observed intensification of extreme precipitation. Geophys. Res. Lett., 47, e2019GL086875, https://doi.org/10.1029/2019GL086875.
Park, C.-E., and Coauthors, 2018: Keeping global warming within 1.5°C constrains emergence of aridification. Nat. Climate Change, 8, 70–74, https://doi.org/10.1038/s41558-017-0034-4.
Pfahl, S., P. A. O’Gorman, and E. M. Fischer, 2017: Understanding the regional pattern of projected future changes in extreme precipitation. Nat. Climate Change, 7, 423–427, https://doi.org/10.1038/nclimate3287.
Piao, S., and Coauthors, 2020: Characteristics, drivers and feedbacks of global greening. Nat. Rev. Earth Environ., 1, 14–27, https://doi.org/10.1038/s43017-019-0001-x.
Pokhrel, Y., and Coauthors, 2021: Global terrestrial water storage and drought severity under climate change. Nat. Climate Change, 11, 226–233, https://doi.org/10.1038/s41558-020-00972-w.
Polson, D., M. Bollasina, G. C. Hegerl, and L. J. Wilcox, 2014: Decreased monsoon precipitation in the Northern Hemisphere due to anthropogenic aerosols. Geophys. Res. Lett., 41, 6023–6029, https://doi.org/10.1002/2014GL060811.
Qian, C., and X. Zhang, 2015: Human influences on changes in the temperature seasonality in mid- to high-latitude land areas. J. Climate, 28, 5908–5921, https://doi.org/10.1175/JCLI-D-14-00821.1.
Reynolds, J. F., and Coauthors, 2007: Global desertification: Building a science for dryland development. Science, 316, 847–851, https://doi.org/10.1126/science.1131634.
Riahi, K., and Coauthors, 2017: The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview. Global Environ. Change, 42, 153–168, https://doi.org/10.1016/j.gloenvcha.2016.05.009.
Rodell, M., and Coauthors, 2004: The Global Land Data Assimilation System. Bull. Amer. Meteor. Soc., 85, 381–394, https://doi.org/10.1175/BAMS-85-3-381.
Samir, K. C., and W. Lutz, 2017: The human core of the shared socioeconomic pathways: Population scenarios by age, sex and level of education for all countries to 2100. Global Environ. Change, 42, 181–192, https://doi.org/10.1016/j.gloenvcha.2014.06.004.
Scheff, J., and D. M. W. Frierson, 2014: Scaling potential evapotranspiration with greenhouse warming. J. Climate, 27, 1539–1558, https://doi.org/10.1175/JCLI-D-13-00233.1.
Schlaepfer, D. R., and Coauthors, 2017: Climate change reduces extent of temperate drylands and intensifies drought in deep soils. Nat. Commun., 8, 14196, https://doi.org/10.1038/ncomms14196.
Sheffield, J., G. Goteti, and E. F. Wood, 2006: Development of a 50-year high-resolution global dataset of meteorological forcings for land surface modeling. J. Climate, 19, 3088–3111, https://doi.org/10.1175/JCLI3790.1.
Skeie, R. B., T. K. Berntsen, G. Myhre, K. Tanaka, M. M. Kvalevåg, and C. R. Hoyle, 2011: Anthropogenic radiative forcing time series from pre-industrial times until 2010. Atmos. Chem. Phys., 11, 11 827–11 857, https://doi.org/10.5194/acp-11-11827-2011.
Song, C., Y. Sheng, S. Zhan, J. Wang, L. Ke, and K. Liu, 2020: Impact of amplified evaporation due to lake expansion on the water budget across the inner Tibetan Plateau. Int. J. Climatol., 40, 2091–2105, https://doi.org/10.1002/joc.6320.
Stavi, I., and Coauthors, 2021: Food security among dryland pastoralists and agropastoralists: The climate, land-use change, and population dynamics nexus. Anthropocene Rev., https://doi.org/10.1177/20530196211007512, in press.
Sun, Q., C. Miao, Q. Duan, H. Ashouri, S. Sorooshian, and K.-L. Hsu, 2018: A review of global precipitation data sets: Data sources, estimation, and intercomparisons. Rev. Geophys., 56, 79–107, https://doi.org/10.1002/2017RG000574.
Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1.
Touma, D., S. Stevenson, F. Lehner, and S. Coats, 2021: Human-driven greenhouse gas and aerosol emissions cause distinct regional impacts on extreme fire weather. Nat. Commun., 12, 212, https://doi.org/10.1038/s41467-020-20570-w.
Ukkola, A. M., I. C. Prentice, T. F. Keenan, A. I. J. M. Van Dijk, N. R. Viney, R. B. Myneni, and J. Bi, 2016: Reduced streamflow in water-stressed climates consistent with CO2 effects on vegetation. Nat. Climate Change, 6, 75–78, https://doi.org/10.1038/nclimate2831.
Vicente-Serrano, S. M., S. M. Quiring, M. Peña-Gallardo, S. Yuan, and F. Domínguez-Castro, 2020: A review of environmental droughts: Increased risk under global warming? Earth-Sci. Rev., 201, 102953, https://doi.org/10.1016/j.earscirev.2019.102953.
Wilcox, L. J., E. J. Highwood, and N. J. Dunstone, 2013: The influence of anthropogenic aerosol on multi-decadal variations of historical global climate. Environ. Res. Lett., 8, 24033, https://doi.org/10.1088/1748-9326/8/2/024033.
Willmott, C. J., and K. Matsuura, 2018: Terrestrial air temperature and precipitation: 1900–2017 Gridded Monthly Time Series (V 5.01), accessed 2020, http://climate.geog.udel.edu/∼climate/html_pages/Global2017/README.GlobalTsT2017.html.
Wu, P., R. Wood, J. Ridley, and J. Lowe, 2010: Temporary acceleration of the hydrological cycle in response to a CO2 rampdown. Geophys. Res. Lett., 37, L12705, https://doi.org/10.1029/2010GL043730.
Wu, P., N. Christidis, and P. Stott, 2013: Anthropogenic impact on Earth’s hydrological cycle. Nat. Climate Change, 3, 807–810, https://doi.org/10.1038/nclimate1932.
Yang, Y., M. L. Roderick, S. Zhang, T. R. McVicar, and R. J. Donohue, 2019: Hydrologic implications of vegetation response to elevated CO2 in climate projections. Nat. Climate Change, 9, 44–48, https://doi.org/10.1038/s41558-018-0361-0.
Yao, J., H. Liu, J. Huang, Z. Gao, G. Wang, D. Li, H. Yu, and X. Chen, 2020: Accelerated dryland expansion regulates future variability in dryland gross primary production. Nat. Commun., 11, 1665, https://doi.org/10.1038/s41467-020-15515-2.
Zhang, C., Y. Yang, D. Yang, and X. Wu, 2021: Multidimensional assessment of global dryland changes under future warming in climate projections. J. Hydrol., 592, 125618, https://doi.org/10.1016/j.jhydrol.2020.125618.
Zhang, X., N. Chen, J. Li, Z. Chen, and D. Niyogi, 2017: Multi-sensor integrated framework and index for agricultural drought monitoring. Remote Sens. Environ., 188, 141–163, https://doi.org/10.1016/j.rse.2016.10.045.
Zhang, X., and Coauthors, 2019: Urban drought challenge to 2030 sustainable development goals. Sci. Total Environ., 693, 133536, https://doi.org/10.1016/j.scitotenv.2019.07.342.
Zhu, Z., and Coauthors, 2016: Greening of the Earth and its drivers. Nat. Climate Change, 6, 791–795, https://doi.org/10.1038/nclimate3004.