Maximum likelihood estimation of natural mortality and quantification of temperature effects on catchability of brown tiger prawn (Penaeus esculentus) in Moreton Bay (Australia) using logbook data

It is common to model the dynamics of fisheries using natural and fishing mortality rates estimated independently using two separate analyses. Fishing mortality is routinely estimated from widely available logbook data, whereas natural mortality estimations have often required more specific, less frequently available, data. However, in the case of the fishery for brown tiger prawn (Penaeus esculentus) in Moreton Bay, both fishing and natural mortality rates have been estimated from logbook data. The present work extended the fishing mortality model to incorporate an eco-physiological response of tiger prawn to temperature, and allowed recruitment timing to vary from year to year. These ecological characteristics of the dynamics of this fishery were ignored in the separate model that estimated natural mortality. Therefore, we propose to estimate both natural and fishing mortality rates within a single model using a consistent set of hypotheses. This approach was applied to Moreton Bay brown tiger prawn data collected between 1990 and 2010. Natural mortality was estimated by maximum likelihood to be equal to 0.032 ± 0.002 week−1, approximately 30% lower than the fixed value used in previous models of this fishery (0.045 week−1).

Improved estimation of size-transition matrices using tag-recapture data

We derive a new method for determining size-transition matrices (STMs) that eliminates probabilities of negative growth and accounts for individual variability. STMs are an important part of size-structured models, which are used in the stock assessment of aquatic species. The elements of STMs represent the probability of growth from one size class to another, given a time step. The growth increment over this time step can be modelled with a variety of methods, but when a population construct is assumed for the underlying growth model, the resulting STM may contain entries that predict negative growth. To solve this problem, we use a maximum likelihood method that incorporates individual variability in the asymptotic length, relative age at tagging, and measurement error to obtain von Bertalanffy growth model parameter estimates. The statistical moments for the future length given an individual’s previous length measurement and time at liberty are then derived. We moment match the true conditional distributions with skewed-normal distributions and use these to accurately estimate the elements of the STMs. The method is investigated with simulated tag–recapture data and tag–recapture data gathered from the Australian eastern king prawn (Melicertus plebejus).

Effects of severity of threshing damage on seed quality of Gatton panic Panicum maximum)

Different degrees of severity of threshing were imposed during combine-harvesting of seed of Gatton panic, a cultivar of Panicum maximum , to determine effects of degree of threshing damage on subsequent properties of seed. Threshing cylinder peripheral speeds and concave clearances covering the normal range employed commercially were varied experimentally in the harvest of 2 crops grown in north Queensland. Harvested seed was dried and cleaned, then stored under ambient conditions. The extent of physical damage was measured, and samples were tested at intervals for viability, germination, dormancy and seedling emergence from soil in a glasshouse and in the field over the 2 seasons following harvest. Physical damage increased as peripheral rotor speed rose and (though less markedly) as concave clearance was reduced. As the level of damage increased, viability was progressively reduced, life expectancy was shortened, and dormancy was broken. When the consequences were measured as seedling emergence from soil, the adverse effects on viability tended to cancel out the benefits of dormancy-breaking, leaving few net differences attributable to the degree of threshing severity. We concluded that there would be no value in trying to manipulate the quality of seed produced for normal commercial use through choice of cylinder settings, but that deliberate light or heavy threshing could benefit special-purpose seed, destined, respectively, for long-term storage or immediate use.

Assessing the relationship between shire winter crop yield and seasonal variability of the MODIS NDVI and EVI images

