The Surface Ocean CO₂ Atlas (SOCAT) is a synthesis of quality-controlled fCO₂ (fugacity of carbon dioxide) values for the global surface oceans and coastal seas with regular updates. Version 3 of SOCAT has 14.7 million fCO₂ values from 3646 data sets covering the years 1957 to 2014. This latest version has an additional 4.6 million fCO₂ values relative to version 2 and extends the record from 2011 to 2014. Version 3 also significantly increases the data availability for 2005 to 2013. SOCAT has an average of approximately 1.2 million surface water fCO₂ values per year for the years 2006 to 2012. Quality and documentation of the data has improved. A new feature is the data set quality control (QC) flag of E for data from alternative sensors and platforms. The accuracy of surface water fCO₂ has been defined for all data set QC flags. Automated range checking has been carried out for all data sets during their upload into SOCAT. The upgrade of the interactive Data Set Viewer (previously known as the Cruise Data Viewer) allows better interrogation of the SOCAT data collection and rapid creation of high-quality figures for scientific presentations. Automated data upload has been launched for version 4 and will enable more frequent SOCAT releases in the future. High-profile scientific applications of SOCAT include quantification of the ocean sink for atmospheric carbon dioxide and its long-term variation, detection of ocean acidification, as well as evaluation of coupled-climate and ocean-only biogeochemical models. Users of SOCAT data products are urged to acknowledge the contribution of data providers, as stated in the SOCAT Fair Data Use Statement. This ESSD (Earth System Science Data) "living data" publication documents the methods and data sets used for the assembly of this new version of the SOCAT data collection and compares these with those used for earlier versions of the data collection (Pfeil et al., 2013; Sabine et al., 2013; Bakker et al., 2014). Individual data set files, included in the synthesis product, can be downloaded here: doi:10.1594/PANGAEA.849770. The gridded products are available here: doi:10.3334/CDIAC/OTG.SOCAT_V3_GRID.

Oceanic dissolved inorganic carbon (T_C) is the largest pool of carbon that substantially interacts with the atmosphere on human timescales. Oceanic T_C is increasing through uptake of anthropogenic carbon dioxide (CO₂), and seawater pH is decreasing as a consequence. Both the exchange of CO₂ between the ocean and atmosphere and the pH response are governed by a set of parameters that interact through chemical equilibria, collectively known as the marine carbonate system. To investigate these processes, at least two of the marine carbonate system's parameters are typically measured – most commonly, two from T_C, total alkalinity (A_T), pH, and seawater CO₂ fugacity ($f_{{\mathrm{CO}}_{\mathrm{2}}}$ ; or its partial pressure, $p_{{\mathrm{CO}}_{\mathrm{2}}}$ , or its dry-air mole fraction, $x_{{\mathrm{CO}}_{\mathrm{2}}}$ ) – from which the remaining parameters can be calculated and the equilibrium state of seawater solved. Several software tools exist to carry out these calculations, but no fully functional and rigorously validated tool written in Python, a popular scientific programming language, was previously available. Here, we present PyCO2SYS, a Python package intended to fill this capability gap. We describe the elements of PyCO2SYS that have been inherited from the existing CO2SYS family of software and explain subsequent adjustments and improvements. For example, PyCO2SYS uses automatic differentiation to solve the marine carbonate system and calculate chemical buffer factors, ensuring that the effect of every modelled solute and reaction is accurately included in all its results. We validate PyCO2SYS with internal consistency tests and comparisons against other software, showing that PyCO2SYS produces results that are either virtually identical or different for known reasons, with the differences negligible for all practical purposes. We discuss insights that guided the development of PyCO2SYS: for example, the fact that the marine carbonate system cannot be unambiguously solved from certain pairs of parameters. Finally, we consider potential future developments to PyCO2SYS and discuss the outlook for this and other software for solving the marine carbonate system. The code for PyCO2SYS is distributed via GitHub (https://github.com/mvdh7/PyCO2SYS, last access: 23 December 2021) under the GNU General Public License v3, archived on Zenodo (Humphreys et al., 2021), and documented online (https://pyco2sys.readthedocs.io/en/latest/, last access: 23 December 2021).

