Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snow contribution”

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 19 records
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.

The remoteness and extreme conditions of the Arctic make it a very difficult environment to investigate. In these polar regions covered by sea ice, the wind is relatively strong due to the absence of obstructions and redistributes a large part of the deposited snow mass, which complicates estimates for precipitation hardly distinguishable from blowing or drifting snow. Moreover, the snow mass balance in the sea ice system is still poorly understood, notably due to the complex structure of its surface. Quantitatively assessing the snow distribution on sea ice and its connection to the sea ice surface features is an important step to remove the snow mass balance uncertainties (i.e., snow transport contribution) in the Arctic environment. In this work we introduce snowBedFoam 1.0., a physics-based snow transport model implemented in the open-source fluid dynamics software OpenFOAM. We combine the numerical simulations with terrestrial laser scan observations of surface dynamics to simulate snow deposition in a MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate) sea ice domain with a complicated structure typical for pressure ridges. The results demonstrate that a large fraction of snow accumulates in their vicinity, which compares favorably against scanner measurements. However, the approximations imposed by the numerical framework, together with potential measurement errors (precipitation), give rise to quantitative inaccuracies, which should be addressed in future work. The modeling of snow distribution on sea ice should help to better constrain precipitation estimates and more generally assess and predict snow and ice dynamics in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Continuous snow depth and temperature measurements from dense network of above-ground distributed temperature profiling systems from 2021-09-23 to 2024-08-23, Seward Peninsula, Alaska

The dataset contains temperature measurements from distributed temperature profiling (DTP) systems (Dafflon et al., 2022; Wielandt et al., 2022; Wang et al., 2024a; Fiolleau et al., 2024) deployed vertically above the ground surface at a large number of locations from 2021 to 2024. The research is designed to improve understanding of the local heterogeneity in snow depth and snow thermal insulation dynamics, as well as their interactions in a discontinuous permafrost region (Wang et al., 2025). The DTP systems were deployed at 96 locations in a watershed along the Nome-Teller road at mile marker 27 (T27) and at 54 locations on a hillslope along the Kougarok road at mile marker 64 (K64) in the Seward Peninsula, Alaska. The probe location information is stored in Probe_locations_T27.csv and Probe_locations_K64.csv. Temperature measurements were recorded at 15-minute intervals using high-precision digital sensors (accuracy: ±0.1°C, resolution: 0.0078°C). The temperature probes, either 1.4 m or 1.6 m long, contain sensors spaced every 5 cm or 10 cm along their length. The temperature data are stored in compressed files following the format: DTP_snow_air_temperature_(site)_(start)_(end).zip, where site is either T27 or K64, and start and end represent the time series period. Within each ZIP file, individual CSV files are named by probe ID and contain temperature records at different heights above the ground surface.This dataset also includes derived snow depth time series over three snow seasons, estimated from temperature measurements. Snow depth was estimated by identifying the consecutive sensor pair that exhibited the largest drop in high-frequency temperature fluctuations (detailed in the methods). These data are stored in: Snow_depths_flags_(site)_(start)_(end).csv, which includes snow depth time series and corresponding quality flags (defined in the methods) from different probes. Additionally, the dataset includes derived metrics and supporting measurements at selected locations over two snow seasons, contributing to the manuscript of Wang et al., 2025. These locations were chosen based on the availability of high-quality snow depth time series during both seasons. The additional data include: (1) Air temperature proxies measured from the top sensors on the pole when they were not buried by snow, stored in Air_temperature_proxies_(site)_(start)_(end).csv (2) Ground interface temperature, recorded at 3 cm above the ground, stored in Ground_interface_temperature_(site)_(start)_(end).csv (3) Site characteristics, including vegetation height, elevation, and the topographic position index (TPI) within a 50 m radius, stored in Selected_probe_locations_gps_vegheight_tpi_elevation_(site).csv. These metrics were derived from 1 m resolution summer LiDAR-based digital elevation models and digital surface models from Singhania et al., 2023, DOI:10.5440/1832016. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv.This dataset is an updated version of a previous archive (Wang et al., 2024b, DOI: 10.15485/2475020), incorporating multiple seasons and improved snow depth estimation. Please note that due to large amount of information present in this dataset, many specificities associated with the acquisition of snow temperature, air temperature proxy and estimation of snow depth, and the future archiving of additional datasets on the soil temperature, thaw depth and soil characteristics at these locations, the author would welcome being contacted by people planning to use this dataset.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort 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 OSTI2023

