A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains

An analytical model of orographic gravity wave drag due to sheared flow past elliptical mountains is developed. The model extends the domain of applicability of the well-known Phillips model to wind profiles that vary relatively slowly in the vertical, so that they may be treated using a WKB approximation. The model illustrates how linear processes associated with wind profile shear and curvature affect the drag force exerted by the airflow on mountains, and how it is crucial to extend the WKB approximation to second order in the small perturbation parameter for these effects to be taken into account. For the simplest wind profiles, the normalized drag depends only on the Richardson number, Ri, of the flow at the surface and on the aspect ratio, γ, of the mountain. For a linear wind profile, the drag decreases as Ri decreases, and this variation is faster when the wind is across the mountain than when it is along the mountain. For a wind that rotates with height maintaining its magnitude, the drag generally increases as Ri decreases, by an amount depending on γ and on the incidence angle. The results from WKB theory are compared with exact linear results and also with results from a non-hydrostatic nonlinear numerical model, showing in general encouraging agreement, down to values of Ri of order one.

An analytical model of mountain wave drag for wind profiles with shear and curvature

An analytical model is developed to predict the surface drag exerted by internal gravity waves on an isolated axisymmetric mountain over which there is a stratified flow with a velocity profile that varies relatively slowly with height. The model is linear with respect to the perturbations induced by the mountain, and solves the Taylor–Goldstein equation with variable coefficients using a Wentzel–Kramers–Brillouin (WKB) approximation, formally valid for high Richardson numbers, Ri. The WKB solution is extended to a higher order than in previous studies, enabling a rigorous treatment of the effects of shear and curvature of the wind profile on the surface drag. In the hydrostatic approximation, closed formulas for the drag are derived for generic wind profiles, where the relative magnitude of the corrections to the leading-order drag (valid for a constant wind profile) does not depend on the detailed shape of the orography. The drag is found to vary proportionally to Ri21, decreasing as Ri decreases for a wind that varies linearly with height, and increasing as Ri decreases for a wind that rotates with height maintaining its magnitude. In these two cases the surface drag is predicted to be aligned with the surface wind. When one of the wind components varies linearly with height and the other is constant, the surface drag is misaligned with the surface wind, especially for relatively small Ri. All these results are shown to be in fairly good agreement with numerical simulations of mesoscale nonhydrostatic models, for high and even moderate values of Ri.

Stochastic resonance in a nonlinear model of a rotating, stratified shear flow, with a simple stochastic inertia-gravity wave parameterization

We report on a numerical study of the impact of short, fast inertia-gravity waves on the large-scale, slowly-evolving flow with which they co-exist. A nonlinear quasi-geostrophic numerical model of a stratified shear flow is used to simulate, at reasonably high resolution, the evolution of a large-scale mode which grows due to baroclinic instability and equilibrates at finite amplitude. Ageostrophic inertia-gravity modes are filtered out of the model by construction, but their effects on the balanced flow are incorporated using a simple stochastic parameterization of the potential vorticity anomalies which they induce. The model simulates a rotating, two-layer annulus laboratory experiment, in which we recently observed systematic inertia-gravity wave generation by an evolving, large-scale flow. We find that the impact of the small-amplitude stochastic contribution to the potential vorticity tendency, on the model balanced flow, is generally small, as expected. In certain circumstances, however, the parameterized fast waves can exert a dominant influence. In a flow which is baroclinically-unstable to a range of zonal wavenumbers, and in which there is a close match between the growth rates of the multiple modes, the stochastic waves can strongly affect wavenumber selection. This is illustrated by a flow in which the parameterized fast modes dramatically re-partition the probability-density function for equilibrated large-scale zonal wavenumber. In a second case study, the stochastic perturbations are shown to force spontaneous wavenumber transitions in the large-scale flow, which do not occur in their absence. These phenomena are due to a stochastic resonance effect. They add to the evidence that deterministic parameterizations in general circulation models, of subgrid-scale processes such as gravity wave drag, cannot always adequately capture the full details of the nonlinear interaction.

The climatology of the middle atmosphere in a vertically extended version of the Met Office's climate model. Part I: mean state

The climatology of a stratosphere-resolving version of the Met Office’s climate model is studied and validated against ECMWF reanalysis data. Ensemble integrations are carried out at two different horizontal resolutions. Along with a realistic climatology and annual cycle in zonal mean zonal wind and temperature, several physical effects are noted in the model. The time of final warming of the winter polar vortex is found to descend monotonically in the Southern Hemisphere, as would be expected for purely radiative forcing. In the Northern Hemisphere, however, the time of final warming is driven largely by dynamical effects in the lower stratosphere and radiative effects in the upper stratosphere, leading to the earliest transition to westward winds being seen in the midstratosphere. A realistic annual cycle in stratospheric water vapor concentrations—the tropical “tape recorder”—is captured. Tropical variability in the zonal mean zonal wind is found to be in better agreement with the reanalysis for the model run at higher horizontal resolution because the simulated quasi-biennial oscillation has a more realistic amplitude. Unexpectedly, variability in the extratropics becomes less realistic under increased resolution because of reduced resolved wave drag and increased orographic gravity wave drag. Overall, the differences in climatology between the simulations at high and moderate horizontal resolution are found to be small.

