Su-Zan Yeh

Modelling the growth of crustacean species

Crustaceans play an important role in marine ecosystem and worldwide fisheries. Accurate and quantitative description of growth is crucial in modelling the demographics and fisheries stock assessment. The stepwise growth as a result of the moulting process and the lack of permanent calcified structures make the traditional approaches developed for finfish inappropriate for crustaceans. This study reviews the data collection, quantitative methods and various sources of uncertainty for modelling the growth of crustacean. The methods were reviewed based on their data requirement and mathematical/statistical complexity ranging from simple growth rate, continuous growth curve to stepwise growth curve and size transition matrix. A comparative example of different growth models was illustrated with four selected crustaceans: American lobster (Homarus americanus, Nephropidae), Dungeness crab (Metacarcinus magister, Cancridae), blue crab (Callinectes sapidus, Portunidae) and pronghorn spiny lobster (Panulirus penicillatus, Palinuridae). Model selection and multi-model inference based on the information theory was discussed for modelling the growth components of moutling increment and intermoult period. Various drivers in determination of the growth pattern of crustaceans were evaluated, including life history strategy and environmental stress. The incorporation of environmental factors into the development of growth models and the recent progress in age-determination were highlighted. We discussed possible research needs for better quantification of crustaceans growth.

Habitat suitability analysis and identification of potential fishing grounds for swordfish, Xiphias gladius, in the South Atlantic Ocean

Swordfish, Xiphias gladius, is a highly migratory species of important commercial value and widely distributed in three oceans. Recently, the South Atlantic swordfish captured as by-catch in longline fisheries targeting tunas has contributed greatly to the overall Atlantic swordfishs landing. In this study, we have developed a habitat suitability index (HSI) model to examine the relationships between their spatio-temporal distribution and environmental factors and to identify potential fishing grounds for the swordfish in the South Atlantic Ocean using the Taiwanese distant-water longline fishery data and remote-sensing oceanographic data for 1998–2007. All the environmental factors considered – sea surface temperature (SST), mixed layer depth (MLD), sea surface height anomaly (SSHA), chlorophyll-a concentration (CHA) and ocean bathymetry (BAH) – were highly significant with most of the catch-per-unit-effort (CPUE) variation explained by SST. The most optimum habitat (i.e. hotspot) was found in the areas with SSTs of 27–28°C, SSHAs of −0.05 to 0.05 m, CHAs of 0.1–0.2 mg m−3 and BAHs of −4000 to −4500 m. The arithmetic mean model with five environmental variables was found to be the most appropriate according to the information theory based on the evaluation of different empirical HSI models in combination with different environmental factors. The bimonthly geographic information system maps of the predicted HSI values were cross-validated by the observed CPUE, suggesting that the model can be used as a tool for reliable prediction of potential fishing grounds. Because the distribution and relative abundance of swordfish are sufficiently heterogeneous in space and time, the output of this study could provide a scientific basis for time–area closures based management of this species.

Application of the sex-specific age-structured assessment method for swordfish, Xiphias gladius, in the North Pacific Ocean

Swordfish are known to be sexually dimorphic. However, previous assessments of the status of swordfish in the North Pacific Ocean have ignored his. A sex-specific age-structured assessment model was therefore constructed and fitted to catch, catch-rate and length–frequency data for the wordfish fisheries that operate in the North Pacific Ocean. Except if natural mortality is lower than its “best” estimate, the results indicate that the pawning stock biomass in 2002 was at a high fraction of its unfished level and that the fishing intensity in 2002 was less than FMSY. Therefore, the wordfish stock in the North Pacific Ocean appears to be relatively stable at the current level of exploitation. However, the results of the assessment odel are sensitive to the values for natural mortality and the steepness of the stock-recruitment relationship. Forward projections based on samples rom a Bayesian posterior distribution indicate that there is negligible risk of the stock dropping below 40% of the unfished spawning stock biomass f fleet-aggregated fishing intensity remains at the current level. However, the risk of population depletion is higher if natural mortality is lower han the “best” estimate. 2006 Elsevier B.V. All rights reserved.

Reproductive biology of female bigeye tuna Thunnus obesus in the western Pacific Ocean.

