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 613 records · Page 34

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗

On-Orbit Radiometric Characterization of OLI (Landsat 8) for Applications in Aquatic Remote Sensing

Landsat-8 carries two separate sensors, namely the Operational Land Imager (OLI) and the Thermal Infrared Radiometer Suite (TIRS), that image the earth surface throughout the visible and thermal portions of the spectrum. Compared to Landsat heritage sensors, the OLI has enhanced features, which include its 12-bit radiometric resolution and the addition of a band centered at 443 nm. The dramatically improved data quality/quantity expands existing applications of Landsat imagery in aquatic sciences from the retrieval of bio-geochemical properties, such as near-surface concentrations of chlorophyll-a (CHL) and total suspended solids (TSS), to benthic mapping. This study offers analysis of OLI's absolute radiometric performance over bodies ofwater using benchmark observations, namely the top-of-atmosphere (TOA) ocean color observations and marine in situ radiometric measurements. Sensor-to-sensor comparisons are performed to derive gain factors (g1) from near-concurrent observations in TOA radiance and reflectance domains. The gains in the radiance domain were further validated/adjusted by determining a second set of gains (g2) via analysis of OLI-derived water-leaving radiance, i.e., L(sub w) (gamma), against in situ measurements made either at the Ocean Color AErosol RObotic NETwork (AERONETOC) sites or during field campaigns. The analyses yield the OLI calibration uncertainties that need to be accounted forwhen studying aquatic environments. Itwas found that, for the visible and near-infrared channels, the OLI radiometric responses, on average, are well in agreement (b 2 % discrepancies) with the TOA radiances estimated by ocean color satellites or those predicted by models based onmeasurements of aquatic and atmospheric properties. However, the TOA radiance at the new 443-nm band is found to be, on average, 3.4 % larger than the reference observations. The inter-sensor comparisons in the reflectance domain, however, indicated slightly different results with the OLI responses being low in the blue bands. To enhance the retrieval accuracy of aquatic-science products from OLI datasets, sets of temporally averaged gains (radiance and reflectance) are derived and recommended for use prior to the retrieval of in-water products.

Remote Sensing↗

Big Data and AI at DoE's Legacy Sites - 20546

More than 30 years have passed since DOE started the decommissioning of nuclear weapon complexes and the clean-up of soil and groundwater. All the sites have been collecting and archiving soil and groundwater monitoring datasets; particularly contaminant concentration time-series. These datasets provide unparalleled opportunities to understand the system behavior (including more fundamental hydrological and geochemical processes, the response to various perturbations, the long-term trend and environmental decay rate towards the regulatory limit). This understanding is critical for providing multiple lines of evidences that can support site closure. In this study, we explore the machine learning (ML) and artificial intelligence (AI) applications to the long-term soil and groundwater management at DoE's legacy sites. ML can improve our understanding of the subsurface systems, which is critical for long-term monitoring and management of the sites, while AI can automate or support some of decision-making processes (e.g., anomaly detection, monitoring well placements). The particular focuses are to develop general algorithms to: (1) to identify distinct spatiotemporal patterns and to identify several groups that have similar temporal behaviors, using unsupervised clustering methods, (2) identify the different temporal scales of hydrological responses to climate perturbations by time-series analysis, and (3) reduce the number of monitoring wells by identifying the minimum sufficient number of wells to capture the heterogeneity of the groundwater contaminant plume and concentration distribution, using the Gaussian Process model. We demonstrate our methodology at the Savannah River Site F-Area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗

Processes involved in the formation of magnesian-suite plutonic rocks from the highlands of the Earth's Moon

