Improving a global model from the boundary layer: total turbulent energy and the neutral limit Prandtl number

Model intercomparisons have identified important deficits in the representation of the stable boundary layer by turbulence parametrizations used in current weather and climate models. However, detrimental impacts of more realistic schemes on the large-scale flow have hindered progress in this area. Here we implement a total turbulent energy scheme into the climate model ECHAM6. The total turbulent energy scheme considers the effects of Earth’s rotation and static stability on the turbulence length scale. In contrast to the previously used turbulence scheme, the TTE scheme also implicitly represents entrainment flux in a dry convective boundary layer. Reducing the previously exaggerated surface drag in stable boundary layers indeed causes an increase in southern hemispheric zonal winds and large-scale pressure gradients beyond observed values. These biases can be largely removed by increasing the parametrized orographic drag. Reducing the neutral limit turbulent Prandtl number warms and moistens low-latitude boundary layers and acts to reduce longstanding radiation biases in the stratocumulus regions, the Southern Ocean and the equatorial cold tongue that are common to many climate models.

Comments on 'Balance and the slow quasimanifold: some explicit results'

The concept of slow vortical dynamics and its role in theoretical understanding is central to geophysical fluid dynamics. It leads, for example, to “potential vorticity thinking” (Hoskins et al. 1985). Mathematically, one imagines an invariant manifold within the phase space of solutions, called the slow manifold (Leith 1980; Lorenz 1980), to which the dynamics are constrained. Whether this slow manifold truly exists has been a major subject of inquiry over the past 20 years. It has become clear that an exact slow manifold is an exceptional case, restricted to steady or perhaps temporally periodic flows (Warn 1997). Thus the concept of a “fuzzy slow manifold” (Warn and Ménard 1986) has been suggested. The idea is that nearly slow dynamics will occur in a stochastic layer about the putative slow manifold. The natural question then is, how thick is this layer? In a recent paper, Ford et al. (2000) argue that Lighthill emission—the spontaneous emission of freely propagating acoustic waves by unsteady vortical flows—is applicable to the problem of balance, with the Mach number Ma replaced by the Froude number F, and that it is a fundamental mechanism for this fuzziness. They consider the rotating shallow-water equations and find emission of inertia–gravity waves at O(F2). This is rather surprising at first sight, because several studies of balanced dynamics with the rotating shallow-water equations have gone beyond second order in F, and found only an exponentially small unbalanced component (Warn and Ménard 1986; Lorenz and Krishnamurthy 1987; Bokhove and Shepherd 1996; Wirosoetisno and Shepherd 2000). We have no technical objection to the analysis of Ford et al. (2000), but wish to point out that it depends crucially on R 1, where R is the Rossby number. This condition requires the ratio of the characteristic length scale of the flow L to the Rossby deformation radius LR to go to zero in the limit F → 0. This is the low Froude number scaling of Charney (1963), which, while originally designed for the Tropics, has been argued to be also relevant to mesoscale dynamics (Riley et al. 1981). If L/LR is fixed, however, then F → 0 implies R → 0, which is the standard quasigeostrophic scaling of Charney (1948; see, e.g., Pedlosky 1987). In this limit there is reason to expect the fuzziness of the slow manifold to be “exponentially thin,” and balance to be much more accurate than is consistent with (algebraic) Lighthill emission.

Carbon suboxide: the infrared spectrum from 1800 to 2600 cm-1

The infrared spectrum of carbon suboxide has been recorded from 1800 to 2600 cm−1 at a resolution of 0.003 cm−1. About 7% of the ca. 40 000 lines observed have been assigned and analyzed, belonging to 36 different bands. Most of these are associated with the fundamental ν3, at 2289.80 cm−1, and the combination band ν2 + ν4, at 2386.61 cm−1, each of which give rise to a system of sum bands, difference bands, and hot bands involving the low-wave-number fundamental ν7 at 18 cm−1. A few other tentative assignments are made. The bands have been analyzed for vibrational and rotational constants.

Correlations of control variables in variational data assimilation

