Home Petrogenesis of the Sukadana Basalt based on petrology and whole rock geochemistry, Lampung, Indonesia: Geodynamic significances
Article Open Access

Petrogenesis of the Sukadana Basalt based on petrology and whole rock geochemistry, Lampung, Indonesia: Geodynamic significances

  • Luhut Pardamean Siringoringo EMAIL logo , Benyamin Sapiie , Alfend Rudyawan and I Gusti Bagus Eddy Sucipta
Published/Copyright: November 15, 2023
Become an author with De Gruyter Brill

Abstract

The petrogenesis of Sukadana Basalt remains an enigma till present. Major and trace element data are analyzed from Sukadana Basalt lava, located at East Lampung, Sumatra, to study the processes involved in the petrogenesis of the erupted magmas and the origin of mantle source compositions. The Sukadana Basalt display SiO2 (48.1–52.5 wt%), MgO (5.3–9.3 wt%), TiO2 (1.3–2.6 wt%), P2O5 (0.2–0.6 wt%), and Fe2O3T (8.9–11.3 wt%) contents. The Sukadana Basalt enriched in light rare earth elements with weak negative Eu anomalies (Eu/Eu* = 0.8–1) and show Ocean Island Basalt (OIB)-like characteristics. There are two different petrogenesis groups, namely group A and group B. Group A samples show enrichment of Hf, Pb, K, and Sr and depletion of Nb. Group B samples show enrichment of K, Sr, and depletion of Pb. These differences are closely related to the mechanism of slab roll-back and normal fault activity. This study shows that Sukadana Basalt has Nb = 7.4–29.8 ppm, Nb/U = 18–060.3, and Nb/La = 0.8–1.6. These characteristics were similar to those found in typical Nb-enriched basalts. Geochemical analyses suggest that the Sukadana Basalt have experienced minimal crustal contamination and Olivine plagioclase, clinopyroxene, and magnetite fractional crystallization. The chemical features, together with high incompatible-element ratios, are consistent with low degrees of partial melting of a dominantly a partial melting of garnet–peridotite mantle source. The trace-element patterns suggest a mantle source influenced by an enriched component. The occurrence of OIB-like basalt suggests significant upwelling of the asthenosphere in response to slab roll-back. These processes occured in the above of a Paleo Indo-Australia subducting N–S beneath the southern part of Sumatra.

Graphical abstract

1 Introduction

Indonesia, located in the Sundaland Block, has the highest number of active volcanoes among other countries in the region [1]. Sundaland is the continental extension of the Eurasian plate in Southeast Asia [2]. These volcanoes are formed by the subduction of the Indian, Australian, Pacific, and Philippine oceanic plates beneath Southeast Asia [3,4,5]. The accumulation of magma on the surface creates volcanoes, with magma produced from the partial melting of the oceanic plate at depths of 150–300 km [6]. The subduction of oceanic plates forms the convergent subduction mechanism, including the Andaman–Sumatra–Java trench in the west and south and the Philippine Trench in the east [7]. In Sumatra Island, the subduction mechanism creates the Bukit Barisan volcanic chain.

Unlike the Bukit Barisan, the Sukadana Basalt is an anomalous volcano. The anomaly encompasses a very wide lava distribution, location in the back-arc, and its composition dominated by mafic minerals (Figure 1(a)). According to Dosso et al. [8], its affinity is between Mid Oceanic Ridge Basalt (MORB) and Oceanic Island Basalt (OIB). However, Nishimura et al. [9] and Gasparon [10] suggested that its magma affinity is MORB, while Romeur [11] proposed that it is Enriched Mid Oceanic Ridge Basalt (E-MORB). These differences highlight the unresolved nature of its petrogenesis.

Figure 1 
               (a) The southern of Sumatra Island with location of the Sukadana Basalt (red rectangle). (b) The cinder cones distribution at the Sukadana Basalt with sampling location of this study and Gasparon’s study [10]. (c) The geology of Basalt Sukadana and its surroundings, as outlined in the Tanjungkarang quadrangle regional geology [22].
Figure 1

(a) The southern of Sumatra Island with location of the Sukadana Basalt (red rectangle). (b) The cinder cones distribution at the Sukadana Basalt with sampling location of this study and Gasparon’s study [10]. (c) The geology of Basalt Sukadana and its surroundings, as outlined in the Tanjungkarang quadrangle regional geology [22].

This study presents new geochemical analysis (major and trace element) and petrographic data from 13 samples of Sukadana Basalt (Figure 1(b)) to better understand the degree of alteration and metamorphism, identify the magma processes (such as fractional crystallization and crustal contamination) and determine the mantle source of the basalt, and develop a new tectonic model for the geodynamic evolution of southern Sumatra. This study is very important in understanding the petrogenesis characteristics of basalt that appear in the back arc of the Eurasian plate and the correlation of the appearance of basalt in the back-arc with the tectonic complexity of the Sunda Strait, which is a transitional zone of oblique and orthogonal subduction.

2 Geology background and petrography

The Indo-Australian Oceanic Plate subducts beneath the Eurasian continental plate at a rate of 6–7 cm/year [12,13], forming a 5,600 km magmatism arc from the Andaman Islands to the Banda Arc. The Sunda Strait, a transitional zone between the Sumatra oblique and Java perpendicular subductions, is of particular interest [14,15,16]. It is related to an extensional feature caused by either the northwestward displacement of the southern Sumatra block along the Sumatra Fault [15] or the rotation of Sumatra relative to Java in the Late Cenozoic [17]. The age of the oceanic lithosphere beneath Sumatra is less than 100 million years old, while the oceanic lithosphere beneath Java and the eastern Sunda is older [18]. A volcanic configuration across the Sunda Strait spans from Panaitan Island through Krakatau Island, the Sebesi and Sebuku Islands, the Rajabasa Volcano, to the Sukadana Basalt plateau from SSW to NNE, according to Harjono et al. [19] and Nishimura et al. [9] (Figure 1(a)). The Sukadana Basalt and Bukit Telor Basalt are located in the back-arc region of Sumatra [10], but the existence of this volcanic lineament is still debated, especially its influence on the presence of the Sukadana Basalt.

The age of Sukadana Basalt is quaternary. K–Ar method show that the age of Sukadana Basalt is 0.8 Ma [9], 1.2 Ma [20], and 0.5–1.3 Ma [11]. This age suggests that the presence of the Sukadana Basalt is not much different from the eruption of Mount Rajabasa, the nearest mount to Sukadana Basalt. Mount Rajabasa age is 0.12–0.31 Ma [21]. Figure 1(c) shows the complexity of Lampung tectonics especially surrounding the Sukadana Basalt. The old and young rocks are exposed not far from each other. This condition makes this area is very interesting to study. The Sukadana Basalt is not the only young rock in this area. The Sukadana Basalt overlays pliocene Lampung Formation and Eocene Tarahan Formations which are emitted by centers within the volcanic arc [14,22]. The Sukadana Basalt is located on Lampung Province. The Sukadana Basalt spans an estimated area of 1,000 km2 and is comprised of numerous basaltic flows with thicknesses of up to 2–3 m. There are several hills present, which likely represent the eruptive centers, with elevations exceeding 200 m above sea level. Although the thickness of the basaltic pile can reach up to 200 m in certain areas, it varies locally [14].

The Sukadana Basalt consists of a cinder cone and plain morphology (Figure 1b). The cinder cone shows dimensions about 290–5,200 m (wide) and 111–270 m (height) above sea level. In observation point 1.2, The Sukadana Basalt shows the cave as a part of a lava tube feature with layered lava structure, dark to green weathered color, light grey fresh color, and vesicular structures (Figure 2(a) and (b)). In addition, observation point 3.1 shows cinder cone morphology without lava tube features (Figure 2(c)). The basalt shows a brown-dark grey weathered color, light grey fresh color with vesicular structures. Certain areas of rocks show vesicular structures that are filled by zeolite (Figure 2(d)). Almost all of the research area found cobble-bolder size basalt either at cinder cone or plain morphology (Figure 2(c), (e), and (g)). The bolder also shows vesicular structures with middle to high quantity of vesicles, dark weathered color, and light grey fresh color (Figure 2(f) and (h)). This vesicular structure indicates that magma contains a massive amount of gas, specifically H2O, CO2, and S [23,24].

Figure 2 
               (a) Basalt outcrops at observation point 1.2 show the cave as a part of a lava tube. (b) Layered lava structure with dark to green weathered color. (c) Observation point 3.1 shows cinder cone morphology without lava tube feature. (d) At observation point 3.1, the basalt shows a brown-dark grey weathered color, light grey fresh color with vesicular structures. Certain areas of rocks show vesicular structures that are filled by zeolite. (e) and (g) Almost all of the research area was found cobble-bolder size basalt either at cinder cone or plain morphology. (f) and (h) The bolder also shows vesicular structures with middle to high quantity of vesicles with dark weathered color.
Figure 2

(a) Basalt outcrops at observation point 1.2 show the cave as a part of a lava tube. (b) Layered lava structure with dark to green weathered color. (c) Observation point 3.1 shows cinder cone morphology without lava tube feature. (d) At observation point 3.1, the basalt shows a brown-dark grey weathered color, light grey fresh color with vesicular structures. Certain areas of rocks show vesicular structures that are filled by zeolite. (e) and (g) Almost all of the research area was found cobble-bolder size basalt either at cinder cone or plain morphology. (f) and (h) The bolder also shows vesicular structures with middle to high quantity of vesicles with dark weathered color.

Petrographically, Sukadana Basalt can be divided according to their mineralogy into sub-groups, including basalt, basaltic trachy andesite, and basaltic andesite. The basalt lithologies contain mainly plagioclase (30–52%), olivine (10–27%), volcanic glass (5–50%), pyroxene (19–30%), and iddingsite rim (3–7%) with vesicular structures. The phenocryst includes plagioclase, olivine, piroxene, and opaque mineral. The plagioclase (grain size: 0.05–1.4 mm) is euhedral–subhedral, carlsbad, and carlsbad–albit twinning, some plagioclases show sieve–glomeroporphyritic texture, and labradorite–andesin composition. The olivine (grain size: 0.05–1.5 mm) is euhedral–subhedral; some olivine shows plagioclase inclusion and iddingsite rim. The piroxene (grain size: 0.05–1.4 mm) consists of clinopiroxene (8–20%) and ortopiroxene (3–15%) with subhedral-anhedral crystals (Figure 3(a) and (b)). The opaque mineral (grain size: 0.025–0.25 mm) is subhedral-anhedral. The groundmass includes plagioclase, piroxene, and volcanic glass.

Figure 3 
               Cross-polarized photomicrographs from Sukadana Basalt. (a) and (b) The basalt lithologies contain mainly plagioclase, olivine, volcanic glass, pyroxene, and iddingsite rim with vesicular structures. Some plagioclases show sieve-glomeroporphyritic texture. The olivine shows plagioclase inclusion and iddingsite rim. (c) The basaltic trachy andesite lithologies contain mainly plagioclase, olivine, pyroxene, and iddingsite rim with vesicular structures. Some plagioclases show sieve texture. Some olivines show plagioclase inclusion and iddingsite rim. (d) The basaltic andesite lithologies contain plagioclase, olivine, pyroxene, volcanic glass, and iddingsite rim with vesicular structures. Some plagioclases display a sieve texture. Some olivines show iddingsite rim.
