Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Using Multiple Geophysical Methods to Refine a Stratigraphic Conceptual Site Model at a Nuclear Waste Site

At contaminated sites, there is a critical reliance on conceptual site models of which a key component is the geologic framework model (GFM), which is used to build predictive fate and transport modeling in support of remediation goals. Here, we detail an ongoing study at a nuclear waste site (Hanford Site, Washington, USA) where surface geophysical methods were used to complement existing borehole information and to help site new boreholes to update and refine the GFM. Ground-based electrical resistivity tomography (ERT), time-domain electromagnetics (TEM), and seismic data was acquired on the Hanford Site’s Central Plateau, where the presence of a hydraulically transmissive subsurface feature, or paleochannel, was suspected but uncertain based on contaminant concentrations at sparsely located boreholes. To assess and compare ERT, TEM, and seismic geophysical datasets, co-located data was acquired in an area where the existing GFM was more certain. Surrounding the suspected paleochannel, ERT and seismic refraction tomography showed similar subsurface structures consistent with a channelized feature. While TEM interpretation was limited in the suspected paleochannel area due to a thick resistive top layer and high electromagnetic noise, TEM showed greater utility in a different area where these factors were less prevalent. Based on the results of quasi-3D ERT inversions, we propose updates to the GFM, including identifying a paleochannel. Further, we evaluate the utility of geophysical methods and examine lessons learned that will be used as a basis for continued characterization efforts to support site management decisions and implementation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

River thorium concentrations can record bedrock fracture processes including some triggered by distant seismic events

Fractures are integral to the hydrology and geochemistry of watersheds, but our understanding of fracture dynamics is very limited because of the challenge of monitoring the subsurface. Here we provide evidence that long-term, high-frequency measurements of the river concentration of the ultra-trace element thorium (Th) can provide a signature of bedrock fracture processes spanning neighboring watersheds in Colorado. River Th concentrations show abrupt (subdaily) excursions and biexponential decay with approximately 1-day and 1-week time constants, concentration patterns that are distinct from all other solutes except beryllium and arsenic. The patterns are uncorrelated with daily precipitation records or seasonal trends in atmospheric deposition. Groundwater Th analyses are consistent with bedrock release and dilution upon mixing with river water. Most Th excursions have no seismic signatures that are detectable 50 km from the site, suggesting the Th concentrations can reveal aseismic fracture or fault events. We find, however, a weak statistical correlation between Th and seismic motion caused by distant earthquakes, possibly the first chemical signature of dynamic earthquake triggering, a phenomenon previously identified only through geophysical methods.

58 GEOSCIENCES↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simultaneous inference of plate boundary stresses and mantle rheology using adjoints: large-scale 2-D models

SUMMARY Plate motions are a primary surface constraint on plate and mantle dynamics and rheology, plate boundary stresses and the occurrence of great earthquakes. Within an optimization method, we use plate motion data to better constrain uncertain mantle parameters. For the optimization problem characterizing the maximum a posteriori rheological parameters we derive gradients using adjoints and expressions to approximate the posterior distributions for stresses within plate boundaries. We apply these methods to a 2-D cross section from the western to eastern Pacific, with temperature distributions and fault zone geometries developed primarily from seismic and plate motion data. We find that the best-fitting stress exponent, n, is about 2.8 and the yield stress about 100 MPa or less. The normal stress on the interplate fault zones is about 100 MPa and the shear stresses about 10 MPa or less.

Rudi, Johann (ORCID:0000000265639265)↗
DOE OSTI2021

Failure prediction and estimation of failure parameters

Machine-learning methods and apparatus are disclosed to determine frictional state or other parameters in an earthquake zone or other failing medium, using acoustic emission, seismic waves, or other detectable indicators of microscopic processes. Predictions of future failures are demonstrated in different regimes. A classifier is trained using time series of acoustic emission data along with historic data of frictional state or failure events. In disclosed examples, random forests and gradient boost trees are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to testing or live data in order to assess a frictional state, assess seismic hazard, or make predictions regarding a future failure event. The technology has been developed in a double direct shear apparatus, but can be widely applied to seismic faults, other terrestrial failures, or failures in man-made structures. Variations are disclosed.

Johnson, Paul Allan↗
DOE OSTI2021

Kimberlina 1.2 Velocity Models and Seismic Data

