Anthropogenic warming disrupts intraseasonal monsoon stages and brings dry-get-wetter climate in future East Asia

East Asia will face a skewed monsoon cycle with soaring flood, drought, and weather whiplash risks in a warming climate. In our objective eight-intraseasonal-monsoon-stage framework, we uncover a ‘dry-get-wetter’ paradigm in East Asia, contesting the fallen ‘rich-get-richer’ common belief. On timing, the Mid-summer and Fall periods are stretching at the expense of three delayed, shortened, and weakened winter stages, especially near the end of the twenty-first century. On threats, entire East Asia will experience up to 14–20 more heavy precipitation days during the rainy Spring to Mid-summer stages. Specifically, the Yangtze basin will suffer from an earlier pluvial period with escalating flood risks. Moreover, societal security and ecosystem resilience in the Huai-Yellow basin, South Japan, and the Korean Peninsula will be challenged by more frequent weather whiplash. Under the monsoon-stage framework, a complete moisture budget decomposition sheds light on the causes of a slower precipitation scaling and the ‘dry-get-wetter’ paradigm.


INTRODUCTION
The East Asian monsoon (EAM) is a subtropical monsoon system north of 20˚N featuring a seasonal reversal of winds and a contrast between rainy and dry seasons 1,2 . On a large scale, it is fundamentally driven by seasonally-varying solar radiations, shaped by the land-sea thermal contrast and upper-level westerly jets due to Tibetan Plateau forcing 1,3-6 . In recent decades, growing concerns over future changes of the EAM in a warming climate are warranted given its profound influence on the livelihoods of more than 1.5 billion people, plus the large uncertainties presented by climate model projections of monsoon precipitation, circulations, and seasonality [7][8][9][10][11] .
Response of the EAM's seasonality to future radiative forcing has been analyzed commonly using precipitation-based monsoon indices [10][11][12][13] . Multiple lines of evidence suggested that the future wet season in East Asia seemed stretched with a remarkably delayed retreat [10][11][12][13][14] and at the expense of the dry duration 15 . While most studies focused on seasonal changes in the future [10][11][12][13][14] , there remains a dearth of research on the future evolution of intraseasonal monsoon circulations and rainband featured by the EAM system 1,16 , which are more relevant to climate-and weather-related activities. For this reason, we follow a recently developed circulation-based clustering approach to derive the EAM annual cycle 17 . The monsoon cycle consists of eight intraseasonal monsoon stages, among which the Early-, Mid-, and Late-winter stages manifest in anticyclonic circulations and widespread dry-and-cold air masses in East Asia, which are most dominant in the Mid-winter. Following the Late-winter, low-level southwesterlies become prominent in the south during the Spring and a quasi-stationary rain belt is then formed in South China during the Pre-Meiyu. The summertime southwesterlies peak in the ensuing Meiyu season, alongside the strongest monsoon rain belts lingering over the Yangtze basin and South Japan. The Midsummer sees a continent-wide weakening of summer monsoons followed by the rise of anticyclonic circulations in the Fall, which completes the annual monsoon cycle. A multi-model weighting scheme is proposed to integrate the climate model outputs based on their competence in reproducing the observed EAM annual cycle (see 'Methods'). Using such a framework, projected changes in the monsoon onset, retreat, and precipitation will facilitate more specific and comprehensive policies in optimizing irrigation schemes, sowing and harvesting, reservoir management, and climate adaptations in each intraseasonal stage 18 .
Another moot point lies in the causes of the climate sensitivity of monsoon precipitation to greenhouse warming 9,[19][20][21][22][23][24] . While water vapors in the atmosphere largely follow the Clausius-Clapeyron (C-C) rate (i.e.,~7% K −1 ) 12,19,21,[25][26][27] , future precipitation change is likely to scale at a much slower rate than that, suggesting additional processes at work to offset the moistening of the atmosphere. Scientists certified that the thermodynamic effect could be counteracted by an overall slowdown of atmospheric circulations 9,19,20,24,28,29 , plausibly due to an increased atmospheric static stability in both lower and upper troposphere 14,30,31 . Paradoxical as it may seem, some studies indicated that East Asia might expect both enhanced precipitation and strengthened monsoon circulation manifested by an increase in lower-level and upper-level winds 32 as well as ascending motions 21 . Apparently, further investigation on the attributions to the future precipitation change in East Asia is much needed. Seeing that the monsoon climate in East Asia features pronounced variabilities within seasons, we show that the governing physical processes could be markedly different in well-defined intraseasonal stages, offering additional insights into the causes of precipitation changes under anthropogenic warming, and benefiting climate model evaluation and projection.
Here we propose a framework to assess the future monsoon cycle that fills the knowledge gap in the projected changes of intraseasonal monsoon stages in East Asia under two emission scenarios, the mitigated (i.e., the Representative Concentration Pathway 4.5 (RCP4.5)) and business-as-usual (RCP8.5) scenarios, based on the well-archived Coupled Model Intercomparison Project Phase 5 (CMIP5) model outputs (see 'Data' section). Under this framework, we discover a skewed monsoon cycle and diverse responses of the monsoon precipitation, hydrological extremes and compound events in six East Asian land regions. Through the lens of a complete moisture budget decomposition, the physical processes that govern the future precipitation change in various monsoon stages are pinpointed. Lastly, we discuss the sensitivity of the monsoon precipitation field to greenhouse warming and the unexpectedly 'dry-get-wetter' paradigm in East Asia.

