Discrete carbonate data from Narragansett Bay in 2022 and 2023

Website: https://www.bco-dmo.org/dataset/961940
Data Type: Other Field Results
Version: 1
Version Date: 2025-05-14

Project
» Drivers of ocean acidification in a temperate urbanized estuary undergoing nutrient loading reductions (OA Drivers in an Urbanized Estuary)
ContributorsAffiliationRole
Wang, HongjieUniversity of Rhode Island (URI)Principal Investigator
Davies, AndrewUniversity of Rhode Island (URI)Co-Principal Investigator
Gomes, KristoferUniversity of Rhode Island (URI)Scientist
Gu, ShuaiTexas A&M University (TAMU)Scientist
Stoffel, HeatherUniversity of Rhode Island (URI)Scientist
Ahumada, GeorgiaUniversity of MiamiStudent
Baskind, AbigailUniversity of Rhode Island (URI)Student
York, Amber D.Woods Hole Oceanographic Institution (WHOI BCO-DMO)BCO-DMO Data Manager

Abstract
This dataset contains discrete carbonate data (TA, DIC, pH, Temperature, Salinity), collected as part of the study described below. See the "Related Publications" section for autonomously collected pH data from this study. Study description: The increase in atmospheric carbon dioxide (CO2) over the last 200 years has largely been mitigated by the ocean’s function as a carbon sink. However, this continuous absorption of CO2 by seawater triggers ocean acidification (OA), a process in which water becomes more acidic and more depleted in carbonate ions that are essential for calcifiers. OA is well-studied in open ocean environments; however, understanding the unique manifestation of OA in coastal ecosystems presents myriad challenges due to considerable natural variability resulting from concurrent and sometimes opposing coastal processes--e.g. eutrophication, changing hydrological conditions, heterogeneous biological activity, and complex water mass mixing. This study analyzed high temporal resolution pH data collected during 2022 and 2023 from Narragansett Bay, RI--a mid-sized, urban estuary that since 2005 has undergone a 50% reduction in nitrogen loading with weekly, discrete bottle samples to verify sensor data. We used autonomous data for pH, temperature, salinity, and dissolved oxygen from 4 sensors in Narragansett Bay. The autonomous data spanned over a year from 2022 to mid-2023 and had temporal resolutions between 10 and 15 minutes. The data have been subjected to QA/QC protocols, such that all pH measurements are final and quality-controlled. As well, pH values normalized to 15°C (using PyCO2SYS) are included. All pH values are in total scale. Discrete samples were taken weekly at the Narragansett Bay Long Term Phytoplankton Time Series site and monthly from Greenwich Bay, collocated with 2 of the sensors. Discrete data were analyzed in lab for dissolved inorganic carbon and total alkalinity, and include in situ temperature and salinity.


Coverage

Location: Narragansett Bay, Rhode Island, United States
Spatial Extent: N:41.64 E:-71.3797 S:41.57 W:-71.39
Temporal Extent: 2022-04-28 - 2023-07-28

Methods & Sampling

This dataset contains discrete carbonate data (TA, DIC, pH, Temperature, Salinity). See the "Related Publications" section for autonomously collected pH data also described in this methodology.

