Medium change monitoring using ambient seismic noise and coda wave interferometry: examples from intraplate NE Brazil and the Mid-Atlantic Ridge

This thesis presents and discusses the results of ambient seismic noise correlation for two different environments: intraplate and Mid-Atlantic Ridge. The coda wave interferometry method has also been tested for the intraplate data. Ambient noise correlation is a method that allows to retrieve the structural response between two receivers from ambient noise records, as if one of the station was a virtual source. It has been largely used in seismology to image the subsurface and to monitor structural changes associated mostly with volcanic eruptions and large earthquakes. In the intraplate study, we were able to detect localized structural changes related to a small earthquake swarm, which main event is mR 3.7, North-East of Brazil. We also showed that the 1-bit normalization and spectral whitening result on the loss of waveform details and that the phase auto-correlation, which is amplitude unbiased, seems to be more sensitive and robust for our analysis of a small earthquake swarm. The analysis of 6 months of data using cross-correlations detect clear medium changes soon after the main event while the auto-correlations detect changes essentially after 1 month. It could be explained by fluid pressure redistribution which can be initiated by hydromechanical changes and opened path ways to shallower depth levels due to later occurring earthquakes. In the Mid-Atlantic Ridge study, we investigate structural changes associated with a mb 4.9 earthquake in the region of the Saint Paul transform fault. The data have been recorded by a single broadband seismic station located at less than 200 km from the Mid-Atlantic ridge. The results of the phase auto-correlation for a 5-month period, show a strong co-seismic medium change followed by a relatively fast post-seismic recovery. This medium change is likely related to the damages caused by the earthquake’s ground shaking. The healing process (filling of the new cracks) that lasted 60 days can be decomposed in two phases, a fast recovery (70% in ~30 days) in the early post-seismic stage and a relatively slow recovery later (30% in ~30 days). In the coda wave interferometry study, we monitor temporal changes of the subsurface caused by the small intraplate earthquake swarm mentioned previously. The method was first validated with synthetics data. We were able to detect a change of 2.5% in the source position and a 15% decrease of the scatterers’ amount. Then, from the real data, we observed a rapid decorrelation of the seismic coda after the mR 3.7 seismic event. This indicates a rapid change of the subsurface in the fault’s region induced by the earthquake.

Estudo da sismicidade na barragem do Açu no período de 2012 a 2013

The seismic activity in the Northeastern of Brazil has been a constant target of study, since it is the most active region of the country. However, some areas have their earthquakes related to human action, what means they are induced. The Açu dam is a classic example of reservoir-induced seismicity and it has been the subject of several studies. Recently, after a considerable period of inactivity, the LabSis / UFRN recorded events related to the dam, which led to the installation of a network around the reservoir. From the data provided by this network, it was observed that the seismic activity is related to a new epicental area inside the lake. The epicentral parameters and the focal mechanism were determined. It was found that the events were related to the reactivation of a basement structure in a new seismogenic subvertical fault with NE-SW-striking subparallel to the São Rafael Fault. These results were used in the preparation of a scientific paper, which discussed the relationship between this seismicity with the geology of the region and with the reservoir water level. The paper showed that the diffusion of pore pressure was the main mechanism that controlled the triggering of the induced seismicity at the reservoir.

Estudo da sismicidade na barragem do Açu no período de 2012 a 2013

The seismic activity in the Northeastern of Brazil has been a constant target of study, since it is the most active region of the country. However, some areas have their earthquakes related to human action, what means they are induced. The Açu dam is a classic example of reservoir-induced seismicity and it has been the subject of several studies. Recently, after a considerable period of inactivity, the LabSis / UFRN recorded events related to the dam, which led to the installation of a network around the reservoir. From the data provided by this network, it was observed that the seismic activity is related to a new epicental area inside the lake. The epicentral parameters and the focal mechanism were determined. It was found that the events were related to the reactivation of a basement structure in a new seismogenic subvertical fault with NE-SW-striking subparallel to the São Rafael Fault. These results were used in the preparation of a scientific paper, which discussed the relationship between this seismicity with the geology of the region and with the reservoir water level. The paper showed that the diffusion of pore pressure was the main mechanism that controlled the triggering of the induced seismicity at the reservoir.

Semi-engineered earthquake-resistant structures: one-storey buildings built up with gabion-box walls

