Decadal variability of extreme high temperature in mid- and high-latitude Asia and its associated North Atlantic air–sea interaction

The decadal variability of extreme high temperature (EHT) in the mid and highlatitudes of continental Asia (MHLCA) and associated mechanisms were investigated. Observational analysis indicated that, after removing global warming, the first leading mode of the EHT events showed a meridional dipole pattern and had significant decadal variability. During the periods 1980–1994 and 2012–2019 (1995–2011), EHT events were more frequent and stronger (less frequent and weaker) in western-central Siberia than normal, whereas they were less frequent and weaker (more frequent and stronger) in the wide area to the south of Lake Baikal. Further Observational and CESM (Community Earth System Model) results based analysis suggested that decadal change in air–sea interaction in the North Atlantic play an important role in shaping the decadal variability of EHT events in MHLCA. On decadal timescales, meridional negative–positive–negative sea surface temperature (SST) anomalies in the North Atlantic and their coupled positive North Atlantic Oscillation can trigger stronger wave activity flux and cause a significant anticyclonic–cyclonic teleconnection wave train in the troposphere over the mid and high latitudes of the Eurasian continent. As a result, an anomalous high-pressure center is evident in western-central Siberia, and thereby reducing total cloud cover and causing stronger solar heating. Thus, EHT events become more frequent and stronger. On the contrary, for the areas south of Lake Baikal, EHT events are weakened by corresponding cyclonic circulation anomalies. Additionally, the CAM5 (Community Atmosphere Model version 5) experiments suggested more important driving role of the decadal North Atlantic tripole SST anomalies in the abovementioned processes.


