Evaluating and validating abundance monitoring methods in the absence of populations of known size: review and application to a passive tracking index

Rarely is it possible to obtain absolute numbers in free-ranging populations and although various direct and indirect methods are used to estimate abundance, few are validated against populations of known size. In this paper, we apply grounding, calibration and verification methods, used to validate mathematical models, to methods of estimating relative abundance. To illustrate how this might be done, we consider and evaluate the widely applied passive tracking index (PTI) methodology. Using published data, we examine the rationality of PTI methodology, how conceptually animal activity and abundance are related and how alternative methods are subject to similar biases or produce similar abundance estimates and trends. We then attune the method against populations representing a range of densities likely to be encountered in the field. Finally, we compare PTI trends against a prediction that adjacent populations of the same species will have similar abundance values and trends in activity. We show that while PTI abundance estimates are subject to environmental and behavioural stochasticity peculiar to each species, the PTI method and associated variance estimate showed high probability of detection, high precision of abundance values and, generally, low variability between surveys, and suggest that the PTI method applied using this procedure and for these species provides a sensitive and credible index of abundance. This same or similar validation approach can and should be applied to alternative relative abundance methods in order to demonstrate their credibility and justify their use.

Modelling the distribution and relative abundance of feral camels in the Northern Territory using count data

The objectives of this study were to predict the potential distribution, relative abundance and probability of habitat use by feral camels in southern Northern Territory. Aerial survey data were used to model habitat association. The characteristics of ‘used’ (where camels were observed) v. ‘unused’ (pseudo-absence) sites were compared. Habitat association and abundance were modelled using generalised additive model (GAM) methods. The models predicted habitat suitability and the relative abundance of camels in southern Northern Territory. The habitat suitability maps derived in the present study indicate that camels have suitable habitat in most areas of southern Northern Territory. The index of abundance model identified areas of relatively high camel abundance. Identifying preferred habitats and areas of high abundance can help focus control efforts.

Pleopod Cuticular Morphology as an Index of Moult Stage in the Ornate Rock Lobster, Panulirus ornatus (Fabricius 1789)

Seven discrete stages and substages of moulting in the ornate rock lobster, Panulirus ornatus, have been distinguished by microscopic examination of the cuticle and setae of the pleopods . The diagnostic features and the duration of each of the stages are described. Freezing did not visually alter the tissue features used to identify each moult stage. Pleopod morphology can reliably indicate whether a lobster has moulted within the previous 24 h or is within 72 h of the next ecdysis.

A method for mapping the distribution and density of rabbits and other vertebrate pests in Australia

The European wild rabbit has been considered Australia’s worst vertebrate pest and yet little effort appears to have gone into producing maps of rabbit distribution and density. Mapping the distribution and density of pests is an important step in effective management. A map is essential for estimating the extent of damage caused and for efficiently planning and monitoring the success of pest control operations. This paper describes the use of soil type and point data to prepare a map showing the distribution and density of rabbits in Australia. The potential for the method to be used for mapping other vertebrate pests is explored. The approach used to prepare the map is based on that used for rabbits in Queensland (Berman et al. 1998). An index of rabbit density was determined using the number of Spanish rabbit fleas released per square kilometre for each Soil Map Unit (Atlas of Australian Soils). Spanish rabbit fleas were released into active rabbit warrens at 1606 sites in the early 1990s as an additional vector for myxoma virus and the locations of the releases were recorded using a Global Positioning System (GPS). Releases were predominantly in arid areas but some fleas were released in south east Queensland and the New England Tablelands of New South Wales. The map produced appears to reflect well the distribution and density of rabbits, at least in the areas where Spanish fleas were released. Rabbit pellet counts conducted in 2007 at 54 sites across an area of south east South Australia, south eastern Queensland, and parts of New South Wales (New England Tablelands and south west) in soil Map Units where Spanish fleas were released, provided a preliminary means to ground truth the map. There was a good relationship between mean pellet count score and the index of abundance for soil Map Units. Rabbit pellet counts may allow extension of the map into other parts of Australia where there were no Spanish rabbit fleas released and where there may be no other consistent information on rabbit location and density. The recent Equine Influenza outbreak provided a further test of the value of this mapping method. The distribution and density of domestic horses were mapped to provide estimates of the number of horses in various regions. These estimates were close to the actual numbers of horses subsequently determined from vaccination records and registrations. The soil Map Units are not simply soil types they contain information on landuse and vegetation and the soil classification is relatively localised. These properties make this mapping method useful, not only for rabbits, but also for other species that are not so dependent on soil type for survival.

