Snow Water Equivalent Estimates Research Articles

Advancing terrestrial snow depth monitoring with machine learning and L-band InSAR data: a case study using NASA’s SnowEx 2017 data

Current terrestrial snow depth mapping from space faces challenges in spatial coverage, revisit frequency, and cost. Airborne lidar, although precise, incurs high costs and has limited geographical coverage, thereby necessitating the exploration of alternative, cost-effective methodologies for snow depth estimation. The forthcoming NASA-ISRO Synthetic Aperture Radar (NISAR) mission, with its 12-day global revisit cycle and 1.25 GHz L-band frequency, introduces a promising avenue for cost-effective, large-scale snow depth and snow water equivalent (SWE) estimation using L-band Interferometric SAR (InSAR) capabilities. This study demonstrates InSAR’s potential for snow depth estimation via machine learning. Using 3 m resolution L-band InSAR products over Grand Mesa, Colorado, we compared the performance of three machine learning approaches (XGBoost, ExtraTrees, and Neural Networks) across open, vegetated, and the combined (open + vegetated) datasets using Root Mean Square Error (RMSE), Mean Bias Error (MBE), and R2 metrics. XGBoost emerged as the superior model, with RMSE values of 9.85 cm, 10.46 cm, and 9.88 cm for open, vegetated, and combined regions, respectively. Validation against in situ snow depth measurements resulted in an RMSE of approximately 16 cm, similar to in situ validation of the airborne lidar. Our findings indicate that L-band InSAR, with its ability to penetrate clouds and cover extensive areas, coupled with machine learning, holds promise for enhancing snow depth estimation. This approach, especially with the upcoming NISAR launch, may enable high-resolution (∼10 m) snow depth mapping over extensive areas, provided suitable training data are available, offering a cost-effective approach for snow monitoring. The code and data used in this work are available at https://github.com/cryogars/uavsar-lidar-ml-project.

The effect of hydrological model structure on spring flow forecasts when assimilating a distributed snow product

Accurate forecasting of spring flow is essential for mitigating flood damage and optimizing hydroelectric power production. In northern countries such as Canada, this flow is mainly driven by snowmelt processes. By integrating snowpack data from diverse sources (in situ, remote sensing, and reanalysis) with modelled snow-related state variables through data assimilation (DA), it is possible to leverage both modeling and observations for more accurate spring flow estimates. Accurate estimates of snow water equivalent (SWE) within a heterogeneous snowpack are crucial for optimizing the advantages of snow DA. Here we assess the potential effect of distributed SNOw Data Assimilation System (SNODAS) SWE data on improving spring flow in two hydrological models having distinct inner structures: HSAMI, a lumped model, and HYDROTEL, a distributed model. DA analyses used an ensemble Kalman filter scheme, which was run over a three-year period. The results were then compared with those from an open-loop experiment. We used Nash–Sutcliffe efficiency (NSE) and bias to assess the results and found that SNODAS SWE DA did not improve 1-day spring flow forecasts for the lumped model. In contrast, HYDROTEL produced more accurate 1-day spring flow estimates for all 3 years, improving NSE of the spring flow forecasts from 0.52 to 0.70, 0.32 to 0.68, and 0.39 to 0.67 for the 2014–2017 period. This study demonstrates that incorporating distributed snow data into a distributed hydrological model can improve spring flow forecasts. Given the availability of SNODAS dataset over most Canadian watersheds, this DA framework is aimed to serve as an asset to enhance operational flood forecasting systems across Canada.

Use of multiple reference data sources to cross-validate gridded snow water equivalent products over North America

