SEARCH · Search NASA

Results for “Snow Telemetry (SNOTEL)”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

18 records

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Multitemporal Snow Cover Mapping in Mountainous Terrain for Landsat Climate Data Record Development

A multitemporal method to map snow cover in mountainous terrain is proposed to guide Landsat climate data record (CDR) development. The Landsat image archive including MSS, TM, and ETM+ imagery was used to construct a prototype Landsat snow cover CDR for the interior northwestern United States. Landsat snow cover CDRs are designed to capture snow-covered area (SCA) variability at discrete bi-monthly intervals that correspond to ground-based snow telemetry (SNOTEL) snow-water-equivalent (SWE) measurements. The June 1 bi-monthly interval was selected for initial CDR development, and was based on peak snowmelt timing for this mountainous region. Fifty-four Landsat images from 1975 to 2011 were preprocessed that included image registration, top-of-the-atmosphere (TOA) reflectance conversion, cloud and shadow masking, and topographic normalization. Snow covered pixels were retrieved using the normalized difference snow index (NDSI) and unsupervised classification, and pixels having greater (less) than 50% snow cover were classified presence (absence). A normalized SCA equation was derived to independently estimate SCA given missing image coverage and cloud-shadow contamination. Relative frequency maps of missing pixels were assembled to assess whether systematic biases were embedded within this Landsat CDR. Our results suggest that it is possible to confidently estimate historical bi-monthly SCA from partially cloudy Landsat images. This multitemporal method is intended to guide Landsat CDR development for freshwaterscarce regions of the western US to monitor climate-driven changes in mountain snowpack extent.

Landsat↗

NASA NTRSDate not supplied

Improving Modeled Snow Mass in GEOS with Snow Cover Fraction Assimilation

The NASA Goddard Earth Observing System (GEOS) Forward Processing (FP) global weather analysis often considerably underestimates near-surface air temperature during the spring snow melt season, which can be caused by excessive simulated snow cover. To mitigate such errors, we developed an empirical snow cover fraction (SCF) analysis to assimilate SCF observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) into the GEOS Catchment land surface model. Results obtained with a land-only version of the system indicate that the SCF analysis reduces errors in snow mass and snow depth by ~10% when compared to in situ measurements from the Snow Telemetry (SNOTEL), Canadian historical Snow Water Equivalent dataset (CanSWE), and Global Historical Climatology Network - Daily (GHCN-Daily) networks.

Rolf Reichle↗

DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗

DOE OSTI2023

Evaluation of E3SM land model snow simulations over the western United States

Abstract. Seasonal snow has crucial impacts on climate, ecosystems, and humans, but it is vulnerable to global warming. The land component (ELM) of the Energy Exascale Earth System Model (E3SM) mechanistically simulates snow processes from accumulation, canopy interception, compaction, and snow aging to melt. Although high-quality field measurements, remote sensing snow products, and data assimilation products with high spatio-temporal resolution are available, there has been no systematic evaluation of the snow properties and phenology in ELM. This study comprehensively evaluates ELM snow simulations over the western United States at 0.125∘ resolution during 2001–2019 using the Snow Telemetry (SNOTEL) in situ networks, MODIS remote sensing products (i.e., MCD43 surface albedo product), the spatially and temporally complete (STC) snow-covered area and grain size (MODSCAG) and MODIS dust and radiative forcing in snow (MODDRFS) products (STC-MODSCAG/STC-MODDRFS), and the snow property inversion from remote sensing (SPIReS) product and two data assimilation products of snow water equivalent and snow depth – i.e., University of Arizona (UA) and SNOw Data Assimilation System (SNODAS). Overall the ELM simulations are consistent with the benchmarking datasets and reproduce the spatio-temporal patterns, interannual variability, and elevation gradients for different snow properties including snow cover fraction (fsno), surface albedo (αsur) over snow cover regions, snow water equivalent (SWE), and snow depth (Dsno). However, there are large biases of fsno with dense forest cover and αsur in the Rocky Mountains and Sierra Nevada in winter, compared to the MODIS products. There are large discrepancies of snow albedo, snow grain size, and light-absorbing particle-induced snow albedo reduction between ELM and the MODIS products, attributed to uncertainties in the aerosol forcing data, snow aging processes in ELM, and remote sensing retrievals. Against UA and SNODAS, ELM has a mean bias of −20.7 mm (−35.9 %) and −20.4 mm (−35.5 %), respectively, for spring, and −13.8 mm (−27.8 %) and −10.2 mm (−22.2 %), respectively, for winter. ELM shows a relatively high correlation with SNOTEL SWE, with mean correlation coefficients of 0.69 but negative mean biases of −122.7 mm. Compared to the snow phenology of STC-MODSCAG and SPIReS, ELM shows delayed snow accumulation onset dates by 17.3 and 12.4 d, earlier snow end dates by 35.5 and 26.8 d, and shorter snow durations by 52.9 and 39.5 d, respectively. This study underscores the need for diagnosing model biases and improving ELM representations of snow properties and snow phenology in mountainous areas for more credible simulation and future projection of mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Hydrologic Model Files (PRMS) for Historical Conditions in the East River Watershed, Colorado between 1987-2019

