SEARCH · Search NASA

Results for “EARTH CRUST”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $\$3.2$ per tonne of coal and $\$186$ billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$\$4.4$ per tonne of ash and $\$1.2$ billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2022

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗

DOE OSTI2020

Comprehensive geoneutrino analysis with Borexino

This paper presents a comprehensive geoneutrino measurement using the Borexino detector, located at Laboratori Nazionali del Gran Sasso (LNGS) in Italy. The analysis is the result of 3262.74 days of data between December 2007 and April 2019. The paper describes improved analysis techniques and optimized data selection, which includes enlarged fiducial volume and sophisticated cosmogenic veto. The reported exposure of ( 1.29 ± 0.05 ) × 10 32 protons × year represents an increase by a factor of two over a previous Borexino analysis reported in 2015. By observing 52.6 - 8.6 + 9.4 ( stat ) - 2.1 + 2.7 ( sys ) geoneutrinos (68% interval) from U 238 and Th 232 , a geoneutrino signal of 47.0 - 7.7 + 8.4 ( stat ) - 1.9 + 2.4 ( sys ) TNU with - 17.2 + 18.3 % total precision was obtained. This result assumes the same Th/U mass ratio as found in chondritic CI meteorites but compatible results were found when contributions from U 238 and Th 232 were both fit as free parameters. Antineutrino background from reactors is fit unconstrained and found compatible with the expectations. The null-hypothesis of observing a geoneutrino signal from the mantle is excluded at a 99.0% C.L. when exploiting detailed knowledge of the local crust near the experimental site. Measured mantle signal of 21.2 - 9.0 + 9.5 ( stat ) - 0.9 + 1.1 ( sys ) TNU corresponds to the production of a radiogenic heat of 24.6 - 10.4 + 11.1 TW (68% interval) from U 238 and Th 232 in the mantle. Assuming 18% contribution of K 40 in the mantle and 8.1 - 1.4 + 1.9 TW of total radiogenic heat of the lithosphere, the Borexino estimate of the total radiogenic heat of the Earth is 38.2 - 12.7 + 13.6 TW , which corresponds to the convective Urey ratio of 0.78 - 0.28 + 0.41 . These values are compatible with different geological predictions, however there is a ~ 2.4 σ tension with those Earth models which predict the lowest concentration of heat-producing elements in the mantle. In addition, by constraining the number of expected reactor antineutrino events, the existence of a hypothetical georeactor at the center of the Earth having power greater than 2.4 TW is excluded at 95% C.L. Particular attention is given to the description of all analysis details which should be of interest for the next generation of geoneutrino measurements using liquid scintillator detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗

DOE OSTI2022

Hydrous SiO 2 in subducted oceanic crust and H 2 O transport to the core-mantle boundary

Subduction of oceanic lithosphere transports surface H 2 O into the mantle. Recent studies show that dense SiO 2 in the form of stishovite, an abundant mineral in subducted oceanic crust at depths greater than ~270 km, has the potential to host and transport a considerable amount of H 2 O into the lower mantle, but the H 2 O storage capacity of SiO 2 phases at high pressure and temperature remains uncertain. We investigate the hydration of stishovite and its higher-pressure polymorphs, β-stishovite and seifertite, with in situ X-ray diffraction experiments at high pressures and temperatures. The H 2 O contents in SiO 2 phases are quantified based on observed increases in unit cell volume relative to the anhydrous SiO 2 system. Density functional theory (DFT) computations permit calibration of water content as a function of volume change based on interstitial substitution of H 2 O. Regression of our experimental data indicates an H 2 O storage capacity in stishovite of ~3.5 wt% in the transition zone and shallow lower mantle, decreasing to about 0.8 wt% at the base of the mantle. We find that SiO 2 -bearing subducted oceanic crust can accommodate all the H 2 O in slab lithosphere that survives sub-arc dehydration. Hydration of silica phases in subducted oceanic crust and their unparalleled capacity to host significant amounts of H 2 O even at high mantle temperatures provides a unique mechanism for transport and storage of water in the deepest mantle.

