ISSN 0916-7390

九州大学大学院理学研究院紀要

M E M O I R S
OF THEFACULTYOFSCIENCESKYUSHU UNIVERSITY

Series D Earth and Planetary Sciences
VOLUME XXXⅤ
No. 1

FUKUOKA, JAPANMarch, 2021

Contents

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo:
Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia ………………

Editorial Board
Atsushi Toramaru, Chief Editor
Takeshi Ikeda

All communications relating to this Memoirs should be
addressed to the Dean of the Faculty of Sciences, Kyushu University
Motooka, Fukuoka 819-0395, JAPAN1

Mem. Fac. Sci., Kyushu Univ., Ser. D, Earth ＆ Planet. Sci., Vol. XXXV, No. 1,
pp. 1-17, March 26, 2021

Eruption Types and Conduit Dynamics of Kukusan and Genteng Volcanoes of
the Ijen Volcanic Complex, Indonesia
Takeru Mitsuoka1, Atsushi Toramaru1, Agung Harijoko2, and Haryo Edi Wibowo2

Abstract
Eruption mechanisms of mafic magma are strongly debated, especially the mafic Plinian eruption. Few
examples of mafic Plinian eruptions prevent us from knowing what condition causes rapid magma ascentand degassing processes. Kukusan and Genteng volcanoes, which are monogenetic volcanoes in the Ijen
caldera complex in Indonesia, produced explosive eruptions of mafic magma. Their chemical compositions
are similar, but the volcanoes are believed to have different eruption types. Therefore, these volcanoes werechosen to reveal the factors affecting the eruption types of mafic magma. Outcrops were observed at 31
locations in fieldwork. Geochemistry, petrography, and microtextural analyses of scoria were conducted.
Based on the dispersal area and fragmentation degree of the volcanic deposits, the eruption type of Kukusan
is considered as sub-Plinian, while that of Genteng is considered as Strombolian. The microtextural analysis
showed that the Kukusan scoria had a higher bubble number density (BND), microlite number density(MND), and microlite crystallinity than the Genteng scoria. High BND and MND values suggest rapid
decompression, i.e. rapid magma ascent in the volcanic conduit. High microlite crystallinity indicates highmagma viscosity. The higher ascent rate and higher viscosity of the Kukusan magma suggest coupleddegassing and fragmentation, leading to a more explosive sub-Plinian eruption. In contrast, the lower ascent
rate and lower viscosity of the Genteng magma suggests decoupled degassing, which leads to a Strombolian
eruption.
Keywords: ascent rate; degassing process; scoria; bubble number density; microlite crystallinity

Ⅰ. Introduction
Explosive volcanic eruptions are classified by the dispersal area and fragmentation degree of pyroclastic
deposits (Walker 1973). Generally, mafic magma generates Hawaiian or Strombolian eruptions, and rarely Plinian
eruptions (Kereszturi and Nemeth 2012). However, it has been recently recognized that mafic magma sometimes
produces a Plinian eruption. Some examples of mafic Plinian eruptions include Etna 122 BC in Italy (Sable et al.
2006), Tarawera 1886 in New Zealand (Houghton et al. 2004), and Pleistocene Fontana Lapilli in Nicaragua(Costantini et al. 2010). Determining the mechanisms that cause various types of eruptions is critical for thehazard mitigation of highly explosive eruptions. Three mechanisms are thought to control mafic Plinian eruptions:rapid magma ascent, coupled degassing, and brittle fragmentation (Sable et al. 2006, Houghton and Gonnermann2008, Costantini et al. 2010). Weakly explosive style of mafic magma is believed to be associated with decoupleddegassing, in which bubbles are able to migrate upward within relatively slow ascending melt, whereas highly
explosive style of silicic or mafic magma is believed to be associated with coupled degassing, in which bubbles
are immobile relative to the melt (Sable et al. 2006). However, few examples of mafic Plinian eruptions prevent
us from knowing what condition causes rapid magma ascent and coupled degassing.
Kukusan and Genteng volcanoes form part of the Ijen Caldera Complex in East Java, Indonesia. Kukusan and
Genteng produced explosive eruptions of mafic magma at between 21 and 2 ka (Sitorus 1990, Sundhoro 1990,
1 Department of Earth and Planetary Sciences, Graduate School of Science, Kyushu University, Japan
2 Department of Geological Engineering, Faculty of Engineering, Universitas Gadjah Mada, Indonesia

2

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

Suhendro et al. 2016). Although some researchers studied the Ijen Caldera Complex, they did not map thedispersal area nor conduct a textural analysis of pyroclastic deposits. Drawing a dispersal map of volcanic
deposits will contribute to volcanic hazard assessment. In this study, we aim to determine the eruption types ofKukusan and Genteng volcanoes. Kukusan and Genteng volcanoes are located in the same caldera and have
similar chemical compositions. The eruption types of Kukusan and Genteng are unknown but are thought to bedifferent based on geological map (Sujanto et al. 1988). Therefore, the pyroclastic deposits of Kukusan and
Genteng are herein used to compare their magma properties, which likely affect the eruption types. The textural
analysis of bubbles and crystals provides information on magma ascent, degassing, and fragmentation in a
volcanic conduit. The study of the conduit dynamics of Kukusan can add a case example of mafic Plinian
eruption, as well as describe the conditions that cause rapid magma ascent, coupled/decoupled degassing, andfragmentation. In this study, we aim to understand the factors that have an effect on their explosive eruption ofmafic magma by comparing the cases of Kukusan and Genteng volcanoes.
Section 2 presents the geological settings of the Ijen Caldera Complex. In Section 3, we outline the methods
utilized in this study. Section 4 reports the fieldwork results and provides the laboratory analyses of the Kukusan
and Genteng pyroclastic deposits. Finally, we interpret the obtained results and discuss the relationship between
conduit dynamics and the Kukusan and Genteng eruption types.