Figure 3

Cross-polarized photomicrographs from Sukadana Basalt. (a) and (b) The basalt lithologies contain mainly plagioclase, olivine, volcanic glass, pyroxene, and iddingsite rim with vesicular structures. Some plagioclases show sieve-glomeroporphyritic texture. The olivine shows plagioclase inclusion and iddingsite rim. (c) The basaltic trachy andesite lithologies contain mainly plagioclase, olivine, pyroxene, and iddingsite rim with vesicular structures. Some plagioclases show sieve texture. Some olivines show plagioclase inclusion and iddingsite rim. (d) The basaltic andesite lithologies contain plagioclase, olivine, pyroxene, volcanic glass, and iddingsite rim with vesicular structures. Some plagioclases display a sieve texture. Some olivines show iddingsite rim.

The basaltic trachy andesite lithologies contain mainly plagioclase (40%), olivine (20%), pyroxene (30%), and iddingsite rim (5%) with vesicular structures. The phenocryst includes plagioclase, olivine, piroxene, and opaque mineral. The plagioclase (grain size: 0.075–0.5 mm) is euhedral–subhedral, carlsbad and carlsbad–albit twinning, some plagioclases show sieve texture, and labradorite-andesin composition. The olivine (grain size: 0.25–0.5 mm) is anhedral-subhedral, some olivines show plagioclase inclusion and iddingsite rim. The piroxene (grain size: 0.05–1.4 mm) consists of clinopiroxene (17%) and ortopiroxene (13%) with subhedral–anhedral crystals (Figure 3(c)). The opaque mineral (grain size: 0.025–0.1 mm) is subhedral–anhedral. The groundmass includes olivine and piroxene.

The basaltic andesite lithologies contain mainly plagioclase (25%), olivine (8%), pyroxene (15%), volcanic glass (50%), and iddingsite rim (2%) with vesicular structures. The phenocryst includes plagioclase, olivine, piroxene, and opaque mineral. The plagioclase (grain size: 0.05–0.375 mm) is euhedral–subhedral, carlsbad, and albit twinning; some plagioclases display a sieve texture and labradorite–andesin composition. The olivine (grain size: 0.1–1.5 mm) is anhedral–euhedral; some olivines show iddingsite rim. The piroxene crystals, with a grain size: 0.075–0.125 mm, consist of clinopiroxene (8%) and ortopiroxene (7%) which are euhedral-anhedral in shape (Figure 3(d)). The groundmass includes plagioclase, piroxene, and volcanic glass.

3 Methods

Thirteen samples of fresh basalt were gathered for geochemical and petrographic analysis. Nine of these samples (69%), were obtained from the center of cones while the remaining 31%, or four samples, were gathered from plain morphologies (Figure 1(b)). The petrographic analysis was performed at the Petrology-Mineralogy Laboratory of the Geological Engineering Department of Sumatera Institute of Technology using a polarization microscope. The geochemical analysis, which included both major and trace elements, took place at the Intertek Laboratory in Jakarta, Indonesia. Samples fused using Lithium Metaborate then were detected using XRF (Malvern Panalytical brand) and recalculated to 100%, volatile free for major elements analysis. Trace elements were identified using the four acids digestion-ICP OES (Agilent brand) and MS method.

The measurement was conducted by utilizing predetermined calibration standards, and the resulting data was processed through laboratory information management system software. This software automates the necessary calculations to determine both the percentage of oxide for individual elements and the total percentage of major elements in the sample, including the Loss on Ignition (LOI). All results have passed rigorous quality control, as indicated by the tolerable deviation. For instance, SiO2 from 1.4 is 51% (Table 1). Quality control protocol used NCS DC 73322 standards, which have a value of 51.09% (lower bound: 50.797%; upper bound: 51.103%). The detection limit was varied during analysis for different elements. It was 50 ppm for Al, Ca, P, S; 5 ppm for Cr, Ti; 1 ppm for Cu, Mn, Ni, Sc, V, Zn, As, Ba, Co, Pb, Se; 0.01% for Fe, Al2O3, CaO, Cr2O3, Fe2O3, K2O, MgO, MnO, Na2O, P2O5, SiO2, TiO2; 20 ppm for K, Mg, Na; 0.5 ppm for Be, Zr, Sr; 0.05 ppm for Bi, Cd, In, Re, Ta, Th, U, Lu, Pr, Tb; 0.1 ppm for Ag, Cs, Ga, Ge, Hf, Li, Mo, Nb, Rb, Sb, Sn, Te, W, Y, Ce, Dy, Er, Eu, Gd, Ho, La, Nd, Sm, Tm, Yb; 0.002% for P2O5, S; 0.02 ppm for Tl.

Table 1

Whole-rock geochemistry data of the Sukadana Basalt

No. samples 1.1 1.2 1.3 1.4 1.5 1.6 2.1 2.2 2.3 3.1 3.2 3.3 3.4
Lithology Basalt Basalt Basaltic-andesite Basalt Basalt Basaltic trachy-Andesite Basalt Basaltic-andesite Basaltic-andesite Basalt Basalt Baaltic-andesite Basalt
Major elements (%)
SiO2 48.2 51.5 52.1 51.0 51.9 52.5 51.9 52.1 52.3 50.2 49.1 51.3 48.1
TiO2 1.5 1.4 1.4 1.3 1.4 1.8 2.3 1.4 1.4 1.5 1.8 1.5 2.6
Al2O3 16.7 16.3 16.3 16.3 15.8 17.1 15.1 16.3 16.3 15.5 15.6 16.3 16.2
Fe2O3 11 10.1 9.9 9.8 9.4 8.9 9.8 10.8 9.4 10 11.1 9.5 11.3
MnO 0.2 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
MgO 9.3 7 6.9 7.9 7.5 5.3 7.1 6.4 6.9 7.7 7.4 6.4 7.5
CaO 8.5 8.5 8 8.4 8.7 8.8 8.5 8.3 8.8 8.7 8.5 8.2 7.5
K2O 0.8 0.8 0.6 0.6 1 1.5 1.3 0.8 1 1.3 1.2 1.3 1.6
Na2O 3 3.6 3.3 3.1 3.4 3.9 3.4 3.5 3.4 3.4 3.4 3.6 3
P2O5 0.2 0.2 0.2 0.3 0.3 0.4 0.5 0.2 0.2 0.3 0.5 0.3 0.6
LOI 0.9 0 0.7 0.4 0.3 0 0.1 0.5 0.1 1.1 0.3 0.6 1.5
Total 100 99.4 99.5 99.3 100 100 100 100 100 100 99 99.3 100
Mg# 62.9 58.2 58.3 62 61.7 54.6 59.1 54.6 59.5 60.8 57.1 57.5 57.2
Trace elements (ppm)
V 187 152 178 179 167 182 189 155 179 171 169 209 196
Cr 315 162 193 284 202 50 297 177 156 270 273 245 251
Ni 188 94 108 171 136 42 114 164 101 199 195 194 203
Zn 117 105 121 113 104 102 112 111 96 106 129 107 124
Rb 6.2 9.3 7.9 12.1 17.9 21.3 19.5 10.7 18.2 22.6 15.6 17.5 25.3
Sr 488 501 486 486 595 760 672 542 584 486 539 574 624
Y 17.2 14.2 12.5 14.5 16.6 23.6 26.9 14 16.3 17.9 20.3 24.2 31.8
Zr 104 69 73 80 103 168 216 90 93 119 145 135 255
Nb 14.8 10.9 7.4 7.9 11.9 22.3 27.2 13.4 11.2 16.4 16.8 16.2 29.8
Ba 197 101 80 122 175 269 300 136 153 327 310 350 505
La 11.6 7 5.7 9.5 10.7 17.6 20.5 9.4 10.3 13.3 16.3 17.3 21.2
Pb 3 6 2 2 2 3 2 2 2 3 <1 3 2
Th 2.5 1.5 1.2 2.1 2.5 3.1 2.2 1.9 2.2 2.8 2.9 2.9 3.1
Sc 21 16 17 19 19 21 20 15 19 20 21 20 23
Co 47 41 38 41 39 32 37 41 36 49 50 44 49
Cs 0.1 0.4 0.2 0.4 0.7 0.1 0.5 0.5 0.6 0.3 0.5 0.3 0.4
Ce 24.9 16 13.6 20.4 23.6 37.8 47.7 20.2 22.7 26.7 33.8 30.8 46.7
Nd 14.6 11.4 10.5 13.4 14.2 23.3 35.7 12.2 14.3 14.5 20 19.3 34.5
Sm 3.9 3.2 3.3 3.6 3.8 5.7 8.5 3.5 3.8 3.9 5.2 4.9 8.1
Eu 1.2 1.2 1 1.1 1.1 1.7 2.2 1.1 1.1 1.1 1.3 1.3 2
Gd 3.9 3.7 3 3.6 3.8 5.6 7.4 3.4 3.7 3.9 5.1 4.9 7.4
Tb 0.6 0.5 0.5 0.5 0.5 0.8 1 0.5 0.5 0.6 0.7 0.7 1.2
Ho 0.6 0.6 0.5 0.5 0.6 0.8 1 0.5 0.6 0.6 0.8 0.7 1.2
Er 1.8 1.3 1.3 1.5 1.7 2.3 2.5 1.4 1.7 1.6 2 1.9 2.9
Tm 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.2 0.2 0.2 0.2 0.2 0.4
Yb 1.5 1.1 1 1.1 1.5 1.9 1.9 1.1 1.4 1.4 1.7 1.6 2
Lu 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.2 0.2 0.2 0.2 0.2 0.3
Hf 2.6 2.7 2.2 2.3 2.7 4.1 5.2 2.4 2.5 3.6 4 3.7 7
Ta 0.9 0.7 0.5 0.5 0.8 1.3 1.5 0.9 0.7 1.7 1.8 1.7 2.2
U 0.6 0.4 0.3 0.4 0.6 0.4 0.6 0.5 0.6 0.7 0.8 0.8 0.8
Al 85,800 83,900 85,000 85,700 83,100 89,000 79,200 80,700 82,100 82,300 81,800 89,100 86,200
Ca 59,000 59,500 56,000 59,700 61,500 63,100 59,600 58,000 61,300 63,200 62,100 62,800 53,800
Fe 7.6 6.6 7.0 6.9 6.6 6.1 6.9 7.3 6.4 6.9 7.7 6.8 7.9
K 6,380 5,910 4,580 5,280 8,550 12,500 10,800 6,400 8,280 10,800 9,770 10,900 13,700
Mg 53,500 39,900 39,000 45,800 43,300 30,200 40,000 37,000 39,800 44,000 41,900 37,400 42,500
Mn 1,130 920 924 982 1,030 1,060 966 1,180 944 1,180 1,130 967 1,030
Na 21,500 26,200 23,900 23,000 24,400 28,400 24,900 25,600 25,100 23,700 22,400 24,900 21,500
P 870 880 750 1,170 1,250 1,700 2,120 920 1,020 1,300 2,040 1,360 2,370
S <50 80 50 <50 <50 <50 <50 <50 <50 130 90 90 60
Ti 8,210 7,980 7,930 7,430 8,070 10,400 12,800 7,440 7,830 8,710 8,940 8,370 14,100
Ag <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 0.10 <0.1 0.1
As 2.00 <1 2 2 2 4 3 2 2 3 2 3 1
Bi <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05
Cd 0.3 <0.05 0.2 0.1 0.2 0.3 0.4 0.2 0.1 0.2 0.2 0.3 0.3
Ge 1.4 1.3 1.5 1.3 1.8 1.6 1.6 1.4 1.6 1.3 1.7 1.7 1.9
Ln 0.1 <0.05 0.1 <0.05 0.1 0.1 0.1 0.1 0.1 0.1 <0.05 <0.05 0.1
Mo 0.7 0.9 1.1 0.8 1.3 1.2 4.7 1.4 1.4 1.5 0.7 1.3 1.7
Re <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05 <0.05
Sb 0.2 0.2 0.1 <0.1 <0.1 0.4 0.5 0.2 0.2 0.1 <0.1 0.2 <0.1
Se <1 <1 <1 <1 <1 <1 <1 <1 <1 1 1 1 2
Sn 1.7 2 1.9 1.4 1.9 2.5 2.1 1.4 1.4 1.8 0.8 1.7 1.8
Te <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 <0.1 0.2 <0.1 <0.1 <0.1
W <0.1 0.6 0.2 0.3 2.2 0.6 1.1 <0.1 0.3 1.5 0.6 1 0.8
Li 7.5 7 5.5 6.5 7.2 6.9 6.4 5.4 6.8 5.8 5.1 7.3 5.9
Be 1 0.9 1 0.8 1.1 1.5 1.5 1.1 1 1.1 1.5 1.3 1.7
Cu 61 48 57 63 57 43 31 61 60 67 77 48 63
Zn 117 105 121 113 104 102 112 111 96 106 129 107 124
Ga 18.8 19.8 18.9 18.6 18.5 21.5 20 18.2 18.6 20 19.9 21.3 22.6
Pr 3.2 2.1 2.1 2.8 3.1 5.1 7.7 2.7 3 3.2 4 4 6.3
Dy 3.5 2.9 2.6 3 3.2 4.6 5.4 2.8 3.3 3.5 4.1 4.0 6.3

