Ozone–gravity wave interaction in the upper stratosphere/lower mesosphere
The increase in amplitudes of upward propagating gravity waves (GWs) with height due to decreasing density is usually described by exponential growth. Recent measurements show some evidence that the upper stratospheric/lower mesospheric gravity wave potential energy density (GWPED) increases more strongly during the daytime than during the nighttime. This paper suggests that ozone–gravity wave interaction can principally produce such a phenomenon. The coupling between ozone-photochemistry and temperature is particularly strong in the upper stratosphere where the time–mean ozone mixing ratio decreases with height. Therefore, an initial ascent (or descent) of an air parcel must lead to an increase (or decrease) in ozone and in the heating rate compared to the environment, and, hence, to an amplification of the initial wave perturbation. Standard solutions of upward propagating GWs with linear ozone–temperature coupling are formulated, suggesting amplitude amplifications at a specific level during daytime of 5 % to 15 % for low-frequency GWs (periods ≥4 h), as a function of the intrinsic frequency which decreases if ozone–temperature coupling is included. Subsequently, the cumulative amplification during the upward level-by-level propagation leads to much stronger GW amplitudes at upper mesospheric altitudes, i.e., for single low-frequency GWs, up to a factor of 1.5 to 3 in the temperature perturbations and 3 to 9 in the GWPED increasing from summer low to polar latitudes. Consequently, the mean GWPED of a representative range of mesoscale GWs (horizontal wavelengths between 200 and 1100 km, vertical wavelengths between 3 and 9 km) is stronger by a factor of 1.7 to 3.4 (2 to 50 J kg−1, or 2 % to 50 % in relation to the observed order of 100 J kg−1, assuming initial GW perturbations of 1 to 2 K in the middle stratosphere). Conclusively, the identified process might be an important component in the middle atmospheric circulation, which has not been considered up to now.
Atmospheric gravity waves (GWs), with horizontal wavelengths of 100 to 2000 km, are produced in the troposphere and propagate vertically through the stratosphere and mesosphere, where gravity wave breaking processes are important drivers of the middle atmospheric circulation (e.g., Andrews et al., 1987; Fritts and Alexander, 2003). Usually, upward propagating GWs are described by sinusoidal wave perturbations in a slowly varying background flow of an exponentially growing amplitude with height due to decreasing density (, where H is the scale height). Recently, Baumgarten et al. (2017) found some evidence that the growth of the GW amplitudes between the middle stratosphere and upper mesosphere might be stronger during the daytime than during nighttime. The aim of the present paper is to examine whether ozone–gravity wave interaction can principally produce such an amplification.
Seasonal variations of gravity wave potential energy density (GWPED) have been derived based on satellite data or lidar measurements (e.g., Geller et al., 2013; Ern et al., 2004, 2018; Kaifler et al., 2015; Baumgarten et al., 2017). At summer middle and polar latitudes, the order of the monthly mean GWPED increases from approximately 1 J kg−1 in the middle stratosphere (30–40 km) to 10 J kg−1 in the lower mesosphere (50–60 km) and 100 J kg−1 in the upper mesosphere (80–90 km), with usual initial GW perturbations in the middle stratosphere in the order of about 1 to 2 K, and wave periods primarily between 4 to 10 h (e.g., Kaifler et al., 2015; Baumgarten et al., 2017, 2018; Ern et al., 2018). Generally, the GW sources in the middle stratosphere are weaker, but the relative increase in the GWPED between the middle stratosphere and upper mesosphere are much stronger in summer than in winter, including a less pronounced seasonal cycle in the upper mesosphere than in the levels below. This is primarily due to the seasonal change in critical level filtering of the GWs by the zonal wind (e.g., Kaifler et al., 2015; Ern et al., 2018), but also due to specific GWs generated by convection and propagating towards polar latitudes (Chen et al., 2019), or to additional sources of GWs in the mesosphere, independent of the GWs at lower levels (Reichert et al., 2021). Recently, model simulations with resolved GWs suggested multistep vertical coupling processes, producing such secondary GWs as a result of dissipating primary GWs, which can strongly enhance the GW amplitudes in the upper mesosphere (e.g., Becker and Vadas, 2018; Vadas et al., 2018; Vadas and Becker, 2018). However, the potential role of daytime–nighttime differences in the increase in GW amplitudes with height have been considered only very sparsely up to now.
Baumgarten et al. (2017) derived monthly means of the GWPED from full-day lidar temperature measurements at northern mid-latitudes (54∘ N, 12∘ E), and found a stronger relative increase between 35 and 40 km and between 55 and 60 km for full-day than nighttime observations during summer months, but less pronounced differences during winter. For example, for July, the GWPED at 55–60 km show values of about J m−3 (or 10 J kg−1) for full-day measurements but about 0 J m−3 for nighttime only (or 0.2 but 0.1 J m−3, if the measured temperature fluctuations are vertically filtered for vertical wavelengths Lm<15 km), where the GWPED at 35–40 km remains nearly unchanged, indicating a difference between full-day and nighttime values by a factor of about 2. Generally, measurements of the mesospheric GWPED are much more uncertain during summer than winter months (e.g., Kaifler et al., 2015; Ehard et al., 2015; Baumgarten et al., 2017), and the signal-to-noise ratio of the lidar measurements is not as good during the daytime than during nighttime (e.g., Rüfenacht et al., 2018), which can stimulate some doubt on the reliability of the daytime–nighttime differences derived from these specific measurements. In addition, taking the potential uncertainties of the analyzing methods into account (i.e., the temporal filtering methods used for the measured time series), Baumgarten et al. (2017) speculated that a change in the phase of long periodic waves (e.g., diurnal and semidiurnal tides) could change the filtering conditions for GWs. However, Baumgarten et al. (2017) conclusively assumed that the detected daytime–nighttime differences are of true geophysical origin, where an unequivocal explanation of this phenomenon remained open. Considering also that the full-day observations of Baumgarten et al. (2018) during May 2016 showed pronounced GW activity, particularly at altitudes between 42 and 50 km where the coupling between ozone and temperature is particularly strong, it seems to be worthwhile to examine whether ozone–gravity wave interaction could principally lead to such daytime–nighttime differences in the GW amplitudes. This must then also lead to a potential effect on the differences in the GWPED between polar day and polar night. The examination of the present paper is based on standard equations describing upward propagating GWs in a constant background flow, excluding other processes controlling the GWPED variability, to provide a clear understanding and quantification of the potential effect, which cannot be achieved based on observational data analysis or comprehensive model calculations alone.
The coupling of temperature and ozone is particularly strong in the upper stratosphere due to the short photochemical lifetime of ozone (e.g., Brasseur and Solomon, 1995). Linear relationships for a change in the heating rate due to a change in ozone, and a change in photochemistry due to a change in temperature, were derived from basic theory or satellite observations, and have been introduced in standard equations of stratospheric dynamics to examine the effects on the stratospheric circulation, planetary-scale wave patterns, and equatorial Kelvin waves (Dickinson, 1973; Douglass et al., 1985; Froidevaux et al., 1989; Cordero et al., 1998; Cordero and Nathan, 2000; Nathan and Cordero, 2007; Ward et al., 2000; Gabriel et al., 2011a). Large-scale ozone-dynamic coupling processes also show significant effects in numerical weather prediction (NWP) or general circulation models (GCMs) (Cariolle and Morcrette, 2006; Gabriel et al., 2007, 2011b; Gillet et al., 2009; Waugh et al., 2009; McCormack et al., 2011; Albers et al., 2013). However, possible effects of mesoscale ozone–gravity wave interaction in the upper stratosphere/lower mesosphere (USLM) have not been considered up to now.
The basic idea of the present paper can be summarized as follows: in the USLM, the time–mean ozone mixing ratio μ0(z) decreases with height (). Therefore, at a specific level in the ULSM, an ascending air parcel initially forced by an upward propagating sinusoidal GW pattern (i.e., the wave crest with vertical velocity perturbation ) must lead to an increase () by both transport (because ) and photochemistry (because the temperature-dependent ozone production increases in the case of adiabatic cooling), and, hence, in the heating rate (), comparable to the latent heat release in the troposphere in the case of condensation. Then, the induced perturbation (, where θ is potential temperature) reinforces the initial ascent, where the lapse rate () decreases ( constant), suggesting an effective ozone adiabatic lapse rate in the upper stratosphere comparable to the moist adiabatic lapse rate in the troposphere. Analogously, a descending air parcel (the wave trough where ) leads to a decrease () and a corresponding change (), reinforcing the initial descent. Overall, this process must lead to a significant amplification of the initial GW amplitude at this level, and, hence, to a successive amplification of the amplitude during the upward level-by-level propagation through the ULSM.
In Sect. 2, standard equations for GWs in a zonal mean background flow with and without linearized ozone–temperature coupling are formulated to quantify the amplitude amplification at a specific level (or altitude) and latitude. Then, in Sect. 3, the cumulative amplitude amplification during the propagation through the USLM is derived, based on an idealized approach of the upward level-by-level propagation of GWs with specific horizontal and vertical wavelengths. Section 4 concludes with a summary and discussion.
In the following section, ozone–gravity wave interaction is analyzed based on standard equations describing GWs in a background atmosphere, where the solutions are illustrated for southern summer conditions. The background is prescribed by monthly and zonal mean temperature T0, ozone μ0, and short-wave heating rate Q0 of January 2001 (Fig. 1a–c), derived from a simulation with the high-altitude general circulation and chemistry model HAMMONIA (details of the model are given by Schmidt et al., 2010). The heating rate Q0 (Fig. 1c) is primarily due to the absorption of solar radiation by ozone, and largely agrees with southern summer solar heating rates derived from satellite measurements by Gille and Lyjak (1986) but with somewhat smaller maximum values (in the order of ∼10 %). Figure 1c shows that Q0 is particularly strong in the USLM where (the dashed line in Fig. 1b indicates ). The HAMMONIA model includes 119 layers up to 250 km with increasing vertical resolution between ∼0.7 km in the middle stratosphere and ∼1.4 km in the middle mesosphere, with a horizontal resolution of 3.75∘. In the following section, this grid is used to illustrate the analytic solutions of upward-propagating GWs.
2.1 Amplification of GW amplitudes at a specific level
2.1.1 Basic equations
Following Fritts and Alexander (2003), we consider standard Eqs. (1)–(5) describing GW propagation in a background flow, with linear GW perturbations T′, θ′, u′, v′, w′, p′, and ρ′ (T′ is temperature; ( is potential temperature;, p(z) is pressure; p00=1000 hPa; z is altitude; u′, v′, and w′ are zonal, meridional, and vertical wind perturbations, respectively; and p′ and ρ′ are the perturbations in pressure and density, respectively). Additionally, we include an ozone-dependent heating rate perturbation in the potential temperature equation (Eq. 5) and Eq. (6) for the ozone perturbation μ′ with a temperature-dependent perturbation in ozone photochemistry S′(T′), where and are linear coupling parameters as a function of latitude ϕ and altitude z specified below, is background density, H∼7 km is scale height, ρ00 is a reference value at altitude z0, u0 is a zonal mean background wind, , where and denote the derivations in longitude and latitude, g is the gravity acceleration, and f is the Coriolis parameter; the background shear terms and are neglected because of the Wentzel–Kramers–Brillouin or WKB approximation:
For , the dispersion relation for gravity waves results from Eqs. (1) to (5) by introducing sinusoidal perturbations , where denotes the perturbation quantities, Xa0 the initial amplitude at altitude zs at the lower boundary of the upper stratosphere, the exponential growth of the amplitude due to decreasing density, k1 and l1 the horizontal and meridional wave number, m1<0 the vertical wave number for upward propagating GWs with and vertical wavelength Lm1, and ω1 the frequency (here, the subscript 1 denotes the solutions for ). We focus on horizontal and vertical wavelengths Lh1≥50 km and Lm1≤15 km, where is the horizontal wave number given by , therefore (. Compressibility effects due to the vertical change in background density are excluded assuming , which is valid for vertical wavelengths Lm≤30 km. Then, the dispersion relation for the intrinsic frequency is given for the frequency range , where denotes the Brunt–Vaisala frequency:
2.1.2 Ozone–temperature coupling
For specifying the parameter b, we consider the vertical ascent in the wave crest of an initial sinusoidal GW perturbation, related to an adiabatic cooling term , which leads to an initial ozone perturbation due to the induced increase via transport, and to a change in ozone photochemistry described by (for the descent in the wave trough, the formulations are analogous but with and ). In the USLM region, ozone is very short-lived and approximate in photochemical equilibrium (Brasseur and Solomon, 1995), i.e., for pure oxygen chemistry it is approximately given by
where J2(O2) and J3(O3) are photo-dissociation rates, and cm6 s−1 and cm3 s−1 are chemical reaction rates for ozone production (O+O2+M→O3+M) and ozone loss (O+O3→2O2) (Appendix C of Brasseur and Solomon, 1995; Table 2 of Schmidt et al., 2010). Accordingly, following Brasseur and Solomon (1995), a relative change in ozone O3/O3 due to a change in temperature ΔT is given by
Then, defining and introducing a total temperature change within a background flow described by , the change S′ is given by
which is the right-hand term of Eq. (6). Overall, the initial ascent leads to an increase in ozone via transport, and the related adiabatic cooling to an increase in ozone because of the induced change . Analogously, the initial descent leads to a decrease in ozone via transport and an induced change . The height-dependence of b is specified by considering that the ozone photochemistry of the USLM region is related to the spatial structure of Q0, which is characterized by a Gaussian-type height-dependence centered at the maximum of Q0 and a rapid decrease with latitude in the extratropical winter hemisphere (see Fig. 1c). Therefore, b is multiplied with the normalized factor , where Q00 is the averaged profile of Q0 over the summer hemisphere (, where hz(z)≈1 in the summer upper stratosphere at the altitude where Q0 reach maximum values). A similar approach of Gaussian-type height-dependence in ozone–temperature coupling was successfully used by Gabriel et al. (2011a) to analyze observed planetary-scale waves in the ozone distribution.
Following previous works (e.g., Cordero and Nathan, 2000; Cordero et al., 1998; Nathan and Cordero, 2007; Ward et al., 2000; Gabriel et al., 2011a), the sensitivity of the upper stratospheric heating rate to a change in ozone is approximately described by the linear approach , where is a time-independent linear function. If we assume the same sensitivity for both the slowly varying background and the mesoscale GW perturbation propagating within the background flow, and , we may write . At a specific altitude z or pressure level p(z), we consider a GW perturbation over the vertical scale of a vertical wavelength, Δz=Lm. Then, considering that with , the first-order heating rate perturbation is given by
which is the right-hand side of Eq. (5) when defining a0=τiQ0 and . Except in polar summer regions, the effect of Q′ is limited by the length of daytime (here denoted by τday) in case of large wave periods. Therefore, we set the time increment to τi=τday in the case of τi>τday, which reduces the effect of Q′ during the time period of 24 h (e.g., τi≤12 h over the Equator). Overall, reassuming an initial ascent , the induced increase in ozone at a pressure level p(z) leads to a heating rate perturbation at this level counteracting to the initial adiabatic cooling and therefore reinforcing the initial ascent. Analogously, an initial descent is reinforced by inducing a perturbation .
Note here that the use of Δz=Lm in Eq. (11) provides a suitable measure of the effect of ozone–temperature coupling on the GW amplitudes at a specific level over the vertical distance Lm. It is also possible to set a smaller vertical scale Δz<Lm leading to smaller values at a specific level, where Δz denotes, for example, the distances of a vertical grid used in a numerical model. This modification does not change the effect over the vertical distance Lm but it provides better vertical resolution when calculating the cumulative amplitude amplification during the upward level-by-level propagation, particularly in the case of small vertical wavelengths or small vertical group velocities, as described in the next subsection.
2.1.3 Amplification of GW amplitudes at a specific level
The parameterizations of Q′ and S′ provide a useful modification of the potential temperature tendency when introducing of Eq. (6) into (Eq. 5):
Here, the amplification factor 1+ab (with ab>0) describes the feedback of the GW-induced ozone perturbation to the change in potential temperature, and an ozone adiabatic lapse rate which is – in the USLM region – smaller than because of . Alternatively, we may write:
where . As with the lapse rate, is smaller than because and . If ozone–temperature coupling becomes weak, below and above the USLM region, converges to .
Analogous to the standard solution given above, we introduce sinusoidal GW perturbations of the form in Eqs. (1)–(4) and (13) (here, the subscript 2 denotes the solutions with ozone–gravity wave coupling) which leads to the modified dispersion relation:
where and .
Equation (13) provides an evident measure of the amplification of a GW amplitude at a specific altitude z or pressure level p(z). On the one hand, introducing the same initial adiabatic potential temperature perturbation , either with or without ozone–temperature coupling, leads to . Consistently, introducing the same initial perturbation leads to or . Furthermore, combining and suggests that the amplitude is stronger than by the factor :
Overall, the introduced process of ozone–temperature coupling leads to a decrease in the GW frequency and a corresponding amplification in the GW amplitude described by the factor or . Note that the vertical variations in could affect the increase in amplitude with height, particularly in the summer upper mesosphere. Therefore, is vertically averaged over the USLM region (from 30 to 0.03 hPa, or ∼25 to ∼70 km altitude) to focus on the effects of ozone–gravity wave interaction only. Moreover, note that the relation implies not only a change in amplitude but also a slight change in the relation of horizontal and vertical wave numbers described by , i.e., a slight change in the direction of upward propagating GWs which is perpendicular to the angle α of the phase lines defined by . However, as illustrated in the following section, ozone–gravity wave interaction is particularly relevant for a range of wavelengths and periods where the induced changes in α are very small (for , or wave periods τi>2 h, the change in α is less than 0.0001∘).
2.1.4 Examples of the amplification of GW amplitudes at specific levels
Figure 1d–f shows the factor 1+ab and the quotient for a GW with horizontal and vertical wavelengths Lk=500 km and Lm=5 km, and the quotient for a GW with Lk=800 km and Lm=3 km. In the first example, the factor 1+ab (Fig. 1d) contributes to the amplification of the GW amplitude at a specific level by up to 6 %–8 %, and the overall factor (Fig. 1e) by up to 8 %–12 % (including a decrease in the lapse rate of up to 3 % described by , not shown here). The second example (Fig. 1f) shows that the factor is larger in the case of larger horizontal and smaller vertical wavelengths, reaching amplifications of up to 12 %–14 % (shaded areas denote the latitudinal range where the amplification is reduced due to the length of daytime, i.e., where τi>τday).
For illustration of the induced change in ozone at a specific level (Fig. 2a–d), we assume an initial GW perturbation with exponentially growing amplitude , with an initial temperature amplitude Ta0 of 1 K at zs≈35 km (ps=6.28 hPa) increasing to ∼8 K at z≈65 km (p=0.1 hPa). In the present paper, we formulate the solutions for pressure levels p, i.e., the initial perturbation is alternatively described by , assuming . Introducing the associated perturbation in Eq. (6) leads to , and, considering and , to an initial ozone perturbation ]. For the example of the ascent () shown in Fig. 2, we set , leading to . For Lk=500 km and Lm=5 km, the contributions (related to transport; Fig. 2a) and (related to S′; Fig. 2b) sum up to a total change of to 0.5 ppm (Fig. 2c) or to 10 % (Fig. 2d) in the USLM region, where the feedback to the heating rate is particularly strong.
The related change in the heating rate at a specific level (Fig. 2e) is given by comparing Eq. (5) with and without ozone–temperature coupling. Assuming that the same initial ascent or adiabatic cooling as above leads to (, or, when introducing to (where in case of ), Fig. 2e shows that reach values of 0.15 K h−1 over the tropics and 0.25 K h−1 at southern summer polar latitudes. Then, consistent with Eq. (16), we yield where for the change in the potential temperature perturbation, i.e., changes in temperature of 0.2–0.3 K in the USLM region (Fig. 2f). In summary, analogously considering the corresponding change for the descent, we yield an increase in the amplitude of the oscillating GW pattern at a specific level by up to 5 %–10 % in ozone and 0.2–0.3 K in temperature.
For other initial wavelengths (or associated frequencies), the latitude-height dependence is very similar to those shown in Figs. 1d–f and 2, whereas the magnitude of the amplification factor becomes smaller in the case of increasing vertical and decreasing horizontal wavelengths, or decreasing frequencies, as illustrated in Fig. 3 for an altitude where reach maximum values (1.156 hPa or ≈47 km altitude). Figure 3a shows values of for wave periods of τi>2 h steadily increasing with an increasing initial period up to values between 1.14 and 1.15. This value is limited, on the one side, because of the increasing duration of nighttime with latitude towards equatorial and northern winter regions (denoted by shaded areas), and, on the other side, because of the increasing Coriolis force in southern summer middle and polar regions (i.e., because of ).
Consistently, the amplification factor increases with decreasing vertical and increasing horizontal wavelengths (Fig. 3b and c show examples for 70 and 10∘ S), where the values are limited by the length of daytime in the case of small relations denoting the conditions where τi>τday (Fig. 3c, shaded area). Figure 3 also indicates that the examples with Lk=500 km and Lm=5 km (Figs. 1e; 2) and Lk=800 km and Lm=3 km (Fig. 1f) represent scales where ozone–gravity wave interaction is particularly efficient.
Overall, Figs. 1d–f, 2 and 3 illustrate the amplification of GW amplitudes at a specific level and a specific time. As far as the GWs are continuously propagating upward through several levels where , the amplification will be successively reinforced at each level. This cumulative amplification can lead to much stronger GW amplitudes at upper mesospheric altitudes in the case with ozone–gravity wave interaction than in the case without, as demonstrated in the next subsection.
2.2 Upward propagating GWs in a background flow
2.2.1 Level-by-level amplification of GW amplitudes
In the following section, a solution for the cumulative amplification during the vertical level-by-level propagation is derived, excluding – to a first guess – other effects like small-scale diffusion, wave-breaking processes, interaction of GWs with atmospheric tides, or so-called secondary GWs. Following the Huygens principle, each point of a propagating wave front at a specific level is the source of a new wave at this level, i.e., a single upward propagating GW, which is amplified at a level zj−1, is the initial perturbation amplified at the next level zj. For illustration (Fig. 4a–c), we choose an initial GW with horizontal and vertical wavelengths Lm=500 km and Lm=5 km as above, where the vertical distance between the levels zj−1 and zj is set by the initial vertical wavelength Δz=Lm. First, we focus on polar latitudes during southern polar summer (70∘ S) with daytime conditions only. Thereafter we consider the modification for middle and equatorial latitudes where GWs with weak vertical group velocities propagate through USLM during both daytime and nighttime.
For orientation, Fig. 4a shows the profiles for Lk=500 km and Lm=5 km at 70∘ S (solid), and, for comparison, for Lm=3 (dashed) and Lm=9 km (dotted), indicating the altitude range where ozone–temperature coupling is relevant (note that the depicted distance of pressure levels approximately represents a 5 km distance in altitude). Beginning with a first level at zs≈35 km (6.28 hPa), the wave propagates through eight layers between ≈35 and ≈70 km (0.06 hPa) where the amplification of the amplitude is relevant. At each of these levels, denoted by (j=1, n; here n=8), the amplitude at zj will be amplified by the factor at zj. Starting with an exponentially growing amplitude (where we set Ta(zs)=1 K again), we yield a new amplitude at the level z1, defining a new exponentially growing amplitude . We then yield at the level z2, defining , and so forth. Finally, the amplitude at the level zn in the middle mesosphere is described by
where the product symbol denotes the multiplication with at each level . As mentioned above, the solutions are calculated on pressure levels, i.e., z represents the geopotential height, and the vertical distance Δz between the levels is given by , where is the height-dependent scale height defined by the background. Note that using a constant scale height H0=7 km instead of H(T0) only leads to second-order changes in the cumulative amplitude amplification (the sensitivity test is described below in Sect. 2.2.4), because H(T0) only varies slightly in the USLM region (between ∼7.5 km at summer stratopause altitudes and ∼6.5 km at 70 km).
Figure 4b shows the initial amplitude Ta (blue line) and the series of the successively amplified amplitudes Tμ1, Tμ2, …, Tμn (from light blue towards red line). Figure 4c shows the related series of constant relative values , , …, , starting at the level zj (solid lines) together with the previous values starting at zj−1, multiplied by the factor (dotted lines), illustrating the successively increasing growth of the amplitude during the upward level-by-level propagation. Finally, the amplitudes converge to Tμn(z) when reaching the upper mesosphere, where Tμn(z) is stronger than Ta(z) by a factor of ∼1.47. Figure 4c also shows the fitted relative increase in the amplitude (thick red line) describing the continuous change in the growth rate of the amplitude, where Tμ(z), or Tμ(p), is defined by
with weighting functions and , where p0 is the background pressure and hPa the level of the maximum of (note that the height of this maximum is slightly decreasing from pm≈0.89 hPa over the South Pole to pm≈1.3 hPa over the Equator).
For middle and equatorial latitudes, daytime–nighttime conditions are considered by setting the amplification factor to during daytime but to Fd=1 during nighttime over the vertical wave propagation distance of 1 full day. In detail, we define the parameter , where τ0=24 h and τday is the duration of daytime within 24 h at the latitude ϕ (with Lday=1 during polar summer and Lday=0 at the Equator). Further, considering the vertical group velocity (with initial frequency ωi1 and vertical wavelength m1 as a first guess), the sinusoidal wave propagation structure between the middle stratosphere and middle mesosphere is described by changing periodically between 1 and −1 over one wavelength, where z and zm are given in kilometers and cgz in kilometers per hour, and where Lcgi=1 at the level pm, or altitude zm(pm). Then, the combined parameter separates the vertical propagation distance into daytime and nighttime fractions by defining a constant value Cd=1 in the case of Ld>1 and Cd=0 in the case of Ld≤1, where the factor provides in the case of daytime and Fd=1 in the case of nighttime.
As an example, Fig. 4d shows the profile of the resulting amplification factor Fd at 10∘ S for a GW with Lk=500 km and Lm=5 km as above, with an associated vertical group velocity cgz of about 7 km per 12 h, illustrating that we define , where zj is located in the daytime region (red) but Fd(zj)=1, where zj is located in the nighttime region (blue). The indicated vertical wave propagation distance during daytime increases towards southern summer polar latitudes but decreases towards northern winter polar latitudes. Note that for vertical wavelengths examined in the present paper (Lm≤15 km), a vertical shift of the phase – as defined by the altitude zm in the definition of Lcgz – does not have a significant impact on the cumulative amplification of the GW amplitudes because of the Gaussian-type structure of the profile of , which has been verified by several test calculations with levels other than pm, or altitudes other than zm.
In the following section, the fitted profiles Tμ are used for further examinations with different horizontal and vertical wavelengths, where the vertical level-by-level amplification is calculated by using the distances Δz=ΔzH of the vertical grid of HAMMONIA instead of Δz=Lm. This includes a smaller amplification factor over the vertical distance ΔzH because of the smaller heating rate perturbation (see Eq. 11 and the related discussion). However, the resulting differences in the amplification at a specific level over the vertical distance Lm are nearly the same, except for some small differences of less than 0.5 % due to the different vertical resolution (i.e., ). Additionally, the resulting cumulative amplification in the upper mesosphere remains nearly unchanged (, where nh is the number of the HAMMONIA levels in the USLM), where small differences between Tμnh and Tμn of less than 10 % occur only at middle and equatorial latitudes in the case of small vertical wavelengths (or small vertical group velocities) when considering the vertical propagation during both daytime and nighttime described below.
2.2.2 Cumulative amplitude amplifications for representative examples
Figure 4e illustrates the dependence of the amplitude amplification on the horizontal and vertical wavelengths Lk and Lm at 70∘ S, where it is not affected by nighttime conditions. In comparison to the example of Lk=500 km and Lm=5 km leading to a cumulative amplification of ∼1.47 (solid red line), a larger vertical wavelength of Lm=9 km leads to a smaller value of ∼1.15 (dotted red line), but a smaller vertical wavelength of Lm=3 km leads to a larger value of ∼2.27 (dashed red line), because the induced increase in the ozone perturbation μ′ produces a heating rate perturbation Q′ within a shorter (in the case of Lm=9 km) or larger (in the case of Lm=3 km) time increment τi. For the same reason, the amplification is generally larger if the horizontal wavelength Lk is larger, e.g., in the case of Lk=800 km, the final amplification in the upper mesospheric amplitudes amounts to ∼1.22 for Lm=9 km (dotted purple line), ∼1.63 for Lm=5 km (solid purple line), and ∼2.56 for Lm=3 km (dashed purple line).
The related GWPED (here denoted by E) is derived following Kaifler et al. (2015):
Introducing and N=Nμ, or and N=N0, leads to the case with (Eμ) or without (Ea) ozone–gravity wave interaction. Figure 4f shows the relative amplitudes related to Fig. 4e. In the case of Lk=500 km (red lines), the final amplifications reach values of ∼1.32 for Lm=9 km (dotted), ∼2.17 for Lm=5 km (solid), and ∼5.21 for Lm=3 km (dashed), and in case of Lk=800 km (purple lines) values of ∼1.50 for Lm=9 km (dotted), ∼2.70 for Lm=5 km (solid), and ∼6.62 for Lm=3 km (dashed). Overall, these factors provide a first-order estimate of the effect of ozone–gravity wave coupling at 70∘ S during polar summer, i.e., in case of large horizontal (≥500 km) and small vertical (≤5 km) wavelengths, we find cumulative amplifications in the upper mesosphere in the order of ∼1.5 to ∼2.5 in the temperature perturbations and in the order of ∼3 to ∼7 in the related GWPED.
2.2.3 Cumulative amplitude amplifications depending on latitude
For the GW with Lk=500 km and Lm=5 km, Fig. 5 shows the latitudinal dependence of the cumulative amplifications of the temperature perturbation (indicated by , Fig. 5a) and the related GWPED (indicated by , Fig. 5b). The values decrease from and over southern summer polar latitudes towards and at lower mid-latitudes (40∘ S), and then less rapidly towards and at 20∘ N. Overall, although the amplification of the GW amplitudes decreases rapidly with the decrease in the length of daytime, it is still quite strong at mid-latitudes.
Figure 6 shows the relations (Fig. 6a) and (Fig. 6b) at upper mesospheric levels (0.01 hPa, ∼80 km) for different horizontal and vertical wavelengths as used for Fig. 4e and f. For both Lk=500 km (red) and Lk=800 km (purple), the amplifications of the temperature perturbations and of the related GWPED are strongest for Lm=3 km (dashed lines), at polar latitudes with values between 2.5 and 3 in and between 7 and 9 in , and at middle and equatorial latitudes between 1.5 and 1.8 in and between 2.4 and 3.5 in . These values decrease with increasing vertical wavelength, i.e., when changing Lm=5 km (solid lines) or Lm=9 km (dotted lines) roughly to ∼1.7 or ∼1.25 in and ∼3.0 or ∼1.5 in at polar latitudes, and roughly to ∼1.25 or ∼1.2 in and ∼1.5 or ∼1.25 in at middle and equatorial latitudes. Overall, for the mesoscale GWs with small vertical and large horizontal wavelengths, the cumulative amplifications due to ozone–gravity wave coupling leads to much stronger amplitudes at upper mesospheric altitudes during daytime than during nighttime, in the GW perturbations by a factor between ∼1.5 at summer mid-latitudes and ∼3 for polar daytime conditions, and in the GWPED by a factor between ∼3 at summer mid-latitudes and ∼9 for polar daytime conditions.
Note that vertical momentum flux terms FGW=ρ0 () can be derived from local profiles T′ if the background is known, i.e., by (k m−1) (Ern et al., 2004). Therefore, the amplification of the GW amplitudes must lead to the same amplification of the flux term FGW and, if the GWs do not break at lower levels, of the associated gravity wave drag GWD in the upper mesosphere, suggesting an important effect of ozone–gravity wave interaction on the meridional mass circulation, particularly at polar latitudes. However, more detailed investigations need extensive numerical model simulations with a spectrum of resolved GWs, which is beyond the scope of the present paper.
2.2.4 Sensitivity to varying conditions
In the following section, we estimate the sensitivity of the GW amplitude amplification on nonlinear processes and background conditions which could modulate the first-guess results described above. For example, the decrease in the frequency towards ωi2<ωi1 includes a slight decrease in the vertical group velocity towards cgz2<cgz1, which can additionally strengthen the process of amplitude amplification because the wave propagates somewhat more slowly through the ULSM region. However, this effect is at least 1 order smaller than the first-order process described above, as derived from test calculations including this effect. For example, for Lk=500 km and Lm=5 km, cgz2 is smaller than cgz1 by 15 %–20 % at southern summer polar latitudes and 5 %–10 % at middle and equatorial latitudes. Subsequently, at a specific level, the amplification factor Fd(cgz2) is stronger than Fd(cgz1) by 2 %–3 % at polar latitudes and less than 1 % at middle and equatorial latitudes. Including this change in the successive level-by-level propagation leads to a weak successive increase in the cumulative amplifications by ∼5 % at 1 hPa to ∼10 % at 0.01 hPa at polar summer latitudes, and by only ∼1 % at 1 hPa to ∼2 % at 0.01 hPa at middle and equatorial latitudes.
We also estimate the sensitivity of the amplitude amplification on the ozone background μ0, considering the observed long-term changes in upper stratospheric ozone in the order of up to −8 % per decade (e.g., Sofieva et al., 2017; WMO, 2018), and the uncertainty in the maximum of the heating rate Q0, which is smaller in the used HAMMONIA data in the order of ∼10 % compared to those derived from satellite measurements, as mentioned above. In the case of a 10 % reduction in ozone, the cumulative amplification in the upper mesospheric GW amplitudes is weaker by about 5 % for the example with Lm=5 km and 10 % for Lm=3 km (i.e., at 70∘ S, we yield a cumulative amplification of ∼1.4 to ∼2.25 instead of ∼1.5 to ∼2.5), and the related amplification of the GWPED is weaker by about 10 % for Lm=5 km and 20 % for Lm=3 km (at 70∘ S, a cumulative amplification of ∼2.7 to ∼7.2 instead of ∼3 to ∼9). Analogously, in the case of an increase in Q0 by 10 %, the cumulative amplification is stronger by 5 % or 10 % in the GW amplitudes and by 10 % or 20 % in the related GWPED amplitudes.
Another question arises about the sensitivity of the effect of ozone–gravity wave coupling to atmospheric tides or the diurnal cycle in stratospheric ozone, which are planetary-scale processes changing the background conditions for the propagation of the mesoscale GW perturbations. For example, Schranz et al. (2018) observed stronger amplitudes in upper stratospheric ozone during daytime than during nighttime in the order of 5 % (summer solstice) to 8 % (May). Such a difference would correspond to a change in the cumulative amplification of the upper mesospheric GW amplitudes or GWPED in the order of 5 % to 10 % or 10 % to 20 %, as follows from the sensitivity of the effect on the prescribed long-term change in stratospheric ozone derived above.
Baumgarten and Stober (2019) derived amplitudes of tides in the order of up to 0.5 K in the middle stratosphere (∼35 km) increasing up to 2 K at ∼50 km and ∼4 K at 70 km, which would correspond to a change in the lapse rate in the order of up to 0.1 K km−1, or in the Brunt-Vaisala frequency in the order of 1 %. As follows from Eq. (14), a change in the amplification factor due to a relative change is given by the factor . Therefore, for wavelengths Lk≥500 km and Lm≤5 km, a relative increase (decrease) of 1 % in would lead to a relative decrease (increase) in the amplification factor of up to 0.035 % at stratopause altitudes, which is much less than the effects of the changes in the vertical group velocity or in ozone described above. Moreover, even if a relative change would be much larger (10 %–50 %), it does not change the amplification factor of a specific level by more than 1 %–3 %, and, hence, the cumulative amplification of the GW amplitudes in the upper mesosphere by more than 5 %–10 %.
Assuming exponential growth of the amplitudes () between two levels, the usual approach of a constant scale height (e.g., H∼7 km) instead of a height-dependent scale height can principally lead to significant differences in the GWPED profiles (e.g., Reichert et al., 2021). For estimating the relevance of a change in H on the cumulative amplitude amplification, the solutions are also calculated for an initial GW perturbation with a prescribed scale height H0=7 km instead of , and a related vertical distance instead of (note that H(T0) varies in the USLM region between ∼7.5 km at summer stratopause altitudes and ∼6.5 km at 70 km). Compared to the values shown in Figs. 5 and 6, the cumulative amplification of the upper mesospheric GW amplitudes is weaker by about 5 % (Lm=5 km) to 10 % (Lm=3 km) over the southern summer polar latitudes and weaker by about 1 % (Lm=5 km) to 3 % (Lm=3 km) at summer mid-latitudes. Correspondingly, the related GWPED values are weaker by about 7.5 % (Lm=5 km) to 20 % (Lm=3 km) over the southern summer polar latitudes, and 1.5 % (Lm=5 km) to 5 % (Lm=3 km) at summer mid-latitudes. Overall, these differences are smaller than the first-order effect of ozone–gravity wave coupling by approximately 1 order, where the use of H(T0) instead of H0 at the levels of relevant amplification leads to somewhat stronger amplitude amplifications, particularly over the southern summer polar latitudes, because of the difference between the high background temperatures in the summer stratopause region and the low background temperatures in the summer mesosphere (see Fig. 1a).
2.2.5 Potential effect on mean GW amplitudes
In the following section, the potential effect of ozone–gravity wave interaction is estimated for an average over a representative range of 16 different mesoscale GW events (horizontal wavelengths: 200, 500, 800, and 1100 km, vertical wavelengths: 3, 5, 7, and 9 km; see, for comparison, the amplification factor as function of wavelengths shown in Fig. 3b, c). Although these settings are idealistic, the results provide a first-guess quantification of the potential effect on time–mean GWPED values usually derived from measurements, where several different GWs contribute to the analyzed temperature fluctuations derived from the detected temperature profiles.
Figure 7 illustrates both the relative and absolute changes in the resulting mean upper mesospheric GW temperature amplitudes (Fig. 7a, b) and in the mean GWPED (Fig. 7c, d). The relative increase in the mean temperature amplitude (Fig. 7a, solid red line) is stronger by a factor increasing from about 1.3 (±0.1) at summer low and middle latitudes up to 1.7 (±0.2) at summer polar latitudes (values in brackets denote 1 standard deviation). This corresponds to a stronger increase from about 7 K (±2 K) up to 17.5 K (±4.5 K) in the case of an initial GW perturbation of 1 K in the middle stratosphere (at 6.28 hPa or ≈35 km) (Fig. 7b, solid orange line), and from about 14 K (±4 K) up to 35 K (±9 K) in the case of an initial GW perturbation of 2 K (Fig. 7b, solid purple line).
The relative increase in the mean GWPED (Fig. 7c, solid red line) is stronger by a factor increasing from about 1.7 (±0.2) at summer low and middle latitudes up to 3.4 (±0.8) at summer polar latitudes. This corresponds to a stronger increase in the absolute GWPED values from about 2 J kg−1 (±0.5 J kg−1) at summer low and middle latitudes up to 12 J kg−1 (±3 J kg−1) at summer polar latitudes in the case of an initial GW perturbation of 1 K at 35 km (Fig. 7d, solid orange line), and from about 8 J kg−1 (±2 J kg−1) up to 48 J kg−1 (±0.5 J kg−1) in the case of an initial GW perturbation of 2 K (Fig. 7d, solid purple line).
In summary, we find an absolute increase in the order of 7 to 35 K in the mean GW temperature amplitudes and 2 to 50 J kg−1 in the mean GWPED values, assuming usual initial GW perturbations in the order of 1 to 2 K in the middle stratosphere, where the effect is particularly large during polar daytime conditions. Note that, assuming exponential growth with height only, this potential effect can be much larger in the case of stronger initial amplitudes in the middle stratosphere (the absolute changes in the temperature amplitudes increase linearly and those in the GWPED values quadratically with increasing initial GW perturbations at 35 km) and in specific geographical regions or time periods where primary GWs with large horizontal and small vertical wavelengths are excited (e.g., where Lk≥800 and Lm≤3 km). However, the GWs with very large amplitudes might dissipate by nonlinear wave-breaking processes before reaching the upper mesosphere.
The present paper shows that ozone–gravity wave interaction in the upper stratosphere/lower mesosphere (USLM) leads to a stronger increase in gravity wave (GW) amplitudes with height during daytime than during nighttime, particularly during polar summer. The results include information about both the amplification of the GW amplitudes at a specific level and the cumulative increase in the amplitudes during the upward level-by-level propagation of the wave from middle stratosphere to upper mesosphere.
In a first step, standard equations describing upward propagating GWs with and without linearized ozone–gravity wave coupling are formulated, where an initial sinusoidal GW perturbation in the vertical ozone transport and temperature-dependent ozone photochemistry produces a heating rate perturbation as a function of the initial intrinsic frequency, which determines the duration of the perturbation at a specific level over the distance of the initial vertical wavelength. The solution reveals an amplification of the ascending and descending perturbations of the sinusoidal GW pattern at this level, i.e., a decrease in the intrinsic frequency due to both the induced changes in the lapse rate (or Brunt–Vaisala frequency) and the positive feedback of the coupling on the initial GW perturbation, and an associated increase in the GW amplitude by a factor defined by the relation of the intrinsic frequencies without (ωi1) and with (ωi2) ozone–gravity wave coupling. This amplitude amplification is dependent on the horizontal and vertical wavelengths, Lk and Lm, where the effect is most efficient for GWs with Lk≥500 km and Lm≤5 km, or initial frequencies τi≥4 h, representing mesoscale GWs forced by cyclones or fronts, or by the orography of mountain ridges like the Rocky Mountains, Andes, or Norwegian Caledonides. For southern summer conditions, strongest amplitude amplifications at specific levels of about 5 %–15 % over the perturbation distance of one vertical wavelength are located near the stratopause, with peak values over the Equator and over summer polar latitudes.
In a second step, an analytic approach of the upward level-by-level propagation of the GW perturbations with and without ozone–gravity wave interaction reveals the cumulative amplitude amplification, where the wave is propagating upward with the vertical group velocity defined by the initial GW parameters, and where daytime–nighttime conditions at middle and equatorial latitudes are considered. Representative examples with different initial wavelengths illustrate that the successive increase in both the GW amplitudes and the related gravity wave potential energy density (GWPED) converge to much stronger amplitudes in the upper mesosphere during daytime than during nighttime. This effect strongly decreases with latitude between summer polar and mid-latitudes because of the decrease in the length of daytime, nearly constant at equatorial latitudes, and decreasing again with latitude towards insignificant values in the winter extratropics.
In summary, the strongest impact of ozone–gravity wave interaction is found for wave periods ≥4 h (related to the wavelengths Lk≥500 km and Lm≤5 km), i.e., in a range of wave periods usually observed at summer middle and polar latitudes. For prescribed single GWs with large horizontal wavelengths (500 to 800 km) and small vertical wavelengths (3 to 5 km), the upper mesospheric GW temperature amplitudes are stronger by a factor between 1.25 and 1.75 at summer low and middle latitudes and between 1.5 and 3 at summer polar latitudes, and the corresponding GWPED by a factor between 1.5 and 3.5 and between 3 and 9. For a representative range of 16 different mesoscale GW events (Lk between 200 and 1100 km, Lm between 3 and 9 km), the mean temperature amplitudes are stronger by a factor between 1.3 at summer low and middle latitudes to 1.7 at summer polar latitudes, e.g., stronger by about 7 to 17.5 K (or 14 to 35 K) in the case of an initial GW perturbation of 1 K (or 2 K) in the middle stratosphere (at ∼35 km). The corresponding relative increase in the mean GWPED is stronger by a factor between 1.7 at summer low and middle latitudes and 3.4 at summer polar latitudes, e.g., for the same example as above, stronger by about 2 to 12 J kg−1 (or 8 to 48 J kg−1). These values range in the order between 2 % and 50 % of the observed order of the mean upper-mesospheric GWPED amplitudes (100 J kg−1). These absolute differences can be larger in the case of stronger initial perturbations in the middle stratosphere, or in specific geographical regions or time periods where primary GWs with large horizontal and small vertical wavelengths (e.g., where Lk≥800 km and Lm≤3 km) are excited. However, the GWs with very large amplitudes might dissipate by nonlinear wave-breaking processes before reaching the upper mesosphere. Overall, these values result from an idealistic approach and cannot entirely explain the details of specific measurements. Nevertheless, they provide a first-guess quantification of the potential effect of ozone–gravity wave interaction on the GW amplitudes.
The variety of horizontal and vertical wavelengths used in the present paper are representative of mesoscale GWs in the USLM region. Observations suggest not only characteristic vertical wavelengths of GWs between ∼2–5 km in the lower stratosphere increasing to ∼10–30 km in the upper mesosphere, but also the existence of large vertical wavelengths greater than 10 km in the ULSM region, particularly above convection in equatorial regions or over southernmost Argentina (e.g., Alexander, 1998; McLandress et al., 2000; Fritts and Alexander, 2003; Alexander and Holton, 2014; Hocke et al., 2016; Baumgarten et al., 2018; Reichert et al., 2021). The results of the present paper suggest that the effect of ozone–gravity wave coupling decreases with increasing vertical wavelengths Lm≥9 km, but strongly increases with decreasing vertical wavelengths Lm≤5 km. The latter could lead to more pronounced GW breaking and dissipation processes in the upper stratosphere during daytime than during nighttime, and – subsequently – to more prominent GWs with larger vertical wavelengths of Lm≥5 km, which would be consistent with the observed GW characteristics at these altitudes presented by Baumgarten et al. (2018).
As mentioned in the introduction, the measurements of Baumgarten et al. (2017) show some evidence that the increase in the GWPED values with height is stronger during full-daytime than during nighttime by a factor of about 2, or, roughly assuming a 2:1 relation of daytime and nighttime (16 h daytime and 8 h nighttime) for high summer mid-latitudes, stronger during daytime than during nighttime by a factor of about 2.5. For comparison, the estimated effect of ozone–temperature coupling for these latitudes (factor of 1.7) is somewhat smaller and would lead to an increase in the nighttime GWPED in the order of ∼50 % (0.7:1.5) of the observed increase. Conclusively, although the difference derived by Baumgarten et al. (2017) might be uncertain as mentioned in the introduction, and although the approach of the present paper cannot entirely explain the details of specific local measurements during a specific time period, the comparison confirms that ozone–gravity wave interaction might be able to produce significant daytime–nighttime differences in the GW amplitudes at high summer mid-latitudes.
Current state-of-the-art general circulation models (GCMs) usually use a variety of prescribed tropospheric sources and tuning parameters in the gravity wave drag (GWD) parameterizations, forcing the middle atmospheric circulation (e.g., McLandress et al., 1998; Fritts and Alexander, 2003; Garcia et al., 2017), where the extreme low temperatures observed in the summer upper mesosphere provide an important benchmark for the quality of the upwelling branch and the associated adiabatic cooling produced by the models. Including ozone–gravity wave interaction into the GCMs might lead to a substantial improvement of the used GWDs and the associated processes driving the summer mesospheric circulation, because the related increase in the GWPED must lead to a similar increase in the vertical momentum flux term determining the GWD. However, the incorporation of ozone–gravity wave interaction into a state-of-the-art GCM using a GWD, or into a numerical model with resolved GWs, needs extensive test simulations, which is beyond the scope of the present paper.
In particular, current GCMs indicate significant changes in the time–mean circulation of the upper mesosphere due to the stratospheric ozone loss over Antarctica during southern spring and early summer via the induced changes in the GWD (Smith et al., 2010; Lossow et al., 2012; Lubi et al., 2016). Long-term changes in upper stratospheric ozone of up to −8 % per decade, derived from satellite measurements (e.g., Sofieva et al., 2017; WMO, 2018), could also affect the mesospheric circulation in the stratosphere and mesosphere by modulating the GW amplitudes and, hence, the GWD. Based on the idealized approach of the present paper, we estimate the sensitivity of the amplification of the GW amplitudes in the upper mesosphere on changes in the ozone background μ0 and the ozone-related heating rate Q0(μ0), revealing that, for horizontal and vertical wavelengths Lk≥500 km and Lm≤5 km, a change of ±10 % in μ0 or Q0 results in a change of ±10 % to ±20 % in the upper mesospheric GWPED. Conclusively, the summer mesospheric upwelling might be much more sensitive to the long-term changes in upper stratospheric ozone as has been suggested by the GCMs up to now.
In the approach of the present paper, the variations due to the diurnal cycle in stratospheric ozone and atmospheric tides are excluded to examine the potential effect of ozone–gravity wave interaction as clearly as possible, based on standard equations describing upward propagating GWs in a constant background. On the one hand, these variations can principally modulate the effect of ozone–gravity wave coupling by changing the planetary-scale background conditions for the propagation of the mesoscale GWs. Assuming – to a first order – linear modulations in the background ozone and background lapse rate according to observed diurnal or tidal variations, the sensitivity calculations of the present paper suggest that the related modulations in the amplitude amplification are smaller than the effect of ozone–gravity wave coupling by approximately 1 order. Further test calculations have shown that the use of a height-dependent scale height H(T0) instead of a constant scale height H0 at the levels of relevant amplification leads to stronger amplitude amplifications, particularly over the southern summer polar latitudes, because of the high temperatures in the stratopause region and the very low temperatures in the upper mesosphere, where the related differences are also smaller than the first-order process (e.g., in the GWPED, for vertical wavelengths between Lm=5 km and Lm=3 km, between about 7.5 % to 20 % at summer polar latitudes and less than 5 % at summer mid-latitudes).
On the other hand, short-term fluctuations in the balanced zonal and meridional winds due to atmospheric tides can principally lead to changes in the upward GW propagation characteristics, and, hence, to significant daytime–nighttime differences in the growth of the GW amplitudes with height, including nonlinear feedbacks of the propagating mesoscale GWs to the short-term balanced flow components. Further, multistep vertical coupling processes producing secondary GWs in the mesosphere could depend on daytime–nighttime conditions or tidal variations, which could also produce significant daytime–nighttime differences in the growth of the GW amplitudes with height. Considering the remarkably strong effect of ozone–gravity wave coupling suggested by the present paper, we may speculate that it significantly affects these possible changes in the GW amplitudes due to short-term fluctuations in the balanced winds or multistep vertical coupling. However, an unequivocal quantification of the effects of these processes and the involved nonlinear interactions of the daytime–nighttime differences in the GWPED needs much more investigations, e.g., based on extensive GW resolving model simulations with interactive ozone photochemistry, which is beyond the scope of the present paper.
The results of the present paper might stimulate further daytime–nighttime observations of GW activity, particularly at specific measurement sites where the GWs are usually characterized by specific horizontal and vertical wavelengths, e.g., downwind of specific mountain ridges (east of Rocky Mountains, Southern Andes or Norwegian Caledonides), which could be helpful to better understand how ozone–gravity wave coupling is operating in situ.
Background data and programs visualizing the presented analytic solutions are available upon request from the author.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The author thanks Hauke Schmidt (Max Planck Institute for Meteorology (MPI-Met), Hamburg) for providing HAMMONIA background data. Thanks also to two reviewers for critical comments.
The publication of this article was funded by the Open Access Fund of the Leibniz Association.
This paper was edited by Mathias Palm and reviewed by two anonymous referees.
Albers, J. R., McCormack, J. P., and Nathan, T. R.: Stratospheric ozone and the morphology of the northern hemisphere planetary waveguide, J. Geophys. Res.-Atmos., 118, 563–576, https://doi.org/10.1029/2012JD017937, 2013.
Alexander, M. J.: Interpretations of observed climatological patterns in stratospheric gravity wave variance, J. Geophys. Res., 103, 8627–8640, 1998.
Alexander, M. J. and Holton, J. R.: On the spectrum of vertically propagating gravity waves generated by a transient heat source, Atmos. Chem. Phys., 4, 923–932, https://doi.org/10.5194/acp-4-923-2004, 2004.
Andrews, D. G., Holton, J. R., and Leovy, C. B.: Middle Atmosphere Dynamics, Academic Press, San Diego, California, 489 pp., 1987.
Baumgarten, K. and Stober, G.: On the evaluation of the phase relation between temperature and wind tides based on ground-based measurements and reanalysis data in the middle atmosphere, Ann. Geophys., 37, 581–602, https://doi.org/10.5194/angeo-37-581-2019, 2019.
Baumgarten, K., Gerding, M., and Lübken, F.-J.: Seasonal variation of gravity wave parameters using different filter methods with daylight lidar measurements at midlatitudes, J. Geophys. Res.-Atmos., 122, 2683–2695, https://doi.org/10.1002/2016JD025916, 2017.
Baumgarten, K., Gerding, M., Baumgarten, G., and Lübken, F.-J.: Temporal variability of tidal and gravity waves during a record long 10-day continuous lidar sounding, Atmos. Chem. Phys., 18, 371–384, https://doi.org/10.5194/acp-18-371-2018, 2018.
Becker, E., and Vadas, S. L.: Secondary gravity waves in the winter mesosphere: Results from a high-resolution global circulation model, J. Geophys. Res., 123, 2605–2627, https://doi.org/10.1002/2017JD027460, 2018.
Brasseur, G. and Solomon, S.: Aeronomy of the Middle Atmosphere, D. Reidel Publishing Company, Dordrecht, Netherlands, 445 pp., 1995.
Cariolle, D. and Morcrette, J.-J.: A linearized approach to the radiative budget of the stratosphere: influence of the ozone distribution, Geophys. Res. Lett., 33, L05806, https://doi.org/10.1029/2005GL025597, 2006.
Chen, D., Strube, C., Ern, M., Preusse, P., and Riese, M.: Global analysis for periodic variations in gravity wave squared amplitudes and momentum fluxes in the middle atmosphere, Ann. Geophys., 37, 487–506, https://doi.org/10.5194/angeo-37-487-2019, 2019.
Cordero, E. C. and Nathan, T. R.: The Influence of Wave– and Zonal Mean–Ozone Feedbacks on the Quasi-biennial Oscillation, J. Atmos. Sci., 57, 3426–3442, 2000.
Cordero, E. C., Nathan, T. R., and Echols, R. S.: An analytical study of ozone feedbacks on Kelvin and Rossby-gravity waves: Effects on the QBO, J. Atmos. Sci., 55, 1051–1062, 1998.
Dickinson, R. E.: Method of parameterization for infrared cooling between the altitude of 30 and 70 km, J. Geophys. Res., 78, 4451, https://doi.org/10.1029/JC078i021p04451, 1973.
Douglass, A. R., Rood, R. B., and Stolarski, R. S.: Interpretation of Ozone Temperature Correlations 2. Analysis of SBUV Ozone Data, J. Geophys. Res., 90, 10693–10708, https://doi.org/10.1029/JD090iD06p10693, 1985.
Ehard, B., Kaifler, B., Kaifler, N., and Rapp, M.: Evaluation of methods for gravity wave extraction from middle-atmospheric lidar temperature measurements, Atmos. Meas. Tech., 8, 4645–4655, https://doi.org/10.5194/amt-8-4645-2015, 2015.
Ern, M., Preusse, P., Alexander, M. J. and Warner, C. D.: Absolute values of gravity wave momentum flux derived from satellite data, J. Geophys. Res., 109, D20103, https://doi.org/10.1029/2004JD004752, 2004.
Ern, M., Trinh, Q. T., Preusse, P., Gille, J. C., Mlynczak, M. G., Russell III, J. M., and Riese, M.: GRACILE: a comprehensive climatology of atmospheric gravity wave parameters based on satellite limb soundings, Earth Syst. Sci. Data, 10, 857–892, https://doi.org/10.5194/essd-10-857-2018, 2018.
Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1003, https://doi.org/10.1029/2001RG000106, 2003.
Froidevaux, L., Allen, M., Berman, S., and Daughton, A.: The Mean Ozone Profile and Its Temperature Sensitivity in the Upper Stratosphere and Lower Mesosphere: An Analysis of LIMS Observations, J. Geophys. Res., 94, 6389–6417, 1989.
Gabriel, A., Peters, D. H. W., Kirchner, I., and Graf, H.-F.: Effect of zonally asymmetric ozone on stratospheric temperature and planetary wave propagation, Geophys. Res. Lett., 34, L06807, https://doi.org/10.1029/2006GL028998, 2007.
Gabriel, A., Körnich, H., Lossow, S., Peters, D. H. W., Urban, J., and Murtagh, D.: Zonal asymmetries in middle atmospheric ozone and water vapour derived from Odin satellite data 2001–2010, Atmos. Chem. Phys., 11, 9865–9885, https://doi.org/10.5194/acp-11-9865-2011, 2011a.
Gabriel, A., Schmidt, H., und Peters, D. H. W.: Effects of the 11-year solar cycle on middle atmospheric stationary wave patterns in temperature, ozone, and water vapor, J. Geophys. Res., 116, D23301, https://doi.org/10.1029/2011JD015825, 2011b.
Garcia, R. R, Smith, A., Kinnison, D., de La Cámara, Á., and Murphy, D. J.: Modification of the Gravity Wave Parameterization in the Whole Atmosphere Community Climate Model: Motivation and Results, J. Atmos. Sci., 74, 275–291, https://doi.org/10.1175/JAS-D-16-0104.1, 2017.
Geller, M., Alexander, M. J., Love, P., Bacmeister, J., Ern, M., Hertzog, A., Manzini, E., Preusse, P., Sato, K., Scaife, A., and Zhou, T.: A Comparison between Gravity Wave Momentum Fluxes in Observations and Climate Models, J. Climate, 26, 6383–6405, https://doi.org/10.1175/JCLI-D-12-00545.1, 2013.
Gille, J. C. and Lyjak, L. V.: Radiative Heating and Cooling Rates in the Middle Atmosphere, J. Atmos. Sci., 43, 2215–2229, 1986.
Gillett, N. P., Scinocca, J. F., Plummer, D. A., and Reader, M. C.: Sensitivity of climate to dynamically-consistent zonal asymmetries in ozone, Geophys. Res. Lett., 36, L10809, https://doi.org/10.1029/2009GL037246, 2009.
Hocke, K., Lainer, M., Moreira, L., Hagen, J., Fernandez Vidal, S., and Schranz, F.: Atmospheric inertia-gravity waves retrieved from level-2 data of the satellite microwave limb sounder Aura/MLS, Ann. Geophys., 34, 781–788, https://doi.org/10.5194/angeo-34-781-2016, 2016.
Kaifler, B., Lübken, F.-J., Höffner, J., Morris, R. J., and Viehl, T. P.: Lidar observations of gravity wave activity in the middle atmosphere over Davis (69∘S, 78∘E), Antarctica, J. Geophys. Res.-Atmos., 120, 4506–4521, https://doi.org/10.1002/2014JD022879, 2015.
Lossow, S., McLandress, C., and Shepherd, T. G.: Influence of the Antarctic ozone hole on the polar mesopause region as simulated by the Canadian Middle Atmosphere Model, J. Atmos. Sol.-Terr. Phy., 74, 111–123, https://doi.org/10.1016/j.jastp.2011.10.010, 2012.
Lubis, S. W., Omrani, N.-E., Matthes, K., and Wahl, S.: Impact of the Antarctic Ozone Hole on the Vertical Coupling of the Stratosphere–Mesosphere–Lower Thermosphere System, J. Atmos. Sci., 73, 2509–2528, https://doi.org/10.1175/JAS-D-15-0189.1, 2016.
McCormack, J. P., Nathan, T. R., and Cordero, E. C.: The effect of zonally asymmetric ozone heating on the Northern Hemisphere winter polar stratosphere, Geophys. Res. Lett., 38, L03802, https://doi.org/10.1029/2010GL045937, 2011.
McLandress, C.: On the importance of gravity waves in the middle atmosphere and their parameterization in general circulation models, J. Atmos. Sol.-Terr. Phy., 60, 1357–1383, 1998.
McLandress, C., Alexander, M. J., and Wu, D. L.: Microwave Limb Sounder observations of gravity waves in the stratosphere: A climatology and interpretation, J. Geophys. Res., 105, 11947–11967, 2000.
Nathan, T. R. and Cordero, E. C.: An ozone-modified refractive index for vertically propagating planetary waves, J. Geophys. Res., 112, D02105, https://doi.org/10.1029/2006JD007357, 2007.
Reichert, R., Kaifler, B., Kaifler, N., Dörnbrack, A., Rapp, M., and Hormaechea, J. L., High-cadence lidar observations of middle atmospheric temperature and gravity waves at the Southern Andes hot spot, J. Geoph. Res.-Atmos., 126, e2021JD034683, doi.org/10.1029/2021JD034683, 2021.
Rüfenacht, R., Baumgarten, G., Hildebrand, J., Schranz, F., Matthias, V., Stober, G., Lübken, F.-J., and Kämpfer, N.: Intercomparison of middle-atmospheric wind in observations and models, Atmos. Meas. Tech., 11, 1971–1987, https://doi.org/10.5194/amt-11-1971-2018, 2018.
Schmidt, H., Brasseur, G. P., and Giorgetta, M. A.: Solar cycle signal in a general circulation and chemistry model with internally generated quasi-biennial oscillation, J. Geophys. Res., 115, D00I14, https://doi.org/10.1029/2009JD012542, 2010.
Schranz, F., Fernandez, S., Kämpfer, N., and Palm, M.: Diurnal variation in middle-atmospheric ozone observed by ground-based microwave radiometry at Ny-Ålesund over 1 year, Atmos. Chem. Phys., 18, 4113–4130, https://doi.org/10.5194/acp-18-4113-2018, 2018.
Smith, A., Garcia, R. R. Marsh, D. R., Kinnison, D. E., and Richter, J. H.: Simulations of the response of mesospheric circulation and temperature to the Antarctic ozone hole, Geophys. Res. Lett., 37, L22803, https://doi.org/10.1029/2010GL045255, 2010.
Sofieva, V. F., Kyrölä, E., Laine, M., Tamminen, J., Degenstein, D., Bourassa, A., Roth, C., Zawada, D., Weber, M., Rozanov, A., Rahpoe, N., Stiller, G., Laeng, A., von Clarmann, T., Walker, K. A., Sheese, P., Hubert, D., van Roozendael, M., Zehner, C., Damadeo, R., Zawodny, J., Kramarova, N., and Bhartia, P. K.: Merged SAGE II, Ozone_cci and OMPS ozone profile dataset and evaluation of ozone trends in the stratosphere, Atmos. Chem. Phys., 17, 12533–12552, https://doi.org/10.5194/acp-17-12533-2017, 2017.
Vadas, S. L. and Becker, E.: Numerical modeling of the excitation, propagation, and dissipation of primary and secondary gravity waves during wintertime at McMurdo Station in the Antarctic, J. Geophys. Res.-Atmos., 123, 9326–9369, https://doi.org/10.1029/2017JD027974, 2018.
Vadas, S. L., Zhao, J., Chu, X., and Becker, E..: The excitation of secondary gravity waves from local body forces: Theory and observation, J. Geophys. Res.-Atmos., 123, 9296–9325, https://doi.org/10.1029/2017JD027970, 2018.
Ward, W. E., Oberheide, J., Riese, M., Preusse, P., and Offermann, D.: Pla netary wave two signatures in CRISTA 2 ozone and temperature data, in Atmospheric Science Across the Stratopause, edited by: Siskind, D. E., Eckermann, S. D., and Summers, M. E., Geophys. Monogr., 103, 319–325, 2000.
Waugh, D. W., Oman, L., Newman, P. A., Stolarski, R. S., Pawson, S., Nielsen, J. E., and Perlwitz, J.: Effect of zonal asymmetries in stratospheric ozone on simulated Southern Hemisphere climate trends, Geophys. Res. Lett., 36, L18701, https://doi.org/10.1029/2009GL040419, 2009.
WMO (World Meteorological Organization): Scientific Assessment of Ozone Depletion: 2018, Global Ozone Research and Monitoring Project, Report No. 58, Geneva, Switzerland, 588 pp., 2018.