Abrupt summer warming and changes in temperature extremes over Northeast Asia since the mid-1990s: Drivers and physical processes

This study investigated the drivers and physical processes for the abrupt decadal summer surface warming and increases in hot temperature extremes that occurred over Northeast Asia in the mid-1990s. Observations indicate an abrupt increase in summer mean surface air temperature (SAT) over Northeast Asia since the mid-1990s. Accompanying this abrupt surface warming, significant changes in some temperature extremes, characterized by increases in summer mean daily maximum temperature (Tmax), daily minimum temperature (Tmin), annual hottest day temperature (TXx), and annual warmest night temperature (TNx) were observed. There were also increases in the frequency of summer days (SU) and tropical nights (TR). Atmospheric general circulation model experiments forced by changes in sea surface temperature (SST)/sea ice extent (SIE), anthropogenic greenhouse gas (GHG) concentrations, and anthropogenic aerosol (AA) forcing, relative to the period 1964–93, reproduced the general patterns of observed summer mean SAT changes and associated changes in temperature extremes, although the abrupt decrease in precipitation since the mid-1990s was not simulated. Additional model experiments with different forcings indicated that changes in SST/SIE explained 76% of the area-averaged summer mean surface warming signal over Northeast Asia, while the direct impact of changes in GHG and AA explained the remaining 24% of the surface warming signal. Analysis of physical processes indicated that the direct impact of the changes in AA (through aerosol–radiation and aerosol–cloud interactions), mainly related to the reduction of AA precursor emissions over Europe, played a dominant role in the increase in TXx and a similarly important role as SST/SIE changes in the increase in the frequency of SU over Northeast Asia via AA-induced coupled atmosphere–land surface and cloud feedbacks, rather than through a direct impact of AA changes on cloud condensation nuclei. The modelling results also imply that the abrupt summer surface warming and increases in hot temperature extremes over Northeast Asia since the mid-1990s will probably sustain in the next few decades as GHG concentrations continue to increase and AA precursor emissions over both North America and Europe continue to decrease.


Introduction
During the 20th century the global averaged surface air temperature (SAT) showed an unequivocal warming due to anthropogenic as well as natural causes (e.g., Kosaka and Xie, 2013;Trenberth et al., 2014;Steinman et al., 2015).Distinguishing between forced and natural variability and understanding the drivers are an essential step in developing ond half of the 20th century, the EASM circulation underwent a weakening tendency, which was characterized by a decreasing lower-tropospheric southwesterly wind and an increasing sea level pressure (SLP) over East Asia.Associated with the weakened EASM circulation were increases in precipitation over the middle and lower reaches of the Yangtze River valley and decreases over northern China, exhibiting a so-called "southern-flood-northern-drought" (SFND) pattern (e.g., Ding et al., 2008Ding et al., , 2009;;Zhao et al., 2010;Huang et al., 2013;Wang et al., 2013, Zhao et al., 2014, Han et al., 2015;Ueda et al., 2015).Also accompanying the weakening EASM circulation were changes in seasonal mean SAT and temperature extremes over China and a large domain of Northeast Asia, and an increase in the frequency of dry spells and droughts (Ren et al., 2010;Wei and Chen, 2011;Li et al., 2012, Qi andWang, 2012;Wang et al., 2012;Wen et al., 2013;Chen and Lu, 2014;Gao et al., 2014aGao et al., , 2014b;;Qian et al., 2014;Song et al., 2014;Sun et al., 2014;Yang et al., 2014).
Previous attribution studies suggest that the changes in SSTs, GHG concentrations, and AA precursor emissions through aerosol-radiation and aerosol-cloud interactions are all important factors for the SFND pattern and associated drought over North and northeastern China (e.g., Kwon et al., 2007;Feng and Hu, 2008;Kosaka and Xie, 2013;Andrews, 2014;Gao et al., 2014aGao et al., , 2014b;;Kamae et al., 2014aKamae et al., , 2014b;;Qian et al., 2014;Song et al., 2014;Trenberth et al., 2014;Urabe and Maeda, 2014;Folini and Wild, 2015;Ueda et al., 2015;Wilcox et al., 2015;Zhang and Zhou, 2015).Previous studies have also revealed that decadal climate variability over East Asia is associated with low-frequency SST variability over the Pacific Ocean (e.g., Kwon et al., 2007, Kosaka andXie, 2013;Kamae et al., 2014a;Trenberth et al., 2014;Ueda et al., 2015) and Atlantic Ocean (Feng and Hu, 2008;Kamae et al., 2014a;Qian et al., 2014) via atmospheric teleconnections.Using AGCM simulations, Andrews (2014) found that changes in SST/SIE contributed about 60% of the summer land warming trends during 1979-2008 in the Northern Hemisphere, while the direct radiative impact of changes in anthropogenic forcings contributed about 40% and was an important factor for reproducing the recent Northern Hemisphere continental warming.Kamae et al. (2014a) demonstrated that the direct warming effect due to anthropogenic forcing and indirect warming effect via increased SST are important for the recent land surface warming and increased frequency of hot summers over the midlatitude Northern Hemisphere.
East Asia is under both local and non-local influences of AA (e.g., Cowan and Cai, 2011;Dong et al., 2016a).The study of Cowan and Cai (2011) indicated that Asian aerosols lead to a weakening of the EASM through aerosol-radiation and aerosol-cloud interactions (e.g., Twomey, 1977;Rosenfeld et al., 2008).The addition of non-Asian aerosols leads to further suppression of monsoon rainfall over Asia, highlighting the importance of the non-Asian aerosols in changes in Asian summer monsoon circulation and rainfall across Asia.Dong et al. (2016a) indicated that aerosol emissions over Asia and Europe both influence the thermal contrast between Asian land and the western North Pacific, leading to a weakening EASM, but through very different processes.Due to the implementation of air quality legislation, AA precursor emissions over both Europe and North America have continuously declined since the 1980s (Smith et al., 2011;Kühn et al., 2014).In contrast, over the same period, aerosol emissions in China and India have continuously increased (Lu et al., 2011).
The SFND pattern and warming over Northeast Asia have emerged particularly strongly since the mid-1990s (Kwon et al., 2007;Chen and Lu, 2014;Han et al., 2015;Ueda et al., 2015).This is roughly coincident with rapid increases in aerosol emissions from China and India and decreases in aerosol emissions from North America and Europe.Although previous studies have demonstrated the effects of recent decadal SST variation on decadal change in summer precipitation and surface warming over East Asia (e.g., Feng and Hu, 2008;Han et al., 2015;Ueda et al., 2015), a comprehensive and quantitative understanding of the drivers and physical processes governing the SFND pattern, and the associated changes in temperature extremes over Northeast Asia, is still lacking.
The main aims in carrying out this study were to understand what has caused the observed rapid changes in SAT and increases in hot temperature extremes over Northeast Asia since the mid-1990s.In particular, we sought to quantify the relative roles of changes in (1) SST/SIE, (2) anthropogenic GHG forcing through its direct impact, and (3) AA forcing through aerosol-radiation and aerosol-cloud interactions, in shaping these regional changes in surface warming and temperature extremes, by performing a set of numerical experiments using the atmospheric component of HadGEM3.It is important to stress that we did not address the anthropogenic contribution of changes in GHG concentrations and AA precursor emissions to observed SST/SIE changes; rather, we considered the SST/SIE changes as an independent forcing factor.
The structure of the paper is as follows: Section 2 describes the observed changes since the mid-1990s over Northeast Asia.In section 3, the model experiments and model climatology are described briefly.Section 4 presents the simulated changes in response to different forcings.Section 5 elucidates the physical processes involved in the responses.Conclusions are given in section 6.

