Drivers behind the summer 2010 wave train leading to Russian heatwave and Pakistan flooding

Summer 2010 saw two simultaneous extremes linked by an atmospheric wave train: a record-breaking heatwave in Russia and severe floods in Pakistan. Here, we study this wave event using a large ensemble climate model experiment. First, we show that the circulation in 2010 reflected a recurrent wave train connecting the heatwave and flooding events. Second, we show that the occurrence of the wave train is favored by three drivers: (1) 2010 sea surface temperature anomalies increase the probability of this wave train by a factor 2-to-4 relative to the model’s climatology, (2) early-summer soil moisture deficit in Russia not only increases the probability of local heatwaves, but also enhances rainfall extremes over Pakistan by forcing an atmospheric wave response, and (3) high-latitude land warming favors wave-train occurrence and therefore rainfall and heat extremes. These findings highlight the complexity and synergistic interactions between different drivers, reconciling some seemingly contradictory results from previous studies.


INTRODUCTION
In summer 2010, a strongly meandering jet stream connected two simultaneous extreme weather events, the Russian heatwave, and the Pakistan flood, leading to severe societal impacts. By the end of July, beginning of August 2010, a series of heavy rainfall outpours hit Pakistan and the western Himalayan foothills region, causing the flooding of one-fifth of the country, severe damage to infrastructure, and taking the lives of 1700 people 1,2 . Simultaneously, western Russia was hit by a severe heatwave with anomalously warm temperatures from the end of June until mid-August 1 . The heatwave was accompanied by drought conditions and forest fires and may have caused the death of~15,000-55,000 people 1,3 . These events were characterized by a semi-permanent blocking anticyclone accompanied by air subsidence, adiabatic warming, reduced cloudiness, and soil drying 1,[4][5][6][7] . Below normal soil moisture limits evapotranspiration, reducing cloudiness and further exacerbating land warming 1,5 . Over Pakistan, a combination of a mid-tropospheric trough bringing relatively dry and cold air from the mid-latitudes, high sea surface temperatures (SST) over the Arabian Sea, and La Niña conditions in the tropical Pacific, allows the advection of moisture into the region 1,8 . This situation led to 3 weeks of intense rainfall over the region, with rainfall anomalies 300% higher than the climatology for the same period 1,2,7 . These two extreme events were connected by an unusually wavy jet stream 1 . Such wave trains can have a strong impact on surface weather conditions and weather extremes, and in particular on heatwaves, with possible consequences for simultaneous crop failures in the northern mid-latitudes 9 .
Thermodynamic warming of the atmosphere intensifies heatwaves and rainfall extremes 10,11 . Extreme rainfall events have increased in recent decades and model simulations project further increases due to thermodynamics, although uncertainties are larger for rainfall than for temperature changes [11][12][13][14][15] . The frequency of heat extremes over Europe has increased in recent decades at a pace that cannot be explained by atmospheric thermodynamics alone [16][17][18] . The effects of global warming on dynamics, including summer wave trains and blocking, are not well understood 9,[19][20][21][22][23][24] . Blockings and atmospheric Rossby waves could contribute to changes in the probability and persistence of extreme warm events in summer 21,25,26 . Some studies suggest that the frequency of occurrence of the mid-latitude wave pattern that characterized summer 2010 and led to the Russian heatwave has greatly increased since 1990 21 , though results are sensitive to the metric used and the analyzed region 27,28 .
SST, especially in the tropics, have a strong influence in shaping atmospheric waves in the North Atlantic/European [29][30][31] and North American sector in boreal summer 32 . SST and tropical convection can act as a source of circumglobal wave trains, thus influencing the mid-latitude circulation 29,[33][34][35][36] . In summer 2010, La Niña-like SST anomalies co-occurred with a warmer than the usual Indian Ocean and a negative Indian Ocean Dipole phase 37 . Model simulations of the 2010 Pakistan flood event show that these SST anomalies are instrumental in causing the circulation anomalies that led to the anomalous rainfall over Pakistan and northern India 8 . Moreover, the combined dynamical and thermodynamic response to anthropogenic-induced trends in SST patterns, shown in atmospheric models, increases the probability of occurrence of the Russian heatwave 38 .
Recent work suggests that high-latitude land warming can also affect the development and persistence of atmospheric wave patterns in boreal summer. Quasi-resonant amplification (QRA) refers to the amplification of synoptic-scale quasi-stationary free waves due to the presence of a mid-latitude waveguide in the zonal wind field 22 . QRA has been linked to the occurrence of amplified planetary waves and some recent heatwaves in the midlatitudes, including the Russian heatwave in 2010 39 . Anthropogenic climate change and in particular enhanced warming over high-latitude land areas might alter the characteristics of these waves by two possible mechanisms: (i) weakening of jet stream and storm tracks, which could, in turn, make wave trains more persistent 19,25 ; and (ii) favoring double jet states characterized by a confined sub-tropical jet that acts as an efficient waveguide 39,40 . These mechanisms may create more favorable conditions for the formation of waveguides and amplified wave patterns and thus for extreme weather 19,40 .
Soil moisture anomalies can act as a source of Rossby waves and thus represent another surface forcing that can impact largescale dynamics and induce or exacerbate heatwaves in the midlatitudes 41 . In northern mid-latitudes, model experiments have shown that dry soils can induce terrestrial heating anomalies which in turn can lead to a stronger wave response [41][42][43] . Soil moisture changes related to global warming could act to exacerbate the impact of heatwave events 44 . Summer soil moisture in extratropical regions is projected to decrease with climate change and such regional changes can substantially affect the climate in both hemispheres, with stronger and more robust impacts on surface temperatures than rainfall 45,46 . Future changes in soil moisture patterns may further exacerbate extreme events characterized by atmospheric circulation patterns similar to those observed in 2010 47,48 .
Here, we study the influence of SST patterns (together with global warming), high-latitude land warming, and early summer soil moisture conditions by quantifying their respective contribution to the 2010 Russian and Pakistan extremes. To do so, we use a very large ensemble of general circulation model simulations forced with observed SST anomalies and radiative forcings provided by the weather@home/climateprediction.com project. We use the global climate model of weather@home and a oneway nested 50 km regional climate model covering the South Asia region (region boundaries denoted in Supplementary Fig. 1). The regional model domain is much larger than the study area ensuring that any boundary effects in the region of interest are minimal. We quantify the relative contribution of these mechanisms to the probability of occurrence of the joint events, reconciling earlier findings on these interactions.

