Modeling hydrothermal system: deriving observables and hydrothermal instability in volcanic and non-volcanic setting

Hydrothermal fluids are a fundamental resource for understanding and monitoring volcanic and non-volcanic systems. This thesis is focused on the study of hydrothermal system through numerical modeling with the geothermal simulator TOUGH2. Several simulations are presented, and geophysical and geochemical observables, arising from fluids circulation, are analyzed in detail throughout the thesis. In a volcanic setting, fluids feeding fumaroles and hot spring may play a key role in the hazard evaluation. The evolution of the fluids circulation is caused by a strong interaction between magmatic and hydrothermal systems. A simultaneous analysis of different geophysical and geochemical observables is a sound approach for interpreting monitored data and to infer a consistent conceptual model. Analyzed observables are ground displacement, gravity changes, electrical conductivity, amount, composition and temperature of the emitted gases at surface, and extent of degassing area. Results highlight the different temporal response of the considered observables, as well as the different radial pattern of variation. However, magnitude, temporal response and radial pattern of these signals depend not only on the evolution of fluid circulation, but a main role is played by the considered rock properties. Numerical simulations highlight differences that arise from the assumption of different permeabilities, for both homogeneous and heterogeneous systems. Rock properties affect hydrothermal fluid circulation, controlling both the range of variation and the temporal evolution of the observable signals. Low temperature fumaroles and low discharge rate may be affected by atmospheric conditions. Detailed parametric simulations were performed, aimed to understand the effects of system properties, such as permeability and gas reservoir overpressure, on diffuse degassing when air temperature and barometric pressure changes are applied to the ground surface. Hydrothermal circulation, however, is not only a characteristic of volcanic system. Hot fluids may be involved in several mankind problems, such as studies on geothermal engineering, nuclear waste propagation in porous medium, and Geological Carbon Sequestration (GCS). The current concept for large-scale GCS is the direct injection of supercritical carbon dioxide into deep geological formations which typically contain brine. Upward displacement of such brine from deep reservoirs driven by pressure increases resulting from carbon dioxide injection may occur through abandoned wells, permeable faults or permeable channels. Brine intrusion into aquifers may degrade groundwater resources. Numerical results show that pressure rise drives dense water up to the conduits, and does not necessarily result in continuous flow. Rather, overpressure leads to new hydrostatic equilibrium if fluids are initially density stratified. If warm and salty fluid does not cool passing through the conduit, an oscillatory solution is then possible. Parameter studies delineate steady-state (static) and oscillatory solutions.

Potential fields study of Marsili basin and Palinuro volcanic complex

In this thesis Marsili back-arc basin and Palinuro Volcanic Complex (Southern Tyrrhenian Sea) have been investigated by using magnetic, bathymetric and gravimetric data. A new velocity model of opening of the Marsili basin has been proposed, highlighting the transition from the horizontal spreading of the back-arc to the vertical accretion of the Marsili seamount. Introducing gravity data, Marsili's internal structure has been modeled and a huge portion of the volcano with low density and vanishing magnetization has been detected. Forward modeling of Palinuro Volcanic Complex showed as Palinuro represents the shallowest evidence of a deep tectonic discontinuity and the possible transition domain between the oceanic crust of Marsili Basin and the continental crust related to the Appenninic chain.

Geochronology, geochemistry and isotopic compositions of the volcanic rocks from oceanic (Hawaii) and continental (Eifel) intra-plate environments

