Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Snow”

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characteristics and effects of aerosols during blowing snow events in the central Arctic

Sea salt aerosol (SSaer) significantly impacts aerosol-radiation and aerosol-cloud interactions, and sublimated blowing snow is hypothesized to be an important SSaer source in polar regions. Understanding blowing snow and other wind-sourced aerosols’ climate relevant properties is needed, especially during winter when Arctic amplification is greatest. However, most of our understanding of blowing snow SSaer comes from modeling studies, and direct observations are sparse. Additionally, SSaer can originate from multiple sources, making it difficult to disentangle emission processes. Here, we present comprehensive observations of wind-sourced aerosol during blowing snow events from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic. High wind speed strongly enhances total aerosol number, submicron sodium chloride mass, cloud condensation nuclei concentrations, and scattering coefficients. Generally, the relative response of aerosol properties to wind speed enhancement is strongest in fall when Arctic aerosol concentrations are lowest. Blowing snow events showed similar aerosol and environmental properties across events, apart from occasions with high snow age (>6 days since last snowfall). Coarse-mode number concentrations (>1 μm) are better explained by variability in wind speed averaged over 12-h air mass back trajectories arriving at the MOSAiC site compared to local, instantaneous wind speed, suggesting the importance of regional transport and consideration of air mass history for wind-driven aerosol production. These MOSAiC observations provide new insights into wind-driven aerosol in the central Arctic and may help validate modeling studies and improve model parameterizations particularly for aerosol direct and indirect radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Snow ALbedo eVOlution (SALVO) Campaign Qualitative Characterization of Snowpack Stratigraphy from April - June, 2024 in Utqiagivk, AK

The physical structure of the snowpack is impacted by and impacts the exchange of energy with the atmosphere. Qualitative characterization of snowpack stratigraphy in snow pits provides context for what physical processes are changing the snow (e.g., melt, wind-packing, hoar formation). 36 snowpits were dug as part of the 2024 Snow Albedo Evolution (SALVO) field campaign. Snow grain stratigraphy was recorded using standard International Association of Cryospheric Sciences (IACS) snow grain classification ( https://cryosphericsciences.org/publications/snow-classification/) and hand hardness tests (in order of least hard to hardest: F (fist), 4F (four fingers), 1F (one finger), P (pencil), K (knife), I (ice)). These snowpits were generally co-located with library site albedo measurements and are given the same index (Date, Site, Location, Position) when appropriate (see broadband and spectral albedo data archived from the SALVO II project).

Depth to bottom of snow layer↗
DOE OSTI2020

Dust dominates high-altitude snow darkening and melt over high-mountain Asia

Westerly-driven, long-range transportation of dust particles in elevated aerosol layers (EALs) is a persistent phenomenon during spring and summer over the Indian subcontinent. During the snow accumulation season, EALs transport substantial dust to the snow-covered slopes of high-mountain Asia (HMA). Here we use multiple satellite-based estimates to demonstrate a robust physical association between the EALs and dust-induced snow darkening over HMA. Results from a fully-coupled atmosphere-chemistry-snow model support these observations, revealing a signature of increasing dust-induced snow darkening with surface elevation over HMA that peaks near 4500 m. Moreover, the influence of dust on snow darkening is greater than that of black carbon above 4000 m. Our findings suggest a discernable role of dust in the observed spatial heterogeneity of snowmelt and snowline trends over HMA and highlight an increasing contribution of dust to snowmelt as the snowline rises with warming.

Elevated aerosol layers (EAL), Light absorbing par↗
DOE OSTI2022

Modeling Spatial Distribution of Snow Water Equivalent by Combining Meteorological and Satellite Data with Lidar Maps

Abstract An accurate characterization of the water content of snowpack, or snow water equivalent (SWE), is necessary to quantify water availability and constrain hydrologic and land surface models. Recently, airborne observations (e.g., lidar) have emerged as a promising method to accurately quantify SWE at high resolutions (scales of ∼100 m and finer). However, the frequency of these observations is very low, typically once or twice per season in the Rocky Mountains of Colorado. Here, we present a machine learning framework that is based on random forests to model temporally sparse lidar-derived SWE, enabling estimation of SWE at unmapped time points. We approximated the physical processes governing snow accumulation and melt as well as snow characteristics by obtaining 15 different variables from gridded estimates of precipitation, temperature, surface reflectance, elevation, and canopy. Results showed that, in the Rocky Mountains of Colorado, our framework is capable of modeling SWE with a higher accuracy when compared with estimates generated by the Snow Data Assimilation System (SNODAS). The mean value of the coefficient of determination R 2 using our approach was 0.57, and the root-mean-square error (RMSE) was 13 cm, which was a significant improvement over SNODAS (mean R 2 = 0.13; RMSE = 20 cm). We explored the relative importance of the input variables and observed that, at the spatial resolution of 800 m, meteorological variables are more important drivers of predictive accuracy than surface variables that characterize the properties of snow on the ground. This research provides a framework to expand the applicability of lidar-derived SWE to unmapped time points. Significance Statement Snowpack is the main source of freshwater for close to 2 billion people globally and needs to be estimated accurately. Mountainous snowpack is highly variable and is challenging to quantify. Recently, lidar technology has been employed to observe snow in great detail, but it is costly and can only be used sparingly. To counter that, we use machine learning to estimate snowpack when lidar data are not available. We approximate the processes that govern snowpack by incorporating meteorological and satellite data. We found that variables associated with precipitation and temperature have more predictive power than variables that characterize snowpack properties. Our work helps to improve snowpack estimation, which is critical for sustainable management of water resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Decline in Seasonal Snow during a Projected 20-Year Dry Spell

Snowpack loss in midlatitude mountains is ubiquitously projected by Earth system models, though the magnitudes, persistence, and time horizons of decline vary. Using daily downscaled hydroclimate and snow projections, we examine changes in snow seasonality across the U.S. Pacific Southwest region during a simulated severe 20-year dry spell in the 21st century (2051–2070) developed as part of the 4th California Climate Change Assessment to provide a “stress test” for water resources. Across California’s mountains, substantial declines (30–100% loss) in median peak annual snow water equivalent accompany changes in snow seasonality throughout the region compared to the historic period. We find that 80% of historic seasonal snowpacks transition to ephemeral conditions. Subsetting empirical-statistical wildfire projections for California by snow seasonality transition regions indicates a two-to-four-fold increase in the area burned, consistent with recent observations of high elevation wildfires following extended drought conditions. By analyzing six of the major California snow-fed river systems, we demonstrate snowpack reductions and seasonality transitions result in concomitant declines in annual runoff (47–58% of historical values). The negative impacts to statewide water supply reliability by the projected dry spell will likely be magnified by changes in snowpack seasonality and increased wildfire activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

[Retracted] Research on Application Experience Design of Ice and Snow Sports Equipment Based on Bee Colony Model

Sports equipment is the key to the smooth development of ice and snow sports. With the rapid development of social economy and the improvement of people’s living standards, the demand for ice and snow sports equipment is increasing day by day. This article presents an improved method based on the chaos theory and the bee colony algorithm to quantify the application experience design of ice and snow sports equipment and reduce the influence of uncertain factors on the design results. First, the chaos theory can establish the dataset of application experience design and analyze the discreteness of the set. According to the bee colony algorithm, the dataset is divided into several groups, and each group obtains the best application experience design by using the design optimization strategy. Finally, the results are mixed to obtain the final experience design results. Through MATLAB simulation analysis and verification, the improved bee colony model can improve the accuracy of application experience design of ice and snow sports equipment in an uncertain environment, shorten the overall design time, and meet the requirements of application experience design of different ice and snow sports equipment. Therefore, the model proposed in this paper is suitable for the application experience design of ice and snow sports equipment.

Li, Yuanjing (ORCID:000000018276647X)↗
DOE OSTI2024

LSTM-Based Data Integration to Improve Snow Water Equivalent Prediction and Diagnose Error Sources

Accurate prediction of snow water equivalent (SWE) can be valuable for water resource managers. Recently, deep learning methods such as long short-term memory (LSTM) have exhibited high accuracy in simulating hydrologic variables and can integrate lagged observations to improve prediction, but their benefits were not clear for SWE simulations. Here we tested an LSTM network with data integration (DI) for SWE in the western United States to integrate 30-day-lagged or 7-day-lagged observations of either SWE or satellite-observed snow cover fraction (SCF) to improve future predictions. SCF proved beneficial only for shallow-snow sites during snowmelt, while lagged SWE integration significantly improved prediction accuracy for both shallow- and deep-snow sites. The median Nash–Sutcliffe model efficiency coefficient (NSE) in temporal testing improved from 0.92 to 0.97 with 30-day-lagged SWE integration, and root-mean-square error (RMSE) and the difference between estimated and observed peak SWE values d max were reduced by 41% and 57%, respectively. DI effectively mitigated accumulated model and forcing errors that would otherwise be persistent. Moreover, by applying DI to different observations (30-day-lagged, 7-day-lagged), we revealed the spatial distribution of errors with different persistent lengths. For example, integrating 30-day-lagged SWE was ineffective for ephemeral snow sites in the southwestern United States, but significantly reduced monthly-scale biases for regions with stable seasonal snowpack such as high-elevation sites in California. These biases are likely attributable to large interannual variability in snowfall or site-specific snow redistribution patterns that can accumulate to impactful levels over time for nonephemeral sites. These results set up benchmark levels and provide guidance for future model improvement strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Determination of Ground Subsidence Around Snow Fences in the Arctic Region

In this study, we analyzed the effects of snow cover changes caused by snow fences (SFs) installed in 2017 in the Alaskan tundra to examine ground subsidence. Digital surface model data obtained through LiDAR-based remote sensing in 2019 and 2022, combined with a field survey in 2021, revealed approximately 0.2 m of ground subsidence around the SF. To investigate the relationship between SF-induced snow cover changes and ground subsidence, geophysical methods, electrical resistivity tomography (ERT) and ground-penetrating radar (GPR), were applied in 2023 to analyze subsurface characteristics. The increased snow cover due to the SF-enhanced insulation, delaying the penetration of winter cold into the subsurface. This delay caused subsurface temperatures to decrease more slowly, melting the upper permafrost and increasing the thickness of the active layer. ERT and GPR surveys well delineated the boundary between the active layer and permafrost, confirming that the increased snow cover thickened the active layer. This thickening led to the melting of pore ice, causing water runoff and ground compaction, which resulted in subsidence. The runoff also formed channels flowing eastward over the SF. This study highlights how changes in snow cover can influence active layer properties, leading to localized environmental changes and ground subsidence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The impact of forest-controlled snow variability on late-season streamflow varies by climatic region and forest structure

Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Snow-corrected vegetation indices for improved gross primary productivity assessment in North American evergreen forests

North American evergreen forests cover large areas and influence the global carbon cycle. Satellite remote sensing has been used to track the phenology of ecosystem photosynthesis of these forests by detecting variation in vegetation optical properties associated with physiological and structural features, and most of these methods have been closely tied to vegetation greenness. However, in evergreens, the application of satellite data to monitor photosynthetic phenology is often limited by the lack of sensitivity of greenness-based indices. In this study, we identified 47 evergreen forest flux sites in North America that had MODIS observation overlapping with the flux tower records. We then calculated four vegetation indices using MODIS MAIAC data (MCD19A1), including NDVI, CCI, NIRv, and kNDVI, for the 47 flux sites and evaluated relationships between gross primary productivity (GPP) and vegetation indices across the North American evergreen forests. Our results showed that snow had substantial effects on the performance of all vegetation indices in tracking GPP phenology, particularly in the early spring when rapid changes occurred to both GPP and snow cover. Furthermore, different vegetation indices were affected differently, indicating contradictory and confounding effects of snow on these indices. After correcting for the snow effects, both CCI and NIRv performed well in tracking GPP phenology, albeit for different reasons. CCI is sensitive to seasonal changes in the relative levels of chlorophyll and carotenoid pigments, which are closely tied to GPP phenology in evergreens. NIRv is sensitive to the absorbed photosynthetically active radiation and to the contribution of deciduous components to the overall optical properties. We also found that correlations between GPP and vegetation indices varied among ecoregions and climate classes. In general, regions with pronounced seasonal GPP patterns had stronger correlations between GPP and greenness-based indices than regions with weaker seasonal GPP patterns. These biome differences were less pronounced for CCI. The snow artifacts and complementary vegetation index effects reported here should be considered in any large-scale studies of GPP using reflectance-based indices from optical satellites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Identifying snowfall elevation patterns by assimilating satellite-based snow depth retrievals

Precipitation in mountain regions is highly variable and poorly measured, posing important challenges to water resource management. Traditional methods to estimate precipitation include in-situ gauges, Doppler weather radars, satellite radars and radiometers, numerical modeling and reanalysis products. Each of these methods is unable to adequately capture complex orographic precipitation. Here, we propose a novel approach to characterize orographic snowfall over mountain regions. We use a particle batch smoother to leverage satellite information from Sentinel-1 derived snow depth retrievals and to correct various gridded precipitation products. This novel approach is tested using a simple snow model for an alpine basin located in Trentino Alto Adige, Italy. Here, we quantify the precipitation biases across the basin and found that the assimilation method (i) corrects for snowfall biases and uncertainties, (ii) leads to cumulative snowfall elevation patterns that are consistent across precipitation products, and (iii) results in overall improved basin-wide snow variables (snow depth and snow cover area) and basin streamflow estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

iButton snow-ground interface temperature measurements in Los Alamos, New Mexico from 2023-2024

Snow/ground interface temperature measurements were collected at two sites in Los Alamos, New Mexico. Data were collected from November 29, 2023 to April 8, 2024 using iButton Link DS1921G-F5# Thermochron miniature temperature sensors (https://www.ibuttonlink.com/products/ds1921g). These sensors are a cost-efficient way to collect snowpack temperatures at a higher spatial resolution than what is normally achieved. iButton data were collected every 3 hours from a total of 19 iButtons. iButtons were placed in pairs, with one iButton placed at the ground surface and another buried 1 - 5 cm below the ground surface. One buried iButton did not successfully collect data, and therefore was excluded from this dataset. Data were collected throughout the snow cover season so that snowpack characteristics could be derived using the temperature data. Specifically, this dataset was used as a validation source for a novel machine learning approach to estimating snow depth (see related publication). Sensors were placed in areas with bare ground or minimal grass coverage, located away from any large vegetation. At Site A (TA51), manual snow depths were collected as validation data. These measurements were taken next to iButtons periodically throughout the winter, and notes on other precipitation types were also recorded. At Site B (TA6 Meteorological Station), a nearby sensor collected snow depths throughout the winter. This dataset contains one *.csv file of snow/ground interface temperatures at two sites, one *.csv file of manually collected snow depths, and one *.kml file of sensor locations. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using Machine Learning Algorithm to Detect Blowing Snow and Fog in Antarctica Based on Ceilometer and Surface Meteorology Systems

Blowing snow is a common weather phenomenon in Antarctica and plays an important role in the water vapor cycle and ice sheet mass balance. Although it has a significant impact on the climate of Antarctica, people do not know much about this process. Fog events are difficult to distinguish from blowing snow events using existing detection algorithms by a ceilometer. In this study, based on ceilometer, the meteorological parameters observed by surface meteorology systems are further combined to detect blowing snow and fog using the AdaBoost algorithm. The weather phenomena recorded by human observers are ‘true’. The dataset is collected from 1 January 2016 to 31 December 2016 at the AWARE site. Among them, three-quarters of the data are used as the training set and the rest of the data as the testing set. The classification accuracy of the proposed algorithm for the testing set is about 94%. Compared with the Loeb method, the proposed algorithm can detect 89.12% of blowing snow events and 76.10% of fog events, while the Loeb method can only identify 64.29% of blowing snow events and 31.87% of fog events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Radar-based Snow Property Profile Retrievals for NSA

The product contains estimated snow microphysical properties retrieved primarily from vertically-profiling radar observations using a priori constraints that are dependent on atmospheric state. The product provides data at 1-minute time resolution, including snowfall rate at the surface and vertical profiles at 300-meter resolution of snow particle size distribution parameters and snow water contents. Estimated uncertainties for the snowfall rates and snow water contents are also provided. The product is obtained using an optimal estimation retrieval in which radar forward model and a priori state characteristics are based on expected snow particle microphysical and radar scattering properties evaluated from field experiment data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impact of Seasonal Snow-Cover Change on the Observed and Simulated State of the Atmospheric Boundary Layer in a High-Altitude Mountain Valley

The structure and evolution of the atmospheric boundary layer (ABL) under clear-sky fair weather conditions over mountainous terrain is dominated by the diurnal cycle of the surface energy balance and thus strongly depends on surface snow cover. We use data from three passive ground-based infrared spectrometers deployed in the East River Valley in Colorado's Rocky Mountains to investigate the response of the thermal ABL structure to changes in surface energy balance during the seasonal transition from low to high snow cover. Temperature profiles were retrieved from the infrared radiances using the optimal estimation physical retrieval Tropospheric Remotely Observed Profiling via Optimal Estimation. A nocturnal surface inversion formed in the valley during clear-sky days, which was subsequently mixed out during daytime with the development of a convective boundary layer when snow cover was low. Over high snow cover, a very shallow convective boundary layer formed, above which the inversion persisted through the daytime hours. We compare these observations to NOAA's operational High-Resolution-Rapid-Refresh model and find large warm biases on clear-sky days resulting from the model's inability to form strong nocturnal inversions and to maintain the stable stratification in the valley during daytime when there was snow on the ground. We suggest several factors contributing to the large model errors. These are (a) the inability of the model to represent well-developed thermally driven flows likely due to the too coarse horizontal grid spacing (3 km), (b) too much convective mixing during daytime, and (c) too strong vertical coupling between the valley atmosphere and the free troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Divergent responses of historic rain-on-snow flood extremes to a warmer climate

Global warming is altering flood risks induced by rain-on-snow events. However, decision-makers lack guidance on how rain-on-snow induced extreme floods could be altered with warming. Here, storyline analyses using a kilometer-scale land surface model reveal diverse responses of four historically-impactful, decision-relevant rain-on-snow induced extreme flood events over the contiguous U.S. to warming, due to alterations in their water budgets. For the 2017-Feb California floods, runoff first increases and then decreases with warming, peaking under the +3 K scenario, while runoff of the 2017-Jan California floods increases monotonically by ~53%/K. Contrastingly, runoff of the 1996-Jan Mid-Atlantic floods decreases gradually with warming. Despite these differences, warming generally shifts flood-generating regimes along elevation profiles. High elevations could experience notably increased runoff, while low elevations encounter a shift from rain-on-snow-driven to rainfall-dominated runoff. These findings underscore the need for flood control planning to quantify region- and elevation-specific changes in rain-on-snow events in a warmer climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Code Description for "Brief Communication: Monitoring snow depth using small, cheap, and easy-to-deploy ground surface temperature sensors"

Temporally continuous snow depth estimates are vital for understanding changing snow patterns and impacts on permafrost in the Arctic. We train a random forest machine learning model to predict snow depth from variability in ground surface temperature. To our knowledge, this is the first time that small ground surface temperature sensors have been used to estimate snow depth. The model performs well at sites where the model was trained and at pan-arctic evaluation sites (RMSE <= 0.15 m). Small temperature sensors are cheap and easy-to-deploy, so this technique enables spatially distributed and temporally continuous snowpack monitoring to an extent previously infeasible. The model is flexible and can be applied to datasets retroactively to retrieve snow depth estimates at additional sites. This code package includes a *.joblib file of the trained random forest model and a *.ipynb file showing how to clean input data, train the random forest model, and apply the model.

Bachand, Claire↗