Geochemistry and Structure of Krakatoa Volcano in the Sunda Strait, Indonesia

geosciences

Article Geochemistry and Structure of Krakatoa Volcano in the Sunda Strait, Indonesia

Mirzam Abdurrachman 1,* ID , Sri Widiyantoro 2, Bambang Priadi 1 and Tauﬁk Ismail 3 1 Department of Geology, Bandung Institute of Technology, Bandung 40132, Indonesia; [email protected] 2 Global Geophysics Research Group, Bandung Institute of Technology, Bandung 40132, Indonesia; [email protected] 3 Department of Geology, Mineral Technology College of Indonesia, Bandung 40263, Indonesia; tauﬁ[email protected] * Correspondence: [email protected]; Tel.: +62-877-2237-9635

Received: 12 February 2018; Accepted: 23 March 2018; Published: 26 March 2018 Abstract: The violent eruption of Krakatoa Volcano located in the Sunda Strait, Indonesia, in 1883 represents one of the deadliest eruptions in human civilization. Although lots of data have been reported, the trajectory of the subducted slab and the upper mantle structure beneath this volcano are still rather poorly known. We combined geochemical data, major, trace and rare earth elements with seismic tomograms to characterize the deep structure of Krakatoa Volcano at the junction of Sumatra and Java subduction systems. Geochemical data are in agreement with the partial melting of mantle wedge in these subduction systems, based on previous studies, and this conclusion is also supported by inferences from P-wave tomographic model. Whereas, the tomographic image of S-wave suggests that subducted slab has been intruded by hot material of mantle upwelling. The presence of both partial melting of mantle wedge and mantle upwelling in the upper mantle might be caused by the thinning of subducted slab beneath Krakatoa Volcano.

Keywords: geochemistry; seismic tomogram; Krakatoa Volcano; Indonesia

1. Introduction The nature of magma’s genesis in subduction zone settings can be determined by whole rock geochemistry and isotopic compositions, even though whole-rock geochemistry and isotopic compositions are not always indicative of actual magma source characteristics, as they have commonly undergone through complicated differentiation processes such as fractional crystallization, magma mixing and also crustal assimilation. Furthermore, the geochemical variations of the mantle wedge, subducted oceanic crust, subducted sediments or melts and slab conﬁguration also play an important role in magma genesis [1–5]. Previous geological studies on Javan volcanoes including Krakatoa Volcano have mainly focused on magma genesis [6–10]. These studies suggest that magma beneath Krakatoa was originally derived from partial melting of the mantle wedge in the subduction system, with minor contamination from upper crustal and Java sediment melts. Here, we focus on whole-rock geochemistry and Sr isotopic data from 1883, 1981 and 2002 eruption products and combine them with seismic tomograms to determine magma source characteristics. We note that regional seismic tomographic studies of the mantle structure include the deep subduction beneath Krakatoa Volcano [11–14]. This study will therefore explore the detailed slab structure, in particular, in the upper mantle beneath Krakatoa. Along with the whole-rock geochemistry and Sr isotopic data, we present P- and S-wave seismic images produced by means of improved tomographic imaging model technique [11] in order to propose a new magma genesis model beneath the volcano.

Geosciences 2018, 8, 111; doi:10.3390/geosciences8040111 www.mdpi.com/journal/geosciences Geosciences 2018, 8, 111 2 of 10

