# Transitions between explosive and effusive phases during the cataclysmic 2010 eruption of Merapi volcano, Java, Indonesia

Katie Preece, Ralf Gertisser, Jenni Barclay, Sylvain J. Charbonnier, Jean-Christophe Komorowski, Richard A. Herd

► To cite this version:

Katie Preece, Ralf Gertisser, Jenni Barclay, Sylvain J. Charbonnier, Jean-Christophe Komorowski, et al. Transitions between explosive and effusive phases during the cataclysmic 2010 eruption of Merapi volcano, Java, Indonesia. Bulletin of Volcanology, 2016, 78, 16 pp. 10.1007/s00445-016-1046-z. insu-03581313

HAL Id: insu-03581313

https://insu.hal.science/insu-03581313

Submitted on 19 Feb 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Transitions between explosive and effusive phases during the cataclysmic 2010 eruption of Merapi volcano, Java, Indonesia

Katie Preece¹ Ⓒ · Ralf Gertisser² · Jenni Barclay¹ · Sylvain J. Charbonnier³ · Jean-Christophe Komorowski⁴ · Richard A. Herd¹

Received: 11 November 2015 / Accepted: 27 June 2016 / Published online: 18 July 2016
© The Author(s) 2016. This article is published with open access at Springerlink.com

**Abstract** Transitions between explosive and effusive activity are commonly observed during dome-forming eruptions and may be linked to factors such as magma influx, ascent rate and degassing. However, the interplay between these factors is complex and the resulting eruptive behaviour often unpredict- able. This paper focuses on the driving forces behind the ex- plosive and effusive activity during the well-documented 2010 eruption of Merapi, the volcano's largest eruption since 1872. Time-controlled samples were collected from the 2010 deposits, linked to eruption stage and style of activity. These include scoria and pumice from the initial explosions, dense and scoriaceous dome samples formed via effusive activity, as well as scoria and pumice samples deposited during subplinian column collapse. Quantitative textural analysis of groundmass feldspar microlites, including measurements of areal number density, mean microlite size, crystal aspect ratio, groundmass crystallinity and crystal size distribution analysis,

reveal that shallow pre- and syn-eruptive magmatic processes acted to govern the changing behaviour during the eruption. High-An (up to ~80 mol% An) microlites from early erupted samples reveal that the eruption was likely preceded by an influx of hotter or more mafic magma. Transitions between explosive and effusive activity in 2010 were driven primarily by the dynamics of magma ascent in the conduit, with degassing and crystallisation acting via feedback mechanisms, resulting in cycles of effusive and explosive activity. Explosivity during the 2010 eruption was enhanced by the presence of a 'plug' of cooled magma within the shallow magma plumbing system, which acted to hinder degassing, leading to overpressure prior to initial explosive activity.

**Keywords** Microlites · CSD · Merapi · Explosive eruption · Effusive eruption

Editorial responsibility: K.V. Cashman

**Electronic supplementary material** The online version of this article (doi:10.1007/s00445-016-1046-z) contains supplementary material, which is available to authorized users.

Katie Preece
K.Preece@uea.ac.uk

¹ School of Environmental Sciences, University of East Anglia, Norwich NR4 7TJ, UK

² School of Physical and Geographical Sciences, Keele University, Keele, Staffordshire ST5 5BG, UK

³ School of Geosciences, University of South Florida, Tampa, FL 33620-5201, USA

⁴ Equipe Systèmes Volcaniques, Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR CNRS 7154, 1 rue Jussieu, Paris 75238 Cedex 05, France

## Introduction

Transitions between explosive and effusive activity are com- monly observed at many subduction zone volcanoes. Worldwide, it has been estimated that 95 % of dome eruptions are associated with an explosive component (Newhall and Melson 1983; Ogburn et al. 2015). Dome-forming eruptions have transitioned between effusive behaviour and explosive Vulcanian and/or (sub)Plinian activity at Mount St. Helens, Montserrat, Pinatubo and Guagua Pichincha to name just a few (e.g. Hammer et al. 1999; Cashman and McConnell 2005; Clarke et al. 2007; Wright et al. 2007; Komorowski et al. 2010). Understanding the driving forces behind eruptive style and factors that influence changes in activity is therefore crucial for hazard assessment and monitoring efforts. However, a complex interplay of factors acts in each system to determine eruptive style, with previous results often
reaching seemingly paradoxical conclusions regarding how these factors influence eruption dynamics. For example, recharge with hotter magma has previously been linked to remobilisation of magma, triggering and intensifying the onset of eruption and affecting ascent rate at volcanoes such as Soufrière Hills, Montserrat (Murphy et al. 2000) and El Reventador (Ridolfi et al. 2008). In contrast, at Quizapu, Chile (Ruprecht and Bachmann 2010) and Mount Hood, USA (Koleszar et al. 2012), recharge with hotter magma has been linked to a decrease in explosivity, via the reduction of bulk viscosity, which facilitates degassing and inhibits fragmentation. Magmatic influx to a reservoir may in turn affect magma ascent dynamics, influencing both the rate of magma ascent and whether ascent is sustained or pulsatory in its nature (Wolf and Eichelberger 1997; Scandone et al. 2007). Magma ascent rate governs, and is governed by, volatile exsolution during decompression, with fast ascent often linked to closed-system degassing, resulting in increased potential for explosive activity, and slower ascent often linked to open-system degassing and effusive activity (e.g. Jaupart and Allègre 1991; Villemant et al. 2008). At shallow depths, degassing and crystallisation of groundmass microlites influence magma rheology and may also contribute to regulating the style of eruption (e.g. Melnik and Sparks 1999, 2005; Barmin et al. 2002; Degruyter et al. 2012). Extensive, shallow-level, degassing-induced crystallisation can effectively seal the conduit, resulting in pressure build-up and an explosive eruption, if pressure exceeds a critical threshold related to magma strength (e.g. Stix et al. 1993; Sparks 1997; Melnik and Sparks 1999, 2005; Clarke et al. 2007; Wright et al. 2007; Komorowski et al. 2010; Burgisser et al. 2011). Therefore, a complex interplay between deep and shallow processes related to magma storage, influx, ascent and the rheological changes caused by crystallisation and degassing is known to influence the final eruptive behaviour, often via complex feedback mechanisms.

This paper investigates what caused transitions between explosive and effusive activity during the VEI 4 eruption of Merapi volcano (Indonesia) in 2010, via textural and compositional analysis of feldspar groundmass microlites. Previous studies have used textural and compositional analyses of experimental and natural samples to shed light on the processes and consequences of devolatilisation and crystallisation of groundmass microlites during magma ascent (e.g. Cashman 1988, 1992; Hammer et al. 1999, 2000; Hammer and Rutherford 2002). However, most studies of experimental and natural samples have focussed on a particular eruptive behaviour of a volcano, with comparatively less attention paid to the transitions between eruptive styles (Cashman and McConnell 2005; Castro and Gardener 2008; Martel 2012). The physical changes experienced by magma during ascent have been interpreted via analysis of the textures and compositions of microlites that have grown in response to those

changes. Differences in eruptive style have been attributed to magma supply rate (Cashman and McConnell 2005), as well as exsolved volatile pressure in the conduit (Castro and Gardner 2008). At Merapi, investigation of lava samples from the 1986–1995 domes, indicated a correlation between effusion rate and microlite number density, with similar microlite textures throughout that period, suggesting crystallisation conditions are cyclic (Hammer et al. 2000). A detailed study of microlites in lava dome samples extruded throughout the 2006 Merapi eruption revealed that factors such as magma ascent path, source depth and near-surface degassing conditions are recorded in the final groundmass texture of Merapi rocks (Preece et al. 2013). New insights into the driving forces behind transitions in eruptive style can be gained by careful textural analysis of chronologically controlled natural samples linked to a well-documented eruption. In this paper, a well-constrained sample suite representing different stages of the 2010 eruption is linked to eruption chronology, allowing the petrological and textural features to be interpreted in light of eruptive style and monitoring data. In particular, shallow magmatic processes that occur in the edifice and conduit are investigated via quantitative textural and compositional analysis of feldspar microlites. We hypothesise that initiation and cessation of the 2010 eruption were predominantly driven by magma flux at depth, with transitions between explosive and effusive activity in 2010 primarily controlled by the dynamics of magma ascent in the shallow conduit. Degassing and degassing-induced crystallisation during magma ascent in the conduit played a crucial role in influencing eruptive behaviour, via complex feedback mechanisms resulting in cycles of explosive and effusive activity. The 2010 eruption was the largest at Merapi since 1872, and hence, it is the first time that a larger than ‘normal’ eruption has been well-monitored at Merapi, with seismic, ground deformation and gas emission data available (Surono et al. 2012; Budi-Santoso et al. 2013). This eruption therefore provides an ideal case study in which to examine the driving forces behind sudden transitions in eruptive styles during a single, relatively short-lived eruption.

## 2010 eruption chronology, deposits and sampling

### Eruption chronology