Quantifying the Contributions of Aerosol- and Snow-Produced ClNO 2 through Observations and 1D Modeling

Nitryl chloride (ClNO 2 ) is a radical reservoir that forms and accumulates in the nocturnal atmospheric boundary layer influenced by combustion emissions and chloride (e.g., sea salt and road salt). Upon sunrise, ClNO 2 rapidly photolyzes to generate highly reactive chlorine radicals (Cl • ) that affect the air quality by generating secondary air pollutants. Recent studies have shown road salt aerosols and saline snowpack to be sources of ClNO 2 in the wintertime urban environment; however, the quantitative contributions of each chloride source are not known. In this study, we examine the vertically resolved contributions of aerosol particles and saline snowpack as sources of ClNO 2 by using an observationally constrained snow–atmosphere coupled one-dimensional model applied to wintertime Kalamazoo, Michigan, U.S. Model simulations show that ClNO 2 emitted from urban snowpack can be vertically transported throughout the entire atmospheric boundary layer and can be a significant source of ClNO 2 , contributing up to ~60% of the ClNO 2 budget near the surface. Here, modeled snowpack ClNO 2 emission rates were 6 (±7) times higher than the observationally derived emission rates, suggesting that not all snow chloride is available for reaction. ClNO 2 production from both aerosol particles and snow emissions are required to best simulate the observed surface-level ClNO 2 . Using the traditional bulk parameterization for ClNO 2 produced from particles significantly overestimated ClNO 2 due to the assumption of having equivalent dinitrogen pentoxide (N 2 O 5 ) uptake and chloride availability for the entire particle population. In comparison, the chemically resolved surface area-based parameterization slightly underestimated the observations and had uncertainties deriving from ClNO 2 production from residential wood burning particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗
DOE OSTI2022

Where does the dust deposited over the Sierra Nevada snow come from?

Abstract. Mineral dust contributes up to one-half of surface aerosol loading in spring over the southwestern United States, posing an environmental challenge that threatens human health and the ecosystem. Using self-organizing map (SOM) analysis with dust deposition and flux data from WRF-Chem and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), we identify four typical dust transport patterns across the Sierra Nevada, associated with the mesoscale winds, Sierra barrier jet (SBJ), North Pacific High (NPH), and long-range cross-Pacific westerlies, respectively. We find that dust emitted from the Central Valley is persistently transported eastward, while dust from the Mojave Desert and Great Basin influences the Sierra Nevada during mesoscale transport occurring mostly in winter and early spring. Asian dust reaching the mountain range comes either from the west through straight isobars (cross-Pacific transport) or from the north in the presence of the NPH. Extensive dust depositions are found on the west slope of the mountain, contributed by Central Valley emissions and cross-Pacific remote transport. In particular, the SBJ-related transport produces deposition through landfalling atmospheric rivers, whose frequency might increase in a warming climate.

Huang, Huilin (ORCID:0000000273286738)↗
DOE OSTI2023

Wind redistribution of snow impacts the Ka- and Ku-band radar signatures of Arctic sea ice

