Origin of Archean Continental Crust in the Tarim Craton: Progresses and Issues
-
摘要:
塔里木克拉通是中国三大古老陆块之一，但由于大面积沉积覆盖，其古老基底的形成演化研究程度较低。然而，近年来在塔里木盆地周缘的库鲁克塔格、敦煌、北阿尔金、铁克里克等地，以及塔里木盆地基底钻井岩心中，均发现了太古宙岩石，表明其可能普遍存在太古宙基底。笔者对塔里木克拉通太古宙岩石的研究历史和最新进展进行了简要总结，对太古宙大陆地壳的形成时间、机制和动力学背景进行了讨论，并指出了未来研究的方向。结果表明，塔里木克拉通的太古宙大陆地壳形成似乎具有区域差异性，北部的库鲁克塔格、敦煌和北阿尔金地区广泛发育新太古代岩浆作用，峰期为～2.5 Ga和～2.7 Ga，北阿尔金地区～3.7 Ga岩石的发现为塔里木克拉通始太古代陆核的存在提供了可靠证据，而西南部的铁克里克地区和盆地基底以中太古代（3.2～2.8 Ga）地壳生长和再造为特征，目前尚未发现新太古代岩石。地球化学、热力学模拟和锆石氧逸度–湿度计研究表明，太古宙大陆地壳可能是不同源岩在不同深度（压力）通过水致熔融产生的，形成于俯冲相关构造背景，而早期板块构造自始太古代以来就已运行。太古宙大陆地壳物质组成的确定、变质–变形的识别、岩浆形成物理化学条件的厘定等方向仍是未来塔里木克拉通太古宙地质研究的重点。
Abstract:The Tarim Craton is one of the three ancient continental blocks in China, but the formation and evolution of its ancient basement have been poorly studied due to extensive sedimentary cover. However, in recent years, Archean rocks have been found in the Kuruktag, Dunhuang, North Altyn Tagh, Tiekelike areas on the periphery of the Tarim Basin, as well as in the drill core from its basement, indicating that there may be a widespread Archean basement. In this paper, the study history and recent progress of Archean rocks in the Tarim Craton are summarized, the formation time, mechanism and geodynamics of Archean continental crust are discussed, and the future research direction is pointed out. The results show that the formation of the Archean continental crust in the Tarim Craton appears to have regional differences. Neoarchean magmatism was widely developed in the Kuruktag, Dunhuang and North Altyn areas, with peaks of ～2.5 Ga and ～2.7 Ga. The discovery of ～3.7 Ga rocks in the North Altyn Tagh area provides reliable evidence for the existence of an Eoarchean continental nucleus in the Tarim Craton. The Tiekelike area and basin basement in the southwest Tarim are characterized by Mesoarchean (3.2～2.8 Ga) crustal growth and reworking, and no Neoarchean rocks have been found. Geochemistry, thermodynamic modelling and zircon oxybarometer-hygrometer indicate that the Archean continental crust might have been produced by water-induced melting of different source rocks at different depths (pressures) and formed in subduction-related tectonic settings, indicating that early plate tectonics have been in operation since the Eoarchean. The elucidation of the components of the Archean continental crust, the identification of metamorphism and deformation, and the determination of the physical and chemical conditions of magma formation are still the focus of future studies of the Archean geology in the Tarim Craton.
-
Keywords:
- Archean /
- continental crust /
- plate tectonics /
- zircon /
- TTG gneisses /
- Tarim craton
-
滑坡是一种破坏性事件，往往对人类生命、财产和生存环境构成严重威胁，严重制约着人类的可持续发展。据统计，发展中国家因滑坡造成的死亡率高达95%，造成的经济损失占国民生产总值的1%～2%（Youssef et al., 2022）。作为发展中国家，中国山区占国土面积的70%，地质环境条件复杂、构造活动强烈、工程建设规模大，滑坡多发且分布较广（刘光旭等，2014）。
中国西北部的黄土覆盖区生态极其脆弱，土壤结构松散，降水集中，水土流失严重，已成为全国最容易发生滑坡的地区（孙萍萍等，2022）。此外，滑坡对当地社会经济的发展和人类安全构成极大威胁（韩玲等，2019）。因此，滑坡预测对该区域滑坡的防治具有重要意义。然而，传统的野外调查难以满足滑坡预测的要求。随着地理信息技术（GIS）的发展，滑坡易发性分区的研究结合了许多统计学方法，即统计指数（Abedini et al., 2018）、逻辑回归（Thai et al., 2018）、确定性因素（Li et al., 2017）、层次分析法（张向营等，2018）、频率比（Abedini et al., 2018）、熵权法（Jaafari et al., 2014）和证据权法。虽然统计学方法和GIS技术可以很好地结合起来，更适合大面积研究，但滑坡易发性分区是一种典型的复杂非线性问题，所以统计学方法得到的结果往往精度较低。因此，越来越多的机器学习方法被引入，如人工神经网络（谢振华等，2017）、核逻辑回归（Chen et al., 2020）、随机森林（Behnia et al., 2018）、支持向量机（武雪玲等，2016；韩玲等，2019） 、朴素贝叶斯（Sameen et al., 2020）和决策树（Kim et al., 2018）。机器学习模型的性能通常取决于训练数据的质量和数量，并且对建模过程中的参数十分敏感。随机划分的数据、不同的研究规模和不同的数据分辨率都会对训练数据产生影响（Bazi et al., 2021）。此外，对于哪种方法最适合研究区开展滑坡易发性分区，目前还没有明确的结论。
近年来，混合模型的出现为滑坡易发性分区研究提供了更多的选择，也可探索出更加科学合理的结果。这些模型包括频率密度耦合核函数逻辑回归模型（FR-KLR）（张庭瑜等，2020）、确定性系数耦合径向基神经网络模型（CF-RBFNN）（Li et al., 2017）等。此外，也有学者利用分形理论对滑坡进行研究，但目前研究的主要方向是用分形维数来描述滑坡的空间分布，几乎没有对滑坡易发性分区进行过研究。
由于评估方法会对滑坡易发性分区结果产生很大影响，探索和研究分区模型是一项重要的工作，并且将分形维数与IOE和SVM融合的研究鲜有报道。因此，本研究基于分形维数和两种原始模型（IOE、SVM），分别构建两种混合模型（F-IOE、F-SVM）对宝鸡市北部典型黄土覆盖区进行滑坡易发性分区研究。研究结果可为研究区以及类似区域的灾害防治提供理论依据，并且相关图件的发布也可为当地的防灾减灾工作提供决策支持。
1. 研究区域
研究区位于陕西省宝鸡市北部典型黄土覆盖区，面积约
2451 km2，地理坐标介于E 107°00′01″～107°45′08″和N 34°29′22″～34°50′38″。由南向北横跨黄土台地和黄土高原两大地貌类型。黄土台地海拔750～950 m，相对高差小于200 m，黄土高原地区海拔900～1700 m，相对高差100～300 m。研究区域内发育一条水系，即千河，水流方向为WN—ES，最终流入渭河。研究区气候类型属暖温带大陆性气候，年平均气温和月平均气温分别为12 ℃和25 ℃，年平均降水量为653 mm，降雨期主要集中在6～9月，近十年平均月最大降雨量为103 mm。由于黄土遇水易塌陷，故该区域容易发生地质灾害。研究区横跨黄土高原、鄂尔多斯盆地和秦岭造山带，由南向北跨越3个构造单元，区域内发育7条断裂带，其中F1-F5断裂走向为NW—SE，F6-F7断裂走向为W—E。根据《陕西省区域地质志》来看，研究区未发生破坏性地震（陕西省地质矿产局，1989）。
2. 数据应用
2.1 滑坡编录图
包含滑坡数量、位置和规模信息的滑坡编录图是滑坡易发性分区研究的基础。根据野外调查、前期地质资料、1980～2017年的历史滑坡数据以及航空遥感数据，共提取地质灾害179处，其中堆积层滑坡166处、基岩滑坡13处。滑坡周界面积的最大、最小和平均尺寸分别为2.2×104 m2、1.1×103 m2、8.3×103 m2。由于近90%的滑坡周界面积不超过
10000 m2，因此采用质心法将多边形转化为点，然后将这些点映射到研究区的滑坡编录图中（图1）。![]()
图 1 研究区位置图（a）和滑坡编录图（b）Figure 1. （a）Location of study area and（b） landslide inventory map2.2 训练和测试数据集的制备
为了建立滑坡易发性分区模型并获得滑坡易发性分区图（LSM），训练和测试数据集的制备是重要且必要的步骤。由于建立滑坡易发性模型必须同时存在正样本和负样本，本研究将179个滑坡样本全部作为正样本。此外，根据地质资料和野外观测数据分析，研究区滑坡一般不会发生在坡度小于5°的区域 （张茂省等，2019），因此将这些区域定义为非滑坡区（Pandey et al., 2020）。在非滑坡区随机生成相同数量的非滑坡点作为负样本，然后将179个滑坡点以70%与30%的比例随机分为两组，同时将179个非滑坡点按照相同的比例分为两组。最后，将70%（250个）的样本用来构建模型，30%（108个）的样本用于模型测试。
由于单一的模型会因数据的随机性而出现过拟合的现象（Mahdadi et al., 2018），因此为了减少随机选择对滑坡易发性分区模型性能的影响，将所有滑坡点以7∶3的比例随机划分30次，并对每组训练数据集进行10折交叉验证。保留最稳定的训练数据集检测多重共线性问题、确定滑坡影响因子和建立滑坡易发性分区模型。
2.3 滑坡影响因子分析
滑坡的发生受区域地质环境、水文、气象和构造运动等一系列条件因素的影响，影响因子的选择和分类将直接影响滑坡易发性分区的准确性（Hong et al., 2017）。本研究在实地调查和以往研究的基础上，共考虑了12个内部和外部的影响因子，包括坡度、坡向、高程、降雨量、平面曲率、距河流的距离、土地利用类型、距道路的距离、距断层的距离、归一化植被指数（NDVI）、地层岩性和剖面曲率。
滑坡与坡度的关系密切，坡度直接影响坡内应力分布及有效面。研究区的坡度信息由分辨率为30 m的DEM数据计算得出，并以5°的等间隔划分为6类（图2a）。坡向被认为是另一个重要的因素，它通过影响土壤的蒸发量、植被的生长和侵蚀来影响边坡的稳定性。研究区的坡向信息也从DEM中提取，并重分为9类（图2b）。同样地，从DEM中提取高程数据，按照
1000 m的间隔重分为7类（图2c）。![]()
图 2 各影响因子图（a）.坡度；（b）.坡向；（c）.高程；（d）.距河流的距离；（e）.距断层的距离；（f）.距道路的距离；（g）.降雨量；（h）.地层岩性；（i）.NDVI；（j）.土地利用类型；（k）.平面曲率；（l）.剖面曲率Figure 2. Map of various Conditioning factors许多学者发现河流侵蚀在滑坡发生中起重要作用（Kumar et al., 2017）。在本研究中，距河流的距离被认为是重要的影响因素之一，并以100 m的等间隔生成了5个不同的缓冲区（图2d）。构造作用也是影响滑坡发生的重要因素。构造作用不仅是单一滑坡发生的必要条件，也是区域滑坡发生的直接控制因素。因此，以距断层的距离作为影响因子，利用ArcGIS软件以
1000 m的相等间隔生成5个不同的缓冲区（图2e）。此外，研究发现道路建设需要开挖大量的边坡，这很可能在研究区产生潜在的滑坡。因此，滑坡易发性分区的另一个条件因素是距道路的距离，以100 m的等间隔生成了4个不同的道路缓冲区（图2f）。根据以往的文献，降雨量对滑坡的稳定性有很大影响（Bui et al., 2017）。本研究使用的降雨数据提取自当地气象台提供的2007～2017年的年平均降雨量。将降雨数据按照150 mm的等间隔分为5类（图2g）。地层岩性决定了机械强度、耐候性和应力分布，进而影响边坡的稳定性和侵蚀强度。根据研究区（表1）的内容，将地层岩性重新划分为6组（图2h）。表 1 研究区地层岩性单元表Table 1. Lithological units of study area类别 地质年代 编码 主要地层岩性 A 新进纪 Q4 砂，砾石，黄土 更新纪 Q3 黄土，砾石 B 上新统 N2 砂土 中新统 N1 石英砂，黏土 C 晚白垩纪 K1 泥岩，砂质泥岩，泥质砂岩 D 早侏罗纪 J3 块状聚集物，粘胶岩，粉砂质泥岩 中侏罗纪 J2 石质，泥岩，粉质泥岩 晚侏罗纪 J1 粉砂岩，煤层 E 早三叠纪 T3 泥岩，页岩，煤层 中三叠纪 T2 中细砂岩，粉砂岩，泥岩 晚三叠纪 T1 石质，细砂岩，粉砂岩，砂质泥岩 F 二叠纪 P 砂质泥岩，细砂岩，粉砂岩 植被对滑坡的发生也有显著影响（Arabameri et al., 2018）。NDVI一直被用作反映地表植被覆盖程度的重要指标，本研究利用ENVI软件，基于GF-2多光谱遥感影像计算NDVI，采用自然断点法（图2i）将NDVI分为5类。土地利用类型通常被认为是滑坡发生的基本影响因素之一，本研究利用ArcGIS软件将土地利用类型分为5组（图2j）。此外，基于DEM数据计算出平面曲率和剖面曲率，并分别使用自然断点法（图2k、图2l）分为5个部分。最后，为了保持数据分辨率的统一，基于ArcGIS软件，将这12个影响因子图全部重采样为30 m空间分辨率。
3. 研究方法
本研究主要包括4个部分：第1部分是基础数据的制备，包括滑坡编录图的构建、滑坡影响因子的分析以及训练和测试数据集的准备；第2部分涉及滑坡影响因子的潜在多重共线性问题的判断、滑坡影响因子的优化以及分形维数的计算；第3部分主要建立了4种滑坡易发性分区模型，即熵权模型（IOE）、基于分形维数的熵权模型（F-IOE）、支持向量机模型（SVM）和基于分形维数的支持向量机模型（F-SVM）。然后根据上述4种模型生成研究区的滑坡易发性分区图；第4部分是利用统计学指标评价滑坡易发性模型的性能，并根据ROC特征曲线对这4种模型的泛化性进行比较。接下来简要介绍本研究中采用的方法。
3.1 多重共线性
多重共线性是指由于影响因子之间存在高度相关性，模型被扭曲或难以估计（Islam et al., 2011）。为了检测多重共线性的存在，本研究引入了方差膨胀因子（VIF）。VIF表示影响因子之间存在潜在多重共线性问题时的方差与不存在潜在多重共线性问题时的方差之比，VIF的倒数是公差（TOL）。通常，VIF值与影响因子之间的多重共线性问题的强度成正比。当VIF＞10且TOL＜0.1时，认为数据集中存在潜在的多重共线性问题。
3.2 影响因子优选
在滑坡易发性分区制图中，并非所有影响因子在建模中都具有相同的预测能力，并且可能存在噪声，这将使结果产生误差。因此，可以计算信息增益率（IG）来定量反映每个影响因子对滑坡易发性模型的贡献，IG值越高，表明滑坡易发性映射的预测能力越高（Chen et al., 2017）。如果IG值等于0，这意味着影响因子没有预测能力，应放弃。IG值计算公式如下：
$$ Entropy(D) = - \sum\limits_{k = 1}^{\left| y \right|} {{p_k}} \log _2^{{p_k}} $$ (1) 其中：D是训练数据集；熵（D）表示训练数据集的熵；y是D中的物种数；$ p_{k} $表示D中k类的比例。然后使用表示影响因子之一的s将训练数据集划分为Dv（v=1，2，3，···，m），并通过等式（2）计算增益（D，s）。
$$ Gain(D,s) = Entropy(D) - \sum\limits_{v - 1}^{\left| m \right|} {\dfrac{{\left| {{D_v}} \right|}}{D}} Entropy({D_v}) $$ (2) 调节因子的IG值计算为：
$$ IG(D,s) = - \frac{{Gain(D,a)}}{{IV(a)}} $$ (3) 其中：IV（s）可通过等式（4）获得。
$$ IV(s) = - \sum\limits_{v = 1}^m {\frac{{\left| {{D_m}} \right|}}{{\left| D \right|}}} {\log _2}_{}^{\tfrac{{\left| {{D_m}} \right|}}{{\left| D \right|}}} $$ (4) 3.3 分形维数
分形维数是分形理论中用来定量描述分形特征和几何复杂性的参数（Shuren et al., 2005）。本研究引入分形维数来表示滑坡的聚类程度和滑坡系统的复杂性。在当前研究中，使用ArcGIS软件将研究区域划分为边长为r的n个网格。然后将r值减少一半。在本研究中，r值分别设定为
10000 m、5000 m、2500 m、1250 m、625 m和312.5 m。计算滑坡点位于不同边长的每个网格中的网格数（Nr），分形维数计算示意图见图3。然后，使用最小二乘法拟合直线，如下所示：$$ \ln {N_r} = a + {f_j}\ln r $$ (5) 其中：$ a $是1个参数；$ {f_j} $表示分形维数（Bouboulis et al., 2006）。最后，通过盒计数技术测量每个影响因子的分形维数，并将其作为滑坡易发性建模的输入数据。
3.4 熵权模型（IOE）和基于分形维数的熵权模型（F-IOE）
IOE作为二元线性统计分析模型被广泛应用（Li et al., 2019）。原始IOE模型通过计算滑坡的频率比来计算每个影响因子的权重，以建立滑坡易发性模型。具体计算过程如下：
$$ F{R_{ij}} = \frac{{{b_{ij}}}}{{{a_{ij}}}} $$ (6) $$ {P_{ij}} = \frac{{F{R_{ij}}}}{{\displaystyle \sum\nolimits_{j = 1}^{{N_j}} {F{R_{ij}}} }} $$ (7) $$ {H_j} = - \sum\limits_{i = 1}^{{N_j}} {{P_{ij}}} {\log _2}{P_{ij}},j = 1,2,3,...,n $$ (8) $$ {H_{j\max }} = {\log _2}{N_j} $$ (9) $$ {I_j} = \frac{{{H_{j\max }} - {H_j}}}{{{H_{j\max }}}} $$ (10) $$ {W_j} = {I_j} \times F{R_{ij}} $$ (11) 其中：$ {a_{ij}} $和$ {b_{ij}} $分别代表像元百分比和滑坡百分比；$ {N_j} $是影响因子的分类数；FRij是滑坡发生的频率；$ {H_j} $和$ {H_{j\max }} $是熵值，$ {W_j} $代表影响因子的综合评分。对于新的混合模型（F-IOE），使用每个因子的分形维数$ {f_j} $代替$ {H_j} $作为建模数据。影响因子的最终权重由F-$ {W_j} $表示（贾俊等，2023）。
3.5 支持向量机模型（SVM）和基于分形维数的支持向量机模型（F-SVM）
支持向量机（SVM）被认为是解决二值分类问题的有效监督学习算法（Wang et al., 2021）。支持向量机模型的基本理论是在特征空间中找到最佳分离超平面，该超平面可以最大化训练数据集中滑坡和非滑坡数据之间的样本间隔（Balogun et al., 2021），可以使用以下数学公式表示：
$$ P = \frac{1}{2}{\left\| w \right\|^2} $$ (12) 约束条件为：
$$ y{}_i\left( {\left( {w \times {x_i}} \right) + b} \right) \geqslant 1 $$ (13) 其中：$ b $为标量；$ {\left\| w \right\|^2} $是超平面法向量的范数。引入拉格朗日乘数法则求极值（$ {\lambda _i} $），生成辅助函数如下：
$$ L = \frac{1}{2}{\left\| w \right\|^2} - \sum\limits_{i - 1}^n {{\lambda _i}} \left( {{y_i}\left( {\left( {w \times {x_i}} \right) + b} \right) - 1} \right) $$ (14) 对于不可分离的情况，引入松弛变量$ \xi_{t} $（i=1，2，···n），等式（12）和等式（13）可以用如下代替：
$$ {y_i}\left( {\left( {w \times {x_i}} \right) + b} \right) \geqslant 1 - {\xi _i} $$ (15) $$ L = \frac{1}{2}{\left\| w \right\|^2} + C\sum\limits_{i = 1}^n {{\xi _i}} \text{，} \left( {C = - \frac{1}{{vn}}} \right) $$ (16) 其中：C为约束系数；$ v $∈（0,1]代表误分类罚值。此外，非线性决策边界可以使用核函数$ K\left(x_{i}, x_{j}\right) $运算。本研究采用径向基函数（RBF）进行滑坡易发性建模，其数学公式如下：
$$ K\left( {{x_i},{x_j}} \right) = \exp \left( { - \delta {{\left\| {{x_i} - {x_j}} \right\|}^2}} \right) \text{，} \delta ＞ 0 $$ (17) 其中：δ 说明了高斯核函数的宽度。
同样，使用分形维数$ {f_j} $替换每个影响因子的原始数值数据作为输入数据，然后应用等式（12）在高维空间中使用超平面对输入数据进行分类。最后，构建一种新的混合模型（F-SVM）。
3.6 统计学指标
3.6.1 PPR，NPR和ACC
本研究采用3种统计指标，即阳性预测率（Positive Prediction Rate, PPR）、阴性预测率（Negative Prediction Rate, NPR）和准确性（Accuracy, ACC）来评估滑坡易发性制图的结果。采用PPR评价滑坡易发性模型对滑坡的预测能力，采用NPR评价滑坡易发性模型对非滑坡的预测能力。采用ACC测量滑坡易发性分区的准确性。这3个统计指标的计算过程如下：
$$ PPR = \frac{{TP}}{{TP + FP}} $$ (18) $$ NPR = \frac{{TN}}{{TN + FN}} $$ (19) $$ ACC = \frac{{TP + TN}}{{TP + TN + FP + FN}} $$ (20) 其中：TP全名为真阳性；TN全名为真阴性，表示正确分类的像素数；FP的全名为假阳性；FN的全名为假阴性，表示未正确分类的像素数（杨光等，2019）。一般来说，ACC值越大，滑坡易发性分区的精度越高（Pham et al., 2021）。
3.6.2 受试者工作特征曲线（ROC）
对4种模型（IOE、F-IOE、SVM和F-SVM）获得的分类结果进行比较评估是一个重要的部分（张文龙等，2023）。在本研究中，引入了ROC曲线来比较这4种滑坡易发性模型，揭示了使用合成法的灵敏度和特异度值的关系。它为连续变量设置几个不同的临界值，以计算一系列灵敏度和特异度，然后绘制以灵敏度作为Y轴1-特异度作为X轴的曲线。曲线下面积（简称AUC）用于比较本研究中每个模型的性能。AUC值范围为0~1。如果AUC＝1，则表明滑坡和非滑坡的分类完全正确(Lombardo et al., 2020)。相反，如果AUC＝0，则分类完全错误。计算公式如下：
$$ \text { 灵敏度 }=\frac{T P}{T P+F N} $$ (21) $$ 1-\text { 特异度 }=1-\frac{T N}{T N+F P} $$ (22) $$ A U C=\frac{(\Sigma T P+\Sigma T N)}{P+N} $$ (23) 其中：P和N分别表示研究区域内滑坡和非滑坡的总数。
4. 结果
4.1 多重共线性诊断结果
在本研究中，VIF和TOL被用于检测影响因子之间潜在的多重共线性问题。从计算结果（表2）可以看出，VIF的最大值为NDVI（VIF=1.675），公差的最小值为土地利用类型（公差=0.650）。所有影响因子的VIF和TOL不在潜在的多重共线性问题阈值内（VIF＞10，TOL＜0.1），表明这12个影响因子中不存在多重共线性问题。
表 2 影响因子的VIF和TOL表Table 2. Variance inflation factors (VIF) and tolerances of each conditioning factor影响因子 TOL VIF 坡度 0.934 1.071 坡向 0.926 1.080 高程 0.656 1.525 距河流的距离 0.908 1.101 距道路的距离 0.877 1.141 距断层的距离 0.916 1.092 NDVI 0.597 1.675 土地利用类型 0.650 1.538 地层岩性 0.814 1.228 降雨量 0.814 1.229 平面曲率 0.912 1.096 剖面曲率 0.925 1.082 4.2 影响因子优选
本研究通过10折交叉验证计算平均信息增益率（AM）来选择影响因子。从图4可以看出，所有影响因子的AM中，高程的最高（0.598），其次是坡度和坡向的相同（0.299），距河流的距离（0.253），降雨量（0.230），土地利用类型（0.210），距道路的距离（0.184），距断层的距离和NDVI的相同（0.115），剖面曲率（0.046）和地层岩性（0.023）。然而，平面曲率的AM值为0，这表明平面曲率因子对建模的贡献为0，可能会造成干扰。因此，在后续建模和评估中排除平面曲率。
![]()
图 4 各影响因子的平均信息增益率图Figure 4. Average Information Gain ratio of each conditioning factor4.3 分形维数的计算
根据分形理论，为了计算影响因子的分形维数，研究区滑坡的空间分布必须满足分形特征（自相关性）。根据上述方法，使用公式（5）获得研究区内滑坡的空间分形维数。从图5可以看出，研究区滑坡的空间分维数为
1.9032 ，R2为0.9649 ，这表明滑坡的空间分布具有明显的分形特征和显著的自相关性。因此，继续计算每个影响因子的空间分维是有意义的，每个因子的分维数$ {f_j} $见表3。表 3 影响因子与滑坡的空间关系表Table 3. Spatial relationship between influencing factors and landslides影响因子 等级 fj FRij Pij Hj Hjmax Ij Wj F-Wj 坡度（°） ＜5 0.0772 0.3035 0.0367 2.2857 2.5850 0.1158 0.1595 0.1201 5~10 0.1906 0.8709 0.1053 10~15 0.2285 0.9124 0.1104 15~20 0.2372 1.1518 0.1393 20~25 0.1682 1.9427 0.2350 ＞25 0.2284 3.0864 0.3733 坡向 水平 0.0000 0.0000 0.0000 2.9004 3.1699 0.0850 0.0732 0.0946 北 0.0459 0.6097 0.0787 东北 0.0317 0.3564 0.0460 东 0.2239 0.9802 0.1266 东南 0.1567 1.0364 0.1338 南 0.1429 1.1387 0.1470 西南 0.1990 1.0215 0.1319 西 0.2260 1.6248 0.2098 西北 0.1251 0.9771 0.1262 高程（m） ＜850 0.5803 2.9085 0.4024 2.4666 2.8074 0.1214 0.1253 0.1279 850~950 0.3346 1.1898 0.1646 950~ 1050 0.2131 0.7202 0.0996 1050 ~1150 0.1009 0.5549 0.0768 1150 ~1250 0.0942 0.5738 0.0794 1250 ~1350 0.1606 0.8630 0.1194 ＞ 1350 0.2092 0.4174 0.0578 距河流的距离（m） ＜200 0.3998 2.2738 0.4703 2.0154 2.3219 0.1320 0.1277 0.1009 200~400 0.2220 0.8336 0.1724 400~600 0.0539 0.4464 0.0923 600~800 0.0706 0.8641 0.1787 ＞800 0.1942 0.4171 0.0863 距断层的距离（m） ＜2000 0.4358 1.2216 0.2275 2.2883 2.3219 0.0145 0.0155 0.1263 2000 ~4000 0.2725 1.2343 0.2299 4000 ~6000 0.3286 1.1496 0.2141 6000 ~8000 0.3992 1.1191 0.2084 ＞ 8000 0.2951 0.6445 0.1200 距道路的距离（m） ＜100 0.4866 1.1476 0.5302 1.3723 2.0000 0.3139 0.1698 0.1030 100~200 0.1407 0.7709 0.3562 200~300 0.1334 0.2458 0.1130 ＞300 0.0000 0.0000 0.0000 4.4 基于IOE和F-IOE模型的滑坡易发性分区
使用IOE和F-IOE模型的最终计算数据见表3。如结果所示，每个影响因子的IOE值和分形维数具有基本相同的趋势。最大$ {W_j} $（
0.3646 ）属于NDVI，其次是土地利用类型（0.2462 ），距道路的距离（0.1698 ），最小$ {W_j} $（0.0155 ）属于距断层的距离。最大F-$ {W_j} $（0.1388 ）属于剖面曲率，其次是NDVI（0.1337 ），降雨量（0.1319 ），最小F-$ {W_j} $（0.0946 ）属于坡向。从以上数据可以看出，NDVI在IOE和FIOE模型的建模过程中具有很高的权重。造成这种现象的原因是研究区水土流失严重，降雨不足。基于IOE和F-IOE模型计算滑坡易发性指数（LSI）。LSI值在0到1之间，越接近1，滑坡发生的概率越高，反之亦然。最后，将LSIIOE和LSIF-IOE划分为5个区间，以生成滑坡易发性图（LSM），同时也通过自然断点法分为5类区域：极低易发区（
0.0940 ～0.2084 ）、（0.1047 ～0.1945 ）；低易发区（0.2084 ～0.2646 ）、（0.1945 ～0.2783 ）；中易发区（0.2646 ～0.3228 ）、（0.2783 ～0.3834 ）；高易发区（0.3228 ～0.3893 ）、（0.3834 ～0.5022 ）；极高易发区（0.3893 ～0.6244 ）、（0.5022 ～0.7214 ）（图6a、图6b）。![]()
图 6 滑坡易发性图（a）.IOE 模型（LSM）；（b）.F-IOE 模型（LSM）；（c）.SVM 模型（LSM）；（d）.F-SVM 模型（LSM）Figure 6. Landslide susceptibility map using4.5 基于支持向量机和F-SVM模型的滑坡易发性分析
在建立基于RBF核函数的支持向量机和F-SVM模型之前，C和δ的确定对于建模至关重要。为了获得准确且稳定的结果，基于Python环境，利用10折交叉验证和网格搜索法，确定了SVM和F-SVM的C和δ分别为（256，
0.1250 ）和（512，0.0165 ）。然后计算SVM和F-SVM模型的LSI，输出范围为0～1，越接近1，滑坡发生的概率越高，反之亦然。最后，将LSISVM和LSIF-SVM划分为5个区间，生成LSM，并通过自然断点法将其分为5类：极低易发区（0.0145 ～0.2459 ）、（0.0740 ～0.3061 ）；低易发区（0.2459 ～0.3695 ）、（0.3061 ～0.3962 ）；中易发区（0.3695 ～0.5161 ）、（0.3962 ～0.4967 ）；高易发区（0.5161 ～0.6974 ）、（0.4967 ～0.6802 ）；极高易发区（0.6974 ～0.9983 ）、（0.6802 ～0.9574 ）（图6c、图6d）。4.6 结果的评估和比较
4.6.1 模型性能评估
基于训练数据集，利用3种统计学指标（PPR、NPR和ACC）对滑坡易发性模型的性能进行评估，结果见表4。从PPR和NPR的计算结果可以看出，F-SVM模型的最高值分别为88.72%和94.02%。说明在本研究中，F-SVM模型对滑坡是否发生的预测能力最强。从ACC计算结果来看，F-SVM模型的值也最高（91.20%），其次是SVM模型（89.20%）、F-IOE模型（87.60%）和IOE模型（84.80%）。
表 4 模型性能评价统计指标计算结果表Table 4. Calculation result of statistical indicators for model performance evaluation指标 模型 IOE SVM F-IOE F-SVM 真阳性 108 113 110 118 真阴性 104 110 109 110 假阳性 21 15 16 15 假阴性 17 12 15 7 PPR（%） 83.72 88.28 87.30 88.72 NPR（%） 85.95 90.16 87.90 94.02 ACC（%） 84.80 89.00 87.60 91.20 4.6.2 模型测试
使用测试数据集和3种统计指标完成了模型测试，计算结果如表5所示。从PPR和NPR的计算结果可以看出，F-SVM模型的值最高，分别为92.73%和94.34%。
表 5 模型测试评价统计指标计算结果表Table 5. Calculation result of statistical indicators for model validation evaluation指标 模型 IOE SVM F-IOE F-SVM 真阳性 48 52 110 118 真阴性 47 48 109 110 假阳性 7 6 16 15 假阴性 6 2 15 7 PPR（%） 87.27 89.66 86.21 92.73 NPR（%） 88.68 96.00 92.00 94.34 ACC（%） 87.96 92.59 88.89 93.52 这表明，在目前的研究中，F-SVM模型在这4种模型中对滑坡和非滑坡发生的分类性能最好。对于ACC指标，F-SVM模型的值也最高（93.52%），其次是SVM模型（92.59%）、F-IOE模型（88.89%）和IOE模型（87.96%）。总的来说，这4种模型的性能对于研究区域的滑坡易发性分区相对可靠。此外，根据统计测试方法的结果，F-SVM模型在所有指标上得分最高，具有更好的分类能力。在本研究中，两种混合模型的性能都优于其原始模型。
4.6.3 模型泛化性对比
本研究通过绘制ROC曲线，计算AUC值来对比滑坡易发性模型泛化性。基于训练和测试数据集计算的AUC值分别如图7a和图7b所示。从基于训练数据集的ROC曲线来看，F-SVM模型的AUC值最高（
0.8527 ），其次是SVM模型（0.8153 ）、IOE模型（0.8057 ）和F-IOE模型（0.8054 ）。从基于测试数据集的ROC曲线来看，F-SVM模型的AUC值最高（0.9761 ），其次是F-IOE模型（0.8591 ）、SVM模型（0.7946 ）和IOE模型（0.7434 ）。结果表明，F-SVM在训练和测试数据集中表现出最强的泛化性，并且测试数据集中的AUC值明显高于其他模型。![]()
图 7 滑坡易发性模型的ROC曲线图（a）训练数据集；（b）测试数据集Figure 7. ROC curves for landslide susceptibility models5. 讨论
在本研究中，基于滑坡编录图，在陕西省宝鸡市北部黄土覆盖区应用了4种滑坡易发性模型，即IOE、F-IOE、SVM、F-SVM模型，并进行了滑坡易发性分区。结果表明，两种混合模型（F-IOE、F-SVM）对于滑坡具有更好的分类能力。虽然许多学者已经使用分形维数从微观和宏观上研究了滑坡，但基本上没有使用分形维数作为输入数据建立混合模型的研究。由于分形维数反映了滑坡更真实的空间分布，分形维数表示的不确定性不等于随机性，分形维数与影响因子之间不存在线性关系。因此，通过使用分形维数作为输入数据构建混合模型，可以实现比普通混合模型更好的效果。另一方面，分形维数可以作为未来研究中影响因子优选的指标。在滑坡易发性建模过程中，尺度和数据分辨率的差异可能会产生一些不确定性。在区域滑坡易发性分区中，通常需要选择一系列在时空尺度上具有不同变化率的影响因子。例如，降雨量和坡度最有可能发生变化，但地层岩性变化的可能性最小。在大尺度上，影响因子的空间总体特征往往更受关注，相反，在小尺度下，影响因子在空间中的整体特征通常被视为约束条件，而细节信息则更受关注。随着空间尺度的增加，变化率高的因子所带来的影响将减少甚至消除，而变化率低的因子所带来的影响将保持甚至更突出。这些原因都将导致滑坡易发性分区的结果具有不确定性，因此，笔者将更加关注尺度变化对滑坡易发性分区的影响。
从模型性能的评估结果来看，用分形维数构建的两种混合模型（F-IOE、F-SVM）的PPR、NPR和ACC值高于其原始模型（IOE、SVM）。同样，从模型测试的结果来看，两种混合模型的上述3个统计学指标也高于其原始模型。
从图7b可以看出，两种混合模型的AUC值分别为
0.8591 （F-IOE）和0.9761 （F-SVM），显著高于原始模型。因此，构建以分形维数为输入数据的混合模型来研究滑坡易发性比原始模型更有优势。然而，研究区域不足以反映整个方法的普遍适用性，需要结合更多的区域进行研究。6. 结论
（1）本研究基于179个滑坡样本点，随机选择70%的滑坡点进行训练，其余30%的滑坡点用于测试。选择了12个影响因子，包括坡度、坡向、高程、距河流的距离、距道路的距离、距断层的距离、NDVI、土地利用类型、地层岩性、降雨量、平面曲率和剖面曲率。多重共线性诊断结果表明，影响因子之间不存在潜在的多重共线性问题，然后计算影响因子的信息增益率，排除了平面曲率因素的影响。
（2）构建了两种原始分类模型，即IOE和SVM模型。同时，计算每个影响因子的分形维数，并将其作为输入数据来构建两种混合模型，即F-IOE和F-SVM模型。最终使用自然断点法将滑坡易发性分区图分为5类：极低易发区、低易发区、中易发区、高易发区和极高易发区。最后，使用3种统计学指标（PPR、NPR和ACC）来评估模型的性能，并使用ROC曲线来比较这4种模型的泛化性。
（3）4种模型的ACC值都高于80%，这表明所有模型的分类结果都可靠，并且两种混合模型的ACC值都高于其原始模型。此外，从基于测试数据集的ROC曲线来看，结果相似，F-IOE模型和F-SVM模型的AUC值分别为
0.8591 和0.9761 ，高于IOE模型（0.7434 ）和SVM模型（0.7946 ）。（4）本研究的结果表明混合模型在滑坡易发性分区中的可靠性高，可以减少野外滑坡详细调查的工作量，为更精准、有效的滑坡灾害预测预警提供数据支撑，同时，相关图件的发布也可为当地的防灾减灾工作提供决策支持。
-
![]()
图 1 塔里木克拉通的位置及前寒武纪岩石分布图 （据Ge et al., 2018修改）
Figure 1. Location and distribution of Precambrian rocks of the Tarim craton
![]()
图 3 敦煌地区太古宙岩石分布图 （据Si et al., 2022修改）
Figure 3. Distribution of Archean rocks in the Dunhuang area
![]()
图 4 北阿尔金地区太古宙岩石分布图 （据Ge et al., 2018, 2020修改）
Figure 4. Distribution of Archean rocks in the North Altyn Tagh area
![]()
图 5 赫罗斯坦杂岩太古宙岩石分布图 （据Ge et al., 2022修改）
Figure 5. Distribution of Archean rocks in the Heluositan Complex
![]()
图 9 热力学–地球化学模拟结果与岩石成分的对比图 （据Ge et al., 2022修改）
Figure 9. Comparison of thermodynamic–geochemical modelling results and rock compositions
![]()
图 6 塔里木克拉通早前寒武纪构造–热事件对比图 （据Ge et al., 2022修改）
Figure 6. Comparison of Early Precambrian tectonothermal events in the Tarim Craton
![]()
图 7 塔里木克拉通早前寒武纪锆石Hf同位素组成与地壳演化图 （据Ge et al., 2022修改）
Figure 7. Early Precambrian zircon Hf isotope composition and crustal evolution of the Tarim craton
![]()
图 8 太古宙TTG成因的热力学–地球化学模拟图 （据Ge et al., 2022修改）
Figure 8. Thermodynamic–geochemical modelling for the formation of TTG melts
![]()
图 10 锆石氧逸度–湿度计的原理和应用 （据Ge et al., 2023修改）
a．锆石氧逸度–湿度计的原理，以Bishop凝灰岩为例，插图展示了水含量计算结果的准确性；b. 锆石氧逸度–湿度计给出的太古宙花岗质岩浆的氧逸度和水含量及其与其他岩浆的对比
Figure 10. Principle and application of the zircon oxybarometer–hygrometer
![]()
图 11 太古宙花岗质岩浆氧逸度和水含量及锆石Hf同位素随时间的长期演化（据Ge et al., 2023修改）
Figure 11. Secular change in magma oxygen fugacity, H2O content and zircon Hf isotopic composition of Archean granitoids
![]()
图 12 太古宙大陆起源的岛弧俯冲–拼贴模式图（据Ge et al., 2022修改）
Figure 12. A subduction–arc accretion model for the origin of continents
-
邓兴梁, 舒良树, 朱文斌, 等 . 新疆兴地断裂带前寒武纪构造-岩浆-变形作用特征及其年龄[J]. 岩石学报,2008 ,24 （12 ）:2800 −2808 .DENG Xingliang, SHU Liangshu, ZHU Wenbin, et al . Precambrian tectonism, deformation and geochronology of igneous rocks in the Xingdi fault zone, Xinjiang[J]. Acta Petrologica Sinica,2008 ,24 （12 ）:2800 −2808 .翟明国, 赵磊, 祝禧艳, 等 . 早期大陆与板块构造启动——前沿热点介绍与展望[J]. 岩石学报,2020 ,36 （8 ）:2249 −2275 . doi: 10.18654/1000-0569/2020.08.01ZHAI Mingguo, ZHAO Lei, ZHU Xiyan, et al . Review and overview for the frontier hotspot: Early continents and start of plate tectonics[J]. Acta Petrologica Sinica,2020 ,36 （8 ）:2249 −2275 . doi: 10.18654/1000-0569/2020.08.01董昕, 张泽明, 唐伟 . 塔里木克拉通北缘的前寒武纪构造热事件——新疆库尔勒铁门关高级变质岩的锆石U-Pb年代学限定[J]. 岩石学报,2011 ,27 （1 ）:47 −58 .DONG Xin, ZHANG Zeming, TANG Wei . Precambrian tectono-thermal events of the northern margin of the Tarim Craton: Constraints of zircon U-Pb chronology from high-grade metamorphic rocks of the Korla, Xinjiang[J]. Acta Petrologica Sinica,2011 ,27 （1 ）:47 −58 .高振家, 陈晋镳, 彭昌文, 等. 新疆北部前寒武系[M]. 北京: 地质出版社, 1993: 171. GAO Zhenjia, CHEN Jinbiao, PENG Changwen., et al. Precambrian Geology of the North Xinjing[M]. Beijing: Geological Press, 1993: 171.
辜平阳, 计文化, 陈锐明, 等 . 塔里木地块东南缘新太古代安南坝石英闪长片麻岩的成因及其对地壳演化的启示[J]. 地球科学,2020 ,45 （9 ）:3268 −3281 .GU Pingyang, JI Wenhua, CHEN Ruiming, et al . Petrogenesis of Neoarchean Ananba Quartz Diorite Gneiss in Southeastern Margin of Tarim: Implications for Crustal Evolution[J]. Earth Science,2020 ,45 （9 ）:3268 −3281 .辜平阳, 徐学义, 何世平, 等 . 塔里木东南缘安南坝地区约 2.5 Ga 花岗闪长质片麻岩的发现及岩石成因[J]. 地质通报,2019 ,38 （5 ）:834 −844 .GU Pingyang, XU Xueyi, HE Shiping, et al . Ca. 2.5 Ga granodioritic gneiss in Annanba area of southeastern Tarim and its petrogenesis[J]. Geological Bulletin of China,2019 ,38 （5 ）:834 −844 .郭新成, 郑玉壮, 高军, 等 . 新疆西昆仑中太古界古陆核的确定及地质意义[J]. 地质论评,2013 ,59 （3 ）:401 −412 . doi: 10.3969/j.issn.0371-5736.2013.03.001GUO Xincheng, ZHENG Yuzhuang, GAO Jun, et al . Determination and Geological Significance of the Mesoarchean Craton in Western Kunlun Mountains, Xinjiang, China[J]. Geological Review,2013 ,59 （3 ）:401 −412 . doi: 10.3969/j.issn.0371-5736.2013.03.001郭瑞清, 秦切, 邹明煜, 等 . 新疆库鲁克塔格西段辉长岩脉年代学、岩石地球化学特征及构造意义[J]. 西北地质,2018 ,51 （4 ）:70 −81 .GUO Ruiqing, QIN Qie, ZOU Mingyu, et al . Geochronology, Petrogeochemical Characteristics and Tectonic Significance of Gabbro Dike from Western Quruqtagh in Xinjiang[J]. Northwestern Geology,2018 ,51 （4 ）:70 −81 .郭召杰, 张志诚, 刘树文, 等 . 塔里木克拉通早前寒武纪基底层序与组合: 颗粒锆石U-Pb年龄新证据[J]. 岩石学报,2003 ,19 （3 ）:537 −542 . doi: 10.3969/j.issn.1000-0569.2003.03.020GUO Zhaojie, ZHANG Zhicheng, LIU Shuwen, et al . U-Pb geochronological evidence for the early Precambrian complex of the Tarim Craton, NW China[J]. Acta Petrologica Sinica,2003 ,19 （3 ）:537 −542 . doi: 10.3969/j.issn.1000-0569.2003.03.020胡霭琴, 罗杰斯 . 新疆塔里木北缘首次发现33亿年的岩石[J]. 科学通报,1992 , （7 ）:627 −630 .HU Aiqin, Rogers G . Discovery of 3.3 Ga Archean Rocks in North Tarim Block of Xinjiang, Western China[J]. Science Bulletin,1992 , （7 ）:627 −630 .胡霭琴, 韦刚健 . 塔里木盆地北缘新太古代辛格尔灰色片麻岩形成时代问题[J]. 地质学报,2006 ,80 （1 ）:126 −134 . doi: 10.3321/j.issn:0001-5717.2006.01.014HU Aiqin, WEI Gangjian . On the Age of the Neo-Archean Qingir Gray Gneisses from the Northern Tarim Basin, Xinjiang, China[J]. Acta Geologica Sinica,2006 ,80 （1 ）:126 −134 . doi: 10.3321/j.issn:0001-5717.2006.01.014黎敦朋, 李新林, 周小康, 等 . 塔里木西南缘新太古代变质辉长岩脉的锆石SHRIMP U-Pb定年及其地质意义[J]. 中国地质,2007 ,34 （2 ）:262 −269 . doi: 10.3969/j.issn.1000-3657.2007.02.007LI Dunpeng, LI Xinlin, ZHOU Xiaokang, et al . SHRIMP U-Pb Zircon Dating of Neoarchean Metagabbro Dikes On the Southwestern Margin of the Tarim Plate and its Significance[J]. Geology in China,2007 ,34 （2 ）:262 −269 . doi: 10.3969/j.issn.1000-3657.2007.02.007李惠民, 陆松年, 郑健康, 等 . 阿尔金山东端花岗片麻岩中3.6 Ga锆石的地质意义[J]. 矿物岩石地球化学通报,2001 ,20 （4 ）:259 −262 . doi: 10.3969/j.issn.1007-2802.2001.04.016LI Huimin, LU Songnian, ZHENG Jiangkang, et al . Dating of 3.6 Ga Zircons in Granite-Gneiss from the Eastern Altyn Mountains and Its Geological Significance[J]. Bulletin of Mineralogy, Petrology and Geochemistry,2001 ,20 （4 ）:259 −262 . doi: 10.3969/j.issn.1007-2802.2001.04.016李晓剑, 王毅, 李慧莉, 等 . 新元古代陆缘岩浆弧——塔里木盆地巴楚隆起的基底: 来自钻井岩芯的最新证据[J]. 岩石学报,2018 ,34 （7 ）:2140 −2164 .LI Xiaojian, WANG Yi, LI Huili, et al . Bachu uplift in the central Tarim Basin based on Neoproterozoic continental arc: New lines of evidence from drilled andesite and dacite[J]. Acta Petrologica Sinica,2018 ,34 （7 ）:2140 −2164 .陆松年 . 新疆库鲁克塔格元古宙地质演化[A]. 中国地质科学院天津地质矿产研究所文集[C].1992 ,26−27 :279 −292 .LU Songnian . Geological evolution of Proterozoic in Kuruktage, Xinjiang[A]. Tianjin Institute of Geology and Mineral Resources, Chinese Academy of Geological Sciences[C].1992 ,26−27 :279 −292 .陆松年, 袁桂邦 . 阿尔金山阿克塔什塔格早前寒武纪岩浆活动的年代学证据[J]. 地质学报,2003 ,77 （1 ）:61 −68 . doi: 10.3321/j.issn:0001-5717.2003.01.008LU Songnian, YUAN Guibang . Geochronology of early Precambrian magmatic activities in Aketashitage, east Altyn Tagh[J]. Acta Geologica Sinica,2003 ,77 （1 ）:61 −68 . doi: 10.3321/j.issn:0001-5717.2003.01.008梅华林, 于海峰, 李铨, 等 . 甘肃北山地区首次发现榴辉岩和古元古花岗质岩石[J]. 科学通报,1998 ,21 （19 ）:2105 −2111 . doi: 10.3321/j.issn:0023-074X.1998.19.022MEI Hualin, YU Haifeng, LI Quan, et al . Archean tonalite in the Dunhuang, Gansu Provience: age from the U-Pb sigle zircon and Nd isotope[J]. Progress in Precambrian Research,1998 ,21 （19 ）:2105 −2111 . doi: 10.3321/j.issn:0023-074X.1998.19.022邬光辉, 张承泽, 汪海, 等 . 塔里木盆地中部塔参1井花岗闪长岩的锆石SHRIMP U-Pb年龄[J]. 地质通报,2009 ,28 （5 ）:568 −571 .WU Guanghui, ZHANG Chengze, WANG Hai, et al . Zircon SHRIMP U-Pb age of granodiorite of the Tacan 1 well in the cen-tral Tarim basin, China[J]. Geologcal Bulletin of China,2009 ,28 （5 ）:568 −571 .辛后田, 刘永顺, 罗照华, 等 . 塔里木盆地东南缘阿克塔什塔格地区新太古代陆壳增生: 米兰岩群和 TTG 片麻岩的地球化学及年代学约束[J]. 地学前缘,2013 ,20 （1 ）:240 −259 .XIN Houtian, LIU Yongshun, LUO Zhaohua, et al . The growth of Archean continental crust in Aqtashtagh Area of Southeast Tarim, China: Constraints from petrochemistry and chronology about Milan Group and TTG-gneiss[J]. Earth Science Frontiers,2013 ,20 （1 ）:240 −259 .辛后田, 罗照华, 刘永顺, 等 . 塔里木东南缘阿克塔什塔格地区古元古代壳源碳酸岩的特征及其地质意义[J]. 地学前缘,2012 ,19 （6 ）:167 −178 .XIN Houtian, LUO Zhaohua, LIU Yongshun, et al . Geological features and significance of Paleoproterozoic carbonatite of crust origin in Aqtashtagh area of southeast Tarim Basin, China[J]. Earth Science Frontiers,2012 ,19 （6 ）:167 −178 .辛后田, 赵凤清, 罗照华, 等 . 塔里木盆地东南缘阿克塔什塔格地区古元古代精细年代格架的建立及其地质意义[J]. 地质学报,2011 ,85 （12 ）:1977 −1993 .XIN Houtian, ZHAO Fengqing, LUO Zhaohua, et al . Determination of the Paleoproterozoic Geochronological Framework in Aqtashtagh Area in Southeastern Tarim, China, and Its Geological Significance[J]. Acta Geologica Sinica,2011 ,85 （12 ）:1977 −1993 .新疆维吾尔自治区地质矿产局. 新疆维吾尔自治区区域地质志[M]. 北京: 地质出版社, 1993: 762. 叶现韬, 张传林 . 阿尔金北缘新太古代 TTG 片麻岩的成因及其构造意义[J]. 岩石学报,2020 ,36 （11 ）:3397 −3413 . doi: 10.18654/1000-0569/2020.11.09YE Xiantao, ZHANG Chuanlin . Petrogenesis and tectonic implications of the Neoarchean TTG gneiss in the North Altyn Tagh area, southeastern Tarim Craton[J]. Acta Petrologica Sinica,2020 ,36 （11 ）:3397 −3413 . doi: 10.18654/1000-0569/2020.11.09张传林, 王中刚, 沈加林, 等 . 西昆仑山阿卡孜岩体锆石SHRIMP定年及其地球化学特征[J]. 岩石学报,2003 ,19 （3 ）:523 −529 . doi: 10.3969/j.issn.1000-0569.2003.03.018ZHANG Chuanlin, WANG Zhonggang, SHEN Jialin, et al . Zircon SHRIMP dating and geochemistry characteristics of Akazi rock mass of Western Kunlun[J]. Acta Petrologica Sinica,2003 ,19 （3 ）:523 −529 . doi: 10.3969/j.issn.1000-0569.2003.03.018赵国春, 张国伟 . 大陆的起源[J]. 地质学报,2021 ,95 （1 ）:1 −19 . doi: 10.1111/1755-6724.14621ZHAO Guochun, ZHANG Guowei . Origin of continents[J]. Acta Geologica Sinica,2021 ,95 （1 ）:1 −19 . doi: 10.1111/1755-6724.14621赵燕, 第五春荣, 敖文昊, 等 . 敦煌地块发现～3.06 Ga花岗闪长质片麻岩[J]. 科学通报,2015 ,60 （1 ）:75 −87 .ZHAO Yan, DIWU Chunrong, AO Wenhao, et al . Ca. 3.06 Ga granodioritic gneiss in Dunhuang block[J]. 科学通报,2015 ,60 （1 ）:75 −87 .赵燕, 第五春荣, 孙勇, 等 . 甘肃敦煌水峡口地区前寒武纪岩石的锆石 U-Pb 年龄, Hf 同位素组成及其地质意义[J]. 岩石学报,2013 ,29 （5 ）:1698 −1712 .ZHAO Yan, DIWU Chunrong, SUN Yong, et al . Zircon geochronology and Lu-Hf isotope compositions for Precambrian rocks of the Dunhuang complex in Shuixiakou area, Gansu Province[J]. Acta Petrologica Sinica,2013 ,29 （5 ）:1698 −1712 .郑永飞 . 太古宙地质与板块构造: 观察与解释[J]. 中国科学: 地球科学,2024 ,54 （1 ）:1 −30 .ZHENG Yongfei . Plate tectonics in the Archean: Observations versus interpretations[J]. Science China Earth Sciences,2024 ,54 （1 ）:1 −30 .朱文斌, 林和丰, 葛荣峰, 等 . 塔里木克拉通北缘库鲁克塔格地块太古宙基底组成与地壳演化[J]. 地质学报,2022 ,96 （9 ）:3084 −3101 .ZHU Wenbin, LIN Hefeng, GE Rongfeng, et al . Archean basement composition and crustal evolution of the Kuluketage block in the northern margin of the Tarim Craton[J]. Acta Geologica Sinica,2022 ,96 （9 ）:3084 −3101 .Aarons S M, Johnson A C, Rader S T . Forming Earth's Continental Crust: A Nontraditional Stable Isotope Perspective[J]. Elements,2021 ,17 （6 ）:413 −418 . doi: 10.2138/gselements.17.6.413Aarons S M, Reimink J R, Greber N D, et al . Titanium isotopes constrain a magmatic transition at the Hadean-Archean boundary in the Acasta Gneiss Complex[J]. Science Advances,2020 ,6 （50 ）:eabc9959 . doi: 10.1126/sciadv.abc9959Armstrong R L . A model for the evolution of strontium and lead isotopes in a dynamic earth[J]. Reviews of Geophysics,1968 ,6 （2 ）:175 −199 . doi: 10.1029/RG006i002p00175Armstrong R L . The Persistent Myth of Crustal Growth[J]. Australian Journal of Earth Sciences,1991 ,38 （5 ）:613 −630 . doi: 10.1080/08120099108727995Arndt N T . The Formation and Evolution of the Continental Crust[J]. Geochemical Perspectives,2013 ,3 （2 ）:405 −533 .Atherton M P, Petford N . Generation of sodium-rich magmas from newly underplated basaltic crust[J]. Nature,1993 ,362 （6416 ）:144 −146 . doi: 10.1038/362144a0Aulbach S, Stagno V . Evidence for a reducing Archean ambient mantle and its effects on the carbon cycle[J]. Geology,2016 ,44 （9 ）:751 −754 . doi: 10.1130/G38070.1Barker F, Arth J G . Generation of trondhjemitic-tonalitic liquids and Archean bimodal trondhjemite-basalt suites[J]. Geology,1976 ,4 （10 ）:596 −600 . doi: 10.1130/0091-7613(1976)4<596:GOTLAA>2.0.CO;2Bédard J H . A catalytic delamination-driven model for coupled genesis of Archaean crust and sub-continental lithospheric mantle[J]. Geochimica et Cosmochimica Acta,2006 ,70 （5 ）:1188 −1214 . doi: 10.1016/j.gca.2005.11.008Bédard J H . Stagnant lids and mantle overturns: Implications for Archaean tectonics, magmagenesis, crustal growth, mantle evolution, and the start of plate tectonics[J]. Geoscience Frontiers,2018 ,9 （1 ）:19 −49 . doi: 10.1016/j.gsf.2017.01.005Belousova E A, Kostitsyn Y A, Griffin W L, et al . The growth of the continental crust: Constraints from zircon Hf-isotope data[J]. Lithos,2010 ,119 （3−4 ）:457 −466 . doi: 10.1016/j.lithos.2010.07.024Cai Z, Jiao C, He B, et al. Archean–Paleoproterozoic tectonothermal events in the central Tarim Block: constraints from granitic gneisses revealed by deep drilling wells [J]. Precambrian Research, 2020: 105776.
Cai Z, Xu Z, Yu S, et al . Neoarchean magmatism and implications for crustal growth and evolution of the Kuluketage region, northeastern Tarim Craton[J]. Precambrian Research,2018 ,304 :156 −170 . doi: 10.1016/j.precamres.2017.11.016Cawood P A, Chowdhury P, Mulder J A, et al. Secular evolution of continents and the Earth system [J]. Reviews of Geophysics, 2023: e2022RG000789.
Collins W J, Murphy J B, Johnson T E, et al . Critical role of water in the formation of continental crust[J]. Nature Geoscience,2020 ,13 :331 −338 . doi: 10.1038/s41561-020-0573-6Condie K C . Episodic Continental Growth and Supercontinents: A Mantle Avalanche Connection?[J]. Earth and Planetary Science Letters,1998 ,163 （1−4 ）:97 −108 . doi: 10.1016/S0012-821X(98)00178-2Condie K C. How to make a continent: thirty-five years of TTG research. In Dilek Y, Furnes H (eds). Evolution of Archean Crust and Early Life[M]. Springer, 2014: 179−193.
Condie K C, Kröner A . The building blocks of continental crust: Evidence for a major change in the tectonic setting of continental growth at the end of the Archean[J]. Gondwana Research,2013 ,23 （2 ）:394 −402 . doi: 10.1016/j.gr.2011.09.011Dhuime B, Hawkesworth C J, Cawood P A, et al . A change in the geodynamics of continental growth 3 billion years ago[J]. Science,2012 ,335 （6074 ）:1334 −1336 . doi: 10.1126/science.1216066Dong C, Ge R, Liu S, et al . Multiple episodes of early Precambrian magmatism and tectonism in the Tarim Craton: A North China connection[J]. Lithos,2022 ,430 :106883 .Drummond M S, Defant M J . A model for trondhjemite-tonalite-dacite genesis and crustal growth via slab melting - Archean to modern comparisons[J]. Journal of Geophysical Research-Solid Earth and Planets,1990 ,95 （B13 ）:21503 −21521 . doi: 10.1029/JB095iB13p21503Feng L, Lin S, Davis D W, et al . Dunhuang Tectonic Belt in northwestern China as a part of the Central Asian Orogenic Belt: Structural and U-Pb geochronological evidence[J]. Tectonophysics,2018 ,747−748 :281 −297 . doi: 10.1016/j.tecto.2018.09.008Fischer R, Gerya T . Regimes of subduction and lithospheric dynamics in the Precambrian: 3D thermomechanical modelling[J]. Gondwana Research,2016 ,37 :53 −70 . doi: 10.1016/j.gr.2016.06.002Ganne J, Feng X . Primary magmas and mantle temperatures through time[J]. Geochemistry, Geophysics, Geosystems,2017 ,18 （3 ）:872 −888 . doi: 10.1002/2016GC006787Gao L, Liu S, Cawood P A, et al . Oxidation of Archean upper mantle caused by crustal recycling[J]. Nature Communications,2022 ,13 （1 ）:3283 . doi: 10.1038/s41467-022-30886-4Ge R F, Zhu W B, Wilde S A, et al . Archean magmatism and crustal evolution in the northern Tarim Craton: Insights from zircon U–Pb–Hf–O isotopes and geochemistry of ～2.7 Ga orthogneiss and amphibolite in the Korla Complex[J]. Precambrian Research,2014a ,252 :145 −165 . doi: 10.1016/j.precamres.2014.07.019Ge R F, Zhu W B, Wilde S A, et al . Zircon U–Pb–Lu–Hf–O isotopic evidence for ≥3.5 Ga crustal growth, reworking and differentiation in the northern Tarim Craton[J]. Precambrian Research,2014b ,249 :115 −128 . doi: 10.1016/j.precamres.2014.05.004Ge R F, Zhu W B, Wilde S A, et al . Synchronous crustal growth and reworking recorded in late Paleoproterozoic granitoids in the northern Tarim craton: In situ zircon U-Pb-Hf-O isotopic and geochemical constraints and tectonic implications[J]. Geological Society of America Bulletin,2015 ,127 （5−6 ）:781 −803 . doi: 10.1130/B31050.1Ge R F, Zhu W B, Wu H L, et al . Timing and mechanisms of multiple episodes of migmatization in the Korla Complex, northern Tarim Craton, NW China: Constraints from zircon U–Pb–Lu–Hf isotopes and implications for crustal growth[J]. Precambrian Research,2013 ,231 :136 −156 . doi: 10.1016/j.precamres.2013.03.005Ge R, Wilde S A, Kemp A I S, et al . Generation of Eoarchean continental crust from altered mafic rocks derived from a chondritic mantle: The ∼3.72 Ga Aktash gneisses, Tarim Craton (NW China)[J]. Earth and Planetary Science Letters,2020 ,538 :116225 . doi: 10.1016/j.jpgl.2020.116225Ge R, Wilde S A, Zhu W, et al . Formation and evolution of Archean continental crust: A thermodynamic–geochemical perspective of granitoids from the Tarim Craton, NW China[J]. Earth-Science Reviews,2022 ,234 :104219 . doi: 10.1016/j.earscirev.2022.104219Ge R, Wilde S A, Zhu W, et al . Earth’s early continental crust formed from wet and oxidizing arc magmas[J]. Nature,2023 ,623 :334 −339 . doi: 10.1038/s41586-023-06552-0Ge R, Zhu W, Wilde S A, et al . Remnants of Eoarchean continental crust derived from a subducted proto-arc[J]. Science Advances,2018 ,4 （2 ）:aao3159 . doi: 10.1126/sciadv.aao3159Gehrels G E, Yin A, Wang X . Magmatic history of the northeastern Tibetan Plateau[J]. Journal of Geophysical Research,2003 ,108 （B9 ）:2423 .Green E, White R W, Diener J, et al . Activity‐composition relations for the calculation of partial melting equilibria in metabasic rocks[J]. Journal of Metamorphic Geology,2016 ,34 （9 ）:845 −869 . doi: 10.1111/jmg.12211Guo Z J, Yin A, Robinson A, et al . Geochronology and Geochemistry of Deep-Drill-Core Samples From the Basement of the Central Tarim Basin[J]. Journal of Asian Earth Sciences,2005 ,25 （1 ）:45 −56 . doi: 10.1016/j.jseaes.2004.01.016Harrison T M. Hadean Earth[M]. Springer, 2020:1−291.
Hastie A R, Fitton J G, Bromiley G D, et al . The origin of Earth's first continents and the onset of plate tectonics[J]. Geology,2016 ,44 （10 ）:855 −858 . doi: 10.1130/G38226.1Hawkesworth C J, Kemp A I S . Evolution of the continental crust[J]. Nature,2006 ,443 （7113 ）:811 −817 . doi: 10.1038/nature05191Herzberg C, Condie K, Korenaga J . Thermal history of the Earth and its petrological expression[J]. Earth and Planetary Science Letters,2010 ,292 （1 ）:79 −88 .Hoffmann J E, Zhang C, Nagel T. The formation of Tonalites–Trondjhemite–Granodiorites in Early continental crust. In Earth's oldest rocks, Van Kranendonk M J, Bennett V, Hoffmann E, Eds.; Elsevier: 2019: 133−168.
Holland T J B, Powell R . An improved and extended internally consistent thermodynamic dataset for phases of petrological interest, involving a new equation of state for solids[J]. Journal of Metamorphic Geology,2011 ,29 （3 ）:333 −383 . doi: 10.1111/j.1525-1314.2010.00923.xHolland T J, Green E C, Powell R . Melting of peridotites through to granites: a simple thermodynamic model in the system KNCFMASHTOCr[J]. Journal of Petrology,2018 ,59 （5 ）:881 −900 . doi: 10.1093/petrology/egy048Jagoutz O, Schmidt M W, Enggist A, et al . TTG-type plutonic rocks formed in a modern arc batholith by hydrous fractionation in the lower arc crust[J]. Contributions to Mineralogy and Petrology,2013 ,166 （4 ）:1099 −1118 . doi: 10.1007/s00410-013-0911-4Jahn B, Glikson A Y, Peucat J J, et al . REE geochemistry and isotopic data of Archean silicic volcanics and granitoids from the Pilbara Block, Western Australia: implications for the early crustal evolution[J]. Geochimica et Cosmochimica Acta,1981 ,45 （9 ）:1633 −1652 . doi: 10.1016/S0016-7037(81)80002-6Johnson T E, Brown M, Gardiner N J, et al . Earth's first stable continents did not form by subduction[J]. Nature,2017 ,543 :239 −242 . doi: 10.1038/nature21383Johnson T E, Kirkland C L, Lu Y, et al . Giant impacts and the origin and evolution of continents[J]. Nature,2022 ,608 （7922 ）:330 −335 . doi: 10.1038/s41586-022-04956-yKamber B S . The evolving nature of terrestrial crust from the Hadean, through the Archaean, into the Proterozoic[J]. Precambrian Research,2015 ,258 :48 −82 . doi: 10.1016/j.precamres.2014.12.007Kemp A I S, Hawkesworth C J, Paterson B A, et al . Episodic growth of the Gondwana supercontinent from hafnium and oxygen isotopes in zircon[J]. Nature,2006 ,439 （7076 ）:580 −583 . doi: 10.1038/nature04505Kendrick J, Duguet M, Yakymchuk C . Diversification of Archean tonalite - trondhjemite - granodiorite suites in a mushy middle crust[J]. Geology,2022 ,50 （1 ）:76 −80 . doi: 10.1130/G49287.1Kleinhanns I C, Kramers J D, Kamber B S . Importance of water for Archaean granitoid petrology: a comparative study of TTG and potassic granitoids from Barberton Mountain Land, South Africa[J]. Contributions to Mineralogy and Petrology,2003 ,145 （3 ）:377 −389 . doi: 10.1007/s00410-003-0459-9Korenaga J . Hadean geodynamics and the nature of early continental crust[J]. Precambrian Research,2021 ,359 :106178 . doi: 10.1016/j.precamres.2021.106178Laurent O, Björnsen J, Wotzlaw J, et al . Earth's earliest granitoids are crystal-rich magma reservoirs tapped by silicic eruptions[J]. Nature Geoscience,2020 ,13 （2 ）:163 −169 . doi: 10.1038/s41561-019-0520-6Lee C A, Luffi P, Chin E J . Building and Destroying Continental Mantle[J]. Annual Review of Earth and Planetary Sciences,2011 ,39 （1 ）:59 −90 . doi: 10.1146/annurev-earth-040610-133505Li W, Costa F . A thermodynamic model for F-Cl-OH partitioning between silicate melts and apatite including non-ideal mixing with application to constraining melt volatile budgets[J]. Geochimica et Cosmochimica Acta,2020 ,269 :203 −222 . doi: 10.1016/j.gca.2019.10.035Liou P, Guo J . Generation of Archaean TTG Gneisses Through Amphibole-Dominated Fractionation[J]. Journal of Geophysical Research: Solid Earth,2019 ,124 （4 ）:3605 −3619 . doi: 10.1029/2018JB017024Long X P, Yuan C, Sun M, et al . Archean crustal evolution of the northern Tarim Craton, NW China: zircon U-Pb and Hf isotopic constraints[J]. Precambrian Research,2010 ,180 （3−4 ）:272 −284 . doi: 10.1016/j.precamres.2010.05.001Long X P, Yuan C, Sun M, et al . The discovery of the oldest rocks in the Kuluketage area and its geological implications[J]. Science in China Series D: Earth Sciences,2011 ,54 （3 ）:342 −348 . doi: 10.1007/s11430-010-4156-zLong X P, Yuan C, Sun M, et al. New geochemical and combined zircon U–Pb and Lu–Hf isotopic data of orthogneisses in the northern Altyn Tagh, northern margin of the Tibetan plateau: Implication for Archean evolution of the Dunhuang Block and crust formation in NW China [J]. Lithos, 2014, 200–201: 418 - 431.
Long X, Wilde S A, Yuan C, et al . Provenance and depositional age of Paleoproterozoic metasedimentary rocks in the Kuluketage Block, northern Tarim Craton: Implications for tectonic setting and crustal growth[J]. Precambrian Research,2015 ,260 :76 −90 . doi: 10.1016/j.precamres.2015.01.008Loucks R R, Fiorentini M L, Henríquez G J . New magmatic oxybarometer using trace elements in zircon[J]. Journal of Petrology,2020 ,61 （3 ）:egaa34 .Lu S N, Li H K, Zhang C L, et al . Geological and geochronological evidence for the Precambrian evolution of the Tarim Craton and surrounding continental fragments[J]. Precambrian Research,2008 ,160 （1−2 ）:94 −107 . doi: 10.1016/j.precamres.2007.04.025Martin H . Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas[J]. Geology,1986 ,14 （9 ）:753 −756 . doi: 10.1130/0091-7613(1986)14<753:EOSAGG>2.0.CO;2Martin H, Moyen J, Guitreau M, et al . Why Archaean TTG cannot be generated by MORB melting in subduction zones[J]. Lithos,2014 ,198−199 :1 −13 .Moore W B, Webb A A G . Heat-pipe earth[J]. Nature,2013 ,501 （7468 ）:501 −505 . doi: 10.1038/nature12473Moyen J F, Stevens G . Experimental constraints on TTG petrogenesis: implications for Archean geodynamics[J]. Geophysical Monograph-American Geophysical Union,2006 ,164 :149 −175 .Moyen J . The composite Archaean grey gneisses: Petrological significance, and evidence for a non-unique tectonic setting for Archaean crustal growth[J]. Lithos,2011 ,123 （1−4 ）:21 −36 . doi: 10.1016/j.lithos.2010.09.015Moyen J . Archean granitoids: classification, petrology, geochemistry and origin[J]. Geological Society, London, Special Publications,2020 ,489 （1 ）:15 −49 . doi: 10.1144/SP489-2018-34Moyen J, Martin H . Forty years of TTG research[J]. Lithos,2012 ,148 :312 −336 . doi: 10.1016/j.lithos.2012.06.010Palin R M, White R W, Green E C R . Partial melting of metabasic rocks and the generation of tonalitic–trondhjemitic–granodioritic (TTG) crust in the Archaean: constraints from phase equilibrium modelling[J]. Precambrian Research,2016 ,287 :73 −90 . doi: 10.1016/j.precamres.2016.11.001Pearson D G, Scott J M, Liu J, et al . Deep continental roots and cratons[J]. Nature,2021 ,596 （7871 ）:199 −210 . doi: 10.1038/s41586-021-03600-5Pourteau A, Doucet L S, Blereau E R, et al . TTG generation by fluid-fluxed crustal melting: Direct evidence from the Proterozoic Georgetown Inlier, NE Australia[J]. Earth and Planetary Science Letters,2020 ,550 :116548 . doi: 10.1016/j.jpgl.2020.116548Rapp R P, Shimizu N, Norman M D . Growth of early continental crust by partial melting of eclogite[J]. Nature,2003 ,425 （6958 ）:605 −609 . doi: 10.1038/nature02031Rollinson H . Do all Archaean TTG rock compositions represent former melts?[J]. Precambrian Research,2021 ,367 :106448 . doi: 10.1016/j.precamres.2021.106448Roman A, Arndt N . Differentiated Archean oceanic crust: Its thermal structure, mechanical stability and a test of the sagduction hypothesis[J]. Geochimica et Cosmochimica Acta,2020 ,278 :65 −77 . doi: 10.1016/j.gca.2019.07.009Rozel A B, Golabek G J, Jain C, et al . Continental crust formation on early Earth controlled by intrusive magmatism[J]. Nature,2017 ,545 （7654 ）:332 −335 . doi: 10.1038/nature22042Rudnick R L . Making continental crust[J]. Nature,1995 ,378 （6557 ）:571 −578 . doi: 10.1038/378571a0Rudnick R L, Gao S. Composition of the continental crust[A]. In: Rudnick R L (ed). Treatise on geochemistry[M]. Elsevier, 2014, 4: 1−51.
Shi M, Hou Q, Wu C, et al . Paleozoic Sanweishan arc in the northern Dunhuang region, NW China: The Dunhuang block is a Phanerozoic orogen, not a Precambrian block[J]. Journal of Asian Earth Sciences,2020 ,194 :103954 . doi: 10.1016/j.jseaes.2019.103954Shu L S, Deng X L, Zhu W B, et al . Precambrian tectonic evolution of the Tarim Block, NW China: New geochronological insights from the Quruqtagh domain[J]. Journal of Asian Earth Sciences,2011 ,42 （5 ）:774 −790 . doi: 10.1016/j.jseaes.2010.08.018Si Y, Ge R, Zhou T, et al. Decoupling of metamorphic zircon U-Pb ages and P-T paths in the Dunhuang metamorphic complex, northwestern China [J]. Precambrian Research, 2022, 379: 106783.
Sizova E, Gerya T, Brown M, et al . Subduction styles in the Precambrian: Insight from numerical experiments[J]. Lithos,2010 ,116 （3−4 ）:209 −229 . doi: 10.1016/j.lithos.2009.05.028Sizova E, Gerya T, Stüwe K, et al . Generation of felsic crust in the Archean: A geodynamic modeling perspective[J]. Precambrian Research,2015 ,271 :198 −224 . doi: 10.1016/j.precamres.2015.10.005Smithies R H, Lu Y, Johnson T E, et al . No evidence for high-pressure melting of Earth's crust in the Archean[J]. Nature Communications,2019 ,10 （1 ）:1 −12 . doi: 10.1038/s41467-018-07882-8Smithies R H, Lu Y, Kirkland C L, et al . Oxygen isotopes trace the origins of Earth’s earliest continental crust[J]. Nature,2021 ,592 （7852 ）:70 −75 . doi: 10.1038/s41586-021-03337-1Smythe D J, Brenan J M . Magmatic oxygen fugacity estimated using zircon-melt partitioning of cerium[J]. Earth and Planetary Science Letters,2016 ,453 :260 −266 . doi: 10.1016/j.jpgl.2016.08.013Sobolev A V, Asafov E V, Gurenko A A, et al . Komatiites reveal a hydrous Archaean deep-mantle reservoir[J]. Nature,2016 ,531 （7596 ）:628 −632 . doi: 10.1038/nature17152Sobolev A V, Asafov E V, Gurenko A A, et al . Deep hydrous mantle reservoir provides evidence for crustal recycling before 3.3 billion years ago[J]. Nature,2019 ,571 （7766 ）:555 −559 . doi: 10.1038/s41586-019-1399-5Sun X, Li X, Lei R, et al . Paleoproterozoic crustal evolution of the Tarim Craton, NW China: Constraints from geochronology and geochemistry of orthogneisses and granitic veins in the Xingdi region of the Quruqtagh Block[J]. Precambrian Research,2023 ,399 :107247 . doi: 10.1016/j.precamres.2023.107247Tang M, Wang X, Shu X, et al . Hafnium isotopic heterogeneity in zircons from granitic rocks: Geochemical evaluation and modeling of “zircon effect” in crustal anatexis[J]. Earth and Planetary Science Letters,2014 ,389 :188 −199 . doi: 10.1016/j.jpgl.2013.12.036Taylor S R, McLennan S M . The geochemical evolution of the continental crust[J]. Reviews of Geophysics,1995 ,33 （2 ）:241 −265 . doi: 10.1029/95RG00262Taylor S R, McLennan S. Planetary crusts: their composition, origin and evolution[M]. Cambridge: Cambridge University Press, 2009, 378.
Van Hunen J, Moyen J . Archean Subduction: Fact or Fiction?[J]. Annual Review of Earth and Planetary Sciences,2012 ,40 （1 ）:195 −219 . doi: 10.1146/annurev-earth-042711-105255Van Hunen J, van den Berg A P . Plate tectonics on the early Earth: Limitations imposed by strength and buoyancy of subducted lithosphere[J]. Lithos,2008 ,103 （1−2 ）:217 −235 . doi: 10.1016/j.lithos.2007.09.016Van Kranendonk M J, Bennett V C, Hoffmann J E. Earth's oldest rocks[M]. Elsevier, 2019, 1078.
Vervoort J D, Kemp A I S . Clarifying the zircon Hf isotope record of crust–mantle evolution[J]. Chemical Geology,2016 ,425 :65 −75 . doi: 10.1016/j.chemgeo.2016.01.023Wang C, Liu L, Wang Y, et al . Recognition and tectonic implications of an extensive Neoproterozoic volcano-sedimentary rift basin along the southwestern margin of the Tarim Craton, northwestern China[J]. Precambrian Research,2015 ,257 :65 −82 . doi: 10.1016/j.precamres.2014.11.022Wang C, Wang Y, Liu L, et al . The Paleoproterozoic magmatic–metamorphic events and cover sediments of the Tiekelik Belt and their tectonic implications for the southern margin of the Tarim Craton, northwestern China[J]. Precambrian Research,2014 ,254 :210 −225 . doi: 10.1016/j.precamres.2014.08.018Wang H Y, Chen H, Zhang Q W, et al . Tectonic mélange records the Silurian–Devonian subduction-metamorphic process of the southern Dunhuang terrane, southernmost Central Asian Orogenic Belt[J]. Geology,2017 ,45 （5 ）:427 −430 . doi: 10.1130/G38834.1Wu Z, Song J, Zhao G, et al . Water-Induced Mantle Overturns Leading to the Origins of Archean Continents and Subcontinental Lithospheric Mantle[J]. Geophysical Research Letters,2023 ,50 （22 ）:e2023GL105178 . doi: 10.1029/2023GL105178Xia X, Cui Z, Li W, et al. Zircon water content: Reference material development and simultaneous measurement with oxygen isotope by SIMS [J]. Journal of Analytical Atomic Spectrometry, 2019.
Xiang H, Connolly J A D . GeoPS: an interactive visual computing tool for thermodynamic modeling of phase equilibria[J]. Journal of Metamorphic Geology,2022 ,40 :243 −255 . doi: 10.1111/jmg.12626Xu Z Q, He B Z, Zhang C L, et al . Tectonic framework and crustal evolution of the Precambrian basement of the Tarim Block in NW China: new geochronological evidence from deep drilling samples[J]. Precambrian Research,2013 ,235 :150 −162 . doi: 10.1016/j.precamres.2013.06.001Yang H, Wu G, Kusky T M, et al . Paleoproterozoic assembly of the North and South Tarim terranes: New insights from deep seismic profiles and Precambrian granite cores[J]. Precambrian Research,2018 ,305 :151 −165 . doi: 10.1016/j.precamres.2017.11.015Ye X, Zhang C, Santosh M, et al . Growth and evolution of Precambrian continental crust in the southwestern Tarim terrane: New evidence from the ca. 1.4 Ga A-type granites and Paleoproterozoic intrusive complex[J]. Precambrian Research,2016 ,275 :18 −34 . doi: 10.1016/j.precamres.2015.12.017Zhang C L, Li H K, Santosh M, et al . Precambrian evolution and cratonization of the Tarim Block, NW China: Petrology, geochemistry, Nd-isotopes and U–Pb zircon geochronology from Archaean gabbro-TTG–potassic granite suite and Paleoproterozoic metamorphic belt[J]. Journal of Asian Earth Sciences,2012 ,47 :5 −20 . doi: 10.1016/j.jseaes.2011.05.018Zhang C L, Li Z X, Li X H, et al . An early Paleoproterozoic high-K intrusive complex in southwestern Tarim block, NW China: Age, geochemistry, and tectonic implications[J]. Gondwana Research,2007 ,12 （1−2 ）:101 −112 . doi: 10.1016/j.gr.2006.10.006Zhang C, Ye X, Zou H, et al . Neoproterozoic sedimentary basin evolution in southwestern Tarim, NW China: New evidence from field observations, detrital zircon U–Pb ages and Hf isotope compositions[J]. Precambrian Research,2016 ,280 :31 −45 . doi: 10.1016/j.precamres.2016.04.011Zhang C, Zou H, Santosh M, et al . Is the Precambrian basement of the Tarim Craton in NW China composed of discrete terranes?[J]. Precambrian Research,2014 ,254 :226 −244 . doi: 10.1016/j.precamres.2014.08.006Zhang J X, Yu S Y, Gong J H, et al . The latest Neoarchean–Paleoproterozoic evolution of the Dunhuang block, eastern Tarim craton, northwestern China: Evidence from zircon U–Pb dating and Hf isotopic analyses[J]. Precambrian Research,2013 ,226 :21 −42 . doi: 10.1016/j.precamres.2012.11.014Zhang Q, Zhao L, Zhou D, et al . No evidence of supracrustal recycling in Si-O isotopes of Earth’s oldest rocks 4 Ga ago[J]. Science Advances,2023 ,9 （26 ）:eadf693 .Zhao Y, Sun Y, Diwu C, et al . The Dunhuang block is a Paleozoic orogenic belt and part of the Central Asian Orogenic Belt (CAOB), NW China[J]. Gondwana Research,2016 ,30 :207 −223 . doi: 10.1016/j.gr.2015.08.012Zhao Y, Sun Y, Yan J, et al . The Archean-Paleoproterozoic crustal evolution in the Dunhuang region, NW China: Constraints from zircon U–Pb geochronology and in situ Hf isotopes[J]. Precambrian Research,2015 ,271 :83 −97 . doi: 10.1016/j.precamres.2015.10.002Zhu R, Zhao G, Xiao W, et al. Origin, Accretion and Reworking of Continents [J]. Reviews of Geophysics, 2021: e2019RG000689.
Zong K Q, Liu Y S, Zhang Z M, et al . The generation and evolution of Archean continental crust in the Dunhuang block, northeastern Tarim craton, northwestern China[J]. Precambrian Research,2013 ,235 :251 −263 . doi: 10.1016/j.precamres.2013.07.002 -
期刊类型引用(3)
1. 李光明，杨玉飞，唐亚明，王小浩，尹春旺，冯凡，周永恒. 数据驱动模型评价滑坡易发性的对比研究：以黄河中游流域为例. 西北地质. 2025(02): 51-65 . 
本站查看
2. 王新刚，王道正，王家鼎，黄强兵，胡胜，连宝琴，辜超颖. 黄土高原泥流灾变机理研究进展与发展趋势. 西北地质. 2025(02): 1-15+345 . 本站查看
3. 王家政，辛鹏，曹生鸿，李先臣，郭丹丹. 降雨诱发黄土滑坡–泥流的规律及阈值曲线研究. 西北地质. 2025(02): 102-110 . 本站查看
其他类型引用(0)
下载:
计量
- 文章访问数: 351
- HTML全文浏览量: 44
- PDF下载量: 151
- 被引次数: 3