Australian researchers have been developing robust yield estimation models, based mainly on the crop growth response to water availability during the crop season. However, knowledge of spatial distribution of yields within and across the production regions can be improved by the use of remote sensing techniques. Images of Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation indices, available since 1999, have the potential to contribute to crop yield estimation. The objective of this study was to analyse the relationship between winter crop yields and the spectral information available in MODIS vegetation index images at the shire level. The study was carried out in the Jondaryan and Pittsworth shires, Queensland , Australia . Five years (2000 to 2004) of 250m resolution, 16-day composite of MODIS Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) images were used during the winter crop season (April to November). Seasonal variability of the profiles of the vegetation index images for each crop season using different regions of interest (cropping mask) were displayed and analysed. Correlation analysis between wheat and barley yield data and MODIS image values were also conducted. The results showed high seasonal variability in the NDVI and EVI profiles, and the EVI values were consistently lower than those of the NDVI. The highest image values were observed in 2003 (in contrast to 2004), and were associated with rainfall amount and distribution. The seasonal variability of the profiles was similar in both shires, with minimum values in June and maximum values at the end of August. NDVI and EVI images showed sensitivity to seasonal variability of the vegetation and exhibited good association (e.g. r = 0.84, r = 0.77) with winter crop yields.

Ultraviolet radiation effects on fruit surface respiration and chlorophyll fluorescence

High-value fruit crops are exposed to a range of environmental conditions that can reduce fruit quality. Solar injury (SI) or sunburn is a common disorder in tropical, sub-tropical, and temperate climates and is related to: 1) high fruit surface temperature; 2) high visible light intensity; and, 3) ultraviolet radiation (UV). Positional changes in fruit that are caused by increased weight or abrupt changes that result from summer pruning, limb breakage, or other damage to the canopy can expose fruit to high solar radiation levels, increased fruit surface temperatures, and increased UV exposure that are higher than the conditions to which they are adapted. In our studies, we examined the effects of high fruit surface temperature, saturating photosynthetically-active radiation (PAR), and short-term UV exposure on chlorophyll fluorescence, respiration, and photosynthesis of fruit peel tissues from tropical and temperate fruit in a simulation of these acute environmental changes. All tropical fruits (citrus, macadamia, avocado, pineapple, and custard apple) and the apple cultivars 'Gala', 'Gold Rush', and 'Granny Smith' increased dark respiration (A0) when exposed to UV, suggesting that UV repair mechanisms were induced. The maximum quantum efficiency of photosystem II (Fv/Fm) and the quantum efficiency of photosystem II (ΦII) were unaffected, indicating no adverse effects on photosystem II (PSII). In contrast, 'Braeburn' apple had a reduced Fv/Fm with no increase in A0 on all sampling dates. There was a consistent pattern in all studies. When Fv/Fm was unaffected by UV treatment, A0 increased significantly. Conversely, when Fv/Fm was reduced by UV treatment, then A0 was unaffected. The pattern suggests that when UV repair mechanisms are effective, PSII is adequately protected, and that this protection occurs at the cost of higher respiration. However, when the UV repair mechanisms are ineffective, not only is PSII damaged, but there is additional short-term damage to the repair mechanisms, indicated by a lack of respiration to provide energy.

Soil sorption-desorption of phosphorus from piggery effluent compared with inorganic sources

The leaching of phosphorus (P) within soils can be a limiting consideration for the sustainable operation of intensive livestock enterprises. Sorption curves are widely used to assist estimation of P retention, though the effect of effluent constituents on their accuracy is not well understood. We conducted a series of P-sorption-desorption batch experiments with an Oxic Haplustalf (soil 1), Haplusterts (soils 2 and 3), and a Natrustalf (soil 4). Phosphorus sources included effluent, orthophosphate-P in a matrix replicating the effluent's salt constituents (the reference solution), and an orthophosphate-P solution. Treated soils were incubated for up to 193 days before sequential desorption extraction. Effluent constituents, probably the organic or particulate components, temporarily increased the vulnerability of sorbed-P to desorption. The increase in vulnerability was removed by 2-113 days of incubation (25 degrees C). Despite vigorous extraction for 20 consecutive days, some P sorbed as part of the treatments of soils 1 and 2 was not desorbed. The increased vulnerability due to effluent constituents lasted a maximum of about one cropping season and, for all other treatments, adsorption curves overestimated vulnerability to desorption. Therefore, adsorption curves provide a conservative estimate of vulnerability to desorption where effluent is used in continued crop production in these soils.

Estimation of light interception in research environments: A joint approach using directional light sensors and 3D virtual plants applied to sunflower (Helianthus annuus) and Arabidopsis thaliana in natural and artificial conditions