Inter-hemispheric mesospheric coupling in a comprehensive middle atmosphere model

Observations of noctilucent clouds have revealed a surprising coupling between the winter stratosphere and the summer polar mesopause region. In spite of the great distance involved, this inter-hemispheric link has been suggested to be the principal reason for both the year-to-year variability and the hemispheric differences in the frequency of occurrence of these high-altitude clouds. In this study, we investigate the dynamical influence of the winter stratosphere on the summer mesosphere using simulations from the vertically extended version of the Canadian Middle Atmosphere Model (CMAM). We find that for both Northern and Southern Hemispheres, variability in the summer polar mesopause region from one year to another can be traced back to the planetary-wave flux entering the winter stratosphere. The teleconnection pattern is the same for both positive and negative wave-flux anomalies. Using a composite analysis to isolate the events, it is argued that the mechanism for interhemispheric coupling is a feedback between summer mesosphere gravity-wave drag (GWD) and zonal wind, which is induced by an anomaly in mesospheric cross-equatorial flow, the latter arising from the anomaly in winter hemisphere GWD induced by the anomaly in stratospheric conditions.

Vertical propagation of information in a middle atmosphere data assimilation system by gravity-wave drag feedbacks

The mesospheric response to the 2002 Antarctic Stratospheric Sudden Warming (SSW) is analysed using the Canadian Middle Atmosphere Model Data Assimilation System (CMAM-DAS), where it represents a vertical propagation of information from the observations into the data-free mesosphere. The CMAM-DAS simulates a cooling in the lowest part of the mesosphere which is accomplished by resolved motions, but which is extended to the mid- to upper mesosphere by the response of the model's non-orographic gravity-wave drag parameterization to the change in zonal winds. The basic mechanism is that elucidated by Holton consisting of a net eastward wave-drag anomaly in the mesosphere during the SSW, although in this case there is a net upwelling in the polar mesosphere. Since the zonal-mean mesospheric response is shown to be predictable, this demonstrates that variations in the mesospheric state can be slaved to the lower atmosphere through gravity-wave drag.

Resonant gravity-wave drag enhancement in linear stratified flow over mountains

High-drag states produced in stratified flow over a 2D ridge and an axisymmetric mountain are investigated using a linear, hydrostatic, analytical model. A wind profile is assumed where the background velocity is constant up to a height z1 and then decreases linearly, and the internal gravity-wave solutions are calculated exactly. In flow over a 2D ridge, the normalized surface drag is given by a closed-form analytical expression, while in flow over an axisymmetric mountain it is given by an expression involving a simple 1D integral. The drag is found to depend on two dimensionless parameters: a dimensionless height formed with z_1, and the Richardson number, Ri, in the shear layer. The drag oscillates as z_1 increases, with a period of half the hydrostatic vertical wavelength of the gravity waves. The amplitude of this modulation increases as Ri decreases. This behaviour is due to wave reflection at z_1. Drag maxima correspond to constructive interference of the upward- and downward-propagating waves in the region z < z_1, while drag minima correspond to destructive interference. The reflection coefficient at the interface z = z_1 increases as Ri decreases. The critical level, z_c, plays no role in the drag amplification. A preliminary numerical treatment of nonlinear effects is presented, where z_c appears to become more relevant, and flow over a 2D ridge qualitatively changes its character. But these effects, and their connection with linear theory, still need to be better understood.

The effect of wind shear and curvature on the gravity wave drag produced by a ridge

The analytical model proposed by Teixeira, Miranda, and Valente is modified to calculate the gravity wave drag exerted by a stratified flow over a 2D mountain ridge. The drag is found to be more strongly affected by the vertical variation of the background velocity than for an axisymmetric mountain. In the hydrostatic approximation, the corrections to the drag due to this effect do not depend on the detailed shape of the ridge as long as this is exactly 2D. Besides the drag, all the perturbed quantities of the flow at the surface, including the pressure, may be calculated analytically.

Is missing orographic gravity wave drag near 60°S the cause of the stratospheric zonal wind biases in chemistry–climate models?

