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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A differential emissivity imaging technique for measuring hydrometeor mass and type

Abstract. The Differential Emissivity Imaging Disdrometer (DEID) is a new evaporation-based optical and thermal instrument designed to measure the mass, size, density and type of individual hydrometeors as well as their bulk properties. Hydrometeor spatial dimensions are measured on a heated metal plate using an infrared camera by exploiting the much higher thermal emissivity of water compared with metal. As a melted hydrometeor evaporates, its mass can be directly related to the loss of heat from the hotplate assuming energy conservation across the hydrometeor. The heat loss required to evaporate a hydrometeor is found to be independent of environmental conditions including ambient wind velocity, moisture level and temperature. The difference in heat loss for snow vs. rain for a given mass offers a method for discriminating precipitation phase. The DEID measures hydrometeors at sampling frequencies of up to 1 Hz with masses and effective diameters greater than 1 µg and 200 µm, respectively, determined by the size of the hotplate and the thermal camera specifications. Measurable snow water equivalent (SWE) precipitation rates range from 0.001 to 200 mm h−1, as validated against a standard weighing bucket. Preliminary field experiment measurements of snow and rain from the winters of 2019 and 2020 provided continuous automated measurements of precipitation rate, snow density and visibility. Measured hydrometeor size distributions agree well with canonical results described in the literature.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Seasonal forecasting skill for the High Mountain Asia region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorologically relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks to months is an area of active research and development. NASA's Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1–3-month lead times in the HMA region, including a portion of the Indian subcontinent, during the retrospective forecast period, 1981–2016. To assess forecast skill, we evaluate 2 m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to the similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from an anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system's ability to forecast HMA hydrometeorology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗
DOE OSTI2025

Canopy structure modulates the sensitivity of subalpine forest stands to interannual snowpack and precipitation variability

A declining spring snowpack is expected to have widespread effects on montane and subalpine forests in western North America and across the globe. The way that tree water demands respond to this change will have important impacts on forest health and downstream water subsidies. Here, we present data from a network of sap velocity sensors and xylem water isotope measurements from three common tree species (Picea engelmannii, Abies lasiocarpa and Populus tremuloides) across a hillslope transect in a subalpine watershed in the Upper Colorado River basin. We use these data to compare tree- and stand-level responses to the historically high spring snowpack but low summer rainfall of 2019 against the low spring snowpack but high summer rainfall amounts of 2021 and 2022. From the sap velocity data, we found that only 40 % of the trees showed an increase in cumulative transpiration in response to the large snowpack year (2019), illustrating the absence of a common response to interannual spring snowpack variability. The trees that increased water use during the year with the large spring snowpack were all found in dense canopy stands – irrespective of species – while trees in open-canopy stands were more reliant on summer rains and, thus, more active during the years with modest snow and higher summer rain amounts. Using the sap velocity data along with supporting measurements of soil moisture and snow depth, we propose three mechanisms that lead to stand density modulating the tree-level response to changing seasonality of precipitation: Topographically mediated convergence zones have consistent access to recharge from snowmelt which supports denser stands with high water demands that are more reliant and sensitive to changing snow. Interception of summer rain in dense stands reduces the throughfall of summer rain to surface soils, limiting the sensitivity of the dense stands to changes in summer rain. Shading in dense stands allows the snowpack to persist deeper into the growing season, providing high local reliance on snow during the fore-summer (early-summer) drought period. Combining data generated from natural gradients in stand density, like this experiment, with results from controlled forest-thinning experiments can be used to develop a better understanding of the responses of forested ecosystems to futures with reduced spring snowpack.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fine-scale vegetation composition and structure shape spatiotemporal variation in surface albedo across a low Arctic tundra landscape

The unprecedented rate of warming in the Arctic is driving changes in the structure and composition of tundra vegetation. Increases in deciduous tall shrub cover, height, and density are of particular concern, as these changes alter local surface albedo in ways that could amplify effects on the regional surface energy budget (SEB). Despite this importance, significant uncertainties remain in understanding the interplay between fine-scale vegetation patterns and emergent albedo dynamics across space and time. Here, we address these uncertainties by (1) quantifying spatiotemporal variation in surface shortwave albedo and (2) determining the relative influence of fine-scale vegetation composition, structure, and environmental conditions on albedo across a representative low-Arctic tundra landscape on Alaska’s Seward Peninsula. To do this, we synthesized multi-scale, multi-platform remote sensing observations, including a novel Landsat-derived albedo time series, a fine-scale map of Arctic plant functional type (PFT) fractional cover, and airborne LiDAR estimates of canopy height and topography. We show that there are substantial reductions in winter albedo for pixels dominated by tall, woody PFTs (28.13%) relative to pixels dominated by non-woody vegetation, but almost no change in summer albedo (3% increase). Further, we identified a unimodal trend in the relationship between canopy height and the timing of the springtime transition from high (snowy) to low (leafy) albedo (peak at 5.5 m), possibly because of competing ‘snow-fence’ and ‘protrusion’ snow-shrub interactions. To explore the primary drivers of albedo, we constructed a random forest model and found that canopy height and the fractional cover of woody PFTs were as- or more important predictors of winter albedo than topographic features. These findings provide strong evidence for the impacts of local vegetation characteristics on regional surface albedo, highlighting the need for better quantification of snow-shrub interactions to accurately predict the Arctic’s SEB under future environmental change.

