Validation of PigBal model predictions for pig manure production

PigBal is a mass balance model that uses pig diet, digestibility and production data to predict the manure solids and nutrients produced by pig herds. It has been widely used for designing piggery effluent treatment systems and sustainable reuse areas at Australian piggeries. More recently, PigBal has also been used to estimate piggery volatile solids production for assessing greenhouse gas emissions for statutory reporting purposes by government, and for evaluating the energy potential from anaerobic digestion of pig effluent. This paper has compared PigBal predictions of manure total, volatile, and fixed solids, and nitrogen (N), phosphorus (P) and potassium (K), with manure production data generated in a replicated trial, which involved collecting manure from pigs housed in metabolic pens. Predictions of total, volatile, and fixed solids and K in the excreted manure were relatively good (combined diet R2 ≥ 0.79, modelling efficiency (EF) ≥ 0.70) whereas predictions of N and P, were generally less accurate (combined diet R2 0.56 and 0.66, EF 0.19 and –0.22, respectively). PigBal generally under-predicted lower N values while over-predicting higher values, and generally over-predicted manure P production for all diets. The most likely causes for this less accurate performance were ammonium-N volatilisation losses between manure excretion and sample analysis, and the inability of PigBal to account for higher rates of P uptake by pigs fed diets containing phytase. The outcomes of this research suggest that there is a need for further investigation and model development to enhance PigBal’s capabilities for more accurately assessing nutrient loads. However, PigBal’s satisfactory performance in predicting solids excretion demonstrates that it is suitable for assessing the methane component of greenhouse gas emission and the energy potential from anaerobic digestion of volatile solids in piggery effluent. The apparent overestimation of N and P excretion may result in conservative nutrient application rates to land and the over-prediction of the nitrous oxide component of greenhouse gas emissions.

Validation of PigBal model predictions for pig manure production

PigBal is a mass balance model that uses pig diet, digestibility and production data to predict the manure solids and nutrients produced by pig herds. It has been widely used for designing piggery effluent treatment systems and sustainable reuse areas at Australian piggeries. More recently, PigBal has also been used to estimate piggery volatile solids production for assessing greenhouse gas emissions for statutory reporting purposes by government, and for evaluating the energy potential from anaerobic digestion of pig effluent. This paper has compared PigBal predictions of manure total, volatile, and fixed solids, and nitrogen (N), phosphorus (P) and potassium (K), with manure production data generated in a replicated trial, which involved collecting manure from pigs housed in metabolic pens. Predictions of total, volatile, and fixed solids and K in the excreted manure were relatively good (combined diet R2 ≥ 0.79, modelling efficiency (EF) ≥ 0.70) whereas predictions of N and P, were generally less accurate (combined diet R2 0.56 and 0.66, EF 0.19 and -0.22, respectively). PigBal generally under-predicted lower N values while over-predicting higher values, and generally over-predicted manure P production for all diets. The most likely causes for this less accurate performance were ammonium-N volatilisation losses between manure excretion and sample analysis, and the inability of PigBal to account for higher rates of P uptake by pigs fed diets containing phytase. The outcomes of this research suggest that there is a need for further investigation and model development to enhance PigBal's capabilities for more accurately assessing nutrient loads. However, PigBal's satisfactory performance in predicting solids excretion demonstrates that it is suitable for assessing the methane component of greenhouse gas emission and the energy potential from anaerobic digestion of volatile solids in piggery effluent. The apparent overestimation of N and P excretion may result in conservative nutrient application rates to land and the over-prediction of the nitrous oxide component of greenhouse gas emissions. © CSIRO 2016.

Predicting mortality of phosphine-resistant adults of Sitophilus oryzae (L) (Coleoptera: Curculionidae) in relation to changing phosphine concentration

