Baseline responses of Alphitobius diaperinus (Coleoptera : Tenebrionidae) to cyfluthrin and detection of strong resistance in field populations in eastern Australia

Resistance to cyfluthrin in broiler farm populations of lesser mealworm, Alphitobius diaperinus (Panzer) (Coleoptera: Tenebrionidae), in eastern Australia was suspected to have contributed to recent control failures. In 2000-2001, beetles from 11 broiler farms were tested for resistance by comparing them to an insecticide-susceptible reference population by using topical application. Resistance was detected in almost all beetle populations (up to 22 times the susceptible at the LC50), especially in southeastern Queensland where more cyfluthrin applications had been made. Two from outside southeastern Queensland were found to be susceptible. Dose-mortality data generated from the reference population over a range of cyflutbrin concentrations showed that 0.0007% cyfluthrin at a LC99.9 level could be used as a convenient dose to discriminate between susceptible and resistant populations. Using this discriminating concentration, from 2001 to 2005, the susceptibilities of 18 field populations were determined. Of these, 11 did not exhibit complete mortality at the discriminating concentration (mortality range 2.8-97.7%), and in general, cyfluthrin resistance was directly related to the numbers of cyfluthrin applications. As in the full study, populations outside of southeastern Queensland were found to have lower levels of resistance or were susceptible. One population from an intensively farmed broiler area in southeastern Queensland exhibited low mortality despite having no known exposure to cyfluthrin. Comparisons of LC50 values of three broiler populations and a susceptible population, collected in 2000 and 2001 and recollected in 2004 and 2005 indicated that values from the three broiler populations had increased over this time for all populations. The continued use of cyfluthrin for control of A. diaperinus in eastern Australia is currently under consideration.

Comparison of bioassay responses to the potential fungal biopesticide Metarhizium anisopliae in Rhipicephalus(Boophilus) microplus and Lucilia cuprina.

Quantal response bioassays were conducted with cattle ticks and sheep blowflies with three different isolates of Metarhizium anisopliae and different methods of inoculation. Ticks were either topically dosed with 2 mu l or immersed in the conidial preparations. Blowflies were either topically dosed with 2 mu l of the conidial preparation or fed on conidia mixed with sugar. Probit analyses were carried out on the mortality data to compare the virulence of these isolates to ticks and blowflies and look for indications of different virulence mechanisms employed by M. anisopliae isolates when invading these hosts. One isolate (ARIM16) showed high virulence to both hosts killing 95% of ticks after 2 days and 88 (+/- 2)% of blowflies after 4 days. Strikingly different mortality patterns indicated that virulence is dependent on different mechanisms in ticks and blowflies. The pattern of mortality seen with ticks suggested that the number of conidia adhering per unit area of the cuticle was more important for rapid tick death than the total number of conidia contacting the entire tick surface. Blowflies fed conidia mixed with food died rapidly after an initial lag phase regardless of dose.

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

Prospects for predicting insect mortality in relation to changing phosphine concentrations

Typically, in bag-stack or silo fumigations the concentration of phosphine is not constant, and yet most of what is known about phosphine efficacy against grain insects comes from studies with fixed concentrations. Indeed, where changing concentration experiments have been performed, researchers have been unable to explain observed efficacy on the basis of data from fixed concentrations. The ability to predict insect mortality in relation to changing phosphine concentrations would facilitate the development of effective fumigation protocols. In this paper, we explore the prospects for making such predictions. After reviewing published and new results, we conclude that the commonly used concentration x time (Ct) product is unreliable for this purpose. New results, for a strongly resistant strain of Rhyzopertha dominica from Australia, suggest that the relationship Cnt = k may be useful for predicting mortality of this type of insect in changing concentrations. However, in the case of a strain of Sitophilus oryzae with a type of resistance common in Australian S. oryzae, the relationship Cnt = k proved to be less reliable.

Predicting mortality of phosphine-resistant adults of Sitophilus oryzae (L) (Coleoptera: Curculionidae) in relation to changing phosphine concentration

Adults of a phosphine-resistant strain of Sitophilus oryzae (L) were exposed to constant phosphine concentrations of 0.0035-0.9 mg litre(-1) for periods of between 20 and 168 h at 25 °C, and the effects of time and concentration on mortality were quantified. Adults were also exposed to a series of treatments lasting 48, 72 or 168 h at 25 °C, during which the concentration of phosphine was varied. The aim of this study was to determine whether equations from experiments using constant concentrations could be used to predict the efficacy of changing phosphine concentrations against adults of S oryzae. A probit plane without interaction, in which the logarithms of time (t) and concentration (C) were variables, described the effects of concentration and time on mortality in experiments with constant concentrations. A derived equation of the form C^nt = k gave excellent predictions of toxicity when applied to data from changing concentration experiments. The results suggest that for resistant S oryzae adults there is nothing inherently different between constant and changing concentration regimes, and that data collected from fixed concentrations can be used to develop equations for predicting mortality in fumigations in which phosphine concentration changes. This approach could simplify the prediction of efficacy of typical fumigations in which concentrations tend to rise and then fall over a period of days.

