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.

Oceanic crust and island arc formation in Central Asia during late Neoproterozoic times - evidence from petrological and geochemical studies

Die vorliegende Arbeit behandelt die Entwicklung des 570 Ma alten, neoproterozoischen Agardagh - Tes-Chem Ophioliths (ATCO) in Zentralasien. Dieser Ophiolith liegt südwestlich des Baikalsees (50.5° N, 95° E) und wurde im frühen Stadium der Akkretion des Zentralasiatischen Mobilgürtels auf den nordwestlichen Rand des Tuvinisch-Mongolischen Mikrokontinentes aufgeschoben. Bei dem Zentralasiatische Mobilgürtel handelt es sich um einen riesigen Akkretions-Subduktionskomplex, der heute das größte zusammenhängende Orogen der Erde darstellt. Im Rahmen dieser Arbeit wurden eine Reihe plutonischer und vulkanischer Gesteine, sowie verschiedene Mantelgesteine des ATCO mittels mikroanalytischer und geochemischer Verfahren untersucht (Elektronenstrahlmikrosonde, Ionenstrahlmikrosonde, Spurenelement- und Isotopengeochemie). Die Auswertung dieser Daten ermöglichte die Entwicklung eines geodynamisch-petrologischen Modells zur Entstehung des ATCO. Die vulkanischen Gesteine lassen sich aufgrund ihrer Spurenelement- und Isotopenzusammensetzung in inselbogenbezogene und back-arc Becken bezogene Gesteine (IA-Gesteine und BAB-Gesteine) unterscheiden. Darüber hinaus gibt es eine weitere, nicht eindeutig zuzuordnende Gruppe, die hauptsächlich mafische Gänge umfasst. Der grösste Teil der untersuchen Vulkanite gehört zur Gruppe der IA-Gesteine. Es handelt sich um Al-reiche Basalte und basaltische Andesite, welche aus einem evolvierten Stammmagma mit Mg# 0.60, Cr ~ 180 µg/g und Ni ~ 95 µg/g hauptsächlich durch Klinopyroxenfraktionierung entstanden sind. Das Stammmagma selbst entstand durch Fraktionierung von ca. 12 % Olivin und geringen Anteilen von Cr-Spinell aus einer primären, aus dem Mantel abgeleiteten Schmelze. Die IA-Gesteine haben hohe Konzentrationen an inkompatiblen Spurenelementen (leichte-(L)- Seltenerdelement-(SEE)-Konzentrationen etwa 100-fach chondritisch, chondrit-normierte (La/Yb)c von 14.6 - 5.1), negative Nb-Anomalien (Nb/La = 0.37 - 0.62) und niedrige Zr/Nb Verhältnisse (7 - 14) relativ zu den BAB-Gesteinen. Initiale eNd Werte liegen bei etwa +5.5, initiale Bleiisotopenverhältnisse sind: 206Pb/204Pb = 17.39 - 18.45, 207Pb/204Pb = 15.49 - 15.61, 208Pb/204Pb = 37.06 - 38.05. Die Anreicherung lithophiler inkompatibler Spurenelemente (LILE) in dieser Gruppe ist signifikant (Ba/La = 11 - 130) und zeigt den Einfluss subduzierter Komponenten an. Die BAB-Gesteine repräsentieren Schmelzen, die sehr wahrscheinlich aus der gleichen Mantelquelle wie die IA-Gesteine stammen, aber durch höhere Aufschmelzgrade (8 - 15 %) und ohne den Einfluss subduzierter Komponenten entstanden sind. Sie haben niedrigere Konzentrationen an inkompatiblen Spurenelementen, flache SEE-Muster ((La/Yb)c = 0.6 - 2.4) und höhere initiale eNd Werte zwischen +7.8 und +8.5. Nb Anomalien existieren nicht und Zr/Nb Verhältnisse sind hoch (21 - 48). Um die geochemische Entwicklung der vulkanischen Gesteine des ATCO zu erklären, sind mindestens drei Komponenten erforderlich: (1) eine angereicherte, ozeaninselbasalt-ähnliche Komponente mit hoher Nb Konzentration über ~ 30 µg/g, einem niedrigen Zr/Nb Verhältnis (ca. 6.5), einem niedrigen initialen eNd Wert (um 0), aber mit radiogenen 206Pb/204Pb-, 207Pb/204Pb- und 208Pb/204Pb-Verhältnissen; (2) eine N-MORB ähnliche back-arc Becken Komponente mit flachem SEE-Muster und einem hohen initialen eNd Wert von mindestens +8.5, und (3) eine Inselbogen-Komponente aus einer verarmten Mantelquelle, welche durch die abtauchende Platte geochemisch modifiziert wurde. Die geochemische Entstehung der ATCO Vulkanite lässt sich dann am besten durch eine Kombination aus Quellenkontamination, fraktionierte Kristallisation und Magmenmischung erklären. Geodynamisch gesehen entstand der ATCO sehr wahrscheinlich in einem intraozeanischen Inselbogen - back-arc System. Bei den untersuchten Plutoniten handelt es sich um ultramafische Kumulate (Wehrlite und Pyroxenite) sowie um gabbroische Plutonite (Olivin-Gabbros bis Diorite). Die geochemischen Charakteristika der mafischen Plutonite sind deutlich unterschiedlich zu denen der vulkanischen Gesteine, weshalb sie sehr wahrscheinlich ein späteres Entwicklungsstadium des ATCO repräsentieren. Die Spurenelement-Konzentrationen in den Klinopyroxenen der ultramafischen Kumulate sind extrem niedrig, mit etwa 0.1- bis 1-fach chondritischen SEE-Konzentrationen und mit deutlich LSEE-verarmten Mustern ((La/Yb)c = 0.27 - 0.52). Berechnete Gleichgewichtsschmelzen der ultramafischen Kumulate zeigen grosse Ähnlichkeit zu primären boninitischen Schmelzen. Die primären Magmen waren daher boninitischer Zusammensetzung und entstanden in dem durch vorausgegangene Schmelzprozesse stark verarmten Mantelkeil über einer Subduktionszone. Niedrige Spurenelement-Konzentrationen zeigen einen geringen Einfluss der abtauchenden Platte an. Die Spurenelement-Konzentrationen der Gabbros sind ebenfalls niedrig, mit etwa 0.5 - 10-fach chondritischen SEE-Konzentrationen und mit variablen SEE-Mustern ((La/Yb)c = 0.25 - 2.6). Analog zu den Vulkaniten der IA-Gruppe haben alle Gabbros eine negative Nb-Anomalie mit Nb/La = 0.01 - 0.31. Die initialen eNd Werte der Gabbros variieren zwischen +4.8 und +7.1, mit einem Mittelwert von +5.9, und sind damit identisch mit denen der IA-Vulkanite. Bei den untersuchten Mantelgesteinen handelt es sich um teilweise serpentinisierte Dunite und Harzburgite, die alle durch hohe Mg/Si- und niedrige Al/Si-Verhältnisse gekennzeichnet sind. Dies zeigt einen refraktären Charakter an und steht in guter Übereinstimmung mit den hohen Cr-Zahlen (Cr#) der Spinelle (bis zu Cr# = 0.83), auf deren Basis der Aufschmelzgrad der residuellen Mantelgesteine berechnet wurde. Dieser beträgt etwa 25 %. Die geochemische Zusammensetzung und die petrologischen Daten der Ultramafite und Gabbros lassen sich am besten erklären, wenn man für die Entstehung dieser Gesteine einen zweistufigen Prozess annimmt. In einer ersten Stufe entstanden die ultramafischen Kumulate unter hohem Druck in einer Magmenkammer an der Krustenbasis, hauptsächlich durch Klinopyroxen-Fraktionierung. Bei dieser Magmenkammer handelte es sich um ein offenes System, dem von unten laufend neue Schmelze zugeführt wurde, und aus dem im oberen Bereich evolviertere Schmelzen geringerer Dichte entwichen. Diese evolvierten Schmelzen stiegen in flachere krustale Bereiche auf und bildeten dort meist isolierte Intrusionskörper. Diese Intrusionskörper erstarrten ohne Magmen-Nachschub, weshalb petrographisch sehr unterschiedliche Gesteine entstehen konnten. Eine geochemische Modifikation der abkühlenden Schmelzen erfolgte allerdings durch die Assimilation von Nebengestein. Da innerhalb der Gabbros keine signifikante Variation der initalen eNd Werte existiert, handelte es sich bei dem assimilierten Material hauptsächlich um vulkanische Gesteine des ATCO und nicht um ältere, möglicherweise kontinentale Kruste.

