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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Reduction of Marine Magnetic Data for Modeling the Main Field of the Earth

The marine data set archived at the National Geophysical Data Center (NGDC) consists of shipborne surveys conducted by various institutes worldwide. This data set spans four decades (1953, 1958, 1960-1987), and contains almost 13 million total intensity observations. These are often less than 1 km apart. These typically measure seafloor spreading anomalies with amplitudes of several hundred nanotesla (nT) which, since they originate in the crust, interfere with main field modeling. The source for these short wavelength features are confined within the magnetic crust (i.e., sources above the Curie isotherm). The main field, on the other hand, is of much longer wavelengths and originates within the earth's core. It is desirable to extract the long wavelength information from the marine data set for use in modeling the main field. This can be accomplished by averaging the data along the track. In addition, those data which are measured during periods of magnetic disturbance can be identified and eliminated. Thus, it should be possible to create a data set which has worldwide data distribution, spans several decades, is not contaminated with short wavelengths of the crustal field or with magnetic storm noise, and which is limited enough in size to be manageable for the main field modeling. The along track filtering described above has proved to be an effective means of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for main field modeling. Its simplicity and ability to adequately handle varying spatial and sampling constraints has outweighed consideration of more sophisticated approaches. This filtering technique also provides the benefits of smoothing out short wavelength crustal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be used in the least square modeling. A useful data set now exists which spans 1953-1987.

Baldwin, R. T.↗
NASA NTRS1992

Palinspastic reconstruction of southeastern California and southwestern Arizona for the middle Miocene

A paleogeographic reconstruction of southeastern California and southwestern Arizona at 10 Ma was made based on available geologic and geophysical data. Clockwise rotation of 39 deg was reconstructed in the eastern Transverse Ranges, consistent with paleomagnetic data from late Miocene volcanic rocks, and with slip estimates for left-lateral faults within the eastern Transverse Ranges and NW-trending right lateral faults in the Mojave Desert. This domain of rotated rocks is bounded by the Pinto Mountain fault on the north. In the absence of evidence for rotation of the San Bernardino Mountains or for significant right slip faults within the San Bernardino Mountains, the model requires that the late Miocene Pinto Mountain fault become a thrust fault gaining displacement to the west. The Squaw Peak thrust system of Meisling and Weldon may be a western continuation of this fault system. The Sheep Hole fault bounds the rotating domain on the east. East of this fault an array of NW-trending right slip faults and south-trending extensional transfer zones has produced a basin and range physiography while accumulating up to 14 km of right slip. This maximum is significantly less than the 37.5 km of right slip required in this region by a recent reconstruction of the central Mojave Desert. Geologic relations along the southern boundary of the rotating domain are poorly known, but this boundary is interpreted to involve a series of curved strike slip faults and non-coaxial extension, bounded on the southeast by the Mammoth Wash and related faults in the eastern Chocolate Mountains. Available constraints on timing suggest that Quaternary movement on the Pinto Mountain and nearby faults is unrelated to the rotation of the eastern Transverse Ranges, and was preceded by a hiatus during part of Pliocene time which followed the deformation producing the rotation. The reconstructed Clemens Well fault in the Orocopia Mountains, proposed as a major early Miocene strand of the San Andreas fault, projects eastward towards Arizona, where early Miocene rocks and structures are continuous across its trace. The model predicts a 14 deg clockwise rotation and 55 km extension along the present trace of the San Andreas fault during late Miocene and early Pliocene time. Palinspastic reconstructions of the San Andreas system based on this proposed reconstruction may be significantly modified from current models.

Richard, Stephen M.↗
NASA NTRS2023

Testing the Gravitational Redshift With an Inner Solar System Probe: the Veritas Case

The NASA Discovery-class mission VERITAS, selected in June 2021, will be launched towards Venus after 2027. In addition to the science instrumentation that will build global foundational geophysical datasets, VERITAS proposed to conduct a technology demonstration for the Deep Space Atomic Clock (DSAC-2). A first DSAC successfully operated in low-Earth orbit for more than two years, demonstrated the trapped ion atomic clock technology, and established a new level of performance for clocks in space. DSAC-2 would have further improvements in size, power, and performance. It would host a 1× 10 -13 grade USO to produce a frequency output with short-term stability of less than 2× 10 -13 /√ τ (where τ is the averaging time).However, due to funding shortfalls, DSAC-2, had to be canceled. The initially foreseen presence of an atomic clock on board the probe, however, raised the question whether this kind of instrumentation could be useful not only for navigation and time transfer but also for fundamental physics tests. In this work, we consider the DSAC-2 atomic clock and VERITAS mission as a specific example to measure possible discrepancies in the redshift predicted by General Relativity by using an atomic clock onboard an interplanetary spacecraft. In particular we investigate the possibility of measuring possible violations of the Local Lorentz Invariance and Local Position Invariance principles. We perform accurate simulations of the experiment during the VERITAS cruise phase. We consider different parametrizations of the possible violations of the General Relativity, different operational conditions, and several different assumptions on the expected measurement performance. We show that DSAC-2 onboard VERITAS would provide new and improved constraints with respect to the current knowledge. Our analysis shows the scientific value of atomic clocks like DSAC-2 hosted onboard interplanetary spacecraft.

Fabrizio De Marchi↗
NASA NTRS2004

Petrogenesis of Mare Basalts, Mg-Rich Suites and SNC Parent Magmas

The successful models for the internal evolution of the Moon must consider the volume, distribution, timing, composition and, ultimately, the petrogenesis of mare basaltic volcanism. Indeed, given the paucity of geophysical data, the internal state of the Moon in the past can be gleaned only be unraveling the petrogenesis of the various igneous products on the Moon and, particularly, the mare basalts. most useful in constraining the depth and composition of their source region [Delano, 1980] despite having undergone a certain degree of shallow level olivine crystallization.The bulk of the lunar volcanic glass suite can be modeled as the partial melting products of an olivine + orthopyroxene source region deep within the lunar mantle. Ti02 contents vary from 0.2 wt % -1 7.0wt [Shearer and Papike, 1993]. Values that extreme would seem to require a Ti- bearing phase such as ilmenite in the source of the high-Ti (but not in the VLT source) because a source region of primitive LMO olivine and orthopyroxene, even when melted in small degrees cannot account for the observed range of Ti02 compositions. The picritic glasses are undersaturated with respect to ilmenite at all pressures investigated therefore ilmenite must have been consumed during melting, leaving an ilmenite free residue and an undersaturated melt [Delano, 1980, Longhi, 1992, Elkins et al, 2000 among others]. Multi- saturation pressures for the glasses potentially represent the last depths at which the liquids equilibrated with a harzburgite residue before ascending to the surface. These occur at great depths within the lunar mantle. Because the liquids have suffered some amount of crystal fractionation, this is at best a minimum depth. If the melts are mixtures, then it is only an average depth of melting. Multisaturation, nevertheless, is still a strong constraint on source mineralogy, revealing that the generation of the lunar basalts was dominated by melting of olivine and orthopyroxene.

Hess, Paul C.↗
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↗
NASA NTRS2009

Proceedings of the 40th Lunar and Planetary Science Conference

