Snow Water Equivalent Estimates Research Articles

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.

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.

Snowpack relative permittivity and density derived from near‐coincident lidar and ground‐penetrating radar

AbstractDepth‐based and radar‐based remote sensing methods (e.g., lidar, synthetic aperture radar) are promising approaches for remotely measuring snow water equivalent (SWE) at high spatial resolution. These approaches require snow density estimates, obtained from in‐situ measurements or density models, to calculate SWE. However, in‐situ measurements are operationally limited, and few density models have seen extensive evaluation. Here, we combine near‐coincident, lidar‐measured snow depths with ground‐penetrating radar (GPR) two‐way travel times (twt) of snowpack thickness to derive >20 km of relative permittivity estimates from nine dry and two wet snow surveys at Grand Mesa, Cameron Pass, and Ranch Creek, Colorado. We tested three equations for converting dry snow relative permittivity to snow density and found the Kovacs et al. (1995) equation to yield the best comparison with in‐situ measurements (RMSE = 54 kg m−3). Variogram analyses revealed a 19 m median correlation length for relative permittivity and snow density in dry snow, which increased to >30 m in wet conditions. We compared derived densities with estimated densities from several empirical models, the Snow Data Assimilation System (SNODAS), and the physically based iSnobal model. Estimated and derived densities were combined with snow depths and twt to evaluate density model performance within SWE remote sensing methods. The Jonas et al. (2009) empirical model yielded the most accurate SWE from lidar snow depths (RMSE = 51 mm), whereas SNODAS yielded the most accurate SWE from GPR twt (RMSE = 41 mm). Densities from both models generated SWE estimates within ±10% of derived SWE when SWE averaged >400 mm, however, model uncertainty increased to >20% when SWE averaged <300 mm. The development and refinement of density models, particularly in lower SWE conditions, is a high priority to fully realize the potential of SWE remote sensing methods.

Spatially continuous snow depth mapping by aeroplane photogrammetry for annual peak of winter from 2017 to 2021 in open areas

Abstract. Information on snow depth and its spatial distribution is important for numerous applications, including natural hazard management, snow water equivalent estimation for hydropower, the study of the distribution and evolution of flora and fauna, and the validation of snow hydrological models. Due to its heterogeneity and complexity, specific remote sensing tools are required to accurately map the snow depth distribution in Alpine terrain. To cover large areas (>100 km2), airborne laser scanning (ALS) or aerial photogrammetry with large-format cameras is needed. While both systems require piloted aircraft for data acquisition, ALS is typically more expensive than photogrammetry but yields better results in forested terrain. While photogrammetry is slightly cheaper, it is limited due to its dependency on favourable acquisition conditions (weather, light conditions). In this study, we present photogrammetrically processed high-spatial-resolution (0.5 m) annual snow depth maps, recorded during the peak of winter over a 5-year period under different acquisition conditions over a study area around Davos, Switzerland. Compared to previously carried out studies, using the Vexcel UltraCam Eagle Mark 3 (M3) sensor improves the average ground sampling distance to 0.1 m at similar flight altitudes above ground. This allows for very detailed snow depth maps in open areas, calculated by subtracting a snow-off digital terrain model (DTM, acquired with ALS) from the snow-on digital surface models (DSMs) processed from the airborne imagery. Despite challenging acquisition conditions during the recording of the UltraCam images (clouds, shaded areas and fresh snow), 99 % of unforested areas were successfully photogrammetrically reconstructed. We applied masks (high vegetation, settlements, water, glaciers) to increase the reliability of the snow depth calculations. An extensive accuracy assessment was carried out using check points, the comparison to DSMs derived from unpiloted aerial systems and the comparison of snow-free DSM pixels to the ALS DTM. The results show a root mean square error of approximately 0.25 m for the UltraCam X and 0.15 m for the successor, the UltraCam Eagle M3. We developed a consistent and reliable photogrammetric workflow for accurate snow depth distribution mapping over large regions, capable of analysing snow distribution in complex terrain. This enables more detailed investigations on seasonal snow dynamics and can be used for numerous applications related to snow depth distribution, as well as serving as a ground reference for new modelling approaches and satellite-based snow depth mapping.

Exploring the potential of thermal infrared remote sensing to improve a snowpack model through an observing system simulation experiment

