Review Article: Global Monitoring of Snow Water Equivalent Using High-Frequency Radar Remote Sensing

Seasonal snow cover is the largest single component of the cryosphere in areal extent, covering an average of 46 × 106 km2 of Earth\u27s surface (31 % of the land area) each year, and is thus an important expression and driver of the Earth\u27s climate. In recent years, Northern Hemisphere spring snow cover has been declining at about the same rate (∼ −13 % per decade) as Arctic summer sea ice. More than one-sixth of the world\u27s population relies on seasonal snowpack and glaciers for a water supply that is likely to decrease this century. Snow is also a critical component of Earth\u27s cold regions\u27 ecosystems, in which wildlife, vegetation, and snow are strongly interconnected. Snow water equivalent (SWE) describes the quantity of water stored as snow on the land surface and is of fundamental importance to water, energy, and geochemical cycles. Quality global SWE estimates are lacking. Given the vast seasonal extent combined with the spatially variable nature of snow distribution at regional and local scales, surface observations are not able to provide sufficient SWE information. Satellite observations presently cannot provide SWE information at the spatial and temporal resolutions required to address science and high-socio-economic-value applications such as water resource management and streamflow forecasting. In this paper, we review the potential contribution of X- and Ku-band synthetic aperture radar (SAR) for global monitoring of SWE. SAR can image the surface during both day and night regardless of cloud cover, allowing high-frequency revisit at high spatial resolution as demonstrated by missions such as Sentinel-1. The physical basis for estimating SWE from X- and Ku-band radar measurements at local scales is volume scattering by millimeter-scale snow grains. Inference of global snow properties from SAR requires an interdisciplinary approach based on field observations of snow microstructure, physical snow modeling, electromagnetic theory, and retrieval strategies over a range of scales. New field measurement capabilities have enabled significant advances in understanding snow microstructure such as grain size, density, and layering. We describe radar interactions with snow-covered landscapes, the small but rapidly growing number of field datasets used to evaluate retrieval algorithms, the characterization of snowpack properties using radar measurements, and the refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. This review serves to inform the broader snow research, monitoring, and application communities on progress made in recent decades and sets the stage for a new era in SWE remote sensing from SAR measurements

Advancing the Monitoring Capabilities of Mountain Snowpack Fluctuations at Various Spatial and Temporal Scales

Snow is a critical water resource for the western US and many regions across the globe. However, our ability to accurately monitor changes in snow mass from satellite remote sensing, specifically its water equivalent, remains a challenge in mountain regions. No single sensor currently has the ability to directly measure snow water equivalent (SWE) from space at a spatial scale suitable for water supply forecasting in mountain environments. This knowledge gap calls for the innovative use of remote sensing techniques, computational tools, and data science methods to advance our ability to estimate mountain snowpacks across a range of spatial and temporal scales. The goal of this dissertation is to advance our capabilities for understanding snowpack across watershed-relevant spatial and temporal scales. Two research approaches were used to accomplish this goal: quantifying the physiographic controls and sensitivities of hydrologically important snow metrics and progressing our ability to use L-band interferometric synthetic aperture radar (InSAR) to measure SWE changes. First, we quantify the physiographic controls and various snowpack metrics in the Sierra Nevada using a novel gridded SWE reanalysis dataset. Such work demonstrates the complexity of snowpack processes and the need for fine-resolution snowpack information. Next, using L-band Interferometric Synthetic Aperture Radar (InSAR) from the NASA SnowEx campaign, both snow ablation and accumulation are estimated in the Jemez Mountains, NM. The radar-derived retrievals are evaluated utilizing a combination of optical snow-cover data, snow pits, meteorological station data, in situ snow depth sensors, and ground-penetrating radar (GPR). Lastly, we compare multisensor optical-radar approaches for SWE retrievals and find that moderate-resolution legacy satellite products provide sufficient results. The results of this work show that L-band InSAR is a suitable technique for global SWE monitoring when used synergistically with optical SCA data and snowpack modeling. While two distinctive methods are present in this research, they both work towards advancing our ability to understand the dynamics of mountain snowpack

Snow Depth Distribution Patterns and Consistency from Airborne Lidar Time Series

Snow provides fresh meltwater to over a billion people worldwide. Snow dominated watersheds drive western US water supply and are increasingly important as demand depletes reservoir and groundwater recharge capabilities. This motivates our inter- and intra-annual investigation of snow distribution patterns, leveraging the most comprehensive airborne lidar survey (ALS) dataset for snow. Validation results for ALS from both the NASA SnowEx 2017 campaign in Grand Mesa, Colorado and the time series dataset from the Tuolumne River Basin in the Sierra Nevada, in California, are presented. We then assess the consistency in the snow depth patterns for the entire basin (at 20-m resolution) and for subbasin regions (at 3-m resolution) from a collection of 51 ALS that span a six-year period (2013-2018) in the Tuolumne Basin. Strong correlations between ALS from different years near peak SWE confirm that spatial patterns exist between snow seasons. Year-to-year snow depth differs in absolute magnitude, but relative differences are consistent spatially, such that deep and shallow zones occur in the same location. We further show that elevation is the terrain parameter with the largest correlation to snow depth at the basin scale, and we map the expected pattern distribution for periods with similar snow-covered extents. Lastly, we show at a subbasin scale that distribution patterns are more consistent in vegetation-limited areas (bedrock dominated terrain and open meadows) compared to vegetation-rich zones (valley hillslopes and dense canopy cover). The maps of snow patterns and their consistency can be used to determine optimal locations of new long-term monitoring sites, design sampling strategies for future snow surveys, and to improve high resolution snow models

Integrating snow science and wildlife ecology in Arctic-boreal North America

Snow covers Arctic and boreal regions (ABRs) for approximately 9 months of the year, thus snowscapes dominate the form and function of tundra and boreal ecosystems. In recent decades, Arctic warming has changed the snowcover\u27s spatial extent and distribution, as well as its seasonal timing and duration, while also altering the physical characteristics of the snowpack. Understanding the little studied effects of changing snowscapes on its wildlife communities is critical. The goal of this paper is to demonstrate the urgent need for, and suggest an approach for developing, an improved suite of temporally evolving, spatially distributed snow products to help understand how dynamics in snowscape properties impact wildlife, with a specific focus on Alaska and northwestern Canada. Via consideration of existing knowledge of wildlife-snow interactions, currently available snow products for focus region, and results of three case studies, we conclude that improving snow science in the ABR will be best achieved by focusing efforts on developing data-model fusion approaches to produce fit-for-purpose snow products that include, but are not limited to, wildlife ecology. The relative wealth of coordinated in situ measurements, airborne and satellite remote sensing data, and modeling tools being collected and developed as part of NASA\u27s Arctic Boreal Vulnerability Experiment and SnowEx campaigns, for example, provide a data rich environment for developing and testing new remote sensing algorithms and retrievals of snowscape properties

Advancements in Measuring and Modeling the Mechanical and Hydrological Properties of Snow and Firn: Multi-sensor Analysis, Integration, and Algorithm Development

