A zona de contato entre as formações Motuca e Sambaíba, Permo-Triássico da bacia do Parnaíba, regiões de Filadélfia (TO), Riachão (MA) e Loreto (MA)

O intervalo que compreende o final do Paleozóico e início do Mesozóico foi marcado por mudanças globais paleogeográficas e paleoclimáticas, em parte atribuídas a eventos catastróficos. A intensa continentalização do supercontinente Pangéia, com a implantação de extensos desertos, sucedeu os ambientes costeiros-plataformais do início do Permiano. Os registros desses eventos no norte do Brasil são encontrados nas bacias intracratônicas, particularmente na Bacia do Parnaíba, representados pela zona de contato entre as formações Motuca e Sambaíba. A Formação Motuca é constituída predominantemente por pelitos vermelhos laminados com lentes de gipsita, calcita e marga. Na porção leste da Bacia do Parnaíba, as fácies da Formação Motuca tornam-se mais arenosas com a ocorrência expressiva de arenitos com estratificação cruzada sigmoidal. A Formação Sambaíba consiste em arenitos de coloração creme alaranjada com estratificação plano-paralela e estratificação cruzada de médio a grande porte. Em geral, o contato entre as unidades é brusco, representado pela passagem de arenitos finos com laminação cruzada cavalgante e acamamento flaser/wavy da Formação Motuca para arenitos médios com falhas/microfalhas sinsedimentares e laminações convolutas da Formação Sambaíba. Foram individualizadas 14 fácies sedimentares, agrupadas em quatro associações: AF1 – Lacustre raso / Planície de lama (mudflat), AF2 – “Panela” salina (saline pan), AF3 – Lençol de areia e AF4 – Campo de dunas. A AF1 foi depositada dominantemente por processos de decantação em um extenso ambiente lacustre raso de baixa energia, influenciado por influxos esporádicos de areias oriundos de rios efêmeros. Este sistema lacustre foi, provavelmente, influenciado por períodos de contração e expansão, devido às variações das condições climáticas predominantemente áridas. Os mais expressivos períodos de contração ocorreram na porção oeste da Bacia do Parnaíba, representados pelo desenvolvimento de planícies de lama (mudflats) associadas a lagoas efêmeras saturadas em carbonatos e a “panelas” salinas (saline pans- AF2). Os lençóis de areia (AF3) são planícies arenosas extensas, localmente com área úmidas, intensamente retrabalhadas por processos eólicos. A AF4 é interpretada como parte de um erg composto por dunas/draas em zona saturada em areia, com interdunas secas subordinadas. Intervalos deformados lateralmente contínuos por centenas de quilômetros ocorrem na zona de contato entre as formações Motuca e Sambaíba. Pelitos com camadas contorcidas e brechadas (Formação Motuca) e arenitos com falhas/microfalhas sinsedimentares, laminação convoluta e diques de injeção preenchidos por argilitos (Formação Sambaíba) são interpretados como sismitos induzidos por terremotos de alta magnitude (>8 na escala Ritcher). Anomalias geoquímicas de elementos traços como Mn, Cr, Co, Cu e Ni na zona de contato entre as formações, juntamente com a presença de micropartículas de composição metálica na matriz argilosa dos sismitos, corroboram com impactos de meteoritos no limite c, possivelmente do astroblema Riachão.

Reconstituição paleoambiental das formações Motuca e Sambaíba, Permo-Triássico da Bacia do Parnaíba no sudoeste do Estado do Maranhão, Brasil

