适用于壳幔混合物源成因花岗岩中暗色微粒包体的单斜辉石压力计

孙兰兰++杨家喜++周争艳++邹莉菲++马嫒嫒++关豪

摘要：花岗岩的壳幔混合物源成因及花岗岩中暗色微粒包体的成因信息是地质学关注的重要研究对象，但至今还没有适用于花岗岩中暗色微粒包体这种基性矿物组合的压力计。依据Patino Douce模拟壳幔岩浆混合过程的试验岩石学研究资料，发现具有基性矿物组合的暗色微粒包体中单斜辉石的Na、Al含量与压力密切相关，得出了压力与单斜辉石中的Na2O含量、Al2O3含量、硬玉分子摩尔分数、契尔马克辉石分子摩尔分数、硬玉分子与契尔马克辉石分子摩尔分数之和的线性关系，以及压力与斜方辉石中Al2O3含量的线性关系，进而得出适用于壳幔混合物源成因花岗岩中暗色微粒包体的有效地质压力计。这一压力计为了解壳幔混合物源成因花岗岩浆形成的深度提供了研究手段。将地质压力计应用于西秦岭天水地区车河花岗岩体及其中的暗色微粒包体研究，对了解壳幔混合作用形成的花岗岩等深部地质问题具有重要意义。

关键词：花岗岩；壳幔混合成因；暗色微粒包体；试验岩石学；单斜辉石；地质压力计；形成深度

中图分类号：P571；P589.1；P595文献标志码：A

0引言

自1969年Didier等根据花岗岩中的包体提出花岗岩的壳源型（Crustal）、幔源型（Mantle）和混合物源型（Mixed）的成因分类[12]以来，花岗岩的物源成为花岗岩成因问题的重要研究内容[39]。数十年来，由于地球化学、同位素技术、试验岩石学研究的快速发展以及大量的花岗岩实例研究，研究者们认识到：除了部分强过铝花岗岩是纯壳源的，拉斑质和碱性花岗岩是纯幔源的或基本是幔源的，其他分布于大陆造山带中的花岗岩都是混合物源成因的，在其形成过程中或多或少都有幔源物质的加入[10]。Pitcher研究表明，花岗岩的生成不仅是地壳物质再循环作用的结果，也有幔源组分的参与，较常见的是酸性岩浆和基性岩浆的混合[11]。1991年，Castro等提出了花岗岩的S、M和H型分类，其中H型即混合花岗岩（相当于壳幔型花岗岩）[12]。

Depaolo等通过对I型花岗岩和陆壳的研究认为，在壳幔相互作用过程中幔源岩浆提供了热量和物质[1318]。一些I型花岗岩中存在大量可以指示源岩信息的暗色微粒包体，并为揭示陆壳和岩石圈地幔的相互作用提供了重要证据。王涛总结并提出了混合花岗岩的类型，其中包括壳源岩浆和幔源岩浆的混合，如幔源玄武岩浆与酸性岩浆混合，充填于岩浆房的硅质壳源岩浆与渗透到岩浆房中的地幔熔体[19]。

Barbarin在对花岗岩的分类评述中，明确提出壳源（过铝质花岗岩）、幔源（拉斑质、碱性和过碱性花岗岩）及壳幔混合（钙碱性花岗岩）3种成因分类，并将它们与地球动力学环境联系起来[10]。自然界中壳幔混合成因的花岗岩是比较常见的。

Yang等通过锆石Hf同位素研究证明地幔参与了花岗质岩浆的形成[20]，试验岩石学也支持这一观点[21]。同时，一些学者认为地幔所提供的热量可能是花岗岩形成的重要因素。这也使得Petford等认为底侵作用和花岗岩成因有密切的联系[2224]。

玄武岩浆的底侵作用和幔源岩浆与壳源岩浆的混合作用是壳幔作用最直观的表现形式[2526]。根据大量的花岗质岩浆来源于壳幔相互作用这一认识，研究者们普遍认为，大量的软流圈或上地幔上涌的幔源岩浆以底侵方式囤积在地壳底部附近。这种幔源岩浆囤积作用带来的地幔热源引起下部地壳发生大规模变质作用和深熔作用，形成巨量花岗岩及未完全混合的镁铁质微粒包体[25，2728]。幔源岩浆的底侵作用引发了壳幔能量和物质的交换[2930]。

