Water level fluctuations due to earthquakes in Koyna-Warna region,India

Earthquakes cause a variety of hydrological phenomena, including changes in the ground water levels in bore wells. The Koyna region in the peninsular shield of India, hitherto considered stable in terms of seismic activity, has been active since 1967. More recently, the earthquakes have been localized to the newly impounded Warna reservoir, which is located south of Koyna, where a burst of seismicity occurred in 1993. The region continues to remain seismically active even after four decades. Twenty-one bore wells were drilled around the seismic source volume in the region to observe water level changes resulting from earthquake phenomena. Our studies have shown coseismic anomalous water level changes to be associated with the moderate earthquakes of April 25, 1997 and February 11, 1998. Our results show that changes in the ground water level in bore wells are correlated with micro-earthquake activity, both preceding and following moderately sized earthquakes. The results have implications in enhancing our understanding of earthquake mechanisms.     

Thessaloniki–Gerakarou Fault Zone (TGFZ): the western extension of the 1978 Thessaloniki earthquake fault (Northern Greece) and seismic hazard assessment

Active faulting and seismic properties are re-investigated in the eastern precinct of the city of Thessaloniki (Northern Greece), which was seriously affected by two large earthquakes during the 20th century and severe damage was done by the 1759 event. It is suggested that the earthquake fault associated with the occurrence of the latest destructive 1978 Thessaloniki earthquake continues westwards to the 20-km-long Thessaloniki–Gerakarou Fault Zone (TGFZ), which extends from the Gerakarou village to the city of Thessaloniki. This fault zone exhibits a constant dip to the N and is characterised by a complicated geometry comprised of inherited 100°-trending faults that form multi-level branching (tree-like fault geometry) along with NNE- to NE-trending faults. The TGFZ is compatible with the contemporary regional N–S extensional stress field that tends to modify the pre-existing NW–SE tectonic fabric prevailing in the mountainous region of Thessaloniki. Both the 1978 earthquake fault and TGFZ belong to a ca. 65-km-long E–W-trending rupture fault system that runs through the southern part of the Mygdonia graben from the Strymonikos gulf to Thessaloniki. This fault system, here called Thessaloniki–Rentina Fault System (TRFS), consists of two 17–20-km-long left-stepping 100°-trending main fault strands that form underlapping steps bridged by 8–10-km-long ENE–WSW faults. The occurrence of large (M6.0) historical earthquakes (in 620, 677 and 700 A.D.) demonstrates repeated activation, and therefore the possible reactivation of the westernmost segment, the TGFZ, could be a major threat to the city of Thessaloniki. Changes in the Coulomb failure function (ΔCFF) due to the occurrence of the 1978 earthquake calculated out in this paper indicate that the TGFZ has been brought closer to failure, a convincing argument for future seismic hazard along the TGFZ.     

The 2003 earthquake swarm in Anjalankoski, south-eastern Finland

During May 2003 a swarm of 16 earthquakes (M_L = 0.6–2.1) occurred at Anjalankoski, south-eastern Finland. The activity lasted for three weeks, but additional two events were observed at the same location in October 2004. A comparison of the waveforms indicated that the source mechanisms and the hypocentres of the events were nearly identical.A relative earthquake location method was applied to better define the geometry of the cluster and to identify the fault plane associated with the earthquakes. The relocated earthquakes aligned along an ENE–WSW trending zone, with a lateral extent of about 1.0 km by 0.8 km. The relative location and the waveform-modelling of depth sensitive surface wave (Rg) and S-to-P converted body wave (sP) phases indicated that the events were unusually shallow, most likely occurring within the first 2 km of the surface. The revised historical earthquake data confirm that shallow swarm-type seismicity is characteristic to the area.The focal mechanism obtained as a composite solution of the five strongest events corresponds to dip-slip motion along a nearly vertical fault plane (strike = 250°, dip = 80°, rake = 90°). The dip and strike of this nodal plane as well as the relocated hypocentres coincide with an internal intrusion boundary of the Vyborg rapakivi batholith.The events occur under a compressive local stress field, which is explained by large gravitational potential energy differences and ridge-push forces. Pore-pressure changes caused by intrusion of ground water and/or radon gas into the fracture zones are suggested to govern the swarm-type earthquake activity.     

