Search NASA⌕ Search

SEARCH · Search NASA

Results for “Snowmelt”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Watershed response of the Feather River Basin, California, United States of America to future climate changes

Increasing mean annual temperatures under climate change are expected to reduce seasonal snowpack, increase evapotranspiration (ET), and alter summer baseflow in headwater watersheds worldwide. Strong regional variability in hydrologic responses highlights the need for catchment-scale, physically based models to assess future flood risk and water availability. This study examines climate-driven changes in the hydrologic response of the Upper Feather River watershed in the Sierra Nevada Mountains, California. Four priority climate models and two representative concentration pathways (RCP4.5 and RCP8.5) are used to evaluate future hydrologic responses of the watershed. Hydrologic processes are simulated using an objectively calibrated Soil and Water Assessment Tool Plus model for a historical baseline (1986–2005) and a future period 2070–2099), driven by observed and projected precipitation and temperature. Results indicate a declining contribution of snowfall to annual precipitation, with peak snowfall and water yield shifting 1–3 months earlier. Long-term annual maximum flows are projected to increase considerably, whereas low-flow responses are mixed, with both increases and decreases projected by the end of the century. These findings highlight the need for adaptive watershed management to enhance flood protection, water storage, and drought resilience. Future water resource planning should also account for one-to-three-month shifts in peak water yield and surface runoff due to changes in snowmelt timing and a lower snowfall-to-rainfall ratio under climate change. The modeling framework and insights are transferable to other snow-dominated headwater watersheds experiencing climate-driven change.

Tigabu, T [UC Davis]↗
DOE OSTI2020

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Science-integrated Artificial-intelligence for Flooding and precipitation Extremes (SAFE)

A grand challenge in hydrologic science is to understand why signals of climate change and variability, which are often visible in precipitation extremes at aggregate scales, are not consistently observed in the case of extreme flooding. However, a solution to this challenge may prove elusive unless the water cycle is viewed in an integrative manner. Thus, for riverine flooding, while Hortonian (infiltration excess) runoff may have stronger correlation with precipitation extremes and hence perhaps to warming trends or climate oscillators, Dunne (saturation excess) runoff may have a more complex relationships with time series of precipitation and with evaporation and transpiration, but rain-on-snow and snowmelt events may depend on land-surface and atmospheric temperatures. Atmospheric rivers and tropical cyclones lead to precipitation or flooding and are impacted by climate. Flooding assessments need to consider long-term baselines, evolving risk factors, coupled natural-human systems, and novel adaptation such as nature-inspired design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Space and Time Dynamics of Transpiration in the East River Watershed: Biotic and Abiotic Controls

Alpine forests have important impacts on water resources by affecting how much of the precipitation that falls on a watershed is routed to streamflow or returned to the atmosphere as transpiration – so-called “green water”. In many of the high-altitude watersheds that supply critical water resources to the western US, it remains a challenge to predict how much water forests will utilize limiting the ability to predict both short and long term information on downstream water resources. The uncertainty emerges because transpiration varies between tree species, across landscapes and over time in ways that cannot readily be predicted from physically-based models for evaporation. In this project, we utilized a technique called sap flux to measure the rate of water use for three common species of fir, spruce and aspen distributed across a hillslope in the East River Watershed from 2019-2021. We also made ~weekly measurements of the concentration of stable isotopic tracers in transpiration to understand not only how much water was being used but whether it originated from summer rain or snowmelt. Firstly, our results show that the species tend to use similar amounts of total water but they achieve this similar cumulative water flux both by using water at different times of the year and by using different water sources. This shows the importance of including species-level information into models used to make seasonal streamflow predictions. Secondly, our results show that trees lower on the hillslope use older waters all through the summer whereas trees higher on the hillslope rely more heavily on recent summer rain. This means trees lower on the hillslope are less sensitive to year-to-year changes in summer climate. Lastly, we found that sites with dense forest stands displayed the largest amounts of year-to-year changes in transpiration and thus these high localized spots in the watershed are the primary drivers of changes in “green water” use. The results provide an unprecedented spatial and temporal picture of forest water use that confirms a number of fundamental hypotheses on hillslope ecohydrology that had not previously been tested against observations. Predictions of streamflow that will benefit from this work are highly valuable for agriculture management, flood control and ecological restoration efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Science Area 1: Standard Award: Model-Data Fusion to Examine Multiscale Dynamical Controls on Snow Cover and Critical Zone Moisture Inputs (Final Report)

