Fisheries science is one of the disciplines of studying marine and riverine organisms, but it is different from usual biology in that it focuses on how to utilize organisms sustainably as renewable ...resources. In terms of sustainable use of fishery resources, it is important to make necessary decisions as soon as possible. However, large uncertainty involved with the fisheries data usually leads to inconclusive argument, resulting in inadequate decision making and the risk of collapsing the resource without doing possible measures. In the latter half of the 20th century, an emphasis on fisheries management science moved away from optimality to robustness. In order to achieve robustness against various uncertainties, biometric statistical methods utilizing the power of modern high-speed computers were actively incorporated into the assessment and risk management of fishery resources. This paper introduces how biometric statistical methods have been applied in recent years in fisheries science, especially in the context of recent development of the MSY-based fisheries management in Japan.
We present the first quantitative review of the stock status relative to the stock biomass (B) and the exploitation rate (U) that achieved the maximum sustainable yield (MSY) (BMSY and UMSY, ...respectively) for 37 Japanese stocks contributing 61% of the total marine capture production in Japan. BMSY and UMSY were estimated by assuming three types of stock-recruitment (S-R) relationships and an age-structured population model or by applying a surplus production model. The estimated stock status shows that approximately half of the stocks were overfishing (U/UMSY > 1), and approximately half of the stocks were overfished (B/BMSY < 0.5) during 2011–2013. Over the past 15 years, U decreased and B slightly increased on average. The rate of decrease in the U of the stocks managed by the total allowable catch (TAC) was significantly greater than that of the other stocks, providing evidence of the effectiveness of TAC management in Japan. The above statuses and trends were insensitive to the assumption of the S-R relationship. The characteristics of Japanese stocks composed mainly of resources with relatively high natural mortality, i.e. productivity, suggest that Japanese fisheries have great potential of exhibiting a quick recovery and increasing their yield by adjusting the fishing intensity to an appropriate level.
Food contamination caused by radioisotopes released from the Fukushima Dai-ichi nuclear power plant is of great public concern. The contamination risk for food items should be estimated depending on ...the characteristics and geographic environments of each item. However, evaluating current and future risk for food items is generally difficult because of small sample sizes, high detection limits, and insufficient survey periods. We evaluated the risk for aquatic food items exceeding a threshold of the radioactive cesium in each species and location using a statistical model. Here we show that the overall contamination risk for aquatic food items is very low. Some freshwater biota, however, are still highly contaminated, particularly in Fukushima. Highly contaminated fish generally tend to have large body size and high trophic levels.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
Pacific saury (Cololabis saira) has a 2‐year life span and age‐1 fish migrates from the central and western North Pacific to Japanese waters from summer to winter. To understand interannual abundance ...variability of Pacific saury in the North Pacific, we examined the extended Japanese standardized catch per unit effort (esCPUE) with generalized linear models (GLMs) during 1982–2018. Explanatory variables included proxies of total and non‐traditional fishing effort (Ettl and Ent) and environmental factors. Each environmental variable was averaged for years t ‐1 and t, corresponding to esCPUE in year t, as esCPUE represented two‐year classes. Ettl and Ent were averaged for years t ‐ 2 and t‐1, assuming fishing of year t affects biomass in year t + 1. We compared four GLMs with different explanatory variables: (a) environmental factors only, (b) fishing effort (without interactions) only, (c) environmental factors and fishing effort (without interactions) and (d) environmental factors and fishing effort (with interactions, i.e., time‐varying effects). Explanatory variables of the best GLM in terms of AICc included the following: Ent with different effects between 1982–2015 and 2016–2018, sea surface water temperature (SST) of the Kuroshio Recirculation Area (KRA) in winter, North Pacific Gyre Oscillation (NPGO) in winter, Southern Oscillation Index in winter, and the biomass of Japanese sardine. The NPGO, a proxy of zooplankton abundance, and KRA‐SST may affect the early growth and survival rates of Pacific saury. We hypothesized the abundance variability of Pacific saury was driven by environmental factors and fishing, the latter of which more adversely affected since 2016.
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, UILJ, UKNU, UL, UM, UPUK
Abstract
Nonlinear phenomena are universal in ecology. However, their inference and prediction are generally difficult because of autocorrelation and outliers. A traditional least squares method for ...parameter estimation is capable of improving short-term prediction by estimating autocorrelation, whereas it has weakness to outliers and consequently worse long-term prediction. In contrast, a traditional robust regression approach, such as the least absolute deviations method, alleviates the influence of outliers and has potentially better long-term prediction, whereas it makes accurately estimating autocorrelation difficult and possibly leads to worse short-term prediction. We propose a new robust regression approach that estimates autocorrelation accurately and reduces the influence of outliers. We then compare the new method with the conventional least squares and least absolute deviations methods by using simulated data and real ecological data. Simulations and analysis of real data demonstrate that the new method generally has better long-term and short-term prediction ability for nonlinear estimation problems using spawner–recruitment data. The new method provides nearly unbiased autocorrelation even for highly contaminated simulated data with extreme outliers, whereas other methods fail to estimate autocorrelation accurately.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
The abundance and recruitment of anchovy Engraulis spp. and sardine Sardinops spp. alternate in a synchronized way across the Pacific. Convergent cross mapping (CCM) indicated that climate change ...drives the alternation of the two species in the California Current System. However, climate indices patterns in the western North Pacific contrast with those in the eastern North Pacific, despite synchronous species alternations occurring. Therefore, it is of great interest to clarify whether climate change, or any other factors, affects the population dynamics of Japanese anchovy and Japanese sardine in the western North Pacific. Using CCM, we tested whether climate change and interspecific interactions affect the population dynamics of these two species. We found that climate change affected recruitment, and we clarified the spatiotemporal pattern of this effect. This result supports the existing hypotheses that population dynamics are regulated by climate change in the western North Pacific. The present study also detected interspecific interactions between sardine and anchovy, which might promote the alternation of the two species, and has not been reported in other regions.
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, UILJ, UKNU, UL, UM, UPUK
Standardized catch per unit effort (CPUE) is a fundamental component of fishery stock assessment. In multispecies fisheries, catchability can differ depending on which species is being targeted, and ...so the yearly trend extracted from the standardized CPUE is likely to be biased. We have, therefore, developed a method for predicting the unobserved variable related to targeted species from among multispecies composition data using a mixture regression model for the transformed residuals. In contrast with traditional methods, the proposed method predicts the target variable in CPUE standardization without removing a subset of the data. Keeping the entire data set avoids information loss, and so CPUE standardization with the predicted target variable should yield an unbiased estimate of the yearly trend. Simple simulation tests demonstrate that our method outperforms traditional methods. We illustrate the use of our method by applying it to CPUE data on arabesque greenling (Pleurogrammus azonus) caught in multispecies trawl fisheries in Hokkaido, Japan.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
A simple age verification method is presented for centrum edge analysis (CEA) of elasmobranch species. In this method, a binomial model is linked with a von Mises distribution for circular data, ...taking the characteristics of the CEA data into account. The periodicity of growth band pairs is categorized as no cycle, an annual cycle, or a biannual cycle. Three models are then constructed according to different periodicities. We use the Akaike information criterion (AIC) to determine which model is the best. The models were applied to the shortfin mako shark (
Isurus oxyrinchus
) data collected in the North Pacific to identify the best periodicity model. The AIC best-fitting model was one involving an annual cycle. The general performance of the method was evaluated using simulated data of various sample sizes, formation times, and durations of growth band pairs. The simulation trials showed that the performance of the method was satisfactory with moderate sample sizes. This method should improve the accuracy of age determination and could be applied to all species that have periodic growth band pairs.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Catchability plays a central role in fisheries stock assessment. Since catchability often varies with time depending on population density, environmental factors, and anthropogenic effects, assuming ...constant catchability in population models can lead to biased abundance estimates. Here we present a novel way to simultaneously estimate time-varying catchability and abundance by integrating a short-term (month-based) removal method and a long-term (year-based) age-structured population dynamics model. We applied this approach to commercial fishery data for a Japanese pufferfish (Takifugu rubripes) population and found that the models with time-varying catchability greatly outperformed the models with constant catchability in terms of predictive ability and model consistency. The temporal variation in catchability was parsimoniously predicted by fishing effort and population size, indicating the existence of effort- and density-dependent catchability. Our approach, integrating population dynamics at different timescales, will help avoid inadvertent overexploitation and contribute to sustainable harvesting by enhancing the estimation accuracy of time-varying catchability and abundance.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract Tuned virtual population analyses are widely used for fisheries stock assessments. However, accurately estimating abundances and fishing mortality coefficients in the terminal year using ...tuned virtual population analyses is generally difficult, particularly when there is a limited number of available abundance indices. We propose a new method of integrating the tuned virtual population analyses with a ridge regression approach. In our method, penalization in the ridge regression is applied to the age-specific fishing mortalities in the terminal year, and the penalty parameter is automatically selected by minimizing the retrospective bias. Therefore, our method is able to simultaneously obtain a stable estimation of fishing mortality coefficients in the terminal year and reduce retrospective bias. Simulation tests based on the northern Japan Sea stock of walleye pollock (Gadus chalcogrammus) in the Sea of Japan demonstrated that this method yielded less biased estimates of abundances and avoided overestimations of fishing mortality coefficients in the terminal year. In addition, despite limited abundance indices, our method can perform reliable abundance estimations even under hyperstability and hyperdepletion conditions.