Gravity Wave Momentum Flux Research Articles

The Mesoscale Gravity Wave Response to the 2022 Tonga Volcanic Eruption: AIRS and MLS Satellite Observations and Source Backtracing

AbstractOn 15 January 2022, the Hunga Tonga–Hunga Ha'apai volcano erupted violently. This exceptional event excited a manifold of atmospheric waves. Here, we focus on the mesoscale part of the wave spectrum. About 8.5 hr after the eruption a strong atmospheric gravity wave (GW) was observed in the stratosphere by the satellite instruments Atmospheric Infrared Sounder (AIRS) and Microwave Limb Sounder (MLS) in the vicinity of Tonga. By ray‐tracing, we confirm the eruption as the source of this GW event. We determine the wave characteristics of the GW in terms of horizontal and vertical wavelengths and GW momentum flux. The strength of the GW is compared to the usual Southern Hemisphere flux values during this week. The event is comparable to the strongest convective events considering MLS, and exceptionally strong considering AIRS, which observes faster waves only.

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

AbstractDue to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large‐scale wind profiles from the ERA‐Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

Propagation paths and source distributions of resolved gravity waves in ECMWF-IFS analysis fields around the southern polar night jet

Abstract. In the southern winter polar stratosphere, the distribution of gravity wave momentum flux in many state-of-the-art climate simulations is inconsistent with long-time satellite and superpressure balloon observations around 60∘ S. Recent studies hint that a lateral shift between prominent gravity wave sources in the tropospheric mid-latitudes and the location where gravity wave activity is present in the stratosphere causes at least part of the discrepancy. This lateral shift cannot be represented by the column-based gravity wave drag parameterisations used in most general circulation models. However, recent high-resolution analysis and re-analysis products of the European Centre for Medium-Range Weather Forecasts Integrated Forecast System (ECMWF-IFS) show good agreement with the observations and allow for a detailed investigation of resolved gravity waves, their sources, and propagation paths. In this paper, we identify resolved gravity waves in the ECMWF-IFS analyses for a case of high gravity wave activity in the lower stratosphere using small-volume sinusoidal fits to characterise these gravity waves. The 3D wave vector together with perturbation amplitudes, wave frequency, and a fully described background atmosphere are then used to initialise the Gravity Wave Regional or Global Ray Tracer (GROGRAT) gravity wave ray tracer and follow the gravity waves backwards from the stratosphere. Finally, we check for the indication of source processes on the path of each ray and, thus, quantitatively attribute gravity waves to sources that are represented within the model. We find that stratospheric gravity waves are indeed subject to far (>1000 km) lateral displacement from their sources, which already take place at low altitudes (<20 km). Various source processes can be linked to waves within stratospheric gravity wave (GW) patterns, such as the orography equatorward of 50∘ S and non-orographic sources above the Southern Ocean. These findings may explain why superpressure balloons observe enhanced gravity wave momentum fluxes in the lower stratosphere over the Southern Ocean despite an apparent lack of sources at this latitude. Our results also support the need to improve gravity wave parameterisations to account for meridional propagation.

Numerical Simulation of Mountain Waves over the Southern Andes. Part II: Momentum Fluxes and Wave–Mean-Flow Interactions

AbstractA companion paper by Lund et al. employed a compressible model to describe the evolution of mountain waves arising due to increasing flow with time over the southern Andes, their breaking, secondary gravity waves and acoustic waves arising from these dynamics, and their local responses. This paper describes the mountain wave, secondary gravity wave, and acoustic wave vertical fluxes of horizontal momentum, and the local and large-scale three-dimensional responses to gravity breaking and wave–mean-flow interactions accompanying this event. Mountain wave and secondary gravity wave momentum fluxes and deposition vary strongly in space and time due to variable large-scale winds and spatially localized mountain wave and secondary gravity wave responses. Mountain wave instabilities accompanying breaking induce strong, local, largely zonal forcing. Secondary gravity waves arising from mountain wave breaking also interact strongly with large-scale winds at altitudes above ~80 km. Together, these mountain wave and secondary gravity wave interactions reveal systematic gravity wave–mean-flow interactions having implications for both mean and tidal forcing and feedbacks. Acoustic waves likewise achieve large momentum fluxes, but typically imply significant responses only at much higher altitudes.

