Search

The teleconnection between the quasi-biennial oscillation (QBO) and the Arctic stratospheric polar vortex, or the Holton–Tan (HT) relationship, may change in a warmer climate or one with stratospheric aerosol intervention (SAI) compared to the present-day climate (PDC). Our results from an Earth system model indicate that, under both global warming (based on RCP8.5 emission scenario) and SAI scenarios, the HT relationship weakens in early winter (November–December), although it is closer to PDC under SAI than under the RCP8.5 scenario. In contrast, the HT relationship in the middle to late winter period (January–February) does not change considerably in response to either RCP8.5 or SAI scenarios compared to PDC. While the weakening of the HT relationship under the RCP8.5 scenario is likely due to the weaker QBO wind amplitudes at the Equator, another physical mechanism must be responsible for the weaker HT relationship under SAI scenarios, since the amplitude of the QBO wind is comparable to the PDC. The strength of the polar vortex does not change under the RCP8.5 scenario compared to PDC, but it becomes stronger under SAI; we attribute the weakening of the HT relationship under SAI to a stronger polar vortex. In general, the changes in the HT relationship cannot be explained by changes to the critical line; the changes in the residual circulation (particularly due to the gravity wave contributions) are important in explaining the changes in the HT relationship under RCP8.5 and SAI scenarios.

We analyse sporadic E (ES) layer occurrence rates (OR) obtained from ionospheric GPS radio occultation measurements by the FORMOSAT-3/COSMIC constellation. Maximum OR are seen at 95–105 km altitude. Midlatitude ES layers are mainly due to wind shear in the presence of tides, and the strongest signals are the migrating diurnal and semidiurnal components. Especially in the Southern Hemisphere, nonmigrating components such as a diurnal westward wave 2 and a semidiurnal westward wave 1 are also visible, especially at higher latitudes. Near the equator, a strong diurnal eastward wavenumber 3 component and a semidiurnal eastward wavenumber 2 component occur in summer and autumn. Terdiurnal and quarterdiurnal components are weaker than the diurnal and semidiurnal ones.

Solar EUV radiation is the dominant driver for upper atmosphere ionization. Ionospheric variations that affect radio signal propagation and thus affect technical systems such as satellite-based positioning systems. One significant time scale for the solar variability is the solar 27-day rotation period that causes a corresponding response in ionospheric observables like the height-dependent electron density (Ne) or the integrated total electron content (TEC). To enhance our understanding of the processes within the ionosphere we investigate a combination of observations and physical models, namely the Coupled Thermosphere Ionosphere Plasmasphere electrodynamics (CTIPe) model and the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIE-GCM), are analyzed to identify differences in Ne, TEC data and ionized oxygen (O+ and O2+). Furthermore, the model results are compared to ground-based ionosonde Ne measurements with regard to the spatial and temporal response of the ionosphere to the 27-day solar rotation period. Modeled Ne correlates strongly with the observed Ne at mid-latitudes, but at low-latitudes the modeled TEC distribution follows the geomagnetic coordinates more strictly when compared to the observational data. Local TEC and F2 layer peak Ne are well represented by CTIPe, whereas TIE-GCM represents the global TEC and F2 layer peak height well.

Winds at the mesosphere and lower thermosphere have been measured by partial reflection radars and specular meteor radars for almost two decades (2004-2022) over Germany and Norway (i.e., middle and high latitudes, respectively). Continuous wind measurements during the mentioned period are important to understand their long-term behavior. The zonal mean wind climatology displays an eastward wind during the winter months and a westward summer jet below ~85km at middle latitudes (~90km at high latitudes). Above the mentioned height, an eastward wind jet is observed. In the meridional wind component, the southward summer wind displays amplitudes between 4 and 5 times less intense than the westward jet. We studied the intensity of the summer wind components, the long-term variability and the possible connection to external forcing (i.e. El Niño-Southern Oscillation, and quasi-biennial oscillation, solar activity and geomagnetic activity). Analyzing the summer winds for low and high geomagnetic activity classified with the Ap index, there is a significant difference between both cases suggesting disturbances in the wind due to high geomagnetic activity. The long-term study shows significant trends at middle latitudes in the monthly summer values of the westward summer jet. As a consequence of the increase in the westward wind, a decrease in the southward component is observed at the same latitudes. While at high latitudes the eastward jet shows a decreasing velocity during July.

