A combined geophysical and lithological study on eruptive history and Quaternary lacustrine stratigraphy of a maar in Leizhou Peninsula,China

Cong Chen,Zhuo Zheng,2* ,Li-Feng Zeng,Fan Xiao,Li-Ping Tian and Kang-You Huang,2

Abstract As the second common type of volcanic vent on Earth, maar-diatreme volcanoes and their post-eruptive lacustrine sediments are a main focus of volcanology, palaeolimnology, palaeoclimatology and palaeontology. A number of maar-type volcanoes have been found in Leizhou Peninsula, South China, but little is known about their eruption processes and detailed stratigraphy of the post-eruptive sediments. We present a combined geophysical and geological analysis to study the eruptive history and post-eruptive sediment stratigraphy of a large maar, the elliptical (1.8×3.0 km2) Jiudouyang (JDY) maar. The lacustrine stratigraphy revealed by drilling cores shows that the JDY maar lake has three major stages of evolution: (i) deep-lake sedimentary environment characterized by high autochthonous diatom productivity; (ii) shallow lake to swamp with very low water levels, characterized by a high total organic carbon (TOC) and abundant wood fragments; and, (iii) intermittent shallow lake and alluvial deposits composed of clay minerals and sand. The electrical resistivity tomography (ERT) values and lithological features are highly consistent, which clearly reveal the presence of ca. 50 m thick lacustrine sediments, directly underlain by a ca. 70 m thick basaltic lava rather than diatreme breccia in the crater. This infill sequence implies an alternation of eruption style from phreatomagmatic to Strombolian and/or lava flow, due to high magma flux and ascent rate of the Hainan Plume during the middle Pleistocene. The ERT data also reveal the initial phreatomagmatic crater floor at ca. 120 m depth. The initial crater had a large diameter/depth ratio (ca. 17), with an elongated shape (major axis to minor axis=0.6), implying possible lateral vent migration during the eruption. A significant erosion under tropical weathering condition during the last few hundred thousand years, accounted for the large size of the maar crater. The study provides insights into the eruptive history and post-eruptive evolution of a large maar, as well as the spatial distribution of the lacustrine sediments.

Keywords: Maar crater,Eruptive history, Infill sequence, Lake evolution, Lacustrine stratigraphy, Leizhou Peninsula

1 Introduction

Maar-diatreme volcanoes are the second common type of volcanic vent on Earth (Lorenz 2007). These volcanoes and the post-eruptive lacustrine sediments are a main focus of the volcanology, palaeolimnology, palaeoclimatology and palaeontology. The morphology of maar craters is governed by the sub-surface phreatomagmatic explosion (Oms et al. 2015). Lorenz (1986) suggested that their diameter is determined by syn-eruptive subsidence and a large maar crater is accompanied by a large sub-surface diatreme, with a crater diameter to diatreme depth ratio of 1:1. However, as noted by Valentine and White (2012), crater diameter is influenced by numerous factors, including magma flux, vent migration, ejection of material and subsidence of failed crater walls. In addition, craters’ diameter to depth ratios reportedly increase with age due to widening and reduction in heights by erosion (e.g., Gen?alio?lu-Ku?cu et al. 2007)and/or vent migration during eruption (e.g., Jordan et al.2013). These processes are complex, but the morphological characteristics of maar craters and the eruption products are keys to indicating the volcanic process and associated eruption environments (Graettinger 2018;Xian et al. 2018).

Most studies on maar-diatreme systems have focused on cases where the craters have an ejecta ring with limited erosion, or sub-surface structures are exposed,which can therefore be relatively easily examined (White and Ross 2011; and references therein). However, subsurface structures of maar-diatreme volcanoes with extremely weathered and eroded ejecta rings can only be modeled or investigated using geophysical approaches and boreholes (e.g., Schulz et al. 2005; Mrlina et al.2009; Gebhardt et al. 2011; Blaikie et al. 2014). An efficient and convenient geophysical method that has been widely used to image internal structures, substrates and post-volcanic stratigraphy of maar-diatreme volcanoes, is the electrical resistivity tomography (ERT) (e.g., Brunner et al. 1999; Oms et al. 2015). ERT surveying, in combination with other geophysical methods and geological techniques, can accurately characterize uppermost parts of maars to help efforts to: elucidate effects of the preeruptive substrate (Martín-Serrano et al. 2009); determine the distribution of post-eruptive sediments and the depocentre for palaeoclimatic drilling (Flechsig et al.2015); and locate the sub-surface structures to reveal the eruptive history(Bolós et al. 2012).

