Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tectonics”

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

Shaping Peptide Assemblies Using Multifaceted Cyclic Tectons

Constructing distinct biomacromolecular assemblies typically necessitates target-specific selection and engineering of building blocks alongside optimization of assembly conditions. The challenge lies in achieving diverse morphological outcomes using simple, shared modules under identical conditions, a hallmark of natural systems that remains elusive in synthetic approaches. Here, in this work, we present a molecular scaffold-based strategy to instruct the coassembly of the same set of peptides into a variety of nanostructures across multiple dimensions. We create trifaceted cyclic scaffolds to manipulate two pairs of dimeric coiled-coil peptides prior to coassembly. These scaffolds, with addressable and orthogonal modules, allow controlled exposure of their cohesive faces, directing the formation of nanotriangles and fibrillar and lamellar assemblies. By tuning interpeptide arrangements that dictate scaffold geometry, we construct nonstraight fibrils with tunable curvature, which are rarely observed before. Notably, these scaffolds exhibit plasticity in adapting the sizes and orientations of cohesive faces to different assembly morphologies. The resultant nanostructures are consistent with the design and simulation results, demonstrating the reliability and predictability of this approach. Multifaceted cyclic scaffolds bridge the intellectual and physical gaps between building peptides and assemblies, holding promise for endowing various existing assembly systems with high tunability and versatility.

assembly morphology↗
DOE OSTI2020

Nowcasting Earthquakes in Southern California With Machine Learning: Bursts, Swarms, and Aftershocks May Be Related to Levels of Regional Tectonic Stress

Abstract Seismic bursts in Southern California are sequences of small earthquakes strongly clustered in space and time and include seismic swarms and aftershock sequences. A readily observable property of these events, the radius of gyration ( R G ), allows us to connect the bursts to the temporal occurrence of the largest M ≥ 7 earthquakes in California since 1984. In the Southern California earthquake catalog, we identify hundreds of these potentially coherent space‐time structures in a region defined by a circle of radius 600 km around Los Angeles. We compute R G for each cluster then filter them to identify those bursts with large numbers of events closely clustered in space, which we call “compact” bursts. Our basic assumption is that these compact bursts reflect the dynamics associated with large earthquakes. Once we have filtered the burst catalog, we apply an exponential moving average to construct a time series for the Southern California region. We observe that the R G of these bursts systematically decreases prior to large earthquakes, in a process that we might term “radial localization.” The R G then rapidly increases during an aftershock sequence, and a new cycle of “radial localization” then begins. These time series display cycles of recharge and discharge reminiscent of seismic stress accumulation and release in the elastic rebound process. The complex burst dynamics we observe are evidently a property of the region as a whole, rather than being associated with individual faults. This new method allows us to improve earthquake nowcasting, which is a technique to evaluate the current state of hazard in a seismically active region.

58 GEOSCIENCES↗
DOE OSTI2020

Heat Flow on the U.S. Beaufort Margin, Arctic Ocean: Implications for Ocean Warming, Methane Hydrate Stability, and Regional Tectonics

The Arctic Ocean represents one of the last geological frontiers on Earth. The formation of the Arctic Ocean remains unclear, and although it is well understood that the Arctic surface temperature is warming at a rate approximately twice as fast as the rest of the Earth, it is unclear how deep Arctic Ocean temperatures are changing. Understanding whether deep Arctic Ocean water is warming is important, since it can lead to the breakdown and release of huge quantities of frozen methane molecules trapped below the Arctic seafloor, which can increase ocean acidity and destabilize the seafloor. To understand both Arctic Ocean temperature change and geologic evolution, we collected temperature measurements at 103 sites. These measurements tell us how temperature increases with depth below the seafloor and can be used to understand both ocean temperature change and regional geology. Analysis of our data shows that at depths of 300–900 m below sea level, the Arctic Ocean has been warming steadily for perhaps several decades—nearly twice as deep as previous studies suggest. At ocean depths greater than 1,000 m, our analysis also reveals surprisingly high temperature increases with depth in the seafloor. The cause of these significant increases is unclear.

58 GEOSCIENCES↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2023

Evaluating the Aftershock Duration of Induced Earthquakes

As the carbon sequestration community prepares to scale up the number and size of commercial operations, the need for tools and methods to assess and mitigate risks associated with these operations becomes increasingly important. One outstanding question is whether aftershocks of induced events decay quickly after injection operations cease or if aftershock activity persists for hundreds of years before returning to background levels more akin to tectonic events (Stein and Liu, 2009). Appropriate estimates of the aftershock duration impact several operational management decisions including mitigation strategies and post-injection monitoring for seismic activity. It is hypothesized that induced earthquake rates may diminish more quickly after injection is stopped, owing to higher stressing rates from injected fluids. Alternatively, it is plausible that only the first event in the sequence is induced by increased fluid overpressures, whereas subsequent events (e.g., aftershocks) respond to the stored tectonic stresses and static and dynamic stress changes due to the mainshock (Keranen et al., 2013). Further, if the aftershock duration can be linked to stressing rates due to injection, then it follows that operational strategies to reduce seismic hazard by reducing injection rates or volumes may be successful. However, if aftershocks of induced events are relieving stored tectonic stresses, then altering injection volumes may not alleviate ongoing seismic activity. Furthermore, knowledge of an aftershock duration could aid in the determination of an appropriate post-injection monitoring period for ongoing seismicity, which is a factor in overall operational costs. In this study, we model induced seismicity sequences in Oklahoma with a coupled Coulomb rate–state earthquake rate model (Dieterich, 1994; Kroll et al., 2017) to estimate aftershocks durations. Results for the current study indicate that elevated rates of aftershock activity following induced mainshocks return to background seismicity rates in less than five years, contrary to the tens to hundreds of years observed for tectonic aftershocks.