Population dynamics of house mice in Queensland grain-growing areas

Context. Irregular plagues of house mice cause high production losses in grain crops in Australia. If plagues can be forecast through broad-scale monitoring or model-based prediction, then mice can be proactively controlled by poison baiting. Aims. To predict mouse plagues in grain crops in Queensland and assess the value of broad-scale monitoring. Methods. Regular trapping of mice at the same sites on the Darling Downs in southern Queensland has been undertaken since 1974. This provides an index of abundance over time that can be related to rainfall, crop yield, winter temperature and past mouse abundance. Other sites have been trapped over a shorter time period elsewhere on the Darling Downs and in central Queensland, allowing a comparison of mouse population dynamics and cross-validation of models predicting mouse abundance. Key results. On the regularly trapped 32-km transect on the Darling Downs, damaging mouse densities occur in 50% of years and a plague in 25% of years, with no detectable increase in mean monthly mouse abundance over the past 35 years. High mouse abundance on this transect is not consistently matched by high abundance in the broader area. Annual maximum mouse abundance in autumn–winter can be predicted (R2 = 57%) from spring mouse abundance and autumn–winter rainfall in the previous year. In central Queensland, mouse dynamics contrast with those on the Darling Downs and lack the distinct annual cycle, with peak abundance occurring in any month outside early spring.Onaverage, damaging mouse densities occur in 1 in 3 years and a plague occurs in 1 in 7 years. The dynamics of mouse populations on two transects ~70 km apart were rarely synchronous. Autumn–winter rainfall can indicate mouse abundance in some seasons (R2 = ~52%). Conclusion. Early warning of mouse plague formation in Queensland grain crops from regional models should trigger farm-based monitoring. This can be incorporated with rainfall into a simple model predicting future abundance that will determine any need for mouse control. Implications. A model-based warning of a possible mouse plague can highlight the need for local monitoring of mouse activity, which in turn could trigger poison baiting to prevent further mouse build-up.

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.

Factors affecting localized abundance and distribution of lesser mealworm in earth-floor broiler houses in subtropical Australia.

Factors that influence the localized abundance and distribution of lesser mealworm, Alphitobius diaperinus (Panzer), in litter of two compacted earth-floor broiler houses in subtropical Australia were studied using various experimental manipulations. Numbers of lesser mealworms substantially increased inside caged areas and under uncaged empty feed pans placed in open areas of the houses. These populations were found to be localized and independent of chicken-feed, manure, and high beetle populations that normally occur under existing feed pans. Substantial horizontal movement of larvae to under feed pans was recorded. Placing metal barriers around these pans significantly restricted this movement. In almost all treatments, lesser mealworms typically peaked in numbers during the middle of the flock time. This temporal pattern of abundance also was observed under pans within barriers, where relatively low insect numbers occurred, but it was not observed in uncaged open areas (where chickens had complete access). It is likely that larvae do not establish in open areas, but fluctuate in numbers as they either move to refuges away from chickens or suffer high rates of mortality. In these refuges, larvae peak in numbers and then leave the litter environment to pupate in the earth floor before the end of the flock time. This behavior might be exploited for management of lesser mealworm by targeting applications of control agents.

Female Reproductive Biology and Spawning Periodicity of Two Species of King Prawns, Penaeus longistylus Kubo and Penaeus latisulcatus Kishinouye, from Queensland's East Coast Fishery

In the coastal region of central Queensland female red-spot king prawns, P. longistylus, and the western or blue-leg king prawns, P. latisulcatus, had high mean ovary weights and high proportions of advanced ovary development during the winter months of July and August of 1985 and 1986. On the basis of insemination, both species began copulating at the size of 26-27 mm CL, but P. longistylus matured and spawned at a smaller size than P. latisulcatus. Abundance of P. longistylus was generally three to four times greater than that of P. latisulcatus but the latter was subject to greater variation in abundance. Low mean ovary weight and low proportions of females with advanced ovaries were associated with the maximum mean bottom sea-water temperature (28.5ºC) for both species. Population fecundity indices indicated that peaks in yolk or egg production (a) displayed a similar pattern for both species, (b) varied in timing from year to year for both species and (c) were strongly influenced by abundance. Generally, sample estimates of abundance and commercial catch rates (CPUE) showed similar trends. Differences between the two may have been due to changes in targeted commercial effort in this multi-species fishery.

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.

Will biological control of Lantana camara ever succeed? Patterns, processes & prospects

