Search NASA⌕ Search

SEARCH · Search NASA

Results for “tectonic processes”

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.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Gordon Research Conference Technical Report

This grant supported travel and conference fees for participants in the 2022 Gordon Research Conference (GRC) and Gordon Research Seminar (GRS) on Rock Deformation. The theme of the conference was “Combining Laboratory Measurements with Observational Constraints to Understand Tectonic Processes”, which focused on how to scale up observations from laboratory experiments to tectonic processes. Dan Faulkner (Univ. of Liverpool) was the GRC chair and Heather Savage (UC Santa Cruz) was the vice chair. The GRS was run by Carolyn Tewksbury-Christle and Taka Kanaya. The majority of the funds from this grant supported early career researchers.

58 GEOSCIENCES↗
DOE OSTI2024

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗
DOE OSTI2026

Corundum Discovered by SuperCam and the Perseverance Rover at Jezero Crater, Mars

Minerals provide clues to the type of geologic environments in which they form. The primary types of minerals that have been identified on Mars indicate a history dominated by volcanism followed by interaction with water. However, other minerals are sometimes found. Here we present results from the SuperCam instrument on the Perseverance rover in which the mineral corundum (α-Al2O3) has been identified using a technique called time-resolved luminescence spectroscopy. Corundum has a unique light-emitting (luminescence) feature based on two well-known peaks as well as the time it takes for this light to fade after being excited by a laser, which is referred to as the lifetime. Corundum forms in specific environments that are enriched in aluminum and depleted in silicon and have been subjected to high temperature and pressure. Such places include deep subsurface magmas or shallower locations that have been exposed to high temperature through tectonic processes or meteorite impacts. Given their small size, association with a mineral called plagioclase, and location on the heavily impacted crater rim, we propose that these corundum minerals formed through impact processes.

58 GEOSCIENCES↗
DOE OSTI2025

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗
DOE OSTI2023

Toward Global Regional Seismic Moment Tensor Inversion with Three-Dimensional Earth Models for Nuclear Explosion Monitoring with Sparse Networks: Demonstration of Reciprocity for Strain Greens Tensor Database Simulation with Salvus

Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).

58 GEOSCIENCES↗
DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2023

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

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

Geochemistry & Geophysics↗
DOE OSTI2024

A Complete Machine-Learning-Based Workflow to Illuminate Earthquake Processes

Under this grant we developed, tested, and made available, machine learning models to improve the tasks in the earthquake monitoring workflow (Figure 1). We implemented these models as a part of an end-to-end workflow for seismic network processing and demonstrated, in a variety of settings, that these methods generalize and that they result in dramatically more comprehensive earthquake catalogs. These catalogs illuminate earthquake processes in detail and to an extent that had previously not been possible, and they do so for both tectonic seismicity and seismicity induced by fluid injection related to unconventional hydrocarbon development. This report summarizes the results from the 15 publications that resulted from this grant. Those contributes are divided into: (1) the development of specific tasks related to monitoring (7 publications), (2) the organization of those tasks into workflows for seismic monitoring (2 publications), and (3) applications of those workflows to data (6 publications).

58 GEOSCIENCES↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2024

EQSIM: Exascale Predictions of Earthquake Effects on Critical Infrastructure

The great “San Francisco” earthquake of 1906 is one of the most recognized, and sobering, demonstrations of the havoc that can be caused by the sudden and violent movement of Earth’s tectonic plates. The estimated 7.9-magnitude quake and subsequent fires decimated the major metropolis and surrounding areas: buildings turned to ruins, hundreds of thousands of people left homeless, and a death toll exceeding 3,000. Today, as evidenced by the catastrophic 7.8-magnitude earthquake that struck Turkey in February 2023, these events still present a significant danger to life and economic security. To mitigate the potential devastation of future earthquakes and better prepare for these inevitable events, researchers are turning to high-performance computers to simulate the underlying geophysical processes and accurately quantify associated risks to critical infrastructure.

42 ENGINEERING↗
DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗
DOE OSTI2026

Constructing a High‐Resolution Aftershock Catalog for the 2017 Mw 8.2 Tehuantepec Earthquake Sequence Using a Machine Learning–Based Workflow

