Stable water isotopes in the MITgcm

We present the first results of the implementation of stable water isotopes in the Massachusetts Institute of Technology general circulation model (MITgcm). The model is forced with the isotopic content of precipitation and water vapor from an atmospheric general circulation model (NCAR IsoCAM), while the fractionation during evaporation is treated explicitly in the MITgcm. Results of the equilibrium simulation under pre-industrial conditions are compared to observational data and measurements of plankton tow records (the oxygen isotopic composition of planktic foraminiferal calcite). The broad patterns and magnitude of the stable water isotopes in annual mean seawater are well captured in the model, both at the sea surface as well as in the deep ocean. However, the surface water in the Arctic Ocean is not depleted enough, due to the absence of highly depleted precipitation and snowfall. A model-data mismatch is also recognizable in the isotopic composition of the seawater-salinity relationship in midlatitudes that is mainly caused by the coarse grid resolution. Deep-ocean characteristics of the vertical water mass distribution in the Atlantic Ocean closely resemble observational data. The reconstructed delta O-18(c) at the sea surface shows a good agreement with measurements. However, the model-data fit is weaker when individual species are considered and deviations are most likely attributable to the habitat depth of the foraminifera. Overall, the newly developed stable water isotope package opens wide prospects for long-term simulations in a paleoclimatic context.