Materiali Ceramici per Edilizia con Funzionalità Fotocatalitica

This thesis wad aimed at the study and application of titanium dioxide photocatalytic activity on ceramic materials. As a matter of fact, photocatalysis is a very promising method to face most of the problems connected with the increasing environmental pollution. Furthermore, titanium dioxide, in its anatase crystallographic phase, is one of the most investigated photocatalytic material and results to be perfectly compatible with silicate body mixes. That goal was pursued by two different strategies: 1. the addition to a body mix used for heavy clay products of several titania powders, with different mean crystallite size, surface area, morphology and anatase/rutile ratio and a titania nanosuspension as well. The titania addition followed two procedures: bulk and spray addition over the ceramic samples surface. Titania was added in two different percentages: 2.5 and 7.5 wt.% in both of the methods. The ceramic samples were then fired at three maximum temperatures: 900, 950 and 1000 °C. Afterwards, the photocatalytic activity of the prepared ceramic samples was evaluated by following the degradation of an organic compound in aqueous medium, under UV radiation. The influence of titania morphological characteristics on the photoactivity of the fired materials was studied by means of XRD and SEM observations. The ceramic samples, sprayed with a slip containing 7.5 wt.% of titania powder and fired at 900 °C, have the best photoactivity, with a complete photo-decomposition of the organic compound. At 1000 °C no sample acted as a photocatalyst due to the anatase-to-rutile phase transformation and to the reaction between titania and calcium and iron oxides in the raw materials. 2. The second one foresaw the synthesis of TiO2-SiO2 solid solutions, using the following stoichiometry: Ti1-xSixO2 where x = 0, 0.1, 0.3 and 0.5 atoms per formula unit (apfu). The mixtures were then fired following two thermal cycles, each with three maximum temperatures. The effect of SiO2 addition into the TiO2 crystal structure and, consequently, on its photocatalytic activity when fired at high temperature, was thoroughly investigated by means of XRD, XPS, FE-SEM, TEM and BET analysis. The photoactivity of the prepared powders was assessed both in gas and liquid phase. Subsequently, the TiO2-SiO2 solid solutions, previously fired at 900 °C, were sprayed over the ceramic samples surface in the percentage of 7.5 wt.%. The prepared ceramic samples were fired at 900 and 1000 °C. The photocatalytic activity of the ceramic samples was evaluated in liquid phase. Unfortunately, that samples did not show any appreciable photoactivity. In fact, samples fired at 900 °C showed a pretty low photoactivity, while the one fired at 1000 °C showed no photoactivity at all. This was explained by the excessive coarsening of titania particles. To summarise, titania particle size, more than its crystalline phase, seems to have a relevant role in the photocatalytic activity of the ceramic samples.

Rapporti tra tettonica e clima durante l'esumazione nelle Alpi Centrali. Evidenze dalla termocronologia e dall'analisi strutturale lungo il traforo ferroviario del Sempione