2. Regional Geology Indonesia is situated at the boundaries of three major plates on the southern margin of the Eurasian Plate, those are: the Eurasian, Indo-Australian and Pacific-Philippine Plates [15] (Figure1). The Indo-Australian Plate moves northward, and is being subducted beneath the Eurasian Plate. The age of the subducted plate varies laterally. It is relatively young beneath Sumatra, i.e., around 40 Ma, where the subduction is highly oblique. In contrast, below Krakatoa, the age of the subducted plate is around 80 Ma, and it becomes significantly older, reaching up to 130 Ma in Eastern Java [16,17], where the convergence direction is almost perpendicular. The lateral variation of nature and age of the subducted plate influences the style of deformation and subducted angle along Sumatra to Java [18], which may produce a thinning in the subducted slab between Sumatra and Java [11], as mirrored by the abrupt change in the seismicity in the Sunda Strait [11]. The convergence rate increases from Sumatra to the East [19], i.e., about 6.7 to 7 cm/year [20,21]. The subduction zone is well defined by the volcanic arc, which stretches about 6000 km from Sumatra to the Molucca Sea, and by deep seismicity that extends down to a depth of 600 km [13,14,21–23].

Figure 1. Tectonic map of Indonesia showing the interaction of 3 major plates, modiﬁed from [15]. Red triangle depicts Krakatoa Volcano. Contour of subducted slab is taken from [12,15].

Krakatoa is an active volcano located at the junction of Sumatra and the Java Islands. This volcano marks the changing of the subduction system from oblique along Sumatra to perpendicular along Java (Figure1). The Krakatoa Volcanic complex is situated approximately 141 km from the Java Trench in the southwest of the volcano, whereas the Benioff zone is at a depth of approximately 120 km beneath the volcano [12,14].

3. Geology of Krakatoa The Krakatoa volcanic complex consists of Rakata, Sertung, Panjang and Anak Krakatoa Islands (Figure2). The ﬁrst three islands are remnants of the caldera collapse of the supposedly ancient Geosciences 2018, 8, 111 3 of 10

Krakatoa, and Rakata Island itself is a volcano that grew up together with Danan and Parbuwatan volcanoes before the cataclysmic eruption of 1883. Following [24], we divided the stratigraphy of Krakatoa into 4 eruption products (Figure2), and a new detailed stratigraphy for Anak Krakatoa is shown in Figure3, from old to young as follows:

1. Pre-historic Krakatoa; these products were erupted around 416 AD and consist of 2 dacitic lava flows, 1 pyroclastic flow and 1 pyroclastic fall. 2. Young Krakatoa; these products were erupted around 1200 AD and consist of 3 eruption centers: Rakata, Danan and Parbuwatan. During this period, the products were dominated by lava. 3. 1883; period of the destruction of the Rakata, Danan and Parbuwatan volcanoes, characterized by the 1883 caldera formation, and produced typical eruptive products. This rock unit was distributed on all three islands—Rakata, Panjang and Sertung—and was composed of pumice pyroclastic flows, minor pyroclastic falls and surge deposits. 4. Anak Krakatoa; since its birth in 1927 to 2017, Anak Krakatoa has erupted at least thirty-two times, showing a combination of explosive and effusive activities. Anak Krakatoa, with an elevation of about 300 m above sea level, is a volcano island, lying in the center of the Krakatoa volcanic complex. It is constructed of alternating layers of 18 lava flows and 18 pyroclastic deposits that have been built since 1927.

Figure 2. Geological map of the Krakatoa complex. See Figure3 for explanation of colors. Geosciences 2018, 8, 111 4 of 10

Figure 3. Stratigraphy of the Krakatoa complex. Sample positions indicated by red star for 1883, blue star for 1981 and green star for 2002. Samples are taken from [6,7].

4. Data and Method In this study, we used representative whole-rock geochemical data (major, trace and rare earth elements, as well as the Sr isotope) from the 1883, 1981 and 2002 eruption products taken from [6,7]; see Table1 and Figure3 for stratigraphic position of the samples. The geochemical data analysis involved basaltic andesite (AK-1) and andesite (DR12653) from Anak Krakatoa, with the other one (L2835) being andesite from the 1883 catastrophic eruption. The samples were processed using XRF and INAA for major, trace and rare earth elements. The Sr isotope was measured dynamically on a single collector Finningan MAT 261. We then combined our geochemical data interpretation results with P- and S-wave tomographic models to enrich our interpretation. We used travel-time residuals taken from an updated global data set by [25]. This data set covers the period 1964–1999 and includes nearly 100,000 global events recorded at a subset of about 7000 seismographic stations worldwide. We employed the same P- and S-ray data coverage to minimize bias or differences in the resulting models that may come from Geosciences 2018, 8, 111 5 of 10 different ray paths. The tomographic models were derived using a non-linear approach, which was implemented through an iterative technique to update 3-D model by employing 3-D ray tracing [11]. As in [11], here we employed a cellular representation of mantle structure by discretizing the entire mantle by a uniform grid of constant velocity cells of 5◦ by 5◦, but in the study area, we used a ﬁner grid of 0.5◦ by 0.5◦, allowing the resolution of relatively small-scale structures.

Table 1. Major, trace and rare earth elements, as well as Sr isotopes, from 1883, 1981 and 2002 AD eruption products. Data source 2002 from [6], 1883 and 1981 from [7]. Note: n.d means not determined.

Sampel L2835 DR12653 AK-1 Age 1883 AD 1981 AD 2002 AD SiO2 (wt %) 57.86 61.45 55.71 TiO2 1.18 1.02 1.06 87Sr/86Sr 0.704393 0.704428 0.704417 Rb (ppm) 59 40 23 Ba 321 275 181 Th 7453 5380 3034 Nb 6 5 4 K 1.73 1.21 0.78 La 21.2 19.9 15.02 Ce 49.4 45.1 35.4 Sr 296 324 383 Nd 27.7 25.7 21.5 Zr 185 184 123 Sm 6.48 6.35 5.29 Ti 0.71 0.61 0.65 Tb 1.09 1.1 0.86 Y 42 44 33 Yb 4.06 4.18 3.06 Lu 0.62 0.62 0.49 Sc 24 21 n.d

5. Results and Discussion

The SiO2 wt % contents for eruption products are 57.68% (1883 AD), 61.45% (1981 AD) and 55.71% (2002 AD), belonging to the medium- to high-K series. In addition, the eruption products are characterized by low TiO2 wt % contents (1.18% for 1883 AD, 1.02% for 1981 AD and 1.02% for 2002 AD). According to [26], volcanic rocks with TiO2 wt % contents <1.4% are typically magma related to the partial melting of the mantle wedge induced by the traditional fluids in the subduction system. The 87Sr/86Sr isotope ratios of Krakatoa coverage in the narrow range are from 0.704393 (1883 AD) to 0.704428 (1981 AD). This range is quite low compared to Javan volcanic rocks, but it is relatively high compared to those of other volcanoes of the region of similar depth to the Benioff zone (e.g., Cikuray, Galunggung, Guntur) [8]. Therefore, the 87Sr/86Sr isotope ratios of Krakatoa are unlikely to reflect pure mantle, indicating a low degree of crustal assimilation [6]. Table1 shows major, trace and rare earth elements, and Sr isotope from the 1883, 1981 and 2002 eruption products. These data were plotted in a La/Yb-Sc/Ni diagram, as proposed by [27], see Figure4, and into a TiO 2-MnO × 10-P2O5 × 10 diagram, as proposed by [28], see Figure5. As shown in Figures4 and5, these volcanic rocks are clearly arc magmas related to the partial melting of the mantle wedge in the subduction zone. The normalized diagram (spidergram) of the primitive mantle (Figure6) shows that eruptive products of all samples are similar to Honshu medium-K basalt, showing highly incompatible elements (left-hand side), which are more enriched than moderately incompatible elements. Moreover, it is markedly depleted in some high field strength elements, notably niobium (Nb). These features are common to many supra subduction zone magmas in mature island arcs and convergent margins. They are believed to reflect intra-crustal magma fractionation Geosciences 2018, 8, 111 6 of 10 superimposed on the signature of melting mantle wedge peridotite metasomatized by slab fluids [29] or residual titenate phase, most likely rutile in the slab during dehydration [30,31].

