Convolutive independent component analysis for processing massive datasets: a case study at Campi Flegrei (Italy)

A novel procedure is proposed to analyse continuous seismic signal on hourly scales to have a prompt discrimination among the different sources. Specifically, this approach is applied to a massive dataset recorded at Campi Flegrei caldera during the year 2006 when a swarm of volcano-tectonic earthquakes occurred. The convolutive independent component analysis is adopted to obtain a clear separation among meteo-marine microseism, anthropogenic noise, hydrothermal tremor in the absence of volcano-tectonic activity, whereas in non-stationary conditions a contribution connected to the corner frequency of the earthquakes emerges. A coarse-grained variable to be monitored continuously is introduced, i.e. the frequency associated with the maximum amplitude of the power spectral density of the deconvolutive independent components. That parameter is sensitive to the variation in the frequency bands of interest (e.g. that corresponding to the corner frequencies of volcano-tectonic events) and can be used as marker of the insurgence of seismic activity.

     

Seismic source parameters of large earthquake in Vrancea

The seismic source parameters for five large Vrancea earthquakes are calculated : corner frequency, focal radius, seismic moment and stress drop.The multiplicity character of some events is put in evidence, using spectral and time domain analysis.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

A geophysical multi-parametric analysis of hydrothermal activity at Dallol,Ethiopia

During December 2003, three seismic stations were installed close to the hornitos of the hydrothermal system at Dallol, complemented by radiometer and infrasonic measurements. A combined geophysical data set was collected for about three days. During this period thermal, seismic and acoustic records indicate the presence of two regimes characterized by a different energy distribution in frequency. Few volcano-tectonic events appear superimposed to the continuous hydrothermal tremor. The continuous data indicate variable shallow processes most likely related with variations in temperature and degassing processes within the shallow geothermal system. This alternation of low and high regimes shows significant similarities with other volcanic systems of different nature, although at Dallol the transition is more evident in the thermal than in the seismic and acoustic data.     

Double seismic zone beneath the middle of Hokkaido, Japan, in the southwestern side of the Kurile Arc

The vertical section of microearthquakes, determined accurately by using the Hokkaido University network, shows two dipping zones (the double seismic zone) 25–30 km apart in the depth range of 80–150 km beneath the middle of Hokkaido in the southwestern side of the Kurile arc. Hypocentral distribution of large earthquakes (m_b > 4) based on the ISC (International Seismological Centre) bulletin also shows the double seismic zone beneath the same region. The hypocentral distribution indicates that the frequency of events occurring in the lower zone is four times greater than that in the upper zone. The difference in seismic activity between the two zones beneath Hokkaido is in contrast with the region beneath northeastern Honshu in the northeastern Japan arc.Composite focal mechanisms of microearthquakes and individual mechanisms of large events mainly characterize the down-dip extension for the lower zone as is observed beneath northeastern Honshu. For the upper zone, however, the stress field is rather complex and not necessarily similar to that beneath northeastern Honshu. This may be considered to indicate the influence of slab contortion or transformation in the Hokkaido corner between the Kurile and the northeastern Japan arcs.     

Source parameters and scaling relations for small earthquakes in the Kachchh region of Gujarat,India

The scaling relationships for stress drop and corner frequency with respect to magnitude have been worked out using 159 accelerograms from 34 small earthquakes (M _w 3.3–4.9) in the Kachchh region of Gujarat. The 318 spectra of P and S waves have been analyzed for this purpose. The average ratio of P- to S-wave corner frequency is found to be 1.19 suggestive of higher corner frequency for P wave as compared to that for S wave. The seismic moments estimated from P waves, M ₀(P), range from 1.98 × 10¹⁴ N m to 1.60 × 10¹⁶ N m and those from S waves, M ₀(S), range from 1.02 × 10¹⁴ N m to 3.4 × 10¹⁶ N m with an average ratio, M ₀(P)/M ₀(S), of 1.11. The total seismic energy varies from 1.83 × 10¹⁰ J to 2.84 × 10¹³ J. The estimated stress drop values do not depend on earthquake size significantly and lie in the range 30–120 bars for most of the events. A linear regression analysis between the estimated seismic moment (M ₀) and corner frequency (f _c) gives the scaling relation M ₀ f _c ³  = 7.6 × 10¹⁶ N m/s³. The proposed scaling laws are found to be consistent with similar scaling relations obtained in other seismically active regions of the world. Such an investigation should prove useful in seismic hazard and risk-related studies of the region. The relations developed in this study may be useful for the seismic hazard studies in the region.     

