Tectônica rúptil e sismicidade na área de inundação do açude do Castamhão (CE): implicações para o risco sísmico

This study describes brittle deformation and seismicity in the Castanhão Dam region, Ceará State, Brazil. This reservoir will include a hidroeletric power plant and will store about 6,7 billions m3 of water. Five main litostratigraphic unit were identified in the region: gneissic-migmatitic basement, metavolcanosedimentary sequence, granitoid plutons of Brasiliano age, Mesozoic basaltic dike swarm, and Cenozoic fluvial terraces of the Jaguaribe river. The region has experienced several faulting events that occurred at different crustal levels. Faults formed at depths less than about 12 km present left-lateral movement and are associated with epidote and quartz infillings. Faults formed at depths less than 7 km are mainly strike-slip present cataclastic rocks, fault breccia and gouge. Both fault groups form mainly NE-trendind lineaments and represent reactivation of ductile shear zones or new formed faults that cut across existing structures. Seismically-induced liquefaction fractures take place in Cenozoic terraces and indicate paleoearthquakes that may have reached at leat 6,8 MS. In short, this work indicate that the level of paleoseismicity is much greater than one observed in the instrumental record. Several faults are favourably oriented for reactivation and induced seismicity should be expected after the Castanhão Dam impoudment

Estudo do potencial da técnica de nanossísmica para o monitoramento de hidrofraturamento em reservatórios de hidrocarbonetos

Hydraulic fracturing is an operation in which pressurised fluid is injected in the geological formation surrounding the producing well to create new permeable paths for hydrocarbons. The injection of such fluids in the reservoir induces seismic events. The measurement of this reservoir stimulation can be made by location these induced microseismic events. However, microseismic monitoring is an expensive operation because the acquisition and data interpretation system using in this monitoring rely on high signal-to-noise ratios (SNR). In general, the sensors are deployed in a monitoring well near the treated well and can make a microseismic monitoring quite an expensive operation. In this dissertation we propose the application of a new method for recording and location of microseismic events called nanoseismic monitoring (Joswig, 2006). In this new method, a continuous recording is performed and the interpreter can separate events from noise using sonograms. This new method also allows the location of seismic sources even when P and S phases onsets are not clear like in situations of 0 dB SNR. The clear technical advantage of this new method is also economically advantageous since the sensors can potentially be installed on the surface rather than in observation well. In this dissertation field tests with controlled sources were made. In the first test small explosives using fire works at 28 m (slant distances) were detected yealding magnitudes between -2.4 ≤ ML ≤ -1.6.. In a second test, we monitored perforation shots in a producing oil field. In this second test, one perforation shot was located with slant distances of 861 m and magnitude 2.4 ML. Data from the tests allow us to say that the method has potential to be used in the oil industry to monitor hydrofracture

Estudoda sismicidade na região de Sobral - CE em 2008

The town of Sobral is located at the northwest part of the Ceará State, 250 km away from its capital, Fortaleza. In January 2008, an intense seismic activity began near Sobral with one event with magnitude 4,2mb on May 21. Since the start of its seismic activity, all events were recorded by the SBBR station (located on EMBRAPA Caprinos Farm), which operates in the region since August 2007. After this event, monitoring the seismic activity was carried out with the deployment of a local three component digital seismographic network, from June 5, 2008 until September 24. Initially, this network was composed of six seismographic stations. Later additional five stations were deployed until August 26 2008. This local network detected approximately 2,800 earthquakes. In this study we analyzed 581 earthquakes recorded by at least three stations for hypocentral and focal mechanism determination, and to contribute to a better explanation of the seismicity which in this region. To determine the hypocenters, we used a half-space model, with vP = 6,00 km/s and vP/vS = 1,71. From the hypocentral determination, it was revealed an active seismic zone with depth ranging between 1 and 8 km, 6 km long in E - W direction. The determination of fault planes and focal mechanism was obtained using the programs FPFIT and PLAN, which allowed comparison between their respective results in order to obtain more accurate results. A set of 24 earthquakes were selected to determine fault using PLAN planes and focal mechanisms using FPFIT. With the aid of detailed map of hypocenters this set, it was possible to identify three structures. Therefore, the set of 24 earthquakes were divided into three subsets. The type of mechanism was predominantly strike-slip with a dextral direction. Although the region has two tectonic structures near the site of the study area: the Café- Ipueiras Fault (normal fault) and the Sobral-Pedro II Lineament (dextral strike-slip fault) it was not possible to correlate the seismicity founded with those structures

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)