The Eifel volcanism is part of the Central European Volcanic Province (CEVP) and is located in the Rhenish Massif, close to the Rhine and Leine Grabens. The Quaternary Eifel volcanism appears to be related to a mantle plume activity. However, the causes of the Tertiary Hocheifel volcanism remain debated. We present geochronological, geochemical and isotope data to assess the geotectonic settings of the Tertiary Eifel volcanism. Based on 40Ar/39Ar dating, we were able to identify two periods in the Hocheifel activity: from 43.6 to 39.0 Ma and from 37.5 to 35.0 Ma. We also show that the pre-rifting volcanism in the northernmost Upper Rhine Graben (59 to 47 Ma) closely precede the Hocheifel volcanic activity. In addition, the volcanism propagates from south to north within the older phase of the Hocheifel activity. At the time of Hocheifel volcanism, the tectonic activity in the Hocheifel was controlled by stress field conditions identical to those of the Upper Rhine Graben. Therefore, magma generation in the Hocheifel appears to be caused by decompression due to Middle to Late Eocene extension. Our geochemical data indicate that the Hocheifel magmas were produced by partial melting of a garnet peridotite at 75-90 km depth. We also show that crustal contamination is minor although the magmas erupted through a relatively thick continental lithosphere. Sr, Nd and Pb isotopic compositions suggest that the source of the Hocheifel magmas is a mixing between depleted FOZO or HIMU-like material and enriched EM2-like material. The Tertiary Hocheifel and the Quaternary Eifel lavas appear to have a common enriched end-member. However, the other sources are likely to be distinct. In addition, the Hocheifel lavas share a depleted component with the other Tertiary CEVP lavas. Although the Tertiary Hocheifel and the Quaternary Eifel lavas appear to originate from different sources, the potential involvement of a FOZO-like component would indicate the contribution of deep mantle material. Thus, on the basis of the geochemical and isotope data, we cannot rule out the involvement of plume-type material in the Hocheifel magmas. The Ko’olau Scientific Drilling Project (KSDP) has been initiated in order to evaluate the long-term evolution of Ko’olau volcano and obtain information about the Hawaiian mantle plume. High precision Pb triple spike data, as well as Sr and Nd isotope data on KSDP lavas and Honolulu Volcanics (HVS) reveal compositional source variations during Ko’olau growth. Pb isotopic compositions indicate that, at least, three Pb end-members are present in Ko’olau lavas. Changes in the contributions of each component are recorded in the Pb, Sr and Nd isotopes stratigraphy. The radiogenic component is present, at variable proportion, in all three stages of Ko’olau growth. It shows affinities with the least radiogenic “Kea-lo8” lavas present in Mauna Kea. The first unradiogenic component was present in the main-shield stage of Ko’olau growth but its contribution decreased with time. It has EM1 type characteristics and corresponds to the “Ko’olau” component of Hawaiian mantle plume. The second unradiogenic end-member, so far only sampled by Honololu lavas, has isotopic characteristics similar to those of a depleted mantle. However, they are different from those of the recent Pacific lithosphere (EPR MORB) indicating that the HVS are not derived from MORB-related source. We suggest, instead, that the HVS result from melting of a plume material. Thus the evolution of a single Hawaiian volcano records the geochemical and isotopic changes within the Hawaiian plume.

The bursting behavior of gas slugs: laboratory and analytical insights into Strombolian volcanic eruptions.

The aim of this thesis is to study how explosive behavior and geophysical signals in a volcanic conduit are related to the development of overpressure in slug-driven eruptions. A first suite of laboratory experiments of gas slugs ascending in analogue conduits was performed. Slugs ascended into a range of analogue liquids and conduit diameters to allow proper scaling to the natural volcanoes. The geometrical variation of the slug in response to the explored variables was parameterised. Volume of gas slug and rheology of the liquid phase revealed the key parameters in controlling slug overpressure at bursting. Founded on these results, a theoretical model to calculate burst overpressure for slug-driven eruptions was developed. The dimensionless approach adopted allowed to apply the model to predict bursting pressure of slugs at Stromboli. Comparison of predicted values with measured data from Stromboli volcano showed that the model can explain the entire spectrum of observed eruptive styles at Stromboli – from low-energy puffing, through normal Strombolian eruptions, up to paroxysmal explosions – as manifestations of a single underlying physical process. Finally, another suite of laboratory experiments was performed to observe oscillatory pressure and forces variations generated during the expansion and bursting of gas slugs ascending in a conduit. Two end-member boundary conditions were imposed at the base of the pipe, simulating slug ascent in closed base (zero magma flux) and open base (constant flux) conduit. At the top of the pipe, a range of boundary conditions that are relevant at a volcanic vent were imposed, going from open to plugged vent. The results obtained illustrate that a change in boundary conditions in the conduit concur to affect the dynamic of slug expansion and burst: an upward flux at the base of the conduit attenuates the magnitude of the pressure transients, while a rheological stiffening in the top-most region of conduit changes dramatically the magnitude of the observed pressure transients, favoring a sudden, and more energetic pressure release into the overlying atmosphere. Finally, a discussion on the implication of changing boundary on the oscillatory processes generated at the volcanic scale is also given.

Reconstruction of Etna recent flank eruptions for assessing diversion barrier project