周珣若结合文献[3]对花岗岩中暗色微粒包体与岩浆混合作用的研究，认为包体在研究岩浆混合作用过程中有十分重要的意义[31]。王德滋等认为，镁铁质岩浆与长英质岩浆之间的混合作用是导致壳幔混源花岗岩类形成的主要机制[32]。暗色、细粒且具火成岩结构的岩石包体是指示岩浆混合作用存在的可靠证据。近些年来，对花岗岩中暗色微粒包体的研究也越来越受到重视[3335]。

Patino Douce的试验模拟了壳幔岩浆的混合过程[36]，混合岩浆早期结晶分异产生了基性矿物组分。通过对Patino Douce的试验资料分析发现，基性矿物组分中单斜辉石的Na、Al含量（质量分数，下同）与压力紧密相关，并由此推导出了壳幔混合成因花岗岩中暗色微粒包体（Mafic Microgranular Enclave，MME）的单斜辉石压力计。这一压力计将为壳幔混合物源成因花岗岩浆及其中暗色微粒包体的形成深度研究提供一项新工具。

1壳幔混合物源成因花岗岩中暗色微粒包体的单斜辉石压力计

单斜辉石是地壳和上地幔许多基性、超基性岩的重要矿物成分，它们的化学组成变化与其温度、压力和总的化学组成变化关系密切。通过研究这些变量如何控制辉石的化学组成，地质学家就可能由共存矿物的简单化学分析来推断这些岩石的形成条件[37]。从这些基性组分中单斜辉石的化学组成方面入手，可以推算出一种单斜辉石压力计，用以研究壳幔混合物源成因花岗岩及其中的暗色微粒包体。

对于花岗岩的壳幔混合物源成因问题，很早就有学者在试验岩石学方面进行过试验模拟研究。1995年，Patino Douce在温度保持在1 000 ℃的条件下，分别在压力（p）为5×105、7×105、10×105、12×105、15×105 Pa的条件下用等比例的玄武岩和黑云母片麻岩全部熔融混合，观察并研究了不同条件下矿物结晶情况和化学成分变化情况（表1）[36]。该试验模拟了壳源岩浆和幔源岩浆的混合过程，从试验结果基性矿物组合数据中发现，矿物组合单斜辉石中的Na、Al含量变化总是与压力变化关系密切。

根据Patino Douce的试验产物中基性矿物组合数据，以6个氧原子为基，分别计算出单斜辉石中各组分的含量（为了对比研究，也挑选了试验数据斜方辉石中的Al进行分析）和辉石各端元组分（表1）。他们是钙质单斜辉石Ca（Mg，Fe）Si2O6、镁铁辉石（Mg，Fe）2Si2O6、硬玉NaAlSi2O6、契尔马克辉石CaAl[AlSi]O6。根据计算结果，可以发现单斜辉石中的Na、Al含量与压力表现出较好的线性相关关系，并在此基础上分别得出了压力与单斜辉石中的Na2O含量、Al2O3含量、硬玉分子摩尔分数（XJd）、契尔马克辉石分子摩尔分数（XTsch）、硬玉分子与契尔马克辉石分子摩尔分数之和（XJd+XTsch）的线性关系，以及压力与斜方辉石中的Al2O3含量的线性关系，同时计算出这些线性关系的相关系数和相对误差（图1、表2）。

2讨论

2.1几种压力计的对比

从表2中辉石的压力计经验公式及相关参数可看出：压力与单斜辉石中Na、Al有关组分表现出良好的线性相关关系，具有较小的相对误差，其中尤其以Al2O3含量对应的压力计和契尔马克辉石对应的压力计线性关系最好。Putirka等研究认为通常

单斜辉石中的硬玉组分对压力是比较敏感的[3840]，结合前人研究成果及一些花岗岩实例研究，单斜辉石中与Na有关（Na2O含量、XJd）的压力计使用效果可能比与Al有关（Al2O3含量、XTsch）的压力计使用效果更好。对于与Al相关的压力计需要更进一步研究。斜方辉石压力计相比单斜辉石压力计线性关系略差，但在无单斜辉石存在的情况下，斜方辉石压力计也是有参考使用价值的。

需要注意的是，尽管得出了以单斜辉石为主的压力计，但其使用存在着限制条件，即该压力计是应用于暗色包体中含有一定量Na、Al的单斜辉石（Al离子数大于Na），而这些暗色包体是在壳幔岩浆混合过程中形成的，即该压力计在有混合岩浆结晶分异过程条件下才有效。另外，该试验是在特定条件下实现的，且试验数据不够多，对压力计的使用还是存在一定影响的。综合看来，该压力计还是有较大的应用价值。

