Description: Published empirical algorithms for oceanic total alkalinity (TA) and dissolved inorganic carbon (DIC) are used with monthly sea surface salinity (SSS) and temperature (SST) derived from satellite (SMOS, Aquarius, SST CCI) and interpolated in situ (CORA) measurements and climatological (WOA) ancillary data to produce monthly maps of TA and DIC at one degree spatial resolution. Earth system model TA and DIC (HADGEM2-ES) are also included. Results are compared with in situ (GLODAPv2) TA and DIC and results analysed in five regions (global, Greater Caribbean, Amazon plume, Amazon plume with in situ SSS 〈 35 and Bay of Bengal). Results are presented in three versions, denoted by 'X' in the lists below: using all available data (X = ''); excluding data with bathymetry 〈 500m (X = 'Depth500'); excluding data with both bathymetry 〈 500m and distance from nearest coast 〈 300 km (X = 'Depth500Dist300'). Datasets S1 to S5 are .csv lists of matchups in each region - date and location, in situ TA and DIC measurements and estimated uncertainties, all input datasets, estimates of TA and DIC from all outputs, and the best available output estimates of TA and DIC for each matchup. S1_GlobalAlgorithmMatchupsX.csv S2_GreaterCaribbeanAlgorithmMatchupsX.csv S3_AmazonPlumeAlgorithmMatchupsX.csv S4_AmazonPlumeLowSAlgorithmMatchupsX.csv S5_BayOfBengalAlgorithmMatchupsX.csv Datasets S6 to S10 are .csv statistical analyses of the performance of each combination of algorithm and input data - carbonate system variable, algorithm, input datasets used, (MAD, RMSD using all available data, output score, RMSD estimated from output score, output and in situ mean and standard deviation, correlation coefficient), all items in brackets presented both unweighted and weighted, number of matchups, number of potential matchups, matchup coverage, RMSD after subtraction of linear regression, percentage reduction in RMSD due to subtraction of linear regression and weighted score divided by number of matchups). S6_GlobalAlgorithmScoresX.csv S7_GreaterCaribbeanAlgorithmScoresX.csv S8_AmazonPlumeAlgorithmScoresX.csv S9_AmazonPlumeLowSAlgorithmScoresX.csv S10_BayOfBengalAlgorithmScoresX.csv Datasets S11 to S15 are zipped netCDF files containing error analyses of all outputs in each region, including the squared error of each output at each matchup, the weight of each squared error (1/squared uncertainty), weight * squared error, number of matchups available to each output, number of matchups available to each combination of two outputs, (score of each output in a given comparison of two outputs, overall output score and RMSD estimated from output score), all items in the last brackets presented both unweighted and weighted. S11_GlobalSquaredErrorsX.nc S12_GreaterCaribbeanSquaredErrorsX.nc S13_AmazonPlumeSquaredErrorsX.nc S14_AmazonPlumeLowSSquaredErrorsX.nc S15_BayOfBengalSquaredErrorsX.nc Datasets S16 to S20 are zipped netCDF files containing global maps of the mean and standard deviation of each of: in situ data; output data; output data - in situ data and number of matchups. Regional files show the same maps, but only including data within the region. S16_GlobalmapsX.nc S17_GreaterCaribbeanmapsX.nc S18_AmazonPlumemapsX.nc S19_AmazonPlumeLowSmapsX.nc S20_BayOfBengalmapsX.nc Datasets S21 and S22 are .csv files containing the effect on estimated RMSD of excluding various combinations of algorithms and/or inputs for TA and DIC in each region. For a given variable and region, the first line shows the algorithm, input data sources, estimated RMSD and bias of the output with lowest estimated RMSD. Subsequent lines show the effect of excluding combinations of algorithms and/or inputs, ordered first by the number of algorithms/inputs excluded (fewest first), then by effect on lowest estimated RMSD. So the first line(s) consist of the effects of excluding the best algorithm and each of the input sources to that algorithm, most important first. Each line consists of the item excluded, ratio of resulting estimated RMSD to original estimated RMSD, resulting bias and number of items excluded. Some exclusions are equivalent, for instance exclusion of WOA nitrate (the only nitrate source) is equivalent to excluding all algorithms using nitrate. Dataset S21 contains a comprehensive list of all possible exclusions, and so is rather hard to read and interpret. To mitigate this, Dataset S22 contains only those exclusion sets with effect greater than 1% and at least 0.1% greater than any subset of its exclusions. S21_importancesX.csv S22_importances2X.csv Dataset S23 is a .csv file containing like-for-like comparisons of RMSD between TA and DIC in each region. Bear in mind that the RMSD shown here is not the same as the estimated RMSD (RMSDe) shown elsewhere. S23_TA_DICcomparisonX.csv