Observed changes over Northeast Asia since
the mid-1990s

Observed changes
Following the work of Chen and Lu (2014), a number of temperature and precipitation indices averaged in a rectangular region (40 • -50 • N, 90 • -130 • E) over Northeast Asia during the summer (June-September, JJAS) or annual mean were calculated.Figure 1 shows the temporal evolution of indices of SAT, precipitation, and a number of temperature extreme anomalies, relative to the climatology, averaged over the whole time series.These time series clearly exhibit decadal to multidecadal variability, with cooling and wet conditions, from the 1960s to the 1980s, which might have been associated with North Atlantic cooling (e.g., Liu and Chiang, 2012).One of the most important features revealed is the abrupt surface warming since the mid-1990s (Fig. 1a).This is in line with the findings of Chen and Lu (2014), who demonstrated that the abrupt point for the decadal shift was 1993/94.We therefore compared the changes during the recent 16 years (1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011) relative to the early period of 1964-93 and found that the area-averaged SAT changes over Northeast Asia ranged from 1.21 • C to 1.27 • C, from four datasets.Accompanying the abrupt increase in the summer mean SAT was a decrease in local precipitation (Fig. 1b), and the area mean decreases ranged from −0.22 to −0.18 mm d −1 , about 10% of the climatological precipitation (e.g., Chen and Lu, 2014;Wilcox et al., 2015).
The rapid summer mean surface warming and drying over Northeast Asia since the mid-1990s accompanied rapid changes in temperature extremes (Figs.1c-e).The changes in these temperature extremes indicate a rapid increase in summer mean Tmax (daily maximum temperature), Tmin (daily minimum temperature), TXx (annual hottest day temperature), and TNx (annual warmest night temperature).However, the change in the diurnal temperature range (DTR) was relatively small (Fig. 1c). Figure 1e also shows a rapid increase (by 14 and 6 days, respectively) in the frequency of summer days (SU; annual number of days when Tmax > 25 • C) and tropical nights [TR; annual number of days when Tmin > 20 • C; see Donat et al. (2013)] over Northeast Asia.
The spatial patterns of changes in SAT and precipitation between the recent period and the early period are illustrated in Fig. 2. The SAT change shows a warming everywhere over the Eurasian continent (e.g., Tang et al., 2012), with the largest warming over Northeast Asia.Changes in precipitation are characterized by a decrease of 0.2-0.8mm d −1 over Northeast Asia, and an increase of 0.2-1.6 mm d −1 over South China and the lower reaches of the Yangtze River valley; and this is the SFND pattern-related to the weakening of the EASM-revealed in many previous studies (e.g., Ding et al., 2008Ding et al., , 2009;;Zhao et al., 2010;Wang et al., 2013;Chen and Lu, 2014).
The spatial patterns of changes in some temperature ex- tremes during the recent period relative to the early period are illustrated in Fig. 3.The main features of change are the increase in both Tmax and Tmin since the mid-1990s in most regions over the Eurasian continent.However, the changes in Tmax and Tmin show large spatial variations (Figs.3a and  b).Tmax changes show a large increase over Eastern Europe and Northeast Asia, and a relatively small increase over East China.Meanwhile, Tmin changes show a relatively uniform increase in a large domain over Asia.As a result, the changes in DTR indicate very large spatial variations, with an increase over Eastern Europe, small changes over Northeast Asia (e.g., Shen et al., 2014), and a decrease over East China.The large decrease in DTR over the Tibetan Plateau is not consistent with the changes in Chinese station data (e.g., You et al., 2015), and this indicates large uncertainty in the temperatureextreme changes in this data-sparse region.Changes in TXx and TNx (Figs. 3d and e) between the two periods show a localized maximum (about 1.8 • C) over Northeast Asia.The changes in frequency of SU and TR also show an increase, with a local maximum of about 12-15 days and 2-8 days over Northeast Asia (Figs. 3f and g).What has caused these observed rapid changes in summer mean SAT and increases in temperature extremes over Northeast Asia?Understanding the drivers of these changes is key for predicting climate and extreme events in the next few decades for China and Northeast Asia.By using the atmospheric component of HadGEM3 in a series of model experiments with different forcings, we investigated the drivers and physical processes for these observed changes over Northeast Asia since the mid-1990s.The findings are reported in the following sections.