Abstract. We use snow course and airborne gamma data available over North America to compare the validation of gridded snow water equivalent (SWE) products when evaluated with one reference dataset versus the other. We assess product performance across both non-mountainous and mountainous regions, determining the sensitivity of relative product rankings and absolute performance measures. In non-mountainous areas, product performance is insensitive to the choice of SWE reference dataset (snow course or airborne gamma): the validation statistics (bias, unbiased root mean squared error, and correlation) are consistent with one another. In mountainous areas, the choice of reference dataset has little impact on relative product ranking but a large impact on assessed error magnitudes (bias and unbiased root mean squared error). Further analysis indicates the agreement in non-mountainous regions occurs because the reference SWE estimates themselves agree up to spatial scales of at least 50 km, comparable to the grid spacing of most available SWE products. In mountain areas, there is poor agreement between the reference datasets, even at short distances (< 5 km). We determine that differences in assessed error magnitudes result primarily from the range of SWE magnitudes sampled by each method, although their respective spatiotemporal distribution and elevation differences between the reference measurements and grid centroids also play a role. We use this understanding to produce a combined reference SWE dataset for North America, applicable to future gridded SWE product evaluations and other applications.

Characterization of non-Gaussianity in the snow distributions of various landscapes

Abstract. Seasonal snowpack is an important predictor of the water resources available in the following spring and early-summer melt season. Total basin snow water equivalent (SWE) estimation usually requires a form of statistical analysis that is implicitly built upon the Gaussian framework. However, it is important to characterize the non-Gaussian properties of snow distribution for accurate large-scale SWE estimation based on remotely sensed or sparse ground-based observations. This study quantified non-Gaussianity using sample negentropy; the Kullback–Leibler divergence from the Gaussian distribution for field-observed snow depth data from the North Slope, Alaska; and three representative SWE distributions in the western USA from the Airborne Snow Observatory (ASO). Snowdrifts around lakeshore cliffs and deep gullies can bring moderate non-Gaussianity in the open, lowland tundra of North Slope, Alaska, while the ASO dataset suggests that subalpine forests may effectively suppress the non-Gaussianity of snow distribution. Thus, non-Gaussianity is found in areas with partial snow cover and wind-induced snowdrifts around topographic breaks on slopes and on other steep terrain features. The snowpacks may be considered weakly Gaussian in coastal regions with open tundra in Alaska and alpine and subalpine terrains in the western USA if the land is completely covered by snow. The wind-induced snowdrift effect can potentially be partitioned from the observed snow spatial distribution guided by its Gaussianity.

Retrieval of snow and soil properties for forward radiative transfer modeling of airborne Ku-band SAR to estimate snow water equivalent: the Trail Valley Creek 2018/19 snow experiment

Abstract. Accurate snow information at high spatial and temporal resolution is needed to support climate services, water resource management, and environmental prediction services. However, snow remains the only element of the water cycle without a dedicated Earth observation mission. The snow scientific community has shown that Ku-band radar measurements provide quality snow information with its sensitivity to snow water equivalent and the wet/dry state of snow. With recent developments of tools like the snow micropenetrometer (SMP) to retrieve snow microstructure data in the field and radiative transfer models like the Snow Microwave Radiative Transfer (SMRT) model, it becomes possible to properly characterize the snow and how it translates into radar backscatter measurements. An experiment at Trail Valley Creek (TVC), Northwest Territories, Canada, was conducted during the winter of 2018/19 in order to characterize the impacts of varying snow geophysical properties on Ku-band radar backscatter at a 100 m scale. Airborne Ku-band data were acquired using the University of Massachusetts radar instrument. This study shows that it is possible to calibrate SMP data to retrieve statistical information on snow geophysical properties and properly characterize a representative snowpack at the experiment scale. The tundra snowpack measured during the campaign can be characterize by two layers corresponding to a rounded snow grain layer and a depth hoar layer. Using RADARSAT-2 and TerraSAR-X data, soil background roughness properties were retrieved (msssoil=0.010±0.002), and it was shown that a single value could be used for the entire domain. Microwave snow grain size polydispersity values of 0.74 and 1.11 for rounded and depth hoar snow grains, respectively, were retrieved. Using the geometrical optics surface backscatter model, the retrieved effective soil permittivity increased from C-band (εsoil=2.47) to X-band (εsoil=2.61) and to Ku-band (εsoil=2.77) for the TVC domain. Using the SMRT and the retrieved soil and snow parameterizations, an RMSE of 2.6 dB was obtained between the measured and simulated Ku-band backscatter values when using a global set of parameters for all measured sites. When using a distributed set of soil and snow parameters, the RMSE drops to 0.9 dB. This study thus shows that it is possible to link Ku-band radar backscatter measurements to snow conditions on the ground using a priori knowledge of the snow conditions to retrieve snow water equivalent (SWE) at the 100 m scale.

