the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# Forcing mechanisms of the terdiurnal tide

### Friederike Lilienthal

### Christoph Jacobi

### Christoph Geißler

Using a nonlinear mechanistic global circulation model we analyze the migrating terdiurnal tide in the middle atmosphere with respect to its possible forcing mechanisms, i.e., the absorption of solar radiation in the water vapor and ozone band, nonlinear tidal interactions, and gravity wave–tide interactions. In comparison to the forcing mechanisms of diurnal and semidiurnal tides, these terdiurnal forcings are less well understood and there are contradictory opinions about their respective relevance. In our simulations we remove the wave number 3 pattern for each forcing individually and analyze the remaining tidal wind and temperature fields. We find that the direct solar forcing is dominant and explains most of the migrating terdiurnal tide's amplitude. Nonlinear interactions due to other tides or gravity waves are most important during local winter. Further analyses show that the nonlinear forcings are locally counteracting the solar forcing due to destructive interferences. Therefore, tidal amplitudes can become even larger for simulations with removed nonlinear forcings.

Atmospheric waves such as solar tides play a crucial role in the dynamics of the mesosphere and lower thermosphere (MLT) region. Tides are global-scale oscillations with periods of a solar day (24 h) or its harmonics (12 h, 8 h, etc.). They are mainly the result of absorption of solar radiation in the water vapor (troposphere) and ozone (stratosphere) region. Tidal amplitudes grow with increasing height due to the decrease of density and conservation of energy (Andrews et al., 1987; Chapman and Lindzen, 1970). In the MLT, tides can reach wind amplitudes comparable to the magnitude of the horizontal mean wind.

Due to the fact that diurnal tides (DTs) and semidiurnal tides (SDTs) usually have larger amplitudes than the harmonics of higher wave numbers or higher frequencies, they have attracted more attention in the past and are therefore relatively well understood. However, there are observations of terdiurnal tides (TDTs) showing local amplitudes comparable to those of DTs during some months of the year (Cevolani and Bonelli, 1985; Jacobi, 2012; Reddi et al., 1993; Thayaparan, 1997; Younger et al., 2002). Observations using midlatitude radar measurements show large TDT amplitudes in autumn and early winter (Beldon et al., 2006; Jacobi, 2012). Namboothiri et al. (2004) also obtained slightly larger amplitudes in winter than in summer while Thayaparan (1997) and Jacobi (2012) additionally emphasize the occurrence of TDTs during spring.

Satellite observations have been used to analyze the TDT on a global scale
(Moudden and Forbes, 2013; Pancheva et al., 2013; Smith, 2000; Yue et al., 2013). Yue et al. (2013) presented
TDT wind amplitudes from the Thermosphere Ionosphere Mesosphere Energetics
and Dynamics (TIMED) Doppler Interferometer (TIDI) of more than
16 m s^{−1} at 50^{∘} N or S above 100 km with an additional
peak in the meridional component at about 82 km between 10 and
20^{∘} N. They identified the first symmetric (3,3) mode (peaking at
8 K above the equator and at midlatitudes), using temperatures from
Sounding of the Atmosphere using Broadband Emission Radiometry (SABER). At an
altitude of 90 km, Moudden and Forbes (2013) found the largest amplitudes above
the equator during equinoxes (6–8 K), and also at 60^{∘} N during
May (7 K) and at 60^{∘} S during October (5 K) using 10 years of SABER temperature data.

Modeling studies of the TDT are mainly concerned with the analysis of forcing
mechanisms (Akmaev, 2001; Du and Ward, 2010; Huang et al., 2007; Smith and Ortland, 2001). This was motivated
by the idea that TDTs are not only the consequence of diurnal solar heating
but are additionally excited by nonlinear interactions between DTs and SDTs
(Glass and Fellous, 1975; Teitelbaum et al., 1989). The theory for these nonlinear
interactions has been outlined by Teitelbaum and Vial (1991) and later by
Beard et al. (1999). They state that the period of a child wave *P*_{3} resulting
from nonlinear interaction is linked to the periods of the parent waves *P*_{1}
and *P*_{2} through $\frac{\mathrm{1}}{{P}_{\mathrm{3}}}=\frac{\mathrm{1}}{{P}_{\mathrm{1}}}+\frac{\mathrm{1}}{{P}_{\mathrm{2}}}$. The same holds
for the wave numbers. If we consider such a pure nonlinear TDT which is only a
result of the interaction between DT and SDT, this means that the wavelength
relation between these tides must be

where *λ*_{DT}, *λ*_{SDT}, and
*λ*_{TDT} are the vertical wavelengths of the DT, SDT and TDT,
respectively. However, it should be noted that, in a real atmosphere with
unknown contributions of different forcings, this criteria is only sufficient
but not necessary to prove the existence of nonlinear interactions. For
example, the wavelengths created by nonlinear interactions may not be
detected if the solar TDT is stronger and is superposed over the nonlinear
TDT. For the same reason, a weak correlation between DT and TDT amplitudes or
between SDT and TDT amplitudes is not necessarily meaningful.

Another possible excitation source is gravity wave–tidal interactions (Huang et al., 2007; Miyahara and Forbes, 1991). More recent simulations (Ribstein and Achatz, 2016) show that details of gravity wave–tidal interactions can change if more comprehensive physics is included but their analysis does not include the TDT.

Teitelbaum et al. (1989) performed the first modeling study on the nonlinear forcing of the TDT and they concluded that the nonlinear interactions and the direct solar forcing lead to comparable terdiurnal amplitudes. Smith and Ortland (2001) used a nonlinear model with specified DT and SDT fields at the lower boundary. They switched off the terdiurnal solar component on the one hand and removed the direct solar forcing of SDTs on the other hand. As a result, they found that the solar forcing is dominant at middle and high latitudes while nonlinear interactions mainly contribute at low latitudes. A similar approach was applied by Akmaev (2001). They stated that the heating due to absorption of solar radiation in the ozone region is the main source for TDTs, while a noticeable nonlinear contribution is only seen during equinoxes. Huang et al. (2007) used a fully nonlinear tidal model with specified diurnal and semidiurnal thermotidal heating. In this model, the occurrence of TDT amplitudes was only possible due to nonlinear interactions, and they were significant in the MLT. Another model study about TDT forcing mechanisms was performed by Du and Ward (2010). They analyzed model output from the extended Canadian Middle Atmosphere Model (CMAM) with self-consistent tides due to radiative heating, convective processes and latent heat release. They performed a correlation analysis of DTs and SDTs with TDTs on a seasonal and short-term scale. They concluded that nonlinear interactions are unlikely to be the source of the migrating TDT and that solar heating is the major source. However, Du and Ward (2010) do not exclude the possibility of nonlinear interactions. They suggest a Hough mode decomposition of the TDT, similar to the analysis of Smith and Ortland (2001). This procedure allows the conclusion to which degree a local forcing actually results in a propagating tidal mode.

To summarize, there are only few modeling studies which address the forcing mechanisms of TDTs, and they do not provide a consistent perspective. Nonlinear interactions seem to play a tangible role in TDT forcing but to what extent is heavily under debate. To shed more light on this matter we have used a nonlinear global circulation model to explore this issue. To this end we performed model simulations with simultaneous nonlinear and solar terdiurnal forcing. Additional model experiments were undertaken, each with one of the forcing mechanisms switched off, in order to analyze TDT amplitudes due to each forcing, separately.

The paper is arranged as follows: the model and the numerical experiments are described in Sect. 2. Section 3 presents the results of the simulations, starting with an overview on the climatology of the reference TDT in the model. The second part of this section describes the TDTs that are obtained when certain forcings are removed. Finally, in Sect. 4 the results from Sect. 3 are discussed and summarized.

