skip to main content

Title: Non-linearity in stock–recruitment relationships of Atlantic cod: insights from a multi-model approach
Abstract The stock–recruitment relationship is the basis of any stock prediction and thus fundamental for fishery management. Traditional parametric stock–recruitment models often poorly fit empirical data, nevertheless they are still the rule in fish stock assessment procedures. We here apply a multi-model approach to predict recruitment of 20 Atlantic cod (Gadus morhua) stocks as a function of adult biomass and environmental variables. We compare the traditional Ricker model with two non-parametric approaches: (i) the stochastic cusp model from catastrophe theory and (ii) multivariate simplex projections, based on attractor state-space reconstruction. We show that the performance of each model is contingent on the historical dynamics of individual stocks, and that stocks which experienced abrupt and state-dependent dynamics are best modelled using non-parametric approaches. These dynamics are pervasive in Western stocks highlighting a geographical distinction between cod stocks, which have implications for their recovery potential. Furthermore, the addition of environmental variables always improved the models’ predictive power indicating that they should be considered in stock assessment and management routines. Using our multi-model approach, we demonstrate that we should be more flexible when modelling recruitment and tailor our approaches to the dynamical properties of each individual stock.
; ; ; ; ; ;
Griffith, Gary
Award ID(s):
1655203 1660584
Publication Date:
Journal Name:
ICES Journal of Marine Science
Page Range or eLocation-ID:
1492 to 1502
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract From hillslope to small catchment scales (< 50 km 2 ), soil carbon management and mitigation policies rely on estimates and projections of soil organic carbon (SOC) stocks. Here we apply a process-based modeling approach that parameterizes the MIcrobial-MIneral Carbon Stabilization (MIMICS) model with SOC measurements and remotely sensed environmental data from the Reynolds Creek Experimental Watershed in SW Idaho, USA. Calibrating model parameters reduced error between simulated and observed SOC stocks by 25%, relative to the initial parameter estimates and better captured local gradients in climate and productivity. The calibrated parameter ensemble was used to produce spatially continuous, high-resolution (10 m 2 ) estimates of stocks and associated uncertainties of litter, microbial biomass, particulate, and protected SOC pools across the complex landscape. Subsequent projections of SOC response to idealized environmental disturbances illustrate the spatial complexity of potential SOC vulnerabilities across the watershed. Parametric uncertainty generated physicochemically protected soil C stocks that varied by a mean factor of 4.4 × across individual locations in the watershed and a − 14.9 to + 20.4% range in potential SOC stock response to idealized disturbances, illustrating the need for additional measurements of soil carbon fractions and their turnover time to improve confidence in the MIMICS simulations of SOCmore »dynamics.« less
  2. In practice, lead time demand (LTD) can be non-standard: skewed, multi-modal or highly variable; factors that compromise the validity of the classic approaches for setting safety stock levels. Motivated by encountering this problem at our industry partner, we develop an approach for setting safety stock levels using the bootstrap, a widely-used statistical procedure. We extend prior research that has used the bootstrap for quantile estimation to address the multi-parameter estimation of safety stocks. We develop a multivariate central limit theorem for the bootstrap mean and bootstrap quantile -- components of the safety stock calculation -- highlighting why the generalization of these bootstrap methods is critical for inventory management. These results provide a theoretical underpinning for the bootstrap estimator of safety stock and permit the construction of confidence intervals for safety stock estimates, allowing decision makers to understand the reliability with which the desired service level will be achieved. Building on our theoretical results, and supported by numerical experiments, we provide insights on the behavior of the bootstrap for various LTD distributions, which our results demonstrate are critical when employing the bootstrap method. Implementation results with our industry partner indicate our approach is quite effective in setting safety stock levels.
  3. Abstract

    Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatiotemporal scale. We propose a spectral error model that accounts for the neighboring space–time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency–wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space–time linear filtering term. The random errors are represented as a nonstationary additive noise. The spectral error model is applied to the IMERG multisatellite precipitation product, and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge–radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass (attenuating the fine-scale variability). While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variancemore »at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from different sources shall not be expected to automatically have statistically independent errors.

    Significance Statement

    Satellite precipitation products are nowadays widely used for climate and environmental research, water management, risk analysis, and decision support at the local, regional, and global scales. For all these applications, knowledge about the accuracy of the products is critical for their usability. However, products are not systematically provided with a quantitative measure of the uncertainty associated with each individual estimate. Various parametric error models have been proposed for uncertainty quantification, mostly assuming that the uncertainty is only a function of the precipitation intensity at the pixel and time of interest. By projecting satellite precipitation fields and their retrieval errors into the Fourier frequency–wavenumber domain, we show that we can explicitly take into account the neighboring space–time multiscale dynamics of precipitation and compute a scale-dependent uncertainty.

    « less
  4. Accurate age data are essential for reliable fish stock assessment. Yet many stocks suffer from inconsistencies in age interpretation. A new approach to obtain age makes use of the chemical composition of otoliths. This study validates the periodicity of recurrent patterns in 25 Mg, 31 P, 34 K, 55 Mn, 63 Cu, 64 Zn, 66 Zn, 85 Rb, 88 Sr, 138 Ba, and 208 Pb in Baltic cod (Gadus morhua) otoliths from tag–recapture and known-age samples. Otolith P concentrations showed the highest consistency in seasonality over the years, with minima co-occurring with otolith winter zones in the known-age otoliths and in late winter – early spring when water temperatures are coldest in tagged cod . The timing of minima differs between stocks, occurring around February in western Baltic cod and 1 month later in eastern Baltic cod; seasonal maxima are also stock-specific, occurring in August and October, respectively. The amplitude in P is larger in faster-growing western compared with eastern Baltic cod. Seasonal patterns with minima in winter – late spring were also evident in Mg and Mn, but less consistent over time and fish size than P. Chronological patterns in P, and to a lesser extent Mg and Mn, may have themore »potential to supplement traditional age estimation or to guide the visual identification of translucent and opaque otolith patterns used in traditional age estimation.« less
  5. Understanding population dynamics is essential for achieving sustainable and productive fisheries. However, estimating recruitment in a stock assessment model involves the challenging task of identifying a self-sustaining population, which often includes representing complex geographic structure. A review of several case studies demonstrated that alternative stock assessment models can influence estimates of recruitment. Incorporating spatial population structure and connectivity into stock assessment models changed the perception of recruit- ment events for a wide diversity of fisheries, but the degree to which estimates were impacted depended on movement rates and relative stock sizes. For multiple population components, estimates of strong recruitment events and the productivity of smaller population units were often more sensitive to connectivity assumptions. Simulation testing, conditioned on these case studies, suggested that accurately accounting for population structure, either in management unit definitions or stock assessment model structure, improved recruitment estimates. An understanding of movement dynamics improved estimation of connected sub-populations. The challenge of representing geographic structure in stock assessment emphasizes the importance of defining self- sustaining management units to justify a unit-stock assumption.