Gravity Wave Parameters Research Articles

Fluctuations of radio occultation signals in sounding the Earth's atmosphere

Abstract. We discuss the relationships that link the observed fluctuation spectra of the amplitude and phase of signals used for the radio occultation sounding of the Earth's atmosphere, with the spectra of atmospheric inhomogeneities. Our analysis employs the approximation of the phase screen and of weak fluctuations. We make our estimates for the following characteristic inhomogeneity types: (1) the isotropic Kolmogorov turbulence and (2) the anisotropic saturated internal gravity waves. We obtain the expressions for the variances of the amplitude and phase fluctuations of radio occultation signals as well as their estimates for the typical parameters of inhomogeneity models. From the GPS/MET observations, we evaluate the spectra of the amplitude and phase fluctuations in the altitude interval from 4 to 25 km in the middle and polar latitudes. As indicated by theoretical and experimental estimates, the main contribution into the radio signal fluctuations comes from the internal gravity waves. The influence of the Kolmogorov turbulence is negligible. We derive simple relationships that link the parameters of internal gravity waves and the statistical characteristics of the radio signal fluctuations. These results may serve as the basis for the global monitoring of the wave activity in the stratosphere and upper troposphere.

First Results From the Ionospheric Extension of WACCM‐X During the Deep Solar Minimum Year of 2008

AbstractNew ionosphere and electrodynamics modules have been incorporated in the thermosphere and ionosphere eXtension of the Whole Atmosphere Community Climate Model (WACCM‐X), in order to self‐consistently simulate the coupled atmosphere‐ionosphere system. The first specified dynamics WACCM‐X v.2.0 results are compared with several data sets, and with the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIE‐GCM), during the deep solar minimum year. Comparisons with Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite of temperature and zonal wind in the lower thermosphere show that WACCM‐X reproduces the seasonal variability of tides remarkably well, including the migrating diurnal and semidiurnal components and the nonmigrating diurnal eastward propagating zonal wavenumber 3 component. There is overall agreement between WACCM‐X, TIE‐GCM, and vertical drifts observed by the Communication/Navigation Outage Forecast System (C/NOFS) satellite over the magnetic equator, but apparent discrepancies also exist. Both model results are dominated by diurnal variations, while C/NOFS observed vertical plasma drifts exhibit strong temporal variations. The climatological features of ionospheric peak densities and heights (NmF2 and hmF2) from WACCM‐X are in general agreement with the results derived from Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) data, although the WACCM‐X predicted NmF2 values are smaller, and the equatorial ionization anomaly crests are closer to the magnetic equator compared to COSMIC and ionosonde observations. This may result from the excessive mixing in the lower thermosphere due to the gravity wave parameterization. These data‐model comparisons demonstrate that WACCM‐X can capture the dynamic behavior of the coupled atmosphere and ionosphere in a climatological sense.

Gravity Wave Dynamics in a Mesospheric Inversion Layer: 2. Instabilities, Turbulence, Fluxes, and Mixing

AbstractA companion paper by Fritts, Laughman, et al. (2017) employed an anelastic numerical model to explore the dynamics of gravity waves (GWs) encountering a mesospheric inversion layer (MIL) having a moderate static stability enhancement and a layer of weaker static stability above. That study revealed that MIL responses, including GW transmission, reflection, and instabilities, are sensitive functions of GW parameters. This paper expands on two of the Fritts, Laughman, et al. (2017) simulations to examine GW instability dynamics and turbulence in the MIL; forcing of the mean wind and stability environments by GW, instability, and turbulence fluxes; and associated heat and momentum transports. These direct numerical simulations resolve turbulence inertial‐range scales and yield the following results: GW breaking and turbulence in the MIL occur below where they would otherwise, due to enhancements of GW amplitudes and shears in the MIL. 2‐D GW and instability heat and momentum fluxes are ~20–30 times larger than 3‐D instability and turbulence fluxes. Mean fields are driven largely by 2‐D GW and instability dynamics rather than 3‐D instabilities and turbulence. 2‐D and 3‐D heat fluxes in regions of strong turbulence yield small departures from initial T(z) and N2(z) profiles, hence do not yield nearly adiabatic “mixed” layers. Our MIL results are consistent with the relation between the turbulent vertical velocity variance and energy dissipation rate proposed by Weinstock (1981) for the limited intervals evaluated.

