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 1,045 records · Page 58Linked to original sources

USGS2016

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗

USGS2020

ShakeAlert Earthquake Early Warning System Performance During the 2019 Ridgecrest Earthquake Sequence

During July 2019, a sequence of earthquakes including a Mw6.4 foreshock and a Mw7.1 mainshock occurred near Ridgecrest, California. ShakeAlert, the U.S. Geological Survey (USGS) ShakeAlert public Earthquake Early Warning (EEW) system being developed for the U.S. West Coast, was operational during this time, though public alerting was only available within LA County. ShakeAlert created alert messages for many of the earthquakes, including the two largest events, and for many of the larger aftershocks. In this study, we dissect log files and replay data through the system to reconstruct the sequence of events and analyze the performance of the system during that time period. While the system performed reasonably well overall, the sequence also revealed various issues and short comings that will be addressed in impending and future system upgrades, with most parts of the system working as they should. ShakeAlert correctly detected and rapidly characterized both the Mw6.4 and Mw7.1 earthquakes within 6.9 s of their origin times and created alert messages that were available to ShakeAlert’s pilot users. No public alerts were sent out by the ShakeAlertLA cellphone app (the only publicly available alerting method at the time) because the predicted shaking for LA County was below the app’s alerting threshold of MMI 4.0. For the Mw6.4 event this was accurate. For the Mw7.1 event, public alerts for LA County were warranted, but ShakeAlert underpredicted the shaking levels because both the point-source and the finite-fault algorithms underestimated the magnitude of the earthquake by 0.8 units. A number of software and hardware issues that were responsible for the magnitude underestimation of the mainshock have been identified and will be addressed in future ShakeAlert releases. We also analyze the hypothetical alerting performance of ShakeAlert had public alerting been available throughout southern California with a lower alerting threshold of 〖MMI〗_alert=2.5MMI 2.5. We find that, despite the magnitude underestimation, ShakeAlert could have provided timely warnings to a large fraction of affected sites, including some of the near-epicentral sites with high ground motion intensities.

California↗

USGS2015

Hundreds of earthquakes per day: The 2014 Guthrie, Oklahoma, Earthquake Sequence

A remarkable increase in seismic activity in Oklahoma since 2009 has been shown to correlate closely with enhanced hydrocarbon extraction and associated wastewater disposal; 99% of this recent Oklahoma earthquake activity has occurred within 15 km of a call II injection well (Ellsworth, 2013). In response to this increase in seismic activity, the U.S. Geological Survey (USGS) partnered with the Oklahoma Geological Survey (OGS) to exchange waveform data from permanent and temporary seismic stations to improve the cataloging of earthquake source parameters for a broad region of north-central Oklahoma. For a particularly persistent earthquake sequence near Guthrie, Oklahoma, a subspace detection method is applied to data from nearby seismic stations. This approach documents the occurrence of hundreds of readily detectable, highly similar, earthquakes per day, with rates occasionally exceeding 1000 earthquakes per day. Time-varying changes in b-value appear episodic, suggesting a correlation with periods of reversible fault weakening and associated failure.

Oklahoman↗

USGS2017

A wideband magnetoresistive sensor for monitoring dynamic fault slip in laboratory fault friction experiments

A non-contact, wideband method of sensing dynamic fault slip in laboratory geophysical experiments employs an inexpensive magnetoresistive sensor, a small neodymium rare earth magnet, and user built application-specific wideband signal conditioning. The magnetoresistive sensor generates a voltage proportional to the changing angles of magnetic flux lines, generated by differential motion or rotation of the near-by magnet, through the sensor. The performance of an array of these sensors compares favorably to other conventional position sensing methods employed at multiple locations along a 2 m long × 0.4 m deep laboratory strike-slip fault. For these magnetoresistive sensors, the lack of resonance signals commonly encountered with cantilever-type position sensor mounting, the wide band response (DC to ≈ 100 kHz) that exceeds the capabilities of many traditional position sensors, and the small space required on the sample, make them attractive options for capturing high speed fault slip measurements in these laboratory experiments. An unanticipated observation of this study is the apparent sensitivity of this sensor to high frequency electomagnetic signals associated with fault rupture and (or) rupture propagation, which may offer new insights into the physics of earthquake faulting.