用小波变换讨论地震中的多普勒现象

地震中多普勒效应可以确定地震的破裂面等,说明对多普勒效应的研究有实际意义,但目前确定地震中是否存在多普勒效应的方法并不成熟。在研究多普勒效应空间分布规律的基础上,提出用小波变换确定地震中是否存在多普勒效应的方法。选择位于汶川地震断层滑动前方的若干台站对台站最初时段的地震记录进行小波变换时,发现随着震中距的增加,小波谱高频幅值明显大于小波谱低频幅值;说明虽然存在介质对地震波的吸收衰减作用,但多普勒效应的存在仍使得小波谱高频幅值增大。选择与汶川地震断层垂直方位的若干台站对台站最初时段的地震记录进行小波变换时,发现随着震中距的增加,小波谱高频幅值迅速降低,震中距大到一定程度后低频部分的小波谱幅值会明显大于高频部分的小波幅值;说明在该方向上,介质对地震波的高频吸收衰减起主要作用,没有发生多普勒效应。     

Dynamics of geophysical medium from power spectral density of microseisms before and after earthquakes: case study of Bureya massif,Amur region

The time–frequency parameters of weak earthquakes and microseisms are studied. The qualitative and quantitative relationships of the power spectral density of seismic waves are established; these relationships vary in terms of frequency spectrum for the areas of the northeastern framing of the Tan Lu fault system and Bureya massif. The reason for these differences could be the influence of the crustal geological structure near the observations points. Resulting from the time–frequency analysis of weak earthquakes, two sites located in the Tan Lu fault zone show an increase in the power spectrum at frequencies of 1 to 5 Hz and the resonant excitation of the medium at high frequencies (12.5–35 Hz) for the area of the Bureya hydroelectric power station during the travelling of seismic waves from weak earthquakes. A longer attenuation of the power spectral density of seismic waves at high frequencies is noted, and this may occur due to resonant excitation of the medium and the influence of the dam on the geological medium. On the contrary, this effect was not observed at the second site located near Lake Udyl. It is shown that the increase in power can be attributed to the interaction between seismic waves and spatial inhomogeneities in the Earth’s crust.     

Local seismic hazard assessment in areas of weak to moderate seismicity—case study from Eastern Germany

Generally the seismic hazard of an area of interest is considered independent of time. However, its seismic risk or vulnerability, respectively, increases with the population and developing state of economy of the area. Therefore, many areas of moderate seismic hazard gain increasing importance with respect to seismic hazard and risk analysis. However, these areas mostly have a weak earthquake database, i.e., they are characterised by relative low seismicity and uncertain information concerning historical earthquakes. In a case study for Eastern Thuringia (Germany), acting as example for similar places in the world, seismic hazard is estimated using the probabilistic approach. Because of the lack of earthquakes occurring in the recent past, mainly historical earthquakes have to be used. But for these the actual earthquake sources or active faults, needed for the analysis, are imprecisely known. Therefore, the earthquake locations are represented by areal sources, a common practice. The definition of these sources is performed carefully, because their geometrical shape and size (apart from the earthquake occurrence model) influence the results significantly. Using analysis tools such as density maps of earthquake epicentres, seismic strain and energy release support this. Oversizing of areal sources leads to underestimation of seismic hazard and should therefore be avoided. Large location errors of historical earthquakes on the other hand are represented by several alternative areal sources with final superimposition of the different results. In a very similar way information known from macroseismic observations interpreted as source rather than as site effects are taken into account in order to achieve a seismic hazard assessment as realistic as possible. In very local cases the meaning of source effects exceeds those of site effects very likely. The influence of attenuation parameter variations on the result of estimated local seismic hazard is relatively low. Generally, the results obtained by the seismic hazard assessment coincide well with macroseismic observations from the thoroughly investigated largest earthquake in the region.     

