Magmatic-tectonic response of the South China Craton to the Paleo-Pacific subduction during the Triassic: a new viewpoint based on Well NK-1*

Chan WANG, Hao LIU, Gang LI, Zhen SUN, Tingting GONG, Li MIAO,Xiaowei ZHU, Yunying ZHANG, Weihai XU,**, Wen YAN,6,**

1 Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Chinese Academy of Sciences,Guangzhou 510301, China

2 Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China

3 Key Laboratory of Tectonic Controlled Mineralization and Oil Reservoir of Ministry of Natural Resources, Chengdu University of Technology, Chengdu 610059, China

4 School of Earth Sciences and Resources, China University of Geosciences (Beijing), Beijing 100083, China

5 College of Earth Science, Chengdu University of Technology, Chengdu 610059, China

6 University of Chinese Academy of Sciences, Beijing 100049, China

Abstract The Nansha Block (NB) is one of the blocks separated from the southern margin of the South China Craton (SCC) by the western Pacific subduction, which contains rich information of geodynamic and tectonic transformation.To reveal the essence of western Paleo-Pacific subduction during the Triassic period,Well NK-1 in this block was selected for petrographic study, and published research data from other cooperative teams were compared.A double-cycle pattern of basic to neutral magmatic volcanism was established, and 36 lithological rhythmic layers and representative cryptoexplosive breccia facies and welded tuff bands were identified.Combined with a reanalysis of published geochronological data, geochemical elements, and isotope geochemistry, we found that the rock assemblages could be divided into an intermediate-acid dacite (DA) series (SiO2＞65%) and basaltic (BA) series (Co＜40 μg/g), which was formed during the early Late Triassic ((218.6±3.2)-(217.9±3.5) Ma).BA exhibits obvious calc-alkaline island-arc magmatic properties: (87Sr/86Sr)i ratio ranging 0.703 77-0.711 18 (average: 0.706 45), 147Sm/144Nd ratio ranging 0.119-0.193 (average: 0.168), and chondrite-normalized rare earth element (REE) curves being flat,while DA exhibits remarkable characteristics of subducted island-arc andesitic magma: (87Sr/86Sr)i ratio(0.709 39-0.711 29; average: 0.710 35),εNd(t) value (-6.2--4.8; average: -5.6) and εHf(t) value (-2.9--1.7,average: -2.2) show obvious crust-mantle mixing characteristics.BA and DA reveal typical characteristics of island-arc magma systems and type Ⅱ enriched mantle (EM-Ⅱ) magma.BA magma was likely resulted from the process whereby the continental crust frontal accretionary wedge was driven by the Paleo-Pacific slab subduction into the deep and began to melt, resulting in a large amount of melt (fluid) joined the asthenosphere on the side of the continental margin.In contrast, DA magma was likely resulted from the process whereby the plate front was forced to bend with increasing subduction distance, which triggered the upwelling of the asthenosphere near the continent and subsequently led to the partial melting of the lithospheric mantle and lower crust due to continuous underplating.The lithospheric thinning environment in the study area at the end of Triassic created suitable conditions for the separation between the NB and SCC, which provided an opportunity for the formation of the early intracontinental rift during the later expansion of the South China Sea (SCS).

Keyword: Paleo-Pacific tectonic domain; tectonic evolution; rift volcanic rock; Nansha Block; South China Craton

1 INTRODUCTION

The Pangea formed on a global scale with the merging of Gondwana between Laurasia during the Permian Period and was surrounded by an ocean known as the Panthalassa or the Paleo-Pacific (Gao et al., 2017; Yan et al., 2017; Hara et al., 2018; Xu et al., 2018; Luo et al., 2021).Later, the Paleo-Pacific gradually expanded and subducted along the eastern margin of Eurasia, resulting in large-scale intracontinental deformation and vast magmatic provinces (Xu et al., 2017; Liu et al., 2019a, b; Sun et al., 2020; Zhu et al., 2020); a cordillera-type linear mountain system extending thousands of kilometers formed on the western margin of the North and South American continents (Zhao et al.,2007, 2016; Li et al., 2018a).These two different subduction responses and their deep genetic mechanisms have aroused great interest from geologists (Zhang et al., 2021).The South China Craton (SCC), as a part of Eurasia, has experienced a series of important geological events (such as the Indosinian movement, Yanshanian movement,destruction of the SCC, and formation of large-scale metal mineral resources in the eastern SCC) under the influence of subduction of the Paleo-Pacific Plate and its superimposed effects (Hong et al.,2017; Li and Zhou, 2018; Hou et al., 2019; Sun et al., 2019; Liu et al., 2021a; Shu et al., 2021; Xie et al., 2021).The orogenic mode and dynamic mechanism of the southeastern margin of the SCC during the Permo-Triassic Indosinian movement have always been controversial.Some scholars think that the Hercynian-Indosinian movement was caused by the closure of the Jinshajiang-Ailaoshan-Song Ma branch ocean basin and the subsequent collision between the Indochina Plate and the SCC(Holloway, 1981; Hall, 2012; Ding et al., 2013; Jian et al., 2019; Chen et al., 2020).Other scholars assert that no collisional orogeny in South China occurred during the Permo-Triassic period but that a remote intraplate collisional orogeny was caused by the closure of the Paleo-Tethys Ocean between the Indochina Plate and the Sibumasu Plate (Jahn, 2010;Shao et al., 2015; Li et al., 2018b; Duan et al.,2020).Other scholars have proposed that tectonic movement was caused by the subduction of the Paleo-Pacific Plate along the southeastern margin of the SCC toward the northwest since the Permian(Wu and Suppe, 2018; Zhang et al., 2018).

The South China Sea (SCS) is a unique marginal sea located at the junction among the Eurasian Plate,Indo-Australian Plate, and Pacific Plate (Yi et al.,2021).Therefore, its formation and evolutionary history have always been a hotspot in the study of ocean-continental transition zones.The Nansha Block (NB), one of the many blocks in the SCS, was a product of the SCC detachment influenced by the Paleo-Pacific subduction mechanism (Luo et al.,2021).Geophysical data suggest that the NB has a Proterozoic folded continental basement that displays affinities with the SCC, which has been confirmed by Pb isotope data for the rocks (Zartman et al.,1988; Yan et al., 2011; Huang et al., 2013; Miao et al.,2021; Wei et al., 2022).Moreover, the NB separated from the Indochina-South China paleocontinent and drifted southeast and southward since the Late Cretaceous; these events are considered to have been caused by the far-field effect of the collision between the Indian Ocean Plate and the Eurasian Plate.Seafloor spreading was halted when the Australian Plate front subducted under the Sunda continent during the middle Miocene (Li et al.,2017; Li and Zhou, 2018).However, no direct petrological evidence has been reported for the above hypothesis, greatly hindering the in-depth study of the tectonic evolution of the NB.For example,when did the NB split from the SCC? What was the dissociation process? What role does the NB played in the formation and evolution of the SCS? Many issues require further exploration and analysis.