Sensors↗

USGS2014

Breaking the oceanic lithosphere of a subducting slab: the 2013 Khash, Iran earthquake

[1] Large intermediate depth, intraslab normal faulting earthquakes are a common, dangerous, but poorly understood phenomenon in subduction zones owing to a paucity of near field geophysical observations. Seismological and high quality geodetic observations of the 2013 M w 7.7 Khash, Iran earthquake reveal that at least half of the oceanic lithosphere, including the mantle and entire crust, ruptured in a single earthquake, confirming with unprecedented resolution that large earthquakes can nucleate in and rupture through the oceanic mantle. A rupture width of at least 55 km is required to explain both InSAR observations and teleseismic waveforms, with the majority of slip occurring in the oceanic mantle. Combining our well-constrained earthquake slip distributions with the causative fault orientation and geometry of the local subduction zone, we hypothesize that the Khash earthquake likely occurred as the combined result of slab bending forces and dehydration of hydrous minerals along a preexisting fault formed prior to subduction.

Geophysical Research Letters↗

USGS2004

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗

USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Holocene surface rupture history of an active forearc fault redefines seismic hazard in Southwestern British Columbia, Canada

Characterizing the hazard associated with Quaternary‐active faults in the forearc crust of the northern Cascadia subduction zone has proven challenging due to historically low rates of seismicity, late Quaternary glacial scouring, and dense vegetation that often obscures fault‐related geomorphic features. We couple lidar topography with paleoseismic trenching across the Leech River Fault on southern Vancouver Island to produce the first detailed surface rupture history of an onland forearc fault in British Columbia, Canada. The results indicate that this fault produced three surface‐rupturing earthquakes in the last ∼9 kyr and is therefore capable of producing large (Mw>6) earthquakes in the future. We provide new constraints on the fault's length (∼130 km) and Holocene slip rate (≥0.2–0.3 mm/year) that, together with the earthquake ages, should be incorporated into new seismic hazard assessments and building code practices relevant to urban centers in southwestern British Columbia (Canada) and northwestern Washington State (United States).

British Columbia↗

USGS1996

Frictional constraints on crustal faulting

We consider how variations in fault frictional properties affect the phenomenology of earthquake faulting. In particular, we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes. We model these variations using a rate- and state-dependent friction law, where we differentiate velocity-weakening behavior into two fields: the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening. Similarly, we differentiate velocity-strengthening behavior into two fields: the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening. The strong seismic field comprises the seismic slip concentrations, or asperities. The two “intermediate” fields, weak seismic and compliant, have frictional velocity dependences that are close to velocity neutral: these fields modulate both the tectonic loading and the dynamic rupture process. During the interseismic period, the weak seismic and compliant regions slip aseismically, while the strong seismic regions remain locked, evolving into stress concentrations that fail only in main shocks. The weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically. This “mixed” frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance. The model also provides a mechanism for rupture arrest: dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas, producing broad areas of accelerated afterslip. Aftershocks occur on both the weak seismic and compliant areas around a fault, but most of the stress is diffused through aseismic slip. Rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed. We test this frictional model by comparing the seismicity and the coseismic slip for the 1966 Parkfield, 1979 Coyote Lake, and 1984 Morgan Hill earthquakes. The interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip: these regions are interpreted as either weak seismic or compliant, depending on whether or not they manifest interevent seismicity.

Journal of Geophysical Research B: Solid Earth↗

USGS1995

San Andreas fault zone drilling project: scientific objectives and technological challenges

