Observed correlations between aftershock spatial distribution and earthquake fault lengths

We observe the spatial distributions of the magnitude of aftershocks following the six earthquakes of focal depth shallower than 20 km with magnitude more than 5.0 from 1983 to 1987 in Japan. The upper limit of the aftershock magnitude is examined as a function of the distance from mainshock hypocentre. The observed spatial distributions of the upper limit are bimodal, with a tendency of the upper limit to decrease as the distance from mainshock hypocentre increases. Moreover, we observe the correlations between the aftershock spatial distribution and earthquake fault length. We focus on the largest aftershocks in each of two aftershock sequences constituting the bimodal distribution. The distances of the two largest aftershocks from the mainshock hypocentre are equal to the fault lengths of shallow earthquakes in Japan and to the maximum earthquake fault lengths.     

Aftershock activity of Bhuj earthquake of january 26th, 2001

Following a large-sized Bhuj earthquake (M _s = 7.6) of January 26th, 2001, a small aperture 4-station temporary local network was deployed, in the epicentral area, for a period of about three weeks and resulted in the recording of more than 1800 aftershocks (-0.07 ≤M _L <5.0). Preliminary locations of epicenters of 297 aftershocks (2.0 ≤M _L <5.0) have brought out a dense cluster of aftershock activity, the center of which falls 20 km NW of Bhachau. Epicentral locations of after-shocks encompass a surface area of about 50 × 40 km² that seems to indicate the surface projection of the rupture area associated with the earthquake. The distribution of aftershock activity above magnitude 3, shows that aftershocks are nonuniformly distributed and are aligned in the north, northwest and northeast directions. The epicenter of the mainshock falls on the southern edge of the delineated zone of aftershock activity and the maximum clustering of activity occurs in close proximity of the mainshock. Well-constrained focal depths of 122 aftershocks show that 89% of the aftershocks occurred at depths ranging between 6 and 25 km and only 7% and 4% aftershocks occur at depths less than 5 and more than 25 km respectively. The Gutenberg-Richter (GR) relationship, logN = 4.52 - 0.89M_L, is fitted to the aftershock data (1.0<-M _L<5.0) and theb-value of 0.89 has been estimated for the aftershock activity.     

Aftershocks of the june 20, 1978, Greece earthquake: A multimode faulting sequence

A 10-station portable seismograph network was deployed in northern Greece to study aftershocks of the magnitude (m_b) 6.4 earthquake of June 20, 1978. The main shock occurred (in a graben) about 25 km northeast of the city of Thessaloniki and caused an east-west zone of surface rupturing 14 km long that splayed to 7 km wide at the west end. The hypocenters for 116 aftershocks in the magnitude range from 2.5 to 4.5 were determined. The epicenters for these events cover an area 30 km (east-west) by 18 km (north-south), and focal depths ranges from 4 to 12 km. Most of the aftershocks in the east half of the aftershock zone are north of the surface rupture and north of the graben. Those in the west half are located within the boundaries of the graben. Composite focalmechanism solutions for selected aftershocks indicate reactivation of geologically mapped normal faults in the area. Also, strike-slip and dip-slip faults that splay off the western end of the zone of surface ruptures may have been activated.The epicenters for four large (M 4.8) foreshocks and the main shock were relocated using the method of joint epicenter determination. Collectively, those five epicenters form an arcuate pattern convex southward, that is north of and 5 km distant from the surface rupturing. The 5-km separation, along with a focal depth of 8 km (average aftershock depth) or 16 km (NEIS main-shock depth), implies that the fault plane dips northward 58° or 73°, respectively. A preferred nodal-plane dip of 36° was determined by B.C. Papazachos and his colleagues in 1979 from a focal-mechanism solution for the main shock. If this dip is valid for the causal fault and that fault projects to the zone of surface rupturing, a decrease of dip with depth is required.     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

Improved earthquake locations in the Koyna-Warna seismic zone

