Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

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 829 records · Page 46

Groundwater Remediation in a Flood plain Aquifer at Shiprock, New Mexico - 16097

The Shiprock, New Mexico, Disposal Site is the location of a former uranium and vanadium ore-processing mill that operated from 1954 to 1968. During milling years, tailings leachate and other mill-related liquids migrated into underlying sediments, resulting in the contamination of groundwater beneath the river terrace containing the mill and a nearby flood plain alluvial aquifer adjacent to the San Juan River. The US DOE Office of Legacy Management manages the site and is implementing the current groundwater remedy for the flood plain area, consisting of natural flushing supplemented by extraction of alluvial groundwater to accelerate reduction in contaminant concentrations. The alluvial aquifer contains elevated concentrations of sulfate, uranium, and nitrate derived largely from the historical eastward migration of groundwater through terrace alluvium and underlying bedrock to the flood plain. The bedrock formation beneath the entire Shiprock site is Mancos shale, which forms an escarpment that separates the terrace from the flood plain. Contamination in the flood plain alluvial aquifer is influenced by a dynamic flow system that changes seasonally and from year to year. Simulations with a flood plain flow model based on monitored groundwater levels and coincident river-surface elevations show that background groundwater flow, unaffected by remediation pumping, is characterized by distinct hydraulic conditions for three separate periods each year. During months of low river flow in winter and early spring, river losses to the aquifer create a large hypo-rheic zone in the south half of the flood plain. Simultaneously, recharge of continually flowing surface water emptying onto the north half of the flood plain from Bob Lee Wash creates divergent, radial flow from the wash outlet and diverts northwestward-migrating contaminant plumes to the northeast where they discharge to the river. Subsequent bank storage processes in May and June from high snow melt runoff change background flow directions, increasing the spread of contamination in the alluvium. From mid-summer through fall, river flows are typically similar in value to those observed prior to the onset of snow melt, but groundwater elevations in the aquifer are noticeably lower than water levels observed in winter and early spring and may induce river losses to the aquifer. Several different contaminant sources have impacted the alluvial aquifer. Much of the contamination originated as tailings leachate and raffinate wastewater that infiltrated the terrace subsurface and eventually reached the flood plain at the escarpment. Additional sources of contamination in the aquifer included discharge of mill effluent to Bob Lee Wash and a pond at the base of the escarpment, and aqueous mobilization of solid-phase contamination in the form of windblown deposits and contaminants adsorbed to alluvial sediments. The degree to which remnant, secondary-source processes, such as leaching of solid-phase deposits and discharge of terrace groundwater across the escarpment, contribute to existing contaminant plumes is uncertain. The flood plain remediation system consists of two groundwater extraction wells, pumping from two horizontal wells placed in trenches near the base of the escarpment, a sump collecting discharges from seeps in the escarpment wall, and discharge of collected water to an evaporation pond on the terrace. Flow-model simulations indicate that most contaminated areas are captured by remediation pumping, resulting in several zones of freshwater between the river and the escarpment. Despite the highly transient nature of groundwater flow beneath the flood plain and the varied contaminant sources, initiation of remediation pumping in 2003 and coordinated operation of the system's various components have successfully removed large amounts of contaminant mass and helped prevent discharge of contaminated groundwater to the river. The percentage reductions of dissolved sulfate, uranium, and nitrate mass in the aquifer estimated for the first 9 years of remediation range from 40% to 72%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Boat-Dragged Riverbed and Mid-Column Water Quality along the Hanford Reach in the Columbia River Shorelines

This dataset contains water quality measurements used to identify sites of hydrologic exchange (i.e.hydrologic exchange flows; HEFs) in the Columbia River near Richland, WA by looking for anomalies in temperature and electrical conductivity (EC) along the bed of the river in February, June, July, August, and November 2017. These are ultimately the ‘ends’ of HEFs as they are locations of subsurface inflow to the river. This dataset contains dissolved oxygen (DO) measurements at 1-min intervals; temperature, EC, depth, and GPS coordinates at 1-second intervals; and N-NO3 at 2-min intervals along the river bed and mid-column water. DO measurements are included in this dataset for completeness but were not used in the analysis of HEFs due to the relatively long time interval of collection.

54 ENVIRONMENTAL SCIENCES↗

An Integrated Feasibility Study of Reservoir Thermal Energy Storage (RTES) in Portland, OR, USA

