A simple column model to explore anticipated problems in variational assimilation of satellite observations

We investigate a simplified form of variational data assimilation in a fully nonlinear framework with the aim of extracting dynamical development information from a sequence of observations over time. Information on the vertical wind profile, w(z ), and profiles of temperature, T (z , t), and total water content, qt (z , t), as functions of height, z , and time, t, are converted to brightness temperatures at a single horizontal location by defining a two-dimensional (vertical and time) variational assimilation testbed. The profiles of T and qt are updated using a vertical advection scheme. A basic cloud scheme is used to obtain the fractional cloud amount and, when combined with the temperature field, this information is converted into a brightness temperature, using a simple radiative transfer scheme. It is shown that our model exhibits realistic behaviour with regard to the prediction of cloud, but the effects of nonlinearity become non-negligible in the variational data assimilation algorithm. A careful analysis of the application of the data assimilation scheme to this nonlinear problem is presented, the salient difficulties are highlighted, and suggestions for further developments are discussed.

A metabolic system-wide characterisation of the pig: a model for human physiology

The pig is a single-stomached omnivorous mammal and is an important model of human disease and nutrition. As such, it is necessary to establish a metabolic framework from which pathology-based variation can be compared. Here, a combination of one and two-dimensional (1)H and (13)C nuclear magnetic resonance spectroscopy (NMR) and high-resolution magic angle spinning (HR-MAS) NMR was used to provide a systems overview of porcine metabolism via characterisation of the urine, serum, liver and kidney metabolomes. The metabolites observed in each of these biological compartments were found to be qualitatively comparable to the metabolic signature of the same biological matrices in humans and rodents. The data were modelled using a combination of principal components analysis and Venn diagram mapping. Urine represented the most metabolically distinct biological compartment studied, with a relatively greater number of NMR detectable metabolites present, many of which are implicated in gut-microbial co-metabolic processes. The major inter-species differences observed were in the phase II conjugation of extra-genomic metabolites; the pig was observed to conjugate p-cresol, a gut microbial metabolite of tyrosine, with glucuronide rather than sulfate as seen in man. These observations are important to note when considering the translatability of experimental data derived from porcine models.

Identification of metabolites in human hepatic bile using 800 MHz 1H NMR spectroscopy, HPLC-NMR/MS and UPLC-MS

The first application of high field NMR spectroscopy (800 MHz for 1H observation) to human hepatic bile (as opposed to gall bladder bile) is reported. The bile sample used for detailed investigation was from a donor liver with mild fat infiltration, collected during organ retrieval prior to transplantation. In addition, to focus on the detection of bile acids in particular, a bile extract was analysed by 800 MHz 1H NMR spectroscopy, HPLC-NMR/MS and UPLC-MS. In the whole bile sample, 40 compounds have been assigned with the aid of two-dimensional 1H–1H TOCSY and 1H–13C HSQC spectra. These include phosphatidylcholine, 14 amino acids, 10 organic acids, 4 carbohydrates and polyols (glucose, glucuronate, glycerol and myo-inositol), choline, phosphocholine, betaine, trimethylamine-N-oxide and other small molecules. An initial NMR-based assessment of the concentration range of some key metabolites has been made. Some observed chemical shifts differ from expected database values, probably due to a difference in bulk diamagnetic susceptibility. The NMR spectra of the whole extract gave identification of the major bile acids (cholic, deoxycholic and chenodeoxycholic), but the glycine and taurine conjugates of a given bile acid could not be distinguished. However, this was achieved by HPLC-NMR/MS, which enabled the separation and identification of ten conjugated bile acids with relative abundances varying from approximately 0.1% (taurolithocholic acid) to 34.0% (glycocholic acid), of which, only the five most abundant acids could be detected by NMR, including the isomers glycodeoxycholic acid and glycochenodeoxycholic acid, which are difficult to distinguish by conventional LC-MS analysis. In a separate experiment, the use of UPLC-MS allowed the detection and identification of 13 bile acids. This work has shown the complementary potential of NMR spectroscopy, MS and hyphenated NMR/MS for elucidating the complex metabolic profile of human hepatic bile. This will be useful baseline information in ongoing studies of liver excretory function and organ transplantation.

Symmetry breaking, mixing, instability, and low frequency variability in a minimal Lorenz-like system

