Calibrating flood inundation models: identifying and addressing uncertainty in satellite observed flood extents

Improvements in the resolution of satellite imagery have enabled extraction of water surface elevations at the margins of the flood. Comparison between modelled and observed water surface elevations provides a new means for calibrating and validating flood inundation models, however the uncertainty in this observed data has yet to be addressed. Here a flood inundation model is calibrated using a probabilistic treatment of the observed data. A LiDAR guided snake algorithm is used to determine an outline of a flood event in 2006 on the River Dee, North Wales, UK, using a 12.5m ERS-1 image. Points at approximately 100m intervals along this outline are selected, and the water surface elevation recorded as the LiDAR DEM elevation at each point. With a planar water surface from the gauged upstream to downstream water elevations as an approximation, the water surface elevations at points along this flooded extent are compared to their ‘expected’ value. The pattern of errors between the two show a roughly normal distribution, however when plotted against coordinates there is obvious spatial autocorrelation. The source of this spatial dependency is investigated by comparing errors to the slope gradient and aspect of the LiDAR DEM. A LISFLOOD-FP model of the flood event is set-up to investigate the effect of observed data uncertainty on the calibration of flood inundation models. Multiple simulations are run using different combinations of friction parameters, from which the optimum parameter set will be selected. For each simulation a T-test is used to quantify the fit between modelled and observed water surface elevations. The points chosen for use in this T-test are selected based on their error. The criteria for selection enables evaluation of the sensitivity of the choice of optimum parameter set to uncertainty in the observed data. This work explores the observed data in detail and highlights possible causes of error. The identification of significant error (RMSE = 0.8m) between approximate expected and actual observed elevations from the remotely sensed data emphasises the limitations of using this data in a deterministic manner within the calibration process. These limitations are addressed by developing a new probabilistic approach to using the observed data.

Enhancement of the aerosol direct radiative effect by semi-volatile aerosol components: airborne measurements in North-Western Europe

A case study of atmospheric aerosol measurements exploring the impact of the vertical distribution of aerosol chemical composition upon the radiative budget in North-Western Europe is presented. Sub-micron aerosol chemical composition was measured by an Aerodyne Aerosol Mass Spectrometer (AMS) on both an airborne platform and a ground-based site at Cabauw in the Netherlands. The examined period in May 2008 was characterised by enhanced pollution loadings in North-Western Europe and was dominated by ammonium nitrate and Organic Matter (OM). Both ammonium nitrate and OM were observed to increase with altitude in the atmospheric boundary layer. This is primarily attributed to partitioning of semi-volatile gas phase species to the particle phase at reduced temperature and enhanced relative humidity. Increased ammonium nitrate concentrations in particular were found to strongly increase the ambient scattering potential of the aerosol burden, which was a consequence of the large amount of associated water as well as the enhanced mass. During particularly polluted conditions, increases in aerosol optical depth of 50–100% were estimated to occur due to the observed increase in secondary aerosol mass and associated water uptake. Furthermore, the single scattering albedo was also shown to increase with height in the boundary layer. These enhancements combined to increase the negative direct aerosol radiative forcing by close to a factor of two at the median percentile level. Such increases have major ramifications for regional climate predictions as semi-volatile components are often not included in aerosol models. The results presented here provide an ideal opportunity to test regional and global representations of both the aerosol vertical distribution and subsequent impacts in North-Western Europe. North-Western Europe can be viewed as an analogue for the possible future air quality over other polluted regions of the Northern Hemisphere, where substantial reductions in sulphur dioxide emissions have yet to occur. Anticipated reductions in sulphur dioxide in polluted regions will result in an increase in the availability of ammonia to form ammonium nitrate as opposed to ammonium sulphate. This will be most important where intensive agricultural practises occur. Our observations over North-Western Europe, a region where sulphur dioxide emissions have already been reduced, indicate that failure to include the semi-volatile behaviour of ammonium nitrate will result in significant errors in predicted aerosol direct radiative forcing. Such errors will be particularly significant on regional scales.

The effect of incremental replacement of wheat with soya hulls in diets for Jersey cows on lactational performance, diet digestibility and feeding behaviour