In regions with long cold overcast winters and sunny summers, Deep Direct-Use (DDU) can be coupled with Reservoir Thermal Energy Storage (RTES) technology to take advantage of pre-existing subsurface permeability to save summer heat for later use during cold seasons. Many aquifers worldwide are underlain by permeable regions (reservoirs) containing brackish or saline groundwater that has limited beneficial use due to poor water quality. We investigate the utility of these relatively deep, slow flowing reservoirs for RTES by conducting an integrated feasibility study in the Portland Basin, Oregon, USA, developing methods and obtaining results that can be widely applied to RTES systems elsewhere. As a case study, we have conducted an economic and social cost-benefit analysis for the Oregon Health and Science University (OHSU), a teaching hospital that is recognized as critical infrastructure in the Portland Metropolitan Area. Our investigation covers key factors that influence feasibility including 1) the geologic framework, 2) heat and fluid flow modeling, 3) capital and maintenance costs, 4) the regulatory framework, and 5) operational risks. By pairing a model of building seasonal heat demand with an integrated model of RTES resource supply, we determine that the most important factors that influence RTES efficacy in the study area are operational schedule, well spacing, the amount of summer heat stored (in our model, a function of solar array size), and longevity of the system. Generally, heat recovery efficiency increases as the reservoir and surrounding rocks warm, making RTES more economical with time. Selecting a base-case scenario, we estimate a levelized cost of heat (LCOH) to compare with other sources of heating available to OHSU and find that it is comparable to unsubsidized solar and nuclear, but more expensive than natural gas. Additional benefits of RTES include energy resiliency in the event that conventional energy supplies are disrupted (e.g., natural disaster) and a reduction in fossil fuel consumption resulting in a smaller carbon footprint. Key risks include reservoir heterogeneity and a possible reduction in permeability through time due to scaling (mineral precipitation). Lastly, a map of thermal energy storage capacity for the Portland Basin yields a total of 43,400 GWh, suggesting tremendous potential for RTES in the Portland Metropolitan Area.

14 SOLAR ENERGY↗

3D Geologic Framework Modelling of the Los Alamos National Laboratory Site and Pajarito Plateau: Integrating a realistic 3D fault network and modelling subsurface relationships in a sparsely sampled and complex geologic region

The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.

58 GEOSCIENCES↗

Thermal Characterization and Exhumation of Northwest San Juan Basin Area, NM

The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.

02 PETROLEUM↗

Evaluation of Stored Reducing Capacity to Optimize Operations and Maintenance for In Situ Biological Treatment of Cr(VI) in Groundwater - 18534

Injection of organic carbon substrates to stimulate in situ biologically mediated reduction of dissolved hexavalent chromium (Cr[VI]) from groundwater by reduction of this oxidized form to the essentially insoluble trivalent form (Cr[III]) is an established remedial technology. The key to successful implementation of in situ biological Cr(VI) reduction is the distribution of electrons from the oxidation of the organic carbon substrates throughout the saturated subsurface to establish robust In Situ Reactive Zones (IRZs). Once IRZs are established, reduction of Cr(VI) is often observed to persist over relatively long time frames (months to years) after the injected organic carbon substrates are degraded or flux of non-amended water flows through the IRZ. The longevity of Cr(VI) treatment is attributed to the establishment of residual-reducing capacity created by the storage of electrons in ferrous iron and sulfide-bearing minerals present in the aquifer matrix, and microbial biomass generated by the injection of organic carbon substrates. In this paper, more than 10 years of data collected from operations and maintenance of IRZ systems for remediation of a large-scale Cr(VI) plume were evaluated to examine the factors that influence the longevity of residual-reducing capacity. The durations that residual-reducing capacity has resulted in sustained Cr(VI) treatment were variable, ranging from months to more than 6 years between carbon amendment events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Detecting dynamic causal inference in nonlinear two-phase fracture flow