The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2022 is an update of the previous version, GLODAPv2.2021 (Lauvset et al., 2021). The major changes are as follows: data from 96 new cruises were added, data coverage was extended until 2021, and for the first time we performed secondary quality control on all sulfur hexafluoride (SF₆) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF₆ data, which are now subjected to secondary quality control, and carbon data measured on board the RV Knorr in the Indian Ocean in 1994–1995 which are now adjusted using certified reference material (CRM) measurements made at the time. GLODAPv2.2022 includes measurements from almost 1.4 million water samples from the global oceans collected on 1085 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, chlorofluorocarbon-11 (CFC-11), CFC-12, CFC-113, CCl₄, and SF₆) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 96 new cruises were derived by comparing those data with the data from the 989 quality-controlled cruises in the GLODAPv2.2021 data product using crossover analysis. SF₆ data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO₂) chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg⁻¹ in dissolved inorganic carbon, 4 µmol kg⁻¹ in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO₂ fugacity (fCO₂), were not subjected to bias comparison or adjustments. The original data, their documentation, and DOI codes are available at the Ocean Carbon and Acidification Data System of NOAA NCEI (https://www.ncei.noaa.gov/access/ocean-carbon-acidification-data-system/oceans/GLODAPv2_2022/, last access: 15 August 2022). This site also provides access to the merged data product, which is provided as a single global file and as four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under https://doi.org/10.25921/1f4w-0t92 (Lauvset et al., 2022). These bias-adjusted product files also include significant ancillary and approximated data, which were obtained by interpolation of, or calculation from, measured data. This living data update documents the GLODAPv2.2022 methods and provides a broad overview of the secondary quality control procedures and results.

The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface to bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2023 is an update of the previous version, GLODAPv2.2022 (Lauvset et al., 2022). The major changes are as follows: data from 23 new cruises were added. In addition, a number of changes were made to the data included in GLODAPv2.2022. GLODAPv2.2023 includes measurements from more than 1.4 million water samples from the global oceans collected on 1108 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, chlorofluorocarbon-11 (CFC-11), CFC-12, CFC-113, CCl₄, and SF₆) have undergone extensive quality control with a focus on the systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 23 new cruises were derived by comparing those data with the data from the 1085 quality-controlled cruises in the GLODAPv2.2022 data product using crossover analysis. SF₆ data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO₂), chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data-handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg⁻¹ in dissolved inorganic carbon, 4 µmol kg⁻¹ in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO₂ fugacity (fCO₂), were not subjected to bias comparison or adjustments. The original data, their documentation, and DOI codes are available at the Ocean Carbon and Acidification Data System of NOAA National Centers for Environmental Information (NCEI), which also provides access to the merged data product. This is provided as a single global file and as four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under https://doi.org/10.25921/zyrq-ht66 (Lauvset et al., 2023). These bias-adjusted product files also include significant ancillary and approximated data, which were obtained by interpolation of, or calculation from, measured data. This living data update documents the GLODAPv2.2023 methods and provides a broad overview of the secondary quality control procedures and results.

The ocean holds a large reservoir of carbon dioxide (CO2), and mitigates climate change through uptake of anthropogenic CO2. Fluxes of CO2 between the atmosphere and surface ocean are regulated by a number of physical and biogeochemical processes, resulting in a spatiotemporally heterogeneous CO2 distribution. Determining the influence of each individual process is useful for interpreting marine carbonate system observations, and is also necessary to investigate how changes in these drivers could affect air-sea CO2 exchange. Biogeochemical processes exert an influence primarily through modifying seawater dissolved inorganic carbon (CT) and total alkalinity (AT), thus changing the seawater partial pressure of CO2 (psw). Here, we propose a novel conceptual framework through which the size of the CO2 source or sink generated by any biogeochemical process, denoted Φ, can be evaluated. This is based on the ‘isocapnic quotient’ (Q), which defines the trajectory through (AT,CT) phase space for which there is no change in psw. We discuss the limitations and uncertainties inherent in this technique, which are negligible for most practical purposes, and its links with existing, related approaches. We investigate the effect on Φ of spatiotemporal heterogeneity in Q in the present day surface ocean for several key biogeochemical processes. This leads the magnitude of the CO2 source or sink generated by processes that modify AT to vary spatiotemporally. Finally, we consider how the strength of each process as a CO2 source or sink may change in a warmer, higher-CO2 future ocean