Quantifying uncertainty in critical loads: (B) acidity mass balance critical loads on a sensitive site

This paper reports an uncertainty analysis of critical loads for acid deposition for a site in southern England, using the Steady State Mass Balance Model. The uncertainty bounds, distribution type and correlation structure for each of the 18 input parameters was considered explicitly, and overall uncertainty estimated by Monte Carlo methods. Estimates of deposition uncertainty were made from measured data and an atmospheric dispersion model, and hence the uncertainty in exceedance could also be calculated. The uncertainties of the calculated critical loads were generally much lower than those of the input parameters due to a "compensation of errors" mechanism - coefficients of variation ranged from 13% for CLmaxN to 37% for CL(A). With 1990 deposition, the probability that the critical load was exceeded was > 0.99; to reduce this probability to 0.50, a 63% reduction in deposition is required; to 0.05, an 82% reduction. With 1997 deposition, which was lower than that in 1990, exceedance probabilities declined and uncertainties in exceedance narrowed as deposition uncertainty had less effect. The parameters contributing most to the uncertainty in critical loads were weathering rates, base cation uptake rates, and choice of critical chemical value, indicating possible research priorities. However, the different critical load parameters were to some extent sensitive to different input parameters. The application of such probabilistic results to environmental regulation is discussed.

The cumulus-capped boundary layer. I: Modelling transports in the cloud layer

Scalar-flux budgets have been obtained from large-eddy simulations (LESs) of the cumulus-capped boundary layer. Parametrizations of the terms in the budgets are discussed, and two parametrizations for the transport term in the cloud layer are proposed. It is shown that these lead to two models for scalar transports by shallow cumulus convection. One is equivalent to the subsidence detrainment form of convective tendencies obtained from mass-flux parametrizations of cumulus convection. The second is a flux-gradient relationship that is similar in form to the non-local parametrizations of turbulent transports in the dry-convective boundary layer. Using the fluxes of liquid-water potential temperature and total water content from the LES, it is shown that both models are reasonable diagnostic relations between fluxes and the vertical gradients of the mean fields. The LESs used in this study are for steady-state convection and it is possible to treat the fluxes of conserved thermodynamic variables as independent, and ignore the effects of condensation. It is argued that a parametrization of cumulus transports in a model of the cumulus-capped boundary layer should also include an explicit representation of condensation. A simple parametrization of the liquid-water flux in terms of conserved variables is also derived.

Numerical simulation of the 12 May 1997 CME Event: The role of magnetic reconnection

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the 12 May 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance (≈1R ) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

An integrative model of amino acid metabolism in the liver of the lactating dairy cow

The objective of this work was to construct a dynamic model of hepatic amino acid metabolism in the lactating dairy cow that could be parameterized using net flow data from in vivo experiments. The model considers 22 amino acids, ammonia, urea, and 13 energetic metabolites, and was parameterized using a steady-state balance model and two in vivo, net flow experiments conducted with mid-lactation dairy cows. Extracellular flows were derived directly from the observed data. An optimization routine was used to derive nine intracellular flows. The resulting dynamic model was found to be stable across a range of inputs suggesting that it can be perturbed and applied to other physiological states. Although nitrogen was generally in balance, leucine was in slight deficit compared to predicted needs for export protein synthesis, suggesting that an alternative source of leucine (e.g. peptides) was utilized. Simulations of varying glucagon concentrations indicated that an additional 5 mol/d of glucose could be synthesized at the reference substrate concentrations and blood flows. The increased glucose production was supported by increased removal from blood of lactate, glutamate, aspartate, alanine, asparagine, and glutamine. As glucose Output increased, ketone body and acetate release increased while CO2 release declined. The pattern of amino acids appearing in hepatic vein blood was affected by changes in amino acid concentration in portal vein blood, portal blood flow rate and glucagon concentration, with methionine and phenylalanine being the most affected of essential amino acids. Experimental evidence is insufficient to determine whether essential amino acids are affected by varying gluconeogenic demands. (C) 2004 Published by Elsevier Ltd.

Electrochemical determination of plant-derived leuco-indigo after chemical reduction by glucose

Electrochemical determination of redox active dye species is demonstrated in indigo samples contaminated with high levels of organic and inorganic impurities. The use of a hydrodynamic electrode system based on a vibrating probe (250 Hz, 200 mu m lateral amplitude) allows time-independent diffusion controlled signals to be enhanced and reliable concentration data to be obtained under steady state conditions at relatively fast scan rates up to 4 V s-1In this work the indigo content of a complex plant-derived indigo sample (dye content typically 30%) is determined after indigo is reduced by addition of glucose in aqueous 0.2 M NaOH. The soluble leuco-indigo is measured by its oxidation response at a vibrating electrode. The vibrating electrode, which consisted of a laterally vibrating 500 mu m diameter gold disc, is calibrated with Fe(CN)(6) 3-/4- in 0.1 M KCl and employed for indigo determination at 55, 65, and 75 C in 0.2 M NaOH. Determinations of the indigo content of 25 different samples of plant-derived indigo are compared with those obtained by conventional spectrophotometry. This comparison suggests a significant improvement by the electrochemical method, which appears to be less sensitive to impurities.

