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

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

Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part I: gravity wave forcing

Parameterization schemes for the drag due to atmospheric gravity waves are discussed and compared in the context of a simple one-dimensional model of the quasi-biennial oscillation (QBO). A number of fundamental issues are examined in detail, with the goal of providing a better understanding of the mechanism by which gravity wave drag can produce an equatorial zonal wind oscillation. The gravity wave–driven QBOs are compared with those obtained from a parameterization of equatorial planetary waves. In all gravity wave cases, it is seen that the inclusion of vertical diffusion is crucial for the descent of the shear zones and the development of the QBO. An important difference between the schemes for the two types of waves is that in the case of equatorial planetary waves, vertical diffusion is needed only at the lowest levels, while for the gravity wave drag schemes it must be included at all levels. The question of whether there is downward propagation of influence in the simulated QBOs is addressed. In the gravity wave drag schemes, the evolution of the wind at a given level depends on the wind above, as well as on the wind below. This is in contrast to the parameterization for the equatorial planetary waves in which there is downward propagation of phase only. The stability of a zero-wind initial state is examined, and it is determined that a small perturbation to such a state will amplify with time to the extent that a zonal wind oscillation is permitted.

Orographic drag associated with lee waves trapped at an inversion

The drag produced by 2D orographic gravity waves trapped at a temperature inversion and waves propagating in the stably stratified layer existing above are explicitly calculated using linear theory, for a two-layer atmosphere with neutral static stability near the surface, mimicking a well-mixed boundary layer. For realistic values of the flow parameters, trapped lee wave drag, which is given by a closed analytical expression, is comparable to propagating wave drag, especially in moderately to strongly non-hydrostatic conditions. In resonant flow, both drag components substantially exceed the single-layer hydrostatic drag estimate used in most parametrization schemes. Both drag components are optimally amplified for a relatively low-level inversion and Froude numbers Fr ≈ 1. While propagating wave drag is maximized for approximately hydrostatic flow, trapped lee wave drag is maximized for l_2 a = O(1) (where l_2 is the Scorer parameter in the stable layer and a is the mountain width). This roughly happens when the horizontal scale of trapped lee waves matches that of the mountain slope. The drag behavior as a function of Fr for l_2 H = 0.5 (where H is the inversion height) and different values of l2a shows good agreement with numerical simulations. Regions of parameter space with high trapped lee wave drag correlate reasonably well with those where lee wave rotors were found to occur in previous nonlinear numerical simulations including frictional effects. This suggests that trapped lee wave drag, besides giving a relevant contribution to low-level drag exerted on the atmosphere, may also be useful to diagnose lee rotor formation.

Dynamical control of the mesosphere by orographic and non-orographic gravity wave drag during the extended northern winters of 2006 and 2009

A version of the Canadian Middle Atmosphere Model (CMAM) that is nudged toward reanalysis data up to 1 hPa is used to examine the impacts of parameterized orographic and non-orographic gravity wave drag (OGWD and NGWD) on the zonal-mean circulation of the mesosphere during the extended northern winters of 2006 and 2009 when there were two large stratospheric sudden warmings. The simulations are compared to Aura Microwave Limb Sounder (MLS) observations of mesospheric temperature, carbon monoxide (CO) and derived zonal winds. The control simulation, which uses both OGWD and NGWD, is shown to be in good agreement with MLS. The impacts of OGWD and NGWD are assessed using simulations in which those sources of wave drag are removed. In the absence of OGWD the mesospheric zonal winds in the months preceding the warmings are too strong, causing increased mesospheric NGWD, which drives excessive downwelling, resulting in overly large lower mesospheric values of CO prior to the warming. NGWD is found to be most important following the warmings when the underlying westerlies are too weak to allow much vertical propagation of the orographic gravity waves to the mesosphere. NGWD is primarily responsible for driving the circulation that results in the descent of CO from the thermosphere following the warmings. Zonal mean mesospheric winds and temperatures in all simulations are shown to be strongly constrained by (i.e. slaved to) the stratosphere. Finally, it is demonstrated that the responses to OGWD and NGWD are non-additive due to their dependence and influence on the background winds and temperatures.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

The physics of orographic gravity wave drag

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.

An integrable evolution equation for surface waves in deep water

In order to describe the dynamics of monochromatic surface waves in deep water, we derive a nonlinear and dispersive system of equations for the free surface elevation and the free surface velocity from the Euler equations in infinite depth. From it, and using a multiscale perturbative method, an asymptotic model for small wave steepness ratio is derived. The model is shown to be completely integrable. The Lax pair, the first conserved quantities as well as the symmetries are exhibited. Theoretical and numerical studies reveal that it supports periodic progressive Stokes waves which peak and break in finite time. Comparison between the limiting wave solution of the asymptotic model and classical results is performed.

Infragravity waves in the nearshore zone /

Includes bibliographical references: (p. 48-55).

Aircraft hazard assessment from a clear-air radar and meteorological tower study of gravity wave events /

"June 1983."

The effect of waves on the ideal fluid forces acting on a slender axisymmetric body having three degrees of freedom /

Mode of access: Internet.

Idealized modeling of the role of stability and shear on mesoscale gravity wave evolution

