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

Abstract 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.

[1]  Janet G. Hering,et al.  Principles and Applications of Aquatic Chemistry , 1993 .

[2]  J.-F. Minster,et al.  Tracers in the Sea , 1982 .

[3]  Edward A. Boyle,et al.  Decoupling of iron and phosphate in the global ocean , 2005 .

[4]  William H. Press,et al.  Numerical Recipes in Fortran 77: The Art of Scientific Computing 2nd Editionn - Volume 1 of Fortran Numerical Recipes , 1992 .

[5]  Reiner Schlitzer,et al.  Carbon export fluxes in the Southern Ocean: results from inverse modeling and comparison with satellite based estimates , 2002 .

[6]  E. Boyle,et al.  Is AOU a good measure of respiration in the oceans? , 2004, Geophysical Research Letters.

[7]  D. Wolf-Gladrow,et al.  CO2 in Seawater: Equilibrium, Kinetics, Isotopes , 2001 .

[8]  R. Wanninkhof.  Relationship between wind speed and gas exchange over the ocean , 1992 .

[9]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[10]  Thomas Kaminski,et al.  Recipes for adjoint code construction , 1998, TOMS.

[11]  R. Weiss.  Carbon dioxide in water and seawater: the solubility of a non-ideal gas , 1974 .

[12]  M. Follows,et al.  The solubility pump of carbon in the subtropical gyre of the North Atlantic , 1996 .

[13]  Wallace S. Broecker,et al.  Atmospheric response to deep-sea injections of fossil-fuel carbon dioxide , 1979 .

[14]  C. Sweeney,et al.  Global sea-air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects , 2002 .

[15]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[16]  John Marshall,et al.  Evaluating carbon sequestration efficiency in an ocean circulation model by adjoint sensitivity analysis , 2004 .

[17]  William H. Press,et al.  Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing , 1992 .

[18]  F. Millero.  Thermodynamics of the carbon dioxide system in the oceans , 1995 .

[19]  S. Gorshkov,et al.  World ocean atlas , 1976 .

[20]  Andrew G. Dickson,et al.  Handbook of methods for the analysis of the various parameters of the carbon dioxide system in sea water. Version 2 , 1994 .

[21]  A. Dickson.  An exact definition of total alkalinity and a procedure for the estimation of alkalinity and total inorganic carbon from titration data , 1981 .

[22]  L. Perelman,et al.  Hydrostatic, quasi‐hydrostatic, and nonhydrostatic ocean modeling , 1997 .

[23]  P. Gent,et al.  Isopycnal mixing in ocean circulation models , 1990 .

[24]  Stephanie Dutkiewicz,et al.  Interactions of the iron and phosphorus cycles: A three‐dimensional model study , 2005 .