Evidence for a fluid flow triggered spatio-temporal migration of seismicity in the 2001 M<Subscript>w</Subscript> 7.7 Bhuj earthquake region,Gujarat, India,during 2001–2013

We studied the variations in spatial and temporal clustering of earthquake activity (during 2001–2013) in the Kachchh seismic zone, Gujarat, India, by precisely relocating 3478 events using a joint hypocentral determination (JHD) relocation technique, and high-quality arrival times of 21032 P- and 20870 S-waves. Temporal disposition of estimated station corrections of P- and S-waves suggests that the fluid flow in the causative fault zone of the 2001 Bhuj mainshock increased during 2001–2010, while it reduced during 2011–2013, due to the healing process associated with the perturbed Kachchh fault zone. We also estimated the isotropic seismic diffusivities from epicentral growth patterns, which are found to be much lower than those observed for reservoir-induced seismicity sites in the world. Finally, we analysed the spatial and temporal evolution of this earthquake sequence by solving the diffusion equation of pore-pressure relaxation caused by co- and post-seismic stress changes associated with earthquakes. The value of the isotropic diffusivity is estimated to be 100 m²/s for the Kachchh rift zone. This gives a higher permeability (after a lapse time of 14 years from the occurrence of the 2001 Bhuj mainshock) in comparison to those observed for other intraplate regions in the world. Our results suggest that the observed spatio-temporal migration of seismicity is consistent with the shallow (meteoric water circulation at 0–10 km depths) and deeper (metamorphic fluid and volatile CO₂ circulation at 10–40 km depths) fluid flows in the permeable and fractured causative fault zone of the 2001 Bhuj earthquake.     

A media-based assessment of damage and ground motions from the january 26th, 2001<Emphasis Type="Italic">M</Emphasis> 7.6 Bhuj,India earthquake