The thermosphere-ionosphere system changes significantly on various temporal scales due to the forcings from solar and geomagnetic processes, and the lower atmosphere. The 27-day variation caused by solar rotation is one of the most important modulating factors in the ionosphere. A robust feature in this context is the ionospheric lag of about 1-2 days in ionospheric parameters such as total electron content (TEC) and F2 layer peak electron density with respect to solar variations at the 27-day solar rotation period. Here, the ionospheric TEC provided by the International GNSS Service (IGS) and the Coupled Thermosphere Ionosphere Plasmasphere electrodynamics (CTIPe) model were used to understand the delayed ionospheric response and the underlying physics. The O/N2 measurements from the imaging ultraviolet spectrograph Global-Scale Observations of the Limb and Disk (GOLD), the Global Ultraviolet Imager (GUVI), were also analyzed during the 2019-2021 period of low solar activity. The comparative study shows that the model successfully reproduces the delayed response of the ionosphere during low solar activity. The observed and modeled O/N2 ratio was found to be positively correlated with the solar EUV proxy (GOLD QEUV), with a lag of about 2 days, indicating a contribution to the ionospheric lag in TEC.Furthermore, the CTIPe model simulations show that the ionospheric time delay is significantly affected by various physical processes such as diffusion, photodissociation, solar and geomagnetic activities, and wave dynamics.

This study investigates the response of the semidiurnal tide (SDT) to the 2013 major sudden stratospheric warming (SSW) event using meteor radar wind observations and mechanistic tidal model simulations. In the model, the background atmosphere is constrained to meteorological fields from the Navy Global Environmental Model - High Altitude analysis system. The solar (thermal) and lunar (gravitational) SDT components are forced by incorporating hourly global temperature tendency fields from the ERA5 forecast model, and by specifying the M2 and N2 lunar gravitational potentials, respectively. The simulated SDT response is compared against meteor wind observations from the CMOR (43.3◦N, 80.8◦W), Collm (51.3◦N, 13.0◦E), and Kiruna (67.5◦N, 20.1◦E) radars, showing close agreement with the observed amplitude and phase variability. Numerical experiments investigate the individual roles of the solar and lunar SDT components in shaping the net SDT response. Further experiments isolate the impact of changing propagation conditions through the zonal mean background atmosphere, non-linear wave-wave interactions, and the SSW-induced stratospheric ozone redistribution. Results indicate that between 80-97 km altitude in the northern hemisphere mid-to-high latitudes the net SDT response is driven by the solar SDT component, which itself is shaped by changing propagation conditions through the zonal mean background atmosphere and by non-linear wave-wave interactions. In addition, it is demonstrated that as a result of the rapidly varying solar SDT during the SSW the contribution of the lunar SDT to the total measured tidal field can be significantly overestimated.