2.2单斜辉石压力计的地质意义

底侵作用使壳幔边界聚集了玄武质岩浆，下地壳的增长是通过底侵的镁铁质火山岩物质实现的，壳幔边界的长英质岩浆在底侵作用过程中也起到了重要作用[41]。放射性稳定同位素的一些证据也表明，底侵作用不仅起到了将地幔热源带到下地壳的作用，而且在此过程中幔源物质和壳源物质发生了广泛的化学反应[4143]。

Patino Douce在试验中模拟幔源岩浆和壳源岩浆的混合，通过结晶分异产生了长英质熔体和镁铁质矿物组合，同位素和微量元素特征也反映出这些生成物具有幔源和壳源组分的特点[36]。在壳幔边界，通过底侵作用幔源岩浆和壳源岩浆发生热量传递和化学反应，混合岩浆经过早期结晶分异生成基性矿物组分，并跟随剩余岩浆一起上侵。花岗岩和其中的暗色微粒包体在这样的环境下产生，通过对这些基性微粒包体的研究可以反映出混合岩浆早期结晶分异时的环境条件。Patino Douce试验反映出基性矿物组分中的Na、Al含量与压力的密切相关关系，对于指示混合岩浆早期结晶分异的深度有着重要意义。

2.3应用实例

西秦岭天水地区车河花岗岩体寄主岩石为似斑状二长花岗岩，暗色微粒包体为二长闪长岩。暗色微粒包体内发育他形单斜辉石，这些单斜辉石颗粒被角闪石包裹。电子探针数据见表3中20组数据。

通过对比研究发现包体内的单颗粒单斜辉石内Na2O含量分布表现出明显差异，近似表现出环带的特征。以TSP3/3样品中11号点的电子探针数据[图2（b）、表3]为例，辉石颗粒中心部位（TSP3/3116）Na2O含量相对最高，达到1.21%，向边部呈现出降低的趋势。根据对应的压力计及深度计算（按0.275×105 Pa·km-1）大致估算出该单斜辉石的形成深度为46～76 km，形成深度部分已经接近壳幔边界。辉石中Al2O3含量、XJd、XTsch分布也表现出这种环带的特征。这种环带特征表明包体中的辉石是在一个压力逐渐降低的过程中形成的，即这些辉石是在岩浆不断上升侵位的过程中形成的。根据暗色微粒包体中单斜辉石的矿物学及地球化学特征和Patino Douce的试验可以推测早期分异结晶过程：高温基性岩浆与较低温的酸性岩浆混合，致使混合后的高温基性岩浆快速结晶分异，以单斜辉石为代表的偏基性矿物组分在此阶段开始结晶，而混合结晶后的剩余岩浆继续上侵进入上地壳。暗色微粒包体中的单斜辉石在压力逐渐降低的过程中形成，据此可推测出这些单斜辉石在开始结晶后随着基性岩浆继续上侵。

通过相关压力计的应用，车河花岗岩暗色微粒包体的相关计算结果（表4，按0.275×105 Pa·km-1计算）部分表现出了不一致性，可能有如下原因：

（1）推测与契尔马克辉石等有关的压力计受温度影响比较大，与Al有关的压力计计算出的压力和深度也会表现出较大误差。相比之下，与Na有关的压力计计算结果更可靠，对于与Al有关的压力计需要更进一步的研究。

（2）针对车河花岗岩体及其中的暗色包体这一实例而言，暗色包体中的单斜辉石颗粒周围多被角闪石包裹，导致单斜辉石边部Na含量偏高，这种包裹现象对数据及计算结果或多或少都有影响。

（3）本文提出了一种针对壳幔混合物源成因花岗岩中暗色微粒包体的研究思路。这个思路也存在着诸多问题，比如试验数据太少，试验条件大多固定受限，一些与压力计有关的平衡反应缺少相关试验的研究基础等。

3结语

（1）单斜辉石压力计是针对壳幔混合物源成因花岗岩中暗色微粒包体提出的，单斜辉石中含有适量的Na、Al（Al离子数大于Na）。这些压力计表现出良好的线性关系，在实际研究中应用是有意义的。该系列的地质压力计是在试验岩石学的基础上提出的，而自然条件下壳幔岩浆混合和结晶分异比较复

