Eruption products of the 1883 eruption of Krakatau and their final settlement ANNALS OF GEOPHYSICS, 58, 2, 2015, S0220; doi:10.4401/ag-6529 S0220 Eruption products of the 1883 eruption of Krakatau and their final settlement Izumi Yokoyama The Japan Academy, Ueno Park, Tokyo, Japan ABSTRACT Firstly the volume of pyroclastic ejecta during the 1883 eruption of Krakatau is re-examined. To revise the volume of flow deposits, the au- thor basically follows Verbeek’s observation while to estimate the fall de- posits, as the last resort, the author assumes that volume ratios fall/flow are common to similar caldera eruptions, and the ratios determined by the caldera- forming eruptions of Novarupta and Pinatubo are applied to the Krakatau eruption. Verbeek’s estimation of the total volume of ejecta, 12 km3 is revised to 19 km3. This is significantly different from the volume of disrupted volcano edifice, 8 km3. Such a result does not support the pred- ecessors’ hypothesis that calderas are formed by collapses of volcano edi- fices into magma reservoirs in replacement of the total ejecta. Through the discussion on the volume estimation of volcanic ejecta on and around Krakatau, the author recognizes that such estimation should be originally very difficult to attain enough accuracy. Much importance of “caldera deposits” to post-eruption settlements of the ejecta is emphasized. In re- lation to caldera formation, mechanical stability of a cavity in the crust is discussed. Lastly, upon the basis of subsurface structure, especially caldera deposits, a structural image of Krakatau caldera is presented. 1. Introduction We have experienced caldera-forming eruptions four times since the 19th century, Tambora, Sumbawa in 1815, Krakatau, Sunda Straits in 1883, Novarupta, Alaska in 1912, and Pinatubo, Luzon in 1991. The Tamb- ora eruption was exceptionally violent but not pur- posefully observed. Sigurdsson and Carey [1989] carried out field studies there and clarified that the eruption was a complex series of events for about 10 days devel- oping its explosive phases. This eruption was rather ex- ceptional in its sequences while the other three were observed more or less at scientific levels of those days and have been well studied until now. In the present paper, the 1883 eruption of Krakatau shall be discussed partly in the light of the last two caldera-forming erup- tions although there exist many calderas formed in geo- logical ages. The 1883 Krakatau eruption has been studied by several distinguished researchers of the 20th century. Some of them have kept a conception that calderas were formed by collapses of volcano edifices into magma reservoirs in replacement of erupted magmatic mate- rial by reason of the paucity of lithic material around the calderas, and this conception presupposes equality of volume between discharged magmatic ejecta and disrupted volcano edifices. In order to verify the equal- ity, we have to estimate the volume with reasonably high accuracy not in an order of magnitude: generally such estimation is impossible even now. The difficulties in studies of the 1883 eruption of Krakatau are as follows: (a) The eruption occurred in the midst of the sea and was accompanied with the dark of huge explosion clouds, and under no instrumental observations except a few accidental ones: a tide gauge and a pressure meter at Batavia located at about 160 km distance from the volcano. Consequently the sequence of this eruption or caldera formation has been not fully clarified. (b) At that time, surveys of pyroclastic ejecta under the sea had been inexperienced. Around Krakatau Is- lands, the new eruption deposits suffered from effects of the Sunda Current: it was important to make a sur- vey as soon as possible after the eruption to minimize post-depositional erosions. Verbeek [1885] may have pre- ferred an immediate survey to an accurate one. At present, an overall estimation of pyroclastic de- posits in the sea area produced by the 1883 eruption is no more possible because of post-depositional erosions for more than 100 years, to say nothing of the land de- posits. Thus, only Verbeek’s observational data are uniquely available but have not been re-examined. We should recognize utility of the past estimation and pro- mote the discussion on the 1883 eruption. The topographic changes of Krakatau and its sur- Article history Received March 19, 2014; accepted April 8, 2015. Subject classification: Krakatau caldera, Volume of the ejecta, Pyroclastic-fall and -flow, Caldera deposit, Magma reservoir. rounding islands, surveyed immediately after the 1883 eruption by Judd [1888, fig. 11] are shown in Figure 1. These changes were caused by deposition of volcanic ejecta and explosive destructions. According to Judd, many of the topographic discrepancies at the southern limit of Krakatau were due to the imperfection of the old survey. This suggests difficulty in re-occupying the survey spots on the sea and in volume estimation of the submarine ejecta by using such charts. After Verbeek [1885], several authors have repeated volume-estimation of pyroclastic ejecta from the 1883 eruption of Krakatau by their own manners. Their re- sults shall be summarized in Section 3.3. About one month after the eruption, Verbeek started surveys of the volcanic ejecta in the sea by the techniques of those days. When we take the effects of post-depositional erosions in the Sunda Straits into consideration, we should appreciate the merit of Verbeek’s early survey though his method was not necessarily complete. 2. Volume estimation of pyroclastic-flow deposits of the 1883 eruption 2.1. Verbeek’s volume estimation of pyroclastic ejecta Verbeek [1885] did not explicitly distinguish fall de- posits from flow deposits. To estimate the total amount of ejecta around Krakatau Islands, he divided the area into two districts at and around the newly formed caldera; one is a circular area of 15 km in radius (Figure 2a) and the other is a wider area surrounding the for- mer and extending toward the N and the SW to the maximum distance of about 1300 km (Figure 2b). In the former, Pieces I ~ IV, and Rings a ~ c, he com- pared the two bathymetric charts surveyed before and after the eruption, and in the latter, Rings d, e, f, and h, YOKOYAMA 2 Figure 1. Chart of Krakatau and the surrounding islands, from the chart surveyed immediately after the eruption [from Judd 1888, fig. 11]. The shaded areas show the form of the islands drawn in the old chart. Dotted lines show sandbanks and lines of breakers. “Danan” is added by the present author. Figure 2. Division maps of Verbeek’s survey of ejecta in and around the Sunda Straits. (a) Survey divisions of the inner area of Krakatau Is- lands: a simplified map of Verbeek’s figure [1885, fig. 4]. Pieces I ~ IV are within a circle of 15 km radius around Krakatau Islands, and Rings b, c and d successively continue to Ring a. ⓢ and ⓒ denote Steers and Calmeyer reefs, respectively. A hatched area is the caldera depres- sion; Verbeek’s method is not applicable there. (b) Survey divisions in the outer area of Krakatau Islands. The innermost circle (r = 50 km) includes Rings a, b, c and Pieces I ~ IV shown in Figure 2a. (a) (b) 3 and Rectangle g, Triangles k and l, he roughly pre- sumed thicknesses of the new deposits on the bases of the reports from the inhabitants and navigating ships in these areas. He actually surveyed the depth changes within ra- dial distance of 15 km from the center of the caldera (Pieces I ~ IV). On the other hand, he did not survey in the area outside of Piece IV, but assigned presumed thickness for each Division as shown in Table 1. In the table, the thickness of deposits discontinu- ously decreases (0.2 to 0.02 m) between Divisions c and d (50 km from the center). Under such pre- arranged condition, the present author assumes that the deposits within a circular area of 50 km in radius (or Pieces I ~ IV and Rings a ~ c) are mainly pyro- clastic-flow deposits including a small quantity of py- roclastic-falls while in the outer area, almost all the deposits are pyroclastic-falls. Such distinction is one of the last resorts as far as we follow Verbeek’s estima- tions, but may be not much unreasonable. In the fol- lowing, Verbeek’s results on pyroclastic-flow obtained only in a circular area of 15 km in radius (Pieces I ~ IV) shall be re-examined, and the fall deposits in other Rings (a ~ c) remain the same as Verbeek’s estimation. And the deposits in Divisions (d ~ l) are regarded as pyroclastic-falls, but the total volume of the fall deposits shall be re-estimated by analogical method with the other caldera-forming eruptions on land (Novarupta and Pinatubo) in Section 3.2. 2.2. Revision of Verbeek’s method to estimate pyroclas- tic-flow deposits Verbeek [1885] estimated volume of the ejecta by comparing the depths of the two charts of this area; namely, the Blommendal chart of 1874 and the Hydro- graaf chart after the 1883 eruption as shown in Figures 3a and 3b, where the isobaths are added every 25 m by the present author. Verbeek presented the sea depths before and after the eruption together on a sheet as shown in Verbeek’s Atlas [1885, fig. 4]. The present au- KRAKATAU ERUPTION PRODUCTS Division a b c d e f g h k l Thickness (in meter) 1 0.3 0.2 0.02 0.05 0.01 0.001 0.005 0.005 0.005 Table 1. Thickness of pyroclasic deposits in Division a ~ l shown in Figures 2a and 2b after Verbeek [1885, pp. 133 and 134]. Figure 3. British Admiralty Charts surveyed before and after the 1883 eruption of Krakatau. In both the figures, the isobaths are added every 25 m by the present author. (a) Chart surveyed in 1854. Spot depths are in fathom (1 fathom = 1.83 m). ⓢ and ⓒ denote the locations of Steers and Calmeyer reefs, respectively. P Perbuwatan, D Danan, R Rakata. (b) Chart surveyed in 1886. The central region collapsed to a depth of 250 m at the maximum. (a) (b) thor does not utilize this figure but compare both the depths in meshes of 1 km square as shown in Figure 4. 2.2.1. Plumbings in the 19th century Verbeek [1885] measured thicknesses of the ma- rine deposits of the Krakatau eruption by comparing the two bathymetric charts around Krakatau Islands, surveyed before and after the eruption; He relied on ba- thymetries by British and Dutch navies. Originally he had to directly measure thicknesses of the deposits, but such measurements require time, and hence he may have wished to obtain the result as early as possible to reduce the effects of post-depositional erosion by the current of the straits. This method is valid in case of thick deposits, several meters thick at least. For successful comparisons of two bathymetric charts surveyed at different periods, accurate re-locations of previous meas- uring points should be very desirable. In 1880s, re-loca- tions were rather doubtful, especially in remote areas. We should take appropriate errors into consideration when we treat with Verbeek’s results. The above com- parison of bathymetric charts is not applicable to the area where the seabed was destructed by the eruption. It is concluded that his results contain some errors due to incomplete methods and inaccurate charts, particu- larly one surveyed before the eruption. Verbeek [1885] did not describe any details of his bathymetries. Here, bathymettic methods in the 19th century shall be reviewed: echo-soundings were first adopted after 1920. In the 19th century, soundings at shallow sea were made usually by using a cable with a plumb of heavier than 7 pounds (3.2 kg). The accuracy of plumbings depends on verticality of plumb cables that is disturbed by the current thereabouts. Around Krakatau a few weeks after its eruption, the plumb of 3.2 kg weight may have penetrated the upper part of ejecta deposits, approximately 0.5 m in depth. By this method, it may have been rather difficult to measure the thickness of the deposits with high accuracy. At present, we can measure it directly by SCUBA core-samplings as used by Mandeville et al. [1996] or by seismic profilers and side-scan SONAR as used by Camus et al. [1991]. In those days, ships usually determined their posi- tions by sextants, and its accuracy depended on visibil- ity of landmarks around the ships. As shown in Figure 1, all Krakatau Islands drastically changed their shapes missing many remarkable topographies those were pre- viously favourable targets of positioning. In sum, the overall errors in the comparison of the two surveys around Krakatau Islands consisted of those of sound- ings and those of positioning, and may have reached possibly more than a few meters. Here, we will follow Verbeek’s method to re-estimate the volume of the 1883 ejecta, i.e. comparisons of the two charts, before and after the eruption. His data were published in his book of 1885, and fortunately we can effectively refer various data of the 1883 eruption and the related prob- lems to a book compiled by Simkin and Fiske [1983]. 2.3. Revised volume of the pyroclastic-flow deposits Verbeek may have obtained the differences be- tween the two soundings on the charts. In Figure 2a, an area hatched by Verbeek [1885] was depressed by the eruption. This is the most important part in our dis- cussion, but the resurveys could not be realized because the bottom topographies changed by the explosions. He grouped the central area of deposits into four Pieces (I ~ IV), assigning thicknesses of the deposits for each Piece as shown in Table 2. In the present paper, to examine Verbeek’s estima- tions, the author follows his method in principle and use similar two charts surveyed in 1854 and 1886, be- fore and after the 1883 eruption, by British Admiralty and referred by Wharton [1888, plts. XXXII and XXXIII]. The charts are reproduced in Figures 3a and 3b where the depths are in fathom (1 fathom = 1.83 m) and the isobaths are added every 25 m by the present author. Strictly speaking, the two charts are not the same as Verbeek’s ones. On this point, Judd [1888] YOKOYAMA 4 Figure 4. Distribution of thickness of volcanic ejecta, newly de- duced from Figures 3a and 3b. Dotted lines show sandbanks and lines of breakers. The broken line in the central part contours an isobath of 150 m after the 1883 eruption. A circle of 15 km radius corresponds to that of Verbeek’s figure (Figure 2a). A mesh is 1 km × 1 km. ○—: < 0 (thin deposits); No numerals: 0~5 m (deposits of a few meter thick); 1: 6~15 m, roughly around 10 m thick; 2: 16~25 m, roughly around 20 m thick 20 m; 3: 26~35 m, roughly around 30 m thick; 4: > 35 m. 5 stated that the English chart was surveyed in 1854 and the Dutch one in 1874, but surveying officers of both the navies made additions and corrections to both the charts from time to time. Therefore, these two charts should be actually the same as Verbeek’s ones. The old chart has less numbers of plumbings and the depths in the two charts (Figures 3a and 3b) do not always correspond with each other. Several depths in the old chart were omitted from the new chart proba- bly because of their uncertainty. It is difficult to draw the isobaths automatically. To improve scatterings and to get smooth distribution of the sea depths, the author newly contoured isobaths every 25 m (about 14 fath- oms) by interpolations. Of course, these isobaths have some ambiguity, and drawings of every 10 m contours are impossible because of paucity and ambiguity of the measurements. In order to estimate the depth changes objectively, both the charts (Figures 3a and 3b) are divided into meshes of 1 km square, and the median value in each mesh is determined as the average depth, and then the corresponding depths of the two charts are compared in each mesh. Thus, the changes in depth of the sea bot- tom due to the new deposits or the bottom displacements are represented in Figure 4 where mark ○— indicates a change smaller than ±5 m, numbers 1, 2, 3, and 4 indi- cate changes of around 10, 20, 30, and 40 m, respectively. As mentioned above, around Krakatau Islands, the pyroclastic-flow deposits are distinguished from the py- roclastic-fall deposits at a boundary of 50 km radius from the center of the explosion. The flow-area is di- vided into 4 Pieces I ~ IV and 3 Rings a, b and c as shown in Figure 2a. Verbeek’s estimations at each divi- sion are revised as follows: Pieces I ~ IV: the volume of deposits in the 4 Pieces is estimated on the base of Figure 4 and the results are tabulated in Table 2; Piece II is not uniform in thick- ness of the deposits, and divided into two by thickness. Piece IV includes the depression of the caldera and a triangular part (hatched part in Figure 2a), of which area is about 38 km2. We cannot apply Verbeek’s method to this depression. Nevertheless Verbeek [1885] assigned 6 m as the thickness of the ejecta in all Piece IV. He did not pay any special attention to the depression that, in the present writer’s opinion, should be the most im- portant part of the Krakatau area after the eruption. We cannot exclude the possibility that many lithic frag- ments deposit there. Its bottom of 200 m deep covers an area of about 11.8 km2 in Figure 5. In Table 2, Ver- beek divided Piece IV into three parts, all 6 m thick de- posits. The area of the caldera depression (hatched area in Figure 2a) is 38 km2 and here the author assumes thickness of the caldera deposits to be 50 m. This is not so overestimated considering that the maximum thick- ness of the deposits is 40 m in Figure 4. Thus Verbeek’s estimates are revised and modified as shown in the last column of Table 2. The revised volume of Pieces I ~ IV is converted into dense rock equivalent (DRE) as 10.8 km3 × 2.1/2.6 = 8.7 km3 DRE. In the table, Piece IV covers the caldera where there should be “caldera deposits” composed of fall- backs, i.e. juvenile ejecta and lithic fragments: it is im- KRAKATAU ERUPTION PRODUCTS Unit Surface (km2) Thickness (m) Volume (km3) Piece I 54 20 1.08 Piece II 45 198 30→ 20 30→ 10 0.90 1.98 Piece III 36 40→ 35 1.26 Piece IV 30 38 (caldera area) 306 6→ 20 6→ 50 6→ 10 0.60 1.9 3.06 Total Σ 707 10.8 ± 3 (bulk) 8.7 (DRE) Ring a (r =15~22.5 km) 883 1 0.883 Ring b (r =22.5~40 km) 3434 0.3 1.030 Ring c (r = 40~50 km ) 2826 0.2 0.565 Total Σ 7143 2.48 ± 0.5 (bulk) 2.0 (DRE) Grand total ΣΣ 7850 10.7 (DRE) Table 2. Revised values of Verbeek’s estimation in an area of 50 km in radius from the centre of the explosion. The deposits are mostly py- roclastic flows. Underlined values are revised ones. “caldera area” in Piece IV is the same as a hatched area in Figure 2a. possible for us to know their partitions at present. In the future, we shall get knowledge of the deposits by drillings at the spot. The caldera deposits shall be re- ferred to in Section 4.1. 2.3.1. Errors of Verbeek’s volume estimation Pieces I ~ III: Verbeek assigned thickness of the ejecta to 20, 30 and 40 m for the three Pieces, respec- tively. In other words, he measured these depths with a unit of 10 m and the error is ± 5 m. Thus, the error in estimation of the volume should be: (Area of Pieces I ~ III) × (Error in thickness) = (Error in volume): (54+243+36) km2 × (± 5 m) = ± 1.7 km3 ≈± 2 km3. Hence the total volume of Pieces I ~ III is expressed as (5.2 km3 ± 2 km3) in bulk. This error originates from only errors in measurement of the thickness, and varies with the locality. 2.3.2. Comments on the pyroclastic ejecta in the innermost region (a) The central area among the three islands was depressed remarkably reaching deeper than 200 m b.s.l. To this caldera area, Verbeek’s method is not applica- ble because the old depths were lost. Nevertheless, Ver- beek [1885] assigned thickness 6 m to the all areas of Piece IV including the collapsed part, or the caldera de- pression in Figure 2a. The present author cannot un- derstand the grounds for estimation of 6 m thick, and anew assigns thickness 50 m to the central area consid- ering the maximum thickness of the deposits shown in Figure 4, as mentioned above. Verbeek did not take thought of caldera deposits. (b) Danan has been assumed to be the main vent of the caldera forming explosions by some authors. However, the center of the existing depression is more probable as the explosion center because Danan was lo- cated outside the depression in Figure 6, and erupted as a parasite of pre-Krakatau volcano [Yokoyama 2014] before and during the 1883 climactic explosions. (c) A large amount of ejecta was deposited roughly 40 m thick at the maximum off the NW coast of Ver- laten (Sertung) Island in Figures 1 and 4. Stehn [1929] actually investigated the same deposition on the island. This means that the pyroclastic flows deposited on these areas after the large tsunami passed over there. (d) Scarcely the sea depths changed at the S of Rakata Island within the limits of surveying errors. This means that the large explosions were directed mainly north- and northwestward. (e) Reefs of Steers and Calmeyer: before the 1883 eruption, the two reefs were below the sea level as shown in Figure 3a, and after the eruption, they were observed above the sea between Lang and Sebesi Is- lands as in Figures 1, 3b and 4. These are located in Piece II as marked by letters and in Figures 2a&C&S YOKOYAMA 6 Figure 5. A part of Indonesian Chart around Krakatau Islands sur- veyed in 1970 by Kepala Dinas Hidro-Oseanografi. Spot depths are in meter. The isobaths of 200 and 250 m in the central part are drawn by the present author. Sertung and Rakata Kcl are another names of Verlaten and Lang, respectively. ⓢ and ⓒ near the cross- ing of the longitude and latitude lines denote Steers and Calmeyer reefs, respectively. A Anak Krakatau, B Bootsmans Rots. Figure 6. Topographic and bathymetric chart on and around Krakatau: a part of Indonesian Chart around Krakatau Islands sur- veyed in 1940. The isobaths are drawn every 50 m and the contour intervals every 20 m. AK Anak Krakatau, R Rakata, D Danan, P Per- buwatan. A circle of broken line shows the assumed caldera rim. 7 and 5. The changes in depth around these reefs were not so much, or less than 25 m. Originally these reefs were below the sea level because, on a chart (scale 1/500 000, 1970) shown in Figure 5, Steers and Calmeyer are 4.3 m and 6 m below the sea level, respectively. Ac- cording to Judd [1888], the ejecta deposited temporar- ily on both the reefs during the eruption, had been completely washed away in a few months by the action of the sea waves. In Atlas of Verbeek [1885, fig. 4] or Simkin and Fiske [1983, p. 180, plt. 4], Steers stood 3 m above sea level and Calmeyer did 6.5 m. Then, the ejecta are supposed to have deposited during the erup- tion on these reefs, 7.3 m and 12.5 m thick, respectively. These thicknesses are agreeable to those given in Figure 4. This means that the large tsunami on August 27 tra- versed this area before the pyroclastic ejecta deposited on these reefs. In other words, first the gigantic explo- sions occurred causing the large tsunami and strong air- waves, and later the pyroclastic-flows deposited on and around these reefs. Also this suggests importance of post-depositional erosions in the present discussion. Rings a, b and c (Figure 2a): we unwillingly cannot re-examine the data because Verbeek [1885] left the unique thickness of the deposits of 1, 0.3 and 0.2 m, re- spectively. These thickness were not actually measured but probably may have been reported by the inhabi- tants on the islands. Eventually, adding the volume of Rings a, b and c (2.0 km3 DRE) to the volume of Pieces I ~ IV, the revised volume of the pyroclastic-flow de- posits amounts to 10.7 km3 DRE as shown in Table 2. In the revised results, it is noticeable that about 80 % of the flow-deposits remain within 15 km distance from the center of the eruptions. 2.3.3. Supplements to Verbeek’s estimation of the flow deposits (a) Volume of disrupted parts of Rakata Island: volume of disrupted volcanic edifices and sea bottom by 1883 eruption was not estimated by Verbeek, but is also important to draw a balance sheet between the vol- ume of volcanic ejecta and disrupted volcanic edifices. Topographic maps of Rakata Island before the 1883 eruption were poor, and there were the three peaks on the island: Rakata (809 m), Danan (449 m) and Per- buwatan (120 m). By the eruption, the northern 2 ⁄ 3 of the island disappeared and Rakata decreased its height to 803 m. The volume of disrupted parts of Rakata Is- land above the sea level is calculated on the old and new topographic maps (Figures 3a and 6) as about 2.9 km3. (b) Volume of topographic depression beneath the sea level: the 1883 eruption formed a depression between Rakata and Verlaten Islands in addition to disruption of 2/3 Rakata Island, and almost simultaneously these de- pressions were filled with mixture of the ejecta and sea water during the volcanic paroxysms, and after the 1883 eruption, a depression of funnel shape with maximum depth 250 m was left among the three islands (Figure 5). This is Krakatau caldera. The topographic depres- sion under the sea level as of 1940 is shown in Figure 6 where probable caldera boundary is shown with a cir- cle of broken line. Intensive explosions of Danan volcano may have caused eastward extension of the central de- pression. Volume of the topographic depression is esti- mated at 5.1 km3 on the charts of Figures 3a and 6. This is “caldera depression”. In the case of Krakatau, the vol- ume of disrupted parts of Rakata Island above sea level (about 2.9 km3) is not included in the count. (c) Deposits on the three islands in the central area: Verbeek [1885] did not pay particular attention to vol- ume of the ejecta on the three islands. He assigned thickness of the ejecta on the three islands extrapolat- ing the depth changes measured at the adjacent sea in Figure 2a where the average thickness of Pieces II, III and IV covering the three islands is about 20 m in Table 2. Considering that thickness of the ejecta on the three islands is about 20 m on the average based on Stehn’s re- port [1929], we may assume the volume of deposits on the three islands as about 0.7 km3 DRE. This value was eventually included in Verbeek’s calculation. 3. Volume of pyroclastic-fall deposits of the 1883 eruption 3.1. Verbeek’s estimation of pyroclastic-fall deposits Originally volume estimation of fall deposits from volcanoes is problematic, especially at volcanoes in the midst of the sea. In the discussion of the 1883 eruption of Krakatau, we expect a large amount of pyroclastic-falls considering the anomalous “bigness” of the eruption [Simkin and Fiske 1983, p. 23] such as (1) a large volume of flow deposit, (2) energy coupled into the airwaves and tsunami, (3) the maximum eruptive column height: 50 km, (4) long duration of climatic eruption: about 23 hours. As already discussed in Section 2.1, according to Verbeek’s data given in Table 1, the ejecta in the area from Ring d to Triangle l, covering the sea and lands of 5 × 105 km2 in area (Figure 2b), are deemed as pyro- clastic-falls. This area is the outside of 50 km radius from the center of the caldera. He assigned 0.001~0.02 m to thickness of the deposits at each division (Table 1) as far as 1400 km distance, and obtained the total amount of the ejecta in this area as 3.6 km3 bulk (Table 4). This is extremely rough estimation in comparison with that of flow deposits. His data based on the reports from the inhabitants and navigating ships at scattered and distant spots around the volcano. KRAKATAU ERUPTION PRODUCTS To convert the bulk volume to DRE volume, we need density of fall deposits. Sarna-Woicicki et al. [1981, fig. 339] obtained density distribution of the May 18, 1980, fall deposits from Mt. St. Helens along axis of downwind lobe: the density of air-fall deposits varies between 0.5~0.1×103 kg/m3 with the distance from the volcano. Here, we adopt density 0.3×103 kg/m3 for the fall deposits. Thus Verbeek’s result of pyroclastic-fall de- posits is modified to 0.4 km3 DRE (Table 4) and the total volume of the flows and falls amounts to 12.2 km3 DRE. 3.2. Estimation of fall-deposits after the analogy of caldera-forming eruptions at Novarupta and Pinatubo vol- canoes When we try to revise the above Verbeek’s esti- mation, it is clear that we can no more collect system- atic data of the pyroclastic falls of the 1883 eruption. Here, to re-estimate the fall deposits, as the last resort, we assume that volume of pyroclastic-fall deposits is proportional to that of pyroclastic-flow deposits in the caldera forming eruptions provided that these erup- tions would occur by similar mechanisms and processes. The volume of pyroclastic-fall deposits of the 1883 Krakatau eruption shall be speculated after the analogy with well-studied caldera-forming eruptions, Novarupta, Alaska in 1912 and Pinatubo, the Philippines in 1991. These two volcanoes are similar to Krakatau in explo- sivity or magnitude of eruptions and properties of magmas (andesite ~ dacite ~ rhyolite), and both the two calderas are filled with lithic-rich caldera deposits. Their eruptions were later than Krakatau, and have been studied much more with modern knowledge and meth- ods, and furthermore both they are on land: the last is absolutely advantageous to estimate volume of the py- roclastic ejecta. There are many papers discussing their eruptions and transport processes of eruption deposits in detail. Some particulars of the three eruptions shall be comparatively discussed and are listed in Table 3. Novarupta: Hildreth [1983] estimated volume of the ash-flow tuff at 11~15 km3 bulk, that is converted to about 8.2 km3 DRE, and Fierstein and Hildreth [1992, appx. 1] did that of the tephra falls at 6.8 km3 DRE. Hilddreth [1983] assumed 5% of the total ash falls and flows to be lithic material that derived from the former ground. This lithic material fills the caldera depression that was formed by explosive excavations. And further he assumed the caldera deposits to be about 1 km3 in volume, a half of the original ground. Thus, the origi- nal ground was approximated to an inverted cone (fun- nel-shape) of 2 km deep and 2 km in top diameter. In other words, about a half of the caldera deposits is rea- sonably presumed to be lithic fallbacks at Novarupta. Pinatubo: the climatic eruption culminated with lithic-rich facies in the pyroclastic-flow deposits and re- sulted in formation of a caldera of 2.5 km in diameter. Scott et al. [1996] estimated volume of the pyroclastic- flow deposits issued in various directions from the crater at 2.1~3.3 km3 DRE by the valley-cross-section method. YOKOYAMA 8 Volcano Krakatau 1883 Novarupta∗ 1912 Pinatubo∗∗ 1991 Rocks andesite ryolite, andesite, dacite dacite (a) Magma volume of pyroclastic- flow deposit (km3 DRE) 10.7 8.2 2.1~3.3 (b) Magma volume of pyroclastic- fall deposit (km3 DRE) 10.7 × 0.8 = 8.6 6.8 1.6~2.0 (b) / (a): Fall-flow ratio ∗∗∗ 0.8 0.7 Total volume of ejecta (km3 DRE) 19.3 14~16 Mean 15 3.7~5.3 Mean 4.5 Volume of disrupted volcanic edifices (km3) 3 (a.s.l.) 5 2.5 + 1.6 = 4.1 incl. dilation Volume of topog. depression (km3) 5 (b.s.l.) 0.7 1.1 Dia. of caldera surface (km) 5 2 2.5 Duration of climatic eruption (hour) 23 11~16 3~5 Max. column height (km) 50 26 34~40 Table 3. Comparison among the three caldera-forming eruptions. (∗): after Hildreth [1983] and Fierstein and Hidreth [1992]; (∗∗): after Pal- adio-Melosantos et al. [1996] and Scott et al. [1996]; (∗∗∗): (b) / (a) is assumed to be 0.8. 9 Paladio-Melosantos et al. [1996] distinguished tephra layers into A, B, C and D. C-layer is the most volumi- nous, 98% of the total volume, amounting to 1.6~2.0 km3 DRE. Scott et al. [1996] also made a comparative discussion between Pinatubo and Novarupta eruptions, and independently assumed the caldera structure of Pinatubo to be funnel-shape of a large depth (6~8 km b.s.l.) and full of lithic material: it is noticeable that the concept of “caldera deposits” was entertained in the discussions of both Novarupta and Pinatubo calderas. Here we refer to the two large eruptions to check their relationships, if any, between volume of pyro- clastic-flow deposits and pumice-fall deposits. If these are correlated, we may apply the correlation to the Krakatau eruption. In Table 3, the ratios of volume DRE of pyroclastic-fall to that of pyroclastic-flow (fall/flow ratio) in the eruptions of Novarupta and Pinatubo are roughly 0.7~0.8. Such similar ratios may be not always accidental. If the ratios are approximately common to the eruptions of similar magmas and magnitudes, we may use the ratio to estimate volume of pumice-fall de- posits of the third eruption of which volume of pyro- clastic-flow is known. We apply the ratio 0.8 to the eruption of Krakatau. Then volume of the Krakatau pyroclastic-fall is approximately estimated at 8.6 km3 DRE that is comparable to that of Novarupta (6.8 km3). The fall deposits of Krakatau may have been distributed over the wider area than Novarupta considering the dif- ference in their column heights. The total volume of the ejecta from Krakatau is revised to (10.7+8.6) km3 = 19.3 km3 DRE as shown in Table 4. In the table, original Verbeek’s estimation of the flow and fall deposits are converted to 11.8 km3 DRE and 0.4 km3 DRE, respectively. Assuming that the pres- ent re-estimation is correct, the difference is mainly due to underestimation of the pyroclastic-fall deposits. The present re-estimation of pyroclastic-fall deposits is ap- proximate but may be more reasonable than Verbeek’s estimate because the latter depended on scanty data. Verbeek [1885] commented on the errors of his esti- mation (18.2 km3 bulk in total) as “The figures every- where used for thickness can be hardly be considered too large, but is most likely too small. The possible in- accuracy of this calculation, I think, is not more than 2 or at most 3 km3.” Strictly speaking, we cannot argue the errors in Verbeek’s estimation because we cannot evaluate the accuracy of the two charts used in his cal- culations. The present result, 19.3 km3 at the minimum, seems not so unreasonable comparing with 15 km3 of the 1912 Novarupta eruption in Table 3. 3.3. Volume estimations of eruption products in the 1833 Krakatau eruption by various authors After Verbeek, volume estimations of eruption products in the 1883 Krakatau eruption have been car- ried out, sometimes as a by-product of fundamental re- searches. Here, only the results in the volume estimation of pyroclastic-flows and -falls, and disrupted parts, by various authors shall be mentioned: (1) Verbeek [1885] first carried out re-plumbing about one month after the 1883 eruption. This is a merit of his work if we consider post-depositional erosions. He did not discriminate between flows and falls, but in the present paper, these two are classified according to the distance from the eruption center (Table 1). Thus his estimation of the total volume of the ejecta amounts to 18.2 km3 bulk or 12.2 km3 DRE (Table 4). (2) Archibald [1888] estimated volume of ash falls collecting reports of the dust falling on 8 ships travers- ing the Indian Ocean assuming its thickness to be 0.005 m and the area represented by these ships to be 2.8×106 km2 (1700 km × 1700 km; cf. Figure 2b). Then he got the total amount of ash fallout as 14.4 km3 bulk, or ×0.3/2.6 = 1.7 km3 DRE. (3) Yokoyama [1981] determined volume of the total ejecta as 101 km3 in order of magnitude based on Verbeek’s estimation (18.2 km3 bulk) and in considera- tion of its ambiguity. Also he deduced volume of “caldera deposit” as 9.3 km3 bulk from gravity anomalies, and later Yokoyama [2014] revised it to 13±3 km3 bulk by new knowledge on Vp/Vs ratios [Jaxybulatov et al. 2011]. (4) Self and Rampino [1981] carried out field in- vestigation on Krakatau Islands to correlate the eruption sequence with the stratigraphy of pyroclastic deposits. They evaluated the total volume of eruption deposits as 18~21 km3 bulk or 9~10 km3 DRE, made up as fol- KRAKATAU ERUPTION PRODUCTS Verbeek [1885] This paper Bullk DRE (km3) DRE (km3) Pyroclastic-flow (2.5+ 12.1) × (2.1/2.6) = 11.8 Revised 10.7 Pyroclastic-fall 3.6 × (0.3/2.6) = 0.4 ∗Newly estimated 8.6 Sum 18.2 12.2 19.3 Table 4. Revised volume of the pyroclastic deposits on and around Krakatau Islands; (∗) analogical estimation in Section 3.2. lows: 12 km3 bulk of pyroclastic deposits and 8.5 km3 bulk of ash fall deposits. (5) Sigurdsson et al. [1991] collected SCUBA cores of the sediments in the area between Krakatau Islands and Sebesi Island, mainly around Steers and Calmeyer reefs, and newly constructed a bathymetric chart based on their 1990 survey, and a contour map of the changes in bathymetry. They followed Verbeek’s manner and deter- mined finally the total volume of the eruption products as 9.0 km3 DRE, and caldera depression as 8.9 km3 bulk. (6) Mandeville et al. [1996] studied the stratigraphy of submarine pyroclastic deposits within a 15 km radius of the caldera collecting land-based and SCUBA cores, and estimated the volume of flow deposits as 12.4 km3 DRE and fall deposits as 2.3 km3 DRE: the total vol- ume amounts to 14.7 km3 DRE. (7) In the present paper, Verbeek’s estimate of vol- ume of the flow deposits is revised to 10.7 km3 DRE fol- lowing his procedures. The volume of pyroclastic fall is assumed to be 80% of that of flow deposits on the anal- ogy of the two caldera-forming eruptions: it is 8.6 km3 DRE. Finally the total amount of pyroclastic material is obtained as 19.3 km3 DRE. And the volume of dis- rupted volcano edifice is reckoned at 8 km3 bulk. In a summary, the total volume of the ejecta from the 1883 Krakatau eruption calculated by various au- thors ranges over 9.0, 9.5, 12.2, 14.7 and 19.3 km3 DRE, as summarized above. The results are rather diverse be- cause the above authors based on the different assump- tions and the different data. After all, it may be concluded that the total volume of the Krakatau pyroclastic prod- ucts including the flows and falls, widely ranges 9~19 km3 DRE, being roughly comparable to Novarupta and larger than Pinatubo, and that the volume of disrupted volcano edifice is roughly 8~9 km3 bulk according to some authors. Now we may say that we must be con- tent with the volume estimation with accuracy of order of magnitude in the discussion of ejecta volume from the 1883 Krakatau eruption. In no distant future, we shall be able to estimate the volume of volcanic ejecta as mass flux at the level of volcano vents by remote sensing techniques in place of the present manner of searching the ejecta around volcanoes after eruptions. 4. Subsurface structure of Krakatau caldera and set- tlement of the eruptive products 4.1. Caldera deposits of Krakatau caldera Generally caldera depression is defined as a vacant space formed by disruption of volcano edifices during violent explosions of pyroclastic ejecta. In the case of Krakatau caldera, the caldera depression was formed by disruption of 2/3 parts of pre-existed Krakatau Is- land and excavation of the seabed to a depth of 250 m. The present author [1981, 2014] discussed subsurface structure of the caldera by the gravity anomalies ob- served there, and assumed “caldera deposits” of low density immediately beneath the caldera depression, and estimated their volume at roughly 13 km3 bulk. Van Bemmelen [1929, 1949] expressed his concep- tion that disrupted volcano edifices collapsed into magma reservoir in replacement of the discharged material by magmatic eruptions. To accept this idea, the following two assumptions should be satisfied: (1) the volume of disrupted volcano edifice (“caldera depression”) should be approximately equal to that of discharged magma, and (2) lithics of the volcano edifice scarcely remain at the surface around the caldera: in the following, we dis- cuss the problems on Krakatau caldera. Firstly, to verify assumption (1), both the volume should be similar in the exact numerals, not in the order of magnitude. The total volume of pyroclastic mate- rial of the Krakatau eruption is estimated at 19.3 km3 DRE in Table 4, and that of disrupted volcano edifice is reckoned as (2.9+5.1) = 8 km3 in Section 2.3: the ratio of (volume of ejecta) / (volume of disrupted volcanic edifice) is (19.3 / 8). Here, we check the cases of No- varupta and Pinatubo calderas: in Table 3, the ratios are (15/5.7) and (4.5/5.2), respectively. In short, Van Bem- melen’s idea is narrowly valid in the Pinatubo eruption though clearly invalid in the Krakatau and the No- varupta eruptions. Hence we may conclude that the as- sumption proves not always to hold good. Secondly, as for assumption (2), Verbeek [1885] re- ported that the ejecta mainly consisted of newly erupted magmatic material and fragments of the old volcano edifices, and that the latter volume amounted to 5% of the total. His estimation 5% may have been very rough admitting that the estimation method has not been es- tablished even now. Stehn [1929] observed a lithic part of 16 m thick of 40 m columnar section at the NW part of Lang Island. We can find other higher percentages in the recent literatures: 13 % at the shore of Lang Is- land [Mandeville et al. 1996]. Simkin and Fiske [1983, p. 135, fig. 32] introduced the investigation report carried out by U.S. steamship Juniata around Krakatau Islands within 2 weeks of the catastrophic eruption: the crew prepared the first post-paroxysm map of Krakatau, dated September 10, 1883, in which we see several small islands formed of the fragmental deposits near Verlaten Island. Verbeek started his survey at the end of Sep- tember and did not register such large amount of de- posits. The delay of 2 weeks in the start of the surveys may be rather significant from the viewpoint of after- depositional erosions. Hence, the present author doubts validity of 5 % reported by Verbeek [1885]. YOKOYAMA 10 11 We may assume the caldera deposits to be com- posed of juvenile pyroclastic ejecta and crushed pre- caldera rocks (lithic material): all these are fallbacks, and the latter may have been due to explosivity of the caldera-forming eruption. Admitting the volume equiv- alence between the caldera depression and the total erupted magmas, it is necessary to locate the magma reservoir and to explain the mechanism how the reser- voir cavity was filled with the disrupted material. The validity of Van Bemmelen’s assumptions should de- pend on the possibility of settlement of the eruption products into magma reservoirs. 4.2. Magma reservoir beneath Krakatau volcano and mechanical stability of a cavity in the earth crust Williams [1941] supported Van Bemmelen’s idea and established the 1883 Krakatau eruption as the typ- ical example in collapses of a pre-existing volcanic edi- fice that was disrupted by large-scale eruption of silicic magma derived from a shallow magma reservoir. “Shallow magma reservoirs” may be a rather diffi- cult conception in the 21st century except “magma stops” between magma reservoirs and vents at some basaltic volcanoes. Beneath Krakatau caldera, searching for liq- uid zones, Harujono et al. [1989] detected several zones of strong attenuation of seismic S-waves and Jaxybula- tov et al. [2011] determined the distribution of high Vp/Vs ratio of seismic waves; Yokoyama [2014] as- sumed a probable magma reservoir at a zone of about 10 km deep, characterized commonly by S-wave atten- uation and high Vp/Vs ratio as shown in Figure 7. Here, mechanical stability of a magma reservoir with a cavity produced by magma ejection shall be ex- amined: a cavity in the deep crust cannot exist against lithostatic pressure. If we assume compressive strength of the shallow crust to be 100 MPa, a cavity cannot exist at depths deeper than about 3 km in average. In other words, a magma reservoir at a depth deeper than 3 km should be crushed according as the magma is discharged from the reservoir. It is doubtful how the disrupted vol- canic edifices can reach such magma reservoirs during formation of calderas. In the present case, the cataclysm is assumed to start from a magma reservoir, not from a surface volcano. When a magma reservoir discharges magma partly or totally through a conduit leading to the volcano, the reservoir should be crushed partly or totally under crustal load when it is located at deeper than 3 km. Here, for the sake of simplicity, we assume a spherical magma reservoir at a depth of 10 km, and concentric crushing force F to act at the center of the reservoir. Timoshenko and Goodier [1970] discussed “Force at a point in an in- finite solid”, and their results are approximately appli- cable to the present problem. Beginning from Laplace’s equation, they finally gave vertical stress v acting along the axis r = 0 of cylindrical co-ordinates (r, i, z) where the origin is the center of the cavity, as (1) Assuming Poisson’s ratio o = 0.25, we get: vz,r=0 = 0.186 F/z 2 (2) Qualitatively speaking, the reservoir should be crushed by lithostatic pressure according as the magma is discharged, and the effect of crushes of the reservoir (crushing force F) decreases rapidly in proportion to z2. In other words, the crushing effect does not reach the sur- face when the reservoir is located at a sufficient depth. If a magma reservoir were located at shallow depths, shallower than 5 km, the crustal strength is much weaker than the deeper parts and cover rocks of the reservoir may be easily exploded. On Krakatau caldera, Yokoyama [2014] assumed depth of the caldera deposits to be about 2 km (Figure 7) and hence a shallow reservoir should be rather near to the caldera bottom. Anyway disrupted volcano edifices cannot settle at a magma reservoir in replacement of the discharged magma but may settle as caldera deposits. 8 1 1 2 3F z,z r 0 2v r o o= - - += ^ ^ h h" , KRAKATAU ERUPTION PRODUCTS Figure 7. Assumptive caldera deposits and magma reservoir pro- jected on a vertical plane along B—B’ line in Figure 6. Dotted part schematically shows caldera deposits. Distribution of the S-wave at- tenuations (striped parts) was obtained by Harjono et al. [1989] and that of the Vp/Vs ratios by Jaxybulatov et al. [2011]. A colored cir- cle on the assumed magma reservoir shows the total amount of the pyroclastic ejecta (19 km3 DRE) discharged by the 1883 eruption. 4.3. A model of Krakatau caldera viewed from its sub- surface structure Williams [1941] discussed general features of calderas and their origin but did not refer to their sub- surface structure. Thus, he had no idea about caldera de- posits. Yokoyama [1963] studied subsurface structure of some calderas and classified them into two: high-grav- ity- anomaly type and low-gravity-anomaly type. The gravity anomalies are due to density of the caldera de- posits which have been verified by drillings at some calderas on land in Japan and Mexico [Yokoyama 2014]. At Krakatau caldera, Yokoyama [1981] analyzed low gravity anomalies observed on the islands and later Deplus et al. [1995] discussed the inner structure of the caldera analyzing gravity anomalies observed in a larger area around Krakatau Islands. Yokoyama [2014] inter- preted the gravity anomalies to be due to mainly the caldera deposits of low density and estimated their vol- ume at 13 ± 3 km3. The Krakatau caldera depression (5.1 km3) is about 200 m deep b.s.l. and full of seawater. Beneath the caldera bottom, there exist caldera deposits that are composed of fallbacks of exploded materials (lithics and juveniles) and are soaked in seawater. The watery deposits are verified by anomalous values of Vp/Vs ratio as high as 2.2 roughly around Krakatau caldera determined by Jaxybulatov et al. [2011]. On the other hand, Harjono et al. [1989] detected several seis- mic S-wave attenuating zones beneath the Krakatau islands. A vertical configuration of Krakatau caldera ob- tained by Yokoyama [2014] is reproduced in Figure 7 where the most probable zone of magma reservoir must be endowed with both the seismic-wave charac- ters in common. This model has been achieved by con- tributions from several authors as mentioned above. The total volume of the ejecta discharged by the 1883 eruption (19.3 km3 DRE) is projected as a sphere on the assumed magma reservoir. 5. Concluding remarks The estimation of volume of pyroclatic ejecta from the 1883 eruption of Krakatau originally carried out by Verbeek [1885] is revised in this paper principally based on his data. As a result, his estimation of the flow deposits proved to be acceptable while the fall deposits to be much underestimated. Its revision in the present paper is based on the assumption that the volume ratio, fall/flow is similar among caldera-forming eruptions, deduced from the 1912 eruption of Novarupta and the 1991 eruption of Pinatubo. As a result, the volume of pyroclastic fall is evaluated to be similar to that of the flow deposits. Originally volume estimation of volcanic ejecta with reasonable accuracy is very difficult, especially on island volcanoes. In the case of the 1883 Krakatau eruption, we must be content to get results with accuracy of the order of magnitude, and so Van Bemmelen’s assump- tion that the dense rock equivalent of the total ejecta derived from a caldera-forming eruption is the same as the volume of the resulting caldera depression, proved not necessarily to hold good. As for the final settlement of the disrupted volcanic materials, some papers have proposed that the materials settled in magma reservoirs in replacement of the discharged magma by reason of volumetric equality between both the materials. How- ever, the equality is not always approved as discussed in the present paper and is not the necessary condition for caldera formations. On the other hand, the caldera de- posits of Krakatau caldera are confirmed as well as some other calderas, and its magma reservoir is con- jectured to be located at a depth of about 10 km from a geophysical standpoint. The paucity of lithic frag- ments among the Krakatau ejecta can be semi-quanti- tatively attributed to their settlement into the caldera deposits. The eruption products of the 1883 Krakatau eruption are discussed on their origin, volume and set- tlement on the base of the two examples of caldera- forming eruptions that occurred under our eyes. If the existence and the origin of “caldera deposits” are approved, we may define “caldera” of low-gravity- anomaly type, from the viewpoint of subsurface struc- ture as a large crater, such as larger than 2 km in diameter, with flat bottom, underlain by caldera deposits and ac- companied with a large quantity of pyroclastic deposits around the volcano. Acknowledgements. The above discussion is indebted to the data originally provided by Verbeek [1885]. The present author had a chance to read this priceless book at the Volcanological Survey of Indonesia, Bandung, and thanks Dr. S. Siswowidjoyo and Prof. A. Sudradjat for their kind cooperation. The author’s knowledge on plumbing in old times owes to Dr. M. Tsuchide of the Hydrographic Department of Japan, to whom the author is grateful. Critical and helpful comments by anonymous referees are gratefully appreciated. References Archibald, E.D. (1888). General analysis of the connec- tion between the unusual meteorological phenom- ena of 1883-6, and the eruptions of Krakatoa in May and August, 1883, In: G.J. Symons (ed.), The eruption of Krakatoa, and subsequent phenomena, 426-463. Bemmelen, R.W., van (1929). Het Caldera probleem, De Mijningenieur, 10, 101-112. Bemmelen, R.W., van (1949). Geology of Indonesia, Chapter 3: Volcanism, A. Active volcanism, The Hague: Govt. Printing Office, 188-224. Camus, G., M. Gloaguen, A. Provost and P.M. Vincent (1991). Emplacement of a debris avalanche during the 1883 eruption of Krakatau (Sunda straits, In- YOKOYAMA 12 13 donesia), XVII Pacific Sci. Cong. (Honolulu, May 27 - June 2, 1991), 3 pp. and 6 figs. Deplus, C.S., S. Bonvalot, D. Dahrin, M. Diament, H. Harjono and J. Dubois (1995). Inner structure of the Krakatau volcanic complex (Indonesia) from grav- ity and bathymetry data, J. Volcanol. Geotherm. Res., 64, 23-52. Fierstein, J., and W. Hildreth (1992). The Plinian erup- tions of 1912 at Novarupta, Katmai National Park, Alaska, B. Volcanol., 54, 646-684. Harjono, H., M. Diament, L. Nouaili and J. Dubois (1989). Detection of magma bodies beneath Krakatau volcano (Indonesia) from anomalous shear waves, J. Volcanol. Geotherm. Res., 39, 335-348. Hildreth, W. (1983). The compositionally zoned erup- tion of 1912 in the Valley of Ten Thousand Smokes, Katmai National Park, Alaska, J. Volcanol. Geot- herm. Res., 18, 1-56. Jaxybulatov, K., I. Koulakov, M.I. Seht, K. Klinge, C. Re- ichert, B. Dahren and V.R. Troll (2011). Evidence for high fluid/melt content beneath Krakatau volcano (Indonesia) from local earthquake tomography, J. Volcanol. Geotherm. Res., 206, 96-105. Judd, J.W. (1888). On the volcanic phenomena of the eruption, and distribution of the ejected materials, In: G.J. Symons (ed.), The eruption of Krakatoa, and subsequent phenomena, 1-56. Mandeville, C.W., S. Carey and H. Sigurdosson (1996). Sedimentology of the Krakatau 1883 submarine py- roclastic deposits, B. Volcanol., 57, 512-529. Paladio-Melosantos, M.L.O., R.U. Solidum, W.E. Scott, R.B. Quiambao, J.V. Umbal, K.S. Rodolfo, B.S. Tubianosa, P.J. Delosreyes, R.A. Alonso and H.B. Ruelo (1996). Tephra falls of the 1991 eruptions of Mount Pinatubo, In: C.G. Newhall and R.S. Punong- bayan (eds.), Fire and mud, eruptions and lahars of Mount Pinatubo, Philippines, 513-535. Sarna-Woicicki, A.M., S. Shipley, R.B. Waitt Jr., D. Dzurisin and S.H. Wood (1981). Areal distribution, thickness, mass, volume, and grain size of air-fall ash from the six major eruptions of 1980, In: The 1980 eruptions of Mount St. Helens, Washington, Geo- logical Survey Professional Paper, 1250, 577-600. Scott, W.E., R.P. Hoblitt, R.C. Torres, S. Self, M.M.L. Martinez and T. Nillos Jr. (1996). Pyroclastic flows of the June 15, 1991, climactic eruption of Mount Pinatubo, In: C.G. Newhall and R.S. Punongbayan (eds.), Fire and mud, eruptions and lahars of Mount Pinatubo, Philippines, 545-570. Self, S., and M.R. Rampino (1981). The 1883 eruption of Krakatau, Nature, 294, 699-704. Sigurdsson, H., and S. Carey (1989) Plinian and co-ign- imbrite tephra fall from the 1815 eruption of Tamb- ora volcano, B. Volcanol., 51, 242-270. Sigurdsson, H., S. Carey, C. Mandeville and S. Bronto (1991). Pyroclastic flows of the1883 Krakatau erup- tion, EOS, 72, 377-392. Simkin, T., and R.S. Fiskes (1983). Krakatau 1883: the volcanic eruption and its effects, Smithsonian Insti- tution Press, Washington D.C., 464 pp. Stehn, Ch.E. (1929). The geology and volcanism of the Krakatau group. Krakatau: Fourth Pacific Science Congress, Part I, 1-55. Timoshenko, S.P., and J.N. Goodier (1970). Theory of elasticity, McGraw-Hill Book Co., New York, 567 pp. Verbeek, R.D.M. (1885). Krakatau, Landsrukkerij, Bata- via, 546 pp. Wharton, W.J.L. (1888). On the seismic sea waves caused by the eruption of Krakatoa, August 26th and 27th, 1883, In: G.J. Symons (ed.), The eruption of Kraka- toa, and subsequent phenomena, 89-150. Williams, H. (1941). Calderas and their origin, Univ. California Pub., Bulletin of the Department of Ge- ological Sciences, 25, University of California Press, 239-346. Yokoyama, I. (1963). Structure of calderas and gravity anomaly, B. Volcanol., 26, 67-72. Yokoyama, I. (1981). A geophysical interpretation of the 1883 Krakatau eruption, J. Volcanol. Geotherm. Res., 9, 359-378. Yokoyama, I. (2014). Krakatau caldera deposits: revis- ited and verification by geophysical means, Annals of Geophysics, 57 (5), S0541. Corresponding author: Izumi Yokoyama, The Japan Academy, Ueno Park, Tokyo, Japan; email: iyokoyama@joy.ocn.ne.jp. © 2015 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. KRAKATAU ERUPTION PRODUCTS << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice