A spatially explicit multi-competitor coexistence model of penaeid (shrimp) distribution on the Australian Great Barrier Reef

A spatially explicit multi-competitor coexistence model was developed for meta-populations of prawns (shrimp) occupying habitat patches across the Great Barrier Reef, where dispersal was localised and dispersal rates varied between species. Prawns were modelled as individuals moving to and from patches or cells according to pre-set decision rules. The landscape was simulated as a matrix of cells with each cell having a spatially explicit survival index for each species. Mixed species prawn assemblages moved over this simplified spatially explicit landscape. A low level of chronic random environmental disturbance was assumed (cyclone and tropical storm damage) with additional acute spatially confined disturbance due to commercial trawling, modelled as an increase in mortality affecting inter-specific competition. The general form of the results was for increased disturbance to favour good-colonising "generalist" species at the expense of good-competitor "specialists". Increasing fishing mortality (local patch extinctions) combined with poor colonising ability resulted in low equilibrium abundance for even the best competitor, while in the same circumstances the poorest competitor but best coloniser could have the highest equilibrium abundance. This mimics the switch from high-value prawn species to lower-value prawn species as trawl effort increases, reflected in historic catch and effort logbook data and reported anecdotaly from the north Queensland trawl fleet. To match the observed distribution and behaviour of prawn assemblages, a combination inter-species competition, a spatially explicit landscape, and a defined pattern of disturbance (trawling) was required. Modelling this combination could simulate not only general trends in spatial distribution of each of prawn species but also localised concentrations observed in the survey data

Simulating crop-parasitic weed interactions using APSIM: Model evaluation and application

The parasitic weed Orobanche crenata inflicts major damage on faba bean, lentil, pea and other crops in Mediterranean environments. The development of methods to control O. crenata is to a large extent hampered by the complexity of host-parasite systems. Using a model of host-parasite interactions can help to explain and understand this intricacy. This paper reports on the evaluation and application of a model simulating host-parasite competition as affected by environment and management that was implemented in the framework of the Agricultural Production Systems Simulator (APSIM). Model-predicted faba bean and O. crenata growth and development were evaluated against independent data. The APSIM-Fababean and -Parasite modules displayed a good capability to reproduce effects of pedoclimatic conditions, faba bean sowing date and O. crenata infestation on host-parasite competition. The r(2) values throughout exceeded 0.84 (RMSD: 5.36 days) for phenological, 0.85 (RMSD: 223.00 g m(-2)) for host growth and 0.78 (RMSD: 99.82 g m(-2)) for parasite growth parameters. Inaccuracies of simulated faba bean root growth that caused some bias of predicted parasite number and host yield loss may be dealt with by more flexibly simulating vertical root distribution. The model was applied in simulation experiments to determine optimum sowing windows for infected and non-infected faba bean in Mediterranean environments. Simulation results proved realistic and testified to the capability of APSIM to contribute to the development of tactical approaches in parasitic weed control.

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.

A hidden Markov model approach for determining vessel activity from vessel monitoring system data.

Many fisheries worldwide have adopted vessel monitoring systems (VMS) for compliance purposes. An added benefit of these systems is that they collect a large amount of data on vessel locations at very fine spatial and temporal scales. This data can provide a wealth of information for stock assessment, research, and management. However, since most VMS implementations record vessel location at set time intervals with no regard to vessel activity, some methodology is required to determine which data records correspond to fishing activity. This paper describes a probabilistic approach, based on hidden Markov models (HMMs), to determine vessel activity. A HMM provides a natural framework for the problem and, by definition, models the intrinsic temporal correlation of the data. The paper describes the general approach that was developed and presents an example of this approach applied to the Queensland trawl fishery off the coast of eastern Australia. Finally, a simulation experiment is presented that compares the misallocation rates of the HMM approach with other approaches.

A Bayesian hurdle model for analysis of an insect resistance monitoring database

Motivated by the analysis of the Australian Grain Insect Resistance Database (AGIRD), we develop a Bayesian hurdle modelling approach to assess trends in strong resistance of stored grain insects to phosphine over time. The binary response variable from AGIRD indicating presence or absence of strong resistance is characterized by a majority of absence observations and the hurdle model is a two step approach that is useful when analyzing such a binary response dataset. The proposed hurdle model utilizes Bayesian classification trees to firstly identify covariates and covariate levels pertaining to possible presence or absence of strong resistance. Secondly, generalized additive models (GAMs) with spike and slab priors for variable selection are fitted to the subset of the dataset identified from the Bayesian classification tree indicating possibility of presence of strong resistance. From the GAM we assess trends, biosecurity issues and site specific variables influencing the presence of strong resistance using a variable selection approach. The proposed Bayesian hurdle model is compared to its frequentist counterpart, and also to a naive Bayesian approach which fits a GAM to the entire dataset. The Bayesian hurdle model has the benefit of providing a set of good trees for use in the first step and appears to provide enough flexibility to represent the influence of variables on strong resistance compared to the frequentist model, but also captures the subtle changes in the trend that are missed by the frequentist and naive Bayesian models. © 2014 Springer Science+Business Media New York.

The Ethics of a Co-regulatory Model for Farm Animal Welfare Research

Standards for farm animal welfare are variously managed at a national level by government-led regulatory control, by consumer-led welfare economics and co-regulated control in a partnership between industry and government. In the latter case the control of research to support animal welfare standards by the relevant industry body may lead to a conflict of interest on the part of researchers, who are dependent on industry for continued research funding. We examine this dilemma by reviewing two case studies of research published under an Australian co-regulated control system. Evidence of unsupported conclusions that are favourable to industry is provided, suggesting that researchers do experience a conflict of interest that may influence the integrity of the research. Alternative models for the management of research are discussed, including the establishment of an independent research management body for animal welfare because of its public good status and the use of public money derived from taxation, with representation from government, industry, consumers, and advocacy groups.

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.