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.     

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.     

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.     

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.     

Recognition of Strong Earthquake Prone Areas in the Altai–Sayan–Baikal Region

The locations of areas prone to strong earthquakes (M ≥ 6.0) in the Altai–Sayan–Baikal region are determined. Based on a scheme of morphostructural zoning of the region and by using the CORA-3 pattern recognition algorithm, all intersections of morphostructural lineaments are separated into two classes: the highly seismic intersections in the vicinities of which strong earthquakes can occur and low seismic in the vicinities of which only earthquakes with M < 6.0 are possible. Recognition was performed for the vectors the components of which were measured values of the geological–geophysical characteristics describing the respective intersection. The result obtained allows the zones of high seismic hazard to be identified more reliably 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.     

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.     

Computation of sliding displacements of seawalls under earthquake conditions

To understand the serviceability aspects of seawalls, it is essential to study the permanent displacements of seawalls that occur during the earthquakes. Studies in the existing literature have concentrated on displacements of retaining walls with dry backfills; to the authors’ observation there is no specific analytical investigation devoted to the earthquake-induced displacements of retaining walls with submerged backfills. This paper focuses on sliding displacements of gravity type seawall retaining a submerged backfill under active earth pressure condition during the earthquakes. The threshold seismic acceleration coefficients required for initiation of sliding and the amount of sliding displacement due to seismic loading are calculated by adopting Newmark’s sliding block method. One of the prime features of the study is the estimation of seismic inertia forces in the submerged soil and wall applying the modified pseudo-dynamic method. The comparison of the results obtained using the proposed analytical formulation with the existing literature found to be in good agreement. A comprehensive parametric study has been conducted to understand the effects of different parameters such as seismic horizontal and vertical acceleration coefficients, soil and wall friction angles, width of the wall, wall inclination and excess pore water pressure ratio.

     

Occurrence Apparent Velocities for Identification and Quantification of Space–Time Clustering Precursory to a Large Earthquake. Application to Large (M > 7.0) Earthquakes in Southern California and Northern Baja California

Space–time seismic clusters, localized bursts of seismic activity, are a feature of background seismicity before the occurrence of large earthquakes, a feature that agrees with observations of diminishing Gutenberg–Richter b-value, fractal dimension, and entropy, and is therefore suggestive of high stress. However, identification and quantification of these space–time clusters, particularly when they are small, is not an easy task and requires a priori assumptions. A novel method for space–time cluster identification, based on an extension of the concept of apparent velocities, is proposed because space–time clusters in the background seismicity have a particular signature in the apparent velocity domain. The contents of histogram peaks due to clusters in the apparent velocity histogram can be used to quantify the cluster activity compared with null hypothesis levels. Identification of the earthquakes corresponding to the apparent velocities in the peaks allows identification of cluster activity in time and space. Apparent velocity peaks do appear in real catalog data for southern California and northern Baja California before the Landers 1992 M = 7.3, Hector Mine 1999 M = 7.1, El Mayor-Cucapah 2010 M = 7.2, and Ridgecrest 2019 M = 7.1 earthquakes, and they appear only within 15 to 25 years before the occurrence of large earthquakes. They are not observed either long before the large earthquakes or after them, and hence could be related to high local states of stress and be of value as a possible precursory observable.

     

Modeling of source parameters and moment tensors of local earthquakes occurring in the eastern Indian shield

