Using harvest statistics to monitor temporal variation in kangaroo density and harvest rate

Management of the commercial harvest of kangaroos relies on quotas set annually as a proportion of regular estimates of population size. Surveys to generate these estimates are expensive and, in the larger states, logistically difficult; a cheaper alternative is desirable. Rainfall is a disappointingly poor predictor of kangaroo rate of increase in many areas, but harvest statistics (sex ratio, carcass weight, skin size and animals shot per unit time) potentially offer cost-effective indirect monitoring of population abundance (and therefore trend) and status (i.e. under-or overharvest). Furthermore, because harvest data are collected continuously and throughout the harvested areas, they offer the promise of more intensive and more representative coverage of harvest areas than aerial surveys do. To be useful, harvest statistics would need to have a close and known relationship with either population size or harvest rate. We assessed this using longterm (11-22 years) data for three kangaroo species (Macropus rufus, M. giganteus and M. fuliginosus) and common wallaroos (M. robustus) across South Australia, New South Wales and Queensland. Regional variation in kangaroo body size, population composition, shooter efficiency and selectivity required separate analyses in different regions. Two approaches were taken. First, monthly harvest statistics were modelled as a function of a number of explanatory variables, including kangaroo density, harvest rate and rainfall. Second, density and harvest rate were modelled as a function of harvest statistics. Both approaches incorporated a correlated error structure. Many but not all regions had relationships with sufficient precision to be useful for indirect monitoring. However, there was no single relationship that could be applied across an entire state or across species. Combined with rainfall-driven population models and applied at a regional level, these relationships could be used to reduce the frequency of aerial surveys without compromising decisions about harvest management.

Stout Whiting Fishery Summary : Commercial Quota Setting for 2017

This fishery assessment report describes the commercial stout whiting fishery operation along Australia’s east coast between Sandy Cape and the Queensland-New South Wales border. The fishery is identified by a T4 symbol. This study follows methods applied in (O'Neill & Leigh, 2016a) and extends the results of that study by using the latest data available up to end of March 2016. The fishery statistics reported herein are for fishing years 1991 to 2016. This study analysed stout whiting catch rates from both Queensland and New South Wales (NSW) for all vessels, areas and fishing gears. The 2016 catch rate index from Queensland and NSW waters was 0.86. This means that the 2016 catch rate index was 86% of the mean standardised catch rate. Results showed that there was a stable trend in catch rates from 2012 to 2016, as in the previous study (O'Neill & Leigh, 2016a), with the 2015 and 2014 catch rates 85% of the mean catch rate. The fish-length frequency and age-length-otolith data were translated using two models which showed: • Where patterns of fish age-abundance were estimated from the fish-length frequency and age-length data, there were slightly decreased estimated measures of fish survival at 38% for 2014, compared to fish survival estimates in 2013 at 40%. The 2014 and 2015 estimated age structure was dominated by 1+ and 2+ old fished, with a slightly higher frequency of age 2 - 3 fish for 2015. • Where only the age-length data were used, estimates showed that from 2011 to 2014 the survival index increased. The estimated survival index increased from 35% in 2013 to 64% in 2014, indicating stronger survival of fish as they recruited and aged. Together the stout whiting catch rate and survival indicators showed the recent fishery harvests were sustainable. Since 1997, T4 management (Stout Whiting Fishery) is centred on annual assessments of total allowable commercial catch (TACC). The TACC is assessed before the start of each fishing year using statistical assessment methodologies, namely evaluation of trends in fish catch rates and catch-at-age frequencies measured against management reference points. The TACC has been under-caught in many years. For setting the 2017 T4 stout whiting TACC, the calculations covered a range of settings to account for the variance in the data and provide options for quota change. The overall (averaged) results suggested: • The procedure where the quota was adjusted based on previous TACC setting in year 2016 gave a recommended TACC for 2017 of between 1100 and 1130 t. • The procedure that focussed directly on optimising the average harvest to match target reference points gave a recommended TACC for 2017 of between 860 and 890 t. Use of these estimates to set TACC will depend on management and industry aims for the fishery.

Stock assessment of the Queensland-New South Wales tailor fishery (Pomatomus saltatrix)