Identifying dynamic causal inference involved in flow and transport processes in complex fractured-porous media is generally a challenging task, because nonlinear and chaotic variables may be positively coupled or correlated for some periods of time, but can then become spontaneously decoupled or non-correlated. In his 2002 paper (Faybishenko, 2002), the author performed a nonlinear dynamical and chaotic analysis of time-series data obtained from the fracture flow experiment conducted by Persoff and Pruess (1995), and, based on the visual examination of time series data, hypothesized that the observed pressure oscillations at both inlet and outlet edges of the fracture result from a superposition of both forward and return waves of pressure propagation through the fracture. In the current paper, the author explores an application of a combination of methods for detecting nonlinear chaotic dynamics behavior along with the multivariate Granger Causality (G-causality) time series test. Based on the G-causality test, the author infers that his hypothesis is correct, and presents a causation loop diagram of the spatial-temporal distribution of gas, liquid, and capillary pressures measured at the inlet and outlet of the fracture. The causal modeling approach can be used for the analysis of other hydrological processes, for example, infiltration and pumping tests in heterogeneous subsurface media, and climatic processes, for example, to find correlations between various meteorological parameters, such as temperature, solar radiation, barometric pressure, etc.

54 ENVIRONMENTAL SCIENCES↗

Software for Support of Groundwater Contaminant Fate and Transport Analysis – 13345

Groundwater contamination by hexavalent chromium and other nuclear reactor operation-related contaminants has resulted in the need for groundwater remedial actions within the Hanford Site reactor areas (the Hanford Site 100 Area). The large geographic extent of the resultant contaminant plumes requires an extensive level of understanding of the aquifer structure, characteristics, and configuration to support assessment and design of remedial alternatives within the former 100-D, 100-H, and 100-K reactor areas. The authors have prepared two- and three-dimensional depictions of the key subsurface geologic structures at two Hanford Site reactor operable units (100-K and 100-D/H). These depictions, prepared using commercial-off-the-shelf (COTS) visualization software, provide a basis for expanding the understanding of groundwater contaminant migration pathways, including identification of geologically defined preferential groundwater flow pathways. These identified preferential flow pathways support the conceptual site model and help explain both historical and current contaminant distribution and transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗

Earth Source Heat: A Cascaded Systems Approach to DDU of Geothermal Energy on the Cornell Campus

This report documents the completion of the DOE study DE-EE0008103: Earth Source Heat: A Cascaded Systems Approach to DDU of Geothermal Energy on the Cornell Campus. This study involved a comprehensive evaluation of the potential for Earth Source Heat (ESH), Cornell’s specific application of Deep Direct Use (DDU) geothermal energy, to create a viable renewable source of thermal energy for its Ithaca, NY campus district heating system. The study included stochastic modeling that married each of two potential subsurface resources (two specific geological reservoirs) to documented campus heating operations based on hourly campus heating profiles and proposed integrated equipment controls (variable speed pumps and heat pumps). The outputs of the subsurface modeling were used as inputs for a Cornell-developed surface heat use Excel-based modeling program (the program was named “Modeled ENergy Use”, or “MEnU”). This allowed Cornell to compute realistic annual totals of extractable and useable heat for the campus, including details like heat losses, minimum summer demand, and controllable building temperature settings. The Cornell team then evaluated the costs and benefits of each modeled scenario. The primary evaluation was a standard “single bottom line” economic calculation for the Levelized Cost of Heat (LCOH). This represents the annualized cost (per unit heat) for Cornell. Additional valuations were also computed based on the benefit to the environment and to the regional economy (i.e., economic benefits external to campus). The results of this study demonstrate that if suitable reservoir flow can be attained, Earth Source Heat would be a viable technology to supply Cornell’s district heating system. This viability is based on the following criteria: (1) The modeled system produced a total useable heat output that exceeded the minimum annual campus heat load determined at the onset of the project (i.e., offsetting at least 20% of the annual campus thermal load). The output range is substantial depending on both the subsurface resource and the surface applications; at the high end, our modeled solutions produced up to ~70% of existing campus heat load with a single well pair producing 70 kg/s with integrated high-temperature heat pumps, and (2) Modeling demonstrated that multiple geothermal reservoirs could provide economically viable results when heat pumps were strategically integrated, a design that has been proven in at least one European installation already. In this context, “economically viable” means that the LCOH for the project is less than the regional commercial price of heat that would be generated using natural gas (the most common and cheapest fossil alternative in our area). Based on a single-bottom-line (LCOH) valuation, our results show that economic success is likely for any reservoir deeper than about 2.25 km where a sufficient permeability exists (or can be created with stimulation or similar engineering enhancement) and reasonable production temperatures sustained over time. The LCOH values of $\$ 4.77$ to $\$ 6.46$ per MMBtu are all lower than the comparable cost of providing heat energy with a natural gas boiler for a facility in the Northeast U.S. that obtains natural gas at commercial rates (Commercial heat costs would include at least $\$ 9.41$ in fuel costs, representing total boiler/transfer efficiency of 85% using gas at $\$ 8.00$ per MMBtu; LCOH would be higher if capital equipment and other operating costs were included). In addition to LCOH, Cornell’s scope also included evaluation of the Environmental value and Regional Economic value of each scenario. When considering the environmental and regional economic value of such a project, the overall “triple bottom line” cost (i.e., the LCOH less the LCOH ENV and LCOH REG values) is even lower and in some cases negative; i.e., in some cases a project achieving the listed performance would have higher value to the environment and to the regional economy than the total project costs over the 30-year project timeframe.

