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.

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.

Sediment and pore water data of five cores from the continental slope off Uruguay

Assessing frequency and extent of mass movement at continental margins is crucial to evaluate risks for offshore constructions and coastal areas. A multidisciplinary approach including geophysical, sedimentological, geotechnical, and geochemical methods was applied to investigate multistage mass transport deposits (MTDs) off Uruguay, on top of which no surficial hemipelagic drape was detected based on echosounder data. Nonsteady state pore water conditions are evidenced by a distinct gradient change in the sulfate (SO4**2-) profile at 2.8 m depth. A sharp sedimentological contact at 2.43 m coincides with an abrupt downward increase in shear strength from approx. 10 to >20 kPa. This boundary is interpreted as a paleosurface (and top of an older MTD) that has recently been covered by a sediment package during a younger landslide event. This youngest MTD supposedly originated from an upslope position and carried its initial pore water signature downward. The kink in the SO4**2- profile approx. 35 cm below the sedimentological and geotechnical contact indicates that bioirrigation affected the paleosurface before deposition of the youngest MTD. Based on modeling of the diffusive re-equilibration of SO4**2- the age of the most recent MTD is estimated to be <30 years. The mass movement was possibly related to an earthquake in 1988 (approx. 70 km southwest of the core location). Probabilistic slope stability back analysis of general landslide structures in the study area reveals that slope failure initiation requires additional ground accelerations. Therefore, we consider the earthquake as a reasonable trigger if additional weakening processes (e.g., erosion by previous retrogressive failure events or excess pore pressures) preconditioned the slope for failure. Our study reveals the necessity of multidisciplinary approaches to accurately recognize and date recent slope failures in complex settings such as the investigated area.

Métodos heurísticos para un problema multicriterio de distribución de ayuda humanitaria

Large scale disasters, such as the one caused by the Typhoon Haiyan, which devastated portions of the Philippines in 2013, or the catastrophic 2010 Haiti earthquake, which caused major damage in Port-au-Prince and other settlements in the region, have massive and lasting effects on populations. Nowadays, disasters can be considered as a consequence of inappropriately managed risk. These risks are the product of hazards and vulnerability, which refers to the extent to which a community can be affected by the impact of a hazard. In this way, developing countries, due to their greater vulnerability, suffer the highest costs when a disaster occurs. Disaster relief is a challenge for politics, economies, and societies worldwide. Humanitarian organizations face multiple decision problems when responding to disasters. In particular, once a disaster strikes, the distribution of humanitarian aid to the population affected is one of the most fundamental operations in what is called humanitarian logistics. This term is defined as the process of planning, implementing and controlling the effcient, cost-effective ow and storage of goods and materials as well as related information, from the point of origin to the point of consumption, for the purpose of meeting the end bene- ciaries' requirements and alleviate the suffering of vulnerable people, [the Humanitarian Logistics Conference, 2004 (Fritz Institute)]. During the last decade there has been an increasing interest in the OR/MS community in studying this topic, pointing out the similarities and differences between humanitarian and business logistics, and developing models suited to handle the special characteristics of these problems. Several authors have pointed out that traditional logistic objectives, such as minimizing operation cost, are not the most relevant goals in humanitarian operations. Other factors, such as the time of operation, or the design of safe and equitable distribution plans, come to the front, and new models and algorithms are needed to cope with these special features. Up to six attributes related to the distribution plan are considered in our multi-criteria approach. Even though there are usually simple ways to measure the cost of an operation, the evaluation of some other attributes such as security or equity is not easy. As a result, several attribute measures are proposed and developed, focusing on different aspects of the solutions. Furthermore, when metaheuristic solution methods are used, considering non linear objective functions does not increase the complexity of the algorithms significantly, and thus more accurate measures can be utilized...

Experimental investigation into amplitude-dependent modal properties of an eleven-span motorway bridge

The authors would like to thank their supporters. New Zealand Earthquake Commission (EQC) Research Foundation provided financial support for experimental work (Grant No. UNI/578). New Zealand Transport Agency (NZTA) provided access to the bridge. Piotr Omenzetter’s work within the LRF Centre for Safety and Reliability Engineering at the University of Aberdeen is supported by Lloyd’s Register Foundation. The Foundation helps to protect life and property by supporting engineering-related education, public engagement and the application of research. Ge-Wei Chen’s doctoral study is supported by China Scholarship Council (CSC) (Grant No. 2011637065).

Forced and ambient vibration testing of full scale bridges : Project No. UNI/578

Preprint

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

Seismic Performance of Hybrid Fiber Reinforced Polymer-Concrete Pier Frame Systems

