# Snow depth variability in the Northern Hemisphere mountains observed from space

Hans Lievens, Matthias Demuzere, Hans-Peter Marshall, Rolf Reichle, Ludovic Brucker, Isis Brangers, Patricia de Rosnay, Marie Dumont, Manuela Girotto, Walter Immerzeel, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Hans Lievens, Matthias Demuzere, Hans-Peter Marshall, Rolf Reichle, Ludovic Brucker, et al. Snow depth variability in the Northern Hemisphere mountains observed from space. Nature Communications, 2019, 10 (1), pp.4629. 10.1038/s41467-019-12566-y. meteo-03657921

HAL Id: meteo-03657921

https://meteofrance.hal.science/meteo-03657921

Submitted on 3 May 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Snow depth variability in the Northern Hemisphere mountains observed from space

Hans Lievens<fnref n="1" />,²*, Matthias Demuzere<fnref n="1" />,²,³, Hans-Peter Marshall⁴,⁵, Rolf H. Reichle⁶, Ludovic Brucker⁶,⁷, Isis Brangers<fnref n="1" />, Patricia de Rosnay⁸, Marie Dumont⁹, Manuela Girotto⁶,⁷,¹⁰, Walter W. Immerzeel¹¹, Tobias Jonas¹², Edward J. Kim⁶, Inka Koch¹³, Christoph Marty<fnref n="1" />⁲, Tuomo Saloranta¹⁴, Johannes Schöber¹⁵ & Gabrielle J.M. De Lannoy<fnref n="1" />

Accurate snow depth observations are critical to assess water resources. More than a billion people rely on water from snow, most of which originates in the Northern Hemisphere mountain ranges. Yet, remote sensing observations of mountain snow depth are still lacking at the large scale. Here, we show the ability of Sentinel-1 to map snow depth in the Northern Hemisphere mountains at 1 km² resolution using an empirical change detection approach. An evaluation with measurements from ~4000 sites and reanalysis data demonstrates that the Sentinel-1 retrievals capture the spatial variability between and within mountain ranges, as well as their inter-annual differences. This is showcased with the contrasting snow depths between 2017 and 2018 in the US Sierra Nevada and European Alps. With Sentinel-1 continuity ensured until 2030 and likely beyond, these findings lay a foundation for quantifying the long-term vulnerability of mountain snow-water resources to climate change.

<footnote n="1">¹Department of Earth and Environmental Sciences, KU Leuven, Heverlee, Belgium. ²Laboratory of Hydrology and Water Management, Ghent University, Ghent, Belgium. ³Department of Geography, Ruhr-University Bochum, Bochum, Germany. ⁴Department of Geosciences, Boise State University, Boise, ID, USA. ⁵Remote Sensing and GIS Center, U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, NH, USA. ⁶NASA Goddard Space Flight Center, Greenbelt, MD, USA. ⁷GESTAR, Universities Space Research Association, Columbia, MD, USA. ⁸European Centre for Medium-Range Weather Forecasts, Reading, UK. ⁹Université Grenoble Alpes, Université de Toulouse, Météo-France, Grenoble, France, CNRS, CNRM, Centre d'Etudes de la Neige, Grenoble, France. ¹⁰Environmental Science and Policy Management Department, University of California, Berkeley, CA, USA. ¹¹Department of Physical Geography, Utrecht University, Utrecht, The Netherlands. ¹²WSL Institute for Snow and Avalanche Research SLF, Davos, Switzerland. ¹³International Centre for Integrated Mountain Development, Kathmandu, Nepal. ¹⁴Hydrology Department, Norwegian Water Resources and Energy Directorate NVE, Oslo, Norway. ¹⁵TIWAG, Tiroler Wasserkraft AG, Innsbruck, Austria. *email: hans.lievens@kuleuven.be</footnote>
Snow has a large-scale cooling effect on our planet by reflecting most of the incoming solar radiation and by dissipating energy during seasonal melt¹,². Water resources from snowmelt, of which the majority originates in mountain ranges, are indispensable. For instance, they provide drinking water to over a billion people³,⁴, supply 3/4th of the crop production in the western US⁵, generate hydro-electric power worldwide⁴, and sustain the vast urban settlements in Himalayan watersheds⁶⁻⁸. Despite this importance, we lack quantitative estimates of how much snow is stored in each mountain range on Earth⁹. The current estimates from the interpolation of local measurements are unrealistic where measurements are sparse⁹, whereas estimates from numerical weather prediction systems are poor due to the large uncertainty in mountain snowfall⁴.

Remote sensing provides essential snow depth information, yet currently operational observations have shortcomings⁹. Passive microwave observations¹⁰,¹¹ typically exclude mountains because their coarse (~25 km) footprints cannot resolve the spatial variability, and the observations saturate in deep snow (>0.8 m depth)¹²,¹³. Airborne lidar systems⁵,¹⁴ are likely the most accurate to date, but are limited to targeted mountain areas and favorable weather conditions. Some alternative methods have shown promise at the local scale. These include structure-from-motion¹⁵,¹⁶, constructing snow depth from overlapping image pairs, X-/K-band (8–18 GHz) Synthetic Aperture Radar (SAR)¹⁷⁻¹⁹, retrieving snow properties based on the backscattered energy from an illuminated snowpack, and SAR interferometry²⁰,²¹, measuring changes in the radar signal phase from refraction in the snow. A community effort to compare some of the above-mentioned and more experimental remote sensing techniques is currently ongoing within the NASA SnowEx activity that should ultimately lead to a space mission designed for snow²². New and robust satellite observations are critically needed to fill the mountain-snow observation gap⁴.

In this study, we demonstrate the ability of the ESA and Copernicus Sentinel-1 constellation (two satellites: 1A and 1B) to map snow depth across the Northern Hemisphere mountains. Presently, Sentinel-1 is the only SAR mission providing high-resolution backscatter measurements (at C-band; 5.4 GHz) with a revisit time of 6 days suitable for snow monitoring. Given the strong absorption of C-band microwave radiation by wet snow, monitoring applications have so far focused on the mapping of melt²³. The use of C-band backscatter for estimating snow depth (or mass, related to depth by the density) has long been swept aside after early satellite measurements had shown a limited sensitivity²⁴,²⁵. However, these studies were mostly surveying shallow snow outside mountain environments and, more importantly, were limited to backscatter measurements in copolarization. Cross-polarized backscatter measurements were to-date only investigated at the local scale using tower installations, with strongly contradicting results²⁶,²⁷. Here, we demonstrate the value of including cross-polarized backscatter measurements from C-band satellite to retrieve snow depth in mountainous areas at the large scale.

The Sentinel-1 snow depth retrievals over the Northern Hemisphere mountain ranges allow for near-real-time monitoring that complements current snow tracking systems from the World Meteorological Organization (WMO). Unlike currently available space-borne passive microwave measurements, Sentinel-1 measurements are suitable for retrieving snow depth in mountainous environments owing to their high resolution and the herein demonstrated sensitivity to deep snow. We provide clear evidence of the value of the empirical Sentinel-1 retrievals through a comparison against point-scale snow depth measurements at ~4000 sites and coarse-scale global reanalysis data. Over the US Sierra Nevada and the European Alps, snow depths of two

consecutive winters are stratified by elevation to highlight inter-annual differences. While resembling the large-scale patterns of the reanalysis data, Sentinel-1 more accurately reveals the spatial detail and the elevation profile of snow depth. Finally, we estimate total snow volumes for the top 100 snowiest mountains in the Northern Hemisphere to support their water resources monitoring.

## Results