The semiannual oscillation (SAO) in the tropical middle atmosphere and its gravity wave driving in reanalyses and satellite observations

Abstract. Gravity waves play a significant role in driving the semiannual oscillation (SAO) of the zonal wind in the tropics. However, detailed knowledge of this forcing is missing, and direct estimates from global observations of gravity waves are sparse. For the period 2002–2018, we investigate the SAO in four different reanalyses: ERA-Interim, JRA-55, ERA-5, and MERRA-2. Comparison with the SPARC zonal wind climatology and quasi-geostrophic winds derived from Microwave Limb Sounder (MLS) and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite observations show that the reanalyses reproduce some basic features of the SAO. However, there are also large differences, depending on the model setup. Particularly, MERRA-2 seems to benefit from dedicated tuning of the gravity wave drag parameterization and assimilation of MLS observations. To study the interaction of gravity waves with the background wind, absolute values of gravity wave momentum fluxes and a proxy for absolute gravity wave drag derived from SABER satellite observations are compared with different wind data sets: the SPARC wind climatology; data sets combining ERA-Interim at low altitudes and MLS or SABER quasi-geostrophic winds at high altitudes; and data sets that combine ERA-Interim, SABER quasi-geostrophic winds, and direct wind observations by the TIMED Doppler Interferometer (TIDI). In the lower and middle mesosphere the SABER absolute gravity wave drag proxy correlates well with positive vertical gradients of the background wind, indicating that gravity waves contribute mainly to the driving of the SAO eastward wind phases and their downward propagation with time. At altitudes 75–85 km, the SABER absolute gravity wave drag proxy correlates better with absolute values of the background wind, suggesting a more direct forcing of the SAO winds by gravity wave amplitude saturation. Above about 80 km SABER gravity wave drag is mainly governed by tides rather than by the SAO. The reanalyses reproduce some basic features of the SAO gravity wave driving: all reanalyses show stronger gravity wave driving of the SAO eastward phase in the stratopause region. For the higher-top models ERA-5 and MERRA-2, this is also the case in the lower mesosphere. However, all reanalyses are limited by model-inherent damping in the upper model levels, leading to unrealistic features near the model top. Our analysis of the SABER and reanalysis gravity wave drag suggests that the magnitude of SAO gravity wave forcing is often too weak in the free-running general circulation models; therefore, a more realistic representation is needed.

Observation of Gravity Waves at the Tropical Tropopause Using Superpressure Balloons

AbstractTropical gravity wave activity is investigated using measurements of momentum fluxes gathered during Strateole‐2 superpressure balloon flights. The data set consists of eight balloon flights performed in the deep tropics from November 2019 to February 2020. The flights lasted for 2–3 months each, and in‐situ meteorological data were collected every 30 s. The relation between gravity waves and deep convection is investigated using geostationary satellite data from the NOAA/NCEP GPM_MERGIR satellite data product, at 1 h resolution. The amplitude of gravity wave momentum fluxes shows a clear dependence on the distance to the nearest convective system, with a strong decay as distance to convection increases. The largest momentum‐flux values (>mPa) are only found less than 200 km away from deep convection. The sensitivity of the wave flux to distance from convection is stronger for high frequency gravity waves (periods shorter than 60 min). Lower frequency waves tend to a non‐zero, background value away from convection, supporting some background value in gravity‐wave drag parameterizations. On the other hand, the wide range of momentum flux values observed close to the convection emphasizes the intermittent nature of the gravity‐wave source. The large scale variation of gravity‐wave intermittency within the equatorial belt is also studied. The results highlight spatial variations of gravity wave activity, with the highest momentum flux recorded over land.

