Fact-checked by Grok 2 weeks ago

Focal mechanism

A focal mechanism describes the geometry of faulting and slip direction that generates seismic waves during an earthquake, encapsulating the orientation of the ruptured fault plane and the relative motion across it.^[1]^[2] Seismologists derive it primarily from the pattern of initial P-wave polarities—upward (compressional) or downward (dilatational)—recorded at distant stations, which delineate two nodal planes separating quadrants of opposite motion.^[3]^[4] This double-couple model is visualized in a stereographic projection called a "beach ball" diagram, where shaded regions denote compression and unshaded areas dilation, aiding inference of the earthquake's causative fault type.^[3]^[4] Principal end-member types include strike-slip (horizontal shear), normal (extensional dip-slip), and thrust or reverse (contractile dip-slip) mechanisms, each reflecting distinct tectonic stresses such as plate boundary interactions.^[4]^[5] Beyond basic fault characterization, focal mechanisms inform broader geodynamic interpretations, including regional stress orientations and seismic hazard assessment, though ambiguity between conjugate nodal planes necessitates auxiliary data like aftershock distributions for resolution.^[2]^[1]

Fundamentals

Definition and Physical Basis

A focal mechanism describes the geometric properties of the fault rupture during an earthquake, specifically the orientation of the fault plane and the direction of slip along it. It represents the seismic source as a point of equivalent forces, capturing how seismic waves radiate from the hypocenter based on the fault's strike, dip, and rake angles. This characterization arises from analysis of first-motion polarities and waveform data from seismometers, providing insight into the underlying tectonic processes.^[1]^[3] The physical basis for focal mechanisms lies in the double-couple model of earthquake faulting, which simulates the shear dislocation on a fault plane without net linear momentum change. In this model, the rupture is equivalent to two orthogonal pairs of equal and opposite forces: one pair acting to "open" the fault and the complementary pair to "close" it, reflecting the sudden release of stored elastic strain energy from tectonic loading. This formulation ensures compliance with conservation laws, as a single force couple would imply an unbalanced explosion or implosion, inconsistent with observed seismic radiation patterns. The model assumes a point source approximation valid for distances much larger than the fault dimensions, with the seismic moment tensor mathematically encoding these force distributions.^[6]^[7] For pure shear faulting, the moment tensor is traceless and symmetric, with its principal axes defining the P (compressional), T (tensional), and intermediate (N) axes that delineate the radiation pattern's quadrants. Compressional first motions occur in the P-axis quadrants due to material convergence, while dilatational motions prevail in T-axis quadrants from divergence, enabling the distinction between fault types such as strike-slip, normal, or thrust. Deviatoric components dominate in natural earthquakes, as volumetric changes (explosive or implosive) are minimal compared to shear, supported by empirical observations of nearly pure double-couple solutions in tectonic events.^[8]^[6]

Double-Couple Source Model

The double-couple (DC) source model approximates the seismic radiation from tectonic earthquakes as a point-source equivalent of shear slip along a fault plane, consisting of two perpendicular force couples that balance to produce no net force, torque, or volume change.^[6] This model arises from the displacement discontinuity across the fault, where the far-field elastic waves mimic the pattern generated by forces of magnitude F separated by distance d, with seismic moment M_0 = Fd.^[9] Unlike single-force or isotropic models, the DC configuration aligns with the absence of overall mass displacement or dilatation in shear-dominated events, as confirmed by matching observed first-motion polarities and waveform amplitudes in global seismicity catalogs.^[8] Mathematically, the DC source is embodied in the deviatoric component of the seismic moment tensor \mathbf{M}, a symmetric 3x3 matrix with trace zero (M_{xx} + M_{yy} + M_{zz} = 0), where the principal axes define the fault-normal compression (P-axis), fault-parallel extension (T-axis), and intermediate null (B-axis) directions.^[7] The radiation pattern for P-waves features four quadrants of alternating compressional and dilatational first motions, while S-waves exhibit a four-lobed pattern, both invariant under 180° rotation about the source axes due to the model's symmetry.^[10] This traceless tensor excludes isotropic (explosive/implosive) contributions, which would imply net volume change and produce uniform P-wave expansion or contraction; empirical inversions of tectonic events routinely yield near-zero isotropic components (typically <5% of total moment), supporting DC dominance over alternatives like compensated linear vector dipoles (CLVD) in standard faulting.^[11]^[12] The model's validity stems from its equivalence to the point-source limit of finite-fault rupture under elastic rebound theory, where prestress accumulation and sudden release on a planar discontinuity generate pure shear without requiring tensile cracking or fluid involvement, as evidenced by dynamic simulations matching broadband seismograms from events like the 1994 Northridge earthquake (moment tensor deviatoric purity >95%).^[6] Deviations to non-DC components, observed in <10% of cataloged earthquakes (e.g., volcanic or induced seismicity), often reflect modeling artifacts, near-source heterogeneity, or atypical processes like tensile failure, but do not invalidate the DC as the baseline for crustal shear sources.^[13] Quantitative fits using DC-constrained inversions, such as those minimizing waveform misfits across global networks, achieve variance reductions of 80-90% for magnitudes M_w > 5, underscoring its causal fidelity to fault mechanics over isotropic or volumetric alternatives.^[8]^[1]

Historical Development

Early Determination Techniques

The earliest techniques for determining focal mechanisms relied on analyzing the polarities of first-arriving P-waves recorded at distant seismograph stations. These methods, pioneered in the 1920s, interpreted the initial ground motion—either compressional (upward deflection) or dilatational (downward deflection)—to infer the orientation of the earthquake source. Percy Byerly introduced a foundational approach in 1928, modeling the source as a single force couple to explain the observed "fling" motion across a fault, though subsequent refinements recognized the double-couple nature of shear faulting for consistency with moment conservation.^[14]^[15] Byerly's method involved projecting station azimuths and takeoff angles onto a focal sphere centered at the earthquake hypocenter, then plotting compressional (+) and dilatational (-) polarities to identify boundaries separating regions of opposite motion. These boundaries formed two orthogonal nodal planes—one interpreted as the fault plane and the other as an auxiliary plane—defining four alternating quadrants of compression and dilatation consistent with a double-couple radiation pattern. Manual plotting on stereographic nets, such as the Wulff or equal-area projection, was used to fit the data, requiring at least 10-20 reliable polarities with good azimuthal coverage for resolution. Early applications, like Byerly's 1938 solution for the July 6, 1934, northern California offshore earthquake, demonstrated strike-slip faulting, validating the technique against geological evidence.^[15]^[16] Challenges in these pre-digital era methods included uncertainties in epicentral location (typically fixed a priori), ray path distortions from heterogeneous Earth structure, and ambiguous polarities due to noise or near-station effects, often leading to non-unique solutions where either nodal plane could be the actual fault. Independent efforts in Japan and Europe during the 1930s paralleled Byerly's work, using similar first-motion data but limited by sparse global station networks, primarily yielding solutions for large-magnitude events (M > 6) detectable teleseismically. Refinements, such as tables extending Byerly's techniques to deep-focus earthquakes by Hodgson and Storey in 1953, improved applicability but retained the labor-intensive graphical fitting process.^[17]^[16]