In estimating the likelihood of an earthquake hazard for a seismically active region, information on the geometry of the potential source is important in quantifying the seismic hazard. The damage from an earthquake varies spatially and is governed by the fault geometry and lithology. As earthquake damage is amplified by guided seismic waves along fault zones, it is important to delineate the disposition of the fault zones by precisely determined hypocentral parameters. We used the double difference (DD) algorithm to relocate earthquakes in the Koyna-Warna seismic zone (KWSZ) region, with the P- and S-wave catalog data from relative arrival time pairs constituting the input. A significant improvement in the hypocentral estimates was achieved, with the epicentral errors <30 m and focal depth errors <75 m i.e. errors have been significantly reduced by an order of magnitude from the parameters determined by HYPO71. The earthquake activity defines three different fault segments. The seismogenic volume is shallower in the south by 3 km, with seismicity in the north extending to a depth of 11 km while in the south the deepest seismicity observed is at a depth of 8 km. By resolving the structure of seismicity in greater detail, we address the salient issues related to the seismotectonics of this region.     

Seismogenic Tectonics and Dynamics of the 2011 Ms5.9 Yingjiang Earthquake in Yunnan,China

In the southern South–North Seismic Zone, China, seismic activity in the Yingjiang area of western Yunnan increased from December 2010, and eventually a destructive earthquake of Ms5.9 occurred near Yingjiang town on 10 March 2011. The focal mechanism and hypocenter location of the mainshock suggest that the Dayingjiang Fault was the site of the mainshock rupture. However, most of foreshocks and all aftershocks recorded by a portable seismic array located close to the mainshock occurred along the N–S-striking Sudian Fault, indicating that this fault had an important influence on these shocks. Coulomb stress calculations show that three strong(magnitude ≥5.0) earthquakes that occurred in the study region in 2008 increased the coulomb stress along the plane parallel to the Dayingjiang Fault. This supports the Dayingjiang Fault, and not the Sudian Fault, as the seismogenic fault of the 2011 Ms5.9 Yingjiang earthquake. The strong earthquakes in 2008 also increased the Coulomb stress at depths of ≤5 km along the entire Sudian Fault, and by doing so increased the shallow seismic activity along the fault. This explains why the foreshocks and aftershocks of the 2011 Yingjiang earthquake were located mostly on the Sudian Fault where it cuts the shallow crust. The earthquakes at the intersection of the Sudian and Dayingjiang faults are distributed mainly along a belt that dips to the southeast at ~40°, suggesting that the Dayingjiang Fault in the mainshock area also dips to the southeast at ~40°.     

The MW = 6.3 2003 Lefkada earthquake (Greece) and induced stress transfer changes

A large earthquake of magnitude M_W = 6.3 occurred on 14 August 2003 NW of the Lefkada Island, which is situated at the Ionian Sea (western Greece). The source parameters of this event are determined using body-wave modeling. The focal depth was found equal to 9 km, the constrained focal mechanism revealed dextral strike–slip motion (φ = 15°, Δ = 80° and λ = 170°), the duration of the source time function was 8 s and the seismic moment 2.9 × 10²⁵ dyn cm. The earthquake occurred close to the northern end of the Kefallinia transform fault, where the 1994 moderate event and its aftershock sequence were also located. The epicentral distribution of the 2003 aftershock sequence revealed the existence of two clusters. The first one is located close to the epicentral area of the mainshock, while the second southern, close to the northwestern coast of the Kefallinia Island. A gap of seismicity is observed between the two clusters. The length of the activated zone is approximately 60 km. The analysis of data revealed that the northern cluster is directly related to the mainshock, while the southern one was triggered by stress transfer caused by the main event.     

Tomographic models and seismotectonics of the Reggio Emilia region, Italy

The aim of this study is to define the Vp and Vp/Vs structure of the fault zone ruptured by the M_L 5.1 earthquake of October 15, 1996 which occurred near Reggio Emilia (central-northern Italy). A 1-month-long seismic sequence followed the mainshock and occurred in a small region along the outer border of the northern Apenninic belt, at depth ranging between 10 and 17 km. P- and S-wave arrival times from 304 aftershocks recorded by two local dense seismic arrays installed in the epicentral region have been inverted to obtain one- and three-dimensional velocity models by using state of the art local earthquake tomographic techniques. Velocity models and aftershock relocation help us to infer the seismotectonic of the region. Earthquakes originated along a NW-dipping backthrust of a NE-trending main thrust, composing the western part of the broad Ferrara Arc. A main high Vp and high Vp/Vs region delineates a pop-up structure in the center of the area. The high Vp/Vs within the pop-up structure supports the presence of a zone with increased pore pressure. The hypocentral depth of both mainshock and aftershocks is greater than those usually found for the main seismogenic regions of the Apenninic belt. P-wave velocity values in the seismogenic area, obtained by tomography, are compatible with rocks of the Mesozoic cover and suggest that seismicity occurred within the Mesozoic units stack at present by compressional tectonics.     

