Despite the importance of the Southern Ocean carbon sink, its response to future atmospheric CO2 perturbations and warming remains highly uncertain. In this study, we use six state-of-the-art Earth system models to assess the response of Southern Ocean air-sea CO2 fluxes (FCO2) to a rapid atmospheric forcing increase and subsequent negative emissions in an idealized carbon dioxide removal reversibility experiment. We find that during positive emissions, the region north of the Polar Front only takes up atmospheric CO2 for 30-50 years before reaching equilibrium; surface stratification and reduction of CO2 solubility with warming diminishes ocean CO2 uptake in this region. In contrast, south of the Polar Front, the upper ocean continues to take up CO2 until the end of positive emissions at 140 years. Sea-ice loss and the accumulation of anthropogenic dissolved inorganic carbon in the upper ocean reduce the upwelling-driven seasonal CO2 outgassing, leading to a stronger Antarctic CO2 sink. CO2 removal triggers a CO2 uptake reduction that slowly converts the Southern Ocean into a CO2 source which persists for at least 50 years post-mitigation. Furthermore, we find that model sensitivity to atmospheric perturbation is closely linked to seasonal FCO2 dynamics. Specifically, models with a thermally dominated pCO2 seasonal cycle exhibit nearly twice the sensitivity to atmospheric perturbations compared to non-thermal models. Our findings further emphasize the necessity of accurate model representation of the seasonal CO2 dynamics for appropriately simulating the future Southern Ocean carbon sink.

The ocean is a major carbon sink and takes up 25-30% of the anthropogenically emitted CO2. A state-of-the-art method to quantify this sink are global ocean biogeochemistry models (GOBMs) but their simulated CO2 uptake differs between models and is systematically lower than estimates based on statistical methods using surface ocean pCO2 and interior ocean measurements. Here, we provide an in-depth evaluation of ocean carbon sink estimates from 1980 to 2018 from a GOBM ensemble. As sources of inter-model differences and ensemble-mean biases our study identifies the (i) model set-up, such as the length of the spin-up, the starting date of the simulation, and carbon fluxes from rivers and into sediments, (ii) the ocean circulation, such as Atlantic Meridional Overturning Circulation and Southern Ocean mode and intermediate water formation, and (iii) the oceanic buffer capacity. Our analysis suggests that the late starting date and biases in the ocean circulation cause a too low anthropogenic CO2 uptake across the GOBM ensemble. Surface ocean biogeochemistry biases might also cause simulated anthropogenic fluxes to be too low but the current set-up prevents a robust assessment. For simulations of the ocean carbon sink, we recommend in the short-term to (1) start simulations in 1765, when atmospheric CO2 started to increase, (2) conduct a sufficiently long spin-up such that the GOBMs reach steady-state, and (3) provide key metrics for circulation, biogeochemistry, and the land-ocean interface. In the long-term, we recommend improving the representation of these metrics in the GOBMs.