Transition to Quantitative Solutions

The initial approaches to focal mechanism determination, pioneered in the early 20th century, relied on qualitative graphical techniques that plotted P-wave first-motion polarities (compressions and dilatations) on equal-area stereographic projections, such as Wulff or Schmidt nets, to delineate nodal planes separating the four quadrants of predicted motions.^[18] These manual methods, exemplified by Percy Willis's work in the 1920s and Percy Byerly's systematic application in the 1930s, involved trial-and-error fitting to identify double-couple orientations consistent with observations, but they were inherently subjective, limited by sparse data, and prone to ambiguity due to the trade-off between fault and auxiliary planes.^[19] The transition to quantitative solutions accelerated in the late 1960s with the emergence of digital computers, enabling numerical optimization algorithms to objectively minimize misfits between observed and synthetic polarity patterns. Early programs, such as that developed by Udías and Baumann in 1968–1969, combined P- and S-wave data in grid-search or least-squares inversions to compute best-fit strike, dip, rake, and scalar moment, providing uncertainty estimates and resolving ambiguities more reliably than graphical trial-and-error.^[20] This shift incorporated additional constraints like amplitude ratios (e.g., P/S or SH/P), reducing reliance on polarities alone and allowing systematic processing of larger station datasets.^[21] By the 1970s, these computational methods facilitated the compilation of global focal mechanism catalogs, marking a departure from ad hoc solutions toward reproducible, data-driven analyses that integrated theoretical double-couple models with empirical observations, though challenges like velocity model uncertainties persisted.^[18]

Methods of Determination

Polarity-Based Approaches

Polarity-based approaches to focal mechanism determination utilize the observed first-motion polarities of P-waves recorded at seismic stations to constrain the orientation of the fault and auxiliary planes. In the double-couple model, the P-wave radiation pattern partitions the focal sphere into alternating compressional (+) and dilatational (-) quadrants, bounded by nodal planes corresponding to the fault plane and its orthogonal auxiliary plane.^[22] The polarity at a station depends on its location relative to these quadrants, determined by the azimuthal takeoff angle and the station's epicentral distance, which influences the ray's emergence angle.^[23] The classical method, developed by Perry Byerly in the 1930s, involves graphically plotting station positions on a stereographic projection of the focal sphere and iteratively adjusting nodal plane orientations to achieve consistency between observed and predicted polarities, allowing for a small percentage of inconsistent readings due to noise or scattering.^[15] Byerly's approach required manual construction of great circles separating clusters of + and - polarities, often using equal-area or Wulff nets to visualize the distribution.^[24] This graphical technique was limited to well-recorded events with azimuthal coverage exceeding 180 degrees to resolve ambiguity between fault and auxiliary planes.^[25] Quantitative polarity-based methods emerged in the mid-20th century, employing grid-search algorithms to systematically test possible double-couple orientations and minimize the number of polarity misfits, typically scored via the percentage of consistent observations or a weighted least-squares metric accounting for takeoff angle uncertainties.^[26] These approaches often incorporate error ellipsoids for station locations and assume a radially symmetric Earth model for ray tracing, with solutions considered reliable if over 80% of polarities match the predicted pattern.^[27] Limitations include sensitivity to near-nodal recordings, where weak signals lead to ambiguous polarities, and assumptions of isotropic media, which can introduce biases in heterogeneous crust.^[23] Contemporary implementations integrate polarity data with auxiliary constraints like S/P amplitude ratios to enhance resolution, particularly for shallow events or sparse networks.^[22] Automated polarity picking using convolutional neural networks has improved efficiency and accuracy for microseismicity, enabling rapid focal mechanism solutions by preprocessing waveforms to classify initial motions as up, down, or unclear with probabilities exceeding 90% for high signal-to-noise ratios.^[28] Such deep learning tools, trained on labeled datasets from regional arrays, reduce human bias and handle impulsive onsets better than traditional STA/LTA detectors, though they require validation against manual reads for novel tectonic regimes.^[29]

Amplitude and Waveform Methods

Amplitude methods for determining focal mechanisms utilize ratios of seismic wave amplitudes, such as S/P or SH/P, derived from observed seismograms at recording stations. These ratios capture the directional dependence of the seismic radiation pattern from a double-couple source, providing constraints that complement first-motion polarities by reducing ambiguities in fault orientation.^[26] For small earthquakes, S/P amplitude ratios are measured directly from waveform envelopes, then inverted alongside polarity data to estimate strike, dip, and rake parameters, with applications demonstrated in regions like southern California where they refine solutions for magnitudes as low as M 2.0.^[30] SH/P ratios, in particular, prove effective in regional networks for selecting the preferred mechanism among multiple polarity-consistent candidates, as their relative insensitivity to isotropic path effects enhances resolution of nodal plane trade-offs.^[31] Despite their utility, amplitude ratio methods are susceptible to biases from near-source effects like directivity, site-specific amplification, and heterogeneous attenuation, which can systematically skew solutions toward thrust or strike-slip orientations if uncorrected.^[32] Studies of induced seismicity highlight that while P- and S-wave amplitudes can be inverted jointly, preprocessing for instrument response and windowing around peak arrivals is essential to avoid contamination from scattered phases.^[26] Empirical corrections, such as those based on 3D velocity models, mitigate these issues but require dense station coverage for validation.^[33] Waveform methods advance beyond ratios by fitting the full time series of seismic traces, incorporating both amplitude and phase information across P, S, and surface waves to invert for the moment tensor. High-frequency full waveform analysis suits small local events (M < 3) in sparse arrays, where grid-search or least-squares optimization matches synthetic seismograms—computed via modal summation or finite-difference modeling—to data, resolving mechanisms unattainable with polarities alone due to cycle-skipping avoidance in broadband fitting.^[34] This approach exploits the sensitivity of waveform shape to source depth and orientation, as validated in oil reservoir monitoring where it yields stable double-couple solutions with rake uncertainties under 20 degrees.^[35] Advanced implementations employ full-waveform inversion (FWI) frameworks, minimizing L2 misfit between observed and simulated waveforms through Bayesian inference or gradient-based optimization, often retrieving isotropic and compensated linear vector dipole components alongside the double-couple.^[36] For moderate events like the 2018 M7.1 Anchorage earthquake, joint inversions of broadband seismograms and GPS data via cut-and-paste techniques confirm mechanisms with scalar moments accurate to within 10%, though computational demands limit real-time use without simplified Green's functions.^[37] Waveform methods outperform polarity-based ones in constraining non-standard sources, such as in volcanic or induced seismicity, but demand accurate velocity models to prevent trade-offs between structure and source parameters.^[38]

Advanced Computational Techniques