The 20 March 1992 South Aegean,Greece, earthquake (Ms= 5.3): possible anomalous effects

The earthquake (M_s= 5.3) of 20 March 1992 and its aftershocks, which occurred near the volcanic island complex of Milos, South Aegean, Greece, are studied on the basis of filed observations and instrumental data. The mainshock caused some building damage, the maximum intensity of VI+ (MM) being assigned to Triovasalos, Milos. Ground cracks, liquefaction in soil, landslides and rockfalls were observed in Milos. Liquefaction took place at an apparently anomalously long epicentral distance (D= 12 km) and is associated with unusually small earthquake magnitude. Abnormal animal behaviour was reported no longer than twelve hours before the mainshock. The b-value (= 1.02) of the G–R relation for the aftershock sequence, the exponentially decreasing number of aftershocks with time, and the difference (= 0.5) in magnitude between the mainshock and its largest aftershock imply that the origin of these earthquakes is tectonic and not associated with the volcanic field of Milos.     

The 1988 Tennant Creek,northern territory,earthquakes: A synthesis

Three large earthquakes with surface‐wave magnitudes 6.3–6.7 on 22 January 1988 were associated with 32 km of surface faulting on two main scarps 30 km southwest of Tennant Creek in the Northern Territory. These events provide an excellent opportunity to study the mechanics of midplate earthquakes because of the abundance of geological and geophysical data in the area, the proximity of the Warramunga seismic array and the ease of access to the fault zone. The 1988 earthquakes were located in the North Australian Craton in an area that had no history of moderate or large earthquakes before 1986. Additionally, no smaller earthquakes from the fault zone were identified at the Warramunga array, which is situated only 30 km from the nearest scarp, between the 1965 installation of the array and 1986. The main shocks were preceded by a swarm of moderatesized (magnitude 4–5) earthquakes in January 1987 and many smaller aftershocks throughout 1987. Careful relocation of all teleseismically recorded earthquakes from the fault zone shows that the 1987 activity was concentrated in an area only 6 km across in the gap between the two main fault scarps. The main shocks also nucleated in the centre of the fault zone near the 1987 activity. Field observations of scarp morphology indicate that the scarp is divided into three segments, each showing primarily reverse faulting. However, whereas the western and eastern segments show movement of the southern block over the northern, the central scarp segment shows the opposite, with the northern block thrust over the southern block.

Analysis of the first arrival times at Warramunga suggests that the three main shocks were associated with the western, central and eastern scarp segments, respectively. The locations of aftershocks determined using data from temporary seismograph arrays in the epicentral area define three inclined zones of activity that are interpreted as fault planes. In the western and eastern portions of the aftershock zone, these concentrations of activity dip to the south at 45° and 35°, respectively, but in the central section the aftershock zone dips to the north at 55°. Focal mechanisms derived from modelling broadband teleseismic data show thrust and oblique thrust faulting for the three main shocks. The first event ruptured unilaterally up and to the northwest on the westernmost fault segment, while the third main shock ruptured horizontally to the southeast. Modelling of repeat levelling data from the epicentral area requires at least three distinct fault planes, with the eastern and western planes dipping to the south and the central plane dipping to the north. The combination of scarp morphology, aftershock distribution and elevation data makes a strong case for rupture of fault planes in conjugate orientation during the 22 January 1988 Tennant Creek earthquakes. More than 20000 aftershocks have been recorded at Warramunga and activity continues to the present‐day with occasional shocks felt in the town of Tennant Creek and some recent off‐fault aftershocks located directly under the Warramunga seismic array. Stratigraphic relationships exposed in trenches excavated across the scarps suggest that during the Quaternary, a large earthquake ruptured the surface along one segment of the 1988 scarps.     