15 GEOTHERMAL ENERGY↗

A conceptual geochemical model of the geothermal system at Surprise Valley, CA

Characterizing the geothermal system at Surprise Valley (SV), northeastern California, is important for determining the sustainability of the energy resource, and mitigating hazards associated with hydrothermal eruptions that last occurred in 1951. Previous geochemical studies of the area attempted to reconcile different hot spring compositions on the western and eastern sides of the valley using scenarios of dilution, equilibration at low temperatures, surface evaporation, and differences in rock type along flow paths. These models were primarily supported using classical geothermometry methods, and generally assumed that fluids in the Lake City mud volcano area on the western side of the valley best reflect the composition of a deep geothermal fluid. In this contribution, we address controls on hot spring compositions using a different suite of geochemical tools, including optimized multicomponent geochemistry (GeoT) models, hot spring fluid major and trace element measurements, mineralogical observations, and stable isotope measurements of hot spring fluids and precipitated carbonates. We synthesize the results into a conceptual geochemical model of the Surprise Valley geothermal system, and show that high-temperature (quartz, Na/K, Na/K/Ca) classical geothermometers fail to predict maximum subsurface temperatures because fluids re-equilibrated at progressively lower temperatures during outflow, including in the Lake City area. We propose a model where hot spring fluids originate as a mixture between a deep thermal brine and modern meteoric fluids, with a seasonally variable mixing ratio. The deep brine has deuterium values at least 3 to 4‰ lighter than any known groundwater or high-elevation snow previously measured in and adjacent to SV, suggesting it was recharged during the Pleistocene when meteoric fluids had lower deuterium values. The deuterium values and compositional characteristics of the deep brine have only been identified in thermal springs and groundwater samples collected in proximity to structures that transmit thermal fluids, suggesting the brine may be thermal in nature. On the western side of the valley at the Lake City mud volcano, the deep brine-meteoric water mixture subsequently boils in the shallow subsurface, precipitates calcite, and re-equilibrates at about 130 °C. On the eastern side of the valley, meteoric fluid mixes to a greater extent with the deep brine, cools conductively without boiling, and the composition is modified as dissolved elements are sequestered by secondary minerals that form along the cooling and outflow path at temperatures <130 °C. Re-equilibration of geothermal fluids at lower temperatures during outflow explains why subsurface temperature estimates based on classical geothermometry methods are highly variable, and fail to agree with temperature estimates based on dissolved sulfate-oxygen isotopes and results of classical and multicomponent geothermometry applied to reconstructed deep well fluids. Furthermore, the proposed model is compatible with the idea suggested by others that thermal fluids on the western and eastern side of the valley have a common source, and supports the hypothesis that low temperature re-equilibration during west to east flow is the major control on hot spring fluid compositions, rather than dilution, evaporation, or differences in rock type.

15 GEOTHERMAL ENERGY↗

Deep Unsaturated Zone Contributions to Carbon Cycling in Semiarid Environments

Understanding terrestrial carbon cycling has relied primarily on studies of topsoils that are typically characterized to depths shallower than 0.5 m. At a semiarid site instrumented down to 7 m, here we measured seasonal- and depth-resolved carbon inventories and fluxes and groundwater and unsaturated zone flow rates. Measurements showed that ~30% of the CO 2 efflux to the atmosphere (60% in winter) originates from below 1 m, contrary to predictions of less than 1% by Earth System Model land modules. Respiration from deeper roots and deeper microbial communities is supported by favorable subsurface temperatures, moisture, and oxygen availability. Below 1 m, dissolved organic carbon fluxes from the overlying soil and C from deep roots and exudates are expected to be important in sustaining microbial respiration. Finally, because these conditions are characteristic of semiarid climate regions, we contend that Earth System Model land modules should incorporate such deeper soil processes to improve CO 2 flux predictions.

