Search USGS⌕ Search

Geology topics

Research about Indonesia

Source-linked reports with geographic coverage including Indonesia.

At least 19 recordsLinked to original sources
USGS2026

Seismic characteristics of the April 2024 eruption of Ruang Volcano, North Sulawesi, Indonesia

The two most recent eruptions of Ruang volcano, in 2002 and 2024, have been large (VEI 4) with very short precursory periods, 3 and 6 days, respectively. Background seismicity at Ruang volcano is generally low with 0–2 volcano-tectonic (VT) earthquakes per day. The 2024 eruptions were preceded by a notable increase in VT earthquakes from 11 to 15 April with a rapid escalation from 16 to 17 April in the form of increased VT swarms, and the appearance of low frequency earthquakes and short bursts of VT drumbeats shortly before the eruption. Similar VT seismic unrest was recorded in March 2015, May 2016, and April 2022, although none of these episodes culminated in an eruption. A comparison of seismic energy release between April 2022 and April 2024 shows distinct differences. In 2024, cumulative seismic energy reached its peak within 5 days, just before the onset of eruption on 17 April. In contrast, during the 2022-episode, energy release followed a more gradual increase over 8 days peaking on 16 April, then declining gradually without leading to eruption. Notably, low frequency (LF) and drumbeat earthquakes were absent in non-eruptive unrest episodes. The transition to proximal LFs in 2024 likely signals pressure fluctuations, intensified magma degassing, and conduit formation, acting as an important precursor to eruption. Additionally, the occurrence of short bursts of drumbeat seismicity likely reflects the onset of the shallow magma ascent and was an important short-term precursor to the eruption. Finally, we examine the potential influence of regional tectonic earthquakes on the escalation of precursory unrest in 2024.

North Sulawesi, Ruang volcano↗
USGS2024

Constraining magma storage conditions of the Toba magmatic system: A plagioclase and amphibole perspective

Silicic magma reservoirs are responsible for producing the largest explosive eruptions in the geologic record. Petrologic and geochronological data provide evidence for these systems spending substantial periods of time (10 4 –10 5 yrs) within the upper crust prior to eruption; however, the long-term thermochemical evolution of these systems is not fully understood, as existing petrologic data make it challenging to quantify the time interval a magmatic system has spent at certain temperatures, or its “thermal history”. Here, we investigate the 74 ka Youngest Toba Tuff (YTT), one of the largest explosive eruptions in the geologic record, to better constrain the long-term thermal evolution of its magmatic system. We combine forward models of Sr diffusion in plagioclase and hornblende, mineral thermometry, and pre-existing trace-element evolution models to quantify the thermochemical evolution of the YTT magmatic system. We find that plagioclase crystals record decades to centuries of storage at temperatures 750 C, while hornblende records up to 6200 years at the same temperatures. Hornblende crystallizes at temperatures around 820 C and adjusting our diffusion modeling to this temperature results in no more than 900 years at initial crystallization conditions. Combined with previous trace-element modeling work, these results indicate that although there was chemical diversity for long durations in the YTT magma system sufficient to produce unique composition eruptive products, the entire system was experiencing a relatively similar thermal history that did not allow for large bodies of eruptible magma to be present for long periods ( 10 2 –10 3 years). Rather, we suggest that magmas within the YTT magmatic system were stored for long durations at thermal conditions where they were uneruptible and only remobilized within a few centuries prior to eruption.

Sumatra↗
USGS2023

High-resolution InSAR reveals localized pre-eruptive deformation inside the crater of Agung Volcano, Indonesia

During a volcanic crisis, high-rate, localized deformation can indicate magma close to the surface, with important implications for eruption forecasting. However, only a few such examples have been reported, because frequent, dense monitoring is needed. High-resolution Synthetic Aperture Radar (SAR) is capable of achieving <1 m spatial resolution and sub-weekly revisit times, but is under-used. Here we use high-resolution satellite SAR imagery from COSMO-SkyMed, TerraSAR-X, and Sentinel-1 to detect intra-crater uplift preceding the November 2017 onset of eruptive activity at Agung, Indonesia. Processing the SAR imagery with an up-to-date, accurate, high-resolution digital elevation model was crucial for preventing aliasing of the deformation signal and for accurate georeferencing. We show that >15 cm of line-of-sight shortening occurred over a 400-by-400 m area on the crater floor in September-October 2017, accompanying a deep seismic swarm and flank dyke intrusion. We attribute the deformation to the pressurization of a shallow (<200 m deep) hydrothermal system by the injection of magmatic gases and fluids. We also observe a second pulse of intra-crater deformation of 3–5 cm within 4 days to 11 hr prior to the first phreatomagmatic eruption, which is consistent with interaction between the hydrothermal system and the ascending magma. This phreatomagmatic eruption created the central pathway used during the final stages of magma ascent. Our observations have important implications for understanding unrest and eruption forecasting, and demonstrate the potential of monitoring with high-resolution SAR.