Nearly all chemistry–climate models (CCMs) have a systematic bias of a delayed springtime breakdown of the Southern Hemisphere (SH) stratospheric polar vortex, implying insufficient stratospheric wave drag. In this study the Canadian Middle Atmosphere Model (CMAM) and the CMAM Data Assimilation System (CMAM-DAS) are used to investigate the cause of this bias. Zonal wind analysis increments from CMAMDAS reveal systematic negative values in the stratosphere near 608S in winter and early spring. These are interpreted as indicating a bias in the model physics, namely, missing gravity wave drag (GWD). The negative analysis increments remain at a nearly constant height during winter and descend as the vortex weakens, much like orographic GWD. This region is also where current orographic GWD parameterizations have a gap in wave drag, which is suggested to be unrealistic because of missing effects in those parameterizations. These findings motivate a pair of free-runningCMAMsimulations to assess the impact of extra orographicGWDat 608S. The control simulation exhibits the cold-pole bias and delayed vortex breakdown seen in the CCMs. In the simulation with extra GWD, the cold-pole bias is significantly reduced and the vortex breaks down earlier. Changes in resolved wave drag in the stratosphere also occur in response to the extra GWD, which reduce stratospheric SH polar-cap temperature biases in late spring and early summer. Reducing the dynamical biases, however, results in degraded Antarctic column ozone. This suggests that CCMs that obtain realistic column ozone in the presence of an overly strong and persistent vortex may be doing so through compensating errors.

Influence of the Antarctic ozone hole on the polar mesopause region as simulated by the Canadian Middle Atmosphere Model

It is well established that variations in polar stratospheric winds can affect mesospheric temperatures through changes in the filtering of gravity wave fluxes, which drive a residual circulation in the mesosphere. The Canadian Middle Atmosphere Model(CMAM) is used to examine this vertical coupling mechanism in the context of the mesospheric response to the Antarctic ozone hole. It is found that the response differs significantly between late spring and early summer, because of a changing balance between the competing effects of parametrised gravity wavedrag (GWD)and changes in resolved wave drag local to the mesosphere. In late spring, the strengthened stratospheric westerlies arising from the ozone hole lead to reduced eastward GWD in the mesosphere and a warming of the polar mesosphere, just as in the well known mesospheric response to sudden stratospheric warmings, but with an opposite sign.In early summer, with easterly flow revailing over most of the polar stratosphere,the strengthened easterly wind shear within the mesosphere arising from the west ward GWD anomaly induces a positive resolved wave drag anomaly through baroclinic instability. The polar cooling induced by this process completely dominates the upper mesospheric response to the ozone hole in early summer. Consequences for the past and future evolution of noctilucent clouds are discussed

Chemistry–climate model simulations of twenty-first century stratospheric climate and circulation changes

The response of stratospheric climate and circulation to increasing amounts of greenhouse gases (GHGs) and ozone recovery in the twenty-first century is analyzed in simulations of 11 chemistry–climate models using near-identical forcings and experimental setup. In addition to an overall global cooling of the stratosphere in the simulations (0.59 6 0.07 K decade21 at 10 hPa), ozone recovery causes a warming of the Southern Hemisphere polar lower stratosphere in summer with enhanced cooling above. The rate of warming correlates with the rate of ozone recovery projected by the models and, on average, changes from 0.8 to 0.48 Kdecade21 at 100 hPa as the rate of recovery declines from the first to the second half of the century. In the winter northern polar lower stratosphere the increased radiative cooling from the growing abundance of GHGs is, in most models, balanced by adiabatic warming from stronger polar downwelling. In the Antarctic lower stratosphere the models simulate an increase in low temperature extremes required for polar stratospheric cloud (PSC) formation, but the positive trend is decreasing over the twenty-first century in all models. In the Arctic, none of the models simulates a statistically significant increase in Arctic PSCs throughout the twenty-first century. The subtropical jets accelerate in response to climate change and the ozone recovery produces awestward acceleration of the lower-stratosphericwind over theAntarctic during summer, though this response is sensitive to the rate of recovery projected by the models. There is a strengthening of the Brewer–Dobson circulation throughout the depth of the stratosphere, which reduces the mean age of air nearly everywhere at a rate of about 0.05 yr decade21 in those models with this diagnostic. On average, the annual mean tropical upwelling in the lower stratosphere (;70 hPa) increases by almost 2% decade21, with 59% of this trend forced by the parameterized orographic gravity wave drag in the models. This is a consequence of the eastward acceleration of the subtropical jets, which increases the upward flux of (parameterized) momentum reaching the lower stratosphere in these latitudes.

Climatology and Forcing of the Quasi-Biennial Oscillation in the MAECHAM5 Model