Arctic↗
DOE OSTI2022

Tracking snowmelt during hydrological surface processes using a distributed hydrological model in a mesoscale basin on the Tibetan Plateau

We report that mountain snowpack is an important water resource for the high altitude and latitude regions where the terrain is complex. However, the snowmelt pathway and its actual contribution to streamflow and soil moisture are rarely reported and remain unclear in such regions. To fill in this knowledge gap, we incorporate a snowmelt pathway tracking algorithm to a high-resolution physics-based distributed-hydrology-soil-vegetation model (DHSVM), to track snowmelt movement and quantify snowmelt contributions in the surface hydrologic processes. A simple reservoir operation scheme is also incorporated in the model. The modified model is applied to a dammed meso-scale watershed in the northeast Tibetan Plateau, China to study the snow and reservoir effects. The results show that annual snow contribution to soil moisture (SC-SM) and snow contribution to streamflow (SC-S) significantly decrease over 1965-2019. At a monthly scale, SC-SM has the largest amplitude at the top soil layer and its peak in the deeper layer lags behind the upper layer, and mean monthly SC-S at all stations show bimodal distributions corresponding to snowfall season. Reservoir regulation has subtle impacts (≤2.0%) on SC-S. If the current climate change rate continues, monthly and annual streamflow at the outlet will decrease primarily due to snowpack reduction. To mitigate climate change impacts, better water resource management is needed in this watershed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of El Niño‐Southern Oscillation and Madden‐Julian Oscillation on the US Puget Sound Regional Hydroclimate

El Niño-Southern Oscillation (ENSO) and Madden-Julian Oscillation (MJO) are two major modes of climate variability with global hydroclimate impacts. However, their impacts often depend on the local climate and geography, resulting in large regional differences. In this study, we examined the connection of ENSO and MJO to the hydroclimate conditions and extremes in the Puget Sound (PS) basin located in the US Pacific Northwest coast. The results indicate that ENSO significantly modulates the cold season temperature and temperature-mediated hydrologic processes. El Niño cold seasons feature less snow accumulation and intensified surface runoff, even if the precipitation amount is similar to La Niña cold seasons. Therefore, El Niño causes more snow drought (in the form of compound dry and warm snow drought) and shifts the surface runoff seasonality by reducing runoff in the subsequent warm season. MJO phases 6–7 trigger more extreme precipitation, temperature, snowmelt, and runoff in the PS region at 0–9–day lags, and such connections are robust regardless of how the ENSO signals are removed. Meanwhile, MJO modulates large-scale extreme weather systems (e.g., atmospheric rivers) with significant enhancement during phases 6–7. ENSO impacts have intensified in the 2001–2020 period, whereas MJO impacts showed some phase shift in this period. This study reveals ENSO and MJO phases 6–7 as useful predictors of the PS hydroclimate anomalies/extremes at seasonal and daily scales, respectively. Utilizing these findings holds the potential to improve regional water resources prediction and management.

ENSO↗
DOE OSTI2022

Spatiotemporal evolution of melt ponds on Arctic sea ice

Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. 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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Lake-Effect Snowstorm Events and Associated Snowfall Totals Integrated from NOAA Storm Reports, ERA5, and HRRR for the Laurentian Great Lakes (1997–2024)

Lake-effect snowstorms are localized, impactful winter weather phenomena that can generate substantial snowfall totals and pose significant challenges for forecasting, transportation, and regional infrastructure. To support the analysis and modeling of these events, this dataset compiles observational reports of lake-effect snowstorms alongside corresponding snowfall estimates derived from gridded atmospheric datasets. The observational component of the data originates from the National Weather Service (NWS) winter storm report, subset to lake-effect snow event type, covering 1997–2024. For each lake-effect snow event, this data provides the impacted county, event start and end datetimes at an hourly resolution, as well as relevant storm narratives. The complementary reanalysis-derived data is sourced from European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5) and High-Resolution Rapid Refresh (HRRR) gridded data. For both gridded datasets, the maximum total snowfall (in units mm) was extracted, constrained by the county and datetimes specified by the observational report. ERA5 data covers the entire observational period (1997–2024), whereas HRRR data is only available from November 2016 – December 2024. Three CSV files are provided here: (1) the observational lake-effect snow event report, (2) ERA5 maximum snowfall detections for each event, and (3) HRRR maximum snowfall detections for each event. Relevant data from the observational files, such as impacted state and county, event datetimes, and event IDs, were included for convenience. Users can inspect and visualize the data using tools such as Microsoft Excel and Python pandas/matplotlib packages. This dataset may support a variety of applications, including climatological analyses of lake-effect snowfall, evaluation of snowfall representation in atmospheric datasets and numerical weather prediction models, and the development of machine learning approaches for detecting or predicting lake-effect snowfall events.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION > SOLID↗
DOE OSTI2020

