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

Introduction The unexpected December 26, 2004, off the west coast of northern Sumatra, Indonesia, MW=9.0 earthquake, which caused devastating tsunami around the Indian Ocean, reminds seis-mologists of the difficulty of earthquake forecast and/or prediction. For seismologists this earth-quake is almost completely unexpected, because there was neither forecasting (which means the estimation of the future earthquake rate as a function of location, time, and magnitude) nor predic-tion (forecasti…     

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.     

印尼苏门答腊西8.7级地震的几个特点及由其引发的一些思考

分析总结了印尼苏门答腊西8．7级地震的几个特点：(1)板缘特大地震;(2)引发海啸,受灾范围特别大,成灾的瞬间突发性却不如内陆地震明显;(3)在远处出现“湖面波动”等同震现象突出;(4)震前全球地震活动出现异常图像;(5)安达曼弧一带与川滇地区强震活动存在相关性。并由此提出了：尽快加紧地震海啸预警机制研究及其技术系统的建立;加强地震预测预警机制研究;要重点研究川滇地区的地震危险性;强化防震减灾科普宣传等建议。     

The August 27, 2008, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 6.3, Kultuk earthquake effects in the Nearepicenter zone: Macroseismic survey results

In this work, the macroseismic effects of the Kultuk earthquake (M _W = 6.3), which took place on August 27, 2008 in the southwestern closure of the Baikal Lake, are under consideration. The intensity of shocks in inhabited localities located in the epicentral zone reached 7–8 points on the MSK-64 scale. The earthquake was named after the local settlement of Kultuk, which was the mostly damaged area by the earthquake. The considered seismic event caused significant material damage (about 250 million rubles according to preliminary estimates). In inhabited localities of Southern Pribaikalie thousands of cases of damaged ovens and chimneys were registered. Some buildings were highly damaged and, accordingly, they are unfit for further use. The earthquake was the cause of numerous rockslides, rockfalls, and landslides on steep natural and artificial slopes. A macroseismic survey allowed us to establish the asymmetrical distribution of the intensity of shocks relative to the epicenter.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> and <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript> comparison for earthquakes in Greece for the period 1969 to 2007

We use 576 earthquakes of magnitude, M _w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M _w, and local magnitude, M _L, for earthquakes in Greece and the surrounding regions. We compare M _w to M _L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M _L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M _w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M _L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M _w and M _L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M _L’s. We propose scaling relations to convert the M _L values reported from the two centers to M _w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.     

Gravity variation before Kunlun mountain pass western <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 earthquake

The relation between the gravity variation features and M _S=8.1 earthquake in Qinghai-Xizang monitoring area is analyzed preliminarily, by using spatial dynamic variation results of regional gravity field from absolute gravity and relative gravity observation in 1998 and 2000. The results show that: 1) M _S=8.1 earthquake in Kulun mountain pass western occurred in the gravity variation high gradient near gravity’s high negative variation; 2) The main tectonic deformation and energy accumulation before M _S=8.1 earthquake are distributed at south side of the epicenter; 3) The range of gravity’s high negative variation at east of the M _S=8.1 earthquake epicenter relatively coincides with that rupture region according to field geology investigation; 4) Gravity variation distribution in high negative value region is just consistent with the second shear strain’s high value region of strain field obtained from GPS observation.     

GPS time-series and its response to <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 Kunlunshan earthquake

In this paper, observation data in 25 GPS reference stations of China have been analyzed by calculating GPS position coordinate time-series with GIPSY. Result shows there is an obvious trend variation in such time-series. The trend variations of time series along the longitude and latitude coordinate reflect the motion of each position in the global-plate, in which the trend variation in the vertical direction reveals some large-scale construction information or reflects the local movement around the positions. The analysis also shows that such time-series have a variation cycle of nearly 1.02 a, but the reason still remains to be further studied. At the end of this paper, response of the time-series of M _S=8.1 Kunlunshan earthquake was analyzed, and the seismogenic process of M _S=8.1 Kunlunshan earthquake, according to the time proceeding and the feature of anomaly, was divided into 3 phases—changes in blocks with forces, strain accumulation, quick accumulation and slow release of energy. At the initial stage of seismogenic process of M _S=8.1 earthquake and at the imminent earthquake, coseismic process as well as during the post earthquake recovery, anomaly in vertical direction is always in a majority. The anomalous movement in vertical direction at the initial stage resulted in a blocking between faults, while at the middle stage of seismogenic process, the differential movement between blocks are in a majority, which is the major reason causing energy accumulating at the blocking stage of faults.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>9 Tohoku earthquake in 2011 in Japan: precursors uncovered by natural time analysis

