Predictive modelling to aid the regional-scale management of a vertebrate pest

Extensive resources are allocated to managing vertebrate pests, yet spatial understanding of pest threats, and how they respond to management, is limited at the regional scale where much decision-making is undertaken. We provide regional-scale spatial models and management guidance for European rabbits (Oryctolagus cuniculus) in a 260,791 km(2) region in Australia by determining habitat suitability, habitat susceptibility and the effects of the primary rabbit management options (barrier fence, shooting and baiting and warren ripping) or changing predation or disease control levels. A participatory modelling approach was used to develop a Bayesian network which captured the main drivers of suitability and spread, which in turn was linked spatially to develop high resolution risk maps. Policy-makers, rabbit managers and technical experts were responsible for defining the questions the model needed to address, and for subsequently developing and parameterising the model. Habitat suitability was determined by conditions required for warren-building and by above-ground requirements, such as food and harbour, and habitat susceptibility by the distance from current distributions, habitat suitability, and the costs of traversing habitats of different quality. At least one-third of the region had a high probability of being highly suitable (support high rabbit densities), with the model supported by validation. Habitat susceptibility was largely restricted by the current known rabbit distribution. Warren ripping was the most effective control option as warrens were considered essential for rabbit persistence. The anticipated increase in disease resistance was predicted to increase the probability of moderately suitable habitat becoming highly suitable, but not increase the at-risk area. We demonstrate that it is possible to build spatial models to guide regional-level management of vertebrate pests which use the best available knowledge and capture fine spatial-scale processes.

Spatial model of site index based on Y-ray spectrometry and a digital elevation model for two Pinus species in Tuan Toolara State Forest, Queensland, Australia

Site index prediction models are an important aid for forest management and planning activities. This paper introduces a multiple regression model for spatially mapping and comparing site indices for two Pinus species (Pinus elliottii Engelm. and Queensland hybrid, a P. elliottii x Pinus caribaea Morelet hybrid) based on independent variables derived from two major sources: g-ray spectrometry (potassium (K), thorium (Th), and uranium (U)) and a digital elevation model (elevation, slope, curvature, hillshade, flow accumulation, and distance to streams). In addition, interpolated rainfall was tested. Species were coded as a dichotomous dummy variable; interaction effects between species and the g-ray spectrometric and geomorphologic variables were considered. The model explained up to 60% of the variance of site index and the standard error of estimate was 1.9 m. Uranium, elevation, distance to streams, thorium, and flow accumulation significantly correlate to the spatial variation of the site index of both species, and hillshade, curvature, elevation and slope accounted for the extra variability of one species over the other. The predicted site indices varied between 20.0 and 27.3 m for P. elliottii, and between 23.1 and 33.1 m for Queensland hybrid; the advantage of Queensland hybrid over P. elliottii ranged from 1.8 to 6.8 m, with the mean at 4.0 m. This compartment-based prediction and comparison study provides not only an overview of forest productivity of the whole plantation area studied but also a management tool at compartment scale.

The effect of data sources and quality on the predictive capacity of CLIMEX models: an assessment of Teleonemia scrupulosa and Octotoma scabripennis for the biocontrol of Lantana camara in Australia

Understanding the effects of different types and quality of data on bioclimatic modeling predictions is vital to ascertaining the value of existing models, and to improving future models. Bioclimatic models were constructed using the CLIMEX program, using different data types – seasonal dynamics, geographic (overseas) distribution, and a combination of the two – for two biological control agents for the major weed Lantana camara L. in Australia. The models for one agent, Teleonemia scrupulosa Stål (Hemiptera:Tingidae) were based on a higher quality and quantity of data than the models for the other agent, Octotoma scabripennis Guérin-Méneville (Coleoptera: Chrysomelidae). Predictions of the geographic distribution for Australia showed that T. scrupulosa models exhibited greater accuracy with a progressive improvement from seasonal dynamics data, to the model based on overseas distribution, and finally the model combining the two data types. In contrast, O. scabripennis models were of low accuracy, and showed no clear trends across the various model types. These case studies demonstrate the importance of high quality data for developing models, and of supplementing distributional data with species seasonal dynamics data wherever possible. Seasonal dynamics data allows the modeller to focus on the species response to climatic trends, while distributional data enables easier fitting of stress parameters by restricting the species envelope to the described distribution. It is apparent that CLIMEX models based on low quality seasonal dynamics data, together with a small quantity of distributional data, are of minimal value in predicting the spatial extent of species distribution.

Reassessing spatial and temporal dynamics of kangaroo populations

