Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 19 records

Effect of tectonic processes on biosphere–geosphere feedbacks across a convergent margin

The subsurface is among Earth’s largest biomes, but the extent to which microbial communities vary across tectonic plate boundaries or interact with subduction-scale geological processes remains unknown. Here we compare bacterial community composition with deep-subsurface geochemistry from 21 hot springs across the Costa Rican convergent margin. We find that cation and anion compositions of the springs reflect the dip angle and position of the underlying tectonic structure and also correlate with the bacterial community. Co-occurring microbial cliques related to cultured chemolithoautotrophs that use the reverse tricarboxylic acid cycle (rTCA) as well as abundances of metagenomic rTCA genes correlate with concentrations of slab-volatilized carbon. This, combined with carbon isotope evidence, suggests that fixation of slab-derived CO2 into biomass may support a chemolithoautotrophy-based subsurface ecosystem. We calculate that this forearc subsurface biosphere could sequester 1.4 × 10 9 to 1.4 × 10 10 mol of carbon per year, which would decrease estimates of the total carbon delivered to the mantle by 2 to 22%. Furthermore, based on the observed correlations, we suggest that distribution and composition of the subsurface bacterial community are probably affected by deep tectonic processes across the Costa Rican convergent margin and that, by sequestering carbon volatilized during subduction, these chemolithoautotrophic communities could in turn impact the geosphere.

54 ENVIRONMENTAL SCIENCES↗

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↗

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗

The Role of Premagmatic Rifting in Shaping a Volcanic Continental Margin: An Example From the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2-D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along-strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post-rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn-rift structures, and the SDR there extend ~200 km seaward of the hinge line. Early post-rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn-rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

58 GEOSCIENCES↗

Monitoring Seismic Velocity Changes Across the San Jacinto Fault Using Train‐Generated Seismic Tremors

Abstract Microseismic noise has been used for seismic velocity monitoring. However, such signals are dominated by low‐frequency surface waves that are not ideal for detecting changes associated with small tectonic processes. Here we show that it is possible to extract stable, high‐frequency body waves using seismic tremors generated by freight trains. Such body waves allow us to focus on small velocity perturbations in the crust with high spatial resolution. We report on 10 years of seismic velocity temporal changes at the San Jacinto Fault. We observe and map a two‐month‐long episode of velocity changes with complex spatial distribution and interpret the velocity perturbation as produced by a previously undocumented slow‐slip event. We verify the hypothesis through numerical simulations and locate this event along a fault segment believed to be locked. Such a slow‐slip event stresses its surroundings and may trigger a major earthquake on a fault section approaching failure.

58 GEOSCIENCES↗

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↗

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↗

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↗

Evolution of Pore Structure and Permeability of Rocks Under Hydrothermal Conditions (Final Report)

The physical and transport properties of porous rocks can be altered by a variety of diagenetic, metamorphic, and tectonic processes, and the changes that result are of critical importance to such industrial applications as resource recovery, carbon dioxide sequestration, and waste isolation in geologic formations. The associated interrelationships between rocks, pore fluids and deformation are also key to understanding many natural processes, including dynamic metamorphism, fault mechanics, fault stability and pressure solution creep. In this project, we investigated the changes of permeability and pore geometry owing to inelastic deformation by solution-transfer, brittle fracturing and dislocation creep. In particular, we studied the coupling between pore fluids and deformation in fluid-filled quartz, calcite, mudstones and ultramafic rocks and examine the effects of coupled mechanical and chemical processes on the evolution of porosity and permeability under hydrothermal conditions. The investigations used a combination of techniques, that included triaxial laboratory experiments; uniaxial compressive loading with in situ observations of microstructure; permeability under hydrostatic compression; observations of microstructure using optical and electron microscopes, and micro-CT imaging; and numerical calculations. Laboratory experiments were designed to provide mechanical and transport data under conditions that isolate the particular mechanisms responsible for the changes. The data obtained were used to quantify changes in surface roughness, porosity, pore dimensions and their spatial fluctuations. The results of the experiments and data from image analyses were compared to the results of network, finite-difference and other numerical models to verify the validity of experimentally established relations between permeability and other rock properties. A bibliography for the period from 2015-2019 is appended below. New results obtained in the final year and the unfunded extension period are detailed below.

58 GEOSCIENCES↗

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↗

An empirical calibration of the serpentine-water oxygen isotope fractionation at T = 25–100°C

