Search NASA⌕ Search

SEARCH · Search NASA

Results for “Aquifer characterization”

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.

124 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sampling in Long-Screened Wells: Issues, Misconceptions, and Solutions

The issues associated with long-screened wells (LSWs) (and open boreholes) at contaminated sites are well documented in the groundwater literature but are still not fully appreciated in practice. As established in seminal and review papers going back over three decades, the interpretation of sampling results from LSWs is challenging in the presence of vertical hydraulic gradients and borehole flow; furthermore, LSWs allow for vertical redistribution of contamination between aquifer layers. Acknowledgment of these issues has led to the development of new technologies and well designs to enable discrete-zone monitoring (DZM), yet LSWs remain common for many reasons, for example, as multipurpose wells, for geophysical logging, and (or) as legacy installations. Despite the literature on LSWs and despite the adoption of DZM at many sites, the use of LSWs persists and the challenges of interpreting sampling results from LSWs remain. In this issue paper, we provide a conceptual overview of the problems posed by LSWs and review existing literature and past work to improve the interpretation of sampling in LSWs. We draw on experience from previous studies at the Hanford Site in eastern WA, USA, and use synthetic examples to illustrate key concepts and challenges for interpretation. A recently published analytical modeling framework is used to develop illustrative synthetic examples and demonstrate a workflow for building scientific intuition to understand issues around interpreting samples from LSWs, which is critical to effective characterization and groundwater remediation at sites with LSWs.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

CMIP6 Historical Simulations (1850‐2014) with GISS‐E2.1

Simulations of the CMIP6 historical period 1850‐‐2014, characterized by the emergence of anthropogenic climate drivers like greenhouse gases, are presented for different configurations of the NASA Goddard Institute for Space Studies (GISS) Earth System ModelE2.1. The GISS‐E2.1 ensembles are more sensitive to forcing than their CMIP5 predecessors (GISS‐E2), but warm less during recent decades due to reduced total forcing. This forcing reduction is attributed to an increase in longwave opacity in pre‐industrial simulations, resulting in an atmosphere less sensitive to further increases in opacity that result from rising greenhouse gas concentrations. This demonstrates the importance of the base climatology to forcing and forced climate trends. Most model versions match observed temperature trends since 1979 from the ocean to the stratosphere. The choice of ocean model is important to the transient climate response, as found in CMIP5 GISS‐E2: the model that more efficiently exports heat to the deep ocean shows a smaller rise in tropospheric temperature. Model sea level rise over the historical period is traced to excessive drawdown of aquifers to meet irrigation demand with a smaller contribution from thermal expansion. This shows how fully coupled models can provide indirect observational constraints upon forcing, in this case, constraining irrigation rates with observed sea level changes. The overall agreement of GISS‐E2.1 with observed trends is familiar from evaluation of its predecessors, as is the conclusion that these trends are almost entirely anthropogenic in origin.

NASA Goddard Institute for Space Studies (GISS) Ea↗
NASA NTRS1998

Gale Crater: An Amazonian Impact Crater Lake at the Plateau/Plain Boundary

