The representation of alkalinity and the carbonate pump from CMIP5 to CMIP6 ESMs and implications for the ocean carbon cycle

Abstract. Ocean alkalinity is critical to the uptake of atmospheric carbon in surface waters and provides buffering capacity towards associated acidification. However, unlike dissolved inorganic carbon (DIC), alkalinity is not directly impacted by anthropogenic carbon emissions. Within the context of projections of future ocean carbon uptake and potential ecosystem impacts, especially through Coupled Model Intercomparison Projects (CMIPs), the representation of alkalinity and the main driver of its distribution in the ocean interior, the calcium carbonate cycle, have often been overlooked. Here we track the changes from CMIP5 to CMIP6 with respect to the Earth system model (ESM) representation of alkalinity and the carbonate pump which depletes the surface ocean in alkalinity through biological production of calcium carbonate, and releases it at depth through export and dissolution. We report a significant improvement in the representation of alkalinity in CMIP6 ESMs relative to those in CMIP5. This improvement can be explained in part by an increase in calcium carbonate (CaCO₃) production for some ESMs, which redistributes alkalinity at the surface and strengthens its vertical gradient in the water column. We were able to constrain a PIC export estimate of 51–70 Tmol yr^-1 at 100 m for the ESMs to match the observed vertical gradient of alkalinity. Biases in the vertical profile of DIC have also significantly decreased, especially with the enhancement of the carbonate pump, but the representation of the saturation horizons has slightly worsened in contrast. Reviewing the representation of the CaCO₃ cycle across CMIP5/6, we find a substantial range of parameterizations. While all biogeochemical models currently represent pelagic calcification, they do so implicitly, and they do not represent benthic calcification. In addition, most models simulate marine calcite but not aragonite. In CMIP6 certain model groups have increased the complexity of simulated CaCO₃ production, sinking, dissolution and sedimentation. However, this is insufficient to explain the overall improvement in the alkalinity representation, which is therefore likely a result of improved marine biogeochemistry model tuning or ad hoc parameterizations. We find differences in the way ocean alkalinity is initialized that lead to offsets of up to 1 % in the global alkalinity inventory of certain models. These initialization biases should be addressed in future CMIPs by adopting accurate unit conversions. Although modelers aim to balance the global alkalinity budget in ESMs in order to limit drift in ocean carbon uptake under preindustrial conditions, varying assumptions in the closure of the budget have the potential to influence projections of future carbon uptake. For instance, in many models, carbonate production, dissolution and burial are independent of the seawater saturation state, and when considered, the range of sensitivities is substantial. As such, the future impact of ocean acidification on the carbonate pump, and in turn ocean carbon uptake, is potentially underestimated in current ESMs and insufficiently constrained.