Search NASASearch

SEARCH · Search NASA

Results for “Subsurface 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.

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

At least 127 records · Page 7

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This 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. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES

Acid Sulfate Alteration on Mars

A variety of mineralogical and geochemical indicators for aqueous alteration on Mars have been identified by a combination of surface and orbital robotic missions, telescopic observations, characterization of Martian meteorites, and laboratory and terrestrial analog studies. Acid sulfate alteration has been identified at all three landing sites visited by NASA rover missions (Spirit, Opportunity, and Curiosity). Spirit landed in Gusev crater in 2004 and discovered Fe-sulfates and materials that have been extensively leached by acid sulfate solutions. Opportunity landing on the plains of Meridiani Planum also in 2004 where the rover encountered large abundances of jarosite and hematite in sedimentary rocks. Curiosity landed in Gale crater in 2012 and has characterized fluvial, deltaic, and lacustrine sediments. Jarosite and hematite were discovered in some of the lacustrine sediments. The high elemental abundance of sulfur in surface materials is obvious evidence that sulfate has played a major role in aqueous processes at all landing sites on Mars. The sulfate-rich outcrop at Meridiani Planum has an SO3 content of up to 25 wt.%. The interiors of rocks and outcrops on the Columbia Hills within Gusev crater have up to 8 wt.% SO3. Soils at both sites generally have between 5 to 14 wt.% SO3, and several soils in Gusev crater contain around 30 wt.% SO3. After normalization of major element compositions to a SO3-free basis, the bulk compositions of these materials are basaltic, with a few exceptions in Gusev crater and in lacustrine mudstones in Gale crater. These observations suggest that materials encountered by the rovers were derived from basaltic precursors by acid sulfate alteration under nearly isochemical conditions (i.e., minimal leaching). There are several cases, however, where acid sulfate alteration minerals (jarosite and hematite) formed in open hydrologic systems, e.g., in Gale crater lacustrine mudstones. Several hypotheses have been suggested for the aqueous formation of sulfate-bearing phases under acidic conditions on the surface of Mars including (1) sulfuric acid weathering of basaltic materials; (2) oxidative weathering of ultramafic igneous rocks containing sulfides; (3) acid fog weathering of basaltic materials, and (4) near-neutral pH subsurface solutions rich in Fe2(+) that were rapidly oxidized to Fe3(+), which produced excess acidity as iron was oxidized on exposure to O2 or photo-oxidized by ultraviolet radiation at the martian surface. Next, we briefly describe evidence for these hypothesis.

Ming, D. W.

The Stop-and-Go Mechanism: Towards an Integrated Approach to Model Seismicity, Outgassing, Deformation, and Thermal Unrest at Active Volcanoes

Connecting the geophysical and geochemical signals recorded at and above the surface of volcanoes with source mechanisms is fundamental to understand transitions from quiescence to eruption, and to integrate the behavior of volcanoes with their regional seismotectonic and hydrological context. Some of these signals include shallow volcanic tremor, a long-lasting (minutes-to-years) ground vibration detected in volcanic areas during unrest; volcanic outgassing, which shows intriguing periodicities over multiple scales, as revealed by space-borne instruments and by ground-based, high-frequency (~1 Hz), sampling methods; inflation and deflation of volcanic edifices, commonly associated with subsurface volume changes; and large-scale (from a few to tens of km2), low-temperature, thermal anomalies, a newly-discovered signal that has been observed to emerge on volcanic flanks from months-to-years prior to gas- and magma-driven eruptions. Traditionally, these geophysical and geochemical signals have been modelled independently, which limits our interpretation of subsurface processes and thus our assessment of unrest. In this work, we show through lumped-parameter models, numerical simulations, and preliminary laboratory experiments that many of the signals recorded around volcanoes can be integrated into a common “stop-and-go” mechanism. In particular, we found that they can arise spontaneously when taking into account the feedbacks between the temporary accumulation of gas and/or magma in the crust (“stop”), and the diffusive transfer of gas and heat towards the surface (“go”). For example, we find that the “stop-and-go” mechanism can explain: (i) the emergence of monochromatic, broadband, and harmonic tremor; (ii) the periodic components commonly observed in volcanic outgassing time series; (iii) the link between large-scale thermal anomalies and shallow hydrothermal systems; (iv) the reported time lags between large-scale thermal anomalies and deformation; and (v) the breathing of potentially destructive volcanic calderas. The “stop-and-go” mechanism provides a consistent and realistic framework to link multiple geophysical and geochemical signals with the processes leading to volcanic unrest and eruption.