Towards a more “scale‐aware” orographic gravity wave drag parametrization: Description and initial testing

AbstractMany current orographic gravity wave drag parametrizations employed in numerical weather prediction and climate models assume idealised elliptical subgrid orography. As a result, they do not fully capture the spectrum of subgrid orographic scales and do not exhibit the desired behaviour when the model grid length is varied. This article describes the motivation, formulation, and testing of a new orographic gravity wave drag parametrization in the Met Office Unified Model that is more “scale‐aware”. The scheme circumvents the need to make assumptions about the shape and distribution of the subgrid orography, moving away from elliptical mountains, through the use of Fourier transforms. The reduction in gravity wave surface stress when orographic flow blocking is present is also accounted for. The result is a new scheme that represents subgrid mountain scales more faithfully and consistently, rather than assuming a single mountain scale for each grid box. Tests using both regional and global simulations, ranging in grid spacing from 2–130 km, demonstrate that the scale‐aware scheme provides a more constant total, resolved plus parametrized, orographic gravity wave momentum flux, when the grid length is varied over most regions of the globe. Global five‐day forecasts also show a significant reduction in the too‐strong zonal wind bias within the Northern and Southern Hemisphere upper troposphere to lower stratosphere at grid spacings of 130 and 40 km when the scale‐aware scheme is employed. Some discussion on future developments for the scheme is also presented, noting that improvements to the orographic flow‐blocking drag component are required before the scheme can be implemented operationally.

Stratospheric gravity waves over the mountainous island of South Georgia: testing a high-resolution dynamical model with 3-D satellite observations and radiosondes

Abstract. Atmospheric gravity waves (GWs) play an important role in atmospheric dynamics but accurately representing them in general circulation models (GCMs) is challenging. This is especially true for orographic GWs generated by wind flow over small mountainous islands in the Southern Ocean. Currently, these islands lie in the “grey zone” of global model resolution, where they are neither fully resolved nor fully parameterised. It is expected that as GCMs approach the spatial resolution of current high-resolution local-area models, small-island GW sources may be resolved without the need for parameterisations. But how realistic are the resolved GWs in these high-resolution simulations compared to observations? Here, we test a high-resolution (1.5 km horizontal grid, 118 vertical levels) local-area configuration of the Met Office Unified Model over the mountainous island of South Georgia (54∘ S, 36∘ W), running without GW parameterisations. The island's orography is well resolved in the model, and real-time boundary conditions are used for two time periods during July 2013 and June–July 2015. We compare simulated GWs in the model to coincident 3-D satellite observations from the Atmospheric Infrared Sounder (AIRS) on board Aqua. By carefully sampling the model using the AIRS resolution and measurement footprints (denoted as model sampled as AIRS hereafter), we present the first like-for-like comparison of simulated and observed 3-D GW amplitudes, wavelengths and directional GW momentum flux (GWMF) over the island using a 3-D S-transform method. We find that the timing, magnitude and direction of simulated GWMF over South Georgia are in good general agreement with observations, once the AIRS sampling and resolution are applied to the model. Area-averaged zonal GWMF during these 2 months is westward at around 5.3 and 5.6 mPa in AIRS and model sampled as AIRS datasets respectively, but values directly over the island can exceed 50 mPa. However, up to 35 % of the total GWMF in AIRS is actually found upwind of the island compared to only 17 % in the model sampled as AIRS, suggesting that non-orographic GWs observed by AIRS may be underestimated in our model configuration. Meridional GWMF results show a small northward bias (∼20 %) in the model sampled as AIRS that may correspond to a southward wind bias compared to coincident radiosonde measurements. Finally, we present one example of large-amplitude (T′≈15–20 K at 45 km altitude) GWs at short horizontal wavelengths (λH≈30–40 km) directly over the island in AIRS measurements that show excellent agreement with the model sampled as AIRS. This suggests that orographic GWs in the full-resolution model with T′≈45 K and λH≈30–40 km can occur in reality. Our study demonstrates that not only can high-resolution local-area models simulate realistic stratospheric GWs over small mountainous islands but the application of satellite sampling and resolution to these models can also be a highly effective method for their validation.

