Reply on CC1

Manual measurement of snow water equivalent (SWE) is still important today for several applications such as hydrological model validation. This measurement can be performed with different types of snow tube sampler or by a snow pit. Although these methods have been performed for several decades, there is an apparent lack of information required to have a consensus regarding the best reference for “true” SWE. We define and estimate the uncertainty and measurement error of different methods of snow pits and snow samplers. Analysis was based upon measurements taken over five consecutive winters (2016–2020) from the same flat and open area. This study compares two snow pit methods and three snow samplers. In addition to including the Standard Federal sampler (SFS), this study documents the first use of two new large diameter samplers, the Hydro-Québec sampler (HQS) and Université Laval sampler (ULS). Large diameter samplers had lowest uncertainty (2.6 to 4.0 %). Snow pit methods had higher uncertainty due to instruments (7.1 to 11.4 %), close to that of the SFS (mean = 10.4 %). Given its larger collected snow volume for estimating SWE and its lower uncertainty, we posit that ULS represents the most appropriate method of reference for “true” SWE. By considering ULS as the reference in calculating mean bias error (MBE), different snow pit methods overestimated SWE by 16.6 to 26.2 %, which was much higher than SFS (8.4 %). This study suggests that large diameter samplers are the best method for estimating “true” SWE.

Manual measurement of snow water equivalent (SWE) is still important today for several applications such as hydrological model validation. This measurement can be performed with different types of snow tube sampler or by a snow pit. Although these methods have been performed for several decades, there is an apparent lack of information required to have a consensus regarding the best reference for “true” SWE. We define and estimate the uncertainty and measurement error of different methods of snow pits and snow samplers. Analysis was based upon measurements taken over five consecutive winters (2016–2020) from the same flat and open area. This study compares two snow pit methods and three snow samplers. In addition to including the Standard Federal sampler (SFS), this study documents the first use of two new large diameter samplers, the Hydro-Québec sampler (HQS) and Université Laval sampler (ULS). Large diameter samplers had lowest uncertainty (2.6 to 4.0 %). Snow pit methods had higher uncertainty due to instruments (7.1 to 11.4 %), close to that of the SFS (mean = 10.4 %). Given its larger collected snow volume for estimating SWE and its lower uncertainty, we posit that ULS represents the most appropriate method of reference for “true” SWE. By considering ULS as the reference in calculating mean bias error (MBE), different snow pit methods overestimated SWE by 16.6 to 26.2 %, which was much higher than SFS (8.4 %). This study suggests that large diameter samplers are the best method for estimating “true” SWE.

Area-based measurements of snow water equivalent (SWE) are important for understanding earth system processes such as glacier mass balance, winter hydrological storage in drainage basins and ground thermal regimes. Remote sensing techniques are ideally suited for wide-scale area-based mapping with the most commonly used technique to measure SWE being passive-microwave, which is limited to coarse spatial resolutions of 25 km or greater, and to areas without significant topographic variation. Passive-microwave also has a negative bias for large SWE. Repeat-pass synthetic aperture radar interferometry (InSAR) as an alternate technique allows measurement of SWE change at much higher spatial resolution. However, it has not been widely adopted because: (1) the phase unwrapping problem has not been robustly addressed, especially for interferograms with poor coherence and; (2) SWE change maps scaled directly from repeat-pass interferograms are not an absolute measurement but contain unknown offsets for each contiguous coherent area. We develop and test a novel method for repeat-pass InSAR based dry-snow SWE estimation that exploits the sensitivity of the dry-snow refraction-induced InSAR phase to topographic variations. The method robustly estimates absolute SWE change at spatial resolutions of < 1 km, without the need for phase unwrapping. We derive a quantitative signal model for this new SWE change estimator and identify the relevant sources of bias. The method is demonstrated using both simulated SWE distributions and a 9-year RADARSAT-2 spotlight-mode dataset near Inuvik, NWT, Canada. SWE results are compared to in situ snow survey measurements and estimates from ERA5 reanalysis. Our method performs well in high-relief areas and in areas with high SWE (> 150 mm), thus providing complementary coverage to other passive- and active-microwave based SWE estimation methods. Further, our method has the advantage of requiring only a single wavelength band and thus can utilize existing spaceborne synthetic aperture radar systems. In application, a first order analysis of SWE trends within three drainage basins suggests that differences between basin-level accumulations are a function of major landcover types, and that re-vegetation following a forest-tundra fire that occurred over 50 years ago continues to affect the spatial distribution of SWE accumulation in the study area.

