Exploring the role of snow metamorphism on the isotopic composition of the surface snow at EastGRIP

Abstract. Stable water isotopes from polar ice cores are invaluable high-resolution climate proxy records. Recent studies have aimed to improve knowledge of how the climate signal is stored in the water isotope record by addressing the influence of post-depositional processes on the surface snow isotopic composition. In this study, the relationship between changes in surface snow microstructure after precipitation/deposition events and water isotopes is explored using measurements of snow specific surface area (SSA). Continuous daily SSA measurements from the East Greenland Ice Core Project site (EastGRIP) situated in the accumulation zone of the Greenland Ice Sheet during the summer seasons of 2017, 2018 and 2019 are used to develop an empirical decay model to describe events of rapid decrease in SSA, driven predominantly by vapour diffusion in the pore space and atmospheric vapour exchange. The SSA decay model is described by the exponential equation SSA(t) = (SSA₀ −26.8) e^−0.54t + 26.8. The model performance is optimal for daily mean values of surface temperature in the range 0 °C to −25 °C and wind speed < 6 m s⁻¹. The findings from the SSA analysis are used to explore the influence of surface snow metamorphism on altering the isotopic composition of surface snow. It is found that rapid SSA decay events correspond to decreases in d-excess over a 2-day period in 72 % of the samples. Detailed studies using Empirical Orthogonal Function (EOF) analysis revealed a coherence between the dominant mode of variance of SSA and d-excess during periods of low spatial variability of surface snow over the sampling transect, suggesting that processes driving change in SSA also influence d-excess. Our findings highlight the need for future studies to decouple the processes driving surface snow metamorphism in order to quantify the fractionation effect of individual processes on the snow isotopic composition.