Search

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 (SF6) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF6 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, CCl4, and SF6) 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. SF6 data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO2) 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−1 in dissolved inorganic carbon, 4 µmol kg−1 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 CO2 fugacity (fCO2), 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.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the first time, an approach is shown to reconcile the difference in our ELUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of 10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for 2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLAND was 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. The global atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2020, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 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.2021 is an update of the previous version, GLODAPv2.2020 (Olsen et al., 2020). The major changes are as follows: data from 43 new cruises were added, data coverage was extended until 2020, all data with missing temperatures were removed, and a digital object identifier (DOI) was included for each cruise in the product files. In addition, a number of minor corrections to GLODAPv2.2020 data were performed. GLODAPv2.2021 includes measurements from more than 1.3 million water samples from the global oceans collected on 989 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) 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 updated to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For this annual update, adjustments for the 43 new cruises were derived by comparing those data with the data from the 946 quality controlled cruises in the GLODAPv2.2020 data product using crossover analysis. Comparisons to estimates of nutrients and ocean CO2 chemistry based on empirical algorithms provided additional context for adjustment decisions in this version. The adjustments 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 with to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 4 µmol kg−1 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 CO2 fugacity (fCO2), were not subjected to bias comparison or adjustments. The original data, their documentation, and DOI codes are available at the Ocean Carbon Data System of NOAA NCEI (https://www.ncei.noaa.gov/access/ocean-carbon-data-system/oceans/GLODAPv2_2021/, last access: 7 July 2021). 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/ttgq-n825 (Lauvset et al., 2021). These bias-adjusted product files also include significant ancillary and approximated data and can be accessed via https://www.glodap.info (last access: 29 June 2021). These were obtained by interpolation of, or calculation from, measured data. This living data update documents the GLODAPv2.2021 methods and provides a broad overview of the secondary quality control procedures and results.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ±  0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budget imbalance BIM of −0.1 GtC yr−1 indicating a near balance between estimated sources and sinks over the last decade. For the year 2019 alone, the growth in EFOS was only about 0.1 % with fossil emissions increasing to 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was 5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminary data for 2020, accounting for the COVID-19-induced changes in emissions, suggest a decrease in EFOS relative to 2019 of about −7 % (median estimate) based on individual estimates from four studies of −6 %, −7 %, −7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. Comparison of estimates from diverse approaches and observations shows (1) no consensus in the mean and trend in land-use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent discrepancy between the different methods for the ocean sink outside the tropics, particularly in the Southern Ocean. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Friedlingstein et al., 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2020 is an update of the previous version, GLODAPv2.2019. The major changes are data from 106 new cruises added, extension of time coverage to 2019, and the inclusion of available (also for historical cruises) discrete fugacity of CO2 (fCO2) values in the merged product files. GLODAPv2.2020 now includes measurements from more than 1.2 million water samples from the global oceans collected on 946 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) 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 updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 106 new cruises with the data from the 840 quality-controlled cruises of the GLODAPv2.2019 data product using crossover analysis. Comparisons to empirical algorithm estimates provided additional context for adjustment decisions; this is new to this version. The adjustments 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−1 in dissolved inorganic carbon, 4 µmol kg−1 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 fCO2, were not subjected to bias comparison or adjustments. The original data and their documentation and DOI codes are available at the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2020/, last access: 20 June 2020). 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/2c8h-sa89 (Olsen et al., 2020). These bias-adjusted product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This living data update documents the GLODAPv2.2020 methods and provides a broad overview of the secondary quality control procedures and results.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use change (ELUC), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), EFF was 9.5±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history, ELUC was 1.5±0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019).

The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface to bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of water samples. This update of GLODAPv2, v2.2019, adds data from 116 cruises to the previous version, extending its coverage in time from 2013 to 2017, while also adding some data from prior years. GLODAPv2.2019 includes measurements from more than 1.1 million water samples from the global oceans collected on 840 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control, especially systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 116 new cruises with the data from the 724 quality-controlled cruises of the GLODAPv2 data product. They correct for errors related to measurement, calibration, and data handling practices, taking into account any known or likely time trends or variations. 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−1 in dissolved inorganic carbon, 4 µmol kg−1 in total alkalinity, 0.01–0.02 in pH, and 5 % in the halogenated transient tracers. The compilation also includes data for several other variables, such as isotopic tracers. These were not subjected to bias comparison or adjustments. The original data, their documentation and DOI codes are available in the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2019/, last access: 17 September 2019). 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/xnme-wr20 (Olsen et al., 2019). The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2.2019 methods and provides a broad overview of the secondary quality control procedures and results.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use and land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was 9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1, SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of 0.5 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in EFF was about 1.6 % and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017, ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 405.0±0.1 ppm averaged over 2017. For 2018, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959–2017, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018, 2016, 2015a, b, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2018.

Abstract Free-ocean CO2 enrichment (FOCE) experiments have been deployed in marine ecosystems to manipulate carbonate system conditions to those predicted in future oceans. We investigated whether the pH/carbonate chemistry of extremely cold polar waters can be manipulated in an ecologically relevant way, to represent conditions under future atmospheric CO2 levels, in an in-situ FOCE experiment in Antarctica. We examined spatial and temporal variation in local ambient carbonate chemistry at hourly intervals at two sites between December and February and compared these with experimental conditions. We successfully maintained a mean pH offset in acidified benthic chambers of −0.38 (±0.07) from ambient for approximately 8 weeks. Local diel and seasonal fluctuations in ambient pH were duplicated in the FOCE system. Large temporal variability in acidified chambers resulted from system stoppages. The mean pH, Ωarag and fCO2 values in the acidified chambers were 7.688 ± 0.079, 0.62 ± 0.13 and 912 ± 150 µatm, respectively. Variation in ambient pH appeared to be mainly driven by salinity and biological production and ranged from 8.019 to 8.192 with significant spatio-temporal variation. This experiment demonstrates the utility of FOCE systems to create conditions expected in future oceans that represent ecologically relevant variation, even under polar conditions.