RESULTS
The 2010 Russian and Pakistan extremes A circumglobal wave, visible in observed meridional winds and geopotential height data (Fig. 1e, f), connects the extreme temperature in Russia (Fig. 1a) with rainfall anomalies in Pakistan (Fig. 1c). Figure 1b shows the probability density function (PDF) of standardized surface atmospheric temperature (SAT) anomalies over WRussia based on the observed 1987-2015 climatology (consistent with the available weather@home climatology). The observed mean value for this region is 293.1 K (=20.1°C), with an interannual standard deviation (s.d.) of 1.7 K. Figure 1d shows the PDF of standardized anomalies of observed Pakistan Rainfall for the 1987-2015 period. The mean value is 2.7 mm/day, with an s.d. of 0.9 mm/day. Both WRussia SAT and Pakistan Rainfall indices have an outlier in 2010, with both extremes exceeding 3 s.d. above the mean. Figure 2 shows a schematic of the model experiment setup together with their corresponding labels. The 2010 large ensemble provided by weather@home is labeled as Glob2010 for the global model and as Reg2010 for the regional model. The corresponding 1987-2015 model climatologies are labeled as GlobClim and RegClim, respectively. See the "Methods" section for further information on the model experiment setup and the definitions of indices.
Our model experiment shows the relative contribution of drivers to the extreme SAT and Rainfall of summer 2010 by conditioning on different potential drivers (Fig. 3). In GlobClim, concurrent events occur twice as frequently as expected for independent variables, supporting the significance of the role of an atmospheric wave train for their concurrent occurrence. In both Glob2010 and Reg2010 these events have clearly larger probability to occur, both individually and in concurrent mode. In 2010, exceedance probabilities increase approximately by a factor of 2 for WRussia SAT index and by a factor 4 for both the Pakistan Rainfall index and the concurrent event. Conditioning on drier soils (Glob2010 | soils and Reg2010|soilM) over western Russia in June and high-latitude land warming (Glob2010 | T65N and Reg2010 | T65N) in July-August, further enhances these probabilities by a factor of~2.5 for WRussia SAT index to a factor of~5 for the Pakistan Rainfall index.
A recurrent wave connecting the extremes Figure 4 shows the ensemble mean composites for SAT over Eurasia (Fig. 4a), Z300 (Fig. 4g), and V300 (Fig. 4h) from Glob2010, and for rainfall over the Indian subcontinent from both Glob2010 and Reg2010 ( Fig. 4c and e, respectively). A pronounced atmospheric wave-5 pattern and associated warm anomalies in Russia (shifted eastward as compared to the observed anomalies, black box in Figs. 1a and 4a), and increased rainfall over Pakistan up tõ 3 mm/day higher than the climatology are shown. Rainfall from Reg2010 ( Fig. 4e) shows the same pattern and magnitude as for Glob2010 (Fig. 4c), though with a more detailed small-scale structure close to the Himalayan Mountain range. In all three cases, the 2010 PDFs are shifted towards higher values compared to the climatological PDF (Fig. 4b, d, f). The PDFs are estimated from 649 temperature values for both Glob2010 and Reg2010, and 4930 (2871) precipitation values for GlobClim (RegClim), i.e. one value for each member and/or year. Temperature and precipitation data are standardized based on the mean and standard deviations of GlobClim and RegClim, respectively. GlobClim has an average SAT of 296.4 K, i.e. 3.3 K above the observed climatological mean for WRussia in the period 25.07-08.08, and an s.d. = 2.4 K. GlobClim (RegClim) has a mean rainfall value of 3.50 mm/ day (3.36 mm/day) and an s.d. of 0.59 mm/day (0.53 mm/day). To account for the model bias with respect to observations, we standardize the model's output by subtracting the mean and dividing by standard deviation. The observed SST pattern used as boundary conditions for the 2010 model experiments can be seen in Supplementary Fig. 2.
The number of extreme events for the WRussia SAT index increases from 10% (i.e. 90th percentile) to~17% in Glob2010, while for the Pakistan Rainfall index it increases from 10% to~34% in Glob2010 (43% in Reg2010), i.e. by about a factor 3-4. The percentage of ensemble members exceeding the 2010 threshold in the climatology is shown together with the percentage of ensemble members exceeding the 2010 threshold in the Glob2010 (or Reg2010) ensemble (see Fig. 4b, d, f). While the WRussia SAT index never reaches the observed (standardized) 2010 threshold, the Pakistan Rainfall indices jump from a 0.2% of occurrence in the climatology to a 0.8% and 1.4% in Glob2010 and Reg2010, respectively. Thus, there is a 4-to-7-fold increase in extreme rainfall events in this region linked to 2010 forcing conditions. Thus, the combination of 2010 SST patterns and radiative forcings strongly influences the appearance of these two interconnected extremes. All pairs of PDFs shown in Fig. 4b, d, and f are significantly different from each other using a Kolmogorov-Smirnov test and a significance level of α = 0.05 (see Supplementary Table 1).
The WRussia SAT 2010 event was so extreme that neither the climatology of the model nor the 2010 ensemble can reproduce values exceeding its threshold (Fig. 4b). This implies that the occurrence probability for such extreme temperatures in the W@H 2010 ensemble is lower than 0.2% (or a one in the 500-year event) indicating that even under present-day GHG forcing the event was extremely rare. These estimated return times are higher than climate. However, their definition was based on the monthly mean July temperatures whereas we focus on a shorter period of 2 weeks coinciding with the most extreme temperature anomalies 5 and rainfall extremes over Pakistan. . This shows the importance of 2010-SST patterns in forcing the circumglobal wave train behind the heat-and rainfall extremes. However, 2010 anomalies cannot be deterministically explained only by the 2010 SST field and radiative forcings. Next to these forcings, other effects from soil moisture variability, land-atmosphere dynamics, and model bias are likely to play a role and may determine the shift in the phase position of the atmospheric wave pattern. Figure 5 shows that a recurrent circumglobal wave pattern connects the Russia and Pakistan events, and therefore extremes in these regions are correlated 50 . Here, we define a wave pattern as recurrent when its spatial (wave number and phase position) and temporal (duration) characteristics are similar among pronounced wave-events in different ensemble members, similarly to the definition in Kornhuber et al. (2019) 50 . Note that anomalies for the sub-selections of Glob2010 and Reg2010 are compared to the ensemble mean of Glob2010 and Reg2010, and not compared to GlobClim or RegClim. We do this because we are interested in how different (pre)conditions within 2010 affect the occurrence of extremes, beyond what's driven by SSTs (as shown in Fig. 4). Figure 5b shows how the percentage of concurrent events (see the "Methods" section for the definition) changes between the climatology of the model and the 2010 run, with ã 4-fold increase (~8-9%) of 2010 concurrent events with respect to the climatology of the model (~2%). It is important to notice that, already in GlobClim, this percentage is~2%, thus twice as large as expected for independent extremes (0.1 × 0.1 = 1%). Note that we can calculate this value only for the global model ensemble, as the regional rainfall data for the climatological run come from a different ensemble of model simulations (batch) of initial conditions (see Supplementary Discussion 1 for further information). For concurrent events selected in Glob2010 and Reg2010 ensembles, this percentage jumps to~8% and~9%, respectively, showing how exceptional 2010 was. Composites of concurrent events (Fig. 5a, c-f) show patterns largely resembling those seen in the observations (Fig. 1) for the Atlantic-Eurasian sector. The same composites for SAT, Z300, and V300 but obtained with Rainfall from Glob2010 (Fig. 5c) are shown in the Supplementary Information (Supplementary Fig. 4). SAT over western Russia shows values up to~4 K higher than the 2010 ensemble mean (Fig. 5a), while rainfall over Pakistan shows values between 2 and 4 mm/day higher than the 2010 ensemble mean (Fig. 5c). In Z300 and V300 fields, a pronounced wave over the Eurasian continent is shown, with positive anomalies for Z300 centered over western Russia and strongest V300 anomalies, showing southerly winds to the west and northerly winds to the east of this positive Z300 anomaly (Fig. 5e, f).
The effect of high-latitude land warming and dry soils By conditioning on high-latitude land warming in the 2010 ensemble (Glob2010 | T65N or Reg2010 | T65N, see the "Methods"   section), the probability of occurrence of both extreme events further increases (Fig. 6). As for Fig. 5, anomalies are deviations from the Glob2010 (or Reg2010) ensemble. While we select for zonal-mean land warming that peaks around 65°N, two main warming anomalies appear south of that latitude (~50°N), i.e. over central Canada and from Eastern Europe to central Asia (Fig. 6a). WRussia SAT extremes increase from~17% for Glob2010 to~26% in Glob2010 | T65N, thus with an absolute increase of~10% (Fig.  6c). Pakistan rainfall shows an increase of the percentage of extreme wet events from~43% in in Reg2010 to~54% in Reg2010 | T65N (Fig. 6e). The corresponding composites of SAT for Glob2010 | T65N (exceeding the 90th percentile) and also exceeding the 90th WRussia SAT index (GlobClim) percentile (Fig. 6b) and of Rainfall for Reg2010 | T65N (exceeding the 90th percentile) and also exceeding the 90th Pakistan Rainfall index (RegClim) percentile (Fig. 6d) show positive SAT anomalies over western Russia paired with enhanced rainfall over the Himalayan foothills and the Indus basin. The PDF pairs shown in Fig. 6c and e are significantly different using the Kolmogorov-Smirnov test (see Supplementary Table 1).
Drier than usual early season (June) soil moisture conditions increase the probability of both 2010 extremes (Fig. 7). Figure 7a shows composites of June soil moisture for Glob2010|soilM compared to the Glob2010 ensemble mean. Composites for SAT (Fig. 7b) and Pakistan rainfall (Fig. 7d) from ensemble members with drier soils show a higher percentage of extremes (Fig. 7c, e). WRussia SAT extremes go from~17% of Glob2010 to~25% in Glob2010|soilM. The percentage of Pakistan Rainfall extremes increases from~43% in Reg2010 to~48% in Reg2010|soilM. This shows that dry soils in western Russia increase the probability of extreme heat over WRussia (i.e. its endemic region) but also influence extreme rainfall in Pakistan by strengthening the connecting wave train (see Supplementary Fig. 5). The effect on extreme rainfall in Pakistan is similar for Glob2010 ( Supplementary  Fig. 6) though the increase in the far tail is somewhat smaller probably due to the lower resolution. The higher resolution of the regional model improves topographic representation and convective activity and is, therefore, better suited for very extreme rainfall events. The PDF pairs are shown in Fig. 7c and e are significantly different using the Kolmogorov-Smirnov test (see Supplementary Table 1). Finally, we analyze how the concurrent extreme is influenced by either dry soils in June or contemporaneous high-latitude land warming. Figure 8 shows the PDFs for the 2010 events that are characterized by high-latitude land warming signals or dry June soil moisture conditions. The probability of concurrent extremes (shown in the upper-right gray-shaded corner) greatly increases in the 2010 ensemble to 8.3%, or about a 5-fold increase. When we condition on (contemporaneous) high-latitude land warming or dry soil moisture conditions (in June), the percentage of concurrent events increases to 12.3% and 10.9%, respectively, in Glob2010. This thus represents a 30-50% increase if compared to the full 2010 ensemble. Similar results are obtained for Reg2010 ( Supplementary Fig. 11).