The earliest evolution of the Moon likely included the formation of a magma ocean and the subsequent development of anorthositic flotation cumulates. This primary anorthositic crust was then intruded by mafic magmas which crystallized to form the lunar highlands magnesian suite. The present study is a compilation of petrologic, mineral-chemical, and geochemical information on all pristine magnesian-suite plutonic rocks and the interpretation of this data in light of 18 'new' samples. Of these 18 clasts taken from Apollo 14 breccias, 12 are probably pristine and include four dunites, two norites, four troctolites, and two anorthosites. Radiogenic isotopic whole rock data also are reported for one of the 'probably pristine' anorthositic troctolites, sample 14303,347. The relatively low Rb content and high Sm and Nd abundances of 14303,347 suggest that this cumulate rock was derived from a parental magma which had these chemical characteristics. Trace element, isotopic, and mineral-chemical data are used to interpret the total highlands magnesian suite as crustal precipitates of a primitive KREEP (possessing a K-, rare earth element (REE)-, and P-enriched chemical signature) basalt magma. This KREEP basalt was created by the mixing of ascending ultramafic melts from the lunar interior with urKREEP (the late, K-, REE-, and P-enriched residuum of the lunar magma ocean). A few samples of the magnesian suite with extremely elevated large-ion lithophile elements (5-10x other magnesian-suite rocks) cannot be explained by this model or any other model of autometasomatism, equilibrium crystallization, or 'local melt-pocket equilibrium' without recourse to an extremely large-ion lithophile element-enriched parent liquid. It is difficult to generate parental liquids which are 2-4 x higher in the REE than average lunar KREEP, unless the liquids are the basic complement of a liquid-liquid pair, i.e., the so-called 'REEP-fraction,' from the silicate liquid immiscibility of urKREEP. Scarce age information on lunar rocks suggests that magnesian-suite magmatism was initiated at progressively more recent time from the northeast to the southwest on the lunar nearside from 4.45 to 4.25 Ga.

Snyder, Gregory A.↗

Multi-Anvil Experimentation Applied to Planetary Differentiation

Planets undergo differentiation that includes segregation of metal from silicate at high temperatures and pressures ranging from deep planetary core pressures (>300 GPa for Earth)to very shallow conditions of asteroids (<100 MPa).The multi-anvil solid media apparatus accesses the middle part of this range from 3 to 30 GPa –pressure relevant to the interior of Mercury, Venus, Earth, Earth’s Moon, and Mars. Early planets are thought to have experienced high temperatures from a combination of heat sources including radioactive decay, gravitational and accretional heating, and impact processes. These heating events led to melting of mantles and cores, thus requiring an understanding of solid-liquid equilibria in metal-silicate systems. In 2006 we established a multi-anvil facility at NASA-JSC combining an 880 ton press and a Kawai/Walker type module from Rockland Research. Our high PT work has been greatly facilitated by use of the COMPRES multi-anvil assemblies (1). Our recent work has included studies of element partitioning between liquid metal and liquid silicate(e.g., 2), as well as between minerals and melts(e.g., 3), both of which have led to better constraints on the timing and conditions of planetary differentiation(e.g. 4). Several examples involving sustained efforts will be summarized below and for the presentation. The distribution of siderophile (iron-loving) elements between core and mantle is controlled by metal-silicate equilibrium across a wide range of pressures. Therefore, experimentation across this pressure range helps to calibrate elemental partitioning models that can be applied to planets, and used to predict mantle chemistry and composition during planetary differentiation. Our studies have focused on a wide range of siderophile elements (refractory Ni, Co, W, Mo; volatile P, Ga, Cu, Sn, Sb; highly siderophile Au, Pd) that have constrained partitioning, valence, and isotopic fractionation, and applied to Earth, Moon, and Mars. When molten mantles (magma oceans) cool enough to initiate crystallization, the solids precipitate at depth and in large planets this involves high pressure phases like garnet, majorite, akimotoite, and ringwoodite. As these solids precipitate they can segregate from liquid by density contrasts, thus causing elemental fractionation which can be used to decipher timing of differentiation. Mineral/melt and metal/silicate equilibria in our lab have helped to better understand high pressure fractionation of isotopic parent/daughter pairs Hf/W, Mn/Cr, Pd/Ag, Pt/Os, Re/Os, and U/Pb, and their application to Earth, Moon and Mars. There remains great potential for multi-anvil experimentation to shed light on many pressure–dependent aspects of planetary evolution such as core formation, high pressure phase equilibria, redox equilibria, and volatile evolution and storage. References1. Leinenweber, K., et al.(2012) American Mineralogist,97, 353–368. 2.Righter, K., et al.(2020) Geochem. Persp. Lett.15, 1-6.3. Righter, K., et al.(2020) Met. Planet. Sci 55, 2741-2757.4. Righter, K., et al.(2020)Earth and Planetary Science Letters,552, 116590.

