Machine-learning methods are evaluated to study the intriguing and debated topic of discrimination among different tectonic environments using geochemical and isotopic data. Volcanic rocks ...characterized by a whole geochemical signature of major elements (SiO
2
, TiO
2
, Al
2
O
3
, Fe
2
O
3T
, CaO, MgO, Na
2
O, K
2
O), selected trace elements (Sr, Ba, Rb, Zr, Nb, La, Ce, Nd, Hf, Sm, Gd, Y, Yb, Lu, Ta, Th) and isotopes (
206
Pb/
204
Pb,
207
Pb/
204
Pb,
208
Pb/
204
Pb,
87
Sr/
86
Sr and
143
Nd/
144
Nd) have been extracted from open-access and comprehensive petrological databases (i.e., PetDB and GEOROC). The obtained dataset has been analyzed using support vector machines, a set of supervised machine-learning methods, which are considered particularly powerful in classification problems. Results from the application of the machine-learning methods show that the combined use of major, trace elements and isotopes allows associating the geochemical composition of rocks to the relative tectonic setting with high classification scores (93 %, on average). The lowest scores are recorded from volcanic rocks deriving from back-arc basins (65 %). All the other tectonic settings display higher classification scores, with oceanic islands reaching values up to 99 %. Results of this study could have a significant impact in other petrological studies potentially opening new perspectives for petrologists and geochemists. Other examples of applications include the development of more robust geothermometers and geobarometers and the recognition of volcanic sources for tephra layers in tephro-chronological studies.
We review pre-eruptive dynamics and evidence of open-system behavior in the volcanic plumbing system beneath Campi Flegrei Caldera, together with estimates of magma residence time, magma ascent, and ...mixing-to-eruption timescales. In detail, we compile pre- and syn-eruptive dynamics reported in the literature for (a) the Campanian Ignimbrite ~ 40 ka, (b) the Neapolitan Yellow Tuff (~ 15 ka), and (c) the recent activity within the Phlegrean area. We first summarize geochemical and textural evidence (e.g., magma mixing, crystal disequilibria, vertical zonings, and isotopic records) of open-system behavior for the pyroclasts erupted in the last 40 ky at Campi Flegrei Caldera. We show that the fingerprint of open-system dynamics is ubiquitous in the deposits associated with the volcanic activity at the Campi Flegrei Caldera in the last 40 ky. Then, we describe the results of geophysical and petrological investigations that allow us to hypothesize the structure of the magma feeding system. We point to a trans-crustal magmatic feeding system characterized by a main storage reservoir hosted at ~ 9 km that feeds and interacts with shallow reservoirs, mainly placed at 2–4 km. Finally, we define a scenario depicting pre-eruptive dynamics of a possible future eruption and provide new constraints on timescales of magma ascent with a physical model based on magma-driven ascending dyke theory. Results show that considerably fast ascent velocities (i.e., of the order of m/s) can be easily achieved for eruptions fed by both shallow (i.e., 3–4 km) and deep (i.e., ~ 9 km) reservoirs. Comparing the results from experimental and numerical methods, it emerges that mixing-to-eruption timescales occurring at shallow reservoirs could be on the order of minutes to hours. Finally, we highlight the volcanological implications of our timescale estimates for magma ascent and mixing to eruption. In particular, explosive eruptions could begin with little physical ‘warning’, of the order of days to months. In this case, the onset of volatile saturation might provide pre-eruptive indicators.
Graphical Abstract
Increasing the spatial resolution of Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) is a challenge in microanalysis of geological samples. Smaller sizes for the laser beam ...will allow for (1) high resolution determination of trace element compositions, (2) accurate estimation of crystal/melt partition coefficients, (3) detailed characterization of diffusion profiles, and (4) analysis of fine volcanic glasses. Here, we report about the figures of merit for LA-ICP Quadrupole MS down to a spatial resolution of 5 μm. This study highlights the possibility to achieve suitable limits of detection, accuracy and precision for geological samples even at spatial resolutions of the order of 5 μm. At a beam size of 15 μm, precision (measured as one sigma) and accuracy (expressed as relative deviation from the reference value) are of the order of 10%. At a smaller beam size of 8 um, precision decreases to 15% for concentration above 1.7 μg g–1. As the beam size is decreased to ∼5 μm, precision declines to about 15% and 20% for concentrations above 10 μg g
–1
using
42
Ca and
29
Si as internal standard, respectively. Accuracy is better or equal to 10% and 13% at beam sizes of 15 and 10 μm, respectively. When the spatial resolution is increased to 8 μm, accuracy remains better than 15% and 20% for
42
Ca and
29
Si as internal standard, respectively. We employed such high-resolution techniques to volcanic glasses in ash particles of the 2010 Eyjafjallajökull eruption. Our results are well consistent with the previously reported data obtained at lower spatial resolution, supporting the reliability of the method.
Characterization, correlation and provenance determination of tephra samples in sedimentary sections (tephrochronological studies) are powerful tools for establishing ages of depositional events, ...volcanic eruptions, and tephra dispersion. Despite the large literature and the advancements in this research field, the univocal attribution of tephra deposits to specific volcanic sources remains too often elusive. In this contribution, we test the application of a machine learning technique named Support Vector Machine to attempt shedding new light upon tephra deposits related to one of the most complex and debated volcanic regions on Earth: the Pliocene-Pleistocene magmatism in Italy. The machine learning algorithm was trained using one of the most comprehensive global petrological databases (GEOROC); 17 chemical elements including major (SiO2, TiO2, Al2O3, Fe2O3T, CaO, MgO, MnO, Na2O, K2O, P2O5) and selected trace (Sr, Ba, Rb, Zr, Nb, La, Ce) elements were chosen as input parameters. We first show the ability of support vector machines in discriminating among different Pliocene-Pleistocene volcanic provinces in Italy and then apply the same methodology to determine the volcanic source of tephra samples occurring in the Caio outcrop, an Early Pleistocene sedimentary section located in Central Italy. Our results show that: 1) support vector machines can successfully resolve high-dimensional tephrochronological problems overcoming the intrinsic limitation of two- and three-dimensional discrimination diagrams; 2) support vector machines can discriminate among different volcanic provinces in complex magmatic regions; 3) in the specific case study, support vector machines indicate that the most probable source for the investigated tephra samples is the so-called Roman Magmatic Province. These results have strong geochronological and geodynamical implications suggesting new age constraints (1.4 Ma instead of 0.8 Ma) for the starting of the volcanic activity in the Roman Magmatic Province.
•Machine learning techniques have been applied for tephrochronological studies.•Support vector machines discriminate among different volcanic provinces.•Support vector machines show a Roman geochemical affinity for the studied tephra.•Results antedate the starting activity of the Roman Magmatic Province to ca. 1.4 Ma.
The mixing of magmas is a common phenomenon in explosive eruptions. Concentration variance is a useful metric of this process and its decay (CVD) with time is an inevitable consequence during the ...progress of magma mixing. In order to calibrate this petrological/volcanological clock we have performed a time-series of high temperature experiments of magma mixing. The results of these experiments demonstrate that compositional variance decays exponentially with time. With this calibration the CVD rate (CVD-R) becomes a new geochronometer for the time lapse from initiation of mixing to eruption. The resultant novel technique is fully independent of the typically unknown advective history of mixing - a notorious uncertainty which plagues the application of many diffusional analyses of magmatic history. Using the calibrated CVD-R technique we have obtained mingling-to-eruption times for three explosive volcanic eruptions from Campi Flegrei (Italy) in the range of tens of minutes. These in turn imply ascent velocities of 5-8 meters per second. We anticipate the routine application of the CVD-R geochronometer to the eruptive products of active volcanoes in future in order to constrain typical "mixing to eruption" time lapses such that monitoring activities can be targeted at relevant timescales and signals during volcanic unrest.
Plinian-type eruptions are extremely hazardous, producing pyroclastic fallout and flows extending many kilometres from the vent. The most commonly invoked eruption trigger for Plinian-type eruptions ...is the intrusion of fresh magma, generally associated with precursory ground deformation and seismicity days/weeks before eruption. Closed-system internal triggering has also been proposed, such as protracted crystallisation of magma, which can produce a build-up of exsolved volatiles and thus pressurise the system prior to eruption. On 22–23 April 2015 Calbuco volcano, Chile, produced a sub-Plinian eruption with <3 h seismic precursory activity and no clear deformation signals in the preceding months. Here, we show that petrological and geochemical evidence do not support a hypothesis of eruption triggering due to pre-eruptive intrusion of fresh magma, but instead are consistent with an internal trigger. We found that basaltic andesitic magma was stored at depths between 8 and 12 km (i.e. 230–320 MPa) beneath Calbuco volcano before the 2015 eruption. The stored magma had an initial temperature of 900–950 °C, was water-saturated (5.5–6.5 wt% H2O) and formed phenocrysts of titanomagnetite, orthopyroxene, clinopyroxene and plagioclase cores (An78–93). Gradual cooling of the magma chamber produced thermal gradients and magma convection, evidenced by plagioclase overgrowth rims (An58–77) and blocky microlites (25–250 μm). Our interpretation is that this continuing crystallisation induced second boiling and an over-pressurisation of the system, leading to the rapid onset of the 2015 eruption. Petrological and geochemical evidence therefore shows that a closed-system magma chamber can evolve into a highly explosive eruption with very little precursory warning, posing a challenge for current volcano monitoring paradigms. We propose that internal triggering should be carefully considered as a mechanism for unexpected sub-Plinian eruptions, prompting a potential revision of existing hazard management strategies.
•Pre-eruptive conditions of the 2015 Calbuco eruption have been constrained.•Sub-Plinian eruptions can be triggered by over-pressure due to crystallisation.•Sub-Plinian eruptions can be internally triggered in a few hours.
Magma mixing is common in the Earth. Understanding the dynamics of the mixing process is necessary for dealing with the likely consequences of mixing events in the petrogenesis of igneous rocks and ...the physics of volcanic eruptive triggers. Here, a new apparatus has been developed in order to perform chaotic mixing experiments in systems of melts with high viscosity contrast. The apparatus consists of an outer and an inner cylinder, which can be independently rotated at finite strains to generate chaotic streamlines. The two cylinder axes are offset. Experiments have been performed for ca. 2 h, at 1,400°C under laminar fluid dynamic conditions (
Re
~ 10
−7
). Two end-member silicate melt compositions were synthesized: (1) a peralkaline haplogranite and (2) a haplobasalt. The viscosity ratio between these two melts was of the order of 10
3
. Optical analysis of post-experimental samples reveals a complex pattern of mingled filaments forming a scale-invariant (i.e. fractal) distribution down to the μm-scale, as commonly observed in natural samples. This is due to the development in space and time of stretching and folding of the two melts. Chemical analysis shows strong non-linear correlations in inter-elemental plots. The original end-member compositions have nearly entirely disappeared from the filaments. The generation of thin layers of widely compositionally contrasting interfaces strongly enhances chemical diffusion producing a remarkable modulation of compositional fields over a short-length scale. Notably, diffusive fractionation generates highly heterogeneous pockets of melt, in which depletion or enrichment of chemical elements occur, depending on their potential to spread via chemical diffusion within the magma mixing system. Results presented in this work offer new insights into the complexity of processes expected to be operating during magma mixing and may have important petrological implications. In particular: (1) it is shown that, in contrast with current thinking, rheologically contrasting magmas can mix (i.e. with large proportions of felsic magmas and high viscosity ratios), thus extending significantly the spectrum of geological conditions under which magma mixing processes can occur efficiently; (2) the mixing process cannot be modeled using the classical linear two-end-member mixing model; and (3) the chemical compositions on short-length scales represent snapshots within the process of mixing and therefore may not reflect the final composition of the magmatic system. This study implies that microanalysis on short-length scales may provide misleading information on the parental composition of magmas.
The Euganean Hills (NE Italy) magmatic district represents the final volcanic activity of the Veneto Volcanic Province. Alkaline to subalkaline magmatic suite dominated by intermediate to felsic ...volcanic rocks characterises the latest volcanic activity of the Euganean Hills. Magmatic (intrusive and volcanic) enclaves are common in Euganean Hills trachytes. We used the ability of clinopyroxene to record variations of P, T, and fO2 to reconstruct the geological history of the volcanic enclaves and trachytic host. Despite similar major and trace elements composition, clinopyroxene from host is higher in Ca and Na (and Fe3+) and lower in Mg than enclaves and is slightly enriched in trace elements but with the same pattern distribution. Minor differences in geochemistry and crystal structure of clinopyroxene from enclaves and trachytic host suggest similar parental magmas that differs by small degrees of fractional crystallisation. Clinopyroxene geobarometry performed combining X-ray diffraction with mineral geochemistry for volcanic enclaves–trachytic host combined with amphibole geobarometry for intrusive enclaves and crystal mushes points to a crystallisation pressure range between 4.8–2.0 kbars. Our data support the model of a complex system of magma chambers at intermediate to shallow crustal level where mafic magma accumulated, evolved by fractionation processes and mixed.
Experiments have been performed to determine the effect of deformation on degassing of bubble-bearing melts. Cylindrical specimens of phonolitic composition, initial water content of 1.5 wt.% and ...2 vol.% bubbles, have been deformed in simple-shear (torsional configuration) in an internally heated Paterson-type pressure vessel at temperatures of 798–848 K, 100–180 MPa confining pressure and different final strains. Micro-structural analyses of the samples before and after deformation have been performed in two and three dimensions using optical microscopy, a nanotomography machine and synchrotron tomography. The water content of the glasses before and after deformation has been measured using Fourier Transform Infrared Spectroscopy (FTIR). In samples strained up to a total of γ ∼ 2 the bubbles record accurately the total strain, whereas at higher strains (γ ∼ 10) the bubbles become very flattened and elongate in the direction of shear. The residual water content of the glasses remains constant up to a strain of γ ∼ 2 and then decreases to about 0.2 wt.% at γ ∼ 10. Results show that strain enhances bubble coalescence and degassing even at low bubble volume-fractions. Noticeably, deformation produced a strongly water under-saturated melt. This suggests that degassing may occur at great depths in the volcanic conduit and may force the magma to become super-cooled early during ascent to the Earth’s surface potentially contributing to the genesis of obsidian.