A singular vector perspective of 4D-Var: Filtering and interpolation

Four-dimensional variational data assimilation (4D-Var) combines the information from a time sequence of observations with the model dynamics and a background state to produce an analysis. In this paper, a new mathematical insight into the behaviour of 4D-Var is gained from an extension of concepts that are used to assess the qualitative information content of observations in satellite retrievals. It is shown that the 4D-Var analysis increments can be written as a linear combination of the singular vectors of a matrix which is a function of both the observational and the forecast model systems. This formulation is used to consider the filtering and interpolating aspects of 4D-Var using idealized case-studies based on a simple model of baroclinic instability. The results of the 4D-Var case-studies exhibit the reconstruction of the state in unobserved regions as a consequence of the interpolation of observations through time. The results also exhibit the filtering of components with small spatial scales that correspond to noise, and the filtering of structures in unobserved regions. The singular vector perspective gives a very clear view of this filtering and interpolating by the 4D-Var algorithm and shows that the appropriate specification of the a priori statistics is vital to extract the largest possible amount of useful information from the observations. Copyright © 2005 Royal Meteorological Society

Comparison of energy evaluation systems and a mechanistic model for milk production by dairy cattle offered fresh grass-based diets

Grass-based diets are of increasing social-economic importance in dairy cattle farming, but their low supply of glucogenic nutrients may limit the production of milk. Current evaluation systems that assess the energy supply and requirements are based on metabolisable energy (ME) or net energy (NE). These systems do not consider the characteristics of the energy delivering nutrients. In contrast, mechanistic models take into account the site of digestion, the type of nutrient absorbed and the type of nutrient required for production of milk constituents, and may therefore give a better prediction of supply and requirement of nutrients. The objective of the present study is to compare the ability of three energy evaluation systems, viz. the Dutch NE system, the agricultural and food research council (AFRC) ME system, and the feed into milk (FIM) ME system, and of a mechanistic model based on Dijkstra et al. [Simulation of digestion in cattle fed sugar cane: prediction of nutrient supply for milk production with locally available supplements. J. Agric. Sci., Cambridge 127, 247-60] and Mills et al. [A mechanistic model of whole-tract digestion and methanogenesis in the lactating dairy cow: model development, evaluation and application. J. Anim. Sci. 79, 1584-97] to predict the feed value of grass-based diets for milk production. The dataset for evaluation consists of 41 treatments of grass-based diets (at least 0.75 g ryegrass/g diet on DM basis). For each model, the predicted energy or nutrient supply, based on observed intake, was compared with predicted requirement based on observed performance. Assessment of the error of energy or nutrient supply relative to requirement is made by calculation of mean square prediction error (MSPE) and by concordance correlation coefficient (CCC). All energy evaluation systems predicted energy requirement to be lower (6-11%) than energy supply. The root MSPE (expressed as a proportion of the supply) was lowest for the mechanistic model (0.061), followed by the Dutch NE system (0.082), FIM ME system (0.097) and AFRCME system(0.118). For the energy evaluation systems, the error due to overall bias of prediction dominated the MSPE, whereas for the mechanistic model, proportionally 0.76 of MSPE was due to random variation. CCC analysis confirmed the higher accuracy and precision of the mechanistic model compared with energy evaluation systems. The error of prediction was positively related to grass protein content for the Dutch NE system, and was also positively related to grass DMI level for all models. In conclusion, current energy evaluation systems overestimate energy supply relative to energy requirement on grass-based diets for dairy cattle. The mechanistic model predicted glucogenic nutrients to limit performance of dairy cattle on grass-based diets, and proved to be more accurate and precise than the energy systems. The mechanistic model could be improved by allowing glucose maintenance and utilization requirements parameters to be variable. (C) 2007 Elsevier B.V. All rights reserved.

An evaluation of a Bayesian method of dose escalation based on bivariate binary responses

Recently, various approaches have been suggested for dose escalation studies based on observations of both undesirable events and evidence of therapeutic benefit. This article concerns a Bayesian approach to dose escalation that requires the user to make numerous design decisions relating to the number of doses to make available, the choice of the prior distribution, the imposition of safety constraints and stopping rules, and the criteria by which the design is to be optimized. Results are presented of a substantial simulation study conducted to investigate the influence of some of these factors on the safety and the accuracy of the procedure with a view toward providing general guidance for investigators conducting such studies. The Bayesian procedures evaluated use logistic regression to model the two responses, which are both assumed to be binary. The simulation study is based on features of a recently completed study of a compound with potential benefit to patients suffering from inflammatory diseases of the lung.

