Spatial variation of seismicity parameters across India and adjoining areas

An attempt has been made to quantify the variability in the seismic activity rate across the whole of India and adjoining areas (0–45°N and 60–105°E) using earthquake database compiled from various sources. Both historical and instrumental data were compiled and the complete catalog of Indian earthquakes till 2010 has been prepared. Region-specific earthquake magnitude scaling relations correlating different magnitude scales were achieved to develop a homogenous earthquake catalog for the region in unified moment magnitude scale. The dependent events (75.3%) in the raw catalog have been removed and the effect of aftershocks on the variation of b value has been quantified. The study area was divided into 2,025 grid points (1°91°) and the spatial variation of the seismicity across the region have been analyzed considering all the events within 300 km radius from each grid point. A significant decrease in seismic b value was seen when declustered catalog was used which illustrates that a larger proportion of dependent events in the earthquake catalog are related to lower magnitude events. A list of 203,448 earth- quakes (including aftershocks and foreshocks) occurred in the region covering the period from 250 B.C. to 2010 A.D. with all available details is uploaded in the website http://www.civil.iisc.ernet.in/*sreevals/resource.htm.

Maximum magnitude estimation considering the regional rupture character

The main objective of the paper is to develop a new method to estimate the maximum magnitude (M (max)) considering the regional rupture character. The proposed method has been explained in detail and examined for both intraplate and active regions. Seismotectonic data has been collected for both the regions, and seismic study area (SSA) map was generated for radii of 150, 300, and 500 km. The regional rupture character was established by considering percentage fault rupture (PFR), which is the ratio of subsurface rupture length (RLD) to total fault length (TFL). PFR is used to arrive RLD and is further used for the estimation of maximum magnitude for each seismic source. Maximum magnitude for both the regions was estimated and compared with the existing methods for determining M (max) values. The proposed method gives similar M (max) value irrespective of SSA radius and seismicity. Further seismicity parameters such as magnitude of completeness (M (c) ), ``a'' and ``aEuro parts per thousand b `` parameters and maximum observed magnitude (M (max) (obs) ) were determined for each SSA and used to estimate M (max) by considering all the existing methods. It is observed from the study that existing deterministic and probabilistic M (max) estimation methods are sensitive to SSA radius, M (c) , a and b parameters and M (max) (obs) values. However, M (max) determined from the proposed method is a function of rupture character instead of the seismicity parameters. It was also observed that intraplate region has less PFR when compared to active seismic region.

Seismic characterization and dynamic site response of a municipal solid waste landfill in Bangalore, India

Seismic design of landfills requires an understanding of the dynamic properties of municipal solid waste (MSW) and the dynamic site response of landfill waste during seismic events. The dynamic response of the Mavallipura landfill situated in Bangalore, India, is investigated using field measurements, laboratory studies and recorded ground motions from the intraplate region. The dynamic shear modulus values for the MSW were established on the basis of field measurements of shear wave velocities. Cyclic triaxial testing was performed on reconstituted MSW samples and the shear modulus reduction and damping characteristics of MSW were studied. Ten ground motions were selected based on regional seismicity and site response parameters have been obtained considering one-dimensional non-linear analysis in the DEEPSOIL program. The surface spectral response varied from 0.6 to 2g and persisted only for a period of 1s for most of the ground motions. The maximum peak ground acceleration (PGA) obtained was 0.5g and the minimum and maximum amplifications are 1.35 and 4.05. Amplification of the base acceleration was observed at the top surface of the landfill underlined by a composite soil layer and bedrock for all ground motions. Dynamic seismic properties with amplification and site response parameters for MSW landfill in Bangalore, India, are presented in this paper. This study shows that MSW has less shear stiffness and more amplification due to loose filling and damping, which need to be accounted for seismic design of MSW landfills in India.

Seismic characterization and dynamic site response of a municipal solid waste landfill in Bangalore, India