RESULTS
A skewed EAM annual cycle Knowledge of the projected changes in the intraseasonal monsoon stages is vital to agrarian activity, freshwater security, and disaster adaptations in East Asia. Interestingly, timings and lengths of the Mid-summer, Fall, and the winter stages are highly susceptible to radiative forcings, while the temporal features of the rainy stages will barely change ( Fig. 1 and Supplementary Fig.  1). Under the high-end scenario, we expect a significantly late retreat of both the Mid-summer and the Fall over the course of the 21st century (Fig. 1b). The late retreat contributes to prolonged Mid-summer (4.4-6.3 days longer) and Fall (~2.4 days longer) in the mid-(2041-2070) and long-term (2071-2100) future (Fig. 1c). However, changes in their onset dates appear less coherent, exemplified by an earlier onset of the Mid-summer by 2.8 days but a delayed Fall onset by 3.8 days in the future (Fig. 1a). These changes collectively explain a larger growth (+15.3%) in Midsummer than Fall (+6.4%) in terms of the stage length (Fig. 1c), coinciding well with the lengthened rainy season as reported previously 10,[12][13][14] .
As surface warming can be amplified via the positive snow/icealbedo feedback 33 , it is not surprising to expect significant changes in the winter stages. Specifically, all the three winter stages will consistently delay their commencements by 3.6-6.5 days on average in the late-21st century (Fig. 1a), while the Early-and Mid-winter will retreat later by 4.3 and 5.8 days, respectively (Fig. 1b). Collectively, the Early-and Mid-winter will delay significantly without much variation in the stage days, whereas the Late-winter will shorten evidently by 2.6-4.6 days in the mid-and long-term future (Fig. 1c). Notably, the Pre-Meiyu, and the Early-and Late-winter exhibit a mild but significant reduction in their stage persistence (see 'Methods'), suggesting more intermittent Pre-Meiyu and winter monsoon seasons than before (Fig. 1d).
Taken together, the monsoon cycle in East Asia is likely to skew under a warming climate, as manifested in a markedly shortened transition from the cold to the warm season but a leisurely lengthened transition in reverse (i.e., the Fall stage). This finding locates one of the systematic changes in the future EAM cycle regardless of the high-emission or mitigated scenario (Supplementary Fig. 1).