58 GEOSCIENCES↗

DOE OSTI2024

Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)

Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.

58 GEOSCIENCES↗

DOE OSTI2024

Xenon compatibility in magmatic processes: Hadean to current contexts

Xenon (Xe) behaviour in petrological processes, albeit essential to constrain mantle ingassing and degassing models, is elusive due to its volatile nature, and lack of direct investigation at the pressures (P) and temperatures (T) relevant to magma formation and crystallisation at depth. Xenon stands out amongst noble gases due to its unique reactivity with silicates of the lower crust and upper mantle, which could at least partially explain that published mineral/melt partitioning coefficients span up to six orders of magnitude. We report partition coefficients of Xe using in situ X-ray fluorescence at high P and T, and mass spectrometry analyses. Xenon is found to be moderately incompatible in anorthite-clinopyroxene mix in equilibrium with basalt (partition coefficient value of 0.16 ± 0.06), and compatible in olivine in equilibrium with basalt (partition coefficients in the range 88 ± 22 to 302 ± 46). While Xe is, thus, concentrated in basaltic melts coexisting with crystallising pyroxenes and feldspars, it is strongly retained in olivine at depth. Consequently, Xe originally contained in solid Earth has been preferentially retained at depth throughout Earth’s history, from the magma ocean stages to present day partial mantle melting processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Warming decreases desert ecosystem functioning by altering biocrusts in drylands

Warming and precipitation fluctuations are changing desert ecosystems in global drylands. However, the effects of climate change on keystone species such as cryptogamic biocrust in drylands remain relatively under-investigated, even though biocrusts play a vital role in desert ecosystems. Here we conducted a long-term experiment (14 years) to simulate the responses of two main types of biocrusts to warming coupled with reduction in precipitation that was achieved by open-top chambers (OTCs) to simulate the predicted warming and precipitation decreasing under climate change scenario. We also conduct a structural analysis to evaluate the resulting changes in desert ecosystem functioning (carbon and nitrogen cycling). Neither warming and corresponding rainfall reduction treatments had a negative effect on lichen species richness, but both treatments reduced lichen cover and biomass. The negative effects of warming on moss-dominated crusts were much greater than those on lichen-dominated crusts. Although mosses and lichens had varying degree responses to warming, the loss of mosses and decreased lichen cover and biomass, as well as the shortening of the wet time, resulted in a reduction in carbon and nitrogen fixation, soil enzyme activity and water-holding capacity of biocrusts and topsoil. These impacts collectively change the water balance of drylands and weaken the hydrological and biogeochemical function of biocrusts. Synthesis and applications: Results from this long-term experiment suggest that the ecosystem C and N cycling and water balance of global drylands may be highly impacted by climate change, in part because of the response of biocrusts, which contribute an important implication for both dryland restoration and earth system dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Dryland feedbacks to future climate change: how species mortality and replacement will affect coupled biogeochemical cycles and energy balance (Final Technical Report)

The overall goal of the dryland climate manipulation experiment is to improve our understanding and ability to forecast dryland responses to increasing temperatures and altered precipitation regimes. Drylands cover over 40% of Earth’s terrestrial surface and are predicted to increase in size 11-23% by 2100. Yet our understanding of this vast biome, which exchanges enormous amounts of CO 2 and energy with the atmosphere, is exceedingly poor. This project aims to transform and improve that understanding in order to reduce uncertainty and increase the confidence with which we can make global predictions of future climate. We focus on both above- and belowground processes, and explore temperature controls over critical aspects of carbon and nutrient cycling for dryland plants, soil, and microbes. We have 60 climate manipulation plots where plants (from above the canopy) and soils are warmed using infrared lamps and where precipitation has been altered using targeted watering regimes. We have measured myriad aspects of plant physiology, morphology, carbon exchange, energy exchange, and phenology, as well as analogous climate-induced changes to biological soil crusts and soil microbial communities. Our goals are to investigate temperature responses and acclimation potential of dryland plant and soil processes. Concurrent soil incubation experiments complement the in-situ soil warming component and enable more controlled mechanistic investigations of temperature response on microbial function. Automated soil CO 2 efflux chambers within the plots provide phenomenal opportunity to improve our understanding of dryland CO 2 exchange with the atmosphere and are coupled with plant CO 2 exchange measurements. The project has been exceptionally successful, with 35 peer-reviewed publications (and more in review), including in high-profile journals such as Nature Climate Change and PNAS. We have given over fifty presentations on the work and it has been featured in many high-profile media outlets, such as Discover Magazine and on NPR.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗

DOE OSTI2023

Molybdenum in basalt-hosted seafloor hydrothermal systems: Experimental, theoretical, and field sampling approaches

Here, seafloor hydrothermal vents represent potential sources of Mo and other biologically relevant transition metals to the global ocean, complementing continental runoff. Here, we use a combination of experimental, theoretical, and field-sampling approaches to investigate the behavior of Mo in basalt-hosted seafloor hydrothermal systems to provide insight into the processes controlling Mo concentrations in hydrothermal fluids and to derive estimates of vent fluid Mo concentrations and fluxes. Results of this study demonstrate that reaction fluids generated from 350°C, 500 bar hydrothermal basalt alteration experiments contain 775–801 nmol/kg Mo and are thus comparable to a recently collected time series of natural seafloor vent fluids that contained 200–220 nmol/kg Mo at 302°C and 29–30 nmol/kg at 281–282°C (Evans et al., 2023). Synchrotron-based analyses of experimentally altered basalt produced in this study and additional natural samples of altered oceanic crust originally collected from Pito Deep Rift reveal the presence of Mo-rich particles consistent with trace molybdenite. Comparisons of Mo:Cu ratios in natural vent fluids and near-vent sediment trap samples from Main Endeavour Field indicate that vent fluid Mo is readily incorporated into buoyant plume particles and advected out of the near-vent field, analogous to previous mass balance studies of Cu in this region. Thermodynamic calculations of molybdenite solubility in the context of mineral-buffered hydrothermal fluids and comparisons with natural and experimental hydrothermal fluids suggest that high-temperature vent fluids contain 30–1500 nmol/kg Mo. While a minor component of the modern Mo budget, hydrothermal Mo fluxes are estimated to have constituted 0.3–200× the contemporaneous continental weathering fluxes prior to the ~2.4 Ga ago “Great Oxidation Event” and widespread oxidative continental weathering. Overall, identification of hydrothermal vents as a source of Mo-rich plume particles with potential for dispersal into the wider marine environment has significant implications for hypotheses regarding the co-evolution of Life and Earth’s environments, specifically the form and availability of Mo in anoxic Archean-Eon oceans, where Mo-dependent enzymatic pathways are thought to have emerged and subsequently evolved.

58 GEOSCIENCES↗

DOE OSTI2021

Surface-wave tomography of the Emeishan large igneous province (China): Magma storage system, hidden hotspot track, and its impact on the Capitanian mass extinction

Large igneous provinces (LIPs) are commonly associated with mass extinctions. However, the precise relations between LIPs and their impacts on biodiversity is enigmatic, given that they can be asynchronous. It has been proposed that the environmental impacts are primarily related to sill emplacement. Therefore, the structure of LIPs' magma storage system is critical because it dictates the occurrence and timing of mass extinction. We use surface-wave tomography to image the lithosphere under the Permian Emeishan large igneous province (ELIP) in southwestern China. We find a northeast-trending zone of high shear-wave velocity (Vs) and negative radial anisotropy (Vsv > Vsh; v and h are vertically and horizontally polarized S waves, respectively) in the crust and lithosphere. We rule out the possibilities of rifting or orogenesis to explain these seismic characteristics and interpret the seismic anomaly as a mafic-ultramafic, dike-dominated magma storage system of the ELIP. We further propose that the anomaly represents a hidden hotspot track that was emplaced before the ELIP eruption. A zone of higher velocity but less-negative radial anisotropy, on the hotspot track but to the northeast of the eruption center in the Panxi region, reflects an elevated proportion of sills emplaced at the incipient stage of the ELIP. Liberation of poisonous gases by the early sill intrusions explains why the mid-Capitanian global biota crisis preceded the peak ELIP eruption by 2–3 m.y.

