Chronology and geochemical characteristics of quartz diorite in Gegen' aobao, Dong Ujimqin Banner, Inner Mongolia, and its geological significance
-
摘要:
内蒙古东乌旗地区新近发现的葛根敖包铅锌矿,成矿与赋矿岩体石英闪长岩关系密切,对石英闪长岩进行系统的LAICP-MS锆石U-Pb定年及锆石Hf同位素分析,并对石英闪长岩进行主量、微量和稀土元素分析。锆石U-Pb定年结果显示,石英闪长岩2个样品的锆石U-Pb年龄分别为299±1Ma和301±1Ma,均为晚石炭世。岩石主量、微量及稀土元素分析显示,石英闪长岩具有典型的富硅、富碱特征,Na2O/K2O值为1.14~1.93,平均1.39,明显富钠高钾,属高钾钙碱性系列;铝饱和指数A/CNK值为0.79~0.95,为准铝质或铝不饱和花岗岩类;轻、重稀土元素分馏显著,球粒陨石标准化稀土元素配分模式表现为明显的右倾,显示岩浆源区可能有石榴子石残留;并有微弱负Eu异常,δEu=0.80~0.91,δCe=0.89~0.96,均靠近1,指示源区无斜长石残留;同时具低Sr、高Yb特征;并指示高场强元素U、La、Nd、Zr相对富集,而Th、Nb、Sr、P、Ti相对亏损,且大离子亲石元素Rb、Ba、K相对富集,Nb/Ta值低,指示岩浆源区可能为幔源为主。锆石Hf同位素分析显示,εHf(t)值全部为正值,+10.8~+13.6(全部大于10,较高),TDM2介于440~629Ma之间,相对集中,明显大于其锆石U-Pb年龄,其年轻的TDM2年龄和较高的εHf(t)值表明,岩浆中有相当大的幔源物质的贡献。结合区域构造演化,认为葛根敖包石英闪长岩岩浆源区可能为亏损地幔物质部分熔融并在上升过程中受到少量壳源物质的混染,为活动陆缘弧环境的产物。
-
关键词:
- 内蒙古东乌旗 /
- 葛根敖包 /
- LA-ICP-MS锆石U-Pb年龄 /
- 地球化学 /
- Hf同位素
Abstract:The newly found Gegen'aobao lead-zinc deposit is within the quartz diorite. Because the exploration is insufficient, the metal resources are unclear. As a newly found deposit, there is no previous research data, and hence systematic isotope dating must be done to determine the age of metallogenesis for the deposit. The U-Pb dating shows that the quartz diorite has the ages of 299±1~301±1Ma. Major element geochemistry shows that the quartz diorite is characterized by high SiO2 and K2O, whereas trace element geochemistry shows that the quartz diorite is a "right-inclined" shape in the chondrite-normalized REE patterns. The rocks are rich in large ion lithophile elements (K, Rb, Ba) and LREEs and depleted in high field strength elements (Th, Nb, Sr, P, Ti) and HREE, with the loss of high field strength elements (HFSEs) such as Th, Nb, P and Ti. The quartz diorite is of high-K calc-alkaline series and has a little negative Eu anomaly (Eu/Eu*=0.80~0.91), with low A/CNK (0.79~0.95). These characteristics shows that the Gegenaobao quartz diorite is mantle source granite. The in situ zircon Hf isotopic data are in agreement with U-Pb age data of zircon grains, and the 176Hf/177Hf isotopic ratio of the quartz diorite is 0.282906~0.282984, with εHf(t) values >0, indicating that the source region of the quartz diorite magma was depleted mantle or newly created lower continental crust beneath the study area. Based on geochemistry, isotopes and regional tectonic evolution history, the authors hold that the quartz diorite was derived from a depleted mantle, which was probably related to an active continental margin arc environment.
-
Keywords:
- Gegen'aobao, Inner Mongolia /
- LA-ICP-MS zircon U-Pb age /
- geochemistry /
- Hf isotope
-
随着城市化进程的加快与人口数量迅速增多,生态环境变化对农业生产的影响已成为全球范围重点关注的问题之一[1]。据《全国生态环境十年变化(2000—2010年)调查评估报告》调查结果表明,全国生态安全形势极为严峻,生态保护与发展矛盾突出,主要表现在农业生产与开发导致的水土流失、土地沙化等问题加剧。中国作为以农业生产为国民经济基础的国家,农业生产的合理布局及可持续发展一直受到政府部门高度重视。在《中共中央、国务院关于建立国土空间规划体系并监督实施的若干意见》中指出“要在资源环境承载能力和国土空间开发适宜性评价的基础上,科学有序统筹布局生态、农业、城镇等功能空间”。
土地生产及农业适宜评估是农业区域开发与规划的重要步骤,这一领域国外研究多与GIS技术有机结合且多注重模型构建及单一作物的GIS应用[2-3]。Stork等[4]基于GIS方法在区域尺度上利用供水条件及土壤质量两大因素来评价薪柴生产潜力区域。Kumar等[5]采用频率比(FR)、层次分析法(AHP),结合遥感(RS)和地理信息系统(GIS)多重手段共同探索潜在的农业用地,对制定研究区可持续农业管理计划具有极好作用。Javad等[6]基于遥感(RS)和GIS技术,采用层次分析法和物质元素技术来评估和绘制大麦的种植布局,实验结果表明层次分析法和物质元素技术混合可使模型准确度更高。而国内研究偏重于应用,早期多以与农业生产相关的土壤、气候、地形等单一因素或单一植被展开研究。如土壤重金属元素的来源辨识、对生态环境污染风险的评价研究[7-9],或土壤养分含量动态变化跟踪及土培对策研究[10-12]。土地利用方式引起的土壤理化性质及生态效应变化等方面也有研究[13-15]。此外,在降水量变化特征及自然灾害对农业生产影响方面的研究也较为广泛[15-18]。现阶段随着国土空间规划相关政策提出,“资源环境承载能力与国土空间开发适宜性”(简称“双评价”)研究也逐渐开展[19-22]。总体看,已有研究大多具有较强的针对性和单一性,而关于不同时间尺度下的区域农业生产环境变化开展的农业适宜性评价鲜有报道。
本文运用ArcGIS的空间分析方法,采用南通市近年来气候条件变化、土壤肥力与土壤重金属污染情况等评价因子加权求和的方式,对南通市县市地区的农业生产适宜性进行分析,为1:50万精度及大比例尺图件关于农业生产要素的区域分析及动态变化影响分析提供参考方法,为推进南通市及更小尺度范围的地区了解区域农业生产环境变化并制定相应环境改善政策、开展农作物产业布局研究、识别农业生产合理区、国土空间规划方案制定提供依据支持。
1. 研究区概况
南通市地处江苏省东南部,处于沿海经济带与长江经济带T字形交汇中心,富集“黄金海岸”与“黄金水道”两大优越的空间地理位置资源,是长江流域重要的港口城市之一,也是中国沿海地区土地资源最丰富的地区之一。南通市凭借其独特的地理位置与丰富的自然资源已成为国内沿海地区农业优质发展的典型地区。南通市属长江三角洲冲积平原,全境地势平坦,多数地区海拔在2~6 m之间。全市现有耕地44.87×104 m2,以水田和旱地为主,当地农业产品以粮棉为主,盛产棉花和特种经济作物,耕地产出率高于全省平均水平[9]。南通市土壤多为水稻土、潮土及滨海盐土,整体土壤肥力较高,呈碱性土质[10]。
已有调查表明,当地年平均气温达到15~17℃,年平均日照时数在1700~2200 h之间,年平均降水量为1000~1100 mm,属于典型的气候温和、光照充足、雨水丰沛的农业适宜区,地处海陆过渡带(滨江临海),受季风影响明显,也是典型的自然灾害频发区。下辖海安县、如东县、如皋市、启东市、海门市,崇川区、港闸区、通州区,考虑到获取数据的重合性,本次研究将崇川区与港闸区合并为南通市辖区。
2. 数据来源与评价体系
2.1 数据来源
本次研究所用数据主要为空间数据和统计数据,其中,空间数据包含南通市行政区划图、南通市基础地理图;统计数据包括积温、年降水量、日照时数、土壤养分、土壤重金属元素含量。积温采用当地一年内日平均气温≥10℃的持续期间日平均气温的总和作为年积温数据;年降水量、日照时数来源于《南通市统计年鉴》(2013—2017);土壤肥力数据、土壤重金属数据来源于土壤肥力普查数据及巩万合等[8, 10]关于南通市土壤肥力及重金属的相关研究。具体如表 1所示。
表 1 数据来源Table 1. Data sources数据名称 数据来源 数据类型 积温 天气网历史天气数据的日平均温度 统计数据 年降水量 《南通市统计年鉴》(2013—2017) 统计数据 日照时数 《南通市统计年鉴》(2013—2017) 统计数据 土壤养分数据 全国第二次土壤普查数据 统计数据 土壤重金属数据 全国第二次土壤普查数据、巩万合等[8, 10]关于南通市土壤肥力及重金属 统计数据 南通市行政区划图 地理信息共享平台(https://www.tianditu.gov.cn) 空间数据 南通市基础地理图 地理信息共享平台(https://www.tianditu.gov.cn) 空间数据 2.2 评价体系构建
农业生产过程极易受多种因素影响,因此结合南通市生态环境特性与农业生产模式,选取积温、年降水量、日照时数、土壤养分(有机质、P、K)、土壤重金属(Cd、Cr、Hg、Pb、As)11个评价指标因子构建农业生产适应性评价体系,如表 2所示。这些评价因子对农业生产影响的范围及强度存在差异,结合南通市实际情况,采用德尔菲法[23-25]将所有指标对当地农业影响的重要程度逐一评分,最终确定各指标权重。在此基础上参考“双评价指南”,并与专家评定结合,共同划定同一因子的适宜程度等级。
表 2 农业生产适应性评价体系Table 2. Evaluation system of agricultural production adaptability评价指标 指标因子 权重 分值 分类 评价指标 指标因子 权重 分值 分类 气候 积温/℃ 0.15 1 < 5000 土壤 土壤磷含量/(mg·kg-1) 0.1 2 50~70 2 5000~5500 3 70~90 3 5500~5700 4 ≥90 4 ≥5700 Cd/(mg·kg-1) 0.05 1 ≥0.2 年降水量/mm 0.1 1 < 1000 2 0.1~0.2 2 1000~1500 3 < 0.1 3 ≥1500 Cr/(mg·kg-1) 0.05 1 ≥60 日照时数/h 0.1 1 1000~1400 2 50~60 2 1400~2200 3 < 50 3 2200~3000 Hg/(mg·kg-1) 0.05 1 ≥0.12 土壤 土壤有机质/% 0.1 1 < 0.6 2 0.1~0.12 2 0.6~1 3 < 0.1 3 1~2 Pb/(mg·kg-1) 0.05 1 ≥25 4 2~3 2 20~25 土壤钾含量/(mg·kg-1) 0.1 1 < 10 3 15~20 2 10~20 4 < 15 3 20~30 As/(mg·kg-1) 0.05 1 ≥8 4 ≥30 2 6~8 土壤钾含量/(mg·kg-1) 0.1 1 < 50 3 < 6 评价单元的综合得分是确定土地开发建设适宜性等级的基础,本次研究通过ArcGIS的空间统计方法,采用11个适宜性评价因子加权求和的数据模型来构建农业生产适宜性指数,计算公式如下所示:
Ti=n∑i=1Wi∗Ci (1) 式中:T代表南通市农业生产适宜性指数及综合得分,i=1,…,n代表第i个评价因子;W表示第i个评价因子的权重;C表示对应的评价因子等级评分。
基于ArcGIS软件平台,以南通市行政区划图为本底,所属县级区域为基本单位,对采集获取的各项数据进行空间叠加分析与属性空间化,构成多因子属性值的综合评价图,进而对南通市农业适宜性程度进行综合评价,划分为农业最适宜区、农业适宜区、农业较适宜区、农业一般适宜区[26-27]。
3. 评价过程与结果
3.1 气候变化
了解区域积温、降水演变规律、区域日照时数,可为农业环境条件筛选适宜种植的作物类型提供参考。
3.1.1 积温年际变化特征
分析积温年际变化有利于对该地区未来的农业生产气候条件做出合理预判。本研究对南通市近5年(2014—2018年)的积温年际变化进行统计,如图 1所示,近5年南通市不同区域积温主要分布在5358.5~5826.5℃之间,年际变化波动幅度较大,随年份递增积温整体呈波动上升趋势,这与王涛等[16]的研究结果一致。对比同一年份南通市各区域之间的积温变化,发现不同地区的积温变化特征基本随年份呈现统一变化规律。经统计,近5年南通市不同区域积温均高于5300℃,产生的热量满足了大部分农作物的需求,属于农业生产的适宜区域。
3.1.2 区域积温空间分布特征
区域积温图可通过直观展示积温具体情况来挖掘其空间分布特征。利用ArcGIS软件平台,以县级区域为单位叠加空间属性将南通市积温分布划分为4类积温区域(综合各地近5年平均积温)。如图 2所示,积温较高区域的温度集中在5563~5712℃之间,主要分布在南通市辖区、通州、海门、启东4个地区,占据了全市47.18%的面积,这部分地区纬度偏低、人口稠密、城市建设用地布局紧张、临江沿海工业聚集程度拥挤,长年处于南通市高积温区域。如东市、海安市区域积温在5494~5563℃之间,处于积温梯度中部,如皋地区积温低于5494℃,处于南通市低积温区域。
3.1.3 降水量与日照时数时间分布特征
图 3为南通市年均降水量与年日照时数的年际变化情况(2013—2017年),从图中历年降水数据可知,年降水量在2013—2016年基本呈逐年递增趋势,在2016年达到近年降水量最高值1921.4 mm,到2017年降水量陡减,同比上年减少了39.8%,表明南通市年降水量起伏波动变化大,极易受洪涝、干旱等气象灾害影响,因此应重点预测降水受灾风险区,提前做好防护措施,避免农业生产与城市基础设施遭受破坏。年日照时数的变化特征与年降水量完全相反,两者之间的变化规律呈负相关。
降水量与日照时数的月分布对农作物生长影响、农业环境变化、灾害风险预测提供分析依据。图 4为南通市降水量月分布情况(2013—2017年月平均数据),通过对降水量的分析,发现南通市全年降水大多集中于5~11月之间,占据年降水量的76.2%~86.3%,其中南通地区的降水峰值主要集中在6月中下旬至7月上中旬的梅雨期,因此应在梅雨期前夕做好重点防御工作。月日照时数按时间分布来看,6~7月的光照时数明显偏低,在光照强度与时间不足的条件下农作物难以进行有效的光合作用,进而导致农作物的减产与停产。此外,梅雨期降水集中,空气湿度与土壤湿度大,极易诱发病虫害[21]。
3.1.5 气候变化对南通市农业生产影响分析
农业对气候变化的反应敏感而脆弱,气候变化还会导致农作物生长过程发生改变[28]。南通市现代农业以小麦、大麦、稻谷等粮食作物与大豆、花生、油菜、棉花等经济作物作为主要农作物,农作物产业分布具备如下的地域特点,小麦、大麦、稻谷等粮食作物多分布在海安、如东、如皋等地区,花生、油菜、棉花等经济作物主要集中于启东、海门等沿海地区。南通市主要农作物的产业空间分布在一定程度上符合农作物对气候条件的要求。
根据南通市积温空间分布、降水量与日照特征对照分析,发现高积温气候适宜区域与经济作物分布区域高度重合,高积温区受沿海区位影响有降水量大、日照时数偏小等特点,可见对光热要求较高的油菜、棉花等作物极易受气候变化影响。对于处于粮食主产区、积温梯度中部的如东与海安来讲,气候变化特征与强度可能会对喜长日照的小麦产生影响,相对于喜高温、多湿、短日照的稻谷,气候变化对其影响较小。
总体看,气候变化对不同农作物的影响强度与范围差异性较大,且难以把控气候演变规律,因此可通过调整农田种植结构与农业布局模式、选育良种、加强气候灾害防控等措施来适应气候变化[29]。
3.2 土壤养分
3.2.1 土壤养分年际变化特征
图 5为南通市养分年际变化情况,随年份增加,全氮与有效磷呈递增趋势,土壤全氮含量浮动趋于稳定,至2012年达到12.5 mg/kg,增幅不足27%。有效磷的起伏变化较大,1982—2012年土壤有效磷增幅超过390%,至2012年含量达到20.6 mg/kg。土壤全氮与有效磷的逐年递增主要来自于有机肥的施用,与之相比,土壤速效钾的年际浮动极特殊,自1982—2012年经历了先减少后增加再减少的连续起伏过程。根据全国第二次土壤普查土壤养分分级标准,南通市土壤全氮与有效磷含量均满足土壤养分标准的二级要求,速效钾含量仅达到四级标准要求。
3.2.2 土壤养分评价
利用南通市各区域土壤养分空间分析,将土壤养分划为养分丰富、养分充足、养分较充足及养分基本充足4级(图 6)。海安、如东因土壤养分条件优越,列为土壤养分丰富区,通州、启东养分含量相对有所下滑,划为土壤养分充足地区,市辖区与海门市划为养分较充足区,与其他区域相比,如皋耕地面积相对偏少且土壤养分含量略有不足,可划分为养分基本充足地区。
从农业角度分析,海安与如东是养分丰富的粮食主产区,满足小麦、水稻等作物对土壤养分的需求,可进一步通过调整农业布局增加土壤养分丰富区域的利用率。通州、启东地区土壤养分充足,可达到油菜、蔬菜、棉花等经济作物的养分要求。
3.3 土壤重金属污染
3.3.1 土壤重金属污染统计分析
土壤重金属是评价耕地土壤环境质量的一项重要指标。表 3为南通市土壤重金属Hg、Pb、Cd、As和Cr的元素含量,根据已有研究发现[7-9],江苏省海岸带主要为碱性土壤,pH多集中在7.5~8.5之间,因此参考《土壤环境质量标准》(GB15618—1995)二级标准中土壤pH>7.5的标准值。
依据土壤环境质量标准,南通市土壤重金属含量均符合标准要求。通过对比,南通市各重金属平均含量与江苏省海岸带土壤重金属含量相似度较高,其中Cd元素含量均超过了江苏省土壤背景值,表明南通市及海岸带重金属污染来源可能相似,南通市是重要的粮食基地,有机肥、农药等大量施用增加了土壤养分,也造成了农用地土壤Cd等其他重金属富集。此外,南通市土壤中Hg元素的平均含量较高,这部分Hg的污染源主要来自工业三废与煤炭燃烧,土壤中的Hg元素富集会导致土壤肥力降低,妨碍农作物根系生长,通过食物链进入人体甚至危害人类健康。近年来,南通市第二产业发展迅速,已占据南通GDP比重的48.4%,可见工业发展对农业生产环境的影响极明显。其他重金属元素含量均低于江苏省土壤背景值。
3.3.2 土壤重金属污染评价
图 7为南通市各地区重金属含量,表 4为各地区综合污染指数,根据已有研究对土壤重金属的内梅罗指数进行计算与评价[7]。整体看,南通市各地区的综合污染指数均在0.26~0.4之间,污染程度处于安全阶段。从综合指数的分布看,重金属的分布也存在一定规律,启东的重金属含量最高,通州、海门、南通市辖区3个地区的指数相近,处于0.3~0.32之间,如皋、海安、如东的综合指数较低,可见重金属污染综合指数基本按照南东—北西方向由大到小呈阶梯状分布。沿海地带是工业发展、进出口集散地、人口密集的重要位置,因此由人类活动引发的重金属富集程度较内陆地区有所抬升。总体看,南通市各地区重金属污染对农业生产影响较小。
表 4 南通市各地区单因子指数与内梅罗综合污染指数Table 4. Index of single factor and comprehensive pollution index in different areas of Nantong区域 单因子指数 内梅罗综合指数 Cd Cr Pb As Hg 如皋 0.10 0.13 0.05 0.29 0.09 0.26 海安 0.11 0.13 0.04 0.28 0.09 0.26 启东 0.18 0.17 0.08 0.51 0.09 0.40 如东 0.18 0.15 0.05 0.30 0.09 0.27 通州 0.10 0.14 0.06 0.39 0.14 0.32 海门 0.17 0.16 0.05 0.34 0.10 0.30 南通 0.15 0.14 0.06 0.37 0.12 0.31 3.3.3 土壤重金属相关性分析
多种重金属元素之间的相关性分析基本反映了元素之间的关联情况,有助于辨识重金属污染来源[31]。表 5为土壤重金属之间的相关系数。Cd-Cr、Cr-As之间的Pearson相关性系数为0.792、0.758,p < 0.05,正相关呈显著,Pb-As之间的相关性系数为0.966,通过p < 0.01的检验,说明Pb与As之间的元素相关性很大,元素之间的同源性可能较高。其他元素之间的相关性略低。
表 5 土壤重金属之间的相关系数Table 5. Correlation coefficient among different heavy metals in soil项目 Cd Cr Pb As Hg Cd 1 Cr 0.792* 1 Pb 0.41 0.673 1 As 0.412 0.758* 0.966** 1 Hg -0.361 -0.163 0.214 0.143 1 注:*代表显著性相关(p < 0.05),**代表极显著相关(p < 0.01) 3.4 江苏海岸带滩涂资源与农业利用
土壤盐渍化是海岸带普遍存在的现象之一,海岸带大部分盐渍土主要集中在沿海滩涂,然而滩涂土壤的含盐量在空间分布上的不均与季节变化上的巨大差异导致盐分数据难以利用。海岸滩涂湿地实际是一种具有较强的气候调节、净化环境等功能的土地后备资源,随着土地资源的日渐紧张,大部分滩涂资源已被开发利用,且围垦后超过60%的土地作为耕地资源进行利用[32]。据了解,现今江苏省有超过2/3的海岸带处于淤长状态,每年新淤达到1333 km2[33]。虽然这些利用的滩涂存在含盐量过高的情况,然而随着新淤面积的不断增加,以及逐渐远离海岸与地表径流的冲洗,出露滩涂盐分随之下降。如东县作为南通市的滩涂大县,海平面以上滩涂面积达8×104 hm2,占南通市滩涂面积的50%,具备了持续生长的后备耕地资源[34]。目前改善土壤盐渍化的方法有很多,种植耐盐作物、水稻,添加土壤改良剂,进行秸秆覆盖、植被覆盖、调整耕作方式等人为措施均可达到降盐、控盐的作用[35],由此可见,如东县等沿海地区具有较强的农业发展潜力。
4. 适宜性区域分析
以农业生态环境为背景,按照南通市区域特点,参考相关行业标准依据专家经验评判确定指标权重及分级参数,建立以气候、土壤为主要指标的评价体系(详见2.2部分)。利用ArcGIS空间统计分析、叠加分析等技术为依托,通过对各指标综合分级空间计算统计(自然间断点分级法)与专家评定相结合的方式展开对南通市各县区农业生产适宜性综合性评价。由于南通市各区县气候条件相对适宜,土壤环境指标均在安全范围内,因此本研究讨论了各区县农业生产环境的相对适宜程度。启东的气候条件及土壤养分条件尚可,而土壤重金属单因子指数及内梅罗指数相对偏高且分值较低设置为一般适宜区域(四级适宜区),可种植如棉花这类重金属耐受性较好的经济作物,建立农产品加工生产基地,也可充分利用区位优势发展水产养殖业。如皋区位条件及连续性耕作,导致了该地气候及土壤养分等指标的综合分值偏低,也划分为一般适宜区范围,区内大多为砂壤土、中壤土等土质也可种植山药、马铃薯及四季蔬菜等。南通各县区中海安的综合分值最高,主要是土壤养分及土壤环境安全等指标级别较高所致,可适合大多数农作物种植。如东县与海安县情况类似,与南通市辖区划分为农业适宜生产区域也称为二级适宜区,可大规模种植四季蔬菜,因其低洼的地势条件及水源条件也发展水稻种植及水产养殖。而通州市、海门市各指标参数较为均匀,按综合分值划分为三级适宜区,因其优越的地理位置可作为上海市农产品供应基地开发,可大力发展形成蔬菜水果、水稻等需求农产品的种植基地及农业商品加工基地。
5. 结论
(1) 从气候的逐年变化情况看,南通市积温基本呈逐年递增趋势,高积温区多分布在纬度较低地区,如南通市辖区、通州、海门、启东,占据了全市面积的47.18%。降水量年际起伏波动较大,变化复杂,且多与年日照时数关系密切,其变化规律与日照时数呈负相关。
(2) 从土壤养分与质量的逐年变化看,土壤全氮与有效磷呈递增趋势,据全国第二次土壤普查土壤养分分级标准,南通市土壤全氮与有效磷处于土壤养分二级水平,速效钾含量偏低仅达到土壤养分四级要求。从养分富集程度的空间分布看,海安、如东为养分丰富区,通州、启东为养分充足地区,市辖区与海门市划为养分较充足区,如皋为养分基本充足区。
(3) 参照《土壤环境质量标准》,南通市土壤重金属含量均符合标准要求,全市综合污染指数在0.26~0.4之间,处于安全阶段。此外,相关性分析发现Pb与As之间相关性系数较高,元素同源性可能性较高。
(4) 根据农业生产适应性评价,将适宜区划分为4级,海安地区的农业生产环境达到最适宜程度(一级适宜区),南通市辖区与如东为二级农业生产适宜区,通州、海门划分为三级适宜区,启东、如皋为四级适宜区。
致谢: 野外工作得到内蒙古盛旺矿业股份有限公司全体员工的大力支持,样品测试分析环节先后得到中国冶金地质总局第一地质勘查院测试中心、中国地质科学院矿产资源研究所国家重点实验室等相关单位和个人的大力支持与帮助,在此一并表示诚挚的感谢。 -
表 1 内蒙古达亚纳石英闪长岩LA-ICP-MS锆石U-Th-Pb测试数据
Table 1 LA-ICP-MS zircon U-Th-Pb isotopic data for the quartz diorite at Dayana, Inner Mongolia
样品-测试点号 Th/ U/ Th/U 207PB/235U 207PB/235U 206PB/238U 206PB/238U 206PB/238U 206PB/238U 10-6 10-6 比值 1σ 比值 1σ 年龄/Ma 1σ ZKI3 ZK13-1973.1810.61.200.41000.01060.03750.00042372ZK13-2 208.8 165.8 1.26 0.3433 0.0119 0.0485 0.0005 306 3 ZK13-3 188.2 135.8 1.39 0.3365 0.0125 0.0486 0.0005 306 3 ZK13-41101.8459.82.400.46940.02560.05070.00103196ZK13-5 373.8 659.1 0.57 0.3378 0.0075 0.0474 0.0003 298 2 ZK13-6 644.2 391.6 1.64 0.3378 0.0083 0.0485 0.0004 305 2 ZKI3-7 467.7 245.1 1.91 0.3525 0.0105 0.0473 0.0004 298 3 ZK13-8196.3170.51.150.63620.03540.04730.00092986ZK13-9 1054.7 602.3 1.75 0.3432 0.0071 0.0474 0.0003 298 2 ZK13-10 1303.2 744.1 1.75 0.3395 0.0077 0.0475 0.0004 299 2 ZK13-113030.8927.63.272.58130.07620.06710.00094195ZK13-12714.8331.72.150.65900.01980.05220.00063284ZK13-13 457.1 602.0 0.76 0.3480 0.0068 0.0477 0.0004 300 2 ZK13-14 1407.2 584.2 2.41 0.3459 0.00X2 0.0483 0.0005 304 3 ZK13-15 885.2 781.4 1.13 0.3348 0.0064 0.0477 0.0004 301 3 ZK13-I6 699.9 422.9 1.65 0.3362 0.0077 0.Q480 0.0005 302 3 ZK13-17332.2145.92.280.52280.03380.05020.00113166ZK13-18 97.7 89.6 1.09 0.3709 0.0303 0.0488 0.0013 307 8 ZK13-19 303.7 294.8 1.03 0.3578 0.0107 0.0488 0.0006 307 3 ZK13-20 1313.4 626.1 2.10 0.3603 0.0077 0.0486 0.0005 306 3 ZK400 ZK400-1 108.9 83.3 1.31 0.3648 0.0233 0.0483 0.0009 304 6 ZK400-2 124.1 70.5 1.76 0.3593 0.0526 0.0473 0.0027 298 16 ZK400-3 133.3 121.3 1.10 0.3476 0.0107 0.0466 0.0004 293 2 ZK400-4 220.7 105.9 2.08 0.3548 0.0343 0.0487 0.0014 307 9 ZK400-5 83.2 70.4 1.18 0.3488 0.0239 0.0476 O.OOOS 300 5 ZK400-689.8390.50.233.18300.03880.21260.001812439ZK400-7101.9341.20.301.42390.01700.14950.00098985ZK400-8 152.5 152.6 1.00 0.3370 0.0116 0.0480 0.0007 302 4 ZK400-942.9234.50.181.55370.02860.14420.00118686ZK400-10 190.2 104.9 1.81 0.3315 0.0117 0.0468 0.0005 295 3 ZK400-1145.227.91.620.39460.06500.04440.001928011ZK400-12 65.8 48.3 1.36 0.3567 0.0177 0.0469 0.0007 296 4 ZK400-13 228.5 192.4 1.19 0.3476 0.0091 0.0476 0.0004 300 3 ZK400-14 137.4 103.2 1.33 0.3473 0.0119 0.0478 0.0005 301 3 ZK400-15 133.4 150.3 0.89 0.3533 0.0119 0.0482 0.0005 303 3 ZK400-16 98.8 57.6 1.72 0.3547 0.0199 0.0471 0.0008 297 5 ZK400-17435.5200.42.171.04750.02370.04870.00053063ZK400-18 158.0 78.7 2.01 0.3382 0.0131 0.0476 0.0005 300 3 ZK400-19 191.3 93.4 2.05 0.3405 0.0192 0.0478 0.0008 301 5 ZK400-20 448.0 337.3 1.33 0.3383 0.0112 0.0471 0.0005 297 3 注:带删除线数据表示测试数据有问题,不可利用 表 2 葛根敖包石英闪长岩主量、微量和稀土元素分析数据
Table 2 Major, trace elements and REE composition of the quartz diorite at the Gegen'aobao
元素 GGAB-ZK301-243 GGAB-ZK301-332 GGAB-ZK001-400 GGAB-ZK001-177 GGAB-ZK001-13 GGAB-ZK008-205 GGAB-ZK008-210 GGAB-ZK008-394 SiO2 65.01 64.95 64.37 65.84 60.13 64.43 64.62 63.15 Al2O3 15.93 15.81 16.14 16.11 15.30 15.85 15.24 15.80 Fe2O3 1.47 1.05 1.16 1.18 2.18 1.29 1.21 1.62 CaO 2.74 3.11 3.22 2.67 5.37 3.78 3.85 3.39 MgO 1.39 1.36 1.22 1.09 2.07 1.55 1.46 1.67 K2O 3.88 2.77 3.88 3.80 3.19 3.13 3.10 3.24 Na2O 4.62 5.34 4.43 4.91 3.71 4.45 4.86 4.64 MnO 0.13 0.14 0.15 0.12 0.26 0.13 0.21 0.34 P2O5 0.24 0.28 0.27 0.22 0.28 0.27 0.26 0.27 TiO2 0.99 1.03 0.98 0.94 1.05 0.93 0.95 0.95 FeO 3.22 3.10 3.52 3.10 4.23 3.81 3.81 3.46 烧失量 1.00 1.63 1.14 0.73 2.38 0.99 0.35 1.45 总计 100.61 100.57 100.46 100.71 100.14 100.60 99.92 99.97 A/CNK 0.95 0.91 0.93 0.95 0.79 0.90 0.83 0.91 Na2O+K2O 8.49 8.11 8.30 8.71 6.90 7.58 7.96 7.88 N2O/K2O 1.19 1.93 1.14 1.29 1.16 1.42 1.57 1.43 La 20.31 15.71 19.53 18.95 20.87 24.56 22.17 27.54 Ce 51.09 41.67 47.51 49.62 45.77 54.30 49.27 57.47 Pr 8.74 7.17 8.03 8.66 7.53 8.42 7.87 8.64 Nd 36.49 31.19 33.44 36.26 29.89 34.10 32.06 33.54 Sm 8.31 7.38 7.64 8.18 6.76 7.24 7.10 7.13 Eu 2.25 2.17 2.13 2.08 1.87 1.97 1.90 2.00 Gd 8.21 6.91 7.03 7.50 6.19 6.90 6.53 6.71 Tb 1.21 1.07 1.07 1.18 0.98 1.07 1.02 0.98 Dy 6.98 5.98 6.40 6.74 5.50 5.99 5.64 5.61 Ho 1.35 1.21 1.21 1.36 1.10 1.19 1.12 1.11 Er 3.88 3.44 3.44 3.89 3.07 3.27 3.19 3.00 Tm 0.57 0.50 0.55 0.59 0.46 0.50 0.47 0.45 Yb 3.65 3.23 3.37 3.86 3.09 3.23 3.14 3.04 Lu 0.54 0.45 0.50 0.57 0.44 0.50 0.46 0.46 SREE 153.58 128.06 141.85 149.42 133.52 153.23 141.94 157.68 LREE 127.19 105.28 118.27 123.74 112.69 130.59 120.36 136.33 HREE 26.38 22.78 23.58 25.68 20.83 22.64 21.57 21.35 LREE/HREE 4.82 4.62 5.02 4.82 5.41 5.77 5.58 6.38 LaN/YbN 3.99 3.49 4.15 3.52 4.85 5.46 5.07 6.51 δEu 0.82 0.91 0.87 0.80 0.87 0.84 0.84 0.87 δCe 0.94 0.96 0.93 0.95 0.89 0.92 0.91 0.91 Li 10.3 17.0 11.1 16.2 14.8 17.4 10.7 15.1 Be 2.61 1.77 2.04 2.67 1.47 1.89 1.94 1.95 V 64.93 69.24 69.13 62.65 127.80 83.89 77.59 76.44 Cr 26.01 18.07 22.87 19.09 22.13 16.92 17.19 16.09 Co 5.01 5.09 4.92 4.10 10.0 7.70 5.83 7.86 Ni 11.85 2.28 1.79 1.73 5.28 2.51 2.38 3.26 Cu 9.68 14.1 35.2 8.40 31.1 9.73 13.5 15.0 Zn 49.67 70.33 115.90 77.70 171.10 63.98 68.49 66.88 Ga 19.73 18.76 18.75 18.87 19.92 18.68 18.94 18.63 Rb 77.65 55.95 86.50 83.34 67.68 66.27 60.49 67.22 Sr 404.80 368.30 288.60 393.30 432.10 461.20 435.20 490.80 Y 29.14 25.46 25.84 28.74 23.08 25.51 24.13 24.22 Zr 247.70 232.08 238.77 263.69 150.60 216.58 202.76 208.08 Nb 14.01 12.73 11.26 13.52 9.72 10.88 10.97 11.01 Cs 2.25 1.96 1.65 7.92 2.21 6.83 9.78 2.54 Hf 5.31 4.91 4.98 7.00 4.43 5.55 5.28 4.78 Ba 977.70 806.50 773.10 978.30 620.30 747.40 765.00 873.80 Ta 1.11 0.90 0.82 0.92 0.74 0.74 0.76 0.81 Pb 12.27 12.48 40.01 16.87 60.56 10.11 14.95 14.04 Th 9.57 7.85 8.16 9.18 6.42 7.01 6.82 7.01 U 2.95 2.33 2.65 3.13 2.19 2.53 2.38 2.30 Nb/Ta 12.62 14.11 13.71 14.66 13.19 14.72 14.49 13.64 Rb/Sr 0.19 0.15 0.30 0.21 0.16 0.14 0.14 0.14 注:主量元素含量单位为%,微量和稀土元素含量单位为10-6 表 3 葛根敖包石英闪长岩锆石Hf同位素分析数据
Table 3 Zircon Lu-Hf isotopic results for the quartz diorite at the Gegen'aobao deposit
测点号 t/Ma 176Yb/177Hf 176Lu/177Hf 176Hf/177Hf 2σ εHf(0) εHf(t) TdM1 TdM2 fu-Hf ZK13-02 306 0.057253 0.001956 0.282948 0.000018 6.2 12.6 442 517 -0.94 ZK13-03 306 0.050751 0.001746 0.282940 0.000018 5.9 12.3 451 533 -0.95 ZK13-05 298 0.070979 0.002803 0.282949 0.000022 6.3 12.3 451 530 -0.92 ZK13-06 305 0.088795 0.003118 0.282906 0.000020 4.7 10.8 520 629 -0.91 ZK13-07 298 0.134357 0.004649 0.282922 0.000024 5.3 10.9 519 616 -0.86 ZK13-09 298 0.094177 0.003372 0.282925 0.000025 5.4 11.3 495 593 -0.90 ZK13-10 299 0.086160 0.003107 0.282935 0.000018 5.7 11.7 477 567 -0.91 ZK13-13 300 0.050320 0.002030 0.282939 0.000018 5.9 12.1 456 542 -0.94 ZK13-14 304 0.055622 0.002038 0.282960 0.000017 6.6 12.9 425 493 -0.94 ZK13-15 301 0.035758 0.001459 0.282914 0.000013 5.0 11.3 485 592 -0.96 ZK13-16 302 0.037999 0.001492 0.282942 0.000015 6.0 12.4 445 528 -0.96 ZK13-18 307 0.072069 0.002644 0.282946 0.000020 6.1 12.4 454 531 -0.92 ZK13-19 307 0.063193 0.002407 0.282916 0.000017 5.1 11.4 495 596 -0.93 ZK13-20 306 0.035327 0.001413 0.282925 0.000013 5.4 11.8 469 564 -0.96 ZK400-01 304 0.034900 0.001294 0.282961 0.000019 6.7 13.1 415 480 -0.96 ZK400-02 298 0.048393 0.001648 0.282920 0.000021 5.2 11.5 478 581 -0.95 ZK400-03 293 0.029959 0.001137 0.282960 0.000019 6.7 12.9 415 487 -0.97 ZK400-04 307 0.037745 0.001360 0.282906 0.000021 4.7 11.2 495 605 -0.96 ZK400-05 300 0.031210 0.001135 0.282907 0.000016 4.8 11.1 491 604 -0.97 ZK400-08 302 0.035079 0.001330 0.282914 0.000016 5.0 11.4 483 589 -0.96 ZK400-10 295 0.038373 0.001500 0.282961 0.000018 6.7 12.9 418 489 -0.95 ZK400-12 296 0.055173 0.002088 0.282985 0.000017 7.5 13.6 389 440 -0.94 ZK400-13 300 0.067177 0.002427 0.282955 0.000017 6.5 12.6 437 511 -0.93 ZK400-14 301 0.048535 0.001808 0.282961 0.000017 6.7 12.9 421 489 -0.95 ZK400-16 297 0.071620 0.002692 0.282984 0.000022 7.5 13.5 398 451 -0.92 ZK400-18 300 0.058848 0.002119 0.282945 0.000020 6.1 12.3 448 530 -0.94 ZK400-19 301 0.053619 0.001883 0.282936 0.000018 5.8 12.0 459 548 -0.94 ZK400-20 297 0.065200 0.002324 0.282922 0.000019 5.3 11.4 485 587 -0.93 -
Wang T, Guo L, Zhang L, et al. Timing and evolution of JurassicCretaceous granitoid magmatisms in the Mongol-Okhotsk belt and adjacent areas, NE Asia:implications for transition from contractional crustal thickening to extensional thinning and geodynamic settings[J]. Journal of Asian Earth Sciences, 2015, 97:365-392. doi: 10.1016/j.jseaes.2014.10.005
Xu B. The Central Asian Orogenic Belt in northern China:Preface. Journal of Asian Earth Sci., 2015, 97:179-182. doi: 10.1016/j.jseaes.2014.11.019
Xu B, Song S, Nie F. The Central Asian Orogenic Belt in northern China:Preface[J]. Journal of Asian Earth Sciences, 2015, 97:179-182. doi: 10.1016/j.jseaes.2014.11.019
Xu B, Zhao P, Wang Y Y, et al. The pre-Devonian tectonic framework of Xing'an-Mongolia orogenic belt (XMOB) in north China. Journal of Asian Earth Sciences. 2015c, 97:183-196. doi: 10.1016/j.jseaes.2014.07.020
Si H J, Bagas L, Peng H, et al. Muscovite Ar-Ar, molybdenite ReOs and zircon U-Pb ages and Sr-Nd-Hf isotopes of the highly fractionated granite-related Shamai tungsten deposit in East Inner Mongolia, China:implications for the timing of mineralization and ore genesis[J]. Lithos, 2015, in publication. https://www.sciencedirect.com/science/article/pii/S0024493715004739
Mao J W, Zhang Z C, Zhang Z H, et al. Re-Os isotopic dating of molybdenites in the Xiaoliugou W (Mo) deposit in the Northern Qilian Mountains and its geological significance[J]. Geochimca et Cosmochimica Acta, 1999, 63(11/12):1815-1818. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=eb7204df38315b7dd7fcab6bc9a0e8bf
Mao J W, Wang Y T, Zhang Z H, et al. Geodynamic settings of Mesozoic large-scale mineralization in North China and adjacent areas[J]. Science China Earth Sciences, 2003, 46(8):839-851. http://cn.bing.com/academic/profile?id=0ed80e14c84d8dcb8f89cfe69ac6e89e&encoded=0&v=paper_preview&mkt=zh-cn
聂凤军, 胡朋, 江思宏, 等.中蒙边境沙麦-玉古兹尔地区钨和钨(钼)矿床地质特征, 形成时代和成因机理[J].地球学报, 2010, 3:383-394. http://d.old.wanfangdata.com.cn/Periodical/dqxb201003012 佘宏全, 李进文, 向安平, 等.大兴安岭中北段原岩锆石U-Pb测年及其与区域构造演化关系[J].岩石学报, 2012, 28(2):571-594. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201202018 向安平, 杨郧城, 李贵涛, 等.黑龙江多宝山斑岩Cu-Mo矿床成岩成矿时代研究[J].矿床地质, 2012, 6:1237-1248. doi: 10.3969/j.issn.0258-7106.2012.06.009 向安平, 佘宏全, 陈毓川, 等.内蒙古红花尔基钨钼矿云英岩化白云母Ar-Ar定年及其地质意义[J].岩矿测试, 2016, 1:108-116. http://d.old.wanfangdata.com.cn/Periodical/ykcs201601018 Li W, Hu C, Zhong R, et al. U-Pb, 39Ar/40Ar geochronology of themetamorphosed volcanic rocks of the Bainaimiao Group in centralInner Mongolia and its implications for ore genesis and geodynamic setting[J]. Asian Earth Journal of Sciences, 2015, 97:251-259. doi: 10.1016/j.jseaes.2014.06.007
Wu G, Chen Y C, Sun F Y, et al. Geochronology, geochemistry, and Sr-Nd-Hf isotopes of the early Paleozoic igneous rocks in the Duobaoshan area, NE China, and their geological significance[J]. Journal of Asian Earth Sciences, 2015, 97:229-250. doi: 10.1016/j.jseaes.2014.07.031
Wang F, Xu W L, Xu Y G, et al. Late Triassic bimodaligneous rocks in eastern Heilongjiang Province, NE China:implications for theinitiation of subduction of the Paleo-Pacific Plate beneath Eurasia[J]. Journal of Asian Earth Sciences, 2015, 97:406-423. doi: 10.1016/j.jseaes.2014.05.025
Wang Z Z, Han B F, Feng L X, et al Geochronology, geochemistry andorigins of the Paleozoic-Triassic plutons in the Langshan area, western InnerMongolia, China[J]. Journal of Asian Earth Sciences, 2015, 97:337-351. doi: 10.1016/j.jseaes.2014.08.005
Wang Z W, Pei F P, Xu W L, et al. Geochronology andgeochemistry of Late Devonian and early Carboniferous igneous rocks of centralJilin Province, NE China:implications for the tectonic evolution of the eastern Central Asian Orogenic Belt[J]. Asian Earth Sci., 2015, 97:260-278. doi: 10.1016/j.jseaes.2014.06.028
Kravchinsky V A, Konstantinow K M, Cogne J P. Palaeomagnetic study of Vendian and Early Cambrian rocks of South Siberia and Central Mongolia:was the Siberian platform assembled at this time?[J]. Precambrian Research, 2001, 110:61-92. doi: 10.1016/S0301-9268(01)00181-4
Cogne J P, Kravchinsky V A, Halim N, et al. Late Jurassic-Early Cretaceous closure of the Mongol-Okhoysk Ocean demonstrated by new Mesozoic palaeomagnetic results from the Trans-Baikal area (SE Siberian)[J]. Geophysical Journal International, 2005, 163:823-832. http://cn.bing.com/academic/profile?id=5826d47d1b5ba9424da841a6150ab675&encoded=0&v=paper_preview&mkt=zh-cn
Metelkin D V, Gordienko I V, Klimuk V S. Paleomagnetism of Upper Jurassic basalts from Transbaikalia:new data on the time of closure of the Mongol-Okhotsk Ocean and Mesozoic intraplate tectonics of Central Asia[J]. Russian Geology and Geophysics, 2007, 48:825-834. doi: 10.1016/j.rgg.2007.09.004
宋彪, 张玉海, 刘敦一.微量原位分析仪器SHRIMP的产生与锆石同位素地质年代学[J].质谱学报, 2002, (1):58-62. doi: 10.3969/j.issn.1004-2997.2002.01.011 Liu Y S, Hu Z C, Gao S, et al. In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard[J]. Chemical Geology, 2008, 257:34-43. doi: 10.1016/j.chemgeo.2008.08.004
Liu Y S, Gao S, Hu Z C, et al. Continental and oceanic crust recycling-induced melt-peridotite interactions in the TransNorth China Orogen:U-Pb dating, Hf isotopes and trace elements in zircons of mantle xenoliths[J]. Journal of Petrology, 2010, 51:537-571. doi: 10.1093/petrology/egp082
Liu Y S, Hu Z C, Zong K Q, et al. Reappraisement and refinement of zircon U-Pb isotope and trace element analyses by LA-ICP-MS[J]. Chinese Science Bulletin, 2010, 55:1535-1546. doi: 10.1007/s11434-010-3052-4
侯可军, 李延河, 田有荣. LA-MC-ICP-MS锆石微区原位U-Pb定年技术[J].矿床地质, 2009, 4:481-492. doi: 10.3969/j.issn.0258-7106.2009.04.010 Sláma J, Kosler J, Condon D J.Plesovice zircon-A new natural reference material for U-Pb and Hf isotopic microanalysis[J]. Chemical Geology, 2008, 249:1-35. doi: 10.1016/j.chemgeo.2007.11.005
侯可军. LA-MC-ICP-MS锆石Hf同位素的分析方法及地质应用[J].岩石学报, 2007, 23(10):2595-2604 doi: 10.3969/j.issn.1000-0569.2007.10.025 Morel M L A, Nebel O, Nebel-Jacobsen Y J. Hafnium isotope characterization of the GJ-1 zircon reference material by solution and laser-ablation MC-ICPMS[J]. Chemical Geology, 2008, 255:231-235. doi: 10.1016/j.chemgeo.2008.06.040
HofmannAW. Chemical differentiation of the earth:The relationship between mantle, continental crust, and oceanic crust[J]. Earth Planetary Science Letters, 1988, 90(3):297-314. doi: 10.1016/0012-821X(88)90132-X
Kelemen P B, H Angh K, Greenem A R. One view of thegeochemistry of subduction-related magmat ic arcs, with anemphasison primitive andesite and lowercrust[J]. Treatise On Geochemistry, 2003, 3:593-659. https://www.sciencedirect.com/science/article/pii/B0080437516030358
Zhou J B, Wang B, Wilde S A. Geochemistry and U-Pb zircon dating of the Toudaoqiao blueschists inthe Great Xing an Range, northeast China, and tectonic implications[J]. Journal of Asian Earth Sciences, 2015, 97:197-210. doi: 10.1016/j.jseaes.2014.07.011
Han G Q, Liu Y J, Neubauer F. U-Pb age and Hf isotopic data of detrital zircons from the Devonian andCarboniferous sandstones in Yimin area, NE China:New evidences to thecollision timing between the Xing'an and Erguna blocks in eastern segment of Central Asian Orogenic Belt[J]. Journal of Asian Earth Sciences, 2015, 97:211-228. doi: 10.1016/j.jseaes.2014.08.006
Martin H, Smithies R H, Rapp R, et al. An overview of adakite, tonalite-trondhjemite-granodiorite (TTG), and sanukitoid:Relationships and some implications for crustal evolution[J]. Lithos, 2005, 79(1/2):1-24. http://cn.bing.com/academic/profile?id=7371dfa5724c0893c0b1241c8ad6e418&encoded=0&v=paper_preview&mkt=zh-cn
陈斌, 赵国春, Wi lde S A.内蒙古苏尼特左旗南两类花岗岩同位素年代学及其构造意义[J].地质论评, 2001, 47(4):361-367. doi: 10.3321/j.issn:0371-5736.2001.04.005 鲍庆中, 张长捷, 吴之理, 等.内蒙古白音高勒地区石炭纪石英闪长岩SHRIM P锆石U-Pb年代学及其意义[J].吉林大学学报(地球科学版), 2007, 37(1):15-23. http://www.cnki.com.cn/Article/CJFDTotal-CCDZ200701002.htm 鲍庆中, 张长捷, 吴之理, 等.内蒙古东南部晚古生代裂谷区花岗质岩石锆石SHRIMP U-Pb定年及其地质意义[J].中国地质, 2007, 34(5):790-798. doi: 10.3969/j.issn.1000-3657.2007.05.005 刘建峰, 迟效国, 张兴洲, 等.内蒙古西乌旗南部石炭纪石英闪长岩地球化学特征及其构造意义[J].地质学报, 2009, 3:365-376. http://d.old.wanfangdata.com.cn/Periodical/dizhixb200903006 Zheng Y F, Zhao Z F, Wu F B, et al. Zircon U-Pb age, Hf and O isotope constrains on protolith origin of ultrahigh-pressure eclogite and gneiss in the Dabie orogeny[J]. Chem. Geol., 2006, 231:135-158. doi: 10.1016/j.chemgeo.2006.01.005
郑永飞, 陈仁旭, 张少兵, 等.大别山超高压榴辉岩和花岗片麻岩中锆石Lu-Hf同位素研究[J].岩石学报, 2007, 2:317-330. http://d.old.wanfangdata.com.cn/Periodical/ysxb98200702012