DISCUSSION
Between late July and early August, an atmospheric wave train connected the western Russia heatwave with the Pakistan flood.
Using a very large model ensemble, we show that this wave can be reproduced by the weather@home model. Already in the model's climatology, this recurrent wave connects the two extremes, making the probability of having a concurrent extreme (i.e. surface temperatures over western Russia and rainfall over Pakistan both exceeding the threshold of the 90th percentile simultaneously) twice as likely compared to assuming independence between the two events. We quantify the influence of different previously proposed drivers 8,40,45  conditions in late spring and early summer over Russia represent a potential driver for enhancing both local surface temperatures and remote rainfall over Pakistan. Consistent with that, the probability of concurrent extremes becomes four times higher in the 2010 ensemble and grows even further when high latitude land warming or dry soils are considered (to about a 7-fold increase). This thus captures the combined, potentially nonlinear, effect of both 2010 SST and radiative forcings plus anomalous soil moisture conditions or high-latitude land warming conditions. Thus, the influence of soil moisture and high-latitude land warming, as determined here, is specific for the 2010 SST background state, and cannot necessarily be generalized to the entire climatology. Similar conclusions regarding the roles of SST, high-latitude warming, and soil moisture, can be drawn both for single and concurrent extremes, highlighting the important role of the wave train in those events. A characteristic feature of 2010 was the La Niña-like SST anomaly in the tropical Pacific. Warm SST anomalies in the Indian Ocean have also contributed to increased low-level moisture content, possibly fueling rainfall activity in the Indian summer monsoon (ISM) region 51 . La Niña-like SST patterns affect ISM rainfall by shifting the region of strong convection related to the Walker circulation cell closer to Southeast Asia, enhancing the ISM circulation cell 52,53 . The most active center of ISM rainfall is usually located over central India, in a region where the monsoon trough is located (~25°N, 80°E) 54,55 . However, Fig. 4 shows that in the 2010 ensemble mean, rainfall is displaced towards the Indus river basin, instead of over the monsoon trough region, as seen in observations (Fig. 2). Hence, the 2010 Pakistan rainfall anomalies cannot be understood only as a deterministic response to 2010 SST patterns and radiative forcings. Complementary analysis is performed on GlobClim for La Niña and El Niño years. Composites of La Niña minus El Niño years show that selecting high-latitude land warming or enhanced Pakistan Rainfall higher than 90th quantile of the full GlobClim ensemble leads to higher than usual SAT anomalies over western Russia and enhanced rainfall anomalies in Pakistan during La Niña years (see Supplementary  Figs. 12 and 13).
Extreme rainfall events on the western Himalayan foothills, such as the 2010 Pakistan flood, have been linked to a combination of a mid-tropospheric trough and associated dry and colder air intruding into the mid-to-low troposphere in the Indian subcontinent from higher latitudes combined with local orographic features 56 . Here, though we do not analyze the V-wind component in the lower troposphere, both observations and model simulations show a trough in 300 hPa geopotential heights located northwest of the Pakistan region (Figs. 1e and 5e, respectively). This trough is part of a mid-latitude wave train stretching from the eastern North Atlantic through Europe and western Russia towards the Indian subcontinent. Thus, model simulations can reproduce the atmospheric circulation pattern linking these two extreme events, and their concurrent probability is higher than that of independent events.
Previous studies showed that soil moisture can act as a Rossby wave source 43 , influencing the severity of temperature extremes in the mid-latitudes 57,58 . Here, we show that local soil moisture conditions can not only increase the probability of local heat extremes but also that of remote rainfall extremes, in this case, Pakistan. The link from drier-than-usual soils in western Russia in June to higher-than-usual temperatures in the same region is straightforward: low soil moisture values reduce evapotranspiration and cloudiness, further enhancing the warming 5 . The link between high temperatures over western Russia and enhanced rainfall over Pakistan can be explained by the connecting wave train. Soil moisture feedbacks thus likely strengthen the blocking over Russia by heating the atmospheric column 4,41 , generating a downstream response of enhanced cold air advection southwards towards Pakistan, which is instrumental in producing massive rainfall with a lead time of~one month (see the "Methods" section). We analyze the soil moisture effect with one month lead time is, contemporaneously, the blocking over Russia could intensify due to a large-scale dynamical process (e.g. highlatitude land warming), and lead to both intensified soil drying locally and southward cold air advection toward Pakistan. In this case, although the two extreme events would result as correlated, one would not have caused the other. Finally, taking into account high-latitude land warming in the Northern Hemisphere as defined by Mann et al. (2018) 40 , reveals that increased temperatures in high-latitudes also increase the probability of occurrence of both extremes. However, the exact causal chain here is less clear. On one hand, high-latitude land warming may be detected as a consequence of the wave pattern that transports additional heat northwards. Alternatively, high-latitude warming itself might change the characteristics of the wave pattern. For example, while the wave-train might be mostly driven by SSTs and local soilmoisture anomalies, its persistence might increase due to enhanced warming in high latitudes. More research is needed to disentangle these aspects.
How the potential drivers of weather extremes over Eurasia, as identified in this study, may change in the future is still an open question. The Indian summer monsoon precipitation is projected to increase under global warming 12,59 . The ENSO response to climate change is uncertain: while observations suggest a La Niñalike warming trend in central-western Pacific SST 37,60 , some studies indicate that in the long run, an El Niño-like pattern might dominate 61,62 . However, a recent study suggests that the La Niña-like trend might in fact be GHG forced and the El Niño-like tendency shown by future projections may be bias-induced 63 . Global warming may also intensify the summer drying of soils and the corresponding modulation of heatwaves and potentially wave trains [44][45][46] . Finally, high-latitude land warming is projected to increase under climate change 64 and thus may contribute to heat extremes in mid-latitude regions and to extreme rainfall events over the western Himalayan region. The concurrent nature of these two extremes, as shown in our experiment, might imply enhanced societal risks in vulnerable regions 9 .
The weather@home models have shown its suitability to conduct climate simulations and attribution studies on several weather extremes such as the 2010 Russian heatwave, Australian wildfires in 2018, extreme rainfall events in Bangladesh, the 2017 Fig. 4 Weather@home 2010 ensemble mean. Panel a: ensemble mean of SAT anomaly for the period 24.07-08.08 for Glob2010. Note that the apparent date mismatch relative to observations stems from the fact that weather@home uses a 360 days calendar. Panel b: PDF for standardized WRussia SAT index for GlobClim (solid black line) and for Glob2010 (solid orange line). Using standardized indices helps to correct for model biases for a better comparison with observation. Panel c: Same as for panel a but for Rainfall anomalies over the Indian subcontinent for Glob2010. Panel d: PDF for standardized Pakistan Rainfall index for GlobClim (solid gray line) and Glob2010 (solid blue line). Panels e and f: Same as panels c and d for Rainfall anomalies in Reg2010 and RegClim. Vertical lines in panels b, d, and f show the climatological 90th percentile threshold (black line) and the 2010 threshold (red line) together with the percentage of ensemble members exceeding each threshold in the climatology (upper value) and in the 2010 ensemble (lower value). Panel g: ensemble mean for Z300 anomalies from Glob2010. Panel h: Same as for Panel g but for V300 anomalies from Glob2010. Anomalies are calculated as the deviation from GlobClim and RegClim ensemble means. The stippling represents significant grid points at α = 0.05 and is calculated applying a Student-t Test. The smoothing of the curve is done using a Gaussian kernel to produce continuous density estimates.
Chinese heatwave, and real-time extreme weather event attribution 49,[65][66][67][68] . Despite the presence of a warm bias in mid-latitudes, the global model (HadAM3P) is reliable in most regions, and in terms of year-to-year variability in global temperature over land 69,70 . Moreover, in this experiment, weather@home, is able to reproduce the circulation pattern that characterizes these extreme events, thus giving additional confidence in the performance of the model.
In conclusion, we have shown that the weather@home model can reproduce the atmospheric circulation pattern which led to the G.Di Capua et al.  Panel a: composite of Glob2010 SAT over Eurasia for the period 24.07-08.08 (anomalies from Glob2010 ensemble mean for the same period) obtained by selecting those ensemble members that show both WRussia SAT (Glob2010) and Pakistan Rainfall (Reg2010) indices above the 90th percentile, i.e. the concurrent events. Panel b: PDF for WRussia SAT index in Glob2010 (solid orange line, top). PDF for Pakistan Rainfall index in Glob2010 (solid blue line, middle right) and scatter plot (middle left) of Pakistan Rainfall index on WRussia SAT for GlobClim (gray dots) and for Glob2010 ensemble (orange dots). PDF for Pakistan Rainfall index in Reg2010 (solid blue line, bottom right) and scatter plot (bottom left) of Pakistan Rainfall index on WRussia SAT for Reg2010 (orange dots). Black lines show the 90th percentile both in the PDFs and in the scatter plot and highlight the concurrent events in the top right gray-shaded panels in each scatter plot. Panels c and d: Same as for Panel a but showing composites of Glob2010 and Reg2010 Rainfall, respectively. Panel e: composite of Glob2010 Z300 for the period 24.07-08.08 (anomalies from Glob2010 ensemble mean for the same period) obtained by selecting those ensemble members that show both Glob2010 WRussia SAT and Reg2010 Pakistan Rainfall indices above the 90th percentile. Panel f: Same as for Panel e but for V300. The stippling represents significant grid points at α = 0.05 and is calculated applying a Student-t-test. The smoothing of the curve is done using a Gaussian kernel to produce continuous density estimates.
concurrent Russia-Pakistan extreme in summer 2010 and we have identified the influence of different drivers of those extremes. The climatological run of our modeling study shows that the probability of getting a concurrent event is almost twice as much as if the two events were independent. This probability increases 5-fold when the 2010 year is simulated, highlighting the importance of the 2010 La Niña-like SST. Moreover, contemporaneous high-latitude land warming further increases the probability both of single events and of concurrent events. Drier than normal local soil moisture features over western Russia are found to increase the probability of the extreme event with a time lag of approximately one month. Thus, the combination of SST anomalies, high-latitude land warming, and low soil moisture, favors the occurrence of the atmospheric wave train leading to the concurrent Russian heatwave and Pakistan flooding. These findings highlight the complexity and dependencies behind such persistent wave trains and high-impact extremes within our highly interconnected and dynamical climate system.

