Hydrovolcanic activity on a continental shelf inferred from the subsurface diatreme- and crater-filling deposits of Jeju Island, Korea

Jeju Island comprises numerous tuff rings and tuff cones and their reworked deposits in the subsurface, which formed on the ca. 120-m-deep Yellow Sea continental shelf under the fluctuating Quaternary sea levels. Tens of meter-thick and massive deposits were found by chance during groundwater drilling at three sites. These deposits are interpreted as either syn-eruptive diatreme-filling deposits or post-eruptive crater-filling deposits, both of hydromagmatic volcanoes. The diatremes were cut into shelf sediment, 70 to 250 m thick, and developed generally within it below the Quaternary sea levels. Abundant external water was therefore available for explosive magma-water interactions at shallow levels. The diatreme deposit in one core shows some features attributable to extreme wetness or water saturation of the diatreme fill, such as the matrix support of larger clasts, meager vertical changes in matrix content, and an absence of features related to particle adhesion. Fluidally shaped clasts with delicate reentrant margins in the core suggest minimal particle abrasion and breakage in a water-saturated and highly fluid slurry of tephra and water that was probably filling a shallow bowl-like diatreme, which is distinguished from both phreatomagmatic and kimberlite diatremes. The diatreme deposits in other cores comprise blocky and angular clasts in a sideromelane ash matrix, suggesting phreatomagmatic explosions at a deeper level. One of the cores contains collapsed deposits of thinly stratified tuff emplaced by pyroclastic surges, indicating that the diatreme is associated with an emergent tuff ring. Both Surtseyan and phreatomagmatic eruptions are therefore interpreted to have occurred on the shelf under the controls of fluctuating Quaternary sea levels. The subsurface diatremes suggest that there can be a variety of diatremes with different sizes, shapes, and material characteristics beneath the craters of hydromagmatic volcanoes, including not only maars but also tuff rings and tuff cones.


