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.

Stock structure of Grey Mackerel, Scomberomorus semifasciatus (Pisces: Scombridae) across northern Australia, based on otolith stable isotope chemistry

The stable isotopes of delta O-18 and delta C-13 in sagittal otolith carbonates were used to determine the stock structure of Grey Mackerel, Scomberomorus semifasciatus. Otoliths were collected from Grey Mackerel at ten locations representing much of their distributional and fisheries range across northern Australia from 2005 to 2007. Across this broad range (similar to 6500 km), fish from four broad locations-Western Australia (S1), Northern Territory and Gulf of Carpentaria (S2, S3, S4, S5, S6, S7), Queensland east coast mid and north sites (S8, S9) and Queensland east coast south site (S10)-had stable isotope values that were significantly different indicating stock separation. Otolith stable isotopes differed more between locations than among years within a location, indicating temporal stability across years. The spatial separation of these populations indicates a complex stock structure across northern Australia. Stocks of S. semifasciatus appear to be associated with large coastal embayments. These results indicate that optimal fisheries management may require a review of the current spatial arrangements, particularly in relation to the evidence of shared stocks in the Gulf of Carpentaria. Furthermore, as the population of S. semifasciatus in Western Australia exhibited high spatial separation from those at all the other locations examined, further research activities should focus on investigating additional locations within Western Australia for an enhanced determination of stock delineation. From the issue entitled "Proceedings of the 4th International Otolith Symposium, 24-28 August 2009, Monterey, California"

Evaluation of aerial photography for predicting trends in structural attributes of Australian woodland including comparison with ground-based monitoring data

The accurate assessment of trends in the woody structure of savannas has important implications for greenhouse accounting and land-use industries such as pastoralism. Two recent assessments of live woody biomass change from north-east Australian eucalypt woodland between the 1980s and 1990s present divergent results. The first estimate is derived from a network of permanent monitoring plots and the second from woody cover assessments from aerial photography. The differences between the studies are reviewed and include sample density, spatial scale and design. Further analyses targeting potential biases in the indirect aerial photography technique are conducted including a comparison of basal area estimates derived from 28 permanent monitoring sites with basal area estimates derived by the aerial photography technique. It is concluded that the effect of photo-scale; or the failure to include appropriate back-transformation of biomass estimates in the aerial photography study are not likely to have contributed significantly to the discrepancy. However, temporal changes in the structure of woodlands, for example, woodlands maturing from many smaller trees to fewer larger trees or seasonal changes, which affect the relationship between cover and basal area could impact on the detection of trends using the aerial photography technique. It is also possible that issues concerning photo-quality may bias assessments through time, and that the limited sample of the permanent monitoring network may inadequately represent change at regional scales

Soil properties and crop yields in a dryland Vertisol sown with cotton-based crop rotations

Information on the effects of growing cotton (Gossypium hirsutum L.)-based crop rotations on soil quality of dryland Vertisols is sparse. The objective of this study was to quantify the effects of growing cereal and leguminous crops in rotation with dryland cotton on physical and chemical properties of a grey Vertisol near Warra, SE Queensland, Australia. The experimental treatments, selected after consultations with local cotton growers, were continuous cotton (T1), cotton-sorghum (Sorghum bicolor (L.) Moench.) (T2), cotton-wheat (Triticum aestivum L.) double cropped (T3), cotton-chickpea (Cicer arietinum L.) double cropped followed by wheat (T4) and cotton-wheat (T5). From 1993 to 1996 land preparation was by chisel ploughing to about 0.2 m followed by two to four cultivations with a Gyral tyne cultivator. Thereafter all crops were sown with zero tillage except for cultivation with a chisel plough to about 0.07-0.1 m after cotton picking to control heliothis moth pupae. Soil was sampled from 1996 to 2004 and physical (air-filled porosity of oven-dried soil, an indicator of soil compaction; plastic limit; linear shrinkage; dispersion index) and chemical (pH in 0.01 M CaCl2, organic carbon, exchangeable Ca, Mg, K and Na contents) properties measured. Crop rotation affected soil properties only with respect to exchangeable Na content and air-filled porosity. In the surface 0.15 m during 2000 and 2001 lowest air-filled porosity occurred with T1 (average of 34.6 m3/100 m3) and the highest with T3 (average of 38.9 m3/100 m3). Air-filled porosity decreased in the same depth between 1997 and 1998 from 45.0 to 36.1 m3/100 m3, presumably due to smearing and compaction caused by shallow cultivation in wet soil. In the subsoil, T1 and T2 frequently had lower air-filled porosity values in comparison with T3, T4 and T5, particularly during the early stages of the experiment, although values under T1 increased subsequently. In general, compaction was less under rotations which included a wheat crop (T3, T4, T5). For example, average air-filled porosity (in m3/100 m3) in the 0.15-0.30 m depth from 1996 to 1999 was 19.8 with both T1 and T2, and 21.2 with T3, 21.1 with T4 and 21.5 with T5. From 2000 to 2004, average air-filled porosity (in m3/100 m3) in the same depth was 21.3 with T1, 19.0 with T2, 19.8 with T3, 20.0 with T4 and 20.5 with T5. The rotation which included chickpea (T4) resulted in the lowest exchangeable Na content, although differences among rotations were small. Where only a cereal crop with a fibrous root system was sown in rotation with cotton (T2, T3, T5) linear shrinkage in the 0.45-0.60 m depth was lower than in rotations, which included tap-rooted crops such as chickpea (T4) or continuous cotton (T1). Dispersion index and organic carbon decreased, and plastic limit increased with time. Soil organic carbon stocks decreased at a rate of 1.2 Mg/ha/year. Lowest average cotton lint yield occurred with T2 (0.54 Mg/ha) and highest wheat yield with T3 (2.8 Mg/ha). Rotations which include a wheat crop are more likely to result in better soil structure and cotton lint yield than cotton-sorghum or continuous cotton.