On the Momentum Flux of Vertically-Propagating Orographic Gravity Waves Excited in Nonhydrostatic Flow over Three-dimensional Orography

AbstractThis work studies nonhydrostatic effects (NHE) on the momentum flux of orographic gravity waves (OGWs) forced by isolated three-dimensional orography. Based on linear wave theory, an asymptotic expression for low horizonal Froude number ( where (U, V) is the mean horizontal wind, γ and a are the orography anisotropy and half-width and N is the buoyancy frequency) is derived for the gravity wave momentum flux (GWMF) of vertically-propagating waves. According to this asymptotic solution, which is quite accurate for any value of Fr, NHE can be divided into two terms (NHE1 and NHE2). The first term contains the high-frequency parts of the wave spectrum that are often mistaken as hydrostatic waves, and only depends on Fr. The second term arises from the difference between the dispersion relationships of hydrostatic and nonhydrostatic OGWs. Having an additional dependency on the horizontal wind direction and orography anisotropy, this term can change the GWMF direction. Examination of NHE for OGWs forced by both circular and elliptical orography reveals that the GWMF is reduced as Fr increases, at a faster rate than for two-dimensional OGWs forced by a ridge. At low Fr, the GWMF reduction is mostly attributed to the NHE2 term, whereas the NHE1 term starts to dominate above about Fr = 0.4. The behavior of NHE is mainly determined by Fr, while horizontal wind direction and orography anisotropy play a minor role. Implications of the asymptotic GWMF expression for the parameterization of nonhydrostatic OGWs in high-resolution and/or variable-resolution models are discussed.

Toward Transient Subgrid-Scale Gravity Wave Representation in Atmospheric Models. Part I: Propagation Model Including Nondissipative Wave–Mean-Flow Interactions

AbstractCurrent gravity wave (GW) parameterization (GWP) schemes are using the steady-state assumption, in which an instantaneous balance between GWs and mean flow is postulated, thereby neglecting transient, nondissipative interactions between the GW field and the resolved flow. These schemes rely exclusively on wave dissipation, by GW breaking or near critical layers, as a mechanism leading to forcing of the mean flow. In a transient GWP, without the steady-state assumption, nondissipative wave–mean-flow interactions are enabled as an additional mechanism. Idealized studies have shown that this is potentially important, and therefore the transient GWP Multiscale Gravity Wave Model (MS-GWaM) has been implemented into a state-of-the-art weather and climate model. In this implementation, MS-GWaM leads to a zonal-mean circulation that agrees well with observations and increases GW momentum-flux intermittency as compared with steady-state GWPs, bringing it into better agreement with superpressure balloon observations. Transient effects taken into account by MS-GWaM are shown to make a difference even on monthly time scales: in comparison with steady-state GWPs momentum fluxes in the lower stratosphere are increased and the amount of missing drag at Southern Hemispheric high latitudes is decreased to a modest but nonnegligible extent. An analysis of the contribution of different wavelengths to the GW signal in MS-GWaM suggests that small-scale GWs play an important role down to horizontal and vertical wavelengths of 50 km (or even smaller) and 200 m, respectively.

Effects of Latitude-Dependent Gravity Wave Source Variations on the Middle and Upper Atmosphere

