Geothermal investigations in an area of induced seismic activity, northern São Paulo State, Brazil

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

Seismic activity triggered by water wells in the Parana Basin, Brazil

Triggered seismicity is commonly associated with deep water reservoirs or injection wells where water is injected at high pressure into the reservoir rock. However, earth tremors related solely to the opening of groundwater wells are extremely rare. Here we present a clear case of seismicity induced by pore-pressure changes following the drilling of water wells that exploit a confined aquifer in the intracratonic Parana Basin of southeastern Brazil. Since 2004, shallow seismic activity, with magnitudes up to 2.9 and intensities V MM, has been observed near deep wells (120-200 m) that were drilled in early 2003 near the town of Bebedouro. The wells were drilled for irrigation purposes, cross a sandstone layer about 60-80 m thick and extract water from a confined aquifer in fractured zones between basalt flow layers. Seismic activity, mainly event swarms, has occurred yearly since 2004, mostly during the rainy season when the wells are not pumped. During the dry season when the wells are pumped almost continuously, the activity is very low. A seismographic network, installed in March 2005, has located more than 2000 microearthquakes. The events are less than 1 km deep (mostly within the 0.5 km thick basalt layer) and cover an area roughly 1.5 km x 5 km across. The seismicity generally starts in a small area and expands to larger distances with an equivalent hydraulic diffusivity ranging from 0.06 to 0.6 m(2)/s. Geophysical and geothermal logging of several wells in the area showed that water from the shallow sandstone aquifer enters the well at the top and usually forms waterfalls. The waterfalls flow down the sides of the wells and feed the confined, fractured aquifer in the basalt layer at the bottom. Two seismic areas are observed: the main area surrounds several wells that are pumped continuously during the dry season, and a second area near another well (about 10 km from the first area) that is not used for irrigation and not pumped regularly. The main area displays cyclic annual activity, but the second area does not. We explain the earthquake swarms as being triggered by pore pressure diffusion in the fractured basalt layer due to additional pressure from the newly connected surface aquifer. This reaches critically prestressed areas up to a few kilometers away from the wells. During periods of continuous pumping, the reduction of pore pressure in the confined aquifer stops the seismic activity. Our study suggests that this kind of activity may be more common than previously thought and implies that many other cases of small tremors associated with the drilling of water wells may have gone unnoticed.

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.

Seismic Activity Rates in the Iberian Peninsula

Evaluating the seismic hazard requires establishing a distribution of the seismic activity rate, irrespective of the methodology used in the evaluation. In practice, how that activity rate is established tends to be the main difference between the various evaluation methods. The traditional procedure relies on a seismogenic zonation and the Gutenberg-Richter (GR) hypothesis. Competing zonations are often compared looking only at the geometry of the zones, but the resulting activity rate is affected by both geometry and the values assigned to the GR parameters. Contour plots can be used for conducting more meaningful comparisons, providing the GR parameters are suitably normalised. More recent approaches for establishing the seismic activity rate forego the use of zones and GR statistics and special attention is paid here to such procedures. The paper presents comparisons between the local activity rates that result for the complete Iberian Peninsula using kernel estimators as well as two seismogenic zonations. It is concluded that the smooth variation of the seismic activity rate produced by zoneless methods is more realistic than the stepwise changes associated with zoned approaches; moreover, the choice of zonation often has a stronger influence on the results than its fairly subjective origin would warrant. It is also observed that the activity rate derived from the kernel approach, related with the GR parameter “a”, is qualitatively consistent with the epicentres in the catalogue. Finally, when comparing alternative zonations it is not just their geometry but the distribution of activity rate that should be compared.

An approach to palaeoseismicity in the Olkiluoto (sea) area during the early Holocene