Ⅱ. Geological settings
Java Island is located in the subduction of the Indian-Australian plate beneath the Eurasian plate (Hamilton1979). The Ijen Volcanic Complex is located in the easternmost part of Java Island (Fig. 1) within the Quaternary
volcanic front of the Sunda arc. The volcanic complex features a caldera rim on the northern boundary, five
stratovolcanoes on the southern boundary, and several post-caldera monogenetic volcanoes within the intracaldera plane (Caudron et al. 2015). The Ijen caldera’s diameter is 18 km in the east-west direction and 15 km in
the north-south direction. The caldera formation occurred between 294 and 50 ka (Sitorus 1990). The post-caldera
eruption formed over 22 separate vents. The vents are separated geographically into two groups: the caldera rim
group and the intra-caldera group (Handley et al. 2007). Post-caldera volcanic deposits represent various rocktypes: lava, pyroclastic flow, pyroclastic air fall and surge deposits (Sitorus 1990). Kukusan and Genteng aremonogenetic volcanoes located inside the Ijen caldera. Kukusan is approximately 4 km away from Genteng.
Based on the classification by Handley et al. (2007), Kukusan and Genteng are the intra-caldera volcanoes, which
erupted between 21 and 2 ka (Sitorus 1990).

Fig. 1 East Java and Bali regional map. The research area is demarcated by the red box and includes the Ijen Caldera Complex and Mt. Ijen.

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

3

Ⅲ. Methods
Fieldwork was conducted two times in 2018 and once in 2019. Pyroclastic stratigraphy was recognized at 31
locations. The observed units were correlated to draw an isopach map of the pyroclastic deposits. The samples
were collected from six locations. The bulk chemical compositions of the major elements were analyzed using a
Primus IV X-ray fluorescence spectrometer (Rigaku Corporation). In slightly weathered rock samples, theweathered part was removed before the analysis. Grain size distribution (GSD) was measured in the laboratory ofPetrology and Volcanology, Kyushu University with sieves ranging from 32 to 1/2 mm (øGSD from -5 to 1 on theøGSD scale of Walker 1971). We manually selected the 8-16 mm juvenile clasts for the component analysis. The
color, morphology, vesicularity, and bubble size of juveniles were observed by the naked eye. Several 8-16 mm
scoria clasts were randomly selected for bulk density analysis. The scoria volume was measured using a 3-D
scanner (Solutionix D700), whereas an electric scale (Shimadzu UW420S) was used to estimate the scoria
weight. Bulk density of each scoria was calculated using its weight and volume. A few representative clasts fromeach volcano were selected: maximum density, minimum density, and modal density. Then, several polished thinsections were produced from the clasts for the petrographic and microtextural analyses.
The microscopic textures were described using electron microscope images. The backscattered electron (BSE)
images were acquired at three magnifications: 30×, 100×, and 1000×. The 30× and 100× images were captured
using a TM3030Plus scanning electron microscope (SEM) with an accelerating voltage of 15 kV. The 1000×
images were taken using a JEOL JXA-8530F field emission electron probe micro analyzer (FE-EPMA) with aprove current of 6 nA and an accelerating voltage of 15 kV. The outlines of the crystals and bubbles were traced
manually using Adobe Illustrator (version 24.3). The areas of crystal and bubble were measured using Image J(version 1.52k). Coalesced bubbles were separated into individual bubbles by manual tracing (Fig. 2).
Images at a magnification of 30× were used to observe phenocryst assemblages, vesicularity, and bubble-freephenocryst crystallinity. We took 10-16 images per thin section, and thereby made a large mosaic image using
Adobe Illustrator. Crystals with an area over 0.01 mm2 were considered phenocryst. The fractions of crystal
phase, vesicularity, and phenocryst crystallinity in 2-D were calculated using the following equations:

(1)

(2)

(3)

where image area is the width × length of the BSE image.
Images at a magnification of 100× were used to calculate mean bubble size and bubble number density(BND). Two images of each scoria were taken in different parts of a thin section. The mean bubble size and BNDof each scoria were obtained by averaging the measurements from these two images. BND in 2-D (NAbub) wasconverted to 3-D (NVbub) using the following formula (Toramaru 2019):(4)

where d is the mean equivalent bubble diameter, and Øbub is the vesicularity in 2-D image. The correction of1-Øbub was applied to represent magma conditions before bubble formation in the ascending melt. Although somestudies have provided improved formulas to solve the stereological problems (Sahagian and Proussevitch 1998,
Higgins 2000, 2002), we favored a much simpler one.
Images at a magnification of 1000× were used to calculate the mean microlite size, microlite number density

4

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo(MND), microlite crystallinity, and microlite aspect ratio. Thick sections of each scoria were made forobservation with EPMA. Two images of each scoria were taken in different parts. The microlite quantitative dataof each scoria were obtained by averaging the measurements from these two images.

Fig. 2 Part of a 100× BSE image showing an example of coalescence bubbles. Such bubbles were manually traced.

Ⅳ. Results
4.1. Stratigraphy of pyroclastic deposits and isopach maps
The pyroclastic deposits of Kukusan volcano were observed at 30 locations (Fig. 3a). Most of these locations
are roadcut outcrops represented by scoria, lithic material, and ash. At location I05, the stratigraphic successions(from bottom to top) consist of massive black coarse ash, white fine ash, and two layers of scoria with lithics (Fig.
4a). Both scoria layers are black to brown, loose-packing, and well-sorted. The white fine ash is widely dispersedin the study area. Due to its distinct color (white), we used it as a key layer to correlate stratigraphic columns. The
white ash is a product of the Plinian eruption; however, it was derived from an unknown sourced outside of the
Ijen volcanic complex (Sundhoro 1990). Similar scoria and white ash layers were observed at locations I13, I17(Fig. 4b), and I15. The scoria layer thickness increases closer to the Kukusan vent: scoria layer thickness equals
46, 65, 105, and 120 cm at locations I05, I13, I17, and I15, respectively. The scoria layer was also observed inside
the Kukusan crater at location I89. The Kukusan scoria dispersed not only to the northern side of the vent, butalso to the eastern and north-eastern sides. At locations I73 and I99, the scoria layer is found between two layersof white fine ash (Fig. 4c). The maximum size of the Kukusan scoria was 10 cm at location I13. At somelocations, the scoria layer contains cauliflower bombs. The observed scoria clasts are bubble-rich. These bubbles
are small (<2 mm), and not elongated. The lithics are red-orange in color. The maximum size of lithics reaches
5.6 cm at location I18. The Kukusan scoria samples were collected at locations I13, I15, I17, I89, and I99. At
these locations, the scoria layer was subdivided into two layers (A and C) from bottom to top (Fig. 4a).
The pyroclastic deposits of Genteng volcano were observed only at one location. The outcrop is located on theflank of the Genteng vent at location I43. The stratigraphic successions (140 cm thick) consists of two layers ofscoria with rare lithics (Fig. 4d). The scoria layers are black to brown, loose-packing, and well-sorted. The
maximum Genteng scoria size is 23 cm. Some scoria clasts are slightly weathered on their surfaces. Most scoria
clasts are bubble-rich. These bubbles are relatively larger (up to 1 cm), and not elongated. The scoria layer was
subdivided into two layers (A and B) from bottom to top (Fig. 4d). Layer A is a mixture of ash and scoria, whilelayer B consists scoria only.
The dispersal area of Kukusan scoria is wide compared with that of Genteng (refer to isopach map, Fig. 3a).
The Genteng scoria can be found solely on the flank of Genteng volcano. Based on the geomorphology andgeological map (Sujanto et al. 1988), we drew a dispersal map for the Genteng scoria (Fig. 3b).

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