Kimberlina 1.2 Velocity model and synthetic seismic data, produced in collaboration of teams at the National Energy Technology Laboratory, Los Alamos National Laboratory, and Lawrence Livermore National Laboratory through the National Risk Assessment Partnership. Data is associated with the following publication: Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019. The Kimberlina 1.2 Velocity models were produced by Zan Wang, Robert Dilmore, William Harbert, and Lianjie Huang at NETL. The following citations are directly related to the creation of the velocity models: Wang, Z. Harbert, W., Dilmore, R., Huang, L. Modeling of time-lapse seismic monitoring using CO2 leakage simulations for a model CO2 storage site with realistic geology: Application in assessment of early leak-detection capabilities. International Journal of Greenhouse Gas Control. V. 76, September 2018, Pages 39-52. https://doi.org/10.1016/j.ijggc.2018.06.011 Wang, Z., Dilmore, R., Harbert, W. Inferring CO2 saturation from synthetic surface seismic and downhole monitoring data using machine learning for leakage detection at CO2 sequestration sites. International Journal of Greenhouse Gas Control, V. 100, September 2020. https://doi.org/10.1016/j.ijggc.2020.103115 The velocity models were built based on the Kimberlina 1.2 aquifer impact data which is associated with the following publications: Buscheck, T.A., Mansoor, K., Yang, X., Wainwright, H., and Carroll, S. (2019). Downhole pressure and chemical monitoring for CO2 and brine leak detection in aquifers above a CO2 storage reservoir. International Journal of Greenhouse Gas Control. 91. 102812. 10.1016/j.ijggc.2019.102812. Xianjin Yang, Thomas A. Buscheck, Kayyum Mansoor, Zan Wang, Kai Gao, Lianjie Huang, Delphine Appriou, Susan A. Carroll, Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, Volume 90, 2019, 102803, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2019.102803 The synthetic seismic data was produced by Youzuo Lin and team at LANL, and are associated with the following citations: Jordan, P. D., and J. L. Wagoner. Characterizing Construction of Existing Wells to a CO2 Storage Target: The Kimberlina Site, California. Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019.

CO2 Leakage↗
NASA NTRS2019

Combined Qualification Vibration Testing and Fixed Base Modal Testing Utilizing a Fixed Based Correction Method

Vibration testing spaceflight hardware is a vital, but time consuming and expensive endeavor. Traditionally modal tests are performed at the component, subassembly, or system level, preferably free-free with mass loaded interfaces or fixed base on a seismic mass to identify the fundamental structural dynamic (modal) characteristics. Vibration tests are then traditionally performed on single-axis slip tables at qualification levels that envelope the maximum predicted flight environment and workmanship in order to verify the spaceflight hardware can survive its flight environment. These two tests currently require two significantly different test setups, facilities, and ultimately reconfiguration of the spaceflight hardware. The vision of this research is to show how traditional fixed-base modal testing can be accomplished using vibration qualification testing facilities, which not only streamlines testing and reduces test costs, but also opens up the possibility of performing modal testing to untraditionally high excitation levels that provide for test-correlated finite element models to be more representative of the spaceflight hardware's response in a flight environment. This paper documents the first steps towards this vision, which is the comparison of modal parameters identified from a traditional fixed-based modal test performed on a modal floor and those obtained by utilizing a fixed based correction method with a large single-axis electrodynamic shaker driving a slip table supplemented with additional small portable shakers driving on the slip table and test article. To show robustness of this approach, the test article chosen is a simple linear weldment, whose mass, size, and modal parameters couple well with the dynamics of the shaker/slip table. This paper will show that all dynamics due to the shaker/slip table were successfully removed resulting in true fixed-base modal parameters, including modal damping, being successfully extracted from a traditional style base-shake vibration test setup.

Base-Shake↗
NASA NTRS2019

Combined Qualification Vibration Testing and Fixed Base Modal Testing Utilizing a Fixed Based Correction Method