From the wide range of possibilities, we propose an instrument capable of measuring annual changes in global Spectral Outgoing Radiation (SORa) from the entire Earth's surface between 200 nm and 1100 nm with a stability of 0.1 Wm-2 over a period of one solar cycle or beyond. Photomultiplier tubes (PMTs) as detectors provide data with a cadence of one second and high dynamic range. Based on Total Solar Radiation TSI(t) data with a stability of 0.01 Wm-2 per year, Spectral Solar Irradiance SSI(t) can be derived and normalized to ΣSSI(t)=TSI(t) for using the Sun as a referenced radiation source supported by solar modeling. Calibrated by SSI(t), a set of 12 spectrometers with 60 PMTs in total and 16 photometers simultaneously detect SOR(t). This database can also be provided to calibrate other space instruments to allow improved comparison of results. In previous missions in space, it has already been shown that the spectrometer design can detect both solar and terrestrial radiation with high dynamic range. The established measurement technique compensates for degradation through repeated calibration. The instrument also enables the determination of the global green Earth coverage and its annual changes by measuring chlorophyll absorption from 350 nm to 490 nm and 620 nm to 690 nm and green backscatter from 500 nm to 600 nm. Mapping the Earth will also make it possible to track annual local changes in green coverage and to assess the impact of different climate policies and climate engineering actions. Another aspect is the derivation of a correction parameter for the Earth Energy Imbalance derived from changes in green areas. Data evaluation can also include determining further parameters such as the Normalized Difference Vegetation Index, the Enhanced Vegetation Index, and the Global Leaf Area Index.

The Arctic has warmed much more than the global mean during past decades. The lapse-rate feedback (LRF) has been identified as large contributor to the Arctic amplification (AA) of climate change. This particular feedback arises from the vertically non-uniform warming of the troposphere, which in the Arctic emerges as strong near-surface, and muted free-tropospheric warming. Stable stratification and meridional energy transport are two characteristic processes that are evoked as causes for this vertical warming structure. Our aim is to constrain these governing processes by making use of detailed observations in combination with the large climate model ensemble of the 6th Coupled Model Intercomparison Project (CMIP6). We build on the result that CMIP6 models show a large scatter in Arctic LRF and AA, which are positively correlated for the historical period 1951–2014. Thereby, we present process-oriented constraints by linking characteristics of the current climate to historical climate simulations. In particular, we compare a large consortium of present-day observations to co-located model data from subsets with weak and strong simulated AA and Arctic LRF in the past. Our results firstly suggest that local Arctic processes mediating the lower thermodynamic structure of the atmosphere are more realistically depicted in climate models with weak Arctic LRF and AA (CMIP6/w) in the past. In particular, CMIP6/w models show stronger inversions at the end of the simulation period (2014) for boreal fall and winter, which is more consistent with the observations. This result is based on radiosonde observations from the year-long MOSAiC expedition in the central Arctic, together with long-term radio soundings at the Utqiaǵvik site in Alaska, USA, and dropsonde measurements from aircraft campaigns in the Fram Strait. Secondly, remote influences that can further mediate the warming structure in the free troposphere are more realistically represented by models with strong simulated Arctic LRF and AA (CMIP6/s) in the past. In particular, CMIP6/s models systemically simulate a stronger Arctic energy transport convergence in the present climate for boreal fall and winter, which is more consistent with reanalysis results. Locally, we find links between changes in transport pathways and vertical warming structures that favor a positive LRF in the CMIP6/s simulations. This hints to the mediating influence of advection on the Arctic LRF. We emphasise that one major attempt of this work is to give insights in different perspectives on the Arctic LRF. We present a variety of contributions from a large collaborative research consortium to ultimately find synergy among them in support of advancing our understanding of the Arctic LRF.

The teleconnection between the Quasi-Biennial Oscillation (QBO) and the Arctic stratospheric polar vortex, or the Holton-Tan (HT) relationship may change in a warmer climate or one with stratospheric aerosol intervention (SAI) as compared to present day climate (PDC). Our results from an Earth system model indicate that, under both global warming (based on RCP8.5 emission scenario) and SAI scenarios, the HT relationship weakens, although it is closer to PDC under SAI than under the RCP8.5 scenario. Such weakening of the HT relationship is more pronounced in early winter (Nov–Dec) compared to the mid-late winter period (Jan–Feb). While the high-latitude responses of temperature to the QBO anomalies are statistically significant under PDC, the responses are not statistically significant in the RCP8.5 and SAI scenarios. While the weakening of the HT relationship under RCP8.5 scenario is likely due to the weaker QBO wind amplitudes at the equator, another physical mechanism must be responsible for the weaker HT relationship under SAI scenario, since the amplitude of the QBO wind is comparable to the PDC. The strength of the polar vortex does not change under the RCP8.5 scenario compared to PDC, but it becomes stronger under SAI; we attribute the weakening of the HT relationship under SAI to such stronger polar vortex. In general, the changes in the HT relationship cannot be solely explained by changes to the critical line; the changes in the residual circulation (particularly due to the gravity wave contributions) are important too in explaining the changes in the HT relationship under RCP8.5 and SAI scenarios.

