Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗
NASA NTRS1991

Geochemistry and petrogenesis of proterozoic mafic dykes in north Kerala, southwestern Indian Shield - Preliminary results

Mafic dyke intrusions occur in three distinct orientations ( NNW-SSE and NW-SE and NE-SW) in north Kerala regions, the southwestern part of the Indian Shield. Dating of two NNW-SSE trending dykes by K-Ar method has yielded Middle Proterozoic ages (ca. 1660 Ma and ca. 1420 Ma respectively). Our initial geochemical results on these dyke rocks (0.65-0.15 wt pct K2O, 0.37-0.38 wt pct P2O5, 3.30-1.00 wt pct TiO2, 11-1 p.p.m. Rb, 250-90 p.p.m. Sr, 230-40 p.p.m. Ba, 160-40 p.p.m. Zr, and 30-10 x chondrite rare earth elemental abundances) indicate a transitional character between abyssal and plateau tholeiites. Petrogenetic modelling suggests that the dyke compositions have been derived by different degrees of partial melting of a heterogenous source mantle. The mantle sources with accessory amphibole and/or garnet, variably enriched in LREE and LIL elements, are compatible with the observed geochemical data.

Radhakrishna, T.↗
DOE OSTI2022

Quantifying pH buffering capacity in acidic, organic-rich Arctic soils: Measurable proxies and implications for soil carbon degradation

Dynamic pH change promoted by biogeochemical reactions in Arctic tundra soils can be a major control on the production and release of CO 2 and CH 4 , which contribute to rising global temperatures. Large quantities of soil organic matter (SOM) in these soils are susceptible to microbial decomposition, leading to pH changes during permafrost thaw. Soil pH buffering capacity (β) modulates the extent of pH change but has not been thoroughly studied and represented in predictive ecosystem scale biogeochemical models in Arctic tundra soils. In this study, we generated titration curves for 21 acidic tundra soils from three Arctic sites across northern Alaska, United States of America. Geochemical and hydrological soil properties were evaluated, and correlations with β were developed. Strong correlations between β and both gravimetric water content (Θ g ) (R 2 = 0.847, p < 0.001) and soil water retention (SWR) (R 2 = 0.849, p = 0.001) indicate that the ability of soil to retain water could be associated with its buffering properties. Correlations between β and soil organic carbon (SOC) and cation exchange capacity (CEC) were also explored, and relationships to SWR are discussed. These correlations were then used with existing soil databases reporting SOC, CEC, and SWR to estimate β across Alaska soils. We further demonstrated the quantitative relationships between β and the simulated rates of biogeochemical reactions and show that lower β leads to higher soil pH and more CH 4 production. Overall, our study provides simple proxies for β in Arctic soils and highlights the importance and implications of representing soil buffering in predictive models, thereby enabling quantitative coupling between pH dynamics associated with biogeochemical reactions. Integrating β into predictive models of Arctic biogeochemical cycling may reduce model uncertainty and further our understanding of permafrost SOM degradation accelerated by warming.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Comparative Magma Oceanography

