Abstract [en]

Magma plumbing systems represent the physical framework of magma transport and storage from the source region in the mantle, through the crust, until reaching the surface in a volcanic eruption. Characterising the different aspects of magma plumbing, in particular the distribution of magma storage zones throughout the crust, is of key importance to better understand the behaviour of individual volcanoes. In particular, shallow crustal magma storage and associated magma-crust interaction processes could potentially explain some of the worlds most unpredictable and explosive volcanoes. This thesis studies magma plumbing architecture in the Sunda Arc (Indonesia), and the North Atlantic Igneous Province, based on elemental and isotope geochemistry, and derived petrological modelling.

In this study, I have employed petrological models, so called geothermobarometers, to calculate pressures and temperatures (P-T) of crustal magma storage. Geothermobarometers are calibrated thermodynamic formulations based on the composition of magmatic minerals and their co-existing melt as a function of the P-T conditions of crystallisation. Using the calculated P-T estimates, I was able to derive the depth of magma storage, and thereby reconstruct the architecture of magma storage systems. A number of different geothermobarometers based on different mineral phases, including plagioclase, clinopyroxene and olivine, were used for this purpose,

The geothermobarometric modelling was combined with additional elemental and isotope geochemical analyses, as well as collaborations with geophysical investigations. These additional approaches were used to corroborate the findings of the geothermobarometric modelling, and also to model and quantify magma-crust interaction processes that take place during crustal magma storage, such as assimilation of crustal lithologies into the magmatic system.

The findings of this thesis build upon the growing body of evidence in support of the prevalence of shallow magma storage in different volcanic settings worldwide. This realisation is relevant to volcano monitoring and hazard mitigation worldwide.

Abstract [en]

Within the KRAKMON project for multiparameter monitoring of Anak Krakatau volcano (Indonesia), a network of temporary stations was installed on the islands of the Krakatau complex as well as in the surrounding areas of the Sunda Strait, Sumatra and Java. The network was operated from June 2005 until January 2006. More than 700 local events were recorded during this experiment, and travel times from these events were used to perform a tomographic inversion for P and S velocities and for the Vp/Vs ratio. In this study, special attention was paid to the validation of the computed model based on different tests, such as inversion of independent data subsets and synthetic modeling. Although the network configuration and the distribution of the events are not favorable for high-quality tomographic imaging, we have obtained some important and robust features which give information about sources of volcanic activity in the Krakatau complex. The most interesting feature of this study is a zone of high Vp/Vs ratio beneath the Krakatau complex. At depths down to 4 km depth we observe anticorrelation of higher P- and lower S-velocities that leads to Vp/Vs ratio higher than 2. This is a probable indicator of the presence of partially molten and/or with high fluid content material with a composition corresponding to deeper layers. It is important that the anomaly of high Vp/Vs ratio beneath the Krakatau complex appears to be separated in two parts at a depth of 5-6 km. This fits to results of geobarometric analysis that presume the existence of several levels of magma chambers beneath Anak Krakatau.

Keyword

Seismic tomography, Krakatau, Volcanism, Magma chamber

National Category

Earth and Related Environmental Sciences

Research subject

Earth Science with specialization in Mineral Chemistry, Petrology and Tectonics

Abstract [en]

Kelut volcano, East Java, is an active volcanic complex hosting a summit crater lake that has been the source of some of Indonesia’s most destructive lahars. In November 2007, an effusive eruption lasting approximately 7 months led to the formation of a 260-m-high and 400-m-wide lava dome that displaced most of the crater lake. The 2007–2008 Kelut dome comprises crystal-rich basaltic andesite with a texturally complex crystal cargo of strongly zoned and in part resorbed plagioclase (An47–94), orthopyroxene (En64–72, Fs24–32, Wo2–4), clinopyroxene (En40–48, Fs14–19, Wo34–46), Ti-magnetite (Usp16–34) and trace amounts of apatite, as well as ubiquitous glomerocrysts of varying magmatic mineral assemblages. In addition, the notable occurrence of magmatic and crustal xenoliths (meta-basalts, amphibole-bearing cumulates, and skarn-type calc-silicates and meta-volcaniclastic rocks) is a distinct feature of the dome. New petrographical, whole rock major and trace element data, mineral chemistry as well as oxygen isotope data for both whole rocks and minerals indicate a complex regime of magma-mixing, decompression-driven resorption, degassing and crystallisation and crustal assimilation within the Kelut plumbing system prior to extrusion of the dome. Detailed investigation of plagioclase textures alongside crystal size distribution analyses provide evidence for magma mixing as a major pre-eruptive process that blends multiple crystal cargoes together. Distinct magma storage zones are postulated, with a deeper zone at lower crustal levels or near the crust-mantle boundary (>15 km depth), a second zone at mid-crustal levels (~10 km depth) and several magma storage zones distributed throughout the uppermost crust (<10 km depth). Plagioclase-melt and amphibole hygrometry indicate magmatic H2O contents ranging from ~8.1 to 8.6 wt.% in the lower crustal system to ~1.5 to 3.3 wt.% in the mid to upper crust. Pyroxene and plagioclase δ18O values range from 5.4 to 6.7 ‰, and 6.5 to 7.6 ‰, respectively. A single whole rock analysis of the 2007–2008 dome lava gave a δ18O value of 7.6 ‰, whereas meta-basaltic and calc-silicate xenoliths are characterised by δ18O values of 6.2 and 10.3 ‰, respectively. Magmatic δ18O values calculated from individual pyroxene and plagioclase analyses range from 5.7 to 7.0 ‰, and 6.2 to 7.4 ‰, respectively. This range in O-isotopic compositions is explained by crystallisation of pyroxenes in the lower to mid-crust, where crustal contamination is either absent or masked by assimilation of material having similar δ18O values to the ascending melts. This population is mixed with isotopically distinct plagioclase and pyroxenes that crystallised from a more contaminated magma in the upper crustal system. Binary bulk mixing models suggest that shallow-level, recycled volcaniclastic sedimentary rocks together with calc-silicates and/or limestones are the most likely contaminants of the 2007–2008 Kelut magma, with the volcaniclastic sediments being dominant.