利用地震资料和天文周期分析的方法开展中期地震预报

地震是一种复杂的自然现象,其成因及发展的物理过程,目前尚在探索之中。地震预报在世界上还是一个尚未解决的难题。但是,地震预报并不是不可能的。地震既然是一种自然现象,必然有其特定的规律。     

卫星热红外地震短临预测方法研究

文章扼要介绍卫星热红外亮温异常的基本原理,分析建立地震短临预测遥感信息模型,通过11 a试验性预报实践和几个地震的案例,如南海东沙地震,台湾集集大地震,震例还有内蒙包头地震,美国华盛顿州奥林匹亚地震及印尼苏门答腊巨震等。从作为震兆的亮温异常提升到热应力场,对预测地震震中更加可靠和更有说服力,说明此项技术具有实际应用前景。     

Earthquake prediction: a political problem?

Earthquake prediction is an area of research of great scientific and public fascination. The reason for this is not only that earthquakes can cause extremely large numbers of fatalities in a short time, but also because earthquakes can have a large social and economic impact on society. Earthquake prediction in the sense of making deterministic predictions about the place, time, and magnitude of earthquakes may very well be fundamentally impossible. However, based on a variety of data, earth scientists can make statements about the probability that earthquakes with a certain size will occur in a certain region over a specified time period. In this context one speaks of “earthquake forecasting.” A number of methods to achieve this are presented. However, it is not obvious how society should respond to these forecasts. It is shown that there is a fundamental dilemma for decision makers that statements of scientists concerning earthquake occurrence either contain very specific information but are very uncertain, or contain very general information but are very certain. Earthquake hazard can to a large extent be reduced by formulating and enforcing appropriate building codes. However, given the fact that the majority of the population that is threatened by earthquakes is living in the third world, it is clear that this cannot easily be realized. For these reasons, earthquake prediction is not only a scientific problem: it also has a complex political dimension.     

Statistical analysis and long-term prediction of seismicity for linear zones

The model of the Poisson point process is too vague for earthquake locations in space and time: earthquakes tend to cluster in middle distances and to repulse in large ones. The Poisson point model with variable density makes it possible to describe the tendency for clustering but does not reveal the periodicity of clusters. The author proposes the point-process model where locations of points are determined not by densities of point distribution, but by densities of interpoint differences distribution. In the model, a latent periodicity is revealed and used for prediction of a point process. In 1983, the point-process model prediction was made for the Kuril Islands for 1983–1987 and two signs of danger in time and location were determined. Then they were confirmed by strong earth-quakes. In 1989, a similar prediction was made for North Armenia. The Spitak earthquake in 1988 is clearly seen from the data of previous earthquakes.     

Variation of singularity of earthquake-size distribution with respect to tectonic regime

Frequency-size relation of earthquakes in a region can be approximated by the Gutenberg-Richter law(GR). This power-law model involves two parameters: a-value measuring seismic activity or earthquake productivity, and b-value describing the relation between frequencies of small and large earthquakes.The spatial and temporal variations of these two parameters, especially the b-value, have been substantially investigated. For example, it has been shown that b-value depends inversely on differential stress. The b-value has also been utilized as earthquake precursor in large earthquake prediction.However, the physical meaning and properties of b-value including its value range still remain as an open fundamental question. We explore the property of b-value from frequency-size GR model in a new form which relates average energy release and probability of large earthquakes. Based on this new form of GR relation the b-value can be related to the singularity index(1-2/3 b) of fractal energy-probability power-law model. This model as applied to the global database of earthquakes with size M ≥ 5 from 1964 to 2015 indicates a systematic increase of singularity from earthquakes occurring on mid-ocean ridges, to those in subduction zones and in collision zones.     

大地震临震预测的研究进展

短临地震预报尤其是临震预报 ,是当今国内外公认的世界性科学难点。文章综述了中国科学家从多学科交叉和内外因耦合的科学途径 ,已在大地震临震信号方面获得了次声波异常、地应力突跳、虎皮鹦鹉跳跃异常、地电和地磁异常、卫星红外异常、重力高频脉冲等临震预测的手段。笔者通过典型分析和普查对比各种地球物理因子和各种天文因子可能影响的整体研究 ,找到了大地震临震的主要触发因子———以月亮为主的非经典引潮力共振的异常叠加。在此基础上 ,采取内外因耦合的途径和方法 ,在 6a多来联合进行的临震预测内部试验表明 ,地震三要素均基本正确的成功率已达40 %。最后 ,还对临震预报问题提出了一些看法。     

