National Oceanic and
Atmospheric Administration
United States Department of Commerce


FY 2023

A multi-stressor data product for marine heatwave, hypoxia, and ocean acidification research, including calculated inorganic carbon parameters from the southern Salish Sea and northern California Current System from 2008-02-04 to 2018-10-19

Alin, S.R., J. Newton, R.A. Feely, D. Greeley, J. Herndon, and A. Kozyr

NCEI Accession 0283266, NOAA National Centers for Environmental Information, Ocean Carbon and Acidification Data System (OCADS),, doi: 10.25921/5g29-q841, View dataset online (2023)

This NCEI Accession contains a data product including calculated CO2 system parameters based on the Salish cruise data product ( We used R seacarb function carb to calculate the most commonly used derived carbonate system parameters, including pH, partial pressures and fugacities of carbon dioxide at in situ temperatures and pressures (pCO2insitu, and fCO2insitu, respectively), and aragonite and calcite saturation states (OmegaAragonite and OmegaCalcite, respectively) (Gattuso et al., 2023). Input parameters from the Salish cruise compiled data set (Alin et al. 2021) comprised dissolved inorganic carbon (DIC_UMOL_KG), total alkalinity (TA_UMOL_KG), phosphate (PHOSPHATE_UMOL_KG), and silicate (SILICATE_UMOL_KG) content values from bottle samples analyzed in the laboratory, along with CTD measurements of temperature (CTDTMP_DEG_C_ITS90), salinity (CTDSAL_PSS78), and pressure (CTDPRS_DBAR). Within seacarb, we used the TEOS-10 thermodynamic seawater equations (IOC, SCOR, and IAPSO, 2010). We adopted the total scale for pH (pHT), the Uppstrom (1974) formulation for deriving total boron concentration from salinity, the seacarb default option for Kf (Perez and Fraga, 1987 for temperatures above 9 °C; Dickson and Goyet, 1994 for those below), and the Dickson (Dickson, 1990) option for Ks (following results of Orr et al., 2015). All input content data were first divided by 106 to convert from µmol kg^–1 to mol kg^–1, and pressure (dbar) was divided by 10 to convert to bar, to conform with the default units of seacarb. For equilibrium constants (K1 and K2), we provide calculated values using both the Lueker et al. (2000) and the Waters et al. (2014) dissociation constants. The Lueker constants (for salinity ranges of 19–43 and temperature ranges of 2–35°C) facilitate comparison with publications arising from West Coast Ocean Acidification (WCOA) cruise data sets (, whereas the Waters constants (for salinity ranges of 1–50 and temperature ranges of 0–50°C) allow users working with more brackish salinities to compare their results directly to those in the Salish cruise data product. All references above are included in the seacarb documentation. Alin et al. (2023b) describe the magnitude of differences in calculated values for the Salish cruise data product using the two different sets of dissociation constants. The data product is available on the Index page accessed by clicking the Database Files link above: The file “SalishCruises_2008to2018_MeasCalcParams_NCEIdataProduct_09262023.csv” contains 3971 complete records of DIC, TA, T, S, O2, and nutrient measurements with the highest quality QC flags and includes calculated values for the carbonate system parameters described above.This effort was conducted in support of the estuarine and coastal monitoring and research objectives of the Washington Ocean Acidification Center (WOAC), the Northwest Association of Networked Ocean Observing Systems, the U.S. National Oceanic and Atmospheric Administration's Pacific Marine Environmental Laboratory, and the U.S. National Oceanic and Atmospheric Administration's Ocean Acidification Program and conforms to climate-quality monitoring guidelines of the Global Ocean Acidification Observing Network ( For any questions about appropriate use or limitations of the data set, please contact Drs. Simone Alin and Jan Newton at email addresses above.

Feature Publications | Outstanding Scientific Publications

Contact Sandra Bigley |