Population structure of South African and Australian Pyrenophora teres isolates.

There are two recognized forms of the disease net blotch of barley: the net form caused by Pyrenophora teres f. teres (PTT) and the spot form caused by P. teres f. maculata (PTM). In this study, amplified fragment length polymorphism analysis was used to investigate the genetic diversity and population structure of 60 PTT and 64 PTM isolates collected across Australia (66 isolates) and in the south-western Cape of South Africa (58 isolates). For comparison, P. tritici-repentis, Exserohilum rostratum and Bipolaris sorokiniana samples were also included in the analyses. Both distance-and model-based cluster analyses separated the PTT and PTM isolates into two strongly divergent genetic groups. Significant variation was observed both among the South African and Australian populations of PTT and PTM and among sampling locations for the PTT samples. Results suggest that sexual reproduction between the two forms is unlikely and that reproduction within the PTT and PTM groups occurs mainly asexually.

Net radiation dynamics: performance of 20 daily net radiation models as related to model structure and intricacy in two climates.

We compared daily net radiation (Rn) estimates from 19 methods with the ASCE-EWRI Rn estimates in two climates: Clay Center, Nebraska (sub-humid) and Davis, California (semi-arid) for the calendar year. The performances of all 20 methods, including the ASCE-EWRI Rn method, were then evaluated against Rn data measured over a non-stressed maize canopy during two growing seasons in 2005 and 2006 at Clay Center. Methods differ in terms of inputs, structure, and equation intricacy. Most methods differ in estimating the cloudiness factor, emissivity (e), and calculating net longwave radiation (Rnl). All methods use albedo (a) of 0.23 for a reference grass/alfalfa surface. When comparing the performance of all 20 Rn methods with measured Rn, we hypothesized that the a values for grass/alfalfa and non-stressed maize canopy were similar enough to only cause minor differences in Rn and grass- and alfalfa-reference evapotranspiration (ETo and ETr) estimates. The measured seasonal average a for the maize canopy was 0.19 in both years. Using a = 0.19 instead of a = 0.23 resulted in 6% overestimation of Rn. Using a = 0.19 instead of a = 0.23 for ETo and ETr estimations, the 6% difference in Rn translated to only 4% and 3% differences in ETo and ETr, respectively, supporting the validity of our hypothesis. Most methods had good correlations with the ASCE-EWRI Rn (r2 > 0.95). The root mean square difference (RMSD) was less than 2 MJ m-2 d-1 between 12 methods and the ASCE-EWRI Rn at Clay Center and between 14 methods and the ASCE-EWRI Rn at Davis. The performance of some methods showed variations between the two climates. In general, r2 values were higher for the semi-arid climate than for the sub-humid climate. Methods that use dynamic e as a function of mean air temperature performed better in both climates than those that calculate e using actual vapor pressure. The ASCE-EWRI-estimated Rn values had one of the best agreements with the measured Rn (r2 = 0.93, RMSD = 1.44 MJ m-2 d-1), and estimates were within 7% of the measured Rn. The Rn estimates from six methods, including the ASCE-EWRI, were not significantly different from measured Rn. Most methods underestimated measured Rn by 6% to 23%. Some of the differences between measured and estimated Rn were attributed to the poor estimation of Rnl. We conducted sensitivity analyses to evaluate the effect of Rnl on Rn, ETo, and ETr. The Rnl effect on Rn was linear and strong, but its effect on ETo and ETr was subsidiary. Results suggest that the Rn data measured over green vegetation (e.g., irrigated maize canopy) can be an alternative Rn data source for ET estimations when measured Rn data over the reference surface are not available. In the absence of measured Rn, another alternative would be using one of the Rn models that we analyzed when all the input variables are not available to solve the ASCE-EWRI Rn equation. Our results can be used to provide practical information on which method to select based on data availability for reliable estimates of daily Rn in climates similar to Clay Center and Davis.

Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (FST = 0· 09 − 0· 49; P < 0· 05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.

Negligible evidence for regional genetic population structure for two shark species Rhizoprionodon acutus (Rüppell, 1837) and Sphyrna lewini (Griffith & Smith, 1834) with contrasting biology

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n = 237) and the milk shark (Rhizoprionodon acutus, n = 207) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751–0.903, respectively; microsatellite loci, 0.038–0.047 respectively). Our results support the spatially homogeneous monitoring and management plan for shark species in Queensland, Australia.

