The Cryosphere (Aug 2024)
Post-depositional modification on seasonal-to-interannual timescales alters the deuterium-excess signals in summer snow layers in Greenland
Abstract
We document the isotopic evolution of near-surface snow at the East Greenland Ice Core Project (EastGRIP) ice core site in northeast Greenland using a time-resolved array of 1 m deep isotope (δ18O, δD) profiles. The snow profiles were taken from May–August during the 2017–2019 summer seasons. An age–depth model was developed and applied to each profile, mitigating the impacts of stratigraphic noise on isotope signals. Significant changes in deuterium excess (d) are observed in surface snow and near-surface snow as the snow ages. Decreases in d of up to 5 ‰ occur during summer seasons after deposition during two of the three summer seasons observed. The d always experiences a 3 ‰–5 ‰ increase after aging 1 year in the snow due to a broadening of the autumn d maximum. Models of idealized scenarios coupled with prior work indicate that the summertime post-depositional changes in d (Δd) can be explained by a combination of surface sublimation, forced ventilation of the near-surface snow down to 20–30 cm, and isotope-gradient-driven diffusion throughout the column. The interannual Δd is also partly explained with isotope-gradient-driven diffusion, but other mechanisms are at work that leave a bias in the d record. Thus, d does not just carry information about source-region conditions and transport history as is commonly assumed, but also integrates local conditions into summer snow layers as the snow ages through metamorphic processes. Finally, we observe a dramatic increase in the seasonal isotope-to-temperature sensitivity, which can be explained solely by isotope-gradient-driven diffusion. Our results are dependent on the site characteristics (e.g., wind, temperature, accumulation rate, snow properties) but indicate that more process-based research is necessary to understand water isotopes as climate proxies. Recommendations for monitoring and physical modeling are given, with special attention to the d parameter.