Detecting snowfall events over the Arctic using optical and microwave satellite measurements

Abstract. The precipitation over the Arctic region is a difficult quantity to determine with high accuracy, as the in situ observation network is sparse, and current climate models, atmospheric reanalyses, and direct satellite-based precipitation observations suffer from diverse difficulties that hinder the correct assessment of precipitation. We undertake a proof-of-concept investigation into how accurately optical satellite observations, namely Sentinel-2 surface-reflectance-based grain-size-connected specific surface area of snow (SSA), and microwave-based snow water equivalent (SWE) estimates can detect snowfall over the Arctic. In addition to the satellite data, we also include ERA5-Land SWE data to support the analysis. Here, we chose a limited area (a circle of 100 km radius around Luosto radar located in Northern Finland) and a short time period (covering March 2018) to test these data sources and their usability in this precipitation assessment problem. We classified differences between observations independently for SSA and SWE and compared the results to the radar-based snowfall information. These initial results are promising. Situations with snowfall are classified with high recalls, 64 % for the satellite-based SWE, 77 % for ERA5-Land-based SWE, and around 90 % for SSA compared to radar-based data. Cases without snowfall are more difficult to classify correctly using satellite-based data. The recall values are 34 % for satellite-based SWE and vary from almost 60 % to over 70 % for SSA. SWE from ERA5-Land has the highest recall value for cases without snowfall, 80 %. These results indicate that optical and microwave-based satellite observations can be used to detect snowfall events over the Arctic.

Machine Learning Enhanced by Feature Engineering for Estimating Snow Water Equivalent

This study compares the calculation of snow water equivalent (SWE) using machine learning algorithms with the conventional degree-day method. The study uses machine learning techniques such as LASSO, Random Forest, Support Vector Machines, and CatBoost. It proposes an innovative use of feature engineering (FE) to improve the accuracy and robustness of SWE predictions by machine learning intended for interpolation, extrapolation, or imputation of missing data. The performance of machine learning approaches is evaluated against the traditional degree-day method for predicting SWE. The study emphasizes and demonstrates gains when modeling is enhanced by transforming basic, raw data through feature engineering. The results, verified in a case study from the mountainous region of Slovakia, suggest that machine learning, particularly CatBoost with feature engineering, shows better results in SWE estimation in comparison with the degree-day method, although the authors present a refined application of the degree-day method by utilizing genetic algorithms. Nevertheless, the study finds that the degree-day method achieved accuracy with a Nash–Sutcliffe coefficient of efficiency NSE = 0.59, while the CatBoost technique enhanced with the proposed FE achieved an accuracy NSE = 0.86. The results of this research contribute to refining snow hydrology modeling and optimizing SWE prediction for improved decision-making in snow-dominated regions.

Bootstrap Methods for Bias-Correcting Probability Distribution Parameters Characterizing Extreme Snow Accumulations

Accurately quantifying the threat of collapse due to the weight of settled snow on the roof of a structure is crucial for ensuring structural safety. This quantification relies upon direct measurements of the snow water equivalent (SWE) of settled snow, though most weather stations in the United States only measure snow depth. The absence of direct load measurements necessitates the use of modeled estimates of SWE, which often results in the underestimation of the scale/variance parameter of the distribution of annual maximum SWE. This paper introduces a novel bias correction method that employs a bootstrap technique with regression-based models to calibrate the variance parameter of the distribution. The efficacy of this approach is demonstrated on real and simulated datasets. The findings reveal varied levels of success, with the efficacy of the proposed approach being inherently dependent on the quality of the selected regression-based model. These findings demonstrate that integrating our approach with a suitable regression-based model can produce unbiased or nearly unbiased annual maximum SWE distribution parameters in the absence of direct SWE measurements.