This thesis studies the static and seismic behavior of simple structures made with gabion box walls. The analysis was performed considering a one-story building with standard dimensions in plan (6m x 5m) and a lightweight timber roof. The main focus of the present investigation is to find the principals aspects of the seismic behavior of a one story building made with gabion box walls, in order to prevent a failure due to seismic actions and in this way help to reduce the seismic risk of developing countries where this natural disaster have a significant intensity. Regarding the gabion box wall, it has been performed some calculations and analysis in order to understand the static and dynamic behavior. From the static point of view, it has been performed a verification of the normal stress computing the normal stress that arrives at the base of the gabion wall and the corresponding capacity of the ground. Moreover, regarding the seismic analysis, it has been studied the in-plane and out-of-plane behavior. The most critical aspect was discovered to be the out-of-plane behavior, for which have been developed models considering the “rigid- no tension model” for masonry, finding a kinematically admissible multiplier that will create a collapse mechanism for the structure. Furthermore, it has been performed a FEM and DEM models to find the maximum displacement at the center of the wall, maximum tension stresses needed for calculating the steel connectors for joining consecutive gabions and the dimensions (length of the wall and distance between orthogonal walls or buttresses) of a geometrical configuration for the standard modulus of the structure, in order to ensure an adequate safety margin for earthquakes with a PGA around 0.4-0.5g. Using the results obtained before, it has been created some rules of thumb, that have to be satisfy in order to ensure a good behavior of these structure.

Seismic Retrofitting of a Historic Structure with Limited Construction Documentation

Many buildings constructed during the middle of the 20th century were constructed with criteria that fall short of current requirements. Although shortcomings are possible in all aspects of the design, the inadequacies in terms of seismic design present a more pressing issue to human life. This risk has been seen in various earthquakes that have struck Italy recently, and subsequently, the codes have been altered to account for this underestimated danger. Structures built after these changes remain at risk and must be retrofitted depending on their use. This report centers around the Giovanni Michelucci Institute of Mathematics at the University of Bologna and the work required to modify the building so that it can withstand 60% of the current design requirements. The goal of this particular report is to verify the previous reports written in Italian and present an accurate analysis along with intervention suggestions for this particular building. The work began with an investigation into the previous sources and work to find out how the structure had been interpreted. After understanding the building, corrections were made where required, and the failing elements were organized graphically to more easily show where the building needed the most work. Once the critical zones were mapped, remediation techniques were tested on the top floor, and the modeling techniques and effects of the interventions were presented to assist in further work on the structure.

Event layer chronology from Lake Suigetsu (Japan, SG06 Project) for the last 40 kyr

Event layers in lake sediments are indicators of past extreme events, mostly the results of floods or earthquakes. Detailed characterisation of the layers allows the discrimination of the sedimentation processes involved, such as surface runoff, landslides or subaqueous slope failures. These processes can then be interpreted in terms of their triggering mechanisms. Here we present a 40 kyr event layer chronology from Lake Suigetsu, Japan. The event layers were characterised using a multi-proxy approach, employing light microscopy and µXRF for microfacies analysis. The vast majority of event layers in Lake Suigetsu was produced by flood events (362 out of 369), allowing the construction of the first long-term, quantitative (with respect to recurrence) and well dated flood chronology from the region. The flood layer frequency shows a high variability over the last 40 kyr, and it appears that extreme precipitation events were decoupled from the average long-term precipitation. For instance, the flood layer frequency is highest in the Glacial at around 25 kyr BP, at which time Japan was experiencing a generally cold and dry climate. Other cold episodes, such as Heinrich Event 1 or the Late Glacial stadial, show a low flood layer frequency. Both observations together exclude a simple, straightforward relationship with average precipitation and temperature. We argue that, especially during Glacial times, changes in typhoon genesis/typhoon tracks are the most likely control on the flood layer frequency, rather than changes in the monsoon front or snow melts. Spectral analysis of the flood chronology revealed periodic variations on centennial and millennial time scales, with 220 yr, 450 yr and a 2000 yr cyclicity most pronounced. However, the flood layer frequency appears to have not only been influenced by climate changes, but also by changes in erosion rates due to, for instance, earthquakes.

Activity concentrations of radionuclides in 5 short cores taken in 2007-2008 measured by ICP-MS and HPGe gamma spectrometry