Agung Volcano, Bali↗
USGS2023

Investigating spatio-temporal variability of initial 230Th/232Th in intertidal corals

One of the key factors in obtaining precise and accurate 230 Th ages of corals, especially for corals with ages less than a few thousand years, is the correction for non-radiogenic 230 Th based on an initial 230 Th/ 232 Th value ( 230 Th/ 232 Th 0 ). Studies that consider coral 230 Th/ 232 Th 0 values in intertidal environments are limited, and it is in these environments that corals have Th concentrations 100–1000 times greater than open-ocean corals. Here we present 66 reconstructed U–Th isochrons of modern and Holocene Porites lutea and lobata corals from throughout the Sumatran forearc islands (SFI). Mean square of weighted deviates ( MSWD ) is used to evaluate the variabilities of 28 location-specific 230 Th/ 232 Th0 values before the regional mean. Results show no significant spatial difference between the means of 3.1 ± 3.2 × 10 −6 for the northern SFI and 4.7 ± 1.3 × 10 −6 for the southern SFI. Using the weighted mean of 4.3 ± 2.5 × 10 −6 for all islands has the advantage that no data need be subjectively excluded in the calculation of reliable 230 Th ages, where a local initial value is unknown. There were no clear centennial to millennial time-scale changes in 230 Th/ 232 Th 0 in the past 2500 years. This study, therefore, suggests a reliable value of initial 230 Th/ 232 Th ratio in intertidal environments.

Sumatran forearc islands↗
USGS2023

Regional crustal structure of Indonesia from receiver functions

Characterizing the crustal structure of Indonesia is important to gain a better understanding of its geodynamic evolution and improve seismic hazard assessments in the area. However, a unified crustal model of the entire Indonesian region and its surroundings is lacking. We present new maps of crustal thickness and bulk V p /V s ratio in Indonesia and the surrounding area that are obtained using P-wave receiver functions at 36 seismic stations from several permanent regional networks. The measured crustal thickness varies from ∼24 km to ∼38 km. The thickest crust, ∼38 km, is beneath Flores Island, southern Maluku, and neighboring northernmost Australia, whereas the thinnest crust, ∼24 km, is found under eastern Malaysia. Thus, crustal thickness varies by ∼14 km (from ∼24 km to ∼38 km) despite the small changes in elevation at the measurement points. The V p /V s ratios are 1.79 &#xB1; "> ± 0.11, with high values (>1.85) found along the Banda-Sunda arc-trench system. We attribute these high values to: (1) the presence of mafic island arc and oceanic crust and (2) partial melting within this volcanic region, which causes a larger decrease in S-wave velocities compared with P-wave velocities. The comparison of the seismic properties of Indonesian island arc crust, particularly the V p /V s ratio, with laboratory measurements and the petrology of the exhumed Talkeetna island arc, Alaska, allows us to infer the crustal composition of Indonesian island arc crust.

Tectonophysics↗
USGS2021

Repeating earthquakes during multiple phases of unrest and eruption at Mount Agung, Bali, Indonesia, 2017