We use the nonlinear Middle and Upper Atmosphere Model (MUAM) to investigate the forcing mechanisms of tides with wave number 3. MUAM is a 3-D mechanistic model based on the COMMA-LIM (Cologne Model of the Middle Atmosphere – Leipzig Institute for Meteorology) model, which is described in detail by Fröhlich et al. (2003a, b). The more recent version of the model, MUAM, is documented by Pogoreltsev (2007); Pogoreltsev et al. (2007) and Lilienthal et al. (2017). MUAM extends from the surface (1000 hPa) to the lower thermosphere while the zonal mean temperatures in the lower 30 km (i.e., at the lower boundary and 10 height levels above) are nudged towards monthly mean ERA-Interim reanalyses of zonal mean temperature (Dee et al., 2011; ERA-Interim, 2018). Note that this only influences the zonal mean, while waves can still develop unaffected by the nudging. The background winds can freely develop in the model and are only indirectly influenced via the zonal mean temperature nudging. In the present version, there is no additional lower boundary forcing. We perform ensemble simulations for each experiment by using 11 different years (2000–2010) of monthly mean reanalysis input data, e.g., our results for January are the average of 11 simulations, nudged with 11 different years of January reanalysis data. In contrast to MUAM model experiments performed by Pogoreltsev et al. (2007) or by Jacobi et al. (2015), stationary planetary waves at the lower boundary are not explicitly forced for these model experiments in order to avoid coupling between stationary planetary waves and tides. This is important because an additional secondary coupling with planetary waves leads to a more complex situation with a more complicated quantification of the individual forcing effects.

The model has a horizontal resolution of 5×5.625^{∘} and a
vertical resolution of 2.842 km in logarithmic pressure height with a
constant scale height of *H*=7 km.

Gravity waves are calculated by an updated Lindzen-type parameterization
(Jakobs et al., 1986; Lindzen, 1981) as described by Fröhlich et al. (2003b) and
Jacobi et al. (2006). Due to the fact that this parameterization does not
account for ionospheric effects, it is coupled with a modified
parameterization after Yiǧit et al. (2008), connected via the eddy diffusion
coefficient which is calculated in the Lindzen scheme and then transferred to
the Yiǧit scheme. Gravity waves with phase speeds of 5 to
30 m s^{−1} are handled by the linear Lindzen-type scheme while the
Yiǧit scheme is restricted to phase speeds of 35 to 105 m s^{−1}.
Therefore, the intrinsic phase speeds of the waves in the Yiǧit scheme
are larger than those in the Lindzen-type scheme, so that they reach their
breaking levels at higher altitudes where the amplitude is larger. As a
result, the Lindzen-type parameterization essentially affects the
stratosphere and mesosphere and the Yiǧit parameterization mainly takes
effect in the thermosphere. Overlaps between both parameterizations are small
and the forcing terms due to gravity waves are summed in the tendency
equation of the model. Further parameterizations of solar and infrared
radiation as well as several ionospheric effects such as Rayleigh friction,
Lorentz force and ion drag are included.

MUAM experiments analyzing TDTs have been performed by Fytterer et al. (2014) who compared the simulated TDT wind shear with global lower ionospheric sporadic E occurrence rates. Additionally, Krug et al. (2015) presented a seasonal climatology of the migrating TDTs based on MUAM simulations.

In the configuration used here, the model incorporates a spin-up of 120 model days. Within that time, zonal mean heating rates (no tides) are building up a background climatology. In the subsequent 90 model days, heating rates are allowed to be zonally variable and tides start to propagate, gradually increasing in time. The heating rates are fully introduced after model day 154. In this model version, the sun's zenith angle does not account for day to day variations and refers to the first day of the respective month. The solar elevation angle, however, includes a diurnal cycle to account for tidal forcing. The last 30 model days are analyzed and presented here. They represent the mean state of the respective months with an equilibrium of background winds and temperature. Tidal amplitudes remain almost constant and show only small day-to-day variations. Note that the nudging in the troposphere and lower stratosphere is still active during that period and the model is not running completely freely at any time. This, however, does not influence the tides because the nudging only influences the zonal mean temperature. The background climatology for zonal wind, meridional wind, and temperature during solstice (January) and equinox (April) conditions is given in Fig. 1 (for details see Sect. 3.1). This simulation does not include any modifications of the tides and therefore serves as a reference, named REF in the following (see also Table 1).

Within the model there are three mechanisms that may excite TDTs: solar heating, nonlinear interactions between tides, and gravity wave–tidal interactions. The first, the diurnal variation of solar heating rates, creates atmospheric tides self-consistently. This mechanism is known to be the most important factor for the forcing of DTs and SDTs (Andrews et al., 1987). The second mechanism is related to nonlinear interactions between different tides. Following Beard et al. (1999), the interaction between a DT and a SDT can lead to the forcing of a TDT. The last source included in MUAM are gravity waves. Miyahara and Forbes (1991) have shown that an interaction between gravity waves and the DT can excite a TDT. Trinh et al. (2018) observed a longitudinal variation of gravity wave activity in the tropical MLT region that may also be caused by gravity wave–tidal interaction.

In order to separate these different mechanisms we analyze the wave number 3 component of the respective forcing and perform separate model runs in which one of the terdiurnal forcing mechanisms is removed at each model time step for each latitude/altitude. We do not consider the temporal dimension for this analysis because wave number spectra prove that TDTs in the model are strongest for wave number 3 (migrating TDTs) and negligible for other wave numbers (non-migrating TDTs, not shown here). This is because non-migrating tides are usually excited by orographic sources, latent heat release, or other geographically fixed effects (Andrews et al., 1987). Note that atmospheric gases such as water vapor or ozone are only included as zonal means which is different from other versions of MUAM (Ermakova et al., 2017). Therefore, we usually refer to the migrating TDT here. The following results are obtained from five ensemble simulations in total, eliminating each forcings separately (NO_SOL, NO_NLIN, and NO_GW), allowing all forcings (REF) and eliminating all forcings (CTRL). An overview is given in Table 1.

Note that the background (monthly mean zonal mean) circulation is not significantly altered when TDT forcings are removed (not shown here). Differences amount to not more than the actual standard deviations in the REF simulation (Fig. 1). Therefore, the influence of a removed wave number 3 forcing is comparable to the year-to-year variation of the background state and propagation conditions for tides remain similar.

The parameterization of solar heating in the middle atmosphere is calculated
following Strobel (1978). It considers heating due to the most important
gases such as water vapor, carbon dioxide, ozone, oxygen, and nitrogen.
Following this, the zonal mean ozone fields up to 50 km altitude are taken from the
Stratosphere-troposphere Processes And their Role in Climate project
(Randel and Wu, 2007; SPARC, 2018). Above 50 km, the ozone mixing ratio
decreases exponentially. The second ozone maximum near 90 km is not
included. In contrast to Jacobi et al. (2015), we restrict our simulations to
ozone data of the year 2005 because we do not intend to perform a trend
analysis. The volume mixing ratio for carbon dioxide has been chosen
according to measurements from Mauna Loa Observatory, also for the year
2005 (NOAA ESRL Global Monitoring Division, 2018; Thoning et al., 1989). Chemical
heating due to recombination of O and O_{2} (Riese et al., 1994) and heating
due to extreme ultra violet radiation (EUV) are included. This is described
in more detail by Fröhlich et al. (2003a).

In the NO_SOL simulation, the total heating rate of all heating contributions is analyzed using a Fourier transform to separate the tidal components. For the analysis of the forcing mechanism we subtract the wave number 3 amplitude from the total heating for each time step and each latitude/altitude, separately. The result of this simulation is a wave number 3 tide that is only due to nonlinear interactions and gravity wave effects.