Maar-diatreme structures can host long-lived and climate-sensitive lakes, which are significant in palaeoclimatic and palaeoenvironmental researches (White and Ross 2011; Chu and Liu 2018; Wu et al. 2019). Maar lake basins are endorheic basins with relatively small drainage areas, where the major input is from meteoric precipitation and/or groundwater. This makes them sensitive to local and regional climate change (e.g., Mingram et al. 2004). Moreover, maar lakes are especially suitable for preserving long-term continuous and high temporal resolution (seasonal to annual) records since their morphology (e.g., small diameter and deep depth)favors the rapid and continuous accumulation of sediments (e.g., Zolitschka 1992; Ortiz et al. 2013). The annually laminated sediments of maar lakes, like ice cores and tree rings, can provide an independently reliable chronology for palaeoenvironmental reconstruction (e.g.,Brauer et al. 1999; Chu et al. 2008). Finally, their sediments are of high value for multi-proxy studies, including palaeontological indicators (e.g., pollen and diatoms)and geophysical (e.g., magnetic properties) and geochemical parameters (e.g., major and trace elements,stable isotopes of oxygen and carbon), of palaeoclimatic and palaeoenvironmental changes (Chu and Liu 2018;Wu et al. 2019).

The volcanic field distributed from Leizhou Peninsula to the northern part of Hainan Island is the largest Quaternary basaltic province in southern China (Ho et al.2000). In respect of tectonic evolution, the formation of this volcanic field is a result of regional faults caused by Cenozoic opening of the South China Sea and the upwelling of the Hainan Plume (Yu et al. 2018). There are more than 150 volcanoes, including stratovolcanoes,shield volcanoes and maar-diatreme volcanoes (Huang et al. 1993). In particular, there are several known maars with large average diameter, such as Tianyang (> 3 km),Qingtongyang (> 2 km), Jiudouyang (> 2 km), and Huguangyan (> 2 km). According to the recently compiled Maar Volcano Location and Shape (MaarVLS)database, maars with > 2 km diameter only account for 12% of the total (29 out of 239), most of which are located in mid- to high-latitudes (Graettinger 2018). In recent decades, a number of important advances in palaeoenvironmental understanding have been obtained from analyses of the maar lacustrine sediments in Leizhou Peninsula and Hainan Island (e.g., Zheng and Lei 1999; Liu et al. 2000; Zheng et al. 2003; Yancheva et al.2007; Chu et al. 2017). However, the eruptive history and lake-peat evolution of the large maars in this lowlatitude coastal area have received little attention (Sun and Fan 2005), and there is no information on the thickness and spatial distributions of the post-volcanic lacustrine sediments in many dry maar lakes.

Thus, in the study presented here, the Jiudouyang(JDY) maar in the middle of Leizhou Peninsula was subjected to electrical resistivity tomography (ERT) imaging,and a core was drilled in the centre of the maar crater.Unfortunately, little geological information about eruptive products is available because of extreme weathering in tropical hot and humid climate. However, the ERT data and the core lithology provide accurate indications of eruptive history and lake-peat evolution of the maar volcano, and the thickness and spatial distributions of the post-volcanic lacustrine sediments. Moreover,the results may also extend a general understanding of the evolution of maar-diatreme volcanoes, superficial lacustrine sediments and Quaternary palaeoenvironments.

2 Geological setting

Volcanic rocks, including basaltic and pyroclastic rocks covers a large area around the South China Sea (Fig. 1a).According to40Ar/39Ar and K/Ar age dating results, the major volcanic eruptions began in the Oligocene, peaked in the middle Pleistocene, and finally ended in the Holocene (Huang et al. 1993; Ho et al. 2000; Yu et al. 2018).The geochemical composition of the basalts suggests that they consist of a mixture of upwelling mantle materials and terrestrial sediments, possibly originating from the opening of the South China Sea associated with lithospheric stretching followed by magma rising and consequent volcanic activities (Tu et al. 1991). An alternative explanation, based on geophysical evidences, involving the action of a mantle plume beneath Leizhou-Hainan region, has also been proposed (Lebedev et al.2000; Wei and Chen 2016). The Hainan Plume contributed to continental breakup and seafloor spreading since ca. 23.8 Ma (Yu et al. 2018). In addition, isotope studies of Holocene basalts suggest that the magma upwelling from the lower mantle was slow (< 1 cm/year), indicating that the Hainan Plume was a dying plume (Zou and Fan 2010). The tectonic evolution of Leizhou Peninsula was mostly controlled by NW-SE and W-E trending faults, and the volcanism was aligned with these faults(Fig. 1a) (Flower et al. 1992; Huang et al. 1993; Huang and Chai 1994). The dominant eruptive types in Leizhou-Hainan region included Strombolian, Hawaiian and phreatomagmatic eruption(Huang et al. 1993).