Atmospheric gravity waves (GWs) are generated in the lower atmosphere by various weather phenomena. They propagate upward, carry energy and momentum to higher altitudes, and appreciably influence the general circulation upon depositing them in the middle and upper atmosphere. We use a three-dimensional first-principle general circulation model (GCM) with implemented nonlinear whole atmosphere GW parameterization to study the global climatology of wave activity and produced effects at altitudes up to the upper thermosphere. The numerical experiments were guided by the GW momentum fluxes and temperature variances as measured in 2010 by the SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) instrument onboard NASA’s TIMED (Thermosphere Ionosphere Mesosphere Energetics Dynamics) satellite. This includes the latitudinal dependence and magnitude of GW activity in the lower stratosphere for the boreal summer season. The modeling results were compared to the SABER temperature and total absolute momentum flux and Upper Atmosphere Research Satellite (UARS) data in the mesosphere and lower thermosphere. Simulations suggest that, in order to reproduce the observed circulation and wave activity in the middle atmosphere, GW fluxes that are smaller than observed fluxes have to be used at the source level in the lower atmosphere. This is because observations contain a broader spectrum of GWs, while parameterizations capture only a portion relevant to the middle and upper atmosphere dynamics. Accounting for the latitudinal variations of the source appreciably improves simulations.

Effective Vertical Diffusion by Atmospheric Gravity Waves

AbstractQuantification of heat and constituent transport by gravity waves (GWs) in global models is challenging due to limited model resolutions. Current parameterization schemes suffer from oversimplification and often underestimate the transport rate. In this study, a new approach is explored to quantify the effective vertical eddy diffusion by using a high‐resolution Whole Atmosphere Community Climate Model (WACCM) simulation based on scale invariance. The WACCM simulation can partially resolve the mesoscale GW spectrum down to 250 km horizontal wavelength, and the heat flux and the effective vertical eddy diffusion by these waves are calculated directly. The effective vertical diffusion by the smaller‐scale, unresolved waves, is then deduced based on scale invariance, following the method outlined by H.‐L. Liu (2019) in quantifying GW momentum flux and forcing. The effective vertical diffusion obtained is generally larger than that obtained from parameterizations, and is comparable with that derived from observations in the mesosphere and lower thermosphere region.

Probing the Analytical Cancellation Factor of Short Scale Gravity Waves Using Na Lidar and Nightglow Data from the Andes Lidar Observatory

The cancellation factor (CF) is a model for the ratio between gravity wave perturbations in the nightglow intensity to those in the ambient temperature. The CF model allows us to estimate the momentum and energy flux of gravity waves seen in nightglow images, as well as the divergence of these fluxes due to waves propagating through the mesosphere and lower thermosphere region, where the nightglow and the Na layers are located. This study uses a set of wind/temperature Na lidar data and zenith nightglow image observations of the OH and O(1S) emissions to test and validate the CF model from the experimental perspective. The dataset analyzed was obtained during campaigns carried out at the Andes Lidar Observatory (ALO), Chile, in 2015, 2016, and 2017. The modeled CF was compared with observed CF values calculated using the ratio of wave amplitude in nightglow images to that seen in lidar temperatures for vertically propagating waves. We show that, in general, the modeled CF underestimates the observed CF results. However, the O(1S) emission line has better agreement with respect to the modeled value due to its supposedly simpler nightglow photochemistry. In contrast, the observed CF for the OH emission deviates by a factor of two from the modeled CF asymptotic value.

Seasonal Variations of High-Frequency Gravity Wave Momentum Fluxes and Their Forcing toward Zonal Winds in the Mesosphere and Lower Thermosphere over Langfang, China (39.4° N, 116.7° E)

Meteor radar data collected over Langfang, China (39.4° N, 116.7° E) were used to estimate the momentum flux of short-period (less than 2 h) gravity waves (GWs) in the mesosphere and lower thermosphere (MLT), using the Hocking (2005) analysis technique. Seasonal variations in GW momentum flux exhibited annual oscillation (AO), semiannual oscillation (SAO), and quasi-4-month oscillation. Quantitative estimations of GW forcing toward the mean zonal flow were provided using the determined GW momentum flux. The mean flow acceleration estimated from the divergence of this flux was compared with the observed acceleration of zonal winds displaying SAO and quasi-4-month oscillations. These comparisons were used to analyze the contribution of zonal momentum fluxes of SAO and quasi-4-month oscillations to zonal winds. The estimated acceleration from high-frequency GWs was in the same direction as the observed acceleration of zonal winds for quasi-4-month oscillation winds, with GWs contributing more than 69%. The estimated acceleration due to Coriolis forces to the zonal wind was studied; the findings were opposite to the estimated acceleration of high-frequency GWs for quasi-4-month oscillation winds. The significance of this study lies in estimating and quantifying the contribution of the GW momentum fluxes to zonal winds with quasi-4-month periods over mid-latitude regions for the first time.