During the latter part of the twentieth century and early twenty-first century, eruptions at Merapi have mainly consisted of effusive lava dome growth and subsequent gravitational collapse to generate small volume pyroclastic density currents (PDCs) (VEI 1–3). In stark contrast, the 2010 eruption (VEI 4) began with explosions rather than the effusion of a lava dome and displayed remarkably high extrusion rates (up to $35 \text{ m}^3 \text{ s}^{-1}$) when dome growth took place subsequently (Surono et al. 2012; Pallister et al. 2013; Ratdomopurbo et al. 2013). Activity in 2010 included
several explosive stages generating a series of laterally directed dome explosions ('blasts'), a subplinian eruption column and maximum PDC runout distances of ~16 km, more than twice those generated in the effusive dome-forming eruption in 2006 (Charbonnier and Gertisser 2008; Charbonnier et al. 2013; Komorowski et al. 2013).

The 2010 eruption chronology (Fig. 1) has previously been documented in detail (e.g. Surono et al. 2012; Pallister et al. 2013; Charbonnier et al. 2013; Komorowski et al. 2013). Komorowski et al. (2013) recognised eight stages of the 2010 eruption, which will be referred to throughout this paper: stage 1 (31 October 2009–26 October 2010), unrest and magmatic intrusion; stage 2 (26 October 2010), initial explosions; stage 3 (29 October–4 November 2010), recurrent rapid dome growth (up to $25 \text{ m}^3 \text{ s}^{-1}$) and destruction; stage 4 (5 November 2010), paroxysmal dome explosions and collapse; stage 5 (5 November 2010), retrogressive summit collapse; stage 6 (5 November 2010), subplinian convective fountain collapse; stage 7 (5–8 November 2010), rapid dome growth ($35 \text{ m}^3 \text{ s}^{-1}$) with alternating

effusive and explosive activity; and stage 8 (8–23 November 2010), declining ash venting and degassing.

## Eruption deposits and sampling

The 2010 eruption deposits were reported in detail by Charbonnier et al. (2013), Cronin et al. (2013), Komorowski et al. (2013) and Preece (2014). The deposits are briefly described and documented here with the aim of relating petrological and textural analysis to eruption chronology and dynamics.

Stage 1 is defined as a period of unrest and magmatic intrusion prior to eruption. Therefore, no deposits or samples are associated with stage 1.

The first deposits of the 2010 eruption were produced during explosive activity on 26 October. These mark the onset of the eruption and the beginning of stage 2. Stage 2 deposits consist of dilute PDC (surge) deposits (Komorowski et al. 2013), valley-confined block-and-ash flows (BAFs) (Charbonnier et al. 2013) and pumice-rich PDCs (Drignon et al. 2014).

<figure><img src="image_2.png" /><figcaption>**Fig. 1** a Map of Java, showing the locality of Merapi volcano. b Map of localities for samples used in textural analysis and the distribution of the deposits erupted on 26 October (stage 2) and on 5 November (stages 4 and 6) plotted onto a shaded relief DEM (after Komorowski et al. 2013).</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Map coordinate system is WGS 84 with decimal degrees. c Eruptive timeline of the 2010 eruption, with stages based on Komorowski et al. (2013) showing expanded section describing the paroxysmal eruptive activity of 5 November and showing samples used for textural analysis</figcaption></figure>
(Fig. 1). Typically, each surge unit consists of two layers; the lower layer composed of massive, grey-coloured or ‘salt and pepper’ coarse ash to fine lapilli. The upper layer is usually a brown-orange or grey-coloured fine ash, which is often stratified. Componentry analysis of clasts and matrix (Charbonnier et al. 2013) reveals that the dominant clast component (~90%) of the stage 2 deposits is fresh, vesicular scoriaceous material (herein referred to as ‘stage 2 scoria’). Stage 2 deposits also contain conspicuous juvenile white pumice clasts (‘stage 2 pumice’), light-coloured dense fragments as well as non-juvenile, hydrothermally altered and accidental lithics (see Online Resource 1 for representative logs).

Stage 3 consisted of repeated, rapid (up to 25 m³ s⁻¹) dome growth and recurrent destruction by explosions and collapses, emplacing dilute PDC (surge) deposits (Komorowski et al. 2013) and valley-confined BAFs (Charbonnier et al. 2013).

The paroxysmal stage 4 activity consisted of five laterally directed dome explosions over a period of 11 min which destroyed the rapidly emplaced lava dome (Komorowski et al. 2013). Stage 4 deposits, formed during these cataclysmic dome explosions, include valley-confined BAFs and associated overbank deposits which were generated via the breakout of confined flows onto interfluve areas, as documented by Charbonnier et al. (2013), as well as high-energy, dilute PDC (surge) or ‘blast’ deposits (Komorowski et al. 2013) (Fig. 1). The total volume of the 2010 deposits is ~36.3 × 10⁶ m³, with >70% of this volume deposited during stage 4 (Charbonnier et al. 2013). Field observations, componentry analysis (Charbonnier et al. 2013) and clast density measurements (Komorowski et al. 2013) reveal bimodal dome density and vesicularity. Stage 4 deposits are essentially monolithological, in that their componentry is dominated (>80%) by dark grey to black, dense (~2.5 g cm⁻³; Komorowski et al. 2013) fragments of the fast-growing dome, which was destroyed by the stage 4 explosions (herein referred to as ‘stage 4 dense dome’ clasts). As well as these dense dome fragments, dark grey to black scoriaceous dome fragments also occur within stage 4 deposits (herein referred to as ‘stage 4 scoriaceous dome’ clasts). This lower density (as low as 1.9 g cm⁻³; Komorowski et al. 2013) scoriaceous dome material reflects the existence of dome regions that had not entirely degassed and produced vesicular juvenile dome fragments. Light grey, dense, crystalline material (similar to the light grey, dense lithology found in stage 2) is found as abundant angular inclusions ranging in size from millimetres to centimetres within the juvenile dome material (herein referred to as ‘light grey inclusions’). Occasionally, this lithology also forms wavy, diffuse bands through the darker dome material. Large, and in some cases, prismatically jointed blocks (up to several metres in diameter) of this inclusion material have also been found loose within the stage 4 PDC deposits (Online Resource 1).

Retrogressive summit dome collapse (Komorowski et al. 2013) during stage 5 (Fig. 1) produced BAFs consisting of

variable lithologies, including older, non-juvenile dome clasts, dense and scoriaceous 2010 dome clasts, light grey dense clasts as well as variable non-juvenile lithics. The deposits often have a distinctive reddish-pink colour. The eruptive lull at the end of stage 5 allowed for the deposition of a layer consisting of orange-pink-coloured fine ash with abundant accretionary lapilli, forming a distinctive marker horizon across Merapi's southern flanks (Online Resource 1).

Stage 6 activity consisted of ash venting and recurrent subplinian convective column collapse, generating PDCs. Stage 6 deposits related to convective column collapse are rich in scoriaceous or pumiceous clasts (Fig. 1 and Online Resource 1). In particular, a stage 6 PDC rich in grey scoriaceous clasts (herein referred to as ‘stage 6 grey scoria’) was sampled in the Kali Gendol and on the Kali Gendol interfluve areas (Fig. 1). Scattered on the surface of the scoria-rich flow deposits, abundant white pumice clasts occur (herein referred to as ‘stage 6 white pumice’) that are also interpreted to be associated with the stage 6 column collapse events (Komorowski et al. 2013). Clast density of pumice is 1.25–1.48 ± 0.2 g cm⁻³ (Komorowski et al. 2013).

The 2010 deposits were characterised during three field seasons (February 2011, July 2011, August 2013). Stratigraphic logs (Online Resource 1) and componentry analysis (Charbonnier et al. 2013) enabled units and lithology types to be linked to the eruptive stages of Komorowski et al. (2013) (Fig. 1). Based on field observations and componentry, >200 samples were collected from logged sites, related to eruption stage and various style of activity. To assess inter-clast chemical variability, samples were analysed by XRF to determine whole-rock major element compositions and by electron microprobe to determine mineral and glass compositions (see Online Resource 2 for samples analysed and Preece (2014) for full sample dataset). Whole-rock compositions are similar for all analysed 2010 juvenile material (54.7–55.7 wt% SiO₂), with light grey inclusion material being slightly less evolved (52.6–54.4 wt% SiO₂) (Online Resource 3). Juvenile whole-rock compositions are similar to the effusive 2006 eruption (Preece et al. 2013), as well as previous dome-forming eruptions (Gertisser et al. 2012). To assess inter-clast textural variability, ~50 thin sections were studied, with a sub-set of these samples chosen for feldspar microlite textural and compositional analysis, including stage 2 scoria and pumice produced during the initial explosions on 26 October, stage 4 dense and scoriaceous samples from the rapidly grown and destroyed stage 3 lava dome, light grey inclusion material recovered from stage 4 PDC deposits, interpreted to be derived from a pre-eruption conduit plug (Preece 2014; Gertisser et al. 2015) as well as grey scoria and white pumice from stage 6 subplinian column collapse. This analysis therefore provides the means to compare groundmass microlite textures produced via effusive and explosive activity that occurred during the main stages of the 2010 eruption.
## Methods