On the stability of populations of mammals, birds, fish and insects

A key concern for conservation biologists is whether populations of plants and animals are likely to fluctuate widely in number or remain relatively stable around some steady-state value. In our study of 634 populations of mammals, birds, fish and insects, we find that most can be expected to remain stable despite year to year fluctuations caused by environmental factors. Mean return rates were generally around one but were higher in insects (1.09 +/- 0.02 SE) and declined with body size in mammals. In general, this is good news for conservation, as stable populations are less likely to go extinct. However, the lower return rates of the large mammals may make them more vulnerable to extinction. Our estimates of return rates were generally well below the threshold for chaos, which makes it unlikely that chaotic dynamics occur in natural populations - one of ecology's key unanswered questions.

Bone morphogenetic proteins (BMP)-4,-6, and-7 potently suppress basal and luteinizing hormone-induced androgen production by bovine theca interna cells in primary culture: Could ovarian hyperandrogenic dysfunction be caused by a defect in thecal BMP signaling?

We reported recently that bovine theca interna cells in primary culture express several type-I and type-II receptors for bone morphogenetic proteins (BMPs). The same cells express at least two potential ligands for these receptors (BMP-4 and - 7), whereas bovine granulosa cells and oocytes express BMP-6. Therefore, BMPs of intrafollicular origin may exert autocrine/paracrine actions to modulate theca cell function. Here we report that BMP-4, - 6, and - 7 potently suppress both basal ( P < 0.0001; respective IC50 values, 0.78, 0.30, and 1.50 ng/ml) and LH-induced ( P < 0.0001; respective IC50 values, 5.00, 0.55, and 4.55 ng/ml) androgen production by bovine theca cells while having only a moderate effect on progesterone production and cell number. Semiquantitative RT-PCR showed that all three BMPs markedly reduced steady-state levels of mRNA for P450c17. Levels of mRNA encoding steroidogenic acute regulatory protein, P450scc, and 3 beta-hydroxysteroid dehydrogenase were also reduced but to a much lesser extent. Immunocytochemistry confirmed a marked reduction in cellular content of P450c17 protein after BMP treatment ( P < 0.001). Exposure to BMPs led to cellular accumulation of phosphorylated Smad1, but not Smad2, confirming that the receptors signal via a Smad1 pathway. The specificity of the BMP response was further explored by coincubating cells with BMPs and several potential BMP antagonists, chordin, gremlin, and follistatin. Gremlin and chordin were found to be effective antagonists of BMP-4 and - 7, respectively, and the observation that both antagonists enhanced ( P < 0.01) androgen production in the absence of exogenous BMP suggests an autocrine/paracrine role for theca-derived BMP- 4 and - 7 in modulating androgen production. Collectively, these data indicate that an intrafollicular BMP signaling pathway contributes to the negative regulation of thecal androgen production and that ovarian hyperandrogenic dysfunction could be a result of a defective autoregulatory pathway involving thecal BMP signaling.

Aldose reductase and the importance of experimental design

Kinetic studies on the AR (aldose reductase) protein have shown that it does not behave as a classical enzyme in relation to ring aldose sugars. As with non-enzymatic glycation reactions, there is probably a free radical element involved derived from monosaccharide autoxidation. in the case of AR, there is free radical oxidation of NADPH by autoxidizing monosaccharides, which is enhanced in the presence of the NADPH-binding protein. Thus any assay for AR based on the oxidation of NADPH in the presence of autoxidizing monosaccharides is invalid, and tissue AR measurements based on this method are also invalid, and should be reassessed. AR exhibits broad specificity for both hydrophilic and hydrophobic aldehydes that suggests that the protein may be involved in detoxification. The last thing we would want to do is to inhibit it. ARIs (AR inhibitors) have a number of actions in the cell which are not specific, and which do not involve them binding to AR. These include peroxy-radical scavenging and effects of metal ion chelation. The AR/ARI story emphasizes the importance of correct experimental design in all biocatalytic experiments. Developing the use of Bayesian utility functions, we have used a systematic method to identify the optimum experimental designs for a number of kinetic model data sets. This has led to the identification of trends between kinetic model types, sets of design rules and the key conclusion that such designs should be based on some prior knowledge of K-m and/or the kinetic model. We suggest an optimal and iterative method for selecting features of the design such as the substrate range, number of measurements and choice of intermediate points. The final design collects data suitable for accurate modelling and analysis and minimizes the error in the parameters estimated, and is suitable for simple or complex steady-state models.