Soaring risks of weather extremes and whiplash
Previous studies revealed a model unanimity on intensified precipitation extremes 12 and prolonged dry spells in East Asia under anthropogenic warming 13 . Here we revisit the projected likelihood of dry and heavy precipitation days (see 'Methods' for definitions) in six East Asian land regions throughout the annual cycle of monsoon. The projected changes in hydrological extremes and weather swings, as shown next, can be markedly different across regions and monsoon stages.
In general, East Asia will face soaring flood risks during the rainy season over the 21st century under both scenarios. We estimate on average 14-20 more heavy precipitation days in all land regions from the Spring to the Mid-summer in the long-term future (Fig. 2a). Remarkably, the Mid-summer alone contributes 6.4-12 more heavy precipitation days to the rainy season in all regions except the mid-lower Yangtze River basin (MLYRB), leading to a 50-200% increase in the likelihood of downpour for the stage. Such a heightened flood risk, to some extent, results from a lengthened Mid-summer under a warmed climate (see 'A skewed EAM annual cycle' section).
Floods could also catch people off guard with early arrivals under greenhouse warming. The pluvial period in the MLYRB is likely to advance, given 6.6 more heavy precipitation days in the Spring during the long-term future of the high-emission scenario ( Fig. 2a(2)). The advanced rainy season is attributable to a northward shift of future springtime rainbands since the mid-21 st century ( Fig. 3a(1)-(4)). Worse still, the enhanced spring rain in the MLYRB will, conceivably, compound with the thawing of snow and potential outbursts from ice-dammed glacial lakes in the upstream Tibetan Plateau 34 , causing more devastating and earlier floods in the basin than ever before. Likewise, concerning regions in the north (i.e., the Huai-Yellow River basin (HYRB), North China and the Korean Peninsula), Meiyu rain belts extend significantly to the north (Fig. 3c(1)-(4)). On the bright side, an earlier and lengthened rainy season for northern regions might be conducive to local agriculture. However, flood control measures will have to be exercised earlier than their historical flood season (i.e., the Midsummer) 1,17,35 to prepare for advanced flooding under a warming climate.
It may be unexpected that stronger northeasterly monsoons during the Fall and the Early-winter are tied to greenhouse warming ( Fig. 3e, f(3), (4)), giving rise to 1-5 more dry days in all regions near the end of the century (Fig. 2b). As such, East Asia will lean towards water scarcity, heatwaves and bushfires when transitioning to the cold season under a warmer climate. While suppressed precipitation over the oceans, South Japan and South China will propagate throughout the winter season, regions in the north (i.e., the HYRB, North China, and the Korean Peninsula) may expect more snowfall or cold rains given 3-6 fewer dry days in the Mid-or Late-winter ( Fig. 2b (3), (4) and (6)). This result reveals a dichotomy between a 'dry-get-drier' and a 'dry-get-wetter' winter climate in different parts of East Asia.
With the projected rise in both dry and wet extremes, we pay special attention to the weather whiplash, which denotes an abrupt shift between downpours and droughts within a short period of time (see 'Methods'). Such a compound extreme could trigger more ruinous impacts on the societies and ecosystems than individual weather extremes, yet it was rarely examined in climate change studies. In the long-term future of the highemission scenario, almost entire East Asia will experience more frequent weather whiplash than before throughout the year (Fig.  2c). For the HYRB, South Japan and the Korean Peninsula that are historically prone to weather whiplash 36,37 , they are projected to suffer more weather swings during the Spring, Pre-Meiyu, Fall and the Early-winter.
It is worth mentioning that weather extremes under the highend scenario tend to increase monotonically from near-term to long-term future ( Supplementary Fig. 2), while the mitigated scenario shares similar trends in hydrological risks at a smaller magnitude ( Supplementary Figs. 3 and 4). This signifies a big difference in the future hydrological risks contingent on our efforts to reduce carbon emissions.

