The transition of a cold winter snowpack to one that is ripe and contributing to runoff is crucial to gauge for water resource management, but is highly variable in space and time. Snow surface melt/freeze cycles, associated with diurnal fluctuations in radiative inputs, are hallmarks of this transition. C-band synthetic aperture radar (SAR) reliably detects meltwater in the snowpack. Sentinel-1 (S1) C-band SAR offers consistent acquisition patterns that allow for diurnal investigations of melting snow. We used over 50 snow pit observations from 2020 in Grand Mesa, Colorado, USA, to track temperature and wetness in the snowpack as a function of depth and time during snowpack phases of warming, ripening, and runoff. We also ran the physically-based SnowModel, which provided a spatially and temporally continuous independent indication of snowpack conditions. Snowpack phases were identified and corroborated by comparing field measurements with SnowModel outputs. Knowledge of snowpack warming, ripening, and runoff phases was used to interpret diurnal changes in S1 backscatter values. Both field measurements and SnowModel simulations suggested that S1 SAR was not sensitive to the initial snowpack warming phase on Grand Mesa. In the ripening and runoff phases, the diurnal cycle in S1 SAR co-polarized backscatter was affected by both surface melt/freeze as well as the conditions of the snowpack underneath (ripening or ripe). The ripening phase was associated with significant increases in morning backscatter values, likely due to volume scattering from surface melt/freeze crusts, as well as significant decreases in evening backscatter values associated with snowmelt. During the runoff phase, both morning and evening backscatter decreased compared to reference values. These unique S1 diurnal signatures, and their interpretations using field measurements and SnowModel outputs, highlight the capacities and limitations of S1 SAR to understand snow surface states and bulk phases, which may offer runoff forecasting or energy balance model validation or parameterization, especially useful in remote or sparsely-gauged alpine basins.

Remote sensing of snow-covered area (SCA) can be binary (indicating the presence/absence of snow cover at each pixel) or fractional (indicating the fraction of each pixel covered by snow). Fractional SCA mapping provides more information than binary SCA, but is more difficult to implement and may not be feasible with all types of remote sensing data. The utility of fractional SCA mapping relative to binary SCA mapping varies with the intended application as well as by spatial resolution, temporal resolution and period of interest, and climate. We quantified the frequency of occurrence of partially snow-covered (mixed) pixels at spatial resolutions between 1 m and 500 m over five dates at two study areas in the western U.S., using 0.5 m binary SCA maps derived from high spatial resolution imagery aggregated to fractional SCA at coarser spatial resolutions. In addition, we used in situ monitoring to estimate the frequency of partially snow-covered conditions for the period September 2013–August 2014 at 10 60-m grid cell footprints at two study areas with continental snow climates. Results from the image analysis indicate that at 40 m, slightly above the nominal spatial resolution of Landsat, mixed pixels accounted for 25%–93% of total pixels, while at 500 m, the nominal spatial resolution of MODIS bands used for snow cover mapping, mixed pixels accounted for 67%–100% of total pixels. Mixed pixels occurred more commonly at the continental snow climate site than at the maritime snow climate site. The in situ data indicate that some snow cover was present between 186 and 303 days, and partial snow cover conditions occurred on 10%–98% of days with snow cover. Four sites remained partially snow-free throughout most of the winter and spring, while six sites were entirely snow covered throughout most or all of the winter and spring. Within 60 m grid cells, the late spring/summer transition from snow-covered to snow-free conditions lasted 17–56 days and averaged 37 days. Our results suggest that mixed snow-covered snow-free pixels are common at the spatial resolutions imaged by both the Landsat and MODIS sensors. This highlights the additional information available from fractional SCA products and suggests fractional SCA can provide a major advantage for hydrological and climatological monitoring and modeling, particularly when accurate representation of the spatial distribution of snow cover is critical.

