Constructing a new individual-based model of phosphine resistance in lesser grain borer (Rhyzopertha dominica): do we need to include two loci rather than one?

In this article, we describe and compare two individual-based models constructed to investigate how genetic factors influence the development of phosphine resistance in lesser grain borer (R. dominica). One model is based on the simplifying assumption that resistance is conferred by alleles at a single locus, while the other is based on the more realistic assumption that resistance is conferred by alleles at two separate loci. We simulated the population dynamic of R. dominica in the absence of phosphine fumigation, and under high and low dose phosphine treatments, and found important differences between the predictions of the two models in all three cases. In the absence of fumigation, starting from the same initial frequencies of genotypes, the two models tended to different stable frequencies, although both reached Hardy-Weinberg equilibrium. The one-locus model exaggerated the equilibrium proportion of strongly resistant beetles by 3.6 times, compared to the aggregated predictions of the two-locus model. Under a low dose treatment the one-locus model overestimated the proportion of strongly resistant individuals within the population and underestimated the total population numbers compared to the two-locus model. These results show the importance of basing resistance evolution models on realistic genetics and that using oversimplified one-locus models to develop pest control strategies runs the risk of not correctly identifying tactics to minimise the incidence of pest infestation.

APSIM – Evolution towards a new generation of agricultural systems simulation

Agricultural systems models worldwide are increasingly being used to explore options and solutions for the food security, climate change adaptation and mitigation and carbon trading problem domains. APSIM (Agricultural Production Systems sIMulator) is one such model that continues to be applied and adapted to this challenging research agenda. From its inception twenty years ago, APSIM has evolved into a framework containing many of the key models required to explore changes in agricultural landscapes with capability ranging from simulation of gene expression through to multi-field farms and beyond. Keating et al. (2003) described many of the fundamental attributes of APSIM in detail. Much has changed in the last decade, and the APSIM community has been exploring novel scientific domains and utilising software developments in social media, web and mobile applications to provide simulation tools adapted to new demands. This paper updates the earlier work by Keating et al. (2003) and chronicles the changing external challenges and opportunities being placed on APSIM during the last decade. It also explores and discusses how APSIM has been evolving to a “next generation” framework with improved features and capabilities that allow its use in many diverse topics.

Estimating daily methane production in individual cattle with irregular feed intake patterns from short-term methane emission measurements

Spot measurements of methane emission rate (n = 18 700) by 24 Angus steers fed mixed rations from GrowSafe feeders were made over 3- to 6-min periods by a GreenFeed emission monitoring (GEM) unit. The data were analysed to estimate daily methane production (DMP; g/day) and derived methane yield (MY; g/kg dry matter intake (DMI)). A one-compartment dose model of spot emission rate v. time since the preceding meal was compared with the models of Wood (1967) and Dijkstra et al. (1997) and the average of spot measures. Fitted values for DMP were calculated from the area under the curves. Two methods of relating methane and feed intakes were then studied: the classical calculation of MY as DMP/DMI (kg/day); and a novel method of estimating DMP from time and size of preceding meals using either the data for only the two meals preceding a spot measurement, or all meals for 3 days prior. Two approaches were also used to estimate DMP from spot measurements: fitting of splines on a 'per-animal per-day' basis and an alternate approach of modelling DMP after each feed event by least squares (using Solver), summing (for each animal) the contributions from each feed event by best-fitting a one-compartment model. Time since the preceding meal was of limited value in estimating DMP. Even when the meal sizes and time intervals between a spot measurement and all feeding events in the previous 72 h were assessed, only 16.9% of the variance in spot emission rate measured by GEM was explained by this feeding information. While using the preceding meal alone gave a biased (underestimate) of DMP, allowing for a longer feed history removed this bias. A power analysis taking into account the sources of variation in DMP indicated that to obtain an estimate of DMP with a 95% confidence interval within 5% of the observed 64 days mean of spot measures would require 40 animals measured over 45 days (two spot measurements per day) or 30 animals measured over 55 days. These numbers suggest that spot measurements could be made in association with feed efficiency tests made over 70 days. Spot measurements of enteric emissions can be used to define DMP but the number of animals and samples are larger than are needed when day-long measures are made.

