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.     

Long-term hydrochemical earthquake precursor studies at the Koyna-Warna reservoir site in western India

Koyna-Warna region of western India is an active seismic zone due to the Reservoir Triggered Seismicity (RTS). Earthquake precursor studies are carried out monitoring hydrochemical and stable isotope signatures in the groundwater from 15 bore wells since January 2005, for more than 12 years (January 2005 to February 2017). Depth of these boreholes ranges from 100 to 250 m. Cyclic or temporal variation in hydrochemistry is observed in few sensitive wells in Koyna region. The Govare well in Koyna is found to be most sensitive and the observed hydrochemical cycle is closely associated with local earthquakes of M > 5. The earthquakes M <5 occurring either in Warna cluster or close to the observation wells, did not generate hydrochemical precursory changes. The increase in hydrochemistry is hypothesized as mixing of two aquifer waters with different hydrochemistry. It is noted that a precursory hydrochemical cycle is observed during first quarter of 2015, but no earthquake M > 5.0 occurred till date. The cyclic changes in hydrochemistry, however, indicate on-going earthquake processes and an impending earthquake of M > 5 in the region.     

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.     

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.     

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.     

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.     

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:

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.     

蠕滑断裂带岩石组成和构造特征分析:以龙门山灌县-安县断裂带为例

断裂蠕滑可以连续释放部分构造应力,但仍可能造成重大的地质灾害,甚至具有发生大地震的可能性。断层岩是断裂作用中的直接产物,其物质组成和内部构造可为揭示断裂带滑移机制提供关键信息。2008年M_w 7.9汶川地震中破裂的龙门山灌县-安县断裂带具有蠕滑性质,是探究大陆内部蠕滑断裂滑移机制的最佳案例。本文以龙门山灌县-安县断裂带地表探槽和深部钻孔的断层岩为研究对象,通过碎屑统计、X射线粉末衍射矿物分析、光学显微镜和扫描电镜观测,结果显示该断裂带断层泥碎屑含量和颗粒大小均小于断层角砾岩,其粘土矿物含量高达50%以上,且断层岩中普遍发育粘土-碎屑组构以及拖尾构造、似S-C组构等多种压溶构造。综合分析发现压溶作用、低摩擦系数物质以及颗粒滑移对灌县-安县断裂带的蠕滑变形都发挥着重要作用,并且三者相辅相成,因此认为灌县-安县断裂带的蠕滑过程主要是压溶作用和摩擦-颗粒滑移机制共同作用,该认识可更好地了解地震周期并为区域防震减灾提供科学依据。

     

龙门山汶川-茂县断裂带的岩石磁学特征及其地震作用的指示意义

汶川-茂县断裂带是龙门山后山断裂,是松潘-甘孜褶皱带与龙门山断裂带之间的边界断层,然而,在2008年M_W 7.9级汶川大地震中并没有发生破裂。同时,汶川-茂县断裂带在汶川地震之前是否发生过历史大地震仍缺少确切的证据。这不仅制约着汶川地震发生机制的认识,而且还影响对龙门山形成演化过程的理解。因此,确定汶川-茂县断裂带的断裂作用环境对于认识龙门山断裂带的地震发生机制至关重要。断裂岩的岩石磁学可以有效地揭示断裂带的物理和化学环境。本文以汶川-茂县断裂带北部地表露头的断裂岩为研究对象,通过岩石磁学研究,并结合显微结构观察和地球化学分析,探讨汶川-茂县断裂带的断裂作用环境。断层泥和断层角砾岩最大磁化率值分别约为围岩的30倍和15倍,具有高磁化率值特征。断层泥的主要载磁矿物为磁铁矿、磁黄铁矿、针铁矿;断层角砾岩的主要载磁矿物为磁黄铁矿、针铁矿。断层泥的高磁化率值异常的主要原因是围岩中的含铁矿物在地震摩擦热和流体作用下新生成磁铁矿和磁黄铁矿。断层角砾岩的高磁化率值异常是围岩在含有大量硫元素的低温热液流体作用下生成了磁黄铁矿。大量针铁矿指示了震后期还原性的低温热液流体作用。断裂岩的高磁化率值异常和新生磁铁矿指示了汶川-茂县断裂带曾经发生了摩擦热温度>500℃的大地震活动,发震和震后期均为还原性的硫元素含量较高的低温热液流体环境。

     

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 的狭窄滑动带内,以发育新鲜的灰色断层泥为特征,这些断层泥是地震断层快速滑动过程中岩石—流体相互作用的结果。     

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.     

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

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

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...     

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.     

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.