Natural radionuclides and man-made 137Cs were analyzed in five short sediment cores taken in northern part of the Gulf of Eilat (Gulf of Aqaba) in order to provide information on sedimentation and mixing rates and sediment sources. The maximum estimates of sedimentation rates based on excess 210Pb were found to vary between 0.105 ± 0.020 and 0.35 ± 0.23 cm · year**-1. Even the lowest estimates are significantly higher than those expected from dust deposition, suggesting other sources and processes being responsible for most of the allochthonous material accumulation, including periodical floods following heavy rain events, internal erosion or triggers, like earthquakes. In 137Cs depth profiles no 1963 related nuclear weapon test maxima were found; instead, the activities decrease monotonically, suggesting that a major process leading to radionuclides' depth distribution might be mixing. The mixing rates calculated from 137Cs, excess 210Pb and excess 228Th reach values up to 2.18 ± 0.69 cm**2 · year**-1.

(Table 1) Experimental details and results

Factor-of-safety analyses of submarine slope failure depend critically on the shear strength of the slope material, which is often evaluated with residual strength values and for normally consolidated sediments. Here, we report on direct measurements of both shear strength and cohesion for a quartz-clay mixture over a wide range of overconsolidation ratios (OCRs). For normally consolidated sediment at low stresses, cohesion is the dominant source of shear strength compared to friction. Significant increases in peak shear strength occur for OCR > 4, and the primary source of this strength increase is due to increased cohesion, rather than friction. The proportion of added shear strength due to cohesion depends log-linearly on the OCR. We show that at shallow depths where OCR values can be high, overconsolidated clays can be stronger than pure or nearly pure quartz sediments, which are cohesionless under near-surface conditions. Our data also suggest that areas which have experienced significant unroofing due to previous mass movements are less likely to experience subsequent failure at shallow depths due to increased peak strength, and if failure occurs it is expected to be deeper where the OCR is lower. In seismically active areas, this is one potential explanation for the general observation of lower slope failure recurrence compared to rates expected from triggering due to local earthquakes.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

GPS inferred velocity and strain rate fields in eastern Canada

