Reconciling the Climate and Ozone Response to the 1257 CE Mount Samalas Eruption

Total Page:16

File Type:pdf, Size:1020Kb

Reconciling the Climate and Ozone Response to the 1257 CE Mount Samalas Eruption Reconciling the climate and ozone response to the 1257 CE Mount Samalas eruption David C. Wadea, Céline M. Vidalb, N. Luke Abrahama,c, Sandip Dhomsed, Paul T. Griffithsa,c, James Keeblea,c, Graham Mannd, Lauren Marshalla, Anja Schmidta,b, and Alexander T. Archibalda,c,1 aCentre for Atmospheric Science, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK; bDepartment of Geography, Downing Place, Cambridge, CB2 3EN, UK; cNational Centre for Atmospheric Science, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK; dSchool of Earth and Environment, University of Leeds, Leeds, LS2 9JT, UK This manuscript was compiled on September 1, 2020 1 The 1257 CE eruption of Mount Samalas (Indonesia) is the source of ies have investigated the climate response to the 1257 Mt. 21 2 the largest stratospheric injection of volcanic gases in the Common Samalas eruption. However, when compared with observa- 22 3 Era. Sulfur dioxide emissions produced sulfate aerosols that cooled tions or proxy-records (11–13), most previous simulations 23 4 Earth’s climate with a range of impacts on society. The co-emission overestimate the surface cooling response to Mt. Samalas, and 24 5 of halogenated species has also been speculated to have led to wide- other large eruptions (14). Despite the 1257 eruption of Mt. 25 6 scale ozone depletion. Samalas emitting 1.8× more SO2 than the 1815 eruption of 26 7 Here we present simulations from HadGEM3-ES, the first fully- Mt. Tambora, the Northern Hemisphere temperature response, 27 8 coupled Earth system model, with interactive atmospheric chemistry reconstructed by tree-ring maximum latewood density (MXD) 28 9 and a microphysical treatment of sulfate aerosol, used to assess the methods, is relatively similar (7, 8). This has been a puzzle 29 10 chemical and climate impacts from the injection of sulfur and halo- for the climate modeling community. Timmreck et al. (15) 30 11 gen species into the stratosphere as a result of the Mt. Samalas and Stoffel et al. (8) were able to simulate a broad agreement 31 12 eruption. Whilst our model simulations support a surface air temper- between their Atmosphere Ocean General Circulation Model 32 13 ature response to the eruption of the order of -1°C, performing well (AO-GCM) simulations and tree-ring based reconstructions, 33 14 against multiple reconstructions of surface temperature from tree- but Stoffel et al. have shown that the surface temperature 34 15 ring records, we find little evidence to support significant injections response was very sensitive to the assumed eruption month 35 16 of halogens into the stratosphere. and the height of injection (8), whilst Timmreck et al. found 36 17 Including modest fractions of the halogen emissions reported from that the size of the prescribed aerosol effective radii were very 37 18 Mt. Samalas leads to significant impacts on the composition of the important (15). 38 19 atmosphere and on surface temperature. As little as 20% of the halo- Both the changes in aerosol number concentrations and 39 20 gen inventory from Mt. Samalas reaching the stratosphere would re- size and gas-phase composition (chemistry) are important 40 21 sult in catastrophic ozone depletion; extending the surface cooling aspects of the role of large injections of sulfur (like those 41 22 caused by the eruption. However, based on available proxy-records from Mt. Samalas) on climate (16, 17). Changes in aerosol 42 23 of surface temperature changes, our model results support only very properties and changes in chemistry are strongly coupled. 43 24 minor fractions (1%) of the halogen inventory reaching the strato- Large eruptions, like Mt. Samalas, can inject their volatile 44 25 sphere and suggest that further constraints are needed to fully re- gases into the stratosphere. In this region, volcanic sulfate 45 26 solve the issue. Samalas | Climate | Ozone | Modeling Volcanic Impacts Significance Statement 1 The eruption of Mount Samalas in 1257 (Indonesia),DRAFT with The Mount Samalas eruption in 1257, one of the largest ex- 2 Volcanic Explosivity Index (VEI) 7, was identified by Lavigne plosive volcanic eruptions in the Common Era, has proven 3 et al. (1) as the source of the largest sulfate spike of the a complex case for climate models – which have generally 4 Common Era in both Greenland and Antarctica ice cores overestimated the climate response compared with proxy data. 5 (2, 3). Analyses of erupted products and their melt inclusions Here we perform the first Earth system model simulations of 6 suggests that the eruption released 158±12 megatons (Tg) of the impacts of the Mount Samalas eruption using a range of 7 sulfur dioxide (SO2), together with 227±18 Tg of chlorine (Cl) SO2 and halogen emission scenarios. Reported halogen emis- 8 and up to 1.3 Tg of bromine (Br) (4); making it responsible sions are considerable from the eruption but using our model 9 for the emissions of the greatest quantity of volcanic gases into simulations and reconstructed climate response we can rule 10 the stratosphere in the Common Era (4–6). out all but minor halogen emissions reaching the stratosphere. 11 Tree-ring proxy records and historical narratives point to- Including a minor fraction of the halogen inventory reaching the 12 ward a strong surface cooling in the Northern Hemisphere stratosphere captures the observed "muted" climate response 13 following the eruption, in the summer of 1258 (1, 7, 8). Con- but results in significant ozone depletion with implications for 14 temporary narrative sources suggest an unusually cloudy, rainy UV exposure and human health. 15 and cold summer in Europe in 1258 (1, 7). The climatic impact 16 of the eruption is suspected to have amplified famines and D.C.W., C.V. and A.T.A. conceived the study and A.S. contributed to the data interpretation. D.C.W., N.L.A., P.T.G., J.K., L.M., G.M. and S.D. performed model development. D.C.W. performed model 17 political turmoil in Europe and Japan, with the most severe and proxy-data analysis. D.C.W., C.V., A.T.A. and A.S. wrote the manuscript. 18 socio-economic consequences reported in England (7,9, 10). The authors declare no conflict of interest. 19 Given the magnitude of the eruption and the scale of the 1 20 societal impacts a number of previous climate modeling stud- To whom correspondence should be addressed. E-mail: [email protected] www.pnas.org/cgi/doi/10.1073/pnas.XXXXXXXXXX PNAS | September 1, 2020 | vol. XXX | no. XX | 1–9 Table 1. Summary of experiments performed to reconcile the climate 46 aerosols (which scatter back sunlight and cool the planet (18)) and ozone response to the Mt. Samalas eruption. Nsims refers to the 47 are formed from the oxidation of SO in the gas phase by the 2 number of ensemble members performed for each set of simulations. 48 hydroxyl radical (OH). Volcanic sulfate aerosols grow through In all simulations the same day of year and altitude is used to in- 49 microphysical processes of condensation and coagulation (19), ject the volatile gases from the eruption of Mt. Samalas. BOTH-SO2 50 with the size of the aerosols being important for their climate refers to the SO2 only experiments, which are averaged in further analysis. Percentages indicate the fraction of emissions compared 51 effects (20). There is an important self-limiting effect of sulfur to those estimated in (4) 52 oxidation in the stratosphere (17, 21). As the size of the SO2 53 injected into the stratosphere increases, the ability of the OH . Ensemble Nsims Stratospheric Injection 54 to react with SO2 decreases; slowing down the rate of aerosol SO2 / Tg(S) HCl / Tg HBr / Tg 55 production, causing sulfate to be produced over a longer time- HI-HAL (20%) 6 142.2 46.68 0.263 56 period and leading to a longer-lived forcing. In addition, the LO-HAL (1%) 6 94.8 2.33 0.013 57 reduction in OH can cause levels of ozone to increase in the BOTH-SO2 HI-SO2 (90%) 6 142.2 0.00 0.000 58 lower stratosphere (21) where it can act as a climate warming BOTH-SO2 LO-SO2 (60%) 6 94.8 0.00 0.000 59 agent. 60 One aspect that affects our understanding of the effects of 61 the Mt. Samalas eruption, and potentially other future erup- 62 tions, is in the role of co-emitted halogens. The impacts of As a prior study showed best model agreement for modeled 107 63 volcanic halogens on stratospheric ozone were first discussed surface temperature with reconstructed surface temperature 108 64 by Stolarski and Cicerone (22). When halogens enter the from tree-rings was with an eruption occurring between May 109 st 65 stratosphere they contribute to the catalytic destruction of and July (8), a 1 June eruption was selected for our sim- 110 66 ozone (23) and lead to commensurate impacts on the composi- ulations at the latitude and longitude corresponding to the 111 67 tion and chemistry of the troposphere. Tie and Brasseur (24), location of Mt. Samalas (8.5 °S, 116.3 °E). The gases listed in 112 68 showed that there is significant sensitivity in the response of Table1 were injected between 19 km and 34 km altitude (see 113 69 stratospheric ozone following eruptions the size of 1991 Mt. Methods). 114 70 Pinatubo eruption to the background chlorine levels. Kutterolf For each of the emissions scenarios six simulations were 115 71 et al. (25) provided renewed interest in the topic of volcanic performed (ensemble members) to investigate the role of inter- 116 72 halogens and several studies have recently been published us- nal variability.
Recommended publications
  • 435 the Date of the Minoan Santorini Eruption
    THE DATE OF THE MINOAN SANTORINI ERUPTION: QUANTIFYING THE “OFFSET” Felix Höflmayer Deutsches Archäologisches Institut, Orient-Abteilung, Peter-Lenné-Strasse 32, 14195 Berlin, Germany. Email: [email protected]. ABSTRACT. Despite many recent attempts to settle the dispute concerning the absolute date of the Minoan Santorini erup- tion, there are still differences between some archaeologists and scientists on the absolute dates and the reliability of radio- carbon dating. The recent publication of over 200 new 14C dates for dynastic Egypt rules out a major flaw in the historical chronology of Egypt and proves the reliability of 14C dating in the Nile Valley. Therefore, the student of Aegean archaeology and eastern Mediterranean interconnections is still confronted with an archaeologically based conventional, or “low,” chro- nology and a 14C-backed “high” chronology. New 14C determinations from different sites of the Aegean support the high chronology for the Late Minoan (LM) IA, while recent re-evaluation of LM IB determinations are slightly higher but more or less in agreement with archaeological estimations. The present contribution reviews archaeological and scientific data for the LM IA period and argues that a reduced (~30 to 50 yr) offset between archaeological and 14C dates for the Minoan San- torini eruption may be possible, thus offering new perspectives for potential solutions for this problem. INTRODUCTION: “HIGH” AND “LOW” CHRONOLOGY AND RECENT DEVELOPMENTS Since the mid-1980s, a date for the Minoan eruption of Santorini 100 to 150 yr higher than previ- ously argued was suggested by radiocarbon dating and climatic events recorded in ice cores and tree rings from North America, Europe, and the eastern Mediterranean, placing the end of the Late Minoan (LM) IA period in the second half of the 17th century BCE (Betancourt 1987; Manning 1988, 1989a,b, 1990a,b, 1999, 2007, 2009; Michael and Betancourt 1988a,b; Friedrich et al.
    [Show full text]
  • And S-Wave Velocity Structures and the in Uence Of
    P- and S-wave Velocity Structures and the Inuence of Volcanic Activities in the East Java Area from Seismic Tomography Syawaldin Ridha Department of Physics, Universitas Brawijaya, Indonesia Sukir Maryanto ( [email protected] ) Universitas Brawijaya https://orcid.org/0000-0002-1882-6818 Agustya A. Martha Meteorological, Climatological, and Geophysical Agency, Indonesia Vanisa Syahra Department of Physics, Universitas Brawijaya, Indonesia Muhajir Anshori Meteorological, Climatological, and Geophysical Agency, Indonesia Pepen Supendi Meteorological, Climatological, and Geophysics Agency, Indonesia Sri Widiyantoro Bandung Institute of Technology: Institut Teknologi Bandung Research Letter Keywords: seismic tomography, East Java, Vp/Vs, Vp, Vs, partial melting Posted Date: May 6th, 2021 DOI: https://doi.org/10.21203/rs.3.rs-438689/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/24 Abstract Indonesia is one of the most interesting targets for seismic tomographic studies due to its tectonic complexity. The subduction zone was formed when the Indian oceanic plate was subducted beneath the Eurasian continental plate. This activity caused the formation of volcanoes along the Sunda Arc, including the area of East Java. In this study, we aim to identify the inuence of volcanic activities which extends from the west to the east of East Java. We used the data of 1,383 earthquakes, recorded by the 22 stations of the Indonesia Tsunami Early Warning System (InaTEWS) seismic network. We relocated the earthquakes and conducted a tomographic study using SIMULPS12. We then explored the anomalies of P- and S-wave velocities and Vp/Vs ratio. The low-velocity zone was observed in the volcanic area related to the partial melting zone or magma chamber with high Vp/Vs.
    [Show full text]
  • The Minoan Eruption of Santorini, Greece
    The Minoan eruption of Santorini, Greece A. BOND & R. S. J. SPARKS SUMMARY The Minoan eruption of Santorini produced logical contrasts in the mass flows which pro- the following sequence of deposits- a plinian duced them. Grain size analyses show wide pumice fall deposit, interbedded surtseyan-type ranges in the lithic contents of the different ash fall and base surge deposits, mud-flow types of deposit: ignimbrite (35--6o~), mud- deposits and ignimbrite interbedded with very flows (2o-3o ~) and the pyroclastic fall and coarse, well-sorted flood deposits. The variation base surge deposits (4-x 5 %). The ignimbrite is of thickness and grain size in the plinian deposit enriched in crystals, complemented by deple- indicates a vent x km west of Thera town. The tion in fine air-fall ash beds that interstratify base surges and surtseyan-type activity is in- with the ignimbrite. The gas velocity of the terpreted as the result of sea water entering the plinian phase is estimated as 55o m/s, the erup- magma chamber. The poorly sorted mud-flow tion column height as greater than 2o km deposits and ignimbrite are distinguished on and it is shown that only particles of 2 mm their grain size, temperature and morphological could have reached Minoan Crete. characteristics, which indicate substantial rheo- I N T H ~. L AT E B R ON Z F. A o E a paroxysmal eruption took place on Santorini Vol- cano, referred to as the Minoan eruption after the Minoan civilisation which inhabited the island at that time. The eruption produced a great volume of pumice and ash and resulted in the formation of the present day caldera which measures x 1- 5 × 8 km and probably had a catastrophic effect on the people living in the southern Aegean.
    [Show full text]
  • Volcanic Eruption Impacts Student Worksheet
    Volcanic Eruption Impacts Student Worksheet Explosive and Effusive Volcanoes The type of volcanic eruption is largely determined by magma composition. Flux-mediated melting at subduction zones creates a felsic magma with high levels of carbon dioxide and water. These dissolved gases explode during eruption. Effusive volcanoes have a hotter, more mafic magma with lower levels of dissolved gas, allowing them to erupt more calmly (effusive eruption). Sinabung (Indonesia) Mount Sinabung is a stratovolcano located 40 km from the Lake Toba supervolcano in North Sumatra. It lies along the Sunda Arc, where the Indo-Australian plate subducts beneath the Sunda and Burma plates. After 1200 years of dormancy, Sinabung began erupting intermittently in 2010. Major eruptions have occurred regularly since November 2013. In November and December 2015, ash plumes reached 6 – 11 km in height on multiple occasions. Pyroclastic flows and ashfall blanketed the region in January 2014 and lava flows travelled down the south flank, advancing 2.5 km by April 2014. Pyroclastic flows in February 2014 killed 17 people in a town 3 km from the vent. In June 2015, ash falls affected areas 10 – 15 km from the summit on many occasions. A lahar in May 2016, caused fatalities in a village 20 km from Sinabung. Pyroclastic flows occurred frequently throughout 2016 and 2017 Eruption of Sinabung 6 October 2016 Major eruptions occurred in 2018 and 2019. In (Y Ginsu, public domain) February 2018, an eruption destroyed a lava dome of 1.6 million cubic metres. At least 10 pyroclastic flows extended up to 4.9 km and an ash plume rose more than 16 km in altitude.
    [Show full text]
  • Apatite Sulfur Isotope Ratios in the 1257 Samalas Eruption (Indonesia)
    Goldschmidt2020 Abstract Apatite sulfur isotope ratios in the 1257 Samalas eruption (Indonesia) R. ECONOMOS1*, Y. JACKSON1, S. DING2, A. FIEGE3, C.VIDAL4, I. PRATOMO5, A. T. HERTWIG6, M.-A. LONGPRÉ2 1Southern Methodist University, Dallas, TX, USA (*[email protected]) 2CUNY Queens College, Flushing, NY, USA 3American Museum of Natural History, New York, NY, USA 4University of Cambridge, Cambridge, UK 5Geological Museum, Bandung, Indonesia 6UCLA, Los Angeles, CA, USA The 1257 eruption of Mt. Samalas in Indonesia produced sulfate anomalies in bi-polar ice cores that are ~2 times larger than those of the 1815 eruption of neighboring Tambora volcano, despite the similar magnitude of both eruptions1. The build-up of such a large volume of eruptible sulfur is likely to be related to pre-eruptive degassing and magma redox conditions. Information about these processes can be preserved in apatite crystals, which integrate sulfur as a trace element at concentrations that allow for δ34S isotope ratio characterization in-situ via Secondary Ionization Mass Spectrometry2. We investigated apatite crystals occuring as inclusions in plagioclase (Pl) crystals and as apatite microphenocrysts in contact with matrix glass, from trachydacitic pumices of both the climactic 1257 eruption and an earlier 2550 B.P. event3. 34 The lowest δ S(CDT) value observed in all four sample groups is 8.5‰, which we interpret as representative of evolved magmas entering the sub-volcanic system. Pl-hosted apatite 34 inclusions from the 1257 eruption range up to a δ S(CDT) of 11‰, while microphenocrysts reach 16‰, consistent with inclusions capturing an earlier stage of magma evolution.
    [Show full text]
  • Title Characteristics of Seismicity Distribution Along the Sunda Arc
    Characteristics of Seismicity Distribution along the Sunda Arc: Title Some New Observations Author(s) GHOSE, Ranajit; OIKE, Kazuo Bulletin of the Disaster Prevention Research Institute (1988), Citation 38(2): 29-48 Issue Date 1988-06 URL http://hdl.handle.net/2433/124954 Right Type Departmental Bulletin Paper Textversion publisher Kyoto University Bull. Disas. Prey. Res. Inst., Kyoto Univ., Vol. 38, Part 2, No. 332, June, 1988 29 Characteristics of Seismicity Distribution along the Sunda Arc: Some New Observations By Ranajit GHOSEand Kazuo OIKE (Manuscript received March 7, 1988) Abstract Spatio-temporal variations of earthquake activity along the Sunda arc were investigated. We prepared a strain release map for this century. Adjacent to the zones of high strain release, presence of seismically quiet zones was noted. A careful inspection of the depth distribution of the earthquakes revealed that in the eastern Sunda arc, possibly there exists a zone of scarce seismicity at an interme- diate depth. We discussed the probable implications. We also analysed the patterns of temporal distributions of earthquakes at the three different seismotectonic provinces of the Sunda arc—Sumatra, Java, and the Lesser Sunda Islands. We could clearly see that, although the causative geodynamic situations for seismicity vary significantly in space along the length of the arc, the period of increase or decrease in seismicity is largely space invariant. The locally differing levels of seismicity are superposed on the common background of long period seismicity fluctuation. Finally, clustering of seismicity at some patches along the Sunda arc was studied with respect to the altimetric gravity anomaly data. We noted some apparent conformities.
    [Show full text]
  • SR 53(7) 28-29.Pdf
    SHORT FEATURE NIKHILANAND PANIGRAHY Mt. Tambora Volcano with its caldera after eruption HE Sun is our default light option. But what happens when But the effect was not just limited to Indonesia. The Tit does not shine? destruction was spread far and wide. Even European countries In the year 1815, the repercussion of a very unusual event could not be spared of its evil consequences. The atmosphere in was felt worldwide for a fairly long period. The incident was the West was covered by the volcanic ash of Tambora, as a result a volcanic eruption in Mount Tambora. In Indonesia, there is of which the sun-rays could not reach the surface of the Earth. an island Sumbawa in its peninsula. This is a part of the Sunda Due to reduction of solar warmth, heavy snowfall and fatal frost islands that forms a segment of Sunda Arc. This is recognised as was found even during June to August, 1816. There was intense a string of volcanic islands. cold and the situation turned worse as famine-conditions In Mount Tambora, there exists a stratovolcano. It contains developed in European and North American countries. lava, pumice, volcanic ash, other materials and different gases. Many people from England started to rush to the lake side In many cases, a volcano looks like a cone, as the hot liquid lava of Geneva, Switzerland for warm climatic conditions. Such a emerging from it cannot fl ow to long distances away from the discouraging state of affair and unfair weather continued for vent or opening of the volcano, due to large viscosity.
    [Show full text]
  • Dating the Santorini Eruption
    Zurich Open Repository and Archive University of Zurich Main Library Strickhofstrasse 39 CH-8057 Zurich www.zora.uzh.ch Year: 2014 The olive-branch dating of the Santorini eruption Cherubini, Paolo ; Humbel, Turi ; Beeckman, Hans ; Gärtner, Holger ; Mannes, David ; Pearson, Charlotte ; Schoch, Werner ; Tognetti, Roberto ; Lev-Yadun, Simcha Abstract: The date of the volcanic eruption of Santorini that caused extensive damage toMinoan Crete has been controversial since the 1980s. Some have placed the event in the late seventeenth century BC. Others have made the case for a younger date of around 1500 BC. A recent contribution to that controversy has been the dating of an olive tree branch preserved within the volcanic ash fall on Santorini. In this debate feature Paolo Cherubini and colleagues argue that the olive tree dating (which supports the older chronology) is unreliable on a number of grounds. There follows a response from the authors of that dating, and comments from other specialists, with a closing reply from Cherubini and his team. DOI: https://doi.org/10.1017/s0003598x00050365 Posted at the Zurich Open Repository and Archive, University of Zurich ZORA URL: https://doi.org/10.5167/uzh-154340 Journal Article Published Version Originally published at: Cherubini, Paolo; Humbel, Turi; Beeckman, Hans; Gärtner, Holger; Mannes, David; Pearson, Charlotte; Schoch, Werner; Tognetti, Roberto; Lev-Yadun, Simcha (2014). The olive-branch dating of the Santorini eruption. Antiquity, 88(339):267-273. DOI: https://doi.org/10.1017/s0003598x00050365 Bronze Age catastrophe and modern controversy: dating the Santorini eruption The date of the volcanic eruption of Santorini that caused extensive damage to Minoan Crete has been controversial since the 1980s.
    [Show full text]
  • The Theran Eruption and Minoan Palatial Collapse: New Interpretations Gained from Modelling the Maritime Network
    The Theran eruption and Minoan palatial collapse: new interpretations gained from modelling the maritime network. Carl Knappett, Department of Art, University of Toronto, Ray Rivers, Department of Physics, Imperial College London, Tim Evans, Department of Physics, Imperial College London, UK. What was the effect on Late Minoan civilisation of the catastrophic destruction of Akrotiri on Thera (Santorini) by volcanic eruption? Not much, according to the evidence for continuing prosperity on Crete. But the authors mobilise their ingenious mathematical model (published in Antiquity 82, 1009-24), this time to show that the effects of removing a major port of call could have impacted after an interval, as increased costs of transport gradually led to ever fewer routes and eventual economic collapse. Keywords: Bronze Age, Aegean, Minoan, Crete, Thera, Santorini, Akrotiri, maritime communications, network analysis Introduction We seek in this paper to provide a novel perspective on the possible causes for the demise of Cretan Bronze Age palatial society c. 1500 BC using a mathematical model developed from a previous study (Knappett et al 2008). Various explanations have been proposed for the collapse - a tsunami generated by the Theran eruption destroying the Minoan fleet, an invasion by Mycenaean mainlanders, or internal socio-political turmoil and unrest (Rehak & Younger 2001). Today many scholars would see many such factors, natural and social, combining in some way. And although most do see some role for the Theran eruption, the fact that it took place some 50-100 years before the collapse makes it difficult to envisage what that role might have been. One persuasive argument sees indirect though insidious effects, both economic (such as ashfall polluting water supply and compromising crop yields in east Crete) and social, such as ideological uncertainty (Driessen & Macdonald 1997: 89-98).
    [Show full text]
  • Seismic Activity Around and Under Krakatau Volcano, Sunda Arc: Constraints to the Source Region of Island Arc Volcanics
    SEISMIC ACTIVITY AROUND AND UNDER KRAKATAU VOLCANO, SUNDA ARC: CONSTRAINTS TO THE SOURCE REGION OF ISLAND ARC VOLCANICS ALEŠ ŠPI ČÁK , VÁCLAV HANUŠ AND JIŘÍ VAN ĚK Geophysical Institute, Academy of Sciences of the Czech Republic, Prague, Czech Republic * ABSTRACT There is general agreement that calc-alkaline volcanic rocks at convergent plate margins are genetically related to the process of subduction (Ringwood, 1974; Maaloe and Petersen, 1981; Hawkesworth et al., 1997). However, opinions on the mode and site of generation of primary magma for island arc volcanism differ substantially. The site of generation of calc-alkaline magma is thought to be either in the mantle wedge (Plank and Langmuir, 1988; McCulloch and Gamble, 1991) or in the subducting slab (White and Dupré, 1986; Defant and Drummond, 1990; Edwards et al., 1993; Ryan and Langmuir, 1993). We present seismological evidence in favour of the latter concept. A distinctive seismicity pattern around and under the Krakatau volcano was identified during systematic studies of the SE Asian convergent plate margins by means of global seismological data. A column-like cluster of events, probably associated with the dynamics of the volcano, is clearly separated from the events in the Wadati-Benioff zone. The accuracy of hypocentral determinations of the events of the cluster does not differ from the accuracy of the events belonging to the subducting slab. The depths of the cluster events vary from very shallow to about 100 km without any apparent discontinuity. On the other hand, there is a pronounced aseismic gap in the Wadati-Benioff zone directly beneath the volcano at depths between 100-150 km.
    [Show full text]
  • Petrogenesis of Rinjani Post-1257-Caldera-Forming-Eruption Lava Flows
    Indonesian Journal on Geoscience Vol. 3 No. 2 August 2016: 107-126 INDONESIAN JOURNAL ON GEOSCIENCE Geological Agency Ministry of Energy and Mineral Resources Journal homepage: hp://ijog.geologi.esdm.go.id ISSN 2355-9314, e-ISSN 2355-9306 Petrogenesis of Rinjani Post-1257-Caldera-Forming-Eruption Lava Flows Heryadi Rachmat1,2, Mega Fatimah Rosana1, A. Djumarma Wirakusumah3, and Gamma Abdul Jabbar4 1Faculty of Geology, Padjadjaran University Jln. Raya Bandung - Sumedang Km. 21, Jatinangor, Sumedang, Indonesia 2Geological Agency Jln. Diponegoro No. 57, Bandung, Indonesia 3Energy and Mineral Institute Jln. Gajah Mada, Karangboyo, Cepu, Kabupaten Blora, Indonesia 4Hokkaido University, Kita 10, Nishi 8, Sapporo, Japan Corresponding author: [email protected] Manuscript received: March 7, 2016; revised: May 17, 2016; approved: June 29, 2016; available online: August 2, 2016 Abstract - After the catastrophic 1257 caldera-forming eruption, a new chapter of Old Rinjani volcanic activity began with the appearance of Rombongan and Barujari Volcanoes within the caldera. However, no published petrogenetic study focuses mainly on these products. The Rombongan eruption in 1944 and Barujari eruptions in pre-1944, 1966, 1994, 2004, and 2009 produced basaltic andesite pyroclastic materials and lava flows. A total of thirty-one samples were analyzed, including six samples for each period of eruption except from 2004 (only one sample). The samples were used for petrography, whole-rock geochemistry, and trace and rare earth element analyses. The Rombongan and Barujari lavas are composed of calc-alkaline and high K calc-alkaline porphyritic basaltic andesite. The magma shows narrow variation of SiO2 content that implies small changes during its generation.
    [Show full text]
  • Reconstruction of the 2018 Anak Krakatau Collapse Using Planetscope Imaging and Numerical Modeling
    Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2020 Reconstruction of the 2018 Anak Krakatau collapse using PlanetScope imaging and numerical modeling Davide Saviano Michigan Technological University, [email protected] Copyright 2020 Davide Saviano Recommended Citation Saviano, Davide, "Reconstruction of the 2018 Anak Krakatau collapse using PlanetScope imaging and numerical modeling", Open Access Master's Thesis, Michigan Technological University, 2020. https://doi.org/10.37099/mtu.dc.etdr/989 Follow this and additional works at: https://digitalcommons.mtu.edu/etdr Part of the Geology Commons, Geomorphology Commons, Tectonics and Structure Commons, and the Volcanology Commons RECONSTRUCTION OF THE 2018 ANAK KRAKATAU COLLAPSE USING PLANETSCOPE IMAGING AND NUMERICAL MODELING By Davide Saviano A THESIS Submitted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE In Geology MICHIGAN TECHNOLOGICAL UNIVERSITY 2020 © 2020 Davide Saviano This thesis has been approved in partial fulfillment of the requirements for the Degree of MASTER OF SCIENCE in Geology. Department of Geological & Mining Engineering & Sciences Thesis Co-Advisor: Dr. Simon A. Carn Thesis Co-Advisor: Dr. Gianluca Groppelli Committee Member: Dr. Roohollah R. Askari Department Chair: Dr. John S. Gierke Table of Contents Abstract ............................................................................................................................... v 1 Introduction ..............................................................................................................
    [Show full text]