(U-Th)/He and fission-track analyses of apatite along deep-seated tunnels crossing high-relief mountain ranges offer the opportunity to investigate climate and tectonic forcing on the topographic evolution. In this study, the thermochronologic analysis of a large set of samples collected in the Simplon railway tunnel (western-central Alps; Italy and Switzerland) and along its surface trace, coupled with kinematic and structural analysis of major fault zones intersecting the tunnel, constrains the phenomena controlling the topographic and structural evolution, during the latest stage of exhumation of the Simplon Massif, and the timing in which they operated. The study area is located at the western margin of the Lepontine metamorphic dome where a complex nappe-stack pertaining to the Penninic and Ultrahelvetic domains experienced a fast exhumation from the latest Oligocene onward. The exhumation was mainly accommodated by a west-dipping low-angle detachment (the Simplon Fault Zone) which is located just 8 km to the west of the tunnel. However, along the section itself several faults related to two principal phases both with important dip-slip kinematics have been detected. Cooling rates derived from our thermocronological data vary from about 10 °C/Ma at about 10 Ma to about 35 °C/Ma in the last 5 Ma. Such increase in the cooling rate corresponds to the most important climatic change recorded in the northern hemisphere in the last 10 Ma, i.e. the shift to wetter conditions at the end of the Messinian salinity crisis and the inception of glacial cycles in the northern hemisphere. In addition, (U-Th)/He and fission-track age patterns lack of important correlation with the topography suggesting that the present-day relief morphology is the result of recent erosional dynamics. More in details, the (U-Th)/He tunnel ages show an impressive uniformity at 2 Ma, whereas cooling rates calculated at 1 Ma increase towards the two major valleys. This indicates a focusing of erosive processes in the valleys which led to the shaping of present-day topography. Structural analysis documents the presence of two phases of brittle deformation postdating the metamorphic phases in the area. The first one is directly related to the last phase of activity along the Simplon Fault Zone and is characterized by extension towards SO and vertical shortening. The young one is characterized by extension towards NO and horizontal shortening in a along the NE-SO direction. Structures related to the first phase of brittle deformation generate important variations in the older ages' dataset, until 3 Ma, suggesting that tectonics controlled rocks exhumation up to that age. Structures related to the second phase generate some variations also in the younger age dataset, highlighting the activity of faults bordering the massif and suggesting a continuous activity also after 2 Ma. However, most of (U-Th)/He tunnel ages, varying slightly around 2 Ma, document that the Simplon area has experienced primarily erosional exhumation in this time span. In conclusion, all our data suggest that in the central Italian Alps the climatic signal gradually overrode the tectonic effects after about 5 Ma, as a consequence of the climatic instability started at end of Messinian salinity crisis and improved by the onset of glaciations in the northern hemisphere.

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Tectonics and kinematics of curved mountain belts: examples from the Andes and the Alps

Curved mountain belts have always fascinated geologists and geophysicists because of their peculiar structural setting and geodynamic mechanisms of formation. The need of studying orogenic bends arises from the numerous questions to which geologists and geophysicists have tried to answer to during the last two decades, such as: what are the mechanisms governing orogenic bends formation? Why do they form? Do they develop in particular geological conditions? And if so, what are the most favorable conditions? What are their relationships with the deformational history of the belt? Why is the shape of arcuate orogens in many parts of the Earth so different? What are the factors controlling the shape of orogenic bends? Paleomagnetism demonstrated to be one of the most effective techniques in order to document the deformation of a curved belt through the determination of vertical axis rotations. In fact, the pattern of rotations within a curved belt can reveal the occurrence of a bending, and its timing. Nevertheless, paleomagnetic data alone are not sufficient to constrain the tectonic evolution of a curved belt. Usually, structural analysis integrates paleomagnetic data, in defining the kinematics of a belt through kinematic indicators on brittle fault planes (i.e., slickensides, mineral fibers growth, SC-structures). My research program has been focused on the study of curved mountain belts through paleomagnetism, in order to define their kinematics, timing, and mechanisms of formation. Structural analysis, performed only in some regions, supported and integrated paleomagnetic data. In particular, three arcuate orogenic systems have been investigated: the Western Alpine Arc (NW Italy), the Bolivian Orocline (Central Andes, NW Argentina), and the Patagonian Orocline (Tierra del Fuego, southern Argentina). The bending of the Western Alpine Arc has been investigated so far using different approaches, though few based on reliable paleomagnetic data. Results from our paleomagnetic study carried out in the Tertiary Piedmont Basin, located on top of Alpine nappes, indicate that the Western Alpine Arc is a primary bend that has been subsequently tightened by further ~50° during Aquitanian-Serravallian times (23-12 Ma). This mid-Miocene oroclinal bending, superimposing onto a pre-existing Eocene nonrotational arc, is the result of a composite geodynamic mechanism, where slab rollback, mantle flows, and rotating thrust emplacement are intimately linked. Relying on our paleomagnetic and structural evidence, the Bolivian Orocline can be considered as a progressive bend, whose formation has been driven by the along-strike gradient of crustal shortening. The documented clockwise rotations up to 45° are compatible with a secondary-bending type mechanism occurring after Eocene-Oligocene times (30-40 Ma), and their nature is probably related to the widespread shearing taking place between zones of differential shortening. Since ~15 Ma ago, the activity of N-S left-lateral strike-slip faults in the Eastern Cordillera at the border with the Altiplano-Puna plateau induced up to ~40° counterclockwise rotations along the fault zone, locally annulling the regional clockwise rotation. We proposed that mid-Miocene strike-slip activity developed in response of a compressive stress (related to body forces) at the plateau margins, caused by the progressive lateral (southward) growth of the Altiplano-Puna plateau, laterally spreading from the overthickened crustal region of the salient apex. The growth of plateaux by lateral spreading seems to be a mechanism common to other major plateaux in the Earth (i.e., Tibetan plateau). Results from the Patagonian Orocline represent the first reliable constraint to the timing of bending in the southern tip of South America. They indicate that the Patagonian Orocline did not undergo any significant rotation since early Eocene times (~50 Ma), implying that it may be considered either a primary bend, or an orocline formed during the late Cretaceous-early Eocene deformation phase. This result has important implications on the opening of the Drake Passage at ~32 Ma, since it is definitely not related to the formation of the Patagonian orocline, but the sole consequence of the Scotia plate spreading. Finally, relying on the results and implications from the study of the Western Alpine Arc, the Bolivian Orocline, and the Patagonian Orocline, general conclusions on curved mountain belt formation have been inferred.

Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