Self-assembly of peptide nanotubes in an organic solvent

The self-assembly of a modified fragment of the amyloid beta peptide, based on sequence A beta(16-20), KLVFF, extended to give AAKLVFF is studied in methanol. Self-assembly into peptide nanotubes is observed, as confirmed by electron microscopy and small-angle X-ray scattering. The secondary structure of the peptide is probed by FTIR and circular dichroism, and UV/visible spectroscopy provides evidence for the important role of aromatic interactions between phenylalanine residues in driving beta-sheet self-assembly. The beta-sheets wrap helically to form the nanotubes, the nanotube wall comprising four wrapped beta-sheets. At higher concentration, the peptide nanotubes form a nematic phase that exhibits spontaneous flow alignment as observed by small-angle neutron scattering.

Characterisation of blends based on hydroxyethylcellulose and maleic acid-alt-methyl vinyl ether

Hydrophilic polymeric films based on blends of hydroxyethylcellulose and maleic acid-co-methyl vinyl ether were produced by casting from aqueous solutions. The physicochemical properties of the blends have been assessed using Fourier transform infrared spectroscopy, thermal gravimetric analysis, differential scanning calorimetry, dielectric spectroscopy, etc. The pristine films exhibit complete miscibility due to the formation of intermacromolecular hydrogen bonding. The thermal treatment of the blend films leads to cross-linking via intermacromolecular esterification and anhydride formation. The cross-linked materials are able to swell in water and their swelling degree can be easily controlled by temperature and thermal treatment time. The formation of the crosslinks is apparent in the dynamic properties of the blends as observed through the mechanical relaxation and dielectric relaxation spectra. The dielectric characteristics of the material are influenced by the effects of change in the local structure of the blend on the ionic conduction processes and the rate of dipolar relaxation. Separation of these processes is attempted using the dielectric modulus method. Significant deviations from a simple additive rule of mixing on the activation energy are observed consistent with hydrogen bonding and crosslinking of the matrix. This paper indicates a method for the creation of films with good mechanical and physical characteristics by exposing the blends to a relatively mild thermal treatment.

Alignment of a model amyloid peptide fragment in bulk and at a solid surface

The alignment of model amyloid peptide YYKLVFFC is investigated in bulk and at a solid surface using a range of spectroscopic methods employing polarized radiation. The peptide is based on a core sequence of the amyloid beta (A beta) peptide, KLVFF. The attached tyrosine and cysteine units are exploited to yield information on alignment and possible formation of disulfide or dityrosine links. Polarized Raman spectroscopy on aligned stalks provides information on tyrosine orientation, which complements data from linear dichroism (LD) on aqueous solutions subjected to shear in a Couette cell. LD provides a detailed picture of alignment of peptide strands and aromatic residues and was also used to probe the kinetics of self-assembly. This suggests initial association of phenylalanine residues, followed by subsequent registry of strands and orientation of tyrosine residues. X-ray diffraction (XRD) data from aligned stalks is used to extract orientational order parameters from the 0.48 nm reflection in the cross-beta pattern, from which an orientational distribution function is obtained. X-ray diffraction on solutions subject to capillary flow confirmed orientation in situ at the level of the cross-beta pattern. The information on fibril and tyrosine orientation from polarized Raman spectroscopy is compared with results from NEXAFS experiments on samples prepared as films on silicon. This indicates fibrils are aligned parallel to the surface, with phenyl ring normals perpendicular to the surface. Possible disulfide bridging leading to peptide dimer formation was excluded by Raman spectroscopy, whereas dityrosine formation was probed by fluorescence experiments and was found not to occur except under alkaline conditions. Congo red binding was found not to influence the cross-beta XRD pattern.

Self-assembly of amphiphilic peptides

The self-assembly of amphiphilic peptides is reviewed. The review covers surfactant-like peptides with amphiphilicity arising from the sequence of natural amino acids, and also peptide amphiphiles (PAs) in which lipid chains are attached to hydrophilic peptide sequences containing charged residues. The influence of the secondary structure on the self-assembled structure and vice versa is discussed. For surfactant-like peptides structures including fibrils, nanotubes, micelles and vesicles have been reported. A particularly common motif for PAs is beta-sheet based fibrils, although other structures have been observed. In these structures, the peptide epitope is presented at the surface of the nanostructure, providing remarkable bioactivity. Recent discoveries of potential, and actual, applications of these materials in biomedicine and bionanotechnology are discussed.

