Search NASA⌕ Search

SEARCH · Search NASA

Results for “moment tensor”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evolution of tensile and shear cracking in crystalline rocks under compression

Here we performed laboratory-scale unconfined compression experiments on Barre granite specimen with a single pre-existing flaw to study the evolution of tensile and shear cracking that occurs during different stress levels. Moment tensor inversion of acoustic emission (AE) events in combination with the nonelastic strain component obtained through the 2D digital image correlation (DIC) technique are used to track the changes in the source mechanisms of the stress-induced cracks. The mode of deformation computed from the image based strain profiles, enables visual comparison of the nucleation, growth and interaction of the tensile and shear cracks with the microseismic source mechanism observed by moment tensor inversion of the AE events. Comparing the results obtained from the two techniques on the same dataset demonstrates a quantitative correlation between the acoustic and visual observations in terms of the cracking mechanisms obtained at different stages of cracking from crack initiation at the flaw tips to the failure of the rock specimen.

58 GEOSCIENCES↗
DOE OSTI2021

Systematic Generation of Moment Invariant Bases for 2D and 3D Tensor Fields

Moment invariants have been successfully applied to pattern detection tasks in 2D and 3D scalar, vector, and matrix valued data. However so far no flexible basis of invariants exists, i.e., no set that is optimal in the sense that it is complete and independent for every input pattern. Here, we prove that a basis of moment invariants can be generated that consists of tensor contractions of not more than two different moment tensors each under the conjecture of the set of all possible tensor contractions to be complete. This result allows us to derive the first generator algorithm that produces flexible bases of moment invariants with respect to orthogonal transformations by selecting a single non-zero moment to pair with all others in these two-factor products. Since at least one non-zero moment can be found in every non-zero pattern, this approach always generates a complete set of descriptors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2020

Directional Cooling-Induced Fracturing Westerly Granite Test Results

Directional Cooling-Induced Fracturing (DCIF) experiments were conducted on a short, cylindrical sample of Westerly granite (diameter = 4 inches, height ~ 2 inches). Liquid nitrogen was poured in a copper cup attached to the top of the sample, and the resulting acoustic emissions (AEs) and temperature changes on the surface of the sample were monitored. The obtained AEs were used to determine the microcracking source locations and amplitude, and the associated moment tensors. Included in this submission is an animation of the AEs, a graphic displaying the temperature changes, and the measured data.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Analysis of the 2021 Milford, Utah earthquake swarm: Enhanced earthquake catalog and migration patterns

An earthquake swarm occurred in Spring 2021 in south-central Utah near the town of Milford. The University of Utah Seismograph Stations located 125 earthquakes between March 19 and May 10 with magnitudes ranging from 0.5 to 3.2. We implement a matched-filter technique in order to identify additional earthquakes that went undetected during the routine network location. The 125 network-located earthquakes are used as templates and are cross-correlated with continuous data for the dates Feb 17—June 10. This time period corresponds to approximately 1 month before the earthquake swarm began through 1 month after it ended. For the matched-filter analysis, we rely heavily on station FOR1, which is located within 5 km of most template events. Four other stations within 20–30 km of most template events provide a supplement to the closest station. The matched-filter implementation results in the detection of over 600 earthquakes in addition to the original 125 catalog events. This is one of the largest swarms ever recorded in Utah, and no previous large swarms have been recorded in this location. We use HypoDD to obtain relative double-difference locations of the catalog events. Both routine locations and HypoDD relocations of the catalog events suggest a fault dipping west beneath the mountain range, opposite of typical Basin and Range normal faults that dip beneath alluvial valleys. Moment tensors for the largest five events show normal faulting consistent with the west-dipping fault seen in the seismicity. Hydrothermal features in the area, including a geothermal power plant, suggest that fluids may be a contributing factor to the earthquake swarm triggering. We examine the role of fluids by exploring bounds on diffusion parameters and investigating spatial migration characteristics of the swarm seismicity. We conclude that this swarm is the result of heterogenous stress conditions in a prefractured region.