5Fig. 3 (a) Outcrop location and dispersal map of the Kukusan scoria. (b) Outcrop location and dispersal map of the Genteng scoria.
Black squares—outcrop locations; the number above the black square—location number; the number below the black square—thickness of scoria layer [cm]; black triangles—volcano vent locations; solid line—an isopach; dashed line—an estimated isopach.

Fig. 4 Photographys of outcrops with volcanic deposits. The length of a measurement stick is 1 m. (a) Kukusan volcanic deposits at
location I05. Two scoria layers are underlain by white ash. (b) Kukusan volcanic deposits at location I17. (c) Kukusan volcanic
deposits at location I99. (d) Genteng volcanic deposits at location I43.

6

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

4.2. Whole rock geochemistry
The SiO2 content of the Kukusan scoria (52.58 ± 0.5 wt.%) is slightly higher than that published previously(50.23 wt.%) by Sundhoro (1990). Similarly, the SiO2 content of the Genteng scoria (54.89 ± 0.5 wt.%) is slightly
higher than that published previously (53.35 wt.%) by Suhendro et al. (2016). Both the Kukusan and the Genteng
scoriae are classified as basaltic andesites (Table 1). Although the chemical composition of the Kukusan magmaslightly differs in the SiO2, Fe2O3, and MgO contents from that of Genteng, the overall chemical compositions of
the two volcanoes are similar.

Table 1 Bulk chemical compositions of the Kukusan and Genteng scoriae. Unit: wt.%.

4.3. Grain size
The grain size distribution (GSD) ranges from -6 to 2 on the øGSD scale in most scoria samples. The Kukusan
GSD peaks at -4 to -2 on the øGSD scale, while Genteng GSD peaks at -6 to -4 on the øGSD scale (Fig. 5a). The
deviation (σ, defined by Walker 1971) among the Kukusan clasts (σ = 0.9-2.0) is lower than that of Genteng clasts(σ = 2.4), meaning that the Kukusan clasts are smaller and better sorted than those of Genteng.

4.4 Component analysis
The Kukusan and Genteng pyroclastic deposits are represented by scoria and lithics. The Kukusan scoria is
typically black in color, blocky and ragged in morphology, and is highly vesicular. The bubble size ranges fromsmall (several millimeters) to tiny (cannot be seen with the naked eyes). The Genteng scoria is black to brown incolor (appears brown on the surface due to weathering), blocky and ragged in morphology, and is highly
vesicular. The bubble size is predominantly small (several millimeters). Thus, the Kukusan and Genteng scoriae
can be distinguished from each other based on differences in dominant bubble size. The Kukusan lithic fragments
are likely to wall rocks of the volcanic conduit. The Kukusan scoria layer contains 2%-14 % of lithic fragments(Fig. 5b) compared with ~1 % in the Genteng scoria layer.

4.5. Bulk density
The Kukusan scoria at location I89A shows a bimodal density distribution (Fig. 5c). Its maximum density andminimum density are 0.96 and 1.77 g/cm3, respectively. The Kukusan density distribution peaks at 1.0-1.1 g/cm3
and 1.4-1.5 g/cm3. On the contrary, the Genteng scoria at location I43B presents a unimodal density distribution(Fig. 5c). Its maximum density and minimum density are 0.95 and 2.06 g/cm3, respectively. The Genteng density
distribution peaks at 1.2-1.3 g/cm3.

4.6. Vesicularity and crystallinity
The scoriae of the maximum density, minimum density, and modal density, which was a peak of density
distribution, were used for vesicularity and crystallinity observations (termed vesicular, dense, and mode,
respectively). Because of the bimodal density distribution of the Kukusan scoria, the scoria corresponding to thesmaller density peak was termed modeV, whereas that corresponding to the larger density peak was termed
modeD (Fig. 5c).
The Kukusan and Genteng scoria vesicularities range 40-70% and 21-64%, respectively (Table 2). The
phenocryst crystallinities of the Kukusan scoria and Genteng scoriae are low and constitute less than 6% and10%, respectively.

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

7

Fig. 5 (a) Grain size distribution of pyroclastic deposits. Md and σ stand for the median diameter and deviation, respectively (Walker
1971). (b) Ratio of scoria and lithic materials in the 8-16 mm clasts. The number above indicates how many samples wereanalyzed in each sampling location. (c) Bulk density distribution of scoria. The number above indicates how many samples wereanalyzed in each sampling location.

Table 2 Density, vesicularity, and phenocryst assemblages of the Kukusan and Genteng scoriae. Density classification: vesicular(minimum density), modeV (smaller density peak of the density distribution), modeD (larger density peak of the density
distribution), and dense (maximum density).

8

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

