Geology Reports⌕ Search

Geology topics

Research about Indonesia

Source-linked reports with geographic coverage including Indonesia.

At least 37 records · Page 2Linked to original sources

Construction of probabilistic event trees for eruption forecasting at Sinabung volcano, Indonesia 2013–14

Eruptions of Sinabung volcano , Indonesia have been ongoing since 2013. Since that time, the character of eruptions has changed, from phreatic to phreatomagmatic to magmatic explosive eruptions, and from production of a lava dome that collapsed to a subsequent thick lava flow that slowly ceased to be active, and later, to a new lava dome. As the eruption progressed, event trees were constructed to forecast eruptive behavior six times, with forecast windows that ranged from 2 weeks to 1 year: November 7–10, December 12–14, and December 27, 2013; and January 9–10, May 13, and October 7, 2014. These event trees were successful in helping to frame the forecast scenarios, to collate current monitoring information, and to document outstanding questions and unknowns. The highest probability forecasts closely matched outcomes of eruption size (including extrusion of the first dome), production of pyroclastic density currents, and pyroclastic density current runout distances. Events assigned low probabilities also occurred, including total collapse of the lava dome in January 2014 and production of a small blast pyroclastic density current in February 2014.

Sinabung volcano↗

Monitoring, forecasting collapse events, and mapping pyroclastic deposits at Sinabung volcano with satellite imagery

During the ongoing (2013–present) eruption of Sinabung volcano, north Sumatra, we have routinely used a variety of satellite remote sensing data to observe and forecast lava dome and lava flow collapse events, to map the resulting pyroclastic deposits, and to estimate effusion rates. In this paper, we focus on the first two years of the current eruption (September 2013–December 2015), and we summarize major events in 2016. We divide the eruption into 5 major phases: 1) phreatomagmatic (July 2013–18 December 2013), 2) first dome growth and collapse (18 December 2013–10 January 2014), 3) lava-flow (10 January 2014–mid-September 2014), 4) second lava dome and collapse (mid-September 2014–July 2015), 5) lava dome collapse and ash explosion phase (August 2015–present). Throughout the eruption, remotely sensed information has been instrumental in assessing the stability of the lava dome and flow and to forecast collapse events that produce pyroclastic density currents (PDCs: block-and-ash flows, co-ignimbrite surges, and blasts). Forecasts based on remote sensing data in combination with seismic, geodetic and gas-monitoring data have also helped inform decisions related to alert levels and evacuations. Relatively unusual aspects of the Sinabung eruption include the transition from dome to flow morphology (phase 2 to phase 3 transition) and the frequent occurrence during phase 3 of collapses from the lava flow-front and flow-margins—collapses that produced extensive pyroclastic density currents. By analogy to the well-known “Merapi type” collapses and pyroclastic deposits, we propose that lava flow-front and flow-margin collapses with associated PDCs be known as “Sinabung type.” Although detailed study of deposits has not been possible due to continuing hazards, our observations suggest that the transition from lava dome to lava flow and the occurrence of flow-front and flow-margin collapses reflect a particular combination of lava viscosity and steepness of slope. Our observations also show clear evidence of at least one slope-parallel high-velocity and dilute PDC (a “blast”) that emanated from a lava-margin collapse site 500 m downslope from the vent. This 1 February 2014 blast downed and singed a forest out to at least 3.9 km from the collapse site and killed 16 people. We also use a combination of field and remotely sensed data to map the distribution of Sinabung deposits. We estimate eruptive volumes and extrusion rates by combining sequential measurements of lava surface and pyroclastic flow areas with thickness estimates derived from simple geometric assumptions, oblique photographs and Digital Elevation Models (DEMs) derived from remotely sensed data. Our estimates of short-term effusion rates vary widely on a daily to weekly basis, from <1 to >20 m 3 s −1 . In a few cases, periods of increased extrusion precede lava flow-front collapses by a few days to a week, suggesting delays in transmittance of effusion pulses as lava moves from vent to flow front. We find that, as of 1 January 2016, the total area of deposits is 10 7 m 2 , and their approximate deposit volume is about 0.3 km 3 , equivalent to 0.2 km 3 Dense Rock Equivalent (DRE). We anticipate that our deposit maps will be valuable in the future as a framework for the study of the magmatic and textural evolution of eruptive products through time.

Sinabung volcano↗

Globally detected volcanic lightning and umbrella dynamics during the 2014 eruption of Kelud, Indonesia

