The East Asian Summer Monsoon Response to Global Warming in a High Resolution Coupled Model: Mean and Extremes

Current climate models still have considerable biases in the simulation of the East Asian summer monsoon (EASM), which in turn reduces their reliability of monsoon projections under global warming. We hypothesize that a higher-resolution coupled climate model with atmospheric and oceanic components at horizontal resolutions of 0.25° and 0.1°, respectively, will better capture regional details and extremes of the EASM. Present-day (PD), 2 × CO2 and 4 × CO2 simulations are conducted with the Community Earth System Model (CESM1.2.2) to evaluate PD simulation performance and quantify future changes. Indeed, our PD simulation well reproduces the climatological seasonal mean and intra-seasonal northward advancement of the monsoon rainband, as well as climate extremes. Compared with the PD simulation, the perturbed CO2 experiments show an intensified EASM response to CO2-induced warming. We find that the precipitation increases of the Meiyu-Baiu-Changma band are caused by comparable contributions from the dynamical and thermodynamical components in 2 × CO2, while they are more driven by the thermodynamical component in 4 × CO2 due to stronger upper atmospheric stability. The regional changes in the probability distribution of the temperature show that extreme temperatures warm faster than the most often temperatures, increasing the skewness. Fitting extreme precipitation values with a generalized Pareto distribution model reveals that they increase significantly in 4 × CO2. Changes of temperature extremes scale with the CO2 concentrations over the monsoon domain but not for precipitation extreme changes. The 99th percentile of precipitation over the monsoon region increases at a super Clausius-Clapeyron rate, ~ 8% K–1, which is mainly caused by increased moisture transport through anomalous southerly winds.


Introduction
The East Asian summer monsoon (EASM) is one of the most vibrant climate system components, encompassing a large spatial coverage and a wide range of time scales. Its seasonal evolution is characterized by the northwestward movement of a zonally elongated Meiyu-Baiu-Changma rainband strongly modulated by the western Pacific subtropical high (WPSH) and the dynamical effect of the Tibetan Plateau, providing water resources for one third of the world's population (e.g., Lau and Ploshay 2009;Son et al. 2019Son et al. , 2020Oh and Ha 2020;Son and Seo 2022). Furthermore, climate extremes (e.g., heatwaves and floods) are closely linked to variations in the EASM, which pose severe threats to human lives and property Zhou and Zou 2014;Wilcox et al. 2015;Lee et al. 2017). Online ISSN 1976-7951 Print ISSN 1976 Korean Meteorological Society Numerical models have proven to be powerful tools in reproducing the basic structure of the monsoon circulation and precipitation (Lau and Ploshay 2009;Xin et al. 2020;Sperber et al. 2013). However, some model biases, such as the monsoon precipitation distribution (Huang et al. 2013), seasonal evolution (Zou and Zhou 2015), and associated climate extremes (Chen et al. 2021) on different time scales, have persisted for decades. These biases possibly result from the unrealistic representation of physical processes (e.g., convection, Liu et al., 2021) and of complex terrain effects (Boos and Hurley 2013;Lee et al. 2015). Previous studies suggested that increasing the model resolution can considerably improve the monsoon performance in both atmosphereonly (e.g., Kitoh and Kusunoki, 2008) and coupled models (e.g., Yao et al. 2017;Chen et al. 2021).
Most global monsoon systems are projected to have a weakened circulation and enhanced precipitation in the Northern Hemisphere under global warming . Besides the humidity increase, greenhouse gas warming generally affects the monsoon dynamically in two contrasting ways (Seo et al. 2013;Lee and Wang 2014;Wang et al. 2020; Moon and Ha, 2020). On one hand, due to relatively larger temperature increases in the upper atmosphere, increased atmospheric stability suppresses convection and precipitation. On the other hand, increasing greenhouse gases induce a land-sea contrast and can act oppositely on different time scales (Shaw and Voigt 2015;Chen and Bordoni 2016). The fast responses due to a reduction of radiative forcing in the absence of sea-surface temperature (SST) changes cause land to warm faster than the ocean, intensifying the WPSH and enhancing the inland moisture transport and subsequent convergence and precipitation increases (Li et al. 2019). On longer time scales, the subsequent SST warming weakens the land-sea thermal contrast, resulting in decreases in monsoon rainfall and circulation (Chen and Bordoni 2016). Despite its important role, the dynamic response of the EASM to CO 2 -induced warming is still largely uncertain, poorly understood and warrants further investigation Zhou et al. 2018). Chen et al. (2018) reported an underestimation of future Meiyu rainfall changes under global warming in coarse resolution models due to a too weak southwesterly response associated with an eastward shift of WPSH. They claimed that increasing model resolution to resolve complex terrain and parameterizations is essential to reduce uncertainties in the projection of EASM changes. Moreover, higher resolution models are also found superior to relatively lower resolution models in simulating precipitation extremes (Iles et al. 2020;Chen et al. 2021) as well as their projections in responses to CO 2 forcing (van der Wiel et al., 2016). This conclusion is also applied to models including air-sea interactions compared to atmosphere-only simulations (Fischer et al. 2018;Kim et al. 2021).
In this study, we used data from coupled experiments with one of the highest atmospheric and oceanic resolutions performed so far, to evaluate its performance on EASM and their responses to different CO 2 concentration levels across various time scales. We show that these high resolution coupled simulations produce better performances and projections of the monsoon, as well as advance our understanding of physical mechanism in CO 2 -forced monsoon changes. The model experimental design and analysis method are documented in Section 2. Results are described in Section 3. Finally, main conclusions are given in Section 4.

