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.

Use of solvents and environmental friendly materials for applications in Green Chemistry

The present Thesis studies three alternative solvent groups as sustainable replacement of traditional organic solvents. Some aspects of fluorinated solvents, supercritical fluids and ionic liquids, have been analysed with a critical approach and their effective “greenness” has been evaluated from the points of view of the synthesis, the properties and the applications. In particular, the attention has been put on the environmental and human health issues, evaluating the eco-toxicity, the toxicity and the persistence, to underline that applicability and sustainability are subjects with equal importance. The “green” features of fluorous solvents and supercritical fluids are almost well-established; in particular supercritical carbon dioxide (scCO2) is probably the “greenest” solvent among the alternative solvent systems developed in the last years, enabling to combine numerous advantages both from the point of view of industrial/technological applications and eco-compatibility. In the Thesis the analysis of these two classes of alternative solvents has been mainly focused on their applicability, rather than the evaluation of their environmental impact. Specifically they have been evaluated as alternative media for non-aqueous biocatalysis. For this purpose, the hydrophobic ion pairing (HIP), which allows solubilising enzymes in apolar solvents by an ion pairing between the protein and a surfactant, has been investigated as effective enzymatic derivatisation technique to improve the catalytic activity under homogeneous conditions in non conventional media. The results showed that the complex enzyme-surfactant was much more active both in fluorous solvents and in supercritical carbon dioxide than the native form of the enzyme. Ionic liquids, especially imidazolium salts, have been proposed some years ago as “fully green” alternative solvents; however this epithet does not take into account several “brown” aspects such as their synthesis from petro-chemical starting materials, their considerable eco-toxicity, toxicity and resistance to biodegradation, and the difficulty of clearly outline applications in which ionic liquids are really more advantageous than traditional solvents. For all of these reasons in this Thesis a critical analysis of ionic liquids has been focused on three main topics: i) alternative synthesis by introducing structural moieties which could reduce the toxicity of the most known liquid salts, and by using starting materials from renewable resources; ii) on the evaluation of their environmental impact through eco-toxicological tests (Daphnia magna and Vibrio fischeri acute toxicity tests, and algal growth inhibition), toxicity tests (MTT test, AChE inhibition and LDH release tests) and fate and rate of aerobic biodegradation in soil and water; iii) and on the demonstration of their effectiveness as reaction media in organo-catalysis and as extractive solvents in the recovery of vegetable oil from terrestrial and aquatic biomass. The results about eco-toxicity tests with Daphnia magna, Vibrio fischeri and algae, and toxicity assay using cultured cell lines, clearly indicate that the difference in toxicity between alkyl and oxygenated cations relies in differences of polarity, according to the general trend of decreasing toxicity by decreasing the lipophilicity. Independently by the biological approach in fact, all the results are in agreement, showing a lower toxicity for compounds with oxygenated lateral chains than for those having purely alkyl lateral chains. These findings indicate that an appropriate choice of cation and anion structures is important not only to design the IL with improved and suitable chemico-physical properties but also to obtain safer and eco-friendly ILs. Moreover there is a clear indication that the composition of the abiotic environment has to be taken into account when the toxicity of ILs in various biological test systems is analysed, because, for example, the data reported in the Thesis indicate a significant influence of salinity variations on algal toxicity. Aerobic biodegradation of four imidazolium ionic liquids, two alkylated and two oxygenated, in soil was evaluated for the first time. Alkyl ionic liquids were shown to be biodegradable over the 6 months test period, and in contrast no significant mineralisation was observed with oxygenated derivatives. A different result was observed in the aerobic biodegradation of alkylated and oxygenated pyridinium ionic liquids in water because all the ionic liquids were almost completely degraded after 10 days, independently by the number of oxygen in the lateral chain of the cation. The synthesis of new ionic liquids by using renewable feedstock as starting materials, has been developed through the synthesis of furan-based ion pairs from furfural. The new ammonium salts were synthesised in very good yields, good purity of the products and wide versatility, combining low melting points with high decomposition temperatures and reduced viscosities. Regarding the possible applications as surfactants and biocides, furan-based salts could be a valuable alternative to benzyltributylammonium salts and benzalkonium chloride that are produced from non-renewable resources. A new procedure for the allylation of ketones and aldehydes with tetraallyltin in ionic liquids was developed. The reaction afforded high yields both in sulfonate-containing ILs and in ILs without sulfonate upon addition of a small amount of sulfonic acid. The checked reaction resulted in peculiar chemoselectivity favouring aliphatic substrates towards aromatic ketones and good stereoselectivity in the allylation of levoglucosenone. Finally ILs-based systems could be easily and successfully recycled, making the described procedure environmentally benign. The potential role of switchable polarity solvents as a green technology for the extraction of vegetable oil from terrestrial and aquatic biomass has been investigated. The extraction efficiency of terrestrial biomass rich in triacylglycerols, as soy bean flakes and sunflower seeds, was comparable to those of traditional organic solvents, being the yield of vegetable oils recovery very similar. Switchable polarity solvents as been also exploited for the first time in the extraction of hydrocarbons from the microalga Botryococcus braunii, demonstrating the efficiency of the process for the extraction of both dried microalgal biomass and directly of the aqueous growth medium. The switchable polarity solvents exhibited better extraction efficiency than conventional solvents, both with dried and liquid samples. This is an important issue considering that the harvest and the dewatering of algal biomass have a large impact on overall costs and energy balance.

