A comprehensive study of the 26th December 2004 Sumatra earthquake: possible implications on Earth rotation and investigations on the coseismic and postseismic stress diffusion associated with the seismic rupture

In this work a multidisciplinary study of the December 26th, 2004 Sumatra earthquake has been carried out. We have investigated both the effect of the earthquake on the Earth rotation and the stress field variations associated with the seismic event. In the first part of the work we have quantified the effects of a water mass redistribution associated with the propagation of a tsunami wave on the Earth’s pole path and on the length-of-day (LOD) and applied our modeling results to the tsunami following the 2004 giant Sumatra earthquake. We compared the result of our simulations on the instantaneous rotational axis variations with some preliminary instrumental evidences on the pole path perturbation (which has not been confirmed yet) registered just after the occurrence of the earthquake, which showed a step-like discontinuity that cannot be attributed to the effect of a seismic dislocation. Our results show that the perturbation induced by the tsunami on the instantaneous rotational pole is characterized by a step-like discontinuity, which is compatible with the observations but its magnitude turns out to be almost one hundred times smaller than the detected one. The LOD variation induced by the water mass redistribution turns out to be not significant because the total effect is smaller than current measurements uncertainties. In the second part of this work of thesis we modeled the coseismic and postseismic stress evolution following the Sumatra earthquake. By means of a semi-analytical, viscoelastic, spherical model of global postseismic deformation and a numerical finite-element approach, we performed an analysis of the stress diffusion following the earthquake in the near and far field of the mainshock source. We evaluated the stress changes due to the Sumatra earthquake by projecting the Coulomb stress over the sequence of aftershocks taken from various catalogues in a time window spanning about two years and finally analyzed the spatio-temporal pattern. The analysis performed with the semi-analytical and the finite-element modeling gives a complex picture of the stress diffusion, in the area under study, after the Sumatra earthquake. We believe that the results obtained with the analytical method suffer heavily for the restrictions imposed, on the hypocentral depths of the aftershocks, in order to obtain the convergence of the harmonic series of the stress components. On the contrary we imposed no constraints on the numerical method so we expect that the results obtained give a more realistic description of the stress variations pattern.

Estimating crustal deformation from geodetic data

By the end of the 19th century, geodesy has contributed greatly to the knowledge of regional tectonics and fault movement through its ability to measure, at sub-centimetre precision, the relative positions of points on the Earth’s surface. Nowadays the systematic analysis of geodetic measurements in active deformation regions represents therefore one of the most important tool in the study of crustal deformation over different temporal scales [e.g., Dixon, 1991]. This dissertation focuses on motion that can be observed geodetically with classical terrestrial position measurements, particularly triangulation and leveling observations. The work is divided into two sections: an overview of the principal methods for estimating longterm accumulation of elastic strain from terrestrial observations, and an overview of the principal methods for rigorously inverting surface coseismic deformation fields for source geometry with tests on synthetic deformation data sets and applications in two different tectonically active regions of the Italian peninsula. For the long-term accumulation of elastic strain analysis, triangulation data were available from a geodetic network across the Messina Straits area (southern Italy) for the period 1971 – 2004. From resulting angle changes, the shear strain rates as well as the orientation of the principal axes of the strain rate tensor were estimated. The computed average annual shear strain rates for the time period between 1971 and 2004 are γ˙1 = 113.89 ± 54.96 nanostrain/yr and γ˙2 = -23.38 ± 48.71 nanostrain/yr, with the orientation of the most extensional strain (θ) at N140.80° ± 19.55°E. These results suggests that the first-order strain field of the area is dominated by extension in the direction perpendicular to the trend of the Straits, sustaining the hypothesis that the Messina Straits could represents an area of active concentrated deformation. The orientation of θ agree well with GPS deformation estimates, calculated over shorter time interval, and is consistent with previous preliminary GPS estimates [D’Agostino and Selvaggi, 2004; Serpelloni et al., 2005] and is also similar to the direction of the 1908 (MW 7.1) earthquake slip vector [e.g., Boschi et al., 1989; Valensise and Pantosti, 1992; Pino et al., 2000; Amoruso et al., 2002]. Thus, the measured strain rate can be attributed to an active extension across the Messina Straits, corresponding to a relative extension rate ranges between < 1mm/yr and up to ~ 2 mm/yr, within the portion of the Straits covered by the triangulation network. These results are consistent with the hypothesis that the Messina Straits is an important active geological boundary between the Sicilian and the Calabrian domains and support previous preliminary GPS-based estimates of strain rates across the Straits, which show that the active deformation is distributed along a greater area. Finally, the preliminary dislocation modelling has shown that, although the current geodetic measurements do not resolve the geometry of the dislocation models, they solve well the rate of interseismic strain accumulation across the Messina Straits and give useful information about the locking the depth of the shear zone. Geodetic data, triangulation and leveling measurements of the 1976 Friuli (NE Italy) earthquake, were available for the inversion of coseismic source parameters. From observed angle and elevation changes, the source parameters of the seismic sequence were estimated in a join inversion using an algorithm called “simulated annealing”. The computed optimal uniform–slip elastic dislocation model consists of a 30° north-dipping shallow (depth 1.30 ± 0.75 km) fault plane with azimuth of 273° and accommodating reverse dextral slip of about 1.8 m. The hypocentral location and inferred fault plane of the main event are then consistent with the activation of Periadriatic overthrusts or other related thrust faults as the Gemona- Kobarid thrust. Then, the geodetic data set exclude the source solution of Aoudia et al. [2000], Peruzza et al. [2002] and Poli et al. [2002] that considers the Susans-Tricesimo thrust as the May 6 event. The best-fit source model is then more consistent with the solution of Pondrelli et al. [2001], which proposed the activation of other thrusts located more to the North of the Susans-Tricesimo thrust, probably on Periadriatic related thrust faults. The main characteristics of the leveling and triangulation data are then fit by the optimal single fault model, that is, these results are consistent with a first-order rupture process characterized by a progressive rupture of a single fault system. A single uniform-slip fault model seems to not reproduce some minor complexities of the observations, and some residual signals that are not modelled by the optimal single-fault plane solution, were observed. In fact, the single fault plane model does not reproduce some minor features of the leveling deformation field along the route 36 south of the main uplift peak, that is, a second fault seems to be necessary to reproduce these residual signals. By assuming movements along some mapped thrust located southward of the inferred optimal single-plane solution, the residual signal has been successfully modelled. In summary, the inversion results presented in this Thesis, are consistent with the activation of some Periadriatic related thrust for the main events of the sequence, and with a minor importance of the southward thrust systems of the middle Tagliamento plain.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Imaging of crustal scatterers using multiple seismic arrays