Mineral and glass major and minor element compositions were determined using Cameca SX-100 electron microprobes at the University of Cambridge, The Open University and the University of Edinburgh. Feldspar microlites were analysed using a beam diameter of 1, 5 or 10 µm, a 15–20-kV accelerating voltage and a 10–20 nA beam current. Groundmass glass was analysed using a 5–10-µm beam size, a 15-kV accelerating voltage and a 4–10-nA beam current. Na was always counted first to minimise alkali migration. Natural silicate minerals were used as primary standards to calibrate instruments and as secondary in-run standards to monitor precision and accuracy during mineral analyses. Rhyolitic glasses, including KE-12, a peralkaline obsidian from Kenya, were used as secondary in-run standards to monitor glass analyses. Detection limits were ~100–500 ppm for the elements analysed and analytical uncertainties, based on repeat analyses of natural mineral and rhyolitic glass standards, were ~1–3% (relative).

Textural analysis of feldspar microlites was carried out using the methods of Preece et al. (2013). In brief, multiple back-scattered electron (BSE) images of continuous areas of groundmass were acquired with a JEOL JSM 5900 LV scanning electron microscope (SEM) at the University of East Anglia, using an accelerating voltage of 20 kV and a working distance of 10 mm at ×1000 or ×1500 magnification. Images were digitally stitched together, and a minimum of ~600 feldspar microlites in each sample were outlined manually using Adobe Illustrator. Measurements of feldspar microlite size and abundance, sample vesicularity and the area occupied by other mineral phases (i.e. pyroxene, Fe-Ti oxide and apatite microlites) were made using ImageJ. Crystal length and width measurements were converted to 3D crystal habits using CSDslice (Morgan and Jerram 2006) to produce a characteristic aspect ratio (short, *S*; intermediate, *I*; long axis, *L*) for each sample. CSDCorrections (version 1.39) (Higgins 2000) was utilised to calculate crystal size distributions (CSDs), using the ellipsoid minor axis for samples with acicular crystals and the ellipsoid major axis for those with rectangular prism aspect ratios, as recommended in CSDCorrections as it reflects the most likely 2D intersections for these shapes. For all samples, the rock fabric was massive and the crystal aspect ratio calculated with CSDslice was used. Crystal roundness was estimated to be 0.1 (on a scale from 0 to 1, where 0 is a euhedral block and 1 is an ellipsoid). The area percent of vesicles was used to correct for sample vesicularity in the calculation of crystal population density. The number of logarithmic size intervals, or bins, was set at five per decade so that each bin is 1.6 times the size of the next smallest bin, and any bins containing less than five crystals were removed from the CSD as they are not precise, as suggested by Higgins (2000).

## Results

### Groundmass feldspar microlite textures and compositions