volcanoes

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport

New dimensions in satellite hydrology

Consideration of the use of remote sensing technology applied from satellites to obtain information for the rapid and continuing assessment of the hydrologic cycle. A detailed account is given of the hydrological information made available through the activities of the ERTS-1 satellite, an experimental satellite entirely devoted to earth resources observations, and the NOAA-2 satellite, a high-resolution operational environmental satellite. Following a description of the satellites and their payloads, it is shown how with their aid information can be obtained regarding atmospheric moisture, surface water and snow cover, glaciers, potential flood situations, and subsurface water fluctuations. In addition, the use of the ERTS-1 and NOAA-2 satellites in watershed characterization and modeling and in monitoring water quality is discussed.

Rango, A.

Climate forcing controls on carbon terrestrial fluxes during shale weathering

Climate influences near-surface biogeochemical processes and thereby determines the partitioning of carbon dioxide (CO 2 ) in shale, and yet the controls on carbon (C) weathering fluxes remain poorly constrained. Using a dataset that characterizes biogeochemical responses to climate forcing in shale regolith, we implement a numerical model that describes the effects of water infiltration events, gas exchange, and temperature fluctuations on soil respiration and mineral weathering at a seasonal timescale. Our modeling approach allows us to quantitatively disentangle the controls of transient climate forcing and biogeochemical mechanisms on C partitioning. We find that ~3% of soil CO 2 (1.02 mol C/m 2 /y) is exported to the subsurface during large infiltration events. Here, net atmospheric CO 2 drawdown primarily occurs during spring snowmelt, governs the aqueous C exports (61%), and exceeds the CO 2 flux generated by pyrite and petrogenic organic matter oxidation (~0.2 mol C/m 2 /y). We show that shale CO 2 consumption results from the temporal coupling between soil microbial respiration and carbonate weathering. This coupling is driven by the impacts of hydrologic fluctuations on fresh organic matter availability and CO 2 transport to the weathering front. Diffusion-limited transport of gases under transient hydrological conditions exerts an important control on CO 2 (g) egress patterns and thus must be considered when inferring soil CO 2 drawdown from the gas phase composition. Our findings emphasize the importance of seasonal climate forcing in shaping the net contribution of shale weathering to terrestrial C fluxes and suggest that warmer conditions could reduce the potential for shale weathering to act as a CO 2 sink.

58 GEOSCIENCES

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

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

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite

An Idea for an Active Seismic Experiment on Mars in 2008

The detection of liquid water is of prime interest and should have deep implications in the understanding of the Martian hydrological cycle and also in exobiology. In the frame of the 2007 joint CNES-NASA mission to Mars, a set of 4 NETLANDERS developed by an European consortium is expected to be launched in June 2007. We propose to use a second spacecraft going or landing to Mars to release near one of the Netlander a series of artificial metallic meteorites, in order to perform an active seismic experiment providing a seismic profile of the crust and subsurface.

Lognonne, Ph.

Terrestrial Water Storage

