Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,153 records · Page 64

Exploring the Use of Non‐Invasive Drone‐Based Ground‐Penetrating Radar (GPR) to Characterize Biogenic Gas Dynamics in Subtropical Peat Soils

Peat soils are a critical component of the global carbon cycle as natural producers of biogenic greenhouse gases (e.g., methane and carbon dioxide) that accumulate within the soil and are released to the atmosphere. Previous studies have showed the ability of ground-based minimally-invasive geophysical methods such as ground-penetrating radar (GPR) to characterize carbon dynamics in peat soils. However, ground-based GPR is limited by scale of measurement and soil disturbance potentially altering gas releases during deployment. Here, we explore the potential of drone-based GPR for identification of hot spots and hot moments of gas accumulation and release in subtropical soils. Here, we collected drone-based GPR data sets across two grids (∼17,500 m 2 ) in the Everglades during January (dry season), September, and November (wet season) of 2023 to characterize peat thickness and seasonal variability of gas content. Results show that drone-based GPR is effective and efficient for: (a) capturing the temporal variation of in situ biogenic gas content in peat soils with changes between 1% and 25 % volumetric gas content over repeatable grids; (b) inferring a total peat thickness between 0.8 and 1.2 m; and (c) estimating flux releases of 63 and 135 mg CH 4 m −2 day −1 for specific locations and periods that are strikingly consistent with our coincident gas trap measurements. This work also indicates that (a) spatial distribution of gas content in the Everglades is strongly controlled by landscape morphology such as ridges and sloughs and (b) the temporal variation of gas content is seasonal with increased gas production during the wet season.

54 ENVIRONMENTAL SCIENCES↗

The Large Impact Process Inferred from the Geology of Lunar Multiring Basins

The study of the geology of multiring impact basins on the Moon over the past ten years has given us a rudimentary understanding of how these large structures have formed and evolved on the Moon and other bodies. Two-ring basins on the Moon begin to form at diameters of about 300 km; the transition diameter at which more than two rings appear is uncertain, but it appears to be between 400 and 500 km in diameter. Inner rings tend to be made up of clusters or aligned segments of massifs and are arranged into a crudely concentric pattern; scarp-like elements may or may not be present. Outer rings are much more scarp-like and massifs are rare to absent. Basins display textured deposits, interpreted as ejecta, extending roughly an apparent basin radius exterior to the main topographic rim. Ejecta may have various morphologies, ranging from wormy and hummocky deposits to knobby surfaces; the causes of these variations are not known, but may be related to the energy regime in which the ejecta are deposited. Outside the limits of the textured ejecta are found both fields of satellitic craters (secondaries) and light plains deposits. Impact melt sheets are observed on the floors of relatively unflooded basins. Samples of impact melts from lunar basins have basaltic major-element chemistry, characterized by K, rare-earth elements (REE), P, and other trace elements of varying concentration (KREEP); ages are between 3.8 and 3.9 Ga. These lithologies cannot be produced through the fusion of known pristine (plutonic) rock types, suggesting the occurrence of unknown lithologies within the Moon. These melts were probably generated at middle to lower crustal levels. Ejecta compositions, preservation of pre-basin topography, and deposit morphologies all indicate that the excavation cavity of multiring basins is between about 0.4 and 0.6 times the diameter of the apparent crater diameter. Basin depths of excavation can be inferred from the composition of basin ejecta. A variety of mechanisms has been proposed to account for the formation of basin rings but none of them are entirely plausible. Mechanisms can be divided into two broad groups: (1) forcible uplift due to fluidization of the target; (2) concentric, brittle, fracturing and failure of the target, on regional (megaterraces) to global scales (lithospheric fracturing). Most basin rings are spaced at a constant factor on all planets. Evidence supports divergent ringforming models, so it may be that the ring-locating mechanism differs from the ring-forming mechanism. Thus, large-scale crustal foundering (megaterracing) could occur along concentric zones of weakness created by some type of resonant wave mechanism (fluidization and uplift); such immediate crustal adjustment could then be followed by long-term adjustment of the fractured lithosphere.

Spudis, Paul D.↗