Seismicity and source parameters of moderate earthquakes in Sikkim Himalaya

In this study, we accurately relocate 360 earthquakes in the Sikkim Himalaya through the application of the double-difference algorithm to 4?years of data accrued from a eleven-station broadband seismic network. The analysis brings out two major clusters of seismicity??one located in between the main central thrust (MCT) and the main boundary thrust (MBT) and the other in the northwest region of Sikkim that is site to the devastating Mw6.9 earthquake of September 18, 2011. Keeping in view the limitations imposed by the Nyquist frequency of our data (10?Hz), we select 9 moderate size earthquakes (5.3????Ml????4) for the estimation of source parameters. Analysis of shear wave spectra of these earthquakes yields seismic moments in the range of 7.95?×?10²¹ dyne-cm to 6.31?×?10²³ dyne-cm and corner frequencies in the range of 1.8?C6.25?Hz. Smaller seismic moments obtained in Sikkim when compared with the rest of the Himalaya vindicates the lower seismicity levels in the region. Interestingly, it is observed that most of the events having larger seismic moment occur between MBT and MCT lending credence to our observation that this is the most active portion of Sikkim Himalaya. The estimates of stress drop and source radius range from 48 to 389?bar and 0.225 to 0.781?km, respectively. Stress drops do not seem to correlate with the scalar seismic moments affirming the view that stress drop is independent over a wide moment range. While the continental collision scenario can be invoked as a reason to explain a predominance of low stress drops in the Himalayan region, those with relatively higher stress drops in Sikkim Himalaya could be attributed to their affinity with strike-slip source mechanisms. Least square regression of the scalar seismic moment (M ₀) and local magnitude (Ml) results in a relation LogM ₀?=?(1.56?±?0.05)Ml?+?(8.55?±?0.12) while that between moment magnitude (M _w ) and local magnitude as M _w ?=?(0.92?±?0.04)Ml?+?(0.14?±?0.06). These relations could serve as useful inputs for the assessment of earthquake hazard in this seismically active region of Himalaya.     

全球大地震时空分布与动力学机制的初步研究

文中根据最新的全球大震资料划分了 10 0a来全球特大地震的时空演化过程 ,表现为环太平洋带与近纬向活动带两种强震分布图像 ,以 2 0a左右的时间段交替出现。逆冲和引张型地震释放了全球强震能量的主要部分。逆冲型地震主要发生在环太平洋地震带上 ,在喜马拉雅碰撞边界和巽他弧也有发生。引张型地震主要出现在洋脊和裂谷地区。而走滑型地震主要分布在东亚大陆内部、北美板块西边界和加勒比板块的周边、西南太平洋带、以及地中海 -喜马拉雅带部分区域 ,呈近纬向分布。地幔对流的垂直运动与水平运动分量的比率 ,决定了逆冲、引张与走滑型强震的能量释放份额。地幔运动在地表体现在逆冲、引张与走滑型强震空间分布的区域差异性。这种差异在时间上的演化形成了全球强震活动环太平带与近纬向分布图像的交替出现。     

Errors in expected human losses due to incorrect seismic hazard estimates

Seismic hazard maps are constructed by extrapolating from the frequency of small earthquakes, the annual probability of large, infrequent, earthquakes. Combining the potential contribution from all seismically active volumes, one calculates the peak ground acceleration with a probability to be exceeded by 10?% in 50?years at any given point. The consequential risk, the losses to be expected, derives from the damage the calculated shaking causes to buildings, and the impact on occupants due to collapsing structures. We show that the numbers of fatalities in recent disastrous earthquakes were underestimated by the world seismic hazard maps by approximately two to three orders of magnitude. Thus, seismic hazard maps based on the standard method cannot be used to estimate the risk to which the population is exposed due to large earthquakes.     

京津冀地区历史地震事件时空特征研究*