The Jiudouyang (JDY) maar is situated in the middle of Leizhou Peninsula (Fig. 1b). The altitude of the present-day crater floor is 58 m a.s.l. (above sea level),and the average height of the ejecta ring is ca. 20 m. The ejecta ring is asymmetrical, being lowest in the north,and the internal slopes are steeper in the east than in the west (Fig. 2). The elliptical crater has been infilled with lacustrine sediments and covers an area of ca. 4.5 km2, with major and minor axes of ca. 3.0 km and ca.1.8 km, respectively (Fig. 2). According to the classification scheme of Gevrek and Kazanci (2000), the JDY maar is “l(fā)arge” (diameter > 1000 m). However, we have less knowledge of the JDY maar than other large maars(e.g., Tianyang and Huguangyan maars) in Leizhou Peninsula, and the age of this maar-diatreme volcano is poorly constrained. Liu (1999) inferred that it erupted in the middle Pleistocene, based on the contact relationship of basalts and country rocks. Recently, two conflicting K-Ar ages (0.81 Ma and 1.08 Ma, respectively) were reported from basalts at ca. 50 m depth of the crater(Wang and Lu 2019). The country rocks are dominated by Quaternary sedimentary deposits, including the Zhanjiang Formation (littoral clay layers, ca. 0.8 Ma) and the Beihai Formation (alluvial sand-gravel layers, ca. 0.55 Ma) (Huang et al. 1993). Nevertheless, the JDY maar crater potentially accumulated lacustrine sediments covering several hundred thousand years and provides an important terrestrial palaeoclimatic record on glacialinterglacial timescale from tropical East Asia. In the past years, two earthquakes occurred near the JDY maar, with Richter magnitudes of 3.1 and 3.6, respectively,highlighting a need for greater understanding of the neo-tectonic activity in the region.

Fig. 1 a Tectonic sketch map of volcanic areas around the South China Sea (after Yu et al. 2018); b Simple map of Leizhou Peninsula, showing locations of the Jiudouyang maar(red star)and other large maars(red circles)

Fig. 2 Orthophoto map of the Jiudouyang maar, showing locations of the JDY15B core (red star) and the two crossed profiles (white lines). The yellow shade indicates the present-day ejecta ring,which has an artificial outlet in northeast

3 Methods

3.1 Electrical resistivity tomography (ERT) data acquisition

Multi-electrode ERT was applied to obtain high-resolution 2-D resistivity data on the JDY maar in January 2015. Data were acquired with a new WGMD-9 Multi-electrode Resistivity Imaging System produced by the Benteng Digital Control Technical Institute(Chongqing, China). We used a Wenner array with 120 electrodes (10 m spacing) connected to an 1190 m-long cable. This array provides a good signal-tonoise ratio and high sensitivity to vertical resistivity variations with smaller acquisition times (Dahlin and Zhou 2004; Martorana et al. 2009). Thus,it offered opportunities for a relatively precise correlation of the ERT results with core stratigraphy and relatively accurate determination of the thickness of crater-filling units. Two profiles were obtained, one along the minor axis (Profile 1)and the other along the major axis(Profile 2)of the crater,with five segments in total (Fig. 2; Table 1). However, the ejecta ring was not covered by these profiles due to the obstruction by the Jiudouyang village and railway. In our survey, the contact resistances were kept low (maximumvalues < 5 kΩ) and the current intensity varied between 1586 mA and 131 mA.The maximum investigation depth was approximately 160-200 m, depending on the electrode spacing and configuration(Dahlin and Zhou 2004).

Table 1 Details of the electrical resistivity tomography profiles

The resistivity data were inverted with RES2DINV software (Loke 2002), which applies a model based on the smoothness-constrained least-squares method that allows adjustment of the damping factor (deGroot-Hedlin and Constable 1990; Sasaki 1992). All except one of the inversions were supported by < 10% RMS (root mean square) error (indicating the difference between calculated and measured apparent resistivity) after five iterations (Table 1). The segments were concatenated to obtain longer profiles. However, the segment lm20150201a had a high RMS error (18.3%), owing to superficial strong resistivity contrasts resulting from road construction. Hence, this segment was excluded from the concatenation.

