Ensemble yield simulations: crop and climate uncertainties, sensitivity to temperature and genotypic adaptation to climate change

Estimates of the response of crops to climate change rarely quantify the uncertainty inherent in the simulation of both climate and crops. We present a crop simulation ensemble for a location in India, perturbing the response of both crop and climate under both baseline (12 720 simulations) and doubled-CO2 (171720 simulations) climates. Some simulations used parameter values representing genotypic adaptation to mean temperature change. Firstly, observed and simulated yields in the baseline climate were compared. Secondly, the response of yield to changes in mean temperature was examined and compared to that found in the literature. No consistent response to temperature change was found across studies. Thirdly, the relative contribution of uncertainty in crop and climate simulation to the total uncertainty in projected yield changes was examined. In simulations without genotypic adaptation, most of the uncertainty came from the climate model parameters. Comparison with the simulations with genotypic adaptation and with a previous study suggested that the relatively low crop parameter uncertainty derives from the observational constraints on the crop parameters used in this study. Fourthly, the simulations were used, together with an observed dataset and a simple analysis of crop cardinal temperatures and thermal time, to estimate the potential for adaptation using existing cultivars. The results suggest that the germplasm for complete adaptation of groundnut cultivation in western India to a doubled-CO2 environment may not exist. In conjunction with analyses of germplasm and local management

Use of a crop model ensemble to quantify CO2 stimulation of water-stressed and well-watered crops

Increased atmospheric concentrations of carbon dioxide (CO2) will benefit the yield of most crops. Two free air CO2 enrichment (FACE) meta-analyses have shown increases in yield of between 0 and 73% for C3 crops. Despite this large range, few crop modelling studies quantify the uncertainty inherent in the parameterisation of crop growth and development. We present a novel perturbed-parameter method of crop model simulation, which uses some constraints from observations, that does this. The model used is the groundnut (i.e. peanut; Arachis hypogaea L.) version of the general large-area model for annual crops (GLAM). The conclusions are of relevance to C3 crops in general. The increases in yield simulated by GLAM for doubled CO2 were between 16 and 62%. The difference in mean percentage increase between well-watered and water-stressed simulations was 6.8. These results were compared to FACE and controlled environment studies, and to sensitivity tests on two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., Bell, M.J., 1995. A peanut simulation model. I. Model development and testing. Agron. J. 87, 1085-1093]. The relationship between CO2 and water stress in the experiments and in the models was examined. From a physiological perspective, water-stressed crops are expected to show greater CO2 stimulation than well-watered crops. This expectation has been cited in literature. However, this result is not seen consistently in either the FACE studies or in the crop models. In contrast, leaf-level models of assimilation do consistently show this result. An analysis of the evidence from these models and from the data suggests that scale (canopy versus leaf), model calibration, and model complexity are factors in determining the sign and magnitude of the interaction between CO2 and water stress. We conclude from our study that the statement that 'water-stressed crops show greater CO2 stimulation than well-watered crops' cannot be held to be universally true. We also conclude, preliminarily, that the relationship between water stress and assimilation varies with scale. Accordingly, we provide some suggestions on how studies of a similar nature, using crop models of a range of complexity, could contribute further to understanding the roles of model calibration, model complexity and scale. (C) 2008 Elsevier B.V. All rights reserved.

Identification of a serine protease involved with the regulation of adrenal growth

