Katie Preece - Academia.edu (original) (raw)

Papers by Katie Preece

Research paper thumbnail of The 2010 Pyroclastic Density Currents of Merapi Volcano, Central Java, Indonesia

ABSTRACT The 2010 pyroclastic density currents (PDC) at Merapi present a rare opportunity to coll... more ABSTRACT The 2010 pyroclastic density currents (PDC) at Merapi present a rare opportunity to collect a uniquely detailed dataset of the source, extent, lateral variations and impact of various PDC deposits on a densely populated area. Using traditional volcanological field-based methods and multi-temporal dataset of high-resolution satellite imagery, a total of 23 PDC events have been recognized, including 5 main channeled flows, 15 overbank flows derived from overspill and re-channelization of the main PDCs into adjacent tributaries and two main surge events. The 2010 PDC deposits covered an area of ~22.3 km2, unequally distributed between valley-filling (6.9%), overbank (22.4%) and ash-cloud surge deposits (71.7%). Their total estimated volume is ~36.3×106 m3, with ~50.2% of this volume accounting for valley-filling deposits, 39.3% for overbank deposits and 10.5% for ash-cloud surge deposits. The internal architecture and facies variations of the 2010 PDC deposit were investigated using data collected from 30 stratigraphic sections measured after one rainy season of erosion. The results show that complex, local-scale variations in flow dynamics and deposit architectures are apparent and that the main factors that control the propagation of the main flows and their potential hazards for overbanking were driven by: (1) the rapid emplacement of several voluminous PDCs, associated with the steady infilling of the receiving landscape after the two first phases of the eruption; (2) longitudinal changes in channel capacity following increased sinuosity in the valley and decreased containment space; and (3) the effects of varying generation mechanisms (gravitational dome collapse, vertical or lateral dome explosions and column-collapse) and source materials involved during individual PDC forming events. Integration of these data into numerical simulations of the 3-5 November channeled and overbank PDCs using two well-established geophysical mass flow models, Titan2D and Volcflow, allow us to evaluate the ability of these models to reproduce the main features of the natural deposits and some of the flow overbanking processes observed in the field. Using such a multi-techniques approach, the dataset obtained in this study is considered not only instrumental for characterizing PDCs and related hazards at Merapi, but will allow comparisons with similar events at other volcanoes around the globe.

Research paper thumbnail of Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

Contributions to Mineralogy and Petrology, 2014

clasts resulting from sub-plinian convective column collapse. This paper presents geochemical dat... more clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxenehosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H 2 O, CO 2 ) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have reequilibrated during shallower storage and/or ascent, at depths of ~0.6-9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform "buffered" Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km 3 in order to produce the mass of SO 2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing.

Research paper thumbnail of GPR Facies Analysis of Block-and-ash Flows at Merapi Volcano, Central Java, Indonesia

Merapi, a 2911-m-high basaltic andesite volcanic complex in Central Java is one of the most frequ... more Merapi, a 2911-m-high basaltic andesite volcanic complex in Central Java is one of the most frequently erupting volcanoes in Indonesia and best known for its nearly persistent volcanic activity characterized by the extrusion of viscous lava domes and collapse of these domes to produce block-and-ash flows (BAFs). During the most recent eruptive episode in 2006, BAFs affected the densely populated

Research paper thumbnail of Block-and-ash flow dynamics and hazards at Merapi volcano (Java, Indonesia): a ground penetrating radar approach

Merapi, a composite volcano in Central Java, is best known for its nearly persistent volcanic act... more Merapi, a composite volcano in Central Java, is best known for its nearly persistent volcanic activity characterized by the extrusion of viscous lava domes and collapse of these domes to produce block-and-ash flows (BAFs). BAFs generally follow narrow valleys resulting in distally extensive, yet spatially confined, deposits having significant depths, even at large distances from the source. As observed during

Research paper thumbnail of Overbank block-and-ash flow deposits and the impact of valley-derived, unconfined flows on populated areas at Merapi volcano, Java, Indonesia

Natural Hazards, 2012

Merapi, an andesitic volcanic complex in Central Java, is one of the most frequently erupting vol... more Merapi, an andesitic volcanic complex in Central Java, is one of the most frequently erupting volcanoes in Indonesia and poses a permanent threat to the surrounding population of over 1 million people. With frequently recurring volcanic activity, the sixty or so reported eruptions since the mid-1500s have caused *7,000 fatalities and destroyed numerous villages in the region. In June 2006, an eruption affected the densely populated area on the volcano's southern and south-eastern flanks for the first time in almost a century. The resultant block-and-ash flows (BAFs) travelled down an incised river valley (Kali Gendol) to a distance of 7 km from the source, breaking out of the main channel at four main locations. Unconfined (overbank) BAFs were generated, which covered the interfluve regions on either side of the main valley and buried buildings and other infrastructure features in the village of Kaliadem, situated on the western bank of the Gendol valley *5 km from the summit of Merapi. Using traditional volcanological field-based methods and non-invasive, high-resolution ground-penetrating radar techniques, the morphology and internal architecture of these overbank deposits were studied in detail in order to evaluate the destructive impact of these flows in a local context. The results show that complex, local-scale variations in flow dynamics and deposit architectures are apparent and that BAFs are capable of transporting significant numbers of large blocks ([1-2 m) out of the valley confines. We propose a conceptual model for the escape of these channelised BAFs onto the interfluvial terrace at Kaliadem and show, through a stratigraphic analysis of the pyroclastic successions underlying the village and adjacent areas on the volcano's southern flank, that the area has been affected repeatedly by overbank BAFs and explosive eruptions over the past few 100 years (and more).