The ESA Hera Mission: Detailed Investigation of the NASA Dart Impact Outcome and Characterization of the Binary Asteroid Didymos

The Hera mission is in development in the Space Safety Program of the European Space Agency (ESA) in collaboration with JAXA for launch in October 2024. It will rendezvous with the binary asteroid (65803) Didymos in early 2027 and over the following 6 months will investigate the properties of Didymos and its small moon Dirmorphos, including the effects of the successful impact of the NASA DART probe. DART was launched on 24 November 2021 at 06:21 UTC and successfully impacted at about 6.1 km/s on the 160 meter-size Dimorphos on 26 September 2022 at 23:14 UTC. Two weeks before impact, it deployed the Light Italian Cubesat for Imaging of Asteroids (LICIACube) that provided images during the few minutes following the impact. A campaign of observations from Earth and space provided distant images of the event, which were used to determine a reduction of 33 min in the orbital period of Dimorphos around Didymos. The Hera mission will investigate the outcome of the DART impact in detail. In particular, it will determine the effects of the impact on Dimorphos’ surface, e.g., whether it produced a crater, and if so its properties, and/or whether it led to global reshaping of Dimorphos. It will accurately measure the mass of Dimorphos, which will allow quantification of the momentum enhancement factor, which is a crucial parameter to evaluate the efficiency of the kinetic impactor technique. Hera will also determine in detail the physical, thermal and compositional properties of the asteroid, including first time measurements of the internal properties, which have a critical influence on the impact outcome. With Hera and DART together, we will have a fully documented deflection test and impact experiment that will allow us to verify our numerical impact models, to be able to extrapolate this deflection test to other asteroids and to improve our understanding of the complex geological properties and response of small asteroids. Hera consists of a mother spacecraft and two cubesats. It will perform the first rendezvous with a binary asteroid, arriving at Didymos about 4 years after DART. Hera will not only contribute significantly to the planetary defense effort by documenting entirely the outcome of the DART impact, it will also greatly advance our understanding of processes in the very low gravity regime of a small asteroidv (Dimorphos) as well as our knowledge of binary asteroids and asteroid geophysics. The mission development is progressing nominally, with the Critical Design Review having been passed in Nov. 2022. The various working groups of the Hera Science Team are working intensively to support this development and be best prepared for the surprises that Hera will offer during its visit to Didymos in 2027.

Asteroid deflection↗

The Effect of Indian Ocean Temperature on the Pacific Trade Winds and ENSO

Abstract A notable shift in the El Niño‐Southern Oscillation (ENSO) has been observed in the early 21st century, characterized by an increased prevalence of Central Pacific (CP) events and strengthened Pacific trade winds. This shift may be attributed to the warming tropical Indian Ocean (TIO). To investigate this, we conduct perturbation experiments using the Insitut Pierre Simon Laplace climate model and nudge TIO surface temperatures to induce warming or cooling effects. Our findings reveal that TIO warming (or cooling) leads to amplified (weakened) mean trade winds and surface warming (cooling) in the Pacific region. Surprisingly, ENSO variability increases in both TIO cooling and warming scenarios. This result is linked to stronger positive feedbacks and a less stable Bjerknes index for either TIO forcing. Additionally, we find that TIO warming leads to more frequent CP events, meridional widening of wind anomalies, and broadening of the ENSO power spectrum toward lower frequencies.

Geology↗

Multi-Temporal Predictive Modelling of Sorghum Biomass Using UAV-Based Hyperspectral and LiDAR Data

