The effect of a well-resolved stratosphere on East Asian winter climate

This study performs a comparison between models with and without well-resolved stratosphere on the simulation of the East Asian winter climate. The results reveal that the Fifth Coupled Model Inter-comparison Project (CMIP5) models with model top above the stratopause have a better simulation of the distribution of surface air temperature, sea level pressure and precipitation than the models with a low-top below the stratopause. The discrepancy of the East Asian winter climate between high-top (HT) and low-top (LT) CMIP5 models is also evident in the future projection under higher (RCP85) and midrange (RCP45) emission scenarios. The HT models present about 1.3 and 1.7 °C higher surface air temperature in East Asian region under RCP45 and RCP 85 scenarios by the end of this century than that of the LT models. Therefore, the insufficient representation of the stratosphere may lead to underestimation of the anthropogenic global warming in regional scale and hence have the potential to lead to insufficient response action and mitigation measures.


Introduction
The stratosphere is conventionally believed to be passively influenced by the tropospheric circulation anomalies. This upward effect has been well documented in observational and theoretical studies of stratospheric phenomena such as stratospheric sudden warming (e.g., Holton and Mass 1976;Labitzke 1982;Matsuno 1971) and quasi-biennial oscillation (e.g., Baldwin et al. 2001;Holton and Lindzen 1972;Lindzen and Holton 1968;Plumb 1977). However, studies in the recent decades have revealed that stratosphere can play an important role in the near surface climate change. For example, data diagnosis showed that the stratosphere zonal-mean zonal wind anomalies can propagate poleward and downward to the troposphere consecutively from the stratopause region (Christiansen 2001;Kodera 1995;Kuroda and Kodera 1999). It was further illustrated that the anomalous Arctic Oscillation (AO) signal, which is the dominant mode of the geopotential fields north of 20°N, appears firstly in the upper stratosphere, and propagates all the way down to the troposphere and can have significant influence on the surface extended-range weather and shortterm climate (Baldwin andDunkerton, 1999, 2001;Baldwin et al., 2003b;Moritz et al., 2002). It's believed that "at its maximum, the magnitude and geographic scale of this influence may be comparable to that of El Niño-Southern Oscillation" (Baldwin et al. 2003a).
The influence from the stratosphere has been verified in numerical studies. For example, in a long-term (260 years pre-industrial control run) simulation performed with the CMCC-CMS climate model including a state-of-the-art Abstract This study performs a comparison between models with and without well-resolved stratosphere on the simulation of the East Asian winter climate. The results reveal that the Fifth Coupled Model Inter-comparison Project (CMIP5) models with model top above the stratopause have a better simulation of the distribution of surface air temperature, sea level pressure and precipitation than the models with a low-top below the stratopause. The discrepancy of the East Asian winter climate between high-top (HT) and low-top (LT) CMIP5 models is also evident in the future projection under higher (RCP85) and midrange (RCP45) emission scenarios. The HT models present about 1.3 and 1.7 °C higher surface air temperature in East Asian region under RCP45 and RCP 85 scenarios by the end of this century than that of the LT models. Therefore, the insufficient representation of the stratosphere may lead to underestimation of the anthropogenic global warming in regional scale and hence have the potential to lead to insufficient response action and mitigation measures.
stratosphere, Manzini et al. (2012) showed that the sea level pressure, surface temperature, sea ice cover and the Atlantic meridional overturning circulation anomalies in the Northern Hemispheric mid and high latitudes can be traced back to the long lasting stratospheric vortex anomalies in the inter-decadal time scales. Karpechko and Manzini (2012) found that the early winter (November-December) stratospheric changes in response to the increased CO 2 concentration could propagate downward and have the largest influence on the surface during mid and late winter (January-March). The association between the strength of stratospheric polar vortex (SPV) and East Asian winter monsoon was discovered in several studies (Wei and Bao 2012;Wei et al. 2015). The second leading mode of the surface air temperature in East Asia (EA), with opposite changes between the northern EA and the southern EA, is highly correlated with the strength of stratospheric polar vortex. And the diagnosis of multiple reanalysis datasets shows that the relationship between the SPV and the second EAWM (East Asian winter monsoon) mode increased since the late 1980s (Wei et al. 2015).
Climate models have become one of the primary tools in climate research as they help to improve our understanding on mechanism, prediction and adaption to climate change. Significant progress has been made in incorporating more realistic physical and chemical processes in the models. Meanwhile, current models are constructed on a finer spatial grid. However, the model progress in the stratosphere seems relatively slow. In the models which were used in the fourth climate change assessment report (AR4) of the Intergovernmental Panel on Climate Change (IPCC), half of the model have model top beneath 2 hPa (upper stratosphere around 48.5 km). These models show a significant cold bias in their upper levels (~10 hPa) compared to the reanalysis data, possibly due to low model top or insufficient stratospheric levels (Cordero and Forster 2006). In the most recent fifth assessment report (AR5) of the IPCC, 54 models were used for the assessment, among which only half (27) models have full stratosphere with model lid top above 1 hPa. These low-top models have weaker polar vortex variability, and reduced planetary wave activity in stratosphere, leading to much less frequency of major sudden stratospheric warming events in the stratosphere (Charlton-Perez et al. 2013;Osprey et al. 2013). The Climate-system Historical Forecast Project using 16 coupled models (8 high-top models and 8 low-top models) shows that the high-top models have a more realistic stratospheric response to El Nino and the QBO, and have a possible stratospheric pathway that enhance the wintertime prediction skill over high latitudes and North Atlantic region (Butler et al. 2016). The comparison study by Scaife et al. (2012) showed that the inclusion of full stratosphere in the Hadley Centre Global Environmental Model (HadGEM) and the ECHO-G Middle Atmosphere Model (EGMAM) could lead to the changes in regional wind and rainfall patterns, which are large enough to significantly alter the regional climate change projection. Hence, it is crucial to incorporate a fully resolved stratosphere in the models to get better near surface climate change simulations.
For the densely populated and economic booming East Asia region, climate variation of even moderate amplitude can have profound impact. There have been several evaluation studies with focus on the model reproduction of East Asian summer and winter monsoon climate (e.g., Boo et al. 2011;Gong et al. 2014;Wei and Bao 2012;Wei et al. 2014). It was shown that in the winter season, the current state-of-the-art CMIP5 models have the ability to reproduce circulation in EA region and progress had been made to relieve the surface cold bias, which was observed in the CMIP3 model ensemble. Meanwhile, CMIP5 models show more model consistence in most EAWM parameters (Wei et al. 2014, hereafter WEI14). However, it's still not clear whether the inclusion of a full stratosphere can make a better simulation of the climate in the East Asian region. And it's also necessary to reevaluate the future climate projection in this region, considering that with or without a full stratosphere may lead to significant climate change difference, which can be large enough to alter the regional climate change projection.
In this study, we examine the performance of the CMIP5 models with a well-resolved stratosphere in East Asian region, and compare with the models with a low model top. We also evaluate the effect of a well-resolved stratosphere on the regional climate change projection. As the stratosphere is mainly coupled with troposphere and have largest climate effect in winter season (Baldwin et al. 2003a, b), the study in this paper focuses mainly in the winter season. The layout of this paper is as follows. In Sect. 2, we will discuss the models and observation-based data used. In Sect. 3, CMIP5 EAWM climatology will be presented and comparison will be made between models with and without full stratosphere. Section 4 presents the discrepancy of the EAWM parameters/indices between high-top and lowtop CMIP5 Models. In Sect. 5, the regional climate change projection to the end of this century will be presented. We will present discussions and summaries in Sect. 6.