The reproductive biology of female bigeye tuna Thunnus obesus was assessed by examining 888 fish (ranging from 84·9 to 174·4 cm fork length, LF ) caught by Taiwanese offshore longliners in the western Pacific Ocean from November 1997 to November 1998 and November to December 1999 and 258 gonad samples from these fish. The overall sex ratio of the catch during the sampling differed significantly from 0·5, but males were predominant in sizes >140 cm LF . Reproductive activity (assessed by histology), a gonado-somatic index, and the size-frequency distributions of whole oocytes indicated that spawning occurred throughout the year and the major spawning season appeared to be from February to September. The estimated sizes at 50% maturity (LF50 ) of females was 102·85 cm (95% c.i.: 90·79-110·21 cm) and the smallest mature female was 99·7 cm LF . They are multiple spawners and oocytes develop asynchronously. The proportion of mature (0·63) and reproductively active (0·70) females with ovaries containing postovulatory follicles indicated that they spawn almost daily. Batch fecundity for 15 females with the most advanced oocytes (>730 µm) ranged from 0·84 to 8·56 million eggs (mean ± s.d. = 3·06 ± 2·09). The relationships between batch fecundity (FB , in millions of eggs) and LF (cm) and round mass (MR , kg) were FB=9·91×10-14LF6·38 (r(2)  = 0·84) and FB=8·89×10-4MR2·05 (r(2)  = 0·80), respectively. The parameters estimated in this study are key information for stock assessments of T. obesus in the western Pacific Ocean and will contribute to the conservation and sustainable yield of this species.

Spatio-temporal distributions of tuna species and potential habitats in the Western and Central Pacific Ocean derived from multi-satellite data

In this study, potential tuna habitats in the Western and Central Pacific Ocean (WCPO) area from 1998 to 2005 were identified using a combination of multi-satellite oceanographic data related to sea-surface temperature (SST), chlorophyll-a concentration (Chl-a) and sea-surface height anomaly (SSHA). Histogram analyses of tuna catch per unit effort (CPUE) by the Taiwanese purse seine fishery in relation to the satellite oceanographic data were used to determine tunas preferred ranges of the three satellite-derived oceanographic parameters. The highest CPUE corresponded to areas where SST ranged from 29.0–29.5°C, Chl-a ranged from 0.0–0.1 mg m−3 and SSHA ranged from 0.0–5.0 cm. Local areas within the WCPO with similar satellite-derived oceanographic parameters were assumed to be potential tuna habitat, and a spatially explicit model was developed to determine their locations. It was found that the potential tuna habitats were distributed in warm pools, in the vicinity of the western cold tongue in the equatorial Pacific. Consequently, by using these potential habitats, tuna-fishery management and protection of habitats of tuna species can be improved, and these can possibly form the basis of fishing-ground forecasting in the future.

Incorporating habitat preference into the stock assessment and management of blue marlin (Makaira nigricans) in the Pacific Ocean

Stock assessments that include a spatial component or relate population dynamics to environmental conditions can be considered one way of implementing an ecosystem approach to fisheries. A spatially-structured population dynamics model that takes account of habitat preference is developed and then applied to Pacific blue marlin (Makaira nigricans), as they prefer certain habitats and migrate seasonally. The model is fitted to fishery catch-rate and size data, along with information on the relative density of the population over space derived from a habitat preference model fitted to oceanographic and biological variables. Results show that blue marlin are more abundant in tropical waters, and females account for most of the biomass. Assessments that allow for environmental factors, movement dynamics and sexual dimorphism indicate that this population is in an over-exploited state, with current spawning stock biomass below the level corresponding to maximum sustainable yield (SMSY) and current fishing mortality exceeding that needed to achieve MSY (FMSY). A risk analysis based on samples from a Bayesian posterior distribution suggests that the population will remain above SMSY after 20 years if exploitation rates are below the level corresponding to FMSY.

Incorporating uncertainty into the estimation of biological reference points for a spiny lobster (Panulirus penicillatus) fishery

Abstract The status of a fishery is often defined as the probability of fishing mortality rate exceeding a perilous level for long‐term sustainability. Lobster stock assessments are often subject to large uncertainty in input data and high levels of natural variability in lobster life history processes, which calls for incorporating uncertainty associated with both indicator and management reference points in an evaluation of biological risk of overfishing. Using a Monte Carlo simulation approach, we evaluated the impacts of uncertainty in modelling on the determination of the status of the Taitung spiny lobster (Panulirus penicillatus) fishery (Taiwan), which has not been quantitatively determined despite its commercial importance. The commonly used biological reference points derived from the per recruit model (F 0.1 the fishing mortality rate where the slope of the curve of yield‐per‐recruit model is 10% of the maximum slope and F 4Q%, the fishing mortality rate that reduces the expected egg production for a cohort of female lobsters to 40% of that produced in the absence of a fishery of the egg‐per‐recruit model) were influenced by uncertainties associated with lobster life history and fishery parameters. A large uncertainty in the current fishing mortality rate (F cnr) and estimates of biological reference points (F BRPs) increased the uncertainty in determining the risk of overexploitation throughout the confidence levels of the stochastic decision‐making framework. This simulation study suggests that the target reference point of F 40% is less sensitive to the input parameters’ uncertainty than F 0.1 We suggest a further evaluation of other F‐based references points and development of biomass‐based reference points before final selection and implementation for the management of the Taitung lobster fishery.