Triassic arc magmatic rocks, as important products of Paleo-Pacific subduction, play an important role in the reconstruction of the Paleo-Pacific subduction model on the southeastern margin of the SCC (Xie et al., 2021; Yang et al., 2021).In this study, we selected core samples from Well NK-1, which is a very important well drilled through volcanic rocks in the NB, southern SCS, as the research subject and extracted important geological information regarding the apparent nearly 1 100-m thick Triassic volcanic rocks.This paper presents a comprehensive study based on the detailed analysis of petrological,petrographic, and mineralogical characteristics of volcanic rocks in Well NK-1 and the sample test data obtained in previous work.The aims are to establish a model of the magmatic eruption cycle,and characterized the magma source area, the genetic mechanism of volcanic rocks and the tectonic setting of magmatic dynamics, and finally to reveal the tectonic evolution process of the NB.This study promoted the understanding of the evolution of the Paleo-Pacific tectonic domain on the eastern margin of Eurasia, which effectively constrains the formation age of volcanic rocks in Well NK-1 and exploration of the affinity between the NB and SCC.

2 GEOLOGICAL SETTING

The SCS, one of the largest marginal seas in the western Pacific Ocean, is located at the junction of the three largest plates in the world at present: the Pacific Plate on the east, the Eurasian Plate on the north-northwest, and the Indo-Australian Plate on the south-southwest (Xu et al., 2012; Yi et al.,2021).The multiple extensional basins formed on the eastern margin of the SCS are closely related to the expansion of the SCS that began at ～34 Ma and stopped at ～16 Ma (Chen et al., 1992; Liu et al.,1997; Huang et al., 2013).Following the cessation of seafloor spreading, the intraplate oceanic island basalt associated with the Hainan mantle plume covered large areas of the SCS and its adjacent areas(e.g., Yan et al., 2011, 2017; Shao et al., 2015; Sun et al., 2020).The NB, located in the southwestern SCS, covers an area of approximately 60×104km2.Its northwestern boundary is a residual spreading ridge in the middle of the southwestern basin and a fault extending to the southwest.Its southern boundary is the Lupal fault.Its eastern boundary is the Zhongnan-Lile fault (Yan et al., 2011; Miao et al., 2021; Yi et al., 2021).There were four tectonic stages in the NB in the pre-Middle Mesozoic geology: (1) the Archean to the Paleoproterozoic(4 000-1 600 Ma) crystalline basement-formation,(2) the Mesoproterozoic to the Neoproterozoic(1 600-541 Ma) rifting, (3) the early Paleozoic(541-359 Ma) intracontinental orogeny, and (4) the early Mesozoic (251-200 Ma) Indosinian orogeny(Wang and Li, 2001, 2003; Zou et al., 2020, 2021;Liu et al., 2021b; Wang et al., 2021).Subsequently,the NB separated from the SCC and complex continental-ocean transition processes occurred from the middle of the Mesozoic to the Cenozoic Era.

3 PETROLOGY

This study is the first attempt to present highresolution detailed observations and rock type divisions of the 975.3-2 020.2-m core column in Well NK-1,as well as the microscopic identification of 85 representative rock thin sections.

Optical thin section observations were conducted,and orthogonally polarized photomicrographs were obtained with a microscope (Carl Zeiss Axio Scope A1, Germany).All preparation of optical thin sections was completed at Guangzhou Tuoyan Testing Technology Co., Ltd., and laboratory identification was completed in the optical section identification room of the National Geology Experimental Teaching Demonstration Center of Chengdu University of Technology.

The bulk geochemistry, zircon U-Pb chronology,and Sr-Nd-Pd-Hf isotope data in this study were obtained from Miao et al.(2021) and Wei et al.(2022) (Supplementary Tables S1-S5).

3.1 Core sample characteristics

Well NK-1 is located on the Meiji Reef, an oval coral atoll in the central Nansha Islands 240 km away from southern Palawan Island (Fig.1).The total drilling depth was 2 020 m, and the core recovery percentage was 91%.A volcanic rock section with a thickness of 1 021.3 m was drilled(Fig.2), and its color changed clearly from dark to light from bottom to top (some of the core segments are interbedded) at macro level.The volcanic rock segment underlies gray-black quartz sandstone and overlies gray-white mixed deposits of clastic/carbonate/reef sediments.According to the division model of volcanic eruption cycles and subcycles, the first cycle of mafic→neutral magmatic volcanism is evident from 2 017.5 m to 1 290.2 m, and the second cycle of mafic→intermediate-felsic magmatic volcanism appears from 1 290.2 m to 996.2 m.The first cycle can be further divided into a mafic magma eruption subcycle (2 017.5-1 963.2 m)→mafic magma overflow subcycle (1 963.2-1 469.2 m)→neutral magma overflow subcycle (1 469.2-1 290.6 m).The second cycle is further divided into a mafic magma overflow subcycle (1 290.6-1 131.5 m)→neutral-mafic magma overflow subcycle (1 131.5-1 077.8 m)→intermediatefelsic magma overflow subcycle (1 077.8-996.2 m).Pelitic tuff at 1 263.7-1 267.4 m, lenticular welded tuff bands at 1 275.2 m, and tabular volcanic gravel at the top of volcanic rocks in the second cycle are also present, suggesting that the magma was a distinct underwater eruption (Wu et al., 2019).This combination of multiple rhythms and numerous cycles is the comprehensive result of multiple magma eruptions and evolution (Gong et al., 2021b).

Fig.1 Survey of the research area and drilling location(after Miao et al., 2021)

3.2 Microscopic characteristics

The petrographic characteristics are newly reported in this paper.Based on the core catalog and detailed identification of 53 rock slices, the first and second cycle characteristics of lithologic facies are clear, at the boundary of approximately 1 400 m; the first cycle consists of basaltic andesite and a combination of trachyte-trachyandesite (Fig.2).The second cycle contains mainly a set of alterations of pyroxene-bearing basalt, basaltic trachyandesitetrachybasalt, and trachyandesite-trachydacite with local tuff lava assemblages.The lavas all have porphyritic structures, and welded tuff or tuff lava textures can be seen in tuff lavas.

The phenocrysts are mainly composed of feldspar and dark minerals, with a small amount of quartz in the second cycle.

Andesine: This mineral is found in all types of lavas and tuff lavas; grains are mostly idiomorphichypidiomorphic and tabular and often develop alteration and twin crystals.In the late stage, weak to moderate saussuritization, epidotization, potassic alteration, and clayization occurred.The direction of plagioclase twin-crystal binding in the upper cycle is inclined in the (010) direction; in addition,alkaline feldspar coronas can be seen at the margins(Fig.3a), or grains are embedded with each other and have a distributed glomeroporphyritic texture(Fig.3b).Moreover, some grains in trachyte and trachybasalt display twin-crystal deformation and sand bell-type structures.

Dark minerals: These minerals are common in trachybasalt and trachyandesite in the two cycles,which have no complete and fresh grains; however,the shapes of early euhedral columns (similar to amphibole or alkaline pyroxene, such as aegirineaugite and aegirine) are retained, and the particle sizes are within 1 mm.Many inclusions, such as apatite, are common, and some particles occasionally retain dark edges.The alteration characteristics show that the stages are in the order of hightemperature actinolitization→lepidolitization (chlorite of the high-temperature phase)+epidotization+silicification and lower-temperature chlorite and carbonation.

Alkali feldspar: These grains are distributed in trachyandesite and trachyte and are mainly perthite;orthoclase retains the sanidine in shape.The perthite has developed a variety of perthitic textures and types; for example, evidence of sanidine can be seen in Carlsbad twins and cleavage in the direction of oblique crossing (001) (Fig.3c).Inlaid metasomatic plagioclase can be observed in the first cycle,although metasomatism is weak overall, and zoisitization is relatively strong.In the second cycle,a few alkaline feldspars are independent minerals,while the remainder are metasomatic plagioclases with “trim shapes” and “patch shapes”, and the main alteration is kaolinization.