Anatomy, phylogeny, and palaeoecology of the basal thalattosuchians (Mesoeucrocodylia) from the Liassic of Central Europe

In the marine Jurassic deposits of Europe, a group of marine crocodilians, the Thalattosuchia, belongs to the frequently found reptiles. Thalattosuchia are widely spread in Central Europe from the Jurassic to Lower Cretaceous, and some taxa are also distributed worldwide. The task of the work was to examine all taxa known from the Liassic of Europe. The most frequently known taxa Steneosaurus bollensis and Pelagosaurus typus are anatomically revised. New discoveries at the skull of Pelagosaurus typus e.g., the fact of a partly paired frontal are described by means of computed tomography investigations. In addition, juvenile specimens of this taxon are studied in detail for the first time. The rarely occurring taxon Platysuchus multiscrobiculatus is anatomically described in detail for the first time. It shows both in the skull and in the postcranial material morphological differences to Steneosaurus bollensis and Pelagosaurus typus. Thus Pl. multiscrobiculatus possesses, e.g., an ilium with a deeper acetabulum and a femur with a distinctly flexed femoral head. A juvenile specimen of Pl. multiscrobiculatus is now discovered and is described in parts for the first time, too. Furthermore, Steneosaurus gracilirostris and Steneosaurus brevior known from Lower Jurassic deposits of England are examined and in parts revised. In this work, Steneosaurs brevior is discovered with one specimen from the Upper Liassic of Holzmaden, Germany for the first time, and provides new evidence for the palaeobiogeographical distribution of the taxon. Because of the high number of investigated specimens, it is possible to study ontogenetic development from juvenile to adult stage in Steneosaurus bollensis, Pelagosaurus typus, and Platysuchus multiscrobiculatus. Biometric data are collected from thalattosuchians and extant crocodilians (e.g. Gavialis gangeticus) to investigate intraspecific variation, ontogenetic development, and taxa differentiation. The skulls of Platysuchus multiscrobiculatus and Steneosaurus bollensis are reconstructed three-dimensionally as wax models. The skull reconstructions form the basis of the jaw muscle restoration of Steneosaurus bollensis in connection with comparative studies at extant crocodilians. By means of functional morphologic analysis of the jaw musculature, the dentition, and the locomotor system of S. bollensis, possible conclusions are drawn for the prey options and the hunting behaviour. To clarify the relationships within the Thalattosuchia, a computer-based cladistic phylogenetic in-group analyse of 25 Thalattosuchia taxa is performed. For the analysis, following Thalattosuchia taxa are studied likewise at original material for comparisons: Metriorhynchus superciliosus, Metriorhynchus hastifer, Metriorhynchus leedsi, Geosaurus gracilis, Geosaurus giganteus, Teleidosaurus calvadosi, Teleidosaurus gaudryi, Teleosaurus cadomensis, Teleosaurus geoffroyi, Steneosaurus priscus, Steneosaurus edwardsi, Steneosaurus heberti, Steneosaurus leedsi, Steneosaurus boutilieri, Steneosaurus megarhinus, Steneosaurus obtusidens, and Machimosaurus hugii. The phylogenetic in-group analyse based on 115 characters, reveals a sister-group relationship of the monophyletic Teleosauridae and monophyletic Metriorhynchidae. Within the groups, some taxa are probably paraphyletic. The taxon Pelagosaurus typus is nested inside the Teleosauridae and not outside or within the Metriorhynchidae, as many authors suggested it so far. Based on these results, a tentative palaeobiogeographical-evolutionary scenario is developed.

