In summer 2022, the Yangtze River basin in China experienced an exceptional hot event. Its intensity reached the strongest level since complete meteorological observations began in 1961 (Figs. 1b,e). Due to reduced precipitation and prolonged high temperatures, the Yangtze River basin experienced an unusually dry summer (Figs. 1a,e); the dry conditions peaked around mid- to late August. The severity of this compound event, which occurred simultaneously with heat and dryness, also broke historical records (Fig. 1e). This compound event resulted in rivers and reservoirs drying up, forest fires, and consecutive electricity loads exceeding historical extremes, causing damage to ecosystems with severe impacts on the regional socioeconomics (www.mem.gov.cn/xw/yjglbgzdt/202209/t20220917_422674.shtml). The objective of this study is to investigate the impact of anthropogenic and natural factors on events similar to the case of 2022 in the Yangtze River basin.

The summer 2022 compound event in its different aspects. Spatial patterns of (a) SPI and (b) STI in East China to illustrate the event, where the black box denotes the Yangtze River basin. (c) The 500-hPa geopotential height (contours; dagpm) and anomalies (shading; dagpm), and the 588-dagpm contours (black solid line) and climatology (green line) for August 2022. (d) The scatterplot of domain-averaged SPI and STI over the Yangtze River basin; the year 2022 is shown as a red triangle, and the top and right sides represent the distribution (gray bar) and estimated PDF (black dashed line) for SPI and STI, respectively. The blue and red solid points represent the 2022 event. The PDF is smoothed by a kernel density function with 100 equally spaced points. (e) Time series of domain-averaged STI, SPI and SDHI from 1961 to 2022. (f) The return period of compound events (solid line) based on an optimal copula function (Frank copula).
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1

