Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical constraints”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of Salinity on Ground Ice Distribution Across an Arctic Coastal Polygonal Tundra Environment

The heterogeneous distribution of ground ice in the Arctic is a key driver of uneven ground subsidence as permafrost thaws, significantly impacting infrastructure and surface/subsurface hydrology. These topographic and hydrological changes contribute to major uncertainties in energy and carbon fluxes and storage in a warming Arctic. This study aims to improve our understanding of the controls on ground ice and organic matter distribution within the top 3 m of permafrost in coastal polygonal tundra near Utqiagvik, Alaska. To this end, we apply a neural network approach to bulk density distributions derived from nondestructive X-ray tomography of soil cores, trained with laboratory analyses, to improve the resolution and spatial coverage of estimates of dry bulk density, ice content, and organic matter content. In addition, we use capacitively coupled geophysical imaging to map soil electrical conductivity and salinity variations. The results show that sedimentary deposits from ocean transgressions, along with subsequent ice wedge polygon geomorphological processes, jointly influence the distribution of ice content at various scales. The impact of the latter decreases with depth, whereas the influence of salinity and sedimentary history increases. Although the controls on the distribution of soil organic matter content (g/cm 3 ) remain unclear, the pronounced heterogeneity in bulk density strongly influences its calculation from laboratory mass fraction measurements (g/g). From a methodological perspective, the interdependencies among soil components and the need for increased data coverage underscore the value of high-resolution density measurements, such as using X-ray tomography. Overall, this study emphasizes the importance of considering salinity constraints on ice content distribution in coastal permafrost regions. The results are expected to aid in the development of data products and process representations in geomorphological and ecosystem models.

