Abstract. The Southern Ocean (44° S–75° S) plays a critical role in the global carbon cycle, yet remains one of the most poorly sampled ocean regions. Different approaches have been used to estimate sea-air CO2 fluxes in this region: synthesis of surface ocean observations, ocean biogeochemical models, and atmospheric and ocean inversions. As part of the RECCAP (REgional Carbon Cycle Assessment and Processes) project, we combine these different approaches to quantify and assess the magnitude and variability in Southern Ocean sea-air CO2 fluxes between 1990–2009. Using all models and inversions (26), the integrated median annual sea-air CO2 flux of −0.42 ± 0.07 Pg C yr−1 for the 44° S–75° S region is consistent with the −0.27 ± 0.13 Pg C yr−1 calculated using surface observations. The circumpolar region south of 58° S has a small net annual flux (model and inversion median: −0.04 ± 0.07 Pg C yr−1 and observations: +0.04 ± 0.02 Pg C yr−1), with most of the net annual flux located in the 44° S to 58° S circumpolar band (model and inversion median: −0.36 ± 0.09 Pg C yr−1 and observations: −0.35 ± 0.09 Pg C yr−1). Seasonally, in the 44° S–58° S region, the median of 5 ocean biogeochemical models captures the observed sea-air CO2 flux seasonal cycle, while the median of 11 atmospheric inversions shows little seasonal change in the net flux. South of 58° S, neither atmospheric inversions nor ocean biogeochemical models reproduce the phase and amplitude of the observed seasonal sea-air CO2 flux, particularly in the Austral Winter. Importantly, no individual atmospheric inversion or ocean biogeochemical model is capable of reproducing both the observed annual mean uptake and the observed seasonal cycle. This raises concerns about projecting future changes in Southern Ocean CO2 fluxes. The median interannual variability from atmospheric inversions and ocean biogeochemical models is substantial in the Southern Ocean; up to 25% of the annual mean flux with 25% of this inter-annual variability attributed to the region south of 58° S. Trends in the net CO2 flux from the inversions and models are not statistically different from the expected increase of –0.05 Pg C yr−1 decade−1 due to increasing atmospheric CO2 concentrations. However, resolving long term trends is difficult due to the large interannual variability and short time frame (1990–2009) of this study.;
