SEARCH · Search NASA

Results for “Surface hydrology”

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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

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

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

ENSO↗

DOE OSTI2023

Identifying hydrologic signatures associated with streamflow depletion caused by groundwater pumping

Abstract Groundwater pumping can reduce streamflow in nearby waterways (‘streamflow depletion’), a process which must be accounted for in integrated management of surface and groundwater resources. However, causal identification of streamflow depletion from hydrographs alone is challenging because pumping impacts are masked by other drivers of hydrologic variability. To identify potential indicators of streamflow depletion, we used synthetic hydrographs and an analytical streamflow depletion model to assess potential pumping impacts on specific hydrograph characteristics (‘hydrologic signatures’) for 215 streamgages spanning the conterminous United States (CONUS). We found that streamflow depletion commonly impacts signatures associated with seasonal and annual low flows and low flow recessions. The largest impacts occurred during dry years, suggesting streamflow depletion may be evident in dry years even where impacts are unmeasurable in wet years. Random forest models indicated that streamflow depletion could significantly impact Annual, Summer, and Fall signatures in most streams. Our finding that multiple hydrologic signatures are consistently responsive to streamflow depletion across CONUS suggests that the underlying hydrological processes linking pumping to streamflow reductions are consistent across diverse settings, information that will aid in identifying indicators of streamflow depletion from streamflow hydrographs.

Lapides, Dana A.↗

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 OSTI2022

Surface parameters and bedrock properties covary across a mountainous watershed: Insights from machine learning and geophysics

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

58 GEOSCIENCES↗

DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗

DOE OSTI2023

Hydrobiogeochemical Controls on the Delivery of Dissolved Organic Matter to Boreal Headwater Streams

Understanding controls on the delivery of dissolved organic matter (DOM) from terrestrial to aquatic systems is key for constraining carbon balances and fluxes in boreal headwater catchments. The largest export of DOM generally occurs during intense periods of water delivery to the landscape (e.g., rain events); however, the timing and magnitude of DOM is complicated by geomorphological, hydrometeorological, and biogeochemical variability. To better assess mechanisms controlling the delivery of DOM to boreal streams, DOM was investigated in two morphologically distinct catchments of an experimental forest watershed in western Newfoundland, one dominated by low relief wetlands and the other by steep hillslopes. The DOM was compared during two fall storms of similar magnitude but contrasting moisture conditions during the autumn transition. Variable concentrations and optical character highlighted antecedent conditions are important for the delivery of DOM to boreal streams. Concentration-discharge relationships with stream hydrograph separation analysis suggest preferential flowpaths through shallow mineral horizons as a key pathway for delivery of DOM in the hillslope dominated catchment compared to rapid input from near-stream wetlands in the low relief catchment. Loss of DOM via preferential pathways suggests an additional mechanism for surface soil carbon loss not directly informed by lateral flow and soil hydrology, which has important implications for our conceptualization and modeling of carbon fluxes in boreal systems. These results suggest watershed scale responses of DOM in boreal headwater catchments are complex, and better linkages between catchment scale hydrology and landscape hydrobiogeochemistry are required to constrain terrestrial to aquatic carbon fluxes in boreal landscapes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Investigating the Impact of Irrigation on Malaria Vector Larval Habitats and Transmission Using a Hydrology‐Based Model

Abstract A combination of accelerated population growth and severe droughts has created pressure on food security and driven the development of irrigation schemes across sub‐Saharan Africa. Irrigation has been associated with increased malaria risk, but risk prediction remains difficult due to the heterogeneity of irrigation and the environment. While investigating transmission dynamics is helpful, malaria models cannot be applied directly in irrigated regions as they typically rely only on rainfall as a source of water to quantify larval habitats. By coupling a hydrologic model with an agent‐based malaria model for a sugarcane plantation site in Arjo, Ethiopia, we demonstrated how incorporating hydrologic processes to estimate larval habitats can affect malaria transmission. Using the coupled model, we then examined the impact of an existing irrigation scheme on malaria transmission dynamics. The inclusion of hydrologic processes increased the variability of larval habitat area by around two‐fold and resulted in reduction in malaria transmission by 60%. In addition, irrigation increased all habitat types in the dry season by up to 7.4 times. It converted temporary and semi‐permanent habitats to permanent habitats during the rainy season, which grew by about 24%. Consequently, malaria transmission was sustained all‐year round and intensified during the main transmission season, with the peak shifted forward by around 1 month. Lastly, we evaluated the spatiotemporal distribution of adult vectors under the effect of irrigation by resolving habitat heterogeneity. These findings could help larval source management by identifying transmission hotspots and prioritizing resources for malaria elimination planning.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Hydraulic Conductivity Measurements, Utqiagvik (Barrow), Alaska, 2014