Seismic design of landfills requires an understanding of the dynamic properties of municipal solid waste (MSW) and the dynamic site response of landfill waste during seismic events. The dynamic response of the Mavallipura landfill situated in Bangalore, India, is investigated using field measurements, laboratory studies and recorded ground motions from the intraplate region. The dynamic shear modulus values for the MSW were established on the basis of field measurements of shear wave velocities. Cyclic triaxial testing was performed on reconstituted MSW samples and the shear modulus reduction and damping characteristics of MSW were studied. Ten ground motions were selected based on regional seismicity and site response parameters have been obtained considering one-dimensional non-linear analysis in the DEEPSOIL program. The surface spectral response varied from 0.6 to 2g and persisted only for a period of 1s for most of the ground motions. The maximum peak ground acceleration (PGA) obtained was 0.5g and the minimum and maximum amplifications are 1.35 and 4.05. Amplification of the base acceleration was observed at the top surface of the landfill underlined by a composite soil layer and bedrock for all ground motions. Dynamic seismic properties with amplification and site response parameters for MSW landfill in Bangalore, India, are presented in this paper. This study shows that MSW has less shear stiffness and more amplification due to loose filling and damping, which need to be accounted for seismic design of MSW landfills in India.

The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data

The stress release model, a stochastic version of the elastic-rebound theory, is applied to the historical earthquake data from three strong earthquake-prone regions of China, including North China, Southwest China, and the Taiwan seismic regions. The results show that the seismicity along a plate boundary (Taiwan) is more active than in intraplate regions (North and Southwest China). The degree of predictability or regularity of seismic events in these seismic regions, based on both the Akaike information criterion (AIC) and fitted sensitivity parameters, follows the order Taiwan, Southwest China, and North China, which is further identified by numerical simulations. (c) 2004 Elsevier Ltd. All rights reserved.

Oceanic intraplate volcanoes exposed: Example from seamounts accreted in Panama

Two Paleogene ocean islands are exposed in the Azuero Peninsula, west Panama, within sequences accreted in the early-Middle Eocene. A multidisciplinary approach involving litho-logic mapping, paleontological age determinations, and petrological study allows reconstruction of the stratigraphy and magmatic evolution of one of these intraplate oceanic volcanoes. From base to top, the volcano's structure comprises submarine basaltic lava flows locally interlayered with hemipelagic sediments, basaltic breccias, shallow-water limestones, and subaerial basaltic lava. Gabbros and basaltic dikes were emplaced along a rift zone of the island. Geochemical trends of basaltic lavas include decreased Mg# {[Mg/(Mg + Fe)] * 100} and, with time, increased incompatible element contents thought to be representative of many poorly documented intraplate volcanoes in the Pacific. Our results show that, in addition to deep drilling, the roots of oceanic islands can be explored through studies of accreted and subaerially exhumed oceanic sequences.

Characteristics of near surface crustal deformation associated with shield seismicity

The purpose of the present study is to understand the surface deformation associated with the Killari and Wadakkancheri earthquake and to examine if there are any evidence of occurrence of paleo-earthquakes in this region or its vicinity. The present study is an attempt to characterize active tectonic structures from two areas within penisular India: the sites of 1993 Killari (Latur) (Mb 6.3) and 1994 Wadakkancheri (M 4.3) earthquakes in the Precambrian shield. The main objectives of the study are to isolate structures related to active tectonism, constraint the style of near – surface deformation and identify previous events by interpreting the deformational features. The study indicates the existence of a NW-SE trending pre-existing fault, passing through the epicentral area of the 1993 Killari earthquake. It presents the salient features obtained during the field investigations in and around the rupture zone. Details of mapping of the scrap, trenching, and shallow drilling are discussed here. It presents the geologic and tectonic settings of the Wadakkancheri area and the local seismicity; interpretation of remote sensing data and a detailed geomorphic analysis. Quantitative geomorphic analysis around the epicenter of the Wadakkancheri earthquake indicates suitable neotectonic rejuvenation. Evaluation of remote sensing data shows distinct linear features including the presence of potentially active WNW-ESE trending fault within the Precambrian shear zone. The study concludes that the earthquakes in the shield area are mostly associated with discrete faults that are developed in association with the preexisting shear zones or structurally weak zones

Seismic behavior of capacity designed masonry walls in low seismicity regions