4 Result

All the major element analysis was redetermined to 100 wt% with no volatile material, as illustrated in Table 1. Most of the samples are petrographically unaltered with only slight signs of alteration, demonstrated by their generally low LOI values (<2 wt%). This value indicates that the original chemical composition of the basaltic melts that were formed is reliable [25,26]. Table 1 presents the compositions of both major and trace elements in the samples. The whole-rock SiO2 content of the basaltic samples varies between 48.1 and 52.5 wt% with low Na2O + K2O content (∼3.7–5.4 wt%). All the samples have MgO content greater than 5.2 wt%, with Mg# ranging from 54.6 to 62.9. The LOI values are between 0 and 1.5 wt% (Table 1). The rocks are classified as basalt, basaltic andesite, and basaltic trachy andesite based on their SiO2 vs total alkali (Na2O + K2O) diagram (Figure 4(a)) [27]. This result is consistent with previous studies from Romeur [11], Gasparon [10], and Zulkarnain [28]. The samples are plot in the subalkaline series on the SiO2 vs Na2O + K2O diagram, except for one which plots in the alkaline series (Figure 4(a)) [29]. Most of the samples belong to the sodic series, but one plots in the postasic series on a K2O vs Na2O discrimination diagram (Figure 4(b)) [30]. The samples have high TiO2 (1.3–2.6 wt%), moderate to high Al2O3 (15.1–17.1 wt%), Fe2O3(t) (total Fe oxide, assuming all Fe is trivalent) = 8.9–11.3 wt%, and P2O5 = 0.2–0.6 wt%.

Figure 4 
               (a) Total alkali-silica classification diagram [27] for Sukadana Basalt with integration from the previous study. The straight blue line represents the boundary between alkaline and sub-alkaline basaltic series after Irvine and Baragar [29]. (b) Na2O vs K2O diagram [30].
Figure 4

(a) Total alkali-silica classification diagram [27] for Sukadana Basalt with integration from the previous study. The straight blue line represents the boundary between alkaline and sub-alkaline basaltic series after Irvine and Baragar [29]. (b) Na2O vs K2O diagram [30].

The Nb vs Ti diagram shows that there are two highly distinct groups (Figure 5). This difference is related to the variation in petrogenesis between these two groups. Group A indicates low values of Nb and Ti, while Group B shows high values of Nb and Ti (Nb > 22; Ti > 10,000). The Rare Earth Element (REE) patterns of Groups A and B (Figure 6) are normalized to primitive mantle (PM) [31] and Chondrite (ChN) [31] identical and similar to the OIB pattern. However, upon closer identification of individual elements, it is observed that the Nb and Ti elements in Group A tend to be depleted.

Figure 5 
               The Nb vs Ti diagram shows two variations of petrogenesis.
Figure 5

The Nb vs Ti diagram shows two variations of petrogenesis.

Figure 6 
               (a) Group of samples with depleted Nb. (b) Group of rock samples with stable Nb. Both groups of REE are normalized to ChN and primitive mantle [31]. OIB, N-MORB, and E-MORB from Sun and McDonough [34]. OIB = Oceanic Island Basalts, E-MORB = Enriched Mid Oceanic Ridge Basalt, N-MORB = Normal Mid Oceanic Ridge Basalt.
Figure 6

(a) Group of samples with depleted Nb. (b) Group of rock samples with stable Nb. Both groups of REE are normalized to ChN and primitive mantle [31]. OIB, N-MORB, and E-MORB from Sun and McDonough [34]. OIB = Oceanic Island Basalts, E-MORB = Enriched Mid Oceanic Ridge Basalt, N-MORB = Normal Mid Oceanic Ridge Basalt.

The patterns are characterized by a slight enrichment of light rare earth elements (LREE) with (La/Yb)N = 3.9–7.2, and The LREE content of the sample is approximately 20–100 times chN, while the HREE content is about 8–18 times chN. There is a weakly negative Eu anomaly present, with Eu/Eu* ratios from 0.8 to 1. In general, the spider diagram of group A normalized to primitive mantle [31], the samples show enrichment of Hf, Pb, K, and Sr and depletion of Nb (Figure 6(a)). On the spider diagram of group B normalized to primitive mantle [31], the rock samples show enrichment of K, Sr, and depletion of Pb (Figure 6(b)).

5 Discussion

5.1 Crustal contamination

When hot basaltic magma passes through continental crust, certain elements from the crust can greatly impact the composition of the resulting rocks [32]. But this can be avoided if magma rises quickly through the lithosphere or crust [33].

The Th/Ta ratio in the Sukadana Basalt varies between 1.4 and 3.8, indicating insignificant influence from crustal contamination in its formation. This interpretation is based on a comparison of the ratio of the original mantle (Th/Ta = 2.3) and the continental crust (Th/Ta = 10) as reported by Sun and McDonough [34]. Moreover, Nb/La ratios (<1.0), which are used as an indicator of crustal influence on magma [35], range from 0.8 to 1.6 in the Sukadana Basalt, showing crustal contamination but not significant.

In Figure 7(a) and (b), which show Nb/Th vs Nb/La [36] and Rb/La vs Th diagrams [37] respectively, the samples in the Sukadana Basalt are close to the mantle composition and not precise in the continental crust region. This indicates insignificant contamination from the upper crust (UC). The Sukadana Basalt has a Zr/Nb ratio that ranges from 6.3 to 10.1. This ratio is higher than what is typically found in OIB mantle sources (3.2–5), but lower than the ratios found in continental crust (16.2) and primitive mantle (14.8). It is also much lower than the ratio found in normal-MORB (30) [38,39], which suggests that significant crustal contamination did not occur. In the results, the composition of Sukadana Basalt is more likely influenced by magma dynamics at its source and partial melting in the mantle region. An integration of Ce/Pb and La/Nb values can be utilized to evaluate the extent of crustal contamination and mantle mixing processes [40,41,42]. Basalts with Ce/Pb ratios in the range of ∼20–30 are improbable to have incorporated crustal material [41]. La/Nb ratios of 0.6–1.2 and Ce/Pb ratios of 2.7– > 33.8 plot in the mantle region on a La/Nb-Ce/Pb diagram (Figure 7(c)).

Figure 7 
                  (a) Nb/La vs Nb/Th diagram [36]. (b) Th vs Rb/La diagram [37]. (c) La/Nb vs Ce/Pb [40]. Values of N-MORB, E-MORB, and OIB are from Sun and McDonough [34], uncontaminated basalts zone from Rooney et al. [41]. Values of primitive mantle (PM) are from McDonough and Sun [31]. Values for the UC and lower crust (LC) are from Wedepohl [82]. LC: Lower Crust; UC: Upper Crust; N-MORB: Normal-Mid Oceanic Ridge Basalt; E-MORB: Enriched-Mid Oceanic Ridge Basalt; OIB: Oceanic Island Basalt; PM: Primitive Mantle.
Figure 7

(a) Nb/La vs Nb/Th diagram [36]. (b) Th vs Rb/La diagram [37]. (c) La/Nb vs Ce/Pb [40]. Values of N-MORB, E-MORB, and OIB are from Sun and McDonough [34], uncontaminated basalts zone from Rooney et al. [41]. Values of primitive mantle (PM) are from McDonough and Sun [31]. Values for the UC and lower crust (LC) are from Wedepohl [82]. LC: Lower Crust; UC: Upper Crust; N-MORB: Normal-Mid Oceanic Ridge Basalt; E-MORB: Enriched-Mid Oceanic Ridge Basalt; OIB: Oceanic Island Basalt; PM: Primitive Mantle.

Nevertheless, only 25% of the samples from this and previous studies plot in the contaminated mantle, characterized by low Ce/Pb ratios. The Nb/Th vs Nb/La and Rb/La vs Th diagrams reveal that the influence of crustal material on the production of high Pb content in the Sukadana Basalt is insignificant. The absence of crustal xenoliths in all samples also provides the conclusion of minimal crustal contamination. Overall, The Sukadana Basalt magmatism had moved through the Sumatra continental crust with minimal crustal contamination. This is supported by the fact that the thickness of the Sumatra crust ranges from 27 to 35 km [43] and the position of Sukadana Basalt is in the back arc. This results in relatively little contact between magma and the continental crust.

5.2 Fractional crystallization

To determine the dominant magma process, whether it is partial melting or fractional crystallization, we used the La vs La/Yb diagram (Figure 8) [44]. This diagram displays that the samples are located between the lines of partial melting and fractional crystallization, indicating that the trace elements of the Sukadana Basalt magma were controlled by both processes during its evolution. The weakly negative Eu anomalies (Eu/Eu* = 0.8–1) in the samples suggest that they likely underwent minor plagioclase fractionation, either prior to eruption or during the magma’s ascent and evolution [45]. In addition to La vs La/Yb diagram, another diagram that can explain the processes of fractional crystallization and partial melting during magma evolution is the Harker diagram [46,47]. This diagram shows the correlation between major elements. Referring to the Harker diagram, Sukadana Basalt magma underwent clinopyroxene and magnetite fractionation. This is indicated by the positive correlation between Fe2O3 and TiO2 with MgO (Figure 9). The diagram also indicates the accumulation of olivine and chromite (Fe oxides) and Ti oxide minerals. Moreover, the positive correlation between CaO and MgO (Figure 9) indicates the presence of clinopyroxene fractionation. Additionally, this correlation shows that basaltic andesite is a product of basalt fractional crystallization. The negative correlation between SiO2, Al2O3, and MgO indicates the accumulation of clinopyroxene.

