Nonlinear stability of continuously stratified quasi-geostrophic flow

Nonlinear stability theorems analogous to Arnol'd's second stability theorem are established for continuously stratified quasi-geostrophic flow with general nonlinear boundary conditions in a vertically and horizontally confined domain. Both the standard quasi-geostrophic model and the modified quasi-geostrophic model (incorporating effects of hydrostatic compressibility) are treated. The results establish explicit upper bounds on the disturbance energy, the disturbance potential enstrophy, and the disturbance available potential energy on the horizontal boundaries, in terms of the initial disturbance fields. Nonlinear stability in the sense of Liapunov is also established.

Nonlinear stability of Eady's model

A nonlinear stability theorem is established for Eady's model of baroclinic flow. In particular, the Eady basic state is shown to be nonlinearly stable (for arbitrary shear) provided (Δz)/(Δy) > 2(5)^1/2f/(πN),where Δz is the height of the domain, Δy the channel width, f the Coriolis parameter, and N the buoyancy frequency. When this criterion is satisfied, explicit bounds can be derived on the disturbance potential enstrophy, the disturbance energy, and the disturbance available potential energy on the rigid lids, which are expressed in terms of the initial disturbance fields. The disturbances are completely general (with nonzero potential vorticity) and are not assumed to be of small amplitude. The results may be regarded as an extension of Arnol'd's second nonlinear stability theorem to continuously stratified quasigeostrophic baroclinic flow.

Nonlinear saturation of baroclinic instability: part II: continuously stratified fluid

Rigorous upper bounds are derived that limit the finite-amplitude growth of arbitrary nonzonal disturbances to an unstable baroclinic zonal flow in a continuously stratified, quasi-geostrophic, semi-infinite fluid. Bounds are obtained bath on the depth-integrated eddy potential enstrophy and on the eddy available potential energy (APE) at the ground. The method used to derive the bounds is essentially analogous to that used in Part I of this study for the two-layer model: it relies on the existence of a nonlinear Liapunov (normed) stability theorem, which is a finite-amplitude generalization of the Charney-Stern theorem. As in Part I, the bounds are valid both for conservative (unforced, inviscid) flow, as well as for forced-dissipative flow when the dissipation is proportional to the potential vorticity in the interior, and to the potential temperature at the ground. The character of the results depends on the dimensionless external parameter γ = f02ξ/β0N2H, where ξ is the maximum vertical shear of the zonal wind, H is the density scale height, and the other symbols have their usual meaning. When γ ≫ 1, corresponding to “deep” unstable modes (vertical scale ≈H), the bound on the eddy potential enstrophy is just the total potential enstrophy in the system; but when γ≪1, corresponding to ‘shallow’ unstable modes (vertical scale ≈γH), the eddy potential enstrophy can be bounded well below the total amount available in the system. In neither case can the bound on the eddy APE prevent a complete neutralization of the surface temperature gradient which is in accord with numerical experience. For the special case of the Charney model of baroclinic instability, and in the limit of infinitesimal initial eddy disturbance amplitude, the bound states that the dimensionless eddy potential enstrophy cannot exceed (γ + 1)2/24&gamma2h when γ ≥ 1, or 1/6;&gammah when γ ≤ 1; here h = HN/f0L is the dimensionless scale height and L is the width of the channel. These bounds are very similar to (though of course generally larger than) ad hoc estimates based on baroclinic-adjustment arguments. The possibility of using these kinds of bounds for eddy-amplitude closure in a transient-eddy parameterization scheme is also discussed.

Nonlinear wave-activity conservation laws and Hamiltonian structure for the two-dimensional anelastic equations

Exact, finite-amplitude, local wave-activity conservation laws are derived for disturbances to steady flows in the context of the two-dimensional anelastic equations. The conservation laws are expressed entirely in terms of Eulerian quantities, and have the property that, in the limit of a small-amplitude, slowly varying, monochromatic wave train, the wave-activity density A and flux F, when averaged over phase, satisfy F = cgA where cg is the group velocity of the waves. For nonparallel steady flows, the only conserved wave activity is a form of disturbance pseudoenergy; when the steady flow is parallel, there is in addition a conservation law for the disturbance pseudomomentum. The above results are obtained not only for isentropic background states (which give the so-called “deep form” of the anelastic equations), but also for arbitrary background potential-temperature profiles θ0(z) so long as the variation in θ0(z) over the depth of the fluid is small compared with θ0 itself. The Hamiltonian structure of the equations is established in both cases, and its symmetry properties discussed. An expression for available potential energy is also derived that, for the case of a stably stratified background state (i.e., dθ0/dz > 0), is locally positive definite; the expression is valid for fully three-dimensional flow. The counterparts to results for the two-dimensional Boussinesq equations are also noted.

