Search NASA⌕ Search

SEARCH · Search NASA

Results for “Watershed transport”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Model output from simulations of manganese-carbon interactions in temperate forest soil profiles

This archive contains model output, code, and scripts for simulations of coupled manganese-carbon cycling in temperate forest soil profiles. These model results were generated as part of a study investigating how manganese availability influences soil organic carbon stocks and demonstrating a new model framework for coupling carbon and manganese cycling. The simulations were in support of a manuscript: "Modeling interactive effects of manganese bioavailability, nitrogen deposition, and warming on soil carbon storage." The study addresses the research questions: How does Mn bioavailability, as driven by subsurface mineral properties, pH, and redox status, affect temperate forest soil organic carbon and litter carbon stocks?How is the relationship between Mn bioavailability and carbon cycling affected by changes in temperature and nitrogen deposition?"Model simulations were conducted in a reactive transport modeling framework using PFLOTRAN coupled to python. Multiple model simulations testing different Mn-bearing mineral solubilities, hydrological patterns, nitrogen deposition rates, and temperatures are included. Soil properties, including total and exchangeable Mn concentrations, are based on values reported for the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), a temperate forested watershed in central Pennsylvania, U.S.A where Mn cycling through vegetation has been documented.File formats include netCDF (.nc), python script (.py), shell script (.sh), plain text PFLOTRAN input file (.in), and plain text PFLOTRAN database file (.dat), and gzipped tar archive (tar.gz).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Seasonal mixing from intermittent flow drives concentration-discharge behaviour in a stream affected by coal mine drainage

Abandoned mining operations continue to severely degrade many ecosystems worldwide by releasing acidic water and/or heavy metals into surface and groundwater. Contaminant concentrations in affected streams vary with discharge in patterns that reflect both geochemical reactions and variable mixing of contaminated and non-contaminated waters. However, controls on concentration-discharge (C-Q) patterns remain unclear, particularly for constituents that experience changing solubility across redox and pH gradients. Understanding the C-Q behaviour of contaminants aids in predicting both downstream transport and effects on aquatic life under variable flow. Here, we examined the C-Q behaviours of non-reactive (Na, K, Ca, Mg, Cl - ) and reactive (Fe, Mn, Al, H + , SO 4 2- ) solutes in a stream contaminated with acid mine drainage in northeastern Ohio, USA. Concentration-discharge patterns at the watershed outlet primarily reflected mixing of contaminated baseflow with intermittent inputs of high pH water draining from a passive limestone treatment system into the stream. The treatment system acted as an ephemeral tributary that mitigated contamination in the stream by diluting solutes, raising pH, and driving metal precipitation, but only when flow was present during wet seasons. Consequently, AMD-derived reactive solutes (H + , Fe, Mn, Al) decreased with increasing stream discharge while relatively conservative solutes (e.g., Ca, Mg, K, Na) decreased only slightly or were chemostatic. This study highlights both the unique C-Q patterns of reactive solutes when compared to those of non-reactive solutes and the potential for intermittent streams to control C-Q behaviour in headwater catchments.

58 GEOSCIENCES↗
DOE OSTI2023

Model Data for the Mesh Convergence Study Demonstrating Benefits of Mixed-polyhedral Mesh in Integrated Hydrology Simulations

This archived model data is related to a study introducing a unique method that employs a stream-aligned mixed-polyhedral mesh to effectively and accurately represent river valleys, stream corridors, and narrow engineered channels in integrated hydrology simulations. The study finds that utilizing stream-aligned mixed-polyhedral meshes in integrated hydrology simulations achieves accuracy on par with a finely refined TIN-based mesh while markedly diminishing computational costs. This archive contains scripts and data files needed to generate the ATS model input, including mesh and ATS input files, for all mesh scenarios using the Watershed Workflow package. Additionally, this archive also provides key outputs from the model simulations that are used in the analysis and post-processing scripts to reproduce figures in the manuscript. The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts. The input files for the ATS model, open-source integrated hydrology, and transport model are in XML format and can be edited in any commonly used text editors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulation of high-frequency dissolved oxygen dynamics in a shallow estuary, the Corsica River, Chesapeake Bay