La vallée du fleuve Saint-Laurent, dans l’est du Canada, est l’une des régions sismiques les plus actives dans l’est de l’Amérique du Nord et est caractérisée par de nombreux tremblements de terre intraplaques. Après la rotation rigide de la plaque tectonique, l’ajustement isostatique glaciaire est de loin la plus grande source de signal géophysique dans l’est du Canada. Les déformations et les vitesses de déformation de la croûte terrestre de cette région ont été étudiées en utilisant plus de 14 ans d’observations (9 ans en moyenne) de 112 stations GPS fonctionnant en continu. Le champ de vitesse a été obtenu à partir de séries temporelles de coordonnées GPS quotidiennes nettoyées en appliquant un modèle combiné utilisant une pondération par moindres carrés. Les vitesses ont été estimées avec des modèles de bruit qui incluent les corrélations temporelles des séries temporelles des coordonnées tridimensionnelles. Le champ de vitesse horizontale montre la rotation antihoraire de la plaque nord-américaine avec une vitesse moyenne de 16,8±0,7 mm/an dans un modèle sans rotation nette (no-net-rotation) par rapport à l’ITRF2008. Le champ de vitesse verticale confirme un soulèvement dû à l’ajustement isostatique glaciaire partout dans l’est du Canada avec un taux maximal de 13,7±1,2 mm/an et un affaissement vers le sud, principalement au nord des États-Unis, avec un taux typique de −1 à −2 mm/an et un taux minimum de −2,7±1,4 mm/an. Le comportement du bruit des séries temporelles des coordonnées GPS tridimensionnelles a été analysé en utilisant une analyse spectrale et la méthode du maximum de vraisemblance pour tester cinq modèles de bruit: loi de puissance; bruit blanc; bruit blanc et bruit de scintillation; bruit blanc et marche aléatoire; bruit blanc, bruit de scintillation et marche aléatoire. Les résultats montrent que la combinaison bruit blanc et bruit de scintillation est le meilleur modèle pour décrire la partie stochastique des séries temporelles. Les amplitudes de tous les modèles de bruit sont plus faibles dans la direction nord et plus grandes dans la direction verticale. Les amplitudes du bruit blanc sont à peu près égales à travers la zone d’étude et sont donc surpassées, dans toutes les directions, par le bruit de scintillation et de marche aléatoire. Le modèle de bruit de scintillation augmente l’incertitude des vitesses estimées par un facteur de 5 à 38 par rapport au modèle de bruit blanc. Les vitesses estimées de tous les modèles de bruit sont statistiquement cohérentes. Les paramètres estimés du pôle eulérien de rotation pour cette région sont légèrement, mais significativement, différents de la rotation globale de la plaque nord-américaine. Cette différence reflète potentiellement les contraintes locales dans cette région sismique et les contraintes causées par la différence des vitesses intraplaques entre les deux rives du fleuve Saint-Laurent. La déformation de la croûte terrestre de la région a été étudiée en utilisant la méthode de collocation par moindres carrés. Les vitesses horizontales interpolées montrent un mouvement cohérent spatialement: soit un mouvement radial vers l’extérieur pour les centres de soulèvement maximal au nord et un mouvement radial vers l’intérieur pour les centres d’affaissement maximal au sud, avec une vitesse typique de 1 à 1,6±0,4 mm/an. Cependant, ce modèle devient plus complexe près des marges des anciennes zones glaciaires. Basées selon leurs directions, les vitesses horizontales intraplaques peuvent être divisées en trois zones distinctes. Cela confirme les conclusions d’autres chercheurs sur l’existence de trois dômes de glace dans la région d’étude avant le dernier maximum glaciaire. Une corrélation spatiale est observée entre les zones de vitesses horizontales intraplaques de magnitude plus élevée et les zones sismiques le long du fleuve Saint-Laurent. Les vitesses verticales ont ensuite été interpolées pour modéliser la déformation verticale. Le modèle montre un taux de soulèvement maximal de 15,6 mm/an au sud-est de la baie d’Hudson et un taux d’affaissement typique de 1 à 2 mm/an au sud, principalement dans le nord des États-Unis. Le long du fleuve Saint-Laurent, les mouvements horizontaux et verticaux sont cohérents spatialement. Il y a un déplacement vers le sud-est d’une magnitude d’environ 1,3 mm/an et un soulèvement moyen de 3,1 mm/an par rapport à la plaque l’Amérique du Nord. Le taux de déformation verticale est d’environ 2,4 fois plus grand que le taux de déformation horizontale intraplaque. Les résultats de l’analyse de déformation montrent l’état actuel de déformation dans l’est du Canada sous la forme d’une expansion dans la partie nord (la zone se soulève) et d’une compression dans la partie sud (la zone s’affaisse). Les taux de rotation sont en moyenne de 0,011°/Ma. Nous avons observé une compression NNO-SSE avec un taux de 3.6 à 8.1 nstrain/an dans la zone sismique du Bas-Saint-Laurent. Dans la zone sismique de Charlevoix, une expansion avec un taux de 3,0 à 7,1 nstrain/an est orientée ENE-OSO. Dans la zone sismique de l’Ouest du Québec, la déformation a un mécanisme de cisaillement avec un taux de compression de 1,0 à 5,1 nstrain/an et un taux d’expansion de 1.6 à 4.1 nstrain/an. Ces mesures sont conformes, au premier ordre, avec les modèles d’ajustement isostatique glaciaire et avec la contrainte de compression horizontale maximale du projet World Stress Map, obtenue à partir de la théorie des mécanismes focaux (focal mechanism method).

Simplified design method for energy dissipating devices in retrofitting of seismically isolated bridges

Abstract: Highway bridges have great values in a country because in case of any natural disaster they may serve as lines to save people’s lives. Being vulnerable under significant seismic loads, different methods can be considered to design resistant highway bridges and rehabilitate the existing ones. In this study, base isolation has been considered as one efficient method in this regards which in some cases reduces significantly the seismic load effects on the structure. By reducing the ductility demand on the structure without a notable increase of strength, the structure is designed to remain elastic under seismic loads. The problem associated with the isolated bridges, especially with elastomeric bearings, can be their excessive displacements under service and seismic loads. This can defy the purpose of using elastomeric bearings for small to medium span typical bridges where expansion joints and clearances may result in significant increase of initial and maintenance cost. Thus, supplementing the structure with dampers with some stiffness can serve as a solution which in turn, however, may increase the structure base shear. The main objective of this thesis is to provide a simplified method for the evaluation of optimal parameters for dampers in isolated bridges. Firstly, performing a parametric study, some directions are given for the use of simple isolation devices such as elastomeric bearings to rehabilitate existing bridges with high importance. Parameters like geometry of the bridge, code provisions and the type of soil on which the structure is constructed have been introduced to a typical two span bridge. It is concluded that the stiffness of the substructure, soil type and special provisions in the code can determine the employment of base isolation for retrofitting of bridges. Secondly, based on the elastic response coefficient of isolated bridges, a simplified design method of dampers for seismically isolated regular highway bridges has been presented in this study. By setting objectives for reduction of displacement and base shear variation, the required stiffness and damping of a hysteretic damper can be determined. By modelling a typical two span bridge, numerical analyses have followed to verify the effectiveness of the method. The method has been used to identify equivalent linear parameters and subsequently, nonlinear parameters of hysteretic damper for various designated scenarios of displacement and base shear requirements. Comparison of the results of the nonlinear numerical model without damper and with damper has shown that the method is sufficiently accurate. Finally, an innovative and simple hysteretic steel damper was designed. Five specimens were fabricated from two steel grades and were tested accompanying a real scale elastomeric isolator in the structural laboratory of the Université de Sherbrooke. The test procedure was to characterize the specimens by cyclic displacement controlled tests and subsequently to test them by real-time dynamic substructuring (RTDS) method. The test results were then used to establish a numerical model of the system which went through nonlinear time history analyses under several earthquakes. The outcome of the experimental and numerical showed an acceptable conformity with the simplified method.

