Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave direction”

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 829 records · Page 46

Unitarity methods in AdS/CFT

We develop a systematic unitarity method for loop-level AdS scattering amplitudes, dual to non-planar CFT correlators, from both bulk and boundary perspectives. We identify cut operators acting on bulk amplitudes that put virtual lines on shell, and show how the conformal partial wave decomposition of the amplitudes may be efficiently computed by gluing lower-loop amplitudes. A central role is played by the double discontinuity of the amplitude, which has a direct relation to these cuts. We then exhibit a precise, intuitive map between the diagrammatic approach in the bulk using cutting and gluing, and the algebraic, holographic unitarity method of [1] that constructs the non-planar correlator from planar CFT data. Our analysis focuses mostly on four-point, one-loop diagrams — we compute cuts of the scalar bubble, triangle and box, as well as some one-particle reducible diagrams — in addition to the five-point tree and four-point double-ladder. Analogies with S-matrix unitarity methods are drawn throughout.

79 ASTRONOMY AND ASTROPHYSICS↗

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of hollow cathode plasma turbulence using coherent Thomson scattering

Turbulence in hollow cathodes used for space propulsion is believed to play an important role in anomalous electron transport and ion heating. In this work, the implementation of coherent Thomson scattering to identify and characterize MHz-frequency ion acoustic turbulence and kHz-frequency oscillations in the plume of a hollow cathode is achieved. In the presence of a background magnetic field of a Hall thruster, a number of unstable modes are observed. A directive ion acoustic mode propagating predominantly within a restricted angle around the magnetic field is found, exhibiting an energy scaling with wavenumber k of the form k -5.2±0.58 , which differs from the classic Kadomtsev k -3 scaling for unmagnetized conditions. Bi-directional ion acoustic mode fluctuations propagating over a range of angles with respect to the magnetic field have been measured, possibly signifying the existence of a large-amplitude plasma wave, similar to the Buneman instability. Finally, electron density fluctuations in the kHz-frequency range, a possible consequence of drift-driven instabilities in the plane perpendicular to the magnetic field, have also been identified. Furthermore, these results not only are an indication of the diversity of wave types that exist in hollow cathode plumes but also point to the key role played by the presence of, and the configuration of, the magnetic field in their appearance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accelerating astronomical and cosmological inference with preconditioned Monte Carlo

ABSTRACT We introduce preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilizes a Normalizing Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison, respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster, respectively, than nested sampling (NS). We found that in higher dimensional applications, the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimized low-depth quantum circuits for molecular electronic structure using a separable-pair approximation

We present a classically tractable model that leads to optimized low-depth quantum circuits leveraging separable-pair approximations. The obtained circuits are well suited as a baseline circuit for emerging quantum hardware and can, in the long term, provide significantly improved initial states for quantum algorithms. The associated wave functions can be represented with linear memory requirement, which allows classical optimization of the circuits and naturally defines a minimum benchmark for quantum algorithms. In this work we employ directly determined pair-natural orbitals within a basis-set-free approach. This leads to accurate representation of the one- and many-body parts for weakly correlated systems and we explicitly illustrate how the model can be integrated into other quantum algorithms for stronger correlated systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Effects of Landau quantization on neutrino emission and absorption

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as 10 15 G and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as 10 16 G. In the presence of a strong magnetic field, Leinson and Perez showed and Baiko and Yakovlev confirmed that the energy levels of electrons and protons are quantized and the direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields 𝐵 ≳ 10 16 G, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. Here, we present a fully relativistic calculation of the direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wave functions present and a new semianalytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at subsaturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

astrophysical electromagnetic fields↗

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗

Hybrid Simulation and Quasi-linear Theory of Bi-Kappa Proton Instabilities