Full-waveform inversion (FWI) represents a sophisticated computational approach for determining focal mechanisms by minimizing the misfit between observed seismograms and synthetic waveforms generated from parameterized source models, enabling the retrieval of moment tensors for events as small as magnitude 0.8.^[39] This technique integrates broadband data across frequency bands, incorporating path effects and site responses, which enhances accuracy over polarity-only methods, particularly in complex velocity structures.^[34] For instance, multiparameter FWI simultaneously inverts for event location, focal mechanism, and moment magnitude, as demonstrated in applications to microseismic datasets where traditional methods fail due to low signal-to-noise ratios.^[40] Deep learning algorithms, such as the Focal Mechanism Network (FMNet), leverage convolutional neural networks trained on three-component waveform data to estimate focal mechanisms in real-time, achieving errors below 20 degrees for strike, dip, and rake angles in regional networks.^[41] These models process full waveforms directly, bypassing explicit phase picking, and outperform P-wave first-motion approaches by incorporating amplitude and phase information, with validation on datasets from the 2019 Ridgecrest sequence showing robust performance for magnitudes above 3.0.^[41] Extensions include multitask networks that jointly predict phase arrivals, polarities, and mechanisms, facilitating end-to-end earthquake monitoring pipelines.^[42] Probabilistic and grid-search algorithms, like the REFOC method, combine P- and S-wave polarities with inter-event amplitude ratios to resolve non-double-couple components and improve stability in sparse networks, as tested on California seismicity data yielding variance reductions of up to 30% compared to legacy inversions.^[22] Open-source tools such as SKHASH implement hash-based grid searches over fault parameters, incorporating S/P amplitude ratios and machine learning constraints for rapid, automated solutions suitable for real-time catalogs.^[43] These techniques collectively advance focal mechanism resolution by exploiting computational power for waveform fitting and statistical inference, though they require high-quality velocity models and dense arrays to mitigate trade-offs between source depth and mechanism orientation.^[44]

Graphical and Mathematical Representations

Beachball Diagrams

Beachball diagrams provide a visual representation of an earthquake's focal mechanism by projecting the lower-hemisphere stereographic view of the P-wave radiation pattern onto a two-dimensional circle.^[45] The diagram divides the focal sphere into four quadrants separated by two orthogonal nodal planes, corresponding to the fault plane and its auxiliary plane in the double-couple model.^[4] Shaded (typically black) regions indicate areas of compressional first-motion polarity for P-waves (upward particle motion), while unshaded (white) regions denote dilatational polarity (downward motion), based on empirical observations from seismograms.^[1] This pattern arises from the constructive and destructive interference of seismic waves emitted from the slipping fault, reflecting the geometry of slip relative to the observer's takeoff angle and azimuth.^[46] The construction begins with polarity data from multiple seismic stations, where the takeoff angle (from source to receiver) determines whether a station records a compressional or dilatational onset. These polarities are plotted on the focal sphere, and the great circles of zero amplitude (nodal planes) are fitted to separate positive and negative quadrants. The resulting diagram compresses the lower hemisphere into a circle, with the center representing the downward vertical and the periphery the horizontal plane; lines curving across the diagram trace the projections of the nodal planes.^[3] Standard conventions ensure consistency: for instance, in strike-slip events, alternating black and white quadrants form a characteristic "X" or crossed pattern, while normal faulting shows larger white lower quadrants and thrust faulting the opposite.^[47] Oblique mechanisms produce irregular quadrant sizes, complicating plane identification without auxiliary data like aftershock distributions.^[46] Interpreting beachballs requires distinguishing the actual fault plane, often aided by the relative motion vectors implied by the double-couple: compression along the P-axis (where all quadrants converge to black) and dilation along the T-axis. The auxiliary plane bisects the acute angle between P and T axes, while the fault plane aligns with the slip direction. Uncertainties in polarity readings or non-double-couple components can distort the pattern, necessitating full moment tensor inversions for precision. These diagrams, despite their simplicity, effectively summarize fault orientation and type, aiding tectonic reconstructions when compiled from regional events.^[1]^[45]

Moment Tensor Formalism

The seismic moment tensor provides a mathematical description of an earthquake source as a point-source approximation, represented by a symmetric 3×3 tensor \mathbf{M} with components M_{ij} (where i,j = 1,2,3 corresponding to orthogonal axes, often east-north-up).^[48]^[49] This tensor quantifies the equivalent body forces at the hypocenter that generate seismic waves, with the far-field displacement u_k proportional to the second time derivative \ddot{M}_{ij} n_i n_j / r, where \mathbf{n} is the unit vector from source to receiver, r is distance, and the proportionality involves medium density and wave speed. The symmetry (M_{ij} = M_{ji}) arises from angular momentum conservation, yielding six independent components: three diagonal (M_{xx}, M_{yy}, M_{zz}) representing volumetric changes or extensions, and three off-diagonal (M_{xy}, M_{xz}, M_{yz}) capturing shear couples.^[50] In the general case, the tensor decomposes into an isotropic component M_{iso} = (M_{xx} + M_{yy} + M_{zz})/3, indicating net volume change (positive for expansion, negative for implosion), and a deviatoric component with zero trace, which includes shear faulting.^[49]^[50] The deviatoric part further splits into a double-couple (DC) mechanism, modeling slip on a planar fault without volume change, and a compensated linear vector dipole (CLVD) for tensile or compressive cracking orthogonal to the DC.^[5] For pure double-couple sources, typical of tectonic earthquakes, the tensor takes the form M_{ij} = M_0 (\nu_i \delta_j + \nu_j \delta_i), where M_0 is the scalar seismic moment, \boldsymbol{\nu} is the fault-normal vector, and \boldsymbol{\delta} is the slip vector in the auxiliary coordinate system.^[48] The scalar moment M_0 measures source size as M_0 = \mu A D, with \mu shear modulus, A fault area, and D average slip, linking the tensor norm to physical rupture parameters via M_0 = \sqrt{\frac{1}{2} \sum M_{ij}^2} for deviatoric tensors.^[51] This formalism extends beyond ideal shear dislocations to non-double-couple sources, such as explosions (dominant isotropic) or volcanic events, by allowing full tensor inversion from waveform data, though resolution depends on network geometry and signal-to-noise ratio.^[5]^[52] Eigenvalue decomposition of \mathbf{M} yields principal axes P (compression), N (null), and T (tension), facilitating interpretation of source orientation independent of geographic coordinates.^[50] Uncertainties in components arise from trade-offs, particularly between isotropic and CLVD for near-vertical faults, requiring constraints like zero-trace assumption for tectonic events.^[52]

Interpretations and Classifications

Standard Fault Types