Long-running datasets from aerial surveys of kangaroos (Macropus giganteus, Macropus [uliginosus, Macropus robustus and Macropus rufus) across Queensland, New South Wales and South Australia have been analysed, seeking better predictors of rates of increase which would allow aerial surveys to be undertaken less frequently than annually. Early models of changes in kangaroo numbers in response to rainfall had shown great promise, but much variability. We used normalised difference vegetation index (NDVI) instead, reasoning that changes in pasture condition would provide a better predictor than rainfall. However, except at a fine scale, NDVI proved no better; although two linked periods of rainfall proved useful predictors of rates of increase, this was only in some areas for some species. The good correlations reported in earlier studies were a consequence of data dominated by large droughtinduced adult mortality, whereas over a longer time frame and where changes between years are less dramatic, juvenile survival has the strongest influence on dynamics. Further, harvesting, density dependence and competition with domestic stock are additional and important influences and it is now clear that kangaroo movement has a greater influence on population dynamics than had been assumed. Accordingly, previous conclusions about kangaroo populations as simple systems driven by rainfall need to be reassessed. Examination of this large dataset has permitted descriptions of shifts in distribution of three species across eastern Australia, changes in dispersion in response to rainfall, and an evaluation of using harvest statistics as an index of density and harvest rate. These results have been combined into a risk assessment and decision theory framework to identify optimal monitoring strategies.

Spatial variation in life history reveals insight into connectivity and geographic population structure of a tropical estuarine teleost: king threadfin, Polydactylus macrochir

Understanding the life history of exploited fish species is not only critical in developing stock assessments and productivity models, but has a dual function in the delineation of connectivity and geographical population structure. In this study, patterns in growth and length and age at sex change of Polydactylus macrochir, an ecologically and economically important protandrous estuarine teleost, were examined to provide preliminary information on the species' connectivity and geographic structure across northern Australia. Considerable variation in life history parameters was observed among the 18 locations sampled. Both unconstrained and constrained (t(0) = 0) estimates of von Bertalanffy growth function parameters differed significantly among all neighbouring locations with the exception of two locations in Queensland's east coast and two in Queensland's Gulf of Carpentaria waters, respectively. Comparisons of back-calculated length-at-age 2 provided additional evidence for growth differences among some locations, but were not significantly different among locations in the south-eastern Gulf of Carpentaria or on Queensland's east coast. The length and age at sex change differed markedly among locations, with fish from the east coast of Australia changing sex from males to females at significantly greater lengths and ages than elsewhere. Sex change occurred earliest at locations within Queensland's Gulf of Carpentaria, where a large proportion of small, young females were recorded. The observed differences suggest that P. macrochir likely form a number of geographically and/or reproductively distinct groups in Australian waters and suggest that future studies examining connectivity and geographic population structure of estuarine fishes will likely benefit from the inclusion of comparisons of life history parameters. (C) 2012 Elsevier B.V. All rights reserved.

Allocasuarina tree hosts determine the spatial distribution of hypogeous fungal sporocarps in three tropical Australian sclerophyll forests

Across three tropical Australian sclerophyll forest types, site-specific environmental variables could explain the distribution of both quantity (abundance and biomass) and richness (genus and species) of hypogeous fungi sporocarps. Quantity was significantly higher in the Allocasuarina forest sites that had high soil nitrogen but low phosphorous. Three genera of hypogeous fungi were found exclusively in Allocasuarina forest sites including Gummiglobus, Labyrinthomyces and Octaviania, as were some species of Castoreum, Chondrogaster, Endogone, Hysterangium and Russula. However, the forest types did not all group according to site-scale variables and subsequently the taxonomic assemblages were not significantly different between the three forest types. At site scale, significant negative relationships were found between phosphorous concentration and the quantity of hypogeous fungi sporocarps. Using a multivariate information theoretic approach, there were other more plausible models to explain the patterns of sporocarp richness. Both the mean number of fungal genera and species increased with the number of Allocasuarina stems, at the same time decreasing with the number of Eucalyptus stems. The optimal conditions for promoting hypogeous fungi sporocarp quantity and sporocarp richness appear to be related to the presence and abundance of Allocasuarina (Casuarinaceae) host trees. Allocasuarina tree species may have a higher host receptivity for ectomycorrhizal hypogeous fungi species that provide an important food resource for Australian mycophagous animals.

Spatial impact of projected changes in rainfall and temperature on wheat yields in Australia

Climate projections over the next two to four decades indicate that most of Australia’s wheat-belt is likely to become warmer and drier. Here we used a shire scale, dynamic stress-index model that accounts for the impacts of rainfall and temperature on wheat yield, and a range of climate change projections from global circulation models to spatially estimate yield changes assuming no adaptation and no CO2 fertilisation effects. We modelled five scenarios, a baseline climate (climatology, 1901–2007), and two emission scenarios (“low” and “high” CO2) for two time horizons, namely 2020 and 2050. The potential benefits from CO2 fertilisation were analysed separately using a point level functional simulation model. Irrespective of the emissions scenario, the 2020 projection showed negligible changes in the modelled yield relative to baseline climate, both using the shire or functional point scale models. For the 2050-high emissions scenario, changes in modelled yield relative to the baseline ranged from −5 % to +6 % across most of Western Australia, parts of Victoria and southern New South Wales, and from −5 to −30 % in northern NSW, Queensland and the drier environments of Victoria, South Australia and in-land Western Australia. Taking into account CO2 fertilisation effects across a North–south transect through eastern Australia cancelled most of the yield reductions associated with increased temperatures and reduced rainfall by 2020, and attenuated the expected yield reductions by 2050.