O intervalo que compreende o final do Paleozoico e início do Mesozoico foi marcado por mudanças globais paleogeográficas e paleoclimáticas, em parte atribuídas a eventos catastróficos. A intensa continentalização do supercontinente Pangeia, com a implantação de extensos desertos, sucedeu os ambientes costeiro-plataformais do início do Permiano. Os registros desses eventos no norte do Brasil são encontrados nas bacias intracratônicas, particularmente na sucessão Permotriássica da Bacia do Parnaíba. A análise de fácies e estratigráfica de afloramentos desta sucessão permitiu a individualização de 14 fácies sedimentares agrupadas em 4 associações de fácies (AF): AF1 e AF2, relacionadas aos depósitos da Formação Motuca, e AF3 e AF4, representativas da base da Formação Sambaíba. A AF1 - Lacustre raso/Mudflat consiste em pelitos vermelhos laminados com lentes de gipsita, calcita e marga, além de lobos de arenitos sigmoidais. A AF2 - Saline pan é constituída por corpos lenticulares de gipso laminado, gipso nodular e gipsarenito, sobrepostos por pelitos esverdeados com nódulos de dolomita e palygorskita. A AF3 - Lençol de areia e AF4 - Campo de dunas são formadas, respectivamente, por arenitos de coloração creme alaranjada com estratificação plano-paralela e estratificação cruzada de médio a grande porte. Destaca-se o registro de intervalos deformados lateralmente contínuos por centenas de quilômetros na zona de contato entre as formações Motuca e Sambaíba. Nestes, ocorrem pelitos com camadas contorcidas e brechadas (Formação Motuca) e arenitos com falhas/microfalhas sinsedimentares, laminação convoluta e diques de injeção preenchidos por argilitos (Formação Sambaíba), interpretados como sismitos induzidos por terremotos de alta magnitude (> 8 na escala Richter).

Mapeamento de ocorrências de eventos perigosos no município de Santos SP: integração com dados de vulnerabilidade socioambiental

This study is based on the mapping of occurrences associated with social vulnerability and natural risks, which refer to the resilience of populations and territories, regarding natural hazards associated with the functioning of natural systems (eg, earthquakes, flood, mass movements). According to UNISDR (2014), the state of São Paulo is a reference in working with Urban Resilience and Disaster in the Resilient Cities Campaign (2014), considering the high investment in Risk Areas Mapping and Public Education Campaigns implemented by the Civil Defense. Thus, this study aims to mapping the occurrence of events related to Tree Falls, Erosion, Landslides, Irregular Housing, Rocky Blocks Falls, Wall Falls, Unroofing and Irregular Construction, attended by the Civil Defense in the city of Santos, São Paulo State, from 2011 to 2014. Thereafter, correlated analyzes to the Environmental Vulnerabilities were generated. The Environmental Vulnerability databases used in this dissertation compose the results of the CNPQ Project - Environmental Vulnerability Mapping of the State of Sao Paulo - Brazil: a methodological contribution of Freitas (2013) and Bortolettoet al (2014), with information collected from the 2010 Census (IBGE, 2010), on a census sectors scale. The adopted methodological procedure involves document analysis followed by data integration in Geographic Information System, through algorithms analysis and mapscrossing.The results obtained in Maps of Social and Environmental Vulnerability Occurrences presented areas of High or Very High Vulnerability. The main variables obtained with such characteristics are Irregular Housing, Landslides and Rocky Blocks Falls, which was associated with hilly terrain formations, with slopes above 30%. To the areas of Medium, Low and Very Low Vulnerability were associated the variables Tree Falls, Wall Falls, Erosion, Unroofing and Irregular Construction, which are spatially distributed without an...

Vulnerabilidade às inundações em Teresina, Piauí

Pós-graduação em Geografia - IGCE

Vulnerabilidade às inundações em Teresina, Piauí

Pós-graduação em Geografia - IGCE

TDEM survey in an area of seismicity induced by water wells in Parana sedimentary basin, Northern Sao Paulo State, Brazil

This article presents TDEM results from an area with recent induced shallow seismicity. The purpose was to do a geoelectrical mapping of sedimentary and fractured basaltic aquifers for better understanding of the hydrogeologic setting. The study area is in the Parana basin where flood basalts are overlain by sedimentary units near the city of Bebedouro, northern Sao Paulo State, Brazil. 86 TDEM soundings were acquired in an area of 90 km(2) in the Andes and Botafogo study areas. The soundings were chosen next to wells for calibration, and also along profiles crossing the seismically active areas. 1-D interpretation results showed the general geoelectrical stratigraphy of this part of the Parana basin. The upper geoelectrical layer is the shallow sedimentary aquifer (Adamantina formation) with less than 80 m thickness. The second geoelectrical layer contains the upper basalts of the Serra Geral formation at about 60-80 m depths. A saturated fractured basalt zone between 100 and 300 m depths was identifiable on various TDEM soundings. This depth range corresponds to the range of hypocentral depths for more than 3000 micro-earthquakes in this area. The lower basalt layer was estimated to lie between 400 and 650 m depth. The deepest geoelectrical layer detected by various TDEM soundings corresponds to the Botucatu sandstone (Guarani aquifer). Results suggest that the high-discharge wells are located in the fractured zone in the middle basalt of the Serra Geral formation. There is a good correlation between seismically active areas, high discharge wells (>190 m(3)/h), and fracture zones in the middle basalt. The results reinforce the hypothesis that the shallow seismic activity in the Bebedouro region is being triggered by high rates of groundwater withdrawal. (C) 2012 Elsevier B.V. All rights reserved.