**Sentinel-1 backscatter signatures over snow.** The Sentinel-1 constellation routinely illuminates the land surface with C-band radiation and measures the backscatter ($\sigma^0$) in co-polarization and cross-polarization, i.e., vertical-vertical (vv) and vertical-horizontal (vh) transmit-receive, respectively. The total backscatter ($\sigma_{pq}^{\text{tot}}$) for transmit-receive polarizations p and q, originating from a snowpack without vegetation, can be approximated by a four-component model²⁸:

$$ \sigma_{pq}^{\text{tot}} = \sigma_{pq}^{\text{air-snow}} + \sigma_{pq}^{\text{snowvol}} + \sigma_{pq}^{\text{snowvol-grnd}} + e^{(-2\tau_p / \cos \theta)} \cdot \sigma_{pq}^{\text{grnd}} \quad (1) $$

with $\sigma_{pq}^{\text{air-snow}}$ the backscatter from the air-snow interface, $\sigma_{pq}^{\text{snowvol}}$ the backscatter from the snow volume, $\sigma_{pq}^{\text{snowvol-grnd}}$ the higher-order interactions between the snow volume and the ground, and $\sigma_{pq}^{\text{grnd}}$ the backscatter from the ground, attenuated by the snow-pack through $e^{(-2\tau_p / \cos \theta)}$. The amount of attenuation depends on the optical thickness of the snow ($\tau_p$) and the radar incidence angle ($\theta$). Commonly, $\sigma_{pq}^{\text{air-snow}}$ and $\sigma_{pq}^{\text{snowvol-grnd}}$ are neglected²⁹ and thus not further considered here.

Sentinel-1 $\sigma^0$ measurements are processed at 1 km² resolution over the Northern Hemisphere mountains³⁰ (Fig. 1a), using standard processing techniques (Fig. 2 and Methods), and their temporal signatures are compared to those of in situ snow depth measurements. Figure 3 illustrates the temporal evolution in $\sigma^0$ and snow depth measurements at 4 representative sites, along with corresponding snow presence observations at 1 km² resolution from the Interactive Multi-sensor Snow and Ice Mapping System (IMS)³¹. The co-polarized $\sigma_{vv}^0$ measurements (Fig. 3, left column) show little variation throughout the winter, due to the limited absorption or scattering by dry snow in vv-polarization at C-band. In contrast, a sharp (~5 dB) decrease in $\sigma_{vv}^0$ is observed with the wetting and depletion of the snow through spring. This confirms previous studies that demonstrated the limited sensitivity of C-band co-polarized satellite measurements to changes in dry snow and the large absorption and reflection of the signal by wet snow²³⁻²⁵. It also confirms previous findings that C-band vv-polarized measurements are suitable for snowmelt mapping²³. The increase in $\sigma_{vv}^0$ towards the end of the snowmelt period in late spring is caused by the decreasing surface area covered by wet snow that contributes to absorption³². Simultaneously, surface thaw and vegetation green-up in snow-free areas can increase scattering³². The evolution of $\sigma_{vv}^0$ during the snow season is dominated by the fourth term in Eq. (1): the backscatter from the ground surface, which is attenuated by the snowpack only when the snow is wet. The second term, representing the scattering in vv-polarization within a dry snow volume, is generally too small to have a noticeable impact at C-band.

The Sentinel-1 cross-polarized $\sigma_{vh}^0$ measurements (Fig. 3, middle column) gradually increase with the accumulation of (dry) snow during the winter. This suggests an increasing depolarization of the incoming v-polarized signal through anisotropic or multiple scattering on ice crystals, bonds or clusters of ice crystals, or inhomogeneities within the snow volume¹⁷,¹⁹,²⁶,³³⁻³⁵. As the thickness of the snowpack increases,
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Study domain and data. **a** The study domain corresponds to the Northern Hemisphere mountains, north of 20°N. Colors indicate the number (N; dimensionless) of Sentinel-1 backscatter measurements for the period September 2016 through August 2018. N is calculated after combining the ascending and descending Sentinel-1A and Sentinel-1B measurements at the daily timescale, and is limited to mountain areas observed by Sentinel-1 in vv-polarization and vh-polarization. Measurements are on average available every ~2 days in Europe and every ~4 days to ~2 weeks in North America and Asia. **b** In situ snow depth measurement sites within the study domain. The sources of the in situ snow depth measurements are given in Table 1. The majority of the 4175 sites is located in the western US and Europe. Measurements in Asia are more scarce. Colors highlight 11 focus mountain ranges, which were selected based on the availability of in situ measurements and based on differences in geographic location, climate, distance to the coast, elevation, and ruggedness, i.e., representing a broad variety of snow conditions</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Flowchart of the Sentinel-1 processing. **a** The standard backscatter processing techniques are performed using the Google Earth Engine Python's api. **b** Further processing steps, and **c** the application of the retrieval algorithm are performed offline</figcaption></figure>

the path length of the radar signal through the snow increases, increasing the opportunities for scattering¹⁹. Note that this anisotropic and multiple scattering in theory also impacts the backscatter in co-polarization ($\sigma_{vv}^0$), but with a negligible

contribution compared to that of the ground surface. During snow ablation, $\sigma_{vh}^0$ decreases considerably, which we hypothesize is caused by the absorption and reflection of the signal by wet snow, and a decreasing amount of snow volume scattering in a shallowing snowpack. To explain the variations in $\sigma_{vh}^0$, both the second and fourth terms in Eq. (1) are needed. The snow volume scattering component in vh-polarization is no longer negligible compared to the contribution of the ground surface.

The cross-polarization ratio $\sigma_{vh}^0/\sigma_{vv}^0$ (in linear scale, converted to dB) shows overall a stronger correlation with snow depth than $\sigma_{vh}^0$ (Fig. 3, middle and right columns). Taking the ratio may partially eliminate the effects of temporal changes in the ground surface, vegetation, or snow conditions, which similarly impact both co- and cross-polarization. In autumn and winter, $\sigma_{vh}^0/\sigma_{vv}^0$ increases due to the increased amount of volume scattering $\sigma_{vh}^0$ for approximately constant surface scattering $\sigma_{vv}^0$. During spring melt, the ratio decreases because of the relatively higher decrease of $\sigma_{vh}^0$ compared to that of $\sigma_{vv}^0$. Changes in snow properties between dry and wet conditions, such as in snow microstructure and liquid water content, potentially modify the proportionality of volume versus surface scattering, and therefore the sensitivity of the ratio to snow depth. Thereby, the sensitivity to snow depth for wet snow is much more uncertain. Strong melt events, with associated high liquid water contents, may cause fluctuations in $\sigma_{vh}^0/\sigma_{vv}^0$ as the signal is more strongly reflected and absorbed by wet snow layers²⁷. Although this effect may be less severe at C-band than at higher frequencies and may occasionally be alleviated by the refreezing of snow²⁷ prior to the early-morning (6 a.m.) and evening (6 p.m.) overpass times of Sentinel-1, it is likely to have a
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Backscatter and snow depth time series at 4 sites. The left column shows the co-polarized backscatter ($\sigma_{vv}^0$; in dB), the middle column the cross-polarized backscatter ($\sigma_{vh}^0$; in dB), and the right column the cross-polarization ratio ($\sigma_{vh}^0/\sigma_{vv}^0$; in dB). Snow cover presence (dimensionless) from IMS is shown in gray. Site details, including mountain range and latitude/longitude/elevation: a Omineca Mountains, British Columbia, Canada (57.53°N/–126.73°E/1331 m), b Absaroka Range, Wyoming, USA (45.22°N/–110.24°E/2697 m), c San Juan Mountains, Colorado, USA (37.71°N/–107.51°E/3535 m), d European Alps, Austria (46.87°N/10.71°E/2464 m). The ranges of the backscatter (dB) and snow depth (m) measurements have been adjusted between sites to improve visualization</figcaption></figure>