受环太平洋地震带影响,华北平原地区地震频发,尤其是处于中国首都经济圈的京津冀地区的地震事件备受关注。通过对历史文献资料及地震台网记录中的地震事件统计、分析,重建该地区地震事件历史并获取其潜在的空间分布特征及时间规律,对未来地震事件的早期预警具有重要参考意义。分析结果表明,公元前231年至公元2018年期间京津冀地区发生的1044起地震事件中,以有感地震和中强地震为主,小地震、强烈地震以及大地震发生频次较低。地震记录完整性分析结果表明,除小地震外,其他等级地震记录自公元1400年以来基本完整。在空间分布上,京津冀地区历史地震呈“T”字形分布,沿1条北西—南东走向地震带和1条北东—南西走向地震带分布。在时间上,京津冀地区地震事件呈现出阶段性的变化,在公元1480—1680年间以及1950年以来2个时间段内较为活跃,发生频率较高,频谱分析结果进一步表明地震记录存在45年的复发周期。在月际尺度上,地震事件同样存在季节性差异且多发于夏秋季节,同时地震密集区域在年内呈现出自西向东迁移的现象。最后,根据历史地震事件发生的时间规律,在未来一段时间内京津冀地区仍将处于地震活跃期,存在发生强震的风险。     

Mud volcano eruptions and earthquakes in the Northern Apennines and Sicily, Italy

The relations between earthquakes and the eruption of mud volcanoes have been investigated at the Pede–Apennine margin of the Northern Apennines and in Sicily. Some of these volcanoes experienced eruptions or increased activity in connection with historical seismic events, showing a good correlation with established thresholds of hydrological response (liquefaction) to earthquakes. However, the majority of eruptions have been documented to be independent of seismic activity, being mud volcanoes often not activated even when the earthquakes were of suitable magnitude and the epicentre at the proper distance for the triggering. This behaviour suggests that paroxysmal activity of mud volcanoes depends upon the reaching of a specific critical state dictated by internal fluid pressure, and implies that the strain caused by the passage of seismic waves can activate only mud volcanoes in near-critical conditions (i.e., close to the eruption). Seismogenic faults, such as the Pede–Apennine thrust, often structurally control the fluid reservoirs of mud volcanoes, which are frequently located at the core of thrust-related folds. Such an intimate link enables mud volcanoes to represent features potentially suitable for recording perturbations associated with the past and ongoing tectonic activity of the controlling fault system.     

KINEMATICS OF QUATERNARY EXTENSION IN THE PAMIR

KINEMATICS OF QUATERNARY EXTENSION IN THE PAMIR1　BurtmanVS ,MolnarP .GeologicalandgeophysicalevidencefordeepsubductionofcontinentalcrustbeneaththePamir[J].SpecPapGeolSocAm ,1993,2 81:76 . 2　StreckerMR ,FrischW ,HamburgerMW ,etal.QuaternarydeformationintheEasternPamirs ,TadzhikistanandKyrgyzs tan[J].Tectonics,1995 ,14(5 ) :10 6 1～ 10 79.…     

龙门山断层破裂的频率非平稳特性

为研究地震过程中的频率非平稳特性,对近年来龙门山断层发生的两次大地震:汶川大地震和芦山大地震的近断层地震记录进行频谱分析。结果表明:相对于芦山地震有较大走向滑动分量的汶川地震,大多数位于汶川地震断层滑动前方的台站接收到更高的频率成分,位于断层滑动后方的台站接收到的地震波频率较低;尽管芦山地震断层相对汶川地震有较小的走向滑动分量,但仍然可以得出与汶川地震相同的结论,不同的是虽然芦山地震沿断层面向上方向分量大,但是其同一台站东西、南北、竖直三方向分量记录幅值相当。把芦山地震三分量记录变换到走向和沿断层面向上方向,证实了沿断层面向上方向高频成分更丰富。在断层滑动前方接收到的地震波频率较高,在断层滑动后方接收到的频率较低,这正是多普勒效应影响的结果。由于多普勒效应的客观存在,其对频率非平稳特性的影响与震源、传播路径和场地效应一样具有普遍性;所以,工程场地接收到的地震波的频率不仅取决于震源、传播路径、场地效应,还取决于断层滑动速度(多普勒效应)。     