This paper reviews the precursory phenomena of the 2011 M _W9 Tohoku earthquake in Japan that emerge solely when we analyze the seismicity data in a new time domain termed natural time. If we do not consider this analysis, important precursory changes cannot be identified and hence are missed. Natural time analysis has the privilege that enables the introduction of an order parameter of seismicity. In this frame, we find that the fluctuations of this parameter exhibit an unprecedented characteristic change, i.e., an evident minimum, approximately two months before Tohoku earthquake, which strikingly is almost simultaneous with unique anomalous geomagnetic field variations recorded mainly on the z component. This is consistent with our finding that such a characteristic change in seismicity appears when a seismic electric signal (SES) activity of the VAN method (from the initials of Varotsos, Alexopoulos, Nomicos) initiates, and provides a direct confirmation of the physical interconnection between SES and seismicity.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

The earthquake of July 22, 2011 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 4.5) in a low-seismicity area of the Argun region

The paper considers the Argun earthquake of July 22, 2011 (M _w = 4.5), which occurred in the Argun River valley in a low-seismicity territory in China. The focal parameters of the earthquake (depth of the hypocenter, moment magnitude, scalar seismic moment, and focal mechanism) were determined by calculating the seismic moment tensor from the amplitude spectra of surface waves and the data on the signs of the first arrivals of body waves at regional stations. The solution of the focal mechanism makes it possible to assume a relationship between the earthquake focus and a fault with a northeastern strike bordering the southeastern side of the Argun Basin (in Chinese territory). The Argun earthquake was felt in Russia with an intensity of II–III to V at the epicentral distances up to 255 km. The intensity of shaking did not exceed values suggested by new GSZ-2012 and GSZ-2014 seismic zoning maps of Russian territory. Nevertheless, the question on the possible occurrence of stronger earthquakes in the studied region remains open.     

Study on a gravity precursor mode of the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>=7.0 Lijiang earthquake

To study the seismogeny process or the precursory behavior of the 1976 M _S=7.0 Lijiang earthquake, we analyze the repeat gravity data with high precision from the Western Yunnan Earthquake Prediction Experiment Area (WYEPEA) and the related results of geology and geophysics survey in this paper. Considering the gross errors caused by observation data and model difference, we have firstly inverted the slip distribution of the main active faults with time based on the robust Bayesian least squares estimation and multi-fault dislocation model. The results show that the slip changes of the faults with time from 1990 to 1997 obviously reflect the preparation process of the Lijiang earthquake. The images of main precursor mode have the characteristic of main shock-after shock type, which is agreement to the model of coupling movement between crust density and crust deformation (DD mode of coupling movement).     

The 2014 Kefalonia Doublet (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>6.1 and <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>6.0), Central Ionian Islands,Greece: Seismotectonic Implications along the Kefalonia Transform Fault Zone

The 2014 Kefalonia earthquake sequence started on 26 January with the first main shock (M_W6.1) and aftershock activity extending over 35 km, much longer than expected from the causative fault segment. The second main shock (M_W6.0) occurred on 3 February on an adjacent fault segment, where the aftershock distribution was remarkably sparse, evidently encouraged by stress transfer of the first main shock. The aftershocks from the regional catalog were relocated using a 7-layer velocity model and station residuals, and their distribution evidenced two adjacent fault segments striking almost N-S and dipping to the east, in full agreement with the centroid moment tensor solutions, constituting segments of the Kefalonia Transform Fault (KTF). The KTF is bounded to the north by oblique parallel smaller fault segments, linking KTF with its northward continuation, the Lefkada Fault.     

Estimating the source parameters for the M<Subscript><Emphasis Type="Italic">W</Emphasis></Subscript> 7.6 April 20, 2006 Olyutorskii earthquake from long period <Emphasis Type="Italic">P</Emphasis>-wave records made by the worldwide network

