Abstract
In December 2007, an Mw2.2 earthquake occurred in a gabbroic dike at 3.3 km depth in a deep gold mine in South Africa. The fore- and aftershock activity was analyzed in an effort to understand the preparation and generation processes of earthquakes; these findings have already been published. The present paper focuses on the stress state in the source region of the mainshock. A 90-m-long borehole across the mainshock fault was drilled ~ 1.5 years after the mainshock and logged by an optical televiewer. The fault can be identified by severe damage to the borehole wall at the point where the borehole intersected the aftershock cluster. Except for a 10-m section in the hanging wall of the fault, borehole cores were fully recovered. Borehole breakout (BO) and core disking (CD) were found to occur. Two stress measurement techniques [Deformation Rate Analysis (DRA) and Diametrical Core Deformation Analysis (DCDA)] were applied to the borehole cores. By combining their results with occurrence criteria for BO and CD, the principal stress state in the source region of the mainshock was determined. The principal directions in the hanging wall of the fault were nearly identical to the virgin stress state, while it was significantly disturbed in the footwall. The vertical stresses were 106 MPa and 40 MPa in the hanging wall and footwall, respectively. The significant difference in the vertical stress between the two sides of the fault can be explained by the stress redistribution associated with a nonuniform slip on a nonplanar fault.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
To understand the mechanical properties of a fault, such as its strength and stability against a slip, investigations of the stress state around the fault are vitally important. For example, Zoback et al. (1987) estimated the stress state around the San Andreas fault, central California, and found that the orientation of the maximum horizontal compression (SHmax) is nearly normal to the fault. Sato et al. (2003) and Yamamoto and Yabe (2001) measured in situ stress near the Nojima fault, central Japan, on which the 1995 Southern Hyogo Prefecture Earthquake (M7.3) occurred. They showed that the shear stress acting on the fault was very small. These results suggest that some faults in the Earth’s crust are weak. On the other hand, Scholz (2000) claimed that the fault-normal compression in central California is a result of active folding in this area. He found that the SHmax direction in southern California has an angle of 30°–60° to the San Andreas fault, implying that the fault is strong. Yamashita et al. (2004) showed that the small shear stress on the Nojima fault was a result of the stress change associated with the coseismic slip. By correcting for the coseismic stress change, they found that the average ratio of the shear stress to the effective normal stress on the Nojima fault was ~ 0.6 before the earthquake. However, their estimate is strongly dependent on their assumptions regarding the pore pressure. Further, Yabe et al. (2008) detected localized stress perturbations close to both sides of the Chelungpu fault, central Taiwan, on which the 1999 Chi-Chi earthquake (Mw7.6) occurred. Thus, despite its importance, the issue of fault strength remains unresolved. Stress measurements around a fault are important, not only to determine the strength of the fault, but also to evaluate the likelihood of earthquake occurrence on the fault (e.g. Yabe & Omura, 2011; Yabe et al., 2010).
In deep gold mines in South Africa, a thin, tabular ore body (locally known as a “reef”) is excavated at depths larger than 3 km. Stresses around the mined-out stope are disturbed by several tens of MPa within a few hundred meters ahead of the stope (Chamber of Mines Research Organization, 1988). Nonlocal stress redistribution induced by relatively large-scale mining may cause relatively large earthquakes of up to M ~ 3 (strong enough to damage nearby mine workings) on preexisting weak planes or discontinuities, such as faults or dike contacts (Gay et al., 1984; Gibowicz & Kijko, 1994; Ziegler et al., 2015). The induced earthquakes occur within a few hundred meters of tunnels or stopes to which mine workers commute every day. Therefore, deep mines provide good opportunities to directly observe source faults of seismic events (e.g. Gay & Ortlepp, 1979; Heesakkers et al., 2011a, b; Mngadi et al., 2019). Further, it is possible to deploy an observation network very close to or within the source regions of the induced earthquakes (e.g. Ogasawara et al., 2002a, b; Reches et al., 2006). The closeness between the observation network and the seismic sources offers opportunities to record less distorted seismic signals (e.g. Yamada et al., 2005, 2007), and to reach the source fault itself with a short drill hole (e.g. Yabe et al., 2013).
A local seismic network (~ 40 m span) consisting of eight acoustic emission (AE) sensors and two strain meters was installed along a tunnel at 3.3 km depth in Mponeng mine (Fig. 1), which is located ~ 70 km west of Johannesburg, in June 2006 by the JAGUARS (Japanese-German Underground Acoustic Emission Research in South Africa) project (Nakatani et al., 2008). The tunnel was ~ 90 m below a gold reef that was planned to be mined in 2006–2007 and passed through a 30-m-thick gabbroic dike, called the Pink and Green dike or, simply, the PG dike. Because the dike does not contain gold, it was left unmined as a dip pillar. The overburden weight above the stope concentrated on the PG dike in the JAGUARS site.
On 27 December 2007, ~ 6 months after the deployment of the JAGUARS network, an Mw2.2 earthquake (mainshock) occurred in the JAGUARS site. The hypocenter of the mainshock was located in the PG dike and ~ 30 m above the JAGUARS network (Naoi et al., 2011; Yabe et al., 2009). Its foreshock and aftershock activity was investigated in detail by Yabe et al. (2009, 2015) and Naoi et al. (2011).
Hofmann et al. (2012) used a homogeneous boundary element model (BEM) to investigate whether this event was induced by stress transfer due to nearby mining. Ziegler et al. (2015) developed a more realistic three-dimensional finite element model (FEM) to evaluate the Coulomb failure stress changes (ΔCFS) on the mainshock fault associated with mining excavation. Both concluded that the Mw2.2 earthquake was induced by mining activity. They also showed that the mainshock fault was the optimally oriented plane under the stress state produced by mining. Since no stress measurement had been carried out in Mponeng mine, Hofmann et al. (2012) and Ziegler et al. (2015) adopted the far-field stress state in Tau Tona mine (neighboring Mponeng mine) estimated by Lucier et al. (2009) as the initial stress state of the models. Therefore, more direct investigation of the stress state in the source area was essential to validate their findings.
Yabe et al. (2013) drilled a 90-m-long hole (Y-hole, hereafter) with an inclination of 22° penetrating the source region of the mainshock ~ 1.5 years after its occurrence. The Y-hole intersected the source fault of the mainshock at ~ 8 m south of the hypocenter. The aftershock distribution along the Y-hole showed a concentration at ~ 70 m from the borehole collar (Fig. 2), suggesting that the borehole intersected the source fault there. The borehole core samples were recovered from the Y-hole, except for a 10-m section in the hanging wall on the source fault. Core disking (CD) and borehole breakout (BO) were found to occur intermittently. Yabe et al. (2013) were able to constrain the principal stress state around the source fault based on the occurrence criteria for CD (Kaga et al., 2003; Li & Schmitt, 1998; Matsuki et al., 2004) and BO (Barton et al., 1988; Zoback et al., 1985). The slip direction of the mainshock was also used as a constraint. The calculated principal directions were roughly consistent with the focal mechanism solution of the mainshock. However, estimation errors of the principal directions were large. Only the upper and the lower limit of the principal stress magnitudes could be estimated.
In the present study, we seek further constraints on the stress states. We apply Deformation Rate Analysis (DRA; Yamamoto, 2009; Yamamoto et al., 1990) and Diametrical Core Deformation Analysis (DCDA; Funato & Ito, 2017) techniques to the core samples. We then combine the DRA and DCDA measurements with the CD and BO occurrence criteria and the slip direction of the mainshock to determine the principal stress state in the source region more precisely.
2 Data and Methods
2.1 Overview of the Mainshock
Since the mainshock occurred during the end-of-year shutdown, more than 20,000 micro-earthquakes down to M ~ −4.4 (Kwiatek et al., 2010) were detected under low-noise environment by the JAGUARS network within the 150 h following the mainshock (Naoi et al., 2011). More than 13,000 of them formed a distinct planar cluster (~ 100 m × 80 m) on which the mainshock hypocenter lay. The aftershock plane has a strike subparallel to that of the PG dike and dip of ~ 60° to the northeast, cutting the dike diagonally (Fig. 1). Naoi et al. (2011) determined the CMT solution of the mainshock from seismic waveforms recorded by an ordinary seismic network operated by the mine. The mainshock had a normal fault mechanism, which is typical for mining-induced earthquakes in deep mines in South Africa. Because one of the nodal planes of the CMT solution coincided with the aftershock plane, the aftershock plane is thought to delineate the mainshock rupture plane.
Ogasawara et al. (2014) monitored the strain accumulation in the PG dike and in the host rock to the east of the dike from March 2007 using Ishii strain meters. The strain accumulation for ~ 9 months to the occurrence of the mainshock showed subvertical contraction and subhorizontal NE–SW extension in both the dike and the host rock. The amounts of the strain accumulation in the dike and the host rock were almost the same. This means that an increase in the subvertical compression associated with nearby mining was larger in the dike than in the host rock, because Young’s modulus of the dike is larger than that of the host rock (Table 1).
Based on the spatial and temporal evolution of the foreshock activity, Yabe et al. (2015) deduced the following mainshock nucleation process. The foreshocks concentrated in three clusters, which were interpreted to represent potential nucleation patches. Subsequent numerical modeling of the stress perturbation due to mining showed that two of the clusters were in a high shear stress zone on the mainshock fault (Hofmann et al., 2012). The foreshock activity in one of the nucleation patches in the high stress zone started to accelerate ~ 7 days before the occurrence of the mainshock, suggesting that this nucleation patch had reached a critical state. However, probably due to strength heterogeneity on the fault, the dynamic expansion of the critical patch was suppressed, although it continued to expand quasi-statically along the high stress zone. When the two nucleation patches in the high stress zone coalesced, the mainshock was finally initiated.
2.2 Stress Measurement Techniques
2.2.1 Diametrical Core Deformation Analysis (DCDA)
The DCDA method was proposed by Funato and Ito (2017). This technique measures the azimuthal variation in the diameter of a borehole core to estimate the differential stress in the plane normal to the borehole. The authors assumed that the borehole core expands elastically in response to stress relief when it is cut from rock mass by boring. The azimuthal variation in the diameter of the expanding core can be expressed by
where \(d(\theta )\) is the core diameter at an azimuth \(\theta\), \(d_{{\rm max}}\) and \(d_{{\rm min}}\) are the maximum and the minimum diameter of borehole core, respectively, and \(\alpha\) is the azimuth of \(d_{{\rm max}}\). The maximum and the minimum diameter of the core sample can be related to the differential stress in the plane normal to the borehole by
where \(\sigma_{x}\) and \(\sigma_{y}\) are the maximum and the minimum stress magnitudes, respectively, in the plane normal to the borehole. E and \({\nu }\) are Young’s modulus and Poisson’s ratio, respectively. d0 represents the diameter of the borehole core before it expands in response to the stress relief. Since it is difficult to know d0 in most cases, Funato and Ito (2017) assumed that \(d_{{\rm min}} - d_{0} \ll d_{0}\) and that \(\frac{{d_{{\rm max}} - d_{{\rm min}} }}{{d_{0} }}\) could be replaced by \(\frac{{d_{{\rm max}} - d_{{\rm min}} }}{{d_{{\rm min}} }}\).
Since borehole cores used in this study were not oriented, we measured the diameter of each borehole core every 2° from an arbitrary reference orientation along four or more profile lines separated by ~ 2 cm. The azimuthal variation in diameter can be disturbed for various reasons, such as dents on the core surface, inhomogeneous elasticity at a vein, and grooving during drilling. Therefore, consistency among phases of the azimuthal variations along profile lines on individual borehole cores was examined by the naked eye. If the out-of-phase variation was obtained along only one profile line, and other profile lines provided consistent diameter variations, data along the out-of-phase profile line were excluded from calculation of the average azimuthal variation. If the out-of-phase variations were found along two or more profile lines on a borehole core, the borehole core was not used for the DCDA.
The average azimuthal variation in each borehole core was calculated from the consistent azimuthal variations in its diameter. When a diameter on a profile line at an azimuth deviated by more than twice the standard deviation from the average diameter at the azimuth, the datum was regarded as an outlier. The average variation was then recalculated by excluding the outliers. The theoretical curve [Eq. (2)] was fitted to thus-obtained average azimuthal variation in core diameter to estimate \(d_{{\rm max}}\) and \(d_{{\rm min}}\).
2.2.2 Deformation Rate Analysis (DRA)
The DRA method was proposed by Yamamoto et al. (1990) and Yamamoto (2009) to recover the stress memory of a rock from hysteresis of the stress–strain curve under cyclic uniaxial loading. The hysteresis is evaluated by the strain difference function \(\Delta \varepsilon_{i, j} (\sigma )\) defined by
where \(\sigma\) is the applied axial stress and \(\varepsilon_{i} (\sigma )\) is the axial strain for the i-th loading. The stress memory is revealed by the bending of this function. In this study, five cyclic loading–unloading tests were conducted to obtain the strain difference functions. Further detail of the DRA is described in previous studies (Hsieh & Dight, 2016; Sato et al., 2003; Wang et al., 2012; Yabe & Omura, 2011; Yabe et al., 2008, 2010; Yamamoto, 2009; Yamamoto et al., 1990;). Further, de Andrade Penido (2021) demonstrated an advantage of simultaneous application of DCDA and DRA.
2.2.3 Borehole Breakout (BO) and Core Disking (CD)
When a borehole is drilled in a rock mass under compressional stresses, the compressional stresses are enhanced on the borehole wall due to stress concentration, while a tensile stress is induced ahead of the borehole bottom. If the enhanced compression is larger than the uniaxial compression strength (UCS) of the rock mass, borehole breakout (BO) takes place on the borehole wall in the direction of the minimum compression in the plane normal to the borehole (\(\sigma_{y}\)). In the case where the magnitude of tensile stress ahead of the borehole end reaches the uniaxial tensile strength (UTS) of the rock mass, core disking (CD) occurs. Therefore, the occurrence of BO and CD provides clues to constrain stresses, if the UCS and the UTS are known.
The occurrence criterion of BO is given by Zoback et al. (1985) and Barton et al. (1988) as the following:
where P is pore pressure. Because water is drained from rock mass in gold mines in South Africa before mining operations commence, it is assumed that \(P = 0\) in this study. The occurrence of BO was visually determined by inspection of the borehole wall by an optical televiewer. Because of the low resolution of the televiewer (Fig. 3a), the orientation of BO could be measured with precision of ~ 10°.
The occurrence criteria of CD are expressed by
where the coefficients A, B and D are functions of the disk thickness, \(\sigma_{z}\) is the normal stress parallel to the borehole and \(\sigma_{m}\) is the average of \(\sigma_{x} ,\; \sigma_{y}\) and \(\sigma_{z}\) (Kaga et al., 2003; Matsuki et al., 2004). In this study, the thicknesses of ~ 10 disks in each section of continuous occurrence of the CD were averaged to determine the coefficients A, B and D.
3 Results
Because the DCDA is a nondestructive and time-efficient technique, we used this technique as the first option to measure the differential stress in the plane normal to the Y-hole. When CD occurred, \(\sigma_{z}\) could be constrained from the occurrence criteria of CD. However, CD did not occur in the dike-footwall section, except at a calcite vein with a thickness of a few decimeters. We applied the DRA to borehole cores recovered from the dike-footwall section to evaluate \(\sigma_{z}\) in this section. The differential stress in the plane normal to the Y-hole in the dike-footwall section was also estimated by the DRA to confirm consistency between measurements by the DCDA and the DRA.
3.1 Borehole Breakout and Core Disking
BOs in the Y-hole and CDs of its cores were investigated by Yabe et al. (2013, Fig. 2). Their results are briefly summarized below.
The contact between the PG dike and the host rock was identified at a distance of 54.2 m from the Y-hole collar. Borehole cores were not recovered within a distance of 70–80 m (just after passing the fault). The UCS and the UTS were measured by uniaxial compression (test piece size was 14 mm in diameter and 36 mm in length) and Brazilian tests (42 or 27 mm in diameter and 21.5 or 16.5 mm in thickness), respectively. These strength values of the host rock, the dike-footwall and the dike-hanging wall, as well as their elastic moduli evaluated from the deformation during the UCS test, are summarized in Table 1 (see also Supplementary Information for data for each test piece). These values are used to constrain stresses based on the occurrence criteria of BO and CD.
Figure 4a shows the orientation of BOs and exfoliations with respect to the distance from the mainshock fault along the Y-hole. The orientation is measured counterclockwise from the horizontal by looking from the borehole collar. Exfoliation is revealed by continuous paired light damage or series of paired dents on the opposite sides of the borehole wall (Fig. 3b). Because the orientations of exfoliation are consistent with those of BO (Fig. 4a), we considered that the exfoliation is an initial stage of BO under a condition that the tangential stress concentration on the borehole wall is almost equal to the UCS of the rock mass. Therefore, the occurrence of the exfoliation may give tighter constraints on the stress magnitude than that of BO, because it can occur only when equality between both sides of Eq. (4) are satisfied. Further, the orientation of the exfoliation represents the orientation of \(\sigma_{y}\).
Along the Y-hole, BO and the exfoliation were found at several places over the host rock and the dike-hanging wall sections. In the dike-footwall section, BO and the exfoliation occurred at limited points near the fault and the contact. Except for proximity to the fault, the orientations of BO and the exfoliation are about −20° irrespective of the rock types. The orientation seems to rotate from 10° to −30° with an increase in the distance from the fault from −0.75 m (in the dike-footwall) to 3.4 m, (in the dike-hanging wall), probably reflecting such stress disturbance localized to the mainshock as reported by Yabe et al. (2008).
CD occurred on the borehole cores in the host rock section and in the dike-hanging wall section, while there was no CD in the dike-footwall section, except at the calcite veins. The average thickness of core disks in the dike-hanging wall section showed a systematic variation with distance from the fault (Fig. 4b). That is, it increases from 10.6 to 16.7 mm with an increase in the distance from the fault from 6 to 11.6 m. This suggests that the stresses rapidly decrease with an increase in the distance from the fault.
In the host rock section, it was observed that damages accumulate along evenly spaced annulus lines on the periphery of a core (Fig. 3c, d), while the core was not split. Li and Schmitt (1998) and Matsuki et al. (2004), for example, demonstrated that tensile stress around the borehole end induced by drilling is the largest at the periphery of the root of the core stub. The damage accumulation should represent nucleation of CD (tensile fracture along a plane normal to the borehole) or indication of CD under a condition that the tensile stress at the periphery of the root of the core stub is slightly less than the UTS. Therefore, the average separation of the damage accumulation lines is also plotted in Fig. 4b. Although they are within error ranges of each other, a tendency for the distance between damage accumulation lines to be greater than the thickness of CDs in the host rock is consistent with the abovementioned interpretation of the peripheral damage accumulation as the indication of CD.
3.2 Differential Stress Measured by Diametrical Core Deformation Analysis
Figure 5 shows examples of DCDA data in individual sections. DCDA was applied to 21 pieces of the borehole core, among which seven pieces were recovered from the host rock section and 13 pieces were from the dike-footwall section. Because the first 10 m of borehole cores in the dike-hanging wall section was severely damaged, only one piece could be measured by DCDA in this section. However, because nine pieces showed incoherent azimuthal variations in diameters, the \(d_{{\rm max}}\) and \(d_{{\rm min}}\) were obtained only from 12 pieces (four from the host rock, seven from the dike-footwall, one from the dike-hanging wall). When an anisotropy of rock mass from which a borehole core was recovered is significant, the differential stress estimation by DCDA is degraded. Figure 6 shows azimuthal variations in Young’s modulus and Poisson’s ratio in the plane normal to the Y-hole obtained from stress–strain curves during the UCS test. It should be noted that the azimuth is measured from an arbitrary reference, which may be different from the reference orientation of DCDA, since the Y-hole core is not oriented. When a sinusoidal function is fitted, the largest anisotropy of Young’s modulus and Poisson’s ratio is estimated to be 14% and 22%, respectively, in the dike-hanging wall. Our preliminary calculation (in preparation) indicates that an anisotropy of ~ 15% would cause ~ 10% error of the DCDA measurement. Further, variations in these elastic moduli at individual azimuths are similar to or larger than the magnitudes of anisotropy, implying that the estimated magnitudes of anisotropy are apparent. Therefore, we consider that an error of the differential stress estimated by DCDA is less than 10%.
Spatial variation in the differential stress determined by DCDA is shown in Fig. 4c with respect to the distance from the mainshock fault along the Y-hole. The differential stresses are 20–30 MPa in the host rock section. It is also 20–30 MPa near the contact in the dike-footwall, while it increases to ~ 100 ± 30 MPa in the dike-footwall section away from the contact. Twin densities of the calcite in the dike-footwall were counted to evaluate the maximum differential stress to be ~ 95 MPa at a point close to the fault (Janssen, personal communication; Rybacki et al., 2013). This is consistent with the differential stress measured by DCDA. The differential stress in the dike-hanging wall section is reduced to ~ 80 MPa, though this is within the error of the differential stress in the dike-footwall.
3.3 Stresses Measured by Deformation Rate Analysis
Normal stresses parallel to the Y-hole axis (\(\sigma_{z}\)) and in four directions (every 45°) in the plane normal to the Y-hole were measured by DRA using core samples recovered from the dike-footwall section at 59 m (Y59) and 68 m (Y68) distances from the hole collar of the Y-hole. Figure 7 shows examples of the strain difference functions obtained from three specimens of Y59 parallel to the Y-hole. Clear bending at 50–70 MPa is commonly seen on these strain difference functions. Though less clear than these cases, other strain difference functions also bend at similar stresses. By taking average of the bending stresses among the strain difference functions, the magnitude of \(\sigma_{z}\) at Y59 was determined to be 61.3 ± 6.4 MPa.
Figure 8 shows the strain difference functions obtained from specimens in the plane normal to the Y-hole at Y59. Azimuthal variations in normal stress in the plane normal to the Y-hole at Y59 and Y68 are plotted in Fig. 9. The magnitude of σz at Y59 is also shown in Fig. 9a. The azimuthal variation in normal stress can be expressed by \({\sigma } ( \theta ) = \sigma_{i} + 0.5\sigma_{d} {\text{cos}} ( {2 ( {\theta - \phi } )} )\), where \(\theta\) is azimuth from a reference direction and \(\phi\) is azimuth of the maximum compression. \(\sigma_{i} =0.5( \sigma_{x} + \sigma_{y})\) and \(\sigma_{d} = \sigma_{x} - \sigma_{y}\) are the isotropic and the differential stress, respectively. By fitting this function to the normal stresses determined by DRA, the isotropic and differential stresses at Y59 and Y68 are estimated to be 70–80 MPa and ~ 70 MPa, respectively. As shown in Fig. 4c, the differential stresses estimated by DRA and DCDA are consistent with each other, suggesting that the two techniques have similar reliability in the estimation of the differential stress.
3.4 Retrieved Principal Stress State
The principal stress state in the source region of the mainshock was estimated by the following grid-search procedure; (1) a set of three magnitudes of the principal stress (\(\sigma_{1} \ge \sigma_{2} \ge \sigma_{3}\)) and three rotation angles of the principal direction (\(\phi_{1}\), \(\phi_{2}\), \(\phi_{3}\); see Fig. 10 for the definition) is given. The given principal magnitudes were in a range from 0 to 300 MPa, at 5 MPa intervals. Because of their symmetry, the given principal directions were in a range from 0° to 180°, at 5° intervals. (2) Magnitudes of \(\sigma_{x}\), \(\sigma_{y}\) and \(\sigma_{z}\) and orientation of \(\sigma_{y}\) under the stress state given in (1) are calculated. (3) It is tested if these stresses are consistent with DCDA (\(\sigma_{x} - \sigma_{y}\)) and DRA (\(\sigma_{z}\), \(\sigma_{x} + \sigma_{y}\)) measurements within errors. Further, it is also examined if these stresses satisfy the occurrence criteria of BO and/or CD within errors of the strength measurements. The orientation of σy is compared with that of BO on the Y-hole wall to confirm that they are consistent within the BO orientation error of 20°. The coseismic slip direction of the mainshock is also considered. That is, misfit between direction of the maximum shear stress resolved on the mainshock fault and the slip direction should not be larger than 45°. (4) If the given principal stress state satisfies all conditions tested in (3), a unity score is added to individual bins of the principal stress magnitudes and those of the principal directions. (5) The principal stress magnitudes and the principal directions with the largest score are regarded as the most robust estimation of the principal stress, because a higher score of a parameter means that its adaptability to the conditions tested in (3) is less affected by values of the other parameters.
Stress measurement methods and conditions to be satisfied to constrain the principal stress in individual sections along the Y-hole are listed in Table 2. Errors given in Table 2 are set to twice the standard deviation of individual measurements. As BO did not occur in the dike-footwall section, the stress concentration on borehole wall (\(3\sigma_{x} - \sigma_{y}\)) in this section is required to be less than UCS. Strictly speaking, at a point where CD occurred, DCDA or DRA measurements cannot be applied. Therefore, the conditions to reproduce DCDA and DRA measurements and those to satisfy the occurrence criteria of CD are not necessarily required to be satisfied at the same time. However, for simplicity, it was assumed that the stress state is uniform in each section along the Y-hole and that all conditions should be satisfied at the same time.
Figure 11 shows the score distributions of individual parameters of the principal stress. Scores of all parameters in the host rock and the dike-hanging wall sections showed a single peak distribution. The most robust estimations could be uniquely determined in these sections. However, in the dike-footwall section, while the principal stress magnitudes and the third rotation angle (\(\phi_{3}\)) could be determined uniquely, scores of the first (\(\phi_{1}\)) and the second (\(\phi_{2}\)) rotations show two peaks. That is, \(\phi_{1}\) = 7.5° or 97.5°, and \(\phi_{2}\) = 27.5° or 112.5°. This may be because no information of \(\sigma_{y}\) orientation was obtained in this section. Four candidates of the principal axes are plotted on the lower hemisphere in Fig. 12. Background colors of this figure represent scores counted every 10° × 10° bins in azimuth and dip of the principal axis and normalized by its maximum for each axis. The normalized score of the maximum stress axis distributes in a broader region from (0, 30) to (270, 60) in azimuth (measured clockwise from north) and dip, while those of the intermediate and the minimum stress show a single peak at about (200, 0) and (100, 30), respectively. The combination of rotation angles with which the principal axes come close to the respective high score regions is (\(\phi_{1}\), \(\phi_{2}\), \(\phi_{3}\)) = (97.5°, 27.5°, 137.5°). Therefore, this combination of rotation angles is adopted as the robust estimation in the dike-footwall region.
The principal stress state thus obtained in each section is shown in Fig. 13, in which those within a half-width of the score distribution are also plotted as a proxy of the estimation errors.
The axes of \(\sigma_{1}\), \(\sigma_{2}\) and \(\sigma_{3}\) in the host rock section are well constrained to be about (318, 73), (119, 17) and (211, 5), respectively, in azimuth and dip. The \(\sigma_{3}\)-axis is almost parallel to the Y-hole. The most robust estimations of the principal stress magnitudes are \(\sigma_{1}\) = 102.5 MPa (98–110 MPa), \(\sigma_{2}\) = 87.5 MPa (75–96 MPa) and \(\sigma_{3}\) = 2.5 MPa (≤ 8 MPa). Here, the stress magnitude ranges in parentheses indicate the stress values where scores were larger than half the maximum score. For the minimum principal stress, since the largest score was at the lowest given magnitude (2.5 MPa), only its upper limit could be identified from score distribution.
The principal directions of \(\sigma_{1}\), \(\sigma_{2}\) and \(\sigma_{3}\) in the dike-footwall section are (353, 28), (240, 37) and (109, 41), respectively, being significantly different from those in the host rock and the dike-hanging wall sections. One of the most significant features is that dip of \(\sigma_{1}\)-axis is slightly smaller than those of \(\sigma_{2}\)- and \(\sigma_{3}\)-axes. This means that stress regime in this section is revers faulting, rather than the normal faulting regime that is typical in deep mines in South Africa. The principal stress magnitudes are \(\sigma_{1}\) = 97.5 MPa (93–109 MPa), \(\sigma_{2}\) = 52.5 MPa (20–66 MPa) and \(\sigma_{3}\) = 2.5 MPa (≤ 8 MPa).
The principal directions of \(\sigma_{1}\), \(\sigma_{2}\) and \(\sigma_{3}\) in the dike-hanging-wall section are (333, 58), (100, 21) and (200, 23), respectively. These are nearly identical to those in the host rock within error. However, the principal stress magnitudes are \(\sigma_{1}\) = 137.5 MPa (130–145 MPa), \(\sigma_{2}\) = 57.5 MPa (48–67 MPa) and \(\sigma_{3}\) = 2.5 MPa (≤ 17 MPa). That is, \(\sigma_{1}\) is ~ 30% larger and \(\sigma_{2}\) is ~ 30% smaller than those in the host rock.
4 Discussion
4.1 Relation to the Mw2.2 Earthquake
The P-, T- and B-axes of the CMT solution of the mainshock determined by Naoi et al. (2011) are projected on the lower hemisphere with the principal stress axes in the host rock section in Fig. 14a. P- and T-axes agree well with the maximum and the minimum stress axes, respectively, in this section. Since we used the slip direction on the fault calculated from the CMT solution to constrain the principal stress, this agreement could be artifact. To test this concern, we carried out the grid-search without the constraint condition of the slip direction. The obtained stress state is shown in Fig. 14b. While the error ranges of the principal stress axes become larger, both magnitudes and directions of the principal stress are almost unchanged. Therefore, the good agreement of P- and T-axes with \(\sigma_{1}\)- and \(\sigma_{2}\)-axes, respectively, is reliable. Shear stress along the slip direction of the mainshock relative to normal stress acting on the fault takes the largest scores in a range of 0.5–0.6, being consistent with a typical friction coefficient of dry rocks (Byerlee, 1978). This suggests that the mainshock occurred on the preferably oriented fault, as was pointed out by Ziegler et al. (2015). Further, it implies that our assumption that the fault was dry and \(P = 0\) in Eq. (4) was suitable.
Since BO did not occur in the dike-footwall section, the only constraint condition on the principal direction in this section is the slip direction of the mainshock. Figure 14c shows the score distribution of the \(\sigma_{1}\)-axis in this section obtained with the constraint of the slip direction. The high-score area is limited around the P-axis. When the constraint of the slip direction is omitted, the high-score area shows a zonal distribution (Fig. 14d). Therefore, the constraint of the slip direction plays an important role in constraining the principal direction in this section of the borehole where there is no information of stress orientation in the plane normal to the borehole.
4.2 Relation to Far-Field Stress State
Figure 15 compares the principal stress state in the host rock (quartzite) section determined in this study with those obtained by Lucier et al. (2009) and Ogasawara et al. (2013). Lucier et al. (2009) modeled the far-field stress state of Tau Tona mine (neighboring Mponeng mine) based on the distribution of borehole failures (BO and drilling induced tensile fracture) along multiple boreholes in quartzite. They concluded that the far-field stress state is a normal faulting regime as shown in Fig. 15b and that gradients of the principal stress magnitude are \(\sigma_{1}\) = 27 MPa/km, \(\sigma_{2}\) = 26 MPa/km and \(\sigma_{3}\) = 13 MPa/km. The principal stress magnitude in Fig. 15b is calculated from these gradients and depth of the JAGUARS site (3.3 km). Ogasawara et al. (2013) measured in situ stress at 3.4 km depth in Mponeng mine by a stress relief technique (Compact Conical-ended Borehole Overcoring, CCBO; Sugawara & Obara, 1999). As their measurement was carried out in quartzite at a site less affected by mining, their result can be proxy of the unaltered virgin stress state in this mine. The principal direction and magnitude of the virgin stress are shown in Fig. 15c. Considering estimation errors of ± 10–20°, the principal directions are consistent among these studies. However, the magnitudes of \(\sigma_{1}\) estimated in this study is ~ 10 MPa larger than those obtained by Lucier et al. (2009) and Ogasawara et al. (2013). On the other hand, the magnitude of \(\sigma_{3}\) estimated in this study is much smaller than that by Lucier et al. (2009) and Ogasawara et al. (2013). This is caused by the alteration of the stress state by mining and the earthquake, as discussed in the next section.
4.3 Alteration of Stress State by Mining and Faulting
4.3.1 Stress Transfer by Mining
As stated earlier, Ogasawara et al. (2014) observed subvertical contraction and subhorizontal NE–SW extension in both of the host rock and the PG dike for 9 months before the mainshock. While the amplitudes of these deformations in the host rock and the dike were almost equal to each other, stress changes were estimated to be larger in the dike because of a larger Young’s modulus.
As the \(\sigma_{1}\)-axis of the virgin stress state is subvertical, it should be enhanced by the increase in subvertical compression by mining. Further, the NE–SW orientation of the extension induced by mining is nearly parallel to the \(\sigma_{3}\)-axis of the virgin stress state. Hence, the magnitude of \(\sigma_{3}\) should be reduced. Thus, the larger magnitudes of \(\sigma_{1}\) and smaller magnitudes of \(\sigma_{3}\) estimated in the host rock and the dike-hanging wall sections than those of the virgin stress may be caused by the mining induced stress transfer. However, since the estimated stress state in the dike-footwall significantly differs from the virgin stress state, a stress alteration mechanism in addition to the stress transfer by mining is needed.
The vertical stress (qv in MPa) in gold mines in South Africa is expressed by a standard stress model (Jager & Ryder, 1999);
where h is the depth from ground surface in meters. As the depth of the JAGUARS site is ~ 3300 m, the expected vertical stress is 89 MPa. The vertical stresses in the host rock, the dike-footwall and the dike-hanging wall sections estimated in this study are \(\sigma_{{\rm v}}\) = 100 MPa, 40 MPa and 106 MPa, respectively, and 81 MPa on average. As uncertainty of our estimation is several or a dozen MPa, the average vertical stress is consistent with the standard model. The spatial variation in the vertical stress is induced partly by the stress transfer by mining, as stated above. However, since the distance from the Y-hole to the nearest mining area was ~ 90 m and both sides of the PG dike were excavated, such a large difference as 66 MPa between the vertical stress in the dike-footwall and the dike-hanging wall sections within a short distance of ~ 20 m is unlikely attributable to the stress transfer by mining. This significant contrast of the vertical stresses should be a result of the stress redistribution by the mainshock.
4.3.2 Stress Redistribution by the Mainshock
The aftershock distribution on a cross-section normal to the fault strike at the intersection of Y-hole with the mainshock fault delineates a nonplanar geometry of the fault (Fig. 16). While the aftershocks shallower then the mainshock hypocenter align along a steeply (82°) dipping plane with a width of ~ 35 m, those deeper than the hypocenter align along a moderately (58°) dipping plane with a width of ~ 25 m. As the intersection of Y-hole with the fault coincides with this bending of the fault, the non-planarity of the fault should asymmetrically disturb the stress states in the hanging wall and the footwall.
The seismic moment (\(M_{0}\)) is related to the moment magnitude (\(M_{w}\)) by \(\log M_{0} = 1.5 M_{w} + 9.1\). The average slip (d) on a fault can be estimated from the seismic moment by \(M_{0} = \mu Sd\), where \(\mu\) and \(S\) are the rigidity and the fault area, respectively. Taking into account a fault length of 100 m estimated from the aftershock distribution (Fig. 1b) and the rigidity of ~ 34 GPa calculated from Young’s modulus and the Poisson ratio of the PG dike (Yabe et al., 2013, Table 1), the average slip of the mainshock is estimated as ~ 10 mm when the slip is uniform over the two fault segments of the steeply and the moderately dipping planes.
To evaluate the stress redistribution in close proximity to the fault, we need detailed information regarding the coseismic slip distribution on the fault. However, seismic waveforms recorded by the mine-wide network were not of sufficiently high frequency to resolve the slip distribution of such a small earthquake as Mw2.2. Waveforms recorded by the JAGUARS network were saturated. Therefore, we assume that the slip on each of the two fault segments is rather uniform with a rake of −90° (normal fault). Further, for simplicity, only vertical stress disturbance is evaluated.
Hofmann et al. (2012) numerically evaluated a distribution of shear stress increases on the mainshock fault due to mining to suggest that slip in the upper half of the fault was much larger than that in the lower half. Therefore, we assumed three models of slips on the upper and the lower fault segments to be (a) 10 mm and 10 mm, (b) 12 mm and 9 mm and (c) 15 mm and 5 mm, respectively. The slips on the segments in individual models were adjusted to increase their difference between the segments from (a) to (c), but to conserve the seismic moment. Stress disturbance by the mainshock has been evaluated using Coulomb 3.3 (Lin & Stein, 2004; Toda et al., 2005), in which deformation of rock mass associated with slips on faults is calculated using Green’s functions of Okada (1992).
As shown in Fig. 16, all three models predict a decrease in the vertical stress by several tens of MPa in a region that the Y-hole passes in the footwall. This agrees with our estimations of vertical stress as small as \(\sigma_{{\rm v}}\) = 40 MPa in the dike-footwall. In contrast, the extent of a region with an increase in the vertical stress by more than a few MPa in the hanging wall shows significant variation among models. That is, the region where the vertical stress is significantly increased becomes wider, as the difference between slips on the upper and the lower fault segment is larger. The observation that the continuous occurrence of BO is limited to 10 m from the fault (Fig. 2) suggests that the extent of the region of significantly increased vertical stress is 10 m. Among the three models tested, model (b) shows a good agreement with this observation. These comparisons of the model prediction on the disturbance of the vertical stress and the observations imply that the reduction in the vertical stress in the dike-footwall was mostly caused by the non-planarity of the fault, while the significant enhancement of the vertical stress in the dike-hanging wall requires a nonuniform slip on the nonplanar fault.
5 Conclusions
Borehole cores were recovered from a borehole intersecting the source fault of an Mw2.2 earthquake that occurred in a gabbroic dike at ~ 3.3 km depth in a deep mine in South Africa. They were analyzed to constrain the principal stress state in the source region of the Mw2.2 earthquake. The source region was divided into three sections along the borehole from the hole collar: the host rock of the dike, the footwall in the dike (dike-footwall) and the hanging wall in the dike (dike-hanging wall). The differential stress in the plane normal to the borehole in each segment was measured by the DCDA method. The isotropic component of the stress in the plane normal to the borehole and the normal stress along the borehole were determined by the DRA method. These demonstrate that the DCDA and the DRA are useful techniques to measure stresses at a depth as great as ~ 3.3 km. Furthermore, they can contribute to the mitigation of the risks of rock bursting in deep mines.
In addition to DCDA and DRA, criteria for borehole breakout and core disking and the coseismic slip direction were taken into account to constrain the principal stress state in each section. The maximum stress directions in the host rock and the dike-hanging wall were subvertical, being almost identical to that in an area where the effect of mining was minimal. These principal directions agree well with the focal mechanism solution of the Mw2.2 earthquake. The average vertical stress among the three sections was consistent with that expected at the depth of studied site from a standard model for the mining district.
It was found that the vertical stresses in the footwall was reduced and that it was enhanced in the hanging wall of the fault of the Mw2.2 earthquake. The former seems attributable to the stress redistribution due to a coseismic slip on a nonplanar fault, which was delineated by the precisely determined aftershock distribution. The latter cannot be reproduced by a uniform slip on the nonplanar fault, but requires a nonuniform slip suggested by the numerically evaluated shear stress distribution, in addition to the non-planarity of the fault.
Data availability
Please contact the corresponding author.
Code availability
No special code was used.
References
Barton, C., Zoback, M., & Burns, K. (1988). In situ stress orientation and magnitude at the Fenton Geothermal Site, New Mexico, determined from wellbore breakouts. Geophysical Research Letter, 15, 467–470. https://doi.org/10.1029/GL015i005p00467
Byerlee, J. D. (1978). Friction of rocks. Pure and Applied Geophysics, 116, 615–626.
Chamber of Mines Research Organization (COMRO) (1988). An Industry Guide to /Methods of Ameliorating the Hazards of Rockfalls and Rockbursts (Chamber of Mines of South Africa).
de Andrade Penido, H., Funato, A., Metsugi, H., Torres, V. F. N., Sotomayor, J. M. G., Dight, P., de Figueiredo, R. P., de Assis, A. P., & Guimaraes, A. (2021). Application of the HF, DRA and DCDA technologies for in situ stress determination in Iron Quadrangle rock masses. Geomechanics and Geoengineering. https://doi.org/10.1080/17486025.2021.1955162
Funato, A., & Ito, T. (2017). A new method of diametrical core deformation analysis for in-situ stress measurements. International Journal of Rock Mechanics and Mining Sciences, 91, 112–118. https://doi.org/10.1016/j.ijrmms.2016.11.002
Gay, N., & Ortlepp, D. (1979). Anatomy of a mining-induced fault zone. Geological Society of America Bulletin, 90, 47–58. https://doi.org/10.1130/0016-7606(1979)90
Gay, N., Spencer, D., van Wyk, J., & van der Heever, P. (1984). The control of geological and mining parameters on seismicity in the Klerksdorp gold mining district. In N. C. Gay & E. H. Wainwright (Eds.), Proc. 1st International Congress on Rockbursts and Seismicity in Mines (pp. 107–120). The South African Institute of Mining and Metallurgy.
Gibowicz, S., & Kijko, A. (1994). An introduction to mining seismology (1st ed.). Academic Press.
Heesakkers, V., Murphy, S., Lockner, D. A., & Reches, Z. (2011a). Earthquake rupture at focal depth, Part I: Structure and rupture of the Pretorius fault, TauTona Mine, South Africa. Pure and Applied Geophysics, 168, 2395–2425. https://doi.org/10.1007/s00024-011-0354-7
Heesakkers, V., Murphy, S., Lockner, D. A., & Reches, Z. (2011b). Earthquake rupture at focal depth, Part II: Mechanism of the 2004 M2.2 earthquake along the Pretorius fault, TauTona mine, South Africa. Pure and Applied Geophysics, 168, 2427–2449. https://doi.org/10.1007/s00024-011-0355-6
Hofmann, G., Ogasawara, H., Katsura, T., & Roberts, D. (2012). An attempt to constrain the stress and strength of a dyke that accommodated a ML2.1 seismic event. Proceedings of the Second Southern Hemisphere International Rock Mechanics Symposium, 436–450.
Hsieh, A., & Dight, P. (2016). The desirable and undesirable effects on stress reconstruction using the deformation rate analysis (DRA). Paper presented at the ISRM International Symposium on In-Situ Rock Stress, Tampere, May 2016.
Jager, A., & Ryder, J. (1999). A handbook on rock engineering practice for tabular hard rock mines. The Safety in Mines Research Advisory Committee (SIMRAC).
Kaga, N., Matsuki, K., & Sakaguchi, K. (2003). The in situ stress states associated with core discing estimated by analysis of principal tensile stress. International Journal of Rock Mechanics and Mining Sciences, 40, 653–665. https://doi.org/10.1016/S1365-1609(03)00057-1
Kwiatek, G., Plenkers, K., Nakatani, M., Yabe, Y., Dresen, G., JAGUARS-Group. (2010). Frequency-magnitude characteristics down to magnitude −4.4 for induced seismicity recorded at Mponeng gold mine, South Africa. Bulletin of the Seismological Society of America, 100, 1165–1173. https://doi.org/10.1785/0120090277
Li, Y., & Schmitt, D. (1998). Drilling-induced core fractures and in situ stress. Journal of Geophysical Research, 103, 5225–5239. https://doi.org/10.1029/97JB02333
Lin, J., & Stein, R. S. (2004). Stress triggering in thrust and subduction earthquakes, and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults. Journal of Geophysical Research, 109, B02303. https://doi.org/10.1029/2003JB002607
Lucier, A. M., Zoback, M. D., Heesakkers, V., Reches, Z., & Murphy, S. K. (2009). Constraining the far-field in situ stress state near a deep South African gold mine. International Journal of Rock Mechanics and Mining Sciences, 46, 555–567. https://doi.org/10.1016/j.ijrmms.2008.09.005
Matsuki, K., Kaga, N., Yokoyama, T., & Tsuda, N. (2004). Determination of three dimensional in situ stress from core discing based on analysis of principal tensile stress. International Journal of Rock Mechanics and Mining Sciences, 41, 1167–1190. https://doi.org/10.1016/j.ijrmms.2004.05.002
Mngadi, S. B., Durrheim, R. J., Manzi, M. S. D., Ogasawara, H., Yabe, Y., Yilmaz, H., Wechsler, N., Van Aswegen, G., Roberts, D., Ward, A. K., Naoi, M., Moriya, H., Nakatani, M., Ishida, A., SATREPS Team, ICDP DSeis Team. (2019). Integration of underground mapping, petrology, and high-resolution microseismicity analysis to characterise weak geotechnical zones in deep South African gold mines. International Journal of Rock Mechanics and Mining Sciences, 114, 79–91. https://doi.org/10.1016/j.ijrmms.2018.10.003
Nakatani, M., Yabe, Y., Philipp, J., Morema, G., Stanchits, S., Dresen, G., JAGUARS-Group. (2008). Acoustic emission measurements in a deep gold mine in South Africa—Project overview and some typical waveforms. Seismological Research Letter, 79, 311.
Naoi, M., Nakatani, M., Yabe, Y., Kwiatek, G., Igarashi, T., & Plenkers, K. (2011). Twenty thousand aftershocks of a very small (M2) earthquake and their relation to the mainshock rupture and geological structures. Bulletin of the Seismological Society of America, 101, 2300–2407. https://doi.org/10.1785/0120100346
Ogasawara, H., Ando, M., & Yagitani, T. (2002a). Seismogenic process monitoring. Balkema.
Ogasawara, H., Kato, H., Hofmann, G., Roberts, D., Piper, P., Clements, T., Yabe, Y., & Sakaguchi, K. (2013). In-situ stress measurements to constrain stress and strength near seismic faults in deep level South African gold mines. Proceeding of 6th International Symposium on In-Situ Rock Stress, Sendai, 614–623.
Ogasawara, H., Katsura, T., Hofmann, G., Yabe, Y., Ishii, H., Roberts, D., Nakao, S., Okubo, M., Ward, A., & Kawakata, H. (2014). The Ishii strainmeters to in-situ monitor rock mass response to mining in South African gold mines. In The Southern African Institute of Mining and Metallurgy (Eds.), Proceedings of the 6th South African Rock Engineering Symposium, 21–33.
Ogasawara, H., Sato, S., Nishii, S., Sumitomo, N., Ishii, H., Iio, Y., Nakao, S., Ando, M., Takano, M., Nagai, N., Ohkura, T., Kawakata, H., Satoh, T., Kusunose, K., Cho, A., Mendecki, A. J., Cichowicz, A., Green, R. W. E., & Kataka, M. O. (2002b). Semi-controlled seismogenic experiments in South African deep gold mines. The Journal of the South African Institute of Mining and Metallurgy, 2002, 243–250.
Okada, Y. (1992). Internal deformation due to shear and tensile faults in a half-space. Bulletin of the Seismological Society of America, 82, 1018–1040.
Reches, Z., the DAFSAM and NELSAM teams. (2006). Building a natural earthquake laboratory at focal depth (DAFSAM-NELSAM project, South Afirca). Scientific Drilling, 3, 30–33. https://doi.org/10.2204/iodp.sd.3.06.2006.
Rybacki, E., Evans, B., Janssen, C., Wirth, R., & Dresen, G. (2013). Influence of stress, temperature, and strain on calcite twins constrained by deformation experiments. Tectonophysics, 601, 20–36. https://doi.org/10.1016/j.tecto.2013.04.021
Sato, N., Yabe, Y., Yamamoto, K., & Ito, H. (2003). In situ stresses near the Nojima fault estimated by deformation rate analysis. Zisin, 2(56), 157–169.
Scholz, C. (2000). Evidence for a strong San Andreas fault. Geology, 28, 163–166. https://doi.org/10.1130/0091-7613(2000)28
Sugawara, K., & Obara, Y. (1999). Draft ISRM suggested method for in-situ stress measurement using the compact conical-ended borehole overcoring (CCBO) technique. International Journal of Rock Mechanics and Mining Science, 36, 307–322.
Toda, S., Stein, R. S., Richards-Dinger, K., & Bozkurt, S. (2005). Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer. Journal of Geophysical Research, 110, B05S16. https://doi.org/10.1029/2004JB003415
Wang, H. J., Dyskin, A. V., Hsieh, A., & Dight, P. (2012). The mechanism of the deformation memory effect and the deformation rate analysis in layered rock in the low stress region. Computers and Geotechnics, 44, 83–92. https://doi.org/10.1080/17486025.2021.1955162
Yabe, Y., Nakatani, M., Naoi, M., Iida, T., Satoh, T., Durrheim, R., Hofmann, G., Roberts, D., Yilmaz, H., Morema, G., & Ogasawara, H. (2013). Estimation of the stress state around the fault source of a Mw 2.2 earthquake in a deep gold mine in South Africa based on borehole breakout and core discing. In T. Ito (Eds.), Proceedings of the 6th International Symposium on In-Situ Rock Stress, 604–613.
Yabe, Y., Nakatani, M., Naoi, M., Philipp, J., Janssen, C., Watanabe, T., Katsura, T., Kawakata, H., Georg, D., & Ogasawara, H. (2015). Nucleation process of an M2 earthquake in a deep gold mine in South Africa inferred from on-fault foreshock activity. Journal of Geophysical Research, 120, 5574–5594. https://doi.org/10.1002/2014JB011680
Yabe, Y., & Omura, K. (2011). In-situ stress at a site close proximity to the Gofukuji fault, central Japan, measured using drilling cores. Island Arc, 20, 160–173. https://doi.org/10.1111/j.1440-1738.2010.00756.x
Yabe, Y., Philipp, J., Nakatani, M., Morema, G., Naoi, M., Kawakata, H., Igarashi, T., Dresen, G., Ogasawara, H., JAGUARS-Group. (2009). Observation of numerous aftershocks of an Mw1.9 earthquake with an AE network installed in a deep gold mine in South Africa. Earth, Planets and Space, 61, e49–e52. https://doi.org/10.1186/BF03352963
Yabe, Y., Song, S., & Wang, C. (2008). In-situ stress at the northern portion of the Chelungpu fault, Taiwan, estimated on boring cores recovered from a 2-km-deep hole of TCDP. Earth, Planets and Space, 60, 809–819. https://doi.org/10.1186/BF03352832
Yabe, Y., Yamamoto, K., Sato, N., & Omura, K. (2010). Comparison of stress state around the Atera fault, central Japan, estimated using core samples and by improved hydraulic fracture tests. Earth, Planets and Space, 62, 257–268. https://doi.org/10.5047/eps.2009.12.002
Yamada, T., Mori, J. J., Ide, S., Abercrombie, R. E., Kawakata, H., Nakatani, M., Iio, Y., & Ogasawara, H. (2007). Stress drops and radiated seismic energies of microearthquakes in a South African gold mine. Journal of Geophysical Research, 112, B03305. https://doi.org/10.1029/2006JB004553
Yamada, T., Mori, J. J., Ide, S., Kawakata, H., Iio, Y., & Ogasawara, H. (2005). Radiation efficiency and apparent stress of small earthquake in a South African gold mine. Journal of Geophysical Research, 101, B01305. https://doi.org/10.1029/2004JB003221
Yamamoto, K. (2009). A theory of rock core-based methods for in-situ stress measurement. Earth, Planets and Space, 61, 1143–1161. https://doi.org/10.1186/BF03352966
Yamamoto, K., Kuwahara, Y., Kato, N., & Hirasawa, T. (1990). Deformation rate analysis: A new method for in situ stress estimation from inelastic deformation of rock samples under uniaxial compressions. Tohoku Geophysical Journal, 33, 127–147.
Yamamoto, K., & Yabe, Y. (2001). Stresses at sites close to the Nojima fault measured from core samples. Island Arc, 10, 266–281. https://doi.org/10.1111/j.1440-1738.2001.00325.x
Yamashita, F., Fukuyama, E., & Omura, K. (2004). Estimation of fault strength: Reconstruction of stress before the 1995 Kobe earthquake. Science, 306, 261–263. https://doi.org/10.1126/science.1101771
Ziegler, M., Reiter, K., Heidbach, O., Zang, A., Kwiatek, G., Stromeyer, D., Dahm, T., Dresen, G., & Hofmann, G. (2015). Mining-induced stress transfer and its relation to a Mw1.9 seismic event in an ultra-deep South African gold mine. Pure and Applied Geophysics, 172, 2557–3257. https://doi.org/10.1007/s00024-015-1033-x
Zoback, M., Moos, D., Mastin, L., & Andersson, R. (1985). Well bore breakouts and in-situ stress. Journal of Geophysical Research, 90, 5523–5530. https://doi.org/10.1029/JB090iB07p05523
Zoback, M., Zoback, M., Mount, V., Suppe, J., Eaton, J., Healy, J., Oppenheimer, D., Reasenberg, P., Jones, L., Raleigh, C., Wong, I., Scotti, O., & Wentworth, C. (1987). New evidence on the state of stress on the San- Andreas fault system. Science, 238, 1105–1111. https://doi.org/10.1126/science.238.4830.1105
Acknowledgements
The authors thank AngloGold Ashanti for kindly providing the research site in Mponeng mine. The authors thank an anonymous reviewer whose constructive comments and suggestions were useful to improve the manuscript. This study was supported by ICDP-DSeis project, JSPS-NRF bilateral research project, SATREPS, JSPS Core-to-Core program, Grant-in-Aid for Scientific Research, Earthquake Research Institute cooperative research program, Tohoku University’s 21st Century COE program, Ritsumeikan University, the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan, under its Earthquake and Volcano Hazards Observation and Research Program, and the DST-NRF South African Research Chair Initiative (SARChI).
Funding
ICDP-DSeis project, JSPS-NRF bilateral research project, SATREPS, JSPS Core-to-Core program, Grant-in-Aid for Scientific Research, Earthquake Research Institute cooperative research program, Tohoku University’s 21st Century COE program, Ritsumeikan University, the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan, under its Earthquake and Volcano Hazards Observation and Research Program, the DST-NRF South African Research Chair Initiative (SARChI).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
Not applicable.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yabe, Y., Abe, S., Hofmann, G. et al. Stress State in the Source Region of Mw2.2 Earthquake in a Deep Gold Mine in South Africa Determined from Borehole Cores. Pure Appl. Geophys. 179, 1679–1700 (2022). https://doi.org/10.1007/s00024-022-02999-w
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00024-022-02999-w