pressure↗

From Lake to River: Documenting an Environmental Transition Across the Jura/Knockfarril Hill Members Boundary in the Glen Torridon Region of Gale Crater (Mars)

Between January 2019 and January 2021, the Mars Science Laboratory team explored the Glen Torridon (GT) region in Gale crater (Mars), known for its orbital detection of clay minerals. Mastcam, Mars Hand Lens Imager, and ChemCam data are used in an integrated sedimentological and geochemical study to characterize the Jura member of the upper Murray formation and the Knockfarril Hill member of the overlying Carolyn Shoemaker formation in northern GT. The studied strata show a progressive transition represented by interfingering beds of fine–grained, recessive mudstones of the Jura member and coarser–grained, cross–stratified sandstones attributed to the Knockfarril Hill member. Whereas the former are interpreted as lacustrine deposits, the latter are interpreted as predominantly fluvial deposits. The geochemical composition seen by the ChemCam instrument show K 2 O–rich mudstones (~1–2 wt.%) versus MgO–rich sandstones (>6 wt.%), relative to the average composition of the underlying Murray formation. We document consistent sedimentary and geochemical data sets showing that low–energy mudstones of the Jura member are associated with the K–rich endmember, and that high–energy cross–stratified sandstones of the Knockfarril Hill member are associated with the Mg–rich endmember, regardless of stratigraphic position. The Jura to Knockfarril Hill transition therefore marks a significant paleoenvironmental change, where a long–lived and comparatively quiescent lacustrine setting progressively changes into a more energetic fluvial setting, as a consequence of shoreline regression due to either increased sediment supply or lake–level drop.

58 GEOSCIENCES↗

Surface Reactivity Analysis of the Crude Oil–Brine–Limestone Interface for a Comprehensive Understanding of the Low-Salinity Waterflooding Mechanism

Low-salinity waterflooding (LSWF) has proven to improve oil recovery in carbonate formations through rock wettability alteration, although the underlying mechanism remains elusive. Multivalent ionic exchange and calcite dissolution have usually been investigated using geochemical analysis in secondary coreflooding. In this work, coreflooding, in tertiary mode, coupled with a surface reactivity analysis approach was employed to investigate the interplay of wettability alteration mechanisms such as mineral dissolution, electrostatic bond attraction, and the effect of pH at in situ conditions. Improved oil recovery (IOR) in tertiary mode observed by coreflooding in Indiana limestone rocks showed an ionic strength dependence, that is, reducing brine ionic strength resulted in an increase in oil recovery. Coreflooding results showed that the seawater and low-salinity brines deprived of Mg 2+ ions resulted in the lowest IOR in tertiary mode, indicating the significance of Mg 2+ on IOR in limestone rocks. Similar results were observed through the contact angle measurement showing the limestone rock wettability state dependence on ionic strength and the effect of Mg2+ ions. Surface reactivity analysis showed an increase in solution pH, Ca 2+ and Mg 2+ ions concentration in the effluent solution from the coreflooding in tertiary mode using low salinity brines (about 40 and 20% increase in the effluent composition for Ca 2+ and Mg 2+ , respectively). These changes in solution composition were used to calculate the in situ oil–brine and rock–brine zeta potential using a validated surface complexation model, showing the changes of zeta potential as brine is injected into limestone rocks. The results show that using seawater-like brine in tertiary mode resulted in no mineral dissolution or ionic exchange. However, improved oil recovery (IOR) using such seawater-like brine was due to wettability alteration caused by reduced electrostatic bond attraction associated with Mg 2+ ions [from 2.6 × 10 –13 (mol/m 2 ) 2 for formation water salinity to 1.5 × 10 –13 (mol/m 2 ) 2 for seawater salinity]. Using low-salinity brines in tertiary mode improved oil recovery by mineral dissolution, resulting in oil desorption and an increase in solution pH. Finally, the increase in solution pH also resulted in reduced electrostatic bond attraction which lead to rock wettability alteration using low-salinity brines.