Retrieving characteristics of inertia gravity wave parameters with least uncertainties using the hodograph method

Abstract. We have analyzed wind velocities measured with high-resolution Global Positioning System (GPS) radiosondes which have been flown continuously from Hyderabad for 120 h with an interval of 6 h. Hodograph method has been used to retrieve inertia gravity waves (IGWs) parameters. Background winds are removed from the time series by detrending, whereas polynomials of different orders are removed to obtain the fluctuations from individual profiles. A Butterworth filter is used to extract the monochromatic IGW component. Another filter finite impulse response (FIR1) is tried in a similar manner to test the effects of filters in estimating IGW characteristics. Results reveal that the fluctuation profiles differ with the change in polynomial orders, but the IGW parameters remain same when a Butterworth filter is chosen to extract the monochromatic wave component. The FIR1 filter produces results with a broader range. The direction of wave propagation can be confirmed with additional temperature information.

Numerical simulation of wave interactions during sudden stratospheric warming

Parameterizations of normal atmospheric modes (NAMs) and orographic gravity waves (OGWs) are implemented into the mechanistic general circulation model of the middle and upper atmosphere (MUA). Numerical experiments of sudden stratospheric warming (SSW) events are performed for climatological conditions typical for January and February using meteorological reanalysis data from the UK MET Office in the MUA model averaged over the years 1992–2011 with the easterly phase of quasi-biennial oscillation (QBO). The simulation shows that an increase in the OGW amplitudes occurs at altitudes higher than 30 km in the Northern Hemisphere after SSW. The OGW amplitudes have maximums at altitudes of about 50 km over the North American and European mountain systems before and during SSW, as well as over the Himalayas after SSW. At high latitudes of the Northern Hemisphere, significant (up to 50–70%) variations in the amplitudes of stationary planetary waves (SPWs) are observed during and after the SSW. Westward travelling NAMs have local amplitude maximums not only in the Northern Hemisphere, but also in the Southern Hemisphere, where there are waveguides for the propagation of these modes. Calculated variations of SPW and NAM amplitudes correspond to changes in the mean temperature and wind fields, as well as the Eliassen-Palm flux and atmospheric refractive index for the planetary waves, during SSW. Including OGW thermal and dynamical effects leads to an increase in amplitude (by 30–70%) of almost all SPWs before and during SSW and to a decrease (up to 20–100%) after the SSW at middle and high latitudes of the Northern Hemisphere.

Momentum Flux of Convective Gravity Waves Derived from an Offline Gravity Wave Parameterization. Part I: Spatiotemporal Variations at Source Level

Abstract Spatiotemporal variations in momentum flux spectra of convective gravity waves (CGWs) at the source level (cloud top), including nonlinear forcing effects, are examined based on calculations using an offline version of CGW parameterization and global reanalysis data for a period of 32 years (1979–2010). The cloud-top momentum flux (CTMF) is not solely proportional to the convective heating rate but is affected by the wave-filtering and resonance factor and background stability and temperature underlying the convection. Consequently, the primary peak of CTMF is in the winter hemisphere midlatitudes, associated with storm tracks, where a secondary peak of convective heating exists, whereas the secondary peak of CTMF appears in the summer hemisphere tropics and intertropical convergence zone (ITCZ), where the primary peak of convective heating exists. The magnitude of CTMF fluctuates largely with 1-yr and 1-day periods in major CTMF regions. At low latitudes and Pacific storm-track regions, a 6-month period is also significant, and the decadal cycle appears in the southern Andes. The equatorial eastern Pacific region exhibits a substantial interannual to decadal scale of variabilities. The correlation between convective heating and the CTMF is relatively lower in the equatorial region than in other regions. The CTMF in 10°N–10°S during the period of the pre-Concordiasi campaign approximately follows a lognormal distribution but with a slight underestimation in the tail of the probability density function. In Part II, the momentum flux and drag of CGW in the stratosphere will be examined.