Postharvest treatment effects on 'B74' mango fruit lenticel discolouration after irradiation

Lenticel discolouration (LD) is a common disorder of mango fruit around the world. It results in poor appearance and disappointment of consumers. LD is exacerbated by treatment of mango fruit with gamma irradiation for insect disinfestation. The issue is problematic on the relatively new mango cultivar 'B74' and may represent an oxidative browning process. With a view to reducing irradiationinduced LD on 'B74', postharvest wax (one and three layers; 75% carnauba wax) and antioxidant (100 mM ascorbic acid, 100 mM calcium chloride, 10, 50 and 100 mM calcium ascorbate) dip treatments were investigated. Treatment of green mature fruit with three layers of wax prior to exposure to 557 Gy gamma irradiation reduced LD by 40% relative to the non-waxed control. However, the fruit failed to ripen properly as evidenced by delayed skin colour change, retarded softening and increased skin browning as compared to the controls and fruit coated with one layer of wax. Treatment with one layer of wax did not reduce LD. Mechanistically, the responses suggest that air exchange plays a pivotal role in LD. A lowered oxygen concentration in the lenticels may reduce the disorder after irradiation treatment. Postharvest treatments with the various antioxidants failed to reduce LD. Rather, all antioxidant treatments at the test concentrations, except calcium chloride, significantly increased skin browning.

Spatial model of site index based on Y-ray spectrometry and a digital elevation model for two Pinus species in Tuan Toolara State Forest, Queensland, Australia

Site index prediction models are an important aid for forest management and planning activities. This paper introduces a multiple regression model for spatially mapping and comparing site indices for two Pinus species (Pinus elliottii Engelm. and Queensland hybrid, a P. elliottii x Pinus caribaea Morelet hybrid) based on independent variables derived from two major sources: g-ray spectrometry (potassium (K), thorium (Th), and uranium (U)) and a digital elevation model (elevation, slope, curvature, hillshade, flow accumulation, and distance to streams). In addition, interpolated rainfall was tested. Species were coded as a dichotomous dummy variable; interaction effects between species and the g-ray spectrometric and geomorphologic variables were considered. The model explained up to 60% of the variance of site index and the standard error of estimate was 1.9 m. Uranium, elevation, distance to streams, thorium, and flow accumulation significantly correlate to the spatial variation of the site index of both species, and hillshade, curvature, elevation and slope accounted for the extra variability of one species over the other. The predicted site indices varied between 20.0 and 27.3 m for P. elliottii, and between 23.1 and 33.1 m for Queensland hybrid; the advantage of Queensland hybrid over P. elliottii ranged from 1.8 to 6.8 m, with the mean at 4.0 m. This compartment-based prediction and comparison study provides not only an overview of forest productivity of the whole plantation area studied but also a management tool at compartment scale.

Net radiation dynamics: performance of 20 daily net radiation models as related to model structure and intricacy in two climates.