Diques clásticos na formação corumbataí (p) no nordeste da bacia do paraná, SP: Análise sistemática e significações estratigráficas, sedimentológicas e tectônicas

This paper aims to present the results of systematic survey on clastic dykes in the Corumbataí Formation (Permian), in the northeast region of the State of São Paulo. Besides this, the paper analyses genetic aspects of those features as well as their stratigraphic and sedimentologic implications in terms of geologic evolution of the northeastern Paraná Sedimentary Basin during Permian times. The field works had been developed in 3 main Corumbataí Formation outcrops (2 road cuts and a quarry) supposed to show the most important clastic dikes occurrences in the studied area. Basically, the sedimentary intrusions are formed by fine sand or silt size particles and had penetrated host rocks as near-vertical, centimeter thick, dykes (most common form) or as horizontal sheets, forming clastic sills (subsidiary form), both with variable geometric forms and dimensions. A lot of dyke walls show undulations suggesting pre-diagenetic clastic intrusions, probably near the ancient depositional surface. Almost all intrusions occur in the superior third portion of the Corumbataí Formation and some similar features seem to appear in the adjoining superposed Pirambóia Formation base. In this article the authors defend a seismic origin hypothesis for the clastic intrusions. It is important to mention that clastic intrusions tend to occur linked to expressive seismic events, with magnitude upper to 5. The analysis of isopach maps of the Permian and Mesozoic units of the Paraná Sedimentary Basin in the study area suggests a depositional system changing, from epicontinental sea conditions to shallow platform and, finally, to coastal deserts. Probably, this environmental change was driven by regional uplift accompanied by seismic events. It is possible that ancient seismicity triggered liquefaction processes and the resulting clastic intrusions. In this sense, those clastic features might be properly namedseismites.

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

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

Seasonal Variation in Reception of Fin Whale Calls at Five Geographic Areas in the North Pacific

In late August 1991 scientists at the National Oceanic and Atmospheric Administration’s (NOAA) National Marine Mammal Laboratory (NMML) and Pacific Marine Environmental Laboratory (PMEL) began a pilot study to investigate the capability of hydrophones from the US. Navy’s fixed array system to detect large whales in the North Pacific by passive reception of their calls. PMEL had previously established a direct data link from five bottom-mounted arrays of the Navy SOSUS (Sound Surveillance System), via the Naval Oceanographic Processing Facility (NOPF) at Whidbey Island, Washington, to study low-level seafloor seismicity (Fox et al. 1994). PMEL subsequently provided NMML tapes of SOSUS hydrophone data from which whale calls were analyzed. As in an analogous study conducted in the North Atlantic (Nishimura and Conlon 1994, Clark 1995, Mellinger and Clark 1995), calls attributable to whales were received at each SOSUS site at rates that varied seasonally (Anonymous 1996).

Upper crustal earthquake swarms in São Caetano: reactivation of the Pernambuco shear zone and trending branches in intraplate Brazil