Mesoscale Gravity Waves (MGWs) are large pressure perturbations that form in the presence of a stable layer at the surface either behind Mesoscale Convective Systems (MCSs) in summer or over warm frontal surfaces behind elevated convection in winter. MGWs are associated with damaging winds, moderate to heavy precipitation, and occasional heat bursts at the surface. The forcing mechanism for MGWs in this study is hypothesized to be evaporative cooling occurring behind a convective line. This evaporatively-cooled air generates a downdraft that then depresses the surface-based stable layer and causes pressure decreases, strong wind speeds and MGW genesis. Using the Weather Research and Forecast Model (WRF) version 3.0, evaporative cooling is simulated using an imposed cold thermal. Sensitivity studies examine the response of MGW structure to different thermal and shear profiles where the strength and depth of the inversion are varied, as well as the amount of wind shear. MGWs are characterized in terms of response variables, such as wind speed perturbations (U'), temperature perturbations (T'), pressure perturbations (P'), potential temperature perturbations (Θ'), and the correlation coefficient (R) between U' and P'. Regime Diagrams portray the response of MGW to the above variables in order to better understand the formation, causes, and intensity of MGWs. The results of this study indicate that shallow, weak surface layers coupled with deep, neutral layers above favor the formation of waves of elevation. Conversely, deep strong surface layers coupled with deep, neutral layers above favor the formation of waves of depression. This is also the type of atmospheric setup that tends to produce substantial surface heating at the surface.

Beryllium abundances along the evolutionary sequence of the open cluster IC 4651-A new test for hydrodynamical stellar models

Context. Previous analyses of lithium abundances in main sequence and red giant stars have revealed the action of mixing mechanisms other than convection in stellar interiors. Beryllium abundances in stars with Li abundance determinations can offer valuable complementary information on the nature of these mechanisms. Aims. Our aim is to derive Be abundances along the whole evolutionary sequence of an open cluster. We focus on the well-studied open cluster IC 4651. These Be abundances are used with previously determined Li abundances, in the same sample stars, to investigate the mixing mechanisms in a range of stellar masses and evolutionary stages. Methods. Atmospheric parameters were adopted from a previous abundance analysis by the same authors. New Be abundances have been determined from high-resolution, high signal-to-noise UVES spectra using spectrum synthesis and model atmospheres. The careful synthetic modeling of the Be lines region is used to calculate reliable abundances in rapidly rotating stars. The observed behavior of Be and Li is compared to theoretical predictions from stellar models including rotation-induced mixing, internal gravity waves, atomic diffusion, and thermohaline mixing. Results. Beryllium is detected in all the main sequence and turn-off sample stars, both slow- and fast-rotating stars, including the Li-dip stars, but is not detected in the red giants. Confirming previous results, we find that the Li dip is also a Be dip, although the depletion of Be is more modest than for Li in the corresponding effective temperature range. For post-main-sequence stars, the Be dilution starts earlier within the Hertzsprung gap than expected from classical predictions, as does the Li dilution. A clear dispersion in the Be abundances is also observed. Theoretical stellar models including the hydrodynamical transport processes mentioned above are able to reproduce all the observed features well. These results show a good theoretical understanding of the Li and Be behavior along the color-magnitude diagram of this intermediate-age cluster for stars more massive than 1.2 M(circle dot).

CMB in a box: Causal structure and the Fourier-Bessel expansion

This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads to position-space solutions which only depend on the sources of anisotropies inside the past light cone of the observer. This foretold manifestation of causality in position (real) space happens order by order in a series expansion in powers of the visibility gamma = e(-mu), where mu is the optical depth to Thomson scattering. We show that the contributions of order gamma(N) to the cosmic microwave background (CMB) anisotropies are regulated by spacetime window functions which have support only inside the past light cone of the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields (including the temperature and polarization momenta) is an alternative to the usual Fourier basis as a framework to compute the anisotropies. The viability of the Fourier-Bessel series for treating the CMB is a consequence of the fact that the visibility function becomes exponentially small at redshifts z >> 10(3), effectively cutting off the past light cone and introducing a finite radius inside which initial conditions can affect physical observables measured at our position (x) over right arrow = 0 and time t(0). Hence, for each multipole l there is a discrete tower of momenta k(il) (not a continuum) which can affect physical observables, with the smallest momenta being k(1l) similar to l. The Fourier-Bessel modes take into account precisely the information from the sources of anisotropies that propagates from the initial value surface to the point of observation-no more, no less. We also show that the physical observables (the temperature and polarization maps), and hence the angular power spectra, are unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal scheme with which one can compute CMB anisotropies.

Associations involving delays (particularly long delays) between certain weather parameters and geomagnetic activity

Four sunspot-minimum periods (1963-1966, 1971-1977, 1983-1987 and 1992-1997) have been examined for the results which are presented. Using several different weather parameters, tropospheric gravity waves, enhanced cold fronts and two rainfall data sets in Eastern Australia, associations at reasonably high levels of significance have been found with enhanced geomagnetic activity (EGA). Statistically this EGA involved either short delays of several days or long delays of about 20 days. The geomagnetic parameters used were (a) the AE index (b) the hourly H component for a number of stations and (c) the daily K-P-sum value. The K-P-sum analyses have shown that the EGA associated with the delays form part of four or five cycles of recurrent geomagnetic activity for 27-day periodicities. Furthermore statistically two recurrent cycles are found to exist concurrently, one apparently related to the short delays and the other to the long delays. Periodicities of 13.5 days are created because the two sets are displaced from each other by approximately this interval. A brief reference is made to the 13.5 periodicity known to exist for geomagnetic activity and the evidence in the literature for active regions on the sun to be displaced by 180 degrees of solar longitude.