Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 739 records · Page 41

Continental-Scale Geographic Trends in Barometric-Pumping Efficiency Potential: A North American Case Study

Barometric pumping is a gas transport mechanism that has important implications for many applications involving subsurface gas seepage processes. Here, this study provides the first continental-scale analysis of barometric-pumping efficiency potential based on meteorology. We quantified the barometric-pumping efficiency potential at 1,257 locations across the continental US and Canada. The results provide continental-scale geographic dependencies of barometric-pumping efficiency potential, indicating a significant correlation with latitude and a nonlinear dependence on longitude. The analysis also indicates that variability in barometric-pumping efficiency potential decreases with distance from the coast and as elevation increases. Locations far from the coastline are more likely to have upper mid-range potentials, while higher elevation locations are more likely to have low potentials. The highest barometric-pumping efficiency potentials are mostly found around the Gulf of St. Lawrence around 50°N. Locations along the Atlantic coast exhibit large-scale variations in potentials with a clear increasing trend with latitude.

58 GEOSCIENCES↗

Distribution of Infiltration in the 216-U-10 and 216-B-3 Pond Systems 1944-1997

A modeling system to simulate fate and transport of contaminant plumes in groundwater beneath the Central Plateau to support the Cumulative Impact Evaluation (CIE) and Composite Analysis (CA) is being developed. The CIE will evaluate the cumulative effect of multiple Central Plateau radiological and chemical sources on groundwater quality (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document ). This will enable cleanup decisions to be evaluated in the context of all Central Plateau waste sites and existing groundwater conditions. The CA evaluates the cumulative effect of all sources of radioactive materials on radiological dose to future members of the public (DOE/RL-2018-60, Annual Status Report (FY 2018): Composite Analysis of Low Level Waste Disposal in the Central Plateau of the Hanford Site ). Fate and transport modeling will be used to forecast future groundwater concentrations for both the CIE and CA. The CIE/CA subsurface modeling system is divided into source (inventory), vadose zone, and saturated zone facets. Source inventory is a necessary input to the vadose zone facet, and releases from the vadose zone to groundwater are inputs to the saturated zone facet. The vadose zone facet simulates historical and future releases from the vadose zone to groundwater for both liquid waste and solid waste disposal sites. Effluent release volumes and radionuclide inventory estimates for liquid waste sites are contained in ECFHANFORD- 17-0079, Hanford Soil Inventory Model (SIM-v2) Calculated Radionuclide Inventory of Direct Liquid Discharges to Soil in the Hanford Site’s 200 Areas , herein referred to as SIM-v2. The liquid waste sites in SIM-v2 include the 216-B-3 Pond near 200 East Area and the 216-U-10 Pond in 200 West Area. The 216-B-3 Pond system consisted of a main pond and 3 expansion lobes that were operational during overlapping time periods. It also included influent ditches, but these are not listed as liquid waste sources in SIM-v2 and are not included in the CIE/CA modeling. The total effluent volume released to the pond system was measured, but the partitioning of this volume to the various expansion lobes was not determined and is not accounted for in SIM-v2. Regarding the 216-U-10 Pond, it received water from several influent ditches. The pond and the influent ditches are listed as liquid waste sites in SIM-v2 and are included as sources in the CIE/CA modeling. Because water in the ditches flowed into the pond, only a portion of the effluent volumes and contaminant inventories assigned to the ditches in SIM-v2 infiltrated from the ditches. Thus, partitioning of infiltration and contaminant inventories between the ditches and the pond is not fully accounted for in SIM-v2. Further, some of the water in the 216-U-10 Pond entered overflow ditches and the amount of this overflow was not recorded and is not accounted for in SIM-v2. Thus, the purpose of this Environmental Calculation File (ECF) is to partition the infiltration of wastewater and contaminant inventory between the main pond and each expansion lobe of the 216-B-3 Pond system, and to partition the infiltration of wastewater and contaminant inventory between the ditches and pond of the 216-U-10 Pond system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Lateral and Vertical Heterogeneity of Thorium in the Procellarum KREEP Terrane: As Reflected in the Ejecta Deposits of Post-Imbrium Craters

The Procellarum KREEP Terrane displays the highest concentrations of Th on the Moon. However, locations of elevated Th in this region appear to be random. As observed in the 5 deg per pixel equal-area Th data, and made more evident in the preliminary 2 deg data, Th is enhanced around the craters Aristillus, Aristarchus, Kepler, Mairan, the Apennine Bench formation, and the Fra Mauro region, while noticeably and unexpectedly lower in other locations (e.g., Archimedes, Copernicus, Eratosthenes, and Plato). We have examined the composition of the materials present in these regions with the goal of understanding the patchy nature to the distribution of Th and ultimately to decipher the geologic processes that have concentrated the Th. At present time, the published resolution of the Lunar Prospector Th gamma-ray data is low (5 deg per pixel), but this will soon be superceded by significantly higher-resolution data (2 deg per pixel). Even at this improved resolution, however, it is difficult to resolve the units that are the major source of Th. In an attempt to circumvent this problem, we employ the higher-resolution Clementine multispectral data for those regions mentioned above. We use the UV-VIS-derived compositional information and the spectral properties of craters, and their ejecta as drill holes through the mare-basalt surface to investigate the thickness and composition of underlying material. With this information we attempt to piece together the stratigraphy and geologic history of the Imbrium-Procellanim region. We processed the five-band multispectral data from the Clementine Mission (415, 750,900,950, and 1000nm) using ISIS software and calibration parameters developed by the USGS, Flagstaff, Arizona. Final image mosaics are in equal-area sinusoidal projection, and have a resolution of 250 m/pixel. Using the method of we produced maps of FeO and Ti02 composition. Here we examine the Th, FeO, and Ti02 composition and spectral properties of the craters discussed above and their ejecta, with the goal of describing the materials they excavate. One interpretation for the origin of the high-Th material is that subsurface KREEPy materials have been excavated by impact craters. The material excavated may be either volcanic KREEP (e,g., Apennine Bench Formation), KREEPy impact-melt breccia formed by the Imbrium impact (e.g., Fra Mauro Formation), or other KREEP-rich crustal material. Determining which type of material is responsible for the elevated Th and its extent is important to understanding the premare and possibly the prebasin stratigraphy of the Imbrium-Procellarum Region. Merging the 5 deg. Th data with the shaded relief map, we observe that the highest Th concentrations are not related to pre-Imbrium upper crustal materials. The Apennines, Alpes, and Caucasus Mountains represent the pre-Imbrian highlands material and do not express concentrations of Th, FeO, and TiO2 as high as the most Th-fich materials exposed within the Procellarum KREEP Terrane. We observe that, in general, these massifs contain 10-14 wt% FeO and 4-7 ppm Th. Determining whether the Th signal is from KREEP basalts or KREEPy impact-melt breccias cannot be done with the Clementine data because the two rock types are compositionally and mineralogically too similar (e.g., the Th-rich, mafic impact-melt breccias in the Apollo sample collection are dominated by a KREEP-basalt like component. Mapping-the distribution and sizes of craters and whether they display elevated Th concentrations or not, should reveal the depth and thickness of the KREEP-rich materials, and whether they are ubiquitous (i.e., impact-melt breccia) or more randomly distributed; this might be taken as an indicator of localized KREEP-basalt flows. Within the southeastern region of the Imbrium basin, there are two Th hot spots. The first is associated with the crater Aristillus, and the latter with the Apennine Bench Formation. Adjacent to these two hot spots are craters with a lower Th signature: Archimedes and Autolycus. We observe in the ejecta of Aristillus, a region of significantly lower FeO (10-14 wt%) relative to the surrounding mare basalt. The crater Autolycus, 50 km to the south, did not excavate similar low-FeO material. We suggest that the lower-FeO material in the ejecta of Aristillus corresponds to Th-rich material; the FeO content observed in Aristillus ejecta is comparable to that of KREEP basalt or mafic impact melt breccia (10-12 wt% FeO). We determine that this low FeO, Th-rich material is volcanic KREEP, as opposed to Imbrium impact melt, on the basis that the low-Fe material is exposed more prominently in ejecta in the northern portion of Aristillus. Our assumption is that if the layer underlying Aristillus was continuous, a more widespread and uniform low-Fe signature wouldbe observed in the ejecta deposit. Archimedes, 110 km southwest of Aristillus, impacted the northern portion of the Apennine Bench prior to the eruption of KREEP basalt. Archimedes rim material is not as enriched in Th as the Apennine Bench, and there are differences between the two in FeO concentration and in their continuum slope. Archimedes exhibits a much steeperor "redder" continuum slope than the Apennine Bench. This steepslope suggests the presence of glassy material. The glassy material is concentrated around an unnamed crater on the southern rim of Archimedes (4.5W, 28.2N) and along the northern rim of Archimedes. We suggest two possibilities, or a combination of the two, to explain the low-Th signal from Archimedes: (1) The Apennine Bench prior to KREEP basalt eruption was lower in Th (4-7 ppm, e.g., similar to the Apennine massifs) and KREEP basalts are absent in the rim of Archimedes; or (2) the glassy (possibly pyroclastic) material layering the rim of the Archimedes, dilutes any high-Th material present with low-Th material. (Additional information is contained in the original)

Gillis, J. J.↗

Mud Volcanoes - Analogs to Martian Cones and Domes (by the Thousands!)

Mud volcanoes are mounds formed by low temperature slurries of gas, liquid, sediments and rock that erupt to the surface from depths of meters to kilometers. They are common on Earth, with estimates of thousands onshore and tens of thousands offshore. Mud volcanoes occur in basins with rapidly-deposited accumulations of fine-grained sediments. Such settings are ideal for concentration and preservation of organic materials, and mud volcanoes typically occur in sedimentary basins that are rich in organic biosignatures. Domes and cones, cited as possible mud volcanoes by previous authors, are common on the northern plains of Mars. Our analysis of selected regions in southern Acidalia Planitia has revealed over 18,000 such features, and we estimate that more than 40,000 occur across the area. These domes and cones strongly resemble terrestrial mud volcanoes in size, shape, morphology, associated flow structures and geologic setting. Geologic and mineralogic arguments rule out alternative formation mechanisms involving lava, ice and impacts. We are studying terrestrial mud volcanoes from onshore and submarine locations. The largest concentration of onshore features is in Azerbaijan, near the western edge of the Caspian Sea. These features are typically hundreds of meters to several kilometers in diameter, and tens to hundreds of meters in height. Satellite images show spatial densities of 20 to 40 eruptive centers per 1000 square km. Many of the features remain active, and fresh mud flows as long as several kilometers are common. A large field of submarine mud volcanoes is located in the Gulf of Cadiz, off the Atlantic coasts of Morocco and Spain. High-resolution sonar bathymetry reveals numerous km-scale mud volcanoes, hundreds of meters in height. Seismic profiles demonstrate that the mud erupts from depths of several hundred meters. These submarine mud volcanoes are the closest morphologic analogs yet found to the features in Acidalia Planitia. We are also conducting laboratory analyses of surface samples collected from mud volcanoes in Azerbaijan, Taiwan and Japan. X-ray diffraction, visible / near infrared reflectance spectroscopy and Raman spectroscopy show that the samples are dominated by mixed-layer smectite clays, along with quartz, calcite and pyrite. Thin section analysis by optical and scanning electron microscopy confirms the mineral identifications. These samples also contain chemical and morphological biosignatures, including common microfossils, with evidence of partial replacement by pyrite. The bulk samples contain approximately 1 wt% total organic carbon and 0.4 mg / gm volatile hydrocarbons. The thousands of features in Acidalia Planitia cited as analogous to terrestrial mud volcanoes clearly represent an important element in the sedimentary record of Mars. Their location, in the distal depocenter for massive Hesperian-age floods, suggests that they contain fine-grained sediments from a large catchment area in the martian highlands. We have proposed these features as a new class of exploration target that can provide access to minimally-altered material from significant depth. By analogy to terrestrial mud volcanoes, these features may also be excellent sites for the sampling martian organics and subsurface microbial life, if such exist or ever existed.

Allen, Carlton C.↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no or some addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4 rocket engine fired on 1m2 area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon. The resulting computed values of plume impingement surface temperature, stagnation pressure, gas velocity, shear stress and heat flux were then matched as closely as possible in the Earth’s atmosphere in a sub-scale rocket engine GO2/CH4 test which was provided by Masten Space Systems in Mojave, California. The rocket engine was mounted on a test stand with vertical translation capabilities so that the landing operations of a lander could be simulated (Figure 1).The pad materials test coupons were placed at the surface of a large bin containing simulated lunar basalt regolith and subjected to a test firing as shown in Figure 2. The results of this testing will be presented with related findings, analysis and discussion.

Plume Surface Interaction↗

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗

Modeling CO 2 plume migration using an invasion-percolation approach that includes dissolution

Geological carbon dioxide (CO 2 ) sequestration has been proposed as a viable technique to decrease effective emissions of CO 2 into the atmosphere. However, the security of this sequestration is tied to our understanding of the long-term migration of CO 2 in subsurface. The dissolution of CO 2 in the reservoir brine is one of the main long-term trapping mechanisms. However, the assumptions used in large-scale reservoir simulations usually lead to an overestimation of the dissolution volume. We propose a modified approach based on the macroscopic invasion-percolation (MIP) theory that allows the dissolution of CO 2 into brine. We used a high-resolution geological model to compare the Darcy-, modified MIP-, and classic MIP-based simulation results. We observed a significant shrinkage in the nonaqueous plume volume when dissolution is considered during the MIP simulation. In the case of Darcy-based simulation, the plume was completely trapped inside the reservoir with limited migration even after a thousand-year simulation. On the other hand, the majority of the plume migrated out of the simulated reservoir in the case of MIP. Our approach provides more realistic estimation of the dissolution volume and nonaqueous plume extent while leveraging the computational efficiency enjoyed by MIP.

54 ENVIRONMENTAL SCIENCES↗

Low-Temperature Geothermal Geospatial Datasets: An Example from Alaska

This project is a component of a broader effort focused on geothermal heating and cooling (GHC) with the aim of illustrating the numerous benefits of incorporating GHC and geothermal heat exchange (GHX) into community energy planning and national decarbonization strategies. To better assist private sector investment, it is currently necessary to define and assess the potential of low-temperature geothermal resources. For shallow GHC/GHX fields, there is no formal compilation of subsurface characteristics shared among industry practitioners that can improve system design and operations. Alaska is specifically noted in this work, because heretofore, it has not received a similar focus in geothermal potential evaluations as the contiguous United States. The methodology consists of leveraging relevant data to generate a baseline geospatial dataset of low-temperature resources (less than 150 degrees C) to compare and analyze information accessible to anyone trying to understand the potential of GHC/GHX and small-scale low-temperature geothermal power in Alaska (e.g., energy modelers, communities, planners, and policymakers). Importantly, this project identifies data related to (1) the evaluation of GHC/GHX in the shallow subsurface, and (2) the evaluation of low-temperature geothermal resource availability. Additionally, data is being compiled to assess repurposing of oil and gas wells to contribute co-produced fluids toward the geothermal direct use and heating and cooling resource potential. In this work we identified new data from three different datasets of isolated geothermal systems in Alaska and bottom-hole temperature data from oil and gas wells that can be leveraged for evaluation of low-temperature geothermal resource potential. The goal of this project is to facilitate future deployment of GHC/GHX analysis and community-led programs and update the low-temperature geothermal resources assessment of Alaska. A better understanding of shallow potential for GHX will improve design and operations of highly efficient GHC systems. The deployment and impact that can be achieved for low-temperature geothermal resources will contribute to decarbonization goals and facilitate widespread electrification by shaving and shifting grid loads.

15 GEOTHERMAL ENERGY↗

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗

Sapflow and xylem water isotopes from Snodgrass Mountain, East River Watershed, Colorado USA

This dataset includes sapflux and stable water isotopes of soil water and xylem water for aspen, fir and spruce trees along the Snodgrass Mountain transect in the East River Watershed, Colorado USA. The purpose of generating this dataset was to understand: (1) the total flux of water being used by trees in the watershed, (2) separate the component of transpiration that was derived from recent precipitation vs. older water such as winter precipitation of groundwater and (3) understand how total water use and water sources for the trees varies between species and position on a hillslope. The data were collected from May 2019 until October 2022. The sap flux data were collected using ICT SFM1 sensors and are presented in both units of cm h^-1 and as kg h^-1 by multiplying the sap flux by the sapwood area of the tree. All sap flow data has been been corrected using estimates of wounding diameter, water content of wood and sap wood depth. The xylem water isotope data were collected approximately weekly from each of the trees instrumented with sap flux. The stems were collected and the water extracted using classic cryogenic methods. Isotope measurements were done using a Picarro 2140i analyzer. We provide an estimate of the Seasonal Origin Index for each measurements following Allen et al., 2019 (10.5194/hess-23-1199-2019) where values of -1 equate to trees relying on winter precipitation and +1 tree relying on summer precipitation. We also provide the simultaneous sap flux for each isotope measurement when this data was available. This is an update to an earlier data repository with the same name that only included 2019 data. The new dataset was posted in January 2023 and is inclusive of the original 2019 data but now includes 2021 and 2022 data. Please note there was an error in the units of transpiration in the original dataset. It was incorrectly listed as mm h-1 when it should have been kg h-1. This change was made in May 2024.

54 ENVIRONMENTAL SCIENCES↗

NASAs Evolvable Mars Campaign: Mars Moons Robotic Precursor

Human exploration missions to the moons of Mars are being considered within NASA's Evolvable Mars Campaign (EMC) as an intermediate step for eventual human exploration and pioneering of the surface of Mars. A range of mission architectures is being evaluated in which human crews would explore one or both moons for as little as 14 days or for as long as 500 days with a variety of orbital and surface habitation and mobility options being considered. Relatively little is known about the orbital, surface, or subsurface characteristics of either moon. This makes them interesting but challenging destinations for human exploration missions during which crewmembers must be able to effectively conduct scientific exploration without being exposed to undue risks due to radiation, dust, micrometeoroids, or other hazards. A robotic precursor mission to one or both moons will be required to provide data necessary for the design and operation of subsequent human systems and for the identification and prioritization of scientific exploration objectives. This paper identifies and discusses considerations for the design of such a precursor mission based on current human mission architectures. Objectives of a Mars' moon precursor in support of human missions are expected to include: 1) identifying hazards on the surface and the orbital environment at up to 50-km distant retrograde orbits; 2) collecting data on physical characteristics for planning of detailed human proximity and surface operations; 3) performing remote sensing and in situ science investigations to refine and focus future human scientific activities; and 4) prospecting for in situ resource utilization. These precursor objectives can be met through a combination or remote sensing (orbital) and in-situ (surface) measurements. Analysis of spacecraft downlink signals using radio science techniques would measure the moon's mass, mass distribution, and gravity field, which will be necessary to enable trajectory planning. Laser altimetry would precisely measure the moon's shape and improve the accuracy of radio science measurements. A telescopic imaging camera would map the moon at submeter resolution and photograph selected areas of interest at subcentimeter resolution and a visible and near-infrared (0.4-3.0 mm) imaging spectrograph would produce a global map of mineral composition variations at a resolution of tens of meters and maps of selected areas of interest at meter resolution. Additional remote sensing capabilities could include a thermal infrared imager (heat flow, thermal inertia, and grain size distributions), a gamma-ray and neutron detector (atomic composition), a ground-penetrating radar (internal structure), and a magnetometer and Langmuir probe (magnetic properties and plasma field). Once on the surface of Phobos or Deimos, necessary instrumentation would include a penetrometer (regolith compressive strength), a motion-imagery camera (to observe the penetrometer tests before, during, and after contact), a dust-adhesion witness plate and camera (dust levitation), a microimager (dust particle sizes and shapes), and an alpha-proton-X-ray, X-ray fluorescence, Mossbauer, or Raman spectrometer (atomic and mineral composition of surface materials) and an optional temperature probe (regolith thermal properties). A variety of robotic mission design options to enable both orbital and surface measurements are being considered that include fully integrated and modular approaches. In-situ measurements from at least one surface location would be required, with additional measurement locations possible through use of multiple landers, through propulsive relocation of a single lander, or through electromechanical surface translation by a walking or hopping lander vehicle, which could also serve to evaluate such mobility capabilities for subsequent human missions. Preliminary orbital analysis suggests that remote sensing would likely be performed while in a distant retrograde orbit around the target moon. Mission design options to enable characterization of both Mars’ moons in a single mission are also being studied.

Gernhardt, Michael L.↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimating Watershed Subsurface Permeability From Stream Discharge Data using Deep Neural Networks, Frontiers in Earth Science: Dataset

This data package contains watershed modeling inputs and outputs as well as deep neural networks training and testing results used in "Estimating Watershed Subsurface Permeability From Stream Discharge Data using Deep Neural Networks" (Cromwell et al., 2021). We train various deep neural network models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models. For detailed information regarding watershed model description and DNNs setup, please refer to Cromwell et al., 2021. The deep_learning.zip file contains the inputs and outputs for the different deep neural networks. The watershed_model.zip file contains the inputs and outputs for watershed simulation using ATS. The figures.zip file contains raw figures and their corresponding scripts. For a more detailed description, please see the README.md file within each compressed file.

54 ENVIRONMENTAL SCIENCES↗

An Empirical Algorithm to Map Perennial Firn Aquifers and Ice Slabs Within the Greenland Ice Sheet Using Satellite L-band Microwave Radiometry

Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature () imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

Julie Z Miller↗

Radiation and Nuclear Technology in Planetary Cave Environments

Many technological and environmental challenges must be resolved to enable successful lunar habitation and exploration, and to maximize scientific return. This presentation is intended as a preliminary discussion on three such areas of interest with the intent to identify likely areas of mutual benefit for collaboration and data sharing between the space nuclear and radiological disciplines with the planetary cave community. First is the previously identified and studied possible application of lunar lava tubes or pits to reduce the solar and cosmic radiation environment burden on crew health and hardware reliability. This is an area of continued interest, and should be kept in the forefront of discussion during the selection of sites for potential exploration or long-term habitation. In theory, the prospect of using existing morphology presents an enticing opportunity to reduce the requirements on landed mass or construction while also reducing dose or fluence of harmful natural radiation. However, such discussions should also include the practical implications of relocating habitation hardware and personnel from a landing site to a subsurface location. Inherent risks associated with landing in proximity to the relevant terrain must also be considered and weighed against those imposed by the natural radiation environment. Second is a discussion on the means by which a long-term habitation module or base of operations may be powered throughout the lunar day-night cycle. Fission surface power presents an opportunity to establish round-the-clock power in any lunar surface environment, including permanently shadowed regions where solar panels cannot operate, or in any other non-polar region where solar radiation is available for no more than two weeks per four-week cycle. Among the drawbacks of fission surface power is the need to either land and co-locate a heavy shield to minimize radiological consequences to personnel or equipment, land and operate construction equipment capable of restructuring the in-situ regolith to provide appropriate shielding, or to make use of existing topology features (e.g. craters or pits) to serve as pre-formed radiation barriers. A realistic assessment of the practicality of this third option should involve the selection of candidate features from existing surveys, and assess the effectiveness of the approach using modern radiation transport methodologies. Consideration must be made for the needs to reject waste heat in the thermal power conversion process, which typically requires the use of large area radiators. Energy must be exchanged from the reactor and power conversion system to such radiators, and the impacts on thermal efficiency and secondary scatter of nuclear radiation must be considered. Further consideration should be made for the implications of placement of a power source such that it does not sacrifice valuable scientific opportunity. Placement of a reactor is likely to thereafter prevent personnel access within that feature. Third is a discussion on the use of portable radioisotope power systems (RPS) within confined spaces, both with respect to the dose effects of emitted nuclear radiation and also performance associated with emitted thermal radiation. Radioisotope power systems provide a unique capability to power objects with no reliance upon solar radiance. However, their useable power production relies upon the flow of heat from ‘hot side’ (decaying radioisotope) to ‘cold side’ (radiators emitting heat to space). In the case of a confined volume in vacuum such as a lunar lava tube, that radiated energy will absorb into the wall, which is already above the ‘cold background’ temperature of dark space, and then gradually rises in temperature. That temperature rise, and its implications on availability of useable electrical power, will depend upon factors such as enclosure volume, thermal emissivity/absorptivity, and thermal conductivity through the depth of absorbing media. The advantages of nuclear technology may prove to enable unique scientific exploration opportunities, both above the lunar surface and below. In turn, the unique advantages of lunar lava tubes or pits may present unique opportunities to enhance the application of these technologies and to mitigate the effects of radiation from both natural and technological sources. A discussion on these pairings may prove to be valuable to all, in our endeavor to explore our neighboring worlds.

Radiation↗

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

The large thermal gradients and high subsurface temperatures of the western region of the U.S. hold great potential for the implementation of enhanced geothermal systems (EGS). The development of these potential EGS resources requires stimulation of the reservoir to enhance permeability and it has been widely reported that a substantial amount of water will be required should conventional hydraulic stimulation be used. This presents a huge challenge and a high risk to the geothermal development because the water stress1 in these areas is already high or extremely high. The use of foam, a gas/liquid mixture predominantly composed of gas, in fracturing is considered and explored in this project as a potential approach to address water concerns with hydraulic stimulation in the development of EGS. This project, led by Oak Ridge National Laboratory (ORNL) in collaboration with Temple University, was awarded in an open lab call in 2018, and was part of the DOE GTO waterless stimulation initiative. The goal of the project was to demonstrate the feasibility of foam fracturing for EGS development through two primary tasks: Task 1: Laboratory study of the effectiveness of foam fracturing for representative geological materials, including cyclic pressurization using foam (led by ORNL) and Task 2: High temperature foam material selection and characterization (led by Temple University). In FY19, ORNL finished the critical review on serval issues associated with foam fracturing and the implementation of the proposed tasks in a lab study (Wang, et al., 2019), and completed the foam fracturing testing using cement as a model material (Wang, et al, 2020a). The work at ORNL was geared up to develop a brand-new foam testing system in FY20. The purchase of main components for the new system was finished in the first half of the FY20. The assembly of the foam testing system and foam fracturing testing were completed in the second half of the FY20 (Wang, et al., 2021a). Task 1 required the development of a test system which can be used to perform hydraulic fracturing of geological specimens with both water and foamed liquids at pressure up to 6,000 psi (41.4 MPa). The system possesses several capabilities that conventional injection systems lack for hydraulic fracturing. In addition to its ability to generate foam with controlled quality, it is capable of cycling pressure levels between specified values up to frequencies of 50 Hz. The latter capability was developed to evaluate the hypothesis that cyclic loading of samples would produce enhanced fracturing. The system consists of two sections: one for foam generation and another for foam injection. The foam is generated through separate control and pressurization of liquid and gas phases with controlled flow rates. The injection section is equipped with a low-flow Coriolis flow that monitors the density of foam to ensure the injection is in the range of target foam quality2. Experimental results of foam fracturing are reported for cylindrical granite specimens using water and aqueous N2 foam as the fracturing fluids. All experiments were performed for unconfined conditions. The effects of injection mode (i.e., monotonic vs cyclic pressurization) on breakdown pressure and failure response sample were investigated using water alone as a fracturing fluid and foams with a range of compositions. It was found that in the case of monotonic injection, the breakdown pressure of granite specimens tended to be slightly higher when fracturing with foam. Additionally, with a foam quality of 90%, the water use can be reduced by 50 to 84%, depending on hole size. On the other hand, it was observed that the breakdown pressure can be brought down to 70% of the monotonic breakdown pressure by using low cycle fatigue. Finally, discussions are presented regarding injectivity and water use reduction.

15 GEOTHERMAL ENERGY↗