Model and model experiments
Since the mid-1990s, there have been changes in the potential drivers of climate over Northeast Asia.Global-mean SST has warmed relative to the early period of 1964-93 (Fig. 4a).SST anomalies are characterized by regional patterns associated with a positive phase of the AMO (Atlantic Multidecadal Oscillation) and a negative phase of the PDO (Pacific Decadal Oscillation) (e.g., Kosaka and Xie, 2013;Kamae et al., 2014a;Steinman et al., 2015).There have been significant increases in anthropogenic GHG concentrations (11% increase in CO 2 , 18% increase in CH 4 , and 6% increase in N 2 O) and significant changes in AA precursor emissions.The changes in annual mean sulfur dioxide emissions-the most important aerosol precursor-show decreases over Europe and North America and increases over Asia (Fig. 4b).
A set of climate model experiments were run to assess the extent that the observed rapid increases in SAT and changes in temperature extremes over Northeast Asia since the mid-1990s can be reproduced when the model is forced by changes in SST/SIE, GHG and AA forcings, and to quantify the relative roles of individual forcing factors.The model used was HadGEM3-A (Hewitt et al., 2011), with a horizontal resolution of 1.875 • longitude by 1.25 • latitude and 85 vertical levels.The model includes an interactive tropospheric chemistry scheme with eight aerosol species and considers both aerosol-radiation and aerosol-cloud interactions.A detailed description of the tropospheric chemistry scheme was documented in Bellouin et al. (2013).The experiments, summarized in Table 1, were carried out and documented in detail in Dong et al. (2016b).The recent recovery of Sahel rainfall amounts was investigated in Dong and Sutton (2015), while the rapid summer surface warming and changes in temperature extremes since the mid-1990s over Western Europe were investigated in Dong et al. (2016b), using the same set of experiments.The response to a particular forcing was diagnosed by the mean difference between a pair of model experiments (one that included and one that excluded that forcing), based on the last 25 years of each experiment.The response to changes in all forcings was estimated by the mean difference between experiments ALL and CON-TROL, and the response to changes in SST/SIE by the difference between SSTONLY and CONTROL.The response to changes in AA forcing was estimated by the mean difference  between experiments ALL and SSTGHG, while the response to changes in GHG by the difference between SSTGHG and SSTONLY, with the assumption that the responses to individual forcings combined linearly.Seasonal mean changes for summer (JJAS) were the main focus of study.The statistical significance of summer mean changes and the 90% confidence intervals of the some area-averaged changes over Northeast Asia, in both observations and model experiments, in response to different forcings, were assessed using the twotailed Student's t-test.

Model-simulated climatology
This section compares the model-simulated climatological features of CONTROL with observed features (Fig. 5).Figures 5a and b show the SAT based on NCEP-NCAR reanalysis data (Kalnay et al., 1996) and observed precipitation based on GPCP (Adler et al., 2003), whereas Figs.5c and d show the corresponding model-simulated SAT and precipitation.
In observations, the SAT pattern was characterized by strong meridional temperature gradients around 45 • N over East Asia (Fig. 5a).The model reproduced the large-scale pattern of SAT well (Fig. 5c).Observed precipitation showed local maxima over the Western Ghats, Indochina, and the western North Pacific (Fig. 5b).The model simulation showed relatively weak precipitation over the Western Ghats, central India, and a large area over the western Bay of Bengal (Fig. 5d), in comparison with observations.However, the model overestimated precipitation over the western equatorial Indian Ocean and the Himalayan foothills (Fig. 5d).Convection tends to develop preferentially over these regions due to the orographic impact as the low-level southwesterly flow hits the Himalayan foothills, and the large moisture and heat availability over the equatorial Indian Ocean.The overestimation of precipitation over the equatorial Indian Ocean appears to be an inherent feature of the Met Office Unified Model (Martin et al., 2010), and is also present in other versions of the Met Office Hadley Centre Global Environmental Model (Levine and Turner, 2012;Guo et al., 2013, Dong et al., 2016a).Despite these deficiencies, the simulated SAT and precipitation over Northeast Asia compared fairly well with observations, giving credibility for using the model to inves-tigate the response of regional climate over Northeast Asia to changes in different forcings.

