Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected. (C) 2013 Elsevier B.V. All rights reserved.

Mixture of time-dependent growth models with an application to blue swimmer crab length-frequency data

Understanding how aquatic species grow is fundamental in fisheries because stock assessment often relies on growth dependent statistical models. Length-frequency-based methods become important when more applicable data for growth model estimation are either not available or very expensive. In this article, we develop a new framework for growth estimation from length-frequency data using a generalized von Bertalanffy growth model (VBGM) framework that allows for time-dependent covariates to be incorporated. A finite mixture of normal distributions is used to model the length-frequency cohorts of each month with the means constrained to follow a VBGM. The variances of the finite mixture components are constrained to be a function of mean length, reducing the number of parameters and allowing for an estimate of the variance at any length. To optimize the likelihood, we use a minorization–maximization (MM) algorithm with a Nelder–Mead sub-step. This work was motivated by the decline in catches of the blue swimmer crab (BSC) (Portunus armatus) off the east coast of Queensland, Australia. We test the method with a simulation study and then apply it to the BSC fishery data.

Modeling chickpea growth and development: Nitrogen accumulation and use

Quantitative information regarding nitrogen (N) accumulation and its distribution to leaves, stems and grains under varying environmental and growth conditions are limited for chickpea (Cicer arietinum L.). The information is required for the development of crop growth models and also for assessment of the contribution of chickpea to N balances in cropping systems. Accordingly, these processes were quantified in chickpea under different environmental and growth conditions (still without water or N deficit) using four field experiments and 1325 N measurements. N concentration ([N]) in green leaves was 50 mg g-1 up to beginning of seed growth, and then it declined linearly to 30 mg g-1 at the end of seed growth phase. [N] in senesced leaves was 12 mg g-1. Stem [N] decreased from 30 mg g-1 early in the season to 8 mg g-1 in senesced stems at maturity. Pod [N] was constant (35 mg g-1), but grain [N] decreased from 60 mg g-1 early in seed growth to 43 mg g-1 at maturity. Total N accumulation ranged between 9 and 30 g m-2. N accumulation was closely linked to biomass accumulation until maturity. N accumulation efficiency (N accumulation relative to biomass accumulation) was 0.033 g g-1 where total biomass was -2 and during early growth period, but it decreased to 0.0176 g g-1 during the later growth period when total biomass was >218 g m-2. During vegetative growth (up to first-pod), 58% of N was partitioned to leaves and 42% to stems. Depending on growth conditions, 37-72% of leaf N and 12-56% of stem N was remobilized to the grains. The parameter estimates and functions obtained in this study can be used in chickpea simulation models to simulate N accumulation and distribution.

Stochastic growth of the Eastern King Prawn (melicertus plebejus (hess, 1865)) harvested off Eastern Australia

Stochastic growth models were fitted to length-increment data of eastern king prawns, Melicertus plebejus (Hess, 1865), tagged across eastern Australia. The estimated growth parameters and growth transition matrix are for each sex representative of the species' geographical distribution. Our study explicitly displays the stochastic nature of prawn growth. Capturing length-increment growth heterogeneity for short-lived exploited species such as prawns that cannot be readily aged is essential for length-based modelling and improved management.

Short-term responses of leaf growth rate to water deficit scale up to whole-plant and crop levels: An integrated modelling approach in maize.

Physiological and genetic studies of leaf growth often focus on short-term responses, leaving a gap to whole-plant models that predict biomass accumulation, transpiration and yield at crop scale. To bridge this gap, we developed a model that combines an existing model of leaf 6 expansion in response to short-term environmental variations with a model coordinating the development of all leaves of a plant. The latter was based on: (1) rates of leaf initiation, appearance and end of elongation measured in field experiments; and (2) the hypothesis of an independence of the growth between leaves. The resulting whole-plant leaf model was integrated into the generic crop model APSIM which provided dynamic feedback of environmental conditions to the leaf model and allowed simulation of crop growth at canopy level. The model was tested in 12 field situations with contrasting temperature, evaporative demand and soil water status. In observed and simulated data, high evaporative demand reduced leaf area at the whole-plant level, and short water deficits affected only leaves developing during the stress, either visible or still hidden in the whorl. The model adequately simulated whole-plant profiles of leaf area with a single set of parameters that applied to the same hybrid in all experiments. It was also suitable to predict biomass accumulation and yield of a similar hybrid grown in different conditions. This model extends to field conditions existing knowledge of the environmental controls of leaf elongation, and can be used to simulate how their genetic controls flow through to yield.