Introduction
In the past 100 years, the global surface temperature has increased by nearly 1 °C, as reported by the intergovernmental panel on climate change fifth assessment report (IPCC 2013).Particularly, from the middle of the 20th century to the present day, global warming has been accelerating (Xu et al. 2018;Cheng et al. 2019), and with it has come more extreme events such as extreme high temperature (EHT) events, extreme heavy precipitation, and exceptionally strong tropical cyclones, occurring in many parts of the world (Rahmstorf and Coumou 2011;Coumou and Rahmstorf 2012;Donat et al. 2013).More seriously, since the beginning of this century, record-breaking monthly mean temperatures have increased worldwide (Coumou et al. 2013).At the same time, record-breaking heat waves struck Europe in 2003 (Schär et al. 2004), Greece in 2007 (Founda and Giannakopoulos 2009), Russia in 2010 (Barriopedro et al. 2011), eastern China in 2013(CMA 2014), southern China in 2018 (Deng et al. 2020), and Europe again in 2019 (Xu et al. 2020), to name just the main examples.These EHT events caused major damage to ecosystems and society.Taking the 2003 European heat wave as an example, it culminated in 14,800 heat-related deaths in France alone, and also caused heavy casualties in other European centuries such as Italy, Spain, Portugal and the United Kingdom (Bouchama 2004).In the recent 2013 China heat wave, reported heatrelated illnesses exceeded 5700 (Gu et al. 2016).Therefore, issues related to EHT have received widespread attention.
Numerous studies suggest that anthropogenic global warming is one of the most important factors resulting in higher frequency and more severe EHT on the global scale (Griffiths et al. 2005;Christidis et al. 2011;Morak et al. 2013;Dong et al. 2016;Kim et al. 2016;Dong et al. 2017).At regional scales, "detection and attribution" studies also indicate the important role played by human activities in increasing the frequency and severity of EHT events (Zwiers et al. 2011;Christidis and Stott 2016;Yin et al. 2017;Chen and Dong 2019;Seong et al. 2021).In the future, in a world that has warmed by 1.5 °C, 2.0 °C or more, human activities will cause more higher-intensity EHT events on both global and regional scales (Chen and Sun 2015;Fischer and Knutti 2015).
In fact, regional EHT events are the result of the combination of natural internal variability of the climate system and anthropogenic emissions of greenhouse gases (Schär et al. 2004;Meehl et al. 2007;Shiogama et al. 2014;Zhou et al. 2014;Sun et al. 2016).Along with mean warming, enhanced variability can also result in more frequent and intense EHT events, and even record-breaking heatwaves (Barriopedro et al. 2011).Thus, natural internal variability also plays an important role in causing regional EHT events (e.g., Wang et al. 2013;Diao et al. 2018).Horton et al. (2015) suggested that both dynamic and thermodynamic effects can contribute to observed changes in temperature extremes.At most locations, warm extremes are associated with positive 500-hPa geopotential height and sea level pressure anomalies slightly downstream and negative anomalies farther upstream (Loikith and Broccoli 2012;Xu et al. 2019b).EHT events are often accompanied by specific atmospheric circulation patterns, less precipitation, and nearby prolonged hot surface conditions (e.g., Fischer et al. 2007;Lau andNath 2012, 2014).In inland area, anomalous circulations can advect hot and dry air, and thereby causing the regional EHT events.An anomalous high-pressure center and anticyclone near the surface can also lead to clear sky and thereby greater solar heating, which together are believed to be the main physical driver of hot extremes (e.g., Meehl and Tebaldi 2004;Tomczyk and Bednorz 2016;Cowan et al. 2017;Luo and Lau 2017;Adams et al. 2021).In addition, the descending branch within the column of the anomalous anticyclone can enhance local adiabatic heating, which also contributes to localized EHT events (Li et al. 2021).Within these processes, soil moisture-temperature feedback is usually involved (Stefanon et al. 2014).
Anomalous sea surface temperature (SST) is one root cause that can lead to large-scale circulation anomalies or local anticyclonic conditions, and ultimately to EHT events.Based on observational and modeling results, Arblaster and Alexander (2012) suggested that extreme maximum temperatures are significantly cooler over Australia, southern Asia, Canada and South Africa during strong La Niña events, as compared to El Niño events, and significantly warmer over the contiguous United States and southern South America.Additionally, tropical Pacific and Indian Ocean warm SST anomalies can influence EHT events over the southern part of China (Hu et al. 2013;Wang et al. 2014;Deng et al. 2019).Even decadal SST anomalies, such as the Pacific Decadal Oscillation and Atlantic Multidecadal Oscillation, have been shown to be closely related to EHT events-specifically, in North China (e.g., Zhang et al. 2020a;Zhu et al. 2020;Yang et al. 2021).
During the period of accelerated global warming in the past few decades, the Arctic region has shown the most significant warming owing to the so-called Arctic amplification effect (Screen and Simmonds 2010).In fact, the mid and high latitudes of continental Asia (hereafter abbreviated to MHLCA), which is adjacent to the Arctic, also show a stronger warming rate compared to other regions Fig. 1 in Xie et al. (2015); and logically, such marked increases in mean surface air temperature (SAT) must lead to similarly large increases in the frequency and severity of EHT events.Moreover, under the future warming scenarios we can expect with unchecked or insufficiently constrained global warming, these increases in the frequency and severity of EHTs event will cause greater impacts on local ecosystems and human society (IPCC 2021).In such a region with the most significant future warming, will natural variability also have an impact on EHT events and superimpose on warming to cause more severe extreme events?To better understand the likely future changes in EHT events in MHLCA, it is particularly important to clarify the mechanism of EHT variation in this area during recent decades.In this study, we will remove the influence of global warming and explore the internal variability factors affecting the decadal variations of extreme high temperature in this region based on 40-year daily observation data and our numerical simulation data.
The rest of this paper is organized as follows: The data, numerical model and methods used in this study are described in Sect. 2. In Sect.3, the observed spatiotemporal characteristics of EHT in MHLCA, as well as the observed and simulated mechanisms, are investigated.Conclusions are presented and discussed in Sect. 4.