The Southern Ocean provides a vital service by absorbing about one-sixth of humankind's annual emissions of CO2. This comes with a cost – an increase in ocean acidity that is expected to have negative impacts on ocean ecosystems. The reduced ability of phytoplankton and zooplankton to precipitate carbonate shells is a clearly identified risk. The impact depends on the significance of these organisms in Southern Ocean ecosystems, but there is very little information on their abundance or distribution. To quantify their presence, we used coulometric measurement of particulate inorganic carbonate (PIC) on particles filtered from surface seawater into two size fractions: 50–1000 µm to capture foraminifera (the most important biogenic carbonate-forming zooplankton) and 1–50 µm to capture coccolithophores (the most important biogenic carbonate-forming phytoplankton). Ancillary measurements of biogenic silica (BSi) and particulate organic carbon (POC) provided context, as estimates of the biomass of diatoms (the highest biomass phytoplankton in polar waters) and total microbial biomass, respectively. Results for nine transects from Australia to Antarctica in 2008–2015 showed low levels of PIC compared to Northern Hemisphere polar waters. Coccolithophores slightly exceeded the biomass of diatoms in subantarctic waters, but their abundance decreased more than 30-fold poleward, while diatom abundances increased, so that on a molar basis PIC was only 1 % of BSi in Antarctic waters. This limited importance of coccolithophores in the Southern Ocean is further emphasized in terms of their associated POC, representing less than 1 % of total POC in Antarctic waters and less than 10 % in subantarctic waters. NASA satellite ocean-colour-based PIC estimates were in reasonable agreement with the shipboard results in subantarctic waters but greatly overestimated PIC in Antarctic waters. Contrastingly, the NASA Ocean Biogeochemical Model (NOBM) shows coccolithophores as overly restricted to subtropical and northern subantarctic waters. The cause of the strong southward decrease in PIC abundance in the Southern Ocean is not yet clear. The poleward decrease in pH is small, and while calcite saturation decreases strongly southward, it remains well above saturation ( > 2). Nitrate and phosphate variations would predict a poleward increase. Temperature and competition with diatoms for limiting iron appear likely to be important. While the future trajectory of coccolithophore distributions remains uncertain, their current low abundances suggest small impacts on overall Southern Ocean pelagic ecology.

Biogeochemical change in the water masses of the Southern Ocean, south of Tasmania, was assessed for the 16-year period between 1995 and 2011 using data from four summer repeats of the WOCE–JGOFS–CLIVAR–GO-SHIP (Key et al., 2015; Olsen et al., 2016) SR03 hydrographic section (at ∼ 140° E). Changes in temperature, salinity, oxygen, and nutrients were used to disentangle the effect of solubility, biology, circulation and anthropogenic carbon (CANT) uptake on the variability of dissolved inorganic carbon (DIC) for eight water mass layers defined by neutral surfaces (γn). CANT was estimated using an improved back-calculation method. Warming (∼ 0.0352 ± 0.0170 °C yr−1) of Subtropical Central Water (STCW) and Antarctic Surface Water (AASW) layers decreased their gas solubility, and accordingly DIC concentrations increased less rapidly than expected from equilibration with rising atmospheric CO2 (∼ 0.86 ± 0.16 µmol kg−1 yr−1 versus ∼ 1 ± 0.12 µmol kg−1 yr−1). An increase in apparent oxygen utilisation (AOU) occurred in these layers due to either remineralisation of organic matter or intensification of upwelling. The range of estimates for the increases in CANT were 0.71 ± 0.08 to 0.93 ± 0.08 µmol kg−1 yr−1 for STCW and 0.35 ± 0.14 to 0.65 ± 0.21 µmol kg−1 yr−1 for AASW, with the lower values in each water mass obtained by assigning all the AOU change to remineralisation. DIC increases in the Sub-Antarctic Mode Water (SAMW, 1.10 ± 0.14 µmol kg−1 yr−1) and Antarctic Intermediate Water (AAIW, 0.40 ± 0.15 µmol kg−1 yr−1) layers were similar to the calculated CANT trends. For SAMW, the CANT increase tracked rising atmospheric CO2. As a consequence of the general DIC increase, decreases in total pH (pHT) and aragonite saturation (ΩAr) were found in most water masses, with the upper ocean and the SAMW layer presenting the largest trends for pHT decrease (∼ −0.0031 ± 0.0004 yr−1). DIC increases in deep and bottom layers (∼ 0.24 ± 0.04 µmol kg−1 yr−1) resulted from the advection of old deep waters to resupply increased upwelling, as corroborated by increasing silicate (∼ 0.21 ± 0.07 µmol kg−1 yr−1), which also reached the upper layers near the Antarctic Divergence (∼ 0.36 ± 0.06 µmol kg−1 yr−1) and was accompanied by an increase in salinity. The observed changes in DIC over the 16-year span caused a shoaling (∼ 340 m) of the aragonite saturation depth (ASD, ΩAr = 1) within Upper Circumpolar Deep Water that followed the upwelling path of this layer. From all our results, we conclude a scenario of increased transport of deep waters into the section and enhanced upwelling at high latitudes for the period between 1995 and 2011 linked to strong westerly winds. Although enhanced upwelling lowered the capacity of the AASW layer to uptake atmospheric CO2, it did not limit that of the newly forming SAMW and AAIW, which exhibited CANT storage rates (∼ 0.41 ± 0.20 mol m−2 yr−1) twice that of the upper layers.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2006–2015), EFF was 9.3 ± 0.5 GtC yr−1, ELUC 1.0 ± 0.5 GtC yr−1, GATM 4.5 ± 0.1 GtC yr−1, SOCEAN 2.6 ± 0.5 GtC yr−1, and SLAND 3.1 ± 0.9 GtC yr−1. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr−1 that took place during 2006–2015. Also, for 2015, ELUC was 1.3 ± 0.5 GtC yr−1, GATM was 6.3 ± 0.2 GtC yr−1, SOCEAN was 3.0 ± 0.5 GtC yr−1, and SLAND was 1.9 ± 0.9 GtC yr−1. GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565 ± 55 GtC (2075 ± 205 GtCO2) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2016).