Figure 8 
                  La vs La/Yb diagram [44].
Figure 8

La vs La/Yb diagram [44].

Figure 9 
                  Harker diagram.
Figure 9

Harker diagram.

The previous explanation indicates that half of the Sukadana Basalt magma process underwent fractional crystallization. This is also supported by unequal granular phenocrystic basalts consisting of pyroxene, plagioclase, and olivine phenocrysts, which suggests that the parental magma underwent fractional crystallization before eruption [48]. This is further confirmed by the differences in Mg# values between the Sukadana Basalt (54.6–62.9) and the original basaltic magma (65–70) [49] or higher than 73 [50,51]. The distinctive SiO2 (48.1–52.5 wt%), MgO (5.3–9.3 wt%), and total Fe2O3 (8.9–11.3 wt%) content also support the fractionation nature of the basalt. In addition, the relatively low Ni (42–203 ppm) and Cr (50–315 ppm) content in the Sukadana Basalt compared to primary mafic magmas (Ni > 400 ppm; Cr > 1,000 ppm) [50,51] indicates that the parental magma underwent fractional crystallization of olivine and clinopyroxene during its evolution.

5.3 Partial melting

All of the samples studied display slightly high LREE/HREE ratios, with (La/Yb)N values ranging from 3.9 to 7.2, indicating that they were derived from melts produced at low degrees of melting. In addition, the effect of partial melting is well demonstrated by a model illustrated on the Gd/Yb vs La/Yb diagram (Figure 10(a)) [52]. The curves Garnet 4% and Garnet 8% indicate the contents of garnet in the source, according to Halliday et al. [53]. The diagram shows 4–6% degree of partial melting of garnet peridotite (4–8%). La/Yb vs Dy/Yb diagram [54,55] also indicates that garnet peridotite is a source region of Sukadana Basalt besides low degrees of partial melting (16–20%) (Figure 10(b)). In extensional settings, the asthenosphere experiences decompression melting at a relatively shallow depth of less than 80 km.

Figure 10 
                  (a) The degree of partial melting shows 4−6% of garnet peridotite (4–8%) for the Sukadana Basalt, as shown in the Gd/Yb vs La/Yb diagram [52]. The curves Garnet 4% and Garnet 8% indicate the contents of garnet in the source, according to Halliday et al. [53]. (b) In the La/Yb vs Dy/Yb diagram, garnet peridotite is also a source region of the Sukadana Basalt. the diagram from Thirlwall et al. [54] and Bogaard and Worner [55]. This diagram indicates the low-degree melting (16–20%).
Figure 10

(a) The degree of partial melting shows 4−6% of garnet peridotite (4–8%) for the Sukadana Basalt, as shown in the Gd/Yb vs La/Yb diagram [52]. The curves Garnet 4% and Garnet 8% indicate the contents of garnet in the source, according to Halliday et al. [53]. (b) In the La/Yb vs Dy/Yb diagram, garnet peridotite is also a source region of the Sukadana Basalt. the diagram from Thirlwall et al. [54] and Bogaard and Worner [55]. This diagram indicates the low-degree melting (16–20%).

5.4 Mantle source

In the plot of Th/Yb vs Ta/Yb [56] for determining the tectonic setting of basaltic rocks, the Sukadana Basalt is classified as Within Plate Volcanic Zones and Within Plate Basalt (WPB) based on trace elements (Figure 11(a)). Samples with codes 2.1, 3.1, 3.2, 3.3, and 3.4 are located in the WPB zone. More specifically, samples with codes 3.1, 3.2, 3.3, and 3.4 are situated in the OIB region. The presence of samples with codes 2.1 and 3.4 in the WPB zone is consistent with the enrichment of Nb in the spider diagram explained earlier.

Figure 11 
                  (a) Ta/Yb vs Th/Yb diagram [56]. (b) Sukadana Basalt lavas in the Zr/Y vs Nb/Y diagram, relative to the mantle compositional components (grey square) and fields of rocks from various tectonic settings [83]. (c) Sukadana Basalt lavas in the Nb/Th vs Zr/Nb diagram after Weaver [39]. (d) Ti/100 vs Zr vs 3Y diagram [84]. (e) 2 Nb vs Zr/4 vs Y diagram [85]; (f) Hf/3 vs Th vs Ta diagram [86]. ARC: Arc-related basalts; DEP: deep depleted mantle; DM: depleted mantle; EM1 and EM2: enriched mantle sources; EN: enriched component; HIMU: high U/Pb mantle source; REC: recycled component; OPB: Oceanic Plateau Basalt, AFC: assimilation-fractional crystallization curve; EMB: enriched MORB; NMB: normal MORB; OIB: Oceanic Island Basalt; PM: Primitive Mantle from Sun and McDonough [34]; UC: Upper Crust from Taylor and McLennan [37]; ACM: Active Continental Margin; WPVZ: Within Plate Volcanic Zones; WPB: Within Plate Basalt, SC: Subduction Components.
Figure 11

(a) Ta/Yb vs Th/Yb diagram [56]. (b) Sukadana Basalt lavas in the Zr/Y vs Nb/Y diagram, relative to the mantle compositional components (grey square) and fields of rocks from various tectonic settings [83]. (c) Sukadana Basalt lavas in the Nb/Th vs Zr/Nb diagram after Weaver [39]. (d) Ti/100 vs Zr vs 3Y diagram [84]. (e) 2 Nb vs Zr/4 vs Y diagram [85]; (f) Hf/3 vs Th vs Ta diagram [86]. ARC: Arc-related basalts; DEP: deep depleted mantle; DM: depleted mantle; EM1 and EM2: enriched mantle sources; EN: enriched component; HIMU: high U/Pb mantle source; REC: recycled component; OPB: Oceanic Plateau Basalt, AFC: assimilation-fractional crystallization curve; EMB: enriched MORB; NMB: normal MORB; OIB: Oceanic Island Basalt; PM: Primitive Mantle from Sun and McDonough [34]; UC: Upper Crust from Taylor and McLennan [37]; ACM: Active Continental Margin; WPVZ: Within Plate Volcanic Zones; WPB: Within Plate Basalt, SC: Subduction Components.

Plotting of the samples in the Zr/Y vs Nb/Y diagram and Nb/Th vs Zr/Nb diagram (Figure 11(b) and (c)) suggests that they are derived from both Oceanic Plateau Basalt (OPB) and OIB-like mantle with slight influence from arc material. This result is also compatible with the other tectonic discrimination diagrams (Figure 11(d)–(f)). The Within Plate Volcanic Zones and ARC indicate the presence of slab-derived components in the Sukadana Basalt. In this study, the components are not significant. In addition, Oceanic Plateau Basalt shows diagnostic geochemical characteristics of E-MORB (i.e., rocks with Nb/La higher than 0.8; Nb-Ta depletion) to Transitional-Mid Oceanic Ridge Basalt [57] indicate extensional tectonic setting and shallow asthenosphere magma source [58,59]. The Sukadana Basalt fulfills all these characteristics. Based on this description, it can be inferred that there has been magma mixing between magma derived from the lower mantle (enriched) and the upper mantle (depleted) before extensional deformation occurred.

The isotope data from Gasparon [10], including 143Nd/144Nd vs 87Sr/85Sr, 87Sr/86Sr vs 206Pb/204Pb, 207Pb/204Pb vs 206Pb/204Pb, and 208Pb/204Pb vs 206Pb/204Pb, confirm that The Sukadana Basalt has characteristics of an OIB. Additionally, it is interesting to note that the magma source shows similarities with the oceanic crust of the South China Sea. This indicates the presence of magma flow from the South China Sea region, which later became part of the Sukadana Basalt magma (Figure 12(a)–(d)).

Figure 12 
                  (a) 143Nd/144Nd vs 87Sr/86Sr. (b) 87Sr/86Sr vs 206Pb/204Pb. (c) 207Pb/204Pb vs 206Pb/204Pb. (d) 208Pb/204Pb vs 206Pb/204Pb for the Sukadana Basalt from Lampung, Indonesia [10]. Data for the South China Sea are from the literature [87,88,89,90], OIB from Castillo [91] and for Indian Ocean-type MORB from Mahoney et al. [92], Sumatran Arc Granitoids and Sumatran Quaternary volcanics from Gasparon [10].
Figure 12

(a) 143Nd/144Nd vs 87Sr/86Sr. (b) 87Sr/86Sr vs 206Pb/204Pb. (c) 207Pb/204Pb vs 206Pb/204Pb. (d) 208Pb/204Pb vs 206Pb/204Pb for the Sukadana Basalt from Lampung, Indonesia [10]. Data for the South China Sea are from the literature [87,88,89,90], OIB from Castillo [91] and for Indian Ocean-type MORB from Mahoney et al. [92], Sumatran Arc Granitoids and Sumatran Quaternary volcanics from Gasparon [10].

Furthermore, the data suggest that the Sukadana Basalt magma has been contaminated by the continental crust of Sumatra (Sumatran Arc Granitoids). This is reasonable because, during its ascent to the Earth’s surface, it would have interacted with the continental crust of Sumatra, which is composed of granitic rocks. However, the whole rock geochemistry data indicate that this contamination is not significant.

The isotopic data also indicate that there is no influence from the magma forming the Indian-Australia oceanic crust and the magma from the quaternary volcanic arc on the composition of the Sukadana Basalt. The absence of the influence from the magma forming the Indian-Australia oceanic crust suggests the absence of slab window and slab tear around the Sukadana Basalt.

The Nb content as an indicator to determine the origin of the magma [44,45,60] shows that the Nb values range from 7.4 to 29.8 ppm, with an average of 15.9 ppm (Table 1). Two of the samples have particularly high Nb content with values exceeding 20 ppm (27.2 and 29.8 ppm). This anomaly has been observed in other studies as well [10,11]. In general, the Sukadana Basalt is considered to be a typical Nb-enriched basaltic rock, with relatively high Nb content (>6 ppm), (Nb/La)PM ratios >0.5 and Nb/U ratios >10, as well as enrichments in LREEs, LILEs, and HFSEs [61,62,63]. The high Nb values can be caused by two type sources [44,45] (1) OIB source mantle (similarly enriched with OIB-like isotopic and trace element signatures) and (2) partial melting of metasomatized enriched mantle wedge (similar to Normal Mid Oceanic Ridge Basalt like isotopes and OIB-like trace elements). The origin of an enriched mantle source (the first type) may result from either the flow of asthenospheric mantle through slab windows [64] or from the mixing of enriched (plume-type) and depleted components inside the mantle wedge [60]. The Sukadana Basalt exhibits high Nb/U ratios (18–60.3) with indication OIB mantle source [65,66] and is also supported by isotope data (Figure 12).

5.5 Geodynamic significance