杂，所得的压力值难免存在一定误差，实际研究中要根据具体情况分析选择合适的对应压力计使用。

（2）将单斜辉石压力计与花岗岩实例研究结合，对揭示壳幔岩浆混合过程有重要意义。地幔与下地壳岩浆在壳幔边界混合，早期结晶分异的基性矿物组分开始结晶后随剩余岩浆上侵，形成了花岗岩中的暗色微粒包体。实例中单斜辉石组分环带分布特征反映了这一压力降低上侵的过程，据此可以约束壳幔岩浆混合后开始结晶分异的深度。

（3）单斜辉石压力计为了解壳幔混合成因花岗岩的形成深度等问题提供了新的研究手段和思路，但仍存在较多问题，如试验数据偏少，试验条件限制，一些平衡反应试验研究基础的缺乏等。对于缺乏相关试验基础的与契尔马克辉石有关的压力计有待进一步的研究。

参考文献：

References：

[1]DIDIER J，LAMEYRE J.Les Granites du Massif Central Francais：Etude Comparee des Leucogranites et Granodiorites[J].Contributions to Mineralogy and Petrology，1969，24（3）：219238.

[2]DIDIER J，DUTHOU J L，LAMEYRE J.Mantle and Crustal Granites：Genetic Classification of Orogenic Granites and the Nature of Their Enclaves[J].Journal of Volcanology and Geothermal Research，1982，14（1/2）：125132.

[3]刘昌实，朱金初，沈渭洲，等.华南花岗岩物源成因特征与陆壳演化[J].大地构造与成矿学，1990，14（2）：125138.

LIU Changshi，ZHU Jinchu，SHEN Weizhou，et al.Evolution of Continental Crust and Material Source （Genesis） of Granitoids in S.China[J].Geotectonica et Metallogenia，1990，14（2）：125138.

[4]BARBARIN B.Mafic Magmatic Enclaves and Mafic Rocks Associated with Some Granitoids of the Central Sierra Nevada Batholith，California：Nature，Origin，and Relations with the Hosts[J].Lithos，2005，80（1/2/3/4）：155177.

[5]YANG J H，WU F Y，CHUNG S L，et al.Multiple Sources for the Origin of Granites：Geochemical and Nd/Sr Isotopic Evidence from the Gudaoling Granite and Its Mafic Enclaves，Northeast China[J].Geochimica et Cosmochimica Acta，2004，68（21）：44694483.

[6]LAMEYRE J，BOWDEN P.Plutonic Rock Types Series：Discrimination of Various Granitoid Series and Related Rocks[J].Journal of Volcanology and Gethermal Research，1982，14（1/2）：169186.

[7]CHAPPELL B W，WHITE A J R.Two Contrasting Granite Types[J].Pacific Geology，1974，8：173174.

[8]MANIAR P D，PICCOLI P M.Tectonic Discrimination of Granitoids[J].Geological Society of America Bulletin，1989，101（5）：635643.

[9]ROSSI P，CHEVREMONT P.Classification des Associations Magmatiques Granitoides[J].Geochronique，1987，21：1418.

[10]BARBARIN B.A Review of the Relationships Between Granitoid Types，Their Origins and Their Geodynamic Environments[J].Lithos，1999，46（3）：605626.

[11]PITCHER W S.The Nature and Origin of Granite[M].Berlin：Springer，1997.

[12]CASTRO A，MORENOVENTAS I， DE LAROSA J D.Htype （Hybrid） Granitoids： A Proposed Revision of the Granitetype Classification and Nomenclature[J].Earthscience Reviews，1991，31（3/4）：237253.

[13] DEPAOLO D J.Neodymium Isotopes in the Colorado Front Range and Crustmantle Evolution in the Proterozoic[J].Nature，1981，291：193196.

[14]COLLINS W J.Evaluation of Petrogenetic Models for Lachlan Fold Belt Granitoids：Implications for Crustal Architecture and Tectonic Models[J].Australian Journal of Earth Science，1998，45（4）：483500.

[15]COLLINS W J，RICHARDS S R，HEALY B E，et al.Origin of Heterogeneous Mafic Enclaves by Twostage Hybridisation in Magma Conduits （Dykes） Below and in Granitic Magma Chambers[J].Geological Society of America Special Papers，2000，350：2745.

