A low magnitude seismic sequence near Isernia (Molise,Central Italy) in January 1986

Following the increase in seismic activity which occurred near Isernia (Molise, Central Italy) in January 1986, a digital seismic network of four stations with three-component, short-period seismometers, was installed in the area by the Osservatorio Vesuviano. The temporary network had an average station spacing of about 8–10 km and, in combination with permanent local seismic stations, allowed the accurate determination of earthquake locations during an operating period of about one month. Among the 1517 detected earthquakes, 170 events were selected with standard errors on epicentre and depth not greater than respectively 0.5 and 1.5 km. The most frequent focal depths ranged between 4 and 8 km, while the epicentres distribution covered a small area NE of Isernia of about 10 km². A main rupture zone could not be clearly identified from the spatial distribution of the earthquakes, suggesting a rupture mechanism involving heterogeneous materials. The activity was characterized by low energy levels, the largest earthquake, on January 18, 1986, havingM _D =4.0. The time sequence of events and pattern of seismic energy release revealed a strong temporal clustering of events, similar to the behaviour commonly associated with seismic swarms.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 John Wiley & Sons, Ltd.     

Refined Locations of The Swarm Earthquakes In The Nový Kostel Focal Zone and Spatial Distribution of The January 1997 Swarm In Western Bohemia,Czech Republic

In order to improve the accuracy of the spatial distribution of earthquake foci in the principal Novy Kostel focal zone, refined focal locations of about 1500 micro-earthquakes of the 1991 – 1997 period were determined using the relative Master-Event location method. To estimate the reduction in the scatter of located hypocentres, the results were compared with those obtained by routinely used FASTHYPO method and cluster analysis (the nearest-neighbour method) was applied to the located foci to evaluate the spatial distribution of the foci. Based on the results of refined location and of the cluster analysis, a concept of seismic energy release in space and time in the main focal zone was developed. Especially the January 1997 earthquake swarm was studied in detail: 946 events were located with the Master-Event location method, and the dimensions and geometry of focal clusters were determined. Type analysis was applied to waveforms to divide approximately 800 located events into eight multiplet groups to each of which a characteristic source mechanism was assigned. The spatial distribution of the foci as well as of the eight types of source mechanisms was revealed in this way and also the planes fitting the clusters of foci with two predominant source mechanisms were determined fairly well.     

Seismicity characteristics of a potentially active Quaternary volcano: The Tatun Volcano Group, northern Taiwan

The Tatun Volcano Group (TVG) is located at the northern tip of Taiwan, near the capital Taipei and close to two nuclear power plants. Because of lack of any activity in historical times it has been classified as an extinct volcano, even though more recent studies suggest that TVG might have been active during the last 20 ka. In May 2003 a seismic monitoring project at the TVG area was initiated by deploying eight three-component seismic stations some of them equipped with both short-period and broadband sensors. During the 18 months observation period local seismicity mainly consisted of high frequency earthquakes either occurring as isolated events, or as a continuous sequence in the form of spasmodic bursts. Mixed and low frequency events were also present during the same period, even though they occurred only rarely. Arrival times from events with clear P-/S-wave phases were inverted in order to obtain a minimum 1D velocity model with station corrections. Probabilistic nonlinear earthquake locations were calculated for all these events using the newly derived velocity model. Most high frequency seismicity appeared to be concentrated near the areas of hydrothermal activity, forming tight clusters at depths shallower than 4 km. Relative locations, calculated using the double-difference method and utilising catalogue and cross-correlation differential traveltimes, showed insignificant differences when compared to the nonlinear probabilistic locations. In general, seismicity in the TVG area seems to be primarily driven by circulation of hydrothermal fluids as indicated by the occurrence of spasmodic bursts, mixed/low frequency events and a b-value (1.17 ± 0.1) higher than in any other part of Taiwan. These observations, that are similar to those reported in other dormant Quaternary volcanoes, indicate that a magma chamber may still exist beneath TVG and that a future eruption or period of unrest should not be considered unlikely.     

Calibration of the Regional Crustal Waveguide and the Retrieval of Source Parameters Using Waveform Modeling

v--vRegional crustal waveguide calibration is essential to the retrieval of source parameters and the location of smaller (M < 4.8) seismic events. This path calibration of regional seismic phases is strongly dependent on the accuracy of hypocentral locations of calibration (or master) events. This information can be difficult to obtain, especially for smaller events. Generally, explosion or quarry blast generated travel-time data with known locations and origin times are useful for developing the path calibration parameters, but in many regions such data sets are scanty or do not exist. We present a method which is useful for regional path calibration independent of such data, i.e. with earthquakes, which is applicable for events down to M_w = 4 and which has successfully been applied in India, central Asia, western Mediterranean, North Africa, Tibet and the former Soviet Union. These studies suggest that reliably determining depth is essential to establishing accurate epicentral location and origin time for events. We find that the error in source depth does not necessarily trade-off only with the origin time for events with poor azimuthal coverage, but with the horizontal location as well, thus resulting in poor epicentral locations. For example, hypocenters for some events in central Asia were found to move from their fixed-depth locations by about 20 km. Such errors in location and depth will propagate into path calibration parameters, particularly with respect to travel times. The modeling of teleseismic depth phases (pP, sP) yields accurate depths for earthquakes down to magnitude M_w = 4.7. This M_w threshold can be lowered to four if regional seismograms are used in conjunction with a calibrated velocity structure model to determine depth, with the relative amplitude of the P_nl waves to the surface waves and the interaction of regional sPmP and pPmP phases being good indicators of event depths. We also found that for deep events a seismic phase which follows an S-wave path to the surface and becomes critical, developing a head wave by S to P conversion is also indicative of depth. The detailed characteristic of this phase is controlled by the crustal waveguide. The key to calibrating regionalized crustal velocity structure is to determine depths for a set of master events by applying the above methods and then by modeling characteristic features that are recorded on the regional waveforms. The regionalization scheme can also incorporate mixed-path crustal waveguide models for cases in which seismic waves traverse two or more distinctly different crustal structures. We also demonstrate that once depths are established, we need only two-stations travel-time data to obtain reliable epicentral locations using a new adaptive grid-search technique which yields locations similar to those determined using travel-time data from local seismic networks with better azimuthal coverage.     

廊固凹陷深部剪切破裂构造的地震学证据

基于区域地震台网观测数据,采用近震波形反演方法,确定2018年2月12日河北永清M4.3地震的最佳双力偶源震源机制解为:节面Ⅰ走向297°,倾角58°,滑动角-32°;节面Ⅱ走向45°,倾角63°,滑动角-144°;是一个略带正断分量的右旋走滑地震.结合近震转换波测定主震的震源深度在19km附近.地震序列的双差定位结果显示:永清地震序列震中呈北东向窄带展布,表明此次地震主要向北东向破裂;深度集中分布在17～19km,整体形态近于铅直,显示发震断裂具有走向北东、倾向南东、倾角陡立的特征,与节面Ⅱ的性质比较吻合,推测节面Ⅱ为发震断层面.将发震断层面参数与震源区附近断裂性质进行对比分析,形成了关于廊固凹陷附近区域地震构造的一些认识:(1)推测永清地震的发震构造不是地壳浅部发育的先存正断裂,而是震源区下方一条地壳尺度的深断裂,该深断裂为新生断裂,具有右旋走滑正断性质,倾角陡峭、近于直立、宽度较大,向上与夏垫断裂相通.(2)综合震源区附近多条深地震反射剖面探测结果,推测永清地震的发震断裂与新夏垫断裂同属一条断裂,称为:新夏垫深断裂.该断裂从夏垫向西南方向延伸至文安,并可能与霸县—束鹿—邯郸断裂带相联系,总长度超过150km.(3)基于2006年文安M5.1地震与2018年永清M4.3地震在震源机制上的相似性及震源位置上的关联性,结合区域构造条件,认为两次地震的发震构造均为新夏垫深断裂.(4)根据研究区几次显著地震的震源深度分布特征,参考区域断层构造、电性结构和流变学模型,推测活化克拉通块体新生断裂的脆韧性转换界面深度在15km附近.     

天山中段地区重复地震时空特征研究

基于2009—2017年新疆区域数字地震台网记录的地震波形数据,利用波形互相关技术及主事件定位方法识别并重新定位了新疆天山中段及其周缘的重复地震。以波形互相关系数0.9作为阈值来确定研究区的重复地震事件,统计结果显示3万零181个事件中的1万1 618个为重复地震事件,这些重复地震事件组成了2395组重复地震对和重复地震丛,占总事件数的38.5%。根据重复地震重定位前后地震对之间距离的统计结果推测,该区域的系统定位误差约为5—10 km。进一步结合该区域最新的震源分类结果对不同震源类型重复地震的时空分布特征予以分析,结果显示:重复矿山爆破事件在空间上呈丛集性分布,且其中的93.6%发生于白天,同时呈现季节性发生模式,即爆破多发生于夏季,而冬季较少;而重复构造地震在空间上大多沿断层分布,24小时内呈随机分布的特征,且研究时段内每个月的活动水平相对平稳;重复诱发地震成丛分布于靠近油气田和水库的区域,但其中部分诱发地震的位置与构造地震重叠,发震时间特征与构造地震相似,为随机分布。      

Detection of earthquake clusters on the basis of waveform similarity: An application in the monferrato region (Piedmont, Italy)

In August 2000 and July 2001 two seismic sequences, characterized by mainshocks with Ml (local magnitude) respectively 5.1 and 4.8, occurred in the Monferrato region (Italy). The regional seismic network of North-western Italy (RSNI) recorded more than 250 foreshocks and aftershocks. The routine locations, obtained from the Hypoellipse code, show a seismic activity concentrated in a circular area, of a radius of about 15 km, located near Acqui Terme, and randomly distributed in depth. Location errors, due to an uneven azimuthal station distribution of the regional seismic network, prevent recognition of the geometry of the active zone. Waveform analysis revealed the presence of several multiplets. In order to discriminate and successively relocate them, an automatic cross-correlation procedure was applied. Normalized cross-correlation matrixes, for the RSNI stations which recorded almost 90% of considered events, on the basis of a signal to noise ratio analysis, were computed using only S wave time windows. The use of a relocation procedure, based on the double-difference method which incorporates ordinary absolute travel-time measurements and/or cross-correlation differential travel-times, allowed us to overcome the constraints of the uneven distribution of stations giving a quite different frame of seismicity. The improved locations showed that the seismic activity is mainly arranged along a NE-SW oriented volume, at a depth range of 8–20 km, involving the basement crystalline units. This orientation is confirmed by the analysis of the focal mechanisms: most focal solutions show a strike slip component with one of the nodal planes consistent with the main orientation of the seismic events.     

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java,Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults. In this study, double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java. To infer the geometry of the structure beneath the study area, precise earthquake hypocenter determination was first performed before tomographic imaging. For this, earthquake waveform data were extracted from the regional Meteorological, Climatological, Geophysical Agency (BMKG) network of Indonesia from South Sumatra to Central Java. The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected, the key features being events of magnitude > 3, azimuthal gap < 210° and number of phases > 8. A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters. The hypocenter locations were then relocated using double-difference tomography (tomoDD). A significant reduction of travel-time (root mean square basis) and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java. Double-difference tomography was found to give a clear velocity structure, especially beneath the volcanic arc area, i.e., under Mt Anak Krakatau, Mt Salak and the mountains complex in the southern part of West Java. Low velocity anomalies for the P and S waves as well as the v_P/v_S ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.     

Source modeling of the Hsingtai,China earthquakes of March 1966

The Hsingtai, China earthquakes of March 1966 were a series of destructive earthquakes associated with the Shu-lu graben. Five strong shocks of M_s ≥ 6 occurred within a period of less than a month, the largest of which was M_s 7.2. Body and surface waves over the period range from several to 100 s have been modeled for the four largest events using synthetic seismograms in the time domain and spectral analysis in the frequency domain. Data from ground deformation, local geology, regional seismic network, and teleseismic joint epicenter determination have also been used to constrain the source model and the rupture process.The fault mechanism of the Hsingtai sequence was mainly strike-slip with a small component of normal dip-slip. The strikes of the four largest shocks range from ～ N26° to 30°E, approximately along strike of the major faults of the Shu-lu graben and the aftershock distribution. The source mechanisms can be explained with a NNW-SSE extensional stress and a NEE-SWW compressional stress acting in the area. The major shocks all had focal depths ～ 10 km.The four largest shocks in the sequence were characterized by a relatively simple and smooth dislocation time history. The durations of the far-field source time functions ranged from 3.5 to 5 s, while the rise times were all ～ 1 s. The seismic moments of the four largest earthquakes ranged from 1.43 × 10²⁵ to 1.51 × 10²⁶ dyne cm^?1. The fault sizes of the four events were very close. Assuming circular faults, the diameters of the four events were determined to be between 10 and 14 km. Stress drops varied from ～ 52 to 194 bars. A trend of increasing stress drop with earthquake size was observed.A survey of stress drop determinations for 15 major intraplate earthquakes shows that on the average the magnitude of stress drop of oceanic intraplate earthquakes and passive continental margin events is higher (～ 200 to several hundred bars) than that of continental intraplate earthquakes (～ 100 bars or less).     

Triggering of repeated earthquakes

Based on the analysis of the world’s earthquakes with magnitudes M ≥ 6.5 for 1960–2013, it is shown that they cause global-scale coherent seismic oscillations which most distinctly manifest themselves in the period interval of 4–6 min during 1–3 days after the event. After these earthquakes, a repeated shock has an increased probability to occur in different seismically active regions located as far away as a few thousand km from the previous event, i.e., a remote interaction of seismic events takes place. The number of the repeated shocks N(t) decreases with time, which characterizes the memory of the lithosphere about the impact that has occurred. The time decay N(t) can be approximated by the linear, exponential, and powerlaw dependences. No distinct correlation between the spatial locations of the initial and repeated earthquakes is revealed. The probable triggering mechanisms of the remote interaction between the earthquakes are discussed. Surface seismic waves traveling several times around the Earth’s, coherent oscillations, and global source are the most preferable candidates. This may lead to the accumulation and coalescence of ruptures in the highly stressed or weakened domains of a seismically active region, which increases the probability of a repeated earthquake.     

Seismogenic stress field beneath Mt. Etna (South Italy) and possible relationships with volcano-tectonic features

Shallow shear-type seismic activity occurring beneath the Etna volcano during 1990–1995 has been analysed for hypocenter locations, focal mechanisms and stress tensor inversion. The results have been examined jointly with Electronic Distance Measurements and tiltmeter data collected in the same period and reported in the literature. Significant seismicity located in the upper 10 km was found to be confined to the time intervals in which ground deformation data indicated inflation of the volcano edifice (e.g., the periods preceding the December 1991–March 1993 and August 1995–March 1996 eruptive phases). The shocks mostly occurred in a sector approximately centered on the crater area and elongated in the East–West direction. The causative seismogenic stress shows a low-dip East–West orientation of σ₁. In agreement with existing knowledge on relationships between local fault systems and magma uprise processes, the shallow seismicity in question is tentatively explained as being due to lateral compression by magma inside a nearly North–South system. The volcano deflation phase revealed by Electronic Distance Measurements and tilt data during the 1991–1993 major eruption was not accompanied by any significant shear-type shallow event. Below the depth of 10 km, the North–South prevailing orientation of σ₁ reflects the dominant role of the regional stress.     

Earthquake mechanism and predictability shown by a laboratory fault

Slip events generated in a laboratory fault model consisting of a circulinear chain of eight spring-connected blocks of approximately equal weight elastically driven to slide on a frictional surface are studied. It is found that most of the input strain energy is released by a relatively few large events, which are approximately time predictable. A large event tends to roughen stress distribution along the fault, whereas the subsequent smaller events tend to smooth the stress distribution and prepare a condition of simultaneous criticality for the occurrence of the next large event. The frequency-size distribution resembles the Gutenberg-Richter relation for earthquakes, except for a falloff for the largest events due to the finite energy-storage capacity of the fault system. Slip distributions, in different events are commonly dissimilar. Stress drop, slip velocity, and rupture velocity all tend to increase with event size. Rupture-initiation locations are usually not close to the maximum-slip locations.     

Relocation of Swarm Events in the June 1987 Earthquake Sequence from the Lake Aswan Area in Egypt

--We have examined the digital waveform data and relocated a number of events within the June 1987 earthquake swarm, which occurred beneath the northern part of Lake Aswan, 70 km southwest of the Aswan High Dam in Egypt. This swarm occurred between June 17th and 19th with a maximum magnitude event of "M_L"=3.5.¶Cross correlation between a chosen master and the analyzed events has been carried out on seismograms from stations of the Aswan network. The cross correlation demonstrates the presence of a difference in both the P wave ((t_p) and the S wave ((t_s) arrival times at each station in the network relative to the arrival times of the master event at the same stations. (t_p ranges between т.15 and 0.11 second, while (t_s ranges between т.17 and 0.11 second.¶The primary interpretation is that the se time differences represent an error in the manually picking arrival times. Then, (t_p and (t_s values for each event result from a change in the hypocentral parameters from those of the master event, assuming the P- and S-wave velocity distribution remains constant during the swarm activity. This interpretation enables us to determine the relative distribution of hypocenters with respect to the hypocentral location of the master event. We present the results from a swarm of 9 events demonstrating they originate from a nearly unique location, rather than the zone identified from the preliminary locations which used manually picked onset times.     

A layered model of the upper crust in the Aigion region of Greece,inferred from arrival times of the 2001 earthquake sequence

The western part of the Gulf of Corinth attracts attention due to its seismically active fault system and considerable seismic hazard. Detailed studies of the seismic activity of the region have been carried out especially as part of the so-called Corinth Rift Laboratory (CRL) Project. For standard earthquake locations, the CRL uses the HYPO algorithm and a special structural model that is composed of homogeneous layers (Rigo et al., 1996). This model was derived from a passive seismic experiment in a broader area around the western part of the Gulf. A significant part of the seismic activity is concentrated close to the town of Aigion, which was damaged by a strong earthquake in 1995. A sequence of smaller earthquakes occurred to the south of this town in the year 2001. In the present paper, we have used this sequence to derive an improved structural model for the region in the vicinity of the town of Aigion. This new model is based on the minimization of travel-time residuals. In particular, we used arrival times from a subset of 88 events recorded by at least 5 stations of the southern part of the CRL network, had magnitudes of over 2.3, and were recorded at the nearest station (station code AIO). A variant of the method of conjugate gradients has been used for this purpose. In comparison with the model derived by Rigo et al. (1996), the new model is characterized by a higher ν_P/ν_S velocity ratio and by higher velocities to a depth of about 7 km. The new model was derived with the aim to get more accurate locations of future events in the vicinity of the town of Aigion.     

Tectonic analysis of mine tremor mechanisms from the Upper Silesian Coal Basin

Fault network of the Upper Silesian Coal Basin (USCB) is built of sets of strike-slip, oblique-slip and dip-slip faults. It is a typical product of force couple which acts evenly with the parallel of latitude, causing horizontal and anti-clockwise movement of rock-mass. Earlier research of focal mechanisms of mine tremors, using a standard fault plane solution, has shown that some events are related to tectonic directions in main structural units of the USCB. An attempt was undertaken to analyze the records of mine tremors from the period 1992–1994 in the selected coal fields. The digital records of about 200 mine tremors with energy larger than 1×10⁴ J (M _L >1.23) were analyzed with SMT software for seismic moment tensor inversion. The decomposition of seismic moment tensor of mine tremors was segmented into isotropic (I) part, compensated linear vector dipole (CLVD) part and double-couple (DC) part. The DC part is prevalent (up to 70%) in the majority of quakes from the central region of the USCB. A group of mine tremors with large I element (up to 50%) can also be observed. The spatial orientation of the fault and auxiliary planes were obtained from the computations for the seismic moment DC part. Study of the DC part of the seismic moment tensor made it possible for us to separate the group of events which might be acknowledged to have their origin in unstable energy release on surfaces of faults forming a regional structural pattern. The possible influence of the Cainozoic tectonic history of the USCB on the recent shape of stress field is discussed.     

Some aspects of earthquake occurrence to the southeast of Cuba

Fifty-two months of data of small earthquakes in an area of 6000 km² to the southeast of Cuba were analyzed in order to determine the space distribution of the events and their temporal behaviour.It is shown that an earthquake swarm to the south of Guantanamo Bay can be regarded as a sequence of foreshocks preceding the earthquake of December 26, 1972. It was concluded that the swarm is distinguished by a clustering of earthquakes whose frequency distribution with time showed a bay form and that the main shock followed a quiet period during which the frequency of events decreased to its original low level.     

The Y. D. and climate abrupt events in the early and middle Holocene: Stalagmite oxygen isotope record from Maolan,Guizhou, China

Changes of climate, especially the abrupt changesof precipitation and temperature lead frequently tograve effect in human productive and living activity.Since the 1950s, results of monitoring indicate thatatmosphere CO2 concentration rises by 1.4×10?6g/gevery year. The atmospheric greenhouse effect leadsprobably to global warming since the 1940s, whichbrings furthermore special attention of scientists andrelevant government organizations in various coun-tries. In the near 20 years, the …     

Seismicity induced by gas production: II. Lithology correlated events,induced stresses and deformation

A 3D relocation technique permits precise locations of induced earthquakes. Geostatistical processing using the data of 87 boreholes provides the basis of a precise 3D structure, with a dome geometry. Conventional laboratory mechanical tests performed on deep rock samples (1000 m to 5000 m) define the rock properties at depths similar to those of the seismic events (1<M _L<4.2) that range from 1 to 7 km.In the studied period, most (85%), of the events were located above the gas reservoir, with very few located in the reservoir itself. Because the production parameters (50 MPa depletion of the gas pressure reservoir) are homogeneous throughout the gas field, the lateral inhomogeneity of the seismic rupture locations are a consequence of variations in the rheological response of the dome to the deformation induced by gas production.Here a ratio of two is found between the elastic modulus of the seismic rock matrix and the elastic modulus of the aseismic rock matrix. The contrast in strength is at least as great, if not greater. Repeated measured surface deformations involve the whole structure. Spatial and temporal deformations indicate that aseismic deformation is quantitatively the main process of this structural deformation. The heterogeneous stress pattern inferred fromP-axes of induced earthquakes disagrees with the tectonic regional stress field. The radial distribution ofP-axes towards the gas reservoir probably reflects the production induced deformation. The inferred deformation of the dome occurs in response to weak induced stresses.