Linking spatial stock dynamics and economics: Evaluation of indicators and fishery management for the travelling eastern king prawn (Melicertus plebejus)

Reduced economic circumstances havemoved management goals towards higher profit, rather than maximum sustainable yields in several Australian fisheries. The eastern king prawn is one such fishery, for which we have developed new methodology for stock dynamics, calculation of model-based and data-based reference points and management strategy evaluation. The fishery is notable for the northward movement of prawns in eastern Australian waters, from the State jurisdiction of New South Wales to that of Queensland, as they grow to spawning size, so that vessels fishing in the northern deeper waters harvest more large prawns. Bioeconomic fishing data were standardized for calibrating a length-structured spatial operating model. Model simulations identified that reduced boat numbers and fishing effort could improve profitability while retaining viable fishing in each jurisdiction. Simulations also identified catch rate levels that were effective for monitoring in simple within-year effort-control rules. However, favourable performance of catch rate indicators was achieved only when a meaningful upper limit was placed on total allowed fishing effort. Themethods and findings will allow improved measures for monitoring fisheries and inform decision makers on the uncertainty and assumptions affecting economic indicators.

Latitudinal and seasonal effects on growth of the Australian eastern king prawn (Melicertus plebejus)

The growth of the Australian eastern king prawn (Melicertus plebejus) is understood in greater detail by quantifying the latitudinal effect. The latitudinal effect is the change in the species' growth rate during migration. Mark-recapture data (N = 1635, latitude 22.21 degrees S-34.00 degrees S) presents northerly movement of the eastern king prawn, with New South Wales prawns showing substantial average movement of 140 km (standard deviation: 176 km) north. A generalized von Bertalanffy growth model framework is used to incorporate the latitudinal effect together with the canonical seasonal effect. Applying this method to eastern king prawn mark-recapture data guarantees consistent estimates for the latitudinal and seasonal effects. For M. plebejus, it was found that growth rate peaks on 25 and 29 January for males and females, respectively; is at a minimum on 27 and 31 July, respectively; and that the shape parameter, k (per year), changes by -0.0236 and -0.0556 every 1 degree of latitude south increase for males and females, respectively.

Growth curves with time-dependent explanatory variables

A class of growth models incorporating time-dependent factors and stochastic perturbations are introduced. The proposed model includes the existing growth models used in fisheries as special cases. Particular attention is given to growth of a population (in average weight or length) from which observations are taken randomly each time and the analysis of tag-recapture data. Two real data sets are used for illustration: (a) to estimate the seasonal effect and population density effect on growth of farmed prawn (Penaeus monodon) from weight data and (b) to assess the effect of tagging on growth of barramundi (Lates calcarifer)

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

A maximum-likelihood method for estimating natural mortality and catchability coefficient from catch-and-effort data

A simple stochastic model of a fish population subject to natural and fishing mortalities is described. The fishing effort is assumed to vary over different periods but to be constant within each period. A maximum-likelihood approach is developed for estimating natural mortality (M) and the catchability coefficient (q) simultaneously from catch-and-effort data. If there is not enough contrast in the data to provide reliable estimates of both M and q, as is often the case in practice, the method can be used to obtain the best possible values of q for a range of possible values of M. These techniques are illustrated with tiger prawn (Penaeus semisulcatus) data from the Northern Prawn Fishery of Australia.

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.

Temporal and Spatial Patterns in Recruitment of Three Penaeid Prawns in Moreton Bay, Queensland, Australia

The size at recruitment, temporal and spatial distribution, and abiotic factors influencing abundance of three commercially important species of penaeid prawns in the sublittoral trawl grounds of Moreton Bay (Queensland, Australia) were compared. Metapenaeus bennettae and Penaeus plebejus recruit to the trawl grounds at sizes which are relatively small (14-15 mm carapace length, CL) and below that at which prawns are selected for, and retained, in the fleet's cod-ends. In contrast, Penaeus esculenlus recruit at the relatively large size of 27 mm CL from February to May, well above the size ranges selected for. Recruitment of M. bennettae extends over several months, September-October and February March, and was thus likely to be bi-annual, while the recruitment period of P. plebejus was distinct, peaking in October-November each year. Size classes of M . bennettae were the most spatially stratified of the three species. Catch rates of recruits were negatively correlated with depth for all three species, and were also negatively correlated with salinity for M. bennettae.

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.

Isolation and expression analysis of multiple isoforms of putative farnesoic acid O-methyltransferase in several crustacean species

Farnesoic acid O-methyltransferase (FaMeT) is the enzyme responsible for the conversion of farnesoic acid (FA) to methyl farnesoate (MF) in the final step of MF synthesis. Multiple isoforms of putative FaMeT were isolated from six crustacean species belonging to the families Portunidae, Penaeidae, Scyllaridae and Parastacidae. The portunid crabs Portunus pelagicus and Scylla serrata code for three forms: short, intermediate and long. Two isoforms (short and long) were isolated from the penaeid prawns Penaeus monodon and Fenneropenaeus merguiensis. Two isoforms were also identified in the scyllarid Thenus orientalis and parastacid Cherax quadricarinatus. Putative FaMeT sequences were also amplified from the genomic DNA of P. pelagicus and compared to the putative FaMeT transcripts expressed. Each putative FaMeT cDNA isoform was represented in the genomic DNA, indicative of a multi-gene family. Various tissues from P. pelagicus were individually screened for putative FaMeT expression using PCR and fragment analysis. Each tissue type expressed all three isoforms of putative FaMeT irrespective of sex or moult stage. Protein domain analysis revealed the presence of a deduced casein kinase II phosphorylation site present only in the long isoform of putative FaMeT.

New Modified Dipping Method Using 4-hexylresorcinol for Preventing Blackspot Formation in Prawns.

Industry acceptance of alternative sulphite treatments for preventing blackspot in prawns is limited. This study demonstrated that dipping prawns (Penaeus monodon and Fenneropenaeus merguiensis) in a less concentrated solution for longer times provides better control of blackspot formation than the recommended 2 minute dip in 50 mg/L 4-hexylresorcinol. Protection improved as modified dipping times and solution concentrations increased. Blackspot protection provided by most of the modified dip methods was more effective than a 1 minute dip in 1% sodium metabisulphite. The recommended 2 minute dip method should be limited to 125 kg when used on consecutive 25-kg batches of prawns. Yields increased by less than 10% even after 72-hour exposure in the modified dip. The new dipping method meets most international residue standards when applying the modified treatments evaluated.

Effects of secondary crops on bacterial growth and nitrogen removal in shrimp farm effluent treatment systems

Secondary crops provide a means of assimilating some effluent nitrogen from eutrophic shrimp farm settlement ponds. However, a more important role may be their stimulation of beneficial bacterial nitrogen removal processes. In this study, bacterial biomass, growth and nitrogen removal capacity were quantified in shrimp farm effluent treatment systems containing vertical artificial substrates and either the banana shrimp Penaeus merguiensis (de Man) or the grey mullet, Mugil cephalus L. Banana shrimp were found to actively graze biofilm on the artificial substrates and significantly reduced bacterial biomass relative to a control (24.5 ± 5.6mgCm−2 and 39.2 ± 8.7mgCm−2, respectively). Bacterial volumetric growth rates, however, were significantly increased in the presence of the shrimp relative to the control 45.2±11.3mgCm−2 per day and 22.0±4.3mgCm−2 per day, respectively). Specific growth rate, or growth rate per cell, of bacteria was therefore appreciably stimulated by the banana shrimp. Nitrate assimilation was found to be significantly higher on grazed substrate biofilm relative to the control (223±54 mgNm−2 per day and 126±36 mg Nm−2 per day, respectively), suggesting that increased bacterial growth rate does relate to enhanced nitrogen uptake. Regulated banana shrimp feeding activity therefore can increase the rate of newbacterial biomass production and also the capacity for bacterial effluent nitrogen assimilation. Mullet had a negligible influence on the biofilm associated with the artificial substrate but reduced sediment bacterial biomass (224 ± 92 mgCm−2) relative to undisturbed sediment (650 ± 254 mgCm−2). Net, or volumetric bacterial growth in the sediment was similar in treatments with and without mullet, suggesting that the growth rate per cell of bacteria in grazed sediments was enhanced. Similar rates of dissolved nitrogen mineralisation werefound in sediments with and without mullet but nitrificationwas reduced. Presence of mullet increased water column suspended solids concentrations, water column bacterial growth and dissolved nutrient uptake. This study has shown that secondary crops, particularly banana shrimp, can play a stimulatory role in the bacterial processing of effluent nitrogen in eutrophic shrimp effluent treatment systems.