Model-simulated changes in response to different forcings
The spatial distributions of changes in summer mean SAT and precipitation in response to different forcings, together and individually, are given in Fig. 6.The model responses to changes in SST/SIE and anthropogenic forcings (Fig. 6a), relative to CONTROL, reproduced the general patterns of observed changes in summer mean SAT over the Eurasian continent, but the simulated warming was more uniform than observed (Fig. 2a).The response in the model of surface warming to changes in all forcings was predominantly due to SST/SIE changes, which were responsible for the general warming found over the low-to-middle latitude Eurasian continent, with a magnitude of 0.4 • C-1.2 • C (Fig. 6c), in which the recent warming of the North Atlantic might have played an important role (e.g., Feng and Hu, 2008;Steinman et al.,  2015).The changes in SAT due to aerosol forcing manifested as a band (40 • -55 • N) of positive anomalies extending from Europe to Northeast Asia, with a magnitude of 0.2 • C-0.8 • C (Fig. 6g); whereas, the direct impact of changes in GHG concentrations induced more localized warming over Europe, with weak warming over Northeast Asia (Fig. 6e).These results indicate that the direct impact of changes in anthropogenic GHG and AA changes plays an important role in midlatitude warming over the Eurasian continent, which is consistent with the recent studies of Kamae et al. (2014a) and Andrews (2014).
In response to changes in all forcings, the model simulated an increase in precipitation over Yangtze River basin, but did not reproduce the drying over Northeast Asia (Fig. 6b).Sperber et al. (2013) showed that current climate models still have biases in simulating Asian summer monsoon circulation and associated precipitation.The results presented in this paper, and those of Christidis et al. (2013) and of Wilcox et al. (2015), suggest that these biases need to be reduced in order for robust attribution of decadal change in precipitation over Asia to be obtained.Changes in GHG concentrations lead to wet conditions over East Asia and Northeast Asia (e.g., Song et al., 2014), while changes in AA precursor emissions lead to opposite changes (e.g., Wilcox et al., 2015).
Associated with a band of surface warming that resulted from AA changes was a band of reduced precipitation extending from Europe to Northeast Asia (Figs. 6g and h), implying a role of changes in non-local aerosol emissions for SAT and precipitation over Northeast Asia.This will be discussed in more detail in section 5.
The responses of seasonal mean Tmax, Tmin and DTR to changes in different forcings, together and individually, are illustrated in Fig. 7.The model-simulated changes in response to all forcings showed an increase in both Tmax and Tmin over the Eurasian continent, in which the changes in Tmax were large over Eastern Europe but the Tmin changes showed a more or less uniform increase over the midlatitudes (Figs.7a and b).As a result, DTR changes showed an increase over Eastern Europe and a weak change over Northeast Asia (Fig. 7c), similar to the observed change (Fig. 3c), albeit the enhancement in DTR over Eastern Europe was overestimated in the model responses to changes in all forcings.Changes in SST/SIE led to an increase in both Tmax and Tmin over the Eurasian continent, with the increase in Tmin being larger than the increase in Tmax over Northeast Asia, and this led to a decrease in local DTR (Figs. 7d-f).The direct impacts of changes in GHG and AA forcings both led to more localized changes in Tmax and Tmin, with the GHG changes   (d-f) response to changes in SST/SIE (SSTONLY minus CONTROL); (g-i) response to the changes in GHG concentrations (SSTGHG minus SSTONLY); and (j-l) response to changes in AA precursor emissions (ALL minus SSTGHG).Thick black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.
inducing a weak cooling in Tmax, a weak warming in Tmin, and, therefore, a decrease in DTR, over Northeast Asia (Figs. 7g-i).The decrease in Tmax was mainly due to the increase in cloud cover locally, and this will be discussed in more detail in section 5.2.In contrast, the direct impacts of changes in aerosols led to an increase in Tmax over regions extending from Eastern Europe to Northeast Asia, while aerosols induced little change in Tmin.As a result, the changes in aerosols induced a significant increase in DTR over Eastern Europe and Northeast Asia.Therefore, weak changes in DTR over Northeast Asia in response to changes in all forcings is due to the compensating impacts from changes in GHG and SST/SIE (reducing DTR) and the impacts from changes in aerosols (enhancing DTR).
The spatial distributions of the responses in the frequencies of SU and TR are illustrated in Fig. 8.In response to changes in all forcings (Figs.8a and b), the model reproduced both the pattern and magnitude of the frequency increases in SU and TR shown in observations over Northeast Asia (Figs. 3f and g). Figure 8 also indicates that changes in SST/SIE played an equally important role as changes in aerosols in the increase in the frequency of SU (Figs. 8c, e and g), while the increase in the frequency of TR was predominantly due to changes in SST/SIE (Figs. 8d, f and h).Consistent with increased local aerosol emissions, the frequency of SU and TR decreased over South, Southeast and East Asia (Figs. 8g and  h).
A number of area-averaged changes in summer seasonal mean SAT and seasonal or annual changes in temperature extremes over Northeast Asia, for both observations and modelsimulated responses, are illustrated in Figs. 9 and 10.Quantitatively, the model-simulated changes in response to changes in all forcings reproduced changes in summer seasonal areaaveraged SAT, Tmax and Tmin realistically, and this led to very small changes in the model-simulated DTR, similar to the observed change (Fig. 9a).The area-averaged increases in the frequency of SU and TR in response to changes in all forcings were also similar to observed changes, albeit weaker (Fig. 10a).In response to changes in all forcings, the model produced an area-averaged summer mean SAT change of 1.07 • C over Northeast Asia, which was very close to the observed changes of 1.21 • C to 1.27 • C. Furthermore, most of the surface warming signal (76% change in SAT) was simulated in response to changes in SST/SIE.If the responses to individual forcings combined linearly, this implies that 21% of the SAT change was the response to changes in AA precursor emissions, and an additional 3% surface warming was the response to changes in GHG concentrations (Fig. 9b).Both changes in SST/SIE and AA led to an increase in Tmax, while the Tmin increase was predominantly due to changes in SST/SIE (Fig. 9b).The model results further suggest that the direct impact of changes in AA precursor emissions played a dominant role in the increase in TXx (Fig. 9b), and an equally important role as SST/SIE changes in the increase in the frequency of SU (Fig. 10b).Meanwhile, the increases in TNx and the frequency of TR were mainly due to the recent warm-ing of the ocean (Figs.9b and 10b).
However, Fig. 9a indicates that the model response to changes in all forcings underestimated the changes in TXx and TNx.This underestimation might have been related to the weak land surface feedback in the model simulation, where the model response showed a weak change in summer mean precipitation (Fig. 6b) and therefore a small change in soil moisture (figure not shown), while observations indicate there was a drying trend over Northeast Asia (Fig. 2b).Previous work (e.g., Seneviratne et al., 2014) points out the role of land surface feedback associated with drying on the change in hot temperature extremes, but this land surface feedback over Northeast Asia in response to changes in all forcings was weak due to the compensating impacts that resulted from the changes in AA, as well as those from changes in SST/SIE and GHGs (Figs. 6d, f and h).
The significant role of changes in AA precursor emissions for changes in both summer mean SAT and temperature extremes over Northeast Asia is interesting.Given the fact that local AA precursor emissions over Asia are increasing and European and North American emissions are decreasing, which regional change of AA precursor emissions was responsible for those changes seen in the model experiment over Northeast Asia?The summer mean SAT changes in re-  sponse to changes in AA precursor emissions shown in Fig. 6 indicate a band of positive SAT anomalies and reduced precipitation extending from Europe to Northeast Asia in the midlatitudes.Changes in AA precursor emissions also induced a band of increase in the seasonal mean Tmax, DTR (Figs. 7j and l), TXx (figure not shown), and the frequency of SU (Fig. 8g).These results imply that changes over Northeast Asia in response to changes in AA precursor emissions are the downstream development of the European response to reduced emissions over Europe.They demonstrate a role of non-local AA changes for the summer mean climate over the Asian continent (e.g., Dong et al., 2016a) and associated changes in temperature extremes.The detailed physical processes responsible for the responses to changes in different forcings are discussed in next section.
5. Physical processes responsible for the model-simulated changes induced by different forcings

