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.

Prospects for predicting insect mortality in relation to changing phosphine concentrations

Typically, in bag-stack or silo fumigations the concentration of phosphine is not constant, and yet most of what is known about phosphine efficacy against grain insects comes from studies with fixed concentrations. Indeed, where changing concentration experiments have been performed, researchers have been unable to explain observed efficacy on the basis of data from fixed concentrations. The ability to predict insect mortality in relation to changing phosphine concentrations would facilitate the development of effective fumigation protocols. In this paper, we explore the prospects for making such predictions. After reviewing published and new results, we conclude that the commonly used concentration x time (Ct) product is unreliable for this purpose. New results, for a strongly resistant strain of Rhyzopertha dominica from Australia, suggest that the relationship Cnt = k may be useful for predicting mortality of this type of insect in changing concentrations. However, in the case of a strain of Sitophilus oryzae with a type of resistance common in Australian S. oryzae, the relationship Cnt = k proved to be less reliable.

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.

Runoff, sediment loss and water quality from forest roads in a southeast Queensland coastal plain Pinus plantation

Runoff and sediment loss from forest roads were monitored for a two-year period in a Pinus plantation in southeast Queensland. Two classes of road were investigated: a gravelled road, which is used as a primary daily haulage route for the logging area, and an ungravelled road, which provides the main access route for individual logging compartments and is intensively used as a haulage route only during the harvest of these areas (approximately every 30 years). Both roads were subjected to routine traffic loads and maintenance during the study. Surface runoff in response to natural rainfall was measured and samples taken for the determination of sediment and nutrient (total nitrogen, total phosphorus, dissolved organic carbon and total iron) loads from each road. Results revealed that the mean runoff coefficient (runoff depth/rainfall depth) was consistently higher from the gravelled road plot with 0.57, as compared to the ungravelled road with 0.38. Total sediment loss over the two-year period was greatest from the gravelled road plot at 5.7 t km−1 compared to the ungravelled road plot with 3.9 t km−1. Suspended solids contributed 86% of the total sediment loss from the gravelled road, and 72% from the ungravelled road over the two years. Nitrogen loads from the two roads were both relatively constant throughout the study, and averaged 5.2 and 2.9 kg km−1 from the gravelled and ungravelled road, respectively. Mean annual phosphorus loads were 0.6 kg km−1 from the gravelled road and 0.2 kg km−1 from the ungravelled road. Organic carbon and total iron loads increased in the second year of the study, which was a much wetter year, and are thought to reflect the breakdown of organic matter in roadside drains and increased sediment generation, respectively. When road and drain maintenance (grading) was performed runoff and sediment loss were increased from both road types. Additionally, the breakdown of the gravel road base due to high traffic intensity during wet conditions resulted in the formation of deep (10 cm) ruts which increased erosion. The Water Erosion Prediction Project (WEPP):Road model was used to compare predicted to observed runoff and sediment loss from the two road classes investigated. For individual rainfall events, WEPP:Road predicted output showed strong agreement with observed values of runoff and sediment loss. WEPP:Road predictions for annual sediment loss from the entire forestry road network in the study area also showed reasonable agreement with the extrapolated observed values.

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.

Seed persistence in the field may be predicted by laboratory-controlled aging

Weed management is complicated by the presence of soil seed banks. The complexity of soil-seed interactions means that seed persistence in the field is often difficult to measure, let alone predict. Field trials, although accurate in their context, are time-consuming and expensive to conduct for individual species. Some ex situ techniques for estimating seed life expectancy have been proposed, but these fail to simulate the environmental complexity of the field. Also, it has been questioned whether techniques such as the controlled aging test (CAT) are useful indicators of field persistence. This study aimed to test the validity of the standard CAT (seed aging at 45 C and 60% relative humidity) in use at the Royal Botanic Gardens, Kew, U.K., for predicting field seed-persistence. Comparison of seed persistence and CAT data for 27 northwest European species suggested a significant positive correlation of 0.31. Subsequently, 13 species of emerging and common weeds of Queensland were assessed for their seed longevity using the CAT. The seed longevity data of these species in the CAT were linked with field seed-persistence data according to three broad seed-persistence categories: <1 yr, 1 to 3 yr, and >3 yr. We discuss the scope for using the CAT as a tool for rapid assignment of species to these categories. There is a need for further studies that compare predictions of seed persistence based on the CAT with seed persistence in the field for a larger range of species and environments.

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.

The accuracy of varietal selection using factor analytic models for multi-environment plant breeding trials

Modeling of cultivar x trial effects for multienvironment trials (METs) within a mixed model framework is now common practice in many plant breeding programs. The factor analytic (FA) model is a parsimonious form used to approximate the fully unstructured form of the genetic variance-covariance matrix in the model for MET data. In this study, we demonstrate that the FA model is generally the model of best fit across a range of data sets taken from early generation trials in a breeding program. In addition, we demonstrate the superiority of the FA model in achieving the most common aim of METs, namely the selection of superior genotypes. Selection is achieved using best linear unbiased predictions (BLUPs) of cultivar effects at each environment, considered either individually or as a weighted average across environments. In practice, empirical BLUPs (E-BLUPs) of cultivar effects must be used instead of BLUPs since variance parameters in the model must be estimated rather than assumed known. While the optimal properties of minimum mean squared error of prediction (MSEP) and maximum correlation between true and predicted effects possessed by BLUPs do not hold for E-BLUPs, a simulation study shows that E-BLUPs perform well in terms of MSEP.