Geochemical constraints on blueschist- and amphibolite-facies rocks of the Central Cordillera of Colombia: the Andean Barragan region

Subduction zones are one of the most characteristic features of planet Earth. Convergent plate junctions exert enormous influence on the formation and recycling of continental crust, and they are also responsible for major mineral resources and earthquakes, which are of crucial importance for society. A subduction-related geologic unit containing high-pressure rocks occurs in the Barragan area (Valle del Cauca Department) on the western flank of the Central Cordillera of the Colombian Andes. Blueschists and amphibolites, serpentinized meta-ultramafic rocks, graphite-chlorite-muscovite-quartz schists, protocataclasites, and graphite-chlorite-andalusite-andesine-garnet-muscovite +/- titanite schists are exposed in this region. In spite of the petrotectonic importance of blueschists, the high-pressure metamorphism of the Central Cordillera of Colombia has been rarely studied. New geochemical data indicate that protoliths of the blueschist- and amphibolite-facies rocks possessed normal mid-ocean ridge basalt bulk compositions. Ar-40/Ar-39 geochronology for a metapelite rock associated with the blueschists shows a plateau age of similar to 120 million years. We suggest that high-P/T conditions were present from similar to 150 to 125 Ma, depending on the model of generation and exhumation considered.

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.

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 use of remote sensing imagery for monitoring recovery in the aftermath of a natural disaster: the Hurricane Katrina Case

Although Recovery is often defined as the less studied and documented phase of the Emergency Management Cycle, a wide literature is available for describing characteristics and sub-phases of this process. Previous works do not allow to gain an overall perspective because of a lack of systematic consistent monitoring of recovery utilizing advanced technologies such as remote sensing and GIS technologies. Taking into consideration the key role of Remote Sensing in Response and Damage Assessment, this thesis is aimed to verify the appropriateness of such advanced monitoring techniques to detect recovery advancements over time, with close attention to the main characteristics of the study event: Hurricane Katrina storm surge. Based on multi-source, multi-sensor and multi-temporal data, the post-Katrina recovery was analysed using both a qualitative and a quantitative approach. The first phase was dedicated to the investigation of the relation between urban types, damage and recovery state, referring to geographical and technological parameters. Damage and recovery scales were proposed to review critical observations on remarkable surge- induced effects on various typologies of structures, analyzed at a per-building level. This wide-ranging investigation allowed a new understanding of the distinctive features of the recovery process. A quantitative analysis was employed to develop methodological procedures suited to recognize and monitor distribution, timing and characteristics of recovery activities in the study area. Promising results, gained by applying supervised classification algorithms to detect localization and distribution of blue tarp, have proved that this methodology may help the analyst in the detection and monitoring of recovery activities in areas that have been affected by medium damage. The study found that Mahalanobis Distance was the classifier which provided the most accurate results, in localising blue roofs with 93.7% of blue roof classified correctly and a producer accuracy of 70%. It was seen to be the classifier least sensitive to spectral signature alteration. The application of the dissimilarity textural classification to satellite imagery has demonstrated the suitability of this technique for the detection of debris distribution and for the monitoring of demolition and reconstruction activities in the study area. Linking these geographically extensive techniques with expert per-building interpretation of advanced-technology ground surveys provides a multi-faceted view of the physical recovery process. Remote sensing and GIS technologies combined to advanced ground survey approach provides extremely valuable capability in Recovery activities monitoring and may constitute a technical basis to lead aid organization and local government in the Recovery management.

Shallow and deep deformation in northern Apennines region using seismological data