58 GEOSCIENCES↗
DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Closure theory for high-collisionality multi-ion plasmas

A general formalism is developed to construct and solve a system of linearized moment equations for parallel and perpendicular closures in high-collisionality plasmas. It is applicable for multiple ion species with arbitrary masses, temperatures, charges, and densities. The convergence of closure coefficients is evaluated by increasing the number of moments from 2 to 32 for scalar, vector, and rank-2 tensor moments. As an example, the complete set of closure coefficients for a deuterium-carbon plasma over the entire Hall parameter range is presented. Furthermore, the closure coefficients at various temperature ratios show that the one-temperature closure coefficients can differ significantly from the two-temperature coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Pn magnitude scale m b ( Pn ) for earthquakes along the equatorial Mid-Atlantic Ridge

We developed a short-period Pn magnitude scale m b (Pn) for earthquakes along the equatorial Mid-Atlantic Ridge. Due to low signal-to-noise ratios, teleseismic body wave magnitude and long-period surface wave magnitude cannot be confidently determined for small earthquakes of m b < 4. Local magnitude scales are also not useful for these events because the oceanic environment does not allow the propagation of crustal phases. However, regional high-frequency Pn waves from these small- to moderate-size (m b 3–6) earthquakes are well recorded in the equatorial Atlantic region and can be used to assign magnitudes. We measured over 2041 Pn peak amplitudes on vertical records from about 21 stations in northeastern Brazil and 11 stations in western Africa in the distance range of 700–3700 km. We analysed data from 189 events from the global centroid moment tensor catalogue to tie our m b (Pn) scale to M w so that seismic moments can be readily estimated. Pn arrivals show apparent group velocity between 7.9 km s –1 at short ranges (~1000 km) and up to 9.1 km s –1 at 3500 km. The measured peak amplitudes have a frequency between 0.8 and 3 Hz at 1000–1800 km, but at greater distances, 1800–3700 km, they show a remarkably consistent frequency of about 0.8 Hz. The peak amplitude attenuates at a higher rate at short distances (~0.65 magnitude units between 700 and 2000 km) but attenuates at a lower rate at long distances (~0.35 magnitude units between 2000 and 3700 km). The low rate of amplitude decay with distance and nearly constant frequency content of the peak amplitudes suggest that Pn waves propagate efficiently in the lower part of the upper mantle in the equatorial Atlantic Ocean basins. Furthermore, these are important attributes of oceanic Pn waves that can be used to assign magnitude for small- to moderate-size earthquakes in the equatorial mid-Atlantic region. The estimated station corrections correlate well with upper mantle low-velocity anomalies, especially in Brazil.

58 GEOSCIENCES↗
NASA NTRS2013

Source Parameter Inversion for Recent Great Earthquakes from a Decade-long Observation of Global Gravity Fields

We quantify gravity changes after great earthquakes present within the 10 year long time series of monthly Gravity Recovery and Climate Experiment (GRACE) gravity fields. Using spherical harmonic normal-mode formulation, the respective source parameters of moment tensor and double-couple were estimated. For the 2004 Sumatra-Andaman earthquake, the gravity data indicate a composite moment of 1.2x10(exp 23)Nm with a dip of 10deg, in agreement with the estimate obtained at ultralong seismic periods. For the 2010 Maule earthquake, the GRACE solutions range from 2.0 to 2.7x10(exp 22)Nm for dips of 12deg-24deg and centroid depths within the lower crust. For the 2011 Tohoku-Oki earthquake, the estimated scalar moments range from 4.1 to 6.1x10(exp 22)Nm, with dips of 9deg-19deg and centroid depths within the lower crust. For the 2012 Indian Ocean strike-slip earthquakes, the gravity data delineate a composite moment of 1.9x10(exp 22)Nm regardless of the centroid depth, comparing favorably with the total moment of the main ruptures and aftershocks. The smallest event we successfully analyzed with GRACE was the 2007 Bengkulu earthquake with M(sub 0) approx. 5.0x10(exp 21)Nm. We found that the gravity data constrain the focal mechanism with the centroid only within the upper and lower crustal layers for thrust events. Deeper sources (i.e., in the upper mantle) could not reproduce the gravity observation as the larger rigidity and bulk modulus at mantle depths inhibit the interior from changing its volume, thus reducing the negative gravity component. Focal mechanisms and seismic moments obtained in this study represent the behavior of the sources on temporal and spatial scales exceeding the seismic and geodetic spectrum.