58 GEOSCIENCES↗
DOE OSTI2023

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗
DOE OSTI2025

Climatic controls on erosional efficiency vary with lithology across the Himalaya

Although climate can strongly influence erosional efficiency (i.e., erosion rate for a given topography), demonstrating its impact in tectonically active areas has been challenging due to other confounding controlling factors, such as lithology. Here, we show that 10 Be-derived erosion rates and efficiencies in the Himalayan orogen exhibit distinct relationships with climatic factors depending on lithology. We compile 173 10 Be-derived, basin-averaged erosion rates across the orogen, including 12 newly measured rates from the Dibang and Lohit valleys in the easternmost Himalaya, regions characterized by high precipitation magnitudes and variability. We group basins based on lithologies separated by orogen-scale thrust faults and quantify erosional efficiency coefficients based on the relationships between erosion rates and topographic metrics. Our results show that erosion rates and erosional efficiency from sedimentary and metasedimentary rocks along the Himalayan range front display a positive, nonlinear correlation with climatic factors, such as the number of extreme rainfall events and mean annual precipitation rates. In contrast, erosion rates from crystalline lithologies in the hanging wall of the Main Central thrust show a strong correlation with fluvial topography, whereas erosional efficiency shows no statistically significant correlation with climatic factors. Rapid erosion rates and high erosional efficiencies in the eastern Himalayan range front are likely driven by extreme precipitation on tectonically active, steep slopes composed of mechanically weak metasedimentary rocks. Our findings highlight the importance of the interplay between controlling factors, which include tectonics, lithology, and climate, that drive surface erosion and influence the topographic evolution of orogenic systems.

Channel steepness↗
DOE OSTI2021

Improved lithospheric attenuation structure of the Arabian Peninsula through the use of national network data

We characterize the attenuation structure of the Arabian Peninsula through the measurement of regional phase amplitudes. High-resolution is achieved by combining stations from global networks with national network data through the cooperative effort of several countries in the region, including Saudi Arabia, Oman, Iraq, and Kuwait. The result is an improved attenuation model of the crust and upper mantle for a broad frequency band that extends from 0.5 to 10 Hz. The observed attenuation is in accordance with various elements of earth structure, including plate boundary type, style of tectonism, thermo-tectonic age, and temperature. Emerging features from the model include details in the structure along the Red Sea, and improved imaging of the southern Arabian Peninsula extending north from the Gulf of Aden. Finally, the resulting attenuation model can be employed for better magnitude estimates, in isolating tectonic and structural features, and in characterizing strong ground motion in the Arabian Peninsula.

58 GEOSCIENCES↗
DOE OSTI2022

Subduction, Underplating, and Return Flow Recorded in the Cycladic Blueschist Unit Exposed on Syros, Greece

Exhumed high-pressure/low-temperature (HP/LT) metamorphic rocks provide insights into deep (~20–70 km) subduction interface dynamics. On Syros Island (Cyclades, Greece), the Cycladic Blueschist Unit preserves blueschist-to-eclogite facies oceanic- and continental-affinity rocks that record the structural and thermal evolution linked to Eocene subduction. Despite decades of research, the metamorphic and deformation history (P-T-D) and timing of subduction and exhumation are matters of ongoing discussion. We suggest that Syros comprises three coherent tectonic slices and that each slice underwent subduction, underplating, and syn-subduction return flow along similar P-T trajectories, but at progressively younger times. Subduction and exhumation are distinguished by lineations and ductile fold axis orientations, and are kinematically consistent with previous studies that document top-to-the-S-SW shear (prograde-to-peak subduction), top-to-the-NE shear (blueschist facies exhumation), and then E-W coaxial stretching (greenschist facies exhumation). Amphibole zonations record cooling during decompression, indicating return flow above a cold slab. Multi-mineral Rb-Sr isochrons and compiled metamorphic geochronology show that the three slices record distinct stages of peak subduction (53–52, ~50, and 45 Ma) that young with structural depth. Retrograde blueschist and greenschist facies fabrics span ~50–40 and ~43–20 Ma, respectively, and also young with structural depth. Synthesized data sets support a revised tectonic framework for Syros, involving subduction of structurally distinct coherent slices and simultaneous return flow of previously accreted tectonic slices in the subduction channel shear zone. Distributed, ductile, dominantly coaxial return flow in an Eocene-Oligocene subduction channel proceeded at rates of ~1.5–5 mm/yr and accommodated ~80% of the total exhumation of this HP/LT complex.