Gale is a 140-km diameter impact crater located at the plateau/plain boundary in the Aeolis Northeast subquadrangle of Mars (5S/223W). The crater is bordered in the northward direction by the Elysium Basin, and in eastward direction by Hesperian channels and the Aeolis Mensae 2. The crater displays a rim with two distinct erosion stages: (a) though eroded, the south rim of Gale has an apparent crest line visible from the north to the southwest (b) the west and northwest rims are characterized by a strong erosion that, in some places, partially destroyed the rampart, leaving remnant pits embayed in smooth-like deposits. The same type of deposits is observed north, outside Gale, it also borders the Aeolis Mensae, covers the bottom of the plateau scarp, and the crater floor. The central part of Gale shows a 6400 km2 subround and asymmetrical deposit: (a) the south part is composed of smooth material, (b) the north part shows spectacular terraces, streamlines, and channels. The transition between the two parts of the deposit is characterized by a scarp ranging from 200 to 2000 in high. The highest point of the scarp is at the center of the crater, and probably corresponds to a central peak. Gale crater does not show a major channel directly inflowing. However, several large fluvi systems are bordering the crater, and could be at the origin of the flooding of the crater, or have contributed to. One fluvial system is entering the crater by the southwest rim but cannot be accounted alone for the volume of sediment deposited in the crater. This channel erodes the crater floor deposit, and ends in a irregular-shaped and dark albedo feature. Gale crater shows the morphology of a crater filled during sedimentation episodes, and then eroded Part of the lower sediment deposition contained in Gale might be ancient and not only aqueous in origin. According to the regional geologic history, the sedimentary deposit could be a mixture of aeolian and pyroclastic material, and aqueous sedimentary material that can originate both from drainage of the regional subsurface aquifer, and/or from surface flood. The central deposit shows three main levels: (a) the current crater floor (north of Gale), (b) an ancient level about 200 rn higher (south of Gale), and (c) the massive terraced deposits. A crater statistics on the 15,400 kM2 area of the crater floor and deposit [3,41 gave: 259+/-112.4 craters, most of them partly embayed in the sedimentary deposit, and all inferior to 5-km diameter. For superimposed crater population only, the result is 194+/-112. The deduced relative ages ranges from Early to Middle Amazonian. The population of craters are comparable for the three levels, implying that the last sedimentation/erosion episode on Gale was recent and affected the whole crater. The streamlined morphology of the border of the deposit, the layering, the channels, and the terraces are compatible with a significant fluvio-lacustrine history of the site. Multiple levels may suggest different episodes, but the common statistical age of the three levels shows that the last episode involved the whole crater. The origin of the lake water in Gale may have varied in time. Three major contributions have been proposed: (a) the drainage of the regional underground aquifer by Gale crater over an area of 110-km radius around the crater which would have provided approximately 1,600 cubic km of water, (b), surface drainage entering Gale by the south and north rims. In the south, a 250-km long system originates in the cratered uplands in a Noachian crater material plain (Nc), and crosses Hesperian and Amazonian crater material plains (AHc) northward [1]. Several fluvial systems originate in the Aeolis Mensae, east of Gale. They may had two functions in time: to recharge, the underground aquifer in the region of Gale, and to supply surface water in the crater by overspilling the northern rim, and (c) surface floods that originated from the rising of the water level in the Elysium Basin. According to the Amazonian age of Gale's floor, and the erosion direction in the crater, a flood from Elysium Basin is the most likely event to explain the material observed in Gale, and the formation of the last lake. This last flood may have been important enough to flood the central deposit up to about 1400 m above the crater floor, leaving two islands (non stream lined features) at the center of the deposit. Terrace spacing suggests a regular drop of the lake level in time. Fractures in terraces perpendicular to the shoreline can be interpreted either as: (a) the result of the drainage systems during the waning of the lake, or (b) traces of the pressure of an ice-covered sheet associated with subglacial drainage. The presence of a lake of such volume during the Amazonian period is one more evidence that water was still active on Mars relatively recently. Gale crater offers the rare opportunity to unveil a key-period of the martian history. The Amazonian might proved not as cold and dry as previously thought. The presence of large lakes and basins (Elysium Basin is large as the Mediterranean Sea), reinforces the model of an extensive water activity during the Amazonian that has still to be understood in the context of an assumed cooling and drying planet. The sediments and rocks that were left of this period in Gale keep the record of the climatic conditions of the Amazonian and the clues that are missing to understand the climatic evolution of Mars. In addition, Gale crater presents the advantage to be located at the plateau/plain boundary, which has never been studied and contains information about the two main martian geological units. As a conclusion, we propose a table that summarizes the worthiness of a mission in Gale crater, and the expected science return relative to the objectives to be met by the Surveyor Program. Additional information contained in the original.

Cabrol, N. A.↗
DOE OSTI2021

The Co-Transport of PFAS and Cr(VI) in porous media

PFAS and Cr are present at some sites as co-contaminants. Here, the objective of this research was to investigate the co-transport behavior of per- and polyfluoroalkyl substances (PFAS) and hexavalent chromium (Cr(VI)) in porous media. Miscible-displacement experiments were conducted using two soils and an aquifer sediment with different geochemical properties. Perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) were employed as model PFAS. The retardation of PFOS was decreased in the presence of Cr(VI). Conversely, the transport and retardation of PFOA was not affected by the presence of Cr(VI). The reduction of PFOS retardation caused by Cr(VI) is likely due to sorption competition for both organic-carbon and inorganic (metal-oxides and clay minerals) domains. The relative contributions of the three soil constituents to PFOS sorption and the potential for competition between PFOS and Cr(VI) is a function of the geochemical composition of the porous media (i.e., organic carbon, metal-oxides and clay minerals). The PFAS had minimal impact on the retention and transport of Cr(VI). To our knowledge, the results presented herein represent the first reported data for PFOS and Cr(VI) co-transport in porous media. The results of this study indicate that the presence of Cr(VI) has the potential to increase the migration potential of PFOS in soil and groundwater, which should be considered when characterizing electroplating facilities, leather tanning facilities, and other co-contaminated sites.

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

SUBTASK 1.6 – BASIN ELECTRIC CARBON STORAGE RESEARCH PROJECT: NOVEL MONITORING TECHNIQUES