In 2017, Mount Agung produced a small (VEI 2) eruption that was preceded by an energetic volcano-tectonic (VT) swarm (>800 earthquakes per day up to M4.9) and two months of declining activity. The period of decreased seismic activity complicated forecasting efforts for scientists monitoring the volcano. We examine the time history of earthquake families at Mount Agung in search of additional insight into the temporal changes in the shallow crust prior to eruption. Specifically, we analyze the period of declining seismic activity about five weeks prior to the eruption when forecasting uncertainty was greatest. We use REDPy (Hotovec-Ellis and Jeffries, 2016) to build a catalog of 6,508 earthquakes from 18 October 2017–15 February 2018 and group them into families of repeating earthquakes based on waveform similarity using a cross-correlation coefficient threshold of 0.8. We show that the evolution of earthquake families provides evidence that Mount Agung was progressing toward eruption even though overall earthquake rates and seismic-energy-release declined. We find that earthquake families that dominated seismicity during the beginning of the crisis ceased near the onset of tremor on 12 November 2017. Then, earthquake families took on characteristics commonly observed during effusive phases of eruptions on 15 November—a full six days before the first phreatomagmatic eruption on 21 November 2017 and a full ten days before the actual onset of lava effusion on 25 November 2017. We interpret the transitions in seismicity as the manifestation of a three-phase physical model including an Intrusion Phase, a Transition Phase, and a Eruptive Phase. During the Intrusion Phase, seismicity was dominated by VT earthquakes with a relatively high percentage of repeaters (59%) grouped into numerous (65) simultaneous families. During the Eruptive Phase, seismicity included both VT and low frequency earthquakes that grouped into relatively long-lived families despite a low overall percentage of repeaters (14%). The Transition Phase exhibited characteristics of earthquake families between the Intrusion Phase and Eruptive Phase. We conclude that the time history of earthquake families provides insight into the evolution of the stress distribution in the volcanic edifice, the development of the volcanic conduit, and seismogenesis of magma effusion. Finally, we discuss the role that repeating earthquakes could play in real-time monitoring at restless volcanoes. Our work suggests eruption forecasts can be improved by incorporating automatic processing codes to assist seismologists during sustained periods of high earthquake rates, even at sparsely monitored volcanoes.

Bali, Mount Agung↗
USGS2020

Local magnitude, coda magnitude, and radiated energy of volcanic tectonic earthquakes from October 2010 to December 2011 at Sinabung volcano, Indonesia

In August 2010, Sinabung volcano began erupting after more than a thousand years of dormancy. Following several weeks of phreatic eruptions, the eruptions ceased and Sinabung entered what became an inter-eruptive period of dominantly seismic unrest. While standard equations for understanding the size of an earthquake (local magnitude (M L ), coda magnitude (M C ), and seismic energy release (E R )) have long been developed, it is best practice to fine tune these relations for a given region and period of study to more accurately describe seismicity and to directly compare it with other volcanic systems. More accurate descriptions of magnitudes and energy release are vital to accurate volcanic eruption forecasting and evaluation of seismic and volcanic risk. In this study, we use high-frequency volcano-tectonic (VT) earthquakes recorded on a temporary three-component network installed between October 2010 and December 2011 in the region around Sinabung volcano to better constrain the seismic parameters of and better understand this previously unstudied volcano. We determine region-specific formulas for M L , M C , and E R as follows: M L = log 10 A + 1.1252 log 10 r + 0.0280 r − 2.5427 , M C = 0.7764 log 10 t c o d a + 0.0676 r − 0.7185 , a n d log 10 ( E R ) = 1.5720 M L + 11.5258 , where A , r , and t coda are maximum amplitude on a Wood-Anderson seismogram, hypocentral distance (km), and the coda duration (s), respectively. Constants in the M L equation have physically interpretable meanings. The constant for the geometrical spreading term (log 10 r term) equals one for perfect spherical spreading of the waveform. Our value is greater than one and thus suggests that wavefronts spread at a slightly different rate than for simple spherical spreading. The constant for the attenuation term ( r term) is consistent with locally mapped attenuative deposits (limestones and tuffs) and previous 3D tomographic results. Our M C equation differs from a previous study, likely because different data in a different time period were used. Earthquake hypocenters are consistent with those located in previous tomographic studies, and we interpret the earthquakes in this study as distal VT earthquakes induced by continued magmatic intrusion at Sinabung over the period of October 2010–December 2011.

Sinabung volcano↗
USGS2020

Relocated aftershocks and background seismicity in eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences

High seismicity rates in eastern Indonesia occur due to the complex interaction of several tectonic plates which resulted in two deadly, destructive earthquake sequences that occurred in Lombok Island and the city of Palu, Sulawesi in 2018. The first sequence began in July with an M w 6.4 event near Lombok, culminating in an M w 7.0 event 8 d later. This was then followed by a nearby M w 6.9 event 12 d later. Approximately 1000 km to the northeast, a separate sequence began several weeks later near Palu where an M w 7.5 event occurred that triggered a tsunami. In this study, we present hypocentre relocations for both earthquake sequences as well as all other regional earthquakes in eastern Indonesia. The relocations were performed using a teleseismic double-difference relocation method and arrival times for P and S waves from stations at local, regional, and teleseismic distances. The catalogue and phase data were taken from the Agency for Meteorology, Climatology and Geophysics (BMKG) of Indonesia and the International Seismological Centre (ISC) for the period of April 2009 through November 2018. The relocated catalogue provides an improved view of seismicity in eastern Indonesia over the study period, sharpening locations and interpretations of seismogenic features throughout the region. In the Lombok area, the relocated earthquakes clearly show a backarc thrust to the north of the Sunda-Banda Arc transition zone. The relocated aftershocks show that the destructive M w 7.0 and M w 6.9 earthquakes of the Lombok sequence ruptured two different regions: The M w 7.0 earthquake propagated westward, whereas the M w 6.9 earthquake propagated eastward. The entire sequence of Lombok earthquakes was most likely started by the M w 6.4 event as the initial event or foreshock, which then triggered backarc thrusts on both sides. Several weeks later and far to the northeast, the M w 7.5 Palu earthquake occurred along the Palu-Koro Fault, filling a seismic gap that had not ruptured in an M w 6.0 event or larger since at least 1900. The distribution of aftershocks indicates that the northern part of the Palu-Koro Fault has lower relative seismicity rates than the southern part at shallow depths, and that off fault aftershocks are mostly located to the east of the Palu-Koro Fault.

Lombok, Palu↗
USGS2020

An experimental investigation of interaction between andesite and hyperacidic volcanic lake water

Alteration in magmatic-hydrothermal systems leads to distinct changes in rock texture and mineralogy, and a strong redistribution of elements between fluid and rock. Here, we experimentally interacted andesite scoria with hyperacidic, high-sulfidation style fluids from Kawah Ijen volcano (Indonesia) at 25 and 100˚C, seeking to reproduce the textures observed in natural samples from this volcano, and to understand the element fluxes that accompany alteration. The susceptibility to alteration in the experiments is Cu-Fe-sulphide > calcic plagioclase > pyroxene > titano-magnetite > sodic plagioclase, with complete preservation of glass. Silicate minerals alter to opaline silica, and gypsum, barite and a Zr-phase precipitate. The selective alteration of the phenocryst minerals results in a preferential release of compatible elements, as the glass is the main incompatible element host. The experiments reproduce the alteration textures of the natural samples, including the preservation of glass, but the predicted compatible over incompatible element enrichment in the alteration element flux is not observed in the natural setting. This suggests that alteration at Kawah Ijen is dominated by lithologies that lack abundant glass, in particular lava flows where the glass has devitrified, despite these lava flows having a lower surface area compared to scoria.

Kawah Ijen↗
USGS2019

The 2017-19 activity at Mount Agung in Bali (Indonesia): Intense unrest, monitoring, crisis response, evacuation, and eruption

After 53 years of quiescence, Mount Agung awoke in August 2017, with intense seismicity, measurable ground deformation, and thermal anomalies in the summit crater. Although the seismic unrest peaked in late September and early October, the volcano did not start erupting until 21 November. The most intense explosive eruptions with accompanying rapid lava effusion occurred between 25 and 29 November. Smaller infrequent explosions and extrusions continue through the present (June 2019). The delay between intense unrest and eruption caused considerable challenges to emergency responders, local and national governmental agencies, and the population of Bali near the volcano, including over 140,000 evacuees. This paper provides an overview of the volcanic activity at Mount Agung from the viewpoint of the volcano observatory and other scientists responding to the volcanic crisis. We discuss the volcanic activity as well as key data streams used to track it. We provide evidence that magma intruded into the mid-crust in early 2017, and again in August of that year, prior to intrusion of an inferred dike between Mount Agung and Batur Caldera that initiated an earthquake swarm in late September. We summarize efforts to forecast the behavior of the volcano, to quantify exclusion zones for evacuations, and to work with emergency responders and other government agencies to make decisions during a complex and tense volcanic crisis.

Scientific Reports↗
USGS2019

The Orange Tuff: A Late Pleistocene tephra-fall deposit emplaced by a VEI 5 silicic Plinian eruption in West Java, Indonesia