Gravity Wave Spectra in the Lower Stratosphere Diagnosed from Project Loon Balloon Trajectories.

Project Loon has been launching super-pressure balloons since January 2013 to provide worldwide Internet coverage. These balloons typically fly between 18-21 km and provide measurements of winds and pressure fluctuations in the lower stratosphere. We divide 1,560 Loon flights into 3,405 two-day segments for gravity wave analysis. We derive the kinetic energy spectrum from the horizontal balloon motion and estimate the temperature perturbation spectrum (proportional to the potential energy spectrum) from the pressure variations. We fit the temperature (and kinetic energy) data to the functional form where ω is the wave frequency, ω o is daily frequency, T' o is the base temperature amplitude and α is the slope. Both the kinetic energy and temperature spectra show -1.9±0.2 power-law dependence in the intrinsic frequency window 3 - 50 cycles/day. The temperature spectrum slope is weakly anti-correlated with the base temperature amplitude. We also find that the wave base temperature distribution is highly skewed. The average tropical modal temperature is 0.77 K. The highest amplitude waves occur over the mountainous regions, the tropics, and the high southern latitudes. Temperature amplitudes show little height variation over our 18-21 km domain. Our results are consistent with other limited super-pressure balloon analyses. The modal temperature is higher than the temperature currently used in Lagrangian model gravity wave parameterizations.

Near InfraRed Imaging Spectrograph (NIRIS) for ground-based mesospheric OH(6-2) and $$\hbox {O}_{2}$$ O 2 (0-1) intensity and temperature measurements

This paper describes the development of a new Near InfraRed Imaging Spectrograph (NIRIS) which is capable of simultaneous measurements of OH(6-2) Meinel and $$\hbox {O}_{2}$$ (0-1) atmospheric band nightglow emission intensities. In this spectrographic technique, rotational line ratios are obtained to derive temperatures corresponding to the emission altitudes of 87 and 94 km. NIRIS has been commissioned for continuous operation from optical aeronomy observatory, Gurushikhar, Mount Abu ( $$24.6^{\circ }\hbox {N}$$ , $$72.8^{\circ }\hbox {E}$$ ) since January 2013. NIRIS uses a diffraction grating of 1200 lines $$\hbox {mm}^{-1}$$ and 1024 $$\times $$ 1024 pixels thermoelectrically cooled CCD camera and has a large field-of-view (FOV) of $$80^{\circ }$$ along the slit orientation. The data analysis methodology adopted for the derivation of mesospheric temperatures is also described in detail. The observed NIRIS temperatures show good correspondence with satellite (SABER) derived temperatures and exhibit both tidal and gravity waves (GW) like features. From the time taken for phase propagation in the emission intensities between these two altitudes, vertical phase speed of gravity waves, $$c_{z}$$ , is calculated and along with the coherent GW time period ‘ $$\tau $$ ’, the vertical wavelength, $$\lambda _{z}$$ , is obtained. Using large FOV observations from NIRIS, the meridional wavelengths, $$\lambda _{y}$$ , are also calculated. We have used one year of data to study the possible cause(s) for the occurrences of mesospheric temperature inversions (MTIs). From the statistics obtained for 234 nights, it appears that in situ chemical heating is mainly responsible for the observed MTIs than the vertical propagation of the waves. Thus, this paper describes a novel near infrared imaging spectrograph, its working principle, data analysis method for deriving OH and $$\hbox {O}_{2}$$ emission intensities and the corresponding rotational temperatures at these altitudes, derivation of gravity wave parameters ( $$\tau $$ , $$c_{z}$$ , $$\lambda _{z}$$ , and $$\lambda _{y})$$ , and results on the statistical study of MTIs that exist in the earth’s mesospheric altitudes.

On the influence of zonal gravity wave distributions on the Southern Hemisphere winter circulation

