The Role of Underlying Boundary Forcing in Shaping the Recent Decadal Change of Persistent Anomalous Activity over the Ural Mountains

Observational analyses demonstrate that the Ural persistent positive height anomaly event(PAE) experienced a decadal increase around the year 2000, exhibiting a southward displacement afterward. These decadal variations are related to a large-scale circulation shift over the Eurasian Continent. The effects of underlying sea ice and sea surface temperature(SST) anomalies on the Ural PAE and the related atmospheric circulation were explored by Atmospheric Model Intercomparison Project(AMIP) experiments from the Coupled Model Intercomparison Project Phase 6 and by sensitivity experiments using the Atmospheric General Circulation Model(AGCM). The AMIP experiment results suggest that the underlying sea ice and SST anomalies play important roles. The individual contributions of sea ice loss in the Barents-Kara Seas and the SST anomalies linked to the phase transition of the Pacific Decadal Oscillation(PDO) and Atlantic Multidecadal Oscillation(AMO) are further investigated by AGCM sensitivity experiments isolating the respective forcings.The sea ice decline in Barents-Kara Seas triggers an atmospheric wave train over the Eurasian mid-to-high latitudes with positive anomalies over the Urals, favoring the occurrence of Ural PAEs. The shift in the PDO to its negative phase triggers a wave train propagating downstream from the North Pacific. One positive anomaly lobe of the wave train is located over the Ural Mountains and increases the PAE there. The negative-to-positive transition of the AMO phase since the late-1990s causes positive 500-h Pa height anomalies south of the Ural Mountains, which promote a southward shift of Ural PAE.

Local Torrential Rainfall Event within a Mei-Yu Season Mesoscale Convective System:Importance of Back-Building Processes

An extreme rainfall event occurred over Hangzhou,China,during the afternoon hours on 24 June 2013.This event occurred under suitable synoptic conditions and the maximum 4-h cumulative rainfall amount was over 150 mm.This rainfall event had two major rainbands.One was caused by a quasi-stationary convective line,and the other by a backbuilding convective line related to the interaction of the outflow boundary from the first rainband and an existing low-level mesoscale convergence line associated with a mei-yu frontal system.The rainfall event lasted 4 h,while the back-building process occurred in 2 h when the extreme rainfall center formed.So far,few studies have examined the back-building processes in the mei-yu season that are caused by the interaction of a mesoscale convergence line and a convective cold pool.The two rainbands are successfully reproduced by the Weather Research and Forecasting(WRF)model with fourlevel,two-way interactive nesting.In the model,new cells repeatedly occur at the west side of older cells,and the backbuilding process occurs in an environment with large CAPE,a low LFC,and plenty of water vapor.Outflows from older cells enhance the low-level convergence that forces new cells.High precipitation efficiency of the back-building training cells leads to accumulated precipitation of over 150 mm.Sensitivity experiments without evaporation of rainwater show that the convective cold pool plays an important role in the organization of the back-building process in the current extreme precipitation case.

Distribution and Formation Causes of PM_(2.5) and O_(3) Double High Pollution Events in China during 2013–20

Fine particulate matter(PM_(2.5))and ozone(O_(3))double high pollution(DHP)events have occurred frequently over China in recent years,but their causes are not completely clear.In this study,the spatiotemporal distribution of DHP events in China during 2013–20 is analyzed.The synoptic types affecting DHP events are identified with the Lamb–Jenkinson circulation classification method.The meteorological and chemical causes of DHP events controlled by the main synoptic types are further investigated.Results show that DHP events(1655 in total for China during 2013–20)mainly occur over the North China Plain,Yangtze River Delta,Pearl River Delta,Sichuan Basin,and Central China.The occurrence frequency increases by 5.1%during 2013–15,and then decreases by 56.1%during 2015–20.The main circulation types of DHP events are“cyclone”and“anticyclone”,accounting for over 40%of all DHP events over five main polluted regions in China,followed by southerly or easterly flat airflow types,like“southeast”,“southwest”,and“east”.Compared with non-DHP events,DHP events are characterized by static or weak wind,high temperature(20.9℃ versus 23.1℃)and low humidity(70.0%versus 64.9%).The diurnal cycles of meteorological conditions cause PM_(2.5)(0300–1200 LST,Local Standard Time=UTC+8 hours)and O_(3)(1500–2100 LST)to exceed the national standards at different periods of the DHP day.Three pollutant conversion indices further indicate the rapid secondary conversions during DHP events,and thus the concentrations of NO_(2),SO_(2) and volatile organic compounds decrease by 13.1%,4.7%and 4.4%,respectively.The results of this study can be informative for future decisions on the management of DHP events.