Measuring prairie snow water equivalent with combined UAV-borne gamma spectrometry and lidar

Abstract. Despite decades of effort, there remains an inability to measure snow water equivalent (SWE) at high spatial resolutions using remote sensing. Passive gamma ray spectrometry is one of the only well-established methods to reliably remotely sense SWE, but airborne applications to date have been limited to observing kilometre-scale areal averages. Noting the increasing capabilities of unoccupied aerial vehicles (UAVs) and miniaturization of passive gamma ray spectrometers, this study tested the ability of a UAV-borne gamma spectrometer and concomitant UAV-borne lidar to quantify the spatial variability of SWE at high spatial resolutions. Gamma and lidar observations from a UAV (UAV-gamma and UAV-lidar) were collected over two seasons from shallow, wind-blown, prairie snowpacks in Saskatchewan, Canada, with validation data collected from manual snow depth and density observations. A fine-resolution (0.25 m) reference dataset of SWE, to test UAV-gamma methods, was developed from UAV-lidar snow depth and snow survey snow density observations. The ability of UAV-gamma to resolve the areal average and spatial variability of SWE was promising with appropriate flight characteristics. Survey flights flown at a velocity of 5 m s−1, altitude of 15 m, and line spacing of 15 m were unable to capture the average or spatial variability of SWE within the uncertainty of the reference dataset. Slower, lower, and denser flight lines at a velocity of 4 m s−1, altitude of 8 m, and line spacing of 8 m were able to successfully observe areal average SWE and its variability at spatial resolutions greater than 22.5 m. Using a combination of UAV-based gamma SWE and UAV-based lidar snow depth improved the spatial representation of SWE substantially and permitted estimation of SWE at a spatial resolution 0.25 m with a ± 14.3 mm error relative to the reference SWE dataset. UAV-borne gamma spectrometry to estimate SWE is a promising and novel technique that has the potential to improve the measurement of shallow prairie snowpacks, and when combined with UAV-borne lidar snow depths, can provide fine-resolution, high-accuracy estimates of prairie SWE. Research on optimal hardware, data processing, and interpolation techniques is called for to further improve this remote sensing product and explore its application in other environments.

Quantifying weather-induced unreliable public transportation service in cold regions under future climate model scenarios

Climate change, particularly in cold regions, significantly challenges public transportation systems. This study conducts a comprehensive analysis of weather patterns and public transit reliability in the context of climate change impacts. Leveraging advanced modeling techniques, including a ridge regression model for snow water equivalent data estimation and a long short-term memory (LSTM) based on recurrent neural network, the study aims to assess the reliability trends of the rapid transit system under various climate scenarios. The findings reveal that climate change in general increases weather-related delays in the Toronto transit system. The number of short delays decreased accordingly due to changes in winter temperatures but exacerbated long delays as the number of weather extremes increased. The LSTM model performed effectively in predicting delays, especially for the rapid transit system sensitive to weather variations. This study emphasizes the need for robust planning and interventions to increase the resilience of transit systems against climate change and highlights the importance of the integration of climate and extreme weather considerations into transportation management.

Quantifying regional variability of machine-learning-based snow water equivalent estimates across the Western United States