confounding impact on the C-band sensitivity to snow depth. Similarly, the backscatter signals are impacted by successive melt-refreeze cycles that can modify the microstructure and stratigraphy of the snowpack. The primary objective of this study is to map the snow depth for dry snow conditions. However, we do not exclude wet snow conditions from the analysis and assess the performance of the retrievals throughout the snow season.

**Sentinel-1 snow depth retrievals.** The retrievals of ~weekly snow depth at 1 km² resolution for September 2016 through August 2018 over the Northern Hemisphere mountains are based on the temporal changes in the Sentinel-1 backscatter polarization ratio ($\sigma_{vh}^0/\sigma_{vv}^0$) and scaled to the range of snow depth measurements at in situ sites (Fig. 2), as discussed in the Methods section. The temporal and spatial resolutions of the ~weekly 1 km² retrievals meet the requirements for watershed-scale applications in mountain regions⁴, and with the Sentinel-1 follow-on missions (1C and 1D) multi-decadal data will be provided for trend analysis. Figure 4 illustrates the Sentinel-1 snow depth averaged over February 2018. We select the month of February as it corresponds with the climatological maximum snow-cover extent according to the WMO Global Cryosphere Watch. The insets illustrate details in snow depth variability in the western United States and Hindu-Kush Himalaya (Fig. 4b, c).

An assessment based on measurements at 4175 locations (Fig. 1b, Table 1 and Methods) reveals the effectiveness of the Sentinel-1 retrievals to capture the overall spatial snow depth

variability. For February 2018, the spatial correlation ($R_s$) with monthly-averaged snow depths from all in situ measurement sites is 0.76 (other months in Table 2). Both the variability between and within mountain ranges is accurately captured. For instance, when averaging snow depths for 11 selected mountain ranges (Fig. 1b), the $R_s$ between these ranges is 0.92. Within these mountain ranges, the $R_s$ is 0.72 on average. The latter value is lower due to the limited number of sites in some ranges and the potentially large differences in spatial representativeness. More specifically, measurement sites suitable for instrumentation are typically in relatively flat areas and may not represent the snow conditions of the surrounding slopes³⁶ that are included in the satellite footprint. The spatial variability in snow depth characterized by Sentinel-1 C-band SAR data is a clear improvement over currently available global estimates from reanalysis data, such as MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2³⁷; selected because it does not assimilate snow depth measurements and thus provides an independent comparison). Corresponding $R_s$ values for MERRA-2 in February 2018 are 0.36, 0.67, and 0.29 for the overall, in-between and within mountain-range variability, respectively.

Figure 5 displays the weekly evolution in snow depth over two consecutive winters averaged over all Northern Hemisphere mountain ranges and over the 11 focus ranges. Across the Northern Hemisphere (Fig. 5a), the retrievals accurately characterize the accumulation of (dry) snow, whereas a slight underestimation is noticeable during snow ablation from March onwards. This underestimation is likely caused by wet snow
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Monthly snow depth retrieved by Sentinel-1. **a** Snow depth (m) over the Northern Hemisphere mountains. **b** Enlargement of a part of the western US and Canada. **c** As in panel **b**, but for a part of the Hindu-Kush Himalaya region. The -weekly 1 km² resolution Sentinel-1 retrievals are averaged over the month of February, 2018</figcaption></figure>

<table><caption>**Table 1 Description of the in situ measurements**</caption><thead><tr><th>Region</th><th>Network name</th><th>Provider</th><th>N sites</th><th>Median elevation</th><th>References</th></tr></thead><tbody><tr><td>Global</td><td>SYNOP</td><td>ECMWF</td><td>894</td><td>432</td><td>51</td></tr><tr><td>Global</td><td>GHCN-D</td><td>NOAA</td><td>3303</td><td>1024</td><td>59,60</td></tr><tr><td>USA</td><td>SNOTEL</td><td>USDA-NRCS</td><td>733</td><td>2409</td><td>-</td></tr><tr><td>USA</td><td>SNOLITE</td><td>USDA-NRCS</td><td>16</td><td>1982</td><td>-</td></tr><tr><td>USA & Canada</td><td>Snow course/aerial marker</td><td>USDA-NRCS</td><td>615</td><td>1981</td><td>-</td></tr><tr><td>USA & Canada</td><td>Cooperator snow sensors</td><td>USDA-NRCS</td><td>32</td><td>1591</td><td>-</td></tr><tr><td>Canada</td><td>ASWS</td><td>British Columbia Government</td><td>50</td><td>1510</td><td>-</td></tr><tr><td>Nepal</td><td>-</td><td>ICIMOD, DHM, KU, Utrecht University, NVE Norway</td><td>9</td><td>4888</td><td>61-63</td></tr><tr><td>France</td><td>-</td><td>Météo France</td><td>129</td><td>1713</td><td>-</td></tr><tr><td>Austria</td><td>-</td><td>ZAMG</td><td>209</td><td>864</td><td>-</td></tr><tr><td>Austria</td><td>-</td><td>TIWAG</td><td>9</td><td>1920</td><td>64,65</td></tr><tr><td>Switzerland</td><td>BEOB, IMIS</td><td>WSL SLF</td><td>407</td><td>1560</td><td>66</td></tr></tbody></table>

Details of the global and regional in situ networks, including the regional coverage, network name (if applicable), data provider, number (N; dimensionless) of sites, median elevation (m), and references (if applicable). The number of sites is reported before averaging into 1 km² EASE-2 grid cells and the total thus exceeds 4175

<table><caption>**Table 2 Correlations between measured and retrieved snow depths**</caption><thead><tr><th>Month</th><th>Overall R<sub>s</sub> (-)</th><th>R<sub>s</sub> (-) between mountain ranges</th><th>R<sub>s</sub> (-) within mountain ranges</th></tr></thead><tbody><tr><td>December</td><td>0.67 (&lt;0.01)</td><td>0.96 (&lt;0.01)</td><td>0.58 (&lt;0.01)</td></tr><tr><td>January</td><td>0.72 (&lt;0.01)</td><td>0.93 (&lt;0.01)</td><td>0.60 (0.02)</td></tr><tr><td>February</td><td>0.76 (&lt;0.01)</td><td>0.92 (&lt;0.01)</td><td>0.72 (&lt;0.01)</td></tr><tr><td>March</td><td>0.69 (&lt;0.01)</td><td>0.85 (&lt;0.01)</td><td>0.67 (&lt;0.01)</td></tr><tr><td>April</td><td>0.72 (&lt;0.01)</td><td>0.84 (&lt;0.01)</td><td>0.60 (0.04)</td></tr></tbody></table>

Spatial correlations (R<sub>s</sub>, dimensionless) and their statistical significance level (p; in parentheses) are shown for monthly-averaged Sentinel-1 and in situ snow depth from December 2017 through April 2018. To calculate the overall R<sub>s</sub>, all available measurements across the Northern Hemisphere are used. The R<sub>s</sub> between mountain ranges is calculated from snow depths averaged per mountain range for the 11 ranges shown in Fig. 1b. The R<sub>s</sub> within mountain ranges is the average of the spatial correlations within each of these ranges
<figure><img src="image_8.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>i</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>j</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>k</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>l</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Fig. 5** Time series of snow depth measurements and retrievals. **a** The average weekly snow depth (m) over all measurement sites and corresponding Sentinel-1 grid cells across the Northern Hemisphere. **b–l** As in **a**, but averaged over the focus mountain ranges (Fig. 1b). The performance metrics shown in the top left of each figure include the temporal correlation ($R_t$; dimensionless), the February spatial correlation ($R_{s2}$; dimensionless), and the number of sites ($N$; dimensionless) used to calculate the metrics</figcaption></figure>