Starting from the classical Saltzman two-dimensional convection equations, we derive via a severe spectral truncation a minimal 10 ODE system which includes the thermal effect of viscous dissipation. Neglecting this process leads to a dynamical system which includes a decoupled generalized Lorenz system. The consideration of this process breaks an important symmetry and couples the dynamics of fast and slow variables, with the ensuing modifications to the structural properties of the attractor and of the spectral features. When the relevant nondimensional number (Eckert number Ec) is different from zero, an additional time scale of O(Ec−1) is introduced in the system, as shown with standard multiscale analysis and made clear by several numerical evidences. Moreover, the system is ergodic and hyperbolic, the slow variables feature long-term memory with 1/f3/2 power spectra, and the fast variables feature amplitude modulation. Increasing the strength of the thermal-viscous feedback has a stabilizing effect, as both the metric entropy and the Kaplan-Yorke attractor dimension decrease monotonically with Ec. The analyzed system features very rich dynamics: it overcomes some of the limitations of the Lorenz system and might have prototypical value in relevant processes in complex systems dynamics, such as the interaction between slow and fast variables, the presence of long-term memory, and the associated extreme value statistics. This analysis shows how neglecting the coupling of slow and fast variables only on the basis of scale analysis can be catastrophic. In fact, this leads to spurious invariances that affect essential dynamical properties (ergodicity, hyperbolicity) and that cause the model losing ability in describing intrinsically multiscale processes.

Identification of components of the Sigma B Regulon in Listeria monocytogenes that contribute to acid and salt tolerance

Sigma B (σB) is an alternative sigma factor that controls the transcriptional response to stress in Listeria monocytogenes and is also known to play a role in the virulence of this human pathogen. In the present study we investigated the impact of a sigB deletion on the proteome of L. monocytogenes grown in a chemically defined medium both in the presence and in the absence of osmotic stress (0.5 M NaCl). Two new phenotypes associated with the sigB deletion were identified using this medium. (i) Unexpectedly, the strain with the ΔsigB deletion was found to grow faster than the parent strain in the growth medium, but only when 0.5 M NaCl was present. This phenomenon was independent of the carbon source provided in the medium. (ii) The ΔsigB mutant was found to have unusual Gram staining properties compared to the parent, suggesting that σB contributes to the maintenance of an intact cell wall. A proteomic analysis was performed by two-dimensional gel electrophoresis, using cells growing in the exponential and stationary phases. Overall, 11 proteins were found to be differentially expressed in the wild type and the ΔsigB mutant; 10 of these proteins were expressed at lower levels in the mutant, and 1 was overexpressed in the mutant. All 11 proteins were identified by tandem mass spectrometry, and putative functions were assigned based on homology to proteins from other bacteria. Five proteins had putative functions related to carbon utilization (Lmo0539, Lmo0783, Lmo0913, Lmo1830, and Lmo2696), while three proteins were similar to proteins whose functions are unknown but that are known to be stress inducible (Lmo0796, Lmo2391, and Lmo2748). To gain further insight into the role of σB in L. monocytogenes, we deleted the genes encoding four of the proteins, lmo0796, lmo0913, lmo2391, and lmo2748. Phenotypic characterization of the mutants revealed that Lmo2748 plays a role in osmotolerance, while Lmo0796, Lmo0913, and Lmo2391 were all implicated in acid stress tolerance to various degrees. Invasion assays performed with Caco-2 cells indicated that none of the four genes was required for mammalian cell invasion. Microscopic analysis suggested that loss of Lmo2748 might contribute to the cell wall defect observed in the ΔsigB mutant. Overall, this study highlighted two new phenotypes associated with the loss of σB. It also demonstrated clear roles for σB in both osmotic and low-pH stress tolerance and identified specific components of the σB regulon that contribute to the responses observed.

Proteomic analyses of a Listeria monocytogenes mutant lacking σB identify new components of the σB regulon and highlight a role for σB in the utilization of glycerol

In Listeria monocytogenes the alternative sigma factor σB plays important roles in both virulence and stress tolerance. In this study a proteomic approach was used to define components of the σB regulon in L. monocytogenes 10403S (serotype 1/2a). Using two-dimensional gel electrophoresis and the recently developed isobaric tags for relative and absolute quantitation technique, the protein expression profiles of the wild type and an isogenic ΔsigB deletion strain were compared. Overall, this study identified 38 proteins whose expression was σB dependent; 17 of these proteins were found to require the presence of σB for full expression, while 21 were expressed at a higher level in the ΔsigB mutant background. The data obtained with the two proteomic approaches showed limited overlap (four proteins were identified by both methods), a finding that highlights the complementarity of the two technologies. Overall, the proteomic data reaffirmed a role for σB in the general stress response and highlighted a probable role for σB in metabolism, especially in the utilization of alternative carbon sources. Proteomic and physiological data revealed the involvement of σB in glycerol metabolism. Five newly identified members of the σB regulon were shown to be under direct regulation of σB using reverse transcription-PCR (RT-PCR), while random amplification of cDNA ends-PCR was used to map four σB-dependent promoters upstream from lmo0796, lmo1830, lmo2391, and lmo2695. Using RT-PCR analysis of known and newly identified σB-dependent genes, as well as proteomic analyses, σB was shown to play a major role in the stationary phase of growth in complex media.