The Energy & Environmental Research Center (EERC) conducted baseline activities associated with an applied research project at Basin Electric Power Cooperative’s (Basin’s) carbon capture and storage (CCS) site in Beulah, North Dakota, to establish novel carbon storage-monitoring techniques as commercial methods under Cooperative Agreement No. DE-FE0024233, Subtask 1.6. The following report summarizes the baseline activities performed and briefly describes the subsequent (operational monitoring) activities that have been proposed to the U.S. Department of Energy (DOE) as part of the overall project to develop and demonstrate novel monitoring techniques at North America’s largest permitted CCS operation. Dakota Gasification Company (DGC), a wholly owned subsidiary of Basin, owns and operates the Great Plains Synfuels Plant (GPSP) approximately 5 miles northwest of the town of Beulah, North Dakota (Figure 1). In 2023, DGC received approval from the North Dakota Industrial Commission (NDIC) to develop a storage facility on-site for injecting a stream of carbon dioxide (CO2) captured from GPSP. DGC will transport the captured CO2 stream with approximately 6.8 miles of transmission lines that extend north of GPSP and inject >1 million tonnes (MMt) of CO2 annually (>1 MMt/yr) over a 12-year period with up to six underground injection control (UIC) Class VI-compliant injection wells completed in the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer underlying GPSP. The Broom Creek Formation lies approximately 5900 feet (ft) below ground surface (bgs) at GPSP. The commercial scale (i.e., >1 MMt/yr) of DGC’s permitted carbon storage project is ideal for developing and testing the novel monitoring techniques included within Subtask 1.6. The goals of this project are to demonstrate 1) the cost-effectiveness of novel monitoring technologies included as part of this research, 2) technology capability for tracking the CO2 plume and/or associated pressure response in the subsurface and monitoring out-of-zone migration, and 3) compliance with UIC Class VI program requirements. The research activities proposed for the overall project include 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; 5) advanced wellbore-monitoring methods; 6) deployment of an AIM monitoring network; 7) EM monitoring of CO2 with real-time data processing; 8) continued seasonal drone-based surveillance studies; 9) seismic monitoring with passive and active surveys; and 10) wellbore monitoring with nuclear magnetic resonance (NMR) for near-surface characterization. Completion of Activities 1.0–5.0 (baseline activities) are described in this report. Upon authorization of funding by DOE, the EERC will initiate Activities 6.0– 10.0 (operational monitoring activities). Current state-of-the-art (SOA) carbon storage-monitoring techniques require countless labor hours dedicated to the acquisition of data. Once data are gathered, these SOA techniques often rely on commercial facilities to process raw data from the field. However, it is anticipated that next-generation monitoring techniques, such as those being demonstrated, will lower acquisition footprints, be less operationally intensive, and improve data acquisition efficiencies. These new techniques are more conducive to the application of machine learning, artificial intelligence, and automation, thus providing a pathway for integration into active control systems, informing site operability, and improving the integration of data for future CCS projects across the United States. Additionally, reclaimed and active mining lands are present within the project site, creating a unique opportunity to demonstrate the effectiveness of remote sensing and surface-based geophysics monitoring techniques at similar project sites that may include disturbed, unconsolidated, or actively excavated near-surface environments. The efforts included in the overall project will produce necessary designs, learnings, and data acquired during the baseline and operational monitoring periods that are necessary for time-lapse demonstration and validation of the described monitoring techniques. In addition, it is anticipated that the monitoring technologies included in this study will be compliant with UIC Class VI requirements to enable the potential for implementation at other CCS sites across the United States.

42 ENGINEERING↗
DOE OSTI2024

Workflow for Developing and Operating Subsurface Hydrogen Storage Facilities in Porous Reservoirs

