The Moho reflectivity from ambient seismic noise autocorrelations beneath the Western Alps

The lower crustal structure beneath the Western Alps -- including the Moho -- bears the signature of past and present geodynamic processes. It has been the subject of many studies until now. However, its current knowledge still leaves significant open questions. In order to derive new information, independent from previous determinations, here I wish to address this topic using a different method --- ambient seismic noise autocorrelation --- that is for the first time applied to reveal Moho depth in the Western Alps. Moho reflections are identified by picking reflectivity changes in ambient seismic noise autocorrelations. The seismic data is retrieved from more than 200 broadband seismic stations, from the China--Italy--France Alps (CIFALPS) linear seismic network, and from a subset of the AlpArray Seismic Network (AASN). The automatically-picked reflectivity changes along the CIFALPS transect in the southwestern Alps show the best results in the 0.5--1 Hz frequency band. The autocorrelation reflectivity profile of the CIFALPS transect shows a steeper subduction profile,~55 to ~70 km, of the European Plate underneath the Adriatic Plate. The dense spacing of the CIFALPS network facilitates the detection of lateral continuity of crustal structure, and of the Ivrea mantle wedge reaching shallow crustal depths in the southwestern Alps. The data of the AASN stations are filtered in the 0.4--1 and 0.5--1 Hz frequency bands. Although the majority of the stations give the same Moho depth for the different frequency bands, the few stations with different Moho depths shows the care that has to be taken when choosing the frequency band for filtering the autocorrelation stacks. The new Moho depth maps by using the AASN stations are a compilation of the first and second picked reflectivity changes. The results show the complex crust-mantle structure with clear differences between the northwestern and southwestern Alps.

Seismic vulnerability of buildings and non-structural elements: insights from the 2011 Lorca earthquake

This PhD dissertation presents a profound study of the vulnerability of buildings and non-structural elements stemming from the investigation of the Mw 5.2 Lorca 2011 earthquake; which constitutes one of the most significant earthquakes in Spain. It left nine fatalities due to falling debris from reinforced concrete buildings, 394 injured and material damage valued at 800 million euros. Within this framework, the most relevant initiatives concerning the vulnerability of buildings and the exposure of Lorca are studied. This work revealed two lines of research: the elaboration of a rational method to determine the adequacy of a specific fragility curve for the particular seismic risk study of a region; and the relevance of researching the seismic performance of non-structural elements. As a consequence, firstly, a method to assess and select fragility curves for seismic risk studies from the catalogue of those available in the literature is elaborated and calibrated by means of a case study. The said methodology is based on a multidimensional index and provides a ranking that classifies the curves in terms of adequacy. Its results for the case of Lorca led to the elaboration of new fragility curves for unreinforced masonry buildings. Moreover, a simplified method to account for the unpredictable directionality of the seism in the creation of fragility curves is contributed. Secondly, the characterisation of the seismic capacity and demand of the non-structural elements that caused most of the human losses is studied. Concerning the capacity, an analytical approach derived from theoretical considerations to characterise the complete out-of-plane seismic response curve of unreinforced masonry cantilever walls is provided; as well as a simplified and more practical trilinear version of it. Concerning the demand, several methods for characterising the Floor Response Spectra of reinforced concrete buildings are tested through case studies.

High-resolution stratigraphy and seismic site characterization of late Quaternary paleovalley systems (Adriatic coastal plain, Italy)