conditions partly reflecting and absorbing the radar signal. However, the underestimation in spring does not occur in all mountain ranges. For example, the temporal evolution in the US Sierra Nevada (Fig. 5d) during 2017–2018, including the spring melt phase, is accurate. Similarly, no underestimation is noticeable in the San Juan Mountains (Colorado, US; Fig. 5e) and Southern Scandinavian mountains (Norway; Fig. 5k). Future research shall therefore investigate under which precise conditions wet snow could cause underestimation. In some ranges, such as the Canadian Rocky Mountains (Fig. 5b) and US Sierra Nevada (Fig. 5d), the evolution of snow depth in the winter of 2016–2017 is less accurate than in 2017–2018, due to the smaller number of Sentinel-1 observations.

Figure 6 shows histograms of performance metrics that evaluate the correspondence between time series of weekly snow depth measurements and Sentinel-1 retrievals across the Northern Hemisphere mountain area, from September 2016 through August 2018. Temporal correlation coefficients ($R_t$; gray

histogram in Fig. 6a) most frequently range between 0.8 and 0.9, have a mean of 0.77 (with significance level *p* < 0.01), and are distributed with a negative skew. When zero snow depths are excluded from the analysis (white histogram), the associated reduction in the snow depth range decreases the mean correlation ($R_t$ = 0.65; *p* < 0.01). The lowest temporal correlations occur in sites with shallow snow depth and short or intermittent snow cover. The mean-absolute-error (MAE) histogram (Fig. 6b) has a positive skew, with values most frequently below 0.1 m and a mean of 0.18 m. When we exclude zero values from the analysis, the average MAE increases to 0.31 m. The bias (Fig. 6c) remains within ± 0.1 m for most of the sites, has a mean of −0.01 m and a nearly normal distribution. Large errors (MAE and bias) can again originate from the differences in representativeness between point-scale measurements and grid-scale retrievals.

For reference, the most accurate snow depth estimates today are likely provided by airborne lidar, offering decimeter-scale accuracy in surface height (either snow-on or snow-off) over local
<figure><img src="image_20.png" /><figcaption>**Fig. 6** Histograms of retrieval performance metrics. **a–c** Histograms of temporal correlation ($R_t$; dimensionless), mean-absolute-error (MAE; m), and mean of retrievals minus measurements (bias; m), calculated from time series of weekly snow depth measurements and corresponding Sentinel-1 retrievals, for September 2016 through August 2018. Gray histograms are calculated with all reported measurements, including zero values (Incl. SD = 0); white histograms exclude zero values (Excl. SD = 0), with SD the snow depth (m). Metrics indicate the mean value ($\mu$) and the number of sites ($N$; dimensionless)</figcaption></figure>

areas¹⁴. As an example, lidar snow depth estimates from the Airborne Snow Observatory have an MAE of ~0.08 m (at the 15 m × 15 m scale) with respect to manual in situ measurements over a relatively flat area near Tioga Pass in the Sierra Nevada, California, USA⁵. The Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) snow depth retrievals, available at global-scale, have correlations of ~0.25 and ~0.41 with estimates from the Snow Data Assimilation System (SNODAS) and WMO in situ measurements¹², and an MAE of 0.20 m against WMO in situ measurements¹⁰. However, a direct comparison with the Sentinel-1 retrieval performance is hampered by the limitation of the AMSR-E retrievals to areas without complex topography and with snow depths typically below 0.8 m.

**Inter-annual differences stratified by elevation.** Characterizing the inter-annual change in snow depth and how it varies with elevation is critical to assess the vulnerability of mountain systems to climate change. For snow cover, there is considerable evidence that global warming reduces the area and duration, thereby reinforcing climate change³⁸. Recent studies indicate that this impact may be stratified by elevation, with projections of decreased snow cover at low or medium elevation, but no clear trend at high elevation⁴,¹³,³⁹,⁴⁰. Whether snow depth (and mass) are similarly affected remains uncertain, particularly in mountain areas where adequate observations are missing⁴,⁹. The most accurate trend estimates of mountain snow depth are currently provided by model reanalysis, although with only moderate consistency between different reanalysis products⁴¹,⁴². Unfortunately, within the complex topography of a mountain range, reanalysis data are too coarse for a stratification by elevation to be meaningful.

Figure 7 illustrates the benefit of Sentinel-1 over the European Alps and the US Sierra Nevada for two consecutive years with contrasting snow conditions. In the European Alps, snow depth in February 2018 greatly exceeded that of February 2017. The Alps were hit by several episodes of extreme snowfall in January 2018, caused by a low-pressure area over the western Mediterranean that brought moist air northwards and resulted in the anomalously high snow depths. The Sierra Nevada featured exceptionally deep snow in February 2017, caused by a series of atmospheric river events⁴³,⁴⁴, whereas the snow depth in

February 2018 was relatively low. Over both the Alps and Sierra Nevada, similar large-scale patterns in snow depth differences (February 2018 minus 2017) are seen in the Sentinel-1 retrievals (Fig. 7a, d) and MERRA-2 reanalysis data (Fig. 7b, e). But the 1 km² Sentinel-1 retrievals provide a much more detailed representation, showing a gradual increase with elevation that closely follows that of the measurements and is poorly represented in MERRA-2 (Fig. 7c, f). Averaged over the 11 selected mountain ranges, the $R_s$ between snow depth measurements and elevation is 0.58. A similar $R_s$ is observed for the Sentinel-1 retrievals (0.50), whereas the MERRA-2 data are overall much less constrained by elevation ($R_s$ of 0.12).

**Mountain range snow volumes.** Worldwide, mountain snow-packs seasonally store and release substantial amounts of water⁴⁵. Observations of the total snow volume (depth × area) in a mountain range can be used to estimate this water storage as a mass, if snow density is known. At the global scale, snow density can for instance be estimated from information on the snow climate class and mechanical compaction, depending on snow depth and day of the year⁴⁶,⁴⁷. Figure 8 shows the snow volume (km³) retrievals for February 2018 from Sentinel-1 for the top 100 snowiest mountain ranges in the Northern Hemisphere, excluding glaciated areas⁴⁸. The largest snow volume (~380 km³) by far is in the Coast Mountains in western Canada, featuring a deep snowpack over a vast area. Given the maritime snow climate of the Coast Mountains, the February snow density is ~310 kg m⁻³. Correspondingly, the first-order estimate of the total water mass in this mountain range is ~11.8 × 10¹⁰ ton (or equally, ~400 kg m⁻²). Averaging the relatively few (48) measurements in this area results in a volume of 260 km³ (black cross in Fig. 8), considerably below the 380 km³ area-wide snow volume retrieval from Sentinel-1 (bar). This significantly lower estimate is likely due to the locations of the in situ measurements. About 75% of the sites are located in the southern half of the mountain range in a warmer climate. Moreover, 70% are located below the mean elevation (964 m) of the mountain range, with the highest site (1835 m) much below the highest mountain peak (4019 m). Confidence in the area-wide Sentinel-1 snow volume retrievals is supported by the close agreement between the cross-masked Sentinel-1 retrievals (i.e., averaged only over grid cells that include measurement sites; purple dots in Fig. 8) and the measurements (black crosses).
<figure><img src="image_21.png" /><figcaption>**Fig. 7** Inter-annual snow depth differences stratified by elevation. **a** Snow depth differences (February 2018 minus February 2017; in m) according to 1 km² Sentinel-1 retrievals over the European Alps. **b** As in a, but for coarse-resolution (0.625° × 0.5°) MERRA-2 reanalysis data. **c** Snow depth (m) of February 2017 and February 2018 stratified by elevation (m), for in situ measurements, Sentinel-1 retrievals and MERRA-2 reanalysis data (subsampled to the 1 km² grid). **d–f** As in panels **a–c**, but for the US Sierra Nevada mountain range. Source data for c,f are provided as a Source Data file</figcaption></figure>