We compiled available news and internet accounts of damage and other effects from the 26th January, 2001, Bhuj earthquake, and interpreted them to obtain modified Mercalli intensities at over 200 locations throughout the Indian subcontinent. These values are used to map the intensity distribution using a simple mathematical interpolation method. The maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated towards the western edge of the inferred fault, consistent with western directivity. Significant sedimentinduced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium such as mud flats and salt pans. In addition we use fault rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0–1000 km. We then convert the predicted hard rock ground motion parameters to MMI using a relationship (derived from internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1–2 units than those estimated from news accounts. This discrepancy is generally consistent with the expected effect of sediment response, but it could also reflect other factors such as a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. Our modeling results also suggest, however, that the Bhuj earthquake generated more high-frequency shaking than is expected for earthquakes of similar magnitude in California, and may therefore have been especially damaging.     

Post-seismic deformation associated with the 2001 Bhuj earthquake

We report the results of GPS measurements of post-seismic deformation due to the 2001 Bhuj earthquake in the Kachchh region, western India. The estimated horizontal velocity vectors in ITRF05 are in the range of 48?C49?mm/year in N46?C50°E. The observed velocity at the Gandhinagar permanent site, a far off site from the earthquake source region and probably unaffected by the post-seismic deformation, is 49?±?1?mm/year in N47°E, which is consistent with the predicted motion of Indian plate at Gandhinagar. At other sites in the source region, transient post-seismic deformation is found to be low; it attenuated rapidly within 3?C4?years of the earthquake and is much low now. Our results support the idea that mantle rheology is weak in the region.     

Postseismic relaxation due to Bhuj earthquake on January 26, 2001: possible mechanisms and processes

Earthquakes cause static stress perturbations in the nearby crust and mantle. Obeying rheological laws, this stress relaxes in a time frame of months to years with the spatial extent of few km to hundreds of km. While postseismic relaxation associated with major inter-plate earthquakes is well established, there have been few opportunities to explore its occurrence following intraplate earthquakes. The M _w 7.6 Bhuj earthquake on January 26, 2001 in western India is considered to be an intraplate event and provided a unique opportunity to examine post-earthquake relaxation processes sufficiently away from plate boundaries. To study the characteristics of transient postseismic deformation, six Global Positioning System campaigns were made at 14 sites. The postseismic transients were delineated after removing plate motions from the position time series. Postseismic deformation has been observed at all the sites in the study area. During 2001?C2007, the site closest to the epicenter exhibited postseismic deformation of about 30 and 25?mm in the north and east components, respectively. Time series of the NS and EW components of the postseismic transients can be fitted to both logarithmic and exponential functions. Close to the epicenter, the logarithmic function fits well to the initial transient, and an exponential function fits well to the later phases. The remaining sites (located east and west of the epicentral region) exhibited significantly diminished north?Csouth relaxation. Rapidly decaying afterslip and poroelastic mechanisms seem to be responsible for postseismic relaxation in the vicinity of epicenter during the initial period subsequent to the Bhuj earthquake. Postseismic relaxation by viscoelastic flow below the seismogenic zone seems to affect displacements across the entire Bhuj region. This paper presents the characteristics of postseismic transients and deformation processes in the scenario of the highly heterogeneous crust in the Bhuj region.     

Triggered/migrated seismicity due to the 2001 M w7.7 Bhuj earthquake, Western India

Paper describes triggered seismicity to 200?km distance and for a decade due to the 2001 M _w7.7 Bhuj earthquake. The Kachchh region is seismically one of the most active intraplate regions of the World due to the occurrence of two large earthquakes 1819 (M _w7.8) and 2001 (M _w7.7). Though, it has high hazard but was known to have low seismicity in view of the occurrence of fewer smaller shocks. However, the status seems to have changed after 2001. Besides the strong aftershock activity for over a decade, seismicity has spread to nearby faults in Kachchh peninsula and at several places southward for 200?km distance in Saurashtra peninsula. Beyond the rupture zone of the 2001 Bhuj earthquake, more than 40 mainshocks of M _w?~?3?C5 have occurred at 20 different locations, which is unusual. The increased seismicity is inferred to be caused by stress perturbation due to the 2001 Bhuj earthquake by viscoelastic process. In Saurashtra, over and above the viscoelastic stress increase, the transient stress increase by water table rise in monsoons seems to be affecting the timing of mainshocks and associated sequences of earthquakes.     

Slip distribution and source parameters of the 20 July 2017 Bodrum-Kos earthquake (Mw6.6) from GPS observations

Greek-Turkish boundary near the cities Kos and Bodrum has been shaken on July 20, 2017 by a Mw6.6 earthquake. The mainshock is located offshore and did not generate an on-land surface rupture. Analyzing pre- and post-earthquake continuous/survey-type static GPS observations, we investigated co-seismic surface displacements at 20 sites to characterize source parameters and slip-distribution of the mainshock. Fault plane solutions as well as co-seismic slip distribution have been acquired through the inversion of co-seismic GPS displacements modeling the event as elastic dislocations in a half space. Fault plane solution shows a southward dipping normal-type fault segment extending a depth down to ~12 km, which remains within the brittle upper crust. Results from the distributed slip inversion show that the mainshock activated a ~65 km fault section, which has three high slip patches, namely western, central and eastern patches, where the coseismic slips reach up to 13, 26, and 5 cm, respectively. This slip pattern indicates that the pre-earthquake coupling, which is storing the slip deficit, occurred on these three patches.     

Statistical parameters of Bhuj earthquake sequence of January 26th, 2001

An intraplate earthquake of magnitude (M _c) 6.9 (Anon 2001a) struck Bhuj and the adjoining region of Kachchh in Gujarat on January 26th, 2001 at about 0316 hrs (GMT) and was followed by a number of aftershocks. The epicentre of this earthquake was located at 23.4‡N and 70.28‡E close to the Kachchh mainland fault. The intensity observed around the epicenter was X on the MSK scale. A study of 531 aftershocks, in the magnitude range of 3.0–5.7, recorded at Vadodara Seismological Observatory till March 31st, 2001 has been carried out and various statistical parameters calculated. The total energy released during the study period is calculated to be 8.2 × 10¹⁴ joule. Sudden occurrence of the main shock without any foreshock in the same tectonic system is a unique feature of this sequence. Theb- value (0.86), value of M₀-M₁ (1.2), high M₁/M₀ (0.89) and high value of the decay constanth (0.91), all support the tectonic origin of the present study.     

3-D seismic structure of the Kachchh,Gujarat, and its implications for the earthquake hazard mitigation

Several pieces of studies on the January 26, 2001, Bhuj earthquake (Mw 7.6) revealed that the mainshock was triggered on the hidden unmapped fault in the western part of Indian stable continental region that caused a huge loss in the entire Kachchh rift basin of Gujarat, India. Occurrences of infrequent earthquakes of Mw 7.6 due to existence of hidden and unmapped faults on the surface have become one of the key issues for geoscientific research, which need to be addressed for evolving plausible earthquake hazard mitigation model. In this study, we have carried out a detailed autopsy of the 2001 Bhuj earthquake source zone by applying three-dimensional (3-D) local earthquake tomography (LET) method to a completely new data set consisting of 576 local earthquakes recorded between November 2006 and April 2009 by a seismic network consisting of 22 numbers of three-component broadband digital seismograph stations. In the present study, a total of 7560 arrival times of P-wave (3820) and S-wave (3740) recorded at least 4 seismograph stations were inverted to assimilate 3-D P-wave velocity (Vp), S-wave velocity (Vs), and Poisson’s ratio (σ) structures beneath the 2001 Bhuj earthquake source zone for reliable interpretation of the imaged anomalies and its bearing on earthquake hazard of the region. The source zone is located near the triple junction formed by juxtapositions of three Indian, Arabian, and Iranian tectonic plates that might have facilitated the process of brittle failure at a depth of 25 km beneath the KRB, Gujarat, which caused a gigantic loss to both property and persons of the region. There may be several hidden seismogenic faults around the epicentral zone of the 2001 Bhuj earthquake in the area, which are detectable using 3-D tomography to minimize earthquake hazard for a region. We infer that the use of detailed 3-D seismic tomography may offer potential information on hidden and unmapped faults beneath the plate interior to unravel the genesis of such big damaging earthquakes. This study may help in evolving a comprehensive earthquake risk mitigation model for regions of analogous geotectonic settings, elsewhere in the world.     

Tectonics and crustal structures related to Bhuj earthquake of January 26, 2001: based on gravity and magnetic surveys constrained from seismic and seismological studies

Gravity and magnetic data of the Kachchh basin and surrounding regions have delineated major E–W and NW–SE oriented lineaments and faults, which are even extending up to plate boundaries in the north Arabian Sea and western boundary of the Indian plate, respectively. The epicentral zone of Bhuj earthquake and its aftershocks is located over the junction of Rann of Kachchh and median uplifts viz. Kachchh mainland and Wagad uplifts, which are separated by thrust faults. Gravity data with constraints from the results of the seismic studies along a profile suggest that the basement is uplifted towards the north along thrust faults dipping 40–60° south. Similarly gravity and magnetic modeling along a profile across Wagad uplift suggest south dipping (50–60°) basement contacts separating rocks of high susceptibility and density towards the north. One of these contacts coincides with the fault plane of the Bhuj earthquake as inferred from seismological studies and its projection on the surface coincides with the E–W oriented north Wagad thrust fault. A circular gravity high in contact with the fault in northern part of the Wagad uplift along with high amplitude magnetic anomaly suggests plug type mafic intrusive in this region. Several such gravity anomalies are observed over the island belt in the Rann of Kachchh indicating their association with mafic intrusions. The contact of these intrusives with the country rock demarcates shallow crustal inhomogeneities, which provides excellent sites for the accumulation of regional stress. A regional gravity anomaly map based on the concept of isostasy presents two centers of gravity lows of −11 to −13 mGal (10⁻⁵ m/s²) representing mass deficiency in the epicentral region. Their best-fit model constrained from the receiver function analysis and seismic refraction studies suggest crustal root of 7–8 km (deep crustal inhomogeneity) under them for a standard density contrast of −400 kg/m³. It is, therefore, suggested that significant amount of stress get concentrated in this region due to (a) buoyant crustal root, (b) regional stress due to plate tectonic forces, and (c) mafic intrusives as stress concentrators and the same might be responsible for the frequent and large magnitude earthquakes in this region including the Bhuj earthquake of January 26, 2001.     

Seismogenesis of the uninterrupted occurrence of the aftershock activity in the 2001 Bhuj earthquake zone, Gujarat, India, during 2001�C2010

Seismogenesis of aftershocks occurring in the Kachchh seismic zone for more than last 10?years is investigated through modeling of fractal dimensions, b-value, seismic velocities, stress inversion, and Coulomb failure stresses, using aftershock data of the 2001 Bhuj earthquake. Three-dimensional mapping of b-values, fractal dimensions, and seismic velocities clearly delineate an area of high b-, D-, and Vp/Vs ratio values at 15?C35?km depth below the main rupture zone (MRZ) of the 2001 mainshock, which is attributed to higher material heterogeneities in the vicinity of the MRZ or deep fluid enrichment due to the release of aqueous fluid/volatile CO₂ from the eclogitisation of the olivine-rich lower crustal rocks. We notice that several aftershocks are occurred near the contacts between high (mafic brittle rocks) and low velocity regions while many of the aftershocks including the 2001 Bhuj mainshock are occurred in the zones of low velocity (low dVp, low dVs and large Vp/Vs) in the 15?C35?km depth range, which are inferred to be the fractured rock matrixes filled with aqueous fluid or volatiles containing CO₂. Further support for this model comes from the presence of hydrous eclogitic layer at sub-lithospheric depths (34?C42?km). The depth-wise stress inversions using the P- and T-axes data of the focal mechanisms reveal an increase in heterogeneity (i.e., misfit) with an almost N?CS ??₁ orientation up to 30?km depth. Then, the misfit decreases to a minimum value in the 30?C40?km depth range, where a 60^o rotation in the ??₁ orientation is also noticed that can be explained in terms of the fluid enrichment in that particular layer. The modeling of Coulomb failure stress changes (??CFS) considering three tectonic faults [i.e., NWF, GF, and Allah bund fault (ABF)] and the slip distribution of the 2001 mainshock on NWF could successfully explain the occurrences of moderate size events (during 2006?C2008) in terms of increase in positive ??CFS on GF and ABF. In a nutshell, we propose that the fluid-filled mafic intrusives are acting as stress accentuators below the Kachchh seismic zone, which generate crustal earthquakes while the uninterrupted occurrence of aftershocks is triggered by stress transfer and aqueous fluid or volatile CO₂ flow mechanisms. Further, our results on the 3-D crustal seismic velocity structure, focal mechanisms, and b-value mapping will form key inputs for understanding wave propagation and earthquake hazard-related risk associated with the Kachchh basin.     

Groundwater radon anomalies associated with earthquakes

Earthquake-related changes in groundwater radon have been detected at a sensitive observation site located right on a major active fault in Northeast Japan. A time-series analysis based on Bayesian statistics was successfully applied to remove background variations from the observed radon data, enabling us to examine the earthquake-related changes in detail.

We set a simple criterion of amplitude and duration for an anomaly observed in our radon data; we define an anomaly as a radon change that kept its level beyond 2σ (a standard deviation over the whole observation period) during a period longer than one day. We have observed 20 radon anomalies that satisfied this criterion from January 1984 to December 1988. Most of these anomalies have turned out to be related to large earthquakes that occurred in East Japan and its surrounding area; we have identified 12 post-seismic and 2-pre-seismic radon anomalies out of a total of 30 earthquakes with magnitude M 6.0 and hypocentral distance D 1000 km.

The typical pattern of the post-seismic anomalies is a radon decrease which started just after an earthquake, lasting for periods ranging from a few days to more than one week. The amplitude of the post-seismic anomalies depends on both magnitude and hypocentral distance, and can, in general, be expressed by a simple magnitude-distance relationships.

A possible pre-seismic anomaly was observed about one week before the largest earthquake that occurred in this region during the observation period (March 6, 1984; M = 7.9, D = 1000 km). Another possible pre-seismic anomaly was observed about three days before two nearby large earthquakes that occurred at almost the same place in a time interval of 53 min (February 6, 1987; M = 6.4 and M = 6.7, D = 130 km).     

1668年郯城8 1/2 级地震构造应力场分析

1668年郯城8 1/2 级地震,发震断层南起郯城窑上北到莒县土岭,全长为130 km,由5条北北东走向的活断层段组成。郯城地震断层南段沿沂沭断裂带内的F₂断裂分布,倾向南东东,倾角为30°～60°。北段紧邻F₁断裂分布,倾向不稳定,倾角较陡(多为70°以上)。南段表现为右行逆冲或逆右行的运动性质,北段则以右行走滑为主。郯城地震断层南、北两段均发育断层泥带、断层角砾带和碎裂带,南段总宽度为几米到十几米,北段总宽度为几十米到近百米,局部发育多条断层泥带。郯城地震断层的排列方式及其几何学特征表明:为老断层复活,而非新生断层。通过断层擦痕的反演同震应力场显示:北段为北东东-南西西向挤压应力场,南段为北东-南西向的挤压应力场,该地震是发生在区域性挤压应力场状态下。这种应力场空间变化可能是地震断层几何学空间变化导致的。其同震应力场与该地区现代区域应力场是一致的,这说明郯城地震并未造成震后应力场调整或震后应力场调整时间较短,未影响到现今应力场。     

Internal geometry of reactivated and non-reactivated sandblow craters related to 2001 Bhuj earthquake,India: A modern analogue for interpreting paleosandblow craters

The liquefaction attributes and crater geometry related to 2001 Bhuj earthquake has been reconstructed by trenching along large known craters formed near Umedpar in Kachchh. The study characterises the liquefied sediments in a large reactivated crater and distinguishes it from a non-reactivated crater located nearby. These characteristics can help in the interpretation of large paleocraters formed as a result of earthquake induced liquefaction.     

Why does the co-seismic slip of the 1999 Chi-Chi (Taiwan) earthquake increase progressively northwestward on the plane of rupture?

The Chi-Chi 1999 earthquake ruptured the out-of-sequence Chelungpu Thrust Fault (CTF) in the fold-and-thrust belt in Western Central Taiwan. An important feature of this rupture is that the calculated slip increases approximately linearly in the SE–NW convergence plate direction from very little at its deeper edge to a maximum near the surface. We propose here a new explanation for this co-seismic slip distribution based on the study of both stress and displacement over the long-term as well as over a seismic cycle. Over the last 0.5 My, the convergence rate in the mountain front belt is accommodated by the frontal Changhua Fault (Ch.F), the CTF and the Shuangtung Fault (Sh.F). Based on previously published balanced cross sections, we estimate that the long-term slip of the Ch.F and of the CTF accommodate 5–30% and 30–55% of the convergence rate, respectively. This long-term partitioning of the convergence rate and the modeling of inter-seismic and post-seismic displacements suggest that the peculiar linear co-seismic slip distribution is accounted for by a combination of the effect of the obliquity of the CTF to the direction of inter-seismic loading, and of increasing aseismic creep on the deeper part of the Ch.F and CTF. Many previous interpretations of this slip distribution have been done including the effects of material properties, lubrication, site effect, fault geometry and dynamic waves. The importance of these processes with respect to the effects proposed here is still unknown. Taking into account the dip angle of the CTF, asperity dynamic models have been proposed to explain the general features of co-seismic slip distribution. In particular, recent works show the importance of heterogeneous spatial distribution of stress prior to the Chi-Chi earthquake. Our analysis of seismicity shows that previous large historic earthquakes cannot explain the amplitude of this heterogeneity. Based on our approach, we rather think that the high stress in the northern part of the CTF proposed by Oglesby and Day [Oglesby, D.D., Day, S.M., 2001. Fault geometry and the dynamics of the 1999 Chi-Chi (Taiwan) earthquake. Bull. Seismol. Soc. Am. 91, 1099–1111] reflects the latitudinal variation of inter-seismic coupling due to the obliquity of the CTF.     

Study of the epicentral trends and depth sections for aftershocks of the 26th january 2001, Bhuj earthquake in Western India

The Geological Survey of India (GSI) established a twelve-station temporary microearthquake (MEQ) network to monitor the aftershocks in the epicenter area of the Bhuj earthquake (M _w7.5) of 26th January 2001. The main shock occurred in the Kutch rift basin with the epicenter to the north of Bhachao village, at an estimated depth of 25 km (IMD). About 3000 aftershocks (M _d ≥ 1.0), were recorded by the GSI network over a monitoring period of about two and half months from 29th January 2001 to 15th April 2001. About 800 aftershocks (M _d ≥ 2.0) are located in this study. The epicenters are clustered in an area 60 km × 30 km, between 23.3‡N and 23.6‡N and 70‡E and 70.6‡E. The main shock epicenter is also located within this zone. Two major aftershock trends are observed; one in the NE direction and other in the NW direction. Out of these two trends, the NE trend was more pronounced with depth. The major NE-SW trend is parallel to the Anjar-Rapar lineament. The other trend along NW-SE is parallel to the Bhachao lineament. The aftershocks at a shallower depth (<10km) are aligned only along the NW-SE direction. The depth slice at 10 km to 20 km shows both the NE-SW trend and the NW-SE trend. At greater depth (20 km–38 km) the NE-SW trend becomes more predominant. This observation suggests that the major rupture of the main shock took place at a depth level more than 20 km; it propagated along the NE-SW direction, and a conjugate rupture followed the NW-SE direction. A N-S depth section of the aftershocks shows that some aftershocks are clustered at shallower depth ≤ 10 km, but intense activity is observed at 15–38 km depth. There is almost an aseismic layer at 10–15 km depth. The activity is sparse below 38 km. The estimated depth of the main shock at 25 km is consistent with the cluster of maximum number of the aftershocks at 20–38 km. A NW-SE depth section of the aftershocks, perpendicular to the major NE-SW trend, indicates a SE dipping plane and a NE-SW depth section across the NW-SE trend shows a SW dipping plane. The epicentral map of the stronger aftershocksM ≥ 4.0 shows a prominent NE trend. Stronger aftershocks have followed the major rupture trend of the main shock. The depth section of these stronger aftershocks reveals that it occurred in the depth range of 20 to 38 km, and corroborates with a south dipping seismogenic plane.     

Estimation of static stress changes after the 2001 Bhuj earthquake: Implications towards the northward spatial migration of the seismic activity in Kachchh,Gujarat

Spatial-temporal patterns of aftershocks of the 2001 Mw7.7 Bhuj earthquake during 2001–2008 reveal a northward spatial migration of seismic activity in the Kachchh seismic zone, which could be related with the loading stresses caused by the continued occurrences of aftershocks on the north Wagad fault (NWF), the causative fault of the 2001-mainshock. Aiming at explaining the observed northward migration of activity, we modelled the Coulomb failure stress change (DCFS) produced by the 2001-mainshock, the 2006 Mw5.6 Gedi fault (GF) and the 2007 Mw4.5 Allah bund fault (ABF) events on optimally oriented plane. A strong correlation between occurrences of earthquakes and regions of increased DCFS is obtained on the associated three faults i.e. NWF, ABF and GF. Predicted DCFS on the GF increased by 0.9 MPa at 3 km depth, where the 7^th March 2006 Mw5.6 event occurred, whereas predicted DCFS on the ABF increased by 0.07 MPa at 30 km depth, where the 15^th December 2007 Mw4.5 event occurred. Focal mechanism solutions of three events on the ABF have been estimated using the iterative inversion of broadband data from 5–10 stations, which are also constrained by the first P-motion data from 8–12 stations. These focal mechanism solutions for the ABF events reveal a dominant reverse movement with a strike-slip component along a preferred northwest or northeast dipping plane (∼50–70°). Focal mechanisms of the events on all the three fault zones reveal an N-S oriented P- axis or maximum principal stress in the region, which agrees with the prevailing N-S compression over the Indian plate. It is apparent that the northward migration of the static stress changes from the NWF, resulting from the occurrence 2001 Bhuj mainshock, might have caused the occurrence of the events on the GF and ABF during 2006–08.     

Response changes of some wells in the mainland subsurface fluid monitoring network of China, due to the September 21, 1999, Ms7.6 Chi-Chi Earthquake

About 60 hydrologic changes in response to the Chi-Chi earthquake with Ms7.6 on September 21, 1999, occurred in 52 wells, including groundwater level, temperature, discharge rate, well pressure and radon, etc., in the subsurface fluid monitoring network. These response changes were mainly co-seismic, but some pre- and post-earthquake changes occurred mainly within 5 days before and after the Chi-Chi earthquake. The response changes of different wells clustering in different tectonic areas showed different features. These changes are distributed in five areas named as A, B, C, D and E. The response changes in A area with short hypo-central distance (less than 550 km) were mainly pre-earthquake changes occurring more than 5 days before the event. Those in area B (in Huanan tectonic block) and C (in Huabei tectonic block) were mainly co-seismic changes. The hypo-central distance is about 1100–1280 and 800–1160 km, respectively. These changes were high-frequency water-level oscillations induced by seismic waves and accompanied by prominent and permanent water-level jumps and drops. There are also some post-seismic changes including discharge rate and water radon and well pressure changes in area C. Those in area D in the Yanshan tectonic block were mainly co-seismic and post-seismic changes including water level, water temperature, and water radon concentration, etc., showing prominent and permanent water-level jumps and drops and rising concentrations of water radon. The hypo-central distance is about 1750–2060 km. Those in Area E were mainly co-seismic changes showing prominent and permanent water-level jump. The hypo-central distance is about 1810–2120 km. Three moderate earthquakes occurred in area D and one strong earthquake occurred in area E 4 months after the Chi-Chi earthquake. The different features of the response changes might be caused by the changes of local hydrologic conditions (like permeability) induced by seismic waves. On the other hand, these response changes might indicate the near-critical conditions in the area where the response changes clustered. Such changes might be understood by the crustal buckling hypothesis. It is thought that the response changes might be a kind of precursor that implies elevated earthquake risk in the region.     

Co-seismic Faults and Geological Hazards and Incidence of Active Fault of Wenchuan Ms 8.0 Earthquake,Sichuan,China

There are two co-seismic faults which developed when the Wenchuan earthquake happened. One occurred along the active fault zone in the central Longmen Mts.and the other in the front of Longmen Mts.The length of which is more than 270 km and about 80 km respectively.The co-seismic fault shows a reverse flexure belt with strike of N45°-60°E in the ground,which caused uplift at its northwest side and subsidence at the southeast.The fault face dips to the northwest with a dip angle ranging from 50°to 60°.The...     

Evaluation of the GPS Precise Point Positioning technique during the 21 July 2017 Kos-Bodrum (East Aegean Sea) Mw 6.6 earthquake

Precise Point Positioning (PPP) algorithms have been widely used in the Global Positioning System (GPS)-based applications. A PPP technique with a single receiver provides effective solutions where accurate absolute positioning is required. This paper provides the performance assessment of GPS PPP for detecting the displacements caused by an earthquake. For this purpose, the earthquake that occurred on 21 July 2017 at Kos-Bodrum with the impact of Mw 6.6 was investigated by analyzing the data of the permanent GPS stations located around the related region with the PPP technique. The location distances of these GPS stations range from 10 to 89 km to the epicenter of this earthquake. GPS data provided from seven permanent stations from the Continuously Operating Reference Stations-Turkey (CORS-TR) and local Bodrum CORS networks were processed to determine the co-seismic displacements during the earthquake. The data of these stations for days of year (DOYs) 200, 201, 202, and 203 were analyzed with post-process static PPP and kinematic PPP methods. GIPSY-OASIS II v6.4 was used for processing the data and all of the solutions were performed in the ITRF2008 reference frame. Two strategies were followed on the post-process static solutions. In the first strategy, 4-day data with 24-h observations were separately analyzed day by day. In the second strategy, the 24-h data were divided into 3-h duration, which is the minimum duration for optimum PPP solutions, and then the analyses were performed. When the displacements between DOYs 200 and 203 are considered in the 24-h data analysis, significant displacements have been observed through northwest direction in the northern stations whereas MUG1 is excluded. Moreover, there is significant displacement through the southeast direction in the station DATC located in the south of the epicenter. When the 3-h solutions are examined, displacements, especially on n and e directions, are observed starting from the solutions, which include Mw 6.6 earthquake. According to the kinematic PPP solutions, the effects of the Mw 6.6 earthquake can be seen clearly in the stations DATC, ORTA, TRKB, and YALI. Considering all outcomes, the PPP technique with both static and kinematic solutions provides effective results for detecting the displacements during the earthquake.     

Co-seismic displacements associated to the Molise (Southern Italy) earthquake sequence of October–November 2002 inferred from GPS measurements

The 2002 earthquake sequence of October 31 and November 1 (main shocks Mw = 5.7) struck an area of the Molise region in Southern Italy. In this paper we analyzed the co-seismic deformation related to the Molise seismic sequence, inferred from GPS data collected before and after the earthquake, that ruptured a rather deep portion of crust releasing a moderate amount of seismic energy with no surface rupture. The GPS data have been reduced using two different processing strategies and softwares (Bernese and GIPSY) to have an increased control over the result accuracy, since the expected surface displacements induced by the Molise earthquake are in the order of the GPS reliability. The surface deformations obtained from the two approaches are statistically equivalent and show a displacement field consistent with the expected deformation mechanism and with no rupture at the surface. In order to relate this observation with the seismic source, an elastic modeling of fault dislocation rupture has been performed using seismological parameters as constraints to the model input and comparing calculated surface displacements with the observed ones. The sum of the seismic moments (8.9 × 10¹⁷ Nm) of the two main events have been used as a constraint for the size and amount of slip on the model fault while its geometry has been constrained using the focal mechanisms and aftershocks locations. Since the two main shocks exhibit the same fault parameters (strike of the plane, dip and co-seismic slip), we modelled a single square fault, size of 15 km × 15 km, assumed to accommodate the whole rupture of both events of the seismic sequence. A vertical E–W trending fault (strike = 266°) has been modeled, with a horizontal slip of 120 mm. Sensitivity tests have been performed to infer the slip distribution at depth. The comparison between GPS observations and displacement vectors predicted by the dislocation model is consistent with a source fault placed between 5 and 20 km of depth with a constant pure right-lateral strike-slip in agreement with fault slip distribution analyses using seismological information. The GPS strain field obtained doesn't require a geodetic moment release larger than the one inferred from the seismological information ruling out significant post-seismic deformation or geodetic deformation released at frequencies not detectable by seismic instruments. The Molise sequence has a critical seismotectonic significance because it occurred in an area where no historical seismicity or seismogenic faults are reported. The focal location of the sequence and the strike-slip kinematics of main shocks allow to distinguish it from the shallower and extensional seismicity of the southern Apennines being more likely related to the decoupling of the southern Adriatic block from the northern one.