We are leading a new international initiative to conduct scientific drilling within the San Andreas fault zone at depths of up to 10 km. This project is motivated by the need to understand the physical and chemical processes operating within the fault zone and to answer fundamental questions about earthquake generation along major plate-boundary faults. Through a comprehensive program of coring, fluid sampling, downhole measurements, laboratory experimentation, and long-term monitoring, we hope to obtain critical information on the structure, composition, mechanical behavior and physical state of the San Andreas fault system at depths comparable to the nucleation zones of great earthquakes. The drilling, sampling and observational requirements needed to ensure the success of this project are stringent. These include: 1) drilling stable vertical holes to depths of about 9 km in fractured rock at temperatures of up to 300°C; 2) continuous coring and completion of inclined holes branched off these vertical boreholes to intersect the fault at depths of 3, 6, and 9 km; 3) conducting sophisticated borehole geophysical measurements and fluid/rock sampling at high temperatures and pressures; and 4) instrumenting some or all of these inclined core holes for continuous monitoring of earthquake activity, fluid pressure, deformation and other parameters for periods of up to several decades. For all of these tasks, because of the overpressured clay-rich formations anticipated within the fault zone at depth, we expect to encounter difficult drilling, coring and hole-completion conditions in the region of greatest scientific interest.

California↗

USGS2022

The Mw 5.1 9 August 2020, Sparta earthquake, North Carolina: The first documented seismic surface rupture in the eastern United States

At 8:07 a.m. EDT on 9 Aug. 2020 a Mw 5.1 earthquake located ~3 km south of Sparta, North Carolina, USA, shook much of the eastern United States, producing the first documented surface rupture due to faulting east of the New Madrid seismic zone. The co-seismic surface rupture was identified along a 2-km-long traceable zone of predominantly reverse displacement, with folding and flexure generating a scarp averaging 8–10-cm-high with a maximum observed height of ~25 cm. Widespread deformation south of the main surface rupture includes cm-dm–long and mm-cm–wide fissures. Two trenches excavated across the surface rupture reveal that this earthquake propagated to the surface along a preexisting structure in the shallow bedrock, which had not been previously identified as an active fault. Surface ruptures by faulting are rarely reported for M <6 earthquakes, and hence the Sparta earthquake provides an opportunity to improve seismic hazard knowledge associated with these moderate events. Furthermore, this earthquake occurred in a very low strain rate intraplate setting, where earthquake surface deformation, regardless of magnitude, is sparse in time and rare to observe and characterize.

North Carolina↗

USGS2007

Seismicity associated with the Sumatra-Andaman Islands earthquake of 26 December 2004

The U.S. Geological Survey/National Earthquake Information Center (USGS/ NEIC) had computed origins for 5000 earthquakes in the Sumatra-Andaman Islands region in the first 36 weeks after the Sumatra-Andaman Islands mainshock of 26 December 2004. The cataloging of earthquakes of mb (USGS) 5.1 and larger is essentially complete for the time period except for the first half-day following the 26 December mainshock, a period of about two hours following the Nias earthquake of 28 March 2005, and occasionally during the Andaman Sea swarm of 26-30 January 2005. Moderate and larger (mb ???5.5) aftershocks are absent from most of the deep interplate thrust faults of the segments of the Sumatra-Andaman Islands subduction zone on which the 26 December mainshock occurred, which probably reflects nearly complete release of elastic strain on the seismogenic interplate-thrust during the mainshock. An exceptional thrust-fault source offshore of Banda Aceh may represent a segment of the interplate thrust that was bypassed during the mainshock. The 26 December mainshock triggered a high level of aftershock activity near the axis of the Sunda trench and the leading edge of the overthrust Burma plate. Much near-trench activity is intraplate activity within the subducting plate, but some shallow-focus, near-trench, reverse-fault earthquakes may represent an unusual seismogenic release of interplate compressional stress near the tip of the overriding plate. The interplate-thrust Nias earthquake of 28 March 2005, in contrast to the 26 December aftershock sequence, was followed by many interplate-thrust aftershocks along the length of its inferred rupture zone.

Bulletin of the Seismological Society of America↗

USGS2013

Enriquillo–Plantain Garden fault zone in Jamaica: paleoseismology and seismic hazard

