Stout Whiting Fishery Summary : Queensland Total Allowable Catch for 2016

In Queensland, stout whiting are fished by Danish seine and fish otter-trawl methods between Sandy Cape and the Queensland-New South Wales border. The fishery is currently identified by a T4 symbol and is operated by two primary quota holders. Since 1997, T4 management has been informed by annual stock assessments in order to determine a total allowable commercial catch (TACC) quota. The TACC is assessed before the start of each fishing year using statistical methodologies. This includes evaluation of trends in fish catch-rates and catch-at-age frequencies against management reference points. The T4 stout whiting TACC for 2014 was adjusted down to 1150 t as a result of elevated estimates of fishing mortality and remained unchanged in 2015 (2013 TACC = 1350 t quota). Two T4 vessels fished for stout whiting in the 2015 fishing year, harvesting 663 t from Queensland waters. Annual T4 landings of stout whiting averaged about 713 t for the fishing years 2013–2015, with a maximum harvest in the last 10 fishing years of 1140 t and a maximum historical harvest of 2400 t in the 1995. Stout whiting catch rates from both Queensland and New South Wales were analysed for all vessels, areas and fishing gears. The 2015 catch rate index was equal to 0.85, down 15% compared to the 2010–2015 fishing year average (reference point =1). The stout whiting fish length and otolith weight frequencies indicated larger and older fish in the calendar year 2014. This data was translated to show improved measures of fish survival at about 38% per year and near the reference point of about 41%. Together, the stout whiting catch rate and survival indicators show the fishery was sustainable. Earlier population modelling conducted for the year 2013 also suggested the stock was sustainable, but the estimate was only marginally above the biomass for maximum sustainable yield. Irrespective, reasons for reduced catch rates should be examined further and interpreted with precaution, particularly given the TACC has been under-caught in many years. For setting of the 2016 TACC, alternate analyses and reference points were compared to address data uncertainties and provide options for quota change. The results were dependent on the stock indicator and harvest procedure used. Uncertainty in all TACC estimates should be considered as they were sensitive to the data inputs and assumptions. For the 2016 T4 fishing year, upper levels of harvest should be limited to 1000–1100 t following procedure equation 1, with target levels of harvest at 750–850 t for procedure equation 2. Use of these estimates to set TACC will depend on management and industry intentions.

Statistical forecasting of the Australian macadamia crop

Three types of forecasts of the total Australian production of macadamia nuts (t nut-in-shell) have been produced early each year since 2001. The first is a long-term forecast, based on the expected production from the tree census data held by the Australian Macadamia Society, suitably scaled up for missing data and assumed new plantings each year. These long-term forecasts range out to 10 years in the future, and form a basis for industry and market planning. Secondly, a statistical adjustment (termed the climate-adjusted forecast) is made annually for the coming crop. As the name suggests, climatic influences are the dominant factors in this adjustment process, however, other terms such as bienniality of bearing, prices and orchard aging are also incorporated. Thirdly, industry personnel are surveyed early each year, with their estimates integrated into a growers and pest-scouts forecast. Initially conducted on a 'whole-country' basis, these models are now constructed separately for the six main production regions of Australia, with these being combined for national totals. Ensembles or suites of step-forward regression models using biologically-relevant variables have been the major statistical method adopted, however, developing methodologies such as nearest-neighbour techniques, general additive models and random forests are continually being evaluated in parallel. The overall error rates average 14% for the climate forecasts, and 12% for the growers' forecasts. These compare with 7.8% for USDA almond forecasts (based on extensive early-crop sampling) and 6.8% for coconut forecasts in Sri Lanka. However, our somewhatdisappointing results were mainly due to a series of poor crops attributed to human reasons, which have now been factored into the models. Notably, the 2012 and 2013 forecasts averaged 7.8 and 4.9% errors, respectively. Future models should also show continuing improvement, as more data-years become available.

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%.

Towards responsible native fish stocking: identifying management concerns and appropriate research methodologies

A major outcome of this project has been the identification and prioritisation of the major management issues related to the ecological impacts of fish stocking and the elucidation of appropriate research methodologies that can be used to investigate these issues. This information is paramount to development of the relevant research projects that will lead to stocking activities aligned with world’s best practice, a requisite for ecologically sustainable recreational freshwater fisheries. In order to quantify the major management issues allied to the sustainability of freshwater fish stocking, stakeholders from around Australia were identified and sent a questionnaire to determine which particular issues they regarded as important. These stakeholders included fisheries managers or researchers from Federal, Territory and State jurisdictions although others, including representatives from environment and conservation agencies and peak recreational fishing and stocking groups were also invited to give their opinions. The survey was completed in late 2007 and the results analysed to give a prioritized list of key management issues relating to the impacts of native fish stocking activities. In the analysis, issues which received high priority rankings were flagged as potential topics for discussion at a future expert workshop. Identified high priority issues fell into the following core areas: marking techniques, genetics, population dynamics, introduction of pathogens and exotic biological material and ecological, biological and conservation issues. The next planned outcome, determination of the most appropriate methodologies to address these core issues in research projects, was addressed through the outputs of an expert workshop held in early 2008. Participants at this workshop agreed on a range of methodologies for addressing priority sustainability issues and decided under what circumstances that these methodologies should be employed.

