-
岩石圈深部探测与研究已经成为固体地球科学发展的前沿之一，而了解深部物质组成、分布和时空演化是探索地球深部组成及动力学的关键。以往这方面的研究多集中于局部和点上的物质组成与循环，而对深部物质组成时空分布等的研究偏少，影响到我们对区域、造山带乃至整个地球尺度的深部物质架构、演变及动力学的完整认识，这主要与缺乏有效深部物质架构的探测技术有关。长期以来，深部探测主要依赖地球物理探测和深部钻探技术（吴功建等，1998; 刘光鼎，2002; 滕吉文，2004; 杨文采等，2004; 金振民等，2004; 许志琴等，2005; 董树文等，2012，2013，2014; Dong Shuwen et al.，2013; Gao Rui et al.，2016; 谢和平等，2017）。钻探可以直接观察深部物质等，但是钻探深度有限，目前最深也只有1万多米，对全面了解地球深部相差甚远。另外，深钻工程成本高，且只是“一孔之见”，资料积累十分有限，难以开展区域构造尺度详细探测。地球物理是最常见的深部探测技术，可以非常好地探测现今地球深部的结构构造，但是难以直接揭示深部物质的属性——特别是随时间演化特征，只能结合其他方法共同约束深部物质的某些特征。这类方法的主要依据为:假设地表出露的来自深部的某些岩石（如深部变质岩、捕虏体等），获取其物理参数（包括温压条件和年龄信息），再应用地球物理技术，推测深部相应物性特征的范围，推断深部物质特征，如常应用于不同比例尺的三维地质填图的物探技术和岩石圈深部探测技术。目前，尚缺乏不同区域尺度的深部物质探测技术。
-
岩浆岩及其携带的深源岩石包体（Enclave）或捕虏体（Xenolith）是来自深部的“使者”，以其作为“岩石探针”（莫宣学等，2011），可以打开探索深部岩石圈岩浆作用和物质特征信息的窗口（邓晋福，1994; Griffin et al.，1998; Condie，1999; 张宏福等，2007; 莫宣学等，2011，2019）。
-
20世纪80~90年代，国内外很多学者对其进行了卓有成效的探索，并被广泛应用于阐述岩石成因与区域演化方面（郑建平等，2021）。近年来，随着分析测试技术的不断发展和提高，发表的岩浆岩相关岩石的测试数据急剧增长。通过岩石学与同位素年代学、地球化学、矿床学等传统学科的紧密交叉融合，岩石探针方面的研究也迅猛发展。目前，国际上以澳大利亚“地球化学演化与大陆成矿”国家研究中心（GEMOC）、“核-壳流体系统”国家卓越中心（CCFS）的研究成果为代表。这些研究在地幔岩石地球化学、岩石圈四维填图和深部岩石圈演化等方面取得了系统性的成果（Griffin et al.，2004; Begg et al.，2009; Belousova et al.，2009; Xu Bo et al.，2021a）。例如，研发了微量元素和同位素分析测试技术，可以系统揭示和量化地幔交代作用（Saunders et al.，2015; Liptai et al.，2017; Kilgore et al.，2020）; 开发出 GLITTER 软件系统，处理同位素定年与微量元素数据（Griffin et al.，2008）以及Terrane Chron®地球化学填图新技术。特别是，O'Reilly和Griffin团队提出了岩石圈3D、4D填图概念（O'Reilly et al.，1996，2001），也是一种以捕虏体为基础的方法（Xenolith-based methodology），其将地球化学、地球物理学、构造地质学和年代学等多种学科综合在一起，分析岩石圈不同深部的性质（O'Reilly et al.，1996）; 或者是利用地幔物质资料、地球物理数据和地球动力学模型，确定大陆岩石圈地幔的成分、架构（Architecture）及演变的一种综合方法。前人将这些思路和方法主要应用于非洲、澳洲和全球的太古宙克拉通岩石圈地幔的研究（Begg et al.，2009; Belousova et al.，2009; Griffin et al.，2009，2013）。通过地球物理手段（如地震波各向异性），追踪上地幔不均一性和岩石圈不连续性，揭示克拉通岩石圈热结构和演化（Özaydin et al.，2021，2022; Afonso et al.，2022）; 在基性岩/包体及同位素填图（岩石探针）基础上，结合地震等地球物理手段，提出并实践了岩石圈3D-4D填图（Belousova et al.，2009; Griffin et al.，2013; Mole et al.，2014; Xu Bo et al.，2021a）。此外，澳大利亚地球科学局、西澳大利亚地质调查局等机构开始了全澳尺度的中酸性岩石Nd/Hf同位素填图（Champion et al.，2016; Lu Yongjun et al.，2021a，2021b），用于示踪克拉通边界、克拉通地壳架构及指导找矿等。
-
但是，这些研究仍处于探索阶段，还存在以下5个问题:① 这些研究主要依赖基性岩及地幔包体/捕虏体，即“以捕虏体为基础的方法”，但这类包体/捕虏体出露有限，且地表最发育的中酸性岩浆岩（如花岗岩等）未充分利用，因此，该方法难以获得更多的深部物质信息; ② 上述3D填图多用于探讨地幔物质架构，对地壳深部物质架构的填图和研究极少，未实现整个地球圈层特别是地壳圈层的物质探测; ③ 这些研究集于大区域的物质（地幔）架构示踪，小范围、高精细，特别是造山带和矿集区范围的示踪研究少; ④ 缺少多元同位素联合示踪及新兴同位素方法的应用，在大区域尺度上，地球化学研究与地球物理、实验岩石学研究结合不紧密; ⑤ 缺少不同方法之间的对比研究，多元同位素之间的关联性，包括传统放射性同位素与新兴稳定同位素关联性不清楚，尚未形成系统的深部物质架构探测的技术方法体系。
-
上述问题的存在制约了研究者更全面、更精细地获取地球深部物质的信息。因而，岩石圈深部探测，特别是深部物质探测技术亟需发展和完善，以便总结出一套能与地球物理探测相匹配的深部物质架构探测技术，提升对地球深部物质分布与演化的认识。因此，需要充分利用地表广泛出露的来自深部不同圈层的各类岩浆岩，排除其他干扰因素，示踪深部源区物质特征。特别是，不局限于某些点上的仅限于岩浆岩成因与物源特征研究，而是更关注深部物质在横向上（水平）、纵向上（垂向）分布的宏观特征，即三维物质架构，甚至四维演变。同时，需要探索研究适用于地壳，特别是造山带地壳的区域性深部物质架构的示踪和探测技术方法，加强岩浆岩Nd、Hf等同位素方法示踪结果的异同及其关联性研究，加强区域同位素填图示踪结果与地球物理探测、实验模拟结果的对应关系研究，从而建立一套揭示三维岩石圈及地球深部物质架构探究方法体系。这样，有望实现地球物理探测与岩石探针、同位素填图的深部物质探测紧密结合，以及岩石圈深部结构探测与深部物质探测的结合，改变长期以来地球深部探测研究缺少深部物质组成架构研究的现状，推动全方位、全要素的岩石圈深部探测工作。
-
本文在前人研究的基础上，重点探索通过“岩石探针”及岩浆岩多元同位素区域填图研究，结合地球物理探测、实验模拟研究，探测深部物质架构及其演变，并探索总结揭示三维岩石圈物质组成架构的方法体系框架，为开展深部物质探测及深部物质演化等研究提供新的思路和方法。
-
1 深部物质架构示踪的方法体系框架
-
近年来，通过同位素填图探测深部物质架构的研究日益增长，应用领域亦不断扩展，如探测深部岩石圈组成与演化、划分不同地球化学省和块体、刻画陆壳生长与造山带发育模式、分析区域成矿规律（王涛等，2018; 侯增谦等，2018; 莫宣学，2019）。依据已有的研究，综合中亚造山带、青藏高原、中央造山带和华北克-扬子克拉通等典型大地构造单元以及美国、俄罗斯、澳大利亚较成熟的深部物质示踪和探测研究实例，结合地球物理和实验模拟技术，笔者初步提出以岩石探针+同位素示踪为核心的、结合地球物理和实验模拟的、揭示三维岩石圈物质架构的方法体系框架（图1，表1）。表1所列前四类方法属于物质示踪探测技术，依据研究的对象和采用的方法，主要包括岩石探针、多元同位素填图、捕获/继承锆石填图、新兴同位素示踪与填图四种方法。
-
图1 以岩石探针和同位素填图为核心探索岩石圈三维物质架构的方法体系原理示意图
-
Fig.1 Schematic diagram of the methodological system for deciphering3-D material architecture of the lithosphere, with rock probe and multi-proxy isotopic mapping as the core
-
1.1 岩石探针
-
岩浆岩及其携带的深源岩石包体或捕虏体可以直接带来深部岩石圈组成结构的信息，包括物理化学条件与参数，以及岩石圈随时间演化过程。特别是穿越不同圈层的物质样品（岩石、矿物等），是窥探地球深部的“超深钻”和“望远镜”，提供了透视深部物质及其演变和揭示地球深部地质作用和成矿效应等深部过程的直观证据，也是地球深部演化过程的重要物质记录。因此，它们被称作探测地球深部的“探针”（Lithoprobe或Rock probe）和“窗口”，可以应用于探测深部物质组成、架构，重建板块构造过程中的大洋与大陆演化历史，揭示壳幔相互作用、地球深部高级变质作用和成矿过程等方面（莫宣学，2011）。
-
由于岩浆岩分布广泛，可以获得更多的深部信息。因此，除充分利用岩石捕虏体外，还可进一步针对广泛分布的岩浆岩开展多种测试分析，获得岩浆岩岩相学和地球化学数据，探测深部信息，特别是开展区域深部物质分布的示踪、探测，以便与地球物理探测紧密结合。
-
笔者认为，岩石探针方法可以泛指利用岩浆岩及其携带的深源岩石包体或捕虏体，开展多种手段的测试研究获取深部信息的方法; 可以理解为广义的岩石探针。而对其集中利用岩浆岩包体或捕虏体，开展研究探索深部物质类型和物化参数特征的方法，可以理解为狭义的岩石探针或岩相学探针，即 “包体或捕虏体方法”（Xenolith-based methodology）。岩石探针方法的研究应用主要包括以下两个方面。
-
1.1.1 深源岩石包体或捕虏体岩相学研究
-
岩浆岩深源岩石包体或捕虏体是深部岩石样品的直接代表，主要来源于地球深部的火山或岩浆所携带的深部岩石圈碎块，是深部物质等探测的“金钉子”。地幔岩捕虏体及捕虏晶直接带来地幔的物质信息，通过对其中矿物特征及成分、形成温度和压力条件等的深入研究，可以确定岩石圈的厚度、组成、物理化学性质及其长期的演化特征，是示踪地幔物质与物质状态的常用方法，故称之为“以捕虏体为基础的方法”（Xenolith-based methodology; O'Reilly et al.，1996）。
-
依据包体或捕虏体来源地幔、地壳等不同深度和物源，可以分别探测地壳和地幔的深部物质组成、物理化学属性等。例如，在一些古老克拉通，利用众多的地幔岩捕虏体探索大陆岩石圈地幔的特性。这方面国内外有大量研究实例（Griffin et al.，2009）。在华北克拉通，利用这类捕虏体结合地球物理反演出深部的结构和状态，取得了很好的效果（郑建平等，2021）。在一些造山带，如青藏高原（如Zhang Zeming et al.，2014），依据来自下地壳的麻粒岩、斜长角闪质捕虏体，直接确定了下地壳的组成、状态及深度。
-
此外，除上述岩浆期结晶形成的主要矿物和副矿物外，岩浆在形成和上升中还携带有继承性或捕虏来源的矿物，如捕获锆石、继承锆石（核）等。继承锆石直接带来了岩浆源区的信息。捕获锆石较复杂，可以来自岩浆源区，也可以来自岩浆岩围岩。对于深成岩而言，基本可以了解岩体定位层次之下的物质，但对于火山岩而言，不排除来自上地壳及地表围岩的可能。近年来，随着测试手段的不断进步，有学者还开展了区分岩浆锆石、独居石、磷灰石与继承锆石、继承独居石和继承磷灰石等相关信息的深入研究（曾令森等，2012; Zhang Jianjun et al.，2015，2017; 刘志超等，2020; Ding Huixia et al.，2021）。
-
1.1.2 岩浆岩岩相学（岩石学和矿物学）研究
-
除了深源岩石捕虏体外，对岩浆岩及其组合开展系统岩相学和岩石成因研究，排除其他因素干扰，也可以推断源区特征，如源区基本的岩石类型或组成等; 再结合年代学、地球化学（包括同位素等）和区域地质背景（如流变学、变质作用等）研究，可以示踪岩浆岩的源岩性质与演变，从而示踪地球深部物质特征。一般而言，通过中酸性岩浆岩研究，可以推测下地壳的物质特征，包括可能的岩石类型; 通过镁铁质-超镁铁质岩石类型、组合研究，结合深源捕虏体和高压矿物及地球化学示踪，可以示踪地幔（岩石圈或软流圈地幔）组成。
-
随着微区原位测试分析技术的发展与应用，矿物学及其原位地球化学研究可以提供更多的深部信息。很多岩浆岩的造岩矿物和副矿物都含有深部的信息。
-
可见，岩石探针提供了直接窥探深部信息的途径，具有明显的优势。莫宣学（2011）提出，岩石探针至少可以提供以下深部信息:① 壳幔物质组成与结构; ② 壳幔的热结构和热状态; ③ 地壳及岩石圈厚度及其空间变化; ④ 软流圈顶面埋深、温压、物质状态、流体或熔浆含量; ⑤ 壳幔氧化—还原状态; ⑥ 深部流体特征; ⑦ 可以提供壳、幔上述各种性质和参数随时间的变化，从而反演壳幔深部过程，这个优点是其他深部探测方法所不具备的。
-
应用岩石探针方法示踪深部物质和深部过程，前人已取得很好的成果。在中国东部特别是华北克拉通，研究者通过对岩浆岩和深源包体的研究推断出了岩石圈不同深度的深部物质组成及结构信息（池际尚等，1988; 邓晋福等，2003; 高山等，2009）。近期，郑建平等（2021）依据宽角反射/折射地震剖面约束的NCC地壳结构模型，结合深部地壳捕虏体的系统研究，多维度地限定了华北地块的地壳物理、化学结构特征，揭示了该区地壳精细的物理化学结构，进而探索其形成演化过程。
-
但是，应用岩石探针探测深部物质也有一定的局限性。因为，深源岩石包体或捕虏体在地表出露比较少，其应用受到很大限制。岩浆岩分布广泛，但是，仅依据岩相学研究直接推断深部源区物质分布状态有难度。开展分析测试研究可以获得更多信息，但仅一个点的研究，也具有时间和空间上的局限性。因此，岩石探针的研究应用还需要与其他物质探测技术，如同位素填图、岩浆岩捕获锆石填图及地球物理方法相结合。
-
1.2 多元同位素填图
-
岩浆岩同位素研究可以示踪源区的某些特征，例如深部物质富集与亏损、地壳物质的新老信息及演化过程。特别是，这方面的研究可以充分利用岩浆岩分布广泛的优势，开展区域性的同位素示踪填图工作，可以了解中下地壳、岩石圈及地幔各类物质目前的时空分布特征。国内外学者开展了这方面的探索（王涛等，2018）。近年来，各类同位素数据（如锆石Hf-O同位素）大量涌现，加之GIS成图技术的进步与推广，国内外出现了很多的依据大量同位素数据，特别是其时空变化特征了解深部物质的分布及物质架构，解决一些重大地质问题。同位素填图形成了有效的、相对独立的研究方法。
-
1.2.1 主要方法
-
1.2.1.1 放射性同位素示踪填图
-
放射性同位素一般以初始（Primitive）到演化（Evolved）的同位素特征来示踪初始或新生（Juvenile）到改造或再造（Reworked）的地壳特征，或亏损到富集的地幔特征，具有一定的时间和物质演变属性。我们将放射性同位素示踪填图分为两类，一类是全岩放射性同位素的数据分析与填图; 另一类是单矿物的方向性同位素的数据分析与填图。
-
全岩Nd同位素填图:全岩Sm-Nd同位素体系是较早应用于岩石地球化学领域相对成熟的方法，可客观反映样品的总体特征，更有利于开展区域同位素填图和对比。利用Nd同位素开展区域变化分析，可以解决很多区域地质地质问题（Hong Dawei et al.，2004; Zhou Xinmin et al.，2006; Wang Tao et al.，2009; Yang Qidi et al.，2017; 王涛等，2018）。但应该注意是，其同位素值是多元物质混合的结果，例如，除记录来自源区的信息外，全岩Nd同位素还记录了上升及定位过程中不同岩石圈层次的围岩信息。故往往全岩Nd同位素反映的信息较锆石Hf同位素更富集，且模式年龄偏老。此外，Nd模式年龄的计算还受到计算过程中参数取值等因素影响，有时误差较大。
-
全岩Sr、Pb等多元同位素填图:随着测试技术的提高，很多其他放射性同位素体系得以广泛应用，且矿物原位微区同位素分析得以进一步普及，为开展区域多元同位素填图奠定了基础。例如，全岩Sr同位素应用较早，常常与Nd同位素一起示踪岩浆岩源区的特征。但全岩Sr同位素体系易受到岩浆分异和流体作用的影响，常常出现异常值。近年来，斜长石、磷灰石等矿物原位Sr同位素分析技术的开发，为解决这一问题提供了新的思路。如Xu Bo et al.（2022）研制了磷灰石的Sr同位素标准物，开展了青藏高原的Sr同位素填图研究（详见后述）。
-
Pb同位素较早应用于区域同位素示踪分析。朱炳泉等（1995，1998）还利用Pb同位素等地球化学区域分析，提出了地球化学激变带，也是区域同位素填图的典范（朱炳泉，2001）。但受制于严格的实验分析条件（如对实验环境的严苛要求），Pb同位素长期未能广泛应用于物质示踪研究。近年来，随着国内实验室条件的普遍提升，关键区带岩浆岩Pb同位素数据量不断增加，有望开展岩浆岩Pb同位素填图实践。由于大陆地壳中Pb同位素丰度远高于亏损地幔和洋中脊玄武岩，如果古老地壳中加入亏损地幔属性的物质，其Pb同位素不会明显变化; 但是反之，会使得地壳中Pb同位素发生显著变化（Sun et al.，1989; 郑永飞等，1998）。根据这一特性，Pb同位素体系在示踪年轻地壳为主的区域内零星分布的古老物质方面有望取得较好的效果。例如，在Sr-Nd同位素填图的基础上，徐盛林等（2022）在西准噶尔地区开展了岩浆岩Pb同位素填图，将西准地区分为高值[（206Pb/204Pb）i＞18.30]、中值[（206Pb/204Pb）i=18.00~18.30]、低值[（206Pb/204Pb）i＜18.00]三个Pb同位素省。
-
锆石原位Hf同位素填图:锆石Hf同位素更稳定，能反映锆石不同部位的原位特征，有利于示踪卷入岩浆事件的年轻和古老物源。开展区域锆石Hf同位素填图需注意:① 一个样品的Hf同位素值变化极大，需要取代表性数值（如算数平均值、中位数、加权平均值）; ② 测试锆石或锆石部位及最终数据结果受人为因素影响，常丢失古老的物质信息，这是造成锆石Hf与全岩Nd同位素参数间解耦的原因之一。因此，在同位素填图中，最好将Nd、Hf两者结合起来分析。
-
近年来，岩浆岩锆石Hf同位素数据量呈爆发式增长趋势，Hf同位素填图的实例也日趋增多。在开展区域岩浆岩Hf同位素填图时，一般以样品基本单位，以每个样品中自结晶锆石的平均值/中位数值代表该样品的Hf同位素特征。
-
应该注意的是，在很多样品中，岩浆期自结晶锆石（不包括捕获锆石、后期热事件中同位素重置的锆石等）的Hf同位素可以呈现一个较宽泛的分布范围（可达十余个epsilon单位），且会有少数锆石测点结果显著地偏离主体测试结果。此时，必然涉及两个问题:是否需要剔除异常值?以什么方式判别异常值?这取决于填图或示踪研究的目的。若是区域Hf同位素填图，特别是需要了解岩浆事件主要的物源特征，需要寻找样品代表性的特征值，可以考虑岩石学、地球化学的基本原理，选用箱线图（Box plot）的方法进行异常值的筛查和剔除。而涉及锆石Hf同位素相关参数的平均值或中位数值时，仅“正常数值范围内”自结晶锆石的相关参数参与计算。若需要示踪了解深部古老物质，则不仅需要特别关注老的继承锆石或捕获锆石，自结晶锆石中“异常点”的影响也不应忽略，以期和Nd等全岩同位素（含有古老物质信息）对比。需要说明的是，此处所指的异常值是统计学意义上的，被标记为异常值的锆石Hf同位素数据，并非没有岩石学和地球化学研究意义。在具体的岩体和火山岩研究中，少数离群值的出现，可能对应某些特殊的深部物质和深部过程，是值得研究的问题。但是在涉及大区域、海量样品的锆石Hf同位素数据时，需要先去除这些统计学意义上的异常值，否则会对样品层面总结的平均值、中位数、均方差等参数施加过大的影响。
-
除岩浆自结晶锆石外，岩浆岩中另一类锆石的Hf同位素信息在实际研究中常被忽略，即岩浆岩中古老的继承/捕获锆石。若需要示踪了解深部古老物质或探索锆石Hf同位素和全岩Nd等同位素体系（含有古老物质信息）的对应关系时，捕获锆石中的同位素信息应予以充分考虑。在区域尺度同位素填图中，还可进一步开展捕获锆石填图（详见后述）。
-
其他单矿物同位素。随着微区原位测试分析技术的发展与应用，矿物学及其原位地球化学研究可以提供更多的深部信息。很多岩浆岩的造岩矿物和副矿物都含有深部的信息。例如，磷灰石作为岩浆岩中常见的副矿物，广泛分布在各类火成岩中，是探究不同体系同位素对标的理想矿物。而且，磷灰石相较于锆石含有更多的微量元素，并且晶格中存在与矿床相关的挥发分元素，如Cl、S; 在演化程度较低的岩浆系统中，对岩浆过程非常敏感。因此，磷灰石中的主量和微量元素及Sr-O同位素数据可以对锆石Hf-O同位素信息进行补充。通过建立不同磷灰石Sr-O同位素的标准物质（Xu Bo et al.，2022），厘定磷灰石Sr-O同位素和挥发分（如Cl、S）含量的地球化学指标，与锆石Hf-O同位素相结合，可以更好地认识岩石成因和成矿过程（Xu Bo et al.，2022）。例如，Xu Bo et al.（2020）通过对多个天然磷灰石样品采用红外-拉曼光谱、X射线衍射、激光烧蚀电感耦合等离子体质谱（LA-ICP-MS）、电子探针（EMPA）、显微X射线荧光光谱（micro-XRF）分析技术及LA-MC-ICP-MS和同步辐射技术（micro-XANES），厘定出不同产地的磷灰石矿物学特征，为探讨天然磷灰石的化学成分、同位素特征和氧化还原状态奠定了基础。通过对特提斯成矿带的含矿斑岩和贫矿岩浆系统内磷灰石的矿物学地球化学研究，揭示碰撞环境下成矿斑岩中磷灰石具有更高的Cl和S含量，指示斑岩体具有较高的成矿潜力。由于Cl和S是成矿流体的重要成分，因此可以揭示挥发分在斑岩系统中迁移和富集的重要过程。结合碰撞斑岩型矿床的87Sr/86Sr同位素特征和模拟计算，对碰撞环境中斑岩系统的挥发分建立了双阶段演化模型，明确了这些挥发分的来源、演化过程及斑岩系统的成矿特征（Xu Bo et al.，2021b）。
-
1.2.1.2 稳定同位素示踪和填图
-
相对于放射性同位素，O-Li等稳定同位素在示踪地表到地幔的熔/流体与矿物直接的相互作用过程中有着独特的优势。O同位素可以用来鉴别高温或低温条件下水-岩相互作用中再循环陆表组分的贡献。尽管洋壳的不同部位具有不同的O同位素组成，但是具有相似的放射性同位素组成。因此，利用O同位素与放射性同位素（如Hf同位素）的相关性，可以判断俯冲的板片和岩石圈之间的流体作用。
-
原位锆石O同位素是近年来一个新兴的分析手段。相对于其他稳定同位素体系，锆石O同位素具有分析快速的特点，有望广泛应用于区域性同位素填图中。原位锆石O同位素组分主要受地表过程的水-岩相互作用影响。而锆石是一种高温、难熔、化学性质稳定的副矿物，受后期岩浆结晶和变质作用的影响较小，可以很好地保存原始岩浆中锆石的O同位素组分（Valley et al.，1994，2005; King et al.，1998）。正常幔源岩浆结晶出来的锆石具有非常一致的δ18O值（5.3‰±0.6‰），该比值受岩浆分异的影响很小（Valley et al.，1998），而经历表壳过程的物质，其O同位素值显著增加。因此，O同位素可以用来鉴别高温或低温条件下的水-岩相互作用的含水组分（陆表组分或幔源），对于分析陆表物质（低温/高温下的水岩反应）的贡献具有重要的意义。尤其在一些典型的新生洋内弧背景，所有的物质都具有相似的亏损的Sr-Nd-Hf同位素特征，接近亏损地幔值（Chen Bin et al.，2005; Geng Hongyan et al.，2009; Tang Gongjian et al.，2012）。这种情况下，传统的放射性Sr-Nd-Hf同位素很难识别岩浆源区组分，尤其是循环的表壳岩组分。最近，有研究者对西准噶尔、东准噶尔等地区的花岗岩开展了锆石O同位素研究（Tang Gongjian et al.，2019; Yin Jiyuan et al.，2020; Zhang Yunying et al.，2022）。很多碱性花岗岩含有非常高的锆石O同位素值（δ18O＞8‰），表明源区含有大量的表壳岩。通过Hf-O同位素模拟估算，其表壳岩的组分可达50%（Tang Gongjian et al.，2019; Yin Jiyuan et al.，2020）。在Nd-Hf填图的基础上，尹继元等（待发表）开展了O同位素区域分布研究及O同位素填图，进一步鉴别出经历了地表物质循环的深部物质的时空分布特征。
-
此外，全岩Li-B稳定同位素可以进一步甄别陆表到地幔之间的流体作用（俯冲板片或海水或沉积物），如相对地幔或海水，轻的Li-B同位素指示了岩浆源区陆壳沉积物的加入。因此，可以通过全岩Li-B等稳定同位素剖面研究，进一步揭示俯冲板片相关的沉积物、流体或熔体作用。如O-Li-B-Mg同位素研究表明，沉积物质在不同的地壳深度下对松潘-甘孜造山带准同期（213~198 Ma）不同花岗岩类型均具有重要贡献（Li Shan et al.，2021）。其中，S型特征花岗岩主要源自中地壳中变沉积物的部分熔融，沉积物贡献最多; 其他花岗岩类型的源区均有沉积物的贡献，且主要源自更深层地壳中变沉积物，但相对较少，因而显示I型和A型特征。
-
随着测试技术的提高，一些新兴同位素测试数据也将急剧增长，也可能发展为新的同位素示踪填图技术。
-
1.2.2 多元同位素解耦、示踪的关联性及联合填图
-
1.2.2.1 同位素解耦及示踪结果的关联性
-
同类的同位素具有相同和相似的示踪原理。例如，Sm-Nd和Lu-Hf等都属于放射性同位素。它们都以初始（Primitive）到演化（Evolved）的同位素特征来示踪初始或新生（Juvenile）到改造或再造（Reworked）的地壳特征，或亏损到富集的地幔特征。因此，有可能相互关联，联合示踪。实现多元同位素联合示踪填图，实现需要研究和了解不同放射性同位素体系方法的相关性与解耦性。在实际研究中，不同放射性同位素体系都具有一定的相关性，但也常出现一定程度的解耦现象。Nd-Hf同位素解耦引起的源区示踪不一致，导致不同学者对花岗岩类代表的构造意义的解释存在明显差异，例如阿尔泰和天山等（Yu Yang et al.，2017; Kong Xiangye et al.，2018; Luo Qun et al.，2019; Huang He et al.，2020）。关于可能导致Nd-Hf同位素发生解耦的成因已有不少研究，可能的因素包括:① 矿物效应:锆石（如Vervoort et al.，2000）和石榴子石效应（Patchett et al.，2004; Schmitz et al.，2004），即部分矿物由于赋存母子体元素的差异，在地质过程的作用下会导致元素异常富集; ② 地幔交代作用:熔体/流体富轻-中稀土元素（包括 Sm、Nd），贫重稀土元素Lu和Hf，对被交代地幔的Sm-Nd同位素体系的影响强于Lu-Hf同位素体系导致Nd-Hf同位素组成解耦（Choi et al.，2012）和软流圈硅酸盐熔体的层离渗透作用（Bizimis et al.，2004; Nowell et al.，2004）等; ③ 地幔中的大洋沉积物组分再循环作用（Garçon et al.，2013）; ④ 地表风化作用（Vervoort et al.，2011）。此外，还可能是因为锆石矿物Hf同位素信息提取不全（Guitreau et al.，2014; Wan Yusheng et al.，2015），即样品岩浆结晶锆石反映的Hf同位素组成，有时可能不（完全）代表其初始岩浆Hf或寄主岩石的Hf同位素特征（吴福元等，2007）。由于Nd同位素主要来自于样品全岩组成，而Hf同位素主要来自于单颗粒锆石，且锆石稳定性高，即使在后期高温事件中依然可保持良好的封闭性，在岩浆作用过程中再循环地壳物质的参与将产生大量继承锆石。
-
在实际研究中，需要了解和评估这些可能的解耦成因。笔者对北疆地区的研究显示，阿尔泰志留纪—泥盆纪岩浆岩普遍出现Nd-Hf同位素解耦的现象，但其成因仍存在争议，目前的解释有:可能继承自其源区特征，也可能由于部分熔融过程中特殊的矿物效应导致（Huang He et al.，2020; Zhang et al.，2020）。而在西准噶尔—西天山地区，Nd-Hf同位素在联合示踪深部物质方面具有较好的一致性和有效性，对应性关系为:长英质岩石εHf（t）锆石=1.0772εNd（t）全岩+5.6916，以及铁镁质岩石的对应关系为:εHf（t）锆石=1. 0739εNd（t）全岩+7.3715（图2）。
-
1.2.2.2 多元同位素体系关系（Nd、Hf等）及联合示踪填图
-
了解一个地区特定的岩石类型不同同位素体系的关联性，可以将多元同位素示踪结果相互校正，互为约束，找到可能联合示踪的依据。例如，可以明确以一种同位素为主，如以εNd（t）=0为基准，判别为年轻和古老物质，其放射性同位素的参数（如锆石Hf同位素）值相应地调整（如εNd（t）=+6）。此外，当数据量有限时，可以用Nd、Hf联合示踪填图，以获得相对客观的物质架构示踪填图结果。
-
在同位素填图中，可以以公认的、更广泛应用的Nd同位素为基础，刻画同位素省的特征。例如，依据中亚造山带的研究，以Nd同位素为特征，可将同位素省分为六类:Ⅰ:极度初始（Highly or extremely primitive，εNd（t）＞+6）、Ⅱ:初始（Primitive; εNd（t）=+6~+4）、Ⅲ:轻微初始（Slightly primitive，+4~0）; Ⅳ:轻微演变（Slightly evolved，εNd（t）=0~-4）、Ⅴ: 演变的（Evolved，εNd（t）=-4~-10）、Ⅵ:强烈演变（Highly evolved，εNd（t）＜-10）。对地壳而言，这些同位素省分别示踪和刻画了地壳的类型和特征:非常年轻（Highly juvenile）、年轻（Juvenile）、略微年轻（Slightly juvenile）、略微再造（Slightly reworked）、再造（Reworked）、强烈再造（Highly reworked）。
-
以此为基础，依据其他同位素如Hf同位素与Nd同位素的关系，可以换算参数，共同刻画这些同位素省的特征，如相应的六类同位素省的Hf同位素特征是:Ⅰ:极度初始（Highly or extremely primitive，εHf（t）＞+12）、Ⅱ:初始（Primitive，εHf（t）=+12~+10）、Ⅲ:轻微初始（Slightly primitive，εHf（t）+10~+6）:Ⅳ:轻微演变（Slightly evolved，εHf（t）=+6~+2）、 Ⅴ:演变的（Evolved，εHf（t）=+2~-4）、Ⅵ类:强烈演变（Highly evolved，εHf（t）＜-4）。例如，王涛等（2020）对北疆基性岩开展同位素填图，示踪地幔的物质特征及其分布，应用Nd、Hf两种数据（将εNd（t）减小6，位置等同于εNd（t））联合示踪填图，获得了较好的效果（图3）。
-
1.2.3 同位素填图有关问题及讨论
-
1.2.3.1 需要注意的问题
-
开展同位素填图时，一是要考虑样品的时间（结晶年龄）; 二是要考虑到样品岩石形成的深度范围。应用εNd、εHf参数开始同位素填图的时候，需要确定一定的时限，如针对每个时期的岩浆岩计算这些参数，并开展区域对比分析。同时用于同位素填图的样品年龄区间越小越好，不能将样品形成年龄差异很大的参数结果简单放在一起开展同位素填图。不过，一般而言，较小时间区间（如100~50 Ma）内样品年龄的差异对这些参数的影响很有限，如中生代（200~145 Ma），可以开展整个中生代岩浆岩的εNd、εHf填图。如有必要，其他时代的岩浆岩可以将样品年龄换算到这个时间，来参与填图，以便利用更多的样品。
-
图2 西天山-西准噶尔典型地区长英质岩石（a）和基性岩石（b）全岩Nd同位素与锆石Hf同位素参数的对应关系图解（数据来自Huang He et al.，2020及项目组未发表资料）
-
Fig.2 Diagram of whole-rock εNd (t) value vs. average εHf (t) value of individual felsic (a) and mafic rocks (b) in the West Junggar and West Tianshan, demonstrating a linear relationship between the two radiogenic isotope systems (based on data from Huang He et al., 2020 and unpublished data)
-
图3 北疆阿尔泰-西准噶尔-天山全岩Nd、锆石Hf和捕获/继承锆石廊带填图及剖面对比（西准噶尔为年轻地壳，阿尔泰、天山为较古老地壳，与地球物理在30~40 km深度的面波Vs探测结果和参数可对比）
-
Fig.3 Cross sections of whole-rock Nd isotopic, zircon Hf isotopic, and ages of inherited/xenocryst zircon grains of magmatic rocks from the Altai-West Junggar-West Tianshan regions (this results together indicate the dominance of juvenile crust in West Junggar and dominance of ancient crust in Altai and Tianshan regions)
-
资料来源:（a）—42.5 km深度面波Vs速度插值结果（据鲍学伟等未发表数据）;（b）—北疆地区Nd同位素填图结果（据王涛等，2022）;（c）—42.5 km深度面波Vs速度插值结果生成的剖面;（d）—生成剖面图的岩浆岩全岩Nd同位素;（e）—捕获锆石Hf同位素（据Huang He et al.，2020及项目组未发表资料）
-
Data sources: (a) —the interpolation result of Vs at the depth of 42.5 km (based on unpublished data by Bao Xuewei et al.) ; (b) —Nd isotopic mapping result of granitoids from northern Xinjiang (after Wang Tao et al., 2022) ; (c) —cross section of the interpolation result of Vs; (d) —cross sections of whole-rock Nd isotopic; (e) —xenolith zircon Hf isotopic results of granitoids in northern Xinjiang (based on data from Huang He et al., 2020 and unpublished data)
-
应用Nd、Hf模式年龄填图的时候，则可以不考虑样品的时限，所有同类岩性的样品均可以用于填图。因为，同一个源区，不管什么时间发生熔融形成岩浆，模式年龄都是一样的。应注意的是，模式年龄误差很大，需要依据与相对应的εNd、εHf参数做相关性分析，将相关性不好的模式年龄排除，挑选出数据质量高、误差小的模式年龄开展填图。
-
此外，开展同位素填图，特别是开展区域分析的时候，要考虑样品的岩性及其可能来源的深度。不能将岩性差别极大的样品，放在一张同位素图上分析。例如，花岗质样品反映的是中下地壳的组成; 而基性岩一般反映的是地幔的组成。这两类岩石的同位素参数不能在一张图上分析区域间的差别。
-
总之，开展同位素填图，对比不同区域的深部物质时候，最好用同一期、同一类岩浆岩，这样才具有可比性，以便开展区域分析和对比研究。若样品有限，且判定早期深部物质架构有没有被后期事件大量改造破坏时，可以将不同期次的岩浆岩按统一的一个时期的年龄来计算，回归到同一个时期，获得该年龄范围内的完善参数，可以用于同位素填图和分析对比。
-
1.2.3.2 常见同位素图表达方法及其优缺点
-
目前，区域 Nd 同位素填图表达方式可分为两类:等值线图和岩体同位素属性填图。同位素等值线图直接依据同位素数值点，应用软件做出等值线图。该方法简便，但存在一定不足。岩体同位素属性图是用同位素参数将岩体属性化，用不同的颜色表达岩体的同位素值，如εNd（t）值图或模式年龄值（Wang Tao et al.，2009; 王涛等，2020）。这种岩体的同位素图类似传统地质图中的岩浆岩图（差异的色度表征不同的岩体时代），但颜色映射的是岩体不同的同位素属性。在此基础上，可以划分不同的同位素省，并据此划分出地壳省，如在阿尔泰、大兴安岭及邻区等地区开展的花岗岩Nd同位素填图（Wang Tao et al.，2009; Yang Qidi et al.，2017）。该方法更客观地反映了实际情况，但需要海量的基础地质信息支撑（如中-大比例尺地质图中各个岩体的大小和分布），且成分变化复杂的岩浆岩Nd同位素值需要取平均值。对于差别过大的同一个岩体，可能需要解体为不同岩体，赋予不同的同位素属性。
-
1.2.3.3 存在问题
-
在实际应用中，需要了解不同同位素填图的方法特点、适用性及在确定三维岩石圈架构中的存在问题。① 等值线同位素图简单，但可能将问题简单化，其研究有模式化的趋向。这种方法很大程度上取决于样品数据的多少和分布均匀，数据点不够或不均匀，可能夸大某些样品数值所代表的空间范围，给出的结果不够客观、全面，而且往往忽视了地质背景，如在生成同位素模式年龄等值线时常缺少构造线或构造单元边界约束。因此，仅依赖“数据”容易得出片面的结论。此外，等值线图的统计方法较多（如克里金法、经验贝叶斯法、反距离权重法），每种方法也可设定不同的条件参数，在样品非等间距、非等密度分布时，各种方法产生的等值线图结果略有差异。② 不同同位素示踪方法（如全岩Sr-Nd和锆石Lu-Hf同位素）给出的判别参数不一定一致，如常将εNd（t）= 0或εHf（t）=0作为判别新生和古老物质的参数，但是，这两个参数是不一致的，在估算的年轻地壳比例差别很大。这就需要了解不同方法的关联性，进行校对。如上文所述，在北疆地区长英质岩石εHf（t）=+6，相当于εNd（t）=0。③ 在实际应用中，不能将不同同位素研究对象和不同测试方法的结果，简单地进行等值线图对比，会导致解释的不真实性，如一个地区的基性岩同位素数据不能和另一个地区的花岗岩数据一起进行填图分析。④ 在进行同位素参数分析填图中，不能忽略岩浆岩形成的时间。如前所述，有时由于不同时代岩浆的时间差别不大，不足以影响同位素参数的计算; 但有些情况下，时代间隔太长，将同位素参数放在一张图上分析，缺乏合理性。可以考虑研究的具体问题，用统一的时间换算同位素比值。
-
此外，经常遇到的一个问题是岩浆成因复杂性对同位素示踪和填图的影响，如地壳混染、岩浆混合等。同位素示踪和填图结果本身是客观的，关键的问题是如何认识和解释这些结果。例如，中酸性岩石中年轻的（高正εHf和εNd）同位素特征既可能来源于年轻地壳，也可能反映了岩浆混合的结果，单纯地通过同位素值无法区分这两个过程。对于这一问题，一方面，无论是否存在岩浆混合，同位素填图的结果都反映了新生物质的存在，区别在于新生物质是在岩浆事件前就加入到深部源区（源区混合），或是伴随岩浆过程幔源岩浆进入到岩石圈中。另一方面，如需精确解剖岩浆过程，还可以开展其他研究加以识别，如基于岩相学研究判断是否有岩浆混合; 利用矿物原位同位素（Sr、O）研究，可以揭示是否有地壳混染等。将岩浆岩同位素填图与矿物原位同位素填图相结合，可以更客地分析应用同位素填图结果。
-
每种方法的优势和不足，如何联合示踪填图，相互作用和约束，总结Nd、Hf等同位素填图方法的优缺点及其相关关系，是需要深入研究的重要问题。
-
1.3 岩浆岩捕获/继承锆石信息填图
-
1.3.1 基本原理与方法
-
区域岩浆岩捕获/继承锆石信息集成（或者说微区信息填图）是探索深部古老物质组成的新思路和有效手段（Zhang Jianjun et al.，2015，2017）。锆石是一种广泛存在于各类岩浆岩中的副矿物，其物理化学性质稳定，是多种同位素体系理想的测试对象，为研究岩浆事件期次和源区特征等提供了重要依据。锆石U-Pb定年是最成熟和最成功的同位素年代学体系（Spencer et al.，2016），锆石O氧同位素可以用来鉴别高温或低温条件下水-岩相互作用的含水组分（陆表组分或幔源）的贡献（Valley et al.，2003; Spencer et al.，2014），而锆石Hf同位素很好地指示了岩浆源区和壳幔分异过程（吴福元等，2007; Kemp et al.，2007; Jeon et al.，2018）。锆石的微量元素可以用来计算其母岩浆的成分特征，约束其氧逸度和形成温压条件（Hoskin et al.，2000; Belousova et al.，2002; 吴元保等，2004; 赵振华，2010; 赵志丹等，2018; 邹心宇等，2021; Lee et al.，2021）。
-
应该注意的是，岩浆岩中除直接在岩浆形成的自结晶锆石外，还存在捕获自岩浆源区、岩浆通道之前岩浆事件或者浅部围岩中的继承或残留锆石及捕虏晶（如pre-magmatic zircons; Miller et al.，2007）。这些捕获/继承锆石的U-Pb年龄明显老于岩浆结晶年龄，因而在岩浆岩年龄计算中不予考虑或常被忽视。但是，越来越多的研究显示，这些捕获/继承锆石能提供揭示“隐藏”在大陆地壳深部的岩浆事件及物质组成的重要信息（Smyth et al.，2007; Siebel et al.，2009; Stern et al.，2010; Buys et al.，2014; Liu Dong et al.，2014）。
-
由于研究对象和拟解决的科学问题不同，一个或多个样品中发现的明显老于岩浆岩结晶年龄的那些捕获或继承性锆石，多被认为难以准确解释，往往只有U-Pb年龄，缺少微量元素和同位素信息。实际上，通过建立数据库，开展区域性捕获、继承性锆石记录信息的综合集成可以揭示某些隐藏的规律。如笔者依据大量数据的统计分析，特别是给出每个样品的空间位置，开展区域填图。对比这些捕获继承锆石的年龄信息，在某一时间切片内，不同构造单元岩浆事件中的捕获继承锆石年龄（期次）、区域性最古老年龄继承锆石的空间分布和Hf同位素资料的统计分析和集成总结，结合区域其他学科资料，综合分析这些新发现的地质意义。在增生型中阿拉善、阿尔泰等地区，笔者应用该方法统计分析了区域岩浆岩中的捕获、继承锆石记录的信息，示踪了区域新老物质分布的特征，结果与Nd、Hf同位素区域填图结果一致，也与地球物理探测结构结果接近（Zhang Jianjun et al.，2015，2017）。因而，笔者提出，随着资料的积累，岩浆岩中的捕获/继承锆石信息填图方法，可以作为岩石探针及同位素示踪深部地壳古老物质组成的一个重要补充（Zhang Jianjun et al.，2017）。
-
该思路和研究方法已在很多研究中被广泛接受，对岩浆岩中识别出的捕获/继承锆石的深入研究很快在更多的地方展开工作。例如，在碰撞型造山带（如Liu Dong et al.，2014）、华北克拉通南部（Qi Nan et al.，2021; Wang Peng et al.，2022）、华南（Wang Kai et al.，2020; Zhang Shitao et al.，2021; Xiang Lu et al.，2022）都有新的相关继承锆石信息填图为重大地质问题提供辅助约束的研究实例。这说明，岩浆岩捕获继承锆石微区信息填图是示踪深部物质的有效方法。
-
1.3.2 功能及优势
-
捕获、继承锆石直接提供了来自深部物质的信息，如年龄信息、物质成分（如Hf、O等同位素）信息。该方法可以解决同位素示踪填图过程中遇到的问题。一些情况下，仅依据Nd-Hf同位素特征变化，很难断定深部物质的构造属性及其边界。例如，在克拉通边界示踪出较年轻的深部物质，其是古老基底改造年轻化的结果?还是增生的较年轻的块体?这些问题仅依赖Nd、Hf等同位素示踪难以解答。而岩浆岩捕获/继承锆石统计分析和填图可以为解决该问题提供直接证据。在一些关键地区，特别是Nd-Hf解耦地区，研究捕获锆石并开展填图更重要。
-
总之，开展不同构造单元重要区带大剖面及关键地区岩浆岩区域（包括廊带）同位素（如全岩Nd、锆石Hf）填图，包括捕虏体、捕获锆石填图，建立相关数据库，可以示踪不同构造单元、不同深度深部物质结构异同，划分大地构造单元，了解不同构造单元深部物质组成特征，结合地球物理探测，建立深部物质组成结构三维模型; 结合岩浆岩年代学研究，构建四维物质结构及演化模型。这是实现从地球物理结构探测到岩石地球化学物质探测的重要思路和途径。
-
1.3.3 存在问题
-
在数据整理阶段，应认真核对原始出处的数据，如一些文献中捕获/继承锆石只包含于U-Pb年龄测试结果表中，未在谐和图中展示。应详细整理和记录岩浆岩的岩石类型、样品经纬度、测试数据、测试单位、挑选锆石单位、测试方法、测试发表时间，锆石阴极发光/背散射图像结构特征（长宽比、核、幔、边，包裹体、裂隙、环带特征等）、数据处理软件等信息，原始文献对捕获或继承锆石原始数据的解释（锆石成因研究，如推测为捕获围岩或之前岩浆事件、或继承源区等）、有无相应的更多测试数据（如微量原始或Hf、O同位素数据）。这些信息在入库之后，在确定数据准确性、重现性和区域对比分析其地质意义时可发挥重要作用。
-
在数据分析阶段，有时捕获、继承锆石的成因和来源无法判断，在未开展进一步工作的情况下，有时难以直接区分。因而，在区域性集成分析时，要注意依据相关信息数据，先梳理出或大胆假设出两个模型端元，如假设锆石都是继承锆石，其成因反映了岩浆源区的信息（如源区古老物质成分信息）; 或者它们均为岩浆上升过程中捕获岩浆通道和围岩中的锆石，代表围岩混染等，反映了深部存在古老基底组分等，然后再小心求证，与其他地质资料结合，相互论证分析，并联系原始资料作者，交流、求证相应的捕获或继承锆石的结构特征及其他相关资料，为解释发现区域上的规律提供依据。
-
在详细的研究过程中，针对聚焦的科学问题，还应适当开展实测，开展重现性验证，并补充新的捕获/继承锆石微区数据，如前人可能缺少的同位素和微量元素等方面的数据，并对寄主岩石开展深入的岩石学研究，为深入研究捕获继承锆石成因和来源，以及与岩浆自结晶锆石信息对比，为揭示深部物质组成特征提供新依据。
-
1.4 新兴同位素示踪方法
-
1.4.1 基本原理和方法
-
传统的Sr-Nd-Hf等放射性同位素体系会随时间演化发生变化，在示踪年轻物质来源或者多期次叠加改造地质体的物质来源方面难度较大。相比之下，新兴的Ca、Mg、Fe、Ti、Zn等金属稳定同位素的分馏只受地质过程影响，不受时间演化影响。这一特性使其在深部物质示踪方面具有独特的优势。近年来，Ca、Mg、Fe、Ti、Zn等金属稳定同位素技术得到重视和发展，并应用于岩浆岩物源示踪分析。例如，Sun Jian et al.（2021）重新建立了精准高效的Ca同位素MC-ICP-MS测试技术，对全球不同地区的代表性火成碳酸岩样品开展了系统的Ca同位素地球化学研究，表明火成碳酸岩可以直接来源于正常的地幔，而不需要地表循环的碳酸盐物质参与。笔者在松潘甘孜三叠纪花岗岩带中，利用Mg、Li、B、Sr、Nd同位素分析，确定了沉积物源在不同类型花岗岩中所起的作用（Li Shan et al.，2021）。
-
通过对基性岩或中酸性岩浆岩全岩或单矿物的新兴金属稳定同位素分析测试，可以对岩浆岩源区的性质进行示踪制约。不同于放射性Sr-Nd-Hf同位素，可开展新兴稳定同位素体系的深部物质示踪，在厘定端元储库“参照系”的同位素组成的同时，还需要查明岩浆演化等过程的金属稳定同位素行为。因此，针对新兴同位素体系的深部物质示踪，需要构建专门的理论方法体系。目前，用于源区示踪的端元储库同位素组成“参照系”有待完善，岩浆演化等地质过程的新兴同位素分馏行为有待查明。
-
1.4.2 功能及优势
-
相较于Sr-Nd-Hf等放射性同位素体系在示踪源区新老物质组分方面的优势，新兴金属稳定同位素体系在深部物质示踪方面的主要功能包括:① 示踪深部的再循环地表物质。其基本原理是:地幔物质的Mg、Zn、Ca等同位素组成较均一，而地表碳酸盐物质与地幔物质的金属同位素组成有较大差异，再循环物质的加入会改变地幔物质的金属稳定同位素组成。例如，Li Shuguang et al.（2017）和Liu Sheng'ao et al.（2016）在中国东部发现大面积的新生代玄武岩Mg、Zn同位素“异常”，指示再循环碳酸盐加入到地幔。② 示踪深部源区的氧逸度状态。其基本原理是:变价金属元素（如Fe）的同位素分馏明显受控于氧逸度状态，不同氧逸度条件下部分熔融形成的熔体具有不同价态金属元素的比例（如Fe3+/Fe2+），因而具有不同的金属稳定同位素组成。例如，中国东部的部分碱性玄武岩和部分洋岛玄武岩的Fe同位素组成异常富集铁的重同位素，指示来源于氧逸度较高的源区（Konter et al.，2016; He Yongsheng et al.，2019）。③ 示踪长英质岩浆岩的构造背景属性。其基本原理是:不同构造背景（如地幔柱背景和岛弧背景）的岩浆体系由于性质不同，演化过程铁钛氧化物的结晶顺序不同，引起的Ti同位素分馏程度也不同，导致地幔柱背景和岛弧背景分异演化形成的中酸性岩浆岩Ti同位素组成明显不同（Deng Zhengbin et al.，2019）。由于新兴同位素在深部物质示踪方面目前仍处于探索阶段，其他功能有待进一步发掘。
-
不同的源区物质来源可能同时导致传统Sr-Nd-Hf同位素和新兴金属稳定同位素的协同变化。例如，Nebel et al.（2019）发现，Pitcairn洋岛玄武岩的Fe同位素和Sr-Nd-Pb同位素有很好的线性关系，认为该玄武岩的源区是两个端元的物质混合。笔者对山西的紫金山岩体开展了初步研究，也发现Fe同位素与Sr-Nd同位素有较好的相关性（孙剑等，未发表）。
-
1.4.3 存在问题
-
部分新兴同位素的分析测试技术仍不成熟，在数据的准确度、精度度及分析效率上需要进一步改进、完善。例如，不同学者对火成碳酸岩的Ca同位素分析数据结果差别很大，直接导致对火成碳酸岩来源认识上的巨大差异（Amsellem et al.，2020; Sun Jian et al.，2021; Banerjee et al.，2021）。
-
控制稳定同位素变化的因素不仅包括源区的同位素组成，而且岩浆-流体演化过程也可能产生较大的同位素分馏效应。
-
1.5 深部物质约束的地球物理方法
-
地球物理方法是地球深部探测的最主要和最常用的方法之一，主要包括重（重力）、磁（航磁、区域高精度重磁等）、电（激发极化、频谱激电、大地电磁）、震（天然地震层析成像、接收函数成像等; 人工源地震的深地震反射、宽角反射/折射地震等）。这些方法可以获得地球深部的物理性质和有关物性参数（O'Reilly et al.，2006; 吕庆田等，2015）。
-
1.5.1 主要方法
-
1.5.1.1 地震方法
-
利用人工地震和天然地震，为地球做“CT”是地球岩石圈深部探测的主要手段，目前已有长期的研究和成熟的技术。地震方法可以获取深部物质（如大陆岩石圈地幔）的密度和弹性性质在水平和垂向方向上的差异（Snyder et al.，2004）; 可以构建现今岩石圈的大范围物理结构（如速度、厚度）等模型。
-
近年来，噪声成像技术是最具革命性的地震成像技术，特别适合探测地壳上地幔顶部精细波速结构。与传统地震成像方法相比，背景噪声成像方法具有以下优点:① 场源丰富，不依赖于震源位置和震源机制信息，背景噪声场无处不在、无时不有; ② 拥有大量的空间采样方向; ③ 分辨率高，横向分辨率主要依赖于台站密度，地震台网密集区域，成像结果具有较高分辨能力，从而可以通过提升台站分布密度来提高成像质量; ④ 能够约束浅部速度结构，当噪声传播路径较短时，可以提取短周期（＜20 s）的面波及体波信号，对浅层结构具有更好的垂向分辨能力。最新噪声面波成像显示，准噶尔盆地中下地壳横波速度为高速异常，反映其基底铁镁质成分较高，可能为洋壳性质（鲍学伟项目组未发表资料）。
-
1.5.1.2 大地电磁
-
深部电导率主要受流体、熔体及电导率高的矿物（如金云母等）控制，而流体、熔体等的产生与构造演化及变形过程的关系十分密切（Xu Yixian et al.，2020）。大地电磁成像通过揭示深部电导率异常，可以为研究区域地质演化及构造变形过程提供较好的约束。例如，前人在西准噶尔进行了三维大地电磁成像，发现120~220 km深度范围存在一北东向的地幔高导体，主要由金云母等矿物导致，该高导体可能指示了该区在古生代大洋俯冲和大陆增生过程中残留的大洋板块（Xu Yixian et al.，2020）。
-
1.5.1.3 航磁
-
航磁指通过观测由岩石、矿体等探测对象的磁性差异而引起的磁异常特征，对探测区的磁性物质分布规律、地质构造等开展研究的一种探测方法。结合地质、岩矿石磁性资料，根据观测的磁异常特征判断引起磁异常的地质体性质，确定其空间展布特征，圈定断裂构造和侵入岩体，推测地质构造特征（朱小三等，2018; Zhu Xiaosan et al.，2019，2022a; Zhu Xiaosan and Lu Minjie，2021）。例如，区域航磁图显示，准噶尔盆地具有高的磁异常，这一特征也支持其深部基底为洋壳的认识。
-
1.5.1.4 重力
-
重力测量是利用组成地壳的各种地层、岩体、矿体的密度差异引起的重力变化而进行探测的一种地球物理方法。通过对重力数据进行校正和处理，消除各类误差，从叠加场中分离或突出目标地质体的异常场，使隐含在数据中的有用信息更易于被识别、比较，以及用于定性和定量解释。通过重力探测方法可以推断地区不同级别和规模的地质构造、划分构造单元，追索、圈定与围岩有明显差异的隐伏、半隐伏岩体或岩层，探测覆盖区基岩地质、沉积岩系各密度界面的起伏及沉积盆地范围，研究深大断裂的展布及火山结构，研究地壳、上地幔结构及地壳厚度变化（Zhu Xiaosan et al.，2022b）。例如，准噶尔盆地具有相对高的重力异常，这一特征与其深部基底为高密度的铁镁质洋壳较一致。此外，重力和地形的综合建模常用于揭示不同岩石圈弹性厚度的变化。
-
地球物理学家利用地震探测等方法，构建现今深部地壳的大范围物理结构（如速度、厚度）等模型，但空间分辨率通常较低，且难以对特定地点/小尺度的地壳结构进行高精度的约束。
-
1.5.2 实际应用
-
1.5.2.1 与实验模拟方法结合反演深部物质性质
-
这是最常用的研究应用。以往开展的不同尺度（岩石圈、造山带、区域）三维地质研究都是应用的这种方法。不同类型岩石的高温高压实验可以获取其密度、地震波速、泊松比、弹性模量等物性参数（Mavko et al.，2003）。但不同岩性的地震波速、密度等物性参数常存在一定程度的重叠，同时岩石的物理性质受到岩石矿物、化学组成及其所处深度和温度的影响，导致岩性识别困难。因此，这些参数与地球物理探测获取的物性参数融合是反演岩石圈架构的理想途径。
-
首先，通过地震波速、密度、热流、拉梅系数的物理参数的递进式分析，发现反演岩石组成的不确定性逐步降低，不同岩性的特征逐渐得到区分（杨立强等，待刊）。在弹性力学中，拉梅第一常数λ和拉梅第二常数μ是应变-应力关系中出现的两个与物质组成性质相关的参数（Feng Kang et al.，1981），其中拉梅第一常数λ没有对应的物理意义，但其对岩石性质的区分度较高（Ji Shaocheng et al.，2002）。而当拉梅系数乘以密度时，不同岩性之间存在明显的区分度，可以作为识别岩性的重要依据。如杨立强等（待刊）在胶东开展了很好地探索研究。综合利用波速、密度、热流等多种物性参数递进式约束岩石组成模型，对岩石物理实验结果进行温压校正，实现与地球物理探测的与岩石组成原位对标，为探索岩石圈组成架构提供新的思路与探索。
-
1.5.2.2 与同位素填图等的对应关系及联合示踪
-
由于同位素填图等深部物质示踪广泛开展，其结果与地球物理探测结果的对比研究成为新的探索方向。笔者在中亚北疆、青藏高原西南三江、扬子克拉通等地区，初步开展了岩石探针+多元同位素填图等多种地球物理探测结果对比研究，显示了很好的对应关系。
-
在北疆，Nd、Hf同位素、岩浆捕获锆石等多元同位素填图揭示的年轻地壳区域，特别是在西准噶尔，显示最年轻的区域。鲍学伟等（未发表资料）从200多个地震震台站数据获得高精度面波Vsv速度模型，包括不同深度切片的模型。研究结果显示，在30~40 km的模型与同位素填图揭示的物质分布型式一致，参数具有很好正相关关系，在西准噶尔地区显示为高速特征（最高值），对应岩浆岩揭示的30~40 km的源区（原来岩浆源区为40~50 km，目前岩体抬升剥蚀了约7 km）（图3）。这就找到了两种方法直接对比的同一个目标地质体。两种方法解释的物质属性也一致，即同位素填图和岩石探针解释为年轻的基性岩; Vsv揭示为镁铁质岩石为主的物质。除地球物理揭示的西准地壳高速特征外，该地区还具有高的电导率（高导体）、高的航磁异常（高磁）和高的布格重力异常，共同指示了基性岩下地壳-亏损的年轻上地幔。因此，联合反演30~40 km的深度为年轻的基性下地壳。具体的对应参数见图3。需要强调的是，这是少有的同一个深部、同一个圈层物质（下地壳）同位素填图等物质探测与地球物理探测直接对标（而不是不定深部的间接对标）的研究实例。
-
在北疆地区，综合上述研究，初步鉴别出以下几种地壳物质类型:① 年轻洋内弧物质+残留洋壳/新生底侵下地壳:该类物质主要分布在西准噶尔，具有高面波Vs，中高航磁异常、高布格重力异常; εHf（t）＞=+10，εNd（t）＞+4; 岩石类型上主要为镁铁质—超镁铁质斜长角闪岩-榴辉岩。② 微陆块中的古老基底:该类物质集中分布在伊犁地块南北两侧，具有中等面波Vs，低航磁异常、低布格重力异常; εHf（t）＜+5，εNd（t）＜0; 岩石成分上主要为闪长质-花岗质片麻岩。③ 微陆块内部再活化的年轻物质:该类物质集中分布在伊犁地块中部阿吾拉勒山—乌孙山一带，具有中高Vsv，强高航磁异常、中高布格重力异常; εHf（t）＞+5，εNd（t）＞0; 岩石成分上为镁铁质-超镁铁质。④ 变质沉积岩+变火成岩:该类物质集中分布在中国阿尔泰山，具有中等Vsv，中等航磁异常、中低布格重力异常; εHf（t）=0~+5; εNd（t）=-4~+4; 岩石成分上主要为闪长质-长英质片麻岩。
-
在拉萨地块，侯增谦等（2020）提出，拉萨地体的深反射地震剖面揭示（Lu Zhanwu et al.，2022）的20 km界面之下的中下地壳发育三个不同的反射区，即与南拉萨地体对应的透明弱反射区、与中萨拉地体对应的北倾弧形反射区、与北拉萨地体对应的伴有局部近水平的弱反射区。三者分别对应中酸性岩浆岩Hf同位素填图揭示的新生地壳块、中拉萨古老地壳块和北萨拉新生地壳块（图7a; Hou Zengqian et al.，2015）。这也是同一个深部的同一个圈层物质即下地壳的对应关系。
-
在三江地区，Hf同位素填图结果与地球物理Vs成像结果具有很好对应（Xu Bo et al.，2021a）:中酸性岩石的εHf＞0的区域，与地幔负Vs异常、Mg#≈89，共同约束年轻岩石圈地幔; εHf＜0区域，对应正Vs异常、Mg#＞89区域，揭示古老岩石圈。这是一种间接的对应。可以解释为中新生代软流圈物质（岩浆-热流）上涌，形成较为热的岩石圈地幔（负Vs），地幔物质（岩浆）上涌导致下地壳的年轻化。在冈底斯带，也具有这种对应关系（Xu Bo et al.，2021a）
-
综上，笔者提出应用地球物理探测和物质探测结合揭示深部物质属性可以总结为两个途径。一类是传统常用的，即首先通过对比地球物理手段获取的物性参数与地表出露的来自深部的岩石（捕虏体）实测的参数或实验模拟的参数，寻找对应性，从而以地球物理手段推测深部的岩石物性甚至物质性质，包括岩石类型等。这是以往应用地球物理和岩石探针结合探索深部物质属性的常用方法。该思路可应用于不同的尺度，大到岩石圈尺度，小到矿区的三维地质填图。例如，在华北克拉通（Gao shan et al.，2004; 高山等，2009）、秦岭-大别造山带（Gao Shan et al.，1992; 高山等，1999）、青藏高原（许志琴等，2006; 赵文津等，2008; 滕吉文等，2012，2019; 杨文采等，2015; Gao Rui et al.，2016; 王绪本等，2017; Lu Zhanwu et al.，2022）等关键地区都获得到了深部地壳、地幔结构。杨立强等（待刊）还进一步探索了该方向研究新的思路和方法。以上的方法需要依赖岩石物理实验结果或需要地表直接获得深部的样品（地幔包体、下地壳包体等）。但是，在很多地区，特别是造山带，很难获得直接样品信息和实验模拟资料。而岩石探针和同位素填图可以从物质示踪和填图的角度获得深部物质信息。这些资料如何与地球物理结合，是一个新的值得探索的问题和方向。
-
笔者初步研究显示，岩石探针和同位素填图结果与地球物理探测结果有很好的对应性。笔者总结为两类对应性。一类是直接对应，即中酸性侵入岩（主要来自地壳）同位素填图示踪的深度是目前30~40 km，即岩浆起源深度（一般40~50 km）-岩浆定位深度。而地球物理手段（如面波、航磁）反演的大约这个深部的结果就属于直接对应关系。另一类是间接对应，即中酸性侵入岩（主要来自地壳）同位素填图示踪的下地壳结果与地球物理探测的地幔-软流圈结果对应。例如，在青藏高原和三江，Hf同位素填图揭示的年轻地壳区域对应地震速度低速区，是软流圈物质上涌与地幔熔融导致下地壳熔融及其新生地壳形成（Xu Bo et al.，2021a）。
-
1.6 实验模拟技术约束
-
1.6.1 基本原理与方法
-
岩石部分熔融实验及其模拟计算是深入了解不同同位素体系间关系及是否解耦的重要途径，也是地球化学等物质示踪与地球物理示踪对标的桥梁。具体而言，一是开展直接的部分熔融实验，针对放射性同位素和稳定同位素，分别通过实测微量元素分配系数和同位素组成获取相关同位素分馏系数; 二是建立基于部分熔融实验的数据库，利用已有的相平衡关系和微量元素数据，针对放射性同位素体系开展模拟计算。
-
1.6.2 部分熔融实验
-
针对具体的地质样品开展直接的部分熔融实验来实测同位素的不平衡和分馏，无疑是最直接的观察手段，将为构建深部物质示踪理论方法体系提供有利支撑。然而，由于放射性同位素的特殊性，在实验的时间尺度内不会出现明显的同位素分馏，因此主要通过测定不同残留矿物-熔体之间相关元素（如Rb/Sr、Sm/Nd和Lu/Hf）的分配系数来预测熔体未来的放射性同位素比值的变化（Iles et al.，2018）。然而，近年来新兴的稳定同位素体系的分馏行为是可以实验测定的。例如，通过实验，可以模拟Fe同位素平衡分馏行为——选取洋中脊玄武岩（MORB）作为初始成分，利用立式筒式炉合成玄武岩玻璃和对应的金属线圈，并控制实验的条件和监测实时变化。利用MC-ICP-MS分别测定实验产物中硅酸盐玻璃和金属线圈的Fe同位素组成，通过数据拟合可以获得硅酸盐熔体中Fe3+和Fe2+之间的Fe同位素平衡分馏系数。该分馏系数有助于厘清地幔部分熔融过程中的Fe同位素分馏行为，进一步构建地幔橄榄岩源区和演化的铁同位素示踪体系，为示踪地幔岩浆运移提供了新思路（侯通等，未发表资料）。
-
1.6.3 实验数据模拟
-
针对放射性同位素体系，考虑到实验研究的局限性和已发表的大量实验数据，笔者建立了实验测定的各类岩浆体系的矿物-熔体元素（主要是Sm/Nd和Lu/Hf）分配系数数据库，用来预测部分熔融过程中放射性同位素行为。
-
具体来说，岩石圈内岩石发生部分熔融的过程中，由于熔体的快速形成和抽取，与残余矿物间扩散速度相对缓慢，往往造成熔体与残余体间难以达到化学和同位素平衡（Zhang et al.，2020）。部分熔融的源岩往往具有同位素不均一性，这是由于岩石形成时各个矿物相的Rb/Sr、Sm/Nd、Lu/Hf等同位素母子体初始比值不同，在源岩形成到部分熔融这段时间内，不同矿物相会演化出不同的同位素组成（Wolf et al.，2019）。同位素组成不均一的源岩部分熔融时，如果参与部分熔融反应的各反应物之间的比例与源岩各矿物相之间的比例不同（即非实比熔融），就会形成同位素组成与源岩迥异的熔体和残余体。因此，即使是同一源岩，当发生多阶段部分熔融（熔融反应不同）时，不同阶段的熔体携带残余体的相互混合，同样会造成花岗岩岩体中同位素的解耦现象（Zeng Lingsen et al.，2005a，2005b）。
-
通过汇总目前国内外的相关实验数据，建立数据库，可以通过数值建模来预测放射性同位素的系统性变化。这种方法有望揭示Nd、Hf同位素深部物质架构示踪结果与参数的关联性。未来可以以该数据库为基础，分别在增生造山带、碰撞造山带和克拉通三种环境下进行实际测试应用，并补充相关实验，为查明不同同位素体系间直接的关联性，以及与地球物理观察结果的对应性提供依据。特别是对于已知矿物组成的岩石及其所处的物理化学条件（或高温高压实验的产物），可以计算其波速及变化量。对于只知道全岩成分的岩石，可以采用PerpleX、pMelt、GeoPS等热力学软件对其组成展开计算模拟，获取源区岩石成分和熔融比例后，进一步计算其地球物理参数。
-
1.7 多元多维深部物质示踪数据集成分析
-
1.7.1 基本原理思路方法
-
由以上分析可以看出，深部物质示踪，特别是大区域的同位素填图、地球物理探测、实验模拟计算等都是建立在巨量的数据基础上，加之不同学科多维的数据分析和集成，需要统一的分析平台。因此，在大数据时代，建立各个学科数据库及其集成分析特别是智能化分析平台至关重要。建立数据库，实现数据共享，并搭建智能化分析平台是今后实现多尺度岩石圈三维物质架构分析的必要途径。
-
笔者在前期需求分析的基础上，开展了数据库概念模型设计，利用Powerdesigner16.5绘制多元同位素数据库概念模型（即E-R图），确定了数据库实体、数据集项、实体之间的关系及实体之间的集成关系。利用Powerdesigner提供的转换工具，依照E-R图，按选定的关系数据模型转换成对应的逻辑模型。在完成数据库逻辑模型的基础上，利用Powerdesigner转换工具，选定使用的数据库（MySQL 8.0）自动生成物理模型。经过上述三个步骤之间的依次迭代完善，构建较成熟的数据库物理模型，确定所有的表和列，定义外键用于确定表之间的关系，最终实现数据在数据库中的存放。
-
基于开源的Django web框架及bootstrap响应式页面布局工具，开发了项目数据共享平台。数据查询页面使用高德地图JavaScript API 2.0，初步实现了数据的框选查询（包括矩形、多边形、圆形）、显示数据点聚合、数据保存下载（Excel格式）等功能。同时，平台集成了开源Dtale可视化库，初步搭建了在线数据可视化应用，目前实现了基本功能，个性化需求还需要通过代码裁剪定制完成。D-Tale存储数据使用了全局变量，通常是最优的，数据在内存中不需要任何序列化/反序列化。但将D-Tale部署在多用户处理请求的web服务器上时，每个用户都有一个单独的python进程和一组单独的全局变量，需要用户配置Redis或者Shelve解决数据存储和全局状态问题。
-
1.7.2 若干实例
-
笔者尝试开展了北疆地球物理、地质和地球化学数据的集成分析。地震数据包括4列13020行数据（鲍学伟项目组未发表资料）; 地球化学数据集包括主量、微量和同位素地球化学（Huang He et al.，2020及黄河项目组未发表资料）。参与分析的数据有105列912行。为了将地震数据和地球化学数据融合，首先对地震数据集进行数学插值，为地球化学数据集增加深度为10 km和17.5 km对应的速度值两个数据列。在新生成的地球化学数据集上，进行模型训练，分类目标（Target）为构造单元，特征值（Feature）包括2列地震数据速度值及其他地球化学元素、年龄等数据的61列，训练的模型构造单元识别准确率为95%。笔者还进一步发现地震数据速度值对构造单元预测结果有较大影响，具体如何解释还需要收集更多数据做进一步分析。
-
这些初步分析结果显示，不同学科数据集中在同一个时空范围内，有望从不同学科、多种类型数据集成分析，获得更客观的认识。
-
2 深部物质架构的技术方法体系联合示踪的初步应用
-
笔者选取大陆岩石圈典型的三个大地构造单元，即中亚增生造山带（北疆）、青藏高原碰撞造山带（拉萨地块）和华北克拉通（破坏区、稳定区），开展了上述方法体系的初步应用，很好地联合反演和揭示了这些地区的三维岩石圈架构。
-
2.1 中亚增生造山带地区（北疆）三维岩石圈深部物质架构
-
2.1.1 北疆增生造山带岩石圈物质架构
-
应用上述深部物质架构的示踪方法体系，特别是区域多元同位素填图，揭示了北疆复杂的深部物质组成架构。阿尔泰深部物质组成具有中部较老、南侧较新的结构。西准噶尔深部地壳物质也以新生地壳组成为主（Nd模式年龄为TDM2为0.92~0.13 Ga; Hf模式年龄为TDM2为0.92~0.26 Ga）。在地壳生长方式上，阿尔泰造山带以陆壳水平生长为主，而准噶尔造山带以垂向生长为主。在构造划分依据上，额尔齐斯断裂带应为阿尔泰与准噶尔块体的构造单元分界线（Song Peng et al.，2019; 王涛等，2020）。
-
在新疆西天山及邻区，Hf同位素填图揭示了同一微陆块内部复杂的新老地壳分布，以及周期性的地壳生长/再循环模式（Huang He et al.，2020）。在伊犁地块南北缘靠近缝合带的区域，分布有古老地壳，其Hf模式年龄为1.4~1.0 Ga。而在块体中心阿吾拉勒山一带，大量发育年轻地壳，Hf模式年龄普遍低于1.0 Ga，局部低至0.6 Ga。从时间演化的角度，现有资料显示，伊犁地块北缘的岛弧岩浆作用在~350 Ma发生了从向陆迁移到向洋迁移的转换，即大洋俯冲（北天山洋）的方式从前进型俯冲转换为后撤型俯冲，且350 Ma至石炭纪末期的板片回卷/后撤伴随着显著的上覆岩石圈伸展、弧后盆地发育。在伊犁地块南缘，岩浆岩时空分布模式记录了在~420 Ma和~350 Ma两次古洋盆（南天山洋）俯冲方式从前进型到后撤型的转换，且泥盆纪末期—早石炭世的板片后撤同样导致弧后盆地的发育。
-
在北疆增生造山带应用多元同位素填图，全岩Nd和锆石Hf同位素体系得到了较一致的示踪结果（见前文1.2.2节，图2）。且同位素填图示踪的物质架构，与多种地球物理手段示踪的岩石圈物性结构与参数有着良好的对应关系（见前文1.5节，图4）。这进一步证实了同位素填图、地球物理探测等多种手段联合示踪岩石圈深部物质架构，在增生造山带具有一致性和有效性。
-
2.1.2 揭示覆盖区、盆地基底属性:以准噶尔盆地为例
-
准噶尔基底的性质一直存在争议。近期，G.Begg等国外研究者认为，准噶尔盆地基底为太古宙古老基底（Wang Jian et al.，2020; G.Begg在DEEP2021研讨会特邀报告）。笔者应用本文提出的方法体系对该区深部物质进行了联合示踪。东、西准噶尔（盆地两侧的造山带）岩浆岩Nd、Hf、Sr同位素填图和捕获锆石信息填图均未发现该区存在古老基底的信息。前已述及，基于200多个地震台站数据的面波数据，同样揭示30~40 km的深度为一高速体。在西准噶尔地区，30~40 km深度上面波示踪的结果与同位素填图示踪的年轻地壳在空间上高度耦合（图3）。基于搜集的地球物理资料，准噶尔盆地基本上都是强磁异常和重力高异常区，且在不同深度延拓的磁异常图和重力异常图（图5）上，准噶尔盆地的磁异常和重力异常特征基本变化很小（Zhu Xiaosan et al.，2022c）。综上所述，推测准噶尔盆地基底很可能为具高重、高磁特性的蛇绿混杂岩所组成（Yang Gaoxue et al.，2012）。此外，准噶尔盆地地壳厚度基本上在40 km左右，甚至更小（图5）。综上所述，应用多元同位素填图和多种地球物理方法联合示踪，结果一致揭示，准噶尔基底是由年轻的基性下地壳组成，而未示踪出大量的古老的基底。
-
图4 深部物质架构示踪方法体系联合反演的北疆岩石圈三维物质架构
-
Fig.4 3-D lithospheric material architecture of northern Xinjiang deciphered based on the methodological system in the present study
-
图5 新疆准噶尔—天山地区航磁、重力布格异常图及向上延拓图以及莫霍面深度图
-
Fig.5 Distribution of upward continuation processed aeromagnetic ΔT and Bouguer gravity anomalies with different heights, and Moho depth for the studied area of the Tianshan Mountains and Junggar basin in the Xinjiang
-
（a）~（c）—分别为延拓高度为5、30、50 km的航磁异常图;（d）~（f）—分别为延拓高度为5、40、60 km的重力布格异常图;（g）—莫霍面深度图
-
(a) ~ (c) —Magnetic anomalies with upward continuation processed maps with different heights of 5 km, 30 km and 50 km; (d) ~ (f) —bouguer gravity anomalies with upward continuation processed maps with different heights of 5 km, 40 km and 60 km; (g) —distribution of Moho depth
-
同时，笔者在盆地钻孔样品收集的Nd同位素也与造山带一致，佐证了盆地基底为年轻地壳，而不是古老物质。因此，这个实例也为解决覆盖区（盆地）基底的问题提供了技术方法。
-
2.2 青藏高原碰撞造山带（冈底斯-三江）三维岩石圈物质架构
-
在青藏高原东部的云南“三江”地区，锆石Hf同位素填图、幔源基性岩全岩地球化学成分、地球物理Vs成像和热化学模拟的联合示踪，揭示出古老与年轻地幔岩石圈水平并置的架构特征（Xu Bo et al.，2021a）。具有低εHf（t）值的前寒武纪岩石圈主要分布于腾冲-保山地体（εHf（t）=-7.2，=1624 Ma）、中咱-印度地体（εHf（t）=-11.5，=1979 Ma）以及中缅马苏地体（εHf（t）=-9.7，=1861 Ma），并显示出更高的地幔Mg#值（大于91）和地震波波速（VS正异常）。以高εHf（t）值为特征的年轻岩石圈则分布于金沙江-哀牢山缝合带附近的东羌塘-思茅地体（εHf（t）=-0.4，=1229 Ma），云南西部的碱性杂岩带（εHf（t）=+0.5，=1083 Ma），以及冕宁-德昌的碳酸岩-碱性杂岩带（εHf（t）=-0.3，=1150 Ma）。这些地区的岩石圈地幔Mg#≈89，同时具有更高的温度、熔融程度，以及更低的地震波波速（Vs负异常）。这些资料表明，中新生代软流圈物质的上涌，改造了上覆三江地区的岩石圈结构，形成了古老与年轻物质并置的架构特征。此外，年轻岩石圈的分布范围与新生代斑岩Cu-Mo-Au和与碳酸岩相关的稀土矿床的分布相吻合，证明岩石圈深部物质组成结构的探测对矿产勘探具有重要的指导意义。
-
穿越青藏高原南部的深反射地震探测精细地揭示了拉萨地体地壳的结构（Lu Zhanwu et al.，2022），其结果与岩浆岩锆石Hf同位素填图有着很好的对应关系（Hou Zengqian et al.，2015; 侯增谦等，2020）。锆石Hf同位素填图揭示出的以年轻地壳物质为主的南拉萨地体与透明弱反射区对应，古老的中拉萨地体则与北倾弧形反射区，而年轻的北拉萨地体亦对应的伴有局部近水平的弱反射区（图6、图7），揭示出三个亚地体之间为逆冲叠覆的空间关系（侯增谦等，2020）。结合区域岩浆岩的岩石学研究可知，北拉萨地体的新生地壳通常被认为与中生代班公湖-怒江洋双向俯冲引起的大规模弧岩浆作用密切相关，因此该地体的逆冲叠置可能发生自新生代的印亚大陆碰撞以来; 而南拉萨地体的新生地壳则由新特提斯洋俯冲阶段和印亚大陆同碰撞阶段的幔源岩浆持续注入而成，表明其逆冲叠覆可能发生于之后的晚碰撞期（侯增谦等，2020）。
-
图6 深部物质架构示踪方法体系联合反演的青藏高原岩石圈三维物质架构
-
Fig.6 3-D lithospheric material architecture of Tibet Plateau deciphered based on the methodological system in the present study
-
图7 拉萨地体Hf同位素填图结果（a）与基于反射深地震资料的综合地质解释（b）（据侯增谦等，2020; Lu Zhanwu et al.，2022）
-
Fig.7 Hf isotope mapping result (a) and possible geological interpretation based on reflect seismic profile across the southern Tibet (b) (modified after Hou Zengqian et al., 2020; Lu Zhanwu et al., 2022)
-
2.3 华北克拉通三维岩石圈物质架构
-
2.3.1 华北克拉通西部稳定区岩石圈物质架构
-
华北克拉通是世界上最古老的克拉通之一，在地表和深部都存在大于3.8 Ga的地壳物质（Liu Dunyi et al.，1992; Ma Qiang et al.，2020）。其广泛出露的前寒武纪变质基底是探索克拉通深部物质组成架构及壳-幔物质循环至关重要的研究对象。通过对华北克拉通中部恒山杂岩体的一系列太古宙—古元古代石榴子石基性麻粒岩、斜长角闪岩及围岩TTG片麻岩开展全岩主量、微量元素、矿物化学和锆石U-Pb-Hf同位素研究，可以约束华北克拉通古老变质基底的变质历史、重建中-下地壳结构。地壳结构从上至下依次为:中地壳（约19 km）的片麻质岩石，下地壳上部（约19~24 km）的角闪岩相岩石，下地壳中部（约24~27 km）二辉麻粒岩和下部（约46 km）的石榴子石基性麻粒岩（图8）。层状地壳受到古元古代热事件的强烈改造，同时显示随着深度的增加，地壳愈加年轻，而成分更亏损的趋势。这种趋势在华北克拉通北缘、南缘及东缘等地区都有所展现。
-
2.3.2 华北克拉通东部破坏区岩石圈物质架构
-
与全球古老克拉通相比，华北克拉通遭受了最明显的破坏（Carson et al.，2005; 朱日祥等，2011），这与华北克拉通所处的构造位置和周边板块的多期俯冲密不可分。晚中生代古太平洋俯冲板块的回转、海沟后撤引起的大陆岩石圈强烈伸展是华北克拉通岩石圈巨厚减薄的主导因素（Tang Yanjie et al.，2013; 朱日祥等，2019）。因此，有效识别周缘多向俯冲对克拉通岩石圈改造的时空-物质组成架构是揭示破坏区岩石圈组成架构的关键。
-
课题组选择处于华北克拉通破坏区的胶东半岛开展了Sr-Nd-Hf-O-Pb多元同位素填图。以中酸性岩代表地壳，基性岩代表地幔，开展了不同深源岩浆岩的Nd-Hf同位素填图（图9）。研究结果显示，胶东半岛地壳与地幔整体耦合，但是在胶莱盆地局部出现壳-幔同位素的解耦，这可能与其强烈的幔源岩浆底侵改造有关，显示胶莱盆地可能受到更强烈的破坏作用。而中酸性岩浆岩的Sr-Nd-Hf-O-Pb多元同位素填图（图9）结果显示，胶莱盆地相对于胶北隆起与苏鲁造山带有更多的新生地壳物质组成，显示了不均一的地壳组成; 锆石Hf-O与全岩Nd-Pb同位素结果基本一致。而胶东半岛地球物理反演获取的岩石圈物质架构显示，胶北隆起和胶莱盆地总体具有克拉通岩石圈层块特征，中地壳之下被强烈改造; 胶北隆起中-下地壳显示的低速体可能与胶西北中生代大规模岩浆-热液-金成矿事件有关; 而苏鲁造山带相对胶北隆起和胶莱盆地，岩石圈结构更复杂，显示被强烈改造的古老造山带物质组成架构。这与同位素填图的结果具有很好的对应性，揭示了同位素填图与地球物理探测对于探测岩石圈架构的潜力。
-
图8 华北克拉通恒山地区地壳结构组成模型图（据Zhang Xiahui et al.，2021）
-
Fig.8 Crustal architecture of the Hengshan region of North China Craton (modified after Zhang Xiahui et al., 2021)
-
图9 基于深部物质架构示踪方法体系联合反演胶东半岛岩石圈三维物质架构（据杨立强等，待刊）
-
Fig.9 3-D lithospheric architecture of Jiaodong Peninsula deciphered based multi-methods of the methodological system introduced in the present study (after Yang Liqiang et al., for publication)
-
3 小结
-
长期以来，深部探测主要依赖地球物理和钻探而缺少岩石圈深部物质探测和研究技术方法。对深部物质的了解主要通过两种途径:① 地球物理推测方法，即依据地表出露的来自深部的某些岩石物理参数（包括温压条件和年龄信息），对比地球物理深部探测结果，解释或推测深部相应参数的物质特征。这是最常见的方法，如不同比例尺的三维地质填图。② 捕虏体方法，即利用岩浆岩深源岩石捕虏体，结合地球物理（地震、大地电磁等），开展岩石圈深部探测，即3D、4D岩石圈填图。
-
除此之外，笔者重点探索第三种途径，即充分利用地表出露的来自不同深部的岩浆岩，排除其他因素，通过岩石探针和同位素填图，示踪深部物源（物质）特征，特别是通过填图，由点到面，了解深部物质在横向上（水平）、纵向上（垂向）的分布特征，即三维物质架构，甚至四维演变。在此基础上，综合上述三种途径，总结、构建较完整的一套以岩石探针和同位素填图为核心的揭示岩石圈三维物质组成架构的方法体系。这为加强深部研究中的物质探测提供了思路和技术方法。
-
该方法体系的构建取得了初步成效。例如，总结出多元同位素填图的关联性; 不同地区、不同岩石类型的Nd、Hf同位素显示正相关关系。这些研究结果，一方面证实了多元同位素体系在探测深部物质方面的结果具有一致性和有效性; 另一方面，也显示多元同位素填图的结果可以相互校正、互相标定，开展联合示踪填图。同位素填图结果与地球物理探测结果具有很好的对应性，在局部地区还总结出同位素参数和地球物理参数之间半定量的对应关系。而实验模拟也为了解不同同位素体系的关联性和与地球物理结果的对应性提供了理论依据。
-
该方法体系在应用方面也取得一些重要成果。例如，中亚增生造山带、青藏高原碰撞造山带和华北-扬子克拉通等地区的研究结果显示，通过岩石探针和Sr、Nd、Hf、Pb等多元同位素填图，能很好地揭示这些地区的岩石圈深部物质组成架构。同时，还能解决一些重大问题。例如，准噶尔盆地基底为年轻基性下地壳（或洋壳为主），为年轻的增生造山带。青藏高原的多元同位素填图与多种方法地球物理探测进一步证实，印度板块俯冲及板片撕裂是导致巨型斑岩成矿带发育的驱动力（Hou Zengqian et al.，submitted）。
-
上述研究显示，区域岩浆岩同位素填图与地球物理等技术的结合在探索深部物质架构方面具有巨大潜力，岩石探针及多元同位素填图与地球物理和实验模拟等方法的结合，有望建立深部物质探测技术方法体系，这将有助于全面开展和推动深部探测与深部物质与演变及其动力学研究，其是今后深地重要探索方向。
-
致谢:谨以此文祝贺中国地质学会成立一百周年。感谢莫宣学、金振民、高锐、吴福元、杨经绥、徐义刚、张宏福、李献华、肖文交、邓军院士及董树文、韩宝福、吕庆田教授（研究员）的指导和有益讨论。匿名审稿人为本文修改提出了宝贵意见。本文写作过程中和文言、高雪、王偲瑞、吴欢欢、潘蓓蓓、李海舟等提供了相关协助，在此一并致谢。
-
参考文献
-
Afonso J C, Ben Mansour W, O'Reilly S Y, Griffin W L, Salajeghegh F, Foley S, Begg G, Selway K, Macdonald A, Januszczak N, Fomin I, Nyblade A A, Yang Yingjie. 2022. Thermochemical structure and evolution of cratonic lithosphere in central and southern Africa. Nature Geoscience, 15: 405~410.
-
Amsellem E, Moynier F, Bertrand H, Bouyon A, Mata J, Tappe S, Day J. 2020. Calcium isotopic evidence for the mantle sources of carbonatites. Science Advances, 6(23): a3269.
-
Banerjee A, Chakrabarti R, Simonetti A. 2021. Temporal evolution of δ44/40Ca and 87Sr/86Sr of carbonatites: implications for crustal recycling through time. Geochimica et Cosmochimica Acta, 307: 168~191.
-
Begg G C, Griffin W L, Natapov L M, O'Reilly S Y, Grand S P, O'Neill C J, Hronsky J M A, Poudjom D Y, Swain C J, Deen T, Bowden P. 2009. The lithospheric architecture of Africa: seismic tomography, mantle petrology, and tectonic evolution. Boulder, CO: Geological Society of America, 23~50.
-
Belousova E, Griffin W, O'Reilly S Y, Fisher N. 2002. Igneous zircon: trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143(5): 602~622.
-
Belousova E A, Kostitsyn Y A, Griffin W L, O'Reilly S Y. 2009. Testing models for continental crustal growth: a TerraneChron approach. Geochimica et Cosmochimica Acta, 73(13S): A107.
-
Bizimis M, Sen G, Salters V J. 2004. Hf-Nd isotope decoupling in the oceanic lithosphere: constraints from spinel peridotites from Oahu, Hawaii. Earth and Planetary Science Letters, 217(1-2): 43~58.
-
Buys J, Spandler C, Holm R J, Richards S W. 2014. Remnants of ancient Australia in Vanuatu: implications for crustal evolution in island arcs and tectonic development of the Southwest Pacific. Geology (Boulder), 42(11): 939~942.
-
Carlson R W, Pearson D G, James D E. 2005. Physical, chemical, and chronological characteristics of continental mantle. Reviews of Geophysics, 43(1): G1001.
-
Champion D C, Huston D L. 2016. Radiogenic isotopes, ore deposits and metallogenic terranes: novel approaches based on regional isotopic maps and the mineral systems concept. Ore Geology Reviews, 76: 229~256.
-
Chen Bin, Arakawa Y. 2005. Elemental and Nd-Sr isotopic geochemistry of granitoids from the West Junggar foldbelt(NW China), with implications for Phanerozoic continental growth. Geochimica et Cosmochimica Acta, 69(5): 1307~1320.
-
Chi Jishang. 1988. Study of Cenozoic Basalts and Upper Mantles in Eastern China. Wuhan: China University of Geosciences Press (in Chinese).
-
Choi S H, Mukasa S B. 2012. Lu-Hf and Sm-Nd isotope systematics of Korean spinel peridotites: a case for metasomatically induced Nd-Hf decoupling. Lithos, 154: 263~276.
-
Condie K C. 1999. Mafic crustal xenoliths and the origin of the lower continental crust. Lithos, 46(1): 95~101.
-
Deng Jinfu. 1994. Petrophysicochemistry and petrophysics. Earth Science Frontiers, (Z1): 57~63 (in Chinese with English abstract).
-
Deng Jinfu, Su Shangguo, Zhao Hailing, Mo Xuanxue, Xiao Qinghui, Zhou Su, Liu Gui, Zhao Guochun. 2003. Deep processes of Mesozoic Yanshanian lithosphere thinning in North China. Earth Science Frontiers, 10(3): 41~50 (in Chinese with English abstract).
-
Deng Zhengbin, Chaussidon M, Savage P, Robert F, Pik R, Moynier F. 2019. Titanium isotopes as a tracer for the plume or island arc affinity of felsic rocks. Proceedings of the National Academy of Sciences, 116(4): 1132~1135.
-
Ding Huixia, Zhang Zeming, Kohn J M, Gou Zhengbin. 2021. Timescales of partial melting and melt crystallization in the eastern Himalayan orogen: insights from zircon petrochronology. Geochemistry, Geophysics, Geosystems, 22(4): e2020GC009539.
-
Dong Shuwen, Li Tingdong, Chen Xuanhua, Wei Wenbo, Gao Rui, Lv Qingtian, Yang Jingsui, Wang Xueqiu, Chen Qunce, Shi Yaolin, Huang Danian, Zhou Qi. 2012. Progress of deep exploration in mainland China: a review. Chinese Journal of Geophysics, 55(12): 3884~3901 (in Chinese with English abstract).
-
Dong Shuwen, Gao Rui, Yin An, Guo Tonglou, Zhang Yueqiao, Hu Jianmin, Li Jianhua, Shi Wei, Li Qiusheng. 2013. What drove continued continent-continent convergence after ocean closure? Insights from high-resolution seismic-reflection profiling across the Daba Shan in central China. Geology, 41(6): 671~674.
-
Dong Shuwen, Li Tingdong, Chen Xuanhua, Gao Rui, Lv Qingtian, Shi Yaolin, Huang Danian, Yang Jingsui, Wang Xueqiu, Wei Wenbo, Chen Qunce. 2014. SinaProbe revealed crustal structures, deep processes, and metallogenic background within China continent. Earth Science Frontiers, 21(3): 201~225 (in Chinese with English abstract).
-
Feng Kang, Shi Zhongci. 1981. Mathematical Theory of Elastic Structures. Beijing: Science Press (in Chinese).
-
Gao Rui, Lu Zhanwu, Klemperer S L, Wang Haiyan, Dong Shuwen, Li Wenhui, Li Hongqiang. 2016. Crustal-scale duplexing beneath the Yarlung Zangbo suture in the western Himalaya. Nature Geosciences, 9(7): 555~560.
-
Gao Shan, Zhang Benren, Luo T C, Li Z J, Xie Qian Li, Gu Xiaoming, Zhang H F, Ouyang Jianping, Wang Dongpo, Gao Changling. 1992. Chemical composition of the continental crust in the Qinling orogenic belt and its adjacent North China and Yangtze cratons. Geochimica et Cosmochimica Acta, 56(11): 3933~3950.
-
Gao Shan, Zhang Benren, Jin Zhenmin. 1999. Delamination of the lower crust of the Qinling-Dabie orogenic belt. Scientia Sinica (Terrae), 29(6): 532~541 (in Chinese).
-
Gao Shan, Rudnick R L, Yuan Hongling, Liu Xiaoming, Liu Yongsheng, Xu Wenliang, Ling Wenli, John Ayers, Wang Xuanche, Wang Qinghai. 2004. Recycling lower continental crust in the North China Craton. Nature, 432: 892~897.
-
Gao Shan, Zhang Junfeng, Xu Wenliang, Liu Yongsheng. 2009. Delamination and destruction of the North China Craton. Chinese Science Bulletin, 54(14): 1962~1973 (in Chinese).
-
Garçon M, Chauvel C, France L C, Huyghe P, Lavé J. 2013. Continental sedimentary processes decouple Nd and Hf isotopes. Geochimica et Cosmochimica Acta, 121: 177~195.
-
Geng Hongyan, Sun Min, Yuan Chao, Xiao Wenjiao, Xian Weisheng, Zhao Guochun, Zhang Lifei, Wong Kenny, Wu Fuyuan. 2009. Geochemical, Sr-Nd and zircon U-Pb-Hf isotopic studies of Late Carboniferous magmatism in the West Junggar, Xinjiang: implications for ridge subduction? Chemical Geology, 266(3-4): 364~389.
-
Griffin B, Andi Z, O'Reilly S, Ryan C G. 1998. Phanerozoic evolution of the lithosphere beneath the Sino-Korean craton. In: Flower M F J, Chung S L, Lo C H, Lee T Y, eds. Mantle Dynamics and Plate Interactions in East Asia, 27: 107~126.
-
Griffin W L, O'Reilly S Y, Doyle B J, Pearson N J, Coopersmith H, Kivi K, Malkovets V, Pokhilenko N. 2004. Lithosphere mapping beneath the North American plate. Lithos, 77(1-4): 873~922.
-
Griffin W L, Powell W J, Pearson N J, O'Reilly S Y. 2008. GLITTER: data reduction software for laser ablation ICP-MS. In: Sylvester P, ed. Laser Ablation-ICP-MS in the Earth Sciences. Mineralogical Association of Canada Short Course Series, 40: 204~207.
-
Griffin W L, O'reilly S Y, Afonso J C, Begg G C. 2009. The composition and evolution of lithospheric mantle: a re-evaluation and its tectonic implications. Journal of Petrology, 50(7): 1185~1204.
-
Griffin W L, Begg G C, O'Reilly S Y. 2013. Continental-root control on the genesis of magmatic ore deposits. Nature Geoscience, 6: 905~910.
-
Guitreau M, Blichert Toft J, Mojzsis S J, Roth A S, Bourdon B, Cates N L, Bleeker W. 2014. Lu-Hf isotope systematics of the hadean-Eoarchean Acasta gneiss complex(northwest territories, Canada). Geochimica et Cosmochimica Acta, 135: 251~269.
-
He Yongsheng, Meng Xunan, Ke Shan, Wu Hongjie, Zhu Chuanwei, Teng Fangzhen, Hoefs J, Huang Jian, Yang Wei, Xu Lijuan, Hou Zhenhui, Ren Zhongyuan, Li Shuguang. 2019. A nephelinitic component with unusual δ56Fe in Cenozoic basalts from eastern China and its implications for deep oxygen cycle. Earth and Planetary Science Letters, 512: 175~183.
-
Hong Dawei, Zhang Jisheng, Wang Tao, Wang Shiguang, Xie Xilin. 2004. Continental crustal growth and the supercontinental cycle: evidence from the Central Asian Orogenic Belt. Journal of Asian Earth Sciences, 23(5): 799~813.
-
Hoskin P, Ireland T. 2000. Rare earth element chemistry of zircon and its use as a provenance indicator. Geology, 28(1): 627~630.
-
Hou Zengqian, Duan Lianfeng, Lu Yongjun, Zheng Yuanchuan, Zhu Dicheng, Yang Zhiming, Yang Zhusen, Wang Baodi, Pei Yingru, Zhao Zhidan, Mccuaig T Campbell. 2015. Lithospheric architecture of the Lhasa Terrane and its control on ore deposits in the Himalayan-Tibetan Orogen. Economic Geology and the Bulletin of the Society of Economic Geologists, 110(6): 1541~1575.
-
Hou Zengqian, Wang Tao. 2018. Isotopic mapping and deep material probing(Ⅱ): imaging crustal architecture and its control on mineral systems. Earth Science Frontiers, 25(6): 20~41 (in Chinese with English abstract).
-
Hou Zengqian, Zheng Yuanchuan, Lu Zhanwu, Xu Bo, Wang Changming, Zhang Hongrui. 2020. Growth, thickening and evolution of the thickened crust of the Tibet Plateau. Acta Geological Sinica, 94(10): 2797~2815 (in Chinese with English abstract).
-
Hou Zengqian, Wang Rui, Zhang Haijiang, Zheng Yuanchuan, Jin Sheng, Weinberg Roberto, Xu Bo, Yang Zhiming, Hao Aowei, Gao Lei, Zhang Letian. Tearing of subducted Indian plate drives giant copper deposit formation in Tibet. (Submitted).
-
Huang He, Wang Tao, Tong Ying, Qin Qie, Ma Xuxuan, Yin Jiyuan. 2020. Rejuvenation of ancient micro-continents during accretionary orogenesis: insights from the Yili block and adjacent regions of the SW Central Asian Orogenic Belt. Earth-Science Reviews, 208: 103255.
-
Iles K A, Hergt J M, Woodhead J D. 2018. Modelling isotopic responses to disequilibrium melting in granitic systems. Journal of Petrology, 59(1): 87~114.
-
Jenkyns H C, Forster A, Schouten S, Sinninghe D J S. 2004. High temperatures in the Late Cretaceous Arctic Ocean. Nature, 432(7019): 888~892.
-
Jeon H, Williams I S. 2018. Trace inheritance—clarifying the zircon O-Hf isotopic fingerprint of I-type granite sources: implications for the restite model. Chemical Geology, 476: 456~468.
-
Ji Shaocheng, Wang Qin, Xia Bin. 2002. Handbook of Seismic Properties of Minerals, Rocks and Ores. Montreal Que, Canada: International Press.
-
Jin Zhenmin, Ou Xingong, Xu Haijun, Kern H, Jin Shuyan, Hou Guangjiu. 2004. Elastic wave velocities of 2000m depth at Chinese Continental Scientific: constraints on deep seismic reflection. Acta Petrologica Sinica, 20(1): 81~96 (in Chinese with English abstract).
-
Kemp A, Hawkesworth C, Foster G, Paterson B, Woodhead J, Hergt J, Gray C, Whitehouse M. 2007. Magmatic and crustal differentiation history of granitic rocks from Hf-O isotopes in zircon. Science, 315: 980~983.
-
Kilgore M L, Peslier A H, Brandon A D, Schaffer L A, Morris R V, Graff T G, Agresti D G, O'Reilly S Y, Griffin W L, D. Pearson G, Gangi K, Shaulis B J. 2020. Metasomatic control of hydrogen contents in the layered cratonic mantle lithosphere sampled by Lac de Gras xenoliths in the central Slave craton, Canada. Geochimica et Cosmochimica Acta, 286: 29~53.
-
King E M, Valley J W, Davis D W, Edwards G R. 1998. Oxygen isotope ratios of Archean plutonic zircons from granite-greenstone belts of the Superior province: indicator of magmatic source. Precambrian Research, 92(4): 365~387.
-
Kong Xiangye, Zhang Chen, Liu Dongdong, Shu Jiang, Luo Qun, Zeng Jianhui, Liu Luofu, Luo Long, Shao Hengbo, Liu Dan, Liu Xiaoyu, Wang Xinpeng. 2018. Disequilibrium partial melting of metasediments in subduction zones: evidence from O-Nd-Hf isotopes and trace elements in S-type granites of the Chinese Altai. Lithosphere, 11(1): 149~168.
-
Konter J G, Pietruszka A J, Hanan B B, Finlayson V A, Craddock P R, Jackson M G, Dauphas N. 2016. Unusual δ56Fe values in Samoan rejuvenated lavas generated in the mantle. Earth and Planetary Science Letters, 450: 221~232.
-
Lee R G, Byrne K, Michael A D, Hart C J R, Hollings P, Gleeson S A, Alfaro M. 2021. Using zircon trace element composition to assess porphyry copper potential of the Guichon Creek batholith and Highland Valley Copper deposit, south-central British Columbia. Mineralium Deposita, 56(2): 215~238.
-
Li Shan, Miller C F, Wang Tao, Xiao Wenjiao, Chew D. 2021. Role of sediment in generating contemporaneous, diverse“type”granitoid magmas. Geology, 50(4): 427~431.
-
Li Shuguang, Yang Wei, Ke Shan, Meng Xunan, Tian Hengci, Xu Lijuan, He Yongsheng, Huang Jian, Wang Xuance, Xia Qunke, Sun Weidong, Yang Xiaoyong, Ren Zhongyuan, Wei Haiquan, Liu Yongsheng, Meng Fancong, Yan Jun. 2017. Deep carbon cycles constrained by a large-scale mantle Mg isotope anomaly in eastern China. National Science Review, 4(1): 111~120.
-
Liptai N, Patkó L, Kovács I J, Hidas K, Pintér Z, Jeffries T, Szabó C. 2017. Multiple metasomatism beneath the Nógrád-Gömör Volcanic Field(northern Pannonian basin)revealed by upper mantle peridotite xenoliths. Journal of Petrology, 58(6): 1107~1144.
-
Liu D Y, Nutman A P, Compston W, Wu J S, Shen Q H. 1992. Remnants of ≥3800 Ma crust in the Chinese part of the Sino-Korean craton. Geology, 20(4): 339~342.
-
Liu Dong, Zhao Zhidan, Zhu Dicheng, Niu Yaoling, Mark Harrison T. 2014. Zircon xenocrysts in Tibetan ultrapotassic magmas: Imaging the deep crust through time. Geology (Boulder), 42(1): 43~46.
-
Liu Guangding. 2002. Review and prospect—21st century solid geophysics. Progress in Geophysics, 17(2): 191~197 (in Chinese with English abstract).
-
Liu Sheng'ao, Wang Zezhou, Li Shuguang, Huang Jian, Yang Wei. 2016. Zinc isotope evidence for a large-scale carbonated mantle beneath eastern China. Earth and Planetary Science Letters, 444: 169~178.
-
Liu Zhichao, Liu Xiaochi, Yu Liangjun, Wang Jiangang. 2020. Highly fractionated origin and magmtic-hydrothermal evolution of the Kampa leucogranites in the Tethyan Himalaya. Journal of Nanjing University(Natural Science), 56(6): 800~814 (in Chinese with English abstract).
-
Lu Yongjun, Wingate M T, Champion D C, Smithies R H, Johnson S P, Mole D R, Poujol M, Zhao J, Maas R, Creaser R A. 2021a. Samarium-neodymium isotope map of Western Australia. Geological Survey of Western Australia, digital data layer.
-
Lu Yongjun, Wingate M T, Romano S S, Mole D R, Kirkland C L, Kemp A I S, Belousova E A, Smithies R H, Gessner K, Johnson S P. 2021b. Zircon lutetium-hafnium isotope map of Western Australia. Record-Geological Survey of Western Australia, 13~16.
-
Lu Zhanwu, Guo Xiaoyu, Gao Rui, Murphy M A, Huang Xingfu, Xu Xiao, Li Sanzhong, Li Wenhui, Zhao Junmeng, Li Chunsen, Xiang Bo. 2022. Active construction of southernmost Tibet revealed by deep seismic imaging. Nature Communicatios, 13: 3143.
-
Luo Qun, Zhang Chen, Jiang Shu, Liu Luofu, Liu Dongdong, Kong Xiangye, Liu Xiaoyu, Wang Xinpeng. 2019. Partial melting of oceanic sediments in subduction zones and its contribution to the petrogenesis of peraluminous granites in the Chinese Altai. Geological Magazine, 156(4): 585~604.
-
Lv Qingtian, Liu Zhendong, Dong Shuwen, Yan Jiayong, Zhang Yongqian. 2015. The nature of Yangtze River deep fault zone: evidence from deep seismic data. Chinese Journal Geophysics, 58(12): 4344~4359 (in Chinese with English abstract).
-
Ma Qiang, Xu Yigang, Huang Xiaolong, Zheng Jianping, Ping Xianquan, Xia Xiaoping. 2020. Eoarchean to Paleoproterozoic crustal evolution in the North China Craton: evidence from U-Pb and Hf-O isotopes of zircons from deep-crustal xenoliths. Geochimica et Cosmochimica Acta, 278: 94~109.
-
Mavko G, Mukerji T, Dvorkin J. 2003. The Rock Physics Handbook. Cambridge: Cambridge University Press.
-
Miller J S, Matzel J E, Miller C F, Burgess S D, Miller R B. 2007. Zircon growth and recycling during the assembly of large, composite arc plutons. Journal of Volcanology and Geothermal Research, 167(1-4): 282~299.
-
Mo Xuanxue. 2011. Magma and magamatic/igneous rocks: a lithoprobe into the deep earth and records of the earth's evolution. Chinese Journal of Nature, 33(5): 255~269 (in Chinese with English abstract).
-
Mo Xuanxue. 2019. Magmatism and deep geological process. Earth Science, 44(5): 1487~1493 (in Chinese with English abstract).
-
Mole D R, Fiorentini M L, Thebaud N, Cassidy K F, Mccuaig T C, Kirkland C L, Romano S S, Doublier M P, Belousova E A, Barnes S J, Miller J. 2014. Archean komatiite volcanism controlled by the evolution of early continents. Proceedings of the National Academy of Sciences, 111(28): 10083~10088.
-
Nebel O, Sossi P A, Bénard A, Arculus R J, Yaxley G M, Woodhead J D, Rhodri D D, Ruttor S. 2019. Reconciling petrological and isotopic mixing mechanisms in the Pitcairn mantle plume using stable Fe isotopes. Earth and Planetary Science Letters, 521: 60~67.
-
Nowell G M, Pearson D G, Bell D R, Carlson R W, Smith C B, Kempton P D, Noble S R. 2004. Hf isotope systematics of kimberlites and their megacrysts: new constraints on their source regions. Journal of Petrology, 45(8): 1583~1612.
-
O'Reilly S Y, Griffin W L. 1996. 4-D lithosphere mapping: methodology and examples. Tectonophysics, 262(1-4): 3~18.
-
O'Reilly S Y, Griffin W L, Djomani Y H P, Morgan P. 2001. Are lithospheres forever? Tracking changes in subcontinental lithospheric mantle through time. GSA Today, 11(4): 4~10.
-
O'Reilly S Y, Griffin W L. 2006. Imaging global chemical and thermal heterogeneity in the subcontinental lithospheric mantle with garnets and xenoliths: geophysical implications. Tectonophysics, 416(1-4): 289~309.
-
Özaydin S, Selway K, Griffin W L. 2021. Are xenoliths from southwestern Kaapvaal Craton representative of the broader mantle? Constraints from magnetotelluric modeling. Geophysical Research Letters, 48(11): e2021GL092570.
-
Özaydin S, Selway K, Griffin W L, Moorkamp M. 2022. Probing the southern African lithosphere with magnetotellurics: 2. Linking electrical conductivity, composition, and tectonomagmatic evolution. Journal of Geophysical Research: Solid Earth, 127(3): e2021JB023105.
-
Patchett P J, Vervoort J D, Söderlund U, Salters V J. 2004. Lu-Hf and Sm-Nd isotopic systematics in chondrites and their constraints on the Lu-Hf properties of the Earth. Earth and Planetary Science Letters, 222(1): 29~41.
-
Qi Nan, Zhou Zhenju, Zou Haibo, Chen Yanjing. 2021. 3. 85 Ga continental crust beneath the southern North China Craton: evidence from zircon xenocrysts in Cretaceous granites. Gondwana Research, 91: 277~285.
-
Saunders J E, Pearson N J, O'Reilly S Y, Griffin W L. 2015. Sulfide metasomatism and the mobility of gold in the lithospheric mantle. Chemical Geology, 410: 149~161.
-
Schmitz M D, Vervoort J D, Bowring S A, Patchett P J. 2004. Decoupling of the Lu-Hf and Sm-Nd isotope systems during the evolution of granulitic lower crust beneath southern Africa. Geology, 32(5): 405~408.
-
Siebel W, Schmitt A K, Danišík M, Chen Fukun, Meier S, Weiβ S, Eroğlu S. 2009. Prolonged mantle residence of zircon xenocrysts from the western Eger rift. Nature Geoscience, 2(12): 886~890.
-
Smyth H R, Hamilton P J, Hall R, Kinny P D. 2007. The deep crust beneath island arcs: inherited zircons reveal a Gondwana continental fragment beneath East Java, Indonesia. Earth and Planetary Science Letters, 258(1-2): 269~282.
-
Snyder D B, Rondenay S, Bostock M G, Lockhart G D. 2004. Mapping the mantle lithosphere for diamond potential using teleseismic methods. Lithos, 77(1-4): 859~872.
-
Song Peng, Wang Tao, Tong Ying, Zhang Jianjun, Huang He, Qin Qie. 2019. Contrasting deep crustal compositions between the Altai and East Junggar orogens, SW Central Asian Orogenic Belt: evidence from zircon Hf isotopic mapping. Lithos, 328: 297~311.
-
Spencer C J, Cawood P A, Hawkesworth C J, Raub T D, Prave A R, Roberts N M W. 2014. Proterozoic onset of crustal reworking and collisional tectonics: reappraisal of the zircon oxygen isotope record. Geology(Boulder), 42(5): 451~454.
-
Spencer C J, Kirkland C L, Taylor R J M. 2016. Strategies towards statistically robust interpretations of in situ U-Pb zircon geochronology. Geoscience Frontiers, 7(4): 581~589.
-
Stern R J, Ali K A, Liégeois J P, Johnson P R, Kozdroj W, Kattan F H. 2010. Distribution and significance of pre-Neoproterozoic zircons in juvenile Neoproterozoic igneous rocks of the Arabian-Nubian Shield. American Journal of Science, 310(9): 791~811.
-
Sun Jian, Zhu Xiangkun, Belshaw N S, Chen Wei, Doroshkevich A G, Luo Wenjuan, Song Wenlei, Chen Beibei, Cheng Zhiguo, Li Zhihong, Wang Yang, Kynicky J, Henderson G M. 2021. Ca isotope systematics of carbonatites: insights into carbonatite source and evolution. Geochemical Perspectives Letters, 17(11-15): SI1~30.
-
Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geological Society, London, Special Publications, 42: 313~345.
-
Tang Gongjian, Wang Qiang, Wyman D A, Dan Wei. 2019. Crustal maturation through chemical weathering and crustal recycling revealed by Hf-O-B isotopes. Earth and Planetary Science Letters, 524: 115709.
-
Tang Gongjian, Wang Qiang, Wyman D A, Li Zhengxiang, Xu Yigang, Zhao Zhenhua. 2012. Recycling oceanic crust for continental crustal growth: Sr-Nd-Hf isotope evidence from granitoids in the western Junggar region, NW China. Lithos, 128-131: 73~83.
-
Tang Yanjie, Zhang Hongfu, Ying Jifeng, Su Benxun. 2013. Widespread refertilization of cratonic and circum-cratonic lithospheric mantle. Earth-Science Reviews, 118: 45~68.
-
Teng Jiwen, Zhang Zhongjie, Bai Wuming. 2004. Lithosphere Physics. Beijing: Science Press (in Chinese).
-
Teng Jiwen, Ruan Xiaomin, Zhang Yongqian, Hu Guoze, Yan Yafen. 2012. The stratificational velocity structure of crust and covering strata of upper mantle and the orbit of deep interaquifer substance locus of movement for Tibetan Plateau. Acta Petrologica Sinica, 28(12): 4077~4100 (in Chinese with English abstract).
-
Teng Jiwen, Yang Dinghui, Tian Xiaobo, Xu Tao, Chen Yun, Bai Zhiming, Liang Xiaofeng, Zhang Xi, Wu Jing, Liu Youshan. 2019. Geophysical investigation progresses of the Qinghai-Tibetan Plateau in the past 70 years (in Chinese). Scientia Sinica Terrae, 49(10): 1546~1564 (in Chinese).
-
Valley J W. 2003. Oxygen isotopes in zircon. Reviews in Mineralogy and Geochemistry, 53(1): 343~385.
-
Valley J W, Chiarenzelli J R, McLelland J M. 1994. Oxygen isotope geochemistry of zircon. Earth and Planetary Science Letters, 126(4): 187~206.
-
Valley J W, Kinny P D, Schulze D J, Spicuzza M J. 1998. Zircon megacrysts from kimberlite: oxygen isotope variability among mantle melts. Contributions to Mineralogy and Petrology, 133(1): 1~11.
-
Valley J W, Lackey J S, Cavosie A J, Clechenko C C, Spicuzza M J, Basei M A S, Bindeman I N, Ferreira V P, Sial A N, King E M, Peck W H, Sinha A K, Wei C S. 2005. 4. 4 billion years of crustal maturation: oxygen isotope ratios of magmatic zircon. Contributions to Mineralogy and Petrology, 150(6): 561~580.
-
Vervoort J D, Patchett P J, Albarède F, Blichert T J, Rudnick R, Downes H. 2000. Hf-Nd isotopic evolution of the lower crust. Earth and Planetary Science Letters, 181(1-2): 115~129.
-
Vervoort J D, Plank T, Prytulak J. 2011. The Hf-Nd isotopic composition of marine sediments. Geochimica et Cosmochimica Acta, 75(20): 5903~5926.
-
Wan Yusheng, Ma Mingzhu, Dong Chunyan, Xie Hangqiang, Xie Shiwen, Ren Peng, Liu Dunyi. 2015. Widespread late Neoarchean reworking of Meso-to Paleoarchean continental crust in the Anshan-Benxi area, North China Craton, as documented by U-Pb-Nd-Hf-O isotopes. American Journal of Science, 315(7): 620~670.
-
Wang Jian, Su Yuping, Zheng Jianping, Belousova E A, Griffin W L, Zhou Xiang, Dai Hongkun. 2020. Hidden Eoarchean crust in the southwestern Central Asian Orogenic Belt. Lithos, 360-361: 105437.
-
Wang Kai, Dong Shuwen, Yao Weihua, Zhang Yueqiao, Li Jianhua, Cui Jianjun, Han Baofu. 2020. Xenocrystic/inherited Precambrian zircons entrained within igneous rocks from eastern South China: tracking unexposed ancient crust and implications for Late Paleoproterozoic orogenesis. Gondwana Research, 84: 194~210.
-
Wang Peng, Mao Jinwen, Ye Huishou, Zhao Guochun, He Sheng, Yan Jianming, Wu Shaokui, Wan Limin, Zhao Xijiang. 2022. Zircon xenocryst Hf-O isotopic compositions in the Qiyugou Au orefield: a record of Paleoproterozoic oceanic slab subduction in the Trans-North China Orogen. Precambrian Research, 368: 106499.
-
Wang Tao, Jahn B M, Kovach P V, Tong Ying, Hong Dawei, Han Baofu. 2009. Nd-Sr isotopic mapping of the Chinese Altai and implications for continental growth in the Central Asian Orogenic Belt. Lithos, 110(1-4): 359~372.
-
Wang Tao, Hou Zengqian. 2018. Isotopic mapping and deep material probing(Ⅰ): revealing the compositional evolution of the lithosphere and crustal growth processes. Earth Science Frontiers, 25(6): 1~9 (in Chinese with English abstract).
-
Wang Tao, Huang He, Song Peng, Wu Huanhuan, Zhang Jianjun, Tong Ying, Du Kaiming, Qin Qie. 2020. Studies of crustal growth and deep lithospheric architecture and new issues: exemplified by the Central Asian Orogenic Belt(northern Xinjiang). Earth Science, 45(7): 2326~2344 (in Chinese with English abstract).
-
Wang Xuben, Yu Nian, Gao Song, Luo Wei, Cai Xuelin. 2017. Research progress in research on electrical structure of crust and upper mantle beneath the eastern margin of the Tibetan Plateau. Chinese Journal of Geophysics, 60(6): 2350~2370 (in Chinese with English abstract).
-
Wolf M, Romer R L, Glodny J. 2019. Isotope disequilibrium during partial melting of metasedimentary rocks. Geochimica et Cosmochimica Acta, 257: 163~183.
-
Wu Fuyuan, Li Xianhua, Zheng Yongfei, Gao Shan. 2007. Lu-Hf isotopic systematics and their application in petrology. Acta Petrologica Sinica, 23(2): 185~220 (in Chinese with English abstract).
-
Wu Gongjian. 1998. Key problems and orientation of the lithospheric research. Earth Science Frontiers, 5(1): 99~109 (in Chinese with English abstract).
-
Wu Yuanbao, Zheng Yongfei. 2004. Zircon genetic mineralogy and its constraints on U-Pb age interpretation. Chinese Science Bulletin, 49(16): 1589~1604 (in Chinese without English abstract).
-
Xiang Lu, Zheng Jianping, Zhai Mingguo. 2022. Archean to Paleoproterozoic crustal evolution of the southern Yangtze Block(South China): U-Pb age and Hf-isotope of zircon xenocrysts from the Paleozoic diamondiferous kimberlites. Precambrian Research, 374: 106651.
-
Xie Heping, Gao Feng, Ju Yang, Zhang Ru, Gao Mingzhong, Deng Jianhui. 2016. Novel idea and disruptive technologies for the exploration and research of deep earth. Advanced Engineer Ring Sciences, 49(1): 1~8 (in Chinese with English abstract).
-
Xu Bo, Kou Guanyu, Etschmann B, Liu Daiyue, Brugger J. 2020. Spectroscopic, Raman, EMPA, micro-XRF and micro-XANES analyses of sulphur concentration and oxidation state of natural apatite crystals. Crystals, 10(11): 1032.
-
Xu Bo, Hou Zengqian, Griffin W L, Zheng Yuanchuan, Wang Tao, Guo Zhen, Hou Jue, Santosh M, O'Reilly S Y. 2021a. Cenozoic lithospheric architecture and metallogenesis in southeastern Tibet. Earth-Science Reviews, 214: 103472.
-
Xu Bo, Hou Zengqian, Griffin W L, Lu Yongjun, Belousova E, Xu Jifeng, O'Reilly S Y. 2021b. Recycled volatiles determine fertility of porphyry deposits in collisional settings. American Mineralogist, 106(4): 656~661.
-
Xu Bo, Hou Zengqian, Griffin W L, Yu Jiaxing, Long Tao, Zhao Yi, Wang Tao, Fu Bin, Belousova E, O'Reilly S Y. 2022. Apatite halogens and Sr-O and zircon Hf-O isotopes: recycled volatiles in Jurassic porphyry ore systems in southern Tibet. Chemical Geology, 605: 120924.
-
Xu Shenglin, Ding Weicui, Chen Xuanhua, Li Tingdong, Han Lele, Liu Yong, Ma Feizhou, Wang Ye. 2022. Late Paleozoic crustal composition and growth in West Junggar: evidence from Sr-Nd-Pb isotopic mapping. Earth Science Frontiers, 29: 261~280 (in Chinese with English abstract).
-
Xu Yixian, Yang Bo, Zhang Anqi, Wu Shucheng, Zhu Lupei, Yang Yingjie, Wang Qingyan, Xia Qunke. 2020. Magnetotelluric imaging of a fossil oceanic plate in northwestern Xinjiang, China. Geology, 48(4): 385~389.
-
Xu Zhiqin, Yang Jingsui, Zhang Zeming, Liu Fulai, Yang Wencai, Jin Zhenmin, Wang Rucheng, Luo Liqiang, Huang Li, Dong Hailiang. 2005. Completion and achievement of the Chinese Continental Scientific Drilling(CCSD)Project. Geology in China, 32(2): 177~183 (in Chinese with English abstract).
-
Xu Zhiqin, Yang Jingsui, Li Haibing, Zhang Jianxin, Zeng Linsen, Jiang Mei. 2006. The Qinghai-Tibet plateau and continental dynamics: a review on terrain tectonic, collisional orogenesis, and processes and mechanisms for the rise of plateau. Geology in China, 33(2): 221~238 (in Chinese with English abstract).
-
Yang Gaoxue, Li Yongjun, Gu Pingyang, Yang Baokai, Tong Lili, Zhang Hongwei. 2012. Geochronological and geochemical study of the Darbut Ophiolitic Complex in the West Junggar(NW China): implications for petrogenesis and tectonic evolution. Gondwana Research, 21(4): 1037~1049.
-
Yang Liqiang, He Wenyan, Gao Xue, Li Nan, Qiu Kunfeng, Li Dapeng. 2022. Cratonic three-dimensional material structure tracing method for lithosphere. Acta Petrologica Sinica(For publication).
-
Yang Qidi, Wang Tao, Guo Lei, Tong Ying, Zhang Lei, Zhang Jianjun, Hou Zengqian. 2017. Nd isotopic variation of Paleozoic-Mesozoic granitoids from the Da Hinggan Mountains and adjacent areas, NE Asia: implications for the architecture and growth of continental crust. Lithos, 272: 164~184.
-
Yang Wencai, Yang Wuyang, Jin Zhenmin, Chen Zhenyan. 2004. Seismic fabric of lithosphere in sulu ultra-high pressure metamorphic zone. Scientia Sinica Terrae, 34(4): 307~319 (in Chinese).
-
Yang Wencai, Hou Zunze, Yu Changqing. 2015. Three dimensional density structural of the Tibetan Plateau and crustal mass movement. Chinese Journal of Geophysics, 58(11): 4223~4234 (in Chinese with English abstract).
-
Yin Jiyuan, Chen Wen, Xiao Wenjiao. 2020. Geochronological and geochemical study of mafic-intermediate dykes from the Northern West Junggar, NW China: source, petrogenesis and tectonic implications. Acta Geologica Sinica (English Edition), 94(S1): 78.
-
Yu Yang, Sun Min, Long Xiaoping, Li Pengfei, Zhao Guochun, Kroener A, Broussolle A, Yang Jinhui, Safonova I Y, Seltmann R, Sun Min, Xiao Wenjiao. 2017. Whole rock Nd-Hf isotopic study of I-type and peraluminous granitic rocks from the Chinese Altai: constraints on the nature of the lower crust and tectonic setting. Gondwana Research, 47: 131~141.
-
Zeng Lingsen, Asimow P D, Saleeby J B. 2005a. Coupling of anatectic reactions and dissolution of accessory phases and the Sr and Nd isotope systematics of anatectic melts from a metasedimentary source. Geochimica et Cosmochimica Acta, 69(14): 3671~3682.
-
Zeng Lingsen, Saleeby J B, Asimow P. 2005b. Nd isotope disequilibrium during crustal anatexis: a record from the Goat Ranch migmatite complex, southern Sierra Nevada Batholith, California. Geology(Boulder), 33(1): 53~56.
-
Zeng Lingsen, Chen Jing, Gao Lie, Chen Zhenyu. 2012. The geochemical nature of apatites in high Sr/Y two-mica granites from the North Himalayan Gneiss Domes, southern Tibet. Acta Petrologica Sinica, 28(9): 2981~2993 (in Chinese with English abstract).
-
Zhang Chen, Liu Dongdong, Zhang Xin, Spencer C, Tang Ming, Zeng Jianhui, Jiang Shu, Jolivet M, Kong Xiangye. 2020. Hafnium isotopic disequilibrium during sediment melting and assimilations. Geochemical Perspectives Letters, I2: 34~39.
-
Zhang Jianjun, Wang Tao, Zhang Lei, Tong Ying, Zhang Zhaochong, Shi Xingjun, Guo Lei, Huang He, Yang Qidi, Huang Wei, Zhao Jianxin, Ye Ke, Hou Jiyao. 2015. Tracking deep crust by zircon xenocrysts within igneous rocks from the northern Alxa, China: constraints on the southern boundary of the Central Asian Orogenic Belt. Journal of Asian Earth Sciences, 108: 150~169.
-
Zhang Jianjun, Wang Tao, Tong Ying, Zhang Zhaochong, Huang He, Zhang Lei, Guo Lei, Song Peng, Hou Zengqian. 2017. Tracking deep ancient crustal components by xenocrystic zircons of Paleozoic felsic igneous rocks from the Altai-East Junggar terrane and adjacent regions and its tectonic significance. Геодинамика и тектонофизика, 8(3): 605~607.
-
Zhang Shitao, Ma Qiang, Chen Huayong, Long Xiaoping, Chu Gaobin, Cheng Jiamin, Xia Xiaoping, Li Rucao. 2021. Precambrian crust growth and reworking of the eastern Yangtze Craton: insights from xenocrystic zircons in the lamprophyres from the Middle-Lower Yangtze Belt, China. Precambrian Research, 355: 106121.
-
Zhang Xiahui, Su Yuping, Zheng Jianping, Liu Penglei, Griffin W L, Ping Xianquan, Wang Jian, Zhou Liang. 2021. Metamorphic history and Neoarchean-Paleoproterozoic crustal growth of the central Trans-North China Orogen: evidence from granulite- to amphibolite-facies rocks of the Hengshan complex. Gondwana Research, 93: 162~183.
-
Zhang Yunying, Sun Min, Yin Jiyuan, Yuan Chao, Sun Zhen, Xia Xiaoping. 2022. Maturation of East Junggar oceanic arc related to supracrustal recycling driven by arc-arc collision: Perspectives from zircon Hf-O isotopes. International Journal of Earth Sciences(Geol Rundsch), https: //doi. org/10. 1007/s00531-022-02164-7.
-
Zhang Zeming, Dong Xin, Xiang Hua, He Zhenyu, Liou J G, Sun Weidong, Tatsumi Y, Teng Fangzhen, Yang Xiaoyong, Ling Mingxing. 2014. Metagabbros of the Gangdese Arc root, south Tibet: implications for the growth of continental crust. Geochimica et Cosmochimica Acta, 143: 268~284.
-
Zhao Wenjin, Wu Zhenhan, Shi Danian, Xiong Jiayu, Xue Guangqi, Su Heping, Hu Daogong, Ye Peisheng. 2008. Comprehensive deep profiling of Tibetan Plateau in the INDEPTH Project. Acta Geoscientical Sinica, 29(3): 328~342 (in Chinese with English abstract).
-
Zhao Zhenhua. 2010. Trace element geochemistry of accessory minerals and its applications in petrogenesis and metallogenesis. Earth Science Frontiers, 17(1): 267~286 (in Chinese with English abstract).
-
Zhao Zhidan, Liu Dong, Wang Qing, Zhu Dicheng, Dong Guochen, Zhou Su, Mo Xuanxue. 2018. Zircon trace elements and their use in probing deep processes. Earth Science Frontiers, 25(6): 124~135 (in Chinese with English abstract).
-
Zheng Jianping, Xia Bing, Ping Xianquan, Wei Ying, Tang Huayun, Su Yuping, Ma Qiang. 2021. Rock probes and seismic methods to constrain the structure, composition and evolution of the deep crust beneath North China Block. China Science Bulletin, 66(23): 3018~3031 (in Chinese with English abstract).
-
Zheng Yongfei, Li Shuguang, Chen Jiangfeng. 1998. Chemical geodynamics. Advance in Earth Sciences, 13(2): 121~128 (in Chinese with English abstract).
-
Zhou Xinmin, Sun Tao, Shen Weizhou, Shu Liangshu, Niu Yaoling. 2006. Petrogenesis of Mesozoic granitoids and volcanic rocks in south China: a response to tectonic evolution. Episodes, 29(1): 26~33.
-
Zhu Bingquan. 2001. Geochemical Province and Geochemical Abrupt Change Zone. Beijing: Science Press (in Chinese).
-
Zhu Bingquan, Chang Xiangyang. 1995. Geochemical boundary of Yangtze in South China and its control on the formation of super-large ore deposits. Scientia Sinica, 25(9): 1004~1008 (in Chinese).
-
Zhu Bingquan, Chang Xiangyang, Qiu Huaning, Sun Dazhong. 1998. Proterozoic basement characteristics in geochemical abrupt zones and their relation to the occurrence of super-large ore deposits. Scientia Sinica Terrae, S2: 63~70 (in Chinese).
-
Zhu Bingquan, Chang Xiangyang. 2001. Geochemical province and geochemical boundary. Advance in Earth Sciences, (2): 153~162 (in Chinese).
-
Zhu Rixiang, Chen Ling, Wu Fuyuan, Liu Junlai. 2011. Timing scale and mechanism of the destruction of the North China Craton. Science China Earth Sciences, 54: 789~797 (in Chinese).
-
Zhu Rixiang, Xu Yigang. 2019. The subduction of the west Pacific plate and the destruction of the North China Craton. Science China Earth Sciences, 62: 1340~1350 (in Chinese).
-
Zhu Xiaosan, Lu Minjie, Guan Ye, Zhang Yinghui. 2018. Experiment on revealing the regional metallogenic structure based on the aeromagnetic data in southern Da Hinggan Mountains. Geological Bulletin of China, 37(2-3): 299~313 (in Chinese with English abstract).
-
Zhu Xiaosan, Zheng Hongwei, Lu Minjie. 2019. Lateral variation of aeromagnetic anomaly in South China and its tectonic implications. International Journal of Earth Sciences, 108: 1493~1507.
-
Zhu Xiaosan, Lu Minjie. 2021. Fault structures and magmatic intrusions inferred from magnetic data for the Southern Great Xing'an Range, northern China. Ore Geology Reviews, 135: 104206.
-
Zhu Xiaosan, Wang Liangliang, Zhou Xiwen. 2022a. Structural features of the Jiangshao fault zone inferred from aeromagnetic data for South China and the East China Sea. Tectonophysics, 826: 229252.
-
Zhu Xiaosan, Lu Minjie, Zheng Hongwei. 2022b. Bouguer gravity anomaly in the andean orogenic belt and its dynamic implications for regional tectonic evolution. Acta Geologica Sinica, 96(3): 806~824.
-
Zhu Xiaosan, Wang Tao, Huang He, Zheng Hongwei. 2022c. Seep structures of northern Xinjiang and their tectonic implications: an aeromagnetic study. International Journal of Earth Sciences, in review.
-
Zou Xinyu, Jiang Jilian, Qin Kezhang, Zhang Yigang, Yang Wei, Li Xianhua. 2021. Progress in the principle and application of zircon trace element. Acta Petrologica Sinica, 37(4): 985~999 (in Chinese with English abstract).
-
池际尚. 1988. 中国东部新生代玄武岩及上地幔研究. 武汉: 中国地质大学出版社.
-
邓晋福. 1994. 岩石物理化学与岩石物理学. 地学前缘, (Z1): 57~63.
-
邓晋福, 苏尚国, 赵海玲, 莫宣学, 肖庆辉, 周肃, 刘翠, 赵国春. 2003. 华北地区燕山期岩石圈减薄的深部过程. 地学前缘, 10(3): 41~50.
-
董树文, 李廷栋, 陈宣华, 魏文博, 高锐, 吕庆田, 杨经绥, 王学求, 陈群策, 石耀霖, 黄大年. 2012. 我国深部探测技术与实验研究进展综述. 地球物理学报, 55(12): 3884~3901.
-
董树文, 李廷栋, 高锐, 吕庆田, 魏文博, 杨经绥, 王学求, 陈群策, 石耀霖, 黄大年, 陈宣华. 2013. 我国深部探测技术与实验研究与国际同步. 地球学报, 34(1): 7~23.
-
董树文, 李廷栋, 陈宣华, 高锐, 吕庆田, 石耀霖, 黄大年, 杨经绥, 王学求, 魏文博, 陈群策. 2014. 深部探测揭示中国地壳结构, 深部过程与成矿作用背景. 地学前缘, 21(3): 201~225.
-
高山, 张本仁, 金振民. 1999. 秦岭-大别造山带下地壳拆沉作用. 中国科学, 29(6): 532~541.
-
高山, 章军锋, 许文良, 刘勇胜. 2009. 拆沉作用与华北克拉通破坏. 科学通报, 54(14): 1962~1973.
-
侯增谦, 王涛. 2018. 同位素填图与深部物质探测(Ⅱ): 揭示地壳三维架构与区域成矿规律. 地学前缘, 25(6): 20~41.
-
侯增谦, 郑远川, 卢占武, 许博, 王长明, 张洪瑞. 2020. 青藏高原巨厚地壳生长、加厚与演化. 地质学报, 94(10): 2797~2815.
-
金振民, 欧新功, 徐海军, H Kern, 金淑燕, 侯光久. 2004. 中国大陆科学钻探主孔100~2000 m岩石弹性波速度: 对地震深反射的约束. 岩石学报, 20(1): 81~96.
-
刘光鼎. 2002. 回顾与展望——21世纪的固体地球物理. 地球物理学进展, 17(2): 191~197.
-
刘志超, 刘小驰, 俞良军, 王建刚. 2020. 喜马拉雅康巴淡色花岗岩的高分异成因及岩浆-热液演化特征. 南京大学学报(自然科学版), 56(6): 800~814.
-
吕庆田, 刘振东, 董树文, 严加永, 张永谦. 2015. “长江深断裂带”的构造性质: 深地震反射证据. 地球物理学报, 58(12): 4344~4359.
-
莫宣学. 2011. 岩浆与岩浆岩: 地球深部“探针”与演化记录. 自然杂志, 33(5): 255~269.
-
莫宣学. 2019. 岩浆作用与地球深部过程. 地球科学, 44(5): 1487~1493.
-
滕吉文, 张中杰, 白武明. 2004. 岩石圈物理学. 北京: 科学出版社.
-
滕吉文, 阮小敏, 张永谦, 胡国泽, 闫亚芬. 2012. 青藏高原地壳与上地幔成层速度结构与深部层间物质的运移轨迹. 岩石学报, 28(12): 4077~4100.
-
滕吉文, 杨顶辉, 田小波, 徐涛, 陈赟, 白志明, 梁晓峰, 张晰, 吴晶, 刘有山. 2019. 青藏高原深部地球物理探测70年. 中国科学: 地球科学, 49(10): 1546~1564.
-
王涛, 侯增谦. 2004. 同位素填图与深部物质探测(Ⅰ): 揭示岩石圈组成演变与地壳生长. 地学前缘, 25(6): 1~9.
-
王涛, 侯增谦. 2018. 同位素填图与深部物质探测(Ⅰ): 揭示岩石圈组成演变与地壳生长. 地学前缘, 25(6): 1~19.
-
王涛, 黄河, 宋鹏, 吴欢欢, 张建军, 童英, 杜开明, 秦切. 2020. 地壳生长及深部物质架构研究与问题: 以中亚造山带(北疆地区)为例. 地球科学, 45(7): 2326~2344.
-
王绪本, 余年, 高嵩, 罗威, 蔡学林. 2017. 青藏高原东缘地壳上地幔电性结构研究进展. 地球物理学报, 60(6): 2350~2370.
-
吴福元, 李献华, 郑永飞, 高山. 2007. Lu-Hf同位素体系及其岩石学应用. 岩石学报, 23(2): 185~220.
-
吴功建. 1998. 岩石圈研究的重要问题与研究方向. 地学前缘, 5(1): 99~109.
-
吴元保, 郑永飞. 2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报, 49(16): 1589~1604.
-
谢和平, 高峰, 鞠杨, 张茹, 高明忠, 邓建辉. 2016. 深地科学领域的若干颠覆性技术构想和研究方向. 工程科学与技术, 49(1): 1~8.
-
徐盛林, 丁伟翠, 陈宣华, 李廷栋, 韩乐乐, 刘勇, 马飞宙, 王叶. 2022. 西准噶尔晚古生代地壳组成与生长: 来自 Sr-Nd-Pb 同位素填图的证据. 地学前缘, 29: 261~280.
-
许志琴, 杨经绥, 张泽明, 刘福来, 杨文采, 金振民, 王汝成, 罗立强, 黄力, 董海良. 2005. 中国大陆科学钻探终孔及研究进展. 中国地质, 32(2): 177~183.
-
许志琴, 杨经绥, 李海兵, 张建新, 曾令森, 姜枚. 2006. 青藏高原与大陆动力学——地体拼合, 碰撞造山及高原隆升的深部驱动力. 中国地质, 33(2): 221~238.
-
杨立强, 和文言, 高雪, 王偲瑞, 李楠, 邱昆峰, 张良, 马强, 苏玉平, 李大鹏, 张智宇, 于红. 2022. 克拉通岩石圈三维物质架构示踪方法. 岩石学报, 待刊.
-
杨文采, 杨午阳, 金振民, 程振炎. 2004. 苏鲁超高压变质带岩石圈的地震组构. 中国科学(D 辑: 地球科学), 34(4): 307~319.
-
杨文采, 侯遵泽, 于常青. 2015. 青藏高原地壳的三维密度结构和物质运动. 地球物理学报, 58(11): 4223~4234.
-
曾令森, 陈晶, 高利娥, 陈振宇. 2012. 藏南北喜马拉雅穹窿高Sr/Y二云母花岗岩中磷灰石地球化学特征及其岩石学意义. 岩石学报, 28(9): 2981~2993.
-
张宏福, 杨岳衡. 2007. 华北克拉通东部含金刚石金伯利岩的侵位年龄和Sr-Nd-Hf同位素地球化学特征. 岩石学报, 23(2): 285~294.
-
赵文津, 吴珍汉, 史大年, 熊嘉育, 薛光琦, 宿和平, 胡道功, 叶培盛. 2008. 国际合作 INDEPTH项目横穿青藏高原的深部探测与综合研究. 地球学报, 29(3): 328~342.
-
赵振华. 2010. 副矿物微量元素地球化学特征在成岩成矿作用研究中的应用. 地学前缘, 17(1): 267~286.
-
赵志丹, 刘栋, 王青, 朱弟成, 董国臣, 周肃, 莫宣学. 2018. 锆石微量元素及其揭示的深部过程. 地学前缘, 25(6): 124~135.
-
郑建平, 夏冰, 平先权, 魏颖, 汤华云, 苏玉平, 马强. 2021. 岩石探针和地震探测手段约束华北深部地壳结构组成及演化. 科学通报, 66(23): 3018~3031.
-
郑永飞, 李曙光, 陈江峰. 1998. 化学地球动力学. 地球科学进展, 13(2): 121~128.
-
朱炳泉. 2001. 地球化学省与地球化学急变带. 北京: 科学出版社.
-
朱炳泉, 常向阳. 1995. 华南扬子地球化学边界及其对超大型矿床形成的控制. 中国科学B 辑, 25(9): 1004~1008.
-
朱炳泉, 常向阳, 邱华宁, 孙大中. 1998. 地球化学急变带的元古宙基底特征及其与超大型矿床产出的关系. 中国科学D辑, (S2): 63~70.
-
朱炳泉, 常向阳. 2001. 地球化学省与地球化学边界. 地球科学进展, (2): 153~162.
-
朱日祥, 陈凌, 吴福元, 刘俊来. 2011. 华北克拉通破坏的时间、范围与机制. 中国科学: 地球科学, 41: 583~592.
-
朱日祥, 徐义刚. 2019. 西太平洋板块俯冲与华北克拉通破坏. 中国科学: 地球科学, 49: 1346~1356.
-
朱小三, 卢民杰, 管烨, 张颖慧. 2018. 基于航磁资料揭示大兴安岭南部区域成矿构造. 地质通报, 37(2-3): 299~313.
-
邹心宇, 蒋济莲, 秦克章, 张毅刚, 杨蔚, 李献华. 2021. 锆石微量元素的理论基础及其应用研究进展. 岩石学报, 37(4): 985~999.
-
摘要
长期以来,岩石圈深部探测主要依赖地球物理手段和深部钻探,缺乏深部物质探测技术。对深部物质的了解也主要局限于两种途径:一类是地球物理推测方法,依据地表获得的深部岩石的物性测定,解释或推测深部物质的某些特征;另一类是捕虏体方法(Xenolith-based methodology),直接获取深部物质信息。本文重点探索的第三种途径,即充分利用地表出露的岩浆岩,通过岩石探针和同位素填图,示踪深部物源(物质)特征。特别是利用大数据分析和数字填图,了解深部物质三维架构及四维演变。在此基础上,总结上述三种途径,构建较完整的以岩石探针和同位素填图为核心的揭示岩石圈三维物质组成架构的方法体系。研究显示,在岩相学研究的基础上,通过系统开展Sr、Nd、Hf、Pb等多元同位素示踪填图,结合地球物理资料,可有效揭示岩石圈深部物质组成架构。通过中亚增生造山带(北疆)、青藏高原碰撞造山带(冈底斯-三江)和华北-扬子克拉通三个典型大地构造单元关键地区的实践,显示多元同位素示踪深部物质的一致性和有效性,以及同位素填图结果与地球物理探测结果的对应性。基于这些成果,笔者初步提出了揭示岩石圈三维物质组成架构的方法体系框架。该方法体系具有良好应用前景,有望成为与地球物理探测相结合和匹配的深部物质架构探测技术,为规范开展深部物质架构探测、物质演变过程及深部动力学过程研究提供技术支撑。
Abstract
To date, our understanding concerning materials in the deep lithosphere relies primarily on two ways: ① a geophysics-based methodology (in which geophysical features of deep components are speculated according to those of exposed rock); ② a xenolith-based methodology. In this study, we explore the third way, a methodology that employs rock probe and multi-proxy isotopic mapping as the core. Combined with synthetical petrographic and petrogenetic studies on magmatic rocks derived from multiple sources at different depths, rock probe and multi-proxy isotopic mapping can be a suite of valuable tools in deciphering the 3-D lithospheric material architecture and its 4-D evolutionary history. Based on this, we attempt to incorporate the three methodologies into a novel methodological system for exploring the lithospheric material architecture. This will shed new light on the detection of the deep lithosphere. Previous studies have revealed that, based on petrographic observations, isotopic mapping by integrating Sr, Nd, Hf, Pb and other isotopes, combined with geophysical results, is capable of unraveling the material architecture of the deep lithosphere. These advances demonstrate the consistency and effectiveness of multi-proxy isotopic systems in tracing deep materials and the correspondence between the isotopic mapping and geophysical detection results. Therefore, we propose a methodological framework for deciphering the 3-D material architecture of the lithosphere. The application of this methodology in the Central Asian Orogenic Belt (accretionary orogen), Tibet Plateau (collisional orogen), and the North China and Yangtze cratons has yielded some encouraging advances, indicating the significant application potential of this methodology. The achievements documented in this paper suggest that a novel methodological system has been prepared for deciphering the deep lithospheric materials and elucidating the architectures of crust and mantle.
Keywords
magmatic rock ； rock probe ； isotopic mapping ； deep material ； geophysics ； experimental simulation
