Search

The Southern Ocean forms an important component of the Earth system as a major sink of CO2 and heat. Recent studies based on the Coupled Model Intercomparison Project version 5 (CMIP5) Earth system models (ESMs) show that CMIP5 models disagree on the phasing of the seasonal cycle of the CO2 flux (FCO2/ and compare poorly with available observation products for the Southern Ocean. Because the seasonal cycle is the dominant mode of CO2 variability in the Southern Ocean, its simulation is a rigorous test for models and their long-term projections. Here we examine the competing roles of temperature and dissolved inorganic carbon (DIC) as drivers of the seasonal cycle of pCO2 in the Southern Ocean to explain the mechanistic basis for the seasonal biases in CMIP5 models. We find that despite significant differences in the spatial characteristics of the mean annual fluxes, the intra-model homogeneity in the seasonal cycle of FCO2 is greater than observational products. FCO2 biases in CMIP5 models can be grouped into two main categories, i.e., group-SST and group-DIC. Group-SST models show an exaggeration of the seasonal rates of change of sea surface temperature (SST) in autumn and spring during the cooling and warming peaks. These higher-than-observed rates of change of SST tip the control of the seasonal cycle of pCO2 and FCO2 towards SST and result in a divergence between the observed and modeled seasonal cycles, particularly in the Sub-Antarctic Zone. While almost all analyzed models (9 out of 10) show these SST-driven biases, 3 out of 10 (namely NorESM1-ME, HadGEM-ES and MPI-ESM, collectively the group-DIC models) compensate for the solubility bias because of their overly exaggerated primary production, such that biologically driven DIC changes mainly regulate the seasonal cycle of FCO2. [ABSTRACT FROM AUTHOR]

Resolving and understanding the drivers of variability of CO2 in the Southern Ocean and its potential climate feedback is one of the major scientific challenges of the ocean-climate community. Here we use a regional approach on empirical estimates of pCO2 to understand the role that seasonal variability has in long-term CO2 changes in the Southern Ocean. Machine learning has become the preferred empirical modelling tool to interpolate time- and locationrestricted ship measurements of pCO2. In this study we use an ensemble of three machine-learning products: support vector regression (SVR) and random forest regression (RFR) from Gregor et al. (2017), and the self-organising-map feedforward neural network (SOM-FFN) method from Landschützer et al. (2016). The interpolated estimates of ΔpCO2 are separated into nine regions in the Southern Ocean defined by basin (Indian, Pacific, and Atlantic) and biomes (as defined by Fay and McKinley, 2014a). The regional approach shows that, while there is good agreement in the overall trend of the products, there are periods and regions where the confidence in estimated ΔpCO2 is low due to disagreement between the products. The regional breakdown of the data highlighted the seasonal decoupling of the modes for summer and winter interannual variability. Winter interannual variability had a longer mode of variability compared to summer, which varied on a 4-6-year timescale. We separate the analysis of the ΔpCO2 and its drivers into summer and winter. We find that understanding the variability of ΔpCO2 and its drivers on shorter timescales is critical to resolving the long-term variability of ΔpCO2. Results show that ΔpCO2 is rarely driven by thermodynamics during winter, but rather by mixing and stratification due to the stronger correlation of ΔpCO2 variability with mixed layer depth. Summer pCO2 variability is consistent with chlorophyll a variability, where higher concentrations of chlorophyll a correspond with lower pCO2 concentrations. In regions of low chlorophyll a concentrations, wind stress and sea surface temperature emerged as stronger drivers of ΔpCO2. In summary we propose that sub-decadal variability is explained by summer drivers, while winter variability contributes to the long-term changes associated with the SAM. This approach is a useful framework to assess the drivers of ΔpCO2 but would greatly benefit from improved estimates of ΔpCO2 and a longer time series. [ABSTRACT FROM AUTHOR]

The Southern Ocean accounts for 40% of oceanic CO2 uptake, but the estimates are bound by large uncertainties due to a paucity in observations. Gap-filling empirical methods have been used to good effect to approximate pCO2 from satellite observable variables in other parts of the ocean, but many of these methods are not in agreement in the Southern Ocean. In this study we propose two additional methods that perform well in the Southern Ocean: support vector regression (SVR) and random forest regression (RFR). The methods are used to estimate ΔpCO2 in the Southern Ocean based on SOCAT v3, achieving similar trends to the SOM-FFN method by Landschützer et al. (2014). Results show that the SOM-FFN and RFR approaches have RMSEs of similar magnitude (14.84 and 16.45 µatm, where 1 atm=101 325 Pa) where the SVR method has a larger RMSE (24.40 µatm). However, the larger errors for SVR and RFR are, in part, due to an increase in coastal observations from SOCAT v2 to v3, where the SOM-FFN method used v2 data. The success of both SOM-FFN and RFR depends on the ability to adapt to different modes of variability. The SOM-FFN achieves this by having independent regression models for each cluster, while this flexibility is intrinsic to the RFR method. Analyses of the estimates shows that the SVR and RFR's respective sensitivity and robustness to outliers define the outcome significantly. Further analyses on the methods were performed by using a synthetic dataset to assess the following: which method (RFR or SVR) has the best performance? What is the effect of using time, latitude and longitude as proxy variables on ΔpCO2? What is the impact of the sampling bias in the SOCAT v3 dataset on the estimates? We find that while RFR is indeed better than SVR, the ensemble of the two methods outperforms either one, due to complementary strengths and weaknesses of the methods. Results also show that for the RFR and SVR implementations, it is better to include coordinates as proxy variables as RMSE scores are lowered and the phasing of the seasonal cycle is more accurate. Lastly, we show that there is only a weak bias due to undersampling. The synthetic data provide a useful framework to test methods in regions of sparse data coverage and show potential as a useful tool to evaluate methods in future studies. [ABSTRACT FROM AUTHOR]