Earthquakes associated with the Clark Hill reservoir, South Carolina — A case of induced seismicity

Following the McCormick County, South Carolina, local magnitude 4.3 earthquake on August 2, 1974, continuous seismic activity has been observed in the area. The epicenters are located within 3 km from the Clark Hill reservoir. The shallow seismic activity appears to be related to water level fluctuations in the reservoir, and follows them by about two days. The frequency of earthquakes is related to the rate of change of water level, while the energy release is seen to depend on the water level itself. Pore pressure fluctuations at focal depths in a regime of high tectonic stresses is the suggested cause for the triggering of earthquakes.     

A scaling model for quantification of earthquakes in and near Japan

A simple method is developed to determine seismic moments of earthquakes. The method is qualified through criteria such as simplicity of calculations, coverage of wide magnitude range, and insensitivity to detailed instrumental response. The method is applied to 163 major earthquakes which occurred underneath Japan and the Japan Sea in the time from 1926 to 1977. Magnitudes of these earthquakes, which have been determined by the Japan Meteorological Agency, (M_JMA) cover the range from 4.3 to 7.5. At first, source spectra are analyzed through a very simple way introducing two new parameters: characteristic period T_c and seismic-moment factor M_c. The former is defined as an average value of apparent periods of seismic waves with the maximum trace amplitude at many stations. The latter is an average of products of maximum trace amplitude and its apparent period multiplied by epicentral distance. It is shown that T_c corresponds to the period of the corner frequency of an earthquake and M_c to the seismic-moment density at the period of T_c. A scaling model of earthquake source spectra is presented which satisfies the empirical relations between the surface-wave magnitude M_s and M_JMA, and M_JMA and the body-wave magnitude m_b. Those relations are independent of the Gutenberg and Richter relation between M_s and m_b, because M_JMA is determined from maximum amplitudes of seismic waves with a period of about 4 sec. The static seismic moment of each earthquake can be estimated from calculated M_c using the source spectra of the scaling model. Seismic moments of 18 earthquakes determined by conventional analyses from near- and/or far-field observations are consistent with static seismic moments thus estimated over the range from 2 × 10²³ to 3 × 10²⁷ dyne cm. This shows the potential in practice of the present method, especially in the routine processing of seismic data.     

An overview on the seismic zonation and microzonation studies in India

The present study presents a review on the progressive development of the seismic zonation map of India both from official agencies and also from independent individual studies. The zonation map have been modified and updated regularly with the occurrence of major destructive earthquakes over the years in the Indian subcontinent with the addition of new data. This study discusses the criteria chosen for the progressive zonation and the major earthquakes that were responsible for retrospection of the earlier published maps. The seismic zonation maps of India have also been prepared by various independent workers by adopting different approaches to achieve the purpose of the zonation. Despite the endeavors from various sources to provide a solution for the problem of earthquake hazards in India, there were many limitations on the zonation map as it gives the picture at a regional scale mostly on the bedrock level without addressing the local site conditions. But nevertheless, the seismic zonation map gives basic guidelines for any region to know the hazard scenario and if any city or urban population is under threat from seismic point of view, further site specific seismic microzonation may be carried out. In the International scenario, the Global Seismic Hazard Assessment Program (GSHAP) in 1999 prepared a hazard map for world in terms of peak ground acceleration (PGA) with a 10% probability of exceedance in 50 years, but it turned out to be an underestimation of the hazard parameter when compared with the observed PGA. To tackle the problem of seismic hazards, there was a need to have a detail study on the local site conditions in terms of its geological, geophysical and geotechnical properties. With the advent of better instrumentation and knowledge on the mechanics of earthquakes, it was possible to identify zones of hazards at a local level and this gives rise to the study of seismic microzonation. Seismic microzonation work has been carried out in India in some of the strategic important mega cities and industrial build up that has the potential of being damaged from future earthquakes, as has been shown in the past. Though the microzonation map is not the final output map, as it can still be updated at later stage with more input data, it does provide a more realistic picture on the site specific seismic hazard.