SEARCH · Engineering Papers

Results for “crustal”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Probing Cretaceous-Paleogene crustal thickness in southern Tibet using quartz-zircon chronobarometry

Knowledge of the crustal thickness history of southern Tibet during the India-Asia collision is key to understanding what deformation mechanisms accommodated northward propagating crustal shortening. Thermoisotopic models suggest that a relatively thin margin of southern Asia persisted from ca. 200 to 45 Ma whereas trace element paleodepth proxies are interpreted to indicate a thinning phase from ca. 100 to 65 Ma when the crust reached ∼30 km-thickness. As neither of these methods directly measures crustal thickness, resolution of this conflict awaits development of a method that can. Here, in this study, we place bounds on the thickening history of southern Tibet using coupled quartz-zircon thermobarometry and U-Pb geochronology of granitoid plutons in the Gangdese batholith. We find that crustal thicknesses were as high (or higher) as 60–75 km at 65 Ma, or 10–15 Ma prior to the onset of hard continental collision. Magmatic inflation in the lower crust likely contributed to crustal thickening and thermal weakening of the lower crust, suggesting that pure shear in the lower crust was likely the primary accommodation mechanism for N-S shortening during collision. These data are in partial agreement with results of thermoisotopic models of crustal thickness but contradict thickness histories derived from empirical trace element proxies.

Chronobarometry↗

DOE OSTI2024

Seasonal modulation of crustal seismicity in northeastern Japan driven by snow load

Numerous studies have reported that surface hydrological loading can seasonally modulate seismicity rates at crustal depths. For example, substantial winter snow accumulation occurs across the Japanese Islands, and these snowy regions appear to have seasonally modulated the occurrence of previous large inland earthquakes. Therefore, it is important to investigate the impact of seasonal stress changes on crustal seismicity to deepen our understanding of earthquake generation. Here we constrain seasonal changes in the surface load across northeastern Japan using Global Navigation Satellite System surface displacements and evaluate the potential relationship between temporal trends in inland seismicity and estimated seasonal stress changes. The spatial distribution of the seasonal surface load is consistent with snow depth along the Sea of Japan. The inland seismicity beneath northeastern Japan is modestly modulated by the seasonal stress changes that are induced by the annual snow load. However, this seasonal response is weaker than that in other regions. This weak modulation may be due to the small surface-load-induced stress perturbation relative to the long-term-averaged stressing rate and/or the limited presence of crustal fluids to trigger seismicity in Japan.

58 GEOSCIENCES↗

DOE OSTI2024

Elastic thermobarometry on metapelites across the crustal section of the Famatinian Arc, Argentina

We applied elastic thermobarometry on garnet-bearing migmatites along two transects through the crustal section at Sierra Valle Fértil-La Huerta, Argentina. We performed quartz-in-garnet barometry and zircon-in-garnet thermometry on metapelites from different paleo-depths across the crustal section. Our work recovers entrapment pressures ranging from 240 to 1330 MPa and entrapment temperatures between 691 and 1574 °C. The entrapment conditions are broadly consistent with anticipated pressures and temperatures along the crustal section derived previously using conventional, thermodynamic thermobarometers. The quartz-in-garnet barometer reproduces those conventionally established entrapment conditions when samples only experienced conditions within the alpha-quartz stability field. Raman-derived pressures for samples that experienced beta-quartz reference conditions are commonly much higher than those established by conventional barometry. Samples that preserve compressive (positive) residual pressures best reproduce reference entrapment pressures. Entrapment temperatures show high variability and overestimation of temperature conditions compared to conventional results. These results indicate elastic thermobarometry over- or under-estimates crystallization conditions in rocks crystallized at high temperatures, as is common in the Famatinian Arc deep-crust. We suggest that modeling quartz behavior across the alpha–beta transition may present challenges, as does shape maturation, viscous deformation, and radiation damage in zircon.

58 GEOSCIENCES↗

DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗

DOE OSTI2021

The Global Fingerprint of Modern Ice-Mass Loss on 3-D Crustal Motion

Crustal motion generated by rapid ice-mass loss from Earth's glaciers and ice sheets has previously been considered in Global Navigational Satellite System (GNSS) analyses and numerical models across regions of ice retreat. However, the fingerprint of ice-mass loss is not limited to glaciated areas, but is characterized by a global pattern of 3-D crustal deformation. We compute “far-field” vertical and horizontal deformation rates that occurred in response to early 21st century mass flux from the Greenland and Antarctic Ice Sheets, global glaciers and ice caps, and associated ocean loading. We demonstrate that mass changes in the Greenland Ice Sheet and high latitude glacier systems each generated average crustal motion of 0.1–0.4 mm/yr across much of the Northern Hemisphere, with significant year-to-year variability in magnitude and direction. Horizontal motions associated with ice-mass loss exceed vertical rates in many far-field areas, and both should be considered in future analysis of GNSS measurements.

