Optimising the FAMOUS climate model: inclusion of global carbon cycling

FAMOUS fills an important role in the hierarchy of climate models, both explicitly resolving atmospheric and oceanic dynamics yet being sufficiently computationally efficient that either very long simulations or large ensembles are possible. An improved set of carbon cycle parameters for this model has been found using a perturbed physics ensemble technique. This is an important step towards building the "Earth System" modelling capability of FAMOUS, which is a reduced resolution, and hence faster running, version of the Hadley Centre Climate model, HadCM3. Two separate 100 member perturbed parameter ensembles were performed; one for the land surface and one for the ocean. The land surface scheme was tested against present-day and past representations of vegetation and the ocean ensemble was tested against observations of nitrate. An advantage of using a relatively fast climate model is that a large number of simulations can be run and hence the model parameter space (a large source of climate model uncertainty) can be more thoroughly sampled. This has the associated benefit of being able to assess the sensitivity of model results to changes in each parameter. The climatologies of surface and tropospheric air temperature and precipitation are improved relative to previous versions of FAMOUS. The improved representation of upper atmosphere temperatures is driven by improved ozone concentrations near the tropopause and better upper level winds.

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.

Available potential energy density for a multicomponent Boussinesq fluid with arbitrary nonlinear equation of state

In this paper, the concept of available potential energy (APE) density is extended to a multicomponent Boussinesq fluid with a nonlinear equation of state. As shown by previous studies, the APE density is naturally interpreted as the work against buoyancy forces that a parcel needs to perform to move from a notional reference position at which its buoyancy vanishes to its actual position; because buoyancy can be defined relative to an arbitrary reference state, so can APE density. The concept of APE density is therefore best viewed as defining a class of locally defined energy quantities, each tied to a different reference state, rather than as a single energy variable. An important result, for which a new proof is given, is that the volume integrated APE density always exceeds Lorenz’s globally defined APE, except when the reference state coincides with Lorenz’s adiabatically re-arranged reference state of minimum potential energy. A parcel reference position is systematically defined as a level of neutral buoyancy (LNB): depending on the nature of the fluid and on how the reference state is defined, a parcel may have one, none, or multiple LNB within the fluid. Multiple LNB are only possible for a multicomponent fluid whose density depends on pressure. When no LNB exists within the fluid, a parcel reference position is assigned at the minimum or maximum geopotential height. The class of APE densities thus defined admits local and global balance equations, which all exhibit a conversion with kinetic energy, a production term by boundary buoyancy fluxes, and a dissipation term by internal diffusive effects. Different reference states alter the partition between APE production and dissipation, but neither affect the net conversion between kinetic energy and APE, nor the difference between APE production and dissipation. We argue that the possibility of constructing APE-like budgets based on reference states other than Lorenz’s reference state is more important than has been previously assumed, and we illustrate the feasibility of doing so in the context of an idealised and realistic oceanic example, using as reference states one with constant density and another one defined as the horizontal mean density field; in the latter case, the resulting APE density is found to be a reasonable approximation of the APE density constructed from Lorenz’s reference state, while being computationally cheaper.

Efficient fully nonlinear data assimilation for geophysical fluid dynamics

A potential problem with Ensemble Kalman Filter is the implicit Gaussian assumption at analysis times. Here we explore the performance of a recently proposed fully nonlinear particle filter on a high-dimensional but simplified ocean model, in which the Gaussian assumption is not made. The model simulates the evolution of the vorticity field in time, described by the barotropic vorticity equation, in a highly nonlinear flow regime. While common knowledge is that particle filters are inefficient and need large numbers of model runs to avoid degeneracy, the newly developed particle filter needs only of the order of 10-100 particles on large scale problems. The crucial new ingredient is that the proposal density cannot only be used to ensure all particles end up in high-probability regions of state space as defined by the observations, but also to ensure that most of the particles have similar weights. Using identical twin experiments we found that the ensemble mean follows the truth reliably, and the difference from the truth is captured by the ensemble spread. A rank histogram is used to show that the truth run is indistinguishable from any of the particles, showing statistical consistency of the method.

Protection of dried probiotic bacteria from bile using bile adsorbent resins