Quartz: This mineral is rounded and occasionally found in the upper trachydacite, where it retains an atypical high-temperatureβ-quartz morphology(Fig.3d).

Matrix: A palimpsest vitreous texture (retained arc compression textures) (Fig.3e) can be seen in the trachyandesite, trachyte, and trachydacite matrices,with conversion to trachyte-cryptocrystalline texture after recrystallization, but spherules that are not completely devitrified are commonly retained.The composition consists of feldspar, a few dark minerals,and cryptocrystalline and dusty iron aggregates.The relatively coarse matrix in the basalt-basaltic andesite is composed of feldspar microcrystals(mostly sericite+clayization+epidotization) (Fig.3f) and dark minerals (pyroxene/amphibole in the early stage plus actinolite, biotite, chlorite, muscovite,etc.) with a mostly intergranular-intersertal texture and local pilotaxitic texture.

Fig.2 Lithological column sketch diagram of volcanic rock in Well NK-1

In addition, amygdaloid bodies filled with quartz+epidote and occasionally with calcite are sporadically found in the lower BA.

In some tuffaceous-bearing lavas, a few plastic magma fragments, crystals, and rigid homologous detritus can be found, and tension-shear fractures have developed in the rocks, which are filled with lepidolite+rectangular/deformed quartz, contemporaneous with the alteration of dark minerals.

Among the accessory minerals, the volcanic rocks from Well NK-1 in the SCS generally form a titanomagnetite-apatite-zircon assemblage rich in zircon, apatite, and other minerals with incompatible elements.

Apatite: This mineral is distributed mainly in the two cycles in the form of long columns, acicular columns, hexagonal columns, and hexagonal bipyramids.Some apatite grains contain fine apatite grains or other minerals, particles of metallic minerals and gas and liquid inclusions.

Fig.3 Microscopic characteristics of Well NK-1 volcanic rocks

Titanomagnetite: This mineral is mostly distributed in the lower part, with contents of 0.5-2 vol.% and relatively few fresh particles; the grains are mostly decomposed and partially or completely converted to an aggregation of leucosphenite+sphene, with three groups of residual cleavage remaining.

Zircon: Trace zircon grains were found in the two cycles; they display positive high relief and tertiary interference colors and have magmatic zircon shapes.A small fraction of the grain face is incomplete and may have burst.

Late tensile-shear fractures have developed in the rocks, and most of these fractures are filled with deformed rectangular quartz+lepidolite.

According to the mineral characteristics of the rocks, titanomagnetite and zircon crystallized earlier among the accessory minerals, while apatite may have had a longer crystallization time span.

In summary, the lithographic characteristics show that the rocks in the early magmatic cycle (in the first cycles) crystallized with a large amount of volatile apatite and that the phenocryst mineral assemblage is dominated by plagioclase+pyroxene/hornblende.The magma was calc-alkaline to moderately alkaline and relatively mafic with high volatility.However, the content of K-feldspar among the phenocrysts of the second magmatic cycle obviously increases, clearly indicating that the magma source region underwent a transformation(Wang et al., 2018).

4 DISCUSSION

4.1 Classification of rock types

Based on the first detailed core characteristics analysis of Well NK-1, combined with the geochemical data published by other teams (Miao et al., 2021; Wei et al., 2022), the development of volcanic eruptions producing rhythmic layers and cycle assemblages indicate the continuation of magmatic eruption (Valentine et al., 2014).However,the alternate lithology evolution during this process suggests that the tectonic environment may be complex (Li, 1993).The microscopic characteristics show obvious two-cycle differences, and the characteristics of the first cycle imply that the mineral assemblage of the early phenocrysts was plagioclase+pyroxene+hornblende, with apatite that crystallized early.In general, the characteristics are reflected by calc-alkaline and mafic rocks, while the magma is rich in volatile components.A small number of plagioclase grains underwent orthoclasation/sanidinization in the later period, with polysynthetic twins crossing obliquely in the (010) direction; these features may indicate that the magma changed in the direction of slightly alkaline composition, potassization,and rapid dehydration, which are characteristic of multisource magma or fluid containing molten lava(Le et al., 1986).The observation of the grains in the second cycle shows autogenous sanidine and quartz phenocrysts, and the potassization is obviously weaker than that in the first cycle, which reflects the evolutionary trend of a homogeneous magma enriched in Si and K.Geochemical data were collected from 79 samples within the 1 489-1 225.4-m rhythmic interval spanning the first and second cycles and 1 485-1 686.8-m rhythmic interval covering the first cycle (Supplementary Tables S1-S3).The sample data in the study area fall into the basic basalt and intermediate-acid trachyte/trachyte dacite areas in Fig.4a, and combined with the characteristics of calc-alkaline basalt and high-K and shoshonite series dacite areas in Fig.4b, the volcanic rock cores of Well NK-1 are divided into the intermediate-acid DA and the basic BA in this paper.

4.2 Magma evolution process

4.2.1 BA evolution process

BA from Well NK-1 showed less variable SiO2content (47.09-53.37 wt.%, average 49.94 wt.%)and fluctuating MgO content (4.53-11.47 wt.%,average 6.91 wt.%).The negative correlation between MgO and SiO2and Al2O3, positive correlation between CaO/Al2O3, Ni and Cr, and positive correlation between TFe2O3, TiO2, and P2O5(Fig.5)show that magma in the BA source area had weak fractional crystallization of olivine, pyroxene, titanic oxide, and apatite, which is consistent with the titanomagnetite and apatite observed in the maceral analysis in Section 3.2.The relatively continuous data distribution also indicates that the BA could be originated from a unified parent magma chamber(Liu et al., 2019c).The low TiO2contents(0.60 wt.%-2.10 wt.%; average: 1.15) and light rare earth element (LREE)/heavy rare earth element(HREE)＞1 (0.70-2.84; average: 1.20) of the BA also suggest that the magma in the source region may have undergone a low degree of crystallization differentiation (Pearce and Norry, 1979).

Fig.4 SiO2 (wt.%) vs.Nb/Y diagram (after Winchester and Floyd (1977)) (a) and Th (μg/g) vs.Co (μg/g) diagram (after Stevenson (2018)) (b)

The linear relationship of the BA data is in overall consistent with the trend of partial melting(Fig.6), which indicates that the magma chamber was dominated by partial melting with a degree of fractional crystallization (Pati?o Douce and Beard,1995).This is consistent with the CIPW calculation results for the BA (the differentiation index (DI)value ranges from 21.66 to 52.88, with an average value of 33.90, and the solidification index (SI)value approaches that of the initial magma, i.e.,23.44-45.99, with an average of 31.44).These characteristics also confirm that BA geochemical data have good indicative significance for the recovery of magma source areas (Rapp and Watson,1995).BA data largely fall into the stable melting zone of spinel at ＜80-km depth (Fig.7a), which shows that the magma chamber in the source area was located in the transition region between the upper mantle asthenosphere and subcontinental lithospheric mantle (SCLM).In the partial melting trend diagram (Fig.7a), BA samples plot between the enriched mid-ocean ridge basalt (E-MORB) and primitive mantle (PM) fields, showing an enrichment trend, which is highly consistent with curve (5) in Fig.7b.However, the data points near curve (3) in Fig.7b indicate that the magmatic source region may have also experienced partial melting of spinel+garnet lherzolite (the partial melting degree is 10%-20%).Because the contents of Zr (average: 111.13 μg/g) and Hf (average:2.67 μg/g) in BA are lower than those in the (upper continental crust) UCC (the average values are 190 μg/g and 5.8 μg/g, respectively) (Supplementary Table S2), it is concluded that the magma in the source area was not obviously affected by crustal contamination during the ascent process (Gao et al.,2017).