In many mountain watersheds of the world, seasonal snowpacks play an important role as natural reservoirs of water. Seasonal snowpacks accumulate water during cold, wet winter months that subsequently melts. Downstream communities depend on water from melting seasonal snowpacks to support agricultural, industrial, and municipal water needs. Rapidly melting snowpacks can also present a flooding hazard, particularly if snowpacks melt at rates faster than anticipated and where adequate reservoir capacity is unavailable to buffer river flows associated with melt. The spatial and temporal dynamics of snow accumulation and melt also play an important role in supporting upland ecosystems in mountain landscapes. Snowmelt provides soil moisture that enable terrestrial ecosystem productivity and exert control on soil microorganisms that play important roles in global carbon cycles. Climate warming is gradually decreasing the amount of precipitation in mountain watersheds arriving as snow, presenting potentially profound disruptions to mountain ecosystems, as well as downstream delivery of water. The overarching goal of this project was to understand how interactions between the near-surface atmosphere and surface topography control the input, accumulation, retention, and release of water from mountain snowpacks. Over a 5-year period, this project pursued an approach combining high-resolution regional climate modeling, satellite and airborne remote sensing data, and ground-based observations to develop and analyze benchmark datasets to address overarching science questions and hypotheses. Key products include a continuous, long-term, high spatiotemporal resolution (1 km/1 hr) dataset characterizing key climate variables in the Upper Colorado River Basin. The dataset included historical estimates of precipitation, temperature, humidity, solar and longwave radiation, and wind speeds. Additionally, the project developed a 20+ year long, 30 m spatial, daily temporal multi-sensor dataset characterizing snow presence/absence in the East/Taylor River watersheds in the Central Rocky Mountains of Colorado. The project supported training of 1 postdoctoral scholar, 1 Ph.D. student, and 1 M.S. student.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

NSA Site Science: Use of ARM Observations from Northern Alaska to Evaluate and Improve Prediction Capabilities

The Arctic is warming at a rate nearly double that of the rest of the planet, leading to profound changes in atmospheric, oceanic, and ice processes. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility has played a significant role in Arctic research, operating observatories in Alaska's North Slope for over 25 years. These observatories provide a rich, and wide-reaching dataset that offers insight into atmospheric processes in northern Alaska. This report details the results of a nine-year research project (2015–2024) supported by the DOE Atmospheric Systems Research (ASR) program, that leverages data from ARM’s deployment of observing facilities at Utqiaġvik (known as the North Slope of Alaska, or NSA, site) and Oliktok Point, Alaska. The project was conducted in two phases: - Phase 1 (2015–2019): Focused on understanding key atmospheric processes at Oliktok Point, including cloud formation, high-latitude precipitation, aerosol-cloud interactions, and cloud properties. - Phase 2 (2019–2024): Extended the research to the broader North Slope region, using data from both Oliktok Point and Utqiaġvik. Topics explored included surface energy budgets, atmospheric stability, ice nucleation processes, and microphysics in Arctic clouds. The project resulted in numerous research products, including 50 peer-reviewed publications and dissertations, 169 presentations, and 10 data products. These products cover a variety of topics, including: - Cloud Macro- and Microphysical Properties: Arctic clouds play a crucial role in energy transfer, and accurate representation in models is critical. The study explored cloud transitions, ice crystal shapes, and dual-wavelength radar data to understand ice crystal habits and size distributions. - Aerosol Properties and Processes: The team examined aerosol sources in the Arctic, including industrial emissions and natural sources. Observations showed significant spatial gradients in aerosol concentrations due to human activities and wildfire smoke. The influence of aerosols on cloud formation and the surface energy budget was also assessed. - Aerosol-Cloud Interactions: Research revealed that aerosols might suppress cloud ice production, affecting cloud radiative forcing and precipitation. The impact of local industrial emissions on cloud properties was also investigated. - Contextualizing the North Slope of Alaska in the context of the broader Arctic: To understand broader trends, the project evaluated large-scale circulation patterns and the influence of weather systems on the Arctic. Studies indicated that large-scale processes play a significant role in temperature patterns and the timing of snowmelt. - Advancing ARM Observational and Modeling Capabilities: The project developed new radar data products and advanced measurement techniques, including clutter mitigation and drizzle detection. Uncrewed aerial systems (UAS) and tethered balloon systems (TBS) were deployed to gather detailed atmospheric data. Additionally, the project supported 10 early career scientists, providing training and mentorship to undergraduate interns, graduate students, postdoctoral researchers, and early career researchers. These efforts contributed to the advancement of ARM research capabilities and fostered a new generation of scientists skilled in Arctic atmospheric research. Ultimately, this ASR-supported project has provided valuable insights into Arctic atmospheric processes and their broader climate implications. Recommendations for future work include continuing support for long-term observing at Arctic locations to foster additional research, further exploration of aerosol-cloud interactions and the potential impacts of enhanced industrialization of the Arctic, and expanded use of uncrewed systems to gather data in this remote and harsh environment. Additionally, the data products developed by this work, and the data products developed through the ARM infrastructure, leave a treasure-trove of additional information that should be explored for many years to come to gain additional insight into physical processes in the Arctic atmosphere that drive the rapid changes occurring in at high latitudes and their global impact.