The Surface Ocean CO2 Atlas (SOCAT) is a synthesis of quality-controlled fCO2 (fugacity of carbon dioxide) values for the global surface oceans and coastal seas with regular updates. Version 3 of SOCAT has 14.7 million fCO2 values from 3646 data sets covering the years 1957 to 2014. This latest version has an additional 4.6 million fCO2 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 fCO2 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 fCO2 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.

The ocean has become more acidic over the last 200 years in response increasing atmospheric carbon dioxide (CO2) levels. Documenting how the ocean has changed is critical for assessing how these changes impact marine ecosystems and for the management of marine resources. Here we use present-day ocean carbon observations, from shelf and offshore waters around Australia, combined with neural network mapping of CO2, sea surface temperature, and salinity to estimate the current seasonal and regional distributions of carbonate chemistry (pH and aragonite saturation state). The observed changes in atmospheric CO2 and sea surface temperature (SST) and climatological salinity are then used to reconstruct pH and aragonite saturation state changes over the last 140 years (1870–2013). The comparison with data collected at Integrated Marine Observing System National Reference Station sites located on the shelf around Australia shows that both the mean state and seasonality in the present day are well represented, with the exception of sites such as the Great Barrier Reef. Our reconstruction predicts that since 1870 decrease in aragonite saturation state of 0.48 and of 0.09 in pH has occurred in response to increasing oceanic uptake of atmospheric CO2. Large seasonal variability in pH and aragonite saturation state occur in southwestern Australia driven by ocean dynamics (mixing) and in the Tasman Sea by seasonal warming (in the case of the aragonite saturation state). The seasonal and historical changes in aragonite saturation state and pH have different spatial patterns and suggest that the biological responses to ocean acidification are likely to be non-uniform depending on the relative sensitivity of organisms to shifts in pH and saturation state. This new historical reconstruction provides an important link to biological observations that will help to elucidate the consequences of ocean acidification.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2, and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), EFF was 9.0 ± 0.5 GtC yr−1, ELUC was 0.9 ± 0.5 GtC yr−1, GATM was 4.4 ± 0.1 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 3.0 ± 0.8 GtC yr−1. For the year 2014 alone, EFF grew to 9.8 ± 0.5 GtC yr−1, 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr−1 that took place during 2005–2014. Also, for 2014, ELUC was 1.1 ± 0.5 GtC yr−1, GATM was 3.9 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and SLAND was 4.1 ± 0.9 GtC yr−1. GATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger SLAND for that year. The global atmospheric CO2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in EFF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of EFF and assumed constant ELUC for 2015, cumulative emissions of CO2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO2) for 1870–2015, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2015).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr−1, ELUC 0.9 ± 0.5 GtC yr−1, GATM 4.3 ± 0.1 GtC yr−1, SOCEAN 2.5 ± 0.5 GtC yr−1, and SLAND 2.8 ± 0.8 GtC yr−1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr−1, 2.2% above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and assuming an ELUC of 1.0 ± 0.5 GtC yr−1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr−1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70% from EFF (390 ± 20 GtC) and 30% from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3).

The Surface Ocean CO2 Atlas (SOCAT), an activity of the international marine carbon research community, provides access to synthesis and gridded fCO2 (fugacity of carbon dioxide) products for the surface oceans. Version 2 of SOCAT is an update of the previous release (version 1) with more data (increased from 6.3 million to 10.1 million surface water fCO2 values) and extended data coverage (from 1968–2007 to 1968–2011). The quality control criteria, while identical in both versions, have been applied more strictly in version 2 than in version 1. The SOCAT website (http://www.socat.info/) has links to quality control comments, metadata, individual data set files, and synthesis and gridded data products. Interactive online tools allow visitors to explore the richness of the data. Applications of SOCAT include process studies, quantification of the ocean carbon sink and its spatial, seasonal, year-to-year and longerterm variation, as well as initialisation or validation of ocean carbon models and coupled climate-carbon models. Data coverage Repository-References: Individual data set files and synthesis product: doi:10.1594/PANGAEA.811776 Gridded products: doi:10.3334/CDIAC/OTG.SOCAT_V2_GRID Available at: http://www.socat.info/ Coverage: 79° S to 90° N; 180° W to 180° E Location Name: Global Oceans and Coastal Seas Date/Time Start: 16 November 1968 Date/Time End: 26 December 2011

We investigated the biogeochemistry of low dissolved oxygen high-nitrate (LDOHN) layers forming against the backdrop of several interleaving regional water masses in the eastern Indian Ocean, off northwest Australia adjacent to Ningaloo Reef. These water masses, including the forming Leeuwin Current, have been shown directly to impact the ecological function of Ningaloo Reef and other iconic coastal habitats downstream. Our results indicate that LDOHN layers are formed from multiple subduction events of the Eastern Gyral Current beneath the Leeuwin Current (LC); the LC originates from both the Indonesian Throughflow and tropical Indian Ocean. Density differences of up to 0.025 kg m−3 between the Eastern Gyral Current and the Leeuwin Current produce sharp gradients that can trap high concentrations of particles (measured as low transmission) along the density interfaces. The oxidation of the trapped particulate matter results in local depletion of dissolved oxygen and regeneration of dissolved nitrate (nitrification). We document an associated increase in total dissolved carbon dioxide, which lowers the seawater pH by 0.04 units. Based on isotopic measurements (δ15N and δ18O) of dissolved nitrate, we determine that ~ 40–100% of the nitrate found in LDOHN layers is likely to originate from nitrogen fixation, and that, regionally, the importance of N-fixation in contributing to LDOHN layers is likely to be highest at the surface and offshore.

Ocean acidification is a threat to marine ecosystems globally. In shallow-water systems, however, ocean acidification can be masked by benthic carbon fluxes, depending on community composition, seawater residence time, and the magnitude and balance of net community production (NCP) and calcification (NCC). Here, we examine how six benthic groups from a coral reef environment on Heron Reef (Great Barrier Reef, Australia) contribute to changes in the seawater aragonite saturation state (Ωa). Results of flume studies using intact reef habitats (1.2 m by 0.4 m), showed a hierarchy of responses across groups, depending on CO2 level, time of day and water flow. At low CO2 (350–450 μatm), macroalgae (Chnoospora implexa), turfs and sand elevated Ωa of the flume water by around 0.10 to 1.20 h−1 – normalised to contributions from 1 m2 of benthos to a 1 m deep water column. The rate of Ωa increase in these groups was doubled under acidification (560–700 μatm) and high flow (35 compared to 8 cm s−1). In contrast, branching corals (Acropora aspera) increased Ωa by 0.25 h−1 at ambient CO2 (350–450 μatm) during the day, but reduced Ωa under acidification and high flow. Nighttime changes in Ωa by corals were highly negative (0.6–0.8 h−1) and exacerbated by acidification. Calcifying macroalgae (Halimeda spp.) raised Ωa by day (by around 0.13 h−1), but lowered Ωa by a similar or higher amount at night. Analyses of carbon flux contributions from benthic communities with four different compositions to the reef water carbon chemistry across Heron Reef flat and lagoon indicated that the net lowering of Ωa by coral-dominated areas can to some extent be countered by long water-residence times in neighbouring areas dominated by turfs, macroalgae and carbonate sand.

The Southern Ocean (44–75° S) plays a critical role in the global carbon cycle, yet remains one of the most poorly sampled ocean regions. Different approaches have been used to estimate sea–air CO2 fluxes in this region: synthesis of surface ocean observations, ocean biogeochemical models, and atmospheric and ocean inversions. As part of the RECCAP (REgional Carbon Cycle Assessment and Processes) project, we combine these different approaches to quantify and assess the magnitude and variability in Southern Ocean sea–air CO2 fluxes between 1990–2009. Using all models and inversions (26), the integrated median annual sea–air CO2 flux of −0.42 ± 0.07 Pg C yr−1 for the 44–75° S region, is consistent with the −0.27 ± 0.13 Pg C yr−1 calculated using surface observations. The circumpolar region south of 58° S has a small net annual flux (model and inversion median: −0.04 ± 0.07 Pg C yr−1 and observations: +0.04 ± 0.02 Pg C yr−1), with most of the net annual flux located in the 44 to 58° S circumpolar band (model and inversion median: −0.36 ± 0.09 Pg C yr−1 and observations: −0.35 ± 0.09 Pg C yr−1). Seasonally, in the 44–58° S region, the median of 5 ocean biogeochemical models captures the observed sea–air CO2 flux seasonal cycle, while the median of 11 atmospheric inversions shows little seasonal change in the net flux. South of 58° S, neither atmospheric inversions nor ocean biogeochemical models reproduce the phase and amplitude of the observed seasonal sea–air CO2 flux, particularly in the Austral Winter. Importantly, no individual atmospheric inversion or ocean biogeochemical model is capable of reproducing both the observed annual mean uptake and the observed seasonal cycle. This raises concerns about projecting future changes in Southern Ocean CO2 fluxes. The median interannual variability from atmospheric inversions and ocean biogeochemical models is substantial in the Southern Ocean; up to 25% of the annual mean flux, with 25% of this interannual variability attributed to the region south of 58° S. Resolving long-term trends is difficult due to the large interannual variability and short time frame (1990–2009) of this study; this is particularly evident from the large spread in trends from inversions and ocean biogeochemical models. Nevertheless, in the period 1990–2009 ocean biogeochemical models do show increasing oceanic uptake consistent with the expected increase of −0.05 Pg C yr−1 decade−1. In contrast, atmospheric inversions suggest little change in the strength of the CO2 sink broadly consistent with the results of Le Quéré et al. (2007).

As a response to public demand for a well-documented, quality controlled, publically available, global surface ocean carbon dioxide (CO2) data set, the international marine carbon science community developed the Surface Ocean CO2 Atlas (SOCAT). The first SOCAT product is a collection of 6.3 million quality controlled surface CO2 data from the global oceans and coastal seas, spanning four decades (1968–2007). The SOCAT gridded data presented here is the second data product to come from the SOCAT project. Recognizing that some groups may have trouble working with millions of measurements, the SOCAT gridded product was generated to provide a robust, regularly spaced CO2 fugacity (fCO2) product with minimal spatial and temporal interpolation, which should be easier to work with for many applications. Gridded SOCAT is rich with information that has not been fully explored yet (e.g., regional differences in the seasonal cycles), but also contains biases and limitations that the user needs to recognize and address (e.g., local influences on values in some coastal regions).

A well-documented, publicly available, global data set of surface ocean carbon dioxide (CO2) parameters has been called for by international groups for nearly two decades. The Surface Ocean CO2 Atlas (SOCAT) project was initiated by the international marine carbon science community in 2007 with the aim of providing a comprehensive, publicly available, regularly updated, global data set of marine surface CO2, which had been subject to quality control (QC). Many additional CO2 data, not yet made public via the Carbon Dioxide Information Analysis Center (CDIAC), were retrieved from data originators, public websites and other data centres. All data were put in a uniform format following a strict protocol. Quality control was carried out according to clearly defined criteria. Regional specialists performed the quality control, using state-of-the-art web-based tools, specially developed for accomplishing this global team effort. SOCAT version 1.5 was made public in September 2011 and holds 6.3 million quality controlled surface CO2 data points from the global oceans and coastal seas, spanning four decades (1968–2007). Three types of data products are available: individual cruise files, a merged complete data set and gridded products. With the rapid expansion of marine CO2 data collection and the importance of quantifying net global oceanic CO2 uptake and its changes, sustained data synthesis and data access are priorities.

The roles of iron and light in controlling biomass and primary productivity are clearly established in the Southern Ocean. However, their influence on net community production (NCP) and carbon export remains to be quantified. To improve our understanding of NCP and carbon export production in the Subantarctic Zone (SAZ) and the northern reaches of the Polar Frontal Zone (PFZ), we conducted continuous onboard determinations of NCP as part of the Sub-Antarctic Sensitivity to Environmental Change (SAZ-Sense) study, which occurred in January–February 2007. Biological O2 supersaturation was derived from measuring O2/Ar ratios by equilibrator inlet mass spectrometry. Based on these continuous measurements, NCP during the austral summer 2007 in the Australian SAZ was approximately 43 mmol O2 m−2 d−1. NCP showed significant spatial variability, with larger values near the Subtropical front, and a general southward decrease. For shallower mixed layers (

The CARINA project is aimed at gathering and providing secondary quality control checks on carbon and carbon-relevant hydrographic and geochemical data from cruises all across the Atlantic, Arctic and Southern Ocean. In total the project gathered 188 cruises that were not previously available to the public. Of these 188 cruises, 37 are part of the Southern Ocean. Parameters from the Southern Ocean cruises, including total carbon dioxide (TCO2), total alkalinity, oxygen, nitrate, phosphate and silicate, were examined for cruise-to-cruise consistency. pH and chlorofluorocarbons (CFCs) are also part of the data base, but are not discussed here. This paper focuses on the quality control of the Southern Ocean data from the Pacific sector which consisted of 29 cruises of which 17 were included in a previous synthesis called GLODAP, 11 were new cruises from the CARINA dataset, and one cruise was included in GLODAP but was updated with new data and therefore also included in CARINA. The Pacific sector quality control procedures included crossover analysis between stations and inversion analysis of all crossover data. The GLODAP data were included into the analysis as reference cruises but without applying the GLODAP recommended adjustments so the corrections could be independently verified. The outcome of this effort is an internally consistent, high-quality carbon data set for all cruises, including the reference cruises.

Carbon and carbon-relevant hydrographic and hydrochemical ancillary data from previously not publicly available cruises were retrieved and recently merged to a new data base, CARINA (CARbon IN the Atlantic). The initial North Atlantic project, an international effort for ocean carbon synthesis, was extended to include the Arctic Mediterranean Seas (Arctic Ocean and Nordic Seas) and all three sectors of the Southern Ocean. Of a total of 188 cruises, 37 cruises are part of the Southern Ocean. The present work focuses on data collected in the Indian sector (20° S–70° S; 30° E–150° E). The Southern Indian Ocean dataset covers the period 1992–2004 and includes seasonal repeated observations. Parameters including salinity, dissolved inorganic carbon (TCO2), total alkalinity (TA), oxygen, nitrate, phosphate and silicate were examined for cruise-to-cruise and overall consistency. In addition, data from an existing, quality controlled data base (GLODAP) were introduced in the CARINA analysis to improve data coverage in the Southern Ocean. A global inversion was performed to synthesize the information deduced from objective comparisons of deep measurements (>1500 m) at nearby stations (generally 2 and other properties of the carbon system that can now be used to investigate the natural variability or stability of ocean chemistry and the accumulation of anthropogenic carbon. This data product also offers an important new synthesis of seasonal to decadal observations to validate ocean biogeochemical models in a region where available historical data were very sparse.

We compiled a large data-set from 22 cruises spanning from 1991 to 2003, of the partial pressure of CO2 (pCO2) in surface waters over the continental shelf (CS) and adjacent open ocean (43° to 46° S; 145° to 150° E), south of Tasmania. Climatological seasonal cycles of pCO2 in the CS, the subtropical zone (STZ) and the subAntarctic zone (SAZ) are described and used to determine monthly pCO2 anomalies. These are used in combination with monthly anomalies of sea surface temperature (SST) to investigate inter-annual variations of SST and pCO2. Monthly anomalies of SST (as intense as 2°C) are apparent in the CS, STZ and SAZ, and are indicative of strong inter-annual variability that seems to be related to large-scale coupled atmosphere-ocean oscillations. Anomalies of pCO2 normalized to a constant temperature are negatively related to SST anomalies. A reduced winter-time vertical input of dissolved inorganic carbon (DIC) during phases of positive SST anomalies, related to a poleward shift of westerly winds, and a concomitant local decrease in wind stress is the likely cause of the negative relationship between pCO2 and SST anomalies. The observed pattern is an increase of the sink for atmospheric CO2 associated with positive SST anomalies, although strongly modulated by inter-annual variability of wind speed. Assuming that phases of positive SST anomalies are indicative of the future evolution of regional ocean biogeochemistry under global warming, we show using a purely observational based approach that some provinces of the Southern Ocean could provide a potential negative feedback on increasing atmospheric CO2.

Coral reefs are highly sensitive to ocean acidification due to rising atmospheric CO 2 concentrations. We present 10 years of data (2009-2019) on the long-term trends and sources of variation in the carbon chemistry from two fixed stations in the Australian Great Barrier Reef. Data from the subtropical mid-shelf GBRWIS comprised 3-h instrument records, and those from the tropical coastal NRSYON were monthly seawater samples. Both stations recorded significant variation in seawater CO 2 fugacity (fCO 2 ), attributable to seasonal, daytime, temperature and salinity fluctuations. Superimposed over this variation, fCO 2 progressively increased by > 2.0 ± 0.3 µatm year -1 at both stations. Seawater temperature and salinity also increased throughout the decade, whereas seawater pH and the saturation state of aragonite declined. The decadal upward fCO 2 trend remained significant in temperature- and salinity-normalised data. Indeed, annual fCO 2 minima are now higher than estimated fCO 2 maxima in the early 1960s, with mean fCO 2 now ~ 28% higher than 60 years ago. Our data indicate that carbonate dissolution from the seafloor is currently unable to buffer the Great Barrier Reef against ocean acidification. This is of great concern for the thousands of coral reefs and other diverse marine ecosystems located in this vast continental shelf system.

Ocean warming, ocean acidification and overfishing are major threats to the structure and function of marine ecosystems. Driven by increasing anthropogenic emissions of CO2, ocean warming is leading to global redistribution of marine biota and altered ecosystem dynamics, while ocean acidification threatens the ability of calcifying marine organisms to form skeletons due to decline in saturation state of carbonate Ω and pH. In Tasmania, the interaction between overfishing of sea urchin predators and rapid ocean warming has caused a phase-shift from productive kelp beds to overgrazed sea urchin barren grounds, however potential impacts of ocean acidification on this system have not been considered despite this threat for marine ecosystems globally. Here we use automated loggers and point measures of pH, spanning kelp beds and barren grounds, to reveal that kelp beds have the capacity to locally ameliorate effects of ocean acidification, via photosynthetic drawdown of CO2, compared to unvegetated barren grounds. Based on meta-analysis of anticipated declines in physiological performance of grazing urchins to decreasing pH and assumptions of nil adaptation, future projection of OA across kelp-barrens transition zones reveals that kelp beds could act as important pH refugia, with urchins potentially becoming increasingly challenged at distances >40 m from kelp beds. Using spatially explicit simulation of physicochemical feedbacks between grazing urchins and their kelp prey, we show a stable mosaicked expression of kelp patches to emerge on barren grounds. Depending on the adaptative capacity of sea urchins, future declines in pH appear poised to further alter phase-shift dynamics for reef communities; thus, assessing change in spatial-patterning of reef-scapes may indicate cascading ecological impacts of ocean acidification., Competing Interests: The authors have declared that no competing interests exist.

To improve estimates of the long-term response of the marine carbon system to climate change a better understanding of the seasonal and interannual variability is needed. We use high-frequency multi-year data at three locations identified as climate change hotspots: two sites located close to South Pacific boundary currents and one in the Subantarctic Zone (SAZ). We investigate and identify the main drivers involved in the seasonal an interannual (2012-2016) variability of the carbon system. The seasonal variability at boundary current sites is temporally different and highly controlled by sea surface temperature. Advection processes also play a significant role on the monthly changes of the carbon system at the western boundary current site. The interannual variability at these sites most likely responds to long-term variability in oceanic circulation ultimately related to climatic indices such as the El Niño Southern Oscillation, the Pacific Decadal Oscillation and the Southern Annular Mode (SAM). In the SAZ, advection and entrainment processes drive most of the seasonality, augmented by the action of biological processes in spring. Given the relevance of advection and entrainment processes at SAZ, the interannual variability is most probably modulated by changes in the regional winds linked to the variability of the SAM.

Free-ocean CO 2 enrichment (FOCE) experiments have been deployed in marine ecosystems to manipulate carbonate system conditions to those predicted in future oceans. We investigated whether the pH/carbonate chemistry of extremely cold polar waters can be manipulated in an ecologically relevant way, to represent conditions under future atmospheric CO 2 levels, in an in-situ FOCE experiment in Antarctica. We examined spatial and temporal variation in local ambient carbonate chemistry at hourly intervals at two sites between December and February and compared these with experimental conditions. We successfully maintained a mean pH offset in acidified benthic chambers of -0.38 (±0.07) from ambient for approximately 8 weeks. Local diel and seasonal fluctuations in ambient pH were duplicated in the FOCE system. Large temporal variability in acidified chambers resulted from system stoppages. The mean pH, Ω arag and fCO 2 values in the acidified chambers were 7.688 ± 0.079, 0.62 ± 0.13 and 912 ± 150 µatm, respectively. Variation in ambient pH appeared to be mainly driven by salinity and biological production and ranged from 8.019 to 8.192 with significant spatio-temporal variation. This experiment demonstrates the utility of FOCE systems to create conditions expected in future oceans that represent ecologically relevant variation, even under polar conditions.

Coral reefs are one of the most vulnerable ecosystems to ocean acidification. While our understanding of the potential impacts of ocean acidification on coral reef ecosystems is growing, gaps remain that limit our ability to translate scientific knowledge into management action. To guide solution-based research, we review the current knowledge of ocean acidification impacts on coral reefs alongside management needs and priorities. We use the world's largest continuous reef system, Australia's Great Barrier Reef (GBR), as a case study. We integrate scientific knowledge gained from a variety of approaches (e.g., laboratory studies, field observations, and ecosystem modelling) and scales (e.g., cell, organism, ecosystem) that underpin a systems-level understanding of how ocean acidification is likely to impact the GBR and associated goods and services. We then discuss local and regional management options that may be effective to help mitigate the effects of ocean acidification on the GBR, with likely application to other coral reef systems. We develop a research framework for linking solution-based ocean acidification research to practical management options. The framework assists in identifying effective and cost-efficient options for supporting ecosystem resilience. The framework enables on-the-ground OA management to be the focus, while not losing sight of CO2 mitigation as the ultimate solution., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The Great Barrier Reef (GBR) is founded on reef-building corals. Corals build their exoskeleton with aragonite, but ocean acidification is lowering the aragonite saturation state of seawater (Ωa). The downscaling of ocean acidification projections from global to GBR scales requires the set of regional drivers controlling Ωa to be resolved. Here we use a regional coupled circulation-biogeochemical model and observations to estimate the Ωa experienced by the 3,581 reefs of the GBR, and to apportion the contributions of the hydrological cycle, regional hydrodynamics and metabolism on Ωa variability. We find more detail, and a greater range (1.43), than previously compiled coarse maps of Ωa of the region (0.4), or in observations (1.0). Most of the variability in Ωa is due to processes upstream of the reef in question. As a result, future decline in Ωa is likely to be steeper on the GBR than currently projected by the IPCC assessment report.

Several studies have suggested that the carbon sink in the Southern Ocean-the ocean's strongest region for the uptake of anthropogenic CO2 -has weakened in recent decades. We demonstrated, on the basis of multidecadal analyses of surface ocean CO2 observations, that this weakening trend stopped around 2002, and by 2012, the Southern Ocean had regained its expected strength based on the growth of atmospheric CO2. All three Southern Ocean sectors have contributed to this reinvigoration of the carbon sink, yet differences in the processes between sectors exist, related to a tendency toward a zonally more asymmetric atmospheric circulation. The large decadal variations in the Southern Ocean carbon sink suggest a rather dynamic ocean carbon cycle that varies more in time than previously recognized., (Copyright © 2015, American Association for the Advancement of Science.)

Sustained observations allow for the tracking of change in oceanography and ecosystems, however, these are rare, particularly for the Southern Hemisphere. To address this in part, the Australian Integrated Marine Observing System (IMOS) implemented a network of nine National Reference Stations (NRS). The network builds on one long-term location, where monthly water sampling has been sustained since the 1940s and two others that commenced in the 1950s. In-situ continuously moored sensors and an enhanced monthly water sampling regime now collect more than 50 data streams. Building on sampling for temperature, salinity and nutrients, the network now observes dissolved oxygen, carbon, turbidity, currents, chlorophyll a and both phytoplankton and zooplankton. Additional parameters for studies of ocean acidification and bio-optics are collected at a sub-set of sites and all data is made freely and publically available. Our preliminary results demonstrate increased utility to observe extreme events, such as marine heat waves and coastal flooding; rare events, such as plankton blooms; and have, for the first time, allowed for consistent continental scale sampling and analysis of coastal zooplankton and phytoplankton communities. Independent water sampling allows for cross validation of the deployed sensors for quality control of data that now continuously tracks daily, seasonal and annual variation. The NRS will provide multi-decadal time series, against which more spatially replicated short-term studies can be referenced, models and remote sensing products validated, and improvements made to our understanding of how large-scale, long-term change and variability in the global ocean are affecting Australia's coastal seas and ecosystems. The NRS network provides an example of how a continental scaled observing systems can be developed to collect observations that integrate across physics, chemistry and biology.

Coral reef calcification is predicted to decline as a result of ocean acidification and other anthropogenic stressors. The majority of studies predicting declines based on in situ relationships between environmental parameters and net community calcification rate have been location-specific, preventing accurate predictions for coral reefs globally. In this study, net community calcification and production were measured on a coral reef flat at One Tree Island, Great Barrier Reef, using Lagrangian flow respirometry and slack water methods. Net community calcification, daytime net photosynthesis and nighttime respiration were higher under the flow respirometry method, likely due to increased water flow relative to the slack water method. The two methods also varied in the degrees to which they were influenced by potential measurement uncertainties. The difference in the results from these two commonly used methods implies that some of the location-specific differences in coral reef community metabolism may be due to differences in measurement methods.

Stenothermal polar benthic marine invertebrates are highly sensitive to environmental perturbations but little is known about potential synergistic effects of concurrent ocean warming and acidification on development of their embryos and larvae. We examined the effects of these stressors on development to the calcifying larval stage in the Antarctic sea urchin Sterechinus neumayeri in embryos reared in present and future (2100+) ocean conditions from fertilization. Embryos were reared in 2 temperature (ambient: -1.0 °C, + 2 °C : 1.0 °C) and 3 pH (ambient: pH 8.0, -0.2-0.4 pH units: 7.8,7.6) levels. Principle coordinates analysis on five larval metrics showed a significant effect of temperature and pH on the pattern of growth. Within each temperature, larvae were separated by pH treatment, a pattern primarily influenced by larval arm and body length. Growth was accelerated by temperature with a 20-28% increase in postoral (PO) length at +2 °C across all pH levels. Growth was strongly depressed by reduced pH with a 8-19% decrease in PO length at pH 7.6-7.8 at both temperatures. The boost in growth caused by warming resulted in larvae that were larger than would be observed if acidification was examined in the absence of warming. However, there was no significant interaction between these stressors. The increase in left-right asymmetry and altered body allometry indicated that decreased pH disrupted developmental patterning and acted as a teratogen (agent causing developmental malformation). Decreased developmental success with just a 2 °C warming indicates that development in S. neumayeri is particularly sensitive to increased temperature. Increased temperature also altered larval allometry. Altered body shape impairs swimming and feeding in echinoplutei. In the absence of adaptation, it appears that the larval phase may be a bottleneck for survivorship of S. neumayeri in a changing ocean in a location where poleward migration to escape inhospitable conditions is not possible., (© 2013 Blackwell Publishing Ltd.)

Ocean acidification, via an anthropogenic increase in seawater carbon dioxide (CO2 ), is potentially a major threat to coral reefs and other marine ecosystems. However, our understanding of how natural short-term diurnal CO2 variability in coral reefs influences longer term anthropogenic ocean acidification remains unclear. Here, we combine observed natural carbonate chemistry variability with future carbonate chemistry predictions for a coral reef flat in the Great Barrier Reef based on the RCP8.5 CO2 emissions scenario. Rather than observing a linear increase in reef flat partial pressure of CO2 (pCO2 ) in concert with rising atmospheric concentrations, the inclusion of in situ diurnal variability results in a highly nonlinear threefold amplification of the pCO2 signal by the end of the century. This significant nonlinear amplification of diurnal pCO2 variability occurs as a result of combining natural diurnal biological CO2 metabolism with long-term decreases in seawater buffer capacity, which occurs via increasing anthropogenic CO2 absorption by the ocean. Under the same benthic community composition, the amplification in the variability in pCO2 is likely to lead to exposure to mean maximum daily pCO2 levels of ca. 2100 μatm, with corrosive conditions with respect to aragonite by end-century at our study site. Minimum pCO2 levels will become lower relative to the mean offshore value (ca. threefold increase in the difference between offshore and minimum reef flat pCO2 ) by end-century, leading to a further increase in the pCO2 range that organisms are exposed to. The biological consequences of short-term exposure to these extreme CO2 conditions, coupled with elevated long-term mean CO2 conditions are currently unknown and future laboratory experiments will need to incorporate natural variability to test this. The amplification of pCO2 that we describe here is not unique to our study location, but will occur in all shallow coastal environments where high biological productivity drives large natural variability in carbonate chemistry., (© 2013 Blackwell Publishing Ltd.)

Space competition between corals and seaweeds is an important ecological process underlying coral-reef dynamics. Processes promoting seaweed growth and survival, such as herbivore overfishing and eutrophication, can lead to local reef degradation. Here, we present the case that increasing concentrations of atmospheric CO(2) may be an additional process driving a shift from corals to seaweeds on reefs. Coral (Acropora intermedia) mortality in contact with a common coral-reef seaweed (Lobophora papenfussii) increased two- to threefold between background CO(2) (400 ppm) and highest level projected for late 21st century (1140 ppm). The strong interaction between CO(2) and seaweeds on coral mortality was most likely attributable to a chemical competitive mechanism, as control corals with algal mimics showed no mortality. Our results suggest that coral (Acropora) reefs may become increasingly susceptible to seaweed proliferation under ocean acidification, and processes regulating algal abundance (e.g. herbivory) will play an increasingly important role in maintaining coral abundance., (© 2010 Blackwell Publishing Ltd/CNRS.)

The oxygen (O(2)) concentration in the surface ocean is influenced by biological and physical processes. With concurrent measurements of argon (Ar), which has similar solubility properties as oxygen, we can remove the physical contribution to O(2) supersaturation and determine the biological oxygen supersaturation. Biological O(2) supersaturation in the surface ocean reflects the net metabolic balance between photosynthesis and respiration, i.e., the net community productivity (NCP). We present a new method for continuous shipboard measurements of O(2)/Ar by equilibrator inlet mass spectrometry (EIMS). From these measurements and an appropriate gas exchange parametrization, NCP can be estimated at high spatial and temporal resolution. In the EIMS configuration, seawater from the ship's continuous intake flows through a cartridge enclosing a gas-permeable microporous membrane contactor. Gases in the headspace of the cartridge equilibrate with dissolved gases in the flowing seawater. A fused-silica capillary continuously samples headspace gases, and the O(2)/Ar ratio is measured by mass spectrometry. The ion current measurements on the mass spectrometer reflect the partial pressures of dissolved gases in the water flowing through the equilibrator. Calibration of the O(2)/Ar ion current ratio (32/40) is performed automatically every 2 h by sampling ambient air through a second capillary. A conceptual model demonstrates that the ratio of gases reaching the mass spectrometer is dependent on several parameters, such as the differences in molecular diffusivities and solubilities of the gases. Laboratory experiments and field observations performed by EIMS are discussed. We also present preliminary evidence that other gas measurements, such as N(2)/Ar and pCO(2) measurements, may potentially be performed with EIMS. Finally, we compare the characteristics of the EIMS with the previously described membrane inlet mass spectrometry (MIMS) approach.

Biogeochemical rate processes in the Southern Ocean have an important impact on the global environment. Here, we summarize an extensive set of published and new data that establishes the pattern of gross primary production and net community production over large areas of the Southern Ocean. We compare these rates with model estimates of dissolved iron that is added to surface waters by aerosols. This comparison shows that net community production, which is comparable to export production, is proportional to modeled input of soluble iron in aerosols. Our results strengthen the evidence that the addition of aerosol iron fertilizes export production in the Southern Ocean. The data also show that aerosol iron input particularly enhances gross primary production over the large area of the Southern Ocean downwind of dry continental areas.

Using inorganic carbon measurements from an international survey effort in the 1990s and a tracer-based separation technique, we estimate a global oceanic anthropogenic carbon dioxide (CO2) sink for the period from 1800 to 1994 of 118 +/- 19 petagrams of carbon. The oceanic sink accounts for approximately 48% of the total fossil-fuel and cement-manufacturing emissions, implying that the terrestrial biosphere was a net source of CO2 to the atmosphere of about 39 +/- 28 petagrams of carbon for this period. The current fraction of total anthropogenic CO2 emissions stored in the ocean appears to be about one-third of the long-term potential.

The delta(13)C value of the dissolved inorganic carbon in the surface waters of the Pacific Ocean has decreased by about 0.4 per mil between 1970 and 1990. This decrease has resulted from the uptake of atmospheric CO(2) derived from fossil fuel combustion and deforestation. The net amounts of CO(2) taken up by the oceans and released from the biosphere between 1970 and 1990 have been determined from the changes in three measured values: the concentration of atmospheric CO(2), the delta(13)C of atmospheric CO(2) and the delta(13)C value of dissolved inorganic carbon in the ocean. The calculated average net oceanic CO(2) uptake is 2.1 gigatons of carbon per year. This amount implies that the ocean is the dominant net sink for anthropogenically produced CO(2) and that there has been no significant net CO(2) released from the biosphere during the last 20 years.