In order to separate the nonlinear forcing we modify the nonlinear terms in the tendency equations of the model (Jakobs et al., 1986), i.e., in the advection terms in the zonal (Eq. 2) and meridional (Eq. 3) momentum equations as well as temperature advection (Eq. 4):

where ** v** is the wind vector,

*u*and

*v*are the horizontal wind components,

*w*is the vertical wind component, and

*T*is the temperature.

*a*is Earth's radius,

*ϕ*,

*λ*, and

*z*are latitude, longitude and altitude, respectively, and

*ρ*

_{0}is the reference density at a given height

*z*. Additionally, the adiabatic contribution included in the temperature equation in principle has to be taken into consideration because it includes nonlinear coupling:

with *R* as the gas constant for dry air, *m*^{′} the ratio of molecular weights
at the respective altitude and at 1000 hPa and *c*_{p} is the specific heat
at constant pressure.

Linearizing these equations, i.e., $T\approx \stackrel{\mathrm{\u203e}}{T}+{T}^{\prime}$, $w\approx \stackrel{\mathrm{\u203e}}{w}+{w}^{\prime}$, etc., results in a separation of purely nonlinear (wave–wave) interactions, wave–background interactions and pure background processes. For example, the adiabatic term from Eq. (5) may be written as follows:

and the terms on the right-hand side of Eqs. (2)–(4) are
treated similarly. The last term in the bracket of Eq. (6) describes nonlinear wave–wave interaction. From these terms of wave–wave
interactions we removed the *k*=3 amplitudes analogous to the modification of
the solar heating terms in the NO_SOL simulation. Removing the nonlinear
interactions will result in a combination of solar and gravity wave driven
TDT (Run NO_NLIN).

The simulations NO_SOL and NO_NLIN are very similar to the approach presented by Akmaev (2001) and Smith and Ortland (2001). Additionally, we consider gravity waves for the generation of TDTs. The contributions of both gravity wave routines (the Lindzen-type and the modified Yiǧit parameterization) to the tendency terms can be simply summed up. The total acceleration of the mean flow due to gravity waves is finally subject to a Fourier filtering of wave number 3, similar to the one for the heating rates and the nonlinear terms. As a result, TDTs of solar and nonlinear origin are remaining (NO_GW simulation).

As a control simulation (CTRL), the wave number 3 component of the solar, nonlinear, and gravity wave forcings are removed simultaneously. This is done in order to test the degree to which all sources of TDTs are captured, and whether the model produces further TDTs of either numerical or physical origin.

In the following analysis, we focus on the months January and April to show solstice and equinox conditions. During this time, the TDT in MUAM is most prominent. Results for July and October are similar and therefore they are not shown here.

## 3.1 Reference simulation: TDT climatology

The REF simulation includes solar, nonlinear and gravity wave forcing for all wave numbers. Therefore, it serves as a reference for all the experiments. The following results are given as a mean of the 11 ensemble members, owing to the nudging of reanalysis data for the years 2000–2010 (color shading) with the respective standard deviations (contour lines).

In Fig. 1 we provide a background climatology of the MUAM zonal mean circulation for solstice (Fig. 1a–c) and equinox (Fig. 1d–f) for the parameters zonal wind (a, d), meridional wind (b, e) and temperature (c, f). The color coding denotes the 11 year means, while the standard deviations are given as black contour lines.

Comparing the MUAM climatology with empirical climatologies such as CIRA86 (Fleming et al., 1990), the radar based GEWM (Portnyagin et al., 2004) or the satellite based UARS (Swinbank and Ortland, 2003) we find good agreement but with slightly larger westerly jets and weaker easterly jets during January in MUAM.

We notice that the model produces small year-to-year variations below
100 km in the Southern Hemisphere and south of 30^{∘} N. There, the
standard deviation *σ* is very small, mostly below *σ*(*T*)=1 K
(*σ*(*u*)=2 m s^{−1}, *σ*(*v*)=0.25 m s^{−1}). Model variations
for middle and high latitudes in the Northern Hemisphere are larger with
standard deviations up to *σ*(*T*)=6 K (*σ*(*u*)=12 m s^{−1},
*σ*(*v*)=2 m s^{−1}) during January and *σ*(*T*)=2 K
(*σ*(*u*)=2 m s^{−1}, *σ*(*v*)=0.75 m s^{−1}) for April. Maxima
of the standard deviation are located at about 60^{∘} N. These
variations have their origin in the year-to-year variability of the polar
vortex for which a range of several K, especially during winter, is
realistic. Due to the fact that MUAM assimilates the zonal mean temperature
up to 30 km altitude, this model variability represents a realistic
atmospheric variability, too.

Figures 2 and 3 show the terdiurnal component of all forcing terms that our analysis takes into account, namely solar forcing, nonlinear forcing and forcing due to gravity wave–tide interactions. All forcing terms are scaled by $\mathrm{exp}\mathit{\{}-z(\mathrm{2}H{)}^{-\mathrm{1}}\mathit{\}}$. This factor is associated with the conservation of wave energy which normalizes the wave growth with height due to the decrease in density. Therefore, the figures show the source region of tidal excitation but they do not provide any information about propagation conditions.

Figure 2 refers to thermal parameters including temperature advection (a, b), the nonlinear component of adiabatic heating (c, d), heating due to gravity waves (e, f) and direct solar heating (g, h). Note that the color scales in Fig. 2 are equal but not continuous in order to cover the magnitudes of all forcings while keeping them comparable to each other. For the thermal forcing of the TDT it can be seen that the direct solar forcing dominates in the troposphere and stratosphere. This is because of the strong absorption of solar radiation by tropospheric water vapor and stratospheric ozone. In the mesosphere (80–100 km), nonlinear effects are mainly responsible for the forcing of terdiurnal fluctuations. Due to absorption of EUV radiation, there is again some solar forcing in the lower thermosphere (Fig. 2g, h at about 120 km altitude) that is comparable to nonlinear thermal forcing (Fig. 2a, b). In this region, heating due to gravity wave effects (Fig. 2e, f) plays a major role. The nonlinear adiabatic heating effect (Fig. 2c, d) is weak everywhere compared to the other forcings and will therefore be neglected in our further considerations.

Figure 3 is similar to Fig. 2 but refers to
wind parameters. These are also scaled by the amplitude's growth rate
$\mathrm{exp}\mathit{\{}-z(\mathrm{2}H{)}^{-\mathrm{1}}\mathit{\}}$. The figure shows nonlinear zonal (a, b) and meridional
(c, d) wind advection as well as zonal (e, f) and meridional (g, h) acceleration
due to gravity waves. In the zonal wind, in the troposphere and stratosphere,
the nonlinear forcing is clearly dominating over gravity wave effects. Zonal
gravity wave forcing becomes strong above 100 km. In January, there is an
additional maximum of gravity wave induced terdiurnal forcing
(Fig. 3e) near 80 km between 30 and 60^{∘} N
which cannot be observed in April (Fig. 3f). For meridional
wind patterns, gravity wave forcing only plays a role between 80 and
100 km (Fig. 3g, h), its magnitude being comparable to
those of the advective nonlinear forcing (Fig. 3c, d). In the
stratosphere and mesosphere, nonlinear advection is the most important source
for the meridional component.

Generally, direct solar forcing is weaker during April (Fig. 2h) than during January (Fig. 2g), but most nonlinear forcings (Figs. 2a, b and 3a–d) become stronger in April and are therefore more dominant during an equinox.

