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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Anticipating how rain-on-snow events will change through the 21st century: lessons from the 1997 new year’s flood event

The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5° C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥ 1.7° C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Median Isn’t the Message: soil nutrient hot spots have a disproportionate influence on biogeochemical structure across years, seasons, and depths

Abstract Soil nutrient distribution is heterogeneous in space and time, potentially altering nutrient acquisition by trees and microorganisms. Ecologists have distinguished “hot spots” (HSs) as areas with enhanced and sustained rates of nutrient fluxes relative to the surrounding soil matrix. We evaluated the spatial and temporal patterns in nutrient flux HSs in two mixed-conifer forest soils by repeatedly sampling the soil solution at the same spatial locations (horizontally and vertically) over multiple seasons and years using ion exchange resins incubated in situ. The climate of these forests is Mediterranean, with intense fall rains occurring following summers with little precipitation, and highly variable winter snowfall. Hot spots formed most often for NO 3 − and Na + . Although nutrient HSs often occurred in the same spatial location multiple times, HSs persisted more often for PO 4 3− NH 4 + , and NO 3 − , and were more transient for Ca 2+ , Mg 2+ , and Na + . Sampling year (annual precipitation ranged from 558 to 1223 mm) impacted the occurrence of HSs for most nutrients, but season was only significant for PO 4 3− , NH 4 + , NO 3 − , and Na + , with HSs forming more often after fall rains than after spring snowmelt. The frequency of HSs significantly decreased with soil depth for all nutrients, forming most commonly immediately below the surficial organic horizon. Although HSs accounted for less than 17% of the sampling volume, they were responsible for 56–88% of PO 4 3− , NH 4 + , and NO 3 − resin fluxes. Our results suggest that macronutrient HSs have a disproportional contribution to soil biogeochemical structure, with implications for vegetation nutrient acquisition strategies and biogeochemical models. Graphical abstract

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reducing nitrous oxide emissions and optimizing nitrogen-use efficiency in dryland crop rotations with different nitrogen rates

Recent interests in improving agricultural production while minimizing environmental footprints emphasized the need for research on management strategies that reduce nitrous oxide (N 2 O) emissions and increase nitrogen-use efficiency (NUE) of cropping systems. This study aimed to evaluate N 2 O emissions, annualized crop grain yield, emission factor, and yield-scaled- and NUE-scaled N 2 O emissions under continuous spring wheat (Triticum aestivum L.) (CW) and spring wheat-pea (Pisum sativum L.) (WP) rotations with four N fertilization rates (0, 50, 100, and 150 kg N ha(-1)). The N 2 O fluxes peaked immediately after N fertilization, intense precipitation, and snowmelt, which accounted for 75-85% of the total annual flux. Cumulative N 2 O flux usually increased with increased N fertilization rate in all crop rotations and years. Annualized crop yield and NUE were greater in WP than CW for 0 kg N ha(-1) in all years, but the trend reversed with 100 kg N ha(-1) in 2013 and 2015. Crop yield maximized at 100 kg N ha(-1), but NUE declined linearly with increased N fertilization rate in all crop rotations and years. As N fertilization rate increased, N fertilizer-scaled N2O flux decreased, but NUE-scaled N 2 O flux increased non-linearly in all years, regardless of crop rotations. The yield-scaled N 2 O flux decreased from 0 to 50 kg N ha -1 and then increased with increased N fertilization rate. Because of non-significant difference of N 2 O fluxes between 50 and 100 kg N ha(-1), but increased crop yield, N 2 O emissions can be minimized while dryland crop yields and NUE can be optimized with 100 kg N ha -1 , regardless of crop rotations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Regional patterns in hydrologic response, a new three-component metric for hydrograph analysis and implications for ecohydrology, Northwest Volcanic Aquifer Study Area, USA