Introduction
A diatreme is a conically shaped, tephra-filled volcanic structure beneath the crater of a volcanic edifice, most of which are inferred to be produced by phreatomagmatic eruptions (Lorenz et al. 1970(Lorenz et al. , 2016Lorenz 1986Lorenz , 2003White and Ross 2011;Valentine and White 2012). A complicated set of processes operate to form the diatreme structures and infilling deposits, such as explosive interactions of magma and water or wet sediments, excavation and disturbance of the substrate, forceful upward injection of tephra and gas followed by falling back of debris, collapse of wall rocks and volcanic edifices on the surface, and convective mixing of the materials associated with these processes (Lorenz 1986(Lorenz , 2003White 1991b;Ross and White 2006;Lorenz and Kurszlaukis 2007;White and Ross 2011;Gernon et al. 2013;Ross et al. 2013). Understanding the processes of diatreme formation and deposit emplacement is important not only for unraveling magmawater interactions and related eruption processes but also for mineral exploration (Leckie et al. 1997;Stiefenhofer and Farrow 2004;Zonneveld et al. 2004;Brown et al. 2008b) and analysis of volcaniclastic-filled sedimentary basins (Kwon and Sohn 2008;Son et al. 2012;Lutz et al. 2013).
Although phreatomagmatic eruptions are one of the most common volcanic phenomena on Earth, and many of them are presumed to form diatremes (Fisher and Schmincke 1984;Cas and Wright 1987;Vespermann and Schmincke 2000;Lorenz et al. 2016), case studies of diatremes are greatly outnumbered by those of surficial tephra rings or cones (White and Ross 2011). In addition, non-kimberlite diatremes have been only rarely explored by drilling (Elliott et al. 2015;Jones et al. 2017;Bolós et al. 2021), although they remain mostly in the subsurface. For kimberlite diatremes, conflicting models remain active in spite of decades of research on open pit and underground mines and cores (e.g., Clement 1982;Stiefenhofer and Farrow 2004;Sparks et al. 2007; Wilson and Head 2007a;Wilson and Head 2007b;Cas et al. 2008;Brown et al. 2008a;Porritt et al. 2008a, b;Brown et al. 2009;White and Ross 2011;Kjarsgaard et al. 2022).
Jeju Island, Korea, is a large basaltic shield volcano situated on the ca. 100-m-deep continental shelf of the southeastern Yellow Sea (Fig. 1a). Extensive borehole drilling since the 1960s reveals that the shield-forming lavas on the surface of the island are underlain by about 100 m of hydrovolcanic deposits, which consist of primary volcaniclastic deposits formed in tuff rings and tuff cones and secondary volcaniclastic deposits that were resedimented in either subaerial or submarine environments (Sohn and Park 2004;Sohn et al. 2008;Jeon et al. 2013). These deposits are interpreted to have resulted from either Surtseyan or other phreatomagmatic activity on the Yellow Sea continental shelf under the fluctuating Quaternary sea levels over a million years before the effusion of the shield-forming lavas on the surface .
In this paper, we describe three selected cores from Jeju Island, which contain previously undescribed hydrovolcanic deposits that are tens of meters thick and massive. Such thickbedded and massive deposits are not found in any of the rim beds of the hydromagmatic volcanoes on the surface and are not known to have been generated by Surtseyan to phreatomagmatic or maar-forming eruptions (White and Ross 2011). They are therefore inferred to have accumulated in local topographic depressions, more than tens of meters deep, i.e., in the craters or diatremes of tuff rings, tuff cones, or maars. We presume that these cores fortuitously penetrated the diatremes of subsurface volcanoes, providing a rare opportunity to observe the features of diatreme fills and crater-filling sediments of hydromagmatic volcanoes. This paper describes and interprets the deposit features of the diatreme-and crater-filling sediments of the subsurface hydromagmatic volcanoes and discusses the controls of the fluctuating Quaternary sea levels on the nature of hydrovolcanic activity on a continental shelf.
important agents of sediment deposition, erosion, and resuspension. Extensive fields of tidal sand ridges are distributed in places on the shelf, including the Jeju Island area (Liu et al. 1998;Shinn et al. 2007). The island consists of an elliptical shield volcano, ca. 70 × 30 km 2 in size, 1950 m in height, and dotted with over 300 scoria cones, a few dozens of hydromagmatic volcanoes, and a few lava domes. The lavas have the compositions of mainly alkali basalts to trachytes and subordinately sub-alkali basalts to andesites (Park 1994;Park et al. 1999;Tatsumi et al. 2005). They erupted generally between ca. 1.0 and 0.01 Ma (Koh et al. 2008(Koh et al. , 2019Koh and Park 2010a, b;Yeo et al. 2019). Holocene to historical eruptions are also increasingly reported (Sohn et al. 2015;Marsden et al. 2021). Geochemical and geophysical studies suggest that Jeju Island is neither a hot spot nor a subduction zone or arc volcano but an intraplate volcano that lacks evidence for plume activity and slab metasomatism. The island is interpreted to have resulted Fig. 1 Location and geology of the study area. a Jeju Island is an intraplate volcano situated on the ca. 100-m-deep continental shelf of the Yellow Sea. b Geologic map of Jeju Island with the locations of the studied cores (simplified after Park et al. 2000). c Distribution of subsurface volcanic cones beneath the shield-forming lavas in Jeju Island based on the compilation of core data from ca. 1400 boreholes. The diameter of individual volcanic cones was determined based on the correlation of multiple cores adjacent to each other, which are several tens to a few hundreds of meters apart. The circles, therefore, reflect the actual edifice size of the subsurface volcanoes, although the accuracy of correlation of the subsurface volcanic cones, which are commonly less than a few kilometers in diameter, cannot but be low when the cores are more than a few hundred meters away from each other. The volcanic edifice type (tuff ring, tuff cone, or scoria cone) was interpreted from the lithofacies characteristics of the volcaniclastic deposits  more likely from decompression melting of the upper mantle at the edges of the subducted Pacific and Philippine Sea plates that have been stagnant in the mantle transition zone for millions of years (Brenna et al. 2015;Ward et al. 2021).

Methods and terminology
Borehole drilling has been carried out for over half a century on Jeju Island to develop and monitor groundwater under the supervision of the Water Resources Office of the Jeju Special Self-Governing Province. The boreholes are drilled about 30 cm in diameter to insert a monitoring device. Mesoscopic structures that are observable on ordinary outcrops can thus be observed on the recovered cores. Graphic logging of the cores was made outdoors at the drilling site within days after drilling was completed because it was not possible to transport the huge cores to an indoor storage facility and because the cores are destroyed, lost, or disturbed within days to months. The graphic logging recorded the characteristics of the lavas and volcaniclastic rocks, including sedimentary structure, clast size, shape, composition, color, and vesicularity, matrix grain size, sorting, and inclusion of exotic particles such as shells, woods, and intraformational clasts. Representative core samples were brought to the laboratory and prepared for slabs and thin sections for microscopic observations of the petrographic features. Point counting of clast types and abundance was made for one of the cores using a transparent film with a grid of 10 × 10 cm area drawn at 1 cm intervals on the core.
There is a spectrum of eruption styles that result from explosive interactions of magma and water. One endmember is an eruption that involves explosive interaction of magma with groundwater or wet sediment in the subsurface and produces commonly buoyancy-dominated tephra jets or plumes and pyroclastic currents. The term "phreatomagmatic" is used for such an eruption in this paper, although the term "Taalian" is also used as a synonym (Kokelaar 1986;Sohn 1996). The term "Surtseyan" is used for the other endmember of eruptions, which occur in shallow waters, involve the interaction of magma with a mobile mix or slurry of tephra and surface water filling a shallow funnel-shaped vent, and generate mostly short-lived and inertia-dominated tephra jets (Thorarinsson et al. 1964;Thorarinsson 1967;Kokelaar 1983;Moore 1985). We use the terms "hydrovolcanic" and "hydromagmatic" as broad terms that can encompass both phreatomagmatic and Surtseyan eruptions, as have already been used for a long time (Wentworth 1938;Sheridan and Wohletz 1983;Wohletz and Sheridan 1983;Wohletz and McQueen 1984;Németh and Kósik 2020). For example, we use the term "hydromagmatic volcano" when collectively referring to all different types of volcanoes produced by "hydrovolcanic activity" or "hydrovolcanism," such as tuff rings, tuff cones, and maars.
We use the terms "tuff ring," "tephra ring," and "ejecta ring" synonymously, irrespective of whether the "ring" surrounds a crater excavated below the pre-eruption surface or not and whether the ring is composed of "tuff" or coarsergrained deposits. Some researchers seem to avoid the use of the term "tuff ring" when describing the volcaniclastic deposits around a maar crater, and instead use the terms "tephra ring" or "ejecta ring" probably because of the belief that "tuff rings" and "maars" are distinctly different volcanoes and the term "tuff ring" is inadequate to describe the volcaniclastic deposits around the crater of a "maar." As has already been argued in another paper (Sohn et al. 2021), tuff rings and maars, or maar-diatreme volcanoes, belong to the same category of volcanoes, and the term "tuff ring" can be used to describe the "tephra ring" or "ejecta ring" of a maar.

Subsurface stratigraphy and lithology
The shield-forming lavas of Jeju Island are underlain by basaltic volcaniclastic deposits named the Seoguipo Formation (Koh 1997). Deposition of the formation began at ~ 1.88 Ma and continued until ~ 0.5 Ma (Koh et al. 2013(Koh et al. , 2017a(Koh et al. , b, 2022. The formation has an average thickness of 100 m and occurs throughout the island in the subsurface except in the northeastern margin. The formation is composed mainly of overlapping rim beds of tuff rings and cones with intercalations of volcaniclastic to epiclastic, subaerial to submarine, and rarely fossiliferous deposits that accumulated under the influence of fluctuating sea levels during the Early Pleistocene (Sohn and Park 2004;Sohn et al. 2008). More than a hundred volcanic cones were identified in the formation by compiling the core data from about 1400 boreholes (unpublished data of K.H. Park) (Fig. 1c). The upper surface of the formation is highly rugged because of the protrusion of these volcanic cones. A few of them are only partly buried by the later lava flows and exposed to the surface (Sohn and Park 2005).
The Seoguipo Formation is underlain by a still poorly consolidated and unlithified sedimentary deposit named the U Formation (Koh 1997). The formation has variable thicknesses between 70 and 250 m and is composed of moderately to well-sorted quartzose sand and mud that are devoid of volcanic particles but contain rare mollusk shell fragments (Jeong et al. 2016). Calcareous nannofossils and diatoms in these sediments suggest that the formation was deposited from the Pliocene to the earliest Pleistocene (Yoon et al. 2004). The formation is therefore interpreted to be deposits of the tidal sand ridges and inter-ridge fines that accumulated on the continental shelf before the onset of volcanism in the Jeju Island area (Koh 1997;Sohn and Park 2004). Although temporally unrelated to Jeju volcanism, the formation played an important role in modulating the eruption and growth of the overlying hydromagmatic volcanoes by acting as a "soft" substrate (Sohn and Park 2005) and as a source of the "dirty" coolant for the hydroexplosions (Sohn 1996;White 1996b;Schipper et al. 2011).
The basement rocks beneath the U Formation are composed of Jurassic and Cretaceous granites in the majority of the island (Kim et al. 2002), Cretaceous to Paleogene volcanic rocks (andesitic to rhyolitic lavas and tuffs) in the eastern margin of the island, and quartzose metasedimentary rocks of unknown age in the southwestern part of the island (Go et al. 2017). These rocks occur 155 to 312 m (ca. 250 m on average) below sea level.

General characteristics of cores
In this paper, we describe three selected cores, which contain tens-of-meter-thick intervals of massive deposits (Fig. 2). The Hallim core ( Fig. 3a) was obtained from a site at an altitude of 57 m above sea level in the northwestern part of Jeju Island (Fig. 1b). The 107-m-thick upper portion of the core consists of a succession of meter-thick aa lavas and decimeter-to meter-thick pahoehoe lavas with interlayers of clinker and paleosol (or sedimentary deposits that have undergone pedogenic alteration) (Fig. 2a). The lowermost lava unit has a glassy rim at the base and is underlain by black sand and granule (Fig. 3b), composed of non-vesicular and angular sideromelane or tachylite particles that are interpreted to be hyaloclastite (Fig. 3c). The Seoguipo Formation is 55 m thick (core depth: 125-180 m) and occurs 68 m below sea level. The U Formation occurs 125 m below the sea level (core depth: > 182 m) and is composed of moderately to well-sorted, quartzose fine sand (Fig. 3d). The contact between the Seoguipo and U Formations is transitional over 2 m (core depth: 180-182 m) and characterized by chaotic mixing of basalt and mudstone clasts and wellsorted quartzose sand (Fig. 3e).
The Hankyung core ( Fig. 4a) was obtained from a site at an altitude of 140 m above sea level in the western part of Jeju Island (Fig. 1b). The upper 187-m-thick portion of the core consists of meter-to decameter-thick aa and pahoehoe lavas with rare intercalations of clinker and paleosol (Fig. 4d). Core recovery was poor between the 187-and 277-m interval above a 3-m-thick trachytic lava (Fig. 2b). Black basaltic glassy sand and granule (hyaloclastite) (Fig. 4b, c), massive mudstone with chaotically deformed basaltic sandy interbeds (Fig. 4e), and massive sand deposits were recovered from this interval with uncertain stratigraphic relationships between them (Fig. 2b). A succession of 118-m-thick basaltic volcaniclastic deposits occurs to a depth of 397 m (257 m below sea level). The U Formation at the base of the core is composed of massive and poorly consolidated mud. Partly glassy basalt clasts, inferred to be hyaloclastite or pillow fragment breccia, are intercalated between the Seoguipo and U Formations (Fig. 4g).
The Jongdal core was obtained from a site at an altitude of 198 m above sea level in the eastern part of Jeju Island (Fig. 1b). The lavas here are 273 m thick and consist of meter-to decameter-thick aa and pahoehoe lavas and intervening clinker and paleosol layers (Fig. 2c). The Seoguipo Formation is 93 m thick (core depth: 275 − 368 m; 77 to 170 m below sea level) and underlain by trachytic lava and hyaloclastite (Fig. 4f). The U Formation is not found at this site probably because the core did not reach the formation.

Lithofacies characteristics of the Seoguipo Formation
Three major lithofacies constitute the Seoguipo Formation of the studied cores: (1) massive, generally matrix-supported basaltic lapilli tuff (facies LTm), (2) steeply inclined-bedded deposit of crudely to thinly stratified lapilli tuff and tuff and massive to crudely stratified gravelly sandstone (facies LTib), and (3) massive, generally matrix-supported, volcaniclastic gravelly sandstone (facies GSm). Other volumetrically minor facies includes massive volcaniclastic sandstone (facies Sm), crudely stratified volcaniclastic sandstone (facies Ss), and massive mudstone (facies Mm). Individual cores show overall similar vertical facies changes from massive lapilli tuff (facies LTm) at the base through massive gravelly sandstone (facies GSm) in the middle to finer-grained and thinner-bedded deposits (facies Ss, Sm, and Mm) at the top (Fig. 2). Steeply inclined-bedded deposit (facies LTib) occurs only in the Hankyung core between the massive lapilli tuff and the massive gravelly sandstone facies.  Fig. 2c). It consists of basalt and sedimentary rock (mostly mudstone) fragments supported in a poorly sorted and detrital quartz-rich sideromelane ash matrix (Fig. 5). Basement-derived accidental lithic clasts are not identified in all cores. The basalt clasts are mostly aphyric, fine-to-medium lapilli-sized, and rarely coarse lapilli-to block-sized (grain size terms after White and Houghton 2006). In the Hallim core, the basalt clasts are generally glassy and poorly to moderately vesicular and have irregular shapes of rag, amoeboid, and disk with uneven margins (Fig. 5a, b). Scoriaceous bomb-like clasts with finely vesicular margins and coarsely vesicular interior are also common (Fig. 5c). In the Hankyung and Jongdal cores, the basalt clasts are non-vesicular and mostly angular and blocky in shape (Fig. 5d), although some block-sized clasts in the Jongdal core have partly irregular or fluidally shaped margins (Fig. 5e).

Description
Sedimentary clasts are composed mainly of mudstone and rarely of sandstone ( Fig. 5f, g). Mudstone fragments are fine lapilli-to block-sized, angular to subrounded and irregular in shape, light gray to light brown in color, and internally massive or crudely laminated (Fig. 5g). Some of them contain abundant shell fragments. On the other hand, sandstone clasts have irregular or contorted shapes and diffuse and disseminated margins ( Fig. 5f) and are found only rarely in the lower part of the facies unit.
Point counting of clast types and their abundance in the Hallim core shows that the core does not show meaningful changes in clast type or abundance through the facies unit except for the abundance of sedimentary clasts and shell fragments near the base of the unit (Table 1). Microscopic observations show that the matrix consists predominantly of non-to poorly vesicular sideromelane shards and subordinately of opaque tachylite particles (Fig. 6). The sideromelane particles have generally round and undeformed vesicles and fluidally shaped (Fig. 6a) to curviplanar or subrounded margins (Fig. 6c, d). Detrital quartz grains are ubiquitous. Characteristically green glauconite (glaucony) particles (Odin and Matter 1981;Banerjee et al. 2016) are also common in the quartz-rich matrix (Fig. 6b). Shell fragments are locally abundant near the contact with the U Formation (Fig. 6b).

Interpretation
The angular/blocky to irregular/fluidally shaped basalt clasts in this facies are interpreted to have formed by the mingling of ascending vesiculating to non-vesiculated magma and wet and unconsolidated, and probably liquefied sediment of the U Formation (Kokelaar 1982;Hanson and Wilson 1993;Zimanowski and Buttner 2002). The coarse mixing f Dark gray glassy particles (left) and trachytic lava (right) at the base of the Jongdal core (depth: 408 m). g Massive lapilli tuff at the base of the Seoguipo Formation overlies the massive mudstone of the U Formation in the Hankyung core (depth: 397 − 398 m). Partly glassy basalt clasts between these formations can be hyaloclastite or pillow lava fragments. The pencil for scale is 12 cm long of magma and sediment was probably followed immediately by either brittle fragmentation or ductile tearing apart of magma, similar to the blocky to fluidal peperite-forming processes (Busby-Spera and White 1987;Hanson and Hargrove 1999;Lorenz et al. 2002;McClintock and White 2006;Martin and Németh 2007;McClintock et al. 2008;Befus et al. 2009), fuel-coolant interactions, and explosive expansion of steam or vapor and dispersion of tephra (Zimanowski et al. 1997(Zimanowski et al. , 2015Zimanowski 2001). The complete absence of bedding in the lapilli tuffs, except for a possible sandy interbed in the Jongdal core (depth: 346 m; Fig. 2c), which may be a sedimentary clast larger than the core diameter, suggests that the lapilli tuffs represent an unsorted mixture of the diatreme-filling slurry. Crude stratification or bedding is expected to be found if the lapilli tuffs were emplaced by multiple pyroclastic fall or density current events on the syneruptive crater floor without later disruption (Stiefenhofer and Farrow 2004;Ross and White 2005;Lorenz 2008;Moss et al. 2008;White and Ross 2011;Elliott et al. 2015;Lorenz et al. 2016).
The preservation of ragged to discoid clast shapes and delicate reentrant margins in the Hallim core  suggests that the basalt clasts experienced minimal abrasion and breakage associated with granular mixing in the diatreme after initial magma fragmentation. The preservation of contorted sedimentary clasts with diffuse and disseminated margins (Fig. 5f) and local concentration of shells and sedimentary clasts (Table 1) also suggests incomplete mixing and homogenization of the diatreme-filling debris in addition to negligible particle abrasion and breakage in the diatreme. We, therefore, infer that the Hallim core penetrated the "dead zone" in the diatreme margins, where the materials were unaffected by the violent diatreme processes, such as debris or gas jetting from below or falling back of debris that was explosively ejected by explosions (White and Ross 2011).
We infer that the diatreme-filling materials had the properties of a mobile slurry and contained abundant water-rich or water-saturated matrix. The slurry might have been cohesive if a sufficient amount of mud, which consists of illite, kaolinite, smectite, and chlorite (Jeong et al. 2016), was introduced into the slurry from the U Formation. The lack of clast-supported to openwork texture in the lapilli tuff suggests that the fine-grained materials in the diatreme could not be efficiently winnowed by, e.g., upward-directed debris jets or gas jets. Meager vertical changes in the proportion of the matrix (Table 1) are also interpreted as evidence against the role of gas or steam fluidization in the diatreme. The absence of ash aggregates or ash-coated particles (Fig. 6a) suggests that the diatreme fill was water-saturated, and particle adhesion was inhibited because of the lack of surface tension effects in watersaturated slurry (Go et al. 2017). We, therefore, infer that the diatreme-filling materials could be locally injected upward by gas-and steam-propelled jets of debris near the loci of explosions, but these jets could not elutriate the fine-grained materials in the other part of the diatreme because these materials were wet or even water-saturated.
Vigorous mixing and particle abrasion of the diatremefilling materials are not inferred to have occurred in the diatremes of the Hankyung and Jongdal cores because the clasts in these cores also preserve sharp corners and irregular protrusions along their margins (Fig. 5d, e). Even fragile mudstone clasts appear to have been subject to negligible abrasion and breakage in the Hankyung core (Fig. 5g). The absence of features related to ash aggregation ( Fig. 6b − d) can also be the result of water saturation of the diatreme fills in these cores.

Interpretation
The meter-thick interval of wood-bearing mudstone and shell-bearing sandstone between this facies and the underlying massive lapilli tuff (Fig. 7c, e, f) indicates that this facies formed after the volcanic eruption ended. The unusually steep bedding in this facies, which far exceeds the angle of repose, is interpreted to be due to neither tectonic tilting nor the artificial effect of borehole drilling but to be the result of penetration of the borehole into steeply dipping deposits of primary tuffs and reworked volcaniclastic deposits. Steep bedding here is also not attributed to subsidence of syn-eruptive bedded upper diatreme deposits (Suhr et al. 2006;Lorenz 2008;White and Ross 2011;Delpit et al. 2014;Lorenz et al. 2016) because this facies is post-eruptive. The steep bedding in this facies is therefore attributed to the failure of the tephra beds and overlying sediments that accumulated on the surface around the volcano. A few alternations of steeply inclined-bedded and massive deposits (Fig. 2b) suggest that the failure occurred in multiple mass movement events involving rockfall or rock slide of semi-consolidated materials with minimal internal deformation and debris flow of remolded materials. The wood fragments (Fig. 7b, f) suggest that the failure occurred after vegetation development on the surface years after the cessation of the eruption. Vegetation development on the surface of the volcano at the time of eruption, which commonly lasts only days to months Siebert 1984, 2000), seems unlikely.

Description
The massive gravelly sandstone facies is 22 m thick in the Hallim core (depth: 134 − 155 m; Fig. 2a), 14 m thick in the Hankyung core (depth: 320 − 334 m; Fig. 2b), and 31 m thick in the Jongdal core (depth: 302 − 333 m; Fig. 2c). It is not clear how many depositional units constitute the massive deposits, but the presence of a mudstone interbed (or a mudstone clast larger than the core diameter) in the Hallim core (Fig. 8a), non-uniform distribution of clast types and shell fragments, and vertical alternations of clast-supported to matrix-supported textures suggest that the massive deposits consist of multiple, meter-thick, amalgamated units that have subtly different clast types and contents. The basalt clasts in this facies are either matrix-or clast-supported and show diverse colors, vesicularity, crystallinity, and angularity in the Hallim and Jongdal cores (Fig. 8b, c). However, the massive deposits of the Hankyung core comprise matrix-supported sideromelane lapilli and rare mudstone (or extremely fine-grained tuff) fragments in a sideromelane ash matrix (Fig. 8d). A lapilli tuff clast (Fig. 8e), deformed "ghost" bedding ( Fig. 8f), and wood fragments (Fig. 2b) are also found in this core. Shell fragments are found in all cores, either sparsely scattered in the massive deposits or forming indistinct layers or lenses. Under the microscope, the matrix consists predominantly of yellowish to brownish sideromelane and relatively rare opaque tachylite particles (Fig. 9). Rounded or abraded particles are more commonly found here than in the matrix of the massive lapilli tuff.
A meter-thick massive sandstone (facies Sm) is intercalated between this facies and the underlying massive lapilli tuff (facies LTm) in the Hallim core (depth 155 m; Fig. 2a). A hyaloclastite layer (poorly retrieved as core) lies between this facies and the underlying steeply inclinedbedded deposit (facies LTib) in the Hankyung core (depth 334 − 345 m; Fig. 2b). The contact with the underlying massive lapilli tuff is indistinct and gradational in the Jongdal core (depth: 333 m; Fig. 2c).

Interpretation
The massive and matrix-or clast-supported gravelly deposits are interpreted to have been deposited by either cohesive or friction-dominated debris flows (Nemec and Steel 1984;Shultz 1984;Blair and McPherson 1998;Sohn et al. 1999;Sohn 2000). The diversity of clast types (Fig. 8b), rounding of some clasts (Fig. 8c), and the inclusion of friable mud (or extremely fine-grained tuff) or lapilli tuff clasts (Fig. 8d, e) suggest that the materials for these debris flows originated from diverse sources, including primary volcaniclastic deposits surrounding the crater of a tuff ring/cone and fluvially or wave-reworked volcaniclastic deposits in a volcaniclastic apron around a volcano (Busby-Spera 1988; Karatson and Németh 2001;Allen et al. 2007;Casalbore et al. 2010). Particle rounding in the matrix is not always evident (Fig. 9), probably because of the short transport distance and insignificant reworking of tephra at the depositional site.
The emplacement setting of the debris flows is interpreted to have been subaqueous based on the intercalation of a mud layer (or a mudstone clast wider than the core diameter) in the Hallim core (Fig. 8a), the occurrence of hyaloclastite breccias beneath and above the facies units (Fig. 2), and the upward facies transition into submarine facies (Sm, Ss, and Mm; see below) (Fig. 2). The absence of steep-sided erosion surfaces incised by rill flows or overland flows is also indirect evidence for subaqueous deposition of the deposits . The sparse distribution of shell fragments cannot be the evidence for subaqueous deposition because the shell fragments could have been derived from the U Formation and later incorporated into this facies during the remobilization of volcaniclasts. However, shell-bearing Fig. 5 Clasts of massive lapilli tuff (facies LTm). a Irregular to ragged and vesicular basalt clasts and mudstone clasts (outlined by yellow line) set in a detrital quartz sand-rich ash matrix from the Hallim core (depth: 174 m). The pencil for scale is 15 cm long. b Irregular and vesicular basalt clasts in a detrital quartz sand-rich matrix from the Hallim core (depth: 168 m). c Scoriaceous bomb-like clasts from the Hallim core (depth: 160 m) having irregular and poorer vesicular margins and a highly vesicular interior. d Generally blocky and non-vesicular basalt clasts from the Hankyung core (depth: 370 m). e Non-vesicular and irregularly shaped basalt clasts from the Jongdal core (depth: 364 m). The matrix is composed of massive lapilli tuff, in which basaltic lapilli are non-vesicular and have blocky and equant shapes. f Massive lapilli tuff of the Hallim core, containing an elongated and contorted sand fragment with diffuse margins (depth: 180 m). The pencil for scale is 12 cm long. g Massive lapilli tuff containing large, massive, or laminated mudstone fragments near the base of the Hankyung core (depth: 390 m). The coin in c and e is 2.4 cm in diameter ◂ lenses or layers are interpreted to indicate a submarine environment.

Description
Core recovery was generally very poor above the massive gravelly sandstone facies in all cores (Fig. 2). Our experience suggests that this interval comprised mostly loose and friable deposits of mud, well-sorted sand, and matrix-free gravel or hyaloclastite. The fragmentary cores retrieved from this interval comprise (1) massive or crudely stratified sandstone (facies Sm and Ss) (Fig. 10a − c), (2) massive mudstone (facies Mm) (Fig. 10d), and (3) thinly to crudely stratified lapilli tuff (facies LTs) (Fig. 10a, e). The sandstones are composed of poorly to moderately sorted, fine-to coarsegrained sideromelane sand. Some sandstone layers are rich in mollusk shell fragments (Fig. 10c). The mudstones locally have burrows and commonly deformed sandstone interbeds (Fig. 4e). A stratified basaltic sandstone fragment standing upright is found in the Hallim core (Fig. 10d). The stratified lapilli tuffs show either thin stratification (Fig. 10a) or very crude stratification (Fig. 10e) with either fine ash-rich or fine ash-poor matrix.

Interpretation
The massive or crudely stratified sandstones are interpreted to have been deposited by subaqueous sandy debris flows or concentrated density flows (Mulder and Alexander 2001). Rare intercalations of shell-rich beds (Fig. 10c), common association with massive mudstone facies (Fig. 10d) or hyaloclastite deposits (Fig. 2), and the occurrence of burrows (depth: 282 − 284 m; Hankyung core; Fig. 2b) suggest collectively that these facies were deposited in a subaqueous setting. The massive mudstones also suggest suspension settling of fines in a subaqueous setting. An outsize sandstone clast (Fig. 10d) is interpreted to be rockfall debris from a nearby steep crater wall. The lack of well-developed stratification or cross-stratification is interpreted to be due to the lack of wave and current activity in a wave-protected volcanic embayment or submerged crater. The unusual intercalations of stratified lapilli tuffs are interpreted to have been derived from nearby volcanic centers by pyroclastic currents (White 1991a;Moss et al. 2008;Gernon et al. 2009). It is worth noting that these intercalations are found in the Hallim and Hankyung cores, where multiple volcanic cones are interpreted to be superposed or laterally juxtaposed in the subsurface, whereas such intercalation is absent in the Jongdal core where volcanic cones are sparsely distributed (Fig. 1c).

Eruptive and depositional setting in relation to Quaternary sea levels
Jeju Island is regarded as having been subject to negligible crustal motion, either vertical or horizontal, because the Korean Peninsula and the Yellow Sea areas surrounding Jeju Island have remained largely aseismic and tectonically stable compared with other regions in East Asia, such as the Japanese Arc and continental China (Hamdy et al. 2005;Chen 2007;Jin et al. 2007). A rough estimation of the eruptive and depositional setting of the subsurface volcanoes could therefore be made based on the depths of the Seoguipo Formation in the cores in relation to the Quaternary sea levels (Shackleton 1987;Waelbroeck et al. 2002). The contact between the Seoguipo Formation and the overlying succession of lavas occurs at a depth of 68 m below sea level in the Hallim core, 47 m below sea level in the Hankyung core, and 76 m below sea level in the Jongdal core. The majority of the Seoguipo Formation could therefore have accumulated below sea level except during the glacial maxima (Fig. 11). The facies interpretations given Fig. 6 Photomicrographs of the matrix of massive lapilli tuff (facies LTm). a A large non-vesicular sideromelane clast has a fluidally shaped reentrant margin. The matrix in the lower left comprises sideromelane and detrital quartz grains. The specimen is from the Hallim core (depth: 165 m). b A specimen from the base of the Seoguipo Formation of the Hankyung core (depth: 397 m) contains abundant shell fragments (Sh) in addition to sideromelane (S) and tachylite (T) particles. Characteristically green particles of glauconite (G), which is a common component of shallow marine sandstone, are abundant in the quartzose matrix. c Sideromelane-rich matrix from the Jongdal core (depth: 347 m). The smallest identifiable particles in the interstices are composed of non-to poorly vesicular sideromelane. d A tachylite clast in contact with the sideromelane-rich matrix from the Jongdal core (depth: 340 m) above and the common intercalations of hyaloclastite near the contact between this formation and the overlying lavas (Fig. 2) support this interpretation.
The level of the contact between the Seoguipo and the U Formations in the cores does not indicate the level of the pre-eruption surface because the contact is more likely the lateral contact between the U Formation and the diatremes that were excavated at least tens of meters to hundreds of meters below the pre-eruption surface. The top surface of the massive lapilli tuff (facies LTm) can, however, provide rough estimates of the upper level of the primary diatreme fill or the level of the syn-eruptive crater floor, which are Crudely stratified gravelly sandstone with very steep bedding (depth: 360 − 362 m), containing angular to rounded, crystalline, and heteromictic basalt gravel and a wood fragment. c Massive mudstone and interbedded shell-bearing sandstone above the massive lapilli tuff (facies LTm) and beneath the steeply inclined-bedded deposits (depth: 373 − 374 m). d Massive gravelly sandstone (facies GSm) containing angular to rounded crystalline basalt gravel (depth: 369 m). e Close-up of the shell-bearing sandstone beneath the steeply inclined-bedded deposits (depth: 373 m). f A wood fragment in massive mudstone (depth: 370 m). The pen for scale is 12 cm long. The coin in f is 2.3 cm in diameter 88 m below sea level in the Hallim core, 229 m below sea level in the Hankyung core, and 135 m below sea level in the Jongdal core. The primary diatreme fills are therefore interpreted to have formed mainly below the fluctuating Quaternary sea levels, even during the glacial maxima in the case of the Hankyung and Jongdal cores (Fig. 11). The syn-eruptive crater floors would have thus been flooded with seawater if conditions were met.
No solid evidence is available regarding the diameter, wall angle, and depth of the diatremes. The diatremes are, however, interpreted to belong to volcanic edifices with basal diameters over a few kilometers (Fig. 1c). If these volcanoes are assumed to have had a crater diameter of 1 km and a relatively gentle wall angle near the angle of repose (30°) because of the soft nature of the substrate, the syneruptive crater would have had a vertical dimension of ca. 288 m according to simple trigonometric calculations. The crater is therefore inferred to have been rooted about 200 m below the pre-eruption surface or near the contact between the U Formation and the basement if the height of the crater Fig. 8 Features of massive gravelly sandstone (facies GSm). a A massive mudstone (facies Mm) is intercalated between massive gravelly sandstone beds (Hallim core; depth: 139 − 141 m). Only part of the mudstone bed is inferred to have been retrieved as a core. b A massive gravelly sandstone from the Hallim core (depth: 155 m) consists of a variety of basalt clasts, which vary in color and texture, together with mudstone clasts (red arrows). It overlies a meter-thick massive sandstone (facies Sm). c Massive gravelly sandstone from the Jongdal core (depth: 323 m), comprising basalt clasts of different colors and textures. d Massive and matrix-supported gravelly sandstone from the Hankyung core (depth: 323 m), containing sideromelane lapilli and irregularly shaped mudstone (or extremely fine tuff) clasts in a sideromelane ash matrix. e A lapilli tuff clast is floating in the middle of a massive gravelly sandstone (Hankyung core; depth: 323 m). f A massive gravelly sandstone from the Hankyung core, showing irregularly deformed "ghost" bedding (depth: 330 m) rim, between ca. 50 and 100 m, is subtracted (Fig. 11). The absence of basement-derived accidental clasts in all cores suggests that the diatremes formed mostly within the U Formation and might have had wall angles smaller than the angle of repose, as is the case of some diatremes that formed in soft substrate environments (Lorenz 2003;Auer et al. 2007). Explosive magma-water interaction is therefore interpreted to have occurred mostly above the interface between the shelf sediment (U Formation) and the basement.
Given the known depths of some diatremes, which can be over 2.5 km deep (Lorenz 1973(Lorenz , 1975Lorenz and Haneke 2004), the diatremes in this study are very shallow, probably having a shallow bowl-like geometry. Abundant shallowlevel external water from multiple sources, including the interstitial water in the U Formation and the seawater that overflowed into the crater or permeated into the diatremes underground, could therefore be supplied to the diatremes to drive explosive magma-water interactions at shallow levels. The diatremes are therefore interpreted to have been kept shallow with limited downward penetration of the root zone. If the root zone of the diatremes penetrated deep into the basement rocks and sufficient access of groundwater or seawater to the explosion loci in the diatremes was prevented, the eruption style of these volcanoes might have changed from hydromagmatic to magmatic activity. The absence of deposits related to magmatic activity (lava and scoria) in the studied diatreme fills is consistent with the interpretation that the diatremes remained shallow throughout the eruption with the unexhausted supply of external water into the diatremes.
The absence of bedded upper diatreme deposits in all studied cores is also interpreted to be related to the shallowness of the diatreme structure. As the syn-eruptive crater floor deepens, the crater floor becomes the site of Fig. 9 Photomicrographs of the matrix of massive gravelly sandstone (facies GSm). a A specimen from the Hankyung core (depth: 323 m), composed of rounded tachylite particles and finer-grained sideromelane particles. b A specimen from the Hallim core (depth: 145 m), composed of angular to rounded sideromelane and tachylite particles. c, d A specimen from the Jongdal core (depth: 309 m), composed generally of rounded sideromelane particles pyroclastic deposition from base surges and tephra fall. These primary tephra beds, alternating with reworked tephra beds, result in the bedded upper diatreme deposits in phreatomagmatic diatremes (White and Ross 2011) and can subside to great depths of the diatremes during and after the diatreme formation (Suhr et al. 2006;Delpit et al. 2014). The absence of such bedded deposits in the studied cores suggests that the majority of the tephra ejected from the crater could be expelled out of the crater because of the low height of the crater rim relative to the syn-eruptive crater floor and the wide-angle spreading of tephra jets from the shallow diatremes. The primary diatreme fills could be therefore overlain directly by post-eruptive resedimented volcaniclastic deposits, which accumulated mostly by masswasting processes in a subaqueous setting in post-eruptive craters.

Hydrovolcanic activity on a continental shelf
The upper surface of the U Formation occurs about 120 m below sea level on average (Koh et al. 2017a(Koh et al. , b, 2022, which is almost identical to the sea levels during the glacial maxima in the Quaternary (Shackleton 1987;Waelbroeck et al. 2002). The pre-eruptive surface of the subsurface volcanoes is therefore interpreted to have been submerged in sea water at a depth of about 120 m during the interglacials, almost identical to the pre-eruptive setting of Surtsey in 1963 (Thorarinsson 1967;Kokelaar 1983;Moore 1985;Jackson et al. 2019;McPhie et al. 2020;Moore and Jackson 2020), whereas it was either barely exposed above sea level or submerged in very shallow waters during the glacials. The Yellow Sea continental shelf is therefore considered to have been an optimal place for either Surtseyan or other phreatomagmatic activity because of the low hydrostatic pressure of the water column, even during the highstands of sea level, and the availability of external water for Surtseyan to phreatomagmatic explosions, either on the surface or in the subsurface. The still unconsolidated sediments of the U Formation could also facilitate efficient mixing between ascending magma and sediment. The early-stage volcanism of Jeju Island could thus be dominated by hydrovolcanic activity, which lasted more than a million years between ~ 1.88 and ~ 0.5 Ma (Koh et al. 2013(Koh et al. , 2017a(Koh et al. , b, 2022 until proto-Jeju Island emerged above the fluctuating Quaternary sea levels in the Middle Pleistocene (Sohn and Park 2004;Sohn et al. 2008).
Either Surtseyan or other phreatomagmatic activity is inferred to have occurred on the continental shelf, resulting in either Surtseyan tuff cones or Taalian tuff rings . Judging on the basis of the general nature of hydrovolcanic activity, Surtseyan activity is supposed to have been more prevalent during the interglacials, although Surtseyan activity can occur in the glacials if abundant shallow-level groundwater is available for magma-water interactions near the surface (e.g., Hamilton and Myers 1963;Leat and Thompson 1988;Sohn and Chough 1992). Phreatomagmatic activity appears to have been possible either during the glacials or during the interglacials because explosive magma-water interaction in the subsurface and diatreme formation can occur irrespective of the presence of a water column of ca. 100 m or more above the eruption site (Leys 1983;Elliott et al. 2015). The phreatomagmatic activity associated with the diatreme of the Hankyung core is, however, interpreted to have occurred in the glacials (see below). High-resolution age data within ca. 10,000 years of error are necessary to confirm the relationship between the Quaternary sea levels and the volcano types in the subsurface.

Volcano types on the surface
The eruption styles and the volcano types on the surface associated with the subsurface diatremes are difficult to infer. However, the thinly stratified tuffs constituting the steeply inclined bedded deposits of the Hankyung core ( Fig. 7a) originated most probably from the rim beds of a tuff ring surrounding the diatreme, which were emplaced by pyroclastic surges . The diatreme of the Hankyung core is therefore interpreted to have been generated by phreatomagmatic explosions during a lowstand of sea level in the glacials. If the eruption occurred in the interglacials, the tuff ring could have been submergent and comprised subaqueous wave-or current-influenced facies (e.g., White 1996a). Post-eruptive crater floor appears to have been far below the pre-eruption surface, probably because of the almost complete emptying of the diatreme fill by the final eruptions. The crater could thus accommodate the fairly thick post-eruptive sediments (Fig. 11).
Similar information on the volcanic edifices on the surface is not available from the Hallim and Jongdal cores. However, the contrastingly different textures of juvenile clasts in these cores (Fig. 5) suggest that their eruptive conditions and settings were quite different. The moderately vesicular and commonly fluidally shaped clasts of the Hallim core are interpreted to have resulted from bulk interaction or coarse mixing (Kokelaar 1986) of vesiculating magma with a mobile slurry of tephra and water that were filling a shallow bowl-like or funnel-shaped vent, not far below the pre-eruption surface, whereas the sideromelane ash matrix resulted from molten fuel-coolant interaction (MFCI), which occurred in concert with the bulk interaction.
The non-vesicular and blocky/angular clasts of the Jongdal core, as well as those of the Hankyung core, resulted most likely from phreatomagmatic interaction of non-vesiculated magma and water in the subsurface under higher confining pressure, involving both contact-surface steam explosivity or MFCI and cooling-contraction granulation possibly after bulk interaction of magma and wet sediment (Kokelaar 1986;Houghton and Wilson 1989;Mattsson 2010). We, therefore, interpret that the volcanic edifice associated with the diatreme in the Hallim core was a Surtseyan tuff cone. If the tuff cone was never emergent above sea level because of its eruption in the interglacials, the tuff cone might more likely be described as a Surtla-type tuff cone (Kokelaar and Durant 1983;Cas et al. 1989;Németh and Kósik 2020). On the other hand, the volcanic edifice associated with the diatreme in the Jongdal core is presumed to have been a tuff ring associated with phreatomagmatic explosions in the subsurface.
The ragged or ribbon-like clasts of the Hallim core (Fig. 5a), as well as the scoriaceous bomb-like clasts with well-preserved reentrant margins (Fig. 5c), suggest that these clasts were subject to insignificant breakage and abrasion after initial fragmentation or ductile tearing apart of magma in the diatreme, although there are cases of preservation of such clasts in pyroclastic flows and debris jets, in which particle abrasion and breakage are inferred to be active (Ross and White 2005;McClintock and White 2006;McClintock et al. 2008). The abundance of matrix and the overall matrix-supported texture throughout the core also suggests that elutriation processes by debris jets or gas jets were not active. We infer that these characteristics can be the features of relatively mild clast interactions in water-saturated and highly fluid and mobile diatreme fills. We suggest that the shallow bowl-like crater and the underlying conduit filled with such water-rich slurries can be regarded as another type of diatreme that can be distinguished from the typical "phreatomagmatic diatremes" of White and Ross (2011), let alone the controversial "kimberlite diatremes" (Stiefenhofer and Farrow 2004;Sparks et al. 2007;Wilson and Head 2007a;Wilson and Head 2007b;Cas et al. 2008;Brown et al. 2008a;Porritt et al. 2008aPorritt et al. , b, 2008bBrown et al. 2009;White and Ross 2011). This study of the Jeju-style diatremes and the recent studies of the Surtsey (Jackson et al. 2019;McPhie et al. 2020;Moore and Jackson 2020) show that all types of hydromagmatic volcanoes, including tuff rings, tuff cones, and maars, can have diatremes in the subsurface with different sizes, shapes, and material characteristics. The term "diatreme" thus appears to lose its meaning in classifying or distinguishing between these hydromagmatic volcano types (also see the discussion in Lorenz et al. 2016). Fig. 10 Features of post-eruptive deposits. a Thinly stratified lapilli tuff (facies LTs) overlying massive sandstone (facies Sm) in the Hallim core (depth: 128 − 129 m). b Massive to very crudely stratified sandstone in the Hallim core (depth: 129 − 130 m). c A mollusk shellrich sand bed intercalated in a crudely stratified and poorly sorted sandstone of the Hankyung core (depth: 330 m). d Massive mudstone (facies Mm) containing an intraformational clast of stratified sandstone (Hallim core; depth: 125 − 126 m). e Crudely stratified and poorly sorted lapilli tuff from the Hankyung core (depth: 299 m) ◂ Fig. 11 Reconstruction of the diatreme position relative to the fluctuating Quaternary sea levels. No information is available regarding the diameter, wall angle, and depth of the diatremes. However, the abundance of detrital quartz sand and mudstone clasts in the diatreme deposits and the absence of basement-derived accidental lithic clasts suggest that the diatremes formed mostly within the U Formation. The wall angle of the diatremes, which is depicted as having an angle of repose (30°), can be smaller. Tuff rings/cones outside the craters are omitted for the clarity and simplicity of the diagram

Conclusions
Tens of meter-thick and massive, basaltic volcaniclastic deposits were retrieved by chance from the subsurface of Jeju Island, Korea, during groundwater drilling. They are interpreted to be either syn-eruptive diatreme-filling deposits or post-eruptive crater-filling deposits of hydromagmatic volcanoes, which formed extensively underground almost all of Jeju Island before the effusion of the lavas that built the island's present shield volcano. The pre-eruption surface of these volcanoes was a continental shelf floored by still unconsolidated sandy to muddy sediments with thicknesses between 70 and 250 m. The surface was at a depth of about 120 m below the present sea level, indicating that the surface was submerged in sea water at a depth of about 120 m during the interglacials and was either barely exposed above sea level or submerged in very shallow waters during the glacials. The setting was, therefore, optimal for either Surtseyan or other phreatomagmatic activity and the construction of either emergent or submerged tuff rings and cones.
The primary diatreme fills in three cores occur 88 m below sea level in the Hallim core, 229 m below sea level in the Hankyung core, and 135 m below sea level in the Jongdal core, suggesting that the diatremes formed mainly below the fluctuating Quaternary sea levels and that the syn-eruptive craters could be easily flooded with sea water in addition to the interstitial water seeping from the shelf sediment. The diatreme fills are devoid of basement-derived accidental lithic clasts, suggesting that the diatremes are only a few hundred meters deep and formed mostly within the shelf sediment with wall angles smaller than the angle of repose. Abundant shallow-level external water was therefore available for explosive magma-water interactions at shallow levels. The possibility of deeper explosions and disturbance of the basement rocks without the emission of the rocks to the surface cannot be ruled out. However, the basement lies fairly shallow beneath the pre-eruption surface, covered only by the unconsolidated sediment of the U Formation, which has an average thickness of 150 m (Koh et al. 2017b(Koh et al. , 2022. The accidental lithic clasts from the basement rocks above the optimal depth of phreatomagmatic explosions (Valentine et al. 2014) are therefore inferred to be present in the diatreme fill if there were active explosions and diatreme formation in the basement.
The diatremes of these volcanoes have some features that can be attributed to extreme wetness or water saturation of the diatreme fill, such as matrix support of clasts, meager vertical changes in the matrix content, and the absence of features related to particle adhesion. Preservation of ragged to discoid shapes and delicate reentrant margins of the basalt clasts in the Hallim core is especially noteworthy because these clast textures can be the result of mild clast interactions in a water-saturated and highly fluid slurry of tephra and water that was filling a shallow funnel-shaped vent and conduit, which needs to be distinguished from both phreatomagmatic and kimberlite diatremes. In these Jeju-style diatremes, abrasion and breakage of coarse particles and elutriation of fine particles appear to occur less vigorously or less efficiently than in typical phreatomagmatic diatremes.
The volcanic edifice type on the surface is interpreted to be an emergent Surtseyan-type tuff cone or a submergent Surtla-type tuff cone in the case of the Hallim core. The eruption could have occurred either in the glacials or in the interglacials. The volcanic edifice type associated with the Hankyung core is interpreted to be a tuff ring based on the post-eruptive collapsed masses of thin-bedded tuff emplaced by pyroclastic surges. The tuff ring formed most likely in the glacials upon an emergent surface. No convincing evidence is available for the Jongdal core, but a tuff ring is likely to have formed by a phreatomagmatic eruption. High-resolution age data from these volcanoes within ca. 10,000 years of error are necessary to confirm the relationship between the Quaternary sea levels and the eruption and edifice styles.
This study suggests that a variety of diatremes with different sizes, shapes, and material characteristics can be produced beneath the craters of hydromagmatic volcanoes, including not only maars but also tuff rings and tuff cones. The term "diatreme" thus appears to lose its meaning in classifying and naming hydromagmatic volcanoes because most, if not all, hydromagmatic volcanoes can have diatremes.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.