Based on the Sukadana Basalt geochemical data, the mechanism of Sukadana Basalt magma emergence to the Earth’s surface can be described in two stages as follows:

  1. The first stage (Figure 13(a)), the subduction of the Indo-Australian plate, underwent roll-back phase. This subduction roll-back has been occurring since the Middle Eocene to the Quaternary period [67,68]. Subduction creates extensional deformation and thinning of the crust in the back-arc region [69,70,71,72]. Additionally, the opening of the Sunda Strait due to the movement of the Sumatra Fault causes extra extensional deformation in the back-arc region during the Quaternary period. At the same time, this subduction induces mantle return flow and upwelling magma from the asthenosphere (poloidal flow) in the lower back-arc region. This flow carries magma from the upper mantle to the lower mantle. Subsequently, magma originating from the lower mantle (OIB-type; enriched in LREE) rises to the upper mantle. This process creates an upper mantle dominated by OIB-like magma. Meanwhile, the other magma remains characteristic of OPB or E-MORB due to less influence from the OIB magma. In this stage, the magma underwent decompression melting in the garnet-peridotite source and then rises to the Earth’s surface through normal faults, reacting to Tuff from Quaternary Lampung Formation before appearing to the surface. This causes the magma to have a slight ARC-like composition, indicated by Pb enrichment. The Lampung Formation is older than Sukadana Basalt. Normal faults are the most favorable structures for magma to emerge from the mantle to the Earth’s surface [7379].

  2. The second stage (Figure 13(b)) is the ongoing subduction roll-back creating a wider accommodation space for OIB-type magma. During this stage, extensional deformation becomes more intense due to the opening of the Sunda Strait caused by the movement of the Sumatra Fault. Subduction roll-back results in the enrichment of LREE in magma, as indicated by increasing Nb content, while the intensive deformation reduces the contact time between the magma and continental crust indicated by decreasing Pb content. The process of decompression melting continues during this stage.

Figure 13 
                  The simple illustration of the rolling back subduction on the southern of Sumatra during Quaternary period.
Figure 13

The simple illustration of the rolling back subduction on the southern of Sumatra during Quaternary period.

If we look at Figure 6, K (potassium) enrichment occurs in all samples. The enrichment of K in basaltic magma can be caused by various possibilities [80]. However, based on the geochemical data of Sukadana Basalt, which shows characteristics of OIB, and there is a direct correlation between Nb and K, we argue that the K enrichment originates from the melting of paleo subducted crustal or sedimentary materials (Th/Nb = 0.1–0.3; Ba/Th < 200).

Besides, based on geochemical data as used in this study, morphological characteristics on the earth’s surface do not indicate a very wide uplift. Typically, a mantle plume results in a significant uplift within a region with a diameter ranging from 1,000 to 2,000 km [81], which was not observed in the Sukadana Basalt area. Furthermore, the appearance of Sukadana Basalt in the continental plate is not associated with the mantle plume mechanism.

6 Conclusions

Based on our comprehensive study of petrology and whole-rock geochemistry for the Sukadana Basalt in East Lampung, Sumatra, the following main conclusions can be drawn:

  1. The Sukadana Basalt originates from an OIB-type magma source and consists of two distinct petrogenesis groups, namely Group A and Group B. Group A shows enrichment of Hf, Pb, K, and Sr and depletion of Nb, and group B shows enrichment of K, Sr, and depletion of Pb. The trace elements have characteristics of sub-alkaline basalt affinity within a within-plate setting, minimal influence from subduction, or insignificant crustal contamination.

  2. The Sukadana Basalt magmas resulted from partial melting of the upwelling asthenospheric mantle, induced by the retreat of the paleo India–Australia oceanic plate. Decompression melting of the ascending asthenospheric mantle peridotite occurred at depths <80 km within the garnet stability zone.

  3. The weakening of Sundaland’s continental lithosphere during the Middle Eocene to Quaternary period was caused by the heat and melt flux associated with the uprising asthenosphere. This process facilitated extensive back-arc extension within an intracontinental setting.

  4. The sub-alkaline basaltic melts were transported to shallow crustal depths through two stages of extensional fault development implicated to the insignificant crustal contamination. The formation or re-activation of these faults was triggered by the movement of the Sumatra fault during the Quaternary.

Acknowledgements

This study was financially supported by the Sumatra Institute of Technology as part of the doctoral degree scholarship scheme. We are grateful to Heinz-Günter Stosch for the material and kind discussion, Sukadana society, ITB staff, Google, USGS, QGIS developer, Intertek Indonesia, Lampung Province government, and geological and geophysical students for their kind assistance during the field observation. We also would like to express our appreciation to the journal editor and reviewers for their valuable and insightful feedback. The constructive comments have significantly contributed to the improvement of our manuscript.

  1. Funding information: This study was financially supported by Sumatra Institute of Technology as part of the doctoral degree scholarship scheme.

  2. Author contributions: LPS contributed to the analysis of the data and the conception of the manuscript and revised it. BS, AR, AR, and IGBES contributed to the review of data gathering and supervising the analysis. All authors have read and agreed to the published version of the manuscript.

  3. Conflict of interest: The authors declare no conflict of interest

References

[1] Metcalfe I. Tectonic evolution of Sundaland. Bull Geol Soc Malaysia. 2017;63:27–60. https://gsm.org.my/content.php? id = 54&pid = 702001-101709.10.7186/bgsm63201702Search in Google Scholar

[2] Hall R, Morley CK. Sundaland basins. Geophys Monogr Ser. American Geophysical Union; 2004. p. 55–85.10.1029/149GM04Search in Google Scholar

[3] Pramumijoyo S, Sebrier M. Neogene and quaternary fault kinematics around the Sunda Strait area, Indonesia. J Southeast Asian Earth Sci. 1991;6:137–45.10.1016/0743-9547(91)90106-8Search in Google Scholar

[4] Stein S, Okal EA. Speed and size of the Sumatra earthquake. Nature. 2005;434:581–2. http://www.nature.com/articles/434581a.10.1038/434581aSearch in Google Scholar PubMed

[5] Abdurrachman M, Widiyantoro S, Priadi B, Ismail T. Geochemistry and structure of krakatoa volcano in the Sunda Strait, Indonesia. Geosci. 2018;8:1–10.10.3390/geosciences8040111Search in Google Scholar

[6] Ringwood A. Petrogenesis of intraplate magmas and structure of the upper mantle. Chem Geol. 1990;82:187–207.10.1016/0009-2541(90)90081-HSearch in Google Scholar

[7] Li J, Ding W, Lin J, Xu Y, Kong F, Li S, et al. Dynamic processes of the curved subduction system in Southeast Asia: A review and future perspective. Earth-Science Rev. 2021;217:103647. 10.1016/j.earscirev.2021.103647.Search in Google Scholar

[8] Dosso L, Joron J-L, Maury R, Bougalut H. Isotopic (Sr, Nd) and trace element study of back-arc basalts behind the Sunda arc. Terra Cogn. 1987;7:398.Search in Google Scholar

[9] Nishimura S, Nishida J, Yokoyama T, Hehuwat F. Neo-tectonics of the Strait of Sunda, Indonesia. J Southeast Asian Earth Sci. 1986;1:81–91.10.1016/0743-9547(86)90023-1Search in Google Scholar

[10] Gasparon M. Origin and evolution of mafic volcanics of Sumatra (Indonesia): their mantle sources, and the roles of subducted oceanic sediments and crustal contamination. Australia: Univerity of Tasmania; 1993.Search in Google Scholar

[11] Romeur M. Series magmatiques arc et arriere-arc de la sonde: Nature des sources impliquees (Elements En Trace Et Isotopes Sr-Nd-Pb). Brest, France: L’UNIVERSITE DE BRETAGNE OCCIDENTALE; 1991.Search in Google Scholar

[12] McCaffrey R, Abers G. Orogeny in arc–continent collision: The Banda Arc and western New Guinea. Geology. 1991;19:563–66.10.1130/0091-7613(1991)019<0563:OIACCT>2.3.CO;2Search in Google Scholar

[13] Hall R. Hydrocarbon basins in SE Asia: understanding why they are there. Pet Geosci. 2009;15(2):131–46.10.1144/1354-079309-830Search in Google Scholar

[14] Barber A, Crow M, Milson J. Sumatra: geology, resources and tectonic evolution. United Kingdom: The Geological Society London; 2005.Search in Google Scholar

[15] Huchon P, Le Pichon X. Sunda Strait and Central Sumatra fault. Geology. 1984;12:668–72.10.1130/0091-7613(1984)12<668:SSACSF>2.0.CO;2Search in Google Scholar

[16] Malod JA, Karta K, Beslier MO, Zen MT. From normal to oblique subduction: Tectonic relationships between Java and Sumatra. J Southeast Asian Earth Sci. 1995;12:85–93.10.1016/0743-9547(95)00023-2Search in Google Scholar

[17] Ninkovich D. Late Cenozoic clockwise rotation of Sumatra. Earth Planet Sci Lett. 1976;29:269–75.10.1016/0012-821X(76)90130-8Search in Google Scholar

[18] Gasparon M, Hilton DR, Varne R. EPSL Crustal contamination processes traced by helium isotopes: Examples from the Sunda arc, Indonesia. Earth Planet Sci Lett. 1994;126:15–22.10.1016/0012-821X(94)90239-9Search in Google Scholar

[19] Harjono H, Diament M, Dubois J, Larue M, Zen MT. Seismicity of the Sunda Strait: Evidence for crustal extension and volcanological implications. Tectonics. 1991;10:17–30.10.1029/90TC00285Search in Google Scholar

[20] Soeria-Atmaja R, Maury R, Bougault H, Joron J, Bellon H, Hasanunddin D. Présence de tholeiites d’arrière − arc Quatenariés en Indonésie: Les basaltes de Sukadana (Sud de Sumatra). Réunion des Sci la Terre. Clermont–Ferrand; 1986.Search in Google Scholar

[21] Hasibuan RF, Ohba T, Abdurrachman M, Hoshide T. Temporal variations of petrological characteristics of Tangkil and Rajabasa Volcanic Rocks, Indonesia. Indones J Geosci. 2020;7:135–59. http://ijog.geologi.esdm.go.id/index.php/IJOG/article/view/651.10.17014/ijog.7.2.135-159Search in Google Scholar

[22] Mangga S, Amirudin, Suwarti T, Gafoer S, Sidarto. Peta Geologi Lembar TanjungKarang, Sumatera. Bandung, Indonesia: Geological Survey Center; 1993.Search in Google Scholar

[23] Nukman M, Moeck I. Structural controls on a geothermal system in the Tarutung Basin, north central Sumatra. J Asian Earth Sci. 2013;74:86–96. 10.1016/j.jseaes.2013.06.012.Search in Google Scholar

[24] Winter JD. Principles of Igneous and Metamorphic Petrology. 2nd edn. New York: Cambridge University Press; 2009.Search in Google Scholar

[25] Rosenstengel LM, Hartmann L. Geochemical stratigraphy of lavas and fault-block structures in the Ametista do Sul geode mining district, Paraná volcanic province, southern Brazil. Ore Geol Rev. 2012;5:332–348.10.1016/j.oregeorev.2012.05.003Search in Google Scholar