A model for diffusive and displacive phase transitions in steels

Heat treatment of steels is a process of fundamental importance in tailoring the properties of a material to the desired application; developing a model able to describe such process would allow to predict the microstructure obtained from the treatment and the consequent mechanical properties of the material. A steel, during a heat treatment, can undergo two different kinds of phase transitions [p.t.]: diffusive (second order p.t.) and displacive (first order p.t.); in this thesis, an attempt to describe both in a thermodynamically consistent framework is made; a phase field, diffuse interface model accounting for the coupling between thermal, chemical and mechanical effects is developed, and a way to overcome the difficulties arising from the treatment of the non-local effects (gradient terms) is proposed. The governing equations are the balance of linear momentum equation, the Cahn-Hilliard equation and the balance of internal energy equation. The model is completed with a suitable description of the free energy, from which constitutive relations are drawn. The equations are then cast in a variational form and different numerical techniques are used to deal with the principal features of the model: time-dependency, non-linearity and presence of high order spatial derivatives. Simulations are performed using DOLFIN, a C++ library for the automated solution of partial differential equations by means of the finite element method; results are shown for different test-cases. The analysis is reduced to a two dimensional setting, which is simpler than a three dimensional one, but still meaningful.

Biomolecular studies in Alzheimer’s Disease models: investigations in vitro and in vivo

The Alzheimer’s disease (AD), the most prevalent form of age-related dementia, is a multifactorial and heterogeneous neurodegenerative disease. The molecular mechanisms underlying the pathogenesis of AD are yet largely unknown. However, the etiopathogenesis of AD likely resides in the interaction between genetic and environmental risk factors. Among the different factors that contribute to the pathogenesis of AD, amyloid-beta peptides and the genetic risk factor apoE4 are prominent on the basis of genetic evidence and experimental data. ApoE4 transgenic mice have deficits in spatial learning and memory associated with inflammation and brain atrophy. Evidences suggest that apoE4 is implicated in amyloid-beta accumulation, imbalance of cellular antioxidant system and in apoptotic phenomena. The mechanisms by which apoE4 interacts with other AD risk factors leading to an increased susceptibility to the dementia are still unknown. The aim of this research was to provide new insights into molecular mechanisms of AD neurodegeneration, investigating the effect of amyloid-beta peptides and apoE4 genotype on the modulation of genes and proteins differently involved in cellular processes related to aging and oxidative balance such as PIN1, SIRT1, PSEN1, BDNF, TRX1 and GRX1. In particular, we used human neuroblastoma cells exposed to amyloid-beta or apoE3 and apoE4 proteins at different time-points, and selected brain regions of human apoE3 and apoE4 targeted replacement mice, as in vitro and in vivo models, respectively. All genes and proteins studied in the present investigation are modulated by amyloid-beta and apoE4 in different ways, suggesting their involvement in the neurodegenerative mechanisms underlying the AD. Finally, these proteins might represent novel potential diagnostic and therapeutic targets in AD.