Attributions of changing precipitation under global warming
Inasmuch as the prominent changes in hydrological risks during future monsoon stages, we revisit a fundamental question about the mechanisms that govern the changing precipitation in the framework of the monsoon cycle. Here we gain insights into the contributions of various physical processes from a complete moisture budget decomposition (see 'Methods'). They include the moisture-driven thermodynamic (δTH h and δTH v ), dynamic (δDY), vertical motion (δVM) and the transient eddy (δTE) effects in δ (P-E) during the long-term future under the RCP8.5 scenario, when changes are most salient. Note that the nonlinear and water  Continental precipitation that scales at a slower rate than the C-C scaling 9,11,19,20,25,26,28 is tied to the tug of war between the enhanced moisture effects and attenuated horizontal advection, vertical transport and convergence processes, especially during the warm-season stages ( Fig. 3a-d(5)). Interestingly, the moisture effects along different dimensions exert markedly diverse impacts   on the hydroclimate. For example, the horizontal moisture effect (δTH h ) explains most of the positive δ(P-E) in the MLYRB from the Spring to the Meiyu seasons, suggesting strengthened horizontal moisture transports for frontogenesis (Fig. 4a-c(1)-(3)). The vertical moisture effect (δTH v ), instead, accounts for enhanced precipitation exclusively along the southern coasts and over the Korean Peninsula, implying the role of enhanced moistening of the lower-troposphere in deep convections in the areas. As hinted earlier, the negative δDY, δVM and δTE components supported by a statistically significant model agreement, considerably offset the moisture effects in South China and the MLYRB (Fig. 4a-c(4)-(6)). These effects are commensurate with the weakened southwesterly monsoons and moisture flux convergence in the above regions ( Fig. 3a-c(4)). Notably, the negative δVM entails the tendency of descending motions that suppress the convection, which could result from the top-heavy heating effect 23,30 wherein a warmer upper troposphere leads to an increased atmospheric static stability in response to radiative forcings. In the dry season, δTH h and δTE are projected to be the main contributors to the enhanced precipitation in the north (Fig. 4f-h(2), (6)), justifying the reduced wintertime drought risks as mentioned earlier (Fig. 2b). Again, we report similar patterns of moisture budget components under the mitigated scenario ( Supplementary Fig. 5), suggesting the robustness of the competing effects in both emission scenarios. However, it should be noted that CESM-CAM and inmcm4 models were both omitted in our moisture budget analysis due to data scarcity, which unfortunately affects the moisture budget terms particularly for the Mid-summer stage where we saw a good performance by CESM-CAM. Given only a handful of grid cells with discernible residuals (Fig. 4a-h(7); see 'Methods' for discussion on uncertainties), we demonstrate the usefulness of the moisture budget analysis in attributing the causes of changing precipitation in the future EAM cycle.
A 'dry-get-wetter' paradigm in East Asia The 'rich-get-richer' pattern is likely inapplicable to East Asia land regions; instead, we find evidence of a 'dry-get-wetter' paradigm in East Asia where the historically dry regions will become wetter, which are most evident in the rainy stages. Owing to the intraseasonal migration of monsoon rainbands during the warm season 1,2,38 , monsoon onset and break periods often occur alternately across regions. During the Spring season, a positive precipitation scaling is expected in the north where the springtime climate used to be relatively dry (Fig. 3a(1), (5)). Likewise, a wetter condition appears in the northern MLYRB where it was historically dry during the Pre-Meiyu stage (Fig. 3b(1), (5)). Similar 'dry-get-wetter' paradigms are discernible in places like the HYRB and North China during the Meiyu and winter stages (Fig. 3c, f-h). Recalling the projected moisture budget components, the weakening of monsoon circulations (i.e., negative δDY and δVM) and transient eddies of moisture flux divergence (i.e., a negative δTE) explains the 'wet-not-get-wetter' paradigm in the historically wet regions (e.g., Fig. 4a-c). Conversely, horizontal and vertical moisture effects (i.e., positive δTH h and δTH v ) are the primary causes of the 'dry-get-wetter' paradigm found over land during several intraseasonal stages. Unawareness of such a 'dry-getwetter' paradigm might miss the golden timings of agricultural production and flood adaptations under climate change, potentially putting societal security at risk.
In contrast, the marine climate over the western North Pacific presents a 'rich-get-richer' paradigm 20,39,40 across monsoon stages under greenhouse warming. A 'wet-get-wetter' pattern is evident over oceans from the Pre-Meiyu to the Mid-summer (Fig. 3b-d) mainly due to the horizontal moistening effect (δTH h ) (Fig. 4b-d).
In wintertime, we notice a sudden sign switch of δTH h (i.e., a horizontal drying effect) coupled with weakened monsoon circulations (i.e., a negative δDY), giving rise to an overall 'dryget-drier' pattern over the oceans (Figs. 3f-h, 4f-h).