The countries of Jamaica, Haiti, and the Dominican Republic all straddle the Enriquillo–Plantain Garden fault zone ( EPGFZ), a major left-lateral, strike-slip fault system bounding the Caribbean and North American plates. Past large earthquakes that destroyed the capital cities of Kingston, Jamaica (1692, 1907), and Port-au-Prince, Haiti (1751, 1770), as well as the 2010 Haiti earthquake that killed more than 50,000 people, have heightened awareness of seismic hazards in the northern Caribbean. We present here new geomorphic and paleoseismic information bearing on the location and relative activity of the EPGFZ, which marks the plate boundary in Jamaica. Documentation of a river bank exposure and several trenches indicate that this fault is active and has the potential to cause major destructive earthquakes in Jamaica. The results suggest that the fault has not ruptured the surface in at least 500 yr and possibly as long as 28 ka. The long period of quiescence and subdued geomorphic expression of the EPGFZ indicates that it may only accommodate part of the ∼7–9 mm=yr plate deformation rate measured geodetically and that slip may be partitioned on other undocumented faults. Large uncertainties related to the neotectonic framework of Jamaica remain and more detailed fault characterization studies are necessary to accurately assess seismic hazards.

Bulletin of the Seismological Society of America↗

USGS2016

Kinematic ground motion simulations on rough faults including effects of 3D stochastic velocity perturbations

We describe a methodology for generating kinematic earthquake ruptures for use in 3D ground‐motion simulations over the 0–5 Hz frequency band. Our approach begins by specifying a spatially random slip distribution that has a roughly wavenumber‐squared fall‐off. Given a hypocenter, the rupture speed is specified to average about 75%–80% of the local shear wavespeed and the prescribed slip‐rate function has a Kostrov‐like shape with a fault‐averaged rise time that scales self‐similarly with the seismic moment. Both the rupture time and rise time include significant local perturbations across the fault surface specified by spatially random fields that are partially correlated with the underlying slip distribution. We represent velocity‐strengthening fault zones in the shallow (<5 km) and deep (>15 km) crust by decreasing rupture speed and increasing rise time in these regions. Additional refinements to this approach include the incorporation of geometric perturbations to the fault surface, 3D stochastic correlated perturbations to the P ‐ and S ‐wave velocity structure, and a damage zone surrounding the shallow fault surface characterized by a 30% reduction in seismic velocity. We demonstrate the approach using a suite of simulations for a hypothetical M w 6.45 strike‐slip earthquake embedded in a generalized hard‐rock velocity structure. The simulation results are compared with the median predictions from the 2014 Next Generation Attenuation‐West2 Project ground‐motion prediction equations and show very good agreement over the frequency band 0.1–5 Hz for distances out to 25 km from the fault. Additionally, the newly added features act to reduce the coherency of the radiated higher frequency ( f >1 Hz) ground motions, and homogenize radiation‐pattern effects in this same bandwidth, which move the simulations closer to the statistical characteristics of observed motions as illustrated by comparison with recordings from the 1979 Imperial Valley earthquake.

Bulletin of the Seismological Society of America↗

USGS2023

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗

USGS1998

Dynamic stress changes during earthquake rupture

We assess two competing dynamic interpretations that have been proposed for the short slip durations characteristic of kinematic earthquake models derived by inversion of earthquake waveform and geodetic data. The first interpretation would require a fault constitutive relationship in which rapid dynamic restrengthening of the fault surface occurs after passage of the rupture front, a hypothesized mechanical behavior that has been referred to as “self-healing.” The second interpretation would require sufficient spatial heterogeneity of stress drop to permit rapid equilibration of elastic stresses with the residual dynamic friction level, a condition we refer to as “geometrical constraint.” These interpretations imply contrasting predictions for the time dependence of the fault-plane shear stresses. We compare these predictions with dynamic shear stress changes for the 1992 Landers ( M 7.3), 1994 Northridge ( M 6.7), and 1995 Kobe ( M 6.9) earthquakes. Stress changes are computed from kinematic slip models of these earthquakes, using a finite-difference method. For each event, static stress drop is highly variable spatially, with high stress-drop patches embedded in a background of low, and largely negative, stress drop. The time histories of stress change show predominantly monotonic stress change after passage of the rupture front, settling to a residual level, without significant evidence for dynamic restrengthening. The stress change at the rupture front is usually gradual rather than abrupt, probably reflecting the limited resolution inherent in the underlying kinematic inversions. On the basis of this analysis, as well as recent similar results obtained independently for the Kobe and Morgan Hill earthquakes, we conclude that, at the present time, the self-healing hypothesis is unnecessary to explain earthquake kinematics.

