The potential to narrow uncertainty in projections of stratospheric ozone over the 21st century

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an "ensemble of opportunity" of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21st century, up-to and after the time when ozone concentrations return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

Bayesian estimation of co-integrating thresholds in the term structure of interest rates.

Threshold Error Correction Models are used to analyse the term structure of interest Rates. The paper develops and uses a generalisation of existing models that encompasses both the Band and Equilibrium threshold models of [Balke and Fomby ((1997) Threshold cointegration. Int Econ Rev 38(3):627–645)] and estimates this model using a Bayesian approach. Evidence is found for threshold effects in pairs of longer rates but not in pairs of short rates. The Band threshold model is supported in preference to the Equilibrium model.

Bayesian estimation of willingness-to-pay where respondents mis-report their preferences

We introduce a modified conditional logit model that takes account of uncertainty associated with mis-reporting in revealed preference experiments estimating willingness-to-pay (WTP). Like Hausman et al. [Journal of Econometrics (1988) Vol. 87, pp. 239-269], our model captures the extent and direction of uncertainty by respondents. Using a Bayesian methodology, we apply our model to a choice modelling (CM) data set examining UK consumer preferences for non-pesticide food. We compare the results of our model with the Hausman model. WTP estimates are produced for different groups of consumers and we find that modified estimates of WTP, that take account of mis-reporting, are substantially revised downwards. We find a significant proportion of respondents mis-reporting in favour of the non-pesticide option. Finally, with this data set, Bayes factors suggest that our model is preferred to the Hausman model.

Estimation of the stoichiometry of volatile fatty acid production in the rumen of lactating cows

The purpose of this study was to improve the prediction of the quantity and type of Volatile Fatty Acids (VFA) produced from fermented substrate in the rumen of lactating cows. A model was formulated that describes the conversion of substrate (soluble carbohydrates, starch, hemi-cellulose, cellulose, and protein) into VFA (acetate, propionate, butyrate, and other VFA). Inputs to the model were observed rates of true rumen digestion of substrates, whereas outputs were observed molar proportions of VFA in rumen fluid. A literature survey generated data of 182 diets (96 roughage and 86 concentrate diets). Coefficient values that define the conversion of a specific substrate into VFA were estimated meta-analytically by regression of the model against observed VFA molar proportions using non-linear regression techniques. Coefficient estimates significantly differed for acetate and propionate production in particular, between different types of substrate and between roughage and concentrate diets. Deviations of fitted from observed VFA molar proportions could be attributed to random error for 100%. In addition to regression against observed data, simulation studies were performed to investigate the potential of the estimation method. Fitted coefficient estimates from simulated data sets appeared accurate, as well as fitted rates of VFA production, although the model accounted for only a small fraction (maximally 45%) of the variation in VFA molar proportions. The simulation results showed that the latter result was merely a consequence of the statistical analysis chosen and should not be interpreted as an indication of inaccuracy of coefficient estimates. Deviations between fitted and observed values corresponded to those obtained in simulations. (c) 2005 Elsevier Ltd. All rights reserved.

Probabilistic simulations of crop yield over western India using the DEMETER seasonal hindcast ensembles

Process-based integrated modelling of weather and crop yield over large areas is becoming an important research topic. The production of the DEMETER ensemble hindcasts of weather allows this work to be carried out in a probabilistic framework. In this study, ensembles of crop yield (groundnut, Arachis hypogaea L.) were produced for 10 2.5 degrees x 2.5 degrees grid cells in western India using the DEMETER ensembles and the general large-area model (GLAM) for annual crops. Four key issues are addressed by this study. First, crop model calibration methods for use with weather ensemble data are assessed. Calibration using yield ensembles was more successful than calibration using reanalysis data (the European Centre for Medium-Range Weather Forecasts 40-yr reanalysis, ERA40). Secondly, the potential for probabilistic forecasting of crop failure is examined. The hindcasts show skill in the prediction of crop failure, with more severe failures being more predictable. Thirdly, the use of yield ensemble means to predict interannual variability in crop yield is examined and their skill assessed relative to baseline simulations using ERA40. The accuracy of multi-model yield ensemble means is equal to or greater than the accuracy using ERA40. Fourthly, the impact of two key uncertainties, sowing window and spatial scale, is briefly examined. The impact of uncertainty in the sowing window is greater with ERA40 than with the multi-model yield ensemble mean. Subgrid heterogeneity affects model accuracy: where correlations are low on the grid scale, they may be significantly positive on the subgrid scale. The implications of the results of this study for yield forecasting on seasonal time-scales are as follows. (i) There is the potential for probabilistic forecasting of crop failure (defined by a threshold yield value); forecasting of yield terciles shows less potential. (ii) Any improvement in the skill of climate models has the potential to translate into improved deterministic yield prediction. (iii) Whilst model input uncertainties are important, uncertainty in the sowing window may not require specific modelling. The implications of the results of this study for yield forecasting on multidecadal (climate change) time-scales are as follows. (i) The skill in the ensemble mean suggests that the perturbation, within uncertainty bounds, of crop and climate parameters, could potentially average out some of the errors associated with mean yield prediction. (ii) For a given technology trend, decadal fluctuations in the yield-gap parameter used by GLAM may be relatively small, implying some predictability on those time-scales.

Estimation of milk production from smallholder dairy cattle in the coastal lowlands of Kenya