Seasonal snow-derived water is a critical component of the water supply in the mountains and downstream regions, and the accurate characterization of available water in the form of snow-water-equivalent (SWE), peak SWE, and snowmelt onset are essential inputs for water management efforts. Arising from recent advancements in artificial intelligence (AI) and machine learning (ML), we introduce a large-scale ML SWE model leveraging publicly available data sources and open-source software. The model demonstrates the application of a limited feature space in a relatively simple ML architecture without the need for process-based formulations to effectively estimate spatially continuous SWE at a daily temporal resolution. Beginning with in situ SWE measurements (i.e., SNOTEL), lidar-derived terrain features, and temporal variables, we employ localized feature engineering and optimization via gradient-boosting decision trees to identify regionally unique drivers of snowpack dynamics and use the optimal features to train regionally independent artificial neural networks to estimate regional SWE at a 1 km spatial resolution. The model results yield respectable skill in reconstructed 1 km gridded SWE magnitudes in a hindcast simulation of the 2019 water year that is independent of the training and testing data. Comparing model estimates to over 6200 observations, the model demonstrates a weighted RMSE of 15.4 cm, Kling-Gupta Efficiency metric of 0.86, and a percent bias of 0.71% across 23 snow-influenced regions in the western U.S. The model simulation produces peak SWE estimates within 10 cm for twenty of the twenty-three regions, demonstrating capability in effectively capturing regional snow accumulation processes. The demonstration of low-error ML workflows capable of providing near-real-time, spatially continuous SWE estimates at a high spatial resolution provides proof-of-concept and a foundation to effectively update snow state variables that drive water supply forecasts in snow-dominated regions.

Snow water equivalent retrieved from X- and dual Ku-band scatterometer measurements at Sodankylä using the Markov Chain Monte Carlo method

Abstract. Radar at high frequency is a promising technique for fine-resolution snow water equivalent (SWE) mapping. In this paper, we extend the Bayesian-based Algorithm for SWE Estimation (BASE) from passive to active microwave (AM) application and test it using ground-based backscattering measurements at three frequencies (X and dual Ku bands; 10.2, 13.3, and 16.7 GHz), with VV polarization obtained at a 50° incidence angle from the Nordic Snow Radar Experiment (NoSREx) in Sodankylä, Finland. We assumed only an uninformative prior for snow microstructure, in contrast with an accurate prior required in previous studies. Starting from a biased monthly SWE prior from land surface model simulation, two-layer snow state variables and single-layer soil variables were iterated until their posterior distribution could stably reproduce the observed microwave signals. The observation model is the Microwave Emission Model of Layered Snowpacks 3 and Active (MEMLS3&a) based on the improved Born approximation. Results show that BASE-AM achieved an RMSE of ∼ 10 cm for snow depth and less than 30 mm for SWE, compared with the RMSE of ∼ 20 cm snow depth and ∼ 50 mm SWE from priors. Retrieval errors are significantly larger when BASE-AM is run using a single snow layer. The results support the potential of X- and Ku-band radar for SWE retrieval and show that the role of a precise snow microstructure prior in SWE retrieval may be substituted by an SWE prior from exterior sources.

Bayesian physical–statistical retrieval of snow water equivalent and snow depth from X- and Ku-band synthetic aperture radar – demonstration using airborne SnowSAr in SnowEx'17

Abstract. A physical–statistical framework to estimate snow water equivalent (SWE) and snow depth from synthetic aperture radar (SAR) measurements is presented and applied to four SnowSAR flight-line data sets collected during the SnowEx'2017 field campaign in Grand Mesa, Colorado, USA. The physical (radar) model is used to describe the relationship between snowpack conditions and volume backscatter. The statistical model is a Bayesian inference model that seeks to estimate the joint probability distribution of volume backscatter measurements, snow density and snow depth, and physical model parameters. Prior distributions are derived from multilayer snow hydrology predictions driven by downscaled numerical weather prediction (NWP) forecasts. To reduce the signal-to-noise ratio, SnowSAR measurements at 1 m resolution were upscaled by simple averaging to 30 and 90 m resolution. To reduce the number of physical parameters, the multilayer snowpack is transformed for Bayesian inference into an equivalent one- or two-layer snowpack with the same snow mass and volume backscatter. Successful retrievals meeting NASEM (2018) science requirements are defined by absolute convergence backscatter errors ≤1.2 dB and local SnowSAR incidence angles between 30 and 45∘ for X- and Ku-band VV-pol backscatter measurements and were achieved for 75 % to 87 % of all grassland pixels with SWE up to 0.7 m and snow depth up to 2 m. SWE retrievals compare well with snow pit observations, showing strong skill in deep snow with average absolute SWE residuals of 5 %–7 % (15 %–18 %) for the two-layer (one-layer) retrieval algorithm. Furthermore, the spatial distributions of snow depth retrievals vis-à-vis lidar estimates have Bhattacharya coefficients above 94 % (90 %) for homogeneous grassland pixels at 30 m (90 m resolution), and values up to 76 % in mixed forest and grassland areas, indicating that the retrievals closely capture snowpack spatial variability. Because NWP forecasts are available everywhere, the proposed approach could be applied to SWE and snow depth retrievals from a dedicated global snow mission.