The data package contains model input and output file for the East River watershed (85 km2) located in the headwaters of the Colorado River. The model is the U.S. Geological Survey (USGS) Precipitation-Modeling Runoff System (PRMS) for the historical, or baseline, condition (1987-2019) to assess the efficiency of summer monsoon rains to generate streamflow. The finite difference grid resolution is 100 m with elevations resampled from the USGS National Elevation Dataset. LANDFIRE raster data is used to derive parameters of dominant cover type, summer and winter cover density, canopy interception characteristics for snow and rain and transmission coefficient for shortwave radiation. Climate forcing uses minimum and maximum daily temperature lapse rates defined by the two Snow Telemetry (SNOTEL) stations located in the vicinity of the East River adjusted for aspect. Schofield SNOTEL snowfall is spatially distributed using LiDAR derived snow depth observations from the Airborne Snow Observatory (ASO) flown 4 April 2016. Snow depths are converted to snow water equivalent based on ground surveys and density modeling. Rainfall is spatially distributed using the monthly Parameter-elevation Relationships on Independent Slopes Model 30-year monthly averages. Simulated solar radiation is calibrated to match weather station observations. Model verification of snow water equivalent (SWE) accumulation and ablation relies on the 2018 and 2019 ASO maps and weather station snow depth. The data package contains model input files in text format, model documentation in pfd and a readme.txt file to help users understand the file structure to run the program.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Forest Canopy Density Effects on Snowpack Across the Climate Gradients of the Western United States Mountain Ranges

Limitations inherent to field studies such as limited spatial and temporal coverage has made it very challenging, if not impossible, to explicitly disentangle the effect of canopy density from climate on snowpack dynamics. Hence, our understanding of how variations in canopy density can affect snow processes under diverse climate conditions has been limited. To address this knowledge gap, this study uses a physics-based modeling approach to represent most key processes that drive snow accumulation and melt in the open and forest with varied forest density, across the climate gradients of the Western U.S. mountain ranges, as represented by 228 Snow Telemetry (SNOTEL) locations. Simulations suggest that, under most winter climate conditions, canopy density exerts a strong influence on under-canopy snowpack dynamics and their controlling processes. Changing canopy density can shift the direction of the modeled differential peak snow accumulation and snow disappearance date (SDD) between the open and forest. For example, while simulated snowpack lasts longer in the open than the high-density forest at most wet/warm and dry/cold locations by up to 55 days, about 90% of all locations show a longer snow duration in the low-density forest by up to 27 days as compared to the open. Interannual variability in climate can also change the magnitude and direction of canopy impact on snow processes. By comparing model simulations to published field observations, and we also identify conditions that require additional field observations to improve process-scale understanding of forest-snow interaction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Estimating snow cover from high-resolution satellite imagery by thresholding blue wavelengths