Spatial patterns of hydrologic response were examined for the Northwest Volcanic Aquifer Study Area (NVASA). The utility of established hydrograph-separation methods for assessing hydrologic response in permeable volcanic terranes was assessed and a new three-component metric for hydrograph analysis was developed. The new metric, which partitions streamflow into subcomponents defined by the timescales of hydrologic response (e.g., fast-runoff, intermediate-interflow and slow-baseflow), was used to gain a fundamental understanding of the regional hydrology, investigate sub-regional differences, influencing factors, and ecohydrological implications. The combined effects of NVASA’s physiography, climate and geology create a strongly coupled surface-groundwater system that produces copious baseflow and limited quantities of runoff and interflow. Patterns of hydrologic response are influenced by the type and rate of precipitation and permeability of the underlying geology. Under variable precipitation conditions the hydrologic response of volcanic terranes with similar permeability and subsurface-storage capacity can be significantly different. From a water management and ecohydrology perspective, understanding regional patterns of hydrologic response and sub-regional differences is fundamental. Results indicate that minimum-flow methods provide the most conservative estimate of baseflow and may be the most robust for filtering out snowmelt bias in baseflow estimates. Baseflow contributes ~75% of the perennial streamflow across the NVASA and represents a critical component of the regional water supply that provides critical cold-water habitat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

AGFormer: Adaptive Spatiotemporal graph informed transformer for multi-reservoir inflow forecasting

Accurate reservoir inflow forecasting is crucial for effective water resource management, yet most machine learning models focus on single-reservoir prediction and overlook spatial dependencies among hydrologically connected reservoirs. Here, we propose AGFormer (Adaptive Graph-Informed Transformer), an end-to-end framework that integrates adaptive graph learning with temporal sequence modeling for multi-reservoir inflow forecasting. A shared encoder and graph attention mechanism generate reservoir-specific embeddings, which are then processed by the Transformer-based encoder–decoder for multi-step inflow forecasting. We also introduce a pretraining paradigm to learn robust temporal embeddings from misaligned historical records. Evaluated on 30 reservoirs in the Upper Colorado River Basin, AGFormer achieves superior seven-day-ahead forecasts, with NSE > 0.75 for 20 reservoirs—outperforming Encoder–Decoder LSTM, GCN+LSTM, and Transformer baselines. Adaptive graph learning captures dynamic inter-reservoir dependencies, and feature attribution aligns with snowmelt-driven hydrology. Incorporating forecasted meteorological inputs further enhances accuracy, demonstrating AGFormer’s potential to support reservoir management under dynamic hydrological conditions.

Adaptive graph learning↗
DOE OSTI2022

Simulation of anoxic lenses as exporters of reactivity in alluvial aquifer sediments

Sedimentary interfaces between contrasting hydrogeological facies in alluvial aquifers drive the development of biogeochemical interfaces that influence subsurface and surface water quality. Here, we calibrate a reactive transport model on a series of dual-domain column experiments, where centimeter-scale, low-permeability, organic-rich anoxic lenses are embedded in coarser-grained aquifer material. Simulations explicitly account for C, Fe, and S cycling at the interface between the lenses and the aquifer in water-saturated conditions. Our results highlight the role of fine-grained, organic-rich inclusions not only as sources or sinks for redox-sensitive species, but also as exporters of nutrients that stimulate downgradient biogeochemical cycling. By releasing large amounts of organic carbon into the surrounding aquifer, such lenses drive the development of proximal secondary reduction zones ("halos"), characterized by high microbial activity (e.g., sulfate reduction) and accumulation of reduced reaction products (e.g., iron sulfide). If similar secondary reduction zones develop within the hydraulically conductive domain of an aquifer, they would be highly susceptible to changes in hydrologic conditions, for instance oxygen pulses associated with seasonal snowmelt. Furthermore, our results also emphasize the limitations of relying solely on aqueous species measurements to inform reactivity in systems where fast redox cycling and/or sizeable particulate transport may limit the signature of reactivity in the dissolved phase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Vertical Hydrologic Exchange Flows Control Methane Emissions from Riverbed Sediments

CH 4 emissions from inland waters are highly uncertain in the current global CH 4 budget, especially for streams, rivers, and other lotic systems. Previous studies have attributed the strong spatiotemporal heterogeneity of riverine CH 4 to environmental factors such as sediment type, water level, temperature, or particulate organic carbon abundance through correlation analysis. However, a mechanistic understanding of the basis for such heterogeneity is lacking. Here, in this study, we combine sediment CH 4 data from the Hanford reach of the Columbia River with a biogeochemical-transport model to show that vertical hydrologic exchange flows (VHEFs), driven by the difference between river stage and groundwater level, determine CH 4 flux at the sediment–water interface. CH 4 fluxes show a nonlinear relationship with the magnitude of VHEFs, where high VHEFs introduce O 2 into riverbed sediments, which inhibit CH 4 production and induce CH 4 oxidation, and low VHEFs cause transient reduction in CH 4 flux (relative to production) due to reduced advective CH 4 transport. In addition, VHEFs lead to the hysteresis of temperature rise and CH 4 emissions because high river discharge caused by snowmelt in spring leads to strong downwelling flow that offsets increasing CH 4 production with temperature rise. Our findings reveal how the interplay between in-stream hydrologic flux besides fluvial-wetland connectivity and microbial metabolic pathways that compete with methanogenic pathways can produce complex patterns in CH 4 production and emission in riverbed alluvial sediments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ion Clusters Reveal the Sources, Impacts, and Drivers of Freshwater Salinization

Population growth, land use change, climate change, and natural resource extraction are driving the salinization of freshwater resources worldwide. Reversing these trends will require data-centric approaches that identify salt sources, environmental drivers, and ecosystem responses. In this study, we applied principal component analysis and hierarchical clustering to identify ion covariance patterns, or “ion clusters,” in Broad Run, an urban stream in the Mid-Atlantic United States. These clusters correspond to distinct hydrologic regimes and reveal specific salinization risks: (1) phosphorus pollution mobilized during summer storms (Cluster 1); (2) elevated concentrations of sulfate and bicarbonate during baseflow (Cluster 2), likely reflecting groundwater discharge; and (3) elevated specific conductance and sodium, chloride, and potassium ion concentrations during snowmelt and rain-on-snow events (Cluster 3), driven by deicer and anti-icer wash-off. These ion fingerprints offer a transferable framework for diagnosing salt sources, assessing ecological risk, and identifying management targets. Our findings underscore the need for next-generation stormwater infrastructure and smart growth policies to protect aquatic life in rapidly urbanizing watersheds.

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

Persistence and Plasticity in Conifer Water–Use Strategies

The selective use of seasonal precipitation by vegetation is critical to understanding the residence time and flow path of water in watersheds, yet there are limited datasets to test how climate alters these dynamics. Here, we use measurements of the seasonal cycle of tree ring δ 18 O for two widespread conifer species in the Rocky Mountains of North America to provide a multi-decadal depiction of the seasonal origins of forest water use. The results show that while the conifer tree stands had a dominant preference for use of snowmelt, there were multi-annual periods over the last four decades when use of summer precipitation was preferential. Utilization of summer rain emerged during years with increased snowfall and tree growth, suggesting that summer rain enhanced the transpiration stream only during the periods of highest water use. We hypothesize this could be explained through shallowing of the root profile during wetter periods and/or through the influence of changing water table depths on the residence time of summer precipitation in the root zone. We suggest the tree ring proxy approach used here could be applied in other watersheds to provide critical insight into the temporal dynamics of plant water use that could not be inferred from short measurement campaigns. Furthermore, these data on the seasonal origins of forest water are critical for understanding forest vulnerability to drought, the processes that affect precipitation pathways and residence time in watersheds and the interpretation of tree ring proxy data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Technology Pathways Could Help Drive the U.S. West Coast Grid's Exposure to Hydrometeorological Uncertainty

Abstract Previous studies investigating deep decarbonization of bulk electric power systems and wholesale electricity markets have not sufficiently explored how future grid pathways could affect the grid's vulnerability to hydrometeorological uncertainty on multiple timescales. Here, we employ a grid operations model and a large synthetic weather ensemble to “stress test” a range of future grid pathways for the U.S. West Coast developed by ReEDS, a well‐known capacity planning model. Our results show that gradual changes in the underlying capacity mix from 2020 to 2050 can cause significant “re‐ranking” of weather years in terms of annual wholesale electricity prices (with “good” years becoming bad, and vice versa). Nonetheless, we find the highest and lowest ranking price years in terms of average electricity price remain mostly tied to extremes in hydropower availability (streamflow) and load (summer temperatures), with the strongest sensitivities related to drought. Seasonal dynamics seen today involving spring snowmelt and hot, dry summers remain well‐defined out to 2050. In California, future supply shortfalls in our model are concentrated in the evening and occur mostly during periods of high temperature anomalies in late summer months and in late winter; in the Pacific Northwest, supply shortfalls are much more strongly tied to negative streamflow anomalies. Under our more robust sampling of stationary hydrometeorological uncertainty, we also find that the ratio of dis‐patchable thermal (i.e., natural gas) capacity to wind and solar required to ensure grid reliability can differ significantly from values reported by ReEDS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Exploring Sources of Surface Bias in HRRR Using New York State Mesonet

In recent years, there has been increasing demand for applications of short-term forecasting of renewable energy potential and assessments of the likelihood of extreme weather events using the High-Resolution Rapid Refresh (HRRR) model. Examining the biases in the newest version of HRRR is necessary to promote further model development. Using data from one of the most comprehensive and dense monitoring networks, New York State Mesonet (NYSM), we evaluate the HRRR version 3 meteorological fields for an entire year. In this work, the land-atmosphere-cloud coupling system is evaluated as an integrated whole. We investigate the physical processes influencing the soil hydrological balance and the thermodynamic interactions, from surface fluxes up to the level of boundary layer convection from both temporal (seasonal and diurnal) and spatial perspectives. Results show that the model 2 m temperature and humidity biases are seasonally dependent, with warm and dry bias present during the warm season, and an extreme nocturnal cold bias in winter. The summer warm bias includes both a land-surface-induced bias and a cloud-induced bias. Inaccurate representation of energy partition and soil hydrological process across different land use types as well as a hydrological bias in describing spring snowmelt are identified as the main source of the land-surface-induced bias. A feedback loop linking cloud presence, flux changes, and temperature contributes to the cloud-induced bias. The positive solar radiation bias increases from clear sky to overcast sky conditions. The most significant bias occurs during overcast and thick cloud conditions associated with frontal passage and thunderstorms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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

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

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An Improved Zhang's Dynamic Water Balance Model Using Budyko-Based Snow Representation for Better Streamflow Predictions

Understanding the water balance of a catchment in relation to its regional climate forcings and catchment characteristics is critical for predicting current and future water resources amid changing climate and land cover. This study intends to improve Zhang's monthly water balance model (a physics-based conceptual hydrologic model) that reflects the physical partitioning process of the hydrological cycle at the basin level based on regional climate and catchment characteristics. The existing model does not include snow process and has confronted evident limitations in snow-affected areas, which is a critical aspect since snowmelt water has been a significant source of water resources for many regions, especially in the temperate and frigid zones. We introduce a snow module based on surface energy balance and Budyko-limits on melting and combine it with the existing water balance equations. Moreover, monthly parameterization is applied to the model to better explain the time-varying hydrological characteristics of a catchment. The proposed model involves five different monthly parameters, which determine the physical partitioning process of the hydrological cycle, and they are regionally calibrated and validated under Budyko-type constraints. The model is applied to 1,210 basins across the continental United States (CONUS), and the simulated streamflow is compared to the observed data. The proposed model significantly outperformed the original model, improving the median NSE by 31% (from 0.51 to 0.67) and increasing the number of catchments with an acceptable NSE by 58%. The spatial variability of the basin characteristics across the CONUS is also investigated based on the calibrated parameters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Snow Albedo Feedbacks Enhance Snow Impurity-Induced Radiative Forcing in the Sierra Nevada

This study employs a fully coupled meteorology-chemistry-snow model to investigate the impacts of light-absorbing particles (LAPs) on snow darkening in the Sierra Nevada. After comprehensive evaluation with spatially and temporally complete satellite retrievals, the model shows that LAPs in snow reduce snow albedo by 0.013 (0-0.045) in the Sierra Nevada during the ablation season (April-July), producing a midday mean radiative forcing (RF) of 4.5 W m -2 which increases to 15-22 W m -2 in July. LAPs in snow accelerate snow aging processes and reduce snow cover fraction, which doubles the albedo change and radiative forcing caused by LAPs. The impurity-induced snow darkening effects decrease snow water equivalent and snow depth by 20 mm and 70 mm in June in Sierra Nevada bighorn sheep habitat. The earlier snowmelt reduces root zone soil moisture by 20%, deteriorating the forage productivity and playing a negative role in the survival of bighorn sheep.

54 ENVIRONMENTAL SCIENCES↗