The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropy T ⊥ > T ∥ (where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κ power-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Power Take-Off Design Study for a Small-Scale Oscillating Surge Wave Energy Converter for Powering the Blue Economy Applications

The power take-off (PTO) is an integral part of wave energy conversion, and the design process is nontrivial. Better PTOs, and better processes for selecting and designing PTO architectures for various applications, would benefit devices that assist in powering the blue economy by decreasing time and money spent on PTO design and increasing the overall energy capture performance of these devices. This paper chronicles the selection process of a PTO for a small-scale surge-type wave energy converter (WEC) for the purpose of informing future PTO selection processes. Three PTO architectures are evaluated in WEC-Sim: a hydraulic check valve PTO, a hydraulic active valving PTO, and a directly electrified PTO. Simple models of each PTO are constructed. Because a model for the small-scale device was initially unavailable, the PTOs are simulated on a large-scale device. The results are scaled down using Froude scaling and compared to results from directly simulating a small-scale model. Strong assumptions are made because this work is early in the design stages, and a coarse look at PTO options was desired. Specifically, the effectiveness of controls is investigated, along with the efficiency of energy conversion. However, energy capture is only part of the consideration; there are also logistic concerns to be considered when selecting a PTO. For example, components for large-scale WECs are so large and expensive that it may make sense to custom-build PTO components, but small-scale WECs would benefit from off-the-shelf availability because the cost of customization would be a significant portion of the total capital cost of deployment at a small scale. Submersible, off-the-shelf components are much easier to source for hydraulic PTOs. Because of highly effective controls, efficient energy conversion, and availability of marine-grade components, an active valving hydraulic PTO is selected for this small-scale surge-type WEC.

energy conversion/systems↗

Gamma rays and gravitational waves from inelastic Higgs portal dark matter

We explore a simple and predictive dark matter scenario involving a complex scalar field, ϕ , coupled to the Higgs portal with no additional field content. In the UV, the field possesses a global U ( 1 ) symmetry which is broken by mass terms and Higgs portal interactions. In the mass basis, the complex field splits into a pair of real scalars with a small mass splitting (in analogy to pseudo-Dirac fermions), such that the Higgs portal acquires both diagonal and off-diagonal terms with respect to these eigenstates. In the parameter space where the off-diagonal interaction predominates, this scenario is safe from direct detection constraints. Moreover, this model provides a viable explanation for the longstanding Galactic Center γ -ray excess. Additionally, this model influences the Higgs potential in a way that could facilitate a strong first-order electroweak phase transition in the early Universe, potentially leading to a stochastic gravitational wave background that could fall within the reach of upcoming space-based detectors.

Hooper, Dan [Wisconsin U., Madison] (ORCID:0009000↗

Revisiting constraints on asymmetric dark matter from collapse in white dwarf stars

The runaway collapse phase of a small dark matter cluster inside a white dwarf star encompasses a reversible stage, where heat can be transferred back and forth between nuclear and dark matter. Induced nuclear burning phases are stable and early carbon depletion undermines previous claims of type Ia supernova ignition. Instead, mini black holes are formed at the center of the star that either evaporate or accrete stellar material until a macroscopic sub-Chandrasekhar-mass black hole is formed. In the latter case, a 0.1 to 1 second lasting electromagnetic transient signal can be detected upon ejection of the white dwarf’s potential magnetic field. Binary systems that transmute to black holes and subsequently merge emit gravitational waves. Advanced LIGO should detect one such sub-Chandrasekhar binary black hole inspiral per year, while future Einstein telescope-like facilities will detect thousands per year. The effective spin parameter distribution is peaked at 0.2 and permits future studies to disentangle from primordial sub-Chandrasekhar black holes. Such signatures are compatible with current direct detection constraints, as well as with neutron star constraints in the case of bosonic dark matter, even though they remain in conflict with the fermionic case for part of the parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗

Acoustic Travel-Time Tomography for Wind Energy

New remote sensing technologies for direct observation of the atmospheric flow in the proximity of operating wind turbines are necessary for the advancement of wind turbine design, wind plant control, and the validation of high-fidelity numerical models; fundamental goals of the DOE Wind Energy Technology Office's (WETO's) A2e program. Modern remote sensing technologies for commercial and research applications rely on the backscatter of light or sound waves, which limits spatial and temporal resolution and cannot produce reliable flow estimates near solid bodies or reflective surfaces. Acoustic tomography (AT) is an innovative sensing technology based on fundamental thermodynamic properties of the atmosphere and an inverse algorithm that can reconstruct both velocity and temperature fields within a network of acoustic transducers. AT relies on the direct travel path of acoustic signals between speakers and microphones in the network and can achieve much higher resolutions than conventional backscatter technologies. For wind energy R&D applications, AT represents an extremely economical option over other remote sensing technologies, integrating mature and commonly available acoustic and signal processing hardware. This report is intended as an update on the development of the AT array on the NREL Flatirons Campus (FC) and includes a brief review of the instrumentation, the calibration of equipment, signal processing, and resulting flow field estimates. Designs for future AT systems and applications are detailed in this report, including summaries of new capabilities, engineering and theoretical challenges, and anticipated costs help to outline future directions for AT technology R\&D. The potential for application of AT to utility-scale turbines is also discussed. Investing in the development of the proposed instrument will establish AT as a transformational RS technique that represents a significant advance in measurement instrumentation, enabling new fundamental insights into the structure of the atmospheric boundary layer and supporting applications for wind energy, pollution transport, and resource characterization. AT is also an extremely economical option over other RS instruments, integrating mature and commonly available acoustic and signal processing hardware.

17 WIND ENERGY↗

Power Take-Off Design Study for a Small-Scale Oscillating Surge Wave Energy Converter for Powering the Blue Economy Applications: Preprint

The power take-off (PTO) is an integral part of wave energy conversion, and the design process is nontrivial. Better PTOs, and better processes for selecting and designing PTO architectures for various applications, would benefit devices that assist in powering the blue economy by decreasing time and money spent on PTO design and increasing the overall energy capture performance of these devices. This paper chronicles the selection process of a PTO for a small-scale surge-type wave energy converter (WEC) for the purpose of informing future PTO selection processes. Three PTO architectures are evaluated in WEC-Sim: a hydraulic check valve PTO, a hydraulic active valving PTO, and a directly electrified PTO. Simple models of each PTO are constructed. Because a model for the small-scale device was initially unavailable, the PTOs are simulated on a large-scale device. The results are scaled down using Froude scaling and compared to results from directly simulating a small-scale model. Strong assumptions are made because this work is early in the design stages, and a coarse look at PTO options was desired. Specifically, the effectiveness of controls is investigated, along with the efficiency of energy conversion. However, energy capture is only part of the consideration; there are also logistic concerns to be considered when selecting a PTO. For example, components for large-scale WECs are so large and expensive that it may make sense to custom-build PTO components, but small- scale WECs would benefit from off-the-shelf availability because the cost of customization would be a significant portion of the total capital cost of deployment at a small scale. Submersible, off-the-shelf components are much easier to source for hydraulic PTOs. Because of highly effective controls, efficient energy conversion, and availability of marine-grade components, an active valving hydraulic PTO is selected for this small-scale surge-type WEC.

electric PTO↗

Perpendicular magnetofluid theory for magnetically confined plasmas in the collisionless limit

Note that the magnetically confined plasmas in thermonuclear fusion turn to be collisionless. Some particles are trapped by the inhomogeneous equilibrium field or by the waves, while some others keep circulating. One cannot expect the charged particles to be collectively localized in space by collisions along the magnetic field lines. Consequently, the fluid parameters relating to the parallel motion, such as parallel fluid velocity, become inadequate. In this paper, the perpendic- ular magnetofluid theory for magnetically confined plasmas in the collisionless limit is developed. It revises Braginskii’s transport theory by releasing the collisional dominance assumption. In the perpendicular direction, the particle spatial localization is solo resulted from the strong magnetic field, while the finite Larmor radius (FLR) effects are taken into account in the next order. In the parallel direction, the particle mobility feature is fully retained by constructing only the per- pendicular moments. Therefore, besides the FLR effects till order O($k$ $^{2}_{⊥}$ $ρ$ $^{2}_{i}$ ), the Landau damping, the trapped particle effects, and other parallel kinetic effects are all kept in the current formalism. Here, k ⊥ is the perpendicular wave number and ρ i is the ion Larmor radius. Here, the theory can be used both for the analytical and numerical studies of magnetically confined plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗

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↗