Volcanic lightning shows considerable promise as a monitoring and research tool to characterize explosive eruptions. Its key strengths are rapid and remote detection, because the radio signals produced by lightning can propagate thousands of km at the speed of light. Despite these tantalizing properties, the scientific work on volcanic lightning has only recently started gaining momentum. Much more is needed to understand what lightning reveals about the evolution of an eruption in near-real time. Here we examine the timing and energy release of lightning generated by the eruption of Kelud volcano in Indonesia on 13 February 2014, as detected by the World Wide Lightning Location Network (WWLLN). The eruption column reached at least 26 km above sea level, representing the highest plume since the advent of global lightning networks in the last decade. Therefore, it provides valuable constraints on the electrification of end-member, sustained Plinian columns. We investigate the lightning in context with satellite images, photographs, and other published studies. Results show that the earliest satellite-detected activity was a thermal anomaly at ~15:46 UTC, corresponding to a directed blast at the onset of eruption (and only a few lightning strokes). Following a brief pause, the eruption produced a sustained column and umbrella cloud that spread outward into the tropical stratosphere. Rates of umbrella expansion provide an average mass eruption rate (MER) in the range of 8 × 10 7 –1 × 10 8 kg s −1 . A more nuanced picture emerges from the time-varying MERs (determined between each satellite pass), which show rapid intensification during the first hour of eruption, followed by constant MER for about an hour, and waning toward the end (after ~17:50 UTC). At this stage, decreasing flux into the umbrella cloud coincides with column instability and formation of pyroclastic density currents, as recorded by photos from the ground ~17:45 UTC. We infer that some of the erupted mass partitioned into ground-hugging currents, leading to a lower apparent MER. Interestingly, there is not a 1:1 correlation between lightning intensity and MER over the course of eruption. Stroke rates increase sharply within the first 30–40 min (during rapid intensification of the plume), and then drop below 2 strokes per min once the MER remains constant. This suggests that electrification was controlled by the rate of increase in MER—in other words, the acceleration of particles out of the vent. We also show that lightning stroke-rates and energies are greatest within 50 km of the vent, even when the ash cloud extends >200 km downwind, indicating that lightning was focused in the regions of highest particle concentration and turbulence. Overall, we conclude that abrupt changes in lightning rates are clearly linked to changes in eruption behavior, and that rapid detection could aid monitoring efforts to characterize eruption rates or styles.

Journal of Volcanology and Geothermal Research↗

3-D seismic tomographic study of Sinabung Volcano, Northern Sumatra, Indonesia, during the inter-eruptive period October 2010-July 2013

We estimated 3D P-wave velocity (Vp), S-wave velocity (Vs), and Vp/Vs tomographic structures in and around the Sinabung Volcano area, Northern Sumatra, Indonesia during the intereruptive period between October 2010 and July 2013 using high-quality P- and S-wave phases from > 600 local volcano tectonic (VT) earthquakes and more than 5000 phase picks. We used the iteritive damped-least-squares method of SIMULPS12 to simultaneously invert for the velocity structure and hypocenter relocations. We find that the majority of the VT earthquakes locate beneath Sinabung volcano to depths of about 8 km, with distal hypocenters dominantly to the north and northwest of the summit. We find six anomalous regions with the most notable including an area of low seismicity, high Vp/Vs and low Vp at depths of around 5 km slightly south of the summit. We believe this anomaly most likely represents hot material with some partial melt that may be the source for the 2010 phreatic eruptions and continued unrest throughout the study period. We also find a region below the summit with high seismicity, high Vp, high Vs and low Vp/Vs that we believe likely represents old intrusive material that has been fractured by the processes of the 2010 phreatic eruptions. Other shallow anomalous values of Vp, Vs and Vp/Vs are consistent with the known geology of Sinabung: anomalies of high Vp and Vs likely represent old, well-consolidated rock, and are observed to the SW and NE of Sinabung; while low Vp, low Vs and high Vp/Vs likely represent increased fracture densities, temperatures, or gas contents in geothermal features to the west of the summit. Finally we compare our results to those of the tomographic study during the October to November 2013 eruptive period in Nugraha et al. (this issue) and find evidence that suggests magma rose from the lower region imaged in our study to within a few kilometers of the surface.

Journal of Volcanology and Geothermal Research↗

La Niña-driven flooding in the Indo-Pacific warm pool during the past millennium