58 GEOSCIENCES↗

DOE OSTI2023

Genomics discovery of giant fungal viruses from subsurface oceanic crustal fluids

Abstract The oceanic igneous crust is a vast reservoir for microbial life, dominated by diverse and active bacteria, archaea, and fungi. Archaeal and bacterial viruses were previously detected in oceanic crustal fluids at the Juan de Fuca Ridge (JdFR). Here we report the discovery of two eukaryotic Nucleocytoviricota genomes from the same crustal fluids by sorting and sequencing single virions. Both genomes have a tRNATyr gene with an intron (20 bps) at the canonical position between nucleotide 37 and 38, a common feature in eukaryotic and archaeal tRNA genes with short introns (<100 bps), and fungal genes acquired through horizontal gene transfer (HGT) events. The dominance of Ascomycota fungi as the main eukaryotes in crustal fluids and the evidence for HGT point to these fungi as the putative hosts, making these the first putative fungi-Nucleocytoviricota specific association. Our study suggests active host-viral dynamics for the only eukaryotic group found in the subsurface oceanic crust and raises important questions about the impact of viral infection on the productivity and biogeochemical cycling in this ecosystem.

Bhattacharjee, Ananda S.↗

DOE OSTI2022

Host-region parameters for an adjustable model for crustal earthquakes to facilitate the implementation of the backbone approach to building ground-motion logic trees in probabilistic seismic hazard analysis

The backbone approach to constructing a ground-motion logic tree for probabilistic seismic hazard analysis (PSHA) can address shortcomings in the traditional approach of populating the branches with multiple existing, or potentially modified, ground-motion models (GMMs) by rendering more transparent the relationship between branch weights and the resulting distribution of predicted accelerations. To capture epistemic uncertainty in a tractable manner, there are benefits in building the logic tree through the application of successive adjustments for differences in source, path, and site characteristics between the host region of the selected backbone GMM and the target region for which the PSHA is being conducted. The implementation of this approach is facilitated by selecting a backbone GMM that is amenable to such host-to-target adjustments for individual source, path, and site characteristics. The NGA-West2 GMM of Chiou and Youngs (CY14) has been identified as a highly adaptable model for crustal seismicity that is well suited to such adjustments. Rather than using generic source, path, and site characteristics assumed appropriate for the host region, the final suite of adjusted GMMs for the target region will be better constrained if the host-region parameters are defined specifically on the basis of their compatibility with the CY14 backbone GMM. To this end, making use of a recently developed crustal shear-wave velocity profile consistent with CY14, we present an inversion of the model to estimate the key source and path parameters, namely the stress parameter and the anelastic attenuation. With these outputs, the effort in constructing a ground-motion logic tree for any PSHA dealing with crustal seismicity can be focused primarily on the estimation of the target-region characteristics and their associated uncertainties. The inversion procedure can also be adapted for any application in which different constraints might be relevant.

Stafford, Peter J.↗

DOE OSTI2020

Waveform Modeling and Simulation for Crustal Phases [Slides]

Goals for this project are to: numerically simulate wave propagation in layered crustal models with laterally-varying structure; Investigate seismic phases Pg and Pn; help understand the variability and complexity of these phases in recorded seismic data; and Use wavefield movies and "wavepaths" to suggest improvements for SALSA3D crustal tomography. In this presentation I will outline progress that has been made to date in an effort to modify a community research modeling code ("simulator")into a tool which can be used by GNDD researchers to study the effect of varying crustal velocity on wavefield propagation.

58 GEOSCIENCES↗

DOE OSTI2024

Seasonal Modulation of Crustal Seismicity in Northeastern Japan Driven by Snow Load

58 GEOSCIENCES↗

DOE OSTI2021

Interplay Between Crustal‐Scale Thrusting, High Metamorphic Heating Rates, and the Development of Inverted Thermal‐Metamorphic Gradients: Numerical Models and Examples From the Caledonides of Northern Scotland

Abstract Understanding the processes that produce high prograde metamorphic heating rates and the development of inverted metamorphic sequences in collisional thrust belts remains a fundamental challenge for tectonics and metamorphic petrology. New 2D finite element models of crustal‐scale thrusts with variable slip rates (10, 20, 35, 50 km Myr −1 ) are used to examine how thrust sheet emplacement contributes to these processes. In the models, average prograde heating rates of 31–118 °C Myr −1 are observed in the footwall, with maximum transient heating rates of ∼167 °C Myr −1 occurring at the highest slip rate. Also, thrust sheet emplacement produces an inverted thermal gradient in the model footwall. At slip rates of 10 km Myr −1 , heating magnitudes >200 °C are observed >7 km structurally beneath the thrust plane. At slip rates of 20–50 km Myr −1 , models produce thermal penetration depths of 4–5 km for similar heating magnitudes. Petrologically determined heating rates for thrust footwalls in the Scandian orogenic wedge of northern Scotland mostly yield rates consistent with model results (10–230 °C Myr −1 ). The model‐derived magnitude of inverted thermal gradients is also similar to those indicated in texturally determined deformation temperature transects across the Scandian orogenic wedge (100 °C–180 °C). Combined, these results indicate that high heating rates can be produced by fault slip at typical plate velocities. Additionally, this implies that crustal‐scale thrusts are likely to produce inverted metamorphic sequences in most systems, provided that P‐T conditions during slip allow for metamorphic or recrystallization processes that would preserve evidence for such features.

