Seismicity and fault aseismic deformation caused by fluid injection in decametric in-situ experiments

Seismicity induced by fluid perturbations became an important societal concern since felt earthquakes (M_w up to 6) occurred after anthropogenic activities. In order to mitigate the risks associated with undesired seismicity, as well as to be able to use the micro-seismicity as a probe for in-depth investigation of fluid-driven processes, it is of crucial importance to understand the links between seismicity, fluid pressure and flow. We have developed a series of in-situ, decameter-scale experiments of fault zone reactivation by controlled fluid injection, in order to improve the near-source geophysical and hydromechanical observations. The deployed geophysical monitoring close to the injection allows one to cover the full frequency range of the fault responses from the static deformation to the very high-frequency seismic emissions (up to 4 kHz). Here, we focus on the microseismicity (M_w ∼ –4 to –3) recorded during two fluid injection experiments in low-permeable shale and highly-fractured limestone formations. In both experiments, the spatio-temporal distribution of the seismic events, the energy balance, and the seismic velocity changes of the fractured medium show that most of the deformation does not actually emit seismic signals. The induced deformation is mainly aseismic. Based on these high-resolution multiparametric observations in the near-field, we therefore proposed a new model for injection-induced seismicity: the seismicity is not directly induced by the increasing fluid pressure, but it is rather triggered by the stress perturbations transferred from the aseismic motion caused by the injection.     

Seismicity of the Baikal rift system from regional network observations

In the paper we report the state-of-the-art of seismicity study in the Baikal rift system and the general results obtained. At present, the regional earthquake catalog for fifty years of the permanent instrumental observations consists of over 185,000 events. The spatial distribution of the epicenters, which either gather along well-delineated belts or in discrete swarms is considered in detail for different areas of the rift system. At the same time, the hypocenters are poorly constrained making it difficult to identify the fault geometry. Clustered events like aftershock sequences or earthquake swarms are typical patterns in the region; moreover, aftershocks of M ⩾ 4.7 earthquakes make up a quarter of the whole catalog. The maximum magnitude of earthquakes recorded instrumentally is M_LH7.6 for a strike-slip event in the NE part of the Baikal rift system and M_LH6.8 for a normal fault earthquake in the central part of the rift system (Lake Baikal basin). Predominant movement type is normal faulting on NE striking faults with a left lateral strike-slip component on W–E planes. In conclusion, some shortcomings of the seismic network and data processing are pointed out.     

Present activity and seismogenic potential of a low-angle normal fault system (Città di Castello,Italy): Constraints from surface geology,seismic reflection data and seismicity

We present new constraints on an active low-angle normal fault system in the Città di Castello–Sansepolcro basin (CSB) of the northern Apennines of Italy. New field data from the geological survey of the Carta Geologica d'Italia (CARG project) define the surface geometry of the normal fault system and lead to an interpretation of the CROP 03 deep-crust seismic reflection profile (Castiglion Fiorentino–Urbania segment), with particular attention paid to the geometry of the Plio-Quaternary extensional structures. Surface and sub-surface geological data are integrated with instrumental and historical seismicity in order to define the seismotectonics of the area.Low-angle east-dipping reflectors are the seismic expression of the well-known Altotiberina Fault (AF), a regional extensional detachment on which both east- and west-dipping high-angle faults, bounding the CSB, sole out. The AF breakaway zone is located ～ 10 km west of the CSB. Within the extensional allochthon, synthetic east-dipping planes prevail. Displacement along the AF is ～ 4.5 km, which agrees with the cumulative offset due to its synthetic splays. The evolution of the CSB has mainly been controlled by the east-dipping fault system, at least since Early Pleistocene time; this system is still active and responsible for the seismicity of the area. A low level of seismic activity was recorded instrumentally within the CSB, but several damaging earthquakes have occurred in historical times. The instrumental seismicity and the intensity data points of the largest historical earthquakes (5 events with maximum MCS intensity of IX to IX–X) allow us to propose two main seismogenic structures: the Monte Santa Maria Tiberina (Mmax = 5.9) and Città di Castello (Mmax up to 6.5) normal faults. Both are synthetic splays of the AF detachment, dipping to the NE at moderate (45–50°) to low (25–30°) angles and cutting the upper crust up to the surface. This study suggests that low-angle normal faults (at least with dips of 25–30°) may be seismogenic.     

Seismic behavior in central Taiwan: Response to stress evolution following the 1999 Mw 7.6 Chi-Chi earthquake