Effect of linseed oil and fish oil alone or as an equal mixture of ruminal fatty acid metabolism in growing steers fed maize silage-based diets

Based on the potential benefits to human health there is interest in increasing 18:3n-3, 20:5n-3, 22:6n-6, and cis-9,trans-11 conjugated linoleic acid (CLA) in ruminant foods. Four Aberdeen Angus steers (406 ± 8.2 kg BW) fitted with rumen and duodenal cannulae were used in a 4 x 4 Latin square experiment with 21 d periods to examine the potential of fish oil (FO) and linseed oil (LO) in the diet to increase ruminal outflow of trans-11 18:1 and total n-3 polyunsaturated fatty acids (PUFA) in growing cattle. Treatments consisted of a control diet (60:40; forage:concentrate ratio, on a DM basis, respectively) based on maize silage, or the same basal ration containing 30 g/kg DM of FO, LO or a mixture (1:1, w/w) of FO and LO (LFO). Diets were offered as total mixed rations and fed at a rate of 85 g DM/kg BW0.75/d. Oils had no effect (P = 0.52) on DM intake. Linseed oil had no effect (P > 0.05) on ruminal pH or VFA concentrations, while FO shifted rumen fermentation towards propionate at the expense of acetate. Compared with the control, LO increased (P < 0.05) 18:0, cis 18:1 (Δ9, 12-15), trans 18:1 (Δ4-9, 11-16), trans 18:2, geometric isomers of ∆9,11, ∆11,13, and ∆13,15 CLA, trans-8,cis-10 CLA, trans-10,trans-12 CLA, trans-12,trans-14 CLA, and 18:3n-3 flow at the duodenum. Inclusion of FO in the diet resulted in higher (P < 0.05) flows of cis-9 16:1, trans 16:1 (Δ6-13), cis 18:1 (Δ9, 11, and 13), trans 18:1 (Δ6-15), trans 18:2, 20:5n-3, 22:5n-3, and 22:6n-3, and lowered (P < 0.001) 18:0 at the duodenum relative to the control. For most fatty acids at the duodenum responses to LFO were intermediate of FO and LO. However, LFO resulted in higher (P = 0.04) flows of total trans 18:1 than LO and increased (P < 0.01) trans-6 16:1 and trans-12 18:1 at the duodenum compared with FO or LO. Biohydrogenation of cis-9 18:1 and 18:2n-6 in the rumen was independent of treatment, but both FO and LO increased (P < 0.001) the extent of 18:3n-3 biohydrogenation compared with the control. Ruminal 18:3n-3 biohydrogenation was higher (P < 0.001) for LO and LFO than FO, while biohydrogenation of 20:5n-3 and 22:6n-3 in the rumen was marginally lower (P = 0.05) for LFO than FO. In conclusion, LO and FO at 30 g/kg DM altered the biohydrogenation of unsaturated fatty acids in the rumen causing an increase in the flow of specific intermediates at the duodenum, but the potential of these oils fed alone or as a mixture to increase n-3 PUFA at the duodenum in cattle appears limited.

The structure and evolution of the stratospheric vortex in response to natural forcings

The structure and evolution of the Arctic stratospheric polar vortex is assessed during opposing phases of, primarily, the El Niño–Southern Oscillation (ENSO) and the Quasi-Biennial Oscillation (QBO), but the 11 year solar cycle and winters following large volcanic eruptions are also examined. The analysis is performed by taking 2-D moments of vortex potential vorticity (PV) fields which allow the area and centroid of the vortex to be calculated throughout the ERA-40 reanalysis data set (1958–2002). Composites of these diagnostics for the different phases of the natural forcings are then considered. Statistically significant results are found regarding the structure and evolution of the vortex during, in particular, the ENSO and QBO phases. When compared with the more traditional zonal mean zonal wind diagnostic at 60°N, the moment-based diagnostics are far more robust and contain more information regarding the state of the vortex. The study details, for the first time, a comprehensive sequence of events which map the evolution of the vortex during each of the forcings throughout an extended winter period.

Doll Play narratives about starting school in children of socially anxious mothers, and their relation to subsequent child school-based anxiety

