Species redistribution creates unequal outcomes for multispecies fisheries under projected climate change.

Climate change drives species distribution shifts, affecting the availability of resources people rely upon for food and livelihoods. These impacts are complex, manifest at local scales, and have diverse effects across multiple species. However, for wild capture fisheries, current understanding is dominated by predictions for individual species at coarse spatial scales. We show that species-specific responses to localized environmental changes will alter the collection of co-occurring species within established fishing footprints along the U.S. West Coast. We demonstrate that availability of the most economically valuable, primary target species is highly likely to decline coastwide in response to warming and reduced oxygen concentrations, while availability of the most abundant, secondary target species will potentially increase. A spatial reshuffling of primary and secondary target species suggests regionally heterogeneous opportunities for fishers to adapt by changing where or what they fish. Developing foresight into the collective responses of species at local scales will enable more effective and tangible adaptation pathways for fishing communities.

Recommendations for quantifying and reducing uncertainty in climate projections of species distributions.

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.

Climate vulnerability assessment for Pacific salmon and steelhead in the California Current Large Marine Ecosystem.

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.

Competing tradeoffs between increasing marine mammal predation and fisheries harvest of Chinook salmon.

Many marine mammal predators, particularly pinnipeds, have increased in abundance in recent decades, generating new challenges for balancing human uses with recovery goals via ecosystem-based management. We used a spatio-temporal bioenergetics model of the Northeast Pacific Ocean to quantify how predation by three species of pinnipeds and killer whales (Orcinus orca) on Chinook salmon (Oncorhynchus tshawytscha) has changed since the 1970s along the west coast of North America, and compare these estimates to salmon fisheries. We find that from 1975 to 2015, biomass of Chinook salmon consumed by pinnipeds and killer whales increased from 6,100 to 15,200 metric tons (from 5 to 31.5 million individual salmon). Though there is variation across the regions in our model, overall, killer whales consume the largest biomass of Chinook salmon, but harbor seals (Phoca vitulina) consume the largest number of individuals. The decrease in adult Chinook salmon harvest from 1975-2015 was 16,400 to 9,600 metric tons. Thus, Chinook salmon removals (harvest + consumption) increased in the past 40 years despite catch reductions by fisheries, due to consumption by recovering pinnipeds and endangered killer whales. Long-term management strategies for Chinook salmon will need to consider potential conflicts between rebounding predators or endangered predators and prey.

Risks of ocean acidification in the California Current food web and fisheries: ecosystem model projections.

The benefits and ecosystem services that humans derive from the oceans are threatened by numerous global change stressors, one of which is ocean acidification. Here, we describe the effects of ocean acidification on an upwelling system that already experiences inherently low pH conditions, the California Current. We used an end-to-end ecosystem model (Atlantis), forced by downscaled global climate models and informed by a meta-analysis of the pH sensitivities of local taxa, to investigate the direct and indirect effects of future pH on biomass and fisheries revenues. Our model projects a 0.2-unit drop in pH during the summer upwelling season from 2013 to 2063, which results in wide-ranging magnitudes of effects across guilds and functional groups. The most dramatic direct effects of future pH may be expected on epibenthic invertebrates (crabs, shrimps, benthic grazers, benthic detritivores, bivalves), and strong indirect effects expected on some demersal fish, sharks, and epibenthic invertebrates (Dungeness crab) because they consume species known to be sensitive to changing pH. The model's pelagic community, including marine mammals and seabirds, was much less influenced by future pH. Some functional groups were less affected to changing pH in the model than might be expected from experimental studies in the empirical literature due to high population productivity (e.g., copepods, pteropods). Model results suggest strong effects of reduced pH on nearshore state-managed invertebrate fisheries, but modest effects on the groundfish fishery because individual groundfish species exhibited diverse responses to changing pH. Our results provide a set of projections that generally support and build upon previous findings and set the stage for hypotheses to guide future modeling and experimental analysis on the effects of OA on marine ecosystems and fisheries.

Extending Vulnerability Assessment to Include Life Stages Considerations.

Species are experiencing a suite of novel stressors from anthropogenic activities that have impacts at multiple scales. Vulnerability assessment is one tool to evaluate the likely impacts that these stressors pose to species so that high-vulnerability cases can be identified and prioritized for monitoring, protection, or mitigation. Commonly used semi-quantitative methods lack a framework to explicitly account for differences in exposure to stressors and organism responses across life stages. Here we propose a modification to commonly used spatial vulnerability assessment methods that includes such an approach, using ocean acidification in the California Current as an illustrative case study. Life stage considerations were included by assessing vulnerability of each life stage to ocean acidification and were used to estimate population vulnerability in two ways. We set population vulnerability equal to: (1) the maximum stage vulnerability and (2) a weighted mean across all stages, with weights calculated using Lefkovitch matrix models. Vulnerability was found to vary across life stages for the six species explored in this case study: two krill-Euphausia pacifica and Thysanoessa spinifera, pteropod-Limacina helicina, pink shrimp-Pandalus jordani, Dungeness crab-Metacarcinus magister and Pacific hake-Merluccius productus. The maximum vulnerability estimates ranged from larval to subadult and adult stages with no consistent stage having maximum vulnerability across species. Similarly, integrated vulnerability metrics varied greatly across species. A comparison showed that some species had vulnerabilities that were similar between the two metrics, while other species' vulnerabilities varied substantially between the two metrics. These differences primarily resulted from cases where the most vulnerable stage had a low relative weight. We compare these methods and explore circumstances where each method may be appropriate.