Observational tracking of the 2D structure of coronal mass ejections between the sun and 1 AU

The Solar TErrestrial RElations Observatory (STEREO) provides high cadence and high resolution images of the structure and morphology of coronal mass ejections (CMEs) in the inner heliosphere. CME directions and propagation speeds have often been estimated through the use of time-elongation maps obtained from the STEREO Heliospheric Imager (HI) data. Many of these CMEs have been identified by citizen scientists working within the SolarStormWatch project ( www.solarstormwatch.com ) as they work towards providing robust real-time identification of Earth-directed CMEs. The wide field of view of HI allows scientists to directly observe the two-dimensional (2D) structures, while the relative simplicity of time-elongation analysis means that it can be easily applied to many such events, thereby enabling a much deeper understanding of how CMEs evolve between the Sun and the Earth. For events with certain orientations, both the rear and front edges of the CME can be monitored at varying heliocentric distances (R) between the Sun and 1 AU. Here we take four example events with measurable position angle widths and identified by the citizen scientists. These events were chosen for the clarity of their structure within the HI cameras and their long track lengths in the time-elongation maps. We show a linear dependency with R for the growth of the radial width (W) and the 2D aspect ratio (χ) of these CMEs, which are measured out to ≈ 0.7 AU. We estimated the radial width from a linear best fit for the average of the four CMEs. We obtained the relationships W=0.14R+0.04 for the width and χ=2.5R+0.86 for the aspect ratio (W and R in units of AU).

Simulation of a domestic ground source heat pump system using a three-dimensional numerical borehole heat exchanger model

Common approaches to the simulation of borehole heat exchangers (BHEs) assume heat transfer in circulating fluid and grout to be in a quasi-steady state and ignore fluctuations in fluid temperature due to transport of the fluid around the loop. However, in domestic ground source heat pump (GSHP) systems, the heat pump and circulating pumps switch on and off during a given hour; therefore, the effect of the thermal mass of the circulating fluid and the dynamics of fluid transport through the loop has important implications for system design. This may also be important in commercial systems that are used intermittently. This article presents transient simulation of a domestic GSHP system with a single BHE using a dynamic three-dimensional (3D) numerical BHE model. The results show that delayed response associated with the transit of fluid along the pipe loop is of some significance in moderating swings in temperature during heat pump operation. In addition, when 3D effects are considered, a lower heat transfer rate is predicted during steady operations. These effects could be important when considering heat exchanger design and system control. The results will be used to develop refined two-dimensional models.

Effect of fluoroquinolone exposure on the proteome of Salmonella enterica serovar Typhimurium

Objectives: The physiological response of Salmonella enterica serovar Typhimurium to fluoroquinolone antibiotics was investigated using proteomic methods. Methods: Proteomes were prepared from strain SL1344 following treatment of broth cultures with ciprofloxacin (0.03 and 0.008 mg/L; 2x and 0.5x MIC) and enrofloxacin (0.03 mg/L) and from a multiple antibiotic resistant (MAR) mutant. Protein expression was determined by two-dimensional HPLC-MSn and also after exposure to ciprofloxacin by two-dimensional gel electrophoresis (2D-GE). Results: The number of proteins (mean +/- SD) detected by 2D-GE derived from control cultures of the wild-type strain was significantly (P < 0.05) reduced from 296 +/- 77 to 153 +/- 36 following treatment with ciprofloxacin (0.03 mg/L). Raised expression (P < 0.05) of 17 proteins was also detected, and increases of up to 8-fold (P < 0.0001) were observed for subunits of F1F0-ATP synthase, TolC and Imp. Analysis by two-dimensional HPLC-MSn provided higher proteome coverage with 787 +/- 50 proteins detected, which was reduced (P < 0.005) to 560 +/- 14 by ciprofloxacin (0.03 mg/L). Increased expression of 43 proteins was observed which included those detected by 2D-GE and additionally the efflux pump protein AcrB. The basal expression of the AcrAB/TolC efflux pump was elevated in the MAR mutant compared with the untreated wild-type and augmented following treatment with ciprofloxacin (0.03 mg/L). F1F0-ATP synthase and Imp were only elevated in the mutant when treated with ciprofloxacin. Conclusions: These studies suggest that increased expression of AcrAB/TolC was associated with resistance while other increases, such as in F1F0-ATP synthase and Imp, were a response to fluoroquinolone.