Trophic size‐structure of sailfish Istiophorus platypterus in eastern Taiwan estimated by stable isotope analysis

To examine trophic dynamics over different size classes, an isotopic study of sailfish Istiophorus platypterus life-history stages was carried out. Samples were collected from eastern Taiwan and the South China Sea during April 2009 and February 2012. A total of 263 samples (111-245 cm, lower jaw fork length, LLJFL ) were examined for changes in trophic structure in relation to LLJFL by using stable isotope analysis of carbon (δ(13) C) and nitrogen (δ(15) N). The δ(15) N values for I. platypterus ranged from 7·51 to 14·19‰ (mean ± s.d. = 12·06 ± 1·16‰) and the δ(13) C values ranged from -22·04 to -15·48‰ (mean ± s.d. = -17·62 ± 1·10‰). The δ(15) N values were positively dependent on LLJFL (r(2)  = 0·377), whereas δ(13) C were negatively dependent on LLJFL (r(2)  = 0·063). There were significantly different seasonal changes in nitrogen and carbon isotopic concentration, but no significant differences in concentrations between eastern Taiwan and the South China Sea were reported. The trophic level (TL ) of each LLJFL class was correlated, starting from 2·84 TL for size class I (LLJFL  < 140 cm) and reaching 5·03 TL for size class VI (LLJFL > 221 cm). The mean ± s.d. TL was 4·43 ± 0·19 for all samples. The results reveal that I. platypterus occupies a wide range of trophic levels and different size classes occupy different trophic positions in the pelagic ecosystem.

LENGTH-BASED ESTIMATES OF GROWTH AND NATURAL MORTALITY FOR BLUE MARLIN (Makaira nigricans) IN THE NORTHWEST PACIFIC OCEAN

Blue marlin is an important resource for commercial and recreational fisheries. However, several life history parameters needed for stock assessments of this species are poorly determined. Sex-specific catch-at-size data (eye fork length, EFL) for blue marlin were collected from the Taiwanese offshore longline fishery in the northwest Pacific Ocean, and analyzed using the MULTIFAN, a length-frequency analysis, to derive growth parameters. The von Bertalanffy growth parameters of blue marlin were estimated to be different between sexes. The best model for the females included 11 age-classes in the length-frequency data sets (asymptotic length L∞ = 312.5 cm EFL, growth coefficient k = 0.111 yr^(-1), and theoretical age at zero length t_0 = -2.42 yr), while 9 age-classes were identified for the males (L∞ = 232.8 cm EFL, k = 0.131 yr^(-1), and t_0 = -3.58 yr). Natural mortality rates, based on empirical equations, were estimated to be 0.258 yr^(-1) for the males and 0.213 yr^(-1) for the females, respectively. The estimates of growth parameters and mortality rates derived from this study could be used in stock assessments for blue marlin, and contribute towards the fisheries management of the species.

The Kuroshio variations from satellite-derived sea surface temperature and Argos satellite-tracking Lagrangian drifters

This study aimed to identify the spatial and temporal variations of the Kuroshio using satellite-derived Sea Surface Temperature (SST) data and a total of 585 Argos satellite-tracking Lagrangian drifters during 1987–2006 in the northwestern Pacific Ocean. The threshold of high surface velocity, more than 75 cm s–1, was used to determine the location of the Kuroshio Current and found that its width ranged from 50 to 100 km. Moreover, this study first conducted a synoptic view for three predominant regions of the Kuroshio, including the area east of Taiwan, the area east of the East China Sea (ECS) and the area south of Japan. We found that the Kuroshio flows at a three-step increase in speed from upstream to downstream. The maximum mean surface velocity, 153 cm s–1, occurred to the south of Japan. Furthermore, the higher surface velocity of the Kuroshio was usually during the summertime (from April to September). In addition, the four-cruise CTD (Conductivity, Temperature, Depth)-transect data observed in January 2005 were collected to identify and validate the path of the Kuroshio axes. The results showed that the measure of 17.5°C at a depth of 200 m could be regarded as an indicator of the Kuroshio axis, corresponding with the maximum surface velocities of the Argos drifters. The monthly satellite-derived SST values of the Kuroshio were usually more than 10°C, so it was difficult to find a specific isotherm to express the location of the Kuroshio. However, extensive Kuroshio fronts derived from satellite SST data were found in this study. These Kuroshio fronts, associated with high SST gradients, were usually located on the boundary between the warm water of the Kuroshio and the colder water of the shelf, corresponding with the boundary of high surface velocity of the Argos drifters. Consequently, this study suggests that these Kuroshio fronts can be regarded as an effective indicator of the Kuroshio paths and its main axes.