LDEO pCO2-Residual Method

The ocean reduces human impacts on global climate by absorbing and sequestering CO2 from the atmosphere. To quantify global, time-resolved air-sea CO2 fluxes, surface ocean pCO2 is needed. A common approach for estimating full-coverage pCO2 is to train a machine learning algorithm on sparse in situ pCO2 data and associated physical and biogeochemical observations. Though these associated variables have understood relationships to pCO2, it is often unclear how they drive pCO2 outputs. Here, we make two advances that enhance connections between physical understanding and reconstructed pCO2. First, we apply pre-processing to the pCO2 data to remove the direct effect of temperature. This enhances the biogeochemical/physical component of pCO2 in the target variable and reduces the complexity that the machine learning must disentangle. Second, we demonstrate that the resulting algorithm has physically understandable connections between input data and the output biogeochemical/physical component of pCO2. The final pCO2 reconstruction agrees modestly better with independent data than most other approaches. Uncertainties in the reconstructed pCO2 and impacts on the estimated CO2 fluxes are quantified. Uncertainty in piston velocity drives substantial flux uncertainties in some regions, but does not increase globally integrated estimates of uncertainty in CO2 fluxes from observation-based products. Our reconstructed CO2 fluxes show larger interannual variability than smoother neural network approaches, but a lesser trend since 2005.