Enteric coated oral tablets or capsules can deliver dried live cells directly into the intestine. Previously, we found that a live attenuated bacterial vaccine acquired sensitivity to intestinal bile when dried, raising the possibility that although gastric acid can be bypassed, significant loss of viability might occur on release from an enteric coated oral formulations. Here we demonstrate that some food-grade lyophilised preparations of Lactobacillus casei and Lactobacillus salivarius also show temporary bile sensitivity that can be rapidly reversed by rehydration. To protect dried bacterial cells from temporary bile sensitivity, we propose using bile acid adsorbing resins, such as cholestyramine, which are bile acid binding agents, historically used to lower cholesterol levels. Vcaps™ HPMC capsules alone provided up to 830-fold protection from bile. The inclusion of 50% w/w cholestyramine in Vcaps™ HPMC capsules resulted in release of up to 1700-fold more live Lactobacillus casei into simulated intestinal fluid containing 1% bile, when compared to dried cells added directly to bile. We conclude that delivery of dried live probiotic organisms to the intestine may be improved by providing protection from bile by addition of bile adsorbing resins and the use of HPMC capsules.

The effect of adding stinging nettle (Urtica dioica) haylage to a total mixed ration on performance and rumen function of lactating dairy cows

In vitro studies found that inclusion of dried stinging nettle (Urtica dioica) at 100 mg/g dry matter (DM) increased the pH of a rumen fluid inoculated fermentation buffer by 30% and the effect was persistent for 7 days. Our objective was to evaluate the effects of adding stinging nettle haylage to a total mixed ration on feed intake, eating and rumination activity, rumen pH, milk yield, and milk composition of lactating dairy cows. Six lactating Holstein-Friesian cows were used in a replicated 3 × 3 Latin Square design experiment with 3 treatments and 3 week periods. Treatments were a control (C) high-starch (311 g/kg DM) total mixed ration diet and two treatment diets containing 50 (N5) and 100 (N10) g nettle haylage (DM/kg) as a replacement for ryegrass silage (Lolium perenne). There was an increase (linear, P < 0.010) in the proportion of large particles and a reduction in medium (linear, P = 0.045) and fine particles (linear, P = 0.026) in the diet offered with increasing nettle inclusion. A numerical decrease (linear, P = 0.106) in DM intake (DMI) was observed as nettle inclusion in the diet increased. Milk yield averaged 20.3 kg/day and was not affected by diet. There was a decrease (quadratic, P = 0.01) in the time animals spent ruminating as nettle inclusion in the diet increased, in spite of an increase in the number of boli produced daily for the N5 diet (quadratic, P = 0.031). Animals fed the N10 diet spent less time with a rumen pH below 5.5 (P < 0.05) than cows fed the N5 diet. Averaged over an 8.5 h sampling period, there were no changes in the concentration or proportions of acetate or propionate in the rumen, but feeding nettle haylage reduced the concentrations of n-butyrate (quadratic, P < 0.001), i-butyrate (linear, P < 0.009) and n-caproate (linear, P < 0.003). Milk and fat and protein corrected milk yield were not affected when nettles replaced ryegrass silage in the diet of lactating dairy cows, despite a numerical reduction in feed intake. Rumination activity was reduced by the addition of nettle haylage to the diet, which may reflect differences in fibre structure between the nettle haylage and ryegrass silage fed. Changes observed in rumen pH suggest potential benefits of feeding nettle haylage for reducing rumen acidosis. However, the extent to which these effects were due to the fermentability and structure of the nettle haylage compared to the ryegrass silage fed, or a bioactive component of the nettles, is not certain

Effect of inclusion of lactulose and Lactobacillus plantarum on the intestinal environment and performance of piglets at weaning

The prebiotic lactulose, a probiotic strain of Lactobacillus plantarum (L. plantarum) and a synbiotic combination of these two agents were evaluated as growth promoters in 25–39-day old commercial weaning pigs. Ninety-six weaning pigs were allocated into 32 pens, taking initial weight into account, and distributed into four groups as follows: a control diet (CTR), the same diet supplemented daily with L. plantarum (109 CFU/mL sprayed on top; 20 mL/pig) (LPN); 10 g/kg lactulose (LAC) or a combination of both treatments (SYN). At day 14, eight piglets from each group were euthanized and proximal colon digesta was sampled for luminal pH, short-chain fatty acids (SCFA) and lactic acid concentrations. Deoxyribonucleic acid was extracted from colonic digesta and the microbial community was profiled by terminal restriction fragment length polymorphism analysis (T-RFLP) and qPCR. Blood urea nitrogen (BUN) and acute-phase proteins (Pig-MAP) were measured. Lactulose treatment (LAC) improved feed intake (P<0.05), average daily gain (P<0.01), feed:gain ratio (P<0.05) and reduced BUN (P<0.01). Both, LAC and LPN treatment, decreased the Enterobacteriaceae:Lactobacillus spp. ratio in the colonic luminal contents (P<0.05). Moreover LPN treatment promoted a decrease in the percentage of branched fatty acids (P<0.01) suggesting a reduction in proteolytic microbial activity. Microbial profiling of colonic luminal contents by T-RFLP revealed changes in some microbial species. Terminal restriction fragments (TRFs) compatible with Bifidobacterium thermoacidophilum were more frequently detected in experimental diets compared to CTR (P<0.05). Pigs receiving SYN diet demonstrated the combined positive effects of individual LAC and LPN treatment although we were not able to show a specific increase in the probiotic strain with the inclusion of lactulose. Collectively, these data suggest the combination of lactulose and L. plantarum acts as a complementary synbiotic, but not as a synergistic combination.