Description: This dataset includes an inventory of all published marine turtle stable isotope studies to help inform ecology and conservation research. We conducted an extensive literature search in Scopus, Web of Science, and Google Scholar in June 2018. All peer-reviewed primary research papers were included in the analysis excluding fossil isotope studies. The inventory provides information on each publication. Version number: v1.3 Submission date stamp: 2019-03-21 Contact email: Julia.haywood@exeter.ac.uk

Description: This submission includes the reference data required to perform a complete verification of the FluxEngine v3.0 install. All data are in netCDF-3 format. Note that this dataset is greater than 100 MB in size. FluxEngine is an open source software toolkit for calculating in situ, regional or global gas fluxes between the atmosphere and ocean. It can be used with model, in situ or satellite Earth observation data. A full description of the toolkit is provided in Shutler et al. (2016) and the FluxEngine software can be freely downloaded from GitHub: https://github.com/oceanflux-ghg/FluxEngine Input data are for the year 2010 with the exception of the Takahashi climatology, which is for the year 2000. The following input data are included: * input_data/air_pressure - atmospheric pressure data from the European Centre for Medium-Range Weather Forecasts (ECMWF, www.ecmwf.int) * input_data/ice - fraction ice coverage from the Ocean and Sea Ice Satellite Application Facility (OSI SAF, www.osi-saf.org) * input_data/rain_gpcp - total precipitation data from the National Oceanic and Atmospheric Administration (NOAA) Global Precipitation Climatology Project (v2.2) * input_data/sig_wv_ht - significant wave height data from the GlobWave project (www.globwave.org) * input_data/sigma0 - radar backscatter from the GlobWave project (www.globwave.org) * input_data/windu10 - wind speed at 10 m above sea level from the GlobWave project (www.globwave.org) * input_data/SMOS - ocean salinity data from Centre Abal de Traitement des Données Soil Moisture and Ocean Salinity (CATDS SMOS v01, www.catds.fr) * input_data/SOCATv4 - The Surface Ocean CO₂ Atlas (SOCAT) (Bakker et al., 2016) version 4 reanalysed to a CO₂ climatology using (Goddijn-Murphy et al., 2015) * input_data/SST - skin sea surface temperature (SST) data from the ATSR Reprocessing for Climate (ARC) project. * input_data/sstfnd_Reynolds - sub skin sea surface temperature (Optimally Interpolated SST, OISST project: Reynolds et al., 2007, Banzon et al., 2016) * input_data/takahashi09 - Takahashi CO₂ climatology dataset (Takahashi et al., 2009) Where appropriate input data have been re-gridded from their original resolution into a monthly 1 × 1 degree resolution. The verification process involves comparing a newly generated output (e.g. generated by a local install of FluxEngine) to a reference dataset for each verification scenario. These reference datasets are included in the 'reference_data' directory, and use the following naming convention: * socatv4, takahashi09_pco2 - uses partial pressure of CO₂ (pCO₂) data from SOCATv4 or Takahashi2009 climatologies, respectively * sst, no_gradients - does/does not use sea surface temperature gradients, respectively * salinity - applies a correction for skin layer salinity * N00 - uses (Nightingale et al., 2000) parameterisation for calculating gas transfer velocity * K0, K1, K2, K3 - uses the generic formulation for calculating gas transfer velocity (using only the 0th, 1st, 2nd and 3rd order components, respectively) * takahashi09_all_inputs - The Takahashi et al. (2009) dataset is used for all inputs in order to reproduce the results in Shutler et al. (2016). * FEv1, FEv2, FEv3 - Reference data was generated using FluxEngine version 1, 2 or 3, respectively Configuration files for each validation set are included in the 'configs' directory. These files are well commented and fully specify the input data to be used, data pre-processing, gas transfer velocity parameterisation and the structure of the gas flux calculation to be performed. Acknowledgements: The Surface Ocean CO₂ Atlas (SOCAT) is an international effort, endorsed by the International Ocean Carbon Coordination Project (IOCCP), the Surface Ocean Lower Atmosphere Study (SOLAS) and the Integrated Marine Biosphere Research (IMBeR) program, to deliver a uniformly quality-controlled surface ocean CO₂ database. The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to SOCAT.