Olkiluoto Island is situated in the northern Baltic Sea, near the southwestern coast of Finland, and is the proposed location of a spent nuclear fuel repository. This study examined Holocene palaeoseismicity in the Olkiluoto area and in the surrounding sea areas by computer simulations together with acoustic-seismic, sedimentological and dating methods. The most abundant rock type on the island is migmatic mica gneiss, intruded by tonalites, granodiorites and granites. The surrounding Baltic Sea seabed consists of Palaeoproterozoic crystalline bedrock, which is to a great extent covered by younger Mesoproterozoic sedimentary rocks. The area contains several ancient deep-seated fracture zones that divide it into bedrock blocks. The response of bedrock at the Olkiluoto site was modelled considering four future ice-age scenarios. Each scenario produced shear displacements of fractures with different times of occurrence and varying recovery rates. Generally, the larger the maximum ice load, the larger were the permanent shear displacements. For a basic case, the maximum shear displacements were a few centimetres at the proposed nuclear waste repository level, at proximately 500 m b.s.l. High-resolution, low-frequency echo-sounding was used to examine the Holocene submarine sedimentary structures and possible direct and indirect indicators of palaeoseismic activity in the northern Baltic Sea. Echo-sounding profiles of Holocene submarine sediments revealed slides and slumps, normal faults, debris flows and turbidite-type structures. The profiles also showed pockmarks and other structures related to gas or groundwater seepages, which might be related to fracture zone activation. Evidence of postglacial reactivation in the study area was derived from the spatial occurrence of some of the structures, especial the faults and the seepages, in the vicinity of some old bedrock fracture zones. Palaeoseismic event(s) (a single or several events) in the Olkiluoto area were dated and the palaeoenvironment was characterized using palaeomagnetic, biostratigraphical and lithostratigraphical methods, enhancing the reliability of the chronology. Combined lithostratigraphy, biostratigraphy and palaeomagnetic stratigraphy revealed an age estimation of 10 650 to 10 200 cal. years BP for the palaeoseismic event(s). All Holocene sediment faults in the northern Baltic Sea occur at the same stratigraphical level, the age of which is estimated at 10 700 cal. years BP (9500 radiocarbon years BP). Their movement is suggested to have been triggered by palaeoseismic event(s) when the Late Weichselian ice sheet was retreating from the site and bedrock stresses were released along the bedrock fracture zones. Since no younger or repeated traces of seismic events were found, it corroborates the suggestion that the major seismic activity occurred within a short time during and after the last deglaciation. The origin of the gas/groundwater seepages remains unclear. Their reflections in the echo-sounding profiles imply that part of the gas is derived from the organic-bearing Litorina and modern gyttja clays. However, at least some of the gas is derived from the bedrock. Additional information could be gained by pore water analysis from the pockmarks. Information on postglacial fault activation and possible gas and/or fluid discharges under high hydraulic heads has relevance in evaluating the safety assessment of a planned spent nuclear fuel repository in the region.

A study on seismicity and seismic hazard for Karnataka State

This paper presents a detailed study on the seismic pattern of the state of Karnataka and also quantifies the seismic hazard for the entire state. In the present work, historical and instrumental seismicity data for Karnataka (within 300 km from Karnataka political boundary) were compiled and hazard analysis was done based on this data. Geographically, Karnataka forms a part of peninsular India which is tectonically identified as an intraplate region of Indian plate. Due to the convergent movement of the Indian plate with the Eurasian plate, movements are occurring along major intraplate faults resulting in seismic activity of the region and hence the hazard assessment of this region is very important. Apart from referring to seismotectonic atlas for identifying faults and fractures, major lineaments in the study area were also mapped using satellite data. The earthquake events reported by various national and international agencies were collected until 2009. Declustering of earthquake events was done to remove foreshocks and aftershocks. Seismic hazard analysis was done for the state of Karnataka using both deterministic and probabilistic approaches incorporating logic tree methodology. The peak ground acceleration (PGA) at rock level was evaluated for the entire state considering a grid size of 0.05A degrees x 0.05A degrees. The attenuation relations proposed for stable continental shield region were used in evaluating the seismic hazard with appropriate weightage factors. Response spectra at rock level for important Tier II cities and Bangalore were evaluated. The contour maps showing the spatial variation of PGA values at bedrock are presented in this work.

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 degrees 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.

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.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Development of open models and tools for seismic risk assessment: application to Portugal

The exponential growth of the world population has led to an increase of settlements often located in areas prone to natural disasters, including earthquakes. Consequently, despite the important advances in the field of natural catastrophes modelling and risk mitigation actions, the overall human losses have continued to increase and unprecedented economic losses have been registered. In the research work presented herein, various areas of earthquake engineering and seismology are thoroughly investigated, and a case study application for mainland Portugal is performed. Seismic risk assessment is a critical link in the reduction of casualties and damages due to earthquakes. Recognition of this relation has led to a rapid rise in demand for accurate, reliable and flexible numerical tools and software. In the present work, an open-source platform for seismic hazard and risk assessment is developed. This software is capable of computing the distribution of losses or damage for an earthquake scenario (deterministic event-based) or earthquake losses due to all the possible seismic events that might occur within a region for a given interval of time (probabilistic event-based). This effort has been developed following an open and transparent philosophy and therefore, it is available to any individual or institution. The estimation of the seismic risk depends mainly on three components: seismic hazard, exposure and vulnerability. The latter component assumes special importance, as by intervening with appropriate retrofitting solutions, it may be possible to decrease directly the seismic risk. The employment of analytical methodologies is fundamental in the assessment of structural vulnerability, particularly in regions where post-earthquake building damage might not be available. Several common methodologies are investigated, and conclusions are yielded regarding the method that can provide an optimal balance between accuracy and computational effort. In addition, a simplified approach based on the displacement-based earthquake loss assessment (DBELA) is proposed, which allows for the rapid estimation of fragility curves, considering a wide spectrum of uncertainties. A novel vulnerability model for the reinforced concrete building stock in Portugal is proposed in this work, using statistical information collected from hundreds of real buildings. An analytical approach based on nonlinear time history analysis is adopted and the impact of a set of key parameters investigated, including the damage state criteria and the chosen intensity measure type. A comprehensive review of previous studies that contributed to the understanding of the seismic hazard and risk for Portugal is presented. An existing seismic source model was employed with recently proposed attenuation models to calculate probabilistic seismic hazard throughout the territory. The latter results are combined with information from the 2011 Building Census and the aforementioned vulnerability model to estimate economic loss maps for a return period of 475 years. These losses are disaggregated across the different building typologies and conclusions are yielded regarding the type of construction more vulnerable to seismic activity.

