Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS

At least 361 records · Page 20Linked to original sources

USGS2009

The postseismic response to the 2002 M 7.9 Denali Fault earthquake: Constraints from InSAR 2003-2005

InSAR is particularly sensitive to vertical displacements, which can be important in distinguishing between mechanisms responsible for the postseismic response to large earthquakes (afterslip, viscoelastic relaxation). We produce maps of the surface displacements resulting from the postseismic response to the 2002 Denali Fault earthquake, using data from the Canadian Radarsat-1 satellite from the periods summer 2003, summer 2004 and summer 2005. A peak-to-trough signal of amplitude 4 cm in the satellite line of sight was observed between summer 2003 and summer 2004. By the period between summer 2004 and summer 2005, the displacement rate had dropped below the threshold required for observation with InSAR over a single year. The InSAR observations show that the principal postseismic relaxation process acted at a depth of ∼50 km, equivalent to the top of the mantle. However, the observations are still incapable of distinguishing between distributed (viscoelastic relaxation) and localized (afterslip) deformation. The imposed coseismic stresses are highest in the lower crust and, assuming a Maxwell rheology, a viscosity ratio of at least 5 between lower crust and upper mantle is required to explain the contrast in behaviour. The lowest misfits are produced by mixed models of viscoelastic relaxation in the mantle and shallow afterslip in the upper crust. Profiles perpendicular to the fault show significant asymmetry, which is consistent with differences in rheological structure across the fault.

Geophysical Journal International↗

USGS2003

Landslides and liquefaction triggered by the M 7.9 denali fault earthquake of 3 November 2002

The moment magnitude (M) 7.9 Denali Fault earthquake in Alaska of 3 November 2002 triggered an unusual pattern of landslides and liquefaction effects. The landslides were primarily rock falls and rock slides that ranged in volume from a few cubic meters to the 40 million-cubic-meter rock avalanche that covered much of the McGinnis Glacier. Landslides were concentrated in a narrow zone ???30 km wide that straddled the fault rupture zone over its entire 300 km length. Large rock avalanches all clustered at the western end of the rupture zone where acceleration levels are reported to have been the highest. Liquefaction effects, consisting of sand blows, lateral spreads, and settlement, were widespread within susceptible alluvial deposits extending from Fairbanks eastward several hundred kilometers. The liquefaction effects displayed a pattern of increasing concentration and severity from west to east and extended well beyond the zone of landslides, which is unusual. The contrasting patterns formed by the distributions of landslides and liquefaction effects initially seemed to be inconsistent; however, preliminary analyses of strong-motion records from the earthquake offer a possible explanation for the unusual ground-failure patterns that are related to three subevents that have been discerned from the earthquake records.

GSA Today↗

USGS2007

Stability and uncertainty of finite-fault slip inversions: Application to the 2004 Parkfield, California, earthquake

The 2004 Parkfield, California, earthquake is used to investigate stability and uncertainty aspects of the finite-fault slip inversion problem with different a priori model assumptions. We utilize records from 54 strong ground motion stations and 13 continuous, 1-Hz sampled, geodetic instruments. Two inversion procedures are compared: a linear least-squares subfault-based methodology and a nonlinear global search algorithm. These two methods encompass a wide range of the different approaches that have been used to solve the finite-fault slip inversion problem. For the Parkfield earthquake and the inversion of velocity or displacement waveforms, near-surface related site response (top 100 m, frequencies above 1 Hz) is shown to not significantly affect the solution. Results are also insensitive to selection of slip rate functions with similar duration and to subfault size if proper stabilizing constraints are used. The linear and nonlinear formulations yield consistent results when the same limitations in model parameters are in place and the same inversion norm is used. However, the solution is sensitive to the choice of inversion norm, the bounds on model parameters, such as rake and rupture velocity, and the size of the model fault plane. The geodetic data set for Parkfield gives a slip distribution different from that of the strong-motion data, which may be due to the spatial limitation of the geodetic stations and the bandlimited nature of the strong-motion data. Cross validation and the bootstrap method are used to set limits on the upper bound for rupture velocity and to derive mean slip models and standard deviations in model parameters. This analysis shows that slip on the northwestern half of the Parkfield rupture plane from the inversion of strong-motion data is model dependent and has a greater uncertainty than slip near the hypocenter.

