Induced stress and microseismicity in the 3000 Orebody, Mount Isa

Summary The 3000 Orebody is one of two orebodies in the Deep Copper Mine at Mount Isa, Australia. Owing to concerns about potential shaft pillar instabilities, an integrated seismic system was introduced to monitor seismic activity associated with pillar and country rock deformation. Coupled with numerical modelling of the stress regime, the system may assist in the characterization of rock mass damage resulting from mining, and perhaps the identification of near- and far-field geological structures that affect stope performance. A study was undertaken to quantify the seismicity and to determine potential applications of the seismic technology. The relation between geological structure and seismicity is strong, suggesting good prospects for the use of the system in the ground-control activities noted above. The induction of seismicity, which involves small magnitude events, is associated with reduction of normal stress on planes of weakness, suggesting that stress path may be an important factor in the level of seismicity observed in hard rock mines.     

The Belchatow, Poland, earthquakes of 1979 and 1980 induced by surface mining

Five seismic events occurred between August 1979 and May 1980 in the Belchatow trench area in central Poland, where large brown-coal deposits have been surface-mined since 1976. The three largest shocks had a local magnitudeM_L = 3.5 , seismic momentM_o = 1.5 · 10¹⁴ N-m and source radiusr = 350m , approximately. The tremors had maximum intensity between 5 and 6, and from the magnitude—intensity—depth relation, the focal depth was estimated to be between 1 and 5 km.

The Belchatow tremors are the manifestation of an unusual type of induced seismicity, i.e., seismicity connected with surface mining. The removal of overburden and extensive ground-water withdrawal by the mining operations seem to be direct factors responsible for stress concentration, although preexisting tectonic stress might be the most important indirect factor responsible for the origin of these events.     

Multi-scale imaging of a slow active fault zone: contribution for improved seismic hazard assessment in the Swiss Alpine foreland

Seismic hazard assessment of slow active fault zones is challenging as usually only a few decades of sparse instrumental seismic monitoring is available to characterize seismic activity. Tectonic features linked to the observed seismicity can be mapped by seismic imaging techniques and/or geomorphological and structural evidences. In this study, we investigate a seismic lineament located in the Swiss Alpine foreland, which was discussed in previous work as being related to crustal structures carrying in size the potential of a magnitude M 6 earthquake. New, low-magnitude (?2.0 ≤ M_L ≤ 2.5) earthquake data are used to image the spatial and temporal distribution of seismogenic features in the target area. Quantitative and qualitative analyses are applied to the waveform dataset to better constrain earthquakes distribution and source processes. Potential tectonic features responsible for the observed seismicity are modelled based on new reinterpretations of oil industry seismic profiles and recent field data in the study area. The earthquake and tectonic datasets are then integrated in a 3D model. Spatially, the seismicity correlates over 10–15 km with a N–S oriented sub-vertical fault zone imaged in seismic profiles in the Mesozoic cover units above a major decollement on top of the mechanically more rigid basement and seen in outcrops of Tertiary series east of the city of Fribourg. Observed earthquakes cluster at shallow depth (<4 km) in the sedimentary cover. Given the spatial extend of the observed seismicity, we infer the potential of a moderate size earthquake to be generated on the lineament. However, since the existence of along strike structures in the basement cannot be excluded, a maximum M 6 earthquake cannot be ruled out. Thus, the Fribourg Lineament constitutes a non-negligible source of seismic hazard in the Swiss Alpine foreland.     

Time-dependent analysis of aftershock events and structural impacts on intraplate crustal seismicity of the Van earthquake (Mw 7.1, 23 October 2011), E-Anatolia

The Van earthquake (M _W 7.1, 23 October 2011) in E-Anatolia is typical representative of intraplate earthquakes. Its thrust focal character and aftershock seismicity pattern indicate the most prominent type of compound earthquakes due to its multifractal dynamic complexity and uneven compressional nature, ever seen all over Turkey. Seismicity pattern of aftershocks appears to be invariably complex in its overall characteristics of aligned clustering events. The population and distribution of the aftershock events clearly exhibit spatial variability, clustering-declustering and intermittency, consistent with multifractal scaling. The sequential growth of events during time scale shows multifractal behavior of seismicity in the focal zone. The results indicate that the extensive heterogeneity and time-dependent strength are considered to generate distinct aftershock events. These factors have structural impacts on intraplate seismicity, suggesting multifractal and unstable nature of the Van event. Multifractal seismicity is controlled by complex evolution of crustal-scale faulting, mechanical heterogeneity and seismic deformation anisotropy. Overall seismicity pattern of aftershocks provides the mechanism for strain softening process to explain the principal thrusting event in the Van earthquake. Strain localization with fault weakening controls the seismic characterization of Van earthquake and contributes to explain the anomalous occurrence of aftershocks and intraplate nature of the Van earthquake.     

Tectonic deformation and subsurface active fault trends of El-Faiyum province,Egypt as deduced by seismicity and potential field data

We utilized the seismicity and the potential field data to study the tectonic deformation and to delineate the seismically active subsurface tectonic trends of El-Faiyum area. To accomplish these goals, we analyzed and interpreted the seismicity data, the reduced to pole total magnetic intensity, and the Bouguer anomaly maps. We also used the spatial distribution of the recent seismic events and the focal mechanism to outline the local seismic zones that control the seismicity of the study area and to determine the sense of the motion along the subsurface active faults. The focal mechanism of the recent seismic events and the interpreted subsurface tectonic faults from the potential field data reflect strike-slip movements with normal components along the subsurface active faults. The gravity and magnetic maps show a NE-SW regional trend with low gravity and magnetic values. The NE-SW regional trend extends across the whole area and could be related to the Pelusium Megashear fault. A NE-SW trend with high gravity and magnetic exists at the northern part and coincides with the Kattania Uplift and the basaltic flows in Gabal Qatrani area. The gravity and magnetic maps also reveal several local anomalies with different polarities, amplitudes, and extensions, which reflect anticlinal and synclinal structures on the basement surface. The seismotectonic map, generated by linking the basement structure map and the spatial distribution of the recent earthquake foci, reveals the dominant tectonic trends and the subsurface active faults.     

The Bishop's Castle (UK) Earthquake of 2 April 1990

A large earthquake, by British standards, occurred near Bishop's Castle in the Welsh Borders on 2 April 1990 at 13:46 GMT. This magnitude 5.1 ML event was felt over a wide area of Britain, from Ayrshire in the north to Cornwall in the south, Kent in the east and Dublin in the west. The epicentre was near the village of Clun, 7 km SSW of Bishop's Castle. Damage was minor and limited to the epicentral area, north to Wrexham and in particular Shrewsbury, which suffered most. Results from a macroseismic survey by BGS revealed that the maximum intensity in the epicentral area was 6 MSK. The mainshock had a focal depth of 14.3±4.7 km; however, better located aftershocks further constrained the mid-crustal seismicity to 15±0.2 km in the best cases. The marked lack of aftershocks contrasts with some previous similar magnitude events for intraplate earthquakes in Britain and throughout the world and may represent a large stress drop due to almost total relief of strain energy by the mainshock. The aftershock epicentral distribution shows a preference for an approximately N-S orientation which is consistent with one of the focal planes of the mainshock focal mechanism and suggests that this is the fault plane. Movement on this plane was predominantly strike-slip with a component of thrust and was consistent with a maximum compressive stress axis orientated NW-SE. The NE striking Welsh Borderland Fault System dominates the epicentral area; however, there is no surface fault which can clearly be related to the seismicity.     

Seismotectonics of Portugal and its adjacent area in the Atlantic

New elements on the seismicity of Portugal and new focal-mechanism solutions of earthquakes with epicentres situated off the coast of the Portuguese mainland and in the Azores region are presented. Historical seismicity data show that in the territory of the Portuguese mainland there are active faults that are responsible for earthquakes that have caused important damage and many casualties. However, most of the intraplate earthquakes with epicentres situated in the Portuguese mainland or near the shore are normally of small magnitude and this renders difficult their interpretation in the light of focal mechanisms. A solution for one earthquake, with magnitude 5 and epicentre at the Nazaré submarine canyon, is presented.Southwestwards of Cape St. Vincent there is an important seismic zone responsible for high-magnitude earthquakes such as that of 1 November 1755. This zone is situated in the region where the extension of the Messejana fault into the ocean joins with the Azores-Gibraltar fault.The seismicity of the area situated between the western coast of the Portuguese mainland and the Azores increases approximately along the 15°W meridian, from the latitude of the Azores-Gibraltar fault up to 44°N. Focal mechanisms of earthquakes with epicentres situated along this line show very similar solutions.The interpretation of the focal mechanism solutions of the earthquakes with epicentres situated in the studied area shows that the stress field trends approximately NW-SE. It is assumed that this stress field results from the interaction of the Eurasian and African plates; however, this direction is not maintained in the Azores region.     

Characteristics of three recent earthquake sequences in the Taupo Volcanic Zone, New Zealand

We have analysed three recent earthquake sequences in the northern part of the Taupo Volcanic Zone. A 1998 sequence at Haroharo with a largest event of M_L 4.8, and a 2004 sequence near Lake Rotoehu (largest event M_L 5.4), had normal b-values, and displayed an aftershock decay pattern, with most of the activity within the first few days. In contrast, a 2005 sequence a few tens of kilometres away at Matata (largest event M_L 4.1), had very different characteristics, with a slow development and decay, no tendency for enhanced seismicity after the larger events, and a very high b-value.The focal mechanisms of the Rotoehu and Matata events are normal, and have stress patterns consistent with the geodetically observed extension of the Taupo Volcanic Zone in a northwest–southeast direction. The extensive recent volcanism in the Okataina Volcanic Centre does not seem to have affected the stress pattern in this area.The Rotoehu sequence showed a strong resemblance, particularly in the time distribution of events, to the well-known swarm activity in the Vogtland region on the German/Czech border, in which larger events were followed by a burst of seismicity, as in a normal aftershock sequence. Some of the arguments that have been advanced to explain the Vogtland swarm as seismicity induced by fluid injection apply to Rotoehu, but there is no direct evidence of fluid involvement. The Matata sequence appears to have a continuing trigger mechanism, either a slow injection of fluid, or a slow slip event, in an environment in which opening pore spaces prevent high overpressures developing. The Matata sequence occurred close to the area of the 1987 M_L 6.3 Edgecumbe Earthquake, so exhibiting two extremes of seismic temporal pattern, namely mainshock–aftershock and a swarm with many events of similar magnitude, within a small area.     

The 2012/2014 swarmquakes of Jaen,Spain: a working hypothesis involving hydroseismicity associated with the hydrologic cycle and anthropogenic activity

A total of 2,300 swarmquake-type seisms have been striking the Torreperogil–Sabiote Hill (TSH) in southern Spain (Jaen) from October 2012 to the present, involving a major seismic crisis with six subseries until May 2013 (TSH seismic series, TSHS), followed by seven final series (FS). Two mild premonitory events occurred nearby in 2010–2011. Several evidences suggest cause–effect relationships between water and seismicity at the TSH: (1) The characteristics of seisms with a vast majority of ill-defined shallow/small events, spatial/temporal variations of the better constrained diffuse swarmquakes, increasing seismic velocities preceding stronger events, marked periodicities, abundant hydro-tremors, etc. (2) A series of spatial/temporal relationships between water and seismicity: Time lags between rainfall/reservoir peaks and seismic peaks; spatial links between swarmquakes and several factors (underground hydrogeology, reservoirs, etc.); etc. (3) From May 2013 to January 2014 evolution of rainfall and FS seismicity declined drastically at the same time. We suggest a hydroseismic hypothesis for TSH swarmquakes resulting from a complex series of cumulative factors: (1) A million-years natural setting involving geological, hydrogeological, geomorphological and tectonic elements. (2) A decades-long series of natural and anthropogenic factors are as follows: the 2012/2013 meteorological change, the massive overexploitation of the confined aquifers and alternating impoundment/emptying of the Giribaile reservoir.     

Seismicity and structural investigations of the Romanian Vrancea Region: Evidence for azimuthal variations of p-wave velocity and poisson's ratio

A study of the shallow and intermediate depth seismicity of the Romanian Vrancea region in the period 1964–1981 has been performed. The seismic events have been relocated by a standard location procedure using a regional velocity model. From the temporal and spatial distribution of the seismic activity, aspects of the seismicity before the large March 4, 1977 earthquake are treated, in particular the seismic gap in space and time prior to this event, found by Mârze (1979), which is critically discussed and revised. The concept of the precursor time/magnitude relationships of different authors is applied and its validity to the Vrancea region assessed. The hypocentral distribution shows that the intermediate depth seismic activity is confined to a small volume with dimensions of only some tens of kilometers. The results are interpreted in terms of the tectonics of the region. From an analysis of the travel-time residuals at different local stations, evidence for lateral velocity heterogeneities beneath the region is obtained e.g. a high velocity zone southeastwards of the Carpathian chain. Finally mean ratios, (i.e. Poisson's ratios), for various stations are calculated from P- and S-wave travel times. They show azimuthal variations of up to 6% for stations within the area where the intermediate seismic activity occurs in comparison with the station Focsani, situated eastwards in the Carpathian foredeeps. All these results are compatible with the plate tectonic concept for the Vrancea region, that is the subduction of an oceanic lithospheric slab under the Carpathian mountain arc, giving rise to such a highly active seismic zone.     

A Waveform Similarity Approach to Investigate Seismicity Patterns

the temporal and spatial evolution of seismicity, a simple cross-correlation technique was applied to a waveform data set of western Alps earthquakes. A selection of events recorded by a reference station of the regional seismic network of North-Western Italy (IGG network) within a particular S-P range was used. The adopted method was tested on 380 events recorded in 1995, allowing for the definition of the critical parameters for an optimal identification of earthquake families.The comparison of the individuated families definition with the relevant epicentres map confirms that, on small scales, the epicentre positions of microearthquakes located by a regional network could be strongly influenced by location errors. The analysis was extended to a more complete data set of 2653 events; 1171 of them were grouped in families. The results demonstrated that our single station method is able to discriminate between closely spaced families on the basis of the position and of the radiation properties.     

Moment tensors and micromechanical models

A numerical modelling approach that simulates cracking and failure in rock and the associated seismicity is presented and a technique is described for quantifying the seismic source mechanisms of the modelled events. The modelling approach represents rock as an assemblage of circular particles bonded together at points of contact. The connecting bonds can break under applied stress forming cracks and fractures in the modelled rock. If numerical damping is set to reproduce realistic levels of attenuation, then energy is released when the bonds break and seismic source information can be obtained as damage occurs. A technique is described by which moment tensors and moment magnitudes can be calculated for these simulated seismic events. The technique basically involves integrating around the source and summing the components of force change at the surrounding particle contacts to obtain the elements of the moment tensor matrix. The moment magnitude is then calculated from the eigenvalues of the moment tensor. The modelling approach is tested by simulating a well-controlled experiment in which a tunnel is excavated in highly stressed granite while microseismicity is recorded. The seismicity produced by the model is compared to the actual recorded seismicity underground. The model reproduces the spatial and temporal distribution of seismicity observed around the tunnel and also the magnitudes of the events. A direct comparison between the actual and simulated moment tensors is not performed due to the two-dimensional nature of the model, however, qualitative comparisons are presented and it is shown that the model produces intuitively realistic source mechanisms. The ability to obtain seismic source information from the models provides a unique means for model validation through comparison with actual recorded seismicity. Once it is established that the model is performing in a realistic manner, it can then be used to examine the micromechanics of cracking, failure and the associated seismicity and to help resolve the non-uniqueness of the geophysical interpretation. This is demonstrated by examining in detail the mechanics of one of the modelled seismic events by observation of the time dependence of the moment tensor and by direct examination of the particle motions at the simulated source.     

福建及其邻近地区地震活动节律性研究

本文通过对福建及其邻近地区五百多年来中强地震序列的研究,结果表明,该地区地震活动存在着明显的节律性,具有在节点上震级最大的特征,从而可预测该区下一次可能发生的中强地震的震级范围。     

The Woods Reef (New South Wales) earthquake of 14 November 1990: Focal mechanism derived from amplitude ratios and synthetic seismograms

On the morning of 15 November 1990 local time, Armidale and the area to the west of Armidale was shaken by a magnitude 3.2 earthquake. The epicentre was located at 30.39° S, 150.88° E and the depth of focus at 12 ± 7 km. As the epicentre was close to the Peel Fault an attempt was made to constrain the focal mechanism of this earthquake. The conventional method, which is based on the analysis of P wave polarities, was not applicable because the event was not strong enough. In an alternative method, the amplitudes of various seismic phases recorded at a number of stations well distributed in azimuth were compared with theoretical amplitudes calculated with the reflectivity method for a point shear dislocation in a layered medium. The differences between observed and calculated amplitudes were minimized as a function of fault strike, fault dip and direction of the slip vector. The analysis indicates that none of the possible fault planes had the strike of the Peel Fault. The solution suggests predominantly strike slip motion along two possible, steeply dipping fault planes. The inferred direction of the maximum compressional stress. is east‐west which is in good agreement with other estimates of the stress field for eastern Australia.     

Seismic activity at the western Pyrenean edge

The present-day seismicity at the westernmost part of the Pyrenean domain reported from permanent networks is of low to moderate magnitude. However, it is poorly constrained due to the scarce station coverage of the area. We present new seismic data collected from a temporary network deployed there for 17 months that provides an enhanced image of the seismic activity and its tectonic implications. Our results delineate the westward continuity of the E–W Pyrenean band of seismicity, through the Variscan Basque Massifs along the Leiza Fault, ending up at the Hendaya Fault. This seismicity belt is distributed on a crustal scale, dipping northward to almost 30 km depth. Other relevant seismic events located in the area can be related to the central segment of the Pamplona fault, and to different E–W thrust structures.     

Double-sided subduction with contrasting polarities beneath the Pamir-Hindu Kush: Evidence from focal mechanism solutions and stress field inversion

The Pamir-Hindu Kush region at the western end of the Himalayan-Tibet orogen is one of the most active regions on the globe with strong seismicity and deformation and provides a window to evaluate continental collision linked to two intra-continental subduction zones with different polarities. The seismicity and seismic tomography data show a steep northward subducting slab beneath the Hindu Kush and southward subducting slab under the Pamir. Here, we collect seismic catalogue with 3988 earthquake events to compute seismicity images and waveform data from 926 earthquake events to invert focal mechanism solutions and stress field with a view to characterize the subducting slabs under the Pamir-Hindu Kush region. Our results define two distinct seismic zones: a steep one beneath the Hindu Kush and a broad one beneath the Pamir. Deep and intermediate-depth earthquakes are mainly distributed in the Hindu Kush region which is controlled by thrust faulting, whereas the Pamir is dominated by strike-slip stress regime with shallow and intermediate-depth earthquakes. The area where the maximum principal stress axis is vertical in the southern Pamir corresponds to the location of a high-conductivity low-velocity region that contributes to the seismogenic processes in this region. We interpret the two distinct seismic zones to represent a double-sided subduction system where the Hindu Kush zone represents the northward subduction of the Indian plate, and the Pamir zone shows southward subduction of the Eurasian plate. A transition fault is inferred in the region between the Hindu Kush and the Pamir which regulates the opposing directions of motion of the Indian and Eurasian plates.     

Global-scale analysis of seismic activity prior to 2004 Sumatra-Andaman mega-earthquake

Following our experience in intermediate-term monitoring of seismicity before large earthquakes worldwide and in some regions we apply a similar approach to the analysis of seismicity on global scale several decades before the December 26, 2004 Sumatra–Andaman mega-earthquake. The Earth lithosphere is being considered as a single whole, representing the ultimate level of the complex Earth's hierarchy. The study brings up the following questions: Are there any anomalies of the global seismicity behaviour observed in advance of the mega-earthquake? Do these anomalies correspond to the global scale seismic patterns similar to those detected on the regional scale before great, major, and strong events? In other words, does the Earth lithosphere considered as a single whole show up an approach of the mega-earthquake in the way typical for events on the smaller levels of hierarchy? The results of our investigation favour the positive answer to the questions. Specifically we found that during the decade before the Sumatra–Andaman mega-earthquake the Earth lithosphere reveals, at least in intermediate-term scale, classical symptoms of instability, which can be depicted by known precursory seismic patterns. These are: (i) transformation of frequency–magnitude distribution, (ii) change in the rate of seismic activity, and (iii) depth redistribution of activity. Moreover changes of dependencies between magnitudes of different types are detected at the same time. The observed global scale patterns of collective behaviour of seismicity may indicate the state of criticality of the Earth lithosphere before the Sumatra–Andaman mega-earthquake.     

A new synthesis of Pyrenean seismicity and its tectonic implications

The Pyrenean range, which results from the convergence of the Iberian and Eurasian plates along the North Pyrenean fault, exhibits a permanent seismic activity with moderate magnitude events. From the end of the 1980s, seismic instrumentation has been deployed in the Pyrenees, making now possible the computation of improved seismicity maps. We have gathered all the arrival times published for the period 1989–1996 by the different Spanish and French institutions in charge of the seismic survey of the range, and reprocessed them in an homogeneous way, in order to obtain a coherent seismicity map over the whole range. Particular attention has been paid to the evaluation of the quality of the locations and to the focal depth determinations. The comparison with previous maps of the Pyrenean seismicity reveals significant improvements in both the quality of locations and the threshold of detection. The new seismicity map reveals that the North Pyrenean fault is active only in the western part of the range. In the central and eastern parts, the seismicity involves other tectonic units such as the Maladeta and Canigou granitic massifs, the North Pyrenean Frontal Thrust, the Tet fault and the volcanic units in Catalonia. Despite the short time interval considered, this new seismicity file may be a valuable tool for future tectonic studies.     

Evidence for a fluid flow triggered spatio-temporal migration of seismicity in the 2001 M<Subscript>w</Subscript> 7.7 Bhuj earthquake region,Gujarat, India,during 2001–2013

We studied the variations in spatial and temporal clustering of earthquake activity (during 2001–2013) in the Kachchh seismic zone, Gujarat, India, by precisely relocating 3478 events using a joint hypocentral determination (JHD) relocation technique, and high-quality arrival times of 21032 P- and 20870 S-waves. Temporal disposition of estimated station corrections of P- and S-waves suggests that the fluid flow in the causative fault zone of the 2001 Bhuj mainshock increased during 2001–2010, while it reduced during 2011–2013, due to the healing process associated with the perturbed Kachchh fault zone. We also estimated the isotropic seismic diffusivities from epicentral growth patterns, which are found to be much lower than those observed for reservoir-induced seismicity sites in the world. Finally, we analysed the spatial and temporal evolution of this earthquake sequence by solving the diffusion equation of pore-pressure relaxation caused by co- and post-seismic stress changes associated with earthquakes. The value of the isotropic diffusivity is estimated to be 100 m²/s for the Kachchh rift zone. This gives a higher permeability (after a lapse time of 14 years from the occurrence of the 2001 Bhuj mainshock) in comparison to those observed for other intraplate regions in the world. Our results suggest that the observed spatio-temporal migration of seismicity is consistent with the shallow (meteoric water circulation at 0–10 km depths) and deeper (metamorphic fluid and volatile CO₂ circulation at 10–40 km depths) fluid flows in the permeable and fractured causative fault zone of the 2001 Bhuj earthquake.     

Seismicity pattern in north Sumatra-Great Nicobar region: In search of precursor for the 26 December 2004 earthquake

We analyse the seismicity pattern including b-value in the north Sumatra-Great Nicobar region from 1976 to 2004. The analysis suggests that there were a number of significant, intermediate and short-term precursors before the magnitude 7.6 earthquake of 2 November 2002. However, they were not found to be so prominent prior to the magnitude 9.0 earthquake of 26 December 2004 though downward migration of activity and a 50-day short-term quiescence was observed before the event. The various precursors identified include post-seismic and intermediate-term quiescence of 13 and 10 years respectively, between the 1976 (magnitude 6.3) and 2002 earthquakes with two years (1990–1991) of increase in background seismicity; renewed seismicity, downward migration of seismic activity and foreshocks in 2002, just before the mainshock. Spatial variation in b-value with time indicates precursory changes in the form of high b-value zone near the epicenter preceding the mainshocks of 2004 and 2002 and temporal rise in b-value in the epicentral area before the 2002 earthquake.