Research article 27 Nov 2020
Research article  27 Nov 2020
Gravitywaveperturbed wind shears derived from SABER temperature observations
 ^{1}Henan Engineering Laboratory for Big Data Statistical Analysis and Optimal Control, School of Mathematics and Information Sciences, Henan Normal University, Xinxiang, 453007, China
 ^{2}State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing, 100190, China
 ^{3} School of Astronomy and Space Science, University of the Chinese Academy of Science, Beijing, 100049, China
 ^{4}Catholic University of America, Washington, DC 20064, USA
 ^{5}Atmospheric and Planetary Sciences, Hampton University, Hampton, VA 23668, USA
 ^{6}High Altitude Observatory, National Center for Atmospheric Research, Boulder, CO 80301, USA
 ^{1}Henan Engineering Laboratory for Big Data Statistical Analysis and Optimal Control, School of Mathematics and Information Sciences, Henan Normal University, Xinxiang, 453007, China
 ^{2}State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing, 100190, China
 ^{3} School of Astronomy and Space Science, University of the Chinese Academy of Science, Beijing, 100049, China
 ^{4}Catholic University of America, Washington, DC 20064, USA
 ^{5}Atmospheric and Planetary Sciences, Hampton University, Hampton, VA 23668, USA
 ^{6}High Altitude Observatory, National Center for Atmospheric Research, Boulder, CO 80301, USA
Correspondence: Jiyao Xu (xujy@nssc.ac.cn)
Hide author detailsCorrespondence: Jiyao Xu (xujy@nssc.ac.cn)
Large wind shears around the mesopause region play an important role in atmospheric neutral dynamics and ionospheric electrodynamics. Based on previous observations using sounding rockets, lidars, radars, and model simulations, large shears are mainly attributed to gravity waves (GWs) and modulated by tides (Liu, 2017). Based on the dispersion and polarization relations of linear GWs and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperature data from 2002 to 2019, a method of deriving GWperturbed wind shears is proposed. The zonalmean GWperturbed shears have peaks (13–17 ms^{−1} km^{−1}) at around the mesopause region, i.e., at z = 90–100 km at most latitudes and at z = 80–90 km around the cold summer mesopause. This latitude–height pattern is robust over the 18 years and agrees with model simulations. The magnitudes of the GWperturbed shears exhibit yeartoyear variations and agree with the lidar and sounding rocket observations in a climatological sense but are 60 %–70 % of the model results in the zonalmean sense. The GWperturbed shears are hemispherically asymmetric and have strong annual oscillation (AO) at around 80 km (above 92 km) at the northern (southern) middle and high latitudes. At middle to high latitudes, the peaks of AO shift from winter to summer and then to winter again with increasing height. However, these GWperturbed shears may be overestimated because the GW propagation direction cannot be resolved by the method and may be underestimated due to the observational filter, sampling distance, and cutoff criterion of the vertical wavelength of GWs.
In the mesosphere and lower thermosphere (MLT), large horizontal winds and their vertical shears (or more precisely vertical wind shears, shears for short) have been revealed from more than 500 wind profiles observed by sounding rockets (Larsen, 2002; Larsen and Fesen, 2009) and from groundbased lidar and radar observations (Larsen and Fesen, 2009; Oppenheim et al., 2009, 2014; Yue el al., 2010). The large horizontal winds (≥ 100–200 ms^{−1}) and shears (≥ 40 ms^{−1} km^{−1}) occur in the height range of z∼ 95–115 km at lower and middle latitudes. Using the lidar and falling sphere observations at high latitudes (the Andøya Rocket Range and the ALOMAR – Arctic Lidar Observatory for Middle Atmosphere Research – observatory at 69.3^{∘} N) in July 2002, Fritts et al. (2004) showed large winds and shears at z∼ 85–95 km and ascribed them to strong gravity wave (GW) activity. Most prominent sources of GWs are convection, orography, and jets and fronts in the troposphere as well as spontaneous adjustment and secondary wave generation in the stratosphere (Fritts and Alexander, 2003). Amplitudes and shears increase when the GWs propagate into the lower thermosphere. It is now accepted that the large winds and shears are a common phenomenon in the MLT region (Liu, 2007, 2017). The large winds and shears are associated with tides and GWs and play an important role in forming the middlelatitude sporadic E layers and driving the equatorial electrojet (Mathews, 1998; Hysell et al., 2002; Arras et al., 2009; Haldoupis, 2012; Shinagawa et al., 2017; Arras and Wickert, 2018; Jacobi et al., 2019; Yu et al., 2019), in controlling atmospheric stabilities and the propagation of GWs, and in transporting and mixing tracers locally and/or globally (Fritts et al., 2004; Liu 2007, 2017; Yue et al., 2013; Stevens et al., 2014).
Based on the definition of the Richardson number, $Ri={N}^{\mathrm{2}}/{S}^{\mathrm{2}}$, the atmosphere is dynamically stable when $Ri>\mathrm{1}/\mathrm{4}$ and dynamically unstable when$\phantom{\rule{0.125em}{0ex}}Ri\le \mathrm{1}/\mathrm{4}$. Here ${N}^{\mathrm{2}}=g/\overline{T}\left(d\overline{T}/\mathrm{d}z+g/{c}_{p}\right)$ is the static stability; ${S}^{\mathrm{2}}={\left(\partial u/\mathrm{d}z\right)}^{\mathrm{2}}+{\left(\partial v/\mathrm{d}z\right)}^{\mathrm{2}}$ is the wind shear; $\overline{T}$ is the background temperature; g and c_{p} are the gravitational acceleration and specific heat for dry air at constant pressure, respectively; and u and v are the zonal and meridional winds, respectively. The threshold of dynamic instability is $Ri=\mathrm{1}/\mathrm{4}$, which means that the maximum wind shear allowed by the background static stability should be S=2N. According to S=2N and the simulated temperature, Liu (2007) showed that the maximum wind shears agree well with the observed large wind shears. However, the globalscale models (e.g., TIMEGCM: thermosphere–ionosphere–mesosphere–electrodynamics general circulation model) cannot reproduce the observed large winds and wind shears but can increase the amplitudes of winds and wind shears with finer spatial resolutions (Larsen and Fesen, 2009). This dependence on model resolution implies that the tides and smallscale waves (e.g., GWs) are the cause of large winds and shears seen in the observations. Using a regionscale numerical model, Liu et al. (2014a) showed that the nonlinear interactions between GWs and tides can produce large winds and shears in the MLT region. The tidal phases modulate the peak height of large winds and shears. Using the spectral element version of the National Center for Atmospheric Research (NCAR) Whole Atmosphere Community Climate Model (WACCM) with a horizontal resolution of ∼ 25 km and vertical resolution of 0.1 scale height, Liu (2017) reproduced the large wind shears, which are in good agreement with observations. Through scale separation, Liu (2017) proposed that smallscale waves (with zonal wavenumber > 6), likely GWs, play a dominant role in producing large shears. The highresolution WACCM can resolve GWs with scales longer than ∼ 100 km. Tidal waves make a secondary contribution to the magnitudes of shears but can modulate the shears produced by GWs.
Observations on large winds and shears were made by a sounding rocket and groundbased lidar and radar for a very limited number of locations and hence cannot provide a global morphology. Highresolution GCM (e.g., WACCM) simulations can provide a global picture of large winds and shears but need be validated through global observations (e.g., satellites). Moreover, it is still challenging to study the intraannual and/or the interannual variations of large winds and shears through highresolution GCM simulations due to the computational cost (Liu, 2017, 2019).
Satellite observations provide a good opportunity to study the climatology of global winds. The Wind Imaging Interferometer (WINDII) and the High Resolution Doppler Imager (HRDI) instruments on board the Upper Atmospheric Research Satellite (UARS) provide global observations of winds in the MLT region (McLandress et al, 1996; Zhang et al., 2007; Shepherd et al., 2012). Combining the winds observed by HRDI and data assimilation system, Swinbank and Ortland (2003) developed a climatology that describes the monthly zonalmean zonal winds from the surface to the upper mesosphere. The TIMED Doppler Interferometer (TIDI) instrument on board the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite measures global winds in the MLT region (Killeen et al., 2006; Niciejewski et al., 2006). These observations advanced our knowledge on the global winds in the MLT region. However, the global characteristics of shears are poorly known due to either the limited altitude coverage or the altitude resolution or large noise of the satellite observations. Since large shears play an important role in the atmospheric dynamics in the MLT region and ionospheric E region and since they are likely caused by GWs (Fritts et al., 2004; Liu, 2007, 2017; Yue et al., 2010), it should be possible to derive wind shears by combining the GW theory and GWs derived from other observed physical quantities (e.g., temperature).
The focus of this paper is to propose a method of deriving shears from GW analyses of temperature observations. Temperature profiles measured by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument (Russell et al., 1999) on board the TIMED satellite from 2002 to 2019 are used for this study providing a 18year period. Temperature profiles have been remarkably stable (Mlynczak et al., 2020). This indicates the stability of the SABER instrument is very high over the recent 18 years of measurements. These profiles cover an altitude range of ∼ 15–110 km and latitude range of 53^{∘} S–83^{∘} N or 83^{∘} S–53^{∘} N. The temperature accuracy is 1–3 K from 30 to 80 km and 5–10 K from 90 to 100 km as reported at http://saber.gatsinc.com/ (last access: March 2020) based on Remsberg et al. (2008).
The remainder of this paper is organized as follows. Section 2 presents the method of deriving wind shears induced by GWs and the uncertainties. Section 3 presents the comparisons of the GWperturbed shears with the model and observational results. Then the latitudinal and intraannual variations of the GWperturbed shears are presented in Sect. 4. The limitations of the method and their possible influences on the GWperturbed shears are discussed in Sect. 5. A short summary is given in Sect. 6.
2.1 Theory of deriving GWperturbed shears
The basic idea of deriving GWperturbed shears is the linear GW theory, which includes the dispersion and polarization relations of a monochromatic GW. For conservative propagation and without refraction, a linear monochromatic GW can be described as (Fritts and Alexander, 2003)
where $i=\sqrt{\mathrm{1}}$ is the imaginary unit. The subscript j denotes a monochromatic GW. ${u}_{j}^{\prime}$ and ${v}_{j}^{\prime}$ are the horizontal wind perturbations parallel and perpendicular to the wave vector of the GW, respectively. ${T}_{j}^{\prime}$ and $\overline{T}$ are the perturbation and background temperatures, respectively. ${\stackrel{\mathrm{\u0303}}{u}}_{j}$, ${\stackrel{\mathrm{\u0303}}{v}}_{j}$, and ${\stackrel{\mathrm{\u0303}}{T}}_{j}$ are the amplitudes of ${u}_{j}^{\prime}$, ${v}_{j}^{\prime}$, and ${T}_{j}^{\prime}/\overline{T}$, respectively. ${\mathit{\phi}}_{j}={k}_{j}x+{l}_{j}y+{m}_{j}z{\mathit{\omega}}_{j}t$ is the phase of GW. k_{j}, l_{j}, and m_{j} are the wavenumbers in the horizontal (x,y) and vertical (z) directions, respectively. ω_{j} and t are the groundbased frequency and time.
Based on the polarization of the monochromatic linear GWs with lower and medium frequencies (Fritts and Rastogi, 1985; Fritts and Alexander, 2003), the relations between ${\stackrel{\mathrm{\u0303}}{u}}_{j}$, ${\stackrel{\mathrm{\u0303}}{v}}_{j}$, and ${\stackrel{\mathrm{\u0303}}{T}}_{j}$ can be derived as (Eckermann et al., 1995; Liou et al., 2003; Gubenko et al., 2008)
where ${\widehat{\mathit{\omega}}}_{j}$ and f=2Ωsin ϕ ($\mathrm{\Omega}=\mathrm{7.292}\times {\mathrm{10}}^{\mathrm{5}}$ s^{−1}, ϕ is latitude) are the intrinsic and inertial frequencies, respectively. The wind shear of each monochromatic GW can be written as
Here we only consider the GWs propagating or projecting in the alongtrack direction. In reality, usually a spectrum of GWs is observed, formed by superposition of several monochromatic GWs. For each monochromatic GW, Eqs. (2)–(5) are valid. Here, we use u^{′} to represent observed GWs, which contains multiple monochromatic GWs (e.g., many ${u}_{j}^{\prime}$). These monochromatic GWs may not propagate in the same direction. However, from two adjacent GW profiles, one cannot get the actual horizontal wavenumber (k_{j,hr}) of each monochromatic but can get the projection of the actual wave vector in the alongtrack direction (k_{j,ha}). The inconsistency between k_{j,hr} and k_{j,ha} introduces uncertainties of GWperturbed shears, which are discussed in Sect. 2.3. Thus, u^{′} can be expressed as the vector sum of the projections of all the actual monochromatic GWs in the alongtrack direction and formulated as
In the same way, the component v^{′} of an observed GWs can be expressed as
For the GWperturbed shears of u^{′} and v^{′}, we have
Finally, the magnitude of the GWperturbed shear can be written as
Due to the uncertainties in determining the GWs' propagation direction, only the magnitudes of GWperturbed shears are analyzed in this work.
2.2 GWperturbed shear derived from synthesized GW profiles and validations
To demonstrate the applicability of the theory and the procedure of retrieval the GWperturbed shears, we construct two synthesized temperature perturbation profiles (shown in Fig. 1a). Each profile is the sum of three monochromatic GWs with three different vertical wavelengths and heightdependent amplitudes,
The background temperature is $\overline{T}\left(z\right)=\mathrm{240}$ K. The black line in Fig. 1a shows the ${T}_{\mathrm{w}\mathrm{1}}^{\prime}\left(z\right)$ and is noted as ${T}_{\mathrm{w}\mathrm{1}}^{\prime}$. The profile ${T}_{\mathrm{w}\mathrm{2}}^{\prime}$ (red line in Fig. 1a) has the same amplitudes and vertical wavenumbers as ${T}_{\mathrm{w}\mathrm{1}}^{\prime}$ but the setting ${\mathit{\phi}}_{\mathrm{1}}={\mathit{\phi}}_{\mathrm{2}}={\mathit{\phi}}_{\mathrm{3}}=\mathit{\pi}/\mathrm{2}$. ${T}_{\mathrm{w}\mathrm{1}}^{\prime}\left(z\right)$ and ${T}_{\mathrm{w}\mathrm{2}}^{\prime}\left(z\right)$ are used to represent two adjacent SABER measurements. It should be noted that we set the three monochromatic GWs in ${T}_{\mathrm{w}\mathrm{2}}^{\prime}\left(z\right)$ as having the same phase shift of π∕2 only for the convenience of theoretical representation. In real atmosphere, GWs with different vertical wavelengths may have different horizontal wavelengths and thus different phase shifts. These phase shifts can be calculated by comparing the phases of the two monochromatic GWs with the same vertical wavelength embedded in the two adjacent GW profiles. The phase of each monochromatic GW can be derived through discrete wavelet transformation (in short DWT – discrete wavelet transform; Torrence and Compo, 1998). Then the phase shifts can be used to estimate the horizontal wavenumber (e.g., Preusse et al., 2002; Ern et al., 2004; Alexander et al., 2008, 2018; Wang and Alexander, 2010). For example, if we assume the horizontal distance (Δr) between the two profiles is 300 km, the phase shift of $\mathrm{\Delta}\mathit{\phi}=\mathit{\pi}/\mathrm{2}$ indicates the horizontal wavenumber ${k}_{\mathrm{h}}=\mathrm{\Delta}\mathit{\phi}/\mathrm{\Delta}r=\left(\mathit{\pi}/\mathrm{2}\right)\left(\mathrm{300}\phantom{\rule{0.125em}{0ex}}\mathrm{km}\right)=\mathrm{2}\mathit{\pi}/\left(\mathrm{1200}\phantom{\rule{0.125em}{0ex}}\mathrm{km}\right)$. For the lower and mediumfrequency GWs, the dispersion relation can be simplified as (Fritts and Alexander, 2003)
After we get the horizontal wavenumber k_{j,h} from a GW profile pair, then the intrinsic frequency for ${T}_{j}^{\prime}\left(z\right)$ can be calculated by Eq. (12). It should be noted that the k_{j,h} is the projection of the actual GW's horizontal wavenumber in the alongtrack direction and is underestimated. This underestimates the intrinsic frequency and thus overestimates the wind and shears based on Eqs. (6)–(9) and will be discussed detail in Sect. 2.3. Using Eqs. (2)–(5) and the prescribed amplitude and vertical wavenumber, which will be determined from the satellite observation and described below, for the individual monochromatic GW, we can get ${u}_{j}^{\prime}$ and ${v}_{j}^{\prime}$ and their shears. Then, according to Eqs. (6)–(9), the vector sum of the three monochromatic GWperturbed wind profiles can be obtained and are shown as black lines in Fig. 1b for zonal wind and Fig. 1c for meridional wind. Here we assume the latitude is at 30^{∘} N, which is a typical latitude at middle latitudes. The corresponding shears are also shown as black lines in Fig. 1d–e. Since the amplitude, vertical, and horizontal wavenumbers of GW profile pairs are prescribed and are not derived through DWT, which is a key step of the spectral decomposition method described below, the retrieved winds and shears are referred to as theoretical results. These theoretical results can be used to measure the results obtained from the spectral decomposition method proposed below. Here we assume that ${u}_{j}^{\prime}$ and ${v}_{j}^{\prime}$ are the winds along and cross the orbit track directions, respectively, and may not coincide with the eastward and northward directions. The alongtrack direction is points from the location of ${T}_{\mathrm{w}\mathrm{1}}^{\prime}\left(z\right)$ to that of ${T}_{\mathrm{w}\mathrm{2}}^{\prime}\left(z\right)$. The crosstrack direction is 90^{∘} counterclockwise from the alongtrack direction. However, GWs may not propagate only in the alongtrack direction. The assumption will introduce uncertainties and will be discussed below.
Now we describe the method of retrieving winds and wind shears induced by GW profile pairs, whose amplitudes, vertical, and horizontal wavenumbers are not prescribed but should be evaluated. We name this method the “spectral decomposition method”, since the principle ideas are the following: (1) decomposing an observed GW profile into multiple monochromatic waves, (2) applying linear GW theory on each monochromatic wave to get the monochromatic wind and shear of each wave component, and (3) finding the vector sum of monochromatic winds and shears to get the wind perturbations and shears induced by the observed GW. The detailed application of this method is described as the following three steps.
The first step is to evaluate the amplitude and vertical wavenumbers of each GW profile by the method of DWT such that we can get the heightdependent amplitudes and vertical wavelengths as well as the phase shifts. For the GW profiles of ${T}_{\mathrm{w}\mathrm{1}}^{\prime}\left(z\right)$ and ${T}_{\mathrm{w}\mathrm{2}}^{\prime}\left(z\right)$, then at each height and vertical wavelength (${\mathit{\lambda}}_{z}=\mathrm{2}\mathit{\pi}/m)$, their DWT are ${\widehat{T}}_{\mathrm{w}\mathrm{1}}\left(z,{\mathit{\lambda}}_{z}\right)$ and ${\widehat{T}}_{\mathrm{w}\mathrm{2}}\left(z,\phantom{\rule{0.125em}{0ex}}{\mathit{\lambda}}_{z}\right)$. Then, their cospectrum C_{1,2} is computed as
Here, ${\stackrel{\mathrm{\u0303}}{T}}_{\mathrm{w}\mathrm{1}}=\left{\widehat{T}}_{\mathrm{w}\mathrm{1}}\right$ and ${\stackrel{\mathrm{\u0303}}{T}}_{\mathrm{w}\mathrm{2}}=\left{\widehat{T}}_{\mathrm{w}\mathrm{2}}\right$ are the amplitudes of ${T}_{\mathrm{w}\mathrm{1}}^{\prime}\left(z\right)$ and ${T}_{\mathrm{w}\mathrm{2}}^{\prime}\left(z\right)$; ${\widehat{T}}_{\mathrm{w}\mathrm{2}}^{\ast}$ is the complex conjugation of ${\widehat{T}}_{\mathrm{w}\mathrm{2}}$. The phase shift is calculated by $\mathrm{\Delta}{\mathit{\phi}}_{\mathrm{1},\mathrm{2}}={\mathrm{tan}}^{\mathrm{1}}\left[\mathrm{Im}\left({C}_{\mathrm{1},\mathrm{2}}\right)/Re\left({C}_{\mathrm{1},\mathrm{2}}\right)\right]$. When performing DWT, we restrict the vertical wavelength ranging from ∼ 5 km to ∼ 30 km for a vertical extent of 90 km, which is the height coverage (18–108 km) of the SABER temperature profiles.
The second step is to evaluate the horizontal wavenumber through the phase shift between two adjacent GW profiles. According to the distance between the two adjacent profiles Δr=300 km and phase shift, we can get ${k}_{\mathrm{h}}=\mathrm{\Delta}{\mathit{\phi}}_{\mathrm{1},\mathrm{2}}/\mathrm{\Delta}r=\left(\mathrm{\Delta}{\mathit{\phi}}_{\mathrm{1},\mathrm{2}}\right)/\left(\mathrm{300}\phantom{\rule{0.125em}{0ex}}\mathrm{km}\right)$ (Ern et al., 2004; Alexander, 2008). Then the intrinsic frequency $\widehat{\mathit{\omega}}$ for the component of m can be calculated by Eq. (12). Here we note that the horizontal wavelengths of GWs in real atmosphere may be shorter than 2Δr; only GWs with horizontal wavelengths longer than 2Δr are considered due to the sampling distances and the limbscanning mode of the SABER instrument (Preusse et al., 2002; Ern et al., 2004).
The third step is to calculate the GWperturbed wind (shown as red dashed lines in Fig. 1b–c) and shear (red dashed lines in Fig. 1d–e) by Eqs. (6)–(9). Then we can get the amplitudes of wind shears, which are the modules of $\partial {u}^{\prime}/\partial z$ and $\partial {v}^{\prime}/\partial z$, respectively. Finally, the GWperturbed shear (S) can be calculated by Eq. (10).
A brief summary of the results from the spectral decomposition method (red) and theory (black) is shown in Fig. 1. From Fig. 1, we can see that the GWperturbed winds and shears derived from spectral decomposition method (red dashed lines) agree well with the theoretical results (black solid lines) below 100 km. The bad consistencies occur at around the upper boundary. Thus, we will focus on the results in z = 30–100 km in the following analysis.
2.3 GWperturbed shear derived from SABER GW profiles and uncertainties
A key step to derive the GWperturbed shears is the extraction of GWs from the SABER temperature profile. The extraction methods of GWs from satellite data have been developed by Fetzer and Gille (1994) and improved greatly since (e.g., Preusse et al., 2002; Ern et al., 2004, 2011, 2018; Chen et al., 2019; Alexander et al., 2008, 2018; Wang and Alexander, 2010; Alexander, 2015). We have developed a similar method in our previous studies (Liu et al., 2014b, 2017, 2019), which is summarized here. First, the daily SABER temperature profiles T(z) in a latitude band of 5^{∘} are selected. Second, at each altitude, these selected data are fitted by harmonics with zonal wavenumbers ranging from 0 to 6, which are mainly planetary waves and nonmigrating tides and are removed from T(z) to get the residual temperature T_{r}(z). The component of wavenumber 0 is considered as the zonalmean temperature $\overline{T}\left(z\right)$. Due to the slowing precessing of the SABER measurement (a full cycle is about 120 d), migrating tides will appear as stationary zonal wave patterns if data from ascending and descending nodes are taken separately. The above two steps are applied on both the ascending and descending nodes, respectively, such that it minimizes the influences of migrating tides on the residual temperature T_{r}(z) and thus GW profiles (Preusse et al., 2002; Ern et al., 2004, 2018). Third, DWT (Eq. 13) is applied on each residual profile to get its wavelet transform. When applying DWT, we restrict the vertical wavelengths in the range from ∼ 5 to ∼ 30 km. From these monochromatic waves, we reconstruct a GW profile ${T}_{\mathrm{w}}^{\prime}\left(z\right)$.
A GW profile pair is defined as the two adjacent SABER GW profiles, whose alongtrack distance is less than 400 km. The 400 km criterion is fulfilled only for the shortdistance pairs for the SABER measurement (Fig. 1 of Ern et al., 2011). Then from a given GW profile pair, we can get pairs of GWperturbed winds and shears by the spectral decomposition method. Figure 2 shows the procedure of deriving winds and shears from a GW profile pair at 42.55^{∘} N, 56.48^{∘} E and 44.97^{∘} N, 56.34^{∘} E on 1 January 2018. Figure 2b and c show the GWperturbed winds in the along and crosstrack directions, respectively, through Eqs. (6) and (7). Figure 2d and e show the GWperturbed shears in the along and crosstrack directions, respectively, through Eqs. (8) and (9). Figure 2f shows the GWperturbed shears calculated by Eq. (10). From Fig. 2, we can see that the peak heights of the temperature, wind, and shear are not at the same height due to their phase differences shown by the polarization relations. The winds and shears in the alongtrack direction are larger than those in the crosstrack direction due to the factor of $f/\widehat{\mathit{\omega}}<\mathrm{1}$ in Eq. (3).
The horizontal wavenumber k_{h}, which is derived from a GW profile pair, is in the alongtrack direction. It is smaller than the real wavenumbers, since GWs may not propagate exactly in the alongtrack direction. This will induce uncertainties in deriving the GWperturbed shears, since the GW propagation direction cannot be determined from a GW profile pair (Ern et al., 2004). The uncertainties of the GWperturbed shears are estimated in below.
The relation of the horizontal wavenumbers in the alongtrack direction (k_{ha}) and in the GW propagation direction (k_{hr}) can be expressed as k_{ha}=k_{hr}cos α. The subscripts “a” and “r” denote the alongtrack and real physical quantities, respectively. The angle (α) between the alongtrack direction and GW propagation direction can vary from 0 to 360^{∘}. Here, the angle α is restricted in the range of 0–90^{∘}. This restriction maps the angle of four quadrants into the first quadrant, since only the magnitude of k_{ha} and k_{hr} are considered here. If α≠0, this will induce uncertainties of intrinsic frequencies ($\widehat{\mathit{\omega}}$) and thus GWperturbed shear (S). Here the uncertainty of S (noted by S_{e}) is defined as the ratio between the alongtrack S (noted by S_{a}) and real S (noted by S_{r}). According to Eqs. (10) and (4)–(5), we get
Equation (14) shows that if α=0 or f=0, then S_{e}=1 and the GWperturbed shear is accurate. If α≠0 and $f=\mathrm{2}\mathrm{\Omega}\mathrm{sin}\mathit{\varphi}\ne \mathrm{0}$, then S_{e} > 1 and the total wind shear is overestimated.
Figure 3 shows the dependencies of S_{e} on horizontal (λ_{h}) and vertical wavelengths (λ_{z}), the angle α, and latitude ϕ. From Fig. 3a, we can see that S_{e} increases with the increasing λ_{h} and the decreasing λ_{z}. Figure 3b shows that S_{e} increases with the increasing λ_{h} and the increasing ϕ. Figure 3c shows that S_{e} increases with the increasing λ_{h} and the increasing α. Regarding the relative importance of λ_{z}, λ_{h}, α, and ϕ in changing S_{e}, the angle α is the most important. If we assume that GW propagates in an arbitrary direction, S_{e} is less than 1.2 (indicated by a red contour line) for a large fraction of GWs at $\mathit{\varphi}=\mathrm{30}{}^{\circ}$. The fractions of GWs, which are overestimated by 20 %, increase with the increasing latitudes. This indicates that the method proposed here can be used to estimate GWperturbed shears even though the wave propagation direction cannot be determined from a GW profile pair. When we analyzed the derived total wind shears, the overestimation at high latitudes should be considered.
The large winds and shears are in the MLT region, where both tides and GWs reach large amplitudes and may interact nonlinearly (Fritts and Vincent, 1987; Li et al., 2009; Liu et al., 2014a) and break (Fritts et al., 2004; Liu and Vadas, 2013; Vadas and Liu, 2013). The largescale winds may rotate with height and act as critical levels, which filter out a broadspectral range of GWs. These filtered GWs may break and deposit their momenta in the background atmosphere, which create body force to general secondary GWs. Moreover, the vertical wavelengths of these filtered GWs change rapidly with height and in the nonlinear regime. These nonlinear GWs may also produce large winds and shears around the mesopause regions. However, these nonlinear GWs cannot be described by linear GW theory proposed here. Thus, the GWperturbed shears underestimate the actual shears in the MLT region due to the unrepresented nonlinear GWs.
Due to the uncertainties, the GWperturbed shears will be further examined by comparing with model and observational results. According to the wind shears during 1–10 July simulated by WACCM (Liu, 2017) and observations, we take the GWperturbed shears during May–August of 2018 as an example for the purpose of comparison. The longer date coverage is chosen to include data from both northviewing and southviewing yaw periods and hence to cover a wider latitude range. Figure 5 shows the zonal mean and standard deviations (SDs) of GWperturbed shears and the top 10 % largest shears during the three periods of 2018. The latitude band used for calculating zonal means has a width of 5^{∘} with an overlap of 2.5^{∘}. The latitude coverage is from 52.5^{∘} S to 82.5^{∘} N or from 82.5^{∘} S to 52.5^{∘} N due to the yaw cycle of the SABER measurement. The three periods have centers in July and extend 1 or 2 months such that the results can illustrate the main features during summer.
Three main features can be found in the zonalmean S shown in Fig. 4. Firstly, the maxima of S are ∼ 13–17 ms^{−1} km^{−1} with SDs of ∼ 9–13 ms^{−1} km^{−1} at z∼ 90–100 km (around the mesopause) at all latitudes of the three different time intervals. Especially, the GWperturbed shears reach their maxima at around 70^{∘} S. The maxima S are ∼ 10–13 ms^{−1} km^{−1} with an SD of ∼ 8–10 ms^{−1} km^{−1} at z ∼ 80–90 km at latitudes higher than 40^{∘} N, where it is near the summer mesopause (cf. contour lines of the zonalmean temperature). Secondly, the zonalmean shears have similar latitude–height patterns during the three different time intervals, although they are averaged over different time interval lengths (e.g., 31 d during 1 July–31 July in 2018 and 62 d during 29 June–30 August in 2018). The same is true during the intervals of 6 May–27 June and 29 June–30 August in 2018, respectively, at latitudes of 52.5^{∘} S–52.5^{∘} N. The GWperturbed shears during the continuous two yaw cycles (6 May–27 June and 29 June–30 August) exhibit a smooth extension from 52.5^{∘} S–52.5^{∘} N to higher latitudes, respectively. This shows that GWperturbed shears are a common phenomenon around the mesopause region (Larson, 2002; Fritts et al., 2004; Larsen and Fesen, 2009; Yue el al., 2010; Liu, 2007, 2017). Thirdly, the top 10 % largest S reach maxima of ∼ 30 ms^{−1} km^{−1} around the mesopause. In general, at each height, the GWperturbed shears reach their minima at lower latitudes and reach their maxima at high latitudes.
3.1 Comparisons with model results
The latitude–height patterns of S derived here agree well those simulated by WACCM during 1–10 July (Liu, 2017). Specifically, the GWperturbed shears derived here have peaks at ∼ z= 80–90 km at latitudes higher than 50^{∘} N during 6 May–27 June in 2018. Moreover, there is another peak at ∼ z= 90–100 km of the southern high latitudes during 29 June–30 August in 2018. These peaks agree with the WACCM simulation results that the large shears have peaks at around the mesopause (Xu et al., 2007; Fig. 2a of Liu, 2017). However, the shear peaks derived here are at a slightly lower height than the WACCM results (Liu, 2017). The GWperturbed shears might be influenced by temperature uncertainties of SABER measurements, which are much larger at around the cold summer mesopause (Remsberg et al., 2008). We removed waves with vertical wavelengths shorter than 5 km to minimize the noise introduced by uncertainties of SABER measurements (Ern et al., 2011, 2018).
Compared to the WACCM simulation results, the differences pertain to the magnitudes of the zonalmean S and specifically the following:

