Spatially distributed snow depth, bulk density, and snow water equivalent from ground-based and airborne sensor integration at Grand Mesa, Colorado, USA

Abstract

Abstract. Estimating snow mass in the mountains remains a major challenge for remote-sensing methods. Airborne lidar can retrieve snow depth, and some promising results have recently been obtained from spaceborne platforms, yet density estimates are required to convert snow depth to snow water equivalent (SWE). However, the retrieval of snow bulk density remains unsolved, and limited data are available to evaluate model estimates of density in mountainous terrain. Toward the goal of landscape-scale retrievals of snow density, we estimated bulk density and length-scale variability by combining ground-penetrating radar (GPR) two-way travel-time observations and airborne-lidar snow depths collected during the mid-winter NASA SnowEx 2020 campaign at Grand Mesa, Colorado, USA. Key advancements of our approach include an automated layer-picking method that leverages the GPR reflection coherence and the distributed lidar–GPR-retrieved bulk density with machine learning. The root-mean-square error between the distributed estimates and in situ observations is 11 cm for depth, 27 kg m−3 for density, and 46 mm for SWE. The median relative uncertainty in distributed SWE is 13 %. Interactions between wind, terrain, and vegetation display corroborated controls on bulk density that show model and observation agreement. Knowledge of the spatial patterns and predictors of density is critical for the accurate assessment of SWE and essential snow research applications. The spatially continuous snow density and SWE estimated over approximately 16 km2 may serve as necessary calibration and validation for stepping prospective remote-sensing techniques toward broad-scale SWE retrieval.