Observational data
To characterize the atmospheric conditions related to the Russian heatwave and the Pakistan flood, we use 1.5°× 1.5°surface atmospheric temperature at 2 m (SAT), meridional wind component at 300 hPa (V300) and geopotential

Weather@home/climatepredicition.net
The weather@home system runs on a volunteer distributed computing network provided by the climateprediction.net project. weather@home uses the HadAM3P atmospheric model developed by the UK Met Office Hadley Centre 73,74 . HadAM3P is forced by observed SST, sea ice fractions, well-mixed GHG concentrations, volcanic activity, and solar forcing and has a resolution of 1.25°× 1.875°6 9 . The model is also forced by observed aerosols (SO 2 and atmospheric dimethyl sulfide) but it does not include black carbon. The regional version, HadRM3P, uses HadAM3P as boundary forcing and has a finer spatial resolution (0.44°× 0.44°) for improved simulation of mesoscale features and orographic precipitation. For a detailed description of the weather@home model and for more information on its performance, biases and applications see further references 49,69,70 .

model experiments
The large ensemble consists of 649 ensemble members for both the HadAM3P (from now on referred to as Glob2010) and HadRM3P (referred to as Reg2010) models. Each ensemble member starts in December 2008 and runs till September 2010. From Glob2010, we analyze geopotential height at 300 hPa (Z300), zonal and meridional winds at 300 hPa (U300 and V300, respectively), the surface atmospheric temperature at 1.5 m (SAT), and soil moisture (soilM). Note that ERA-Interim SAT is calculated at 2 m, however, defining SAT at this slightly higher level is not expected to make a material difference. Rainfall is available from both Reg2010 and Glob2010 simulations. All variables are available at daily resolution for both Glob2010 and Reg2010 with the only exception of soilM, which is provided on a monthly time scale. To reduce the noise associated with precipitation data, 3-day running means of both Glob2010 and Reg2010 Rainfall output are taken. The climatology of both regional and global models is calculated over the period 1987-2015 (for a total of 29 years) and hereafter referred to as GlobClim and RegClim. Each year in GlobClim (RegClim) has 170 (99) ensemble members, leading to a total of 4930 (2871) years. Figure 2 shows a schematic of the model experiment setups together with their corresponding labels.
Both Glob2010 and Reg2010 ensembles are forced by observed SST and radiative forcing agents. The pattern of observed SST trends in the past decades, induced by the GHG forcing, resembles the La Niña SST anomaly pattern ( Supplementary Fig. 2) 37,60 . However, our model setup does not allow for a thorough assessment of the effect of GHG and SST forcing separately. In the following, we calculate anomalies by subtracting the Glob2010 (or Reg2010) ensemble average (thus removing the SST/GHG effect) and show both anomalies and ensemble means. Each  (7) with the total number of T-65N and SoilM events in Glob2010 (65 each). The smoothing of the curve is done using a Gaussian kernel to produce continuous density estimates.
Reg2010 simulation is driven by an associated Glob2010 simulation, but the RegClim simulations were driven by different global simulations not part of the GlobClim ensemble used here. Hence, while for Glob2010 and Reg2010 it is possible to select a set of ensemble members in Glob2010 and select the corresponding ensemble members in Reg2010 (or vice versa), this is not possible between GlobClim and RegClim (for more detailed information see Supplementary Discussion 1).