Adults of a phosphine-resistant strain of Sitophilus oryzae (L) were exposed to constant phosphine concentrations of 0.0035-0.9 mg litre(-1) for periods of between 20 and 168 h at 25 °C, and the effects of time and concentration on mortality were quantified. Adults were also exposed to a series of treatments lasting 48, 72 or 168 h at 25 °C, during which the concentration of phosphine was varied. The aim of this study was to determine whether equations from experiments using constant concentrations could be used to predict the efficacy of changing phosphine concentrations against adults of S oryzae. A probit plane without interaction, in which the logarithms of time (t) and concentration (C) were variables, described the effects of concentration and time on mortality in experiments with constant concentrations. A derived equation of the form C^nt = k gave excellent predictions of toxicity when applied to data from changing concentration experiments. The results suggest that for resistant S oryzae adults there is nothing inherently different between constant and changing concentration regimes, and that data collected from fixed concentrations can be used to develop equations for predicting mortality in fumigations in which phosphine concentration changes. This approach could simplify the prediction of efficacy of typical fumigations in which concentrations tend to rise and then fall over a period of days.

Effect of seasonal conditions on fleece weight and its components in merino sheep grazing Mitchell grass pasture in North West Queensland

SUMMARY Seasonal conditions in the pre to post natal period and selected periods before and during wool growth were described using climatic measures and estimates of the quality and quantity of pasture on offer derived from a validated pasture production model (GRASP). The variation in greasy and clean fleece weight, yield, staple length, fibre diameter, neck and side wrinkle score of Merinos grazing Mitchell grass in north west Queensland was explained in terms of these pasture and climatic measures and animal characteristics such as reproductive status, age and skin area. Multiple regression equations predicting clean and greasy fleece weight from the proportion of days in the wool growth period that the green pool in the pasture was less than one kg/ha, the percentage utilisation of the pasture, age, reproductive status and skin area of the ewes explained 87% and 79% of the variation respectively. Equations with similar predictors explained 58-85% of the variation of the other components. The inclusion of pasture conditions in the pre to post natal period did not significantly improve the predictions of the animal’s later performance. 22nd Biennial Conference.

Modelling Leaf Production and Crop Development in Maize (Zea mays L.) After Tassel Initiation Under Diverse Conditions of Temperature and Photoperiod

Prediction of the initiation, appearance and emergence of leaves is critically important to the success of simulation models of crop canopy development and some aspects of crop ontogeny. Data on leaf number and crop ontogeny were collected on five cultivars of maize differing widely in maturity and genetic background grown under natural and extended photoperiods, and planted on seven sowing dates from October 1993 to March 1994 at Gatton, South-east Queensland. The same temperature coefficients were established for crop ontogeny before silking, and the rates of leaf initiation, leaf tip appearance and full leaf expansion, the base, optimum and maximum temperatures for each being 8, 34 and 40 degrees C. After silking, the base temperature for ontogeny was 0 degrees C, but the optimum and maximum temperatures remained unchanged. The rates of leaf initiation, appearance of leaf tips and full leaf expansion varied in a relatively narrow range across sowing times and photoperiod treatments, with average values of 0.040 leaves (degrees Cd)-1, 0.021 leaves (degrees Cd)-1, and 0.019 leaves (degrees Cd)-1, respectively. The relationships developed in this study provided satisfactory predictions of leaf number and crop ontogeny (tassel initiation to silking, emergence to silking and silking to physiological maturity) when assessed using independent data from Gatton (South eastern Queensland), Katherine and Douglas Daly (Northern Territory), Walkamin (North Queensland) and Kununurra (Western Australia).

Diet quality and liveweight gain of steers grazing Leucaena-grass pasture estimated with faecal near infrared reflectance spectroscopy (F.NIRS)