Vibration testing spaceflight hardware is a vital, but time consuming and expensive endeavor. Traditionally modal tests are performed at the component, subassembly, or system level, preferably free-free with mass loaded interfaces or fixed base on a seismic mass to identify the fundamental structural dynamic (modal) characteristics. Vibration tests are then traditionally performed on single-axis slip tables at qualification levels that envelope the maximum predicted flight environment plus 3 dB and workmanship in order to verify the spaceflight hardware can survive its flight environment. These two tests currently require two significantly different test setups, facilities, and ultimately reconfiguration of the spaceflight hardware. The vision of this research is to show how traditional fixed-base modal testing can be accomplished using vibration qualification testing facilities, which not only streamlines testing and reduces test costs, but also opens up the possibility of performing modal testing to untraditionally high excitation levels that provide for test-correlated finite element models to be more representative of the spaceflight hardware's response in a flight environment. This paper documents the first steps towards this vision, which is the comparison of modal parameters identified from a traditional fixed-based modal test performed on a modal floor and those obtained by utilizing a fixed based correction method with a large single-axis electrodynamic shaker driving a slip table supplemented with additional small portable shakers driving on the slip table and test article. To show robustness of this approach, the test article chosen is a simple linear weldment, whose mass, size, and modal parameters couple well with the dynamics of the shaker/slip table. This paper will show that all dynamics due to the shaker/slip table were successfully removed resulting in true fixed-base modal parameters, including modal damping, being successfully extracted from a traditional style base-shake vibration test setup.

Winkel, James P.↗
NASA NTRS1975

A seismic refraction technique used for subsurface investigations at Meteor Crater, Arizona

A seismic refraction technique for interpreting the subsurface shape and velocity distribution of an anomalous surface feature such as an impact crater is described. The method requires the existence of a relatively deep refracting horizon and combines data obtained from both standard shallow refraction spreads and distant offset shots by using the deep refractor as a source of initial arrivals. Results obtained from applying the technique to Meteor crater generally agree with the known structure of the crater deduced by other investigators and provide new data on an extensive fractured zone surrounding the crater. The breccia lens is computed to extend roughly 190 m below the crater floor, about 30 m less than the value deduced from early drilling data. Rocks around the crater are fractured as distant as 900 m from the rim crest and to a depth of at least 800 m beneath the crater floor.

Ackermann, H. D.↗
NASA NTRS1992

Geophysical Methods: an Overview

Geophysics is expected to have a major role in lunar resource assessment when manned systems return to the Moon. Geophysical measurements made from a lunar rover will contribute to a number of key studies: estimating regolith thickness, detection of possible large-diameter lava tubes within maria basalts, detection of possible subsurface ice in polar regions, detection of conductive minerals that formed directly from a melt (orthomagmatic sulfides of Cu, Ni, Co), and mapping lunar geology beneath the regolith. The techniques that can be used are dictated both by objectives and by our abilities to adapt current technology to lunar conditions. Instrument size, weight, power requirements, and freedom from orientation errors are factors we have considered. Among the geophysical methods we believe to be appropriate for a lunar resource assessment are magnetics, including gradiometry, time-domain magnetic induction, ground-penetrating radar, seismic reflection, and gravimetry.

Becker, A.↗
DOE OSTI2020

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗
DOE OSTI2025

STILGAR: Subsurface Models for Graymont Pleasant Gap Mine

The detection, location, and monitoring of underground structures are of great importance to national and global security. Tunnels and voids generate seismic signatures detectable at the surface, but using non-invasive seismic data to image near-surface presents several challenges in real-world applications. In this report, we describe the use of a dense surface seismic deployment to generate subsurface models of the Graymont Pleasant Gap mine - a single-layer mine with a complex structure embedded in a high-velocity P-wave limestone bedrock. Our approach consists of three key methods. We use P-wave arrival times from local blast events to perform a tomography inversion with the tomoTD method, constructing a P-wave velocity model of the subsurface. We model the layer above the mine using Rayleigh wave ellipticity and inversion techniques. We leverage ongoing anthropogenic activities to identify and locate noise sources both on the surface and within the subsurface. With this integrated approach we aim to overcome the challenges and enhance our ability to non-invasively characterize underground structures, contributing to improved seismic monitoring techniques.

58 GEOSCIENCES↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

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↗
NASA NTRS2020

The Influence of Frequency and Duration Of Seismic Ground Motion On The Size Of Triggered Landslides—A Regional View