Tracking (poly)phenol components from raspberries in ileal fluid

The (poly)phenols in ileal fluid after ingestion of raspberries were analysed by targeted and non-targeted LC-MSn approaches. Targeted approaches identified major anthocyanin and ellagitannin components at varying recoveries and with considerable inter-individual variation. Non-targeted LC-MSn analysis using an Orbitrap mass spectrometer gave exact mass MS data which was sifted using a software program to select peaks that changed significantly after supplementation. This method confirmed the recovery of the targeted components but also identified novel raspberry-specific metabolites. Some components (including ellagitannin and previously unidentified proanthocyanidin derivatives) may have arisen from raspberry seeds that survived intact in ileal samples. Other components include potential breakdown products of anthocyanins, unidentified components and phenolic metabolites formed in either the gut epithelia or after absorption into the circulatory system and efflux back into the gut lumen. The possible physiological roles of the ileal metabolites in the large bowel are discussed.

Solution calorimetry as a tool for investigating drug interaction with intestinal fluid

Solution calorimetry offers a reproducible technique for measuring the enthalpy of solution (ΔsolH) of a solute dissolving into a solvent. The ΔsolH of two solutes, propranolol HCl and mannitol were determined in simulated intestinal fluid (SIF) solutions designed to model the fed and fasted states within the gut, and in Hanks’ balanced salt solution (HBSS) of varying pH. The bile salt and lipid within the SIF solutions formed mixed micelles. Both solutes exhibited endothermic reactions in all solvents. The ΔsolH for propranolol HCl in the SIF solutions differed from those in the HBSS and was lower in the fed state than the fasted state SIF solution, revealing an interaction between propranolol and the micellar phase in both SIF solutions. In contrast, for mannitol the ΔsolH was constant in all solutions indicating minimal interaction between mannitol and the micellar phases of the SIF solutions. In this study, solution calorimetry proved to be a simple method for measuring the enthalpy associated with the dissolution of model drugs in complex biological media such as SIF solutions. In addition, the derived power–time curves allowed the time taken for the powdered solutes to form solutions to be estimated.

The inclusion of forage mixtures in the diet of growing dairy heifers: impacts on digestion, energy utilisation, and methane emissions

Intensive farming focusing on monoculture grass species to maximise forage production has led to a reduction in the extent and diversity of species-rich grasslands. However, plant communities with higher species number (richness) are a potential strategy for more sustainable production and mitigation of greenhouse gas (GHG) emissions. Research has indicated the need to understand opportunities that forage mixtures can offer sustainable ruminant production systems. The objective of the two experiments reported here were to evaluate multiple species forage mixtures in comparison to ryegrass-dominant pasture, when conserved or grazed, on digestion, energy utilisation, N excretion, and methane emissions by growing 10–15 month old heifers. Experiment 1 was a 4 × 4 Latin square design with five week periods. Four forage treatments of: (1) ryegrass (control); permanent pasture with perennial ryegrass (Lolium perenne); (2) clover; a ryegrass:red clover (Trifolium pratense) mixture; (3) trefoil; a ryegrass:birdsfoot trefoil (Lotus corniculatus) mixture; and (4) flowers; a ryegrass:wild flower mixture of predominately sorrel (Rumex acetosa), ox-eye daisy (Leucanthemum vulgare), yarrow (Achillea millefolium), knapweed (Centaurea nigra) and ribwort plantain (Plantago lanceolata), were fed as haylages to four dairy heifers. Measurements included digestibility, N excretion, and energy utilisation (including methane emissions measured in respiration chambers). Experiment 2 used 12 different dairy heifers grazing three of the same forage treatments used to make haylage in experiment 1 (ryegrass, clover and flowers) and methane emissions were estimated using the sulphur hexafluoride (SF6) tracer technique. Distribution of ryegrass to other species (dry matter (DM) basis) was approximately 70:30 (clover), 80:20 (trefoil), and 40:60 (flowers) for experiment 1. During the first and second grazing rotations (respectively) in experiment 2, perennial ryegrass accounted for 95 and 98% of DM in ryegrass, and 84 and 52% of DM in clover, with red clover accounting for almost all of the remainder. In the flowers mixture, perennial ryegrass was 52% of the DM in the first grazing rotation and only 30% in the second, with a variety of other flower species occupying the remainder. Across both experiments, compared to the forage mixtures (clover, trefoil and flowers), ryegrass had a higher crude protein (CP) content (P < 0.001, 187 vs. 115 g kg −1 DM) and DM intake (P < 0.05, 9.0 vs. 8.1 kg day −1). Heifers in experiment 1 fed ryegrass, compared to the forage mixtures, had greater total tract digestibility (g kg −1) of DM (DMD; P < 0.008, 713 vs. 641) and CP (CPD, P < 0.001, 699 vs. 475), and used more intake energy (%) for body tissue deposition (P < 0.05, 2.6 vs. −4.9). For both experiments, heifers fed flowers differed the most compared to the ryegrass control for a number of measurements. Compared to ryegrass, flowers had 40% lower CP content (P < 0.001, 113 vs. 187 g kg −1), 18% lower DMD (P < 0.01, 585 vs. 713 g kg −1), 42% lower CPD (P < 0.001, 407 vs. 699 g kg −1), and 10% lower methane yield (P < 0.05, 22.6 vs. 25.1 g kg −1 DM intake). This study has shown inclusion of flowers in forage mixtures resulted in a lower CP concentration, digestibility and intake. These differences were due in part to sward management and maturity at harvest. Further research is needed to determine how best to exploit the potential environmental benefits of forage mixtures in sustainable ruminant production systems.