Model and Observation-Based Data
We use three century-long simulations of the Community Earth System Model (CESM) version 1.2.2. The horizontal resolutions of atmospheric and oceanic components are around 0.25° and 0.1°, respectively, ranking among the finest resolved coupled model simulations conducted up to now. The present-day simulation (PD) is run for 140 years with a CO 2 concentration level of 367 ppm corresponding to the Year 2000 condition in as Small et al. (2014). This level enables us to make comparisons with recent-era climatology but note that differences may potentially come from the constant forcing used when comparing with transient reanalysis. From year 71 of the PD simulation, two new simulations were spun up, where the CO 2 concentrations were increased to 2 times (734 ppm; 2 × CO 2 ) and 4 times (1468 ppm; 4 × CO 2 ), and run for 100 years for each experiment. The details of the model set up have been documented in Chu et al. (2020) and Huang et al. (2021). Only results based on the last 40-years, after the simulations reached near-equilibrium as shown in the timeseries of global surface temperature, radiative forcing at the top of atmosphere (see Fig. S1 in Chu et al., 2020), and intensity of the WPSH and EASM (Fig. S1), are analyzed. In contrast to lower resolution simulations, the PD experiment shows considerable improvements in the annual-mean spatial distribution of precipitation and SST (Chu et al. 2020), which may substantially improve the projected SST warming and land-sea contrast. Besides, the high resolution coupled simulation is also expected to resolve the local terrain effect (e.g., the Maritime Continent; Chang et al., 2016;Yamanaka et al., 2018) well and, thus, to improve the monsoon representation (e.g., Liu et al., 2021). This improved air-sea interaction is critical to simulating the land-sea thermal contrast and subsequent monsoon circulation (Wang et al. 2005). With these experiments, we evaluate how these high resolution coupled simulations replicate the EASM and investigate its response to CO 2 forcing by examining the differences between perturbed CO 2 and the PD simulations.
The seasonal mean in summer (June-August) and intraseasonal features on pentad scales are evaluated against observations at a horizontal resolution of around 25 km. For evaluating the simulation quality of precipitation, we choose the Asian Precipitation-Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APH-RODITE; Yatagai et al. 2012) dataset during 1976to 2015Huffman et al., 2007) during 1998 to 2019, and the ECMWF fifth generation reanalysis (ERA5; Hersbach et al., 2020) during 1981 to 2020 for daily precipitation. The simulated temperature field is assessed against APHRODITE and ERA5 reanalysis, and the latter is also used to evaluate the simulated large-scale circulation.