For most of the other featured mountain ranges in North America, snow volumes compare reasonably well between area-wide and cross-masked Sentinel-1 retrievals and in situ measurements. Mountain ranges with noticeable differences are mostly limited to those with few measurement sites, such as the Mackenzie and Chugach Mountains and the Alaska Range. Similar findings apply to the European ranges, whereas more contrasting results are obtained in Asia. Here, insufficient measurements preclude confident estimates of the total snow volume. A clear example is the Himalaya, where the local variability in snow conditions is huge and the few (6) clustered measurement sites in Langtang Valley, Nepal, do not represent even the cross-masked Sentinel-1 snow volume retrieval, much less the area-wide mountain range Sentinel-1 retrieval. This highlights the need for extending measurement networks in High Mountain Asia, where water resources from snowmelt are critical<sup>8,49</sup>. At the same time, it also stresses the value of the frequent and systematic C-band SAR observations from Sentinel-1 to quantify snow depth in areas where expanding the network is difficult owing to extreme altitude, poor accessibility or safety.

## Discussion

The Sentinel-1 snow depth retrievals are primarily depending on the cross-polarized backscatter measurements. Temporal variations in $\sigma_{vh}^0$ with snow depth evolution have, to our knowledge, not been investigated before with C-band satellite observations. Only few studies deployed tower-mounted radar instruments, with contradictory results. Overall, our study aligns well with the scatterometer measurements over a site in Michigan, USA, revealing an increase in $\sigma_{vh}^0$ with an increase in snow depth<sup>26</sup>.

However, the scatterometer measurements showed signs of saturation for depths exceeding ~60 cm, which is not observed in the Sentinel-1 measurements. Potential causes for this discrepancy include the use of artificial snow in the scatterometer experiments (characterized by a homogeneous layer of snow composed of small, rounded particles and the absence of snow melt-freeze metamorphism), or differences in the spatial support of the measurements (18 m × 30 m for the scatterometer versus 1 km² for the processed Sentinel-1 data). In strong contrast, a (minor) decrease in $\sigma_{vh}^0$ with increasing snow depth was observed from scatterometer measurements over a site in the Swiss Alps<sup>27</sup>. An inverse relationship can occur with site-specific ground, vegetation and snow conditions, if the attenuation of ground scattering by the snowpack is stronger than the scattering contribution from the snowpack<sup>25</sup>. However, the results of the study in the Swiss Alps could also be impacted by the backscatter measurement principle: the total backscatter was calculated by integrating the scattering contributions from a few dominant surfaces, i.e., the snow surface, the ground surface and/or horizontal layers within the snowpack, and was thus not including the multiple scattering<sup>27</sup>.

In agreement with the scatterometer measurements in Michigan<sup>26</sup> and with our Sentinel-1 observations, radiative transfer model simulations generally indicate an increase in $\sigma_{vh}^0$ with snow depth. Recent model developments tend towards a relatively weak dependence on frequency<sup>34</sup>. A critical aspect in this context is the development of state-of-the-art, theoretically-based radiative transfer models, which allow for simulating volume scattering from snow, represented by clustered, non-spherical particles<sup>17,33–35</sup>. This presents a major improvement over conventional solutions for scattering from individual, spherical snow
<figure><img src="image_22.png" /><figcaption>**Fig. 8** Snow volumes for the top 100 snowiest mountain ranges. The snow volumes (km³) are averaged estimates over February 2018 and exclude glaciated areas. Bars represent the area-wide Sentinel-1 estimates averaged over the entire mountain range, black crosses the average of in situ measurements, and purple dots the corresponding average of cross-masked Sentinel-1 estimates, averaged only over grid cells that include in situ measurement sites. The number of in situ measurement sites, if any, is indicated next to the mountain range name. Source data are provided as a Source Data file</figcaption></figure>

particles using Rayleigh theory, assuming that scattering is proportional to the fourth power of frequency and third power of grain size³⁵. Most of these advances in radiative transfer modeling were thus far focusing on microwave measurements at higher frequencies (e.g., X-band and Ku-band). We encourage future modeling efforts to unravel the different mechanisms that cause the demonstrated sensitivity to deep snow at C-band.

To further improve the Sentinel-1 retrieval algorithm, we recommend future tower-mounted radar experiments, for instance to investigate the impacts of snow wetness on $\sigma_{vh}^0/\sigma_{vv}^0$. Such measurements would support a quantitative assessment of the uncertainty in the snow depth retrievals caused by wet snow, which in turn may result in a correction for snow liquid water content in the retrieval algorithm. In fact, previous studies demonstrated the retrieval of the snow wet-dry state²³ (and even liquid water content⁵⁰) directly from C-band SAR observations.

Wet snow maps derived from Sentinel-1 could be provided as auxiliary information (quality flag) with the snow depth retrievals. Similarly, the retrievals may be further refined by accounting for the temporal evolution in snow microstructure¹⁰. For instance, a climatological snow classification with associated snow microstructure and density estimates is available at the global scale⁴⁶,⁴⁷ and could be used to improve the parameterization. Further research is also recommended to investigate the retrieval performance in glaciated areas. Currently, data from the Randalph Glacier Inventory 6.0⁴⁸ are processed and provided with the snow depth retrievals to allow for an optional masking of glaciated areas.

In the future, sensor synergies can be exploited by merging our data over mountain areas with passive microwave observations¹⁰⁻¹² outside these areas, providing global coverage, or with active observations at higher frequency (e.g., X-band or Ku-band) that are
likely more favorable for areas with shallow snow<sup>17-19</sup>. The combination with estimates from SAR interferometry<sup>20,21</sup> could further improve the representation of temporal changes. Finally, the assimilation of the snow depth retrievals into a land surface model could yield improved and continuous (in time and space) global estimates of snow depth and mass that may also improve the initialization of numerical weather and climate models and thus improve predictions<sup>51</sup>. For such applications, the long-term continuity in C-band backscatter measurements with the Sentinel-1 constellation and the RADARSAT Constellation Mission is a strong asset, offering the systematic observations that are required to improve the representation of the cryosphere.

## Methods

**Study domain.** We processed Sentinel-1 observations for the period September 1, 2016 through August 31, 2018 over the Northern Hemisphere mountain areas north of 20°N (Fig. 1a). To delineate mountainous areas, we used the Global Mountain Biodiversity Assessment (GMBA) inventory<sup>30</sup>, providing features such as geographic coordinates, area, name, and bioclimatic region for more than 1000 individual mountain ranges in the world, which are classified based on the combination of elevation and ruggedness.

**In situ measurements.** Daily, point-scale measurements of snow depth were assembled over mountainous areas for the period September 1, 2016 through August 31, 2018. The measurements originate from a variety of sources, including regional, national and global networks (Table 1). Quality control was applied to the measurements by eliminating values higher than twice the 90th-percentile of the time series (after excluding zero values), and eliminating sites with less than 3 reported values during the evaluation period. The measurement locations were projected onto the 1 km² global cylindrical Equal Scalable Earth version 2.0 (EASE-2) grid<sup>52</sup>. Where several locations fell within the same grid cell, their time series measurements were averaged. This resulted in 4175 unique grid cells with in situ measurements (Fig. 1b).

**Sentinel-1 data processing.** The processing of the ground-range detected (GRD) Sentinel-1A and Sentinel-1B Interferometric Wide Swath (IW) σ⁰ measurements in vv-polarization and vh-polarization was performed using Google Earth Engine's Python api. We applied the standard Sentinel-1 processing techniques, including thermal noise correction, radiometric calibration, and range-Doppler terrain correction (Fig. 2). Further, the σ⁰ data with native 5 m × 20 m resolution and 10 m × 10 m grid spacing were resampled (by averaging in linear scale) and projected onto the 1 km² global cylindrical EASE-2 grid. The 1 km² scale of the resampled Sentinel-1 data matches the resolution of the auxiliary information on snow and land cover used in the retrieval algorithm (Section Sentinel-1 snow depth retrieval algorithm). Moreover, the averaging reduces the speckle noise inherent in radar observations and improves the signal-to-noise ratio, which is typically low for the cross-polarized σ⁰ data. Future research will investigate the potential of retrieving snow depths at a higher spatial resolution.