We report the extent and duration of snow cover, a critical component of the hydrologic cycle and the global climate system, is expected to shift dramatically under climate change. Therefore, developing high-resolution assessments of snow cover change is crucial for estimating the impact of changing snow cover on watershed and ecosystems processes in cold regions. Remote sensing tools provide a powerful method for mapping snow-covered area (SCA) across a landscape. The most common method for estimating SCA utilizes the normalized difference snow index (NDSI), which relies on spectral measurements in the shortwave-infrared wavelengths (SWIR). NDSI can effectively estimate catchment- to regional-scale SCA, but it cannot be used to assess fine-scale SCA because of current limitations on the spatial resolution of satellite-derived SWIR measurements. Here, we map SCA using a threshold of blue wavelengths and high-resolution satellite imagery. The thresholding method, which we call the Blue Snow Threshold algorithm (BST), has previously been used with digital camera imagery. We refine and automate the algorithm for use with cloud-free high-resolution satellite imagery and find that the BST can be used to assess fine-scale SCA. For validation, we compared BST-derived estimates of SCA to a) airborne lidar surveys, b) Landsat fractional SCA, and c) snow disappearance dates from Snow Telemetry (SNOTEL) stations. When compared to airborne lidar surveys of SCA, the BST predicted SCA had a range of F-scores between 0.81 and 0.94 in four study areas in California and Colorado. We also found general agreement between SCA and snow disappearance at multiple SNOTEL sites across the western United States. Given the relatively recent availability of high-resolution satellite imagery with spectral measurements in the visible wavelengths but lacking in SWIR, the BST offers a reliable and easy-to-apply tool for examining fine-scale snow-related processes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

Changing Snow Cover and Stream Discharge in the Western United States - Wind River Range, Wyoming

Earlier onset of springtime weather has been documented in the western United States over at least the last 50 years. Because the majority (>70%) of the water supply in the western U.S. comes from snowmelt, analysis of the declining spring snowpack has important implications for the management of water resources. We studied ten years of Moderate-Resolution Imaging Spectroradiometer (MODIS) snow-cover products, 40 years of stream discharge and meteorological station data and 30 years of snow-water equivalent (SWE) SNOw Telemetry (SNOTEL) data in the Wind River Range (WRR), Wyoming. Results show increasing air temperatures for.the 40-year study period. Discharge from streams in WRR drainage basins show lower annual discharge and earlier snowmelt in the decade of the 2000s than in the previous three decades. Changes in streamflow may be related to increasing air temperatures which are probably contributing to a reduction in snow cover, although no trend of either increasingly lower streamflow or earlier snowmelt was observed within the decade of the 2000s. And SWE on 1 April does not show an expected downward trend from 1980 to 2009. The extent of snow cover derived from the lowest-elevation zone of the WRR study area is strongly correlated (r=0.91) with stream discharge on 1 May during the decade of the 2000s. The strong relationship between snow cover and streamflow indicates that MODIS snow-cover maps can be used to improve management of water resources in the drought-prone western U.S.

Hall, Dorothy K.↗

NASA NTRS2018

Quantifying the Early Snowmelt Event of 2015 in the Cascade Mountains, USA by Developing and Validating MODIS-Based Snowmelt Timing Maps

Spring snowmelt serves as the major hydrological contribution to many watersheds of the US West. Since the 1970s the conterminous western USA has seen an earlier arrival of spring snowmelt. The extremely low snowpack and early melt of 2015 in the Cascade Mountains may be a harbinger of winters to come, underscoring the interest in advancements in spring snowmelt monitoring. Target-of-opportunity and point measurements of snowmelt using meteorological stations or stream gauges are common sources of these data, however, there have been few attempts to identify snowmelt timing using remote sensing. In this study, we describe the creation of snowmelt timing maps (STMs) which identify the day of year that each pixel of a remotely sensed image transitions from "snow-covered" to "no snow" during the spring melt season, controlling for cloud coverage and ephemeral spring snow storms. Derived from the 500 meter MODerate-resolution Imaging Spectroradiometer (MODIS) standard snow map, MOD10A2, this new dataset provides annual maps of snowmelt timing, with corresponding maps of cloud interference and interannual variability in snow coverage from 2001-2015. We first show that the STMs agree strongly with in-situ snow telemetry (SNOTEL) meteorological station measurements in terms of snowmelt timing. We then use the STMs to investigate the early snowmelt event of 2015 in the Cascade Mountains, USA, highlighting the protected areas of Mt. Rainier, Crater Lake, and Lassen Volcanic National Parks. In 2015 the Cascade Mountains experienced snowmelt 41 days earlier than the 2001-2015 average, with 25 percent of its land area melting more than 65 days earlier than average. The upper elevations of the Cascade Mountains experienced the greatest snowmelt anomaly. Our results are relevant to land managers and biologists as they plan adaptation strategies for mitigating the effects of climate change throughout temperate mountains.