Large-scale two-dimensional turbulence in the atmosphere

Global FGGE data are used to investigate several aspects of large-scale turbulence in the atmosphere. The approach follows that for two-dimensional, nondivergent turbulent flows which are homogeneous and isotropic on the sphere. Spectra of kinetic energy, enstrophy and available potential energy are obtained for both the stationary and transient parts of the flow. Nonlinear interaction terms and fluxes of energy and enstrophy through wavenumber space are calculated and compared with the theory. A possible method of parameterizing the interactions with unresolved scales is considered. Two rather different flow regimes are found in wavenumber space. The high-wavenumber regime is dominated by the transient components of the flow and exhibits, at least approximately, several of the conditions characterizing homogeneous and isotropic turbulence. This region of wavenumber space also displays some of the features of an enstrophy-cascading inertial subrange. The low-wavenumber region, on the other hand, is dominated by the stationary component of the flow, exhibits marked anisotropy and, in contrast to the high-wavenumber regime, displays a marked change between January and July.

Measurements of NO, NOy, N2O, and O3 during SPURT: implications for transport and chemistry in the lowermost stratosphere

We present measurements of NO, NOy, O3, and N2O within the lowermost stratosphere (LMS) over Europe obtained during the SPURT project. The measurements cover all seasons between November 2001 and July 2003. They span a broad band of latitudes from 30° N to 75° N and a potential temperature range from 290 to 380 K. The measurements represent a comprehensive data set of these tracers and reveal atmospheric transport processes that influence tracer distributions in the LMS. Median mixing ratios of stratospheric tracers in equivalent latitude-potential temperature coordinates show a clear seasonal cycle related to the Brewer-Dobson circulation, with highest values in spring and lowest values in autumn. Vertical tracer profiles show strong gradients at the extratropical tropopause, suggesting that vertical (cross-isentropic) mixing is reduced above the tropopause. Pronounced meridional gradients in the tracer mixing ratios are found on potential temperature surfaces in the LMS. This suggests strongly reduced mixing along isentropes. Concurrent large gradients in static stability in the vertical direction, and of PV in the meridional direction, suggest the presence of a mixing barrier. Seasonal cycles were found in the correlation slopes ΔO3/ΔN2O and ΔNOy/ΔN2O well above the tropopause. Absolute slope values are smallest in spring indicating chemically aged stratospheric air originating from high altitudes and latitudes. Larger values were measured in summer and autumn suggesting that a substantial fraction of air takes a "short-cut" from the tropical tropopause region into the extratropical LMS. The seasonal change in the composition of the LMS has direct implications for the ozone chemistry in this region. Comparisons of measured NO with the critical NO value at which net ozone production changes from negative to positive, imply ozone production up to 20 K above the local tropopause in spring, up to 30 K in summer, and up to 40 K in autumn. Above these heights, and in winter, net ozone production is negative.

Irreversible compressible work and APE dissipation in turbulent stratified fluid

Although it plays a key role in the theory of stratified turbulence, the concept of available potential energy (APE) dissipation has remained until now a rather mysterious quantity, owing to the lack of rigorous result about its irreversible character or energy conversion type. Here, we show by using rigorous energetics considerations rooted in the analysis of the Navier-Stokes for a fully compressible fluid with a nonlinear equation of state that the APE dissipation is an irreversible energy conversion that dissipates kinetic energy into internal energy, exactly as viscous dissipation. These results are established by showing that APE dissipation contributes to the irreversible production of entropy, and by showing that it is a part of the work of expansion/contraction. Our results provide a new interpretation of the entropy budget, that leads to a new exact definition of turbulent effective diffusivity, which generalizes the Osborn-Cox model, as well as a rigorous decomposition of the work of expansion/contraction into reversible and irreversible components. In the context of turbulent mixing associated with parallel shear flow instability, our results suggests that there is no irreversible transfer of horizontal momentum into vertical momentum, as seems to be required when compressible effects are neglected, with potential consequences for the parameterisations of momentum dissipation in the coarse-grained Navier-Stokes equations.

Thermodynamics/dynamics coupling in weakly compressible turbulent stratified fluids

