This data package contains a hybrid surface OA data product that is produced based on three recent observational data products: (a) the Surface Ocean CO2 Atlas (SOCAT, version 2022), (b) the Global Ocean Data Analysis Product version 2 (GLODAPv2, version 2022), and (c) the Coastal Ocean Data Analysis Product in North America (CODAP-NA, version 2021), and 14 Earth System Models from the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The trajectories of ten OA indicators, including fugacity of carbon dioxide, pH on Total Scale, total hydrogen ion content, free hydrogen ion content, carbonate ion content, aragonite saturation state, calcite saturation state, Revelle Factor, total dissolved inorganic carbon content, and total alkalinity content are provided under preindustrial conditions, historical conditions, and future Shared Socioeconomic Pathways: SSP1-19, SSP1-26, SSP2-45, SSP3-70, and SSP5-85 from 1750 to 2100 on a global surface ocean grid. These OA trajectories are improved relative to previous OA data products with respect to data quantity, spatial and temporal coverage, diversity of the underlying data and model simulations, and the provided SSPs over the 21st century.

Dataset Description: Gridded monthly data products of surface ocean acidification indicators from 1998 to 2024 and on a 0.25° by 0.25° spatial grid have been developed for eleven U.S. Large Marine Ecosystems (LMEs) using a machine learning algorithm called random forest regression (RFR). The data products are called RFR-LMEs, and were constructed using observations from the Surface Ocean CO2 Atlas — co-located with surface ocean properties from various satellite, reanalysis, and observational products — with an approach that utilized two types of machine learning algorithms: (1) Gaussian mixture models to cluster the data into subregions with similar environmental variability and (2) RFRs that were trained and applied separately in each cluster to interpolate the observational data in space and time. RFR-LMEs also rely on previously published seawater property estimation routines to obtain the full suite of ocean acidification indicators. The products show a domain-wide carbo n dioxide partial pressure increase of 1.6 ± 0.4 μatm yr−1 and pH decrease of 0.0015 ± 0.0004 yr−1. More information on the creation and validation of RFR-LMEs is available in the following publication: Sharp, J.D., Jiang, L., Carter, B.R., Lavin, P.D., Yoo, H., Cross, S.L., 2024. A mapped dataset of surface ocean acidification indicators in large marine ecosystems of the United States. Scientific Data, 11, 715, 10.1038/s41597-024- 03530-7.

The dataset contains sea-surface temperature, salinity, dissolved inorganic carbon, total alkalinity, pH, and partial pressure of CO2 for the Indian Ocean region (Longitude: 30°E-120°E, Latitude: 30°S-30°N). The data is available from 1980 to 2019 on a monthly time scale. Each of these data variables has a spatial resolution of 1/12°.