Eight Jersey cows were used in two balanced 4 x 4 Latin Squares to investigate the effects of replacement of dietary starch with non-forage fibre on productivity, diet digestibility and feeding behaviour. Total-mixed rations consisted of maize silage, grass silage and a soyabean meal-based concentrate mixture, each at 250g/kg DM, with the remaining 250g consisting of cracked wheat/soya hulls (SH) in the ratios of 250:0, 167:83; 83:167 and 0:250 g, respectively, for treatments SH0, SH83, SH167 and SH250. Starch concentrations were 302, 248, 193 and 140g/kg DM, and NDF concentrations were 316, 355, 394 and 434g/kg DM, for treatments SHO, SH83, SH167 and SH250, respectively. Total eating time increased (p < 0.05) as SH inclusion increased, but total rumination time was unaffected. Digestibility of DM, organic matter and starch declined (p < 0.01) as SH inclusion increased, whilst digestibility of NDF and ADF increased (p < 0.01). Dry-matter intake tended to decline with increasing SH, whilst bodyweight, milk yield and fat and lactose concentrations were unaffected by treatment. Milk protein concentration decreased (p < 0.01) as SH level increased. Feed conversion efficiency improved (p < 0.05) as SH inclusion rose, but it was not possible to determine whether this was due to the increased fibre levels alone, or the favourable effect on rumen fermentation of decreasing starch levels. (c) 2006 Elsevier B.V. All rights reserved.

Model selection for the mixed logit with Bayesian estimation

In this paper, the mixed logit (ML) using Bayesian methods was employed to examine willingness-to-pay (WTP) to consume bread produced with reduced levels of pesticides so as to ameliorate environmental quality, from data generated by a choice experiment. Model comparison used the marginal likelihood, which is preferable for Bayesian model comparison and testing. Models containing constant and random parameters for a number of distributions were considered, along with models in ‘preference space’ and ‘WTP space’ as well as those allowing for misreporting. We found: strong support for the ML estimated in WTP space; little support for fixing the price coefficient a common practice advocated and adopted in the environmental economics literature; and, weak evidence for misreporting.

Effect of replacing grass silage with maize silage in the diet on bovine milk fatty acid composition

Even though extensive research has examined the role of nutrition on milk fat composition, there is less information on the impact of forages on milk fatty acid (FA) composition. In the current study, the effect of replacing grass silage (GS) with maize silage (MS) as part of a total mixed ration on animal performance and milk FA composition was examined using eight multiparous mid-lactation cows in a replicated 4 X 4 Latin square with 28-day experimental periods. Four treatments comprised the stepwise replacement of GS with MS (0, 160, 334 and 500 g/kg dry matter (DM)) in diets containing a 54:46 forage: concentrate ratio on a DM basis. Replacing GS with MS increased (P < 0.001) the DM intake, milk yield and milk protein content. Incremental replacement of GS with MS in the diet enhanced linearly (P < 0.001) the proportions of 6:0-14:0, decreased (P < 0.01) the 16:0 concentrations, but had no effect on the total milk fat saturated fatty acid content. Inclusion of MS altered the distribution of trans-18:1 isomers and enhanced (P < 0.05) total trans monounsaturated fatty acid and total conjugated linoleic acid content. Milk total n-3 polyunsaturated fatty acid (PUFA) content decreased with higher amounts of MS in the diet and n-6 PUFA concentration increased, leading to an elevated n-6: n-3 PUFA ratio. Despite some beneficial changes associated with the replacement of GS with MS, the overall effects on milk FA composition would not be expected to substantially improve long-term human health. However the role of forages on milk fat composition must also be balanced against the increases in total milk and protein yield on diets containing higher proportions of MS.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.

Effect of method of application of a fibrolytic enzyme product on digestive processes and milk production in Holstein-Friesian cows

Four multiparous cows with cannulas in the rumen and proximal duodenum were used in early lactation in a 4 x 4 Latin square experiment to investigate the effect of method of application of a fibrolytic enzyme product on digestive processes and milk production. The cows were given ad libitum a total mixed ration (TMR) composed of 57% (dry matter basis) forage (3:1 corn silage:grass silage) and 43% concentrates. The TMR contained (g/kg dry matter): 274 neutral detergent fiber, 295 starch, 180 crude protein. Treatments were TMR alone or TMR with the enzyme product added (2 kg/1000 kg TMR dry matter) either sprayed on the TMR 1 h before the morning feed (TMR-E), sprayed only on the concentrate the day before feeding (Concs-E), or infused into the rumen for 14 h/d (Rumen-E). There Was no significant effect on either feed intake or milk yield but both were highest on TMR-E. Rumen digestibility of dry matter, organic matter, and starch was unaffected by the enzyme. Digestibility of NDF was lowest on TMR-E in the rumen but highest postruminally. Total Tract digestibility was highest on TMR-E for dry matter, organic matter, and starch but treatment differences were nonsignificant for neutral detergent fiber: Corn silage stover retention time in the rumen was reduced by all enzyme treatments but postruminal transit time vas increased so the decline in total tract retention. time with enzymes was not significant. It is suggested that the tendency for enzymes to reduce particle retention time in the rumen may, by reducing the time available for fibrolysis to occur, at least partly explain the variability in the reported responses to enzyme treatment.

