We evaluate performance of the Catchment Land Surface Model (CLSM) under flood conditions after the assimilation of observations of the terrestrial water storage anomaly (TWSA) from NASA’s Gravity Recovery and Climate Experiment (GRACE). Assimilation offers three key benefits for the viability of GRACE observations to operational applications: (1) near-real time analysis; (2) a downscaling of GRACE’s coarse spatial resolution; and (3) state disaggregation of the vertically-integrated TWSA. We select the 2011 flood event in the Missouri river basin as a case study, and find that assimilation generally made the model wetter in the months preceding flood. We compare model outputs with observations from 14 USGS groundwater wells to assess improvements after assimilation. Finally, we examine disaggregated water storage information to improve the mechanistic understanding of event generation. Validation establishes that assimilation improved the model skill substantially, increasing regional groundwater anomaly correlation from 0.58 to 0.86. For the 2011 flood event in the Missouri river basin, results show that groundwater and snow water equivalent were contributors to pre-event flood potential, providing spatially-distributed early warning information.

This study investigates the effect of projected temperature increases on maritime mountain snowpack in the McKenzie River Basin (MRB; 3041 km²) in the Cascades Mountains of Oregon, USA. We simulated the spatial distribution of snow water equivalent (SWE) in the MRB for the period of 1989–2009 with SnowModel, a spatially-distributed, process-based model (Liston and Elder, 2006b). Simulations were evaluated using point-based measurements of SWE, precipitation, and temperature that showed Nash-Sutcliffe Efficiency coefficients of 0.83, 0.97, and 0.80, respectively. Spatial accuracy was shown to be 82% using snow cover extent from the Landsat Thematic Mapper. The validated model then evaluated the inter- and intra-year sensitivity of basin wide snowpack to projected temperature increases (2 °C) and variability in precipitation (±10%). Results show that a 2 °C increase in temperature would shift the average date of peak snowpack 12 days earlier and decrease basin-wide volumetric snow water storage by 56%. Snowpack between the elevations of 1000 and 2000 m is the most sensitive to increases in temperature. Upper elevations were also affected, but to a lesser degree. Temperature increases are the primary driver of diminished snowpack accumulation, however variability in precipitation produce discernible changes in the timing and volumetric storage of snowpack. The results of this study are regionally relevant as melt water from the MRB's snowpack provides critical water supply for agriculture, ecosystems, and municipalities throughout the region especially in summer when water demand is high. While this research focused on one watershed, it serves as a case study examining the effects of climate change on maritime snow, which comprises 10% of the Earth's seasonal snow cover.

We characterize how regional watersheds function as simple, dynamic systems through a series of hysteresis loops using measurements from NASA's Gravity Recovery and Climate Experiment (GRACE) satellites. These loops illustrate the temporal relationship between runoff and terrestrial water storage in three regional-scale watersheds (> 150 000 km2) of the Columbia River Basin, USA and Canada. The shape and size of the hysteresis loops are controlled by the climate, topography, and geology of the watershed. The direction of the hystereses for the GRACE signals moves in opposite directions from the isolated groundwater hystereses. The subsurface water (soil moisture and groundwater) hystereses more closely resemble the storage-runoff relationship of a soil matrix. While the physical processes underlying these hystereses are inherently complex, the vertical integration of terrestrial water in the GRACE signal encapsulates the processes that govern the non-linear function of regional-scale watersheds. We use this process-based understanding to test how GRACE data can be applied prognostically to predict seasonal runoff (mean Nash-Sutcliffe Efficiency of 0.91) and monthly runoff during the low flow/high demand month of August (mean Nash-Sutcliffe Efficiency of 0.77) in all three watersheds. The global nature of GRACE data allows this same methodology to be applied in other regional-scale studies, and could be particularly useful in regions with minimal data and in trans-boundary watersheds.

Snow is a critical component of the climate system, provides fresh water for millions of people globally, and affects forest and wildlife ecology. Snowy regions are typically data sparse, especially in mountain environments. Remotely-sensed snow cover data are available globally but are challenging to convert into accessible, actionable information. SnowCloudMetrics is a web portal for on-demand production and delivery of snow information including snow cover frequency (SCF) and snow disappearance date (SDD) using Google Earth Engine (GEE). SCF and SDD are computed using the Moderate Resolution Imaging Spectroradiometer (MODIS) Snow Cover Binary 500 m (MOD10A1) product. The SCF and SDD metrics are assessed using 18 years of Snow Telemetry records at more than 750 stations across the Western U.S. SnowCloudMetrics provides users with the capacity to quickly and efficiently generate local-to-global scale snow information. It requires no user-side data storage or computing capacity, and needs little in the way of remote sensing expertise. SnowCloudMetrics allows users to subset by year, watershed, elevation range, political boundary, or user-defined region. Users can explore the snow information via a GEE map interface and, if desired, download scripts for access to tabular and image data in non-proprietary formats for additional analyses. We present global and hemispheric scale examples of SCF and SDD. We also provide a watershed example in the transboundary, snow-dominated Amu Darya Basin. Our approach represents a new, user-driven paradigm for access to snow information. SnowCloudMetrics benefits snow scientists, water resource managers, climate scientists, and snow related industries providing SCF and SDD information tailored to their needs, especially in data sparse regions.

Dynamic natural processes govern snow distribution in mountainous environments throughout the world. Interactions between these different processes create spatially variable patterns of snow depth across a landscape. Variations in accumulation and redistribution occur at a variety of spatial scales, which are well established for moderate mountain terrain. However, spatial patterns of snow depth variability in steep, complex mountain terrain have not been fully explored due to insufficient spatial resolutions of snow depth measurement. Recent advances in uncrewed aerial systems (UASs) and structure from motion (SfM) photogrammetry provide an opportunity to map spatially continuous snow depths at high resolutions in these environments. Using UASs and SfM photogrammetry, we produced 11 snow depth maps at a steep couloir site in the Bridger Range of Montana, USA, during the 2019–2020 winter. We quantified the spatial scales of snow depth variability in this complex mountain terrain at a variety of resolutions over 2 orders of magnitude (0.02 to 20 m) and time steps (4 to 58 d) using variogram analysis in a high-performance computing environment. We found that spatial resolutions greater than 0.5 m do not capture the complete patterns of snow depth spatial variability within complex mountain terrain and that snow depths are autocorrelated within horizontal distances of 15 m at our study site. The results of this research have the potential to reduce uncertainty currently associated with snowpack and snow water resource analysis by documenting and quantifying snow depth variability and snowpack evolution on relatively inaccessible slopes in complex terrain at high spatial and temporal resolutions.