Seismogenic fault reactivation of continental-scale structures has been observed in a few intraplate areas, but its cause is still amatter of debate. The objective of the present study is to analyze two seismic swarms that occurred along the EW-trending Pernambuco ductile shear zone and in a NE-trending branch, in 2007 and 2010 in São Caetano County, Northeastern Brazil.We studied both epicentral areas using a nine- and a seven-station network during 180 and 54 days, respectively. The results indicate that the 2007 swarm correspond to a right-lateral, strike–slip fault with a normal component of slip (strike 74°, dip 60°, and rake−145°) and the 2010 swarmcorresponds to a normal fault (strike 265°, dip 79°, and rake −91°). The former reactivated a NE-trending branch, whereas the latter reactivated the main E-W-trending mylonitic belt of the Pernambuco shear zone. These results are consistent with seismogenic reactivation of this major structure, generated by the present-day EW-trending compression and NS-trending extension, as observed by previous studies. This shear zone was reactivated as rift faults in the Cretaceous during the South America–Africa breakup. However, our study confirms that the basement fabric such as continental-scale ductile shear zones, show evidence of crustal weakness outside areas of previous rifting, and it reveals the potential for large earthquakes along dormant rift segments associated with major basement shear belts.

Testing the influence of far-field topographic forcing on subduction initiation at a passive margin

Despite favourable gravitational instability and ridge-push, elastic and frictional forces prevent subduction initiation fromarising spontaneously at passive margins. Here,we argue that forces arising fromlarge continental topographic gradients are required to initiate subduction at passivemargins. In order to test this hypothesis,we use 2Dnumerical models to assess the influence of the Andean Plateau on stressmagnitudes and deformation patterns at the Brazilian passive margin. The numerical results indicate that “plateau-push” in this region is a necessary additional force to initiate subduction. As the SE Brazilianmargin currently shows no signs of self-sustained subduction, we examined geological and geophysical data to determine if themargin is in the preliminary stages of subduction initiation. The compiled data indicate that the margin is presently undergoing tectonic inversion, which we infer as part of the continental–oceanic overthrusting stage of subduction initiation. We refer to this early subduction stage as the “Brazilian Stage”, which is characterized by N10 kmdeep reverse fault seismicity at themargin, recent topographic uplift on the continental side, thick continental crust at themargin, and bulging on the oceanic side due to loading by the overthrusting continent. The combined results of the numerical simulations and passivemargin analysis indicate that the SE Brazilian margin is a prototype candidate for subduction initiation.

La erupción submarina de 2011-12 en El Hierro (Islas Canarias): cronología de eventos y gestión de la crisis

VIII Congreso geológico de España, Oviedo, 17-19 julio 2012

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

Astronomical site testing in the era of the extremely large telescopes

