Atmospheric internal gravity wave

An atmospheric internal gravity wave is an internal gravity wave propagating through a stably stratified atmosphere. Flow across mountains can generate nearly stationary waves whose vertical propagation is controlled by the wind and buoyancy-frequency profiles.

The intrinsic frequency is the frequency observed in a frame moving with the basic flow. For a wave of frequency $\omega$ and horizontal wavenumber $k$ in a uniform current $U$, it is $\hat{\omega }=\omega -kU$.

A critical level is a height at which the intrinsic frequency of an internal gravity wave vanishes. For finite positive buoyancy frequency, its local vertical wavenumber diverges as $\hat{\omega }\to 0$. The inviscid Taylor–Goldstein equation is singular there; viscosity, diffusion or nonlinear wave breaking can regularize the resulting small scales.

For upward internal gravity waves, let $B=c_{gz}{N}^2{A}_{\zeta }^2/\hat{\omega }$ be constant. Combining $c_{gz}=(N/k)\sin \theta {\cos }^2\theta$ with the monochromatic internal-wave overturning criterion gives ${\cot }^4\theta <B{k}^3\hat{\omega }/(N^3{\sin }^3\theta )$. Near a critical level the onset balance is $\cot \theta \sim (B{k}^3\hat{\omega }/N^3{)}^{1/4}$. This is compatible with $\hat{\omega }=N\cos \theta$; it fixes the overturning point rather than redefining the dispersion relation.

For $S<0$ and constant positive buoyancy frequency, the critical height is $U_0/|S|$. A stationary upward ray approaches it with vertical group velocity proportional to $U^2$, requiring infinite ideal ray travel time. Large gradient Richardson number can preserve the short-wavelength WKB approximation, but parcel-displacement amplitude grows and eventually invalidates linearity.

A radiation condition selects the wave solution whose energy propagates away from its source. For an internal gravity wave, the vertical phase and group velocities have opposite signs, so upward radiation fixes the sign of the vertical wavenumber.

For a slowly switched-on plane internal gravity wave in the half-space $z>0$, the method of multiple scales gives transport of its amplitude at the vertical group velocity. For ${\omega }^2=N^2{k}^2/(k^2+m^2)$, $c_{gz}=-\omega m/(k^2+m^2)$. Initially undisturbed fluid and forcing at the bottom require $c_{gz}>0$. The outgoing radiation condition therefore selects $m<0$ when $\omega >0$, although the vertical phase velocity is downward.

Select outgoing waves as a limit of a problem with weak absorption. With time convention $e^{-i\omega t}$ use $\mathrm{Im}\omega >0$; with $e^{i\omega t}$ use $\mathrm{Im}\omega <0$. The resulting Green function and transform-contour prescriptions determine which poles represent incoming or outgoing waves. This is especially useful when real-axis singularities obscure a Wiener-Hopf factorization.

The wave momentum flux ${\rho }_0\overline{u^{\prime }{w}^{\prime }}$ is the vertical transport of horizontal momentum by correlated velocity perturbations. Its vertical convergence exerts the mean-flow force $-{\rho }_0{\partial }_z\overline{u^{\prime }{w}^{\prime }}$.

The Taylor–Goldstein equation governs linear two-dimensional disturbances of an inviscid, vertically sheared, stably stratified parallel flow. For horizontal phase velocity $c$, base velocity $U(z)$, buoyancy frequency $N(z)$, and horizontal wavenumber $k$,

For positive background velocity $U=U_0{e}^{Mz}$ and buoyancy frequency $N=N_0{e}^{Mz}$, a stationary nonzero-horizontal-wavenumber Taylor–Goldstein equation reduces to the displayed constant-coefficient equation. Its structures are oscillatory, affine or exponential according to the sign of $N_0^2/U_0^2-M^2-k^2$. The disturbance Laplacian is a constant times its streamfunction, but its buoyancy anomaly is $N_0^2{e}^{Mz}\varphi /U_0$. Therefore vorticity self-advection vanishes while buoyancy self-advection produces a nonzero second harmonic for $M\ne 0$. These nontrivial real linear modes are not exact finite-amplitude solutions without additional fields. For $M=0$, both polarizations are constant and these nonlinear self-interactions vanish.