Meteor radars have become widely used instruments to study atmospheric dynamics, particularly in the 70 to 110 km altitude region. These systems have been proven to provide reliable and continuous measurements of horizontal winds in the mesosphere and lower thermosphere. Recently, there have been many attempts to utilize specular and/or transverse scatter meteor measurements to estimate vertical winds and vertical wind variability. In this study we investigate potential biases in vertical wind estimation that are intrinsic to the meteor radar observation geometry and scattering mechanism, and we introduce a mathematical debiasing process to mitigate them. This process makes use of a spatiotemporal Laplace filter, which is based on a generalized Tikhonov regularization. Vertical winds obtained from this retrieval algorithm are compared to UA-ICON model data. This comparison reveals good agreement in the statistical moments of the vertical velocity distributions. Furthermore, we present the first observational indications of a forward scatter wind bias. It appears to be caused by the scattering center's apparent motion along the meteor trajectory when the meteoric plasma column is drifted by the wind. The hypothesis is tested by a radiant mapping of two meteor showers. Finally, we introduce a new retrieval algorithm providing a physically and mathematically sound solution to derive vertical winds and wind variability from multistatic meteor radar networks such as the Nordic Meteor Radar Cluster (NORDIC) and the Chilean Observation Network De meteOr Radars (CONDOR). The new retrieval is called 3DVAR+DIV and includes additional diagnostics such as the horizontal divergence and relative vorticity to ensure a physically consistent solution for all 3D winds in spatially resolved domains. Based on this new algorithm we obtained vertical velocities in the range of w = ± 1–2 m s−1 for most of the analyzed data during 2 years of collection, which is consistent with the values reported from general circulation models (GCMs) for this timescale and spatial resolution.

The temporal and spatial variation of the upper transition height ( H T ) is studied using the measurements of the Incoherent Scatter Radar (ISR) at Arecibo Observatory ( 18.35 ° N , 66.75 ° W ) , Puerto Rico (magnetic latitude 30 ° N ). Observed ion compositions, electron densities, electron and ion temperatures are analysed to understand the diurnal and seasonal variations of H T at low and high solar activity covering solar cycles 23 and 24. The results show that during the low solar activity conditions, the mean H T is very low (around 750 km during the daytime and 500 km at night), while during high solar activity conditions, it reaches about 1600 km. Under very quiet geomagnetic conditions, the mean H T is about 700 km during day time and 500 km during night time. The obtained results are compared with the International Reference Ionosphere (IRI) model predictions using the TBT15 model ion composition option. During low and moderate solar activity conditions, the TBT15 model overestimates about 5–10% at night and 15–20% in the afternoon. Furthermore, from the current investigation using the IRI’s TBT15 and DY85 model ion composition options, it is found that the newer TBT15 model option is more appropriate to investigate H T variations as compared to the DY85 model option.

The Draconid meteor shower shows strong bursts of activity at irregular intervals, with nearly no activity in intervening years. Five outbursts of the Draconid meteor shower were observed with specular meteor radars in Canada and Europe between 1999 and 2018. The outbursts generally lasted between 6 and 8 h, and most were not fully visible at a single geographical site, emphasizing the need for observations at multiple longitudes for short-duration shower outbursts. There is at least a factor of two difference in the peak flux as measured on different radars; the initial trail radius effect is undercorrected for Draconid meteors, which are known to be fragile.