Abstract. A mechanistic global circulation model is used to simulate the Southern Hemisphere stratospheric, mesospheric, and lower thermospheric circulation during austral winter. The model includes a gravity wave (GW) parameterization that is initiated by prescribed 2-D fields of GW parameters in the troposphere. These are based on observations of GW potential energy calculated using GPS radio occultations and show enhanced GW activity east of the Andes and around the Antarctic. In order to detect the influence of an observation-based and thus realistic 2-D GW distribution on the middle atmosphere circulation, we perform model experiments with zonal mean and 2-D GW initialization, and additionally with and without forcing of stationary planetary waves (SPWs) at the lower boundary of the model. As a result, we find additional forcing of SPWs in the stratosphere, a weaker zonal wind jet in the mesosphere, cooling of the mesosphere and warming near the mesopause above the jet. SPW wavenumber 1 (SPW1) amplitudes are generally increased by about 10 % when GWs are introduced being longitudinally dependent. However, at the upper part of the zonal wind jet, SPW1 in zonal wind and GW acceleration are out of phase, which reduces the amplitudes there.

Exploring gravity wave characteristics in 3-D using a novel S-transform technique: AIRS/Aqua measurements over the Southern Andes and Drake Passage

Abstract. Gravity waves (GWs) transport momentum and energy in the atmosphere, exerting a profound influence on the global circulation. Accurately measuring them is thus vital both for understanding the atmosphere and for developing the next generation of weather forecasting and climate prediction models. However, it has proven very difficult to measure the full set of GW parameters from satellite measurements, which are the only suitable observations with global coverage. This is particularly critical at latitudes close to 60° S, where climate models significantly under-represent wave momentum fluxes. Here, we present a novel fully 3-D method for detecting and characterising GWs in the stratosphere. This method is based around a 3-D Stockwell transform, and can be applied retrospectively to existing observed data. This is the first scientific use of this spectral analysis technique. We apply our method to high-resolution 3-D atmospheric temperature data from AIRS/Aqua over the altitude range 20–60 km. Our method allows us to determine a wide range of parameters for each wave detected. These include amplitude, propagation direction, horizontal/vertical wavelength, height/direction-resolved momentum fluxes (MFs), and phase and group velocity vectors. The latter three have not previously been measured from an individual satellite instrument. We demonstrate this method over the region around the Southern Andes and Antarctic Peninsula, the largest known sources of GW MFs near the 60° S belt. Our analyses reveal the presence of strongly intermittent highly directionally focused GWs with very high momentum fluxes (∼ 80–100 mPa or more at 30 km altitude). These waves are closely associated with the mountains rather than the open ocean of the Drake Passage. Measured fluxes are directed orthogonal to both mountain ranges, consistent with an orographic source mechanism, and are largest in winter. Further, our measurements of wave group velocity vectors show clear observational evidence that these waves are strongly focused into the polar night wind jet, and thus may contribute significantly to the missing momentum at these latitudes. These results demonstrate the capabilities of our new method, which provides a powerful tool for delivering the observations required for the next generation of weather and climate models.

Mean-Flow Effects of Thermal Tides in the Mesosphere and Lower Thermosphere

Abstract This study addresses the heat budget of the mesosphere and lower thermosphere with regard to the energy deposition of upward-propagating waves. To this end, the energetics of gravity waves are recapitulated using an anelastic version of the primitive equations. This leads to an expression for the energy deposition of waves that is usually resolved in general circulation models. The energy deposition is shown to be mainly due to the frictional heating and, additionally, due to the negative buoyancy production of wave kinetic energy. The frictional heating includes contributions from horizontal and vertical momentum diffusion, as well as from ion drag. This formalism is applied to analyze results from a mechanistic middle-atmosphere general circulation model that includes energetically consistent parameterizations of diffusion, gravity waves, and ion drag. This paper estimates 1) the wave driving and energy deposition of thermal tides, 2) the model response to the excitation of thermal tides, and 3) the model response to the combined energy deposition by parameterized gravity waves and resolved waves. It is found that thermal tides give rise to a significant energy deposition in the lower thermosphere. The temperature response to thermal tides is positive. It maximizes at polar latitudes in the lower thermosphere as a result of poleward circulation branches that are driven by the predominantly westward Eliassen–Palm flux divergence of the tides. In addition, thermal tides give rise to a downward shift and reduction of the gravity wave drag in the upper mesosphere. Including the energy deposition in the model causes a substantial warming in the upper mesosphere and lower thermosphere.