California↗

USGS2008

Finite-fault analysis of the 2004 Parkfield, California, earthquake using P nl waveforms

We apply a kinematic finite-fault inversion scheme to P nl displacement waveforms recorded at 14 regional stations (Δ<2°) to recover the distribution of coseismic slip for the 2004 Parkfield earthquake using both synthetic Green’s functions (SGFs) calculated for one-dimensional (1D) crustal-velocity models and empirical Green’s functions (EGFs) based on the recordings of a single M w 5.0 aftershock. Slip is modeled on a rectangular fault subdivided into 2×2 km subfaults assuming a constant rupture velocity and a 0.5 sec rise time. A passband filter of 0.1–0.5 Hz is applied to both data and subfault responses prior to waveform inversion. The SGF inversions are performed such that the final seismic moment is consistent with the known magnitude ( M w 6.0) of the earthquake. For these runs, it is difficult to reproduce the entire P nl waveform due to inaccuracies in the assumed crustal structure. Also, the misfit between observed and predicted vertical waveforms is similar in character for different rupture velocities, indicating that neither the rupture velocity nor the exact position of slip sources along the fault can be uniquely identified. The pattern of coseismic slip, however, compares well with independent source models derived using other data types, indicating that the SGF inversion procedure provides a general first-order estimate of the 2004 Parkfield rupture using the vertical P nl records. The best-constrained slip model is obtained using the single-aftershock EGF approach. In this case, the waveforms are very well reproduced for both vertical and horizontal components, suggesting that the method provides a powerful tool for estimating the distribution of coseismic slip using the regional P nl waveforms. The inferred slip model shows a localized patch of high slip (55 cm peak) near the hypocenter and a larger slip area (~50 cm peak) extending between 6 and 20 km to the northwest.

California↗

USGS2011

Testing long-period ground-motion simulations of scenario earthquakes using the Mw 7.2 El Mayor-Cucapah mainshock: Evaluation of finite-fault rupture characterization and 3D seismic velocity models

Using a suite of five hypothetical finite-fault rupture models, we test the ability of long-period ( T >2.0 s) ground-motion simulations of scenario earthquakes to produce waveforms throughout southern California consistent with those recorded during the 4 April 2010 M w 7.2 El Mayor-Cucapah earthquake. The hypothetical ruptures are generated using the methodology proposed by Graves and Pitarka (2010) and require, as inputs, only a general description of the fault location and geometry, event magnitude, and hypocenter, as would be done for a scenario event. For each rupture model, two Southern California Earthquake Center three-dimensional community seismic velocity models (CVM-4m and CVM-H62) are used, resulting in a total of 10 ground-motion simulations, which we compare with recorded ground motions. While the details of the motions vary across the simulations, the median levels match the observed peak ground velocities reasonably well, with the standard deviation of the residuals generally within 50% of the median. Simulations with the CVM-4m model yield somewhat lower variance than those with the CVM-H62 model. Both models tend to overpredict motions in the San Diego region and underpredict motions in the Mojave desert. Within the greater Los Angeles basin, the CVM-4m model generally matches the level of observed motions, whereas the CVM-H62 model tends to overpredict the motions, particularly in the southern portion of the basin. The variance in the peak velocity residuals is lowest for a rupture that has significant shallow slip (<5 km depth), whereas the variance in the residuals is greatest for ruptures with large asperities below 10 km depth. Overall, these results are encouraging and provide confidence in the predictive capabilities of the simulation methodology, while also suggesting some regions in which the seismic velocity models may need improvement.

California↗

USGS2017

Long-term afterslip of the M6.0, 2004 Parkfield, California, earthquake—Implications for forecasting amount and duration of afterslip on other major creeping faults