Abstract. The assimilation of data from Earth observation satellites into numerical models is considered to be the path forward to estimate snow cover distribution in mountain catchments, providing accurate information on the mountainous snow water equivalent (SWE). The land surface temperature (LST) can be observed from space, but its potential to improve SWE simulations remains underexplored. This is likely due to the insufficient temporal or spatial resolution offered by the current thermal infrared (TIR) missions. However, three planned missions will provide global-scale TIR data at much higher spatiotemporal resolution in the coming years. To investigate the value of TIR data to improve SWE estimation, we developed a synthetic data assimilation (DA) experiment at five snow-dominated sites covering a latitudinal gradient in the Northern Hemisphere. We generated synthetic true LST and SWE series by forcing an energy balance snowpack model with the ERA5-Land reanalysis. We used this synthetic true LST to recover the synthetic true SWE from a degraded version of ERA5-Land. We defined different observation scenarios to emulate the revisiting times of Landsat 8 (16 d) and the Thermal infraRed Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA) (3 d) while accounting for cloud cover. We replicated the experiments 100 times at each experimental site to assess the robustness of the assimilation process with respect to cloud cover under both revisiting scenarios. We performed the assimilation using two different approaches: a sequential scheme (particle filter) and a smoother (particle batch smoother). The results show that LST DA using the smoother reduced the normalized root mean square error (nRMSE) of the SWE simulations from 61 % (open loop) to 17 % and 13 % for 16 d revisit and 3 d revisit respectively in the absence of clouds. We found similar but higher nRMSE values by removing observations due to cloud cover but with a substantial increase in the standard deviation of the nRMSE of the replicates, highlighting the importance of revisiting times in the stability of the assimilation performance. The smoother largely outperformed the particle filter algorithm, suggesting that the capability of a smoother to propagate the information along the season is key to exploit LST information for snow modelling. Finally, we have compared the benefit of assimilating LST with synthetic observations of fractional snow cover area (FSCA). LST DA performed better than FSCA DA in all the study sites, suggesting that the information provided by LST is not limited to the duration of the snow season. These results suggest that the LST data assimilation has an underappreciated potential to improve snowpack simulations and highlight the value of upcoming TIR missions to advance snow hydrology.

Assessment of L-band InSAR snow estimation techniques over a shallow, heterogeneous prairie snowpack

Snow water equivalent (SWE) is a critical input for weather, climate, and water resource management models at local to global scales. Despite its importance, global SWE measurements that are accurate, consistent, and at sufficiently high spatiotemporal resolutions are not currently available. L-band interferometric synthetic aperture radar (InSAR) techniques have been used to measure SWE at local to regional scales, and two upcoming L-band SAR satellite missions have renewed interest in these techniques to provide regular SWE measurements at the global scale. However, previous research demonstrating the capabilities of L-band InSAR-SWE measurement has been limited to mountain or tundra snowpack regimes. Here we examine the feasibility of applying the same techniques over a prairie snowpack, which are typically characterized by shallow snow depths (mean snow depth of 0.22 m in this study), exposed agricultural vegetation, and high spatial variability over short distances. Our study area in central Montana, USA (47.060, -109.951) was a validation site for NASA SnowEx 2021, as part of the UAVSAR snow timeseries. Airborne L-band SAR imagery was acquired by the UAVSAR platform while concurrent snow measurements were collected using uncrewed aerial vehicle (UAV)-based LiDAR, UAV-based photogrammetry, and ground-based manual techniques. This validation dataset enables an investigation of the effects of sub-pixel snow cover heterogeneity and exposed agricultural vegetation stubble on SAR data and the resulting SWE estimations. Results based on repeated application of the Kolmogorov–Smirnov test show that UAVSAR VV phase change is sensitive to differences in snow cover but relatively unaffected by differences in agricultural stubble height. However, we did not find similarly definitive results when we used the same phase change data to estimate SWE. Although broad spatial patterns were similar in both LiDAR-derived and InSAR-derived SWE estimates, considerable differences in the two estimates were apparent in areas with large sub-pixel snow depth variability. Our results indicate that additional work is necessary to derive accurate SWE estimates in prairie environments. Regular measurements from L-band SAR satellites will provide an excellent opportunity to refine InSAR-based snow estimation techniques over shallow, heterogeneous snowpacks.

A Multifactor Eigenvector Spatial Filtering-Based Method for Resolution-Enhanced Snow Water Equivalent Estimation in the Western United States