Direct detection of dimethylstannylene and tetramethyldistannene in solution and the gas phase by laser flash photolysis of 1,1-dimethylstannacyclopent-3-enes

The photochemistry of 1,1-dimethyl- and 1,1,3,4-tetramethylstannacyclopent-3-ene (4a and 4b,respectively) has been studied in the gas phase and in hexane solution by steady-state and 193-nm laser flash photolysis methods. Photolysis of the two compounds results in the formation of 1,3-butadiene (from 4a) and 2,3-dimethyl-1,3-butadiene (from 4b) as the major products, suggesting that cycloreversion to yield dimethylstannylene (SnMe2) is the main photodecomposition pathway of these molecules. Indeed, the stannylene has been trapped as the Sn-H insertion product upon photolysis of 4a in hexane containing trimethylstannane. Flash photolysis of 4a in the gas phase affords a transient absorbing in the 450-520nm range that is assigned to SnMe2 by comparison of its spectrum and reactivity to those previously reported from other precursors. Flash photolysis of 4b in hexane solution affords results consistent with the initial formation of SnMe2 (lambda(max) approximate to 500 nm), which decays over similar to 10 mu s to form tetramethyldistannene (5b; lambda(max) approximate to 470 nm). The distannene decays over the next ca. 50 mu s to form at least two other longer-lived species, which are assigned to higher SnMe2 oligomers. Time-dependent DFT calculations support the spectral assignments for SnMe2 and Sn2Me4, and calculations examining the variation in bond dissociation energy with substituent (H, Me, and Ph) in disilenes, digermenes, and distannenes rule out the possibility that dimerization of SnMe2 proceeds reversibly. Addition of methanol leads to reversible reaction with SnMe2 to form a transient absorbing at lambda(max) approximate to 360 nm, which is assigned to the Lewis acid-base complex between SnMe2 and the alcohol.

Toughness & cutting: A new way of determining ductile fracture toughness and strength

A review is given of the mechanics of cutting, ranging from the slicing of thin floppy offcuts (where there is negligible elasticity and no permanent deformation of the offcut) to the machining of ductile metals (where there is severe permanent distortion of the offcut/chip). Materials scientists employ the former conditions to determine the fracture toughness of ‘soft’ solids such as biological materials and foodstuffs. In contrast, traditional analyses of metalcutting are based on plasticity and friction only, and do not incorporate toughness. The machining theories are inadequate in a number of ways but a recent paper has shown that when ductile work of fracture is included many, if not all, of the shortcomings are removed. Support for the new analysis is given by examination of FEM simulations of metalcutting which reveal that a ‘separation criterion’ has to be employed at the tool tip. Some consideration shows that the separation criteria are versions of void-initiation-growth-and-coalescence models employed in ductile fracture mechanics. The new analysis shows that cutting forces for ductile materials depend upon the fracture toughness as well as plasticity and friction, and reveals a simple way of determining both toughness and flow stress from cutting experiments. Examples are given for a wide range of materials including metals, polymers and wood, and comparison is made with the same properties independently determined using conventional testpieces. Because cutting can be steady state, a new way is presented for simultaneously measuring toughness and flow stress at controlled speeds and strain rates.

The use of wireless data communication and body sensing devices to evaluate occupants' comfort in buildings

Physiological parameters measured by an embedded body sensor system were demonstrated to respond to changes of the air temperature in an office environment. The thermal parameters were monitored with the use of a wireless sensor system that made possible to turn any existing room into a field laboratory. Two human subjects were monitored over daily activities and at various steady-state thermal conditions when the air temperature of the room was altered from 22-23°C to 25-28°C. The subjects indicated their thermal feeling on questionnaires. The measured skin temperature was distributed close to the calculated mean skin temperature corresponding to the given activity level. The variation of Galvanic Skin Response (GSR) reflected the evaporative heat loss through the body surfaces and indicated whether sweating occurred on the subjects. Further investigations are needed to fully evaluate the influence of thermal and other factors on the output given by the investigated body sensor system.

Stress relaxation behaviour of frozen sucrose solutions