Modelling Mortality and Reproduction Rates for Management of Sheep Flocks in Northern Australia

By quantifying the effects of climatic variability in the sheep grazing lands of north western and western Queensland, the key biological rates of mortality and reproduction can be predicted for sheep. These rates are essential components of a decision support package which can prove a useful management tool for producers, especially if they can easily obtain the necessary predictors. When the sub-models of the GRAZPLAN ruminant biology process model were re-parameterised from Queensland data along with an empirical equation predicting the probability of ewes mating added, the process model predicted the probability of pregnancy well (86% variation explained). Predicting mortality from GRAZPLAN was less successful but an empirical equation based on relative condition of the animal (a measure based on liveweight), pregnancy status and age explained 78% of the variation in mortalities. A crucial predictor in these models was liveweight which is not often recorded on producer properties. Empirical models based on climatic and pasture conditions estimated from the pasture production model GRASP, predicted marking and mortality rates for Mitchell grass (Astrebla sp.) pastures (81% and 63% of the variation explained). These prediction equations were tested against independent data from producer properties and the model successfully validated for Mitchell grass communities.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Thermal requirements, field mortality and population phenology modelling of Paropsis atomaria Olivier, an emergent pest in subtropical hardwood plantations

Paropsis atomaria is a recently emerged pest of eucalypt plantations in subtropical Australia. Its broad host range of at least 20 eucalypt species and wide geographical distribution provides it the potential to become a serious forestry pest both within Australia and, if accidentally introduced, overseas. Although populations of P. atomaria are genetically similar throughout its range, population dynamics differ between regions. Here, we determine temperature-dependent developmental requirements using beetles sourced from temperate and subtropical zones by calculating lower temperature thresholds, temperature-induced mortality, and day-degree requirements. We combine these data with field mortality estimates of immature life stages to produce a cohort-based model, ParopSys, using DYMEX™ that accurately predicts the timing, duration, and relative abundance of life stages in the field and number of generations in a spring–autumn (September–May) field season. Voltinism was identified as a seasonally plastic trait dependent upon environmental conditions, with two generations observed and predicted in the Australian Capital Territory, and up to four in Queensland. Lower temperature thresholds for development ranged between 4 and 9 °C, and overall development rates did not differ according to beetle origin. Total immature development time (egg–adult) was approximately 769.2 ± S.E. 127.8 DD above a lower temperature threshold of 6.4 ± S.E. 2.6 °C. ParopSys provides a basic tool enabling forest managers to use the number of generations and seasonal fluctuations in abundance of damaging life stages to estimate the pest risk of P. atomaria prior to plantation establishment, and predict the occurrence and duration of damaging life stages in the field. Additionally, by using local climatic data the pest potential of P. atomaria can be estimated to predict the risk of it establishing if accidentally introduced overseas. Improvements to ParopSys’ capability and complexity can be made as more biological data become available.

Time-dependent instantaneous mortality rates from multiple tagging experiments with exact times of release and recovery