The question of whether the Earth ever passed through a magma ocean stage is of considerable interest. Geochemical evidence strongly suggests that the Moon had a magma ocean and the evidence is mounting that the same was true for Mars. Analyses of martian (SNC) meteorites have yielded insights into the differentiation history of Mars, and consequently, it is interesting to compare that planet to the Earth. Three primary features of Mars contrast strongly to those of the Earth: (i) the extremely ancient ages of the martian core, mantle, and crust (about 4.55 b.y.); (ii) the highly depleted nature of the martian mantle; and (iii) the extreme ranges of Nd isotopic compositions that arise within the crust and depleted mantle. The easiest way to explain the ages and diverse isotopic compositions of martian basalts is to postulate that Mars had an early magma ocean. Cumulates of this magma ocean were later remelted to form the SNC meteorite suite and some of these melts assimilated crustal materials enriched in incompatible elements. The REE pattern of the crust assimilated by these SNC magmas was LREE enriched. If this pattern is typical of the crust as a whole, the martian crust is probably similar in composition to melts generated by small degrees of partial melting (about 5%) of a primitive source. Higher degrees of partial melting would cause the crustal LREE pattern to be essentially flat. In the context of a magma ocean model, where large degrees of partial melting presumably prevailed, the crust would have to be dominated by late-stage, LREE-enriched residual liquids. Regardless of the exact physical setting, Nd and W isotopic evidence indicates that martian geochemical reservoirs must have formed early and that they have not been efficiently remixed since. The important point is that in both the Moon and Mars we see evidence of a magma ocean phase and that we recognize it as such. Several lines of theoretical inference point to an early Earth that was also hot and, perhaps, mostly molten. The Giant Impact hypothesis for the origin of the Moon offers a tremendous input of thermal energy and the same could be true for core formation. And current solar system models favor the formation of a limited number of large (about 1000 km) planetesimals that, upon accreting to Earth, would cause great heating, being lesser versions of the Giant Impact. Several lines of geochemical evidence do not favor this hot early Earth scenario. (i) Terrestrial man-tle xenoliths are sometimes nearly chondritic in their major element compositions, suggesting that these rocks have never been much molten. Large degrees of partial melting probably promote differentiation rather than homogenization. (ii) Unlike the case of Mars, the continental crust probably did not form as a highly fractionated residual liquid from a magma ocean (about 99% crystallization), but, rather, formed in multiple steps. [The simplest model for the formation of continental crust is complicated: (a) about 10% melting of a primitive mantle, making basalt; (b) hydrothermal alteration of that basalt, converting it to greenstone; and (c) 10% partial melting of that greenstone, producing tonalite.] This model is reinforced by the recent observation from old (about 4.1 b.y.) zircons that the early crust formed from an undepleted mantle having a chondritic Lu/Hf ratio. (iii) If the mantle were once differentiated by a magma ocean, the mantle xenolith suite requires that it subsequently be homogenized. The Os isotopic compositions of fertile spinel lherzolites place constraints on the timing of that homogenization. The Os isotopic composition of spinel lherzolites approaches that of chondrites and correlates with elements such as Lu and Al. As Lu and Al concentrations approach those of the primitive mantle, Os isotopic compositions approach chondritic. The Re and Os in these xenoliths were probably added as a late veneer. Thus, the mantle that received the late veneer must have been nearly chondritic in terms of its major elements (excluding Fe). If the mantle that the veneer was mixed into was not al-ready homogenized, then Os isotopes should not correlate with incompatible elements such as Al. Consequently, either early differentiation of the mantle did not occur or the homogenization of this differentiation must have occurred before the late veneer was added. The timing of the late veneer is itself uncertain but presumably postdated core formation at about 4.45 b.y. and did not postdate the 3.8-3.9 b.y. late bombardment of the Moon. This timing based on siderophile elements is consistent with the Hf isotopic evidence cited above. If the Earth, Moon and Mars had magma oceans, the Earth subsequently rehomogenized whereas the Moon and Mars did not. The simplest solution to this observation is that homogenization of igneous differentiates was never necessary on Earth, either because the hypothetical magma ocean never occurred or because this event did not produce mantle differentiation.

Jones, J. H.↗
NASA NTRS1986

Support for global science - Remote sensing's challenge

Advances in remote sensing techniques are discussed. The benefits possible to remote sensing with the new Earth Observing System, which is composed of the Space Station and coorbiting and polar satellite platforms, are examined. Current changes in the remote sensing field, which involve a change from an industrial society to an informational society, force technology to high technology with high touch, short term to long term, centralized to decentralized, hierarchies to networks, and either/or to multiple option systems are studied. The explanatory and objective types of analyses for investigating biophysical, geochemical, and socioeconomic processes are described; the procedures include: morphometric, cause and effect, temporal and functional and ecological system analyses, inventory, mapping, monitoring, and modeling.

Estes, J. E.↗
NASA NTRS2004

Petrology of New Stannern-trend Eucrites and Eucrite Genesis

Eucrites are basaltic meteorites of the howardite-eucrite-diogenite (HED) suite that originated on a differentiated asteroid, possibly 4 Vesta. Basaltic eucrites are divided into three subgroups based on composition: main group, Stannern-trend, and Nuevo Laredo-trend. The main group and Nuevo Laredo-trend define a sequence formed by fractional crystallization of pigeonite and plagioclase from primitive parent melts [2, 3]. The Stannern-trend cannot be explained this way, but may rather represent a partial-melt sequence of their parent body. However, this model seems inadequate to explain eucrite siderophile element contents, and it is difficult to develop a single unifying model for petrogenesis of all eucrites. Until recently, there were only four Stannern- trend eucrites. One is an anomalous partial cumulate. There is little geochemical variation among these meteorites, so the Stannern-trend was poorly defined. Geochemical studies have identified four additional eucrites as members of the Stannern-trend; one extends the Stannern-trend closer to the main group [5]. No detailed descriptions of these rocks have been published. In order to better integrate these eucrites into the suite, we have done petrologic study of them. They are: LEW 88010, PCA 82501, PCA 91006 and PCA 91179.