Three drafts of Bos indicus cross steers (initially 178-216 kg) grazed Leucaena-grass pasture [Leucaena leucocephala subspecies glabrata cv. Cunningham with green panic (Panicum maximum cv. trichoglume)] from late winter through to autumn during three consecutive years in the Burnett region of south-east Queensland. Measured daily weight gain (DWGActual) of the steers was generally 0.7-1.1 kg/day during the summer months. Estimated intakes of metabolisable energy and dry matter (DM) were calculated from feeding standards as the intakes required by the steers to grow at the DWGActual. Diet attributes were predicted from near infrared reflectance spectroscopy spectra of faeces (F.NIRS) using established calibration equations appropriate for northern Australian forages. Inclusion of some additional reference samples from cattle consuming Leucaena diets into F.NIRS calibrations based on grass and herbaceous legume-grass pastures improved prediction of the proportion of Leucaena in the diet. Mahalanobis distance values supported the hypothesis that the F.NIRS predictions of diet crude protein concentration and DM digestibility (DMD) were acceptable. F.NIRS indicated that the percentage of Leucaena in the diet varied widely (10-99%). Diet crude protein concentration and DMD were usually high, averaging 12.4 and 62%, respectively, and were related asymptotically to the percentage of Leucaena in the diet (R2 = 0.48 and 0.33, respectively). F.NIRS calibrations for DWG were not satisfactory to predict this variable from an individual faecal sample since the s.e. of prediction were 0.33-0.40 kg/day. Cumulative steer liveweight (LW) predicted from F.NIRS DWG calibrations, which had been previously developed with tropical grass and grass-herbaceous legume pastures, greatly overestimated the measured steer LW; therefore, these calibrations were not useful. Cumulative steer LW predicted from a modified F.NIRS DWG calibration, which included data from the present study, was strongly correlated (R2 = 0.95) with steer LW but overestimated LW by 19-31 kg after 8 months. Additional reference data are needed to develop robust F.NIRS calibrations to encompass the diversity of Leucaena pastures of northern Australia. In conclusion, the experiment demonstrated that F.NIRS could improve understanding of diet quality and nutrient intake of cattle grazing Leucaena-grass pasture, and the relationships between nutrient supply and cattle growth.

NIR model development and robustness in prediction of melon fruit total soluble solids

Near infrared spectroscopy (NIRS) can be used for the on-line, non-invasive assessment of fruit for eating quality attributes such as total soluble solids (TSS). The robustness of multivariate calibration models, based on NIRS in a partial transmittance optical geometry, for the assessment of TSS of intact rockmelons (Cucumis melo) was assessed. The mesocarp TSS was highest around the fruit equator and increased towards the seed cavity. Inner mesocarp TSS levels decreased towards both the proximal and distal ends of the fruit, but more so towards the proximal end. The equatorial region of the fruit was chosen as representative of the fruit for near infrared assessment of TSS. The spectral window for model development was optimised at 695-1045 nm, and the data pre-treatment procedure was optimised to second-derivative absorbance without scatter correction. The 'global' modified partial least squares (MPLS) regression modelling procedure of WINISI (ver. 1.04) was found to be superior with respect to root mean squared error of prediction (RMSEP) and bias for model predictions of TSS across seasons, compared with the 'local' MPLS regression procedure. Updating of the model with samples selected randomly from the independent validation population demonstrated improvement in both RMSEP and bias with addition of approximately 15 samples.

Review: Near infrared spectroscopy of faeces to evaluate the nutrition and physiology of herbivores