Data on catch sizes, catch rates, length-frequency and age composition from the Australian east coast tailor fishery are analysed by three different population dynamic models: a surplus production model, an age-structured model, and a model in which the population is structured by both age and length. The population is found to be very heavily exploited, with its ability to reproduce dependent on the fishery’s incomplete selectivity of one-year-old fish. Estimates of recent harvest rates (proportion of fish available to the fishery that are actually caught in a single year) are over 80%. It is estimated that only 30–50% of one-year-old fish are available to the fishery. Results from the age-length-structured model indicate that both exploitable biomass (total mass of fish selected by the fishery) and egg production have fallen to about half the levels that prevailed in the 1970s, and about 40% of virgin levels. Two-year-old fish appear to have become smaller over the history of the fishery. This is assumed to be due to increased fishing pressure combined with non-selectivity of small one-year-old fish, whereby the one-year-old fish that survive fishing are small and grow into small two-year-old fish the following year. An alternative hypothesis is that the stock has undergone a genetic change towards smaller fish; the true explanation is unknown. The instantaneous natural mortality rate of tailor is hypothesised to be higher than previously thought, with values between 0.8 and 1.3 yr–1 consistent with the models. These values apply only to tailor up to about three years of age, and it is possible that a lower value applies to fish older than three. The analysis finds no evidence that fishing pressure has yet affected recruitment. If a recruitment downturn were to occur, however, under current management and fishing pressure there is a strong chance that the fishery would need a complete closure for several years to recover, and even then recovery would be uncertain. Therefore it is highly desirable to better protect the spawning stock. The major recommendations are • An increase in the minimum size limit from 30cm to 40cm in order to allow most one-year-old fish to spawn, and • An experiment on discard mortality to gauge the proportion of fish between 30cm and 40cm that are likely to survive being caught and released by recreational line fishers (the dominant component of the fishery, currently harvesting roughly 1000t p.a. versus about 200t p.a. from the commercial fishery).

Inferential, non-parametric statistics to assess quality of probabilistic forecast systems