3.2 Core drilling and laboratory measurements

The intersection of the two profiles (20°38′5.64″N,110°2′15.72″E, 58 m a.s.l.) was chosen as the drilling location to optimize the correlation with both profiles(Fig. 2), and a core (designated as JDY15B) was obtained at this point by rotary drilling in June 2015. The core was 49.49 m long with a mean core recovery of 98%,penetrating the entire lacustrine sediments and reaching the fresh volcanic basalts. In the laboratory, the core was split lengthwise into two halves, then visually examined and photographed. Subsamples of the entire core were subsequently obtained by pushing plastic boxes (2 × 2 ×2 cm) into its surface for measurement of magnetic volume-specific susceptibility (κ) using a BartingtonInstrument MS2E sensor in combination with a MS2 control unit. Magnetic susceptibility (κ) is an indicator of ferromagnetic mineral concentration and provides a valuable tool for characterizing lithological changes. To estimate productivity, 94 samples at 50 cm interval were ground to powder and were acid-washed with 10% HCl to remove possible inorganic carbon, and then their total organic carbon (TOC) was measured using a Vairo EL cube elemental analyzer (Elementar, Germany). Six selected samples (Table 2) were measured for major element concentrations using a PANalytical Magix PW2424 X-ray fluorescence (XRF) spectrometer (precision < 5%).These measurements were all finished in the Sun Yatsen University, excepting that major elements were measured in the ALS Minerals Laboratory of Guangzhou.The geological data from the core provided valuable information that facilitated the interpretation of post-eruptive crater-filling processes and ERT results.

Table 2 Major element concentrations (in wt.%) of selected samples from the JDY15B core

4 Results

4.1 Electrical resistivity tomography (ERT)

Both profiles indicate the presence of 40-50 m thick superficial horizontal layers of low resistivity (5-60 Ω·m), despite higher resistivity near the surface due to cultivation and road construction (Fig. 3). The resistivity values are lowest (<10 Ω·m) at depths of 10-20 m. Under these conductive layers, the inversion models indicate that resistivity rapidly increases to more than 100 Ω·m (up to 1000 Ω·m), although the superficial conductive layers extend deeper into the underlying resistive unit in some places. The ERT measurements also indicate that the thickness of superficial conductive layers varies little inside the crater, at least within the constraints imposed by the spatial resolution of ERT methodology.

Profile 1 shows a thick high resistivity layer (>100 Ω·m)under the superficial conductive layers, extending across almost all of the crater’s minor axis, but split by a conductive column,ca.400 m to 500 m from the NWW starting point (Fig. 3a). The resistivity of this layer tends to be lower in the last 40 m, and in the SEE flank it is thinner and overlies conductive layers(<60 Ω·m) below depths of ca.120 m(Fig.3a).

In profile 2 (Fig. 3b), the resistivity model shows a similar overall distribution of resistivity to profile 1, but with some distinct differences. The most striking feature is a larger resistive body (from 200 to 1000 Ω·m) at approximately 45 m depth. The lower part of this large resistive body is ca. 500 m wide (from 800 to 1300 m distance) and vertically rooted beyond the maximum depth of investigation. In the upper part, this resistive body rapidly widens,and extends across nearly the entire crater between 40 m and 120 m depth, but with a concave upper limit in the centre. The rest of profile 2,below 120 m, is occupied by low-slope conductive layers(<60 Ω·m).

4.2 Lithology of the JDY15B core

The JDY15B core extends to 49.49 m, which consists of 47.92 m of lacustrine deposits overlying 1.57 m of pyroclastics and unweathered basaltic lava flow. The lacustrine deposits consist of abundant diatom and clay minerals. The overall succession can be divided into seven units (designated Unit 1 to Unit 7) from the bottom to the top based on the lithology as follows(Fig. 4):

Fig. 3 Electrical resistivity tomography (ERT) models of a Profile 1 and b Profile 2. In Profile 2, the lm20150201a segment is not directly concatenated to the profile to avoid overestimation of resistivity. The location of the JDY15B core is also shown. RMS error means the difference between calculated and measured apparent resistivity

Unit 1 (47.92-49.49 m): fresh basaltic lava flow (Figs.4, 5a) overlain by a thin layer (ca. 30 cm) of pyroclastics.κ values increase abruptly (Fig. 4). Based on the major elements of bulk sample (Table 2), the basaltic lava belongs to quartz tholeiites (Cross et al. 1902), which is one of major types of the Upper Cenozoic basalts in Leizhou Peninsula and Hainan Island(Ho et al. 2000).