We developed an automated approach for mapping persistent ice and snow cover (glaciers and perennial snowfields) from Landsat TM and ETM+ data across a variety of topography, glacier types, and climatic conditions at high latitudes (above ~65°N). Our approach exploits all available Landsat scenes acquired during the late summer (1 August–15 September) over a multi-year period and employs an automated cloud masking algorithm optimized for snow and ice covered mountainous environments. Pixels from individual Landsat scenes were classified as snow/ice covered or snow/ice free based on the Normalized Difference Snow Index (NDSI), and pixels consistently identified as snow/ice covered over a five-year period were classified as persistent ice and snow cover. The same NDSI and ratio of snow/ice-covered days to total days thresholds applied consistently across eight study regions resulted in persistent ice and snow cover maps that agreed closely in most areas with glacier area mapped for the Randolph Glacier Inventory (RGI), with a mean accuracy (agreement with the RGI) of 0.96, a mean precision (user’s accuracy of the snow/ice cover class) of 0.92, a mean recall (producer’s accuracy of the snow/ice cover class) of 0.86, and a mean F-score (a measure that considers both precision and recall) of 0.88. We also compared results from our approach to glacier area mapped from high spatial resolution imagery at four study regions and found similar results. Accuracy was lowest in regions with substantial areas of debris-covered glacier ice, suggesting that manual editing would still be required in these regions to achieve reasonable results. The similarity of our results to those from the RGI as well as glacier area mapped from high spatial resolution imagery suggests it should be possible to apply this approach across large regions to produce updated 30-m resolution maps of persistent ice and snow cover. In the short term, automated PISC maps can be used to rapidly identify areas where substantial changes in glacier area have occurred since the most recent conventional glacier inventories, highlighting areas where updated inventories are most urgently needed. From a longer term perspective, the automated production of PISC maps represents an important step toward fully automated glacier extent monitoring using Landsat or similar sensors.

The surface mass balance (SMB) of West Antarctica is an important glaciological input to understanding polar climate and sea-level rise but with historically poor in situ data coverage. Previous studies demonstrate the utility of frequency-modulated continuous-wave radar to image subsurface layering in ice sheets, providing an additional source of data with which to estimate SMB. Traditional methods, however, require time-intensive manual oversight. Here, we present a probabilistic, fully automated approach to estimate annual SMB and uncertainties from radar echograms using successive peak-finding and weighted neighborhood search algorithms with the Monte Carlo simulations based on annual-layer likelihood scores. We apply this method to ground-based and airborne radar in a 175-km transect of the West Antarctic interior and compare the results to traditional manual methods and independent estimates from firn cores. The method demonstrates an automated estimation of SMB across a range of accumulation rates (100-450-mm water equivalent per year) and layer gradients up to 2 m/km. Based on likelihood-weighted F-scores, automated layer picks have a success rate between 64% and 84.6% compared with manually picked layers for three validation sites dispersed across the region. Comparisons between the automated SMB estimates and independent firn cores show a bias of 24 ± 70-mm water equivalent per year (12% ± 35% water equivalent of the in situ mean accumulation rate) although individual core site biases differ. This new approach permits the fully automated extraction of annual SMB rates and should be broadly and readily applicable to previously collected and ongoing radar data sets across polar regions.

The Bonneville Salt Flats (BSF) of Utah, with a notable history in land-speed racing and resource extraction, presents a dynamic landscape shaped by both natural processes and anthropogenic activities. This study leverages Interferometric Synthetic Aperture Radar (InSAR) techniques to monitor the BSF, focusing on the potential of coherence as an indicator of surface changes. A total of 165 Sentinel-1 A&B coherence images, with perpendicular baselines under 100 meters and temporal baselines of 12, 24, and 36 days, were analyzed using spatially and temporally reduced statistics, along with median absolute deviation variance. The findings demonstrate consistently high coherence (>0.90) during dry conditions, contrasted by significantly lower coherence (<0.40) during periods of precipitation and flooding. Coherence showed moderate reduction (<0.70) in zones of recreational activities (driving & racing), highlighting the relative impact of land-use activities on the salt crust. Coherence observations indicate minimal difference during periods overlapping with and following official racing events compared to periods of racetrack preparation and general recreational use. Zones with the highest variability are linked to disturbances from recreation and resource extraction, while areas with minimal variability correspond to regions with thicker halite crust. The data also revealed that large precipitation events significantly influence coherence, with a clear power-law decrease in mean coherence as precipitation and floodwater accumulation increased. This research not only validates the use of coherence for monitoring environmental evolution and the influences of anthropogenic activities at the BSF but also supports the potential of future InSAR applications for accurate phase assessments during dry periods.