Following the 1999 M_w 7.6 Chi-Chi earthquake, a large amount of seismicity occurred in the Nantou region of central Taiwan. Among the seismic activities, eight M_w ⩾ 5.8 earthquakes took place following the Chi-Chi earthquake, whereas only four earthquakes with comparable magnitudes took place from 1900 to 1998. Since the seismicity rate during the Chi-Chi postseismic period has never returned to the background level, such seismicity activation cannot simply be attributed to modified Omori’s Law decay. In this work, we attempted to associate seismic activities with stress evolution. Based on our work, it appears that the spatial distribution of the consequent seismicity can be associated with increasing coseismic stress. On the contrary, the stress changes imparted by the afterslip; lower crust–upper mantle viscoelastic relaxation; and sequent events resulted in a stress drop in most of the study region. Understanding seismogenic mechanisms in terms of stress evolution would be beneficial to seismic hazard mitigation.     

Low b-value prior to the Indo-Myanmar subduction zone earthquakes and precursory swarm before the May 1995 M 6.3 earthquake

Some 455 events (m_b ⩾ 4.5) in the Indo-Myanmar subduction zone are compiled using the ISC/EHB/NEIC catalogues (1964–2011) for a systematic study of seismic precursors, b-value and swarm activity. Temporal variation of b-value is studied using the maximum likelihood method beside CUSUM algorithm. The b-values vary from 0.95 to 1.4 for the deeper (depth ⩾60 km) earthquakes, and from 0.85 to 1.3 for the shallower (depth <60 km) earthquakes. A sudden drop in the b-value, from 1.4 to 0.9, prior to the occurrence of larger earthquake(s) at the deeper depth is observed. It is also noted that the CUSUM gradient reversed before the occurrence of larger earthquakes. We further examined the seismicity pattern for the period 1988–1995 within a radius of 150 km around the epicentre (latitude: 24.96°N; longitude: 95.30°E) of a deeper event M 6.3 of May 6, 1995 in this subduction zone. A precursory swarm during January 1989 to July 1992 and quiescence during August 1992 to April 1995 are identified before this large earthquake. These observations are encouraging to monitor seismic precursors for the deeper events in this subduction zone.     

Bearing of plate geometry and rheology on shallow-focus mega-thrust seismicity with special reference to 26 December 2004 Sumatra event

We propose a model pertaining to the generation of 26th December 2004 off Sumatra mega-event in the backdrop of other similar type earthquakes along subduction zones around the world. Reconstructions of Benioff trajectories through the hypocenters of historical earthquakes including six mega-earthquakes indicate (i) confinement of hypocenters right within the descending lithosphere, and (ii) natural coincidence of foci of the mega-events around the zones of plate flexing. These observations are discussed in detail with special emphasis on the Sumatra margin considering the role of rheological anomaly across the cross-section of the descending lithosphere; yield strength envelope and residual stress accumulation through time. The intraplate origin of shallow mega-thrust earthquakes allowed us to advocate the ‘zone of flexing’ along the profiles of the subducting plates as nodal area for stress concentration. We propose here that at elevated confining pressure and temperature, loading of unidirectional cyclic stress on time-average bending stress enhanced the material yield strength (i.e., strain-hardening), and leads the semi-brittle portion of the lithosphere into near-brittle condition through rheological transformation. Under subsequent rise in neutral surface and increase in compressive stress field, non-coaxial deformation triggered shear failure on 26th December 2004 preferably at the rheological interface between strain-hardened near-brittle layer and deformed ductile layer within the sub-oceanic mantle.A two-stage fracture mechanism viz. a slow (～1.1 km/s) bilateral initiation in an essentially strain-hardened near-brittle domain and a follow-up very rapid progression (3.3 km/s) in the brittle, crustal domain was mainly involved in the generation of 2004 off Sumatra mega-event. Estimation shows an amount of 3.38 × 10²² to 4.50 × 10²² N m seismic moment (M_o) and 8.95–9.03 moment magnitude (M_w) for the southern part of the 1300 km extended rupture i.e. between the North Andaman to the north and the Sumatra at its south. The study necessitates the reassessment of other shallow-focus mega-thrust earthquakes along the subduction margins around the globe.     

Hierarchy of earthquakes in the Baikal rift system: implications for lithospheric stress

We investigated space-time lithospheric stress patterns of the Baikal rift system according to the hierarchy of earthquakes using mechanisms of 265 K_P ≥ 10 events recorded from 1950 to 1998 and seismic moments of 802 K_P ≥ 11 events from 1968 to 1994. The lithosphere of the region was confirmed to undergo rifting with mostly normal-slip events, while local areas of frequent strike-slip and reverse motions may record stress heterogeneity. The dominance of rifting, although being evident in the stress dynamics, is unstable, which is indicated by increase in strike-slip and reverse motions to as many as normal slip events in the latest 1980s–earliest 1990s. The lithospheric stress patterns inferred from seismic-moment data are generally consistent with those derived from the classical focal mechanism method. The suggested approach of seismic zoning according to earthquake slip geometry may provide a more reliable background for successful mitigation of seismic hazard in the region.     

P,S wave velocity model of the crust and upper most mantle of Albania region

This paper describes the one-dimensional (1D) velocity model computed by VELEST in the SEISAN seismic analysis system, inverting re-picked P-wave and S-wave arrival times recorded during 2002–2006 by the Albanian, Montenegro, Thessalonica and Macedonia seismic networks. The re-picked data yield P-wave and S-wave velocities proved to be more suitable compared to bulletin data for this detailed inversion study. Seismic phases recorded by the Albania seismic network and integrated with data from the Montenegro, Thessalonica and Macedonia networks are used to prepare the Albanian seismic bulletin. Earthquake hypocenters from the Albanian bulletins have also location errors that are negligible for civil protection purposes, large scale seismotectonic analyses and more accurate hypocentral determinations which are necessary for detailed seismotectonic and geodynamic studies.It was noted that the smoothness of the velocity variation increased with depth. A velocity of 5.5 km/s was calculated for the upper crust, 6.1 km/s was calculated for the middle crust and 6.9 km/s was computed for the lower crust. P wave velocity was 7.85 km/s at depth of 50 km and for the upper mantle it is 8.28 km/s. Using the improved velocity model, the earthquakes which occurred in Albania in the past 5 years were able to be relocated, achieving constrained hypocentral determinations for events in Albania. The interpretation of the 1 D velocity models infers interesting features of the deep structure of Albania. These results represent an important step towards more detailed seismotectonic analyses.     

Shear wave attenuation characteristics over the Central India Tectonic Zone and its surroundings

The study area lies between latitude 18–26°N and longitude 73–83°E, and mainly covers the Central India Tectonic Zone (CITZ). The frequency-dependent shear wave quality factor (Q_s) has been estimated over the CITZ and its surroundings using Double Spectral Ratio (DSR) method. We have considered 25 local earthquakes with magnitude (M_L) varies from 3.0 to 4.7 recorded at 11 stations running under national seismic network. The Fast Fourier Transformed (FFT) spectra were computed from the recorded waveform having time-window from onset of S-phase to 1.0 s and for a frequency-band of 0.1–10 Hz. Three different shear wave velocities (i.e., 3.87, 3.39 and 3.96 km/s) were obtained over the study area based on a pair of earthquakes recorded at a pair of stations. The low Q_s values of 51–96 at 1 Hz (i.e., Q_s = 51f^0.49; Q_s = 90f^0.488 and Q_s = 96f^0.53) were found in the area covering the Son–Narmada–Tapti (SONATA) lineament, CITZ, eastern part of the Satpura fold belt, Vindhyan and Gondwana basins, Godavari and Mahanadi grabens, and southern part of Gangetic plain. Intermediate Q_s values of the order of 204–277 (i.e., Q_s = 204f^0.56 and Q_s = 277f^0.55) were noted in the cartonic areas, namely, Bundelkhand, Dharwar-Bhandara and Bastar. While the higher Q_s values of 391–628 at 1 Hz (i.e., Q_s = 391f^0.49, Q_s = 409f^0.48, Q_s = 417f^0.48, Q_s = 500f^0.66, Q_s = 585f^0.65 and Q_s = 628f^0.69) were found in the eastern part of the SONATA, CITZ, and the northeastern part of the Satpura fold belt. The low Q_s values might be attributing to the more heterogeneous SONATA rift system. Low Q_s values further may presumably be associated with lower-level of seismicity and apparently account for higher tectonic stress accumulation over long duration. The long-term accumulated stress is generally released through occasional triggering of moderate magnitude earthquakes in the SONATA zone. Surrounding the SONATA region, the higher Q_s values possibly accounts for a more homogeneous subsurface structure along the SONATA zone.     

The intracratonic Caraíbas–Itacarambi earthquake of December 09, 2007 (4.9 mb), Minas Gerais State,Brazil

On December 9, 2007, a 4.9 m_b earthquake occurred in the middle of the São Francisco Craton, in a region with no known previous activity larger than 4 m_b. This event reached intensity VII MM (Modified Mercalli) causing the first fatal victim in Brazil. The activity had started in May 25, 2007 with a 3.5 magnitude event and continued for several months, motivating the deployment of a local 6-station network. A three week seismic quiescence was observed before the mainshock. Initial absolute hypocenters were calculated with best fitting velocity models and then relative locations were determined with hypoDD. The aftershock distribution indicates a 3 km long rupture for the mainshock. The fault plane solution, based on P-wave polarities and hypocentral trend, indicates a reverse faulting mechanism on a N30°Ε striking plane dipping about 40° to the SE. The rupture depth extends from about 0.3 to 1.2 km only. Despite the shallow depth of the mainshock, no surface feature could be correlated with the fault plane. Aeromagnetic data in the epicentral area show short-wavelength lineaments trending NNE–SSW to NE–SW which we interpret as faults and fractures in the craton basement beneath the surface limestone layer. We propose that the Caraíbas–Itacarambi seismicity is probably associated with reactivation of these basement fractures and faults under the present E–W compressional stress field in this region of the South American Plate.     

Contamination of seismicity catalogs by quarry blasts: An example from İstanbul and its vicinity,northwestern Turkey

Scientists have proposed two fault systems of different ages in the Sea of Marmara: the Thrace-Eski?ehir Fault Zone of Early Miocene–Early Pliocene age and the North Anatolian Fault Zone of Late Pliocene–Recent age. Different seismicity rates and extensions of these faults onto land near ?stanbul have been suggested. One of the reasons for these differences is the contamination of seismicity catalogs by seismic events from quarries operated in ?stanbul and its vicinity, including Gaziosmanpa?a (Cebeci and Kemerburgaz), Çatalca, Ömerli, Gebze, and Hereke.In this study, we investigated waveforms of 179 seismic events (1.8 < Md < 3.0) from the KOERI, NEMC digital database. We determined differences between earthquakes and quarry blasts based on time- and frequency-domain analyses of their seismograms (amplitude peak ratio, power ratio, and spectral amplitude ratio) and used these differences as discriminants. The results of this study indicate that 15% and 85% of the investigated seismic events are earthquakes and quarry blasts, respectively.     

Detailed fault structure of the Tarutung Pull-Apart Basin in Sumatra,Indonesia, derived from local earthquake data

The Tarutung Basin is located at a right step-over in the northern central segment of the dextral strike-slip Sumatran Fault System (SFS). Details of the fault structure along the Tarutung Basin are derived from the relocations of seismicity as well as from focal mechanism and structural geology. The seismicity distribution derived by a 3D inversion for hypocenter relocation is clustered according to a fault-like seismicity distribution. The seismicity is relocated with a double-difference technique (HYPODD) involving the waveform cross-correlations. We used 46,904 and 3191 arrival differences obtained from catalogue data and cross-correlation analysis, respectively. Focal mechanisms of events were analyzed by applying a grid search method (HASH code). Although there is no significant shift of the hypocenters (10.8 m in average) and centroids (167 m in average), the application of the double difference relocation sharpens the earthquake distribution. The earthquake lineation reflects the fault system, the extensional duplex fault system, and the negative flower structure within the Tarutung Basin. The focal mechanisms of events at the edge of the basin are dominantly of strike-slip type representing the dextral strike-slip Sumatran Fault System. The almost north–south striking normal fault events along extensional zones beneath the basin correlate with the maximum principal stress direction which is the direction of the Indo-Australian plate motion. The extensional zones form an en-echelon pattern indicated by the presence of strike-slip faults striking NE–SW to NW–SE events. The detailed characteristics of the fault system derived from the seismological study are also corroborated by structural geology at the surface.     

Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

Seismicity and velocity structures along the south-Alpine thrust front of the Venetian Alps (NE-Italy)

In this paper we show the seismicity and velocity structure of a segment of the Alpine retro-belt front along the continental collision margin of the Venetian Alps (NE Italy). Our goal is to gain insight on the buried structures and deep fault geometry in a “silent” area, i.e., an area with poor instrumental seismicity but high potential for future earthquakes, as indicated by historical earthquakes (1695 Me = 6.7 Asolo and 1936 Ms = 5.8 Bosco del Cansiglio). Local earthquakes recorded by a dense temporary seismic network are used to compute 3-D Vp and Vp/Vs tomographic images, yielding well resolved images of the upper crust underneath the south-Alpine front. We show the presence of two main distinct high Vp S-verging thrust units, the innermost coincides with the piedmont hill and the outermost is buried under a thick pile of sediments in the Po plain.Background seismicity and Vp/Vs anomalies, interpreted as cracked fluid-filled volumes, suggest that the NE portion of the outermost blind thrust and its oblique/lateral ramps may be a zone of high fluid pressure prone to future earthquakes.Three-dimensional focal mechanisms show compressive and transpressive solutions, in agreement with the tectonic setting, stress field maps and geodetic observations. The bulk of the microseismicity is clustered in two different areas, both in correspondence of inherited lateral ramps of the thrust system. Tomographic images highlight the influence of the paleogeographic setting in the tectonic style and seismic activity of the region.     

