Hazard Function Models to Estimate Mortality Rates Affecting Fish Populations with Application to the Sea Mullet (Mugil cephalus) Fishery on the Queensland Coast (Australia)

Fisheries management agencies around the world collect age data for the purpose of assessing the status of natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided unbiased estimations of relevant quantities. The method was applied to the data from the Queensland (Australia) sea mullet (Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate of natural mortality affecting this stock: 0.22±0.08 year −1 .

Are high-impact species predictable? An analysis of naturalised grasses in northern Australia

Predicting which species are likely to cause serious impacts in the future is crucial for targeting management efforts, but the characteristics of such species remain largely unconfirmed. We use data and expert opinion on tropical and subtropical grasses naturalised in Australia since European settlement to identify naturalised and high-impact species and subsequently to test whether high-impact species are predictable. High-impact species for the three main affected sectors (environment, pastoral and agriculture) were determined by assessing evidence against pre-defined criteria. Twenty-one of the 155 naturalised species (14%) were classified as high-impact, including four that affected more than one sector. High-impact species were more likely to have faster spread rates (regions invaded per decade) and to be semi-aquatic. Spread rate was best explained by whether species had been actively spread (as pasture), and time since naturalisation, but may not be explanatory as it was tightly correlated with range size and incidence rate. Giving more weight to minimising the chance of overlooking high-impact species, a priority for biosecurity, meant a wider range of predictors was required to identify high-impact species, and the predictive power of the models was reduced. By-sector analysis of predictors of high impact species was limited by their relative rarity, but showed sector differences, including to the universal predictors (spread rate and habitat) and life history. Furthermore, species causing high impact to agriculture have changed in the past 10 years with changes in farming practice, highlighting the importance of context in determining impact. A rationale for invasion ecology is to improve the prediction and response to future threats. Although our study identifies some universal predictors, it suggests improved prediction will require a far greater emphasis on impact rather than invasiveness, and will need to account for the individual circumstances of affected sectors and the relative rarity of high-impact species.

Predicting survival of western rock lobsters Panulirus cygnus using discriminant analysis of hemolymph parameters taken immediately following simulated handling treatments

Instances of morbidity amongst rock lobsters (Panulirus cygnus) arriving at factories in Western Australia (WA) have been attributed to stress during post-harvest handling. This study used discriminant analysis to determine whether physiological correlates of stress following a period of simulated post-harvest handling had any validity as predictors of future rejection or morbidity of western rock lobsters. Groups of 230 western rock lobsters were stored for 6 h in five environments (submerged/flowing sea water, submerged/re-circulating sea water, humid air, flowing sea water spray, and re-circulated sea water spray). The experiment was conducted in late spring (ambient sea water 22°C), and repeated again in early autumn (ambient sea water 26°C). After 6 h treatment, each lobster was graded for acceptability for live export, numbered, and its hemolymph was sampled. The samples were analysed for a number of physiological and health status parameters. The lobsters were then stored for a week in tanks in the live lobster factory to record mortality. The mortality of lobsters in the factory was associated with earlier deviations in hemolymph parameters as they emerged from the storage treatments. Discriminant analysis (DA) of the hemolymph assays enabled the fate of 80-90% of the lobsters to be correctly categorised within each experiment. However, functions derived from one experiment were less accurate at predicting mortality when applied to the other experiments. One of the reasons for this was the higher mortality and the more severe patho-physiological changes observed in lobsters stored in humid air or sprays at the higher temperature. The analysis identified lactate accumulation during emersion and associated physiological and hemocyte-related effects as a major correlate of mortality. Reducing these deviations, for example by submerged transport, is expected to ensure high levels of survival. None of the indicators tested predicted mortality with total accuracy. The simplest and most accurate means of comparing emersed treatments was to count the mortality afterwards.

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.