For two-dimensional flow under the Boussinesq approximation with $u={\psi }_z$, $w=-{\psi }_x$, a linear mode whose Laplacian is proportional to its streamfunction has zero self-advection of vorticity. Its buoyancy anomaly is $\sigma =a(z)\varphi$, with $a=N^2/(U-c)$. Self-advection of buoyancy vanishes only when the displayed Jacobian is zero. Constant stratification and constant flow make $a$ constant and yield finite-amplitude exact wave solutions. A varying polarization coefficient generally creates a second harmonic, even when the vorticity self-interaction vanishes.

For a regular real mode between rigid boundaries, multiply the Taylor–Goldstein equation by its vertical structure and integrate to get the stationary quotient $k^2=\int ({\ell }^2{\varphi }^2-{\varphi }_z^2)dz/\int {\varphi }^{2}dz$. Along a differentiable dispersion branch, stationarity removes the variation of the eigenfunction, giving $2k=(\partial I/\partial c)dc/dk$. Since ${\partial }_c{\ell }^2=N^2/(U-c{)}^3+{\ell }^2/(U-c)$, the displayed group velocity follows from $d(kc)/dk$. Critical levels and vanishing branch derivatives require a separate limiting analysis.

For a mode with nonreal phase velocity $c$, let $V=U-c$, choose a continuous branch of $V^a$, and put $\hat{w}=V^{a}q$. Under impermeable boundary conditions, multiplying the transformed Taylor–Goldstein equation by $V^{2a}\overline{q}$ and applying integration by parts givesChoosing $a=1/2$ gives the Miles–Howard theorem; choosing $a=1$ makes the real phase velocity of an unstable mode a weighted mean of $U$.

A smooth inviscid stratified parallel flow with gradient Richardson number at least $1/4$ everywhere has no exponentially growing two-dimensional normal modes. Put $a=1/2$ in the power-transformed Taylor–Goldstein energy identity and take its imaginary part:The integral is positive for a nonzero mode when the numerator is nonnegative, forcing $c_i=0$. This is a modal stability theorem, not a prohibition on transient growth. Maslowe's review discusses the theorem and the role of critical layers.

For an unstable Taylor–Goldstein equation mode in a finite channel, take $a=1$ in the power-transformed Taylor–Goldstein energy identity. Its imaginary part givesThus the real phase velocity lies strictly between the extremes of a nonconstant smooth shear profile. Equality would force the regular eigenfunction to vanish on an interval, hence everywhere by uniqueness for its ordinary differential equation.

At a material interface without surface tension, continuity of displacement and of the pressure evaluated on the displaced interface givesThe first follows from the kinematic boundary condition $ik(U-c)\hat{\eta }=\hat{w}$. The second uses $\hat{p}={\rho }_0[(U-c){\hat{w}}^{\prime }-U^{\prime }\hat{w}]/(ik)$ and hydrostatic displacement $[\hat{p}-g\rho \hat{\eta }]=0$. These conditions apply to density, velocity and vorticity jumps within the Boussinesq approximation.

For ${\psi }^{\prime \prime }={\alpha }^2\psi$ on $-1<z<1$ with endpoint values $P=\psi (1)$ and $Q=\psi (-1)$, the interior derivatives are ${\psi }^{\prime }(1)=\alpha [\coth (2\alpha )P-\mathrm{csch}(2\alpha )Q]$ and ${\psi }^{\prime }(-1)=\alpha [\mathrm{csch}(2\alpha )P-\coth (2\alpha )Q]$. Fitting the hyperbolic-function solution proves the map. It turns matching conditions into a finite-dimensional interface system, with exponentially small cross-interface coupling at large $\alpha$.

At an isolated interface with continuous base velocity, a stable density drop $\Delta \rho /2$ and a jump $[U^{\prime }]$ of base shear, decaying disturbances $w\propto e^{-k|z-z_i|}$ obey the displayed dispersion relation. It follows by inserting their derivative jump $[w^{\prime }]=-2kw$ into the jump conditions for stratified inviscid shear flow. Both the density jump and the vorticity jump contribute to the intrinsic phase velocity $c-U_i$.