Next-Generation Intensity-Duration-Frequency Curves for Climate-Resilient Infrastructure Design: Advances and Opportunities

National and international security communities (e.g., U.S. Department of Defense) have shown increasing attention for innovating critical infrastructure and installations due to recurring high-profile flooding events in recent years. The standard infrastructure design approach relies on local precipitation-based intensity-duration-frequency (PREC-IDF) curves that do not account for snow process and assume stationary climate, leading to high failure risk and increased maintenance costs. This paper reviews the recently developed next-generation IDF (NG-IDF) curves that explicitly account for the mechanisms of extreme water available for runoff including rainfall, snowmelt, and rain-on-snow under nonstationary climate. The NG-IDF curve is an enhancement to the PREC-IDF curve and provides a consistent design approach across rain- to snow-dominated regions, which can benefit engineers and planners responsible for designing climate-resilient facilities, federal emergency agencies responsible for the flood insurance program, and local jurisdictions responsible for developing design manuals and approving subsequent infrastructure designs. Further, we discuss the recent advances in climate and hydrologic science communities that have not been translated into actional information in the engineering community. To bridge the gap, we advocate that building climate-resilient infrastructure goes beyond the traditional local design scale where engineers rely on recipe-based methods only; the future hydrologic design is a multi-scale problem and requires closer collaboration between climate scientists, hydrologists, and civil engineers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Variability in observed stable water isotopes in snowpack across a mountainous watershed in Colorado

In this study, isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2020

Sensitivity evaluation of the Kudryavtsev permafrost model

Modeling is an important way to assess current and future permafrost spatial distribution and dynamics, especially in data poor areas like the Arctic region. Here, we evaluate a physics-based analytical model, Kudryavtsev's active layer model, which is widely used because it has relatively few data requirements. This model was recently incorporated into a component modeling toolbox, allowing for coupled modeling of permafrost and geomorphic processes over geological timescales. However, systematic quantitative assessment of the influence of its controlling parameters on permafrost temperature and active layer thickness predictions has not been undertaken before. In this work, we investigate the sensitivity of the Kudryavtsev's active layer model by Monte Carlo simulations to generate probability distributions for input parameters and compare predictions with a comprehensive benchmark dataset of in-situ permafrost observations over entire Alaska. Predicted permafrost surface temperature is highly dependent on mean annual air temperature (r = 0.78 on average), annual temperature amplitude (–0.41), and winter-averaged snow thickness (0.30). Uncertainty of predicted permafrost temperature is relatively small (RMSE = 1 °C), when air temperature and snow depth are well constrained. Similarly, RMSE between simulated and observed ALT at stations is ~0.08 m. However, under given air temperature and snow conditions, soil water content bias can significantly affect modeled active layer thickness (RMSE = 0.1 m or 40% of the observed active layer thickness). If soil water content has a large bias, improvements in other parameters may not significantly improve the active layer predictions of the Kudryavtsev's model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗
DOE OSTI2021

Geophysical Monitoring Shows that Spatial Heterogeneity in Thermohydrological Dynamics Reshapes a Transitional Permafrost System

Climate change is causing rapid changes of Arctic ecosystems. Yet, data needed to unravel complex subsurface processes are very rare. Using geophysical and in-situ sensing, this study closes an observational gap associated with thermohydrological dynamics in discontinuous permafrost systems. It highlights the impact of vegetation and snow thickness distribution on subsurface thermohydrological properties and processes. Large snow accumulation near tall shrubs insulates the ground and allows for rapid and downward heat flow. Thinner snowpack above graminoid results in surficial freezing and prevents water from infiltrating into the subsurface. Analyzing short term disturbances, we found that lateral flow could be a driving factor in talik formation. Inter-annual measurements show that deep permafrost temperatures increased by about 0.2°C over two years. The results, which suggest that snow-vegetation-subsurface processes are tightly coupled, will be useful for improving predictions of Arctic feedback to climate change, including how subsurface thermohydrology influences CO 2 and CH 4 fluxes.

58 GEOSCIENCES↗
DOE OSTI2021

Future Changes in Snowpack, Snowmelt, and Runoff Potential Extremes Over North America

Snowpack and snowmelt-driven extreme events (e.g., floods) have large societal consequences including infrastructure failures. However, it is not well understood how projected changes in the snow-related extremes differ across North America. Using dynamically downscaled regional climate model (RCM) simulations, we found that the magnitudes of extreme snow water equivalent, snowmelt, and runoff potential (RP; snowmelt plus precipitation) decrease by 72%, 73%, and 45%, respectively, over the continental United States and southern Canada but increase by up to 8%, 53%, and 41% in Alaska and northern Canada by the late 21st century. In California and the Pacific Northwest, there is a notable increase in extreme RP by 21% contrary to a decrease in snowmelt by 31% by the late century. These regions could be vulnerable to larger rain-on-snow floods in a warmer climate. Regions with a large variability among RCM ensembles are identified, which require further investigation to reduce the regional uncertainties.

54 ENVIRONMENTAL SCIENCES↗