4.7. Phenocryst assemblage
Phenocrysts in the Kukusan scoria are composed of plagioclase (Pl), olivine (Ol), and iron-titanium oxide(Ox) (Fig. 6a, Fig. 6b). Pl is subhedral and commonly shows albite twin. Ol is euhedral or anhedral, whereas Ox
is euhedral or subhedral. Phenocrysts in the Genteng scoria are composed of Pl, clinopyroxene (Cpx),
orthopyroxene (Opx), Ol, and Ox (Fig. 6c, Fig. 6d). Pl is subhedral or euhedral and commonly shows albite twin.
Cpx, Opx, and Ol are euhedral or anhedral, whereas Ox is euhedral or subhedral. The Pl fraction in the total
phenocrysts in the Kukusan scoria (45%-79%) is lower than that in the Genteng scoria (83%-95%) (Table 2). On
the contrary, the Ol fraction in the total phenocrysts in the Kukusan scoria (17%-46%) is higher than that in the
Genteng scoria (0%-4%).

Fig. 6 Petrographic images of scoria. The black bar is 1 mm. Note the phenocrysts of plagioclase (Pl), olivine (Ol), orthopyroxene(Opx), and Fe-Ti oxide. (a) Petrographic image of the Kukusan scoria under open nicol (location I89A). Multiple bubbles arevisible. (b) Petrographic image of the Kukusan scoria under crossed nicol. Phenocrysts and micro-phenocryst of Pl are visible. (c)Petrographic image of the Genteng scoria under open nicol (location I43B). (d) Petrographic image of the Genteng scoria undercrossed nicol.

4.8. Bubble texture
Based on the qualitative observation of bubbles, the Kukusan and Genteng scoriae are characterized by amicrovesicular texture. Most bubble shapes are rounded and connected, indicating bubble coalescence by someextent. The bubbles are not elongated in any direction. The Kukusan scoria contains smaller bubbles and higher
BND than the Genteng scoria (Fig. 7). The mean bubble size of the Kukusan scoria (0.034-0.041 mm) is slightly
smaller than that of the Genteng scoria (0.049-0.061 mm). The BND of the Kukusan scoria (1.1 × 1013-3.4 × 1013m-3) is one order of magnitude higher than that of the Genteng scoria (3.2 × 1012-4.5 × 1012 m-3) (Table 3).

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

9Fig. 7 (a) Scanned thin section of the Kukusan modeD scoria. The red square indicates the approximate location of (c). (b) Scanned
thin section of the Genteng mode scoria. The red square indicates the approximate location of (d). (c) 100× BSE image of bubbles
in the Kukusan scoria. Black—bubble, gray—glass. Mean bubble size and BND constitute 0.035 mm and 1.9 × 1013 m-3,
respectively. (d) 100× BSE image of bubbles in the Genteng scoria. Mean bubble size and BND constitute 0.052 mm and 4.5 ×
1012 m-3, respectively.

Table 3 Mean bubble size, bubble number density, mean microlite size, microlite number density, microlite crystallinity, andmicrolite aspect ratio of the Kukusan and Genteng scoriae.

4.9. Microlite texture
All studied samples contain abundant microlites. The microlite shows various shapes (Fig. 8). Microlites in
the Kukusan scoria are composed of Pl, Px, and Ox. Pl is a tabular or acicular shape; Px is a skeletal or tabular
shape; and Ox is an equant shape. Microlites in the Genteng scoria are similarly composed of Pl, Px, and Ox. Pl
and Ox have the same shapes as in the Kukusan scoria. Px is octagon or tabular. The mean microlite size of the
Kukusan scoria (3.7-3.8 µm) is smaller than that of the Genteng scoria (5.0-7.9 µm) (Table 3). The MND of the
Kukusan scoria (7.9 × 1010 m-2-9.7 × 1010 m-2) is higher than that of the Genteng scoria (1.0 × 1010 m-2-3.5 × 1010
m-2). The microlite crystallinity of the Kukusan scoria (36%-48 %) is slightly higher than that of the Genteng
scoria (29%-31 %). The mean aspect ratio of the Kukusan scoria (4.3-4.8) is similar to that of the Genteng scoria(3.6-10.2).

10

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

Fig. 8 (a) Scanned thick section of the Kukusan modeD scoria. The red square indicates the approximate location of (c). (b) Scanned
thick section of the Genteng mode scoria. The red square indicates the approximate location of (d). (c) 1000× BSE image of
microlites in the Kukusan scoria. Mean microlite size, MND, and microlite crystallinity are 3.8 µm, 7.9 × 1010 m-2, and 41%,
respectively. (d) 1000× BSE image of microlites in the Genteng scoria. Mean microlite size, MND, and microlite crystallinity are
7.9 µm, 1.0 × 1010 m-2, and 29%, respectively.

Ⅴ. Discussion
5.1. Eruption types
The Kukusan and Genteng scoria deposits are loose-packing, well-sorted, and lack accretionary lapilli,suggesting that scoria is a fall deposit and deposited in a dry eruption (Walker and Croasdale 1971). Eruption
types are classified by dispersal area (D) and fragmentation degree (F) (Walker 1973). Based on the isopach map
constructed herein, the dispersal area of Kukusan is estimated to be 64 km2 by extrapolating the area enclosed by
the 1 cm thickness isopach (Fig. 9), whereas that of Genteng is estimated to be 0.30 km2 by the dispersal area ofGenteng scoria (Fig. 3b). We calculated the fragmentation degree using the weight % of clasts smaller than 1 mm.
Based on the grain size analysis results, the fragmentation degree of Kukusan equals 1-7 wt.% at location I13,
I15, I89 and I99, whereas that of Genteng equals approximately 16 wt.% at location I43. It should be noted thatthe fragmentation degree of the Genteng scoria was determined solely from the samples deposited close to thevent. Based on the dispersal area and fragmentation degree, we concluded that the Kukusan eruption type was
sub-Plinian, while that of Genteng was Strombolian (Fig. 10).

Fig. 9 Area enclosed by the Kukusan isopach against the thickness of the Kukusan scoria deposits. The maximum thickness isopach
is 100 cm. The area enclosed by 1 cm isopach is estimated to be 64 km2 by extrapolation.

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

11

Fig. 10 Diagram showing the classification of eruption types by dispersal and fragmentation (Walker 1973). Kukusan is classified as
the sub-Plinian eruption, while Genteng is classified as the Strombolian eruption.