Extreme precipitation events are one of the most consequential components of climate change for society. The El Niño-Southern Oscillation (ENSO) is the dominant mode of precipitation variability in the tropics and causes severe flooding and drought in many socioeconomically vulnerable regions. It remains unclear how tropical rainfall extremes and ENSO are changing in response to anthropogenic forcing, demanding that we investigate the relationships between precipitation, ENSO, and external forcing in the past. Lake sediment records have provided benchmark records of extreme flood events from the eastern tropical Pacific, where paleofloods have been interpreted to reflect El Niño events during the last millennium. However, the connections between flooding and ENSO variability in this region are uncertain, and the eastern Pacific can only capture precipitation events driven by El Niño, not La Niña. Thus, it is unclear how the ENSO system and tropical rainfall extremes have changed in the recent past. Here, we reconstruct flood events during the past millennium using a lake sediment record from East Java, Indonesia, which can provide insight into flooding driven by La Niña. We detect flood frequency variations in the western tropical Pacific that are highly coherent with records from the eastern part of the basin over the past millennium. Our findings demonstrate that heavy rainfall and flooding occurs more frequently on both sides of the tropical Pacific during periods of warmer Northern Hemisphere mean temperatures, implying that ENSO-driven rainfall extremes could intensify in the near future.

Quaternary Science Reviews↗

Opal-A in glassy pumice, acid alteration, and the 1817 phreatomagmatic eruption at Kawah Ijen (Java), Indonesia

At Kawah Ijen (Indonesia), vigorous SO 2 and HCl degassing sustains a hyperacid lake (pH ~0) and intensely alters the subsurface, producing widespread residual silica and advanced argillic alteration products. In 1817, a VEI 2 phreatomagmatic eruption evacuated the lake, depositing a widespread layer of muddy ash fall, and sending lahars down river drainages. We discovered multiple types of opaline silica in juvenile low-silica dacite pumice and in particles within co-erupted laharic sediments. Most spectacular are opal-replaced phenocrysts of plagioclase and pyroxene adjacent to pristine matrix glass and melt inclusions. Opal-bearing pumice has been found at numerous sites, including where post-eruption infiltration of acid water is unlikely. Through detailed analyses of an initial sampling of 1817 eruption products, we find evidence for multiple origins of opaline materials in pumice and laharic sediments. Evidently, magma encountered acid-altered materials in the subsurface and triggered phreatomagmatic eruptions. Syn-eruptive incorporation of opal-alunite clasts, layered opal, and fragment-filled vesicles of opal and glass, all suggest magma-rock interactions in concert with vesiculation, followed by cooling within minutes. Our experiments at magmatic temperature confirm that the opaline materials would show noticeable degradation in time periods longer than a few tens of minutes. Some glassy laharic sedimentary grains are more andesitic than the main pumice type and may represent older volcanic materials that were altered beneath the lake bottom and were forcefully ejected during the 1817 eruption. A post-eruptive origin remains likely for most of the opal-replaced phenocrysts in pumice. Experiments at 25°C and 100°C reveal that when fresh pumice is bathed in Kawah Ijen hyperacid fluid for 6 weeks, plagioclase is replaced without altering either matrix glass or melt inclusions. Moreover, lack of evidence for high-temperature annealing of the opal suggests that post-eruption alteration of pumice is more likely than pre-eruption envelopment of euhedral opal-replaced phenocrysts in dacitic melt. At Ijen and elsewhere, the ascent of magma into hydrous acid-altered mineral assemblages (e.g., opal, kaolinite, alunite) could induce rapid dehydration of hydrous minerals and amorphous materials, generating considerable steam and contributing to magmatic-hydrothermal and phreatomagmatic explosions.

Kawah Ijen↗

Hypocenter relocation along the Sunda arc in Indonesia, using a 3D seismic velocity model

The tectonics of the Sunda arc region is characterized by the junction of the Eurasian and Indo‐Australian tectonic plates, causing complex dynamics to take place. High‐seismicity rates in the Indonesian region occur due to the interaction between these tectonic plates. The availability of a denser network of seismometers after the earthquakes of M w "> M w 9.1 in 2004 and M w "> M w 8.6 in 2005 supports various seismic studies, one of which regards the precise relocation of the hypocenters. In this study, hypocenter relocation was performed using a teleseismic double‐difference (DD) relocation method (teletomoDD) combining arrival times of P and S waves from stations at local, regional, and teleseismic distances. The catalog data were taken from the Agency of Meteorology, Climatology, and Geophysics (BMKG) of Indonesia, and the International Seismological Centre (ISC) for the time period of April 2009 to May 2015. The 3D seismic‐wave velocity model with a grid size 1 &#xB0; &#xD7; 1 &#xB0; "> 1°×1° was used in the travel‐time calculations. Relocation results show a reduction in travel‐time residuals compared with the initial locations. The relocation results better illuminate subducted slabs and active faults in the region such as the Mentawai back thrust and the outer rise in the subduction zone south of Java. Focal mechanisms from the Global Centroid Moment Tensor catalog are analyzed in conjunction with the relocation results, and our synthesis of the results provides further insight into seismogenesis in the region.

