Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

Seismic activity associated with the subducting motion of the Philippine Sea plate beneath the Kanto district, Japan

The Pacific plate and the Philippine Sea plate overlap and subduct underneath the Kanto region, central Japan, causing complex seismic activities in the upper mantle. In this research, we used a map selection tool with a graphic display to create a data set for earthquakes caused by the subducting motion of the Philippine Sea plate that are easily determined. As a result, we determined that there are at least four earthquake groups present in the upper mantle above the Pacific plate. Major seismic activity (Group 1) has been observed throughout the Kanto region and is considered to originate in the uppermost part of mantle in the subducted Philippine Sea plate, judging from the formation of the focal region and comparison with the 3D structure of seismic velocity. The focal mechanism of these earthquakes is characterized by the down-dip compression. A second earthquake layer characterized by down-dip extension (Group 2), below the earthquakes in this group, is also noted. The focal region for those earthquakes is considered to be located at the lower part of the slab mantle, and the Pacific plate located directly below is considered to influence the activity. Earthquakes located at the shallowest part (Group 3) form a few clusters distributed directly above the Group 1 focal region. Judging from the characteristics of later phases in these earthquakes and comparing against the 3D structure of seismic velocity, the focal regions for the earthquakes are considered to be located near the upper surface of the slab. Another earthquake group (Group 4) originates further below Group 2; it is difficult to consider these earthquakes within a single slab. The seismic activities representing the upper area of the Philippine Sea plate are Group 3. This paper proposes a slab geometry model that is substantially different from conventional models by strictly differentiating the groups.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

Earthquake generation in Cyprus revealed by the evolving stress field

Strong earthquake occurrence (M ≥ 6.0) onshore and offshore the Cyprus Island constitutes significant seismic hazard because they occur close to populated areas. Seismicity is weak south of the Island along the Cyprean Arc and strong events are aligned along the Paphos transform fault and Larnaka thrust fault zone that were already known and the Lemessos thrust fault zone that defined in the present study. By combining the past history of strong (M ≥ 6.0) events and the long-term tectonic loading on these major fault zones, the evolution of the stress field from 1896 until the present is derived. Although uncertainties exist in the location, magnitude and fault geometries of the early earthquakes included in our stress evolutionary model, the resulting stress field provides an explanation of later earthquake triggering. It was evidenced that the locations of all the strong events were preceded by a static stress change that encouraged failure. The current state of the evolved stress field may provide evidence for the future seismic hazard. Areas of positive static stress changes were identified in the southwestern offshore area that can be considered as possible sites of future seismic activity.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

冲击凿岩钎杆的瞬态动力学分析

以ANSYS为平台，按照真实工作状态对凿岩钎杆进行了瞬态动力学分析，并对不同活塞形状下分析结果数据进行了比较，为冲击凿岩系统设计提供了理论依据。     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

Geological and geomorphological characteristics of landslides triggered by the 2004 Mid Niigta prefecture earthquake in Japan

The 2004 Mid Niigta prefecture earthquake (M_JMA 6.8) triggered more than one thousand landslides in the Miocene to Quaternary sedimentary rocks in Japan. The most common landslides were shallow disrupted landslides on steep slopes, which has been common in many previous disastrous earthquakes in the world. The Mid Niigta prefecture earthquake also triggered more than one hundred deep landslides, providing valuable information on the conditions for their occurrence. A field investigation and the interpretation of aerial photographs taken before and after the earthquake suggest that reactivation of existing landslides and undercutting of slopes are the most important factors for deep landslides to be triggered by earthquakes. In addition, planar sliding surfaces seem to be essential for the generation of catastrophic landslides triggered by this earthquake. Planar bedding–parallel sliding surfaces were formed at the boundary between the overlying permeable sandstone and underlying siltstone or along the bedding planes of alternating beds of sandstone and siltstone. Sliding surfaces along the slope-parallel oxidation front were formed in the area of black mudstone. New landslides (rockslide-avalanches) occurred with the sliding surfaces in a several-cm thick tuff interbedded in siltstone. One rockslide-avalanche occurred on a slope where buckling deformation preceded the earthquake. Gentle valley bottom sediments were mobilized in many locations, probably because they were saturated and partial liquefaction had occurred due to the earthquake shaking.     

Mineralogy, geochemistry, and Nd isotope composition of the Rainbow hydrothermal field, Mid-Atlantic Ridge

Petrological, geochemical, and Nd isotopic analyses have been carried out on rock samples from the Rainbow vent field to assess the evolution of the hydrothermal system. The Rainbow vent field is an ultramafic-hosted hydrothermal system located on the Mid-Atlantic Ridge characterized by vigorous high-temperature venting (∼365°C) and unique chemical composition of fluids: high chlorinity, low pH and very high Fe, and rare earth element (REE) contents (Douville et al., Chemical Geology 184:37–48, 2002). Serpentinization has occurred under a low-temperature (<270°C) retrograde regime, later overprinted by a higher temperature sulfide mineralization event. Retrograde serpentinization reactions alone cannot reproduce the reported heat and specific chemical features of Rainbow hydrothermal fluids. The following units were identified within the deposit: (1) nonmineralized serpentinite, (2) mineralized serpentinite—stockwork, (3) steatite, (4) semimassive sulfides, and (5) massive sulfides, which include Cu-rich massive sulfides (up to 28wt% Cu) and Zn-rich massive sulfide chimneys (up to 5wt% Zn). Sulfide mineralization has produced significant changes in the sulfide-bearing rocks including enrichment in transition metals (Cu, Zn, Fe, and Co) and light REE, increase in the Co/Ni ratios comparable to those of mafic Cu-rich volcanic-hosted massive sulfide deposits and different ¹⁴³Nd/¹⁴⁴Nd isotope ratios. Vent fluid chemistry data are indicative of acidic, reducing, and high temperature conditions at the subseafloor reaction zone where fluids undergo phase separation most likely under subcritical conditions (boiling). An explanation for the high chlorinity is not straightforward unless mixing with high salinity brine or direct contribution from a magmatic Cl-rich aqueous fluid is considered. This study adds new data, which, combined with the current knowledge of the Rainbow vent field, brings compelling evidence for the presence, at depth, of a magmatic body, most likely gabbroic, which provides heat and metals to the system. Co/Ni ratios proved to be good tools used to discriminate between rock units, degree of sulfide mineralization, and positioning within the hydrothermal system. Deeper units have Co/Ni <1 and subsurface and surface units have Co/Ni >1.     

Damage Survey of Water Supply Systems and Fragility Curve of PVC Water Pipelines in the Chi–Chi Taiwan Earthquake

Substantial damage to water supply systems, including water delivery pipelines, water treatment plants, reservoirs, and water storage tanks, was reported after the 1999 Chi–Chi Taiwan Earthquake. This paper first summarizes the damage survey and then presents the results of seismic fragility analysis for underground pipelines. Construction blueprints of the water delivery pipelines and repair work orders of 11 townships and cities in the disastrous area were digitized into a Geographical Information System (GIS) for analysis and assessment. With the aid of the GIS system, we found that PVC pipes made up 86% of water delivery pipelines while steel, cast iron, ductile iron, PE and others took the rest. Therefore, this paper focuses on the fragility analysis of PVC pipes. Three different methods were applied to derive the fragility relations between the PVC water pipes having nominal diameters (approximately inner diameters) greater than or equal to 65 mm and earthquake intensity parameters such as peak ground acceleration and peak ground velocity. The results were then examined with those of other countries. The discrepancy between our results and the empirical equation used by HAZUS, an earthquake loss estimation software developed by the Federal Emergency Management Agency was not significant.