Light interception is a major factor influencing plant development and biomass production. Several methods have been proposed to determine this variable, but its calculation remains difficult in artificial environments with heterogeneous light. We propose a method that uses 3D virtual plant modelling and directional light characterisation to estimate light interception in highly heterogeneous light environments such as growth chambers and glasshouses. Intercepted light was estimated by coupling an architectural model and a light model for different genotypes of the rosette species Arabidopsis thaliana (L.) Heynh and a sunflower crop. The model was applied to plants of contrasting architectures, cultivated in isolation or in canopy, in natural or artificial environments, and under contrasting light conditions. The model gave satisfactory results when compared with observed data and enabled calculation of light interception in situations where direct measurements or classical methods were inefficient, such as young crops, isolated plants or artificial conditions. Furthermore, the model revealed that A. thaliana increased its light interception efficiency when shaded. To conclude, the method can be used to calculate intercepted light at organ, plant and plot levels, in natural and artificial environments, and should be useful in the investigation of genotype-environment interactions for plant architecture and light interception efficiency. This paper originates from a presentation at the 5th International Workshop on Functional–Structural Plant Models, Napier, New Zealand, November 2007.

Estimation in a multiplicative mixed model involving a genetic relationship matrix.

Genetic models partitioning additive and non-additive genetic effects for populations tested in replicated multi-environment trials (METs) in a plant breeding program have recently been presented in the literature. For these data, the variance model involves the direct product of a large numerator relationship matrix A, and a complex structure for the genotype by environment interaction effects, generally of a factor analytic (FA) form. With MET data, we expect a high correlation in genotype rankings between environments, leading to non-positive definite covariance matrices. Estimation methods for reduced rank models have been derived for the FA formulation with independent genotypes, and we employ these estimation methods for the more complex case involving the numerator relationship matrix. We examine the performance of differing genetic models for MET data with an embedded pedigree structure, and consider the magnitude of the non-additive variance. The capacity of existing software packages to fit these complex models is largely due to the use of the sparse matrix methodology and the average information algorithm. Here, we present an extension to the standard formulation necessary for estimation with a factor analytic structure across multiple environments.

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.

Influence of plant traits on production costs and profitability of strawberry in southeast Queensland.

In Queensland the subtropical strawberry ( Fragaria * ananassa) breeding program aims to combine traits into novel genotypes that increase production efficiency. The contribution of individual plant traits to cost and income under subtropical Queensland conditions was investigated, with the overall goal of improving the profitability of the industry through the release of new strawberry cultivars. The study involved specifying the production and marketing system using three cultivars of strawberry that are currently widely grown annually in southeast Queensland, developing methods to assess the economic impact of changes to the system, and identifying plant traits that influence outcomes from the system. From May through September P (price; $ punnet -1), V (monthly mass; tonne of fruit on the market) and M (calendar month; i.e. May=5) were found to be related ( r2=0.92) by the function (SE) P=4.741(0.469)-0.001630(0.0005) V-0.226(0.102) M using data from 2006 to 2010 for the Brisbane central market. Both income and cost elements in the gross margin were subject to sensitivity analysis. 'Harvesting' and 'Handling/Packing' 'Groups' of 'Activities' were the major contributors to variable costs (each >20%) in the gross margin analysis. Within the 'Harvesting Group', the 'Picking Activity' contributed most (>80%) with the trait 'display of fruit' having the greatest (33%) influence on the cost of the 'Picking Activity'. Within the 'Handling/Packing Group', the 'Packing Activity' contributed 50% of costs with the traits 'fruit shape', 'fruit size variation' and 'resistance to bruising' having the greatest (12-62%) influence on the cost of the 'Packing Activity'. Non-plant items (e.g. carton purchases) made up the other 50% of the costs within the 'Handling/Packing Group'. When any of the individual traits in the 'Harvesting' and 'Handling/Packing' groups were changed by one unit (on a 1-9 scale) the gross margin changed by up to 1%. Increasing yield increased the gross margin to a maximum (15% above present) at 1320 g plant -1 (94% above present). A 10% redistribution of total yield from September to May increased the gross margin by 23%. Increasing fruit size increased gross margin: a 75% increase in fruit size (to ~30 g) produced a 22% increase in the gross margin. The modified gross margin analysis developed in this study allowed simultaneous estimation of the gross margin for the producer and gross value of the industry. These parameters sometimes move in opposite directions.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Scale- and Context-Dependent Selection of Recreational Harvest Estimation Methods: The Australasian Experience

