Search

Parameterizations of subgrid-scale gravity waves (GWs) in atmospheric models commonly involve the description of the dissipation of GWs. Where they dissipate, GWs have an increased effect on the large-scale flow. Instabilities that trigger wave breaking are an important starting point for the route to dissipation. Possible destabilizing mechanisms are numerous, but the classical vertical static instability is still regarded as a key indicator for the disposition to wave breaking. In this work, we investigate how the horizontal variations associated with a GW could alter the criterion for static instability. To this end, we use an extension of the common parcel displacement method. This three-dimensional static stability analysis predicts a significantly larger range of instability than does the vertical static stability analysis. In this case, the Lindzen-type saturation adjustment to a state of marginal stability is perhaps a less suitable ansatz for the parameterization of the GW breaking. To develop a possible ansatz for the GW dissipation due to three-dimensional instability, we apply the methods of irreversible thermodynamics, which are embedded in the Gibbs formalism of dynamics. In this way, the parameterization does not only satisfy the second law of thermodynamics, but it can also be made consistent with the conservation of energy and further (non-)conservation principles. We develop the parameterization for a discrete spectrum of GW packets. Offline computations of GW drag and dissipative heating rates are performed for two vertical profiles of zonal wind and temperature for summer and winter conditions from CIRA data. The results are compared to benchmarks from the literature.

The ICON model with upper atmosphere extension (UA-ICON) running with the Atmosphere in the Earth System (AES) physics package suffers from a weak stratospheric polar vortex, though this package is intended for simulations of ICON on the seasonal to climatological time scales. For example, UA-ICON (AES) data from model free runs predict an unrealistically high number of SSWs per decade. However, due to the weak polar vortex, the simulated SSWs are of too weak intensity concerning the polar cap averaged temperature and mid-latitude zonal mean, zonal wind anomalies compared to observations. The climatology and variability of the Northern hemispheric (NH) stratospheric winter circulation widely improve when applying ICON and UA-ICON with the physics package for numerical weather prediction (NWP). However, whereas the NWP physics package is implemented technically for UA-ICON, it is not well-tuned for the upper mesosphere/ lower thermosphere (UMLT) region. This concerns, e.g. the temperatures in the summer mesopause regions at both hemispheres, which are much too warm, or the zonal mean wind structures in the lower thermosphere, which do not show the wind reversals from westerly to the easterly zonal wind in winter and vice versa in summer. These problems are caused mainly by an underrepresentation of the non-orographic gravity wave drag (NGWD), implemented with the parametrization based on Warner and McIntyre. The orographic gravity wave drag (OGWD) implemented with a parametrization based on Lott and Miller acts mainly in the stratosphere, but OGWD tuning also influences temperatures in the mesopause., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

In a companion paper, the Multiscale Gravity Wave Model (MS-GWaM) has been introduced and its application to a global model as a transient subgrid-scale parameterization has been described. This paper focuses on the examination of intermittency of gravity waves (GWs) modeled by MS-GWaM. To introduce the variability and intermittency in wave sources, convective GW sources are formulated, using diabatic heating diagnosed by the convection parameterization, and they are coupled to MS-GWaM in addition to a flow-independent source in the extratropics accounting for GWs due neither to convection nor to orography. The probability density function (PDF) and Gini index for GW pseudomomentum fluxes are assessed to investigate the intermittency. Both are similar to those from observations in the lower stratosphere. The intermittency of GWs over tropical convection is quite high and is found not to change much in the vertical direction. In the extratropics, where nonconvective GWs dominate, the intermittency is lower than that in the tropics in the stratosphere and comparable to that in the mesosphere, exhibiting a gradual increase with altitude. The PDFs in these latitudes seem to be close to the lognormal distributions. Effects of transient GW–mean-flow interactions on the simulated GW intermittency are assessed by performing additional simulations using the steady-state assumption in the GW parameterization. The intermittency of parameterized GWs over tropical convection is found to be overestimated by the assumption, whereas in the extratropics it is largely underrepresented. Explanation and discussion of these effects are included.