Research paper thumbnail of Paroxysmal dome explosion during the Merapi 2010 eruption: Processes and facies relationships of associated high-energy pyroclastic density currents

Journal of Volcanology and Geothermal Research, 2013

An 11-minute sequence of laterally-directed explosions and retrogressive collapses on 5 November ... more An 11-minute sequence of laterally-directed explosions and retrogressive collapses on 5 November 2010 at Merapi (Indonesia) destroyed a rapidly-growing dome and generated high-energy pyroclastic density currents (PDCs) spreading over 22 km 2 with a runout of 8.4 km while contemporaneous co-genetic valley-confined PDCs reached 15.5 km. This event formed Stage 4 of the multi-stage 2010 eruption, the most intense eruptive episode at Merapi since 1872. The deposits and the widespread devastating impact of associated high-energy PDCs on trees and buildings show striking similarities with those from historical volcanic blasts (Montagne Pelée, Martinique, Bezymianny, Russia, Mount St. Helens, USA, Soufrière Hills, Montserrat). We provide data from stratigraphic and sedimentologic analyses of 62 sections of the first unequivocal blast-like deposits in Merapi's recent history. We used high resolution satellite imagery to map eruptive units and flow direction from the pattern of extensive tree blowdown. The stratigraphy of Stage 4 consists of three depositional units (U0, U1, U2) that we correlate to the second, third and fourth explosions of the seismic record. Both U1 and U2 show a bi-partite layer stratigraphy consisting each of a lower L1 layer and an upper L2 layer. The lower L1 layer is typically very coarse-grained, fines-poor, poorly-sorted and massive, and was deposited by the erosive waxing flow head. The overlying L2 layer is much finer grained, fines-rich, moderately to well-sorted, with laminar to wavy stratification. L2 was deposited from the waning upper part and wake of the PDC. Field observations indicate that PDC height reached~330 m with an internal velocity of~100 m s −1 within 3 km from the source. The summit's geometry and the terrain morphology formed by a major transversal ridge and a funneling deep canyon strongly focused PDC mass towards a major constriction, thereby limiting the loss of kinetic energy. This favored elevated PDC velocities and high particle concentration, promoted overspilling of PDCs across high ridges into other river valleys, and generated significant dynamic pressures to distances of 6 km that caused total destruction of buildings and the forest. The Merapi 2010 eruption highlights that explosive and gravitational disintegration of a rapidly growing dome can generate devastating high-energy, high-velocity PDCs. This constitutes a credible high impact scenario for future multi-stage eruptions at Merapi and at other volcanoes that pose particular monitoring, crisis response, and risk reduction challenges.

Research paper thumbnail of Magmatic differentiation processes at Merapi Volcano: inclusion petrology and oxygen isotopes

Journal of Volcanology and Geothermal Research, 2013

Indonesian volcano Merapi is one of the most hazardous volcanoes on the planet and is characteris... more Indonesian volcano Merapi is one of the most hazardous volcanoes on the planet and is characterised by periods of active dome growth and intermittent explosive events. Merapi currently degasses continuously through high temperature fumaroles and erupts basaltic-andesite dome lavas and associated block-and-ash-flows that carry a large range of magmatic, coarsely crystalline plutonic, and meta-sedimentary inclusions. These inclusions are useful in order to evaluate magmatic processes that act within Merapi's plumbing system, and to help an assessment of which phenomena could trigger explosive eruptions. With the aid of petrological, textural, and oxygen isotope analysis we record a range of processes during crustal magma storage and transport, including mafic recharge, magma mixing, crystal fractionation, and country rock assimilation. Notably, abundant calc-silicate inclusions (true xenoliths) and elevated δ 18 O values in feldspar phenocrysts from 1994, 1998, 2006, and 2010 Merapi lavas suggest addition of limestone and calc-silicate materials to the Merapi magmas. Together with high δ 13 C values in fumarole gas, crustal additions to mantle and slab-derived magma and volatile sources are likely a steady state process at Merapi. This late crustal input could well represent an eruption trigger due to sudden over-pressurisation of the shallowest parts of the magma storage system independently of magmatic recharge and crystal fractionation. Limited seismic precursors may be associated with this type of eruption trigger, offering a potential explanation for the sometimes erratic behaviour of Merapi during volcanic crises.

Research paper thumbnail of Textural and micro-petrological variations in the eruptive products of the 2006 dome-forming eruption of Merapi volcano, Indonesia: Implications for sub-surface processes

Journal of Volcanology and Geothermal Research, 2013

The interplay between magma ascent, degassing and changing magmatic properties are widely recogni... more The interplay between magma ascent, degassing and changing magmatic properties are widely recognized as critical factors controlling the style of silicic volcanic eruptions. Microlite textures in samples from the prolonged dome-forming eruption of Merapi in 2006 provide a record of changing magmatic ascent conditions and shallow conduit processes throughout the eruption. Analysis of microlite textural parameters, including measurements of areal number density (N A ), mean microlite size, crystal aspect ratio and groundmass crystallinity (φ), combined with the monitoring record and field observations, indicate that magma ascent paths change between continuous ascent at varying rates from a deeper magma storage region, to ascent being temporarily stalled at shallow depths in the latter stages of the eruption, supporting the idea of an ephemeral shallow magma storage region at Merapi. Plagioclase microlite compositions show evidence of decompression-induced degassing, often displaying rims of anorthoclase and more K-rich alkali feldspar (sanidine). Anorthite contents also support the textural data of later erupted magma being temporarily stalled at shallow depths. Crystal size distributions (CSDs) are interpreted to show that both growth-dominated and nucleation-dominated crystallisation regimes existed during the 2006 eruption, resulting from changing conditions of undercooling (ΔT) during variable magma ascent paths. By contrast, microlite textural analysis and feldspar microlite compositions of samples from the fast-growing lava dome of the second phase of the 2010 eruption prior to the cataclysmic events on 5 November indicate faster ascent rates, a crystallisation regime more strongly dominated by nucleation due to high ΔT and interaction of the 2010 magma with more hotter magma from depth.