[16]SNYDER D，TAIT S.The Imprint of Basalt on the Geochemistry of Silicic Magmas[J].Earth and Planetary Science Letters，1998，160（3/4）：433445.

[17]CLYNNE M A.A Complex Magma Mixing Origin for Rock Erupted in 1915，Lassen Peak，California[J].Journal of Petrology，1999，40（1）：105132.

[18]QIN J F，LAI S C，GRAPES R，et al.Geochemical Evidence for Origin of Magma Mixing for the Triassic Monzonitic Granite and Its Enclaves at Mishuling in the Qinling Orogen （Central China）[J].Lithos，2009，112（3/4）：259276.

[19]王涛.花岗岩混合成因研究及大陆动力学意义[J].岩石学报，2000，16（2）：161168.

WANG Tao.Origin of Hybrid Granitoids and the Implications for Continental Dynamics[J].Acta Petrologica Sinica，2000，16（2）：161168.

[20]YANG J H，WU F Y，WILDE S A，et al.Tracing Magma Mixing in Granite Genesis：In Situ UPb Dating and Hfisotope Analysis of Zircons[J].Contributions to Mineralogy and Petrology，2007，153（2）：177190.

[21]PATINO DOUCE A E.What Do Experiments Tell Us About the Relative Contributions of Crust and Mantle to the Origin of Granitic Magmas？ [J].Geological Society，London，Special Publications，1999，168：5575.

[22]PETFORD N，GALLAGHER K.Partial Melting of Mafic （Amphibolitic） Lower Crust by Periodic Influx of Basaltic Magma[J].Earth and Planetary Science Letters，2001，193（3/4）：483499.

[23]赵振华，熊小林，王强，等.新疆北部晚古生代的底侵作用：来自橄榄玄粗岩与埃达克岩的证据[J].地质学报，2007，81（5）：606619.

ZHAO Zhenhua，XIONG Xiaolin，WANG Qiang，et al.Underplating During Late Paleozoic in North Xinjiang：Evidence from Shoshonitic Series Volcanic Rocks and Adakite[J].Acta Geologica Sinica，2007，81（5）：606619.

[24]杜杨松，秦新龙，李铉具.安徽铜陵地区中生代幔源岩浆底侵作用：来自矿物巨晶和岩石包体的证据[J].岩石矿物学杂志，2004，23（2）：109116.

DU Yangsong，QIN Xinlong，LI Xuanju.Mesozoic Mantlederived Magma Underplating in Tongling，Anhui Province：Evidence from Megacrysts and Xenoliths[J].Acta Petrologica et Mineralogica，2004，23（2）：109116.

[25]齐有强，胡瑞忠，刘燊，等.岩浆混合作用研究综述[J].矿物岩石地球化学通报，2008，27（4）：409416.

QI Youqiang，HU Ruizhong，LIU Shen，et al.Review on Magma Mixing and Mingling[J].Bulletin of Mineralogy，Petrology and Geochemistry，2008，27（4）：409416.

[26]王德滋.华南花岗岩研究的回顾与展望[J].高校地质学报，2004，10（3）：305314.

WANG Dezi.The Study of Granitic Rocks in South China：Looking Back and Forward[J].Geological Journal of China Universities，2004，10（3）：305314.

[27]肖庆辉，邱瑞照，邢作云，等.花岗岩成因研究前沿的认识[J].地质论评，2007，53（增）：1727.

XIAO Qinghui，QIU Ruizhao，XING Zuoyun，et al.Major Frontiers on Studies of Granite Formation[J].Geological Review，2007，53（S）：1727.

[28]李昌年.岩浆混合作用及其研究评述[J].地质科技情报，2002，21（4）：4954.

LI Changnian.Comment on the Magma Mixing and Their Research[J].Geological Science and Technology Information，2002，21（4）：4954.

[29]杜杨松.壳幔成矿学初探[J].矿床地质，1999，18（4）：341346.

DU Yangsong.A Preliminary Approach to Crustmantle Metallogeny[J].Mineral Deposits，1999，18（4）： 341346.

[30]ZHANG X H，MAO Q，ZHANG H F，et al.Mafic and Felsic Magma Interaction During the Construction of HighK CalcAlkaline Plutons Within a Metacratonic Passive Margin：The Early Permian Guyang Batholith from the Northern North China Craton[J].Lithos，2011，125（1/2）：569591.