We compared daily net radiation (Rn) estimates from 19 methods with the ASCE-EWRI Rn estimates in two climates: Clay Center, Nebraska (sub-humid) and Davis, California (semi-arid) for the calendar year. The performances of all 20 methods, including the ASCE-EWRI Rn method, were then evaluated against Rn data measured over a non-stressed maize canopy during two growing seasons in 2005 and 2006 at Clay Center. Methods differ in terms of inputs, structure, and equation intricacy. Most methods differ in estimating the cloudiness factor, emissivity (e), and calculating net longwave radiation (Rnl). All methods use albedo (a) of 0.23 for a reference grass/alfalfa surface. When comparing the performance of all 20 Rn methods with measured Rn, we hypothesized that the a values for grass/alfalfa and non-stressed maize canopy were similar enough to only cause minor differences in Rn and grass- and alfalfa-reference evapotranspiration (ETo and ETr) estimates. The measured seasonal average a for the maize canopy was 0.19 in both years. Using a = 0.19 instead of a = 0.23 resulted in 6% overestimation of Rn. Using a = 0.19 instead of a = 0.23 for ETo and ETr estimations, the 6% difference in Rn translated to only 4% and 3% differences in ETo and ETr, respectively, supporting the validity of our hypothesis. Most methods had good correlations with the ASCE-EWRI Rn (r2 > 0.95). The root mean square difference (RMSD) was less than 2 MJ m-2 d-1 between 12 methods and the ASCE-EWRI Rn at Clay Center and between 14 methods and the ASCE-EWRI Rn at Davis. The performance of some methods showed variations between the two climates. In general, r2 values were higher for the semi-arid climate than for the sub-humid climate. Methods that use dynamic e as a function of mean air temperature performed better in both climates than those that calculate e using actual vapor pressure. The ASCE-EWRI-estimated Rn values had one of the best agreements with the measured Rn (r2 = 0.93, RMSD = 1.44 MJ m-2 d-1), and estimates were within 7% of the measured Rn. The Rn estimates from six methods, including the ASCE-EWRI, were not significantly different from measured Rn. Most methods underestimated measured Rn by 6% to 23%. Some of the differences between measured and estimated Rn were attributed to the poor estimation of Rnl. We conducted sensitivity analyses to evaluate the effect of Rnl on Rn, ETo, and ETr. The Rnl effect on Rn was linear and strong, but its effect on ETo and ETr was subsidiary. Results suggest that the Rn data measured over green vegetation (e.g., irrigated maize canopy) can be an alternative Rn data source for ET estimations when measured Rn data over the reference surface are not available. In the absence of measured Rn, another alternative would be using one of the Rn models that we analyzed when all the input variables are not available to solve the ASCE-EWRI Rn equation. Our results can be used to provide practical information on which method to select based on data availability for reliable estimates of daily Rn in climates similar to Clay Center and Davis.

Simple software processes and tests improve the reliability and usefulness of a model.

Models are abstractions of reality that have predetermined limits (often not consciously thought through) on what problem domains the models can be used to explore. These limits are determined by the range of observed data used to construct and validate the model. However, it is important to remember that operating the model beyond these limits, one of the reasons for building the model in the first place, potentially brings unwanted behaviour and thus reduces the usefulness of the model. Our experience with the Agricultural Production Systems Simulator (APSIM), a farming systems model, has led us to adapt techniques from the disciplines of modelling and software development to create a model development process. This process is simple, easy to follow, and brings a much higher level of stability to the development effort, which then delivers a much more useful model. A major part of the process relies on having a range of detailed model tests (unit, simulation, sensibility, validation) that exercise a model at various levels (sub-model, model and simulation). To underline the usefulness of testing, we examine several case studies where simulated output can be compared with simple relationships. For example, output is compared with crop water use efficiency relationships gleaned from the literature to check that the model reproduces the expected function. Similarly, another case study attempts to reproduce generalised hydrological relationships found in the literature. This paper then describes a simple model development process (using version control, automated testing and differencing tools), that will enhance the reliability and usefulness of a model.

Antibody Responses to Sarcoptes scabiei Apolipoprotein in a Porcine Model: Relevance to Immunodiagnosis of Recent Infection

No commercial immunodiagnostic tests for human scabies are currently available, and existing animal tests are not sufficiently sensitive. The recombinant Sarcoptes scabiei apolipoprotein antigen Sar s 14.3 is a promising immunodiagnostic, eliciting high levels of IgE and IgG in infected people. Limited data are available regarding the temporal development of antibodies to Sar s 14.3, an issue of relevance in terms of immunodiagnosis. We utilised a porcine model to prospectively compare specific antibody responses to a primary infestation by ELISA, to Sar s 14.3 and to S. scabiei whole mite antigen extract (WMA). Differences in the antibody profile between antigens were apparent, with Sar s 14.3 responses detected earlier, and declining significantly after peak infestation compared to WMA. Both antigens resulted in >90% diagnostic sensitivity from weeks 8–16 post infestation. These data provide important information on the temporal development of humoral immune responses in scabies and further supports the development of recombinant antigen based immunodiagnostic tests for recent scabies infestations.