Long-duration (seasonal) storage of natural gas (NG), which primarily consists of methane (CH 4 ), has been practiced for more than a hundred years at underground gas storage (UGS) facilities that use depleted hydrocarbon reservoirs, saline aquifers, and salt caverns. To enable hydrogen (H 2 ) to be used as a long-duration, energy-storage medium, similar facilities are envisioned for underground H 2 storage (UHS) of either H 2 or H 2 /NG mixtures. Experience with UGS can be used to guide recommended practices for developing and operating UHS facilities in porous reservoirs. The most important factors (formation/fluid properties and engineering choices) that influence the performance of UHS reservoirs have been identified and quantified in previous studies. These factors and choices influence phenomena that determine the sweep efficiency of the stored working gas. These phenomena include viscous fingering, hysteretic capillary trapping, and gravity override of the working gas, as well as the upconing of nonproductive fluid that determine the sweep efficiency of the stored working gas. This report describes initial recommended-practices and a project-development workflow for UHS facilities that utilize porous reservoirs, based on the current state-of-knowledge about H 2 behavior in the subsurface. The workflow sequentially addresses all aspects of UHS project development, including the identification of H 2 sources and users, site ranking and down-selection, geologic and reservoir-engineering characterization, reservoir design, testing, risk management, commissioning, operations, and monitoring for a UHS facility. The goal is to enable UHS facilities to be developed in an efficient and timely manner, while carefully managing project risks. This workflow is similar to that which has been developed for UGS facilities (see Figure 1 of API, 2022), with the addition of tasks and subtasks specific to H 2 and UHS. The project-development workflow is broken down into three major stages: (1) define the H 2 use case; (2) rank, down-select, and characterize potential, candidate UHS sites; and (3) reservoir design, integrity testing, risk assessment, commissioning, operations, and monitoring for selected UHS sites. Each major stage is further broken down into tasks and subtasks, which are described at a high level. This report also provides more detailed descriptions of all tasks and subtasks that involve reservoir analysis and testing.

08 HYDROGEN↗
NASA NTRS2013

The Origin of Carbon-bearing Volatiles in Surprise Valley Hot Springs in the Great Basin: Carbon Isotope aud Water Chemistry Characterizations

There are numerous hydrothermal fields within the Great Basin of North America, some of which have been exploited for geothermal resources. With methane and other carbon-bearing compounds being observed, in some cases with high concentrations, however, their origins and formation conditions remain unknown. Thus, studying hydrothermal springs in this area provides us an opportunity to expand our knowledge of subsurface (bio)chemical processes that generate organic compounds in hydrothermal systems, and aid in future development and exploration of potential energy resources as well. While isotope measurement has long been used for recognition of their origins, there are several secondary processes that may generate variations in isotopic compositions: oxidation, re-equilibration of methane and other alkanes with CO2, mixing with compounds of other sources, etc. Therefore, in addition to isotopic analysis, other evidence, including water chemistry and rock compositions, are necessary to identify volatile compounds of different sources. Surprise Valley Hot Springs (SVHS, 41º32'N, 120º5'W), located in a typical basin and range province valley in northeastern California, is a terrestrial hydrothermal spring system of the Great Basin. Previous geophysical studies indicated the presence of clay-rich volcanic and sedimentary rocks of Tertiary age beneath the lava flows in late Tertiary and Quaternary. Water and gas samples were collected for a variety of chemical and isotope composition analyses, including in-situ pH, alkalinity, conductivity, oxidation reduction potential (ORP), major and trace elements, and C and H isotope measurements. Fluids issuing from SVHS can be classified as Na-(Cl)-SO4 type, with the major cation and anion being Na+ and SO4 2-, respectively. Thermodynamic calculation using ORP and major element data indicated that sulfate is the most dominant sulfur species, which is consistent with anion analysis results. Aquifer temperatures at depth estimated by both dissolved SiO2 and Na-K-Ca geothermometers are in the range of 125.0 to 135.4 oC, and higher than the values measured at orifices (77.3 to 90.0 oC). CO2 and homologs of straight chain alkanes (C1-C5) were identified in gas samples. Carbon isotope values of alkanes increase with carbon numbers. The 13C fractionation between CO2 and dissolved inorganic carbon suggests they are out of carbon isotope equilibrium. The hypothesis regarding the formation of carbon-bearing compounds in SVHS may involve two processes: 1) Under high heat flow conditions which are caused by regional faulting and crustal extension, original high molecular weight organic compounds (kerogens) in clay-rich rocks decomposed to generate methane and other alkane homologs. 2) The SVHS area is associated with outflow structures, and distant from the heat source. Anaerobic oxidation of methane (AOM) with sulfate at shallow depth (< 90 oC) is suggested as being responsible for the generation of CO2 in SVHS.

Fu, Qi↗
NASA NTRS2013

The Origin of Carbon-bearing Volatiles in Surprise Valley Hot Springs in the Great Basin: Carbon Isotope and Water Chemistry Characterizations