Observation, theory, and intuition all suggest that larger earthquakes should trigger larger landslides. Many factors could contribute to this, including depth-dependent shear strength or non-linearity of ground motion in soils and rock, but we hypothesize that the key characteristics of large earthquakes causing this phenomenon are (in addition to magnitude) the frequency and duration of the strong ground motion. Because of the paucity of site-specific data for detailed analysis, we take a regional approach to this question by analyzing strong-motion records and earthquake-induced landslide (EQIL) inventories from six well-documented earthquakes. Ground motion is characterized using earthquake magnitude and the median durations and frequencies (mean periods) of subsets of strong-motion records relevant to landslide triggering. EQIL inventories are characterized using the median landslide area of the entire inventory as well as the median areas of the largest 1% of the landslides and the largest 10 landslides. We then compare ground-motion characteristics with landslide size statistics to determine possible correlations. Comparisons of all earthquake- and landslide-size statistics show strong positive correlations between landslide size and (1) magnitude, (2) ground-motion duration, and (3) mean period. Although all the ground-motion measures yield highly correlated regressions, mean period appears to be the best overall predictor of landslide size. Landslide modeling using Newmark's sliding-block method also shows that longer mean periods and durations and larger magnitudes correlate strongly with increases in modeled displacements. These results support our hypothesis that increasing period and duration of seismic ground motion are the physical factors driving increased landslide sizes for larger earthquakes. Additional studies including data from a much larger set of earthquakes is needed to confirm the results of this initial study.

Hakan Tanyas↗
DOE OSTI2022

Maximizing the Proppant Carrying and Viscoelastic Properties of the Bakken Hypersaline-Produced Water with High-Viscosity Friction Reducers for Sustainable Applications

Summary The development and production of unconventional reservoirs, such as the Bakken Formation, have become a resolved mystery for operators in North America since the arrival and advancement of horizontal drilling and hydraulic fracturing technologies. As a result, unconventional reservoir assets became the central focus of the oil and gas industry at the state, national, and global levels. The produced water from these activities in the Bakken Formation have high salt contents (110,000–350,000 ppm) total dissolve solids (TDS) and can pose significant challenges to the environment if not treated. Deep injection into disposal wells is the routine method used to get rid of the Bakken produced water. However, there have been some concerns that unrestrained injections, in addition to polluting the groundwater, could potentially lead to seismic activities either at the time of injection or in the near future. To diminish the environmental impacts that may be associated with induced seismicity, including the reduction of the costs of water acquisition, the produced water can be treated and reused in the hydraulic fracturing processes. Also, the treated water could be used for irrigation purposes, for power generation, and coal mining operations. The issues of waste water and residual oil high in TDS are challenges yet to be effectively addressed despite preceding research and studies on advancing produced water technologies. The goal of this study is to explore all applicable ways by which the produced water from the Bakken Formation can serve as a replacement base fluid for use with polymers like the high-viscosity friction reducers (HVFRs) to create hydraulic fracturing fluids that can be stable with reservoir conditions and also be able to minimize environmental impacts and cost of operations. Experimental investigations using the high-salinity produced water from the Bakken Formation with HVFRs were carried out. The studies included a base case that served as a bench mark for comparing the effectiveness of the other scenarios. The results indicate that the Bakken hyper-saline produced water can withstand effect of heavy metals, salinity, hardness and remain stable through different shear rates (66–330 s−1) when treated with higher dosages [4–8 gal/1,000 gal (gpt)] of HVFRs. Filtration and dilution were the only methods used on the Bakken Formation produced wate for this research

Engineering↗
DOE OSTI2022

Ambient Noise Attenuation Tomography Reveals an Asymmetric Damage Zone Across San Jacinto Fault Near Anza, California

Here, we perform seismic attenuation tomography of the shallow structure for the San Jacinto Fault in the Ramona Reservation of southern California. The study uses ambient seismic noise recorded by a linear array of 65 3-C sensors across the fault. We extract amplitude decay information, with uncertainty, from noise interferometry functions. To account for strong heterogeneities in the complex shallow fault zone structure, we apply a frequency-dependent amplitude correction for focusing/defocusing effects using inverted phase velocity maps obtained by solving the transport equation. We then estimate the attenuation structure based on the linear station-triplet method for both Love and Rayleigh waves. The attenuation tomography indicates strong attenuation correlated with known San Jacinto Fault surface traces. The Love wave attenuation tomography reveals an asymmetric damage zone that exists primarily on the fault side with faster seismic velocity, consistent with earthquake ruptures on a fault bimaterial interface with preferred propagation direction to the northwest.

58 GEOSCIENCES↗