Correlazioni tra proprietá reologiche, struttura e processabilitá di blend di LLDPE/LDPE

Polymer blends constitute a valuable way to produce relatively low cost new materials. A still open question concerns the miscibility of polyethylene blends. Deviations from the log-additivity rule of the newtonian viscosity are often taken as a signature of immiscibility of the two components. The aim of this thesis is to characterize the rheological behavior in shear and elongation of five series of LLDPE/LDPE blends whose parent polymers have been chosen with different viscosity and SCB content and length. Synergistic effects have been measured for both zero shear viscosity and melt strength. Both SCB length and viscosity ratio between the components have been found to be key parameters for the miscibility of the pure polymers. In particular the miscibility increases with increasing SCB length and with decreasing the LDPE molecular weight and viscosity. This rheological behavior has significant effects on the processability window of these blends when the uni or biaxial elongational flows are involved. The film blowing is one of the processes for which the synergistic effects above mentioned can be crucial. Small scale experiments of film blowing performed for one of the series of blends has demonstrated that the positive deviation of the melt strength enlarges the processability window. In particular, the bubble stability was found to improve or disappear when the melt strength of the samples increased. The blending of LDPE and LLDPE can even reduce undesired melt flow instability phenomena widening, as a consequence, the processability window in extrusion. One of the series of blends has been characterized by means of capillary rheometry in order to allow a careful morphological analysis of the surface of the extruded polymer jets by means of Scanning Electron Microscopy (SEM) with the aim to detect the very early stages of the small scale melt instabilty at low shear rates (sharksin) and to follow its subsequent evolution as long as the shear rate was increased. With this experimental procedure it was possible to evaluate the shear rate ranges corresponding to different flow regions: smooth extrudate surface (absence of instability), sharkskin (small scale instability produced at the capillary exit), stick-slip transition (instability involving the whole capillary wall) and gross melt fracture (i.e. a large scale "upstream" instability originating from the entrance region of the capillary). A quantitative map was finally worked out using which an assessment of the flow type for a given shear rate and blend composition can be predicted.

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.

Provenienza dei sedimenti arenitici nel bacino di Tracia (eo-oligocene, Turchia nord-occidentale e Grecia nord-orientale)