California↗

USGS2011

Developing seismogenic source models based on geologic fault data

Calculating seismic hazard usually requires input that includes seismicity associated with known faults, historical earthquake catalogs, geodesy, and models of ground shaking. This paper will address the input generally derived from geologic studies that augment the short historical catalog to predict ground shaking at time scales of tens, hundreds, or thousands of years (e.g., SSHAC 1997). A seismogenic source model, terminology we adopt here for a fault source model, includes explicit three-dimensional faults deemed capable of generating ground motions of engineering significance within a specified time frame of interest. In tectonically active regions of the world, such as near plate boundaries, multiple seismic cycles span a few hundred to a few thousand years. In contrast, in less active regions hundreds of kilometers from the nearest plate boundary, seismic cycles generally are thousands to tens of thousands of years long. Therefore, one should include sources having both longer recurrence intervals and possibly older times of most recent rupture in less active regions of the world rather than restricting the model to include only Holocene faults (i.e., those with evidence of large-magnitude earthquakes in the past 11,500 years) as is the practice in tectonically active regions with high deformation rates. During the past 15 years, our institutions independently developed databases to characterize seismogenic sources based on geologic data at a national scale. Our goal here is to compare the content of these two publicly available seismogenic source models compiled for the primary purpose of supporting seismic hazard calculations by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) and the U.S. Geological Survey (USGS); hereinafter we refer to the two seismogenic source models as INGV and USGS, respectively. This comparison is timely because new initiatives are emerging to characterize seismogenic sources at the continental scale (e.g., SHARE in the Euro-Mediterranean, http://www.share-eu.org/; EMME in the Middle East, http://www.emme-gem.org/) and global scale (e.g., GEM, http://www.globalquakemodel.org/; Anonymous 2008). To some extent, each of these efforts is still trying to resolve the level of optimal detail required for this type of compilation. The comparison we provide defines a common standard for consideration by the international community for future regional and global seismogenic source models by identifying the necessary parameters that capture the essence of geological fault data in order to characterize seismogenic sources. In addition, we inform potential users of differences in our usage of common geological/seismological terms to avoid inappropriate use of the data in our models and provide guidance to convert the data from one model to the other (for detailed instructions, see the electronic supplement to this article). Applying our recommendations will permit probabilistic seismic hazard assessment codes to run seamlessly using either seismogenic source input. The USGS and INGV database schema compare well at a first-level inspection. Both databases contain a set of fields representing generalized fault three-dimensional geometry and additional fields that capture the essence of past earthquake occurrences. Nevertheless, there are important differences. When we further analyze supposedly comparable fields, many are defined differently. These differences would cause anomalous results in hazard prediction if one assumes the values are similarly defined. The data, however, can be made fully compatible using simple transformations.

Seismological Research Letters↗

USGS1977

A review of the effects of earthquakes on underground mines

The nature and geography of earthquakes is summarized and the effects of earthquakes on underground structures, due both to the causative process, faulting, and the secondary process, shaking, primarily in the continental United States, is reviewed. The principal sources of information are eyewitness reports, as instrumentally derived seismological data have been available for only about 70 years. Reference is made to the results of instrumented studies in mines and data pertaining to mine bumps and rock bursts in mines. The perceptibility and physical effects of earthquakes in mines and underground workings deserve more systematic study than they have received. Severe damage is inevitable when a mine or tunnel intersects a fault along which movement occurs during an earthquake. Mines in the epicentral region of strong earthquakes, but not transected by fault movement, may suffer severe damage by shaking. Mines outside of the epicentral region are likely to suffer little or no damage from a strong earthquake. Other factors being equal, it appears reasonable that the severity of damage due to shaking would probably be least in a mine located in fresh, highly competent rock; somewhat greater damage would probably be expected in a mine in weathered or less competent rock; and the greatest damage would be expected in a mine located in loose unconsolidated or incompetent rock.

Open-File Report↗