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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Biogenic coal-to-methane conversion can be enhanced with small additions of algal amendment in field-relevant upflow column reactors

Subsurface coal environments, where biogenic coal-to-methane conversion occurs, are difficult to access, resulting in inherent challenges and expenses for in situ experiments. Previous batch reactor studies provided insights into specific processes, pathways, kinetics, and engineering strategies, but field-relevance is restricted due to limited substrate availability or byproduct accumulation that may influence reactions or metabolisms. In this study, continuous-flow column reactors were used to overcome some batch limitations, improve the understanding of in situ conditions, and increase field-relevance for subsurface engineering technology development. The bench-scale reactor system was constructed to investigate the addition of algal amendment for enhancing microbial coal-to-methane conversion previously developed in batch systems. Four reactor columns were packed with coal and inoculated with a microbial consortium from the same Flowers-Goodale coal bed. Two reactors were amended with 13 C-labeled algal amendment on day 0, and two were unamended. On day 61, one previously amended and one previously unamended reactor were re-amended. Produced gases were captured in a gas trap, and CH 4 and CO 2 were quantified. The reactor amended twice produced 1712.6 µmol CH 4 (4.6% as 13 CH 4 ). The reactor amended only on day 0 produced 1485.5 µmol CH 4 (2.6% as 13 CH 4 ). The reactor amended only on day 61 produced 278.9 µmol CH 4 (3.9% as 13 CH 4 ). The reactor with no amendment produced no measurable gases for the duration of the 172-day experiment. Amendment increased the rate of coal-to-methane conversion and total gas production; most of the produced gases were due to coal conversion with only small contributions (<7%) from amendment conversion.

09 BIOMASS FUELS↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Aram Chaos: a Long Lived Subsurface Aqueous Environment with Strong Water Resources Potential for Human Missions on Mars

Aram Chaos, Mars is a crater 280 kilometers in diameter with elevations circa. minus 2 to minus 3 kilometers below datum that provides a compelling landing site for future human explorers as it features multiple scientific regions of interest (ROI) paired with a rich extensible Resource ROI that features poly-hydrated sulfates [1]. The geologic history of Aram Chaos suggests several past episodes of groundwater recharge and infilling by liquid water, ice, and other materials [1-3]. The creation of the fractured region with no known terrestrial equivalent may have been caused by melting of deep ice reservoirs that triggered the collapse of terrain followed by catastrophic water outflows over the region. Aram Chaos is of particular scientific interest because it is hypothesized that the chaotic terrain may be the source of water that contributed to the creation of nearby valleys such as Ares Vallis flowing toward Chryse Planitia. The liquid water was likely sourced as groundwater and therefore represents water derived from a protected subsurface environment making it a compelling astrobiological site [2]. The past history of water is also represented by high concentrations of hematite, Fe-oxyhydroxides, mono-hydrated and poly-hydrated sulfates [1, 2]. Poly-hydrated sulfates are likely to contain abundant water that evolves at temperatures below 500 degrees Centigrade thus conferring Aram Chaos a potentially high value for early in-situ resource utilization (ISRU) [4]. The geologic history also calls for future prospecting of deep ice deposits and possibly liquid water via deep drilling. The most recent stratigraphic units in the central part of Aram Chaos are not fractured, and are part of a dome-shaped formation that features bright, poorly-consolidated material that contains both hydrated sulfates and ferric oxides according to OMEGA (Observatoire pour la Minéralogie, l'Eau, les Glaces et l'Activité) data [5]. These surface material characteristics are preliminary indications of their potential use in civil engineering activities that involve regolith moving and hauling, while further study is needed to assess traverse-ability challenges. The widespread distribution of sulfates is also of interest as a resource for the use of sulfur as a binding compound in regolith-based concrete for constructions. The terrain depressions caused by the rock fracturing events may challenge surface mobility but also suggest the possibility of using such natural features for additional shielding from space radiation and as emplacement of nuclear surface power reactors for the same reason. The high concentration of hematite (up to 16 percent) in some of the smoother recent terrains of the central part of Aram Chaos [2] is a favorable attribute for metal extraction ISRU to create iron-based feedstock for in-situ fabrication of replacement parts or their repairs. Preliminary data on Aram Chaos indicate that it offers a combination of many critical criteria for human missions to the surface of Mars: equatorial region at low Mars Orbiter Laser Altimeter (MOLA), evidence of hydrated minerals over large areas and at high concentrations tied to historic evidence of liquid water over long periods.

Sibille, L.↗
DOE OSTI2019

Quantity and quality of water percolating below the root zone of three biofuel feedstock crop systems

Biofuel demand is on the rise, yet there is limited information available about the likely implications of increased planting of cellulosic biofuel crops on water and nutrient fluxes. To help address this gap, this study quantifies subsurface water, nitrate, and dissolved phosphorus fluxes below feedstock crops managed for biofuel production on naturally well drained soils in southern Wisconsin. Automated equilibrium tension lysimeters were used to measure subsurface (below root zone) fluxes of water and nutrients below continuous corn (CC), monoculture switchgrass (SG), and hybrid poplar (HP) cropping systems for three years. Average water fluxes below SG were nearly double those below CC and HP. Most of this difference can be linked to seasonal events, particularly during late winter and early spring when subsurface water fluxes and increased soil moisture were observed below SG plots, likely due to enhanced macropore flow. Nitrate fluxes were approximately 14 times higher below CC than the 37.4 kg NO 3 -N/ha three year annual average observed for SG; the nitrate fluxes were negligible below HP despite similar nitrogen application rates. The high flux in one of the two CC plots is likely due to excessive manure applications prior to the study period; however, the nitrate loss potential from CC is much higher than the perennial crops. Dissolved reactive phosphorus (DRP) concentrations were very low below all three cropping systems (average of 0.076, 0.083, 0.034 mg/L in SG, CC, and HP respectively). For the cropping systems evaluated, nitrate fluxes were driven by concentration differences in leachate and not influenced by nominal differences in N fertilizer inputs, whereas DRP fluxes were mainly driven by water fluxes and therefore independent of phosphorus inputs and soil test phosphorus. The data in this paper provide a foundation to better understand differences in fluxes below different biofuel cropping systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Can We Probe the Conductivity of the Lithosphere and Upper Mantle Using Satellite Tidal Magnetic Signals?

A few studies convincingly demonstrated that the magnetic fields induced by the lunar semidiurnal (M2) ocean flow can be identified in satellite observations. This result encourages using M2 satellite magnetic data to constrain subsurface electrical conductivity in oceanic regions. Traditional satellite-based induction studies using signals of magnetospheric origin are mostly sensitive to conducting structures because of the inductive coupling between primary and induced sources. In contrast, galvanic coupling from the oceanic tidal signal allows for studying less conductive, shallower structures. We perform global 3-D electromagnetic numerical simulations to investigate the sensitivity of M2 signals to conductivity distributions at different depths. The results of our sensitivity analysis suggest it will be promising to use M2 oceanic signals detected at satellite altitude for probing lithospheric and upper mantle conductivity. Our simulations also suggest that M2 seafloor electric and magnetic field data may provide complementary details to better constrain lithospheric conductivity.

Schnepf, N. R.↗
NASA NTRS2017

Electric Potential and Electric Field Imaging with Dynamic Applications: 2017 Research Award Innovation

The technology and methods for remote quantitative imaging of electrostatic potentials and electrostatic fields in and around objects and in free space is presented. Electric field imaging (EFI) technology may be applied to characterize intrinsic or existing electric potentials and electric fields, or an externally generated electrostatic field may be used for illuminating volumes to be inspected with EFI. The baseline sensor technology (e-Sensor) and its construction, optional electric field generation (quasi-static generator), and current e- Sensor enhancements (ephemeral e-Sensor) are discussed. Critical design elements of current linear and real-time two-dimensional (2D) measurement systems are highlighted, and the development of a three dimensional (3D) EFI system is presented. Demonstrations for structural, electronic, human, and memory applications are shown. Recent work demonstrates that phonons may be used to create and annihilate electric dipoles within structures. Phonon induced dipoles are ephemeral and their polarization, strength, and location may be quantitatively characterized by EFI providing a new subsurface Phonon-EFI imaging technology. Initial results from real-time imaging of combustion and ion flow, and their measurement complications, will be discussed. These new EFI capabilities are demonstrated to characterize electric charge distribution creating a new field of study embracing areas of interest including electrostatic discharge (ESD) mitigation, crime scene forensics, design and materials selection for advanced sensors, combustion science, on-orbit space potential, container inspection, remote characterization of electronic circuits and level of activation, dielectric morphology of structures, tether integrity, organic molecular memory, atmospheric science, and medical diagnostic and treatment efficacy applications such as cardiac polarization wave propagation and electromyography imaging.

Generazio, Ed↗
DOE OSTI2025

Transportability of exogenous microbial community correlates with interwell connectivity in deep aquifers

Subsurface resource engineering operations often utilize continuous injection of externally-sourced water into geological reservoirs for formation pressure maintenance, resource recovery or energy/waste storage. Such injected water generally contains naturally occurring microbes. Little is known, however, about how the injectate microbes transport through geological media as a community, how such transportability is affected by injector-producer connectivity, and whether such knowledge can be utilized for flowpath characterization. In this study, we analyzed daily-to-weekly timeseries microbial community data from the injected- and produced-fluids of a ten-month flow test at a deep, well-characterized engineered aquifer. We found that the injectate microbial community was distinct from the indigenous community at the amplicon sequence variant (ASV) level, and that the transportability of injectate community towards a given producer, quantified by an “nASV-Overlap” metric we propose, had strong and significant positive correlation with known injector-producer connectivities at our site. This suggests that the better the connectivity, the higher the probability for more injectate species to flow through the interwell region and arrive at a producer. Because interwell connectivity is an important yet usually unknown parameter in subsurface resource engineering, such correlation in turn points to nASV-Overlap as a useful indicator of interwell connectivity for aquifer characterization and long-term monitoring. Based on our findings, an nASV-Overlap-based microbial tracing approach was developed for characterizing and monitoring the relative connectivities across multiple producers with a given injector. A side-by-side comparison between the new nASV-Overlap approach and traditional artificial tracer methods is presented, and their respective strengths and limitations are discussed.

Deep biosphere↗
DOE OSTI2020

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗
DOE OSTI2022

Geological activity shapes the microbiome in deep-subsurface aquifers by advection

Subsurface environments host diverse microorganisms in fluid-filled fractures; however, little is known about how geological and hydrological processes shape the subterranean biosphere. Here, we sampled three flowing boreholes weekly for 10 mo in a 1478-m-deep fractured rock aquifer to study the role of fracture activity (defined as seismically or aseismically induced fracture aperture change) and advection on fluid-associated microbial community composition. We found that despite a largely stable deep-subsurface fluid microbiome, drastic community-level shifts occurred after events signifying physical changes in the permeable fracture network. The community-level shifts include the emergence of microbial families from undetected to over 50% relative abundance, as well as the replacement of the community in one borehole by the earlier community from a different borehole. Null-model analysis indicates that the observed spatial and temporal community turnover was primarily driven by stochastic processes (as opposed to deterministic processes). We, therefore, conclude that the observed community-level shifts resulted from the physical transport of distinct microbial communities from other fracture(s) that outpaced environmental selection. Given that geological activity is a major cause of fracture activity and that geological activity is ubiquitous across space and time on Earth, our findings suggest that advection induced by geological activity is a general mechanism shaping the microbial biogeography and diversity in deep-subsurface habitats across the globe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
NASA NTRS1993

Sunspot dynamics

This report describes recent results of our theoretical and observational work on dynamical phenomena in sunspots. The overall goal of this research has been a better understanding of the various oscillatory, transient, and steady motions in a sunspot and their relation to the basic structure of the sunspot. The principal topics of the research reported here are the following: (1) sunspot seismology, i.e., the study of the interaction of solar p-modes with a sunspot as a probe of the subsurface structure of a sunspot; (2) local sources of acoustic waves in the solar photosphere; and (3) siphon flows in isolated magnetic flux tubes and their relation to the photospheric Evershed flow and to intense magnetic elements outside of sunspots.

Thomas, John H.↗
NASA NTRS1992

Thermographic imaging of cracks in thin metal sheets