04 OIL SHALES AND TAR SANDS↗

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mineralogical, Elemental, and Tomographic Reconnaissance Investigation for CLPS (METRIC)

METRIC is a robotic science laboratory that can determine the mineralogy, elemental chemistry, micromorphology, and thermophysical properties of planetary regolith. The METRIC suite comprises METRIC XRD/F, an X-ray diffraction/X-ray fluorescence instrument that can determine the mineralogy and elemental chemistry of regolith samples; METRIC XCT, a micro X-ray computed tomography instrument that can be used to evaluate grain/crystallite sizes and textures; METRIC IRS, an imaging spectrometer mounted on a rover that can determine mineralogy and thermophysical properties at the landing site; and a pneumatic sample collection, processing, distribution system developed by Honeybee Robotics. The payload elements could be deployed on a static lander or a rover. Data returned from the METRIC payload would inform origin, formation, and evolution of rocky planetary bodies. METRIC XRD/F draws on heritage from the CheMin instrument on the Mars Science Laboratory (MSL) Curiosity rover [1], with a few important improvements. Like CheMin, METRIC XRD/F operates in transmission geometry and uses piezoelectric actuators on sample cells in a tuning fork geometry to induce convective grain motion of the regolith to create a randomly oriented powder. MSL CheMin uses an energy-sensitive CCD to collect XRD patterns and XRF spectra simultaneously from the same sample cell, resulting in qualitative XRF data. METRIC XRD/F uses two different sample cells, one optimized for XRD and one optimized for XRF, and a silicon drift detector to detect fluoresced X-rays. This improvement to the XRF capabilities provides quantitative geochemical data of major elements down to Z = 11 and allows for the detection of minor and trace elements that are critical for evaluating geologic evolution of the Moon (e.g., P and Th). Modest improvements to the XRD geometry and hardware allow for better angular resolution and the ability to distinguish between members of the pyroxene group. METRIC XCT uses the same geometry and much of the same hardware as METRIC XRD/F, where a CCD would capture images of a regolith sample in a 3 mm diameter sample tube that is rotated 360° in steps <1°. Image brightness can be used to infer compositional data, where brighter materials indicate a higher Z, much like scanning electron microscopy. Data from METRIC XCT complement those from METRIC XRD/F. Particle size, shape, and texture can provide petrologic and provenance information, whereas vesicle size and morphology in volcanic or impact melt lithologies can inform cooling rates. METRIC IRS is a hyperspectral thermal imager that can be mounted to a lander or rover to provide mineralogical data from the broader landing site and help determine whether the samples analyzed by METRIC XRD/F and XCT are representative. The METRIC IRS spectral range (8–14 μm) and resolution (10.8 cm-1) allow for quantitative mineralogy from modelling Reststrahlen bands of major rock-forming minerals (e.g., silicates, phosphates). Radiance cubes can be processed and modelled with an onboard high-performance computer to determine mineral abundances of plagioclase, high-Ca pyroxene, pigeonite, orthopyroxene, olivine, and glass. Regolith samples can be acquired, processed, and delivered to the X-ray instruments via multiple sample handling systems, but the pneumatic sampling systems developed by Honeybee Robotics [e.g., 2] are best suited for relatively low-cost missions that are being competed for the Moon (e.g., NASA’s Payloads and Research Investigations for the Surface of the Moon program). There are pneumatic sampling systems that collect surface material and other systems that pneumatically drill up to ~1 m below the surface, providing material that has not been space weathered and has not been affected by the lander’s exhaust. [1] Blake, D. F., Vaniman, D., Achilles, C., Anderson, R., Bish, D., et al. (2012). Space Sci. Rev. 170, 341-478. https://doi.org/10.1007/s11214-012-9905-1. [2] Zacny, K., Betts, B., Hedlund, M., Long, P., Gramlich, M., Tura, K., Chu, P., Jacob, A., Garcia, A. (2014). IEEE Aerospace Conference, 3-7 March 2014, Big Sky, MT, U.S.A.