Gardner, K. G.↗
NASA NTRS2023

Boninites as Mercury Lava Analogues: Geochemical and Spectral Measurements from Pillow Lavas on Cyprus Island

In the absence of Mercurian rocks or meteorites in our collections, komatiites and boninites are often proposed as the best analog rocks to Mercury lavas. However, despite previous work on the possible analogy between komatiites and Mercury rocks, similar work has not been done for boninites. In this work, we investigate the whole-rock geochemistry and visible/near-infrared (VNIR) spectroscopy of boninitic material collected at three specific areas of the Troodos Massif (Cyprus island). The objective is to evaluate if collected boninites, these along with other boninites present in the literature, can be analogous to Mercury geochemical terranes. On average, we find an unusually high MgO/SiO 2 ratio (0.68) for the boninites from the Troodos Massif compared with previous boninite analysis. This MgO/SiO 2 value is most closely related to the high-Mg regions of Mercury, while the average Al 2 O 3 /SiO 2 ratio (0.25) is consistent with the Mercurian intermediate terrain and to the Mercury’s largest pyroclastic deposit. In addition, further affinity to the high-Mg regions and the intermediate terrains of Mercury are shown in regard to Si vs. Mg, Si vs. Ca, and Si vs. Fe content for one sample in particular. We then conduct magmatic modeling on this specific sample to provide a possible parental melt composition for analogue Mercurian magmas. In conclusion, we suggest these specific locations on the Troodos Massif in Cyprus as good geochemical analogue sites for the high-Mg regions of Mercury and explain how boninites could be important benchmark samples for the chemical and spectral data expected from the BepiColombo mission.

Boninites↗
DOE OSTI2023

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2014

Phase Equilibrium Experiments on Potential Lunar Core Compositions: Extension of Current Knowledge to Multi-Component (Fe-Ni-Si-S-C) Systems

Numerous geophysical and geochemical studies have suggested the existence of a small metallic lunar core, but the composition of that core is not known. Knowledge of the composition can have a large impact on the thermal evolution of the core, its possible early dynamo creation, and its overall size and fraction of solid and liquid. Thermal models predict that the current temperature at the core-mantle boundary of the Moon is near 1650 K. Re-evaluation of Apollo seismic data has highlighted the need for new data in a broader range of bulk core compositions in the PT range of the lunar core. Geochemical measurements have suggested a more volatile-rich Moon than previously thought. And GRAIL mission data may allow much better constraints on the physical nature of the lunar core. All of these factors have led us to determine new phase equilibria experimental studies in the Fe-Ni-S-C-Si system in the relevant PT range of the lunar core that will help constrain the composition of Moon's core.

Righter, K.↗
DOE OSTI2022

Stability of mineral-organic matter associations under varying biogeochemical conditions

It is well-recognized and documented in the literature that a considerable fraction of the soil organic matter (SOM) is associated with Fe-oxide minerals. Recent studies have also showed that the persistence of various classes of SOM may depend on soil Fe mineralogy. Other laboratory studies with model mineral-associated organic matter (MAOM) moieties, which were primarily focused on ferrihydrite, either coated or coprecipitated with natural organic matter (NOM), have suggested that the mineral phase transformation and adsorption-desorption of SOM is rather complex but appears to be largely controlled by the C/Fe ratio of the MAOM moiety and geochemical conditions. It is also clear from multiple recent studies that the nature or type of SOM (microbial vs. plant-derived) adsorbed on mineral surfaces and the rate and extent of SOM biodegradation depend on the type of mineral, and/or metal ions, that are associated with the SOM. Modeling studies that are primarily centered on achieving a better understanding of the effect of mineral type on SOM interactions with minerals and degradation are also reported. These studies have improved the understanding of the nature and stability of MAOM in soils of varying types and under different conditions, positioning researchers to better predict the properties and stabilities of organo-mineral associations in soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

E3SM Water Cycle Visualization Project (Final Report)