Assessing the impact of distributed snow water equivalent calibration and assimilation of Copernicus snow water equivalent on modelled snow and streamflow performance

AbstractAccuracy of the Copernicus snow water equivalent (SWE) product and the impact of SWE calibration and assimilation on modelled SWE and streamflow was evaluated. Daily snowpack measurements were made at 12 locations from 2016 to 2019 across a 4104 km2 mixed‐forest basin in the Great Lakes region of central Ontario, Canada. Sub‐basin daily SWE calculated from these sites, observed discharge, and lake levels were used to calibrate a hydrologic model developed using the Raven modelling framework. Copernicus SWE was bias corrected during the melt period using mean bias subtraction and was compared to daily basin average SWE calculated from the measured data. Bias corrected Copernicus SWE was assimilated into the models using a range of parameters and the parameterizations from the model calibration. The bias corrected Copernicus product agreed well with measured data and provided a good estimate of mean basin SWE demonstrating that the product shows promise for hydrology applications within the study region. Calibration to spatially distributed SWE substantially improved the basin scale SWE estimate while only slightly degrading the flow simulation demonstrating the value of including SWE in a multi‐objective calibration formulation. The particle filter experiments yielded the best SWE estimation but moderately degraded the flow simulation. The particle filter experiments constrained by the calibrated snow parameters produced similar results to the experiments using the upper and lower bounds indicating that, in this study, model calibration prior to assimilation was not valuable. The calibrated models exhibited varying levels of skill in estimating SWE but demonstrated similar streamflow performance. This indicates that basin outlet streamflow can be accurately estimated using a model with a poor representation of distributed SWE. This may be sufficient for applications where estimating flow is the primary water management objective. However, in applications where understanding the physical processes of snow accumulation, melt and streamflow generation are important, such as assessing the impact of climate change on water resources, accurate representations of SWE are required and can be improved via multi‐objective calibration or data assimilation, as demonstrated in this study.

Combining Hydrological Models and Remote Sensing to Characterize Snowpack Dynamics in High Mountains