Effect of the Initial Vortex Structure on Intensity Change During Eyewall Replacement Cycle of Tropical Cyclones:A Numerical Study

This study investigates the effect of the initial tropical cyclone(TC)vortex structure on the intensity change during the eyewall replacement cycle(ERC)of TCs based on two idealized simulations using the Weather Research and Forecasting(WRF)model.Results show that an initially smaller TC with weaker outer winds experienced a much more drastic intensity change during the ERC than an initially larger TC with stronger outer winds.It is found that an initially larger TC vortex with stronger outer winds favored the development of more active spiral rainbands outside the outer eyewall,which slowed down the contraction and intensification of the outer eyewall and thus prolonged the duration of the concentric eyewall and slow intensity evolution.In contrast,the initially smaller TC with weaker outer winds corresponded to higher inertial stability in the inner core and weaker inertial stability but stronger filamentation outside the outer eyewall.These led to stronger boundary layer inflow,stronger updraft and convection in the outer eyewall,and suppressed convective activity outside the outer eyewall.These resulted in the rapid weakening during the formation of the outer eyewall,followed by a rapid re-intensification of the TC during the ERC.Our study demonstrates that accurate initialization of the TC structure in numerical models is crucial for predicting changes in TC intensity during the ERC.Additionally,monitoring the activity of spiral rainbands outside the outer eyewall can help to improve short-term intensity forecasts for TCs experiencing ERCs.

The Impact of Atmospheric Heat Sources over the Eastern Tibetan Plateau and the Tropical Western Pacific on the Summer Rainfall over the Yangtze-River Basin

The variability of the summer rainfall over China is analyzed using the EOF procedure with a new parameter （namely, mode station variance percentage） based on 1951-2000 summer rainfall data from 160 stations in China. Compared with mode variance friction, the mode station variance percentage not only reveals more localized characteristics of the variability of the summer rainfall, but also helps to distinguish the regions with a high degree of dominant EOF modes representing the analyzed observational variable. The atmospheric circulation diagnostic studies with the NCEP/NCAR reanalysis daily data from 1966 to 2000 show that in summer, abundant （scarce） rainfall in the belt-area from the upper-middle reaches of the Yangtze River northeastward to the Huaihe River basin is linked to strong （weak） heat sources over the eastern Tibetan Plateau, while the abundant （scarce） rainfall in the area to the south of the middle-lower reaches of the Yangtze River is closely linked to the weak （strong） heat sources over the tropical western Pacific.

Identifying Global Monsoon Troughs and Global Atmospheric Centers of Action on a Pentad Scale

Using two datasets of global pentad grid precipitation and global 850 hPa geopotential height during 1979-2007,this study identified global monsoon troughs and global atmospheric centers of action (ACAs) on a pentad scale.The global monsoon troughs consist of planetary-scale monsoon troughs and peninsula-scale monsoon troughs.Forced by seasonal variations in solar radiation,the inter-tropical convergence zones (ITCZs) represent the planetary-scale monsoon troughs,which are active and shift over the tropical North Pacific,the tropical North Atlantic,and the tropical South Indian oceans.The peninsula-scale monsoon troughs are originated from regional land-sea topography and varied with contrasts in seasonal land-sea surface temperatures and precipitation.During the boreal summer,five peninsula-scale troughs and one planetary-scale trough are distributed in the Asia-Northwest Pacific (NWP) region.In total,22 troughs,nine monsoon troughs,and 19 ACAs in the lower troposphere were identified.Relevant ACAs may be useful in constructing regional monsoon and circulation indices.

Performance of a Reconfigured Atmospheric General Circulation Model at Low Resolution

