Maximum magnitudes in aftershock sequences in Taiwan

In this work, Båth’s Law, the b-value in Gutenberg–Richter Law (G–R Law) in the form of the 1/β relationship, and both the a- and b-values in the G–R Law were introduced in order to estimate maximum aftershock magnitudes of earthquake sequences in the Taiwan region. The averaged difference of magnitude between the mainshock and the maximum aftershock is 1.20, and is consistent with Båth’s Law, however, with a large uncertainty. The large uncertainty implies that the difference may result from a variable controlled by other factors, such as the aftershocks number of an earthquake sequence and magnitude threshold for mainshock. With 1/β, since 86% of the earthquake sequences with a M ⩾ 6.0 mainshock follow this relationship, the upper bound of the maximum magnitude can be estimated for an earthquake sequence with a large mainshock. The a- and b-values in the G–R Law was also considered by evaluating maximum aftershock magnitudes. As there are low residuals between the model and the observations, the results suggest that the G–R Law is a good index for maximum aftershock magnitude determinations. In order to evaluate the temporal decays of maximum aftershock magnitudes, modified Omori’s Law was introduced. Using the approaches mentioned above, the maximum magnitudes and the temporal evolution of an earthquake sequence could be modeled. Among them, the model of the G–R Law has the best fit with observations for most of earthquake sequences. It shows its feasibility. The results of this work may benefit seismic hazards mitigation in the form of rapid re-evaluations for short-term seismic hazards immediately following devastating earthquakes.     

Spatial variation of aftershock activity across the rupture zone of the 17 August 1999 Izmit earthquake, Turkey

We examined the spatial variation in the aftershock activity from the 17 August 1999 Izmit, Turkey earthquake. We found that this aftershock sequence is non-uniform both in space and time, aspects that need to be taken into account in any further statistical analysis. Other aspects of this aftershock sequence are similar to other aftershock sequences, namely low b-values and a high degree of spatial variation. We have detected three zones of relatively high b-values, two of which coincide with asperities revealed by previous slip inversion studies. The third zone with an anomalous b-value is located beyond the fault rupture and indicates a weakened fractured zone in the Yalova-Tuzla area. This b-value analysis provided no evidence for any significant difference that may exist between the two sides of the mainshock fault plane.     

Properties of the aftershocks sequences of the 2000 and 2002 earthquakes in Aksehir-Afyon graben,west-central Anatolia,Turkey

The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤  b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days.     

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.     

Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

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.     

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.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

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.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

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.     

Repeating aftershocks of the great 2004 Sumatra and 2005 Nias earthquakes

We investigate repeating aftershocks associated with the great 2004 Sumatra–Andaman (Mw 9.2) and 2005 Nias–Simeulue (Mw 8.6) earthquakes by cross-correlating waveforms recorded by the regional seismographic station PSI and teleseismic stations. We identify 10 and 18 correlated aftershock sequences associated with the great 2004 Sumatra and 2005 Nias earthquakes, respectively. The majority of the correlated aftershock sequences are located near the down-dip end of a large afterslip patch. We determine the precise relative locations of event pairs among these sequences and estimate the source rupture areas. The correlated event pairs identified are appropriately referred to as repeating aftershocks, in that the source rupture areas are comparable and significantly overlap within a sequence. We use the repeating aftershocks to estimate afterslip based on the slip-seismic moment scaling relationship and to infer the temporal decay rate of the recurrence interval. The estimated afterslip resembles that measured from the near-field geodetic data to the first order. The decay rate of repeating aftershocks as a function of lapse time t follows a power-law decay 1/t^p with the exponent p in the range 0.8–1.1. Both types of observations indicate that repeating aftershocks are governed by post-seismic afterslip.     

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.     

Finite-fault slip models for the 15 April 1979 (Mw 7.1) Montenegro earthquake and its strongest aftershock of 24 May 1979 (Mw 6.2)

We revisit the April 1979 Montenegro earthquake sequence to invert for finite-fault slip models for the mainshock of 15 April 1979 (M_w 7.1) and of the strongest aftershock of 24 May 1979 (M_w 6.2) using P, SH and SV waveforms, retrieved from IRIS data center. We also used body waveform modelling inversion to confirm the focal mechanism of the mainshock as a pure thrust mechanism and rule out the existence of considerable strike slip component in the motion. The mainshock occurred along a shallow (depth 7 km), low angle (14°) thrust fault, parallel to the coastline and dipping to the NE. Our preferred slip distribution model for the mainshock indicates that rupture initiated from SE and propagated towards NW, with a speed of 2.0 km/s. Moment was released in a main slip patch, confined in an area of L  50 km × W  23 km. The maximum slip ( 2.7 m) occurred  30 km to the NW of the hypocenter (location of rupture initiation). The average slip is 49 cm and the total moment release over the fault is 4.38e19 Nm. The slip model adequately fits the distribution of the M_w ≥ 4.3 aftershocks, as most of them are located in the regions of the fault plane that did not slip during the mainshock. The 24 May 1979 (M_w 6.2) strongest aftershock occurred  40 km NW of the mainshock. Our preferred slip model for this event showed a characteristic two-lobe pattern, where each lobe is  7.5 × 7.5 km². Rupture initiated in the NW lobe, where the slip obtained its maximum value of 45 cm, very close to the hypocenter, and propagated towards the south-eastern lobe where it reached another maximum value — for this lobe — of 30 cm, approximately 10 km away from the hypocenter. To indirectly validate our slip models we produced synthetic PGV maps (Shake maps) and we compared our predictions with observations of ground shaking from strong motion records. All comparisons were made for rock soil conditions and in general our slip models adequately fit the observations especially at the closest stations where the shaking was considerably stronger. Through the search of the parameter space for our inversions we obtained an optimum location for the mainshock at 42.04°N and 19.21° E and we also observed that better fit to the observations was obtained when the fault was modeled as a blind thrust fault.     

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.     

Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Relative stress variations as determined by b-values from earthquakes in circum-pacific subduction zones