Tsunami characteristics along the Peru-Chile trench: analysis of the 2015 Mw8.3 Illapel, the 2014 Mw8.2 iquique and the 2010 Mw8.8 maule tsunamis in the near-field

Tsunamis occur quite frequently following large magnitude earthquakes along the Chilean coast. Most of these earthquakes occur along the Peru-Chile Trench, one of the most seismically active subduction zones of the world. This study aims to understand better the characteristics of the tsunamis triggered along the Peru-Chile Trench. We investigate the tsunamis induced by the Mw8.3 Illapel, the Mw8.2 Iquique and the Mw8.8 Maule Chilean earthquakes that happened on September 16th, 2015, April 1st, 2014 and February 27th, 2010, respectively. The study involves the relation between the co-seismic deformation and the tsunami generation, the near-field tsunami propagation, and the spectral analysis of the recorded tsunami signals in the near-field. We compare the tsunami characteristics to highlight the possible similarities between the three events and, therefore, attempt to distinguish the specific characteristics of the tsunamis occurring along the Peru-Chile Trench. We find that these three earthquakes present faults with important extensions beneath the continent which result in the generation of tsunamis with short wavelengths, relative to the fault widths involved, and with reduced initial potential energy. In addition, the presence of the Chilean continental margin, that includes the shelf of shallow bathymetry and the continental slope, constrains the tsunami propagation and the coastal impact. All these factors contribute to a concentrated local impact but can, on the other hand, reduce the far-field tsunami effects from earthquakes along Peru-Chile Trench.

Comportamento e reforço de nós viga-pilar

A presente dissertação tem como principal objetivo a análise do comportamento do nó viga-pilar, com a finalidade de dimensionar uma solução de reforço sísmico a aplicar em nós de estruturas antigas de betão armado. Uma parte considerável das construções em betão armado construídas em Portugal até 1970 são muito frágeis relativamente às ações sísmicas. Devido a este fato, os edifícios estão associados a uma conceção desajustada face à regulamentação atual pois à época não estava prevista a consideração de ações sísmicas na execução dos projetos de estruturas. Nesse período, era prática comum na construção o uso de armadura lisa conjugada com uma má pormenorização das armaduras e baixa resistência mecânica do betão. Tal levava a que as estruturas de betão armado construídas em Portugal até à década de setenta, possuírem uma baixa capacidade para resistir a solicitações cíclicas, como as induzidas pelos sismos. Os nós viga-pilar das estruturaras de betão armado são considerados as zonas mais crítica e vulnerável a ações sísmicas. Por isso é frequente em estruturas antigas, a necessidade de as reforçar, para melhorar a performance do edifício durante os sismos. Têm sido estudadas por investigadores várias técnicas e soluções para reforçar os nós viga-pilar de estruturas antigas de betão armado. Serão analisados nesta dissertação os danos observados em dois nós viga-pilar ensaiados experimentalmente, onde são simuladas as ações cíclicas horizontais dos sismos. Com base nestes dados, é proposta e detalhada uma solução de reforço para melhorar o comportamento histerético dessa ligações, de forma a conferir uma maior capacidade resistente e/ou ductilidade.

Building monitoring during the partial implosion stage

The partial collapse of a building in Colombia caused severe damage to its structural components -- An implosion was realized to induce the collapse of 50% of the deteriorated building -- To evaluate the influence of the implosion on the remaining structure, a monitoring survey was realized using triaxial accelerometers -- Time signals associated with ambient, seismic and forced vibration were obtained -- A study of the records in the time and the frequency domain was made -- The analysis of the information allowed determining some structural properties that were useful to calibrate the analytical model of the structure

Projeto de reforço de um edifício de habitação

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Estruturas