As described above, the nonlinear terdiurnal forcing is a result of interactions between the migrating DT and the migrating SDT. These interactions can only take place if both DT and SDT have a considerable amplitude. To test this relation between the different harmonics, the product of DT and SDT amplitudes serves as a proxy for the terdiurnal nonlinear forcing. Due to the fact that the forcing terms in Figure 3 are scaled by the growth rate of the amplitudes, we also scaled the product of DT and SDT amplitudes to show the source region of the possible nonlinear interaction. As an example, Fig. 4 shows the results for temperature (a), zonal wind (b), and meridional wind amplitudes (c) during January.

It can be seen that the scaled product of DT and SDT amplitudes exhibits similar structures to the nonlinear terdiurnal forcing terms in Figs. 2a and 3a, c. For example, the zonal and meridional component (Fig. 4b, c), have regions of enhanced amplitude near 50 km extending from low latitudes poleward to high latitudes and with a minimum over the equator. This is in good agreement with the nonlinear zonal and meridional forcing in Fig. 3a, c. The similarities for the temperature component (Figs. 4a and 2a) are less clear but we want to emphasize that the multiplied amplitudes of DT and SDT only serve as proxy. The pure existence of an overlapping DT and SDT source region does not necessarily induce an interaction.

TDT amplitudes are presented for January (Fig. 5a–c) and
April (Fig. 5d–f). In contrast to the forcing terms, they
are not scaled. Zonal wind amplitudes become stronger in April
(Fig. 5e) compared to January (Fig. 5b)
above 110 km but this is not the case for the temperature and meridional
wind amplitude. Amplitudes at 100 km altitude reach only about 1.5 K
and 4 m s^{−1} (zonal/meridional wind). This is much smaller than
observed by radars
(Beldon et al., 2006; Jacobi, 2012; Namboothiri et al., 2004; Thayaparan, 1997) and by
satellite measurements (Moudden and Forbes, 2013; Pancheva et al., 2013; Yue et al., 2013).
They reported amplitudes of about 5–6 m s^{−1} at 90 km
(Namboothiri et al., 2004; Thayaparan, 1997) during equinoxes and local winter with
a minimum during summer. These radars are located between 40 and 50^{∘} N and in these regions, Fig. 5 also shows larger wind
amplitudes during winter and equinoxes. Beldon et al. (2006) and
Jacobi (2012) observe a maximum larger than 10 m s^{−1}
(95 km) during autumn and early winter and a smaller one during spring. The
absence of a mid-winter maximum can be explained by the location of the
radars (*>*50^{∘} N) which is north of the region with a winter maximum
as can be seen in Fig. 5b, c.

However, considering only the maxima does not give a good comparison between
seasons, and the height-latitudinal structure is more important. Especially
in temperature (Fig. 5a) and zonal wind
(Fig. 5b) we note a double-peak structure in January with
maxima at very low latitudes and a minimum at the equator. This turns into a
triple-peak structure in April (Fig. 5d, e) with maxima
slightly more poleward (30^{∘} N or S) and directly at the equator. This
structure is also visible in SABER measurements reported by
Pancheva et al. (2013) for March and December. In the meridional wind, the
structure of the TDT is not as clear in January (Fig. 5c),
with several maxima between ±60^{∘} , the strongest one appearing
near the equator. In April (Fig. 5f), it has four distinct
peaks with maxima at low and midlatitudes but, in contrast to temperature and
zonal wind, a minimum at the equator. This pattern of superposed maxima on
minima and vice versa between the zonal and meridional wind components is
expected from the wave structure itself.

The standard deviation of tidal amplitudes is relatively small, not more than
10 *%* of the total amplitude. Thus, our results prove to be robust in
structure and strength.

The TDT phases are shown in Fig. 6. At each latitude, the corresponding vertical wavelength can be obtained from the vertical phase gradient. The wavelength is taken as the vertical distance between two points of identical phases. A full span of phases should be covered between these points, and for upward propagating waves, the phase gradient for the determination should be negative. Where the amplitude is large, vertical wavelengths turn out to be longer, i.e., the vertical phase gradients are small. Where the amplitude is small, wavelengths are shorter with larger phase gradients. Thayaparan (1997), Namboothiri et al. (2004) and Jacobi (2012) report a similar relationship with vertical wavelengths being short in summer when the amplitude minimizes. Typically, the wavelengths in Fig. 6 reach 100 km and more. In January (Fig. 5a–c), the structure of phases appears to be more complex while in April (Fig. 5d–f) there are large areas of constant phase, especially at low latitudes.

Figure 7 presents the seasonal cycle of TDT amplitudes at an
altitude of 106 km. Results of satellite data analyses have frequently
been presented at 90 and 110 km
(Moudden and Forbes, 2013; Pancheva et al., 2013; Yue et al., 2013), and therefore we choose an altitude
between these heights. The temperature TDT at this altitude
(Fig. 7a) appears to be strongest during equinoxes near the
equator (3.0 K) and at midlatitudes (30–40^{∘} N or S). The
amplitudes in autumn (2.2 K) are larger than those in spring (1.6 K).
Further maxima are present during local winter at 30–40^{∘} N or S
(2.6 K at Northern Hemisphere and 2.3 K at southern Hemisphere).
For latitudes poleward of 50^{∘} N or S, amplitudes are much lower and peak during summer (*<*1.1 K).

The structure of MUAM temperature amplitudes is generally confirmed by SABER
measurements (Moudden and Forbes, 2013; Pancheva et al., 2013; Yue et al., 2013) who reported
maxima of about 5 K during equinoxes near the equator at 90 km
altitude. Note that this amplitude is almost twice as large as the one
obtained from our model simulations even though the altitude is smaller. For
midlatitudes, Moudden and Forbes (2013) also found maxima during northern winter
(3–4 K at 30–50^{∘} N, 90 km) but not during southern
winter. This is in contrast to the results of Pancheva et al. (2013) and
Yue et al. (2013) who found maxima during equinoxes and local winter in both
hemispheres (110 km altitude) which qualitatively agrees well with our
results at 106 km near 40^{∘} N or S.

Maxima in zonal wind (Fig. 7b) and meridional wind TDT
(Fig. 7c) are also found during local winter at midlatitudes.
They are slightly larger in the Northern Hemisphere (5.9 m s^{−1} in
both components) than in the Southern Hemisphere (4.7 m s^{−1} in both
components). During equinoxes, the maxima are smaller and located close to
the equator (zonal wind only, *<*4.0 m s^{−1}), at low latitudes
(meridional wind only, *<*4.3 m s^{−1}) and at midlatitudes
(zonal and meridional wind *<*3.8 m s^{−1}).

Zonal and meridional amplitudes at midlatitudes (40–50^{∘} N or S)
agree well with TIDI measurements (Yue et al., 2013) showing maxima during
equinoxes at both hemisphere and during southern winter. However, the
northern winter maximum is not seen in the zonal wind analysis by
Yue et al. (2013). Another meridional wind peak is reported by Yue et al. (2013)
near 30^{∘} N during July which can be found in our simulations, as
well. However, amplitudes tend to be underestimated by a factor of about
3–4.

Some differences between model results and satellite measurements may be explained by the orbit of the satellite passing high latitudes less frequently and leading to larger uncertainties at these latitudes. However, this cannot explain the large discrepancies in the magnitude of the TDT. Smaller model amplitudes may be due to processes that are not included in the simulations such as latent heat release.

## 3.2 Separating the forcings

In order to determine the effect of each individual forcing on the amplitude of the TDT we performed the simulations with different forcings switched off, as listed in Table 1.