An 18‐Year Climatology of Directional Stratospheric Gravity Wave Momentum Flux From 3‐D Satellite Observations

AbstractAtmospheric gravity waves (GWs) are key drivers of the atmospheric circulation, but their representation in general circulation models (GCMs) is challenging, leading to significant biases in middle atmospheric circulations. Unresolved GW momentum transport in GCMs must be parameterized, but global directional GW observations are needed to constrain this. Here we present an 18‐year climatology of directional stratospheric GW momentum flux (GWMF) from global AIRS/Aqua 3‐D satellite observations during 2002 to 2019. Striking hemispheric asymmetries are found at high latitudes, including dramatic reductions and reversals of GWMF during sudden stratospheric warmings. During Southern Hemisphere winter, a lateral convergence of GWMF toward 60°S is found that has no Northern Hemisphere counterpart. In the tropics, we find that zonal GWMF in AIRS measurements is strongly modulated by the semiannual oscillation (SAO) but not the quasi‐biennial oscillation (QBO). Our results provide guidance for future GW parameterizations needed to resolve long‐standing biases in GCMs.

Application of Deep Learning to Estimate Atmospheric Gravity Wave Parameters in Reanalysis Data Sets

AbstractGravity waves play an essential role in driving and maintaining global circulation. To understand their contribution in the atmosphere, the accurate reproduction of their distribution is important. Thus, a deep learning approach for the estimation of gravity wave momentum fluxes was proposed, and its performance at 100 hPa was tested using data from low‐resolution zonal and meridional winds, temperature, and specific humidity at 300, 700, and 850 hPa in the Hokkaido region (Japan). To this end, a deep convolutional neural network was trained on 29‐year reanalysis data sets (JRA‐55 and DSJRA‐55), and the final 5‐year data were reserved for evaluation. The results showed that the fine‐scale momentum flux distribution of the gravity waves could be estimated at a reasonable computational cost. Particularly, in winter, when gravity waves are stronger, the median root means square errors (RMSEs) of the maximum momentum flux and the characteristic zonal wavenumber were 0.06–0.13 mPa and 1.0 × 10−5, respectively.

Lidar observations of large-amplitude mountain waves in the stratosphere above Tierra del Fuego, Argentina

Large-amplitude internal gravity waves were observed using Rayleigh lidar temperature soundings above Rio Grande, Argentina (54^circ ; hbox {S}, 68^circ ; hbox {W}), in the period 16–23 June 2018. Temperature perturbations in the upper stratosphere amounted to 80 K peak-to-peak and potential energy densities exceeded 400 J/kg. The measured amplitudes and phase alignments agree well with operational analyses and short-term forecasts of the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF), implying that these quasi-steady gravity waves resulted from the airflow across the Andes. We estimate gravity wave momentum fluxes larger than 100 mPa applying independent methods to both lidar data and IFS model data. These mountain waves deposited momentum at the inner edge of the polar night jet and led to a long-lasting deceleration of the stratospheric flow. The accumulated mountain wave drag affected the stratospheric circulation several thousand kilometers downstream. In the 2018 austral winter, mountain wave events of this magnitude contributed more than 30% of the total potential energy density, signifying their importance by perturbing the stratospheric polar vortex.

Response of the Quasi‐Biennial Oscillation to a warming climate in global climate models