[26] Hartmann LA, Medeiros JTN, Baggio SB, Antunes LM. Controls on prolate and oblate geode geometries in the Veia Alta basalt flow, largest world producer of amethyst, Paraná volcanic province, Brazil. Ore Geol Rev. 2015;66:243–51. https://linkinghub.elsevier.com/retrieve/pii/S0169136814002868.10.1016/j.oregeorev.2014.11.005Search in Google Scholar

[27] Le Bas MJ, Maitre RWL, Streckeisen A, Zanettin B. A chemical classification of volcanic rocks based on the total alkali-silica diagram. J Petrol. 1986;27:745–50.10.1093/petrology/27.3.745Search in Google Scholar

[28] Zulkarnain I. Geochemical evidence of Island-Arc Origin for Sumatra Island; A New Perspective based on Volcanic Rocks in Lampung Province, Indonesia. Indones. J Geosci. 2011;6:213–25.10.17014/ijog.v6i4.128Search in Google Scholar

[29] Irvine TN, Baragar WR. A guide to the classification of the common volcanic rocks. Can J Earth Sci. 1971;8:235–458.10.1139/e71-055Search in Google Scholar

[30] Le Maitre R. Igneous rocks: A classification and glossary of terms. Recomentations of the International Union of Geological Sciences Subcommission on the Systematics of Igneous Rocks. Cambridge: Cambridge University Press; 2002.Search in Google Scholar

[31] McDonough WF, Sun S-s. The composition of the Earth. Chem Geol. 1995;120:223–53. https://linkinghub.elsevier.com/retrieve/pii/0009254194001404.10.1016/0009-2541(94)00140-4Search in Google Scholar

[32] Ashwal LD. Sub-lithospheric mantle sources for overlapping southern African Large Igneous Provinces. South African. J Geol. 2021;124:421–42. https://pubs.geoscienceworld.org/sajg/article/124/2/421/598793/Sub-lithospheric-mantle-sources-for-overlapping.10.25131/sajg.124.0023Search in Google Scholar

[33] O’Hara MJ. Are Ocean Floor Basalts Primary Magma. Nature. 1968;220:683–6. https://www.nature.com/articles/220683a0.10.1038/220683a0Search in Google Scholar

[34] Sun SS, McDonough WF. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. Geol Soc Spec Publ. 1989;42:313–45.10.1144/GSL.SP.1989.042.01.19Search in Google Scholar

[35] Kieffer B, Arndt N, Lapierre H, Bastien F, Bosch D, Pecher A, et al. Flood and shield basalts from Ethiopia: magmas from the African Superswell. J Petrol. 2004;45:793–834.10.1093/petrology/egg112Search in Google Scholar

[36] Zhang C-L, Yang D-S, Wang H-Y, Takahashi Y, Ye H-M. Neoproterozoic mafic-ultramafic layered intrusion in Quruqtagh of northeastern Tarim Block, NW China: Two phases of mafic igneous activity with different mantle sources. Gondwana Res. 2011;19:177–90. https://linkinghub.elsevier.com/retrieve/pii/S1342937X10000729.10.1016/j.gr.2010.03.012Search in Google Scholar

[37] Taylor S, McLennan S. The Continental Crust: Its Composition and Evolution. Carlton: Blackwell Scientific Publication; 1985.Search in Google Scholar

[38] Saunders AD, Norry MJ, Tarney J. Origin of MORB and chemically-depleted mantle reservoirs: Trace element constraints. In: Menzies M, Cox K, editors. In: Ocean Cont Lithosph Similarities Differ. J. Petrol. Oxford University Press; 1988. p. 415–445.10.1093/petrology/Special_Volume.1.415Search in Google Scholar

[39] Weaver BL. The origin of ocean island basalt end-member compositions: trace element and isotopic constraints. Earth Planet Sci Lett. 1991;104:381–97.10.1016/0012-821X(91)90217-6Search in Google Scholar

[40] Barry T. Petrogenesis of cenozoic basalts from Mongolia: Evidence for the role of asthenospheric versus metasomatized lithospheric mantle sources. J Petrol. 2003;44:55–91. 10.1093/petrology/44.1.55.Search in Google Scholar

[41] Rooney T, Furman T, Bastow I, Ayalew D, Yirgu G. Lithospheric modification during crustal extension in the Main Ethiopian Rift. J Geophys Res. 2007;112:B10201. 10.1029/2006JB004916.Search in Google Scholar

[42] Sheldrick TC, Barry TL, Van Hinsbergen DJJ, Kempton PD. Constraining lithospheric removal and asthenospheric input to melts in Central Asia: A geochemical study of Triassic to Cretaceous magmatic rocks in the Gobi Altai (Mongolia). Lithos. 2018;296–299:297–315. https://linkinghub.elsevier.com/retrieve/pii/S0024493717303997.10.1016/j.lithos.2017.11.016Search in Google Scholar

[43] Bora DK, Borah K, Goyal A. Crustal shear-wave velocity structure beneath Sumatra from receiver function modeling. J Asian Earth Sci. 2016;121:127–38. 10.1016/j.jseaes.2016.03.007 Search in Google Scholar

[44] Moradi S, Khaksar T, Nazarinia A, Hussain A. Petrology and geochemistry of Plio-Quaternary high-Nb basalts from Shahr-e-Babak area:Insights into post-collision magmatic processes in the Kerman Cenozoic Magmatic Arc. Geol Acta. 2022;20:1–19. https://revistes.ub.edu/index.php/GEOACTA/article/view/36839.10.1344/GeologicaActa2022.20.8Search in Google Scholar

[45] Chen F, Cui X, Lin S, Wang J, Ren G, Li K, et al. The earliest Neoproterozoic Nb-enriched mafic magmatism indicates subduction tectonics in the southwestern Yangtze Block, South China. Precambrian Res. 2023;384:106938. 10.1016/j.precamres.2022.106938.Search in Google Scholar

[46] Zhang Z, Li S, Wang G, Li X, Wang G, Suo Y, et al. Plume interaction and mantle heterogeneity: A geochemical perspective. Geosci Front. 2020;11:1571–9. https://linkinghub.elsevier.com/retrieve/pii/S1674987120300517.10.1016/j.gsf.2020.02.009Search in Google Scholar

[47] Misra KC. Introduction to Geochemistry Principles and Applications. 1st edn. United Kingdom: Wiley-Blackwell; 2012.Search in Google Scholar

[48] Manikyamba C, Ganguly S, Santosh M, Saha A, Lakshminarayana G. Geochemistry and petrogenesis of Rajahmundry trap basalts of Krishna-Godavari Basin. India Geosci Front. 2015;6:437–51. 10.1016/j.gsf.2014.05.003.Search in Google Scholar

[49] Frey F, Green D, Roy S. Integrated models of basalt petrogenesis: A study of quartz tholeiites to olivine melilitites from South Eastern Australia utilizing geochemical and experimental petrological data. Petrology. 1978;19:463–513.10.1093/petrology/19.3.463Search in Google Scholar

[50] Wilson M. Igneous Petrogenesis: a Global Tectonic Approach. Mineral Mag 1989;53:88.10.1007/978-1-4020-6788-4Search in Google Scholar

[51] Sharma M. Siberian traps. In: Mahoney JJ, Coffin MF, editors. Large Igneous Prov Cont Ocean Planet Flood Volcanism. New York: American Geophysical Union Geophysical Monograph; 1997. p. 273–95.10.1029/GM100p0273Search in Google Scholar

[52] Yokoyama T, Aka FT, Kusakabe M, Nakamura E. Plume–lithosphere interaction beneath MtConstraints from 238U–230Th–226Ra and Sr–Nd–Pb isotope systematics. Geochim Cosmochim Acta. 71, Cameroon volcano, West Africa; 2007. p. 1835–54. https://linkinghub.elsevier.com/retrieve/pii/S0016703707000282.10.1016/j.gca.2007.01.010Search in Google Scholar

[53] Halliday AN, Lee DC, Tommasini S, Davies GR, Paslick CR, Fitton JG, et al. Incompatible trace elements in OIB and MORB and source enrichment in the sub oceanic mantle. Earth Planet Sci Lett. 1995;133:379–395.10.1016/0012-821X(95)00097-VSearch in Google Scholar

[54] Thirlwall F, Upton B, Jenkins C. Interaction between con tinental lithosphere and Iceland plume Sr-Nd-Pb isotope geo − chemistry of tertiary basalts. NE Greenland Petrology. 1994;35:839–79.10.1093/petrology/35.3.839Search in Google Scholar

[55] Bogaard PJ, Worner G. Petrogenesis of Basanitic to Tholeiitic Volcanic Rocks from the Miocene Vogelsberg, Central Germany. J Petrol. 2003;44:569–602. 10.1093/petrology/44.3.569.Search in Google Scholar

[56] Pearce JA. Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust. Lithos. 2008;100:14–48. https://linkinghub.elsevier.com/retrieve/pii/S0024493707001375.10.1016/j.lithos.2007.06.016Search in Google Scholar

[57] Xia L, Li X. Basalt geochemistry as a diagnostic indicator of tectonic setting. Gondwana Res. 2019;65:43–67. 10.1016/j.gr.2018.08.006.Search in Google Scholar

[58] Hamilton WB. The closed upper-mantle circulation of plate tectonics. 2013;30:359–410. 10.1029/GD030p0359.Search in Google Scholar

[59] Anderson DL, Zhang Y-S, Tanimoto T. Plume heads, continental lithosphere, flood basalts and tomography. Vol. 68. Geol Soc London: Spec Publ; 1992. p. 99–124. 10.1144/GSL.SP.1992.068.01.07.Search in Google Scholar

[60] Macpherson CG, Chiang KK, Hall R, Nowell GM, Castillo PR, Thirlwall MF. Plio-Pleistocene intra-plate magmatism from the southern Sulu Arc, Semporna peninsula, Sabah, Borneo: Implications for high-Nb basalt in subduction zones. J Volcanol Geotherm Res. 2010;190:25–38. 10.1016/j.jvolgeores.2009.11.004.Search in Google Scholar

[61] Hastie AR, Mitchell SF, Kerr AC, Minifie MJ, Millar IL. Geochemistry of rare high-Nb basalt lavas: Are they derived from a mantle wedge metasomatised by slab melts? Geochim Cosmochim Acta. 2011;75:5049–72. https://linkinghub.elsevier.com/retrieve/pii/S0016703711003462.10.1016/j.gca.2011.06.018Search in Google Scholar

[62] Liao F-X, Chen N-S, Santosh M, Wang Q-Y, Gong S-L, He C, et al. Paleoproterozoic Nb–enriched meta-gabbros in the Quanji Massif. NW China: Implications for assembly of the Columbia supercontinent Geosci Front. 2018;9:577–90. https://linkinghub.elsevier.com/retrieve/pii/S1674987117301032.10.1016/j.gsf.2017.05.007Search in Google Scholar

[63] Han Q, Peng S. Paleoproterozoic subduction within the Yangtze Craton: Constraints from Nb-enriched mafic dikes in the Kongling complex. Precambrian Res. 2020;340:105634. https://linkinghub.elsevier.com/retrieve/pii/S0301926819300749.10.1016/j.precamres.2020.105634Search in Google Scholar