The Thrace Basin is the largest and thickest Tertiary sedimentary basin of the eastern Balkans region and constitutes an important hydrocarbon province. It is located between the Rhodope-Strandja Massif to the north and west, the Marmara Sea and Biga Peninsula to the south, and the Black Sea to the est. It consists of a complex system of depocenters and uplifts with very articulate paleotopography indicated by abrupt lateral facies variations. Its southeastern margin is widely deformed by the Ganos Fault, a segment of the North Anatolian strike-slip fault system . Most of the Thrace Basin fill ranges from the Eocene to the Late Oligocene. Maximum total thickness, including the Neogene-Quaternary succession, reaches 9.000 meters in a few narrow depocenters. This sedimentary succession consists mainly of basin plain turbiditic deposits with a significant volcaniclastic component which evolves upwards to shelf deposits and continental facies, with deltaic bodies prograding towards the basin center in the Oligocene. This work deals with the provenance of Eocene-Oligocene clastic sediments of the southern and western part of Thrace Basin in Turkey and Greece. Sandstone compositional data (78 gross composition analyses and 40 heavy minerals analyses) were used to understand the change in detrital modes which reflects the provenance and geodinamic evolution of the basin. Samples were collected at six localities, which are from west to est: Gökçeada, Gallipoli and South-Ganos (south of Ganos Fault), Alexandroupolis, Korudağ and North-Ganos (north of Ganos Fault). Petrologic (framework composition and heavy-mineral analyses) and stratigraphic-sedimentologic data, (analysis of sedimentologic facies associations along representative stratigraphic sections, paleocurrents) allowed discrimination of six petrofacies; for each petrofacies the sediment dispersal system was delineated. The Thrace Basin fill is made mainly of lithic arkoses and arkosic litharenites with variable amount of low-grade metamorphic lithics (also ophiolitic), neovolcanic lithics, and carbonate grains (mainly extrabasinal). Picotite is the most widespread heavy mineral in all petrofacies. Petrological data on analyzed successions show a complex sediment dispersal pattern and evolution of the basin, indicating one principal detrital input from a source area located to the south, along both the İzmir-Ankara and Intra-Pontide suture lines, and a possible secondary source area, represented by the Rhodope Massif to the west. A significant portion of the Thrace Basin sediments in the study area were derived from ophiolitic source rocks and from their oceanic cover, whereas epimetamorphic detrital components came from a low-grade crystalline basement. An important penecontemporaneous volcanic component is widespread in late Eocene-Oligocene times, indicating widespread post-collisional (collapse?) volcanism following the closure of the Vardar ocean. Large-scale sediment mass wasting from south to north along the southern margin of the Thrace Basin is indicated (i) in late Eocene time by large olistoliths of ophiolites and penecontemporaneous carbonates, and (ii) in the mid-Oligocene by large volcaniclastic olistoliths. The late Oligocene paleogeographic scenario was characterized by large deltaic bodies prograding northward (Osmancik Formation). This clearly indicates that the southern margin of the basin acted as a major sediment source area throughout its Eocene-Oligocene history. Another major sediment source area is represented by the Rhodope Massif, in particolar the Circum-Rhodopic belt, especially for plutonic and metamorphic rocks. Considering preexisting data on the petrologic composition of Thrace Basin, silicilastic sediments in Greece and Bulgaria (Caracciolo, 2009), a Rhodopian provenance could be considered mostly for areas of the Thrace Basin outside our study area, particularly in the northern-central portions of the basin. In summary, the most important source area for the sediment of Thrace Basin in the study area was represented by the exhumed subduction-accretion complex along the southern margin of the basin (Biga Peninsula and western-central Marmara Sea region). Most measured paleocurrent indicators show an eastward paleoflow but this is most likely the result of gravity flow deflection. This is possible considered a strong control due to the east-west-trending synsedimentary transcurrent faults which cuts the Thrace Basin, generating a series of depocenters and uplifts which deeply influenced sediment dispersal and the areal distribution of paleoenvironments. The Thrace Basin was long interpreted as a forearc basin between a magmatic arc to the north and a subduction-accretion complex to the south, developed in a context of northward subduction. This interpretation was challenged by more recent data emphasizing the lack of a coeval magmatic arc in the north and the interpretation of the chaotic deposit which outcrop south of Ganos Fault as olistoliths and large submarine slumps, derived from the erosion and sedimentary reworking of an older mélange unit located to the south (not as tectonic mélange formed in an accretionary prism). The present study corroborates instead the hypothesis of a post-collisional origin of the Thrace Basin, due to a phase of orogenic collapse, which generated a series of mid-Eocene depocenters all along the İzmir-Ankara suture (following closure of the Vardar-İzmir-Ankara ocean and the ensuing collision); then the slab roll-back of the remnant Pindos ocean played an important role in enhancing subsidence and creating additional accommodation space for sediment deposition.

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.