Computational fluid dynamics modelling of the air movement in an environmental test chamber with a respiring manikin

In recent years, computational fluid dynamics (CFD) has been widely used as a method of simulating airflow and addressing indoor environment problems. The complexity of airflows within the indoor environment would make experimental investigation difficult to undertake and also imposes significant challenges on turbulence modelling for flow prediction. This research examines through CFD visualization how air is distributed within a room. Measurements of air temperature and air velocity have been performed at a number of points in an environmental test chamber with a human occupant. To complement the experimental results, CFD simulations were carried out and the results enabled detailed analysis and visualization of spatial distribution of airflow patterns and the effect of different parameters to be predicted. The results demonstrate the complexity of modelling human exhalation within a ventilated enclosure and shed some light into how to achieve more realistic predictions of the airflow within an occupied enclosure.

Rheological and electrical properties of NaY zeolite electrorheological fluid

The relations between the rheological and electrical properties of NaY zeolite electrorheological fluid and its solid phase are studied. It is found that then exist complex relations between its electrical and theological properties. The temperature spectra of dielectric properties of the fluid under high AC electric field are strongly field strength dependent. The relation between the DC conductivity of the fluid and the exciting electric field is experimentally presented as log sigma =A+BE1/2, when A is a strong function, but B, a very weak function of temperature. The shear stress of the fluid under a fixed electric field and temperature decreases with shear rate. A relaxation time for the adsorbed charges is estimated to be about 0.3 to 6.6 s in the temperature range from 280 to 380 K. The relaxation time qualitatively corresponds to the shear rate at which the shear stress begins to drop. The time dependent leaking current of the ER fluids under DC electric field is also measured. The conductivity increase is mainly caused by the structure evolution of particles. The experimental results can he explained with the calculations of Davis (J. Appl. Phys. 81(1997) pp.1985-1991) and Martin (J. Chem. Phys. 110(1999) pp.4854-4866). It is predicted that the NaY zeolite ER fluid strength would get degraded slowly.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

Structure and viscoelasticity of an electrorheological fluid in oscillatory shear: computer simulation investigation

The three-dimensional molecular dynamics simulation method has been used to study the dynamic responses of an electrorheological (ER) fluid in oscillatory shear. The structure and related viscoelastic behaviour of the fluid are found to be sensitive to the amplitude of the strain. With the increase of the strain amplitude, the structure formed by the particles changes from isolated columns to sheet-like structures which may be perpendicular or parallel to the oscillating direction. Along with the structure evolution, the field-induced moduli decrease significantly with an increase in strain amplitude. The viscoelastic behaviour of the structures obtained in the cases of different strain amplitudes was examined in the linear response regime and an evident structure dependence of the moduli was found. The reason for this lies in the anisotropy of the arrangement of the particles in these structures. Short-range interactions between the particles cannot be neglected in determining the viscoelastic behaviour of ER fluids at small strain amplitude, especially for parallel sheets. The simulation results were compared with available experimental data and good agreement was reached for most of them.