X-ray diffraction↗

The combined effects of Mg 2+ and Sr 2+ incorporation during CaCO 3 precipitation and crystal growth

Calcite (CaCO 3 ) composition and properties are defined by the chemical environment in which CaCO 3 forms. However, a complete understanding of the relationship between aqueous chemistry during calcite precipitation and resulting chemical and physical CaCO 3 properties remains elusive; therefore, we present an investigation into the coupled effects of divalent cations Sr 2+ and Mg 2+ on CaCO 3 precipitation and subsequent crystal growth. Through chemical analysis of the aqueous phases and microscopy of the resulting calcite phases in compliment with density functional theory calculations, we elucidate the relationship between crystal growth and the resulting composition (elemental and isotopic) of calcite. The results of this experimental and modeling work suggest that Mg 2+ and Sr 2+ have cation-specific impacts that inhibit calcite crystal growth, including: (1) Sr 2+ incorporates more readily into calcite than Mg 2+ (D Sr > D Mg ), and increasing [Sr 2+ ] t or [Mg 2+ ] t increases D Sr ; (2) the inclusion of Mg 2+ into structure leads to a reduction in the calcite unit cell volume, whereas Sr 2+ leads to an expansion; (3) the inclusion of both Mg 2+ and Sr 2+ results in a distribution of unit cell impacts based on the relative positions of the Sr 2+ and Mg 2+ in the lattice. These experiments were conducted at saturation indices of CaCO 3 of ~4.1, favoring rapid precipitation. This rapid precipitation resulted in observed Sr isotope fractionation confirming Sr isotopic fractionation is dependent upon the precipitation rate. We further note that the precipitation and growth of calcite favors the incorporation of the lighter 86 Sr isotope over the heavier 87 Sr isotope, regardless of the initial solution conditions, and the degree of fractionation increases with D Sr . In sum, these results demonstrate the impact of solution environment to influence the incorporation behavior and crystal growth behavior of calcite. These factors are important to understand in order to effectively use geochemical signatures resulting from calcite precipitation or dissolution to gain specific information.

58 GEOSCIENCES↗

Formation of ferrihydrite and associated iron corrosion products in permeable reactive barriers of zero-valent iron

Ferrihydrite, which is known to form in the presence of oxygen and to be stabilized by the adsorption of Si, PO4 and SO4, is ubiquitous in the fine-grained fractions of permeable reactive barrier (PRB) samples from the U.S. Coast Guard Support Center (Elizabeth City, NC) and the Denver Federal Center (Lakewood, CO) studied by high-resolution transmission electron microscopy and selected area electron diffraction. The concurrent energy-dispersive X-ray data indicate a strong association between ferrihydrite and metals such as Si, Ca, and Cr. Magnetite, green rust 1, aragonite, calcite, mackinawite, greigite and lepidocrocite were also present, indicative of a geochemical environment that is temporally and spatially heterogeneous. Whereas magnetite, which is known to form due to anaerobic Fe0 corrosion, passivates the Fe0 surface, ferrihydrite precipitation occurs away from the immediate Fe0 surface, forming small (<0.1 microm) discrete clusters. Consequently, Fe0-PRBs may remain effective for a longer period of time in slightly oxidized groundwater systems where ferrihydrite formation occurs compared to oxygen-depleted systems where magnetite passivation occurs. The ubiquitous presence of ferrihydrite suggests that the use of Fe0-PRBs may be extended to applications that require contaminant adsorption rather than, or in addition to, redox-promoted contaminant degradation.