Fisheries managers are becoming increasingly aware of the need to quantify all forms of harvest, including that by recreational fishers. This need has been driven by both a growing recognition of the potential impact that noncommercial fishers can have on exploited resources and the requirement to allocate catch limits between different sectors of the wider fishing community in many jurisdictions. Marine recreational fishers are rarely required to report any of their activity, and some form of survey technique is usually required to estimate levels of recreational catch and effort. In this review, we describe and discuss studies that have attempted to estimate the nature and extent of recreational harvests of marine fishes in New Zealand and Australia over the past 20 years. We compare studies by method to show how circumstances dictate their application and to highlight recent developments that other researchers may find of use. Although there has been some convergence of approach, we suggest that context is an important consideration, and many of the techniques discussed here have been adapted to suit local conditions and to address recognized sources of bias. Much of this experience, along with novel improvements to existing approaches, have been reported only in "gray" literature because of an emphasis on providing estimates for immediate management purposes. This paper brings much of that work together for the first time, and we discuss how others might benefit from our experience.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.

A novel method for the age estimation of Saddletail snapper (Lutjanus malabaricus) using Fourier Transform-near infrared (FT-NIR) spectroscopy

Near infrared (NIR) spectroscopy was investigated as a potential rapid method of estimating fish age from whole otoliths of Saddletail snapper (Lutjanus malabaricus). Whole otoliths from 209 Saddletail snapper were extracted and the NIR spectral characteristics were acquired over a spectral range of 800–2780 nm. Partial least-squares models (PLS) were developed from the diffuse reflectance spectra and reference-validated age estimates (based on traditional sectioned otolith increments) to predict age for independent otolith samples. Predictive models developed for a specific season and geographical location performed poorly against a different season and geographical location. However, overall PLS regression statistics for predicting a combined population incorporating both geographic location and season variables were: coefficient of determination (R2) = 0.94, root mean square error of prediction (RMSEP) = 1.54 for age estimation, indicating that Saddletail age could be predicted within 1.5 increment counts. This level of accuracy suggests the method warrants further development for Saddletail snapper and may have potential for other fish species. A rapid method of fish age estimation could have the potential to reduce greatly both costs of time and materials in the assessment and management of commercial fisheries.

Phenotyping novel stay-green traits to capture genetic variation in senescence dynamics

Stay-green plants retain green leaves longer after anthesis and can have improved yield, particularly under water limitation. As senescence is a dynamic process, genotypes with different senescence patterns may exhibit similar final normalised difference vegetative index (NDVI). By monitoring NDVI from as early as awn emergence to maturity, we demonstrate that analysing senescence dynamics improves insight into genotypic stay-green variation. A senescence evaluation tool was developed to fit a logistic function to NDVI data and used to analyse data from three environments for a wheat (Triticum aestivum L.) population whose lines contrast for stay-green. Key stay-green traits were estimated including, maximum NDVI, senescence rate and a trait integrating NDVI variation after anthesis, as well as the timing from anthesis to onset, midpoint and conclusion of senescence. The integrative trait and the timing to onset and mid-senescence exhibited high positive correlations with yield and a high heritability in the three studied environments. Senescence rate was correlated with yield in some environments, whereas maximum NDVI was associated with yield in a drought-stressed environment. Where resources preclude frequent measurements, we found that NDVI measurements may be restricted to the period of rapid senescence, but caution is required when dealing with lines of different phenology. In contrast, regular monitoring during the whole period after flowering allows the estimation of senescence dynamics traits that may be reliably compared across genotypes and environments. We anticipate that selection for stay-green traits will enhance genetic progress towards high-yielding, stay-green germplasm.