The thesis contributed to the volcanic hazard assessment through the reconstruction of some historical flank eruptions of Etna in order to obtain quantitative data (volumes, effusion rates, etc.) for characterizing the recent effusive activity, quantifying the impact on the territory and defining mitigation actions for reducing the volcanic risk as for example containment barriers. The reconstruction was based on a quantitative approach using data extracted from aerial photographs and topographic maps. The approach allows to obtain the temporal evolution of the lava flow field and estimating the Time Average Discharge Rate (TADR) by dividing the volume emplaced over a given time interval for the corresponding duration. The analysis concerned the 2001, 1981 and 1928 Etna eruptions. The choice of these events is linked to their impact on inhabited areas. The results of the analysis showed an extraordinarily high effusion rate for the 1981 and 1928 eruptions (over 600 m^3/s), unusual for Etna eruptions. For the 1981 Etna eruption an eruptive model was proposed to explain the high discharge rate. The obtained TADRs were used as input data for simulations of the propagation of the lava flows for evaluating different scenarios of volcanic hazard and analyse different mitigation actions against lava flow invasion. It was experienced how numerical simulations could be adopted for evaluating the effectiveness of barrier construction and for supporting their optimal design. In particular, the gabions were proposed as an improvement for the construction of barriers with respect to the earthen barriers. The gabion barriers allow to create easily modular structures reducing the handled volumes and the intervention time. For evaluating operational constrain an experimental test was carried out to test the filling of the gabions with volcanic rock and evaluating their deformation during transport and placement.

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

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

Toward a Magnitude Scale for Low Frequency Seismicity The LP energy for Campi Flegrei caldera

The energy released during a seismic crisis in volcanic areas is strictly related to the physical processes in the volcanic structure. In particular Long Period seismicity, that seems to be related to the oscillation of a fluid-filled crack (Chouet , 1996, Chouet, 2003, McNutt, 2005), can precedes or accompanies an eruption. The present doctoral thesis is focused on the study of the LP seismicity recorded in the Campi Flegrei volcano (Campania, Italy) during the October 2006 crisis. Campi Flegrei Caldera is an active caldera; the combination of an active magmatic system and a dense populated area make the Campi Flegrei a critical volcano. The source dynamic of LP seismicity is thought to be very different from the other kind of seismicity ( Tectonic or Volcano Tectonic): it’s characterized by a time sustained source and a low content in frequency. This features implies that the duration–magnitude, that is commonly used for VT events and sometimes for LPs as well, is unadapted for LP magnitude evaluation. The main goal of this doctoral work was to develop a method for the determination of the magnitude for the LP seismicity; it’s based on the comparison of the energy of VT event and LP event, linking the energy to the VT moment magnitude. So the magnitude of the LP event would be the moment magnitude of a VT event with the same energy of the LP. We applied this method to the LP data-set recorded at Campi Flegrei caldera in 2006, to an LP data-set of Colima volcano recorded in 2005 – 2006 and for an event recorded at Etna volcano. Experimenting this method to lots of waveforms recorded at different volcanoes we tested its easy applicability and consequently its usefulness in the routinely and in the quasi-real time work of a volcanological observatory.

Transform Tectonics and Non-Volcanic Oceanic Islands

Oceanic islands can be divided, according to their origin, in volcanic and tectonic. Volcanic islands are due to excess volcanism. Tectonic islands are mainly formed due to vertical tectonic motions of blocks of oceanic lithosphere along transverse ridges flanking transform faults at slow and ultraslow mid-ocean ridges. Vertical tectonic motions are due to a reorganization of the geometry of the transform plate boundary, with the transition from a transcurrent tectonics to a transtensive and/or transpressive tectonics, with the formation of the transverse ridges. Tectonic islands can be located also at the ridge–transform intersection: in this case the uplift is due by the movement of the long-lived detachment faults located along the flanks of the mid-ocean ridges. The "Vema" paleoisland (equatorial Atlantic) is at the summit of the southern transverse ridge of the Vema transform. It is now 450 m bsl and it is capped by a carbonate platform 500 m-thick, dated by 87Sr/86Sr at 10 Ma. Three tectonic paleoislands are on the summit of the transverse ridge flanking the Romanche megatrasform (equatorial Atlantic). They are now about 1,000 m bsl and they are formed by 300 m-thick carbonate platforms dated by 87Sr/86Sr, between 11 and 6 Ma. The tectonic paleoisland “Atlantis Bank" is located in the South-Western Indian Ridge, along the Atlantis II transform, and it is today 700 m bsl. The only modern example of oceanic tectonics island is the St. Paul Rocks (equatorial Atlantic), located along the St. Paul transform. This archipelago is the top of a peridotitic massif that it is now a left overstep undergoing transpression. Oceanic volcanic islands are characterized by rapid growth and subsequent thermal subsidence and drowning; in contrast, oceanic tectonic islands may have one or more stages of emersion related to vertical tectonic events along the large oceanic fracture zones.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