A VEI 5 dacite eruption emplaced the Orange Tuff about between 34.3 cal kBP and 17.2 cal kBP. Gunung Salak is the unit’s source and the Orange Tuff represents the most recent such eruption from any of the volcanoes southwest of Bogor, Indonesia. The Orange Tuff is the region’s first such documented tephra-fall deposit whose characteristics and phenocryst geochemistry make it readily identifiable over at least 1250 km 2 . Magnetite compositions and temperature and f O 2 estimates inferred from Fe-Ti oxide compositions are particularly useful for identifying the unit. Deposit characteristics suggest that the eruption lasted 1–11 h with mass eruption rates of 1.0–8.3 × 10 8 kg/s and a column height of 31–40 km. The eruption’s column height and the deposit’s 2.5–11 km 3 volume suggest that the unit was dispersed over a much wider area than mapped. The unit is a marker bed throughout its mapped distribution and has potential to be applied over a much broader area as a regional marker bed. The large population and infrastructure proximal to Salak suggest that the unit should be considered in hazards assessments despite its age and the lack of subsequent similar eruptions.

Bulletin of Volcanology↗
USGS2019

Diverse late‐stage crystallization and storage conditions in melt domains from the Youngest Toba Tuff revealed by age and compositional heterogeneity in the last increment of accessory phase growth

The chemical record contained within the nal increment of growth on crystals is utilized to reveal the dynamics and time- scales of magma assembly and storage before eruption of the cataclysmic 2800 km3 Youngest Toba Tu (YTT), Indonesia. In situ U–Th disequilibrium dates and trace element concentrations were obtained via secondary ionization mass spectrometry (SIMS) on unsectioned and unpolished faces of individual zircon and allanite crystals. The six high-silica (> 73 wt% SiO2) pumices from which crystals were derived are among the more evolved and lower crystallinity (< 25 wt%) pumices from the YTT eruption, and likely represent the melt-dominated portion of the magma system. Discrete SIMS measurement cycles were coupled with statistical treatments to detect zircon and allanite surface zoning domains at the ~ 1 μm scale. Coupled r-MELTS and accessory phase saturation modeling indicates that at the granite ternary minimum or ‘eutectoid’ conditions that de ne this portion of the YTT, zircon and allanite crystallization is dependent on and proportionate to major phase crystallization, and is more limited than at pre-eutectoid conditions. A lower proportion of near-eruption zircon surface ages in the comparatively cool and wet YTT relative to other hotter and drier voluminous silicic eruptions could re ect the in u- ence of eutectoid storage conditions on magmatic responses to remobilization-related magmatic recharge.

Youngest Toba Tuff↗
USGS2019

Character of community response to volcanic crises at Sinabung and Kelud volcanoes

Volcano disaster mitigation is a priority in Indonesia's national development plan due to the large number of active volcanoes, high frequency of eruptions and dense population within volcanic hazards zones. As a consequence of these factors and the many decades of experience in mitigation of volcanic risk, Indonesia demonstrates that effective community response is critical to avoid disasters. Understanding of hazards, risks and early-warning systems are the main factors in building an effective community response. The level of understanding of these factors is proportional to a community's independence in taking appropriate actions and is reflected in the community's preparedness. Further, the experience gained from responses to the eruptions of Sinabung and Kelud as well as from other Indonesian volcanoes shows that a local leader plays an important role in mobilizing a community. Such a leader utilizes knowledge of local culture (local wisdom) and an understanding of the community's character to encourage community members to participate and empower themselves. Sharing of information and the formulation of a disaster mitigation plan by government and community also results in a shared commitment to increase participation and community empowerment. Indonesia uses four volcano alert levels to activate community mitigation plans. In order of increasing criticality, these are: Normal, Advisory (Waspada), Watch (Siaga) and Warning (Awas). Implementation of policies, strategies and actions are tied to each level. In addition, mandatory actions by certain stakeholders are required at certain alert levels. These mandatory actions increase the effectiveness of disaster mitigation by both government and community. A comparison of the crisis responses of the Sinabung and Kelud communities is carried out here in order to better understand problems, learn lessons and improve the process of Indonesia's community mitigation policies. Local culture, the start time of community involvement (in normal versus crisis time), degree and level of government involvement, roles of local leaders, local perception of hazards, and political intervention are all factors that influence a community's response at the time of eruption. We stress that an understanding of the “ community character ,” which includes a number of cultural, social and knowledge parameters is essential for effective crisis management. The responses of Sinabung and Kelud illustrate both the key role of community in disaster mitigation and the synergy that can result from close coordination and collaboration between government and community.