Figure 4. La/Yb vs. Sc/Ni diagram modiﬁed after [27].

Figure 5. TiO2 -MnO × 10-P2O5 × 10 diagram modiﬁed after [28]. Geosciences 2018, 8, 111 7 of 10

Figure 6. Normalized diagram to primitive mantle showing typical volcanic rocks from supra subduction zone. Data source: Primitive mantle [32], Honshu [33].

Figure 7. (a) Vertical section of P-wave model across Krakatoa Volcano. Tomogram is plotted as velocity perturbations relative to the ak135 velocity model [34]. The thin line in the cross section indicates the 410-km mantle discontinuity and the arrow depicts hot mantle upwelling. (b) Same as (a), but based on an S-wave model. (c) Index map, showing location of Krakatoa Volcano and A-B cross section line. Note the thinning in the subducted slab in the upper mantle below Krakatoa, which may be due to: (i) intrusion of hot mantle upwelling, (ii) strong slab pull in the study region [11,13], and (iii) the change from oblique subduction in Sumatra to frontal subduction in Java. Geosciences 2018, 8, 111 8 of 10

The P-wave tomographic model shows evidence of cold material, depicted by a high-velocity anomaly beneath Krakatoa Volcano (Figure7a). Figure7b shows the same cross-section, but through an S-wave tomographic model. These images clearly show a low-velocity anomaly related to hot materials of mantle intruding into the cold subducted slab. The TiO2 wt % contents, La/Yb-Sc/Ni, TiO2-MnO × 10-P2O5 × 10 and spidergram show a good agreement with the P-wave tomographic model, where the P-wave data clearly show cold material subducted beneath Krakatoa Volcano (Figure7a). This condition implies a common magma genesis in the subduction system. On the other hand, the S-wave tomographic model shows more clearly the hot material of the mantle upwelling that intruded into the subducted slab beneath this volcano (Figure7b) implying a magma genesis related to upper mantle plume.

6. Conclusions Whole-rock geochemistry data in conjunction with Sr isotopic data have provided new insight into magma genesis beneath Krakatoa Volcano. We found that the data from eruptive products of 1883, 1981 and 2002 AD represent volcanic rocks derived from partial melting of the mantle wedge in the subduction system, as shown by the P-wave tomogram. The presence of partial melting of the mantle wedge and mantle upwelling beneath Krakatoa Volcano may be controlled by the thinning in the subducted slab as depicted clearly by the S-wave tomogram shown in Figure7b, owing to strong slab pull in the eastern part of the Sunda Arc [11,13,18] and which may also be related to mantle upwelling beneath Krakatoa and the changing of the subduction style between Sumatra and Java.

Acknowledgments: We are thankful to the head and staff of BKSDA Bengkulu, Lampung Province, Sumatra for their permission, assistance and help during our field work. Author Contributions: Mirzam Abdurrachman, Sri Widiyantoro, Bambang Priadi, and Taufik Ismail analyzed the data and wrote the paper. Conflicts of Interest: The authors declare no conflict of interest.

References