The ability to precisely constrain the physical–chemical conditions of serpentinization, such as temperature and fluid source, is limited by the accuracy of calibrations for oxygen isotope fractionation between serpentine and water – i.e., 1000 lnα( serpentine-water ) – which disagree by up to 20‰ when extrapolated to T<200°C. In this study, we present a new empirical calibration of 1000 lnα( serpentine-water ) aiming to improve applications of oxygen isotope geochemistry to very low-T serpentinization (T<100°C). We used the high-spatial resolution capabilities of Secondary Ion Mass Spectrometry (SIMS) to analyze oxygen isotope ratios in mineral pairs of calcite + serpentine, quartz + serpentine and talc + serpentine co-crystallized at scales≤50μm in six serpentinite samples from the Samail ophiolite (Oman). SIMS analysis shows that the mineral pairs are relatively homogenous in oxygen isotope ratios with variability in δ 18 O values ≤2‰ (2s). Clumped isotope thermometry and petrological constraints indicate crystallization temperatures from 25 to 100°C for the investigated samples. These independent constraints on temperature allowed us to derive 1000 lnα( serpentine-water ) by combining Δ 18 O( quartz-serpentine ), Δ 18 O( talc-serpentine ) and Δ 18 O( calcite-serpentine ) measured by SIMS in the investigated samples with Δ 18 O( quartz-water ), Δ 18 O( talc-water ) and Δ 18 O( calcite-water ) calibrations available in the literature. Our empirical calibration of 1000 lnα( serpentine-water ) = 1.04±0.20 (2SE) × 10 6 /T 2 (T in K), from T = 25 to 100°C, is within uncertainty of former high-temperature empirical calibrations extrapolated to T<100°C and matches experimental calibrations when extrapolated to T=250°C. In conclusion, the new serpentine-water calibration enables more accurate reconstructions of fluid-rock interactions occurring during low-temperature serpentinization processes in various tectonic settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Mountain biodiversity and ecosystem functions: interplay between geology and contemporary environments

Although biodiversity and ecosystem functions are strongly shaped by contemporary environments, such as climate and local biotic and abiotic attributes, relatively little is known about how they depend on long-term geological processes. Here, along a 3000-m elevational gradient with tectonic faults on the Tibetan Plateau (that is, Galongla Mountain in Medog County, China), we study the joint effects of geological and contemporary environments on biological communities, such as the diversity and community composition of plants and soil bacteria, and ecosystem functions. We find that these biological communities and ecosystem functions generally show consistent elevational breakpoints at 2000–2800 m, which coincide with Indus-Yalu suture zone fault and are similar to the elevational breakpoints of soil bacteria on another mountain range 1000 km away. Mean annual temperature, soil pH and moisture are the primary contemporary determinants of biodiversity and ecosystem functions, which support previous findings. However, compared with the models excluding geological processes, inclusion of geological effects, such as parent rock and weathering, increases 67.9 and 35.9% of the explained variations in plant and bacterial communities, respectively. Such inclusion increases 27.6% of the explained variations in ecosystem functions. The geological processes thus provide additional links to ecosystem properties, which are prominent but show divergent effects on biodiversity and ecosystem functions: parent rock and weathering exert considerable direct effects on biodiversity, whereas indirectly influence ecosystem functions via interactions with biodiversity and contemporary environments. Thus, the integration of geological processes with environmental gradients could enhance our understanding of biodiversity and, ultimately, ecosystem functioning across different climatic zones.

54 ENVIRONMENTAL SCIENCES↗

Lattice modeling of excavation damage in argillaceous clay formations: Influence of deformation and strength anisotropy

This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes.

58 GEOSCIENCES↗

Controlling Factors Analyses of Sediment Dispersal From Source-to-Sink Systems of Lacustrine Rift Basin: A Case Study From Paleogene Shahejie-Dongying Formations of Nanpu Sag, Eastern China

Tectonic activity not only shapes the basic stratigraphic framework of rift basins, but also profoundly affects the sediment dispersal in rift basins. In this study, analyses of heavy mineral assemblages in different periods demonstrate that there are three obvious tectono-sedimentary evolutionary stages (E s 3 –E s 2 , E s 1 –E d 2 , and E d 1 , respectively) in the Paleogene provenance area of Nanpu Sag, and the volume of sand bodies increases from the bottom of the Paleogene Shahejie (E s ) Formation to the top of the Dongying (E d ) Formation in Nanpu Sag. Besides, this study comprehensively utilize the analyses of seismic interpretation, palynology, heavy mineral, and borehole core samples to investigate the controlling factors of sediment dispersal in the rift basin. The assemblages of heavy minerals in different periods reflect the rock composition and catchment area of different provenance areas, and their vertical differences reflect the evolution process of the provenance area and reflect the uplift-denudation process of the provenance area. The results reveal that the synergy of the evolution of tectonic activity and the adjustment of topographical evolution are the main controlling factors of sediment dispersal in Nanpu Sag, while climate change is not the main controlling factor. We conclude that an increased sediment supply rate in the long term reflects the control of tectonic activity on basin topography, rather than climate fluctuations. The differences in morphological modification result in differential sediment dispersal, which is principally related to the differential extrusion of the fault system. The catchment area and provenance distance adjustment is evidenced by the vertical changes of heavy mineral characteristics of single-well and interaction and linkage of boundary faults, and the adjustment of topography evolution. A consideration is that the interaction and linkage of boundary faults and complex subsidence history are multi-directional, and differential evolution of provenance area is universal in lacustrine rift basins, all of this highlights the adjustment of sediment pathways generated by this characteristic of rift basins and emphasizes the importance of controlling factors analyses in understanding differential sediment dispersal that presents in the rift basins. Besides, four sets of sediment dispersal patterns were delineated based on different developmental regions in the rift basin, which are fault segmental point and multi-stage fault terrace, single-stage fault terrace and axial fault valley, axial fault terrace, and paleo-terrace and axial fault valley, respectively. This study has a certain guiding significance for the prediction of the spatial distribution of sand bodies in the rift basin and the exploration of potential oil and gas targets in the rift basin.

Yu, Jifeng↗

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↗

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↗

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↗