Seasonal snowpacks, characterized by their snow water equivalent (SWE), can play a major role in the hydrological cycle of montane environments with months of snow accretion followed by episodes of melt controlling flood risk and water resource availability downstream. Quantifying the temporal and spatial patterns of snowpack accumulation and its subsequent melt and runoff is an internationally significant challenge, particularly within mountainous regions featuring complex terrain with limited or absent observational data. Here we report a new approach to snowpack characterization using open-source global satellite and modelled data products (precipitation and SWE) greatly enhancing the utility of the widely used Soil and Water Assessment Tool (SWAT). The paper focusses on the c. 23,000 km2 Chenab river basin (CRB) in the headwaters of the Indus Basin, globally important because of its large and growing population and increasing water insecurity due to climate change. We used five area-weighted averaged satellite, gridded and reanalysis precipitation datasets: ERA5-Land, CMORPH, TRMM, APHRODITE and CPC UPP. As well as comparison to local weather station data, these were used in SWAT to model streamflow for evaluation against observed streamflow at the basin outlet. ERA5-Land data provided the best streamflow match-ups and was used to infer snowpack (SWE) dynamics at basin and sub-basin scales. Snow reference data were derived from remote sensing and modelled SWE re-analysis products: ULCA-SWE and KRA-SWE, respectively. Beyond conventional auto-calibration and single-variable approaches we undertook multi-variable calibration using R-SWAT to manually adjust snow parameters alongside observed streamflow data. Characterization of basin-wide patterns of snowpack build-up and melt (SWE dynamics) were greatly strengthened using KRA-SWE data accompanied by improved streamflow simulation in sub-basins dominated by seasonal snow cover. UCLA-SWE data also improved SWE estimations using R-SWAT but weakened the performance of simulated streamflow due to under capture of seasonal runoff from permanent snow/ice fields in the CRB. This research highlights the utility and value of remote sensing and modelling data to drive better understanding of snowpack dynamics and their contribution to runoff in the absence of in situ snowpack data in high-altitude environments. An improved understanding of snow-bound water is vital in natural hazard risk assessment and in better managing worldwide water resources in the populous downstream regions of mountain-fed large rivers under threat from climate change.

Combining Differential SAR Interferometry and Copolar Phase Differences for Snow Water Equivalent Estimation

Combining Differential SAR Interferometry and Copolar Phase Differences for Snow Water Equivalent Estimation

Six Consecutive Seasons of High‐Resolution Mountain Snow Depth Maps From Satellite Stereo Imagery

AbstractFine‐scale seasonal snow depth observations can improve estimates of snow water equivalent at critical times of year. Airborne lidar is the current gold standard for snow depth measurement, but it involves high costs and relatively limited coverage. Using very‐high‐resolution satellite stereo images from WorldView‐2, WorldView‐3, and Pléiades‐HR 1A/1B, we produced a six‐year time series (2017–2022) of spatially continuous digital elevation models for an 874 km2 study area over Grand Mesa, Colorado. We generated high‐resolution stereo snow depth maps that capture intra‐ and interannual variability and span multiple anomalous years (58%–158% of median peak SNOTEL snow depth). Comparisons with near‐contemporaneous airborne lidar snow depth measurements showed good agreement, with median offset of −0.13 m, precision of 0.19 m and accuracy of 0.31 m. Our results suggest that satellite stereo can provide snow depth observations with the spatiotemporal coverage needed to improve operational forecast models and inform adaptive management strategies.

Evaluation of snow cover properties in ERA5 and ERA5-Land with several satellite-based datasets in the Northern Hemisphere in spring 1982–2018

Abstract. Seasonal snow cover of the Northern Hemisphere (NH) greatly influences surface energy balance; hydrological cycle; and many human activities, such as tourism and agriculture. Monitoring snow cover at a continental scale is only possible from satellites or using reanalysis data. This study aims to analyze the time series of snow water equivalent (SWE), snow cover extent (SCE), and surface albedo in spring in ERA5 and ERA5-Land reanalysis data and to compare the time series with several satellite-based datasets. As reference data for the SWE intercomparison, we use bias-corrected SnowCCI v1 data for non-mountainous regions and the mean of Brown, MERRA-2, and Crocus v7 datasets for the mountainous regions. For surface albedo, we use the black-sky albedo datasets CLARA-A2 SAL, based on AVHRR data, and MCD43D51, based on MODIS data. Additionally, we use Rutgers and JAXA JASMES SCE products. Our study covers land areas north of 40∘ N and the period between 1982 and 2018 (spring season from March to May). The analysis shows that both ERA5 and ERA5-Land overestimate total NH SWE by 150 % to 200 % compared to the SWE reference data. ERA5-Land shows larger overestimation, which is mostly due to very high SWE values over mountainous regions. The analysis revealed a discontinuity in ERA5 around the year 2004 since adding the Interactive Multisensor Snow and Ice Mapping System (IMS) from the year 2004 onwards considerably improves SWE estimates but makes the trends less reliable. The negative NH SWE trends in ERA5 range from −249 to −236 Gt per decade in spring, which is 2 to 3 times larger than the trends detected by the other datasets (ranging from −124 to −77 Gt per decade). SCE is accurately described in ERA5-Land, whereas ERA5 shows notably larger SCE than the satellite-based datasets. Albedo estimates are more consistent between the datasets, with a slight overestimation in ERA5 and ERA5-Land. The negative trends in SCE and albedo are strongest in May, when the albedo trend varies from −0.011 to −0.006 per decade depending on the dataset. The negative SCE trend detected by ERA5 in May (-1.22×106 km2 per decade) is about twice as large as the trends detected by all other datasets (ranging from −0.66 to -0.50×106 km2 per decade). The analysis also shows that there is a large spatial variability in the trends, which is consistent with other studies.