Background: Child social anxiety is common, and predicts later emotional and academic impairment. Offspring of socially anxious mothers are at increased risk. It is important to establish whether individual vulnerability to disorder can be identified in young children. Method: The responses of 4.5 year-old children of mothers with social phobia (N = 62) and non-anxious mothers (N = 60) were compared, two months before school entry, using a Doll Play (DP) procedure focused on the social challenge of starting school. DP responses were examined in relation to teacher reports of anxious-depressed symptoms and social worries at the end of the child’s first school term. The role of earlier child behavioral inhibition and attachment, assessed at 14 months, was also considered. Results: Compared to children of non-anxious mothers, children of mothers with social phobia were significantly more likely to give anxiously negative responses in their school DP (OR = 2.57). In turn, negative DP predicted teacher reported anxious-depressed and social worry problems. There were no effects of infant behavioral inhibition or attachment. Conclusion: Vulnerability in young children at risk of anxiety can be identified using Doll Play narratives.

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Toward understanding of differences in current cloud retrievals of ARM ground-based measurements

Accurate observations of cloud microphysical properties are needed for evaluating and improving the representation of cloud processes in climate models and better estimate of the Earth radiative budget. However, large differences are found in current cloud products retrieved from ground-based remote sensing measurements using various retrieval algorithms. Understanding the differences is an important step to address uncertainties in the cloud retrievals. In this study, an in-depth analysis of nine existing ground-based cloud retrievals using ARM remote sensing measurements is carried out. We place emphasis on boundary layer overcast clouds and high level ice clouds, which are the focus of many current retrieval development efforts due to their radiative importance and relatively simple structure. Large systematic discrepancies in cloud microphysical properties are found in these two types of clouds among the nine cloud retrieval products, particularly for the cloud liquid and ice particle effective radius. Note that the differences among some retrieval products are even larger than the prescribed uncertainties reported by the retrieval algorithm developers. It is shown that most of these large differences have their roots in the retrieval theoretical bases, assumptions, as well as input and constraint parameters. This study suggests the need to further validate current retrieval theories and assumptions and even the development of new retrieval algorithms with more observations under different cloud regimes.

Effects of oxidised low density lipoprotein on dendritic cells: a possible immunoregulatory component of the atherogenic micro-environment?

Objective: The objective of this study was to explore the relationship between low density lipoprotein (LDL) and dendritic cell (DC) activation, based upon the hypothesis that reactive oxygen species (ROS)-mediated modification of proteins that may be present in local DC microenvironments could be important as mediators of this activation. Although LDL are known to be oxidised in vivo, and taken up by macrophages during atherogenesis; their effect on DC has not been explored previously. Methods: Human DCs were prepared from peripheral blood monocytes using GM-CSF and IL-4. Plasma LDLs were isolated by sequential gradient centrifugation, oxidised in CuSO4, and oxidation arrested to yield mild, moderate and highly oxidised LDL forms. DCs exposed to these LDLs were investigated using combined phenotypic, functional (autologous T cell activation), morphological and viability assays. Results: Highly-oxidised LDL increased DC HLA-DR, CD40 and CD86 expression, corroborated by increased DC-induced T cell proliferation. Both native and oxidised LDL induced prominent DC clustering. However, high concentrations of highly-oxidised LDL inhibited DC function, due to increased DC apoptosis. Conclusions: This study supports the hypothesis that oxidised LDL are capable of triggering the transition from sentinel to messenger DC. Furthermore, the DC clustering–activation–apoptosis sequence in the presence of different LDL forms is consistent with a regulatory DC role in immunopathogenesis of atheroma. A sequence of initial accumulation of DC, increasing LDL oxidation, and DC-induced T cell activation, may explain why local breach of tolerance can occur. Above a threshold level, however, supervening DC apoptosis limits this, contributing instead to the central plaque core.

Sequence analysis and distribution of an IS3-like insertion element isolated from Salmonella enteritidis

The nucleotide sequence of a 3 kb region immediately upstream of the sef operon operon of Salmonella enteritidis was determined. A 1230 base pair insertion sequence which shared sequence identity (> 75%) with members of the IS3 family was revealed. This element, designated IS1230, had almost identical (90% identity) terminal inverted repeats to Escherichia coli IS3 but unlike other IS3-like sequences lacked the two characteristic open reading frames which encode the putative transposase. S. enteritidis possessed only one copy of this insertion sequence although Southern hybridisation analysis of restriction digests of genomic DNA revealed another fragment located in a region different from the sef operon which hybridised weakly which suggested the presence of an IS1230 homologue. The distribution of IS1230 and IS1230-like elements was shown to be widespread amongst salmonellas and the patterns of restriction fragments which hybridised differed significantly between Salmonella serotypes and it is suggested that IS1230 has potential for development as a differential diagnostic tool.