The quasi-biennial oscillation (QBO) in the equatorial zonal wind is an outstanding phenomenon of the atmosphere. The QBO is driven by a broad spectrum of waves excited in the tropical troposphere and modulates transport and mixing of chemical compounds in the whole middle atmosphere. Therefore, the simulation of the QBO in general circulation models and chemistry climate models is an important issue. Here, aspects of the climatology and forcing of a spontaneously occurring QBO in a middle-atmosphere model are evaluated, and its influence on the climate and variability of the tropical middle atmosphere is investigated. Westerly and easterly phases are considered separately, and 40-yr ECMWF Re-Analysis (ERA-40) data are used as a reference where appropriate. It is found that the simulated QBO is realistic in many details. Resolved large-scale waves are particularly important for the westerly phase, while parameterized gravity wave drag is more important for the easterly phase. Advective zonal wind tendencies are important for asymmetries between westerly and easterly phases, as found for the suppression of the easterly phase downward propagation. The simulation of the QBO improves the tropical upwelling and the atmospheric tape recorder compared to a model without a QBO. The semiannual oscillation is simulated realistically only if the QBO is represented. In sensitivity tests, it is found that the simulated QBO is strongly sensitive to changes in the gravity wave sources. The sensitivity to the tested range of horizontal resolutions is small. The stratospheric vertical resolution must be better than 1 km to simulate a realistic QBO.

Sensitivity of simulated climate to conservation of momentum in gravity wave drag parameterization

The Canadian Middle Atmosphere Model is used to examine the sensitivity of simulated climate to conservation of momentum in gravity wave drag parameterization. Momentum conservation requires that the parameterized gravity wave momentum flux at the top of the model be zero and corresponds to the physical boundary condition of no momentum flux at the top of the atmosphere. Allowing momentum flux to escape the model domain violates momentum conservation. Here the impact of momentum conservation in two sets of model simulations is investigated. In the first set, the simulation of present-day climate for two model-lid height configurations, 0.001 and 10 hPa, which are identical below 10 hPa, is considered. The impact of momentum conservation on the climate with the model lid at 0.001 hPa is minimal, which is expected because of the small amount of gravity wave momentum flux reaching 0.001 hPa. When the lid is lowered to 10 hPa and momentum is conserved, there is only a modest impact on the climate in the Northern Hemisphere; however, the Southern Hemisphere climate is more adversely affected by the deflection of resolved waves near the model lid. When momentum is not conserved in the 10-hPa model the climate is further degraded in both hemispheres, particularly in winter at high latitudes, and the impact of momentum conservation extends all the way to the surface. In the second set of simulations, the impact of momentum conservation and model-lid height on the modeled response to ozone depletion in the Southern Hemisphere is considered, and it is found that the response can display significant sensitivity to both factors. In particular, both the lower-stratospheric polar temperature and surface responses are significantly altered when the lid is lowered, with the effect being most severe when momentum is not conserved. The implications with regard to the current round of Intergovernmental Panel on Climate Change model projections are discussed.

Angular momentum conservation and gravity wave drag parameterization: implications for climate models

The robustness of the parameterized gravity wave response to an imposed radiative perturbation in the middle atmosphere is examined. When momentum is conserved and for reasonable gravity wave drag parameters, the response to a polar cooling induces polar downwelling above the region of the imposed cooling, with consequent adiabatic warming. This response is robust to changes in the gravity wave source spectrum, background flow, gravity wave breaking criterion, and model lid height. When momentum is not conserved, either in the formulation or in the implementation of the gravity wave drag parameterization, the response becomes sensitive to the above-mentioned factors—in particular to the model lid height. The spurious response resulting from nonconservation is found to be nonnegligible in terms of the total gravity wave drag–induced downwelling.

Large-scale dynamics of the mesosphere and lower thermosphere: an analysis using the extended Canadian Middle Atmosphere Model

The extended Canadian Middle Atmosphere Model is used to investigate the large-scale dynamics of the mesosphere and lower thermosphere (MLT). It is shown that the 4-day wave is substantially amplified in southern polar winter in the presence of instabilities arising from strong vertical shears in the MLT zonal mean zonal winds brought about by parameterized nonorographic gravity wave drag. A weaker 4-day wave in northern polar winter is attributed to the weaker wind shears that result from weaker parameterized wave drag. The 2-day wave also exhibits a strong dependence on zonal wind shears, in agreement with previous modeling studies. In the equatorial upper mesosphere, the migrating diurnal tide provides most of the resolved westward wave forcing, which varies semiannually in conjunction with the tide itself; resolved forcing by eastward traveling disturbances is dominated by smaller scales. Nonmigrating tides and other planetary-scale waves play only a minor role in the zonal mean zonal momentum budget in the tropics at these heights. Resolved waves are shown to play a significant role in the zonal mean meridional momentum budget in the MLT, impacting significantly on gradient wind balance. Balance fails at low latitudes as a result of a strong Reynolds stress associated with the migrating diurnal tide, an effect which is most pronounced at equinox when the tide is strongest. Resolved and parameterized waves account for most of the imbalance at higher latitudes in summer. This results in the gradient wind underestimating the actual eastward wind reversal by up to 40%.