Despite biocontrol research spanning over 100 years, the hybrid weed, commonly referred to as Lantana camara, is not under adequate control. Host specificity and varietal preference of released agents, climatic suitability of a region for released agents, number of agents introduced and range or area of infestation appear to play a role in limiting biocontrol success. At least one of 41 species of mainly leaf- or flower-feeding insects has been introduced, or spread, to 41 of the 70 countries or regions where lantana occurs. Over half (26) of these species have established, achieving varying levels of herbivory and presumably some degree of control. Accurate taxonomy of the plant and adaptation of potential agents to the host plant are some of the better predictors of at least establishment success. Retrospective analysis of the hosts of introduced biocontrol agents for L. camara show that a greater proportion of agents that were collected from L. camara or Lantana urticifolia established, than agents that were collected from other species of Lantana. Of the introduced agents that had established and were oligophagous, 18 out of 22 established. The proportion of species establishing, declined with the number of species introduced. However, there was no trend when oceanic islands were treated separately from mainland areas and the result is likely an artefact of how introductions have changed over time. A calculated index of the degree of herbivory due to agents known to have caused some damage per country, was not related to land area infested with lantana for mainlands nor for oceanic islands. However, the degree of herbivory is much higher on islands than mainlands. This difference between island and mainland situations may reflect population dynamics in patchy or metapopulation landscapes. Basic systematic studies of the host remain crucial to successful biocontrol, especially of hybrid weeds like L. camara. Potential biocontrol agents should be monophages collected from the most closely related species to the target weed or be phytophages that attack several species of lantana. Suitable agents should be released in the most ideal ecoclimatic area. Since collection of biocontrol agents has been limited to a fraction of the known number of phytophagous species available, biocontrol may be improved by targeting insects that feed on stems and roots, as well as the agents that feed on leaves and flowers.

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.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60-103 days or 94-136 days, and those for AGECL were 95-125 or 141-176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

Benefits of oxygation of subsurface drip-irrigation water for cotton in a Vertosol

Australian cotton (Gossypium hirsutum L.) is predominantly grown on heavy clay soils (Vertosols). Cotton grown on Vertosols often experiences episodes of low oxygen concentration in the root-zone, particularly after irrigation events. In subsurface drip-irrigation (SDI), cotton receives frequent irrigation and sustained wetting fronts are developed in the rhizosphere. This can lead to poor soil diffusion of oxygen, causing temporal and spatial hypoxia. As cotton is sensitive to waterlogging, exposure to this condition can result in a significant yield penalty. Use of aerated water for drip irrigation (‘oxygation’) can ameliorate hypoxia in the wetting front and, therefore, overcome the negative effects of poor soil aeration. The efficacy of oxygation, delivered via SDI to broadacre cotton, was evaluated over seven seasons (2005–06 to 2012–13). Oxygation of irrigation water by Mazzei air-injector produced significantly (P < 0.001) higher yields (200.3 v. 182.7 g m–2) and water-use efficiencies. Averaged over seven years, the yield and gross production water-use index of oxygated cotton exceeded that of the control by 10% and 7%, respectively. The improvements in yields and water-use efficiency in response to oxygation could be ascribed to greater root development and increased light interception by the crop canopies, contributing to enhanced crop physiological performance by ameliorating exposure to hypoxia. Oxygation of SDI contributed to improvements in both yields and water-use efficiency, which may contribute to greater economic feasibility of SDI for broadacre cotton production in Vertosols.

Stock assessment of the Queensland east coast common coral trout (Plectropomus leopardus) fishery