Unit 2 (39.95-47.92 m): dark-greenish diatomite laminae (Figs. 4, 5a). The thickness of the lamination generally varies between sub-millimeter to millimeter scales(Fig. 5a, d). κ and TOC values remain stable, i.e., κ is generally between 0.5×10?5SI and 4×10?5SI and TOC is generally between 2% and 8% (Fig. 4). This unit is diagnosed as deep-lake facies.

Unit 3 (38.89-39.95 m): dark-greenish grey diatomaceous clay containing abundant pyroclastics (fresh without obvious weathering) (Figs. 4, 5a). The volcanic scoriae (Fig. 5c) are characterized by low K2O (0.75%)and high SiO2(51.22%) (Table 2). κ values are higher but TOC values are lower than adjacent units (diatomite laminae; Fig. 4), due to abundant pyroclastics. This unit is also diagnosed as deep-lake facies.

Unit 4 (11.99-38.89 m): dark-greenish grey diatomaceous clay to diatomite, with some varve-like laminae(27.17-30.24 m; Figs. 4, 5a). κ and TOC values remain relatively stable, similar to Unit 2. The laminae layer has relatively lower κ value. This unit is also diagnosed as deep-lake facies.

Unit 5 (10.66-11.99 m): water-saturated peat with wood fragments and organic detritus (Figs. 4, 5a, b),characterized by the lowest κ (0 SI) and the highest TOC (up to 30%) of the entire core. This unit is diagnosed as swamp facies.

Unit 6 (6.18-10.66 m): green-greyish clay minerals(Figs. 4, 5a). κ values generally increase with the depth from ca. 3×10?5SI to ca. 10×10?5SI. TOC values are low (average 0.2%). Concentrations of MgO (1.43%),CaO (0.57%), K2O (0.44%) and Na2O (0.13%) are lower than in the fresh basalt (Unit 1; Table 2). The value of chemical index of alteration (CIA) is 93 (Table 2), indicating extreme chemical weather intensity. This unit is diagnosed as shallow-lake facies.

Fig.4 Lithology,magnetic susceptibility(κ)and total organic carbon(TOC)of the JDY15B core

Fig. 5 Photograph of the JDY15 core. a Slices of lithological Units 1-7, showing with the corresponding depth of each slice in the bottom; b Wood fragments in Unit 5 at 11.35 m depth; c Volcanic scoria in Unit 3 at 39.94 m depth; d Diatomite laminae in Unit 2 at 46.39 m depth, varying between sub-millimeter to millimeter scales

Unit 7 (0-6.18 m): reddish yellow clay minerals and fine quartz and feldspar sands, with occasionally irregular and poorly rounded gravel at 1.96-3.07 m depth(Figs. 4, 5a). Elements of Mg, Ca, K and Na are also extremely depleted. The CIA value is greater than 98 (99 in Table 2), indicating extreme chemical weather intensity.Kaolinite is the dominant clay mineral in regional maar lake sediments under extremely tropical weathering(Zhou et al. 2009; Wang and Lu 2019).This unit has the highest magnetic susceptibility (κ) values (up to 26 × 10?5SI) of the entire core, while TOC values are close to 0%. This unit is diagnosed as alluvial facies.

5 Discussion

5.1 Post-eruptive evolution history

From the facies analysis of the maars in the Eifel area(Germany), Pirrung et al. (2003) developed a classification scheme for sedimentary facies division, including five lithozones representing typical stages of maar development in a maar-diatreme structure. These zones are:A(diatreme breccias), B (non-bedded syn?/post-eruptive collapse breccia), C (graded debris flow deposits and turbidites), D (lacustrine laminated sediments with turbidites and debris flow deposits) and E (swamp or peat deposits). Our core from the centre of the JDY crater included sediments almost completely covering post-eruptive development phases, from deep lake to filled basin (Fig. 4). After the eruptions, the crater was rapidly filled with water from meteoric precipitation and/or groundwater because of heavy summer rainfall.On the basis of lithological study, three major stages of the maar post-eruptive evolution can be drawn according to the changes of sedimentary phase.

