Geoscientific Model Development (Mar 2024)

Modelling water isotopologues (<sup>1</sup>H<sup>2</sup>H<sup>16</sup>O, <sup>1</sup>H<sub>2</sub><sup>17</sup>O) in the coupled numerical climate model iLOVECLIM (version 1.1.5)

  • T. Extier,
  • T. Caley,
  • D. M. Roche,
  • D. M. Roche

DOI
https://doi.org/10.5194/gmd-17-2117-2024
Journal volume & issue
Vol. 17
pp. 2117 – 2139

Abstract

Read online

Stable water isotopes are used to infer changes in the hydrological cycle for different climate periods and various climatic archives. Following previous developments of δ18O in the coupled climate model of intermediate complexity, iLOVECLIM, we present here the implementation of the 1H2H16O and 1H217O water isotopes in the different components of this model and calculate the associated secondary markers deuterium excess (d-excess) and oxygen-17 excess (17O-excess) in the atmosphere and ocean. So far, the latter has only been modelled by the atmospheric model LMDZ4. Results of a 5000-year equilibrium simulation under preindustrial conditions are analysed and compared to observations and several isotope-enabled models for the atmosphere and ocean components. In the atmospheric component, the model correctly reproduces the first-order global distribution of the δ2H and d-excess as observed in the data (R=0.56 for δ2H and 0.36 for d-excess), even if local differences are observed. The model–data correlation is within the range of other water-isotope-enabled general circulation models. The main isotopic effects and the latitudinal gradient are properly modelled, similarly to previous water-isotope-enabled general circulation model simulations, despite a simplified atmospheric component in iLOVECLIM. One exception is observed in Antarctica where the model does not correctly estimate the water isotope composition, a consequence of the non-conservative behaviour of the advection scheme at a very low moisture content. The modelled 17O-excess presents a too-important dispersion of the values in comparison to the observations and is not correctly reproduced in the model, mainly because of the complex processes involved in the 17O-excess isotopic value. For the ocean, the model simulates an adequate isotopic ratio in comparison to the observations, except for local areas such as the surface of the Arabian Sea, a part of the Arctic and the western equatorial Indian Ocean. Data–model evaluation also presents a good match for the δ2H over the entire water column in the Atlantic Ocean, reflecting the influence of the different water masses.