The physics of orographic gravity wave drag
- Department of Meteorology, University of Reading, Reading, UK
The drag and momentum fluxes produced by gravity waves generated in flow over orography are reviewed, focusing on adiabatic conditions without phase transitions or radiation effects, and steady mean incoming flow. The orographic gravity wave drag is first introduced in its simplest possible form, for inviscid, linearized, non-rotating flow with the Boussinesq and hydrostatic approximations, and constant wind and static stability. Subsequently, the contributions made by previous authors (primarily using theory and numerical simulations) to elucidate how the drag is affected by additional physical processes are surveyed. These include the effect of orography anisotropy, vertical wind shear, total and partial critical levels, vertical wave reflection and resonance, non-hydrostatic effects and trapped lee waves, rotation and nonlinearity. Frictional and boundary layer effects are also briefly mentioned. A better understanding of all of these aspects is important for guiding the improvement of drag parametrization schemes.
The atmosphere above the boundary layer is almost invariably characterized by positive static stability, supporting the propagation of internal gravity waves, where buoyancy is the main restoring force. Related waves exist over a wide range of horizontal scales, from ≈1 km, where the waves are non-hydrostatic, to ≈10 km, where they are approximately hydrostatic, through to ≈100 km, where rotation of the Earth becomes important (inertia-gravity waves), up to ≈1000 km, where the variation of the Coriolis parameter with latitude must be taken into account (Rossby-gravity waves) . One of the key permanent sources for these waves is the Earth's orography, although other transient sources exist, for example, thermal contrasts between the continents and oceans, or atmospheric convection .
If these orographically generated waves (hereafter called mountain waves, for simplicity, although the term is often reserved for mesoscale waves) are able to propagate energy away from their sources (i.e., they are not evanescent), they exert a drag force on the mountains that generate them. In turn, by Newton's 3rd law, those mountains must exert an equal and opposite force on the atmosphere. However, while drag on the orography must obviously be exerted at the surface, the reaction force on the atmosphere may be distributed in space, and often acts at high altitudes. For mountain waves whose energy propagates vertically, the levels where the reaction force is exerted are given by Eliassen-Palm's theorem , or its extension to 3D flow by Broad .
Since the horizontal scales of the shortest mountain waves (primarily those unaffected by rotation) are smaller than current grid spacings in global weather and climate prediction models, these waves cannot be represented explicitly and so must be parametrized, because they affect the resolved flow [5, 6]. Perceived overestimations in the intensity of westerly winds in mid-latitudes were initially alleviated by adopting a so-called “envelope orography” [7, 8], whereby the surface elevation in mountainous areas is artificially increased by a certain amount. But a mountain wave drag parametrization was recognized by Palmer et al.  as a more physically satisfactory way of correcting this bias. Since then, increasingly sophisticated drag parametrizations have been developed and incorporated in large-scale models [10–19].
Drag not only decelerates the global atmospheric circulation, but, via thermal wind balance, it also has an impact on temperatures in the stratosphere at high latitudes. The formation of polar stratospheric clouds, which have an important role in ozone depletion, is very sensitive to those temperatures , and therefore to the drag representation. Additionally, it has been shown recently that the Brewer-Dobson circulation in the stratosphere, which ascends at the tropics and descends at the poles, contributing to the transport of ozone and other chemical species, is partially driven by mountain waves, although its primary forcing is believed to be non-orographic gravity waves .
On the other hand, mountain wave drag integrated globally determines the torque exerted on the Earth by the atmosphere and vice-versa, which affects the angular momentum budget of the atmosphere [22, 23]. This budget may be evaluated independently using series of the duration of day [24, 25]. There are known imbalances in the global angular momentum budget as simulated by numerical models, and this may be caused by errors in drag parametrizations [26–28]. The total terrain-induced drag may be split into three components: 1. drag associated with mountain waves, 2. drag associated with low-level blocking caused by stratification effects, and 3. turbulent form drag. There is the possibility, not only that the total drag is incorrectly estimated, but also that errors in each of these 3 components partially compensate.
For all these reasons, a fundamental understanding of mountain wave drag is of the utmost importance for the reliability of the models used to forecast the weather, reproduce the present climate, and especially predict future climate scenarios, on which important political decisions are based.
Several field campaigns to investigate flow over mountains have taken place over the last decades, for example: ALPEX [29, 30], PYREX , MAP [32, 33] and T-REX [34, 35]. These campaigns were generally carried out by large consortia, an exception being the smaller-scale campaign in the island of Madeira reported by Miranda et al. . While they contributed decisively to advance our knowledge about orographic flows, and included pressure measurements that allowed obtaining rough estimates of mountain wave drag, their success in this respect was limited, because of the limited available data sampling and complexity of real-world conditions (see also ). This hampered a more detailed comparison with theory, numerical simulations and lab experiments, which tend to assume much more idealized conditions.
There have been numerous reviews about flow over mountains [38–45], some on mountain wave drag parametrization [46–48], and also a few on the dependence of the resolved drag on numerical model resolution [33, 49], but none focused specifically on how the drag is affected by different physical processes. The present contribution aims to fill that gap, trying to show how accumulated knowledge, primarily from theory and numerical simulations, helps us understand the behavior of mountain wave drag.
Despite having grown rapidly only over the last 30 years, the subject has become too vast to cover entirely in a single review paper. Therefore, I decided to omit here studies focused on the effects of moisture, phase transitions and precipitation, radiation, and mean flow unsteadiness. Numerical aspects of mountain wave simulations, such as resolution dependence, and aspects directly linked with parametrization, are also, to a large extent, overlooked. I apologize in advance for failing to acknowledge several excellent investigations, even in the research areas that are covered here, and also for my probable bias toward theory.
The remainder of this paper comprises: an introduction to mountain wave drag in its simplest form in section 2, followed by an overview of the impact on the drag of various physical effects. Section 3 deals with orography anisotropy, section 4 addresses vertical wind shear, and in sections 5 and 6 the effects of total and partial critical levels are reviewed. Section 7 deals with partial wave reflection and resonance in layered atmospheres, in section 8 non-hydrostatic effects and trapped lee waves are mentioned, followed by rotation in section 9, and specifically nonlinear effects in section 10. To conclude, section 11 presents a brief overview of boundary layer effects. The paper ends with a summary of some open questions and possible future directions in the field.
2. General Framework for Linearized Mountain Waves
To understand how different physical processes affect mountain wave drag it is useful to formulate first the mountain wave problem in its most basic form. Consider the inviscid, adiabatic, non-rotating, steady, linearized equations of motion with the Boussinesq approximation:
Here (U, V) is the mean incoming wind velocity (which is assumed to depend only on height z), (u, v, w) is the velocity perturbation associated with the waves, p is the pressure perturbation, ρ0 is a reference constant density, b0 = gθ/θ0 is the buoyancy (where g is the acceleration of gravity, θ is a potential temperature perturbation and θ0 is a reference constant potential temperature), and N2 = (g/θ0)dΘ/dz is the static stability (where Θ is the mean potential temperature, assumed to depend only on z). These five equations for the five unknowns u, v, w, p and b0 may be combined into a single equation for w,
which must be solved subject to appropriate boundary conditions.
If the waves are forced by an isolated mountain, w (and all other perturbation variables) must decay to zero as x → ±∞ and y → ±∞. Additionally, the flow must follow the terrain elevation at the surface, i.e., in the linearized approximation
where h(x, y) is the surface elevation and (U0, V0) = (U, V)(z=0) is the mean incoming wind velocity at the surface. The waves must also either decay to zero as z → +∞, or their energy must propagate upward (the so-called radiation boundary condition).
Differentiating (1) with respect to x and (2) with respect to y, adding the two equations and using also (5) to eliminate u and v, an equation for p expressed in terms of w is obtained:
In the linearized approximation, the pressure drag force exerted on the mountain is given in general by:
where Dx and Dy are the drag components along the x and y direction, respectively. An important related quantity is the vertical flux of horizontal momentum associated with the waves,
It can be shown using (1, 2, 7), integrating by parts and using also the horizontal boundary conditions at x, y → ±∞ that:
which basically expresses Newton's 3rd law.
If the mountain is isolated and the flow perturbations decay to zero horizontally away from it, h, w and p may be expressed as Fourier integrals, viz.
where ĥ, ŵ and are the corresponding Fourier transforms, , and (k, l) is the horizontal wavenumber vector. When (13) is inserted into (6), that equation becomes:
which is sometimes called the Taylor-Goldstein equation. The lower boundary condition (7) can then be written as:
From (8, 13, 14), it can be shown that:
On the other hand, using (12) and (14), it is possible to express (9) in wavenumber space,
and also (10) as:
where the asterisk denotes complex conjugate. Treating the wave problem using Fourier analysis is particularly convenient, because the criteria determining whether the waves are vertically propagating of evanescent must be specified in wavenumber space.
2.1. Drag for Hydrostatic Flow with Constant Wind and Static Stability
If the flow is hydrostatic and both N and (U, V) = (U0, V0) are constant, all terms within square brackets in (15) vanish except the first one, which becomes constant. The analytical solution to that equation is then
is the vertical wavenumber of the waves and it has been assumed that N2 > 0, hence N > 0 is real. Note that, for hydrostatic flow m is always real, which means that all waves propagate vertically. Equation (21) already implies that the wave energy propagates upward, since the condition for the vertical component of the group velocity to be positive is that m and Uk + Vl have the same sign .
If the solution (20) is inserted into (17) and that equation is used in (18), then
where (21) has also been employed. If, additionally, the orography h(x, y) is assumed to have a height of O(h0) and a width of O(a) (see Figure 1), (k, l) may be made dimensionless using a as (k', l') = (ak, al), and ĥ using both a and h0 as ĥ′ = ĥ/(h0 a2), in terms of which (22) becomes:
where = (U2 + V2)1/2 is the magnitude of the incoming wind velocity and ψ is the angle it makes with the x direction, defined by:
Figure 1. Schematic diagram of mountain waves for a flow with ψ = 0. a and b are the mountain half-widths along x and y (for an axisymmetric mountain b = a) and h0 is the mountain height. Marked in the figure are flow streamlines (thick black lines), the boundary layer top (dashed blue line) and rotors (curved green arrows, see section 8). Note that throughout most of this paper, except in section 11, the boundary layer is neglected and a free-slip boundary condition is assumed at the surface.
In (23) the integral is dimensionless, and its numerical value depends only on the wind direction and shape of the orography. The drag therefore scales as ρ0Nah20 as is well known . In the case of a circular bell-shaped mountain,
the value of the drag in any direction is the same, and from (23), taking ψ = 0 without any loss of generality, the value of the integral is evaluated to be 1/(16π), so the drag is given by:
It should be noted that, for vertically propagating waves, the horizontal wavelength and amplitude of the dominant waves can be estimated from the corresponding orography width a and height h0, respectively. This implies that ã = Na/U quantifies non-hydrostatic effects (the flow is highly non-hydrostatic when ã ≤ O(1) and is hydrostatic when ã ≫ 1), while = Nh0/U quantifies flow nonlinearity (conditions are approximately linear for ≪ 1 and nonlinear when is of O(1) or larger). The effects of rotation can also be estimated in a similar way: if fa/U ≪ 1 (where f is the Coriolis parameter) then rotation is negligible, while for fa/U of O(1) or larger, it substantially affects the waves. In the present section it was assumed from the outset that ≪ 1 and fa/U ≪ 1, and in this subsection in particular also that ã ≫ 1. These assumptions will be relaxed in the following sections.
3. Orography Anisotropy
In general, the orography that generates mountain waves is not axisymmetric. A useful limit, representative of elongated mountain ranges, is that of 2D orography (say, aligned in the y direction), where (1–8), (15–17) and (20, 21) in the preceding section remain valid, provided that ∂n (…)/∂yn = 0 (for any n = 1, 2, 3, …) and l = 0 are imposed on them, since the flow is now independent of y. In that case, the Fourier transforms should be understood as being 1D (only taken along x) instead of 2D. The drag, on the other hand, must be defined per unit length along y, and is given by:
as an integral in physical or in wavenumber space, respectively. Following a procedure similar to that used in the preceding section, it can be shown that this drag scales as ρ0NUh20 and, for a bell-shaped mountain ridge,
(where now ĥ′ = ĥ/(h0a)), the drag takes the form
While the drag on an axisymmetric mountain is aligned with the incoming wind, the drag on a 2D ridge is always perpendicular to the ridge (and only depends on the wind in that direction, as shown by (30)). The considerably more realistic case of flow over a mountain with an elliptical horizontal cross-section, with a generic form h(x, y) = h[(x/a)2 + (y/b)2], was considered by Phillips , where b is the width of the mountain along y (see Figure 1). This kind of geometry is relevant, because it can be used to represent more accurately the anisotropy of real orography. This approach is adopted in each model grid box in the drag parametrizations of weather forecast models at the European Centre for Medium-Range Weather Forecasts (ECMWF)  and UK Meteorological Office . Using (23) it can be shown that, for this geometry, the two components of the drag take the form
γ = a/b is the aspect ratio of the mountain and ĥ is now normalized as ĥ′ = ĥ/(h0ab). G is a coefficient that describes the variation of the surface elevation in the radial direction.
In (31) the drag components are aligned with the main axes of the elliptical mountain. The components aligned with the mean incoming wind and perpendicular to it are 
Graphs of these two quantities normalized by ρ0Nbh20 (D3 and T3, respectively) are presented in Figure 2 as functions of γ and ψ. Clearly, the only situations with γ ≠ 1 where the drag is aligned with the incoming wind are when this wind is aligned with the principal axes of the mountain (ψ = 0o or ψ = 90o). Otherwise, the drag is oblique to both the incoming wind and the principal axes of the mountain.
Figure 2. Normalized drag components parallel (D3) and perpendicular (T3) to the incoming wind as a function of ψ and γ (see curve labels for details). Reproduced from Figure 8 of Phillips . © American Meteorological Society. Used with permission.
Hines  showed that the contribution to the drag given by wavenumbers within a small angular range surrounding the direction perpendicular to the mountain is overwhelmingly dominant for elongated mountains with aspect ratios γ < 1/3. For that reason, he suggested that the drag be approximated as perpendicular to the major axis of the mountain in those cases, and as receiving contributions from two wavenumber directions symmetrical with respect to the incoming wind otherwise. This allowed the derivation of analytical drag expressions (involving an explicit evaluation of approximated forms for (33,34)). These ideas were used by Scinocca and McFarlane  to develop a General Circulation Model (GCM) drag parametrization.
Bauer et al.  and Wells et al.  studied the interaction of orography anisotropy with flow nonlinearity also assuming elliptical mountains, concluding that flow perpendicular to the major axis of the mountain is substantially more nonlinear than flow parallel to that axis. Wells et al. , in particular, found that the dependence of the drag on flow direction is quite well captured by linear theory, but its value is typically underestimated for flow across the mountain (and overestimated for flow along it). Wells et al.  additionally explored the effect of weak rotation (fa/U ≈ 0.1), which was shown to introduce a lateral asymmetry in flow over symmetric mountains, but otherwise did not affect much the pressure distribution or the drag. Epifanio and Durran  further investigated the joint effects of orography anisotropy and nonlinearity, showing that in flow across an elongated 3D mountain the drag may only be approximated accurately by the drag in flow across an infinite 2D ridge for relatively weak nonlinearity, whereas finite anisotropy becomes much more important for strongly nonlinear flow.
4. Incoming Wind Shear
Blumen and McGregor 's pioneering study used linear theory to calculate the drag over a 2D and an approximately axisymmetric mountain, both for an incoming flow with horizontal shear and for a two-layer atmosphere with positive vertical wind shear and lower static stability in the bottom layer, and no shear and higher stability in the top layer. They found that the drag increased relative to its constant-wind and constant-stability reference value due to horizontal shear, was higher for flow over a 2D mountain than for flow over a 3D mountain, and was an oscillating function of the lower layer height. Blumen and McGregor  interpreted this as due to interference between waves transmitted and reflected at the interface between the two layers (see section 7 below), but did not systematically explore the behavior of the drag with the flow parameters.
An important measure of the wind shear in the incoming wind profile is the Richardson number
which quantifies the opposing effects of shear and stratification in destabilizing and stabilizing the flow, respectively. In the cases addressed next Ri > 1/4, which corresponds to dynamically stable (i.e., non-turbulent) flow.
Smith  developed a linear model for lee cyclogenesis, where he assumed a unidirectional wind profile with constant negative shear,
where U0 > 0 and α < 0 are constants, and as a by-product calculated mountain wave drag for hydrostatic, non-rotating flow over a 2D bell-shaped ridge (29), yielding
where, from (38), Ri = N2/α2, and the drag is normalized here by the corresponding unsheared reference value (30). Equation (39) shows that the drag decreases from its unsheared value as Ri decreases, until it vanishes at Ri = 1/4, where the flow is expected to become dynamically unstable. Grisogono  used linear theory conjugated with a WKB approximation  to calculate additional wind profile effects on the surface drag exerted by flow over a 2D ridge, finding that the drag increases for a wind profile with negative curvature at the surface, but decreases instead for positive curvature.
Grubišić and Smolarkiewicz  developed linear and numerical calculations of hydrostatic unidirectional negatively sheared flow (38) over an axisymmetric mountain that extended those of Smith  to 3D orography. They obtained a semi-analytical formula for the surface drag,
(normalized by (26)), which predicts that the drag is independent of the sign of the shear rate (a property that actually also applies to (39)) and that the drag at Ri = 1/4 is non-zero (reducing to Dx/D0x = 4/(3π)), in contrast with the 2D case. Figure 3 (left panel) shows the variation of the drag with Ri from (40) and from numerical simulations by Grubišić and Smolarkiewicz  and Teixeira et al. . The agreement is good down to quite low values of Ri. Grubišić and Smolarkiewicz  also found that nonlinear effects become important at progressively lower values of the dimensionless mountain height as Ri decreases.
Figure 3. Normalized drag as a function of Ri. Left: drag for wind profile (38) as given by (42) (solid line), from (40) (dotted line), and from numerical simulations (symbols). Right: drag for wind profile (41) as given by (43) (solid line), from Grisogono  (dashed line), from the corresponding high-Ri asymptotic limit (dotted line), and from numerical simulations (symbols). Reproduced from Figures 1, 7 of Teixeira et al. . © American Meteorological Society. Used with permission.
For the same kind of wind profile (38), but for flow over a 2D ridge, Wang and Lin  investigated the interaction between shear and nonlinear effects for various flow quantities, using high values of Ri. They found that the drag could exceed by a large amount the linear estimate (39), taking normalized values larger than 2 or 3. For = O(1) and Ri > 300, negative shear was found to enhance the drag more than positive shear, but for higher this relation was reversed. Curiously, Wang and Lin  found that, in contrast with Grubišić and Smolarkiewicz , nonlinear effects become important at lower values of as Ri increases. But the range of Ri and they considered was much higher: Ri > 20 and > 0.8, while in Grubišić and Smolarkiewicz  Ri ≤ 9 and ≤ 0.3, which perhaps explains the difference.
Teixeira et al.  developed a linear, hydrostatic model based on the WKB approximation to evaluate effects on the surface drag of generic, relatively slowly-varying wind profiles. They aimed to explain why in flow with constant shear over a 3D mountain the drag decreases with Ri (as in  or ), whereas for a directional shear flow where the wind turns with height at a constant rate keeping a constant magnitude, expressed as:
(where β is a constant), the drag increases instead as Ri decreases, as originally shown by Miranda  and confirmed by Valente . Teixeira et al.  found that only a 2nd-order WKB approximation is able to capture shear corrections to the surface pressure that are asymmetric with respect to the obstacle, and thus that modify the drag. Figure 3 shows the drag for hydrostatic flow over an axisymmetric mountain for the wind profiles (38) and (41), respectively, which is given by WKB theory  as:
respectively. Although in (41) the wind turns with height, and so does the drag in numerical simulations for sufficiently low Ri, this is not captured by (43). However, the behavior of the drag component aligned with the surface wind (Dx) is reproduced quite accurately. It is easy to show that (40) asymptotically approaches (42) as Ri → ∞.
Teixeira and Miranda  addressed the simpler case of flow over a 2D ridge, where for the linear wind profile (38) the drag was found to vary with Ri as:
Teixeira and Miranda  extended these calculations to flow over mountains with an elliptical horizontal cross section, such as assumed by Phillips , obtaining corrections to the drag due to wind profile effects that are in a form suitable to be implemented in the ECMWF drag parametrization . These corrections, like those in Teixeira et al.  and Teixeira and Miranda , do not depend on the detailed shape of the orography, a property that relies on the hydrostatic assumption. The impact of these corrections on the drag at the global scale and on the integrated torque exerted by mountains on the atmosphere was evaluated using reanalysis meteorological data and orography representative of the real mountain ranges around the world by Miranda et al. . They found that positive corrections to the drag may be important in localized regions such as Antarctica, perhaps contributing to a slight alleviation of the westerly bias exhibited by the drag torque.
More recently, Tang et al.  extended the model of Teixeira et al.  to a non-Boussinesq atmosphere. They showed using hydrostatic linear theory with a 2nd order WKB approximation that, for constant-shear flow over an axisymmetric mountain, the drag behavior depends on the sign of the shear. Figure 4 reproduces their results, for the linear wind profile with directional shear
showing that the drag, which in the Boussinesq approximation is given by:
(as shown originally by Teixeira et al. ), decreases more slowly with Ri when the shear is negative and faster when the shear is positive.
Figure 4. Normalized x (left) and y (right) components of the drag as a function of Ri for the wind profile (45). The solid and dash-dotted lines correspond to non-Boussinesq calculations (with different signs of α – see legend), and the dotted line is the original Boussinesq result (46). Reproduced from Figure 1 of Tang et al.  with kind permission from Springer Science and Business Media.
5. Total Critical Levels
A total critical level (also sometimes called 2D or unidirectional critical level) may be defined as a level in the atmosphere where the Taylor-Goldstein Equation (15) has a singularity for all wavenumbers in the mountain wave spectrum at the same time. This can occur when either (U, V) = 0 in flow over 3D orography, or when U = 0 in a flow which may have directional shear over 2D orography for which l = 0. Eliassen and Palm  have shown that for low-amplitude 2D mountain waves that propagate vertically, the vertical flux of horizontal momentum only varies at these critical levels, where the reaction force to the drag exerted on the mountain acts on the atmosphere. Eliassen-Palm's theorem does not, however, specify in which way this momentum flux varies, or the value of that force.
Booker and Bretherton  were the first to undertake that task. They found that, for linear flow, the momentum flux is always absorbed by critical levels, being attenuated by the factor
where Ric is Ri at the critical level zc. This means that, unless Ric is of O(1), there will be total absorption of the wave momentum flux, by a process which does not depend on viscosity and can be captured in an inviscid framework.
Booker and Bretherton's results were for low-amplitude waves. In the case of larger-amplitude, nonlinear mountain waves, Breeding  showed that total critical levels behave in a similar way to that predicted from linear theory (i.e., leading to nearly total momentum flux absorption) for Ric ≥ 5, but reflect a substantial amount of wave energy when 0.25 < Ri < 1. They estimated the amount of reflection as 35% for Ric = 0.53 and 7% for Ric = 2.12, but were unable to establish a relation between the reflection coefficient and the wave amplitude (i.e., nonlinearity).
As a mountain wave, or an internal gravity wave in general, propagates toward a critical level, the magnitude of the associated horizontal velocity fluctuations increases indefinitely in the inviscid linear approximation . This would necessarily lead to dynamical instability. Through linearized calculations including viscosity, Fritts and Geller  derived criteria for the viscous stabilization of critical levels, whereby this kind of instability is prevented. They concluded that the flow is stabilized if the viscous length scale associated with the critical levels is somewhat larger than the depth of the layer that becomes unstable in inviscid conditions. From this general criterion, they inferred that internal gravity waves should be stabilized by viscosity only at heights above 130 km. They pointed out, however, that if the atmosphere is turbulent, an eddy viscosity replaces the molecular viscosity in the definition of the viscous length scale, and so waves in the troposphere become typically stable at critical levels. Figure 5 shows the variation of the wave momentum flux at a critical level from viscous and inviscid theory.
Figure 5. Profile of the momentum flux in a region containing a total critical level (at ζ = 0, where ζ = (k|α|/ν)1/3 z is height scaled by the viscous length scale, with ν being viscosity) for high Ric. Solid line: Viscous theory, dashed line: inviscid theory. Reproduced from Figure 4 of Fritts and Geller . © American Meteorological Society. Used with permission.
Clark and Peltier  and Scinocca and Peltier  used the concept of wave reflection at a total critical level to explain the occurrence of downslope windstorms and high-drag states in 2D mountain waves. Scinocca and Peltier focused particularly on the effect of the Richardson number at the critical level Ric. (see section 7).
McFarlane  developed a drag parametrization based on results from linear theory, where the effect of wave saturation was taken into account in flows with directional shear. The distribution of stress with height was calculated by prescribing that the wave amplitude be limited due to wave breaking associated with the decay of density with height, or critical levels, which they defined as the levels where the wind becomes zero in the direction of the surface wind. This definition, which is independent of the orography orientation and thus not strictly consistent, nevertheless allows directional shear flows to be treated using the conceptual framework of total critical levels.
Grubišić and Smolarkiewicz  extended the work of Booker and Bretherton  by camputing the wave momentum flux absorption at a critical level for waves generated by a 3D mountain. For a mean incoming flow with the form (38), they obtained the following attenuation factor
where the difference from (47) results from the factor (k2 + l2)/k2. This implies that waves with a wavenumber vector perpendicular to the mean flow (k = 0 and l ≠ 0) are totally absorbed, no matter how low Ric is, whereas the attenuation of waves with wavenumber vectors parallel to the mean flow (l = 0) is similar to that of 2D waves.
Following Lindzen and Tung , but carrying out a more systematic exploration of parameter space, [76, 77] used linear theory to address wave ducting, and its implications for downslope windstorms, in flow over 2D mountains. They treated cases where both Ric > 0.25, and thus the critical level absorbs upward propagating waves in accordance with (47), and cases where Ric < 0.25, where the waves are unattenuated, but instead amplified by the critical level, in what has been termed “overreflection” .
Lott  addressed trapped lee wave absorption at the critical level that exists at the surface due to the no-slip boundary condition. Using a linear model including viscosity, he concluded that absorption of wave momentum at the surface is insignificant for Ric < 0.25 at high Reynolds number Re, leading to long-lived trapped lee waves, whereas it becomes progressively higher as Ric increases for Ric > 0.25. Lott  also found that wave absorption for Ric > 0.25 is always smaller than predicted by both inviscid critical-level theory and simplified models that represent friction as a Rayleigh damping (e.g., ).
6. Partial Critical Levels in Flow Over 3D Mountains
Partial critical levels (also called 3D or directional critical levels) are levels in flow with directional shear over 3D mountains where the Taylor-Goldstein Equation (15) becomes singular, but not for all wavenumbers in the wave spectrum simultaneously. Mathematically, they are defined by the condition
implying that the mean incoming wind (U, V) and the horizontal wavenumber of the waves (k, l) are perpendicular. Unlike total critical levels, which affect the whole wave spectrum at a single well-defined height, the height of partial critical levels depends on the azimuthal angle of the wavenumber. For this reason, these critical levels only exist for 3D orography, which contains wavenumbers with a range of orientations.
Using linear theory, Grimshaw  derived for the first time an expression for the momentum flux absorption at partial critical levels. The attenuation of the momentum flux in the absence of rotation for a linear wind profile expressed by:
where V0 and β are constants, was given by:
Equation (48) can be viewed as a particular case of (51) when the flow is unidirectional (V0 = β = 0).
Shutts  noted that most drag parametrizations assume that the wave momentum flux is associated with a single monochromatic wave, following . He tried to improve that situation by calculating the wave momentum flux in flows with directional shear over 3D mountains using hydrostatic linear theory. He developed a practical approach to compute the momentum flux over anisotropic orography, described in the spectral domain, including a selective critical-level absorption effect. Shutts  and Shutts and Gadian  supported this approach by calculating the momentum flux for idealized flow with constant directional shear over an axisymmetric mountain (a particular case of (50)) and also for a wind that turns with height, (41). They assumed that Ric ≫ 1, and thus that all momentum flux is totally absorbed at critical levels (see also ). Shutts and Gadian  also carried out numerical simulations that showed that the momentum flux (and the drag) can be substantially underestimated by linear theory, and can exit laterally from the domain at high levels due to spreading of the wave pattern away from its source (see Figure 6), even for hydrostatic flow. This has implications for drag parametrizations, which typically adopt a single-column approach.
Figure 6. Normalized momentum fluxes as a function of height for the wind profile (41) at very high Ri. Dashed line: infinite-area linear solution, dash- dotted line: linear solution taking into account the finite domain extent, solid line: numerical simulations. Reproduced from Figure 12 of Shutts and Gadian . Copyright © 1999 Royal Meteorological Society.
This implies that the momentum flux divergence d(Mx, My)/dz, which corresponds to the drag exerted on the atmosphere, is perpendicular to the wind at each height, and that the momentum flux (Mx, My) will not change unless the wind turns. This result places a strong constraint on the momentum flux profiles, which is of key importance for including directional shear effects in drag parametrizations. Broad  also considered the question of whether gravity waves in a directionally sheared flow break or not. This question was prompted by the fact that, before they reach a unidirectional critical level, finite-amplitude waves always overturn, since their amplitude tends to grow without bound before being dissipated. In directionally sheared flows, due to the nature of partial critical levels, this effect would be expected to be weaker. Using linear theory conjugated with ray tracing to clarify this issue, Broad  concluded that wave breaking always occurs. This depletes the momentum flux above wave breaking zones, qualitatively in the same way as critical level absorption in the linear analysis of Shutts and Gadian . However, Broad's flow overturning condition was formulated in Fourier space, which seems questionable.
Doyle and Jiang  reported observations of mountain waves over the south-western French Alps during the MAP campaign. These waves were of relatively small amplitude, and were embedded in a large-scale flow with substantial directional shear. They observed that the amplitude of the waves decreased rapidly with height (by about 50%) in a layer with directional shear where the wind turned by about 90o. This provided the first experimental support for the prediction by linear theory (e.g., [80, 81]) of momentum flux absorption by partial critical levels in flow with directional shear, highlighting the importance of representing this effect in drag parametrizations.
Using hydrostatic linear theory conjugated with a 3rd-order WKB approximation, Teixeira and Miranda  developed calculations of the momentum fluxes associated with mountain waves for generic (albeit slowly-varying) wind profiles over axisymmetric mountains, at values of Ri of O(1). Unlike Shutts  and Shutts and Gadian  (for Ric ≫ 1), Teixeira and Miranda  found that wave absorption at critical levels at low Ri is not total, but rather the attenuation factor is:
where the subscript c denotes values taken at the critical level. This further generalizes (51) to generic wind profiles with relatively high Ri (although actually it has good accuracy down to Ri ≥ 0.5). For the wind profile (50), (53) is the asymptotic limit of (51) for Ri → ∞. Figure 7 shows the normalized momentum flux for the wind profile (41), given by Teixeira and Miranda  as:
at a height where the wind has turned by an angle of 180o, and therefore the momentum flux has been filtered for the whole of the wave spectrum. Teixeira and Miranda  found that the momentum fluxes normalized by the surface drag in the absence of shear, like the surface drag, do not depend on the detailed shape of the orography, and may be expressed as 1D integrals, while the momentum flux divergence, which corresponds to the drag acting on the atmosphere, has a closed analytical form for generic wind profiles.
Figure 7. Normalized momentum flux at the height βz/π = 1 as a function of Ri for the wind profile (41). Solid line: (54), dashed line: Ri ≫ 1 limit (total absorption), symbols: numerical simulations for different values of : circles: 0.01, squares: 0.1, triangles: 0.2, diamonds: 0.5. Reproduced from Figure 5 Teixeira and Miranda of . © American Meteorological Society. Used with permission.
Xu et al.  used linear theory, conjugated with linear wind profiles and orography similar to those considered by Teixeira and Miranda  to show that the wave momentum flux turns with height in the same direction as the mean incoming wind, but by a smaller angle, and the vertical divergence of the momentum flux always vanishes at the surface, pointing perpendicularly to the left (right) of a mean flow that backs (veers) with height. This divergence was found to attain a maximum at a height proportional to the wind speed at the surface and inversely proportional to the shear intensity, of higher magnitude for large wind turning angles and low Ri. These findings are consistent with the results of Teixeira and Miranda , but had not been explicitly stated by them. Using nonlinear numerical simulations of the same flow for = 0.0125, 0.625, 1.25 and 2.5,  showed that, when = 0.625 the momentum fluxes are enhanced by nonlinear effects. For = 1.25, low- and high-drag states are produced for different wind-turning angles, depending on whether flow over the mountain or flow around the mountain (with lee vortices) are preferred.
Very recently, Teixeira and Yu  extended the model of Teixeira and Miranda  to mountains with an elliptical horizontal cross-section, bringing it one step closer to a concrete parametrization proposal. They also developed the necessary theory for the situation where the wind turns non-monotonically with height or by an angle larger than 180o, where there is more than one critical level per wavenumber of the wave spectrum. It is then necessary to split the atmosphere into various layers, and express either the momentum flux or its divergence as a function of the momentum flux at the bottom of each layer.
7. Partial Wave Reflection and Resonance
Since their first systematic recordings , downslope windstorms are known to be associated with high-drag states, where the drag is several times larger than expected from linear theory. High-drag states are thought to be caused by resonance associated with vertical wave reflection with constructive interference at sharp variations in atmospheric parameters.
Peltier and Clark  found that, for constant wind profiles, high-drag states can be explained by processes incorporated in Long's nonlinear wave model of flow over 2D mountains (see section 10). Other nonlinear effects, intrinsically related to the wind variation with height, could only be reproduced by numerical simulations (e.g., highly nonlinear trapped lee waves). Peltier and Clark  introduced the concept of wave amplification by reflection at self-induced (wave-overturning) critical levels, suggesting it as the key mechanism behind high-drag states. In order to investigate this issue further, Clark and Peltier  carried out numerical simulations with an incoming flow of the form
where d is the depth of the shear layer surrounding the critical level. In these simulations, where the Richardson number Ric = N2d2/U20 was kept above 0.25, resonant drag enhancement corresponding to constructive interference between upward and downward propagating waves below zc was observed for:
where n is an integer number and λz = 2πU0/N is the hydrostatic vertical wavelength of internal gravity waves in the constant-wind layer near the surface. Figure 8 displays the variation of the drag with zc/λz, clearly showing the spacing of λz between maxima, which is twice the value expected from linear resonance theory (see below).
Figure 8. Drag obtained from numerical simulations as a function of zc/λz (denoted here by zi/λz). Reproduced from Figure 5 of Clark and Peltier . © American Meteorological Society. Used with permission.
Smith  developed a different “hydraulic” approach, using a hydrostatic version of Long's equation to simulate flow over a 2D mountain for a constant upstream wind profile, with a “dead layer” downstream of the mountain, where the flow is turbulent, well-mixed (and therefore unstratified), and approximately quiescent on average. This model was able to predict the altitude of the well-mixed layer, the strength of the winds at the surface beneath that layer, and the mountain wave drag, which is given by:
where H0 denotes the height of the streamline upstream of the mountain that will divide to contain the well-mixed layer downstream of the mountain, and H1 is the depth of the accelerated layer beneath the well-mixed layer.
Bacmeister and Pierrehumbert  evaluated the merits of the theories of Clark and Peltier  and of Smith  through a number of numerical experiments of nonlinear flow over 2D mountains also with the wind profile (55), reaching the conclusion that Smith's hydraulic theory works better than the theory of Clark and Peltier. They criticized this latter theory for not justifying the phase of the reflected wave that leads to resonance, and found that the height of the critical level for which high-drag states occur depends on , in agreement with Smith's theory (in what was termed “resonance shift”). They additionally showed that high-drag states are possible for low values of , but take a long time to stabilize.
Although Smith's theory is believed to be the most satisfactory to date, other possible explanations for high-drag states have been explored. Durran , Wang and Lin , and Leutbecher  investigated the role played by static stability discontinuities in the establishment of high-drag states and downslope windstorms through wave ducting, and tested the accuracy of their representation by linear theory. Wang and Lin  for flow over 2D mountains and Leutbecher  for flow over an axisymmetric mountain showed that, when N is higher in the bottom layer and lower in the top layer, drag maxima (as well as wind intensity maxima) exist for:
where n is an integer and z1 is the height of the discontinuity in N, while for the opposite distribution of N maxima occur instead for z1 = (n/2) λz. Durran  showed that large drag amplification, and more severe underestimation of it by linear theory (by a factor of 3 or 4), occurs when the static stability is higher in the bottom layer than in the top one. He concluded that in such multi-layer configurations nonlinear effects may become relevant for ≥ 0.3, which is substantially lower than for an atmosphere with constant parameters. For = 0.67, 0.83, and 1, Wang and Lin  confirmed that the location of the drag maxima as a function of z1 is roughly unchanged, but their magnitude increases and the detailed variation of the drag with z1 is modified. In 3D flow, Leutbecher  found that the drag maxima become less pronounced, especially for high values of z1, because of directional wave dispersion.
Scinocca and Peltier  assessed the importance of the Richardson number at the critical level Ric in the establishment of high-drag states for the wind profile (55). They found that high-drag states take longer to achieve as Ric decreases, and may even be prevented for low Ric > 1/4, if is small enough. They also found that the steady-state high-drag states become insensitive to Ric for 0.625 ≤ ≤ 0.85, and attributed resonance shift to the different values of Ric used by previous authors, rather than to nonlinear effects.
Using a piecewise-linear wind profile qualitatively similar to (55), Wang and Lin [76, 77] studied the implications of wave ducting for downslope windstorms and high-drag states. They found that, at least in the linear regime, wave reflections at the shear discontinuities existing below and above the critical level (whose intensity is controlled by Ric) rather than at the critical level itself determine the wave response, with wind and drag maxima at the surface attained when (58) is fulfilled, where now z1 is the height of the lowest shear discontinuity. They explained some of the discrepancies observed between previous studies on high-drag states (namely  and ) by the use of zc instead of z1 as a key flow parameter. Wang and Lin  also showed how nonlinear high-drag states (with a spacing of λz rather than λz/2 as a function of z1) are selected from the linear ducted wave modes through modification of the static stability profile by the wave flow perturbations.
Miranda and Valente  treated the same kind of situation, but for flow over an axisymmetric mountain. concluding that the drag amplification is much more modest than for 2D mountains, and the spacing of high-drag states as a function of the critical level height is consistent with linear theory, rather than with Smith's theory, being (1/2) λz. The occurrence of resonance shift was also seen to be much reduced. As in Leutbecher , these effects are attributable to directional wave dispersion. Using linear theory and for the same wind profile over both 2D and axisymmetric mountains, Teixeira et al.  derived analytical expressions for the drag, the former of which can be written as:
with drag maxima occurring for (58), and being enhanced as Ric decreases. In weakly-nonlinear flow ( = 0.5), high drag-states were found to become much less sensitive to Ric (in agreement with ), and to occur for higher values of z1 (consistently with ) (see Figure 9). While in flow over an axisymmetric mountain, the locations of drag maxima were found to be accurately predicted by linear theory, in 2D flow intermediate maxima tended to disappear, confirming the results of Miranda and Valente . For the same 2D flow as considered by Teixeira et al. , Teixeira and Miranda  showed that wave breaking for Rec < 9/4 behaves differently from wave breaking in an atmosphere with constant parameters, occurring not directly over the mountain, but in regions displaced both vertically and horizontally upstream and downstream.
Figure 9. Drag normalized by (30) as a function of Nz1/(πU0). Left: linear flow. Lines: linear theory, symbols: numerical simulations. Solid line and squares: Ric = 2, dotted line and diamonds: Ric = 1, dashed line and circles: Ric = 0.5. Right: Numerical simulations for = 0.5. Squares: Ric = 2, circles: Ric = 1, triangles: Ric = 0.5. Reproduced from Figures 1, 5 of Teixeira et al. . Copyright © 2005 Royal Meteorological Society.
Teixeira et al.  showed that differences in the drag behavior for wind profiles where the wind varies linearly near the surface and is constant aloft depend (in hydrostatic conditions) on whether critical levels exist, and are located above or below the shear discontinuity. When a large fraction of the wave spectrum has a critical level below the shear discontinuity, the drag behaves relatively similarly as if the lower layer extended indefinitely. Otherwise, the drag undergoes oscillations due to resonant constructive and destructive interference of reflected waves.
Wells and Vosper  identified a wave and drag amplification mechanism by parametric resonance, when the profile of the Scorer parameter
varies sinusoidally with a wavelength of half the basic vertical hydrostatic wavelength. For flow over a 2D ridge, they showed that this could lead to an underestimation of the drag by linear theory by a factor of about 2.5 (see Figure 10). Teixeira et al.  developed a perturbation model to explain this drag amplification, showing that the drag behavior at resonance depends on the phase of the oscillation of the Scorer parameter, and also crucially on the representation of dissipative processes.
Figure 10. Drag normalized by its reference value for constant l0 as a function of λz/λs, where λs is the wavelength of the vertical oscillation of the Scorer parameter. See legend for details. Reproduced from Figure 9 of Wells and Vosper . © Crown Copyright 2010. Published by John Wiley & Sons, Ltd.
8. Non-Hydrostatic Effects and Trapped Lee Waves
Non-hydrostatic effects (which occur when ã ≲ O(1)) lead to mountain wave dispersion, which tends to weaken the drag in a flow with constant wind and static stability . However, when in a unidirectional flow the Scorer parameter decreases with height, non-hydrostatic effects promote the propagation of trapped lee waves near the surface, which may increase the drag.
Following Wurtele et al. , Keller  used linear theory to evaluate the accuracy of the hydrostatic approximation in flow over a 2D ridge with a linearly increasing incoming wind velocity (i.e., (38) with α > 0) either extending indefinitely, or in a layer with lower static stability, topped by a semi-infinite layer with no shear and higher static stability. He concluded that the momentum flux is substantially reduced by non-hydrostatic effects and its contributions are shifted downstream from the mountain due to wave trapping near the surface. The momentum flux was shown to depend on the lower layer height and Ri within it (as  would confirm later for hydrostatic flow). Sharman and Wurtele  treated essentially the same situation, but for flow with positive or negative shear extending indefinitely over 3D mountains of varying aspect ratio. They found that the drag over elongated mountains perpendicular to the flow does not change much for γ < 1/3. The drag also tended to become lower as the flow became more non-hydrostatic for negative shear while for positive shear and the same ã the drag was higher, a behavior they attributed to the existence of trapped lee waves.
Lott  calculated the drag associated with trapped lee waves using a 2D numerical linear model. Despite the intrinsic non-hydrostaticity of these waves, their associated drag was shown to be reasonably estimated using the hydrostatic reference value. However, Lott also found that the momentum flux profile was very different from that associated with vertically propagating waves, decaying with altitude in the absence of dissipation, because the assumptions behind Eliassen-Palm's theorem are not satisfied, and a substantial fraction of the wave momentum exits the computational domain through the leeward boundary. Broad  explained this variation with height of the momentum flux, showing that it is associated with a decrease of the dynamic pressure across the mountain because of the downstream wave train. He proposed an expression to replace the Eliassen-Palm theorem in those circumstances,
where W is the amplitude of the trapped lee waves far downstream of the obstacle.
In a pioneering study, Bretherton  systematically analyzed the drag produced by trapped lee waves, departing from Bernoulli's equation, but he assumed a rigid-lid upper boundary condition, which is not very realistic in the atmosphere. This limitation was removed by Smith  who generalized the trapped lee wave drag expression derived by Bretherton to an unbounded atmosphere, yielding
where the index of the sum j is over all trapped lee wave modes and kj (j = 1, 2,..) are the corresponding resonant wavenumbers. This expression is quite powerful, as it does not assume any specific form for the atmospheric profiles, other than that the integral in (62) must converge.
For practical applications, it is useful to explore the drag behavior for more specific atmospheric profiles. Tutiś  calculated trapped lee wave solutions in an atmosphere where the Scorer parameter l0 follows a hyperbolic-secant profile. She concluded that linear theory gives good predictions for the wavelength and drag associated with those waves, but the wave amplitude is underestimated. Unfortunately, Tutis used a rigid-lid upper boundary condition in her calculations, which is not very realistic. Adopting a two-layer model with higher Scorer parameter in a lower layer and lower Scorer parameter aloft (following ), Teixeira et al.  explicitly evaluated the drag associated with waves trapped in the lower layer and waves propagating in the upper layer, obtaining the following expression for the former drag component when the wind velocity is constant:
where m1 = (l21 − k2)1/2, n2 = (k2 − l22)1/2, l1 and l2 are the Scorer parameters in the lower and upper layers, respectively, and H is the height of the lower layer. DL normalized by (30), is plotted in Figure 11 (left). Teixeira et al.  considered a similar problem, but for the two-layer atmosphere of Vosper  instead, where the lower layer is neutral (l1 = 0), and there is a temperature inversion between the two layers, represented by a jump in potential temperature, which is quantified by a Froude number Fr. They found that in that case the trapped lee wave drag is due to interfacial waves propagating at the inversion, which have only one mode, being given by:
where kL is the resonant wavenumber of the trapped lee waves. Figure 11 (right) presents (64) normalized by (30) using the static stability of the upper layer.
Figure 11. Left: trapped lee wave drag (here denoted by D2) normalized by (30) as a function of l1H/π. Solid line: l1a = 10, dashed line: l1a = 5, dotted line: l1 a = 2. Reproduced from Figure 6 of Teixeira et al. . Copyright © 2012 Royal Meteorological Society. Right: Drag normalized by (30) as a function of Fr for l2H = 0.5 and l2a = 1. Solid line: total drag, dotted line: internal gravity wave drag, dashed line: trapped lee wave drag, all from theory; symbols: numerical simulations. Reproduced from Figure 9 of Teixeira et al. . © American Meteorological Society. Used with permission.
Interesting effects occur when trapped lee waves generated by more than one mountain interact. Vosper  studied gravity wave drag on two consecutive 3D mountains using a linear numerical model. He concluded that, when non-hydrostatic effects are unimportant, the drag behaves qualitatively as predicted by Grisogono et al.  for flow over 2D ridges (see section 9). When the atmospheric profiles permit the propagation of trapped lee waves, the drag is strongly enhanced when the waves from an upstream mountain interfere constructively with those from a downstream one. Drag extrema can then be produced even when the mountains are quite far apart. Grubišić and Stiperski  pursued a similar kind of investigation on lee wave interference and resoance, but for pairs of 2D ridges, and using numerical simulations. They found that the appearance of constructive or destructive interference between trapped lee waves generated by each mountain depends both on the separation between them and the ratio of their heights. As expected, constructive and destructive interference lead to high- and low-drag states, respectively. Stiperski and Grubišić  added surface friction to the problem addressed by Grubišić and Stiperski , finding that the important parameter to quantify the drag modulation with distance between the two mountains is the ratio of valley width to the wavelength of the trapped lee waves. However, in contrast with inviscid simulations, the wavelength of the velocity perturbation associated with the waves is shorter than the wavelength associated with the pressure perturbation, which is the one that controls drag behavior.
An important phenomenon associated with trapped lee waves is rotors, which are vortices with their axes of rotation aligned with the mountains, where the flow reverses direction at the surface (see Figure 1). These structures were investigated by Doyle and Durran  using 2D high-resolution numerical simulations. They found that increasing the roughness length in their model decreased the drag, because it tended to shift the point of boundary-layer separation upstream, and therfore reduce the extent of the lee-side (negative) pressure anomaly. A significant correlation between trapped lee wave drag, as given by linear theory, and the occurrence of rotors, as described in the regime diagram of Vosper , was found by Teixeira et al. . That link deserves more detailed investigation.
9. Effects of Rotation
9.1. f-Plane Approximation
For orography with larger width, fa/U = O(1), and the effects of the Earth's rotation may influence mountain waves. The extension of (15) to an atmosphere with an f−plane approximation is:
Note that (65) assumes d2U/dz2 = d2V/dz2 = 0, since this is the only choice consistent with steady-state flow. In a rotating atmosphere, thermal-wind balance implies that vertical shear necessarily corresponds to a mean horizontal temperature gradient, so the potential temperature is no longer steady, because of temperature advection, unless the flow is unidirectional. Consequently, it can be shown (see ) that it is only consistent to assume that N is independent of time and of horizontal position if d2U/dz2 = d2V/dz2 = 0.
Jones  used an equation equivalent to (65) for unidirectional flow (V = 0) to investigate the interaction of inertia-gravity waves with critical levels. He noted that two types of critical levels, where (65) becomes singular, exist: similar to those in non-rotating flow, where U = 0 (or in general Uk + Vl = 0), and inertial critical levels, where Uk = ±f (or in general Uk + Vl = ±f). Jones found that in rotating flow it is the wave angular momentum, rather than the linear momentum, that is conserved in the vertical away from critical levels. He also concluded that absorption of momentum by critical levels differs only slightly between cases with and without rotation. This result was confirmed by Grimshaw , provided that f/N ≪ 1 (which is typically satisfied in the atmosphere). Broad 's generalization of the Eliassen-Palm theorem to 3D flow mentioned in section 6 substantiated these ideas further. The expression valid for rotating flow is still given by (52), but with Mx and My redefined as:
where δ and ζ are the fluid parcel (or streamline) displacements along x and y associated with the waves (an analogous result had been derived previously for 2D rotating flow by Eliassen ).
Wurtele et al. [122, 123] revisited the problem addressed by Jones  using numerical simulations of flow over 2D mountains. They found that while linear theory predicts absorption of a monochromatic wave at levels where |Uk| = f, in nonlinear flow there is also partial reflection. They further showed that both the horizontal and vertical velocity perturbations diverge at these inertial critical levels, making the momentum flux diverge as well. However, for waves with a continuous spectrum, this singular behavior is avoided because each wavenumber has a different inertial critical level, and so momentum absorption is well predicted by linear theory. Wurtele et al.  also found that, for positively sheared flow, the drag and wave momentum flux produced by a monochromatic wave may oscillate in time, with no steady state attained. Shutts  explored further a similar problem using linear theory, but for flow over an axisymmetric mountain, finding that, for a monochromatic wave, only the u and v velocity perturbations diverge at the inertial critical level, while w does not (in contrast with ).
Another relevant aspect of rotation is how it affects the surface drag. Using linear theory and considering unsheared hydrostatic flow over a 2D bell-shaped ridge, Smith  found that mountain wave drag decreases as fa/U increases, being reduced to half of its non-rotating value for fa/U = 0.63 (see Figure 12). While Smith's drag formula can only be expressed analytically in terms of Bessel functions, Miranda and James  obtained a closed expression for the drag in entirely similar conditions, but for flow over an axisymmetric bell-shaped mountain (25), namely
which shows that the drag decreases exponentially for large fa/U. On the other hand, Grisogono et al.  used linear theory to examine the effects of rotation in flow over a double bell-shaped mountain. They derived an expression for the drag as a function of the mountain separation and fa/U, where a is the width of the first mountain, noting that, when f ≠ 0, the drag displays slight oscillations as the ridge spacing varies, with a wavelength that decreases as fa/U increases, due to wave interference. In the same vein, Teixeira et al.  derived asymptotic formulae, expressed in terms of elementary functions, for the mountain wave drag in flow over 2D mountains, when non-hydrostatic effects and rotation are important, but treated as relatively weak. They showed that these approximations are especially accurate for the full range of ã and fa/U when there is a well-defined hydrostatic and non-rotating limit (i.e., f/N ≪ 1), which is typically the case.
Figure 12. Drag normalized by (30) as a function of fa/U for flow over a bell-shaped ridge. Reproduced from Figure 1 of Smith . © American Meteorological Society. Used with permission.
Using hydrostatic numerical simulations of flow with uniform N and U over elliptical mountains, Ólafsson and Bougeault  and Ólafsson  separately analyzed the effects of rotation and surface friction on orographic drag. They found that for low the drag is reduced as fa/U increases but for > 1.4 (blocked flow regime), rotation instead increases the drag, which means that rotation may extend the range of usefulness of linear theory. Additionally they showed that in the non-blocked regime ( = O(1)), the mountain exerts a cyclonic couple on the flow, and the drag oscillates in time (due to wave breaking), whereas in the blocked regime (large ) it exerts an anti-cyclonic couple instead. Petersen et al.  investigated how upstream wind direction influences the flow across a large meridionally-oriented elliptical mountain (representative of Greenland), finding that the drag force and the wake pressure deficit are larger when the flow is from the southwest than when it is from the northwest. For a similar situation, Doyle et al.  examined the sensitivity of wave breaking to terrain slope. They confirmed the result of Ólafsson and Bougeault  that the drag is reduced by rotation for = 0.5, but instead increased by rotation for = 1.5.
The interaction of rotation with shear also has consequences for the surface drag. In his theory of lee cyclogenesis, Smith  showed that quasi-geostrophic flow with constant negative shear and a critical level ((38) with U0 > 0 and α < 0) generates a baroclinic standing wave train on the lee of a 2D ridge. He evaluated the drag exerted by this kind of lee wave as:
where k* = f/(Nzc) is the corresponding resonant wavenumber. Following a similar approach, Martin and Lott  addressed the synoptic response to mountain waves generated by the passage of fronts over orography, and the corresponding absorption of wave momentum flux by directional critical levels. They concluded that the drag exerted by the orography on the flow acts to the left of the wind in a cold front case and to the right of the wind in a warm front case, extending the results of Xu et al.  (see section 6) to rotating flow.
9.2. Beta-Plane Approximation
When the horizontal scale of the mountains is even larger, the variation of the Coriolis parameter with latitude becomes important in the dynamics of orographic waves, which are now Rossby-gravity waves, or simply Rossby waves when the stratification becomes of secondary importance. While the effect of rotation in its simplest form decreases the drag (as was seen in the preceding subsection), the beta effect increases it again.
Janowitz  calculated the drag produced by barotropic Rossby waves in a single-layer flow with a rigid-lid upper boundary condition over a generic isolated mountain. He noted that the drag increases with increasing β = df/dy, and also that the effects of stratification are negligible for this upper boundary condition, because for typical values of the atmospheric parameters no internal (trapped) wave modes are excited.
McCartney  studied inertial Taylor columns on a beta plane, using a two-layer model of flow over a cylindrical mountain. He noted that for eastward flow there was an extensive meandering wake downstream of the mountain, composed of a train of Rossby-gravity waves, and explicitly calculated the associated drag. Bannon  did essentially the same, but for single-layer flow over a Gaussian axisymmetric mountain. Thompson and Flierl  revisited the problem of one-layer flow over a cylindrical mountain, finding that for eastward flow the barotropic Rossby wave drag can be much larger than predicted by quasi-geostrophic theory, being given by:
where D00 = ρ0f0UR2H, f0 is a constant reference Coriolis parameter, R is the radius of the mountain, H is the height of the fluid layer and J1 is the Bessel function of the first kind of order 1.
Teixeira and Grisogono  calculated the drag associated with internal Rossby-gravity waves, by considering a continuously stratified atmosphere flowing over an elliptical mountain. They showed that this drag is especially important for mountains that are aligned meridionally, much larger than the drag given by models where f is constant, substantially larger the the drag produced by internal gravity waves (see Figure 13), and comparable to the barotropic Rossby wave drag calculated by Thompson and Flierl , (70). A more explicit link between these two types of drag is currently missing.
Figure 13. Drag normalized by (31) as a function of the mountain width in the zonal direction for γ = 0.25 and latitudes (ϕ) between 30o and 90o (see legend for details). Reproduced from Figure 3 of Teixeira and Grisogono . Copyright © 2008 Royal Meteorological Society.
10. Nonlinear Effects
Nonlinear effects in mountain waves become important roughly when ≳ O(1). Linear theory, as expressed by (15) or (65), is then not valid anymore, and an alternative method of solution (generally, numerical simulations) must be adopted. There are cases where further progress is possible analytically, however.
Long  derived an equation for stratified flow over a 2D obstacle in non-rotating conditions, which is valid for waves of arbitrary amplitude (if the flow remains steady):
where δ = z − z0 is the vertical displacement of a streamline relative to its upstream undisturbed position z0. When U2 ρ is constant to a good approximation, (71) reduces to:
Equation (72) is remarkable because, while describing nonlinear waves, it is linear (and with constant coefficients if N/U is also constant). Long  applied (71) and (72) to cases of fluid bounded above by a rigid lid. Huppert and Miles  and Miles and Huppert  implemented a radiation boundary condition at the upper and downstream limits of the domain, obtaining solutions to (72) for waves generated by mountains with a semi-elliptical vertical cross-section, and arbitrary shape but relatively small amplitude. Miles and Huppert  calculated the drag for a number of orography forms, including a bell-shaped ridge (29), yielding
correct up to second order in , showing that the drag is underestimated by linear theory. They also found that, even in hydrostatic conditions, the drag is not invariant to flow reversal for asymmetric mountains.
Lilly and Klemp  investigated further the effects of terrain shape and asymmetry using a hydrostatic version of Long's equation, namely
They found that the nonlinear lower and upper (radiative) boundary conditions have a strong influence on the wave structure at large amplitudes, and the drag is significantly enhanced for mountains with a gentle windward slope and a steep leeward slope. Miller and Durran  assessed the sensitivity of downslope windstorms to orography asymmetry using numerical simulations for various atmospheric profiles, finding it to be more important in the case of two-layer flow without a critical level or wave breaking than for flows beneath a critical level or beneath breaking waves. Miller and Durran  also calculated the drag in a single-layer atmosphere using Long's equation, showing that although the drag is super-linear when the mountain is symmetric or has a steep lee slope, it is sub-linear when the upstream slope is the steepest one.
The use of numerical simulations proved to be particularly valuable for exploring the highly nonlinear flow regimes where wave breaking and flow blocking or splitting occur (and thus Long's equation is not valid anymore). Clark and Peltier  performed pioneering numerical simulations of flow over a 2D ridge with constant upstream U and N. They found that while for low the wave drag and momentum flux profile agree with the predictions of linear theory, above a certain threshold there is convective instability, the drag is dramatically enhanced, and the momentum flux profile is strongly divergent (in contrast with what is expected from the Eliassen-Palm theorem), due to wave breaking. In simulations of high-drag states using a hydrostatic primitive equation model, Stein  found that the dimensional drag exerted on a 2D ridge varies proportionally to 2 in the quasi-linear regime, proportionally to a power of larger than 2 in the high-drag regime coinciding with wave breaking, and proportionally to 1.3 in the blocked flow regime. Using a similar model setup to address transient upstream flow effects, Garner  showed that, while the upstream wind weakens gradually as increases, there is a sudden increase in the downslope wind speed at the onset of wave breaking, accompanied by drag enhancement, despite the fact that upstream changes in the flow only become important for much higher . Garner  additionally improved Smith 's formula for the drag in high-drag states, obtaining
where α = 2.2, or
where β = 1/8, when upstream effects become important.
While trying to clarify whether breaking mountain waves decelerate the local mean flow, Durran  showed that zones of flow deceleration in the vicinity of the mountain are compensated by zones of flow acceleration. Thus he suggested that information both about the horizontal and vertical flux of wave momentum need to be parametrized, and the momentum fluxes associated with trapped lee waves should be considered a potentially important source of low-level drag (as was seen in section 8). Mayr and Gohm  revisited the problem of flow over two successive 2D mountains, for = 0.05, 0.51, 1.01, and 1.52, showing that the differing steepness of the upwind and downwind slopes caused by the separating valley strengthens nonlinear effects, and for wide enough valleys two wave breaking regions can form above both peaks, which causes a sharp increase of the drag (by more than a factor of 3 above the linear estimate).
Highly nonlinear flow over 3D mountains (where no equivalent to Long's equation has been derived) was also simulated numerically by various authors. Smolarkiewicz and Rotunno  noted that the conditions for flow reversal on the windward side of a 3D obstacle in inviscid non-rotating flow for ≈ 1 are predicted with some accuracy by linear theory  in the case of axisymmetric mountain, but less so for mountains with different horizontal aspect ratios. Smolarkiewicz and Rotunno  also calculated the drag for an axisymmetric mountain, showing that a maximum several times higher than the linear estimate occurs at ≈ 2. Miranda and James  studied more systematically the drag dependence on (Figure 14, left), finding that the flow is well described by linear theory for < 0.5, for 1 < < 2 it enters the wave-breaking regime, with quasi-periodic oscillations of the drag and an average maximum attained for ≈ 1.3, and for 2 < < 5 it is characterized by streamline splitting around the mountain, a vortex pair on its lee side, and decreasing drag. They concluded that the use of 2D drag estimates may severely overestimate the drag when the flow is in this last regime. Ólafsson and Bougeault  carried out analogous simulations of flow over an elliptical mountain of aspect ratio 5 for 0.500 < < 6.818, showing that for low the linear predictions of flow stagnation aloft (leading to wave breaking) and of stagnation at the surface (leading to flow splitting)  are accurate, but for higher wave breaking occurs mostly at the edges of the mountain, not over its axis of symmetry. Using an essentially similar setup, Bauer et al.  explored more systematically the dependence of the flow on the horizontal aspect ratio of the mountain, finding that, when γ decreases, the threshold for the onset of wave breaking, vortex formation and windward stagnation decreases. The interval of where high-drag states occur, which is centered around ≈ 1, was seen to widen as γ decreases (see Figure 14, right). For a similar flow, Doyle et al.  examined the sensitivity of wave breaking to terrain slope, concluding that that the normalized drag shows a stronger dependence on terrain slope for = 0.5 than for = 1.5, and becomes saturated beyond = 1.5, probably because this dimensionless mountain height is sufficiently high to produce wave breaking regardless of the slope.
Figure 14. Left: drag in flow over an axisymmetric bell-shaped mountain normalized by its linear hydrostatic value (solid line) and by its linear non-hydrostatic value (dashed line) as a function of −1. Reproduced from Figure 1 of Miranda and James . Copyright © 1992 Royal Meteorological Society. Right: normalized drag in flow over an elliptical Gaussian mountain for different β = 1/γ (see legend) as a function of (here denoted by Hm). Reproduced from Figure 15 of Bauer et al. . © American Meteorological Society. Used with permission.
For flows in the splitting regime, the laboratory experiments of strongly stratified flow past axisymmetric obstacles of Vosper et al.  revealed that the wave amplitude decreases, while the drag coefficient increases as increases, reaching values between 2.8 and 5.4 times larger than in the neutral flow limit (depending on obstacle shape) when > 4. A major contribution to the total drag (about 25%) resulted from a separated wake with vortex shedding, rather than from gravity waves.
Eckermann et al.  evaluated the dependence of the momentum fluxes on for flow over 3D obstacles, essentially confirming the findings of Ólafsson and Bougeault  and Bauer et al.  for the surface drag. They additionally showed that the drag vacillates for flow across an elongated obstacle in the 0.7 ≤ < 3 range, due to cyclical buildup and breakdown of wave activity, with the amplitude of these vacillations decreasing as increases. At > 2 − 3, Eckermann et al.  noted that the drag decreased smoothly proportionally to −1/3 to values lower than 1, and the dimensional drag approached a constant value at high , a constraint they suggested implementing in drag parametrizations.
Concerning the variation of the wave momentum flux with height due to wave breaking, Lindzen  formulated the concept of supersaturation of vertically propagating internal gravity waves, which has proved influential in drag parametrization schemes . The underlying idea is that gravity waves, including mountain waves, tend to increase their amplitude as they propagate upward in the atmosphere, due to the decrease of density with height (a non-Boussinesq effect). In practice, when waves reach a saturated state which corresponds to the onset of wave breaking, their amplitude must be limited. This determines the wave momentum flux distribution (and consequently the drag force on the atmosphere) for finite-amplitude waves. Unfortunately, Lindzen's treatment was limited to monochromatic waves. Kim and Mahrt  reported measurements of the momentum flux associated with mountain waves in wave breaking regions from the ALPEX campaign. They attributed the momentum flux divergence in one of the observations to directional wind shear. Using ideas from , they then generalized the wave saturation theory of Lindzen  to vertically varying mean flows, finding that the wave momentum fluxes from the extended model show good agreement with observations from ALPEX.
11. Boundary Layer Effects
Most of the results described above (except essentially those on lee wave rotors in section 8) assumed frictionless flow and a free-slip boundary condition at the surface. A few studies specifically addressing boundary layer effects are mentioned next.
Richard et al.  investigated the effect of friction in downslope windstorms over 2D orography using a hydrostatic numerical model with a turbulent kinetic energy parametrization, finding that it delays the onset of strong surface winds and prevents the downstream propagation of the zone of maximum wind speed. The drag in the corresponding high-states was seen to be established over a longer time interval, with a somewhat lower steady value than in inviscid flow (see Figure 15). Ólafsson and Bougeault  similarly found that in flow over 3D elliptical mountains the drag is substantially reduced by surface friction for < 3, with a suppression of wave breaking and high-drag states (as Jiang et al.  corroborated later). Peng and Thompson  further investigated this effect using a non-hydrostatic numerical model. They hypothesized that the observed reduction in the mountain wave amplitude is due to reduction in the slope of the boundary layer top, and that flow over the mountain outside the boundary layer, including its wave drag, may be reproduced by replacing the terrain elevation with the boundary layer elevation. They found that this effect becomes less pronounced for high because the boundary layer depth occupies then a lower fraction of the relevant flow depth. This was later confirmed by Jiang et al. .
Figure 15. Time evolution of the drag for simulations of flow over 2D orography using free-slip and no-slip lower boundary conditions (see curve labels details). Reproduced from Figure 4 of Richard et al. . © American Meteorological Society. Used with permission.
Belcher and Wood  and Athanassiadou  investigated the relation between form and wave drag in stably stratified flow over 2D hills, using the linear model of Hunt et al.  with the added effect of stable stratification. They showed that the form drag first increases with stability due to increased shear in the boundary layer, but for higher stability it decreases due to boundary layer thinning. The wave drag, which for weak stratification is smaller than the form drag, increases with stratification, eventually becoming dominant. Belcher and Wood  concluded that the reduction of the wave drag relative to its inviscid value is due to the downstream shift that turbulence in the boundary layer imposes on the pressure distribution outside the boundary layer.
Vosper and Brown  investigated the effect of small-scale hills on orographic gravity wave drag and flow blocking produced by 2D and 3D mountains using numerical simulations. They found that corrugations can significantly reduce the amplitude of the mountain waves generated by the broader mountain, and suppress the unsteadiness of the wake. This was seen to be associated with a significant reduction of the total drag compared to the sum of the contributions from the two orography scales, by a factor of more than 30%. Steeneveld et al.  showed that a (realistic) parametrization of drag produced by small-scale orography as partly gravity wave drag instead of drag entirely due to turbulence improved various performance indicators in a forecast model, namely: cyclone filling rates and boundary layer depths, particularly in quite stable conditions.
Nappo and Chimonas  used linear theory to study momentum flux absorption at a critical level within a stably stratified boundary layer, comparing it to frictional drag. They found that gravity wave drag could be comparable, or even larger than frictional drag, and hence should be taken into account in boundary layer drag parametrizations. Grisogono  used linear theory conjugated with a WKB approximation to evaluate the dissipative effect of turbulence in the boundary layer (represented through a constant eddy diffusivity) on the momentum flux profile. He found that the momentum flux decreases exponentially with height due to this effect.
Jiang et al.  investigated the absorption of trapped lee waves due to boundary layers, using linear theory and numerical simulations, finding that they decay exponentially with downstream distance. This decay rate was seen to increase with surface roughness (higher momentum flux), but to decrease as the heat flux increases. Trapped lee wave absorption was therefore found to be maximized for stable boundary layers (in accordance with ). Smith  developed a hydrostatic 3D linear model to represent the absorption of mountain waves by boundary layers, using a bulk approach including a Rayleigh damping coefficient to represent friction. He showed that variations in the boundary layer thickness reduce the mountain wave amplitude (in agreement with ), the pressure drag, and even more severely the wave momentum flux, confirming that part of the pressure drag is due to boundary layer friction.
Using a mesoscale model with boundary layer parametrizations and analytical models, Jiang et al.  showed that the boundary layer tends to reduce the wave drag by up to 60% and the momentum flux above the boundary layer by up to 80%, and slightly delay the onset of wave breaking. They further showed that these effects are stronger over rougher surfaces and for slower boundary layer flow. Jiang and Doyle  found that mountain waves, and in particular the wave drag, exhibit substantial diurnal variation in response to changes in boundary layer depth and stability. During daytime, a convective boundary layer weakens the mountain waves and reduces the surface drag by up to 90%, whereas during nighttime, when the boundary layer is shallow, the drag may increase to several times its free-slip value. This latter phenomenon, which might be due to waves propagating in thin inversion layers existing near the surface (cf. [111, 112]), still awaits a deeper understanding.
12. Concluding Remarks
As was seen in this review, the physical processes affecting mountain wave drag are so numerous, varied and complex, that it is hopeless to include all of them in drag parametrizations. However, it is still useful to study these processes in order to better understand them fundamentally, assess their practical relevance and also the feasibility of parametrizing them. The following list of open questions is by no means exhaustive, being largely dictated by personal preference. It is geared toward a fundamental understanding of mountain wave drag, and the development of simple models useful for parametrization.
Concerning the effect of wind profile shear on the surface drag, results obtained using the WKB approximation [50, 66, 67] have the limitation that they are expressed locally in terms of derivatives of the wind velocity at the surface. The drag behavior is quite sensitive to the particular height where these derivatives are evaluated (which must lie outside the boundary layer, since the theory on which they are based is inviscid). It would be useful to develop an approach where the drag was instead expressed in terms of some integral property of the wind profile, and also where faster variations of the atmospheric parameters, which lead to wave reflections and other non-local effects, could be accommodated.
Critical levels, both total and partial, still need to be better understood, namely the circumstances under which they change from absorbing to reflecting. Perhaps this could be clarified, at least for flow over 2D orography, through the use of weakly-nonlinear theory, more specifically adopting Long's equation including de effect of wind shear. This approach would probably also be able to shed some light on the behavior of nonlinear downslope windstorms, in particular the weakening of intermediate drag maxima in 2D flows with critical levels . Although weakly nonlinear theory has already been applied to 3D problems , it would be interesting to develop a model similar to that of Smith  for downslope windstorms over 3D mountains, something that has not been done until now.
Although a generic expression for trapped lee wave drag analogous to (62) has been derived by Gregory et al. , the dependence of trapped lee wave drag on input flow parameters for simple atmospheric profiles, as investigated by Teixeira et al. ,  for 2D moutains, has not been extended to 3D orograpy, but that is necessary since in real mountains 3D effects are important. It would also be useful to extend the calculations of trapped lee wave drag to more realistic atmospheric profiles, where the windspeed and static stability are not limited to be piecewise constant, as in Teixeira et al. [111, 112] but are, for example, linear, as in Keller  or Shutts , or even more general.
Despite the fact that the drag produced by mountain waves affected by the Earth's rotation is not as relevant for drag parametrization (because it is fully resolved), the predictions of this drag from simple models that account for barotropic Rossby waves  or internal Rossby-gravity waves  have not been verified against numerical simulations. A connection between the first approach (where the drag does not depend on stratification) and the second one (which is an extension of the theory of smaller-scale mountain waves to a beta-plane) needs to be established. It would be useful to know how accurately these two disparate approaches represent the drag in realistic flows.
Finally, there is still much to learn about intrinsically nonlinear effects, such as wave breaking in directionally sheared flows. While the linear analysis of Broad  has shed some light on this issue, further clarification about when and how wave energy is dissipated in such flows over 3D mountains of realistic height is crucial for the improvement of drag parametrization schemes. Currently, the representation of the momentum flux depletion with height that forces the deceleration of the large-scale atmospheric circulation is based on a monochromatic wave formulation (which does not represent the effects of critical level filtering by directional wind shear). An approach akin those suggested by Shutts and Gadian , Teixeira and Miranda  and Teixeira and Yu , or an extension of it to higher wave amplitudes, would improve consistency with current treatments of the surface drag in parameterizations (e.g., ).
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author acknowledges the financial support of the Portuguese Science Foundation (FCT) under Grant PTDC/CTE-ATM/122501/2010 and of the European Commission, under Marie Curie Career Integration Grant GLIMFLO, contract PCIG13-GA-2013-618016.
7. Wallace JM, Tibaldi S, Simmons AJ. Reduction of systematic forecast errors in the ECMWF model through the introduction of an envelope orography. Q J R Meteor Soc. (1983) 109:683–717. doi: 10.1002/qj.49710946202
9. Palmer TN, Shutts GJ, Swinbank R. Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization. Q J R Meteor Soc. (1986) 112:1001–39. doi: 10.1002/qj.49711247406
11. McFarlane NA. The effect of orographically excited gravity-wave drag on the general circulation of the lower stratosphere and troposphere. J Atmos Sci. (1987) 44:1775–1800. doi: 10.1175/1520-0469(1987)044<1775:TEOOEG>2.0.CO;2
13. Kim YJ, Arakawa A. Improvement of orographic gravity wave parameterization using a mesoscale gravity wave model. J Atmos Sci. (1995) 52:1875–902. doi: 10.1175/1520-0469(1995)052<1875:IOOGWP>2.0.CO;2
16. Gregory D, Shutts GJ, Mitchell JR. A new gravity-wave-drag scheme incorporating anisotropic orography and low-level wave breaking. Impact upon the climate of the UK Meteorological Office Unified Model. Q J R Meteor Soc. (1998) 124:463–93. doi: 10.1002/qj.49712454606
26. Madden RA, Speth P. Estimates of Atmospheric Angular momentum, friction, and mountain torques during 1987-1988. J Atmos Sci. (1995) 52:3681–94. doi: 10.1175/1520-0469(1995)052<3681:EOAAMF>2.0.CO;2
34. Grubišić V, Doyle JD, Kuettner J, Mobbs S, Smith RB, Whiteman CD, et al. The terrain-induced rotor experiment: a field campaign overview including observational highlights. Bull Am Meteor Soc. (2008) 89:1513–33. doi: 10.1175/2008BAMS2487.1
35. Doyle JD, Grubišić V, Brown WOJ, Wekker SFJD, Dornbrack A, Jiang QF, et al. Observations and numerical simulations of subrotor vortices during T-REX. J Atmos Sci. (2009) 66:1229–1249. doi: 10.1175/2008JAS2933.1
46. Kim YJ, Eckermann SD, Chun HY. An overview of the past, present and future of gravity-wave drag parametrization for numerical climate and weather prediction models - Survey article. Atmos–Ocean (2003) 41:65–98. doi: 10.3137/ao.410105
48. Lott F, Millet C. The representation of gravity waves in atmospheric general circulation models (CGMs). In: Pinchon AL, Blanc E, Hauchecorne A, editors. Infrasound Monitoring for Atmospheric Studies. Springer (2009). p. 685–699.
50. Teixeira MAC, Miranda PMA, Valente MA. An analytical model of mountain wave drag for wind profiles with shear and curvature. J Atmos Sci. (2004) 61:1040–54. doi: 10.1175/1520-0469(2004)061<1040:AAMOMW>2.0.CO;2
51. Phillips DS. Analytical surface pressure and drag for linear hydrostatic flow over three-dimensional elliptical mountains. J Atmos Sci. (1984) 41:1073–84. doi: 10.1175/1520-0469(1984)041<1073:ASPADF>2.0.CO;2
54. Hines CO. A modeling of atmospheric gravity waves and wave drag generated by isotropic and anisotropic terrain. J Atmos Sci. (1988) 45:309–322. doi: 10.1175/1520-0469(1988)045<0309:AMOAGW>2.0.CO;2
56. Bauer MH, Mayr GJ, Vergeiner I, Pichler H. Strongly nonlinear flow over and around a three-dimensional mountain as a function of the horizontal aspect ratio. J Atmos Sci. (2000) 57:3971–3991. doi: 10.1175/1520-0469(2001)058<3971:SNFOAA>2.0.CO;2
69. Tang JY, Tang J, Wang Y. Analytical investigation on 3D non-Boussinesq mountain wave drag for wind profiles with vertical variations. Appl Math Mech. (2007) 28:317–25. doi: 10.1007/s10483-007-0305-z
74. Scinocca JF, Peltier WR. On the Richardson number dependence of nonlinear critical-layer flow over localized topography. J Atmos Sci. (1991) 48:1560–72. doi: 10.1175/1520-0469(1991)048<1560:OTRNDO>2.0.CO;2
76. Wang TA, Lin YL. Wave ducting in a stratified shear flow over two-dimensional mountains. Part I: general linear criteria. J Atmos Sci. (1999) 56:412–36. doi: 10.1175/1520-0469(1999)056<0412:WDIASS>2.0.CO;2
77. Wang TA, Lin YL. Wave ducting in a stratified shear flow over a two-dimensional mountain. Part II: implications for the development of high-drag states for severe downslope windstorms. J Atmos Sci. (1999) 56:437–52. doi: 10.1175/1520-0469(1999)056<0437:WDIASS>2.0.CO;2
80. Shutts G. Gravity-wave drag parametrization over complex terrain: the effect of critical-level absorption in directional wind-shear. Q J R Meteor Soc. (1995) 121:1005–21. doi: 10.1002/qj.49712152504
87. Xu X, Xue M, Wang Y. Gravity wave momentum flux in directional shear flows over three-dimensional mountains: linear and nonlinear numerical solutions as compared to linear analytical solutions. J Geophys Res. (2013) 118:7670–81. doi: 10.1002/jgrd.50471
88. Teixeira MAC, Yu CL. The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains. Eur J Mech B Fluids (2014). 47:16–31. doi: 10.1016/j.euromechflu.2014.02.004
90. Peltier WR, Clark TL. The evolution and stability of finite-amplitude mountain waves. Part II: surface wave drag and severe downslope windstorms. J Atmos Sci. (1979) 36:1498–1529. doi: 10.1175/1520-0469(1979)036<1498:TEASOF>2.0.CO;2
93. Durran DR. Another look at downslope windstorms. Part I: the development of analogous to supercritical flow in an infinitely deep, continuously stratified fluid. J Atmos Sci. (1986) 43:2527–43. doi: 10.1175/1520-0469(1986)043<2527:ALADWP>2.0.CO;2
100. Teixeira MAC, Miranda PMA, Argaín JL. Mountain waves in two-layer sheared flows: critical level effects, wave reflection and drag enhancement. J Atmos Sci. (2008) 65:1912–26. doi: 10.1175/2007JAS2577.1
102. Wurtele MG, Sharman RD, Keller TL. Analysis and simulations of a troposphere-stratosphere gravity wave model. Part I. J Atmos Sci. (1987) 44:3269–81. doi: 10.1175/1520-0469(1987)044<3269:AASOAT>2.0.CO;2
111. Teixeira MAC, Argaín JL, Miranda PMA. Drag produced by trapped lee waves and propagating mountain and propagating mountain waves in a two-layer atmosphere. Q J R Meteor Soc. (2013) 139:964–81. doi: 10.1002/qj.2008
119. Inverarity GW, Shutts GJ. A general, linearized vertical structure equation for the vertical velocity: properties, scalings and special cases. Q J R Meteor Soc. (2000) 126:2709–24. doi: 10.1002/qj.49712656906
122. Wurtele MG, Datta A, Sharman RD. The propagation of gravity-inertia waves and lee waves under a critical level. J Atmos Sci. (1996) 53:1505–23. doi: 10.1175/1520-0469(1996)053<1505:TPOGWA>2.0.CO;2
130. Petersen GN, Kristjansson JE, Ólafsson H. The effect of upstream wind direction on atmospheric flow in the vicinity of a large mountain. Q J R Meteor Soc. (2005) 131:1113–28. doi: 10.1256/qj.04.01
136. Thompson L, Flierl GR. Barotropic flow over finite isolated topography: steady solutions on the beta-plane and the initial value problem. J Fluid Mech. (1993) 250:553–86. doi: 10.1017/S0022112093001569
144. Stein J. Investigation of the regime diagram of hydrostatic flow over a mountain with a primitive equation model. Part I: Two-dimensional flows. Mon Wea Rev. (1992) 120:2962–76. doi: 10.1175/1520-0493(1992)120<2962:IOTRDO>2.0.CO;2
148. Smolarkiewicz PK, Rotunno R. Low Froude number flow past three-dimensional obstacles. Part II: upwind flow reversal zone. J Atmos Sci. (1990) 47:1498–1511. doi: 10.1175/1520-0469(1990)047<1498:LFNFPT>2.0.CO;2
152. Eckermann S, Lindeman J, Broutman D, Ma J, Boybeyi Z. Momentum fluxes of gravity waves generated by variable Froude number flow over three-dimensional obstacles. J Atmos Sci. (2010) 67:2260–78. doi: 10.1175/2010JAS3375.1
162. Steeneveld GJ, Holtslag AAM, Nappo CJ, de Wiel BJHV, Mahrt L. Exploring the possible role of small-scale terrain drag on stable boundary layers over land. J Appl Meteor Climatol. (2008) 47:2518–30. doi: 10.1175/2008JAMC1816.1
Keywords: mountain waves, drag, momentum flux, stratification, wind shear, coriolis force, nonlinearity, non-hydrostatic effects
Citation: Teixeira MAC (2014) The physics of orographic gravity wave drag. Front. Phys. 2:43. doi: 10.3389/fphy.2014.00043
Received: 12 March 2014; Accepted: 19 June 2014;
Published online: 11 July 2014.
Edited by:Gert-Jan Steeneveld, Wageningen University, Netherlands
Reviewed by:Simon Henri Chabrillat, Belgian Institute for Space Aeronomy, Belgium
Ivana Stiperski, University of Innsbruck, Austria
Copyright © 2014 Teixeira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Miguel A. C. Teixeira, Department of Meteorology, University of Reading, Earley Gate, PO Box 243, Reading RG6 6BB, UK e-mail: firstname.lastname@example.org