Instantaneous natural mortality rates and a nonparametric hunting mortality function are estimated from a multiple-year tagging experiment with arbitrary, time-dependent fishing or hunting mortality. Our theory allows animals to be tagged over a range of times in each year, and to take time to mix into the population. Animals are recovered by hunting or fishing, and death events from natural causes occur but are not observed. We combine a long-standing approach based on yearly totals, described by Brownie et al. (1985, Statistical Inference from Band Recovery Data: A Handbook, Second edition, United States Fish and Wildlife Service, Washington, Resource Publication, 156), with an exact-time-of-recovery approach originated by Hearn, Sandland and Hampton (1987, Journal du Conseil International pour l'Exploration de la Mer, 43, 107-117), who modeled times at liberty without regard to time of tagging. Our model allows for exact times of release and recovery, incomplete reporting of recoveries, and potential tag shedding. We apply our methods to data on the heavily exploited southern bluefin tuna (Thunnus maccoyii).

Prediction of mortality and conception rates of beef breeding cattle in northern Australia

In current simulation packages for the management of extensive beef-cattle enterprises, the relationships for the key biological rates (namely conception and mortality) are quite rudimentary. To better estimate these relationships, cohort-level data covering 17 100 cow-years from six sites across northern Australia were collated and analysed. Further validation data, from 7200 cow-years, were then used to test these relationships. Analytical problems included incomplete and non-standardised data, considerable levels of correlation among the 'independent' variables, and the close similarity of alternate possible models. In addition to formal statistical analyses of these data, the theoretical equations for predicting mortality and conception rates in the current simulation models were reviewed, and then reparameterised and recalibrated where appropriate. The final models explained up to 80% of the variation in the data. These are now proposed as more accurate and useful models to be used in the prediction of biological rates in simulation studies for northern Australia. © The State of Queensland (through the Department of Agriculture, Fisheries and Forestry) 2012. © CSIRO.

Factors associated with calf mortality in tropically adapted beef breeds managed in extensive Australian production systems

Data from 9296 calves born to 2078 dams over 9 years across five sites were used to investigate factors associated with calf mortality for tropically adapted breeds (Brahman and Tropical Composite) recorded in extensive production systems, using multivariate logistic regression. The average calf mortality pre-weaning was 9.5% of calves born, varying from 1.5% to 41% across all sites and years. In total, 67% of calves that died did so within a week of their birth, with cause of death most frequently recorded as unknown. The major factors significantly (P < 0.05) associated with mortality for potentially large numbers of calves included the specific production environment represented by site-year, low calf birthweight (more so than high birthweight) and horn status at branding. Almost all calf deaths post-branding (assessed from n = 8348 calves) occurred in calves that were dehorned, totalling 2.1% of dehorned calves and 15.9% of all calf deaths recorded. Breed effects on calf mortality were primarily the result of breed differences in calf birthweight and, to a lesser extent, large teat size of cows; however, differences in other breed characteristics could be important. Twin births and calves assisted at birth had a very high risk of mortality, but <1% of calves were twins and few calves were assisted at birth. Conversely, it could not be established how many calves would have benefitted from assistance at birth. Cow age group and outcome from the previous season were also associated with current calf mortality; maiden or young cows (<4 years old) had increased calf losses overall. More mature cows with a previous outcome of calf loss were also more likely to have another calf loss in the subsequent year, and this should be considered for culling decisions. Closer attention to the management of younger cows is warranted to improve calf survival.

Bycatch and strandings programs as ecological indicators for data-limited cetaceans

An integrated approach of using strandings and bycatch data may provide an indicator of long-term trends for data-limited cetaceans. Strandings programs can give a faithful representation of the species composition of cetacean assemblages, while standardised bycatch rates can provide a measure of relative abundance. Comparing the two datasets may also facilitate managing impacts by understanding which species, sex or sizes are the most vulnerable to interactions with fisheries gear. Here we apply this approach to two long-term datasets in East Australia, bycatch in the Queensland Shark Control Program QSCP, 1992–2012) and strandings in the Queensland Marine Wildlife Strandings and Mortality Program StrandNet, 1996–2012). Short-beaked common dolphins, Delphinus delphis, were markedly more frequent in bycatch than in the strandings dataset, suggesting that they are more prone to being incidentally caught than other cetacean species in the region. The reverse was true for humpback whales, Megaptera novaeangliae, bottlenose dolphins, Tursiops spp.; and species predominantly found in offshore waters. QSCP bycatch was strongly skewed towards females for short-beaked common dolphins, and towards smaller sizes for Australian humpback dolphins, Sousa sahulensis. Overall, both datasets demonstrated similar seasonality and a similar long-term increase from 1996 until 2008. Analysis on a species-by-species basis was then used to explore potential explanations for long-term trends, which ranged from a recovering stock (humpback whales) to a shift in habitat use (short-beaked common dolphins).

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 .

Bycatch and strandings programs as ecological indicators for data-limited cetaceans

An integrated approach of using strandings and bycatch data may provide an indicator of long-term trends for data-limited cetaceans. Strandings programs can give a faithful representation of the species composition of cetacean assemblages, while standardised bycatch rates can provide a measure of relative abundance. Comparing the two datasets may also facilitate managing impacts by understanding which species, sex or sizes are the most vulnerable to interactions with fisheries gear. Here we apply this approach to two long-term datasets in East Australia, bycatch in the Queensland Shark Control Program (QSCP, 1992–2012) and strandings in the Queensland Marine Wildlife Strandings and Mortality Program (StrandNet, 1996–2012). Short-beaked common dolphins, Delphinus delphis, were markedly more frequent in bycatch than in the strandings dataset, suggesting that they are more prone to being incidentally caught than other cetacean species in the region. The reverse was true for humpback whales, Megaptera novaeangliae, bottlenose dolphins, Tursiops spp.; and species predominantly found in offshore waters. QSCP bycatch was strongly skewed towards females for short-beaked common dolphins, and towards smaller sizes for Australian humpback dolphins, Sousa sahulensis. Overall, both datasets demonstrated similar seasonality and a similar long-term increase from 1996 until 2008. Analysis on a species-by-species basis was then used to explore potential explanations for long-term trends, which ranged from a recovering stock (humpback whales) to a shift in habitat use (short-beaked common dolphins).