NO_SOL represents a TDT that is only due to nonlinear and gravity wave
effects because wave number 3 direct solar heating is removed in the whole
model domain. Therefore, possible sources of this wave are nonlinear
interactions between other tides, i.e., between the DT and the SDT, and
gravity wave–tide interactions only. The resulting amplitudes and phases are
shown in Figs. 8 and 9. As
expected, the amplitudes are strongly reduced. However, they are not
completely extinguished. In all parameters there is a clear maximum at
northern midlatitudes (about 60^{∘} N) during January reaching 4 K ±0.6 K (temperature), 5 m s${}^{-\mathrm{1}}\pm \mathrm{1.2}$ m s^{−1} (zonal
wind) and 4 m s${}^{-\mathrm{1}}\pm \mathrm{0.8}$ m s^{−1} (meridional wind) above
120 km. In the zonal wind component there is a secondary maximum at about
30^{∘} N as well. During April, the maxima are shifted towards the
equator with amplitudes similar to those in January. This indicates that
secondary terdiurnal forcing is most evident during local winter as is
confirmed from the annual cycle of the NO_SOL simulation (not shown here).
TDT phases from this simulation (Fig. 9) are much more
irregular in comparison to the REF simulation (Fig. 6) and
show vertical wavelengths shorter than 50 km for those latitudes where TDT
amplitudes are strong.

The simulation NO_NLIN only includes direct solar forcing and gravity
wave–tide interactions. Therefore, it does not include nonlinear
interactions. Figure 10 shows the mean amplitude
differences between the NO_NLIN and REF ensembles where red (blue) colors
denote larger (smaller) amplitudes in NO_NLIN. Welch's *t* test was applied
and areas with *α**<*0.01 are hatched. It turns out that
decreased amplitudes are not the only consequence of the removed nonlinear
forcing since there are also areas where the amplitude has increased. This
result occurs primarily during January in all parameters. The strongest
increase of about 3 K (3 m s^{−1}) is located where the REF
amplitude reaches its maximum. There, the amplitude in the NO_NLIN
simulation is about 25 % larger compared to the REF simulation. Another
large red area is located at about 60^{∘} N at an altitude of
110 km. In this area the amplitudes in the REF simulation are relatively
small (not more than 2 K and 2 m s^{−1}), but the differences
between the simulations are similar so that the NO_NLIN amplitudes are twice
as strong as the REF amplitudes.

In April only weak enhancements of about 0.5–1.5 K (0.5–2 m s^{−1})
appear for individual grid points and these are not located in the areas of
larger amplitudes. Generally, the negative amplitude differences dominate and
areas of positive change are negligible.

We do not show the phases of the NO_NLIN simulation and the NO_GW simulation here because both of these simulations still include the solar forcing which dominates the other remaining forcing. As a result, the phases are almost identical to those shown in Fig. 6 for the REF simulation.

In order to investigate the reason for the positive difference in amplitude
it is useful to compare phase shifts Δ*ϕ* between the NO_NLIN TDT
(with solar and gravity wave forcing) and the NO_SOL TDT (with nonlinear and
gravity wave forcing). The gravity wave forcing appears in both simulations
and therefore the phase shift between the tides associated with these
simulations can be mainly attributed to the phase shift between a pure solar
wave and a pure nonlinear wave. The differences in the background wind and
therefore tidal propagation conditions between the simulations are small. For
120^{∘} *<*Δ*ϕ**<*240^{∘}destructive interference occurs and
leads to a decrease in amplitude for the case of superposition.

Figure 11 shows the amplitude differences as presented in
Fig. 10 but now scaled by the growth rate (factor
$\mathrm{exp}\mathit{\{}-z(\mathrm{2}H{)}^{-\mathrm{1}}\mathit{\}}$) to show the source of the positive amplitude
differences. Here, the hatched areas show regions of destructive interference
(120^{∘} *<*Δ*ϕ**<*240^{∘}) between the phases of NO_NLIN and
NO_SOL occur. It is clear that the red areas and the destructive
interferences match almost perfectly for both January and April conditions,
and for all parameters.

Figure 12 shows the mean amplitude differences between the
NO_GW and REF ensembles. For this simulation, positive amplitude differences
occur at several heights/latitudes, when the gravity wave–tide interaction as
a forcing of TDTs is removed. In this case destructive interference seems to
be more independent of the season as can be seen in January and April.
However, the regions where the zonal wind amplitude has increased are rather
small. This increase is most apparent around 60^{∘} N and 110 km
altitude during January. Note that this area is positive for both the
meridional and zonal wind and also appears in the NO_NLIN simulation
(Fig. 10). For the temperature and meridional wind
component we find, as in Fig. 10, that amplitudes in
regions with strong REF amplitudes are enhanced even when the wave number 3
gravity wave forcing is removed. Furthermore, the amplitude changes in
Fig. 12 are larger during April compared to January. This
is consistent with the larger TDT reference amplitudes occurring during
April. Generally, all amplitude differences are stronger for NO_GW than for
NO_NLIN.

The CTRL simulation provides a measure of TDT amplitudes due to effects that
have not been considered in the previous simulations. The presence of regions
of significant amplitude indicates that there still exist other sources in
the model. Figure 13 shows the TDT amplitudes for the CTRL
simulation. Note that the scale is different from Fig. 5
so that the smaller magnitudes can be seen. The structure of this remaining
tide is not completely irregular. However, the amplitudes are small with
maximum values below 1 K, 1.2 m s^{−1} (zonal wind) and
1.4 m s^{−1} (meridional wind). During January, maxima are located in
the Northern Hemisphere at low and midlatitudes and during April at the
equator (temperature) and at southern low and midlatitudes (wind).

The results of our REF simulation present a climatology and structure of the
TDT that generally agrees with observations and earlier model studies. MUAM
produces relatively small amplitudes for the TDT, e.g., 5 m s^{−1} for
the zonal wind component at 106 km altitude during winter or
12 m s^{−1} at an altitude of 120 km during April. In fact, it is an
ongoing issue that numerical models tend to underestimate tides, at least for
some regions or seasons (Pokhotelov et al., 2018; Smith, 2012).

In contrast to reports by Cevolani and Bonelli (1985); Reddi et al. (1993); Thayaparan (1997) or Yue et al. (2013) the TDT in our simulations does not attain the amplitude of a typical DT or SDT. However, this property of the TDT was mainly reported for short temporal scales of only few days which are not represented in MUAM.

MUAM simulations show strongest wind amplitudes at midlatitudes
(30–50^{∘} N) during winter with smaller maxima during spring and
autumn. This is in accordance with radar measurements at these latitudes
(Namboothiri et al., 2004; Thayaparan, 1997) who observed amplitudes of at
least 5 m s^{−1} during the whole year except during summer. At
slightly higher latitudes (50–60^{∘} N), the winter maxima disappear
and those near the equinoxes become more important as reported by,
e.g., Beldon et al. (2006) or Jacobi (2012). There are also agreements with
satellite analyses by Moudden and Forbes (2013), e.g., during equinoxes maxima
appear at the equator and at midlatitudes. However, Moudden and Forbes (2013)
observe those maxima to be more poleward (at about 60^{∘} N or S) than we
do (30–40^{∘} N or S in MUAM). They also find that winter maxima are
located about 30–40^{∘} N while poleward of 55^{∘} the maxima
appear during summer.