Six individual ice cores were collected from the Barrow Environmental Observatory in Barrow, Alaska, in May of 2013 as part of the Next Generation Ecosystem Experiment (NGEE). Each core was drilled at a different location to varying depths. After drilling, the cores were stored in coolers packed with dry ice and flown to Lawrence Berkeley National Laboratory (LBNL) in Berkeley, CA. 3-dimensional images of the cores were constructed using medical X-ray computed tomography (CT) scanner at 120kV. Hydraulic conductivity samples were extracted from these cores at LBNL Richmond Field Station in Richmond, CA, in February 2014 by cutting 5 to 8 inch segments using a chop saw. Samples were packed individually and stored at -20C freezing temperatures to minimize any changes in structure or loss of ice content prior to analysis. Hydraulic conductivity was determined through falling head tests using a permeameter [ELE International, Model #: K-770B] (Appendix A). Samples were placed in a latex membrane via a membrane stretcher while frozen. Use of a membrane stretcher made the membranes easier to secure and minimized contact with the sample. A clear polycarbonate sleeve, fabricated with a stainless steel ring at the bottom to keep the sleeve from floating, was placed around the sample inside the permeameter to minimize deformation during analysis. The permeameter was filled with water and 1.0 PSI of air was applied for confining pressure during sample defrost. Outflow valves were left open to allow for incremental thawing and samples were left to thaw for approximately 12 hours. After approximately 12 hours of thaw, initial falling head tests were performed. When the flow was significantly too fast or too slow, the analysis was stopped and the burette size was adjusted accordingly (i.e. a larger diameter burette was used for flows that were faster than desired or a smaller diameter burette was used for flows that were slower than desired). Two to four measurements were collected on each sample and collection stopped when the applied head load exceeded 25% change from the original load. Analyses were performed between 2 to 3 times for each sample. The final hydraulic conductivity calculations were computed using methodology of Das et al., 1985.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Airborne Measurements of Scale‐Dependent Latent Heat Flux Impacted by Water Vapor and Vertical Velocity Over Heterogeneous Land Surfaces During the CHEESEHEAD19 Campaign

Abstract The water vapor transport associated with latent heat flux (LE) in the planetary boundary layer (PBL) is critical for the atmospheric hydrological cycle, radiation balance, and cloud formation. The spatiotemporal variability of LE and water vapor mixing ratio ( r v ) are poorly understood due to the scale‐dependent and nonlinear atmospheric transport responses to land surface heterogeneity. Here, airborne in situ measurements with the wavelet technique are utilized to investigate scale‐dependent relationships among LE, vertical velocity ( w ) variance (), and r v variance () over a heterogeneous surface during the Chequamegon Heterogeneous Ecosystem Energy‐balance Study Enabled by a High‐density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign. Our findings reveal distinct scale distributions of LE, , and at 100 m height, with a majority scale range of 120 m–4 km in LE, 32 m–2 km in , and 200 m–8 km in . The scales are classified into three scale ranges, the turbulent scale (8–200 m), large‐eddy scale (200 m–2 km), and mesoscale (2–8 km) to evaluate scale‐resolved LE contributed by and . The large‐eddy scale in PBL contributes over 70% of the monthly mean total LE with equal parts (50%) of contributions from and . The monthly temporal variations mainly come from the first two major contributing classified scales in LE, , and . These results confirm the dominant role of the large‐eddy scale in the PBL in the vertical moisture transport from the surface to the PBL, while the mesoscale is shown to contribute an additional ∼20%. This analysis complements published scale‐dependent LE variations, which lack detailed scale‐dependent vertical velocity and moisture information.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Hydrologic response of artificially drained agricultural watersheds: insights from high-resolution integrated surface/subsurface simulations

Artificial drainage systems comprising subsurface networks of perforated pipes (tile drains) and engineered surface ditches are widely used to remove excess water from poorly drained agricultural regions. Artificial drainage lowers the water table by design but also has inadvertent effects on the watershed-scale hydrologic response with important implications for flood risk and nutrient exports. We investigated the effects of tile drains on watershed-scale hydrologic response in the Portage River, OH, Watershed using a high-resolution physically based integrated surface/subsurface hydrology model with recently developed capabilities to represent artificial drainage. Tile drains were found to enhance streamflow during times of low flow, generally consistent with previous studies. Streamflow flashiness was found to have a non-monotonic dependence on tile spacing with a minimum at intermediate spacings (∼50 m). Flashiness and the event hydrographs for small tile spacing were similar to the situation with no tiles, but flow paths from farm to stream were very different for those two end member cases, emphasizing the limitations of the stream hydrograph in characterizing hydrologic response. For typical tile spacings, peak flow can either be enhanced or attenuated by the presence of tiles, depending on the size of the event and the antecedent meteorological conditions. Tiles enhance peak flow when the event is below a threshold of ∼25 mm or when events arrive in dry conditions. Peak flow is reduced by tiles when events are large and arrive in conditions that are not overly dry. The dependence on event size and antecedent conditions is explained by differences in available storage and flow paths to the streams. These results provide additional insights into how tile drainage modulates event-scale hydrologic response, an important control on flood generation mechanisms and nutrient exports.

agricultural watershed↗

DOE OSTI2022

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Abstract. Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (tens of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) – Daily Surface Weather and Climatological Summaries (Daymet), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30 % of the watershed area). However, the 4 km-resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using the temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), and sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications for model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗

DOE OSTI2023

Integrating fine root diameter and watershed mapping to characterize rhizosphere hydrology

Root morphology and soil hydraulic characteristics were integrated using watershed distance mapping to show water distribution and uptake across the plant-soil interface. Poplar (Populus deltoides, P. trichocarpa), maize (Zea mays), juniper (Juniperus virginiana), grape (Vitis rotundifolia) and maple (Acer saccharum) seedlings were grown in sand, after which root diameter and soil water dynamics were assessed via sequential neutron radiography. Three local soil regions (root-soil interface or edge, rhizosphere, bulk soil) were classified based on both radial distance from the root surface and diameter of the nearest root, from which changes in water content and distribution were characterized using digital image processing. Water content dynamics across the rhizosphere showed two different species-independent processes: a consistently elevated water content at the root-soil edge interface which increased with root diameter, and hysteresis as the rhizosphere transitioned to bulk soil (~0.5 cm from the root), independent of root diameter. Water uptake per unit root surface area declined exponentially with root diameter, independent of species. Results highlight the species-independent hydrologic characteristics of the rhizosphere and the potential for evaluating them in a local spatially connected soil context. Avenues for improved integration of soil and root characteristics are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Geophysical and Environmental Monitoring Data, and Subsurface Flow Modelling Results for Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO

This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