Firn Density Distribution and Annual Snow Water Equivalent Estimates from Ground Penetrating Radar

Accurate estimation of glacier volume-to-mass conversion relies on a thorough understanding of firn density, both in-depth and over time. Ground-penetrating radar (GPR) serves as a suitable geophysical tool to trace internal reflection horizons (IRHs) and estimate the physical properties of different layers. Our goal is to characterize the IRHs as annual layers and ascertain the spatial firn density-depth profile in the accumulation zone of the Aletsch glacier. The process involves identifying IRHs from radargrams and iteratively selecting the annual layers by excluding unreasonable layer structures. For an accurate estimation of firn density distribution, it is necessary to derive the velocity-depth profile of electromagnetic waves within the firn zone. The common mid-point (CMP) method was applied to track the velocity distribution within the firn body. Additionally, a method was introduced to estimate the velocity-depth profile for longer GPR profiles by backtracking the calculated velocity from the CMP gather. To validate IRHs as annual firn layers, we utilized annual accumulation measurements at a nearby stake for Snow Water Equivalent (SWE) estimation. The resulting firn density-depth profile was compared to different firn densification models, considering regional meteorological information. This approach enables us to determine a reliable density-depth function for bulk SWE computations. The study also addresses uncertainties associated with selecting IRHs as annual layers and enhances the application of local volume-to-mass estimates.