Earthquake source parameters and crustal Q are being estimated simultaneously through the inversion of S-wave displacement spectra from three-component recordings of ten local cratonic intraplate earthquakes from 3-6 broadband stations in the eastern Indian shield, wherein, an iterative Levenberg-Marquardt inversion technique is used. The estimated seismic moment (M_o) and source radii (r) vary from 7.4 x 10¹² to 7.1 x 10¹⁴ N-m and 144.2 to 211.3 m, respectively, while estimated stress drops (Δσ) and multiplicative factor (E_mo) values range from 0.11 to 4.13 MPa and 1.33 to 2.16, respectively. The corner frequencies range from 6.23 to 8.62 Hz while moment magnitudes vary from 2.44 to 3.57. The radiated seismic energy and apparent stresses range from 8.3 x 10⁶ to 2.0 x 10¹⁰ Joules and 0.06 to 0.94 MPa, respectively, wherein the estimated corner frequencies and seismic moment satisfy the relation Mo ∞ f _c ^–(3+ε) for ε = 12.7. Thus, the source scaling of these events clearly deviates from the self-similarity i.e. f^–3. Estimated Zuniga parameters reveal that all selected events satisfy the partial stress drop model, which is in good agreement with the global observations. Our estimated crustal S-wave quality factors vary from 1091 to 4926 with an average of 3006, suggesting a less heterogeneous crustal structure underlying the study region.We also perform moment tensor inversion of five selected local events using ISOLA software, which reveals that the dominant deformation mode for the eastern Indian shield is left-lateral strike slip motion with minor normal dip-slip component on an almost vertical plane. This observation suggests that neotectonic vertical movements might have played a key role in generating these earthquakes. Our modeling also depicts that the seismically mildly active Singhbhum shear zone and Eastern Ghats mobile belt are characterized by the left-lateral strike motion while two events in the Chotanagpur half graben belt suggest a normal dip-slip motion along a south dipping plane. A north-south orientation of P-axis is found to be dominant in the area, which is consistent with the prevailing north–south compression over the Indian plate.     

Inelastic displacement demand of RC buildings subjected to earthquakes generated by intermediate-depth Vrancea seismic source

Evaluating inelastic displacement demand of structures exposed to seismic hazard is required for the design of new buildings as well as for seismic risk assessment of existing structures. Most of the buildings are designed to withstand strong earthquakes by responding in the nonlinear range. Having special parts of the structure designed to develop a stable hysteretic behaviour allows the structure to deform in order to accommodate the displacement demand imposed by strong ground motions. This paper is centred on finding a correspondence between the maximum elastic and inelastic displacement responses of the single degree of freedom (SDOF) systems subjected to earthquakes generated by Vrancea seismic source. Vrancea intermediate-depth earthquakes are responsible for the seismic hazard throughout Romanian territory. They have distinctive features, such as large displacement demand and large predominant periods, which makes Romania a special seismic environment. Using a database of Romanian and Japanese strong ground motions generated by intermediate-depth earthquakes and performing nonlinear dynamic analysis on the SDOF oscillators following the Takeda model, this study estimates the inelastic to elastic displacement ratio of reinforced concrete systems. Soil conditions, epicentral distance and magnitude influence on inelastic response is analysed using constant ductility response spectra. The main findings of the study are: the local increase of the inelastic to elastic displacement ratio for type C soil (Eurocode 8 classification) for large magnitude earthquakes and the significant effect of soil conditions on the inelastic response of the SDOF systems. The inelastic amplification was evaluated using a functional form depending on system ductility, soil conditions and earthquake magnitude.

     

South Australian earthquakes, 1980–92

The seismicity of South Australia over the period 1980–92 is presented as a follow‐up to earlier studies. The South Australian seismic network has undergone a significant expansion in the last decade, and with it an increase in the number and precision of located earthquakes. The distribution of recent seismic activity is similar to the historical pattern of earthquakes and the previous instrumental seismicity maps, all of which show the three main areas as being the Flinders‐Mt Lofty Ranges, Eyre Peninsula, and the southeast. The one notable exception in the recent study is the presence of earthquake activity in the Musgrave Block, a previously aseismic region. Intensity characteristics are reported for earthquakes that were sufficiently widely felt. Fault plane solutions for three Flinders Ranges earthquakes (previously unpublished) are also presented; the focal mechanisms are consistent with predominant northeast‐southwest compression. The seismic moment method was used to estimate the seismic risk for the major population centres in terms of probability of exceedance of seismic intensity within a given period. These estimates are based on the recurrence parameters and intensity attenuation function for the region. The results place Adelaide close to the AS2121 ‐ 1979 Earthquake Code Zone I/Zone 2 boundary.     

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.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Seismic response of high concrete face rockfill dams subjected to non-uniform input motion