The 40th Lunar and Planetary Science Conference included sessions on: Phoenix: Exploration of the Martian Arctic; Origin and Early Evolution of the Moon; Comet Wild 2: Mineralogy and More; Astrobiology: Meteorites, Microbes, Hydrous Habitats, and Irradiated Ices; Phoenix: Soil, Chemistry, and Habitability; Planetary Differentiation; Presolar Grains: Structures and Origins; SPECIAL SESSION: Venus Atmosphere: Venus Express and Future Missions; Mars Polar Caps: Past and Present; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part I; 5 Early Nebula Processes and Models; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Cosmic Gymnasts; Mars: Ground Ice and Climate Change; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part II; Chondrite Parent-Body Processes; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Salubrious Surfaces; SNC Meteorites; Ancient Martian Crust: Primary Mineralogy and Aqueous Alteration; SPECIAL SESSION: Messenger at Mercury: A Global Perspective on the Innermost Planet; CAIs and Chondrules: Records of Early Solar System Processes; Small Bodies: Shapes of Things to Come; Sulfur on Mars: Rocks, Soils, and Cycling Processes; Mercury: Evolution and Tectonics; Venus Geology, Volcanism, Tectonics, and Resurfacing; Asteroid-Meteorite Connections; Impacts I: Models and Experiments; Solar Wind and Genesis: Measurements and Interpretation; Mars: Aqueous Processes; Magmatic Volatiles and Eruptive Conditions of Lunar Basalts; Comparative Planetology; Interstellar Matter: Origins and Relationships; Impacts II: Craters and Ejecta Mars: Tectonics and Dynamics; Mars Analogs I: Geological; Exploring the Diversity of Lunar Lithologies with Sample Analyses and Remote Sensing; Chondrite Accretion and Early History; Science Instruments for the Mars Science Lander; . Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Volcanism; Early Solar System Chronology; Seek Out and Explore: Upcoming and Future Missions; Mars: Early History and Impact Processes; Mars Analogs II: Chemical and Spectral; Achondrites and their Parent Bodies; and Planning for Future Exploration of the Moon The poster sessions were: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1; LRO and LCROSS; Geophysical Analysis of the Lunar Surface and Interior; Remote Observation and Geologic Mapping of the Lunar Surface; Lunar Spectroscopy; Venus Geology, Geophysics, Mapping, and Sampling; Planetary Differentiation; Bunburra and Buzzard Coulee: Recent Meteorite Falls; Meteorites: Terrestrial History; CAIs and Chondrules: Records of Early Solar System Processes; Volatile and Organic Compounds in Chondrites; Crashing Chondrites: Impact, Shock, and Melting; Ureilite Studies; Petrology and Mineralogy of the SNC Meteorites; Martian Meteorites; Phoenix Landing Site: Perchlorate and Other Tasty Treats; Mars Polar Atmospheres and Climate Modeling; Mars Polar Investigations; Mars Near-Surface Ice; Mars: A Volatile-Rich Planet; Mars: Geochemistry and Alteration Processes; Martian Phyllosilicates: Identification, Formation, and Alteration; Astrobiology; Instrument Concepts, Systems, and Probes for Investigating Rocks and Regolith; Seeing is Believing: UV, VIS, IR, X- and Gamma-Ray Camera and Spectrometer Instruments; Up Close and Personal: In Situ Analysis with Laser-Induced Breakdown Spectroscopy and Mass Spectrometry; Jupiter and Inscrutable Io; Tantalizing Titan; Enigmatic Enceladus and Intriguing Iapetus; Icy Satellites: Cryptic Craters; Icy Satellites: Gelid Geology/Geophysics; Icy Satellites: Cool Chemistry and Spectacular Spectroscopy; Asteroids and Comets; Comet Wild 2: Mineralogy and More; Hypervelocity Impacts: Stardust Models, LDEF, and ISPE; Presolar Grains; Early Nebular Processes: Models and Isotopes; Solar Wind and Genesis: Measurements and Interpretation; Education and Public Outreach; Mercury; Pursuing Lunar Exploration; Sources and Eruptionf Lunar Basalts; Chemical and Physical Properties of the Lunar Regolith; Lunar Dust and Transient Surface Phenomena; Lunar Databases and Data Restoration; Meteoritic Samples of the Moon; Chondrites, Their Clasts, and Alteration; Achondrites: Primitive and Not So Primitive; Iron Meteorites; Meteorite Methodology; Antarctic Micrometeorites; HEDs and Vesta; Dust Formation and Transformation; Interstellar Organic Matter; Early Solar System Chronology; Comparative Planetology; Impacts I: Models and Experiments; Impacts II: Craters and Ejecta; Mars: Volcanism; Mars: Tectonics and Dynamics; Martian Stratigraphy: Understanding the Geologic History of Mars Through the Sedimentary Rock Record; Mars: Valleys and Valley Networks; Mars: Aqueous Processes in Valles Marineris and the Southern Highlands; Mars: Aqueous Geomorphology; Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Remote Sensing; Mars: Geologic Mapping, Photogrammetry, and Cratering; Martian Mineralogy: Constraints from Missions and Laboratory Investigations; Mars Analogs: Chemical and Physical; Mars Analogs: Sulfates and Sulfides; Missions: Approaches, Architectures, Analogs, and Actualities; Not Just Skin Deep: Electron Microscopy, Heat Flow, Radar, and Seismology Instruments and Planetary Data Systems, Techniques, and Interpretation.

Source record↗
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↗
DOE OSTI2022

Size Ranges of Magnetic Domain States in Tetrataenite