Area-based measurements of snow water equivalent (SWE) are important for understanding earth system processes such as glacier mass balance, winter hydrological storage in drainage basins and ground thermal regimes. Remote sensing techniques are ideally suited for wide-scale area-based mapping with the most commonly used technique to measure SWE being passive-microwave, which is limited to coarse spatial resolutions of 25 km or greater, and to areas without significant topographic variation. Passive-microwave also has a negative bias for large SWE. Repeat-pass synthetic aperture radar interferometry (InSAR) as an alternate technique allows measurement of SWE change at much higher spatial resolution. However, it has not been widely adopted because: (1) the phase unwrapping problem has not been robustly addressed, especially for interferograms with poor coherence and; (2) SWE change maps scaled directly from repeat-pass interferograms are not an absolute measurement but contain unknown offsets for each contiguous coherent area. We develop and test a novel method for repeat-pass InSAR based dry-snow SWE estimation that exploits the sensitivity of the dry-snow refraction-induced InSAR phase to topographic variations. The method robustly estimates absolute SWE change at spatial resolutions of < 1 km, without the need for phase unwrapping. We derive a quantitative signal model for this new SWE change estimator and identify the relevant sources of bias. The method is demonstrated using both simulated SWE distributions and a 9-year RADARSAT-2 spotlight-mode dataset near Inuvik, NWT, Canada. SWE results are compared to in situ snow survey measurements and estimates from ERA5 reanalysis. Our method performs well in high-relief areas and in areas with high SWE (> 150 mm), thus providing complementary coverage to other passive- and active-microwave based SWE estimation methods. Further, our method has the advantage of requiring only a single wavelength band and thus can utilize existing spaceborne synthetic aperture radar systems. In application, a first order analysis of SWE trends within three drainage basins suggests that differences between basin-level accumulations are a function of major landcover types, and that re-vegetation following a forest-tundra fire that occurred over 50 years ago continues to affect the spatial distribution of SWE accumulation in the study area.

Abstract. During the World Meteorological Organization (WMO) Solid Precipitation Intercomparison Experiment (SPICE), automated measurements of snow water equivalent (SWE) were made at the Sodankylä (Finland), Weissfluhjoch (Switzerland) and Caribou Creek (Canada) SPICE sites during the northern hemispheric winters of 2013/14 and 2014/15. Supplementary intercomparison measurements were made at Fortress Mountain (Kananaskis, Canada) during the 2013/14 winter. The objectives of this analysis are to compare automated SWE measurements with a reference, comment on their performance and, where possible, to make recommendations on how to best use the instruments and interpret their measurements. Sodankylä, Caribou Creek and Fortress Mountain hosted a Campbell Scientific CS725 passive gamma radiation SWE sensor. Sodankylä and Weissfluhjoch hosted a Sommer Messtechnik SSG1000 snow scale. The CS725 operating principle is based on measuring the attenuation of soil emitted gamma radiation by the snowpack and relating the attenuation to SWE. The SSG1000 measures the mass of the overlying snowpack directly by using a weighing platform and load cell. Manual SWE measurements were obtained at the intercomparison sites on a bi-weekly basis over the accumulation–ablation periods using bulk density samplers. These manual measurements are considered to be the reference for the intercomparison. Results from Sodankylä and Caribou Creek showed that the CS725 generally overestimates SWE as compared to manual measurements by roughly 30–35 % with correlations (r2) as high as 0.99 for Sodankylä and 0.90 for Caribou Creek. The RMSE varied from 30 to 43 mm water equivalent (mm w.e.) and from 18 to 25 mm w.e. at Sodankylä and Caribou Creek, which had respective SWE maximums of approximately 200 and 120 mm w.e. The correlation at Fortress Mountain was 0.94 (RMSE of 48 mm w.e. with a maximum SWE of approximately 650 mm w.e.) with no systematic overestimation. The SSG1000 snow scale, having a different measurement principle, agreed quite closely with the manual measurements at Sodankylä and Weissfluhjoch throughout the periods when data were available (r2 as high as 0.99 and RMSE from 8 to 24 mm w.e. at Sodankylä and from 56 to 59 mm w.e. at Weissfluhjoch, where maximum SWE was approximately 850 mm w.e.). When the SSG1000 was compared to the CS725 at Sodankylä, the agreement was linear until the start of ablation when the positive bias in the CS725 increases substantially relative to the SSG1000. Since both Caribou Creek and Sodankylä have sandy soil, water from the snowpack readily infiltrates into the soil during melt, even if the soil is frozen. However, the CS725 does not differentiate this water from the unmelted snow. This issue can be identified, at least during the late spring ablation period, with soil moisture and temperature observations like those measured at Caribou Creek. With a less permeable soil and surface runoff, the increase in the instrument bias during ablation is not as significant, as shown by the Fortress Mountain intercomparison.

