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.

Fishery biology and management of Protonibea diacanthus (Sciaenidae) aggregations in far Northern Cape York Peninsula waters

The sciaenid Protonibea diacanthus is a large, long-lived predatory fish of inshore northern Australian waters, which forms annual aggregations that are fished extensively by traditional (subsistence) and recreational fishers. There are now widespread concerns that the resource is being overexploited. Indigenous fishers of the Cape York Northern Peninsula Area (NPA) relate that large adult fish (up to 1500 mm total length (TL)) made up the bulk of the catch from the sciaenid aggregations until about 1994. In contrast, sexually mature P. diacanthus comprised only a small component (12 fish out of 270=4.4%) examined in a 1999–2000 sampling programme that was biased towards the largest individuals available. At 790 mm TL, the minimum size at first maturity for female P. diacanthus in this study is much smaller than the 920 mm TL reported previously in Queensland waters. Developing ovaries were observed in specimens sampled from sciaenid aggregations which formed in NPA waters between May and September 2000. However, no fish with ripe or spent gonads were found in the study, so the current timing and location of the spawning season for P. diacanthus in the region remain unknown. Food items observed in the analysis of the diet of P. diacanthus from the NPA included a variety of teleosts and invertebrates. The range of animal taxa represented in the prey items support the description of an ‘opportunistic predator’ attributed to the species. In our sampling, the stomach contents of fish caught during the time of the aggregation events did not differ from those observed at other times of the year. A total of 114 P. diacanthus were tagged and released at aggregation sites during the study period, and 3 fish (2.6%) were subsequently recaptured. The low rate of tag returns from the wild stock tagging programme, both in this study (2.6%) and from recreational fisher tag/release programmes for the sciaenid elsewhere in Queensland (6.5%), were not explained by tag loss nor mortality, given the high retention rate of tags and the zero mortality seen in tank trials. In response to the biological findings from this study, indigenous community councils of the NPA imposed a 2-year fishing moratorium for P. diacanthus. Surveys at aggregation sites in 2002 and 2003 established that much larger fish (mean size 103.5 cm TL) were again present on the grounds, albeit in very low numbers. These recent preliminary results highlight the critical need for continued monitoring and management of the P. diacanthus fishery in the NPA, if prospects for resource recovery are to be realised. The NPA initiative has provided a rare opportunity to negotiate a co-management strategy, based on scientific data and traditional knowledge, for the recovery of a cultural and economically significant fished resource.

Queensland Coastal Wetland Resource Investigation of the Bowen Region: Cape Upstart to Gloucester Island

The wetland resources of the Queensland coastline have been mapped as a baseline dataset for Marine Protected Area investigation and particularly Fish Habitat Area (FHA) declaration, Ramsar site nomination and continued monitoring of these important fish habitats. This report summarises the results of the mapping undertaken in the Bowen region from the East Coast of Cape Upstart (Abbot Bay) to Gloucester Island (encompassing Edgecumbe Bay). The study was undertaken in order to: 1. document and map the coastal wetland communities within the Bowen region; 2. document levels of existing disturbance to and protection of the wetlands; 3. examine existing recreational and commercial fisheries in the region; and 4. evaluate the significance of the coastal wetlands in the region. Dataset URL Link: Queensland Coastal Wetlands Resources Mapping data. [Dataset]

Resource Assessment of the Tidal Wetland Vegetation of Western Cape York Peninsula, North Queensland, Report to Ocean Rescue 2000

A project to allow the resource assessment of tidal wetland vegetation of western Cape York Peninsula, in north Queensland, was undertaken as part of the longterm assessment of the coastal fisheries resources of Queensland. The project incorporated a littoral invertebrate fauna component. Extending from May 1993 to December 1994, fieldwork was undertaken in May 1993, November 1993 and April 1994. The aims of this project were to: • obtain baseline information on the distribution of marine plants of western Cape York Peninsula; • commence a preliminary assessment of the littoral invertebrate fauna and their habitat requirements with a view to extending knowledge of their biogeographic affinities; • perform biogeographic classification of the tidal wetlands at a meso and local scale for marine conservation planning; • evaluate the conservation values of the areas investigated from the viewpoint of fisheries productivity and as habitat for important/threatened species. Dataset URL Link: Queensland Coastal Wetlands Resources Mapping data. [Dataset]

An investigation of migration and possible stock structuring by the stout whiting, Sillago robusta, in southern Queensland waters, and its impact on managing the fishery. Final Report on the Pilot Program

A highly polymorphic genetic locus of Stout Whiting was examined for evidence of geographical subdivision amongst samples collected from three locales in southern Queensland waters. Statistical indicators of subdivision were not significantly different from zero, suggesting that it is unlikely that the Stout Whiting resource in southern Queensland is genetically subdivided into separate stocks. It is recommended that the full-scale genetic program not proceed and that the resource be managed as a single stock.

Development of a Gulf community based natural resource monitoring program.

Development of a Gulf community based natural resource monitoring program, with sawfish as an initial focus.

Targeting resource investments to achieve sediment reduction and improved Great Barrier Reef health

Concerns about excessive sediment loads entering the Great Barrier Reef (GBR) lagoon in Australia have led to a focus on improving ground cover in grazing lands. Ground cover has been identified as an important factor in reducing sediment loads, but improving ground cover has been difficult for reef stakeholders in major catchments of the GBR. To provide better information an optimising linear programming model based on paddock scale information in conjunction with land type mapping was developed for the Fitzroy, the largest of the GBR catchments. This identifies at a catchment scale which land types allow the most sediment reduction to be achieved at least cost. The results suggest that from the five land types modelled, the lower productivity land types present the cheapest option for sediment reductions. The study allows more informed decision making for natural resource management organisations to target investments. The analysis highlights the importance of efficient allocation of natural resource management funds in achieving sediment reductions through targeted land type investments. © 2012.

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.

