Abstract
<jats:p>Abstract. Constraining the Southern Ocean carbon sink is vital for understanding global carbon cycling and climate. Monthly 1°x1° gridded products of surface seawater partial pressure of carbon dioxide (pCO2) are routinely used to detect change in Southern Ocean carbon cycling. Despite Southern Ocean sampling being sparse and skewed to summer months, the uncertainty arising from incomplete sampling at the monthly 1°x1° scale is unquantified to date, and it remains unclear how this pCO2 sampling error affects estimates of monthly gridded pCO2 products and derived air-sea CO2 fluxes. Here, we quantify the Southern Ocean pCO2 sampling error using synthetic observations from an eddy-permitting ocean-sea ice-biogeochemistry configuration of the Energy Exascale Earth System Model. Using instantaneous snapshots of the model state at the times and locations of real-world ship and float observations, we reconstruct gridded monthly 1°x1° pCO2 datasets over 2014–2023 and compare them with the simulated model pCO2 output, isolating the error due to incomplete observational coverage at the monthly 1°x1° scale. We find that, on average, pCO2 derived from synthetic observations is 13–14 μatm lower than the modeled monthly mean, with individual grid cells exhibiting uncertainties up to an order of magnitude larger, primarily caused by undersampling during high-wind events in the Southern Ocean. Our results imply that the rate of oceanic CO2 uptake derived from available pCO2 observations is biased high. Since the pCO2 uncertainty due to incomplete observational coverage exceeds other sources of uncertainty in gridded pCO2 data products, e.g., measurement error, it represents a substantial, previously overlooked contribution that should be propagated through gap-filled pCO2 data products and estimates of the Southern Ocean and global ocean carbon sink.</jats:p>