The first stage of the JDY crater is the transformation into a deep-water maar lake(from Unit 2 to Unit 4) after the eruptions. The big and deep maar lake was a suitable environment for diatom colonization. The most frequent diatom species are Melosira, Stephanodiscus, Navicula and Cyclotella in the lacustrine sediments (Wang et al.1995). The chemical composition of bulk sample in this stage shows a relatively high SiO2(ca. 66 wt.%), which is a result of abundant diatoms. The geochemical analysis of most diatom shells shows that SiO2content is more than 90%, with minor Fe2O3, TiO2, etc. (Feng et al.1995). Major element compositions of the volcanic scoriae at 39.9 m depth are identical to the basaltic lava flow at 48.6 m depth and the Late Cenozoic basalts in Leizhou Peninsula (Huang et al. 1993; Ho et al. 2000; Wang and Lu 2019), characterized by low K2O, low P2O5and high SiO2(Table 2). This implies that the pyroclastics randomly found in the lacustrine sediments (Unit 3) are derived primarily from nearby eruption rather than a distal source. Variations of κ and TOC are small in magnitude (Fig. 4), which is a function of the highly internal productivity of the deep lake environment. During this period, the sediments were primarily controlled by internal lake processes and characterized by high accumulation of autochthonous biogenic silicates (diatom) and organic matter, while allochthonous materials from the crater rim and dust were secondary.

The second stage was a period close to filling up when the lake accommodation for sedimentation reduced and the water level shoaled. Consequently, the lake changed into the shallow condition and finally became a swamp environment. The corresponding sediment layer (Unit 5)is characterized by a high TOC, abundant wood fragments (Fig. 5b) and the extremely low κ (Fig. 4). Organic matter diagenesis may induce the dissolution of ironbearing minerals via iron reduction and consequently,lower κ values of sediments (Sheng et al. 2019). Few diatoms in this stage possibly indicate a high humic acid content. It must be noted that post-eruptive sedimentation does not simply end with swamp or peat deposits,as suggested by Pirrung et al. (2003). Sedimentary environment changes within a maar crater are generally controlled by the crater morphology and water level. The latter is largely affected by regional hydroclimate changes, i.e., the balance between precipitation and evaporation. Therefore, the crater potentially experiences shallow-lake-swamp-alluvial stages repeatedly associated with regional hydroclimate changes, lower water level during the drier glacial and higher water level during the wetter interglacial, till complete filling. Similarly, the Tianyang maar, another large maar, has experienced at least three sedimentary cycles in Leizhou Peninsula (Wang et al. 2012). The present-day annual groundwater table fluctuation can be as high as 3 m(Zhang et al. 2016). Furthermore, Leizhou Peninsula is located on the frontier of the East Asian summer monsoon, the intensity of which substantially varied during the past glacial-interglacial cycles (Beck et al. 2018) and consequently could influence the water depth and sedimentary environment of the maar craters.

The third stage in our study site from the depth 10.66 m upward shows a high content of clay minerals and sands with much higher magnetic susceptibility (Unit 6-Unit 7; Fig. 4). The depletion of Mg, Ca, K and Na elements and the high CIA (93-99) indicate a strong chemical weathering condition (Table 2). The sedimentary environment becomes an intermittent shallow lake.The uppermost red clay with gravel is derived from the intermittent local alluvial deposits, in which kaolinite is the dominant clay mineral.

5.2 Correlation between the ERT results and core stratigraphy

The resistivity of sediments and rocks depends primarily on the resistivity and content of the pore fluid and the nature of the pore space (Sasaki 1992). Alterations of different lithological units with the depth make ERT suitable to investigate the volcanic structure and sediment sequence of the Jiudouyang maar.The low resistivity of the superficial layers in the two profiles (Fig. 3) is clearly indicative of lacustrine or swamp/peatland sediments, and the JDY15B core provides details of their corresponding lithological composition (Fig. 6). Unit 7 has resistivity values of 20-60 Ω·m, and consists of clay and fine sand with occasional gravel. Unit 6 is composed largely of clay minerals, with 10-20 Ω·m of resistivity.Unit 5 has the lowest resistivity values (< 10 Ω·m),resulting from the high content of water in peat deposits. The diatomaceous clay or diatomite of Unit 2-Unit 4 has relatively large resistivity (10-100 Ω·m), depending on pore water content and degree of compaction. Between depths of 40-50 m, the resistivity values rapidly increase from 60 Ω·m to 100 Ω·m,and the resistivity of the basaltic lava flow (Unit 1) exceeds 100 (up to 1000) Ω·m.