The stress relaxation behaviour of two frozen sucrose solutions (7% and 19%) during indentation in the temperature range of -20C to -40C were investigated. The stress relaxation is similar to that of pure polycrystalline ice, which is controlled by steady-state creep. The steady state creep rate exponent, m, of 7% and 19% sucrose solutions lies between 2.3 and 3.6. The steady state creep rate constant, B, of 19% sucrose solution is greater than that of 7% sucrose solution. It is suggested that the steady-state creep rate exponent m depends on contributions from the proportions of favourably oriented grains, unfavourably oriented grains and grain boundaries to creep and that these components depend on the value of internal stress which is related to the hardness of samples at the different testing temperatures. The steady-state creep rate constant B depends on the mobility of dislocations in sucrose solutions which, in turn, depends on the temperature and the concentration of sucrose.

Comparison of small and large deformation rheological properties of wheat dough and gluten

The rheological properties of dough and gluten are important for end-use quality of flour but there is a lack of knowledge of the relationships between fundamental and empirical tests and how they relate to flour composition and gluten quality. Dough and gluten from six breadmaking wheat qualities were subjected to a range of rheological tests. Fundamental (small-deformation) rheological characterizations (dynamic oscillatory shear and creep recovery) were performed on gluten to avoid the nonlinear influence of the starch component, whereas large deformation tests were conducted on both dough and gluten. A number of variables from the various curves were considered and subjected to a principal component analysis (PCA) to get an overview of relationships between the various variables. The first component represented variability in protein quality, associated with elasticity and tenacity in large deformation (large positive loadings for resistance to extension and initial slope of dough and gluten extension curves recorded by the SMS/Kieffer dough and gluten extensibility rig, and the tenacity and strain hardening index of dough measured by the Dobraszczyk/Roberts dough inflation system), the elastic character of the hydrated gluten proteins (large positive loading for elastic modulus [G'], large negative loadings for tan delta and steady state compliance [J(e)(0)]), the presence of high molecular weight glutenin subunits (HMW-GS) 5+10 vs. 2+12, and a size distribution of glutenin polymers shifted toward the high-end range. The second principal component was associated with flour protein content. Certain rheological data were influenced by protein content in addition to protein quality (area under dough extension curves and dough inflation curves [W]). The approach made it possible to bridge the gap between fundamental rheological properties, empirical measurements of physical properties, protein composition, and size distribution. The interpretation of this study gave indications of the molecular basis for differences in breadmaking performance.

Numerical study of 2D heat transfer in a scraped surface heat exchanger

A numerical study of fluid mechanics and heat transfer in a scraped surface heat exchanger with non-Newtonian power law fluids is undertaken. Numerical results are generated for 2D steady-state conditions using finite element methods. The effect of blade design and material properties, and especially the independent effects of shear thinning and heat thinning on the flow and heat transfer, are studied. The results show that the gaps at the root of the blades, where the blades are connected to the inner cylinder, remove the stagnation points, reduce the net force on the blades and shift the location of the central stagnation point. The shear thinning property of the fluid reduces the local viscous dissipation close to the singularity corners, i.e. near the tip of the blades, and as a result the local fluid temperature is regulated. The heat thinning effect is greatest for Newtonian fluids where the viscous dissipation and the local temperature are highest at the tip of the blades. Where comparison is possible, very good agreement is found between the numerical results and the available data. Aspects of scraped surface heat exchanger design are assessed in the light of the results. (C) 2003 Elsevier Ltd. All rights reserved.

Biodiversity of human faecal bacteria isolated from phytic acid enriched chemostat fermenters

Background: Myo-inositol hexaphosphate (IP6) or phytic acid is found mostly in cereals and legumes and is thought to possess anti-carcinogenic properties. Aim: To isolate and identify faecal bacteria capable of phytic acid metabolism and to assess the effectiveness of prebiotics (dietary oligosaccharides, metabolised by selective colonic bacteria) in preserving the integrity of phytic acid. Methods: Faecal samples from three volunteers were used in continuous culture experiments under varying conditions of pH, substrate concentration and dilution rates, seventy three different isolates cultured at steady state were then screened for phytic acid metabolism and identified through partial sequencing of their 16S rRNA genes (16S ribosomal ribonucleic acid). Utilisation of phytic acid was also assessed in a continuous culture system enriched with prebiotic fructooligosaccharides (FOS). Results: Bacteroides spp., Clostridium spp. and facultatively anaerobic bacteria generally appeared to maintain viable counts in the presence of phytic acid. Bifidobacterium spp. and Lactobacillus spp. appeared less able to maintain viable counts in the presence of phytic acid. These results were confirmed by an increase in viable counts of Bacteroides spp., Clostridium spp. and a decrease in viable counts of Bifidobacterium spp. and Lactobacillus spp. once phytic acid was introduced to a FOS enriched continuous culture. Conclusions: The phytate metabolising biodiversity from the human large intestine does not appear to encompass major bacterial genera associated with beneficial or benign health effects (e.g. Lactobacillus spp. and Bifidobacterium spp).