High-throughput phenotyping using high spatial, spectral, and temporal resolution remote sensing (RS) data has become a critical part of the plant breeding chain focused on reducing the time and cost of the selection process for the “best” genotypes with respect to the trait(s) of interest. In this paper, the potential of accurate and reliable sorghum biomass prediction using visible and near infrared (VNIR) and short-wave infrared (SWIR) hyperspectral data as well as light detection and ranging (LiDAR) data acquired by sensors mounted on UAV platforms is investigated. Predictive models are developed using classical regression-based machine learning methods for nine experiments conducted during the 2017 and 2018 growing seasons at the Agronomy Center for Research and Education (ACRE) at Purdue University, Indiana, USA. The impact of the regression method, data source, timing of RS and field-based biomass reference data acquisition, and the number of samples on the prediction results are investigated. R2 values for end-of-season biomass ranged from 0.64 to 0.89 for different experiments when features from all the data sources were included. Geometry-based features derived from the LiDAR point cloud to characterize plant structure and chemistry-based features extracted from hyperspectral data provided the most accurate predictions. Evaluation of the impact of the time of data acquisition during the growing season on the prediction results indicated that although the most accurate and reliable predictions of final biomass were achieved using remotely sensed data from mid-season to end-of-season, predictions in mid-season provided adequate results to differentiate between promising varieties for selection. The analysis of variance (ANOVA) of the accuracies of the predictive models showed that both the data source and regression method are important factors for a reliable prediction; however, the data source was more important with 69% significance, versus 28% significance for the regression method.

09 BIOMASS FUELS↗

Using Multiple Geophysical Methods to Refine a Stratigraphic Conceptual Site Model at a Nuclear Waste Site

At contaminated sites, there is a critical reliance on conceptual site models of which a key component is the geologic framework model (GFM), which is used to build predictive fate and transport modeling in support of remediation goals. Here, we detail an ongoing study at a nuclear waste site (Hanford Site, Washington, USA) where surface geophysical methods were used to complement existing borehole information and to help site new boreholes to update and refine the GFM. Ground-based electrical resistivity tomography (ERT), time-domain electromagnetics (TEM), and seismic data was acquired on the Hanford Site’s Central Plateau, where the presence of a hydraulically transmissive subsurface feature, or paleochannel, was suspected but uncertain based on contaminant concentrations at sparsely located boreholes. To assess and compare ERT, TEM, and seismic geophysical datasets, co-located data was acquired in an area where the existing GFM was more certain. Surrounding the suspected paleochannel, ERT and seismic refraction tomography showed similar subsurface structures consistent with a channelized feature. While TEM interpretation was limited in the suspected paleochannel area due to a thick resistive top layer and high electromagnetic noise, TEM showed greater utility in a different area where these factors were less prevalent. Based on the results of quasi-3D ERT inversions, we propose updates to the GFM, including identifying a paleochannel. Further, we evaluate the utility of geophysical methods and examine lessons learned that will be used as a basis for continued characterization efforts to support site management decisions and implementation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

The Meteoritic Component in Impact Deposits

This proposal requested support for a broad-based research program designed to understand the chemical and mineralogical record of accretion of extraterrestrial matter to the Earth. The primary goal of this research is to study the accretion history of the Earth, to understand how this accretion history reflects the long-term flux of comets, asteroids, and dust in the inner solar system and how this flux is related to the geological and biological history of the Earth. This goal is approached by seeking out the most significant projects that can be attacked utilizing the expertise of the PI and potential collaborators. The greatest expertise of the PI is the analysis of meteoritic components in terrestrial sediments. This proposal identifies three primary areas of research, involving impact events in the early Archean (3.2 Ga), the late Eocene (35 Ma) and the late Pliocene (2 Ma). In the early Archean we investigate sediments that contain the oldest recorded impacts on Earth. These are thick spherule beds, three of which were deposited within 20 m.y. If these are impact deposits the flux of objects to Earth at this time was much greater than predicted by current models. Earlier work used Cr isotopes to prove that one of these contain extraterrestrial matter, from a projectile with Cr isotopes similar to CV chondrites. We planned to expand this work to other spherule beds and to search for additional evidence of other impact events. With samples from D. Lowe (Stanford Univ.) the PI proposed to screen samples for high Ir and Cr so that appropriate samples can be provided to A. Shukolyukov for Cr-isotopic analyses. This work was expected to provide evidence that at least one interval in the early Archean was a period of intense bombardment and to characterize the composition of objects accreted. The late Eocene is also a period of intense bombardment with multiple spherule deposits and two large craters. Farley et al. (1998) demonstrated an increased (3)He flux to marine sediments that was attributed to an increase in interplanetary dust due to a shower of comets invading the inner solar system. We planned to detect a change in the Cr-isotopic composition in the flux of fine-grained extraterrestrial matter to ocean sediments. This would provide evidence for the comet shower hypothesis. We planned attempt to locate late Eocene impact deposits in a new high resolution, hi-latitude site that had the potential for excellent preservation and new information of the sources and effects of these impacts. Although most impacts on Earth occur in deep-ocean basins, only one such event is known - the late Pliocene impact of the Eltanin asteroid. The ejecta from this impact includes Ir- rich impact melt, spherules, and actual meteorites from the km-sized mesosiderite asteroid. Through a study of the ejecta, we can learn about the formation, distribution, alteration and preservation of Ir-rich deposits. We can also learn about meteorite survival during hypervelocity impacts and we can study pieces of a km-sized object to learn more about the mesosiderite parent body. Analyses of sediment cores and geophysical exploration of the impact site can further our understanding of the processes involved in deep-ocean impacts and potential effects on the terrestrial climate and biosphere. We planned analyses of this ejecta, a search for ejecta 5000 km distant from the impact area, and a new oceanographic expedition to study and sample the impact site. In addition to these three specific areas of research, the PI planned to remain flexible and available to exploit new opportunities presented by new discoveries, and to engage in new collaborations if other researchers require his expertise to develop new projects that fit within the objectives of the overall research program.