Paleoclimate simulations usually require model runs over a very long time. The fast integration version of a state-of-the-art general circulation model （GCM）, which shares the same physical and dynamical processes but with reduced horizontal resolution and increased time step, is usually developed. In this study, we configure a fast version of an atmospheric GCM （AGCM）, the Grid Atmospheric Model of IAP/LASG （Institute of Atmospheric Physics/State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics）, at low resolution （GAMIL-L, hereafter）, and compare the simulation results with the NCEP/NCAR reanalysis and other data to examine its performance. GAMIL-L, which is derived from the original GAMIL, is a finite difference AGCM with 72 × 40 grids in longitude and latitude and 26 vertical levels. To validate the simulated climatology and variability, two runs were achieved. One was a 60-year control run with fixed climatological monthly sea surface temperature （SST） forcing, and the other was a 50-yr （1950-2000） integration with observational time-varying monthly SST forcing. Comparisons between these two cases and the reanalysis, including intra-seasonal and inter-annual variability are also presented. In addition, the differences between GAMIL-L and the original version of GAMIL are also investigated.The results show that GAMIL-L can capture most of the large-scale dynamical features of the atmosphere, especially in the tropics and mid latitudes, although a few deficiencies exist, such as the underestimated Hadley cell and thereby the weak strength of the Asia summer monsoon. However, the simulated mean states over high latitudes, especially over the polar regions, are not acceptable. Apart from dynamics, the thermodynamic features mainly depend upon the physical parameterization schemes. Since the physical package of GAMIL-L is exactly the same as the original high-resolution version of GAMIL, in which the NCAR Community Atmosphere Model （CAM2） physical package was used, there are only small differences between them in the precipitation and temperature fields. Because our goal is to develop a fast-running AGCM and employ it in the coupled climate system model of IAP/LASG for paleoclimate studies such as ENSO and Australia-Asia monsoon, particular attention has been paid to the model performances in the tropics. More model validations, such as those ran for the Southern Oscillation and South Asia monsoon, indicate that GAMIL-L is reasonably competent and valuable in this regard.

An Observational Analysis of the Oceanic and Atmospheric Structure of Global-Scale Multi-decadal Variability

The aim of the present study was to identify multi-decadal variability （MDV） relative to the current centennial global warming trend in available observation data.The centennial global wanning trend was first identified in the global mean surface temperature （STgm） data.The MDV was identified based on three sets of climate variables,including sea surface temperature （SST）,ocean temperature from the surface to 700 m,and the NCEP and ERA40 reanalysis datasets,respectively.All variables were detrended and low-pass filtered.Through three independent EOF analyses of the filtered variables,all results consistently showed two dominant modes,with their respective temporal variability resembling the Pacific Decadal Oscillation/Inter-decadal Pacific Oscillation （PDO/IPO） and the Atlantic Multi-decadal Oscillation （AMO）.The spatial structure of the PDO-like oscillation is characterized by an ENSO-like structure and hemispheric symmetric features.The structure associated with the AMO-like oscillation exhibits hemispheric asymmetric features with anomalous warm air over Eurasia and warm SST in the Atlantic and Pacific basin north of 10°S,and cold SST over the southern oceans.The Pacific and Atlantic MDV in upper-ocean temperature suggest that they are mutually linked.We also found that the PDO-like and AMO-like oscillations are almost equally important in global-scale MDV by EOF analyses.In the period 1975-2005,the evolution of the two oscillations has given rise to strong temperature trends and has contributed almost half of the STgm warming.Hereon,in the next decade,the two oscillations are expected to slow down the global warming trends.

Atmospheric Profile Retrieval with AIRS Data and Validation at the ARM CART Site

The physical retrieval algorithm of atmospheric temperature and moisture distribution from the Atmospheric InfraRed Sounder （AIRS） radiances is presented. The retrieval algorithm is applied to AIRS clear-sky radiance measurements. The algorithm employs a statistical retrieval followed by a subsequent nonlinear physical retrieval. The regression coefficients for the statistical retrieval are derived from a dataset of global radiosonde observations （RAOBs） comprising atmospheric temperature, moisture, and ozone profiles. Evaluation of the retrieved profiles is performed by a comparison with RAOBs from the Atmospheric Radiation Measurement （ARM） Program Cloud And Radiation Testbed （CART） in Oklahoma, U. S. A.. Comparisons show that the physically-based AIRS retrievals agree with the RAOBs from the ARM CART site with a Root Mean Square Error （RMSE） of 1K on average for temperature profiles above 850 hPa, and approximately 10% on average for relative humidity profiles. With its improved spectral resolution, AIRS depicts more detailed structure than the current Geostationary Operational Environmental Satellite （GOES） sounder when comparing AIRS sounding retrievals with the operational GOES sounding products.