The adrenal cortex is a dynamic organ in which the cells of the outer cortex continually divide. It is well known that this cellular proliferation is dependent on constant stimulation from peptides derived from the ACTH precursor pro-opiomelanocortin (POMC) because disruption of pituitary corticotroph function results in rapid atrophy of the gland. Previous results from our laboratory have suggested that the adrenal mitogen is a fragment derived from the N-terminal of POMC not containing the gamma-MSH sequence. Because such a peptide is not generated during processing of POMC in the pituitary, we proposed that the mitogen is generated from circulating pro-gamma-MSH by an adrenal protease. Using degenerate oligonucleotides, we identified a secreted serine protease expressed by the adrenal gland that we named adrenal secretory protease (ASP). In the adrenal cortex, expression of ASP is limited to the outer zona glomerulosa/fasciculata, the region where cortical cells are believed to be derived, and is significantly up-regulated during compensatory growth. Y1 adrenocortical cells transfected with a vector expressing an antisense RNA (and thus having reduced levels of endogenous ASP) were found to grow slower than sense controls while also losing their ability to utilize exogenous pro-gamma-MSH in the media supporting a role for ASP in adrenal growth. Digestion of an N-POMC peptide substrate encompassing the residues around the dibasic cleavage site at positions 49/50 with affinity-purified ASP showed cleavage not to occur at the dibasic site but two residues downstream leading us to propose the identity of the adrenal mitogen to be N-POMC (1-52).

A molecular and structural characterization of senescing Arabidopsis siliques and comparison of transcriptional profiles with senescing petals and leaves

Senescence of plant organs is a genetically controlled process that regulates cell death to facilitate nutrient recovery and recycling, and frequently precedes, or is concomitant with, ripening of reproductive structures. In Arabidopsis thaliana, the seeds are contained within a silique, which is itself a photosynthetic organ in the early stages of development and undergoes a programme of senescence prior to dehiscence. A transcriptional analysis of the silique wall was undertaken to identify changes in gene expression during senescence and to correlate these events with ultrastructural changes. The study revealed that the most highly up-regulated genes in senescing silique wall tissues encoded seed storage proteins, and the significance of this finding is discussed. Global transcription profiles of senescing siliques were compared with those from senescing Arabidopsis leaf or petal tissues using microarray datasets and metabolic pathway analysis software (MapMan). In all three tissues, members of NAC and WRKY transcription factor families were up-regulated, but components of the shikimate and cell-wall biosynthetic pathways were down-regulated during senescence. Expression of genes encoding ethylene biosynthesis and action showed more similarity between senescing siliques and petals than between senescing siliques and leaves. Genes involved in autophagy were highly expressed in the late stages of death of all plant tissues studied, but not always during the preceding remobilization phase of senescence. Analyses showed that, during senescence, silique wall tissues exhibited more transcriptional features in common with petals than with leaves. The shared and distinct regulatory events associated with senescence in the three organs are evaluated and discussed.

A randomised clinical trial to assess the effect of total enteral and total parenteral nutritional support on metabolic, inflammatory and oxidative markers in patients with predicted severe acute pancreatitis (APACHE II >= 66)

Background: Total enteral nutrition (TEN) within 48 h of admission has recently been shown to be safe and efficacious as part of the management of severe acute pancreatitis. Our aim was to ascertain the safety of immediate TEN in these patients and the effect of TEN on systemic inflammation, psychological state, oxidative stress, plasma glutamine levels and endotoxaemia. Methods: Patients admitted with predicted severe acute pancreatitis (APACHE II score 15) were randomised to total enteral (TEN; n = 8) or total parenteral nutrition (TPN; n = 9). Measurements of systemic inflammation (C-reactive protein), fatigue ( visual analogue scale), oxidative stress ( plasma thiobarbituric acid- reactive substances), plasma glutamine and anti-endotoxin IgG and IgM antibody concentrations were made on admission and repeated on days 3 and 7 thereafter. Clinical progress was monitored using APACHE II score. Organ failure and complications were recorded. Results: All patients tolerated the feeding regime well with few nutrition-related complications. Fatigue improved in both groups but more rapidly in the TEN group. Oxidative stress was high on admission and rose by similar amounts in both groups. Plasma glutamine concentrations did not change significantly in either group. In the TPN group, 3 patients developed respiratory failure and 3 developed non-respiratory single organ failure. There were no such complications in the TEN group. Hospital stay was shorter in the TEN group [ 7 (4-14) vs. 10 (7-26) days; p = 0.05] as was time to passing flatus and time to opening bowels [1 (0-2) vs. 2 (1-5) days; p = 0.01]. The cost of TEN was considerably less than of TPN. Conclusion: Immediate institution of nutritional support in the form of TEN is safe in predicted severe acute pancreatitis. It is as safe and as efficacious as TPN and may be beneficial in the clinical course of this disease. Copyright (C) 2003 S. Karger AG, Basel and IAP.