Analysis Methods
To disentangle the contribution of circulation and humidity changes accounting for precipitation changes, the moisture budget is analyzed by decomposing the moisture flux convergence (MFC) as in Seager et al. (2010) andD'Agostino et al. (2019). The changes of precipitation minus evaporation (P -E) should approximately be equal to changes of the vertically integrated MFC: Here ρ w is the water density, g is the gravitational constant, p s is the surface pressure, p t is the top pressure (i.e., 10 hPa), u is the horizontal wind, q is the specific humidity, and S is the surface quantity related to moisture transport by surface flow, which is usually small and not further discussed here. The overbar denotes the monthly mean, and the prime represents daily anomalies relative to the monthly mean. Thus, the first and second terms on the right-hand side represent the monthly moisture transport changes by the dynamical (DY) and thermodynamical (TH) components, respectively. The third term denotes the contribution from transient eddies on time scales less than one month. We checked that the quadratic term (e.g., ∫ p s p t ∇ • u q dp ) is small and, similar to Oh et al. (2018), we neglected it.
The monsoon onset and withdrawal dates are defined as the first and last pentad of the rainy season, respectively, whose rainfall amount exceeds and drops below the January mean precipitation by 5 mm/day over the monsoon domain (Wang and Ho 2002). The monsoon domain refers to the region where the differences between maximum pentad precipitation and January mean precipitation are larger than 5 mm/day during the extended summer (May-September).
Two extreme climate indices, the 95 th and 99 th percentile of daily near-surface air temperature (T95 and T99) and precipitation (R95 and R99) are analyzed in this study. The generalized Pareto distribution (GPD) is used to model the exceedances over a precipitation threshold (Coles 2001). Comparing the shape and scale parameters of the distribution in different experiments helps us to identify whether the changes of extreme precipitation are significant or not. The significance of the spatial distribution of the differences between increased CO 2 runs and the PD simulation is estimated by the Student's t test. The non-parametric Kolmogorov test is adopted to assess the significance of the probability density distribution differences between the CO 2 perturbed and PD simulations. Figure 1 shows the distribution of the climatological summer mean precipitation and 850-hPa wind in observations and the PD simulation. The precipitation pattern is characterized by its maxima on the windward side of the high terrain continent over the western edge of India, Indochina and the nearby Bay of Bengal, and the Philippines, along the Himalayas, and over the western tropical Pacific, South Korea, and southern Japan. The PD simulation overestimates the precipitation over the Himalaya as does the ERA5 reanalysis. This is consistent with previous CESM simulations (e.g., Palazzi et al. 2015). Over the tropical regions, the simulated westerly winds are much weaker over the Bay of Bengal and the South China Sea compared with those in the ERA5 reanalysis. As a result, the precipitation deficit appears over the west coast of Indochina and the Philippines and surrounding oceans to their west, while excessive precipitation occurs in the western Bay of Bengal and the South China Sea, exhibiting a zonal dipolar pattern of biases. Chen et al. (2019) reported a similar finding in a multiyear 3-day hindcast simulation with a well-constrained synoptic scale circulation, and thus attributed the biases to local circulation errors associated with fast physical processes (e.g., convection and the diurnal cycle). Another possible reason could be that a resolution of 25 km still cannot resolve local terrain peaks, leading to a weak blocking effect and less precipitation on the windward side (e.g. the Arakan Mountains in Indochina; Wu et al. 2014). These dry biases near the western coast of Indochina and the Philippines generate weaker low-level convergence and westerly winds as a positive feedback. Consequently, the tropical westerlies easily turn northward toward South China, showing a stronger EASM circulation and precipitation in the PD than the observations. Meanwhile, the WPSH is stronger in PD than ERA5 as manifested by the excessively strong tropical easterly winds reaching at around 125°E compared to 130°E in the latter, resulting in dry biases over the western tropical Pacific and a stronger monsoon circulation. We note that most of the coupled models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive tend to underestimate the EASM precipitation associated with a weaker WPSH and monsoon circulation than observations Fig. S7 in Wilcox et al., 2020).Whereas discrepancies are present in the PD, the model well replicates the distribution of the observed seasonal mean precipitation and 850-hPa wind, with spatial correlation coefficients of precipitation above 0.70 (statistically significant at the 1% level). Therefore, our high resolution coupled simulation can reproduce well the seasonal mean features of the EASM circulation and precipitation.