Assuming a relation of “b” to stress state, the possibility of globe-wide stress variation and transmission was investigated. The NOAA earthquake data file served to determine the temporal change in “b” of log N = a − bM from 1963 to 1975.Periods of six to eight years are observed in the b-values (stress pattern) for most circum-Pacific areas (South America, Tonga, Kermadec, New Hebrides, Kamchatka and Eastern Aleutians).In the Kurils, fore- and aftershock sequences of large earthquakes seem to mask any characteristic global pattern that might exist. These sequences exhibit low b-values (high stress) through the time of foreshocks and early stages of aftershocks, followed by rapid increase in b-values (decrease in stress).Use of a worldwide earthquake data file clearly yields less resolution of the temporal “b” variation than the use of local network studies published by other authors.Incidental to the study, 1124 earthquakes of the NOAA data file yield the M_s − m_b relations: M_s = 1.16m_b − 0.835 for 4.5 m_b 6 and: log₁₀M_s = 0.1432m_b − 0.0629 form_b > 6 with correlative coefficients of 0.994 and 0.992 respectively.     

Comparing different models of aftershock rate decay: The role of catalog incompleteness in the first times after main shock

We evaluated the efficiency of various models in describing the time decay of aftershock rate of 47 simple sequences occurred in California (37) from 1933 to 2004 and in Italy (10) from 1976 to 2004. We compared the models by the corrected Akaike Information Criterion (AIC_c) and the Bayesian Information Criterion (BIC), both based on the log-likelihood function and also including a penalty term that takes into account the number of independent observations and of free parameters of each model. These criteria follow two different approaches (probabilistic and Bayesian respectively) well covering the wide spectra of current views on model comparison. To evaluate the role of catalog incompleteness in the first times after the main shock, we compared the performance of different models by varying the starting time T_s and the minimum magnitude threshold M_min for each sequence. We found that Omori-type models including parameter c are preferable to those not including it, only for short T_s and low M_min while the latters generally perform better than the formers for T_s longer than a few hours and M_min larger than the main shock magnitude M_m minus 3 units. For T_s > 1 day or M_min > M_m − 2.5, only about 15% of the sequences still give a preference to models including c. This clearly indicates that a value of parameter c different from zero does not represent a general property of aftershock sequences in California and Italy but it is very likely induced in most cases by catalog incompleteness in the first times after the main shock. We also considered other models of aftershock decay proposed in the literature: the Stretched Exponential Law in two forms (including and not including a time shift) and the band Limited Power Law (LPL). We found that such models perform worse than the Modified Omori Model (MOM) and other Omori-type models for the large majority of sequences, although for LPL, the relatively short duration of the analyzed sequences (one year) might also contribute to its poor performance. Our analysis demonstrates that the MOM with c kept fixed to 0 represent the better choice for the modeling (and the forecasting) of simple sequence behavior in California and Italy.     

Static and dynamic fault parameters of the Saitama earthquake of July 1, 1968

The source mechanism of the Saitama earthquake (36.07°N,139.40°E, M_s = 5.4) of July 1, 1968, is studied on the basis of P-wave first motion, aftershock, long-period surface-wave data and low-magnification long-period seismograms recorded in the nearfield. A precise location of the aftershocks is made using P and S—P time data obtained by a micro-earthquake observatory network. The synthetic near-field seismograms based on the Haskell model are directly compared with the observed near-field seismograms for wave form and amplitude to determine the dynamic fault parameters. The results obtained are as follows: source geometry, reverse dip slip with considerable right-lateral strike-slip component; dip direction, N6°E; dip angle 30°; fault dimension, 10 × 6 km²; rupture velocity, 3.4 km/sec in the direction S30°E; average dislocation, 92 cm; average dislocation velocity, 92 cm/sec; seismic moment, 1.9 · 10²⁵ dyn-cm; stress drop, 100 bar. The effective stress is about the same as the stress drop. For major earthquakes in the Japanese Islands, the dislocation velocity, .D, is found to be proportional to the stress drop, σ. This relation can be expressed by .D - (β/μ)σ, where β is the shear velocity and μ is the rigidity. This result has an importance in engineering seismology because the stress drop scales the seismic motion in the vicinity of an earthquake fault.