There are numerous hydrothermal fields within the Great Basin of North America, some of which have been exploited for geothermal resources. With methane and other carbon-bearing compounds being observed, in some cases with high concentrations, however, their origins and formation conditions remain unknown. Thus, studying hydrothermal springs in this area provides us an opportunity to expand our knowledge of subsurface (bio)chemical processes that generate organic compounds in hydrothermal systems, and aid in future development and exploration of potential energy resources as well. While isotope measurement has long been used for recognition of their origins, there are several secondary processes that may generate variations in isotopic compositions: oxidation, re-equilibration of methane and other alkanes with CO2, mixing with compounds of other sources, etc. Therefore, in addition to isotopic analysis, other evidence, including water chemistry and rock compositions, are necessary to identify volatile compounds of different sources. Surprise Valley Hot Springs (SVHS, 41 deg 32'N, 120 deg 5'W), located in a typical basin and range province valley in northeastern California, is a terrestrial hydrothermal spring system of the Great Basin. Previous geophysical studies indicated the presence of clay-rich volcanic and sedimentary rocks of Tertiary age beneath the lava flows in late Tertiary and Quaternary. Water and gas samples were collected for a variety of chemical and isotope composition analyses, including in-situ pH, alkalinity, conductivity, oxidation reduction potential (ORP), major and trace elements, and C and H isotope measurements. Fluids issuing from SVHS can be classified as Na-(Cl)-SO4 type, with the major cation and anion being Na+ and SO4(2-), respectively. Thermodynamic calculation using ORP and major element data indicated that sulfate is the most dominant sulfur species, which is consistent with anion analysis results. Aquifer temperatures at depth estimated by both dissolved SiO2 and Na-K-Ca geothermometers are in the range of 125.0 to 135.4 C, and higher than the values measured at orifices (77.3 to 90.0 C). CO2 and homologs of straight chain alkanes (C1-C5) were identified in gas samples. Carbon isotope values of alkanes increase with carbon numbers. The C-13 fractionation between CO2 and dissolved inorganic carbon suggests they are out of carbon isotope equilibrium. The hypothesis regarding the formation of carbon-bearing compounds in SVHS may involve two processes: 1) Under high heat flow conditions which are caused by regional faulting and crustal extension, original high molecular weight organic compounds (kerogens) in clay-rich rocks decomposed to generate methane and other alkane homologs. 2) The SVHS area is associated with outflow structures, and distant from the heat source. Anaerobic oxidation of methane (AOM) with sulfate at shallow depth (< 90 C) is suggested as being responsible for the generation of CO2 in SVHS.

Fu, Qi↗
DOE OSTI2020

Development of a Conditioned System-Level Groundwater Model to Evaluate Long-Term Groundwater Impacts within a Performance Assessment - 20115

A preliminary performance assessment (PA) of single-shell tank Waste Management Area (WMA) A-AX located at the U.S. Department of Energy (DOE)'s Hanford Site in southeastern Washington is being conducted to satisfy the requirements of DOE Order 435.1 [1] as it relates to closure of the single-shell radioactive waste tanks in the WMA A-AX tank farms. A PA assesses the fate, transport, and impacts of radionuclides within a low-level radioactive waste disposal facility in its assumed closure configuration and the subsequent potential doses to humans over a 1,000-year compliance period and a 10,000-year performance evaluation period. The WMA A-AX preliminary PA evaluation is structured around the complementary use of process-level and system-level models to calculate the facility performance against established DOE Order 435.1 [1] performance objectives. Process-level models are those that represent a detailed phenomenological representation of processes of concern in the PA. Process models typically only represent one or a few of the components of the PA, such as groundwater flow and transport, and must be integrated with other modeling elements to perform PA calculations. System-level models are those that are abstracted from the process models, retaining the essential features of the process model, while allowing integration of all aspects of the PA in a single modeling framework. System-level models are often characterized by coarser numerical discretization, lower dimensionality, or other similar simplifications compared to the process-level model. Traditionally, a three-dimensional (3-D) process-level model is utilized primarily to evaluate the long-term impact on groundwater and the potential doses to individuals who consume contaminated groundwater. System-level models are also utilized to evaluate the groundwater pathway in PAs. These models typically have reduced dimensionality (1-D) and are conditioned utilizing flow fields (Darcy fluxes) and moisture content distributions that are abstracted from the process level models. The abstraction approach assures that the flow field in both models is consistent for a specific set of input parameters for flow, differing only in the discretization and dimensionality of the two models. The preliminary WMA A-AX PA incorporates a detailed representation of the geological system and hydraulic properties within the 3-D model STOMP{sup C} numerical code so that the effects of relevant features and processes on water flow and radionuclide transport in the subsurface can be evaluated. The complementary system-level model is developed utilizing the GoldSim{sup C} code to implement a simplified, 1-D equivalent model to represent the groundwater pathway. Contaminant transport through the vadose zone and unconfined aquifer for Tc-99 and I-129 were evaluated in each of the models. Adjustments to the saturated portion of the GoldSim{sup C} 1-D model were required to mimic the dispersion effect captured with the 3-D STOMP{sup C} model. Once this conditioning was conducted, highly similar results for the transport of Tc-99 and I-129 were achieved at the point of calculation, located 100 m downgradient from the WMA A-AX fenceline. These radionuclides represent elements that are regarded as primary dose drivers in the PA groundwater pathway analysis. The high degree of conformance between the two models suggests that the use of the equivalent 1-D system model is suitable for evaluating the full suite of radionuclides that will be released from the tank sources within WMA A-AX over the 1,000-year compliance period and over the 10,000-year evaluation period. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗
DOE OSTI2022