Projecting future expansion of invasive species: Comparing and improving methodologies for species distribution modeling

Modeling the distributions of species, especially of invasive species in non-native ranges, involves multiple challenges. Here, we developed some novel approaches to species distribution modeling aimed at reducing the influences of such challenges and improving the realism of projections. We estimated species-environment relationships with four modeling methods run with multiple scenarios of (1) sources of occurrences and geographically isolated background ranges for absences, (2) approaches to drawing background (absence) points, and (3) alternate sets of predictor variables. We further tested various quantitative metrics of model evaluation against biological insight. Model projections were very sensitive to the choice of training dataset. Model accuracy was much improved by using a global dataset for model training, rather than restricting data input to the species’ native range. AUC score was a poor metric for model evaluation and, if used alone, was not a useful criterion for assessing model performance. Projections away from the sampled space (i.e. into areas of potential future invasion) were very different depending on the modeling methods used, raising questions about the reliability of ensemble projections. Generalized linear models gave very unrealistic projections far away from the training region. Models that efficiently fit the dominant pattern, but exclude highly local patterns in the dataset and capture interactions as they appear in data (e.g. boosted regression trees), improved generalization of the models. Biological knowledge of the species and its distribution was important in refining choices about the best set of projections. A post-hoc test conducted on a new Partenium dataset from Nepal validated excellent predictive performance of our “best” model. We showed that vast stretches of currently uninvaded geographic areas on multiple continents harbor highly suitable habitats for Parthenium hysterophorus L. (Asteraceae; parthenium). However, discrepancies between model predictions and parthenium invasion in Australia indicate successful management for this globally significant weed. This article is protected by copyright. All rights reserved.

Projecting future expansion of invasive species: Comparing and improving methodologies for species distribution modeling

Modeling the distributions of species, especially of invasive species in non-native ranges, involves multiple challenges. Here, we developed some novel approaches to species distribution modeling aimed at reducing the influences of such challenges and improving the realism of projections. We estimated species-environment relationships with four modeling methods run with multiple scenarios of (1) sources of occurrences and geographically isolated background ranges for absences, (2) approaches to drawing background (absence) points, and (3) alternate sets of predictor variables. We further tested various quantitative metrics of model evaluation against biological insight. Model projections were very sensitive to the choice of training dataset. Model accuracy was much improved by using a global dataset for model training, rather than restricting data input to the species’ native range. AUC score was a poor metric for model evaluation and, if used alone, was not a useful criterion for assessing model performance. Projections away from the sampled space (i.e. into areas of potential future invasion) were very different depending on the modeling methods used, raising questions about the reliability of ensemble projections. Generalized linear models gave very unrealistic projections far away from the training region. Models that efficiently fit the dominant pattern, but exclude highly local patterns in the dataset and capture interactions as they appear in data (e.g. boosted regression trees), improved generalization of the models. Biological knowledge of the species and its distribution was important in refining choices about the best set of projections. A post-hoc test conducted on a new Partenium dataset from Nepal validated excellent predictive performance of our “best” model. We showed that vast stretches of currently uninvaded geographic areas on multiple continents harbor highly suitable habitats for Parthenium hysterophorus L. (Asteraceae; parthenium). However, discrepancies between model predictions and parthenium invasion in Australia indicate successful management for this globally significant weed. This article is protected by copyright. All rights reserved.

Statistical methods for analysis of multi-harvest data from perennial pasture variety selection trials

Variety selection in perennial pasture crops involves identifying best varieties from data collected from multiple harvest times in field trials. For accurate selection, the statistical methods for analysing such data need to account for the spatial and temporal correlation typically present. This paper provides an approach for analysing multi-harvest data from variety selection trials in which there may be a large number of harvest times. Methods are presented for modelling the variety by harvest effects while accounting for the spatial and temporal correlation between observations. These methods provide an improvement in model fit compared to separate analyses for each harvest, and provide insight into variety by harvest interactions. The approach is illustrated using two traits from a lucerne variety selection trial. The proposed method provides variety predictions allowing for the natural sources of variation and correlation in multi-harvest data.