Current gravity wave (GW) parameterization (GWP) schemes are using the steady-state assumption, in which an instantaneous balance between GWs and mean flow is postulated, thereby neglecting transient, nondissipative interactions between the GW field and the resolved flow. These schemes rely exclusively on wave dissipation, by GW breaking or near critical layers, as a mechanism leading to forcing of the mean flow. In a transient GWP, without the steady-state assumption, nondissipative wave–mean-flow interactions are enabled as an additional mechanism. Idealized studies have shown that this is potentially important, and therefore the transient GWP Multiscale Gravity Wave Model (MS-GWaM) has been implemented into a state-of-the-art weather and climate model. In this implementation, MS-GWaM leads to a zonal-mean circulation that agrees well with observations and increases GW momentum-flux intermittency as compared with steady-state GWPs, bringing it into better agreement with superpressure balloon observations. Transient effects taken into account by MS-GWaM are shown to make a difference even on monthly time scales: in comparison with steady-state GWPs momentum fluxes in the lower stratosphere are increased and the amount of missing drag at Southern Hemispheric high latitudes is decreased to a modest but nonnegligible extent. An analysis of the contribution of different wavelengths to the GW signal in MS-GWaM suggests that small-scale GWs play an important role down to horizontal and vertical wavelengths of 50 km (or even smaller) and 200 m, respectively.

A virtual meteor radar system based on the upper-atmosphere extension of the high-resolution ICOsahedral Non-hydrostatic general circulation model is constructed to validate multistatic specular meteor radar (SMR) analyses. The virtual radar system examines the validity of mean winds and gradients estimation techniques used in multistatic SMRs. The study is motivated by unexpected mean values and tide-like features recently observed in the vertical velocities estimated from multistatic SMRs at different latitudes in the mesosphere and lower thermosphere. The proposed analysis confirms multistatic SMR systems' excellent capability to measure the horizontal mean wind components and gradient terms. It is also found that multistatic SMRs can estimate mean vertical winds if they have an amplitude greater than ±2 m/s. Due to the smoothing inherent to the model results, these results should be treated as lower bounds to the error incurred using real data. Hourly variability in vertical velocity estimates up to ±1–2 m/s in the observed vertical winds are due to contamination by small-scale horizontal structures in the horizontal winds.

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.

Specular meteor radars (SMRs) are a major ground-based instrument to study the mesosphere and the lower thermosphere (MLT) dynamics. The recently developed multi-static approach of SMRs allows maximising the number of measurements from different viewing angles, hence enabling the estimation of horizontal wind fields and their second-order statistics (power spectrum, momentum fluxes). We have installed the operational versions of these techniques in Germany, Peru and Argentina, called SIMONe (Spread-spectrum Interferometric Multistatic meteor radar Observing Network) systems. Here, we present a validation study of multi-static meteor radar analysis by using virtual radar systems on the upper-atmosphere extension of the ICOsahedral Non-hydrostatic (UA-ICON) general circulation model with a horizontal grid spacing of 5 km. This particular study is focusing on the estimates of gradients and vertical velocities with these multi-static systems.

The aim of the presented work is to improve the parameterization of subgrid-scale gravity wave (GW) effects on the resolved flow in atmospheric models in a large altitude range from the upper troposphere to the mesopause (~85km). State of the art GW parameterization schemes are using the steady-state approximation for the wave field and therefore assume an instantaneous GW propagation neglecting direct interactions between the GW field and the resolved flow within the whole altitude range mentioned above. As such, these schemes rely on dissipative processes - GW breaking and critical layer filtering - as the only mechanism to accelerate/decelerate the resolved flow. In contrast to this, by dropping the steady-state assumption a contribution to the mean-flow forcing emerges in the form of direct GW-mean-flow interactions. Several idealized studies show that, besides dissipative effects, direct GW-mean-flow interactions contribute to GW dynamics in an important extent (Bölöni et al., 2016, J. Atmos. Sci.}, 73, 4833-4852, Wilhelm et al., 2018, J. Atmos. Sci., 75, 2257-2280, Wei et al., 2019, J. Atmos. Sci., 76, 2715-2738). This motivates the implementation of a transient GW model (MS-GWaM: Multi Scale Gravity Wave Model) to UA-ICON, the upper atmosphere version of ICON (Borchert et al., 2019, Geosci. Model Dev., 12, 3541-3569) which does not rely on the steady-state assumption and thus includes direct GW-mean-flow interactions. As a reference and a representative of currently available GW parameterization schemes a steady-state version of MS-GWaM (ST-MS-GWaM) has been implemented to UA-ICON as well, which shares the treatment of all possible components (wave sources and wave saturation scheme) with the transient MS-GWaM scheme and differs from it "only" in the treatment of propagation, i.e. excluding direct GW-mean-flow interactions and thus transience. Both MS-GWaM and ST-MS-GWaM reproduce the observed wind and temperature climatology (e.g. URAP data: Swinbank, R. and D. A. Ortland, 2003, J. Geophys. Res., 108, D19, 4615) reasonably well, but the transient propagation makes a robust difference in the circulation in perpetual runs. The transient propagation in MS-GWaM substantially contributes to an increase of the GW intermittency in the whole altitude range, giving a better comparison with super-pressure balloon observations (e.g. Hertzog et al., 2012, J. Atmos. Sci., 69, 3433-3448), whereas the lack of transience prevents any occurrence of higher GW momentum flux values than the launch magnitude itself. This is explained by the fact that the direct GW-mean-flow interactions involve a highly transient evolution of the wave action and the vertical group velocity, which often leads to increased pseudo-momentum fluxes as compared to the launch magnitude.