Comparison of different automatic methods for estimating snow water equivalent

Direct measurement of snow water equivalent (SWE) in snow-dominated mountainous areas is difficult, thus its prediction is essential for water resources management in such areas. In addition, because of nonlinear trend of snow spatial distribution and the multiple influencing factors concerning the SWE spatial distribution, statistical models are not usually able to present acceptable results. Therefore, applicable methods that are able to predict nonlinear trends are necessary. In this research, to predict SWE, the Sohrevard Watershed located in northwest of Iran was selected as the case study. Database was collected, and the required maps were derived. Snow depth (SD) at 150 points with two sampling patterns including systematic random sampling and Latin hypercube sampling (LHS), and snow density at 18 points were randomly measured, and then SWE was calculated. SWE was predicted using artificial neural network (ANN), adaptive neuro-fuzzy inference system (ANFIS) and regression methods. The results showed that the performance of ANN and ANFIS models with two sampling patterns were observed better than the regression method. Moreover, based on most of the efficiency criteria, the efficiency of ANN, ANFIS and regression methods under LHS pattern were observed higher than the systematic random sampling pattern. However, there were no significant differences between the two methods of ANN and ANFIS in SWE prediction. Data of both two sampling patterns had the highest sensitivity to the elevation. In addition, the LHS and the systematic random sampling patterns had the least sensitivity to the profile curvature and plan curvature, respectively.

The assessment of snow water equivalent (SWE) is crucial for hydrological studies in glaciated catchments to quantify accumulation and ablation of a seasonal snow cover for both, the glaciated and non-glaciated terrain. Presently, the majority of the SWE assessment on glaciated terrain relies on manual measurements once or a few times per year, given the limited techniques for continuous SWE monitoring and the challenging conditions in a high mountain environment. Cosmic Ray Neutron Sensors (CRNS) offer to overcome these limitations providing sub-daily SWE estimates derived from neutron counts.This study employs a CRNS installed on an Alpine glacier (Hintereisferner, Austria) over three years, complemented with an additional CRNS for one winter roughly 300 m lower in elevation along the glacier&amp;#8217;s central flow line. Comparing CRNS outputs with frequent manual SWE measurements, the results demonstrate close agreement in SWE and snow density. The CRNS were found to be remarkably resilient in harsh conditions, providing nearly continuous hourly data over the measurement period. Additionally, the study evaluates at the CRNS locations the performance of two snow models, which might be considered as end members of model complexity &amp;#8211; SNOWPACK and &amp;#916;SNOW in its latest version. While SNOWPACK, with its physically-based approach, yields the best results, &amp;#916;SNOW stands out for its simplicity, requiring only daily snow depth observations as input and performs almost as well as SNOWPACK in terms of mean absolute SWE error.&amp;#160;Comparing the SWE measurements with winter precipitation from weighing gauges distributed in the catchment gives interesting details of precipitation gradients with elevation. Precipitation gradients interpolated from non-glaciated to glaciated terrain are considerably higher than on non-glaciated terrain only. On the latter, empirical bulk correction factors are frequently applied, which might fail on glaciers due to their different topographic setting. This highlights the need of separate treatment of snow on glaciers in hydrological models for correct SWE representation across the catchment.&amp;#160;

Peak snowpack in the western USA has decreased 21% since the mid‐20th century, and these trends are forecasted to continue over the next century. In water‐limited systems, forest productivity during the growing season is assumed to be linked with snowpack during the previous winter, but this linkage has proven not to be universal in studies across varying ecosystems. We compared peak Normalized Difference Vegetation Index (NDVI, a proxy for vegetation productivity) at relatively fine scales (30‐m pixel) to point measurements of snow water equivalent (SWE) at 169 locations scattered among the northern Intermountain West for the 25‐yr period between 1993 and 2017. We further compared NDVI to additional environmental variables including topographical variables (elevation, aspect, and slope), and climatic variables (Palmer Drought Severity Index, snowmelt date, and snowmelt length). We hypothesized that peak SWE would be a strong predictor of peak NDVI, but the strength of the SWE–NDVI relationships would vary spatially in our topographically complex study area. We observed weak relationships between snowpack and vegetation productivity, with only 10% of the locations displaying statistically significant correlations between annual peak SWE and peak NDVI. Furthermore, we determined that the response of NDVI to snowpack and environmental variables differed among subregions of our study area. For example, we observed strong SWE–NDVI associations in the continental‐most portion of our study area and weak correlations in other subregions of our study area. Our results highlight that in regions characterized by complex topography, (1) environmental drivers of NDVI are highly variable and (2) fine‐resolution studies more accurately capture the spatial heterogeneity of NDVI compared to course‐resolution studies. However, one consequence of fine‐resolution studies may be that microclimate effects at each site swamp the signal of predictor‐ecological response relationships, whereas these relationships may be more discernable in larger scale remote sensing studies.

