The seasonal mountain snowpack of the western U.S. (WUS) is a key water resource to millions of people and an important component of the regional climate system. Impurities at the snow surface can affect snowmelt timing and rate through snow radiative forcing (RF), resulting in earlier streamflow, snow disappearance, and less water availability in dry months. Predicting the locations, timing, and intensity of impurities is challenging, and little is known concerning whether snow RF has changed over recent decades. Here we analyzed the relative magnitude and spatio-temporal variability of snow RF across the WUS at three spatial scales (pixel, watershed, regional) using remotely sensed RF from spatially and temporally complete (STC) MODIS data sets (STC-MODSCAG/MODDRFS) from 2001-2022. To quantify snow RF impacts, we calculated a pixel-integrated metric over each snowmelt season (March 1st - June 30th) in all 22 years. We tested for long-term trend significance with the Mann-Kendall test and trend magnitude with Theil-Sen's slope. Mean snow RF was highest in the Upper Colorado region, but notable in less-studied regions, including the Great Basin and Pacific Northwest. Watersheds with high snow RF also tended to have high spatial and temporal variability in RF, and these tended to be near arid regions. Snow RF trends were largely absent; only a small percent of mountain ecoregions (0.03-8%) had significant trends, and these were typically decreasing trends. All mountain ecoregions exhibited a net decline in snow RF. While the spatial extent of significant RF trends was minimal, we found declining trends most frequently in the Sierra Nevada, North Cascades, and Canadian Rockies, and increasing trends in the Idaho Batholith. This study establishes a two-decade chronology of snow impurities in the WUS, helping inform where and when RF impacts on snowmelt may need to be considered in hydrologic models and regional hydroclimate studies.