Simulating pasture growth rates in Australian and New Zealand grazing systems

DairyMod, EcoMod, and the SGS Pasture Model are mechanistic biophysical models developed to explore scenarios in grazing systems. The aim of this manuscript was to test the ability of the models to simulate net herbage accumulation rates of ryegrass-based pastures across a range of environments and pasture management systems in Australia and New Zealand. Measured monthly net herbage accumulation rate and accumulated yield data were collated from ten grazing system experiments at eight sites ranging from cool temperate to subtropical environments. The local climate, soil, pasture species, and management (N fertiliser, irrigation, and grazing or cutting pattern) were described in the model for each site, and net herbage accumulation rates modelled. The model adequately simulated the monthly net herbage accumulation rates across the range of environments, based on the summary statistics and observed patterns of seasonal growth, particularly when the variability in measured herbage accumulation rates was taken into account. Agreement between modelled and observed growth rates was more accurate and precise in temperate than in subtropical environments, and in winter and summer than in autumn and spring. Similarly, agreement between predicted and observed accumulated yields was more accurate than monthly net herbage accumulation. Different temperature parameters were used to describe the growth of perennial ryegrass cultivars and annual ryegrass; these differences were in line with observed growth patterns and breeding objectives. Results are discussed in the context of the difficulties in measuring pasture growth rates and model limitations.

Models for the field-based toxicity of copper and zinc salts to wheat in 11 Australian soils and comparison to laboratory-based models

Laboratory-based relationships that model the phytotoxicity of metals using soil properties have been developed. This paper presents the first field-based phytotoxicity relationships. Wheat(Triticum aestivum L) was grown at 11 Australian field sites at which soil was spiked with copper (Cu) and zinc (Zn) salts. Toxicity was measured as inhibition of plant growth at 8 weeks and grain yield at harvest. The added Cu and Zn EC10 values for both endpoints ranged from approximately 3 to 4760 mg/kg. There were no relationships between field-based 8-week biomass and grain yield toxicity values for either metal. Cu toxicity was best modelled using pH and organic carbon content while Zn toxicity was best modelled using pH and the cation exchange capacity. The best relationships estimated toxicity within a factor of two of measured values. Laboratory-based phytotoxicity relationships could not accurately predict field-based phytotoxicity responses.

Genetic variation and responses in reproductive performance of sows in lines selected for growth rate under restricted feeding

The objective of this study was to examine genetic changes in reproduction traits in sows (total number born (TNB), number born alive (NBA), average piglet birth weight (ABW) and number of piglets weaned (NW), body weight prior to mating (MW), gestation length (GL) and daily food intake during lactation (DFI)) in lines of Large White pigs divergently selected over 4 years for high and low post-weaning growth rate on a restricted ration. Heritabilities and repeatabilities of the reproduction traits were also determined. The analyses were carried out on 913 litter records using average information-restricted maximum likelihood method applied to single trait animal models. Estimates of heritability for most traits were small, except for ABW (0·33) and MW (0·35). Estimates of repeatability were slightly higher than those of heritability for TNB, NBA and NW, but they were almost identical for ABW, MW, GL and DFI. After 4 years of selection, the high growth line sows had significantly heavier body weight prior to mating and produced significantly more piglets born alive with heavier average birth weight than the low line sows. There were, however, no statistical differences between the selected lines in TNB or NW. The lower food intake of high relative to low line sows during lactation was not significant, indicating that daily food intake differences found between grower pigs in the high and low lines (2·71 v. 2·76 kg/day, s.e.d. 0·024) on ad libitum feeding were not fully expressed in lactating sows. It is concluded that selection for growth rate on the restricted ration resulted in beneficial effects on important measures of reproductive performance of the sows.

Parameter estimation and sensitivity analysis of fat deposition models in beef steers using acslXtreme.