AbstractWe compare the response of the Quasi‐Biennial Oscillation (QBO) to a warming climate in eleven atmosphere general circulation models that performed time‐slice simulations for present‐day, doubled, and quadrupled CO2 climates. No consistency was found among the models for the QBO period response, with the period decreasing by 8 months in some models and lengthening by up to 13 months in others in the doubled CO2 simulations. In the quadrupled CO2 simulations, a reduction in QBO period of 14 months was found in some models, whereas in several others the tropical oscillation no longer resembled the present‐day QBO, although it could still be identified in the deseasonalized zonal mean zonal wind timeseries. In contrast, all the models projected a decrease in the QBO amplitude in a warmer climate with the largest relative decrease near 60 hPa. In simulations with doubled and quadrupled CO2, the multi‐model mean QBO amplitudes decreased by 36 and 51%, respectively. Across the models the differences in the QBO period response were most strongly related to how the gravity wave momentum flux entering the stratosphere and tropical vertical residual velocity responded to the increases in CO2 amounts. Likewise it was found that the robust decrease in QBO amplitudes was correlated across the models to changes in vertical residual velocity, parametrized gravity wave momentum fluxes, and to some degree the resolved upward wave flux. We argue that uncertainty in the representation of the parameterized gravity waves is the most likely cause of the spread among the eleven models in the QBO's response to climate change.

Sampling Errors in Observed Gravity Wave Momentum Fluxes from Vertical and Tilted Profiles

Observations from radiosondes or from vertically pointing remote sensing profilers are often used to estimate the vertical flux of momentum due to gravity waves. For planar, monochromatic waves, these vertically integrated fluxes are equal to the phase averaged flux and equivalent to the horizontal averaging used to deduce momentum flux from aircraft data or in numerical models. Using a simple analytical solution for two-dimensional hydrostatic gravity waves over an isolated ridge, it is shown that this equivalence does not hold for mountain waves. For a vertical profile, the vertically integrated flux estimate is proportional to the horizontally integrated flux and decays with increasing distance of the profile location from the mountain. For tilted profiles, such as those obtained from radiosonde ascents, there is a further sampling error that increases as the trajectory extends beyond the localised wave field. The same sampling issues are seen when the effects of the Coriolis force on the gravity waves are taken into account. The conclusion of this work is that caution must be taken when using radiosondes or other vertical profiles to deduce mountain wave momentum fluxes.

What can we learn about orographic drag parametrisation from high‐resolution models? A case study over the Rocky Mountains

AbstractComprehensive high‐resolution numerical weather prediction models provide a virtual laboratory for modelling the atmospheric flow over complex mountain ranges. In this study, global and regional simulations with horizontal grid spacing ranging from 2 to 32 km, focused over the northern Rocky Mountains, are used to assess the orographic blocking and gravity wave drag parametrisations employed in the Met Office Unified Model (UM) and the European Centre for Medium‐Range Weather Forecasts Integrated Forecasting System (IFS). The total, resolved and parametrised drag components in coarse‐resolution simulations are compared with those in high‐resolution simulations, in which the orographic drag processes are better resolved. The total surface stresses and gravity wave momentum fluxes in the free atmosphere of the global 16 km UM and IFS simulations are shown to compare well with 2 km regional simulations in terms of variability and mean. While the total gravity wave momentum flux is somewhat underestimated by both global models, its vertical distribution is well captured. The “seamlessness” of the parametrisation scheme is then assessed by comparing the total orographic stress – and its components – across several horizontal resolutions of the UM. The surface stress remains relatively constant across resolutions, such that the reduction in resolved orographic stress at coarser resolutions is compensated for by an almost equivalent increase in parametrised orographic stress. However, the parametrised orographic gravity wave momentum flux in the free atmosphere remains almost constant with resolution, failing to compensate for the lack of resolved flux at coarse resolutions. This leads to an underestimation of the total gravity wave drag at coarser resolutions. Further analysis suggests that this underestimation is due to the monochromatic wave assumption made by the gravity wave drag parametrisation scheme.