Thigpen, J. Ryan↗

DOE OSTI2024

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Challenges and Potential of Waveform Modeling for Crustal Scale Predictions

Waveform modeling is crucial to improving our understanding of observed seismograms. Forward simulation of wavefields provides quantitative methods of testing interactions between complicated source functions and the propagation medium. Here, we discuss three experiments designed to improve under standing of high frequency seismic wave propagation. First, we compare observed and predicted travel times of crustal phases for a set of real observed earthquakes with calculations and synthetic seismograms. Second, we estimate the frequency content of a series of nearly co-located earthquakes of varying magnitude for which we have a relatively well- known 1D velocity model. Third, we apply stochastic perturbations on top of a 3D tomographic model and qualitatively assess how those variations map to differences in the seismograms. While different in scope and aim, these three vignettes illustrate the current state of crustal scale waveform modeling and the potential for future studies to better constrain the structure of the crust.

58 GEOSCIENCES↗

DOE OSTI2024

Advancing the Limits of InSAR to Detect Crustal Displacement from Low-Magnitude Earthquakes through Deep Learning

Detecting surface deformation associated with low-magnitude (M w ≤ 5) seismicity using interferometric synthetic aperture radar (InSAR) is challenging due to the subtlety of the signal and the often challenging imaging environments. However, low-magnitude earthquakes are potential precursors to larger seismic events, and thus characterizing the crustal displacement associated with them is crucial for regional seismic hazard assessment. We combine InSAR time-series techniques with a Deep Learning (DL) autoencoder denoiser to detect the magnitude and extent of crustal deformation from the M w = 3.4 Gallina, New Mexico earthquake that occurred on 30 July 2020. Although InSAR alone cannot detect event-related deformation from such a low-magnitude seismic event, application of the DL method reveals maximum displacements as small as (±2.5 mm) in the vicinity of both the fault and earthquake epicenter without prior knowledge of the fault system. This finding improves small-scale displacement discernment with InSAR by an order of magnitude relative to previous studies. We additionally estimate best-fitting fault parameters associated with the observed deformation. The application of the DL technique unlocks the potential for low-magnitude earthquake studies, providing new insights into local fault geometries and potential risks from higher-magnitude earthquakes. This technique also permits low-magnitude event monitoring in areas where seismic networks are sparse, allowing for the possibility of global fault deformation monitoring.

58 GEOSCIENCES↗

DOE OSTI2020

Crustal Structure of the Mesopotamian Plain, East of Iraq

The crustal structure of Iraq was investigated through analyzing teleseismic data from 12 new seismic stations. Three seismic stations are located within the Zagros Fold-Thrust Belt, and nine are within the Mesopotamian Plain. Joint inversion of P wave receiver function and Rayleigh wave dispersion data were employed to resolve S wave velocity structure models beneath each station. Combining these models with available Moho depths within and around the study area reveals that the Moho depth is smoothly increasing from the Arabian platform toward the Zagros Mountains. Further, an exceptional deeper root is observed near the eastern edge of the Mesopotamian Plain where sedimentary pile is the thickest. This root was interpreted as a structure inherited from the successive Mesozoic rifting of the NE Arabian platform enhanced by sedimentary loading and progressive subsidence. A low-velocity uppermost mantle was observed beneath the Arabian Foredeep attesting the existence of a low-strength lithospheric mantle beneath the region southwest from the Zagros deformation front. The weak rheology of the uppermost mantle may have allowed local sinking of the crust and deepening of the Moho boundary due to vertical loading and the Late Miocene lateral contraction.

58 GEOSCIENCES↗

DOE OSTI2021

Experimental Investigation of Mercury's Magma Ocean Viscosity: Implications for the Formation of Mercury's Cumulate Mantle, Its Subsequent Dynamic Evolution, and Crustal Petrogenesis