5.2. Comparison with other mafic explosive eruptions
Mafic magma is characterized by various types of explosive eruption type: Hawaiian (Kilauea volcano, USA,1984-1986), Strombolian (Stromboli volcano, Italy, 2002; Genteng), sub-Plinian (Izu-oshima, Japan, 1986;
Kukusan), and Plinian (Etna volcano, Italy, 122 BC). The magmatic compositions of these volcanoes are similar
in terms of the SiO2 content (50-55 wt.%). However, their eruption types differ. It permits us to ask a question,“What is the main factor that influences the eruption type of these volcanoes”?
We compared the BND, MND, and microlite crystallinity of the Strombolian and (sub-)Plinian eruptions
using the available data. The BND of Kukusan is similar to that of Etna and Izu-Oshima, and one order ofmagnitude higher than that of Genteng (Fig. 11). The MND of Kukusan is lower than that of Etna, similar to thatof Izu-Oshima and higher than that of Genteng (Fig. 12). The microlite crystallinity of Kukusan is much lower
than that of Etna, but similar to that of Izu-Oshima and slightly higher than that of Genteng. Based on these facts,
we suggest that BND, MND, and microlite crystallinity correlate positively with the column height and represent
the factors that are related to the eruption types (from Strombolian to Plinian).

Fig. 11 Comparison of bubble number density for explosive eruption against the magma SiO2 content. The dashed line is the contourof decompression rate with a water content of 4 wt.%. (Toramaru 2006). References: (1) this study, (2) Mangan and Cashman(1996), (3) Lautze and Houghton (2005, 2007, 2008), (4) Toramaru (2006), (5) Sable et al. (2006), and (6) Criswell (1987), andKlug and Cashman (1994).

12

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

Fig. 12 Comparison of microlite number density for explosive eruption against microlite crystallinity (bubble-free). Symbol size
indicates the average length of microlite. References: (1) this study, (2) Sable et al. (2006), and (3) Toramaru et al. (2008).

5.3 Ascent rate
The ascent rate can be regarded as proportional to the decompression rate, assuming steady-state flow
conditions. The decompression rate is calculated as a function of SiO2 and water contents, and BND (Toramaru
2006). The SiO2 content of Kukusan is similar to that of Genteng with the difference being approximately 3 wt.%.
The initial water content of Kukusan and Genteng is unknown: however, it is unlikely to exert a significant
influence on the decompression rate (Toramaru 2006). Therefore, we deem that BND influences thedecompression rate the most. With the BND of the Kukusan scoria up to one order of magnitude higher than thatof Genteng scoria, the Kukusan decompression rate is up to one order of magnitude higher than that of Genteng(Fig. 11).
Microlite crystallization is induced by magma cooling or decompression (Toramaru 2019). Magmatemperature is considered almost constant when magma ascends rapidly in the conduit. A rate of water exsolution
from a melt, which is proportional to a decompression rate, is calculated as a function of SiO2 and water contents,and MND, assuming that the water exsolution proceeds at equilibrium (Toramaru et al. 2008). Since pyroxene
nucleates homogeneously more than plagioclase (Toramaru 2019), we used pyroxene microlite to calculate MND.
The water exsolution rate (dCw/dt) was calculated using the following formula:(5)

where NVmic is the MND of pyroxene; CSi is the SiO2 content (wt.%); Cw is the initial water content (wt.%).
NVmic was calculated using the following formula:(6)

where d is the equivalent mean microlite size. NVmic may have small error because the microlite shape is not
typically spherical but rather tabular or square. The pyroxene MND of Kukusan (1.5 × 10 13 mm -3) isapproximately one order of magnitude higher than that of Genteng (9.9 × 1011 mm-3). The initial water contents ofEtna and Stromboli are up to 4 wt.% and 3.4 wt.%, respectively (Houghton and Gonnermann 2008). Since theinitial water contents of Kukusan and Genteng were not measured, we assumed a water content of 3 wt. %. The

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

13

water exsolution rate of Kukusan (9 × 10-1 wt.% /s) is approximately one order of magnitude higher than that ofGenteng (4 × 10-2 wt.% /s). Similar results were obtained for other water contents (1-4 wt.%), suggesting that thedecompression rate of Kukusan is much higher than that of Genteng.
The number density of bubbles and microlites demonstrated that the ascent rate of the Kukusan magma ismuch higher than that of the Genteng magma. There are a few factors to affect an eruption explosivity, i.e.
column height (e.g., conduit geometry) (Cassidy et al. 2018). However, a high ascent rate is one of the mostimportant factors that increase the column height. It is suggested that the higher BND and MND values of the
Kukusan scoria are related to the cause to induce the highly explosive eruption of Kukusan volcano.

5.4. Viscosity of melt containing microlites
The rheology of magma is essential when considering conduit dynamics. The viscosity of the pure melt phase(η0) is studied well as Newtonian fluid (e.g., Hui and Zhang 2007) and it can be estimated from temperature (T),water content (Cw) and geochemical composition of the melt. When the magma contains microlites with a rangeof aspect ratios, it behaves as a non-Newtonian fluid so that the magma viscosity increases. The relative viscosity(ῆ), which is the microlite-bearing magma viscosity normalized by the viscosity of the pure melt phase, can be
estimated from crystallinity (Ømic) and crystal aspect ratio (rp) (Mader et al. 2013). The aspect ratio significantly
affects a maximum packing fraction of crystal (ØMmic).
In order to assume the temperatures in magma chambers, the chemical compositions of Kukusan and Genteng
were applied to Rhyolites-MELTS (v.1.2.0) to calculate liquidus temperature. Assuming a water content of 3
wt.% and constant pressure of 1 kbar (corresponds to a depth of approximately 3 km), the Kukusan magma
liquidus equals 1084°C. The magma temperature does not change significantly (1091°C) assuming a watercontent of 3 wt.% and constant pressure of 2 kbar (corresponds to a depth of approximately 7 km). Under thesame conditions of initial water content (3 wt.%) and pressure (1 kbar) as the Kukusan magma, the Genteng
magma liquidus equals 1077 °C. In both volcanoes, the magma temperature was assumed to be less than 1100°C.
As phenocrysts are thought to be crystallized in the magma chamber, the temperature in the conduit tends to belower than the liquidus.
We estimated the viscosity of Kukusan and Genteng melt-phase magma (η0) using the chemical compositionsof the Kukusan and Genteng melts for two cases of magma temperatures at 1100°C and 1000°C. The chemicalcomposition of each melt was calculated from the chemical composition of the bulk rock and crystallinity ofphenocrysts. For example, the SiO2 content of each melt was determined using the following formula:(7)