Wind-driven redistribution of snow on sea ice alters its topography and microstructure, yet the impact of these processes on radar signatures is poorly understood. Here, we examine the effects of snow redistribution over Arctic sea ice on radar waveforms and backscatter signatures obtained from a surface-based, fully polarimetric Ka- and Ku-band radar at incidence angles between 0° (nadir) and 50°. Two wind events in November 2019 during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition are evaluated. During both events, changes in Ka- and Ku-band radar waveforms and backscatter coefficients at nadir are observed, coincident with surface topography changes measured by a terrestrial laser scanner. At both frequencies, redistribution caused snow densification at the surface and the uppermost layers, increasing the scattering at the air–snow interface at nadir and its prevalence as the dominant radar scattering surface. The waveform data also detected the presence of previous air–snow interfaces, buried beneath newly deposited snow. The additional scattering from previous air–snow interfaces could therefore affect the range retrieved from Ka- and Ku-band satellite altimeters. With increasing incidence angles, the relative scattering contribution of the air–snow interface decreases, and the snow–sea ice interface scattering increases. Relative to pre-wind event conditions, azimuthally averaged backscatter at nadir during the wind events increases by up to 8 dB (Ka-band) and 5 dB (Ku-band). Results show substantial backscatter variability within the scan area at all incidence angles and polarizations, in response to increasing wind speed and changes in wind direction. Our results show that snow redistribution and wind compaction need to be accounted for to interpret airborne and satellite radar measurements of snow-covered sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Summer atmospheric circulation over Greenland in response to Arctic amplification and diminished spring snow cover

Abstract The exceptional atmospheric conditions that have accelerated Greenland Ice Sheet mass loss in recent decades have been repeatedly recognized as a possible dynamical response to Arctic amplification. Here, we present evidence of two potentially synergistic mechanisms linking high-latitude warming to the observed increase in Greenland blocking. Consistent with a prominent hypothesis associating Arctic amplification and persistent weather extremes, we show that the summer atmospheric circulation over the North Atlantic has become wavier and link this wavier flow to more prevalent Greenland blocking. While a concomitant decline in terrestrial snow cover has likely contributed to this mechanism by further amplifying warming at high latitudes, we also show that there is a direct stationary Rossby wave response to low spring North American snow cover that enforces an anomalous anticyclone over Greenland, thus helping to anchor the ridge over Greenland in this wavier atmospheric state.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Contribution of Tibetan Plateau ecosystems to local and remote precipitation through moisture recycling

The ecosystems of the Tibetan Plateau (TP) provide multiple important ecosystem services that benefit both local populations and those beyond, such as through climate regulation services on precipitation for East Asia and China. However, the precipitation regulation service of the TP ecosystems for supplying moisture and maintaining precipitation is yet to be evaluated. In this study, we used the moisture recycling framework and a moisture tracking model to quantify the precipitation regulation services of TP ecosystems for their contribution to precipitation. We found TP ecosystems contributed substantially to local and downwind precipitation, with a declined contribution from 221 mm/year for the TP and neighboring areas, to <10 mm/year for eastern China and other surrounding countries. Among ecosystem types, grassland contributed most to precipitation, followed by barren and snow lands, forests, and shrublands. In terms of seasonality, precipitation contribution from TP ecosystems was greater in summer months than in non-summer months for western China, while the opposite was true for eastern China – though the magnitude was much smaller. Over the past two decades, the significant evapotranspiration increases in TP translated to a widespread increase in precipitation contribution for TP and downwind beneficiary regions from 2000 to 2020. Our study provides a quantitative way to understand the precipitation regulation services of TP ecosystems through moisture recycling, substantiating their key role to safeguard precipitation and the water cycle for downwind regions– effectively acting as an ecological security shelter that could be perceived by the public.

Li, Yan↗
DOE OSTI2023

Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