Paleomagnetic studies of meteorites provide unique constraints on the evolution of magnetic fields in the early solar system. These studies rely on the identification of magnetic minerals that can retain stable magnetizations over ≳4.5 billion years (Ga). The ferromagnetic mineral tetrataenite (γ''-Fe 0.5 Ni 0.5 ) is found in iron, stony-iron and chondrite meteorite groups. Nanoscale intergrowths of tetrataenite have been shown to carry records of paleomagnetic fields, although the effect of magnetostatic interactions on their magnetic remanence acquisition remains to be fully understood. Tetrataenite can also occur as isolated, non-interacting, nanoscale grains in many meteorite groups, although the paleomagnetic potential of these grains is particularly poorly understood. Here, we aim to improve our understanding of tetrataenite magnetization to refine our knowledge of existing paleomagnetic analyses and broaden the spectrum of meteorite groups that can be used for future paleomagnetic studies. We present the results of analytical calculations and micromagnetic modeling of isolated tetrataenite grains with various geometries. We find that tetrataenite forms a stable single domain state at grain lengths between 6 and ~160 nm dependent on its elongation. It also possesses a magnetization resistant to viscous remagnetization over the lifetime of the solar system at 293 K. At larger grain sizes, tetrataenite's lowest energy state is a lamellar two-domain state, stable at Ga-scale timescales. Unlike many other magnetic minerals, tetrataenite does not form a single-vortex domain state due to its large uniaxial anisotropy. Our results show that single domain and two-domain tetrataenite grains carry an extremely stable magnetization and therefore are promising for paleomagnetic studies.

58 GEOSCIENCES↗
DOE OSTI2025

Revealing the Moon's Taurus‐Littrow Landslide via Integrated Analysis of Pristine Apollo 17 Soil Core 73001/2

The “light mantle” deposit at the base of South Massif in the Moon's Taurus‐Littrow Valley was a primary science target for the Apollo 17 exploration. The possibility that it was a landslide triggered by ejecta from Tycho Crater is critical for establishing the age of Tycho and constraining recent lunar impact chronology; however, the mechanism of emplacement of the deposit has recently been questioned. The newly opened 73001/73002 double drive tube from Station 3 sampled 70.6 cm deep into the regolith and represents the first stratigraphic section of an extraterrestrial landslide deposit returned to Earth. Here we provide an overview of the stratigraphy of the 73001/73002 core based on top to bottom variations revealed by coordinated laboratory analyses and explore constraints on the emplacement of the light mantle deposit. Briefly, the upper ∼10 cm of 73002 contains a disturbed zone from space weathering and emplacement of ejecta from a nearby crater that excavated and ejected basaltic material. Below 10 cm is a nearly uniform unit of immature regolith. These data support a single event for the emplacement of the deposit at this location, followed by weathering and mixing of materials from nearby crater ejecta in the upper 10 cm. Slight variations in chemistry and clast components may reflect the relative stratigraphy of the South Massif slope, with material toward the bottom of 73001 originating from lower slopes and material from higher up in the core representing regolith from higher up the South Massif slopes.

58 GEOSCIENCES↗
DOE OSTI2026

Learning earthquake ground motions via conditional generative modeling

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

geophysics↗
DOE OSTI2024

On the predictability of turbulent fluxes from land: PLUMBER2 MIP experimental description and preliminary results

Accurate representation of the turbulent exchange of carbon, water, and heat between the land surface and the atmosphere is critical for modelling global energy, water, and carbon cycles in both future climate projections and weather forecasts. Evaluation of models' ability to do this is performed in a wide range of simulation environments, often without explicit consideration of the degree of observational constraint or uncertainty and typically without quantification of benchmark performance expectations. We describe a Model Intercomparison Project (MIP) that attempts to resolve these shortcomings, comparing the surface turbulent heat flux predictions of around 20 different land models provided with in situ meteorological forcing evaluated with measured surface fluxes using quality-controlled data from 170 eddy-covariance-based flux tower sites. Predictions from seven out-of-sample empirical models are used to quantify the information available to land models in their forcing data and so the potential for land model performance improvement. Sites with unusual behaviour, complicated processes, poor data quality, or uncommon flux magnitude are more difficult to predict for both mechanistic and empirical models, providing a means of fairer assessment of land model performance. When examining observational uncertainty, model performance does not appear to improve in low-turbulence periods or with energy-balance-corrected flux tower data, and indeed some results raise questions about whether the energy balance correction process itself is appropriate. In all cases the results are broadly consistent, with simple out-of-sample empirical models, including linear regression, comfortably outperforming mechanistic land models. In all but two cases, latent heat flux and net ecosystem exchange of CO 2 are better predicted by land models than sensible heat flux, despite it seeming to have fewer physical controlling processes. Land models that are implemented in Earth system models also appear to perform notably better than stand-alone ecosystem (including demographic) models, at least in terms of the fluxes examined here. The approach we outline enables isolation of the locations and conditions under which model developers can know that a land model can improve, allowing information pathways and discrete parameterisations in models to be identified and targeted for future model development.

54 ENVIRONMENTAL SCIENCES↗