HTML
--> --> -->To handle these types of issues, we use the approach of clustering as described in Clark et al. (2021) in which daily June, July, and August fields of sea level pressure from the ECMWF ERA5 (Hersbach et al., 2020) reanalysis are combined into eight collections (hereafter “clusters”) using the k-means clustering technique. With new data now available from 2019 and 2020, the clusters of Clark et al. (2021) have been updated, using ERA5 output from 1979 onwards. Throughout this article, fields from ERA5 will be used as a representation of the observations of 2020 due to their spatial coverage and corresponding associated fields, for which gridded observations (of moisture fluxes, for example) are not yet readily available.
To examine how summers like 2020 might play out in a future warming world, we use model simulations from two ensembles of fully coupled ocean-atmosphere models. The first (denoted PPE), consists of a collection of 15 perturbed physics variants of the Met Office HadGEM3-GC3.05 model (Yamazaki et al., 2021) driven by observed greenhouse gas concentrations up to 2005 and following the (Relative Concentration Pathway) RCP8.5 (Moss et al., 2008) scenario to 2100. The second ensemble consists of simulations contributed to the CMIP6 (Eyring et al., 2016) archive, following a (Shared Socioeconomic Pathway) SSP5-8.5 (O'Neill et al., 2016) pathway, very similar to that of RCP8.5. Both are representative of a world of future rapid economic development and convergence. CMIP6 models, whose simulations were used for the analysis (based on data available at the time of doing the analysis) are shown in Table 1. Models were divided according to their resolution to gain an insight into whether their resolution had any effect on results. A single realization from each model was used.
Horizontal resolution | Model names | Total number of models | |
Sea level pressure (SLP) | Greater than approximately 250 km | CNRM-CM6-1-HR, NorESM2-MM, EC-Earth3, GFDL-CM4, CMCC-CM2-SR5, TaiESM1, MRI-ESM2-0, BCC-CSM2-MR, CESM2-WACCM, CNRM-CM6-1, CNRM-ESM2-1, GFDL-ESM4, CESM2, HadGEM3-GC31-MM, EC-Earth3-Veg | 15 |
Approximately 250 km | NorESM2-LM, INM-CM5-0, MPI-ESM1-2-LR, IPSL-CM6A-LR, FGOALS-g3, NESM3, MIROC6, ACCESS-CM2,ACCESS-ESM1-5,UKESM1-0-LL | 10 | |
Precipitation | Greater than approximately 250 km | As for SLP, but without CMCC-CM2-SR5, TaiESM1, CESM2-WACCM, CESM2, HadGEM3-GC31-MM, EC-Earth3-Veg | 9 |
Approximately 250 km | As for SLP, but without INM-CM5-0 | 9 |
Table1. CMIP6 models contributing to Fig. 6 and supplementary Figs. S6 and S7.
3.1. Circulation and moisture transport anomalies directly affecting China
Before performing any circulation clustering, we briefly examine some aspects of the 2020 summer, shown in Fig. 1. A westward zonal wind anomaly (compared to its climate, shown in Fig. S1e of the Electronic supplementary material, ESM) at 850 hPa was seen over the Philippines (Fig. 1b), the southern half of the South China Sea and Indochina, with the anomaly streamlines approximately parallel to the pressure anomalies. This wind anomaly matches the westward anomaly of the zonal water vapor flux seen in Fig. 1d, which would have transported moisture westwards from the South China Sea towards Indochina. The westward extent of this can be seen in Fig. S1f in the ESM. The anomalous northward water vapor flux (Fig. 1c) suggests that this anomalous moisture then passed into southern China, in addition to that which climatologically (see Fig. S1c) passes directly from the South China Sea. Over southern China, on the northern side of the surface pressure anomaly over the South China Sea, the eastward zonal wind was greater than usual (Fig. 1b), This eastward wind anomaly appears to have also manifested itself in the eastward water vapor flux (Fig. 1d) which would have resulted in anomalous moisture passing zonally along the mei-yu front, enhancing the rainfall over parts of eastern China under the front.Figure1. (a) ERA5 rainfall and (c), (d) vertically integrated northward and eastward water vapor flux (color shading) of 15 June to 15 August 2020 as a fraction of the corresponding 1979 to 2020 climatology. (b) Streamlines of 850 hPa wind speed anomalies (from 1979 to 2020 climatology) on the same days as in panels (a), (c), and (d). Green contours show equivalent sea level pressure anomalies (0.4 hPa intervals). Rainfall in (a) is masked (grey) where the climatology is less than 1.0 mm d?1. Pressure contours are masked (to avoid noise) in grid boxes with an altitude greater than 3000 m above sea level.
2
3.2. Initial cluster analysis of 2020
To compare the circulation of the 2020 summer to that of other years, we use the clustering approach described in section two. For reference, Fig. 2 shows the eight clusters of daily ERA sea level pressure from all June, July, and August (JJA) days from 1979 to 2020. They are very similar to those generated, in the absence of data, from 2019 and 2020 in Clark et al. (2021). A chart showing how the cluster occurrences have been distributed throughout the period is given in Fig. 3a. A rich variety of circulation types is seen in almost all years, clearly illustrating the complexity of issues mentioned early in section two. There are, however, some patterns. Clusters 1, 2, 3, and 4, for example, typically occur earlier in summer than clusters 5, 6. Clusters 7 and 8 typically occur late in summer. In some years, there are also highly unusual sequences of days when specific clusters persisted much longer than normal. Examples of these include cluster 5 in 1985, and cluster 7 in June 2004 and August 1998, 2010, and 2014. The June 2004 occurrences of cluster 7 appear to have been exceptional, given that this cluster is very rare in June.Figure2. ERA5 rainfall averaged over JJA (1979 to 2020) days in each cluster, as a fraction of overall (regardless of cluster) JJA climate. A value of 2.0, for example, shows rainfall twice the climate. A mask (grey) has been applied where JJA rainfall climate is less than 1.0 mm d?1. Dashed and solid contour lines show negative and positive anomalies, respectively, relative to the mean ERA5 sea level surface pressure anomaly (0.4 hPa intervals) in each cluster.
Figure3. (a) Identities of the ERA5 sea level pressure clusters of Fig. 2 on each June, July, and August day of each 1979 to 2020 year. (b) Statistics (green and brown) of the frequency of each cluster in JJA of each year, and in 2020 JJA (red). Trend values show the linear trend per 10 years, rounded to the nearest whole number.
Considering the summer of 2020, some details emerge quite easily. Clusters 5 and 6, for example, are almost entirely absent, even during August, when they typically occur. The difference between 2020 and the previous two years (2018 and 2019) in early August is striking. An unusual persistence of cluster 4 in late July 2020 is also shown. Cluster 4 was also more frequent than normal in both June and August. Figure 3b allows the frequency of each cluster in 2020 to be compared with the overall mean (which can be considered as a climatology) and the linearly calculated trend of the 42 year from 1979 to 2020. It shows that clusters 1, 2, and 4 were unusually frequent in 2020; for clusters 2 and 4, but not cluster 1, this matched the sign of the trend of their frequency.
To examine transitions between clusters, Fig. S2 in the ESM gives bar graphs showing the prevalence of clusters of the day immediately following days of those clusters (1 to 4) which were most frequent in 2020. Climatologically, individual days of a cluster are most likely to be immediately followed by a day of the same cluster. This pattern also occurred in 2020. In 2020 however, this occurred even more often than normal for days of clusters 1 and 4 (see Figs. S2a, d in the ESM). For these two clusters, though, the pattern was stronger in other years (2010 and 2005 respectively) when approximately 80% of the days of these two clusters were followed by an additional day of the same cluster. When days of 2020 did transition to a different cluster they tended to evolve almost exclusively into other types characterized by positive surface pressure anomalies over the South China Sea. There are hints that a weak cyclical pattern also occurred in 2020. When days of cluster 1 were not immediately followed by an additional day of cluster 1, they most often evolved into at least a day of either cluster 4 or 2, typically characterized by either strengthening or weakening, respectively, of the positive pressure anomalies over the South China Sea (Figs. 2b, d). Days of cluster 4 subsequently often evolved into those of cluster 2 and days of cluster 2 tended to evolve into days of either clusters 1 or 3. Interestingly, all of those which evolved into cluster 3 subsequently evolved into cluster 1.
2
3.3. Reconstructing 2020 rainfall from cluster rainfall and the role of 850 hPa zonal flow
In addition to indicating the mean sea level pressure in each cluster, Fig. 2 also provides the mean rainfall in each cluster, expressed as a fraction of the overall mean of all JJA days from the 42-years spanning 1979 to 2020. The spatial patterns of rainfall are clearly related to the magnitude of the pressure anomalies, for example, over the West Pacific between 20°N and 30°N in clusters 2 to 6 (see Figs. 2b–f). Although each cluster is, by definition, a collection of several hundred differing days grouped together solely by the sea level pressure patterns, it is interesting to examine how well the 2020 rainfall shown in Fig. 1a can be reconstructed, using a simple weighted average of the rainfall means of Fig. 2, where the weighting is given according to the frequencies of the clusters in 2020. However, although the reconstructed rainfall has similarities with that of Fig. 1a (see Fig. S3 in the ESM), with a moderate spatial correlation (0.31 and 0.43 depending on whether the rainfall used from each cluster was expressed in absolute terms or as a fraction of the overall climate), there are two very large differences. Firstly, the spatial variability of the reconstructed rainfall is much smaller, resulting from the use of the means of the cluster-dependent rainfall rather than information from the full distribution of days within each cluster. Secondly, the north-south dipole of the wetter and drier than climatology conditions, seen in Fig 1a, is also displaced. It is not surprising that none of the clusters of Fig. 2 show a spatial pattern of rainfall deficit in southeast China as pronounced as seen in Fig. 1a, considering how unusual the 2020 summer was. To further investigate this, Fig. 4 provides maps of the rainfall and 850 hPa zonal wind, during the 2020 summer in each of the four most frequent clusters of 2020, relative to what usually happens in each cluster. In all four, an enhancement of the characteristics of the rainfall in eastern China, shown in Figs. 2a–d, was seen (Fig. 4 left column). Immediately to the north and south of these rainfall enhancements, dipoles of opposing enhancement of the zonal wind were also seen (Fig. 4 right column), suggesting a propensity for additional (i.e. greater than typically expected), local scale cyclonicity. This appeared particularly prevalent in clusters 1 and 3 where meridional gradients of the zonal wind enhancements were remarkably sharp. As explained by Chou et al. (1990), the anomalous westerly wind, immediately to the south of the rainfall, would have favored further ascent and decent anomalies on the north and south sides of the low-level jet, respectively, enhancing the corresponding anomalies of rainfall.Figure4. ERA5 rainfall and 850 hPa zonal wind averaged over the days spanning 15 June to 15 August 2020 in each cluster. The rainfall and wind are shown as a fraction and anomaly respectively, of the cluster specific rainfall, compiled from all relevant 1979 to 2020 JJA days.
2
3.4. Grouping of the clusters
In light of the information from the above sections concerning the roles of clusters 1, 2, 3, and 4 in 2020, we now group the clusters for additional analysis. This can be done statistically using dendrograms applied to the pressure anomalies of the whole domain shown in Fig. 2. However, for examining the 2020 summer an approach focused on eastern China, where the rainfall anomalies were greatest, appeared to be more beneficial. The proximity of the West Pacific subtropical anticyclone to the region impacted is already known to have a large influence on the rainfall in the region (Rodwell and Hoskins, 2001). As highlighted by Lu (2001), Kim et al. (2002), and Preethi et al. (2017), westward extensions of the anticyclone strongly favor greater than normal rainfall over East Asia and the Korean peninsula with moisture transported from the South China Sea and West Pacific. Although primarily a feature at 500 hPa, the variability of the anticyclone also manifests itself through surface pressure anomalies over the South China Sea. For examining the 2020 summer in East Asia, it is thus appropriate to focus on this region for grouping the clusters. In clusters 1, 2, 3, and 4 (denoted Cluster Group CG1234), the surface pressure in the region is greater than normal (see Figs. 2a–d). In clusters 5 and 6 (denoted CG56), the pressure is less than normal (see Figs. 2e, f). The remaining clusters, 7 and 8 (denoted CG78), show no significant anomalies in the region but are characterized by positive anomalies over northwest China. By grouping the clusters in this way, the information of Fig. 3a becomes much clearer, as shown in Fig. 5a, with the 2020 summer dominated by days of the CG1234 cluster group. Of the 62 days, from 15 June to 15 August 2020, 59 of them were representative of this group of clusters and generally corresponded to the time when most of the rainfall occurred in eastern China. Furthermore, the second group (CG56) barely occurred at all in JJA 2020. From Fig. 5b, the cluster group occurrences during 2020, relative to the means, appear unrelated to the trends of the past 42 years.Figure5. As Fig. 3 but of the clusters of Fig. 2 grouped as follows: CG1234: clusters 1, 2, 3, and 4, CG56: Clusters 5 and 6, CG78: Clusters 7 and 8. SCS: South China Sea
Supporting the grouping of the clusters, Fig. S4 in the ESM shows anomalies of the 500 hPa geopotential height in the three groups. The position of the climatological WPSH is clearly westwardly displaced and stronger than usual during the days of the CG1234 group. During days of the CG56 group, it is weaker than normal and eastwardly displaced.
An alternative option of simply producing three clusters instead of eight clusters when clustering the daily sea level pressure fields was also attempted. However, this approach produced clusters with little spatial diversity (see Fig. S5 in the ESM). The distinctive negative pressure and associated positive rainfall anomalies over the South China Sea for example, in clusters 5 and 6 of the original 8 clusters, visible in Figs. 2e, f are difficult to identify in Fig. S5. The lack of diversity consequently limits their usefulness in examining extreme events such as those of the 2020 summer in eastern Asia.
2
3.5. Distributions of cluster group frequencies
Figures 6a, c, and e give distributions of the frequency of each of the three groups of clusters discussed above for the JJA period of each year. From the green bars of Fig. 6a, for example, it can be seen that during each of the 11 different years, there were between 50 and 55 days of the CG1234 anomaly group according to ERA5, roughly equating to between 54% and 59% of the 92 days in each annual period. Furthermore, in every year at least a third of the 92 days, were of the same group of clusters. In 2020, 73 of the 92 days (or 79%), according to ERA5 (the red dot in Fig. 6a), were of the group. In contrast, only two of the ERA5 JJA days of 2020 (Fig. 6b) were of the group of the two CG56 clusters, a group which typically occurs on 20 to 25 of the 92 days in JJA of each year. The group of clusters of high pressure anomalies over northwest China (CG78) typically occurs outside of the 15 June to 15 August period, both before and after. Therefore, their occurrence appears to have had little relevance during the mid-summer wet period of 2020.Figure6. (a), (c), and (e) Distributions of the number of JJA days in individual years, in the three groups of clusters (CG1234, CG56, and CG78), indicated in the row titles. Vertical green bars show ERA5 distribution. The red dot shows the number of days in the ERA5 reanalysis, of each cluster in 2020. PPE (see Yamazaki et al.) lines are of an RCP8.5 Relative Concentration Pathway. CMIP6 lines are of the SSP5-8.5 Social Economic Pathway. PPE and CMIP6 lines shown are ensemble means of distributions (distributions calculated separately for each simulation). (b), (d), and (f) Rainfall in the three cluster groups of (a), (c), and (e) as a fraction of overall climatological rainfall. A mask (grey) has been applied where JJA rainfall climate is less than 1.0 mm d?1. Green contours show equivalent sea level pressure anomalies (0.4 hPa intervals).
The light blue and pink lines of Figs. 6a, c, and e are the equivalent distributions from the PPE and CMIP6 (using models of all resolutions) ensembles described in section two, for the period 1979–2020. To produce these, each day of each simulation during the JJA period was compared (using Euclidean distance) to the centroids of the original eight clusters (shown in Fig. 2) to determine the most likely cluster of each day. The clusters of these days were then grouped in the same way as was done for those of the ERA5 reanalysis. The PPE and CMIP6 ensembles both appear to reproduce the shape and location of the ERA5 distributions remarkably well, especially when one considers that the 2006 to 2020 part of the simulations is forced according to a future scenario, rather than by observations. Given the observed impacts of the 2020 summer, the future frequency of days of the CG1234 group is of interest. With this in mind, we note that the prevalence of years with larger numbers of days of this group is, according to the PPE, greater for the future 2058 to 2099 time slice. For the CMIP6 models, the future prevalence appears unchanged. Additional plots (not shown), of the CMIP6 models according to their resolution (shown in Table 1) yield very similar results. A caveat here, however, is that only a single realization from each CMIP6 model was used for the analysis. The corresponding future rainfall of this time slice is discussed later.
The right-hand column of Fig. 6 gives the ERA5 rainfall for each of the three groups of clusters as a fraction of the climatological (i.e. regardless of cluster) rainfall . The first (Fig. 6b) has wetter than normal conditions over central and eastern China. The fractions, however, are no greater than 1.2. Much greater fractions are seen (Fig. 6d), for example, in the second group (CG56), over the west Pacific, associated with the cyclonic pressure anomaly. Interestingly, over eastern China, however, rainfall deficits, compared to climate are seen for this second group. In most years, the occurrence of days of this group of clusters would thus be expected to provide respite in eastern China from the rainfall occurring during days of the first (CG1234 anomaly) group. The absence of days in the second group in 2020, especially in July, would have reduced the likelihood of such respites occurring. Mean days of the third group (of CG78 clusters) are slightly drier than normal(Fig. 6f). However, these occurred in early June, before the main rainfall period in the regions most impacted by flooding, and then again in late August, after the rainfall had moved north. Their usefulness as a respite was thus very limited.
-->
4.1. Summers with circulation similar to that of 2020
We now examine how the prevalence of the CG1234 anomaly group is related to the rainfall and water vapor fluxes of the 2020 summer. We do this by focusing on the five years (1995, 1998, 2003, 2010, and 2013) in which the frequency of this group of clusters was most similar to 2020. The means of these are given in Fig. 7. The precipitation and moisture fluxes are expressed as a fraction (Fig. 7a, c, and e) of the overall climate (regardless of cluster) and are remarkably similar to those of Fig. 1. In particular, the spatial water vapor flux patterns show very clearly how the frequency of the CG1234 group could potentially be used as a proxy for the fluxes. Furthermore, in these five years, the northward and eastward water vapor fluxes over Indochina and southern China, respectively, are even greater than their own cluster-specific climate (see Figs. 7d, f), suggesting that the flux magnitudes are themselves related to the circulation cluster group frequency.Figure7. ERA5 precipitation and water vapor fluxes for days spanning 15 June to 15 August for cluster group CG1234 of clusters of surface pressure greater than normal over the South China Sea, averaged over the five years (1995, 1998, 2003, 2010, and 2013) in which the frequency of this group of clusters was most similar to that of 2020. AAF: As a fraction. Masks have been applied where the climates used in each plot were less than 1.0 mm d?1 in (a) and (b); and 0.05 kg m?1 s?1 in (c–f). Green contours show the corresponding surface pressure anomaly (0.4 hPa intervals).
2
4.2. Future precipitation in years of circulation similar to that of 2020
In a warmer world, the atmosphere, due to the Clausius-Clapeyron relation, has a greater capacity to hold and transport water vapor from moisture sources to destination regions (Held and Soden, 2006). In the absence of constraints on the moisture sources, rainfall in these destination regions is thus likely to increase. Due to this, and the impacts of the 2020 floods, we now examine how summer rainfall in years with a circulation similar to that of 2020 may change in a future time slice (2058 to 2099). This is useful since we have already seen indications of an increase in the frequency of these types of summers in the future in the PPE ensemble (Fig. 6a).The top row of Fig. 8 gives the mean rainfall of the five years in the 2058 to 2099 PPE time slice in which the frequency of the CG1234 cluster group is the most frequent. The values are expressed as fractions of the equivalent rainfall from the 1979 to 2020 PPE time slice. The identities of the years were determined independently for each simulation. These five years aretherefore indicative of how the conditions that were observed in 2020 summer could change in the future. Results are given as lower, median, and upper percentile estimates from across the PPE ensemble. Of the three estimates, confidence will naturally be greater in the median. However, the lower and especially the upper estimates should not be ignored because of their potential impacts should they occur. The median estimate (Fig. 8b) strongly suggests that future 2020-circulation-type summers will be wetter than their corresponding recent-past equivalents. From the upper estimate, there is a low (10%) risk of many locations even experiencing 1.5 times as much rainfall as in 2020. Considering that the values in Figs. 8a–c are expressed as fractions relative to the rainfall in years of circulation similar to the extremely wet 2020 summer, suggests that adapting to such increases could potentially be very challenging.
Figure8. Fractional anomalies of HadGEM3-PPE rainfall during CG1234 days of clusters of surface pressure greater than normal over the South China Sea between 15 June and 15 August for the five years of 2058 to 2099 in which these clusters were most frequent. (a) to (c) Lower, median, and upper percentile estimates (from across the PPE ensemble) as a fraction of the mean (calculated separately for each PPE simulation) of CG1234 days from the equivalent five years of 1979 to 2020 in the PPE. (d) to (f) Lower, median, and upper percentile estimates (from across the PPE ensemble) as a fraction of the mean of all days (calculated separately for each PPE simulation) of 1979 to 2020. Masks have been applied where the 1979 to 2020 values were less than 1.0 mm d?1.
The bottom row of Fig. 8 shows the future rainfall as a fraction of the overall 1979 to 2020 model climatology (calculated independently for each simulation). This allows an approximate visual comparison with Figs. 1a, 6b, and 7a. In almost all regions of China, rainfall is likely to exceed the overall 1979 to 2020 climatology during these years. A resulting expansion of the area impacted by the type of rainfall event which occurred in 2020 thus looks likely, presenting a further additional challenge.
A caveat, however, arises from how the estimates from across the ensemble are calculated. To produce these, the future rainfall fractions (of the 1979 to 2020 period) were ranked separately for each grid box. Spatial coherence between regions is thus not guaranteed. The simulation contributing to the greatest increase in one region, may not be the same elsewhere. Additionally, the PPE ensemble is small, and its range of responses may be skewed as a result of the parameters perturbed. As mentioned in section two, the PPE is also built from a single model. The range of estimates is thus likely to be narrower than that which might be obtained from an ensemble of PPEs built from multiple models. To quantify such estimates is beyond the scope of this study. However, it is possible to get an approximate, initial estimate from the inclusion of a single realization from different models. Figures S6 and S7 in the ESM give the equivalent estimates of Fig. 8, but from the CMIP6 models (see section two for their identities). The increases of the median and upper rainfall estimates are smaller than those from the PPE but their likelihood over eastern and northeastern China appears robust. The CMIP6 model resolution has a small effect (from comparing Figs. S6 and S7) with models of higher resolution having slightly larger increases.