Accurate snow water equivalent (SWE) products are vital for monitoring hydrological processes and managing water resources effectively. However, the coarse spatial resolution (typically at 25 km from passive microwave remote sensing images) of the existing SWE products cannot meet the needs of explicit hydrological modeling. Linear regression ignores the spatial autocorrelation (SA) in the variables, and the measure of SA in the data assimilation algorithm is not explicit. This study develops a Resolution-enhanced Multifactor Eigenvector Spatial Filtering (RM-ESF) method to estimate daily SWE in the western United States based on a 6.25 km enhanced-resolution passive microwave record. The RM-ESF method is based on a brightness temperature gradience algorithm, incorporating not only factors including geolocation, environmental, topographical, and snow features but also eigenvectors generated from a spatial weights matrix to take SA into account. The results indicate that the SWE estimation from the RM-ESF method obviously outperforms other SWE products given its overall highest correlation coefficient (0.72) and lowest RMSE (56.70 mm) and MAE (43.88 mm), compared with the AMSR2 (0.33, 131.38 mm, and 115.45 mm), GlobSnow3 (0.50, 100.03 mm, and 83.58 mm), NCA-LDAS (0.48, 98.80 mm, and 81.94 mm), and ERA5 (0.65, 67.33 mm, and 51.82 mm), respectively. The RM-ESF model considers SA effectively and estimates SWE at a resolution of 6.25 km, which provides a feasible and efficient approach for SWE estimation with higher precision and finer spatial resolution.

Estimation of the Snow Water Equivalent Using Muon Scattering Radiography

AbstractDespite the important hydrological and ecological implications of the snowpack, its real time monitoring remains challenging. This is particularly relevant in relation to the Snow Water Equivalent (SWE), as the available technologies which measure it, exhibit a number of limitations that hinder their operational implementation. In this work, we explore the potential of a new technology, Muon Scattering Radiography, to infer the SWE. We coupled snowpack simulations generated by the SNOWPACK model, with a muon scattering simulation program based on GEANT4. The SWE is modeled as a function of the muon scattering distributions. Predictions of the SWE along the year are provided showing a root‐mean‐square error of 12 mm for 5 hr continuous measurements. We also performed laboratory measurements using ice samples, confirming the SWE estimation capabilities and the potential of the technique to operate as a SWE monitoring tool.

Estimation of snow bulk density and snow water equivalent on the Tibetan Plateau using snow cover duration and snow depth

Study regionTibetan Plateau (TP), China. Study focusSnow water equivalent (SWE) measurements is few compared with snow depth (SD) but combining with snow bulk density (SBD) can convert SD to SWE. No robust model has been established to estimate SBD and SWE for the TP. In this study, we aims to develop a model capable of obtaining SBD and SWE on the TP based on a significant dataset of simultaneously measured SD and SWE over nearly four decades. New hydrological insights for the regionWe found unique seasonal characteristics and elevation differences in the SBD on the TP. The newly introduced factor snow cover duration (SCD) can reflect the beginning and status of each snowfall-snowmelt cycle in this type of discontinuous snow region. The developed artificial neural network (ANN) model based on two-factor of SCD and SD can capture these unique characteristics, with minimal errors and the best performance in simulating SBD and then estimating SWE on the TP. Our method improves the underestimation of the widely used regional average density during the heavy snow period, enables a reasonable conversion of the large amount of historical and growing SDs into critical SWEs, and thus improves our ability to assess snow water resources on the TP.

NH-SWE: Northern Hemisphere Snow Water Equivalent dataset based on in situ snow depth time series