The Davis Growth Model (a dynamic steer growth model encompassing 4 fat deposition models) is currently being used by the phenotypic prediction program of the Cooperative Research Centre (CRC) for Beef Genetic Technologies to predict P8 fat (mm) in beef cattle to assist beef producers meet market specifications. The concepts of cellular hyperplasia and hypertrophy are integral components of the Davis Growth Model. The net synthesis of total body fat (kg) is calculated from the net energy available after accounting tor energy needs for maintenance and protein synthesis. Total body fat (kg) is then partitioned into 4 fat depots (intermuscular, intramuscular, subcutaneous, and visceral). This paper reports on the parameter estimation and sensitivity analysis of the DNA (deoxyribonucleic acid) logistic growth equations and the fat deposition first-order differential equations in the Davis Growth Model using acslXtreme (Hunstville, AL, USA, Xcellon). The DNA and fat deposition parameter coefficients were found to be important determinants of model function; the DNA parameter coefficients with days on feed >100 days and the fat deposition parameter coefficients for all days on feed. The generalized NL2SOL optimization algorithm had the fastest processing time and the minimum number of objective function evaluations when estimating the 4 fat deposition parameter coefficients with 2 observed values (initial and final fat). The subcutaneous fat parameter coefficient did indicate a metabolic difference for frame sizes. The results look promising and the prototype Davis Growth Model has the potential to assist the beef industry meet market specifications.

A novel breeding programme for improved growth in barramundi Lates calcarifer (Bloch) using foundation stock from progeny-tested parents

Rapid genetic gains for growth in barramundi ( Lates calcarifer) appear achievable by starting a breeding programme using foundation stock from progeny tested broodstock. The potential gains of this novel breeding design were investigated using biologically feasible scenarios tested with computer simulation models. The design involves the production of a large number of full-sib families using artificial mating which are compared in common growout conditions. The estimated breeding values of their paternal parents are calculated using a binomial probit analysis to assess their suitability as foundation broodstock. The programme can theoretically yield faster rates of genetic gain compared to other breeding programmes for aquaculture species. Assuming a heritability of 0.25 for growth, foundation broodstock evaluated in two years had breeding values for faster growth ranging from 21% to 51% depending on the genetic diversity of stock under evaluation. As a comparison it will take between nine and twenty-two years to identify broodstock with similar breeding values in a contemporary barramundi breeding programme.

Modeling population growth and site specific control of the invasive Lantana camara L. (Verbenaceae) under differing fire regimes

It is at the population level that an invasion either fails or succeeds. Lantana camara L. (Verbenaceae) is a weed of great significance in Queensland Australia and globally but its whole life-history ecology is poorly known. Here we used 3 years of field data across four land use types (farm, hoop pine plantation and two open eucalyptus forests, including one with a triennial fire regime) to parameterise the weed’s vital rates and develop size-structured matrix models. Lantana camara in its re-colonization phase, as observed in the recently cleared hoop pine plantation, was projected to increase more rapidly (annual growth rate, λ = 3.80) than at the other three sites (λ 1.88–2.71). Elasticity analyses indicated that growth contributed more (64.6 %) to λ than fecundity (18.5 %) or survival (15.5 %), while across size groups, the contribution was of the order: juvenile (19–27 %) ≥ seed (17–28 %) ≥ seedling (16–25 %) > small adult (4–26 %) ≥ medium adult (7–20 %) > large adult (0–20 %). From a control perspective it is difficult to determine a single weak point in the life cycle of lantana that might be exploited to reduce growth below a sustaining rate. The triennial fire regime applied did not alter the population elasticity structure nor resulted in local control of the weed. However, simulations showed that, except for the farm population, periodic burning could work within 4–10 years for control of the weed, but fire frequency should increase to at least once every 2 years. For the farm, site-specific control may be achieved by 15 years if the biennial fire frequency is tempered with increased burning intensity.

Long-term variation of tropical rock lobster panulirus ornatus (decapoda, palinuridae) growth inTorres Strait, Australia

This analysis of all carapace length measurements collected between 1989 and 2009, during scientific surveys, describes the variation of tropical rock lobster, Panulirus ornatus, somatic growth in Torres Strait. Multiple models of carapace length frequency distributions were compared by maximum likelihood to determine which hypotheses were most supported by the data. The best model assumed sex and cohort-specific Von Bertalanffy's parameters. These estimates are consistent with results derived from tagging data collected in the 1980s and provide new information on parameters' uncertainty. In the past two decades, growth rates have fluctuated inter-annually without displaying any distinctive trend. Associated uncertainties are large, suggesting that sampling will need to be intensified in order to detect an effect of climate change.

Estimates of genetic parameters for growth and wood properties in Eucalyptus pellita F. Muell. to support tree breeding in Vietnam