Generally, the resistivity values of profiles 1 and 2 are both consistent well with the lithology of the JDY15B core (Fig. 6), which validates data acquisition and inversion. Interestingly,the peat deposits(Unit 5) and basaltic lava flow (Unit 1) can be clearly discerned in the ERT imaging profiles (Fig. 6). Similar results have been obtained from analyses of maar lake sediments in the Baruth maar and Messel Pit in Germany (Schulz et al.2005), and the Mytina maar in the Czech Republic(Flechsig et al. 2015). These findings suggest that ERT surveys can elucidate the spatial distribution of posteruptive lacustrine sediments in maar craters and help to determine depocentres of drilling cores for palaeoenvironmental research, despite some inevitable error due to the limited spatial resolution of the ERT measurements (Dahlin and Zhou 2004). Our ERT survey of the present JDY crater indicates that the thickness of the post-eruptive sediments is relatively constant (40-50 m)across the crater (Fig. 3).

5.3 Subsurface structures and eruptive history

Fig.6 Correlation between the JDY15B core lithology and the ERT results,revealing coherent lithological characteristics

The resistive layer (> 100 Ω·m and up to 1000 Ω·m) detected below 50 m depth correlates with lava flow directly overlain by the post-eruptive deep lake sediments(diatomite) (Fig. 6). This demonstrates that the eruption was complex rather than a simple type, ending with a Strombolian type and/or lava flow after an initial phreatomagmatic eruption. The ERT results indicate that the lava flow was ca. 70 m thick (at ca. 50-120 m depth). In contrast, the infill following a simple phreatomagmatic eruption is often mainly composed of thick syn?/posteruptive diatreme breccia and subsequently lacustrine deposits (White and Ross 2011; Kereszturi et al. 2014).A previous study suggested that syn?/post-eruptive diatreme breccia are poorly distinguished from lacustrine deposits by ERT measurements (Flechsig et al. 2015).

Alternation of eruption styles is common in monogenetic eruptions, which usually start with a phreatomagmatic eruption and then change into magmatic explosive or effusive activity (Kereszturi et al. 2014). Variations in eruption styles occur due to changes in the balance between magma and underground water through exhaustion of water and/or increases in the rising magma.These changes may be driven by both internal (e.g.,magma flux and ascent rate) and external factors (e.g.,regional tectonic setting and substrate conditions) (Kereszturi et al. 2011; Kshirsagar et al. 2016). In the JDY maar, it is difficult to obtain reliable details about the eruptive products, owing to extreme weathering. According to Huang et al. (1993), the pre-eruptive substrates were dominated by unconsolidated Quaternary sediments, including the Zhanjiang Formation (littoral clay layers) and/or the Beihai Formation (alluvial sandgravel layers), which were saturated with groundwater and were favorable for the onset of the phreatomagmatic eruption (Fig. 7a). However, since volcanic activities peaked in the middle Pleistocene (Huang et al. 1993; Ho et al. 2000), the magma flux and ascent rate of the Hainan Plume were likely high in this period, and eventually, the unconsolidated substrates would not be able to supply sufficient water to fuel phreatomagmatic eruption further. Therefore, the style of eruption shifted to the Strombolian type when groundwater was exhausted and/or the water/magma ratio decreased (Fig. 7c). Toward the end, the system shifted to ‘dry’ effusive activity,resulting in lava flow occupying almost the whole maar crater producing a lava lake (Fig. 7d).

An alternative explanation of the infilling sequence indicated by the core is a basaltic flood eruption (i.e., lava flow) followed by a phreatomagmatic eruption. We consider that this is not the situation of the Jiudouyang maar. As the JDY15B core revealed,the basaltic lava flow was overlain by a thin layer (ca. 30 cm) of pyroclastics and then deep-lake sediments (diatomite). Phreatomagmatic eruption(s) should produce thicker proximal pyroclastic layer than 30 cm, considering the large diameter of the JDY maar. Instead, the thin pyroclastic layer was derived from syn?/post-eruptive fall debris from crater walls. However, a possibility of external origin of lava flow could not be excluded. Since volcanic activities peaked in the middle Pleistocene in Leizhou Peninsula(Huang et al. 1993; Ho et al. 2000), there might be several types of volcanic eruptions (e.g., phreatomagmatic eruption, Strombolian eruption and lava flood eruption)in a small region. In other words, lava flow in the JDY crater might originate from a nearby eruption rather than itself. To clarify this possibility, however, detailed geological and geophysical investigation covering a larger area is needed in further research. Based on available information from ERT results and JDY core lithology, we prefer to alter eruption styles from phreatomagmatic eruption to magmatic explosive or effusive activity in the Jiudouyang maar.