Research paper thumbnail of The pre-eruptive magma plumbing system of the 2007–2008 dome-forming eruption of Kelut volcano, East Java, Indonesia

Contributions to Mineralogy and Petrology, 2013

Kelut volcano, East Java, is an active volcanic complex hosting a summit crater lake that has bee... more 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 (An 47-94 ), orthopyroxene (En 64-72 , Fs 24-32 , Wo 2-4 ), clinopyroxene (En 40-48 , Fs 14-19 , Wo 34-46 ), Ti-magnetite (Usp 16-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 H 2 O contents ranging from *8.1 to 8.6 wt.% in the lower crustal system to *1.5 to Communicated by J. Hoefs.

Research paper thumbnail of Uranium-Series Isotopic Constraints on Recent Changes in the Eruptive Behaviour of Merapi Volcano, Java, Indonesia

ABSTRACT Merapi volcano (Central Java) is one of the most active and deadly volcanoes in Indonesi... more ABSTRACT Merapi volcano (Central Java) is one of the most active and deadly volcanoes in Indonesia. The 2010 eruption was the volcano's largest eruption since 1872 and erupted much more violently than expected. Prior to 2010, volcanic activity at Merapi was characterised by several months of slow dome growth punctuated by gravitational dome failures, generating small-volume pyroclastic density currents (Merapi-type nuées ardentes). The unforeseen, large-magnitude events in 2010 were different in many respects: pyroclastic density currents travelled > 15 km beyond the summit causing widespread devastation in proximal areas on Merapi's south flank and ash emissions from sustained eruption columns resulted in ash fall tens of kilometres away from the volcano. The 2010 events have proved that Merapi's relatively small dome-forming activity can be interrupted at relatively short notice by larger explosive eruptions, which appear more common in the geological record. We present new geochemical and Uranium-series isotope data for the volcanic products of both the 2006 and 2010 eruptions at Merapi to investigate the driving forces behind this unusual explosive behaviour and their timescales. An improved knowledge of these processes and of changes in the pre-eruptive magma system has important implications for the assessment of hazards and risks from future eruptive activity at Merapi.

Research paper thumbnail of Merapi (Java, Indonesia): anatomy of a killer volcano

Geology Today, 2011

Merapi is Indonesia's most dangerous volcano with a history of deadly eruptions. Over the past tw... more Merapi is Indonesia's most dangerous volcano with a history of deadly eruptions. Over the past two centuries, the volcanic activity has been dominated by prolonged periods of lava dome growth and intermittent gravitational or explosive dome failures to produce pyroclastic flows every few years. Explosive eruptions, such as in 2010, have occurred occasionally during this period, but were more common in pre-historical time, during which a collapse of the western sector of the volcano occurred at least once. Variations in magma supply from depth, magma ascent rates and the degassing behaviour during ascent are thought to be important factors that control whether Merapi erupts effusively or explosively. A combination of subsurface processes operating at relatively shallow depth inside the volcano, including complex conduit processes and the release of carbon dioxide into the magmatic system through assimilation of carbonate crustal rocks, may result in unpredictable explosive behaviour during periods of dome growth. Pyroclastic flows generated by gravitational or explosive lava dome collapses and subsequent lahars remain the most likely immediate hazards near the volcano, although the possibility of more violent eruptions that affect areas farther away from the volcano cannot be fully discounted. In order to improve hazard assessment during future volcanic crises at Merapi, we consider it crucial to improve our understanding of the processes operating in the volcano's plumbing system and their surface manifestations, to generate accurate hazard zonation maps that make use of numerical mass flow models on a realistic digital terrain model, and to utilize probabilistic information on eruption recurrence and inundation areas.

Research paper thumbnail of Magma evolution and volatile behaviour prior to and during the Plinian Lower Pumice 2 eruption, Santorini, Greece

The Plinian Lower Pumice 2 (LP2) eruption (c. 180 ka) was the first major caldera-forming eruptio... more The Plinian Lower Pumice 2 (LP2) eruption (c. 180 ka) was the first major caldera-forming eruption of the Santorini volcanic complex (Greece). The eruption shows some striking similarities to the caldera-forming Late Bronze Age (Minoan) eruption in terms of field, petrological and geochemical characteristics of its eruptive products, which are discussed here to reveal the storage conditions of the LP2

Research paper thumbnail of Using Ground-Penetrating Radar to Unravel the Dynamics and Hazards of Block-and-ash Flows at Merapi Volcano, Indonesia

Assessing the dynamics and depositional processes of block-and-ash flows (BAFs) can be extremely ... more Assessing the dynamics and depositional processes of block-and-ash flows (BAFs) can be extremely difficult when based on traditional field-based investigations alone, as these are often complicated by a combination of poor exposure, rapid lateral facies variations of associated deposits and unknown paleotopography. However, when combined with other techniques such as ground-penetrating radar (GPR), a method that allows the internal architecture

Research paper thumbnail of GPR facies analysis of block-and-ash flows, Merapi volcano, Central Java (Indonesia): new insights into deposit architecture and emplacement mechanisms

Research paper thumbnail of Evaluation of the impact of the 2010 pyroclastic density currents at Merapi volcano from high-resolution satellite imagery, field investigations and numerical simulations

Research paper thumbnail of The 2010 Eruption of Merapi Volcano, Java, Indonesia: Petrological Insights into Magma Dynamics and Eruptive Behaviour

ABSTRACT The violent eruption of Merapi volcano (Central Java) that started on 26 October 2010 wa... more ABSTRACT The violent eruption of Merapi volcano (Central Java) that started on 26 October 2010 was the volcano's largest since 1872 and the deadliest event since 1930. Before 2010, Merapi's more recent (historical) eruptive activity was repeatedly characterised by periods of slow lava dome extrusion punctuated by gravitational dome failures, generating small-volume pyroclastic density currents (PDCs) with runout distances of typically less than 10 km. The unforeseen, large-magnitude events in 2010 were unusual in many respects: (1) the eruption was short-lived and started with an explosive phase that was not preceded by a lava dome at the surface; (2) between 31 October and 4 November, a lava dome appeared and grew rapidly within the summit crater, exceeding growth rates observed at the peak of the penultimate eruption in 2006 by a factor of ~ 22; (3) during the most vigorous eruptive phase on 5 November, at least one PDC travelled more than 15 km (more than twice the distance of the largest flows in 2006) beyond the summit along the Gendol river valley, causing widespread devastation on Merapi's south flank; (4) in a late phase of the eruption, pumice-rich PDCs were produced, forming a thin veneer on top of the deposits of the largest PDCs from 5 November; (5) ash emissions from sustained eruption columns resulted in ash fall tens of kilometres from the volcano, affecting, amongst other areas, the volcano's western slopes and the city of Yogyakarta ~ 25 km to the south; and (6) the total deposited volume in 2010, based on provisional estimates, may have been ~ 10 times higher than that of other recent eruptions. Here we report and present new geochemical, Sr-Nd-O isotope and U-series data for the volcanic products (lava dome fragments, magmatic inclusions, scoria, pumice and ash) from various stages of the 2010 eruption of Merapi. These data are discussed in the context of other recent to historical, typically less explosive, dome-forming eruptions to elucidate the driving forces behind the unusual explosive behaviour in 2010 and their timescales. The 2010 events highlight that dome extrusion and relatively small, prolonged dome-forming eruptions at Merapi can be interrupted at relatively short notice by larger and more vigorous eruption phases or eruptions, which appear more common in the geological record. An improved knowledge of these processes and of changes in the pre-eruptive magma system has important implications for the assessment of hazards and risks from future eruptive activity at Merapi.

Research paper thumbnail of The Plinian Lower Pumice 2 eruption, Santorini, Greece: Magma evolution and volatile behaviour

Journal of Volcanology and Geothermal Research, 2009

The Plinian Lower Pumice 2 (LP2) eruption (172 ka) was one of the first major caldera-forming eru... more The Plinian Lower Pumice 2 (LP2) eruption (172 ka) was one of the first major caldera-forming eruptions of the Santorini volcanic complex (Greece). The eruption shows some striking similarities to the caldera-forming Late Bronze Age (Minoan) eruption in terms of field, petrological and geochemical characteristics of its eruptive products, which are used to reveal the storage conditions of the LP2 magmas, pre-eruptive magmatic processes and the behaviour and degassing of volatiles prior to and during eruption.The LP2 eruption comprises four, predominantly rhyodacitic eruptive units (LP2-A, B, C, D). The lowermost unit of the Plinian LP2 deposits (LP2-A) consists of a basal phreatomagmatic bed (LP2-A1), which is overlain by three discrete pumice fall deposits (LP2-A2-1, A2-2, A3), the most prominent of which (LP2-A3) contains abundant, quench-textured scoriae that range in composition from basalt to basaltic andesite. The eruption proceeded with the deposition of pumice-rich pyroclastic flows (LP2-B) characterised by a lower, stratified and cross-bedded ignimbrite (LP2-B1) that may grade into a massive, non-welded ignimbrite (LP2-B2), a lithic-rich pumiceous breccia (LP2-C) and a co-ignimbrite lithic lag breccia (LP2-D).The main volume of rhyodacitic magma, which formed by fractionation of olivine, clinopyroxene, orthopyroxene, plagioclase, amphibole, Fe–Ti oxides, pyrrhotite and apatite from basaltic parental magmas and assimilation of crustal rocks, was held at mid-crustal levels (≤ 16 km depth), magmatic temperatures of 831 ± 12 °C and an oxygen fugacity slightly above the fayalite–magnetite–quartz (FMQ) oxygen buffer. Injection of ∼ 200 °C hotter mafic magma into the rhyodacitic reservoir and subsequent mingling and minor hybridisation with the resident magma helped to remobilise the rhyodacitic host magma and determined the final compositional range of the erupted products.Melt inclusion data show that sulphur concentrations were reduced to < 270 ppm in the rhyodacite, primarily due to partitioning of sulphur into pyrrhotite or, depending on temperature, a FeS-rich melt during magmatic differentiation at oxygen fugacities around the FMQ oxygen buffer. Sulphur concentrations in groundmass glasses of the LP2 pumices suggest that ∼ 43% of the remaining sulphur was released into the atmosphere during the LP2 eruption, the climatic effects of which are considered minor when compared to eruptions of more oxidised silicic arc magmas. Chlorine remained dissolved in the melt during magmatic differentiation prior to and during the LP2 eruption, indicating that chlorine emissions to the atmosphere were negligible.

Research paper thumbnail of The 2010 Pyroclastic Density Currents of Merapi Volcano, Central Java, Indonesia

ABSTRACT The 2010 pyroclastic density currents (PDC) at Merapi present a rare opportunity to coll... more ABSTRACT The 2010 pyroclastic density currents (PDC) at Merapi present a rare opportunity to collect a uniquely detailed dataset of the source, extent, lateral variations and impact of various PDC deposits on a densely populated area. Using traditional volcanological field-based methods and multi-temporal dataset of high-resolution satellite imagery, a total of 23 PDC events have been recognized, including 5 main channeled flows, 15 overbank flows derived from overspill and re-channelization of the main PDCs into adjacent tributaries and two main surge events. The 2010 PDC deposits covered an area of ~22.3 km2, unequally distributed between valley-filling (6.9%), overbank (22.4%) and ash-cloud surge deposits (71.7%). Their total estimated volume is ~36.3×106 m3, with ~50.2% of this volume accounting for valley-filling deposits, 39.3% for overbank deposits and 10.5% for ash-cloud surge deposits. The internal architecture and facies variations of the 2010 PDC deposit were investigated using data collected from 30 stratigraphic sections measured after one rainy season of erosion. The results show that complex, local-scale variations in flow dynamics and deposit architectures are apparent and that the main factors that control the propagation of the main flows and their potential hazards for overbanking were driven by: (1) the rapid emplacement of several voluminous PDCs, associated with the steady infilling of the receiving landscape after the two first phases of the eruption; (2) longitudinal changes in channel capacity following increased sinuosity in the valley and decreased containment space; and (3) the effects of varying generation mechanisms (gravitational dome collapse, vertical or lateral dome explosions and column-collapse) and source materials involved during individual PDC forming events. Integration of these data into numerical simulations of the 3-5 November channeled and overbank PDCs using two well-established geophysical mass flow models, Titan2D and Volcflow, allow us to evaluate the ability of these models to reproduce the main features of the natural deposits and some of the flow overbanking processes observed in the field. Using such a multi-techniques approach, the dataset obtained in this study is considered not only instrumental for characterizing PDCs and related hazards at Merapi, but will allow comparisons with similar events at other volcanoes around the globe.

Research paper thumbnail of Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

Contributions to Mineralogy and Petrology, 2014

clasts resulting from sub-plinian convective column collapse. This paper presents geochemical dat... more clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxenehosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H 2 O, CO 2 ) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have reequilibrated during shallower storage and/or ascent, at depths of ~0.6-9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform "buffered" Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km 3 in order to produce the mass of SO 2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing.