Analysis of the seismic response of high CFRDs under non-uniform ground motion input is conducted using a novel non-uniform input motion calculation method combined with nonlinear FEM. The non-uniform input motion calculation method and its basic assumption are validated. The response of CFRDs under uniform and non-uniform input is compared to discuss the necessity of conducting seismic analysis of high CFRDs under realistic non-uniform ground motion input. When the acceleration at the surface of the free field for dynamic simulations with uniform and non-uniform input is kept consistent, the seismic response of CFRDs under non-uniform input is in general significantly smaller, while the dynamic tensile stress around the edges of the concrete face slab is greater. The simulation results suggest that non-uniformity of the ground motion input has important effects on the seismic response of high CFRDs and should be considered in the seismic design of CFRDs. The influence of the incident angle of seismic waves is also investigated, with results indicating that the influence is waveform dependent, while being frequency independent.

     

Strong historical earthquakes in the northwestern Issyk Kul’ basin (northern Tien Shan)

The northern Tien Shan is the northern front of the Himalayan mountain belt, which resulted from the collision between the Indian and Eurasian Plates. This region encompasses the most active seismic zones of the orogen, which generated the strongest (M > 8) earthquakes. Since there are scarcely any written accounts, the only way to trace back strong earthquakes is the paleoseismologic method. Since 1984 we have been studying the northwestern Issyk Kul’ basin, where there are differently directed anticlines, which constitute the Kungei meganticline. Here, several active tectonic structures (faults, folds) are located, whose development was accompanied by strong earthquakes. Our field studies of 2008 in the Iiri-Taldybulak Valley, along the adyrs (foothills) of the Kungei-Ala-Too Range, revealed two unknown historical earthquakes. The first one, which occurred along the southern rupture in the late 7th century A.D., gave rise to a seismic scarp; the latter broke through the river floodplain and a tash-koro (ancient settlement). The second one, which occurred along the northern rupture in the late 9th century A.D., increased the height of the seismic scarp, existing on the Early Holocene and older terraces. Note that this region already records a strong seismic event around 500 A.D. Archeologic data have revealed one more strong earthquake, which took place in the 14th century A.D. Note that the above-mentioned strong seismic events are coeval with the decline of the nomadic cultures (Wusun, Turkic, Mogul) in the northern Tien Shan and Zhetysu (Semirech’e).     

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.     

Modelling Geophysical Precursors to the Prehistoric c. AD1305 Kaharoa Rhyolite Eruption of Tarawera Volcano, New Zealand

The 1305 Kaharoa rhyolite eruptive episode is the largest volcanic event(4 km³ magma) to have occurred in New Zealand during the last 1000 years. Proximal areas were devastated by pyroclastic flows, and tephra fell over much of the northern North Island. No eyewitness observations are recorded, but ejecta analyses show that the rhyolite eruptions were primed and triggered by basalt intrusions. This key finding, combined with observations of similar modern eruptions, has allowed construction of a conceptual scenario of the seismic and other activity that likely preceded the Kaharoa episode.The precursory scenario begins at -5 years (before the first eruption). Rising basalt magma intrusions generate deep long-period earthquakes in the lower crust, before intersecting and heating a rhyolite magma body at 6 km depth beneath Tarawera. By -1 year, increased heat flux from the rhyolite magma body had raised temperatures and pressures in the overlying hydrothermal system; generating shallow long-period earthquakes and increased heat flow at the surface. At -2 months, shallow volcano-tectonic earthquake activity intensified, driven by inflation of the rhyolite magma body, with magmatic gas appearing in fumarole discharges. Rapidly accelerating seismicity, ground deformation and surface heat flow occurred in the last few weeks and days, before the initial vent-opening explosions intensified into major plinian eruptions.Effectiveness of the present volcano monitoring system at Tarawera can be evaluated against this scenario. The precursory seismic activity, including the critical deep long-period earthquakes, would be recorded but not accurately located. Similarly, the existing ground deformation monitoring systems would detect early magma chamber inflation, but discrimination from the background tectonic tilting signal would be difficult. Continuous telemetering of geodetic data from existing and additional instruments would be required for any useful monitoring of rapid ground deformation in the final precursory phases.