Response induced by changes in SST/SIE
The spatial distributions of the summer mean responses of some key variables induced by changes in SST/SIE are shown in Fig. 11, while the area-averaged changes over Northeast Asia in surface radiation, surface turbulent heat fluxes and some other key variables in response to changes in different forcings are given in Fig. 12.In response to SST/SIE changes, the most significant changes were the increased water vapor in the atmosphere over the ocean and Eurasian continent (Fig. 11a).The increased water vapor led to an increase in clear-sky net surface longwave (LW) radiation (Fig. 11b), which tended to warm the surface.However, Fig. 11c indicates that the net surface LW change was small, and this small net LW anomaly (Fig. 11c) likely reflected a balance between increased emissions from the warmer surface and increased absorption of radiation emitted by the warmer and moister atmosphere.In addition, the reduction in medium-level cloud (Fig. 11d) related to the increased stability over land in the midlatitudes (e.g., Dong et al., 2009) induced positive anomalies in the shortwave (SW) cloud radiative effect (CRE) (Fig. 11e).The reduction in the upward surface sensible heat flux over Northeast Asia (Fig. 11f), related to surface air warming more than the surface skin temperature (Fig. 12b), also contributed to the surface warming.In summary, it was the increased water vapor in the atmosphere over the Eurasian continent in the midlatitudes in response to SST/SIE changes that induced the surface warming over Northeast Asia, with fur- ther positive SW feedbacks through reduced cloud cover and decreased upward sensible heat flux related to the increase in the difference between the surface air and skin temperatures.This role of reduced cloud cover on summer surface warming over Eastern Europe agrees with the studies of Tang and Leng (2012) and Tang et al. (2012), who highlighted the role of decreased cloud cover in the recent observed changes in European hot summers.Our study further suggests that changes in SST/SIE might be partially responsible for the reduced summer cloud cover over Europe during recent decades.