The source parameters of the M _W = 7.6 Olyutorskii earthquake were estimated using the moments of the slip rate function with degrees 1 and 2. The moments were estimated from broadband P-wave records at 52 stations of the worldwide network. The first step was to find a function S(t) for each station; this function is an apparent source time function, i.e., the P-wave slip as radiated by the source toward a station under consideration. The method of empirical Green’s functions was used to estimate S(t). The next step was to calculate the moments of S(t) of degrees 1 and 2 over time and to set up relevant equations to be solved by least squares for the unknown source moments. The horizontal linear source was used as a nonparametric model for calculating the source moments. Haskell’s parametric model was used for further interpretation of the source moments. The resulting estimates are as follows: the source centroid was 13–25 km southwest of the epicenter, the source was 105–120 km long, the source strike was 222°–228°, the rupture velocity was 2.7–3.0 km/s, and the total radiation duration was 24–27 s. These estimates indicate a bilateral rupture dominated by a southwestward sense of rupture propagation. The source characteristics are consistent with the aftershock area geometry and with the focal mechanism, as well as with surface breakage as observed by geologists in the field.     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

Seismic damage to structures in the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>6.5 Ludian earthquake

On 3 August 2014, the Ludian earthquake struck northwest Yunnan Province with a surface wave magnitude of 6.5. This moderate earthquake unexpectedly caused high fatalities and great economic loss. Four strong motion stations were located in the areas with intensity V, VI, VII and IX, near the epicentre. The characteristics of the ground motion are discussed herein, including 1) ground motion was strong at a period of less than 1.4 s, which covered the natural vibration period of a large number of structures; and 2) the release energy was concentrated geographically. Based on materials collected during emergency building inspections, the damage patterns of adobe, masonry, timber frame and reinforced concrete (RC) frame structures in areas with different intensities are summarised. Earthquake damage matrices of local buildings are also given for fragility evaluation and earthquake damage prediction. It is found that the collapse ratios of RC frame and confined masonry structures based on the new design code are significantly lower than non-seismic buildings. However, the RC frame structures still failed to achieve the ‘strong column, weak beam’ design target. Traditional timber frame structures with a light infill wall showed good aseismic performance.     

The seismotectonic and macroseismic features of the Wenchuan (Sichuan) earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript> = 8.0) of May 12, 2008

A disastrous earthquake with a magnitude M _S = 8.0 (M _W = 7.9), in China called “the 5.12 Wenchuan earthquake,” occurred on May 12, 2008, in Sichuan province on the border between the Sino-Tibetan Mountains and the Sichuan depression. The instrumental epicenter was registered in the southeastern part of Wenchuan county, and the hypocenter depth was 14 km. As the strongest and most destructive earthquake within mainland China, it caused numerous human losses and destruction of buildings and infrastructure. The seismic effect from the main shock and aftershocks was felt in many counties, towns, and villages, though Sichuan province suffered the most. The maximum intensity of the shocks was estimated at 11 degrees, according to the Chinese macroseismic scale. In the process of source opening, from the southern part of Wenchuan county to the vicinities of Quingchuan, a seismic fault system with a total length up to 240 km out-cropped on the earth’s surface, confined to the Longmenshan fault belt. The seismic fault system disturbed the original ground, resulting in the collapse or damage to various constructions, such as buildings, homes, bridges, roads, etc. Fault offsets had a dextral strike-slip and thrust kinematic combination. The earthquake generated several tens of thousands of landslides, rockfalls, and debris flows. Many dammed ponds appeared in the epicentral zone due to the activation of landslides. Thus, the geological effects turned out to be the most destructive factor in this case. At the same time, the seismic intensity of surface shaking was abnormally low even in direct proximity to the seismic fault system. Usually it was no more than 7–8 degrees. This macroseismic phenomenon may turn out to be rather typical for many major earthquakes.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

Application of Fisher method to discriminating earthquakes and explosions using criterion <Emphasis Type="Italic">m</Emphasis><Subscript>b</Subscript>/<Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>

We try to give a quantitative and global discrimination function by studying m _b/M _S data using Fisher method that is a kind of pattern recognition methods. The reliability of the function is also analyzed. The results show that this criterion works well and has a global feature, which can be used as first-level filtering criterions in event identification. The quantitative and linear discrimination function makes it possible to identify events automatically and achieve the goal to react the events quickly and effectively.     

Rupture behavior and deformation localization of the Kunlunshan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>7.8) and their tectonic implications

Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M _w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M _w =6.8, M _w =6.2 and M _w ⩽=7.8 events, respectively. The M _w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)