Experimental Analysis and Numerical Simulation of Quench in Superconducting HTS Tapes and Coils

The quench characteristics of second generation (2 G) YBCO Coated Conductor (CC) tapes are of fundamental importance for the design and safe operation of superconducting cables and magnets based on this material. Their ability to transport high current densities at high temperature, up to 77 K, and at very high fields, over 20 T, together with the increasing knowledge in their manufacturing, which is reducing their cost, are pushing the use of this innovative material in numerous system applications, from high field magnets for research to motors and generators as well as for cables. The aim of this Ph. D. thesis is the experimental analysis and numerical simulations of quench in superconducting HTS tapes and coils. A measurements facility for the characterization of superconducting tapes and coils was designed, assembled and tested. The facility consist of a cryostat, a cryocooler, a vacuum system, resistive and superconducting current leads and signal feedthrough. Moreover, the data acquisition system and the software for critical current and quench measurements were developed. A 2D model was developed using the finite element code COMSOL Multiphysics R . The problem of modeling the high aspect ratio of the tape is tackled by multiplying the tape thickness by a constant factor, compensating the heat and electrical balance equations by introducing a material anisotropy. The model was then validated both with the results of a 1D quench model based on a non-linear electric circuit coupled to a thermal model of the tape, to literature measurements and to critical current and quench measurements made in the cryogenic facility. Finally the model was extended to the study of coils and windings with the definition of the tape and stack homogenized properties. The procedure allows the definition of a multi-scale hierarchical model, able to simulate the windings with different degrees of detail.

Synthesis and electronic properties of luminescent silicon nanocrystals and copper indium sulphide quantum dots

In the last decades, nanomaterials, and in particular semiconducting nanoparticles (or quantum dots), have gained increasing attention due to their controllable optical properties and potential applications. Silicon nanoparticles (also called silicon nanocrystals, SiNCs) have been extensively studied in the last years, due to their physical and chemical properties which render them a valid alternative to conventional quantum dots. During my PhD studies I have planned new synthetical routes to obtain SiNCs functionalised with molecules which could ameliorate the properties of the nanoparticle. However, this was certainly challenging, because SiNCs are very susceptible to many reagents and conditions that are often used in organic synthesis. They can be irreversibly quenched in the presence of alkalis, they can be damaged in the presence of oxidants, they can modify their optical properties in the presence of many nitrogen-containing compounds, metal complexes or simple organic molecules. If their surface is not well-passivated, the oxygen can introduce defect states, or they can aggregate and precipitate in several solvents. Therefore, I was able to functionalise SiNCs with different ligands: chromophores, amines, carboxylic acids, poly(ethylene)glycol, even ameliorating functionalisation strategies that already existed. This thesis will collect the experimental procedures used to synthesize silicon nanocrystals, the strategies adopted to functionalise effectively the nanoparticle with different types of organic molecules, and the characterisation of their surface, physical properties and luminescence (mostly photogenerated, but also electrochemigenerated). I also spent a period of 7 months in Leeds (UK), where I managed to learn how to synthesize other cadmium-free quantum dots made of copper, indium and sulphur (CIS QDs). During my last year of PhD, I focused on their functionalisation by ligand exchange techniques, yielding the first example of light-harvesting antenna based on those quantum dots. Part of this thesis is dedicated to them.