Estimating snow mechanical properties – such as elastic modulus, stiffness, and strength – is important for understanding how effectively a vehicle can travel over snow-covered terrain. Vehicle instrumentation data and observations of the snowpack are valuable for improving the estimates of winter vehicle performance. Combining in-situ and remotely-sensed snow observations, driver input, and vehicle performance sensors requires several techniques of data integration. I explored correlations between measurements spanning from millimeter to meter scales, beginning with the SnowMicroPenetrometer (SMP) and instruments applied to snow that were designed for measuring the load bearing capacity and the compressive and shear strengths of roads and soils. The spatial distribution of snow’s mechanical properties is still largely unknown. From this initial work, I determined that snow density remains a useful proxy for snowpack strength. To measure snow density, I applied multi-sensor electromagnetic methods. Using spatially distributed snowpack, terrain, and vegetation information developed in the subsequent chapters, I developed an over-snow vehicle performance model. To measure the vehicle performance, I joined driver and vehicle data in the coined Normalized Difference Mobility Index (NDMI). Then, I applied regression methods to distribute NDMI from spatial snow, terrain, and vegetation properties. Mobility prediction is useful for the strategic advancement of warfighting in cold regions. The security of water resources is climatologically inequitable and water stress causes international conflict. Water resources derived from snow are essential for modern societies in climates where snow is the predominant source of precipitation, such as the western United States. Snow water equivalent (SWE) is a critical parameter for yearly water supply forecasting and can be calculated by multiplying the snow depth by the snow density. In this work, I combined high-spatial resolution light detection and ranging (LiDAR) measured snow depths with ground-penetrating radar (GPR) measurements of two-way travel-time (TWT) to solve for snow density. Then using LiDAR derived terrain and vegetation features as predictors in a multiple linear regression, the density observations are distributed across the SnowEx 2020 study area at Grand Mesa, Colorado. The modeled density resolved detailed patterns that agree with the known interactions of snow with wind, terrain, and vegetation. The integration of radar and LiDAR sensors shows promise as a technique for estimating SWE across entire river basins and evaluating observational- or physics-based snow-density models. Accurate estimation of SWE is a means of water security. In our changing climate, snow and ice mass are being permanently lost from the cryosphere. Mass balance is an indicator of the (in)stability of glaciers and ice sheets. Surface mass balance (SMB) may be estimated by multiplying the thickness of any annual snowpack layer by its density. Though, unlike applications in seasonal snowpack, the ages of annual firn layers are unknown. To estimate SMB, I modeled the firn depth, density, and age using empirical and numerical approaches. The annual SMB history shows cyclical patterns representing the combination of atmospheric, oceanic, and anthropogenic climate forcing, which may serve as evaluation or assimilation data in climate model retrievals of SMB. The advancements made using the SMP, multi-channel GPR arrays, and airborne LiDAR and radar within this dissertation have made it possible to spatially estimate the snow depth, density, and water equivalent in seasonal snow, glaciers, and ice sheets. Open access, process automation, repeatability, and accuracy were key design parameters of the analyses and algorithms developed within this work. The many different campaigns, objectives, and outcomes composing this research documented the successes and limitations of multi-sensor estimation techniques for a broad range of cryosphere applications

Evaluation of snow depth retrievals from ICESat-2 using airborne laser-scanning data

The unprecedented precision of the altimetry satellite ICESat-2 and the increasing availability of high-resolution elevation datasets open new opportunities to measure snow depth in mountains, a critical variable for ecosystems and water resources monitoring. We retrieved snow depth over the upper Tuolumne basin (California, USA) for three years by differencing ICESat-2 ATL06 snow-on elevations and various snow-off elevation sources, including ATL06 and external digital elevation models. Snow depth derived from ATL06 data only (snow-on and snow-off) provided a poor temporal and spatial coverage, limiting its utility. However, using airborne lidar or satellite photogrammetry elevation models as snow-off elevation source yielded an accuracy of ~0.2 m (bias), a precision of ~0.5 m for low slopes and ~1.2 m for steeper areas, compared to eight reference airborne lidar snow depth maps. The snow depth derived from ICESat-2 ATL06 will help address the challenge of measuring the snow depth in unmonitored mountainous areas.</p

Estimating The Spatial Distribution of Snow Water Equivalent Using in situ and Remote Sensing Observations

