diff --git a/404.html b/404.html index 2a115817..63edca0c 100644 --- a/404.html +++ b/404.html @@ -23,7 +23,7 @@ - + @@ -54,9 +54,10 @@ + - + @@ -278,13 +279,13 @@
Description: Seasonal bottom temperatures on the Northeast Continental Shelf between 1959 and 2022 in a 1/12° grid.
Indicator category: Published Methods, Synthesis of published information
-Found in: State of the Ecosystem - Gulf of Maine & Georges Bank (2023); State of the Ecosystem - Mid-Atlantic Bight (2023)
-Contributor(s): Hubert du Pontavice, Vincent Saba, Zhuomin Chen
-Data steward: Hubert du Pontavice, hubert.dupontavice@noaa.gov
-Point of contact: Hubert du Pontavice, hubert.dupontavice@noaa.gov
+Found in: State of the Ecosystem - Gulf of Maine & Georges Bank (2023+); State of the Ecosystem - Mid-Atlantic Bight (2023+)
+Contributor(s): Joe Caracappa, Hubert du Pontavice, Vincent Saba, Zhuomin Chen
+Data steward: Joe Caracappa, joseph.caracappa@noaa.gov
+Point of contact: Joe Caracappa, joseph.caracappa@noaa.gov
Public availability statement: Source data are NOT publicly available. Please email hubert.dupontavice@noaa.gov for further information and queries of bottom temperature source data.
Four ocean models were used to get high-resolution daily bottom temperature on the NEUS between 1959 and 2022.
-For the period between 1959 and 1992, we used daily ocean bottom temperature from the long-term (1958–2007) high-resolution numerical simulation of the Northwest Atlantic Ocean in the Regional Ocean Modelling System (ROMS), a split-explicit, free-surface, terrain-following, hydrostatic, primitive equation model (Shchepetkin and McWilliams (2005)). The model domain covers the Northwest Atlantic Ocean with ~7km horizontal resolution and 40 vertical terrain- following layers. A detailed description of ROMS-NWA can be found in Chen et al. (2018).
+For the period between 1959 and 1992, we used daily ocean bottom temperature from the long-term (1958–2007) high-resolution numerical simulation of the Northwest Atlantic Ocean in the Regional Ocean Modelling System (ROMS), a split-explicit, free-surface, terrain-following, hydrostatic, primitive equation model (Shchepetkin and McWilliams (2005)). The model domain covers the Northwest Atlantic Ocean with ~7km horizontal resolution and 40 vertical terrain- following layers. A detailed description of ROMS-NWA can be found in Chen et al. (2018a).
For the period between 1992 and 2020, the daily bottom temperature outputs from the GLORYS12v1 ocean reanalysis product were used. GLORYS12v1 is a global ocean, eddy-resolving, and data assimilated hindcast from Mercator Ocean (European Union-Copernicus Marine Service, 2018; Fernandez and Lellouche2018; Jean-Michel et al. (2021a)) with 1/12 degree horizontal resolution and 50 vertical levels. The base ocean model is the Nucleus for European Modelling of the Ocean 3.1 (NEMO 3.1; Madec, 2016) driven at the surface by the European Centre for the Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis (Dee et al. (2011)). Remotely sensed and in situ observations are jointly assimilated by means of a reduced-order Kalman filter.
-For the year 2021, we used daily bottom temperature from the Operational Mercator global ocean analysis and forecast system (GLO12v3 called PSY4V3R1 in “A High-Resolution Ocean Bottom Temperature Product for the Northeast u.s. Continental Shelf Marine Ecosystem” (2023) and Lellouche et al. (2018)). GLO12v3 is a global ocean, eddy-resolving, monitoring forecasting system (Lellouche et al. (2018)) with the same ocean model grid (1/12 degree horizontal resolution and 50 vertical levels) and has many similarities with GLORYS12v1. Remotely sensed and in situ observations are also jointly assimilated by means of a reduced-order Kalman filter.
-For the year 2022, we used GLO12v4 which is a revised and updated version of GLO12v3 (European Union-Copernicus Marine Service, 2016). The general model structure is similar to GLO12v3 with some changes in model configuration, parameterizations, relaxations to avoid spurious drifts, river inputs, atmospheric fluxes and data assimilation (more detail in https://data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_PHY_001_024/description)
+For the year 2021 and 2022, we used GLO12v4 which is a revised and updated version of GLO12v3 (European Union-Copernicus Marine Service, 2016). The general model structure is similar to GLO12v3 with some changes in model configuration, parameterizations, relaxations to avoid spurious drifts, river inputs, atmospheric fluxes and data assimilation (more detail in https://data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_PHY_001_024/description)
Description: Chesapeake Bay Seasonal SST Anomalies
Found in: State of the Ecosystem - Mid-Atlantic (2021+)
Indicator category: Database pull with analysis
@@ -929,32 +930,25 @@Public availability statement: Source data are publicly available.
Public availability statement: Source data are publicly available here.
Data for Chesapeake Bay seasonal sea surface temperature (SST) anomalies were derived from the NOAA Multi-satellite AVHRR SST data set, available from NOAA CoastWatch East Coast Regional Node. The data set is a composite of overpasses from all operational satellites currently flying the Advanced Very High Resolution Radiometer (AVHRR) instrument. SST is derived using the Operational Non-linear Multichannel SST Algorithm (X. Li, Pichel, Maturi, et al. (2001), X. Li, Pichel, Clemente-Colón, et al. (2001)). Both daytime and nighttime overpasses are composited into daily and then seasonal SST products. The data extend from 2008 to present, and provide a 1.25 km x 1.25 km grid of SST measurements.
+Data for the Chesapeake Bay seasonal sea surface temperature (SST) anomalies are derived from the AVHRR and VIIRS Multi-Sensor Composite Sea Surface Temperature data set, available from the NOAA CoastWatch East Coast Regional Node (https://eastcoast.coastwatch.noaa.gov/). The data set is a composite of overpasses from two instruments: the Advanced Very High Resolution Radiometer (AVHRR) instrument on the European MetOp satellites, and the Visible Infrared Imaging Radiometer Suite (VIIRS) on the S-NPP and NOAA satellites (starting with NOAA-20 and follow-on NOAA satellites). Data from all current operational satellites are used in order to increase geographic coverage on a per-day basis. SST is derived using the Advanced Clear-Sky Processor for Oceans (ACSPO) processing system for consistency across instruments (Petrenko et al., 2014; Petrenko et al., 2016). Only nighttime overpasses are incorporated into the composite, i.e. the data do not represent daytime solar heating of the water surface. The data extend from December 2006 to the present. The AVHRR 1 km spatial resolution data, and the VIIRS 750 m spatial resolution data, are co-gridded to an 830 m spatial grid.
+More information about the AVHRR and VIIRS Multi-Sensor Composite SST data set is available at: +https://eastcoast.coastwatch.noaa.gov/cw_avhrr-viirs_sst.php
Anomaly maps of SST are generated by creating long-term ‘climatological’ seasonal average SST for the years from 2008 to the year immediately prior to the current year (max(Year) - 1)
. The reference period serves as a benchmark for comparing current observations. The current-year seasonal SST is then subtracted from the long-term seasonal average. Seasons for Chesapeake Bay are Dec-Feb (winter), Mar-May (spring), Jun-Aug (summer), and Sep-Nov (fall).
Individual nighttime overpasses from all instruments on the current operational satellites are composited into daily gridded scenes. The daily gridded scenes are then averaged seasonally. A long-term ‘climatological’ seasonal average SST is generated for each season for the years from 2007 to the year immediately prior to the current year (max(Year) - 1). This reference period serves as a benchmark for comparing current observations. The current-year seasonal SST is then subtracted from the long-term seasonal average for all grid cells to create the anomaly map. Seasons for Chesapeake Bay are Dec-Feb (winter), Mar-May (spring), Jun-Aug (summer), and Sep-Nov (fall).
Code for processing Chesapeake Bay temperature data can be found here.
catalog link -<https://noaa-edab.github.io/catalog/ch_bay_temp.html>
+<https://noaa-edab.github.io/catalog/ches_bay_sst.html>The first step was to define the Cold Pool domain, which is typically located within the MAB and the southern flank of Georges Bank (Chen et al. (2018); Robert W. Houghton et al. (1982); Lentz (2017a)). Here, we delineated a spatial domain covering the management area of the SNEMA yellowtail flounder (since this method was initially developed to study the Cold Pool impact on yellowtail flounder recruitment) comprising the MAB and in the SNE shelf between the 20 and 200 m isobaths (Chen et al. (2018); Chen and Curchitser (2020)). We restricted the time period from June (to match the start of the settlement period; SULLIVAN, COWEN, and STEVES (2005)) to September (which is the average end date of the Cold Pool (calendar day 269) estimated by Chen and Curchitser (2020). The Cold Pool domain was defined as the area, wherein average bottom temperature was cooler than 10°C between June and September from 1959 to 2022. We then developed the three Cold Pool indices using bottom temperature from ocean models.
+The first step was to define the Cold Pool domain, which is typically located within the MAB and the southern flank of Georges Bank (Chen et al. (2018b); Robert W. Houghton et al. (1982); Lentz (2017a)). Here, we delineated a spatial domain covering the management area of the SNEMA yellowtail flounder (since this method was initially developed to study the Cold Pool impact on yellowtail flounder recruitment) comprising the MAB and in the SNE shelf between the 20 and 200 m isobaths (Chen et al. (2018b); Chen and Curchitser (2020)). We restricted the time period from June (to match the start of the settlement period; SULLIVAN, COWEN, and STEVES (2005)) to September (which is the average end date of the Cold Pool (calendar day 269) estimated by Chen and Curchitser (2020). The Cold Pool domain was defined as the area, wherein average bottom temperature was cooler than 10°C between June and September from 1959 to 2022. We then developed the three Cold Pool indices using bottom temperature from ocean models.
The spatial component of the Cold Pool and the habitat provided by the cold pool was calculated using the Spatial Extent Index (Model_SEI). The Model_SEI is estimated by the number of cells where bottom temperature remains below 10C for at least 2 months between June and September.
The Bottom temperature data is the average ROMS-NWA bottom temperature over the decade \[d\] in the grid cell \[i\]. All above methods duPontavice et al. (2022).
Bottom temperature from Glorys reanalysis and Global Ocean Physics Analysis were not being processed.
-Bottom temperature from ROMS-NWA (used for the period 1959-1992) were bias-corrected. Previous studies that focused on the ROMS-NWA-based Cold Pool highlighted strong and consistent warm bias in bottom temperature of about 1.5C during the stratified seasons over the period of 1958-2007 (Chen et al. (2018); Chen and Curchitser (2020)). In order to bias-correct bottom temperature from ROMS-NWA, we used the monthly climatologies of observed bottom temperature from the Northwest Atlantic Ocean regional climatology (NWARC) over decadal periods from 1955 to 1994. The NWARC provides high resolution (1/10° grids) of quality-controlled in situ ocean temperature based on a large volume of observed temperature data (Seidov, Baranova, Johnson, et al. (2016), Seidov, Baranova, Boyer, et al. (2016)) (https://www.ncei.noaa.gov/products/northwest-atlantic-regional-climatology). The first step was to re-grid the ROMS-NWA to obtain bottom temperature over the same 1/10° grid as the NWARC. Then, a monthly bias was calculated in each grid cell and for each decade (1955–1964, 1965–1974, 1975–1984, 1985–1994) in the MAB and in the SNE shelf:
+Bottom temperature from ROMS-NWA (used for the period 1959-1992) were bias-corrected. Previous studies that focused on the ROMS-NWA-based Cold Pool highlighted strong and consistent warm bias in bottom temperature of about 1.5C during the stratified seasons over the period of 1958-2007 (Chen et al. (2018b); Chen and Curchitser (2020)). In order to bias-correct bottom temperature from ROMS-NWA, we used the monthly climatologies of observed bottom temperature from the Northwest Atlantic Ocean regional climatology (NWARC) over decadal periods from 1955 to 1994. The NWARC provides high resolution (1/10° grids) of quality-controlled in situ ocean temperature based on a large volume of observed temperature data (Seidov, Baranova, Johnson, et al. (2016), Seidov, Baranova, Boyer, et al. (2016)) (https://www.ncei.noaa.gov/products/northwest-atlantic-regional-climatology). The first step was to re-grid the ROMS-NWA to obtain bottom temperature over the same 1/10° grid as the NWARC. Then, a monthly bias was calculated in each grid cell and for each decade (1955–1964, 1965–1974, 1975–1984, 1985–1994) in the MAB and in the SNE shelf:
\[{BIAS}_{i,\ d}=\ T_{i,d}^{Climatology}\ -\ {\bar{T}}_{i,\ d}^{ROMS-NWA}\]
where \[T_{i,d}^{Climatology}\] is the NWARC bottom temperature in the grid cell i for the decade d and \[{\bar{T}}_{i,\ d}^{ROMS-NWA}\] is the average ROMS-NWA bottom temperature over the decade d in the grid cell i. All above methods duPontavice et al. (2022).
Depth-averaged spatial temperature is calculated based on the daily Cold Pool dataset, which is quantified following Chen et al. (2018).
+Depth-averaged spatial temperature is calculated based on the daily Cold Pool dataset, which is quantified following Chen et al. (2018b).
The spatial distribution and abundance of cyst cells from the annual survey are used to drive an ecosystem forecast model for the Gulf of Maine (Anderson et al. (2005), Y. Li et al. (2009), Y. Li et al. (2020), McGillicuddy et al. (2011)). The model also includes many other inputs of dynamic oceanographic data such as currents, temperature, and nutrients. Operational Harmful Algal Bloom forecast is served online at https://coastalscience.noaa.gov/research/stressor-impacts-mitigation/hab-forecasts/gulf-of-maine-alexandrium-catenella-predictive-models/.
+The spatial distribution and abundance of cyst cells from the annual survey are used to drive an ecosystem forecast model for the Gulf of Maine (Anderson et al. (2005), Li et al. (2009), Li et al. (2020), McGillicuddy et al. (2011)). The model also includes many other inputs of dynamic oceanographic data such as currents, temperature, and nutrients. Operational Harmful Algal Bloom forecast is served online at https://coastalscience.noaa.gov/research/stressor-impacts-mitigation/hab-forecasts/gulf-of-maine-alexandrium-catenella-predictive-models/.
Description: Harbor Porpoise and Gray Seal Indicator
+Description: Harbor Porpoise Bycatch Indicator
Found in: State of the Ecosystem - Gulf of Maine & Georges Bank (2018+), State of the Ecosystem - Mid-Atlantic (2018+)
Indicator category: Synthesis of published information; Published methods
-Contributor(s): Christopher D. Orphandies, Debra Palka
+Contributor(s): Kristin Precoda, Christopher D. Orphanides, Debra Palka
Data steward: Debra Palka debra.palka@noaa.gov
Point of contact: Debra Palka debra.palka@noaa.gov
Public availability statement: Source data are available in public stock assessment reports.
@@ -931,7 +932,7 @@Reported harbor porpoise bycatch estimates and potential biological removal levels can be found in publicly available documents; detailed in Marine Mammal Protection Stock Assessments. More detailed documentation as to the methods employed can be found in NOAA Fisheries Northeast Fisheries Science Center (NEFSC) Center Reference Documents (CRDs) found on the NEFSC publications page.
+Reported harbor porpoise bycatch estimates and potential biological removal levels can be found in publicly available documents; detailed in Marine Mammal Protection Stock Assessments. More detailed documentation as to the methods employed can be found in NOAA Fisheries Northeast Fisheries Science Center (NEFSC) Center Reference Documents (CRDs) found on the NEFSC publications page.
Bycatch estimates as found in stock assessment reports were plotted along with confidence intervals. The confidence intervals were calculated from published CVs assuming a normal distribution (\(\sigma = \mu CV\); \(CI = \bar{x} \pm \sigma * 1.96\)).
-Data were analyzed and formatted for inclusion in the ecodata
R package using the R code found here, Harbor Porpoise and Gray Seal.
Data were analyzed and formatted for inclusion in the ecodata
R package using the R code found here, Harbor Porpoise.
catalog link -https://noaa-edab.github.io/catalog/harborporpoise.html -https://noaa-edab.github.io/catalog/grayseal.html
+https://noaa-edab.github.io/catalog/harborporpoise.htmlFound in: State of the Ecosystem - Gulf of Maine & Georges Bank (2020(Different Methods), 2021+), State of the Ecosystem - Mid-Atlantic (2020(Different Methods), 2021+)
Indicator category: Database pull with analysis
Contributor(s): Heather Baertlein, Jackie Wilson, George Silva, Jennifer Cudney
-Data steward: Kimberly Bastille
+Data steward: Jennifer Cudney jennifer.cudney@noaa.gov
Point of contact: Jennifer Cudney jennifer.cudney@noaa.gov
Public availability statement: Source data are NOT publicly available.
Data from eDealer database (https://www.fisheries.noaa.gov/atlantic-highly-migratory-species/atlantic-highly-migratory-species-dealer-reporting) and Bluefin Tuna Dealer reports on SAFIS (https://www.accsp.org/what-we-do/safis/). The eDealer data were supplemented with ACCSP records, GulfFIN records, and vessel logbook catches for which no dealer reports were submitted.
+Data from eDealer database (https://www.fisheries.noaa.gov/atlantic-highly-migratory-species/atlantic-highly-migratory-species-dealer-reporting) and Bluefin Tuna Dealer reports on SAFIS (https://www.accsp.org/what-we-do/safis/). The eDealer data were supplemented with ACCSP records, GulfFIN records, and pelagic and coastal fisheries vessel logbook catches reported to SEFSC for which no dealer reports were submitted.
Data were processed for Fisheries of the United States and then aggregated by regions to avoid confidentiality issues. Data of Atlantic shark, swordfish, bigeye tuna, albacore tuna, yellowfin tuna and skipjack tuna were initially extracted from our eDealer database. Additional landings of these HMS not in eDealer were found in ACCSP, GulfFIN, and the SEFSC Atlantic HMS vessel logbook records. Bluefin tuna landings data from the Bluefin Tuna Dealer reports in SAFIS were also extracted and combined with the eDealer data for other HMS .
-Procedures of quality assurance were conducted. Duplicate records were removed from the data. This may occur from multiple submissions of reports by the same dealer. It may also occur when two or more dealers report the same landings in “Packing” situations. While most vessels immediately sell their catch to the dealer at their port of landing, some vessels sell their catch to a dealer(s) in another location. Transport to alternate locations requires processing of the fish to preserve quality. This processing activity is done by the dealer at the port of landing and is referred to as “Packing”. Differences in federal and state definitions of who is considered the “dealer” of the product, and thus ultimately responsible for submitting the landings report, often results in multiple reports being created for the same landings. These duplicate reports need to be accounted for when summarizing the data to reflect accurate landings. Therefore, searches for duplicate reports of the same landing were conducted and eliminated prior to summarizing the data for the Fisheries of the United States.
-All reported landings were converted to live weights using conversion ratios appropriate for the species/species group and reported grade of the product. Shark fins were not reported to live weight as these weights are included in the converted whole weight of the reported shark landing.
-States, where the landings occurred, were grouped into ‘ecological production units’ (EPUs), as defined by GARFO staff. “New England”, or NE, includes Maine, New Hampshire and Massachusetts, as well as landings from Canada. The “Mid-Atlantic Bight”, or MAB, includes states from Rhode Island to North -Carolina.
+Data, from 2015-2022, were processed for Fisheries of the United States and then aggregated by regions to avoid confidentiality issues. Data of Atlantic shark, swordfish, bigeye tuna, albacore tuna, yellowfin tuna and skipjack tuna were initially extracted from our eDealer database and the Caribbean Commercial Vessel logbook database. Additional landings of these HMS not in this dataset were collected from ACCSP, GulfFIN, and the SEFSC Atlantic HMS vessel logbook databases. Bluefin tuna landings data from the Bluefin Tuna Dealer reports in SAFIS were also extracted and combined with the eDealer data for other HMS.
+Procedures of quality assurance were conducted. Duplicate records were removed from the data. This may occur from multiple submissions of reports by the same dealer. It may also occur when two or more dealers report the same landings in “Packing” situations. While most vessels immediately sell their catch to the dealer at their port of landing, some vessels sell their catch to a dealer(s) in another location. Transport to alternate locations requires processing of the fish to preserve quality. This processing activity is done by the dealer at the port of landing and is referred to as “Packing”. Differences in federal and state reporting requirements and definitions of who is considered the “dealer” of the product, and thus ultimately responsible for submitting the landings report, may result in multiple reports being created for the same landings. These duplicate reports need to be accounted for when summarizing the data to reflect accurate landings prior to summarizing the data for analyses, including the Fisheries of the United States.
+All reported landings were converted to live weights using conversion ratios appropriate for the species/species group and reported grade of the product. Shark fins, shark heads, shark tails, and shark bellies were not reported to live weight as these weights are included in the converted whole weight of the reported shark landing.
+States, where the landings occurred, were grouped into ‘ecological production units’ (EPUs), as defined by GARFO staff. “New England”, or NE, includes Maine, New Hampshire and Massachusetts, as well as landings from Canada. The “Mid-Atlantic Bight”, or MAB, includes states from Rhode Island to North Carolina. Landings in states outside of these EPUs were excluded from further summaries.
Seven HMS Management Groups represent 26 highly migratory species in the dataset. HMS Management Groups may include a single species or a group of species. HMS groups include “Bluefin Tuna”, “BAYS”, “Swordfish”, “Large Coastal Sharks”, “Small Coastal Sharks”, “Pelagic Sharks”, “Smoothhound Sharks”. “BAYS” includes bigeye, albacore, yellowfin and skipjack tunas. “Large Coastal Sharks” includes blacktip, bull, great hammerhead, scalloped hammerhead, smooth hammerhead, lemon, nurse, sandbar, silky, spinner, and tiger sharks. “Small Coastal Sharks” includes Atlantic sharpnose, blacknose, bonnethead, finetooth sharks. “Pelagic Sharks” includes blue, porbeagle, shortfin mako, and thresher sharks. “Smoothhound Sharks” includes smooth dogfish shark.
-Price per pound was used to determine the ex-vessel value. For landings with prices per pound reported as “N/A”, 0, $0.01 or left blank, average prices were calculated for each species and state. Those averages replaced the missing values to determine landings revenue. Revenue from sales to the aquarium trade were also excluded to avoid extreme values associated with shipping live specimens.
-High migratory landings include 26 species of tunas, sharks and swordfish.
-Data were processed and analyzed using SAS and Microsoft Excel pivot tables. -The count of dealers and vessels in each regional species/management group sum was used to determine if a sufficient number of records were available to make the data public or if it needed to be marked as confidential. Additionally, ratios of landings reported by dealers/fishermen were compared in each regional species/managment group sum to determine if any one entity contributed more than 2/3 of the total which would require it being marked as confidential.
+Price per pound was used to determine the ex-vessel value. For landings with a reported disposition of “Food” and prices per pound reported as “N/A”, 0, $0.01 or left blank, average prices were calculated for each species and state. Those average prices replaced the missing values to determine landings revenue. Revenue from sales to the aquarium trade were also excluded to avoid extreme values associated with shipping live specimens.
Highly migratory species landings include 26 species of tunas, sharks and swordfish. Data were processed and analyzed using SAS and Microsoft Excel pivot tables. The count of entities represented was used to determine if a sufficient number of records were available to make the data public or if it needed to be marked as confidential; generally, data was marked confidential if it did not meet the “rule of three” (i.e., at least three unique entities represented).
HMS landings data were formatted for inclusion in the ecodata
R package using the R code found here.
22 March 2024
+5 April 2024
Point of contact: Richard Pace, richard.pace@noaa.gov
+Contributor(s): Debra Palka
+Data steward: Debra Palka, debra.palka@noaa.gov
+Point of contact: Daniel Linden, daniel.linden@noaa.gov
Public availability statement: Source data are available from the New England Aquarium upon request. Derived data are available here.
The North Atlantic right whale abundance estimates were taken from a published document (see Pace, Corkeron, and Kraus 2017), except for the most recent 2016 and 2017 estimates. Abundance estimates from 2016 and 2017 were taken from the 2016 National Oceanographic and Atmospheric Administration marine mammal stock assessment (Hayes et al. 2017) and an unpublished 2017 stock assessment.
-Calves birth estimates are taken from a published report (Pettis, Pace, and Hamilton 2019) put out yearly by the North American Right Whale Consortium.
+The North Atlantic right whale abundance estimates were taken from a published document (see Linden_pop_2023?).
+Calves birth estimates are available in Pace, Corkeron, and Kraus (2017), with more recent years shared here.
Analysis for right whale abundance estimates is provided by Pace, Corkeron, and Kraus (2017), and code can be found in the supplemental materials.
+Analysis for right whale abundance estimates is based on methods by Pace, Corkeron, and Kraus (2017), as documented most recently by (Linden_pop_2023?). Data and code can be found in the following Github repository: NEFSC/PSD-NARW_popsize.
Time series of right whale and calf abundance estimates were formatted for inclusion in the ecodata
R package using this R code.
catalog link -https://noaa-edab.github.io/catalog/narw.html
Point of contact: Grace Saba saba@marine.rutgers.edu
+Contributor(s): Grace Saba, Lori Garzio
+Data steward: Grace Saba saba@marine.rutgers.edu, Lori Garzio lgarzio@marine.rutgers.edu
+Point of contact: Grace Saba saba@marine.rutgers.edu, Lori Garzio lgarzio@marine.rutgers.edu
Public availability statement: Source data is available to the public (see Data Sources).
The New England Fishery Management Council (NEFMC) and Mid-Atlantic Fishery Management Council (MAFMC) have recently requested regional Ocean Acidification (OA) information in the State of the Ecosystem reports. The work included in the State of the Ecosystem 2021 report, seasonal dynamics of pH in shelf waters in the Mid-Atlantic, was synthesized from Wright-Fairbanks et al. (2020). The maps included in the State of the Ecosystem 2022 reports include a plot of bottom pH in summer over the entire U.S. Northeast Shelf (2007-2021), and glider-based pH profiles during summer 2021 in both the Mid-Atlantic Bight (MAB) and the Gulf of Maine. The plots in the State of the Ecosystem 2023 reports included maps of bottom summer aragonite saturation (ΩArag, or omega) over the entire U.S. Northeast Shelf (2007-2022) and locations where summer bottom ΩArag reached lab-derived sensitivity levels of designated target species.
The products developed for the State of the Ecosystem 2024 reports include: plots of the seasonal progression (Spring-Fall 2023) of oceanographic properties (including temperature, chlorophyll, dissolved oxygen, pH, and aragonite saturation state) on the New Jersey coastal shelf; plots summarizing a multi-stressor event in the Mid-Atlantic during summer 2023; static and animated maps of summer-time bottom pH and aragonite saturation state on the U.S. Northeast Shelf (2007-2023); and maps of locations where species sensitivity levels for aragonite saturation state were reached in bottom water during the summer (2007-2023).
-Products from all State of the Ecosystem reports to date were developed using openly accessible, quality-controlled data from vessel-based discrete samples and glider-based measurements (see Data Sources), and data from published laboratory-based experimental studies (see Plotting).
+Products from all State of the Ecosystem reports to date were developed using openly accessible, quality-controlled data from vessel-based discrete samples and glider-based measurements (see Data Sources), and data from published laboratory-based experimental studies (see Plotting)
Glider-based observations of pH (and other variables including temperature, salinity, chlorophyll-a, and dissolved oxygen) began in the southern MAB region in May 2018 (Saba et al. 2019), and seasonal glider pH missions thereafter began in February 2019 (Wright-Fairbanks et al. 2020; although no deployments occurred in 2020 as a result of the COVID pandemic). Simultaneous measurements from the glider’s pH, temperature, and salinity sensors enable the derivation of total alkalinity and calculation of other carbonate system parameters including aragonite saturation state (ΩArag). The glider pH observation program expanded spatially and temporally, with additional deployments in the northern MAB (New York Bight) and the Gulf of Maine, starting in February 2021. A typical glider mission runs for about 4 weeks, covers 500 km, and collects data though the full water column.
-For the 2023 glider data cross-section plots (2023 seasonal progression of oceanographic properties on the New Jersey coastal shelf and plots summarizing a multi-stressor event in the Mid-Atlantic during summer 2023): -Full-resolution delayed-mode glider datasets from six deployments (Spring-Fall 2023) were downloaded from the RUCOOL Glider ERDDAP server. These datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata can be found here.
-For the summer 2007-2023 bottom-water pH and aragonite saturation state data synthesis (static and animated maps of summer-time bottom pH and aragonite saturation state on the U.S. Northeast Shelf [2007-2023]; maps of summer-time bottom locations where species sensitivity levels for aragonite saturation state were reached [2007-2023]): -Full-resolution delayed-mode glider datasets from seven deployments were downloaded from the RUCOOL Glider ERDDAP server and the IOOS Glider DAC ERDDAP server (SBU01 2022-2023 deployments). These datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata for each glider deployment used in the synthesis can be found here. Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from glider-based (and vessel-based, see below) measurements can be found here. -Vessel-based data were mined from the Coastal Ocean Data Analysis Product in North America, version v2021 (Jiang et al. 2021). This data product synthesizes two decades of quality-controlled inorganic carbon system parameters (including pH, total alkalinity, dissolved inorganic carbon) along with other physical and chemical parameters (temperature, salinity, dissolved oxygen, nutrients) collected from the North American continental shelves. Additionally, four recent vessel-based datasets that were not included in CODAP-NA (Jiang et al. 2021) were included in this synthesis. These datasets were collected during more recent NOAA NEFSC Ecosystem Monitoring (EcoMon) surveys: June 2019 (Cruise ID HB1902), August 2019 (Cruise ID GU1902), August 2021 (Cruise ID PC2104), June 2022 (Cruise ID HB2204). These datasets were downloaded via the NCEI Ocean Carbon and Acidification Data Portal. Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from vessel-based (and glider-based, see above) measurements can be found here.
+For the 2023 glider data cross-section plots (2023 seasonal progression of oceanographic properties on the New Jersey coastal shelf and plots summarizing a multi-stressor event in the Mid-Atlantic during summer 2023): Full-resolution delayed-mode glider datasets from six deployments (Spring-Fall 2023) were downloaded from the RUCOOL Glider ERDDAP server. These datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata can be found here.
+For the summer 2007-2023 bottom-water pH and aragonite saturation state data synthesis (static and animated maps of summer-time bottom pH and aragonite saturation state on the U.S. Northeast Shelf [2007-2023]; maps of summer-time bottom locations where species sensitivity levels for aragonite saturation state were reached [2007-2023]): Full-resolution delayed-mode glider datasets from seven deployments were downloaded from the RUCOOL Glider ERDDAP server and the IOOS Glider DAC ERDDAP server (SBU01 2022-2023 deployments). These datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata for each glider deployment used in the synthesis can be found here. Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from glider-based (and vessel-based, see below) measurements can be found here. Vessel-based data were mined from the Coastal Ocean Data Analysis Product in North America, version v2021 (Jiang et al. 2021). This data product synthesizes two decades of quality-controlled inorganic carbon system parameters (including pH, total alkalinity, dissolved inorganic carbon) along with other physical and chemical parameters (temperature, salinity, dissolved oxygen, nutrients) collected from the North American continental shelves. Additionally, four recent vessel-based datasets that were not included in CODAP-NA (Jiang et al. 2021) were included in this synthesis. These datasets were collected during more recent NOAA NEFSC Ecosystem Monitoring (EcoMon) surveys: June 2019 (Cruise ID HB1902), August 2019 (Cruise ID GU1902), August 2021 (Cruise ID PC2104), June 2022 (Cruise ID HB2204). These datasets were downloaded via the NCEI Ocean Carbon and Acidification Data Portal. Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from vessel-based (and glider-based, see above) measurements can be found here.
Glider data were processed and quality-controlled by software technician Lori Garzio at Rutgers University.
-CODAP-NA data were accessed and downloaded on October 14, 2021.
+CODAP-NA, Version 2021 data were accessed and downloaded on October 14, 2021.
EcoMon datasets were accessed and downloaded on October 13, 2022 (Cruise IDs HB1902, GU1902, PC2104) and November 17, 2023 (Cruise ID HB2204).
For processing and quality-control procedures of glider-based data, see Wright-Fairbanks et al. (2020). For the 2023 glider data cross-section plots (2023 seasonal progression of oceanographic properties on the New Jersey coastal shelf and plots summarizing a multi-stressor event in the Mid-Atlantic during summer 2023): datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata can be found here.
+For processing and quality-control procedures of glider-based data, see Wright-Fairbanks et al. (2020). For the 2023 glider data cross-section plots (2023 seasonal progression of oceanographic properties on the New Jersey coastal shelf and plots summarizing a multi-stressor event in the Mid-Atlantic during summer 2023): datasets were QC’d and pH was calculated from raw variables (where applicable) using python code, and the final processed NetCDF datasets containing all relevant metadata can be found here.
Data from CODAP-NA were filtered temporally to include only those collected during summer months (June-August) and were spatially limited to the U.S. Northeast Shelf. The resulting datasets included those from major vessel-based campaigns (East Coast Ocean Acidification, ECOA I and II cruises 2015 and 2018; The Gulf of Mexico and East Coast Carbon cruises, GOMECC 2007 and 2012; EcoMon 2012-2013, 2015-2019).
For vessel-based datasets, when ΩArag was unavailable it was calculated using PyCO2SYS (Humphreys et al. 2020) with inputs of pressure, temperature, salinity, total alkalinity, and pH.
For MAB glider datasets, total alkalinity was calculated from salinity using a linear relationship determined from in situ water sampling data taken during glider deployment and recovery in addition to ship-based water samples (Wright-Fairbanks et al. 2020). For the Gulf of Maine glider dataset, total alkalinity was calculated from temperature and salinity using Table 3 Equation IV in McGarry et al. (2021). Calculations for ΩArag were then conducted using PyCO2SYS (Humphreys et al. 2020) with inputs of pressure, temperature, salinity, total alkalinity, and pH.
-Bottom water values were defined as the median of the measurements (or calculated ΩArag values) within the deepest 1m of a glider profile or, for vessel-based measurements, the deepest measurement of a vertical CTD/Rosette cast where water samples were collected, for profiles deeper than 10m. In order to validate whether the deepest depth was at or near the bottom, the sampling depth was compared to water column depth (when provided) or water depths extracted from a GEBCO bathymetry grid based on the sample collection coordinates. Any glider profiles/vessel-based casts with the deepest measurement shallower than the bottom 20% of total water column depth were removed. This allowed for a sliding scale instead of providing a strict cut off (e.g., 1 m above the bottom). Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from glider- and vessel-based measurements can be found here.
+Bottom water values were defined as the median of the measurements (or calculated ΩArag values) within the deepest 1m of a glider profile or, for vessel-based measurements, the deepest measurement of a vertical CTD/Rosette cast where water samples were collected, for profiles deeper than 10m. In order to validate whether the deepest depth was at or near the bottom, the sampling depth was compared to water column depth (when provided) or water depths extracted from a GEBCO bathymetry grid based on the sample collection coordinates. Any glider profiles/vessel-based casts with the deepest measurement shallower than the bottom 20% of total water column depth were removed. This allowed for a sliding scale instead of providing a strict cut off (e.g., 1 m above the bottom). Resulting data files containing combined summer-time bottom pH and aragonite saturation state data from glider- and vessel-based measurements can be found here.
A set of plots was constructed for the 2024 State of the Ecosystem reports:
Code for data manipulation and plotting can be found here: https://github.com/lgarzio/cinar-soe.
-catalog link -https://noaa-edab.github.io/catalog/ocean_acidification.html
+Data processing code can be found on Github here, and all data files use in these analyses and syntheses can be found here.
+catalog link https://noaa-edab.github.io/catalog/ocean_acidification.html
Two data frames are in the stocksmart
package, stockAssessmentData
and stockAssessmentSummary
.
In stockAssessmentData
we have time series. Columns are StockName, Stockid, Assessmentid, Year, Value, Metric, Description, Units, AssessmentYear, Jurisdiction, FMP, CommonName, ScientificName, ITIS, AssessmentType, StockArea, RegionalEcosystem and the reported metrics are Catch, Fmort, Recruitment, Abundance, Index.
In stockAssessmentSummary
we have assessment metadata. Columns are Stock ID, Stock Name, Jurisdiction, FMP, Science Center, Regional Ecosystem, FSSI Stock?, ITIS Taxon Serial Number, Scientific Name, Common Name, Stock Area, Assessment ID, Assessment Year, Assessment Month, Last Data Year, Review Result, Assessment Model, Model Version, Lead Lab, Citation, Final Assessment Report 1, Final Assessment Report 2, Point of Contact, Life History Data, Abundance Data, Catch Data, Assessment Level, Assessment Frequency, Model Category, Catch Input Data, Abundance Input Data, Biological Input Data, Ecosystem Linkage, Composition Input Data, F Year, Estimated F, F Unit, F Basis, Flimit, Flimit Basis, Fmsy, Fmsy Basis, F/Flimit, F/Fmsy, Ftarget, Ftarget Basis, F/Ftarget, B Year, Estimated B, B Unit, B Basis, Blimit, Blimit Basis, Bmsy, Bmsy Basis, B/Blimit, B/Bmsy, MSY, MSY Unit, Assessment Type.
For 2021-2023, stocksmart
was updated with all current assessments, so data extraction was simply:
Year-specific naming conventions for assess and decoder files were dropped in 2021 to facilitate future data updates.
In 2020, assessment summary data were extracted from stockAssessmentSummary
for 2019 and prior records, and the 2020 assessments results were added from the preliminary results provided by Jeff Vieser.