Metals/chemistry↗

Prebiotic Evolution of Nitrogen Compounds

Support from this four year grant has funded our research on two general problems. One involves attempts to model the abiotic formation of simple source compounds for functional biomolecules, their concentration from dilute state in the hydrosphere and, in several cases, surface induced reactions to form precursor monomers for bioactive end products (refs. 1-5). Because of the pervasiveness and antiquity of phosphate based biochemistry and the catalytic activity of RNA we have exploring the hypothesis of an RNA World as an early stage in the emergence of life. This concept is now rather generally considered, but has been questioned due to the earlier lack of an experimentally demonstrated successful scheme for the spontaneous formation of ribose phosphate, the key backbone molecule in RNA. That impediment has now been removed. This has been achieved by demonstrating probable sources of activated (condensed) highly soluble and strongly sorbed phosphates in nature (Refs. 1,2) and effective condensation of aldehyde phosphates to form ribose phosphate in high yield (ref.6), thereby placing the RNA World concept on a somewhat safer experimental footing. Like all work in this field these experiments are oversimplifications that largely ignore competing side reactions with other compounds expected to be present. None the less our choice of experimental conditions aim at selective processes that eliminate interfering reactions. We have also sought to narrow the credibility gap by simulating geophysically and geochemically plausible conditions surrounding the putative prebiotic reactions.

Arrhenius, G.↗

BONCAT-Live for isolation and cultivation of active environmental bacteria

In diverse environments, microbes drive a myriad of processes, from geochemical and nutrient cycling to interspecies interactions, including associations with plants and animals. Their physiological state is dynamic and impacted by abiotic and biotic conditions, responding to environmental fluctuations by changes in cellular metabolism, according to their genetic potential. Molecular, cellular, and genomic approaches can identify and measure microbial responses and adaptation to environmental changes in native communities. However, isolating individual microbial cells that respond to specific changes for cultivation has been difficult. To address this, we implemented a novel bacterial isolation approach (BONCAT-Live) by integrating bio-orthogonal non-canonical amino acid tagging (BONCAT) in diverse native communities, with isolation and cultivation of cells responding to specific stimuli, at different time scales. In frozen Arctic permafrost samples, we identified and isolated dormant bacteria that become active after thawing under native or nutrient-enriched conditions. From the Populus tree rhizosphere, we isolated strains that thrive under high concentrations of root exudates that act as defense compounds and nutrients. In the human microbiome, we identified and isolated bacteria that rapidly proliferated when exposed to metabolites provided by the host or other co-occurring microbes. Further characterization of isolated bacterial strains will provide opportunities for in-depth determination of how these microbes adapt to changes in their environments, individually and as part of model communities.

Analytical Methods↗

Convection in three dimensions with surface plates - Generation of toroidal flow

This work presents numerical calculations of mantle convection that incorporate some of the basic observational constraints imposed by plate tectonics. The model is three-dimensional and includes surface plates; it allows plate velocity to change dynamically according to the forces which result from convection. It is shown that plates are an effective means of introducing a toroidal component into the flow field. After initial transients the plate motion is nearly parallel to transform faults and in the direction that tends to minimize the toroidal flow field. The toroidal field decays with depth from its value at the surface; the poloidal field is relatively constant throughout the layer but falls off slightly at the top and bottom boundaries. Layered viscosity increasing with depth causes the toroidal field to decay more rapidly, effectively confining it to the upper, low-viscosity layer. The effect of viscosity layering on the poloidal field is relatively small, which is attributed to its generation by temperature variations distributed throughout the system. The generation of toroidal flow by surface plates would seem to account for the observed nearly equal energy of toroidal and poloidal fields of plate motions on the earth. A low-viscosity region in the upper mantle will cause the toroidal flow to decay significantly before reaching the lower mantle. The resulting concentration of toroidal flow in the upper mantle may result in more thorough mixing there and account for some of the geochemical and isotopic differences proposed to exist between the upper and lower mantles.