Structural Variation of an Atmospheric Heat Source over the Qinghai-Xizang Plateau and Its Influence on Precipitation in Northwest China

NCEP/NCAR reanalysis data and a 47-year precipitation dataset are utilized to analyze the relationship between an atmospheric heat source （hereafter called 〈 Q1 〉） over the Qinghai-Xizang Plateau （QXP） and its surrounding area and precipitation in northwest China. Our main conclusions are as follows： （1） The horizontal distribution of 〈 Q1 〉 and its changing trend are dramatic over QXP in the summer. There are three strong centers of 〈 Q1 〉 over the south side of QXP with obvious differences in the amount of yearly precipitation and the number of heat sinks predominate in the arid and semi-arid regions of northwest China （NWC）, beside the northern QXP with an obvious higher intensity in years with less precipitation. （2） In the summer, the variation of the heat source＇s vertical structure is obviously different between greater and lesser precipitation years in eastern northwest China （ENWC）. The narrow heat sink belt forms between the northeast QXP and the southwestern part of Lake Baikal. In July and August of greater precipitation years, the heating center of the eastern QXP stays nearly over 35°N, and at 400 hPa of the eastern QXP, the strong upward motion of the heating center constructs a closed secondary vertical circulation cell over the northeast QXP （40~ 46~N）, which is propitious to add precipitation over the ENWC. Otherwise, the heating center shifts to the south of 30°N and disappears in July and August of lesser precipitation years, an opposite secondary circulation cell forms over the northeast QXP, which is a disadvantage for precipitation. Meanwhile, the secondary circulation cell in years with more or less precipitation over the ENWC is also related to the heat source over the Lake Baikal. （3） The vertical structure of the heat source over the western QXP has obvious differences between greater and lesser precipitation years in western northwest China in June and July. The strong/weak heat source over the western QXP produces relatively strong/weak ascending motion and correspondingly constructs a secondary circulation cell in lesser/greater precipitation years.

Seasonal evolution of the dominant modes of the Eurasian snowpack and atmospheric circulation from autumn to the subsequent spring and the associated surface heat budget

This study investigates the seasonal evolution of the dominant modes of the Eurasian snowpack and atmospheric circulation from autumn to the subsequent spring using snow water equivalent （SWE）, snow cover frequency （SCF）, and 500 hPa geopotential height data. It is found that the Eurasian SWE/SCF and circulation dominant modes are stably coupled from autumn to the subsequent spring.The temporal coherence of the seasonal evolution of the dominant modes is examined.The seasonal evolution of the Eurasian circulation and SWE dominant modes exhibit good coherence, whereas the evolution of the Eurasian SCF dominant mode is incoherent during the autumn-winter transition season. This incoherence is associated with a sign-change in the SCF anomalies in Europe during the autumn-winter transition season, which is related to the wind anomalies over Europe. In addition, the surface heat budget associated with the Eurasian SWE/SCF and circulation dominant modes is analyzed. The sensible heat flux （SHF） related to the wind-induced thermal advection dominates the surface heat budget from autumn to the subsequent spring, with the largest effect during winter. The surface net shortwave radiation is mainly modulated by snow cover rather than cloud cover, which is estimated to be as important as, or likely superior to, the SHF for the surface heat budget during spring.The NCEP-NCAR surface heat flux reanalysis data demonstrate a consistency with the SWE/SCF and air temperature observational data, indicating a good capability for snow-atmosphere interaction analysis.

ATMOSPHERIC BOUNDARY LAYER CONCEPT MODEL OF THE PEARL RIVER DELTA AND ITS APPLICATION

Based on the geographical circumstance, climate and the boundary layer meteorology features of the Pearl River Delta, a boundary layer concept model of the Pearl River Delta was built. The concept model consists of four fundamental factors that affect the boundary layer meteorology of the Pearl River Delta and can convincingly explain the reason of the air quality change in the Pearl River Delta. The model can be used to the diffusion capability analysis, the air pollution potential forecasting or haze forecasting, etc.

RELATIONSHIP BETWEEN THE 30 TO 60 DAY OSCILLATION OF ATMOSPHERIC HEAT SOURCE AND THE DROUGHT AND FLOOD EVENTS IN JUNE IN THE SOUTH OF CHINA