In traditional and geophysical fluid dynamics, it is common to describe stratified turbulent fluid flows with low Mach number and small relative density variations by means of the incompressible Boussinesq approximation. Although such an approximation is often interpreted as decoupling the thermodynamics from the dynamics, this paper reviews recent results and derive new ones that show that the reality is actually more subtle and complex when diabatic effects and a nonlinear equation of state are retained. Such an analysis reveals indeed: (1) that the compressible work of expansion/contraction remains of comparable importance as the mechanical energy conversions in contrast to what is usually assumed; (2) in a Boussinesq fluid, compressible effects occur in the guise of changes in gravitational potential energy due to density changes. This makes it possible to construct a fully consistent description of the thermodynamics of incompressible fluids for an arbitrary nonlinear equation of state; (3) rigorous methods based on using the available potential energy and potential enthalpy budgets can be used to quantify the work of expansion/contraction B in steady and transient flows, which reveals that B is predominantly controlled by molecular diffusive effects, and act as a significant sink of kinetic energy.

A GCSS model intercomparison for a tropical squall line observed during toga-coare. II: Intercomparison of single-column models and a cloud-resolving model

This paper presents single-column model (SCM) simulations of a tropical squall-line case observed during the Coupled Ocean-Atmosphere Response Experiment of the Tropical Ocean/Global Atmosphere Programme. This case-study was part of an international model intercomparison project organized by Working Group 4 ‘Precipitating Convective Cloud Systems’ of the GEWEX (Global Energy and Water-cycle Experiment) Cloud System Study. Eight SCM groups using different deep-convection parametrizations participated in this project. The SCMs were forced by temperature and moisture tendencies that had been computed from a reference cloud-resolving model (CRM) simulation using open boundary conditions. The comparison of the SCM results with the reference CRM simulation provided insight into the ability of current convection and cloud schemes to represent organized convection. The CRM results enabled a detailed evaluation of the SCMs in terms of the thermodynamic structure and the convective mass flux of the system, the latter being closely related to the surface convective precipitation. It is shown that the SCMs could reproduce reasonably well the time evolution of the surface convective and stratiform precipitation, the convective mass flux, and the thermodynamic structure of the squall-line system. The thermodynamic structure simulated by the SCMs depended on how the models partitioned the precipitation between convective and stratiform. However, structural differences persisted in the thermodynamic profiles simulated by the SCMs and the CRM. These differences could be attributed to the fact that the total mass flux used to compute the SCM forcing differed from the convective mass flux. The SCMs could not adequately represent these organized mesoscale circulations and the microphysicallradiative forcing associated with the stratiform region. This issue is generally known as the ‘scale-interaction’ problem that can only be properly addressed in fully three-dimensional simulations. Sensitivity simulations run by several groups showed that the time evolution of the surface convective precipitation was considerably smoothed when the convective closure was based on convective available potential energy instead of moisture convergence. Finally, additional SCM simulations without using a convection parametrization indicated that the impact of a convection parametrization in forced SCM runs was more visible in the moisture profiles than in the temperature profiles because convective transport was particularly important in the moisture budget.

Weak intermolecular interactions in an ionically bound molecular adsorbate: cyclopentadienyl=Cu(111)

The dissociative adsorption of cyclopentadiene (C5H6) on Cu(111) yields a cyclopentadienyl (Cp) species with strongly anionic characteristics. The Cp potential energy surface and frictional coupling to the substrate are determined from measurements of dynamics of the molecule together with density functional calculations. The molecule is shown to occupy degenerate threefold adsorption sites and molecular motion is characterized by a low diffusional energy barrier of 40 +/- 3 meV with strong frictional dissipation. Repulsive dipole-dipole interactions are not detected despite charge transfer from substrate to adsorbate.

Mode-specific chemisorption of CH4 on Pt{110}-(1 x 2) explored by first-principles molecular dynamics

The chemisorption of CH4 on Pt{110}-(1 x 2) has been studied by vibrational analysis of the reaction pathway defined by the potential energy surface and, in time reversal, by first-principles molecular dynamics simulations of CH4 associative desorption, with the electronic structure treated explicitly using density functional theory. We find that the symmetric stretch vibration ν1 is strongly coupled to the reaction coordinate; our results therefore provide a firm theoretical basis for recently reported state-resolved reactivity measurements, which show that excitation of the ν1 normal mode is the most efficient way to enhance the reaction probability

Estimating Lorenz's reference state in an ocean with a nonlinear equation of state for seawater