Spectral-Element and Adjoint 3D Full-Wave Tomography for the Lithosphere of Central Italy

The primary objective of this thesis is to obtain a better understanding of the 3D velocity structure of the lithosphere in central Italy. To this end, I adopted the Spectral-Element Method to perform accurate numerical simulations of the complex wavefields generated by the 2009 Mw 6.3 L’Aquila event and by its foreshocks and aftershocks together with some additional events within our target region. For the mainshock, the source was represented by a finite fault and different models for central Italy, both 1D and 3D, were tested. Surface topography, attenuation and Moho discontinuity were also accounted for. Three-component synthetic waveforms were compared to the corresponding recorded data. The results of these analyses show that 3D models, including all the known structural heterogeneities in the region, are essential to accurately reproduce waveform propagation. They allow to capture features of the seismograms, mainly related to topography or to low wavespeed areas, and, combined with a finite fault model, result into a favorable match between data and synthetics for frequencies up to ~0.5 Hz. We also obtained peak ground velocity maps, that provide valuable information for seismic hazard assessment. The remaining differences between data and synthetics led us to take advantage of SEM combined with an adjoint method to iteratively improve the available 3D structure model for central Italy. A total of 63 events and 52 stations in the region were considered. We performed five iterations of the tomographic inversion, by calculating the misfit function gradient - necessary for the model update - from adjoint sensitivity kernels, constructed using only two simulations for each event. Our last updated model features a reduced traveltime misfit function and improved agreement between data and synthetics, although further iterations, as well as refined source solutions, are necessary to obtain a new reference 3D model for central Italy tomography.

Integrating new and traditional approaches for the estimate of slip-rates of active faults: examples from the Mw 6.3, 2009 L’Aquila earthquake area, Central Italy