Evaluation of Present-Day Simulation
The intra-seasonal evolution of the EASM is accompanied by a northward advancement of an elongated rain belt following the northward movement of the WPSH from early June to mid-July (Ding 1992;Wang and Ho 2002;Chu et al. 2017). During the 31 May-9 June period, the heavy rainband extends from South China to the southern tip of Japan ( Fig. 2a and e). This rain belt intensifies and reaches the Yangtze River basin and southern Japan in the following two pentads associated with enhanced WPSH and southerly winds coming from the eastern Bay of Bengal and the South China Sea ( Fig. 2b and f). By the end of June, the intense precipitation zone is situated at central China stretching northeastward to central Japan together with the northward shift of the WPSH ( Fig. 2c and g). Then the WPSH moves further north but with slightly weakened southwesterly winds on its northwestern flank during early July ( Fig. 2d and h). The rainfall belt covers South Korea and most of Japan and becomes weaker. The model well reproduces the slanted and elongated rainfall belt and its stepwise northward advancement associated with the WPSH northward shift, although the monsoon arrival date is about 2 pentads earlier than observed ( Fig. 2i-l). This finding is similar to that reported in Lau and Ploshay (2009) using an atmospheric general circulation model at a resolution of around 50 km. However, our PD simulation produces a much closer match in magnitude and pattern of rainfall to observations than their results. Inspection of the low-level circulation shows that the earlier onset of the rainy season in the simulation is likely induced by the too strong southerly wind over East Asia. The excessively powerful monsoon circulation originates mainly from the excessively strong southerly wind coming from the Bay of Bengal, the South China Sea, and the western tropical Pacific associated with dry biases there. The precipitation and low-level circulation biases on the intra-seasonal scale resemble those on the seasonal mean scale suggesting a common cause for these biases as discussed before.
Extreme weather and climate phenomena are closely linked with variations of the large-scale circulation patterns over Asia during the monsoon season (e.g., Chen and Zhai 2014;Luo and Lau 2017). The patterns of T95 and T99 are very similar, with higher temperature over the land than the ocean (Fig. 3a, b, d and e), particularly over northern India and central eastern China where a strong monsoon circulation occurs. The PD experiment generally overestimates the temperature extremes, especially over warm areas ( Fig. 3c and f). Despite this overestimation, the model well captures the pattern of temperature extremes over Asia with a pattern correlation of about 0.95 (statistically significant Year 140. The numbers in the parenthesis denote the pattern correlation of precipitation between observations and PD at the 1% level). The pattern of precipitation extremes overall follows the seasonal mean precipitation distribution (Fig. 4). There is a wide spread in the observed precipitation extremes, with TRMM as the highest and APHRO-DITE as the lowest in terms of intensity, consistent with Kim et al. (2019). Their differences can be more than 50 mm/ day for the R99 values over the windward regions, such as the Bay of Bengal, the South China Sea, and the western tropical Pacific ( Fig. 4e and f). Previous studies found that TRMM generally overestimates heavy precipitation (Zhao and Yatagai 2014) while APHRODITE underestimates the precipitation extremes (Lai et al. 2020). The simulated magnitude of precipitation extremes falls into the range of the observations (Fig. 4d and h). The R95 values in the PD simulation are near the lower bound of the observations (i.e., APHRODITE) while R99 is closer to the upper bound (i.e., TRMM). The pattern correlations between simulations and various observations are above 0.55 (significant at the 1% level), suggesting the model can reproduce well the synoptic features associated with the monsoon circulation. Overall, the high resolution PD simulation has a good fidelity of the monsoon precipitation and circulation on seasonal, intraseasonal, and synoptic scales, which enables us to examine their responses to CO 2 forcing with more confidence. Figure 5 shows the responses of precipitation and the lowlevel circulation to doubling and quadrupling of the CO 2 concentration. In the following analysis, we also include the Indian monsoon region which shows a weaker monsoon circulation in responses to CO 2 forcing, consistent with a generally projected weakening of the global monsoon circulation . This weakening of the South Asian monsoon circulation due to CO 2 forcing mainly serves as a comparison to the enhancement of the EASM. By doubling the CO 2 concentration level, the circulation changes over the northern subtropics are characterized by anomalous anticyclonic circulation centered over southern India and the western Pacific and a trough centered at Indochina. This pattern weakens the climatological trough over eastern India while it reinforces the climatological southerly wind over East Asia, suggesting a weakened South Asian monsoon circulation and an enhanced EASM circulation. Under the influence of southwesterly oceanic flow, southwestern India and southern China experience a precipitation increase by about 20%. Moreover, the precipitation increases are on the northern side of the climatological rainfall band over eastern China, indicating a northward shift of the rainband Fig. 2 The spatial distribution of climatological precipitation (mm/day) in TRMM for the (a) period 1 (P1; 31 May-9 Jun), (b) P2 (10 Jun-19 Jun), (c) P3 (20 Jun-29 Jun), and (d) P4(30 Jun-9 Jul). (e-h), (i-l) Same as (a-d) but for precipitation and 850-hPa wind (m/s) based on ERA5 and PD, respectively due to enhanced southerly anomalies ( Fig. 1 and 5a). These anomalous southwesterly flows over western India and East Asia are mainly attributed to the CO 2 -induced differential horizontal warming patterns manifested by a warmer Northern than Southern Hemisphere and warmer land than ocean (Lee and Wang 2014;Lau and Kim 2017;Wang et al. 2020). Consequently, there are anomalously strong cross-equatorial easterly flows over the southern Indian Ocean because of the interhemispheric temperature contrast, which then turn to southwesterlies towards the Indian subcontinent. Over East Asia, warmer land causes low-pressure anomalies to extend from eastern China to the northwestern Pacific across Korea and Japan (Fig. 6d). By contrast, relatively colder oceaninduced high-pressure anomalies appear over the western subtropical Pacific, showing a stronger WPSH in response to CO 2 forcing (Lee and Wang 2014). The sea-level pressure gradient between land and ocean leads to intensified southwesterly wind prevailing over the East China Sea and the northwestern Pacific (Fig. 5a). To the north, there is an anomalous cyclonic circulation centered over Japan in conjunction with the southward movement of the jet over the Pacific sector (Fig. 6a). Note the meridional shift of the jet stream due to greenhouse warming is of considerable uncertainty in climate models, mainly depending on the stratospheric cooling and Arctic surface warming (Shaw et al. 2016;Yim et al. 2016). In response to quadrupling the CO 2 concentration, the overall patterns of changes in precipitation and low-level circulation are very similar to that induced by doubling CO 2 , except that the magnitude is much larger, and more significant changes are shown (Fig. 5b). Larger interhemispheric temperature gradients and land-sea pressure gradients (Fig. 6h) lead to much stronger southerly wind anomalies over East Asia reaching northeastern China in comparison to south of 30°N in the 2 × CO 2 experiment. More moisture is transported to East China resulting in precipitation increases, showing an even stronger monsoon circulation in the 4 × CO 2 than the 2 × CO 2 simulation. Compared to 2 × CO 2 , the jet stream is much weaker (Fig. 6a and  e), which is also conducive to stronger low-level southerly wind anomalies. Note the increased precipitation is situated at the location of the climatological rainfall band over eastern China ( Fig. 1 and 5b), suggesting an intensification in the 4 × CO 2 in contrast to a northward movement of Meiyu in the 2 × CO 2 , albeit with stronger southerly winds in the former. This may indicate different mechanisms in driving the rainfall enhancement, which will be discussed below.