Stock structure of exploited shark species in north-eastern Australia

The project has provided management and other stakeholders with information necessary to make informed decisions about the management of four of the key exploited shark species caught in the Queensland inshore net fishery and northern New South Wales line fishery. The project has determined that spatial management of milk sharks within Queensland, and scalloped hammerhead, common black tip and Australian black tip sharks within Queensland and New South Wales is appropriate. The project has determined that both black tip shark species are likely to require co-operative management arrangements between Queensland and New South Wales. For scalloped hammerheads separate stocks between the two jurisdictions were identified from the fisheriesdependent samples, however genetic exchange across borders is likely to be facilitated by movement of adult females and perhaps larger males to a lesser extent. This information will greatly assist compliance with the Commonwealth Environment Protection and Biodiversity Conservation Act (1999) for shark fisheries in north-eastern Australia by providing the necessary basis for robust assessment of the status of stocks of the study species, thereby helping to deliver their sustainable harvest. It also helps to achieve objectives of the Australian National Shark Plan. The project provides the appropriate spatial framework for future monitoring and assessment of the study species. This is at a time when shark fisheries are receiving close attention from all sectors and when monitoring programs are being implemented, aimed at better assessment of stock status. This project has provided the crucial information for developing an appropriate monitoring design as well as the necessary basis for making statements about stock status. The project has addressed research priorities identified by the Queensland Fisheries Research Advisory Board, Great Barrier Reef Marine Park Authority and Queensland Fisheries. Previously management has assumed a single stock for each species on the east coast of Queensland, and management of shark fisheries in New South Wales (NSW) and Queensland has been independent of one another. The project has been able to enhance and develop links between research, management and industry. Strong positive relationships with commercial fishers were crucial in the collection of samples throughout the study area and fisheries managers were part of the project team throughout the study period. During the project the study area was extended to include both Queensland and NSW waters, creating mutualistic and positive links between the States’ research and management agencies. Extension of project results included management representatives from NSW and Queensland, as well as the Northern Territory where similar shark fisheries operate and similar species are targeted. The project was able to provide significant human capital development opportunities providing considerable value to the project outcomes. Use of vertebral microchemistry and life history characteristics as stock determination methods provided material for two PhD students based at James Cook University: Ron Schroeder, vertebral chemistry; and Alastair Harry, life history characteristic. The project has developed novel research methods that have great capacity for future application, including: • Development of a simple and rapid genetic diagnostic tool (RT-HRM-PCR assay) for differentiating among the black tip shark species, for which no simple morphological identifier exists; and • Development of laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) methods for analysing and interpreting microchemical composition of shark vertebrae. The study has provided further confirmation of the effectiveness of using a holistic approach in stock structure studies and justifies investment into such studies.

CattleTickBase: An integrated Internet-based bioinformatics resource for Rhipicephalus (Boophilus) microplus

The Rhipicephalus micro plus genome is large and complex in structure, making it difficult to assemble a genome sequence and costly to resource the required bioinformatics. In light of this, a consortium of international collaborators was formed to pool resources to begin sequencing this genome. We have acquired and assembled genomic DNA into contigs that represent over 1.8 Gigabase pairs of DNA from gene-enriched regions of the R. micro plus genome. We also have several datasets containing transcript sequences from a number of gene expression experiments conducted by the consortium. A web-based resource was developed to enable the scientific community to access our datasets and conduct analysis through a web-based bioinformatics environment called YABI. The collective bioinformatics resource is termed CattleTickBase. Our consortium has acquired genomic and transcriptomic sequence data at approximately 0.9X coverage of the gene-coding regions of the R. microplus genome. The YABI tool will facilitate access and manipulation of cattle tick genome sequence data as the genome sequencing of R. microplus proceeds. During this process the CattleTickBase resource will continue to be updated. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology Inc.

Extremely Low Microsatellite Diversity but Distinct Population Structure in a Long-Lived Threatened Species, the Australian Lungfish Neoceratodus forsteri (Dipnoi)

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F_ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.