The maxima of the zonalmean S. Figure 2a of Liu (2017) showed that the maxima of the zonalmean S are 20–40 ms^{−1} km^{−1} near the mesopause at latitudes higher than 50^{∘} N and 50^{∘} S and are ∼ 20–25 ms^{−1} km^{−1} at latitudes lower than 50^{∘} N and 50^{∘} S. This indicates that the GWperturbed zonalmean S profiles derived here are about 70 % of those of the WACCM simulation results.

The top 10 % largest S. The simulated top 10 % largest S profiles (e.g., the minima of the total 10 % largest S) are larger than 50 ms^{−1} km^{−1} at high latitudes and 35 ms^{−1} km^{−1} over the Equator (Fig. 5 of Liu, 2017). Thus, the top 10 % largest S derived here (with maxima of ∼ 30 ms^{−1} km^{−1}) is about 60 % of the simulated results.

The SDs. The simulated SDs reach their maxima of ∼ 1.4 times of the shears (Fig. 2b of Liu, 2017). The ratio of 1.4 indicates that the maxima of the simulated SD is about 40 ms^{−1} km^{−1}, which is larger than those derived from observations (∼ 11–13 ms^{−1} km^{−1}, see the middle row of Fig. 4) at latitudes lower than 50^{∘} N and 50^{∘} S.