Aftershock Propagation Characteristics During the First Three Hours Following the 26 December 2004 Sumatra-Andaman Earthquake

Within three hours of the mainshock rupture of the 26 December 2004 Sumatra-Andaman earthquake, 45 aftershocks occurred that are distributed all along the mega-thrust fault plane and also along the West Andaman fault. Seven of these aftershocks struck sequentially and unilaterally from the mainshock in the south towards north within 2h 9m 50.76s indicating an overall rate of aftershock propagation to the tune of 167 meters/sec. Seismic moment calculated from fault parameters gives a value of 1.2 × 10³⁰ dyne cm. Three separate fault segments are identified from distribution of aftershocks with propagation rates 330, 250 and 85 meters/sec in the southern, central and northern segments. These 7 unilaterally propagating shocks along the mega-thrust are probably not aftershocks of the mainshock rather these are sequentially triggered shocks each rupturing a small segment of the fault. Location of the mainshock and several aftershocks are guided by several lithospheric hinge faults identified previously.     

Recent activity and seismotectonics of the Eastern Pyrenees

Results from a recent earthquake in the Eastern Pyrenees are presented and the seismotectonics of the region is analyzed from the presently available data. On 26 September 1984 an earthquake (M_L = 4.4) took place in the area of the historical destructive earthquake of 1428. Several portable stations installed in the epicentral area to record aftershocks permitted of defining a precise location at 42°19.2′N, 2°10.2′E and 5 km depth. A maximum felt intensity of V (MSK) is obtained from macroseismic data. The epicentral location lies within a block bounded by E-W-trending structures and the focal solution shows right-lateral shearing with a NW-SE pressure axis.The seismicity in the Eastern Pyrenees shows a complex pattern which can be associated with both E-W fractures and NE-SW fault systems. Focal solutions of another two recent earthquakes of M_L ~ 4, with differences in horizontal pressure axis, are also discussed.     

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.     

断层旋性与地震危险性

文中以“平行同旋走滑断层减震”的观点论证了兰州、北京、昆明这些位于强震活动区的大城市今后百年内不会发生 6 .5级以上地震。以“平行异旋走滑断层加震”的观点解释了西南棱形块体北边界和南边界在发生大震方面相互促进的现象。对于由构造分段求震级来说 ,在遇到不同的横交断层作为分段点时 ,还需考虑将来发震时始破裂点的位置以及发震断层的旋性 ,不然就会造成对震级估计不足 ,继而成为抗震建设的潜在不安全因素。在主震后为了预报余震的强度 ,可应用物理学中的科里奥利力理论 ,应用时必须知道断层类型和旋性。对于走滑断层来说 ,左旋余震弱 ,右旋余震强。例如 1997年藏北玛尼 7.5级地震 ,余震仅为 5 .3级 ,震级偏小 ,因主震为左旋的缘故 ;1976年唐山 7.8级地震 ,余震可达 7.1级 ,因主震是右旋。对于逆断层来说 ,上盘错动方向在当地子午面左侧者余震强度大 ,在右侧者余震强度弱。据此讨论了 1999年台湾南投 7.6级大震余震强度达 7.1级是因为主震为逆断层 ,上盘向西错动。     

Aftershock Propagation Characteristics During the First Three Hours Following the 26 December 2004 Sumatra-Andaman Earthquake

Within three hours of the mainshock rupture of the 26 December 2004 Sumatra-Andaman earthquake, 45 aftershocks occurred that are distributed all along the mega-thrust fault plane and also along the West Andaman fault. Seven of these aftershocks struck sequentially and unilaterally from the mainshock in the south towards north within 2h 9m 50.76s indicating an overall rate of aftershock propagation to the tune of 167 meters/sec. Seismic moment calculated from fault parameters gives a value of 1.2 × 10³⁰ dyne cm. Three separate fault segments are identified from distribution of aftershocks with propagation rates 330, 250 and 85 meters/sec in the southern, central and northern segments. These 7 unilaterally propagating shocks along the mega-thrust are probably not aftershocks of the mainshock rather these are sequentially triggered shocks each rupturing a small segment of the fault. Location of the mainshock and several aftershocks are guided by several lithospheric hinge faults identified previously.     

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.     

