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.

A retrospective evaluation of sustainable yields for Australia's Northern Prawn Fishery

The fundamental aim in fisheries management is to determine an optimal fishing effort for sustainably harvesting from a replenishable resource. The current management objective of Australia's Northern Prawn Fishery is to maximize the long-term net economic return following Australian government policy, resulting in an average recent catch of tiger prawn species of about 1,250 tons only. However, the maximum sustainable catch stated from different studies is around 3,000-4,700 tons. We also evaluated the net profit assuming that there was no buyback scheme in 2005 and the fishing fleet was kept at 89 vessels since 2005 and concluded that 40% more catch on average (2006-2009) and an additional total profit of A$ 17 million ( excluding crew cost) could have been gained in addition to the many millions of dollars of savings in the buyback scheme. These findings have great implications for future management in Australia and elsewhere because there is a grave concern of overfishing worldwide.

Analysing commercial catch and effort data from a penaeid trawl fishery: A comparison of linear models, mixed models, and generalised estimating equations approaches

Statistical methods are often used to analyse commercial catch and effort data to provide standardised fishing effort and/or a relative index of fish abundance for input into stock assessment models. Achieving reliable results has proved difficult in Australia's Northern Prawn Fishery (NPF), due to a combination of such factors as the biological characteristics of the animals, some aspects of the fleet dynamics, and the changes in fishing technology. For this set of data, we compared four modelling approaches (linear models, mixed models, generalised estimating equations, and generalised linear models) with respect to the outcomes of the standardised fishing effort or the relative index of abundance. We also varied the number and form of vessel covariates in the models. Within a subset of data from this fishery, modelling correlation structures did not alter the conclusions from simpler statistical models. The random-effects models also yielded similar results. This is because the estimators are all consistent even if the correlation structure is mis-specified, and the data set is very large. However, the standard errors from different models differed, suggesting that different methods have different statistical efficiency. We suggest that there is value in modelling the variance function and the correlation structure, to make valid and efficient statistical inferences and gain insight into the data. We found that fishing power was separable from the indices of prawn abundance only when we offset the impact of vessel characteristics at assumed values from external sources. This may be due to the large degree of confounding within the data, and the extreme temporal changes in certain aspects of individual vessels, the fleet and the fleet dynamics.

Subsampling multi-species trawl catches from tropical northern Australia: Does it matter which part of the catch is sampled?

Subsampling is a common method for estimating the abundance of species in trawl catches. However, the accuracy of subsampling in representing the total catch has not been assessed. To estimate one possible source of bias due to subsampling, we tested whether the position on trawler sorting trays from which subsamples were taken affected their ability to represent species in catches. This was done by sorting catches into 10 kg subsamples and comparing subsamples taken from different positions on the sorting tray. Comparisons were made after species were grouped into three categories of abundance, either 'rare', 'common' or 'abundant'. A generalised linear model analysis showed that taking subsamples from different positions on the sorting tray had no major effect on estimating the total numbers or weights of fish or invertebrates, or the total number of fish or invertebrate taxa, recorded in each position. Some individual taxa showed differences between positions on the sorting tray (11.5% of taxa ina three-position design; 25% in a five-position design). But consistent and meaningful patterns in the position of these taxa on the sorting tray could only be seen for the pony fish Leiognathus moretoniensis and the saucer scallop Amusium pleuronectes. Because most bycatch laxa are well mixed throughout the catch, subsamples can be taken from any position on trawler sorting trays without introducing bias.

A generalized estimating equations approach for analysis of the impact of new technology on a trawl fishery

The article describes a generalized estimating equations approach that was used to investigate the impact of technology on vessel performance in a trawl fishery during 1988-96, while accounting for spatial and temporal correlations in the catch-effort data. Robust estimation of parameters in the presence of several levels of clustering depended more on the choice of cluster definition than on the choice of correlation structure within the cluster. Models with smaller cluster sizes produced stable results, while models with larger cluster sizes, that may have had complex within-cluster correlation structures and that had within-cluster covariates, produced estimates sensitive to the correlation structure. The preferred model arising from this dataset assumed that catches from a vessel were correlated in the same years and the same areas, but independent in different years and areas. The model that assumed catches from a vessel were correlated in all years and areas, equivalent to a random effects term for vessel, produced spurious results. This was an unexpected finding that highlighted the need to adopt a systematic strategy for modelling. The article proposes a modelling strategy of selecting the best cluster definition first, and the working correlation structure (within clusters) second. The article discusses the selection and interpretation of the model in the light of background knowledge of the data and utility of the model, and the potential for this modelling approach to apply in similar statistical situations.