Development of the snow- and ice-accounting routine (SIAR)

Abstract. Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across nine trenches) collected over two winters at Trail Valley Creek, NWT, Canada, was applied in synthetic radiative transfer experiments. This allowed for robust assessment of the impact of estimation accuracy of unknown snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability in total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths of layer thickness were less than a metre for all layers. Depth hoar was consistently ∼30 % of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and specific surface area (SSA) of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7 % under the median value was needed to accurately retrieve SWE. In shallow snowpacks &lt;0.6 m, depth hoar SSA estimates of ±5 %–10 % around the optimal retrieval SSA allowed SWE retrievals within a tolerance of ±30 mm. Where snowpacks were deeper than ∼30 cm, accurate values of representative SSA for depth hoar became critical as retrieval errors were exceeded if the median depth hoar SSA was applied.

Abstract. We analyze snow water equivalent (SWE) measurements from a cosmic ray sensor (CRS) on the lower accumulation area of Mera Glacier (central Himalaya, Nepal) between November 2019 and November 2021. The CRS aligned well with field observations and revealed accumulation in pre-monsoon and monsoon observations, followed by ablation in post-monsoon and winter observations. COSIPY simulations suggest significant surface melting, water percolation, and refreezing within the snowpack, consistent with CRS observations yet liable to be missed by surface mass balance surveys. We conclude that CRS can be used to determine mass fluxes in various climatic settings, but the interpretation of the total changes in SWE needs complementary measurements and model analysis to determine the share of specific mass fluxes, such as melt and refreezing.

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra‐taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.

The current study evaluates the impacts of various sources of uncertainty involved in hydrologic modeling on parameter behavior and regionalization utilizing different Bayesian likelihood functions and the Differential Evolution Adaptive Metropolis (DREAM) algorithm. The developed likelihood functions differ in their underlying assumptions and treatment of error sources. We apply the developed method to a snow accumulation and ablation model (National Weather Service SNOW17) and generate parameter ensembles to predict snow water equivalent (SWE). Observational data include precipitation and air temperature forcing along with SWE measurements from 24 sites with diverse hydroclimatic characteristics. A multiple linear regression model is used to construct regionalization relationships between model parameters and site characteristics. Results indicate that model structural uncertainty has the largest influence on SNOW17 parameter behavior. Precipitation uncertainty is the second largest source of uncertainty, showing greater impact at wetter sites. Measurement uncertainty in SWE tends to have little impact on the final model parameters and resulting SWE predictions. Considering all sources of uncertainty, parameters related to air temperature and snowfall fraction exhibit the strongest correlations to site characteristics. Parameters related to the length of the melting period also show high correlation to site characteristics. Finally, model structural uncertainty and precipitation uncertainty dramatically alter parameter regionalization relationships in comparison to cases where only uncertainty in model parameters or output measurements is considered. Our results demonstrate that accurate treatment of forcing, parameter, model structural, and calibration data errors is critical for deriving robust regionalization relationships.

In this report we summarize the outcome of our research project Near real-time observations of snow water equivalent for SIOS on Svalbard – (SWESOS), which was funded by Svalbard Integrated Arctic Earth Observing System (SIOS). We evaluate a new automated monitoring technique for measuring snow water equivalent (SWE) using a passive gamma ray sensor at the Bayelva site on Svalbard. Our goal is to provide high quality, continuous, near-real time data for the observation of SWE. From the new automated measurements we derive the essential climate variables of snow cover (snow height, SWE, spatial coverage) and dynamics (timing of onset and melt, wetting) and provide a unique dataset by linking existing snow datasets with the new SWE measurements. This dataset will be essential to develop and calibrate snow, permafrost and hydrologic models in this data scarce region.