DISCUSSION
The present study contributes to an improved understanding of the projected changes in the intraseasonal monsoon stages that have once been a missing piece in the EAM research. In the following, we summarize several robust and systematic changes in a warming EAM climate and discuss the societal implications and the possible mechanisms at work.
In terms of monsoon seasonality, we unveil a skewed EAM annual cycle manifest in a prolonged rainy season at the cost of delayed, shortened, and weakened winter stages, which is most evident in the late-21st century under the high-end scenario. Such a lengthened rainy season originates from the extension of the Mid-summer and the Fall stages, inferring a longer transition from the warm to the cold season than the reverse way. While a stretched rainy season may favor sowing and harvesting, it also accompanies a greater and earlier flood risk in the entire East Asian continent, with alarmingly 14-20 more heavy precipitation days on average throughout the warm season in the worst-case scenario. However, future hydrological risks would be significantly reduced given substantial mitigations to reduce carbon emissions as soon as possible. Forward-looking adaptations are also required to cope with the soaring flood potential, especially during the Spring and the Mid-summer. We pay special attention to the MLYRB where floods are projected to come earlier with more severe impacts, given the enhanced springtime rainfall compounded with the upstream snow melting and the potential outburst of the ice-dammed glacial lakes 34 . Apart from a continent-wide increase of precipitation, the Mid-summer flood risk will be heightened by the antecedent wetness of the soil due to the previous rainy stages, causing entire East Asia to be more prone to floods. The projected rise of drought risks when transitioning to the prolonged Fall and the delayed Early-winter could potentially disrupt harvesting, wintertime water security, and ecosystem services. More snowfall or cold rains may occur in the HYRB, North China, and the Korean Peninsula despite a shortened winter. Notably, more frequent weather whiplash is expected in East Asia, especially over the HYRB, South Japan and the Korean Peninsula during the stages outside the peak summer and winter seasons, which could cause more severe societal impacts than weather extremes individually.
We present evidence of the rarely-documented 'dry-get-wetter' paradigm in East Asia, demonstrating the complexity of a subtropical monsoon system like the EAM, in which the commonly believed 'rich-get-richer' mechanism 39,40 could not universally apply to the region 24 . The complete moisture budget decomposition unravels the role of the horizontal and vertical moisture effects (δTH h and δTH v ) in the 'dry-get-wetter' regions during the warm-season stages, while the weakening of the dynamic effect (δDY), transient eddies of moisture convergence (δTE) and the vertical motions (δVM) elucidate the 'wet-not-get-wetter' paradigm in the historically wet regions. Broadly in tune with the previous studies 9,11,19,20,28 , we further show that the competition between the enhanced moisture effects (both horizontally and vertically) and the attenuated horizontal advection, vertical motions and convergence explains slower precipitation scaling than the C-C relation, as evident in the warm-season stages.
One possible mechanism for the skewed monsoon cycle could link with the change of the westerly jet stream. During the prolonged Mid-summer in a warmer climate, we observe a stronger westerly jet lingered on the southern Tibetan Plateau due to an enhanced meridional temperature gradient in the midtroposphere by the thermal wind relation (Supplementary Fig. 6d). The intensified jet stream in the south could facilitate the frontogenesis and convective initiations in East Asia 41,42 , and may thereby sustain the Mid-summer monsoon circulations ( Fig.  3d(1)). In contrast, faster surface warming near 40˚N, as evident in the Mid-summer and the Fall, attenuates the westerlies and induces a cyclonic circulation anomaly at the upper level ( Supplementary Fig. 6d, e). This atmospheric response may suppress the development of the wintertime Siberian high and partly explain the postponement of the winter stages (Fig. 1a, c). The cause of the weakened Asian monsoon circulations, as revealed from the moisture budget decomposition, could stem from the enhanced condensational heating in the upper troposphere in the tropics 9,13 , which consequently weakens the planetary-scale meridional thermal gradient and thereby the monsoon circulations 43,44 . The physical processes triggering the skewed monsoon cycle and the weakened monsoon circulations warrant further verification.
It should be noted that the multi-model weighting scheme proposed in this study provides a technical alternative to investigate the projected monsoon climate. It may be improved by considering further the mutual independence among the models 45 . While the systematic biases in the models may be reduced by the multi-model scheme and temporal differences in climate periods, the extent to which the biases impact the results, on both the mean state and the intraseasonal variability, is worth in-depth investigation.
Findings from this work provide policymakers, farmers, water managers and other stakeholders in East Asia, with practical information about the projected changes of the intraseasonal monsoon stages that are more pertinent to agricultural activities, climate policies and adaptations, land-use planning and water management. This study lays the foundation of our future work on the intraseasonal EAM climate under the Shared Socioeconomic Pathways (SSP) scenarios by the next generation of climate models.