Journal of Volcanology and Geothermal Research↗
USGS2019

Using a process-based model of pre-eruptive seismic patterns to forecast evolving eruptive styles at Sinabung Volcano, Indonesia

Most volcanoes worldwide are not monitored in real-time; for those that are, patterns of pre-eruptive earthquakes coupled with conceptual models of magma ascent enable short-term forecasting of eruption onset. Basic event locations, characterization of background seismicity, and recognition of changes in earthquake types and energy release are most important to successful eruption forecasting. During renewed activity at Sinabung volcano, Indonesia, this approach was used by the Center for Volcanology and Geological Hazards Mitigation (CVGHM) and the USGS Volcano Disaster Assistance Program to forecast eruption onset, identify changes in eruptive styles and raise or lower alert levels and extend or contract evacuation zones. After > 400 years of quiescence, Sinabung began erupting in August 2010. The volcano was unmonitored at the onset of these eruptions, which were phreatic, but soon after a monitoring network was installed by CVGHM. Increasing swarms of high-frequency volcano tectonic (VT) earthquakes were used to forecast continuing phreatic eruptions. Volcanic activity decreased in mid-September 2010, while additional intrusions at depth (inferred from continued distal VT swarms) continued through September 2013, when explosive phreatic eruptions recurred. Explosive eruptions were forecast based on increases in real-time seismic amplitude measurement (RSAM) and VT seismicity. Seismicity changed markedly in late November and early December 2013 with the occurrence of deep earthquakes and an overall transition from low-frequency (LF) dominated and irregular (in time and magnitude) earthquakes to more regular LF and hybrid seismicity – a transition that accompanied the continued rise, eventual emergence and growth of a lava dome in the summit crater. This lava dome was first observed on 18 December. In late December 2013 to early January 2014, the eruptive style changed again as additional ascending magma deformed the summit and the dome grew beyond the capacity of the summit crater, resulting in the en masse collapse of the lava dome (2 Mm3) on 11 January and the largest pyroclastic flow to date. The collapse was forecast on the basis of a several order of magnitude increase in RSAM, continued strong distal VT seismicity, an increase in proximal seismicity, and large-scale observed deformation of the summit area. Similarly, a later collapse of a second summit lava dome on 1 February 2014 was forecast on the basis of increased distal seismicity. Here, we demonstrate how a process-based volcano seismicity model was used in combination with real-time data to forecast the time and magnitude of eruptions, as well as changes in eruption style.

Sinabung Volcano↗
USGS2019

Overview of the eruptions of Sinabung Volcano, 2010 and 2013-present, and details of the 2013 phreatomagmatic phase

A small phreatic eruption of Sinabung Volcano, North Sumatra on August 28, 2010, at 18:30 local time marked the first eruption in the past ~1200 years. The eruption took place from two small vents in the south crater area. Explosions and ash emissions from these vents generated multiple ash plumes that reached altitudes of up to 5 km during early- to mid-September. By the end of September, only low level steam plumes were visible and the alert level was reduced from Level 4 (highest) to Level 3. The 2010 eruption effectively ended at this time. Beginning two days after the initial 2010 eruption, activity of the eruption has been monitoring continuously by a telemetered seismic network surrounding the volcano and by remotely sensed observations. This monitoring system was supplemented with a near-field continuous GPS network, beginning in February 2011. Persistent fumarolic emissions continued for almost 3 years following the 2010 eruption, before a new eruption began on 15 September 2013. This eruption continues to the present. The ongoing eruption is divided into 5 major phases: 1) phreatomagmatic phase (July 2013 - 18 December 2013); 2) first dome and collapse phase with pyroclastic density currents (PDCs; block-and-ash flows and related surges) to south (18 December 2013 - 10 January 2014); 3) lava-flow and collapse phase (10 January 2014 - mid-September 2014); 4) second lava dome and collapse phase with PDCs to south (mid-September 2014 - July 2015); 5) lava dome collapse and ash explosion phase with PDCs to southeast and east (August 2015 - present). The volcano erupted intermittently during the early phreatomagmatic phase with small vertical ash explosions. Then the eruption became increasingly vigorous with more repetitive and intense vertical ash explosions during late October through November. The first small pyroclastic density currents (PDCs) began on November 1. These pyroclastic flows descended the southeastern flank to a distance of 2 km.