Research paper thumbnail of GPR Facies Analysis of Block-and-ash Flows at Merapi Volcano, Central Java, Indonesia

Merapi, a 2911-m-high basaltic andesite volcanic complex in Central Java is one of the most frequ... more Merapi, a 2911-m-high basaltic andesite volcanic complex in Central Java is one of the most frequently erupting volcanoes in Indonesia and best known for its nearly persistent volcanic activity characterized by the extrusion of viscous lava domes and collapse of these domes to produce block-and-ash flows (BAFs). During the most recent eruptive episode in 2006, BAFs affected the densely populated

Research paper thumbnail of Block-and-ash flow dynamics and hazards at Merapi volcano (Java, Indonesia): a ground penetrating radar approach

Merapi, a composite volcano in Central Java, is best known for its nearly persistent volcanic act... more Merapi, a composite volcano in Central Java, is best known for its nearly persistent volcanic activity characterized by the extrusion of viscous lava domes and collapse of these domes to produce block-and-ash flows (BAFs). BAFs generally follow narrow valleys resulting in distally extensive, yet spatially confined, deposits having significant depths, even at large distances from the source. As observed during

Research paper thumbnail of Overbank block-and-ash flow deposits and the impact of valley-derived, unconfined flows on populated areas at Merapi volcano, Java, Indonesia