Machine Learning-Assisted High-Temperature Reservoir Thermal Energy Storage Optimization: Numerical Modeling and Machine Learning Input and Output Files

This data set includes the numerical modeling input files and output files used to synthesize data, and the reduced-order machine learning models trained from the synthesized data for reservoir thermal energy storage site identification. In this study, a machine-learning-assisted computational framework is presented to identify High-Temperature Reservoir Thermal Energy Storage (HT-RTES) site with optimal performance metrics by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems. All the simulations and the neural network model were done by Idaho National Laboratory. A detailed description of the work was reported in publication linked below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Evaluation of the Change in Uranium Mobility in Sediments from the Hanford 300-FF-5 Stage B Polyphosphate Field Injection

The purpose of this study was to evaluate the change in uranium mobility in the periodically re-wetted zone (PRZ) and upper unconfined aquifer of the Hanford Site’s 300-FF-5 operable unit as a result of Stage B polyphosphate (PO 4 ) injections. These 300 Area sediments are beneath or near the former North Process Pond, which has previously been shown to contain metatorbernite [Cu(UO 2 ) 2 (PO 4 ) 2 •8H 2 O] and uranophane [Ca(UO 2 ) 2 (SiO 3 OH) 2 •5H 2 O], along with aqueous and adsorbed uranium (U), and small amounts of U in calcite as a result of acidic disposal with a significant number of co-contaminants. Previous laboratory-scale studies have identified multiple mechanisms that decrease uranium mobility with polyphosphate injection: (a) precipitation of a low-solubility uranium-phosphate mineral (autunite); (b) precipitation of a Ca-phosphate solid that incorporates U; (c) precipitation of Ca-phosphate mineral (apatite group minerals) that coats existing uranium phases (i.e., U in calcite, U in Fe-oxides), reducing uranium leaching; and (d) aqueous Ca-U-carbonate species adsorption onto apatite. Given the significant amount of apatite precipitation, other U phases such as carbonates could also coprecipitate with apatite. To quantify differences in solubility, leachability, and phases of uranium between pre- and post-PO 4 -treated sediments, different types of analyses were conducted, including analysis of uranium in (a) a carbonate solution extraction from sediment over 1000 h (41.67 days), (b) six sequential liquid extractions from sediment over approximately 1 week, (c) leaching of uranium from sediments in one dimensional (1-D) column experiments over 2 months, and (d) solid phase characterization of select sediments. An additional extraction and solid phase measurements were conducted to evaluate phosphate precipitate formation. Metals were also analyzed in sequential liquid extractions as indicators of (a) amendment injection delivery, (b) co-contaminant movement, and (c) changes in solid phases. Aqueous and solid phase carbonate measurements were also conducted to evaluate changes in carbonates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Technical note: Using long short-term memory models to fill data gaps in hydrological monitoring networks

Abstract. Quantifying the spatiotemporal dynamics in subsurface hydrological flows over a long time window usually employs a network of monitoring wells. However, such observations are often spatially sparse with potential temporal gaps due to poor quality or instrument failure. In this study, we explore the ability of recurrent neural networks to fill gaps in a spatially distributed time-series dataset. We use a well network that monitors the dynamic and heterogeneous hydrologic exchanges between the Columbia River and its adjacent groundwater aquifer at the U.S. Department of Energy's Hanford site. This 10-year-long dataset contains hourly temperature, specific conductance, and groundwater table elevation measurements from 42 wells with gaps of various lengths. We employ a long short-term memory (LSTM) model to capture the temporal variations in the observed system behaviors needed for gap filling. The performance of the LSTM-based gap-filling method was evaluated against a traditional autoregressive integrated moving average (ARIMA) method in terms of error statistics and accuracy in capturing the temporal patterns of river corridor wells with various dynamics signatures. Our study demonstrates that the ARIMA models yield better average error statistics, although they tend to have larger errors during time windows with abrupt changes or high-frequency (daily and subdaily) variations. The LSTM-based models excel in capturing both high-frequency and low-frequency (monthly and seasonal) dynamics. However, the inclusion of high-frequency fluctuations may also lead to overly dynamic predictions in time windows that lack such fluctuations. The LSTM can take advantage of the spatial information from neighboring wells to improve the gap-filling accuracy, especially for long gaps in system states that vary at subdaily scales. While LSTM models require substantial training data and have limited extrapolation power beyond the conditions represented in the training data, they afford great flexibility to account for the spatial correlations, temporal correlations, and nonlinearity in data without a priori assumptions. Thus, LSTMs provide effective alternatives to fill in data gaps in spatially distributed time-series observations characterized by multiple dominant frequencies of variability, which are essential for advancing our understanding of dynamic complex systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Contributions of GRACE to Climate Monitoring