HAWAIIAN SKIRT: an F-box gene that regulates organ fusion and growth in arabidopsis

A fast neutron-mutagenized population of Arabidopsis ( Arabidopsis thaliana) Columbia-0 wild-type plants was screened for floral phenotypes and a novel mutant, termed hawaiian skirt ( hws), was identified that failed to shed its reproductive organs. The mutation is the consequence of a 28 bp deletion that introduces a premature amber termination codon into the open reading frame of a putative F-box protein ( At3g61590). The most striking anatomical characteristic of hws plants is seen in flowers where individual sepals are fused along the lower part of their margins. Crossing of the abscission marker, Pro(PGAZAT):beta-glucuronidase, into the mutant reveals that while floral organs are retained it is not the consequence of a failure of abscission zone cells to differentiate. Anatomical analysis indicates that the fusion of sepal margins precludes shedding even though abscission, albeit delayed, does occur. Spatial and temporal characterization, using Pro(HWS):beta-glucuronidase or Pro(HWS):green fluorescent protein fusions, has identified HWS expression to be restricted to the stele and lateral root cap, cotyledonary margins, tip of the stigma, pollen, abscission zones, and developing seeds. Comparative phenotypic analyses performed on the hws mutant, Columbia-0 wild type, and Pro(35S):HWS ectopically expressing lines has revealed that loss of HWS results in greater growth of both aerial and below-ground organs while overexpressing the gene brings about a converse effect. These observations are consistent with HWS playing an important role in regulating plant growth and development.

The effect of postnatal depression on mother-infant interaction, infant response to the still-face perturbation, and performance on an instrumental learning task

A representative community sample of primiparous depressed women and a nondepressed control group were assessed while in interaction with their infants at 2 months postpartum. At 3 months, infants were assessed on the Still-face perturbation of face to face interaction, and a subsample completed an Instrumental Learning paradigm. Compared to nondepressed women, depressed mothers' interactions were both less contingent and less affectively attuned to infant behavior. Postnatal depression did not adversely affect the infant's performance in either the Still-face perturbation or the Instrumental Learning assessment. Maternal responsiveness in interactions at 2 months predicted the infant's performance in the Instrumental Learning assessment but not in the Still-face perturbation. The implications of these findings for theories of infant cognitive and emotional development are discussed.

Response of porcine intestinal in vitro organ culture tissues following exposure to Lactobacillus plantarum JC1 and Salmonella Typhimurium SL1344.

The development of novel intervention strategies for the control of zoonoses caused by bacteria such as Salmonella spp. in livestock requires appropriate experimental models to assess their suitability. Here, a novel porcine intestinal in vitro organ culture (IVOC) model utilizing cell crown (CC) technology (CCIVOC) (Scaffdex) was developed. The CCIVOC model was employed to investigate the characteristics of association of S. enterica serovar Typhimurium strain SL1344 with porcine intestinal tissue following exposure to a Lactobacillus plantarum strain. The association of bacteria to host cells was examined by light microscopy and electron microscopy (EM) after appropriate treatments and staining, while changes in the proteome of porcine jejunal tissues were investigated using quantitative label-free proteomics. Exposure of porcine intestinal mucosal tissues to L. plantarum JC1 did not reduce the numbers of S. Typhimurium bacteria associating to the tissues but was associated with significant (P < 0.005) reductions in the percentages of areas of intestinal IVOC tissues giving positive staining results for acidic mucins. Conversely, the quantity of neutrally charged mucins present within the goblet cells of the IVOC tissues increased significantly (P < 0.05). In addition, tubulin- was expressed at high levels following inoculation of jejunal IVOC tissues with L. plantarum. Although L. plantarum JC1 did not reduce the association of S. Typhimurium strain SL1344 to the jejunal IVOC tissues, detection of increased acidic mucin secretion, host cytoskeletal rearrangements, and proteins involved in the porcine immune response demonstrated that this strain of L. plantarum may contribute to protecting the pig from infections by S. Typhimurium or other pathogens.