This data package contains 1x1 degree coastal climatologies (mean fields of oceanographic variables on a regular geographic grid at specific depths) of pH on the Total Scale, total hydrogen ion content, free hydrogen ion content, carbonate ion content, aragonite saturation state, calcite saturation state, total dissolved inorganic carbon content, and total alkalinity content from surface to 500 meters water depth on North American ocean margins. These climatologies were developed with the World Ocean Atlas (WOA) gridding technologies of the NOAA National Centers for Environmental Information (NCEI), based on the recently released Coastal Ocean Data Analysis Product in North America (CODAP-NA, https://doi.org/10.5194/essd-13-2777-2021), along with the Global Ocean Data Analysis Product version 2 (GLODAPv2, version 2021, https://doi.org/10.5194/essd-13-5565-2021). The relevant variables were adjusted to the year of 2010 based on the algorithms as developed by Carter et al. (2021, https://doi.org/10.1002/lom3.10461) before the gridding. This data package contains a total of 8 NetCDF files, one for each of the variable. It is recommended to use the objectively analyzed mean fields (with "_an" postfix) for each variable.

This dataset consists of year 2002 and preindustrial (pi) OA metric values and their uncertainties (u), total OA metric changes (d) due to anthropogenic carbon accumulation to the year 2002 and the component of those changes caused by carbonate system nonlinearities (n), with associated uncertainties provided. Uncertainties were estimated using a 1000 iteration Monte Carlo simulation. Data from the upper 2000 m of the GLODAPv2.2016b mapped data product (https://doi.org/10.3334/cdiac/otg.ndp093_glodapv2), described in Lauvset et al., 2016 (https://doi.org/10.5194/essd-8-325-2016), and from the preformed properties product of Carter et al., 2021 (https://doi.org/10.1029/2020GB006623) were used to make these estimates. Calculation details are described in Fassbender et al., 2023 (https://doi.org/10.1029/2023GB007843). Year 2002 aragonite saturation state and pH values, and their uncertainties, are reproduced from the GLODAPv2.2016b mapped data product (https://doi.org/10.7289/v5kw5d97) and are provided here for user convenience with the permission of the original data producer. Version 1.0.

This dataset includes high-frequency (hourly to sub-hourly) coastal acidification time-series data collected during nine deployments in the aforementioned seven estuaries along the US East Coast, US West Coast and Gulf of Mexico from 2015-04-23 to 2020-07-29. These data include water temperature, salinity, partial pressure of carbon dioxide (pCO2) in water, dissolved oxygen (DO) in water, and pH on the total scale. The instruments used to collected these data include Sunburst SAMI-CO2, Pro-Oceanus CO2-Pro CV and a LiCOr LI-820 CO2 gas analyzers for autonomous pCO2 measurements, Sea-Bird SeapHOx and SeaFET instruments for pH measurements, Sea-Bird SeapHOx and Aanderaa Oxygen Optode instruments for DO measurements, and YSI water sensing instrument packages for measurements of conductivity (salinity), temperature and depth. Beginning in 2015, the U.S. Environmental Protection Agency’s (EPA) National Estuary Program (NEP) started a collaboration with partners in seven estuaries along the East Coast (Barnegat Bay; Casco Bay), West Coast (Santa Monica Bay; San Francisco Bay; Tillamook Bay), and the Gulf of Mexico (GOM) Coast (Tampa Bay; Mission-Aransas Estuary) of the United States to expand the use of autonomous monitoring partial pressure of carbon dioxide (pCO2) and pH sensors to evaluate carbonate chemistry in the estuarine environment.

This dataset contains spatial distribution of surface ocean pH (total hydrogen scale), acidity (or hydrogen ion activity, unit: nmol/kg, or 10^-9 mol/kg) and Revelle Factor (a measure of the ocean's buffer capacity, unitless) on a 1x1 degree global grid (Longitude: [20.5:1:379.5], Latitude: [-89.5:1:89.5]) in all 12 months of the years from 1770 to 2100 (1770, 1780, 1790, ..., 2100). This data product is produced by combining a recent observational seawater carbon dioxide (CO2) data product, i.e., the 6th version of the Surface Ocean CO2 Atlas (1991-2018, ~23 million observations), with temporal trends at individual locations of the global ocean from a robust Earth System Model (ESM2M), to provide a high-resolution regionally varying view of global surface ocean pH, acidity, and the Revelle Factor. The climatology extends from the pre-Industrial era (1770 C.E.) to the end of this century under historical atmospheric CO2 concentrations (pre-2005) and the Representative Concentrations Pathways (RCP2.6, RCP4.5, RCP6.0 and RCP8.5, post-2005) of the Intergovernmental Panel on Climate Change (IPCC)’s 5th Assessment Report (AR5). By linking the modeled pH trends to the observed modern pH distribution, the climatology benefits from recent improvements in both model design and observational data coverage, and is likely to provide improved regional OA trajectories than the model output could alone, therefore, will help guide the regional OA adaptation strategies. Revelle Factor is defined as the ratio between the fractional change in pCO2 to the fractional change in dissolved inorganic carbon (DIC). This dataset is available in netCDF format. Some plots and animation files are also available for your presentation purposes. For details of the calculation and gridding method, please refer to Jiang, L.-Q., B. R. Carter, R. A. Feely, S. Lauvset, and A. Olsen (2019), Surface ocean pH and buffer capacity: past, present and future, Nature Scientific Reports, 9:18624, https://doi.org/10.1038/s41598-019-55039-4.

This dataset contains discrete measured and calculated surface seawater properties from samples taken on the first week of every month over the course of one or two days between 11/2021 and 07/2022 (9 months) at the Equatorial Pacific region off the coast of Colombia, around Gorgona Island National Natural Park. The dataset contains also measured monthly sea-surface temperature (SST), sea-surface salinity (SSS), pH in total scale (derived directly from in situ pH measurements) with depth of 1.8 metres along 7 sampling sites in a coast-ocean gradient (1 to 7) with the following geographical coverage: 1. Lat. -78.201138 Lon. 3.0163611; 2. Lat. -78.221416 Lon. 2.9705833; 3. Lat. -78.163861 Lon. 2.9960278; 4. Lat. -78.101555 Lon. 2.9473889; 5. Lat. -78.101555 Lon. 2.9322222; 6. Lat. -78.122865 Lon. 2.8756667; 7. Lat. -77.979865 Lon. 2.7204915. The measurements were taken with two CTD divers (SonTek CastAway-CTD and Eijkelkamp CTD-Diver 11.11.60) and one Niskin bottle (5L). The water samples where taken to a oceanography lab in the Pontifical Xavierian University to be tested for total alkalinity (TA), and to the Autonomous University of Baja California for Dissolved Inorganic Carbon (DIC) measurements. This was done in order to calculate the carbon fluxes on a monthly scale throughout nine months during the 2021-2022 La Niña year for the whole study site, allowing for observations in the spatiotemporal flux behavior variability (shown in dataset as fco2_measured).

This dataset consists of time-resolved reconstructions of ocean interior acidification from 1800 through 1994, 2004, and 2014. The basis of these reconstructions are observation-based estimates of the accumulation of anthropogenic carbon, combined with climatologies of hydrographic and biogeochemical properties in the ocean interior. Acidification trends are determined for several parameters of the marine CO2 system, namely the saturation state of aragonite (Ωarag), the carbonate ion concentration ([CO32-]), the free proton concentration ([H+]), and pH on the total scale (pHT). The underlying anthropogenic carbon concentration (ΔCant), the computed sensitivities of the four marine CO2 system parameters and their absolute state estimates are provided as well. The datasets contain in addition to the standard estimate also 14 sensitivity cases, which are intended to assess the robustness of our acidification estimates to changes in the estimation procedure of ΔCant as well as the climatological distributions of other hydrographic properties. All estimates are provided on a horizontal grid with 1° x 1° resolution and for 28 depth layers from 0 - 3000m. These data provide strong constraints on ocean interior acidification over the industrial era, unravelling in particular its progression since 1994.