Standard fault types in focal mechanisms represent end-member double-couple models of shear slip on planar faults, assuming no isotropic (volumetric) component and classifying motion by the rake angle—the direction of slip relative to fault strike. Pure strike-slip occurs at rake angles near 0° or 180°, normal dip-slip near -90°, and reverse dip-slip near +90°; these correspond to distinct patterns in lower-hemisphere projections (beachballs), where shaded quadrants denote compressional first-motion polarities and unshaded denote dilatational.^[45]^[1] Strike-slip faults exhibit horizontal relative motion parallel to the fault strike on near-vertical planes (dip ≈90°), driven by lateral shear without significant vertical offset. Right-lateral (dextral) slip displaces the opposite block to the right when facing the fault, as on California's San Andreas Fault; left-lateral (sinistral) does the opposite. Beachball diagrams show a cross-pattern of two adjacent shaded and unshaded quadrants, with vertical nodal planes—the actual fault and orthogonal auxiliary—producing symmetric compression along the maximum compressive stress (σ1) direction perpendicular to slip.^[4]^[1]^[45] Normal (extensional) faults involve dip-slip motion where the hanging wall slides downward relative to the footwall along inclined planes (typical dip 45°–70°), accommodating crustal extension in rift zones. The tension axis (T) aligns subhorizontally perpendicular to the fault strike, and the pressure axis (P) subvertically; beachballs display an asymmetrical pattern with opposing shaded quadrants concentrated toward the steeper-dipping nodal plane, reflecting P-wave compressions in deeper lobes.^[4]^[1]^[45] Reverse (compressional) faults, including thrusts, feature upward hanging-wall motion relative to the footwall, with variable dips but often shallow (<30°) for thrusts under horizontal compression. The P-axis aligns subhorizontally parallel to σ1, and T subvertically; beachballs mirror normal faults but inverted, with shaded quadrants emphasizing shallow compressional lobes and the shallower nodal plane as the likely fault.^[4]^[1]^[45] These types assume pure end-members; real mechanisms often include oblique components, and fault-plane ambiguity is resolved via auxiliary data like aftershock lineations aligning with the true rupture plane.^[1]^[45]

Deviations from Double-Couple Assumptions

The double-couple model underlying standard focal mechanisms assumes pure shear deformation on a fault plane, implying no net volume change and equivalence to orthogonal quadrantal force couples. Deviations arise when the seismic source includes non-double-couple (non-DC) components, detectable through full moment tensor inversions that decompose the tensor into double-couple (DC), isotropic (ISO), and compensated linear vector dipole (CLVD) parts, where the DC fraction is quantified as the percentage of the tensor's deviatoric energy attributed to shear slip. Non-DC components, comprising ISO (spherical expansion or contraction) and CLVD (elongated tension or compression without net volume change), typically constitute less than 10% of the moment tensor for most tectonic earthquakes but can exceed 20-50% in induced, mining, or volcanic events.^[53]^[54] Isotropic components represent volumetric changes, with positive ISO indicating explosive sources like nuclear tests or cavity collapses, and negative ISO suggesting implosions such as reservoir-induced seismicity or magma chamber deflation; for instance, analysis of the 2019 Ridgecrest sequence revealed ISO components up to 15% in some foreshocks, potentially linked to pore pressure variations or fluid interactions rather than pure faulting. CLVD components, often vertical in orientation, model tensile crack opening or shear-enhanced dilation, prevalent in volcanic settings where 101 vertical-CLVD events near active volcanoes were identified with CLVD percentages reaching 40%, attributed to dyke propagation or ring-fault mechanics. In mining-induced seismicity, non-DC fractions as high as 30-60% have been documented, reflecting combined shear, tensile, and implosive processes in heterogeneous rock masses.^[55]^[56]^[57] Apparent non-DC components in tectonic catalogs, with mean deviations from pure DC around 5-10%, may partly stem from inversion artifacts such as unmodeled path effects, noise, or suboptimal station coverage rather than true source physics, as synthetic tests show that full moment tensor solutions can amplify minor asymmetries into spurious ISO/CLVD signals unless constrained by double-couple assumptions. Consistency checks across catalogs using varied inversion methods reveal that while some non-DC signals persist (e.g., in swarm seismicity), many align with DC-only solutions upon reanalysis, underscoring the need for multiple datasets to distinguish genuine deviations—such as those in deep subduction zones from phase changes—from methodological biases. Verification through waveform polarity constraints or geodetic data integration enhances reliability, as unconstrained inversions risk overestimating non-DC by factors of 2-3.^[54]^[58]^[59]

Applications in Seismology

Tectonic and Stress Analysis

Focal mechanisms provide direct evidence of the orientation of fault planes and slip vectors during earthquakes, enabling the classification of tectonic regimes at plate boundaries and intraplate settings. In extensional environments, such as mid-ocean ridges or continental rifts, normal faulting dominates, with the principal stress axis σ3 (least compressive) oriented subhorizontally perpendicular to the rift axis, consistent with gravitational unloading driving divergence.^[60] Compressional regimes at subduction zones feature thrust or reverse faults, where σ1 (maximum compressive stress) aligns subhorizontally and converges toward the trench, as exemplified by the 2004 Sumatra-Andaman earthquake's megathrust mechanism with a near-unity rake on a low-angle plane.^[3] Strike-slip faulting prevails along transform boundaries like the San Andreas Fault, where σ1 and σ3 are both subhorizontal and at approximately 45° to the fault strike, reflecting shear-dominated plate motion.^[61] Stress analysis employs inversion of populations of focal mechanisms to resolve regional principal stress orientations and shape ratios, assuming that slip occurs on favorably oriented planes under a uniform stress field as per Anderson's faulting theory. Techniques like the Michael (1984) method minimize misfit between observed slip vectors and those predicted by a stress tensor, handling nodal plane ambiguity by testing both possibilities and yielding σ1, σ2, and σ3 directions with associated uncertainties.^[61] Advanced Bayesian or bootstrapped approaches further quantify errors from heterogeneous stress or non-double-couple components, as applied in southern California where inversions reveal σ1 oriented northwest-southeast, aligning with Pacific-North American plate convergence. These inversions have mapped stress provinces globally, identifying rotations or heterogeneities, such as in the Himalayan arc where local variations deviate from India-Eurasia collision predictions.^[62] In intraplate settings, focal mechanisms reveal inherited or reactivated structures under far-field plate stresses; for instance, in stable cratons, compressional mechanisms indicate distant subduction influences, with inversion results showing σ1 aligned with absolute plate motion.^[63] Such analyses constrain geodynamic models by linking seismicity patterns to lithospheric strength and viscosity, though assumptions of homogeneity limit resolution in complex zones.^[64] Empirical validations against borehole breakouts or GPS data confirm inversion reliability where data density suffices, typically requiring dozens of well-constrained mechanisms for robust axes orientations within 10-20° uncertainty.^[65]

Seismic Hazard and Risk Assessment