Data, numerical model and methods
The daily maximum temperature (Tmax) from the Climate Prediction Center (CPC) global temperature data provided by the National Oceanic and Atmospheric Administration/ Oceanic and Atmospheric Research/Earth System Research Laboratories Physical Sciences Laboratory (https:// psl.noaa.gov/ data/ gridd ed/ data.cpc.globa ltemp.html) was used to calculate the intensity and frequency of EHT for the period 1980-2019.The monthly SAT, which was calculated based on the CPC global temperature data, was used to examine changes in mean SAT in boreal summer (June, July, and August).The observed SST from the hadley center monthly sea surface temperature dataset (Rayner et al. 2003), monthly precipitation in the Climate Research Unit (CRU) data set (Harris et al. 2020), and monthly ERA-Interim data from the European Centre for Medium-Range Weather Forecasts (Dee et al. 2011) are used to investigate related climate change during the period 1980-2019.
The Community Earth System Model (CESM) version 1.2.2, which was developed and is maintained by the National Center for Atmospheric Research, was used to carry out a long-term simulation.This model is a coupled global climate model.It has seven components, including atmosphere, ocean, land, sea ice, land ice, river runoff, and ocean waves.They are coupled by a coupler that coordinates these components' temporal evolution and passes information between them (more detailed information can be found at https:// www.cesm.ucar.edu/).A long-term simulation was performed with the "B1850C5" configuration, which means that the external forcing agents had no year-to-year variations.The values for greenhouse gas concentrations, aerosols, and ozone were determined based on the settings from pre-industrial control simulations in Phase 6 of the CMIP (coupled model intercomparison project), and were fixed at preindustrial levels (1850).The horizontal resolution was set to "f19_g16", which means that the horizontal resolutions of the atmosphere and ocean components were approximately 2° (1.9° in latitude and 2.5° in longitude) and 1°, respectively.The initial states used in the experiment were obtained from the b40_1850_c02c_76jpf files, which contain the initial fields for the pre-industrial simulation provided by the National Center for Atmospheric Research (NCAR) through long-term simulations.In this study, we further continuously integrated the 200 model years.The first 100 model years were used as spin-up to allow the model to reach an equilibrium state, accounting for changes in the computing environment and subtle parameters.We analyzed the final 100 years.Analysis of this long-term simulation can avoid influences from external forcings on the related mechanism and processes.In addition, two simulations using the Community Atmosphere Model version 5 (CAM5), which is the atmospheric component of the CESM, were performed to examine the influence of anomalous SST on EHT in MHLCA.More detailed information can be found in the following simulation result section.
In this study, the EHT intensity was defined as the maximum value of daily maximum temperature in summer of that year.Besides, the EHT frequency was defined as the number of days when the daily maximum temperature in summer of that year was greater than the 90th percentile value in the base period .The 90th percentile value was calculated using the same method as in Bonsal et al. (2001).Take one grid for example, there were 2760 Tmax values (30 years × 92 days in each summer), and they were first sorted in ascending order: X 1 , X 2 ,…,X 2760 .The probability P that a random value is less than or equal to the rank of that value X m is estimated by P =(m − 0.31)/(N + 0.38).N = 2760 and the Tmax representing the 90th percentile was linearly interpolated between the 2484th-ranked value (corresponding to P = 89.98%)and 2485th-ranked value (corresponding to P = 90.01%).In this study, to remove the effects of global warming on the EHT over MHLCA, we detrended the EHT intensity, the EHT frequency and associated observed meteorological variables before analysis.Referring to Zhang et al. (2018), the linear trends were removed for all the observational data to reduce the anthropogenic and other external forcings' influences.For statistical analysis, we mainly used the empirical orthogonal function (EOF) analysis, regression analysis and composite analysis, and the statistical significance of the composite analysis was assessed by applying the Student's t test.

Observed decadal changes in EHT
First, we made an EOF analysis of the detrended anomalies in the frequency and intensity of EHT events during the period 1980-2019.As shown in Fig. 1a and c, the first leading modes of the frequency and intensity of EHT were almost the same, showing meridional dipole modes.These two first leading modes can explain 17.4% and 14.8% of their total variances.Their positive and negative centers were located in western Siberia and the wide area to the south of Lake Baikal, respectively.It suggested that EHT experienced opposite trends in these two areas during the period 1980-2019.Their principal components (PCs) were also consistent during this period (Fig. 1b and d).On the decadal timescale, the two PCs both turned into a negative phase around 1995 and rebounded into a positive phase in 2012.This indicated that the EHT in MHLCA presented significant decadal variabilities.
The result of a seven-year moving t test further confirmed the decadal shift of EHT and the abovementioned jump points (i.e., the years 1995 and 2012, Fig. 2a).Accordingly, we chose the years from 1980 to 1994 and from 2012 to 2019 as the positive phases and the years from 1995 to 2011 as the negative phase.The composite results were consistent with the above EOF analysis.In western Siberia, the intensity of EHT was stronger and its occurrence was more frequent during the positive phases, whereas EHT was weaker and less frequent during the negative phase (Fig. 2b and c).On the contrary, the EHT in the mid-latitude region south of Lake Baikal showed opposite variations.In the following sections, we explored the mechanisms underpinning the decadal variabilities of EHT during the past 40 years.

Anomalous atmospheric circulation and North Atlantic SST
Regional atmospheric circulation was very different between the positive and negative phases.As shown in Fig. 3a, anomalous downward vertical motion can be found over western-central Siberia during the positive phases, whereas anomalous upward vertical motion was mainly present south of Lake Baikal and in some eastern Siberian areas.Correspondingly, significant anticyclonic and cyclonic circulation anomalies were evident in these two regions, respectively (Fig. 3b).Particularly at north of Lake Baikal, sea level pressure anomalies were increased, indicating significant anticyclonic conditions there (Fig. 3c).Due to anomalous vertical motion over MHLCA, the total cloud cover was significantly reduced over western-central Siberia (Fig. 3d).As a result, surface net shortwave radiation was increased in western-central Siberia, and dominanted regional surface heat fluxes (Fig. 4a).Finally, the SAT increased significantly in western-central Siberia (Fig. 3e).In the clear-sky, surface net shortwave radiation didn't show significant changes (Fig. 4b).It suggested that reduced total cloud cover and associated changes in solar radiation were key factors leading to local warming.For some areas south of Lake Baikal, the processes were opposite.Therefore, the anomalous spatial patterns in local circulation and related radiation were consistent with changes in SAT and the decadal variability of EHT events (Fig. 2b, c).This suggested that local circulation anomalies between positive and negative phases was the main reason behind the decadal variability of EHT events in MHLCA.On a larger scale, a significant anticyclonic-cyclonic teleconnection wave train in the troposphere can be seen in the mid and high latitudes of the Northern Hemisphere (Fig. 5).Over the Eurasia, the anomalous wave train is similar with the British-Baikal Corridor (Xu et al. 2019a), and the local circulation anomalies in MHLCA was one part of this teleconnection wave train.Analysis of the wave activity flux (WAF) suggested that this tropospheric teleconnection wave train probably propagated from the North Atlantic to the high-latitude Asian region.And in the source area of the wave train, positive North Atlantic Oscillation (NAO) was evident (Fig. 5a and b).
Therefore, we further examined the changes in North Atlantic SST between the positive and negative phases.As shown in Fig. 6, during the positive phases, the North Atlantic showed a meridional tripole SST anomaly pattern.Significant negative SST anomalies were evident in the tropics and north of 45ºN in the Atlantic, whereas positive SST anomalies lied between in the mid-latitude Atlantic.Such tripole SST pattern was well coupled with overlying positive NAO-like criculation.
Many previous studies have examined the decadal and longer scale North Atlantic SST variations and their role in modulating climate variabilities over the North Atlantic and Eurasian continent (e.g., Kushnir 1994;Schlesinger and Ramankutty 1994;Sutton et al. 2003;Lu et al. 2020).Actually, the tripole anomaly pattern was a dominant mode of North Atlantic SST variabilities after removing the linear trend.As shown in Fig. 7, this first leading mode can explain 28.2% of the total variance.Including annual variabilities, its PC showed that this tripole SST pattern had significant decadal variations during the past 40 years.This shifting also occurred in the years 1995 and 2012, which were the same years as for the decadal changes in EHT in MHLCA, suggesting there was a close relationship between them.
Such decadal North Atlantic SST pattern and anomalous NAO-like pattern were likely to be the root cause of the wave train in the troposphere, and thereby the decadal changes in EHT in MHLCA.

Observed remote influences of North Atlantic SST
To highlight the influences of the meridional tripole SST anomalies on the downstream climate, we used a summer North Atlantic tripole index (NATI) in the following analysis.The NATI was defined as (Fig. 7) These three NATI regions were chosen according to the leading mode (Fig. 7) and decadal change in the North Atlantic SST pattern (Fig. 6).The interannual and decadal variabilities of NATI were consistent with the PC of first leading mode of North Atlantic SST anomalies.The correlation coefficient between them was 0.74 (p < 0.01).
Figure 8a and b showed the wind fields at 300 hPa and 500 hPa regressed upon NATI.An anticyclone-cyclone-anticyclone-cyclone-anticyclone pattern was evident in the troposphere from the North Atlantic to central Siberia.This circulation pattern showed a significant wave-train feature,  the abovementioned wave train from the North Atlantic to Siberia.Consequently, anticyclonic circulation anomalies led to negative total cloud cover anomalies over westerncentral Siberia, thereby resulting in higher SAT and more frequent and stronger EHT events there (Fig. 9).This result confirmed the observed decadal linkage between the air-sea interaction in the North Atlantic and increased EHT events in western-central Siberia.
In addition, compared with composite analysis (Fig. 4), the wave train signals were more pronounced in the regression map.It suggested that the eastward wave train probably has more significant impacts on climate in MHLCA on an interannual timescale.Of course, a variety of other factors (such as anthropogenic and natural external forcings, etc.) in the observation might have influences on the above process, and thus can weaken the associated signals on the decadal timescale.More importantly, the length of observational data also limited the analysis of decadal and longerterm changes.Therefore, parallel analysis using a long-term simulation with no variations in external forcing agents was necessary.

Simulated linkage between North Atlantic climate and EHT
In the long-term CESM simulation, the meridional tripole SST pattern was also the first leading mode of the North Atlantic (Fig. 10a).In its positive phase, positive SST anomalies were evident in the mid-latitude North Atlantic, which were surrounded by negative SST anomalies.Two centers of the negative SST anomalies were located in the Labrador Sea and subtropical North Atlantic, respectively.This simulated tripole SST pattern resembled the observed one.At the same time, its PC and simulated NATI were consistent and showed significant decadal variabilities (Fig. 10b, c).Particularly, in the second half of this longterm simulation, the 9-year running averages of both the PC and NATI were almost always positive values during model years 46-57 and 94-100, whereas they were almost always negative values during model years 58-93.Thus, focusing on these model years, we examined the corresponding differences in circulation and temperature over the Eurasian continent, and especially the middle and high latitudes of Asia, between different decadal phases of the North Atlantic tripole SST pattern.
As shown in Fig. 11, a similar anomalous circulation pattern can be found in the troposphere over the North Atlantic and Eurasian continental region.Positive NAOlike pattern was evident over the North Atlantic.At the same time, a significant anticyclonic-cyclonic teleconnection wave train propagated from over the island of Newfoundland to Europe, and finally to western-central Siberia, which resembled the results based on the observational data (Figs. 5 , 8).The anomalous tropospheric WAF was also the same, demonstrating eastward propagation characteristics of anomalous circulation signals.Therefore, on the decadal timescale in both the simulation and observation, corresponding to the decadal meridional negative-positive-negative SST pattern in the North Atlantic and positive NAO-like pattern, anomalous anticyclonic circulation was evident in the troposphere over westerncentral Siberia.As a result, the simulated total cloud cover decreased, ultimately leading to increased SAT in western-central Siberia (Fig. 12).The abovementioned simulated process was almost the same as the observed one, which further confirmed the close linkage between In fact, a large number of studies have focused on the air-sea interaction process in the North Atlantic.On the interannual scale, the climate variations over the North Atlantic are related to changes in the NAO (Glowienkahense 1990;Hurrell and VanLoon 1997;Rogers 1997;Frankignoul et al. 2001;Hurrell and Deser 2009).Through large-scale air-sea interaction, the interannual tripole SST pattern in the wintertime North Atlantic Ocean was suggested to be driven by the NAO (e.g., Cayan 1992;Battisti et al. 1995;Marshall et al. 2001aMarshall et al. , 2001b;;Zhou et al. 2006).However, on decadal and longer timescales, the ocean was thought to significantly contribute to climate variability (Bjerknes 1964;Sutton and Hodson 2003).In the North Atlantic, the oceanic dynamics also played an important role in generating decadal SST anomalies (Eden and Jung 2001;Marshall et al. 2001aMarshall et al. , 2001b)).In addition, changes in North Atlantic sea ice is also suggested to have an impact on decadal SST variability (Deser et al. 1993).Focusing on the North Atlantic tripole SST pattern in the present study, Häkkinen (2000) indicated the importance of oceanic dynamics in shaping the decadal variation.Gulev et al. (2013) also suggested that, on the decadal timescale, the North Atlantic Ocean can drive the turbulent heat flux and force the atmosphere.Of course, both atmospheric forcing and ocean dynamics may play important roles in the formation of such decadal tripole SST anomaly pattern (Arthun et al. 2021).In this study, to examine the causes in shaping the wave train and decadal change in EHT in MHLCA, the observed summertime turbulent heat fluxes (sensible plus latent heat fluxes) over the Northern Ocean were diagnosed (Fig. 6b).Similarly, tripole turbulent heat fluxes can be seen over the North Atlantic, corresponding well to the tripole SST anomalies between the positive and negative EHT phases.It indicated that the cooling (warming) North Atlantic SST releases less (more) energy to force the overlying atmosphere at the decadal scale, rather than atmosphere controlled the ocean through anomalous turbulent heat fluxes.From observed changes of the indies, we can also see that the decadal variation of the summer NAO also lags behind that of NATI (Fig. 6b).Thus, on the decadal scale, SST variation may play a more important role in regulating the overlying atmosphere.
To further identify the driving effects of North Atlantic tripole SST on the wave train and decadal changes in EHT in MHLCA, two simulations using the CAM5 were performed.The configure set "F_2000_CAM5" was selected for them.It meant that the atmospheric composition was constant in the year 2000.The resolution was 2.5° longitude × 1.9° latitude, with 26 vertical hybrid levels.The first simulation was a control run with the model's modern climatological SST and sea ice (hurrell_sst_ifrac.1 × 1.050606.nc,which is the default setting of SST and sea ice cover climatological data used by CAM to simulate present climate).The second simulation was a sensitivity run, an idealised SST was constructed by imposing the observed differences in As shown in Fig. 13, the tripole-like SAT anomalies were evident in the North Atlantic, which well responded to the SST boundary forcing.At the same time, significant surface warming can be found in western-central Siberia.Similar with processes in the observations and coupled model, an anomalous anticyclonic-cyclonic teleconnection wave train was simulated from the North Atlantic to Siberia.Thus, the local anticyclonic conditions led to the reduced total cloud cover and thereby the surface warming in western-central Siberia (figure not shown).It indicated that the tripole SST anomalies in the North Atlantic can cause lead to decadal change in EHT in MHLCA by causing anomalous atmospheric meridional temperature gradient over the North Atlantic and related large-scale circulation anomalies over Eurasia.In addition, we can find that the cyclone and anticyclone centers of the wave train simulated by CAM5 were weaker and slightly displaced from the observed and CESM simulated positions.It implied that the air-sea coupling process may also play an important role in strengthening the excitation and downstream propagation of wave trains.