EASM Responses to Global Warming
The enhanced moisture transport to inland areas and subsequent convergence can be induced by the increases in either convergence or specific humidity. To further understand their contributions to the precipitation increase over Asia, we analyze the moisture budget following Seager et al. (2010). The changes in MFC are further decomposed into The numbers in the parenthesis denote the pattern correlation between observations and PD the monthly DY via circulation changes, the monthly TH via humidity changes, and the contribution from transient eddies on time scales less than a month as in previous studies (D'Agostino et al. 2019;Seager et al. 2010). The distribution of the P -E changes (Fig. 7a, e) is very similar to that of the precipitation changes (Fig. 5a, b), suggesting a minor role of evaporation changes. In response to doubling CO 2 , there is a decrease of P -E over central India while at the same time there is an increase over southwestern India (Fig. 7a). This decrease over the core monsoon region is consistent with the low-level anticyclone (Fig. 5a) and downward motion anomalies (Fig. 6b). To the east, increased precipitation appears over Indochina and extends from South China to Japan across South Korea (Fig. 7a), corresponding to the anomalous trough over Indochina and intensified WPSH, respectively (Fig. 5a). Inspection of the decomposition shows that the pattern of P-E changes basically follows that of DY ( Fig. 7a and b), indicating the vital role of circulation changes in shaping the pattern of monsoon precipitation changes Lau and Kim 2017). The general negative anomalies in DY (Fig. 7b) mean suppressed convection associated with low-level divergence. This enhanced divergence is mainly due to the increased atmospheric stability with more greenhouse gases (Fig. 6c), which has been discussed in previous studies Lee and Wang 2014). An exception is southern China and mid-latitude inland regions where positive DY anomalies are presented corresponding to the intensified monsoon circulation due to the strengthened sea-level pressure gradient (Fig. 6d). Note the positive anomalies are in the north of the climatological rainband over eastern China, suggesting the driving role of DY in the northward shift of the rainband. Besides, the decreased low-level stability over land is also conducive to ascending anomalies and subsequent low-level convergence (Fig. 6c). The decreased stability due to increased CO 2 forcing is likely induced by the stronger low-level temperature increases over middle-to high-latitude inland areas in boreal summer (see Fig. 10a in Wang et al. 2020). TH generally contributes a relatively uniform pattern of positive anomaly (Fig. 7c) that mainly results from increased specific humidity as temperature increases (Held and Soden 2006). Over the core South Asian monsoon region, the wet anomalies due to TH only partially offset the contribution from DY, leading to dry anomalies there. In the main rainband of the EASM, those two components contribute to comparable wet anomalies. However, previous studies show more important DY effects than TH effects on the EASM precipitation (Oh et al. 2018;Li et al. 2019). Li et al. (2019) found that DY plays a dominant role for precipitation increase for a warming of 2-5 K above the pre-industrial global mean surface temperature. As warming level increases, DY can be two times of TH in generating positive precipitation anomalies. One possible reason accounting for the smaller contribution of DY in this study is the increased stability in the upper atmosphere in response to CO 2 forcing (Fig. S2), which can be measured by the temperature differences between 500 and 200 hPa. The increased upper atmospheric stability inhibits the ascending motion, deep convection, and subsequent low-level convergence, thus reducing the DY contribution. Besides, the wide inter-model spread of responses to warming in DY and TH is found in the CMIP5 (IPCC, 2013) and CMIP6 models , which may also explain the relatively stronger TH contribution in this work. Interestingly, the pattern of contribution from transient eddies is opposite to that from DY, although its magnitude is much smaller (Fig. 7d). In the 4 × CO 2 simulation, the pattern of sea-level pressure changes is very similar to the 2 × CO 2 simulation but intensified by a factor of around 2 ( Fig. 6d and  h). Consequently, the circulation responses are enhanced, characterized by a stronger anticyclone over India and intensified southwesterly over East Asia. (Fig. 5b). The changes in DY are also stronger for 4 × CO 2 than for 2 × CO 2 , but weaker in the rainband over East Asia (Fig. 7b and f). This is because the further enhancement of upper atmospheric stability (Fig. S2) suppresses the ascending anomalies and deep convections (Fig. 6f) and offsets the enhanced monsoon responses in 2 × CO 2 although the lower atmosphere is more unstable over East Asia (Fig. 6g). By contrast, the enhanced upper atmospheric stability in 4 × CO 2 reinforces the weakened monsoon responses in 2 × CO 2 over the South Asian monsoon region. The weakened increase of the rainband due to DY over East Asia in 4 × CO 2 than 2 × CO 2 indicates that there is an optimal CO 2 -induced warming level, which maximizes the influence of DY on precipitation increases over the monsoon front. Not surprisingly, the TH contribution further increases through humidity increases as a response to stronger warming for 4 × CO 2 than 2 × CO 2 ( Fig. 7c and g). The precipitation increases due to TH outweigh the decreases caused by DY ( Fig. 7f and g), leading to an overall increased precipitation over India for 4 × CO 2 (Fig. 7e). Over China, TH accounts for more precipitation increases than DY ( Fig. 7f and g), suggesting a dominant role of the former for 4 × CO 2 . The Meiyu band is intensified over eastern China mainly following TH changes in 4 × CO 2 instead of shifting northward in 2 × CO 2 because DY is suppressed due to further enhanced upper-atmospheric stability in 4 × CO 2 . We further decompose the DY term into advection and convergence components following Han et al. (2019). The advection changes between 4 × CO 2 and 2 × CO 2 show increases over the Meiyu band (Fig. S3b) in conjunction with the intensified southwesterly winds (Fig. 5). The decreased DY effects over eastern China is mainly driven by convergence changes (Fig. S3a) modulated by the strengthened upper atmospheric stability, indicating a coupled TH and DY impact through temperature changes.
The monsoon onset and withdrawal represent the start and end dates of the rainy season, respectively. Their predictions are crucial for agriculture and the economy over Asia but remain largely uncertain under global warming. Some studies suggest an earlier onset and later withdrawal of the EASM Lee and Wang 2014) while other The patterns of onset and retreat dates in PD are similar to those in previous studies using observations (Wang and Ho 2002;Dallmeyer et al., 2015), showing the model's capability in capturing the key monsoon evolution (Fig. 8a and  d). In response to a doubling of the CO 2 concentration, the onset changes are characterized by a dipole pattern, with positive anomalies over South Asia (e.g., Indochina) and negative anomalies over East Asia and the western Pacific (Fig. 8b). This means a later onset in the former and an earlier onset in the latter. Increases of the withdrawal date are shown over the monsoon domain, indicating a late withdrawal (Fig. 8e). The larger increases of the withdrawal than onset date over India and surrounding oceanic areas suggest a longer rainy season there . By contrast, the stronger positive anomalies of onset than withdrawal over Indochina mean a shorter rainy season there. Over East Asia, the lengthened rainy season is attributed to both the earlier onset and delayed withdrawal (Cui et al. 2020;Ha et al. 2020;Moon and Ha, 2020). The generally stronger responses of withdrawal than onset except over Indochina suggest its higher sensitivity to CO 2 -induced warming so that the changes of the rainy season length are mainly driven by withdrawal changes. With quadrupling CO 2 forcing, the magnitudes of changes in both onset and withdrawal dates are more enhanced than those with doubling CO 2 , particularly for withdrawal ( Fig. 8c and f). Over the Yangtze River basin, the earlier onset (~ 4-6 pentads) and delayed withdrawal (~ 8-10 pentads) can prolong the monsoon season by more than 2 months. Note that only withdrawal changes over eastern China and the Arabian Sea in 4 × CO 2 are shown as evidently significant with respect to their interannual Fig. 6 The spatial distribution of summer changes in (a) 300-hPa zonal wind (m/s), (b) 500-hPa vertical velocity (0.01 Pa/s), (c) temperature differences between 850 and 500 hPa (K), and (d) sea-level pressure (hPa) between 2 × CO 2 and PD simulation. (e-h) Same as (a-d) but for 4 × CO 2 simulation. Black dots represent the changes at 0.05 significance level. Blue contours in (a, e) denote the zonal wind climatology in PD.