Natural Hazards, 2012

Merapi, an andesitic volcanic complex in Central Java, is one of the most frequently erupting vol... more Merapi, an andesitic volcanic complex in Central Java, is one of the most frequently erupting volcanoes in Indonesia and poses a permanent threat to the surrounding population of over 1 million people. With frequently recurring volcanic activity, the sixty or so reported eruptions since the mid-1500s have caused *7,000 fatalities and destroyed numerous villages in the region. In June 2006, an eruption affected the densely populated area on the volcano's southern and south-eastern flanks for the first time in almost a century. The resultant block-and-ash flows (BAFs) travelled down an incised river valley (Kali Gendol) to a distance of 7 km from the source, breaking out of the main channel at four main locations. Unconfined (overbank) BAFs were generated, which covered the interfluve regions on either side of the main valley and buried buildings and other infrastructure features in the village of Kaliadem, situated on the western bank of the Gendol valley *5 km from the summit of Merapi. Using traditional volcanological field-based methods and non-invasive, high-resolution ground-penetrating radar techniques, the morphology and internal architecture of these overbank deposits were studied in detail in order to evaluate the destructive impact of these flows in a local context. The results show that complex, local-scale variations in flow dynamics and deposit architectures are apparent and that BAFs are capable of transporting significant numbers of large blocks ([1-2 m) out of the valley confines. We propose a conceptual model for the escape of these channelised BAFs onto the interfluvial terrace at Kaliadem and show, through a stratigraphic analysis of the pyroclastic successions underlying the village and adjacent areas on the volcano's southern flank, that the area has been affected repeatedly by overbank BAFs and explosive eruptions over the past few 100 years (and more).