The impact of global positioning systems and plotters on fishing power in the northern prawn fishery, Australia

The impact of global positioning systems (GPS) and plotter systems on the relative fishing power of the northern prawn fishery fleet on tiger prawns (Penaeus esculentus Haswell, 1879, and P. semisulcatus de Haan, 1850) was investigated from commercial catch data. A generalized linear model was used to account for differences in fishing power between boats and changes in prawn abundance. It was found that boats that used a GPS alone had 4% greater fishing power than boats without a CPS. The addition of a plotter raised the power by 7% over boats without the equipment. For each year between the first to third that a fisher has been working with plotters, there is an additional 2 or 3% increase. It appears that when all boats have a GPS and plotter for at least 3 years, the fishing power of the fleet will increase by 12%. Management controls have reduced the efficiency of each boat and lowered the number of days available to fish, but this may not have been sufficient to counteract the increases. Further limits will be needed to maintain the desired levels of mortality.

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.

Stock-recruitment relationships of the tiger prawns (Penaeus esculentus and Penaeus semisulcatus) in the Australian Northern Prawn Fishery

This paper investigates the stock-recruitment and equilibrium yield dynamics for the two species of tiger prawns (Penaeus esculentus and Penaeus semisulcatus) in Australia's most productive prawn fishery: the Northern Prawn Fishery. Commercial trawl logbooks for 1970-93 and research surveys are used to develop population models for these prawns. A population model that incorporates continuous recruitment is developed. Annual spawning stock and recruitment indices are then estimated from the population model. Spawning stock indices represent the abundance of female prawns that are likely to spawn; recruitment indices represent the abundance of all prawns less than a certain size. The relationships between spawning stock and subsequent recruitment (SRR), between recruitment and subsequent spawning stock (RSR), and between recruitment and commercial catch were estimated through maximum-likelihood models that incorporated autoregressive terms. Yield as a function of fishing effort was estimated by constraining to equilibrium the SRR and RSR. The resulting production model was then used to determine maximum sustainable yield (MSY) and its corresponding fishing effort (f(MSY)). Long-term yield estimates for the two tiger prawn species range between 3700 and 5300 t. The fishing effort at present is close to the level that should produce MSY for both species of tiger prawns. However, current landings, recruitment and spawning stock are below the equilibrium values predicted by the models. This may be because of uncertainty in the spawning stock-recruitment relationships, a change in carrying capacity, biased estimates of fishing effort, unreliable catch statistics, or simplistic assumptions about stock structure. Although our predictions of tiger prawn yields are uncertain, management will soon have to consider new measures to counteract the effects of future increases in fishing effort.

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.

A spatially explicit multi-competitor coexistence model of penaeid (shrimp) distribution on the Australian Great Barrier Reef

A spatially explicit multi-competitor coexistence model was developed for meta-populations of prawns (shrimp) occupying habitat patches across the Great Barrier Reef, where dispersal was localised and dispersal rates varied between species. Prawns were modelled as individuals moving to and from patches or cells according to pre-set decision rules. The landscape was simulated as a matrix of cells with each cell having a spatially explicit survival index for each species. Mixed species prawn assemblages moved over this simplified spatially explicit landscape. A low level of chronic random environmental disturbance was assumed (cyclone and tropical storm damage) with additional acute spatially confined disturbance due to commercial trawling, modelled as an increase in mortality affecting inter-specific competition. The general form of the results was for increased disturbance to favour good-colonising "generalist" species at the expense of good-competitor "specialists". Increasing fishing mortality (local patch extinctions) combined with poor colonising ability resulted in low equilibrium abundance for even the best competitor, while in the same circumstances the poorest competitor but best coloniser could have the highest equilibrium abundance. This mimics the switch from high-value prawn species to lower-value prawn species as trawl effort increases, reflected in historic catch and effort logbook data and reported anecdotaly from the north Queensland trawl fleet. To match the observed distribution and behaviour of prawn assemblages, a combination inter-species competition, a spatially explicit landscape, and a defined pattern of disturbance (trawling) was required. Modelling this combination could simulate not only general trends in spatial distribution of each of prawn species but also localised concentrations observed in the survey data

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.

Additional records of the rock shrimp, Sicyonia australiensis Hanamura & Wadley, 1998 (Decapoda, Sicyoniidae) from eastern Australia, western central Pacific

This paper reports on the collection of S. australiensis from the continental shelf off southern Queensland, easter Australia, in the western Central Pacific, documenting for the first time the occurrence of the species outside of eastern Bass Strait.

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.