Gable, Carl W.↗

Exploring Geothermal Potential of Great Basin Sub-Regions: Preprint

The INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems (INGENIOUS) project aims to discover new, economically viable hidden geothermal systems in the Great Basin region by building on previous work in play fairway analysis and machine learning. A key objective of this project is to develop an exploration workflow to reduce geothermal exploration risks for hidden geothermal systems. A single preliminary play fairway workflow was developed from the assessment of the regional INGENIOUS geological, geophysical, and geochemical datasets. This workflow provided new preliminary predictive geothermal fairway maps for the INGENIOUS study area, which encompasses most of Nevada, western Utah, southern Idaho, southeastern Oregon, and easternmost California. However, a recent study (incorporating machine learning techniques) of a portion of Nevada identified four geologic domains and determined that the relative importance of individual datasets or features as indicators of geothermal potential may differ across these domains. The INGENIOUS study area includes a much larger and more geologically diverse region; therefore, additional geologic domains or sub-regions are expected. To assess the sub-regions in the INGENIOUS study area, principal component analysis and k-means clustering were applied. Preliminary results indicate that the INGENIOUS regional data cluster into groups that relate to different geologic domains in the Great Basin region. These include domains such as the Walker Lane, extensional western Great Basin region, broad lower strain region in the eastern Great Basin of western Utah and eastern Nevada, Quaternary volcanic fields, and the area adjacent to the Snake River Plain. These clusters are assessed to determine the key geologic drivers of the identified clusters. Understanding this variability can provide key insights for the exploration and characterization of hidden geothermal systems in the Great Basin region and could indicate the need to develop multiple geothermal conceptual models and play fairway workflows for the INGENIOUS study area.

GEOTHERMAL ENERGY↗

Exploring Geothermal Potential of Great Basin Sub-Regions

The INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems (INGENIOUS) project aims to discover new, economically viable hidden geothermal systems in the Great Basin region by building on previous work in play fairway analysis and machine learning. A key objective of this project is to develop an exploration workflow to reduce geothermal exploration risks for hidden geothermal systems. A single preliminary play fairway workflow was developed from the assessment of the regional INGENIOUS geological, geophysical, and geochemical datasets. This workflow provided new preliminary predictive geothermal fairway maps for the INGENIOUS study area, which encompasses most of Nevada, western Utah, southern Idaho, southeastern Oregon, and easternmost California. However, a recent study (incorporating machine learning techniques) of a portion of Nevada identified four geologic domains and determined that the relative importance of individual datasets or features as indicators of geothermal potential may differ across these domains. The INGENIOUS study area includes a much larger and more geologically diverse region; therefore, additional geologic domains or sub-regions are expected. To assess the sub-regions in the INGENIOUS study area, principal component analysis and k-means clustering were applied. Preliminary results indicate that the INGENIOUS regional data cluster into groups that relate to different geologic domains in the Great Basin region. These include domains such as the Walker Lane, extensional western Great Basin region, broad lower strain region in the eastern Great Basin of western Utah and eastern Nevada, Quaternary volcanic fields, and the area adjacent to the Snake River Plain. These clusters are assessed to determine the key geologic drivers of the identified clusters. Understanding this variability can provide key insights for the exploration and characterization of hidden geothermal systems in the Great Basin region and could indicate the need to develop multiple geothermal conceptual models and play fairway workflows for the INGENIOUS study area.

exploration↗