54 ENVIRONMENTAL SCIENCES↗

Biostimulation induces syntrophic interactions that impact C, S and N cycling in a sediment microbial community

Stimulation of subsurface microorganisms to induce reductive immobilization of metals is a promising approach for bioremediation, yet the overall microbial community response is typically poorly understood. Here we used community proteogenomics to test the hypothesis that excess input of acetate activates syntrophic interactions among autotrophs and heterotrophs. A flow-through sediment column was incubated in a groundwater well of an acetate-amended aquifer. Genomic sequences from the community recovered during microbial sulfate reduction were used to econstruct, de novo, near-complete genomes for Desulfobacter (Deltaproteobacteria) and relatives of Sulfurovum and Sulfurimonas (Epsilonproteobacteria), and Bacteroidetes. Partial genomes were obtained for Clostridiales (Firmicutes) and Desulfuromonadales-like Deltaproteobacteria. The majority of proteins identified by mass spectrometry corresponded to Desulfobacter-like species, and demonstrate the role of this organism in sulfate reduction (Dsr and APS), nitrogen-fixation (Nif) and acetate oxidation to CO2 during amendment. Results suggest less abundant Desulfuromonadales and Bacteroidetes also actively contributed to CO2 production via the TCA cycle. Proteomic data indicate that sulfide was partially re-oxidized by Epsilonproteobacteria through nitrate-dependent sulfide oxidation (using Nap, Nir, Nos, SQR and Sox), with CO2 fixed using the reverse TCA cycle. Modeling shows that this reaction was thermodynamically possible, and kinetically favorable relative to acetate-dependent denitrification. We conclude that high-levels of carbon amendment aimed to stimulate anaerobic heterotrophy led to carbon fixation in co-dependent chemoautotrophs. These results have implications for understanding complex ecosystem behavior, and show that high levels of organic carbon supplementation can expand the range of microbial functionalities accessible for ecosystem manipulation.

metagenomics/proteomics/subsurface/sediment/syntro↗

2012 Groundwater Monitoring and Inspection Report Gnome-Coach, New Mexico, Site

Gnome-Coach was the site of a 3-kiloton underground nuclear test conducted in 1961. Surface and subsurface contamination resulted from the underground nuclear testing, post-test drilling, and a groundwater tracer test performed at the site. Surface reclamation and remediation began after the underground testing. A Completion Report was prepared, and the State of New Mexico is currently proceeding with a conditional certificate of completion for the surface. Subsurface corrective action activities began in 1972 and have generally consisted of annual sampling and monitoring of wells near the site. In 2008, the annual site inspections were refined to include hydraulic head monitoring and collection of samples from groundwater monitoring wells onsite using the low-flow sampling method. These activities were conducted during this monitoring period on January 18, 2012. Analytical results from this sampling event indicate that concentrations of tritium, strontium-90, and cesium-137 were generally consistent with concentrations from historical sampling events. The exceptions are the decreases in concentrations of strontium-90 in samples from wells USGS-4 and USGS-8, which were more than 2.5 times lower than last year's results. Well USGS-1 provides water for livestock belonging to area ranchers, and a dedicated submersible pump cycles on and off to maintain a constant volume in a nearby water tank. Water levels in wells USGS-4 and USGS-8 respond to the on/off cycling of the water supply pumping from well USGS-1. Well LRL-7 was not sampled in January, and water levels were still increasing when the transducer data were downloaded in September. A seismic reflection survey was also conducted this year. The survey acquired approximately 13.9 miles of seismic reflection data along 7 profiles on and near the site. These activities were conducted from February 23 through March 10, 2012. The site roads, monitoring well heads, and the monument at surface ground zero were in good condition at the time of the site inspection. However, it was reported in September 2012 that the USGS-1 well head had been damaged by a water truck in April 2012.

54 ENVIRONMENTAL SCIENCES↗

Final Technical Report for Department of Energy Award DE-SC0010857, “Towards parameterization of root-rock hydrologic interactions in the Earth System Model”