Kyte, Frank T.↗

Transformational faulting in Mn 2 GeO 4 from olivine to wadsleyite structure: Implications for physical mechanism of deep-focus earthquakes

High-pressure and temperature deformation experiments interfaced with acoustic emission (AE) monitoring have been conducted to study transformational faulting in Mn 2 GeO 4 olivine, which transforms to the β phase, isostructural to wadsleyite. Metastable Mn 2 GeO 4 olivine exhibits a marked embrittlement behavior at temperatures between 800 and 1100 K, emitting numerous AEs. At each temperature, brittle deformation is characterized by a two-stage process: (1) a “preparation” stage with numerous diffusedly located low-magnitude AEs and large b values (>2), and (2) a failure stage where larger-magnitude AEs form a planar distribution with b values about 1. Microstructure analysis reveals extensive kink band development in olivine grains in the recovered samples. Kink band boundaries (KBBs), with a typical thickness of ∼100 nm, are filled with a nanometric β-Mn 2 GeO 4 “gouge”. A dense array of secondary shear localizations is often present within the kink bands, suggesting significant shear deformation therein. The combined observations suggest that faulting in metastable Mn 2 GeO 4 olivine is a self-similar process, from grain-scale to the sample-scale. Both observed embrittlement behavior and the microstructure of metastable Mn 2 GeO 4 olivine are essentially identical to those in Mg 2 GeO 4 olivine we have reported previously, indicating that the physical mechanism of faulting in metastable olivine is insensitive to the specific crystallographic structure of the high-pressure phase. The low b values (about 1) observed in the faulting process in our experiments are similar to those of deep focus earthquakes in cold subduction zones. Our observed mechanism explains deep focus seismicity in cold metastable mantle wedges, provided that the self-similarity assumption holds to geological scales.

58 GEOSCIENCES↗

Improved dryland carbon flux predictions with explicit consideration of water-carbon coupling

Dryland ecosystems are dominant influences on both the trend and interannual variability of the terrestrial carbon sink. Despite their importance, dryland carbon dynamics are not well-characterized by current models. Here, we present DryFlux, an upscaled product built on a dense network of eddy covariance sites in the North American Southwest. To estimate dryland gross primary productivity, we fuse in situ fluxes with remote sensing and meteorological observations using machine learning. DryFlux explicitly accounts for intra-annual variation in water availability, and accurately predicts interannual and seasonal variability in carbon uptake. Applying DryFlux globally indicates existing products may underestimate impacts of large-scale climate patterns on the interannual variability of dryland carbon uptake. We anticipate DryFlux will be an improved benchmark for earth system models in drylands, and prompt a more sensitive accounting of water limitation on the carbon cycle.