Parameter Inversion↗
NASA NTRS1994

The effects of gravity darkening on the ultraviolet continuum polarization produced by circumstellar disks

We investigate the effects of gravity darkening on the UV continuum polarization produced by an axisymmetric disk that surrounds a rapidly rotating star. Although the model is a single scattering approximation, we do include the effects of attenuation (electron scattering plus hydrogen bound-free absorption) by the disk, using an approach similar to that of Sobolev (1963). Because of the gravity darkening of the star and the attenuation within the disk, the radiation field is not axially symmetric about the radius vector. This implies that the polarization source functions are no longer provided by the finite disk depolarization factors of Cassinelli, Nordsieck, & Murison (1987), which are functions of the intensity moments in a spherically symmetric atmosphere. We reformulate the polarization source functions using generalized intensity moment tensors (J, H(sub i), K(sub ij)) that are valid for an arbitrary radiation field and envelope geometry. We find that the polarization source functions are simplest when using intensity moments in the observer's reference frame. On the other hand, the intensity moments are most easily evaluated in the stellar reference frame. Using the rotation transformation properties of the generalized intensity moments, we relate the observer's moments to those evaluated in the stellar reference frame. Our procedure for determining the polarization source functions thus merely involves choosing a set of Euler angles for the coordinate rotations, and then evaluating the associated rotation matrix. The geometrical complications of polarization transfer are thus reduced to obtaining a coordinate rotation matrix.

Bjorkman, J. E.↗
NASA NTRS1994

Dynamical influences on the moment of inertia tensor from lateral viscosity variations inferred from seismic tomographic models

We have investigated the influences of lateral variations of viscosity on the moment of inertia tensor from viscous flows due to the density anomalies in the mantle inferred from seismic tomographic models. The scaling relations between the density and the seismic anomalies is taken as either a constant or a function increasing with depth in accord with the recent high-pressure experimental studies. The viscosity is taken as an exponential function of the 3D density anomaly. In models with an isoviscous background, the effects on the perturbed moment of inertia tensor from the lateral viscosity variations are smaller than those due to variations in the radial viscosity profiles. In mantle models with a background viscosity increasing with depth, the influences of the lateral viscosity variations are significant. The most striking feature in the latter case is that the two off-diagonal elements delta I(sub xz) and delta I(sub yz) in the inertia tensor exhibit greatest sensitivity to lateral variations of the viscosity. While the other elements of the inertia change by only about a few tens of percent in the range of lateral viscosity contrast considered (less than 300), delta I(sub xz) and delta I(sub yz) can vary up to 40 times even with a change in sign, depending on the radial viscosity stratification and the location of the strongest lateral variations. The increase in the velocity-density scaling relation with depth can reduce the influences of the lateral viscosity variations, but it does not change the overall sensitive nature of delta I(sub xz) and delta I(sub yz). This study demonstrates clearly that the lateral viscosity variations, especially in the upper mantle, must be considered in the determination of long-term polar wander, since the variations in the delta I(sub xz) and delta I(sub yz) terms are directly responsible for exciting rotational movements.

Zhang, Shuxia↗
NASA NTRS2004