Experiments with Seasonal Forecasts of ocean conditions for the Northern region of the California Current upwelling system.

Resource managers at the state, federal, and tribal levels make decisions on a weekly to quarterly basis, and fishers operate on a similar timeframe. To determine the potential of a support tool for these efforts, a seasonal forecast system is experimented with here. JISAO's Seasonal Coastal Ocean Prediction of the Ecosystem (J-SCOPE) features dynamical downscaling of regional ocean conditions in Washington and Oregon waters using a combination of a high-resolution regional model with biogeochemistry and forecasts from NOAA's Climate Forecast System (CFS). Model performance and predictability were examined for sea surface temperature (SST), bottom temperature, bottom oxygen, pH, and aragonite saturation state through model hindcasts, reforecast, and forecast comparisons with observations. Results indicate J-SCOPE forecasts have measurable skill on seasonal timescales. Experiments suggest that seasonal forecasting of ocean conditions important for fisheries is possible with the right combination of components. Those components include regional predictability on seasonal timescales of the physical environment from a large-scale model, a high-resolution regional model with biogeochemistry that simulates seasonal conditions in hindcasts, a relationship with local stakeholders, and a real-time observational network. Multiple efforts and approaches in different regions would advance knowledge to provide additional tools to fishers and other stakeholders.

An Integrated Coral Reef Ecosystem Model to Support Resource Management under a Changing Climate.

Millions of people rely on the ecosystem services provided by coral reefs, but sustaining these benefits requires an understanding of how reefs and their biotic communities are affected by local human-induced disturbances and global climate change. Ecosystem-based management that explicitly considers the indirect and cumulative effects of multiple disturbances has been recommended and adopted in policies in many places around the globe. Ecosystem models give insight into complex reef dynamics and their responses to multiple disturbances and are useful tools to support planning and implementation of ecosystem-based management. We adapted the Atlantis Ecosystem Model to incorporate key dynamics for a coral reef ecosystem around Guam in the tropical western Pacific. We used this model to quantify the effects of predicted climate and ocean changes and current levels of current land-based sources of pollution (LBSP) and fishing. We used the following six ecosystem metrics as indicators of ecosystem state, resilience and harvest potential: 1) ratio of calcifying to non-calcifying benthic groups, 2) trophic level of the community, 3) biomass of apex predators, 4) biomass of herbivorous fishes, 5) total biomass of living groups and 6) the end-to-start ratio of exploited fish groups. Simulation tests of the effects of each of the three drivers separately suggest that by mid-century climate change will have the largest overall effect on this suite of ecosystem metrics due to substantial negative effects on coral cover. The effects of fishing were also important, negatively influencing five out of the six metrics. Moreover, LBSP exacerbates this effect for all metrics but not quite as badly as would be expected under additive assumptions, although the magnitude of the effects of LBSP are sensitive to uncertainty associated with primary productivity. Over longer time spans (i.e., 65 year simulations), climate change impacts have a slight positive interaction with other drivers, generally meaning that declines in ecosystem metrics are not as steep as the sum of individual effects of the drivers. These analyses offer one way to quantify impacts and interactions of particular stressors in an ecosystem context and so provide guidance to managers. For example, the model showed that improving water quality, rather than prohibiting fishing, extended the timescales over which corals can maintain high abundance by at least 5-8 years. This result, in turn, provides more scope for corals to adapt or for resilient species to become established and for local and global management efforts to reduce or reverse stressors.

Indirect effects of conservation policies on the coupled human-natural ecosystem of the upper Gulf of California.