Red (blue) in (b, f) indicates descent (ascent)
Korean Meteorological Society variability in PD. The interannual variations of onset and withdrawal dates are large, particularly in the latter (Fig. S4), which is likely because of the large variability of pentad precipitation (e.g., Montes et al., 2021). Interestingly, the duration of the monsoon season over Indochina becomes longer in 4 × CO 2 , in contrast to the shorter one in 2 × CO 2 because of a more delayed retreat for 4 × CO 2 concentration.
Climate extremes are projected to increase significantly under global warming over the Asian monsoon regions (Zhang and Zhou 2019;Kim and Bae 2020;Lee et al., 2021). Given its large population density, Asia is highly vulnerable to climate extremes (e.g., heatwaves and floods). However, regional projections of climate extremes are still rather uncertain in the latest climate models (Seneviratne and Hauser 2020), which may come from insufficient spatial resolution, particularly for precipitation extremes (van der Wiel et al., 2016;Iles et al. 2020). Table 1 shows precipitation extremes pooled over six representative inland monsoon regions in response to CO 2 forcing.
The ranges of the R95 increases are from 1 to 32% for 2 × CO 2 and from 16 to 40% for 4 × CO 2 , which are much wider than that for temperature ( Table 1). The increases of R99 are generally larger than those of R95, suggesting that stronger precipitation extremes tend to be even stronger as CO 2 increases. This means a higher sensitivity of stronger precipitation extremes to CO 2 forcing. On average, the intensification of R95 over surface temperature changes is approximately 5%-6% K -1 while that of R99 is roughly 8% K -1 , with the Clausius-Clapeyron (CC) relation about 6.5% K -1 falling between both (Allen and Ingram 2002). This means that R95 is still not extreme enough to effectively consume all the atmospheric moisture to precipitate out as only the uppermost percentage of extreme changes follows the CC relation assuming marginal impact from Fig. 7 The spatial distribution of summer changes (mm/day) in (a) precipitation minus evaporation (P -E), (b) dynamic term (DY), (c) thermo-dynamic term (TH), (d) transient term between 2 × CO 2 simulation and PD. (e-h) Same as (a-d) but for 4 × CO 2 simulation circulation changes (Pall et al. 2007). By contrast, the super CC scaling of R99 is likely because the anomalous southerly winds bring more moisture via dynamical transport. Unlike the temperature extremes, the regional mean precipitation extreme increases do not scale with CO 2 concentration, which are on average below and above a factor of 2 for R95 and R99, respectively.
To determine whether precipitation extremes will increase significantly under CO 2 -induced warming, the GPD is applied to model the tail of the precipitation distribution exceeding their 99.99 th values (R99.99) at the present day for each region. There are three parameters to determine the GPD: location, shape, and scale (Coles 2001). Larger values of scale and shape mean stronger extreme values, and the location parameter refers to the threshold above which the tail fits the GPD well. We choose R99.99 as the threshold because we cannot reject the GPD as the appropriate distribution of the precipitation data at the 1% significance level for all regions, particularly for larger regions with more samples. In this sense, Fig. 9 shows how intense precipitation extremes respond to CO 2 forcing. Error bars denote the standard deviations of the shape and scale parameters. Fig. 8 The spatial distribution of monsoon onset date (pentad) (a) in PD, (b) changes between 2 × CO 2 simulation and PD, and (c) changes between 4 × CO 2 simulation and PD. (d-f) Same as (a-c) but for monsoon withdrawal date. Black dots denote the region where the responses are larger than one standard deviation of onset/withdrawal dates in the PD experiment Table 1 The average of 95 th and 99 th percentile of nearsurface air temperature (K) and precipitation (%) changes, and surface temperature (K) changes between 2 × CO 2 , 4 × CO 2 and PD simulations over different land regions. The shape parameters are larger in CO 2 perturbed simulations than PD over South China and Japan, suggesting significant increases in precipitation extremes. Similarly, the larger scale parameters over North China indicate stronger extremes there. Over South Korea and Central East China, the stronger precipitation extremes are only seen in 4 × CO 2 but not for 2 × CO 2 . Over India, the scale parameter is larger in 4 × CO 2 compared to PD, but the shape parameter is relatively smaller. Thus, there is apprarent uncertainty whether R99.99 over India will significantly increase or not in response to 4 × CO 2 forcing. Besides precipitation, all the regions experience temperature increases of about 3.5 and 7.0 K on average in 2 × CO 2 and 4 × CO 2 , respectively (Table 1). This means that the temperature increases over the monsoon domain roughly scale with the CO 2 concentration. The increase of T99 is slightly larger than that of T95, showing that stronger extremes warm more in response to CO 2 forcing (Duan et al. 2020). This conclusion is also confirmed by the probability density function of near-surface air temperature for each experiment (Fig. 10). The probability density differences among different experiments are significant at the 1% level according to the Kolmogorov-Smirnov test. In Fig. 10, the vertical dashed lines generally shift to the right faster than the vertical solid lines, indicating that T95 increases faster than the mode of the distribution. In other words, warmer temperatures warm more under CO 2 -induced warming, which will increase the skewness of the temperature distribution. Besides, the mode of the distribution decreases as CO 2 concentration increases, except over Japan, suggesting a wider temperature range with more CO 2 . Note there is skewness in the distribution over South Korea and Japan, which is not present in the observations, implying model biases in reproducing temperature distribution over islands. This bias may be because the model cannot well simulate the anomalous advection of the stationary temperature gradient and the covariance between wind and temperature anomalies, which are critical for temperature skewness . The underlying causes are unclear and beyond the scope of this study.