Análise estática não linear de um edifício em betão armado

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil

Investigating Environmental Determinants of Diet, Physical Activity, and Overweight among Adults in Sao Paulo, Brazil

There is worldwide recognition that the burden of noncommunicable diseases (NCDs) and obesity-related health problems is rapidly increasing in low- and middle-income countries. Environmental determinants of obesity are likely to differ between countries, particularly in those undergoing rapid socioeconomic and nutrition transitions such as Brazil. This study aims to describe some built environment and local food environment variables and to explore their association with the overweight rate and diet and physical activity area-level aggregated indicators of adults living in the city of Sao Paulo, the largest city in Brazil. This formative study includes an ecological analysis of environmental factors associated with overweight across 31 submunicipalities of the city of Sao Paulo using statistical and spatial analyses. Average prevalence of overweight was 41.69% (95% confidence interval 38.74, 44.64), ranging from 27.14% to 60.75% across the submunicipalities. There was a wide geographical variation of both individual diet and physical activity, and indicators of food and built environments, favoring wealthier areas. After controlling for area socioeconomic status, there was a positive correlation between regular fruits and vegetables (FV) intake and density of FV specialized food markets (r = 0.497; p < 0.001), but no relationship between fast-food restaurant density and overweight prevalence was found. A negative association between overweight prevalence and density of parks and public sport facilities was seen (r = -0.527; p < 0.05). Understanding the relationship between local neighborhood environments and increasing rates of poor diet, physical activity, and obesity is essential in countries undergoing rapid economic and urban development, such as Brazil, in order to provide insights for policies to reduce increasing rates of NCDs and food access and health inequalities.

Intraplate earthquake swarm in Belo Jardim, NE Brazil: reactivation of a major Neoproterozoic shear zone (Pernambuco Lineament)

Intraplate earthquakes in stable continental areas have been explained basically by reactivation of pre-existing zones of weakness, stress concentration, or both. Zones of weakness are usually identified as sites of the last major orogeny, provinces of recent alkaline intrusions, or stretched crust in ancient rifts. However, it is difficult to identify specific zones of weakness and intraplate fault zones are not always easily correlated with known geological features. Although Northeastern Brazil is one of the most seismically active areas in the country (magnitudes 5 roughly every 5 yr), with hypocentral depths shallower than similar to 10 km and seismic zones as long as 30-40 km, no clear relationship with the known surface geology can be usually established with confidence, and a clear identification of zones of weakness has not yet been possible. Here we present the first clear case of seismic activity occurring as reactivation of an old structure in Brazil: the Pernambuco Lineament, a major Neoproterozoic shear zone. The 2004 earthquake swarm of Belo Jardim (magnitudes up to 3.1) and the recurrent activities in the nearby towns of Sao Caetano and Caruaru (magnitudes up to 4.0 and 3.8), show that the Pernambuco Lineament is a weak zone. A local seismic network showed that the Belo Jardim swarm of 2004 November occurred by normal faulting on a North dipping, E-W oriented fault plane in close agreement with the E-W trending structures within the Pernambuco Lineament. The Belo Jardim activity was concentrated in a 1.5 km (E-W) by 2 km (downdip) fault area, and average depth of 4.5 km. The nearby Caruaru activity occurs as both strike-slip and normal faulting, also consistent with local structures of the Pernambuco Lineament. The focal mechanisms of Belo Jardim, Caruaru and S. Caetano, indicate E-W compressional and N-S extensional principal stresses. The NS extension of this stress field is larger than that predicted by numerical models such as those of Coblentz & Richardson and we propose that additional factors such as flexural stresses from the nearby Sergipe-Alagoas marginal basin could also affect the current stress field in the Pernambuco Lineament.