用线性预测理论研究地震中长期预报

引言长期以来,人们十分注重对地震的时间、空间序列进行研究,以探索地震发生的规律和进行预报。过去相当多的工作讨论了地震的周期性,地震在空间分布上的填空性及迁移规律问题。引入现代统计方法研究地震的时空序列,就有可能在更一般的条件下对原因和过程都尚未清楚的地震事件进行统计分析,概括出某些数学模型进行外推预测。     

Statistical evaluation of the effect of earthquake with other related factors on landslide susceptibility: using the watershed area of Shihmen reservoir in Taiwan as a case study

Landslide susceptibility evaluation is one of the most important issues in watershed management. After an earthquake, the landslide susceptibility decreases functionally with increases in the distance from the epicenter. Under the same rainfall intensity, landslides are more likely to occur in an area where earthquakes occur more frequently. However, the questions of how much an earthquake should be weighted and how to evaluate the effects of an earthquake still need to be studied. To understand how earthquakes affect rainfall-triggered landslides, the horizontal peak ground acceleration (PGA) data from the Central Weather Bureau Seismic Network is used as the earthquake factor and combined with other factors to determine the weight of earthquakes in landslide susceptibility using logistic regression. The results indicate that the ability of landslide prediction is better when considering the earthquake factor. This study also proved that although there are no co-seismic landslides (after earthquakes) in the study area, the earthquake factor is still required to increase the model accuracy. PGA has been described as a usable factor. In areas with frequent earthquakes and high geological activity, when using historical data to evaluate landslide susceptibility, the earthquake factor should be taken into consideration to prevent errors.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Process Analysis of In-situ Strain during the Ms8.0 Wenchuan Earthquake—Data from the Stress Monitoring Station at Shandan

Abstract: There were huge life and property losses during the Ms8.0 Wenchuan earthquake on May 12, 2008. Strain fluctuation curves were completely recorded at stress observatory stations in the Qinghai-Tibet plateau and its surroundings in the process of the earthquake. This paper introduces the geological background of the Wenchuan earthquake and the profile of in-situ stress monitoring stations. In particular, data of 174 earthquakes (Ms4.0-Ms8.5) were processed and analyzed with various methods, which were recorded at the Shandan station from August 2007 to December 2008. The results were compared with other seismic data, and further analyses were done for the recoded strain seismic waves, co-seismic strain stepovers, pre-earthquake strain valleys, Earth’s free oscillations before and after the earthquake and their physical implications. During the Wenchuan earthquake, the strainmeter recorded a huge extensional strain of 70 seconds, which shows that the Wenchuan earthquake is a rupture process predominated by thrusting. Significant precursory strain anomalies were detected 48 hours, 30 hours, 8 hours and 37 minutes before the earthquake. The anomalies are very high and their forms are very similar to that of the main shock. Similar anomalies can also be found in strain curves of other shocks greater than Ms7.0, indicating that such anomalies are prevalent before a great earthquake. In this paper, it is shown that medium aftershocks (Ms5.5-6.0) can also cause Earth’s free oscillations. Study of free oscillations is of great significance to understand the internal structure of the Earth and focal mechanisms of earthquakes and to recognize slow shocks, thus providing a scientific basis for the prevention and treatment of geological disasters and the prediction of future earthquakes.     

论地震效应与超长期天气预报的关系

大气环流所需能源除太阳辐射为主要能源外 ,地震亦是不容忽视的外在能源。据统计 ,每年全球发生震级M >5以上的地震约 10 0 0次 ,每次地震在震中附近有一个增温区 ,可向大气环流补充能量。特别是环太平洋地震带 ,这里既是地震多发区 ,又是台风多发区 ,地震能影响台风运动并影响大气环流运动。除日食能影响大气环流外 ,地震亦是一个不可忽视的因素。通过对 1993,1999年用日食效应及日食、地震效应数学模拟计算 ,充分证明地震是可以影响大气环流运动的新因素。     

A tectonophysical model of the Baikal seismic zone: testing and implications for medium-term earthquake prediction