Each Sentinel-1A and Sentinel-1B satellite have an exact 12-day repeat cycle, with 175 orbits per cycle. In a given 1 km² EASE-2 grid cell, the Sentinel-1 observations from the different orbits within one repeat cycle have different incidence angles, ranging from 29.1° to 46.0° relative to a flat surface. Over mountain areas, the impact of the incidence angle on σ⁰ may be large, e.g., when a terrain slope is facing towards or away from the sensor line of sight in ascending and descending overpasses. To account for the effect of the incidence angle, we separated the σ⁰ values pertaining to each of the 175 orbits in a 12-day cycle (repeated for every cycle, with identical incidence angles for each of the orbits). The static bias between σ⁰ values from different orbits was then removed by rescaling the mean σ⁰ of each orbit to the overall mean (i.e., of the entire σ⁰ time series, including all orbits) and applying this mean correction to the individual σ⁰ measurements. The accuracy of the incidence-angle normalization increases as more observations become available. Outliers were removed by excluding values that are 3 dB above the 90th-percentile or 3 dB below the 10th-percentile of the time series. Four sub-sets of Sentinel-1 data (i.e., ascending (6 p.m.) and descending (6 a.m.) data from Sentinel-1A and Sentinel-1B) were pre-processed separately and combined into a single Sentinel-1 dataset.

Both Sentinel-1 satellites share the same orbital plane with a 6-day offset, thus the two-satellite constellation offers an exact 6-day repeat cycle. However, the observation frequency during the period considered here varies from as much as ~daily (for certain areas in Europe) to every ~2 weeks (Fig. 1a), depending on the latitude, the availability of ascending and descending orbits, the availability of Sentinel-1B data (reaching full capability in early 2017), and the evolving Sentinel-1 observation acquisition strategy. Prior to 2017, observations can (regionally) be irregular in time, with extended no-data gaps. For instance, relatively few observations were available over parts of the USA and Canada for September-December 2016 (Fig. 3a-c). We did not process any Sentinel-1 data taken before September 2016 owing to the extremely limited coverage outside Europe.

**Sentinel-1 snow depth retrieval algorithm.** The snow depth (m) retrieval algorithm relies on an empirical change detection method applied to the Sentinel-1 measurements of the cross-polarization ratio (σ⁰<sub>vh</sub>/σ⁰<sub>vv</sub>, in dB). As auxiliary input data, it uses 1 km² snow cover (SC; 1 if present or 0 if absent) from IMS and fractional forest cover (of evergreen species) from the 1 km² global consensus land cover dataset<sup>53</sup>. Firstly, a change detection index (hereafter referred to as snow index, SI; in dB) is calculated for each location *i* and time step *t*. This index links the temporal changes in Sentinel-1 σ⁰<sub>vh</sub>/σ⁰<sub>vv</sub> with the accumulation or ablation of snow:

$$ \text{SI}(i, t) = \begin{cases} \max(0, [\text{SI}(i, t-1) + \sigma_{\text{vh}}^0 / \sigma_{\text{vv}}^0 (i, t) - \sigma_{\text{vh}}^0 / \sigma_{\text{vv}}^0 (i, t-1)]) & \text{if } \text{SC}(i, t) = 1 \\ 0 & \text{if } \text{SC}(i, t) = 0 \end{cases} \quad (2) $$

where the maximum operator avoids negative indices. Secondly, the SI is rescaled into snow depth (SD; in m) as:

$$ \text{SD}(i, t) = \left( \frac{a}{1 - b\text{FC}(i)} \right) \text{SI}(i, t) \quad (3) $$

where the parameters *a* (in m dB⁻¹) and *b* (dimensionless) are constant in space and time, and FC(*i*) is the evergreen forest cover fraction (dimensionless). Forests typically attenuate snow backscatter<sup>54</sup>; they scatter part of the incoming energy, thus preventing it from reaching the snow, and part of the signal backscattered from the snow, thus preventing it from reaching the satellite sensor. This attenuation is corrected as in established coarse-scale passive microwave retrievals¹², with the same parameter value *b* = 0.6. The performance of the algorithm can potentially be further improved by optimization of *b*. Moreover, vegetation dynamics can have a very similar impact as snow on backscatter<sup>55</sup>. Therefore, the snow depth retrievals are masked when climatological leaf area index (LAI) from the Moderate Resolution Imaging Spectroradiometer (MOD15A2) exceeds 25% of the dynamic range in the LAI climatology time series. Finally, the snow depth retrievals based on Eq. (3) are smoothed to further reduce the impact of Sentinel-1 observation noise and to reduce short-term, high-magnitude fluctuations¹⁰. This is done using linear inverse distance weighting with a 2 km radius in space and a 10 day radius in time. Depending on user requirements, the smoothing can be adjusted or switched off (e.g., for data assimilation).

The parameter *a* (Eq. (3)), impacting the magnitude of the Sentinel-1 snow depth retrievals, was estimated based on in situ measurements. More specifically, *a* was optimized by minimizing the MAE between the times series of the global average snow depth measurements and corresponding Sentinel-1 retrievals. For the optimization, we excluded the time periods from March to August to limit the impact of wet snow, which is not yet accounted for in the algorithm and would likely increase errors in the retrievals. The in situ measurement locations (Section In situ measurements) were randomly sampled in two subsets with equal size: one for calibration and the other for validation. The calibrated value of *a* is 1.1 m dB⁻¹.

**Validation.** Several performance metrics are calculated to validate the Sentinel-1 snow depth retrievals using in situ measurements. As the retrieval algorithm includes the rescaling parameter *a* (Eq. 3) that is optimized based on a subset of the in situ measurements (i.e., the calibration subset), some metrics are only calculated using the remaining subset of the measurements (i.e., the validation subset) to provide an independent assessment. The spatial Pearson correlation coefficient (*R*<sub>s</sub>, dimensionless) is calculated between monthly-averaged point-scale in situ snow depth measurements and Sentinel-1 retrievals for corresponding grid-cell locations (*N* = 4175). The calculation of *R*<sub>s</sub> uses data from both the calibration and validation subsets, because the calibrated parameter *a* does not impact the metric. Similarly, temporal correlation coefficients (*R*<sub>t</sub>; dimensionless) between time series of in situ snow depth measurements and corresponding Sentinel-1 retrievals are also calculated for sites pertaining to both the calibration and validation subsets for the same reason. However, we included only sites (*N* = 1800) with more than 25 (*i.e.*, once every 4 weeks, on average) non-zero values, since the calculation of *R*<sub>t</sub> is sensitive to the sample size. The MAE and bias (mean of retrievals minus measurements) are calculated only for sites (*N* = 2000) of the validation subset with at least one non-zero value, as these metrics are less dependent on the sample size. This also allows evaluating the performance for sites with occasional snow or few reported values (such as the monthly USDA Aerial Markers).