58 GEOSCIENCES↗

DOE OSTI2023

Equation of State and Spin Crossover of (Al, Fe)–Phase H

The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases δ-AlOOH (delta phase), MgSiO 2 (OH) 2 (phase H), and ε-FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core-mantle boundary. Here we have constrained the equation of state, including the effects of a spin crossover in the Fe 3+ atoms, of (Al, Fe)-phase H: Al 0.84 Fe 3+ 0.07 Mg 0.02 Si 0.06 OOH, using powder X-ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)-phase H and δ-(Al, Fe)OOH. The changes in spin state of Fe 3+ in (Al, Fe)-phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared with δ-(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)-phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)-phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid-lower mantle (~1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)-phase H could reproduce some of the seismic signatures of large, low seismic-velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ-AlOOH)-(MgSiO 2 (OH) 2 )-(ε-FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.

58 GEOSCIENCES↗

DOE OSTI2021

r-Process Radioisotopes from Near-Earth Supernovae and Kilonovae

The astrophysical sites where r-process elements are synthesized remain mysterious: it is clear that neutron star mergers (kilonovae (KNe)) contribute, and some classes of core-collapse supernovae (SNe) are also possible sources of at least the lighter r-process species. The discovery of 60Fe on the Earth and Moon implies that one or more astrophysical explosions have occurred near the Earth within the last few million years, probably SNe. Intriguingly, 244 Pu has now been detected, mostly overlapping with 60 Fe pulses. However, the 244 Pu flux may extend to before 12 Myr ago, pointing to a different origin. Here, motivated by these observations and difficulties for r-process nucleosynthesis in SN models, we propose that ejecta from a KN enriched the giant molecular cloud that gave rise to the Local Bubble, where the Sun resides. Accelerator mass spectrometry (AMS) measurements of 244 Pu and searches for other live isotopes could probe the origins of the r-process and the history of the solar neighborhood, including triggers for mass extinctions, e.g., that at the end of the Devonian epoch, motivating the calculations of the abundances of live r-process radioisotopes produced in SNe and KNe that we present here. Given the presence of 244 Pu, other r-process species such as 93 Zr, 107 Pd, 129 I, 135 Cs, 182 Hf, 236 U, 237 Np, and 247 Cm should be present. Their abundances and well-resolved time histories could distinguish between the SN and KN scenarios, and we discuss prospects for their detection in deep-ocean deposits and the lunar regolith. We show that AMS 129 I measurements in Fe–Mn crusts already constrain a possible nearby KN scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

CRITICAL MINERAL RESOURCE ASSESSMENT OF COALS FROM THE POWDER RIVER BASIN, WYOMING, USA

Coal and associated sediments are potential sources of Rare Earth Elements (REE) and other critical minerals necessary for many consumer and defense-related products and development of energy technologies. For this study, 15 samples were sourced from Dry Fork Mine within the Powder River Basin near Gillette, Wyoming, to determine REE content along with select transition metals (t-metals) of interest in non-combusted coals and associated sediments. An excavator collected samples along a mine highwall, within benched intervals representing a foot to a few feet of coal strata. Homogenized sample subsets were crushed in an aluminum-ceramic vial, then passed through a 106-micron sieve. Sieved samples were pyrolyzed by a LECO TGA-108, which concurrently measured sample moisture, volatile, and ash content. Residual ash was then fused into a glass bead with lithium metaborate at a 6:1 Li-metaborate to sample ratio and then dissolved in 20% nitric acid. The nitric acid solutions containing digested beads were then diluted and passed through a 0.45 micron syringe filter. Filtered samples were analyzed by ICP-MS for REE and select t-metals. Reagent blanks were analyzed to assess for any potential contamination. Total REE plus Y (REE+Y) concentrations across the 15 samples range from 52 ppm to 1015 ppm with an average of 477 ppm on a dry ash basis. When normalized to upper continental crust (UCC; Taylor and McLennan, 1985) values, 13 of the 15 samples show average enrichment factors (EF) for REE+Y greater than 1. EF values indicate M-type enrichment (Sm through Er), as opposed to distinct preferential enrichment of light or heavy REE. T-metals with the highest average concentrations are V, Cu, Ni, and Zr at 215, 243, 123, and 118 ppm, respectively. Multiple t-metals had EF>1; however, Cu, Ge, and Sb have the highest average EF relative to UCC for t-metals, at 9.7, 8.0, and 17.0, respectively.

Copeland, Grant↗

DOE OSTI2020

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗

DOE OSTI2023

Critical Minerals and Rare Earth Elements in Powder River Basin Coal and Associated Sediments, Wyoming, USA

The demand for Rare Earth Elements (REE) and other critical minerals (CM) required for consumer products, defense-related applications, and low-carbon energy technology is increasing rapidly. Unconventional sources of REE/CM, such as coal and associated sediments, could prove very important in building resilient and ethical energy supply chains. Utilization of existing infrastructure and the highly trained energy workforce in traditionally coal-producing regions in emerging REE/CM industries could provide an economic boost to coal communities. The Powder River Basin (PRB) of Wyoming and Montana, USA, is a prime candidate to investigate the feasibility of extracting REE/CM from coal and associated sediments. The PRB hosts thick (>50 ft) coal seams that are mined at the surface and more than 40% of the coal produced in the US comes from the PRB. Geochemical data from multiple locations in PRB coal systems show REE enrichments at the top and bottom margins and at internal partings in coal seams. This trend is exhibited in two cores drilled at Peabody’s North Antelope Rochelle mine in the east central PRB, the largest coal mine in the world. The two cores include the Wyodak Anderson coal zone as well as the over- and underlying shale units. The Wyodak Anderson coal zone is part of the Paleocene Tongue River Member of the Fort Union Formation in which most of the coal resources in the PRB reside. To understand REE/CM enrichment, a total of 188 samples were analyzed for their major and trace element chemistry from both cores, sampled at one foot or smaller intervals. The two cores show distinct REE enrichments in the uppermost and lowermost 3 to 12 feet of the Wyodak Anderson coal zone, as well as enrichments in bounding carbonaceous shale units. Total REE+Y (REY) concentrations as high as 2510 ppm (all concentrations reported on an ash basis), or 15 times average upper continental crust values (Taylor and McLennan, 1995), were identified in the coal. Partings within the coal zone also show relative REE enrichment, with concentrations up to 485 ppm REY. The remaining interior portions of the coal zone contain lower concentrations of REE, resulting in an average REY for all samples from both cores of 288 ppm. The proportion of high-value critical REE (Nd, Eu, Tb, Dy, Er, Y) compared to REY averages 36%, which is higher than the critical REE proportion in average upper continental crust of 33%. This finding highlights the relative enrichment of middle REE compared to light REE, in contrast to many conventional REE deposits. Other critical trace elements that show enrichments above average upper continental crust include Ga, Nb, and V. Within this sample suite, Ga concentrations range from 2.8 to 213 ppm, with an average of 29 ppm; Nb from 2.9 to 166 ppm with an average of 26 ppm; and V from 27 to 1695 ppm with an average of 188 ppm. None of these trace elements exhibit the same distinct pattern of enrichment as REE in the upper and lower bounding layers and internal partings of the coal zone. Major element chemistry indicates CaO concentrations averaging ~25% in the coal samples from both cores, which has important implications for the extractability of REE from coal. Calcium-rich PRB coal has been shown to be more amenable to REE extraction than lower calcium coals by some methods (Stuckman et al., 2019; Taggart et al., 2016), highlighting the importance of REE extractability in addition to REE concentration in assessing the value of potential feedstocks. The original coal resource in the PRB is estimated at 1.16 trillion short tons (Luppens et al., 2015). Thus, the PRB represents an important potential unconventional source for REE and other CM. Ongoing research to fully characterize the REE/CM resource, identify enrichment mechanisms, and further develop and scale extraction technologies is necessary to understand the full potential of this resource.

Phillips, Erin↗