Anomalous Seismicity and Accelerating Moment Release Preceding the 2001 and 2002 Earthquakes in Northern Baja California, Mexico

— An algorithm recently developed by RUNDLE et al. (2002) to find regions of anomalous seismic activity associated with large earthquakes identified the location of an M _w = 5.6 earthquake near Calexico, Mexico. In this paper we analyze the regional seismicity before this event, and a nearby M _w = 5.7 event, using time-to-failure algorithms developed by BOWMAN et al. (1998) and BOWMAN and KING (2001a,b). The former finds the radius of a circular region surrounding the epicenter that optimizes the time-to-failure acceleration of seismic release. The latter optimizes acceleration based on the expected stress accumulation pattern for a dislocation source. Both methods found a period of accelerating seismicity in an optimal region, the size of which agrees with previously proposed scaling relations. This positive result suggests that the Rundle algorithm may provide a useful technique to identify regions of accelerating seismicity, which can then be analyzed using signal optimization time-to-failure techniques.     

The December 26, 2004, off the west coast of northern Sumatra,Indonesia, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>=9.0, earthquake and the critical-point-like model of earthquake preparation

Long-term seismic activity prior to the December 26, 2004, off the west coast of northern Sumatra, Indonesia, M _W=9.0 earthquake was investigated using the Harvard CMT catalogue. It is observed that before this great earthquake, there exists an accelerating moment release (AMR) process with the temporal scale of a quarter century and the spatial scale of 1 500 km. Within this spatial range, the M _W=9.0 event falls into the piece-wise power-law-like frequency-magnitude distribution. Therefore, in the perspective of the critical-point-like model of earthquake preparation, the failure to forecast/predict the approaching and/or the size of this earthquake is not due to the physically intrinsic unpredictability of earthquakes.     

Intermediate-Term Medium-Range Precursory Accelerating Seismicity Prior to the 12 May 2008, Wenchuan Earthquake

In the study of the predictability of great earthquakes in the perspective of seismicity analysis, two issues are presently controversial, and need more detailed studies based on real earthquake cases. The first issue is whether there exists pre-shock accelerating behavior of seismicity which is robust against the changing of spatio-temporal ranges for the sampling of seismic events, and the second is whether such an accelerating behavior is physically associated with an approach to the critical point. To answer these two questions, a retrospective case study was conducted on the 12 May 2008, Wenchuan earthquake, using the local earthquake catalogue in Sichuan and Yunnan Provinces, China, with cutoff magnitude M _L3.0, from 1977 to 2008. The results show that the answer to the first question appears to be ‘yes’; that is, in a finite spatial domain within the last couple of years before the event, clear accelerating seismicity could be observed. The answer to the second question cannot be obtained merely by examining seismicity data. However, detailed analysis of the accelerating behavior reveals a potential spatial correlation between the accelerating region and a known asperity, which might be an evidence for that the observed acceleration may have a geometrical or mechanical rather than statistical origin.     

Earthquakes Along the Passive Margin of Greenland: Evidence for Postglacial Rebound Control

—?An intriguing observation in Greenland is a clear spatial correlation between seismicity and deglaciated areas along passive continental margins, a piece of evidence for earthquake triggering due to postglacial rebound. Another piece of evidence for induced seismicity due to deglaciation derives from earthquake source mechanisms. Sparse, low magnitude seismicity has made it difficult to determine focal mechanisms from Greenland earthquakes. On the basis of two normal faulting events along deglaciated margins and from the spatial distribution of epicenters, earlier investigators suggested that the earthquakes of Greenland are due to postglacial rebound. This interpretation is tested here by using more recent data. Broadband waveforms of teleseismic P waves from the August 10, 1993 (m _b = 5.4) and October 14, 1998 (m _b = 5.1) earthquakes have been inverted for moment tensors and source parameters. Both mechanisms indicate normal faulting with small strike-slip components: the 1993 event, strike = 348.9°, dip = 41.0°, rake =?56.3°, focal depth = 11?km, seismic moment = 1.03?×?10²⁴ dyne-cm, and M _w = 5.3; the 1998 event, strike = 61.6°, dip = 58.0°, rake =?95.5°, focal depth = 5?km, seismic moment = 5.72?×?10²³ dyne-cm, and M _w = 5.1. These and the two prior events support the theory that the shallow part of the lithosphere beneath the deglaciated margins is under horizontal extension. The observed stress field can be explained as flexural stresses due to removal of ice loads and surface loads by glacial erosion. These local extensional stresses are further enhanced by the spreading stress of continental crust and reactivate preexisting faults. Earthquake characteristics observed from Greenland suggest that the dominant seismogenic stresses are from postglacial rebound and spreading of the continental lithosphere.     