where i represents the phenocrysts of Pl, Opx, Cpx, Ol, and Ox. Assuming a water content of 3.0 wt.% and amagma temperature of 1100°C, the melt viscosity of Kukusan was estimated as log η0 = 1.9, while that ofGenteng was estimated as log η0 = 2.0. Assuming a water content of 3.0 wt.% and a magma temperature of 1000
°C, the melt viscosity of Kukusan was estimated as log η0 = 2.7, while that of Genteng was estimated as log η0 =2.8. With similar magma temperatures, the Kukusan and Genteng melt viscosities do not differ significantly.
Because the phenocryst crystallinities of Kukusan and Genteng scoriae are low (< 10%), phenocryst may barely
affect the melt viscosity.
The relative viscosities (ῆ) of Kukusan and Genteng microlite-bearing magmas were estimated using the Ømicand rp values (Mader et al. 2013). The Kukusan magma behaves as a Herschel-Bulkley fluid because the Ømic/ØMmic ratio exceeds 0.9 (Table 4). Assuming that the strain rate is 1 s-1 and the yield stress is 0 Pa, the relative
Kukusan viscosity is the second or third power of ten (log ῆ = 2.1-3.1). When the yield stress is greater than 0 Pa,
the relative viscosity increases. On the contrary, the Genteng magma behaves as a power-law fluid, because theØmic/ØMmic ratio is 0.6 (Table 4). The relative Genteng viscosity is less than ten (log ῆ = 0.8). The effective magma
viscosity (η), which is microlite-bearing magma viscosity, is assumed to be a multiplication of the pure melt
phase viscosity (η0) and the relative viscosity (ῆ). The effective viscosity (η) of Kukusan is 104-106 Pa s, whereasthat of Genteng is 103-104 Pa s. The microlite crystallinity of the Kukusan scoria (36%-48%) is slightly higher
than that of the Genteng scoria (29%-31%). However, the effective viscosities in these two volcanoes differ
dramatically. If Kukusan and Genteng volcanoes had the same temperature and initial water content, the effective
viscosity of Kukusan should be up to two orders of magnitude higher than that of Genteng.

14

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

Table 4 Relative viscosity of the Kukusan and Genteng microlite-bearing magmas estimated following Mader et al. (2013) andusing the microlite crystallinity and aspect ratio.

5.5 Degassing process
We analyze the relationship between degassing in the ascending magma and eruption types by comparing the
sub-Plinian eruption in Kukusan and the Strombolian eruption in Genteng. Fig. 13a illustrates a vertical crosssection of a conduit containing a bubble and some microlites during the upward transport of magma. Bubble
ascent and expansion tends to be restricted by the high viscosity of the melt. To determine the mobility of bubbles
in ascending the melt, it is crucial to examine the bubble ascent velocity relative to the magma ascent velocity.
When the bubble velocity is high, decoupled degassing may occur. Conversely, when the bubble velocity is nearlyzero, coupled degassing occurs. Assuming a spherical shape of bubbles and a steady-state flow of the melt, thebubble velocity is thought to be the Stokes velocity. The Stokes velocity is negatively proportional to magma
viscosity. Because the effective viscosity of the Kukusan magma is higher, the Stokes velocity should be smaller.
Therefore, the Stokes velocity in Kukusan is estimated to be low. In addition, since the velocity of magma ascentis quite high, the time available for bubble migration relative to melt is short. Hence, the bubbles may become
more immobile in the ascending melt. Consequently, the Kukusan magma undergoes coupled degassing,producing the sub-Plinian eruption (Fig. 13b). Instead, the effective viscosity of the Genteng magma wasevaluated to be relatively low. In addition, the mean Genteng bubble size is larger than that of Kukusan. As the
Stokes velocity is negatively proportional to magma viscosity and is positively proportional to the square of the
bubble size, it is higher in the Genteng magma than in Kukusan one. With the lower ascent rate of the Genteng
magma, the bubbles migrate upwards easier in the Genteng magma. Consequently, the Genteng magma undergoes
decoupled degassing and loses its explosiveness, producing the Strombolian eruption.

5.6. Fragmentation
The sub-Plinian and Plinian eruptions of mafic magma are thought to be associated with brittle fragmentation(Sable et al. 2006, Houghton and Gonnermann 2008). A crystallization experiment showed that the criteria offragmentation are given as an initial magma temperature of less than 1100°C, and a crystal content of over 30 %,which corresponds to a magma viscosity of approximately 105 Pa s (Arzilli et al. 2019). Applying this criteria to
our results, it is found that the Kukusan magma satisfies these conditions: initial magma temperature less than
1100°C; crystal content 36-48%; and magma viscosity 104-106 Pa s. Even though the initial magma temperatureof Genteng is similar to that of Kukusan, the crystallinity of Genteng is lower than that of Kukusan. Thus, theGenteng effective viscosity is low and is probably insufficient to cause fragmentation. Therefore, the Genteng
magma could not reach fragmentation, unlike the Kukusan magma.

Ⅵ. Conclusions
・Based on the dispersal area and fragmentation degree of scoria, we argue that the eruption types of Kukusan
and Genteng volcanoes are sub-Plinian and Strombolian, respectively.
・The microtextural analysis shows that the Kukusan scoria has a higher bubble number density, microlite
number density, and microlite crystallinity than the Genteng scoria. Assuming the same water contents, theascent rate and effective viscosity of the Kukusan magma are higher than those of the Genteng magma.
・We ascertain that the eruption type of basaltic andesite magma is related to its rapid magma ascent and magma
viscosity. The higher ascent rate and viscosity of the Kukusan magma cause coupled degassing andfragmentation, leading to a more explosive sub-Plinian eruption. Conversely, the lower ascent rate andviscosity of the Genteng magma cause decoupled degassing, leading to the Strombolian eruption.

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

15

Fig. 13 (a) Vertical cross-section of a conduit (schematic) containing a bubble, some microlites, and melt. High microlite
crystallinity increases the melt viscosity and suspends the bubble ascent. (b) Sketch of the Kukusan and Genteng conduits. Higher
ascent rate and higher effective viscosity of the Kukusan magma cause bubble immobility relative to the melt. The effective
viscosity is high so that reaches magma fragmentation. Lower ascent rate and lower effective viscosity of the Genteng magma
make bubbles ascend easier.

