MS-GWaM: A Three-Dimensional Transient Gravity Wave Parametrization for Atmospheric Models

Georg S. Voelker aInstitute for Atmosphere and Environment, Goethe University Frankfurt, Frankfurt, Germany

Search for other papers by Georg S. Voelker in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0003-3658-8515
,
Gergely Bölöni bDeutscher Wetterdienst, Offenbach am Main, Germany

Search for other papers by Gergely Bölöni in
Current site
Google Scholar
PubMed
Close
,
Young-Ha Kim aInstitute for Atmosphere and Environment, Goethe University Frankfurt, Frankfurt, Germany

Search for other papers by Young-Ha Kim in
Current site
Google Scholar
PubMed
Close
,
Günther Zängl bDeutscher Wetterdienst, Offenbach am Main, Germany

Search for other papers by Günther Zängl in
Current site
Google Scholar
PubMed
Close
, and
Ulrich Achatz aInstitute for Atmosphere and Environment, Goethe University Frankfurt, Frankfurt, Germany

Search for other papers by Ulrich Achatz in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Parameterizations for internal gravity waves in atmospheric models are traditionally subject to a number of simplifications. Most notably, they rely on both neglecting wave propagation and advection in the horizontal direction (single-column assumption) and an instantaneous balance in the vertical direction (steady-state assumption). While these simplifications are well justified to cover some essential dynamic effects and keep the computational effort small, it has been shown that both mechanisms are potentially significant. In particular, the recently introduced Multiscale Gravity Wave Model (MS-GWaM) successfully applied ray-tracing methods in a novel type of transient but columnar internal gravity wave parameterization (MS-GWaM-1D). We extend this concept to a three-dimensional version of the parameterization (MS-GWaM-3D) to simulate subgrid-scale nonorographic internal gravity waves. The resulting global wave model—implemented into the weather forecast and climate code Icosahedral Nonhydrostatic (ICON)—contains three-dimensional transient propagation with accurate flux calculations, a latitude-dependent background source, and convectively generated waves. MS-GWaM-3D helps reproduce expected temperature and wind patterns in the mesopause region in the climatological zonal mean state and thus proves a viable internal gravity wave (IGW) parameterization. Analyzing the global wave action budget, we find that horizontal wave propagation is as important as vertical wave propagation. The corresponding wave refraction includes previously missing but well-known effects such as wave refraction into the polar jet streams. On a global scale, three-dimensional wave refraction leads to a horizontal flow-dependent redistribution of waves such that the structures of the zonal mean wave drag and consequently the zonal mean winds are modified.

© 2024 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Georg S. Voelker, voelker@iau.uni-frankfurt.de

Abstract

Parameterizations for internal gravity waves in atmospheric models are traditionally subject to a number of simplifications. Most notably, they rely on both neglecting wave propagation and advection in the horizontal direction (single-column assumption) and an instantaneous balance in the vertical direction (steady-state assumption). While these simplifications are well justified to cover some essential dynamic effects and keep the computational effort small, it has been shown that both mechanisms are potentially significant. In particular, the recently introduced Multiscale Gravity Wave Model (MS-GWaM) successfully applied ray-tracing methods in a novel type of transient but columnar internal gravity wave parameterization (MS-GWaM-1D). We extend this concept to a three-dimensional version of the parameterization (MS-GWaM-3D) to simulate subgrid-scale nonorographic internal gravity waves. The resulting global wave model—implemented into the weather forecast and climate code Icosahedral Nonhydrostatic (ICON)—contains three-dimensional transient propagation with accurate flux calculations, a latitude-dependent background source, and convectively generated waves. MS-GWaM-3D helps reproduce expected temperature and wind patterns in the mesopause region in the climatological zonal mean state and thus proves a viable internal gravity wave (IGW) parameterization. Analyzing the global wave action budget, we find that horizontal wave propagation is as important as vertical wave propagation. The corresponding wave refraction includes previously missing but well-known effects such as wave refraction into the polar jet streams. On a global scale, three-dimensional wave refraction leads to a horizontal flow-dependent redistribution of waves such that the structures of the zonal mean wave drag and consequently the zonal mean winds are modified.

© 2024 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Georg S. Voelker, voelker@iau.uni-frankfurt.de

1. Introduction

Internal gravity waves (IGWs) play a significant role in distributing energy and momentum throughout the atmosphere. Being generated through the perturbation of balanced flow states in stratified environments (e.g., convection, flow over topography, flow instabilities, and wave turbulence), they may propagate over large distances and deposit their energy, generate turbulence, and modify the mean flow dynamics in general far away from their source (Fritts and Alexander 2003; Alexander et al. 2010; Sutherland 2010; Nappo 2013; Williams et al. 2017; Sutherland et al. 2019; Achatz 2022). While propagating, they may exchange energy among themselves, transiently interact with mean flows, or influence the transport of chemical species. Even though their generation regions are mostly located in the troposphere, their effects are strongest in the middle atmosphere (Lindzen 1981; Kim et al. 2003; Sigmond and Scinocca 2010). However, these major effects may couple to and thus impact, for instance, tropospheric weather patterns or climate conditions (Scaife et al. 2005, 2012; Sigmond and Scinocca 2010).

In state-of-the-art general circulation models (GCMs) or numerical weather prediction (NWP) models, gravity waves are only partially resolved and thus need parameterization (e.g., Kim et al. 2003; Holt et al. 2016). Even when entering the global kilometer-resolving regime, their effects are not entirely represented (Kruse et al. 2022; Polichtchouk et al. 2022, 2023). IGW parameterizations commonly rely on WKBJ theories (Bretherton 1966; Grimshaw 1975; Achatz et al. 2017), incorporating some major simplifications (Lindzen 1981; Medvedev and Klaassen 1995; Warner and McIntyre 1996; Hines 1997b,a; Lott and Miller 1997; Alexander and Dunkerton 1999; Scinocca 2003; Orr et al. 2010; Lott and Guez 2013). Most notably, there are three assumptions that are usually made and shall be considered in this work. First, local horizontal homogeneity is assumed in the so-called single-column approximation, i.e., horizontal gradients of the resolved flow are neglected, which leads horizontal wavenumbers to be invariant during propagation. Furthermore, responses of the resolved flow due to the horizontal finiteness of an IGW field are neglected. Second, a steady-state assumption is adopted such that IGWs instantly propagate through the atmosphere, leading to a neglect of all transient propagation effects. Third, it is commonly assumed that the resolved flow is (approximately) in hydrostatic and geostrophic balance, implying a reduced formula for the resolved flow response to IGWs. In the extratropics, this can lead to significant modifications in the IGW forcing when the flow is imbalanced.

Recent investigations suggested that the mentioned effects can have important impacts on modeled flows (Sato et al. 2009; Bölöni et al. 2016; Ehard et al. 2017; Wei et al. 2019). In particular, Bölöni et al. (2016) showed that including the time dependence of gravity wave propagation significantly improves the estimated wave drag in idealized simulations. Sato et al. (2009) found that three-dimensional propagation pathways can impact the structure of the stratosphere and mesosphere using wave-resolving general circulation models. Similarly, Ehard et al. (2017) highlight the importance of horizontal gravity wave propagation for the dynamics of the Antarctic polar night jet by analyzing various observations and reanalyses. Finally, Wei et al. (2019) showed that describing the wave momentum fluxes with a direct formulation can improve the simulated wave impact on the resolved flow. Several studies have successfully relaxed some of the simplifications in model parameterizations (e.g., Muraschko et al. 2015; Wilhelm et al. 2018; Quinn et al. 2020). In particular, Bölöni et al. (2021) and Kim et al. (2021) presented a novel Lagrangian IGW parameterization, the Multiscale Gravity Wave Model (MS-GWaM). It is built on a weakly nonlinear WKBJ theory, including transient wave–mean flow interactions. The model has been implemented in a single-column mode into a state-of-the-art weather forecast and climate code. It will therefore be referred to as MS-GWaM-1D throughout this manuscript. Here, we present MS-GWaM-3D, an extension of MS-GWaM-1D, which models the full 3D transient propagation of IGWs and their forcing of the general (cf. balanced) resolved flow. While both the columnar approximation and assumptions on the resolved flow are relaxed in MS-GWaM-3D, we shall focus on the effects of horizontal propagation in this study. In particular, the modified refraction behavior and resolved-flow response in the statistical mean are investigated. Indeed, we find that many aspects show an increased realism of the wave dynamics. However, our results also underline the challenges left for future studies. Ultimately, we aim to improve the representation of the parameterized IGW processes while keeping simulations numerically efficient.

This paper is structured as follows. First, a brief recapitulation of the underlying WKBJ theory for transient, three-dimensional IGWs (section 2) is presented. It is then followed by the description of the employed ray-tracing techniques (section 3). Last, results from simulations with MS-GWaM-3D (section 4) are visualized and discussed. The manuscript then closes with some concluding remarks on the achievements and challenges of ray-tracing parameterization (section 5).

2. Nonlinear, three-dimensional, and transient internal gravity waves

Although generally three-dimensional and transient, internal gravity waves are commonly parameterized using both single-column and steady-state approximations. Here, we attempt to relax both these assumptions and build a three-dimensional and transient parameterization. An underlying multiscale theory has been presented by Achatz et al. (2010, 2017, 2023), and thus, we restrict ourselves to summarizing the results here for brevity.

a. Wave evolution

We follow the analysis of Bretherton (1966), Grimshaw (1974), Hasha et al. (2008), Achatz et al. (2017), Achatz (2022), and Achatz et al. (2023) and summarize the Bretherton–Grimshaw modulation equations for a quasi-monochromatic wave train on a sphere. In particular, the wave and mean flow dynamics are derived by means of a WKBJ expansion of the compressible Euler equations. We extend the existing theories by the inclusion of the atmospheric density scale height to improve realism in the nonhydrostatic regime. Most notably, the scale height correction plays an important role in the course of vertical wave reflection, where vertical wavenumbers become small with respect to the inverse scale height Γ. We then arrive at the well-known dispersion relation,
ω(r,t)=Ω(r,t,k)=kU±f2(kr2+Γ2)+N2(kλ2+kϕ2)K2,
with the position vector, r = λeλ ϕeϕ + rer, and the wavevector, k = kλeλ + kϕeϕ + krer. Here, λ, ϕ, and r are the standard geographical spherical coordinates and eλ, eϕ, and er are the corresponding unit vectors. Note that the mean-flow velocity U(r, t) = Ueλ + Veϕ, the Coriolis frequency f = f(r), and the buoyancy frequency N = N(r, t) are functions of space and time as indicated. For convenience, we have defined the total wavevector squared as K2 = kk + Γ2, with the scale height correction Γ = Γ(r, t). Finally, we would like the reader to note that the theory is written in spherical coordinates. The corresponding eikonal equations and the group velocities cg for an individual wave component then follow the relations:
(t+cgr)k=k˙=rΩ,cg=kΩ.
Here, we denote the gradients with respect to position and wavenumber by ∇r and ∇k, respectively. Having the application to a ray-tracing scheme in mind, we indicate by a dot the derivative along characteristics, so-called rays, with their positions denoted by r(t) so that
r˙=cg.
Rays follow the local group velocity and hence the propagation of wave energy. In particular, the prognostic equations for the zonal, meridional, and radial wavenumbers (kλ, kϕ, kr) in spherical coordinates read
k˙λ=1rcosϕ[kλU+|kh|22ω^K2(λN2N2f2K2λΓ2)]kλrcgr+kλtanϕrcgϕ,
k˙ϕ=1r[kϕU+|kh|22ω^K2(ϕN2+kr2+Γ2|kh|2ϕf2N2f2K2ϕΓ2)]kϕrcgrkλtanϕrcgλ,
kr˙=[krU+|kh|22ω^K2(rN2N2f2K2rΓ2)]+kλrcgλ+kϕrcgϕ.
Note that the metric terms ensure propagation on great circles in the case of a uniform atmosphere at rest, as laid out in detail by Hasha et al. (2008).1 Note that the representation in spherical coordinates comes with the drawback of a pole problem in the above relations when propagating very close to the pole. Practically, singularities due to division by zero, however, do not pose a serious problem as they only occur in close vicinity of the poles. To avoid potential problems, we do not consider wave generation where |ϕ| > 85° (see section 3c).
Finally, for the locally monochromatic field, its amplitude is governed by the wave action equation,
0=tA+r(cgA)+S,
for the wave action density A in physical space, where S denotes the sources and sinks of wave action, such as the wave saturation as described below. Wave–wave interactions between IGWs and interactions between IGWs and the (subgrid-scale) geostrophic modes (GM) are neglected, and thus, the wave action conservation is valid for any overlapping wave fields, provided their amplitudes are sufficiently weak. Numerous processes, e.g., wave refraction, can lead to rays crossing in physical space at so-called caustics, so that the assumption of local monochromaticity can break down in the course of the integration of the eikonal equations above. To avoid numerical instabilities due to this issue, we introduce the spectral wave action density N=N(r,k,t)=jAj(r,t)δ[kk(j)(r,t)] in phase space (cf. Muraschko et al. 2015; Achatz 2022; and references therein), where the index j indicates any member of a possibly infinitely large and infinitely dense set of locally monochromatic fields that are being superposed. The resulting phase-space wave action conservation then reads
0=tN+r(cgN)+k(k˙N)+S=(t+cgr+k˙k)N+S.
Note that due to the nondivergence of the six-dimensional phase space velocity,
0=rcg+kk˙,
the flux and the transport formulations in Eq. (8) are equivalent. Again, S denotes both sources and the dissipation of phase-space wave action. Thus, the phase-space wave action density N is conserved along wave propagation paths up to wave dissipation and generation. Moreover, the nondivergence [Eq. (9)] of the phase space velocity also implies that any phase space volume following its rays conserves its volume content. As a consequence, rays cannot cross in phase space and caustics do not occur. Summarizing, Eqs. (3)(6) and Eq. (8) form the basis for the forward ray tracing of gravity waves.

b. Wave impact on the mean flow