This thesis is based on the integration of traditional and innovative approaches aimed at improving the normal faults seimogenic identification and characterization, focusing mainly on slip-rate estimate as a measure of the fault activity. The L’Aquila Mw 6.3 April 6, 2009 earthquake causative fault, namely the Paganica - San Demetrio fault system (PSDFS), was used as a test site. We developed a multidisciplinary and scale‐based strategy consisting of paleoseismological investigations, detailed geomorphological and geological field studies, as well as shallow geophysical imaging and an innovative application of physical properties measurements. We produced a detailed geomorphological and geological map of the PSDFS, defining its tectonic style, arrangement, kinematics, extent, geometry and internal complexities. The PSDFS is a 19 km-long tectonic structure, characterized by a complex structural setting and arranged in two main sectors: the Paganica sector to the NW, characterized by a narrow deformation zone, and the San Demetrio sector to SE, where the strain is accommodated by several tectonic structures, exhuming and dissecting a wide Quaternary basin, suggesting the occurrence of strain migration through time. The integration of all the fault displacement data and age constraints (radiocarbon dating, optically stimulated luminescence (OSL) and tephrochronology) helped in calculating an average Quaternary slip-rate representative for the PSDFS of 0.27 - 0.48 mm/yr. On the basis of its length (ca. 20 km) and slip per event (up to 0.8 m) we also estimated a max expected Magnitude of 6.3-6.8 for this fault. All these topics have a significant implication in terms of surface faulting hazard in the area and may contribute also to the understanding of the PSDFS seismic behavior and of the local seismic hazard.

Kinematics of the Sicily and Calabria Subduction System from Elastic Block Modeling of GPS Data

We use data from about 700 GPS stations in the EuroMediterranen region to investigate the present-day behavior of the the Calabrian subduction zone within the Mediterranean-scale plates kinematics and to perform local scale studies about the strain accumulation on active structures. We focus attenction on the Messina Straits and Crati Valley faults where GPS data show extentional velocity gradients of ∼3 mm/yr and ∼2 mm/yr, respectively. We use dislocation model and a non-linear constrained optimization algorithm to invert for fault geometric parameters and slip-rates and evaluate the associated uncertainties adopting a bootstrap approach. Our analysis suggest the presence of two partially locked normal faults. To investigate the impact of elastic strain contributes from other nearby active faults onto the observed velocity gradient we use a block modeling approach. Our models show that the inferred slip-rates on the two analyzed structures are strongly impacted by the assumed locking width of the Calabrian subduction thrust. In order to frame the observed local deformation features within the present- day central Mediterranean kinematics we realyze a statistical analysis testing the indipendent motion (w.r.t. the African and Eurasias plates) of the Adriatic, Cal- abrian and Sicilian blocks. Our preferred model confirms a microplate like behaviour for all the investigated blocks. Within these kinematic boundary conditions we fur- ther investigate the Calabrian Slab interface geometry using a combined approach of block modeling and χ2ν statistic. Almost no information is obtained using only the horizontal GPS velocities that prove to be a not sufficient dataset for a multi-parametric inversion approach. Trying to stronger constrain the slab geometry we estimate the predicted vertical velocities performing suites of forward models of elastic dislocations varying the fault locking depth. Comparison with the observed field suggest a maximum resolved locking depth of 25 km.

Indoor and Outdoor Natural Radioactivity in the Vulsini Volcanic District (Central Italy): Estimation of Doses and Radiological Risks

Terrestrial radioactivity for most individual is the major contributor to the total dose and is mostly provided by 238U, 232Th and 40K radionuclides. In particular indoor radioactivity is principally due to 222Rn, a radioactive noble gas descendent of 238U, second cause of lung cancer after cigarettes smoking. Vulsini Volcanic District is a well known quaternary volcanic area located between the northern Latium and southern Tuscany (Central Italy). It is characterized by an high natural radiation background resulting from the high concentrations of 238U, 232Th and 40K in the volcanic products. In this context, subduction-related metasomatic enrichment of incompatible elements in the mantle source coupled with magma differentiation within the upper crust has given rise to U, Th and K enriched melts. Almost every ancient village and town located in this part of Italy has been built with volcanic rocks pertaining to the Vulsini Volcanic District. The radiological risk of living in this area has been estimated considering separately: a. the risk associated with buildings made of volcanic products and built on volcanic rock substrates b. the risk associated to soil characteristics. The former has been evaluated both using direct 222Rn indoor measurements and simulations of “standard rooms” built with the tuffs and lavas from the Vulsini Volcanic District investigated in this work. The latter has been carried out by using in situ measurements of 222Rn activity in the soil gases. A radon risk map for the Bolsena village has been developed using soil radon measurements integrating geological information. Data of airborne radioactivity in ambient aerosol at two elevated stations in Emilia Romagna (North Italy) under the influence of Fukushima plume have been collected, effective doses have been calculated and an extensive comparison between doses associated with artificial and natural sources in different area have been described and discussed.