Reservoir Triggered Seismicity (RTS) and well water level response in the Koyna–Warna region, India

Water level fluctuations in twenty-one observation wells have been monitored for the last 10 years around the seismically active Koyna–Warna region, western India where earthquakes continue to occur even after four decades of the initiation of the seismic activity in the region. Fourteen of the observation wells act as volume strain meters as their water levels show earth tidal signals. Our analysis suggests three types of response of the well water levels to seismo-tectonic effects, i) one to local earthquakes, ii) to regional and teleseismic events, and iii) to local fluctuations in rock strain on regional scale. We observed five cases of co-seismic step-like well water level changes, of the order of few centimeters in amplitude, related to earthquakes in the magnitude range 4.3 ≤ M ≤ 5.2. All these earthquakes occurred within the network of wells drilled for the study and within 25 km distance of the recording wells. In three cases, drop in well levels preceded co-seismic step-like increases, which may be of premonitory nature. The second type of response is observed to be due to the passing of seismic waves from regional and teleseismic earthquakes like the M 7.7 Bhuj event on January 26, 2001 and the M 9.3 December 26, 2004 Sumatra earthquake. The third type is a well level anomaly of centimeter amplitude coherently occurring in several wells. The anomalies are similar in shape and last for several hours to days.From our studies we conclude that the wells in the network appear to respond to regional strain variations and transient changes due to distant earthquakes. The two factors which are important to co-seismic steps due to local earthquakes are the magnitude and epicentral distance. From the limited number of events we found that all local earthquakes exceeding M ≥ 4.3 have produced co-seismic changes. No such changes were observed for earthquakes below this magnitude threshold.     

Crustal velocity structure beneath the western Andes of Colombian using receiver-function inversion

Analysis of teleseismic records obtained in two broadband seismic stations of three components located on the Andean region of Colombia is presented in this work. The two stations are located at the Western Cordillera (WC), station BOL, and at the Central Cordillera (CC), station PBLA. The analysis of seismograms was performed by inversion of the receiver functions (RF) in order to obtain the crustal velocity structure beneath the receivers. The receiver function is a spectral ratio obtained from teleseismic earthquakes recorded by broadband seismic stations, which allows the calculation of the velocity structure beneath the receiver by removing source effects in the horizontal components of the seismic traces. Data stacking was performed in order to improve signal to noise ratio and then the data was inverted by using two optimization algorithms: a genetic algorithm (GA), and a simulated annealing algorithm (SA). The present work calculates the receiver functions using teleseismic earthquakes at epicentral distances (Δ) ranging between 30° and 90° and recorded at the two stations within the years 2007 and 2009.Delay times between P and PS waves converted at the Moho boundary were used to constrain the velocity structure. The receiver functions at the stations were generated from seismic events within a broad range of back azimuth. Data from gravity and magnetism were also used during the geophysical survey. The depth of the Moho boundary was found to be at 40 km in the WC beneath station BOL and at 43 km in the CC beneath station PBLA. The upper crust, with a thickness of 5 km, is characterized by a shear wave velocity of about 3.0 km s⁻¹; the shallower layers, at approximately 1.0 km, have shear wave velocities between 2.2 and 2.6 km s⁻¹, which corresponds to sediments overlying the upper crust. These observations support the hypothesis of a thickness of the crust at the root of the mountain range to be between 32 and 50 km. The calculated receiver functions were compared with artificial ones generated from the inversion of 48000 models of horizontal layers for each station using a GA and an SA that allowed a satisfactory coverage of all the sample space in order to avoid non-unique solutions. Beneath station BOL a moderate low-velocity zone (LVZ) was found, which was caused by accretionary processes of the ophiolite complex in the WC.     

Role of African–Eurasian plate setting in the felt areas of intermediate‐depth earthquakes: an investigation using crowdsourced data

Greek intermediate‐depth earthquakes, occurring in the subducted plate of the Hellenic Arc, are felt at greater distances than expected, reaching Italy in some cases. We study in detail macroseismic intensity data from intermediate‐depth Italian and Greek earthquakes collected from Internet users who felt the shaking in Italy. The huge amount of data allowed us to outline the felt/not‐felt limit and to find a correspondence between attenuation areas and the presence of asthenospheric material at shallow depths. We show that plate boundaries, known to produce the majority of earthquakes, are, in some specific cases, the boundaries of areas in which earthquakes are felt. The Ionian subducted lithosphere propagates seismic waves with low attenuation over large distances, whereas high‐attenuation zones in Italy, linked to asthenospheric upwelling, limit the propagation, as evidenced also by PGA values. We identify a typical pattern that can be used to recognize intermediate‐depth earthquakes, and to properly locate historical events.     