We present the longest record of surface afterslip on a continental strike‐slip fault for the 2004 M 6.0 Parkfield, California, earthquake, from which we can derive critical information about the duration and predictability of afterslip relevant to urban displacement hazard applications. Surface slip associated with this event occurred entirely postseismically along the interseismically creeping (0.6–1.5 cm/yr) main trace of the San Andreas fault. Using the first year of afterslip data, the program AFTER correctly predicted the cumulative surface afterslip (maximum ∼35 cm) eventually attained. By 1 yr postearthquake, observed afterslip had accumulated to only ∼74% of its modeled final value u f in units of length. The 6‐yr data suggested final slip would be reached everywhere by ∼6–12 yrs. Parkfield’s afterslip lasted much longer (∼6–12 yrs) than afterslip following a 2014 M 6.0 event in Napa, California, where no interseismic creep was known, and its afterslip neared completion (∼97% of u f ) by 1 yr. The uncertainty in u f for the Napa event fell to ≤2 cm in only three months, versus in 2 yrs for the Parkfield event, mostly because duration of the power‐law stage of afterslip at Parkfield is much longer, ∼1000 (493–1666) days versus ∼100 (35–421) days for Napa. Because the urban Hayward fault near San Francisco, California, like the Parkfield section, exhibits interseismic creep in a similar geological regime, significant afterslip might last for up to a decade following an anticipated M ≥6.7 earthquake, potentially delaying postearthquake recovery.

California↗

USGS2021

The impact of 3D finite‐fault information on ground‐motion forecasting for earthquake early warning

We identify aspects of finite‐source parameterization that strongly affect the accuracy of estimated ground motion for earthquake early warning (EEW). EEW systems aim to alert users to impending shaking before it reaches them. The U.S. West Coast EEW system, ShakeAlert, currently uses two algorithms based on seismic data to characterize the earthquake’s location, magnitude, and origin time, treating it as a point or line source. From this information, ShakeAlert calculates shaking intensity and alerts locations where shaking estimates exceed a threshold. Several geodetic EEW algorithms under development would provide 3D finite‐fault information. We investigate conditions under which this information produces sufficiently better intensity estimates to potentially improve alerting. Using scenario crustal and subduction interface sources, we (1) identify the most influential source geometry parameters for an EEW algorithm’s shaking forecast, and (2) assess the intensity alert thresholds and magnitude ranges for which more detailed source characterization affects alert accuracy. We find that alert regions determined using 3D‐source representations of correct magnitude and faulting mechanism are generally more accurate than those obtained using line sources. If a line‐source representation is used and magnitude is calculated from the estimated length, then incorrect length estimates significantly degrade alert region accuracy. In detail, the value of 3D‐source characterization depends on the user’s chosen alert threshold, tectonic regime, and faulting style. For the suite of source models we tested, the error in shaking intensity introduced by incorrect geometry could reach levels comparable to the intrinsic uncertainty in ground‐motion calculations (e.g., 0.5–1.3 modified Mercalli intensity [MMI] units for MMI 4.5) but, especially for crustal sources, was often less. For subduction interface sources, 3D representations substantially improved alert area accuracy compared to line sources, and incorrect geometry parameters were more likely to cause error in calculated shaking intensity that exceeded uncertainties.

Bulletin of the Seismological Society of America↗

USGS2022

Topographic response to simulated Mw 6.5-7.0 earthquakes on the Seattle Fault

We explore the response of ground motions to topography during large crustal fault earthquakes by simulating several magnitude 6.5–7.0 rupture scenarios on the Seattle fault, Washington State. Kinematic simulations are run using a 3D spectral element code and a detailed seismic velocity model for the Puget Sound region. This model includes realistic surface topography and a near‐surface low‐velocity layer; a mesh spacing of ∼30 m at the surface allows modeling of ground motions up to 3 Hz. We simulate 20 earthquake scenarios using different slip distributions and hypocenter locations on a planar fault surface. Results indicate that average ground motions in simulations with and without topography are similar. However, shaking amplification is common at topographic highs, and more than a quarter of all sites experience short‐period (≤2 s) ground‐motion amplification greater than 25%–35%, compared with models without topography. Comparisons of peak ground velocity at the top and bottom of topographic features demonstrate that amplification is sensitive to period, with the greatest amplifications typically manifesting near a topographic feature’s estimated resonance frequency and along azimuths perpendicular to its primary axis of elongation. However, interevent variability in topographic response can be significant, particularly at shorter periods (<1 s). We do not observe a clear relationship between source centroid‐to‐site azimuths and the strength of topographic amplification. Overall, our results suggest that although topographic resonance does influence the average ground motions, other processes (e.g., localized focusing and scattering) also play a significant role in determining topographic response. However, the amount of consistent, significant amplification due to topography suggests that topographic effects should likely be considered in some capacity during seismic hazard studies.