Mercury has a compositionally heterogeneous surface that was produced by different periods of igneous activity during Mercury’s history, perhaps suggesting heterogeneous mantle sources. Furthermore, understanding the structure of Mercury’s mantle formed during the planet’s magma ocean stage could help in developing a petrologic model for Mercury, and thus, its dynamic history in the context of crustal petrogenesis. We present results of falling sphere viscometry experiments on late stage Mercurian magma ocean analogue compositions. Owing to the presence of sulfur on the surface of Mercury, two compositions were tested, one with sulfur and one without. The liquids have viscosities of 0.6-3.9 (sulfur-bearing) and 0.6-10.9 Pa·s (sulfur-free), similar to an andesite at the tested experimental conditions (1.4-6.2 GPa, 1600-2000°C). We present new viscosity models that enable extrapolation beyond the experimental conditions and evaluate grain growth and the potential for crystal entrainment in a cooling, convecting magma ocean. We consider scenarios with and without a graphite flotation crust, which suggests two possible endmember outcomes for Mercury’s mantle structure. With a graphite flotation crust, crystallization of the mantle would be fractional with negatively buoyant minerals sinking to form a stratified cumulate pile according to the crystallization sequence. Without a flotation crust, crystals would largely remain entrained in the convecting liquid during solidification, producing a homogeneous mantle. In the context of these endmember models, the chemically heterogeneous surface could result from dynamical stirring or mixing of a mantle that was initially compositionally stratified, or from different extents of melting of a homogeneous mantle.

58 GEOSCIENCES↗

DOE OSTI2022

Trans-crustal structural control of CO 2 -rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

58 GEOSCIENCES↗

DOE OSTI2020

Modulation of seismic noise near the San Jacinto fault in southern California: origin and observations of the cyclical time dependence and associated crustal properties

SUMMARY We examine the cyclic amplitude variation of seismic noise recorded by continuous three-component broad-band seismic data with durations spanning 91–713 d (2008–2011) from three different networks: Anza seismic network, IDA network and the transportable seismic array. These stations surround the San Jacinto Fault Zone (SJFZ) in southern California. We find the seismic noise amplitudes exhibit a cyclical variation between 0.3 and 7.2 Hz. The high-frequency (≥0.9 Hz) noise variations can be linked to human activity and are not a concern. Our primary interest is signals in the low frequencies (0.3–0.9 Hz), where the seismic noise is modulated by semi-diurnal tidal mode M2. These long-period (low-frequency) variations of seismic noise can be attributed to a temporal change of the ocean waves breaking at the shoreline, driven by ocean tidal loading. We focus on the M2 variation of seismic noise at f = 0.6 Hz, travelling distances of ∼92 km through the crust from offshore California to the inland Anza, California, region. Relative to the shoreline station, data from the inland stations show a phase lag of ∼ –12°, which we attribute to the cyclic change in M2 that can alter crustal seismic attenuation. We also find that for mode M2 at 0.6 Hz, the amplitude variations of the seismic quality factor (Q) depend on azimuth and varies from 0.22 per cent (southeast to northwest) to 1.28 per cent (northeast to southwest) with Q = 25 for Rayleigh waves. We propose the direction dependence of the Q variation at 0.6 Hz reflects the preferred orientation of subfaults parallel to the main faulting defined by the primarily N45°W strike of the SJFZ.

Geochemistry & Geophysics↗

DOE OSTI2022

Crustal and uppermost mantle structure of Cape Verde from ambient noise tomography

SUMMARY We present a seismic ambient noise tomography of the Cape Verde archipelago, located in the Atlantic Ocean, approximately 600 km west of Senegal. We used 38 seismic broad-band stations that continuously recorded for 10 months, in order to construct the first 3-D model of Sv-wave velocities for the crust and uppermost mantle beneath the Cape Verde region. We started by computing phase cross-correlations for vertical component recordings using all possible inter-island station pairs. Next, a time–frequency phase-weighted stack was applied to obtain robust Rayleigh-wave group-velocity dispersion curves in the period band between 10 and 24 s. Group-velocity maps at different periods are obtained by inverting the dispersion curves. We then inverted the group-velocity maps to obtain the 3-D shear wave velocity structure of the crust and uppermost mantle beneath Cape Verde. The final 3-D model extends from 8 km down to 23 km and has a lateral resolution of about 50 km. The crust in the southwestern sector, encompassing Fogo, presents lower S-wave velocities that may be caused by the presence of melt pockets and/or hydrothermal fluids circulation. The uppermost mantle beneath the northwestern sector is characterized by higher S-wave velocities in agreement with previous results obtained from Ps and Sp receiver functions. Those high-velocity anomalies can reflect non-altered crust or remnants of magma chambers or solidified basaltic intrusions, which fed the volcanism in these islands. Our maps revealed the presence of crustal underplating across the entire archipelago, yet stronger beneath the groups Santo Antão—São Vicente—São Nicolau and Fogo—Santiago—Maio.

Carvalho, J. (ORCID:0000000310778994)↗