58 GEOSCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2021

The Indian Ocean, its supra-subduction history, and implications for ophiolites

We report Ophiolite complexes represent fragments of ocean crust and mantle formed at spreading centers and emplaced on land. The setting of their origin, whether at mid-ocean ridges, back-arc basins, or forearc basins has been debated. Geochemical classification of many ophiolite extrusive rocks reflect an approach interpreting their tectonic environment as the same as rocks with similar compositions formed in various modern oceanic settings. This approach has pointed to the formation of many ophiolitic extrusive rocks in a supra-subduction zone (SSZ) environment. Paradoxically, structural and stratigraphic evidence suggests that many apparent SSZ-produced ophiolite complexes are more consistent with mid-ocean ridge settings. Compositions of lavas in the southeastern Indian Ocean resemble those of modern SSZ environments and SSZ ophiolites, although Indian Ocean lavas clearly formed in a mid-ocean ridge setting. These facts suggest that an interpretation of the tectonic environment of ophiolite formation based solely on their geochemistry may be unwarranted. New seismic images revealing extensive Mesozoic subduction zones beneath the southern Indian Ocean provide one mechanism to explain this apparent paradox. Cenozoic mid-ocean-ridge–derived ocean floor throughout the southern Indian Ocean apparently formed above former sites of subduction. Compositional remnants of previously subducted mantle in the upper mantle were involved in generation of mid-ocean ridge lavas. The concept of historical contingency may help resolve the ambiguity on understanding the environment of origin of ophiolites. Many ophiolites with “SSZ” compositions may have formed in a mid-ocean ridge setting such as the southeastern Indian Ocean.

58 GEOSCIENCES↗
DOE OSTI2021

Attention Network Forecasts Time-to-Failure in Laboratory Shear Experiments

Rocks under stress deform by creep mechanisms that include formation and slip on small-scale internal cracks. Intragranular cracks and slip along grain contacts release energy as elastic waves termed acoustic emissions (AE). AEs are thought to contain predictive information that can be used for fault failure forecasting. Here, we present a method using unsupervised classification and an attention network to forecast labquakes using AE waveform features. Our data were generated in a laboratory setting using a biaxial shearing device with granular fault gouge intended to mimic the conditions of tectonic faults. Here, we analyzed the temporal evolution of AEs generated throughout several hundred laboratory earthquake cycles. We used a Conscience Self-Organizing Map (CSOM) to perform topologically ordered vector quantization based on waveform properties. The resulting map was used to interactively cluster AEs. We examined the clusters over time to identify those with predictive ability. Finally, we used a variety of LSTM and attention-based networks to test the predictive power of the AE clusters. By tracking cumulative waveform features over the seismic cycle, the network is able to forecast the time-to-failure (TTF) of lab earthquakes. Our results show that analyzing the data to isolate predictive signals and using a more sophisticated network architecture are key to robustly forecasting labquakes. In the future, this method could be applied on tectonic faults to monitor earthquakes and augment early warning systems.

58 GEOSCIENCES↗
DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2025

The Mystery of Baselevel Controls in the Incision History of the Central Colorado Plateau

Erosion can remain active and changing in landscapes long after tectonic drivers have ceased, potentially due to local-geologic controls, climate changes, or geodynamics. We present new fluvial incision-rate histories and terrain analyses of the Colorado River system through the central Colorado Plateau to understand what has caused the variable erosion across this post-orogenic landscape. Results from new cosmogenic and luminescence dating of fluvial terrace and upland gravel deposits in Glen and Meander Canyons establish incision-rate histories that are marked by an Early-Middle Pleistocene erosion hiatus, followed by ~200 m of rapid incision over the last ~350 kyr. Projection of fluvial topography from above knickzones of the Colorado River drainage system roughly agree with the observed magnitude of recent incision and reflect a common baselevel fall from Pliocene river integration through Grand Canyon, which is still propagating through the drainage. A response-time model indicates that baselevel fall from integration likely took 2–4 Myr to reach the central Colorado Plateau and 100s kyr to travel across the study area, potentially accounting for incision rate changes in the fluvial terrace records of Meander and Glen Canyons. The upstream-migrating incision has likely been partitioned into multiple waves across the landscape due to the local geologic controls of lava damming, salt tectonics, and heterogenous bedrock. As baselevel fall from Pliocene Colorado River integration diffuses upstream, it can only account for perhaps a quarter of the total ~2 km of exhumation in the central Colorado Plateau, demanding an unknown driver for significant erosion in the Pliocene.

58 GEOSCIENCES↗
DOE OSTI2022

Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise

SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.

Silveira, Graça (ORCID:0000000221102554)↗
DOE OSTI2025

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