The structures of S. The structures in Fig. 2b of Liu (2017) have much finer scales than those in Fig. 4. The different spatial scale coverages, which will be discussed in Sect. 5, might be responsible for the smaller values of wind shears derived from observations.
In general, the GWperturbed shears derived here can reproduce the latitude–height pattern and 60 %–70 % of the simulated shear magnitudes in the zonalmean sense.
3.2 Comparisons with observational results
To compare the GWperturbed shears derived here and those observed by groundbased lidar and sounding rockets (Larsen, 2002; Yue et al., 2010), we show in Fig. 5 the profiles of S and their zonal means as well as the top 10 % and 1 % largest S during January and July at around 40^{∘} N and the Equator. January and July are representative months for winter and summer, respectively. The latitudes of the Equator and 40^{∘} N may be representative for low and middle latitudes, where sounding rocket measurements were performed (Fig. 1 of Larsen, 2002; Larsen and Fesen, 2009). The 40^{∘} N is near the latitude of 41^{∘} N, where the Colorado State University (CSU) Na lidar observations were performed (Yue et al., 2010).
Figure 5 shows that the magnitudes of the GWperturbed S profiles increase with height and can be larger than 40 ms^{−1} km^{−1} above 90 km (80 km) during January (July). The height variations and magnitudes of GWperturbed S profiles derived here compare well with the over 400 chemical tracer measurements by sounding rocket (Fig. 10 of Larsen, 2002) below 100 km. Moreover, the height variations and magnitudes of GWperturbed S profiles derived here compare well with the CSU lidar observations below 100 km. Specifically, the magnitudes of S profiles observed by CSU lidar increase with height and have maxima of ∼ 40 ms^{−1} km^{−1} (Fig. 1c–d of Yue et al., 2010). This agrees well with the S profiles shown in Fig. 5. After averaging the shears observed by the CSU lidar during summer and winter months, the zonalmean shears are ∼ 12–17 ms^{−1}km^{−1} with SDs of ∼ 10 ms^{−1} km^{−1} (Fig. 11c of Yue et al., 2010). They are slightly larger than the magnitudes of 10–13 ms^{−1} km^{−1} derived here (the right column of Fig. 5). The zonalmean shears (the right column and upper row of Fig. 5) derived here have similar magnitudes during January and July above 90 km. However, the mean shears observed by the CSU lidar have similar magnitudes in winter and summer at ∼ z=80–105 km. We note that the magnitudes of S derived here are larger in July than in January at the Equator. This might be related to the intraannual variations of S and will be studied in Sect. 4.
A short summary of the above comparisons is below. The GWperturbed shears derived from SABER observations agree with the previous observations and model results in general. This provides a global climatology of large shears around the mesopause region partially based on observations. The magnitudes of S derived from SABER observations are similar to those observed by lidar and sounding rocket but are about 60 %–70 % of the highresolution WACCM results in the zonalmean sense. The difference probably comes from (1) the coarse horizontal samplings of ∼ 250–350 km of satellite observations; (2) the use of only GWs with λ_{z}≥5 km and λ_{h}≥2Δr (Δr is the distance between the two GW profiles in a pair) for deriving shears; and (3) the fact that the observational filter should underestimate the amplitudes of GWs, especially for GWs with shorter vertical wavelengths (Ern et al., 2018), and thus underestimate the GWperturbed shears. This will be further discussed in Sect. 5.
With the advantage of 18year SABER observations, the climatology of the GWperturbed shears can be explored on the aspects of their latitudinal variations and intraannual variations. The four seasons in the Northern Hemisphere (NH) are spring (MAM: March, April, and May), summer (JJA: June, July, and August), autumn (SON: September, October, and November), and winter (DJF: December, January, and February).
4.1 Latitude variations of GWperturbed shears
Figure 6 shows the latitude–height contours of the zonal mean and SDs of the GWperturbed shears and the top 10 % largest shears during four composite seasons. Each composite season is made up by the superpositions of the corresponding seasons from 2002 to 2019. The numbers of profiles used to derive GWperturbed wind shears in each season (the fourth row of Fig. 6) have peaks around 50^{∘} N and 50^{∘} S due to the changes of ascending and descending nodes. The sharp changes of the profile numbers around 50^{∘} N and 50^{∘} S might induce the discontinuity of the latitudinal variations of GWperturbed shears. By further examination of the discontinuities in Fig. 6, we find that the zonal means exhibit more obvious discontinuities at around 50^{∘} N (50^{∘} S) than those at around 50^{∘} S (50^{∘} N) during spring and summer (autumn and winter). This is because there are fewer samplings (the fourth row of Fig. 6) at around 50^{∘} N (50^{∘} S) than those at around 50^{∘} S (50^{∘} N) during spring and summer (autumn and winter).
The hemispheric asymmetry of the sampling is induced by the inconsistency of the date coverages of yaw cycle and season. Consequently, we show in Fig. 7 the zonal means and SDs of GWperturbed shears and the top 10 % largest shears during six composite yaw cycles. The composite yaw cycle is the superposition of all the yaw cycles, which have nearly identical date coverage relative to the beginning of each calendar year. For example, the first yaw cycle of each year covers 25 January–18 March in 2002, 16 January–18 March in 2003, …, and 28 December 2018–26 February 2019. These dates are mainly in January and February, with a few extending to March and December. We label “1228–0318” on the top of the first column of Fig. 7 to note the all the dates (month and day) covered by the first yaw cycle. The results in the other five composite yaw cycles are shown in the same manner. The continuous two composite yaw cycles may have overlaps, with the longest overlap time of about 20 d. Thus, the composite yaw cycle can represent the results during 2 months around the center date of each composite yaw cycle.
From Fig. 6, we can see that the zonal means of the GWperturbed shears (the first row) increase with the increasing height and latitude in general. The peaks are ∼ 10–15 ms^{−1} km^{−1} above 90 km at latitudes of 82.5^{∘} S–50^{∘} N (50^{∘} S–82.5^{∘} N) during spring and summer (autumn and winter). Moreover, the wind shears have peaks at a lower height (z ∼ 80–90 km) and at latitudes of 82.5–50^{∘} S (50–82.5^{∘} N) during autumn and winter (spring and summer). These lower height peaks during spring and autumn (highlighted by blue rectangles) are weaker than those during summer and winter (highlighted by red rectangles). Comparing with Fig. 7 (the second and third columns, 0228–0512 and 0502–0715), we find that the weak peak at z ∼ 80–90 km during spring is contributed by the wind shears during May (highlighted by a red rectangle in Fig. 7 during the yaw cycle of 0502–0715), since there is no peak at similar location during the yaw cycle of 0228–0512. The same is true during autumn, when the weak peak at z∼ 80–90 km is contributed from the wind shears during November (highlighted by a red rectangle in Fig. 7 during the yaw cycle 1031–0114), since there is no peak at a similar location during the yaw cycle of 1228–0318. The stronger peaks during summer and winter in Fig. 6 are contributed from those during the yaw cycles of 0502–0715 and 1031–0114, respectively. The stronger peaks above 90 km and at ∼ z=80–90 km (marked by red rectangles) are both at around the mesopause as referred to the zonalmean temperature (contour lines in the second rows of Figs. 6 and 7).
The SD and the top 10 % largest shears, which are shown in the second and third rows of Figs. 6 and 7, respectively, have similar patterns as that of zonalmean shears. The maxima of the SD and the top 10 % largest shears are, respectively, ∼ 12 and ∼ 30 ms^{−1} km^{−1}, which are slightly less than that shown in Fig. 4. This is because the sampling profiles in Figs. 6 and 7 (composite season or yaw cycle over 18 years) are much larger than those in Fig. 4 (only one yaw cycle in 1 year).
Since the patterns of zonal mean and SDs of the GWperturbed shears and the top 10 % largest shears are similar to each other (as shown in Figs. 4, 6, and 7), only the zonalmean shears during each summer from 2002 to 2019 are shown in Fig. 8. It can be seen that the latitude–height distributions of GWperturbed shears, including the peaks at lower heights (around the mesopause region) of high latitudes, are similar to the 18year mean results shown in Figs. 6 and 7. However, the GWperturbed shear magnitudes (shown in Fig. 8) exhibit yeartoyear variations. For example, at the Southern Hemisphere (SH) high latitudes, the wind shears above 90 km are strongest during 2008 and 2019 and weakest during 2002. At the NH high latitudes, the GWperturbed shears at ∼ z=85–95 km vary by year more greatly and have smaller values, as compared to those at around 80 km.
4.2 Intraannual oscillations of GWperturbed shears
Since the GWperturbed shears are prominent around the mesopause region, their intraannual oscillations will be studied at ∼ z=60–100 km. Figure 9 shows the monthly zonalmean GWperturbed shears at four latitudes bands of the NH and SH from 2002 to 2019. A general feature of time–height variations GWperturbed shears are the annual (AO) and semiannual oscillations (SAO). To quantify the exact amplitudes and phases of AO and SAO, harmonic fitting is applied on the GWperturbed shears. The fitting function has periods of both AO and SAO. Figure 10 shows the amplitudes and phases of both AO and SAO at four latitude bands of the Northern Hemisphere and Southern Hemisphere (SH).
At 50^{∘} N and 50^{∘} S (the first row of Fig. 9 and the first column of Fig. 10), the GWperturbed shears exhibit different height dependencies of AO and SAO. At 50^{∘} N, both AO and SAO reach their maxima at 80 km, while SAO has another peak at 97 km. At z=75–92 km, the AO is dominant and peaks in June. Below 75 km and above 92 km, the AO and SAO are almost equal partitioned. At 50^{∘} S, both AO and SAO reach their maxima at ∼ 81 km, while AO has another peak at 98 km. Above 92 and below 68 km, the AO is dominant and has phase in June. At z=80–90 km, the AO and SAO are almost equally partitioned and have peaks in December and June, respectively. The amplitudes of SAO at 50^{∘} N and 50^{∘} S have similar amplitudes and height variations. However, the amplitude of AO at 50^{∘} N is smaller than that at 50^{∘} S. This makes the GWperturbed shears hemispherically asymmetric.
It should be noted that the phase of AO at 50^{∘} N shifts from December at ∼ 65 km to June at ∼ 75 km and then shifts from June at ∼ 88 km to December at ∼ 100 km, whereas the phase of AO at 50^{∘} S shifts from June at ∼ 70 km to December at ∼ 77 km and then shifts from December at ∼ 85 km to June at ∼ 65 km. In summary, in each hemisphere, the phase of AO shifts from winter below ∼ 65 km (∼ 70 km) to summer at ∼ z=75–88 km (∼ z=77–85 km) and then shifts to winter again above 95 km at 50^{∘} N (50^{∘} S), respectively.
At 35^{∘} N and 35^{∘} S (the second row of Fig. 9 and the second column of Fig. 10), the GWperturbed shears exhibit both AO and SAO. At 35^{∘} N, the amplitudes of AO and SAO vary with height in a similar pattern as those at 50^{∘} N but have smaller values. The phases of AO and SAO are also in June when their amplitudes are prominent at z=75–92 km and then shift to winter above ∼ 92 km. At 35^{∘} S, the amplitude and phase of AO vary with height in a similar pattern as those at 50^{∘} S. The amplitude of SAO is dominant below 90 km with a peak in June.
At 20^{∘} N (the third row of Fig. 9 and the third column of Fig. 10), the amplitudes of AO and SAO exhibit similar height variations as those at 50^{∘} N and 35^{∘} N but have smaller values. The AO and SAO reach their peaks at around 81 km in June and have values of ∼ 0.7 and ∼ 0.5 ms^{−1}km^{−1}, respectively. At 20^{∘} S, the SAO is in the dominant position and has peaks in June at ∼ 92 km and ∼ 67 km. At 5^{∘} N and 5^{∘} S and at z=85–98 km, the AO is in the dominant position and has peak shifting from July to March. At ∼ z=71–80, the SAO is in the dominant position and has peak shifting from March to January.
The AO and SAO of the GWperturbed shears are summarized below. The amplitudes of AO have peaks at around 80 km and decrease with the decreasing latitudes. The phases of AO shift from winter to summer and then to winter again with the increasing height. The amplitudes of SAO decrease with the decreasing latitudes. The phases of SAO are in May and June when the SAOs reach their peaks at 50, 35, and 20^{∘} N and S. At 5^{∘} N and 5^{∘} S, the SAOs shift their phase from March to January in their peak height. The AO and SAO are hemispherically asymmetric. At ∼ z=75–90 km, the AOs (SAOs) are in the dominant position at latitudes higher than 20^{∘} N (20^{∘} S). Above ∼ 90 km, the AO and SAO are almost equally partitioned at 50 and 35^{∘} N, whereas the AO is in the dominant position at 50 and 35^{∘} S.
The GWperturbed shears, which are derived from the projection of actual GWs in the alongtrack direction, are overestimated as compared to the actual GWperturbed shears. The extent of overestimation depends the actual GW propagations, which have seasonal and latitudinal preferences. Thus the seasonal and latitudinal variations of GWperturbed shears may be influenced by the preferred GWs propagation directions to some extent. Comparing with the AO and SAO in the GW square temperature amplitude (GWSTA) and absolute momentum flux (GWMF) presented by Chen et al. (2019), we find that the AO and SAO of GWperturbed shears agree with GWSTA and GWMF on the aspects of phase shifts and hemispheric asymmetry (Figs. 2 and 3 of Chen et al., 2019). However, the heights at which phase shifts occur are different. One reason is that the AO and SAO in the background temperature and static stabilities (Liu et al., 2020). The other reason is that the GWMF is inverse proportion to the vertical wavenumber m, while the GWperturbed shear is proportion to the m. The resulting height of the phase shift of GWMF is at a lower height than that of GWperturbed shears, since m increases with height below z = 90 km (Ern et al., 2018).
To determine the horizontal wavenumbers in the zonal and meridional directions, at least three profiles should be sampled at different locations of the same wave (Wang and Alexander, 2010; Alexander, 2015; Schmidt et al., 2016; Alexander et al., 2018). For the SABER measurement, there are 15 orbits in the ascending and descending nodes, respectively. The nearest distance between two orbits is about 24^{∘}, which is much longer than the horizontal wavelengths of most GWs. This limits our ability to deduce the zonal and meridional horizontal wavenumbers and leads to the uncertainties in deriving the GWperturbed shears. The horizontal wavenumber derived from a GW profile pair is in the alongtrack direction. It is in general smaller than the horizontal wavenumbers of GWs in reality. The angle (α) between the alongtrack direction and real GW propagation direction is the dominant source of the uncertainties. This will overestimate the GWperturbed shears as shown by Eq. (14) and Fig. 3. According to Eq. (14), the influences of α on S_{e} increase with the increasing latitudes due to increasing inertial frequency f. For the extreme case, S_{e}=1 for any α, since f=0 over the Equator. On the other hand, at the high latitudes, GWs propagate mostly in the zonal direction, since their sources are mainly jets or fronts and topography (Fritts and Alexander, 2003; Plougonven and Zhang, 2014). Fortunately, the SABER orbit track intersects with zonal direction at a smaller angle at high latitudes than that at lower latitudes due to the changes of ascending and descending nodes. This reduces the uncertainties of the GWperturbed shears for the zonally propagating waves. Thus, the uncertainties might be smaller than 1.2 as shown in Fig. 3c.
Even with an overestimation of 1.2, the GWperturbed shears derived here are smaller than those of highresolution model simulations by 60 %–70 % (Liu, 2017). The smaller GWperturbed shears might be induced by the following two reasons: (1) the coarse horizontal samplings of satellite observations and (2) only GWs with λ_{h}≥2Δr (Nyquist limit) and λ_{z}≥5 km used to derive shears. The reason for (1) is that the latitude–height variations of SDs in Figs. 4, 6, and 7 are smoother than those in Fig. 2b of Liu (2017). The smallscale variations in Fig. 2b of Liu (2017) might be smoothed out due to the observational filter of the SABER limbsounding pattern (Preusse et al., 2002; Ern et al., 2018).
According to the SABER sampling (Fig. 1 of Ern et al., 2011), the sampling distance of a GW profile pair Δr∼250–350 km limits the resolved GWs with λ_{h}≥500–700 km, whereas the horizontal resolution of WACCM is about 25 km (Liu, 2017); this can resolve GWs with λ_{h}≥50 km according to Nyquist limit, though waves with λ_{h}≤200 km are excessively damped by numerical diffusion in WACCM (H.L. Liu et al., 2014). The longer sampling distances miss the GWs with shorter horizontal wavelengths, which might also contribute the GWperturbed shears and thus reduces the magnitudes of GWperturbed shears. The influence of the horizontal sampling distances on wind shears can be further confirmed by the simulation results presented by Shinagawa et al. (2017), who showed that the longitude–latitude distributions of the zonal wind shears have peak values of 16–18 ms^{−1} km^{−1} at 100 km during summer and winter (Figs. 6 and 7 of their paper). According to climatology, the zonal mean of the zonal wind shears should be much smaller than 16–18 ms^{−1} km^{−1} at 100 km. This magnitude is also smaller than the simulation results presented by Liu (2017). The smaller magnitude of the wind shears might result from the different resolution used by their models. The Groundtotopside model of Atmosphere and Ionosphere for Aeronomy (GAIA) used by Shinagawa et al. (2017) has a grid size of 2.8^{∘} longitude × 2.8^{∘} latitude horizontally and a 0.2 scale height vertically, whereas the WACCM used by H.L. Liu et al. (2014) and Liu (2017) has a quasiuniform horizontal resolution of ∼25 km and a 0.1 scale height vertically.
For reason (2), the cutoff criterion of λ_{z}≥5 km is used here to get a more reliable GW profile through DWT. This cutoff criterion is related to the vertical resolution of the SABER measurement and is the same as that used by Ern et al. (2018) to remove the artificial oscillations. To test the influences of the cutoff criterions on the GWperturbed shears, we perform the same procedure as that described in Sects. 2 and 3 but relax the cutoff criterion to λ_{z}≥3 km. The GWperturbed shears derived with a cutoff criterion of λ_{z}≥3 km are shown in Fig. 11. It can be seen that the latitude–height patterns of the GWperturbed shears are the same as those shown in Fig. 4. However, the maxima of the zonalmean GWperturbed shears increase from 12–17 ms^{−1} km^{−1} for λ_{z}≥5 km to 18–24 ms^{−1} km^{−1} for λ_{z}≥3 km. This illustrates that magnitudes of the GWperturbed shears increase with the decreasing cutoff vertical wavelengths.
The uncertainties in the theory presented in Sect. 2 arise from the fact that we assume $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ except at the Equator. For the mediumfrequency GWs (${\widehat{\mathit{\omega}}}_{j}\gg f$ or $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0})$, we get
The GWperturbed shears are contributed only by $\partial {u}_{j}^{\prime}/\partial z$. Then S derived under the midfrequency assumption of $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$ is less than that derived under the assumption of $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ by a factor of
To fully explore the differences of magnitudes of S derived under the assumptions of $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ and $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$, we show the GWperturbed shears for $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$ in the same manner as those in Figs. 4–9 such that we can judge whether the assumptions need to be made. In the same manner as Fig. 4, we show in Fig. 12 the latitude–height contours of the zonal means of S derived under the assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$. The SDs and the top 10 % largest shears are not shown here, since they have similar patterns as the zonal means of S but have smaller maxima than those shown in Fig. 4. From Fig. 12 we can see that the latitude–height patterns of the S are the same as those shown in Fig. 4. However, the maxima of the zonal means of S decrease from 12–17 ms^{−1} km^{−1} for $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ to 11–15 ms^{−1} km^{−1} for $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$. The maxima of the zonal means of S at latitudes higher than 50^{∘} N are at a higher height for $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$ than that for $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$. In the same manner as Fig. 5, we show in Fig. 13 the profiles of S as well as the top 10 % and 1 % largest S derived under the assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$ during January and July at around 40^{∘} N. From Fig. 13 we can see that the height variations of these profiles are similar to those shown in Fig. 5. However, the magnitudes of these profiles are slightly smaller than those shown in Fig. 5.
In the same manner as Figs. 6–7, we show in Figs. 14–15 the zonal means of S during the four composite seasons and six yaw cycles, respectively, derived under the mediumfrequency assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$. From Figs. 14–15 we can see these results have similar patterns as those for $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ but have slightly smaller magnitudes. Moreover, the maxima of these results at high latitudes of summer hemispheres are at a higher height for $f/{\widehat{\mathit{\omega}}}_{j}\approx \mathrm{0}$ than those for $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$. The peaks circled by blue rectangles are very weak and almost disappear in Figs. 14–15 as compared to those shown in Figs. 6–7. The SDs and the top 10 % largest S are not shown here, since they have similar patterns as the zonal means of S but have smaller maxima than those shown in Figs. 6–7.
In the same manner as Figs. 8–9, we show in Figs. 16–17 the latitude–height contours and time–height contours of the zonal means of S, respectively, derived under the mediumfrequency assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$. Comparing between Figs. 8–9 and 16–17, we can see that the zonal means of S derived under the mediumfrequency assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$ have similar patterns as those for $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ but have slightly smaller magnitudes.
In summary, the GWperturbed shears derived under the assumptions of $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$ and $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$ have similar patterns on the aspects of latitude–height, time–height contours. The magnitudes of the GWperturbed shears derived under the assumption of $f/{\widehat{\mathit{\omega}}}_{j}=\mathrm{0}$ are slightly smaller than those under the assumption of $f/{\widehat{\mathit{\omega}}}_{j}\ne \mathrm{0}$.
Due to the important role the large vertical wind shears play in the dynamics and electrodynamics of the MLT, there is a need for global observation. In response to this need, a method of deriving GWperturbed shears is proposed in this work. The theoretical basis of the method is the dispersion and polarization relations of linear GWs. Data employed are SABER temperature profiles measured over the past 18 years. Based on the method and the data, the global GWperturbed shears are studied over a time span of 18 years.
The GWperturbed shears derived here agree with previous lidar and sounding rocket observations in the aspects of height structures and magnitudes in the climatological sense. Moreover, the GWperturbed shears derived here agree with the highresolution model simulation results in the aspects of latitude–height patterns but have smaller magnitudes in the zonalmean sense. The GWperturbed shears reach their maxima around the mesopause region and increase with the increasing latitudes. At most latitudes and during all seasons, the maxima of GWperturbed shears are at ∼ z=90–100 km. At high latitudes of the summer hemisphere, the maxima of GWperturbed shears are at a lower height (∼ z=80–90 km). This latitude–height pattern of GWperturbed shears is independent of the year. The magnitudes of the GWperturbed shears exhibit yeartoyear variations.
The GWperturbed shears exhibit more prominent AO and SAO at high latitudes than those at lower latitudes. The height variations of the amplitudes AO and SAO are hemispherically asymmetric. The strong AO occurs at around 80 km in the NH and above 92 km in the SH. At middle to high latitudes, the phases of AO shift from winter to summer and then to winter again with the increasing height. The amplitudes of SAO decrease with the decreasing latitudes. The phases of SAO are in May and June when the SAO reaches its peak at middle to high latitudes.
The main limitation of the method is the overestimation of the GWperturbed shears due to the unresolved GW propagation direction by the method. The other limitations, such as the observational filter, long sampling distance, and cutoff criterion of the vertical wavelength, will underestimate the GWperturbed shears. To overcome these limitations, it is necessary to develop new techniques for remote sensing temperatures or winds from space, such as limbimaging techniques allowing one to infer temperatures in three dimensions along the orbital track at high horizontal and vertical resolution.
The SABER data were downloaded from ftp://saber.gatsinc.com/Version2_0/Level2A/ (last access: March 2020). Outputs of the method are available from the corresponding author upon request.
XL and JX designed the study and wrote the paper. JY and HL contributed to the discussion of the results and the preparation of the paper. All authors discussed the results and commented on the paper at all stages.
The authors declare that they have no conflict of interest.
The SABER data were obtained from ftp://saber.gatsinc.com/Version2_0/Level2A/ (last access: March 2020). This work was supported by the National Natural Science Foundation of China (grant nos. 41831073 and 41874182). This work was also supported in part by the Natural Science Foundation of Henan, the Specialized Research Fund, and the Open Research Program of the State Key Laboratory of Space Weather. National Center for Atmospheric Research is a major facility sponsored by the National Science Foundation (cooperative agreement no. 1852977).
This research has been supported by the National Natural Science Foundation of China (grant nos. 41831073 and 41874182) and the National Science Foundation (grant no. 1852977).
This paper was edited by FranzJosef Lübken and reviewed by two anonymous referees.
Alexander, M. J.: Global and seasonal variations in threedimensional gravity wave momentum flux from satellite limbsounding temperatures, Geophys. Res. Lett., 42, 6860–6867, https://doi.org/10.1002/2015GL065234, 2015.
Alexander, M. J., Gille, J., Cavanaugh, C., Coffey, M., Craig, C., Eden, T., Francis, G., Halvorson, C., Hannigan, J., Khosravi, R., Kinnison, D., Lee, H., Massie, S., Nardi, B., Barnett, J., Hepplewhite, C., Lambert, A., and Dean, V.: Global estimates of gravity wave momentum flux from High Resolution Dynamics Limb Sounder observations, J. Geophys. Res., 113, 1–11, https://doi.org/10.1029/2007jd008807, 2008.
Alexander, P., Schmidt, T., and de la Torre, A.: A Method to Determine Gravity Wave Net Momentum Flux, Propagation Direction, and “Real” Wavelengths: A GPS Radio Occultations Soundings Case Study, Earth Space Sci., 5, 222–230, https://doi.org/10.1002/2017EA000342, 2018.
Arras, C. and Wickert, J.: Estimation of ionospheric sporadic E intensities from GPS radio occultation measurements, J. Atmos. Sol.Terr. Phys., 171, 60–63, https://doi.org/10.1016/j.jastp.2017.08.006, 2018.
Arras, C., Jacobi, C., and Wickert, J.: Semidiurnal tidal signature in sporadic E occurrence rates derived from GPS radio occultation measurements at higher midlatitudes, Ann. Geophys., 27, 2555–2563, https://doi.org/10.5194/angeo2725552009, 2009.
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/angeo374872019, 2019. 1 Introduction Back to top
Eckermann, S. D., Hirota, I., and Hocking, W. K.: Gravity wave and equatorial wave morphology of the stratosphere derived from longterm rocket soundings, Q. J. Roy. Meteor. Soc., 121, 149–186, https://doi.org/10.1002/qj.49712152108, 1995.
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.Atmos., 109, 1–17, https://doi.org/10.1029/2004JD004752, 2004.
Ern, M., Preusse, P., Gille, J. C., Hepplewhite, C. L., Mlynczak, M. G., Russell, J. M., and Riese, M.: Implications for atmospheric dynamics derived from global observations of gravity wave momentum flux in stratosphere and mesosphere, J. Geophys. Res.Atmos., 116, 1–24, https://doi.org/10.1029/2011JD015821, 2011.
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/essd108572018, 2018.
Fetzer, E. J., and Gille, J. C.: Gravity Wave Variance in LIMS Temperatures, Part I: Variability and Comparison with Background Winds, J. Atmos. Sci., 51, 2461–2483, https://doi.org/10.1175/15200469(1994)051<2461:GWVILT>2.0.CO;2, 1994.
Fritts, D. C., and Rastogi, P. K.: Convective and dynamical instabilities due to gravity wave motions in the lower and middle atmosphere: Theory and observations, Radio Sci., 20, 1247–1277, https://doi.org/10.1029/RS020i006p01247, 1985.
Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1–64, https://doi.org/10.1029/2001RG000106, 2003.
Fritts, D. C. and Vincent, R. A.: Mesospheric Momentum Flux Studies at Adelaide, Australia: Observations and a Gravity Wave–Tidal Interaction Model, J. Atmos. Sci., 44, 605–619, https://doi.org/10.1175/15200469(1987)044<0605:MMFSAA>2.0.CO;2, 1987.
Fritts, D. C., Williams, B. P., She, C. Y., Vance, J. D., Rap, M., Lübken, F.J., Müllemann, A., Schmidlin, F. J. and Goldberg, R. A.: Observations of extreme temperature and wind gradients near the summer mesopause during the MaCWAVE/MIDAS rocket campaign, Geophys. Res. Lett., 31, L24S06, https://doi.org/10.1029/2003GL019389, 2004.
Gubenko, V. N., Pavelyev, A. G., and Andreev, V. E.: Determination of the intrinsic frequency and other wave parameters from a single vertical temperature or density profile measurement, J. Geophys. Res.Atmos., 113, 1–9, https://doi.org/10.1029/2007JD008920, 2008.
Haldoupis, C.: Midlatitude sporadic E. A typical paradigm of atmosphereionosphere coupling, Space Sci. Rev., 168, 441–461, https://doi.org/10.1007/s1121401197868, 2012.
Hysell, D. L., Chau, J. L., and Fesen, C. G.: Effects of large horizontal winds on the equatorial electrojet, J. Geophys. Res.Space, 107, 1–12, https://doi.org/10.1029/2001JA000217, 2002.
Jacobi, C., Arras, C., Geißler, C., and Lilienthal, F.: Quarterdiurnal signature in sporadic E occurrence rates and comparison with neutral wind shear, Ann. Geophys., 37, 273–288, https://doi.org/10.5194/angeo372732019, 2019.
Killeen, T. L., Wu, Q., Solomon, S. C., Ortland, D. A., Skinner, W. R., Niciejewski, R. J., and Gell, D. A.: TIMED Doppler Interferometer: Overview and recent results, J. Geophys. Res.Space, 111, 1–19, https://doi.org/10.1029/2005JA011484, 2006
Larsen, M. F.: Winds and shears in the mesosphere and lower thermosphere: Results from four decades of chemical release wind measurements, J. Geophys. Res., 107, 1215, https://doi.org/10.1029/2001JA000218, 2002
Larsen, M. F. and Fesen, C. G.: Accuracy issues of the existing thermospheric wind models: can we rely on them in seeking solutions to winddriven problems?, Ann. Geophys., 27, 2277–2284, https://doi.org/10.5194/angeo2722772009, 2009.
Li, T., She, C. Y., Liu, H. L., Yue, J., Nakamura, T., Krueger, D. A., Wu, Q., Dou, X., and Wang, S.: Observation of local tidal variability and instability, along with dissipation of diurnal tidal harmonics in the mesopause region over Fort Collins, Colorado (41^{∘} N, 105^{∘} W), J. Geophys. Res.Atmos., 114, D06106, https://doi.org/10.1029/2008JD011089, 2009.
Liou, Y. A., Pavelyev, A. G., Huang, C. Y., Igarashi, K., Hocke, K., and Yan, S. K.: Analytic method for observation of the gravity waves using radio occultation data, Geophys. Res. Lett., 30, 1–5, https://doi.org/10.1029/2003GL017818, 2003.
Liu, H.L.: On the large wind shear and fast meridional transport above the mesopause, Geophys. Res. Lett., 34, L08815, https://doi.org/10.1029/2006GL028789, 2007.
Liu, H.L.: Large wind shears and their implications for diffusion in regions with enhanced static stability: The mesopause and the tropopause, J. Geophys. Res.Atmos., 122, 9579–9590, https://doi.org/10.1002/2017JD026748, 2017.
Liu, H.L.: Quantifying gravity wave forcing using scale invariance, Nat. Commun., 10, 1–12, https://doi.org/10.1038/ s4146701910527z, 2019.
Liu, H. L. and Vadas, S. L.: Largescale ionospheric disturbances due to the dissipation of convectivelygenerated gravity waves over Brazil, J. Geophys. Res.Space, 118, 2419–2427, https://doi.org/10.1002/jgra.50244, 2013.
Liu, H.L., McInerney, J. M., Santos, S., Lauritzen, P. H., Taylor, M. A, and Pedatella, N. M.: Gravity waves simulated by highresolution Whole Atmosphere Community Climate Model, Geophys. Res. Lett., 41, 9106–9112, https://doi.org/10.1002/2014GL062468, 2014.
Liu, X., Xu, J., Yue, J., Liu, H.L., and Yuan, W.: Large winds and wind shears caused by the nonlinear interactions between gravity waves and tidal backgrounds in the mesosphere and lower thermosphere, J. Geophys. Res.Space, 119, 7698–7708, https://doi.org/10.1002/2014JA020221, 2014a.
Liu, X., Yue, J., Xu, J., Wang, L., Yuan, W., Russell, J. M., and Hervig, M. E.: Gravity wave variations in the polar stratosphere and mesosphere from SOFIE/AIM temperature observations, J. Geophys. Res.Atmos., 119, 7368–7381, https://doi.org/10.1002/2013JD021439, 2014b.
Liu, X., Yue, J., Xu, J., Garcia, R. R., Russell, J. M., Mlynczak, M., Wu, D., and Nakamura, T.: Variations of global gravity waves derived from 14 years of SABER temperature observations, J. Geophys. Res.Atmos., 122, 6231–6249, https://doi.org/10.1002/2017JD0266042017, 2017.
Liu, X., Xu, J., Yue, J., Vadas, S. L., and Becker, E.: Orographic primary and secondary gravity waves in the middle atmosphere from 16year SABER observations, Geophys. Res. Lett., 46, 4512–4522, https://doi.org/10.1029/2019GL082256, 2019.
Liu, X., Xu, J. Y., and Yue, J.: Global static stability and its relation to gravity waves in the middle atmosphere, Earth Planet. Phys., 4, 1–9, https://doi.org/10.26464/epp2020047, 2020.
Mathews, J. D.: Sporadic E: current views and recent progress, J. Atmos. Sol.Terr. Phy., 60, 413–435, https://doi.org/10.1016/S13646826(97)000436, 1998.
McLandress, C., Shepherd, G. G., Solheim, B. H., Burrage, M. D., Hays, P. B., and Skinner, W. R.: Combined mesosphere/thermosphere winds using WINDII and HRDI data from the Upper Atmosphere Research Satellite, J. Geophys. Res.Atmos., 101, 10441–10453, https://doi.org/10.1029/95JD01706, 1996
Mlynczak, M. G., Daniels, T., Hunt, L. A., Yue, J., Marshall, B. T., Russell, J. M., Remsberg, E. E., Tansock, J., Esplin, R., Jensen, M., Shumway, A., Gordley, L., and Yee, J. H.: Radiometric Stability of the SABER Instrument, Earth Space Sci., 7, 1–8, https://doi.org/10.1029/2019EA001011, 2020.
Niciejewski, R., Wu, Q., Skinner, W., Gell, D., Cooper, M., Marshall, A., Killeen, T., Solomon, S., and Ortland, D.: TIMED doppler interferometer on the thermosphere ionosphere mesosphere energetics and dynamics satellite: Data product overview, J. Geophys. Res.Space, 111(A11), 1–20. https://doi.org/10.1029/2005JA011513, 2006
Oppenheim, M. M., Sugar, G., Slowey, N. O., Bass, E., Chau, J. L., and Close, S.: Remote sensing lower thermosphere wind profiles using nonspecular meteor echoes, Geophys. Res. Lett., 36, 1–5, https://doi.org/10.1029/2009GL037353, 2009.
Oppenheim, M. M., Arredondo, S., and Sugar, G.: Intense winds and shears in the equatorial lower thermosphere measured by highresolution nonspecular meteor radar, J. Geophys. Res.Space, 119, 2178–2186, https://doi.org/10.1002/2013 JA019272, 2014.
Plougonven, R. and Zhang, F: Internal gravity waves from atmospheric jets and fronts, Rev. Geophys., 52, 33–76, https://doi.org/10.1002/2012RG000419, 2014.
Preusse, P., Dörnbrack, A., Eckermann, S. D., Riese, M., Schaeler, B., Bacmeister, J. T., Broutman, D., and Grossmann, K. U.: Spacebased measurements of stratospheric mountain waves by CRISTA 1. Sensitivity, analysis method, and a case study, J. Geophys. Res.Atmos., 107, 8178, https://doi.org/10.1029/2001JD000699, 2002.
Remsberg, E. E., Marshall, B. T., GarciaComas, M., Krueger, D., Lingenfelser, G. S., MartinTorres, J., Mlynczak, M. G., Russell III, J. M., Smith, A. K., Zhao, Y., Brown, C., Gordley, L. L., LopezGonzalez, M. J., LopezPuertas, M., She, C.Y., Taylor, M. J., and Thompson, R. E.: Assessment of the quality of the version 1.07 temperatureversuspressure profiles of the middle atmosphere from TIMED/SABER, J. Geophys. Res.Atmos., 113, 1–27, https://doi.org/10.1029/2008JD010013, 2008.
Russell III, J. M., Mlynczak, M. G., Gordley, L. L., Tansock Jr., J. J., and Esplin, R. W.: Overview of the SABER experiment and preliminary calibration results, Proc. SPIE, 3756, 277–288, https://doi.org/10.1117/12.366382, 1999.
Schmidt, T., Alexander, P., and de la Torre, A.: Stratospheric gravity wave momentum flux from radio occultations, J. Geophys. Res.Atmos., 121, 4443–4467, https://doi.org/10.1002/2015JD024135, 2016.
Shepherd, G. G., Thuillier, G., Cho, Y., Duboin, M., Evans, W. F. J., Gault, W. A., Hersom, C., Kendall, D. J. W., Lathuillère, C., Lowe, R. P., Mcdade, I. C., Rochon, Y. J., Shepherd, M. G., Solheim, B. H., Wang, D., and Ward, W. E.: The Wind Imaging Interferometer (WINDII) on the Upper Atmosphere Research Satellite: A 20 year perspective, Rev. Geophys., 50, 1–38, https://doi.org/10.1029/2012RG000390, 2012.
Shinagawa, H., Miyoshi, Y., Jin, H., and Fujiwara, H.: Global distribution of neutral wind shear associated with sporadic E layers derived from GAIA, J. Geophys. Res.Space, 122, 4450–4465, https://doi.org/10.1002/2016JA023778, 2017.
Stevens, M. H., Lossow, S., Siskind, D. E., Meier, R. R., Randall, C. E., Russell, J. M., Urban, J., and Murtagh, D.: Space shuttle exhaust plumes in the lower thermosphere: Advective transport and diffusive spreading, J. Atmos. Sol.Terr. Phy., 108, 50–60, https://doi.org/10.1016/j.jastp.2013.12.004, 2014.
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, 1–9, https://doi.org/10.1029/2002jd003135, 2003.
Torrence, C. and Compo, G. P.: A Practical Guide to Wavelet Analysis, B. Am. Meteorol. Soc., 79, 61–78, https://doi.org/10.1175/15200477(1998)079<0061:APGTWA>2.0.CO;2, 1998.
Vadas, S. L. and Liu, H. L.: Numerical modeling of the largescale neutral and plasma responses to the body forces created by the dissipation of gravity waves from 6 h of deep convection in Brazil, J. Geophys. Res.Space., 118, 2593–2617, https://doi.org/10.1002/jgra.50249, 2013.
Wang, L. and Alexander, M. J.: Global estimates of gravity wave parameters from GPS radio occultation temperature data, J. Geophys. Res.Atmos., 115, 1–12, https://doi.org/10.1029/2010JD013860, 2010.
Xu, J., Liu, H.L., Yuan, W., Smith, A. K., Roble, R. G., Mertens, C. J., Roble, R. G., Mertens, C. G., Russell III, I. M., and Mlynczak, M. G.: Mesopause structure from Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED)/Sounding of the Atmosphere Using Broadband Emission Radiometry (SABER) observations. J. Geophys. Res.Atmos., 112, D09102, https://doi.org/10.1029/2006JD007711, 2007.
Yu, B., Xue, X., Yue, X., Yang, C., Yu, C., Dou, X., Ning, B., and Hu, L.: The global climatology of the intensity of the ionospheric sporadic E layer, Atmos. Chem. Phys., 19, 4139–4151, https://doi.org/10.5194/acp1941392019, 2019.
Yue, J., She, C.Y., and Liu, H.L.: Large wind shears and stabilities in the mesopause region observed by Na windtemperature lidar at midlatitude, J. Geophys. Res.Atmos., 115, A10307, https://doi.org/10.1029/2009JA014864, 2010.
Yue, J., Liu, H. L., Meier, R. R., Chang, L., Gu, S. Y., and Russell, J.: On the fast zonal transport of the STS121 space shuttle exhaust plume in the lower thermosphere, J. Atmos. Sol.Terr. Phy., 94, 19–27, https://doi.org/10.1016/j.jastp.2012.12.017, 2013.
Zhang, S. P., Mclandress, C., and Shepherd, G. G.: Satellite observations of mean winds and tides in the lower thermosphere: 2. Wind Imaging Interferometer monthly winds for 1992 and 1993, J. Geophys. Res.Atmos., 112, 1–14, https://doi.org/10.1029/2007JD008457, 2007.
 Abstract
 Introduction
 Method of deriving GWperturbed wind shears and validations
 Comparisons of GWperturbed shears with model and observational results
 Climatology of GWperturbed shears
 Discussion
 Summary
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Method of deriving GWperturbed wind shears and validations
 Comparisons of GWperturbed shears with model and observational results
 Climatology of GWperturbed shears
 Discussion
 Summary
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References