Variational data assimilation systems for numerical weather prediction rely on a transformation of model variables to a set of control variables that are assumed to be uncorrelated. Most implementations of this transformation are based on the assumption that the balanced part of the flow can be represented by the vorticity. However, this assumption is likely to break down in dynamical regimes characterized by low Burger number. It has recently been proposed that a variable transformation based on potential vorticity should lead to control variables that are uncorrelated over a wider range of regimes. In this paper we test the assumption that a transform based on vorticity and one based on potential vorticity produce an uncorrelated set of control variables. Using a shallow-water model we calculate the correlations between the transformed variables in the different methods. We show that the control variables resulting from a vorticity-based transformation may retain large correlations in some dynamical regimes, whereas a potential vorticity based transformation successfully produces a set of uncorrelated control variables. Calculations of spatial correlations show that the benefit of the potential vorticity transformation is linked to its ability to capture more accurately the balanced component of the flow.

Stability of stationary solutions of the forced Navier-Stokes equations on the two-torus

We study the linear and nonlinear stability of stationary solutions of the forced two-dimensional Navier-Stokes equations on the domain [0,2π]x[0,2π/α], where α ϵ(0,1], with doubly periodic boundary conditions. For the linear problem we employ the classical energy{enstrophy argument to derive some fundamental properties of unstable eigenmodes. From this it is shown that forces of pure χ2-modes having wavelengths greater than 2π do not give rise to linear instability of the corresponding primary stationary solutions. For the nonlinear problem, we prove the equivalence of nonlinear stability with respect to the energy and enstrophy norms. This equivalence is then applied to derive optimal conditions for nonlinear stability, including both the high-and low-Reynolds-number limits.

Infrared dynamics of decaying two-dimensional turbulence governed by the Charney-Hasegawa-Mima equation

The low wave number range of decaying turbulence governed by the Charney-Hasegawa-Mima (CHM) equation is examined theoretically and by direct numerical simulation. Here, the low wave number range is defined as values of the wave number k below the wave number kE corresponding to the peak of the energy spectrum, or alternatively the centroid wave number of the energy spectrum. The energy spectrum in the low wave number range in the infrared regime (k →0) is theoretically derived to be E(k) ∼k5, using a quasinormal Markovianized model of the CHM equation. This result is verified by direct numerical simulation of the CHM equation. The wave number triads (k,p,q) responsible for the formation of the low wave number spectrum are also examined. It is found that the energy flux Π(k) for k< kE can be entirely expressed by Π(-)(k), which is the total net input of energy to wave numbers numbers p, q > k. Furthermore, the contribution of nonlocal triad interactions to the energy flux is found to be predominant in the range log (k/kE)<-0.5, where the nonlocal interactions are defined to be those triad interactions for which the ratio of the largest leg of the triad to the smallest leg is larger than four. ©2001 The Physical Society of Japan

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.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low-shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional-wave velocity is strongly correlated with the shear-wave velocity in this suite of rocks. The ratio V-P/V-S varies significantly with change of both pore-fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot-flow and squirt-flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt-flow of fluid between the pores of various geometries. The contribution of the Biot-flow loss mechanism to the total loss is negligible. The compressional-wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low-shale, high-porosity sandstone reservoirs.

Low density lipoprotein undergoes oxidation within lysosomes in cells

The oxidised low density lipoprotein (LDL) hypothesis of atherosclerosis proposes that LDL undergoes oxidation in the interstitial fluid of the arterial wall. We have shown that aggregated (vortexed) nonoxidised LDL was taken up by J774 mouse macrophages and human monocyte-derived macrophages and oxidised intracellularly, as assessed by the microscopic detection of ceroid, an advanced lipid oxidation product. Confocal microscopy showed that the ceroid was located in the lysosomes. To confirm these findings, J774 macrophages were incubated with acetylated LDL, which is internalised rapidly to lysosomes, and then incubated (chase incubation) in the absence of any LDL. The intracellular levels of oxysterols, measured by HPLC, increased during the chase incubation period, showing that LDL must have been oxidised inside the cells. Furthermore, we found that this oxidative modification was inhibited by lipid-soluble antioxidants, an iron chelator taken up by fluid-phase pinocytosis and the lysosomotropic drug chloroquine, which increases the pH of lysosomes. The results indicate that LDL oxidation can occur intracellularly, most probably within lysosomes.

Nitrogen fertilizer and seed rate effects on Hagberg failing number of hybrid wheats and their parents are associated with alpha-amylase activity, grain cavity size and dormancy

Field experiments were carried out to assess the effects of nitrogen fertilization and seed rate on the Hagberg falling number (HFN) of commercial wheat hybrids and their parents. Applying nitrogen (200 kg N ha(-1)) increased HFN in two successive years. The HFN of the hybrid Hyno Esta was lower than either of its parents (Estica and Audace), particularly when nitrogen was not applied. Treatment effects on HFN were negatively associated with a-amylase activity. Phadebas grain blotting suggested two populations of grains with different types of a-amylase activity: Estica appeared to have a high proportion of grains with low levels of late maturity endosperm a-amylase activity (LMEA); Audace had a few grains showing high levels of germination amylase; and the hybrid, Hyno Esta, combined the sources from both parents to show heterosis for a-amylase activity. Applying nitrogen reduced both apparent LMEA and germination amylase. The effects on LMEA were associated with the size and disruption of the grain cavity, which was greater in Hyno Esta and Estica and in zero-nitrogen treatments. External grain morphology failed to explain much of the variation in LMEA and cavity size, but there was a close negative correlation between cavity size and protein content. Applying nitrogen increased post-harvest dormancy of the grain. Dormancy was greatest in Estica and least in Audace. It is proposed that effects of seed rate, genotype and nitrogen fertilizer on HFN are mediated through factors affecting the size and disruption of the grain cavity and therefore LMEA, and through factors affecting dormancy and therefore germination amylase. (c) 2004 Society of Chemical Industry.

Long-term effectiveness of sowing high and low diversity seed mixtures to enhance plant community development on ex-arable fields

Questions: How is succession on ex-arable land affected by sowing high and low diversity mixtures of grassland species as compared to natural succession? How long do effects persist? Location: Experimental plots installed in the Czech Republic, The Netherlands, Spain, Sweden and the United Kingdom. Methods: The experiment was established on ex-arable land, with five blocks, each containing three 10 m x 10 m experiment tal plots: natural colonization, a low- (four species) and high-diversity (15 species) seed mixture. Species composition and biomass was followed for eight years. Results: The sown plants considerably affected the whole successional pathway and the effects persisted during the whole eight year period. Whilst the proportion of sown species (characterized by their cover) increased during the study period, the number of sown species started to decrease from the third season onwards. Sowing caused suppression of natural colonizing species, and the sown plots had more biomass. These effects were on average larger in the high diversity mixtures. However, the low diversity replicate sown with the mixture that produced the largest biomass or largest suppression of natural colonizers fell within the range recorded at the five replicates of the high diversity plots. The natural colonization plots usually had the highest total species richness and lowest productivity at the end of the observation period. Conclusions: The effect of sowing demonstrated dispersal limitation as a factor controlling the rate of early secondary succession. Diversity was important primarily for its 'insurance effect': the high diversity mixtures were always able to compensate for the failure of some species.

Alternative methodologies - stretching the in vitro box

Current gas-based in vitro evaluation systems are extremely powerful research techniques. However they have the potential to generate a great deal more than simple fermentation dynamics. Details from four experiments are presented in which adaptation, and novel application, of an in vitro system allowed widely differing objectives to be examined. In the first two studies, complement methodologies were utilised. In such assays, an activity or outcome is inferred through the occurrence of a secondary event rather than by direct observation. Using an N-deficient incubation medium, the increase in starch fermentation, when supplemented with individual amino acids (i.e., known level of N) relative to that of urea (i.e., known quantity and N availability), provided an estimate of their microbial utilisation. Due to the low level of response observed with some arnino acids (notably methionine and lysine), it was concluded, that they may not need to be offered in a rumen-inert form to escape rumen microbial degradation. In another experiment, the extent to which degradation of plant cell wall components was inhibited by lipid supplementation was evaluated using fermentation gas release profiles of washed hay. The different responses due to lipid source and level of inclusion suggested that the degree of rumen protection required to ameliorate this depression was supplement dependent. That in vitro inocula differ in their microbial composition is of little interest per se, as long as the outcome is the same (i.e., that similar substrates are degraded at comparable rates and end-product release is equivalent). However where a microbial population is deficient in a particular activity, increasing the level of inoculation will have no benefit. Estimates of hydrolytic activity were obtained by examining fermentation kinetics of specific substrates. A number of studies identified a fundamental difference between rumen fluid and faecal inocula, with the latter having a lower fibrolytic activity, which could not be completely attributed to microbial numbers. The majority of forage maize is offered as an ensiled feed, however most of the information on which decisions such as choice of variety, crop management and harvesting date are made is based on fresh crop measurements. As such, an attempt was made to estimate ensiled maize quality from an in vitro analysis of the fresh crop. Fermentation profiles and chemical analysis confirmed changes in crop composition over the growing season, and loss of labile carbohydrates during ensiling. In addition, examination of degradation residues allowed metabolizable energy (ME) contents to be estimated. Due to difficulties associated with starch analysis, the observation that this parameter could be predicted by difference (together with an assumed degradability), allowed an estimate of ensiled maize ME to be developed from fresh material. In addition, the contribution of the main carbohydrates towards ME showed the importance of delaying harvest until maximum starch content has been achieved. (c) 2005 Elsevier B.V. All rights reserved.

Low molecular weight heparin significantly reduces embolisation after carotid endarterectomy - a randomised controller trial

Objectives The administration of unfractionated heparin (UFH) prior to carotid clamping during carotid endarterectomy (CEA) transiently increases the platelet aggregation response to arachidonic acid (AA) despite the use of aspirin. We hypothesized that this phenomenon might be reduced by using low molecular weight heparin (LMWH) resulting in fewer emboli in the early post-operative period. Methods 183 aspirinated patients undergoing CEA were randomised to 5000 IU UFH (n = 91) or 2500 IU LMWH (dalteparin, n = 92) prior to carotid clamping. End-points were: transcranial Doppler (TCD) measurement of embolisation, effect on bleeding and platelet aggregation to AA and adenosine 5′-diphosphate (ADP). Results Patients randomised to UFH had twice the odds of experiencing a higher number of emboli in the first 3 h after CEA, than those randomised to LMWH (p = 0.04). This was not associated with increased bleeding (mean time from flow restoration to operation end: 23 min (UFH) vs. 24 min (LMWH), p = 0.18). Platelet aggregation to AA increased significantly following heparinisation, but was unaffected by heparin type (p = 0.90). The platelets of patients randomised to LMWH exhibited significantly lower aggregation to ADP compared to UFH (p < 0.0001). Conclusions Intravenous LMWH is associated with a significant reduction in post-operative embolisation without increased bleeding. The higher rate of embolisation seen with UFH may be mediated by increased platelet aggregation to ADP, rather than to AA.

Low density lipoprotein undergoes oxidation within lysosomes in cells

The oxidized low density lipoprotein (LDL) hypothesis of atherosclerosis proposes that LDL undergoes oxidation in the interstitial fluid of the arterial wall. We have shown that aggregated (vortexed) nonoxidized LDL was taken up by J774 mouse macrophages and human monocyte-derived macrophages and oxidized intracellularly, as assessed by the microscopic detection of ceroid, an advanced lipid oxidation product. Confocal microscopy showed that the ceroid was located in the lysosomes. To confirm these findings, J774 macrophages were incubated with acetylated LDL, which is internalized rapidly to lysosomes, and then incubated (chase incubation) in the absence of any LDL. The intracellular levels of oxysterols, measured by HPLC, increased during the chase incubation period, showing that LDL must have been oxidized inside the cells. Furthermore, we found that this oxidative modification was inhibited by lipid-soluble antioxidants, an iron chelator taken up by fluid-phase pinocytosis and the lysosomotropic drug chloroquine, which increases the pH of lysosomes. The results indicate that LDL oxidation can occur intracellularly, most probably within lysosomes.