Fig. 7 Eruptive history of the Jiudouyang maar in the middle Leizhou Peninsula: a Onset of phreatomagmatic eruption; b Lateral vent migration;c Strombolian eruption;d Effusive activity(lava flow)

The ERT results indicate that the floor of the initial maar crater is at depth of ca. 120 m (Fig. 3). The JDY maar is 2.4 km in average diameter and the present height of its ejecta ring is ca. 20 m. Thus, it now has a diameter/depth ratio of ca. 17, significantly larger than younger maars (ca. 70-10 ka) with diameter/depth ratios ranging from 3 to 7, according to Ross et al. (2011). We suggest that the large diameter/depth ratio of the older crater of the JDY maar (dating from the middle Pleistocene) may be mainly due to widening and shallowing of the crater ring, resulting from a long-term erosion in the tropical climate, as manifested in present-day surface processes. On the other hand, given its elongated shape(major axis to minor axis = 0.6), its large diameter/depth ratio may also be partly due to lateral vent migrations(Fig. 7b), although the deeper diatreme structures could not be determined in our ERT investigation. Analog experiments suggest that craters with one explosion epicentre tend to be circular and likely reach 560 m in diameter after tens of explosions and associated collapse,while large and complex craters require lateral vent migration (Sonder et al. 2015; Valentine et al. 2015). Furthermore, multiple coalesced maar craters are also present in Huguangyan and Qingtongyang maars in Leizhou Peninsula, indicating that lateral vent migration commonly occurred during the formation of large maars in the region,in accordance with observations of many large maars around the world (e.g., Németh et al. 2001; Auer et al.2007;Jordan et al.2013;Amin and Valentine 2017).

6 Conclusions

The combination of geophysical (ERT) technique and geological (stratigraphy and lithology) data analysis allowed us to outline the volcanic structure and spatial distribution of post-eruptive deposits in Jiudouyang maar lake, which is one of the largest maar-diatreme volcanoes in Leizhou Peninsula. Post-eruptive evolution of the maar crater had three major stages: (i) deep lake environment, characterized by high autogenetic diatom productivity; (ii) shallow lake to swamp environment,with very low water levels characterized by very high TOC values and abundant wood fragments; and, (iii)intermittent shallow lake and alluvial deposits, mainly composed of clay minerals. The correlation between ERT results and sediment lithology revealed a coherent lithological succession. The organic-rich peat layer had the lowest resistivity values (<10 Ω·m). The lacustrine silty sediments had higher resistivity values but generally less than 100 Ω·m, while the resistivity values of basaltic lava exceeded 100 Ω·m (up to 1000 Ω·m). The estimated thickness of the post-eruptive sediment overlying the basaltic lave flow was ca.50 m.

The style of eruption shifted to Strombolian or‘dry’effusive activity due to high magma flux and ascent rate,producing a ca. 70 m-thick basaltic lava flow. An initial phreatomagmatic eruption emplaced in unconsolidated Quaternary sediments produced a bowl-like crater, with a floor at a depth of ca. 120 m. The present-day large diameter to depth ratio and elongated shape of the crater implied possible lateral vent migration during phreatomagmatic eruption and long-term erosion.

Our research extended understandings of the eruptive history and evolution of large maars, and the spatial distribution of post-eruptive sediments in Leizhou Peninsula, and provided new insights into studies of the Cenozoic volcanic activity and the Quaternary palaeoenvironment changes.

Acknowledgements

We thank Hui-Wei Zheng, Mei-Ling Man, Shi-Fu Li and Hua-Cong Huang in School of Earth Sciences and Engineering, Sun Yat-sen University for their help in the fieldwork. We also thank the local people for their selfless help during the fieldwork. We are grateful to the Editor-in-Chief Prof. Zeng-Zhao Feng and two anonymous reviewers for their valuable comments and suggestions, which have substantially improved the quality of the original manuscript.

Authors’ contributions

ZZ, KYH, and CC designed the research; CC and ZZ prepared the manuscript,which was edited by all of the co-authors; CC, LFZ, and LPT did the fieldwork and laboratory analysis; FX provided helpful and constructive advice in the fieldwork and data analysis. All authors read and approved the final paper.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2016 YFA0600500), the National Natural Science Foundation of China(Grant Nos. 41472143 and 41661144003), and the Guangdong Province Introduction of Innovative R&D Team of Geological Processes and Natural Disasters around the South China Sea (Grant No. 2016ZT06N331).

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Received: 5 June 2020 Accepted: 2 December 2020