High-latitude land warming
From the Glob2010 and Reg2010 ensembles, we make two sub-selections of ensemble members that meet specific characteristics. The first subselection consists of those ensemble members showing enhanced warming over land in high-latitudes (~65°N) in the Northern Hemisphere. We define high-latitude land warming (abbreviated as T65N) following Mann et al. (2018) that introduces a metric based on surface air temperature to identify QRA events. Observed QRA events are diagnosed by peaks in zonal SAT profiles around~65°N, labeled as climatological QRA fingerprint (see Supplementary Fig. 14). Here, we calculate the contemporaneous (lead zero) correlation between the QRA fingerprint and the zonal SAT profile from the Glob2010 ensemble. Those ensemble members that are characterized by high correlation values (above the 90th percentile = r 90th = 0.72, where r min = −0.61 and r max = 0.90) are selected for the "high-latitude land warming" sub-selection (Glob2010 | T65N and Reg2010 | T65N).

Soil moisture
The second sub-selection of Glob2010 and Reg2010 ensemble members is applied using local soil moisture characteristics over western Russia. First, we calculate the monthly soil moisture (soilM) index by averaging over the region of the western Russia heatwave (45°-65°N, 25°-60°E) from the Glob2010 ensemble. Then, ensemble members that show soilM values lower than the 10th percentile in the month of June are included in these sub-selections as Glob2010|soilM and Reg2010|soilM. Thus, the soil moisture signal is detected one month in advance (one month lead) with respect to the target variable, which is detected at the end of July/ beginning of August.