snowmelt timing↗

DOE OSTI2024

Using Temporal Deep Learning Models to Estimate Daily Snow Water Equivalent Over the Rocky Mountains

Abstract In this study we construct and compare three different deep learning (DL) models for estimating daily snow water equivalent (SWE) from high‐resolution gridded meteorological fields over the Rocky Mountain region. To train the DL models, Snow Telemetry (SNOTEL) station‐based SWE observations are used as the prediction target. All DL models produce higher median Nash‐Sutcliffe Efficiency (NSE) values than a conceptual SWE model and interpolated gridded data sets, although mean squared errors also tend to be higher. Sensitivity of the SWE prediction to the model's input variables is analyzed using an explainable artificial intelligence (XAI) method, yielding insight into the physical relationships learned by the models. This method reveals the dominant role precipitation and temperature play in snowpack dynamics. In applying our models to estimate SWE throughout the Rocky Mountains, an extrapolation problem arises since the statistical properties of SWE (e.g., annual maximum) and geographical properties of individual grid points (e.g., elevation) differ from the training data. This problem is solved by normalizing the SWE with its historical maximum value to alleviate extrapolation for all tested DL models. Our work shows that the DL models are promising tools for estimating SWE, and sufficiently capture relevant physical relationships to make them useful for spatial and temporal extrapolation of SWE values.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Automated Cloud Based Long Short-Term Memory Neural Network Based SWE Prediction

Snow derived water is a critical component of the US water supply. Measurements of the Snow Water Equivalent (SWE) and associated predictions of peak SWE and snowmelt onset are essential inputs for water management efforts. This paper aims to develop an integrated framework for real-time data ingestion, estimation, prediction and visualization of SWE based on daily snow datasets. In particular, we develop a data-driven approach for estimating and predicting SWE dynamics using the Long Short-Term Memory neural network (LSTM) method. Our approach uses historical datasets (precipitation, air temperature, SWE, and snow thickness) collected at NRCS Snow Telemetry (SNOTEL) stations to train the LSTM network and current year data to predict SWE behavior. The performance of our prediction was compared for different prediction dates and prediction training datasets. Our results suggest that the proposed LSTM network can be an efficient tool for forecasting the SWE timeseries, as well as Peak SWE and snowmelt timing. Results showed that the window size impacts the model performance (where the Nash Sutcliffe efficiency (NSE) ranged from 0.96 to 0.85 and the Rooted Mean Square Error (RMSE) ranged from 0.038 to 0.07) with an optimum number that should be calibrated for different stations and climate conditions. In addition, by implementing the LSTM prediction capability in a cloud based site-monitoring platform, we automate model-data integration. By making the data accessible through a graphical web interface and an underlying API which exposes both training and prediction capabilities. The associated results can be made easily accessible to a broad range of stakeholders.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Evaluating 3 decades of precipitation in the Upper Colorado River basin from a high-resolution regional climate model