Stochastic model for simulating maize yield

Maize is one of the most important crops in the world. The products generated from this crop are largely used in the starch industry, the animal and human nutrition sector, and biomass energy production and refineries. For these reasons, there is much interest in figuring the potential grain yield of maize genotypes in relation to the environment in which they will be grown, as the productivity directly affects agribusiness or farm profitability. Questions like these can be investigated with ecophysiological crop models, which can be organized according to different philosophies and structures. The main objective of this work is to conceptualize a stochastic model for predicting maize grain yield and productivity under different conditions of water supply while considering the uncertainties of daily climate data. Therefore, one focus is to explain the model construction in detail, and the other is to present some results in light of the philosophy adopted. A deterministic model was built as the basis for the stochastic model. The former performed well in terms of the curve shape of the above-ground dry matter over time as well as the grain yield under full and moderate water deficit conditions. Through the use of a triangular distribution for the harvest index and a bivariate normal distribution of the averaged daily solar radiation and air temperature, the stochastic model satisfactorily simulated grain productivity, i.e., it was found that 10,604 kg ha(-1) is the most likely grain productivity, very similar to the productivity simulated by the deterministic model and for the real conditions based on a field experiment. © 2012 American Society of Agricultural and Biological Engineers.

Odour, dust and non-methane volatile organic-compound emissions from tunnel-ventilated layer-chicken sheds: a case study of two farms

An observational study was undertaken to measure odour and dust (PM10 and PM2.5) emission rates and identify non-methane volatile organic compounds (NMVOCs) and odorants in the exhaust air from two tunnel-ventilated layer-chicken sheds that were configured with multi-tiered cages and manure belts. The study sites were located in south-eastern Queensland and the West Gippsland region of Victoria, Australia. Samples were collected in summer and winter on sequential days across the manure-belt cleaning cycle. Odour emissions ranged from 58 to 512 ou/s per 1000 birds (0.03-0.27 ou/s.kg) and dust emission rates ranged 0.014-0.184 mg/s per 1000 birds for PM10 and 0.001-0.190 mg/s per 1000 birds for PM2.5. Twenty NMVOCs were identified, including three that were also identified as odorants using thermal desorption-gas chromatography-mass spectrometry/olfactometry analysis. Odour emission rates were observed to vary with the amount of manure accumulation on the manure belts, being lowest 2-4 days after removing manure. Odour emission rates were also observed to vary with diurnal and seasonal changes in ventilation rate. Dust emissions were observed to increase with ventilation rate but not with manure accumulation. Some NMVOCs were identified at both farms and in different seasons whereas others were observed only at one farm or in one season, indicating that odorant composition was influenced by farm-specific practices and season.

Evaluation and histological examination of a Campylobacter fetus subsp. venerealis small animal infection model

Bovine genital campylobacteriosis (BGC), caused by Campylobacter fetus subsp. venerealis, is associated with production losses in cattle worldwide. This study aimed to develop a reliable BGC guinea pig model to facilitate future studies of pathogenicity, abortion mechanisms and vaccine efficacy. Seven groups of five pregnant guinea pigs (1 control per group) were inoculated with one of three strains via intra-peritoneal (IP) or intra-vaginal routes. Samples were examined using culture, PCR and histology. Abortions ranged from 0 to 100 and re-isolation of causative bacteria from sampled sites varied with strain, dose of bacteria and time to abortion. Histology indicated metritis and placentitis, suggesting that the bacteria induce inflammation, placental detachment and subsequent abortion. Variation of virulence between strains was observed and determined by culture and abortion rates. IP administration of C. fetus subsp. venerealis to pregnant guinea pigs is a promising small animal model for the investigation of BGC abortion.