Terrestrial water storage can be defined as the summation of all water on the land surface and in the subsurface. It includes surface soil moisture, root zone soil moisture, groundwater, snow,ice, water stored in the vegetation, river and lake water. Terrestrial water storage (TWS) changes have been observed by the Gravity Recovery and Climate Experiment (GRACE) mission since 2002. GRACE has provided an unprecedented view of the terrestrial water storage variations at large scales. Extremes in water storage often are associated with droughts and flooding events because they are driven by the surplus or deficit of water. Few hydrologic observing networks yield sufficient data for comprehensive monitoring of changes in the total amount of water stored in a region. GRACE observations have helped to fill this gap. This book chapter is divided into 4 sections.Section 1 provides an overview of the GRACE mission. Section 2 reviews the terrestrial water storage solutions that are available. Section 3 describes regional to global applications of GRACE for monitoring extremes in the terrestrial water storage, including droughts and flooding events. Conclusions and future directions are reported in Section 4.

Terrestrial Water Storage

Water resources and hydrology of Mars

The surface of Mars has been extensively modified by a large variety of water erosional and depositional processes. Although liquid water is presently unstable on the planet's surface, in its cold, hyperarid climate, there is abundant geomorphological evidence of past fluvial valley development multiple episodes of catastrophic flooding, periglacial landforms, ice-related permafrost, lake deposits, eroded impact craters and possible glacial landforms throughout much of Mars' geological history. The amount of water required to form such features is estimated to be equivalent to a planet-wide layer approximately 50 meters deep. Some of this water undoubtedly was removed from the planet by atmospheric escape processes, but much probably remains in the subsurface of Mars. Jakosky summarized the present partitioning of water on Mars, expressed as an average global depth, as follows: in the polar caps, 30 meters; in the megaregolith, 500 to 1000 meters; structurally bound in clays, 10 meters; and in high latitude regolith, a few meters. However, most of this water is probably in the form of ice, except in anomalous areas of possible near surface liquid water, and in regions where hydrothermal systems are still active. The best locations for prospecting are those areas where water or ice is sufficiently concentrated at shallow enough depths to make it feasible to pump out or mine.

Baker, V. R.

AIRSAR Data for Geological and Geomorphological Mapping in the Great Sandy Desert and Pilbara Regions of Western Australia

Enhancements of AIRSAR data have demonstrated the benefits of synthetic aperture radar (SAR) for revealing an additional and mich higher level of information about the composition of the terrain than enhancements f either SPOT-PAN or Landsat TM data. With appropriate image processing techniques, surface and near surface geological structures, hydrological systems (both current and ancient) and landform features, have been evidenced in a diverse range of landscapes. In the Great Sandy Desert region where spectral variability is minimal, radar's sensitivity to the micromorphology of sparse exposures of subcrop and lag gravels has provided a new insight into the region's geological framework, its landforms, and their evolution. In the Pilbara region, advanced processing of AIRSAR data to unmix the backscatter between and within the three frequencies of data has highlighted subsurface extensions of greenstone lithologies below sand cover and morphological evidence of past flow conditions under former climate regimes. On the basis of these observations, it is recommend that radar remote sensing technology involving the use of high resolution, polarimetric data be seriously considered as a viable tool for exploration in erosional and depositional environments located within Australia's mineral and oil-prospective provinces.

Tapley, Ian J.

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY

Antecedent Hydrologic Conditions Reflected in Stream Lithium Isotope Ratios During Storms

Antecedent hydrological conditions are recorded through the evolution of dissolved lithium isotope signatures (δ 7 Li) by juxtaposing two storm events in an upland watershed subject to a Mediterranean climate. Discharge and δ 7 Li are negatively correlated in both events,but mean δ 7 Li ratios and associated ranges of variation are distinct between them. We apply a previously developed reactive transport model (RTM) for the site to these event-scale flow perturbations, but observed shifts in stream δ 7 Li are not reproduced. To reconcile the stability of the subsurface solute weathering profile with our observations of dynamic stream δ 7 Li signatures, we couple the RTM to a distribution of fluid transit times that evolve based on storm hydrographs. The approach guides appropriate flux-weighting of fluid from the RTM over a range of flow path lengths, or equivalently fluid residence times. This flux-weighted RTM approach accurately reproduces dynamic storm δ 7 Li-discharge patterns distinguished by the antecedent conditions of the watershed.

58 GEOSCIENCES

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

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic