Simulation of Seismicity in the Block-structure Model of Italy and its Surroundings

The numerical block-model of the lithosphere dynamics is used to simulate seismicity in Italy and its surroundings, based on the available structural and geodynamics information. The purpose of the study is to understand which are the tectonic processes that control the main features of the observed seismicity and the kinematics of the region. The influence of the rheology of the fault systems is studied as well. The model we use differs from other modeling approaches in that it simulates earthquakes and hence it possibly relates to seismicity and geodynamics. The model provides an effective capability to include the set of documented constraints supplied by widely available earthquake catalogs. This is done by means of the comparison of the GR relation, of the focal mechanisms and of the space distribution for observed and computed seismicity. The region is modeled as a system of perfectly rigid blocks, separated by infinitely thin fault planes, in viscoelastic interaction between themselves and with the underlying medium. The movement of the boundary blocks and of the underlying medium determines the motion of the blocks. The synthetic seismicity obtained with the defined block-model is similar to the observed one for the most seismically active areas. A linear frequency-magnitude (FM) relation (Gutenberg-Richter law) is obtained for synthetic earthquakes; the slope (b-value) of the FM plot appears larger for the synthetic seismicity than for the observed one. Nevertheless, the b-value is essentially larger in northern and central Italy than that in southern Italy, both in the model and in the observations. The analysis of the source mechanisms of the synthetic earthquakes shows a good agreement with the observations. In the model normal faulting is typical for the Apennines, the eastern edge of Sicily and the Calabrian arc, while reverse faulting takes place at the northwestern boundary of the Adriatic Sea, in the southern Alps and along the eastern edge of the Adria, along the Dinarides. The model correctly reproduces the extension zone along the Apennines and the contraction zone along the northwestern boundary of the Adriatic Sea; the counter-clockwise rotation of the Adria is mimed. The resulting movements of the blocks are in overall agreement with GPS (Global Positioning System) observations. The results of the modeling experiments suggest that the main features of dynamics and seismicity in the central Mediterranean region cannot be satisfactorily explained as a consequence of Africa and Eurasia convergence only; the passive subduction in the Calabrian arc and the different rheology of faults are essential as well.     

Formation of the superdeep South Caspian basin: subsidence driven by phase change in continental crust

The large hydrocarbon basin of South Caspian is filled with sediments reaching a thickness of 20–25 km. The sediments overlie a 10–18 km thick high-velocity basement which is often interpreted as oceanic crust. This interpretation is, however, inconsistent with rapid major subsidence in Pliocene-Pleistocene time and deposition of 10 km of sediments because the subsidence of crust produced in spreading ridges normally occurs at decreasing rates. Furthermore, filling a basin upon a 10–18 km thick oceanic crust would require twice less sediments. Subsidence as in the South Caspian, of ≥20 km, can be provided by phase change of gabbro to dense eclogite in a 25–30 km thick lower crust. Eclogites which are denser than the mantle and have nearly mantle P velocities but a chemistry of continental crust may occur beneath the Moho in the South Caspian where consolidated crust totals a thickness of 40–50 km. The high subsidence rates in the Pliocene-Pleistocene may be attributed to the effect of active fluids infiltrated from the asthenosphere to catalyze the gabbro-eclogite transition. Subsidence of this kind is typical of large petroleum provinces. According to some interpretations, historic seismicity with 30–70 km focal depths in a 100 km wide zone (beneath the Apsheron-Balkhan sill and north of it) has been associated with the initiation of subduction under the Middle Caspian. The consolidated lithosphere of deep continental sedimentary basins being denser than the asthenosphere, can, in principle, subduct into the latter, while the overlying sediments can be delaminated and folded. Yet, subduction in the South Caspian basin is incompatible with the only 5–10 km shortening of sediments in the Apsheron-Balkhan sill and south of it and with the patterns of earthquake foci that show no alignment like in a Benioff zone and have mostly extension mechanisms.     

Improving automatic location procedure by waveform similarity analysis: An application in the South Western Alps (Italy)

The accuracy of automatic procedures for locating earthquakes is influenced by several factors such as errors in picking seismic phases, network geometry, modeling errors and velocity model uncertainties. The main purpose of this work is to improve the performances of the automatic procedure employed for the “quasi-real-time” location of seismic events in North Western Italy by developing a procedure based on a waveform similarity analysis and by using only one seismic station.To detect “earthquake families” a cross-correlation technique was applied to a data set of seismic waveforms recorded in the period 1985-2002, in a small test area (1600 km²) located in the South Western Alps (Italy). Normalized cross-correlation matrices were calculated using about 2700 seismic events, selected on the basis of the signal to noise ratio, manually picked and located by using the Hypoellipse code. The waveform similarity analysis, based on the bridging technique, allowed grouping about 65% of the selected events into 80 earthquake families (multiplets) located inside the area considered. For each earthquake family a master event is selected, manually re-picked and re-located by using Hypoellipse code. Having chosen a reference station (STV) on the basis of the completeness of the available data set, an automatic procedure has been developed with the aim of cross-correlating new seismic recordings (automatically picked) to the waveforms of the events belonging to the detected families. If the new event is proved to belong to a family (on the basis of the cross-correlation values), its hypocenter co-ordinates are defined by the location of the master event of the associated family. The performance of the proposed procedure is tested and demonstrated using a data set of 104 selected earthquakes recorded in the period January 2003-June 2004 and located in the test area. The automatic procedure is able to locate, associating events with the multiplets detected by the waveform similarity analysis, about 50% of the test events, almost independently of the accuracy of the automatic phase picker and without the biasing of the network geometry and of the velocity model uncertainties.     

Seismic doublets and multiplets at Polish coal and copper mines

The following criteria for selection of doublets at Polish coal mines were accepted: the difference in magnitude (based on seismic moment) of two events not larger than 0.15, the distance between their hypocenters not greater than 150 m, and the time interval between their occurrence not longer than 10 days. Similarly, the criteria for seismic events at copper mines are: the difference in magnitude not exceeding 0.15, the distance not greater than 200 m, and the time interval not longer than 20 days. Seismic events from the Wujek and Ziemowit coal mines that occurred between 1993 and 1995, and seismic events from the Polkowice copper mine that occurred between 1994 and 1996 and from the Rudna copper mine that occurred between 1994 and 2004 were considered. Their source parameters and focal mechanisms were known in most cases from previous studies. Altogether 108 seismic pairs from coal mines and 118 pairs from copper mines were found, forming doublets, triplets and quadruplets, within the magnitude range from 0.7 to 3.5. The distance and time intervals between two events forming pairs are not dependent on magnitude of these events. The focal mechanism of seismic events forming pairs is similar in over 60 percent of pairs at coal mines and in about one third of pairs at copper mines. Spatial distributions of doublets in particular sections of coal and copper mines display dominant linear trends, characteristic for a given area, which are often in conformity with the direction of nodal planes determined by fault plane solution of one or both the events forming a doublet. In such cases, the rupture plane can be discriminated among the nodal planes.     

Seismic Hazard in Terms of Spectral Accelerations and Uniform Hazard Spectra in Northern Algeria

Seismic hazard in terms of spectral acceleration (SA) has been estimated for the first time in northern Algeria. For this purpose, we have used the spatially-smoothed seismicity approach. The present paper is intended to be a continuation of previous work in which we have evaluated the seismic hazard in terms of peak ground acceleration (PGA) using the same methodology. To perform these evaluations, four complete and Poissonian seismic models have been used. One of them considers earthquakes with magnitudes above M_S 6.5 in the last 300 years, that is, the most energetic seismicity in the region. Firstly, seismic hazard maps in terms of SA, at periods of 0.1, 0.2, 0.3, 0.4, 0.5, 1.0, 1.5 and 2.0 sec, with 39.3% and 10% probability of exceedance in 50 years, have been obtained. Therefore, uniform hazard spectra (UHS) are computed and examined in detail for twelve of the most industrial and populated cities in northern Algeria. All the reported results in this study are for rock soil and 5% of damping. It is noteworthy that, in the seismic hazard maps as well as in the UHS plots, we observe maximum SA values in the central area of the Tell. The higher values are reached in the Chleff region (previously El Asnam), specifically around the location of the destructive earthquakes of September 9, 1954 (M_S 6.8), and October 10, 1980 (M_S 7.3). These maximum values, 0.4 g and 1.0 g, are associated with periods of about 0.2 and 0.3 sec for return periods of 100 and 475 years, respectively.     

The Applicability of Modern Methods of Earthquake Location

We compare traditional methods of seismic event location, based on phase pick data and analysis of events one-at-a-time, with a modern method based on cross-correlation measurements and joint analysis of numerous events. In application to four different regions representing different types of seismicity and monitored with networks of different station density, we present preliminary results indicating what fraction of seismic events may be amenable to analysis with modern methods. The latter can supply locations ten to a hundred times more precise than traditional methods. Since good locations of seismic sources are needed as the starting point for so many user communities, and potentially can be provided due to current improvements in easily-accessible computational capability, we advocate wide-scale application of modern methods in the routine production of bulletins of seismicity. This effort requires access to waveform archives from well-calibrated stations that have long operated at the same location.     

强震孕育过程中区域地震活动演化与活动地块的关系

通过对中国大陆及邻区4个Ⅰ级活动地块中1972—2003年全部地震资料的系统分析,认为强震中期阶段典型的地震活动演化图像是中等地震聚集活动,震前震中区附近的地震密度和有效速率增高。文中用中等地震增长率来刻划强震孕育过程中中等地震活动的加速增长,结果表明:中等地震聚集活动密度的增长率能够较好地反映强震前中期异常变化特征,对各活动地块的分析表明华北地块中等地震增长率达到180%以上时,未来1～3a中等地震异常聚集区及其附近会发生强震;而西域、青藏和滇缅地块中等地震增长率达到250%以上时,应注意中等地震异常聚集区及其附近强震的发生。因此,中等地震增长率可以作为中期预报的定量化指标     

依兰-伊通断裂北段地震活动性研究

依据地质构造背景，应用测震学的基本方法，讨论依兰-伊通断裂北段的地震活动特征，由地震活动在时空强显示的一些规律性，认为近期该区域中强地震活动可能处于较弱势状态。