Understanding shallow water biogeochemical dynamics is a challenge in coastal regions, due to the presence of highly variable land-water interface fluxes, tight coupling with sediment processes, tidal dynamics, and diurnal variability in biogeochemical processes. While the deployment of continuous monitoring devices has improved our understanding of high-frequency (12 - 24 hours) variability and spatial heterogeneity in shallow regions, mechanistic modeling of these dynamics has lagged behind conceptual and empirical models. The inherent complexity of shallow water systems is represented in the Corsica River estuary, a small basin within the Chesapeake Bay ecosystem, where abundant monitoring data have been collected from long-term monitoring stations, continuous monitoring sensors, synoptic sensor surveys, and measurements of sediment-water fluxes. A state-of-the-art modeling system, the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), was applied to the Corsica domain with a high-resolution grid and nutrient loads from the most recent version of the Chesapeake Bay watershed model. The Corsica SCHISM model reproduced observed high-frequency variability in dissolved oxygen, as well as seasonal variability in chlorophyll-a and sediment-water fluxes. Time-series signal analyses using Empirical Model Decomposition and spectral analysis revealed that the diurnal and M2 tide frequencies are the dominant high-frequency modes and physical transport contributes a larger share to dissolved oxygen budgets than biogeochemical processes on an hourly time scale. Heterogeneity and patchiness in dissolved oxygen resulting from phytoplankton distributions and geometry-driven eddies amplify the physical transport effect, and on longer time scales oxygen is controlled more by photosynthesis and respiration. Our simulation demonstrates that interactions among physical and biological dynamics generate complex high-frequency variability in water quality and non-linear reposes to nutrient loading and environmental forcing in shallow water systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Gap-filled meteorological data (2011-2020) and modeled potential evapotranspiration data from the KCOMTCRE2 WeatherUnderground weather station, from the East River Watershed, Colorado

This dataset is a gap-filled meteorological dataset (years 2011-2020) that includes modeled potential evapotranspiration from a station near the pumphouse on the East River (KCOMTCRE2 WeatherUnderground ) near the top of the ridge at the Pumphouse PLM wells. These datasets are useful as model inputs. The KOMTCRE2 station is located on Mt. Crested Butte in Colorado. More information can be found at: https://www.wunderground.com/dashboard/pws/KCOMTCRE2/table/2020-12-7/2020-12-7/daily. The purpose of this dataset is to use as input to hydrological and flow-reactive transport models as meteorological forcing. One attached file includes metadata files with definitions for all variables and units. Subsequent files include the same time series dataset, but aggregated and gap-filled at different time-stamp intervals. These include 5 minute time-series, 15 minute time-series, 30 minute time-series, 1 hour time-series, daily time-series, and hourly and daily modeled evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

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

Community-Engaged Modeling of Urban Flood Adaptation Pathways

Climate change is intensifying the hydrologic cycle, leading to more frequent and severe rainfall-driven (pluvial) flooding in urban areas. In the mid-Atlantic US cities, aging and under-designed stormwater infrastructure is increasingly strained by these events, resulting in recurring damage to property and disruptions to transportation networks. In this study, we combine community engagement with hydrologic modeling to develop and evaluate potential urban flood adaptation strategies. Over a three-year period, local technical experts and community representatives met regularly to discuss flooding concerns, identify priorities, and co-develop adaptation strategies. These discussions informed the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed, the focus location of this study. The flooding model integrates complex surface and subsurface stormwater infrastructure data, local expert knowledge, and community insights. We simulate stakeholder-prioritized adaptations, such as green and gray infrastructure strategies. Model results demonstrate that enhanced infrastructure maintenance is the most effective adaptation for reducing flood depths, but has varied effects across the watershed, and can increase flooding in some locations. Spatially concentrated greening provides limited benefit to the watershed as a whole, but moderate benefit in community priority areas. Together, these adaptations have the potential to reduce flood depths by as much as 58% in some locations, greatly reducing property damage and mobility impacts, primary concerns of stakeholders. Future work will implement robust optimization tools to search for adaptations which meet stakeholder objectives and perform highly under varied future climate conditions. This work contributes to the expanding literature on collaborative modeling and demonstrates that community-engaged approaches can enhance model credibility and generate more actionable insights for communities seeking to strengthen climate resilience.