The goal of the project is to understand how plants survive droughts, as the decimation of transpiration could shift the surface energy balance from latent heat to sensible heat, leading to warming of the lower atmosphere and amplification of drought. The hypothesis we investigated is that there is a store of moisture in the weathered bedrock below the organic soil mantle, so that plants that have roots deep enough to access this moisture reservoir could survive drought. We developed a new stochastic parameterization of hydraulic conductivity that introduces heterogeneity into the traditional formulation and captures the preferential flow through the weathered bedrock (Vrettas and Fung, 2015). With the new parameterization in the Richards equation, we succeeded in reproducing the fluctuations of the water table in seven well locations over six years. We also and carried out a series of model experiments that explore how subsurface properties impact evapotranspiration (ET) in a Mediterranean climate where a significant portion of ET is observed to take place in the dry and sunny summer when the precipitation is insufficient to meet the demand. Our results show that hydraulic redistribution is important for sustaining ET in the dry seasons when the vertical gradient in water potential is large. The results highlight the importance of lithology, species composition and root function for ET, especially under dry conditions.

58 GEOSCIENCES↗

Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation

Long-term seismic monitoring of carbon capture and storage projects is needed to verify that the injected gas is safely stored in the subsurface until permanence can be assured. Conventional surface seismic monitoring techniques are usually expensive, require highly invasive surface operations, and need significant time investments on the part of personnel for both the field effort and processing the acquired data. For these reasons, permanent reservoir monitoring technologies are preferred, as they can offer a cost-effective solution for long-term monitoring. As part of the monitoring program of the Archer Daniels Midland’s large-scale injection of CO 2 in Decatur, Illinois, USA, a continuous seismic monitoring array was installed using a combination of surface orbital vibrator (SOV) sources and fiber-optic cables for distributed acoustic sensing (DAS) acquisition with the objective to build a continuous monitoring array. The aim of the presented project was to build a monitoring array and platform that integrates real-time seismic data with conventional data streams and provides continuous data analysis using dynamic computational models to deliver a comprehensive real-time assessment of subsurface conditions. It is in this context that the Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation project was proposed with the objective to develop an integrated architecture that utilizes a permanent seismic monitoring network, combines the real-time geophysical and process data with reservoir flow and geomechanical models to create a comprehensive monitoring, visualization, and control system that delivers critical information for process surveillance and optimization.

54 ENVIRONMENTAL SCIENCES↗

Two- and Three-Dimensional Depiction of Subsurface Geology Using Commercial Software for Support of Groundwater Contaminant Fate and Transport Analysis - 13345

Groundwater contamination by hexavalent chromium and other nuclear reactor operation-related contaminants has resulted in the need for groundwater remedial actions within the Hanford Site reactor areas (the Hanford Site 100 Area). The large geographic extent of the resultant contaminant plumes requires an extensive level of understanding of the aquifer structure, characteristics, and configuration to support assessment and design of remedial alternatives within the former 100-D, 100-H, and 100-K reactor areas. The authors have prepared two- and three-dimensional depictions of the key subsurface geologic structures at two Hanford Site reactor operable units (100-K and 100-D/H). These depictions, prepared using commercial-off-the-shelf (COTS) visualization software, provide a basis for expanding the understanding of groundwater contaminant migration pathways, including identification of geologically-defined preferential groundwater flow pathways. These identified preferential flow pathways support the conceptual site model and help explain both historical and current contaminant distribution and transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Inferring subsurface heterogeneity from push-drift tracer tests

In this study, we consider the late-time tailing in a tracer test performed with a push-drift methodology (i.e., quasi-radial injection followed by drift under natural gradient). Numerical simulations of such tests are performed on 1000 multi-Gaussian 2-D log-hydraulic conductivity field realizations of varying heterogeneity, each under eight distinct mean flow directions. The ensemble pdfs of solute return times are found to exhibit power law tails for each considered variance of the log-hydraulic conductivity field, σ$2\atop{ln K}$. The tail exponent is found to relate straightforwardly to σ$2\atop{ln K}$ and, within the parameter space we explored, to be independent of push-phase pumping rate, pumping duration, and local-scale dispersivity. We conjecture that individual push-drift tracer tests in wells with screened intervals much greater than the vertical correlation length of the aquifer will exhibit quasi-ergodicity and that their tail exponent may be used to infer σ$2\atop{ln K}$. Lastly, we calibrate a predictive relationship of this sort from our Monte Carlo study, and apply it to data from a push-drift test performed at a site of approximately known heterogeneity—closely matching the existing best estimate of heterogeneity.

58 GEOSCIENCES↗