The first tectonophysical model of the Baikal seismic zone represents a separate complex region of the lithosphere. It has a pinnate structure with a backbone belt of current deformation, which is a concentrator of largest earthquakes, and branching, repeatedly reactivated large and small faults. In its vertical section, the seismic zone is tree-like, the stem and the branches being faults of different size ranks which can generate earthquakes when reactivated. The real-time short-period fault motions and the respective seismicity occurring at a certain time and in certain places are triggered by strain waves, which disturb the metastable state of the faulted lithosphere subject to regional stress. The modeling work includes developing general requirements for tectonophysical models of continental rifts and special methods for identifying the faults that become active within short historic time spans, as well as techniques for locating potential events in space and time in specific active faults. The methods and model testing for medium-term earthquake prediction are described by the example of the well-documented Baikal seismic zone, which is the most active part of the Baikal rift system. The tectonophysical model for the Baikal zone is statistically supported by field data, and this allows estimating the velocities and periods of strain waves for different zone segments and faults, with implications for nearest-future earthquake prediction.     

Premonitory crustal deformations, strains and seismotectonic features (b-values) preceding Koyna earthquakes

There have been instances of premonitory variations in tilts, displacements, strains, telluric current, seismomagnetic effects, seismic velocities ( V_p, V_s) and their ratio (V_p/V_s), b-values, radon emission, etc. preceding large and moderate earthquakes, especially in areas near epicentres and along faults and other weak zones. Intensity and duration (T) of these premonitory quantities are very much dependent on magnitude (M) of the seismic event. Hence, these quantities may be utilised for prediction of an incoming seismic event well in advance of the actual earthquake. In the recent past, tilts, strain in deep underground rock and crustal displacements have been observed in the Koyna earthquake region over a decade covering pre- and postearthquake periods; and these observations confirm their reliability for qualitative as well as quantitative premonitory indices. Tilt began to change significantly one to two years before the Koyna earthquake of December 10, 1967, of magnitude 7.0. Sudden changes in ground tilt measured in a watertube tiltmeter accompanied an earthquake of magnitude 5.2 on October 17, 1973 and in other smaller earthquakes in the Koyna region, though premonitory changes in tilt preceding smaller earthquakes were not so much in evidence. However, changes in strains in deep underground rock were observed in smaller earthquakes of magnitude 4.0 and above. Furthermore, as a very large number of earthquakes (M = 1–7.0) were recorded in the extensive seismic net in the Koyna earthquake region during 1963–1975, precise b-value variations as computed from the above data, could reveal indirectly the state of crustal (tectonic) strain variations in the earthquake focal region and consequently act as a powerful premonitory index, especially for the significant Koyna earthquakes of December 10, 1967 (M = 7.0) and October 17, 1973 (M = 5.2). The widespread geodetic and magnetic levelling observations covering the pre- and postearthquake periods indicate significant vertical and horizontal crustal displacements, possibly accompanied by large-scale migration of underground magma during the large seismic event of December 10, 1967 in the Koyna region (M = 7.0). Duration (T) of premonitory changes in tilt, strains, etc., is generally governed by the equation of the type logT = A + BM (A and B are statistically determined coefficients). Similar other instances of premonitory evidences are also observed in micro-earthquakes (M = − 1 to 2) due to activation of a fault caused by nearby reservoir water-level fluctuations.     

A method for estimating earthquake occurrence probability using first- and multiple-order Markov chain models

It is valuable in earthquake prediction to determine the occurrence probability of major earthquakes by making use of data obtained from precursory phenomena up to the time of the evaluation. In this study, the time evolution of the state determined by earthquakes and precursory phenomena was modelled using Markov chains. Various probabilities suitable for earthquake prediction were derived from the transition probability of the Markov chain with a chosen length of memory time. As an example, earthquake sequence records for northern China, which covered a period of about two thousand years, were examined and the results were also obtained from modern scientific observations of the radon anomaly which covered a period of about 10 years.Assuming moderate-sized earthquakes of two magnitude ranges 4 3/4 M 5 3/4, 6 M 6 3/4 as precursors to large earthquakes of a magnitude range, 7 M 8 1/2, transition probabilities were calculated for a time interval from 2 to 40 yr. The results showed that the precursory time of moderate-sized earthquakes is mainly distributed in a time span of around several years and that the earthquake occurrence probability is considerably large when the precursory earthquakes occur successively. Furthermore, it was shown that a larger moderate-sized earthquake (4 3/4 M 5 3/4) is a more effective precursor than an entire moderate-sized earthquake (4 3/4 M 6 3/4).Second, a multiple precursor case was tried by means of simulation based on the radon anomaly data obtained during a limited observation period. Simultaneous occurrence of two precursors makes the earthquake occurrence probability increase by 1.5–2.0 for a reasonable choice of a mean recurrence interval of the radon anomaly compared with the case where only a moderate-sized earthquake was treated as the precursor. However, the probability is much the same if the average recurrence interval of the radon anomaly is the same through time, including the preparatory period before the earthquake.     

