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; Watson et al., 2020). 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., 2019). 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., 2019) 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).