Washington↗

USGS2025

Examining the role of elevated and sustained strain in dynamically triggering earthquakes on the Anza section of the San Jacinto fault

Microearthquakes can be dynamically triggered in southern California by remote earthquakes. However, directly connecting dynamic triggering mechanisms with observational data remains challenging. One proposed failure mechanism suggests that both the amplitude and duration of cyclic fatigue caused by the passing seismic wave contribute to triggering occurrence. Here, we measure dynamic strains recorded by borehole strainmeters in the Anza section of the San Jacinto fault zone from 710 earthquakes that occurred over 300 km away between 2008 and 2017 to systematically investigate the role of elevated and sustained strain in controlling dynamic triggering. We design a suite of tests to evaluate whether specific amplitude thresholds and durations of strain can predict dynamic triggering cases. We further test whether the peak dynamic strain (PDS) can predict triggering occurrence in combination with the strain amplitude and duration. Based on these tests, there is no strain amplitude–duration threshold that can distinguish triggering occurrence in Anza. Dynamic triggering is more likely to occur if a remote earthquake causes a PDS above 100 nanostrain, though many cases were triggered at smaller PDSs. The lack of clear correlation between triggering and characteristics of the dynamic strain field suggests that the tested features of the incoming waves do not determine triggering occurrence and local fault conditions and slip processes are more important in controlling dynamic triggering in Anza.

California↗

USGS2025

High-frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington↗

USGS2023

Shallow fault slip of the 2020 M5.1 Sparta, North Carolina, earthquake

The 2020 M 5.1 Sparta, North Carolina, earthquake is the largest in the eastern United States since the 2011 M 5.8 Mineral, Virginia, earthquake and produced a ∼2.5‐km‐long surface rupture, unusual for an event of this magnitude. A geological field study conducted soon after the event indicates oblique slip along a east‐southeast‐trending fault with a consistently observed thrust component. My analysis of regional seismic waveforms, Interferometric Synthetic Aperture Radar, and Global Positioning System survey data yields a compact shallow rupture extending from Earth’s surface down‐dip to the southwest over a ∼3 km fault length. The inferred kinematic rupture is primarily toward the up‐dip and eastward along‐strike directions and has predominantly thrust motion in the west, transitioning to roughly equal thrust and left‐lateral strike‐slip motion in the east. No normal faulting component, as proposed in an earlier geophysical study, is necessary to explain the data. The prevalence of only dip‐slip motions observed at Earth’s surface may demand slip partitioning between dip slip and lateral motions at depth.

North Carolina, Virginia, West Virginia↗

USGS2008

The Hayward Fault— Is it due for a repeat of the powerful 1868 earthquake?

On October 21, 1868, a magnitude 6.8 earthquake struck the San Francisco Bay region. Although the region was then sparsely populated, this quake on the Hayward Fault was one of the most destructive in California's history. Recent studies show that such powerful Hayward Fault quakes have repeatedly jolted the region in the past. U.S. Geological Survey (USGS) scientists describe this fault as a tectonic time bomb, due anytime for another magnitude 6.8 to 7.0 earthquake. Because such a quake could cause hundreds of deaths, leave thousands homeless, and devastate the region's economy, the USGS and other organizations are working together with new urgency to help prepare Bay Area communities for this certain future quake.

California↗

USGS1978