Potential precipitation-related feedback
In the past two decades, more compound hot and dry extremes have been reported worldwide (e.g., AghaKouchak et al. 2014;Miralles et al. 2014;Wu et al. 2021).Enhanced land-atmosphere coupling, associated with persistent soil moisture deficits, appeared to intensify surface warming and anticyclonic circulation anomalies, and thereby leading to stronger EHT events (Zhang et al. 2020b).In this positive feedback process, a reduction in precipitation was also one of the important processes (Miralles et al. 2019).
Focusing on the present study, in both the observation and simulation, decreased precipitation was evident over the central Siberian region, whereas increased precipitation was apparent over the south of Lake Baikal (Fig. 14).In the Siberia, the decreased precipitation was likely caused by the decreased SST in the tropical and high-latitude Atlantic (Sun et al. 2015;Nicolì et al. 2020).At the same time, this precipitation anomaly pattern was partially consistent with anomalous EHT events, suggesting that precipitation-related feedback was possibly involved in the decadal variability of EHT events in MHLCA.

Conclusion and discussion
In this study, the observed spatiotemporal characteristics of EHT events in MHLCA (the mid and high latitudes of continental Asia) were investigated.After removing the linear trend, significant decadal variation of EHT events could be found in the observation.During the periods 1980-1994 and 2012-2019, EHT events were more frequent and stronger in western-central Siberia than normal, ) between positive and negative phases.Areas that exceed the 95% confidence level are dotted whereas they were less frequent and weaker in the wide area south of Lake Baikal.The pattern of anomalous EHT events was found to show a meridional mode, and during the period 1995-2011 this EHT event pattern reversed.
Further analysis indicated that the decadal North Atlantic tripole SST anomaly pattern and associated sea-air interaction played important roles in shaping the abovementioned decadal variability of EHT events in MHLCA.During a positive NATI phase (i.e., negative-positive-negative SST anomaly phase), stronger eastward propagation of WAF (wave activity flux) can cause a significant anticyclonic-cyclonic teleconnection wave train in the troposphere over the mid and high latitudes of the Eurasian continent.Anomalous anticyclonic circulation was evident over western-central Siberia, which led to an anomalous high-pressure center and reduced total cloud cover; and because of stronger solar radiation heating, EHT events became more frequent and stronger in this region.In addition, anomalous anticyclonic circulation also caused an anomalous cyclonic condition over the wide area to the south of Lake Baikal, which reduced the frequency and weakened the severity of EHT events there.Results based on a long-term CESM simulation confirmed the important role of North Atlantic air-sea interaction in abovementioned processes.Further analysis based on the observed turbulent heat fluxes and CAM5 simulations suggested that, at the decadal scale, the tripole SST was likely to play a more important driving role.
It was noted that there were also significant positive SST anomalies in the subtropical Pacific between positive and negative EHT phases (Fig. 6a).Analysis based on observations cannot infer whether the positive SST anomalies contribute to the warming in western-central Siberia.It was interesting that CAM simulations produce similar warming in the eastern subtropical Pacific (Fig. 13a).It implied that the North Atlantic tripole SST anomalies might be important factor leading to positive SST anomalies in the subtropical Pacific through atmospheric bridges, which needs further investigation.
Recently, Hong et al. ( 2022) also investigated variations of EHT events in MHLCA, with a particular focus on interannual-scale physical processes.The Polar-Eurasian teleconnection pattern, NAO and Pacific-Japan/East Asia-Pacific pattern were all suggested to be the atmospheric factors in regulating the EHT events in MHLCA on the interannual scale.The related local physical processes were consistent with the present study.At the same time, they further pointed out the contribution of Atlantic multidecadal oscillation (AMO) and global warming to the longterm trend of EHT events in the region.In this study, however, we focused more on the impacts of decadal variability on extreme events in MHLCA, and indicated the roles of the decadal North Atlantic tripole SST anomaly and its coupled atmospheric circulation.
In this study, the analyses based on observations were detrended, and the external forcing factors such as greenhouse gases in the long-term simulations were kept constant.This means that the study highlighted the impact of decadal-scale internal variability in the climate system on MHLCA extreme weather events.Numerous studies have projected EHT events/heatwaves to become more frequent and stronger in the future (e.g., Zhou et al. 2016;Bador et al. 2017;Chen and Dong 2021;Rasmijn et al. 2018;Yin et al. 2020).Based on an early model result, Meehl and Tebaldi (2004) suggested that the European and North American heat-wave-related atmospheric circulation pattern will be intensified by ongoing increases in greenhouse gas concentrations, which in turn will lead to more severe heat waves.
In the MHLCA region, warming is more pronounced due to the Arctic amplification effect.However, it remains unclear whether the coupling process of the North Atlantic at the decadal scale will change under a more intense warming background in the future, and whether the large-scale circulation anomalies generated from the North Atlantic coupled processes will be affected by global warming.Firstly, the increase in global average temperature can significantly contribute to the enhancement of EHT events (e.g., Dong et al. 2017).Against the background of future warming, it is necessary to further investigate whether the circulation system related to decadal-scale internal variability of the North Atlantic will strengthen and further increase the intensity and frequency of EHT using the multimodel results from the latest two phases (5 and 6) of CMIP.tific and Technological Innovation Capacity Improvement Project of Chengdu University of Information Technology (KYQN202202).