Source Parameters and Scaling Relations for Reservoir Induced Seismicity in the Longtan Reservoir Area

We investigate the relationship between the impoundment and seismicity in the Longtan reservoir, southwestern China and find evidence that the seismicity was reservoir induced. After the reservoir impoundment, a pronounced increase in seismicity was observed in five clusters mainly concentrated in the areas where few earthquakes had occurred before the first filling. The observed induced seismicity shows a strong correlation with the filling cycles. The activity levels in the five clusters are different due to differences in the structures and permeabilities of the faults. Source parameters for 1,616 earthquakes with M _L 0.1–4.2 recorded by 24 fixed and temporary stations deployed around the reservoir were calculated after applying corrections for geometrical spreading, frequency-dependent Q, and site effects. The static stress drop and apparent stress in this area both appear to increase with increasing seismic moment over the entire magnitude range. Our results show that reservoir induced earthquakes have ten times lower average stress drop than natural tectonic earthquakes. These results may indicate that the reservoir induced seismicity can occur with a lower tectonic stress due to the high pore pressures of the underground medium, and that the effect of the water decreases the coefficient of friction.     

川东南长宁地区地震活动及b值演化特征

利用全国统一目录和流动台站目录,研究了四川盆地东南部长宁地区的地震活动特征和b值的空间分布特征.研究结果显示,长宁地区的地震活动在时间上呈现明显的分段特征,地震活动在2015年后明显增强;在空间上,长宁地区的地震活动主要集中在以28.3°N为界限的南、北2个地区,对于这2个区域的b值演化,计算结果显示出不同的分段特征....     

2009年7月24日西藏尼玛西南MS5.6地震研究

主要研究2009年7月24日西藏尼玛西南MS5.6地震的基本参数、地震序列特征、震源参数、发震构造等;利用震中附近600km范围内台站测定参数研究地震的震源机制解,与哈佛大学给出的震源机制解较一致,且与通过现场考察的发震断层走向具有一致性。研究认为本次地震发生在冈底斯山—拉萨块体内部,断裂为NNW向,主要受张应力作用产生左旋走滑正断层活动。此外还分析了震前地震学条带异常特征,结束表明,震前1年出现NW向条带非常显著,研究结论为该地区今后地震预测提供科学依据。     

Benioff Strain Release Before Earthquakes in China: Accelerating or Not?