Response induced by changes in GHG forcing
The spatial distributions of the summer mean responses of a number of key variables induced by changes in GHG forcing are given in Fig. 13.The main characteristics of surface energy changes over land in the midlatitudes of 30 • -50 • N indicate an increase in the net surface LW radiation (Fig. 13a).This increase was partly due to the increase in downward surface clear LW radiation (Fig. 13b; as expected for an increase in the Greenhouse Effect) and partly due to positive LW feedback resulting from increased water vapor in the atmosphere (Fig. 13c).In summary, it was the increased Greenhouse Effect that induced surface warming over Northeast Asia, with positive feedbacks from increased water vapor and reduced upward sensible heat flux (Fig. 13d), in response to GHG changes.The reduction in the upward sensible heat flux was the result of surface air warming being more than the surface skin temperature (Fig. 12b).The changes in surface SW radiation were closely related to changes in the SW CRE (Figs. 13e and f) as a consequence of increased total cloud cover (Fig. 12b) and water vapor (Fig. 13c), which tended to cool the surface over Northeast Asia.

Response induced by changes in AA precursor emissions
The spatial distributions of the summer mean responses of a number of key variables induced by changes in AA precursor emissions are illustrated in Fig. 14.Changes in sulfate aerosol optical depth (AOD) (Fig. 14a) indicate a decrease over Europe and an increase over South, Southeast, East Asia, and part of Northeast Asia, closely related to sulfur dioxide emission changes (e.g., Fig. 4b).Note that changes in AA in the model experiments included changes in sulfate, soot, biomass, and also other species (e.g., Bellouin et al., 2013).However, changes in soot AOD and biomass AOD (Figs. 14h  and i) were much smaller than changes in sulfate AOD, suggesting that the responses to AA changes were predominantly due to changes in sulfate.
Local decreases in AOD over Europe led to significant regional surface warming over Western and Eastern Europe (Fig. 6g), through aerosol-radiation and aerosol-cloud interactions (e.g., Nabat et al., 2014;Dong et al., 2016aDong et al., , 2016b)).This in turn led to reduced precipitation (Fig. 6h) and a decrease in soil moisture over Europe (Fig. 14b).Associated with changes over Europe was the downstream extension of a significant decrease in water vapor in the atmosphere from Eastern Europe to Northeast Asia (Fig. 14c), resulting from decreased evaporation (Fig. 14g) related to the drying of soil (Fig. 14b).The reduction in column-integrated water vapor induced a decrease in atmospheric heating by SW absorption (e.g., Mitchell et al., 1987, Dong et al., 2016b).This led to an increase in surface net clear-sky SW radiation, and therefore surface net SW radiation (Fig. 14d).The decrease in columnintegrated water vapor also induced a reduction in the LW radiative cooling of the atmosphere to the surface, leading to a decrease in net surface LW (Fig. 12a).The reduction in cloud cover associated with the decreased water vapor in the atmosphere led to a positive change in the SW CRE and a negative change in the LW CRE (Figs. 14e and f).This gave rise to larger changes in net surface SW and LW radiations than clear-sky net surface SW and LW radiations (Fig. 12a).
The increased surface net SW radiation and reduced upward latent heat flux (Figs.14d and g) led to a warming in SAT over the midlatitudes of 30 • -50 • N.These modelling results suggest that decreases in AA precursor emissions over Europe might also be partially responsible for the reduced summer cloud cover during recent decades over Eastern Europe (e.g., Tang and Leng, 2012;Tang et al., 2012), while landatmosphere interactions downstream drive responses of temperature and temperature extremes over Northeast Asia (Dai et al., 1999;Seneviratne et al., 2010;Mueller and Seneviratne, 2012;Schubert et al., 2014).The reduction in soil moisture led to reduced evapotranspiration and therefore resulted in decreased upward latent heat flux and increased upward sensible flux (Figs.12a and 14g).The increased upward sensible heat flux in turn warmed the air, leading to an increase in evaporative demand, which tended to dry the soil further.This positive feedback between land surface and atmosphere led to a decrease in cloud cover and an increase in surface SW radiation, which again enhanced land surface drying.This positive feedback loop illustrates the importance of land surface processes in driving the regional responses of surface temperature and hot temperature extremes over Northeast Asia (Dai et al., 1999;Seneviratne et al., 2010; Mueller and Seneviratne, 2012;Schubert et al., 2014), and the resultant increases in Tmax, TXx and the frequency of SU, as well as the enhancement of DTR (Fig. 9b).In summary, the surface warming over Northeast Asia in response to changes in AA precursor emissions resulted from the increased net SW radiation due to the reduction in cloud cover and positive SW CRE related to tropospheric and surface drying associated with a coupled atmosphere-land response, rather than a direct impact of AA changes on cloud condensation nuclei and AOD.Also, the reduced upward latent heat flux tended to warm the surface and feed back positively to the reduction in cloud cover.The reduction in cloud cover in turn induced increases in Tmax, TXx and the frequency of SU.