Acknowledgements
The authors are grateful to T. Miyamoto, K. Shimada, and Sariyanto for helping with sample analysis. We
thank T. Ikeda for his critical comments to clarify this work. We thank M. Nishiwaki and G. N. R. Bunga Naenfor their fieldwork help. We also thank M. Ohashi for his helpful review of the manuscript and relevant
comments.

References
Arzilli, F., La Spina, G., Burton, M. R., Polacci, M., Le Gall, N., Hartley, M. E., Di Genova, D., Cai, B., Vo, N. T.,
Bamber, E. C., Nonni, S., Atwood, R., Llewellin, E. W., Brooker, R. A., Mader, H. M. and Lee, P. D. (2019)
Magma fragmentation in highly explosive basaltic eruptions induced by rapid crystallization. Nature
Geoscience, 12(12), 1023–1028.
Cassidy, M., Manga, M., Cashman, K. and Bachmann, O. (2018) Controls on explosive-effusive volcanic
eruption styles. Nature Communications, 9(1), 2839.

16

Takeru Mitsuoka, Atsushi Toramaru, Agung Harijoko, and Haryo Edi Wibowo

Caudron, C., Syahbana, D. K., Lecocq, T., Van Hinsberg, V., McCausland, W., Triantafyllou, A., Camelbeeck, T.,
Bernard, A. and Surono (2015) Kawah Ijen volcanic activity: a review. Bulletin of Volcanology, 77(3), 16.
Costantini, L., Houghton, B. F. and Bonadonna, C. (2010) Constraints on eruption dynamics of basaltic explosive
activity derived from chemical and microtextural study: The example of the Fontana Lapilli Plinian eruption,
Nicaragua. Journal of Volcanology and Geothermal Research, 189(3–4), 207–224.
Criswell, C. W. (1987) Chronology and pyroclastic stratigraphy of the May 18, 1980, eruption of Mount St.
Helens, Washington. Journal of Geophysical Research, 92(10), 237–266.
Hamilton, W. B. (1979) Tectonics of the Indonesian region. Washington, DC, 345pp.
Handley, H. K., Macpherson, C. G., Davidson, J. P., Berlo, K. and Lowry, D. (2007) Constraining fluid andsediment contributions to subduction-related magmatism in Indonesia: Ijen Volcanic Complex. Journal of
Petrology, 48(6), 1155–1183.
Higgins, M. D. (2000) Measurement of crystal size distributions. American Mineralogist, 85(9), 1105–1116.
Higgins, M. D. (2002) Closure in crystal size distributions (CSD), verification of CSD calculations, and thesignificance of CSD fans. American Mineralogist, 87(1), 171–175.
Houghton, B. F., Wilson, C. J. N., Del Carlo, P., Coltelli, M., Sable, J. E. and Carey, R. (2004) The influence ofconduit processes on changes in style of basaltic Plinian eruptions: Tarawera 1886 and Etna 122 BC. Journalof Volcanology and Geothermal Research, 137(1-3 SPEC. ISS.), 1–14.
Houghton, B. F. and Gonnermann, H. M. (2008) Basaltic explosive volcanism: Constraints from deposits andmodels. Chemie der Erde, 68(2), 117–140.
Hui, H. and Zhang, Y. (2007) Toward a general viscosity equation for natural anhydrous and hydrous silicate
melts. Geochimica et Cosmochimica Acta, 71(2), 403–416.
Kereszturi, G. and Nemeth, K. (2012) Monogenetic Basaltic Volcanoes: Genetic Classification, Growth,
Geomorphology and Degradation. Updates in Volcanology - New Advances in Understanding Volcanic
Systems, 16–17.
Klug, C. and Cashman, K. V. (1994) Vesiculation of May 18, 1980, Mount St. Helens magma. Geology, 22(5),468–472.
Lautze, N. C. and Houghton, B. F. (2005) Physical mingling of magma and complex eruption dynamics in theshallow conduit at Stromboli volcano, Italy. Geology, 33(5), 425–428.
Lautze, N. C. and Houghton, B. F. (2007) Linking variable explosion style and magma textures during 2002 at
Stromboli volcano, Italy. Bulletin of Volcanology, 69(4), 445–460.
Lautze, N. C. and Houghton, B. F. (2008) Single explosions at Stromboli in 2002: Use of clast microtextures tomap physical diversity across a fragmentation zone. Journal of Volcanology and Geothermal Research,
170(3–4), 262–268.
Mader, H. M., Llewellin, E. W. and Mueller, S. P. (2013) The rheology of two-phase magmas: A review andanalysis. Journal of Volcanology and Geothermal Research, 257, 135–158.
Mangan, M. T. and Cashman, K. V (1996) Joumalof volcanology and geothermal research The structure ofbasaltic scoria and reticulite and inferences for vesiculation, foam formation, and fragmentation in lava
fountains. Journal of Volcanology and Geothermal Research, 73(1–2), 1–18.
Sable, J. E., Houghton, B. F., Del Carlo, P. and Coltelli, M. (2006) Changing conditions of magma ascent andfragmentation during the Etna 122 BC basaltic Plinian eruption: Evidence from clast microtextures. Journal ofVolcanology and Geothermal Research, 158(3–4), 333–354.
Sahagian, D. L. and Proussevitch, A. A. (1998) 3D particle size distributions from 2D observations: Stereology
for natural applications. Journal of Volcanology and Geothermal Research, 84(3–4), 173–196.
Sitorus, K. (1990) Volcanic stratigraphy and geochemistry of the Idjen Caldera Complex, East Java, Indonesia.
Unpublished MSc thesis, University of Wellington, New Zealand.
Suhendro, I., Harijoko, A. and Bunga Naen, G. N. R. (2016) Karakteristik Batuan Hasil Gunung Api Dalam
Kaldera (Intra Caldera) Ijen, Desa Kalianyar, Kecamatan Sempol, Kabupaten Bondowoso. Unpublished
Boucher thesis, Gadjah Mada Universitiy, Indonesia.
Sujanto, Syarifuddin, M. Z. and Sitorus, K. (1988) Geological Map of the Ijen Caldera Complex, East Java.
unpublished.
Sundhoro, H. (1990) A study of the stratigraphy, volcanology and geochemistry of proclastic rocks from the Ijen
caldera complex, East Java, Indonesia. Unpublished MSc thesis, University of Wellington, New Zealand.
Toramaru, A. (2006) BND (bubble number density) decompression rate meter for explosive volcanic eruptions.
Journal of Volcanology and Geothermal Research, 154(3–4), 303–316.
Toramaru, A., Noguchi, S., Oyoshihara, S. and Tsune, A. (2008) MND(microlite number density) water
exsolution rate meter. Journal of Volcanology and Geothermal Research, 175(1–2), 156–167.
Toramaru, A. (2019) Vesicularion and Crystallization of Magma: Fundamentals of the Volcanic Eruption Process.