Detailed study of the seismotectonic setting in the Lucanian Apennines and surrounding areas (Southern Italy)

In this research work I analyzed the instrumental seismicity of Southern Italy in the area including the Lucanian Apennines and Bradano foredeep, making use of the most recent seismological database available so far. I examined the seismicity occurred during the period between 2001 and 2006, considering 514 events with magnitudes M ≥ 2.0. In the first part of the work, P- and S-wave arrival times, recorded by the Italian National Seismic Network (RSNC) operated by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), were re-picked along with those of the SAPTEX temporary array (2001–2004). For some events located in the Upper Val d'Agri, I also used data from the Eni-Agip oil company seismic network. I computed the VP/VS ratio obtaining a value of 1.83 and I carried out an analysis for the one-dimensional (1D) velocity model that approximates the seismic structure of the study area. After this preliminary analysis, making use of the records obtained in the SeSCAL experiment, I incremented the database by handpicking new arrival times. My final dataset consists of 15,666 P- and 9228 S-arrival times associated to 1047 earthquakes with magnitude ML ≥ 1.5. I computed 162 fault-plane solutions and composite focal mechanisms for closely located events. I investigated stress field orientation inverting focal mechanism belonging to the Lucanian Apennine and the Pollino Range, both areas characterized by more concentrated background seismicity. Moreover, I applied the double difference technique (DD) to improve the earthquake locations. Considering these results and different datasets available in the literature, I carried out a detailed analysis of single sub-areas and of a swarm (November 2008) recorded by SeSCAL array. The relocated seismicity appears more concentrated within the upper crust and it is mostly clustered along the Lucanian Apennine chain. In particular, two well-defined clusters were located in the Potentino and in the Abriola-Pietrapertosa sector (central Lucanian region). Their hypocentral depths are slightly deeper than those observed beneath the chain. I suggest that these two seismic features are representative of the transition from the inner portion of the chain with NE-SW extension to the external margin characterized by dextral strike-slip kinematics. In the easternmost part of the study area, below the Bradano foredeep and the Apulia foreland, the seismicity is generally deeper and more scattered and is associated to the Murge uplift and to the small structures present in the area. I also observed a small structure NE-SW oriented in the Abriola-Pietrapertosa area (activated with a swarm in November 2008) that could be considered to act as a barrier to the propagation of a potential rupture of an active NW-SE striking faults system. Focal mechanisms computed in this study are in large part normal and strike-slip solutions and their tensional axes (T-axes) have a generalized NE-SW orientation. Thanks to denser coverage of seismic stations and the detailed analysis, this study is a further contribution to the comprehension of the seismogenesis and state of stress of the Southern Apennines region, giving important contributions to seismotectonic zoning and seismic hazard assessment.

Kinematic description of the rupture from strong motion data: strategies for a robust inversion