Ultrashallow seismic experiment on a trenched section of the Chelunpu fault zone, Taiwan

Ultrashallow P-wave seismic reflection experiments were conducted at a model test site and in a trenched shallow fault zone along the Chelunpu fault line. The field layout was designed to have the shallowest undistorted reflection from about 1 m depth with 0.5 m vertical resolution. The smallest group interval tested in this study was 0.05 m with a 0.25 ms sample interval, which can avoid spatial aliasing of ground roll if the target is very shallow and the velocities are low. Data processing was designed to be simple but consistent. As the ultrashallow reflections may be contaminated with high-amplitude coherent noise in many aspects, first break muting and surgical muting were performed on each file as detailed as possible, and f–k filtering was applied mainly for the purpose of attenuating the aliasing energy and back-scattered noise. Data acquired in this study show that the low P-wave velocities (< 200 m/s) and high dominant frequencies (120–200 Hz) of near-surface layers may have a potential vertical resolution of 0.4 m or even better.Comparing the test profile with the ground-penetrating radar (GPR) control profile of the same test site and correlating the results obtained from the study site with those of the geologic cross-section of the trench, this experiment demonstrates the possibility of using seismic methods in investigating shallow structures at depths of less than a few meters with vertical resolution comparable to the GPR technique.     

Improved seismic velocity reference model from local earthquake data in Northwestern Italy

Resolution and reliability estimates of results obtained by seismic tomography strongly depend on the reference model. Inadequate initial reference models may severely distort tomographic images or introduce artefacts that lead to misinterpretations of the results. Reference models are usually obtained by means of a priori near-surface geological information or by geophysical information derived by controlled-source seismology.
Starting from the idea that a reference model must approximate the weighted average of data selected for the three-dimensional (3D) inversion, one-dimensional (1D) model for Northwestern Italy is derived that is able to minimize mean of RMS of a set of well-locatable earthquakes, by computing a solution of the coupled hypocentre 1D velocity problem.
Such a model, termed the Minimum 1D model, can be used both as an initial reference model for 3D inversion and as a reference velocity model for high-quality routine earthquake location.     

Temporal variation in the geometry of a strike–slip fault zone: Examples from the Dead Sea Transform

The location of the active fault strands along the Dead Sea Transform fault zone (DST) changed through time. In the western margins of Dead Sea basin, the early activity began a few kilometers west of the preset shores and moved toward the center of the basin in four stages. Similar centerward migration of faulting is apparent in the Hula Valley north of the Sea of Galilee as well as in the Negev and the Sinai Peninsula. In the Arava Valley, seismic surveys reveal a series of buried inactive basins whereas the current active strand is on their eastern margins. In the central Arava the centerward migration of activity was followed by outward migration with Pleistocene faulting along NNE-trending faults nearly 50 km west of the center. Largely the faulting along the DST, which began in the early–middle Miocene over a wide zone of up to 50 km, became localized by the end of the Miocene. The subsidence of fault-controlled basins, which were active in the early stage, stopped at the end of the Miocene. Later during the Plio-Pleistocene new faults were formed in the Negev west of the main transform. They indicate that another cycle has begun with the widening of the fault zone. It is suggested that the localization of faulting goes on as long as there is no change in the stress field. The stresses change because the geometry of the plates must change as they move, and consequently the localization stage ends. The fault zone is rearranged, becomes wide, and a new localization stage begins as slip accumulates. It is hypothesized that alternating periods of widening and narrowing correlate to changes of the plate boundaries, manifest in different Euler poles.     

Long-term earthquake prediction in the North Pacific seismic zone based on the time and magnitude predictable model

Earthquake recurrence intervals for large and great shallow mainshocks in 12 seismogenic sources along the North Pacific seismic zone (Alaska-Aleutians-Kamchatka-Kuril Islands) have been estimated and used for the determination of the following relations:

Tectonic stress field in the epicentral zone of the Latur earthquake of 1993

In situ stress measurements by hydraulic fracturing were carried out in the 617 m deep borehole specially drilled in the epicentral zone of the 1993 Latur earthquake for the purpose of research. The stress measurements carried out at 592 m depth in this borehole are the deepest of all such measurements made so far in the Indian shield. The maximum and minimum principal horizontal stresses (S _H max andS _h min) have been derived from the hydrofracture data using the classical method. TheS _{H max} andS _h min are found to be 16.5 and 9.6 MPa at 373 m depth, and 25.0 and 14.1 MPa at 592 m depth, indicating that the vertical gradients ofS _hmax andS _hmin in the epicentral zone are 39 MPa/km and 21 MPa/km respectively. The principal horizontal stresses in the epicentral zone are comparable with those at Hyderabad and 30% higher than in most other comparable intra-continental regions. Analysis of the results indicate that the stresses in the focal region of the 1993 Latur earthquake have not undergone any significant change following its occurrence and this is in agreement with a similar inference drawn from the seismic data analysis. It appears that the Latur earthquake was caused due to rupturing of the overpressured fault segment at the base of the seismogenic zone.     

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

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

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.     

Geochemistry of soil gas in the seismic fault zone produced by the Wenchuan Ms 8.0 earthquake, southwestern China

The spatio-temporal variations of soil gas in the seismic fault zone produced by the 12 May 2008 Wenchuan Ms 8.0 earthquake were investigated based on the field measurements of soil gas concentrations after the main shock. Concentrations of He, H₂, CO₂, CH₄, O₂, N₂, Rn, and Hg in soil gas were measured in the field at eight short profiles across the seismic rupture zone in June and December 2008 and July 2009. Soil-gas concentrations of more than 800 sampling sites were obtained. The data showed that the magnitudes of the He and H₂ anomalies of three surveys declined significantly with decreasing strength of the aftershocks with time. The maximum concentrations of He and H₂ (40 and 279.4 ppm, respectively) were found in three replicates at the south part of the rupture zone close to the epicenter. The spatio-temporal variations of CO₂, Rn, and Hg concentrations differed obviously between the north and south parts of the fault zone. The maximum He and H₂ concentrations in Jun 2008 occurred near the parts of the rupture zone where vertical displacements were larger. The anomalies of He, H₂, CO₂, Rn, and Hg concentrations could be related to the variation in the regional stress field and the aftershock activity.     

四川汶川5.12大地震同震滑动断层泥的发现及意义

2008年汶川8.0级地震沿龙门山断裂带内的映秀—北川断裂和灌县—安县断裂产生了近300 km的同震地表破裂带。震后地质科学考察发现地表变形以逆冲为主,并伴有右旋走滑。地震地表破裂带大多沿古生代碳质泥岩、页岩和三叠系煤系地层内的滑动面出露地表,这些软弱地层为地震破裂带冲到地表提供了超低摩擦滑动带。我们发现在同震垂直和水平位错达6m左右的地表破裂带,地震的同震滑动发生在厚度约0.5~2cm 的狭窄滑动带内,以发育新鲜的灰色断层泥为特征,这些断层泥是地震断层快速滑动过程中岩石—流体相互作用的结果。     

Devonian-Cretaceous repeated subsidence and uplift along the Teisseyre-Tornquist zone in SE Poland — Insight from seismic data interpretation