The presence of cracks significantly decreases the structural integrity of thin metal sheets used in aerospace applications. Thermographic detection of surface temperature variations due to these cracks is possible after external heating. An approximate line source of heat is used to produce an inplane flow of heat in the sheet. A crack in the sheet perturbs the inplane flow of heat and can be seen in an image of the surface temperature of the sheet. An effective technique for locating these perturbations is presented which reduces the surface temperature image to an image of variations in the inplane heat flow. This technique is shown to greatly increase the detectability of the cracks. This thermographic method has advantages over other techniques in that it is able to remotely inspect a large area in a short period of time. The effectiveness of this technique depends on the shape, position and orientation of the heat source with respect to the cracks as well as the extent to which the crack perturbs the surface heat flow. The relationship between these parameters and the variation in the heat flow is determined both by experimental and computational techniques. Experimental data is presented for through-the-thickness, subsurface and surface EDM notches. Data for through-the-thickness fatigue cracks are also presented.

Cramer, K. E.↗
NASA NTRS2001

The Evolution of the Martian Hydrosphere: Implications for the Fate of a Primordial Ocean and the Current State of the Northern Plains

In this paper we consider the hydraulic and thermal conditions that gave rise to the elevated source regions of the Late Hesperian outflow channels and explore their implications for the evolution of the Martian hydrosphere. We find that if the outflow channel floodwaters were derived from a subpermafrost aquifer, then it implies that, throughout the planet's first billion years of evolution, as much as one third of its surface was covered by standing bodies of water and ice. Following the development of the global dichotomy, the bulk of this water would have existed as an ice-covered ocean in the northern plains. We demonstrate that the progressive crustal assimilation of this early surface reservoir of H2O (punctuated by possible episodes of less extensive flooding) was a natural consequence of the planet's subsequent climatic and geothermal evolution-potentially cycling the equivalent of a km-deep global ocean of water through the atmosphere and subsurface every approx. 10(exp 9) years. In response to the long-term decline in planetary heat flow, the progressive coldtrapping of H2O into the growing cryosphere is expected to have significantly depleted the original inventory of groundwater-a development that could well explain the apparent decline in outflow channel activity observed during the Amazonian. Although primarily a theoretical analysis, our findings appear remarkably consistent with the geomorphic and topographic evidence that Mars once possessed a primordial ocean and that a substantial relic of that body continues to survive as massive ice deposits within the northern plains. Confirmation of the presence of such deposits, combined with the potential detection of a global-scale groundwater system, would provide persuasive support for the validity of this analysis.

Clifford, Stephen M.↗
NASA NTRS2012

Development of a Deep-Penetrating, Compact Geothermal Heat Flow System for Robotic Lunar Geophysical Missions

Geothermal heat flow measurements are a high priority for the future lunar geophysical network missions recommended by the latest Decadal Survey of the National Academy. Geothermal heat flow is obtained as a product of two separate measurements of geothermal gradient and thermal conductivity of the regolith/soil interval penetrated by the instrument. The Apollo 15 and 17 astronauts deployed their heat flow probes down to 1.4-m and 2.3-m depths, respectively, using a rotary-percussive drill. However, recent studies show that the heat flow instrument for a lunar mission should be capable of excavating a ~3-m deep hole to avoid the effect of potential long-term changes of the surface thermal environment. For a future robotic geophysical mission, a system that utilizes a rotary/percussive drill would far exceed the limited payload and power capacities of the lander/rover. Therefore, we are currently developing a more compact heat flow system that is capable of 3-m penetration. Because the grains of lunar regolith are cohesive and densely packed, the previously proposed lightweight, internal hammering systems (the so-called moles ) are not likely to achieve the desired deep penetration. The excavation system for our new heat flow instrumentation utilizes a stem which winds out of a pneumatically driven reel and pushes its conical tip into the regolith. Simultaneously, gas jets, emitted from the cone tip, loosen and blow away the soil. Lab tests have demonstrated that this proboscis system has much greater excavation capability than a mole-based heat flow system, while it weighs about the same. Thermal sensors are attached along the stem and at the tip of the penetrating cone. Thermal conductivity is measured at the cone tip with a short (1- to 1.5-cm long) needle sensor containing a resistance temperature detector (RTD) and a heater wire. When it is inserted into the soil, the heater is activated. Thermal conductivity of the soil is obtained from the rate of temperature increase during the heating. By stopping during the excavation, it is possible to measure thermal conductivities at different depths. The gas jets are turned off when the penetrating cone reaches the target depth. Then, the stem pushes the needle sensor into the undisturbed soil at the bottom of the hole and carries out a thermal conductivity measurement. When the measurement is complete, the system resumes excavation. RTDs, placed along the stem at short (approx 30 cm) intervals, will monitor long-term temperature stability of the subsurface. Temperature in the shallow subsurface would fluctuate with the diurnal, annual, and precession cycles of the Moon. These thermal waves penetrate to different depths into the regolith. Longterm monitoring of the subsurface temperature would allow us to accurately delineate these cyclic signals and separate them from the signal associated with the outward flow of the Moon s endogenic heat. Further, temperature toward bottom of the 3-m hole should be fairly stable after the heat generated during the excavation dissipates into the surrounding soil. The geothermal gradient may be determined reliably from temperature measurements at the RTDs near the bottom. In order to minimize the heat conduction along the stem from affecting the geothermal gradient measurements, we plan to use low-conductive materials for the stem and develop a mechanism to achieve close coupling between the RTDs and the wall of the excavated hole.

Nagihara, Seiichi↗
NASA NTRS1993

The thumbprint terrain: What will Mars Observer tell us?

Some of the more puzzling features seen on Mars are those known as curvilinear features or 'thumb print' terrain, types of patterned ground found in the northern plains. The thumbprint terrain, named for its resemblance to the lines of a human thumbprint is found on what appears to be level, relatively low-lying ground near the crustal dichotomy boundary. It is often found near the termini of large channels. There are three types of thumbprint terrain in the classification of Rossbacher and Judson. The first consists of ridges and depressions about a kilometer wide, separated by a few kilometers, and with an apparent relief of 40-100 m. The second consists of steep-sided, flat-floored depressions about a kilometer wide, separated by a few kilometers, and with an apparent relief of 10-290 m. The third consists only of albedo markings, with a typical scale of about a kilometer. Many models were proposed to explain the origin of this terrain. It was suggested that it was caused by lava flows, removal of debris mantles, glaciers, or karst. However, the most popular models at present involve the action of subsurface ice to form such thermokarst features as striped ground, solifluction lobes, and/or linear, ice-cored ridges. There are several instruments on the Mars Observer spacecraft that will be able to provide us with information useful in distinguishing between these models.

Schaefer, M. W.↗
DOE OSTI2019

Radiokrypton unveils dual moisture sources of a deep desert aquifer

We report that in arid regions, groundwater is a vital resource that can also provide a long-term record of the regional water cycle. However, the use of groundwater as a paleoclimate proxy has been limited by the complex hydrology and the lack of appropriate chronometers to determine the recharge time without complication. Applying 81 Kr, a long-lived radioisotope tracer, we investigate the paleohydroclimate and subsurface water storage properties of the Nubian Sandstone Aquifer in the Negev Desert, Israel. Based on the spatial distributions of stable isotopes and the abundance of 81 Kr, we resolve subsurface mixing and identify two distinct moisture sources of the recharge: one recent (<38 ky ago) from the Mediterranean and the other 361 ± 30 ky ago from the tropical Atlantic, both of which occurred under conditions of low orbital eccentricity comparable to that of the present. The recent recharge provided by the moisture from Mediterranean cyclones can be attributed to the southward shift of the storm track during the Last Glacial Maximum, and the earlier recharge can be attributed to moisture from the Atlantic delivered as tropical plumes under a climate colder than the present. Furthermore, the residence time of the latter reveals that tectonically active terrain can store groundwater for an unexpectedly long period, likely due to strongly attenuated groundwater flow across the fault zones. Lastly with this tracer, groundwater can now serve as a direct record of paleoprecipitation over land and of subsurface water storage from the mid-Pleistocene and onward.

58 GEOSCIENCES↗
DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗
NASA NTRS1991

Strain accommodation beneath structures on Mars

A recent review of tectonic features on Mars shows that most of their subsurface structures can be confidently extended only a few kilometers deep (exceptions are rifts, in which bounding normal faults penetrate the entire brittle lithosphere, with ductile flow at deeper levels). Nevertheless, a variety of estimates of elastic lithosphere thickness and application of accepted failure criteria under likely conditions on Mars suggest a brittle lithosphere that is many tens of kilometers thick. This raises the question of how the strain (extension or shortening) accommodated by grabens and wrinkle ridges within the upper few kilometers is being accommodated at deeper levels in the lithosphere. Herein, the nonrift tectonic features present on Mars are briefly reviewed, along with their likely subsurface structures, and some inferences and implications are presented for behavior of the deeper lithosphere.

Golombek, M. P.↗