The quality of astronomical sites is the first step to be considered to have the best performances from the telescopes. In particular, the efficiency of large telescopes in UV, IR, radio etc. is critically dependent on atmospheric transparency. It is well known that the random optical effects induced on the light propagation by turbulent atmosphere also limit telescope’s performances. Nowadays, clear appears the importance to correlate the main atmospheric physical parameters with the optical quality reachable by large aperture telescopes. The sky quality evaluation improved with the introduction of new techniques, new instrumentations and with the understanding of the link between the meteorological (or synoptical parameters and the observational conditions thanks to the application of the theories of electromagnetic waves propagation in turbulent medias: what we actually call astroclimatology. At the present the site campaigns are evolved and are performed using the classical scheme of optical seeing properties, meteorological parameters, sky transparency, sky darkness and cloudiness. New concept are added and are related to the geophysical properties such as seismicity, microseismicity, local variability of the climate, atmospheric conditions related to the ground optical turbulence and ground wind regimes, aerosol presence, use of satellite data. The purpose of this project is to provide reliable methods to analyze the atmospheric properties that affect ground-based optical astronomical observations and to correlate them with the main atmospheric parameters generating turbulence and affecting the photometric accuracy. The first part of the research concerns the analysis and interpretation of longand short-time scale meteorological data at two of the most important astronomical sites located in very different environments: the Paranal Observatory in the Atacama Desert (Chile), and the Observatorio del Roque de Los Muchachos(ORM) located in La Palma (Canary Islands, Spain). The optical properties of airborne dust at ORM have been investigated collecting outdoor data using a ground-based dust monitor. Because of its dryness, Paranal is a suitable observatory for near-IR observations, thus the extinction properties in the spectral range 1.00-2.30 um have been investigated using an empirical method. Furthermore, this PhD research has been developed using several turbulence profilers in the selection of the site for the European Extremely Large Telescope(E-ELT). During the campaigns the properties of the turbulence at different heights at Paranal and in the sites located in northern Chile and Argentina have been studied. This given the possibility to characterize the surface layer turbulence at Paranal and its connection with local meteorological conditions.

Pushover analysis of an existing reinforced concrete bridge:Jamboree Road Overcrossing in Irvine, California

The work for the present thesis started in California, during my semester as an exchange student overseas. California is known worldwide for its seismicity and its effort in the earthquake engineering research field. For this reason, I immediately found interesting the Structural Dynamics Professor, Maria Q. Feng's proposal, to work on a pushover analysis of the existing Jamboree Road Overcrossing bridge. Concrete is a popular building material in California, and for the most part, it serves its functions well. However, concrete is inherently brittle and performs poorly during earthquakes if not reinforced properly. The San Fernando Earthquake of 1971 dramatically demonstrated this characteristic. Shortly thereafter, code writers revised the design provisions for new concrete buildings so to provide adequate ductility to resist strong ground shaking. There remain, nonetheless, millions of square feet of non-ductile concrete buildings in California. The purpose of this work is to perform a Pushover Analysis and compare the results with those of a Nonlinear Time-History Analysis of an existing bridge, located in Southern California. The analyses have been executed through the software OpenSees, the Open System for Earthquake Engineering Simulation. The bridge Jamboree Road Overcrossing is classified as a Standard Ordinary Bridge. In fact, the JRO is a typical three-span continuous cast-in-place prestressed post-tension box-girder. The total length of the bridge is 366 ft., and the height of the two bents are respectively 26,41 ft. and 28,41 ft.. Both the Pushover Analysis and the Nonlinear Time-History Analysis require the use of a model that takes into account for the nonlinearities of the system. In fact, in order to execute nonlinear analyses of highway bridges it is essential to incorporate an accurate model of the material behavior. It has been observed that, after the occurrence of destructive earthquakes, one of the most damaged elements on highway bridges is a column. To evaluate the performance of bridge columns during seismic events an adequate model of the column must be incorporated. Part of the work of the present thesis is, in fact, dedicated to the modeling of bents. Different types of nonlinear element have been studied and modeled, with emphasis on the plasticity zone length determination and location. Furthermore, different models for concrete and steel materials have been considered, and the selection of the parameters that define the constitutive laws of the different materials have been accurate. The work is structured into four chapters, to follow a brief overview of the content. The first chapter introduces the concepts related to capacity design, as the actual philosophy of seismic design. Furthermore, nonlinear analyses both static, pushover, and dynamic, time-history, are presented. The final paragraph concludes with a short description on how to determine the seismic demand at a specific site, according to the latest design criteria in California. The second chapter deals with the formulation of force-based finite elements and the issues regarding the objectivity of the response in nonlinear field. Both concentrated and distributed plasticity elements are discussed into detail. The third chapter presents the existing structure, the software used OpenSees, and the modeling assumptions and issues. The creation of the nonlinear model represents a central part in this work. Nonlinear material constitutive laws, for concrete and reinforcing steel, are discussed into detail; as well as the different scenarios employed in the columns modeling. Finally, the results of the pushover analysis are presented in chapter four. Capacity curves are examined for the different model scenarios used, and failure modes of concrete and steel are discussed. Capacity curve is converted into capacity spectrum and intersected with the design spectrum. In the last paragraph, the results of nonlinear time-history analyses are compared to those of pushover analysis.

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.