[64] Castillo PR. Origin of the adakite-high-Nb basalt association and its implications for postsubduction magmatism in Baja California, Mexico. Geol Soc Am Bull. 2008;120:451–62. https://pubs.geoscienceworld.org/gsabulletin/article/120/3-4/451-462/2266.10.1130/B26166.1Search in Google Scholar

[65] Kepezhinskas P, Defant MJ, Drummond MS. Progressive enrichment of island arc mantle by melt-peridotite interaction inferred from Kamchatka xenoliths. Geochim Cosmochim Acta. 1996;60:1217–29. https://linkinghub.elsevier.com/retrieve/pii/0016703796000014.10.1016/0016-7037(96)00001-4Search in Google Scholar

[66] Hofmann AW. Sampling Mantle Heterogeneity through Oceanic Basalts: Isotopes and Trace Elements. Treatise on Geochemistry. USA: Elsevier; 2007. p. 1–44. https://linkinghub.elsevier.com/retrieve/pii/B008043751602123X.10.1016/B0-08-043751-6/02123-XSearch in Google Scholar

[67] Pubellier M, Morley CK. The basins of Sundaland (SE Asia): Evolution and boundary conditions. Mar Pet Geol. 2014;58:555–78. 10.1016/j.marpetgeo.2013.11.019.Search in Google Scholar

[68] Lai Y-M, Chung S-L, Ghani AA, Murtadha S, Lee H-Y, Chu M-F. Mid-Miocene volcanic migration in the westernmost Sunda arc induced by India-Eurasia collision. Geology. 2021;49:713–7. https://pubs.geoscienceworld.org/gsa/geology/article/49/6/713/595184/Mid-Miocene-volcanic-migration-in-the-westernmost.10.1130/G48568.1Search in Google Scholar

[69] Chen Z, Schellart WP, Strak V, Duarte JC. Does subduction-induced mantle flow drive backarc extension. Earth Planet Sci Lett. 2016;441:200–10. 10.1016/j.epsl.2016.02.027.Search in Google Scholar

[70] Schellart WP, Moresi L. A new driving mechanism for backarc extension and backarc shortening through slab sinking induced toroidal and poloidal mantle flow: Results from dynamic subduction models with an overriding plate. J Geophys Res Solid Earth. 2013;118:3221–48.10.1002/jgrb.50173Search in Google Scholar

[71] Xue K, Schellart WP, Strak V. Overriding plate deformation and topography during slab rollback and slab rollover: insights from subduction experiments. Tectonics. 2022;41:1–19.10.1029/2021TC007089Search in Google Scholar

[72] Feng G, DIlek Y, Niu X, Liu F, Yang J. Geochemistry and geochronology of OIB-type, Early Jurassic magmatism in the Zhangguangcai range, NE China, as a result of continental back-arc extension. Geol Mag. 2021;158:143–57.10.1017/S0016756818000705Search in Google Scholar

[73] Ayalew D, Jung S, Romer RL, Garbe-Schönberg D. Trace element systematics and Nd, Sr and Pb isotopes of Pliocene flood basalt magmas (Ethiopian rift): A case for Afar plume-lithosphere interaction. Chem Geol. 2018;493:172–88.10.1016/j.chemgeo.2018.05.037Search in Google Scholar

[74] Faccenna C, Becker TW, Lallemand S, Lagabrielle Y, Funiciello F, Piromallo C. Subduction-triggered magmatic pulses: A new class of plumes. Earth Planet Sci Lett. 2010;299:54–68. 10.1016/j.epsl.2010.08.012.Search in Google Scholar

[75] Zi JW, Haines PW, Wang XC, Jourdan F, Rasmussen B, Halverson GP, et al. Pyroxene 40Ar/39Ar dating of basalt and applications to large igneous provinces and precambrian stratigraphic correlations. J Geophys Res Solid Earth. 2019;124:8313–30.10.1029/2019JB017713Search in Google Scholar

[76] Shahraki M. Dynamics of mantle circulation and convection: The signatures in the satellite derived gravity fields. Germany: Johann Wolfgang Goethe University; 2013.Search in Google Scholar

[77] Yan Q, Shi X, Metcalfe I, Liu S, Xu T, Kornkanitnan N, et al. Hainan mantle plume produced late Cenozoic basaltic rocks in Thailand, Southeast Asia. Sci. Rep. Nature.com; 2018. https://www.nature.com/articles/s41598-018-20712-7.10.1038/s41598-018-20712-7Search in Google Scholar PubMed PubMed Central

[78] Wang Y, Santosh M, Luo Z, Hao J. Large igneous provinces linked to supercontinent assembly. J Geodyn. 2015;85:1–10. 10.1016/j.jog.2014.12.001.Search in Google Scholar

[79] De Souza ZS, Vasconcelos PM, Knesel KM, da Silveira Dias LG, Roesner EH, Cordeiro de Farias PR, et al. The tectonic evolution of Cenozoic extensional basins, northeast Brazil: Geochronological constraints from continental basalt 40Ar/39Ar ages. J South Am Earth Sci. 2013;48:159–72. 10.1016/j.jsames.2013.09.008.Search in Google Scholar

[80] Gupta AK. Origin of potassium-rich silica-deficient igneous rocks. India: Springer; 201510.1007/978-81-322-2083-1Search in Google Scholar

[81] Ritter JRR, Christensen UR. Mantle plumes: A Multidiscip Approach. Heidelberg, Germany: Springer; 2007.10.1007/978-3-540-68046-8Search in Google Scholar

[82] Hans Wedepohl K. The composition of the continental crust. Geochim Cosmochim Acta. 1995;59:1217–32. https://linkinghub.elsevier.com/retrieve/pii/0016703795000382.10.1016/0016-7037(95)00038-2Search in Google Scholar

[83] Condie KC. High field strength element ratios in Archean basalts: A window to evolving sources of mantle plumes? Lithos. 2005;79:491–504.10.1016/j.lithos.2004.09.014Search in Google Scholar

[84] Pearce JA, Cann JR. Tectonic setting of basic volcanic rocks determined using trace element analyses. Earth Planet Sci Lett. 1973;19:290–300. https://linkinghub.elsevier.com/retrieve/pii/0012821X73901295.10.1016/0012-821X(73)90129-5Search in Google Scholar

[85] Meschede M. A method of discriminating between different types of mid-ocean ridge basalts and continental tholeiites with the Nb□1bZr□1bY diagram. Chem Geol. 1986;56:207–18. https://linkinghub.elsevier.com/retrieve/pii/0009254186900045.10.1016/0009-2541(86)90004-5Search in Google Scholar

[86] Wood CA. Morphometric analysis of cinder cone degradation. J Volcanol Geotherm Res. 1980;8:137–60.10.1016/0377-0273(80)90101-8Search in Google Scholar

[87] Yan Q, Shi X, Wang K, Bu W, Xiao L. Major element, trace element, and Sr, Nd and Pb isotope studies of Cenozoic basalts from the South China Sea. Sci China Ser D Earth Sci. Vol. 51. 2008. p. 550–66. 10.1007/s11430-008-0026-3.Search in Google Scholar

[88] Tu K, Flower MFJ, Carlson RW, Xie G, Chen C-Y, Zhang M. Magmatism in the South China Basin. Chem Geol. 1992;97:47–63. https://linkinghub.elsevier.com/retrieve/pii/000925419290135R.10.1016/0009-2541(92)90135-RSearch in Google Scholar

[89] Yan Q, Shi X, Castillo PR. The late Mesozoic–Cenozoic tectonic evolution of the South China Sea: A petrologic perspective. J Asian Earth Sci 2014(85):178–201. https://linkinghub.elsevier.com/retrieve/pii/S1367912014000510.10.1016/j.jseaes.2014.02.005Search in Google Scholar

[90] Yan Q, Castillo P, Shi X, Wang L, Liao L, Ren J. Geochemistry and petrogenesis of volcanic rocks from Daimao Seamount (South China Sea) and their tectonic implications. Lithos. 2015;218–219:117–26. https://linkinghub.elsevier.com/retrieve/pii/S0024493715000092.10.1016/j.lithos.2014.12.023Search in Google Scholar

[91] Castillo P. The Dupal anomaly as a trace of the upwelling lower mantle. Nature. 1988;336:667–70. https://www.nature.com/articles/336667a0.10.1038/336667a0Search in Google Scholar

[92] Mahoney JJ, Natland JH, White WM, Poreda R, Bloomer SH, Fisher RL, et al. Isotopic and geochemical provinces of the western Indian Ocean Spreading Centers. J Geophys Res Solid Earth. 1989;94:4033–52. 10.1029/JB094iB04p04033.Search in Google Scholar