We use point-scale in situ measurements for evaluating the performance of the Sentinel-1 snow depth retrievals. However, this likely provides a conservative estimate of the true retrieval performance. Especially in mountain areas, the point-scale snow measurements do not necessarily resemble the 1 km² grid-cell average snow conditions represented by Sentinel-1. The local variability in conditions can be large due to differences in elevation, slope and aspect, as well as wind and vegetation impacts on snow distribution<sup>56</sup>. In situ sites are preferentially located in relatively flat and non-forested terrain that is often not representative of the large variations in slope and forest cover in the surrounding area<sup>36</sup>. Consequently, in situ measurements often poorly represent the snow accumulation and melt rates on nearby slopes<sup>9,57</sup>. Ideally, the validation of the Sentinel-1 retrievals would therefore be performed using data that can be processed at the matching scale, such as gridded estimates from airborne lidar<sup>5,14</sup>, regional model simulations<sup>41</sup> or high-resolution model reanalysis<sup>58</sup>. However, such information is not available at the
global scale, which limits this study to the comparison with in situ measurements, despite the representativeness differences.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

**Data availability**

The Sentinel-1 snow depth retrievals are available online at https://ees.kuleuven.be/project/c-snow. The source data for Figs. 1a, 4, 5, 7a, d correspond to the Sentinel-1 snow depth retrievals provided through the above-mentioned website. The source data underlying Figs. 7c, f and 8 are provided as a Source Data file.

**Code availability**

The source code is available from the corresponding author upon reasonable request.

Received: 5 April 2019; Accepted: 18 September 2019;

Published online: 11 October 2019

**References**

1. Marks, D. & Dozier, J. Climate and energy exchange at the snow surface in the Alpine Region of the Sierra Nevada: 2. Snow cover energy balance. *Water Resour. Res.* **28**, 3043–3054 (1992).

2. Henderson, G. R., Peings, Y., Furtado, J. C. & Kushner, P. J. Snow-atmosphere coupling in the Northern Hemisphere. *Nat. Clim. Change* **8**, 954–963 (2018).

3. Barnett, T. P., Adam, J. C. & Lettenmaier, D. P. Potential impacts of a warming climate on water availability in snow-dominated regions. *Nature* **438**, 303–309 (2005).

4. Bormann, K. J., Brown, R. D., Derksen, C. & Painter, T. H. Estimating snow-cover trends from space. *Nat. Clim. Change* **8**, 924–928 (2018).

5. Painter, T. H. et al. The Airborne Snow Observatory: fusion of scanning lidar, imaging spectrometer, and physically-based modeling for mapping snow water equivalent and snow albedo. *Remote Sens. Environ.* **184**, 139–152 (2016).

6. Immerzeel, W. W., Droogers, P., de Jong, S. M. & Bierkens, M. F. P. Large-scale monitoring of snow cover and runoff simulation in Himalayan river basins using remote sensing. *Remote Sens. Environ.* **113**, 40–49 (2009).

7. Immerzeel, W. W., van Beeck, L. P. H. & Bierkens, M. F. P. Climate change will affect the Asian water towers. *Science* **328**, 1382–1385 (2010).

8. Best, J. Anthropogenic stresses on the world's big rivers. *Nat. Geosci.* **12**, 7–21 (2019).

9. Dozier, J., Bair, E. H. & Davis, R. E. Estimating the spatial distribution of snow water equivalent in the world's mountains. *WIREs Water* **3**, 461–474 (2016).

10. Kelly, R. E., Chang, A. T., Tsang, L. & Foster, J. L. A prototype AMSR-E global snow area and snow depth algorithm. *IEEE Trans. Geosci. Remote Sens.* **41**, 230–242 (2003).

11. Takala, M. et al. Estimating Northern Hemisphere snow water equivalent for climate research through assimilation of spaceborne radiometer data and ground-based measurements. *Remote Sens. Environ.* **115**, 3517–3529 (2011).

12. Tedesco, M. & Narvekar, P. S. Assessment of the NASA AMSR-E SWE product. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **3**, 141–159 (2010).

13. Smith, T. & Bookhagen, B. Changes in seasonal snow water equivalent distribution in High Mountain Asia (1987 to 2009). *Sci. Adv.* **17**, e1701550 (2018).

14. Deems, J. S., Painter, T. H. & Finnegan, D. C. Lidar measurements of snow depth: a review. *J. Glaciol.* **59**, 467–479 (2013).

15. Nolan, M., Larsen, C. & Sturm, M. Mapping snow depth from manned aircraft on landscape scales at centimeter resolution using structure-from-motion photogrammetry. *Cryosphere* **9**, 1445–1463 (2015).

16. Marti, R. et al. Mapping snow depth in open alpine terrain from stereo satellite imagery. *Cryosphere* **10**, 1361–1380 (2016).

17. Yueh, S. H. et al. Airborne Ku-band polarimetric radar remote sensing of terrestrial snow cover. *IEEE Trans. Geosci. Remote Sens.* **47**, 3347–3364 (2009).

18. Rott, H. et al. Cold regions hydrology high-resolution observatory for snow and cold land processes. *Proc. IEEE* **98**, 752–765 (2010).

19. King, J. et al. Spatio-temporal influence of tundra snow properties on Ku-band (17.2 GHz) backscatter. *J. Glaciol.* **61**, 267–279 (2015).

20. Guneriussen, T., Høgda, K. A., Johnsen, H. & Lauknes, I. InSAR for estimation of changes in snow water equivalent of dry snow. *IEEE Trans. Geosci. Remote Sens.* **39**, 2101–2108 (2001).

21. Leins, S., Wiesmann, A., Lemmetyinen, J. & Hajnsek, I. Snow water equivalent of dry snow measured by differential interferometry. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **8**, 3773–3790 (2015).

22. Kim, E. et al. NASA's SnowEx campaign: observing seasonal snow in a forested environment. In *Proc. IGARSS*, 1388–1390 (IEEE, 2017).

23. Nagler, T., Roth, H., Ripper, E., Bippus, G. & Hetzenecker, M. Advancements for snowmelt monitoring by means of Sentinel-1 SAR. *Remote Sens.* **8**, 348 (2016).

24. Bernier, M. et al. Determination of snow water equivalent using RADARSAT SAR data in eastern Canada. *Hydrol. Proc.* **13**, 3041–3051 (1999).

25. Shi, J. & Dozier, J. Estimation of snow water equivalence using SIR-C/X-SAR, part II: Inferring snow depth and particle size. *IEEE Trans. Geosci. Remote Sens.* **38**, 2475–2488 (2000).

26. Kendra, J. R., Sarabandi, K. & Ulaby, F. T. Radar measurements of snow: Experiment and analysis. *IEEE Trans. Geosci. Remote Sens.* **36**, 864–879 (1998).

27. Strozzi, T. & Mätzler, C. Backscattering measurements of Alpine snow covers at 5.3 and 35 GHz. *IEEE Trans. Geosci. Remote Sens.* **36**, 838–848 (1998).

28. Ulaby, F. T., Moore, R. K. & Fung, A. K. Microwave Remote Sensing, Active and Passive. Vol. II. Radar Remote Sensing and Surface Scattering and Emission Theory. (Artech House, Norwood, 1982).

29. Zhu, J. et al. Forward and inverse radar modeling of terrestrial snow using SnowSAR data. *IEEE Trans. Geosci. Remote Sens.* **56**, 7122–7132 (2018).

30. Körner, C. et al. A global inventory of mountains for bio-geographical applications. *Alp. Bot.* **127**, 1–15 (2017).

31. National Ice Center. IMS daily Northern Hemisphere snow and ice analysis at 1 km, 4 km, and 24 km resolutions. (National Snow and Ice Data Center, Digital media, Boulder) (2008).

32. Nghiem, S. V. & Tsai, W.-Y. Global snow cover monitoring with spaceborne Ku-band scatterometer. *IEEE Trans. Geosci. Remote Sens.* **39**, 2118–2134 (2001).

33. Du, J., Shi, J. & Rott, H. Comparison between a multi-scattering and multi-layer snow scattering model and its parameterized snow backscattering model. *Remote Sens. Environ.* **114**, 1089–1098 (2010).