A study of fault plane mechanisms preceding the Thanksgiving Day, 1974 earthquake at Hollister, California

Fault plane solutions were obtained for 80 microearthquakes in a 900 square km area centered on the Busch fault during a one year period preceding the Thanksgiving Day, 1974 Hollister earthquake. Strike slip motion occurred in 65 of the quakes, mainly at active spots on the Sargent, Castro, Busch, and Calaveras faults. There are 11 thrust-type and 4 normal-type solutions. No discernable rotation of compressional axes was observed to precede the Thanksgiving Day earthquake.

California↗

USGS1995

Breaks in Pavement and Pipes as Indicators of Range-Front Faulting Resulting from the 1989 Loma Prieta Earthquake near the Southwest Margin of the Santa Clara Valley, California

Damage to pavement and near-surface utility pipes, caused by the October 17, 1989, Loma Prieta earthquake, provide indicators for ground deformation in a 663 km2 area near the southwest margin of the Santa Clara Valley, California. The spatial distribution of 1284 sites of such damage documents the extent and distribution of detectable ground deformation. Damage was concentrated in four zones, three of which are near previously mapped faults. The zone through Los Gatos showed the highest concentration of damage, as well as evidence for pre- and post-earthquake deformation. Damage along the foot of the Santa Cruz Mountains reflected shortening that is consistent with movement along reverse faults in the region and with the hypothesis that tectonic strain is distributed widely across numerous faults in the California Coast Ranges.

Open-File Report↗

USGS2012

Earthquakes and faults in southern California (1970-2010)

The map depicts both active and inactive faults and earthquakes magnitude 1.5 to 7.3 in southern California (1970–2010). The bathymetry was generated from digital files from the California Department of Fish And Game, Marine Region, Coastal Bathymetry Project. Elevation data are from the U.S. Geological Survey National Elevation Database. Landsat satellite image is from fourteen Landsat 5 Thematic Mapper scenes collected between 2009 and 2010. Fault data are reproduced with permission from 2006 California Geological Survey and U.S. Geological Survey data. The earthquake data are from the U.S. Geological Survey National Earthquake Information Center.

California↗

USGS2023

Using corrected and imputed polarity measurements to improve focal mechanisms in a regional earthquake catalog near the Mt. Lewis Fault Zone, California

We utilized relative polarity measurements and machine learning techniques to better resolve focal mechanisms and stress orientations considering a catalog of ∼29,000 relocated earthquakes that occurred during 1984–2021 in the southeastern San Francisco Bay Area. Earthquake focal mechanisms are commonly produced using P wave first motion polarities, which traditionally requires events to be well-recorded across a seismic network with good focal sphere coverage. We adapted recently developed approaches that are less dependent on high signal-to-noise records and exploit similar waveforms to produce relative polarity and amplitude measurements between earthquake pairs. These techniques were previously only applied on localized earthquake sequences, and we further developed these approaches so that they can be utilized for regional catalogs. We validated or corrected manually identified polarities by performing polarity consensuses using earthquake pairs. Missing and unreliable polarity measurements were imputed using iterative random forests, an unsupervised ensemble machine learning method. Relative P and S wave amplitude measurements were made between earthquakes, constraining S / P ratios for low signal-to-noise waveforms. Using these techniques, we were able to reduce focal mechanism uncertainties by an average of ∼13° and produced well-constrained focal mechanisms for ∼6 times as many earthquakes than those produced using only the traditionally derived polarities. We performed stress inversions using the focal mechanisms by grouping the focal mechanism results into a quadtree structure. Our stress results are consistent with previous work, albeit at a higher spatial resolution, and demonstrate these techniques can aid our understanding of fault structures and kinematics in more detail than was previously possible.

California↗

USGS2002

Paleoseismic event dating and the conditional probability of large earthquakes on the southern San Andreas fault, California