Focal mechanisms play a critical role in seismic hazard assessment by providing constraints on fault orientation, slip type, and rupture characteristics, which are essential for modeling earthquake sources in probabilistic seismic hazard analysis (PSHA). In PSHA, these mechanisms inform the characterization of seismogenic sources, including fault geometries and maximum magnitudes, enabling more accurate estimation of ground motion exceedance probabilities. For instance, strike-slip mechanisms, common along transform faults, differ in radiation patterns and directivity effects compared to thrust mechanisms in subduction zones, influencing the selection of ground motion prediction equations (GMPEs) that account for style-of-faulting variability.^[66]^[67] Incorporating focal mechanism data into PSHA deaggregation allows identification of controlling earthquakes, where the resulting design spectrum is tied to specific source parameters like rake angle and fault strike, potentially altering hazard levels by up to 20-30% in regions with heterogeneous stress fields. This approach extends traditional PSHA by introducing a priori focal mechanism distributions, reducing epistemic uncertainty in source models, particularly in areas with limited historical data. Real-time focal mechanism solutions further enhance operational hazard monitoring, as rapid assessment of slip type aids in forecasting aftershock sequences and evaluating tsunami potential for shallow thrust events.^[66]^[38] In seismic risk assessment, focal mechanisms contribute to scenario-based simulations that integrate hazard with exposure and vulnerability, such as determining amplified shaking on hanging walls during reverse faulting or forward directivity in unilateral ruptures along strike-slip faults. Regional studies, like those in subduction-prone areas, use cataloged mechanisms to define zonation parameters, including predominant fault types and depths, which refine loss estimates for infrastructure. For example, in Jamaica's PSHA, zones are characterized by dominant normal or strike-slip mechanisms to parameterize recurrence and attenuation, directly impacting building code provisions and insurance modeling. Limitations arise from non-double-couple components in complex sources, which may bias risk if not resolved, underscoring the need for high-quality, diverse datasets from global catalogs.^[68]^[38]

Monitoring Induced and Volcanic Seismicity

Focal mechanisms are essential for characterizing induced seismicity from activities like hydraulic fracturing, wastewater injection, and geothermal stimulation, as they delineate fault plane orientations, slip directions, and stress perturbations that inform operational risks. In the Permian Basin of South Texas, an enhanced catalog of over 1,000 focal mechanisms for induced earthquakes between 2017 and 2022 revealed dominant strike-slip and reverse faulting aligned with regional extensional stresses, facilitating improved mapping of reactivated faults and rupture kinematics.^[69] Advanced methods integrating borehole data and in situ stress estimates have constrained first-motion polarities for small-magnitude events (M < 2), showing shear-dominated double-couple solutions consistent with poroelastic responses to injection pressures rather than tensile opening.^[70] These analyses enable real-time traffic-light systems for mitigating seismicity, as demonstrated in monitoring programs where non-double-couple components indicate volumetric changes from fluid pathways.^[26] In volcanic environments, focal mechanisms of volcano-tectonic (VT) earthquakes—typically double-couple types reflecting brittle shear failure—serve as proxies for subsurface stress regimes and magma-induced faulting, aiding eruption forecasting. At Campi Flegrei caldera, Italy, focal mechanisms from 2012–2019 swarms predominantly featured normal faulting on high-angle planes (strike ~300°–330°, dip ~60°–80°), correlating with extensional tectonics and unrest phases.^[71] Temporal variations in these mechanisms, such as shifts from strike-slip to normal during precursory activity, have been linked to stress perturbations from magma ascent, as evidenced at Mount Etna where 2001 eruption sequences showed evolving P- and T-axes indicative of localized extension.^[72] Clustering of moment tensors further classifies event families, distinguishing VT shear events from hybrid or long-period signals with isotropic components that suggest fluid excitation or crack resonance, enhancing multi-parametric monitoring networks.^[73] Such determinations, often derived from dense seismic arrays, reveal volcanic stress fields influenced by edifice loading and plate tectonics, with reliability improved by waveform fitting over polarity alone.^[74]

Limitations and Uncertainties

Sources of Error in Solutions

Focal mechanism solutions are inherently non-unique due to the symmetry of the double-couple model, which produces two orthogonal nodal planes—one representing the actual fault and the other an auxiliary plane—without additional geologic or geodetic constraints to distinguish them.^[75] This ambiguity persists in standard inversions unless resolved by relative relocation techniques or waveform fitting that favor one plane based on rupture directivity.^[76] Data quality introduces significant errors, particularly low signal-to-noise ratios from sources like oceanic microseisms, wind, or cultural noise, which degrade polarity readings and waveform fits.^[77] For small-magnitude events (M < 3), limited high-quality recordings exacerbate this, often resulting in poorly constrained solutions due to insufficient amplitude or phase information.^[33] Azimuthal gaps in station coverage further amplify uncertainties by poorly sampling the radiation pattern, leading to biased parameter estimates such as strike or rake angles.^[78] Modeling assumptions contribute to discrepancies, including inaccuracies in crustal velocity models that distort Green's functions and introduce biases in moment tensor components; for instance, a 5% perturbation in P- or S-wave velocities can shift solutions by tens of degrees. Source location errors, especially in hypocentral depth, compound these issues by misaligning observed and synthetic seismograms, with joint inversions sometimes mitigating but not eliminating the effect.^[26] Amplitude ratio methods, such as S/P ratios, can bias inversions due to unmodeled distance decay (e.g., ratios dropping from ~0.55 to ~0.15 over 50–120 km for M1 events) and magnitude-dependent attenuation, worsening ~20% of solutions despite average improvements.^[79]

Ongoing Debates on Source Complexity

One persistent debate in focal mechanism analysis concerns the prevalence and physical interpretation of non-double-couple (non-DC) components, which deviate from the standard shear-faulting model assumed in double-couple (DC) representations. These components, including isotropic (indicating volume changes) and compensated linear vector dipole (CLVD) elements, have been reported in up to 10-20% of cataloged events, particularly in induced seismicity, volcanic settings, and mining tremors, but their detection is hindered by resolution limits in moment tensor inversions. Critics argue that apparent non-DC signatures often arise as artifacts from velocity model errors, inadequate station coverage, or unmodeled anisotropy, rather than genuine tensile cracking or fluid-driven processes.^[80]^[81]^[82] Proponents of physical non-DC mechanisms cite cases like induced earthquakes in fluid-injection zones, where non-DC fractions correlate with pore-pressure changes and hybrid faulting, suggesting mechanisms beyond pure shear slip. For instance, analysis of events in the Permian Basin revealed non-DC components linked to fault reactivation influenced by hydraulic fracturing, challenging the DC-only paradigm for hazard modeling. However, comparative studies across global catalogs show high variability in non-DC estimates—ranging from negligible to significant—due to differences in inversion methods, waveform frequency bands, and prior constraints on tensor decomposition, underscoring unresolved trade-offs between data fit and model parsimony.^[83]^[84]^[13] A related contention involves source time function complexity, where multi-sub-event ruptures or directivity effects mimic non-DC signals in low-frequency inversions, complicating distinctions between true complexity and observational biases. Recent work on the 2024 Noto Peninsula earthquake highlighted discrepancies between first-motion polarities (favoring DC strike-slip) and full-waveform inversions (yielding non-DC elements), attributed potentially to unmodeled heterogeneity or actual implosive components from fluid migration. This fuels arguments for waveform-polarity hybrid methods to better isolate non-DC reliability, yet empirical tests indicate that even high-quality datasets yield unstable decompositions without broad azimuthal coverage.^[85]^[86] In tectonic earthquakes, the rarity of robust non-DC detections—fewer than 5% in well-constrained catalogs—supports the dominance of DC models for large events, but microseismicity studies reveal frequent CLVD-like patterns potentially tied to fault zone anisotropy or afterslip. Debates persist on whether advancing full moment tensor (including single-force terms) analyses will validate non-DC as routine indicators of rupture physics or expose them as predominantly methodological artifacts, with implications for refining seismic hazard maps that assume simple faulting. Ongoing efforts emphasize synthetic benchmarks and machine-learning aids for decomposition stability, though consensus remains elusive amid trade-offs in computational feasibility and physical interpretability.^[87]^[84]