METHODS Data
The fifth generation European Centre for Medium-Range Weather Forecast (ECMWF) atmospheric reanalysis (ERA5) at 1.5°gird resolution during 1979-2005 is adopted as a benchmark dataset to validate climate models' performance. We analyze outputs from 13 medium-to-high spatial resolutions climate models from CMIP5 46 (Supplementary Table 1) for their greater potentials in reproducing the EAM climate than the coarser-resolution models 47 48 , which performs the RCP8.5 projection as the CMIP5 models. The multi-ensembles mean for each model is computed for further analysis. As the spatial resolution is different across models, all the model output fields are linearly interpolated into 1.5˚gird resolution for a direct comparison to the benchmark fields.

Derivation of EAM annual cycle and the related definitions
We follow the approach developed by Dai et al. (2021) 17 to construct the EAM annual cycle, which employs a 1 × 8 Self-organizing Map (SOM) 49 to derive eight distinct monsoon stages in the EAM domain (109.5-141˚E, 19.5-45˚N) based on the daily climatology of 850-hPa horizontal winds smoothed by the 5-day-moving-mean. The clustering is trained with 850-hPa winds in the historical run  in each climate model, and the results are then evaluated against that from the training with the ERA5 reanalysis. Readers are referred to Dai et al. 17 for more details of the approach.
The stage onset date is defined as the starting day of the first five consecutive days of the corresponding stage in the year after the stage assignment 17 . Likewise, the retreat date is the last day of the last five consecutive days of the corresponding stage. Given the purely data-driven definitions, a stage's retreat date could sometimes occur later than the onset date of the ensuing stage. To avoid double-counting the days in our analysis, we investigate the number of days of a stage in one annual cycle so as to infer its duration. The 5-day persistence is defined as the conditional probability that if the first day belongs to a stage, the ensuing four days belong to the same stage. The change of the metric infers whether the stage becomes more persistent or intermittent.

Multi-model weighting scheme
Three indices considering the temporal and spatial features of the monsoon stages are designed to filter the well-performed models in reproducing the EAM annual cycle. For each monsoon stage, we compute two ratios of the time overlapped to either the observed or the model's stage period (O 1 and O 2 ) to evaluate the overall temporal matching, namely, where T o and T m is a set of calendar days of the observed and model's stage period, respectively, and ∩ denotes the intersection. The third index is the spatial pattern correlation coefficient R to evaluate the spatial similarity of the climatological 850-hPa wind fields in the corresponding observed stage and that in the model simulation. For each monsoon stage, we then select well-performed models given all three indices (i.e., O 1 , O 2 , and R) greater than 0.5. The rationale of this preprocessing is to remove poorly performed models that could undermine the confidence of the analysis. Detailed results of the model selection are listed in Supplementary Table 2.
Here, we propose a weighting scheme in terms of the model's potential in simulating the temporal and spatial features of the monsoon stage ( Supplementary Fig. 7a). Specifically, for each monsoon stage, the weighting w i assigned to a selected model i is defined as where D i denotes the Euclidean distance in climatological 850-hPa wind fields between the observed stage and the simulated stage from model i, and n denotes the number of selected models for the monsoon stage. Note that we scaled D i from 0 to 1 to make it comparable to O 2 prior to the calculation of weightings. Only O 2 is considered in Eq. (3) as it is directly proportional to O 1 (Eqs. (1) and (2)). Results using the weighting scheme largely outperform individual models based on the stages' temporal characteristics (i.e., onset, retreat and number of days), the 850-hPa winds, precipitation and temperature patterns ( Supplementary Fig. 7b-h). Albeit with a comparable performance of the multi-model weighting and the arithmetic averaging schemes, the weighting scheme is still adopted for analysis as it is based on more reliable models ( Supplementary Fig. 7b, f).

Hydrological risk analysis
For each East Asian land region, dry days (heavy precipitation days) are obtained by the number of days with the regional-total precipitation lower (higher) than the 10th (90th) percentile of the historical values  during each monsoon stage for each climate model. We compute the likelihood by dividing the dry or heavy precipitation days by the number of days of a stage. Regarding the weather whiplash (i.e., an abrupt shift between two opposite weather extremes), it is defined as a dry day and a heavy precipitation day that occur one after another within five days. The total number of weather whiplash in a monsoon stage is determined by the period between its onset date and that of the ensuing stage. The multimodel weighting scheme is applied to the likelihood of heavy precipitation days, dry days, and the frequency of occurrence of weather whiplash deduced from each model (Fig. 2).