Research paper thumbnail of Paroxysmal dome explosion during the Merapi 2010 eruption: Processes and facies relationships of associated high-energy pyroclastic density currents

Journal of Volcanology and Geothermal Research, 2013

An 11-minute sequence of laterally-directed explosions and retrogressive collapses on 5 November ... more An 11-minute sequence of laterally-directed explosions and retrogressive collapses on 5 November 2010 at Merapi (Indonesia) destroyed a rapidly-growing dome and generated high-energy pyroclastic density currents (PDCs) spreading over 22 km 2 with a runout of 8.4 km while contemporaneous co-genetic valley-confined PDCs reached 15.5 km. This event formed Stage 4 of the multi-stage 2010 eruption, the most intense eruptive episode at Merapi since 1872. The deposits and the widespread devastating impact of associated high-energy PDCs on trees and buildings show striking similarities with those from historical volcanic blasts (Montagne Pelée, Martinique, Bezymianny, Russia, Mount St. Helens, USA, Soufrière Hills, Montserrat). We provide data from stratigraphic and sedimentologic analyses of 62 sections of the first unequivocal blast-like deposits in Merapi's recent history. We used high resolution satellite imagery to map eruptive units and flow direction from the pattern of extensive tree blowdown. The stratigraphy of Stage 4 consists of three depositional units (U0, U1, U2) that we correlate to the second, third and fourth explosions of the seismic record. Both U1 and U2 show a bi-partite layer stratigraphy consisting each of a lower L1 layer and an upper L2 layer. The lower L1 layer is typically very coarse-grained, fines-poor, poorly-sorted and massive, and was deposited by the erosive waxing flow head. The overlying L2 layer is much finer grained, fines-rich, moderately to well-sorted, with laminar to wavy stratification. L2 was deposited from the waning upper part and wake of the PDC. Field observations indicate that PDC height reached~330 m with an internal velocity of~100 m s −1 within 3 km from the source. The summit's geometry and the terrain morphology formed by a major transversal ridge and a funneling deep canyon strongly focused PDC mass towards a major constriction, thereby limiting the loss of kinetic energy. This favored elevated PDC velocities and high particle concentration, promoted overspilling of PDCs across high ridges into other river valleys, and generated significant dynamic pressures to distances of 6 km that caused total destruction of buildings and the forest. The Merapi 2010 eruption highlights that explosive and gravitational disintegration of a rapidly growing dome can generate devastating high-energy, high-velocity PDCs. This constitutes a credible high impact scenario for future multi-stage eruptions at Merapi and at other volcanoes that pose particular monitoring, crisis response, and risk reduction challenges.

Research paper thumbnail of Magmatic differentiation processes at Merapi Volcano: inclusion petrology and oxygen isotopes

Journal of Volcanology and Geothermal Research, 2013

Indonesian volcano Merapi is one of the most hazardous volcanoes on the planet and is characteris... more Indonesian volcano Merapi is one of the most hazardous volcanoes on the planet and is characterised by periods of active dome growth and intermittent explosive events. Merapi currently degasses continuously through high temperature fumaroles and erupts basaltic-andesite dome lavas and associated block-and-ash-flows that carry a large range of magmatic, coarsely crystalline plutonic, and meta-sedimentary inclusions. These inclusions are useful in order to evaluate magmatic processes that act within Merapi's plumbing system, and to help an assessment of which phenomena could trigger explosive eruptions. With the aid of petrological, textural, and oxygen isotope analysis we record a range of processes during crustal magma storage and transport, including mafic recharge, magma mixing, crystal fractionation, and country rock assimilation. Notably, abundant calc-silicate inclusions (true xenoliths) and elevated δ 18 O values in feldspar phenocrysts from 1994, 1998, 2006, and 2010 Merapi lavas suggest addition of limestone and calc-silicate materials to the Merapi magmas. Together with high δ 13 C values in fumarole gas, crustal additions to mantle and slab-derived magma and volatile sources are likely a steady state process at Merapi. This late crustal input could well represent an eruption trigger due to sudden over-pressurisation of the shallowest parts of the magma storage system independently of magmatic recharge and crystal fractionation. Limited seismic precursors may be associated with this type of eruption trigger, offering a potential explanation for the sometimes erratic behaviour of Merapi during volcanic crises.

Research paper thumbnail of Textural and micro-petrological variations in the eruptive products of the 2006 dome-forming eruption of Merapi volcano, Indonesia: Implications for sub-surface processes

Journal of Volcanology and Geothermal Research, 2013

The interplay between magma ascent, degassing and changing magmatic properties are widely recogni... more The interplay between magma ascent, degassing and changing magmatic properties are widely recognized as critical factors controlling the style of silicic volcanic eruptions. Microlite textures in samples from the prolonged dome-forming eruption of Merapi in 2006 provide a record of changing magmatic ascent conditions and shallow conduit processes throughout the eruption. Analysis of microlite textural parameters, including measurements of areal number density (N A ), mean microlite size, crystal aspect ratio and groundmass crystallinity (φ), combined with the monitoring record and field observations, indicate that magma ascent paths change between continuous ascent at varying rates from a deeper magma storage region, to ascent being temporarily stalled at shallow depths in the latter stages of the eruption, supporting the idea of an ephemeral shallow magma storage region at Merapi. Plagioclase microlite compositions show evidence of decompression-induced degassing, often displaying rims of anorthoclase and more K-rich alkali feldspar (sanidine). Anorthite contents also support the textural data of later erupted magma being temporarily stalled at shallow depths. Crystal size distributions (CSDs) are interpreted to show that both growth-dominated and nucleation-dominated crystallisation regimes existed during the 2006 eruption, resulting from changing conditions of undercooling (ΔT) during variable magma ascent paths. By contrast, microlite textural analysis and feldspar microlite compositions of samples from the fast-growing lava dome of the second phase of the 2010 eruption prior to the cataclysmic events on 5 November indicate faster ascent rates, a crystallisation regime more strongly dominated by nucleation due to high ΔT and interaction of the 2010 magma with more hotter magma from depth.