The Mount View (MV) and Quonset Point (QP) sensors are maintained by the Narragansett Bay Fixed Station Monitoring Network (NBFSMN). Each sensor measures pH, along with temperature, salinity, and DO in the surface, every 15 and were deployed only from spring to fall.  Complete documentation provided by the Rhode Island Department of Environmental Management (RIDEM, 2020; https://dem.ri.gov/sites/g/files/xkgbur861/files/programs/benviron/water/quality/surfwq/pdfs/nbfsmn.pdf).

The Conanicut Hydrocat 620 and Potowomut Hydrocat 720 sensors are maintained by the Rhode Island Consortium for Coastal Ecology Assessment, Innovation, and Modeling (RI C-AIM; https://data-explorer.riddc.brown.edu/dataset/buoy-telemetry/summary) and Rhode Island Sea Grant. Each sensor measures pH, along with temperature, salinity, and DO, every 15 minutes approximately 1 meter below the surface (RI C-AIM, 2024). The hydrocat sensors were deployed year-round and were recovered every 2 to 4 weeks, weather permitting, for maintenance at which point the flow path and the conductivity cell are flushed with 1% Triton detergent followed by vigorous flushing with DI water.

Discrete bottle samples were collected to verify the sensor data. Samples for the southern region of the bay (i.e. near QP and Conanicut buoy) were collected weekly as part of the Narragansett Bay Long Term Phytoplankton Time Series (PLT), which is located just off the Conanicut buoy. The time series typically collects samples on Monday mornings at approximately 7:30am, barring weeks when inclement weather conditions delayed or canceled sampling. Samples for the northern region of the bay (i.e. Greenwich Bay near MV and Potowomut buoy) were collected approximately once a month. After collection, samples were poisoned with 100 uM  of saturated mercuric chloride solution and stored in the refrigerator until analysis for DIC and TA. DIC and TA were measured according to Dickson et al. (2007) using the Apollo SciTech Model AS-C6L Dissolved Inorganic Carbon Analyzer and the Apollo SciTech Model AS-ALK3 Total Alkalinity Titrator. Instruments were calibrated to Certified Reference Material (CRM) from Scripps Institute of Oceanography at room temperature (21 - 22° C). For DIC analysis, a sample of CRM was run prior to and at the end of sample analysis for quality control. TA analysis was calibrated to CRM, and either a sample of CRM or aged open ocean water was run at beginning of daily analysis and end of daily analysis for quality control. Lab-based measurements carry a ±0.2% uncertainty for alkalinity and ±0.1% uncertainty for DIC.

Instruments:
(additional context to accompany the list in the "Instruments" section)

The Mount View (MV; EXO V2) and Quonset Point (QP; 6600EDS) sensors are YSI brand multi-parameter sondes maintained by the Narragansett Bay Fixed Station Monitoring Network (NBFSMN). Both sondes were equipped to measure temperature, salinity, dissolved oxygen, pH [NBS], depth, and chlorophyll. The surface equipment was housed in a tube on a monitoring buoy station with the instrument depth about 1 meter from surface. All NBFSMN data were subject to quality assurance measures including verification of calibrations and consistency among multiple instruments, corrections for sensor drift and biases due to biofouling, removal of outliers, and interpolation across selected intervals of missing data, in accordance with the NBFSMN’s EPA approved Quality Assurance Project Plan (RIDEM, 2020). To keep consistency among instruments, servicing protocols consist of swapping the deployed instruments with newly calibrated instruments on a 2-week interval to minimize biofouling. All data are verified through a three-point comparison: data from the retrieved sonde are compared to the newly calibrated sonde, as well as an independent profiling sonde, all at the deployment depth. Outliers and data errors are removed based on criteria set in the NBFSMN’s Quality Assurance Project Plan. Data correction and gaps in coverage on average affect up to 6% of the record at an individual station for an annual assessment. Corrections were filled using linear interpolation and/or sensor offsets following Quality Assurance and Quality Control (QA/QC) protocols detailed in the NBFSMN’s Quality Assurance Project Plan.

The Conanicut Hydrocat 620 and Potowomut Hydrocat 720 sensors were equipped with Seabird Hydrocat-EP CTD sensors (model number# HC-EP.1011S), equipped with a Seabird Hydrocat-EP pH Module (model#802532). Hydrocat-EPs collected multiple parameters including temperature, salinity, pH, dissolved oxygen. The hydrocat sensors were deployed year-round and were recovered every 2 to 4 weeks, weather permitting, for maintenance at which point the flow path and the conductivity cell are flushed with 1% Triton detergent followed by vigorous flushing with DI water. Following the cleaning, conductivity checks with conductivity standards (1413 μS,10,000 μS) ensures proper function, as well as a three point pH calibration of the pH Module (pH 4,7,10). pH calibrations were performed in Seabird's UCI software using their pH calibration wizard. Sensors were sent back to Seabird Scientific for yearly manufacturer calibrations. QA/QC was conducted in R using the R package OCE (Kelley & Richards, 2024), following QA/QC tests established by OOIO (NSF Ocean Observatories Initiative, 2012) including, a stuck value test, despiking, and a range test against global values for measured parameters provided by OOIO and local values in Narragansett Bay provided by the Narragansett Bay Long-Term Plankton Time Series (https://web.uri.edu/gso/research/plankton/).

DIC and TA in discrete samples were measured according to Dickson et al. (2007) using the Apollo SciTech Model AS-C6L Dissolved Inorganic Carbon Analyzer and the Apollo SciTech Model AS-ALK3 Total Alkalinity Titrator. Instruments were calibrated to Certified Reference Material (CRM) from Scripps Institute of Oceanography at room temperature (21 - 22° C). For DIC analysis, a sample of CRM was run prior to and at the end of sample analysis for quality control. TA analysis was calibrated to CRM, and either a sample of CRM or aged open ocean water was run at beginning of daily analysis and end of daily analysis for quality control. Lab-based measurements carry a ±0.2% uncertainty for alkalinity and ±0.1% uncertainty for DIC.

Temperature and salinity in the discrete carbonate data were measured with the Seabird instruments (analogous sites in the autonomous dataset).  Seabird Hydrocat-EP CTD sensors (model number# HC-EP.1011S), equipped with a Seabird Hydrocat-EP pH Module (model#802532). Hydrocat-EPs collected multiple parameters, including temperature, salinity, pH, and dissolved oxygen. 


Data Processing Description

This dataset contains discrete carbonate data (TA, DIC, pH, Temperature, Salinity). See the "Related Publications" section for autonomously collected pH data also described in this methodology.

We performed a gross range test on sensor data to eliminate unrealistic data (i.e. pH measurements less than 5 or greater than 10, or DO less than 1 mg L-1) that indicate biofouling or other sensor malfunction. We then removed obvious outliers--any data points that were more than 3 standard deviations greater or less than the sensor's annual mean or any data points that were more than 1.5 standard deviations greater than or less than the sensor's 24-hour moving mean. 

Autonomous observations were compared to discrete samples. Occasionally, periods of sensor data diverged from bottle sample data with consistent, identifiable bias--for instance, a consistent underestimation or overestimation of pH by a value less than 0.1 pH unit or a dynamic bias less than 0.1 pH unit that changes with temperature or time. In the case of QP, we corrected this by adjusting all the sensor data points by the mean difference between sensor pH and bottle pH; specifically, we lowered QP pH by 0.084, the mean amount by which QP overestimated pH. Conanicut Hydrocat 620 overestimated pH in January-March 2023 when temperatures are lowest. We applied a more dynamic correction as a third order polynomial function of temperature, such that sensor pH was adjusted to be lower at colder temperatures.  Conanicut Hydrocat 620 and QP sensors were compared to weekly bottle samples collected as part of the Narragansett Bay Long Term Phytoplankton Time Series, whereas Potowomut Hydrocat 720 and MV sensors were compared to bottle samples collected approximately monthly at the same location at the Potowomut sensor.

pH measurements from all sensors were originally on NBS scale but converted to total scale using the Python PyCO2SYS package (Humphreys et al., 2022).


BCO-DMO Processing Description

* The data table within the submitted file "DiscreteData_Baskind2024.csv" was imported into the BCO-DMO data system for this dataset.  Table will appear as Data File: 961940_v1_discrete_carbonate-narr-bay.csv (along with other download format options).

* Column names adjusted to conform to BCO-DMO naming conventions designed to support broad re-use by a variety of research tools and scripting languages. [Only numbers, letters, and underscores.  Can not start with a number]

* BCO-DMO requires all columns have names. After a conversation with data submitter, the un-named first column with a row id was dropped from the table.

* Additional site_lat, site_lon columns were added to the data table from site information included in provided metadata for the locations.


[ table of contents | back to top ]

Data Files

File
961940_v1_discrete_carbonate-narr-bay.csv
(Comma Separated Values (.csv), 12.33 KB)
MD5:3a7831cf0d40ab4089a02a0203cefd28
Primary data file for dataset ID 961940, version 1

[ table of contents | back to top ]

Related Publications

Baskind A, Ahumada G, Gomes K, Stoffel H, Gu S, Davies A and Wang H (2025) Unraveling natural carbonate variability in Narragansett Bay, RI using multiple high temporal resolution pH time series. Front. Mar. Sci. 12:1552350. doi: 10.3389/fmars.2025.1552350
Results
Dickson, A.G.; Sabine, C.L. and Christian, J.R. (eds) (2007) Guide to best practices for ocean CO2 measurement. Sidney, British Columbia, North Pacific Marine Science Organization, 191pp. (PICES Special Publication 3; IOCCP Report 8). DOI: https://doi.org/10.25607/OBP-1342
Methods
Humphreys, M. P., Lewis, E. R., Sharp, J. D., & Pierrot, D. (2022). PyCO2SYS v1.8: marine carbonate system calculations in Python. Geoscientific Model Development, 15(1), 15–43. https://doi.org/10.5194/gmd-15-15-2022
Software
Kelley, D. E., Richards, C., & Layton, C. (2022). oce: an R package for Oceanographic Analysis. Journal of Open Source Software, 7(71), 3594. https://doi.org/10.21105/joss.03594
Software
Kelley, D., & Richards, C. (2007). oce: Analysis of Oceanographic Data [dataset]. In CRAN: Contributed Packages. The R Foundation. https://doi.org/10.32614/cran.package.oce https://doi.org/10.32614/CRAN.package.oce
Software
Kelley, D., & Richards, C. (2024). oce: Analysis of oceanographic data (Version 1.8-4) [R package]. Retrieved from https://dankelley.github.io/oce/
Software
McDougall, T. J., & Barker, P. M. (2011). Getting started with TEOS-10 and the Gibbs Seawater (GSW) oceanographic toolbox. Scor/iapso WG, 127(532), 1-28. ISBN: 978-0-646-55621-5
Methods
NSF Ocean Observatories Initiative. (2012). Data product specification for pH of seawater (Version 1-00). Retrieved from https://oceanobservatories.org/wp-content/uploads/2015/10/1341-00510_Data_Product_Spec_PHWATER_OOI.pdf
Methods
RIDEM. (2020). Quality assurance project plan: Narragansett Bay fixed site water quality monitoring network seasonal monitoring. Rhode Island Department of Environmental Managament. Retrieved from https://dem.ri.gov/sites/g/files/xkgbur861/files/programs/benviron/water/quality/surfwq/pdfs/nbfsmn.pdf
Methods
TEOS-10 developers (2017). GSW-Python — gsw 3.6.16.post1 documentation. (n.d.-c). Retrieved from https://teos-10.github.io/GSW-Python/
Software
University of Rhode Island (2025). The Narragansett Bay Long-Term Plankton Time Series. Retrieved from https://web.uri.edu/gso/research/plankton/
Methods

[ table of contents | back to top ]

Related Datasets

IsRelatedTo
Baskind, A., Ahumada, G., Gomes, K., Stoffel, H., Gu, S., Davies, A., Wang, H. (2025) Autonomously collected pH data from Narragansett Bay in 2022 and 2023. Biological and Chemical Oceanography Data Management Office (BCO-DMO). (Version 1) Version Date 2025-05-14 doi:10.26008/1912/bco-dmo.961920.1 [view at BCO-DMO]
Relationship Description: Data collected as part of the same study to be published in Baskind et al. (2025, doi: 10.3389/fmars.2025.1552350).

[ table of contents | back to top ]

Parameters

ParameterDescriptionUnits
Sample

The name of the sample and bottle number

unitless
ISO_DateTime_UTC

Date and time of measurement in UTC time zone

unitless
Location

Sample location (PLT for Narragansett Bay Long Term Phytoplankton Time Series site and GB for Greenwich Bay)

unitless
site_lat

Sample location latitude

decimal degrees
site_lon

Sample location longitude

decimal degrees
depth_category

depth of sample (surface or bottom)

unitless
Salinity

In situ salinity

Practical Salinity Units (PSU)
In_Situ_Temperature

In situ temperature

degrees Celsius (degC)
DIC

measured dissolved inorganic carbon

micromoles per kilogram (umol/kg)
TA

measured total alkalinity

micromoles per kilogram (umol/kg)
pH_Total

pH in total scale calculated from DIC, TA, temperature, and salinity using PyCO2SYS

total pH scale
pH_NBS

pH in NBS scale calculated from DIC, TA, temperature, and salinity using PyCO2SYS

NBS scale
pH_norm_12pt5deg

pH in total scale calculated from DIC, TALK, temperature, and salinity using PyCO2SYS and normalized to 12.5°C, the approximate annual mean temperature of Narragansett Bay

total pH scale


[ table of contents | back to top ]

Instruments

Dataset-specific Instrument Name
Generic Instrument Name
Apollo SciTech AS-C6L Dissolved Inorganic Carbon (DIC) analyzer
Dataset-specific Description
DIC and TA in discrete samples were measured according to Dickson et al. (2007) using the Apollo SciTech Model AS-C6L Dissolved Inorganic Carbon Analyzer and the Apollo SciTech Model AS-ALK3 Total Alkalinity Titrator. 
Generic Instrument Description
An instrument designed for the analysis of dissolved inorganic carbon in samples from various aquatic environments. It comprises of a laser-based CO2 detector (LI-7815), a digital syringe pump, a mass flow controller, CO2 stripping reactor, an electronic cooling system and a computer communication assembly (RS-485, USB). The AS-C6L supersedes the earlier AS-C3 model, which used non-dispersive infra-red CO2 detection (LI-7000, discontinued). The AS-C6L improves on the AS-C3 by incorporating a multi-sampler of one set of standards plus 8 samples, and uses improved Apollo SciTech software. The AS-C6L is suitable for use in either shipboard or land-based laboratories. It maintains a precision of +/-0.1 % for seawater (or +/-2 umol/kg), enables sample volumes ranging from 0.5 - 3.5 ml per analysis, and an analytical rate of approximately 3 minutes.

Dataset-specific Instrument Name
Generic Instrument Name
Apollo SciTech Model AS-ALK3 total alkalinity titrator
Dataset-specific Description
DIC and TA in discrete samples were measured according to Dickson et al. (2007) using the Apollo SciTech Model AS-C6L Dissolved Inorganic Carbon Analyzer and the Apollo SciTech Model AS-ALK3 Total Alkalinity Titrator.
Generic Instrument Description
An automated acid-base titrator for use in aquatic carbon dioxide parameter analysis. The titrator provides standardisation and sample analysis, using the Gran titration procedure for alkalinity determination of seawater and brackish waters. It is designed for both shipboard and land based laboratory use. The precision of the instrument is 0.1 percent or higher, and sample volumes may range from 10-25 ml. Titraton takes approximately 8 minutes per sample, and the repeatability is within plus or minus 1-2 micromoles per kg.

Dataset-specific Instrument Name
Seabird Hydrocat-EP CTD sensors (model number# HC-EP.1011S)
Generic Instrument Name
CTD Sea-Bird
Dataset-specific Description
Temperature and salinity in the discrete carbonate data were measured with the Seabird instruments (analogous sites in the autonomous dataset).  Seabird Hydrocat-EP CTD sensors (model number# HC-EP.1011S), equipped with a Seabird Hydrocat-EP pH Module (model#802532). Hydrocat-EPs collected multiple parameters, including temperature, salinity, pH, dissolved oxygen. 
Generic Instrument Description
A Conductivity, Temperature, Depth (CTD) sensor package from SeaBird Electronics. This instrument designation is used when specific make and model are not known or when a more specific term is not available in the BCO-DMO vocabulary. Refer to the dataset-specific metadata for more information about the specific CTD used. More information from: http://www.seabird.com/

Dataset-specific Instrument Name
Seabird Hydrocat-EP pH Module (model#802532)
Generic Instrument Name
Water Quality Multiprobe
Dataset-specific Description
Temperature and salinity in the discrete carbonate data were measured with the Seabird instruments (analogous sites in the autonomous dataset).  Seabird Hydrocat-EP CTD sensors (model number# HC-EP.1011S), equipped with a Seabird Hydrocat-EP pH Module (model#802532). Hydrocat-EPs collected multiple parameters including temperature, salinity, pH, dissolved oxygen. 
Generic Instrument Description
An instrument which measures multiple water quality parameters based on the sensor configuration.


[ table of contents | back to top ]

Project Information

Drivers of ocean acidification in a temperate urbanized estuary undergoing nutrient loading reductions (OA Drivers in an Urbanized Estuary)

Coverage: Narragansett Bay, Rhode Island


NSF Award Abstract:
About one-third of human-produced carbon dioxide (CO2) has dissolved in the ocean, significantly slowing down climate change and global warming. However, as more CO2 dissolves in the sea, the ocean water becomes acidified because CO2 is a weak acid. We call this change ocean acidification (OA). OA has a significant and lasting impact on marine life, as acidified water causes shell dissolution for numerous marine species. US society is already grappling with the adverse effects of OA. OA is predicted to lead to economic losses of $400 million annually by 2100 for the US shellfish industry. Therefore, it is urgent to understand how OA is changing and the potential approaches to managing it in coastal waters. Presently, we do not have a good understanding of OA change at either the large or the small scale because so many competing processes impact OA.

One of the most dramatic processes that affect ocean acidification (OA) is the extra nutrient loading from rivers due to human inputs. While moderate nutrient levels are necessary to support marine life, nutrient levels that are too high can cause oceanic plankton blooms. As plankton die, they sink to the seafloor, where they are decomposed by bacteria and speed up OA where shellfish live. In recent decades, management agencies have been working on improving water quality by reducing the nutrient input to coastal waters, which has the potential to reduce the OA even though OA is not the initial motivation. To our knowledge, there has yet to be targeted research to understand how OA responds to aggressive nutrient reduction. The average nutrient input to Narragansett Bay? one of the most rapidly warming estuaries in the US?has been reduced over one-third in recent decades because of intentional nutrient reduction. Therefore, Narragansett Bay can be considered a natural laboratory that can help in better understanding the OA change that other coastal areas may soon experience resulting from environmental management actions and climate change. This project will leverage the existing water quality monitoring network, collect new data, and utilize a coastal biogeochemical model to quantify OA change and the mechanisms driving its change in Narragansett Bay, RI. Through engagement with the public and local shellfish growers, we will increase their awareness of OA and enhance their ability to carry out ongoing or new aquaculture operations in the near and long term. As more US states reduce nutrient input to improve coastal water quality, this project will allow management agencies nationwide to make next-generation refinements to best limit OA for sustainable fisheries and wildlife in the face of climate change.

This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria



[ table of contents | back to top ]

Funding

Funding SourceAward
NSF Division of Ocean Sciences (NSF OCE)

[ table of contents | back to top ]