Revisiting Southern Hemisphere polar stratospheric temperature trends in WACCM: The role of dynamical forcing

AbstractThe latest version of the Whole Atmosphere Community Climate Model (WACCM), which includes a new chemistry scheme and an updated parameterization of orographic gravity waves, produces temperature trends in the Antarctic lower stratosphere in excellent agreement with radiosonde observations for 1969–1998 as regards magnitude, location, timing, and persistence. The maximum trend, reached in November at 100 hPa, is −4.4 ± 2.8 K decade−1, which is a third smaller than the largest trend in the previous version of WACCM. Comparison with a simulation without the updated orographic gravity wave parameterization, together with analysis of the model's thermodynamic budget, reveals that the reduced trend is due to the effects of a stronger Brewer‐Dobson circulation in the new simulations, which warms the polar cap. The effects are both direct (a trend in adiabatic warming in late spring) and indirect (a smaller trend in ozone, hence a smaller reduction in shortwave heating, due to the warmer environment).

Influence of parameterized small‐scale gravity waves on the migrating diurnal tide in Earth's thermosphere

AbstractEffects of subgrid‐scale gravity waves (GWs) on the diurnal migrating tides are investigated from the mesosphere to the upper thermosphere for September equinox conditions, using a general circulation model coupled with the extended spectral nonlinear GW parameterization of Yiğit et al. (). Simulations with GW effects cut off above the turbopause and included in the entire thermosphere have been conducted. GWs appreciably impact the mean circulation and cool the thermosphere down by up to 12–18%. GWs significantly affect the winds modulated by the diurnal migrating tide, in particular, in the low‐latitude mesosphere and lower thermosphere and in the high‐latitude thermosphere. These effects depend on the mutual correlation of the diurnal phases of the GW forcing and tides: GWs can either enhance or reduce the tidal amplitude. In the low‐latitude MLT, the correlation between the direction of the deposited GW momentum and the tidal phase is positive due to propagation of a broad spectrum of GW harmonics through the alternating winds. In the Northern Hemisphere high‐latitude thermosphere, GWs act against the tide due to an anticorrelation of tidal wind and GW momentum, while in the Southern high‐latitudes they weakly enhance the tidal amplitude via a combination of a partial correlation of phases and GW‐induced changes of the circulation. The variable nature of GW effects on the thermal tide can be captured in GCMs provided that a GW parameterization (1) considers a broad spectrum of harmonics, (2) properly describes their propagation, and (3) correctly accounts for the physics of wave breaking/saturation.

Seasonal variation of gravity wave parameters using different filter methods with daylight lidar measurements at midlatitudes

AbstractThe daylight‐capable Rayleigh‐Mie‐Raman (RMR) lidar at the midlatitude station in Kühlungsborn (54°N, 12°E) is in operation since 2010. The RMR lidar system is used to investigate different fractions of atmospheric waves, like gravity waves (GW) and thermal tides (with diurnal, semidiurnal, and terdiurnal components) at day and night. About 6150 h of data have been acquired until 2015. The general challenge for GW observations is the separation of different wave contributions from the observed superposition of GW, tides, or even longer periodic waves. Unfiltered lidar data always include such a superposition. We applied a Butterworth filter to separate GW and tides by vertical wavelength with a cutoff wavelength of 15 km and by observed periods with a cutoff period of 8 h. GW activity and characteristics are derived in an altitude range between 30 and 70 km. The retrieved vertically filtered temperature deviations contain GW with small vertical wavelengths over a broad range of periods, while only a small range of periods is included in the temporally filtered temperature deviations. We observe an annual variation of the wave activity for unfiltered and vertically filtered data, which is caused from tides and inertia gravity waves. In contrast to that, filtering in time leads to a weak semiannual variation for gravity waves with periods of 4–8 h, especially in higher altitudes. During summer, these waves have the half of the total amount of the potential energy budget compared to the unfiltered data. This shows the importance of waves with periods smaller than 8 h.