The instability of the Northern Hemisphere polar vortex is mainly caused by the breaking of planetary-scale (Rossby) waves (RWs). However, gravity waves (GWs) may also play an important role in polar vortex preconditioning before breakdown events. Moreover, GWs affect dynamics in the stratosphere by altering the upward propagation of RWs at short time scales and therefore indirectly influence polar vortex stability. Due to the coarse spatial resolution of global chemistry-climate models, current efforts in climate research rely on simulations where the majority of the GW spectrum is parameterized. In the present study, we apply a recently developed method for detecting strong orographic gravity wave (OGW) drag events in the lower stratosphere above the Himalayas. For this, we use a specified dynamics simulation of the chemistry-climate model CMAM (Canadian Middle Atmosphere Model) spanning the period 1979–2010. We show that strong OGW drag events above the Himalayas are associated with anomalously increased upward RW propagation in the stratosphere. This, in turn, is associated with an increase of the refractive index in the mid-latitude lower stratosphere, a region where the OGW drag dominates. Our results also illustrate that OGW strong events have the potential to alter ozone variability via changes of mixing in the surf zone and advection from lower latitudes. Altogether, we detail a preconditioning process of the polar vortex morphology by GWs above the Himalayas and how this relates to the proximity to polar vortex breakdown.

When orographic gravity waves (OGWs) break, they dissipate their momentum and energy and thereby influence the thermal and dynamical structure of the atmosphere. This OGW forcing mainly takes place in the middle atmosphere. It is zonally asymmetric and strongly intermittent. So-called OGW hotspot regions have been shown to exert a large impact on the total wave forcing, in particular in the lower stratosphere (LS). Motivated by this we investigate the asymmetrical distribution of the three-dimensional OGW drag (OGWD) for selected hotspot regions in the specified dynamics simulation of the chemistry-climate model CMAM (Canadian Middle Atmosphere Model) for the period 1979–2010. As an evaluation, we first compare zonal mean OGW fluxes and GW drag (GWD) of the model simulation with observations and reanalyses in the northern hemisphere. We find an overestimation of GW momentum fluxes and GWD in the model's LS, presumably attributable to the GW parameterizations which are tuned to correctly represent the dynamics of the southern hemisphere. In the following, we define three hotspot regions which are of particular interest for OGW studies, namely the Himalayas, the Rocky Mountains and East Asia. The GW drags in these hotspot regions emerge as strongly intermittent, a result that can also quantitatively be corroborated with observational studies. Moreover, a peak-detection algorithm is applied to capture the intermittent and zonally asymmetric character of OGWs breaking in the LS and to assess composites for the three hotspot regions. This shows that LS peak OGW events can have opposing effects on the upper stratosphere and mesosphere depending on the hotspot region. Our analysis constitutes a new method for studying the intermittency of OGWs, thereby facilitating a new possibility to assess the effect of particular OGW hotspot regions on middle atmospheric dynamics.

We used a nonlinear mechanistic global circulation model to analyze the migrating quarterdiurnal tide (QDT) in the middle atmosphere with focus on its possible forcing mechanisms: the absorption of solar radiation by ozone and water vapor, nonlinear tidal interactions, and gravity wave–tide interactions. We show a climatology of the QDT amplitudes, and we examine the contribution of the different forcing mechanisms to the QDT amplitude. To this end, we first extracted the QDT from the model tendency terms and then removed the respective QDT contribution from the different tendency terms. We find that the solar forcing mechanism is the most important one for the QDT; however, the nonlinear and gravity wave forcing mechanisms also play a role in autumn and winter, particularly at lower and middle latitudes in the mesosphere and lower thermosphere. Furthermore, destructive interference between the individual forcing mechanisms is observed. Therefore, tidal amplitudes become even larger in simulations with the nonlinear or gravity wave forcing mechanisms removed.

