Smith, James A., Mercedes Pozo Buil, Barbara A Muhling, Desiree Tommasi, Stephanie Brodie, Timothy H Frawley, Jerome Fiechter, Stefan Koenigstein, Amber Himes-Cornell, Michael A Alexander, Steven J Bograd, Nathalí Cordero-Quirós, Larry B Crowder, Enrique N Curchitser, Stephanie J Green, Natasha A Hardy, Alan C Haynie, Elliot L Hazen, Kirstin Holsman, Gwendal Le Fol, Nerea Lezama-Ochoa, Ryan R Rykaczewski, Charles A Stock, Stephen Stohs, Jonathan Sweeney, Heather Welch, and Michael G Jacox, February 2023: Projecting climate change impacts from physics to fisheries: A view from three California Current fisheries. Progress in Oceanography, 211, 102973, DOI:10.1016/j.pocean.2023.102973. Abstract
Motivated by a need for climate-informed living marine resource management, increased emphasis has been placed on regional end-to-end modeling frameworks designed to project climate impacts on marine ecosystems and evaluate the efficacy of potential management strategies under changing conditions. The ‘Future Seas’ project was initiated with a focus on three fisheries (Pacific sardine, swordfish, and albacore tuna) in the California Current System (CCS). This work leverages a suite of climate, ocean, ecosystem, and economic models to project physical, ecological, and socio-economic change, evaluate management strategies, and quantify uncertainty in model projections. Here we describe the components of the modeling framework, considerations underlying choices made in model development, engagement with stakeholders, and key physical, ecological, and socio-economic results to date, including projections to 2100. Our broad aims are to (i) synthesize a large body of climate and fisheries research that has been conducted, and continues, under the Future Seas umbrella, and (ii) provide insight and recommendations to those pursuing similar efforts for other applications and in other regions. In general, our results indicate that all three species will likely shift their distributions (predominantly poleward) in the future, which impacts accessibility to fishing fleets, spatial management, and quota allocation. For similar integrative climate-to-fisheries projections, we recommend attention is given to: recognizing potential biases arising from differences between the climate products used for ecological model fitting and those used for model projection; how sources of projection uncertainty are prioritized, incorporated, and communicated; and quantitatively linking scenarios – especially socio-economic scenarios – with climate and ecological projections.
Taboada, Fernando G., Jong-Yeon Park, Barbara A Muhling, Desiree Tommasi, Kisei R Tanaka, Ryan R Rykaczewski, Charles A Stock, and Jorge L Sarmiento, March 2023: Anticipating fluctuations of bigeye tuna in the Pacific Ocean from three-dimensional ocean biogeochemistry. Journal of Applied Ecology, 60(3), DOI:10.1111/1365-2664.14346463-479. Abstract
1) Subseasonal to decadal ocean forecasting can make significant contributions to achieving effective management of living marine resources in a changing ocean. Most applications rely on indirect proxies, however, often measured at the ocean surface and lacking a direct mechanistic link to the dynamics of marine populations.
2) Here, we take advantage of three-dimensional, dynamical reconstructions and forecasts of ocean biogeochemistry based on a global Earth system model to hindcast and assess the capacity to anticipate fluctuations in the dynamics of bigeye tuna (Thunnus obesus Lowe) in the Pacific Ocean during the last six decades. We reconstructed spatial patterns in catch per unit effort (CPUE) through the combination of physiological indices capturing both habitat preferences and physiological tolerance limits in bigeye tuna.
3) Our analyses revealed a sequence of four distinct regimes characterized by changes in the zonal distribution and average CPUE of bigeye tuna in the Pacific Ocean. Habitat models accounting for basin-wide fluctuations in the thermal structure and oxygen concentration throughout the water column captured interannual fluctuations in CPUE and regime switches that models based solely on surface information were unable to reproduce. Decade-long forecast experiments further suggested that forecasts of three-dimensional biogeochemical information might enable anticipation of fluctuations in bigeye tuna several years ahead.
4) Synthesis and applications. Together, our results reveal the impact of variability of biogeochemical conditions in the ocean interior on the dynamics of bigeye tuna on the Pacific Ocean, raising concerns about the future impact of ocean warming and deoxygenation. The results also lend support to incorporating subsurface biogeochemical information into ecological forecasts to implement efficient dynamic management strategies and promote the sustainable use of marine living resources.
Efforts to manage living marine resources (LMRs) under climate change need projections of future ocean conditions, yet most global climate models (GCMs) poorly represent critical coastal habitats. GCM utility for LMR applications will increase with higher spatial resolution but obstacles including computational and data storage costs, obstinate regional biases, and formulations prioritizing global robustness over regional skill will persist. Downscaling can help address GCM limitations, but significant improvements are needed to robustly support LMR science and management. We synthesize past ocean downscaling efforts to suggest a protocol to achieve this goal. The protocol emphasizes LMR-driven design to ensure delivery of decision-relevant information. It prioritizes ensembles of downscaled projections spanning the range of ocean futures with durations long enough to capture climate change signals. This demands judicious resolution refinement, with pragmatic consideration for LMR-essential ocean features superseding theoretical investigation. Statistical downscaling can complement dynamical approaches in building these ensembles. Inconsistent use of bias correction indicates a need for objective best practices. Application of the suggested protocol should yield regional ocean projections that, with effective dissemination and translation to decision-relevant analytics, can robustly support LMR science and management under climate change.
Hervieux, G, Michael A Alexander, Charles A Stock, Michael G Jacox, Kathleen Pegion, E Becker, Frederic Castruccio, and Desiree Tommasi, December 2019: More reliable coastal SST forecasts from the North American multimodel ensemble. Climate Dynamics, 53(12), DOI:10.1007/s00382-017-3652-7. Abstract
The skill of monthly sea surface temperature (SST) anomaly predictions for large marine ecosystems (LMEs) in coastal regions of the United States and Canada is assessed using simulations from the climate models in the North American Multimodel Ensemble (NMME). The forecasts based on the full ensemble are generally more skillful than predictions from even the best single model. The improvement in skill is particularly noteworthy for probability forecasts that categorize SST anomalies into upper (warm) and lower (cold) terciles. The ensemble provides a better estimate of the full range of forecast values than any individual model, thereby correcting for the systematic over-confidence (under-dispersion) of predictions from an individual model. Probability forecasts, including tercile predictions from the NMME, are used frequently in seasonal forecasts for atmospheric variables and may have many uses in marine resource management.
Jacox, Michael G., Desiree Tommasi, Michael A Alexander, G Hervieux, and Charles A Stock, July 2019: Predicting the evolution of the 2014-16 California Current System marine heatwave from an ensemble of coupled global climate forecasts. Frontiers in Marine Science, 6(497), DOI:10.3389/fmars.2019.00497. Abstract
Throughout 2014-2016, the California Current System (CCS) was characterized by large and persistent sea surface temperature anomalies (SSTa), which were accompanied by widespread ecological and socioeconomic consequences that have been documented extensively in the scientific literature and in the popular press. This marine heatwave and others have resulted in a heightened awareness of their potential impacts and prompted questions about if and when they may be predictable. Here, we use output from an ensemble of global climate forecast systems to document which aspects of the 2014-2016 CCS heatwave were predictable and how forecast skill, or lack thereof, relates to mechanisms driving the heatwave’s evolution. We focus on four prominent SSTa changes within the 2014-2016 period: (i) the initial onset of anomalous warming in early 2014, (ii) a second rapid SSTa increase in late 2014, (iii) a sharp reduction and subsequent return of warm SSTa in mid-2015, and (iv) another anomalous warming event in early 2016. Models exhibited clear forecast skill for the first and last of these fluctuations, but not the two in the middle. Taken together with the state of knowledge on the dominant forcing mechanisms of this heatwave, our results suggest that CCS SSTa forecast skill derives from predictable evolution of pre-existing SSTa to the west (as in early 2014) and the south (as in early 2016), while the inability of models to forecast wind-driven SSTa in late 2014 and mid-2015 is consistent with the lack of a moderate or strong El Niño or La Niña event preceding those periods. The multi-model mean forecast consistently outperformed a damped persistence forecast, especially during the period of largest SSTa, and skillful CCS forecasts were generally associated with accurate representation of large-scale dynamics. Additionally, a large forecast ensemble (85 members) indicated elevated probabilities for observed SSTa extremes even when ensemble mean forecasts exhibited limited skill. Our results suggest that different types or aspects of marine heat waves are more or less predictable depending on the forcing mechanisms at play, and events that are consistent with predictable ocean responses could inform ecosystem-based management of the ocean.
Muhling, Barbara A., Carlos F Gaitán, Charles A Stock, Vincent S Saba, Desiree Tommasi, and Keith W Dixon, March 2018: Potential Salinity and Temperature Futures for the Chesapeake Bay Using a Statistical Downscaling Spatial Disaggregation Framework. Estuaries and Coasts, 41(2), DOI:10.1007/s12237-017-0280-8. Abstract
Estuaries are productive and ecologically important ecosystems, incorporating environmental drivers from watersheds, rivers, and the coastal ocean. Climate change has potential to modify the physical properties of estuaries, with impacts on resident organisms. However, projections from general circulation models (GCMs) are generally too coarse to resolve important estuarine processes. Here, we statistically downscaled near-surface air temperature and precipitation projections to the scale of the Chesapeake Bay watershed and estuary. These variables were linked to Susquehanna River streamflow using a water balance model and finally to spatially resolved Chesapeake Bay surface temperature and salinity using statistical model trees. The low computational cost of this approach allowed rapid assessment of projected changes from four GCMs spanning a range of potential futures under a high CO2 emission scenario, for four different downscaling methods. Choice of GCM contributed strongly to the spread in projections, but choice of downscaling method was also influential in the warmest models. Models projected a ~2–5.5 °C increase in surface water temperatures in the Chesapeake Bay by the end of the century. Projections of salinity were more uncertain and spatially complex. Models showing increases in winter-spring streamflow generated freshening in the Upper Bay and tributaries, while models with decreased streamflow produced salinity increases. Changes to the Chesapeake Bay environment have implications for fish and invertebrate habitats, as well as migration, spawning phenology, recruitment, and occurrence of pathogens. Our results underline a potentially expanded role of statistical downscaling to complement dynamical approaches in assessing climate change impacts in dynamically challenging estuaries.
Morse, R E., K D Friedland, Desiree Tommasi, Charles A Stock, and Janet A Nye, January 2017: Distinct zooplankton regime shift patterns across ecoregions of the U.S. Northeast continental shelf Large Marine Ecosystem. Journal of Marine Systems, 165, DOI:10.1016/j.jmarsys.2016.09.011. Abstract
We investigated regime shifts in seasonal zooplankton communities of the Northeast continental shelf Large Marine Ecosystem (NES) and its subcomponent ecoregions over a multi-decadal period (1977–2013). Our cross ecoregion analysis shows that regime shifts in different ecoregions often exhibited very distinct characteristics, emphasizing more granular fluctuations in NES plankton communities relative to previous work. Shifts early in the time series generally reflected an increase in abundance levels. The response of zooplankton abundance within fall communities was more similar among ecoregions than for spring communities. The Gulf of Maine exhibited highly distinct patterns from other ecoregions, with regime shifts identified in the early 1980s, early 2000s, and mid-2000s for spring communities. Regime shifts were identified in the early to mid-1990s for the NES, Georges Bank, and the Mid-Atlantic Bight ecoregions, while the fall communities experienced shifts in the early 1990s and late 1980s for the NES and Georges Bank, but in the late 1990s in the Mid-Atlantic Bight. A constrained correspondence analysis of zooplankton community against local and basin-scale climatological indices suggests that water temperature, stratification, and the Atlantic multidecadal oscillation (AMO) were the predominant factors in driving the zooplankton community composition.
Tommasi, Desiree, Charles A Stock, A J Hobday, R Methot, Isaac C Kaplan, J P Eveson, Kirstin Holsman, Timothy J Miller, Sarah K Gaichas, Marion Gehlen, A Pershing, Gabriel A Vecchi, Rym Msadek, Thomas L Delworth, C M Eakin, Melissa A Haltuch, Roland Séférian, C M Spillman, J R Hartog, Samantha A Siedlecki, Jameal F Samhouri, Barbara A Muhling, R G Asch, M Pinsky, Vincent S Saba, Sarah B Kapnick, and Carlos F Gaitán, et al., March 2017: Managing living marine resources in a dynamic environment: The role of seasonal to decadal climate forecasts. Progress in Oceanography, 152, DOI:10.1016/j.pocean.2016.12.011. Abstract
Recent developments in global dynamical climate prediction systems have allowed for skillful predictions of climate variables relevant to living marine resources (LMRs) at a scale useful to understanding and managing LMRs. Such predictions present opportunities for improved LMR management and industry operations, as well as new research avenues in fisheries science. LMRs respond to climate variability via changes in physiology and behavior. For species and systems where climate-fisheries links are well established, forecasted LMR responses can lead to anticipatory and more effective decisions, benefitting both managers and stakeholders. Here, we provide an overview of climate prediction systems and advances in seasonal to decadal prediction of marine-resource relevant environmental variables. We then describe a range of climate-sensitive LMR decisions that can be taken at lead-times of months to decades, before highlighting a range of pioneering case studies using climate predictions to inform LMR decisions. The success of these case studies suggests that many additional applications are possible. Progress, however, is limited by observational and modeling challenges. Priority developments include strengthening of the mechanistic linkages between climate and marine resource responses, development of LMR models able to explicitly represent such responses, integration of climate driven LMR dynamics in the multi-driver context within which marine resources exist, and improved prediction of ecosystem-relevant variables at the fine regional scales at which most marine resource decisions are made. While there are fundamental limits to predictability, continued advances in these areas have considerable potential to make LMR managers and industry decision more resilient to climate variability and help sustain valuable resources. Concerted dialog between scientists, LMR managers and industry is essential to realizing this potential.
Tommasi, Desiree, Charles A Stock, Kathleen Pegion, and Gabriel A Vecchi, et al., March 2017: Improved management of small pelagic fisheries through seasonal climate prediction. Ecological Applications, 27(2), DOI:10.1002/eap.1458. Abstract
Populations of small pelagic fish are strongly influenced by climate. The inability of managers to anticipate environment-driven fluctuations in stock productivity or distribution can lead to overfishing and stock collapses, inflexible management regulations inducing shifts in the functional response to human predators, lost opportunities to harvest populations, bankruptcies in the fishing industry, and loss of resilience in the human food supply. Recent advances in dynamical global climate prediction systems allow for sea surface temperature (SST) anomaly predictions at a seasonal scale over many shelf ecosystems. Here we assess the utility of SST predictions at this “fishery relevant” scale to inform management, using Pacific sardine as a case study. The value of SST anomaly predictions to management was quantified under four harvest guidelines (HGs) differing in their level of integration of SST data and predictions. The HG that incorporated stock biomass forecasts informed by skillful SST predictions led to increases in stock biomass and yield, and reductions in the probability of yield and biomass falling below socioeconomic or ecologically acceptable levels. However, to mitigate the risk of collapse in the event of an erroneous forecast, it was important to combine such forecast-informed harvest controls with additional harvest restrictions at low biomass.
Decisions made by fishers and fisheries managers are informed by climate and fisheries observations that now often span more than 50 years. Multi-annual climate forecasts could further inform such decisions if they were skillful in predicting future conditions relative to the 50-year scope of past variability. We demonstrate that an existing multi-annual prediction system skillfully forecasts the probability of next year, the next 1–3 years, and the next 1–10 years being warmer or cooler than the 50-year average at the surface in coastal ecosystems. Probabilistic forecasts of upper and lower seas surface temperature (SST) terciles over the next 3 or 10 years from the GFDL CM 2.1 10-member ensemble global prediction system showed significant improvements in skill over the use of a 50-year climatology for most Large Marine Ecosystems (LMEs) in the North Atlantic, the western Pacific, and Indian oceans. Through a comparison of the forecast skill of initialized and uninitialized hindcasts, we demonstrate that this skill is largely due to the predictable signature of radiative forcing changes over the 50-year timescale rather than prediction of evolving modes of climate variability. North Atlantic LMEs stood out as the only coastal regions where initialization significantly contributed to SST prediction skill at the 1 to 10 year scale.
Friedland, K D., R T Leaf, J Kane, Desiree Tommasi, R G Asch, N Rebuck, R Ji, S I Large, Charles A Stock, and Vincent S Saba, July 2015: Spring bloom dynamics and zooplankton biomass response on the US Northeast Continental Shelf. Continental Shelf Research, 102, DOI:10.1016/j.csr.2015.04.005. Abstract
The spring phytoplankton bloom on the US Northeast Continental Shelf is a feature of the ecosystem production cycle that varies annually in timing, spatial extent, and magnitude. To quantify this variability, we analyzed remotely-sensed ocean color data at two spatial scales, one based on ecologically defined sub-units of the ecosystem (production units) and the other on a regular grid (0.5°). Five units were defined: Gulf of Maine East and West, Georges Bank, and Middle Atlantic Bight North and South. The units averaged 47×103 km2 in size. The initiation and termination of the spring bloom were determined using change-point analysis with constraints on what was identified as a bloom based on climatological bloom patterns. A discrete spring bloom was detected in most years over much of the western Gulf of Maine production unit. However, bloom frequency declined in the eastern Gulf of Maine and transitioned to frequencies as low as 50% along the southern flank of the Georges Bank production unit. Detectable spring blooms were episodic in the Middle Atlantic Bight production units. In the western Gulf of Maine, bloom duration was inversely related to bloom start day; thus, early blooms tended to be longer lasting and larger magnitude blooms. We view this as a phenological mismatch between bloom timing and the “top-down” grazing pressure that terminates a bloom. Estimates of secondary production were available from plankton surveys that provided spring indices of zooplankton biovolume. Winter chlorophyll biomass had little effect on spring zooplankton biovolume, whereas spring chlorophyll biomass had mixed effects on biovolume. There was evidence of a “bottom up” response seen on Georges Bank where spring zooplankton biovolume was positively correlated with the concentration of chlorophyll. However, in the western Gulf of Maine, biovolume was uncorrelated with chlorophyll concentration, but was positively correlated with bloom start and negatively correlated with magnitude. This observation is consistent with both a “top-down” mechanism of control of the bloom and a “bottom-up” effect of bloom timing on zooplankton grazing. Our inability to form a consistent model of these relationships across adjacent systems underscores the need for further research.
Kearney, Kelly A., Desiree Tommasi, and Charles A Stock, September 2015: Simulated ecosystem response to volcanic iron fertilization in the subarctic Pacific ocean. Fisheries Oceanography, 24(5), DOI:10.1111/fog.12118. Abstract
The eruption of the Kasatochi volcano in August 2008 stimulated an anomalously high phytoplankton bloom in the otherwise iron-limited subarctic Pacific ocean. It has been proposed that this increased production may have been responsible for record returns of some Pacific salmon stocks in the following years. Here, we investigate the potential effect of volcanic-induced iron fertilization on the entire ecosystem, from phytoplankton through to top predators, using a fully-coupled end-to-end ecosystem model. Our simulations indicate that the volcanic iron fertilization could only stimulate modest increases, at most 10%, in the standing stock biomass of upper trophic level species, including fisheries targets such as Pacific salmon. Propagation of energy to higher trophic levels depends on the timing of the eruption, with more efficient crustaceous zooplankton pathways being favored earlier in the growing season and less-efficient gelatinous zooplankton pathways dominating during later months. However, effects were of modest magnitude for all eruption timings, and the strong level of connectivity within the food web makes the preferential stimulation of a single salmon stock implausible. This adds additional support to evidence suggesting that the Kasatochi eruption did not play a large role in subsequent high salmon returns and questions the value of much smaller-scale artificial fertilization for fisheries. Indeed, the onset of macronutrient limitation coupled with the highly-connected nature of the food web exert strong controls on the fisheries response to even complete removal of iron limitation in the subarctic Pacific.
Sea surface temperature (SST) anomalies are often both leading indicators and important drivers of marine resource fluctuations. Assessment of the skill of SST anomaly forecasts within coastal ecosystems accounting for the majority of global fish yields, however, has been minimal. This reflects coarse global forecast system resolution and past emphasis on the predictability of ocean basin-scale SST variations. This paper assesses monthly to inter-annual SST anomaly predictions in coastal “Large Marine Ecosystems” (LMEs). We begin with an analysis of 7 well-observed LMEs adjacent to the United States and then examine how mechanisms responsible for prediction skill in these systems are reflected in predictions for LMEs globally. Historical SST anomaly estimates from the 1/4o daily Optimal Interpolation Sea Surface Temperature reanalysis (OISST.v2) were first found to be highly consistent with in-situ measurements for 6 of the 7 U.S. LMEs. Thirty years of retrospective forecasts from climate forecast systems developed at NOAA’s Geophysical Fluid Dynamics Laboratory (CM2.5-FLOR) and the National Center for Environmental Prediction (CFSv2) were then assessed against OISST.v2. Forecast skill varied widely by LME, initialization month, and lead but there were many cases of high skill that also exceeded that of a persistence forecast, some at leads greater than 6 months. Mechanisms underlying skill above persistence included accurate simulation of a) seasonal transitions between less predictable locally generated and more predictable basin-scale SST variability; b) seasonal transitions between different basin-scale influences; c) propagation of SST anomalies across seasons through sea ice; and d) re-emergence of previous anomalies upon the breakdown of summer stratification. Globally, significant skill above persistence across many tropical systems arises via mechanisms a) and b). Combinations of all four mechanisms contribute to less prevalent but nonetheless significant skill in extratropical systems. While continued refinement of global climate forecast systems and observations are needed to improve coastal SST anomaly prediction and extend predictions to other ecosystem relevant variables (e.g., salinity), present skill warrants close examination of forecasts for marine resource applications.
Tommasi, Desiree, Janet A Nye, and Charles A Stock, et al., July 2015: Effect of Environmental Conditions on Juvenile Recruitment of Alewife (Alosa pseudoharengus) and Blueback Herring (A. aestivalis) in Freshwater: A Coastwide Perspective. Canadian Journal of Fisheries and Aquatic Sciences, 72(7), DOI:10.1139/cjfas-2014-0259. Abstract
The abundance of alewife (Alosa pseudoharengus) and blueback herring (Alosa aestivalis) has declined throughout their range, and there are increasing concerns about their conservation status. Because of their diadromous life history, variability in rates of survival in freshwater can affect overall recruitment. The objective of our study was to assess how river temperature and flow influence young of the year (YOY) river herring recruitment in the Northeast US. Observations of adult and juvenile fish in five rivers were used to construct spawner-YOY recruits models; these rivers were chosen because of the length of the time series (>15 years) and the paired observations of spawners and juveniles. An environmentally-explicit stock recruitment model explained a substantial fraction (41 to 80%) of the variance in YOY abundance, depending on river system. Our approach allowed for a preliminary discussion of potential mechanisms, which need to be further substantiated by focused field and laboratory studies. Early summer river flow and river temperature had the greatest influence indicating the importance of conditions in nursery habitats. In certain systems, spring or fall conditions were also important determinants of survival suggesting additional effects of the environment on spawning of adults and juvenile egress from freshwater nursery habitats.