Visual inspection of the groundmass in texturally representative samples reveals clear differences between each lithology type (Fig. 2). Textural analysis of feldspar microlites, the most abundant groundmass phase, reveals variations in crystal size, population density and morphology between different lithologies (Fig. 3 and Table 1). Microlite abundance, expressed as areal number density ($N_A$), describes the number of crystals per unit area (mm²) on a vesicle- and phenocryst-free basis. The highest number densities are in samples from the stage 4 dome, especially the dense dome samples (up to 63,962 mm⁻²). The lowest number densities are recorded in samples of stage 2 and 6 pumice (minimum 3230 mm⁻²). For most lithologies, mean crystal area ranges between ~4–15 µm², although the average crystal size microlites in the white pumice samples is larger (28 and 49 µm²). Feldspar microlite crystallinity ($\phi$) is the fraction of groundmass area that is occupied by feldspar microlites and, following the calculations of Hammer et al. (2000), is expressed on a vesicle-, phenocryst- and other crystal phase-free area, meaning that the remaining area is composed of only groundmass glass, thereby taking into account the area of liquid available for late-stage crystallisation. It is a function of both microlite abundance and size, as many small crystals may produce the same crystallinity as few large crystals. Feldspar microlite crystallinity ranges from 0.05 to 0.32 in 2010 samples (Table 1). The morphology of the crystals calculated using CSDslice (Morgan and Jerram 2006) is generally revealed to be acicular, apart from those in the stage 2 scoria and the stage 6 white pumice, which are classed as having a rectangular prism habit. Short/long (*S*/*L*) axis ratios range from 0.10 to 0.67, with acicular shaped crystals having a low *S*/*L* and more blocky shaped crystals having a higher *S*/*L*. The stage 6 white pumice microlites have the highest *S*/*L*, with both sampled populations calculated at 0.67 (Fig. 3). In summary, the 2010 products have either many small microlites (e.g. stage 4 dense dome) or few large microlites, such as in the stage 6 white pumice samples. The relatively few and large microlites of the white pumice have blocky morphologies, whereas the abundant small microlites in the dome samples are acicular in shape. Light grey inclusion material has a highly crystalline groundmass, which has not been analysed quantitatively because many microlites are touching, making discrimination of separate crystals difficult (Fig. 2g, h). However, qualitative observations show that the highly crystalline groundmass is composed of tabular or equant looking crystals, with low proportions of groundmass glass. The remaining glass often appears speckled, with heterogeneities of light and dark areas in backscatter images. Within the light grey inclusions, biotite (Fig. 2g) and a crystalline $SiO_2$ phase, cristobalite, are also present (Fig. 2h). Platy cristobalite crystals fill small vesicles and are pervasive within the groundmass,
<figure><img src="image_4.png" /><figcaption>**Fig. 2** a–1 Representative BSE images of 2010 lithologies showing feldspar microlites (*fsp*), pyroxene microlites (*pyx*), Fe-Ti oxide microlites (*ox*), groundmass glass (*g*) and vesicles (*ves*). Note the presence of biotite (*bio*) and cristobalite (*cr*) and the speckled</figcaption></figure>

appearance of the groundmass in the light grey inclusion material. All images were taken at ×1500 magnification except those of stage 2 and stage 6 pumice samples taken at ×1000 magnification and *g* taken at ×400 magnification

often with ‘fish-scale’ cracked morphology or a microbotryoidal texture. Biotite has not previously been observed in Merapi products before 2010. Costa et al. (2013) note the presence of biotite in 2010 samples, with Preece (2014)

proposing its presence is predominantly restricted to the light grey inclusion material.

The two-dimensional intersection measurements as outlined above were converted into 3D data using *CSDSlice*
<figure><img src="image_5.png" /><figcaption>**Fig. 2** (continued)</figcaption></figure>

(Morgan and Jerram 2006) and CSDCorrections (Higgins 2000) to produce CSDs, as described in Preece et al. (2013)

(Fig. 4). All of the CSDs display curved trends, which may be divided into two or three segments (Fig. 4). The decrease in

**Fig. 3** Variation in textural parameters in 2010 products, plotted with those from 2006 dome products (Preece et al. 2013) for comparison. **a** Areal feldspar microlite number density (N<sub>A</sub> mm<sup>-2</sup>) vs. mean crystal area (μm<sup>2</sup>), with solid lines indicating modelled microlite crystallinity at a particular N<sub>A</sub> and crystal size, labelled with percentage values. **b** Areal feldspar microlite number density (N<sub>A</sub> mm<sup>-2</sup>) vs. feldspar microlite crystallinity (φ) including data from other recent (1986–1995) effusive dome-forming eruptions (Hammer et al. 2000). **c** Areal feldspar microlite number density (N<sub>A</sub> mm<sup>-2</sup>) vs. short/long axis (S/L). **d** Mean crystal area (μm<sup>2</sup>) vs. S/L

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 1** Textural variation of 2010 feldspar microlites</caption><thead><tr><th rowspan="2">Stage</th><th rowspan="2">Lithology</th><th rowspan="2">Sample</th><th colspan="3">2D areal measurements</th><th colspan="3">3D volumetric measurements</th></tr><tr><th>n<sup>a</sup></th><th>Total area (μm<sup>2</sup>)</th><th>N<sub>A</sub> (mm<sup>-2</sup>)</th><th>Area % fsp (vesicle corrected)</th><th>φ</th><th>Mean crystal area (μm<sup>2</sup>)</th><th>Groundmass glass (area %)</th><th>Vesicle area %</th><th>S/L</th><th>Calculated aspect ratio (R<sup>2</sup>)</th><th>Nv (mm<sup>-3</sup>)</th><th>Vol. from inter-area (%)</th></tr></thead><tbody><tr><td>2</td><td>S2S</td><td>M11-28a</td><td>716</td><td>265,942</td><td>10,396</td><td>11.0</td><td>0.12</td><td>10.58</td><td>81.38</td><td>74.8</td><td>0.56</td><td>1.00:1.25:1.80 (0.8883)</td><td>3,050,000</td><td>11.3</td></tr><tr><td>2</td><td>S2P</td><td>M13-26</td><td>696</td><td>375,000</td><td>5704</td><td>8.2</td><td>0.08</td><td>14.63</td><td>89.34</td><td>66.14</td><td>0.56</td><td>1.00:1.50:1.80 (0.7626)</td><td>1,740,000</td><td>7.9</td></tr><tr><td>2</td><td>S2P</td><td>M13-44</td><td>604</td><td>406,000</td><td>3230</td><td>4.8</td><td>0.05</td><td>15.05</td><td>89.30</td><td>54.48</td><td>0.36</td><td>1.00:1.40:2.80 (0.7546)</td><td>1,060,000</td><td>4.9</td></tr><tr><td>4</td><td>DD</td><td>M11-15</td><td>678</td><td>10,682</td><td>63,471</td><td>27.8</td><td>0.29</td><td>4.40</td><td>67.73</td><td>N/A</td><td>0.20</td><td>1.00:1.60:5.00 (0.7245)</td><td>17,500,000</td><td>27.9</td></tr><tr><td>4</td><td>DD</td><td>M11-01</td><td>790</td><td>12,351</td><td>63,962</td><td>25.8</td><td>0.28</td><td>4.03</td><td>66.43</td><td>N/A</td><td>0.10</td><td>1.00:1.50:10.00 (0.6928)</td><td>8,880,000</td><td>25.7</td></tr><tr><td>4</td><td>SD</td><td>M11-02</td><td>736</td><td>15,612</td><td>50,150</td><td>23.9</td><td>0.25</td><td>4.76</td><td>70.70</td><td>6.0</td><td>0.10</td><td>1.00:3.20:10.00 (0.5101)</td><td>10,900,000</td><td>23.9</td></tr><tr><td>6</td><td>GS</td><td>M11-75</td><td>717</td><td>27,568</td><td>26,009</td><td>15.2</td><td>0.16</td><td>5.84</td><td>82.04</td><td>N/A</td><td>0.10</td><td>1.00:1.50:10.00 (0.7869)</td><td>2,850,000</td><td>14.4</td></tr><tr><td>6</td><td>GS</td><td>M11-51</td><td>779</td><td>18,082</td><td>53,338</td><td>24.5</td><td>0.25</td><td>4.59</td><td>72.58</td><td>18.2</td><td>0.18</td><td>1.00:1.40:5.50 (0.7512)</td><td>13,600,000</td><td>24.2</td></tr><tr><td>6</td><td>WP</td><td>M11-50</td><td>686</td><td>330,746</td><td>6111</td><td>29.6</td><td>0.32</td><td>48.52</td><td>62.97</td><td>64.1</td><td>0.67</td><td>1.00:1.15:1.50 (0.8533)</td><td>655,000</td><td>28.0</td></tr><tr><td>6</td><td>WP</td><td>M11-55</td><td>645</td><td>297,422</td><td>5653</td><td>16.0</td><td>0.16</td><td>28.35</td><td>83.01</td><td>59.7</td><td>0.67</td><td>1.00:1.10:1.50 (0.8407)</td><td>853,000</td><td>15.2</td></tr></tbody></table>

Lithology types are as follows: S2S stage 2 scoria, S2P stage 2 pumice, DD dome (dense), SD dome (scoriaceous), GS grey scoria, WP white pumice. Area percent feldspar has been vesicle corrected where necessary

<sup>a</sup>Number of crystals analysed

population density at the smallest sizes, resulting in a left-hand 'downturn' is considered to result from conversion errors in the CSD software, resulting from a lower probability of intersecting smaller crystal sizes (Cashman 1988; Brugger and Hammer 2010a). The y-axis intercept (n°), calculated using the steepest part of the CSD, not including the left-hand downturn, is indicative of the final nucleation density that ranges from 18.84 to 23.37 mm⁻⁴, with stage 6 white pumice samples having the lowest intercept values and the dense dome fragments having the highest. This is consistent with the N<sub>A</sub> data calculated, which also shows that the dense dome samples have the highest crystal number density and the white pumice has the lowest.

Feldspar microlite compositions (Fig. 5) in the 2010 products are highly variable, ranging from An<sub>1</sub>Ab<sub>41</sub>Or<sub>58</sub> to An<sub>84</sub>Ab<sub>16</sub>Or<sub>&lt;1</sub> (Online Resource 4). Microlites from the stage 2 scoria range between An<sub>26</sub>Ab<sub>64</sub>Or<sub>10</sub> and An<sub>80</sub>Ab<sub>19</sub>Or<sub>&lt;1</sub> and those from stage 2 pumice between An<sub>33</sub>Ab<sub>59</sub>Or<sub>8</sub> and An<sub>50</sub>Ab<sub>47</sub>Or<sub>3</sub>. Microlites from dome samples range in composition from An<sub>36</sub>Ab<sub>53</sub>Or<sub>11</sub> to An<sub>76</sub>Ab<sub>21</sub>Or<sub>3</sub>, although nearly 60 % of the measured microlites have &gt;An<sub>60</sub> (Fig. 5). Stage 6 grey scoria microlites have a similar overall range in composition between An<sub>31</sub>Ab<sub>59</sub>Or<sub>10</sub> and An<sub>66</sub>Ab<sub>32</sub>Or<sub>2</sub>, but, in contrast to the dome microlites, only ~15 % of those from the grey scoria are &gt;An<sub>60</sub> (Fig. 5). Stage 6 white pumice microlites are generally more albitic (An<sub>12</sub>Ab<sub>61</sub>Or<sub>27</sub> and An<sub>59</sub>Ab<sub>39</sub>Or<sub>2</sub>), with 90 % of those microlites analysed having ≤An<sub>40</sub>. Microlites from the light grey inclusion material show the widest overall range (An<sub>1</sub>Ab<sub>41</sub>Or<sub>58</sub> to An<sub>84</sub>Ab<sub>16</sub>Or<sub>&lt;1</sub>), although more than 70 % of these crystals are alkali feldspars with &gt;Or<sub>20</sub> (Fig. 5). Microlite compositions are shown with isothermal sections of the dry ternary solvus, calculated using SOLVCALC (Wen and Nekvasil 1994) at 0.5 kbar (Fig. 5). The dry solvus position is not sensitive to pressure changes between 0 and 3 km, the region of expected degassing in the conduit, so depth only controls feldspar composition in terms of the extent of H₂O degassing (Hammer et al. 2000). The microlite compositions with respect to the dry solvus do not necessarily imply that microlite crystallisation occurred at different temperatures during cooling between ~600 and 1000 °C but rather that changes in H₂O activity changed the position of the liquidus and, therefore, the composition of the crystallising feldspar, with enrichment in alkalis due to degassing and the compositional evolution of the residual melt phase.

## Discussion

### Textural and compositional evidence for microlite growth and magma ascent conditions

Groundmass microlite textural differences are due to differences in undercooling (ΔT) during crystallisation, which may occur either during ascent or during dome residence. The 2010 dome, which was destroyed on 5 November (stage 4 dense

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 4** Crystal size distribution (CSD) plots of feldspar microlites from the 2010 eruption. **a** All 2010 CSDs. **b–k** Individual CSDs for each sample (**coloured solid line**), with a regression line fitted to the steepest</figcaption></figure>

part of the CSD curve (**black dashed line**), values for the slope of the regressed line, the y-axis intercept or nucleation density (n°) and the R² value. n denotes the number of crystals analysed in each sample

and scoriaceous dome samples), began to extrude on 1 November (Pallister et al. 2013). Therefore, these 2010 dome clasts resided for a maximum of 4 days within the dome. These dome samples have φ and groundmass glass area values similar to those of the stage 6 grey scoria and white pumice, which erupted explosively and were never part of a dome. This suggests that these microlite textures predominantly reflect conditions during magma ascent rather than microlite growth during dome residence, as also observed in 2006 dome samples (Preece et al. 2013).

Stage 2 scoria and pumice have low N<sub>A</sub> and φ and are closer to intersecting the origin (Fig. 3b), suggesting crystallisation was dominated by nucleation rather than by growth, as also observed in other nucleation-dominated systems (e.g. Pinatubo; Hammer et al. 1999). In stage 4 dome samples, high microlite number densities, small crystal sizes

and acicular morphologies also point to nucleation-dominated crystallisation at high ΔT (e.g. Lofgren 1980). The extrusion rate of the dome that produced the stage 4 samples was up to 25 m³ s⁻¹ (Pallister et al. 2013). In comparison, this is nearly 8 times greater than peak extrusion rates at Merapi in 2006 (Ratdomopurbo et al. 2013) and 78 times greater than in 1995 (Hammer et al. 2000). Compared to previous dome eruptions at Merapi, high extrusion rates and high ascent rates of the 2010 dome magma resulted in high degrees of ΔT related to degassing during fast decompression in the conduit, causing nucleation-dominated crystallisation and yielding textures with many, small acicular crystals. High ascent rates are also confirmed by the lack of reaction rims on many of the amphibole phenocrysts and microphenocrysts in the dome samples (Preece 2014). Values of N<sub>A</sub> are higher for a given φ compared to the 2006 and other recent dome-forming eruptions
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Feldspar microlite compositions plotted in the feldspar ternary (An-Ab-Or) diagram display differences in the compositions between the different 2010 lithologies; 2006 microlites (Preece et al. 2013) are</figcaption></figure>

shown for comparison. *Curves* represent isothermal sections of the dry ternary solvus at 600, 800 and 1000 °C, calculated using SOLVCALC (Wen and Nekvasil 1994)

(Fig. 3b), indicating that crystallisation was dominated by high levels of nucleation in comparison to other recent eruptions. Stage 6 white pumice microlites were presumably also formed during high rates of magma ascent in the subplinian stage of the 5 November paroxysmal activity and yet have different groundmass textures to the dome, with fewer, larger equant microlites (Fig. 3). At very high $\Delta T$, often associated with explosive eruptive behaviour, crystallisation is retarded due to an increase in viscosity with H$_2$O loss (e.g. Hess and Dingwell 1996) and due to a lag time between conditions of supersaturation and crystallisation. Nucleation kinetics are more sluggish than those of gas exsolution, attributed to the time necessary for structural reorganisation of the melt (Hammer et al. 1999; Brugger and Hammer 2010b), leading to a nucleation delay of between a few hours to days, depending on $\Delta T$ (Couch et al. 2003a, b; Brugger and Hammer 2010b; Martel 2012). In the stage 6 white pumice samples, nucleation densities are low for a given $\phi$ compared to other 2010 samples (Fig. 3), indicating that crystallisation proceeded by growth of existing crystals rather than nucleation of new ones. Stage 6 white pumice microlites with large size, low number density and equant morphologies probably crystallised at lower $\Delta T$ under growth-dominated conditions, before the rapid ascent to the surface. The final rapid ascent did not result in a nucleation event and formation of many

small crystals, as seen in the dome samples, because the ascent was sufficiently rapid to hinder crystallisation.

Groundmass microlite textures in stage 6 grey scoria samples are intermediate between those observed in the dome clasts and those in the white pumice, although, in many respects, they are more closely similar to the dome samples and fall on the same nucleation-dominated trend (Fig. 3b). The crystal number densities are more variable than for other lithologies, varying from high number densities, similar to those of the dome, to densities of ~26,000 mm$^{-2}$ intermediate between those of the dome and white pumice (Fig. 3). Small mean crystal sizes and low aspect ratios reveal that crystals are small and acicular, as in the dome samples. High number densities of small, acicular crystals indicate crystallisation at conditions of high $\Delta T$, although not high enough to retard nucleation as recorded by stage 6 white pumice samples. Although quantitative textural analysis was not carried out on the light grey inclusions in the 2010 lava dome, high groundmass crystallinity, low proportions of groundmass glass and the presence of cristobalite suggest that the light grey inclusion material has spent a prolonged period of time crystallising at shallow depths (<50 MPa) within the magmatic system (Couch et al. 2003b). The presence of cristobalite is indicative of extensive late-stage vapour-phase crystallisation as previously observed, for example, in the Mount St. Helens
cryptodome in 1980 and subsequent dome rocks (Hoblitt and Harmon 1993; Pallister et al. 2008). The inhomogeneous groundmass glass is similar to the glass observed from Mount St. Helens (Cashman 1988, 1992), Santiaguito (Scott et al. 2012) and Soufrière Hills, Montserrat (Couch et al. 2003b). This texture has been interpreted as a phase separation and devitrification of the glass during shallow storage or slow extrusion (Cashman 1992; Scott 2012). Another groundmass phase observed in the light grey inclusions is biotite. This is extremely rare at Merapi and has not been recorded in previous eruptive products prior to 2010. Biotite may have a magmatic origin but possibly may also originate from hydrothermal alteration of basaltic andesite or may be derived from metasomatic skarn material (e.g. Panigrahi et al. 2008; Afshooni et al. 2013; Balassone et al. 2013). Biotite composition, especially in terms of its FeO, MgO, MnO, TiO₂ and Al₂O₃ contents, is indicative of its mode of formation. For example, the mole fraction of Mg in the octahedral site (X_Mg) has been used to distinguish between igneous biotite and that which formed in hydrothermal alteration zones of porphyry deposits, with the hydrothermal biotite containing higher X_Mg and less FeO (Selby and Nesbitt 2000; Panigrahi et al. 2008). At Vesuvius for example, distinct mica compositions have been identified, depending on whether they are magmatic or associated with hydrothermal or metasomatic material (Balassone et al. 2013). When compared to the Vesuvius mica compositions, the Merapi biotite are most similar to the biotite of magmatic origin, as they are Al-, Ti- and Fe-rich compared to the skarn and hydrothermal types. As the light grey inclusion material has a highly crystalline groundmass, the biotite potentially formed in the latest magmatic stages from the presumably more evolved residual melt.

In order to further elucidate magmatic ascent conditions, melt H₂O contents were calculated with the plagioclase-liquid hygrometer of Waters and Lange (2015), using the microlite and groundmass glass compositions (Online Resource 2) and a temperature of 900 °C as inputs. The temperature is based on minimum estimates from amphibole thermometry (Costa et al. 2013; Preece 2014). The range in calculated melt H₂O contents and the corresponding pressures, calculated using Papale et al. (2006), is as follows: stage 2 scoria, 2.4–3.1 wt% H₂O (413–624 bars); stage 2 pumice, 2.3–2.6 wt% H₂O (388–470 bars); stage 6 scoria, 2.1–2.8 wt% H₂O (334–526 bars); stage 6 pumice, 2.2–3.1 wt% H₂O (378–642 bars); and light grey inclusions, 1.1–2.3 wt% H₂O (150–429 bars). Calculations were not performed on stage 4 samples, as groundmass glass patches were too small to analyse accurately. The similar H₂O content and pressure range for all juvenile 2010 material correspond to microlite crystallisation depths of 1.2–2.4 km, assuming a crustal density of 2800 kg m⁻³ (Costa et al. 2013) and shallower depths of 0.6–1.5 km for the microlites within the light grey inclusions. This suggests that magma decompression and resultant microlite crystallisation

originated from similar depths throughout the eruption. However, these results must be taken with some caution, as degassing-induced crystallisation may not occur under full equilibrium conditions.

High-An (up to An₈₀) compositions of feldspar microlites in stage 4 dome clasts and stage 2 scoria can result from either higher water pressures or higher temperatures compared to the other, lower An microlites, as experiments show that An content increases with increasing temperature and P(H₂O) (Couch et al. 2003a; Martel 2012). However, hygrometry results suggest that stage 2 microlites were forming at similar P(H₂O) compared to all other stages. As such, it is likely that the high-An microlites in the initial explosions and subsequent dome formed via the interaction with an influx of hotter magma, the presence of which is consistent with previous petrological work (Costa et al. 2013; Preece et al. 2014) and monitoring data from the 2010 eruption (Surono et al. 2012). This interpretation is consistent with high-An microlites documented at Mount Pelée (Martel et al. 2006) and Soufrière Hills volcano (Humphreys et al. 2009), which are thought to be inherited from mafic magma. The stage 6 white pumice microlites grew at low ΔT under growth-dominated conditions, suggesting the magma stalled in the conduit for a period of time, at depths of ~1.4–2.4 km. The very fast final ascent and magma fragmentation before the explosive activity only occurred within the last kilometres of ascent, rather than by very fast ascent from greater depths. Therefore, stage 6 may have been caused by rapid decompression and fragmentation driven by unloading, leading to explosive behaviour. Anorthoclase and more K-rich feldspar microlites that are predominant in the light grey inclusion material suggest that crystallisation occurred either at low P(H₂O) and/or lower temperatures. Hygrometry results are consistent with low P(H₂O), suggesting crystallisation occurred at depths of 0.6–1.5 km, shallower than other 2010 products and in accord with this material originating from a plug in the shallow conduit (Preece 2014; Gertisser et al. 2015).

## Crystal size distribution analysis

Crystal size distributions of the samples are concave-upwards curves, reflecting either shape variability of the sample (Castro et al. 2003) or progressive changes in crystal growth and nucleation rates as a function of changing ΔT in the system. Growth of crystals is dominant in the lower conduit, and nucleation begins to dominate in the upper conduit, leading to crystallisation of smaller crystals. Therefore, different microlite sizes reflect magmatic conditions at different depths (Melnik et al. 2011). As crystal growth rate is not constant over time, due to changes in ΔT during ascent, a single growth rate cannot be accurately used to solve crystallisation times. However, average crystallisation rates may be calculated using known time constraints based on observations of the eruption chronology. For example, during 1–26 October
2010, an increase in seismic activity and summit deformation marked an ‘intrusive phase’ believed to reflect the movement of magma to shallower (<5 km below the summit) regions in the volcano magma plumbing system (Budi-Santoso et al. 2013), likely inducing microlite crystallisation. Based on this observation, it is possible to calculate an average crystallisation rate for the first microlites grown in the stage 2 scoria and pumice, erupted on 26 October. Using the equation $Slope = -1/G\tau$, where G is growth rate and τ is crystallisation time (Marsh 1988), average crystallisation times can be calculated using the shallowest segment of the CSD (representing the largest microlites) and a crystallisation period of 26 days (Fig. 4). This reveals a crystal growth rate of $9.5 \times 10^{-9}$ mm s⁻¹ for the large stage 2 scoria microlites and slower rates of $2.2–3.5 \times 10^{-9}$ mm s⁻¹ for stage 2 pumice microlites. These calculated rates are close to the low end of the range of growth rates ($10^{-6}$ to $10^{-8}$ mm s⁻¹) deemed appropriate for syn-eruptive plagioclase crystallisation based on experimental data (Brugger and Hammer 2010b). The fact that the calculated growth rates are at the lower end of this range may be expected for these early microlites that formed deeper within the conduit during growth-dominated crystallisation conditions. Crystals formed later during ascent, at shallower levels within the conduit, may have faster growth rates. The $N_A-\phi$ relationship (Fig. 6) suggests higher eruption intensity and possibly higher decompression rates, leading to faster growth rates for microlites in stage 4 dome and stage 6 grey scoria samples. In comparison, samples from the 1991 Pinatubo eruption have a higher average feldspar microlite growth rate of $1.7 \times 10^{-7}$ mm s⁻¹ (Hammer et al. 1999), consistent with higher eruption intensity compared to Merapi in 2010. This is also in accord with the $N_A-\phi$ plot (Fig. 6), where the Merapi stage 4 dome and stage 6 grey scoria samples plot within the region between the Pinatubo and previous Merapi effusive dome products, in agreement with the Merapi 2010 eruptive style also being of an intermediate intensity between

Pinatubo and previous Merapi effusive dome-forming eruptions.

## Transitions between explosive and effusive eruptive behaviour

The 2010 eruption was the largest at Merapi since 1872 and displayed transitions between both explosive and effusive behaviour. Whole-rock compositions remained constant throughout the duration of the eruption and are similar to previous effusive dome-forming eruptions, indicating bulk magmatic composition cannot be a major factor in the change in behaviour. Instead, the initiation and cessation of the 2010 eruption was predominantly controlled by magma flux from depth, with shallow conduit processes modulating eruptive style, as evidenced by microlites.

The initial explosion at the onset of the 2010 eruption on 26 October was not preceded by dome extrusion, as was the case for the 5 November paroxysmal explosions. On 26 October, it is probable that there was an influx of hotter, potentially more mafic magma (Fig. 7a), as evidenced by high-An microlites in early products, as well as other petrological indications (e.g. Costa et al. 2013; Preece et al. 2014). Gas accumulation due to magmatic influx and rapid ascent (evidenced by low $\phi$ and $N_A$ in initial explosive products) was not balanced by release through permeable conduit walls and fractures to the surface, leading to a build-up of overpressure in a closed system (Fig. 7a, b). The gas-charged, rapidly risen magma erupted as stage 2 scoria and pumice (Komorowski et al. 2013; Preece 2014; Drignon et al. 2014). In addition, the light grey dense material is interpreted to be a plug of cooled rigid magma that resided within the shallow system before being partially reheated and remobilised during eruption (Preece 2014; Gertisser et al. 2015). This material has an extensively crystallised groundmass, with low volumes of residual melt (glass) and abundant cristobalite. Cristobalite is often

<figure><img src="image_10.png" /><figcaption>**Fig. 6** $N_A-\phi$ relationship of 2010 Merapi samples compared to those from the lower-intensity Merapi 2006 dome-forming eruption (Preece et al. 2013) and older Merapi domes (1986–1995) (Hammer et al. 2000), as well as samples from the higher-intensity pre-climatic stages of the Pinatubo 1991 eruption (Hammer et al. 1999)</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Fig. 7** Schematic diagram of the Merapi conduit during the 2010 eruption. **a** Early October (stage 1): rapid magma ascent caused by input of hotter or more mafic magma, with the presence of a 'plug' in the conduit, causing gas overpressure. **b** This resulted in the initial explosions beginning on 26 October (stage 2), which opened the system, allowing for rapid dome extrusion between 26 October and 4 November (stage 3). **c** Fast magma ascent and open-system degassing promoted rapid nucleation and crystallisation with permeable gas loss during rapid dome extrusion (stage 3). **d** Permeable gas loss, vesicle</figcaption></figure>

collapse and formation of a dense dome, effectively sealed the conduit, producing conditions favourable to explosive activity on 5 November (stage 4) followed by **e** retrogressive gravitational collapse of old domes at the summit (stage 5). **f** Closed-system degassing sustained explosive behaviour, generating a convective column, which collapsed to produce the scoria- and pumice-rich PDCs on 5 November (stage 6). **g** A return to open-system degassing was marked by effusive activity and Strombolian jetting of gas and tephra between 6 and 8 November (stage 7). Lithology abbreviations as in Table 1

precipitated within vesicles via vapour transport in the upper-most conduit, and the extent of cristobalite mineralisation may reduce porosity and permeability (Horwell et al. 2013). It is therefore possible that cristobalite crystallisation within the light grey dense inclusion material contributed to sealing gas escape pathways prior to the 2010 eruption, augmenting gas overpressure. The remnant 2006 dome may also have been a factor in sealing gas escape pathways at the onset of eruption in 2010 (Fig. 7a). Seismic features support the fact that explosions during this initial stage of the eruption were driven by gas pressure in the conduit (Jousset et al. 2013). Volcano tectonic (VT), multi-phase (MP) and very long period (VLP) earthquakes, as well as unprecedented rates of summit deformation, were all detected prior to the eruption, attributed to fluid movement and pressure build-up in the conduit (Surono et al. 2012; Jousset et al. 2013; Budi-Santoso et al. 2013).

Once the juvenile magma had disrupted the 'plug' during stage 2 (Fig. 7b), rapid extrusion of a lava dome became possible (Fig. 7c). Rapid magma ascent and extrusion of a

dense dome, evidenced by microlite textures and the presence of amphiboles without breakdown rims in stage 4 products, as well as by monitoring data (Surono et al. 2012; Pallister et al. 2013), are interpreted to be a contributing factor to the paroxysmal behaviour in 2010. Rapid magma ascent and degassing resulted in high degrees of $\Delta T$, causing a rapid nucleation event and nucleation-dominated crystallisation of many microlites at shallow levels within the conduit (leading to high $N_A$ values for stage 4 dome samples) (Fig. 7c). It is probable that degassing and crystallisation caused rheological stiffening of magma and increased viscosity (e.g. Sparks 1997; Melnik and Sparks 1999, 2005). This could have resulted in increased gas overpressure, if the pressure increase due to crystallisation could not be balanced via gas flow and escape to the exterior (Sparks 1997). Rapid nucleation and crystallisation concentrated the remaining gas into smaller amounts of residual melt, so that more gas bubbles formed as a result of supersaturation, promoting higher levels of bubble connectivity (Sparks 1997; Clarke et al. 2007).
High bubble connectivity promotes permeable gas loss, resulting in vesicle collapse and formation of a dense lava dome (stage 4 dense dome samples) (Fig. 7d), which acted to effectively seal the conduit, forming conditions favourable to short-lived, directed dome explosions driven by conduit processes akin to Vulcanian-like explosions (Clarke et al. 2007) (Fig. 7d). Open-system degassing and crystallisation of the dense dome therefore effectively sealed the conduit, which, when combined with high magma flux rates in the conduit, led to an increase in pressurisation before the 5 November paroxysmal explosions (Fig. 7d). Sealing of the conduit is also indicated by the decrease in SO₂ emissions during rapid dome growth at the beginning of November, compared to the explosive stages of the eruption (Surono et al. 2012).

Effective sealing of the conduit by the dome consequently ‘closed’ the system, so that magma ascending below the dome proceeded to degas, with little or no gas escape from the system (Fig. 7d). The light grey ‘plug’ material may also have contributed to this, as it is found as inclusions within the dense dome lava and as prismatically jointed clasts within the stage 4 PDC deposit. Stage 4 explosions and stage 5 retrogressive summit collapse (Fig. 7e), acted to open the system prior to stage 6. Evidence of closed-system degassing is seen in melt inclusions from stage 6 grey scoria and white pumice clasts (Preece et al. 2014), probably linked to fast magma ascent rates. Closed-system degassing and fragmentation of the magma sustained explosive behaviour, generating a pulsating subplinian convective column, which collapsed to produce the scoria- and pumice-rich PDCs during stage 6 of the 2010 eruption (Fig. 7f). The cataclysmic stages of the eruption on 5 November therefore started with repeated short directed explosions (‘blasts’), driven by conduit conditions similar to Vulcanian eruptions (stage 4), before a subplinian convective column was generated by the rising magma in an open conduit (stage 6). The following day, on 6 November, effusive activity returned, with even higher dome extrusion rates than previously (35 m³ s⁻¹), accompanied by Strombolian (i.e. open-vent) jetting of gas and tephra in the morning of 6 November (Surono et al. 2012; Pallister et al. 2013) (Fig. 7g). This stage represents the transition from explosive (Vulcanian/subplinian) activity back to effusive activity and a return to open-system degassing. Presumably, the rapid dome extrusion that occurred on 6 November could potentially have caused further large explosions by sealing the conduit as occurred during the paroxysmal stage. Instead, observed Strombolian-style gas and tephra fountaining suggests that the preceding stage 6 explosive activity ‘re-opened’ the system, creating a pathway for gas release. Dome growth ceased by 8 November and was followed by dome subsidence and small gas emissions (Surono et al. 2012; Pallister et al. 2013). This indicates that waning magmatic flux inhibited further explosive behaviour.

## Conclusions

In summary, the beginning and end of the 2010 eruption were controlled by magma flux from depth, but eruptive style and transitions between explosive and effusive activity were regulated by shallow conduit processes. Degassing and degassing-induced crystallisation were controlled by magma ascent processes in the conduit, which in turn were responsible for changes in eruptive style. A large influx of hotter, possibly mafic magma from depth, triggered gas overpressure and faster magma ascent rates compared to those observed in 2006. The presence of a ‘plug’ in the conduit, coupled with high ascent rates led to closed-system conditions, which pressurised the system and led to explosions on 26 October. Initial explosions temporarily ‘opened’ the system, allowing ascending magma to degas more freely, and extrude as a lava dome, representing the transition back to effusive activity. Open-system degassing led to a rapid nucleation and crystallisation event and extrusion of dense, degassed lava. This effectively resealed the system, and closed-system degassing of the ascending magma began to prevail, leading to a build-up of gas overpressure until the cataclysmic dome explosions and subplinian stage on 5 November. These explosions then enabled open-conduit conditions, facilitating the transition back to effusive activity, allowing for rapid extrusion and emplacement of a dome the following day. By this stage, the eruption was waning and so further explosive activity did not occur due to decreasing magma flux from depth.

Transitions between explosive and effusive activity in 2010 were driven primarily by the dynamics of magma ascent in the conduit, with degassing and crystallisation acting via feedback mechanisms, resulting in cycles of effusive and explosive activity. For example, explosive activity on 26 October and 5 November acted to ‘open’ the system and allow for efficient degassing, causing the transition to effusive activity. Effusive activity at the beginning of November, promoted the transition back to explosive activity, with rapid dome growth effectively sealing the system. In other words, explosive activity enabled subsequent effusive activity, and effusive activity provoked explosive activity.

The 2010 eruption demonstrates the capacity for dome-forming periods at Merapi to switch to explosive behaviour and that future eruptions may begin explosively, with little warning time and without initial dome growth. Petrological evidence suggests that pre-cursors to a large eruption in the future may include a deep pulse of gas-rich magma and a build-up of gas overpressure or rapid dome extrusion. Warning signs that these processes are occurring may be monitored remotely via seismic networks, satellite monitoring, gas flux and ground deformation measurements. The findings of this study are pertinent to other dome-forming volcanoes, where the dynamics of magma influx, magma ascent, degassing and crystallisation may also play an important role in determining eruptive style.
**Acknowledgments** The authors are grateful to Aurélie Germa (University of South Florida), Susanna Jenkins (University of Bristol) and colleagues in Indonesia from the BBPTK, including N. Cholik, for assistance and discussions during several field seasons. Bertrand Lézé (UEA) is thanked for technical advice and assistance with XRF and SEM analysis. Chiara Petrone and Iris Buisman (University of Cambridge), Andy Tindle (The Open University) and Chris Hayward (University of Edinburgh) are thanked for technical assistance with electron microprobe analysis. We are grateful for the thorough and constructive comments from reviewers Julia Hammer and Thomas Giachetti, as well as Editor Kathy Cashman, which helped to improve and clarify this manuscript. This work has been supported by the Natural Environment Research Council (NERC) through Urgency grant NE/029927/1. KP was funded by a NERC studentship, grant number NE/H524506/1.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## References

Afshooni SZ, Mirnejad H, Esmaeily D, Asadi Haroni H (2013) Mineral chemistry of hydrothermal biotite from the Kahang porphyry copper deposit (NE Isfahan), Central Province of Iran. Ore Geol Rev 54: 214-232

Balassone G, Scordari F, Lacalamita M, Schingaro E, Mormone A, Piochi M, Petti C, Mondillo N (2013) Trioctahedral micas in xenolithic ejecta from recent volcanism of the Somma-Vesuvius (Italy): crystal chemistry and genetic inferences. Lithos 160-161:84-97

Barmin A, Melnik O, Sparks RSJ (2002) Periodic behaviour in lava dome eruptions. Earth Planet Sci Lett 199:173-184

Brugger CR, Hammer JE (2010a) Crystal size distribution analysis of plagioclase in experimentally decompressed hydrous rhyodacite magma. Earth Planet Sci Lett 300: 246-254

Brugger CR, Hammer JE (2010b) Crystallization kinetics in continuous decompression experiments: implications for interpreting natural magma ascent processes. J Petrol 51: 1941-1965

Budi-Santoso A, Lesage P, Dwiyono S, Sumarti S, Subandriyo, Surono, Jousset P, Metaxian J-P (2013) Analysis of the seismic activity associated with the 2010 eruption of Merapi Volcano, Java. J Volcanol Geotherm Res 261:153-170

Burgisser A, Arbaret L, Druitt TH, Giachetti T (2011) Pre-explosive conduit conditions of the 1997 Vulcanian explosions at Soufrière Hills Volcano, Montserrat: II. Overpressure and depth distributions. J Volcanol Geotherm Res 199:193-205

Cashman KV (1988) Crystallization of Mount St. Helens 1980-1986 dacite: a quantitative textural approach. Bull Volcanol 50:194-209

Cashman KV (1992) Groundmass crystallization of Mount St. Helens dacite, 1980-1986: a tool for interpreting shallow magmatic processes. Contrib Mineral Petrol 109:431-449

Cashman KV, McConnell SM (2005) Multiple levels of magma storage during the 1980 summer eruptions of Mount St. Helens, WA. Bull Volcanol 68:57-75

Castro JM, Gardner JE (2008) Did magma ascent rate control the explosive-effusive transition at the Inyo volcanic chain, California? Geology 36:279-282

Castro JM, Cashman KV, Manga M (2003) A technique for measuring 3D crystal-size distributions of prismatic microlites in obsidian. Am Mineral 88:1230-1240

Charbonnier SJ, Gertisser R (2008) Field observations and surface characteristics of pristine block-and-ash flow deposits from the 2006 eruption of Merapi volcano, Java, Indonesia. J Volcanol Geotherm Res 177:971-982

Charbonnier SJ, Germa AM, Connor CB, Gertisser R, Preece K, Komorowski J-C, Lavigne F, Dixon TH, Connor LJ (2013) Evaluation of the impact of the 2010 pyroclastic density currents at Merapi volcano from high-resolution satellite imagery analysis, field investigations and numerical simulations. J Volcanol Geotherm Res 261:295-315

Clarke AB, Stephens S, Teasdale R, Sparks RSJ, Diller K (2007) Petrologic constraints on the decompression history of magma prior to Vulcanian explosions at the Soufrière Hills Volcano, Montserrat. J Volcanol Geotherm Res 161:261-274

Costa F, Andreastuti S, de Maisonneuve CB, Pallister JS (2013) Petrological insights into the storage conditions, and magmatic processes that yielded the centennial 2010 Merapi explosive eruption. J Volcanol Geotherm Res 261:209-235

Couch S, Harford CL, Sparks RSJ, Carroll MR (2003a) Experimental constraints on the conditions of formation of highly calcic plagioclase microlites at the Soufrière Hills volcano, Montserrat. J Petrol 44: 1455-1475

Couch S, Sparks RSJ Carroll MR (2003b) The kinetics of degassing-induced crystallization at Soufrière Hills volcano, Montserrat. J Petrol 44: 1477-1502

Cronin SJ, Lube G, Dayudi DS, Sumarti S, Subrandiyo S, Surono (2013) Insights into the October-November 210 Gunung Merapi eruption (Central Java, Indonesia) from the stratigraphy, volume and characteristics of its pyroclastic deposits. J Volcanol Geotherm Res 261: 244-259

Degruyter W, Bachmann O, Burgisser A, Manga M (2012) The effects of outgassing on the transition between effusive and explosive silicic eruptions. Earth Planet Sci Lett 349-350:161-170

Drignon M, Arbaret L, Burgisser A, Komorowski J-C, Martel C, Raditya P (2014) Transition from effusive to explosive activity during lava dome eruption: the example of the 2010 of Merapi Volcano (Java, Indonesia). AGU Fall Meeting 2014:V51B-4759

Gertisser R, Charbonnier SJ, Keller J, Quidelleur X (2012) The geological evolution of Merapi volcano, Indonesia. Bull Volcanol 74: 1213-1233

Gertisser R, Handley H, Reagan M, Preece K, Berlo K, Barclay J, Herd R (2015) Timescales of magma transfer, degassing and crustal assimilation at Merapi volcano, Indonesia. 25th Annual V.M. Goldschmidt Conference, Prague, Czech Republic, 2015.

Hammer JE, Rutherford MJ (2002) An experimental study of the kinetics of decompression-induced crystallization in silicic melt. J Geophys Res 107:1-23

Hammer JE, Cashman KV, Hoblitt RP, Newman S (1999) Degassing and microlite crystallization during pre-climactic events of the 1991 eruption of Mt. Pinatubo, Philippines. Bull Volcanol 60:355-380

Hammer JE, Cashman KV, Voight B (2000) Magmatic processes revealed by textural and compositional trends in Merapi dome lavas. J Volcanol Geotherm Res 100:165-192

Hess K-U, Dingwell DB (1996) Viscosities of hydrous leucogranite melts; a non-Arrhenian model. Am Mineral 81:1297-1300

Higgins MD (2000) Measurement of crystal size distributions. Am Mineral 85:1105-1116

Hoblitt RP, Harmon RS (1993) Bimodal density distribution of cryptodome dacite from the 1980 eruption of Mount St. Helens, Washington. Bull Volcanol 55:421-437

Horwell CJ, Williamson BJ, Llewellin EW, Damby DE, Le Blond JS (2013) The nature and formation of cristobalite at the Soufrière
Hills volcano, Montserrat: implications for the petrology and stability of silicic lava domes. Bull Volcanol 75:1–19

Humphreys MCS, Christopher T, Hards V (2009) Microlite transfer by dis-aggregation of mafic inclusion following magma mixing at Soufrière Hills volcano, Montserrat. Contrib Mineral Petrol 157:609–624

Jaupart C, Allègre CJ (1991) Gas content, eruption rate and instabilities of eruption regime in silicic volcanoes. Earth Planet Sci Lett 102:413–429

Jousset P, Budi-Santoso A, Jolly AD, Boichu M, Surono, Dwiyono S, Sunarti S, Hidayati S, Thierry P (2013) Signs of magma ascent in LP and VLP seismic events and link to degassing: an example from the 2010 explosive eruption at Merapi volcano, Indonesia. J Volcanol Geotherm Res 261:171–192

Koleszar AM, Kent AJR, Wallace PJ, Scott WE (2012) Controls on long-term explosivity at andesitic arc volcanoes: insights from Mount Hood, Oregon. J Volcanol Geotherm Res 219:220–21

Komorowski JC, Legendre Y, Christopher T, Bernstein M, Stewart R, Joseph E, Fournier N, Chardot L, Finizola A, Wadge G, Syers R, Williams C, Bass V (2010) Insights into processes and deposits of hazardous vulcanian explosions at Soufrière Hills Volcano during 2008 and 2009 (Montserrat, West Indies). Geophys Res Lett 37:L00E19

Komorowski JC, Jenkins S, Baxter PJ, Picquout A, Lavigne F, Charbonnier S, Gertisser R, Preece K, Cholik N, Budi-Santoso A, Surono (2013) Paroxysmal dome explosion during the Merapi 2010 eruption: processes and facies relationships of associated high-energy pyroclastic density currents. J Volcanol Geotherm Res 261:260–294

Lofgren G (1980) Experimental studies on the dynamic crystallization of silicate melts. In: Hargreaves RB (ed) Physics of magmatic processes. Princeton University Press, Princeton, pp 487–551

Marsh BD (1988) Crystal size distribution (CSD) in rocks and the kinetics and dynamics of crystallization I. Theory. Contrib Mineral Petrol 99: 277–291

Martel C (2012) Eruption dynamics inferred from microlite crystallization experiments: application to Plinian and dome-forming eruptions of Mt. Pelée (Martinique, Lesser Antilles). J Petrol 53:699–725

Martel C, Radadi Ali A, Poussineau S, Gourgaud A, Pichavant M (2006) Basalt-inherited microlites in silicic magmas: evidence from Mount Pelée (Martinique, French West Indies). Geology 34:905–908

Melnik O, Sparks RSJ (1999) Nonlinear dynamics of lava dome extrusion. Nature 402:37–41

Melnik O, Sparks RSJ (2005) Controls on conduit magma flow dynamics during lava dome building eruptions. J Geophys Res 110:B02209

Melnik OE, Blundy JD, Rust AC, Muir DD (2011) Subvolcanic plumbing systems imaged through crystal size distributions. Geology 39:403–406

Morgan DJ, Jerram DA (2006) On estimating crystal shape for crystal size distribution analysis. J Volcanol Geotherm Res 154:1–7

Murphy MD, Sparks RSJ, Barclay J, Carroll MR, Brewer TS (2000) Remobilization of andesite magma by intrusion of mafic magma at the Soufrière Hills volcano, Montserrat, West Indies. J Petrol 41:21–42

Newhall CG, Melson WG (1983) Explosive activity associated with the growth of volcanic domes. J Volcanol Geotherm Res 17:111–131

Ogburn SE, Loughlin SC, Calder ES (2015) The association of lava dome growth with major explosive activity (VEI ≥4): DomeHaz, a global dataset. Bull Volcanol 77:40

Pallister JS, Thornber CR, Cashman KV, Clynne MA, Lowers HA, Mandeville CW, Brownfield IK, Meeker GP (2008) Petrology of the 2004–2006 Mount St. Helens lava dome—implications for magmatic plumbing and eruption triggering. US Geol Surv Prof Pap 1750:647–702

Pallister JS, Schneider DJ, Griswold JP, Keeler RH, Burton WC, Noyles C, Newhall CG, Ratdomopurbo A (2013) Merapi 2010 eruption— chronology and extrusion rates monitored with satellite radar and used in eruption forecasting. J Volcanol Geotherm Res 261:144–152

Panigrahi MK, Naik RK, Pandit D, Misra KC (2008) Reconstructing physico-chemical parameters of hydrothermal mineralization of copper at the Malankjhand deposit, India, from mineral chemistry of biotite, chlorite and epidote. Geochem J 42:443–460

Papale P, Moretti R, Barbato D (2006) The compositional dependence of the saturation surface of H₂O + CO₂ fluids in silicate melts. Chem Geol 229:78–95

Preece K, Barclay J, Gertisser R, Herd RA (2013) Textural and micro-petrological variations in the eruptive products of the 2006 dome-forming eruption of Merapi volcano, Indonesia: implications for sub-surface processes. J Volcanol Geotherm Res 261:98–120

Preece KJ (2014) Transitions between effusive and explosive activity at Merapi volcano, Indonesia: a volcanological and petrological study of the 2006 and 2010 eruptions. Ph.D. Thesis, University of East Anglia

Preece K, Gertisser R, Barclay J, Berlo K, Herd RA, Edinburgh Ion Microprobe Facility (2014) Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia. Contributions to Mineralogy and Petrology 168:1061

Ratdomopurbo A, Beauducel F, Subandriyo J, Agung Nandaka IGM, Newhall CG, Suharna, Sayudi DS, Suparwaka H, Sunarta S (2013) Overview of the 2006 eruption of Mt. Merapi. J Volcanol Geotherm Res 261:97–97

Ridolfi F, Puerini M, Renzulli A, Menna M, Toulkeridis T (2008) The magmatic feeding system of El Reventador volcano (Sub-Andean zone, Ecuador) constrained by texture, mineralogy and thermobarometry of the 2002 erupted products. J Volcanol Geotherm Res 176:94–106

Ruprecht P, Bachmann O (2010) Pre-eruptive reheating during magma mixing at Quizapu volcano and the implications for the explosive-ness of silicic are volcanoes. Geology 38:919–922

Scandone R, Cashman KV, Malone SD (2007) Magma supply, magma ascent and style of volcanic eruptions. Earth Planet Sci Lett 253: 513–529

Scott JAJ (2012) Origin and evolution of the Santiaguito lava dome complex, Guatemala. PhD Thesis, University of Oxford, UK

Scott JAJ, Mather TA, Pyle DM, Rose WI, Chigna G (2012) The magmatic plumbing system beneath Santiaguito Volcano, Guatemala. J Volcanol Geotherm Res 237:238–54–68

Selby D, Nesbitt BE (2000) Chemical composition of biotite from the Casino porphyry Cu-Au-Mo mineralization, Yukon, Canada; evaluation of magmatic and hydrothermal fluid chemistry. Chem Geol 171:77–93

Sparks RSJ (1997) Causes and consequences of pressurisation in lava dome eruptions. Earth Planet Sci Lett 150:177–189

Stix J, Zapata JA, Calvache M, Cortés GP, Fischer TP, Gómez D, Narvaez L, Ordoñez M, Ortega A, Torres R, Williams SN (1993) A model of degassing at Galeras Volcano, Colombia, 1988–1993. Geology 21: 963–967

Surono, Jousset P, Pallister J, Boichu M, Buongiorno MF, Budisantoso A, Costa F, Andreastuti S, Prata F, Schneider D, Clarisse L, Humaida H, Sumarti S, Bignami C, Griswold J, Carn S, Oppenheimer C, Lavigne F (2012) The 2010 explosive eruption of Java's Merapi volcano—a '100-year' event. J Volcanol Geotherm Res 241:242–121–135

Villemant B, Mouatt J, Michel A (2008) Andesitic magma degassing investigated through H₂O vapour-melt partitioning of halogens at Soufrière Hills volcano, Montserrat, (Lesser Antilles). Earth Planet Sci Lett 269:212–229

Waters LE, Lange RA (2015) An updated calibration of the plagioclase-liquid hygrometer-thermometer applicable to basalts through rhyolites. Am Mineral 100:2172–2184

Wen S, Nekvasil H (1994) SOLVCALC: an interactive graphics program package for calculating the ternary feldspar solvus and for two-feldspar geothermometry. Comput Geosci 20:1025–1040

Wolf KJ, Eichelberger JC (1997) Syneruptive mixing, degassing, and crystallization at Redoubt Volcano, eruption of December, 1989 to May 1990. J Volcanol Geotherm Res 75:19–37

Wright HMN, Cashman KV, Rosi M, Cioni R (2007) Breadcrust bombs as indicators of Vulcanian eruption dynamics at Guagua Pichincha volcano, Ecuador. Bull Volcanol 69:281–300