Abstract
The Afro-Asian summer monsoon (AfroASM) sustains billions of people living in many developing countries covering West Africa and Asia, vulnerable to climate change. Future increase in AfroASM precipitation has been projected by current state-of-the-art climate models, but large inter-model spread exists. Here we show that the projection spread is related to present-day interhemispheric thermal contrast (ITC). Based on 30 models from the Coupled Model Intercomparison Project Phase 6, we find models with a larger ITC trend during 1981–2014 tend to project a greater precipitation increase. Since most models overestimate present-day ITC trends, emergent constraint indicates precipitation increase in constrained projection is reduced to 70% of the raw projection, with the largest reduction in West Africa (49%). The land area experiencing significant increases of precipitation (runoff) is 57% (66%) of the raw projection. Smaller increases of precipitation will likely reduce flooding risk, while posing a challenge to future water resources management.
Similar content being viewed by others
Introduction
The Afro-Asian monsoon, a major component of the global monsoon system, consists of West African, South Asian and East Asian monsoon regions1,2,3,4,5. Projection of Afro-Asian summer monsoon (AfroASM) precipitation, which would greatly affect local freshwater resources for billions of people3,6,7, thus is crucial to climate change adaptation and mitigation activities. Unfortunately, the projections of climate models show a large spread8,9,10,11,12, hampering the assessment of regional climate change. To yield a more reliable future projection, understanding and reducing the uncertainty are of urgent need.
Great efforts have been devoted to the understanding of the source of uncertainties in monsoon precipitation changes. Multi-model ensemble (MME) projections of CMIP (Coupled Model Intercomparison Project) models suggest that AfroASM precipitation would increase by approximately 8% and 14% under medium and high emission scenarios in the long-term projection, with large inter-model spread of 1%~14% and 3%~25%, respectively4,11,13,14,15,16. The thermodynamic process related to increases in atmospheric moisture enhances precipitation robustly across models, while the dynamic process associated with changes of circulation contributes to large inter-model uncertainty7,11. Total uncertainty in future projection comes from distinct socioeconomic scenarios, stochastic internal variability and different model structures8,17,18,19,20,21. The model uncertainty explains more than 70% of total uncertainties of AfroASM precipitation changes in CMIP5 models in the long-term projection8,22. The uncertainty can even offset the reliability of monsoon precipitation projection, and hamper the use of this information for policy making8,23.
In recent years, the emergent constraint technique, which is based on the physical link between a modeled but observable variable in the present day and a projected variable in the future climate system, has been developed to reduce the projection uncertainty and improve the reliability of future projection24,25,26,27,28. So far, many factors, such as the SST over cold-tongue regions and the convection over western Pacific, have been used to constrain the CMIP5 projection of summer precipitation over the East Asian, South Asian and West African monsoon regions9,10,29,30,31.
In the ongoing CMIP6, while global land monsoon precipitation is projected to increase in the long term under high emission scenarios, there exists large uncertainty at regional scales, in particular over the AfroASM regions11. Recent studies reported the connection between the uncertainty of AfroASM precipitation changes and the increase of interhemispheric thermal contrast in the projection4,12,13,32. But how to constrain the projection and reduce the spread remains unknown. Given the fact that the monsoon precipitation in West Africa and Asia shows in-phase changes due to the modulation of interhemispheric thermal contrast (ITC) and sea surface temperature (SST) variation of North Atlantic on millennial1,33, centurial4,12,32, and decadal timescales3,34, we hypothesize that the spread of AfroASM precipitation can potentially be constrained by the large scale interhemispheric or land-sea thermal contrast. We examine this hypothesis by using the output of the new Scenario Model Intercomparison Project in CMIP635.
In this work, by constraining the spread and biases of ITC in a hierarchical statistical framework, precipitation increase in the constrained projection (0.57 ± 0.38 mm day−1, constrained projection with ±1σ across models) is about 70% of the raw projection, with the largest reduction in the West African monsoon region. About 10% of the inter-model uncertainty in future precipitation changes is reduced. Given that the emergent constraint improves the reliability in AfroASM precipitation projections, we further investigate the impacts of the constrained projection on the potential water availability. The fractions of land area that will experience a significant increase of precipitation and potential water availability are about 57% and 66% of the raw projection, respectively.
Results
Dominant uncertainty of AfroASM precipitation projection
CMIP6 models project a general increase in AfroASM precipitation under the high-emission scenario, Shared Socioeconomic Pathways (SSP) 5–8.5 (2050~2099), except for part of the West African monsoon region (Fig. 1; see Methods and Supplementary Table S1). The SSP5-8.5 scenario is a fossil-fuel development pathway, in which the anthropogenic radiative forcing will increase by 8.5 W m−2 at the end of 21st century36,37,38. The regional average of the AfroASM precipitation increase is 14% relative to the baseline (1965~2014), with a large inter-model spread (1%~27% for the 5th–95th ensemble range; Fig. 1a). The signal-to-noise ratio (SNR), defined as the ratio between the ensemble mean and inter-model standard deviation (STD) of the projected changes, is less than 1.5 over 90% of the Afro-Asian monsoon regions, demonstrating large inter-model uncertainties in the projected AfroASM precipitation (Fig. 1c).
To reveal the sources of model spread, we conduct an inter-model empirical orthogonal function (EOF) analysis on the projected changes of AfroASM precipitation (Fig. 2; see Methods). The leading principal component (PC1) accounts for 26% of inter-model variance (Fig. 2a and Supplementary Fig. S1). The leading uncertainty mode exhibits a spatially consistent increase of precipitation over the AfroASM domain, with a systematic enhancement of monsoon circulation from West Africa, through Indian Peninsular to East Asia (Fig. 2a and Supplementary Fig. S2). To exclude that the above pattern may be dominated by strong diversity in mean precipitation and spatial variability across model, we further scale each model prior to taking the inter-model EOF analysis (see Methods), and obtain similar patterns compared with that in Fig. 2a (not shown). Hence, the synchronized precipitation changes across the entire AfroASM regions in the model spread imply that a large-scale controlling factor may play a dominant role.
Physical linkage between present-day bias and future projection uncertainty
The large-scale monsoon circulation is apparently driven by the thermal contrast between Northern Hemisphere (NH) and Southern Hemisphere (SH) due to moist static energy gradients associated with the seasonal swing of solar incidence39,40,41,42. Given the driving mechanism of monsoon, to understand the leading mode of projection uncertainty, we focus on the warming contrast between NH and SH that drives large-scale monsoon circulation. We regress surface temperature warming in 2050–2099 across models onto the normalized PC1 (Fig. 2b). A robust “NH warmer than SH” pattern is obtained, suggesting that models with larger increases of NH-SH ITC tend to project a wetter AfroASM, which is consistent with the basic driving mechanism of monsoon. In addition, the inter-model spread of NH-SH ITC in future projections correlates with the trend over the historical period (Fig. 2c). This pattern indicates that a model with a larger ITC trend in the present-day climate will project a greater increase of AfroASM precipitation in the future, as shown in the corresponding EOF1 (Fig. 2a). In addition, a remarkable warming anomaly related to PC1 in the historical period is seen over the Southern Ocean, which may be associated with the model biases in sea ice coverage43.
Why does the significant inter-model correlation between the present-day ITC trend and future AfroASM precipitation exist? A greater NH warming than SH is inherent to global warming given a smaller heat capacity due to a larger land area fraction in the NH and the Arctic amplification44,45,46. Hence, the inter-model scatter of NH-SH ITC is associated with that of global mean warming rate under a specific radiative forcing, i.e., the equilibrium climate sensitivity (ECS) which measures warming magnitude under doubled CO2 concentration relative to the pre-industrial period (Fig. 3a; see Methods). Models with a larger ECS show a larger NH-SH ITC in both the historical and future periods (Fig. 3b). A larger increase of ITC would lead to a stronger PC1 and thereby project more precipitation over AfroASM regions (Fig. 3c and Supplementary Fig. S1). The underlying physical mechanism is that a model with larger increase of ITC induces a stronger enhancement of low-level cross-equatorial flow over North Atlantic Ocean, Somalia and South China Sea (Supplementary Figs. S3a and S3b). The pattern of low-level cross-equatorial flow regressed onto the PC1 across models closely resembles that regressed onto the projected ITC (Supplementary Figs. S3c and S3d). Hence a larger projected increase of ITC would induce a stronger low-level cross-equatorial flow and thereby more moisture transport, finally resulting in more increase of AfroASM precipitation.
The underlying mechanism provides a solid physical basis to the observational constraint. Thus, we can constrain the future inter-model uncertainty of AfroASM precipitation based on the present-day observed ITC trend.
Constrained projection of AfroASM precipitation
The relationship between present-day warming patterns and projected spread across models allows the emergent constraint on the AfroASM precipitation, using multiple observation datasets. To measure the model’s fidelity in simulating the present-day observed climate, we define an ITC pattern index (ITCI). For each model, the ITCI is produced by projecting the present-day trend pattern of surface temperature onto the warming pattern associated with the inter-model PC1 shown in Fig. 2c. Using the above projecting pattern, the ITCI is defined as the difference between NH and SH (see Methods). The ITCI can well explain the leading mode of model uncertainty in projected AfroASM precipitation, as evinced by the significant correlation coefficient with the area mean of AfroASM precipitation (r = 0.58, p < 0.01; Fig. 4a) and the PC1 (r = 0.61, p < 0.01; Fig. 4b), respectively.
Based on the relationship between present-day climate (X) and projected PC1 (Y), we constrain the PC1 using a linear fit: \(Y=\bar{Y}+\rho (X-\bar{X})\), where ρ is the corrected regression coefficient, and \(\bar{X}\) and \(\bar{Y}\) represent the present-day and projected multi-model ensemble (MME), respectively (see Methods). To constrain the PC1, we firstly calculate the mean ITCI based on four observational datasets (vertical red dashed line in Fig. 4). The values of ITCI simulated by ~70% models are larger than the observations, indicating a systematic overestimation of the present-day ITC by CMIP6 models. The models with a larger ECS, such as CanESM5 and CanESM5-CanOE, simulate a stronger ITCI (Fig. 3a and Supplementary Table S3). The constrained AfroASM area-mean precipitation is 10 ± 6% (0.57 ± 0.38 mm day−1, constrained projection with ±1σ across models after emergent constraint). The constrained value of PC1 is −0.60 ± 0.80.
Since the regional precipitation changes are informative and crucial for climate change adaptation activities, we further constrain the spatial patterns in the projection. Based on the observational constraint of PC1 value and PC1-related patterns, we attain the corrected projections of precipitation and circulation patterns (Fig. 5; see Methods). The constrained projection indicates an increase of AfroASM precipitation by 10% in 2050–099 relative to 1965–2014 (Fig. 5a), which is ~70% of that of the raw (viz, unconstrained) projection. Correspondingly, the constrained monsoon circulation change is weaker than the raw multi-model ensemble (Supplementary Fig. S4). Locally, the constrained increases over West African, East Asian, and South Asian monsoon regions are 7%, 8%, and 12% respectively (Fig. 5a). The strongest reduction is seen over the West African monsoon region, where the constrained projection is only 49% of the raw MME projection, while over the East Asian monsoon region the constrained projection is 70% of the raw MME (Fig. 5b). The model uncertainty is also reduced after emergent constraint. Probability density function (PDF) of constrained PC1 is narrower than the original one, with a reduction of variance by 37% (Fig. 5c). Considering the explained variances of PC1 (26%), the total variance is reduced by ~10% (37% × 26%).
Impacts on the potential water availability
The AfroASM region holds a high density of population. More monsoon precipitation is expected to increase the potential water availability, which is mirrored in the runoff47,48, while the associated intense monsoon precipitation will also lead to flood and landslide49,50,51,52. The projected increase in monsoon precipitation under global warming is expected to partly offset the drying tendency since the 1950s (Fig. 1a). The smaller increase of precipitation in the constrained projection will reduce the increased potential water availability as expected from the raw projection, meanwhile the possible disasters related to heavy rainfall and floods will reduce accordingly. Here, we further estimate the impact of emergent constraint on the change in areal extent of precipitation and potential water availability.
To quantify the impact on the areal extent of precipitation, we examine the land area fraction that experiences a significant increase of precipitation (Fig. 6; see Methods). The fraction with a significant increase of precipitation is 24% in the constrained projection, only 57% of the raw projection. Regionally, in the constrained projection, the land area fraction in the East Asian monsoon region is only 37% of the raw projection, while in the West African and South Asian monsoon regions, the corresponding results are 50% and 69%, respectively.
Based on the significantly positive correlation between precipitation and runoff (Supplementary Fig. S5), we further quantify the changes of potential water availability in the constrained projection (Fig. 6; see Methods). About 27% land area in the AfroASM region will witness a significant increase of potential water availability in the constrained projection, which is only 66% of that of the raw projection. Regionally, the constrained land area fraction in the West African monsoon region is only 55% of the raw projection, while in the South Asia and East Asian monsoon regions, the corresponding result is 71% and 76%, respectively. Hence, less land area in AfroASM will experience a significant increase of precipitation and runoff in the constrained projection. These imply that the characteristics of precipitation change in the future will be milder than the raw projection.
Discussion
Our emergent constraint on the projection of monsoon precipitation in the AfroASM regions is based on SSP5-8.5 scenario. We extend our analysis to the medium emission scenario of SSP2-4.5 and get similar conclusion (Supplementary Figs. S6 and S7). In addition, to examine the robustness of the emergent constraint, we check the inter-model correlation coefficient between ITCI and PC1 by using different subsets of model and including randomized outliers, and come to similar conclusion (Supplementary Fig. S8). The independence of the results on the model ensemble and the future emission scenarios confirms the robustness of the conclusions.
Our results reveal that the raw projection overestimates the increase of precipitation in the AfroASM region. While the constrained projection of the increase in AfroASM precipitation is 70% of the raw projection in the context of regional average, the effects of emergent constraint on the changes of precipitation and water availability are more pronounced at regional scales. The projection of precipitation constrained by the observation is 49% (70%) of the raw projection in the West Africa (East Asia) monsoon region, even reduced by 70% (50%) more widely over sub-regional or local scales. The land fraction that will experience a significant increase of precipitation is 50% (37%) of that of the raw projection in West Africa (East Asia) monsoon region. The change of precipitation is echoed in the runoff as an indicator of potential water availability and the risk of flood. The smaller increase of potential water availability than the raw projection may pose a challenge to climate change adaptation and mitigation activities related to water management and food security53,54, although a smaller than expected increase in rainfall will also reduce the risk of extreme precipitation and flooding.
Given that the inter-model uncertainty of global mean warming is closely associated with the inter-model spread of ECS55,56, with normalizing by the global mean warming, the precipitation response (viz, hydrological sensitivity) still shows a remarkable spread across models10,12,57. A recent study reported that the projected uncertainty of hydrological sensitivity over the AfroASM regions is related to the projected uncertainty of ITC and land-sea thermal contrast12. Since the inter-model uncertainties of both hydrological sensitivity and global mean surface air temperature (GSAT) are significantly correlated with the ITC (Supplementary Fig. S9), we further constrain the hydrological sensitivity and GSAT separately based on the intermodel relationship between projected uncertainty and present-day biases (see Methods). The results based on constraining hydrological sensitivity and GSAT separately are consistent with that based on constraining the precipitation changes directly. The constraining effects from hydrological sensitivity and GSAT reduce the raw precipitation projection by 13% and 21%, respectively. Thus, the GSAT warming plays a dominant role in the emergent constraint on precipitation changes, while the contribution from hydrological sensitivity should not be neglected.
Methods
Observations
To correct the present-day model biases and constrain the future projections, multiple observational datasets are used. The monthly gridded observational surface temperature datasets used are (1) Berkeley Earth Surface Temperature (BEST)58, (2) Cowtan and Way version 259, (3) NASA Goddard Institute for Space Studies Surface Temperature version 5 (GISTEMP v5)60, and (4) NOAA Global Surface Temperature version 5.0 (NOAAGlobalTemp v5)61,62. The observational surface temperature datasets have been regridded on 2.5° × 2.5°.
Precipitation datasets are taken from (1) Climatic Research Unit (CRU) Time-Series version 4.0263, (2) Global Precipitation Climatology Centre version 7 (GPCC v7)64, (3) Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP v1201)65, and (4) the Global Precipitation Climatology Project (GPCP v2.3)66. Details of the datasets are provided in Supplementary Table 2.
Model simulations
We use the monthly output from 30 CMIP6 models (Supplementary Table S1) in historical simulation, and future projection under SSP5-8.5 scenario35,37. Given the recent trends in decarbonization, SSP5-8.5 is a highly unlikely scenario38. Hence, the monthly data from 18 CMIP6 models in future projection under median emission scenario, SSP2-4.5, are also used to verify the independence of emission scenarios. The first available realization for each model is used to give equal weight to each model. To examine the contribution from internal variability, we use the output of 29 CMIP6 models from the piControl simulations, in which all external forcings are held constant at their 1850 levels (Supplementary Table S1). All the data is re-gridded to 2.5° × 2.5° grids using first-order conservative interpolation, except for the circulation patterns which is re-gridded using bilinear interpolation.
To represent future projection of the summer monsoon rainfall, we focus on the changes in seasonal mean precipitation from June to September (JJAS) in 2050–2099 relative to the mean in 1965–2014.
Afro-Asian monsoon
The Afro-Asian monsoon region is defined as the land monsoon area over Eurasian continent and North Africa3,67 (Fig. 1). The land monsoon domain is defined as the land area where the precipitation difference between the local summer and winter is larger than 2.0 mm day−1, and local summer precipitation exceeds 55% of the annual total precipitation67, based on the climatological mean of CMAP65 and GPCP66. Local summer is defined as May to September for the Northern Hemisphere (NH). The Afro-Asian monsoon consists of 3 regional monsoons, including East Asia, South Asia, and North Africa monsoon (Fig. 1b, c contour).
A strong monsoon circulation is marked by a strengthening of the vertical zonal wind shear68,69. To quantify the response of AfroASM circulation, we define a circulation index as the vertical shear of zonal winds between 850 and 200 hPa averaged in a zone stretching from North Atlantic eastward to the Philippines (0–20°N, 30°W–120°E).
Inter-model empirical orthogonal function (EOF) analysis
The leading modes of inter-model uncertainty in the summer precipitation projection over Afro-Asian monsoon region is obtained by applying the typical EOF method to model-spatial dimension:
in which ∆ denotes projected changes, m denotes model number, n denotes the spatial area, and num is the mode number. Prime represents the deviation from the multi-model ensemble (MME). PCs are normalized here. The inter-model EOF method has been successfully used for the East Asia monsoon29,31, tropical ocean SST70,71,72, and extra-tropical oceans SST73.
The leading principal component (PC1) of AfroASM precipitation accounts for 26% of the total intermodal variance (Fig. 2).
Scaling individual models
To confirm the results are not dependent on the strength of individual models’ hydrological cycles, the mean precipitation changes over tropics and global for each model have been removed from the original precipitation changes, respectively, and then the precipitation changes have been normalized by the corresponding spatial standard deviation. We take the inter-model EOF analysis for the scaled precipitation changes. The results show that the patterns of projected changes of precipitation and low-level circulation regressed onto the scaled PC1 across models closely resemble that in Fig. 2a. Thus, we only present the results without scaling in the paper.
The definition of interhemispheric thermal contrast (ITC) and ITC pattern index
The future increase of AfroASM precipitation is closely associated with the projected NH-SH thermal contrast, with correlation coefficient higher than 0.7 across models12,13,32. To constrain the future projection using the present-day observation, we select two key regions which represent the NH-SH thermal contrast to define ITC and ITC pattern index (ITCI).
The ITC is defined as the difference of the area-averaged surface temperature between (20°N~50°N, 0~360°) and (20°S~50°S, 0~360°). To represent the pattern of ITC, a pattern index is produced by projecting surface temperature trend onto the pattern associated with PC1 (TPC1) shown in Fig. 2c.
The period 1981–2014 is chosen to calculate the pattern indices and constrain the projection. The warming trend of ITC in 1981–2014 has been dominated by the response to greenhouse gases, with no significant trend in aerosol cooling74,75,76,77,78.
For the ITCI of each model, to clearly reflect the present-day warming pattern, the historical warming trend of surface temperature in each model (THist) is projected onto the inter-model warming trend shown in Fig. 2c (TPC1) in Northern Hemisphere (NH; 20°N~50°N, 0~360°) and Southern Hemisphere (SH; 20°S~50°S, 0~360°), respectively, following Chen31:
where 〈〉 denotes area mean.
To calculate the index in observation, THist is derived from the four observational surface temperature datasets (Supplementary Table S2).
Equilibrium climate sensitivity
To investigate the source of model biases of the present-day ITC trend, we use the model’s equilibrium climate sensitivity (ECS). The ECS is represented by the effective climate sensitivity which is estimated by regressing the net top-of-atmosphere radiance against the global mean surface air temperature changes in the first 150 years of the CO2 quadruples experiment75,78,79,80. The ECS of most models in this study is derived from the Table 2 in Meehl80, except for that of CanESM5-CanOE which is derived from Swart81, and FGOALS-g3 which is derived from Zhou82 and Li83.
The contribution from the unforced internal variability
To quantify the impact from the unforced internal variability to the ITCI, we calculate the ITCI based on the piControl simulations of 29 CMIP6 models. We calculate the trend of randomly selected continuous 34-year period and repeat this process over 1000 times to obtain 1000 synthetic members from each model. The contribution of the internal variability is measure by the variance across different synthetic members \(({{\sigma }}_{{internal}}^{2})\), which is approximately 7.8 × 10−4 K2 34 yr−2, compared to 3.9 × 10−3 K2 34 yr−2 of present-day ITCI across 30 models (\({\sigma }_{{intermodel}}^{2}\)). The contribution from the internal variability accounts for 20% based on the variance ratio (\({\sigma }_{{internal}}^{2}/{\sigma }_{{intermodel}}^{2}\)). The range of unforced internal variability in Fig. 4 (light gray shading) is represented by ±1σ across different synthetic members (\({\sigma }_{{internal}}^{2}\)).
Hierarchical statistical framework for emergent constraint
To constrain the projected AfroASM precipitation, we use the hierarchical emergent constraint (HEC) framework proposed by Bowman24. The HEC framework accounts for both the correlation between future and present-day climate, and the precision in the observational datasets, compared with the classical emergent constraint24,31.
In the HEC framework, we establish a link between future climate change (Y) and present-day climate (X) to constrain Y. The emergent constraint is based on the linear regression between Y and X obtained from climate models:
where ρi is the regression coefficient. \(\bar{X}\) and \(\bar{Y}\) are the multi-model ensemble mean of X and Y, respectively. Y is AfroASM precipitation changes or PC1, and X is ITCI in Eq. (2).
Since we use the observation in current climate (XO) to constrain Y, the uncertainty in the observations should be considered. Under the Gaussian assumptions which relates the observations to current climate24, the signal-noise ratio (SNR) in the observation is the ratio between the variance across models (\({\sigma }_{X}^{2}\)) and observational datasets (\({\sigma }_{O}^{2}\)):
The regression coefficient is multiplied by \(\frac{1}{1+{SN}{R}^{-1}}\) to correct ρi31. If the SNR is much larger than 1 (SNR » 1), the effect of correction can be neglected. In our analysis, SNR of ITCI is 110.
Based on the Eqs. (3) and (4), the constrained results and variance of future climate change \(\bar{{Y}_{C}}\) can be expressed as:
Where ρ is the corrected regression coefficient, i.e., \(\frac{{\rho }_{i}}{1+{SN}{R}^{-1}}\); r is the corrected correlation coefficient between X and Y, i.e., \(\frac{{r}_{i}^{2}}{1+{SN}{R}^{-1}}\), and the ri is the original correlation coefficient between X and Y.
Based on the Eq. (6), the relative variance reduction (\(1-\frac{{\sigma }_{{Y}_{C}}^{2}}{{\sigma }_{Y}^{2}}\)) from the HEC framework is \(\frac{{r}_{i}^{2}}{1+{SN}{R}^{-1}}\). The total reduced variance (TRV) after constraining the PC1 can be expressed as the weighting on the corresponding explained variances of PC1 (PCV1):
To constrain the projection results of individual model, we correct the PC1 of individual model using the \(\bar{{X}_{O}}\). Base on the relationship of emergent constraint, the constrained PC1 of each model (\(\bar{{Y}_{m,C}}\)) can be expressed as:
where Xm and Ym is the ITCI and the PC1 of each model, respectively. The inter-model standard deviation of the constrained PC1 is 0.80.
Corrected multi-model mean projection
The corrected PC1 are estimated by the emergent constraint using the observed ITCI. Since the pattern scaling has been shown to work robustly for seasonal averages of precipitation22, the projection of the AfroASM precipitation can be corrected based on EOF reconstruction, following Eq. (1):
where subscript “O” denotes the corrected PC1 constrained by observation, \(\triangle \overline{{\Pr }}\) represents the multi-model ensemble mean of precipitation, and \({{\Pr }}_{{PC}1}^{{\prime} }\) represents the PC1-related pattern of precipitation changes shown in Fig. 2a. The wind fields at 850 hPa can be corrected in a similar way of Eq. (9), but the \({{\Pr }}_{{PC}1}^{{\prime} }\) term need to be replaced by the regression coefficients related to PC1.
Spatial aggregated probability density function
To measure the area with a certain change, we calculate the spatial aggregated probability density function (PDF) in the projected changes of precipitation. According to the latitude-dependent area, the grid points falling in each bin of the PDF have been weighted. Hence, the spatial distribution is an aggregated of all grid area satisfying the conditional sampling. The PDF is derived from the nonparametric assessment of the PDF. The spatial PDF is proposed by Fischer84 and successfully used for the detection of extreme climate events85,86,87,88.
Land fraction that experiences a significant increase of precipitation
The significant increase in the spatial aggregated PDF over AfroASM and three submonsoon regions is defined as the increase exceeds the inter-model standard deviation over AfroASM and three submonsoon regions, respectively. The area that experiences significant increases is aggregated spatially to represent the area which witness a significant increase of precipitation. Fraction is calculated with respect to the total area.
Impact of precipitation change on runoff projection
For each model, we calculate the projected changes in summer mean runoff for the period of 2050–2099 relative to the baseline. For the entire AfroASM and each submonsoon regions, there is a strong linear correlation between runoff changes and precipitation changes across models (Supplementary Fig. S5). This relationship enables a constraint on future runoff changes by using observationally constrained precipitation changes. Hence, following previous studies48, the constrained changes of runoff (\(\triangle {R}_{{constrained}}\)) are derived as following equation:
where ∆Prconstrained denotes the constrained change of precipitation based on Eq. (10), and k denotes the regression coefficient between changes of runoff and precipitation, and b denotes intercept.
Constrained projection of hydrological sensitivity and GSAT
The hydrological sensitivity is defined as the precipitation response normalized by the GSAT warming in 2050–2099 relative to 1965–2014. Since the inter-model spread of hydrological sensitivity is closely related to that of ITC (Supplementary Fig. S9a), which is consistent with recent study12, we constrain the projection of hydrological sensitivity over AfroASM region based on Eq. (5) and Eq. (8). The constrained response of hydrological sensitivity is only 87% of that of the raw projection.
We constrained the projected GSAT warming using the observed GSAT trend in 1981–2014, following Tokarska75 and Lee45. The constrained GSAT warming is 2.62 ± 0.57 K, weaker than the raw projection (3.30 ± 0.78 K) under SSP5-8.5.
Data availability
The data that support the findings of this study are freely available. CMIP6 model data are from the Earth System Grid Federation [https://esgf-node.llnl.gov/search/cmip6/]. Observational temperature BEST is from the Berkeley Earth [http://berkeleyearth.org/data-new/], Cowtan and Way v2 is from the University of York and the University of Ottawa [https://www-users.york.ac.uk/~kdc3/papers/coverage2013/series.html], GISTEMP is from the NASA GISS [https://data.giss.nasa.gov/gistemp/]. Observational precipitation CRU is from the University of East Anglia [http://badc.nerc.ac.uk/data/cru/]. NOAAGlobalTemp v5, GPCC v7, CMAP v1201 and GPCP v2.2 are provided by the NOAA/OAR/ESRL, PSD, Boulder, CO, USA [https://psl.noaa.gov/data/gridded/].
Code availability
The data in this study is analyzed with NCAR Command Language (NCL; http://www.ncl.ucar.edu/). The relevant codes in this work are available, upon request, from the corresponding author T. Z.
References
Stager, J. C., Ryves, D. B., Chase, B. M. & Pausata, F. S. R. Catastrophic drought in the Afro-Asian monsoon region during Heinrich event 1. Science 331, 1299–1302 (2011).
Ji, J., Petit-Maire, N. & Yan, Z. The last 1000 years: climatic change in arid Asia and Africa. Glob. Planet. Change 7, 203–210 (1993).
Li, Y., Ding, Y. & Li, W. Interdecadal variability of the Afro-Asian summer monsoon system. Adv. Atmos. Sci. 34, 833–846 (2017).
He, C., Li, T. & Zhou, W. Drier North American monsoon in contrast to Asian-African monsoon under global warming. J. Clim. 1–47, https://doi.org/10.1175/jcli-d-20-0189.1 (2020).
A. Cherchi, A. Turner. IPCC. Annex V: Monsoons. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds. Masson-Delmotte, V. et al.) (Cambridge University Press, 2021). https://www.ipcc.ch/report/ar6/wg1/downloads/report/IPCC_AR6_WGI_Annex_V.pdf.
Zhang, W., Zhou, T., Zou, L., Zhang, L. & Chen, X. Reduced exposure to extreme precipitation from 0.5 °C less warming in global land monsoon regions. Nat. Commun. 9, 1–8 (2018).
Wang, B. et al. Monsoon climate change assessment. Bull. Am. Meteorol. Soc. https://doi.org/10.1175/BAMS-D-19-0335.1 (2020).
Zhou, T., Lu, J., Zhang, W. & Chen, Z. The sources of uncertainty in the projection of global land monsoon precipitation. Geophys. Res. Lett. 47, https://doi.org/10.1029/2020gl088415 (2020).
Park, J. Y., Bader, J. & Matei, D. Northern-hemispheric differential warming is the key to understanding the discrepancies in the projected Sahel rainfall. Nat. Commun. 6, 1–8 (2015).
Li, G., Xie, S.-P., He, C. & Chen, Z. Western Pacific emergent constraint lowers projected increase in Indian summer monsoon rainfall. Nat. Clim. Chang. 7, 708–712 (2017).
Chen, Z. et al. Global land monsoon precipitation changes in CMIP6 projections. Geophys. Res. Lett. 47, https://doi.org/10.1029/2019gl086902 (2020).
Cao, J. et al. Sources of the inter‐model spread in projected global monsoon hydrological sensitivity. Geophys. Res. Lett. 47, https://doi.org/10.1029/2020gl089560 (2020).
Wang, B., Jin, C. & Liu, J. Understanding future change of global monsoon projected by CMIP6 models. J. Clim. 33, 6471–6489 (2020).
Jin, C., Wang, B. & Liu, J. Future changes and controlling factors of the eight regional monsoons projected by CMIP6 models. J. Clim. 33, 9307–9326 (2020).
Ha, K.-J., Moon, S., Timmermann, A. & Kim, D. Future changes of summer monsoon characteristics and evaporative demand over Asia in CMIP6 simulations. Geophys. Res. Lett. https://doi.org/10.1029/2020GL087492 (2020).
Moon, S. & Ha, K.-J. Future changes in monsoon duration and precipitation using CMIP6. npj Clim. Atmos. Sci. 3, https://doi.org/10.1038/s41612-020-00151-w (2020).
Cox, P. & Stephenson, D. A changing climate for prediction. Science 317, 207–208 (2007).
Tebaldi, C. & Knutti, R. The use of the multi-model ensemble in probabilistic climate projections. Science 365, 2053–2075 (2007).
Hawkins, E. & Sutton, R. The potential to narrow uncertainty in regional climate predictions. Bull. Am. Meteorol. Soc. 90, 1095–1108 (2009).
Kirtman, B. et al. Near-term climate change: projections and predictability. In Climate Change 2013 the Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds. Stocker, T. F. et al.) vol. 9781107057 953–1028, (Cambridge University Press, 2013).
Brunner, L. et al. Reduced global warming from CMIP6 projections when weighting models by performance and independence. Earth Syst. Dyn. 1–24, https://doi.org/10.5194/esd-2020-23 (2020).
Tebaldi, C. & Arblaster, J. M. Pattern scaling: Its strengths and limitations, and an update on the latest model simulations. Clim. Change 122, 459–471 (2014).
Xie, S.-P. et al. Towards predictive understanding of regional climate change. Nat. Clim. Chang. 5, 921–930 (2015).
Bowman, K. W., Cressie, N., Qu, X. & Hall, A. A hierarchical statistical framework for emergent constraints: application to snow‐albedo feedback. Geophys. Res. Lett. 45, https://doi.org/10.1029/2018gl080082 (2018).
Brient, F. Reducing uncertainties in climate projections with emergent constraints: concepts, examples and prospects. Adv. Atmos. Sci. 37, 1–15 (2020).
Hall, A., Cox, P., Huntingford, C. & Klein, S. Progressing emergent constraints on future climate change. Nat. Clim. Chang. 9, 269–278 (2019).
Caldwell, P. M., Zelinka, M. D. & Klein, S. A. Evaluating emergent constraints on equilibrium climate sensitivity. J. Clim. 31, 3921–3942 (2018).
Klein, S. A. & Hall, A. Emergent constraints for cloud feedbacks. Curr. Clim. Chang. Rep. 1, 276–287 (2015).
Zhou, S., Huang, P., Huang, G. & Hu, K. Leading source and constraint on the systematic spread of the changes in East Asian and western North Pacific summer monsoon. Environ. Res. Lett. 14, 124059 (2019).
Yan, Y., Lu, R. & Li, C. Relationship between the future projections of Sahel rainfall and the simulation biases of present South Asian and Western North Pacific rainfall in summer. J. Clim. 32, 1327–1343 (2019).
Chen, X., Zhou, T., Wu, P., Guo, Z. & Wang, M. Emergent constraints on future projections of the western North Pacific Subtropical High. Nat. Commun. 11, https://doi.org/10.1038/s41467-020-16631-9 (2020).
Endo, H., Kitoh, A. & Ueda, H. A unique feature of the Asian summer monsoon response to global warming: the role of different land–sea thermal contrast change between the lower and upper troposphere. SOLA 14, 57–63 (2018).
Wang, P. X. et al. The global monsoon across timescales: coherent variability of regional monsoons. Clim. Past 10, 2007–2052 (2014).
Feudale, L. & Kucharski, F. A common mode of variability of African and Indian monsoon rainfall at decadal timescale. Clim. Dyn. 41, 243–254 (2013).
Eyring, V. et al. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev. 9, 1937–1958 (2016).
Riahi, K. et al. The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: an overview. Glob. Environ. Chang. 42, 153–168 (2017).
O’Neill, B. C. et al. The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev. 9, 3461–3482 (2016).
Hausfather, Z. & Peters, G. P. Emissions – the ‘business as usual’ story is misleading. Nature 577, 618–620 (2020).
Trenberth, K. E., Stepaniak, D. P. & Caron, J. M. The global monsoon as seen through the divergent atmospheric circulation. J. Clim. 13, 3969–3993 (2000).
Zhou, W. & Xie, S. A hierarchy of idealized monsoons in an intermediate GCM. J. Clim. 31, 9021–9036 (2018).
Hill, S. A. Theories for past and future monsoon rainfall changes. Curr. Clim. Chang. Rep. 5, 160–171 (2019).
Geen, R., Bordoni, S., Battisti, D. S. & Hui, K. L. Monsoons, ITCZs and the concept of the global monsoon. Rev. Geophys. 58, 1–60 (2020).
Shi, X., Chen, X., Dai, Y. & Hu, G. Climate sensitivity and feedbacks of BCC-CSM to idealized CO2 forcing from CMIP5 to CMIP6. J. Meteorol. Res. 34, 865–878 (2020).
Byrne, M. P. & O’Gorman, P. A. Land-ocean warming contrast over a wide range of climates: Convective quasi-equilibrium theory and idealized simulations. J. Clim. 26, 4000–4016 (2013).
Lee, J. Y. et al. Chapter 4: Future global climate: scenario-based projections and near-term information. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds. Masson-Delmotte, V. et al.) (Cambridge University Press, 2021).
Previdi, M., Smith, K. L. & Polvani, L. M. Arctic amplification of climate change: a review of underlying mechanisms. Environ. Res. Lett. 16, https://doi.org/10.1088/1748-9326/ac1c29 (2021).
Schewe, J. et al. Multimodel assessment of water scarcity under climate change. Proc. Natl Acad. Sci. USA 111, 3245–3250 (2014).
Wang, T. et al. Atmospheric dynamic constraints on Tibetan Plateau freshwater under Paris climate targets. Nat. Clim. Chang. 11, 219–225 (2021).
Cho, C., Li, R., Wang, S. Y., Yoon, J. H. & Gillies, R. R. Anthropogenic footprint of climate change in the June 2013 northern India flood. Clim. Dyn. 46, 797–805 (2016).
Kundzewicz, Z. W. et al. Flood risk and its reduction in China. Adv. Water Resour. 130, 37–45 (2019).
Li, C., Chai, Y., Yang, L. & Li, H. Spatio-temporal distribution of flood disasters and analysis of influencing factors in Africa. Nat. Hazards 82, 721–731 (2016).
Ranasinghe, R. et al. Chapter 12: Climate change information for regional impact and for risk assessment. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds. Masson-Delmotte, V. et al.) (Cambridge University Press, 2021).
Trenberth, K. E. & Asrar, G. R. Challenges and opportunities in water cycle research: WCRP Contributions. Surv. Geophys. 35, 515–532 (2014).
Mbow, C., et al. Food security. In Climate Change and Land: an IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems (eds. Shukla, P. R. et al.) 437–550 (2019).
Cox, P. M., Huntingford, C. & Williamson, M. S. Emergent constraint on equilibrium climate sensitivity from global temperature variability. Nature 553, 319–322 (2018).
Sherwood, S. C. et al. An assessment of Earth’s climate sensitivity using multiple lines of evidence. Rev. Geophys. 58, 1–92 (2020).
Allen, M. & Ingram, W. J. Constraints on future changes in climate and the hydrologic cycle. Nature 489, 224–232 (2002).
Rohde, R., Muller, R., Jacobsen, R., Perlmutter, S. & Mosher, S. Berkeley Earth temperature averaging process. Geoinformatics Geostatistics An Overv. 01, https://doi.org/10.4172/2327-4581.1000103 (2013).
Cowtan, K. & Way, R. G. Coverage bias in the HadCRUT4 temperature series and its impact on recent temperature trends. Q. J. R. Meteorol. Soc. 140, 1935–1944 (2014).
Lenssen, N. et al. Improvements in the GISTEMP uncertainty model. J. Geophys. Res. Atmos. 124, 6307–6326 (2019).
Vose, R. S. et al. NOAA’s merged land-ocean surface temperature analysis. Bull. Am. Meteorol. Soc. 93, 1677–1685 (2012).
Zhang, H.-M. et al. Updated temperature data give a sharper view of climate trends. Eos. 100, https://doi.org/10.1029/2019EO128229 (2019).
Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations - the CRU TS3.10 Dataset. Int. J. Climatol. 34, 623–642 (2014).
Schneider, U. et al. GPCC’s new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle. Theor. Appl. Climatol. 115, 15–40 (2013).
Xie, P. & Arkin, P. A. Global precipitation: a 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bull. Am. Meteorol. Soc. 78, 2539–2558 (1997).
Adler, R. F. et al. The version-2 global precipitation climatology project (GPCP) monthly precipitation analysis (1979-present). J. Hydrometeorol. 4, 1147–1167 (2003).
Wang, B., Liu, J., Kim, H.-J., Webster, P. J. & Yim, S.-Y. Recent change of the global monsoon precipitation (1979–2008). Clim. Dyn. 39, 1123–1135 (2012).
Webster, P. J. & Yang, S. Monsoon and ENSO: selectively interactive systems. Q. J. R. Meteorol. Soc. 118, 877–926 (1992).
Wang, B. et al. Northern Hemisphere summer monsoon intensified by mega-El Nino/southern oscillation and Atlantic multidecadal oscillation. Proc. Natl Acad. Sci. USA 110, 5347–5352 (2013).
Li, G. & Xie, S. P. Tropical biases in CMIP5 multimodel ensemble: The excessive equatorial Pacific cold tongue and double ITCZ problems. J. Clim. 27, 1765–1780 (2014).
Li, G. & Xie, S. P. Origins of tropical-wide SST biases in CMIP multi-model ensembles. Geophys. Res. Lett. 39, L22703 (2012).
Li, G., Xie, S.-P. & Du, Y. Monsoon-induced biases of climate models over the tropical Indian Ocean. J. Clim. 28, 3058–3072 (2015).
Wang, C., Zhang, L., Lee, S.-K., Wu, L. & Mechoso, C. R. A global perspective on CMIP5 climate model biases. Nat. Clim. Chang. 4, 201–205 (2014).
Jiménez-de-la-Cuesta, D. & Mauritsen, T. Emergent constraints on Earth’s transient and equilibrium response to doubled CO2 from post-1970s global warming. Nat. Geosci. 12, 902–905 (2019).
Tokarska, K. B. et al. Past warming trend constrains future warming in CMIP6 models. Sci. Adv. 6, https://doi.org/10.1126/sciadv.aaz9549 (2020).
Liang, Y., Gillett, N. P. & Monahan, A. H. Climate model projections of 21st century global warming constrained using the observed warming trend. Geophys. Res. Lett. 47, https://doi.org/10.1029/2019gl086757 (2020).
Nijsse, F. J. M. M., Cox, P. M., Huntingford, C. & Williamson, M. S. Decadal global temperature variability increases strongly with climate sensitivity. Nat. Clim. Chang. 9, 598–601 (2019).
Forster, P. M. et al. Evaluating adjusted forcing and model spread for historical and future scenarios in the CMIP5 generation of climate models. J. Geophys. Res. Atmos. 118, 1139–1150 (2013).
Gregory, J. M. et al. A new method for diagnosing radiative forcing and climate sensitivity. Geophys. Res. Lett. 31, https://doi.org/10.1029/2003gl018747 (2004).
Meehl, G. A. et al. Context for interpreting equilibrium climate sensitivity and transient climate response from the CMIP6 Earth system models. Sci. Adv. 6, eaba1981 (2020).
Swart, N. C. et al. The Canadian Earth System Model version 5 (CanESM5.0.3). Geosci. Model Dev. 12, 4823–4873 (2019).
Zhou, T. et al. Development of climate and earth system models in China: past achievements and new CMIP6 results. J. Meteorol. Res. 34, 1–19 (2020).
Li, L. et al. The flexible global ocean‐atmosphere‐land system model grid‐point version 3 (FGOALS‐g3): description and evaluation. J. Adv. Model. Earth Syst. 12, https://doi.org/10.1029/2019ms002012 (2020).
Fischer, E. M., Beyerle, U. & Knutti, R. Robust spatially aggregated projections of climate extremes. Nat. Clim. Chang. 3, 1033–1038 (2013).
Fischer, E. M. & Knutti, R. Detection of spatially aggregated changes in temperature and precipitation extremes. Geophys. Res. Lett. 41, 547–554 (2014).
Frolicher, T. L., Fischer, E. M. & Gruber, N. Marine heatwaves under global warming. Nature 560, 360–364 (2018).
Zhao, S. & Zhou, T. Are the observed changes in heat extremes associated with a half‐degree warming increment analogues for future projections? Earth’s Futur. 7, 978–992 (2019).
Zhao, S., Zhou, T. & Chen, X. Consistency of extreme temperature changes in China under a historical half-degree warming increment across different reanalysis and observational datasets. Clim. Dyn. 54, 2465–2479 (2020).
Acknowledgements
This study is supported by the National Key Research and Development Program of China (2020YFA0608904), National Natural Science Foundation of China (Grant No. 41988101), and National Key Research and Development Program of China (2018YFA0606501). We thank the World Climate Research Programme’s Working Group on Coupled Modeling, which is responsible for CMIP6, and the climate modeling groups (listed in Supplementary Table S1) for producing and making available their model output (https://esgf-node.llnl.gov/search/cmip6/). We also acknowledge NOAA/OAR/ESRL PSD for providing the observational precipitation datasets (Global Precipitation Climatology Project (GPCP v2.3): https://psl.noaa.gov/data/gridded/data.gpcp.html; Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP v1201): https://psl.noaa.gov/data/gridded/data.cmap.html; Global Precipitation Climatology Centre version 7 (GPCC v7): https://www.esrl.noaa.gov/psd/data/gridded/data.gpcc.html) and the University of East Anglia Climatic Research Unit (CRU) for providing the observational precipitation datasets (https://catalogue.ceda.ac.uk/uuid/10d3e3640f004c578403419aac167d82).
Author information
Authors and Affiliations
Contributions
T.Z. designed the research, provided comments, and revised the manuscript. Z.C. performed the analysis and drafted the manuscript. X.C., W.Z., Lixia Z., Liwei Z. and M.W. provided comments, and helped to organize and revise the draft. All of the co-authors contributed to scientific interpretations and helped to improve the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chen, Z., Zhou, T., Chen, X. et al. Observationally constrained projection of Afro-Asian monsoon precipitation. Nat Commun 13, 2552 (2022). https://doi.org/10.1038/s41467-022-30106-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-30106-z
This article is cited by
-
Historical footprints and future projections of global dust burden from bias-corrected CMIP6 models
npj Climate and Atmospheric Science (2024)
-
Robust future intensification of winter precipitation over the United States
npj Climate and Atmospheric Science (2024)
-
Impact of spring land-surface conditions over the Tibetan Plateau on the early summer Asian monsoon using an AGCM large ensemble
Climate Dynamics (2024)
-
Evaluation and Projection of Extreme High Temperature Indices in Southwestern China Using NEX-GDDP-CMIP6
Journal of Meteorological Research (2024)
-
Future climate change in the Northern Indian Ocean as simulated with a high-resolution regional earth system model
Climate Dynamics (2024)