We introduce a quantitative approach to paleoearthquake dating and apply it to paleoseismic data from the Wrightwood and Pallett Creek sites on the southern San Andreas fault. We illustrate how stratigraphic ordering, sedimentological, and historical data can be used quantitatively in the process of estimating earthquake ages. Calibrated radiocarbon age distributions are used directly from layer dating through recurrence intervals and recurrence probability estimation. The method does not eliminate subjective judgements in event dating, but it does provide a means of systematically and objectively approaching the dating process. Date distributions for the most recent 14 events at Wrightwood are based on sample and contextual evidence in Fumal et al. (2002) and site context and slip history in Weldon et al. (2002). Pallett Creek event and dating descriptions are from published sources. For the five most recent events at Wrightwood, our results are consistent with previously published estimates, with generally comparable or narrower uncertainties. For Pallett Creek, our earthquake date estimates generally overlap with previous results but typically have broader uncertainties. Some event date estimates are very sensitive to details of data interpretation. The historical earthquake in 1857 ruptured the ground at both sites but is not constrained by radiocarbon data. Radiocarbon ages, peat accumulation rates, and historical constraints at Pallett Creek for event X yield a date estimate in the earliest 1800s and preclude a date in the late 1600s. This event is almost certainly the historical 1812 earthquake, as previously concluded by Sieh et al. (1989). This earthquake also produced ground deformation at Wrightwood. All events at Pallett Creek, except for event T, about A.D. 1360, and possibly event I, about A.D. 960, have corresponding events at Wrightwood with some overlap in age ranges. Event T falls during a period of low sedimentation at Wrightwood when conditions were not favorable for recording earthquake evidence. Previously proposed correlations of Pallett Creek X with Wrightwood W3 in the 1690s and Pallett Creek event V with W5 around 1480 (Fumal et al., 1993) appear unlikely after our dating reevaluation. Apparent internal inconsistencies among event, layer, and dating relationships around events R and V identify them as candidates for further investigation at the site. Conditional probabilities of earthquake recurrence were estimated using Poisson, lognormal, and empirical models. The presence of 12 or 13 events at Wrightwood during the same interval that 10 events are reported at Pallett Creek is reflected in mean recurrence intervals of 105 and 135 years, respectively. Average Poisson model 30-year conditional probabilities are about 20% at Pallett Creek and 25% at Wrightwood. The lognormal model conditional probabilities are somewhat higher, about 25% for Pallett Creek and 34% for Wrightwood. Lognormal variance ??ln estimates of 0.76 and 0.70, respectively, imply only weak time predictability. Conditional probabilities of 29% and 46%, respectively, were estimated for an empirical distribution derived from the data alone. Conditional probability uncertainties are dominated by the brevity of the event series; dating uncertainty contributes only secondarily. Wrightwood and Pallett Creek event chronologies both suggest variations in recurrence interval with time, hinting that some form of recurrence rate modulation may be at work, but formal testing shows that neither series is more ordered than might be produced by a Poisson process.

California↗

USGS2003

Detailed fault structure of the 2000 Western Tottori, Japan, earthquake sequence

We investigate the faulting process of the aftershock region of the 2000 western Tottori earthquake (Mw 6.6) by combining aftershock hypocenters and moment tensor solutions. Aftershock locations were precisely determined by the double difference method using P- and S-phase arrival data of the Japan Meteorological Agency unified catalog. By combining the relocated hypocenters and moment tensor solutions of aftershocks by broadband waveform inversion of FREESIA (F-net), we successfully resolved very detailed fault structures activated by the mainshock. The estimated fault model resolves 15 individual fault segments that are consistent with both aftershock distribution and focal mechanism solutions. Rupture in the mainshock was principally confined to the three fault elements in the southern half of the zone, which is also where the earliest aftershocks concentrate. With time, the northern part of the zone becomes activated, which is also reflected in the postseismic deformation field. From the stress tensor analysis of aftershock focal mechanisms, we found a rather uniform stress field in the aftershock region, although fault strikes were scattered. The maximum stress direction is N107??E, which is consistent with the tectonic stress field in this region. In the northern part of the fault, where no slip occurred during the mainshock but postseismic slip was observed, the maximum stress direction of N130??E was possible as an alternative solution of stress tensor inversion.

Tottori↗