Sunda arc↗

High spatio-temporal resolution observations of crater-lake temperatures at Kawah Ijen volcano, East Java, Indonesia

The crater lake of Kawah Ijen volcano, East Java, Indonesia, has displayed large and rapid changes in temperature at point locations during periods of unrest, but measurement techniques employed to-date have not resolved how the lake&rsquo;s thermal regime has evolved over both space and time. We applied a novel approach for mapping and monitoring variations in crater-lake apparent surface (&ldquo;skin&rdquo;) temperatures at high spatial (~32 cm) and temporal (every two minutes) resolution at Kawah Ijen on 18 September 2014. We used a ground-based FLIR T650sc camera with digital and thermal infrared (TIR) sensors from the crater rim to collect (1) a set of visible imagery around the crater during the daytime and (2) a time series of co-located visible and TIR imagery at one location from pre-dawn to daytime. We processed daytime visible imagery with the Structure-from-Motion photogrammetric method to create a digital elevation model onto which the time series of TIR imagery was orthorectified and georeferenced. Lake apparent skin temperatures typically ranged from ~21 to 33oC. At two locations, apparent skin temperatures were ~ 4 and 7 oC less than in-situ lake temperature measurements at 1.5 and 5 m depth, respectively. These differences, as well as the large spatio-temporal variations observed in skin temperatures, were likely largely associated with atmospheric effects such as evaporative cooling of the lake surface and infrared absorption by water vapor and SO2. Calculations based on orthorectified TIR imagery thus yielded underestimates of volcanic heat fluxes into the lake, whereas volcanic heat fluxes estimated based on in-situ temperature measurements (68 to 111 MW) were likely more representative of Kawah Ijen in a quiescent state. The ground-based imaging technique should provide a valuable tool to continuously monitor crater-lake temperatures and contribute insight into the spatio-temporal evolution of these temperatures associated with volcanic activity.

Kawah Ijen volcano↗

New insights into the Kawah Ijen hydrothermal system from geophysical data

Volcanoes with crater lakes and/or extensive hydrothermal systems pose significant challenges with respect to monitoring and forecasting eruptions, but they also provide new opportunities to enhance our understanding of magmatic–hydrothermal processes. Their lakes and hydrothermal systems serve as reservoirs for magmatic heat and fluid emissions, filtering and delaying the surface expressions of magmatic unrest and eruption, yet they also enable sampling and monitoring of geochemical tracers. Here, we describe the outcomes of a highly focused international experimental campaign and workshop carried out at Kawah Ijen volcano, Indonesia, in September 2014, designed to answer fundamental questions about how to improve monitoring and eruption forecasting at wet volcanoes.

Kawah Ijen↗

Stress and mass changes at a “wet” volcano: Example during the 2011–2012 volcanic unrest at Kawah Ijen volcano (Indonesia)

Since 2010, Kawah Ijen volcano has been equipped with seismometers, and its extremely acid volcanic lake has been monitored using temperature and leveling sensors, providing unprecedented time resolution of multiparametric data for an acidic volcanic lake. The nature of stress and mass changes of the volcano is studied by combining seismic analyses and volcanic lake measurements that were made during the strongest unrest ever recorded by the seismic network at Kawah Ijen. The distal VT earthquake swarm that occurred in May 2011 was the precursor of volcanic unrest in October 2011 that caused an increase in shallow earthquakes. The proximal VT earthquakes opened pathways for fluids to ascend by increasing the permeability of the rock matrix. The following months were characterized by two periods of strong heat and mass discharge into the lake and by the initiation of monochromatic tremor (MT) activity when steam/gases interacted with shallow portions of the aquifer. Significant seismic velocity variations, concurrent with water level rises in which water contained a large amount of steam/gas, were associated with the crises, that caused an although the unrest did not affect the shallow hydrothermal system at a large scale. Whereas shallow VT earthquakes likely reflect a magmatic intrusion, MT and relative seismic velocity changes are clearly associated with shallow hydrothermal processes. These results will facilitate the forecast of future crises.