Conclusions and Discussion
In this study, we have assessed the performance of the EASM circulation and precipitation in one of the most finely resolved coupled global model simulations conducted so far, with horizontal resolutions of around 0.25° and 0.1° for the atmospheric and oceanic components, respectively. The PD simulation well replicates the spatial distribution and intensity of monsoon features on seasonal mean, intraseasonal, and synoptic scales albeit with overestimated southerly winds over East Asia. This stronger monsoon is accompanied by weaker subtropical westerly winds associated with dry biases over the western coast of Indochina and the Philippines, and stronger subtropical easterly winds associated with an excessive powerful WPSH. These biases are possibly linked to local circulation errors associated with fast physical processes (e.g., convections and diurnal cycle) and unresolved terrain (Chen et al. 2019;Wu et al. 2014). Spatial resolutions higher than 25 km in atmospheric models may be needed to further improve the monsoon performance. Cloud resolving simulations may be essential to capture the coastal diurnal cycle and land-sea breeze associated with local complex terrain and convections over the Maritime Continent and the Indochina Peninsula (Oouchi et al. 2009;Li et al. 2020;Liu et al., 2021). A caveat is that increasing the horizontal resolution does not necessarily yield improved monsoon circulation and rainfall relative to observations (Lim et al. 2014). The model performance may be degraded as the convection schemes in many global coupled models are originally designed for relatively coarse resolutions (Wang et al. 2022). Thus, scale-aware parameterizations may be necessary and should be adjusted to match up with high resolution simulations (Johnson et al., 2016). The optimized horizontal resolution including relative importance of atmospheric and oceanic resolution (e.g., Sein et al., 2018;de la Vara et al., 2020;Jullien et al., 2020) for improving monsoon performance also warrants further investigation, which is beyond the scope of this study.
Additional CO 2 perturbed experiments, 2 × CO 2 and 4 × CO 2, have been analyzed in order to investigate the responses of the EASM to CO 2 -induced warming by comparing with PD. In responses to doubling the CO 2 concentration, the EASM circulation is characterized by stronger southwesterly anomalies associated with the intensified WPSH (e.g., Zhou and Zou 2010;Dai et al. 2013;Endo et al., 2018) and an anomalous trough over Indochina. This finding is in contrast to a weakened circulation over most monsoon regions under global warming (IPCC, 2013;Li et al. 2019). By contrast, the South Asian monsoon is weakened, characterized by an anomalous anticyclone over the Indian subcontinent opposing to the climatological trough over eastern India. The anomalous southwesterly over northern Arabian Sea and southern China is mainly driven by the sea-level pressure gradient via the CO 2 -induced warmer Northern Hemisphere than Southern Hemisphere Only land grid points are considered. The vertical solid and dashed lines denote the mode and the 95 th percentile of the distribution, respectively and warmer land than ocean, respectively. These anomalous flows bring abundant moisture to southwestern India and southern China, leading to precipitation increases there by 20%. Compared to doubling the CO 2 concentration, the patterns of changes in the monsoon circulation and precipitation are very similar in 4 × CO 2 but with a much stronger magnitude.
To disentangle the relative roles of circulation and humidity changes in driving monsoon precipitation anomalies due to CO 2 forcing, the MFC changes are decomposed into contributions from monthly DY and TH components and transient eddies on time scales less than a month. In 2 × CO 2 , the changes of DY resemble those of P -E, suggesting its dominant role in shaping the precipitation responses Lau and Kim 2017). The DY and TH components contribute to comparable wet anomalies of the rainband over East Asia, with marginal impact from transient eddies. In 4 × CO 2 , TH outweighs DY, leading to general precipitation increases over the monsoon region. Compared to 2 × CO 2 , the decreases in wet anomalies of the rainfall band due to DY is caused by the increased upper-atmospheric stability. This means that the CO 2 -induced temperature changes can in turn modulate the circulation changes and subsequent DY effect, suggesting a coupling between the DY and TH components. Our decomposition of MFC cannot totally separate the contribution from circulation changes and temperature changes although this method initially is used to quantify the relative importance of circulation changes and humidity changes. One possible solution could be to conduct additional CO 2 experiments with the circulation constrained toward that in the PD experiment to disentangle the DY effects using a dynamical nudging method (e.g., Liu et al., 2021).
The monsoon season becomes longer over East Asia in response to doubling the CO 2 concentration due to the advanced onset and delayed withdrawal, particularly in the latter, consistent with previous studies (Lee and Wang 2014;Cui et al. 2020). This increased duration may be related to seasonal-mean precipitation increases under global warming (Goswami and Xavier 2005). The late withdrawal may be associated with higher mean sea-level pressure anomalies in the tropical western Pacific in favour of moisture transport to monsoon regions (Kitoh and Uchiyama 2006). The early onset is suggested to be linked to the convective activity over the tropical western Pacific (Jin and Stan 2019). In 4 × CO 2 , the changes in onset and withdrawal dates are amplified, particularly in the latter, suggesting its higher sensitivity to CO 2 forcing. The monsoon duration can be longer by 2 months than that in PD over the Yangtze River basin.
For precipitation extremes, R95 increases at a rate of 5%-6% K -1 on average over the monsoon region while R99 increases at a larger rate of ~ 8% K -1 . The super CC scaling for R99 is mainly due to more oceanic moisture transport through anomalous southerly winds. On average, the precipitation extreme increases over the monsoon domain are not scaling with the CO 2 concentration levels, with a factor of below and above 2 for R95 and R99, respectively, from 2 × CO 2 to 4 × CO 2 . Fitting the precipitation extremes exceeding the 99.99 th percentile of PD with an extreme value distribution (i.e., GPD) shows significant increases over South China, North China, and Japan in CO 2 perturbed experiments, and over Central East China and South Korea in 4 × CO 2 .
Note that the CO 2 -perturbed simulations do not reach equilibrium completely as reflected by the long-term trend of global mean sea-surface temperature (Fig. S5). Thus, the results presented in this study combine both the fast (i.e., radiative forcing) and slow (i.e., sea-surface temperature) responses to CO 2 forcing. Previous studies revealed the competing effect of fast and slow responses in generating monsoon changes, which involves an enhancement associated with increased land-sea contrast in the former and a weakening accompanied with decreased land-sea contrast in the later (Shaw and Voigt 2015;Chen and Bordoni 2016). The intensified EASM indicates the dominant role of fast responses in our study. Disentangling the contributions of DY and TH in both fast and slow responses needs additional experiments forced by climatological SST as Chen and Bordoni (2016), which will be explored in future studies.
We acknowledge that the robustness of the findings needs to be assessed in other climate models. Nevertheless, this work serves as an example of how to evaluate the performance of monsoon circulation across different time scales in one of the highest resolution global coupled models. Understanding the responses of the EASM to different levels of CO 2 concentrations and the underlying physical mechanism can be beneficial to make better monsoon projection and mitigation strategy to future monsoon changes.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.