On the solution of the carbonate chemistry system in ocean biogeochemistry models

We present a simplified method for solving the local equilibrium carbonate chemistry in numerical ocean biogeochemistry models. Compared to the methods typically used, the scheme is fast, efficient and compact. The accuracy of the solution is dictated by the number of species retained in the expression for alkalinity and there is almost no computational penalty for retaining minor contributions. We demonstrate that this scheme accurately reproduces the results of the commonly used method in the context of a three-dimensional global ocean carbon cycle model. Using this model we also show that neglecting the regional variations in surface dissolved inorganic phosphorus and silicic acid concentrations can lead to significant systematic bias in regional estimates of air–sea carbon fluxes using such models.