Spatiotemporal behavior of the TIGGE medium-range ensemble forecasts

Using the recently-developed mean–variance of logarithms (MVL) diagram, together with the TIGGE archive of medium-range ensemble forecasts from nine different centres, an analysis is presented of the spatiotemporal dynamics of their perturbations, showing how the differences between models and perturbation techniques can explain the shape of their characteristic MVL curves. In particular, a divide is seen between ensembles based on singular vectors or empirical orthogonal functions, and those based on bred vector, Ensemble Transform with Rescaling or Ensemble Kalman Filter techniques. Consideration is also given to the use of the MVL diagram to compare the growth of perturbations within the ensemble with the growth of the forecast error, showing that there is a much closer correspondence for some models than others. Finally, the use of the MVL technique to assist in selecting models for inclusion in a multi-model ensemble is discussed, and an experiment suggested to test its potential in this context.

Assessing the vulnerability of blanket peat to climate change using an ensemble of statistical bioclimatic envelope models

We assessed the vulnerability of blanket peat to climate change in Great Britain using an ensemble of 8 bioclimatic envelope models. We used 4 published models that ranged from simple threshold models, based on total annual precipitation, to Generalised Linear Models (GLMs, based on mean annual temperature). In addition, 4 new models were developed which included measures of water deficit as threshold, classification tree, GLM and generalised additive models (GAM). Models that included measures of both hydrological conditions and maximum temperature provided a better fit to the mapped peat area than models based on hydrological variables alone. Under UKCIP02 projections for high (A1F1) and low (B1) greenhouse gas emission scenarios, 7 out of the 8 models showed a decline in the bioclimatic space associated with blanket peat. Eastern regions (Northumbria, North York Moors, Orkney) were shown to be more vulnerable than higher-altitude, western areas (Highlands, Western Isles and Argyle, Bute and The Trossachs). These results suggest a long-term decline in the distribution of actively growing blanket peat, especially under the high emissions scenario, although it is emphasised that existing peatlands may well persist for decades under a changing climate. Observational data from long-term monitoring and manipulation experiments in combination with process-based models are required to explore the nature and magnitude of climate change impacts on these vulnerable areas more fully.

Initiation of deep convection at marginal instability in an ensemble of mesoscale models: a case-study from COPS

The present study investigates the initiation of precipitating deep convection in an ensemble of convection-resolving mesoscale models. Results of eight different model runs from five non-hydrostatic models are compared for a case of the Convective and Orographically-induced Precipitation Study (COPS). An isolated convective cell initiated east of the Black Forest crest in southwest Germany, although convective available potential energy was only moderate and convective inhibition was high. Measurements revealed that, due to the absence of synoptic forcing, convection was initiated by local processes related to the orography. In particular, the lifting by low-level convergence in the planetary boundary layer is assumed to be the dominant process on that day. The models used different configurations as well as different initial and boundary conditions. By comparing the different model performance with each other and with measurements, the processes which need to be well represented to initiate convection at the right place and time are discussed. Besides an accurate specification of the thermodynamic and kinematic fields, the results highlight the role of boundary-layer convergence features for quantitative precipitation forecasts in mountainous terrain.

Ensemble prediction of transitions of the North Atlantic eddy-driven jet