The 8 September 2017 Mw 8.2 Tehuantepec earthquake was the largest instrumentally recorded normal‐faulting earthquake in Mexico. The mainshock occurred offshore within the Tehuantepec seismic gap, generating >30,000 aftershocks in the following year. We applied an open‐source, machine learning (ML)–assisted workflow to construct a high‐resolution aftershock catalog using data from temporary and permanent seismic networks in southern Mexico. The workflow integrates PhaseNet for phase detection; GaMMA for phase association; and VELEST, HypoInverse, and HypoDD for velocity modeling and relocation. We processed seven months of continuous waveform data from 29 broadband stations, including a temporary rapid‐response deployment that improved station coverage of the offshore rupture zone. To evaluate performance, we compared our results against analyst‐reviewed picks and event locations from the Servicio Sismológico Nacional catalog. The resulting catalog contains 11,374 relocated earthquakes and represents the most comprehensive published dataset for this sequence, incorporating the first full use of the temporary network. Relocated hypocenters show improved depth control and align well with the Slab2.0 subduction geometry, revealing clearer separation between offshore slab events and onshore crustal seismicity. This study demonstrates that combining ML‐based detection with established methods provides a scalable and reproducible approach for constructing high‐quality earthquake catalogs in tectonically complex environments and offers practical guidance for adapting similar workflows to other earthquake sequences.

Garcia, Marc [The University of Texas at El Paso, ↗
DOE OSTI2025

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES↗
DOE OSTI2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics↗
DOE OSTI2026

Long‐term growth and persistence of granitic inselbergs in a semi‐arid cratonic landscape

Granitic inselbergs are enduring landforms that punctuate planation surfaces in tectonically stable settings, yet the mechanisms and timescales of their persistence remain debated. Here, in this study, we present a multinuclide cosmogenic dataset ( 10 Be, 26 Al and 21 Ne) from bedrock, fluvial sediments and planation surfaces in the semi-arid Central Ceará Domain, north-eastern Brazil, to quantify denudation rates and reconstruct the long-term dynamics of residual relief formation. This triple-nuclide approach enables independent validation of denudation contrasts and long-term exposure histories, providing quantitative constraints on inselberg evolution that a single-nuclide approach could not resolve. Bedrock denudation rates (1–4 m Myr −1 ) are systematically lower than basin-wide rates (8–14 m Myr −1 ), revealing a regional pattern of vertical differential denudation. When combined with summit elevations and surface exposure ages of up to 5 Myr, these contrasts suggest that cumulative relief growth has been operating over timescales of tens of millions of years (1–57 Myr). Morphological and structural observations further indicate localized mass wasting and heterogeneous erosion styles, modulated by rock fabric and joint density. Cosmogenic nuclide inventories in detrital sediments reveal partial reworking from older sediment stocks, including surface pebbles with exposure ages up to 5 Myr, which record multistage of burial and re-exposure histories and reflect episodic remobilization in low-connectivity, semi-arid catchments. Rather than static relicts of ancient surfaces, inselbergs in this cratonic setting emerge as dynamic landforms actively sustained by slow but persistent differential denudation. These results provide new constraints on long-term cratonic landscape evolution and underscore the interplay between lithological resistance, structural inheritance and sediment transfer processes in sustaining topographic relief. They also call for a critical reassessment of classical inselberg models, suggesting that hybrid mechanisms, rather than single-process paradigms, more accurately capture the complexity of residual landform development in cratonic interiors.

Cosmogenic nuclides↗
DOE OSTI2025

Evaluation of Station Performance of the Idaho National Laboratory Seismic Monitoring Network Using Network Detection Thresholds

The Idaho National Laboratory (INL) Seismic Monitoring Network is located in eastern Idaho and monitors a portion of the intermountain seismic belt. It has been in place for 50 yr and has undergone several major changes, the most recent of which has been the transition to the Antelope real‐time acquisition system and the implementation of automatic phase picking algorithms to aid in analysis. This study discusses the efforts to evaluate the performance of the INL seismic monitoring network (and other surrounding stations) using the new real‐time acquisition system. The method outlined by Wilson et al. (2021) is used to develop an empirical relationship between the observability of local earthquakes as a function of magnitude and distance. This relationship is used to produce detection thresholds for Pwaves for all stations of interest. The INL seismic network has two main goals: monitor tectonic‐and volcanic‐related events and measure ground motions for input into seismic hazard analysis. Because of these two overall objectives, several seismic stations have been installed near critical facilities and, therefore, are not as quiet as stations that are used primarily for earthquake detection. This is reflected in their detection thresholds, which are much smaller for stations away from facilities. This study shows that the INL Seismic Monitoring Network is able to detect earthquakes near INL facilities with M L > 1.2, with redundancies built in to ensure this sensitivity even if data became unavailable from some stations. This study also shows “holes” in the monitoring network where the detection of smaller earthquakes is highly dependent on sparsely placed seismic stations. In conclusion, the results of this study will be used to govern plans for expansion of earthquake monitoring in Idaho and the surrounding region and to fine‐tune the detection thresholds for individual stations.

58 - GEOSCIENCES↗
DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