The TDT in model simulations by Smith and Ortland (2001) has wind maxima near
50^{∘} N or S in the respective winter hemisphere. This is slightly more
equatorward but generally agrees with our results. They have peak amplitudes
of about 10 m s^{−1}, about twice as large as in MUAM. In our model,
the zonal wind amplitudes at low latitudes are generally weaker than at
midlatitudes. Slightly enhanced amplitudes can be seen at low latitudes in
the summer hemisphere and during equinoxes above the equator. This structure
is similar to the TDT reported by Smith and Ortland (2001) at 97 km and by
Du and Ward (2010) at 95 km but vanishes at higher altitudes (Du and Ward, 2010).
The temperature amplitude in our model has a strong maximum during equinoxes
at the equator and at midlatitudes. This is not seen in earlier simulations
(Du and Ward, 2010) but agrees with observations
(Beldon et al., 2006; Jacobi, 2012; Moudden and Forbes, 2013).

In order to investigate the different generation mechanisms of the TDT we present their respective source regions. In addition to the methods used by, e.g., Akmaev (2001); Huang et al. (2007); Smith and Ortland (2001) or Du and Ward (2010), who focus on direct solar heating and nonlinear interactions between tides only, we also consider gravity wave–tide interactions as suggested by, e.g., Miyahara and Forbes (1991) and Huang et al. (2007). To summarize, the solar forcing is dominant in the troposphere and stratosphere, nonlinear interactions are present in the mesosphere, and the forcing due to gravity waves is only significant in the zonal component above the mesopause. This analysis, however, only allows insight into the local terdiurnal forcing, which does not necessarily result in a propagating tide. As suggested by Du and Ward (2010), a Hough mode decomposition of the forcing terms and of the actual tide could shed more light into this issue. Here, a different approach has been applied to analyze the propagating tides due to different forcing mechanisms. Similar to Akmaev (2001) and Smith and Ortland (2001), we perform further model simulation where possible forcing mechanisms are switched off individually.

Removing the direct terdiurnal solar heating leads to a significant decrease
in amplitude (see Fig. 8) and therefore we conclude that
the solar forcing is the most important and dominant TDT source amongst all
possible mechanisms. With respect to the relevance of the solar forcing, our
results generally agree with earlier simulations by
Akmaev (2001); Smith and Ortland (2001) and Du and Ward (2010). However, the amplitudes in
our simulations associated with the additional forcing mechanisms amount to
several K or m s^{−1}, respectively, at few latitudes and altitudes, which is
about one third to one half of the total amplitude. This gives rise to the
assumption that nonlinear interaction between tides and/or gravity waves
should also be considered as an important forcing. The “leftover
amplitudes”, which include nonlinear and gravity wave induced forcing,
exhibit a maximum at northern low and midlatitudes during January and April
alike and phases for this tide are much more complex than those associated
with solar heating (Figs. 8 and 9).
Huang et al. (2007) also underline the importance of nonlinear interactions but
they obtain pure nonlinear tidal amplitudes up to 15 m s^{−1} and
12 K during equinoxes near 100 km. These amplitudes are much larger
than those in MUAM. In contrast, the simulations by Smith and Ortland (2001) reveal
that nonlinear interactions are weak and only contribute at low latitudes.

Removing the nonlinear tidal interactions leads to an increase in amplitude
for some heights/latitudes during January by up to 2 K (3 m s^{−1}).
Although Smith and Ortland (2001) and Akmaev (2001) used the same procedure to
analyze the solar and nonlinear forcing contribution, they did not observe
this behavior of increased amplitudes. However, Smith et al. (2004) studied the
forcing mechanisms of the quarterdiurnal tide (period of 6 h) and they
have seen a similar feature. They conclude that the nonlinear forcing may
reduce rather than enhance the tide. This can be explained as a result of
destructive interferences between the purely solar forced tide and the
nonlinearly forced tide. Due to the destructive phase shift the waves are
counteracting each other and therefore reduce the amplitude when appearing
together.

Similar results are obtained when the terdiurnal gravity wave–tide interactions are removed but an increase in amplitude in this case is observed for both January and April conditions. Here, the zonal wind component is not affected by this positive amplitude change but temperature and meridional wind are.

This conclusion supports the results of Smith and Ortland (2001) and partly those of
Akmaev (2001) who found some minor nonlinear contributions but assume
the solar forcing to be a major source. While Smith and Ortland (2001) also obtain
largest nonlinear contribution at low and middle latitudes,
Akmaev (2001) point out that nonlinear interactions take place during
equinoxes. However, Akmaev (2001) only analyzed a latitude of
44^{∘} N where amplitudes seem to maximize during equinoxes and
therefore one may conclude that nonlinear interactions generally come into
play where the TDT is large. Therefore, we cannot agree with Du and Ward (2010)
who concluded that nonlinear interactions are negligible. However, we did not
perform a correlation analysis between DTs, SDTs, and TDTs and therefore we
cannot directly compare with their results. Furthermore, our simulations do
not agree with Huang et al. (2007) who obtain very large wind amplitudes over
15 m s^{−1} and temperature amplitudes over 10 K in the MLT region
for TDTs due to nonlinear interactions only. However, they also find
nonlinear amplitude maxima during equinoxes at low and middle latitudes which
is in agreement with our results.

Finally, a control simulation (CTRL) tested the TDT amplitude when all three
forcings considered here are removed simultaneously to check whether there is
a remaining weak forcing that has not yet been considered. Amplitudes for
that simulation are relatively small (*<*0.6 K and *<*1.5 m s^{−1}) but
have a clear structure with maxima at 50^{∘} N or S during local winter.
Rind et al. (2014) have noted that numerical noise can produce regular
signatures like a quasi-biennial oscillation. Therefore, noise cannot be
excluded as a tidal source in the CTRL simulation. Another reasonable TDT
source in our model could be originating from the thermospheric
parameterizations, which include some nonlinear terms. These sources,
however, are likely to be dependent on the model used and it is not likely
that the remaining amplitudes in Fig. 13 have a real
meteorological meaning.

In the future, it would be interesting to analyze non-migrating tides, as well. To do this, we would need to include additional sources such as latent heat release or 3-D ozone and water vapor (Ermakova et al., 2017). As we have seen, gravity waves are a crucial parameter for tidal forcing and they also have a large influence on the background circulation of the middle atmosphere. Therefore, the coupling of two different gravity wave parameterization is going to be replaced by an original whole atmosphere scheme after Yiǧit et al. (2008).

The MUAM model code can be obtained from the corresponding author on request.

FL designed and performed the MUAM model runs. CJ together with FL drafted the first version of the text. CJ and CG contributed to the analysis and interpretation of the results.

The authors declare that they have no conflict of interest.

This research has been funded by Deutsche
Forschungsgemeinschaft under grant JA 836/30-1. SPARC global ozone fields
were provided by William J. Randel (NCAR) through
ftp://sparc-ftp1.ceda.ac.uk/sparc/ref_clim/randel/o3data/ (Randel and Wu, 2007; SPARC, 2018). Mauna Loa
carbon dioxide mixing ratios were provided by NOAA through
ftp://aftp.cmdl.noaa.gov/data/trace_gases/co2/flask/surface/ (Dee et al., 2011; ERA-Interim, 2018).
ERA-Interim data have been provided by ECMWF on
http://apps.ecmwf.int/datasets/data/interim_full_moda/?levtype1/4pl (NOAA ESRL Global Monitoring Division, 2018; Thoning et al., 1989). We
further acknowledge support from the German Research Foundation (DFG) and
Universität Leipzig within the program of Open Access
Publishing.

Edited by: Franz-Josef
Lübken

Reviewed by: two anonymous referees

Akmaev, R.: Seasonal variations of the terdiurnal tide in the mesosphere and lower thermosphere: a model study, Geophys. Res. Lett., 28, 3817–3820, https://doi.org/10.1029/2001GL013002, 2001. a, b, c, d, e, f, g, h, i, j

