In the mountainous headwaters of the Colorado River episodic dust deposition from adjacent arid and disturbed landscapes darkens snow and accelerates snowmelt, impacting basin hydrology. Patterns and impacts across the heterogenous landscape cannot be inferred from current in situ observations. To fill this gap daily remotely sensed retrievals of radiative forcing and contribution to melt were analyzed over the MODIS period of record (2001-2023) to quantify spatiotemporal impacts of snow darkening. Each season radiative forcing magnitudes were lowest in early spring and intensified as snowmelt progressed, with interannual variability in timing and magnitude of peak impact. Over the full record, radiative forcing was elevated in the first decade relative to the last decade. Snowmelt was accelerated in all years and impacts were most intense in the central to southern headwaters. The spatiotemporal patterns motivate further study to understand controls on variability and related perturbations to snow water resources.

With climate warming, many regions are experiencing changes in snow accumulation and persistence. These changes are known to affect streamflow volume, but the magnitude of the effect varies between regions. This research evaluates whether variables derived from remotely sensed snow cover can be used to estimate annual streamflow at the small watershed scale across the western U.S., a region with a wide range of climate types. We compared snow cover variables derived from MODIS, snow persistence (SP), and snow season (SS), to more commonly utilized metrics, snow fraction (fraction of precipitation falling as snow, SF), and peak snow water equivalent (SWE). Each variable represents different information about snow, and this comparison assesses similarities and differences between the snow metrics. Next, we evaluated how two snow variables, SP and SWE, related to annual streamflow (Q) for 119 USGS reference watersheds and examined whether these relationships varied for wet/warm (precipitation surplus) and dry/cold (precipitation deficit) watersheds. Results showed high correlations between all snow variables, but the slopes of these relationships differed between climates, with wet/warm watersheds displaying lower SF and higher SWE for the same SP. In dry/cold watersheds, both SP and SNODAS SWE correlated with Q spatially across all watersheds and over time within individual watersheds. We conclude that SP can be used to map spatial patterns of annual streamflow generation in dry/cold parts of the region. Applying this approach to the Upper Colorado River Basin demonstrates that 50% of streamflow comes from areas >3,000 masl. If the relationship between SP and Q is similar in other dry/cold regions, this approach could be used to estimate annual streamflow in ungauged basins.