Many statistical forecast systems are available to interested users. In order to be useful for decision-making, these systems must be based on evidence of underlying mechanisms. Once causal connections between the mechanism and their statistical manifestation have been firmly established, the forecasts must also provide some quantitative evidence of `quality’. However, the quality of statistical climate forecast systems (forecast quality) is an ill-defined and frequently misunderstood property. Often, providers and users of such forecast systems are unclear about what ‘quality’ entails and how to measure it, leading to confusion and misinformation. Here we present a generic framework to quantify aspects of forecast quality using an inferential approach to calculate nominal significance levels (p-values) that can be obtained either by directly applying non-parametric statistical tests such as Kruskal-Wallis (KW) or Kolmogorov-Smirnov (KS) or by using Monte-Carlo methods (in the case of forecast skill scores). Once converted to p-values, these forecast quality measures provide a means to objectively evaluate and compare temporal and spatial patterns of forecast quality across datasets and forecast systems. Our analysis demonstrates the importance of providing p-values rather than adopting some arbitrarily chosen significance levels such as p < 0.05 or p < 0.01, which is still common practice. This is illustrated by applying non-parametric tests (such as KW and KS) and skill scoring methods (LEPS and RPSS) to the 5-phase Southern Oscillation Index classification system using historical rainfall data from Australia, The Republic of South Africa and India. The selection of quality measures is solely based on their common use and does not constitute endorsement. We found that non-parametric statistical tests can be adequate proxies for skill measures such as LEPS or RPSS. The framework can be implemented anywhere, regardless of dataset, forecast system or quality measure. Eventually such inferential evidence should be complimented by descriptive statistical methods in order to fully assist in operational risk management.

Aspects of reproduction and diet of the Australian endemic skate Dipturus polyommata (Ogilby) (Elasmobranchii: Rajidae), by-catch of a commercial prawn trawl fishery

The Australian endemic skate Dipturus polyommata collected from by-catch of a benthic prawn fishery off southern Queensland was examined to provide information on reproduction and diet. Morphological relationships of total length (LT) to disc width and LT to mass were estimated. Size at birth was estimated at c. 100-110 mm and size at first feeding at c. 105-110 mm LT. Size at 50% maturity (LT50 and 95% CI) was 321 (305-332) and 300 (285-306) mm LT for females and males, respectively. Size at first maturity corresponded to 87.7% of observed maximum size in females (366 mm LT) and 87.5% in males (343 mm L T). Two females, representing 18.2% of mature females sampled in the austral winter were each carrying two egg cases. Descriptions of egg cases are given. Diet described by the index of relative importance as a percentage (%IRI) was predominantly crustacean based with carid shrimps (53.64%) and penaeoid prawns (23.30%) the most significant prey groups. Teleosts (11.72%), gammarid amphipods (5.31%) and mysids (4.72%) were also important to the diet of the species, while a further six prey groups made only a minor contribution to diet (1.31%). An ontogenetic change was evident between the diets of immature and mature skates. Immature animals fed more extensively on carids and amphipods and mature animals on penaeoids, teleosts and mysids.

Spatial and temporal variation and effects of changes in management in discard rates from the commercial reef line fishery of the Great Barrier Reef, Australia

Discarding in commercially exploited fisheries has received considerable attention in the last decade, though only more recently in Australia. The Reef Line fishery (RLF) of the Great Barrier Reef (GBR) in Australia is a large-scale multi-sector, multi-species, highly regulated hook and line fishery with the potential for high levels of discarding. We used a range of data sources to estimate discard rates and discard quantities for the two main target groups of the RLF, the coral trout, Plectropomus spp, and the red throat emperor, Lethrinus miniatus, and investigated possible effects on discarding of recent changes in management of the fishery. Fleet-wide estimates of total annual quantities discarded from 1989 to 2003 were 292-622 t and 33-95 t for coral trout and red throat emperor, respectively. Hypothetical scenarios of high-grading after the introduction of a total allowable commercial catch for coral trout resulted in increases in discard quantities up to 3895 t, while no high-grading still meant 421 t were discarded. Increasing the minimum size limit of red throat emperor from 35 to 38 cm also increased discards to an estimated 103 t. We provide spatially and temporally explicit estimates of discarding for the two most important species in the GBR RLF of Australia to demonstrate the importance of accounting for regional variation in quantification of discarding. Effects of management changes on discarding are also highlighted. This study provides a template for exploring discarding levels for other species in the RLF and elsewhere.

Frequency and precision of aerial surveys for kangaroo management

The appropriate frequency and precision for surveys of wildlife populations represent a trade-off between survey cost and the risk of making suboptimal management decisions because of poor survey data. The commercial harvest of kangaroos is primarily regulated through annual quotas set as proportions of absolute estimates of population size. Stochastic models were used to explore the effects of varying precision, survey frequency and harvest rate on the risk of quasiextinction for an arid-zone and a more mesic-zone kangaroo population. Quasiextinction probability increases in a sigmoidal fashion as survey frequency is reduced. The risk is greater in more arid regions and is highly sensitive to harvest rate. An appropriate management regime involves regular surveys in the major harvest areas where harvest rate can be set close to the maximum sustained yield. Outside these areas, survey frequency can be reduced in relatively mesic areas and reduced in arid regions when combined with lowered harvest rates. Relative to other factors, quasiextinction risk is only affected by survey precision (standard error/mean × 100) when it is >50%, partly reflecting the safety of the strategy of harvesting a proportion of a population estimate.

A comparison of mark-recapture distance-sampling methods applied to aerial surveys of eastern grey kangaroos

Aerial surveys of kangaroos (Macropus spp.) in Queensland are used to make economically important judgements on the levels of viable commercial harvest. Previous analysis methods for aerial kangaroo surveys have used both mark-recapture methodologies and conventional distance-sampling analyses. Conventional distance sampling has the disadvantage that detection is assumed to be perfect on the transect line, while mark-recapture methods are notoriously sensitive to problems with unmodelled heterogeneity in capture probabilities. We introduce three methodologies for combining together mark-recapture and distance-sampling data, aimed at exploiting the strengths of both methodologies and overcoming the weaknesses. Of these methods, two are based on the assumption of full independence between observers in the mark-recapture component, and this appears to introduce more bias in density estimation than it resolves through allowing uncertain trackline detection. Both of these methods give lower density estimates than conventional distance sampling, indicating a clear failure of the independence assumption. The third method, termed point independence, appears to perform very well, giving credible density estimates and good properties in terms of goodness-of-fit and percentage coefficient of variation. Estimated densities of eastern grey kangaroos range from 21 to 36 individuals km-2, with estimated coefficients of variation between 11% and 14% and estimated trackline detection probabilities primarily between 0.7 and 0.9.

Analysis of high yielding maize production - A study based on a commercial crop

This paper reports on the use of APSIM - Maize for retrospective analysis of performance of a high input, high yielding maize crop and analysis of predicted performance of maize grown with high inputs over the long-term (>100 years) for specified scenarios of environmental conditions (temperature and radiation) and agronomic inputs (sowing date, plant population, nitrogen fertiliser and irrigation) at Boort, Victoria, Australia. It uses a high yielding (17 400 kg/ha dry grain, 20 500 kg/ha at 15% water) commercial crop grown in 2004-05 as the basis of the study. Yield for the agronomic and environmental conditions of 2004-05 was predicted accurately, giving confidence that the model could be used for the detailed analyses undertaken. The analysis showed that the yield achieved was close to that possible with the conditions and agronomic inputs of 2004-05. Sowing dates during 21 September to 26 October had little effect on predicted yield, except when combined with reduced temperature. Single year and long-term analyses concluded that a higher plant population (11 plants/m2) is needed to optimise yield, but that slightly lower N and irrigation inputs are appropriate for the plant population used commercially (8.4 plants/m2). Also, compared with changes in agronomic inputs increases in temperature and/or radiation had relatively minor effects, except that reduced temperature reduces predicted yield substantially. This study provides an approach for the use of models for both retrospective analysis of crop performance and assessment of long-term variability of crop yield under a wide range of agronomic and environmental conditions.

Concentrations of constitutive alk(en)ylresorcinols in peel of commercial mango varieties and resistance to postharvest anthracnose

Mature green mango fruits of commercially important varieties were screened to investigate the levels of constitutive antifungal compounds in peel and to assess anthracnose disease after inoculation with Colletotrichum gloeosporioides. High pressure liquid chromatography was used to quantify the levels of 5-n-heptadecenylresorcinol and 5-n-pentadecylresorcinol in the peel extracts. The fruit peel of the varieties ‘Kensington Pride’ and ‘Keitt’ were observed to have the highest levels of both 5-n-heptadecenylresorcinol (107.3-123.7 and 49.9-61.4 μg/g FW, respectively) and 5-n-pentadecylresorcinol (6.32-7.99 and 3.30-6.05 μg/g FW, respectively), and the fruit of the two varieties were found to have some resistance to postharvest anthracnose. The varieties ‘Kent’, ‘R2E2’, ‘Nam Doc Mai’, ‘Calypso’, and ‘Honey Gold’ contained much lower concentrations of resorcinols in their peel and three of these varieties were found to be more susceptible to anthracnose. Concentrations of 5-nheptadecenylresorcinol were significantly lower at the ‘sprung’ and ‘eating ripe’ stages of ripening compared to levels at harvest. Concentrations of 5-n-pentadecylresorcinol did not differ significantly across the three stages of ripening. The levels of these two resorcinols were found to be strongly inter-correlated (P < 0.001, r2 = 0.71), with concentrations of 5-nheptadecenylresorcinol being an average 18 times higher than those of 5-npentadecylresorcinol. At the ‘eating ripe’ stage, significant relationships were observed between the concentrations of each type of alk(en)ylresorcinol and anthracnose lesion areas following postharvest inoculation, P<0.001, r2= 0.69 for 5-n pentadecylresorcinol, and P<0.001, r2= 0.44 for 5-n-heptadecenylresorcinol.

From inferential statistics to climate knowledge

Climate variability and change are risk factors for climate sensitive activities such as agriculture. Managing these risks requires "climate knowledge", i.e. a sound understanding of causes and consequences of climate variability and knowledge of potential management options that are suitable in light of the climatic risks posed. Often such information about prognostic variables (e.g. yield, rainfall, run-off) is provided in probabilistic terms (e.g. via cumulative distribution functions, CDF), whereby the quantitative assessments of these alternative management options is based on such CDFs. Sound statistical approaches are needed in order to assess whether difference between such CDFs are intrinsic features of systems dynamics or chance events (i.e. quantifying evidences against an appropriate null hypothesis). Statistical procedures that rely on such a hypothesis testing framework are referred to as "inferential statistics" in contrast to descriptive statistics (e.g. mean, median, variance of population samples, skill scores). Here we report on the extension of some of the existing inferential techniques that provides more relevant and adequate information for decision making under uncertainty.

Effect of micro-propagation on the health status of strawberry planting material for commercial production of strawberry runners for Queensland

Plant tissue culture has been used for a number of years to produce micropropagated strawberry plants for planting into runner growing beds in the Stanthorpe (Queensland) and Bothwell (Tasmania) regions. This process has allowed the rapid release of new cultivars from the LAWS (Late Autumn, Winter, Spring) breeding program into the current runner production system. Micro-propagation in vitro allows plants to be produced during the autumn and winter months, when mother plants would normally be in a fruit production phase in the field in Queensland. The plants produced are of a high health status when they are planted. The subsequent arrival and build up of various diseases in the runner fields are closely monitored. Using tissue culture for the first generation reduces the time the plants spend in the field by twelve months, reducing disease incidence. To date, any disease outbreak has been successfully managed using early detection and rapid response methods.

Herpesviral haematopoietic necrosis virus (CyHV-2) infection: case studies from commercial goldfish farms

Herpesviral haematopoietic necrosis is a disease of goldfish, Carassius auratus, caused by Cyprinid herpesvirus-2 (CyHV-2) infection. Quantitative PCR was carried out on tissue homogenates from healthy goldfish fingerlings, broodfish, eggs and fry directly sampled from commercial farms, from moribund fish submitted to our laboratory for disease diagnosis, and on naturally-infected CyHV-2 carriers subjected to experimental stress treatments. Healthy fish from 14 of 18 farms were positive with copy numbers ranging from tens to 10(7) copies mu g(-1) DNA extracted from infected fish. Of 118 pools of broodfish tested, 42 were positive. The CyHV-2 was detected in one lot of fry produced from disinfected eggs. Testing of moribund goldfish, in which we could not detect any other pathogens, produced 12 of 30 cases with 10(6)-10(8) copies of CyHV-2 mu g(-1) DNA extracted. Subjecting healthy CyHV-2 carriers to cold shock (22-10 degrees C) but not heat, ammonia or high pH, increased viral copy numbers from mean copy number (+/- SE) of 7.3 +/- 11 to 394 +/- 55 mu g(-1) DNA extracted after 24 h. CyHV-2 is widespread on commercial goldfish farms and outbreaks apparently occur when healthy carriers are subjected to a sharp temperature drop followed by holding at the permissive temperature for the disease.

Quantifying Transmission of Campylobacter jejuni in Commercial Broiler Flocks

Since meat from poultry colonized with Campylobacter spp. is a major cause of bacterial gastroenteritis, human exposure should be reduced by, among other things, prevention of colonization of broiler flocks. To obtain more insight into possible sources of introduction of Campylobacter into broiler flocks, it is essential to estimate the moment that the first bird in a flock is colonized. If the rate of transmission within a flock were known, such an estimate could be determined from the change in the prevalence of colonized birds in a flock over time. The aim of this study was to determine the rate of transmission of Campylobacter using field data gathered for 5 years for Australian broiler flocks. We used unique sampling data for 42 Campylobacter jejuni-colonized flocks and estimated the transmission rate, which is defined as the number of secondary infections caused by one colonized bird per day. The estimate was 2.37 +/- 0.295 infections per infectious bird per day, which implies that in our study population colonized flocks consisting of 20,000 broilers would have an increase in within-flock prevalence to 95% within 4.4 to 7.2 days after colonization of the first broiler. Using Bayesian analysis, the moment of colonization of the first bird in a flock was estimated to be from 21 days of age onward in all flocks in the study. This study provides an important quantitative estimate of the rate of transmission of Campylobacter in broiler flocks, which could be helpful in future studies on the epidemiology of Campylobacter in the field.

Development and testing of an evaluation procedure for commercial manure additive products.

Manure additive products can be used to reduce odour emissions (OE) from livestock farms. The standardised evaluation of these manure additive products under specific farm conditions is important. In this study, the efficacy of a manure additive (WonderTreat(TM), CKLS, Inc., Hong-Kong) was assessed under Australian conditions utilising a combination of laboratory and field-scale evaluation techniques. As a first step, the efficacy of the manure additive was assessed in a laboratory-scale trial using a series of uniformly managed digesters and standard odour, liquor ammonia and hydrogen sulphide concentration measurement procedures. This showed that the addition of WonderTreat(TM) at the 'low dose rate' (LDR) (102.6 g m-2) used during the trial significantly, but only marginally (30%; P = 0.02) reduced the OE rate (mean 13.9 OU m-2 s-1) of anaerobic pig liquor relative to an untreated control (UC) (19.9 OU m-2 s-1). However, the 'high dose rate' (HDR) (205.3 g m-2) also assessed during the trial preformed similarly (19.7 OU m-2 s-1) to the UC. No statistically significant difference in the concentrations of a range of measured water quality variables at the 5% level was observed between the treatments or controls digesters. As a second step, a field-scale assessment of the manure additive was undertaken at a commercial piggery. Two piggery manure lagoons (each with approximately 2500 m2 surface area) were included in the study; one was treated with WonderTreat(TM) while the other was used as control. The efficacy of the treatment was assessed using olfactometric evaluation of odour samples collected from the surface of the pond using a dynamic wind tunnel and ancillary equipment. No statistically significant reduction in OE rate could be demonstrated (P = 0.35), partially due to the limited number of samples taken during the assessment. However, there was a numerical reduction in the average OE rate of the treatment pond (29 OU m-2 s-1 at 1 m s-1) compared to the control lagoon (38 OU m-2 s-1 at 1 m s-1).