Across a horizontal material interface with continuous background mass density but a velocity jump, an internal gravity wave preserves laboratory frequency and horizontal wavenumber. The matching conditions are continuity of displacement $w/\hat{\omega }$ and pressure, not continuity of $w$. For wave amplitudes and upward phase-line angles measured from the vertical ${\theta }_1,{\theta }_2$, define $r={\hat{\omega }}_2/{\hat{\omega }}_1$. Then the transmission coefficient isThe displacement condition is $A_i+A_r=A_t\cos {\theta }_2/(r\cos {\theta }_1)$, while pressure continuity gives $A_i-A_r=r{A}_t\sin {\theta }_2/\sin {\theta }_1$. Adding eliminates the reflected amplitude and gives the coefficient. It assumes nonzero intrinsic frequencies and propagating outgoing branches.

At a material velocity jump with continuous background mass density, vertical displacement $\eta$ and fluid pressure are continuous, while $w=ikU\eta$ need not be. For a stationary component, $p=U^2{\eta }_z$ and the signed upward-wave impedance is $Z=U^{2}m$. Reflection and vertical-displacement transmission are $R=(Z_1-Z_2)/(Z_1+Z_2)$ and $T=2Z_1/(Z_1+Z_2)$.

Two equal stable density jumps of size $\Delta \rho /2$ at $z=\pm h/2$ lie in the global linear shear flow $U=\Delta Uz/h$. Decaying normal modes haveThe jump conditions for stratified inviscid shear flow give the determinant equationWriting $\alpha =kh/2$, $J=g\Delta \rho h/({\rho }_0\Delta U^2)$ and $\tilde{c}=2c/\Delta U$, this isThe two roots for ${\tilde{c}}^2$ are real; one is negative exactly when $2\alpha /(1+e^{-2\alpha })<J<2\alpha /(1-e^{-2\alpha })$. At large $\alpha$, this narrow band centres on $J=2\alpha$, where the isolated counterpropagating interfacial gravity waves have the same zero laboratory speed. This realizes counterpropagating wave instability.

The Scorer parameter is the height-dependent coefficient $l^2$ in the stationary atmospheric-wave equation ${\hat{w}}^{\prime \prime }+(l^2-k^2)\hat{w}=0$. Regions with $l^2>k^2$ support vertically oscillatory disturbances, whereas $l^2<k^2$ makes them vertically evanescent.

For a stationary internal gravity wave, the Taylor–Goldstein equation gives $m^2=N^2/U^2-U^{\prime \prime }/U-k^2$. In a propagating interval with slowly varying $m$, the WKB approximation has amplitude proportional to $m^{-1/2}$ and phase derivative $m$. The conditions $|m^{\prime }|\ll m^2$ and $|m^{\prime \prime }|\ll |m{|}^3$ fail at a regular turning level. With $U>0$, the upward branch has $m>0$ and negative intrinsic frequency; laboratory group velocity is $U(k^2,km)/(k^2+m^2)$ when the curvature term is negligible.

An atmospheric internal gravity wave is vertically trapped when its squared vertical wavenumber changes from positive to negative with altitude. A decrease of the Scorer parameter, caused for example by increasing wind speed or decreasing buoyancy frequency, creates a turning level and an evanescent upper region.

A vertically propagating stationary topographic internal gravity wave meets a turning level when its squared vertical wavenumber changes from positive to negative. The stationary internal-wave WKB solution then fails locally. A simple zero is resolved by an Airy function transition connecting incident and reflected waves below to a decaying evanescent wave above. For negligible wind curvature the condition is $|kU(z_t)|=N(z_t)$, distinct from the zero-intrinsic-frequency critical layer in a shear flow.

For constant buoyancy frequency $N$ and $U=U_0+Sz$ with $S>0$, a stationary component of horizontal wavenumber $k$ turns at $z_t=(N/k-U_0)/S$ if it propagates at the source. Its vertical wavenumber vanishes there; WKB approximation fails and an Airy function transition connects reflection to an upper evanescent wave.