Abstract. Convection-permitting regional climate models (RCMs) have recently become tractable for applications at multi-decadal timescales. These types of models have tremendous utility for water resource studies, but better characterization of precipitation biases is needed, particularly for water-resource-critical mountain regions, where precipitation is highly variable in space, observations are sparse, and the societal water need is great. This study examines 34 years (1987–2020) of RCM precipitation from the Weather Research and Forecasting model (WRF; v3.8.1), using the Climate Forecast System Reanalysis (CFS; CFSv2) initial and lateral boundary conditions and a 1 km × 1 km innermost grid spacing. The RCM is centered over the Upper Colorado River basin, with a focus on the high-elevation, 750 km2 East River watershed (ERW), where a variety of high-impact scientific activities are currently ongoing. Precipitation is compared against point observations (Natural Resources Conservation Service Snow Telemetry or SNOTEL), gridded climate datasets (Newman, Livneh, and PRISM), and Bayesian reconstructions of watershed mean precipitation conditioned on streamflow and high-resolution snow remote-sensing products. We find that the cool-season precipitation percent error between WRF and 23 SNOTEL gauges has a low overall bias (x^ = 0.25 %, s = 13.63 %) and that WRF has a higher percent error during the warm season (x^ = 10.37 %, s = 12.79 %). Warm-season bias manifests as a high number of low-precipitation days, though the low-resolution or SNOTEL gauges limit some of the conclusions that can be drawn. Regional comparisons between WRF precipitation accumulation and three different gridded datasets show differences on the order of ± 20 %, particularly at the highest elevations and in keeping with findings from other studies. We find that WRF agrees slightly better with the Bayesian reconstruction of precipitation in the ERW compared to the gridded precipitation datasets, particularly when changing SNOTEL densities are taken into account. The conclusions are that the RCM reasonably captures orographic precipitation in this region and demonstrates that leveraging additional hydrologic information (streamflow and snow remote-sensing data) improves the ability to characterize biases in RCM precipitation fields. Error characteristics reported in this study are essential for leveraging the RCM model outputs for studies of past and future climates and water resource applications. The methods developed in this study can be applied to other watersheds and model configurations. Hourly 1 km × 1 km precipitation and other meteorological outputs from this dataset are publicly available and suitable for a wide variety of applications.

Geology↗

DOE OSTI2020

Hybrid predictive modeling approach simulated evapotranspiration and ecosystem respiration data

This dataset includes measured data used for developing hybrid-predictive-modeling (HPM) approach and simulated evapotranspiration and ecosystem respiration data across several Fluxnet sites, SNOTEL sites and East River locations (Chen et al., 2020 in review). Fluxnet sites considered in this study include: CA-OAS, CA-OBS, US-NR1, US-TON, US-VAR, US-SRM, US-WHS, US-WKG. Snotel sites Butte (ER-BT), Porphyry Creek (ER-PK) and Schofield Pass (ER-SP) are included as well 16 east river locations with different vegetation types denoted by evergreen forests (EF), deciduous forests (DF), riparian shrublands (RS) and meadow grassland (MS). This data package includes: 1) A summary spread sheet of the sites and acronyms for data variables. 2) Processed Fluxnet datasets that are obtained from Fluxnet.fluxdata.org with HPM estimated evapotranspiration (ET) and ecosystem respiration (RECO). Trained models and metrics for parameters for the Fluxnet sites were added in an update on November 2, 2020; 3) Processed SNOTEL datasets obtained from www.nrcs.usda.gov/snow/ and HPM estimated ET and RECO; 4) Weather data at east river locations obtained from DAYMET (daymet.ornl.gov) as well as HPM estimated ET and RECO. NDVI datasets are obtained from Landsat 5, 7 and 8 on Google Earth Engine platform. Detailed QA/QC and sampling methods of measured data are available at the original data source.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2015

Evaluation of the Snow Simulations from the Community Land Model, Version 4 (CLM4)

This paper evaluates the simulation of snow by the Community Land Model, version 4 (CLM4), the land model component of the Community Earth System Model, version 1.0.4 (CESM1.0.4). CLM4 was run in an offline mode forced with the corrected land-only replay of the Modern-Era Retrospective Analysis for Research and Applications (MERRA-Land) and the output was evaluated for the period from January 2001 to January 2011 over the Northern Hemisphere poleward of 30 deg N. Simulated snow-cover fraction (SCF), snow depth, and snow water equivalent (SWE) were compared against a set of observations including the Moderate Resolution Imaging Spectroradiometer (MODIS) SCF, the Interactive Multisensor Snow and Ice Mapping System (IMS) snow cover, the Canadian Meteorological Centre (CMC) daily snow analysis products, snow depth from the National Weather Service Cooperative Observer (COOP) program, and Snowpack Telemetry (SNOTEL) SWE observations. CLM4 SCF was converted into snow-cover extent (SCE) to compare with MODIS SCE. It showed good agreement, with a correlation coefficient of 0.91 and an average bias of -1.54 x 10(exp 2) sq km. Overall, CLM4 agreed well with IMS snow cover, with the percentage of correctly modeled snow-no snow being 94%. CLM4 snow depth and SWE agreed reasonably well with the CMC product, with the average bias (RMSE) of snow depth and SWE being 0.044m (0.19 m) and -0.010m (0.04 m), respectively. CLM4 underestimated SNOTEL SWE and COOP snow depth. This study demonstrates the need to improve the CLM4 snow estimates and constitutes a benchmark against which improvement of the model through data assimilation can be measured.