58 GEOSCIENCES↗
DOE OSTI2024

The Relationship of Hydrology in the Green River to Colorado Pikeminnow Populations

The Federally listed Colorado Pikeminnow (CP) (Ptychocheilus Lucius) has been declining in the Green River Basin over the last 20 years. Muth et al. (2000) recommended flows and temperatures to benefit endangered fish species downstream of Flaming Gorge Dam on the Green River. The recommendations identified a range of flows during the base flow period (flows that occur after the annual snowmelt runoff period, usually from early summer to the following spring) that were intended to benefit CP. The recommendations indicated that the mean flow for the summer–winter period should be established each year on the basis of anticipated hydrologic conditions, but adjustments could be made if hydrologic conditions change. Flow should gradually decline from peak flow to base flow, with the base flow reached by early to middle summer (depending on hydrologic conditions) and maintained through February. The rate of decline should depend on rates of decline in Yampa River flows and the recommended rates of decline for dam releases. No specific recommendations were made for base flow variation in Reach 1, but variation around the annual mean base flow should be restricted to achieve recommended levels of variation in Reach 2 (middle Green River) (Table 5.5 in Muth et al. 2000). Implementation of the flow and temperature recommendations served as the proposed action for an EIS on Flaming Gorge Dam operations. A Biological Assessment of the impacts on endangered species of implementing the recommendations was prepared and a Biological Opinion that called for development of a study plan to examine the effects of implementation. Projects to address this evaluation were identified in the Green River Study Plan and have been conducted since then. The purpose of the Study Plan was to identify and recommend to the Recovery Program those monitoring or research projects necessary for implementation and evaluation of the flow and temperature recommendations. Those projects included studies to evaluate the anticipated effects of implementing the recommendations (including potential adverse effects identified in the 2006 Biological Opinion) and studies to examine recognized uncertainties of the recommendations. A new recovery program summer base flow study plan was approved in September of 2020 to evaluate stable and elevated summer base flows below Flaming Gorge Dam intended to benefit CP spawning and rearing success. Anecdotal observations suggest CP production may have been higher in years prior to flow restrictions for endangered fish. Based on long-term CP monitoring data in the Green River, Bestgen and Hill (2016) concluded that CP would benefit from higher summer baseflows (1,700-3,000 ft 3 /sec) with flows above or below associated with lower abundance of age-0 CP. Therefore, flexibility in the ROD is currently being utilized to establish higher stable summer base flows in an attempt create more advantageous conditions for spawning and rearing CP. These fish recovery flows are restrictive to WAPA in the summer months, therefore it is necessary to conduct an analysis of historic river data relationships to CP abundance to better understand if the recommended fish flows (elevated and stable base flows) are necessary for fish recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

How Snow Drives the Seasonal Evolution of Land and Sea Surface Albedos in the Alaskan High Arctic: Final Technical Report

The purpose of the project was to observe and quantify temporal variation in snow albedo and snow characteristics across the Arctic coastal landscape as winter transitioned into spring and snowmelt occurred, both on tundra and sea ice. This transition is bounded by fully snow-covered landscapes with broadband albedos of approximately 0.8 and snow-free landscapes with albedos of 0.15 (tundra or ponded sea ice). For these landscapes, we monitored the spring surface characteristics and radiative properties nearly daily at three locations on or near the Department of Energy ARM North Slope of Alaska (NSA) User Facility in Utqiagvik, Alaska: the central NSA Facility (hereafter called ARM), one near NSA-E12 (BEO), and one on the sea ice of Elson Lagoon (ICE) during three melt seasons (2019, 2022, and 2024). The field campaign component of this award was named SALVO ( S now AL bedo E VO lution). Typical field seasons began in mid-April and lasted until mid-June. Main measurements included snow depth (at 1-m intervals), broadband albedo (at 5-m intervals), spectral albedo (at 5- m intervals), and multiple digital images. Orthomosaics were converted into binary images to determine the snow-covered fraction over time across various landscapes. Additional measurements included basic weather data and snow-ground (or ice) interface temperatures. Sky conditions were observed and photographed to help assess albedo values.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparative Study of Hydrology and Icemelt in Three Nepal River Basins Using the Glacio-Hydrological Degree-Day Model (GDM) and Observations From the Advanced Scatterometer (ASCAT)

