Assessing the uncertainty associated with intermittent rainfall fields

[1] In many practical situations where spatial rainfall estimates are needed, rainfall occurs as a spatially intermittent phenomenon. An efficient geostatistical method for rainfall estimation in the case of intermittency has previously been published and comprises the estimation of two independent components: a binary random function for modeling the intermittency and a continuous random function that models the rainfall inside the rainy areas. The final rainfall estimates are obtained as the product of the estimates of these two random functions. However the published approach does not contain a method for estimation of uncertainties. The contribution of this paper is the presentation of the indicator maximum likelihood estimator from which the local conditional distribution of the rainfall value at any location may be derived using an ensemble approach. From the conditional distribution, representations of uncertainty such as the estimation variance and confidence intervals can be obtained. An approximation to the variance can be calculated more simply by assuming rainfall intensity is independent of location within the rainy area. The methodology has been validated using simulated and real rainfall data sets. The results of these case studies show good agreement between predicted uncertainties and measured errors obtained from the validation data.

Global hydrology modelling and uncertainty: running multiple ensembles with a campus grid.

Uncertainties associated with the representation of various physical processes in global climate models (GCMs) mean that, when projections from GCMs are used in climate change impact studies, the uncertainty propagates through to the impact estimates. A complete treatment of this ‘climate model structural uncertainty’ is necessary so that decision-makers are presented with an uncertainty range around the impact estimates. This uncertainty is often underexplored owing to the human and computer processing time required to perform the numerous simulations. Here, we present a 189-member ensemble of global river runoff and water resource stress simulations that adequately address this uncertainty. Following several adaptations and modifications, the ensemble creation time has been reduced from 750 h on a typical single-processor personal computer to 9 h of high-throughput computing on the University of Reading Campus Grid. Here, we outline the changes that had to be made to the hydrological impacts model and to the Campus Grid, and present the main results. We show that, although there is considerable uncertainty in both the magnitude and the sign of regional runoff changes across different GCMs with climate change, there is much less uncertainty in runoff changes for regions that experience large runoff increases (e.g. the high northern latitudes and Central Asia) and large runoff decreases (e.g. the Mediterranean). Furthermore, there is consensus that the percentage of the global population at risk to water resource stress will increase with climate change.

TIGGE: Comparison of the prediction of northern hemisphere extratropical cyclones by different ensemble prediction systems

The Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE) is a World Weather Research Programme project. One of its main objectives is to enhance collaboration on the development of ensemble prediction between operational centers and universities by increasing the availability of ensemble prediction system (EPS) data for research. This study analyzes the prediction of Northern Hemisphere extratropical cyclones by nine different EPSs archived as part of the TIGGE project for the 6-month time period of 1 February 2008–31 July 2008, which included a sample of 774 cyclones. An objective feature tracking method has been used to identify and track the cyclones along the forecast trajectories. Forecast verification statistics have then been produced [using the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis as the truth] for cyclone position, intensity, and propagation speed, showing large differences between the different EPSs. The results show that the ECMWF ensemble mean and control have the highest level of skill for all cyclone properties. The Japanese Meteorological Administration (JMA), the National Centers for Environmental Prediction (NCEP), the Met Office (UKMO), and the Canadian Meteorological Centre (CMC) have 1 day less skill for the position of cyclones throughout the forecast range. The relative performance of the different EPSs remains the same for cyclone intensity except for NCEP, which has larger errors than for position. NCEP, the Centro de Previsão de Tempo e Estudos Climáticos (CPTEC), and the Australian Bureau of Meteorology (BoM) all have faster intensity error growth in the earlier part of the forecast. They are also very underdispersive and significantly underpredict intensities, perhaps due to the comparatively low spatial resolutions of these EPSs not being able to accurately model the tilted structure essential to cyclone growth and decay. There is very little difference between the levels of skill of the ensemble mean and control for cyclone position, but the ensemble mean provides an advantage over the control for all EPSs except CPTEC in cyclone intensity and there is an advantage for propagation speed for all EPSs. ECMWF and JMA have an excellent spread–skill relationship for cyclone position. The EPSs are all much more underdispersive for cyclone intensity and propagation speed than for position, with ECMWF and CMC performing best for intensity and CMC performing best for propagation speed. ECMWF is the only EPS to consistently overpredict cyclone intensity, although the bias is small. BoM, NCEP, UKMO, and CPTEC significantly underpredict intensity and, interestingly, all the EPSs underpredict the propagation speed, that is, the cyclones move too slowly on average in all EPSs.

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