Mountain snowpack is one of the primary surface water sources for about one-sixth of the global population. More than 75% of the total runoff originates from mountain snowpacks in the Western U.S. Snowmelt water recharges reservoirs and aquifers gradually in the melting season, providing vital water supplies for urban and agricultural areas. Therefore, accurately monitoring the spatial and temporal distribution of mountain snowpack &ndash; often measured as snow water equivalent (SWE) &ndash; is crucial for effective water management. While existing SWE estimation approaches remain highly uncertain, particularly when applied over large mountainous regions, the remotely-sensed snow data provide new opportunities to better characterize the spatial distributions of mountain snowpack. This dissertation investigates the approaches that optimally blend satellite, airborne, and ground snow observations to improve (near) real-time SWE estimation over mountainous terrain. The second chapter of this dissertation evaluates the accuracy of existing SWE estimation models in Sierra Nevada California. Five large-scale SWE datasets at fine spatial resolutions (&lt;= 1000 m) are comprehensively validated and compared with the Airborne Snow Observatory (ASO) SWE data in the Tuolumne River Basin (2013-2017), and ground snow pillow and snow course SWE observations across the Sierra Nevada (2004-2014). These SWE datasets include REC-INT, REC-ParBal, a Sierra Nevada SWE reanalysis (REC-DA), and two operational SWE datasets from the Snow Data Assimilation System (SNODAS) and the National Water Model (NWM-SWE), respectively. The results show that the REC-DA overall provides the most accurate SWE estimates across the Sierra Nevada (R2 = 0.87, MAE = 66 mm, PBIAS = 8.3%), followed by the REC-ParBal (R2 = 0.73, MAE = 83 mm, PBIAS = -6.4%), which is the least biased SWE estimates. Generally, SNODAS (R2 = 0.66, MAE = 106 mm, PBIAS = 9.3%) and REC-INT (R2 = 0.61, MAE = 131 mm, PBIAS = -28.3%) exhibit comparable but lower accuracy than the earlier mentioned two datasets, while NWM-SWE (R2 = 0.49, MAE = 142 mm, PBIAS = -25.2%) shows the least accuracy among the five SWE datasets. Given that REC-DA is not applicable in real-time, in the third chapter, a SWE data-fusion framework is developed, which integrates the historical SWE patterns derived from REC-DA into a statistically-based linear regression model (LRM) to estimate SWE in real-time. To investigate the influence of satellite-observed daily mean fractional snow-covered area (DMFSCA) on SWE estimation accuracy, two LRMs are compared: a baseline regression model (LRM-baseline) in which physiographic data and historical SWE patterns are used as independent variables, and an FSCA-informed regression model (LRM-FSCA) in which the DMFSCA from MODIS satellite imagery is included as an additional independent variable. By incorporating DMFSCA, LRM-FSCA outperforms LRM-baseline with improved R2 from 0.54 to 0.60, and reduced PBIAS from 2.6% to 2.2% in snow pillow cross-validation. The improvement in LRM-FSCA&rsquo;s performance is more significant during snow accumulation periods than during the snowmelt seasons. Compared to the ASO SWE, the LRM-FSCA explains 85% of the variance on average, which is at least 21% higher than the operational SNODAS (R2 = 0.64) and NWM-SWE (R2 = 0.33) in comparison. In chapter 4, a SWE bias correction framework (SWE-BCF) is developed that incorporates the ASO SWE and machine learning (ML) algorithms to further improve LRM SWE estimates in real-time. The performance of a wide range of commonly used machine learning algorithms is examined in the SWE-BCF including Gaussian Process Regression (GPR), Support Vector Machine (SVM), Bayesian Regularized Neural Networks (BRNN), Random Forest (RF), and Gradient Boosting Machine (GBM). The results indicate that all ML algorithms are capable of improving LRM-SWE accuracy substantially. While no single model performs significantly better than others, GPR, overall, shows the best performance with a 20% (0.14) increase in mean R2 value, a 31% (51 mm) reduction in mean RMSE, and a 61% (18.0%) reduction in absolute PBIAS compared with the original LRM using ASO SWE data for model validation. RF shows the most robust and stable performance in SWE bias correction with a 10% (0.08) increase in median R2 and a 41% (50 mm) reduction in median RMSE compared with the original LRM.</p

Developing Parameter Constraints for Radar-based SWE Retrievals

