Hostname: page-component-7b9c58cd5d-7g5wt Total loading time: 0 Render date: 2025-03-16T03:07:51.599Z Has data issue: false hasContentIssue false

Provenance and depositional history of the Mesozoic Sanjiang Basin (northeastern China): implications for the uplift history of the northeastern Asian continental margin

Published online by Cambridge University Press:  03 November 2021

Jin-Peng Luan
Affiliation:
Key Laboratory of Ministry of Education for Safe Mining of Deep Metal Mines, Northeastern University, Shenyang 110819, China School of Earth and Environmental Sciences, University of Queensland, Brisbane 4072, Queensland, Australia
Gideon Rosenbaum
Affiliation:
School of Earth and Environmental Sciences, University of Queensland, Brisbane 4072, Queensland, Australia
Jian-Fei Fu*
Affiliation:
Key Laboratory of Ministry of Education for Safe Mining of Deep Metal Mines, Northeastern University, Shenyang 110819, China
*
Author for correspondence: Jian-Fei Fu, Email: fujianfei@mail.neu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

The uplift history of northeastern Asian continental margin has been the subject of much debate. The Sanjiang Basin is an ideal area to investigate the uplift history of this margin. We present new data from sandstone petrology, whole-rock geochemistry, U–Pb geochronology and Hf isotopy of detrital zircons to trace the provenance of the basin and to unravel the uplift history of the northeastern Asian continental margin. We investigated the chemical and plagioclase indices of alteration and index of compositional variability (CIA, PIA and ICV, respectively) in samples from two formations. The geochemical proxies indicate that source rocks were subjected to an intensified weathering process. Based on the high SiO2/Al2O3 values and low K2O, we infer that the two basinal strata experienced silicification and insignificant potassium metasomatism. Sandstone petrography is indicative of low degrees of sedimentary sorting, suggesting a proximal deposition. The provenance fingerprints of light rare earth, high-field-strength and transition metal elements indicate that parts of the provenance included recycled sediments, and that first-cycle sediments were derived mainly from felsic rocks with a minor contribution from intermediate and mafic rocks. The combined detrital zircon U–Pb ages and Hf isotopic results constrain the contributions of different source terranes over time. The provenance results, in combination with seismic profile interpretations from the Sanjiang Basin, suggest that the northeastern Asian continental margin experienced at least three stages of uplift, which were driven by subduction initiation of the Paleo-Pacific Ocean during the Jurassic Period, a plate motion change during the Early Cretaceous Epoch, and a shallowing of the slab dip angle during the Late Cretaceous Epoch.

Type
Original Article
Copyright
© The Author(s), 2021. Published by Cambridge University Press

1. Introduction

Detrital zircon U–Pb geochronology is commonly used for provenance analyses, providing information on sedimentary dispersal, basin evolution, palaeo-drainage systems, surface uplift and the tectonic setting during deposition (Rainbird et al. Reference Rainbird, Hamilton and Young2001; Kosler & Sylvester, Reference Kosler and Sylvester2003; Barth et al. Reference Barth, Wooden, Jacobson and Probst2004; Lawton et al. Reference Lawton, Bradford, Vega, Gehrels and Amato2009; Cawood et al. Reference Cawood, Hawkesworth and Dhuime2012; Shaanan et al. Reference Shaanan, Rosenbaum and Campbell2019). This methodology has been applied extensively in studies that aimed to understand the evolution of Mesozoic–Cenozoic basins in northeastern Asia, including the Sanjiang Basin (Fig. 1a; e.g. Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015, Reference Zhang, Dilek, Chen, Yang and Meng2017). Despite the power of detrital zircon U–Pb dating, the approach has some inherent limitations. Firstly, its ability to detect mafic and ultramafic sources is limited due to the scarcity of zircon grains in such rocks. Secondly, detrital zircon U–Pb ages cannot distinguish between coeval igneous rocks that have been derived from different magmatic sources (Reiners et al. Reference Reiners, Campbell, Nicolescu, Allen, Hourigan, Garver, Mattinson and Cowan2005). Thirdly, the resilience of zircon grains and their ability to survive through multiple episodes of sedimentary recycling and transportation (Bruguier & Lancelet, Reference Bruguier and Lancelet1997; Cherniak & Watson, Reference Cherniak and Watson2000) means that the provenance of older detrital zircon grains (e.g. Neoproterozoic zircons) could be elusive and ambiguous. These limitations make it difficult to decipher the depositional history of the Sanjiang Basin and surface processes (e.g. uplift history) of the northeastern Asian continental margin.

Fig. 1. Topographic and tectonic map of (a) northeastern Asia and (b) northeastern China (after Tao et al. Reference Tao, Niu, Ning, Chen, Grand, Kawakatsu, Tanaka, Obayashi and Ni2014); (c) aeromagnetic anomaly map of northeastern China; and (d) geological map of northeastern China showing the location of the Sanjiang Basin and the ages of igneous rocks (after Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Xu et al. Reference Xu, Pei, Wang, Meng, Ji, Yang and Wang2013).

The Sanjiang Basin is a large Mesozoic sedimentary basin overlying the basement of the Jiamusi Massif in the easternmost Central Asian Orogenic Belt (CAOB; Fig. 1a, b; Zhang et al. Reference Zhang, Chen, Yang, Feng, Wu, Batt, Zhao, Sun, A, Wang and Yang2012). Published detrital zircon U–Pb data indicate that the detrital material in the Sanjiang Basin was likely derived from multiple sources of similar ages (e.g. Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015). While previous detrital zircon ages have provided valuable information on the provenance of some basinal strata (Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015), the provenance and depositional history of the Sanjiang Basin remain controversial, partly because of the inherent limitations of detrital zircon U–Pb geochronology. Understanding the provenance and depositional history of the Sanjiang Basin is crucial for reconstructing the uplift history of the northeastern Asian continental margin. To obtain a more diagnostic signature of the provenance region, complementary data from clastic sedimentary rocks and individual zircon grains must be used. Such data can be obtained, for example, from sandstone petrology (Dickinson & Suczek, Reference Dickinson and Suczek1979; Dickinson, Reference Dickinson and Zuffa1985), whole-rock geochemistry (e.g. Migani et al. Reference Migani, Borghesi and Dinelli2015), trace-element composition of zircon grains (e.g. Belousova et al. Reference Belousova, Griffin, O’Reilly and Fisher2002; Rubatto, Reference Rubatto2002; Campbell et al. Reference Campbell, Rosenbaum, Allen and Mortimer2020), zircon morphology (e.g. Pidgeon et al. Reference Pidgeon, Nemchin and Hitchen1998; Shaanan & Rosenbaum, Reference Shaanan and Rosenbaum2018) and zircon Hf isotope geochemistry (Hawkesworth & Kemp, Reference Hawkesworth and Kemp2006). Sandstone petrology can provide information on the hydrodynamic conditions during sedimentary transport, as well as the tectonic setting of the source regions (Dickinson, Reference Dickinson and Zuffa1985). Whole-rock geochemistry of the sedimentary rocks can provide information on the weathering history, sedimentary recycling and the tectonic setting of the source terranes (Roser & Korsch, Reference Roser and Korsch1986; Armstrong-Altrin et al. Reference Armstrong-Altrin, Lee, Verma and Ramasamy2004). Hf isotope analyses can provide constraints on crustal recycling, with derived model Hf ages constraining the timing of extraction from a mantle reservoir.

In this paper, we present sandstone petrography and whole-rock geochemistry results from the Sanjiang Basin, complemented by detrital zircon U–Pb ages and in situ Hf isotope data. Combined with available data from exploratory drilling and seismic profiles, our results elucidate the provenance and sedimentary history of the Sanjiang Basin and provide a new insight into the uplift history of the northeastern Asian continental margin.

2. Geological setting

The study area in northeastern China belongs to the eastern part of the CAOB and is bounded by the Siberian and North China cratons (Fig. 1a, b; Kröner et al. Reference Kröner, Kovach, Belousova, Hegner, Armstrong, Dolgopolova, Seltmann, Alexeiev, Hoffmann, Wong, Sun, Cai, Wang, Tong, Wilde, Degtyarev and Rytsk2014; Xiao et al. Reference Xiao, Windley, Sun, Li, Huang, Han, Yuan, Sun and Chen2015). The basement geology in this area includes a number of Proterozoic massifs and terranes, such as the Erguna Massif, Xing’an Massif, Songnen Massif, Jiamusi Massif and Nadanhada Terrane (Fig. 1b; Sengör et al. Reference Sengör, Natal’in and Burtman1993; Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Han et al. Reference Han, Zhou, Li and Song2017; Luan et al. Reference Luan, Wang, Xu, Ge, Sorokin, Wang and Guo2017 a, b, Reference Luan, Yu, Yu, Cui and Xu2019; Yang et al. Reference Yang, Ge, Bi, Wang, Tian, Dong and Chen2018). During the Mesozoic Era, the area was subjected to widespread crustal extension, indicated by abundant fault-bounded sedimentary basins, metamorphic core complexes and intraplate magmatism (Wang et al. Reference Wang, Zheng, Zhang, Zeng, Donskaya, Guo and Li2011; Zhang et al. Reference Zhang, Yan and Ji2019). Major graben-hosted basins in the area include the Sanjiang, Songliao, Erlian and Hailar basins (Fig. 1b). The easternmost basin, Sanjiang Basin, is located close to the Circum-Pacific subduction zone and has been regarded as a typical Mesozoic retro-arc foreland basin (Zhang et al. Reference Zhang, Dilek, Chen, Yang and Meng2017).

The Sanjiang Basin is situated along the China–Russia border, with its northern part (in Russia) called Amur Basin (Fig. 1b). The basin is surrounded by Mesozoic orogens (Lesser Xing’an Range, Zhangguangcai Range, Middle Jilin Orogen) and the Nadanhada Terrane (Fig. 1b). The area includes abundant igneous plutons of various ages (Fig. 1c, d), including early Palaeozoic (Wilde et al. Reference Wilde, Zhang and Wu2000, Reference Wilde, Wu and Zhang2003; Bi et al. Reference Bi, Ge, Yang, Zhao, Yu, Zhang, Wang and Tian2014 a; Yang et al. Reference Yang, Ge, Zhao, Bi, Wang, Dong and Xu2017), Permian (Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Bi et al. Reference Bi, Ge, Zhang, Yang and Wang2014 b, Reference Bi, Ge, Yang, Zhao, Xu and Wang2015; Yang et al. Reference Yang, Ge, Zhao, Yu and Zhang2015 b) and Triassic (Yang et al. Reference Yang, Ge, Zhao, Dong, Xu, Wang, Ji and Yu2015 a, b; Guo et al. Reference Guo, Xu, Yu, Wang, Tang and Li2016). A few Neoproterozoic and Early Cretaceous igneous rocks also occur in the region (Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Luan et al. Reference Luan, Xu, Wang, Wang and Guo2017 b, Reference Luan, Yu, Yu, Cui and Xu2019; Yang et al. Reference Yang, Ge, Zhao, Bi, Wang, Dong and Xu2017, Reference Yang, Ge, Bi, Wang, Tian, Dong and Chen2018). The Mashan, Heilongjiang and Yuejinshan complexes are three metamorphic complexes in the study area (Fig. 1d).