The global warming has been observed to be more severe in the Arctic compared to the rest of the world. This enhanced warming i.e. Arctic Amplification is not just the result of local feedback processes in the Arctic. The stratospheric pathways of Arctic-midlatitude linkages and large-scale dynamical processes can contribute to the Arctic Amplification. The polar stratospheric dynamics crucially depends on the atmospheric waves at all scales. The winter polar vortex can be disturbed by gravity waves in the middle atmosphere. To investigate the sensitivity of the polar vortex dynamics, large-scale dynamical processes, and the stratospheric pathways of the Arctic-midlatitude linkages to the modification of gravity wave drag, we conduct sensitivity experiments using the global atmospheric model ICON-NWP (ICOsahedral Nonhydrostatic Model for Numerical Weather Prediction). These sensitivity experiments are performed by imposing a repeated annual cycle of the year 1985 for sea surface temperatures and sea ice as lower boundary conditions and for greenhouse gas concentrations as external forcing. This year is selected as both El-Nino Southern Oscillation and Pacific decadal oscillation were in their neutral phase and no explosive volcanic eruption has occurred. Hence, lower boundary and external forcing conditions in this year can serve as a useful proxy for the multi-year mean condition and an estimate of its internal variability. We performed simulations where in the control simulation the sub-grid parameterization scheme for both orographic and non-orographic gravity wave drags are switched on. The other two experiments are identical to the control simulation except that either orographic or non-orographic gravity wave drags are switched off. Recently, deep learning has extraordinarily progressed our ability to recognize complex patterns in big datasets. Deep neural networks have shown great capabilities to capture the dynamical process of the atmosphere. Applying deep learning algorithms on experiments’ results, the impact of gravity wave drag modifications on large-scale mechanisms of the Arctic Amplification will be analyzed. Special emphasis will be put on the effects of gravity wave drag modifications on the polar vortex dynamics.

In recent years, multistatic specular meteor radars (SMRs) have been introduced to study the Mesosphere and Lower Thermosphere (MLT) dynamics with increasing spatial and temporal resolution. SMRs, compared to other ground-based observations, have the advantage of continuously measuring the region between 80 and 100 km independent of weather, season, or time of day. In this paper, frequency spectra of MLT horizontal winds are explored through observations from a campaign using the SIMONe (Spread-spectrum Interferometric Multistatic meteor radar Observing Network) approach conducted in northern Germany in 2018 (hereafter SIMONe 2018). The 7-day SIMONe 2018 comprised of fourteen multistatic SMR links and allows us to build a substantial database of specular meteor trail events, collecting more than one hundred thousand detections per day within a geographic area of $$\sim$$ ∼ 500 km $$\times$$ × 500 km. We have implemented two methods to obtain the frequency spectra of the horizontal wind components: (1) Mean Wind Estimation (MWE) and (2) Wind field Correlation Function Inversion (WCFI), which utilizes the mean and the covariances of the line of sight velocities, respectively. Monte Carlo simulations of a gravity wave spectral model were implemented to validate and compare both methods. The simulation analyses suggest that the WCFI helps us to capture the energy of smaller scale wind fluctuations than those capture with MWE. Characterization of the spectral slope of the horizontal wind at different MLT altitudes has been conducted on the SIMONe 2018, and it provides evidence that gravity waves with periods smaller than 7 h and greater than 2 h dominate with horizontal structures significantly larger than 500 km. In the future, these analyses can be extended to understand the significance of small-scale fluctuations in the MLT, which were not possible with conventional MWE methods. Graphical Abstract

Data assimilation schemes blend observational data, with limited coverage, with a short term forecast to produce an analysis, which is meant to be the best estimate of the current state of the atmosphere. Appropriately specifying observation error statistics is necessary to obtain an optimal analysis. Observation error can originate from instrument error as well as the error of representation. While representation error is most commonly associated with unresolved scales and processes, this term is often considered to include contributions from pre-processing or quality control and errors associated with the observation operator. With a focus on practical operational implementation, this chapter aims to define the components of observation error, discusses their sources and characteristics, and provides an overview of current methods for estimating observation error statistics. We highlight the implicit assumptions of these methods, as well as their shortcomings. We will detail current operational practice for diagnosing observation error and accounting for correlated observation error. Finally, we provide a practical methodology for using these diagnostics, as well as the associated innovation-based observation impact, to optimize the assimilation of meteor radar observations in the upper atmosphere.

