Search USGS⌕ Search

USGS · 70022759

Energy budgets of mining-induced earthquakes and their interactions with nearby stopes

Abstract

In the early 1960's, N.G.W. Cook, using an underground network of geophones, demonstrated that most Witwatersrand tremors are closely associated with deep level gold mining operations. He also showed that the energy released by the closure of the tabular stopes at depths of the order of 2 km was more than sufficient to account for the mining-induced earthquakes. I report here updated versions of these two results based on more modern underground data acquired in the Witwatersrand gold fields. Firstly, an extensive suite of in situ stress data indicate that the ambient state of crustal stress here is close to the failure state in the absence of mining even though the tectonic setting is thoroughly stable. Mining initially stabilizes the rock mass by reducing the pore fluid pressure from its initial hydrostatic state to nearly zero. The extensive mine excavations, as Cook showed, concentrate the deviatoric stresses, in localized regions of the abutments, back into a failure state resulting in seismicity. Secondly, there appears to be two distinct types of mining-induced earthquakes: the first is strongly coupled to the mining and involves shear failure plus a coseismic volume reduction; the second type is not evidently coupled to any particular mine face, shows purely deviatoric failure, and is presumably caused by more regional changes in the state of stress due to mining. Thirdly, energy budgets for mining induced earthquakes of both types indicate that, of the available released energy, only a few per cent is radiated by the seismic waves with the majority being consumed in overcoming fault friction. Published by Elsevier Science Ltd.In the early 1960's, N.G.W. Cook, using an underground network of geophones, demonstrated that most Witwatersrand tremors are closely associated with deep level gold mining operations. He also showed that the energy released by the closure of the tabular stopes at depths of the order of 2 km was more than sufficient to account for the mining-induced earthquakes. I report here updated versions of these two results based on more modern underground data acquired in the Witwatersrand gold fields. Firstly, an extensive suite of in situ stress data indicate that the ambient state of crustal stress here is close to the failure state in the absence of mining even though the tectonic setting is thoroughly stable. Mining initially stabilizes the rock mass by reducing the pore fluid pressure from its initial hydrostatic state to nearly zero. The extensive mine excavations, as Cook showed, concentrate the deviatoric stresses, in localized regions of the abutments, back into a failure state resulting in seismicity. Secondly, there appears to be two distinct types of mining-induced earthquakes: the first is strongly coupled to the mining and involves shear failure plus a coseismic volume reduction; the second type is not evidently coupled to any particular mine face, shows purely deviatoric failure, and is presumably caused by more regional changes in the state of stress due to mining. Thirdly, energy budgets for mining induced earthquakes of both types indicate that, of the available released energy, only a few per cent is radiated by the seismic waves with the majority being consumed in overcoming fault friction.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Art McGarr. 2000. Energy budgets of mining-induced earthquakes and their interactions with nearby stopes. https://pubs.usgs.gov/publication/70022759

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Improving building footprint extraction using NAIP and 3DEP lidar derived features with deep learning

Accurate building footprint extraction is critical for applications ranging from population estimation to disaster management. Although optical imagery provides detailed spectral information, it often struggles with shadows, occlusions, and background clutter in dense urban environments. Lidar data, by contrast, offer precise elevation and structural attributes but face challenges such as variable point density and noise. This study integrates multispectral imagery from the U.S. Department of Agriculture (USDA) National Agriculture Imagery Program (NAIP) with lidar-derived feature height and intensity from the U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) to improve footprint extraction using a U-Net–based deep learning model. A six-band input stack (RGB, near-infrared, height, intensity) was developed, normalized, and tiled for training and evaluation against Microsoft Global Building Footprints (GBF). Results from the Houston, TX test site show that the six-band model achieved a precision of 0.86, recall of 0.88, F1 score of 0.87, and Intersection-over-Union (IoU) of 0.76, consistently outperforming four-band baselines by reducing false positives while maintaining sensitivity. Predictions on withheld Houston tiles confirmed strong within-region generalization, yielded a precision of 0.78, recall of 0.81, F1 score of 0.79, and IoU of 0.66. Qualitative analysis further revealed limitations stemming from both training label quality and vegetation–building confusion. These findings demonstrate the complementary value of integrating spectral and structural information for robust building footprint extraction and how domain adaptation strategies can be used to enhance cross-regional transferability.

Texas↗

Time-to-depth conversion of seismic-reflection data from eastern Lake Superior and implications for the eastern arm of the Midcontinent Rift

Seismic-reflection data were acquired in the mid 1980s along several lines across eastern Lake Superior by industry and the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) (Fig. 1). The lines form part of a larger network of crossing lines over the entire lake, which can be used to develop three-dimensional geologic models of the Mesoproterozoic Midcontinent Rift that lies below. To better interpret these lines, we developed velocity models to convert seismic reflections versus two-way travel time (TWTT) to reflections versus depth. In addition, the velocity models themselves provide insights into the structure of the Midcontinent Rift by recognizing common velocity ranges for certain rock types (Grauch, 2023).

eastern Lake Superior↗

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