Inversion of low frequency seismo-volcanic events

The inversion of seismo-volcanic events is performed to retrieve the source geometry and to determine volumetric budgets of the source. Such observations have shown to be an important tool for the seismological monitoring of volcanoes. We developed a novel technique for the non-linear constrained inversion of low frequency seismo-volcanic events. Unconstrained linear inversion methods work well when a dense network of broadband seismometers is available. We propose a new constrained inversion technique, which has shown to be efficient also in a reduced network configuration and a low signal-noise ratio. The waveform inversion is performed in the frequency domain, constraining the source mechanism during the event to vary only in its magnitude. The eigenvectors orientation and the eigenvalue ratio are kept constant. This significantly reduces the number of parameters to invert, making the procedure more stable. The method has been tested over a synthetic dataset, reproducing realistic very-long-period (VLP) signals of Stromboli volcano. The information obtained by performing the synthetic tests is used to assess the reliability of the results obtained on a VLP dataset recorded on Stromboli volcano and on a low frequency events recorded at Vesuvius volcano.

Miocene basin development and volcanism along a strike-slip to flat-slab subduction transition: Stratigraphy, geochemistry, and geochronology of the central Wrangell volcanic belt, Yakutat-North America collision zone

New geochronologic, geochemical, sedimentologic, and compositional data from the central Wrangell volcanic belt (WVB) document basin development and volcanism linked to subduction of overthickened oceanic crust to the northern Pacific plate margin. The Frederika Formation and overlying Wrangell Lavas comprise >3 km of sedimentary and volcanic strata exposed in the Wrangell Mountains of south-central Alaska (United States). Measured stratigraphic sections and lithofacies analyses document lithofacies associations that reflect deposition in alluvial-fluvial-lacustrine environments routinely influenced by volcanic eruptions. Expansion of intrabasinal volcanic centers prompted progradation of vent-proximal volcanic aprons across basinal environments. Coal deposits, lacustrine strata, and vertical juxtaposition of basinal to proximal lithofacies indicate active basin subsidence that is attributable to heat flow associated with intrabasinal volcanic centers and extension along intrabasinal normal faults. The orientation of intrabasinal normal faults is consistent with transtensional deformation along the Totschunda-Fairweather fault system. Paleocurrents, compositional provenance, and detrital geochronologic ages link sediment accumulation to erosion of active intrabasinal volcanoes and to a lesser extent Mesozoic igneous sources. Geochemical compositions of interbedded lavas are dominantly calc-alkaline, range from basaltic andesite to rhyolite in composition, and share geochemical characteristics with Pliocene-Quaternary phases of the western WVB linked to subduction-related magmatism. The U/Pb ages of tuffs and Ar-40/Ar-39 ages of lavas indicate that basin development and volcanism commenced by 12.5-11.0 Ma and persisted until at least ca. 5.3 Ma. Eastern sections yield older ages (12.5-9.3 Ma) than western sections (9.6-8.3 Ma). Samples from two western sections yield even younger ages of 5.3 Ma. Integration of new and published stratigraphic, geochronologic, and geochemical data from the entire WVB permits a comprehensive interpretation of basin development and volcanism within a regional tectonic context. We propose a model in which diachronous volcanism and transtensional basin development reflect progressive insertion of a thickened oceanic crustal slab of the Yakutat microplate into the arcuate continental margin of southern Alaska coeval with reported changes in plate motions. Oblique northwestward subduction of a thickened oceanic crustal slab during Oligocene to Middle Miocene time produced transtensional basins and volcanism along the eastern edge of the slab along the Duke River fault in Canada and subduction-related volcanism along the northern edge of the slab near the Yukon-Alaska border. Volcanism and basin development migrated progressively northwestward into eastern Alaska during Middle Miocene through Holocene time, concomitant with a northwestward shift in plate convergence direction and subduction collision of progressively thicker crust against the syntaxial plate margin.