The gravity waves play a crucial role in driving and shaping the middle atmospheric circulation. The Upper-Atmospheric extension of the ICOsahedral Non-hydrostatic (UA-ICON) general circulation model was recently developed with satisfying performances in both idealized test cases and climate simulations, however the sensitivity of the circulation to the parameterized orographic and non-orographic gravity wave drag remains largely unexplored. Using UA-ICON and ICON-NWP, the sensitivity of the dynamics and circulation to both orographic and non-orographic parameterized gravity waves effects are investigated. ICON-NWP stands for the numerical-weather prediction mode of the ICON model (see Zängl et al, 2015, QJRMetSoc), with a model top at about 80 km altitude. The UA-ICON mode differs from ICON-NWP in deep-atmosphere dynamics (instead of shallow-atmosphere dynamics) and upper-atmosphere physics parameterizations being switched on. In addition, the model top is at about 150 km. The sensitivity experiments involve employing repeated annual cycle sea surface temperatures, sea ice, and greenhouse gases under year 1988. This year is selected as both El-Nino southern oscillation and pacific decadal oscillation are in their neutral phase and no explosive volcano eruption has occurred and hence conditions in this year can serve as a useful proxy for the multi-year mean condition and an estimate of its internal variability. For both UA-ICON and ICON-NWP, we perform simulations where in the control (CTL) simulation both orographic and non-orographic gravity wave drags are switched on. The other two experiments are identical to the control simulation except that either orographic (OGWD-off) or b) non-orographic (NGWD-off) gravity wave drags are switched off. The analysis include comparisons between CTL and OGWD-off and NGWD-off simulations and include wave-mean flow interaction diagnostics (Eliassen-Palm flux and its divergence and refractive index of Rossby waves) and mass stream function of the Brewer-Dobson circulation. We also investigate the sudden stratospheric warming frequency and polar vortex morphology in order to understand whether a missing gravity wave forcing can further amplify or curtail the effects of future climate. We present our goal, method as well as first results and discuss possible further analysis.

In the lower ionospheric E region, shallow regions of high electron density are found, which are called sporadic E (ES) layers. ES layers consist of thin clouds of accumulated ions. They occur mainly at middle latitudes, and they are most frequently found during the summer season. ES are generally formed at heights between 90 and 120 km. At midlatitudes, their occurrence can be described through the wind shear theory. According to this theory, ES formation is due to interaction between the metallic ion concentration, the Earth’s magnetic field, and the vertical shear of the neutral wind. Here, we analyze ES occurrence rates (OR) obtained from ionospheric radio occultation measurements by the FORMOSAT-3/COSMIC constellation. To derive information on ES from RO, we use the Signal-to-Noise ratio (SNR) profiles of the GPS L1 phase measurements. If large SNR standard deviation values occur that are concentrated within a layer of less than 10 km thickness, we assume that the respective SNR profile disturbance is owing to an ES layer.Midlatitude ES are found to be mainly connected with a migrating diurnal and semidiurnal component. Especially at high latitudes of the southern hemisphere, nonmigrating components such as a diurnal westward wave 2 and a semidiurnal westward wave 1 are also visible. Near the equator, a strong diurnal eastward wavenumber 3 component and a semidiurnal eastward wavenumber 2 component are found in summer and autumn. Terdiurnal and quarterdiurnal components are weaker than the diurnal and semidiurnal ones. We discuss seasonal and global distributions of migrating and nonmigrating components, and their relation to neutral wind shear derived from ground-based observations and numerical modeling.