For its particular position and the complex geological history, the Northern Apennines has been considered as a natural laboratory to apply several kinds of investigations. By the way, it is complicated to joint all the knowledge about the Northern Apennines in a unique picture that explains the structural and geological emplacement that produced it. The main goal of this thesis is to put together all information on the deformation - in the crust and at depth - of this region and to describe a geodynamical model that takes account of it. To do so, we have analyzed the pattern of deformation in the crust and in the mantle. In both cases the deformation has been studied using always information recovered from earthquakes, although using different techniques. In particular the shallower deformation has been studied using seismic moment tensors information. For our purpose we used the methods described in Arvidsson and Ekstrom (1998) that allowing the use in the inversion of surface waves [and not only of the body waves as the Centroid Moment Tensor (Dziewonski et al., 1981) one] allow to determine seismic source parameters for earthquakes with magnitude as small as 4.0. We applied this tool in the Northern Apennines and through this activity we have built up the Italian CMT dataset (Pondrelli et al., 2006) and the pattern of seismic deformation using the Kostrov (1974) method on a regular grid of 0.25 degree cells. We obtained a map of lateral variations of the pattern of seismic deformation on different layers of depth, taking into account the fact that shallow earthquakes (within 15 km of depth) in the region occur everywhere while most of events with a deeper hypocenter (15-40 km) occur only in the outer part of the belt, on the Adriatic side. For the analysis of the deep deformation, i.e. that occurred in the mantle, we used the anisotropy information characterizing the structure below the Northern Apennines. The anisotropy is an earth properties that in the crust is due to the presence of aligned fluid filled cracks or alternating isotropic layers with different elastic properties while in the mantle the most important cause of seismic anisotropy is the lattice preferred orientation (LPO) of the mantle minerals as the olivine. This last is a highly anisotropic mineral and tends to align its fast crystallographic axes (a-axis) parallel to the astenospheric flow as a response to finite strain induced by geodynamic processes. The seismic anisotropy pattern of a region is measured utilizing the shear wave splitting phenomenon (that is the seismological analogue to optical birefringence). Here, to do so, we apply on teleseismic earthquakes recorded on stations located in the study region, the Sileny and Plomerova (1996) approach. The results are analyzed on the basis of their lateral and vertical variations to better define the earth structure beneath Northern Apennines. We find different anisotropic domains, a Tuscany and an Adria one, with a pattern of seismic anisotropy which laterally varies in a similar way respect to the seismic deformation. Moreover, beneath the Adriatic region the distribution of the splitting parameters is so complex to request an appropriate analysis. Therefore we applied on our data the code of Menke and Levin (2003) which allows to look for different models of structures with multilayer anisotropy. We obtained that the structure beneath the Po Plain is probably even more complicated than expected. On the basis of the results obtained for this thesis, added with those from previous works, we suggest that slab roll-back, which created the Apennines and opened the Tyrrhenian Sea, evolved in the north boundary of Northern Apennines in a different way from its southern part. In particular, the trench retreat developed primarily south of our study region, with an eastward roll-back. In the northern portion of the orogen, after a first stage during which the retreat was perpendicular to the trench, it became oblique with respect to the structure.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.

Improved ground motion intensity measures for reliability-based demand analysis of structures