Data Availability
The daily maximum temperature is openly available from the Climate Prediction Center Global Temperature dataset at https:// psl.noaa.gov/ data/ gridd ed/ data.cpc.globa ltemp.html.The monthly precipitation is openly available from the Climatic Research Unit Time-Series version 4.04 (Harris et al. 2020).The observed SST is openly available from the Hadley Center Monthly Sea Surface Temperature dataset at https:// www.metoffi ce.gov.uk/ hadobs/ hadis st/ (Rayner et al. 2003).Monthly ERA-Interim data is openly available from the European Centre for Medium-Range Weather Forecasts at https:// apps.ecmwf.int/ datas ets/ data/ inter im-mdfa/ levty pe= sfc/ (Dee et al. 2011).The output of CESM (version 1.2.2) long-term simulation and CAM5 simulations are too large to be publicly archived with available resources.Their detailed experimental design has been listed in the Sect. 2 and the simulation results section, respectively.They are also available from the corresponding author on reasonable request.

Fig. 1 a
Fig. 1 a First leaing EOF mode of anomalous summer EHT frequency (unit: days).b Normalized PC time series corresponding to the EOF mode of summer EHT frequency during 1980-2019 (bars)

Fig. 7 aFig. 8
Fig. 7 a First leading EOF mode of observed summer North Atlantic SST anomalies and b its PC during 1980-2019 (bars) along with the nine-year running mean (black line).In a black boxes are used to define the summer NATI, representing the subpolar area (region 1 : 45°-60° N, 50°-15° W), mid-latitude North Atlantic area (region 2 :

Fig. 9
Fig. 9 Regression maps of summer a total cloud cover (unit: %), b SAT (unit: °C), c EHT frequency (unit: days) and d EHT intensity (unit: °C) anomalies upon the summer NATI.The results are based

Fig. 10 a
Fig. 10 a First leading EOF mode of CESM simulated summer North Atlantic SST anomalies and b its PC (color bars) and corresponding 9-year running mean (black line).c Simulated summer NATI and (color bars) and its 9-year running mean (black line)

Fig. 12
Fig. 12 Differences in CESM simulated summer a total cloud cover (unit: %) and b SAT (unit: °C) between positive and negative NATI phases.Areas that exceed the 95% confidence level are dotted