Received: 2023-05-22
Revised: 2023-09-01
Accepted: 2023-09-02
Published Online: 2023-11-15

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Diagenesis and evolution of deep tight reservoirs: A case study of the fourth member of Shahejie Formation (cg: 50.4-42 Ma) in Bozhong Sag
  3. Petrography and mineralogy of the Oligocene flysch in Ionian Zone, Albania: Implications for the evolution of sediment provenance and paleoenvironment
  4. Biostratigraphy of the Late Campanian–Maastrichtian of the Duwi Basin, Red Sea, Egypt
  5. Structural deformation and its implication for hydrocarbon accumulation in the Wuxia fault belt, northwestern Junggar basin, China
  6. Carbonate texture identification using multi-layer perceptron neural network
  7. Metallogenic model of the Hongqiling Cu–Ni sulfide intrusions, Central Asian Orogenic Belt: Insight from long-period magnetotellurics
  8. Assessments of recent Global Geopotential Models based on GPS/levelling and gravity data along coastal zones of Egypt
  9. Accuracy assessment and improvement of SRTM, ASTER, FABDEM, and MERIT DEMs by polynomial and optimization algorithm: A case study (Khuzestan Province, Iran)
  10. Uncertainty assessment of 3D geological models based on spatial diffusion and merging model
  11. Evaluation of dynamic behavior of varved clays from the Warsaw ice-dammed lake, Poland
  12. Impact of AMSU-A and MHS radiances assimilation on Typhoon Megi (2016) forecasting
  13. Contribution to the building of a weather information service for solar panel cleaning operations at Diass plant (Senegal, Western Sahel)
  14. Measuring spatiotemporal accessibility to healthcare with multimodal transport modes in the dynamic traffic environment
  15. Mathematical model for conversion of groundwater flow from confined to unconfined aquifers with power law processes
  16. NSP variation on SWAT with high-resolution data: A case study
  17. Reconstruction of paleoglacial equilibrium-line altitudes during the Last Glacial Maximum in the Diancang Massif, Northwest Yunnan Province, China
  18. A prediction model for Xiangyang Neolithic sites based on a random forest algorithm
  19. Determining the long-term impact area of coastal thermal discharge based on a harmonic model of sea surface temperature
  20. Origin of block accumulations based on the near-surface geophysics
  21. Investigating the limestone quarries as geoheritage sites: Case of Mardin ancient quarry
  22. Population genetics and pedigree geography of Trionychia japonica in the four mountains of Henan Province and the Taihang Mountains
  23. Performance audit evaluation of marine development projects based on SPA and BP neural network model
  24. Study on the Early Cretaceous fluvial-desert sedimentary paleogeography in the Northwest of Ordos Basin
  25. Detecting window line using an improved stacked hourglass network based on new real-world building façade dataset
  26. Automated identification and mapping of geological folds in cross sections
  27. Silicate and carbonate mixed shelf formation and its controlling factors, a case study from the Cambrian Canglangpu formation in Sichuan basin, China
  28. Ground penetrating radar and magnetic gradient distribution approach for subsurface investigation of solution pipes in post-glacial settings
  29. Research on pore structures of fine-grained carbonate reservoirs and their influence on waterflood development
  30. Risk assessment of rain-induced debris flow in the lower reaches of Yajiang River based on GIS and CF coupling models
  31. Multifractal analysis of temporal and spatial characteristics of earthquakes in Eurasian seismic belt
  32. Surface deformation and damage of 2022 (M 6.8) Luding earthquake in China and its tectonic implications
  33. Differential analysis of landscape patterns of land cover products in tropical marine climate zones – A case study in Malaysia
  34. DEM-based analysis of tectonic geomorphologic characteristics and tectonic activity intensity of the Dabanghe River Basin in South China Karst
  35. Distribution, pollution levels, and health risk assessment of heavy metals in groundwater in the main pepper production area of China
  36. Study on soil quality effect of reconstructing by Pisha sandstone and sand soil
  37. Understanding the characteristics of loess strata and quaternary climate changes in Luochuan, Shaanxi Province, China, through core analysis
  38. Dynamic variation of groundwater level and its influencing factors in typical oasis irrigated areas in Northwest China
  39. Creating digital maps for geotechnical characteristics of soil based on GIS technology and remote sensing
  40. Changes in the course of constant loading consolidation in soil with modeled granulometric composition contaminated with petroleum substances
  41. Correlation between the deformation of mineral crystal structures and fault activity: A case study of the Yingxiu-Beichuan fault and the Milin fault
  42. Cognitive characteristics of the Qiang religious culture and its influencing factors in Southwest China
  43. Spatiotemporal variation characteristics analysis of infrastructure iron stock in China based on nighttime light data
  44. Interpretation of aeromagnetic and remote sensing data of Auchi and Idah sheets of the Benin-arm Anambra basin: Implication of mineral resources
  45. Building element recognition with MTL-AINet considering view perspectives
  46. Characteristics of the present crustal deformation in the Tibetan Plateau and its relationship with strong earthquakes
  47. Influence of fractures in tight sandstone oil reservoir on hydrocarbon accumulation: A case study of Yanchang Formation in southeastern Ordos Basin
  48. Nutrient assessment and land reclamation in the Loess hills and Gulch region in the context of gully control
  49. Handling imbalanced data in supervised machine learning for lithological mapping using remote sensing and airborne geophysical data
  50. Spatial variation of soil nutrients and evaluation of cultivated land quality based on field scale
  51. Lignin analysis of sediments from around 2,000 to 1,000 years ago (Jiulong River estuary, southeast China)
  52. Assessing OpenStreetMap roads fitness-for-use for disaster risk assessment in developing countries: The case of Burundi
  53. Transforming text into knowledge graph: Extracting and structuring information from spatial development plans
  54. A symmetrical exponential model of soil temperature in temperate steppe regions of China
  55. A landslide susceptibility assessment method based on auto-encoder improved deep belief network
  56. Numerical simulation analysis of ecological monitoring of small reservoir dam based on maximum entropy algorithm
  57. Morphometry of the cold-climate Bory Stobrawskie Dune Field (SW Poland): Evidence for multi-phase Lateglacial aeolian activity within the European Sand Belt
  58. Adopting a new approach for finding missing people using GIS techniques: A case study in Saudi Arabia’s desert area
  59. Geological earthquake simulations generated by kinematic heterogeneous energy-based method: Self-arrested ruptures and asperity criterion
  60. Semi-automated classification of layered rock slopes using digital elevation model and geological map
  61. Geochemical characteristics of arc fractionated I-type granitoids of eastern Tak Batholith, Thailand
  62. Lithology classification of igneous rocks using C-band and L-band dual-polarization SAR data
  63. Analysis of artificial intelligence approaches to predict the wall deflection induced by deep excavation
  64. Evaluation of the current in situ stress in the middle Permian Maokou Formation in the Longnüsi area of the central Sichuan Basin, China
  65. Utilizing microresistivity image logs to recognize conglomeratic channel architectural elements of Baikouquan Formation in slope of Mahu Sag
  66. Resistivity cutoff of low-resistivity and low-contrast pays in sandstone reservoirs from conventional well logs: A case of Paleogene Enping Formation in A-Oilfield, Pearl River Mouth Basin, South China Sea
  67. Examining the evacuation routes of the sister village program by using the ant colony optimization algorithm
  68. Spatial objects classification using machine learning and spatial walk algorithm
  69. Study on the stabilization mechanism of aeolian sandy soil formation by adding a natural soft rock
  70. Bump feature detection of the road surface based on the Bi-LSTM
  71. The origin and evolution of the ore-forming fluids at the Manondo-Choma gold prospect, Kirk range, southern Malawi
  72. A retrieval model of surface geochemistry composition based on remotely sensed data
  73. Exploring the spatial dynamics of cultural facilities based on multi-source data: A case study of Nanjing’s art institutions
  74. Study of pore-throat structure characteristics and fluid mobility of Chang 7 tight sandstone reservoir in Jiyuan area, Ordos Basin
  75. Study of fracturing fluid re-discharge based on percolation experiments and sampling tests – An example of Fuling shale gas Jiangdong block, China
  76. Impacts of marine cloud brightening scheme on climatic extremes in the Tibetan Plateau
  77. Ecological protection on the West Coast of Taiwan Strait under economic zone construction: A case study of land use in Yueqing
  78. The time-dependent deformation and damage constitutive model of rock based on dynamic disturbance tests
  79. Evaluation of spatial form of rural ecological landscape and vulnerability of water ecological environment based on analytic hierarchy process
  80. Fingerprint of magma mixture in the leucogranites: Spectroscopic and petrochemical approach, Kalebalta-Central Anatolia, Türkiye
  81. Principles of self-calibration and visual effects for digital camera distortion
  82. UAV-based doline mapping in Brazilian karst: A cave heritage protection reconnaissance
  83. Evaluation and low carbon ecological urban–rural planning and construction based on energy planning mechanism
  84. Modified non-local means: A novel denoising approach to process gravity field data
  85. A novel travel route planning method based on an ant colony optimization algorithm
  86. Effect of time-variant NDVI on landside susceptibility: A case study in Quang Ngai province, Vietnam
  87. Regional tectonic uplift indicated by geomorphological parameters in the Bahe River Basin, central China
  88. Computer information technology-based green excavation of tunnels in complex strata and technical decision of deformation control
  89. Spatial evolution of coastal environmental enterprises: An exploration of driving factors in Jiangsu Province
  90. A comparative assessment and geospatial simulation of three hydrological models in urban basins
  91. Aquaculture industry under the blue transformation in Jiangsu, China: Structure evolution and spatial agglomeration
  92. Quantitative and qualitative interpretation of community partitions by map overlaying and calculating the distribution of related geographical features
  93. Numerical investigation of gravity-grouted soil-nail pullout capacity in sand
  94. Analysis of heavy pollution weather in Shenyang City and numerical simulation of main pollutants
  95. Road cut slope stability analysis for static and dynamic (pseudo-static analysis) loading conditions
  96. Forest biomass assessment combining field inventorying and remote sensing data
  97. Late Jurassic Haobugao granites from the southern Great Xing’an Range, NE China: Implications for postcollision extension of the Mongol–Okhotsk Ocean
  98. Petrogenesis of the Sukadana Basalt based on petrology and whole rock geochemistry, Lampung, Indonesia: Geodynamic significances
  99. Numerical study on the group wall effect of nodular diaphragm wall foundation in high-rise buildings
  100. Water resources utilization and tourism environment assessment based on water footprint
  101. Geochemical evaluation of the carbonaceous shale associated with the Permian Mikambeni Formation of the Tuli Basin for potential gas generation, South Africa
  102. Detection and characterization of lineaments using gravity data in the south-west Cameroon zone: Hydrogeological implications
  103. Study on spatial pattern of tourism landscape resources in county cities of Yangtze River Economic Belt
  104. The effect of weathering on drillability of dolomites
  105. Noise masking of near-surface scattering (heterogeneities) on subsurface seismic reflectivity
  106. Query optimization-oriented lateral expansion method of distributed geological borehole database
  107. Petrogenesis of the Morobe Granodiorite and their shoshonitic mafic microgranular enclaves in Maramuni arc, Papua New Guinea
  108. Environmental health risk assessment of urban water sources based on fuzzy set theory
  109. Spatial distribution of urban basic education resources in Shanghai: Accessibility and supply-demand matching evaluation
  110. Spatiotemporal changes in land use and residential satisfaction in the Huai River-Gaoyou Lake Rim area
  111. Walkaway vertical seismic profiling first-arrival traveltime tomography with velocity structure constraints
  112. Study on the evaluation system and risk factor traceability of receiving water body
  113. Predicting copper-polymetallic deposits in Kalatag using the weight of evidence model and novel data sources
  114. Temporal dynamics of green urban areas in Romania. A comparison between spatial and statistical data
  115. Passenger flow forecast of tourist attraction based on MACBL in LBS big data environment
  116. Varying particle size selectivity of soil erosion along a cultivated catena
  117. Relationship between annual soil erosion and surface runoff in Wadi Hanifa sub-basins
  118. Influence of nappe structure on the Carboniferous volcanic reservoir in the middle of the Hongche Fault Zone, Junggar Basin, China
  119. Dynamic analysis of MSE wall subjected to surface vibration loading
  120. Pre-collisional architecture of the European distal margin: Inferences from the high-pressure continental units of central Corsica (France)
  121. The interrelation of natural diversity with tourism in Kosovo
  122. Assessment of geosites as a basis for geotourism development: A case study of the Toplica District, Serbia
  123. IG-YOLOv5-based underwater biological recognition and detection for marine protection
  124. Monitoring drought dynamics using remote sensing-based combined drought index in Ergene Basin, Türkiye
  125. Review Articles
  126. The actual state of the geodetic and cartographic resources and legislation in Poland
  127. Evaluation studies of the new mining projects
  128. Comparison and significance of grain size parameters of the Menyuan loess calculated using different methods
  129. Scientometric analysis of flood forecasting for Asia region and discussion on machine learning methods
  130. Rainfall-induced transportation embankment failure: A review
  131. Rapid Communication
  132. Branch fault discovered in Tangshan fault zone on the Kaiping-Guye boundary, North China
  133. Technical Note
  134. Introducing an intelligent multi-level retrieval method for mineral resource potential evaluation result data
  135. Erratum
  136. Erratum to “Forest cover assessment using remote-sensing techniques in Crete Island, Greece”
  137. Addendum
  138. The relationship between heat flow and seismicity in global tectonically active zones
  139. Commentary
  140. Improved entropy weight methods and their comparisons in evaluating the high-quality development of Qinghai, China
  141. Special Issue: Geoethics 2022 - Part II
  142. Loess and geotourism potential of the Braničevo District (NE Serbia): From overexploitation to paleoclimate interpretation
Downloaded on 8.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/geo-2022-0544/html
Scroll to top button