Eurocode 8 representing a new generation of structural design codes in Europe defines ‎requirements for the design of buildings against earthquake action. In Central and ‎Western Europe, the newly defined earthquake zones and corresponding design ground ‎acceleration values, will lead in many cases to earthquake actions which are remarkably ‎higher than those defined so far by the design codes used until now in Central Europe. ‎ In many cases, the weak points of masonry structures during an earthquake are the corner ‎regions of the walls. Loading of masonry walls by earthquake action leads in most cases ‎to high shear forces. The corresponding bending moment in such a wall typically causes a ‎significant increase of the eccentricity of the normal force in the critical wall cross ‎section. This in turn leads ultimately to a reduction of the size of the compression zone in ‎unreinforced walls and a high concentration of normal stresses and shear stresses in the ‎corner regions. ‎ Corner-Gap-Elements, consisting of a bearing beam located underneath the wall and ‎made of a sufficiently strong material (such as reinforced concrete), reduce the effect of ‎the eccentricity of the normal force and thus restricts the pinching effect of the ‎compression zone. In fact, the deformation can be concentrated in the joint below the ‎bearing beam. According to the principles of the Capacity Design philosophy, the ‎masonry itself is protected from high stresses as a potential cause of brittle failure. ‎ Shaking table tests at the NTU Athens Earthquake Engineering Laboratory have proven ‎the effectiveness of the Corner-Gap-Element. The following presentation will cover the ‎evaluation of various experimental results as well as a numerical modeling of the ‎observed phenomena.‎

Induced seismicity in the Castarthao reservoir, NE Brazil - Preliminary results

The Castanhao reservoir was built in the state of Ceara, a dry region in Northeastern Brazil, to regulate the flow of the Jaguaribe River, for irrigation, and for power generation. It is an earth-filled dam, 60 m high, with a water capacity of 4.5 x 10(9) m(3). The seismicity in the area has been monitored since 1998, with a few interruptions, using one analog or one digital station and, during a few periods, a three-station network. The first earthquakes likely to be induced events were detected in 2003, when the water level was about 20 in high. In early 2004 a very heavy rainfall season quickly filled the reservoir. Shortly after, an increase in the seismic activity occurred and many micro-earthquakes were recorded. We suggest that this activity resulted from an increase in pore pressure due to undrained response. Therefore, we may classify this cluster of microearthquakes as ""initial seismicity."" We deployed a network with four analog stations in the area, following this activity, to determine the epicentral zone. At least three epicentral areas under the reservoir were detected. The spatio-temporal analysis of the available data revealed that the seismicity occurs in clusters and that these were activated at different periods. We identified four sets of faults (N-S-, E-W-, NW-SE-, and NE-SW-oriented), some of which moved in shallow crustal levels and as recently as the Quaternary (1.8 Ma). Under the present-day stress regime, the last two sets moved as strike-slip structures. We suggest a possible correlation between dormant faults and the observed induced seismicity. (c) 2008 Elsevier B.V. All rights reserved.

Early Permian post-glacial bivalve faunas of the Itarare Group, Parana Basin, Brazil: Paleoecology and biocorrelations with South American intraplate basins

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Strain rates estimated by geodetic observations in the Borborema Province, Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

TDEM survey in an area of seismicity induced by water wells in Parana sedimentary basin, Northern Sao Paulo State, Brazil

This article presents TDEM results from an area with recent induced shallow seismicity. The purpose was to do a geoelectrical mapping of sedimentary and fractured basaltic aquifers for better understanding of the hydrogeologic setting. The study area is in the Parana basin where flood basalts are overlain by sedimentary units near the city of Bebedouro, northern Sao Paulo State, Brazil. 86 TDEM soundings were acquired in an area of 90 km(2) in the Andes and Botafogo study areas. The soundings were chosen next to wells for calibration, and also along profiles crossing the seismically active areas. 1-D interpretation results showed the general geoelectrical stratigraphy of this part of the Parana basin. The upper geoelectrical layer is the shallow sedimentary aquifer (Adamantina formation) with less than 80 m thickness. The second geoelectrical layer contains the upper basalts of the Serra Geral formation at about 60-80 m depths. A saturated fractured basalt zone between 100 and 300 m depths was identifiable on various TDEM soundings. This depth range corresponds to the range of hypocentral depths for more than 3000 micro-earthquakes in this area. The lower basalt layer was estimated to lie between 400 and 650 m depth. The deepest geoelectrical layer detected by various TDEM soundings corresponds to the Botucatu sandstone (Guarani aquifer). Results suggest that the high-discharge wells are located in the fractured zone in the middle basalt of the Serra Geral formation. There is a good correlation between seismically active areas, high discharge wells (>190 m(3)/h), and fracture zones in the middle basalt. The results reinforce the hypothesis that the shallow seismic activity in the Bebedouro region is being triggered by high rates of groundwater withdrawal. (C) 2012 Elsevier B.V. All rights reserved.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.