Annual forecasting of the Australian macadamia crop - integrating tree census data with statistical climate-adjustment models

To facilitate marketing and export, the Australian macadamia industry requires accurate crop forecasts. Each year, two levels of crop predictions are produced for this industry. The first is an overall longer-term forecast based on tree census data of growers in the Australian Macadamia Society (AMS). This data set currently accounts for around 70% of total production, and is supplemented by our best estimates of non-AMS orchards. Given these total tree numbers, average yields per tree are needed to complete the long-term forecasts. Yields from regional variety trials were initially used, but were found to be consistently higher than the average yields that growers were obtaining. Hence, a statistical model was developed using growers' historical yields, also taken from the AMS database. This model accounted for the effects of tree age, variety, year, region and tree spacing, and explained 65% of the total variation in the yield per tree data. The second level of crop prediction is an annual climate adjustment of these overall long-term estimates, taking into account the expected effects on production of the previous year's climate. This adjustment is based on relative historical yields, measured as the percentage deviance between expected and actual production. The dominant climatic variables are observed temperature, evaporation, solar radiation and modelled water stress. Initially, a number of alternate statistical models showed good agreement within the historical data, with jack-knife cross-validation R2 values of 96% or better. However, forecasts varied quite widely between these alternate models. Exploratory multivariate analyses and nearest-neighbour methods were used to investigate these differences. For 2001-2003, the overall forecasts were in the right direction (when compared with the long-term expected values), but were over-estimates. In 2004 the forecast was well under the observed production, and in 2005 the revised models produced a forecast within 5.1% of the actual production. Over the first five years of forecasting, the absolute deviance for the climate-adjustment models averaged 10.1%, just outside the targeted objective of 10%.

Probabilistic forecasts of the onset of the north Australian wet season

The amount and timing of early wet-season rainfall are important for the management of many agricultural industries in north Australia. With this in mind, a wet-season onset date is defined based on the accumulation of rainfall to a predefined threshold, starting from 1 September, for each square of a 1° gridded analysis of daily rainfall across the region. Consistent with earlier studies, the interannual variability of the onset dates is shown to be well related to the immediately preceding July-August Southern Oscillation index (SOI). Based on this relationship, a forecast method using logistic regression is developed to predict the probability that onset will occur later than the climatological mean date. This method is expanded to also predict the probabilities that onset will be later than any of a range of threshold dates around the climatological mean. When assessed using cross-validated hindcasts, the skill of the predictions exceeds that of climatological forecasts in the majority of locations in north Australia, especially in the Top End region, Cape York, and central Queensland. At times of strong anomalies in the July-August SOI, the forecasts are reliably emphatic. Furthermore, predictions using tropical Pacific sea surface temperatures (SSTs) as the predictor are also tested. While short-lead (July-August predictor) forecasts are more skillful using the SOI, long-lead (May-June predictor) forecasts are more skillful using Pacific SSTs, indicative of the longer-term memory present in the ocean.

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.

Development of near infrared analysis of faeces to estimate non-grass proportions in diets selected by cattle grazing tropical pastures

Grass (monocots) and non-grass (dicots) proportions in ruminant diets are important nutritionally because the non-grasses are usually higher in nutritive value, particularly protein, than the grasses, especially in tropical pastures. For ruminants grazing tropical pastures where the grasses are C-4 species and most non-grasses are C-3 species, the ratio of C-13/C-12 in diet and faeces, measured as delta C-13 parts per thousand, is proportional to dietary non-grass%. This paper describes the development of a faecal near infrared (NIR) spectroscopy calibration equation for predicting faecal delta C-13 from which dietary grass and non-grass proportions can be calculated. Calibration development used cattle faeces derived from diets containing only C-3 non-grass and C-4 grass components, and a series of expansion and validation steps was employed to develop robustness and predictive reliability. The final calibration equation contained 1637 samples and faecal delta C-13 range (parts per thousand) of [12.27]-[27.65]. Calibration statistics were: standard error of calibration (SEC) of 0.78, standard error of cross-validation (SECV) of 0.80, standard deviation (SD) of reference values of 3.11 and R-2 of 0.94. Validation statistics for the final calibration equation applied to 60 samples were: standard error of prediction (SEP) of 0.87, bias of -0.15, R-2 of 0.92 and RPD of 3.16. The calibration equation was also tested on faeces from diets containing C-4 non-grass species or temperate C-3 grass species. Faecal delta C-13 predictions indicated that the spectral basis of the calibration was not related to C-13/C-12 ratios per se but to consistent differences between grasses and non-grasses in chemical composition and that the differences were modified by photosynthetic pathway. Thus, although the calibration equation could not be used to make valid faecal delta C-13 predictions when the diet contained either C-3 grass or C-4 non-grass, it could be used to make useful estimates of dietary non-grass proportions. It could also be ut :sed to make useful estimates of non-grass in mixed C-3 grass/non-grass diets by applying a modified formula to calculate non-grass from predicted faecal delta C-13. The development of a robust faecal-NIR calibration equation for estimating non-grass proportions in the diets of grazing cattle demonstrated a novel and useful application of NIR spectroscopy in agriculture.