Fig.5 Bivariate plots of selected characteristic components vs.MgO (wt.%) of BA from Well NK-1

Fig.6 Binary diagrams illustrating the Y (μg/g)-Ce (μg/g) (a), La/Sm-La (μg/g) (b), Ba/Zr-Ba (μg/g) (c), and Th/Nb-Th (μg/g)(d) compositions of BA and DA rocks from Well NK-1

Fig.7 (Tb/Yb)N vs.(La/Sm)N diagram (a) (after Wang et al.(2019)) and Sm/Yb vs.La/Sm diagram (b) (after Li et al.(2018a))

4.2.2 DA evolution process

Compared to those of the UCC, DA exhibited a higher LREE/HREE ratio (average: 2.81), greater enrichment of large-ion lithophile elements (LILEs,such as Rb and K) and high-field strength elements(HFSEs) (such as Th, U, Zr, and Hf) and greater depletion of LILEs (such as Sr) (Fig.8).Combined with the higher DI value (average: 79.2) and lower SI value (average: 6.6) calculated by DA CIPW, the andesitic magma chamber in the study area experienced obvious crystallization differentiation(Gardien et al., 1995), which is consistent with the characteristics of the dominant fractional crystallization in the magma chamber of the DA source area, as shown in Fig.6.In the fractional crystallization and assimilation discrimination diagram (Fig.9), DA is dominated by fractional crystallization, but the degree of fractional crystallization is weak (Fig.9d),which indicates that the obvious element geochemical fractionation characteristics of DA should reflect fractional crystallization in the source region.The geochemical characteristics of DA are significantly different from those of BA, and the trace elements and REEs are significantly different, which suggests that the two types of rocks may derive from different magma sources (Iwamori and Nakamura, 2015).

DA samples show good continuity in Harker diagram (Fig.10), which should be the product of magma evolution in the same source area (Wang et al., 2019).The contents of MgO, CaO, Al2O3,TFe2O3, TiO2, P2O5, Sr, and Cr in the DA samples decreased significantly with the increase in SiO2content, at the same time, combined with the obvious depletion of Ba, Eu, Nb, Ta, and Ti in DA.These results indicate that the magma in the DA source area experienced fractional crystallization of mafic minerals (such as pyroxene and amphibole),oxides containing Fe-Ti (such as titanomagnetite),plagioclase, K-feldspar, and apatite (Le and Thompson, 1988).In the Sr-Rb vs.Ba-Sr diagram(Fig.11), DA shows characteristics consistent with the trend line of biotite (Bi), suggesting that the source region may be water-rich and mafic-rich, which is also consistent with the enrichment of biotite minerals reflected in macerals in Section 3.2.

Fig.8 Chondrite-normalized REE distribution pattern of the DA (a); primitive mantle-normalized spider diagram of DA(b), normalized after Sun and McDonough (1989)

Fig.9 The fractional crystallization (FC) and assimilation (A) on the evolution of the magma producing the subvolcanic rocks (after Nicolae and Saccani (2003))

The origin of the parent magma of DA is worth discussing.Based on the discrimination diagram of multiple source areas, most DA samples fall into the overlap area of the partial melting of metasandstone and partial melting of the basic rock field (Fig.12a)and of the A-type granite formed by partial melting of intermediate-basic igneous rocks and biotite gneissic rock-derived melt field (Fig.12b-d).This feature suggests that the magma in the DA source area of the study area has obvious subduction arc properties, which is most likely due to partial melting of continental lithospheric mantle and lower crust (Defant and Drummond, 1990; Li et al.,2018a).Figure 12e-f shows that the DA samples fall nearby sandstone source rock and shale, showing high agreement with sand-rich protolith and clay-poor source regions, which is similar to the characteristics of the Bethanga pluton in the Lachlan Fold belt, Australia.These characteristics suggest that the DA source area has obvious crust-mantle mixing characteristics rather than the partial melting of a single crust-derived material (Stevenson, 2018).Moreover,the geochemical element enrichment characteristics of DA are more obvious than those of UCC.

The basic and intermediate-acidic rock assemblages in Well NK-1 show geochemical characteristics similar to those of island-arc magmas.BA obviously belongs to the island-arc calc-alkaline basalt, while DA shows the characteristics of island-arc andesitic magmas, which originated from different magma chambers (Xu et al., 2017; Wei et al., 2022).The magma in the BA source area is dominated by the partial melting of spinel and garnet lherzolite with weak fractional crystallization (fractional crystallization of mafic minerals, minerals containing Fe-Ti and apatite).Geochemistry shows the trend of the chondrite-normalized REE curve similar to that of PM (Fig.13a) and the PM-normalized curve similar to that of (lower continental crust) LCC (Fig.13b).The geochemical characteristics of DA show that the source area is influenced by crust-mantle interactions, and the magma evolution shows limited crystallization differentiation and weak assimilation and contamination.Fractional crystallization is mainly dominated by the biotite phase with the addition of the fractional crystallization of mafic minerals, minerals containing Fe-Ti, apatite, plagioclase,and K-feldspar.The rock properties of BA and DA suggest that they originated during the mature island-arc stage influenced by the subduction of the western Pacific Ocean and were constrained by complex magmatism under the influence of subduction-related tectonics (Gao et al., 2017;Hernández-Uribe et al., 2020).The interbedded production of BA and DA rocks revealed by petrographic characteristics indicates that the late magma chamber eruption mode was dominated by upwelling accompanied by a double magma system,forming a rhythmic magma eruption environment.

Fig.10 Bivariate plots of selected characteristic components vs.SiO2 (wt.%) of DA from Well NK-1

4.3 Genetic type

4.3.1 Genetic analysis of BA

The BA samples retrieved from Well NK-1 exhibit obvious subduction island-arc calc-alkaline basalt characteristics (Fig.4b), with a flat chondritenormalized REE pattern (Fig.13a).The BA is more depleted of HFSEs such as Th, U, Ta, and Nb, more enriched by LILEs such as Rb and K than LCC based on the PM-normalized spider diagram(Fig.13b), and has various isotope ratios of (87Sr/86Sr)i(0.703 77-0.711 18, 0.706 45 on average) and147Sm/144Nd (0.119-0.193, 0.168 on average)(Supplementary Table S5).These characteristics indicate that the magma of the BA source area is relatively enriched in incompatible elements, and the material composition of the source area is relatively complex, with the addition of crustderived characteristic material (Blichert-Toft and Albarède, 1999; Iwamori and Nakamura, 2015).

Fig.11 Simulated calculation showing the process of fractional crystallization for Sr (μg/g) vs.Rb (μg/g) (a) and Ba (μg/g) vs.Sr (μg/g) (b) (after Pitcher et al., 1985)