Conclusions
This study investigated the drivers and physical processes involved in the abrupt summer mean surface warming over Northeast Asia and associated changes in temperature extremes since the mid-1990s.Firstly, we analyzed observations based on multiple data sets.Then, we carried out a set of numerical experiments using the atmospheric component of HadGEM3-A, to assess and quantify the relative roles of changes in three forcing factors [(1) SST/SIE; (2) anthropogenic GHG forcing through its direct impact; and (3) AA forcing through aerosol-radiation and aerosol-cloud interactions] in governing the changes in summer mean SAT and temperature extremes over Northeast Asia.No attempt was made to assess and quantify the anthropogenic contribution of changes in GHG concentrations and AA precursor emissions to the observed SST/SIE changes.The main findings and physical processes involved can be summarized as follows: (1) Observations indicate that there was an abrupt increase in summer mean SAT over Northeast Asia in the mid-1990s.The area-averaged SAT change over Northeast Asia in JJAS during the recent 16 years from 1996 to 2011, relative to the period 1964-93, ranged from 1.21 • C to 1.27 • C, from four data sets.Accompanying the abrupt increase in the summer mean SAT was a decrease in local precipitation.
(2) Also accompanying the abrupt increase in summer mean SAT were rapid changes in a number of temperature extremes.Specifically, there was a rapid increase in JJAS mean Tmax, Tmin, and in TXx and TNx over Northeast Asia, while the change in DTR was relatively small.
(3) Observational analysis also revealed a rapid increase (by 14 and 6 days, respectively) in the frequency of SU and TR over Northeast Asia since the mid-1990s.
(4) The model experiment forced by changes in SST/SIE and anthropogenic (GHG and AA) forcings together, relative to CONTROL, reproduced the general patterns of observed summer mean SAT changes and changes in temperature extremes over Northeast Asia.Changes in SST/SIE explained 76% of the area-averaged SAT signal over Northeast Asia, with 21% and 3% of the surface warming signal explained by the direct impacts of changes in AA precursor emissions and GHG concentrations, respectively.However, the precipitation responses to changes in all forcings were small, related to the compensating impacts from changes in SST/SIE and AA, thus precluding a definitive attribution regarding recent precipitation changes over Northeast Asia.
(5) The direct impact of the changes in AA precursor emissions, through aerosol-radiation and aerosol-cloud interactions, mainly related to the reduction in emissions over Europe, played a dominant role in the increase in TXx, and an equally important role as SST/SIE changes in the increase in the frequency of SU through the reduction in cloud cover, related to coupled atmosphere-land surface feedbacks, rather than through a direct impact of AA changes on cloud condensation nuclei and AOD.
(6) The relatively small change in DTR induced by changes in all forcings was due to the compensating impact from changes in GHG concentrations and SST/SIE (reducing the DTR related to increased water vapor in the atmosphere and the associated increase in cloud cover) and the impact from changes in aerosols (enhancing the DTR related to decreased water vapor in the atmosphere and the associated reduction in cloud cover).
Whilst each forcing factor led to summer mean surface warming and changes in temperature extremes over the Eurasian continent in the model simulations, the spatial patterns of the responses were different and possessed distinct features in each case.For example, changes in SST/SIE induced a more or less uniform summer mean surface warming, while changes in AA precursor emissions induced a band of warming and changes in temperature extremes in the midlatitudes extending from Europe downstream to Northeast Asia at the surface.The results of this study indicate that reduced cloud cover enhances surface warming in the midlatitudes over the Eurasian continent in response to changes in both AA precursor emissions and SST/SIE, through its impact on SW radiation at the surface.In response to changes in AA, the model results indicate a local positive feedback between surface warming and reduced cloud cover, precipitation, soil moisture, and evaporation.This mechanism is an important factor leading to local increases in summer mean Tmax, DTR, TXx, and the frequency of SU.The local change in cloud cover, in turn, is closely related to change in the hydrological cycle and land surface processes.These results strengthen the conclusion that a realistic representation of cloud and its radiative effect, as well as the hydrological cycle, in climate models, are fundamental to the credibility for using such models to predict responses to a change in anthropogenic climate forcing.In the near term (the next few decades), GHG concentrations will continue to rise and AA precursor emissions over North America and Europe will continue to decrease.The results of this study imply that the rapid summer surface warming and increases in hot temperature extremes over Northeast Asia since the mid-1990s will probably continue, or even enhance, in the near term.

Fig. 2 .
Fig. 2. Spatial patterns of differences in (a) SAT (units: • C) and (b) precipitation (units: mm d −1 ) between the recent period of 1996-2011 and the early period of 1964-93.The thick red boxes highlight Northeast Asia, where some area-averaged indices were calculated.

Fig. 3 .
Fig. 3. Spatial patterns of differences in some temperature extremes between the recent period of 1996-2011 and the early period of 1964-93: (a) Tmax; (b) Tmin; (c) DTR; (d) TXx; (e) TNx; (f) SU; and (g) TR.The units for SAT, Tmax, Tmin, DTR, TXx and TNx are • C. The units for SU and TR are d yr −1 .The thick red boxes highlight Northeast Asia.

Fig. 5 .
Fig. 5. Spatial patterns of the JJAS climatology for (a, c) SAT (units: • C) and (b, d) precipitation (units: mm d −1 ) in (a, b) observations and (c, d) CONTROL.The thick red boxes highlight Northeast Asia.

Fig. 6 .
Fig. 6.Spatial patterns of summer mean (JJAS) SAT (units: • C) (left panels) and precipitation (units: mm d −1 ) (right panels) changes: (a, b) response forced by changes in SST/SIE, GHG concentrations, and AA precursor emissions (ALL minus CONTROL); (c, d) response to changes in SST/SIE (SSTONLY minus CONTROL); (e, f) response to changes in GHG concentrations (SSTGHG minus SSTONLY); and (g, h) response to changes in AA precursor emissions (ALL minus SSTGHG).Thick black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.