The 2002 Molise earthquake sequence: What can we learn about the tectonics of southern Italy?

On October 31, 2002 a M_L = 5.4 earthquake occurred in southern Italy, at the margin between the Apenninic thrust belt (to the west) and the Adriatic plate (to the east). In this area, neither historical event nor seismogenic fault is reported in the literature. In spite of its moderate magnitude, the earthquake caused severe damage in cities close to the epicenter and 27 people, out of a total of 29 casualties, were killed by the collapse of a primary school in S. Giuliano di Puglia. By inverting broadband regional waveforms, we computed moment tensor solutions for 15 events, as small as M_L = 3.5 (M_w = 3.7). The obtained focal mechanisms show pure strike-slip geometry, mainly with focal planes oriented to NS (sinistral) and EW (dextral). In several solutions focal planes are rotated counterclockwise, in particular for later events, occurring west of the mainshock. From the relocated aftershock distribution, we found that the mainshock ruptured along an EW plane, and the fault mechanisms of some aftershocks were not consistent with the mainshock fault plane. The observed stress field, resulting from the stress tensor inversion, shows a maximum principal stress axis with an east–west trend (N83°W), whereas the minimum stress direction is almost N–S. Considering both the aftershock distribution and moment tensor solutions, it appears that several pre-existing faults were activated rather than a single planar fault associated with the mainshock. The finite fault analysis shows a very simple slip distribution with a slow rupture velocity of 1.1 km/s, that could explain the occurrence of a second mainshock about 30 h after. Finally, we attempt to interpret how the Molise sequence is related to the normal faulting system to the west (along the Apennines) and the dextral strike-slip Mattinata fault to the east.     

Fine structure of seismotectonics in western Shillong massif,north east India

The epicentral tract of the great Assam earthquake of 1897 of magnitude 8·7 was monitored for about 6 months using an array of portable seismographs. The observed seismicity pattern shows several diversely-oriented linear trends, some of which either encompass or parallel known geological faults. A vast majority of the recorded micro-earthquakes had estimated focal depths between 8–14 km. The maximum estimated depth was 45 km. On the basis of a seismic velocity model for the region reported recently and these depth estimates we suggest that the rupture zone of the great 1897 earthquake had a depth of 11–12 km under the western half of the Shillong massif. Four composite fault plane solutions define the nature of dislocation in three of the seismic zones. Three of them show oblique thrusting while one shows pure dip slip reverse faulting. The fault plane solutions fit into a regional pattern of a belt of earthquakes extending in NW-SE direction across the north eastern corner of the Bengal basin. The maximum principle stress axis is approximately NS for all the solutions in conformity with the inferred direction of the Indian-EuroAsian plate convergence in the eastern Himalaya.     

Source study and tectonic implications of the 1995 Ventimiglia (border of Italy and France) earthquake (ML=4.7)

We analyse the source process and the aftershock distribution of the April 21, 1995, Ventimiglia, M_L=4.7 earthquake using the records of permanent high dynamic broad-band seismic stations and a temporary network deployed on land and at sea few hours after the earthquake. This event occurred on the western Mediterranean coast, near the border between Italy and France, at a depth of 9 km, at a point where Alpine tectonic units and Late Oligocene extensional structure overlap and are currently undergoing compressional stress. The focal solutions of the mainshock and three aftershocks depict a dominant reverse faulting with an important strike-slip component, which underlines two nodal planes: a NW–SE-dipping north fault and a NE–SW-dipping south fault. We operate a careful re-location of the aftershocks using a master-event technique and data from the temporal network and obtain a predominant NW–SE alignment. Then, we analyse the rupture process using an empirical Green function approach. We find that the mainshock broke a 0.5 to 1 km fault length and that the rupture propagated during 0.1–0.2 s probably in a SE direction. Those two arguments, together with the recent fault trace that exists close to the epicentre, leads us to propose that this event expresses the reactivation of an old transverse NW–SE structure with a dextral movement. This study thus emphasizes the role of inherited, deep-rooted, transcurrent features in the tectonic reactivation of this passive margin. It also underlines the importance of combining short-period and broad-band seismology to better resolve and understand regional tectonic processes in areas of moderate seismic activity and complex geology.