Numerical simulation of dynamic mechanisms of the 2008 Wenchuan Ms8.0 earthquake: implications for earthquake prediction

The sudden and unexpected Wenchuan earthquake (Ms = 8.0) occurred on the Longmen Shan Fault, causing a large number of casualties and huge property loss. Almost no definite precursors were reported prior to this event by Chinese scientists, who made a first successful prediction of the 1975 Haicheng earthquake (M = 7.3) in China. Does the unsuccessful prediction of the Wenchuan earthquake mean earthquake prediction is inherently impossible? In order to answer this question, the paper simulated inter- and co-seismic deformation, and recurrence of strong earthquakes associated with the Longmen Shan listric thrust fault by means of viscoelastic finite element method. The modeling results show that the computed interseismic strain accumulation in the lower crust beneath the Eastern Tibet is much faster than that in the other regions. In particular, the elastic strain energy density rate accumulates very rapid in and around the Longmen Shan fault in the depth above ~25 km that may explain why the great Wenchuan earthquake occurs in the region of such a slow surface deformation rate. The modeled coseismic displacements around the fault are consistent with surface rupture, aftershock distribution, and GPS measurement. Also, the model displays the slip history on the Longmen Shan fault, implying that the average earthquake recurrence interval on the Longmen Shan fault is very long, 3,300 years, which is in good agreement with the observed by paleoseismological investigations and estimates by other methods. Moreover, the model results indicate that the future earthquake could be evaluated based on numerical computation, rather than on precursors or on statistics. Numerical earthquake prediction (NEP) seems to be a promising avenue to a successful prediction, which will play an important part in natural hazard mitigation. NEP is difficult but possible, which needs well supporting.     

Earthquake prediction: 20 years of global experiment

Earthquake professionals have for many decades recognized the benefits to society from reliable earthquake predictions, but uncertainties regarding source initiation, rupture phenomena, and accuracy of both the timing and magnitude of the earthquake occurrence have oftentimes seemed either very difficult or impossible to overcome. The problem is that most of these methods cannot be adequately tested and evaluated either because of (a) lack of a precise definition of “prediction” and/or (b) shortage of data for meaningful statistical verification. This is not the case for the pattern recognition algorithm M8 designed in 1984 for prediction of great, Magnitude 8, earthquakes, hence its name. By 1986, the algorithm was rescaled for applications aimed at smaller magnitude earthquakes, down to M5+ range, and since then it has become a useful tool for systematic monitoring of seismic activity in a number of test seismic regions worldwide. After confirmed predictions of both the 1988 Spitak (Armenia) and the 1989 Loma Prieta (California) earthquakes, a “rigid test” to evaluate the efficiency of the intermediate-term middle-range earthquake prediction technique has been designed. Since 1991, each half-year, the algorithm M8 alone and in combination with its refinement MSc has been applied in a real-time prediction mode to seismicity of the entire Earth, and this test outlines, where possible, the areas in the two approximations where magnitude 8.0+ and 7.5+ earthquakes are most likely to occur before the next update. The results of this truly global 20-year-old experiment are indirect confirmations of the existing common features of both the predictability and the diverse behavior of the Earth’s naturally fractal lithosphere. The statistics achieved to date prove (with confidence above 99 %) rather high efficiency of the M8 and M8-MSc predictions limited to intermediate-term middle- and narrow-range accuracy. These statistics support the following general conclusions—(1) precursory seismic patterns do exist; (2) the size of an area where precursory seismic patterns show up is much larger than that of the source zone of the incipient target earthquake; (3) many precursory seismic patterns appear to be similar, even in regions of fundamentally different tectonic environments; and (4) some precursory seismic patterns are analogous to those in advance of extreme catastrophic events in other complex nonlinear systems (e.g., magnetic storms, solar flares, “starquakes”, etc.)—that are of high importance for further searches of the improved earthquake forecast/prediction algorithms and methods.     

A stochastic model of earthquake occurrence and the accompanying horizontal land deformation

A large-scale earthquake is believed to be associated with a release of strain energy accumulated in the crust, probably by the motion of upper-mantle lithosphere. Such an earthquake mechanism is well simulated by a belt-conveyer model proposed by Utsu (1972). The probability of earthquake occurrence can be estimated on the assumption that the motion of a slider on the belt-conveyer is mathematically formulated as a Markov process.In the probabilistic expressions, the results of Mogi's (1962) rock-fracture experiments are applied to the hazard-rate function of earthquake occurrence. The hazard-rate function has two coefficients, A and B, to be determined by the experiments. It is concluded that, when B is small, a number of small-scale earthquakes occur in the early time after the accumulation of crustal strain energy starts, but that the accumulated strain energy changes catastrophically into a single large-scale earthquake, when B is large.     