As an alternative to transverse spiral or hoop steel reinforcement, fiber reinforced polymers (FRPs) were introduced to the construction industry in the 1980's. The concept of concrete-filled FRP tube (CFFT) has raised great interest amongst researchers in the last decade. FRP tube can act as a pour form, protective jacket, and shear and flexural reinforcement for concrete. However, seismic performance of CFFT bridge substructure has not yet been fully investigated. Experimental work in this study included four two-column bent tests, several component tests and coupon tests. Four 1/6-scale bridge pier frames, consisting of a control reinforced concrete frame (RCF), glass FRP-concrete frame (GFF), carbon FRP-concrete frame (CFF), and hybrid glass/carbon FRP-concrete frame (HFF) were tested under reverse cyclic lateral loading with constant axial loads. Specimen GFF did not show any sign of cracking at a drift ratio as high as 15% with considerable loading capacity, whereas Specimen CFF showed that lowest ductility with similar load capacity as in Specimen GFF. FRP-concrete columns and pier cap beams were then cut from the pier frame specimens, and were tested again in three point flexure under monotonic loading with no axial load. The tests indicated that bonding between FRP and concrete and yielding of steel both affect the flexural strength and ductility of the components. The coupon tests were carried out to establish the tensile strength and elastic modulus of each FRP tube and the FRP mold for the pier cap beam in the two principle directions of loading. A nonlinear analytical model was developed to predict the load-deflection responses of the pier frames. The model was validated against test results. Subsequently, a parametric study was conducted with variables such as frame height to span ratio, steel reinforcement ratio, FRP tube thickness, axial force, and compressive strength of concrete. A typical bridge was also simulated under three different ground acceleration records and damping ratios. Based on the analytical damage index, the RCF bridge was most severely damaged, whereas the GFF bridge only suffered minor repairable damages. Damping ratio was shown to have a pronounced effect on FRP-concrete bridges, just the same as in conventional bridges. This research was part of a multi-university project, which is founded by the National Science Foundation (NSF) Network for Earthquake Engineering Simulation Research (NEESR) program.

Description of manganese deposits from the Papago mine, Saipan, Mariana Islands

Saipan, situated about 15° N. and 146° E., is one of the larger and more southerly of the Mariana Islands. The 15 small islands of this chain are strung along an eastwardly convex ridge for more than 400 miles north to south, midway between Honshu and New Guinea and about 1,200 miles east of the Philippines. Paralleling this ridge 60 to 100 miles further east is a deep submarine trench, beyond which lies the Pacific Basin proper. To the west is the Philippine Sea, generally deeper than 2,000 fathoms. The trench coincides with a zone of negative gravity anomalies, earthquake foci occur at increasing depths westward from it, and silica- and alumina-rich volcanic rocks characterize the emergent island chain itself. The contrast between these features and those of the Pacific Basin proper to the east is held to favor the conclusion that the Mariana island arc and trench define the structural and petrographic front of Asia.

Position measurements at five permanent GPS stations in the Bavarian Alps as part of the Geodetic Alpine Integrated Network (GAIN) of the ALPS-GPSQUAKENET project

In the frame of the transnational ALPS-GPSQUAKENET project, a component of the Alpine Space Programme of the European Community Initiative Programme (CIP) INTERREG III B, the Deutsches Geodätisches Forschungsinstitut (DGFI) in Munich, Germany, installed in 2005 five continuously operating permanent GPS stations located along the northern Alps boundary in Bavaria. The main objective of the ALPS-GPSQUAKENET project was to build-up a high-performance transnational space geodetic network of Global Positioning System (GPS) receivers in the Alpine region (the so-called Geodetic Alpine Integrated Network, GAIN). Data from this network allows for studying crustal deformations in near real-time to monitor Earthquake hazard and improve natural disaster prevention. The five GPS stations operatied by DGFI are mounted on concrete pillars attached to solid rock. The names of the stations are (from west to east) Hochgrat (HGRA), Breitenberg (BREI), Fahrenberg (FAHR), Hochries (HRIE) and Wartsteinkopf (WART). The provided data series start from October 7, 2005. Data are stored with a temporal spacing of 15 seconds in daily RINEX files.

(Table 1) Radiocarbon results from this study

The Indonesian Arc represents the subduction of the Indian-Australian plate beneath Asia. It has been the scene of catastrophic tectonic activity, including the recent 2004 M=9.1 Aceh earthquake and resulting Indian Ocean tsunami. We have dated planktonic forams associated with historic tephras (Tambora, 1815 and Krakatau, 1883) in marine sediment cores to determine radiocarbon reservoir ages for 2 locations along the arc. Our best estimates for 19th century regional reservoir corrections (DeltaR) are +90 ± 40 yr for surface-dwelling species and +220 ± 40 yr for mixed planktic assemblages containing some upper thermocline species, but scatter in the data suggests that past surface reservoir ages may have varied by about ±100 yr. We used the results of this study to investigate a proposed very large AD 535 eruption at or near Krakatau. We find no evidence for ash from such an eruption, and although this is negative evidence, we consider it sufficiently strong to rule out any possibility that one took place.