Fig. 7 .
Fig. 7. Spatial patterns of summer mean responses for Tmax (left column), Tmin (middle column), and DTR (right column) (units: • C): (a-c) response forced by changes in SST/SIE, GHG concentrations, and AA precursor emissions (ALL minus CONTROL);(d-f) response to changes in SST/SIE (SSTONLY minus CONTROL); (g-i) response to the changes in GHG concentrations (SSTGHG minus SSTONLY); and (j-l) response to changes in AA precursor emissions (ALL minus SSTGHG).Thick black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.

Fig. 8 .
Fig. 8. Spatial patterns of responses for SU (left column) and TR (right) (units: d yr −1 ): (a, b) response forced by changes in SST/SIE, GHG concentrations, and AA precursor emissions (ALL minus CONTROL); (c, d) response to changes in SST/SIE (SSTONLY minus CONTROL); (e, f) response to changes in GHG concentrations (SSTGHG minus SSTONLY); and (g, h) response to changes in AA precursor emissions (ALL minus SSTGHG).Thick black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.

Fig. 9 .
Fig. 9. Observed and model-simulated summer seasonal mean (JJAS) changes for SAT, Tmax, Tmin, and DTR, and annual changes in TXx and TNx over Northeast Asia (units: • C): (a) observed changes (based on CRUTS3.2 and HadEX2) and simulated responses to changes in SST/SIE, GHG concentrations, and AA precursor emissions (ALL minus CONTROL) (colored bars indicate the central estimates and whiskers show the 90% confidence intervals, based on a two-tailed Student's t-test); (b) model-simulated changes in response to different forcings [SST & SIE is the response to changes in SST/SIE (SSTONLY minus CONTROL); GHG is the response to GHG concentrations (SSTGHG minus SSTONLY); and AA is the response to changes in AA precursor emissions (ALL minus SSTGHG)].

Fig. 10 .
Fig. 10.(a) Observed (HadEX2) and simulated annual changes over Northeast Asia in response to changes in SST/SIE, GHG concentrations, and AA precursor emissions (ALL minus CONTROL) for the frequency of summer days and tropical nights (units: d yr −1 ; colored bars indicate the central estimates and whiskers show the 90% confidence intervals, based on a two-tailed Student's t-test).(b) Model-simulated changes in response to different forcings [SST & SIE is the response to changes in SST/SIE (SSTONLY minus CONTROL); GHG is the response to GHG concentrations (SST-GHG minus SSTONLY); and AA is the response to changes in AA precursor emissions (ALL minus SSTGHG)].

Fig. 11 .
Fig. 11.Spatial patterns of responses to changes in SST/SIE (SSTONLY minus CONTROL): (a) column-integrated water vapor (units: kg m −2 ); (b) clear-sky surface LW; (c) surface LW; (d) medium-level cloud cover (units: %); (e) SW CRE; and (f) surface sensible heat flux.Radiation and fluxes are in W m −2 , and positive values mean downward.Radiation is the net component.Thin black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.

Fig. 12 .
Fig. 12. Model-simulated summer seasonal mean (JJAS) changes over Northeast Asia for (a) net surface radiation and flux changes (units: W m −2 ; positive values mean downward) and (b) soil moisture (SM; units: kg m −2 ), precipitation (Pr; units: 0.1 mm d −1 ), surface air and skin temperature difference (Ta-Ts; units: 0.1 • C), columnintegrated water vapor (WV; units: kg m −2 ), and cloud fraction (units: %) [SST & SIE is the response to changes in SST/SIE (SSTONLY minus CONTROL); GHG is the response to GHG concentrations (SSTGHG minus SSTONLY); and AA is the response to changes in AA precursor emissions (ALL minus SSTGHG)].The colored bars indicate the central estimates and whiskers show the 90% confidence intervals, based on a two-tailed Student's t-test.

Fig. 13 .
Fig. 13.Spatial patterns of responses to changes in GHG forcing (SSTGHG minus SSTONLY): (a) surface LW radiation; (b) clear-sky surface LW radiation; (c) column-integrated water vapor (units: kg m −2 ); (d) surface sensible heat flux; (e) surface SW radiation; and (f) SW CRE.Radiation and fluxes are in W m −2 , and positive values mean downward.Radiation is the net component.Thin black lines highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.

Fig. 14 .
Fig. 14.Spatial patterns of responses to changes in aerosols (ALL minus SSTGHG): (a) sulfate AOD at 0.55 um; (b) soil moisture (units: kg m −2 ), (c) column-integrated water vapor (units: kg m −2 ); (d) surface SW radiation; (e) total cloud cover (units: %); (f) SW CRE; (g) surface latent heat flux; (h) soot AOD at 0.55 um; and (i) biomass AOD at 0.55 um.Radiation and fluxes are in W m −2 , and positive values mean downward.Radiation is the net component.Thin black lines in (b-g) highlight regions where the changes are statistically significant at the 90% confidence level, based on a two-tailed Student's t-test.The thick red boxes highlight Northeast Asia.
SSTONLY Monthly climatological SST/SIE averaged over the period 1996-2011, with GHG concentrations and AA precursor emissions the same as in CON-TROL.