Co-Seismic Mass Displacement and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) displacements in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field. The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross. The calculation uses the normal mode summation scheme, applied to over twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Centroid Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies, conspiring to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards approx. 140 deg.E, roughly opposite to the observed polar drift direction. Currently, the Gravity Recovery And Climate Experiment (GRACE) is measuring the time-variable gravity to high degree and order with unprecedented accuracy. Our results show that great earthquakes such as the 1960 Chilean or 1964 Alaskan events cause gravitational field changes that are large enough to be detected by GRACE.

Chao, B. F.↗
NASA NTRS1999

Co-Seismic Mass Dislocation and Its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results. The calculation uses the normal mode summation scheme, applied to 15,814 major earthquakes that occurred during 1976-1998, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J(sub 2) and J(sub 22) while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards about 140 degree E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, Benjamin F.↗
DOE OSTI2024

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗
NASA NTRS2002

Co-Seismic Mass Dislocation and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the shaking that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross (1987). The calculation uses the normal mode summation scheme, applied to nearly twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to nudge the Earth rotation pole towards approximately 140 degrees E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, B. F.↗
NASA NTRS2000

Coseismic Excitation of the Earth's Polar Motion

Apart from the "shaking" near the epicenter that is the earthquake, a seismic event creates a permanent field of dislocation in the entire Earth. This redistribution of mass changes (slightly) the Earth's inertia tensor; and the Earth's rotation will change in accordance with the conservation of angular momentum. Similar to this seismic excitation of Earth rotation variations, the same mass redistribution causes (slight) changes in the Earth's gravitational field expressible in terms of changes in the Stokes coefficients of its harmonic expansion. In this paper, we give a historical background of the subject and discuss the related physics; we then compute the geodynamic effects caused by earthquakes based on normal-mode summation scheme. The effects are computed using the centroid moment tensor (CMT) solutions for 15,814 major earthquakes from Jan., 1977, through Feb., 1999, as provided in the Harvard CMT catalog. The computational results further strengthens these findings and conclusions: (i) the strong tendency for earthquakes to make the Earth rounder and more compact (however slightly) continues; (ii) so does the trend in the seismic "nudging" of the rotation pole toward the general direction of approx. 140 E, roughly opposite to that of the observed polar drift, but two orders of magnitude smaller in drift speed.

Chao, B. F.↗
DOE OSTI2021

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗
DOE OSTI2023

Microseismic Monitoring of Laboratory Hydraulic Fracturing Experiments in Granitic Rocks for Different Fracture Propagation Regimes

While hydraulic fracturing (HF) is a widely employed process, the underlying fracturing processes are still heavily contested. The attributes of the HF generated fracture network can exhibit substantial variation when dealing with specific HF propagation regimes encountered in the field. In this study, HF experiments were performed on true-triaxially loaded Barre granite cubes, with microseismic monitoring, to identify and characterize the fracturing mechanisms associated with different viscosity injection fluids. Utilizing fluids with high (oil/1450 cP) and low (water/1 cP) viscosity represented two key HF propagation regimes: viscosity- and toughness-dominated. The experiments conducted with oil involved higher breakdown pressures, larger fluid volumes, and slower fracture propagation speeds. Furthermore, the frequency–magnitude distribution (b value) for all experiments (1.9–2.3) was similar to those encountered for large-scale operations. Slightly larger b values were encountered during the initiation phase (2.4–2.7) relative to the fracture propagation and post-fracturing phases (1.9–2.2). Techniques such as polarity and moment-tensor inversion were utilized to characterize the source mechanisms. For the HF experiments with oil, tensile fractures were most dominant (92%) in the initiation phase compared to fracture propagation and post-fracturing phases (70–75%). Similar tensile fracturing dominancy was not observed with water, attributable to fluid permeation and leak-off. Regardless of the injection fluid or classification criteria employed, tensile fractures were the dominant type consistently, with fewer occurring in water experiments but the specific ratio of crack types varied with different source mechanism criteria employed.

58 GEOSCIENCES↗