Terrestrial snow is an important freshwater reservoir with significant influence on the climate and energy balance. It exhibits natural spatiotemporal variability which has been enhanced by climate change, thus it is important to monitor on a large scale. Active microwave, or radar remote sensing has shown frequency-dependent promise in this regard, however, interpretation remains a challenge. The aim of this thesis was to develop constraints for radar based SWE retrievals which characterize and limit uncertainty with a focus on the underlying physical processes, snowpack stratigraphy, the influence of vegetation, and effects of background scattering. The University of Waterloo Scatterometer (UWScat) was used to make measurements at 9.6 and 17.2 GHz of snow and bare ground in a series of field-based campaigns in Maryhill and Englehart, ON, Grand Mesa, CO (NASA SnowEx campaign, year 1), and Trail Valley Creek, NT. Additional measurements from Tobermory, ON, and Churchill, MB (Canadian Snow and Ice Experiment) were included. The Microwave Emission Model for Layered Snowpacks, Version 3, adapted for backscattering (MEMLS3&a) was used to explore snowpack parameterization and SWE retrieval and the Freeman-Durden three component decomposition (FD3c) was used to leverage the polarimetric response. Physical processes in the snow accumulation environment demonstrated influence on regional snowpack parameterization and constraints in a SWE retrieval context with a single-layer snowpack parameterization for Maryhill, ON and a two-layer snowpack parameterization for Englehart, ON resulting in a retrieval RMSE of 21.9 mm SWE and 24.6 mm SWE, respectively. Use of in situ snow depths improved RMSE to 12.0 mm SWE and 10.9 mm SWE, while accounting for soil scattering effects further improved RMSE by up to 6.3 mm SWE. At sites with vegetation and ice lenses, RMSE improved from 60.4 mm SWE to 21.1 mm SWE when in situ snow depths were used. These results compare favorably with the common accuracy requirement of RMSE ≤ 30 mm and underscore the importance of understanding the driving physical processes in a snow accumulation environment and the utility of their regional manifestation in a SWE retrieval context. A relationship between wind slab thickness and the double-bounce component of the FD3c in a tundra snowpack was introduced for incidence angles ≥ 46° and wind slab thickness ≥ 19 cm. Estimates of wind slab thickness and SWE resulted in an RMSE of 6.0 cm and 5.5 mm, respectively. The increased double-bounce scattering was associated with path length increase within a growing wind slab layer. Signal attenuation in a sub-canopy SWE retrieval was also explored. The volume scattering component of the FD3c yielded similar performance to forest fraction in the retrieval with several distinct advantages including a real-time description of forest condition, accounting for canopy geometry without ancillary information, and providing coincident information on forest canopy in remote locations. Overall, this work demonstrated how physical processes can manifest regional outcomes, it quantified effects of natural inclusions and background scattering on SWE retrievals, it provided a means to constrain wind slab thickness in a tundra environment, and it improved characterization of coniferous forest in a sub-canopy SWE retrieval context. Future work should focus on identifying ice and vegetation conditions prior to SWE retrieval, testing the spatiotemporal validity of the methods developed herein, and finally, improving the integration of snowpack attenuation within retrieval efforts

HUMAN AND CLIMATE IMPACTS ON FLOODING VIA REMOTE SENSING, BIG DATA ANALYTICS, AND MODELING

Over the last 20 years, the amount of streamflow has greatly increased and spring snowmelt floods have occurred more frequently in the north-central U.S. In the Red River of the North Basin (RRB) overlying portions of North Dakota and Minnesota, six of the 13 major floods over the past 100 years have occurred since the late 1990s. Based on numerous previous studies as well as senior flood forecasters’ experiences, recent hydrological changes related to human modifications [e.g. artificial subsurface drainage (SSD) expansion] and climate change are potential causes of notable forecasting failures over the past decade. My dissertation focuses on the operational and scientific gaps in current forecasting models and observational data and provides insights and value to both the practitioner and the research community. First, the current flood forecasting model needs both the location and installation timing of SSD and SSD physics. SSD maps were developed using satellite “big” data and a machine learning technique. Next, using the maps with a land surface model, the impacts of SSD expansion on regional hydrological changes were quantified. In combination with model physics, the inherent uncertainty in the airborne gamma snow survey observations hinders the accurate flood forecasting model. The operational airborne gamma snow water equivalent (SWE) measurements were improved by updating antecedent surface moisture conditions using satellite observations on soil moisture. From a long-term perspective, flood forecasters and state governments need knowledge of historical changes in snowpack and snowmelt to help flood management and to develop strategies to adapt to climate changes. However, historical snowmelt trends have not been quantified in the north-central U.S. due to the limited historical snow data. To overcome this, the current available historical long-term SWE products were evaluated across diverse regions and conditions. Using the most reliable SWE product, a trend analysis quantified the magnitude of change extreme snowpack and melt events over the past 36 years. Collectively, this body of research demonstrates that human and climate impacts, as well as limited and noisy data, cause uncertainties in flood prediction in the great plains, but integrated approaches using remote sensing, big data analytics, and modeling can quantify the hydrological changes and reduce the uncertainties. This dissertation improves the practice of flood forecasting in Red River of the North Basin and advances research in hydrology and snow science