In highly urbanized coastal lowlands, effective site characterization is crucial for assessing seismic risk. It requires a comprehensive stratigraphic analysis of the shallow subsurface, coupled with the precise assessment of the geophysical properties of buried deposits. In this context, late Quaternary paleovalley systems, shallowly buried fluvial incisions formed during the Late Pleistocene sea-level fall and filled during the Holocene sea-level rise, are crucial for understanding seismic amplification due to their soft sediment infill and sharp lithologic contrasts. In this research, we conducted high-resolution stratigraphic analyses of two regions, the Pescara and Manfredonia areas along the Adriatic coastline of Italy, to delineate the geometries and facies architecture of two paleovalley systems. Furthermore, we carried out geophysical investigations to characterize the study areas and perform seismic response analyses. We tested the microtremor-based horizontal-to-vertical spectral ratio as a mapping tool to reconstruct the buried paleovalley geometries. We evaluated the relationship between geological and geophysical data and identified the stratigraphic surfaces responsible for the observed resonances. To perform seismic response analysis of the Pescara paleovalley system, we integrated the stratigraphic framework with microtremor and shear wave velocity measurements. The seismic response analysis highlights strong seismic amplifications in frequency ranges that can interact with a wide variety of building types. Additionally, we explored the applicability of artificial intelligence in performing facies analysis from borehole images. We used a robust dataset of high-resolution digital images from continuous sediment cores of Holocene age to outline a novel, deep-learning-based approach for performing automatic semantic segmentation directly on core images, leveraging the power of convolutional neural networks. We propose an automated model to rapidly characterize sediment cores, reproducing the sedimentologist's interpretation, and providing guidance for stratigraphic correlation and subsurface reconstructions.

Comparison of the seismic performances of a steel frame building equipped with different geometrical configurations of Crescent Shaped Braces

The Crescent Shaped Brace (CSB) is a new simple steel hysteretic device proposed to be used as an enhanced diagonal brace in framed structures. The CSB allows the practical designer to choose the lateral stiffness independently from the yield strength of the device, due to its peculiar ad-hoc shape. In the present thesis, a complete study referring to different CSB configurations has been presented. After the validation of the hysteretic capacities of the Crescent Shaped Braces, the seismic concept of the "enhanced first story isolation" system has been proposed within the PBSD. It relies on the total separation between the Vertical Resisting System (VRS) and the Horizontal Resisting System (HRS) in order to attain a certain objective curve of the structure. An applicative example has been studied following this concept and exploiting the advantages of the CSBs as seismic dissipative devices used for the HRS. Then several geometrical configurations called Single CSB system, Single 2 CSB system, Double CSB system, Coupled CSB system, Coupled with high length CSB system, and the final one was Cross bracing system have been introduced and modelled with SAP2000 and the results have been compared.

Anisotropic seismic tomography with the reversible jump Markov chain Monte Carlo

The established isotropic tomographic models show the features of subduction zones in terms of seismic velocity anomalies, but they are generally subjected to the generation of artifacts due to the lack of anisotropy in forward modelling. There is evidence for the significant influence of seismic anisotropy in the mid-upper mantle, especially for boundary layers like subducting slabs. As consequence, in isotropic models artifacts may be misinterpreted as compositional or thermal heterogeneities. In this thesis project the application of a trans-dimensional Metropolis-Hastings method is investigated in the context of anisotropic seismic tomography. This choice arises as a response to the important limitations introduced by traditional inversion methods which use iterative procedures of optimization of a function object of the inversion. On the basis of a first implementation of the Bayesian sampling algorithm, the code is tested with some cartesian two-dimensional models, and then extended to polar coordinates and dimensions typical of subduction zones, the main focus proposed for this method. Synthetic experiments with increasing complexity are realized to test the performance of the method and the precautions for multiple contexts, taking into account also the possibility to apply seismic ray-tracing iteratively. The code developed is tested mainly for 2D inversions, future extensions will allow the anisotropic inversion of seismological data to provide more realistic imaging of real subduction zones, less subjected to generation of artifacts.

Utilização de métodos acústicos em estudos de dinâmica costeira: exemplo na desembocadura lagunar de cananéia

Em função de suas condições de interface entre águas doces e salinas, desembocaduras estuarinas e lagunares constituem sistemas geomorfológicos altamente complexos e dinâmicos. Como conseqüência da variabilidade espacial e temporal dos fluxos de maré, o leito responde com uma grande variabilidade nas características morfológicas e sedimentares. Neste sentido, é possível relacionar diretamente a circulação de fundo e o transporte sedimentar com as feições submersas geradas. Perfis de ecossondagem, sonar de varredura lateral e sísmica de alta resolução, executados na desembocadura lagunar de Cananéia, revelaram a existência de uma dinâmica de fundo extremamente complexa, caracterizada por marcas onduladas e ondas de areia de alturas métricas. As maiores ondas de areia, localizadas em uma depressão na desembocadura lagunar, apresentam inversão de polaridade em sua assimetria, com a presença de ondas simétricas de grande tamanho no ponto de inversão. Este padrão morfológico não apresenta variação temporal em escala anual, sugerindo a persistência de um padrão de fluxos sobre o leito. Esta dinâmica revela, também, a constância de fluxos convergentes que aparentemente independem das condições de maré enchente ou vazante. Os resultados permitiram o estabelecimento de um primeiro modelo qualitativo de circulação de fundo na área, com aplicações potenciais na navegação e estudos de proteção da costa.