54 ENVIRONMENTAL SCIENCES↗

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 2—Practice and Demonstration

Summary This study continues the work of presenting a novel approach for making petrophysical assessments of tight core samples. This method, the full-immersion pressure-pulse decay (Hannon 2019), involves applying a rapid increase in pressure in a chamber surrounding the entire outer surface area of a cylindrical sample, shutting the system in, and monitoring the pressure decay in the chamber as it reaches a new equilibrium. A precursor article covered the numerical simulator designed to model flow through the sample, demonstrating its performance and accuracy in addition to providing a first-order comparison between the speed and shape of the pressure-decay responses of the full-immersion method with those of other similar transient methods. This study covers the parameter-estimation procedure and experimental verification through a proof-of-concept laboratory investigation. The investigations provided here demonstrate that under appropriate, achievable experimental conditions, the pressure data can be analyzed in such a way that returns an estimate of the porosity and apparent permeabilities both parallel and perpendicular to bedding from a single test performed on a single cylindrical sample. After determining these experimental conditions (the uniqueness window), this report outlines a data-inversion strategy to estimate the petrophysical properties (porosity, horizontal permeability, and vertical permeability) from each test. This strategy is put to the test through comparisons with measurements performed by a commercial core laboratory. A common set of samples recovered from an outcrop of a tight-gas sandstone formation were investigated using the full-immersion method, and their results are compared with those from conventional steady-state measurement procedures performed by the commercial laboratory. Comparisons between petrophysical characterizations of these samples, which had permeabilities between 25 nd and 2.3 μd, demonstrated close agreement in most cases. However, whereas steady-state measurements performed at the professional laboratory required 4 to 5 hours of testing time per measurement of a single permeability, similar assessments using the full-immersion technique, requiring approximately 5 to 10 minutes to complete, returned estimates of the horizontal and vertical permeability simultaneously. Additional analyses are provided to determine principal reasons of discrepancies in instances where agreement was not as strong. Based on lessons learned from these experiences, the report closes with suggestions on areas of improvement in the experimental approach. Once complete, these developments should propel this technology to fill a critical need to determine petrophysical properties (porosity and permeability) of tight rocks in a time-efficient manner and in a way that does not compromise their accuracy.

Energy & Fuels↗

Paleobiology of the Mesoproterozoic Billyakh Group, Anabar Uplift, northern Siberia

Silicified peritidal carbonates of the Mesoproterozoic Kotuikan and Yusmastakh Formations, Anabar Uplift, northeastern Siberia, contain exceptionally well-preserved microfossils. The assemblage is dominated by ellipsoidal akinetes of nostocalean cyanobacteria (Archaeoellipsoides) and problematic spheroidal unicells (Myxococcoides); both are allochthonous and presumably planktonic. The assemblage also includes distinctive mat-forming scytonematacean and entophysalidacean cyanobacteria, diverse short trichomes interpreted as cyanobacterial hormogonia or germinated akinetes, rare longer trichomes, and several types of colonial unicells. Although many taxa in the Kotuikan-Yusmastakh assemblage are long-ranging prokaryotes, the overall character of the assemblage is distinctly Mesoproterozoic, with its major features shared by broadly coeval floras from Canada, China, India, and elsewhere in Siberia. Microfossils also occur in middle to inner shelf shales of the Ust'-Il'ya and lower Kotuikan Formations. Leiosphaerid acritarchs (up to several hundred microns in diameter) characterize this facies. As in other Mesoproterozoic acritarch assemblages, acanthomorphic and other complex forms that typify Neoproterozoic assemblages are absent. The combination in Billyakh assemblages of exceptional preservation and low eukaryotic diversity supports the hypothesis that nucleated organisms diversified markedly near the Mesoproterozoic-Neoproterozoic boundary. The assemblages also demonstrate the antiquity of cyanobacteria capable of cell differentiation and suggest the importance of both changing peritidal substrates and evolving eukaryotes in determining stratigraphic patterns of Proterozoic prokaryotes. The permineralized assemblage contains 33 species belonging to 17 genera. Ten new species or new combinations are proposed: Archaeoellipsoides costatus n. sp., A. elongatus n. comb., A. dolichos n. comb., A. minor n. nom., A. crassus n. comb., A. major n. comb., A. bactroformis n. sp., Veteronostocale medium n. sp., Filiconstrictosus cephalon n. sp., and Partitiofilum yakschinii n. sp.

Non-NASA Center↗

The distribution of ground ice on Mars

A wealth of geologic evidence indicates that subsurface water ice has played an important role in the evolution of Martian landforms. Theoretical models of the stability of ground ice show that in the near-surface regolith ice is currently stable at latitudes poleward of about +/- 40 deg and below a depth of a few centimeters to a meter, with some variations with longitude. If ice is not previously present at a particular location where it is stable, atmospheric water will diffuse into the regolith and condense as ice, driven by the annual subsurface thermal oscillations. The lower boundary of this ice deposit is found to occur at a depth (typically a few meters) where the annual thermal oscillations give way to the geothermal gradient. In the equatorial regions near-surface ice is currently not stable, resulting in the sublimation of any existing ice and subsequent loss to the atmosphere. However, subliming ice might be maintained at a steady-state depth, where diffusion and loss to the atmosphere are balanced by resupply from a possible deeper source of water (either deeper, not yet depleted, ice deposits or ground water). This depth is typically a few tens to hundreds of meters and depends primarily on the surface temperature and the nature of the geothermal gradient, being deeper for a higher surface temperature and a lower geothermal gradient. Such an equatorial deposit is characterized by the regolith ice content being low nearer the surface and increasing with depth in the deposit. Oscillations in the orbit will affect this picture of ground ice in two ways: by causing periodic changes in the pattern of near-surface stability and by producing subsurface thermal waves that may be capable of driving water ice deeper into the regolith.

Mellon, M. T.↗

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗

Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Pierazzo, E.↗

Microbial Lipid and C Isotopic Biosignatures of a Unique Community at Grand Prismatic Spring, Yellowstone National Park

The microbial communities found in modern hot springs are considered analogs to ones that may have existed in hydrothermal systems on the early Earth and possibly Mars. Our goal was to characterize the microbial biosignatures and to assess the preservation of organic matter in the silica-depositing Grand Prismatic Spring in Yellowstone National Park. This study combines 16S rRNA surveys, lipid biomarkers, and C isotopes to query, "Who's there and what are they doing?" On the edge of the approximately 90 m diameter blue vent pool (56.1 C, pH 8.5), a floating green streamer community grew over a benthic pink community. The membrane lipids in the green streamers and pink mat were composed of unusual ester-linked fatty acids, indicating the presence of novel bacterial groups. In particular, we discovered a series of 2-methyl and 2,X-dimethyl phospholipid fatty acids (C18-22). We are attempting to use the 16S rRNA surveys to link these compounds to source organisms. Wax esters, biomarkers for Chloroflexi, were present in both communities, but displayed different profiles. A higher proportion of branched wax esters were found in the green streamers, and were associated with a relatively high concentration of long-chain di- and trienes (C29-31). This suggests that Chloroflexus primarily grew in the green streamers, while a pink mat of Roseiflexus grew on the sinter substrate underneath. Cyanobacterial alkanes were found in the green streamers (n-C17, 7-, 6- and 5-monomethyl-C17, 7,11-dimethyl-C17, n-C19, n-C19:1). We also detected a series of monoalkylglycerylethers and geologically relevant hopanoids in both communities. Carbon isotope analyses indicated that Chloroflexus was growing photoheterotrophically using cyanobacterial photosynthate. Roseiflexus also traditionally grows photoheterotrophically, but the C isotopic signatures of the lipids in the pink mat were approximately 10 %0 lighter than the cyanobacterial and Chloroflexus lipids, indicating a potentially novel metabolic mechanism or prior secondary reworking of substrates before reaching Roseiflexus. This arrangement of Synechococcus-Chloroflexus green streamers floating over a benthic pink community of Roseiflexus is different from the classical laminated Synechococcus-Chloroflexi mats at Octopus Hot Spring.

Jahnke, Linda↗

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