class: right, bottom, my-title, title-slide .title[ # Assessing small pelagic fish trends in space
and time using piscivore diet data ] .subtitle[ ## Bluefish Research Track Assessment Review
Term of Reference 1
7 December 2022 ] .author[ ### Sarah Gaichas
1
, James Gartland
2
, Brian Smith
1
, Elizabeth Ng
3
, Michael Celestino
4
,
Anthony Wood
1
, Katie Drew
5
, Abigail Tyrell
1, 6
, and James Thorson
7
] .institute[ ###
1
NOAA NMFS Northeast Fisheries Science Center, Woods Hole, MA, USA;
2
Virginia Institute of Marine Science, Gloucester Point, VA, USA;
3
University of Washington, Seattle, WA, USA;
4
New Jersey Department of Environmental Protection, Port Republic, NJ, USA;
5
Atlantic States Marine Fisheries Commission, Arlington, VA, USA;
6
Ocean Associates Inc, Arlington, VA, USA;
7
NOAA NMFS Alaska Fisheries Science Center, Seattle, WA, USA
] --- class: top, left # Does prey drive availability of bluefish? ## Bluefish, *Pomatomus saltatrix* .pull-left[ ![Bluefish illustration, credit NOAA Fisheries](https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Bluefish-NOAAFisheries.png) ] .pull-right[ .large[ "... it is perhaps the most ferocious and bloodthirsty fish in the sea, leaving in its wake a trail of dead and mangled mackerel, menhaden, herring, alewives, and other species on which it preys." <a name=cite-collette_bigelow_2002></a>([Collette, et al., 2002](#bib-collette_bigelow_2002)) "From Raritan Bay to Rockaway Inlet, we have had a phenomenal bluefish year with lots of bunker and other bait, ultimately leading to an abundance of bluefish." [Mid-Atlantic Bluefish Fishery Performance Report, 2021](https://www.mafmc.org/s/8_BF-FPR-2021.pdf) ] ] .center[ Can localized predator-prey observations scale to coastwide assessment and management? ] ??? <font size="-1"></font> Changing distribution and abundance of small pelagics may drive changes in predator distributions, affecting predator availability to fisheries and surveys. However, small pelagic fish are difficult to survey directly, so we developed a novel method of assessing small pelagic fish aggregate abundance via predator diet data. We used piscivore diet data collected from multiple bottom trawl surveys within a Vector Autoregressive Spatio-Temporal (VAST) model to assess trends of small pelagics on the Northeast US shelf. The goal was to develop a spatial “forage index” to inform survey and/or fishery availability in the bluefish (Pomatomus saltatrix) stock assessment. Using spring and fall surveys from 1973-2020, 20 small pelagic groups were identified as major bluefish prey using the diet data. Then, predators were grouped by diet similarity to identify 19 piscivore species with the most similar diet to bluefish in the region. Diets from all 20 piscivores were combined for the 20 prey groups at each surveyed location, and the total weight of small pelagic prey per predator stomach at each location was input into a Poisson-link delta model to estimate expected prey mass per predator stomach. Best fit models included spatial and spatio-temporal random effects, with predator mean length, number of predator species, and sea surface temperature as catchability covariates. Spring and fall prey indices were split into inshore and offshore areas to reflect changing prey availability over time in areas available to the recreational fishery and the bottom trawl survey, and also to contribute to regional ecosystem reporting --- ## Bluefish diet in the Northeast US: a mix of managed and unmanaged small pelagics .pull-left[ <img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Squid-Longfin-v1-NOAAFisheries.png" width="50%" /><img src="https://objects.liquidweb.services/images/201909/robert_aguilar,_serc_48650727166_3ce8edc47d_b.jpg" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Butterfish-NOAAFisheries.png" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Scup-NOAAFisheries_0.png" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Herring-Atlantic-NOAAFisheries.png" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Mackerel-Atlantic-NOAAFisheries.png" width="50%" /> ] .pull-right[ <img src="https://diveary.com/img/species/cache/323__jpg.jpeg" width="50%" /><img src="https://www.vims.edu/research/departments/fisheries/programs/juvenile_surveys/netnotes_listing/_photosets/0910_dec_jan/_northernsandlance.jpg" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Hake-Silver-NOAAFisheries.png" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2021-03/squid_illex_nb_w_0.png" width="50%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2021-01/640x427-Atlantic_Menhaden_NB_W.jpg" width="50%" /> ] .right[ Northeast Fisheries Science Center Diet Data Online: https://fwdp.shinyapps.io/tm2020/ ] ??? Using NEFSC bottom trawl survey diet data from 1973-2021, 20 small pelagic groups were identified as major bluefish prey with 10 or more observations (in descending order of observations): Longfin squids (*Doryteuthis* formerly *Loligo* sp.), Anchovy family (Engraulidae), bay anchovy (*Anchoa mitchilli*), Atlantic butterfish, (*Peprilus triachanthus*), Cephalopoda, (*Anchoa hepsetus*), red eye round herring (*Etrumeus teres*), Sandlance (*Ammodytes* sp.), scup (*Stenotomus chrysops*), silver hake (*Merluccius bilinearis*), shortfin squids (*Illex* sp.), Atlantic herring (*Clupea harengus*), Herring family (Clupeidae), Bluefish (*Pomatomus saltatrix*), silver anchovy (*Engraulis eurystole*), longfin inshore squid (*Doryteuthis pealeii*), Atlantic mackerel (*Scomber scombrus*), flatfish (Pleuronectiformes), weakfish (*Cynoscion regalis*), and Atlantic menhaden (*Brevoortia tyrannus*). Prey categories such as fish unidentified, Osteichthyes, and unidentified animal remains were not included in the prey list. Although unidentified fish and Osteichthyes can comprise a significant portion of bluefish stomach contents, we cannot assume that unidentified fish in other predator stomachs represent unidentified fish in bluefish stomachs. Image credits: Striped and bay anchovy photo--Robert Aguilar, Smithsonian Environmental Research Center; redeye round herring photo--https://diveary.com ; sandlance photo--Virginia Institute of Marine Science; all others NOAA Fisheries. --- background-image: url("https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/Ng2022dietvsassess.png"), url("https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/NgPredStomForageB2022.png") background-size: 800px, 400px background-position: right, left ## Fish stomach contents → Atlantic herring biomass estimates <a name=cite-ng_predator_2021></a>([Ng, et al., 2021](https://doi.org/10.1093/icesjms/fsab026)) <!--.pull-left-40[ ![:img Ng paper title](https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/NgPredStomForageB2022.png) ] .pull-right-60[]--> ??? --- ## Vector Autoregressive Spatio-Temporal (VAST) modeling <a name=cite-thorson_comparing_2017></a><a name=cite-thorson_guidance_2019></a>([Thorson, et al., 2017](https://doi.org/10.1093/icesjms/fsw193); [Thorson, 2019](http://www.sciencedirect.com/science/article/pii/S0165783618302820)) VAST models two linear predictors for an index: 1. encounter rate, and 2. positive catch (amount in stomach) A full model for the first linear predictor `\(\rho_1\)` for each observation `\(i\)` can include: * fixed intercepts `\(\beta_1\)` for each category `\(c\)` and time `\(t\)`, * spatial random effects `\(\omega_1\)` for each location `\(s\)` and category, * spatio-temporal random effects `\(\varepsilon_1\)` for each location, category, and time, * fixed vessel effects `\(\eta_1\)` by vessel `\(v\)` and category, and * fixed catchability impacts `\(\lambda_1\)` of covariates `\(Q\)` for each observation and variable `\(k\)`: `$$\rho_1(i) = \beta_1(c_i, t_i) + \omega_1^*(s_i, c_i) + \varepsilon_1^*(s_i, c_i, t_i) + \eta_1(v_i, c_i) + \sum_{k=1}^{n_k} \lambda_1(k) Q(i,k)$$` The full model for the second linear predictor `\(\rho_2\)` has the same structure, estimating `\(\beta_2\)`, `\(\omega_2\)`, `\(\varepsilon_2\)`, `\(\eta_2\)`, and `\(\lambda_2\)` using the observations, categories, locations, times, and covariates. We modeled aggregate small pelagic prey as a single category, and apply a Poisson-link delta model to estimate expected prey mass per predator stomach as in ([Ng, et al., 2021](https://doi.org/10.1093/icesjms/fsab026)). VAST model code and documentation: https://github.com/James-Thorson-NOAA/VAST ??? Spatial and spatio-temporal correlation decay with increasing distance estimated as `\(\kappa\)` in a Matern function with fixed smoothness and geometric anisotropy (directional correlation, optionally estimated by the model). Initial model selection consistently supported the inclusion of spatial and spatio-temporal random effects and anisotropy across all datasets: fall, spring, and annual. --- background-image: url("https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/BluefishOnlyFallData_by_year.png") background-size: 650px background-position: right top ## Bluefish stomachs only? .pull-left[ Due to uneven "sampling" of Atlantic herring by predators, ([Ng, et al., 2021](https://doi.org/10.1093/icesjms/fsab026)) recommended aggregating across predators to improve the diet-based Atlantic herring biomass index. Between 1985-2021 there were: .contrib[ * 25634 survey stations with diet collections. * 22751 survey stations with piscivore diets. + 9027 piscivore stations with bluefish prey; + 40% of piscivore stations have bluefish prey. * 1814 survey stations with bluefish diets. + 905 bluefish stations with bluefish prey; + 50% of bluefish stations have bluefish prey. ] For this index combining multiple small pelagic prey, aggregating across predators most similar to bluefish both increases sample size and reduces sampling variability due to different predator availability to surveys. ] .pull-right[ <!--![Fall bluefish NEFSC](https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/BluefishOnlyFallData_by_year.png) --> .center[ .footnote[ Bluefish diet collection stations, fall Northeast Fisheries Science Center surveys ] ] ] ??? The figure shows bluefish diet collection stations for fall surveys, 1985-2019. NEAMAP survey stations with diet collections for piscivores (n = 3838) had a higher proportion with our defined bluefish prey (n = 2418, 63.0015633%). --- ## Aggregating predators: diet similarity to bluefish .pull-left[ <img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/unnamed-chunk-3-1.png" width="504" /> ] .pull-right[ <img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Bluefish-NOAAFisheries.png" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-07/640x427-Cod-Atlantic-NOAAFisheries.png?itok=kNKcZ7iV" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Shark-SpinyDogfish-NOAAFisheries.png?itok=cdjTG3Hz" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Monkfish-NOAAFisheries.png?itok=7JAgAz-u" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Flounder-Summer-NOAAFisheries.png?itok=II2ii-Qw" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Pollock-Atlantic-NOAAFisheries.png?itok=ZFoDB-Qr" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Hake-Red-NOAAFisheries.png?itok=SyCYEmmm" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Halibut-Atlantic-right-NOAAFisheries.png?itok=uPvRdIBx" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-StripedBass-NOAAFisheries.png?itok=4ZQoQM0S" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2020-12/640x427-Hake_White_NB_W.jpg?itok=yHy_AuTT" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-09/640x427-Hake-Silver-NOAAFisheries.png" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/dam-migration/640x427-cusk.jpg?itok=u0fw0hiv" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2022-08/640x427-Squid-Longfin-v1-NOAAFisheries.png" width="25%" /><img src="https://media.fisheries.noaa.gov/styles/original/s3/2021-03/squid_illex_nb_w_0.png" width="25%" /><img src="https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/weakfishASMFC.png" width="25%" /><img src="https://photolib.noaa.gov/Portals/0//GravityImages/36881/ProportionalFixedWidth/sanc100698384x800x800.jpg" width="25%" /> <img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/unnamed-chunk-5-1.png" width="25%" /><img src="https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/depositphotos_7476236-stock-photo-zeus-faber-fish.jpg" width="25%" /><img src="https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/spanishmackerelASMFC.png" width="25%" /><img src="https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/spottedseatroutASMFC.png" width="25%" /> ] ??? All size classes of 50 fish predators captured in the NEFSC bottom trawl survey were grouped by diet similarity to identify the size classes of piscivore species with the most similar diet to bluefish in the region. Diet similarity analysis was completed using the Schoener similarity index (@schoener_nonsynchronous_1970; B. Smith, pers. comm.), and is available available via [this link on the NEFSC food habits shiny app](https://fwdp.shinyapps.io/tm2020/#4_DIET_OVERLAP_AND_TROPHIC_GUILDS). The working group evaluated several clustering methods to develop the predator list (see [this link with detailed cluster results](https://sgaichas.github.io/bluefishdiet/PreySimilarityUpdate.html)). Predators with highest diet similarity to Bluefish from the NEFSC diet database (1973-2020) include Atlantic cod, Atlantic halibut, buckler dory, cusk, fourspot flounder, goosefish, longfin squid, shortfin squid, pollock, red hake, sea raven, silver hake, spiny dogfish, spotted hake, striped bass, summer flounder, thorny skate, weakfish, and white hake. The NEAMAP survey operates closer to shore than the current NEFSC survey. The NEAMAP dataset includes predators sampled by the NEFSC survey and adds two species, Spanish mackerel and spotted sea trout, not captured by the NEFSC survey offshore but included based on working group expert judgement of prey similarity to bluefish. Predator size classes included are listed in Table 2 of the forage fish index working paper at [this link](https://sgaichas.github.io/bluefishdiet/VASTcovariates_forageindex_WP.html). Image credits: Weakfish and Spanish mackerel-- https://marinefishesofgeorgia.org ; spotted seatrout-- https://fishinginmiami.com ; Sea Raven photo 2/11/2019 11:07:56 AM, Photographer: Andrew J. Martinez, Location: Massachusetts, Stellwagen Bank NMS; all others NOAA Fisheries. --- background-image: url("https://github.com/NOAA-EDAB/presentations/raw/master/docs/EDAB_images/AllPiscivoresFallData_by_year.png") background-size: 600px background-position: right top .pull-left[ ## "Catchability" covariates for aggregate predator samplers at a location Number of predator species → likely to affect *encounter rate* Mean size of predators → likely to affect *amount of prey* ([Ng, et al., 2021](https://doi.org/10.1093/icesjms/fsab026)) Sea surface temperature (SST) → likely to affect predator activity and feeding rate *encounter rate and amount of prey* * Many missing SST measurements for surveys before 1991 * NOAA OI SST V2 High Resolution Dataset <a name=cite-reynolds_daily_2007></a>([Reynolds, et al., 2007](https://journals.ametsoc.org/view/journals/clim/20/22/2007jcli1824.1.xml)) filled gaps Model selection consistently included number of predator species, mean predator size, and SST as covariates using fall, spring, and annual datasets ] .pull-right[ ] .footnote[ All piscivore diet collection stations, fall Northeast Fisheries Science Center surveys → ] ??? Diets from all 22 piscivores (including bluefish) were combined for the 20 forage fish (bluefish prey) groups at each surveyed location, and the mean weight of forage fish per predator stomach at each location was calculated. Data for each station included station ID, year, season, date, latitude, longitude, vessel, mean bluefish prey weight (g), mean piscivore length (cm), number of piscivore species, and sea surface temperature (degrees C). Because approximately 10% of survey stations were missing in-situ sea water temperature measurements, National Oceanic and Atmospheric Administration Optimum Interpolation Sea Surface Temperature (NOAA OI SST) V2 High Resolution Dataset [@reynolds_daily_2007] data provided by the NOAA PSL, Boulder, Colorado, USA, from their website at https://psl.noaa.gov were used to fill gaps. For survey stations with in-situ temperature measurements, the in-situ measurement was retained. For survey stations with missing temperature data, OI SST was substituted for input into VAST models. Models were developed combining all data for the year ("Annual") and with separate data for "Spring" (collection months January - June) and "Fall" (collection months July-December) to align with assumptions used in the bluefish stock assessment. Modeled years included 1985-2021 to align with other data inputs in the bluefish stock assessment. SST is also likely to affect prey distribution, but differently for each prey species. Therefore, SST is not modeled as a density covariate for aggregate small pelagics. --- ## Spatial partitioning: examining small pelagics trends at multiple scales <div class="figure"> <img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/maps-1.png" alt="Maps of key areas for Bluefish assessment indices. The full VAST model grid is shown in brown." width="33%" /><img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/maps-2.png" alt="Maps of key areas for Bluefish assessment indices. The full VAST model grid is shown in brown." width="33%" /><img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/maps-3.png" alt="Maps of key areas for Bluefish assessment indices. The full VAST model grid is shown in brown." width="33%" /> <p class="caption">Maps of key areas for Bluefish assessment indices. The full VAST model grid is shown in brown.</p> </div> Indices for aggregate small pelagics from piscivore stomachs can be calculated for any subset of the full model domain. Bias correction of the resulting indices is then applied <a name=cite-thorson_implementing_2016></a>([Thorson, et al., 2016](https://www.sciencedirect.com/science/article/pii/S0165783615301399)). ??? NEFSC survey strata definitions are built into the VAST `northwest-atlantic` extrapolation grid already. We defined additional new strata to address the recreational inshore-offshore 3 mile boundary. The area within and outside 3 miles of shore was defined using the `sf` R package as a 3 nautical mile (approximated as 5.556 km) buffer from a high resolution coastline from the`rnaturalearth` R package. This buffer was then intersected with the current `FishStatsUtils::northwest_atlantic_grid` built into VAST and saved using code [here](https://github.com/sgaichas/bluefishdiet/blob/main/VASTcovariates_updatedPreds_sst_3mi.Rmd#L49-L94). Then, the new State and Federal waters strata were used to split NEFSC survey strata where applicable, and the new full set of strata were used along with a modified function from `FishStatsUtils::Prepare_NWA_Extrapolation_Data_Fn` to build a custom extrapolation grid for VAST as described in detail [here](https://sgaichas.github.io/bluefishdiet/VASTcovariates_finalmodbiascorrect_3misurvstrat.html). --- background-image: url("https://github.com/sgaichas/bluefishdiet/raw/main/pyindex/allagg_fall_500_lennosst_ALLsplit_biascorrect/ln_density-predicted.png") background-size: 560px background-position: right top ## Results: Fall Forage Index .pull-left[ <div class="figure"> <img src="20221207_BluefishRT_ForageIndex_Gaichas_files/figure-html/WHAMfall-1.png" alt="Time series of VAST estimated fall forage indices for input into the bluefish assessment, 1985-2021" width="504" /> <p class="caption">Time series of VAST estimated fall forage indices for input into the bluefish assessment, 1985-2021</p> </div> VAST estimated Fall forage biomass density → ] .pull-right[ ] ??? --- ## Thank you! ## References .contrib[ <a name=bib-collette_bigelow_2002></a>[Collette, B. B. et al.](#cite-collette_bigelow_2002) (2002). _Bigelow and Schroeder's Fishes of the Gulf of Maine, Third Edition_. 3rd ed. edition. Washington, DC: Smithsonian Books. ISBN: 978-1-56098-951-6. <a name=bib-ng_predator_2021></a>[Ng, E. L. et al.](#cite-ng_predator_2021) (2021). "Predator stomach contents can provide accurate indices of prey biomass". In: _ICES Journal of Marine Science_ 78.3, pp. 1146-1159. ISSN: 1054-3139. DOI: [10.1093/icesjms/fsab026](https://doi.org/10.1093%2Ficesjms%2Ffsab026). URL: [https://doi.org/10.1093/icesjms/fsab026](https://doi.org/10.1093/icesjms/fsab026) (visited on Sep. 01, 2021). <a name=bib-reynolds_daily_2007></a>[Reynolds, R. W. et al.](#cite-reynolds_daily_2007) (2007). "Daily High-Resolution-Blended Analyses for Sea Surface Temperature". EN. In: _Journal of Climate_ 20.22. Publisher: American Meteorological Society Section: Journal of Climate, pp. 5473-5496. ISSN: 0894-8755, 1520-0442. DOI: [10.1175/2007JCLI1824.1](https://doi.org/10.1175%2F2007JCLI1824.1). URL: [https://journals.ametsoc.org/view/journals/clim/20/22/2007jcli1824.1.xml](https://journals.ametsoc.org/view/journals/clim/20/22/2007jcli1824.1.xml) (visited on Aug. 01, 2022). <a name=bib-thorson_guidance_2019></a>[Thorson, J. T.](#cite-thorson_guidance_2019) (2019). "Guidance for decisions using the Vector Autoregressive Spatio-Temporal (VAST) package in stock, ecosystem, habitat and climate assessments". En. In: _Fisheries Research_ 210, pp. 143-161. ISSN: 0165-7836. DOI: [10.1016/j.fishres.2018.10.013](https://doi.org/10.1016%2Fj.fishres.2018.10.013). URL: [http://www.sciencedirect.com/science/article/pii/S0165783618302820](http://www.sciencedirect.com/science/article/pii/S0165783618302820) (visited on Feb. 24, 2020). <a name=bib-thorson_comparing_2017></a>[Thorson, J. T. et al.](#cite-thorson_comparing_2017) (2017). "Comparing estimates of abundance trends and distribution shifts using single- and multispecies models of fishes and biogenic habitat". In: _ICES Journal of Marine Science_ 74.5, pp. 1311-1321. ISSN: 1054-3139. DOI: [10.1093/icesjms/fsw193](https://doi.org/10.1093%2Ficesjms%2Ffsw193). URL: [https://doi.org/10.1093/icesjms/fsw193](https://doi.org/10.1093/icesjms/fsw193) (visited on Nov. 04, 2021). <a name=bib-thorson_implementing_2016></a>[Thorson, J. T. et al.](#cite-thorson_implementing_2016) (2016). "Implementing a generic method for bias correction in statistical models using random effects, with spatial and population dynamics examples". En. In: _Fisheries Research_ 175, pp. 66-74. ISSN: 0165-7836. DOI: [10.1016/j.fishres.2015.11.016](https://doi.org/10.1016%2Fj.fishres.2015.11.016). URL: [https://www.sciencedirect.com/science/article/pii/S0165783615301399](https://www.sciencedirect.com/science/article/pii/S0165783615301399) (visited on Jul. 29, 2022). ] ## Additional resources .contrib[ [Northeast US State of the Ecosystem Reports](https://www.fisheries.noaa.gov/new-england-mid-atlantic/ecosystems/state-ecosystem-reports-northeast-us-shelf) ] .footnote[ Slides available at https://noaa-edab.github.io/presentations Contact: <Sarah.Gaichas@noaa.gov> ]