We present a non linear technique to invert strong motion records with the aim of obtaining the final slip and rupture velocity distributions on the fault plane. In this thesis, the ground motion simulation is obtained evaluating the representation integral in the frequency. The Green’s tractions are computed using the discrete wave-number integration technique that provides the full wave-field in a 1D layered propagation medium. The representation integral is computed through a finite elements technique, based on a Delaunay’s triangulation on the fault plane. The rupture velocity is defined on a coarser regular grid and rupture times are computed by integration of the eikonal equation. For the inversion, the slip distribution is parameterized by 2D overlapping Gaussian functions, which can easily relate the spectrum of the possible solutions with the minimum resolvable wavelength, related to source-station distribution and data processing. The inverse problem is solved by a two-step procedure aimed at separating the computation of the rupture velocity from the evaluation of the slip distribution, the latter being a linear problem, when the rupture velocity is fixed. The non-linear step is solved by optimization of an L2 misfit function between synthetic and real seismograms, and solution is searched by the use of the Neighbourhood Algorithm. The conjugate gradient method is used to solve the linear step instead. The developed methodology has been applied to the M7.2, Iwate Nairiku Miyagi, Japan, earthquake. The estimated magnitude seismic moment is 2.6326 dyne∙cm that corresponds to a moment magnitude MW 6.9 while the mean the rupture velocity is 2.0 km/s. A large slip patch extends from the hypocenter to the southern shallow part of the fault plane. A second relatively large slip patch is found in the northern shallow part. Finally, we gave a quantitative estimation of errors associates with the parameters.

Finite element models of coseismic deformation due to the 2009 L'Aquila (Italy) and 2008 Wenchuan(China) earthquakes

The topic of my Ph.D. thesis is the finite element modeling of coseismic deformation imaged by DInSAR and GPS data. I developed a method to calculate synthetic Green functions with finite element models (FEMs) and then use linear inversion methods to determine the slip distribution on the fault plane. The method is applied to the 2009 L’Aquila Earthquake (Italy) and to the 2008 Wenchuan earthquake (China). I focus on the influence of rheological features of the earth's crust by implementing seismic tomographic data and the influence of topography by implementing Digital Elevation Models (DEM) layers on the FEMs. Results for the L’Aquila earthquake highlight the non-negligible influence of the medium structure: homogeneous and heterogeneous models show discrepancies up to 20% in the fault slip distribution values. Furthermore, in the heterogeneous models a new area of slip appears above the hypocenter. Regarding the 2008 Wenchuan earthquake, the very steep topographic relief of Longmen Shan Range is implemented in my FE model. A large number of DEM layers corresponding to East China is used to achieve the complete coverage of the FE model. My objective was to explore the influence of the topography on the retrieved coseismic slip distribution. The inversion results reveals significant differences between the flat and topographic model. Thus, the flat models frequently adopted are inappropriate to represent the earth surface topographic features and especially in the case of the 2008 Wenchuan earthquake.

Magma injections and destabilization of basaltic volcanoes: a numerical study. Application to La Reunion (Indian ocean, France) and Stromboli (Tyrrhenian sea, Italy)

Most basaltic volcanoes are affected by recurrent lateral instabilities during their evolution. Numerous factors have been shown to be involved in the process of flank destabilization occurring over long periods of time or by instantaneous failures. However, the role of these factors on the mechanical behaviour and stability of volcanic edifices is poorly-constrained as lateral failure usually results from the combined effects of several parameters. Our study focuses on the morphological and structural comparison of two end-member basaltic systems, La Reunion (Indian ocean, France) and Stromboli (southern Tyrrhenian sea, Italy). We showed that despite major differences on their volumes and geodynamic settings, both systems present some similarities as they are characterized by an intense intrusive activity along well-developed rift zones and recurrent phenomena of flank collapse during their evolution. Among the factors of instability, the examples of la Reunion and Stromboli evidence the major contribution of intrusive complexes to volcano growth and destruction as attested by field observations and the monitoring of these active volcanoes. Classical models consider the relationship between vertical intrusions of magma and flank movements along a preexisting sliding surface. A set of published and new field data from Piton des Neiges volcano (La Reunion) allowed us to recognize the role of subhorizontal intrusions in the process of flank instability and to characterize the geometry of both subvertical and subhorizontal intrusions within basaltic edifices. This study compares the results of numerical modelling of the displacements associated with high-angle and low-angle intrusions within basaltic volcanoes. We use a Mixed Boundary Element Method to investigate the mechanical response of an edifice to the injection of magmatic intrusions in different stress fields. Our results indicate that the anisotropy of the stress field favours the slip along the intrusions due to cointrusive shear stress, generating flank-scale displacements of the edifice, especially in the case of subhorizontal intrusions, capable of triggering large-scale flank collapses on basaltic volcanoes. Applications of our theoretical results to real cases of flank displacements on basaltic volcanoes (such as the 2007 eruptive crisis at La Reunion and Stromboli) revealed that the previous model of subvertical intrusions-related collapse is a likely mechanism affecting small-scale steeply-sloping basaltic volcanoes like Stromboli. Furthermore, our field study combined to modelling results confirms the importance of shallow-dipping intrusions in the morpho-structural evolution of large gently-sloping basaltic volcanoes like Piton de la Fournaise, Etna and Kilauea, with particular regards to flank instability, which can cause catastrophic tsunamis.