Research paper thumbnail of The pre-eruptive magma plumbing system of the 2007–2008 dome-forming eruption of Kelut volcano, East Java, Indonesia

Contributions to Mineralogy and Petrology, 2013

Kelut volcano, East Java, is an active volcanic complex hosting a summit crater lake that has bee... more 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 (An 47-94 ), orthopyroxene (En 64-72 , Fs 24-32 , Wo 2-4 ), clinopyroxene (En 40-48 , Fs 14-19 , Wo 34-46 ), Ti-magnetite (Usp 16-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 H 2 O contents ranging from *8.1 to 8.6 wt.% in the lower crustal system to *1.5 to Communicated by J. Hoefs.

Research paper thumbnail of Uranium-Series Isotopic Constraints on Recent Changes in the Eruptive Behaviour of Merapi Volcano, Java, Indonesia

ABSTRACT Merapi volcano (Central Java) is one of the most active and deadly volcanoes in Indonesi... more ABSTRACT Merapi volcano (Central Java) is one of the most active and deadly volcanoes in Indonesia. The 2010 eruption was the volcano&#39;s largest eruption since 1872 and erupted much more violently than expected. Prior to 2010, volcanic activity at Merapi was characterised by several months of slow dome growth punctuated by gravitational dome failures, generating small-volume pyroclastic density currents (Merapi-type nuées ardentes). The unforeseen, large-magnitude events in 2010 were different in many respects: pyroclastic density currents travelled &gt; 15 km beyond the summit causing widespread devastation in proximal areas on Merapi&#39;s south flank and ash emissions from sustained eruption columns resulted in ash fall tens of kilometres away from the volcano. The 2010 events have proved that Merapi&#39;s relatively small dome-forming activity can be interrupted at relatively short notice by larger explosive eruptions, which appear more common in the geological record. We present new geochemical and Uranium-series isotope data for the volcanic products of both the 2006 and 2010 eruptions at Merapi to investigate the driving forces behind this unusual explosive behaviour and their timescales. An improved knowledge of these processes and of changes in the pre-eruptive magma system has important implications for the assessment of hazards and risks from future eruptive activity at Merapi.

Research paper thumbnail of Merapi (Java, Indonesia): anatomy of a killer volcano

Geology Today, 2011

Merapi is Indonesia's most dangerous volcano with a history of deadly eruptions. Over the past tw... more Merapi is Indonesia's most dangerous volcano with a history of deadly eruptions. Over the past two centuries, the volcanic activity has been dominated by prolonged periods of lava dome growth and intermittent gravitational or explosive dome failures to produce pyroclastic flows every few years. Explosive eruptions, such as in 2010, have occurred occasionally during this period, but were more common in pre-historical time, during which a collapse of the western sector of the volcano occurred at least once. Variations in magma supply from depth, magma ascent rates and the degassing behaviour during ascent are thought to be important factors that control whether Merapi erupts effusively or explosively. A combination of subsurface processes operating at relatively shallow depth inside the volcano, including complex conduit processes and the release of carbon dioxide into the magmatic system through assimilation of carbonate crustal rocks, may result in unpredictable explosive behaviour during periods of dome growth. Pyroclastic flows generated by gravitational or explosive lava dome collapses and subsequent lahars remain the most likely immediate hazards near the volcano, although the possibility of more violent eruptions that affect areas farther away from the volcano cannot be fully discounted. In order to improve hazard assessment during future volcanic crises at Merapi, we consider it crucial to improve our understanding of the processes operating in the volcano's plumbing system and their surface manifestations, to generate accurate hazard zonation maps that make use of numerical mass flow models on a realistic digital terrain model, and to utilize probabilistic information on eruption recurrence and inundation areas.

Research paper thumbnail of Magma evolution and volatile behaviour prior to and during the Plinian Lower Pumice 2 eruption, Santorini, Greece

The Plinian Lower Pumice 2 (LP2) eruption (c. 180 ka) was the first major caldera-forming eruptio... more The Plinian Lower Pumice 2 (LP2) eruption (c. 180 ka) was the first major caldera-forming eruption of the Santorini volcanic complex (Greece). The eruption shows some striking similarities to the caldera-forming Late Bronze Age (Minoan) eruption in terms of field, petrological and geochemical characteristics of its eruptive products, which are discussed here to reveal the storage conditions of the LP2

Research paper thumbnail of Using Ground-Penetrating Radar to Unravel the Dynamics and Hazards of Block-and-ash Flows at Merapi Volcano, Indonesia

Assessing the dynamics and depositional processes of block-and-ash flows (BAFs) can be extremely ... more Assessing the dynamics and depositional processes of block-and-ash flows (BAFs) can be extremely difficult when based on traditional field-based investigations alone, as these are often complicated by a combination of poor exposure, rapid lateral facies variations of associated deposits and unknown paleotopography. However, when combined with other techniques such as ground-penetrating radar (GPR), a method that allows the internal architecture

Research paper thumbnail of GPR facies analysis of block-and-ash flows, Merapi volcano, Central Java (Indonesia): new insights into deposit architecture and emplacement mechanisms

Research paper thumbnail of Evaluation of the impact of the 2010 pyroclastic density currents at Merapi volcano from high-resolution satellite imagery, field investigations and numerical simulations

Research paper thumbnail of The 2010 Eruption of Merapi Volcano, Java, Indonesia: Petrological Insights into Magma Dynamics and Eruptive Behaviour

ABSTRACT The violent eruption of Merapi volcano (Central Java) that started on 26 October 2010 wa... more ABSTRACT The violent eruption of Merapi volcano (Central Java) that started on 26 October 2010 was the volcano&#39;s largest since 1872 and the deadliest event since 1930. Before 2010, Merapi&#39;s more recent (historical) eruptive activity was repeatedly characterised by periods of slow lava dome extrusion punctuated by gravitational dome failures, generating small-volume pyroclastic density currents (PDCs) with runout distances of typically less than 10 km. The unforeseen, large-magnitude events in 2010 were unusual in many respects: (1) the eruption was short-lived and started with an explosive phase that was not preceded by a lava dome at the surface; (2) between 31 October and 4 November, a lava dome appeared and grew rapidly within the summit crater, exceeding growth rates observed at the peak of the penultimate eruption in 2006 by a factor of ~ 22; (3) during the most vigorous eruptive phase on 5 November, at least one PDC travelled more than 15 km (more than twice the distance of the largest flows in 2006) beyond the summit along the Gendol river valley, causing widespread devastation on Merapi&#39;s south flank; (4) in a late phase of the eruption, pumice-rich PDCs were produced, forming a thin veneer on top of the deposits of the largest PDCs from 5 November; (5) ash emissions from sustained eruption columns resulted in ash fall tens of kilometres from the volcano, affecting, amongst other areas, the volcano&#39;s western slopes and the city of Yogyakarta ~ 25 km to the south; and (6) the total deposited volume in 2010, based on provisional estimates, may have been ~ 10 times higher than that of other recent eruptions. Here we report and present new geochemical, Sr-Nd-O isotope and U-series data for the volcanic products (lava dome fragments, magmatic inclusions, scoria, pumice and ash) from various stages of the 2010 eruption of Merapi. These data are discussed in the context of other recent to historical, typically less explosive, dome-forming eruptions to elucidate the driving forces behind the unusual explosive behaviour in 2010 and their timescales. The 2010 events highlight that dome extrusion and relatively small, prolonged dome-forming eruptions at Merapi can be interrupted at relatively short notice by larger and more vigorous eruption phases or eruptions, which appear more common in the geological record. An improved knowledge of these processes and of changes in the pre-eruptive magma system has important implications for the assessment of hazards and risks from future eruptive activity at Merapi.

Research paper thumbnail of The Plinian Lower Pumice 2 eruption, Santorini, Greece: Magma evolution and volatile behaviour

Journal of Volcanology and Geothermal Research, 2009

The Plinian Lower Pumice 2 (LP2) eruption (172 ka) was one of the first major caldera-forming eru... more The Plinian Lower Pumice 2 (LP2) eruption (172 ka) was one of the first major caldera-forming eruptions of the Santorini volcanic complex (Greece). The eruption shows some striking similarities to the caldera-forming Late Bronze Age (Minoan) eruption in terms of field, petrological and geochemical characteristics of its eruptive products, which are used to reveal the storage conditions of the LP2 magmas, pre-eruptive magmatic processes and the behaviour and degassing of volatiles prior to and during eruption.The LP2 eruption comprises four, predominantly rhyodacitic eruptive units (LP2-A, B, C, D). The lowermost unit of the Plinian LP2 deposits (LP2-A) consists of a basal phreatomagmatic bed (LP2-A1), which is overlain by three discrete pumice fall deposits (LP2-A2-1, A2-2, A3), the most prominent of which (LP2-A3) contains abundant, quench-textured scoriae that range in composition from basalt to basaltic andesite. The eruption proceeded with the deposition of pumice-rich pyroclastic flows (LP2-B) characterised by a lower, stratified and cross-bedded ignimbrite (LP2-B1) that may grade into a massive, non-welded ignimbrite (LP2-B2), a lithic-rich pumiceous breccia (LP2-C) and a co-ignimbrite lithic lag breccia (LP2-D).The main volume of rhyodacitic magma, which formed by fractionation of olivine, clinopyroxene, orthopyroxene, plagioclase, amphibole, Fe–Ti oxides, pyrrhotite and apatite from basaltic parental magmas and assimilation of crustal rocks, was held at mid-crustal levels (≤ 16 km depth), magmatic temperatures of 831 ± 12 °C and an oxygen fugacity slightly above the fayalite–magnetite–quartz (FMQ) oxygen buffer. Injection of ∼ 200 °C hotter mafic magma into the rhyodacitic reservoir and subsequent mingling and minor hybridisation with the resident magma helped to remobilise the rhyodacitic host magma and determined the final compositional range of the erupted products.Melt inclusion data show that sulphur concentrations were reduced to < 270 ppm in the rhyodacite, primarily due to partitioning of sulphur into pyrrhotite or, depending on temperature, a FeS-rich melt during magmatic differentiation at oxygen fugacities around the FMQ oxygen buffer. Sulphur concentrations in groundmass glasses of the LP2 pumices suggest that ∼ 43% of the remaining sulphur was released into the atmosphere during the LP2 eruption, the climatic effects of which are considered minor when compared to eruptions of more oxidised silicic arc magmas. Chlorine remained dissolved in the melt during magmatic differentiation prior to and during the LP2 eruption, indicating that chlorine emissions to the atmosphere were negligible.