Radar observations from two SKiYMET radars at Collm (51°N, 13°E) and Kazan (56°N, 49°E) during 2016-2017 are used to investigate the longitudinal variability of the mesosphere/lower thermosphere (MLT) wind regime over western and eastern Europe. Both of the meteor radars have similar setups and apply the same analysis procedures to correctly compare MLT parameters and validate the simulated winds. The radar observations confirm the established seasonal variability of the wind distribution, but this distribution is not identical for the two stations. The results show good qualitative agreement with global circulations model predictions by the Middle and Upper Atmosphere Model (MUAM) and the Upper Atmosphere ICOsahedral Non-hydrostatic model (UA-ICON). The MUAM and UA-ICON models well reproduce the main dynamical features, namely the vertical and temporal distributions of the winds observed throughout the year. However, there are also some differences in the longitudinal wind variability of the models and radar observations. Numerical experiments with modified parameterization settings have also been carried out to study the response of the MLT wind circulation to the gravity waves originating from the lower atmosphere. The MUAM model results show that a decrease/increase in the gravity wave intensity at the lower atmosphere leads to an increase/decrease of the mesospheric zonal wind jet extension and the zonal wind reversal.

In the extratropical atmosphere, Rossby waves (RWs) and internal gravity waves (GWs) propagating from the troposphere mediate a coupling with the middle atmosphere by influencing the dynamics herein. In current generation chemistry-climate models (CCMs), GWs are usually smaller than the model resolution and the majority of their spectrum therefore must be parameterized. From observations, we know that GWs are intermittent and asymmetrically distributed around the globe, which holds to some extent also for the parameterized GW drag (GWD) (in particular for orographic GWD (oGWD)). The GW parameterizations in CCMs are usually tuned to mitigate biases in the zonal mean climatology of particular quantities, but the complex interaction of parameterized GWs with the large- scale circulation and resolved waves in the models remains to date poorly understood. This presentation will combine observational evidence, idealized modeling and dynamical analysis of a CCM output to study both the short-term and long-term model response to the oGWD. Our results demonstrate that the oGW-resolved dynamics interaction is a complex two-way process, with the most prominent oGWD impact being the alteration of propagation of planetary-scale Rossby waves on a time-scale of a few days. The conclusions give a novel perspective on the importance of oGWD for the stratospheric polar vortex and atmospheric transport studies outlining potential foci of future research.

The global warming has been observed to be more severe in the Arctic compared to the rest of the world. This enhanced warming i.e. Arctic Amplification is not just the result of local feedback processes in the Arctic. The stratospheric pathways of Arctic-midlatitude linkages and large-scale dynamical processes can contribute to the Arctic Amplification. The polar stratospheric dynamics crucially depends on the atmospheric waves at all scales. The winter polar vortex can be disturbed by gravity waves in the middle atmosphere. To investigate the sensitivity of the polar vortex dynamics, large-scale dynamical processes, and the stratospheric pathways of the Arctic-midlatitude linkages to the modification of gravity wave drag, we conduct sensitivity experiments using the global atmospheric model ICON-NWP (ICOsahedral Nonhydrostatic Model for Numerical Weather Prediction). These sensitivity experiments are performed by imposing a repeated annual cycle of the year 1986 for sea surface temperatures and sea ice as lower boundary conditions and for greenhouse gas concentrations as external forcing. This year is selected as both El-Nino Southern Oscillation and Pacific decadal oscillation were in their neutral phase and no explosive volcanic eruption has occurred. Hence, lower boundary and external forcing conditions in this year can serve as a useful proxy for the multi-year mean condition and an estimate of its internal variability. We performed simulations where in the control simulation the sub-grid parameterization scheme for both orographic and non-orographic gravity wave drags are switched on. The other two experiments are identical to the control simulation except that either orographic or non-orographic gravity wave drags are switched off. Recently, deep learning has extraordinarily progressed our ability to recognize complex patterns in big datasets. Deep neural networks have shown great capabilities to capture the dynamical process of the atmosphere. Applying deep learning algorithms on experiments’ results, the impact of gravity wave drag modifications on large-scale mechanisms of the Arctic Amplification will be analyzed. Special emphasis will be put on the effects of gravity wave drag modifications on the polar vortex dynamics.