On the Dynamical Control of the Mesosphere–Lower Thermosphere by the Lower and Middle Atmosphere

Abstract The NCAR Whole Atmosphere Community Climate Model (WACCM) is used to investigate the dynamical influence of the lower and middle atmosphere on the upper mesosphere and lower thermosphere. In simulations using a methodology adapted from the “specified dynamics” (nudged) version of the model, horizontal winds and temperature over part of the vertical range of the atmosphere are relaxed toward results from a previous simulation that serves as the true simulation, equivalent to meteorological analysis. In the upper mesosphere, the magnitude of the divergence of the constrained simulations from the true simulation depends on the vertical extent and frequency of the data used for nudging the model and grows with altitude. The simulations quantify the error growth of the model dynamical fields when data and forcing terms are known exactly and there are no model biases. The error growth rate and the ultimate discrepancy between the nudged and true fields depend strongly on the method used for representing gravity wave drag. The largest error growth occurs when the gravity wave parameterization uses interactive wave sources that depend on convective activity or fronts. Errors are reduced when the same parameterization is used with smoothly varying specified wave sources. The smallest errors are seen when the parameterized gravity wave drag is replaced by linear Rayleigh friction damping on the wind speed. These comparisons demonstrate the role of gravity waves in transporting the variability of the troposphere into the mesosphere and lower thermosphere.

Sensitivity of the quasi‐biennial oscillation simulated in WACCM to the phase speed spectrum and the settings in an inertial gravity wave parameterization

AbstractThe application of inertial gravity wave parameterization has allowed for the spontaneous generation of quasi‐biennial oscillation (QBO) in the Whole Atmosphere Community Climate Model (WACCM), although there is some mismatch when comparing with observations. The parameterization is based on Lindzen's linear saturation theory, modified to describe inertia‐gravity waves (IGW) by considering the Coriolis effect. In this work, we improve the parameterization by importing a more realistic IGW phase speed spectrum that exhibits a double peak Gaussian distribution calculated from tropical radiosonde observations. A series of numeric simulations are performed to test the sensitivity of QBO‐like oscillation features to the phase speed spectrum and the settings of parameterized IGW. All these simulations are capable of generating equatorial wind oscillations in the stratosphere based on standard spatial resolution settings. Central phase speeds of the “double‐Gaussian parameterization” affect QBO magnitudes and periods, and the momentum flux of IGW determines the acceleration rate of zonal wind. Furthermore, stronger IGW forcing can lead to a propagation of the QBO‐like oscillation to lower altitude. The intermittency factor of the parameterization also prominently affects the QBO period. Stratospheric QBO‐like oscillation with obvious improvements is generated using the new IGW parameterization in a long‐time simulation.

Middle atmosphere dynamical sources of the semiannual oscillation in the thermosphere and ionosphere

AbstractThe strong global semiannual oscillation (SAO) in thermospheric density has been observed for five decades, but definitive knowledge of its source has been elusive. We use the National Center of Atmospheric Research thermosphere‐ionosphere‐mesosphere electrodynamics general circulation model (TIME‐GCM) to study how middle atmospheric dynamics generate the SAO in the thermosphere‐ionosphere (T‐I). The “standard” TIME‐GCM simulates, from first principles, SAOs in thermospheric mass density and ionospheric total electron content that agree well with observed climatological variations. Diagnosis of the globally averaged continuity equation for atomic oxygen ([O]) shows that the T‐I SAO originates in the upper mesosphere, where an SAO in [O] is forced by nonlinear, resolved‐scale variations in the advective, net tidal, and diffusive transport of O. Contrary to earlier hypotheses, TIME‐GCM simulations demonstrate that intra‐annually varying eddy diffusion by breaking gravity waves may not be the primary driver of the T‐I SAO: A pronounced SAO is produced without parameterized gravity waves.

Modification of the Gravity Wave Parameterization in the Whole Atmosphere Community Climate Model: Motivation and Results