Common coral trout Plectropomus leopardus is an iconic fish of the Great Barrier Reef (GBR) and is the most important fish for the commercial fishery there. Most of the catch is exported live to Asia. This stock assessment was undertaken in response to falls in catch sizes and catch rates in recent years, in order to gauge the status of the stock. It is the first stock assessment ever conducted of coral trout on the GBR, and brings together a multitude of different data sources for the first time. The GBR is very large and was divided into a regional structure based on the Bioregions defined by expert committees appointed by the Great Barrier Reef Marine Park Authority (GBRMPA) as part of the 2004 rezoning of the GBR. The regional structure consists of six Regions, from the Far Northern Region in the north to the Swains and Capricorn–Bunker Regions in the south. Regions also closely follow the boundaries between Bioregions. Two of the northern Regions are split into Subregions on the basis of potential changes in fishing intensity between the Subregions; there are nine Subregions altogether, which include four Regions that are not split. Bioregions are split into Subbioregions along the Subregion boundaries. Finally, each Subbioregion is split into a “blue” population which is open to fishing and a “green” population which is closed to fishing. The fishery is unusual in that catch rates as an indicator of abundance of coral trout are heavily influenced by tropical cyclones. After a major cyclone, catch rates fall for two to three years, and rebound after that. This effect is well correlated with the times of occurrence of cyclones, and usually occurs in the same month that the cyclone strikes. However, statistical analyses correlating catch rates with cyclone wind energy did not provide significantly different catch rate trends. Alternative indicators of cyclone strength may explain more of the catch rate decline, and future work should investigate this. Another feature of catch rates is the phenomenon of social learning in coral trout populations, whereby when a population of coral trout is fished, individuals quickly learn not to take bait. Then the catch rate falls sharply even when the population size is still high. The social learning may take place by fish directly observing their fellows being hooked, or perhaps heeding a chemo-sensory cue emitted by fish that are hooked. As part of the assessment, analysis of data from replenishment closures of Boult Reef in the Capricorn–Bunker Region (closed 1983–86) and Bramble Reef in the Townsville Subregion (closed 1992–95) estimated a strong social learning effect. A major data source for the stock assessment was the large collection of underwater visual survey (UVS) data collected by divers who counted the coral trout that they sighted. This allowed estimation of the density of coral trout in the different Bioregions (expressed as a number of fish per hectare). Combined with mapping data of all the 3000 or so reefs making up the GBR, the UVS results provided direct estimates of the population size in each Subbioregion. A regional population dynamic model was developed to account for the intricacies of coral trout population dynamics and catch rates. Because the statistical analysis of catch rates did not attribute much of the decline to tropical cyclones, (and thereby implied “real” declines in biomass), and because in contrast the UVS data indicate relatively stable population sizes, model outputs were unduly influenced by the unlikely hypothesis that falling catch rates are real. The alternative hypothesis that UVS data are closer to the mark and declining catch rates are an artefact of spurious (e.g., cyclone impact) effects is much more probable. Judging by the population size estimates provided by the UVS data, there is no biological problem with the status of coral trout stocks. The estimate of the total number of Plectropomus leopardus on blue zones on the GBR in the mid-1980s (the time of the major UVS series) was 5.34 million legal-sized fish, or about 8400 t exploitable biomass, with an 2 additional 3350 t in green zones (using the current zoning which was introduced on 1 July 2004). For the offshore regions favoured by commercial fishers, the figure was about 4.90 million legal-sized fish in blue zones, or about 7700 t exploitable biomass. There is, however, an economic problem, as indicated by relatively low catch rates and anecdotal information provided by commercial fishers. The costs of fishing the GBR by hook and line (the only method compatible with the GBR’s high conservation status) are high, and commercial fishers are unable to operate profitably when catch rates are depressed (e.g., from a tropical cyclone). The economic problem is compounded by the effect of social learning in coral trout, whereby catch rates fall rapidly if fishers keep returning to the same fishing locations. In response, commercial fishers tend to spread out over the GBR, including the Far Northern and Swains Regions which are far from port and incur higher travel costs. The economic problem provides some logic to a reduction in the TACC. Such a reduction during good times, such as when the fishery is rebounding after a major tropical cyclone, could provide a net benefit to the fishery, as it would provide a margin of stock safety and make the fishery more economically robust by providing higher catch rates during subsequent periods of depressed catches. During hard times when catch rates are low (e.g., shortly after a major tropical cyclone), a change to the TACC would have little effect as even a reduced TACC would not come close to being filled. Quota adjustments based on catch rates should take account of long-term trends in order to mitigate variability and cyclone effects in data.

Mapping the carriage of flaA-restriction fragment length polymorphism Campylobacter genotypes on poultry carcasses through the processing chain and comparison to clinical isolates

Poultry are considered a major source for campylobacteriosis in humans. A total of 1866 Campylobacter spp. isolates collected through the poultry processing chain were typed using flaA-restriction fragment length polymorphism to measure the impact of processing on the genotypes present. Temporally related human clinical isolates (n = 497) were also typed. Isolates were obtained from whole chicken carcass rinses of chickens collected before scalding, after scalding, before immersion chilling, after immersion chilling and after packaging as well as from individual caecal samples. A total of 32 genotypes comprising at least four isolates each were recognised. Simpson's Index of Diversity (D) was calculated for each sampling site within each flock, for each flock as a whole and for the clinical isolates. From caecal collection to after packaging samples the D value did not change in two flocks, decreased in one flock and increased in the fourth flock. Dominant genotypes occurred in each flock but their constitutive percentages changed through processing. There were 23 overlapping genotypes between clinical and chicken isolates. The diversity of Campylobacter is flock dependant and may alter through processing. This study confirms that poultry are a source of campylobacteriosis in the Australian population although other sources may contribute.