We systematically analyzed the Benioff strain release before 65 earthquakes with M_S over 6.0 in China from 1978 to 2003 to investigate the generality of the widely discussed accelerating moment release (AMR) phenomenon before strong and intermediate-strength earthquakes. In this approach, a strong or intermediate-strength earthquake is selected as a ``target earthquake,' and retrospective analysis of seismic activity before the ``target earthquake' is performed. Simple searching area (three circular areas with different radius centered at the epicenter of the ``target earthquake') and unified temporal range (8 years) are taken in the analysis. Justification of whether AMR exists is by both visual inspection and by power-law curve fitting. It is found that more than 3/5 of the earthquakes under consideration exhibit clear pre-shock AMR property, and 1/3 of the events seem to be sensitive to the searching area. AMR behavior shows apparent focal mechanism dependence: 15 out of 17 dip-slip earthquakes with stable moment release characteristics against the changing of searching areas exhibit AMR behavior, while 16 out of 25 strike-slip earthquakes with stable moment release characteristics exhibit AMR behavior.     

Repeated Intermittent Earthquake Cycles in the San Francisco Bay Region

Forecasts of future earthquake hazard in the San Francisco Bay region (SFBR) are dependent on the distribution used for the possible magnitude of future events. Based on the limited observed data, it is not possible to statistically distinguish between many distributions with very different tail behavior. These include the modified and truncated Gutenberg–Richter distributions, and a composite distribution assembled by the Working Group on California Earthquake Probabilities. There is consequent ambiguity in the estimated probability of very large, and hence damaging, events. A related question is whether the energy released in earthquakes is a small or large proportion of the stored energy in the crust, corresponding loosely to the ideas of self-organized criticality, and intermittent criticality, respectively. However, the SFBR has experienced three observed accelerating moment release (AMR) cycles, terminating in the 1868 Hayward, 1906 San Andreas and 1989 Loma Prieta events. A simple stochastic model based on elastic rebound has been shown to be capable of producing repeated AMR cycles in large synthetic catalogs. We propose that such catalogs can provide the basis of a test of a given magnitude distribution, via comparisons between the AMR properties of the real and synthetic data. Our results show that the truncated Gutenberg–Richter distribution produces AMR behavior closest to the observed AMR behavior. The proviso is that the magnitude parameters b and m _max are such that a sequence of large events that suppresses activity for several centuries is unlikely to occur. Repeated simulation from the stochastic model using such distributions produces 30-year hazard estimates at various magnitudes, which are compared with the estimates from the 2003 Working Group on California Earthquake Probabilities.     

Characteristics of seismic activity before the MS8.0 Wenchuan earthquake

The characteristics of spatio-temporal seismicity evolution before the Wenchuan earthquake are studied. The results mainly involve in the trend abnormal features and its relation to the Wenchuan earthquake. The western Chinese mainland and its adjacent area has been in the seismically active period since 2001, while the seismic activity shows the obvious quiescence of M≥?7.0, M≥?6.0 and M?≥5.0 earthquakes in Chinese mainland. A quiescence area with M?≥7.0 has been formed in the middle of the North-South seismic zone since 1988, and the Wenchuan earthquake occurred just within this area. There are a background seismicity gap of M?≥5.0 earthquakes and a seismogenic gap of M_L?≥4.0 earthquakes in the area of Longmenshan fault zone and its vicinity prior to the Wenchuan earthquake. The seismic activity obviously strengthened and a doughnut-shape pattern of M?≥4.6 earthquakes is formed in the middle and southern part of the North-South seismic zone after the 2003 Dayao, Yunnan, earthquake. Sichuan and its vicinity in the middle of the doughnut-shape pattern show abnormal quiescence. At the same time, the seismicity of earthquake swarms is significant and shows heterogeneity in the temporal and spatial process. A swarm gap appears in the M4.6 seismically quiet area, and the Wenchuan earthquake occurred just on the margin of the gap. In addition, in the short term before the Wenchuan earthquake, the quiescence of earthquake with M_L≥?4.0 appears in Qinghai-Tibet block and a seismic belt of M_L?≥3.0 earthquakes, with NW striking and oblique with Longmenshan fault zone, is formed.     

云南11组MS≥6.5强震前h''值异常研究

利用极限时间法对云南省1970年以来Ms≥6．5 11组强地震进行h′值计算,结果显示：11组强震前出现h′值异常,表现为不同程度的大幅(平均减幅达28％)突然减小现象,即地震活动出现异常平静的现象,可成为地震预报意见的关键指标;强震前出现8次地震活动增强一平静一主震的发震模式,占73％,出现3次增强一平静一增强(短期内)一主震的发震模式,占23％;h′值量板图对强震前的地震活动平静和地震活动增强现象展现得较为清晰;通过对11组强震进行预测,预测时间约为5个月。     

Seismicity anomalies before the great earthquake of <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 in the Kunlun Pass and its significance to earthquake prediction

A great earthquake of M _S=8.1 took place in the west of Kunlun Pass on November 14, 2001. The epicenter is located at 36.2°N and 90.9°E. The analysis shows that some main precursory seismic patterns appear before the great earthquake, e.g., seismic gap, seismic band, increased activity, seismicity quiet and swarm activity. The evolution of the seismic patterns before the earthquake of M _S=8.1 exhibits a course very similar to that found for earthquake cases with M _S≥7. The difference is that anomalous seismicity before the earthquake of M _S=8.1 involves in the larger area coverage and higher seismic magnitude. This provides an evidence for recognizing precursor and forecasting of very large earthquake. Finally, we review the rough prediction of the great earthquake and discuss some problems related to the prediction of great earthquakes.     

Changes in source parameters of foreshocks and aftershocks of the 2001 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=6.0 Yajiang,Sichuan, earthquake

In this paper changes in focal mechanisms, parameters of wave spectra, and stress drops for the M _S=5.0 foreshock and M _S=6.0 mainshock in February 2001 in Yajiang County, Sichuan, and seismicity in epicentral region are studied. Comparison of focal mechanisms for the Yajiang earthquakes with distribution patterns of aftershocks, the nodal plane I, striking in the direction of NEN, of the Yajiang M=5.0 event is chosen as the faulting plane; the nodal plane II, striking in the direction of WNW, of the M=6.0 event as the faulting plane. The strikes of the two faulting planes are nearly perpendicular to each other. The level of stress drops in the epicentral region before the occurrence of the M=6.0 earthquake increases, which is consistent with increase of seismicity in the epicentral region. The rate decay of the Yajiang earthquake sequence, changes in wave spectra for foreshocks and aftershocks, and focal mechanisms are complex.     

The December 26, 2004, off the west coast of northern Sumatra, Indonesia, M W=9.0, earthquake and the critical-point-like model of earthquake preparation

Long-term seismic activity prior to the December 26, 2004, off the west coast of northern Sumatra, Indonesia, M _W=9.0 earthquake was investigated using the Harvard CMT catalogue. It is observed that before this great earthquake, there exists an accelerating moment release (AMR) process with the temporal scale of a quarter century and the spatial scale of 1 500 km. Within this spatial range, the M _W=9.0 event falls into the piece-wise power-law-like frequency-magnitude distribution. Therefore, in the perspective of the critical-point-like model of earthquake preparation, the failure to forecast/predict the approaching and/or the size of this earthquake is not due to the physically intrinsic unpredictability of earthquakes. Foundation item: Ministry of Science and Technology Project (2004CB418406). Contribution No. 05FE3010, Institute of Geophysics, China Earthquake Administration.     

Multifractal analysis of three large earthquakes in Chile: Antofagasta 1995, Valparaiso 1985, and Maule 2010

A multifractal analysis of seismicity of three large earthquakes in Chile is made: the Central Zone 1985 (M _W = 8.0), Antofagasta 1995 (M _W = 8.1), and Maule 2010 (M _W = 8.8) earthquakes. The analysis shows that the fractal dimension spectrum D _q decreases with time before an earthquake. This fact suggests that the spatial distribution of seismic events could form a cluster before a main shock.     

Multifractality in Seismicity Spatial Distributions: Significance and Possible Precursory Applications as Found for Two Cases in Different Tectonic Environments

We explore fractal properties of two observed seismicity distributions prior to the 2003 M _w 7.4 Colima, Mexico and 1992 M _w 7.3 Landers, USA earthquakes, together with several mathematical fractal distributions and two non-fractal ones, in order to estimate minimum reliable sample sizes, determine whether fractality for observed seismicity is essentially different from random uniform distributions, and explore the possibility of extracting premonitory information from fractal characteristics of seismicity before large earthquakes. Sample sizes above 800 events for whole catalogs appear to be sufficient to maintain ordered multifractality and to yield dimension estimates that vary smoothly and reliably. Fractal estimates appear to be best for whole catalogs that include aftershocks. The fractal characteristics of spatial distributions of seismicity are essentially different from those of the uniform random distribution, which is the null hypothesis of a non-fractal distribution with minimum information. The fractal dimensions and afractality measures of seismicity distributions change with time and show distinctive behaviors associated with foreshocks and main events, although these behaviors are different for each example. Results suggest the possibility of a priori identification of foreshocks to large earthquakes. A combination of fractal dimension and afractality measures over time may be helpful in large earthquake premonitory studies.     

Asperity Distribution of the 1952 Great Kamchatka Earthquake and its Relation to Future Earthquake Potential in Kamchatka

—The 1952 Kamchatka earthquake is among the largest earthquakes of this century, with an estimated magnitude of M _w = 9.0. We inverted tide gauge records from Japan, North America, the Aleutians, and Hawaii for the asperity distribution. The results show two areas of high slip. The average slip is over 3 m, giving a seismic moment estimate of 155×10²⁰Nm, or M _w = 8.8. The 20th century seismicity of the 1952 rupture zone shows a strong correlation to the asperity distribution, which suggests that the large earthquakes (M > 7) are controlled by the locations of the asperities and that future large earthquakes will also recur in the asperity regions.     

Source Characteristics of the 12 November 1996 Mw 7.7 Peru Subduction Zone Earthquake

—The 12 November 1996 M _w 7.7 Peru subduction zone earthquake occurred off the coast of southern Peru, near the intersection of the South American trench and the highest topographical point of the subducting Nazca Ridge. We model the broadband teleseismic P-waveforms from stations in the Global Seismic Network to constrain the source characteristics of this subduction zone earthquake. We have analyzed the vertical component P-waves for this earthquake to constrain the depth, source complexity, seismic moment and rupture characteristics. The seismic moment determined from the nondiffracted P-waves is 3–5 × 10²⁰ N·m, corresponding to a moment magnitude M _w of 7.6–7.7. The source time function for the 1996 Peru event has three pulses of seismic moment release with a total duration of approximately 45–50 seconds. The largest moment release occurs at approximately 35–40 seconds and is located ～90km southeast of the rupture initiation. Approximately 70% of the seismic moment was released in the third pulse.¶We find that the 1996 event reruptured part of the rupture area of the previous event in 1942. The location of the 1996 earthquake corresponds to a region along the Peru coast with the highest uplift rates of marine terraces. This suggests that the uplift may be due to repeated earthquakes such as the 1996 and 1942 events.     

Stresses and seismicity at the present stage of evolution of the Baikal rift zone lithosphere

The paper presents results of studying stresses and seismicity of the lithosphere in the Baikal rift zone at the present (instrumental) stage of its evolution. These results are well consistent with the theory of self-organization of complex dynamic systems and can be used for the classification of certain features and properties of the Baikal rift seismogenic system studied. Application of the theory of nonlinear dynamic systems to the analysis of stresses and seismicity allowed us to develop a concept of a complex spatiotemporal structure of the stress state of the lithosphere and seismicity in the region. In terms of this concept, the distribution of strong earthquakes in time is explained in terms of bifurcations (catastrophes) of stresses in the rift zone. Extrapolation of the results indicates that a stress catastrophe in lithosphere of the rift zone can take place in the next few years years, which increases the probability of occurrence of strong (M ≈ 7) earthquakes in the Baikal region. A model with bifurcation of triple equilibrium most consistent with the phase image of regional stresses is proposed as a scenario of the stress state evolution in the lithosphere of the Baikal rift zone.     

Estimation of the probability of strongest seismic disasters based on the extreme value theory

This paper presents the review of the experience in applying the approach based on the limiting distributions of the extreme value theory (the generalized Pareto distribution, GPS, and generalized extreme value distribution, GEV) for deriving the distributions of maximal magnitudes and related ground accelerations from the earthquakes on the future time intervals of a given duration. The results of analyzing the global and regional earthquake catalogs and the ground peak accelerations during the earthquakes are described. It is shown that the magnitude of the strongest possible earthquake M _max (and analogous characteristics for other types of data), which is often used in seismic risk assessment, is potentially unstable. We suggest a stable alternative for M _max in the form of quantiles Q _q (τ) of the maximal possible earthquake, which could occur during the future time interval of length τ. The quantity of the characteristic maximal event M _c, which has been introduced in our previous publications, is another helpful robust scalar parameter. All the cases of approximation of the tails of empirical distributions, which were studied in our works, turned out to be finite (bounded); however, the rightmost point of these distributions, M _max, is often poorly detectable and unstable. Therefore, the M _max parameter has a low practical value.