Spangler, Ava [Pennsylvania State University] (ORC↗
DOE OSTI2026

Community-Informed Urban Flood Modeling for Impact Mitigation

Climate change is intensifying the hydrologic cycle, leading to more frequent and severe rainfall-driven (pluvial) flooding in urban areas. In the mid-Atlantic US cities, aging and under-designed stormwater infrastructure is increasingly strained by these events, resulting in recurring damage to property and disruptions to transportation networks. In this study, we combine community engagement with hydrologic modeling to develop and evaluate potential urban flood adaptation strategies. Over a three-year period, local technical experts and community representatives met regularly to discuss flooding concerns, identify priorities, and co-develop adaptation strategies. These discussions informed the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed, the focus location of this study. The flooding model integrates complex surface and subsurface stormwater infrastructure data, local expert knowledge, and community insights. We simulate stakeholder-prioritized adaptations, such as green and gray infrastructure strategies. Model results demonstrate that enhanced infrastructure maintenance is the most effective adaptation for reducing flood depths, but has varied effects across the watershed, and can increase flooding in some locations. Spatially concentrated greening provides limited benefit to the watershed as a whole, but moderate benefit in community priority areas. Together, these adaptations have the potential to reduce flood depths by as much as 58% in some locations, greatly reducing property damage and mobility impacts, primary concerns of stakeholders. Future work will implement robust optimization tools to search for adaptations which meet stakeholder objectives and perform highly under varied future climate conditions. This work contributes to the expanding literature on collaborative modeling and demonstrates that community-engaged approaches can enhance model credibility and generate more actionable insights for communities seeking to strengthen climate resilience.

Baltimore MD↗
DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old‐growth to second‐growth forests in headwater ecosystems

Abstract Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear‐cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old‐growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear‐cut harvesting, from mixed‐species old‐growth to lodgepole pine‐dominated second‐growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3‐year period (2016–2018) at the Fraser Experimental Forest in north‐central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old‐ versus second‐growth forest. Fluorescence spectroscopy showed that the DOM exported from the old‐growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second‐growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old‐growth forest relative to second growth. Old‐growth and second‐growth forests are co‐mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

Fegel, Timothy S.↗
DOE OSTI2021

Modeling the Impact of Riparian Hollows on River Corridor Nitrogen Exports

Recent studies in snowmelt-dominated catchments have documented changes in nitrogen (N) retention over time, such as declines in watershed exports of N, though there is a limited understanding of the controlling processes driving these trends. Working in the mountainous headwater East River Colorado watershed, we explored the effects of riparian hollows as N-cycling hotspots and as important small-scale controls on observed watershed trends. Using a modeling-based approach informed by remote sensing and in situ observations, we simulated the N-retention capacity of riparian hollows with seasonal and yearly hydrobiogeochemical perturbations imposed as drivers. We then implemented a scaling approach to quantify the relative contribution of riparian hollows to the total river corridor N budget. We found that riparian hollows primarily serve as N sinks, with N-transformation rates significantly limited by periods of enhanced groundwater upwelling and promoted at the onset of rainfall events. Given these observed hydrologic controls, we expect that the nitrate ( NO 3 - ) sink capacity of riparian hollows will increase in magnitude with future climatic perturbations, specifically the shift to more frequent rainfall events and fewer snowmelt events, as projected for many mountainous headwater catchments. Our current estimates suggest that while riparian hollows provision ~5–20% of NO 3 - to the river network, they functionally act as inhibitors to upland NO 3 - reaching the stream. Our work linking transient hydrological conditions to numerical biogeochemical simulations is an important step in assessing N-retaining features relative to the watershed N budget and better understanding the role of small-scale features within watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Riverine Plastic Pollution: Field Sampling Protocol and Implementation in U.S. Rivers

Riverine plastic pollution is increasingly being recognized as a serious environmental concern. The flux of plastic from U.S. rivers into the oceans is estimated to be on the order of hundreds of tons per year. Plastic pollution has been detected in all major U.S. rivers but the observations did not measure metrics that could be used for mass flux calculations. The Department of Energy is addressing the problem of riverine plastic pollution with the Waterborne Plastics Resource Assessment and Debris Characterization (WaterPACT) project. Phase one of the WaterPACT project includes the development of a sampling protocol, field sampling events that implement the protocol, and numerical modeling studies. The field sampling will quantify and characterize the plastic pollution and the numerical modeling studies will analyze the movement of the plastic pollution by rivers. This report details the WaterPACT sampling protocol and its implementation in four major U.S. rivers. The WaterPACT sampling protocol was implemented in the Mississippi, Columbia, Delaware, and Los Angeles rivers with at least three sampling events in each river. The rivers represent four types ranging from the large Mississippi with its vast agricultural watershed to the small Los Angeles with its highly engineered urban waterway. The three sampling events covered a variety of discharge conditions from low drought flow to extreme storm events. The data obtained by the WaterPACT sampling of the four rivers will be combined to provide an invaluable dataset to researchers studying the source and chemical composition of plastic pollution and modelers estimating the flux of plastic being transported to and released into the oceans.

13 HYDRO ENERGY↗
DOE OSTI2021

Bed form-induced hyporheic exchange and geochemical hotspots

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old-growth to second-growth forests in headwater ecosystems

Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear-cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old-growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear-cut harvesting, from mixed-species old-growth to lodgepole pine-dominated second-growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3-year period (2016–2018) at the Fraser Experimental Forest in north-central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old- versus second-growth forest. Fluorescence spectroscopy showed that the DOM exported from the old-growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second-growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old-growth forest relative to second growth. Old-growth and second-growth forests are co-mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