Stock assessment models are fitted to abundance-index, fishery catch, and age–length–sex composition data that are estimated from survey and fishery records. Research has developed spatiotemporal ...methods to estimate abundance indices, but there is little research regarding model-based methods to generate age–length–sex composition data. We demonstrate a spatiotemporal approach to generate composition data and a multinomial sample size that approximates the estimated imprecision. A simulation experiment comparing spatiotemporal and design-based methods demonstrates a 32% increase in input sample size for the spatiotemporal estimator. A Stock Synthesis assessment used to manage lingcod (Ophiodon elongatus) in the California Current also shows a 17% increase in sample size and better model fit using the spatiotemporal estimator, resulting in smaller standard errors when estimating spawning biomass. We conclude that spatiotemporal approaches are feasible for estimating both abundance-index and compositional data, thereby providing a unified approach for generating inputs for stock assessments. We hypothesize that spatiotemporal methods will improve statistical efficiency for composition data in many stock assessments and recommend that future research explore the impact of including additional habitat or sampling covariates.
Estimating trends in abundance from fishery catch rates is one of the oldest endeavors in fisheries science. However, many jurisdictions do not analyze fishery catch rates due to concerns that these ...data confound changes in fishing behavior (adjustments in fishing location or gear operation) with trends in abundance. In response, we developed a spatial dynamic factor analysis (SDFA) model that decomposes covariation in multispecies catch rates into components representing spatial variation and fishing behavior. SDFA estimates spatiotemporal variation in fish density for multiple species and accounts for fisher behavior at large spatial scales (i.e., choice of fishing location) while controlling for fisher behavior at fine spatial scales (e.g., daily timing of fishing activity). We first use a multispecies simulation experiment to show that SDFA decreases bias in abundance indices relative to ignoring spatial adjustments and fishing tactics. We then present results for a case study involving petrale sole (Eopsetta jordani) in the California Current, for which SDFA estimates initially stable and then increasing abundance for the period 1986-2003, in accordance with fishery-independent survey and stock assessment estimates.
Rockfish are an important component of West Coast fisheries and California Current food webs, and recruitment (cohort strength) for rockfish populations has long been characterized as highly variable ...for most studied populations. Research efforts and fisheries surveys have long sought to provide greater insights on both the environmental drivers, and the fisheries and ecosystem consequences, of this variability. Here, variability in the temporal and spatial abundance and distribution patterns of young-of-the-year (YOY) rockfishes are described based on midwater trawl surveys conducted throughout the coastal waters of California Current between 2001 and 2019. Results confirm that the abundance of winter-spawning rockfish taxa in particular is highly variable over space and time. Although there is considerable spatial coherence in these relative abundance patterns, there are many years in which abundance patterns are very heterogeneous over the scale of the California Current. Results also confirm that the high abundance levels of YOY rockfish observed during the 2014-2016 large marine heatwave were largely coastwide events. Species association patterns of pelagic YOY for over 20 rockfish taxa in space and time are also described. The overall results will help inform future fisheries-independent surveys, and will improve future indices of recruitment strength used to inform stock assessment models and marine ecosystem status reports.
The ability of management strategies to achieve the fishery management goals are impacted by environmental variation and, therefore, also by global climate change. Management strategies can be ...modified to use environmental data using the "dynamic B sub(0)" concept, and changing the set of years used to define biomass reference points. Two approaches have been developed to apply management strategy evaluation to evaluate the impact of environmental variation on the performance of management strategies. The "mechanistic approach" estimates the relationship between the environment and elements of the population dynamics of the fished species and makes predictions for population trends using the outputs from global climate models. In contrast, the "empirical approach" examines possible broad scenarios without explicitly identifying mechanisms. Many reviewed studies have found that modifying management strategies to include environmental factors does not improve the ability to achieve management goals much, if at all, and only if the manner in which these factors drive the system is well known. As such, until the skill of stock projection models improves, it seems more appropriate to consider the implications of plausible broad forecasts related to how biological parameters may change in the future as a way to assess the robustness of management strategies, rather than attempting specific predictions per se.
Major ecological realignments are already occurring in response to climate change. To be successful, conservation strategies now need to account for geographical patterns in traits sensitive to ...climate change, as well as climate threats to species-level diversity. As part of an effort to provide such information, we conducted a climate vulnerability assessment that included all anadromous Pacific salmon and steelhead (Oncorhynchus spp.) population units listed under the U.S. Endangered Species Act. Using an expert-based scoring system, we ranked 20 attributes for the 28 listed units and 5 additional units. Attributes captured biological sensitivity, or the strength of linkages between each listing unit and the present climate; climate exposure, or the magnitude of projected change in local environmental conditions; and adaptive capacity, or the ability to modify phenotypes to cope with new climatic conditions. Each listing unit was then assigned one of four vulnerability categories. Units ranked most vulnerable overall were Chinook (O. tshawytscha) in the California Central Valley, coho (O. kisutch) in California and southern Oregon, sockeye (O. nerka) in the Snake River Basin, and spring-run Chinook in the interior Columbia and Willamette River Basins. We identified units with similar vulnerability profiles using a hierarchical cluster analysis. Life history characteristics, especially freshwater and estuary residence times, interplayed with gradations in exposure from south to north and from coastal to interior regions to generate landscape-level patterns within each species. Nearly all listing units faced high exposures to projected increases in stream temperature, sea surface temperature, and ocean acidification, but other aspects of exposure peaked in particular regions. Anthropogenic factors, especially migration barriers, habitat degradation, and hatchery influence, have reduced the adaptive capacity of most steelhead and salmon populations. Enhancing adaptive capacity is essential to mitigate for the increasing threat of climate change. Collectively, these results provide a framework to support recovery planning that considers climate impacts on the majority of West Coast anadromous salmonids.
Environmental conditions can have spatially complex effects on the dynamics of marine fish stocks that change across life-history stages. Yet the potential for non-stationary environmental effects ...across multiple dimensions, e.g. space and ontogeny, are rarely considered. In this study, we examined the evidence for spatial and ontogenetic non-stationary temperature effects on Pacific hake Merluccius productus biomass along the west coast of North America. Specifically, we used Bayesian additive models to estimate the effects of temperature on Pacific hake biomass distribution and whether the effects change across space or life-history stage. We found latitudinal differences in the effects of temperature on mature Pacific hake distribution (i.e. age 3 and older); warmer than average subsurface temperatures were associated with higher biomass north of Vancouver Island, but lower biomass offshore of Washington and southern Vancouver Island. In contrast, immature Pacific hake distribution (i.e. age 2) was better explained by a nonlinear temperature effect; cooler than average temperatures were associated with higher biomass coastwide. Together, our results suggest that Pacific hake distribution is driven by interactions between age composition and environmental conditions and highlight the importance of accounting for varying environmental effects across multiple dimensions.
This paper investigates environmental drivers of U.S. West Coast petrale sole (Eopsetta jordani) recruitment as an initial step toward developing an environmental recruitment index that can inform ...the stock assessment in the absence of survey observations of age‐0 and age‐1 fish. First, a conceptual life history approach is used to generate life‐stage‐specific and spatio‐temporally specific mechanistic hypotheses regarding oceanographic variables that likely influence survival at each life stage. Seven life history stages are considered, from female spawner condition through benthic recruitment as observed in the Northwest Fisheries Science Center West Coast Groundfish Bottom Trawl Survey (age‐2 fish). The study area encompasses the region from 40 to 48°N in the California Current Ecosystem. Hypotheses are tested using output from a regional ocean reanalysis model outputs and model selection techniques. Four oceanographic variables explained 73% of the variation in recruitment not accounted for by estimates based exclusively on the spawning stock size. Recruitment deviations were (a) positively correlated with degree days during the female precondition period, (b) positively correlated with mixed‐layer depth during the egg stage, (c) negatively correlated with cross‐shelf transport during the larval stage, and (d) negatively correlated with cross‐shelf transport during the benthic juvenile stage. While multiple mechanisms likely affect petrale sole recruitment at different points during their life history, the strength of the relationship is promising for stock assessment and integrated ecosystem assessment applications.
Species distribution models (SDMs) are a common approach to describing species’ space‐use and spatially‐explicit abundance. With a myriad of model types, methods and parameterization options ...available, it is challenging to make informed decisions about how to build robust SDMs appropriate for a given purpose. One key component of SDM development is the appropriate parameterization of covariates, such as the inclusion of covariates that reflect underlying processes (e.g. abiotic and biotic covariates) and covariates that act as proxies for unobserved processes (e.g. space and time covariates). It is unclear how different SDMs apportion variance among a suite of covariates, and how parameterization decisions influence model accuracy and performance. To examine trade‐offs in covariation parameterization in SDMs, we explore the attribution of spatiotemporal and environmental variation across a suite of SDMs. We first used simulated species distributions with known environmental preferences to compare three types of SDM: a machine learning model (boosted regression tree), a semi‐parametric model (generalized additive model) and a spatiotemporal mixed‐effects model (vector autoregressive spatiotemporal model, VAST). We then applied the same comparative framework to a case study with three fish species (arrowtooth flounder, pacific cod and walleye pollock) in the eastern Bering Sea, USA. Model type and covariate parameterization both had significant effects on model accuracy and performance. We found that including either spatiotemporal or environmental covariates typically reproduced patterns of species distribution and abundance across the three models tested, but model accuracy and performance was maximized when including both spatiotemporal and environmental covariates in the same model framework. Our results reveal trade‐offs in the current generation of SDM tools between accurately estimating species abundance, accurately estimating spatial patterns, and accurately quantifying underlying species–environment relationships. These comparisons between model types and parameterization options can help SDM users better understand sources of model bias and estimate error.
Concurrent declines in demersal fish stock abundances and shifts in long-term average environmental conditions have been well documented in the Pacific. Management advice ignoring environmental ...forcing of recruitment may cause stocks to be over- or under-harvested, so it is important to consider including environmental forcing on recruitment in stock assessment models. Simulation testing is used to determine the statistical power of stock assessments to identify long-term, decadal-scale environmental forcing of recruitment and the ability to estimate management reference points when the duration of the fisheries time series is equal to or less than the period of the environmental cycle. Commonly used assessment methods generally lead to lower total type I (incorrectly rejecting the null hypothesis of no environmental impact) and type II (failing to detect an environmental impact when such a relationship exists) error rates. The promise of integrating environmental data with decadal-scale variability directly into stock assessments via the stock–recruitment relationship is outweighed by the pitfall of high type I error rates that are due to fishing-induced stock declines that coincide with directional environmental change. However, the impact of type I errors can be minimized by choosing an appropriate combination of assessment method and reference point estimators.
Projecting the future distributions of commercially and ecologically important species has become a critical approach for ecosystem managers to strategically anticipate change, but large ...uncertainties in projections limit climate adaptation planning. Although distribution projections are primarily used to understand the scope of potential change—rather than accurately predict specific outcomes—it is nonetheless essential to understand where and why projections can give implausible results and to identify which processes contribute to uncertainty. Here, we use a series of simulated species distributions, an ensemble of 252 species distribution models, and an ensemble of three regional ocean climate projections, to isolate the influences of uncertainty from earth system model spread and from ecological modeling. The simulations encompass marine species with different functional traits and ecological preferences to more broadly address resource manager and fishery stakeholder needs, and provide a simulated true state with which to evaluate projections. We present our results relative to the degree of environmental extrapolation from historical conditions, which helps facilitate interpretation by ecological modelers working in diverse systems. We found uncertainty associated with species distribution models can exceed uncertainty generated from diverging earth system models (up to 70% of total uncertainty by 2100), and that this result was consistent across species traits. Species distribution model uncertainty increased through time and was primarily related to the degree to which models extrapolated into novel environmental conditions but moderated by how well models captured the underlying dynamics driving species distributions. The predictive power of simulated species distribution models remained relatively high in the first 30 years of projections, in alignment with the time period in which stakeholders make strategic decisions based on climate information. By understanding sources of uncertainty, and how they change at different forecast horizons, we provide recommendations for projecting species distribution models under global climate change.
Projecting the future distributions of species has become critical for climate adaption planning, but we don't have a precise understanding of how accurate projections can be without waiting decades for validation. Here, we model simulated species responses under climate change to quantify and understand sources of uncertainty, and how they change at different forecast horizons. We provide recommendations for projecting species distribution models under global climate change.