Reproductive Biology and Spawning Periodicity of Endeavour Shrimps Metapenaeus endeavouri (Schmitt, 1926) and Metapenaeus ensis (de Haan, 1850) from a Central Queensland (Australia) Fishery

Metapenaeus endeavouri and M. ensis from coastal trawl fishing grounds off central Queensland, Australia, have marked seasonal reproductive cycles. Female M. endeavouri grew to a larger size than female M. ensis and occurred over a wider range of sites and depths. Although M. ensis was geographically restricted in distribution to only the shallowest sites it was highly abundant. Mating activity in these open thelycum species, indicated by the presence or absence of a spermatophore, was relatively low and highly seasonal compared with closed thelyeum shrimps. Seasonal variation in spermatophore insemination can be used as an independent technique to study spawning periodicity in open thelycum shrimps. Data strongly suggest an inshore movement of M. endeavouri to mature and spawn. This differs from most concepts of Penaeus species life cycles, but is consistent with the estuarine significance in the life cycle of Metapenaeus species. Monthly population fecundity indices suggest summer spawning for both species, which contrasts with the winter spawning of other shrimps from the same multispecies fishery.

The fishery and reproductive biology of Barking crayfish, Linuparus trigonus (Von Siebold, 1824) along Queenslands east coast

This paper describes the fishery and reproductive biology for Linuparus trigonus obtained from trawl fishermen operating off Queensland’s east coast, Australia. The smallest mature female lobster measured 59.8 mm CL, however, 50% maturity was reached between 80 and 85 mm CL. Brood fecundity (BF) was size dependent and ranged between 19,287 and 100,671 eggs in 32 females from 59.8 to 104.3 mm CL. The relationship was best described by the power equation BF = 0.1107*CL to the power of 2.9241 (r to the power of 2 = 0:74). Egg size ranged from 0.96 to 1.12 mm in diameter (mean = 1:02 (+or-) 0:01 mm). Egg weight and size were independent of lobster size. Length frequencies displayed multi-modal distributions.The percentage of female to male lobsters was relatively stable for small size classes (30 to 70 mm CL; 50.0 to 63.6% females), but female proportions rose markedly between 75 and 90 mm (72.2 to 85.4%) suggesting that at the onset of sexual maturity female growth rates are reduced. In size classes greater than 95 mm, males were numerically dominant. A description of the L. trigonus fishery in Queensland is also detailed.

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

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

Quantifying the effects of bycatch reduction devices in Queensland's (Australia) shallow water eastern king prawn (Penaeus plebejus) trawl fishery

This study presents results from an experimental 10-day research charter that was designed to quantify the effects of a) a turtle excluder device (TED) and b) a radial escape section bycatch reduction device (BRD) and c) both devices together, on prawn and bycatch catch rates in the Queensland shallow water eastern king prawn (Penaeus plebejus) trawl fishery. The bycatch was comprised of 250 taxa, mainly gurnards, whiting, lizard fish, flathead, dragonets, portunid crabs, turretfish and flounders. The observed mean catch rates of bycatch and marketable eastern king prawns from the standard trawl net (i.e., net with no TED or BRD) used during the charter were 11.06 (se 0.90) kg per hectare swept by the trawl gear (ha-1) and 0.94 kg ha-1, respectively. For the range of depths sampled (20.1-90.7 m), bycatch catch rates declined significantly at a rate of 0.14 kg ha-1 for every 1 m increase in depth, while prawn catch rates were unaffected. When both the TED and radial escape section BRD were used together they resulted in a 24% reduction in total bycatch catch rate compared to a standard net, but at a 20% reduction in marketable prawn catch rate. The largest reductions were achieved for stout whiting Sillago robusta (57% reduction) and yellowtail scad Trachurus novaezelandiae (32% reduction). Multidimensional scaling and analysis of similarities revealed that bycatch assemblages differed significantly between depths and latitude, but not between the different combinations of bycatch reduction devices. Despite the lowered prawn catch rates, the reduced bycatch catch rates are promising, particularly for S. robusta which is not permitted to be retained by the prawn trawl fleet and yet experiences considerable incidental fishing mortality, and because it is targeted in a separate licensed commercial fishery.

Seasonal resource availability and use by an endangered tropical mycophagous marsupial

This study highlights the importance of considering how seasonality of rainfall affects availability of resources and consequently species distributions within tropical ecosystems. The endangered northern bettong, Bettongia tropica Wakefield is thought to be restricted to habitats where seasonal availability of hypogeous fungi, their principal food resource, remains high. To test this hypothesis fungal abundance was quantified in the early wet, late wet, early dry and late dry seasons within known bettong habitat. A relationship was found between precipitation and fungal availability, with the abundance of hypogeous fungi being significantly lower in the late dry season. Fungal availability correlated strongly with the seasonal rainfall pattern determined from 74-year monthly means. This contrasts with a previous study where mycophagy, measured by faecal analysis, remained high across seasons presumably because of aseasonal rainfall during that study period. Alloteropsis semialata R.Br. (cockatoo grass) use by bettongs increased significantly during the period of low fungal availability. This suggests that the importance of cockatoo grass as an alternative food resource during annual and extended dry periods has previously been underestimated. With the frequency and intensity of drought expected to increase with global climate change, these findings have significant implications for bettong management. The important and possibly equivalent dependence of B. tropica on both hypogeous fungi and A. semialata helps to explain their habitat preference and identifies this species as a true ecotonal specialist.

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

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

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

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