Data
The monthly mean data of National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) Global Reanalysis 1 (Kalnay et al. 1996;Kistler et al. 2001) and the NOAA's Precipitation Reconstruction Dataset (Chen et al. 2002) were used in this study, including surface air temperature (SAT), sea level pressure (SLP), zonal wind (u) and meridional wind (v) at 1000 hPa (UV1000), and zonal wind at 200 hPa (U200). Our analyses are focused on the boreal winter. Therefore, the wintertime mean was computed by averaging 3-month periods (December-February (DJF)) in each year. Finally, the climatological wintertime means were calculated by averaging each DJF over the period of model experiments. To compare with CMIP3 and CMIP5 model output, the climatology is made in the period of 1970/71-1999/2000. It should be noted that the use of other observation-based reanalysis, such as the EAR-40 from European Centre for Medium-Range Weather Forecasts, or the Japanese 55-year Reanalysis (JRA-55), leads to the same result. We also used the monthly mean temperature data of 160 stations from China Meteorological Administration for comparison. The same results were obtained (figures not shown).

CMIP5 models and experiments
Model Intercomparison Projects (MIPs) provide the best opportunity to evaluate how the current state-of-the-art climate models characterize the variability of the climate system. The most conspicuous MIP is probably the Coupled Model Intercomparison Project (CMIP), which compare simulations from global coupled climate models with its components describing the atmosphere, ocean, sea ice and land surface, and provide a standard set of model simulations aiming to evaluate how well the models simulate and project the future climate change and understand some of the factors responsible for the differences in model projections (http://cmip-pcmdi.llnl.gov/cmip5/). The simulation results were used in the climate change assessment report of IPCC. Since its inception in 1995, CMIP has now evolved to the fifth phase (CMIP5), which was widely used in the fifth assessment report of the IPCC.
In this study, we adopted 42 couple climate model outputs from the Program for Climate Model Diagnosis and Intercomparison (PCMDI, http://cmip-pcmdi.llnl.gov/). The names of the analyzed models are listed in Table 1, together with the model resolution of the atmosphere. Table 1 is compiled using information that the participating modeling centers provided to the PCMDI (see http://wwwpcmdi.llnl.gov/ipcc/about_ipcc.php for more information about the models). As stratosphere can provide an important source of model improvement, models are divided into two groups based on whether the full stratosphere is resolved. If the model lid top is higher than stratopause (1 hPa), it belongs to the high-top (HT) group; otherwise it is a low-top (LT) model. Accordingly, there are 20 HT models and 22 LT models utilized in this study.
We used long-term experiments of the historical, RCP4.5 and RCP8.5, three core simulations of CMIP5. Historical experiments are forced by observed atmospheric composition changes and including time-evolving land cover, and designed to inspect how well the model can reproduce the main aspects of the characteristics and variations of historical climate. RCPs (Representative Concentration Pathways) experiments are designed as future projection simulations forced by specific concentrations. The RCP4.5 represents additional radiative forcing increase to a level of about 4.5 W/m 2 in 2100 which is a midrange mitigation emission scenario, and the RCP8.5 is a high emission scenario with an additional radiative forcing increasing to about 8.5 W/ m 2 in 2100 (Taylor et al. 2012). The output parameters from CMIP5 models, including surface air temperature, precipitation rate, meridional wind at 1000 hPa, sea level pressure, geopotential height at 500 hPa, zonal wind in at 200 hPa, are used and compared with the observation.
To compare with previous version of CMIP and confirm the model progress, the model output from 24 CGCM, released by CMIP3 and used for the AR4 evaluation, are used in this study. The model names of the CMIP3 models are listed in Table 2. These model results were made available by the IPCC Data Distribution Center (ftp://ftp-esg. ucllnl.org) and the Hadley Centre for Climate Prediction and Research. To be consistent with the reanalysis datasets, all of the CMIP3 and CMIP5 simulation results were regrided to a horizontal grid resolution of 2.5° × 2.5° and focused on the winter (DJF) period.

The model climatology of EAWM
The comparison of CMIP5 historical run, RCP4.5 run and RCP8.5 run between HT and LT models, together with the NCEPR1 reanalysis and CMIP3 20C3 M run, is made for the mean state of winter (DJF) SAT (represented by 0 °C isotherm, Fig. 1a) and SLP (represented by the 1026 hPa isobar, Fig. 1b). Multi-model ensemble (MME) is made for each group. Black contours are for NCEPR1 reanalysis and purple solid contours are for CMIP3 models. The dashed contours are for the low-top models. As pointed out by WEI14, the CMIP5 models are capable of simulating a good and reasonable EAWM feature and great improvement has been obtained in characterizing the regional climatology, compared with the previous CMIP3 models. The 0 °C isotherms in CMIP5 are more northward and close to both the reanalysis and station observation (Fig. 1a of WEI14). Although the CMIP5 LT models have a better SAT simulation than that of the CMIP3 models, the CMIP HT models shows a much greater improvement. The 0 °C isotherm of HT models (solid green contour) are almost overlapping with the NCEPR1 reanalysis (black contour), while   Fig. 1 The comparison of CMIP5 historical run (green, 1970-1999 climatology), RCP45 run (blue, 2070-2099 climatology) and RCP85 run (red, 2070-2099 climatology) between HT models (solid contours) and LT models (dashed contours), together with CMIP3 20C3 M run (purple) and NCEPR1 reanalysis (black, 1970-1999 climatology) for EAWM fields: a the mean position of surface temperature (SAT) 0 °C isotherms; b sea level pressure represented by 1026-hPa isobars the cold bias is still obvious in the CMIP5 LT models. As the 0 °C line define the appropriate southern boundary of the snowfield, a better simulation and northward withdraw of the SAT field help to better characterize the hydrological cycle and radiation energy budget in EA region. The EAWM is characterized by a southward extension of surface cold air and the Siberian High downstream of the Tibetan Plateau over eastern China. The simulation of the SLP (Fig. 1b) shows that the CMIP3 MME is featured with too southward extension of the Siberian High over eastern China. The 1026-hPa isobar of CMIP3 is almost two degree south of the NCEPR1 reanalysis. The similar results can be obtained if utilizing another isobar such as 1024 or 1028 hPa (figures not shown). In CMIP5 LT model historical run (the green dashed line), the too southward extension of the Siberian High is maintained, indicating these models' deficiency in characterizing the sea level pressure distribution. While the southern boundary of the SLP 1026-hPa isobar in the CMIP5 HT historical MME is very close to observation, indicating the superiority of the Fig. 2 The DJF precipitation rate (unit: mm/day): a NOAA's Precipitation Reconstruction data (1970-1999 climatology); b CMIP3 MME (1970-1999 climatology); c the differences between CMIP5 models (historical run) and CMIP3 models (20C3 M run); d the differences between CMIP5 HT and LT models in historical run (1970-1999 climatology); e the differences between CMIP5 HT and LT models in RCP45 run (the 2070-2099 climatology); and f the differences between CMIP5 HT and LT models in RCP85 run (the 2070-2099 climatology). The units are mm/ day in all panels. Dots signify statistically significant regions above the 5% level according to the Student's t test HT models in characterizing the southern boundary of the Siberian High. However, the SLP contracts too westward in HT models, compared with the LT models and NCEPR1 reanalysis.
The DJF precipitation distribution in East Asia shows most rainfall in the southeast coastal region, and decreases northwestward and northward. The rain belt above 2 mm/ day extends northeastward from southeast China to Japan, with its eastern margin extending into the North Pacific storm track. The maximum rainfall centers above 3 mm/day over southwest Japan and central Japan. As WEI14 pointed out, the winter precipitation is overestimated by both CMIP3 and CMIP5 MME over most East Asian region, and the CMIP5 makes great improvement with less rainfall over the East Asia oceanic region of South China Sea, and Sea of Japan (Fig. 2c). The difference between CMIP5 HT and LT historical climatology (Fig. 2d) indicates that HT models have less rainfall over the South China Sea, East China Sea, and the ocean region south of Japan. Therefore, the precipitation amount and pattern are better simulated in the CMIP5 HT models. The overall improvement of the CMIP5 in simulation of precipitation was mainly caused by the improvement of the HT models.

Discrepancy of the EAWM between high-top and low-top CMIP5 models
In order to assess the performance of CMIP5 models with and without full stratosphere in charactering EAWM, we examined in several key parameters that are usually used as the strength or state of the EAWM. These parameters are listed in Table 3. In Fig. 3, we plotted the above indices in NCEPR1, historical climatology (1970), RCP45 projection (2070-2099), and RCP85 projection (2070-2099 for HT (red box) and LT (blue box) models. Systematic discrepancy can be observed between HT and LT models under all scenarios. For TAS over East Asian region (Fig. 3a), the HT models present higher temperature than the LT models, which is about 1.2, 1.3 and 1.5 °C higher than LT models under historical, RCP45 and RCP85. These differences are greater than the standard deviation of the NCEPR1 analysis, which is 0.9 °C. These values are nonignorable, compared with the projected change between the two periods of 1970-1999 and 2070-2099, which are 2.7 and 4.9 °C under RCP45 and RCP85 in HT models, and 2.6 and 4.6 °C under RCP45 and RCP85 in LT models.
As for the SLP averaged over the eastern China (I SLPEC , Fig. 2b), the values of SLP difference between the HT and LT models are much larger than the projected change under all RCPs and in both model groups. This is also true for the strength of Siberian High (I SH , Fig. 2d) and the lowlevel northerly along the southeast China coastal region (I V , Fig. 2f). As for the east-west SLP contrast (I SLPdiff , Fig. 2c), and the 200 hPa jet strength (I jet , Fig. 2e), the discrepancy between HT and LT models are comparable to the their individual projected change. Figure 4 shows the difference of SAT and SLP between HT and LT model under historical, RCP45 and RCP85 scenarios. Under all scenarios, the HT models generally simulated higher SAT and lower SLP over the East Asian region than the LT models. The SAT difference is above 1 °C over most land region. The maximum SAT difference is above 2 degree over Northeast China, Mongolia, and southwest China in historical simulation. It's interesting that the significant and large value regions are almost identical under all the scenarios, suggesting that the mechanism which causing the model discrepancy remains the same. However, the difference is larger over southern China under RCP85 than the RCP45, implying the possible amplifying of the difference between HT and LT by global warming.
The historical SLP difference shows consistent dropping between HT and LT models, while over the Siberian High (SH) center region (40°-60°N, 80°-120°E) the SLP I v The average low-level meridional wind along the EA coastal region The 1000 hPa northerly averaged for the domain: 15°-30°N, 115°-130°E Ji et al. (1997) difference value is small, indicating a northward contraction of the Siberian High in the HT models. While under RCPs, the SLP difference between HT and LT groups shows negative anomalies over the SH center region and East Asia, with the maximum center over North China.

The EAWM projection
The multi-model ensemble mean SAT (relative to the 1970-1999 base period) over the East Asian region (20°-60°N, 100°-140°E) through the 21st century is estimated with CMIP5 HT and LT models under the RCP45 (Fig. 5a) and RCP85 (Fig. 5b) scenarios. The projected SAT in 2100 is 3.2 °C and 1.9 °C respectively for HT and LT under RCP45, and is 6.3 and 4.6 °C respectively for HT and LT under RCP85. Through the 21st century, the average difference between HT and LT MME is 1.2 and 1.3 °C under RCP45 and RCP85, respectively. As the projection of regional and global temperature is usually based on all models regardless of the model stratosphere configuration, this kind of ensemble will get a MME projection around (c) (f) Fig. 4 The differences of surface air temperature (left panels) and sea level pressure (right panels) between CMIP5 HT and LT models under historical (a and d, upper panels, 1970-1999 climatology), RCP45 (b and e, middle panels, 2070-2099 climatology) and RCP85 (c and f, bottom panels, 2070-2099 climatology). Dots signify statistically significant regions above the 5% level according to the Student's t test the middle of HT and LT, which is 2.5 and 5.4 °C under RCP45 and RCP85, respectively. Therefore, the conventional MME results in a relatively lower temperature increase as compared with the full stratosphere-resolved HT models. This lower temperature increase will lead to an underestimation of the SAT change of around 0.7 °C and 0.9 °C under the global warming RCP45 and RCP85 scenarios, respectively. It's worth noting that the projected SAT change of the EA region in both HT and LT groups is comparable to the global SAT change under RCP45 [MME around 1.85 °C (Acollins et al. 2013)], while is much higher than the global SAT change under RCP85 [MME around 4.15 °C (Acollins et al. 2013)]. Therefore, it seems that the high emission of RCP85 will be amplified by the seasonal and regional effect. The warming can be extreme in some specific region and season, which worth further study.
The most unanimous change under the global warming scenario in various EAWM parameters (Fig. 3) is the surface air temperature (Fig. 3a) and high-level jet strength (Fig. 3e). In both HT and LT model group, the EA jet stream shows increase tendency, depending on the scenario with higher RCP forcing leading to stronger jet stream. Although the jet stream is considered an important index charactering the state of EAWM, with stronger jet stream corresponding to stronger EAWM, the EAWM doesn't show any systematic strengthening. Another critical component of the EAWM, the Siberian High, shows pronounced SLP decrease in the Siberian High center (Fig. 3d) under RCP45 and RCP85 in the HT models. Meanwhile, the 1000 hPa northerly along the EA coastal region, averaged for the domain of 15°-30°N, 115°-130°E, is featured with slightly decrease under RCP85 in the HT models. The zonal SLP gradient (Fig. 3c) doesn't have any discernible change under RCP85 and RCP45 in the HT models. We will discuss in the next section that this kind of jet increase is the response of anthropogenic global warming, instead of a local monsoon system change.

Discussion and conclusions
In this study, we compared the simulation and projection of the EAWM by two subensembles of the full CMIP5 ensemble, one containing 22 models with full stratosphere and the other containing 20 models with low model top below the stratopause. In the EA region, although the simulation of the EAWM by the two ensembles is similar and better than that in the CMIP3 models, the subensemble with the full stratosphere (HT models) has a better simulation of the temperature (the cold bias is removed), precipitation (the oceanic rainfall surplus is reduced), and the southward extension of the Siberian High (the too southward extension is eliminated). Therefore, the better simulation of EA regional winter climate depends largely on the inclusion of a full stratosphere in the coupled model.
With a full stratosphere, the boundary effect in the stratosphere can be alleviated and a real stratospheric polar vortex can be obtained. Studies revealed that, with high-top model, the common cold bias of the stratospheric polar vortex is corrected and the strength of the polar jet is weakened to be comparable to observation (Charlton-Perez et al. 2013). Figure 6 shows the difference of geopotential height at various vertical levels between HT and LT model under The differences of winter (DJF) geopotential height at a 30 hPa, b 100 hPa, c 850 hPa, and d sea level pressure between CMIP5 HT and LT models under historical run. The winter climatol-goy is based on the data between 1970 and 1999. Dots signify statistically significant regions above the 5% level according to the Student's t test the historical run. The meridional geopotential height gradient is decreased around the stratospheric polar vortex in the MME of HT models, with height increase in the polar cap and significant height decrease in the middle latitudes and the sub-polar region. This north-south geopotential height anomaly makes a negative North Annular Model (NAM) pattern and has a deep quasi-barotropic influence through the height of stratosphere and troposphere. At the lower levels of 850 hPa and the sea level, the polar region pressure is increased, while the pressure around the middle latitude is decreased. Meanwhile, at the middle stratosphere level (Fig. 6a), a wavy pattern of wavenumber 1 is apparent. And at the low tropospheric level, the SLP over Siberia, East Asia, North America, and North Atlantic is decreased, making a wavy pattern of wavenumber 2. The SLP decreases over the Siberia and East Asia region. The circulation change favors the weakening of the EAWM and therefore leads to the SAT increase in the East Asian region. The zonal-mean zonal wind difference (Fig. 7a) between HT and LT models shows zonal wind reduction around the polar throughout the stratosphere and troposphere, while wind increase in the upper troposphere and lower stratosphere around 35°N, centering on the subtropical tropopause. This wind change may be associated with the change of planetary wave in the stratosphere-troposphere system. The HT ensemble has a stronger stratospheric wave-guide than the LT ensemble, while its tropospheric wave-guide is slightly weaker than the LT models. According to Chen et al. (2003), this kind of dipole structure in Eliassen-Palm (EP) flux (Andrews et al. 1987;Edmon et al. 1980) anomalies could lead to the weakening of the polar night jet and the strengthening of the subtropical jet. The polar jet and subtropical jet change under the RCPs have the similar pattern, indicating that the dynamical mechanism of including a full stratosphere remains the same. However, it should bet noted that model's sensitivity to increase of greenhouse gases may increase the subtropical jet. The greenhouse gas-associated global warming causes a mixing of tropospheric warming and stratospheric cooling. Since the tropopause has an average height of ~17 km in the tropics and only ~9 km at the poles, both the temperature and pressure gradients increases near the upper troposphere and lower stratosphere in the subtropical regions, leading to zonal wind acceleration from upper troposphere to middle stratosphere there. The subtropical jet increase of the HT models indicates that the models with a full stratosphere are more capable of capturing the GHG-associated global change signal. Accordingly, the subtropical jet gets stronger in the HT models, especially under the high emission RCP85 scenario.
Currently, the models with full stratosphere are usually featured with higher model top above stratopause, as well as non-orographic gravity waves (NOGW). In this study, 17 of the 20 HT models considered NOGW, and only 5 of the 22 LT models have the NOGW. Meanwhile, most (12 of the 20 models) of the HT models considered the carbon cycle and natural aerosols and therefore belongs to Earth system models (ESMs), while half (11 of the 22 models) of the LT models are ESMs. Therefore, the inclusion of the NOGW and upgrade to ESMs may also help to improve the simulation of stratospheric polar vortex and the EAWM.
The projection under RCP scenarios shows significant discrepancy between HT and LT models. Under the (a) (b) (c) Fig. 7 The zonal-mean zonal wind (shadings) and Eliassen-Palm (EP) flux (vectors) between CMIP5 HT and LT models under a historical (1970-1999 climatology), RCP45 (2070RCP45 ( -2099 and RCP85 (2070RCP85 ( -2099. Dots signify statistically sig-nificant regions above the 5% level according to the Student's t test for the zonal-mean zonal wind. The EP fluxes are scaled by the inverse of the square root of the air density and the unit is kg s −2 midrange emission RCP45 scenario, the projected SAT change in 2100 over the EA region is 3.2 and 1.9 °C for HT and LT models, respectively. And it's 6.3 °C (HT) and 4.6 °C (LT) in 2100 under the higher emission RCP85 scenario. Therefore, the conventional MME using both HT and LT models will lead to lower temperature increase than the full stratosphere-resolved models ensemble. As climate models have now become one of the primary tools in climate change assessment and projection, which was used as the base for mitigation measure and adaption policy. The systemic underestimation of the temperature change in the LT climate models have the potential to lead to insufficient response action to the anthropogenic global warming in the next decades. Although our study focuses on the EA region in the cold season, it's worth to note that the stratosphere configuration in the climate model can have significant influence in other seasons and regions. The involvement of a full stratosphere results in a better simulation of the global-scale stratospheric polar vortex, which can have global influences. And the amplitude of global change under RCP scenarios may be different from current evaluation. This will be discussed in another study.