Andrews, D. G., Holton, J. R., and Leovy, C. B.: Middle Atmosphere Dynamics, Academic Press Inc. Ltd., London, 1987. a, b, c

Beard, A. G., Mitchell, N. J., Williams, P. J. S., and Kunitake, M.: Non-linear interactions between tides and planetary waves resulting in periodic tidal variability, J. Atmos. Sol.-Terr. Phys., 61, 363–376, https://doi.org/10.1016/S1364-6826(99)00003-6, 1999. a, b

Beldon, C., Muller, H., and Mitchell, N.: The 8-hour tide in the mesosphere and lower thermosphere over the UK, 1988–2004, J. Atmos. Sol.-Terr. Phys., 68, 655–668, https://doi.org/10.1016/j.jastp.2005.10.004, 2006. a, b, c, d, e

Cevolani, G. and Bonelli, P.: Tidal activity in the middle atmosphere, Il Nuovo Cimento C, 8, 461–490, https://doi.org/10.1007/BF02582675, 1985. a, b

Chapman, S. and Lindzen, R.: Atmospheric tides – thermal and gravitational, p. 200, ix, D. Reidel Publishing Company, Dordrecht, Holland, 1970. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a, b

Du, J. and Ward, W.: Terdiurnal tide in the extended Canadian Middle Atmospheric Model (CMAM), J. Geophys. Res., 115, D24106, https://doi.org/10.1029/2010JD014479, 2010. a, b, c, d, e, f, g, h, i, j

ERA-Interim: Monthly mean temperature and geopotential fields on pressure levels 1979-date; Eurpoean Reanalysis Interim, available at: http://apps.ecmwf.int/datasets/data/interim-full-moda/levtype%=pl, last access: 13 August 2018. a, b

Ermakova, T. S., Statnaya, I. A., Fedulina, I. N., Suvorova, E. V., and Pogoreltsev, A. I.: Three-dimensional semi-empirical climate model of water vapor distribution and its implementation to the radiation module of the middle and upper atmosphere model, Russ. Meteorol. Hydro., 42, 594–600, https://doi.org/10.3103/S1068373917090060, 2017. a, b

Fleming, E. L., Chandra, S., Barnett, J., and Corney, M.: Zonal mean temperature, pressure, zonal wind and geopotential height as functions of latitude, Adv. Space Res., 10, 11–59, https://doi.org/10.1016/0273-1177(90)90386-E, 1990. a

Fröhlich, K., Pogoreltsev, A., and Jacobi, C.: The 48 Layer COMMA-LIM Model: Model description, new Aspects, and Climatology, Rep. Inst. Meteorol. Univ. Leipzig, 161–189, available at: http://nbn-resolving.de/urn:nbn:de:bsz:15-qucosa-217766 (last access: 30 October 2018), 2003a. a, b

Fröhlich, K., Pogoreltsev, A., and Jacobi, C.: Numerical simulation of tides, Rossby and Kelvin waves with the COMMA-LIM model, Adv. Space Res., 32, 863–868, https://doi.org/10.1016/S0273-1177(03)00416-2, 2003b. a, b

Fytterer, T., Arras, C., Hoffmann, P., and Jacobi, C.: Global distribution of the migrating terdiurnal tide seen in sporadic E occurrence frequencies obtained from GPS radio occultations, Earth Planets Space, 66, 1–9, https://doi.org/10.1186/1880-5981-66-79, 2014. a

Glass, M. and Fellous, J. L.: The eight-hour (terdiurnal) component of atmospheric tides, Space Res., 15, 191–197, 1975. a

Huang, C., Zhang, S., and Yi, F.: A numerical study on amplitude characteristics of the terdiurnal tide excited by nonlinear interaction between the diurnal and semidiurnal tides, Earth Planets Space, 59, 183–191, https://doi.org/10.1186/BF03353094, 2007. a, b, c, d, e, f, g

Jacobi, C.: 6 year mean prevailing winds and tides measured by VHF meteor
radar over Collm (51.3^{∘} N, 13.0^{∘} E), J. Atmos. Sol.-Terr. Phys., 78–79,
8–18, https://doi.org/10.1016/j.jastp.2011.04.010, 2012. a, b, c, d, e, f, g, h

Jacobi, C., Fröhlich, K., and Pogoreltsev, A.: Quasi two-day-wave modulation of gravity wave flux and consequences for the planetary wave propagation in a simple circulation model, J. Atmos. Sol.-Terr. Phys., 68, 283–292, https://doi.org/10.1016/j.jastp.2005.01.017, 2006. a

Jacobi, C., Lilienthal, F., Geißler, C., and Krug, A.: Long-term
variability of mid-latitude mesosphere-lower thermosphere winds over Collm
(51^{∘} N, 13^{∘} E), J. Atmos. Sol.-Terr. Phys., 136, 174–186,
https://doi.org/10.1016/j.jastp.2015.05.006, SI: Vertical Coupling, 2015. a, b

Jakobs, H. J., Bischof, M., Ebel, A., and Speth, P.: Simulation of gravity wave effects under solstice conditions using a 3-D circulation model of the middle atmosphere, J. Atmos. Sol.-Terr. Phys., 48, 1203–1223, https://doi.org/10.1016/0021-9169(86)90040-1, 1986. a, b

Krug, A., Lilienthal, F., and Jacobi, C.: The terdiurnal tide in the MUAM circulation model, Rep. Inst. Meteorol. Univ. Leipzig, 53, 33–44, available at: http://nbn-resolving.de/urn:nbn:de:bsz:15-qucosa2-166502 (last access: 30 October 2018), 2015. a

Lilienthal, F., Jacobi, C., Schmidt, T., de la Torre, A., and Alexander, P.: On the influence of zonal gravity wave distributions on the Southern Hemisphere winter circulation, Ann. Geophys., 35, 785–798, https://doi.org/10.5194/angeo-35-785-2017, 2017. a

Lindzen, R. S.: Turbulence and stress owing to gravity wave and tidal breakdown, J. Geophys. Res.-Oceans, 86, 9707–9714, https://doi.org/10.1029/JC086iC10p09707, 1981. a

Miyahara, S. and Forbes, J. M.: Interactions between Gravity Waves and the Diurnal Tide in the Mesosphere and Lower Thermosphere, J. Met. Soc. Jap. Ser. II, 69, 523–531, https://doi.org/10.2151/jmsj1965.69.5_523, 1991. a, b, c

Moudden, Y. and Forbes, J. M.: A decade-long climatology of terdiurnal tides using TIMED/SABER observations, J. Geophys. Res.-Space, 118, 4534–4550, https://doi.org/10.1002/jgra.50273, 2013. a, b, c, d, e, f, g, h, i

Namboothiri, S., Kishore, P., Murayama, Y., and Igarashi, K.: MF radar
observations of terdiurnal tide in the mesosphere and lower thermosphere at
Wakkanai (45.4^{∘} N, 141.7^{∘} E), Japan, J. Atmos. Sol.-Terr. Phys., 66,
241–250, https://doi.org/10.1016/j.jastp.2003.09.010, 2004. a, b, c, d, e

NOAA ESRL Global Monitoring Division: Atmospheric Carbon Dioxide Dry Air Mole Fractions from quasi-continuous measurements at Mauna Loa, Hawaii, National Oceanic and Atmospheric Administration (NOAA), Earth System Research Laboratory (ESRL), Global Monitoring Division (GMD), Boulder, Colorado, USA, available at: https://www.esrl.noaa.gov/gmd/ccgg/trends/data.html (last access: 30 October 2018), compiled by Thoning, K. W., Kitzis, D. R., and Crotwell, A., Version 2015-12, updated annually, 2018. a, b