With the highest albedo of the land surface, snow plays a vital role in Earth's surface energy budget and water cycle. Snow albedo is primarily controlled by snow grain properties (e.g., size and shape) and light-absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) Land Model (ELM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate, and Koch snowflake) and internal mixing of dust–snow, and it systematically evaluates the impacts on the surface energy budget and water cycle over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC–snow and dust–snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation (ELM_Control) with the default configurations of spherical snow grain shape, internal mixing of BC–snow, external mixing of dust–snow, and without TOP as well as the ELM simulation with new model features (ELM_New) can both capture the overall snow distribution reasonably. Additionally, ELM_New overall shows smaller biases in snow cover fraction than ELM_Control in spring when snowmelt is important for water management. The estimated LAP-induced RF in ELM_New ranges from 0 to 19.3 W m -2 with the area-weighted average value of 1.5 W m -2 that is comparable to the reported values in existing studies. The Koch snowflake shape, among other non-spherical shapes, shows the largest difference from the spherical shape in spring when snow processes related to the surface energy budget and water cycle have high importance. The impacts of the mixing state of LAP in snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP–snow can lead to larger snow albedo reduction and snowmelt, which further affect the surface energy budget and water cycle. The individual contributions of non-spherical snow shape, mixing state of LAP–snow, and local topography impacts on the snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and nonlinear interactions among the factors. Overall, the changes in net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W m -2 and -29.7 to 12.2 W m -2 , respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

58 GEOSCIENCES↗
DOE OSTI2022

Atmospheric particle abundance and sea salt aerosol observations in the springtime Arctic: a focus on blowing snow and leads

Abstract. Sea salt aerosols play an important role in the radiation budget and atmospheric composition over the Arctic, where the climate is rapidly changing. Previous observational studies have shown that Arctic sea ice leads are an important source of sea salt aerosols, and modeling efforts have also proposed blowing snow sublimation as a source. In this study, size-resolved atmospheric particle number concentrations and chemical composition were measured at the Arctic coastal tundra site of Utqiaġvik, Alaska, during spring (3 April–7 May 2016). Blowing snow conditions were observed during 25 % of the 5-week study period and were overpredicted by a commonly used blowing snow parameterization based solely on wind speed and temperature. Throughout the study, open leads were present locally. During periods when blowing snow was observed, significant increases in the number concentrations of 0.01–0.06 µm particles (factor of 6, on average) and 0.06–0.3 µm particles (67 %, on average) and a significant decrease (82 %, on average) in 1–4 µm particles were observed compared to low wind speed periods. These size distribution changes were likely caused by the generation of ultrafine particles from leads and/or blowing snow, with scavenging of supermicron particles by blowing snow. At elevated wind speeds, both submicron and supermicron sodium and chloride mass concentrations were enhanced, consistent with wind-dependent local sea salt aerosol production. At moderate wind speeds below the threshold for blowing snow as well as during observed blowing snow, individual sea spray aerosol particles were measured. These individual salt particles were enriched in calcium relative to sodium in seawater due to the binding of this divalent cation with organic matter in the sea surface microlayer and subsequent enrichment during seawater bubble bursting. The chemical composition of the surface snowpack also showed contributions from sea spray aerosol deposition. Overall, these results show the contribution of sea spray aerosol production from leads on both aerosols and the surface snowpack. Therefore, if blowing snow sublimation contributed to the observed sea salt aerosol, the snow being sublimated would have been impacted by sea spray aerosol deposition rather than upward brine migration through the snowpack. Sea spray aerosol production from leads is expected to increase, with thinning and fracturing of sea ice in the rapidly warming Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Shrubs Strongly Influence Snow Properties in Two Subarctic Watersheds

Understanding changes in snow distribution in permafrost ecosystems is fundamental to predicting their response to future climate change. The expansion of tall shrubs into tundra ecosystems can trap snow and insulate permafrost ecosystems during the winter, but the overall insulation effect is dependent upon many ecosystem properties. To study shrub–snow–ground interactions, small temperature sensors were deployed at two research sites on the Seward Peninsula of Alaska, USA, during the 2019–2020 winter. Snow temperatures were used to extrapolate multiple metrics, including freezing n-factors, the snow insulation effect, snow cover duration, and the length of the snowmelt period. Statistical and spatial analysis showed that shrub patches were a dominant control on all snow metrics. Within shrub patches, average ground temperatures were 2.1°C warmer, snow persisted 50 days longer, snow insulation was double, and a longer, later spring snowmelt period occurred compared to nonshrubby areas. Site-level differences contributed relatively little to variation in snow metrics, indicating that shrub presence is an overarching driver of snow–ground interactions at the locations examined. Shrub expansion, which is anticipated under climate change, will strongly impact future permafrost distribution and Arctic energy, water, and carbon cycles through snow–shrub–ground feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A cleaner snow future mitigates Northern Hemisphere snowpack loss from warming