The NASA/German Gravity Recovery and Climate Experiment (GRACE) was launched in March 2002. Rather than looking downward, GRACE continuously monitors the locations of and precise distance between twin satellites which orbit in tandem about 200 km apart. Variations in mass near Earth's surface cause heterogeneities in its gravity field, which in turn affect the orbits of satellites. Thus scientists can use GRACE data to map Earth's gravity field with enough accuracy to discern month to month changes caused by ocean circulation and redistribution of water stored on and in the land. Other gravitational influences, such as atmospheric circulation, post-glacial rebound, and solid earth movements are either independently determined and removed or are negligible on a monthly to sub-decadal timescale. Despite its coarse spatial (>150,000 sq km at mid-latitudes) and temporal (approx monthly) resolutions, GRACE has enabled significant advancements in the oceanic, hydrologic, and cryospheric science, and has great potential for climate monitoring, because it is the only global observing system able to measure ocean bottom pressures, total terrestrial water storage, and ice mass changes. The best known GRACE results are estimates of Greenland and Antarctic ice sheet loss rates. Previously, scientists had estimated ice mass losses using ground and satellite based altimetry and surface mass balance estimates based on snowfall accumulation and glacier discharge. While such measurements are still very useful for their spatial detail, they are imperfectly correlated with large-scale ice mass changes, due to snow and ice compaction and incomplete spatial coverage. GRACE enables scientists to generate monthly time series of Greenland and Antarctic ice mass, which have confirmed the shrinking of the polar ice sheets, one of the most obvious and indisputable manifestations of climate change. Further, GRACE has located and quantified hot spots of ice loss in southeastern Greenland and western Antarctica. For 2002 to present, the rate of ice mass loss has been 200 to 300 GT/yr in Greenland and 70 to 210 GT/yr in Antarctica, and some scientists are suggesting that the rates are accelerating. Similarly, GRACE has been used to monitor mass changes in alpine glaciers. Tamisiea et al. first characterized glacier melt along the southern coast of Alaska, more recently estimated to be occurring at a rate of 84 GT/yr. Chen et al. estimated that Patagonian glaciers are melting at a rate of 28 GT/yr, and estimated that the high mountains of central Asia lose ice at a rate of 47 GT/yr. Tapley et al. and Wahr et al. presented the first GRACE based estimates of changes in column-integrated terrestrial water storage (TWS; the sum of ground-water, soil moisture, surface waters, snow, ice, and water stored in vegetation) at continental scales. Since then, dozens of studies have shown that GRACE based estimates of regional to continental scale TWS variations agree with independent information, and some innovative uses of GRACE data have been developed. Rodell et al. (2004) and Swenson and Wahr (2006) demonstrated that by combining GRACE derived terrestrial water storage changes with observations of precipitation and runoff in a river basin scale water budget, it was possible to produce new estimates of evapotranspiration and atmospheric moisture convergence, essential climate variables that are difficult to estimate accurately. Similarly, GRACE has been used to constrain estimates of global river discharge and the contribution of changes in TWS to sea level rise. Crowley et al. observed a negative correlation between interannual TWS anomalies in the Amazon and the Congo River basin. Yeh et al. and Rodell et al. estimated regionally averaged groundwater storage variations based on GRACE and auxiliary observations. Rodell et al. and Tiwari et al. applied that method to quantify massive groundwater depletion in northern India caused by over reliance on aquifers for irration, and Famiglietti et al. found a similar situation in California's Central Valley. Zaitchik et al. and Lo et al. described approaches to use GRACE to constrain hydrological models, enabling integration of GRACE data with other observations and achieving much higher spatial and temporal resolutions than GRACE alone. Such approaches are now supporting applications including drought and water resources monitoring. Oceanography has likewise benefitted from the independent nature of GRACE observations. One application is measurement of the mass component of sea level rise, which complements radar altimetry and in situ measurements. GRACE also measures ocean bottom pressures (OBP), which help to refine understanding and modeling of ocean circulation and the ocean's fresh water budget, among other things. For example, Hayakawa et al. showed that GRACE observes OBP patterns absent from the background models of oceanic variability. Morison et al. used GRACE to describe important decadal scale shifts in circulation and an ongoing trend of freshening of the western Arctic, important indicators of climate variability. The research of Song and Zlotnicki and Chambers and Willis on GRACE-derived ocean bottom pressures in the sub-polar gyre led to the discovery of an ENSO teleconnection and a long-term change in OBP in the North Pacific sub-polar gyre that was not predicted by an ocean model. Further, Chambers and Willis were able to identify an internal redistribution of mass between Atlantic and Pacific Oceans lasting at least six years, which was not predicted by ocean models and was the first direct evidence of sustained mass transport from one ocean basin to another on periods longer than a year. Boening et al. observed a record increase in OBP over part of the southeastern Pacific in late 2009 and early 2010, primarily caused by wind stress curl associated with a strong and persistent anticyclone and likely related to the concurrent Central Pacific El Nino. GRACE has far surpassed its 5-year design lifetime, but it will likely succumb to the aging of batteries and instrument systems sometime in the next few years. NASA has begun initial development of a follow-on to GRACE with very similar design, which could launch as soon as 2016 and would provide continuity in the data record while improving resolution slightly. Higher resolution time variable gravity missions are also on the drawing board .