Kawah Ijen↗

Kawah Ijen volcanic activity: A review

Kawah Ijen is a composite volcano located at the easternmost part of Java island in Indonesia and hosts the largest natural acidic lake in the world. We have gathered all available historical reports on Kawah Ijen’s activity since 1770 with the purpose of reviewing the temporal evolution of its activity. Most of these observations and studies have been conducted from a geochemical perspective and in punctuated scientific campaigns. Starting in 1991, the seismic activity and a set of volcanic lake parameters began to be weekly available. We present a database of those measurements that, combined with historical reports, allow us to review each eruption/unrest that occurred during the last two centuries. As of 2010, the volcanic activity is monitored by a new multi-disciplinary network, including digital seismic stations, and lake level and temperature measurements. This detailed monitoring provides an opportunity for better classifying seismic events and forecasting volcanic unrest at Kawah Ijen, but only with the understanding of the characteristics of this volcanic system gained from the historical review presented here.

Kawah Ijen↗

Comparing spatial capture–recapture modeling and nest count methods to estimate orangutan densities in the Wehea Forest, East Kalimantan, Indonesia

Accurate information on the density and abundance of animal populations is essential for understanding species' ecology and for conservation planning, but is difficult to obtain. The endangered orangutan ( Pongo spp.) is an example; due to its elusive behavior and low densities, researchers have relied on methods that convert nest counts to orangutan densities and require substantial effort for reliable results. Camera trapping and spatial capture&ndash;recapture (SCR) models could provide an alternative but have not been used for primates. We compared density estimates calculated using the two methods for orangutans in the Wehea Forest, East Kalimantan, Indonesia. Camera trapping/SCR modeling produced a density estimate of 0.16 &plusmn; 0.09&ndash;0.29 indiv/km 2 , and nest counts produced a density estimate of 1.05 &plusmn; 0.18&ndash;6.01 indiv/km 2 . The large confidence interval of the nest count estimate is probably due to high variance in nest encounter rates, indicating the need for larger sample size and the substantial effort required to produce reliable results using this method. The SCR estimate produced a very low density estimate and had a narrower but still fairly wide confidence interval. This was likely due to unmodeled heterogeneity and small sample size, specifically a low number of individual captures and recaptures. We propose methodological fixes that could address these issues and improve precision. A comparison of the overall costs and benefits of the two methods suggests that camera trapping/SCR modeling can potentially be a useful tool for assessing the densities of orangutans and other elusive primates, and warrant further investigation to determine broad applicability and methodological adjustments needed.

East Kalimantan↗

Coral 13C/12C records of vertical seafloor displacement during megathrust earthquakes west of Sumatra

The recent surge of megathrust earthquakes and tsunami disasters has highlighted the need for a comprehensive understanding of earthquake cycles along convergent plate boundaries. Space geodesy has been used to document recent crustal deformation patterns with unprecedented precision, however the production of long paleogeodetic records of vertical seafloor motion is still a major challenge. Here we show that carbon isotope ratios ( ) in the skeletons of massive Porites corals from west Sumatra record abrupt changes in light exposure resulting from coseismic seafloor displacements. Validation of the method is based on the coral response to uplift (and subsidence) produced by the March 2005 M w 8.6 Nias&ndash;Simeulue earthquake, and uplift further south around Sipora Island during a M&sim;8.4 megathrust earthquake in February 1797. At Nias, the average step-change in coral was 0.6&plusmn;0.1&permil;/m for coseismic displacements of +1.8 m and &minus;0.4 m in 2005. At Sipora, a distinct change in Porites microatoll growth morphology marks coseismic uplift of 0.7 m in 1797. In this shallow water setting, with a steep light attenuation gradient, the step-change in microatoll is 2.3&permil;/m , nearly four times greater than for the Nias Porites . Considering the natural variability in coral skeletal , we show that the lower detection limit of the method is around 0.2 m of vertical seafloor motion. Analysis of vertical displacement for well-documented earthquakes suggests this sensitivity equates to shallow events exceeding M w &sim;7.2 in central megathrust and back-arc thrust fault settings. Our findings indicate that the coral paleogeodesy technique could be applied to convergent tectonic margins throughout the tropical western Pacific and eastern Indian oceans, which host prolific coral reefs, and some of the world's greatest earthquake catastrophes. While our focus here is the link between coral , light exposure and coseismic crustal deformation, the same principles could be used to characterize interseismic strain during earthquake cycles over the last several millennia.

Earth and Planetary Science Letters↗