An improvement in clear-air turbulence forecasting based on spontaneous imbalance theory: the ULTURB algorithm

Recent research has shown that Lighthill–Ford spontaneous gravity wave generation theory, when applied to numerical model data, can help predict areas of clear-air turbulence. It is hypothesized that this is the case because spontaneously generated atmospheric gravity waves may initiate turbulence by locally modifying the stability and wind shear. As an improvement on the original research, this paper describes the creation of an ‘operational’ algorithm (ULTURB) with three modifications to the original method: (1) extending the altitude range for which the method is effective downward to the top of the boundary layer, (2) adding turbulent kinetic energy production from the environment to the locally produced turbulent kinetic energy production, and, (3) transforming turbulent kinetic energy dissipation to eddy dissipation rate, the turbulence metric becoming the worldwide ‘standard’. In a comparison of ULTURB with the original method and with the Graphical Turbulence Guidance second version (GTG2) automated procedure for forecasting mid- and upper-level aircraft turbulence ULTURB performed better for all turbulence intensities. Since ULTURB, unlike GTG2, is founded on a self-consistent dynamical theory, it may offer forecasters better insight into the causes of the clear-air turbulence and may ultimately enhance its predictability.

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

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

Linear criteria for gravity-wave breaking in resonant stratified flow over a ridge

Using linear theory, it is shown that, in resonant flow over a 2D mountain ridge, such as exists when a layer of uniform wind is topped by an environmental critical level, the conditions for internal gravity-wave breaking are different from those determined in previous studies for non-resonant flows. For Richardson numbers in the shear layer not exceeding 2.25, two zones of flow overturning exist, respectively below and downstream and above and upstream of the expected locations. Flow overturning occurs for values of the dimensionless height of the ridge smaller than those required for a uniform wind profile. These results may have implications for the physical understanding of high-drag states.

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.

Measuring motion with kinematically redundant accelerometer arrays: theory, simulation and implementation

This work presents two schemes of measuring the linear and angular kinematics of a rigid body using a kinematically redundant array of triple-axis accelerometers with potential applications in biomechanics. A novel angular velocity estimation algorithm is proposed and evaluated that can compensate for angular velocity errors using measurements of the direction of gravity. Analysis and discussion of optimal sensor array characteristics are provided. A damped 2 axis pendulum was used to excite all 6 DoF of the a suspended accelerometer array through determined complex motion and is the basis of both simulation and experimental studies. The relationship between accuracy and sensor redundancy is investigated for arrays of up to 100 triple axis (300 accelerometer axes) accelerometers in simulation and 10 equivalent sensors (30 accelerometer axes) in the laboratory test rig. The paper also reports on the sensor calibration techniques and hardware implementation.

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.

Instabilities in thermal gravity with a cosmological constant

It is shown that in quantum gravity at finite temperature, the effective potential evaluated in the tadpole approximation can have a local minimum below a certain critical temperature. However, when the leading higher order thermal loop corrections are included, one finds that no static solution exists at high temperature. (C) 2008 Elsevier B.V. All rights reserved.

Fermion perturbations in string theory black holes

In this paper we study fermion perturbations in four-dimensional black holes of string theory, obtained either from a non-extreme configuration of three intersecting five-branes with a boost along the common string or from a non-extreme intersecting system of two two-branes and two five-branes. The Dirac equation for the massless neutrino field, after conformal re-scaling of the metric, is written as a wave equation suitable to study the time evolution of the perturbation. We perform a numerical integration of the evolution equation, and with the aid of Prony fitting of the time-domain profile, we calculate the complex frequencies that dominate the quasinormal ringing stage, and also determine these quantities by the semi-analytical sixth-order WKB method. We also find numerically the decay factor of fermion fields at very late times, and show that the falloff is identical to those showing for massless fields in other four-dimensional black hole spacetimes.

Quadratic gravity in (2+1)D with a topological Chern-Simons term

Three-dimensional quadratic gravity, unlike general relativity in (2+1)D, is dynamically nontrivial and has a well behaved nonrelativistic potential. Here we analyse the changes that occur when a topological Chem-Simons term is added to this theory. It is found that the harmless massive scalar mode of the latter gives rise to a troublesome massive spin-0 ghost, while the massive spin-2 ghost is replaced by two massive physical particles both of spin 2. We also found that light deflection does not have the 'wrong sign' such as in the framework of three-dimensional quadratic gravity.

Can one cure the nommitarity of three-dimensional quadratic gravity using a topological Chern-Simons term?

Quadratic gravity in (2+1)D is nonunitarity at the tree level. When a topological Chern-Simons term is added to this theory, the harmless massive scalar mode of the former gives rise to a troublesome massive spin-0 ghost, while the massive spin-2 ghost is replaced by two massive physical particles both of spin-2. Therefore, unlike what it is claimed in the literature, quadratic Chern-Simons gravity in (2+1)D is nonunitary at the tree level.

Consistent modified gravity: dark energy, acceleration and the absence of cosmic doomsday

We discuss modified gravity which includes negative and positive powers of curvature and provides gravitational dark energy. It is shown that in GR plus a term containing a negative power of curvature, cosmic speed-up may be achieved while the effective phantom phase (with w less than -1) follows when such a term contains a fractional positive power of curvature. Minimal coupling with matter makes the situation more interesting: even 1/R theory coupled with the usual ideal fluid may describe the (effective phantom) dark energy. The account of the R(2) term (consistent modified gravity) may help to escape cosmic doomsday.

Prediction of R plus R-2 gravity for the deflection of a photon passing close to the Sun

The cross-section for the scattering of a photon by the Sun's gravitational field, treated as an external field, is computed in the framework of R + R-2 gravity. Using this result, we found that for a photon just grazing the Sun's surface the deflection is 1.75 which is exactly the same as that given by Einstein's theory. An explanation for this pseudo-paradox is provided.

Propagator, tree-level unitarity and effective nonrelativistic potential for three-dimensional higher-derivative gravity augmented by a gravitational Chern-Simons term

An algorithm for computing the propagator for three-dimensional quadratic gravity with a gravitational Chern-Simons term, based on an extension of the three-dimensional Barnes-Rivers operators, is proposed. A systematic study of the tree-level unitarity of this theory is developed and its agreement with Newton's law is investigated by computing the effective nonrelativistic potential. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Quadratic gravity theories in 2+1 dimensions and the topological Chern-Simons term

The addition of a topological Chern-Simons term to three-dimensional higher-derivative gravity is not a good therapy to cure the nonunitarity of the aforementioned theory. Moreover, R+R-2 gravity in (2+1)D, which is unitary at the tree level, becomes tree-level nonunitary when it is augmented by the abovementioned topological term. Therefore, unlike what is claimed in the literature, topological higher-derivative gravity in (2+1)D is not tree-level unitary and neither is topological three-dimensional R+R-2 gravity.