The pulsations of the B5IVe star HD 181231 observed with CoRoT and ground-based spectroscopy

Context. HD 181231 is a B5IVe star, which has been observed with the CoRoT satellite during similar to 5 consecutive months and simultaneously from the ground in spectroscopy and spectropolarimetry. Aims. By analysing these data, we aim to detect and characterize as many pulsation frequencies as possible, to search for the presence of beating effects possibly at the origin of the Be phenomenon. Our results will also provide a basis for seismic modelling. Methods. The fundamental parameters of the star are determined from spectral fitting and from the study of the circumstellar emission. The CoRoT photometric data and ground-based spectroscopy are analysed using several Fourier techniques: CLEAN-NG, PASPER, and TISAFT, as well as a time-frequency technique. A search for a magnetic field is performed by applying the LSD technique to the spectropolarimetric data. Results. We find that HD 181231 is a B5IVe star seen with an inclination of similar to 45 degrees. No magnetic field is detected in its photosphere. We detect at least 10 independent significant frequencies of variations among the 54 detected frequencies, interpreted in terms of non-radial pulsation modes and rotation. Two longer-term variations are also detected: one at similar to 14 days resulting from a beating effect between the two main frequencies of short-term variations, the other at similar to 116 days due either to a beating of frequencies or to a zonal pulsation mode. Conclusions. Our analysis of the CoRoT light curve and ground-based spectroscopic data of HD 181231 has led to the determination of the fundamental and pulsational parameters of the star, including beating effects. This will allow a precise seismic modelling of this star.

The B0.5IVe CoRoT target HD 49330 I. Photometric analysis from CoRoT data

Context. Be stars undergo outbursts producing a circumstellar disk from the ejected material. The beating of non-radial pulsations has been put forward as a possible mechanism of ejection. Aims. We analyze the pulsational behavior of the early B0.5IVe star HD 49330 observed during the first CoRoT long run towards the Galactical anticenter (LRA1). This Be star is located close to the lower edge of the beta Cephei instability strip in the HR diagram and showed a 0.03 mag outburst during the CoRoT observations. It is thus an ideal case for testing the aforementioned hypothesis. Methods. We analyze the CoRoT light curve of HD 49330 using Fourier methods and non-linear least square fitting. Results. In this star, we find pulsation modes typical of beta Cep stars (p modes) and SPB stars (g modes) with amplitude variations along the run directly correlated with the outburst. These results provide new clues about the origin of the Be phenomenon as well as strong constraints on the seismic modelling of Be stars.

Regional scale analysis of landform configuration with base-level (isobase) maps

Base-level maps (or ""isobase maps"", as originally defined by Filosofov, 1960), express a relationship between valley order and topography. The base-level map can be seen as a ""simplified"" version of the original topographic surface, from which the ""noise"" of the low-order stream erosion was removed. This method is able to identify areas with possible tectonic influence even within lithologically uniform domains. Base-level maps have been recently applied in semi-detail scale (e.g., 1:50 000 or larger) morphotectonic analysis. In this paper, we present an evaluation of the method's applicability in regional-scale analysis (e.g., 1:250 000 or smaller). A test area was selected in northern Brazil, at the lower course of the Araguaia and Tocantins rivers. The drainage network extracted from SRTM30_PLUS DEMs with spatial resolution of approximately 900 m was visually compared with available topographic maps and considered to be compatible with a 1:1,000 000 scale. Regarding the interpretation of regional-scale morphostructures, the map constructed with 2nd and 3rd-order valleys was considered to present the best results. Some of the interpreted base-level anomalies correspond to important shear zones and geological contacts present in the 1:5 000 000 Geological Map of South America. Others have no correspondence with mapped Precambrian structures and are considered to represent younger, probably neotectonic, features. A strong E-W orientation of the base-level lines over the inflexion of the Araguaia and Tocantins rivers, suggest a major drainage capture. A N-S topographic swath profile over the Tocantins and Araguaia rivers reveals a topographic pattern which, allied with seismic data showing a roughly N-S direction of extension in the area, lead us to interpret this lineament as an E-W, southward-dipping normal fault. There is also a good visual correspondence between the base-level lineaments and geophysical anomalies. A NW-SE lineament in the southeast of the study area partially corresponds to the northern border of the Mosquito lava field, of Jurassic age, and a NW-SE lineament traced in the northeastern sector of the study area can be interpreted as the Picos-Santa Ines lineament, identifiable in geophysical maps but with little expression in hypsometric or topographic maps.