Abstract The current standard version of the Whole Atmosphere Community Climate Model (WACCM) simulates Southern Hemisphere winter and spring temperatures that are too cold compared with observations. This “cold-pole bias” leads to unrealistically low ozone column amounts in Antarctic spring. Here, the cold-pole problem is addressed by introducing additional mechanical forcing of the circulation via parameterized gravity waves. Insofar as observational guidance is ambiguous regarding the gravity waves that might be important in the Southern Hemisphere stratosphere, the impact of increasing the forcing by orographic gravity waves was investigated. This reduces the strength of the Antarctic polar vortex in WACCM, bringing it into closer agreement with observations, and accelerates the Brewer–Dobson circulation in the polar stratosphere, which warms the polar cap and improves substantially the simulation of Antarctic temperature. These improvements are achieved without degrading the performance of the model in the Northern Hemisphere stratosphere or in the mesosphere and lower thermosphere of either hemisphere. It is shown, finally, that other approaches that enhance gravity wave forcing can also reduce the cold-pole bias such that careful examination of observational evidence and model performance will be required to establish which gravity wave sources are dominant in the real atmosphere. This is especially important because a “downward control” analysis of these results suggests that the improvement of the cold-pole bias itself is not very sensitive to the details of how gravity wave drag is altered.

Impacts of SABER CO2‐based eddy diffusion coefficients in the lower thermosphere on the ionosphere/thermosphere

AbstractThis work estimates global‐mean Kzz using Sounding of the Atmosphere using Broadband Emission Radiometry/Thermosphere‐Ionosphere‐Mesosphere Energetics and Dynamics monthly global‐mean CO2 profiles and a one‐dimensional transport model. It is then specified as a lower boundary into the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIE‐GCM). Results first show that global‐mean CO2 in the mesosphere and lower thermosphere region has annual and semiannual oscillations (AO and SAO) with maxima during solstice seasons along with a primary maximum in boreal summer. Our calculated AO and SAO in global‐mean CO2 are then modeled by AO and SAO in global‐mean Kzz. It is then shown that our estimated global‐mean Kzz is lower in magnitude than the suggested global‐mean Kzz from Qian et al. (2009) that can model the observed AO and SAO in the ionosphere/thermosphere (IT) region. However, our estimated global‐mean Kzz is similar in magnitude with recent suggestions of global‐mean Kzz in models with explicit gravity wave parameterization. Our work therefore concludes that global‐mean Kzz from global‐mean CO2 profiles cannot model the observed AO and SAO in the IT region because our estimated global‐mean Kzz may only be representing eddy diffusion due to gravity wave breaking. The difference between our estimated global‐mean Kzz and the global‐mean Kzz from Qian et al. (2009) thus represents diffusion and mixing from other nongravity wave sources not directly accounted for in the TIE‐GCM lower boundary conditions. These other sources may well be the more dominant lower atmospheric forcing behind the AO and SAO in the IT region.

Global Distribution of Gravity Wave Sources and Fields in the Martian Atmosphere during Equinox and Solstice Inferred from a High-Resolution General Circulation Model

Abstract Results of simulations with a new high-resolution Martian general circulation model (MGCM) (T106 spectral resolution, or ~67-km horizontal grid size) have been analyzed to reveal global distributions of gravity waves (GWs) during the solstice and equinox periods. They show that shorter-scale harmonics progressively dominate with height, and the body force per unit mass (drag) they impose on the larger-scale flow increases. Mean magnitudes of the drag in the middle atmosphere are tens of meters per second per sol, while instantaneously they can reach thousands of meters per second per sol. Inclusion of small-scale GW harmonics results in an attenuation of the wind jets in the middle atmosphere and in the tendency of their reversal. GW energy in the troposphere due to the shortest-scale harmonics is concentrated in the low latitudes for both seasons and is in a good agreement with observations. The vertical fluxes of wave horizontal momentum are directed mainly against the larger-scale wind. Orographically generated GWs contribute significantly to the total energy of small-scale disturbances and to the drag created by the latter. These waves strongly decay with height, and thus the nonorographic GWs of tropospheric origin dominate near the mesopause. The results of this study can be used to better constrain and validate GW parameterizations in MGCMs.