An assessment of the water supply and its seasonal and annual changes over the century in the High Mountain Asia (HMA) region is of increasing interest due to its potential impact on one-sixth of the global population. In order to understand the changing hydrology and snow and ice melt, we used remotely sensed Advanced Scatterometer (ASCAT) observations of glacier melt (GM) and a distributed and gridded Glacio-hydrological Degree-day Model (GDM) in three river basins: Tamor, Trishuli and Marsyangdi. The GDM-estimated contribution of snowmelt, icemelt, rainfall and baseflow in river flows is found to be most accurate in the Trishuli River basin, with Nash-Sutcliffe efficiency (NSE) between the estimated and observed discharges of 0.81 and volume differences of -0.5%, and reasonably accurate in the Tamor River basin, with NSE of 0.69 and volume difference of -7.51%. Similarly, NSE of 0.81 and volume difference of 4.64% in Marsyangdi River basin. We find strong similarities in the timing of glacier melting using the GDM and from observations from the ASCAT GM, determining the seasonal start of glacier melting to within 6 days on average. In all basins ASCAT GM observes melting at higher elevations relative to GDM, average of 5,328 m a.s.l. Systematic differences in glacier melting area determined by modeling and satellite observations indicate ASCAT may have suboptimal resolution, view geometry and/or polarimetry for delineating glacier melting at the process-scale in complex topography, especially in the ablation zone. This is the first step in examining the remote sensing products that could potentially be incorporated into hydrologic models to increase the accuracy of the hydrologic flow as well as the ability to estimate river discharge in other basins with limited data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The understudied winter: Evidence of how precipitation differences affect stream metabolism in a headwater

Climate change is causing pronounced shifts during winter in the US, including shortening the snow season, reducing snowpack, and altering the timing and volume of snowmelt-related runoff. These changes in winter precipitation patterns affect in-stream freeze-thaw cycles, including ice and snow cover, and can trigger direct and indirect effects on in-stream physical, chemical, and biological processes in ~60% of river basins in the Northern Hemisphere. We used high-resolution, multi-parameter data collected in a headwater stream and its local environment (climate and soil) to determine interannual variability in physical, chemical, and biological signals in a montane stream during the winter of an El Niño and a La Niña year. We observed ~77% greater snow accumulation during the El Niño year, which caused the formation of an ice dam that shifted the system from a primarily lotic to a lentic environment. Water chemistry and stream metabolism parameters varied widely between years. They featured anoxic conditions lasting over a month, with no observable gross primary production (GPP) occurring under the ice and snow cover in the El Niño year. In contrast, dissolved oxygen and GPP remained relatively high during the winter months of the La Niña year. These redox and metabolic changes driven by changes in winter precipitation have significant implications for water chemistry and biological functioning beyond the winter. Our study suggests that as snow accumulation and hydrologic conditions shift during the winter due to climate change, hot-spots and hot-moments for biogeochemical processing may be reduced, with implications for the downstream movement of nutrients and transported materials.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2024

Historical and projected future runoff over the Mekong River basin