Commonly used farm disinfectants can select for mutant Salmonella enterica serovar Typhimurium with decreased susceptibility to biocides and antibiotics without compromising virulence

Objectives: To determine if one passage of Salmonella enterica serovar Typhimurium in the presence of farm disinfectants selected for mutants with decreased susceptibility to disinfectants and/or antibiotics. Methods: Eight Salmonella Typhimurium strains including field isolates and laboratory mutants were exposed to either a tar oil phenol (PFD) disinfectant, an oxidizing compound disinfectant (OXC), an aldehyde based disinfectant (ABD) or a dairy sterilizer disinfectant (based on quaternary ammonium biocide) in agar. The susceptibility of mutants obtained after disinfectant exposure to antibiotics and disinfectants was determined as was the accumulation of norfloxacin. The proteome of SL1344 after exposure to PFD and OXC was analysed using two-dimensional liquid chromatography mass spectrometry. Results: Strains with either acrB or tolC inactivated were more susceptible to most disinfectants than other strains. The majority (3/5) of mutants recovered after disinfectant exposure required statistically significantly longer exposure times to disinfectants than their parent strains to generate a 5 log kill. Small decreases in antibiotic susceptibility were observed but no mutants were multiply antibiotic-resistant (MAR). Notably exposure to ABD decreased susceptibility to ciprofloxacin in some strains. Mutants with increased disinfectant tolerance were able to survive and persist in chicks as well as in parent strains. Analysis of proteomes revealed significantly increased expression of the AcrAB-TolC efflux system after PFD exposure. Conclusions: Data presented demonstrate that efflux pumps are required for intrinsic resistance to some disinfectants and that exposure to disinfectants can induce expression of the AcrAB-TolC efflux system, but that single exposure was insufficient to select for MAR strains.

The unsteady flow of a weakly compressible fluid in a thin porous layer II: three-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a three-dimensional layer, composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Numerical solution of this three-dimensional evolution problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l, a situation which occurs frequently in the application to oil and gas reservoir recovery and which leads to significant stiffness in the numerical problem. Under the assumption that $\epsilon\propto h/l\ll 1$, we show that, to leading order in $\epsilon$, the pressure field varies only in the horizontal directions away from the wells (the outer region). We construct asymptotic expansions in $\epsilon$ in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive expressions for all significant process quantities. The only computations required are for the solution of non-stiff linear, elliptic, two-dimensional boundary-value, and eigenvalue problems. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the layer, $\epsilon$, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighbourhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer III: three-dimensional computations

We describe a novel method for determining the pressure and velocity fields for a weakly compressible fluid flowing in a thin three-dimensional layer composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Our approach uses the method of matched asymptotic expansions to derive expressions for all significant process quantities, the computation of which requires only the solution of linear, elliptic, two-dimensional boundary value and eigenvalue problems. In this article, we provide full implementation details and present numerical results demonstrating the efficiency and accuracy of our scheme.

Isomap nonlinear dimensionality reduction and bimodality of Asian monsoon convection

It is known that the empirical orthogonal function method is unable to detect possible nonlinear structure in climate data. Here, isometric feature mapping (Isomap), as a tool for nonlinear dimensionality reduction, is applied to 1958–2001 ERA-40 sea-level pressure anomalies to study nonlinearity of the Asian summer monsoon intraseasonal variability. Using the leading two Isomap time series, the probability density function is shown to be bimodal. A two-dimensional bivariate Gaussian mixture model is then applied to identify the monsoon phases, the obtained regimes representing enhanced and suppressed phases, respectively. The relationship with the large-scale seasonal mean monsoon indicates that the frequency of monsoon regime occurrence is significantly perturbed in agreement with conceptual ideas, with preference for enhanced convection on intraseasonal time scales during large-scale strong monsoons. Trend analysis suggests a shift in concentration of monsoon convection, with less emphasis on South Asia and more on the East China Sea.