This paper develops a risk table to facilitate incorporation of additional information into the fisheries stock assessment and management process. The risk table is designed to evaluate unanticipated ...ecosystem and environmental impacts on marine resources that may require a rapid management response. The risk table is a standardized framework to document concerns about the assessment model, population dynamics, and the ecosystem/environment that are not explicitly addressed within the stock assessment model. A scoring procedure is used to evaluate the severity of the concern. These concerns can then be evaluated in support for or against a reduction from the maximum Acceptable Biological Catch while providing reviewers and stakeholders transparent documentation of the concerns. The risk table was applied successfully to several stocks on a trial basis during the 2018 groundfish assessment cycle for the North Pacific Fishery Management Council, and will be used for all full groundfish assessments in 2019. Rapid changes in climate are likely for Alaska marine ecosystems in coming decades, and these changes are not entirely predicable. Therefore, we avocate that the risk table approach should be included in the suite of management tools used to address the effects of climate change on Alaska marine resources.
Catch quotas for walleye pollock Gadus chalcogrammus, the dominant species in the groundfish fishery off Alaska, are set by applying harvest control rules to annual estimates of spawning stock ...biomass (SSB) from age-structured stock assessments. Adult walleye pollock abundance and maturity status have been monitored in early spring in Shelikof Strait in the Gulf of Alaska for almost three decades. The sampling strategy for maturity status is largely characterized as targeted, albeit opportunistic, sampling of trawl tows made during hydroacoustic surveys. Trawl sampling during pre-spawning biomass surveys, which do not adequately account for spatial patterns in the distribution of immature and mature fish, can bias estimated maturity ogives from which SSB is calculated. Utilizing these maturity data, we developed mixed-effects generalized additive models to examine spatial and temporal patterns in walleye pollock maturity and the influence of these patterns on estimates of SSB. Current stock assessment practice is to estimate SSB as the product of annual estimates of numbers at age, weight at age, and mean maturity at age for 1983-present. In practice, we found this strategy to be conservative for a time period from 2003-2013 as, on average, it underestimates SSB by a 4.7 to 11.9% difference when compared to our estimates of SSB that account for spatial structure or both temporal and spatial structure. Inclusion of spatially explicit information for walleye pollock maturity has implications for understanding stock reproductive biology and thus the setting of sustainable harvest rates used to manage this valuable fishery.
Some of the longest and most comprehensive marine ecosystem monitoring programs were established in the Gulf of Alaska following the environmental disaster of the Exxon Valdez oil spill over 30 years ...ago. These monitoring programs have been successful in assessing recovery from oil spill impacts, and their continuation decades later has now provided an unparalleled assessment of ecosystem responses to another newly emerging global threat, marine heatwaves. The 2014-2016 northeast Pacific marine heatwave (PMH) in the Gulf of Alaska was the longest lasting heatwave globally over the past decade, with some cooling, but also continued warm conditions through 2019. Our analysis of 187 time series from primary production to commercial fisheries and nearshore intertidal to offshore oceanic domains demonstrate abrupt changes across trophic levels, with many responses persisting up to at least 5 years after the onset of the heatwave. Furthermore, our suite of metrics showed novel community-level groupings relative to at least a decade prior to the heatwave. Given anticipated increases in marine heatwaves under current climate projections, it remains uncertain when or if the Gulf of Alaska ecosystem will return to a pre-PMH state.
Abstract
Climate-driven changes in the timing of spawning or migration can affect the availability of fish to surveys designed to monitor their abundance, complicating efforts to assess stock status ...and sustainably manage fisheries. From 2017 to 2019, trends in biomass estimates from four surveys used to monitor Gulf of Alaska pollock diverged. These conflicting trends increased uncertainty in the stock assessment and occurred during a time of rapid environmental change. We hypothesized that changes in spawn timing affected availability of pollock to a winter survey that targets pre-spawning aggregations. To test this, we reconstructed relative spawn timing using two independent data sources: spring larval surveys and observations of spawning state in mature female pollock. We found that changes in spawn timing relative to survey timing explained a significant portion of recent and historical discrepancies between survey and model estimates of biomass. We then incorporated measures of spawn timing/survey timing mismatch as catchability covariates in an enhanced state-space stock assessment model. Including spawn timing-based catchability covariates significantly improved the model fit to survey data and provided a mechanistic explanation for recent survey discrepancies.
Information on the annual variability in abundance and growth of juvenile groundfish can be useful for predicting fisheries stocks, but is often poorly known owing to difficulties in sampling fish in ...their first year of life. In the Western Gulf of Alaska (WGoA) and Eastern Bering Sea (EBS) ecosystems, three species of puffin (tufted and horned puffin, Fratercula cirrhata, Fratercula corniculata, and rhinoceros auklet, Cerorhinca monocerata, Alcidae), regularly prey upon (i.e., “sample”) age‐0 groundfish, including walleye pollock (Gadus chalcogramma, Gadidae) and Pacific cod (Gadus microcephalus, Gadidae). Here, we test the hypothesis that integrating puffin dietary data with walleye pollock stock assessment data provides information useful for fisheries management, including indices of interannual variation in age‐0 abundance and growth. To test this hypothesis, we conducted cross‐correlation and regression analyses of puffin‐based indices and spawning stock biomass (SSB) for the WGoA and EBS walleye pollock stocks. For the WGoA, SSB leads the abundance of age‐0 fish in the puffin diet, indicating that puffins sample the downstream production of the WGoA spawning stock. By contrast, the abundance and growth of age‐0 fish sampled by puffins lead SSB for the EBS stock by 1–3 years, indicating that the puffin diet proxies incoming year class strength for this stock. Our study indicates connectivity between the WGoA and EBS walleye pollock stocks. Integration of non‐traditional data sources, such as seabird diet data, with stock assessment data appears useful to inform information gaps important for managing US fisheries in the North Pacific.
Predation can have substantial and long-term effects on the population dynamics of ecologically important prey. Diverse predator assemblages, however, may produce stabilizing (i.e., portfolio) ...effects on prey mortality when consumption varies asynchronously among predators. We calculated spatiotemporal variation in predation on a dominant forage species to quantify synchrony and portfolio effects in a food web context and better understand diversity–stability relationships in a large marine ecosystem that has undergone considerable changes in community composition. We selected Walleye Pollock (Gadus chalcogrammus) as our case study because they support some of the largest, most valuable commercial fisheries in the world and serve as essential prey for an array of economically and culturally important species. Thus, there are sufficient data for Pollock with which to test ecological theories in an empirical setting. Spatially explicit predation indices accounted for annual variation in predator biomass, bioenergetics-based rations, and age-specific proportions of Pollock consumed by a suite of groundfishes in the Gulf of Alaska (1990–2015). We found that Arrowtooth Flounder (Atheresthes stomias) was, by far, the dominant Pollock predator (proportional consumption: 0.74 ± 0.14). We also found synchronous trends in consumption among predator species, indicating a lack of portfolio effects at the basin scale. This combination of a single dominant predator and synchronous consumption dynamics suggests strong top-down control over Pollock in the Gulf of Alaska, though the degree of synchrony was highly variable at all spatial scales. Whereas synchrony generally increased in the western subregion, consumption in the central Gulf of Alaska became less synchronous through time. This suggests diminished trophic stability in one area and increased stability in another, thereby emphasizing the importance of spatiotemporal heterogeneity in maintaining food web structure and function. Finally, total Pollock consumption was highly variable (ranging from 1.87 to 7.63 Tg) and often exceeded assessment-based estimates of productivity. We assert that using our holistic and empirically derived predation index as a modifier of assumed constant natural mortality would provide a practical method for incorporating ecological information into single-species stock assessments.
Abstract
As climate stressors are impacting marine ecosystems and fisheries across the world, ecosystem models that incorporate environmental variables are increasingly used to inform ecosystem-based ...fisheries management. The assumptions around the mechanistic links between climate stressors and the biological processes in these models are important, but the implications for model outcomes of which stressors are captured and how they affect modeled biological processes are seldom explored. Using a whole-ecosystem model (Atlantis) for the Gulf of Alaska, we explore the effects of capturing physical (increased temperature) and biogeochemical (decreased low trophic level productivity) climate stressors, and disentangle the effects of each stressor on the productivity of forage fish, groundfish, and fish-eating seabirds. We then test the effects of alternative model specifications of temperature-driven habitat determination and bioenergetics. Increased temperature resulted in increased weight-at-age and higher natural mortality, while decreased productivity resulted in decreased weight-at-age and higher natural mortality. Model specification of temperature dependence of movement and spawning influenced model outcomes, and decoupling these processes from temperature led to overly optimistic biomass predictions. As the use of ecosystem models to inform fisheries management becomes more operational, we illustrate that the assumptions around the links between climate stressors and ecological processes influence model outcomes.
Climate variability affects marine ecosystems. The mechanisms relating low-frequency environmental fluctuations (regime shifts) and their impacts on marine populations are not well established for ...most species, but there is an expectation that the accuracy of scientific advice provided to fisheries management would be improved if management strategies incorporated the effects of regime shifts on stocks. Management strategy evaluation is used to examine the impact of regime shifts in average recruitment on the performance of management strategies for the fishery for walleye pollock,
Theragra chalcogramma
, in the Gulf of Alaska. The current and four alternative management strategies are evaluated. The alternatives include management strategies with different definitions of the average recruitment used when calculating management reference points, e.g., a 25-year sliding-window method and a "dynamic B
0
" method. The current management strategy, which ignores the possibility of future regime shifts, kept the spawning biomass higher relative to the target level than the other management strategies and had the lowest risk of fishing mortality exceeding the overfishing limit. The sliding-window management strategy achieved the highest catches and the lowest interannual variation in catch, although at the cost of a higher risk of the fishing mortality exceeding the overfishing limit.
Although highly variable in time and space, predation remains the greatest source of mortality for juvenile and lower trophic-level fishes. As such, predation can have substantial and long-term ...effects on the dynamics of these prey. Gulf of Alaska walleye pollock (Gadus chalcogrammus) has shown considerable variability in biomass over the past four decades. During this same time, the demersal fish community transitioned from being dominated by pollock to a system comprised primarily of upper trophic-level predators. We estimated time-varying predation mortality to better understand its effects on the population dynamics of pollock in this currently “top heavy” system. Our index of predation accounted for spatiotemporal variation in predator biomass, bioenergetics-based rations, and age-specific proportions of pollock consumed (1990–2019). To evaluate population-level impacts of predation, we included an index of removals as part of the stock assessment model. This formulation allowed for non-annual data inputs and included a proportionality constant with which to scale predation. Age-specific natural mortality was allowed to vary according to a penalized random walk. We found that natural mortality ranged from 37% higher to 17% lower than the long-term mean. Resulting estimates of total pollock biomass differed by as much as 37% relative to a model without time-varying natural mortality, though the maximum difference for exploitable biomass was only 14%. Using an empirically-derived predation index to modify constant natural mortality allows stock assessment scientists to evaluate impacts of time-varying predation on assessed species. This approach provides a relatively simple way of incorporating ecological information into single-species stock assessments and may reduce bias compared to conventional models that do not account for changes in predation mortality. Notably, including predation mortality in single-species assessments may help identify inconsistencies in biomass estimates that warrant further consideration.
•Time-varying predation mortality impacts prey population dynamics.•Age-structured stock assessments often assume constant natural mortality, M.•We compared effects of constant M and time-varying predation on pollock biomass.•Constant M produced apparent bias, though it was insufficient to result in overfishing.•Additional research is necessary to estimate uncertainty in time-varying predation.
Hierarchical Bayesian meta-analysis can be a useful method for improving estimation of key parameters for harvested fish populations. In hierarchical models, data from multiple populations are used ...simultaneously to obtain estimates of parameters for individual populations and characterize the variability among populations. Many populations of Pacific rockfishes (
Sebastes
spp.) have declined off the US West Coast since the 1980s, and there is also concern for their conservation in Canada. We develop a hierarchical Bayesian meta-analysis to improve estimates of stock-recruit parameters, characterize management-related parameters (e.g., optimal harvest rate), and address uncertainties in the structural form of the stock-recruit function for Pacific rockfishes. We estimate steepness and optimal harvest rates for 14 populations of Pacific rockfishes under alternative assumptions about the underlying stock-recruit function (Beverton-Holt and Ricker). We provide a posterior predictive distribution of steepness for rockfishes that can be used as a prior in future assessments for similar populations. We also evaluate whether F
40%
is an appropriate proxy for F
MSY
for Pacific rockfishes and show that uncertainty in the natural mortality rate can have a significant effect on management advice derived from meta-analyses of stock-recruit data.