Journal of Volcanology and Geothermal Research↗
USGS2019

Joint 3-D tomographic imaging of Vp, Vs and Vp/Vs and hypocenter relocation at Sinabung volcano, Indonesia from November to December 2013

We conducted travel time tomography using P- and S-wave arrival times of volcanic-tectonic (VT) events that occurred between November and December 2013 to determine the three-dimensional (3D) seismic velocity structure (Vp, Vs, and Vp/Vs) beneath Sinabung volcano, Indonesia in order to delineate geological subsurface structure and to enhance our understanding of the volcanism itself. This was a time period when phreatic explosions became phreatomagmatic and then magma migrated to the surface forming a summit lava dome. We used 4846 VT events with 16,138 P- and 16,138 S-wave arrival time phases recorded by 6 stations for the tomographic inversion. The relocated VTs collapse into three clusters at depths from the surface to sea level, from 2 to 4 km below sea level, and from 5 to 8.5 km below sea level. The tomographic inversion results show three prominent regions of high Vp/Vs (~ 1.8) beneath Sinabung volcano at depths consistent with the relocated earthquake clusters. We interpret these anomalies as intrusives associated with previous eruptions and possibly surrounding the magma conduit, which we cannot resolve with this study. One anomalous region might contain partial melt, at sea level and below the eventual eruption site at the summit. Our results are important for the interpretation of a conceptual model of the “plumbing system” of this hazardous volcano.

Sinabung Volcano↗
USGS2019

Gas and ash emissions associated with the 2010–present activity of Sinabung Volcano, Indonesia

Sinabung Volcano (Sumatra, Indonesia) awoke from over 1200 years of dormancy with multiple phreatic explosions in 2010. After a period of quiescence, Sinabung activity resumed in 2013, producing frequent explosions, lava dome extrusion, and pyroclastic flows from dome collapses, becoming one of the world's most active volcanoes and displacing over 20,000 citizens. This study presents a compilation of the geochemical datasets collected by the Indonesian Center for Volcanology and Geological Hazard Mitigation (CVGHM) from 2010 - current (2016), which provides insights into the evolution of the eruption. Based on observations of SO2 emissions, ash componentry, leachate chemistry, and bulk ash geochemistry, the eruption can be split into five distinct phases. The initial stage of phreatic summit explosions occurred from August - October 2010, during which background SO2 emissions averaged ~550 ± 180 t/d (1 s.d.). An eruptive pause (phase two) starting in October 2010 abruptly ended in September 2013 with a resumption of conduit-clearing eruptions. This third phase had a relatively modest background SO2 emission rate (avg. ~410 ± 275 t/d) and produced ash consisting entirely of accidental ejecta with high S/Cl leachate ratios (up to 30), suggestive of deep-sourced magma and the incorporation of hydrothermal sulfur-bearing phases. The most intense phase of the eruption (phase four) occurred from December 2013 to February 2014, when juvenile magma first reached the surface. This period included dozens of large eruptions per day, high SO2 emission rates (average: 1,120 ± 1,030 t/d, peak: ~3,800 t/d), the onset of lava dome extrusion, and a dramatic drop in S/Cl ash leachates to ratios < 5, all reflecting increased degassing from shallow magma and the clearing out of sulfurous phases from the old hydrothermal system. From late February 2014 through the time of writing (September 2016), Sinabung settled into a relatively steady state of lower activity (phase five). Ash emissions now consist of dominantly juvenile material, and background SO2 emission rates have been progressively decreasing to an average of ~250 - 300 t/d. Starting August 2016, SO2 emissions started being measured in a continuous manner using a network of permanent scanning DOAS instruments. We find that long-term SO2 emission rates have been gradually declining at Sinabung since early 2014, consistent with an apparent decrease in magma supply. Our degassing model suggests that large explosions and pyroclastic flows could continue in the near-term owing to conduit plugging and dome collapses, remaining a major threat until the magma supply rate decreases further and the eruption ends.

Sumatra↗