[31]周珣若.花岗岩混合作用[J].地学前缘，1994，1（1/2）：8797.

ZHOU Xunruo.Hybridization in the Genesis of Granitoids[J].Earth Science Frontiers，1994，1（1/2）：8797.

[32]王德滋，谢磊.岩浆混合作用：来自岩石包体的证据[J].高校地质学报，2008，14（1）：1621.

WANG Dezi，XIE Lei.Magma Mingling：Evidence from Enclaves[J].Geological Journal of China Universities，2008，14（1）：1621.

[33]薛玉山，柳振江，王建平，等.花岗岩暗色微粒包体特征及其研究方向[J].地质找矿论丛，2014，29（1）：2030.

XUE Yushan，LIU Zhenjiang，WANG Jianping，et al.Characteristics of Mafic Microgranule Enclaves（MMEs） in Granitoids and Future Study Prospects[J].Contributions to Geology and Mineral Resources Research，2014，29（1）：2030.

[34]陈广俊，孙丰月，李碧乐，等.东昆仑沟里地区暗色包体及其寄主岩石地球化学特征及成因[J].吉林大学学报：地球科学版，2014，44（3）：892904.

CHEN Guangjun，SUN Fengyue，LI Bile，et al.Geochemistry and Petrogenesis of Gouli Mafic Enclaves and Their Host Rocks in Eastern Kunlun[J].Journal of Jilin University：Earth Science Edition，2014，44（3）：892904.

[35]高永宝，李侃，钱兵，等.东昆仑卡而却卡铜矿区花岗闪长岩及其暗色微粒包体成因：锆石UPb年龄、岩石地球化学及SrNdHf同位素证据[J].中国地质，2015，42（3）：646662.

GAO Yongbao，LI Kan，QIAN Bing，et al.The Genesis of Granodiorites and Dark Enclaves from the Kaerqueka Deposit in East Kunlun Belt：Evidence from Zircon UPb Dating，Geochemistry and SrNdHf Isotopic Compositions[J].Geology in China，2015，42（3）：646662.

[36]PATINO DOUCE A E.Experimental Generation of Hybrid Silicic Melts by Reaction of HighA1 Basalt with Metamorphic Rocks[J].Journal of Geophysical Research，1995，100（B8）：1562315639.

[37]HERZBERG C T.Pyroxene Geothermometry and Geobarometry：Experimental and Thermodynamic Evaluation of Some Subsolidus Phase Relations Involving Pyroxenes in the System CaOMgOAl2O3SiO2[J].Geochimica et Cosmochimica Acta，1978，42（7）：945957.

[38]PUTIRKA K，JOHNSON M，Kinzler R，et al.Thermobarometry of Mafic Igneous Rocks Based on Clinopyroxeneliquid Equilibria，030 kbar[J].Contributions to Mineralogy and Petrology，1996，123（1）：92108.

[39]PUTIRKA K D，MIKAELIAN H，RYERSON F，et al.New Clinopyroxeneliquid Thermobarometers for Mafic，Evolved and Volatilebearing Lava Compositions，with Applications to Lavas from Tibet and the Snake River Plain，Idaho[J].American Mineralogist，2003，88（10）：15421554.

[40]雷敏，吴才来，高前明，等.铜陵地区中酸性侵入岩及其包体的成因和矿物温压计的应用[J].岩石矿物学杂志，2010，29（3）：271288.

LEI Min，WU Cailai，GAO Qianming，et al.Petrogenesis of Intermediateacid Intrusive Rocks and Enclaves in Tongling Area and the Application of Mineral Thermobarometry[J].Acta Petrologica et Mineralogica，2010，29（3）：271288.

[41]KEMPTON P D，HARMON R S.Oxygen Isotope Evidence for Largescale Hybridization of the Lower Crust During Magmatic Underplating[J].Geochimica et Cosmochimica Acta，1992，56（3）：971986.

[42]KEMPTON P D，HARMON R S，HAWKESWORTH C J，et al.Petrology and Geochemistry of Lower Crustal Granulites from the Geronimo Volcanic Field，Southeastern Arizona[J].Geochimica et Cosmochimica Acta，1990，54（12）：34013426.

[43]DOWNES H，KEMPTON P D，BRIOT D，et al.Pb and O Isotope Systematics in Granulite Facies Xenoliths，French Massif Central：Implications for Crustal Processes[J].Earth and Planetary Science Letters，1991，102（3/4）：342357.