The intermittency of gravity waves (GWs) is investigated using Multi-Scale Gravity Wave Model (MS-GWaM) implemented in the upper-atmosphere extension of ICON model. The intermittency of GWs is originated from that of wave sources but altered during propagation of the waves. Conventional GW parametrization (GWP), which diagnoses vertical profiles of GW properties under the steady-state assumption, can take into account the source intermittency if the GWP employs flow-dependent sources, while it cannot present the change of intermittency by transient evolutions of GWs. MS-GWaM is a prognostic model that explicitly solves the evolution of positions of waves (as well as their wavenumbers and amplitudes) in time and thus capable of describing the intermittency change. In order to include the source intermittency and variability, we couple the convective source, as diagnosed by subgrid-scale cumulus parametrization in ICON, to MS-GWaM, based on an analytic formulation of GW response to this source. In addition to this, a spatio-temporally uniform, persistent source is prescribed in the extratropics to take into account other non-orographic sources. Orographic sources are currently not used. The GW intermittency is measured by the Gini index, and is found to be quite high in the tropics, compared to that in the extratropics. In both regions, the index has similar values to those obtained from superpressure balloon observations reported in previous studies. A control experiment is performed using GWP based on the steady-state assumption, but coupled to the same wave sources, to assess the effects of transient modelling using MS-GWaM on the simulated intermittency. From comparison to the control experiment, the intermittency is found to increase largely for GWs from the uniform source but to decrease for convective GWs by the transient modelling.

The source mechanism of inertia–gravity waves (IGWs) observed in numerical simulations of the differentially heated rotating annulus experiment is investigated. The focus is on the wave generation from the balanced part of the flow, a process presumably contributing significantly to the atmospheric IGW field. Direct numerical simulations are performed for an atmosphere-like configuration of the annulus and possible regions of IGW activity are characterised by a Hilbert-transform algorithm. In addition, the flow is separated into a balanced and unbalanced part, assuming the limit of a small Rossby number, and the forcing of IGWs by the balanced part of the flow is derived rigorously. Tangent-linear simulations are then used to identify the part of the IGW signal that is rather due to radiation by the internal balanced flow than to boundary-layer instabilities at the side walls. An idealised fluid set-up without rigid horizontal boundaries is considered as well, to investigate the effect of the identified balanced forcing unmasked by boundary-layer effects. The direct simulations of the realistic and idealised fluid set-ups show a clear baroclinic-wave structure exhibiting a jet–front system similar to its atmospheric counterparts, superimposed by four distinct IGW packets. The subsequent tangent-linear analysis indicates that three wave packets are radiated from the internal flow and a fourth one is probably caused by boundary-layer instabilities. The forcing by the balanced part of the flow is found to play a significant role in the generation of IGWs, so it supplements boundary-layer instabilities as a key factor in the IGW emission in the differentially heated rotating annulus.

How the upper-atmosphere branch of the circulation contributes to and interacts with the circulation of the middle and lower atmosphere is a research area with many open questions. Inertia–gravity waves, for instance, have moved in the focus of research as they are suspected to be key features in driving and shaping the circulation. Numerical atmospheric models are an important pillar for this research. We use the ICOsahedral Non-hydrostatic (ICON) general circulation model, which is a joint development of the Max Planck Institute for Meteorology (MPI-M) and the German Weather Service (DWD), and provides, e.g., local mass conservation, a flexible grid nesting option, and a non-hydrostatic dynamical core formulated on an icosahedral–triangular grid. We extended ICON to the upper atmosphere and present here the two main components of this new configuration named UA-ICON: an extension of the dynamical core from shallow- to deep-atmosphere dynamics and the implementation of an upper-atmosphere physics package. A series of idealized test cases and climatological simulations is performed in order to evaluate the upper-atmosphere extension of ICON.

How the upper-atmosphere branch of the circulation contributes to and interacts with the circulation of the middle and lower atmosphere is a research area with many open questions. Inertia-gravity waves, for instance, have moved in the focus of research as they are suspected to be key features in driving and shaping the circulation. Numerical atmospheric models are an important pillar for this research. We use the ICOsahedral Non-hydrostatic (ICON) general circulation model, which is a joint development of the Max Planck Institute for Meteorology (MPI-M) and the German Weather Service (DWD), and provides, e.g., local mass conservation, a flexible grid nesting option and a non-hydrostatic dynamical core formulated on an icosahedral-triangular grid. We extended ICON to the upper atmosphere and present here the two main components of this new configuration named UA-ICON: an extension of the dynamical core from shallow- to deep-atmosphere dynamics, and the implementation of an upper-atmosphere physics package. A series of test cases and climatological simulations show that UA-ICON performs satisfactorily and is in good agreement with the observed global atmospheric circulation.

A finite-volume model of the classic differentially heated rotating annulus experiment is used to study the spontaneous emission of gravity waves (GWs) from jet stream imbalances, which may be an important source of these waves in the atmosphere and for which no satisfactory parameterisation exists. Experiments were performed using a classic laboratory configuration as well as using a much wider and shallower annulus with a much larger temperature difference between the inner and outer cylinder walls. The latter configuration is more atmosphere-like, in particular since the Brunt–Väisälä frequency is larger than the inertial frequency, resulting in more realistic GW dispersion properties. In both experiments, the model is initialised with a baroclinically unstable axisymmetric state established using a two-dimensional version of the code, and a low-azimuthal-mode baroclinic wave featuring a meandering jet is allowed to develop. Possible regions of GW activity are identified by the horizontal velocity divergence and a modal decomposition of the small-scale structures of the flow. Results indicate GW activity in both annulus configurations close to the inner cylinder wall and within the baroclinic wave. The former is attributable to boundary layer instabilities, while the latter possibly originates in part from spontaneous GW emission from the baroclinic wave.

Atmospheric new particle formation is a general phenomenon observed over coniferous forests. So far nucleation is either parameterised as a function of gaseous sulphuric acid concentration only, which is unable to explain the observed seasonality of nucleation events at different measurement sites, or as a function of sulphuric acid and organic molecules. Here we introduce different nucleation parameters based on the interaction of sulphuric acid and terpene oxidation products and elucidate the individual importance. They include basic trace gas and meteorological measurements such as ozone and water vapour concentrations, temperature (for terpene emission) and UV B radiation as a proxy for OH radical formation. We apply these new parameters to field studies conducted at conducted at Finnish and German measurement sites and compare these to nucleation observations on a daily and annual scale. General agreement was found, although the specific compounds responsible for the nucleation process remain speculative. This can be interpreted as follows: During cooler seasons the emission of biogenic terpenes and the OH availability limits the new particle formation while towards warmer seasons the ratio of ozone and water vapour concentration seems to dominate the general behaviour. Therefore, organics seem to support ambient nucleation besides sulphuric acid or an OH-related compound. Using these nucleation parameters to extrapolate the current conditions to prognosed future concentrations of ozone, water vapour and organic concentrations leads to a significant potential increase in the nucleation event number.

The differentially heated rotating annulus is a widely studied tabletop-size laboratory model of the general mid-latitude atmospheric circulation. The two most relevant factors of cyclogenesis, namely rotation and meridional temperature gradient are quite well captured in this simple arrangement. The radial temperature difference in the cylindrical tank and its rotation rate can be set so that the isothermal surfaces in the bulk tilt, leading to the formation of baroclinic waves. The signatures of these waves at the free water surface have been analyzed via infrared thermography in a wide range of rotation rates (keeping the radial temperature difference constant) and under different initial conditions. In parallel to the laboratory experiments, five groups of the MetStr\"om collaboration have conducted numerical simulations in the same parameter regime using different approaches and solvers, and applying different initial conditions and perturbations. The experimentally and numerically obtained baroclinic wave patterns have been evaluated and compared in terms of their dominant wave modes, spatio-temporal variance properties and drift rates. Thus certain ``benchmarks'' have been created that can later be used as test cases for atmospheric numerical model validation.