34. Xu, X., Tsang, L. & Yueh, S. Electromagnetic models of co/cross polarization of Bicontinuous/DMRT in radar remote sensing of terrestrial snow at X- and Ku-band for CoReH2O and SCLP applications. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **5**, 1024–1032 (2012).

35. Chang, W., Tan, S., Lemmetyinen, J. & Tsang, L. Dense media radiative transfer applied to SnowScat and SnowSAR. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **7**, 3811–3825 (2014).

36. Grünewald, T. & Lehning, M. Are flat-field snow depth measurements representative? A comparison of selected index sites with areal snow depth measurements at the catchment scale. *Hydrol. Process.* **29**, 1717–1728 (2015).

37. Gelaro, R. et al. The modern-era retrospective analysis for research and applications, Version 2 (MERRA-2). *J. Clim.* **30**, 5419–5454 (2017).

38. Brown, R. et al. Snow, water, ice and permafrost in the arctic (SWIPA). *AMAP* 25–64 (2019).

39. Beniston, M. et al. The European mountain cryosphere: a review of its current state, trends, and future challenges. *Cryosphere* **12**, 759–794 (2018).

40. Skaugen, T., Stranden, H. B. & Saloranta, T. Trends in snow water equivalent in Norway (1931–2009). *Hydrol. Res.* **43**, 489–499 (2012).

41. Wrzesien, M. L. et al. A new estimate of North America mountain snow accumulation from regional climate model simulations. *Geophys. Res. Lett.* **45**, 1423–1432 (2018).

42. Mudryk, L. R., Derksen, C., Kushner, P. J. & Brown, R. Characterization of Northern Hemisphere snow water equivalent datasets. *J. Clim.* **28**, 8037–8051 (2015).

43. Wen, Y., Behrangi, A., Chen, H. & Lambrigtsen, B. How well were the early 2017 California atmospheric river precipitation events captured by satellite products and ground-based radars? *Q. J. R. Meteorol. Soc.* **144**, 344–359 (2018).

44. Behrangi, A., Bormann, K. & Painter, T. H. Using the Airborne Snow Observatory to assess remotely sensed snowfall products in the California Sierra Nevada. *Water Resour. Res.* **54**, 7331–7346 (2018).

45. Huning, L. & AghaKouchak, A. Mountain snowpack response to different levels of warming. *PNAS* **115**, 10932–10937 (2018).

46. Sturm, M. & Holmgren, J. A Seasonal snow cover classification system for local to global applications. *J. Clim.* **8**, 1261–1283 (1995).

47. Sturm, M. et al. Estimating snow water equivalent using snow depth data and climate classes. *J. Hydrometeorol.* **11**, 1380–1394 (2010).

48. RGI Consortium. *Randolph Glacier Inventory-a Dataset of Global Glacier Outlines: Version 6.0*. Technical report, Global Land Ice Measurements from Space, Colorado, USA (2017).

49. Bolch, T. et al. *Status change cryosphere Ext. Hindu. Kush Himalaya Reg.* https://doi.org/10.1007/978-3-319-92288-1_7 (2019).

50. Shi, J. & Dozier, J. Inferring snow wetness using C-band data from SIR-C's polarimetric synthetic aperture radar. *IEEE Trans. Geosci. Remote Sens.* **33**, 905–914 (1995).

51. de Rosnay, P., Isaksen, L. & Dahoui, M. Snow data assimilation at ECMWF. *ECMWF Newsl.* **143**, 26–31 (2015).
52. Brodzik, M. J., Billingsley, B., Haran, T., Raup, B. & Savoie, M. H. EASE-Grid 2.0: Incremental but significant improvements for Earth-gridded data sets. *ISPRS Int. J. Geo-Inf.* **1**, 32–45 (2012).

53. Tuanmu, M. N. & Jetz, W. A global 1-km consensus land-cover product for biodiversity and ecosystem modeling. *Glob. Ecol. Biogeogr.* **23**, 1031–1045 (2014).

54. Macelloni, G., Brogione, M., Montomoli, F. & Fontanelli, G. Effect of forests on the retrieval of snow parameters from backscatter measurements. *Eur. J. Remote Sens.* **45**, 121–132 (2012).

55. Vreugdenhil, M. et al. Sensitivity of Sentinel-1 backscatter to vegetation dynamics: an Austrian case study. *Remote Sens.* **10**, 1396 (2018).

56. Blöschl, G. Scaling issues in snow hydrology. *Hydrol. Process.* **13**, 3041–3051 (1999).

57. Meromy, L., Molotch, N. P., Link, T. E., Fassnacht, S. R. & Rice, R. Subgrid variability of snow water equivalent at operational snow stations in the western USA. *Hydrol. Process.* **27**, 2838–2400 (2013).

58. Margulis, S., Cortés, G., Girotto, M. & Durand, M. A Landsat-era Sierra Nevada snow reanalysis (1985–2015). *J. Hydrometeorol.* **17**, 1203–1221 (2016).

59. Menne, M. J. et al. *Global Historical Climatology Network-Daily (GHCN-Daily)*, Version 3.25. NOAA National Climatic Data Center (2012).

60. Menne, M. J., Durre, I., Vose, R. S., Gleason, B. E. & Houston, T. G. An overview of the global historical climatology network-daily database. *J. Atmos. Ocean. Tech.* **29**, 897–910 (2012).

61. Stigter, E. E. et al. Assimilation of snow cover and snow depth into a snow model to estimate snow water equivalent and snowmelt runoff in a Himalayan catchment. *Cyro sphere* **11**, 1647–1664 (2017).

62. Shea, J. M. et al. A comparative high-altitude meteorological analysis from three catchments in the Nepalese Himalaya. *Int. J. Water Resour. D.* **31**, 174–200 (2015).

63. Saloranta, T. et al. A model setup for mapping snow conditions in High-Mountain Himalaya. *Front. Earth Sci.* **7**, https://doi.org/10.3389/feart.2019.00129 (2019).

64. Krajci, P., Kirnbauer, R., Parajka, J., Schöber, J. & Blöschl, G. The Kühtai data set: 25 years of lysimetric, snow pillow, and meteorological measurements. *Water Resour. Res.* **53**, 5158–5165 (2017).

65. Schöber, J., Klebelsberg, A., Schattan, P., Helfricht, K. & Fey, C. Automatic weather station data from Weisssee snow research site (2014–2018). *PANGAEA* https://doi.org/10.1594/PANGAEA.898217 (2019).

66. Marty, C., Tilg, A. & Jonas, T. Recent evidence of large-scale receding snow water equivalents in the European Alps. *J. Hydrometeorol.* **18**, 1021–1031 (2017).

## Acknowledgements

This work was funded through the BELSPO SNOPOST and C-SNOW projects. Part of the study was performed by H.L. at NASA/GSFC, with computational resources provided by the NASA High-End Computing Program through the Center for Climate Simulation.

R.R. and M.G. were supported by the NASA Terrestrial Hydrology program. Sentinel-1A/B data are from the ESA and Copernicus Sentinel Satellites project and were processed using Google Earth Engine.

## Author contributions

H.L. and G.D.L. designed the study. M.De. processed the Sentinel-1 data. H.L., H.M., R.R., L.B., I.B., M.G., E.K. and G.D.L. contributed to the Sentinel-1 snow depth retrievals. Pd.R., M.Du., W.I., T.J., I.K., C.M., T.S., J.S. provided the in situ snow depth measurements. All co-authors contributed to the writing of the paper and the discussion and interpretation of the methods and results.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-019-12566-y.

Correspondence and requests for materials should be addressed to H.L.

Peer review information *Nature Communications* thanks Monique Bernier, Juha Lemmetyinen and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Reprints and permission information is available at http://www.nature.com/reprints

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_23.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019