Based on the NCEP/NCAR reanalysis data and the observed precipitation data in the south of China from 1958 to 2000,the impact of 30 to 60 day oscillation of atmospheric heat sources on the drought and flood events in June in the south of China is discussed.During the flood(drought) events,there exists an anomalous low-frequency anticyclone(cyclone) at the low level of the troposphere over the South China Sea and the northwestern Pacific,accompanied with anomalous low-frequency heat sinks(heat sources),while there exists an anomalous low-frequency cyclone(anticyclone) with anomalous heat sources(sinks) over the area from the south of China to the south of Japan.On average,the phase evolution of the low-frequency in drought events is 7 to 11 days ahead of that in flood events in May to June in the south of China.In flood events,low-frequency heat sources and cyclones are propagated northward from the southern South China Sea,northwestward from the warm pool of the western Pacific and westward from the northwestern Pacific around 140°E,which have very important impact on the abundant rainfall in June in the south of China.However,in drought events,the northward propagations of the low-frequency heat sources and cyclones from the South China Sea and its vicinity are rather late compared with those in flood events,and there is no obvious westward propagation of the heat sources from the northwestern Pacific.The timing of the low-frequency heat source propagation has remarkable impact on the June rainfall in the south of China.

Retrieval of Atmospheric and Oceanic Parameters and the Relevant Numerical Calculation

It is well known that retrieval of parameters is usually ill-posed and highly nonlinear, so parameter retrieval problems are very difficult. There are still many important theoretical issues under research, although great success has been achieved in data assimilation in meteorology and oceanography. This paper reviews the recent research on parameter retrieval, especially that of the authors. First, some concepts and issues of parameter retrieval are introduced and the state-of-the-art parameter retrieval technology in meteorology and oceanography is reviewed briefly, and then atmospheric and oceanic parameters are retrieved using the variational data assimilation method combined with the regularization techniques in four examples： retrieval of the vertical eddy diffusion coefficient; of the turbulivity of the atmospheric boundary layer; of wind from Doppler radar data, and of the physical process parameters. Model parameter retrieval with global and local observations is also introduced.

Application of an Artificial Neural Network for a Direct Estimation of Atmospheric Instability from a Next-Generation Imager

Atmospheric instability information derived from satellites plays an important role in short-term weather forecasting, especially the forecasting of severe convective storms. For the next generation of weather satellites for Korea＇s multi-purpose geostationary satellite program, a new imaging instrument has been developed. Although this imaging instrument is not de- signed to perform full sounding missions and its capability is limited, its multi-spectral infrared channels provide information on vertical sounding. To take full advantage of the observation data from the much improved spatiotemporal resolution of the imager, the feasibility of an artificial neural network approach for the derivation of the atmospheric instability is investigated. The multi-layer perceptron model with a feed-forward and back-propagation training algorithm shows quite a sensitive re- sponse to the selection of the training dataset and model architecture. Through an extensive performance test with a carefully selected training dataset of 7197 independent profiles, the model architectures are selected to be 12, 5000, and 0.3 for the number of hidden nodes, number of epochs, and learning rate, respectively. The selected model gives a mean absolute error, RMSE, and correlation coefficient of 330 J kg-1, 420 J kg-1, and 0.9, respectively. The feasibility is further demonstrated via application of the model to real observation data from a similar instrument that has comparable observation channels with the planned imager.

Evaluation and Improvement of a SVD-Based Empirical Atmospheric Model

An empirical atmospheric model(EAM) based on the singular value decomposition(SVD) method is evaluated using the composite El Ni(?)o/Southern Oscillation(ENSO) patterns of sea surface temperature (SST) and wind anomalies as the target scenario.Two versions of the SVD-based EAM were presented for comparisons.The first version estimates the wind anomalies in response to SST variations based on modes that were calculated from a pair of global wind and SST fields(i.e.,conventional EAM or CEAM).The second version utilizes the same model design but is based on modes that were calculated in a region-wise manner by separating the tropical domain from the remaining extratropical regions(i.e.,region-wise EAM or REAM). Our study shows that,while CEAM has shown successful model performance over some tropical areas, such as the equatorial eastern Pacific(EEP),the western North Pacific(WNP),and the tropical Indian Ocean(TIO),its performance over the North Pacific(NP) seems poor.When REAM is used to estimate the wind anomalies instead of CEAM,a marked improvement over NP readily emerges.Analyses of coupled modes indicate that such an improvement can be attributed to a much stronger coupled variability captured by the first region-wise SVD mode at higher latitudes compared with that captured by the conventional one. The newly proposed way of constructing the EAM(i.e.,REAM) can be very useful in the coupled studies because it gives the model a wider application beyond the commonly accepted tropical domain.