Description: The Surface Ocean CO₂ Atlas (SOCAT) version 2019 dataset (Bakker et al., 2016) is a quality-controlled dataset containing 25.7 million surface ocean gaseous CO₂ measurements collated from thousands of individual submissions. These gaseous CO₂ measurements are typically collected at many different depths (of the order of several metres below the surface) using many different systems, and the sampling depth varies dependent upon the sampling platform and/or setup. Different platforms (e.g. ships of opportunity, research vessels) and systems will collect water samples at different depths, and the sampling depth can even vary dependent upon sea state. Therefore, the collated SOCAT dataset contains high quality data, but these data are all valid for different and inconsistent depths. Therefore the SOCAT provided individual gaseous CO₂ measurements and gridded data are sub-optimal for calculating global or regional atmosphere-ocean gas exchange (and the resultant net CO₂ sinks) and sub-optimal for verifying gas fluxes from (or assimilation into) numerical models. Accurate calculations of CO₂ flux between the atmosphere and oceans require CO₂ concentrations at the top and bottom of the mass boundary layer, the ~100 μm deep layer that forms the interface between the ocean and the atmosphere (Woolf et al., 2016). Ignoring vertical temperature gradients across this very small layer can result in significant biases in the concentration differences and the resulting gas fluxes (e.g. ~5 to 29% underestimate in global net CO₂ sink values, Woolf et al., 2016). It is currently impossible to measure the CO₂ concentrations either side of this very thin layer, but it is possible to calculate the concentrations either side of this layer using the SOCAT data, satellite observations and knowledge of the carbonate system. Therefore to enable the SOCAT data to be optimal for an accurate atmosphere-ocean gas flux calculation, a reanalysis methodology was developed to enable the calculation of the fugacity of CO₂ (fCO₂) for the bottom of the mass boundary layer (termed sub-skin value). The theoretical basis and justification for this is described in detail within Woolf et al., (2016) and the re-analysis methodology is described in detail in (Goddijn-Murphy et al., 2015). The re-analysis calculation exploits paired in situ temperature and fCO₂ measurements in the SOCAT dataset, and uses an Earth observation dataset to provide a depth-consistent (sub-skin) temperature field to which all fugacity data are reanalysed. The outputs provide paired fCO₂ (and partial pressure of CO₂) and temperature data that correspond to a consistent sub-skin layer temperature. These can then be used to accurately calculate concentration differences and atmosphere-ocean CO₂ gas fluxes. This data submission contains a reanalysis of the fugacity of CO₂ (fCO₂) from the SOCAT version 2019 dataset to a consistent sub-skin temperature field. The reanalysis was performed using a tool that is distributed within the FluxEngine open source software toolkit (https://github.com/oceanflux-ghg/FluxEngine) (Shutler et al., 2016; Holding et al., in-review). All data processing and driver scripts are available from the FluxEngine ancillary tools repository https://github.com/oceanflux-ghg/FluxEngineAncillaryTools. The NOAA Optimum Interpolation Sea Surface Temperature (OISST) dataset (Reynolds et al., 2007) was used to provide the climate quality and depth consistent temperature data. The original OISST data were first resampled to provide monthly mean values on a 1º by 1º degree grid. These data were then used as the temperature input for the reanalysis. The resulting reanalysed data are provided as a tab-separated value file (individual data points) and as netCDF-5 file (gridded monthly means). These are the same file formats as provided by SOCAT and analogous to the SOCAT single data point and gridded data. Each row in the tab-separated value file corresponds to a row in the original SOCAT version 2019 dataset. The original SOCAT version 2019 data are included in full, with four additional columns containing the reanalysed data: * T_reynolds - The temperature (in degrees C) taken from the consistent OISST temperature field for the corresponding time and location. * fCO2_reanalysed - The fugacity of CO₂ (in μatm) reanalysed to the consistent surface temperature indicated by T_reynolds. * pCO2_SST - The partial pressure of CO₂ (in μatm) corresponding to the in situ (measured) temperature. * pCO2_reanalysed - The partial pressure of CO₂ (in μatm) reanalysed to the consistent surface temperature indicated by T_reynolds. The netCDF gridded version of the reanalysed dataset contains monthly mean data, binned into a 1º by 1º grid and uses the same units, missing value indicators and time and space resolution as the original SOCAT gridded product to maximise compatibility. The gridding is performed using the SOCAT gridding methodology (Sabine et al. 2013). The implementation of the gridding has been verified by performing the gridding on the original (non-reanalysed) SOCAT data and all results were identical to 8 decimal places. The result of gridding the original SOCAT data are included within these netCDF data, along with additional variables containing the equivalent results for the reanalysed SOCAT data. Statistical sample mean, minimum, maximum, standard deviation and count data for each grid cell are included, with unweighted and cruise-weighted versions (following the convention used by SOCAT). Full meta data are included within the file. Notes: 1. Due to the temporal range of the OISST dataset the reanalysed values are only available from 1981 onwards. Pre-1981 rows contain "NaN" (not-a-number) in the reanalysis columns. 2. The download for this submission is provided as a single .zip file (1.1 GB, uncompressed: 10.7GB) containing two files: SOCATv2019_reanalysed_subskin.tsv (containing every data point, ungridded) and SOCATv2019_reanalysed_subskin.nc (the gridded monthly mean data). How to cite these data: Please cite this PANGAEA submission, the theory (Woolf et al., 2016), the reanalysis methodology (Goddijn-Murphy et al., 2015), the FluxEngine toolbox which was used to perform the reanalysis (Shutler et al., 2016, Holding et al. in review) and the original SOCAT dataset (Bakker et al., 2016) and/or gridded equivalent (Sabine et al., 2013). Acknowledgements: The Surface Ocean CO₂ Atlas (SOCAT) is an international effort, endorsed by the International Ocean Carbon Coordination Project (IOCCP), the Surface Ocean Lower Atmosphere Study (SOLAS) and the Integrated Marine Biosphere Research (IMBeR) program, to deliver a uniformly quality-controlled surface ocean CO₂ database. The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to SOCAT. These data were provided by two Integrated Carbon Observing System (ICOS) European Union (EU) readiness projects, Ringo (grant no. 730944) and BONUS Integral (grant no. 03FO773A).

Description: Ocean surface winds, currents, and waves play a crucial role in exchanges of momentum, energy, heat, freshwater, gases, and other tracers between the ocean, atmosphere, and ice. Despite surface waves being strongly coupled to the upper ocean circulation and the overlying atmosphere, efforts to improve ocean, atmospheric, and wave observations and models have evolved somewhat independently. From an observational point of view, community efforts to bridge this gap have led to proposals for satellite Doppler oceanography mission concepts, which could provide unprecedented measurements of absolute surface velocity and directional wave spectrum at global scales. This paper reviews the present state of observations of surface winds, currents, and waves, and it outlines observational gaps that limit our current understanding of coupled processes that happen at the air-sea-ice interface. A significant challenge for the coming decade of wind, current, and wave observations will come in combining and interpreting measurements from (a) wave-buoys and high-frequency radars in coastal regions, (b) surface drifters and wave-enabled drifters in the open-ocean, marginal ice zones, and wave-current interaction “hot-spots,” and (c) simultaneous measurements of absolute surface currents, ocean surface wind vector, and directional wave spectrum from Doppler satellite sensors.

Description: Shelf seas play an important role in the global carbon cycle, absorbing atmospheric carbon dioxide (CO2) and exporting carbon (C) to the open ocean and sediments. The magnitude of these processes is poorly constrained, because observations are typically interpolated over multiple years. Here, we used 298500 observations of CO2 fugacity (fCO2) from a single year (2015), to estimate the net influx of atmospheric CO2 as 26.2 ± 4.7 Tg C yr−1 over the open NW European shelf. CO2 influx from the atmosphere was dominated by influx during winter as a consequence of high winds, despite a smaller, thermally-driven, air-sea fCO2 gradient compared to the larger, biologically-driven summer gradient. In order to understand this climate regulation service, we constructed a carbon-budget supplemented by data from the literature, where the NW European shelf is treated as a box with carbon entering and leaving the box. This budget showed that net C-burial was a small sink of 1.3 ± 3.1 Tg C yr−1, while CO2 efflux from estuaries to the atmosphere, removed the majority of river C-inputs. In contrast, the input from the Baltic Sea likely contributes to net export via the continental shelf pump and advection (34.4 ± 6.0 Tg C yr−1).

Description: The flow (flux) of climate-critical gases, such as carbon dioxide (CO2), between the ocean and the atmosphere is a fundamental component of our climate and an important driver of the biogeochemical systems within the oceans. Therefore, the accurate calculation of these air–sea gas fluxes is critical if we are to monitor the oceans and assess the impact that these gases are having on Earth's climate and ecosystems. FluxEngine is an open-source software toolbox that allows users to easily perform calculations of air–sea gas fluxes from model, in situ, and Earth observation data. The original development and verification of the toolbox was described in a previous publication. The toolbox has now been considerably updated to allow for its use as a Python library, to enable simplified installation, to ensure verification of its installation, to enable the handling of multiple sparingly soluble gases, and to enable the greatly expanded functionality for supporting in situ dataset analyses. This new functionality for supporting in situ analyses includes user-defined grids, time periods and projections, the ability to reanalyse in situ CO2 data to a common temperature dataset, and the ability to easily calculate gas fluxes using in situ data from drifting buoys, fixed moorings, and research cruises. Here we describe these new capabilities and demonstrate their application through illustrative case studies. The first case study demonstrates the workflow for accurately calculating CO2 fluxes using in situ data from four research cruises from the Surface Ocean CO2 ATlas (SOCAT) database. The second case study calculates air–sea CO2 fluxes using in situ data from a fixed monitoring station in the Baltic Sea. The third case study focuses on nitrous oxide (N2O) and, through a user-defined gas transfer parameterisation, identifies that biological surfactants in the North Atlantic could suppress individual N2O sea–air gas fluxes by up to 13 %. The fourth and final case study illustrates how a dissipation-based gas transfer parameterisation can be implemented and used. The updated version of the toolbox (version 3) and all documentation is now freely available.

Description: The Sea surface KInematics Multiscale monitoring (SKIM) satellite mission is designed to explore ocean surface current and waves. This includes tropical currents, notably the unknown patterns of divergence and their impact on the ocean heat budget near the Equator, monitoring of the emerging Arctic up to 82.5$(\circ)$N. SKIM will also make unprecedented direct measurements of strong currents, from boundary currents to the Antarctic circumpolar current, and their interaction with ocean waves with expected impacts on air-sea fluxes and extreme waves. For the first time, SKIM will directly measure the ocean surface current vector from space. The main instrument on SKIM is a Ka-band conically scanning, multi-beam Doppler radar altimeter/wave scatterometer that includes a state-of-the-art nadir beam comparable to the Poseidon-4 instrument on Sentinel 6. The well proven Doppler pulse-pair technique will give a surface drift velocity representative of the top two meters of the ocean, after subtracting a large wave-induced contribution. Horizontal velocity components will be obtained with an accuracy better than 7 cm/s for horizontal wavelengths larger than 80~km and time resolutions larger than 15 days, with a mean revisit time of 4 days for of 99\% of the global oceans. This will provide unique and innovative measurements that will further our understanding of the transports in the upper ocean layer, permanently distributing heat, carbon, plankton, and plastics. SKIM will also benefit from co-located measurements of water vapor, rain rate, sea ice concentration, and wind vectors provided by the European operational satellite MetOp-SG(B), allowing many joint analyses. SKIM is one of the two candidate satellite missions under development for ESA Earth Explorer 9. The other candidate is the Far infrared Radiation Understanding and Monitoring (FORUM). The final selection will be announced by September 2019, for a launch in the coming decade.

Description: Shelf seas play an important role in the global carbon cycle, absorbing atmospheric carbon dioxide (CO2) and exporting carbon (C) to the open ocean and sediments. The magnitude of these processes is poorly constrained, because observations are typically interpolated over multiple years. Here, we used 298500 observations of CO2 fugacity (fCO2) from a single year (2015), to estimate the net influx of atmospheric CO2 as 26.2 ± 4.7 Tg C yr−1 over the open NW European shelf. CO2 influx from the atmosphere was dominated by influx during winter as a consequence of high winds, despite a smaller, thermallydriven, air-sea fCO 2 gradient compared to the larger, biologically-driven summer gradient. In order to understand this climate regulation service, we constructed a carbon-budget supplemented by data from the literature, where the NW European shelf is treated as a box with carbon entering and leaving the box. This budget showed that net C-burial was a small sink of 1.3 ± 3.1 Tg C yr−1, while CO2 efflux from estuaries to the atmosphere, removed the majority of river C-inputs. In contrast, the input from the Baltic Sea likely contributes to net export via the continental shelf pump and advection (34.4 ± 6.0 Tg C yr−1).

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Boas, A. B. V., Ardhuin, F., Ayet, A., Bourassa, M. A., Brandt, P., Chapron, B., Cornuelle, B. D., Farrar, J. T., Fewings, M. R., Fox-Kemper, B., Gille, S. T., Gommenginger, C., Heimbach, P., Hell, M. C., Li, Q., Mazloff, M. R., Merrifield, S. T., Mouche, A., Rio, M. H., Rodriguez, E., Shutler, J. D., Subramanian, A. C., Terrill, E. J., Tsamados, M., Ubelmann, C., & van Sebille, E. Integrated observations of global surface winds, currents, and waves: Requirements and challenges for the next decade. Frontiers in Marine Science, 6, (2019): 425, doi:10.3389/fmars.2019.00425.

Description: Ocean surface winds, currents, and waves play a crucial role in exchanges of momentum, energy, heat, freshwater, gases, and other tracers between the ocean, atmosphere, and ice. Despite surface waves being strongly coupled to the upper ocean circulation and the overlying atmosphere, efforts to improve ocean, atmospheric, and wave observations and models have evolved somewhat independently. From an observational point of view, community efforts to bridge this gap have led to proposals for satellite Doppler oceanography mission concepts, which could provide unprecedented measurements of absolute surface velocity and directional wave spectrum at global scales. This paper reviews the present state of observations of surface winds, currents, and waves, and it outlines observational gaps that limit our current understanding of coupled processes that happen at the air-sea-ice interface. A significant challenge for the coming decade of wind, current, and wave observations will come in combining and interpreting measurements from (a) wave-buoys and high-frequency radars in coastal regions, (b) surface drifters and wave-enabled drifters in the open-ocean, marginal ice zones, and wave-current interaction “hot-spots,” and (c) simultaneous measurements of absolute surface currents, ocean surface wind vector, and directional wave spectrum from Doppler satellite sensors.

Description: AV was funded by NASA Earth and Space Science Fellowship award number 80NSSC17K0326. MB was funded by NOAA (FundRef number 100007298) through the NGI (grant number 18-NGI3-42). SG was funded by NASA grants NNX16AH67G, NNX14A078G, NNX17AH53G, and 80NSSC19K0059. MT acknowledges support from the Natural Environment Research Council (grant number NE/R000654/1). MT, MR, JS, and EvS were partially funded by the SKIM Mission Science Study (SKIM-SciSoc) project ESA RFP 3-15456/18/NL/CT/gp. AA was supported by DGA grant No D0456JE075 and the French Brittany Regional Council. MF was supported by NASA Ocean Vector Winds Science Team Grant 80NSSC18K1611 and Jet Propulsion Laboratory/CalTech subcontract 1531731. FA, BC, and AM were supported by ESA under the Sea State CCI project, with additional support from CNES and ANR grants for ISblue (ANR-17-EURE-0015) and LabexMER (ANR-10-LABX-19). MZ was funded by NASA (grant number NNX16AH67G).