Array seismology is an useful tool to perform a detailed investigation of the Earth’s interior. Seismic arrays by using the coherence properties of the wavefield are able to extract directivity information and to increase the ratio of the coherent signal amplitude relative to the amplitude of incoherent noise. The Double Beam Method (DBM), developed by Krüger et al. (1993, 1996), is one of the possible applications to perform a refined seismic investigation of the crust and mantle by using seismic arrays. The DBM is based on a combination of source and receiver arrays leading to a further improvement of the signal-to-noise ratio by reducing the error in the location of coherent phases. Previous DBM works have been performed for mantle and core/mantle resolution (Krüger et al., 1993; Scherbaum et al., 1997; Krüger et al., 2001). An implementation of the DBM has been presented at 2D large-scale (Italian data-set for Mw=9.3, Sumatra earthquake) and at 3D crustal-scale as proposed by Rietbrock & Scherbaum (1999), by applying the revised version of Source Scanning Algorithm (SSA; Kao & Shan, 2004). In the 2D application, the rupture front propagation in time has been computed. In 3D application, the study area (20x20x33 km3), the data-set and the source-receiver configurations are related to the KTB-1994 seismic experiment (Jost et al., 1998). We used 60 short-period seismic stations (200-Hz sampling rate, 1-Hz sensors) arranged in 9 small arrays deployed in 2 concentric rings about 1 km (A-arrays) and 5 km (B-array) radius. The coherence values of the scattering points have been computed in the crustal volume, for a finite time-window along all array stations given the hypothesized origin time and source location. The resulting images can be seen as a (relative) joint log-likelihood of any point in the subsurface that have contributed to the full set of observed seismograms.

Sea-Level climate variability in the Mediterranean Sea

Sea-level variability is characterized by multiple interacting factors described in the Fourth Assessment Report (Bindoff et al., 2007) of the Intergovernmental Panel on Climate Change (IPCC) that act over wide spectra of temporal and spatial scales. In Church et al. (2010) sea-level variability and changes are defined as manifestations of climate variability and change. The European Environmental Agency (EEA) defines sea level as one of most important indicators for monitoring climate change, as it integrates the response of different components of the Earths system and is also affected by anthropogenic contributions (EEA, 2011). The balance between the different sea-level contributions represents an important source of uncertainty, involving stochastic processes that are very difficult to describe and understand in detail, to the point that they are defined as an enigma in Munk (2002). Sea-level rate estimates are affected by all these uncertainties, in particular if we look at possible responses to sea-level contributions to future climate. At the regional scale, lateral fluxes also contribute to sea-level variability, adding complexity to sea-level dynamics. The research strategy adopted in this work to approach such an interesting and challenging topic has been to develop an objective methodology to study sea-level variability at different temporal and spatial scales, applicable in each part of the Mediterranean basin in particular, and in the global ocean in general, using all the best calibrated sources of data (for the Mediterranean): in-situ, remote-sensig and numerical models data. The global objective of this work was to achieve a deep understanding of all of the components of the sea-level signal contributing to sea-level variability, tendency and trend and to quantify them.