The effect of wheat straw residue on the emergence and early growth of sugar beet (Beta vulgaris) and oilseed rape (Brassica napus)

The management of straw residue can be a concern in non-inversion tillage systems where straw tends to be incorporated at shallow depths or left on the soil surface. This can lead to poor crop establishment because straw residue can impede or hinder crop emergence and growth. Small container-based experiments were undertaken using varying amounts of wheat straw residue either incorporated or placed oil the soil surface. The effects on (lays to seedling emergence, percentage emergence, seedling dry-weight and soil temperature using sugar beet and oilseed rape were investigated because these crops often follow wheat in a cropping sequence. The position of the straw residue was found to be the primary factor in reducing crop emergence and growth. Increasing the amount of straw residue (from 3.3 t ha(-1) to 6.7 t ha(-1)) did not show any consistent trends in reducing crop emergence or growth. However, in some instances, results indicated that an interaction between the position and the amount of straw residue Occurred particularly when the straw and seed was placed on the soil surface. Straw placed on the soil surface significantly reduced mean day-time soil temperature by approximately 2.5 degrees C compared to no residue. When the seed and straw was placed on the soil Surface a lack of seed-to-soil contact caused a reduction in emergence by approximately 30% because of the restriction in available moisture that limited the ability for seed imbibition. This trend was reversed when the seed was placed in the soil, but with straw residue still on the soil surface, because the surface straw was likely to reduce moisture evaporation and improved seed-to-soil contact that led to rapid emergence. In general, when straw was mixed in or placed on the soil surface along with the seed, sugar beet and oilseed rape emergence and early growth biomass was significantly restricted by approximately 50% compared to no residue. The consequences of placing seed with or near to straw residue have been shown to cause a restriction in crop establishment. In both oilseed tape and sugar beet, this could lead to a reduction in final crop densities, poor, uneven growth and potentially lower yields that could lower financial margins. Therefore, if farmers are planning to use non-inversion tillage methods for crop establishment, the management and removal of straw residue from near or above the seed is considered important for successful crop establishment. (C) 2008 Elsevier B.V. All rights reserved.

Phylogenetic mixture models can reduce node-density artifacts

We investigate the performance of phylogenetic mixture models in reducing a well-known and pervasive artifact of phylogenetic inference known as the node-density effect, comparing them to partitioned analyses of the same data. The node-density effect refers to the tendency for the amount of evolutionary change in longer branches of phylogenies to be underestimated compared to that in regions of the tree where there are more nodes and thus branches are typically shorter. Mixture models allow more than one model of sequence evolution to describe the sites in an alignment without prior knowledge of the evolutionary processes that characterize the data or how they correspond to different sites. If multiple evolutionary patterns are common in sequence evolution, mixture models may be capable of reducing node-density effects by characterizing the evolutionary processes more accurately. In gene-sequence alignments simulated to have heterogeneous patterns of evolution, we find that mixture models can reduce node-density effects to negligible levels or remove them altogether, performing as well as partitioned analyses based on the known simulated patterns. The mixture models achieve this without knowledge of the patterns that generated the data and even in some cases without specifying the full or true model of sequence evolution known to underlie the data. The latter result is especially important in real applications, as the true model of evolution is seldom known. We find the same patterns of results for two real data sets with evidence of complex patterns of sequence evolution: mixture models substantially reduced node-density effects and returned better likelihoods compared to partitioning models specifically fitted to these data. We suggest that the presence of more than one pattern of evolution in the data is a common source of error in phylogenetic inference and that mixture models can often detect these patterns even without prior knowledge of their presence in the data. Routine use of mixture models alongside other approaches to phylogenetic inference may often reveal hidden or unexpected patterns of sequence evolution and can improve phylogenetic inference.

Fitting models for the joint action of two drugs using SAS(R)

Combinations of drugs are increasingly being used for a wide variety of diseases and conditions. A pre-clinical study may allow the investigation of the response at a large number of dose combinations. In determining the response to a drug combination, interest may lie in seeking evidence of synergism, in which the joint action is greater than the actions of the individual drugs, or of antagonism, in which it is less. Two well-known response surface models representing no interaction are Loewe additivity and Bliss independence, and Loewe or Bliss synergism or antagonism is defined relative to these. We illustrate an approach to fitting these models for the case in which the marginal single drug dose-response relationships are represented by four-parameter logistic curves with common upper and lower limits, and where the response variable is normally distributed with a common variance about the dose-response curve. When the dose-response curves are not parallel, the relative potency of the two drugs varies according to the magnitude of the desired effect and the models for Loewe additivity and synergism/antagonism cannot be explicitly expressed. We present an iterative approach to fitting these models without the assumption of parallel dose-response curves. A goodness-of-fit test based on residuals is also described. Implementation using the SAS NLIN procedure is illustrated using data from a pre-clinical study. Copyright © 2007 John Wiley & Sons, Ltd.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Antiviral effects of antisense morpholino oligomers in murine coronavirus infection models

The recent emergence of novel pathogenic human and animal coronaviruses has highlighted the need for antiviral therapies that are effective against a spectrum of these viruses. We have used several strains of murine hepatitis virus (MHV) in cell culture and in vivo in mouse models to investigate the antiviral characteristics of peptide-conjugated antisense phosphorodiamidate morpholino oligomers (P-PMOs). Ten P-PMOs directed against various target sites in the viral genome were tested in cell culture, and one of these (5TERM), which was complementary to the 5' terminus of the genomic RNA, was effective against six strains of MHV. Further studies were carried out with various arginine-rich peptides conjugated to the 5TERM PMO sequence in order to evaluate efficacy and toxicity and thereby select candidates for in vivo testing. In uninfected mice, prolonged P-PMO treatment did not result in weight loss or detectable histopathologic changes. 5TERM P-PMO treatment reduced viral titers in target organs and protected mice against virus-induced tissue damage. Prophylactic 5TERM P-PMO treatment decreased the amount of weight loss associated with infection under most experimental conditions. Treatment also prolonged survival in two lethal challenge models. In some cases of high-dose viral inoculation followed by delayed treatment, 5TERM P-PMO treatment was not protective and increased morbidity in the treated group, suggesting that P-PMO may cause toxic effects in diseased mice that were not apparent in the uninfected animals. However, the strong antiviral effect observed suggests that with further development, P-PMO may provide an effective therapeutic approach against a broad range of coronavirus infections.

A comparison of three different models for estimating relative risk in meta-analysis of clinical trials under unobserved heterogeneity

We focus on the comparison of three statistical models used to estimate the treatment effect in metaanalysis when individually pooled data are available. The models are two conventional models, namely a multi-level and a model based upon an approximate likelihood, and a newly developed model, the profile likelihood model which might be viewed as an extension of the Mantel-Haenszel approach. To exemplify these methods, we use results from a meta-analysis of 22 trials to prevent respiratory tract infections. We show that by using the multi-level approach, in the case of baseline heterogeneity, the number of clusters or components is considerably over-estimated. The approximate and profile likelihood method showed nearly the same pattern for the treatment effect distribution. To provide more evidence two simulation studies are accomplished. The profile likelihood can be considered as a clear alternative to the approximate likelihood model. In the case of strong baseline heterogeneity, the profile likelihood method shows superior behaviour when compared with the multi-level model. Copyright (C) 2006 John Wiley & Sons, Ltd.

The effect of potato variety mixtures on epidemics of late blight, in relation to plot size and level of resistance

Potatoes of a number of varieties of contrasting levels of resistance were planted in pure or mixed stands in four experiments over 3 years. Three experiments compared the late blight severity and progress in mixtures with that in pure stands. Disease on susceptible or moderately resistant varieties typical of those in commercial use was similar in mixtures and pure stands. In 2 of 3 years, there were slight reductions on cv. Sante, which is moderately susceptible, in mixture with cv. Cara, which is moderately resistant. Cara was unaffected by this mixture. Mixtures of an immune or near-immune partner with Cara or Sante substantially reduced disease on the latter. The effect of the size of plots of individual varieties or mixtures on blight severity was compared in two experiments. Larger plots had a greater area under the disease progress curve, but the average rate of disease progress was greater in smaller plots; this may be because most disease progress took place later, under more favourable conditions, in the smaller plots. In one experiment, two planting densities were used. Density had no effect on disease and did not interact with mixture effects. The overall conclusion is that, while mixtures of potato varieties may be desirable for other reasons, they do not offer any improvement on the average of the disease resistance of the components.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.