The Mekong River (MR) crosses the borders and connects six countries, including China, Myanmar, Laos, Thailand, Cambodia, and Vietnam. It provides critical water resources and supports natural and agricultural ecosystems, socioeconomic development, and the livelihoods of the people living in this region. Understanding changes in the runoff of this important international river under projected climate change is critical for water resource management and climate change adaptation planning. However, research on long-term runoff dynamics for the MR and the underlying drivers of runoff variability remains scarce. Here, we analyse historical runoff variations from 1971 to 2020 based on runoff gauge data collected from eight hydrological stations along the MR. With these runoff data, we then evaluate the runoff simulation performance of five global hydrological models (GHMs) forced by four global climate models (GCMs) under the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP). Furthermore, based on the best simulation combination, we quantify the impact of future climate change on river runoff changes in the MR. The result shows that the annual runoff in the MR has not changed significantly in the past 5 decades, while the establishment of dams and reservoirs in the basin visibly affected the annual runoff distribution. The ensemble-averaged result of the Water Global Assessment and Prognosis version 2 (WaterGAP2; i.e. GHM) forced by four GCMs has the best runoff simulation performance. Under Representative Concentration Pathways (RCPs; i.e. RCP2.6, RCP6.0 and RCP8.5), the runoff of the MR is projected to increase significantly (p<0.05); e.g. 3.81 ± 3.47 m 3 s -1 a -1 (9 ± 8 % increase in 100 years) at the upper reach under RCP2.6 and 16.36 ± 12.44 m 3 s -1 a -1 (13 ±10 % increase in 100 years) at the lower reach under RCP6.0. In particular, under the RCP6.0 scenario, the increase in annual runoff is most pronounced in the middle and lower reaches, due to increased precipitation and snowmelt. Under the RCP8.5 scenario, the runoff distribution in different seasons varies obviously, increasing the risk of flooding in the wet season and drought in the dry season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Projecting end-of-century climate extremes and their impacts on the hydrology of a representative California watershed

Abstract. In California, it is essential to understand the evolution of water resources in response to a changing climate to sustain its economy and agriculture and to build resilient communities. Although extreme conditions have characterized the historical hydroclimate of California, climate change will likely intensify hydroclimatic extremes by the end of the century (EoC). However, few studies have investigated the impacts of EoC extremes on watershed hydrology. We use cutting-edge global climate and integrated hydrologic models to simulate EoC extremes and their effects on the water-energy balance. We assess the impacts of projected driest, median, and wettest water years under Representative Concentration Pathway (RCP) 8.5 on the hydrodynamics of the Cosumnes River basin. Substantial changes to annual average temperature (>+2.5 ∘C) and precipitation (>+38 %) will characterize the EoC extreme water years compared to their historical counterparts. A shift in the dominant form of precipitation, mostly in the form of rain, is projected to fall earlier. These changes reduce snowpack by more than 90 %, increase peak surface water and groundwater storages up to 75 % and 23 %, respectively, and drive the timing of peak storage to occur earlier in the year. Because EoC temperatures and soil moisture are high, both potential and actual evapotranspiration (ET) increase. The latter, along with the lack of snowmelt in the warm EoC, causes surface water and groundwater storages to significantly decrease in summer, with groundwater showing the highest rates of decrease. These changes result in more ephemeral EoC streams with more focused flow and increased storage in the mainstem of the river network during the summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Watershed zonation through hillslope clustering for tractably quantifying above- and below-ground watershed heterogeneity and functions

Abstract. In this study, we develop a watershed zonation approach for characterizing watershed organization and functions in a tractable manner by integrating multiple spatial data layers. We hypothesize that (1) a hillslope is an appropriate unit for capturing the watershed-scale heterogeneity of key bedrock-through-canopy properties and for quantifying the co-variability of these properties representing coupled ecohydrological and biogeochemical interactions, (2) remote sensing data layers and clustering methods can be used to identify watershed hillslope zones having the unique distributions of these properties relative to neighboring parcels, and (3) property suites associated with the identified zones can be used to understand zone-based functions, such as response to early snowmelt or drought and solute exports to the river. We demonstrate this concept using unsupervised clustering methods that synthesize airborne remote sensing data (lidar, hyperspectral, and electromagnetic surveys) along with satellite and streamflow data collected in the East River Watershed, Crested Butte, Colorado, USA. Results show that (1) we can define the scale of hillslopes at which the hillslope-averaged metrics can capture the majority of the overall variability in key properties (such as elevation, net potential annual radiation, and peak snow-water equivalent – SWE), (2) elevation and aspect are independent controls on plant and snow signatures, (3) near-surface bedrock electrical resistivity (top 20 m) and geological structures are significantly correlated with surface topography and plan species distribution, and (4) K-means, hierarchical clustering, and Gaussian mixture clustering methods generate similar zonation patterns across the watershed. Using independently collected data, we show that the identified zones provide information about zone-based watershed functions, including foresummer drought sensitivity and river nitrogen exports. The approach is expected to be applicable to other sites and generally useful for guiding the selection of hillslope-experiment locations and informing model parameterization.

58 GEOSCIENCES↗
DOE OSTI2024

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗
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↗