High bycatch of non-target species and species of conservation concern often drives the implementation of fisheries policies. However, species- or fishery-specific policies may lead to indirect consequences, positive or negative, for other species or fisheries. We use an Atlantis ecosystem model of the Northern Gulf of California to evaluate the effects of fisheries policies directed at reducing bycatch of vaquita (Phocoena sinus) on other species of conservation concern, priority target species, and metrics of ecosystem function and structure. Vaquita, a Critically Endangered porpoise endemic to the Upper Gulf of California, are frequently entangled by finfish gillnets and shrimp driftnets. We tested five fishery management scenarios, projected over 30 years (2008 to 2038), directed at vaquita conservation. The scenarios consider progressively larger spatial restrictions for finfish gillnets and shrimp driftnets. The most restrictive scenario resulted in the highest biomass of species of conservation concern; the scenario without any conservation measures in place resulted in the lowest. Vaquita experienced the largest population increase of any functional group; their biomass increased 2.7 times relative to initial (2008) levels under the most restrictive spatial closure scenario. Bycatch of sea lions, sea turtles, and totoaba decreased > 80% in shrimp driftnets and at least 20% in finfish gillnet fleets under spatial management. We found indirect effects on species and ecosystem function and structure as a result of vaquita management actions. Biomass and catch of forage fish declined, which could affect lower-trophic level fisheries, while other species such as skates, rays, and sharks increased in both biomass and catch. When comparing across performance metrics, we found that scenarios that increased ecosystem function and structure resulted in lower economic performance indicators, underscoring the need for management actions that consider ecological and economic tradeoffs as part of the integrated management of the Upper Gulf of California.

Exploring trade-offs between fisheries and conservation of the vaquita porpoise (Phocoena sinus) using an Atlantis ecosystem model.

BACKGROUND: Minimizing fishery bycatch threats might involve trade-offs between maintaining viable populations and economic benefits. Understanding these trade-offs can help managers reconcile conflicting goals. An example is a set of bycatch reduction measures for the Critically Endangered vaquita porpoise (Phocoena sinus), in the Northern Gulf of California, Mexico. The vaquita is an endemic species threatened with extinction by artisanal net bycatch within its limited range; in this area fisheries are the chief source of economic productivity. METHODOLOGY/PRINCIPAL FINDINGS: We analyze trade-offs between conservation of the vaquita and fisheries, using an end-to-end Atlantis ecosystem model for the Northern Gulf of California. Atlantis is a spatially-explicit model intended as a strategic tool to test alternative management strategies. We simulated increasingly restrictive fisheries regulations contained in the vaquita conservation plan: implementing progressively larger spatial management areas that exclude gillnets, shrimp driftnets and introduce a fishing gear that has no vaquita bycatch. We found that only the most extensive spatial management scenarios recovered the vaquita population above the threshold necessary to downlist the species from Critically Endangered. The scenario that excludes existing net gear from the 2008 area of vaquita distribution led to moderate decrease in net present value (US$ 42 million) relative to the best-performing scenario and a two-fold increase in the abundance of adult vaquita over the course of 30 years. CONCLUSIONS/SIGNIFICANCE: Extended spatial management resulted in the highest recovery of the vaquita population. The economic cost of proposed management actions was unequally divided between fishing fleets; the loss of value from finfish gillnet fisheries was never recovered. Our analysis shows that managers will have to confront difficult trade-offs between management scenarios for vaquita conservation.

Impacts of fishing low-trophic level species on marine ecosystems.

Low-trophic level species account for more than 30% of global fisheries production and contribute substantially to global food security. We used a range of ecosystem models to explore the effects of fishing low-trophic level species on marine ecosystems, including marine mammals and seabirds, and on other commercially important species. In five well-studied ecosystems, we found that fishing these species at conventional maximum sustainable yield (MSY) levels can have large impacts on other parts of the ecosystem, particularly when they constitute a high proportion of the biomass in the ecosystem or are highly connected in the food web. Halving exploitation rates would result in much lower impacts on marine ecosystems while still achieving 80% of MSY.

A statistical approach for estimating fish diet compositions from multiple, data sources: Gulf of California case study.

Trophic ecosystem models are one promising tool for providing ecosystem-based management advice. Diet and interaction rate parameters are critical in defining the behavior of these models, and will greatly influence any predictions made in response to management perturbations. However, most trophic ecosystem models must rely on a patchwork of data availability and must contend with knowledge gaps and poor quantification of uncertainty. Here we present a statistical method for combining diet information from field samples and literature to describe trophic relationships at the level of functional groups. In this example, original fieldwork in the northern Gulf of California, Mexico, provides gut content data for targeted and untargeted fish species. The field data are pooled with diet composition information from FishBase, an online data repository. Diet information is averaged across stomachs to represent an average predator, and then the data are bootstrapped to generate likelihood profiles. These are fit to a Dirichlet function, and from the resulting marginal distributions, maximum-likelihood estimates are generated with confidence intervals representing the likely contribution to diet for each predator-prey combination. We characterize trophic linkages into two broad feeding guilds, pelagic and demersal feeders, and explore differentiation within those guilds. We present an abbreviated food web for the northern Gulf of California based on the results of this study. This food web will form the basis of a trophic dynamic model. Compared to the common method of averaging diet compositions across predators, this statistical approach is less influenced by the presence of long tails in the distributions, which correspond to rare feeding events, and is therefore better suited to small data sets.