northern hemisphere↗

DOE OSTI2026

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗

DOE OSTI2026

Data from: "Reply to ‘The challenge of defining effectively-no-snow’"

This repository contains the data and code associated with the paper titled "Reply to ‘The challenge of defining effectively-no-snow’" published in Nature Reviews Earth and Environment, 2026. In this reply, we argue that the 10th percentile of peak SWE (Snow Water Equivalent), which we propose in the original article, can be used as intended given it's a standardized, impact-based benchmark for comparing snow conditions across regions, not as a literal measure of snow absence. We present new evidence with SNOwpack TELemetry (SNOTEL) data showing that years meeting the threshold are overwhelmingly associated with subsequent drought (given United States Drought Monitor conditions), supporting its hydrologic and societal relevance. We conclude that while the distinction between "effectively no snow" and "zero snow" should be clearly communicated, the original definition remains appropriate for assessing impacts on snow-dependent water systems. The file code_nree_ML_reply_2026.Rmd contains the main processing scripts which analyze the SNOTEL data. Data from the US Drought Monitor was downloaded at: https://usdmdataservices using the Get Drought Severity Statistics By Area Percent' option, saved to the *_HUC4_delineated.csv files (Hydrologic Unit Code), which are labeled accordingly. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > TERRESTRIAL HYDROSPHERE > SNOW/ICE↗

NASA NTRS1980

Evolution of snow sensors and the potential of microwave devices in operational telemetry networks

Both manually collected snow-course data and telemetered information from SNOTEL sites throughout the western United States are used to make forecasts of streamflow. Although these conventional methods have proven highly reliable over the years, they still exhibit some shortcomings. Active and passive microwave remote sensing systems recently developed and tests present the potential to eliminate some negative aspects associated with other types of sensors currently in use. Because of their ability to see through the snowpack they offer a unique opportunity to improve snowpack measurement techniques. As yet neither the active nor passive system is developed sufficiently to supplant existing methods in operational telemetry networks. Active microwave sensors are the most promising for the near future and are being used in research and development programs to study snow pillow performance in the SNOTEL system.

Shafer, B. A.↗

NASA NTRS2021

Modeling Snow Ablation over the Western United States Mountains: Patterns and Controlling Factors

When compared with differences in snow accumulation predicted by widely used hydrological models, there is a much greater divergence among otherwise “good” models in their simulation of the snow ablation process. Here, we explore differences in the performance of the Variable Infiltration Capacity model (VIC), Noah land surface model with multiparameterization options (Noah-MP), the Catchment model, and the third-generation Simplified Simple Biosphere model (SiB3) in their ability to reproduce observed snow water equivalent (SWE) during the ablation season at 10 Snowpack Telemetry (SNOTEL) stations over 1992–2012. During the ablation period, net radiation generally has stronger correlations with observed melt rates than does air temperature. Average ablation rates tend to be higher (in both model predictions and observations) at stations with a large accumulation of SWE. The differences in the dates of last snow between models and observations range from several days to approximately a month (on average 5.1 days earlier than in observations). If the surface cover in the models is changed from observed vegetation to bare soil in all of the models, only the melt rate of the VIC model increases. The differences in responses of models to canopy removal are directly related to snowpack energy inputs, which are further affected by different algorithms for surface albedo and energy allocation across the models. We also find that the melt rates become higher in VIC and lower in Noah-MP if the shrub/grass present at the observation sites is switched to trees.

Snow↗