The study of the mechanical energy budget of the oceans using Lorenz available potential energy (APE) theory is based on knowledge of the adiabatically re-arranged Lorenz reference state of minimum potential energy. The compressible and nonlinear character of the equation of state for seawater has been thought to cause the reference state to be ill-defined, casting doubt on the usefulness of APE theory for investigating ocean energetics under realistic conditions. Using a method based on the volume frequency distribution of parcels as a function of temperature and salinity in the context of the seawater Boussinesq approximation, which we illustrate using climatological data, we show that compressibility effects are in fact minor. The reference state can be regarded as a well defined one-dimensional function of depth, which forms a surface in temperature, salinity and density space between the surface and the bottom of the ocean. For a very small proportion of water masses, this surface can be multivalued and water parcels can have up to two statically stable levels in the reference density profile, of which the shallowest is energetically more accessible. Classifying parcels from the surface to the bottom gives a different reference density profile than classifying in the opposite direction. However, this difference is negligible. We show that the reference state obtained by standard sorting methods is equivalent, though computationally more expensive, to the volume frequency distribution approach. The approach we present can be applied systematically and in a computationally efficient manner to investigate the APE budget of the ocean circulation using models or climatological data.

Impact of Interbasin Exchange on the Atlantic Overturning Circulation

The ther mohaline exchange between the Atlantic and the Souther n Ocean is analyzed, using a dataset based on WOCE hydrographic data. It is shown that the salt and heat transports brought about by the South Atlantic subtropical gyre play an essential role in the Atlantic heat and salt budgets. It is found that on average the exported North Atlantic Deep W ater (NADW) is fresher than the retur n flows (basically composed of ther mocline and inter mediate water), indicating that the overtur ning circulation (OC) exports freshwater from the Atlantic. The sensitivity of the OC to interbasin fluxes of heat and salt is studied in a 2 D model, representing the Atlantic between 60 8 N and 30 8 S. The model is forced by mixed boundar y conditions at the sur face, and by realistic fluxes of heat and salt at its 30 8 S boundar y. The model circulation tur ns out to be ver y sensitive to net buoyancy fluxes through the sur face. Both net sur face cooling and net sur face saltening are sources of potential energy and impact positively on the circulation strength. The vertical distributions of the lateral fluxes tend to stabilize the stratification, and, as they extract potential energy from the system, tend to weaken the flow . These results imply that a change in the composition of the NADW retur n transports, whether by a change in the ratio ther mocline/inter mediate water , o r by a change in their ther mohaline characteristics, might influence the Atlantic OC considerably . It is also shown that the circulation is much more sensitive to changes in the shape of the lateral buoyancy flux than to changes in the shape of the sur face buoyancy flux, as the latter does not explicitly impact on the potential energy of the system. It is concluded that interocean fluxes of heat and salt are important for the strength and operation of the Atlantic ther mohaline circulation, and should be correctly represented in models that are used for climate sensitivity studies.

Re-examining the U.K.'s greatest tornado outbreak: forecasting the limited extent of tornadoes along a cold front

On 23 November 1981, a strong cold front swept across the U.K., producing tornadoes from the west to the east coasts. An extensive campaign to collect tornado reports by the Tornado and Storm Research Organisation (TORRO) resulted in 104 reports, the largest U.K. outbreak. The front was simulated with a convection-permitting numerical model down to 200-m horizontal grid spacing to better understand its evolution and meteorological environment. The event was typical of tornadoes in the U.K., with convective available potential energy (CAPE) less than 150 J kg-1, 0-1-km wind shear of 10-20 m s-1, and a narrow cold-frontal rainband forming precipitation cores and gaps. A line of cyclonic absolute vorticity existed along the front, with maxima as large as 0.04 s-1. Some hook-shaped misovortices bore kinematic similarity to supercells. The narrow swath along which the line was tornadic was bounded on the equatorward side by weak vorticity along the line and on the poleward side by zero CAPE, enclosing a region where the environment was otherwise favorable for tornadogenesis. To determine if the 104 tornado reports were plausible, first possible duplicate reports were eliminated, resulting in as few as 58 tornadoes to as many as 90. Second, the number of possible parent misovortices that may have spawned tornadoes is estimated from model output. The number of plausible tornado reports in the 200-m grid-spacing domain was 22 and as many as 44, whereas the model simulation was used to estimate 30 possible parent misovortices within this domain. These results suggest that 90 reports was plausible.

Conservative force fields in non-Gaussian statistics

In this Letter, we determine the kappa-distribution function for a gas in the presence of an external field of force described by a potential U(r). In the case of a dilute gas, we show that the kappa-power law distribution including the potential energy factor term can rigorously be deduced in the framework of kinetic theory with basis on the Vlasov equation. Such a result is significant as a preliminary to the discussion on the role of long range interactions in the Kaniadakis thermostatistics and the underlying kinetic theory. (C) 2008 Elsevier B.V. All rights reserved.