Based on aeromagnetic and seismic data, the Sanjiang Basin can be divided into two sub-basins, separated by the so-called Fujin Uplift (Figs 1c, 2; Zhang et al. Reference Zhang, Chen, Yang, Feng, Wu, Batt, Zhao, Sun, A, Wang and Yang2012, Reference Zhang, Dilek, Chen, Yang and Meng2017). The basin fill includes Mesozoic and Cenozoic terrestrial strata, consisting of volcanic, volcaniclastic, alluvial-fan, fluvial and lacustrine sediments (Fig. 2; Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015). The thickness of the sedimentary succession is variable, reaching a maximum of c. 9 km in the depocentre. Mesozoic strata include, in ascending order, the Suibin, Dongrong, Didao (K1d), Chengzihe (K1ch), Muling (K1m), Dongshan (K1ds), Hailang (K2hl), Qixinghe (K2qx) and Yanwo (K2yw) formations (Fig. 2). Cenozoic strata include the Baoquanling (Eb) and Fujin (Ef) formations. Exploratory drilling identified five unconformities within the Mesozoic–Cenozoic stratigraphic succession (Fig. 2; Zhang et al. Reference Zhang, Chen, Yang, Feng, Wu, Batt, Zhao, Sun, A, Wang and Yang2012).

Fig. 2. Cross-section and stratigraphy of the Sanjiang Basin (Zhang et al. Reference Zhang, Dilek, Chen, Yang and Meng2017). Locations of seismic profiles are shown in Figure 1c.

3. Samples and analytical methods

We measured two sections near Fujin (Fig. 1d). Section I (47° 19ʼ 37ʼʼ N, 132° 28ʼ 18ʼʼ E) is dominated by grey and black mudstone, sandstone and pebbly sandstone, with thin siltstone and silty mudstone layers (Fig. 3). Rocks in Section I contain bivalve and plant fossils, suggesting a paralic depositional environment. Based on the fossil assemblage, we consider that Section I belongs to the Dongrong Formation (Fig. 2). Section II (47° 12ʼ 42ʼʼ N, 132° 13ʼ 02ʼʼ E) consists of green mudstone, grey siltstone and medium sandstone, as well as interbedded grey mudstone and silty mudstone (Fig. 3). Based on the palynological record from this section (Aquiapollenites sp., Cyathidites sp., Cicatricosisporites sp. and Protoconiferus), we consider that Section II belongs to the Yanwo Formation (Fig. 2). The Yanwo Formation is intruded by 56–50 Ma granites (Fig. 4a, b; Li et al. Reference Li, Liu, Li, Zhao, Shi, Li and Yang2019). We collected 26 samples from the two sections for petrographic and whole-rock geochemical analysis. In addition, three samples were used for zircon U–Pb and Hf isotopic analysis. Mesozoic plant fossils (Cladophlebis and Podozamites) were found in two samples from the Dongrong Formation (Fig. 4c, d).

Fig. 3. (a) Detailed geological map of study area; and (b) measured geological sections of the Dongrong and Yanwo formations.

Fig. 4. Photographs showing (a) intrusive relationship; and (b) 50–55 Ma granite. Photomicrograph of (c) Cladophlebis fossil; (d) Podozamites fossil; (e) selected sandstone samples (cross-polarized light); (f) andesitic fragment; (g) felsic fragment; and (h) metamorphosed felsic lithic fragment. Lv – volcanic lithic fragment; Lm – metamorphic lithic fragment.