Abstract. Ground-based datasets of observed snow water equivalent (SWE) are scarce, while gridded SWE estimates from remote-sensing and climate reanalysis are unable to resolve the high spatial variability of snow on the ground. Long-term ground observations of snow depth, in combination with models that can accurately convert snow depth to SWE, can fill this observational gap. Here, we provide a new SWE dataset (NH-SWE) that encompasses 11 071 stations in the Northern Hemisphere (NH) and is available at https://doi.org/10.5281/zenodo.7515603 (Fontrodona-Bach et al., 2023). This new dataset provides daily time series of SWE, varying in length between 1 and 73 years, spanning the period 1950–2022, and covering a wide range of snow climates including many mountainous regions. At each station, observed snow depth was converted to SWE using an established snow-depth-to-SWE conversion model, with excellent model performance using regionalised parameters based on climate variables. The accuracy of the model after parameter regionalisation is comparable to that of the calibrated model. The key advantages and strengths of the regionalised model presented here are its transferability across climates and the high performance in modelling daily SWE dynamics in terms of peak SWE, total snowmelt and duration of the melt season, as assessed here against a comparison model. This dataset is particularly useful for studies that require accurate time series of SWE dynamics, timing of snowmelt onset, and snowmelt totals and duration. It can, for example, be used for climate change impact analyses, water resources assessment and management, validation of remote sensing of snow, hydrological modelling, and snow data assimilation into climate models.

Intercomparison of snow water equivalent products in the Sierra Nevada California using airborne snow observatory data and ground observations

Whereas many independent methods are used to estimate snow water equivalent (SWE) and its spatial distribution and seasonal variability, a need exists for a systematic characterization of inter-model differences at annual, seasonal, and regional scales necessary to quantify the associated uncertainty in these datasets. This study conducts a multi-scale validation and comparison, based on Airborne Snow Observatory data, of five state-of-the-art SWE datasets in the Sierra Nevada, California, including three SWE datasets from retrospective models: an INiTial REConstruction model (REC-INT), an improved REConstruction model based on the ParBal energy balance model (REC-ParBal), and a Sierra Nevada SWE REConstruction with Data Assimilation (REC-DA), and two operational SWE datasets from the U.S. National Weather Service, including the Snow Data Assimilation System (SNODAS) and the National Water Model (NWM-SWE). The results show that REC-DA and REC-ParBal provide the two most accurate estimates of SWE in the snowmelt season, both with small positive biases. REC-DA provides the most accurate spatial distribution of SWE (R2 = 0.87, MAE = 66 mm, PBIAS = 8.3%) at the pixel scale, while REC-ParBal has the least basin-wide PBIAS (R2 = 0.79, MAE = 73 mm, PBIAS = 4.1%) in the snowmelt season. Moreover, REC-DA underestimates peak SWE by −5.8%, while REC-ParBal overestimates it by 7.5%, when compared with the measured peak SWE at snow pillow stations across the Sierra Nevada. The two operational SWE products—SNODAS and NWM-SWE—are less accurate. Furthermore, the inter-model comparison reveals a certain amount of disagreement in snow water storage across time and space between SWE datasets. This study advances our understanding of regional SWE uncertainties and provides critical insights to support future applications of these SWE data products and therefore has broad implications for water resources management and hydrological process studies.

Estimation of Snow Water Equivalent in Semiarid Zone from Data of Global Numerical Models ICON and GFS/NCEP (Case Study of the Selenga River Basin)

Estimation of Snow Water Equivalent in Semiarid Zone from Data of Global Numerical Models ICON and GFS/NCEP (Case Study of the Selenga River Basin)

Snow Height and Snow Water Equivalent Estimation from Snow Cover Fraction Using Sentinel-2 Satellite Images in North Kazakhstan

Climate change's influence on snowpack can significantly affect natural and anthropogenic processes. Water resources and agri-business, which depend on winter precipitation, are highly affected by variations in a snowbank and melting regimes. This paper demonstrates the comparison results of the snowpack thickness estimation in the LLP "North Kazakhstan AES" adopting distinctive techniques (quadratic, exponential, and linear functions) for assessing Fractional Snow Cover (SCF) and demonstrating Snow Water Equivalent (SWE) on the one hand, and in-situ perspective on the other. Between the 26–29 of February 2020, a field measurement was managed on the territory of 25,000 hectares. Accordingly, the thickness of the snowpack was surveyed at 560 points, and its density was measured at 70 points. Applying existing methodologies of SCF computation, it became apparent that the quadratic equation provides more reliable results at RMSE of 0.01 m, followed by linear -0.12 m and exponential -0.13 m methods. This work showed a strong correlation between snow height and SCF, namely the quadratic function in Northern Kazakhstan. Thus, we strongly suggest using Sentinel-2 MSI and the quadratic SCF estimation function for snow cover estimation, further spring flood forecasting, and other hydrological studies. Keywords: snow height, snow cover fraction, normalized-difference snow index, snow water equivalent, Sentinel-2 MSI, North Kazakhstan