The summer 2022 compound event in its different aspects. Spatial patterns of (a) SPI and (b) STI in East China to illustrate the event, where the black box denotes the Yangtze River basin. (c) The 500-hPa geopotential height (contours; dagpm) and anomalies (shading; dagpm), and the 588-dagpm contours (black solid line) and climatology (green line) for August 2022. (d) The scatterplot of domain-averaged SPI and STI over the Yangtze River basin; the year 2022 is shown as a red triangle, and the top and right sides represent the distribution (gray bar) and estimated PDF (black dashed line) for SPI and STI, respectively. The blue and red solid points represent the 2022 event. The PDF is smoothed by a kernel density function with 100 equally spaced points. (e) Time series of domain-averaged STI, SPI and SDHI from 1961 to 2022. (f) The return period of compound events (solid line) based on an optimal copula function (Frank copula).
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1
The summer 2022 compound event in its different aspects. Spatial patterns of (a) SPI and (b) STI in East China to illustrate the event, where the black box denotes the Yangtze River basin. (c) The 500-hPa geopotential height (contours; dagpm) and anomalies (shading; dagpm), and the 588-dagpm contours (black solid line) and climatology (green line) for August 2022. (d) The scatterplot of domain-averaged SPI and STI over the Yangtze River basin; the year 2022 is shown as a red triangle, and the top and right sides represent the distribution (gray bar) and estimated PDF (black dashed line) for SPI and STI, respectively. The blue and red solid points represent the 2022 event. The PDF is smoothed by a kernel density function with 100 equally spaced points. (e) Time series of domain-averaged STI, SPI and SDHI from 1961 to 2022. (f) The return period of compound events (solid line) based on an optimal copula function (Frank copula).
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1
Data and methods
The observation data used in this study consist of monthly surface air temperature and precipitation with a resolution of 0.5° × 0.5° from 1961 to 2022. The data were retrieved from over 2,000 quality-controlled weather stations in China (Wu and Gao 2013), covering the Yangtze River basin (about 105°–123°E and 28°–35°N, black box in Fig. 1a). We use the geopotential height anomalies at 500 hPa relative to 1991–2020 climatology from the NCEP–NCAR reanalysis to describe the circulation background of this event.
We also use data from six CMIP6 models run under both natural and anthropogenic forcings (ALL) and natural-only forcing (NAT), respectively, during 1961–2020 (Eyring et al. 2016). As the 2015–20 segments of NAT simulations used natural forcing from the shared socioeconomic pathway 2-4.5 (SSP2-4.5) scenario, the historical ALL simulations, normally ended in 2014, were simply completed to 2020 with their corresponding SSP2-4.5 scenario runs. Our practice here is identical to Hu et al. (2023). We interpolated monthly precipitation and surface air temperature from different models onto the grid of observation dataset using bilinear interpolation. The basic information of models used is listed in Table S1 in the online supplemental material (https://doi.org/10.1175/BAMS-D-23-0149.2).
In this work, we are also interested in the anthropogenic influence on the likelihood of a hot event occurring while a configuration of dryness (X1 = x1) is fixed. This can be derived from the conditional probability density distribution:
The fraction of attributable risk (FAR = 1 – PNAT/PALL) is used to estimate the effect of anthropogenic influence on the occurrence probability of the 2022 event (Fischer and Knutti 2015; Stott et al. 2016). The terms PNAT and PALL represent the probability of the event under natural forcing and historical forcing, respectively. FAR can be interpreted as the fraction of the likelihood of an event that is attributable to the anthropogenic forcing. A bootstrap resampling was performed 1,000 times to estimate the 90% confidence interval (CI) of FAR as in Christidis et al. (2013).
Results
The majority of the study domain experienced extreme hot and dry conditions in 2022 (Figs. 1a,b). The atmospheric circulation provided favorable background for these conditions. The western Pacific subtropical high (WPSH) was anomalously westward shifted in June, July, and August (Fig. S1). The entire basin was controlled by a strong positive anomalies of geopotential height at 500 hPa (Fig. 1c), corresponding to the peak of dry conditions in August. Such an anomalous circulation pattern prevented cold air in the north from moving southward. Simultaneously, in the region controlled by anomalous high pressure, prevailing descending air inhibits precipitation and enhances warming through subsidence.
The average STI and SPI over the domain show significant correlation (r = –0.28, p < 0.05), resulting from the influences of the land–atmosphere interaction and the circulation anomalies. The main variability of SPI is at the interannual scale. For STI, however, there is a clear decadal variability, characterized by decreasing trends before the 1990s and increasing trends afterward. The domain-averaged STI of 2022 reaches 1.88, breaking the observation record since 1961 (Figs. 1d,e). The PDF deduced from the historical observations indicates that the probability of such a hot event is nearly 0.007 (Fig. 1d), corresponding to a return period of 150 years, with a 90% confidence interval of 59–721 years. The domain-averaged SPI of 2022 is –0.6 (Figs. 1d,e), indicating a severe meteorological drought with a probability of 0.1. Such a dry event is a 1-in-10-yr event (90% confidence interval: 1 in 7–17 years). In terms of a compound event combining both heat and dryness, a decrease of SHDI indicating an increase in the severity of a compound event was observed since the 1990s, and the severity of the 2022 compound event also largely broke the historical record (Fig. 1e). The best copula fit suggests that the return period of the 2022 hot and dry event is 753 years (533–1,151 years for the 90% confidence interval; Fig. 1f). It is to be noted that the longer return period of the compound event, relative to the single event, may be related to the heat–dryness negative correlation (Hao et al. 2022).
The performance of models in simulating STI and SPI over the Yangtze River basin is first assessed before conducting the attribution analysis. The CMIP6 models can well simulate the increases in the severity of compound events since the 1990s with a weak interannual variability (Fig. S2a). The CMIP6 models can also simulate the probability distribution of observed SPI and STI quite well, as shown in Figs. S2b and S2c. The distributions of both STI and SPI in the model and observation cannot be distinguished with the Kolmogorov–Smirnov (K-S) test (p value > 0.05). Then, we compare the probability of STI and SPI events similar to the one in 2022 in ALL and NAT simulations to quantify the impact of anthropogenic factor. For STI, the ALL forcing causes the PDF to generally shift to the right (Fig. 2a). Under NAT forcing, the probability of a heatwave event similar to the 2022 event is very low, estimated at 0.00035. However, under ALL forcing, this probability increases to 0.0087. In other words, human influence has increased the probability of the record-breaking heatwave in 2022 by nearly 25 times and explained 96% (90% CI: 90%, 99%) of the attributable risk for the event (Fig. 2d).

Attribution and projection results. The probability density function (PDF) of domain-averaged (a) STI, (b) SPI, and (c) bivariate return period of compound event for ALL forcing (red line) and NAT forcing (blue line); the selected copula function is rotated Clayton 270° for NAT and rotated and Tawn type 1 90° for ALL, respectively. (d) The best estimation (blue line) and 90% confidence interval (gray shadings) of the fraction of attributable risk (FAR) for SPI, STI, and the compound event. (e) The PDF of domain-averaged SPI, and (f) the bivariate return period of the compound event for historical forcing (blue line) and future SSP2-4.5 forcing (red line); the selected copula function is rotated Clayton 270° for future forcing. The black dashed lines and triangles denote the 2022 event.
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1

Attribution and projection results. The probability density function (PDF) of domain-averaged (a) STI, (b) SPI, and (c) bivariate return period of compound event for ALL forcing (red line) and NAT forcing (blue line); the selected copula function is rotated Clayton 270° for NAT and rotated and Tawn type 1 90° for ALL, respectively. (d) The best estimation (blue line) and 90% confidence interval (gray shadings) of the fraction of attributable risk (FAR) for SPI, STI, and the compound event. (e) The PDF of domain-averaged SPI, and (f) the bivariate return period of the compound event for historical forcing (blue line) and future SSP2-4.5 forcing (red line); the selected copula function is rotated Clayton 270° for future forcing. The black dashed lines and triangles denote the 2022 event.
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1
Attribution and projection results. The probability density function (PDF) of domain-averaged (a) STI, (b) SPI, and (c) bivariate return period of compound event for ALL forcing (red line) and NAT forcing (blue line); the selected copula function is rotated Clayton 270° for NAT and rotated and Tawn type 1 90° for ALL, respectively. (d) The best estimation (blue line) and 90% confidence interval (gray shadings) of the fraction of attributable risk (FAR) for SPI, STI, and the compound event. (e) The PDF of domain-averaged SPI, and (f) the bivariate return period of the compound event for historical forcing (blue line) and future SSP2-4.5 forcing (red line); the selected copula function is rotated Clayton 270° for future forcing. The black dashed lines and triangles denote the 2022 event.
Citation: Bulletin of the American Meteorological Society 104, 11; 10.1175/BAMS-D-23-0149.1
For SPI, anthropogenic forcing causes the entire PDF to shift to the left, indicating that anthropogenic forcing enhances dryness in the Yangtze River basin (Fig. 2b). Under NAT forcing conditions, the probability of a dry event similar to the 2022 event is 0.073. Under ALL forcing, this probability increases to 0.107. The risk of a dry event in the Yangtze River basin reaching or exceeding the intensity of 2022 under historical forcing is nearly 1.5 times larger compared to natural forcing, of which 32% (90% CI: 21%, 42%) can be attributed to human influence (Fig. 2d).
For the compound event (Fig. 2c), the 2022-like compound event has a high recurrence interval of nearly 3,500 years under natural world forcing, whereas it is a much more frequent 1 in 500 years in the actual world. Human activities have increased the likelihood of this record-breaking compound event by 7 times, of which 85.5% (90% CI: 66%, 96%) can be attributed to human activities (Fig. 2d). When compared to FAR for STI, it is evident that the occurrence of this compound event is predominantly driven by heat induced by anthropogenic influence.
It is worthy of note that such conclusions could not be directly extrapolated to a future scenario. Actually, we also performed the same analysis (results shown in Figs. 2e,f) with historical (1961–2020) and SSP2-4.5 (2041–2100) simulations of the same models as we used to deal with the pair NAT/ALL. As expected for STI, the probability for heat increases from 0.0087 in “historical” to 0.97 in “SSP2-4.5.” It is a little surprising to see that for SPI the probability of dryness decreases from 0.107 to 0.038, showing rainfall increases in the last half of the century (Fig. 2e). This is consistent with most other studies investigating future rainfall changes in the Yangtze River basin (Wang et al. 2014), a consequence of an enhanced East Asian summer monsoon (EASM) circulation. We suspect that the drying in ALL (compared to NAT) was mainly provoked by the important aerosol forcing, as demonstrated in Zhang et al. (2020) and Qian et al. (2009). The aerosol effects that contribute to the weakening of the EASM under ALL forcing have also been accounted for in CMIP6 models (Monerie et al. 2022). The likelihood of a compound hot and dry event, similar to the one in 2022, is expected to increase in the Yangtze River basin (Fig. 2f); the return periods from the historical and future simulations show a considerable shift from 1 in 500 years to 1 in 27 years. Nevertheless, these events are projected to be less severe in terms of dryness compared to the historical period, indicating that the compound event is primarily driven by heat, rather than dryness.
Dry conditions are generally favorable for hot events to occur, since the surface energy partition is altered. It is thus interesting to investigate the occurrence probability of a hot event like that of 2022 under a given dryness, and the fractional contribution of human activities to this conditional event. We found that human activities indeed increased the likelihood of a hot event (STI being greater than 1.88) given dry conditions (SPI being –0.6). The probability of this conditional event is 0.0052 for natural forcing. It increases to 0.017 for historical forcing (Fig. S2a). This gives a FAR of 70% (90% CI: 41%, 95.7%), implying that the anthropogenic influence explains 70% of the attributable risk for this conditional event. Additionally, we found that the best estimated FAR of the conditional event increases as STI increases (Fig. S2b), which means that human influence may lead to the occurrence of more intense hot events.
Conclusions
In the summer of 2022, an exceptional hot and dry compound event occurred in the Yangtze River basin in China. We found that anthropogenic effects have increased the occurrence probability of both heat and dryness as independent events, and compound events involving them together. Specifically, a fraction of 32% (FAR) for the individual dryness and a fraction of near 96% for the individual heat is attributable to anthropogenic forcing. For the compound event, the fraction of attributable risk in relation to human activities is 86%. Based on the attribution, the compound event is primarily driven by heat, rather than dryness. This is even more evident for the future, as the probability of dryness occurrence will decrease in this region, while the probability of a compound event will continue to increase.
Acknowledgments.
W. L. was funded by the National Key R&D Program of China (Grant 2022YFF0801704). W. L. and Z. J. were funded by the National Natural Science Foundation of China (42275184). We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6. We thank the climate modeling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access (https://esgf-node.llnl.gov/projects/cmip6/), and the multiple modeling funding agencies who support CMIP6 and ESGF.
References
Chen, Y., and Coauthors, 2019: Anthropogenic warming has substantially increased the likelihood of July 2017–like heat waves over central eastern China [in “Explaining Extreme Events of 2017 from a Climate Perspective”]. Bull. Amer. Meteor. Soc., 100, S91–S95, https://doi.org/10.1175/BAMS-D-18-0087.1.
Chiang, F., P. Greve, O. Mazdiyasni, Y. Wada, and A. AghaKouchak, 2021: A multivariate conditional probability ratio framework for the detection and attribution of compound climate extremes. Geophys. Res. Lett., 48, e2021GL094361, https://doi.org/10.1029/2021GL094361.
Christidis, N., P. A. Stott, A. A. Scaife, A. Arribas, G. S. Jones, D. Copsey, J. R. Knight, and W. J. Tennant, 2013: A new HadGEM3-A-based system for attribution of weather- and climate-related extreme events. J. Climate, 26, 2756–2783, https://doi.org/10.1175/JCLI-D-12-00169.1.
Eyring, V., S. Bony, G. A. Meehl, C. A. Senior, B. Stevens, R. J. Stouffer, and K. E. Taylor, 2016: Overview of the Coupled Model Intercomparison Project phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev., 9, 1937–1958, https://doi.org/10.5194/gmd-9-1937-2016.
Hao, Z., F. Hao, V. P. Singh, and X. Zhang, 2019: Statistical prediction of the severity of compound dry-hot events based on El Niño-Southern Oscillation. J. Hydrol., 572, 243–250, https://doi.org/10.1016/j.jhydrol.2019.03.001.
Hao, Z., and Coauthors, 2022: Compound droughts and hot extremes: Characteristics, drivers, changes, and impacts. Earth-Sci. Rev., 235, 104241, https://doi.org/10.1016/j.earscirev.2022.104241.
Hu, T., Y. Sun, X. Zhang, and D. Wang, 2023: Anthropogenic influence on the 2021 wettest September in northern China. Bull. Amer. Meteor. Soc., 104, E243–E248, https://doi.org/10.1175/BAMS-D-22-0156.1.
Monerie, P. A., L. J. Wilcox, and A. G. Turner, 2022: Effects of anthropogenic aerosol and greenhouse gas emissions on Northern Hemisphere monsoon precipitation: Mechanisms and uncertainty. J. Climate, 35, 2305–2326, https://doi.org/10.1175/JCLI-D-21-0412.1.
Nelsen, R. B., 2007: An Introduction to Copulas. Springer, 272 pp.
Qian, Y., D. Gong, J. Fan, L. R. Leung, R. Bennartz, D. Chen, and W. Wang, 2009: Heavy pollution suppresses light rain in China: Observations and modeling. J. Geophys. Res., 114, D00K02, https://doi.org/10.1029/2008JD011575.
Stott, P. A., and Coauthors, 2016: Attribution of extreme weather and climate-related events. Wiley Interdiscip. Rev.: Climate Change, 7, 23–41, https://doi.org/10.1002/wcc.380.
Vicente-Serrano, S. M., and J. I. López-Moreno, 2005: Hydrological response to different time scales of climatological drought: An evaluation of the standardized precipitation index in a mountainous Mediterranean basin. Hydrol. Earth Syst. Sci., 9, 523–533, https://doi.org/10.5194/hess-9-523-2005.
Wang, B., S.-Y. Yim, J.-Y. Lee, J. Liu, and K.-J. Ha, 2014: Future change of Asian-Australian monsoon under RCP 4.5 anthropogenic warming scenario. Climate Dyn., 42, 83–100, https://doi.org/10.1007/s00382-013-1769-x.
Wu, J., and X. J. Gao, 2013: A gridded daily observation dataset over China region and comparison with the other datasets (in Chinese). Chin. J. Geophys., 56, 1102–1111, https://doi.org/10.6038/cjg20130406.
Zhang, W., and Coauthors, 2020: Anthropogenic influence on 2018 summer persistent heavy rainfall in central western China [in “Explaining Extreme Events of 2018 from a Climate Perspective”]. Bull. Amer. Meteor. Soc., 101, S65–S70, https://doi.org/10.1175/BAMS-D-19-0147.1.