InSAR deformation measurements of the earthquake cycle in transcurrent tectonic domains, analytical and analog modeling

I applied the SBAS-DInSAR method to the Mattinata Fault (MF) (Southern Italy) and to the Doruneh Fault System (DFS) (Central Iran). In the first case, I observed limited internal deformation and determined the right lateral kinematic pattern with a compressional pattern in the northern sector of the fault. Using the Okada model I inverted the observed velocities defining a right lateral strike slip solution for the MF. Even if it fits the data within the uncertainties, the modeled slip rate of 13-15 mm yr-1 seems too high with respect to the geological record. Concerning the Western termination of DFS, SAR data confirms the main left lateral transcurrent kinematics of this fault segment, but reveal a compressional component. My analytical model fits successfully the observed data and quantifies the slip in ~4 mm yr-1 and ~2.5 mm yr-1 of pure horizontal and vertical displacement respectively. The horizontal velocity is compatible with geological record. I applied classic SAR interferometry to the October–December 2008 Balochistan (Central Pakistan) seismic swarm; I discerned the different contributions of the three Mw > 5.7 earthquakes determining fault positions, lengths, widths, depths and slip distributions, constraining the other source parameters using different Global CMT solutions. A well constrained solution has been obtained for the 09/12/2008 aftershock, whereas I tested two possible fault solutions for the 28-29/10/08 mainshocks. It is not possible to favor one of the solutions without independent constraints derived from geological data. Finally I approached the study of the earthquake-cycle in transcurrent tectonic domains using analog modeling, with alimentary gelatins like crust analog material. I successfully joined the study of finite deformation with the earthquake cycle study and sudden dislocation. A lot of seismic cycles were reproduced in which a characteristic earthquake is recognizable in terms of displacement, coseismic velocity and recurrence time.

Fault delineation and stress orientations from the analysis of background, low magnitude seismicity in Southern Apennines (Italy)

The aim of this work was to show that refined analyses of background, low magnitude seismicity allow to delineate the main active faults and to accurately estimate the directions of the regional tectonic stress that characterize the Southern Apennines (Italy), a structurally complex area with high seismic potential. Thanks the presence in the area of an integrated dense and wide dynamic network, was possible to analyzed an high quality microearthquake data-set consisting of 1312 events that occurred from August 2005 to April 2011 by integrating the data recorded at 42 seismic stations of various networks. The refined seismicity location and focal mechanisms well delineate a system of NW-SE striking normal faults along the Apenninic chain and an approximately E-W oriented, strike-slip fault, transversely cutting the belt. The seismicity along the chain does not occur on a single fault but in a volume, delimited by the faults activated during the 1980 Irpinia M 6.9 earthquake, on sub-parallel predominant normal faults. Results show that the recent low magnitude earthquakes belongs to the background seismicity and they are likely generated along the major fault segments activated during the most recent earthquakes, suggesting that they are still active today thirty years after the mainshock occurrences. In this sense, this study gives a new perspective to the application of the high quality records of low magnitude background seismicity for the identification and characterization of active fault systems. The analysis of the stress tensor inversion provides two equivalent models to explain the microearthquake generation along both the NW-SE striking normal faults and the E- W oriented fault with a dominant dextral strike-slip motion, but having different geological interpretations. We suggest that the NW-SE-striking Africa-Eurasia convergence acts in the background of all these structures, playing a primary and unifying role in the seismotectonics of the whole region.