Seasonal snowpack is an essential component in the hydrological cycle and plays a significant role in supplying water resources to downstream users. Yet the snow water equivalent (SWE) in seasonal snowpacks, and its space–time variation, remains highly uncertain, especially over mountainous areas with complex terrain and sparse observations, such as in High Mountain Asia (HMA). In this work, we assessed the spatiotemporal distribution of seasonal SWE, obtained from a new 18-year HMA Snow Reanalysis (HMASR) dataset, as part of the recent NASA High Mountain Asia Team (HiMAT) effort. A Bayesian snow reanalysis scheme previously developed to assimilate satellite-derived fractional snow-covered area (fSCA) products from Landsat and MODIS platforms has been applied to develop the HMASR dataset (at a spatial resolution of 16 arcsec (∼500 m) and daily temporal resolution) over the joint Landsat–MODIS period covering water years (WYs) 2000–2017. Based on the results, the HMA-wide total SWE volume is found to be around 163 km³ on average and ranges from 114 km³ (WY2001) to 227 km³ (WY2005) when assessed over 18 WYs. The most abundant snowpacks are found in the northwestern basins (e.g., Indus, Syr Darya and Amu Darya) that are mainly affected by the westerlies, accounting for around 66 % of total seasonal SWE volume. Seasonal snowpack in HMA is depicted by snow accumulating through October to March and April, typically peaking around April and depleting in July–October, with variations across basins and WYs. When examining the elevational distribution over the HMA domain, seasonal SWE volume peaks at mid-elevations (around 3500 m), with over 50 % of the volume stored above 3500 m. Above-average amounts of precipitation causes significant overall increase in SWE volumes across all elevations, while an increase in air temperature (∼1.5 K) from cooler to normal conditions leads to an redistribution in snow storage from lower elevations to mid-elevations. This work brings new insight into understanding the climatology and variability of seasonal snowpack over HMA, with the regional snow reanalysis constrained by remote-sensing data, providing a new reference dataset for future studies of seasonal snow and how it contributes to the water cycle and climate over the HMA region.

Snow water equivalent (SWE), particularly in mountains regions, has been an elusive hydrologic measurement. We examine the utility of a data assimilation approach to generate space‐time continuous estimates of SWE from more readily available snow depth (SD) measurements. A multitemporal lidar data set provides a unique opportunity to assimilate single SD images and verify posterior estimates against SD images at nonassimilation times. Application over three water years shows significant improvement in the posterior estimates with an average correlation between estimated and measured SD fields of 0.88 compared to 0.52 for prior estimates that do not benefit from the assimilated SD data. We also show that posterior estimates are consistent with independent in situ SWE and streamflow measurements. This work demonstrates that using high‐resolution/high‐accuracy, but infrequent, SD measurements combined with a data assimilation framework could make significant inroads toward the goal of spatially distributed SWE and snowmelt estimates at the global scale.

Water stored in mountain snowpacks (i.e., snow water equivalent, SWE) represents an important but poorly characterized component of the terrestrial water cycle. The Western United States snow reanalysis (WUS–SR) dataset is novel in its combination of spatial resolution (~500 m), spatial extent (31°–49° N; 102°–125° W), and temporal continuity (daily over 1985–2021). WUS–SR is generated using a Bayesian framework with model-based snow estimates updated through the assimilation of cloud-free Landsat fractional snow-covered area observations. Over the WUS, the peak SWE verification with independent in situ measurements show correlation coefficient, mean difference (MD), and root mean squared difference (RMSD) of 0.77, −0.15 m, and 0.28 m, respectively. The effects of forest cover and Landsat image availability on peak SWE are assessed. WUS–SR peak SWE is well correlated (ranging from 0.75 to 0.91) against independent lidar-derived SWE taken near April 1^st, with MD <0.15 m and RMSD <0.38 m. The dataset is useful for characterizing WUS mountain snow storage, and ultimately for improving snow-derived water resources management.