The Teisseyre-Tornquist Zone that separates the East European Craton from the Palaeozoic Platform forms one of the most fundamental lithospheric boundaries in Europe. Devonian to Cretaceous-Paleogene evolution of the SE segment of this zone was analyzed using high-quality seismic reflection data that provided detailed information regarding entire Palaeozoic and Mesozoic sedimentary cover, with particular focus on problems of Late Carboniferous and Late Cretaceous-Paleogene basin inversion and uplift. Two previously proposed models of development and inversion of the Devonian-Carboniferous Lublin Basin seem to only partly explain configuration of this sedimentary basin. A new model includes Late Devonian-Early Carboniferous reverse faulting within the cratonic area NE from the Kock fault zone, possibly first far-field effect of the Variscan orogeny. This was followed by Late Carboniferous inversion of the Lublin Basin. Inversion tectonics was associated with strike-slip movements along the Ursynów-Kazimierz fault zone, and thrusting along the Kock fault zone possibly triggered by deeper strike-slip movements. Late Carboniferous inversion-related deformations along the NE boundary of the Lublin Basin were associated with some degree of ductile (quasi-diapiric) deformation facilitated by thick series of Silurian shales. During Mesozoic extension and development of the Mid-Polish Trough major fault zones within the Lublin Basin remained mostly inactive, and subsidence centre moved to the SW, towards the Nowe Miasto-Zawichost fault zone and further to the SW into the present-day Holy Cross Mts. area. Late Cretaceous-Paleogene inversion of the Mid-Polish Trough and formation of the Mid-Polish Swell was associated with reactivation of inherited deeper fault zones, and included also some strike-slip faulting. The study area provides well-documented example of the foreland plate within which repeated basin inversion related to compressive/transpressive deformations was triggered by active orogenic processes at the plate margin (i.e. Variscan or Carpathian orogeny) and involved important strike-slip reactivation of crustal scale inherited fault zones belonging to the Teisseyre-Tornquist Zone.     

Rupture mechanism of the 1993 Killari earthquake, India: constraints from aftershocks and static stress change

The Killari earthquake of September 29, 1993 (Mw=6.2) in peninsular India triggered several aftershocks that were recorded by a network of 21 stations. We computed the change in regional static stress caused by coseismic slip on the earthquake rupture and correlated it with the aftershocks with a view to constrain some of the rupture parameters of this earthquake. We evaluated the six available estimates of fault plane solutions for this earthquake and concluded that reverse slip on a 42° dipping, N112° trending fault, which extends up to the surface from a depth of 7 km, produces maximum correlation between the increased static stress and aftershock distribution. Our analysis suggests that the majority of coseismic slip occurred on the part of the rupture that lies in the depth range of 3–6.5 km.     

Intervention and assessment of earthquake knowledge at rural schools near the New Madrid seismic zone,USA

Natural Hazards - Possible recurrence of large earthquakes such as the 1811–1812 New Madrid sequence presents a significant hazard in the Central United States. The efficacy of earthquake...     

四川汶川5.12大地震同震滑动断层泥的发现及构造意义

2008年汶川8．0级地震沿龙门山断裂带内的映秀-北川断裂和灌县-安县断裂产生了近300km的同震地表破裂带。震后地质科学考察发现地表变形以逆冲为主,并伴有右旋走滑。地震地表破裂带大多沿古生代碳质泥岩、页岩和三叠系煤系地层内的滑动面出露地表,这些软弱地层为地震破裂带冲到地表提供了超低摩擦滑动带。我们发现在同震垂直和水平位错达6m左右的地袁破裂带,地震的同震滑动发生在厚度约O．5—2cm的狭窄滑动带内,以发育新鲜的灰色断层泥为特征,这些断层泥是地震断层快速滑动过程中岩石-流体相互作用的结果。     

汶川地震公路隧道普通段震害分析及震害机制研究

在汶川地震中隧道普通段出现了严重的震害,为了提高高烈度地震区隧道普通段的抗震性能,通过统计分析汶川地震公路隧道普通段震害调查资料,对隧道普通段进行了震害分析。结果表明,隧道普通段震害主要发生在软岩均一段、软硬围岩交接段软岩部分以及围岩缺陷段。通过三维有限差分数值模拟和现场典型震害分析,研究了公路隧道普通段的震害机制,探明了围岩均一普通段隧道震害机制,即硬岩隧道普通段基本无震害,软岩隧道普通段震害较严重,衬砌结构受地震惯性力影响较大,强制位移次之;探明了围岩软硬交接普通段隧道震害机制,即震害主要发生在软硬围岩交接段软岩部分,震害的主要影响因素为强制位移,地震惯性力次之;探明了围岩缺陷普通段隧道震害机制,即震害主要由地震惯性力造成强制位移影响很小。研究成果可为高烈度地震区公路隧道普通段的抗震设防提供参考。     