An overview of the earth crust under China

We assess the results of a number of deep seismic soundings performed in China over the last few decades, and study the variations in crustal structure in 18 tectonic units comprised of three platforms and 15 fold systems. Thickness data on 344 Mesozoic–Cenozoic sedimentary basins, as well as data on Moho depth are collected in order to discuss the relationship between the thickness of the basins and the average thickness of the consolidated crust in each tectonic unit. The degree of mirror-image symmetry between Mesozoic–Cenozoic sedimentary basins and the uplifting topmost parts of the mantle is herein analyzed using deep geophysical data on sedimentary basins. By applying standard methods of least-squares analysis to both datasets, we have obtained both the average thickness of the consolidated crust and the mirror-image symmetry factor for every platform and tectonic fold system, thereby allowing us to explore the correlation between the depths of the bottom of the sedimentary basins and the top of the uplifting mantle. The thickness of the consolidated crust in China is found to be between 20 and 63 km, following a pattern of gradual thickening from east to west. Expressed in terms of spatial seismicity and the concentration of seismic energy, and according to the sharing-out of earthquake hypocenters in the top 80 km of the earth, the rheology of the area does not appear to suit the widely accepted “jelly–sandwich” model for the continents, which does not seem to be entirely valid in China. The findings on the mirror-image symmetry factor show that this parameter varies mainly in the range ? 0.5 to ? 1.8 compared with each other tectonic unit. Most of the tectonic systems in China appear to be isostatically compensated. Tibet is an exception, in that the symmetry factor is positive for the Gangdise–Nyainqentanglha (1.0) and Himalayan (0.1) fold systems, implying that these tectonic zones are far from being in the isostatic equilibrium of the other regions. We have also analyzed the zoning characteristics through the logarithmic relationship ln R = ln (h/|a|), using the average thickness of the consolidated crust and the absolute value of the symmetry factor. The key finding is that regardless of geographical location, all the values of R = h/|a| for those sedimentary basins in which oil/gas reservoirs have to date been found, fall into the narrow range of values between 19.38 and 37.40. There is some appeal in a possible relationship involving the ratio of crustal thickness to symmetry factor, more so when the results obtained appear to suggest a prognostic tool for exploratory practice in relation to oil/gas reservoirs.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Potential of future seismogenesis in Hebei Province (NE China) due to stress interactions between strong earthquakes

Northeast China, a densely populated area, is affected by intense seismic activity, which includes large events that caused extensive disaster and tremendous loss of life. For contributing to the continuous efforts for seismic hazard assessment, the earthquake potential from the active faults near the cities of Zhangjiakou and Langfang in Hebei Province is examined. We estimate the effect of the coseismic stress changes of strong (M ⩾ 5.0) earthquakes on the major regional active faults, and mapped Coulomb stress change onto these target faults. More importantly our calculations reveal that positive stress changes caused by the largest events of the 1976 Tangshan sequence make the Xiadian and part of Daxing fault, thus considered the most likely sites of the next strong earthquake in the study area. The accumulated static stress changes that reached a value of up to 0.4 bar onto these faults, were subsequently incorporated in earthquake probability estimates for the next 30 years.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Seismicity distribution in the vicinity of the Chile Triple Junction,Aysén Region,southern Chile

The Aysén Region, southern Chile, is the area located at the southern end of the Nazca-South America subduction zone, to the east of the Chile Triple Junction. This region has historically presented low levels of seismicity mostly related to volcanism. Nonetheless, a seismic sequence occurred in 2007, related to the reactivation of the strike-slip Liquiñe-Ofqui Fault System (LOFS), confirmed that this region is not exempt from major seismic activity M ∼ 7. Here we present results from background local seismicity of two years (2004–2005) preceding the sequence of 2007. Event magnitudes range between 0.5 and 3.4 M_L and hypocenters occur at shallow depths, mostly within the upper 10 km of crust, in the overriding South American plate. No events were detected in the area locus of the 2007 sequence, and the Wadati–Benioff (WB) plane is not observable given the lack of subduction inter-plate seismicity in the area. A third of the seismicity is related to Hudson volcano activity, and sparse crustal events can be spatially associated with the trace of the Liquiñe-Ofqui fault, showing the largest detected magnitudes, in particular at the place where the two main branches of the LOFS meet. Other minor sources of seismicity correspond to glacial calving in the terminal zones of glaciers and mining explosions.