In Performance-Based Earthquake Engineering (PBEE), evaluating the seismic performance (or seismic risk) of a structure at a designed site has gained major attention, especially in the past decade. One of the objectives in PBEE is to quantify the seismic reliability of a structure (due to the future random earthquakes) at a site. For that purpose, Probabilistic Seismic Demand Analysis (PSDA) is utilized as a tool to estimate the Mean Annual Frequency (MAF) of exceeding a specified value of a structural Engineering Demand Parameter (EDP). This dissertation focuses mainly on applying an average of a certain number of spectral acceleration ordinates in a certain interval of periods, Sa,avg (T1,…,Tn), as scalar ground motion Intensity Measure (IM) when assessing the seismic performance of inelastic structures. Since the interval of periods where computing Sa,avg is related to the more or less influence of higher vibration modes on the inelastic response, it is appropriate to speak about improved IMs. The results using these improved IMs are compared with a conventional elastic-based scalar IMs (e.g., pseudo spectral acceleration, Sa ( T(¹)), or peak ground acceleration, PGA) and the advanced inelastic-based scalar IM (i.e., inelastic spectral displacement, Sdi). The advantages of applying improved IMs are: (i ) "computability" of the seismic hazard according to traditional Probabilistic Seismic Hazard Analysis (PSHA), because ground motion prediction models are already available for Sa (Ti), and hence it is possibile to employ existing models to assess hazard in terms of Sa,avg, and (ii ) "efficiency" or smaller variability of structural response, which was minimized to assess the optimal range to compute Sa,avg. More work is needed to assess also "sufficiency" and "scaling robustness" desirable properties, which are disregarded in this dissertation. However, for ordinary records (i.e., with no pulse like effects), using the improved IMs is found to be more accurate than using the elastic- and inelastic-based IMs. For structural demands that are dominated by the first mode of vibration, using Sa,avg can be negligible relative to the conventionally-used Sa (T(¹)) and the advanced Sdi. For structural demands with sign.cant higher-mode contribution, an improved scalar IM that incorporates higher modes needs to be utilized. In order to fully understand the influence of the IM on the seismis risk, a simplified closed-form expression for the probability of exceeding a limit state capacity was chosen as a reliability measure under seismic excitations and implemented for Reinforced Concrete (RC) frame structures. This closed-form expression is partuclarly useful for seismic assessment and design of structures, taking into account the uncertainty in the generic variables, structural "demand" and "capacity" as well as the uncertainty in seismic excitations. The assumed framework employs nonlinear Incremental Dynamic Analysis (IDA) procedures in order to estimate variability in the response of the structure (demand) to seismic excitations, conditioned to IM. The estimation of the seismic risk using the simplified closed-form expression is affected by IM, because the final seismic risk is not constant, but with the same order of magnitude. Possible reasons concern the non-linear model assumed, or the insufficiency of the selected IM. Since it is impossibile to state what is the "real" probability of exceeding a limit state looking the total risk, the only way is represented by the optimization of the desirable properties of an IM.

Substructuring approache in state space models for dynamic system parameters identification

In the recent decade, the request for structural health monitoring expertise increased exponentially in the United States. The aging issues that most of the transportation structures are experiencing can put in serious jeopardy the economic system of a region as well as of a country. At the same time, the monitoring of structures is a central topic of discussion in Europe, where the preservation of historical buildings has been addressed over the last four centuries. More recently, various concerns arose about security performance of civil structures after tragic events such the 9/11 or the 2011 Japan earthquake: engineers looks for a design able to resist exceptional loadings due to earthquakes, hurricanes and terrorist attacks. After events of such a kind, the assessment of the remaining life of the structure is at least as important as the initial performance design. Consequently, it appears very clear that the introduction of reliable and accessible damage assessment techniques is crucial for the localization of issues and for a correct and immediate rehabilitation. The System Identification is a branch of the more general Control Theory. In Civil Engineering, this field addresses the techniques needed to find mechanical characteristics as the stiffness or the mass starting from the signals captured by sensors. The objective of the Dynamic Structural Identification (DSI) is to define, starting from experimental measurements, the modal fundamental parameters of a generic structure in order to characterize, via a mathematical model, the dynamic behavior. The knowledge of these parameters is helpful in the Model Updating procedure, that permits to define corrected theoretical models through experimental validation. The main aim of this technique is to minimize the differences between the theoretical model results and in situ measurements of dynamic data. Therefore, the new model becomes a very effective control practice when it comes to rehabilitation of structures or damage assessment. The instrumentation of a whole structure is an unfeasible procedure sometimes because of the high cost involved or, sometimes, because it’s not possible to physically reach each point of the structure. Therefore, numerous scholars have been trying to address this problem. In general two are the main involved methods. Since the limited number of sensors, in a first case, it’s possible to gather time histories only for some locations, then to move the instruments to another location and replay the procedure. Otherwise, if the number of sensors is enough and the structure does not present a complicate geometry, it’s usually sufficient to detect only the principal first modes. This two problems are well presented in the works of Balsamo [1] for the application to a simple system and Jun [2] for the analysis of system with a limited number of sensors. Once the system identification has been carried, it is possible to access the actual system characteristics. A frequent practice is to create an updated FEM model and assess whether the structure fulfills or not the requested functions. Once again the objective of this work is to present a general methodology to analyze big structure using a limited number of instrumentation and at the same time, obtaining the most information about an identified structure without recalling methodologies of difficult interpretation. A general framework of the state space identification procedure via OKID/ERA algorithm is developed and implemented in Matlab. Then, some simple examples are proposed to highlight the principal characteristics and advantage of this methodology. A new algebraic manipulation for a prolific use of substructuring results is developed and implemented.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.