Temperature and rainfall indices
An index for the western Russia heatwave is defined, both in observations and model experiments, by applying a spatial average of SAT over western Russia (45°-65°N, 25°-60°E). Similarly, the Pakistan flood is defined by averaging rainfall over the Pakistan/western Himalayan foothills region (25°-40°N, 65°-85°E). Following Lau and Kim (2011), all indices are temporally averaged in the period 24.07.10-08.08.10 (starting on the 24.07 instead of 25.07 since the model uses a 360 days calendar), obtaining one index in each region for each ensemble member. We will refer to these indices as to WRussia SAT and Pakistan Rainfall indices, respectively. To account for model biases both indices are standardized by removing the climatological mean and dividing by their respective standard deviation derived from the 1987-2015 climatology. Extreme heat events over western Russia are defined as ensemble members where the WRussia SAT index exceeds the 90th (climatological) percentile threshold. Similarly, extreme rainfall events over Pakistan are defined as those ensemble members where the Pakistan Rainfall index exceeds the 90th (climatological) percentile threshold. Concurrent events are defined as those ensemble members where both the WRussia SAT index and the Pakistan rainfall index exceed their respective climatological 90th percentile thresholds simultaneously.

DATA AVAILABILITY
The model output data for the experiment used in this study is freely available at the Centre for Environmental Data Analysis (https://catalog.ceda.ac.uk/uuid/ dae723e78f9f4a0b9712f352f0a0231d). Data can be cited as: Sparrow