Rodell, Matthew↗
DOE OSTI2021

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗
DOE OSTI2020

Pahute Mesa-Oasis Valley Hydrostratigraphic Framework Model for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada (Rev. 0)

A new, revised three-dimensional (3-D) hydrostratigraphic framework model (HFM) for the Central and Western Pahute Mesa Corrective Action Units 101 and 102 was completed in 2019. The initial (Phase I) model was completed in 2002 and was documented in A Hydrostratigraphic Model and Alternatives for the Groundwater Flow and Contaminant Transport Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada (BN, 2002). Subsequent flow and transport modeling revealed the need to collect additional data to reduce uncertainties. The Phase II data collection and characterization effort included extensive well drilling and testing activities, and rebuilding of the Phase I 3-D HFM, resulting in the Phase II HFM (NSTec, 2014). Further drilling (Well ER-20-12), data collection, and interpretive work (e.g., Thirsty Canyon Lineament [Wurtz and Day, 2018]) were incorporated into the Pahute Mesa-Oasis Valley (PM-OV) HFM. The PM-OV HFM incorporates the area of interest from the Phase II HFM and adds additional areas, primarily to the north and west with minor additions to the south and east. The new, expanded area matches the PM-OV groundwater basin as defined in Fenelon et al. (2016) and includes Pahute Mesa, a former nuclear testing area at the Nevada National Security Site; and Oasis Valley, a groundwater discharge area downgradient from contaminant source areas on Pahute Mesa. The data and information necessary to build the PM-OV HFM are provided in Navarro (2019). The PM-OV HFM will be used by hydrologic modelers who are tasked with developing a model to determine how contaminants are transported by groundwater flow in an area of complex geology. The model area is large (more than 6,250 square kilometers) and geologically complex, including Paleozoic- to Mesozoic-age structures (e.g., thrust faults, Thirsty Canyon Lineament), at least seven Tertiary-age calderas, several intrusive bodies, and many relatively recent basin-and-range normal faults. Investigators from the U.S. Geological Survey; National Security Technologies, LLC; National Laboratories; and U.S. Department of Energy contractors have organized the more than 300 volcanic units and other rocks in the model area into 77 hydrostratigraphic units (HSUs). The volcanic rocks are subdivided into 42 aquifers, 11 confining units, and 9 composite units (containing both aquifer and confining unit rocks). The underlying pre-Tertiary rocks were divided into 5 HSUs, comprising 2 aquifers and 3 confining units. The model also includes 1 alluvial aquifer unit, 8 intrusive confining units, and 1 granite confining unit. The model depicts the thickness, extent, and geometric relationships of these HSUs (“layers” in the model) along with all the major structural features that control them, including calderas and faults. The most substantive differences between the Phase II HFM and the PM-OVHFM include (1) addition of the new data and resulting revised interpretations from the additional area encompassed by the PM-OV groundwater basin; (2) incorporation of 16 new faults; (3) 12 boreholes from the northern extension; (4) 2 additional caldera-collapse collars; and (5) correlation of the simulated faulting and hydrostratigraphy in the eastern part of the model domain with the delineations represented in the neighboring Rainier Mesa-Shoshone Mountain HFM-Prime and Yucca Flat-Climax Mine HFM. All PM-OV HFM delineations were based on published and unpublished data from geophysical investigations, surface geologic maps, and drill-hole records. Phase II approached model uncertainty using the model-building software to test model configurations and hydrologic consequences using a range of values for parameters such as fault offset amounts and HSU elevations. The PM-OV HFM expands on this approach to encompass the entire PM-OV groundwater basin.

54 ENVIRONMENTAL SCIENCES↗