Abstract Light-absorbing particles (LAP) deposited on seasonal snowpack can result in snow darkening, earlier snowmelt, and regional climate change. However, their future evolution and contributions to snowpack change relative to global warming remain unclear. Here, using Earth System Model simulations, we project significantly reduced black carbon deposition by 2081-2100, which reduces the December-May average LAP-induced radiative forcing in snow over the Northern Hemisphere from 1.3 Wm −2 during 1995-2014 to 0.65 (SSP126) and 0.49 (SSP585) Wm −2 . We quantify separately the contributions of climate change and LAP evolution on future snowpack and demonstrate that projected LAP changes in snow over the Tibetan Plateau will alleviate future snowpack loss due to climate change by 52.1 ± 8.0% and 8.0 ± 1.1% at the end of the century for the two scenarios, mainly due to reduced black carbon contamination. Our findings highlight a cleaner snow future and its benefits for future water supply from snowmelt especially under the sustainable development pathway of SSP126.

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 OSTI2026

Warming and snow loss increase reliance on old groundwater in a Colorado River headwater

Atmospheric warming is reducing snowpack, with uncertain effects on mountainous streamflow, a crucial water resource. Despite limited historical observations of groundwater–streamflow interactions above 2,500 m, new measurements in the Upper Colorado River headwaters indicate declining groundwater storage that is dated decades to millennia old. Here we use integrated hydrologic modelling spanning water years 2015–2021 to determine whether the loss of old-age groundwater buffers streamflow during low-snow years and whether that loss is exacerbated with warming. Results show that old-groundwater contributions to streams remain relatively steady through time, unlike the more variable contributions from young groundwater. Numerical experiments of increased surface air temperatures (+2.5 °C and +4 °C) increase rain–snow fractions and evapotranspiration and decrease runoff ratio by 2–3% per degree Celsius increase. As streamflow declines with warming, the age of groundwater supporting it gets older, in part owing to intermediate-aged (1–3 year) groundwater declining twice as fast. Simulations show that water table depths at higher elevations (>3,700 m) decline disproportionately and fail to recover even during wet years. These findings suggest altered groundwater–streamflow interactions with warming and snow loss, with implications for water resources.

Siirila-Woodburn, Erica R. [Lawrence Berkeley Nati↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Leveraging Groundwater Dynamics to Improve Predictions of Summer Low-Flow Discharges

Summer streamflow predictions are critical for managing water resources; however, warming-induced shifts from snow to rain regimes impact low-flow predictive models. Additionally, reductions in snowpack drive earlier peak flows and lower summer flows across the western United States increasing reliance on groundwater for maintaining summer streamflow. However, it remains poorly understood how groundwater contributions vary interannually. Here, we quantify recession limb groundwater (RLGW), defined as the proportional groundwater contribution to the stream during the period between peak stream flow and low flow, to predict summer low flows across three diverse western US watersheds. We ask (a) how do snow and rain dynamics influence interannual variations of RLGW contributions and summer low flows?; (b) which watershed attributes impact the effectiveness of RLGW as a predictor of summer low flows? Linear models reveal that RLGW is a strong predictor of low flows across all sites and drastically improves low-flow prediction compared to snow metrics at a rain-dominated site. Results suggest that strength of RLGW control on summer low flows may be mediated by subsurface storage. Subsurface storage can be divided into dynamic (i.e., variability saturated) and deep (i.e., permanently saturated) components, and we hypothesize that interannual variability in dynamic storage contribution to streamflow drives RLGW variability. In systems with a higher proportion of dynamic storage, RLGW is a better predictor of summer low flow because the stream is more responsive to dynamic storage contributions compared to deep-storage-dominated systems. Overall, including RLGW improved low-flow prediction across diverse watersheds.

54 ENVIRONMENTAL SCIENCES↗