Statistical aspects of Parkfield earthquake sequence and Parkfield prediction experiment

This note discusses three interconnected statistical problems concerning the Parkfield sequence of moderate earthquakes and the Parkfield prediction experiment: (a) Is it possible that the quasi-periodic Parkfield sequence of characteristic earthquakes is no uncommon, specific phenomenon (the research hypothesis), but can be explained by a preferential selection from available earthquake catalogs? To this end we formulate the null hypothesis (earthquakes occur according to the Poisson process in time and their size follows the Gutenberg-Richter relation). We test whether the null hypothesis can be rejected as an explanation for the Parkfield sequence. (b) If the null hypothesis cannot be refuted, what is the probability of magnitude m ≥ 6 earthquake occurrence in the Parkfield region? (c) The direct goal of the Parkfield experiment is the registration of precursory phenomena prior to a m6 earthquake. However, in the absence of the characteristic earthquake, can the experiment resolve which of the two competing hypotheses is true in a reasonable time? Statistical analysis is hindered by an insufficiently rigorous definition of the research model and inadequate or ambiguous data. However, we show that the null hypothesis cannot be decisively rejected. The quasi-periodic pattern of intermediate size earthquakes in the Parkfield area is a statistical event likely to occur by chance if it has been preferentially selected from available earthquake catalogs. The observed magnitude-frequency curves for small and intermediate earthquakes in the Parkfield area agree with the theoretical distribution computed on the basis of a modified Gutenberg-Richter law (gamma distribution), using deformation rates for the San Andreas fault. We show that the size distribution of the Parkfield characteristic earthquakes can also be attributed to selection bias. According to the null hypothesis, the yearly probability of a m ≥ 6 earthquake originating in the Parkfield area is less than 1%, signifying that several more decades of observation may be needed before the expected event occurs. By its design, the Parkfield experiment cannot be expected to yield statistically significant conclusions on the validity of the research hypothesis for many decades.     

Increased correlation range of seismicity before large events manifested by earthquake chains

“Earthquake chains” are clusters of moderate-size earthquakes which extend over large distances and are formed by statistically rare pairs of events that are close in space and time (“neighbors”). Earthquake chains are supposed to be precursors of large earthquakes with lead times of a few months. Here we substantiate this hypothesis by mass testing it using a random earthquake catalog. Also, we study stability under variation of parameters and some properties of the chains. We found two invariant parameters: they characterize the spatial and energy scales of earthquake correlation. Both parameters of the chains show good correlation with the magnitudes of the earthquakes they precede. Earthquake chains are known as the first stage of the earthquake prediction algorithm reverse tracing of precursors (RTP) now tested in forward prediction. A discussion of the complete RTP algorithm is outside the scope of this paper, but the results presented here are important to substantiate the RTP approach.     

Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, Japan

We found a characteristic space–time pattern of the tidal triggering effect on earthquake occurrence in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, central Japan, where a large interplate earthquake may be impending. We measured the correlation between the Earth tide and earthquake occurrence using microearthquakes that took place in the Philippine Sea plate for about two decades. For each event, we assigned the tidal phase angle at the origin time by theoretically calculating the tidal shear stress on the fault plane. Based on the distribution of the tidal phase angles, we statistically tested whether they concentrate near some particular angle or not by using Schuster's test. In this test, the result is evaluated by p-value, which represents the significance level to reject the null hypothesis that earthquakes occur randomly irrespective of the tidal phase angle. As a result of analysis, no correlation was found for the data set including all the earthquakes. However, we found a systematic pattern in the temporal variation of the tidal effect; the p-value significantly decreased preceding the occurrence of M ≥ 4.5 earthquakes, and it recovered a high level afterwards. We note that those M ≥ 4.5 earthquakes were considerably larger than the normal background seismicity in the study area. The frequency distribution of tidal phase angles in the pre-event period exhibited a peak at the phase angle where the tidal shear stress is at its maximum to accelerate the fault slip. This indicates that the observed small p-value is a physical consequence of the tidal effect. We also found a distinctive feature in the spatial distribution of p-values. The small p-values appeared just beneath the strongly coupled portion of the plate interface, as inferred from the seismicity rate change in the past few years.