Pancheva, D., Mukhtarov, P., and Smith, A.: Climatology of the migrating terdiurnal tide (TW3) in SABER/ TIMED temperatures, J. Geophys. Res.-Space, 118, 1755–1767, https://doi.org/10.1002/jgra.50207, 2013. a, b, c, d, e, f

Pogoreltsev, A. I.: Generation of normal atmospheric modes by stratospheric vacillations, Izvestiya, Atmospheric and Oceanic Physics, 43, 423–435, https://doi.org/10.1134/S0001433807040044, 2007. a

Pogoreltsev, A. I., Vlasov, A. A., Fröhlich, K., and Jacobi, C.: Planetary waves in coupling the lower and upper atmosphere, J. Atmos. Sol.-Terr. Phys., 69, 2083–2101, https://doi.org/10.1016/j.jastp.2007.05.014, 2007. a, b

Pokhotelov, D., Becker, E., Stober, G., and Chau, J. L.: Seasonal variability of atmospheric tides in the mesosphere and lower thermosphere: meteor radar data and simulations, Ann. Geophys., 36, 825–830, https://doi.org/10.5194/angeo-36-825-2018, 2018. a

Portnyagin, Y., Solovjova, T., Merzlyakov, E., Forbes, J., Palo, S., Ortland, D., Hocking, W., MacDougall, J., Thayaparan, T., Manson, A., Meek, C., Hoffmann, P., Singer, W., Mitchell, N., Pancheva, D., Igarashi, K., Murayama, Y., Jacobi, C., Kuerschner, D., Fahrutdinova, A., Korotyshkin, D., Clark, R., Taylor, M., Franke, S., Fritts, D., Tsuda, T., Nakamura, T., Gurubaran, S., Rajaram, R., Vincent, R., Kovalam, S., Batista, P., Poole, G., Malinga, S., Fraser, G., Murphy, D., Riggin, D., Aso, T., and Tsutsumi, M.: Mesosphere/lower thermosphere prevailing wind model, Adv. Space Res., 34, 1755–1762, https://doi.org/10.1016/j.asr.2003.04.058, 2004. a

Randel, W. J. and Wu, F.: A stratospheric ozone profile data set for 1979–2005: Variability, trends, and comparisons with column ozone data, J. Geophys. Res.-Atmos., 112, d06313, https://doi.org/10.1029/2006JD007339, 2007. a, b

Reddi, C. R., Rajeev, K., and Geetha, R.: Tidal winds in the
radio-meteor region over Trivandrum (8.5^{∘} N, 77^{∘} E), J. Atmos. Terr.
Phys., 55, 1219–1231, https://doi.org/10.1016/0021-9169(93)90049-5, 1993. a, b

Ribstein, B. and Achatz, U.: 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, 121, 8936–8950, https://doi.org/10.1002/2016JA022478, 2016. a

Riese, M., Offermann, D., and Brasseur, G.: Energy released by recombination of atomic oxygen and related species at mesopause heights, J. Geophys. Res.-Atmos., 99, 14585–14593, https://doi.org/10.1029/94JD00356, 1994. a

Rind, D., Jonas, J., Balachandran, K. N., Schmidt, G. A., and Lean, J.: The QBO in two GISS global climate models: 1. Generation of the QBO, J. Geophys. Res.-Atmos., 119, 8798–8824, https://doi.org/10.1002/2014JD021678, 2014. a

Smith, A. and Ortland, D.: Modeling and Analysis of the Structure and Generation of the Terdiurnal Tide, J. Atmos. Sci., 58, 3116–3134, https://doi.org/10.1175/1520-0469(2001)058<3116:MAAOTS>2.0.CO;2, 2001. a, b, c, d, e, f, g, h, i, j, k, l, m

Smith, A., Pancheva, D., and Mitchell, N.: Observations and modeling of the 6-hour tide in the upper mesosphere, J. Geophys. Res., 109, D10105, https://doi.org/10.1029/2003jd004421, 2004. a

Smith, A. K.: Structure of the terdiurnal tide at 95 km, Geophys. Res, Lett., 27, 177–180, https://doi.org/10.1029/1999GL010843, 2000. a

Smith, A. K.: Global Dynamics of the MLT, Surv. Geophys., 33, 1177–1230, https://doi.org/10.1007/s10712-012-9196-9, 2012. a

SPARC: Monthly mean global ozone data set 1979–2005 from Stratosphere-troposphere Processes And their Role in Climate, available at: https://www.sparc-climate.org/data-centre/data-access/referen% ce-climatology/randals-climatologies/updated-ozone-trends, last access: 13 August 2018. a, b

Strobel, D. F.: Parameterization of the atmospheric heating rate from 15 to
120 km due to O_{2} and O_{3} absorption of solar radiation, J. Geophys. Res.-Oceans, 83, 6225–6230, https://doi.org/10.1029/JC083iC12p06225, 1978. a

Swinbank, R. and Ortland, D. A.: Compilation of wind data for the Upper Atmosphere Research Satellite (UARS) Reference Atmosphere Project, J. Geophys. Res.-Atmos., 108, 4615, https://doi.org/10.1029/2002JD003135, 2003. a

Teitelbaum, H. and Vial, F.: On tidal variability induced by nonlinear interaction with planetary waves, J. Geophys. Res.-Space, 96, 14169–14178, https://doi.org/10.1029/91JA01019, 1991. a

Teitelbaum, H., Vial, F., Manson, A., Giraldez, R., and Massebeuf, M.: Non-linear interaction between the diurnal and semidiurnal tides: terdiurnal and diurnal secondary waves, J. Atmos. Terr. Phys., 51, 627–634, https://doi.org/10.1016/0021-9169(89)90061-5, 1989. a, b

Thayaparan, T.: The terdiurnal tide in the mesosphere and lower thermosphere
over London, Canada (43^{∘} N, 81^{∘} W), J. Geophys. Res.-Atmos.,
102, 21695–21708, https://doi.org/10.1029/97JD01839, 1997. a, b, c, d, e, f, g

Thoning, K. W., Tans, P. P., and Komhyr, W. D.: Atmospheric carbon dioxide at Mauna Loa Observatory: 2. Analysis of the NOAA GMCC data, 1974–1985, J. Geophys. Res.-Atmos., 94, 8549–8565, https://doi.org/10.1029/JD094iD06p08549, 1989. a, b

Trinh, Q. T., Ern, M., Doornbos, E., Preusse, P., and Riese, M.: Satellite observations of middle atmosphere–thermosphere vertical coupling by gravity waves, Ann. Geophys., 36, 425–444, https://doi.org/10.5194/angeo-36-425-2018, 2018. a

Yiǧit, W., Aylward, A., and Medvedev, A.: Parameterization of the effects of vertically propagating gravity waves for thermosphere general circulation models: Sensitivity study, J. Geophys. Res., 113, D19106, https://doi.org/10.1029/2008JD010135, 2008. a, b

Younger, P. T., Pancheva, D., Middleton, H. R., and Mitchell, N. J.: The 8-hour tide in the Arctic mesosphere and lower thermosphere, J. Geophys. Res.-Space, 107, 1420, https://doi.org/10.1029/2001JA005086, 2002. a

Yue, J., Xu, J., Chang, L. C., Wu, Q., Liu, H.-L., Lu, X., and Russell, J.: Global structure and seasonal variability of the migrating terdiurnal tide in the mesosphere and lower thermosphere, J. Atmos. Sol.-Terr. Phys., 105–106, 191–198, https://doi.org/10.1016/j.jastp.2013.10.010, 2013. a, b, c, d, e, f, g, h, i, j