郯庐断裂带周边地震活动图像及未来大震危险性分析

对比1668年郯城8.5级地震前的地震围空区图像和现代华东地区地震围空区图像,发现二者很相似,都是环绕郯庐断裂带周围呈椭圆形,长轴走向近南北,面积约39×10~4km~2。涵盖河北、河南、山东、江苏、安徽、上海、黄海地区。空区范围、形态、面积、形成时间近似。另外,自2010年以来华东地区地震前兆台网出现大范围异常,集中在郯庐断裂带附近。因此,郯庐带南段的大震危险概率可能正在增大,需引起进一步关注。     

Geometry of the fault zone of the 2003 Ms = 7.5 Chuya earthquake and associated stress fields, Gorny Altai

The co-seismic deformations produced during the September 27, 2003 Chuya earthquake (Ms = 7.5) that affected the Gorny Altai, Russia, are described and discussed along a 30 km long segment. The co-seismic deformations have manifested themselves both in unconsolidated sediments as R- and R′-shears, extension fractures and contraction structures, and in bedrock as the reactivation of preexisting schistosity zones and individual fractures, as well as development of new ruptures and coarse crushing zones. It has been established that the pattern of earthquake ruptures represents a typical fault zone trending NW–SE with a width reaching 4–5 km and a dextral strike–slip kinematics. The initial stress field that produced the whole structural pattern of co-seismic deformations during the Chuya earthquake, is associated with a transcurrent regime with a NNW–SSE, almost N–S, trending of compressional stress axis (σ₁), and a ENE–WSW, almost E–W, trending of tensional stress axis (σ₃). The state of stress in the newly-formed fault zone is relatively uniform. The local stress variations are expressed in insignificant deviation of σ₁ from N–S to NW–SE or NE–SW, in short-term fluctuations of relative stress values in keeping their spatial orientations, or in a local increase of the plunge angle of the σ₁. The geometry of the fault zone associated with the Chuya earthquake has been compared with the mechanical model of fracturing in large continental fault zones with dextral strike–slip kinematics. It is apparent that the observed fracture pattern corresponds to the late disjunctive stage of faulting when the master fault is not fully developed but its segments are already clearly defined. It has been shown that fracturing in widely different rocks follows the common laws of the deformation of solid bodies, even close to the Earth surface, and with high rates of movements.     

Detailed fault structure of the Tarutung Pull-Apart Basin in Sumatra,Indonesia, derived from local earthquake data

The Tarutung Basin is located at a right step-over in the northern central segment of the dextral strike-slip Sumatran Fault System (SFS). Details of the fault structure along the Tarutung Basin are derived from the relocations of seismicity as well as from focal mechanism and structural geology. The seismicity distribution derived by a 3D inversion for hypocenter relocation is clustered according to a fault-like seismicity distribution. The seismicity is relocated with a double-difference technique (HYPODD) involving the waveform cross-correlations. We used 46,904 and 3191 arrival differences obtained from catalogue data and cross-correlation analysis, respectively. Focal mechanisms of events were analyzed by applying a grid search method (HASH code). Although there is no significant shift of the hypocenters (10.8 m in average) and centroids (167 m in average), the application of the double difference relocation sharpens the earthquake distribution. The earthquake lineation reflects the fault system, the extensional duplex fault system, and the negative flower structure within the Tarutung Basin. The focal mechanisms of events at the edge of the basin are dominantly of strike-slip type representing the dextral strike-slip Sumatran Fault System. The almost north–south striking normal fault events along extensional zones beneath the basin correlate with the maximum principal stress direction which is the direction of the Indo-Australian plate motion. The extensional zones form an en-echelon pattern indicated by the presence of strike-slip faults striking NE–SW to NW–SE events. The detailed characteristics of the fault system derived from the seismological study are also corroborated by structural geology at the surface.