This work quantifies the uncertainty of accumulation-season peak snow water storage in the portions of the midlatitude American Cordillera where snow is a dominant driver of hydrology. This is accomplished through intercomparison of commonly used global and regional products over the Western United States (WUS) and Andes domains, which have similar hydrometeorology but are disparate with respect to the amount of available in situ information. The recently developed WUS Snow Reanalysis (WUS-SR) and Andes Snow Reanalysis (Andes-SR) datasets, which have been extensively verified against in situ measurements, are used as baseline reference datasets in the intercomparison. Relative to WUS-SR climatological peak snow water equivalent (SWE) storage (269 km³), high- and moderate-resolution products (i.e., those with resolutions less than ∼10 km) are in much better agreement (284±14 km³; overestimated by 6 %) compared to low-resolution products (127±54 km³; underestimated by 53 %). In comparison to the Andes-SR peak snow storage (29 km³), all other products show large uncertainty and bias (19±16 km³; underestimated by 34 %). Examination of spatial patterns related to orographic effects showed that only the high- to moderate-resolution Snow Data Assimilation System (SNODAS) and University of Arizona (UA) products show comparable estimates of windward–leeward SWE patterns over a subdomain (Sierra Nevada) of the WUS. Coarser products distribute too much snow on the leeward side in both the Sierra Nevada and Andes, missing orographic and rain shadow patterns that have important hydrological implications. The uncertainty of peak seasonal snow storage is primarily explained by precipitation uncertainty in both the WUS (R²=0.55) and Andes (R²=0.84). Despite using similar forcing inputs, snow storage diverges significantly within the ECMWF Reanalysis v5 (ERA5) (i.e., ERA5 vs. ERA5-Land) products and the Global Land Data Assimilation System (GLDAS) (modeled with Noah, Variable Infiltration Capacity (VIC), and Catchment model) products due to resolution-induced elevation differences and/or differing model process representation related to rain–snow partitioning and accumulation-season snowmelt generation. The availability and use of in situ precipitation and snow measurements (i.e., in WUS) in some products adds value by reducing snow storage uncertainty; however, where such data are limited, i.e., in the Andes, significant biases and uncertainty exist.

Accurate estimates of snow water equivalent (SWE) based on remote sensing have been elusive, particularly in mountain areas. However, there now appears to be some potential for direct satellite-based SWE observations along ground tracks that only cover a portion of a spatial domain (e.g., watershed). Fortunately, spatiotemporally continuous meteorological and surface variables could be leveraged to infer SWE in the gaps between satellite ground tracks. Here, we evaluate statistical and machine learning (ML) approaches to performing track-to-area (TTA) transformations of SWE observations in California's upper Tuolumne River watershed using synthetic data. The synthetic SWE measurements are designed to mimic a potential future P-band Signals of Opportunity (P-SoOP) satellite mission with a (along-track) spatial resolution of about 500 m. We construct relationships between multiple meteorological and surface variables and synthetic SWE observations along observation tracks, and we then extend these relationships to unobserved areas between ground tracks to estimate SWE over the entire watershed. Domain-wide, SWE inferred on 1 April using two synthetic satellite tracks (∼4.5 % basin coverage) led to percent errors of basin-averaged SWE (PEBAS) of 24.5 %, 4.5 % and 6.3 % in an extremely dry water year (WY2015), a normal water year (WY2008) and an extraordinarily wet water year (WY2017), respectively. Assuming a 10 d overpass interval, percent errors of basin-averaged SWE during both snow accumulation and snowmelt seasons were mostly less than 10 %. We employ a feature sensitivity analysis to overcome the black-box nature of ML methods and increase the explainability of the ML results. Our feature sensitivity analysis shows that precipitation is the dominant variable controlling the TTA SWE estimation, followed by net long-wave radiation (NetLong). We find that a modest increase in the accuracy of SWE estimation occurs when more than two ground tracks are leveraged. The accuracy of 1 April SWE estimation is only modestly improved for track repeats more often than about 15 d.