A study was carried out on 92 smallholder farms in Kwale district in Coast Province of Kenya to estimate the milk yield. The effect of concentrate feed supplementation on milk yield was also evaluated. Data were collected during a one-year observational longitudinal study. Analysis was done for 371 observations following 63 calving events. The mean annual milk offtake was estimated at 2021 kg/cow. Forty-nine (77.8%) of the lactating cows were supplemented with concentrate feeds at varying rates of 0.5-3.0 kg/cow per day. Supplementary feeding of lactating cows led to a significantly higher mean daily milk yield compared to non-supplemented cows throughout the year (p<0.05). The mean annual milk offtake from supplemented cows (2195 kg/cow) was 18.6% more than offtake from non-supplemented cows, a difference that was statistically significant (p<0.05). Therefore, supplementary feeding of commercial feed concentrates was a rational management practice. It was also concluded that milk production from smallholder dairy cows in the coastal lowlands of Kenya was comparable to that from similar production systems but lower than national targets.

Cross-sectional estimation of Babesia bovis antibody prevalence in cattle in two contrasting dairying areas in Tanzania

The crude prevalence of antibodies to Babesia bovis infection in cattle was estimated by serology using indirect ELISA during the period January to April, 1999. Sera were obtained from 1395 dairy cattle (of all ages, sexes and breeds) on smallholder farms, the majority being kept under a zero grazing regime. The crude prevalence of antibodies to Babesia bovis was 6 % for Tanga and 12 % for Iringa. The forces of infection based on the age sero-prevalence profile, were estimated at six for Iringa and four for Tanga per 100 cattle years-risk, respectively. Using random effect logistic regression as the analytical method, the factors (variables) of age, source of animals and geographic location were hypothesised to be associated with sero-positivity of Babesia bovis in the two regions.

Investigating the effects of predator removal and habitat management on nest success and breeding population size of a farmland passerine: a case study

Passerines are especially vulnerable to predation at the pre-independence stage. Although the role of nest success in British farmland passerine declines is contentious, improvement in nest success through sympathetic management could play a role in their reversal. Because habitat is known to interact with predation, management options for mitigation will need to consider effects of nest predation. We present results from an observational study of a population of Common Blackbird Turdus merula on a farm which has experienced a range of agri-environment and game-management options, including a period with nest predator control, as a case study to address some of these issues. We used an information theoretic model comparison procedure to look for evidence of interactions between habitat and nest predation, and then asked whether habitat management and nest predator abundances could explain population trends at the site through their effects on nest success. Interactions were detected between measures of predator abundance and habitat variables, and these varied with nest stage - habitat within the vicinity of the nest appeared to be important at the egg stage, and nest-placement characteristics were important at the nestling stage. Although predator control appeared to have a positive influence on Blackbird breeding population size, the non-experimental set-up meant we could not eliminate other potential explanations. Variation in breeding population size did not appear to be influenced by variation in nest success alone. Our study demonstrates that observational data can only go so far in detection of such effects, and we discuss how it might be taken further. Agri-environment and game-management techniques are likely to influence nest predation pressure on farmland passerines, but the patterns, mechanisms and importance to population processes remain not wholly understood.

Investigating the effects of predator removal and habitat management on nest success and breeding population size of a farmland passerine: a case study

Passerines are especially vulnerable to predation at the pre-independence stage. Although the role of nest success in British farmland passerine declines is contentious, improvement in nest success through sympathetic management could play a role in their reversal. Because habitat is known to interact with predation, management options for mitigation will need to consider effects of nest predation. We present results from an observational study of a population of Common Blackbird Turdus merula on a farm which has experienced a range of agri-environment and game-management options, including a period with nest predator control, as a case study to address some of these issues. We used an information theoretic model comparison procedure to look for evidence of interactions between habitat and nest predation, and then asked whether habitat management and nest predator abundances could explain population trends at the site through their effects on nest success. Interactions were detected between measures of predator abundance and habitat variables, and these varied with nest stage - habitat within the vicinity of the nest appeared to be important at the egg stage, and nest-placement characteristics were important at the nestling stage. Although predator control appeared to have a positive influence on Blackbird breeding population size, the non-experimental set-up meant we could not eliminate other potential explanations. Variation in breeding population size did not appear to be influenced by variation in nest success alone. Our study demonstrates that observational data can only go so far in detection of such effects, and we discuss how it might be taken further. Agri-environment and game-management techniques are likely to influence nest predation pressure on farmland passerines, but the patterns, mechanisms and importance to population processes remain not wholly understood.

Functional benefits of predator species diversity depend on prey identity

1. Determining the functional significance of species diversity in natural enemy assemblages is a key step towards prediction of the likely impact of biodiversity loss on natural pest control processes. While the biological control literature contains examples in which increased natural enemy diversity hinders pest control, other studies have highlighted mechanisms where pest suppression is promoted by increased enemy diversity. 2. This study aimed to test whether increased predator species diversity results in higher rates of predation on two key, but contrasting, insect pest species commonly found in the rice ecosystems of south-east Asia. 3. Glasshouse experiments were undertaken in which four life stages of a planthopper (Nilaparvata lugens) and a moth (Marasmia patnalis) were caged with single or three-species combinations of generalist predators. 4. Generally, predation rates of the three-species assemblages exceeded expectation when attacking M. patnalis, but not when attacking N. lugens. In addition, a positive effect of increased predator species richness on overall predation rate was found with M. patnalis but not with N. lugens. 5. The results are consistent with theoretical predictions that morphological and behavioural differentiation among prey life stages promotes functional complementarity among predator species. This indicates that emergent species diversity effects in natural enemy assemblages are context dependent; they depend not only on the characteristics of the predators species, but on the identity of the species on which they prey.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.