Combined GNSS reflectometry–refractometry for automated and continuous in situ surface mass balance estimation on an Antarctic ice shelf

Abstract. Reliable in situ surface mass balance (SMB) estimates in polar regions are scarce due to limited spatial and temporal data availability. This study aims at deriving automated and continuous specific SMB time series for fast-moving parts of ice sheets and shelves (flow velocity > 10 m a−1) by developing a combined global navigation satellite system (GNSS) reflectometry and refractometry (GNSS-RR) method. In situ snow density, snow water equivalent (SWE), and snow deposition or erosion are estimated simultaneously as an average over an area of several square meters and independently on weather conditions. The combined GNSS-RR method is validated and investigated regarding its applicability to a moving, high-latitude ice shelf. A combined GNSS-RR system was therefore installed in November 2021 on the Ekström ice shelf (flow velocity ≈ 150 m a−1) in Dronning Maud Land, Antarctica. The reflected and refracted GNSS observations from the site are post-processed to obtain snow accumulation (deposition and erosion), SWE, and snow density estimates with a 15 min temporal resolution. The results of the first 16 months of data show a high level of agreement with manual and automated reference observations from the same site. Snow accumulation, SWE, and density are derived with uncertainties of around 9 cm, 40 kg m−2 a−1, and 72 kg m−3, respectively. This pilot study forms the basis for extending observational networks with GNSS-RR capabilities, particularly in polar regions. Regional climate models, local snow modeling, and extensive remote sensing data products will profit from calibration and validation based on such in situ time series, especially if many such sensors will be deployed over larger regional scales.

Assimilation of airborne gamma observations provides utility for snow estimation in forested environments

Abstract. An airborne gamma-ray remote-sensing technique provides a strong potential to estimate a reliable snow water equivalent (SWE) in forested environments where typical remote-sensing techniques have large uncertainties. This study explores the utility of assimilating the temporally (up to four measurements during a winter period) and spatially sparse airborne gamma SWE observations into a land surface model (LSM) to improve SWE estimates in forested areas in the northeastern US. Here, we demonstrate that the airborne gamma SWE observations add value to the SWE estimates from the Noah LSM with multiple parameterization options (Noah-MP) via assimilation despite the limited number of measurements. Improvements are witnessed during the snow accumulation period, while reduced skills are seen during the snowmelt period. The efficacy of the gamma data is greater for areas with lower vegetation cover fraction and topographic heterogeneity ranges, and it is still effective at reducing the SWE estimation errors for areas with higher topographic heterogeneity. The gamma SWE data assimilation (DA) also shows a potential to extend the impact of flight-line-based measurements to adjacent areas without observations by employing a localization approach. The localized DA reduces the modeled SWE estimation errors for adjacent grid cells up to 32 km distance from the flight lines. The enhanced performance of the gamma SWE DA is evident when the results are compared to those from assimilating the existing satellite-based SWE retrievals from the Advanced Microwave Scanning Radiometer 2 (AMSR2) for the same locations and time periods. Although there is still room for improvement, particularly for the melting period, this study shows that the gamma SWE DA is a promising method to improve the SWE estimates in forested areas.