TOPOGRAPHIC REGISTERS OF PALEO-VALLEYS ON THE SOUTHEASTERN BRAZILIAN CONTINENTAL SHELF

The evolution of paleo-incised-valleys in the Sao Paulo State region of the southeastern Brazilian continental shelf is presented in this study in relation to the post Last Glacial Maximum (LGM) sea-level rises based on the submarine topography modeled by a detailed Digital Elevation Model and evidences noted in high resolution seismic profiles. The hypothesis that has guided this study is that the set of paleo-valley characteristics (i.e. the fluvial parameters of modern coastal drainage systems, the topographical shape and dimensions of the valleys and of the subsurface channels) may indicate aspects of the relation between the influence of the fluvial and the eustatic variation regime in geomorphological-stratigraphic registers. Models described in the literature sustain the view that faster marine transgressions tend to increase erosion in estuaries, which may explain the lack of registers of paleo-drainage both in topography and the sub-surface in areas with wider shelves. On the other hand, on narrower shelves, with a higher slope angle, the transgression process can preserve, or even enhance, the incised valley registers during shoreface retreat. In the area studied, we observed that the dimensions and form of the continental shelf varies from the northern to the southern part of the area, affecting aspects of the geomorphological registers of the submerged incised valleys.

Parametrical study of masonry walls subjected to in-plane loading through numerical modeling

This paper deals with the numerical assessment of the influence of parameters such as pre-compression level, aspect ratio, vertical and horizontal reinforcement ratios and boundary conditions on the lateral strength of masonry walls under in-plane loading. The numerical study is performed through the software DIANA (R) based on the Finite Element Method. The validation of the numerical model is carried out from a database of available experimental results on masonry walls tested under cyclic lateral loading. Numerical results revealed that boundary conditions play a central role on the lateral behavior of masonry walls under in-plane loading and determine the influence of level of pre-compression as well as the reinforcement ratio on the wall strength. The lateral capacity of walls decreases with the increase of aspect ratio and with the decrease of pre-compression. Vertical steel bars appear to have almost no influence in the shear strength of masonry walls and horizontal reinforcement only increases the lateral strength of masonry walls if the shear response of the walls is determinant for failure, which is directly related to the boundary conditions. (C) 2011 Elsevier Ltd. All rights reserved.

Determination of in situ uniaxial compressive strength of coal seams based on geophysical data

The central issue for pillar design in underground coal mining is the in situ uniaxial compressive strength (sigma (cm)). The paper proposes a new method for estimating in situ uniaxial compressive strength in coal seams based on laboratory strength and P wave propagation velocity. It describes the collection of samples in the Bonito coal seam, Fontanella Mine, southern Brazil, the techniques used for the structural mapping of the coal seam and determination of seismic wave propagation velocity as well as the laboratory procedures used to determine the strength and ultrasonic wave velocity. The results obtained using the new methodology are compared with those from seven other techniques for estimating in situ rock mass uniaxial compressive strength.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Texture alignment in simple shear

We illustrate the flow behaviour of fluids with isotropic and anisotropic microstructure (internal length, layering with bending stiffness) by means of numerical simulations of silo discharge and flow alignment in simple shear. The Cosserat theory is used to provide an internal length in the constitutive model through bending stiffness to describe isotropic microstructure and this theory is coupled to a director theory to add specific orientation of grains to describe anisotropic microstructure. The numerical solution is based on an implicit form of the Material Point Method developed by Moresi et al. [1].