1. Davidson, J.P. Crustal contamination versus subduction zone enrichment: Examples from the Lesser Antilles and implications for mantle source compositions of island arc volcanic rocks. Geochim. Cosmochim. Acta 1987, 51, 2185–2198. [CrossRef] 2. Hildreth, W.; Moorbath, S. Crustal contributions to arc magmatism in the Andes of Central Chile. Contributions Mineral. Petrol. 1988, 98, 455–489. [CrossRef] 3. McCulloch, M.T.; Gamble J., A. Geochemical and geodynamical constraints on subduction zone magmatism. Earth Planet. Sci. Lett. 1991, 102, 358–374. [CrossRef] 4. Hawkesworth, C.J.; Hergt, J.M.; Ellam, R.M.; McDermott, F. Element ﬂuxes associated with subduction related magmatism. Philos. Trans. R. Soc. Lond. A 1991, 335, 393–405. [CrossRef] 5. Pearce, J.A.; Peate, D.W. Tectonic implications of the composition of volcanic arc magmas. Annu. Rev. Earth Planet. Sci. 1995, 24, 251–285. [CrossRef] 6. Gardner, M.F.; Troll, V.R.; Gamble, J.A.; Gertisser, R.; Hart, G.L.; Ellam, R.M.; Harris, C.; Wolff, J.A. Crustal differentiation processes at Krakatau Volcano, Indonesia. J. Petrol. 2013, 54, 149–182. [CrossRef] 7. Turner, S.; Foden, J.U. Th and Ra disequilibria, Sr, Nd, and Pb isotope and trace element variations in Sunda arc lavas: Predominance of a subducted sediment component. Contributions Mineral. Petrol. 2001, 142, 43–57. [CrossRef] 8. Handley, H.K.; Blichert-Toft, J.; Gertisser, R.; Macpherson, C.G.; Turner, S.P.; Zaennudin, A.; Abdurrachman, M. Insights from Pb and O isotopes into along-arc variations in subduction inputs and crustal assimilation for volcanic rocks in Java, Sunda arc, Indonesia. Geochim. Cosmochim. Acta 2014, 139, 205–226. [CrossRef] 9. Abdurrachman, M. Geology and Petrology of Quaternary Papandayan Volcano and Genetic Relationship of Volcanic Rocks from the Triangular Volcanic Complex around Bandung Basin, West Java, Indonesia. Unpublished Ph.D. Thesis, Akita University, Akita, Japan, 2012. Geosciences 2018, 8, 111 9 of 10

10. Abdurrachman, M.; Yamamoto, M. Geochemical variation of Quaternary volcanic rocks in Papandayan area, West Java, Indonesia: A role of crustal component. J. Geol. Soc. Thailand 2013, 40–57. 11. Widiyantoro, S.; Pesicek, J.D.; Thurber, C.H. Subducting slab structure below the eastern Sunda arc inferred from non-linear seismic tomographic imaging. Geol. Soc. Lond. Spec. Publ. 2011, 355, 139–155. [CrossRef] 12. Abdurrachman, M.; Widiyantoro, S.; Priadi, B.; Alim, M.Z.A.; Dewangga, A.H. Proposed new wadati— benioff zone model in Java-Sumatra subduction zone and its tectonic implication. In Proceedings of the Joint Convention Balikpapan 2015 HAGI-IAGI-IAFMI-IATMI, Balikpapan, Indonesia, 5–8 October 2015. 13. Widiyantoro, S.; van der Hilst, R.D. Mantle structure beneath Indonesia inferred from high-resolution tomographic imaging. Geophys. J. Int. 1997, 130, 167–182. [CrossRef] 14. Widiyantoro, S.; van der Hilst, R.D. Structure and evolution of lithospheric slab beneath the Sunda arc, Indonesia. Science 1996, 271, 1566–1570. [CrossRef] 15. Hamilton, W.B. Tectonics of the Indonesian Region. Available online: https://pubs.usgs.gov/pp/1078/ report.pdf (accessed on 10 January 2018). 16. Mueller, R.D.; Roest, W.R.; Royerj, Y.; Gahagan, L.M.; Sclater, J.G. Digital isochrones of the world’s ocean floor. J. Geophys. Res. 1997, 102, 3211–3214. [CrossRef] 17. Syracuse, E.M.; Abers, G.A. Global compilation of variation in slab depth beneath arc volcanoes and implications. Geochem. Geophys. Geasyst. 2006, 7, Q05017. [CrossRef] 18. Cloetingh, S.; Wortel, R. Stress in the Indo-Australian plate. Tectonophysics 1986, 132, 49–67. [CrossRef] 19. Minster, J.B.; Jordan, T.H. Present-day plate motions. J. Geophys. Res. 1978, 83, 5331–5354. [CrossRef] 20. Tregoning, P.; Brunner, F.K.; Bock, Y.; Puntodewo, S.S.O.; McCaffrey, R.; Genrich, J.F.; Calais, E.; Rais, J.; Subarya, C. First geodetic measurement of convergence across the Java trench. Geophys. Res. Lett. 1994, 21, 2135–2138. [CrossRef] 21. Hall, R. Indonesia, Geology. In Encyclopedia of Islands; Gillespie, R., Clague, D.A., Eds.; University of California Press: Oakland, CA, USA, 2008. 22. Ghose, R.; Yoshioka, S.; Oike, K. Three dimensional numerical simulations of the subduction dynamics in the Sunda arc region, Southeast Asia. Tectonophysics 1990, 181, 223–255. [CrossRef] 23. Puspito, N.T.; Shimazaki, K. Mantle structure and seismotectonics of the Sunda and Banda arcs. Tectonophysics 1995, 251, 215–228. [CrossRef] 24. Stehn, C.E. The geology and volcanism of the Krakatau group. In Guidebook for the 4th Pacific Science Congress; University of Michigan: Ann Arbor, MI, USA, 1929; pp. 1–55. 25. Engdahl, E.R.; van der Hilst, R.D.; Buland, R. Global teleseismic earthquake relocation with improved travel times and procedures for depth determination. Bull. Seism. Soc. Am. 1998, 88, 722–743. 26. Metcalf, R.V.; Shervais, J.W. Suprasubduction-Zone Ophiolites: Is There Really an Ophiolite Conundrum? In Special Paper 438: Ophiolites, Arcs, and Batholiths: A Tribute to Cliff Hopson; Geological Society of America: Boulder, CO, USA, 2008; Volume 438, pp. 191–222, ISBN 978-0-8137-2438-6. 27. Bailey, J.C. Geochemical criteria for a refined tectonic discriminant of orogenic andesites. Chem. Geol. 1981, 32, 139–154. [CrossRef]

