Sticky Snow: Combining Snow and Radiative Transfer Models in the Percolation Area of the Greenland Ice Sheet

• Main Author: Kreynen, Dylan (author)
• Other Authors: Lhermitte, S.L.M. (mentor), Li, W. (mentor), Slobbe, D.C. (mentor), Delft University of Technology (degree granting institution)
• Format: Master Thesis
• Language: English
• Published: 2021
• Subjects: microwave; remote sensing; brightness temperature; snow model; Greenland; ice sheet; radiative transfer; GRS; Ice Sheet
• Online Access: http://resolver.tudelft.nl/uuid:14d5d85d-2346-4428-9b40-d2495ef6def9

Subsurface firn processes play a crucial role in ice sheet mass loss mechanisms. On Greenland surface meltwater percolates to deeper layers where porous firn retains it, directly inhibiting runoff. However, secondary effects such as the formation of impermeable ice slabs may indirectly and irreversibly accelerate runoff and with it global sea level rise. Microwave remote sensing offers opportunities to monitor these processes, but due to the simplicity of their underlying snow models retrieval methods fail over areas subject to melt and refreezing - areas where the firn's (in)ability to buffer meltwater is critical. This study presents a new forward model which given initial conditions and atmospheric forcing first solves for the firn state through a full-complexity snow model (SNOWPACK) and then simulates multifrequency brightness temperature (Tb) time series (using radiative transfer model SMRT). As part of a comprehensive sensitivity analysis three ensembles of multi-decade Tb time series (19 and 37 GHz) were modelled for the DYE-2 site in the percolation area of the Greenland Ice Sheet. Model performance based on RMSE w.r.t. independent Tb satellite observations was found to be sensitive to biases introduced in the atmospheric forcing record (with air temperature, precipitation and relative humidity controlling variance) and snow model settings (new snow grain size and albedo settings) and not to initial profile conditions. However, computed RMSEs were high (min. 17.8 K at 37V and 19.4 K at 19V) due to trends in modelled Tb consistently underestimating observed trends when taken over an accumulation season. It is shown that this can only be explained by the constant-with-time stickiness assumption used to link the snow model’s microstructure representation to the sticky hard sphere model employed for the radiative transfer scheme. A seasonal stickiness signal is made evident for the conditions at DYE-2 and linked to its yearly melt-refreeze-accumulation cycle. These results demonstrate that earlier ...