Several sets of Ta and Nb values are close to or lower than the PM (Fig.13b) without obvious crystallization differentiation of the BA source area that analyzed in Section 4.2.1, which indicates that the initial source area is likely a PM or depleted mantle (DM) (Huang et al., 2013; Li et al., 2018a).The main reason for the enrichment in incompatible elements in BA may be related to the addition of oceanic crust subduction fluid and melt rather than the melting of a typical crust-derived accretionary wedge(Gee et al., 2020).TheδCe (δCe=(Cesample/Cechondrite)/values (0.96-1.04, 1.01 on average) of the BA are lower than those of the UCC (1.11 on average), which indicates that a certain amount of marine sediments may have been trapped during the subduction of the oceanic crust (Singh and Johannes, 1996; Li et al., 2012).Previous studies have shown that these BA magma source characteristics can be caused by arc magmas related to subduction, lithospheric mantle inherited from melt or fluid metasomatism, or contamination from felsic crustal materials (Defant and Drummond,1990; Hoskin and Schaltegger, 2003).The contamination of felsic crustal materials was first excluded by the nonsignificant contamination from the BA in Section 4.2.1; however, the obvious enrichment of incompatible elements in BA also excludes the possibility of arc magmas related to subduction.The low (87Sr/86Sr)iratio (0.703 77-0.711 18, 0.706 45 on average) and the large fluctuationεNd(t) values (-5.7-5.5, 2.7 on average) of the BA (Supplementary Table S5) indicate that the dual influence of mantle- and ancient crust-derived material in the magma chamber was probably caused by the transformation of lithospheric mantle by subduction melt (fluid) (Jahn, 2010; Hong et al.,2017).In Fig.14a, BA samples fall within a triangular pattern between the 10% UCC and 50%LCC margins, indicating that the mixed substances of UCC and LCC are added to the source area.Figure 14b shows that BA samples fall into the upper right of the bulk silicate earth (BSE) point.Both series of rocks show an evolutionary trend toward type EM-Ⅱ enriched mantle (EM-Ⅱ), which indicates that the BA source area may not have been affected by a single siliceous and aluminous crustal material.The greater depletion of some elements in BA may not only be the enriched geochemical characteristics of metasomatic lithospheric mantle formation but may also be caused by the magmatic source area originating from PM or DM (Fig.13)(Pati?o Douce and Johnston, 1991).

Fig.12 Elemental discrimination diagrams

Fig.13 The chondrite-normalized REE distribution pattern of the BA (a) and the primitive mantle-normalized spider diagram of the BA (b) (normalized after Sun and McDonough (1989))

Fig.14 Discrimination diagrams of εNd(t) vs.(87Sr/86Sr)i (after Zindler and Hart (1986)) (a) and 143Nd/144Nd vs.(87Sr/86Sr)i(after Zindler and Hart (1986)) (b)

A small number of BA samples were scattered in the mantle plume source region and the normal midocean ridge basalt (N-MORB) and arc magma(ARC) source region boundary, concentrated in the interval between PM and upper crust (UC)(Fig.15a), and showed magma contamination characteristics (Fig.15b).This indicates that the crustal and mantle contamination characteristics of the basic magma chamber in the study area are obvious (Rapp and Watson, 1995).Combined with the fact that the BA (except NK-1-V2-88 and NK-1-V2-79) does not exhibit positive Zr and Hf anomalies, it is suggested that the assimilation and contamination during magma upwelling should be weak, and the element enrichment and crust-mantle contamination characteristics are derived from the magma chamber properties in the source region.The BA has the characteristics of typical subducted island-arc volcanic rocks (Fig.15c).In Fig.15d, BA data fall into the oceanic island-arc range affected by the volcanic arc and the calc-alkaline basalt region transitive to the active continental margin.Combined with the trend consistent with the evolution of the intraoceanic arc (Fig.15e), it is suggested that the BA should be a mature island-arc volcanic rock, and its forming environment is biased to the side of the continental margin of the intraoceanic arc, with obvious continental crust pollution attributes (Fig.15f).It is worth noting that although NK-1-V2-193, NK-1-V2-79, and NK-1-V2-88 samples in BA obviously have high Nb content (up to 29.66 μg/g), their Zr and Hf enrichment trends are not obvious, which is different from previous studies of Nb-rich basalts.The analysis (Fig.16a-b) shows that the overall data fall into the normal arc magma field, showing the characteristics of classical island-arc magmatic rocks, which indicates that the enrichment of incompatible elements in BA is more likely caused by the fluid or melt formed by the continental crust accretionary wedge accompanied by oceanic crust subduction (Zhao et al., 2016; Yu et al., 2020).

Fig.15 Discrimination diagrams of Nb/Y vs.Zr/Y (after Fitton et al.(1997)) (a); Th/Yb vs.Nb/Yb (after Pearce (2008))(b); La/Nb vs.La (μg/g) (after Wang et al.(2019)) (c); Th/Yb vs.Ta/Yb (after Pearce (1982)) (d); Ce (μg/g) vs.Yb(μg/g) (after Stevenson (2018)) (e); and La/Ba vs.La/Nb (after Miller et al.(1999)) (f)

Based on the above characteristics, it is concluded that the BA features obvious properties of subduction zone composition (SZC), and the EM-Ⅱcharacteristics shown by Sr-Nd isotope data indicate that the magma source area may be influenced by crust-mantle contamination (e.g., subducted recirculating sediments, recirculating oceanic crust, and SCLM)(Huang et al., 2013; Gee et al., 2020).According to the petrological characteristics in Section 3.2, the BA has a relatively euhedral to subhedral mafic hydrous amphibole, which suggests that it may be related to the water-bearing fluid derived from the subducted plate.According to the discrimination of Ba/Th-La/Sm (Fig.17a), the BA has a high coupling with subduction plate dehydration and a high oxygen fugacity trend.Combined with the influence of fluid metasomatism shown by the BA in Fig.17b, it can be concluded that the source area of the BA obviously experienced fluid metasomatism after dehydration of subducted plate.Ta*and Nb*can be used as criteria to distinguish the mantle-derived characteristics of basalts.In the Ta*vs.Nb*diagram (Fig.17c), BA samples are mainly located near the PM and below the average continental crust (ACC), showing obvious clustering.One set of data suggests the addition of continental crust material, while the other set of data indicates that the primitive magma should have derived from the PM.Further input of BA data into [Hf/Sm]PM-[Ta/La]PMdiscrimination (Fig.17d) shows that the samples basically plot in the subduction metasomatism field, showing the characteristics of the source area influenced by both subduction melt and fluid.In conclusion, the melt (fluid) produced by the subduction of crust-derived sediments wrapped by oceanic plates may be the main factor leading to water-, volatile matter- and incompatible elementrich conditions in the BA source area.

Fig.16 Plots of Sr/Y vs.Y (μg/g) (after Defant and Drummond (1990)) (a) and La/Yb vs.Yb (μg/g) (after Defant and Drummond (1990)) (b)

Fig.17 Plots of Ba/Th vs.La/Sm (after Labanieh et al.(2012) and Li et al.(2022)) (a); U/Th vs.Th (μg/g) (after Sun et al.(2003)) (b); Ta* vs.Nb* (after Pearce et al.(2005)) (c); and [Hf/Sm]PM vs.[Ta/La]PM (after Pearce et al.(2005)) (d)

4.3.2 Genetic analysis of DA

Fig.18 Nb (μg/g) vs.10000Ga/Al diagram (after Whalen et al.(1987)) (a); (Na2O+K2O)/CaO vs.(Zr+Nb+Ce+Y) (μg/g)diagram (after Whalen et al.(1987)) (b); Rb/Nb vs.Y/Nb diagram (after Eby (1992)) (c) and Sc/Nb vs.Y/Nb diagram (after Eby (1990)) (d)

The chondrite-normalized REE pattern of the DA(Fig.8a) shows a gentle right-leaning pattern, with more enrichment of ΣREE (245.49-376.16 μg/g, with an average value of 284.36 μg/g) than those of UCC, a slight positiveδCe anomaly (1.01-1.07, 1.04 on average), a moderate negativeδEu (δEu=(Eusample/Euchondrite)/anomaly (0.60-0.70, 0.63 on average), and a high LREE/HREE ratio (2.40-2.92, 2.81 on average).The PM-normalized spider diagram (Fig.8b) shows curve characteristics similar to those of the UCC,with significant enrichments in LILEs Rb and K and HFSEs Th, U, Zr, Hf and REEs and depletions in LILEs Ba, Sr, and P and HFSEs Ta, Nb, and Ti,which shows significant geochemical characteristics of siliceous and aluminous crust.DA rocks have the characteristics of A-type granite (Fig.18a-b), which may derive from the alkali-rich magma source area (Pati?o Douce, 1997; Feio et al., 2012).Figure 18c-d shows that DA data basically plot into the A2-type region, suggesting that the magma source area of DA may have formed in a thinning crustal environment.However, there are some differences between DA and the traditional A-type granites in the “wing-type” distribution of the chondrite-normalized REE pattern, enrichment in Yb and K2O and strong negative δEu anomalies.The Nb/U ratio of DA (7.67 on average) is lower than that of MORB (47 on average) and that of the average crust (10 on average); and the Ce/Pb ratio of DA (2.62-8.50, 5.10 on average) is similar to that of ACC, indicating the addition of a certain mantle source material.Combined with the data of DA samples falling in the normal arc magma field(Fig.16a-b) and similar characteristics of island-arc Zhejiang and Fujian granites (Zhang et al., 2010)(characterized by moderate negativeδEu anomalies,Sr content and Yb content varying in 30-400 μg/g and 1-7 μg/g, respectively), it can be concluded that the magma source area of DA is a typical subducted island-arc andesitic magma produced in an environment of crustal thinning (Stevenson, 2018).The high (87Sr/86Sr)iratio of DA (0.709 39-0.711 29,0.710 35 on average), the less variedεNd(t) values(-6.2--4.8, -5.6 on average) andεHf(t) values (-2.9--1.7, -2.2 on average) (Supplementary Table S5), the isotopic values close to the young UC (Fig.14a) and the evolutionary trend toward EM-Ⅱ (Fig.14b)suggest that the magma source area was the product of partial melting of the lower crust, the magma chamber environment was relatively closed, and the element enrichment attribute was the result of the evolution of magma differentiation and crystallization(Li et al., 2009).

According to the analysis of the Harker diagram(Fig.10) and source area discrimination diagram(Fig.12), DA rocks have the characteristics of obvious mafic magma melting or biotite gneissderived melt, and the obvious crystallization differentiation attribute indicates that DA should be derived from the location of the lower crust and lithospheric mantle.Further analysis of the Pb isotope characteristics of DA shows that the data plot close to the EM-Ⅱ region (Fig.19a-b), showing obvious crust-mantle metasomatism.Combined with the EM-Ⅱ characteristics of DA in Fig.14, both suggest that the subduction of the western Pacific resulted in continuous asthenospheric underplating on the continental margin, and partial melting of the lithospheric mantle and lower crust formed the magma chamber in the DA source area.DA is highly similar to the Triassic granites of the SCC in isotopic geochemistry, which indicates an affinity between the Meiji Reef and SCC (Jahn, 2010; Li et al., 2018a; Hou et al., 2019).

Based on the analysis of rock geochemistry and isotopic geochemistry of BA and DA, the overall lithology changed gradually from the bottom of DA to the top of BA and DA alternately, with good continuity in the core feature.This may imply that the continuous fractional crystallization of pyroxene and olivine in the basic magma produced a large amount of silicate unsaturated alkaline melt, resulting in more obvious alkaline characteristics of the upper cycle.

The large fluctuationsεHf(t) values (-1.5-11.5,5.7 on average) and TDM(Ma) ages (808-3 434 Ma,1 415 Ma on average) larger than zircon U-Pb ages in BA, and the relatively stableεHf(t) values(-1.7--2.9, -2.2 on average) and TDM(Ma) age(1 128-1 238 Ma, 1 162 Ma on average) larger than zircon U-Pb age in DA (Supplementary Tables S4-S5),suggesting that a portion of crust-derived accretionary wedge wrapped by subduction in the western Pacific Ocean entered deep, a large number of subduction melt (fluid) formed with the intensification of subduction depth were rapidly added to the asthenosphere, forming enrichment mantle-derived magma chamber (EM-Ⅱ-typed) of the BA in the study area.As asthenospheric upwelling continued to underplate, the DA source magma chamber was formed by partial melting of the lithospheric mantle and lower crust, and the upwelling of magma intensified the gradual thinning of the crust on the continental margin.Subsequently,tectonic fractures appeared, and then the magma chamber in the DA source area was the first to erupt;magma from the deeper BA source area also rapidly erupted to form rocks along this conduit.The rhythmic interbedding revealed by the petrographic analysis also confirmed this view.BA and DA show a trend of decreasingεHf(t) values with the age of TDM(Ma), and the correlation is high in Fig.20a,which also indicates that the BA was influenced by obvious new crust, while DA received more ancient crust-derived materials (Fig.20b).However, the two samples with older model ages in BA indicate the addition of ancient crust-derived materials at the same time (Blichert-Toft and Albarède, 1999).

Fig.19 (207Pb/204Pb)i vs.(206Pb/204Pb)i diagram (after Zindler and Hart (1986)) (a); and (208Pb/204Pb)i vs.(206Pb/204Pb)i diagram(after Zindler and Hart (1986)) (b)

4.4 Tectonic spatio-temporal attribute

In addition to the eruption and overflow facies,the cryptoexplosive breccia facies appeared in the upper 1 315.3-1 292.5 m in the petrographic characteristics of the core samples from Well NK-1 of the Meiji Reef, indicating that the volcanic rock layer was thickly accumulated at this time and rich in fluid or gas components, resulting in magma emplacement to a shallow place and pressure eruption.According to the research results of Valentine et al.(2014), the depth of magma eruption at this time was relatively shallow (＜200 m) and there was some water environment interference.Their results are not only consistent with the microscopy identification results but also suggest that the magma eruption at this time occurred in a relatively open (extension and relaxation) environment.The lenticular welded tuff bands developed at 1 275.2 m, indicate that the volcanic eruption environment in the study area was already in the aquatic environment.The limited potassization,slight silicification, and late carbonation in the firstcycle core samples also indicate that the magma in the source region transformed into alkaline and volatile components, suggesting that the tectonic environment in the study area is only equivalent to the early stage of extrusion to tension transformation (Caricchi et al., 2014).

According to zircon U-Pb geochronology analysis of DA samples from 1 225.4 to 1 684.0 m in cores(Supplementary Table S4; Miao et al., 2021; Wei et al., 2022), the age of volcanic rocks in the Meiji Reef ranges from 218.6±3.2 Ma to 217.9±3.5 Ma.Because DA and BA were interbedded without obvious alteration, it is believed that DA and BA were the products of the magmatic eruption at the same time.The analysis of the tectonic environment of the DA rocks (Figs.21-22) shows that the DA volcanic rocks of the Meiji Reef should have been produced in the initial back-arc region of the subduction zone, and the tectonic environment was the transitional stage from extrusion to extension.

Based on the analysis of the tectonic environment of the BA (Fig.23) and the characteristics of the extensional environment on the terrigenous side shown by DA, it is concluded that the BA-produced environment of the Meiji Reef was a subduction arcback initial extensional basin.Combined with the extrusion-extension transition characteristics of DA in the study area, it is concluded that the overall tectonic environment of the Meiji Reef in the Late Triassic was more likely produced by the initial extensional environment (initial continental margin extensional rift environment) under the relative extrusion environment than the obvious back-arc extension characteristics (Meschede, 1986).This interpretation is consistent with the accordion-type periodic “open-close” dynamic background proposed by Yang et al.(2002), which is dominated by extensional tectonics, indicating that the magma production environment in the study area is a dynamic rhythmic change rather than an example of simple extensional tectonics.This interpretation is also consistent with the complex tectonic and geological background of the ancient continent in South China.

4.5 Mechanism of subduction on the western side of the ancient Pacific Ocean

Fig.20 εHf(t) vs.TDM (Ma) diagram (a) and εHf(t) frequency histogram (b) of volcanic rocks from Well NK-1

Fig.21 TiO2 (wt.%) vs.Al2O3 (wt.%) diagram (after Müller and Groves (2000)) (a); Sr (μg/g) vs.Yb (μg/g) diagram (after Kong et al.(2012)) (b); Ti/100 vs.La vs.Hf×10 ternary diagram (after Müller and Groves (2000)) (c); and Zr×3 vs.Nb×50 vs.Ce/P2O5 ternary diagram (after Müller and Groves (2000)) (d)

The SCC was affected by the Paleo-Pacific subduction during the Triassic period, but the study of subduction arc magmatic rocks during this period is obviously lacking, which makes the mechanism of Paleo-Pacific subduction controversial.The U-Pb isotopic geochronology data of volcanic rocks from Well NK-1 of the Meiji Reef in this paper show that magmatic eruptions occurred in the early Late Triassic, which provides important subduction time nodes and deep geochemical field data for the study of the ocean-continental transition zone between the SCC and the Western Pacific.The BA in the study area shows the magmatic characteristics of subducted island arcs with multisource properties, and subduction components should be added to the mantle wedge(Fig.14).Some previous studies attributed the changes in geochemical dynamics in the subduction zone to the mixture of three components, namely,DM, sediment melt, and aqueous fluid generated by dehydration of altered ocean crust (Pati?o Douce,1996).The enrichment of LILEs (such as Ba and Sr)and LREEs (such as La, Ce, and Nd) in arc magma has been confirmed by rock simulation experiments,as mainly caused by subduction components added to the mantle wedge (Pati?o Douce and Beard,1995).The content of Ba and Sr in BA is significantly higher than that of PM in the study area, indicating that there should be an obvious subduction component in the source area, which is closely related to subduction melt (fluid).According to the core depth (age) and multi-element and isotope geochemical analysis in Fig.24, the zircon U-Pb ages are basically the same (220.5-215 Ma),indicating that the volcanic rocks with a thickness of nearly 1 000 m should have been formed in a simultaneous eruption.In addition, the BA data in Fig.24 show significant fluctuations with the depth of the well column, especially the significantly fluctuating Ti/Ti*curve, which indicates that the BA source area is affected by both recirculating oceanic crust and recirculating sediments.In addition, the ratios of Nd/Sr and Ba/La in BA show the characteristics of sharp fluctuations at the same location (1 400-1 450 m and 1 200-1 300 m),indicating that the subduction material composition decreased first and then intensified during the eruption process.The consistent geochemical data characteristics of DA again indicate a stable and relatively closed environment in the source area.

Fig.22 Distribution of the volcanic rock data retrieved from Well NK-1 in Rb (μg/g) vs.(Yb+Ta) (μg/g) (a) and Rb (μg/g)vs.(Y+Nb) (μg/g) discrimination diagrams of Pearce et al.(1984) (b); R2 vs.R1 diagram (after Batchelor and Bowden (1985) (c) and lg[CaO/(Na2O+K2O)] vs.SiO2 (wt.%) diagram (after Maniar and Piccoli (1989) (d)

The obvious fluctuation data characteristics of the BA suggest the complexity of the source area.The tholeiitic index (THI) of the BA is 0.705,calculated using the THI proposed by Zimmer et al.(2010), which represents the relative depletion of FeO in magma.It should be the calc-alkaline magma series, which is consistent with the calculation range of global island-arc basalt.In the SiO2-TFeO/MgO diagram of Hora et al.(2009) (Fig.25), BA data mostly fall into the tholeiitic region below the CA(calc-alkaline)/TH (tholeiitic) index=1 line and are in the range of medium-Fe to high-Fe transitions,which is consistent with the global data of MORB and backarc basin basalt (BABB) and is in the primary stage of MORB and BABB evolution.However, it is worth noting that the obvious islandarc basalt (IAB) imprint of the BA geochemical properties may not support the true sense of MORB and BABB.This indicates that the BA has calcalkaline island-arc properties and develops into intraplate tholeiite.Combined with the calc-alkaline characteristics of the BA geochemical characteristics shown above, it is believed that it belongs to the late stage of island-arc magmatic development and has obvious mature island-arc characteristics, but it is different from the basalt of the back-arc basin and should be the product of early-stage back-arc extension.The study of global magma data showed that THI formed a very obvious inverse correlation with H2O, and the increase in water content in magma dominated the transition from calc-alkaline to intraplate tholeiitic magma.Furthermore, using the H2O content calculation formulas of Zimmer et al.(2010) and Kelley et al.(2010), the BA magmatic water content is 5.7 wt.% and 5.6 wt.%, respectively.The results are consistent in general with the evolution trend of the liquid lines of descent.The study of Tollan and Hermann (2019) showed that the dehydration effect of the subducted slab could form island-arc magmas with water-rich phases and high oxygen fugacity, which can metasomatize the overlying continental mantle wedge.However, H2O in aqueous solution can further produce H2and O2,which makes the mantle source more oxidized.

Fig.24 The relationship between core depth and age and multi-element and isotopic geochemistry from Well NK-1

Fig.25 SiO2 (wt.%) vs.TFeO/MgO plot (after Arculus (2004))

Based on this analysis and previous studies on the subduction model of the western side of the Paleo-Pacific Ocean, it is concluded that the Meiji Reef underwent obvious Pacific Plate slab subduction during the Late Triassic.A large number of previous studies have found that a large amount of terrigenous accretionary wedge material can be wrapped by oceanic plate subduction into the deep asthenosphere to form recirculating sediment melt(Turner et al., 2000; Yan et al., 2020); however, the oceanic crust can also form recirculating oceanic crust melt (Wang et al., 2002).Both of these materials can be added to the asthenosphere to cause changes in the magma chamber.This further reveals the characteristics of high oxygen fugacity, rich in incompatible elements and the oxidation environment of magma in the BA source area.However, DA in the study area shows obvious geochemical characteristics of high-K volcanic rocks, which further indicates that the magma was produced far from the arc magma center, and the magma source area is relatively stable, which is not the product of intense extrusion collision mode (Raymond et al.,2020; Gong et al., 2021a).Combined with the DA relaxed magmatic-tectonic environment indication and the BA marginal extensional rift environment discrimination, the subduction model in the study area is similar to the low-angle slab subduction model established by predecessors.The new young oceanic crust in the Paleo-Pacific Ocean at this time also provided the prerequisite for low-angle subduction.The relative extension “detachment” tectonic environment is also conducive to the continuous eruption of magma over a large area.Many crustderived materials are wrapped by Paleo-Pacific Ocean slab subduction.When slab subduction reached a certain depth and position, a large amount of melt (fluid) formed by continuous plate dehydration continued to weaken the lithospheric mantle and form EM-Ⅱ mantle characteristics,resulting in the unique complex geochemical properties of the BA in the study area.A large amount of lithospheric mantle and lower crust with the intensification of underplating began to partially melt to form the DA parent magma, and the gradual thinning of the crust formed a large number of tensile fractures.BA and DA magma chambers erupted along these fractures under pressure to form the Late Triassic volcanic rocks of the Meiji Reef.

Ishikawa and Tera (1997) studied a large number of underwater volcanoes in the Chiba arc and found that they may have formed in a rear-arc environment, which mainly produced a series of tholeiites, from calc-alkaline basalt to rhyolite,mainly basaltic andesite, which is clearly different from back-arc basalt.Peate et al.(1997) also pointed out that Vanuatu island-arc magma in the southwestern Pacific Ocean may be of the “reararc” type, which has a low degree of magma melting and is a K-containing volcanic rock series at the farthest end of the arc magma.The geochemical properties of the two series rocks in the study area are well consistent with the above characteristics,and the tectonic spatio-temporal properties also have rear-arc characteristics.Therefore, it is considered that the study area probably belongs to this obvious rear-arc tectonic environment, and the two series of rocks belong to the farthest part of the arc magma.Combined with the characteristics of geochronology, it is concluded that the western margin of the Paleo-Pacific Ocean began to subduct to the ancient land of South China during the Late Triassic, a mature arc magmatic system was formed, and back-arc extensional tectonics gradually formed.This extensional tectonics provided an incentive for the NB to split away from the ancient South China land,and the Late Triassic was very likely an important tectonic time point for detachment.It also indicates that the ancient South China land should have been a relatively relaxed passive continental margin at this time, and the western Pacific subduction mode should have been an Andean type (Pitcher et al., 1985).

The obvious extensional tectonic environment in the study area suggests that the SCC tectonics were not converted to the Paleo-Pacific tectonic domain at this time.Combined with the characteristics that the Changning-Menglian suture belt entered the postcontinental collision extension period and the Qinling-Dabieshan-Sulu collision belt entered the main collision period, it can be concluded that the southeastern margin of South China was still a passive collision during the Late Triassic period,which facilitated the subduction of the western Pacific Plate and resulted in the generation of a large number of back-arc rifts (Wu and Suppe, 2018).This has laid a foundation for the emergence of the detachment of a large number of islands in Southeast Asia (e.g., NB, Zhongsha Block, Mindoro Block, Palawan Block, etc.).During the Late Triassic, the NB was located near the back-arc extensional tectonics on one side of the continental margin (southeast margin of the SCC) under the subduction environment of the western Pacific Ocean and was subjected to obvious slab subduction of the Paleo-Pacific Ocean.

This “staged” slab subduction pattern of the Paleo-Pacific on the SCC indicates that the whole southeastern margin of the SCC was obviously influenced by the Paleo-Pacific tectonic domain during the Late Triassic, and the initiation of this large-scale tectonic domain had a profound impact on the whole SCC.For example, the large-scale mineralization in the SCC during the Late Triassic (high- and lowtemperature W-Sn-Au-Sb-Pb-Zn metallogenic belt located east to west in the SCC during 230-200 Ma),oceanic and continental interaction sedimentary tectonics and Paleo-climatic abrupt change, and other important events occurred during this period.A large number of initial continental margin rifts formed by lithospheric crust thinning caused by Paleo-Pacific Plate slab subduction provide tectonic space for a large number of enriched deposits.This dynamic rhythmic open-close model facilitates the continuous differentiation of magma (conducive to the collection of metallogenic elements) and intrusion.This is the main reason why BA and DA are interbedded in the study area.We believe that the back-arc extensional environment caused by the subduction of the Paleo-Pacific Plate in the study area was an important cause to the separation of the NB from the ancient South China block.The Late Triassic, an important tectonic event node, is an important time point for the formation of the initial intracontinental rift in the SCS.The two series of volcanic rocks in Well NK-1 of the Meiji Reef provide strong support for this theory.However, the multistage tectonic movement after this time point led to the “intermittent” adjustment of the shape and boundary of the NB and finally created the current complex strike-slip and sag pattern (Fig.26).

5 CONCLUSION

This detailed study on genetic mechanism of volcanic rocks in Well NK-1 of the Meiji Reef suggested the following:

(1) The volcanic rocks in Well NK-1 comprise 36 rhythmic layers encompassing an overflow phase and an intermittent eruption-overflow phase, and contain locally extreme cryptoexplosive breccia facies.These rocks can be further classified as a first mafic→neutral magmatic volcanic activity cycle (3 subcycles, in apparent thickness of 727.3 m) and a second mafic→intermediate-felsic magmatic volcanic activity cycle (including 3 subcycles, in apparent thickness of 294.0 m), which exhibits obvious characteristics of alkaline magma.

(2) Combining the geochronological and geochemical data of Well NK-1 published by our team, we can conclude that BA and DA were formed during the Late Triassic between 218.6±3.2 Ma and 217.9±3.5 Ma.The two series of rocks exhibit obvious island-arc magmatic properties and EM-Ⅱ-typed magmatic characteristics.The PM-normalized spider diagram shows that BA exhibited more similar characteristics to those of LCC and was significantly enriched in incompatible elements.However, DA was more consistent with the UCC and enriched in crustderived elements.The BA magma source area is partially melted, which is influenced by a large amount of subducted melt (fluid).The DA magma source area is dominated by fractional crystallization,and assimilation and contamination are very weak.The Sr-Nd-Pb-Hf isotopic characteristics of these two series rocks indicate obvious crust-mantle mixing.The BA source area reflects obvious influences of young oceanic crust and old crustal recycling, while the DA source area mainly indicates the addition of ancient crustal material.

(3) The geochemical characteristics of BA and DA are highly consistent with those of rear-arc magmas, showing enrichment characteristics of incompatible elements higher than those of MORB and lower than those of OIB, as well as enhanced isotopic properties.A large amount of sediment entrained by the Paleo-Pacific Ocean was slab subducted westward into the deep SCC, and continuous plate dehydration caused a large amount of melt (fluid) to join the asthenosphere on the side of the continent.Thus, the EM-Ⅱ-typed source area was formed, where a large number of water-rich,volatile-rich fluids and recirculating sediment melts were continuously enriched.

Combining the Late Triassic tectonic geological background in the study area and the two series of tectonic attribute analysis data in this paper, it is proposed that the southeastern margin of the SCC at this time was still a passive continental margin, the magma eruption environment was a relatively extruded extensional environment, which was a phased “periodic open-close” tectonic system, and a mature island-arc system formed along the southeastern margin of the ancient South China continent during this period.The unique accordion-like periodic open-close movement may have reflected the main tectonic dynamic mechanism of the separation of the NB from the SCC, and thus, the back-arc extensional tectonics may have induced the formation of the initial continental rift in the SCS.

6 DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available from Miao et al.(2021) (https://doi.org/10.1016/j.lithos.2021.106337) and Wei et al.(2022)(https://doi.org/10.1111/ter.12556).

7 ACKNOWLEDGMENT

We thank Xiuquan MIAO from Guangzhou Institute of Geochemistry, Chinese Academy of Sciences and Wu WEI from Institute of Geology and Geophysics, Chinese Academy of Sciences for their collaboration with our team and for their data contribution.We also thank Jianghong DENG and Bin DENG from Chengdu University of Technology for their help and advice in petrology and petrography.We are grateful to the Well NK-1 research team of the Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, for providing the core samples.