Uniform slip rates of the Altyn Tagh and the Kunlun faults likely reflect lateral variation of frictional strength of the faults

The slip rate predicted from geodetic and geomorphological measurements is quite uniform on ~800–1,000 km length of the Altyn Tagh and the Kunlun faults. GPS velocity field documents that tectonic loading on the two active faults changes greatly along their strikes. To explore the mechanical relationship between far‐field tectonic loading and fault slip‐rate accumulation, we built a 3D viscoelastic finite‐element model with fault motion governed by frictional strength of contact interfaces. Based on numerical experiments, it is found that the observed uniform slip rate could reflect lateral variation of frictional strength along fault strike. Variation of predicted frictional coefficient ranges from ~0.13 to ~0.02, verifying that the two active faults must be weak for their mechanical strength. In addition, the good fitness between the relatively weak segment of faults and the location of strong earthquakes suggests that seismic activity along the two faults could be related to their frictional strength.     

Effects of subducting lithosphere on p-wave amplitudes and residuals

In studying the structure of the lithosphere and asthenosphere using the kinematic and dynamic parameters of seismic waves, in any area, the tectonics of focal regions should also be considered. Subduction zones represent large-scale inhomogeneities which affect the propagation of seismic waves both at small and teleseismic epicentral distances.A study on the magnitude corrections of seismic stations in Central Europe and Scandinavia revealed that even in the case of close seismic stations the observed differences of the amplitudes of teleseismic P waves depend on the strike and dip of the sinking plates in the Northwestern Pacific. The smaller the angle between a seismic ray and the bottom of the lithospheric plate, the larger the numerical decrease of the magnitude corrections. A similar dependence was found for the P-wave residuals; the rays propagating more along the slab are accelerated more. Besides that, we have observed systematic differences in the P-wave residuals of about 2 s for foci with different positions within the subducting lithosphere of the Kurile arc in investigating the dependence of the residuals on the epicentral distance in the direction of subduction. The differences in the residuals disappear in this region at an epicentral distance of about 52°.     

The Tornquist Zone, a north east inclining lithospheric transition at the south western margin of the Baltic Shield: Revealed through a nonlinear teleseismic tomographic inversion

From July 1996 to August 1997 the TOR project operated 130 seismographs in North Germany, Denmark and South Sweden, with the aim of collecting signals from local, regional and teleseismic earthquakes. This data set is particularly interesting since the seismic antenna crosses the most significant geological boundary in Europe, the Tornquist Zone, which in the northern part is the border between the Baltic Shield and the younger European lithosphere. Previous studies have shown significant physical changes in the crust and upper mantle across this transition zone, including two independent teleseismic tomographic studies of the TOR data set. But these two studies disagree on the orientation of the slope of the transition. Both studies used an iterative linearized inversion method. We will in this work Preprint submitted to Elsevier Science 27 July 2005 present an inversion based on Bayesian statistics, where the solution space is examined in order to study a very large number of tomographic solutions and to examine the solution uniqueness and uncertainty. The method is applied to measurements of 3345 relative teleseismic P-phase travel times from 48 teleseismic earthquakes with good azimuthal coverage with respect to the great circle arc of the TOR array. We find the lithospheric transition to be a north east inclination of around 30° to 45° off vertical.     

Aseismic deep subduction of the Philippine Sea plate and slab window

We have made great efforts to collect and combine a large number of high-quality data from local earthquakes and teleseismic events recorded by the dense seismic networks in both South Korea and West Japan. This is the first time that a large number of Korean and Japanese seismic data sets are analyzed jointly. As a result, a high-resolution 3-D P-wave velocity model down to 700-km depth is determined, which clearly shows that the Philippine Sea (PHS) plate has subducted aseismically down to ∼460 km depth under the Japan Sea, Tsushima Strait and East China Sea. The aseismic PHS slab is visible in two areas: one is under the Japan Sea off western Honshu, and the other is under East China Sea off western Kyushu. However, the aseismic PHS slab is not visible between the two areas, where a slab window has formed. The slab window is located beneath the center of the present study region where many teleseismic rays crisscross. Detailed synthetic tests were conducted, which indicate that both the aseismic PHS slab and the slab window are robust features. Using the teleseismic data recorded by the Japanese stations alone, the aseismic PHS slab and the slab window were also revealed (Zhao et al., 2012), though the ray paths in the Japanese data set crisscross less well offshore. The slab window may be caused by the subducted Kyushu-Palau Ridge and Kinan Seamount Chain where the PHS slab may be segmented. Hot mantle upwelling is revealed in the big mantle wedge above the Pacific slab under the present study region, which may have facilitated the formation of the PHS slab window. These novel findings may shed new light on the subduction history of the PHS plate and the dynamic evolution of the Japan subduction zone.     

Short-period surface waves across the Iberian Peninsula area

A large number of records of near earthquakes, for the period 1950–1975, have been revised at Toledo Observatory (I.G.C.), looking for clear trains of surface waves.All records used correspond to the Wiechert, WWSS or HGLP seismometers of this station.Multiple filtering techniques have been applied to the selected seismograms to obtain Love and Rayleigh fundamental mode dispersion curves.As most of the paths are of mixed continental and oceanic structure, the first interpretation has been made considering an average model for each set of similar seismic paths. However, for some regions, the separation of the two structures (pure continental and pure oceanic) has been attempted.     

基于时频域极化滤波的勒夫型槽波提取

矿井槽波地震勘探数据中,含有多种类型地震波,包括线性极化的纵波、横波、勒夫型槽波和椭圆极化的瑞利型槽波。由于槽波频散现象严重,不同类型槽波互相重叠,在地震记录上难以区分,合理地进行波场分离是槽波数据处理的重要环节。根据不同类型槽波振动极化特征,提出基于S变换的时频域自适应协方差矩阵极化滤波方法;在时频域中,针对槽波水平双分量信号建立自适应协方差矩阵;利用矩阵对应特征值和特征向量,计算椭圆率和方位角参数来描述槽波极化振动特征;依据极化参数建立极化滤波函数,对槽波地震信号进行波场分离。经过理论合成信号和实际数据测试表明,该方法可以有效地从槽波地震记录中提取出勒夫型槽波,避免其它类型地震波在槽波数据处理过程中产生的干扰,为槽波地震勘探数据处理提供新的技术手段。      

Particle motions in a physical model of shallow angle thrust faulting

A physical model of thrust faulting has been constructed out of elastic foam rubber. The thrust wedge has an angle of 25°., and is forced from the back by a steel-hydraulic system of effectively infinite rigidity, analogous to the force exerted by tectonic plates. The observed particle motions show many features different from those commonly assumed in dislocation models of subduction thrusts. Interface waves associated with fault opening propagate along the thrust plane (from the back) and temporarily decouple the overlying hanging wall plate from the foot wall. Because of the geometry of the reflecting fault boundary, and the free surface (also reflecting), energy is trapped in the wedge and, as a consequence, the particle motions and energy in the hanging wall are much greater than in the foot wall. The interface wave and the energy trapped in the upper plate propagate up the wedge and break out at the toe of the thrust with a spectacular increase in motions. If this model is analogous to the real earth, it suggests that some common assumptions in dislocation modeling may be incorrect. The model may explain apparent missing energy radiated to teleseismic distances, the anomalously low number of empirical Green’s functions needed to model teleseismicp waves, and the evidence of intense shaking on the hanging wall toe of some thrust fault earthquakes. The results support the suggestion that interface waves may explain the lack of high frictional heat generation from thrust faults, and may explain the ‘paradox of large overthrusts’. The results also suggest that the seismic hazard of great subduction thrust earthquakes and continental shallow angle thrust faults might in some cases be seriously underestimated.     

Late-Quaternary Slip Rate and Seismic Activity of the Xianshuihe Fault Zone in Southwest China

The Xianshuihe fault zone is a seismo-genetic fault zone of left-lateral slip in Southwest China. Since 1725, a total of 59 Ms ≥ 5.0 earthquakes have occurred along this fault zone, including 18 Ms 6.0–6.9 and eight Ms ≥ 7.0 earthquakes. The seismic risk of the Xianshuihe fault zone is a large and realistic threat to the western Sichuan economic corridor. Based on previous studies, we carried out field geological survey and remote sensing interpretation in the fault zone. In addition, geophysical surveys, trenching and age-dating were conducted in the key parts to better understand the geometry, spatial distribution and activity of the fault zone. We infer to divide the fault zone into two parts: the northwest part and the southeast part, with total eight segments. Their Late Quaternary slip rates vary in a range of 11.5 mm/a –(3±1) mm/a. The seismic activities of the Xianshuihe fault zone are frequent and strong, periodical, and reoccurred. Combining the spatial and temporal distribution of the historical earthquakes, the seismic hazard of the Xianshuihe fault zone has been predicted by using the relationship between magnitude and frequency of earthquakes caused by different fault segments. The prediction results show that the segment between Daofu and Qianning has a possibility of Ms ≥ 7.0 earthquakes, while the segment between Shimian and Luding is likely to have earthquakes of about Ms 7.0. It is suggested to establish a GPS or In SAR-based real-time monitoring network of surface displacement to cover the Xianshuihe fault zone, and an early warning system of earthquakes and post seismic geohazards to cover the major residential areas.     

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.     

Application of GPS in identifying active fault plane in Western Maharashtra Peninsular shield of India

A broad overview of how GPS can be used to identify the active fault plane along which earthquakes occur in the seismic prone Western Maharashtra Peninsular shield of India is given. Indian Institute of Technology, Bombay in collaboration with Indian Institute of Geomagnetism, Mumbai was actively working in the field of crustal deformation monitoring along Western Maharashtra from year 2004–2006. A regional GPS network had been established for this purpose, periodically observed between May 2004 to May 2006, and analysed. To understand the pattern of seismicity along Western Maharashtra region, the GPS-estimated coordinates and displacements were used as an indicator to estimate and identify regions of high strain rates, the reliability of which was confirmed with real seismic data collected. These strain rates were used as a means to identify the fault plane along which the earthquakes occurred during the period of observation.     

Estimation of the rupture velocity and fault length of the 2004 Sumatra–Andaman earthquake using a dense broadband seismic array in Taiwan

The rupture process of the disastrous Sumatra–Andaman earthquake of 26 December 2004 was analyzed by array processes for teleseismic P-waves recorded by a dense broadband seismic array in Taiwan with epicentral distances of close to 31°. The azimuthal variation from the BATS array center to both ends of the rupture fault is approximately 21°, which is larger than that reported previously for seismic arrays used to image the rupture process of this earthquake, thereby providing a high spatial resolution in studying the source rupture behavior. Two array-processing methods were used to analyze teleseismic P-wave trains. Both analyses were based on data recorded by a broadband network, covering a region of 200 × 400 km, with the aim of evaluating the rupture behavior of the earthquake. Consistent results from both analyses indicate that the earthquake had a rupture duration exceeding 500 s, with major asperities encountered at 80, 260, and 330 s after the initiation of rupturing. We traced the ruptured fault for more than 1200 km from the point of initial rupture. The average rupture velocity was approximately 3.0 km/s and the major northward rupture propagation began at 80 s after the initiation of rupturing.     

Seismic interpretation of crustal structure in the flinders‐mt lofty ranges and gulf regions,South Australia

There appears to be little correlation of earthquake epicentres with known surface geological features in South Australia. Seismic wave travel‐time residuals are used to derive corrections for the velocity and depth parameters for the simple uniform crustal model which approximates to that in South Australia. Local studies of Moho depth in the seismic zone and analysis of travel‐time station corrections from both local earthquake and teleseismic data suggest that lateral and vertical variations in the South Australian crust are small. Data presented in this paper appear to be consistent with a plate tectonic model derived from focal mechanism studies (Stewart & Mount, 1972) for the active South Australian seismic zones.     

The Seismic Project of the National Tsunami Hazard Mitigation Program

In 1997, the Federal Emergency Management Agency (FEMA), National Oceanic and Atmospheric Administration (NOAA), U.S. Geological Survey (USGS), and the five western States of Alaska, California, Hawaii, Oregon, and Washington joined in a partnership called the National Tsunami Hazard Mitigation Program (NTHMP) to enhance the quality and quantity of seismic data provided to the NOAA tsunami warning centers in Alaska and Hawaii. The NTHMP funded a seismic project that now provides the warning centers with real-time seismic data over dedicated communication links and the Internet from regional seismic networks monitoring earthquakes in the five western states, the U.S. National Seismic Network in Colorado, and from domestic and global seismic stations operated by other agencies. The goal of the project is to reduce the time needed to issue a tsunami warning by providing the warning centers with high-dynamic range, broadband waveforms in near real time. An additional goal is to reduce the likelihood of issuing false tsunami warnings by rapidly providing to the warning centers parametric information on earthquakes that could indicate their tsunamigenic potential, such as hypocenters, magnitudes, moment tensors, and shake distribution maps. New or upgraded field instrumentation was installed over a 5-year period at 53 seismic stations in the five western states. Data from these instruments has been integrated into the seismic network utilizing Earthworm software. This network has significantly reduced the time needed to respond to teleseismic and regional earthquakes. Notably, the West Coast/Alaska Tsunami Warning Center responded to the 28 February 2001 Mw 6.8 Nisqually earthquake beneath Olympia, Washington within 2 minutes compared to an average response time of over 10 minutes for the previous 18 years.     

Focal mechanism solutions and dynamic parameters of earthquakes in the eastern Himalayas and northern Burma

Focal-mechanism solutions of four earthquakes in the eastern Himalayas and northern Burma are determined using the first motion of compressional waves. Two possible solutions thus obtained for each event reveal steeply dipping fault with predominantly strike-slip motion. The stress directions inferred from the focal mechanism solutions are interpreted in the light of predictions of the plate tectonics theory, viz., the underthrusting of the Indian plate in the Burma region in an easterly direction.Dynamic parameters (seismic moment, apparent stress and average dislocations) are obtained using the corrected spectra of Love waves. The earthquakes are found to possess low seismic moment and apparent stress values. A comparison of these estimates with values for intraplate earthquakes is given. It is suggested that these earthquake might be a consequence of a nonhomogeneous rupture process.     

Measuring the north–south coseismic displacement component with high‐resolution multi‐aperture InSAR

Earthquake source inversions based on space‐borne Synthetic Aperture Radar interferometry (InSAR) are used extensively. Typically, however, only the line‐of‐sight (LoS) surface displacement component is measured, which is mainly sensitive to the vertical and E–W deformations, although well‐established methods also exist to estimate the flight‐path component, which is highly sensitive to the N–S displacement. With high‐resolution sensors, these techniques are particularly appealing, because accuracies in the order of 3 cm can be achieved, while retaining spatial resolutions between 45 m and a few km, depending on the required level of filtering. We discuss the application to COSMO‐SkyMed SAR imagery of the Spectral Diversity or Multi Aperture Interferometry technique, presenting the first SAR flight‐path displacement field associated with the Mw 6.3, 2009 L'Aquila event (central Apennines). Finally, we observe and characterize a previously unknown misregistration pattern.     

The seismogenic thickness in Italy: constraints on potential magnitude and seismic hazard

The thickness of the seismogenic layer is a key parameter for seismic hazard, since it can be used to constrain the maximum depth of faulting and the potential magnitude. In this study, we compute the seismogenic thickness in the Italian region by defining the lower seismicity cut‐off, using high‐quality hypocentral locations of earthquakes that occurred in the past decade. Along the eastern Alps, the seismogenic thickness is about 12–14 km, laterally homogeneous along the entire south‐verging thrust front. In the Apennines extensional belt, lateral changes in seismogenic thickness are evident, and correlate with changes in the seismic energy released by past earthquakes. The potential magnitude is larger in the southern Apennines where the seismogenic thickness is greater (16–18 km) than in the northern Apennines where it is less (6–10 km) and seismic energy is partially released by the creeping of faults.     

The March 11, 2002 Masafi, United Arab Emirates earthquake: Insights into the seismotectonics of the northern Oman Mountains

A moderate (M 5) earthquake struck the northeastern United Arab Emirates (UAE) and northern Oman on March 11, 2002. The event was felt over a wide area of the northern Emirates and was accompanied by smaller (felt) events before and after the March 11 main shock. The event was large enough to be detected and located by global networks at teleseismic distances. We estimated focal mechanism and depth from broadband complete regional waveform modeling. We report a normal mechanism with a slight right-lateral strike-slip component consistent with the large-scale tectonics. The normal component suggests relaxation of obducted crust of the Semail ophiolite (specifically, the Khor Fakkan Block) while the right-lateral strike-slip component of the mechanism is consistent with shear across the Oman Line. Felt earthquakes are rare in the region, however no regional seismic network exists in the UAE to determine local seismicity. This event offers a unique opportunity to study the active tectonics of the region as well as inform future studies of seismic hazard in the UAE and northern Oman.