28. Mullen, E.D. MnO/TiO2/P2O5: A minor element discriminant for basaltic rocks of oceanic environments and its implication for petrogenesis. Earth Planet. Sci. Lett. 1983, 62, 53–62. [CrossRef] 29. Gill, R. Igneous Rocks and Processes: A Practical Guide; John Wiley & Sons, Wiley-Blackwell: Hoboken, NJ, USA, 2010; Volume 57, ISBN 978-0-632-06377-2. 30. Brenan, J.M.; Shaw, H.F.; Ryerson, F.J.; Phinney, D.L. Mineral-aqueous fluid partitioning of trace elements at 900 ◦C and 2.0 GPa: Constraints on the trace element geochemistry of mantle and deep crustal fluids. Geochim. Cosmochim. Acta 1995, 59, 3331–3350. [CrossRef] 31. Stalder, R.; Foley, S.F.; Brey, G.P.; Horn, I. Mineral-aqueous fluid partitioning of trace elements at 900–1200 ◦C and 3.0–5.7 GPa: New experimental data for garnet, clinopyroxene, and rutile, and implications for mantle metasomatism. Geochim. Cosmochim. Acta 1998, 62, 1781–1801. [CrossRef] 32. Sun, S.S.; McDonough, W.F. Chemical and isotopic systematics of oceanic basalts: Implication for mantle composition and processes. In Magmatism in Ocean Basins; Saunders, A.D., Norry, M.J., Eds.; Geological Society of London Special Publication: London, UK, 1989; Volume 42, pp. 313–345. Geosciences 2018, 8, 111 10 of 10

33. Gust, D.A.; Arculus, R.J.; Kersting, A.B. Aspect of magma sources and processes in the Honshu arc. Can. Mineral. 1997, 35, 347–365. 34. Kennett, B.L.N.; Engdahl, E.R.; Buland, R. Constraints on seismic velocities in the Earth from travel times. Geophys. J. Int. 1995, 122, 108–124. [CrossRef]

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).