Motivated by the need to improve visualizations of Earth's complex water cycle, this team embarked on non-trivial tasks to push from traditional methods of viewing data from simulations as static line graphs and contour plots into new realms with multiple dimensions (three spatial dimensions and the time component) viewable at once. To do this, we chose to feature the extremes in the general circulation of the atmosphere because these Earth system elements are short-lived, but impactful events. We used simulation data produced by the Department of Energy's (DOE) Energy Exascale Earth System Model (E3SM),which is designed to answer DOE science questions using DOE computing facilities. The E3SM project includes three sets of simulation experiments: cryosphere, bio-geochemical cycles, and the water cycle. The water cycle experiment campaign includes a set of simulations designed to understand how the water cycle will change in coming decades under a changing climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Some Opinions on Remote Sensing and Geologic Studies

The principal role of remote sensing data in geologic studies is as a source of geologic information from which meaningful geologic interpretations can be made. Remote sensing data are important in sedimentary basin analysis and other geologic studies as independent and sometimes unique sources of important lithologic and structural information; however, their greatest benefit to exploration-oriented investigations may come when these data are used with other relevant data in a digital database approach to exploration. Modern computer technology facilitates the rapid integration and synthesis of satellite, topographic, gravity, aeromagnetic, geochemical, and other data collected from a given region. Once such data are geometrically registered, they can be digitally processed, within the constraints of defined geologic models, to rapidly identify, and focus further exploration efforts on, target areas that have the greatest potential for success.

Bailey, G. B.↗
NASA NTRSDate not supplied

Mineralogical, Elemental, and Tomographic Reconnaissance Investigation for CLPS (METRIC): A Proposed Mission to the Lunar South Pole-Aitken Basin

The global Lunar Magma Ocean (LMO) hypothesis is based on analyses of ferroan anorthosites (FAN), magnesian plutonic rocks (Mg-suite), and KREEPy (Potassium - K; Rare Earth Elements - REE; Phosphorous - P) material present at the Apollo landing sites and the assumed global distribution of these lithologies. However, orbital spacecraft data over the last two decades have highlighted that the location of the Apollo missions within and near the Procellarum KREEP Terrane (PKT) is compositionally anomalous and not representative of the entire lunar surface and suggests that Apollo samples provide a biased view of the Moon. Further, meteorites, which originate from random areas on the lunar surface, including those not visited by Apollo, Luna, or Chang’E-5, provide complementary datasets that reveal a more complex lunar formation and evolution than the simple LMO onion skin model. While many recent investigations support and refine the LMO model, other recent studies have questioned whether this mechanism alone is responsible for the primordial differentiation of the Moon. Thus, relating the distinct and asymmetrically distributed geochemical terranes on the Moon, especially KREEP which is a predicted global LMO product, to lunar formation and differentiation remains a fundamental goal of lunar science. The South Pole-Aitken (SPA) impact carved a two-and-a-half thousand-kilometer basin into the Moon’s ancient crust, excavating material from as deep as the upper mantle. Orbital remote sensing of the basin indicates compositions with elevated thorium (Th) contents and unique mineralogies that provide a window into early lunar structure and LMO processes. We propose to land a powerful characterization laboratory at Birkeland crater on the far side of the Moon–where Th-rich SPA ejecta have been re-exposed near the surface–in order to address outstanding questions about the evolution of the Moon.

K M Cannon↗
DOE OSTI2022

Emergent Behavior at the Calcite–Water Interface during Reactive Transport in a Simple Microfluidic Channel

Geochemical reactive transport processes in natural mineral-fluid systems may produce a wide array of emergent phenomena that are difficult to predict from basic principles and to reproduce in model systems. In this study, we present experimental results obtained from a simple microfluidic system with which we explored the consequences of reacting the calcite (104) cleavage surface with an acidic Pb-bearing solution (pH = 3.5, [Pb] total = 5 mM) as a function of flow rate. This system is relevant to passive remediation systems for Pb-rich acid mine drainage. We observed periodic banding in the amounts of Pb sorption at flow velocities ≥ 926 μm s -1 , where the band spacing was spatially correlated with the amount of calcite dissolution and the development of micropyramidal topography on the calcite (104) surface. The equivalent coverage of Pb deposited in these Pb-rich bands was at least several monolayers per unit cell, yet there was no evidence for precipitation of any secondary Pb phase implying incorporation of Pb within the near-surface calcite lattice. We also observed spatial variations in nucleation and growth of euhedral secondary Pb-carbonate minerals hydrocerusite and cerussite at flow rates ≤ 278 μm s -1 ,. These findings demonstrate potential for exploiting the rich phenomenology afforded by the interplay among transport phenomena and chemical kinetics in experimental systems designed to yield deeper insights into geochemical self-organization.

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

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