Oltre Oceano: Altri orizzonti del possibile. Epistemiologie di Abya Yala e progettualità esistenziale “Lekil Kuxlejal,Sumak Kawsay, Kajkrasa Ruyina”. Filosofie della Terra per una Pedagogia del Buen Vivir

Questo elaborato propone alcune riflessioni sulla necessità urgente di un nuovo paradigma educativo, mediante la re-organizzazione delle scienze della conoscenza, scienze in parole di Morin, “disgiunte e frazionate, inadeguate ad affrontare problemi che richiedono oggi approcci multidisciplinari”. La sfida: affrontare i nuovi problemi di una convivenza planetaria, attraverso le connessioni del pensiero ecologico, in questo studio asse centrale delle cosmovisioni e della Sapienza ancestrale dei Popoli di AbyaYala (America Latina). Popoli in cui la Vita come orizzonte di Armonia ed Equilibrio si concretizza in pratiche di Vita Quotidiana grazie ad una Pedagogia del BuenVivir, inclusiva e partecipativa, rispettosa della diversità biologica e delle differenze culturali, nonché della Sacralità della Terra e della Vita in tutte le sue manifestazioni. La cornice teorica considerata fa riferimento in modo particolare a: L’Ecologia della Mente (Bateson); Il problematicismo Pedagogico e l’Educazione alla Progettualità Esistenziale (G.M.Bertin, Contini); l’Ecologia dei Saperi e le Epistemologie del Sud (Boaventura di Sousa Santos, sociologo portoghese), in modo da tessere ponti di dialogo fra le diverse discipline, in particolare fra la pedagogia, la geografia, l’antropologia, la filosofia, la sociologia, la letteratura, il diritto e anche con le neuroscienze.

Chemical, biochemical and microbiological indicators to assess soil quality in temperate agro-ecosystems

Soil is a critically important component of the earth’s biosphere. Developing agricultural production systems able to conserve soil quality is essential to guarantee the current and future capacity of soil to provide goods and services. This study investigates the potential of microbial and biochemical parameters to be used as early and sensitive soil quality indicators. Their ability to differentiate plots under contrasting fertilization regimes is evaluated based also on their sensitivity to seasonal fluctuations of environmental conditions and on their relationship with soil chemical parameters. Further, the study addresses some of the critical methodological aspects of microplate-based fluorimetric enzyme assays, in order to optimize assay conditions and evaluate their suitability to be used as a toll to asses soil quality. The study was based on a long-term field experiment established in 1966 in the Po valley (Italy). The soil was cropped with maize (Z. mays L.) and winter wheat (T. aestivum L.) and received no organic fertilization, crop residue or manure, in combination with increasing levels of mineral N fertilizer. The soil microbiota responded to manure amendment increasing it biomass and activity and changing its community composition. Crop residue effect was much more limited. Mineral N fertilization stimulated crop residue mineralization, shifted microbial community composition and influenced N and P cycling enzyme activities. Seasonal fluctuations of environmental factors affected the soil microbiota. However microbial and biochemical parameters seasonality did not hamper the identification of fertilization-induced effects. Soil microbial community abundance, function and composition appeared to be strongly related to soil organic matter content and composition, confirming the close link existing between these soil quality indicators. Microplate-based fluorimetric enzyme assays showed potential to be used as fast and throughput toll to asses soil quality, but required proper optimization of the assay conditions for a precise estimation of enzymes maximum potential activity.

New water use efficiency strategies to cope with climate change

Crop water requirements are important elements for food production, especially in arid and semiarid regions. These regions are experience increasing population growth and less water for agriculture, which amplifies the need for more efficient irrigation. Improved water use efficiency is needed to produce more food while conserving water as a limited natural resource. Evaporation (E) from bare soil and Transpiration (T) from plants is considered a critical part of the global water cycle and, in recent decades, climate change could lead to increased E and T. Because energy is required to break hydrogen bonds and vaporize water, water and energy balances are closely connected. The soil water balance is also linked with water vapour losses to evapotranspiration (ET) that are dependent mainly on energy balance at the Earth’s surface. This work addresses the role of evapotranspiration for water use efficiency by developing a mathematical model that improves the accuracy of crop evapotranspiration calculation; accounting for the effects of weather conditions, e.g., wind speed and humidity, on crop coefficients, which relates crop evapotranspiration to reference evapotranspiration. The ability to partition ET into Evaporation and Transpiration components will help irrigation managers to find ways to improve water use efficiency by decreasing the ratio of evaporation to transpiration. The developed crop coefficient model will improve both irrigation scheduling and water resources planning in response to future climate change, which can improve world food production and water use efficiency in agriculture.