Dafflon, Baptiste [Lawrence Berkeley National Labo↗
DOE OSTI2022

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗
DOE OSTI2020

Inversion of Time-Lapse Seismic Reservoir Monitoring Data Using CycleGAN: A Deep Learning-Based Approach for Estimating Dynamic Reservoir Property Changes

Carbon capture and storage is being pursued globally as a geoengineering measure for reducing the emission of anthropogenic CO 2 the atmosphere. Comprehensive monitoring, verification, and accounting programs must be established for demonstrating the safe storage of injected CO 2 . One of the most commonly deployed monitoring techniques is time-lapse seismic reservoir monitoring (also known as 4-D seismic), which involves comparing 3-D seismic survey data taken at the same study site but over different times. Analyses of 4-D seismic data volumes can help improve the quality of storage reservoir characterization, track the movement of injected CO 2 plume, and identify potential CO 2 spillover/leakage from the storage reservoirblue. However, the derivation of high-resolution CO 2 saturation maps from 4-D seismic data is a highly nonlinear and ill-posed inverse problem, often requiring significant computational effort. In this research, we apply a physics-based deep learning method to facilitate the solution of both the forward and inverse problems in seismic inversion while honoring physical constraints. A cycle generative adversarial neural network (CycleGAN) model is trained to learn the bidirectional functional mappings between the reservoir dynamic property changes and seismic attribute changes, such that both forward and inverse solutions can be obtained efficiently from the trained model. We show that our CycleGAN-based approach not only improves the reliability of 4-D seismic inversion but also expedites the quantitative interpretation. Our deep learning-based workflow is generic and can be readily used for reservoir characterization and reservoir model updates involving the use of 4-D seismic data.

58 GEOSCIENCES↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

An End-To-End Earthquake Detection Method for Joint Phase Picking and Association Using Deep Learning

Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. Here, in this study, we introduce a new end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations to associate and detect earthquakes across a seismic network. We use these sub-networks together with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the neural network model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches. Because our approach preserves information across tasks in the detection pipeline, it has the potential to outperform approaches that do not.

58 GEOSCIENCES↗
DOE OSTI2023

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗
DOE OSTI2021

Connectivity-informed drainage network generation using deep convolution generative adversarial networks

Abstract Stochastic network modeling is often limited by high computational costs to generate a large number of networks enough for meaningful statistical evaluation. In this study, Deep Convolutional Generative Adversarial Networks (DCGANs) were applied to quickly reproduce drainage networks from the already generated network samples without repetitive long modeling of the stochastic network model, Gibb’s model. In particular, we developed a novel connectivity-informed method that converts the drainage network images to the directional information of flow on each node of the drainage network, and then transforms it into multiple binary layers where the connectivity constraints between nodes in the drainage network are stored. DCGANs trained with three different types of training samples were compared; (1) original drainage network images, (2) their corresponding directional information only, and (3) the connectivity-informed directional information. A comparison of generated images demonstrated that the novel connectivity-informed method outperformed the other two methods by training DCGANs more effectively and better reproducing accurate drainage networks due to its compact representation of the network complexity and connectivity. This work highlights that DCGANs can be applicable for high contrast images common in earth and material sciences where the network, fractures, and other high contrast features are important.

42 ENGINEERING↗
DOE OSTI2023

Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models

Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.

Geochemistry & Geophysics↗
DOE OSTI2024

Wasatch Fault Structure from Machine Learning Arrival Times and High-Precision Earthquake Locations

Abstract On 18 March 2020, a magnitude 5.7 earthquake hit the Salt Lake valley in the state of Utah, United States. Using a dense geophone deployment and machine learning (ML), an additional several thousand events were detected and located. Currently, both the mainshock and the majority of the aftershocks are suspected to have occurred on or near a deeper portion of the Salt Lake segment of the Wasatch fault—part of a large range-bounding fault system thought to be capable of generating an Mw 7.2 earthquake. However, a small subset of aftershocks may have occurred on a portion of the more steeply, eastward dipping, and poorly understood West Valley fault. Unfortunately, the catalog locations and lack of focal mechanisms for this subset of aftershocks provide only a crude constraint on the true fault structure. To better illuminate fault structure, we relocate the ML-generated catalog with a range of magnitudes from −2 to 4.6, using: (1) NonLinLoc, a nonlinear location algorithm, (2) source-specific station terms, and (3) waveform coherence. We further compute first-motion focal mechanisms for 68 events. Results of the relocation suggest a simpler, minimally listric Wasatch fault geometry, contrary to what has been previously proposed. We also find that analysis of the focal mechanisms and waveform similarity indicates minimal event similarity throughout the Magna sequence, suggesting a highly complex and heterogeneous rupture zone, as opposed to rupture on a single plane. These findings suggest an increased seismic hazard due to the overall shallowness of the earthquake sequence and highly varied rupture mechanisms.

Geochemistry & Geophysics↗
DOE OSTI2024

Sulfur speciation in dacitic melts using X-ray absorption near-edge structure spectroscopy of the S K -edge (S-XANES): Consideration of radiation-induced changes and the implications for sulfur in natural arc systems

Abstract The synchrotron technique of micro X-ray absorption near-edge structure spectroscopy at the sulfur K-edge (S-XANES) provides a unique opportunity to measure the proportion of different oxidation states of sulfur (S) in silicate glasses. Although applied extensively in the analysis of basaltic silicate glasses, few S-XANES studies have investigated variations in S oxidation states with fO2 in felsic silicate glasses. In addition, no study has systematically compared the S-XANES results obtained from the same samples at different photon flux densities to quantify the relationship between exposure time and changes in S speciation in silicate glass, as has been done for Fe and V. This study evaluates observed differences in S speciation measured in experimentally produced H2O-saturated dacitic glasses over a range of reducing to oxidizing conditions (from log fO2 = ΔFMQ-0.7 to ΔFMQ+3.3; FMQ is the fayalite-magnetite-quartz mineral redox buffer) and equilibrated at 1000 °C and 300 MPa. S-XANES spectra were collected at three different photon flux densities using three microspectroscopy beamlines. As is observed in S-XANES analyses of basaltic silicate glasses, beam-induced changes to the S6+/ΣS are observed as a function of photon flux density and beam exposure time. Our results demonstrate that silicate glasses of dacitic composition undergo beam-induced photo-reduction in samples equilibrated at ΔFMQ > +1.75 and photo-oxidation if equilibrated at ΔFMQ < +1. The time required to observe beam-induced changes in the spectra varies as a function of flux density, and our study establishes an upper photon density limit at ~1.0 × 1012 photons/µm2. The S6+/ΣS calculated from spectra collected below this absorbed photon limit at intermediate flux densities (~1–4 × 109 photons/s per µm2) are affected by beam damage, as no conditions were found to be completely free of beam-induced changes. However, the S6+/ΣS ratios calculated below the limit at intermediate flux densities are consistent with thermodynamic constraints, demonstrating that S6+/ΣS ratios calculated from S-XANES spectra can be considered reliable for estimating the oxygen fugacity. Our results carry important implications for the S budget of felsic magmas and dissolution mechanisms in evolved melts. While our results from all three flux densities show the presence of S4+ dissolved in relatively oxidized (ΔFMQ > +1.75) dacitic glass, even in the spectra exposed to the lowest photon densities, we are unable to rule out the possibility that the S4+ signal is the result of instantaneous X-ray irradiation induced beam damage using S-XANES alone. When our spectra are compared to S-XANES spectra from basaltic silicate glasses, important differences exist in the solubility of S2– and S6+ between dacitic silicate melts, pointing to differences in solubility mechanisms as melt composition changes. This study highlights the need for further investigation into beam damage systematics, presence of S4+, and the solubility mechanisms of different oxidation states of S as silicate melt composition changes.

Geochemistry & Geophysics↗
DOE OSTI2024

A role for subducting clays in the water transportation into the Earth’s lower mantle

Subducting sedimentary layer typically contains water and hydrated clay minerals. The stability of clay minerals under such hydrous subduction environment would therefore constraint the lithology and physical properties of the subducting slab interface. Here we show that pyrophyllite (Al 2 Si 4 O 10 (OH) 2 ), one of the representative clay minerals in the alumina-silica-water (Al 2 O 3 -SiO 2 -H 2 O, ASH) system, breakdowns to contain further hydrated minerals, gibbsite (Al(OH) 3 ) and diaspore (AlO(OH)), when subducts along a water-saturated cold subduction geotherm. Such a hydration breakdown occurs at a depth of ~135 km to uptake water by ~1.8 wt%. Subsequently, dehydration breakdown occurs at ~185 km depth to release back the same amount of water, after which the net crystalline water content is preserved down to ~660 km depth, delivering a net amount of ~5.0 wt% H 2 O in a phase assemblage containing δ-AlOOH and phase Egg (AlSiO 3 (OH)). Our results thus demonstrate the importance of subducting clays to account the delivery of ~22% of water down to the lower mantle.

58 GEOSCIENCES↗
DOE OSTI2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2021

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
DOE OSTI2025

Dust Direct Radiative Effect Including Large Particles and Component Minerals

The direct radiative effect (DRE) of dust aerosols in Earth system models (ESMs) remains highly uncertain, largely due to inadequate representations of particle size distribution (PSD) and mineral composition. Using NASA's Earth Surface Mineral Dust Source Investigation (EMIT) soil mineralogy data and observed PSD in an ESM that resolves dust mineral composition and emitted diameters from 0.1 to 70 μm, we find a near-neutral global dust net DRE (−0.057 W m −2 ), weaker than most previous estimates. Large dust (diameter >10 μm) contributes 30% of the global longwave dust optical depth, providing observational constraints on large-particle abundance, and offsets 20% of the dust shortwave cooling over major source regions. Incorporation of EMIT mineralogy reduces shortwave uncertainty by more than 50%. The remaining uncertainty mainly exists in processes controlling dust abundance, particularly the poorly understood transport of large dust, and longwave optical properties, which require additional observational constraints to more accurately quantify the dust DRE.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2023

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗
DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗
DOE OSTI2023

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