References

[1]
Focal Mechanisms - Pacific Northwest Seismic Network
Seismologists refer to the direction of slip in an earthquake and the orientation of the fault on which it occurs as the focal mechanism.
[2]
Focal Mechanism - an overview | ScienceDirect Topics
Focal mechanisms refer to the methods used to analyze seismic waves radiated from an earthquake to determine the fault's geometry and motion during rupture.
[3]
Earthquake focal mechanisms – life's a beach(ball) - USGS.gov
Jun 24, 2024 · The fault orientation and the slip direction of an earthquake are described by something called a focal mechanism. And focal mechanisms are usually visualized ...
[4]
Focal Mechanisms Explained - IRIS
When an earthquake occurs, seismologists create graphics of focal mechanisms, informally referred to as beach balls,to show the faulting motions that ...
[5]
Berkeley Moment Tensor Project
Types of Faulting These diagrams are called focal mechanism diagrams, or sometimes “beach-ball” diagrams. These circular diagrams shown in the following figure ...<|separator|>
[6]
[PDF] GS 388 Handout: Earthquake Focal Mechanisms 1 DOUBLE COUPLE
The elastic wave radiation from an earthquake can be modeled in two equivalent ways. The first uses the concept of of a point force, i.e. a force applied ...
[7]
[PDF] 12.510 Introduction to Seismology - MIT OpenCourseWare
A quick review of the Moment Tensor and focal mechanisms. For a point ... momentum conservation, a complementary couple balances the double couple forces.<|separator|>
[8]
[PDF] EARTHQUAKE, FOCAL MECHANISM - Northwestern University
University, Evanston, IL, USA. Definition. Focal mechanisms are geometrical or mathematical repre- sentations of the faulting during an earthquake. In very.
[9]
[PDF] 4.02 Seismic Source Theory - Elsevier
The next level of complexity was to introduce a double couple source, a source without resultant force or moment. The physical origin of the double couple model ...
[10]
Double-couple earthquake source: symmetry and rotation
The face/side (up/down or foot/hangingwall) of a focal plane generally is unknown. In such a case, the DC focal mechanism has a C2 symmetry; we call it DC2.
[11]
Earthquake Source Complexity Controls the Frequency ...
Aug 19, 2021 · The spatial patterns of earthquake ground motion amplitudes are commonly represented using a double-couple model that corresponds to shear ...
[12]
Non-double-couple earthquakes 2. Observations
Most studies assume that earthquakes have double-couple (DC) source ... isotropic components, indicating volume changes of either explosive or implosive polarity.
[13]
Non‐Double‐Couple Components of Seismic Source: Method and ...
Apr 11, 2024 · Those unique source components are usually represented by two non-double-couple (non-DC) components in a seismic moment decomposition: an ...
[14]
[PDF] Commonly used earthquake source models - Geologos
Byerly (1928) was the first to solve this problem by presenting a single force couple model to identify the fling motion of the two sides of a fault. An im- ...Missing: Percy | Show results with:Percy
[15]
Fifteen— Development of Fault-Plane Studies for the Mechanism of ...
The first methods to study the source mechanism of earthquakes were based on observations of the compressional or dilatational nature of the first impulse of ...
[16]
Tables extending byerly's fault-plane techniques to earthquakes of ...
Mar 3, 2017 · Byerly has developed a method by which the direction of faulting in a normal-focus earthquake may be determined from the distribution of ...Missing: Percy | Show results with:Percy
[17]
Observatory Seismology "d0e9506"
The Focal Sphere. Independent of Byerly's work, studies of the mechanisms of earthquakes using first-motion data were pursued in Japan and Europe.Missing: Percy | Show results with:Percy
[18]
FMC—Earthquake focal mechanisms data management, cluster and ...
Since the 70s, earthquake focal mechanisms started to be computed in a systematic way and global catalogues of focal mechanisms. Because of the continuous ...Missing: transition | Show results with:transition
[19]
The S-wave project for focal mechanism studies, earthquakes of 1963
Mar 3, 2017 · The solution of a mid-Atlantic earthquake of November 17, 1963 is presented as an example of a focus which clearly conforms to the double couple ...Missing: assisted | Show results with:assisted
[20]
A computer program for focal mechanism determination combining ...
Mar 3, 2017 · A computer program for focal mechanism determination combining P and S wave data Available ... Article history. Received: 10 Jun 1968. First ...
[21]
[PDF] Earthquake Focal Mechanisms from Relative Amplitudes of P, pP ...
Earthquake Focal Mechanism Solutions 1922-1962". Pubi Dom Obs,. Ottawa, 33, No. 1 (1967). 15. A Udias and D Baumann: "A Computer Program for Focal Mechanism.
[22]
A New Focal Mechanism Calculation Algorithm (REFOC) Using Inter ...
Feb 13, 2023 · We present a new focal mechanism calculation method using polarities, S/P amplitude ratios and inter-event P/P and S/S amplitude ratios The ...
[23]
P-wave first-motion polarity determination of waveform data in ...
Nov 29, 2019 · P-wave first-motion polarity is the most useful information in determining the focal mechanisms of earthquakes, particularly for smaller ...
[24]
Remark on Byerly's fault-plane method - GeoScienceWorld
Mar 3, 2017 · In Byerly's fault-plane method, if one circle of radius R 2 is determined, the locus of the centers of the other is a straight line perpendicular to the line ...Missing: Percy | Show results with:Percy
[25]
The Focal Mechanism of Earthquakes - ScienceDirect.com
This chapter discusses the development of Byerly's method of fault-plane determination. The essential features of the methods of focal mechanism study are ...
[26]
Focal mechanism determination for induced seismicity using the ...
In this study, we have developed a new method to determine the focal mechanism for induced seismicity. Three misfit functions are used in our method to measure ...Missing: assisted | Show results with:assisted
[27]
P Wave Arrival Picking and First‐Motion Polarity Determination With ...
Jun 4, 2018 · The first-motion polarities are used to determine earthquake focal mechanisms. We use all of the predicted first-motions from the validation ...Abstract · Introduction · Methods and Results · Discussion
[28]
Simultaneous Seismic Phase Picking and Polarity Determination ...
Jan 19, 2023 · We employ a state‐of‐the‐art neural network infused with an attention mechanism to simultaneously pick arrivals and determine the first‐motion polarity.
[29]
PoViT-UQ: P-wave polarity and arrival time determination using ...
When focal mechanisms are estimated using P-wave first-motion polarities, accurate polarity determination is critical. In recent years, deep-learning-based ...
[30]
Estimating Focal Mechanism of Small Earthquakes Using S/P ...
May 3, 2025 · We propose a method for estimating focal mechanisms using data from a novel fiber-optic sensing technique called distributed acoustic sensing ( ...
[31]
Regional Network focal mechanism determination using 3D velocity ...
Mar 3, 2017 · The SH/P amplitude ratios are found to be very effective as constraints for choosing the best solutions among those consistent with the ...
[32]
How S/P Amplitude Ratio Data Can Bias Earthquake Focal ...
May 30, 2025 · In this study, I characterize several underappreciated issues with S/P amplitude ratio data that are relevant to focal mechanism inversion.
[33]
[PDF] A New Focal Mechanism Calculation Algorithm (REFOC) Using Inter ...
Small earthquakes' focal mechanisms are usually solved using P-wave first-motion polarities and/or S-/P-wave amplitude ratios, which are limited due to the low.
[34]
Focal mechanism determination using high-frequency waveform ...
We present a new method using high-frequency full waveform information to determine the focal mechanisms of small, local earthquakes monitored by a sparse ...
[35]
Focal Mechanism Determination using High Frequency, Full ...
In this research, we use high frequency waveform information to determine the focal mechanisms of small local earthquakes at an oil reservoir.
[36]
Full-waveform earthquake source inversion using simulation-based ...
SUMMARY. This paper presents a novel framework for full-waveform seismic source inversion using simulation-based inference (SBI). Traditional probabilistic.Missing: peer | Show results with:peer
[37]
Focal mechanism inversion of the 2018 MW7.1 Anchorage ...
In this study, high-rate GPS data and near-field broadband seismograms are used in separate and joint inversions by the generalized Cut-and-Paste (gCAP) method ...
[38]
Real-time determination of earthquake focal mechanism via deep ...
The authors here present a deep learning method to determine the source focal mechanism of earthquakes in realtime.Missing: transition | Show results with:transition
[39]
Focal mechanisms and moment magnitudes of micro-earthquakes in ...
This paper documents an investigation on the use of full waveform inversion to retrieve focal mechanisms of 11 micro-earthquakes (Mw 0.8 to 1.4).
[40]
Full-waveform event location and moment tensor inversion for ...
Events are located and focal mechanisms are derived via a multiparameter inversion, which minimizes the differences between synthetic and observed data. The ...
[41]
Real-time determination of earthquake focal mechanism via deep ...
Mar 4, 2021 · The proposed FMNet is a deep-learning-based intelligent algorithm that allows us to estimate the source focal mechanism using data from a ...Missing: graphical | Show results with:graphical
[42]
Towards End-to-End Earthquake Monitoring Using a Multitask Deep ...
Jun 7, 2025 · This extended model, PhaseNet+, simultaneously performs phase arrival-time picking, first-motion polarity determination, and phase association.
[43]
SKHASH: A Python Package for Computing Earthquake Focal ...
Apr 10, 2024 · We introduce a Python package for computing focal mechanism solutions. This algorithm, which we refer to as SKHASH, is largely based on the HASH algorithm ...Introduction · Method · New Features<|separator|>
[44]
Inversion for Focal Mechanisms Using Waveform Envelopes and ...
Nov 27, 2018 · This article investigates an alternative methodology that is more robust with respect to inappropriate velocity models: the inversion of ...
[45]
Focal Mechanisms... or "Beachballs" | U.S. Geological Survey
A focal mechanism, or "beachball", is a graphic symbol that indicates the type of slip that occurs during an earthquake: strike-slip, normal, thrust (reverse), ...
[46]
[PDF] A Primer on Focal Mechanism Solutions for Geologists
The complete charac- terization of an earthquake's focal mechanism provides important information, including the origin time, epicenter loca- tion, focal depth ...
[47]
[PDF] Earthquake Focal Mechanisms
Fault types and “Beach Ball” plots. USGS. Page 19. Example Focal mechanism diagrams on mid- ocean ridges. Stein and Wysession, An Introduction to seismology ...Missing: explanation | Show results with:explanation
[48]
A geometric setting for moment tensors - Oxford Academic
The seismic moment tensor is a mathematical description for earthquake sources. A moment tensor can represent a single earthquake, as in the case of a small ...
[49]
Moment Tensor - an overview | ScienceDirect Topics
The moment tensor is defined as a symmetrical 3 × 3 matrix that describes the source mechanisms of seismic events, including earthquakes and other seismic ...
[50]
Moment tensors - a practical guide - mXrap
Oct 7, 2025 · I will mostly discuss the 3D case; the 2D diagram is just a stereonet projection of the 3D beach ball. The construction of a beach ball diagram ...
[51]
Seismic Moment - an overview | ScienceDirect Topics
This is the definition of seismic moment tensor: It is the stress produced by the inelastic deformation of a body that is elastic everywhere. It should be clear ...
[52]
Seismic moment tensors and estimated uncertainties in southern ...
Mar 24, 2016 · Moment tensors are used to interpret the style of faulting and deformation in active tectonic settings, they are used within seismic wavefield ...Missing: formalism | Show results with:formalism
[53]
Non-double-couple earthquakes. 1. Theory
This paper summarizes the theory of “moment tensor” expansions of equivalent-force systems and analyzes many possible physical non-DC earthquake processes.
[54]
Apparent Non-Double-Couple Components as Artifacts of Moment ...
Apr 4, 2024 · Compilations of earthquake moment tensors from global and regional catalogs find pervasive non-double-couple (NDC) components with a mean deviation from a ...
[55]
Isotropic Source Components of Events in the 2019 Ridgecrest ...
Sep 13, 2021 · We examine source mechanisms of M ≥ 3.0 earthquakes in the 2019 M w 7.1 Ridgecrest sequence, using waveforms from 39 near-fault and regional stations.
[56]
Physical mechanisms for vertical‐CLVD earthquakes at active ...
Mar 4, 2013 · We identified 101 earthquakes with vertical compensated-linear-vector-dipole (vertical-CLVD) focal mechanisms that took place near active ...Introduction · Data and Analysis · Discussion of Physical... · Conclusions
[57]
Source mechanism of mining induced seismic events
Aug 1, 2008 · The double couple (DC) – a superposition of two single couples which balance their moments – is one of them. Being the body force equivalent of ...
[58]
Apparent non-double-couple components as artifacts of moment ...
Apr 4, 2024 · Compilations of earthquake moment tensors from global and regional catalogs find pervasive non-double-couple (NDC) components with a mean ...
[59]
Consistency of Non‐Double‐Couple Components of Seismic ...
Feb 16, 2022 · Earthquake moment tensors can be decomposed into double‐couple components describing slip on planar faults and non‐double‐couple (NDC) ...Origin of NDC Components · Moment Tensor Dataset · Generation of NDC...
[60]
Tectonic Regimes Inferred From Clustering of Focal Mechanisms ...
Dec 9, 2020 · Underlying patterns in focal mechanisms can be revealed by an approach that combines unsupervised learning techniques The prevailing ...
[61]
An improved method for determining the regional stress tensor using ...
Oct 10, 1984 · The method can be used to investigate possible stress inhomogeneities during earthquake sequences by analyzing subsets of the data population.
[62]
An improved method for the determination of the tectonic stress field ...
All conventional stress inversion methods, when applied to earthquake focal mechanism data, suffer from uncertainty as to which plane is the true fault plane.
[63]
Earthquakes focal mechanism and stress field pattern in the ...
The inversion technique scheme is used also in the present study for determining the regional stress field parameters for earthquake focal mechanism solutions ...
[64]
A refined methodology for stress inversions of earthquake focal ...
Nov 15, 2016 · We develop an improved methodology for reliable high-resolution inversions of focal mechanisms to background stress field orientation and stress ratio R in two ...
[65]
Inferences about the local stress field from focal mechanisms
Mar 3, 2017 · Focal mechanisms determined from regional-network earthquake data or aftershock field investigation often contain members ranging from ...
[66]
Influence of Focal Mechanism in Probabilistic Seismic Hazard Analysis
Mar 9, 2017 · An extension of probabilistic seismic hazard analysis is proposed to introduce a priori information about seismic source parameters.
[67]
[PDF] Probabilistic Seismic Hazard Assessment: Subduction zone sources ...
- Fault plane focal mechanism. -Estimate of the maximum magnitude that can occur on the fault. - Estimate of the mean recurrence interval of magnitudes on ...<|control11|><|separator|>
[68]
[PDF] Probabilistic Seismic Hazard Assessment for Jamaica
Probabilistic Seismic Hazard Assessment for Jamaica ... Table 2 Seismogenic zone characteristics. Zone Description. Type. Depth (km) Main focal mechanism.
[69]
An Enhanced Focal Mechanism Catalog of Induced Earthquakes in ...
May 8, 2025 · Earthquake focal mechanisms provide important information for analyzing subsurface fault zone geometries, rupture kinematics, and crustal stress ...
[70]
constraining focal mechanisms of injection induced seismicity using ...
We propose a new method to better constrain first motion focal mechanisms of fluid-induced seismicity by introducing in situ stress information estimated from ...SUMMARY · INTRODUCTION · CONCEPT AND... · APPLICATION TO FIELD DATA
[71]
Focal mechanisms of recent seismicity at Campi Flegrei, Italy
Dec 15, 2019 · Focal mechanisms of VT earthquakes at Campi Flegrei show predominance of normal kinematics on high angle faults.
[72]
Monitoring eruption activity using temporal stress changes at Mount ...
Feb 19, 2016 · Earthquake focal mechanisms represent the fault orientation and slip direction, which are influenced by the stress field. Focal mechanisms of ...
[73]
Seismicity monitoring by cluster analysis of moment tensors
Dec 31, 2013 · We suggest a new clustering approach to classify focal mechanisms from large moment tensor catalogues, with the purpose of automatically identify families of ...
[74]
Earthquake Focal Mechanisms as a Stress Meter of Active Volcanoes
Oct 1, 2022 · This study indicates that focal mechanisms can be used to monitor the stress state of an active volcano. Combining focal mechanisms with other ...
[75]
Inversion of earthquake focal mechanisms to obtain the ...
The determination of the seismotectonic regime based on the inversion of sets of focal mechanisms of earthquakes has been studied by several authors (e.g. ...
[76]
Resolving fault plane ambiguity for small earthquakes
Resolving fault-plane ambiguity from CMT solutions can also reduce the uncertainties in regional stress inversions based on earthquake focal mechanism data ( ...Missing: uniqueness | Show results with:uniqueness
[77]
Mitigating the effect of errors in source parameters on seismic ...
Errors in earthquake parameters can have a variety of causes. Data noise is an important factor and can be caused by oceanic microseisms, wind and local sources ...
[78]
Earthquake focal mechanisms with distributed acoustic sensing - PMC
Jul 13, 2023 · Earthquake focal mechanisms provide critical in-situ insights about the subsurface faulting geometry and stress state.Results · Improved Focal Mechanisms... · Methods
[79]
How S/P Amplitude Ratio Data Can Bias Earthquake Focal Mechanism Estimates | Seismological Research Letters | GeoScienceWorld
### Key Points on How S/P Amplitude Ratios Bias Focal Mechanism Inversions and Other Limitations
[80]
Non-Double-Couple Earthquake Mechanism as an Artifact of the ...
Apr 1, 2010 · A sizable amount of moderate and strong earthquakes exhibit a considerable percentage of non-double-couple (non-DC) components in the mechanism.Introduction · Velocity Profile Mismodeling · Discussion And ConclusionsMissing: debate | Show results with:debate
[81]
When are Non-Double-Couple Components of Seismic Moment ...
Mar 10, 2023 · Thus, constraining the isotropic component also reduces the appearance of spurious CLVD components.
[82]
Resolution of non-double-couple components in the seismic ...
Abstract. We perform a detailed synthetic study on the resolution of non-double-couple (non-DC) components in the seismic moment tensors from short-period.<|separator|>
[83]
Faults and Non‐Double‐Couple Components for Induced ...
Aug 6, 2018 · In this paper, we examine fault-related (double-couple (DC)) and possibly fluid-related (non-double-couple (non-DC)) mechanisms of induced ...
[84]
When are Non-Double-Couple Components of Seismic Moment ...
Mar 11, 2023 · We explore this question by comparing moment tensors (MTs) of earthquakes in three global catalogs, which use different inversion procedures.
[85]
The Impact of Source Time Function Complexity on Stress‐Drop ...
Aug 21, 2024 · We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different ...<|separator|>
[86]
Non-double-couple components of the 2024 Noto earthquake ...
Sep 25, 2025 · Moment tensors (MTs) are an essential tool for representing the source mechanisms of seismic events. Under the point-source approximation, ...
[87]
Non‐Double‐Couple Microearthquakes in the Focal Area of the ...
Feb 10, 2020 · Earthquakes with non-double-couple (NDC) components are indicators of geometric fault complexity, anisotropy, tensile failure, ...Introduction · NDC Microearthquakes · Modeling the NDC Events · DiscussionMissing: debate | Show results with:debate