Near infrared (NIR) spectroscopy, usually in reflectance mode, has been applied to the analysis of faeces to measure the concentrations of constituents such as total N, fibre, tannins and delta C-13. In addition, an unusual and exciting application of faecal NIR [F.NIR] analyses is to directly predict attributes of the diet of herbivores such as crude protein and fibre contents, proportions of plant species and morphological components, diet digestibility and voluntary DM intake. This is an unusual application of NIR spectroscopy insofar as the spectral measurements are made, not on the material of interest [i.e. the diet), but on a derived material (i.e. faeces). Predictions of diet attributes from faecal spectra clearly depend on there being sufficient NIR spectral information in the diet residues present in faeces to describe the diet, although endogenous components of faeces such as undigested debris of micro-organisms from the rumen and Large intestine and secretions into the gastrointestinal tract wilt also contribute spectral information. Spectra of forage and of faeces derived from the forage are generally similar and the observed differences are principally in the spectral regions associated with constituents of forages known to be of low, or of high, digestibility. Some diet components (for example, ureal which are likely to be entirely digested apparently cannot be predicted from faecal NIR spectra because they cannot contribute to faecal spectra except through modifying the microbial and endogenous components. The errors and robustness of F.NIR calibrations to predict the crude protein concentration and digestibility of the diet of herbivores are generally comparable with those to directly predict the same attributes in forage from NIR spectra of the forage. Some attributes of the animal, such as species, gender, pregnancy status and parasite burden have been successfully discriminated into classes based on their faecal NIR spectra. Such discrimination was likely associated with differences in the diet selected and/or differences in the metabolites excreted in the faeces. NIR spectroscopy of faeces has usually involved scanning dried and ground samples in monochromators in the 400-2500nm or 1100-2500nm ranges. Results satisfactory for the purpose have also been reported for dried and ground faeces scanned using a diode array instrument in the 800-1700nm range and for wet faeces and slurries of excreta scanned with monochromators. Chemometric analysis of faecal spectra has generally used the approaches established for forage analysis. The capacity to predict many attributes of the diet, and some aspects of animal physiology, from NIR spectra of faeces is particularly useful to study the quality and quantity of the diet selected by both domestic and feral grazing herbivores and to enhance production and management of both herbivores and their grazing environment.

Seasonal changes in pasture quality and diet selection and their relationship with liveweight gain of steers grazing tropical grass and grass-legume pastures in northern Australia.

The variation in liveweight gain in grazing beef cattle as influenced by pasture type, season and year effects has important economic implications for mixed crop-livestock systems and the ability to better predict such variation would benefit beef producers by providing a guide for decision making. To identify key determinants of liveweight change of Brahman-cross steers grazing subtropical pastures, measurements of pasture quality and quantity, and diet quality in parallel with liveweight were made over two consecutive grazing seasons (48 and 46 weeks, respectively), on mixed Clitoria ternatea/grass, Stylosanthes seabrana/grass and grass swards (grass being a mixture of Bothriochloa insculpta cv. Bisset, Dichanthium sericeum and Panicum maximum var. trichoglume cv. Petrie). Steers grazing the legume-based pastures had the highest growth rate and gained between 64 and 142 kg more than those grazing the grass pastures in under 12 months. Using an exponential model, green leaf mass, green leaf %, adjusted green leaf % (adjusted for inedible woody legume stems), faecal near infrared reflectance spectroscopy predictions of diet crude protein and diet dry matter digestibility, accounted for 77, 74, 80, 63 and 60%, respectively, of the variation in daily weight gain when data were pooled across pasture types and grazing seasons. The standard error of the regressions indicated that 95% prediction intervals were large (+/- 0.42-0.64 kg/head.day) suggesting that derived regression relationships have limited practical application for accurately estimating growth rate. In this study, animal factors, especially compensatory growth effects, appeared to have a major influence on growth rate in relation to pasture and diet attributes. It was concluded that predictions of growth rate based only on pasture or diet attributes are unlikely to be accurate or reliable. Nevertheless, key pasture attributes such as green leaf mass and green leaf% provide a robust indication of what proportion of the potential growth rate of the grazing animals can be achieved.

The role of fruit traits of bird-dispersed plants in invasiveness and weed risk assessment

Aim: Birds play a major role in the dispersal of seeds of many fleshy-fruited invasive plants. The fruits that birds choose to consume are influenced by fruit traits. However, little is known of how the traits of invasive plant fruits contribute to invasiveness or to their use by frugivores. We aim to gain a greater understanding of these relationships to improve invasive plant management. Location: South-east Queensland, Australia. Methods: We measure a variety of fruit morphology, pulp nutrient and phenology traits of a suite of bird-dispersed alien plants. Frugivore richness of these aliens was derived from the literature. Using regressions and multivariate methods, we investigate relationships between fruit traits, frugivore richness and invasiveness. Results: Plant invasiveness was negatively correlated to fruit size, and all highly invasive species had quite similar fruit morphology [smaller fruits, seeds of intermediate size and few (<10) seeds per fruit]. Lower pulp water was the only pulp nutrient trait associated with invasiveness. There were strong positive relationships between the diversity of bird frugivores and plant invasiveness, and in the diversity of bird frugivores in the study region and another part of the plants' alien range. Main conclusions: Our results suggest that weed risk assessments (WRA) and predictions of invasive success for bird-dispersed plants can be improved. Scoring criteria for WRA regarding fruit size would need to be system-specific, depending on the fruit-processing capabilities of local frugivores. Frugivore richness could be quantified in the plant's natural range, its invasive range elsewhere, or predictions made based on functionally similar fruits.

Analysis and Modelling of the Effects of Water Stress on Maize Growth and Yield in Dryland Conditions.

It is essential to provide experimental evidence and reliable predictions of the effects of water stress on crop production in the drier, less predictable environments. A field experiment undertaken in southeast Queensland, Australia with three water regimes (fully irrigated, rainfed and irrigated until late canopy expansion followed by rainfed) was used to compare effects of water stress on crop production in two maize (Zea mays L.) cultivars (Pioneer 34N43 and Pioneer 31H50). Water stress affected growth and yield more in Pioneer 34N43 than in Pioneer 31H50. A crop model APSIM-Maize, after having been calibrated for the two cultivars, was used to simulate maize growth and development under water stress. The predictions on leaf area index (LAI) dynamics, biomass growth and grain yield under rain fed and irrigated followed by rain fed treatments was reasonable, indicating that stress indices used by APSIM-Maize produced appropriate adjustments to crop growth and development in response to water stress. This study shows that Pioneer 31H50 is less sensitive to water stress and thus a preferred cultivar in dryland conditions, and that it is feasible to provide sound predictions and risk assessment for crop production in drier, more variable conditions using the APSIM-Maize model.

The effect of data sources and quality on the predictive capacity of CLIMEX models: an assessment of Teleonemia scrupulosa and Octotoma scabripennis for the biocontrol of Lantana camara in Australia

Understanding the effects of different types and quality of data on bioclimatic modeling predictions is vital to ascertaining the value of existing models, and to improving future models. Bioclimatic models were constructed using the CLIMEX program, using different data types – seasonal dynamics, geographic (overseas) distribution, and a combination of the two – for two biological control agents for the major weed Lantana camara L. in Australia. The models for one agent, Teleonemia scrupulosa Stål (Hemiptera:Tingidae) were based on a higher quality and quantity of data than the models for the other agent, Octotoma scabripennis Guérin-Méneville (Coleoptera: Chrysomelidae). Predictions of the geographic distribution for Australia showed that T. scrupulosa models exhibited greater accuracy with a progressive improvement from seasonal dynamics data, to the model based on overseas distribution, and finally the model combining the two data types. In contrast, O. scabripennis models were of low accuracy, and showed no clear trends across the various model types. These case studies demonstrate the importance of high quality data for developing models, and of supplementing distributional data with species seasonal dynamics data wherever possible. Seasonal dynamics data allows the modeller to focus on the species response to climatic trends, while distributional data enables easier fitting of stress parameters by restricting the species envelope to the described distribution. It is apparent that CLIMEX models based on low quality seasonal dynamics data, together with a small quantity of distributional data, are of minimal value in predicting the spatial extent of species distribution.

Development of a productivity asessment toll for native spotted gum forest on private land based on estimates of forest growth on Crown land.

Reliable estimates of forest productivity are essential for improved predictions of timber yields for the private native spotted gum resource in southern Qld and northern NSW. The aim of this research was to estimate the potential productivity of native spotted gum forests on private land by making use of available inventory data collated from Qld and northern NSW for spotted gum forest on Crown land (i.e. state forests). We measured a range of site-related factors to determine their relative importance in predicting productivity of spotted gum forest. While measures such as stand height and height-diameter relationships are known to be useful predictors of productivity, we aimed to determine productivity for a site where this information was not available. Through estimation of stand growth rates we developed a spotted gum productivity assessment tool (SPAT) for use by landholders and extension officers. We aimed to develop a tool to allow private landholders to see the benefits of maintaining their timber resource. This paper summarises the information used to develop the SPAT with a particular focus on forest growth relationships.

Potential Distribution and Risk Assessment of an Invasive Plant Species: A Case Study of Hymenachne amplexicaulis in Australia

Given the limited resources available for weed management, a strategic approach is required to give the best bang for your buck. The current study incorporates: (1) a model ensemble approach to identify areas of uncertainty and commonality regarding a species invasive potential, (2) current distribution of the invaded species, and (3) connectivity of systems to identify target regions and focus efforts for more effective management. Uncertainty in the prediction of suitable habitat for H. amplexicaulis (study species) in Australia was addressed in an ensemble-forecasting approach to compare distributional scenarios from four models (CLIMATCH; CLIMEX; boosted regression trees [BRT]; maximum entropy [Maxent]). Models were built using subsets of occurrence and environmental data. Catchment risk was determined through incorporating habitat suitability, the current abundance and distribution of H. amplexicaulis, and catchment connectivity. Our results indicate geographic differences between predictions of different approaches. Despite these differences a number of catchments in northern, central, and southern Australia were identified as high risk of invasion or further spread by all models suggesting they should be given priority for the management of H. amplexicaulis. The study also highlighted the utility of ensemble approaches in indentifying areas of uncertainty and commonality regarding the species invasive potential.

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.