NUMERICAL EXPERIMENTS OF EFFECT OF SSTA OVER THE INDIAN OCEAN ON ATMOSPHERIC LOW-FREQUENCY OSCILLATION IN THE EXTRATROPICAL LATITUDE

Numerical experiments on forcing dissipation and heating response of dipole (unipole) are carried out using global spectral models with quasi-geostrophic barotropic vorticity equations. For each experiment model integration is run for 90 days on the condition of three-wave quasi-resonance. The results are given as follows: Under the effects of dipole (unipole) forcing source and basic flow intensity, there exist strong interactions among the three planetary waves and quasi-biweekly and intraseasonal oscillation of the three planetary waves. In the meantime, the changes in the intensity of dipole or unipole forcing source and basic flow have different frequency modulation effects on LFO in the middle and higher latitudes. The results of the stream function field of three quasi-resonant waves evolving with time confirm that the low-frequency oscillation exists in extratropical latitude.

ON THE STUDY OF THE RELATIONSHIPS BETWEEN GUANGZHOU ATMOSPHERIC ENVIRONMENT FACTORS AND THE SARS EPIDEMIC

Based on SARS epidemic data and the corresponding atmospheric data, we used the timescale-partitioning technique, spectrum analysis and correlation analysis to investigate the impacts of the atmospheric environmental factors on the SARS epidemic. Results showed that there were close relations between environmental factors and SARS: The daily probable cases of SARS varied in 3-5 day cycles, much the same as the atmospheric elements did. The variations of the epidemics correlated remarkably with atmospheric elements. So conclusions can be drawn that weather changes have influences on the variations of daily SARS cases. In addition, statistical results showed that cold air activities aggravated the SARS epidemic.

Causes of the record-low Antarctic sea-ice in austral summer 2022

2022年夏季(2021年12月至2022年2月)南极海冰面积达到历史新低,西南极减少最显著.2021年8~10月南半球环状模接近历史最强和7~9月海洋性大陆附近海温显著增暖是两个关键因素.由于平流层臭氧破纪录减少使得前者维持历史最强或接近最强,导致阿蒙森低压(ASL)加深并向西南移动,有利于海冰减少.后者持续到夏季,有利于拉尼娜的发展,通过激发罗斯贝波列加深ASL.在热力上,臭氧减少导致向下净短波辐射增加,引起西南极增暖.此外,净短波辐射-海温-云形成正反馈,与埃克曼输送一起,放大表面增暖,从而促进海冰融化.

Role of Stratospheric Processes in Climate Change: Advances and Challenges

In this review,instead of summarizing all the advances and progress achieved in stratospheric research,the main advances and new developments in stratosphere-troposphere coupling and stratospheric chemistry-climate interactions are summarized,and some outstanding issues and grand challenges are discussed.A consensus has been reached that the stratospheric state is an important source of improving the predictability of the troposphere on sub-seasonal to seasonal(S2S)time scales and beyond.However,applying stratospheric signals in operational S2S forecast models remains a challenge because of model deficiencies and the complexities of the underlying mechanisms of stratosphere-troposphere coupling.Stratospheric chemistry,which controls the magnitude and distribution of many important climate-forcing agents,plays a critical role in global climate change.Convincing evidence has been found that stratospheric ozone depletion and recovery have caused significant tropospheric climate changes,and more recent studies have revealed that stratospheric ozone variations can even exert an impact on SSTs and sea ice.The climatic impacts of stratospheric aerosols and water vapor are also important.Although their quantitative contributions to radiative forcing have been reasonably well quantified,there still exist large uncertainties in their long-term impacts on climate.The advances and new levels of understanding presented in this review suggest that whole-atmosphere interactions need to be considered in future for a better and more thorough understanding of stratosphere-troposphere coupling and its role in climate change.