Eruption Types and Conduit Dynamics of Kukusan and
Genteng Volcanoes of the Ijen Volcanic Complex, Indonesia

17

Tokyo University, 493pp.
Walker, G. P. L. (1971) Grain-Size Characteristics of Pyroclastic Deposits. The Journal of Geology, 79(6), 696–
714.
Walker, G. P. L. (1973) Explosive volcanic eruptions - a new classification sheme. Geologische Rundschau,
62(2), 431–446.
Walker, G. P. L. and Croasdale, R. (1971) Characteristics of some basaltic pyroclastics. Bulletin Volcanologique,35(2), 303–317.

Memoirs of the Faculty of Sciences, Kyushu University
Series D, Earth and Planetary Sciences
This publication is published irregularly. All the back issues are listed below.

Back Issues
Memoirs of the Faculty of Science, Kyushu Imperial University. Series D, Geology.
Vol. 1 No. 1 (July, 1940), No. 2 (Mar., 1941), No. 3 (July, 1942).
Vol. 2 No. 1 (Feb., 1943), No. 2 (Aug.,1944).
Vol. 3 No. 1 (Apr., 1947).
Memoirs of the Faculty of Science, Kyushu University. Series D, Geology.
Vol. 3 No. 2 (Nov., 1949), No. 3 (May.,1952), No. 4 (Dec., 1952).
Vol. 4 No. 1 (June, 1954), No. 2 (July, 1954).
Vol. 5 No. 1 (Aug., 1954), No. 2 (Oct., 1954), No. 3 (Dec., 1955), No. 4 (Jan., 1957).
Vol. 6 No. 1 (Feb., 1957), No. 2 (May.,1957), No. 3 (May.,1958).
Vol. 7 No. 1 (Mar., 1958).
Vol. 8 No. 1 (Mar., 1958), No. 2 (May.,1958), No. 3 (Mar.,1959), No. 4 (May.,1959).
Vol. 9 No. 1 (Dec., 1959), No. 2 (Nov., 1959), No. 3 (Mar.,1960).
Special Issue No. 1 (Nov., 1959), No. 2 (Jan., 1960).
Vol. 10 No. 1 (Dec., 1960), No. 2 (Mar.,1961).
Vol. 11 No. 1 (Mar., 1961), No. 2 (Nov., 1961), No. 3 (Dec., 1961).
Vol. 12 No. 1 (Mar., 1962), No. 2 (June, 1962), No. 3 (June, 1962).
Vol. 13 No. 1 (Feb., 1963).
Vol. 14 No. 1 (Jan., 1963), No. 2 (Mar.,1963), No. 3 (Jan., 1964).
Vol. 15 No. 1 (June, 1964), No. 2 (Mar.,1965).
Vol. 16 No. 1 (Mar., 1965), No. 2 (May.,1965), No. 3 (Nov., 1965).
Vol. 17 No. 1 (Nov., 1965), No. 2 (Dec., 1965), No. 3 (Sept.,1966).
Vol. 18 No. 1 (Feb., 1967), No. 2 (Dec., 1967).
Vol. 19 No. 1 (Jan., 1969), No. 2 (Jan., 1969), No. 3 (Nov., 1969).
Vol. 20 No. 1 (Jan., 1970), No. 2 (Nov., 1970).
Vol. 21 No. 1 (Oct., 1971), No. 2 (Dec., 1972).
Vol. 22 No. 1 (Nov., 1973), No. 2 (Feb., 1975).
Vol. 23 No. 1 (Mar., 1975), No. 2 (Nov., 1975), No. 3 (Feb., 1977).
Vol. 24 No. 1 (Nov., 1978), No. 2 (Nov., 1979), No. 3 (Jan., 1981), No. 4 (Dec.,1981).
Vol. 25 No. 1 (Nov., 1983), No. 2 (Nov., 1984), No. 3 (Feb., 1985).
Vol. 26 No. 1 (Jan., 1986), No. 2 (Jan., 1988), No. 3 (Dec., 1989).
Memoirs of the Faculty of Science, Kyushu University. Series D, Earth and Planetary Sciences.
Vol. 27 No. 1 (Jan., 1991), No. 2 (Feb., 1992).
Vol. 28 No. 1 (Dec., 1993), No. 2 (Dec., 1994).
Vol. 29 No. 1 (Dec., 1995), No. 2 (Dec., 1996).
Vol. 30 No. 1 (Jan., 1998), No. 2 (Dec., 1998), No. 3 (Dec., 1999).
Vol. 31 No. 1 (Dec., 2000), No. 2 (Feb., 2002), No. 3 (Feb., 2005), No. 4 (Feb., 2007).
Vol. 32 No. 1 (Feb., 2008), No. 2 (Mar.,2009), No. 3 (Mar.,2011), No. 4 (Mar.,2012).
Vol. 33 No. 1 (Dec.,2013).
Vol. 34 No. 1 (Nov.,2016).
Vol. 35 No. 1 (Mar.,2021).
Published by
FACULTY OF SCIENCES
Kyushu University
Fukuoka, Japan

令和3年3月26日

発行

編 集 兼
発 行 者

九州大学大学院理学研究院〒812-8581
福岡市東区箱崎6丁目10番1号

編

株 式 会 社 ミドリ 印 刷
福岡県福岡市博多区博多駅南6丁目17-12

集

TYPESETTING PERFORMED AT MIDORI PRINTING CO., LTD., FUKUOKA