To account for the impact of the wave perturbations on the mean flow, we note that the IGW perturbation fluxes enter the mean flow tendencies as
tU|IGW=1ρ¯(ρ¯vu)+fger×ub,
tθ|IGW=uθ,
where er is the radial unit vector, g is Earth’s gravity, and ρ¯ and θ¯ are the reference density and potential temperature (Achatz et al. 2017, 2023). Moreover, v′, b′, and θ′ denote the IGW velocity, buoyancy, and potential temperature perturbations and u′ denotes the horizontal component of v′. The brackets 〈⋅〉 represent the phase average of the wave perturbations. Using the dispersion and polarization relations, they may be expressed as
ρ¯vu=[ω^2c^gkhN+f2(er×c^g)(er×khN)ω^2f2]d3k,
uθ=(θ¯gρ¯c^grfN2ω^2f2er×khN)d3k,
fger×ub=fθ¯er×uθ,
where we denote the horizontal wavevector as kh = kλeλ + kϕeϕ and where c^g=cgU is the intrinsic group velocity. Note that the elastic term (f/g)er×ub and the potential temperature flux divergence ∇ ⋅ 〈uθ′〉 are negligible where the Coriolis frequency f is small with respect to the intrinsic wave frequency |f||ω^|. Thus, the two terms are weak for nonhydrostatic waves in the mid- to high-frequency range. In addition, the momentum flux divergence may then be approximated with the pseudomomentum flux divergence, (ρ¯vu)c^gkhNd3k. Interestingly, the corresponding flux, often termed the Eliassen–Palm flux, also determines the evolution of the mean flow potential vorticity and the geostrophically and hydrostatically balanced mean flow even without any spectral restriction. It is thus often used in IGW parameterizations (e.g., Warner and McIntyre 2001; Orr et al. 2010; Bölöni et al. 2021). For a brief comparison of the two approaches, see appendix A. Recently, Wei et al. (2019) showed that the direct flux formulation more accurately represents the impact of an internal wave packet on the mean flow. Hence, we use the full set of fluxes shown in Eqs. (10)(14).

Combined with the eikonal equations [Eqs. (4)(6)], wave propagation [Eq. (3)], and wave action conservation [Eq. (8)], these relations form a closed prognostic system for transient wave propagation. Note that the equations are energy-conserving wherever both wave sources and dissipation are zero (not shown explicitly).

3. Gravity wave ray tracing

As pointed out in section 2, the subgrid-scale gravity wave dynamics are natively described on the rays of the corresponding gravity wave propagation. To include both the transience and the horizontal wave propagation, one may employ various techniques to obtain a numerical solution. Ray tracing has gained significant attention recently as an efficient and accurate method for simulating wave fields (Marks and Eckermann 1995; Muraschko et al. 2015; Amemiya and Sato 2016; Voelker et al. 2021; Bölöni et al. 2021). Here, we build on the implementations of Bölöni et al. (2021) and Kim et al. (2021) and extend the approach by the horizontal wave propagation and the full flux calculation [Eqs. (10) and (11)]. While we present a broad overview in this section, we would like to refer the reader to appendix A for additional details.

a. Representation of rays as phase-space ray volumes

One of the major challenges for ray tracing is the problem of caustics and the subsequent ambiguity of the notion of wave amplitudes and wave energy of a local spectrum. As laid out above, it is possible to reformulate the wave action conservation as a transport equation in the six-dimensional phase space. The nondivergence of the phase space velocity [Eq. (9)] implies that any phase space volume following its rays conserves its volume content, i.e.,
0=dt(dxdydzdkxdkydkz)=dt(r2cosϕdλdϕdrdkλdkϕdkr),
where (x, y, z) are the Cartesian coordinates and (kx, ky, kz) are the wave vector components in these coordinates. In Cartesian coordinates, one has even
0=xcgx+kxk˙x,
0=ycgy+kyk˙y,
0=zcgz+kzk˙z,
so that
0=dt(dxdkx)=dt(dydky)=dt(dzdkz).
Such a simple splitting relies on the independence of the spatial coordinates and does not exist in spherical coordinates. The approach taken in MS-GWaM, for the time being, is to divide the phase-space content with nonzero wave action density N into finite-size small ray volumes, with extent Δ(λ, ϕ, r, kλ, kϕ, kr). Within the spanned ray volume, the wave action density N is assumed to be constant. The ray-volume content is then propagated according to Eqs. (3)(6) following a central carrier ray. The change in physical extent Δ(λ, ϕ, r) is determined using rays at the faces of the ray volumes carrying identical wave properties (kλ, kϕ, kr) as in the central carrier ray but being exposed to deviations in the background fields. In particular, the ray-volume extent is approximated as locally Cartesian by
Δ(x,y,z)=(rcosϕΔλ,rΔϕ,Δr),
Δ(kx,ky,kz)=Δ(kλ,kϕ,kr),
where the new Δ(kx, ky, kz) is determined by Eq. (19) and then inverted using Eq. (21) to finally obtain the new Δ(kλ, kϕ, kr). Note that this procedure, albeit converging for infinitesimally small ray volumes, implies the following simplification. In a uniform atmosphere at rest, where the group velocity differences between the opposing ray-volume faces are exactly zero, the ray-volume extent in spherical coordinates Δ(λ, ϕ, r) is exactly conserved. This implies, however, that any ray volume propagating poleward must shrink in the tangent-linear extent Δx and expand in spectral width Δkλ (and vice versa for southward propagation). Although this effect may be small in regions sufficiently separated from the pole, we plan to apply an advanced method to remove this potential problem in the next version of MS-GWaM.

b. Coupling of Lagrangian particles to the Eulerian mean flow

The coupling between the large-scale flow and the wave perturbations necessitates linking the Lagrangian wave representation to the Eulerian model grid. For the evaluation of the modulation equations, we linearly interpolate all necessary mean fields and their gradients to the ray-volume center or faces using a first-order Taylor expansion. In particular, the group velocities are calculated at the physical boundaries of the ray volume (see section a in appendix B). This serves a dual purpose: The mean of the group velocities on opposing ray-volume faces may be used to advance the location of the central carrier ray, and the difference between them serves as an estimate for the compression or inflation of the ray volume in physical space. The phase-space extent is then determined as described above.

To account for the impact of waves on the mean flow, it is necessary to differentiate between vertical and horizontal components, considering possibly irregular horizontal grid structures. The vertical gradient of the wave flux, i.e., r(ρ¯wu) in Eq. (10), is computed following Muraschko et al. (2015) and Bölöni et al. (2016, 2021). In particular, the vertical flux associated with a gridcell boundary is integrated over all ray volumes contained in the vertically staggered cell that covers the desired cell boundary (for more details, see Fig. 1a in Bölöni et al. 2021). The gradients are then computed from the centered difference of the vertical cell boundary values. For the horizontal fluxes, the gridcell-averaged values (rather than those at the boundaries) are calculated to avoid problems with complex grid geometries. We integrate all ray volumes whose carrier rays are contained in the chosen grid cell’s volume. To mend errors due to ray volumes overlapping multiple adjacent cells, the obtained flux fields are horizontally smoothed (section d in appendix B). The resulting fields and their divergences [Eqs. (10) and (11)] then contribute to the tendencies provided by the parameterization.

c. Sources and sinks of internal gravity waves