Moisture budget decomposition
To evaluate the extent to which the changes in different processes (e.g., thermodynamic, dynamic, and transient effects) would contribute to the change of precipitation, we start from the vertically integrated atmospheric moisture budget equation 50 , which writes where P is the precipitation, E is the evaporation, w is the precipitable water in the atmosphere, ρ w is the density of water, g is the gravitational acceleration, q is the specific humidity, u is the three-dimensional wind consisting of horizontal wind velocities and pressure velocity ω, and p s denotes surface pressure. Inspired by Seager et al. 51 's derivation, the changes in the moisture budget between two time periods can be written in a complete form: where δDY ¼ À 1 ρ w g Z 1000 hPa 10 hPa ∇ h Á q his δu ð Þ dp; δVM ¼ À 1 ρ w g Z 1000 hPa 10 hPa ∂ ∂p q his δω ð Þ dp; δTE ¼ À 1 ρ w g Z 1000 hPa 10 hPa ∇ Á δðq 0 u 0 Þdp; δNL ¼ À 1 ρ w g Z 1000 hPa 10 hPa ∇ Á δuδq ð Þdp; δST ¼ À 1 ρ w g Z 1000 hPa 10 hPa δ ∂q ∂t dp; and δ Á ð Þ ¼ Á ð Þ fut À Á ð Þ his Á (13) In this context, overbars indicate the climatological stage-means and primes are the departures from the climatological stage-means. δ(·) denotes the difference in the quantity between the future and historical periods. ∇ and ∇ h indicate the three-dimensional and horizontal divergence operators, respectively. δ(P-E) refers to the change in the net flux of water substance at the surface. δTH h and δTH v represents the moisture-driven thermodynamic effect in response to the change in the horizontal and vertical gradients of specific humidity, respectively. δDY is the horizontal dynamic effect, while δVM is the vertical motion term that could be induced by both mechanical and thermal forcings. δTE represents the change in transient eddies of moisture flux convergence, δNL is the nonlinear term, and δST is the water storage term. Since the vertical integration is merely computed from 10 to 1000 hPa, here we preserve the vertical component (i.e., δTH v and δVM) to explain the precipitation change in the convective region 27,39,52 . Note that the pressure-level data underneath the topography were removed from the vertical integral. We omitted the CESM-CAM5 and inmcm4 models in the moisture budget analysis due to data scarcity. For the same token, the surface term in the original derivation 51 was incorporated into the residual term (Res). Large residuals could be found in regions with complex topography (e.g., in mountainous or coastal areas) (Fig. 4), where data with a coarse vertical resolution fail to represent the state of the system 29,50 . As daily evaporation is not available in the CMIP5 model outputs, we assume all days of the month share the same rate of evaporation bounded by the monthly amount. All the terms were first computed for each model and were then interpolated to a common 1.5˚× 1.5˚gridded resolution before performing the multi-model weighted averaging.

Statistical significance
We adopted a two-tailed Student's t-test with a significance level of 0.05 in all analyses present in this article. The t-test was performed with the null hypothesis that annual values of a given 30-year period come from a tdistribution with a mean equal to the historical mean (e.g., Figs. 1-3) with an unknown variance. For Fig. 4, the t-test was conducted on the timemean of a given period from the climate models with the null hypothesis that the mean is zero and variance is unknown.

DATA AVAILABILITY
The meteorological data is retrieved from the ERA5 by the European Center for Medium-Range Weather Forecast (ECMWF) at https://www.ecmwf.int/en/forecasts/ datasets/reanalysis-datasets/era5. The CMIP5 data can be accessed from https://esgfnode.llnl.gov/search/cmip5/. CESM-CAM5 outputs from the CESM Large Ensemble Community Project is available from https://www.cesm.ucar.edu/projects/ community-projects/LENS/data-sets.html. Derived data supporting the findings of this study are available from the corresponding author upon reasonable request.