The consistency of ensemble forecasts from three global medium-range prediction systems with the observed transition behaviour of a three-cluster model of the North Atlantic eddy-driven jet is examined. The three clusters consist of a mid jet cluster taken to represent an undisturbed jet and south and north jet clusters representing southward and northward shifts of the jet. The ensemble forecasts span a period of three extended winters (October–February) from October 2007–February 2010. The mean probabilities of transitions between the clusters calculated from the ensemble forecasts are compared with those calculated from a 23-extended-winter climatology taken from the European Centre for Medium-Range Weather Forecasts 40-Year Re-analysis (ERA40) dataset. No evidence of a drift with increasing lead time of the ensemble forecast transition probabilities towards values inconsistent with the 23-extended-winter climatology is found. The ensemble forecasts of transition probabilities are found to have positive Brier Skill at 15 day lead times. It is found that for the three-extended-winter forecast set, probabilistic forecasts initialized in the north jet cluster are generally less skilful than those initialized in the other clusters. This is consistent with the shorter persistence time-scale of the north jet cluster observed in the ERA40 23-extended-winter climatology. Copyright © 2011 Royal Meteorological Society

Comparing the UK Met Office climate prediction system DePreSys with idealized predictability in the HadCM3 model

The initial condition effect on climate prediction skill over a 2-year hindcast time-scale has been assessed from ensemble HadCM3 climate model runs using anomaly initialization over the period 1990–2001, and making comparisons with runs without initialization (equivalent to climatological conditions), and to anomaly persistence. It is shown that the assimilation improves the prediction skill in the first year globally, and in a number of limited areas out into the second year. Skill in hindcasting surface air temperature anomalies is most marked over ocean areas, and is coincident with areas of high sea surface temperature and ocean heat content skill. Skill improvement over land areas is much more limited but is still detectable in some cases. We found little difference in the skill of hindcasts using three different sets of ocean initial conditions, and we obtained the best results by combining these to form a grand ensemble hindcast set. Results are also compared with the idealized predictability studies of Collins (Clim. Dynam. 2002; 19: 671–692), which used the same model. The maximum lead time for which initialization gives enhanced skill over runs without initialization varies in different regions but is very similar to lead times found in the idealized studies, therefore strongly supporting the process representation in the model as well as its use for operational predictions. The limited 12-year period of the study, however, means that the regional details of model skill should probably be further assessed under a wider range of observational conditions.

State estimation using the particle filter with mode tracking

A particle filter is a data assimilation scheme that employs a fully nonlinear, non-Gaussian analysis step. Unfortunately as the size of the state grows the number of ensemble members required for the particle filter to converge to the true solution increases exponentially. To overcome this Vaswani [Vaswani N. 2008. IEEE Trans Signal Process 56:4583–97] proposed a new method known as mode tracking to improve the efficiency of the particle filter. When mode tracking, the state is split into two subspaces. One subspace is forecast using the particle filter, the other is treated so that its values are set equal to the mode of the marginal pdf. There are many ways to split the state. One hypothesis is that the best results should be obtained from the particle filter with mode tracking when we mode track the maximum number of unimodal dimensions. The aim of this paper is to test this hypothesis using the three dimensional stochastic Lorenz equations with direct observations. It is found that mode tracking the maximum number of unimodal dimensions does not always provide the best result. The best choice of states to mode track depends on the number of particles used and the accuracy and frequency of the observations.

Sensory and instrumental analyses of volatiles generated during the extrusion cooking of oat flours

Three batches of oats were extruded under four combinations of process temperature (150 or 180 °C) and process moisture (14.5 and 18%). Two of the extrudates were evaluated by a sensory panel, and three were analyzed by GC-MS. Maillard reaction products, such as pyrazines, pyrroles, furans, and sulfur-containing compounds, were found in the most severely processed extrudates (high-temperature, low-moisture). These extrudates were also described by the assessors as having toasted cereal attributes. Lipid degradation products, such as alkanals, 2-alkenals, and 2,4-alkadienals, were found at much higher levels in the extrudates of the oat flour that had been debranned. It contained lower protein and fiber levels than the others and showed increased lipase activity. Extrudates from these samples also had significantly lower levels of Maillard reaction products that correlated, in the sensory analysis, with terms such as stale oil and oatmeal. Linoleic acid was added to a fourth oat flour to simulate the result of increased lipase activity, and GC-MS analysis showed both an increase in lipid degradation products and a decrease in Maillard reaction products.