Similar to Bölöni et al. (2021) and Kim et al. (2021), the present implementation of MS-GWaM utilizes two sources of gravity waves as a lower boundary condition (Fig. 1). In particular, a background source in the mid- to high latitudes, accounting, e.g., for IGWs emitted from jets and fronts, as well as a convective source, is employed. The former prescribes at a given launch level a seasonally varying pseudomomentum flux with a spectral shape following the Desaubie spectrum (cf. Orr et al. 2010). At the solstice, it has a latitudinal profile with winter hemispheric fluxes of Fw = 8 mPa, no fluxes for latitudes ϕ ∈ [20°S, 20°N], and summer hemispheric fluxes of Fs = 2.4 mPa. No background waves are launched at latitudes higher than |ϕ| > 85°. The different regimes are connected through smooth transitions. The profile oscillates in time with a sinusoidal yearly cycle. Summarizing, we employ a latitudinal profile α(φ) and a total absolute momentum flux at the lower boundary during the northern summer solstice Fsol,
α(φ)={min[1,max(0,1φ805)],ifφ801,ifφ[50,80)12[1cos(πφ2030)],ifφ[20,50)0,else,
Fsol(φ)=α(φ)Fs+α(φ)Fw,
where φ is the latitude (°N). The absolute momentum flux at the lower boundary condition then reads
Fam=12[Fsol(φ)(1+cosτ)+Fsol(φ)(1cosτ)],
with the nondimensional time difference, τ = 2π(ttsol)/1y, relative to the time of the northern summer solstice tsol. These background waves are launched at every time step at a height corresponding to a pressure of 300 hPa (black, dotted line in Figs. 1a,c) with a momentum flux which is equally divided into four horizontal directions as done by Bölöni et al. (2021). In particular, we construct a ghost layer at the launch altitude within which ray volumes are allowed to propagate vertically only. At every launch interval, that is, every 60 s (equivalent to the model time step), the fraction of the ray volume which has propagated out of the ghost layer is separated and activated. Only then, it will propagate freely according to the modulation equations. Ray volumes which are reflected within the ghost layer are discarded. More technical information on the launch process can be found in the work of Bölöni et al. (2021). This background wave energy source is directionally homogeneous, with a purely latitudinal profile on the 300-hPa pressure surface. Additionally, convectively generated waves are considered with a launch spectrum as described by Kim et al. (2021). In their formulation, the areal fraction of the convective latent heat release ϵq inversely controls the total entropy forcing and thus the launch momentum flux (Kim et al. 2021). Here, we increase the area fraction to εq = 7%, resulting in a reduced-gravity wave flux. For the convective source, we couple MS-GWaM to the convection parameterization and utilize the convection strength and depth calculated within. The ray volumes are then launched at the cloud top and under consideration of the background wind. As an example of the resulting lower boundary conditions, we show the absolute vertical momentum fluxes, defined by
Fam=ρ¯(wu2+wυ2)1/2,
through the surface at a constant height of 16 km for June 1991 (Fig. 1). The thick line presents the median of the zonally averaged fluxes during the month, and the shading around the median shows the range between the 10th and 90th percentiles, highlighting the variability in the fluxes. Even though this height is rather close to the launch altitudes, we observe that the horizontal wave propagation changes the statistics of the vertical fluxes of background waves (Figs. 1a,c). In particular, we note that the gravity wave transience leads to a reduction in the median absolute momentum fluxes and a strong variability with fluxes regularly exceeding the launch flux. This can be explained by a spatially and temporally varying mean flow through which the waves propagate from the launch level to the depicted 16-km altitude. For instance, the waves may encounter reductions and increases in their vertical group velocity, critical layers, or vertical reflection. In addition, the horizontal propagation leads to both an increased dissipation below 16 km and a horizontal redistribution, effectively changing the wave flux statistics. Comparing Figs. 1a and 1b, we thus find a reduced median flux for MS-GWaM-3D and background waves far outside the mid- to high-latitude generation regions. In contrast to the background source, the convectively generated waves show rather similar statistics at this altitude between MS-GWaM-3D and MS-GWaM-1D (Figs. 1b,d). Naturally, their statistics differ from the background waves as convective events are spatially and temporally rather sparse but can lead to high momentum fluxes. Thus, the spatiotemporal median is very low, while the peak momentum fluxes are comparable to the background waves.
Fig. 1.
Fig. 1.

Zonally averaged absolute vertical fluxes of horizontal momentum at 16-km height for (a),(b) MS-GWaM-1D and (c),(d) MS-GWaM-3D for June 1991. The solid lines represent the 50th percentile (median), and the surfaces show the ranges from the 10th to 90th percentiles. Results are separated into (a),(c) background waves and (b),(d) convectively generated waves. The black dotted lines represent the launched absolute momentum fluxes at 300 hPa during the summer solstice.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Finally, we would like to add that the launch fluxes Fw and Fs, as well as the area fraction ϵq, were tuned to achieve a realistic zonal mean mesopause. In particular, both the wind reversal and the summer hemispheric cold pole are expected to occur at approximately 80 km (see discussion of mean model states in section 4a). To avoid high computational costs due to optimization algorithms, short runs were considered, integrating the model for 2 weeks in both seasons. The model results, calculated for a physical range of the three tunable parameters, were then compared, and the best fit was chosen. The resulting directional momentum fluxes compare in order of magnitude with satellite observations at altitudes 30–40 km but differ in horizontal distribution (Jiang et al. 2006; Ern et al. 2018; Hindley et al. 2020; Polichtchouk et al. 2022; direct comparisons not shown). The latter is, however, expected given that MS-GWaM incorporates a flow-independent background source.

The two sources of wave energy are accompanied by a saturation parameterization following Bölöni et al. (2016, 2021), which is based on the idea of Lindzen (1981) but in a spectrally integrated manner. In particular, the threshold for static instability is formulated for the wave action as follows:
R3kr2|kh|2ω^K2Nd3k=P<ρ¯2.
Note that for a quasi-monochromatic wave, this relation reduces to the well-known threshold, |Bw|2<N4/kr2, with Bw being the amplitude of IGW buoyancy. The concept of the integrated density perturbation P is then used in a discrete sense and rewritten as the sum over all physically overlapping ray volumes j at the location rb such that
P=jPj=jkrj2|khj|2ω^jKj2Ajαj.
Here, Aj=NjΔkλΔkϕΔkr and αj are the ray volume’s total wave action and the spatial fraction that the jth ray volume takes of the volume of the cell it is contained in, respectively. Thus, rather than each individual spectral component, we consider a spectrally integrated saturation. The reduction in the individual wave action components Aj to their saturation value Aj is then achieved through setting
NjNj*=Nj[12κKj2Δt],
with the turbulent diffusivity coefficient
κ=max(0,Pρ¯22ΔtjPjKj2).
Saturation is thus triggered where the instability threshold [Eq. (26)] is violated by the constructive interference of locally overlapping internal waves. All wave amplitudes, or equivalently wave action densities, are then set to the corresponding saturation values. Note that Eq. (29) does not imply that the turbulent diffusivity diverges as Δt → ∞. Because at the beginning of each time step Pρ¯/2, the numerator of the ratio is at best growing linearly with Δt, so that κ always ends up finite. This spectral saturation implementation was shown to compare well to wave-resolving simulations in idealized cases by Bölöni et al. (2016). Furthermore, it is planned to couple the saturation scheme to the turbulence parameterization of the underlying atmospheric model, Icosahedral Nonhydrostatic (ICON). An in-depth analysis of the breaking processes and the impact of the spectral saturation on a global scale is hence left to a future study.

4. Atmospheric simulations with UA-ICON and MS-GWaM-3D

Building on Bölöni et al. (2021) and Kim et al. (2021), who introduced MS-GWaM-1D, we choose to work with the ICON model with its upper-atmosphere extension (Zängl et al. 2015; Borchert et al. 2019). Here, we use version 2.6.5-nwp1b with a horizontal resolution of approximately 160 km (model grid R2B04) and the physics packages for the NWP and the upper atmosphere. Note that MS-GWaM replaces the nonorographic gravity wave drag but leaves the operational orographic wave parameterization intact (in particular, Lott and Miller 1997). Through the coupling to the mean wind, the latter thus indirectly interacts with MS-GWaM. The setup has a model top of 150 km with a vertical grid extent of a few 10 m in the boundary layer, 700–1500 m in the stratosphere, and a maximum of approximately 4 km in the lower thermosphere. A sponge layer acts above an altitude of 110 km, which is why we restrict our analysis to altitudes below 100 km. The model is initialized with IFS analysis data below 60 km and with the climatological thermodynamic state at rest above. It is then spun up for a month to exclude adjustment effects in our analysis. These runs are repeated to simulate June and December for the years 1991 through 1998. Additionally, we deploy runs for June and December 1991 with extended diagnostics. In particular, the free parameters controlling the launch fluxes are the convective area fraction for the convectively generated waves and the summer and winter launch amplitudes of the background waves (see section 3c). In general, stronger gravity wave launch fluxes shift the wave saturation to lower altitudes, leading to a lower wind reversal at the mesopause. Comparisons of simulations with the horizontal wind model (HWM) (HWM2014; Drob et al. 2015) and the empirical temperature model NRLMSIS2.1 (MSIS; Emmert et al. 2022) then give a good estimate for the performance of the chosen parameters. The HWM2014 and MSIS are chosen based on their wealth of incorporated data in the middle and upper atmospheres in the time range of our simulations. In particular, they incorporate reanalysis, ground-based datasets, and satellite-based observations (for more details, see Drob et al. 2015; Emmert et al. 2021, 2022). The reference climatologies are then generated for exactly the simulated months. With the help of short simulations spanning a large parameter space, we identified an optimum with an area fraction of εq = 7%, a winter launch flux of Fw = 8 mPa, and a summer launch flux of Fs = 2.4 mPa (for more details, see section 3c).

a. Zonal mean flow

To test the general performance of the gravity wave parameterization, we analyze the zonal-mean zonal wind and temperatures as composites over eight simulations of June and December (Figs. 2 and 3). To highlight the general effect of the parameterizations and visualize the general model performance, we accompany the data from simulations without nonorographic gravity wave parameterization. Naturally, the runs including a gravity wave drag parameterization mend model biases and thus differ significantly from the runs without. The comparison is, therefore, to be understood as a reference for the general behavior of ICON in the present setup and in how far model biases stem from the lack of a gravity wave drag parameterization or may be rooted elsewhere. Additionally, we show the reference climatologies HWM2014 and MSIS. Note that, as described above, both versions of MS-GWaM were tuned such that the mesopause height at z ≈ 80 km agrees with HWM2014. In general, MS-GWaM-1D and MS-GWaM-3D produce physical zonal-mean zonal winds with the expected wind reversals in both the summer and winter hemispheres (Figs. 2b–d,f–h). As one may expect, runs with the same model setup but without any nonorographic parameterization are subject to strong jets and entirely lack the mesopause reversal (Figs. 2a,e).

Fig. 2.
Fig. 2.

Zonal mean zonal winds from composites of runs for (a)–(d) June and (e)–(h) December. The composites correspond to averages for the years 1991 through 1998. Runs from (b),(f) MS-GWaM-1D and (c),(g) MS-GWaM-3D are compared to (a),(e) simulations without any nonorographic IGW parameterization and (d),(h) the HWM2014 climatologic model.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Fig. 3.
Fig. 3.

As in Fig. 2, but for zonally averaged temperatures. For comparison, we show the (d),(h) climatological model MSIS2.1.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Notable differences to the HWM2014 can be seen in the summer hemispheric jets, in the tropical stratosphere, and in the Antarctic winter jet. Comparing the two MS-GWaM results with the run without any nonorographic gravity wave drag parameterization (Figs. 2a,e), we find that the tropical stratospheric jet at z ≈ 50 km height is a robust feature of the current model version, which is present in all three runs. The aforementioned tuning runs indicated that there is little to no influence on the structure by the tuning parameters (not shown). We thus conclude that it is generally independent of the gravity wave parameterization. The weak northern summer jet, however, seems to be related to the transient parameterization in general, as it is present in runs with both MS-GWaM-1D and MS-GWaM-3D (Figs. 2b,c). While we note a similar bias in the work of Bölöni et al. (2021), the runs without any nonorographic gravity wave parameterization do not suffer from this bias (Fig. 2d). Finally, we observe that the slant of the Antarctic winter jet toward the equator is not entirely reproduced by MS-GWaM-3D (Fig. 2c). This is due to a modified zonal and meridional gravity wave drags in the Southern Hemisphere above 60-km altitude in MS-GWaM-3D as opposed to MS-GWaM-1D as will be discussed in section 4c. Possible reasons for the modified wave drag might be improper spectral characteristics of the current background source, errors near the pole due to the representation of the ray-volume extent in spherical coordinate systems (cf. section 3a), or even codependencies of the columnar gravity wave parameterization with other parameterizations such as the employed orographic wave drag parameterization. Moreover, other components of the model, such as the radiation scheme or the modified ozone climatology, are tuned to produce realistic temperature distributions in combination with the standard orographic wave drag. In general, the confidence in the present implementation is high, as MS-GWaM-1D was validated by Bölöni et al. (2021) and Kim et al. (2021) and MS-GWaM-3D was validated against the results of Wei et al. (2019) using an idealized cubic geometry implemented in ICON (not shown). However, a more detailed analysis may be needed to identify the reasons behind the phenomenon and mend the setup accordingly.

Similar to the zonally averaged zonal winds, we find that the zonally averaged temperatures and the corresponding MSIS fields agree well (Fig. 3). In particular, the MS-GWaM reproduces the cold summer pole at the mesopause and the warm winter pole associated with the stratopause (Figs. 3b,c,f,g). These features are naturally not present when the nonorographic gravity wave drag parameterization is switched off (Figs. 3a,e).

In general, the zonal mean model state generally seems to be moderately affected by enabling three-dimensional wave propagation in MS-GWaM. While the runs with the transient parameterizations generally exert some bias in the northern summer hemisphere, the Antarctic winter jet is too strong at altitudes around 50–80 km for MS-GWaM-3D. Runs for December do not contain these systematic biases. This resembles the findings of Bölöni et al. (2021) where a weak dependency of the zonal mean wind on the transience of the gravity wave parameterization was detected. Albeit somewhat counterintuitive, it may not be a surprise, as the spatial and temporal averages represent the quasi-steady state of the model setups. These, in turn, are tuned to approximate a similar zonal mean model state. It was, however, also shown that both the transience and the localization of the gravity waves can have significant impacts on the structure of the horizontal mean flow (e.g., Sato et al. 2009; Šácha et al. 2016; Ehard et al. 2017). For instance, the mean-state results may change if a more realistic spatiotemporal distribution of the wave source spectrum was employed. While the investigation of realistic source distribution is left for future work, in the following sections, we focus on the nature of three-dimensional propagation of parameterized IGWs, examining wave action fluxes and budgets and horizontal structures of IGW momentum fluxes and mean wave drags.

b. Global three-dimensional IGW distribution

Gravity waves are generally three-dimensional phenomena that are present in large parts of the atmosphere. While wave-resolving simulations can provide insights into global wave activity, they typically span short model times due to the limitations of both computational resources and required storage. Additionally, the separation of IGWs from other dynamics and their characterization are generally nontrivial tasks. Being based on a multiscale WKBJ theory (Achatz et al. 2017), MS-GWaM simulates the global behavior of nonorographic gravity waves without the need to explicitly resolve their fast varying phases. As a result, it can be used to predict global distributions of phase-space wave action densities and derived quantities such as momentum fluxes and the resulting wave drag. As an example, we analyze the global wave action budget and augment it with the zonal wave drag to illustrate the resulting effect on the mean flow.

Consider the phase-space integrated wave action conservation [Eq. (8)], averaged in time,
ΔAΔt=hcghA¯+rcgrA¯+S¯,
where we define (cgh,cgr)A=(cgh,cgr)Nd3k, and the overbars denote the temporal average over the time interval Δt. Here, the three right-hand side terms correspond to the temporally averaged wave action tendency contributions from the horizontal wave propagation, the vertical wave propagation, and the wave saturation. Moreover, the left-hand term is small for sufficiently long averaging intervals. For a sufficiently small left-hand side of Eq. (30) while preserving the planetary wave structure, we choose Δt = 10 days and show the predictions of the right-hand side contributions by MS-GWaM-1D and MS-GWaM-3D at altitudes of approximately 25 and 60 km for 10–20 June 1991. While the 25-km level shows the stark differences between the 1D and 3D schemes, the wave dynamics at 60-km altitude highlight the interaction with the Antarctic polar night jet in more detail.

As one may expect, MS-GWaM-1D predicts a balance between vertical wave action propagation and dissipation at either altitude (Figs. 4c,e and 5c,e). Thus, wave breaking and strong mean flow impacts can only occur where IGWs were previously able to vertically propagate through the wind shear of the underlying air column. Negative values for the change in wave action due to vertical propagation are associated with a finite temporal averaging interval Δt. For instance, propagating mean flow structures such as planetary waves may lead to changes in the wave action budget, which adjust only on the time scale of the wave propagation. Although not linearly related, the strongest zonal wave drags occur where the wave dissipation is high. At a rather low altitude of 25 km, this results in a rather uniform distribution of wave dissipation and drag in the Southern Hemisphere, accompanied with intermittent convective events in the Northern Hemisphere (Figs. 4e–g). Running the same case with MS-GWaM-3D, we find that both the contributions due to the horizontal and vertical propagation are of similar magnitude at an altitude of 25 km, such that the resulting balance includes all terms (Figs. 4b,d,f). Near horizontally sheared wind structures, such as the edges of the Antarctic winter jet, the two contributions from the horizontal and vertical propagation form inverted dipoles associated with wave refraction phenomena. In the case of the mentioned Antarctic winter jet, as we will show below, this corresponds to the often observed wave refraction into the jet (e.g., Ehard et al. 2017; Hindley et al. 2020). Interestingly, the wave refraction follows the structure of the resolved planetary waves and ultimately leads to a shift of lower-level wave action dissipation toward regions of strong horizontal shear, as, for instance, found between the summer and winter jets (Fig. 4f). Despite the different launch amplitudes (see section 3c), all contributions have similar magnitudes on both hemispheres. Correspondingly, the horizontal propagation significantly modifies the zonal wave drag at 25 km. While the Southern Hemisphere drag is much reduced, the Northern Hemisphere exhibits comparatively strong drags near planetary wave structures. Interestingly, the wave action budgets change structurally with altitude. At z ≈ 60 km, both the budgets of MS-GWaM-1D and MS-GWaM-3D are mostly controlled by vertical wave action propagation and dissipation (Figs. 5c–f). In contrast to the lower altitude, horizontal wave action propagation is much reduced. However, horizontal propagation through lower altitudes shifts the action dissipation and wave drag near the southern winter jet toward the equator (Figs. 5e–h). We suspect this shift to cause the missing wave drag at higher altitudes near the pole, leading to a biased jet structure in the zonal mean. These results suggest an important influence of the horizontal propagation near the Antarctic winter jet, and thus, we analyze the zonal mean drag structures in the following.

Fig. 4.
Fig. 4.

Wave action tendencies for (left) MS-GWaM-1D and (right) MS-GWaM-3D due to (a),(b) horizontal wave propagation, (c),(d) vertical wave propagation, and (e),(f) wave dissipation. The last row depicts the corresponding zonal wave drag. All tendencies are 10-day averages for 10–20 Jun 1991 at an altitude of z ≈ 25 km. Additionally, contours of the zonal wind in steps of 10 m s−1 are shown.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Fig. 5.
Fig. 5.

As in Fig. 4, but for z ≈ 60 km.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

c. Wave propagation and drag in the vicinity of horizontally sheared jet structures

Expanding on the horizontal refraction behavior found in the horizontal maps (Figs. 4 and 5), we show zonal averages of wave action fluxes (Fig. 6) and the corresponding wave drag (Fig. 7) temporally averaged over free runs for June and December 1991 through 1998. By construction, MS-GWaM-1D only allows for vertical propagation (Figs. 6a,c). The corresponding wave action tendency is thus associated with the vertical wave action flux convergence. The consequently generated wave drag accelerates the zonal wind structures, with maxima occurring at the mesopause (Figs. 7a,c). Note that the maximum wave action flux convergences are not equivalent to the maximum wave drag due to the scaling with the wavevector (cf. section 2). MS-GWaM-3D shows strong horizontal propagation in regions of strong horizontal and vertical shear, such as jet edges (Figs. 6b,d). Notably, in the Southern Hemisphere winter, the wave action flux vector points northward near the South Pole, i.e., into the southern winter jet (Fig. 6b). Consequently, the polar wave action fluxes are reduced at altitudes near the mesopause and the wave drag maximum is weaker south of 75°S (Figs. 6b and 7b). Moreover, the equatorial regions are strongly affected by the horizontal wave propagation (Figs. 7b,d). Here, the wind shear is generally strong, but the mean wave action fluxes are weak due to the intermittency of the convective wave source.

Fig. 6.
Fig. 6.

Wave action changes due to wave propagation in (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D for composites of (a),(b) June and (c),(d) December. The field is overlaid with vectors of wave action fluxes. They are scaled with respect to the plot dimensions and units, and the horizontal direction is enhanced by a factor of 3 for visualization. The arrow transparency represents the wave action flux strength. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments, with dashed and solid lines for negatives and nonnegatives, respectively. See Fig. 7 for the corresponding wave drag.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Fig. 7.
Fig. 7.

Zonally averaged zonal gravity wave drag due to the parameterized, nonorographic IGWs. The results are averaged over composites of (a),(b) June and (c),(d) December 1991–98 and plotted for both (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

To highlight the dynamics near the Antarctic winter jet, we additionally show the zonal and meridional wave drags for both MS-GWaM-1D and MS-GWaM-3D near the jet, averaged for all June (Fig. 8). Here, the reduced zonal wave drag for MS-GWaM-3D, south of 75°S, becomes particularly visible (Figs. 8a,b). Moreover, the maximum zonal wave drag is shifted southward due to the horizontal propagation. Stronger differences can be seen for the generally weaker meridional component (Figs. 8c,d). In particular, it changes sign and exerts a southward drag on the jet when enabling horizontal propagation. Note that this behavior is only observed above the jet. Globally, the meridional drag accelerates the flow from the winter toward the summer hemisphere for both models (not shown). It remains unclear, however, whether the modeled meridional drag is a consequence of the too strong winter jet over the pole or vice versa. For instance, the ray-volume discretization, the interaction with the orographic wave drag parameterization, or the structure of the modeled temperatures may be a cause of the biased wind structures (also compare section 4a and Fig. 2c). Further and more targeted investigation may be needed to bring clarity about this effect.

Fig. 8.
Fig. 8.

Zonally averaged zonal gravity wave drag due to the parameterized, nonorographic IGWs. The results are averaged over composites of (a),(b) June and (c),(d) December 1991–98 and plotted for both (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

Differences in refraction behavior can also be seen in the summer hemispheres. However, the model biases of the northern summer hemispheric jet structures (cf. section 4a) complicate the interpretation of the wave action budgets in the corresponding regions.

One may now ask the question how realistic the modeled meridional wave propagation near the Antarctic winter jet is and how it corresponds to satellite observations. To provide some insight, we show the meridional momentum fluxes [cf. Eq. (12)] from the composites of June for the years 1991 through 1998 at altitudes ≈ 40 km (Fig. 9). Note that the meridional fluxes of MS-GWaM-1D must be interpreted keeping in mind the assumed local horizontal homogeneity. That is, albeit the wave parameters suggest horizontal momentum fluxes, the horizontal propagation (including the wave advection) is set to zero in the columnar setup. The resulting synthetic meridional fluxes seem to be unrelated to the underlying zonal jet structure (Fig. 9a). As waves are launched homogeneously in all directions at the lower boundary (section 3c), the resulting meridional fluxes are a product of the partial wave filtering due to local vertical gradients of both the horizontal wind and the buoyancy frequency. In contrast, MS-GWaM-3D shows meridional convergence of wave momentum fluxes into the polar jet (Fig. 9b). In particular, the meridional fluxes are northward over the Antarctic continent and southward further north. The change in sign, i.e., the convergence zone, roughly coincides with the zonal wind maximum. This convergence effect of internal gravity waves into the polar night jets is well documented from both observations and general circulation models (e.g., Dunkerton and Butchart 1984; Sato et al. 2012; Ehard et al. 2017; Hindley et al. 2020; Gupta et al. 2021). Notably, the structure and amplitudes of the meridional momentum fluxes from MS-GWaM-3D simulations resemble the observations of Hindley et al. (2020; compare their Fig. 3d). Naturally, MS-GWaM cannot reproduce the strong momentum fluxes originating from orographic waves over the Antarctic Peninsula, the Southern Andes, or the islands of the Southern Ocean. The otherwise fairly close fit, however, gives us high confidence in the correctness of MS-GWaM-3D.

Fig. 9.
Fig. 9.

Meridional momentum fluxes in (a) MS-GWaM-1D and (b) MS-GWaM-3D for composites of June from 1991 through 1998 at an altitude of z ≈ 40 km overlaid by the corresponding zonal winds (contours; m s−1). The maps show the typically observed refraction of IGWs into the southern polar night jet and modified wind patterns.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

5. Discussion of achievements and challenges

Although IGWs are an essential part of atmospheric dynamics, their parameterization in general circulation models is typically subject to various simplifications. In particular, the single-column assumption, the steady-state approximation, and the balanced mean flow assumption have usually been made. Here, we present a MS-GWaM-3D, which aims at parameterizing IGWs through Lagrangian ray tracing without the three mentioned assumptions. In general, we find that including 3D transience significantly impacts global IGW propagation patterns and thus horizontal IGW distributions. Correspondingly, the associated wave drag and mean flows are modified.

We compare two distinct situations: For MS-GWaM-1D, where the columnar approximation is applied, the vertical propagation is balanced by wave breaking in the wave action equation (Figs. 4c,e and 5c,e). Only waves that did not encounter critical layers according to their wave properties may propagate to higher altitudes and eventually break. In contrast, MS-GWaM-3D allows for modulated and spatially unconstrained propagation and enables waves to be refracted around wind structures. As a consequence, both horizontal and vertical propagation balance with the wave dissipation (Figs. 4b,d,f and 5b,d,f). The equal order of magnitude in the contributions of horizontal and vertical wave propagation to the wave action budget emphasizes the importance of including horizontal propagation. Moreover, it questions the validity of using columnar methods when, e.g., studying IGW distributions. Both methods do, however, perform as expected in reproducing the cold summer pole and the warm winter pole at altitudes ∼85 and ∼60 km, respectively, in the climatological zonal mean (Fig. 3). The corresponding wind reversals and the mesopause altitudes are reasonably predicted (Fig. 2). We are thus confident that MS-GWaM-3D, as MS-GWaM-1D before it, covers some major effects of IGWs on the mean flow dynamics, rendering MS-GWaM-3D a viable IGW parameterization. The 3D wave propagation is also found to have significant impacts on the variability of the zonal mean flow, such as the quasi-biennial oscillation (QBO), as presented by Kim et al. (2023).

Including the horizontal propagation does, however, also introduce some important differences in the simulated climatology. In particular, the Southern Hemisphere winter jet becomes stronger near the pole at altitudes 60–80 km for MS-GWaM-3D (Fig. 2c). The reason behind this change is, however, not obvious. One possible explanation is the northward refraction of gravity waves near the Antarctic winter jet in MS-GWaM-3D as compared to MS-GWaM-1D. Both the 3D wave action tendencies (dipole structures in Figs. 4d,f) and the zonally averaged wave action fluxes (vectors in Figs. 6a,b) suggest that the 3D modulation causes IGWs to propagate northward and thus relate to a weaker gravity wave drag at high altitudes and similar latitudes (Figs. 8a,b). Additionally, the meridional wave drag near the jet is oppositely directed, that is southward, for MS-GWaM-3D (Figs. 8c,d). It remains unclear, however, whether the wave refraction and ultimately the wave drag cause the shifts of the jets as the two structures a nonlinearly coupled. Possible other causes are the numerical simplifications made in the setup of the ray tracer in the spherical coordinate system or codependencies of the gravity wave parameterization with other parameterizations, for instance, the orographic gravity wave drag. Finally, the jets are sensitive to the structure of the modeled temperature distribution, which, in turn, is dependent on, for instance, the ozone climatology and the radiation scheme. A thorough investigation will be needed, and we hope to answer these questions in an envisioned follow-up study. Finally, we find that the meridional wave propagation in MS-GWaM-3D around the Antarctic polar night jet resembles the observations of Hindley et al. (2020) both in amplitudes and in structure, up to orographic wave sources such as the southern Andes (Fig. 9). This gives us high confidence into the implementations as the convergence into the polar night jet is well documented in a number of studies (e.g., Dunkerton and Butchart 1984; Sato et al. 2012; Ehard et al. 2017; Hindley et al. 2020; Gupta et al. 2021).

There is a long list of effects that may be analyzed in more detail from this starting point. Notably, the role of IGWs in sudden stratospheric warmings (e.g., Baldwin et al. 2021), the accuracy of the final warming date (e.g., De la Cámara et al. 2016; Eichinger et al. 2023), the missing wave drag at 60°S (e.g., Holt et al. 2023), etc., are of interest and are planned for future simulations and analyses. Current investigations include the impact of 3D transient waves on the quasi-biennial oscillation and wave intermittency. Moreover, there are efforts into increasing the realism of the included IGW sources related to jet frontal systems and flow over topography. Finally, code optimization may reduce the computational cost, which amounts to a runtime factor ∼30 with respect to runs without any nonorographic wave drag parameterization (cf. section c in appendix B). With these challenges in mind, we ultimately aim for the application of MS-GWaM-3D in climate simulations and possibly numerical weather predictions.

Furthermore, there are numerous potential improvements that may be added to MS-GWaM in the future, albeit not considered here. Most importantly, MS-GWaM-3D—for the time being—does not include a flow-dependent description of the emission of IGWs from jets and fronts (Charron and Manzini 2002; Richter et al. 2010; De la Cámara and Lott 2015) or flow over topography (Palmer et al. 1986; Bacmeister et al. 1994; Lott and Miller 1997; Xie et al. 2020, 2021; van Niekerk and Vosper 2021; van Niekerk et al. 2023; Eichinger et al. 2023). For the latter reason, MS-GWaM needs augmentation with a subgrid-scale orography parameterization. Moreover, it is not clear how important modulated triadic resonant interactions among IGWs or between subgrid-scale geostrophic modes and IGWs are for atmospheric dynamics, which could, however, be investigated through ray-tracing techniques (cf. Kafiabad et al. 2019; Voelker et al. 2021).

1

Revisiting past studies, we recognize that Ribstein et al. (2015) and Ribstein and Achatz (2016) omit some of the important metric terms so that their results should be interpreted with care.

Acknowledgments.

UA thanks the German Research Foundation (DFG) for partial support through the research unit “Multiscale Dynamics of Gravity Waves” (MS-GWaves; Grants AC 71/8-2, AC 71/9-2, and AC 71/12-2) and CRC 301 “TPChange” (Project 428312742 and Projects B06 “Impact of small-scale dynamics on UTLS transport and mixing” and B07 “Impact of cirrus clouds on tropopause structure”). YHK and UA thank the German Federal Ministry of Education and Research (BMBF) for partial support through the program Role of the Middle Atmosphere in Climate (ROMIC II: QUBICC) and Grant 01LG1905B. UA and GSV thank the German Research Foundation (DFG) for partial support through the CRC 181 “Energy transfers in Atmosphere and Ocean” (Project 274762653 and Projects W01 “Gravity-wave parameterization for the atmosphere” and S02 “Improved Parameterizations and Numerics in Climate Models”). UA is furthermore grateful for the support by Eric and Wendy Schmidt through the Schmidt Futures VESRI “DataWave” project. This work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under Project bb1097.

Data availability statement.

All data used in this publication may be available upon request to the corresponding author.

APPENDIX A

The Gravity Wave Forcing of the Mean Flow and Its Relation to the EP Flux Formulation

As laid out in section 2, MS-GWaM-3D utilizes the perturbation flux formulation as introduced by Wei et al. (2019). The interested reader may note that the resulting mean flow impact differs from the commonly used pseudomomentum flux convergence related to the Eliassen–Palm flux (Eliassen and Palm 1961; Andrews and McIntyre 1976). To highlight the relation to the latter, we would first like to remind the reader about the underlying assumption of the multiscale theory applied here. First, the theory is built on the compressible Euler equations such that nonlinearities not only arise from the advection operators but also from the pressure gradient terms. As MS-GWaM acts as a gravity wave parameterization for the underlying Eulerian model, phase averages are constructed in a Eulerian sense, accordingly. In terms of scale separation, it is assumed that horizontal scales are much larger than vertical scales. Furthermore, the Brunt–Väisälä frequency is assumed to be much larger than the Coriolis frequency. These assumptions then lead to the regime of small Rossby numbers with a leading-order mean flow in hydrostatic and geostrophic balance. Finally, all wave scales are assumed to be much smaller than the synoptic mean flow scales.

Under these conditions, the impact of the internal gravity waves on the mean momentum comprises the momentum flux convergence and an elastic term arising from variations in density and pressure in the compressible atmosphere. In addition, the mean entropy is forced by a wave potential temperature flux convergence. More specifically, the mean flow impact can then be written as
tU|IGW=1ρ¯(ρ¯vu)+fger×ub,
tθ|IGW=uθ,
with the wave fluxes expressed as
ρ¯vu=[ω^2c^gkhN+f2(er×c^g)(er×khN)ω^2f2]d3k,
uθ=[θ¯gρ¯c^grfN2ω^2f2er×khN]d3k,
fger×ub=fθ¯er×uθ.
It is now interesting to additionally consider the quasigeostrophic potential vorticity evolution in the presence of internal gravity waves. The corresponding prognostic equation reads
(t+Ur)P=err×[1ρ¯(rHrG0)],
where P=r2ψ+ρ¯1z(ρ¯f2/N2zψ) is the quasigeostrophic potential vorticity, denoted in terms of the meanflow streamfunction ψ. Note that the latter relates to the leading-order synoptic-scale wind through U = er × ∇rψ and that within geostrophic theory also the thermodynamic fields can be obtained to leading order from the streamfunction. Hence, as long as the mean flow is within the synoptic scaling regime, it is controlled completely by Eq. (A6). Furthermore, we have denoted the horizontal pseudomomentum fluxes H and G defined by the spectral integral over the phase-space pseudomomentum fluxes, G=kc^gNd3k and H=lc^gNd3k. For a full derivation, we refer the reader to the work of Achatz et al. (2017) and Achatz (2022), and a recent summary was given by Achatz et al. (2023). Of interest is now that, if the mean flow is in geostrophic and hydrostatic balance, Eq. (A6) can also be derived from a modified system of mean-flow equations where one replaces
ρ¯vuc^gkhNd3k,
ub=gθ¯uθ0,
and deletes the elastic term. This is the simplification that is most commonly applied in gravity wave parameterizations: The momentum flux convergence in the momentum equation is replaced by the convergence of pseudomomentum flux, the elastic term is neglected, and the entropy flux is set to zero. Under the conditions stated, the mean flow then follows to leading order of the same dynamics as predicted by (A1) and (A2).

Since the zonal means of zonal wind and temperature are to a good approximation in geostrophic and hydrostatic balance, we therefore do not find a statistically significant difference in those fields. However, Wei et al. (2019) showed that the local transient response of the Eulerian phase averaged wind in the presence of a finite amplitude gravity wave packet becomes more accurate when considering the fluxes as denoted in Eqs. (A1) and (A2). We thus implement this formulation but leave the in-depth analysis of the local wind response to an envisioned follow-up study.

APPENDIX B

Details of the Numerical Implementation

To simulate transient internal gravity waves, MS-GWaM uses a spectral extension of WKBJ theory (Bretherton 1966; Grimshaw 1975; Achatz et al. 2017; Achatz 2022; and references therein). A Lagrangian approach is employed to solve the equations numerically (Muraschko et al. 2015; Bölöni et al. 2016, 2021; Kim et al. 2021). Therein, the six-dimensional phase space, spanned by position and wavenumber, is split into small finite-size ray volumes, propagating according to the eikonal equations [Eqs. (3)(6)]. They keep their regular shape of rectangles in (λ, ϕ, r, kλ, kϕ, kr), but their extent in all these six directions is allowed to vary in a volume-preserving manner. The wave action density carried by each ray volume is constant, unless wave breaking leads to the onset of turbulence and hence a decrease in N. These discrete Lagrangian volumes must, however, be connected to the local Eulerian model state to incorporate direct nonlinear interactions between the mean flow and the gravity waves (cf. sections 2 and 3). This coupling poses a number of computational challenges, which are solved as described below.

a. Three-dimensional field interpolation

When calculating the propagation of a finite-size ray volume in an inhomogeneous medium, it is important to consider that even under the assumption that all wavenumbers are constant throughout the volume, the group velocities can differ significantly within the same ray volume. In particular, the scale height correction Γ, the Brunt–Väisälä frequency N, and the mean winds U are the functions of space and differ between the four faces forming with the meridional and zonal boundaries of the ray volume (see crosses in Fig. B1). To account for a corresponding change in the area of the ray volume, we interpolate all needed background fields horizontally to the face center positions and the ray-volume center through linear Taylor expansions as follows:
χ(rj)=χ(rcc)+[rχ]cc(rjrcc),
where χ symbolizes the interpolated field and the subscripts cc and j refer to evaluations at the cell center and the ray-volume face or ray-volume center, respectively. In the vertical, we apply linear interpolations between the grid points. While the tendencies for the wavenumbers are computed with the ray-volume centered values, the interpolated fields at the ray-volume faces are used to compute group velocities, enabling the prediction of the changing extent of the ray volume in physical space [Eqs. (1) and (2)]. In particular, the difference in the group velocities at the center of opposing ray-volume faces predicts the rate of compression along the corresponding direction. The average of the opposing group velocities is then used to predict the change in location along the characteristic [Eq. (3)]. Consequently, the change in the wavenumber extents Δ(kλ, kϕ, kr) is determined following the procedure described in section 3a.

b. Horizontal propagation and code parallelization

Since MS-GWaM propagates Lagrangian ray volumes through space, it cannot easily be parallelized using the traditional MPI infrastructure of the ICON model developed for the synchronization of Eulerian fields. In particular, the ray volumes may propagate over large distances and transfer between parallel MPI domains. To accompany this while avoiding load balancing problems, we associate each ray volume with a parent cell based on its current location and store it in a corresponding grid-based array. As an example, in the shown sketch (Fig. B1), the hatched ray volume is associated with the white-bounded underlying cell. It may, however, have a group velocity such that it propagates into a new parent cell, represented by the arrow and the dark gray shaded triangle, respectively. Therefore, after each propagation step, all cells such as the dark gray shaded triangle search all rays in all neighboring cells (light gray shading) and transfer those that are now located in their own cell area. With this construction, we can localize arrays of ray volumes based on the Eulerian grid column. The integration along the wave characteristics, launch processes, saturation calculations, and wave projection onto the Eulerian grid can then be parallelized using the domain decomposition of the dynamical core with the given model infrastructure. Note that the modulation equations are integrated using a three-step Runge–Kutta low-storage scheme (Williamson 1980).

Fig. B1.
Fig. B1.

Sketch of the horizontal extent of a Lagrangian ray volume propagating in the unstructured triangular grid of the ICON model. The propagating ray volume (hatched box) may be associated with a current parent column (framed white grid cell). Given a horizontal group velocity (black arrow), we may expect it to propagate into a neighboring column (dark gray shaded grid cell). The crosses at the edges of the ray volume mark the locations for which background field data are needed to calculate the propagation of the ray volume.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

c. Ray-volume splitting

In general, the numerical treatment of internal gravity waves as ray volumes is set up to be exact in the limit of infinitesimally small ray volumes and Eulerian grid cells. The actual sizes of both are, however, finite. While the total phase space volume of the rays is constant by construction, they may stretch or compress in either physical space or spectral space. To avoid large numerical artifacts, we therefore split all ray volumes j vertically, which become larger than the local thresholds:
Δrj>max(αshc,1km),
with the local gridcell height hc and a scale factor αs (which is set to 2.5 in this study). The volume is cut into two equally sized ray volumes in the vertical, with the wavevectors identical to the original parent wave. The carrier ray location of each split volume is then set to its center. Similarly, ray volumes are split horizontally where either of the three thresholds is reached:
cosϕjΔλj>αsac/R,Δϕj>αsac/R,
cosϕjΔλjΔϕj>αsac/R2.
Here, ac represents the local cell area and R represents Earth’s mean radius. Splitting is done along the direction of the larger extent.

d. Horizontal flux smoothing

To avoid numerical artifacts based on the employed simplification in the ray volume representation, we smooth all diagnostic outputs and wave fluxes. Each triangular cell combines three vertices, one in each corner (cf. Fig. B1). Each of these vertices is connected to five or six neighboring cells, depending on the location of the triangulated icosahedron. The smoothing algorithm first averages the corresponding fields over these pentagons or hexagons and then combines the three resulting values into one average. Consequently, the weights are distributed such that the center cell is weighted with a factor 3, the direct adjacent cells with one common edge are weighted with a factor 2, and all neighboring cells sharing a vertex with the center cell are weighted with a factor 1. While this procedure generally reduces the amplitude of strong flux peaks, it also reduces the occurrence of grid-scale noise and ensures model stability.

e. Computational performance

Finally, conceptional models need to prove themselves not only through the accurate representation of not-resolved physics but also through competitive performance. In particular, ray-tracing models may suffer from load balancing problems due to the uneven distribution of ray integrations between processors but also the existence of too many rays in general. As described in section b of appendix B, MS-GWaM utilizes a cell-based representation of ray volumes, which is highly parallelizable such that load balancing problems may not inhibit the model performance. Additionally, the total number of ray volumes is restricted by the maximum numbers nbg ≤ 4800 and nconv ≤ 8800 per model-grid column for background and convective gravity waves, respectively. For the lowest impact, the ray volumes with the lowest wave energy are chosen to be removed. Detailed analyses show that the impact of this procedure is indeed small and may be neglected (not shown).

To measure the performance of the model, we compare the runtimes of all 2-month runs without any nonorographic gravity wave parameterization with both the MS-GWaM-1D and MS-GWaM-3D setups (Fig. B2). Additionally, we show the fraction of runtime that is associated with MS-GWaM and normalize all values with the number of simulated days. These measures show that MS-GWaM-1D and MS-GWaM-3D slow down ICON by factors of approximately 12 and 30, respectively. That is, the 3D setup has a performance penalty equal to a factor of approximately 2.5. Still, the gravity wave parameterization is several orders of magnitude faster than wave-resolving simulations—for which Bölöni et al. (2021) estimate a cost factor of about six orders of magnitude—and thus a competitive tool for the estimate and analysis of global IGW effects.

Fig. B2.
Fig. B2.

Performance comparison for UA-ICON at the R2B4 resolution (∼160 km) parallelized over 1920 MPI cores for runs without IGW parameterization (gray), MS-GWaM-1D (black), and MS-GWaM-3D (red). In addition to the total run time of ICON per simulated day, we show the fraction of run time associated with MS-GWaM-1D and MS-GWaM-3D. The median values for the three setups are 25 , 317, and 743 s day−1, and thus, MS-GWaM-1D and MS-GWaM-3D slow down the model by about a factor of 12 and 30, respectively.

Citation: Journal of the Atmospheric Sciences 81, 7; 10.1175/JAS-D-23-0153.1

REFERENCES

  • Achatz, U., 2022: Atmospheric Dynamics. 1st ed. Springer, 554 pp.

  • Achatz, U., R. Klein, and F. Senf, 2010: Gravity waves, scale asymptotics and the pseudo-incompressible equations. J. Fluid Mech., 663, 120147, https://doi.org/10.1017/S0022112010003411.

    • Search Google Scholar
    • Export Citation
  • Achatz, U., B. Ribstein, F. Senf, and R. Klein, 2017: The interaction between synoptic-scale balanced flow and a finite-amplitude mesoscale wave field throughout all atmospheric layers: Weak and moderately strong stratification. Quart. J. Roy. Meteor. Soc., 143, 342361, https://doi.org/10.1002/qj.2926.

    • Search Google Scholar
    • Export Citation
  • Achatz, U., Y.-H. Kim, and G. S. Voelker, 2023: Multi-scale dynamics of the interaction between waves and mean flows: From nonlinear WKB theory to gravity-wave parameterizations in weather and climate models. J. Math. Phys., 64, 111101, https://doi.org/10.1063/5.0165180.

    • Search Google Scholar
    • Export Citation
  • Alexander, M. J., and T. J. Dunkerton, 1999: A spectral parameterization of mean-flow forcing due to breaking gravity waves. J. Atmos. Sci., 56, 41674182, https://doi.org/10.1175/1520-0469(1999)056<4167:ASPOMF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Alexander, M. J., and Coauthors, 2010: Recent developments in gravity-wave effects in climate models and the global distribution of gravity-wave momentum flux from observations and models. Quart. J. Roy. Meteor. Soc., 136, 11031124, https://doi.org/10.1002/qj.637.

    • Search Google Scholar
    • Export Citation
  • Amemiya, A., and K. Sato, 2016: A new gravity wave parameterization including three-dimensional propagation. J. Meteor. Soc. Japan, 94, 237256, https://doi.org/10.2151/jmsj.2016-013.

    • Search Google Scholar
    • Export Citation
  • Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen-Palm relation and the mean zonal acceleration. J. Atmos. Sci., 33, 20312048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Bacmeister, J. T., P. A. Newman, B. L. Gary, and K. R. Chan, 1994: An algorithm for forecasting mountain wave–related turbulence in the stratosphere. Wea. Forecasting, 9, 241253, https://doi.org/10.1175/1520-0434(1994)009<0241:AAFFMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Baldwin, M. P., and Coauthors, 2021: Sudden stratospheric warmings. Rev. Geophys., 59, e2020RG000708, https://doi.org/10.1029/2020RG000708.

    • Search Google Scholar
    • Export Citation
  • Bölöni, G., B. Ribstein, J. Muraschko, C. Sgoff, J. Wei, and U. Achatz, 2016: The interaction between atmospheric gravity waves and large-scale flows: An efficient description beyond the nonacceleration paradigm. J. Atmos. Sci., 73, 48334852, https://doi.org/10.1175/JAS-D-16-0069.1.

    • Search Google Scholar
    • Export Citation
  • Bölöni, G., Y.-H. Kim, S. Borchert, and U. Achatz, 2021: Toward transient subgrid-scale gravity wave representation in atmospheric models. Part I: Propagation model including nondissipative wave–mean-flow interactions. J. Atmos. Sci., 78, 13171338, https://doi.org/10.1175/JAS-D-20-0065.1.

    • Search Google Scholar
    • Export Citation
  • Borchert, S., G. Zhou, M. Baldauf, H. Schmidt, G. Zängl, and D. Reinert, 2019: The upper-atmosphere extension of the ICON general circulation model (version: ua-icon-1.0). Geosci. Model Dev., 12, 35413569, https://doi.org/10.5194/gmd-12-3541-2019.

    • Search Google Scholar
    • Export Citation
  • Bretherton, F. P., 1966: The propagation of groups of internal gravity waves in a shear flow. Quart. J. Roy. Meteor. Soc., 92, 466480, https://doi.org/10.1002/qj.49709239403.

    • Search Google Scholar
    • Export Citation
  • Charron, M., and E. Manzini, 2002: Gravity waves from fronts: Parameterization and middle atmosphere response in a general circulation model. J. Atmos. Sci., 59, 923941, https://doi.org/10.1175/1520-0469(2002)059<0923:GWFFPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • De la Cámara, A., and F. Lott, 2015: A parameterization of gravity waves emitted by fronts and jets. Geophys. Res. Lett., 42, 20712078, https://doi.org/10.1002/2015GL063298.

    • Search Google Scholar
    • Export Citation
  • De la Cámara, A., F. Lott, V. Jewtoukoff, R. Plougonven, and A. Hertzog, 2016: On the gravity wave forcing during the southern stratospheric final warming in LMDZ. J. Atmos. Sci., 73, 32133226, https://doi.org/10.1175/JAS-D-15-0377.1.

    • Search Google Scholar
    • Export Citation
  • Drob, D. P., and Coauthors, 2015: An update to the horizontal wind model (HWM): The quiet time thermosphere. Earth Space Sci., 2, 301319, https://doi.org/10.1002/2014EA000089.

    • Search Google Scholar
    • Export Citation
  • Dunkerton, T. J., and N. Butchart, 1984: Propagation and selective transmission of internal gravity waves in a sudden warming. J. Atmos. Sci., 41, 14431460, https://doi.org/10.1175/1520-0469(1984)041<1443:PASTOI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ehard, B., and Coauthors, 2017: Horizontal propagation of large-amplitude mountain waves into the polar night jet. J. Geophys. Res. Atmos., 122, 14231436, https://doi.org/10.1002/2016JD025621.

    • Search Google Scholar
    • Export Citation
  • Eichinger, R., and Coauthors, 2023: Emulating lateral gravity wave propagation in a global chemistry–climate model (EMAC v2.55.2) through horizontal flux redistribution. Geosci. Model Dev., 16, 55615583, https://doi.org/10.5194/gmd-16-5561-2023.

    • Search Google Scholar
    • Export Citation
  • Eliassen, A., and E. Palm, 1961: On the transfer of energy in stationary mountain waves. Geofys. Publ., 22 (3), 123.

  • Emmert, J. T., and Coauthors, 2021: NRLMSIS 2.0: A whole-atmosphere empirical model of temperature and neutral species densities. Earth Space Sci., 8, e2020EA001321, https://doi.org/10.1029/2020EA001321.

    • Search Google Scholar
    • Export Citation
  • Emmert, J. T., and Coauthors, 2022: NRLMSIS 2.1: An empirical model of nitric oxide incorporated into MSIS. J. Geophys. Res. Space Phys., 127, e2022JA030896, https://doi.org/10.1029/2022JA030896.

    • Search Google Scholar
    • Export Citation
  • Ern, M., Q. T. Trinh, P. Preusse, J. C. Gille, M. G. Mlynczak, J. M. Russell III, and M. Riese, 2018: GRACILE: A comprehensive climatology of atmospheric gravity wave parameters based on satellite limb soundings. Earth Syst. Sci. Data, 10, 857892, https://doi.org/10.5194/essd-10-857-2018.

    • Search Google Scholar
    • Export Citation
  • Fritts, D. C., and M. J. Alexander, 2003: Gravity wave dynamics and effects in the middle atmosphere. Rev. Geophys., 41, 1003, https://doi.org/10.1029/2001RG000106.

    • Search Google Scholar
    • Export Citation
  • Grimshaw, R., 1974: Internal gravity waves in a slowly varying, dissipative medium. Geophys. Fluid Dyn., 6, 131148, https://doi.org/10.1080/03091927409365792.

    • Search Google Scholar
    • Export Citation
  • Grimshaw, R., 1975: Nonlinear internal gravity waves in a rotating fluid. J. Fluid Mech., 71, 497512, https://doi.org/10.1017/S0022112075002704.

    • Search Google Scholar
    • Export Citation
  • Gupta, A., T. Birner, A. Dörnbrack, and I. Polichtchouk, 2021: Importance of gravity wave forcing for springtime southern polar vortex breakdown as revealed by ERA5. Geophys. Res. Lett., 48, e2021GL092762, https://doi.org/10.1029/2021GL092762.

    • Search Google Scholar
    • Export Citation
  • Hasha, A., O. Bühler, and J. Scinocca, 2008: Gravity wave refraction by three-dimensionally varying winds and the global transport of angular momentum. J. Atmos. Sci., 65, 28922906, https://doi.org/10.1175/2007JAS2561.1.

    • Search Google Scholar
    • Export Citation
  • Hindley, N. P., C. J. Wright, L. Hoffmann, T. Moffat-Griffin, and N. J. Mitchell, 2020: An 18-year climatology of directional stratospheric gravity wave momentum flux from 3-D satellite observations. Geophys. Res. Lett., 47, e2020GL089557, https://doi.org/10.1029/2020GL089557.

    • Search Google Scholar
    • Export Citation
  • Hines, C. O., 1997a: Doppler-spread parameterization of gravity-wave momentum deposition in the middle atmosphere. Part 1: Basic formulation. J. Atmos. Sol.-Terr. Phys., 59, 371386, https://doi.org/10.1016/S1364-6826(96)00079-X.

    • Search Google Scholar
    • Export Citation
  • Hines, C. O., 1997b: Doppler-spread parameterization of gravity-wave momentum deposition in the middle atmosphere. Part 2: Broad and quasi monochromatic spectra, and implementation. J. Atmos. Sol.-Terr. Phys., 59, 387400, https://doi.org/10.1016/S1364-6826(96)00080-6.

    • Search Google Scholar
    • Export Citation
  • Holt, L. A., M. J. Alexander, L. Coy, A. Molod, W. Putman, and S. Pawson, 2016: Tropical waves and the quasi-biennial oscillation in a 7-km global climate simulation. J. Atmos. Sci., 73, 37713783, https://doi.org/10.1175/JAS-D-15-0350.1.

    • Search Google Scholar
    • Export Citation
  • Holt, L. A., C. M. Brabec, and M. J. Alexander, 2023: Exploiting high-density zonal-sampling of HIRDLS profiles near 60°S to investigate missing drag in chemistry-climate models. J. Geophys. Res. Atmos., 128, e2022JD037398, https://doi.org/10.1029/2022JD037398.

    • Search Google Scholar
    • Export Citation
  • Jiang, J. H., S. D. Eckermann, D. L. Wu, and D. Y. Wang, 2006: Inter-annual variation of gravity waves in the Arctic and Antarctic winter middle atmosphere. Adv. Space Res., 38, 24182423, https://doi.org/10.1016/j.asr.2005.09.036.

    • Search Google Scholar
    • Export Citation
  • Kafiabad, H. A., M. A. Savva, and J. Vanneste, 2019: Diffusion of inertia-gravity waves by geostrophic turbulence. J. Fluid Mech., 869, R7, https://doi.org/10.1017/jfm.2019.300.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-H., G. Bölöni, S. Borchert, H.-Y. Chun, and U. Achatz, 2021: Toward transient subgrid-scale gravity wave representation in atmospheric models. Part II: Wave intermittency simulated with convective sources. J. Atmos. Sci., 78, 13391357, https://doi.org/10.1175/JAS-D-20-0066.1.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-H., G. S. Voelker, G. Bölöni, G. Zängl, and U. Achatz, 2023: Crucial role of obliquely propagating gravity waves in the quasi-biennial oscillation dynamics. arXiv:2309.15301v2, https://doi.org/10.48550/arXiv.2309.15301.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-J., S. D. Eckermann, and H.-Y. Chun, 2003: An overview of the past, present and future of gravity-wave drag parametrization for numerical climate and weather prediction models. Atmos.–Ocean, 41, 6598, https://doi.org/10.3137/ao.410105.

    • Search Google Scholar
    • Export Citation
  • Kruse, C. G., and Coauthors, 2022: Observed and modeled mountain waves from the surface to the mesosphere near the Drake Passage. J. Atmos. Sci., 79, 909932, https://doi.org/10.1175/JAS-D-21-0252.1.

    • Search Google Scholar
    • Export Citation
  • Lindzen, R. S., 1981: Turbulence and stress owing to gravity wave and tidal breakdown. J. Geophys. Res., 86, 97079714, https://doi.org/10.1029/JC086iC10p09707.

    • Search Google Scholar
    • Export Citation
  • Lott, F., and M. J. Miller, 1997: A new subgrid-scale orographic drag parametrization: Its formulation and testing. Quart. J. Roy. Meteor. Soc., 123, 101127, https://doi.org/10.1002/qj.49712353704.

    • Search Google Scholar
    • Export Citation
  • Lott, F., and L. Guez, 2013: A stochastic parameterization of the gravity waves due to convection and its impact on the equatorial stratosphere. J. Geophys. Res. Atmos., 118, 88978909, https://doi.org/10.1002/jgrd.50705.

    • Search Google Scholar
    • Export Citation
  • Marks, C. J., and S. D. Eckermann, 1995: A three-dimensional nonhydrostatic ray tracing model for gravity waves: Formulation and preliminary results for the middle atmosphere. J. Atmos. Sci., 52, 19591984, https://doi.org/10.1175/1520-0469(1995)052%3C1959:ATDNRT%3E2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Medvedev, A. S., and G. P. Klaassen, 1995: Vertical evolution of gravity wave spectra and the parameterization of associated wave drag. J. Geophys. Res., 100, 25 84125 853, https://doi.org/10.1029/95JD02533.

    • Search Google Scholar
    • Export Citation
  • Muraschko, J., M. D. Fruman, U. Achatz, S. Hickel, and Y. Toledo, 2015: On the application of Wentzel–Kramer–Brillouin theory for the simulation of the weakly nonlinear dynamics of gravity waves. Quart. J. Roy. Meteor. Soc., 141, 676697, https://doi.org/10.1002/qj.2381.

    • Search Google Scholar
    • Export Citation
  • Nappo, C. J., 2013: An Introduction to Atmospheric Gravity Waves. Elsevier Science, 359 pp.

  • Orr, A., P. Bechtold, J. Scinocca, M. Ern, and M. Janiskova, 2010: Improved middle atmosphere climate and forecasts in the ECMWF model through a nonorographic gravity wave drag parameterization. J. Climate, 23, 59055926, https://doi.org/10.1175/2010JCLI3490.1.

    • Search Google Scholar
    • Export Citation
  • Palmer, T. N., G. J. Shutts, and R. Swinbank, 1986: Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization. Quart. J. Roy. Meteor. Soc., 112, 10011039, https://doi.org/10.1002/qj.49711247406.

    • Search Google Scholar
    • Export Citation
  • Polichtchouk, I., N. Wedi, and Y.-H. Kim, 2022: Resolved gravity waves in the tropical stratosphere: Impact of horizontal resolution and deep convection parametrization. Quart. J. Roy. Meteor. Soc., 148, 233251, https://doi.org/10.1002/qj.4202.

    • Search Google Scholar
    • Export Citation
  • Polichtchouk, I., A. V. Niekerk, and N. Wedi, 2023: Resolved gravity waves in the extratropical stratosphere: Effect of horizontal resolution increase from O(10) to O(1) km. J. Atmos. Sci., 80, 473486, https://doi.org/10.1175/JAS-D-22-0138.1.

    • Search Google Scholar
    • Export Citation
  • Quinn, B., C. Eden, and D. Olbers, 2020: Application of the IDEMIX concept for internal gravity waves in the atmosphere. J. Atmos. Sci., 77, 36013618, https://doi.org/10.1175/JAS-D-20-0107.1.

    • Search Google Scholar
    • Export Citation
  • Ribstein, B., and U. Achatz, 2016: The interaction between gravity waves and solar tides in a linear tidal model with a 4-D ray-tracing gravity-wave parameterization. J. Geophys. Res. Space Phys., 121, 89368950, https://doi.org/10.1002/2016JA022478.

    • Search Google Scholar
    • Export Citation
  • Ribstein, B., U. Achatz, and F. Senf, 2015: The interaction between gravity waves and solar tides: Results from 4-D ray tracing coupled to a linear tidal model. J. Geophys. Res. Space Phys., 120, 67956817, https://doi.org/10.1002/2015JA021349.

    • Search Google Scholar
    • Export Citation
  • Richter, J. H., F. Sassi, and R. R. Garcia, 2010: Toward a physically based gravity wave source parameterization in a general circulation model. J. Atmos. Sci., 67, 136156, https://doi.org/10.1175/2009JAS3112.1.

    • Search Google Scholar
    • Export Citation
  • Šácha, P., F. Lilienthal, C. Jacobi, and P. Pisoft, 2016: Influence of the spatial distribution of gravity wave activity on the middle atmospheric dynamics. Atmos. Chem. Phys., 16, 15 75515 775, https://doi.org/10.5194/acp-16-15755-2016.

    • Search Google Scholar
    • Export Citation
  • Sato, K., S. Watanabe, Y. Kawatani, Y. Tomikawa, K. Miyazaki, and M. Takahashi, 2009: On the origins of mesospheric gravity waves. Geophys. Res. Lett., 36, L1980, https://doi.org/10.1029/2009GL039908.

    • Search Google Scholar
    • Export Citation
  • Sato, K., S. Tateno, S. Watanabe, and Y. Kawatani, 2012: Gravity wave characteristics in the Southern Hemisphere revealed by a high-resolution middle-atmosphere general circulation model. J. Atmos. Sci., 69, 13781396, https://doi.org/10.1175/JAS-D-11-0101.1.

    • Search Google Scholar
    • Export Citation
  • Scaife, A. A., J. R. Knight, G. K. Vallis, and C. K. Folland, 2005: A stratospheric influence on the winter NAO and North Atlantic surface climate. Geophys. Res. Lett., 32, L18715, https://doi.org/10.1029/2005GL023226.

    • Search Google Scholar
    • Export Citation
  • Scaife, A. A., and Coauthors, 2012: Climate change projections and stratosphere–troposphere interaction. Climate Dyn., 38, 20892097, https://doi.org/10.1007/s00382-011-1080-7.

    • Search Google Scholar
    • Export Citation
  • Scinocca, J. F., 2003: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models. J. Atmos. Sci., 60, 667682, https://doi.org/10.1175/1520-0469(2003)060<0667:AASNGW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sigmond, M., and J. F. Scinocca, 2010: The influence of the basic state on the Northern Hemisphere circulation response to climate change. J. Climate, 23, 14341446, https://doi.org/10.1175/2009JCLI3167.1.

    • Search Google Scholar
    • Export Citation
  • Sutherland, B. R., 2010: Internal Gravity Waves. Cambridge University Press, 377 pp.

  • Sutherland, B. R., U. Achatz, C.-C. P. Caulfield, and J. M. Klymak, 2019: Recent progress in modeling imbalance in the atmosphere and ocean. Phys. Rev. Fluids, 4, 010501, https://doi.org/10.1103/PhysRevFluids.4.010501.

    • Search Google Scholar
    • Export Citation
  • van Niekerk, A., and S. B. Vosper, 2021: Towards a more “scale-aware” orographic gravity wave drag parametrization: Description and initial testing. Quart. J. Roy. Meteor. Soc., 147, 32433262, https://doi.org/10.1002/qj.4126.

    • Search Google Scholar
    • Export Citation
  • van Niekerk, A., S. B. Vosper, and M. A. C. Teixeira, 2023: Accounting for the three-dimensional nature of mountain waves: Parametrising partial critical-level filtering. Quart. J. Roy. Meteor. Soc., 149, 515536, https://doi.org/10.1002/qj.4421.

    • Search Google Scholar
    • Export Citation
  • Voelker, G. S., T. R. Akylas, and U. Achatz, 2021: An application of WKBJ theory for triad interactions of internal gravity waves in varying background flows. Quart. J. Roy. Meteor. Soc., 147, 11121134, https://doi.org/10.1002/qj.3962.

    • Search Google Scholar
    • Export Citation
  • Warner, C. D., and M. E. McIntyre, 1996: On the propagation and dissipation of gravity wave spectra through a realistic middle atmosphere. J. Atmos. Sci., 53, 32133235, https://doi.org/10.1175/1520-0469(1996)053<3213:OTPADO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Warner, C. D., and M. E. McIntyre, 2001: An ultrasimple spectral parameterization for nonorographic gravity waves. J. Atmos. Sci., 58, 18371857, https://doi.org/10.1175/1520-0469(2001)058<1837:AUSPFN>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Wei, J., G. Bölöni, and U. Achatz, 2019: Efficient modeling of the interaction of mesoscale gravity waves with unbalanced large-scale flows: Pseudomomentum-Flux convergence versus direct approach. J. Atmos. Sci., 76, 27152738, https://doi.org/10.1175/JAS-D-18-0337.1.

    • Search Google Scholar
    • Export Citation
  • Wilhelm, J., T. R. Akylas, G. Bölöni, J. Wei, B. Ribstein, R. Klein, and U. Achatz, 2018: Interactions between mesoscale and submesoscale gravitywaves and their efficient representation in mesoscale-resolving models. J. Atmos. Sci., 75, 22572280, https://doi.org/10.1175/JAS-D-17-0289.1.

    • Search Google Scholar
    • Export Citation
  • Williams, P. D., and Coauthors, 2017: A census of atmospheric variability from seconds to decades. Geophys. Res. Lett., 44, 11 20111 211, https://doi.org/10.1002/2017GL075483.

    • Search Google Scholar
    • Export Citation
  • Williamson, J. H., 1980: Low-storage Runge-Kutta schemes. J. Comput. Phys., 35, 4856, https://doi.org/10.1016/0021-9991(80)90033-9.

  • Xie, J., M. Zhang, Z. Xie, H. Liu, Z. Chai, J. He, and H. Zhang, 2020: An orographic-drag parametrization scheme including orographic anisotropy for all flow directions. J. Adv. Model. Earth Syst., 12, e2019MS001921, https://doi.org/10.1029/2019MS001921.

    • Search Google Scholar
    • Export Citation
  • Xie, J., M. Zhang, Q. Zeng, Z. Xie, H. Liu, Z. Chai, J. He, and H. Zhang, 2021: Implementation of an orographic drag scheme considering orographic anisotropy in all flow directions in the Earth System Model CAS-ESM 2.0. J. Adv. Model. Earth Syst., 13, e2021MS002585, https://doi.org/10.1029/2021MS002585.

    • Search Google Scholar
    • Export Citation
  • Zängl, G., D. Reinert, P. Rípodas, and M. Baldauf, 2015: The ICON (ICOsahedral non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core. Quart. J. Roy. Meteor. Soc., 141, 563579, https://doi.org/10.1002/qj.2378.

    • Search Google Scholar
    • Export Citation
Save
  • Achatz, U., 2022: Atmospheric Dynamics. 1st ed. Springer, 554 pp.

  • Achatz, U., R. Klein, and F. Senf, 2010: Gravity waves, scale asymptotics and the pseudo-incompressible equations. J. Fluid Mech., 663, 120147, https://doi.org/10.1017/S0022112010003411.

    • Search Google Scholar
    • Export Citation
  • Achatz, U., B. Ribstein, F. Senf, and R. Klein, 2017: The interaction between synoptic-scale balanced flow and a finite-amplitude mesoscale wave field throughout all atmospheric layers: Weak and moderately strong stratification. Quart. J. Roy. Meteor. Soc., 143, 342361, https://doi.org/10.1002/qj.2926.

    • Search Google Scholar
    • Export Citation
  • Achatz, U., Y.-H. Kim, and G. S. Voelker, 2023: Multi-scale dynamics of the interaction between waves and mean flows: From nonlinear WKB theory to gravity-wave parameterizations in weather and climate models. J. Math. Phys., 64, 111101, https://doi.org/10.1063/5.0165180.

    • Search Google Scholar
    • Export Citation
  • Alexander, M. J., and T. J. Dunkerton, 1999: A spectral parameterization of mean-flow forcing due to breaking gravity waves. J. Atmos. Sci., 56, 41674182, https://doi.org/10.1175/1520-0469(1999)056<4167:ASPOMF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Alexander, M. J., and Coauthors, 2010: Recent developments in gravity-wave effects in climate models and the global distribution of gravity-wave momentum flux from observations and models. Quart. J. Roy. Meteor. Soc., 136, 11031124, https://doi.org/10.1002/qj.637.

    • Search Google Scholar
    • Export Citation
  • Amemiya, A., and K. Sato, 2016: A new gravity wave parameterization including three-dimensional propagation. J. Meteor. Soc. Japan, 94, 237256, https://doi.org/10.2151/jmsj.2016-013.

    • Search Google Scholar
    • Export Citation
  • Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen-Palm relation and the mean zonal acceleration. J. Atmos. Sci., 33, 20312048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Bacmeister, J. T., P. A. Newman, B. L. Gary, and K. R. Chan, 1994: An algorithm for forecasting mountain wave–related turbulence in the stratosphere. Wea. Forecasting, 9, 241253, https://doi.org/10.1175/1520-0434(1994)009<0241:AAFFMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Baldwin, M. P., and Coauthors, 2021: Sudden stratospheric warmings. Rev. Geophys., 59, e2020RG000708, https://doi.org/10.1029/2020RG000708.

    • Search Google Scholar
    • Export Citation
  • Bölöni, G., B. Ribstein, J. Muraschko, C. Sgoff, J. Wei, and U. Achatz, 2016: The interaction between atmospheric gravity waves and large-scale flows: An efficient description beyond the nonacceleration paradigm. J. Atmos. Sci., 73, 48334852, https://doi.org/10.1175/JAS-D-16-0069.1.

    • Search Google Scholar
    • Export Citation
  • Bölöni, G., Y.-H. Kim, S. Borchert, and U. Achatz, 2021: Toward transient subgrid-scale gravity wave representation in atmospheric models. Part I: Propagation model including nondissipative wave–mean-flow interactions. J. Atmos. Sci., 78, 13171338, https://doi.org/10.1175/JAS-D-20-0065.1.

    • Search Google Scholar
    • Export Citation
  • Borchert, S., G. Zhou, M. Baldauf, H. Schmidt, G. Zängl, and D. Reinert, 2019: The upper-atmosphere extension of the ICON general circulation model (version: ua-icon-1.0). Geosci. Model Dev., 12, 35413569, https://doi.org/10.5194/gmd-12-3541-2019.

    • Search Google Scholar
    • Export Citation
  • Bretherton, F. P., 1966: The propagation of groups of internal gravity waves in a shear flow. Quart. J. Roy. Meteor. Soc., 92, 466480, https://doi.org/10.1002/qj.49709239403.

    • Search Google Scholar
    • Export Citation
  • Charron, M., and E. Manzini, 2002: Gravity waves from fronts: Parameterization and middle atmosphere response in a general circulation model. J. Atmos. Sci., 59, 923941, https://doi.org/10.1175/1520-0469(2002)059<0923:GWFFPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • De la Cámara, A., and F. Lott, 2015: A parameterization of gravity waves emitted by fronts and jets. Geophys. Res. Lett., 42, 20712078, https://doi.org/10.1002/2015GL063298.

    • Search Google Scholar
    • Export Citation
  • De la Cámara, A., F. Lott, V. Jewtoukoff, R. Plougonven, and A. Hertzog, 2016: On the gravity wave forcing during the southern stratospheric final warming in LMDZ. J. Atmos. Sci., 73, 32133226, https://doi.org/10.1175/JAS-D-15-0377.1.

    • Search Google Scholar
    • Export Citation
  • Drob, D. P., and Coauthors, 2015: An update to the horizontal wind model (HWM): The quiet time thermosphere. Earth Space Sci., 2, 301319, https://doi.org/10.1002/2014EA000089.

    • Search Google Scholar
    • Export Citation
  • Dunkerton, T. J., and N. Butchart, 1984: Propagation and selective transmission of internal gravity waves in a sudden warming. J. Atmos. Sci., 41, 14431460, https://doi.org/10.1175/1520-0469(1984)041<1443:PASTOI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ehard, B., and Coauthors, 2017: Horizontal propagation of large-amplitude mountain waves into the polar night jet. J. Geophys. Res. Atmos., 122, 14231436, https://doi.org/10.1002/2016JD025621.

    • Search Google Scholar
    • Export Citation
  • Eichinger, R., and Coauthors, 2023: Emulating lateral gravity wave propagation in a global chemistry–climate model (EMAC v2.55.2) through horizontal flux redistribution. Geosci. Model Dev., 16, 55615583, https://doi.org/10.5194/gmd-16-5561-2023.

    • Search Google Scholar
    • Export Citation
  • Eliassen, A., and E. Palm, 1961: On the transfer of energy in stationary mountain waves. Geofys. Publ., 22 (3), 123.

  • Emmert, J. T., and Coauthors, 2021: NRLMSIS 2.0: A whole-atmosphere empirical model of temperature and neutral species densities. Earth Space Sci., 8, e2020EA001321, https://doi.org/10.1029/2020EA001321.

    • Search Google Scholar
    • Export Citation
  • Emmert, J. T., and Coauthors, 2022: NRLMSIS 2.1: An empirical model of nitric oxide incorporated into MSIS. J. Geophys. Res. Space Phys., 127, e2022JA030896, https://doi.org/10.1029/2022JA030896.

    • Search Google Scholar
    • Export Citation
  • Ern, M., Q. T. Trinh, P. Preusse, J. C. Gille, M. G. Mlynczak, J. M. Russell III, and M. Riese, 2018: GRACILE: A comprehensive climatology of atmospheric gravity wave parameters based on satellite limb soundings. Earth Syst. Sci. Data, 10, 857892, https://doi.org/10.5194/essd-10-857-2018.

    • Search Google Scholar
    • Export Citation
  • Fritts, D. C., and M. J. Alexander, 2003: Gravity wave dynamics and effects in the middle atmosphere. Rev. Geophys., 41, 1003, https://doi.org/10.1029/2001RG000106.

    • Search Google Scholar
    • Export Citation
  • Grimshaw, R., 1974: Internal gravity waves in a slowly varying, dissipative medium. Geophys. Fluid Dyn., 6, 131148, https://doi.org/10.1080/03091927409365792.

    • Search Google Scholar
    • Export Citation
  • Grimshaw, R., 1975: Nonlinear internal gravity waves in a rotating fluid. J. Fluid Mech., 71, 497512, https://doi.org/10.1017/S0022112075002704.

    • Search Google Scholar
    • Export Citation
  • Gupta, A., T. Birner, A. Dörnbrack, and I. Polichtchouk, 2021: Importance of gravity wave forcing for springtime southern polar vortex breakdown as revealed by ERA5. Geophys. Res. Lett., 48, e2021GL092762, https://doi.org/10.1029/2021GL092762.

    • Search Google Scholar
    • Export Citation
  • Hasha, A., O. Bühler, and J. Scinocca, 2008: Gravity wave refraction by three-dimensionally varying winds and the global transport of angular momentum. J. Atmos. Sci., 65, 28922906, https://doi.org/10.1175/2007JAS2561.1.

    • Search Google Scholar
    • Export Citation
  • Hindley, N. P., C. J. Wright, L. Hoffmann, T. Moffat-Griffin, and N. J. Mitchell, 2020: An 18-year climatology of directional stratospheric gravity wave momentum flux from 3-D satellite observations. Geophys. Res. Lett., 47, e2020GL089557, https://doi.org/10.1029/2020GL089557.

    • Search Google Scholar
    • Export Citation
  • Hines, C. O., 1997a: Doppler-spread parameterization of gravity-wave momentum deposition in the middle atmosphere. Part 1: Basic formulation. J. Atmos. Sol.-Terr. Phys., 59, 371386, https://doi.org/10.1016/S1364-6826(96)00079-X.

    • Search Google Scholar
    • Export Citation
  • Hines, C. O., 1997b: Doppler-spread parameterization of gravity-wave momentum deposition in the middle atmosphere. Part 2: Broad and quasi monochromatic spectra, and implementation. J. Atmos. Sol.-Terr. Phys., 59, 387400, https://doi.org/10.1016/S1364-6826(96)00080-6.

    • Search Google Scholar
    • Export Citation
  • Holt, L. A., M. J. Alexander, L. Coy, A. Molod, W. Putman, and S. Pawson, 2016: Tropical waves and the quasi-biennial oscillation in a 7-km global climate simulation. J. Atmos. Sci., 73, 37713783, https://doi.org/10.1175/JAS-D-15-0350.1.

    • Search Google Scholar
    • Export Citation
  • Holt, L. A., C. M. Brabec, and M. J. Alexander, 2023: Exploiting high-density zonal-sampling of HIRDLS profiles near 60°S to investigate missing drag in chemistry-climate models. J. Geophys. Res. Atmos., 128, e2022JD037398, https://doi.org/10.1029/2022JD037398.

    • Search Google Scholar
    • Export Citation
  • Jiang, J. H., S. D. Eckermann, D. L. Wu, and D. Y. Wang, 2006: Inter-annual variation of gravity waves in the Arctic and Antarctic winter middle atmosphere. Adv. Space Res., 38, 24182423, https://doi.org/10.1016/j.asr.2005.09.036.

    • Search Google Scholar
    • Export Citation
  • Kafiabad, H. A., M. A. Savva, and J. Vanneste, 2019: Diffusion of inertia-gravity waves by geostrophic turbulence. J. Fluid Mech., 869, R7, https://doi.org/10.1017/jfm.2019.300.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-H., G. Bölöni, S. Borchert, H.-Y. Chun, and U. Achatz, 2021: Toward transient subgrid-scale gravity wave representation in atmospheric models. Part II: Wave intermittency simulated with convective sources. J. Atmos. Sci., 78, 13391357, https://doi.org/10.1175/JAS-D-20-0066.1.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-H., G. S. Voelker, G. Bölöni, G. Zängl, and U. Achatz, 2023: Crucial role of obliquely propagating gravity waves in the quasi-biennial oscillation dynamics. arXiv:2309.15301v2, https://doi.org/10.48550/arXiv.2309.15301.

    • Search Google Scholar
    • Export Citation
  • Kim, Y.-J., S. D. Eckermann, and H.-Y. Chun, 2003: An overview of the past, present and future of gravity-wave drag parametrization for numerical climate and weather prediction models. Atmos.–Ocean, 41, 6598, https://doi.org/10.3137/ao.410105.

    • Search Google Scholar
    • Export Citation
  • Kruse, C. G., and Coauthors, 2022: Observed and modeled mountain waves from the surface to the mesosphere near the Drake Passage. J. Atmos. Sci., 79, 909932, https://doi.org/10.1175/JAS-D-21-0252.1.

    • Search Google Scholar
    • Export Citation
  • Lindzen, R. S., 1981: Turbulence and stress owing to gravity wave and tidal breakdown. J. Geophys. Res., 86, 97079714, https://doi.org/10.1029/JC086iC10p09707.

    • Search Google Scholar
    • Export Citation
  • Lott, F., and M. J. Miller, 1997: A new subgrid-scale orographic drag parametrization: Its formulation and testing. Quart. J. Roy. Meteor. Soc., 123, 101127, https://doi.org/10.1002/qj.49712353704.

    • Search Google Scholar
    • Export Citation
  • Lott, F., and L. Guez, 2013: A stochastic parameterization of the gravity waves due to convection and its impact on the equatorial stratosphere. J. Geophys. Res. Atmos., 118, 88978909, https://doi.org/10.1002/jgrd.50705.

    • Search Google Scholar
    • Export Citation
  • Marks, C. J., and S. D. Eckermann, 1995: A three-dimensional nonhydrostatic ray tracing model for gravity waves: Formulation and preliminary results for the middle atmosphere. J. Atmos. Sci., 52, 19591984, https://doi.org/10.1175/1520-0469(1995)052%3C1959:ATDNRT%3E2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Medvedev, A. S., and G. P. Klaassen, 1995: Vertical evolution of gravity wave spectra and the parameterization of associated wave drag. J. Geophys. Res., 100, 25 84125 853, https://doi.org/10.1029/95JD02533.

    • Search Google Scholar
    • Export Citation
  • Muraschko, J., M. D. Fruman, U. Achatz, S. Hickel, and Y. Toledo, 2015: On the application of Wentzel–Kramer–Brillouin theory for the simulation of the weakly nonlinear dynamics of gravity waves. Quart. J. Roy. Meteor. Soc., 141, 676697, https://doi.org/10.1002/qj.2381.

    • Search Google Scholar
    • Export Citation
  • Nappo, C. J., 2013: An Introduction to Atmospheric Gravity Waves. Elsevier Science, 359 pp.

  • Orr, A., P. Bechtold, J. Scinocca, M. Ern, and M. Janiskova, 2010: Improved middle atmosphere climate and forecasts in the ECMWF model through a nonorographic gravity wave drag parameterization. J. Climate, 23, 59055926, https://doi.org/10.1175/2010JCLI3490.1.

    • Search Google Scholar
    • Export Citation
  • Palmer, T. N., G. J. Shutts, and R. Swinbank, 1986: Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization. Quart. J. Roy. Meteor. Soc., 112, 10011039, https://doi.org/10.1002/qj.49711247406.

    • Search Google Scholar
    • Export Citation
  • Polichtchouk, I., N. Wedi, and Y.-H. Kim, 2022: Resolved gravity waves in the tropical stratosphere: Impact of horizontal resolution and deep convection parametrization. Quart. J. Roy. Meteor. Soc., 148, 233251, https://doi.org/10.1002/qj.4202.

    • Search Google Scholar
    • Export Citation
  • Polichtchouk, I., A. V. Niekerk, and N. Wedi, 2023: Resolved gravity waves in the extratropical stratosphere: Effect of horizontal resolution increase from O(10) to O(1) km. J. Atmos. Sci., 80, 473486, https://doi.org/10.1175/JAS-D-22-0138.1.

    • Search Google Scholar
    • Export Citation
  • Quinn, B., C. Eden, and D. Olbers, 2020: Application of the IDEMIX concept for internal gravity waves in the atmosphere. J. Atmos. Sci., 77, 36013618, https://doi.org/10.1175/JAS-D-20-0107.1.

    • Search Google Scholar
    • Export Citation
  • Ribstein, B., and U. Achatz, 2016: The interaction between gravity waves and solar tides in a linear tidal model with a 4-D ray-tracing gravity-wave parameterization. J. Geophys. Res. Space Phys., 121, 89368950, https://doi.org/10.1002/2016JA022478.

    • Search Google Scholar
    • Export Citation
  • Ribstein, B., U. Achatz, and F. Senf, 2015: The interaction between gravity waves and solar tides: Results from 4-D ray tracing coupled to a linear tidal model. J. Geophys. Res. Space Phys., 120, 67956817, https://doi.org/10.1002/2015JA021349.

    • Search Google Scholar
    • Export Citation
  • Richter, J. H., F. Sassi, and R. R. Garcia, 2010: Toward a physically based gravity wave source parameterization in a general circulation model. J. Atmos. Sci., 67, 136156, https://doi.org/10.1175/2009JAS3112.1.

    • Search Google Scholar
    • Export Citation
  • Šácha, P., F. Lilienthal, C. Jacobi, and P. Pisoft, 2016: Influence of the spatial distribution of gravity wave activity on the middle atmospheric dynamics. Atmos. Chem. Phys., 16, 15 75515 775, https://doi.org/10.5194/acp-16-15755-2016.

    • Search Google Scholar
    • Export Citation
  • Sato, K., S. Watanabe, Y. Kawatani, Y. Tomikawa, K. Miyazaki, and M. Takahashi, 2009: On the origins of mesospheric gravity waves. Geophys. Res. Lett., 36, L1980, https://doi.org/10.1029/2009GL039908.

    • Search Google Scholar
    • Export Citation
  • Sato, K., S. Tateno, S. Watanabe, and Y. Kawatani, 2012: Gravity wave characteristics in the Southern Hemisphere revealed by a high-resolution middle-atmosphere general circulation model. J. Atmos. Sci., 69, 13781396, https://doi.org/10.1175/JAS-D-11-0101.1.

    • Search Google Scholar
    • Export Citation
  • Scaife, A. A., J. R. Knight, G. K. Vallis, and C. K. Folland, 2005: A stratospheric influence on the winter NAO and North Atlantic surface climate. Geophys. Res. Lett., 32, L18715, https://doi.org/10.1029/2005GL023226.

    • Search Google Scholar
    • Export Citation
  • Scaife, A. A., and Coauthors, 2012: Climate change projections and stratosphere–troposphere interaction. Climate Dyn., 38, 20892097, https://doi.org/10.1007/s00382-011-1080-7.

    • Search Google Scholar
    • Export Citation
  • Scinocca, J. F., 2003: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models. J. Atmos. Sci., 60, 667682, https://doi.org/10.1175/1520-0469(2003)060<0667:AASNGW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sigmond, M., and J. F. Scinocca, 2010: The influence of the basic state on the Northern Hemisphere circulation response to climate change. J. Climate, 23, 14341446, https://doi.org/10.1175/2009JCLI3167.1.

    • Search Google Scholar
    • Export Citation
  • Sutherland, B. R., 2010: Internal Gravity Waves. Cambridge University Press, 377 pp.

  • Sutherland, B. R., U. Achatz, C.-C. P. Caulfield, and J. M. Klymak, 2019: Recent progress in modeling imbalance in the atmosphere and ocean. Phys. Rev. Fluids, 4, 010501, https://doi.org/10.1103/PhysRevFluids.4.010501.

    • Search Google Scholar
    • Export Citation
  • van Niekerk, A., and S. B. Vosper, 2021: Towards a more “scale-aware” orographic gravity wave drag parametrization: Description and initial testing. Quart. J. Roy. Meteor. Soc., 147, 32433262, https://doi.org/10.1002/qj.4126.

    • Search Google Scholar
    • Export Citation
  • van Niekerk, A., S. B. Vosper, and M. A. C. Teixeira, 2023: Accounting for the three-dimensional nature of mountain waves: Parametrising partial critical-level filtering. Quart. J. Roy. Meteor. Soc., 149, 515536, https://doi.org/10.1002/qj.4421.

    • Search Google Scholar
    • Export Citation
  • Voelker, G. S., T. R. Akylas, and U. Achatz, 2021: An application of WKBJ theory for triad interactions of internal gravity waves in varying background flows. Quart. J. Roy. Meteor. Soc., 147, 11121134, https://doi.org/10.1002/qj.3962.

    • Search Google Scholar
    • Export Citation
  • Warner, C. D., and M. E. McIntyre, 1996: On the propagation and dissipation of gravity wave spectra through a realistic middle atmosphere. J. Atmos. Sci., 53, 32133235, https://doi.org/10.1175/1520-0469(1996)053<3213:OTPADO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Warner, C. D., and M. E. McIntyre, 2001: An ultrasimple spectral parameterization for nonorographic gravity waves. J. Atmos. Sci., 58, 18371857, https://doi.org/10.1175/1520-0469(2001)058<1837:AUSPFN>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Wei, J., G. Bölöni, and U. Achatz, 2019: Efficient modeling of the interaction of mesoscale gravity waves with unbalanced large-scale flows: Pseudomomentum-Flux convergence versus direct approach. J. Atmos. Sci., 76, 27152738, https://doi.org/10.1175/JAS-D-18-0337.1.

    • Search Google Scholar
    • Export Citation
  • Wilhelm, J., T. R. Akylas, G. Bölöni, J. Wei, B. Ribstein, R. Klein, and U. Achatz, 2018: Interactions between mesoscale and submesoscale gravitywaves and their efficient representation in mesoscale-resolving models. J. Atmos. Sci., 75, 22572280, https://doi.org/10.1175/JAS-D-17-0289.1.

    • Search Google Scholar
    • Export Citation
  • Williams, P. D., and Coauthors, 2017: A census of atmospheric variability from seconds to decades. Geophys. Res. Lett., 44, 11 20111 211, https://doi.org/10.1002/2017GL075483.

    • Search Google Scholar
    • Export Citation
  • Williamson, J. H., 1980: Low-storage Runge-Kutta schemes. J. Comput. Phys., 35, 4856, https://doi.org/10.1016/0021-9991(80)90033-9.

  • Xie, J., M. Zhang, Z. Xie, H. Liu, Z. Chai, J. He, and H. Zhang, 2020: An orographic-drag parametrization scheme including orographic anisotropy for all flow directions. J. Adv. Model. Earth Syst., 12, e2019MS001921, https://doi.org/10.1029/2019MS001921.

    • Search Google Scholar
    • Export Citation
  • Xie, J., M. Zhang, Q. Zeng, Z. Xie, H. Liu, Z. Chai, J. He, and H. Zhang, 2021: Implementation of an orographic drag scheme considering orographic anisotropy in all flow directions in the Earth System Model CAS-ESM 2.0. J. Adv. Model. Earth Syst., 13, e2021MS002585, https://doi.org/10.1029/2021MS002585.

    • Search Google Scholar
    • Export Citation
  • Zängl, G., D. Reinert, P. Rípodas, and M. Baldauf, 2015: The ICON (ICOsahedral non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core. Quart. J. Roy. Meteor. Soc., 141, 563579, https://doi.org/10.1002/qj.2378.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Zonally averaged absolute vertical fluxes of horizontal momentum at 16-km height for (a),(b) MS-GWaM-1D and (c),(d) MS-GWaM-3D for June 1991. The solid lines represent the 50th percentile (median), and the surfaces show the ranges from the 10th to 90th percentiles. Results are separated into (a),(c) background waves and (b),(d) convectively generated waves. The black dotted lines represent the launched absolute momentum fluxes at 300 hPa during the summer solstice.

  • Fig. 2.

    Zonal mean zonal winds from composites of runs for (a)–(d) June and (e)–(h) December. The composites correspond to averages for the years 1991 through 1998. Runs from (b),(f) MS-GWaM-1D and (c),(g) MS-GWaM-3D are compared to (a),(e) simulations without any nonorographic IGW parameterization and (d),(h) the HWM2014 climatologic model.

  • Fig. 3.

    As in Fig. 2, but for zonally averaged temperatures. For comparison, we show the (d),(h) climatological model MSIS2.1.

  • Fig. 4.

    Wave action tendencies for (left) MS-GWaM-1D and (right) MS-GWaM-3D due to (a),(b) horizontal wave propagation, (c),(d) vertical wave propagation, and (e),(f) wave dissipation. The last row depicts the corresponding zonal wave drag. All tendencies are 10-day averages for 10–20 Jun 1991 at an altitude of z ≈ 25 km. Additionally, contours of the zonal wind in steps of 10 m s−1 are shown.

  • Fig. 5.

    As in Fig. 4, but for z ≈ 60 km.

  • Fig. 6.

    Wave action changes due to wave propagation in (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D for composites of (a),(b) June and (c),(d) December. The field is overlaid with vectors of wave action fluxes. They are scaled with respect to the plot dimensions and units, and the horizontal direction is enhanced by a factor of 3 for visualization. The arrow transparency represents the wave action flux strength. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments, with dashed and solid lines for negatives and nonnegatives, respectively. See Fig. 7 for the corresponding wave drag.

  • Fig. 7.

    Zonally averaged zonal gravity wave drag due to the parameterized, nonorographic IGWs. The results are averaged over composites of (a),(b) June and (c),(d) December 1991–98 and plotted for both (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments.

  • Fig. 8.

    Zonally averaged zonal gravity wave drag due to the parameterized, nonorographic IGWs. The results are averaged over composites of (a),(b) June and (c),(d) December 1991–98 and plotted for both (a),(c) MS-GWaM-1D and (b),(d) MS-GWaM-3D. The contours represent the corresponding zonally averaged zonal wind with 10 m s−1 increments.

  • Fig. 9.

    Meridional momentum fluxes in (a) MS-GWaM-1D and (b) MS-GWaM-3D for composites of June from 1991 through 1998 at an altitude of z ≈ 40 km overlaid by the corresponding zonal winds (contours; m s−1). The maps show the typically observed refraction of IGWs into the southern polar night jet and modified wind patterns.

  • Fig. B1.

    Sketch of the horizontal extent of a Lagrangian ray volume propagating in the unstructured triangular grid of the ICON model. The propagating ray volume (hatched box) may be associated with a current parent column (framed white grid cell). Given a horizontal group velocity (black arrow), we may expect it to propagate into a neighboring column (dark gray shaded grid cell). The crosses at the edges of the ray volume mark the locations for which background field data are needed to calculate the propagation of the ray volume.

  • Fig. B2.

    Performance comparison for UA-ICON at the R2B4 resolution (∼160 km) parallelized over 1920 MPI cores for runs without IGW parameterization (gray), MS-GWaM-1D (black), and MS-GWaM-3D (red). In addition to the total run time of ICON per simulated day, we show the fraction of run time associated with MS-GWaM-1D and MS-GWaM-3D. The median values for the three setups are 25 , 317, and 743 s day−1, and thus, MS-GWaM-1D and MS-GWaM-3D slow down the model by about a factor of 12 and 30, respectively.

All Time Past Year Past 30 Days
Abstract Views 217 176 0
Full Text Views 3608 3594 523
PDF Downloads 435 420 54