Six thin-sections from sandstone samples were selected for petrographic analysis using a microscope with an attached point-counting stage. In each sample, at least 300 points were counted following the Gazzi–Dickinson method (Ingersoll et al. Reference Ingersoll, Fullard, Ford, Grimm, Pickle and Sares1984; Dickinson, Reference Dickinson and Zuffa1985). Recorded components include monocrystalline quartz (Qm), polycrystalline quartz (Qp), plagioclase (P), alkali feldspar (Kf) and lithic fragments (L). Matrix and cement were not counted but estimated using standard charts for visual percentage estimation. The results are presented in online Supplementary Table S1 (available at http://journals.cambridge.org/geo and https://data.mendeley.com/datasets/nfy7gxg38p/1).

Samples for whole-rock analysis were crushed in a c. 200 μm mesh using an agate mill. The samples were then digested in Teflonbombs and measured by X-ray fluorescence (Rigaku RIX 2100 spectrometer), fused-glass disks and an inductively coupled plasma mass spectrometer (ICP-MS) (Agilent 7500a). Standard samples used were BHVO-1 (basalt), BCR-2 (basalt) and AGV-1 (andesite). Analytical results for these standards indicate that the analytical uncertainties are less than 3%, and the analytical precision for major and trace elements is better than 5% and 10%, respectively (Rudnick et al. Reference Rudnick, Gao, Ling, Liu and Mc Donough2004). The whole-rock geochemical analyses were conducted at the State Key Laboratory of Geological Processes and Mineral Resources at the China University of Geosciences (Wuhan, China). Major- and trace-element results are provided in online Supplementary Table S2 (available at http://journals.cambridge.org/geo and https://data.mendeley.com/datasets/nfy7gxg38p/1).

Zircon grains were extracted by standard density and magnetic separation techniques, followed by handpicking under a binocular microscope. Cathodoluminescence (CL) images were obtained to reveal internal structures within the zircon grains. Zircon U–Pb analyses were performed using an Agilent 7500a ICP-MS equipped with a 193 nm laser. The analyses created a crater diameter of 32 μm. For further information on the analytical procedures, see Yuan et al. (Reference Yuan, Gao, Liu, Li, Günther and Wu2004). Data reduction and correction for common Pb were done following the methods by Anderson (Reference Anderson2002), Ludwig (Reference Ludwig2003) and Liu et al. (Reference Liu, Gao, Hu, Gao, Zong and Wang2010). In situ zircon Hf isotope analyses were conducted using a Neptune multi-collector (MC-) ICP-MS. By modifying the X skimmer cone during the experiment, the signal intensities of Yb, Lu and Hf were improved by factors of 5.3, 4.0 and 2.4, respectively. The crater diameter of the single spot ablation was 44 μm. For further information about the operating conditions, see Hu et al. (Reference Hu, Liu, Gao, Hu, Dietiker and Günther2008, Reference Hu, Liu, Gao, Liu, Yang, Zhang, Tong, Lin, Zong, Li, Chen and Zhou2012). Zircon U–Pb dating and Hf isotope analyses were performed at Wuhan SampleSolution Analytical Technology Co., Ltd. (Wuhan, China). Zircon U–Pb ages and Hf isotope results are provided in online Supplementary Tables S3 and Table S4 (available at http://journals.cambridge.org/geo and https://data.mendeley.com/datasets/nfy7gxg38p/1), respectively.

4. Results

4.a. Sandstone petrography

All sandstone samples are arkoses and wackes (Fig. 5; online Supplementary Table S1), composed mainly of monocrystalline quartz, plagioclase and volcanic lithic fragments with minor polycrystalline quartz, K-feldspar and metamorphic lithic fragments (Fig. 4e–h). Quartz grains are moderately sorted and are sub-rounded to sub-angular (Fig. 4e). Feldspar grains are sub-rounded, and some plagioclase grains are altered into sericite (Fig. 4e). Volcanic lithic fragments consist of andesitic (Fig. 4f) and felsic fragments (Fig. 4g), as well as minor amounts of metamorphosed felsic lithic fragments (Fig. 4h). Biotite grains were observed in the wacke samples. Based on the degree of sorting and mineral morphology, the sandstone samples are considered compositionally and texturally immature.

Fig. 5. Lithological classification of siliciclastic rocks (Herron, Reference Herron1988). The grey symbols correspond to samples from Cao et al. (Reference Cao, Guo, Shan, Ma and Sun2015) and Hu et al. (Reference Hu, Meng and Liu2021).

4.b. Whole-rock geochemistry

Samples from the Dongrong and Yanwo formations show variable major-element compositions (online Supplementary Table S2). Medium- and fine-grained wackes and arkoses are characterized by relatively high concentrations of SiO2 (64.2–74.5 wt%), CaO (0.80–5.20 wt%), Na2O (2.28–3.83 wt%) and K2O (3.17–3.93 wt%), and low concentrations of Al2O3 (12.1–14.9 wt%), total Fe2O3 (1.65–3.27 wt%), TiO2 (0.30–0.66 wt%) and P2O5 (0.08–0.19 wt%). In contrast, siltstones and mudstones have relatively low concentrations of SiO2 (40.0–61.7 wt%), total Fe2O3 (1.90–6.50 wt%), MgO (0.70–1.42 wt%), K2O (0.22–1.62 wt%) and MnO (0.01–0.02 wt%); high concentrations of Al2O3 (14.6–24.8 wt%) and Na2O (1.30–3.60 wt%); and variable concentrations of CaO (0.50–2.40 wt%), TiO2 (0.66–1.15 wt%) and P2O5 (0.04–0.31 wt%). Samples 14JP1 and 14JP2, from the Dongrong Formation yielded chemical indices of alteration (CIA) of 43.2–84.8, plagioclase indices of alteration (PIA) of 42.3–85.6, and indices of compositional variability (ICV) of 0.34–1.12 (Fig. 6; online Supplementary Table S2). Sample 14JP3 from the Yanwo Formation yielded CIA, PIA and ICV values of 44.2–83.6, 43.5–84.3 and 0.30–1.05, respectively (Fig. 6; online Supplementary Table S2).

Fig. 6. (a) CIA–ICV diagram for sedimentary rocks from the Sanjiang Basin (Nesbitt & Young, Reference Nesbitt and Young1989; Cox et al. Reference Cox, Lowe and Cullers1995). (b) A–CN–K diagram for sedimentary rocks from the Sanjiang Basin (Fedo et al. Reference Fedo, Nesbitt and Young1995). Symbols as for Figure 5. The grey symbols correspond to samples from Cao et al. (Reference Cao, Guo, Shan, Ma and Sun2015) and Hu et al. (Reference Hu, Meng and Liu2021).

All samples are enriched in light rare earth elements (LREE) and depleted in heavy rare earth elements (HREE), with chondrite-normalized La/Yb (La/YbN) values of 4.8–14.1 (Fig. 7a; online Supplementary Table S2). The samples show various degrees of negative Eu anomaly (Eu/Eu* = 0.30–0.75), and compositional variations in both high-field-strength elements (HFSE) and large-ion-lithophile elements (LILE). Most samples are enriched in HFSE (e.g. Th, Zr and Hf) and depleted in LILE (e.g. Ba, K and Sr), P and Ti. Some samples (e.g. 14JP2 and 14JP3) show variable concentrations of Ba, K, La and Zr (Fig. 7b; online Supplementary Table S2). All samples have relatively low concentrations of transition metals, such as Cr, Ni, Cu and Sc (online Supplementary Table S2).

Fig. 7. (a) Chondrite-normalized REE patterns and (b) primitive mantle (PM) normalized trace-element diagrams. The chondrite and PM values used for normalization are from Boynton (Reference Boynton and Henderson1984) and Sun & McDonough (Reference Sun, McDonough, Saunders and Norry1989), respectively.

4.c. Detrital zircon UPb geochronology and Hf isotope analysis

From each of the three samples selected for U–Pb dating, more than 100 detrital zircon grains were dated. The results yielded 320 concordant ages. Out of these zircon grains (with determined U–Pb ages), 133 grains were analysed for Hf isotopic compositions. The results are presented in online Supplementary Tables S3 and S4. 206Pb/238U ages are used for grains younger than 1000 Ma. For statistical purposes, ages with discordance of > 10% were excluded. Most detrital zircon grains are euhedral–subhedral and display banded structures or oscillatory growth zoning (Fig. 8a–c). They are generally characterized by high Th/U ratios (≥ 0.1), with the exception of four grains that yielded Th/U ratios of 0.050.08. As evidenced by the zircon morphology and the high Th/U ratios (online Supplementary Table S3), most of the detrital zircon grains are interpreted as originating from a magmatic source (Belousova et al. Reference Belousova, Griffin, O’Reilly and Fisher2002; Rubatto, Reference Rubatto2002).

Fig. 8. Representative cathodoluminescence (CL) images of selected zircons and zircon U–Pb concordia diagrams of the Dongrong and Yanwo formations. Red and dotted circles on the zircon grains indicate locations of the analysis spots used for LA-ICP-MS U–Pb dating and Hf isotopic analyses, respectively.

4.c.1. Sample 14JP1-1

A total of 106 analytical spot analyses were conducted on 106 grains from sample 14JP1-1. The results yielded five main concordant 206Pb/238U age populations: 548493 Ma (n = 38), c. 454 Ma (mean square weighted deviation (MSWD) = 0.01, n = 6), 291254 Ma (n = 45), 241227 Ma (n = 4) and 208 ± 2 Ma (MSWD = 0.55, n = 6) (online Supplementary Table S3; Fig. 8d). In addition, seven grains yielded 206Pb/238U ages of 929 ± 15, 895 ± 10, 847 ± 11, 726 ± 8, 369 ± 5, 358 ± 5 and 336 ± 4 Ma. The youngest concordant detrital zircon age is 206 ± 3 Ma.

The ϵHf(t) values are divided into different groups based on the U–Pb ages of the zircon grains (online Supplementary Table S4). Cambrian zircon grains (548493 Ma) yielded ϵHf(t) values between −5.1 and +2.3, with corresponding two-stage model ages (T DM2) of 18241333 Ma (Fig. 9; online Supplementary Table S4). Ordovician zircon grains (c. 455 Ma) yielded ϵHf(t) values of 2.6–3.5, with corresponding T DM2 of 12851207 Ma (Fig. 9; online Supplementary Table S4). Permian zircon grains (291254 Ma) yielded ϵHf(t) values between −2.4 and +4.7, with corresponding T DM2 of 14411238 Ma (Fig. 9; online Supplementary Table S4). Late Triassic (208–206 Ma) zircon grains yielded ϵHf(t) values of between +1.0 and +3.6, with corresponding T DM2 of 11771017 Ma (Fig. 9; online Supplementary Table S4).

Fig. 9. Correlations between Hf isotopic compositions and ages of detrital zircons from the Dongrong and Yanwo formations. J – Jurassic; T – Triassic; P – Permian; C – Carboniferous; D – Devonian; O – Ordovician.

4.c.2. Sample 14JP2-1

A total of 114 analytical spot analyses were conducted on 114 grains (online Supplementary Table S3; Fig. 8). The results yielded four major concordant 206Pb/238U age populations: 548 ± 3 Ma (MSWD = 0.91, n = 15), 524451 Ma (n = 52), 298–257 Ma (n = 26) and 242 ± 3 Ma (MSWD = 2.4, n = 8). Additional minor age populations are: 750 ± 14 Ma (MSWD = 0.12, n = 3), 422 ± 6 Ma (MSWD = 0.01, n = 2), 223 ± 5 Ma (MSWD = 0.16, n = 2) and 201 ± 3 Ma (MSWD = 0.18, n = 3). Three grains yielded 206Pb/238U ages of 886 ± 12, 387 ± 5 and 353 ± 5 Ma. The youngest concordant age is 199 ± 3 Ma.

Neoproterozoic zircon grains dated at c. 886 Ma and 755–747 Ma yielded ϵHf(t) values of c. −4.3 and 4.7–7.0, respectively (Fig. 9; online Supplementary Table S4). Zircon grains dated at c. 548 Ma yielded ϵHf(t) values ranging between −6.3 and +3.0, with corresponding T DM2 of 18971317 Ma (Fig. 9; online Supplementary Table S4). Cambrian (505–499 Ma) zircon grains yielded ϵHf(t) values between −6.5 and −0.9, with corresponding T DM2 of 18821557 Ma (Fig. 9; online Supplementary Table S4). Ordovician (478–456 Ma) zircon grains yielded ϵHf(t) values between −0.3 and +2.1, with corresponding T DM2 of 14871302 Ma (Fig. 9; online Supplementary Table S4). Silurian (c. 421 Ma) zircon grains yielded ϵHf(t) and TDM2 values of 13 and 539 Ma, respectively. The two Devonian zircon grains (387 and 353 Ma) yielded ϵHf(t) values of 11.6 and 5.8, respectively (Fig. 9; online Supplementary Table S4). Permian (297–249 Ma) zircon grains yielded ϵHf(t) values between −3.8 and +4.8, with corresponding T DM2 of 1528–1001 Ma (Fig. 9; online Supplementary Table S4). Triassic (246–199 Ma) zircon grains yielded ϵHf(t) values between −2.4 and +1.6, with corresponding T DM2 of 1431–1169 Ma (Fig. 9; online Supplementary Table S4).

4.c.3. Sample 14JP3-1

A total of 100 analytical spot analyses were conducted on 100 grains. The results yielded seven major concordant 206Pb/238U age populations: 446 ± 5 Ma (MSWD = 1.06, n = 5), 299254 Ma (n = 33), 252–235 Ma (n = 16), 208 ± 2 Ma (MSWD = 0.55, n = 6), 172 ± 4 Ma (MSWD = 1.7, n = 5), 95 ± 4 Ma (MSWD = 2.1, n = 5) and 85 ± 1 Ma (MSWD = 1.1, n = 16) (online Supplementary Table S3; Fig. 8c). Three additional minor age populations are 836 ± 12 Ma (MSWD = 0.33, n = 3), 367 ± 12 Ma (MSWD = 0.33, n = 2) and 121 ± 2 Ma (MSWD = 0.33, n = 3). Six grains yielded 206Pb/238U ages of 898 ± 9, 497 ± 6, 473 ± 6, 307 ± 4, 155 ± 5 and 131 ± 5 Ma. The youngest concordant detrital zircon age is 84 ± 1 Ma.

Neoproterozoic (c. 841 Ma) zircon grains are characterized by a wide variation of ϵHf(t) values (−4.7 to 10.7), and corresponding T DM2 range from 2013 to 1033 Ma (Fig. 9; online Supplementary Table S4). Ordovician–Devonian zircon grains are characterized by a narrower distribution of ϵHf(t) values (between −8.2 and −6.5), except for a single 367 Ma zircon that has an ϵHf(t) value of −0.4 (Fig. 9; online Supplementary Table S4). Permian (299–252 Ma) zircon grains, which are the most dominant zircon age population in sample 13JP3-1, yielded a wide variation of ϵHf(t) values ranging from 3.0 to 9.9, with corresponding T DM2 of 1475–647 Ma (Fig. 9; online Supplementary Table S4). Triassic–Cretaceous zircon grains are characterized by positive ϵHf(t) values. Triassic grains yielded ϵHf(t) values of 2.4–11.4, and their T DM2 ranged from 1440 to 552 Ma. Jurassic and Cretaceous zircons yielded ϵHf(t) values of 4.7–13.3 and 4.5–10.8; their corresponding T DM2 are 917–343 Ma and 875–487 Ma, respectively (Fig. 9; online Supplementary Table S4).

5. Discussion

5.a. Timing of deposition

Our results provide new constraints on the sedimentary history of the Sanjiang Basin. Previous detrital zircon geochronological studies have focused mainly on Cretaceous petroleum- and coal-bearing strata collected from boreholes in the Sanjiang Basin (Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015), but the depositional ages of exposed outcrops have remained poorly constrained. In fact, the two sections studied here, in the Fujin Uplift, have previously been considered as lower Permian strata (HBGMR, 1993). Our geological investigation and stratigraphic correlations indicate that these rocks belong to the Dongrong and Yanwo formations, respectively (Figs 1d, 2, 3). Furthermore, the discovery of Cladophlebis and Podozamites fossil assemblages in the Dongrong Formation (Fig. 4a, b) indicates that the age of the Dongrong Formation must be Mesozoic (Ash, Reference Ash1993; Falcon-Lang, Reference Falcon-Lang2002; Vavrek & Larsson, Reference Vavrek and Larsson2007). Our new detrital zircon data support this conclusion.

Sample 14JP1 from the lower part of the Dongrong Formation yielded a youngest concordant age peak of 208 ± 2 Ma. Sample 14JP2, from the middle part of the same formation, yielded a youngest concordant age peak of 201 ± 3 Ma. These Late Triassic ages provide a maximum depositional age constraint for the Dongrong Formation.

Sample 14JP3, from the bottom of the Yanwo Formation, yielded a youngest concordant age peak of 85 ± 1 Ma. This age provides a maximum depositional age constraint for this formation. The Yanwo Formation is intruded by a granodiorite body dated 54 ± 1 Ma (Fig. 4a, b; Wang et al. Reference Wang, Yang, Ge, Bi, Zhang and Xu2016a; Li et al. Reference Li, Liu, Li, Zhao, Shi, Li and Yang2019). The Yanwo Formation must therefore have been deposited after 85 Ma and before 54 Ma.

5.b. Depositional history

The depositional history of clastic sedimentary rocks typically involves weathering (Nesbitt & Young, Reference Nesbitt and Young1982), hydrodynamic sorting during transport (McLennan et al. Reference McLennan, Hemming, Mcdaniel, Hanson, Johnsson and Basu1993; Cullers, Reference Cullers1994), sedimentary recycling (Long et al. Reference Long, Sun, Yuan, Xiao and Cai2008) and diagenesis (Fedo et al. Reference Fedo, Nesbitt and Young1995, Reference Fedo, Eriksson and Krogstad1996). These processes control the chemical composition of detrital sediments. Using our petrographic and geochemical data, we discuss the effects of the above geological processes during the deposition of the Sanjiang Basin.

5.b.1. Weathering and sedimentary recycling

CIA and PIA are geochemical proxies that can be used to quantify the intensity of weathering in the source region (Nesbitt, Reference Nesbitt1979; Fedo et al. Reference Fedo, Nesbitt and Young1995). CIA values from the Dongrong Formation (43.2–84.8, with higher values in the upper part; online Supplementary Table S2) are generally higher than those found in the upper continental crust (UCC; CIA = 48; Rudnick & Gao, Reference Rudnick, Gao, Holland and Turekian2003), but lower than post-Archean Australian shale (PAAS; CIA = 70; McLennan et al. Reference McLennan, Hemming, Mcdaniel, Hanson, Johnsson and Basu1993). The results indicate that weathering was intensified during the deposition of the Dongrong Formation. PIA values from the same samples range from 42.3 to 85.6, and also show higher values in the upper part of the formation (online Supplementary Table S2), supporting our conclusion. Similarly, variations in the ICV values also indicate that the degree of weathering has increased from bottom to top. ICV values in sample 14JP1, which are generally higher than PAAS (ICV = 0.84; Fig. 6a), are typical values for an early stage of weathering (Fedo et al. Reference Fedo, Nesbitt and Young1995). In contrast, the lower ICV values (Fig. 6a) and relatively higher Al2O3 content (Fig. 6b) in sample 14JP2 from the upper part of the Dongrong Formation suggest an advanced stage of weathering. We therefore conclude that the provenance of the Dongrong Formation was initially weak weathering processes. During the deposition of this formation, an increase in the intensity of weathering processes produced compositionally mature alumina-rich sediments.

Sample 14JP3 from the Yanwo Formation yielded PIA, ICV and CIA values of 43.5–84.3, 0.30–1.05 and 44.2–83.6, respectively (Fig. 6; online Supplementary Table S2). These geochemical proxies suggest that the weathering processes, which affected the sources of the Yanwo Formation, were intensified during the deposition of this formation.

Sedimentary recycling in oxidizing conditions commonly results in fractionation of Th and U, because U4+ is readily oxidized to U6+ during weathering (McLennan & Taylor, Reference McLennan and Taylor1980). The removal of relatively soluble U6+ could cause an increase in the Th/U ratio; the Th/U ratio is therefore a good indicator for sedimentary recycling. The Mo/U molar ratio and Mo enrichment factor (MoEF) of sediments are good indicators of the redox conditions of water (Algeo & Tribovillard, Reference Algeo and Tribovillard2009). Under oxic conditions, Mo enrichment is limited and Mo is present as the stable and largely unreactive molybdate oxyanion (MoO42–). Our samples have lower Mo/U molar ratios than suboxic modern seawater (c. 7.5–7.9) with low MoEF (< 3.2), suggesting oxidizing conditions. The Th/U ratios of Unite A2 and Unite B show a positive correlation with the CIA values (Fig. 10; online Supplementary Table S2). Some of the measured Th/U ratios are higher than the average value for PAAS (4.7; McLennan & Taylor, Reference McLennan and Taylor1980), suggesting that parts of the provenance were derived from recycled sediments. These samples are also enriched in Zr (online Supplementary Table S2), which is another indication for the incorporation of recycled weathered material (recycling causes enrichment in heavy minerals, such as zircon and monazite).

Fig. 10. Th/U–CIA diagrams for the Dongrong and Yanwo formations. Symbols as for Figure 5.

5.b.2. Hydrodynamic sorting

The hydrodynamic conditions can be evaluated by textural maturity with respect to grain sizes and shapes, as well as mineralogical and geochemical compositions (McLennan et al. Reference McLennan, Hemming, Mcdaniel, Hanson, Johnsson and Basu1993). The six medium-grained arkose and wacke samples studied by us are compositionally and texturally immature, reflecting a low degree of sedimentary sorting. Geochemically, the samples are characterized by low contents of Al2O3 and high SiO2/Al2O3, suggesting low proportions of aluminous clay minerals. Furthermore, the low contents of HFSE (Zr and Hf), REE and transition metals (Cr and Sc) suggest that heavy minerals are rare or absent. Given that aluminous clay minerals and heavy minerals typically accumulate by sedimentary sorting, our observations support the suggestion that the sediments experienced a low degree of sedimentary sorting, likely associated with a relatively short distance of transport from the area of provenance.

5.b.3. Diagenesis

The geochemistry of sedimentary rocks can be affected by potassium metasomatism and silicification. Potassium metasomatism in sandstone could involve the conversion of aluminous clay minerals (e.g. kaolinite) to illite and/or conversion of plagioclase to K-feldspar (Fedo et al. Reference Fedo, Nesbitt and Young1995). In mudstone, potassium metasomatism is mainly associated with clay alteration. The latter can be identified petrographically by observing illite either in matrix or as an alteration product of partially weathered plagioclase grains. Such features were not observed in our samples. The role of potassium metasomatism can also be evaluated using the A-CN-K plot, with the expectation that samples subjected to this process would be distributed along a path towards the K2O apex of the triangle. This is not the case in our samples, which are depleted in K2O (Fig. 6b). We therefore conclude that our samples did not experience significant potassium metasomatism.

The recognition of polycrystalline quartz in the matrix (Fig. 4c) is indicative of silicification. Further evidence of silicification is metamorphic lithic fragments, which are characterized by polycrystalline quartz (Fig. 3f). The significantly higher SiO2/Al2O3 values in samples 14JP1-4 and 14JP3-6, in comparison with other samples with the same lithology (online Supplementary Table S2), might also be attributed to the process of silicification.

5.c. Interpretation of geochemical proxies

Major element geochemical proxies (e.g. K2O/Al2O3, SiO2/Al2O3 and A-CN-K) are commonly used to obtain information on the composition of source rocks (Nesbitt, Reference Nesbitt1979; Cox et al. Reference Cox, Lowe and Cullers1995; Fedo et al. Reference Fedo, Nesbitt and Young1995; Hayashi et al. Reference Hayashi, Fujisawa, Holland and Ohmoto1997). However, such proxies might be problematic in the Sanjiang Basin because of the mobility of some elements (e.g. Ca2+, Na+, K+ and silica) by weathering, silicification and leaching. As discussed in the previous section, our samples did not experience significant potassium metasomatism, but some samples (14JP2-2, -3, -6 and 14JP3-1, -7, -8) show evidence of recycled sediments. After ruling out these samples, the geochemical proxies indicate that the provenance was dominated by felsic and intermediate rocks, with a minor component of mafic rocks (Fig. 11a). This conclusion is consistent with the petrography and zircon morphology.

Fig. 11. (a) Discriminant function diagram of major-element provenance (Roser & Korsch, Reference Roser and Korsch1986). (b) Eu/Eu* versus Th/Sc (after McLennan et al. Reference McLennan, Taylor, McCulloch and Maynard1990). (c) Co/Th versus La/Sc (after McLennan et al. Reference McLennan, Taylor and Eriksson1983). (d) La-Th-Sc diagram for tectonic setting discrimination (after Bhatia & Crook, Reference Bhatia and Crook1986). (e) Tectonic setting of sandstone derivation from different types of provenances (after Dickinson, Reference Dickinson and Zuffa1985). (f) Cumulative proportions versus time (in Ma) of detrital zircon ages from the Sanjiang Basin (after Cawood et al. Reference Cawood, Hawkesworth and Dhuime2012). Discriminant Function 1 = (−1.773×TiO2) + (0.607×Al2O3)  +  (0.760×Fe2O3)  +  (−1.500×MgO)  +  (0.616×CaO) + (0.509×Na2O) + (−1.224×K2O) + (−9.090); Discriminant Function 2 = (0.445×TiO2) + (0.070×Al2O3) + (−0.250×Fe2O3) + (−1.142×MgO) + (0.438×CaO) + (1.475×Na2O) + (1.462×K2O) + (−6.861). Sample symbols as for Figure 5. The grey symbols correspond to samples from Cao et al. (Reference Cao, Guo, Shan, Ma and Sun2015) and Hu et al. (Reference Hu, Meng and Liu2021).

Some elements, such as LREE, HFSE, Co and Sc, have extremely low water–rock partition coefficients. As such, they are remarkably stable during weathering, transport and diagenesis (Chaudhuri & Cullers, Reference Chaudhuri and Cullers1979). These elements can be used as a provenance fingerprint. Diagrams of Eu/Eu* versus Th/Sc and Co/Th versus La/Sc show that the provenance areas were associated mainly with felsic rocks, without any mafic contribution (Fig. 11b, c). These results are inconsistent with our previous conclusion that a minor portion of the provenance was made of intermediate and mafic rocks. This discrepancy could be explained by the relatively weak chemical signal of mafic provenance, which means that a slight contribution from a mafic source might not be detectable in the whole-rock geochemistry of the sedimentary rock. Accordingly, the suggestion that the provenance was dominated by felsic rocks with a minor amount of intermediate and mafic rocks is not unreasonable.

The original tectonic setting of the provenance region could be inferred from the La-Th-Sc triangle plot (Fig. 11d; McLennan & Taylor, Reference McLennan and Taylor1991). According to this interpretation, samples 14JP1 and 14JP2 (from the Dongrong Formation) were derived from a continental margin (active or passive) and a continental island arc, respectively. The sandstone compositions indicate that the provenance changed from a transitional continental setting to recycled orogen (Fig. 11e). Indeed, the onset of the Paleo-Pacific plate subduction beneath the continental margin of East Asia during latest Triassic time likely resulted in a change in the source regions (Ma et al. Reference Ma, Zhu, Zhou and Qiao2017). In addition, the overall spread and proportion of detrital zircon ages indicate that the Dongrong Formation formed in a convergent plate margin (Fig. 11f; Cawood et al. Reference Cawood, Hawkesworth and Dhuime2012), which is also consistent with the tectonic evolution of Jiamusi Massif during the Triassic Period.

As demonstrated by the La-Th-Sc triangle plot (Fig. 11d), the provenance of the Yanwo Formation formed in a continental margin. The sandstone compositions indicate that the deposits were derived mainly from a recycled orogen (Fig. 11e). In contrast, detrital zircon data indicate that the Yanwo Formation formed in a collisional plate margin (Fig. 11f). The tectonic evolution of the study area was controlled by W-dipping subduction (present coordinate) since Early Jurassic time, which led to the development of a fold–thrust belt and, ultimately, to the collision between the Songnen and Jiamusi massifs. Accordingly, the formation of the Yanwo Formation could be linked to the subduction of the Paleo-Pacific Ocean and the associated accretionary orogenesis (Zhou & Li, Reference Zhou and Li2017).

5.d. Sedimentary provenance

Our sandstone petrography results are indicative of low degrees of sedimentary sorting, which likely reflect a short distance of sedimentary transport. Potential source regions could therefore be the Jiamusi Massif, Lesser Xing’an and Zhangguangcai ranges, or the Nadanhada Terrane (Fig. 1b). The geochemistry of the sedimentary sequence indicates that parts of the deposits included recycled sediments, and that first-cycle sediments were derived mainly from felsic rocks with a minor contribution from intermediate and mafic rocks. The combination of U–Pb ages and Hf data allows us to evaluate the relative contributions from the different sources.

5.d.1. Provenance of the Dongrong Formation

Zircon U–Pb dating and Hf isotope data indicate that the provenance of samples 14JP1 and 14JP2 included Cambrian, Upper Ordovician, Permian and Triassic igneous rocks, with a minor component of Neoproterozoic, Devonian and Silurian rocks (Fig. 12a). Neoproterozoic and Cambrian detrital zircon grains were likely derived from the interior of the Jiamusi Massif. The minor amount of Upper Cambrian detrital zircon grains could have been derived from the Lesser Xing’an Range (Fig. 12a). Specifically, granitic gneisses and granites dated at c. 895, 751, 541 and 507–498 Ma from the Jiamusi Massif and Lesser Xing’an Range could represent the provenance terranes (Wang et al. Reference Wang, Xu, Pei, Wang and Cao2015; Yang et al. Reference Yang, Ge, Zhao, Bi, Wang, Dong and Xu2017, Reference Yang, Ge, Bi, Wang, Tian, Dong and Chen2018; Luan et al. Reference Luan, Yu, Yu, Cui and Xu2019). We note that the c. 551 Ma detrital zircon grains have similar ϵHf(t) values as inherited zircon grains from granites in the Yuejinshan Complex, which is part of the Nadanhada Terrane (Bi et al. Reference Bi, Ge, Yang, Wang, Dong, Liu and Ji2017).

The source region of the Ordovician (c. 477 and 458–453 Ma) detrital zircon grains could be associated with the Lesser Xing’an and Zhangguangcai ranges (Fig. 12a). Granitoids from the Lesser Xing’an Range, dated c. 476 Ma and c. 457 Ma, likely provided detritus to the Sanjiang Basin (Wang et al. Reference Wang, Xu, Pei, Wang and Guo2016 b). Silurian–Carboniferous detrital zircon grains were also likely derived from the Lesser Xing’an and Zhangguangcai ranges, with a minor contribution from the Middle Jilin Orogen (Fig. 12a). Indeed, zircon grains found in a Devonian (366 ± 5 Ma) granite from the Lesser Xing’an range are characterized by similar ϵHf(t) values as those found in detrital zircon grains from the Dongrong Formation (Guo et al. Reference Guo, Xu, Wang, Wang and Luan2018). In addition, 364–340 Ma rhyolites and granites from the Middle Jilin Orogen might have contributed to the detrital sediments (Wang et al. Reference Wang, Xu, Pei, Wang and Cao2015).

The main source rocks of the Dongrong Formation were Permian (291–257 Ma) and Early Triassic (c. 240 Ma) in age, which could be associated with four potential source regions (Fig. 12a): (1) 300–261 Ma granitoids from the Jiamusi Massif (Yu et al. Reference Yu, Wang, Xu, Gao and Tang2013; Yang et al. Reference Yang, Ge, Zhao, Yu and Zhang2015 b; Bi et al. Reference Bi, Ge, Yang, Wang, Xu, Yang, Xing and Chen2016; Dong et al. Reference Dong, Ge, Yang, Bi, Wang and Xu2017 a, b; Long et al. Reference Long, Xu, Guo, Sun and Luan2020); (2) 293–257 Ma rhyolites, andesites and gabbros, and 245 Ma granitoids from the Lesser Xing’an and Zhangguangcai ranges (Long et al. Reference Long, Xu, Guo, Sun and Luan2019, Reference Long, Xu, Guo, Sun and Luan2020); (3) 263–256 Ma ganitoids, monzonites and gabbros, and c. 244 Ma monzonites and granites from the Middle Jilin Orogen (Cao et al. Reference Cao, Xu, Pei, Wang, Wang and Wang2013; Wang et al. Reference Wang, Xu, Pei, Wang and Cao2015); and (4) 288–269 Ma granitoids and gabbros from the Nadanhada Terrane (Bi et al. Reference Bi, Ge, Yang, Zhao, Yu, Zhang, Wang and Tian2014 a, Reference Bi, Ge, Yang, Zhao, Xu and Wang2015; Sun et al. Reference Sun, Xu, Wilde and Chen2015). Hf compositions of Upper Triassic detrital zircon grains are similar to those found in c. 210 Ma granites and gabbros from the Lesser Xing’an and Zhangguangcai ranges, indicating that the sediments were likely derived from these granites (Long et al. Reference Long, Xu, Guo, Sun and Luan2020).

5.d.2. Provenance of the Yanwo Formation

The provenance of the Yanwo Formation is associated mainly with upper Permian, Lower Triassic, Middle Jurassic and Upper Cretaceous rocks, with a minor contribution of Neoproterozoic, Ordovician and lower Permian rocks (Fig. 12b). Neoproterozoic detrital zircon grains are characterized by two groups with distinct Hf compositions. The first group of detrital zircon grains (837–831 Ma) were possibly derived from c. 840 Ma gneissic granodiorites from the Lesser Xing’an Range (Luan et al. Reference Luan, Xu, Wang, Wang and Guo2017 b). The origin of the second group (c. 841 Ma zircon grains with ϵHf(t) values of 10.7 ± 0.9) remains unclear. Cambrian and Lower Ordovician detrital zircon grains were possibly derived from granitoids in the Lesser Xing’an Range and the Jiamusi Massif (Fig. 12b). The 444 ± 5 Ma detrital zircons have similar Hf compositions as those from the 450 ± 2 Ma granites in the Lesser Xing’an Range (Wang et al. Reference Wang, Xu, Pei, Wang and Cao2015). The 367 ± 11 Ma zircons have similar Hf composition as inherited zircons from the 267 ± 5 Ma rhyolite in the Yuejinshang Complex (Bi et al. Reference Bi, Ge, Yang, Wang, Dong, Liu and Ji2017). Consequently, these granites and rhyolites likely contributed to the sediments of the Yanwo Formation.

Permian and Triassic zircon grains from the Yanwo Formation have similar Hf compositions to those from the Dongrong Formation (Fig. 12b), suggesting that their provenances were similar. The Hf compositions of Jurassic detrital zircons indicate that their provenance was associated with the Middle Jilin Orogen, Nadanhada Terrane and the Zhangguangcai Range. The source rocks might include 175–172 Ma and 164 Ma granitoids and diorites from the Middle Jilin Orogen, c. 174 and 167 Ma rhyolites, andesites and basalts from the Nadanhada Terrane (Wang et al. Reference Wang, Xu, Pei, Wang and Cao2015, Reference Wang, Xu, Xu, Yang, Wu and Sun2017), and c. 177 Ma granitoids from the Zhangguangcai Range (Long et al. Reference Long, Xu, Guo, Sun and Luan2020). The Hf compositions of Cretaceous detrital zircon grains are similar to these from Nadahada Terrane and the Jiamusi Massif (Fig. 12b). Granitoids from the interior of the Jiamusi Massif, dated at c. 117 and 95–89 Ma (Yu et al. Reference Yu, Wang, Xu, Gao and Tang2013), and granitoids from Nadanhada Terrane (c. 133 and 120 Ma), might have contributed to the sediments of the Yanwo Formation (Wang et al. Reference Wang, Xu, Xu, Yang, Wu and Sun2017).

In summary, the evidence suggests that the sediments of the Sanjiang Basin were derived mainly from the Jiamusi Massif and the surrounding orogens and terranes. Source rocks were dominated by felsic rocks with a minor amount of intermediate and mafic rocks.

5.e. Mesozoic uplift history of the northeastern Asian continental margin

During the Mesozoic Era, the northeastern Asian continental margin was an active continental margin, characterized by a fore-arc region, island arc and retro-arc foreland basin systems (Zhang et al. Reference Zhang, Yan and Ji2019). Subduction processes were likely responsible for uplift and exhumation of the continental margin, thus controlling drainage systems and the architecture of sedimentary basins (Zhang et al. Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015; Zhou et al. Reference Zhou, István, Liu, Li and Eynatten2020). Many previous studies have investigated surface processes of the northeastern Asian continental margin using information associated with the distribution of magmatism (e.g. Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Cao et al. Reference Cao, Xu, Pei, Wang, Wang and Wang2013; Song et al. Reference Song, Stepashko and Ren2015; Wang et al. Reference Wang, Xu, Xing, Wang, Zhang, Wu, Sun and Ge2019), contractional structures (e.g. Zhou & Li, Reference Zhou and Li2017; Zheng et al. Reference Zheng, Sun, Wan, Wang, He and Zhang2019) and unconformities within the stratigraphic succession (e.g. Deng et al. Reference Deng, He, Pan and Zhu2013; Song et al. Reference Song, Ren, Stepashko and Li2014). The integrated evidence from the Sanjiang Basin, based on stratigraphic data, provenance studies and seismic reflection profiles (Zhang et al. Reference Zhang, Dilek, Chen, Yang and Meng2017), provides new insights into the uplift history and associated subduction processes.

Crustal-scale shortening and uplift are generally expressed by the presence of changes in provenance, unconformities, positive flower structures, imbricate thrusts and folds (Fig. 2). The presence of a Jurassic unconformity indicates that the study area might have experienced one stage of surface uplift. The Jurassic provenance changes of basinal strata from Lesser Xing’an and Zhangguangcai ranges to Nadanhada Terrane and Middle Jilin Orogen, and the Jurassic detrital zircons of Yanwo Formation, were mainly derived from the Nadanhada Terrane in the eastern of the Sanjiang Basin. The provenance change appears to correspond to surface uplift in the eastern part of the Sanjiang Basin. Based on geochronological and geochemical studies, it has been suggested that the tectonic evolution of the northeastern Asian continental margin was controlled by W-wards subduction (present coordinate) of the Paleo-Pacific Ocean during the Early Jurassic Epoch (190–159 Ma; Wu et al. Reference Wu, Sun, Ge, Zhang, Grant, Wilde and Jahn2011; Xu et al. Reference Xu, Pei, Wang, Meng, Ji, Yang and Wang2013). We propose that the uplift of the continental margin resulted from the initiation of this W-wards subduction.

Positive flower structures and imbricate thrusts are recognized in Lower Cretaceous strata (cross-section A–B; Fig. 2), and folds and imbricate thrusts occur in Upper Cretaceous strata (cross-section C-D; Fig. 2). These contractional structures are consistent with the two major unconformities within the Cretaceous strata (Fig. 2), suggesting that the study area experienced at least two stages of contractional deformation and uplift during the Cretaceous Period. The Cretaceous unconformity in the Sanjiang Basin seems to be expressed by the absence of detrital zircon ages dated 117–144, 133146 and 115145 Ma in the Lower Cretaceous Didao, Chengzihe and Muling formations, respectively (Fig. 13). In comparison, the Upper Cretaceous Hailang and Yanwo formations lack zircon grains dated at 97144, 99120 and 131149 Ma (Fig. 13), respectively. According to the geochronological compatibility and the two major unconformities, we conclude that the northeastern Asian continental margin experienced at least two stages of uplift during the Cretaceous Period (i.e. 144–133 and 120–99 Ma).

Fig. 13. Comparison of probability diagrams of detrital zircon ages for sedimentary sequences from the Sanjiang Basin. Reference data for the Didao, Chengzihe, Muling and Hailang formations from Zhang et al. (Reference Zhang, Chen, Batt, Dilek, Na, Sun, Yang, Meng and Zhao2015).

During latest Jurassic – Early Cretaceous time, the Izanagi Plate (Paleo-Pacific Ocean) underwent a 24° clockwise rotation (Seton et al. Reference Seton, Müller, Zahirovic, Gaina, Torsvik, Shephard, Talsma, Gurnis, Turner, Maus and Chandler2012). This change in the plate motion could have generated shear stresses along the continental margin, which might have led to the uplift during 144–133 Ma (Fig. 14b). Subsequently, during the latest Early Cretaceous Epoch, the subduction dip angle of the Paleo-Pacific slab was reduced (Sun et al. Reference Sun, Chen, Milan, Wilde, Jourdan and Xu2018), enhancing compressional deformation and further surface uplift (Fig. 14c). The two stages of uplift were likely responsible for the formation of the unconformity in the Cretaceous strata and the development of contractional structures (positive flower structures and imbricate thrusts). In conclusion, the northeastern Asian continental margin experienced at least three stages of uplift, possibly driven by subduction initiation (Early Jurassic) and the changes in the plate kinematics (Early and Late Cretaceous).

Fig. 14. Schematic illustration of the Early Jurassic and Cretaceous uplift events in the northeastern Asian continental margin.

6. Conclusions

Based on sandstone petrography, geochemistry of sedimentary rocks, U–Pb ages and Hf isotope of detrital zircons, the following conclusions are drawn.

  1. 1. The Dongrong Formation was deposited after 208–201 Ma. Deposition of the Yanwo Formation occurred after 85 Ma and before 54 Ma.

  2. 2. The two basinal strata were subjected to weathering processes that progressively intensified during the deposition of these formations. The rocks experienced low degrees of sedimentary sorting, silicification and insignificant potassium metasomatism. Some of the mudstones contain recycled sediments.

  3. 3. The provenance of the Dongrong and Yanwo formations was likely associated with the Jiamusi Massif and surrounding orogens and terranes. The source rocks were dominated by felsic rocks with a minor component of mafic rocks.

  4. 4. The northeastern Asian continental margin experienced at least three stages of uplift. The first uplift phase was possibly driven by the initiation of subduction of the Paleo-Pacific Ocean during Early Jurassic time, whereas the second and third phases are interpreted to be associated with a change in plate kinematics during Early Cretaceous time.

Acknowledgments

We thank the editor Stephen Hubbard and two anonymous reviewers for their constructive suggestions and comments. We are grateful to staff of the Wuhan Samplesolution Analytical Technology Co., Ltd. and the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, for their advice and assistance during zircon LA-ICP-MS U–Pb dating, major- and trace-element analyses, and Hf isotope analyses. This work was financially supported by the Fundamental Research Funds for the Central Universities (N2123029) and the Natural Science Foundation of Hebei Province (D2020501002). All data are archived at Mendeley Date (https://data.mendeley.com/datasets/nfy7gxg38p/1).

Conflict of interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Supplementary material

For supplementary material accompanying this paper visit https://doi.org/10.1017/S0016756821000996

References

Algeo, TJ and Tribovillard, N (2009) Environmental analysis of paleoceanographic systems based on molybdenum–uranium covariation. Chemical Geology 268, 211–25.CrossRefGoogle Scholar
Anderson, T (2002) Correction of common Lead in U–Pb analyses that do not report 204Pb. Chemical Geology 192, 5979.CrossRefGoogle Scholar
Armstrong-Altrin, JS, Lee, YI, Verma, SP and Ramasamy, S (2004) Geochemistry of sandstones from the upper miocene Kudankulam Formation, southern India: implications for provenance, weathering, and tectonic setting. Journal of Sedimentary Research 74, 285–97.CrossRefGoogle Scholar
Ash, SR (1993) Two new late Triassic plants from the Petrified forest of Arizona. Journal of Paleontology 47, 4653.Google Scholar
Barth, AP, Wooden, JL, Jacobson, CE and Probst, K (2004) U–Pb geochronology and geochemistry of the McCoy Mountains Formation, southeastern California: a Cretaceous retroarc foreland basin. Geological Society of America Bulletin 116, 142–53.CrossRefGoogle Scholar
Belousova, EA, Griffin, WL, O’Reilly, SY and Fisher, NI (2002) Igneous zircon: trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology 143, 602622.CrossRefGoogle Scholar
Bhatia, MR and Crook, KAW (1986) Trace-element characteristics of graywackes and tectonic setting discrimination of sedimentary basins. Contributions to Mineralogy and Petrology 92, 181–93.CrossRefGoogle Scholar
Bi, JH, Ge, WC, Yang, H, Wang, ZH, Dong, Y, Liu, XW and Ji, Z (2017) Age, petrogenesis, and tectonic setting of the Permian bimodal volcanic rocks in the eastern Jiamusi Massif, NE China. Journal of Asian Earth Sciences 134, 160–75.CrossRefGoogle Scholar
Bi, JH, Ge, WC, Yang, H, Wang, ZH, Xu, WL, Yang, JH, Xing, DH and Chen, HJ (2016) Geochronology and geochemistry of late Carboniferous-middle Permian I- and A-type granites and gabbro-diorites in the eastern Jiamusi Massif, NE China: implications for petrogenesis and tectonic setting. Lithos 266–267, 213–32.CrossRefGoogle Scholar
Bi, JH, Ge, WC, Yang, H, Zhao, GC, Xu, WL and Wang, ZH (2015) Geochronology, geochemistry and zircon Hf isotopes of the Dongfanghong gabbroic complex at the eastern margin of the Jiamusi Massif, NE China: petrogenesis and tectonic implications. Lithos 234–235, 2746.CrossRefGoogle Scholar
Bi, JH, Ge, WC, Yang, H, Zhao, GC, Yu, JJ, Zhang, YL, Wang, ZH and Tian, DX (2014a) Petrogenesis and tectonic implications of early Paleozoic granitic magmatism in the Jiamusi Massif, NE China: geochronological, geochemical and Hf isotopic evidence. Journal of Asian Earth Sciences 96, 308–31.CrossRefGoogle Scholar
Bi, JH, Ge, WC, Zhang, YL, Yang, H and Wang, ZH (2014b) Petrogenesis of Permian Jinshan granitic complex in the eastern Jiamusi Massif and its geological implications. Journal of Earth Sciences and Environment 36, 1631 (in Chinese with English abstract).Google Scholar
Boynton, WV (1984) Geochemistry of the rare earth elements: meteorite studies. In Rare Earth Element Geochemistry (ed Henderson, P), pp. 63114. Amsterdam: Elsevier.CrossRefGoogle Scholar
Bruguier, O and Lancelet, JR (1997) U–Pb dating on single detrital zircon grains from the Triassic Songpan-Ganze flysch (central China): provenance and tectonic correlations. Earth and Planetary Science Letters 152, 217231.CrossRefGoogle Scholar
Campbell, MJ, Rosenbaum, G, Allen, CM and Mortimer, N (2020) Origin of dispersed Permian–Triassic fore-arc basin terranes in New Zealand: Insights from zircon petrochronology. Gondwana Research 78, 210–27.CrossRefGoogle Scholar
Cao, HH, Xu, WL, Pei, FP, Wang, ZW, Wang, F and Wang, ZJ (2013) Zircon U–Pb geochronology and petrogenesis of late Paleozoic-Early Mesozoic intrusive rocks in the eastern segment of the northern margin of the North China Block. Lithos 170–171, 191207.CrossRefGoogle Scholar
Cao, HS, Guo, W, Shan, XL, Ma, L and Sun, PC (2015) Paleolimnological environments and organic accumulation of the Nenjiang Formation in the southern Songliao Basin, China. Oil Shale 32, 524.CrossRefGoogle Scholar
Cawood, PA, Hawkesworth, CJ and Dhuime, B (2012) Detrital zircon record and tectonic setting. Geology 40, 875–78.CrossRefGoogle Scholar
Chaudhuri, S and Cullers, RL (1979) The distribution of rare-earth elements in deeply buried Gulf Coast sediments. Chemical Geology 24, 327–38.CrossRefGoogle Scholar
Cherniak, DJ and Watson, EB (2000) Pb diffusion in zircon. Chemical Geology 172, 524.CrossRefGoogle Scholar
Cox, R, Lowe, DR and Cullers, RL (1995) The influence of sediment recycling and basement composition on evolution of mudrock chemistry in the southwestern United States. Geochimica et Cosmochimica Acta 59, 2919–40.CrossRefGoogle Scholar
Cullers, RL (1994) The chemical signature of source rocks in size fractions of Holocene stream sediment derived from metamorphic rocks in the Wet Mountains region, Colorado, USA. Chemical Geology 113, 327–43.CrossRefGoogle Scholar
Deng, CL, He, HY, Pan, YX and Zhu, RX (2013) Chronology of the terrestrial Upper Cretaceous in the Songliao Basin, northeast Asia. Palaeogeography, Palaeoclimatology, Palaeoecology 385, 4454.CrossRefGoogle Scholar
Dickinson, WR (1985) Interpreting provenance relations from detrital modes of sandstones. In Provenance of Arenites (ed Zuffa, GG), pp. 333–61. Boston: D. Reidel Publishing Company.CrossRefGoogle Scholar
Dickinson, WR and Suczek, CA (1979) Plate tectonics and sandstone compositions. AAPG Bulletin 63, 2167–82.Google Scholar
Dong, Y, Ge, WC, Yang, H, Bi, JH, Wang, ZH and Xu, WL (2017a) Permian tectonic evolution of the Mudanjiang Ocean: evidence from zircon U–Pb–Hf isotopes and geochemistry of a N-S trending granitoid belt in the Jiamusi Massif, NE China. Gondwana Research 49, 147–63.CrossRefGoogle Scholar
Dong, Y, Ge, WC, Yang, H, Xu, WL, Bi, JH and Wang, ZH (2017b) Geochemistry and geochronology of the Late Permian mafic intrusions along the boundary area of Jiamusi and Songnen-Zhangguangcai Range massifs and adjacent regions, northeastern China: petrogenesis and implications for the tectonic evolution of the Mudanjiang Ocean. Tectonophysics 694, 356–67.CrossRefGoogle Scholar
Falcon-Lang, HJ (2002) Terrestrial paleoecology of the Cretaceous (Early Aptian) Cerro Negro Formation, South Shetland Islands, Antarctica: a record of polar vegetation in a volcanic arc environment. Palaios 17, 491506.2.0.CO;2>CrossRefGoogle Scholar
Fedo, CM, Eriksson, KA and Krogstad, EJ (1996) Geochemistry of shales from the Archean (∼3.0 Ga) Buhwa Greenstone Belt, Zimbabwe: implications for provenance and source-area weathering. Geochimica et Cosmochimica Acta 60, 1751–63.CrossRefGoogle Scholar
Fedo, CM, Nesbitt, HW and Young, GM (1995) Unraveling the effects of potassium metasomatism in sedimentary rocks and paleosols, with implications for paleoweathering conditions and provenance. Geology 23, 921–24.2.3.CO;2>CrossRefGoogle Scholar
Guo, P, Xu, WL, Wang, ZW, Wang, F and Luan, JP (2018) Geochronology and geochemistry of Late Devonian–Carboniferous igneous rocks in the Songnen-Zhangguangcai Range Massif, NE China: Constraints on the late Paleozoic tectonic evolution of the easten Central Asian Orogenic Belt. Gondwama Research 57, 119–32.CrossRefGoogle Scholar
Guo, P, Xu, WL, Yu, JJ, Wang, F, Tang, J and Li, Y (2016) Geochronology and geochemistry of Late Triassic bimodal igeous rocks at the eastern margin of Songnen–Zhangguangcai Range Massif, NE China: petrogenesis and tectonic implications. International Geology Review 58, 196215.CrossRefGoogle Scholar
Han, J, Zhou, JB, Li, L and Song, MC (2017) Mesoproterozoic (∼1.4 Ga) A-type gneissic granites in the break-up of Columbia in the Eastern CAOB. Precambrian Research 296, 2038.CrossRefGoogle Scholar
Hawkesworth, CJ and Kemp, AIS (2006) Evolution of the continental crust. Nature 443, 811–17.CrossRefGoogle ScholarPubMed
Hayashi, KI, Fujisawa, H, Holland, HD and Ohmoto, H (1997) Geochemistry of ∼1.9 Ga sedimentary rocks from northeastern Labrador, Canada. Geochimica et Cosmochimica Acta 61, 4115–37.CrossRefGoogle ScholarPubMed
HBGMR (Heilongjiang Bureau of Geology and Mineral Resources) (1993) Regional Geology of Heilongjiang Province. Beijing: Geological Publishing House, pp. 8584 (in Chinese with English abstract).Google Scholar
Herron, MM (1988) Geochemical classification of terrigenous sands and shales from core or log data. Journal of Sedimentary Petrology 58, 820–29.Google Scholar
Hu, F, Meng, QT and Liu, ZJ (2021) Mineralogy and element geochemistry of oil shales in the Lower Cretaceous Qingshankou Formation of the southern Songliao Basin, northeast China: implications of provenance, tectonic setting, and paleoenvironment. ACS Earth Space Chemistry 5, 365380.CrossRefGoogle Scholar
Hu, ZC, Liu, YS, Gao, S, Hu, SH, Dietiker, R and Günther, D (2008) A local aerosol extraction strategy for the determination of the aerosol composition in laser ablation inductively coupled plasma mass spectrometry. Journal of Analytical Atomic Spectrometry 23, 1192–203.CrossRefGoogle Scholar
Hu, ZC, Liu, YS, Gao, S, Liu, WG, Yang, L, Zhang, W, Tong, XR, Lin, L, Zong, KQ, Li, M, Chen, HH and Zhou, L (2012) Improved in situ Hf isotope ratio analysis of zircon using newly designed X skimmer cone and Jet sample cone in combination with the addition of nitrogen by laser ablation multiple collector ICP-MS. Journal of Analytical Atomic Spectrometry 27, 1391–99.CrossRefGoogle Scholar
Ingersoll, RV, Fullard, TF, Ford, RL, Grimm, JP, Pickle, JD and Sares, SW (1984) The effect of grain size on detrital modes: a test of the Gazzi-Dickinson point-counting method. Journal of Sediment Petrology 54, 103–16.Google Scholar
Kosler, J and Sylvester, PJ (2003) Present trends and the future of zircon in geochronology: laser ablation ICPMS Zircon. Reviews in Mineralogy and Geochemistry 53, 243–71.CrossRefGoogle Scholar
Kröner, A, Kovach, V, Belousova, E, Hegner, E, Armstrong, R, Dolgopolova, A, Seltmann, R, Alexeiev, DV, Hoffmann, JE, Wong, J, Sun, M, Cai, K, Wang, T, Tong, Y, Wilde, SA, Degtyarev, KE and Rytsk, E (2014) Reassessment of continental growth during the accretionary history of the Central Asian Orogenic Belt. Gondwana Research 25, 103–25.CrossRefGoogle Scholar
Lawton, TF, Bradford, IA, Vega, FJ, Gehrels, GE and Amato, JM (2009) Provenance of upper cretaceous-Paleogene sandstones in the foreland basin system of the Sierra Madre Oriental, northeastern Mexico, and its bearing on fluvial dispersal systems of the Mexican Laramide Province. Geological Society of America Bulletin 121, 820–36.CrossRefGoogle Scholar
Li, PC, Liu, ZH, Li, SC, Zhao, QY, Shi, Q, Li, CH and Yang, XH (2019) Late Paleocene-early Eocene granitoids in the Jiamusi Massif, NE China: Zircon U–Pb ages, geochemistry and tectonic implications. International Geology Review 61, 116.CrossRefGoogle Scholar
Liu, YS, Gao, S, Hu, ZC, Gao, CG, Zong, KQ and Wang, DB (2010) Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U–Pb dating, Hf isotopes and trace elements in zircons of mantle xenoliths. Journal of Petrology 51, 537–71.CrossRefGoogle Scholar
Long, XP, Sun, M, Yuan, C, Xiao, WJ and Cai, KD (2008) Early Paleozoic sedimentary record of the Chinese Altai: implication for its tectonic evolution. Sedimentary Geology 208, 88100.CrossRefGoogle Scholar
Long, XY, Xu, WL, Guo, P, Sun, CY and Luan, JP (2019) Was Permian magmatism in the eastern Songnen and western Jiamusi massifs, NE China, related to the subduction of the Mudanjiang oceanic plate? Geological Journal 55, 1781–807.CrossRefGoogle Scholar
Long, XY, Xu, WL, Guo, P, Sun, CY and Luan, JP (2020) Opening and closure history of the Mudanjiang Ocean in the eastern Central Asian Orogenic Belt: geochronological and Geochemical constraints from early Mesozoic intrusive rocks. Gondwana Research 84, 111–30.CrossRefGoogle Scholar
Luan, JP, Wang, F, Xu, WL, Ge, WC, Sorokin, AA, Wang, ZW and Guo, P (2017a) Provenance, age, and tectonic implications of Neoproterozoic strata in the Jiamusi Massif: evidence from U–Pb ages and Hf isotope compositions of detrital and magmatic zircons. Precambrian Research 297, 1732.CrossRefGoogle Scholar
Luan, JP, Xu, WL, Wang, F, Wang, ZW and Guo, P (2017b) Age and geochemistry of Neoproterozoic granitoids in the Songnen–Zhangguangcai Range Massif, NE China: Petrogenesis and tectonic implications. Journal of Asian Earth Sciences 148, 265–76.CrossRefGoogle Scholar
Luan, JP, Yu, JJ, Yu, JL, Cui, YC and Xu, WL (2019) Early Neoproterozoic magmatism and associated metamorphism in the Songnen Massif, NE China: Petrogenesis and tectonic implications. Precambrian Research 328, 250–68.CrossRefGoogle Scholar
Ludwig, KR (2003) ISOPLOT 3: A Geochronological Toolkit for Microsoft Excel. Berkeley: Berkeley Geochronology Centre, Special Publication no. 4, 74 pp.Google Scholar
Ma, XH, Zhu, WP, Zhou, ZH and Qiao, SL (2017) Transformation from Paleo-Asian Ocean clousure to Paleo-Pacific subduction: New constraints from granitoids in the eastern Jilin-Heilongjiang Belt, NE China. Journal of Asian Earth Sciences 144, 261–86.CrossRefGoogle Scholar
McLennan, SM, Hemming, S, Mcdaniel, DK and Hanson, GN (1993) Geochemical approaches to sedimentation, provenance, and tectonics. In Processes Controlling the Composition of Clastic Sediments (eds Johnsson, MJ and Basu, A), pp. 240. Geological Society of America, Boulder, Special Paper, 285.Google Scholar
McLennan, SM and Taylor, SR (1980) Th and U in sedimentary rocks: crustal evolution and sedimentary recycling. Nature 285, 621–24.CrossRefGoogle Scholar
McLennan, SM and Taylor, SR (1991) Sedimentary rocks and crustal evolution: tectonic setting and secular trends. Journal of Geology 99, 121.CrossRefGoogle Scholar
McLennan, SM, Taylor, SR and Eriksson, KA (1983) Geochemistry of Archean shales from the Pilbara Supergroup, Western Australia. Geochimica et Cosmochimica Acta 47, 1211–22.CrossRefGoogle Scholar
McLennan, SM, Taylor, SR, McCulloch, MT and Maynard, JB (1990) Geochemical and Nd–Sr isotopic composition of deep-sea turbidites: crustal evolution and plate tectonic associations. Geochimica et Cosmochimica Acta 54, 2015–50.CrossRefGoogle Scholar
Migani, F, Borghesi, F and Dinelli, E (2015) Geochemical characterization of surface sediments from the northern Adriatic wetlands around the Poriver delta. Part I: Bulk composition and relation to local background. Journal of Geochemical Exploration 156, 7288.CrossRefGoogle Scholar
Nesbitt, HW (1979) Mobility and fractionation of rare earth elements during weathering of a granodiorite. Nature 279, 206–10.CrossRefGoogle Scholar
Nesbitt, HW and Young, GM (1982) Early Proterozoic climates and plate motions inferred from major element chemistry of lutites. Nature 299, 715–7.CrossRefGoogle Scholar
Nesbitt, HW and Young, GM (1989) Formation and diagenesis of weathering profiles. Journal of Geology 97, 129–47.CrossRefGoogle Scholar
Pidgeon, RT, Nemchin, AA and Hitchen, GJ (1998) Internal structures of zircons from Archaean granites from the Darling Range batholith: implications for zircon stability and the interpretation of zircon U–Pb ages. Contributions to Mineralogy and Petrology 132, 288–99.CrossRefGoogle Scholar
Rainbird, RH, Hamilton, MA and Young, GM (2001) Detrital zircon geochronology and provenance of the Torridonian, NW Scotland. Journal of the Geological Society of London 158, 1527.CrossRefGoogle Scholar
Reiners, PW, Campbell, IH, Nicolescu, S, Allen, CM, Hourigan, JK, Garver, JI, Mattinson, JM and Cowan, DS (2005) (U–Th)/(He–Pb) double dating of detrital zircons. American Journal of Science 305, 259311.CrossRefGoogle Scholar
Roser, BP and Korsch, RJ (1986) Determination of tectonic setting of sandstone–mudstone suites using SiO2 content and K2O/Na2O ratio. Journal of Geology 94, 635–50.CrossRefGoogle Scholar
Rubatto, D (2002) Zircon trace element geochemistry: partitioning with garnet and the link between U–Pb ages and metamorphism. Chemical Geology 184, 123–38.CrossRefGoogle Scholar
Rudnick, RL and Gao, S (2003) Composition of the continental crust. In Treatise on Geochemistry, volume 3 (eds Holland, HD and Turekian, KK), pp. 164. Amsterdam: Elsevier.Google Scholar
Rudnick, RL, Gao, S, Ling, WL, Liu, YS and Mc Donough, WF (2004) Petrology and geochemistry of spinel peridotite xenoliths from Hannuoba and Qixia, North China Craton. Lithos 77, 609–37.CrossRefGoogle Scholar
Sengör, AMC, Natal’in, BA and Burtman, VS (1993) Evolution of the Altaid tectonic collage and Paleozoic crustal growth in Eurasia. Nature 364, 299307.CrossRefGoogle Scholar
Seton, M, Müller, RD, Zahirovic, S, Gaina, C, Torsvik, T, Shephard, G, Talsma, A, Gurnis, M, Turner, M, Maus, S and Chandler, M (2012) Global continental and ocean basin reconstructions since 200 Ma. Earth-Science Reviews 113, 212–70.CrossRefGoogle Scholar
Shaanan, U and Rosenbaum, G (2018) Detrital zircons as palaeodrainage indicators: insights into southeastern Gondwana from Permian basins in eastern Australia. Basin Research 30, 3647.CrossRefGoogle Scholar
Shaanan, U, Rosenbaum, G and Campbell, MJ (2019) Detrital fingerprint: the use of early Precambrian zircon age spectra as unique identifiers of Phanerozoic terranes. Earth and Planetary Science Letters 506, 97103.CrossRefGoogle Scholar
Song, Y, Ren, JY, Stepashko, A and Li, JG (2014) Post-rift geodynamics of the Songliao Basin, NE China: Origin and significance of T11 (coniacian) unconformity. Tectonophysics 634, 118.CrossRefGoogle Scholar
Song, Y, Stepashko, A and Ren, JY (2015) The Cretaceous climax of compression in Eastern Asia: age of 87–89 Ma (late Turonian/Coniacian), Pacific cause, continental consequences. Cretaceous Research 55, 262–84.CrossRefGoogle Scholar
Sun, M, Chen, H, Milan, LA, Wilde, SA, Jourdan, F and Xu, Y (2018) Continental arc and back-arc migration in eastern NE China: New constraints on Cretaceous Paleo-Pacific subduction and rollback. Tectonics 37, 3893–915.CrossRefGoogle Scholar
Sun, MD, Xu, YG, Wilde, SA and Chen, HL (2015) Provenance of Cretaceous trench slope sediments from the Mesozoic Wandashan Orogen, NE China: Implications for determining ancient drainage systems and tectonics of the Paleo-Pacific. Tectonics 34, 1269–89.CrossRefGoogle Scholar
Sun, SS and McDonough, WF (1989) Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In Magmatism in Ocean Basins (eds Saunders, AD and Norry, MJ), pp. 313345. Geological Society of London, Special Publication no. 42.Google Scholar
Tao, K, Niu, FL, Ning, JY, Chen, YSJ, Grand, S, Kawakatsu, H, Tanaka, S, Obayashi, M and Ni, J (2014) Crustal structure beneath NE China imaged by NECESSArray receiver function data. Earth and Planetary Science Letters 398, 4857.CrossRefGoogle Scholar
Vavrek, MJ and Larsson, HCE (2007) A Late Triassic flora from wast-central Axel Heibeng Island, Nunavut, Canada. Journal of Canada Earth Sciences 44, 1653–59.CrossRefGoogle Scholar
Wang, F, Xu, WL, Li, J, Pei, FP and Cao, HH (2009) Chronology and geochemistry of Early Cretaceous gabbro-diorite in Yantongshan area of Jilin. Global Geology 28, 403413 (in Chinese with English abstract).Google Scholar
Wang, F, Xu, WL, Xing, KC, Wang, YN, Zhang, HH, Wu, W, Sun, CY and Ge, WC (2019) Final closure of the Paleo-Asian Ocean and onset of subduction of Paleo-Pacific Ocean: Constraints from Early Mesozoic magmatism in the central southern Jilin province, NE China. Journal of Geophysical Research: Solid Earth 124, 2601–22.CrossRefGoogle Scholar
Wang, F, Xu, YG, Xu, WL, Yang, L, Wu, W and Sun, CY (2017) Early Jurassic calc-alkaline magmatism in Northeast China: Magmatic response to subduction of the Paleo-Pacific plate beneath the Eurasian continent. Journal of Asian Earth Sciences 143, 249–68.CrossRefGoogle Scholar
Wang, T, Zheng, YD, Zhang, JJ, Zeng, LS, Donskaya, T, Guo, L and Li, JB (2011) Pattern and kinematic polarity of late Mesozoic extension in continental NE Asia: perspectives from metamorphic core complexes. Tectonics 30, 127.CrossRefGoogle Scholar
Wang, ZH, Yang, H, Ge, WC, Bi, JH, Zhang, YL and Xu, WL (2016a) Discovery and geological significance of the Eocene granodiorites in the Sanjiang Basin, NE China: Evidence from zircon U–Pb chronology, geochemistry and Sr-Nd-Hf isotopes. Acta Petrologica Sinica 32, 1823–38.Google Scholar
Wang, ZJ, Xu, WL, Pei, FP, Wang, ZW and Cao, HH (2015) Geochronology and geochemistry of middle Permian-Middle Triassic intrusive rocks from central-eastern Jilin Province, NE China: constraints on the tectonic evolution of the eastern segment of the Paleo-Asian Ocean. Lithos 238, 1325.CrossRefGoogle Scholar
Wang, ZW, Xu, WL, Pei, FP, Wang, F and Guo, P (2016b) Geochronology and geochemistry of early Paleozoic igneous rocks of the Lesser Xing’an Range, NE China: implications for the tectonic evolution of the eastern Central Asian Orogenic Belt. Lithos 26, 144–63.CrossRefGoogle Scholar
Wilde, SA, Wu, FY and Zhang, XZ (2003) Late Pan-African magmatism in Northeastern China: SHRIMP U–Pb zircon evidence for igneous ages from the Mashan Complex. Precambrian Research 122, 311–27.CrossRefGoogle Scholar
Wilde, SA, Zhang, XZ and Wu, FY (2000) Extension of a newly-identified 500 Ma metamorphic terrain in Northeast China: further U–Pb SHRIMP dating of the Mashan Complex, Heilongjiang Province, China. Tectonophysics 328, 115–30.CrossRefGoogle Scholar
Wu, FY, Sun, DY, Ge, WC, Zhang, YB, Grant, ML, Wilde, SA and Jahn, BM (2011) Geochronology of the Phanerozoic granitoids in northeastern China. Journal of Asian Earth Sciences 41, 130.CrossRefGoogle Scholar
Xiao, WJ, Windley, BF, Sun, S, Li, JL, Huang, BC, Han, CM, Yuan, C, Sun, M and Chen, HL (2015) A tale of amalgamation of three Permo-Triassic collage systems in Central Asia: oroclines, sutures, and terminal accretion. Annual Review of Earth and Planetary Sciences 43, 477507.CrossRefGoogle Scholar
Xu, WL, Pei, FP, Wang, F, Meng, E, Ji, WQ, Yang, DB and Wang, W (2013) Spatial-temporal relationships of Mesozoic volcanic rocks in NE China: constraints on tectonic overprinting and transformations between multiple tectonic systems. Journal of Asian Earth Sciences 74, 167–93.CrossRefGoogle Scholar
Yang, H, Ge, WC, Bi, JH, Wang, ZH, Tian, DX, Dong, Y and Chen, HJ (2018) The Neoproterozoic-early Paleozoic evolution of the Jiamusi Block, NE China and its east Gondwana connection: geochenmical and zircon U–Pb–Hf isotopic constraints from the Mashan Complex. Gondwana Research 54, 102–21.CrossRefGoogle Scholar
Yang, H, Ge, WC, Zhao, GC, Bi, JH, Wang, ZH, Dong, Y and Xu, WL (2017) Zircon U–Pb ages and geochemistry of newly discovered Neoproterozoic orthogneisses in the Mishan region, NE China: constraints on the high-grade metamorphism and tectonic affinity of the Jiamusi-Khanka Block. Lithos 268–271, 1631.CrossRefGoogle Scholar
Yang, H, Ge, WC, Zhao, GC, Dong, Y, Xu, WL, Wang, ZH, Ji, Z and Yu, JJ (2015a) Late Triassic intrusive complex in the Jidong region, Jiamusi-Khanka Block, NE China: geochemistry, zircon U–Pb ages, Lu–Hf isotopes, and implications for magma mingling and mixing. Lithos 224–225, 143–59.CrossRefGoogle Scholar
Yang, H, Ge, WC, Zhao, GC, Yu, JJ and Zhang, YL (2015b) Early Permian-Late Triassic granitic magmatism in the Jiamusi-Khanka Massif, eastern segment of the Central Asian Orogenic Belt and its implications. Gondwana Research 27, 1509–33.CrossRefGoogle Scholar
Yu, JJ, Wang, F, Xu, WL, Gao, FH and Tang, J (2013) Late Permian tectonic evolution at the southeastern margin of the Songnen-Zhangguangcai Range Massif, NE China: constraints from geochronology and geochemistry of granitoids. Gondwana Research 24, 635–47.CrossRefGoogle Scholar
Yuan, HL, Gao, S, Liu, XM, Li, HM, Günther, D and Wu, FZ (2004) Accurate U–Pb age and trace element determinations of zircon by laser ablation inductively coupled plasma mass spectrometry. Geostandard Newsletter 28, 353–70.CrossRefGoogle Scholar
Zhang, FQ, Chen, HL, Batt, GE, Dilek, Y, Na, AM, Sun, MD, Yang, SF, Meng, QA and Zhao, XQ (2015) Detrital zircon U–Pb geochronology and stratigraphy of the Cretaceous Sanjiang Basin in NE China: provenance record of an abrupt tectonic switch in the mode and nature of the NE Asian continental margin evolution. Tectonophysics 665, 5878.CrossRefGoogle Scholar
Zhang, FQ, Chen, HL, Yang, SF, Feng, ZQ, Wu, HY, Batt, GE, Zhao, XQ, Sun, MD, A, MN, Wang, SH and Yang, JG (2012) Late Mesozoic–Cenozoic evolution of the Sanjiang Basin in NE China and its tectonic implications for the West Pacific continental margin. Journal of Asian Earth Sciences 49, 287–99.CrossRefGoogle Scholar
Zhang, FQ, Dilek, Y, Chen, HL, Yang, SF and Meng, QA (2017) Structural architecture and stratigraphic record of Late Mesozic sedimentary basins in NE China: Tectonic archives of the Late Cretaceous continental margin evolution in East Asia. Earth-Science Reviews 171, 598620.CrossRefGoogle Scholar
Zhang, KJ, Yan, LL and Ji, C (2019) Switch of NE Asia from extension to contraction at the mid-Cretaceous: A tale of the Okhotsk ocean plateau from initiation by the Perm Anomaly extrusion in the Mongol-Okhotsk ocean? Earth-Science Reviews 198, 102941.CrossRefGoogle Scholar
Zheng, H, Sun, XM, Wan, K, Wang, PJ, He, S and Zhang, XQ (2019) Structure and tectonic evolution of the Late Jurassic-Early Cretaceous Wandashan accretionary complex, NE China. International Geology Review 61, 1738.CrossRefGoogle Scholar
Zhou, JB and Li, L (2017) The Mesozoic accretionary complex in Northeast China: evidence for the accretion history of Paleo-Pacific subduction. Journal of Asian Earth Sciences 145, 91100.CrossRefGoogle Scholar
Zhou, JP, István, D, Liu, YJ, Li, WM and Eynatten, HV (2020) Late Cretaceous-Tertiary tectonic inversion of northeastern Asian continental margin: Insight from the low temperature thermochronology in NE China. Gondwana Research, https://doi.org/10.1016/j.gr.2020.05.017.CrossRefGoogle Scholar
Figure 0

Fig. 1. Topographic and tectonic map of (a) northeastern Asia and (b) northeastern China (after Tao et al.2014); (c) aeromagnetic anomaly map of northeastern China; and (d) geological map of northeastern China showing the location of the Sanjiang Basin and the ages of igneous rocks (after Wu et al.2011; Xu et al.2013).

Figure 1

Fig. 2. Cross-section and stratigraphy of the Sanjiang Basin (Zhang et al.2017). Locations of seismic profiles are shown in Figure 1c.

Figure 2

Fig. 3. (a) Detailed geological map of study area; and (b) measured geological sections of the Dongrong and Yanwo formations.

Figure 3

Fig. 4. Photographs showing (a) intrusive relationship; and (b) 50–55 Ma granite. Photomicrograph of (c) Cladophlebis fossil; (d) Podozamites fossil; (e) selected sandstone samples (cross-polarized light); (f) andesitic fragment; (g) felsic fragment; and (h) metamorphosed felsic lithic fragment. Lv – volcanic lithic fragment; Lm – metamorphic lithic fragment.

Figure 4

Fig. 5. Lithological classification of siliciclastic rocks (Herron, 1988). The grey symbols correspond to samples from Cao et al. (2015) and Hu et al. (2021).

Figure 5

Fig. 6. (a) CIA–ICV diagram for sedimentary rocks from the Sanjiang Basin (Nesbitt & Young, 1989; Cox et al.1995). (b) A–CN–K diagram for sedimentary rocks from the Sanjiang Basin (Fedo et al.1995). Symbols as for Figure 5. The grey symbols correspond to samples from Cao et al. (2015) and Hu et al. (2021).

Figure 6

Fig. 7. (a) Chondrite-normalized REE patterns and (b) primitive mantle (PM) normalized trace-element diagrams. The chondrite and PM values used for normalization are from Boynton (1984) and Sun & McDonough (1989), respectively.

Figure 7

Fig. 8. Representative cathodoluminescence (CL) images of selected zircons and zircon U–Pb concordia diagrams of the Dongrong and Yanwo formations. Red and dotted circles on the zircon grains indicate locations of the analysis spots used for LA-ICP-MS U–Pb dating and Hf isotopic analyses, respectively.

Figure 8

Fig. 9. Correlations between Hf isotopic compositions and ages of detrital zircons from the Dongrong and Yanwo formations. J – Jurassic; T – Triassic; P – Permian; C – Carboniferous; D – Devonian; O – Ordovician.

Figure 9

Fig. 10. Th/U–CIA diagrams for the Dongrong and Yanwo formations. Symbols as for Figure 5.

Figure 10

Fig. 11. (a) Discriminant function diagram of major-element provenance (Roser & Korsch, 1986). (b) Eu/Eu* versus Th/Sc (after McLennan et al.1990). (c) Co/Th versus La/Sc (after McLennan et al.1983). (d) La-Th-Sc diagram for tectonic setting discrimination (after Bhatia & Crook, 1986). (e) Tectonic setting of sandstone derivation from different types of provenances (after Dickinson, 1985). (f) Cumulative proportions versus time (in Ma) of detrital zircon ages from the Sanjiang Basin (after Cawood et al.2012). Discriminant Function 1 = (−1.773×TiO2) + (0.607×Al2O3)  +  (0.760×Fe2O3)  +  (−1.500×MgO)  +  (0.616×CaO) + (0.509×Na2O) + (−1.224×K2O) + (−9.090); Discriminant Function 2 = (0.445×TiO2) + (0.070×Al2O3) + (−0.250×Fe2O3) + (−1.142×MgO) + (0.438×CaO) + (1.475×Na2O) + (1.462×K2O) + (−6.861). Sample symbols as for Figure 5. The grey symbols correspond to samples from Cao et al. (2015) and Hu et al. (2021).

Figure 11

Fig. 12. (a) Comparison of detrital zircon ϵHf(t) and zircon Hf isotope compositions of igneous rocks from the study area. Reference data are from Bi et al. (2014a, b, 2015, 2016, 2017), Cao et al. (2013), Guo et al. (2016, 2018), Long et al. (2019, 2020), Luan et al. (2017b, 2019), Wang et al. (2009, 2015, 2016a, b, 2017, 2019) and Yang et al. (2015a, b, 2017, 2018).

Figure 12

Fig. 13. Comparison of probability diagrams of detrital zircon ages for sedimentary sequences from the Sanjiang Basin. Reference data for the Didao, Chengzihe, Muling and Hailang formations from Zhang et al. (2015).

Figure 13

Fig. 14. Schematic illustration of the Early Jurassic and Cretaceous uplift events in the northeastern Asian continental margin.

Supplementary material: File

Luan et al. supplementary material

Luan et al. supplementary material 1

Download Luan et al. supplementary material(File)
File 208.9 KB
Supplementary material: File

Luan et al. supplementary material

Luan et al. supplementary material 2

Download Luan et al. supplementary material(File)
File 196.5 KB
Supplementary material: File

Luan et al. supplementary material

Luan et al. supplementary material 3

Download Luan et al. supplementary material(File)
File 184.6 KB
Supplementary material: File

Luan et al. supplementary material

Luan et al. supplementary material 4

Download Luan et al. supplementary material(File)
File 995 KB