Key message Eucalyptus pellita demonstrated good growth and wood quality traits in this study, with young plantation grown timber being suitable for both solid and pulp wood products. All traits examined were under moderate levels of genetic control with little genotype by environment interaction when grown on two contrasting sites in Vietnam. Context Eucalyptus pellita currently has a significant role in reforestation in the tropics. Research to support expanded of use of this species is needed: particularly, research to better understand the genetic control of key traits will facilitate the development of genetically improved planting stock. Aims This study aimed to provide estimates of the heritability of diameter at breast height over bark, wood basic density, Kraft pulp yield, modulus of elasticity and microfibril angle, and the genetic correlations among these traits, and understand the importance of genotype by environment interactions in Vietnam. Methods Data for diameter and wood properties were collected from two 10-year-old, open-pollinated progeny trials of E. pellita in Vietnam that evaluated 104 families from six native range and three orchard sources. Wood properties were estimated from wood samples using near-infrared (NIR) spectroscopy. Data were analysed using mixed linear models to estimate genetic parameters (heritability, proportion of variance between seed sources and genetic correlations). Results Variation among the nine sources was small compared to additive variance. Narrow-sense heritability and genetic correlation estimates indicated that simultaneous improvements in most traits could be achieved from selection among and within families as the genetic correlations among traits were either favourable or close to zero. Type B genetic correlations approached one for all traits suggesting that genotype by environment interactions were of little importance. These results support a breeding strategy utilizing a single breeding population advanced by selecting the best individuals across all seed sources. Conclusion Both growth and wood properties have been evaluated. Multi-trait selection for growth and wood property traits will lead to more productive populations of E. pellita both with improved productivity and improved timber and pulp properties.

Integrating pest population models with biophysical crop models to better represent the farming system

Farming systems frameworks such as the Agricultural Production Systems simulator (APSIM) represent fluxes through the soil, plant and atmosphere of the system well, but do not generally consider the biotic constraints that function within the system. We designed a method that allowed population models built in DYMEX to interact with APSIM. The simulator engine component of the DYMEX population-modelling platform was wrapped within an APSIM module allowing it to get and set variable values in other APSIM models running in the simulation. A rust model developed in DYMEX is used to demonstrate how the developing rust population reduces the crop's green leaf area. The success of the linking process is seen in the interaction of the two models and how changes in rust population on the crop's leaves feedback to the APSIM crop modifying the growth and development of the crop's leaf area. This linking of population models to simulate pest populations and biophysical models to simulate crop growth and development increases the complexity of the simulation, but provides a tool to investigate biotic constraints within farming systems and further moves APSIM towards being an agro-ecological framework.

Assessment of functional forms of crop yield loss models of invasive plant species applied in decision support tools and bioeconomic modelling

AbstractObjectives Decision support tools (DSTs) for invasive species management have had limited success in producing convincing results and meeting users' expectations. The problems could be linked to the functional form of model which represents the dynamic relationship between the invasive species and crop yield loss in the DSTs. The objectives of this study were: a) to compile and review the models tested on field experiments and applied to DSTs; and b) to do an empirical evaluation of some popular models and alternatives. Design and methods This study surveyed the literature and documented strengths and weaknesses of the functional forms of yield loss models. Some widely used models (linear, relative yield and hyperbolic models) and two potentially useful models (the double-scaled and density-scaled models) were evaluated for a wide range of weed densities, maximum potential yield loss and maximum yield loss per weed. Results Popular functional forms include hyperbolic, sigmoid, linear, quadratic and inverse models. Many basic models were modified to account for the effect of important factors (weather, tillage and growth stage of crop at weed emergence) influencing weed–crop interaction and to improve prediction accuracy. This limited their applicability for use in DSTs as they became less generalized in nature and often were applicable to a much narrower range of conditions than would be encountered in the use of DSTs. These factors' effects could be better accounted by using other techniques. Among the model empirically assessed, the linear model is a very simple model which appears to work well at sparse weed densities, but it produces unrealistic behaviour at high densities. The relative-yield model exhibits expected behaviour at high densities and high levels of maximum yield loss per weed but probably underestimates yield loss at low to intermediate densities. The hyperbolic model demonstrated reasonable behaviour at lower weed densities, but produced biologically unreasonable behaviour at low rates of loss per weed and high yield loss at the maximum weed density. The density-scaled model is not sensitive to the yield loss at maximum weed density in terms of the number of weeds that will produce a certain proportion of that maximum yield loss. The double-scaled model appeared to produce more robust estimates of the impact of weeds under a wide range of conditions. Conclusions Previously tested functional forms exhibit problems for use in DSTs for crop yield loss modelling. Of the models evaluated, the double-scaled model exhibits desirable qualitative behaviour under most circumstances.