LA-ICP-MS zircon U-Pb dating and geochemical characteris-tics of Baoxing rock mass in Xifeng area of northern Liaoning Province
-
摘要:
宝兴岩体位于华北板块北缘陆缘活动带东段,岩性为花岗闪长岩,LA-ICP-MS锆石206Pb/238U年龄加权平均值为235.0±1.3Ma,属于中三叠世。岩石的K2O/Na2O平均值为1.07,属高钾钙碱性系列;富集大离子亲石元素,亏损高场强元素,具有活动陆缘钙碱性岩系特征;稀土元素总量平均值为147×10-6,稀土元素配分曲线右倾,属具中等负Eu异常的轻稀土元素富集、重稀土元素亏损型。通过岩石学、岩石地球化学及年代学研究,结合邻区构造演化研究,认为中三叠世研究区在伸展作用下岩石圈进行拆沉,从而形成大规模的岩浆活动。
Abstract:Located on the northern margin of North China Block, the Baoxing rock mass is mainly composed of granodiorite. According to LA-ICP-MS zircon U-Pb dating, it was formed in Middle Triassic, with the 206Pb/238U age of 235.0±1.3Ma. The average K2O/Na2O ratio is 1.07, indicating that the rocks belong to high-K calc-alkaline series. The rocks are enriched in LILE and depleted in HFSE, showing typical features of calc-alkaline series on an active continental margin. The average total amount of REE is 147×10-6, and the curve exhibits obvious right-inclined feature on REE distribution pattern, suggesting that the rocks are of the type of LREE enrichment and HREE depletion with medium negative Eu anomaly. According to the petrology, geochemistry and geochronology in combination with the tectonic evolution research on adjacent areas, lithospheric delamination occurred under the extension in the study area, and then caused a large-scale magmatic activity in the Middle Triassic.
-
矿产开发不仅由于矿渣、矸石、尾矿等的堆积,造成对土地的占压,同时由于采空区的形成产生大量的地裂缝。地裂缝的发生对其周围4m范围内的植被生长有较大的影响[1-3]。地裂缝区域的土壤水分常作为衡量其对土壤环境影响的重要指标。地裂缝的发育,改变了区域土体的结构、微环境及其水文循环过程,植被生长受到水分胁迫机理研究受到学者的广泛关注[4-8]。赵红梅等[9-11]对地裂缝一维纵向剖面上的含水率分布特征进行了研究,表明地裂缝对土壤水分的影响,杨泽元等、王文科等[12-14]就风沙滩地区采煤塌陷裂缝对包气带水分运移的影响模型进行了研究,李海欣等[1-2]就煤矿开采沉陷地裂缝对植被覆盖的影响等进行了研究。但是关于地裂缝对土壤水分的影响范围、程度及发生机制研究则较少。明确矿采对区域生态环境的影响、实现地区资源开发同生态环境保护的协调发展、防治土地沙化受到国家和地方政府的重视[15-21]。
宁东煤炭能源基地总面积约3484km2,地处毛乌素沙漠腹地,降水稀少,蒸发强烈,地质环境极其脆弱,土壤水分是制约植被生长及其迹地植被恢复的重要瓶颈。
本次以宁东煤矿区为研究区域,选择发育于不同土壤母质上的地裂缝,进行地裂缝在二维剖面上的土壤含水率时空演变特征及不同土壤母质对地裂缝的影响研究,为宁东矿采区资源开发与生态环境保护的协调发展,以及防治土地沙化提供理论依据。
1. 研究方法
本文以宁东煤矿区的羊场湾煤矿采空区作为研究场所,地理坐标为北纬37°59'19"、东经106°35' 29",海拔1361m,地貌类型为冲洪积平原。植物优势种群以沙蒿、柠条为主,植被覆盖度30%~50%。
采用筛分法和密度计相结合的方法进行试验区土壤颗粒分级,以土壤岩性组分进行土壤命名。土壤剖面自上而下可分为4层:0~0.1m为风积沙,含根系;0.1~0.2m为风积沙;0.2~0.5m为粉土;0.5~ 1.0m为黄土。
选择地形平整、植被生长较好、均一的采空地裂缝区域为研究区,以地裂缝为起点,垂直地裂缝方向的二维剖面土壤水分为研究对象,地裂缝的宽度为2~3m,沿垂直地裂缝的方向开挖长3.6m、宽0.9m、深1.0m的浅井,土壤二维含水率监测平面图见图 1。在其一侧的剖面垂直方向由垂向0.0m、0.1m、0.2m、0.3m、0.5m、0.7m、1.0m土层和沿地裂缝垂直水平方向的0.0m、0.3m、0.6m、0.9m、1.2m、1.5m、1.8m、2.1m、2.4m、2.7m、3.0m、3.3m和3.6m处分别设置土壤含水率数据采集点,土壤二维含水率监测点剖面分布见图 2。数据采集仪器设备为英国Delta-T公司生产的土壤水分、温度、电导率便携式测量仪,测量精度约3.0%,分辨率0.1%。采用电导原理测量土壤介质的体积含水率,每次数据采集完后对浅井进行回填处理,以减小开挖浅井造成的土壤水分散失,数据采集时间为2017年6~10月,每月1次。
为了避免大气降水对观测的影响,每次数据采集时间选择晴朗的天气进行。
收集附近宁夏灵武市气象站资料,包括降水量、蒸发量、气温等气候要素,作为影响因子。数据处理采用GMS软件中的Interpolation模块,作图采用Excel2016、Coreldraw13.0软件。
2. 研究结果
2.1 气象特征
宁东煤矿区地处毛乌素沙漠腹地,气温呈现单峰季节变化特征,年平均气温7.5℃,最高气温29.3℃,出现于植物生长季7月,最低气温-14.6℃,出现于冬季1月。1~7月为持续升温期,7~12月为持续降温期,区内平均气温及降水量、蒸发量随时间变化曲线见图 3。
降水亦呈现单峰季节变化特征,年降水量284mm,生长季5~9月月降水量占年降水量的83%,8月的降水量最大,达72.6mm。年平均蒸发量达1723mm,是降水量的6.1倍。蒸发量的季节动态亦呈单峰变化趋势,其峰值出现点提前于降水2个月。
2.2 含水率时空特征
为了查明土-气界面地表、地裂缝侧向裸露面对土壤水分散失的影响程度,本文绘制了不同月份土壤剖面含水率等值线图,结合气象要素进行分析,以总结土壤剖面含水率时空变化规律,6~10月土壤剖面含水率分布等值线图见图 4—图 8。
受裸露的地裂缝断面水分蒸发的影响,在越靠近地裂缝前缘的部分,其土壤含水率越小,在裸露的地裂缝断面上,其含水率接近残余含水率,水平方向0~1.2m及垂直方向0.4~1.0m的范围属于水分主要散失区,且主要分布在粉土、黄土区段(图 4),表明6月土壤水分主要通过侧面裂缝的裸露面散失,且靠近地表范围基本不受影响,表生植被生长没有来自水分缺失胁迫的影响。6月气温逐步升高,月平均气温20.9℃,植被进入复苏期,根系耗水能力逐步增强,在蒸发作用下,土壤水分略有损耗,但有降雨及时补充。
随着7月气温的骤然升高,土壤水分散失的范围也急剧扩大,在6月的基础上,水平方向大于3.6m及垂直方向大于1.0m的范围属于水分主要散失区(图 5),表明7月受外界气温变化的影响,土壤水分的散失范围扩展到地表,通过地表、侧面裂缝的裸露面散失,表生植被生长受到来自水分缺失胁迫的影响。
风积沙下部的粉土、黄土颗粒较细,渗透性能较差,严重阻碍了大气降雨由风积沙层进入粉土层的入渗速率,且风积沙的渗透性能好,大气降水短时间内滞留在该层,有利于植被根系吸收水分;同时,细颗粒物质的比表面积增加,蒸发面扩大,其水分也更容易蒸发散失。
7月气温快速升高,月平均气温22.7℃,植被进入旺盛生长期,根系耗水能力进一步增强,且在腾发作用下,土壤亏缺大量水分,尽管7月降水量达60.1mm,但依旧无法补充土壤水分的亏缺量。
随着8月气温的略微回落,土壤水分散失的范围也逐步收缩,在7月的基础上,水平方向0~3m及垂直方向0.2~1.0m的范围属于水分主要散失区(图 6),表明8月受外界气温回落的影响,土壤水分的散失范围收缩到地裂缝的裸露面,主要通过侧面裂缝的裸露面散失,表生植被生长略微受到来自水分缺失胁迫的影响。
8月气温略有回落,月平均气温20.8℃,植被进入生长期,根系耗水能力持续增强,但本月降水量增加,达72.6mm,土壤水分得到及时补充。
随着9月气温的持续回落,加之8月的降水量增加,土壤水分散失的范围也逐步收缩,在8月的基础上,水平方向0~1.8m及垂直方向0.2~1.0m的范围属于水分主要散失区(图 7),表明8月受外界气温持续回落的影响,土壤水分的散失范围继续收缩到地裂缝的裸露面,主要通过侧面裂缝的裸露面散失,表生植被生长不会受水分缺失胁迫的影响。
9月气温降低,月平均气温仅15.2℃,植被进入衰退期,根系耗水能力随之下降,土壤水分得到有效补充。
随着10月气温的大幅回落,土壤水分散失的范围也急剧收缩,在9月的基础上,水平方向0~1.5m及垂直方向上接近0.4~1.0m的范围属于水分主要散失区(图 8),表明10月受外界气温大幅回落的影响,土壤水分的散失范围继续收缩到地裂缝的裸露面,主要通过侧面裂缝的裸露面散失,且水平方向的范围也在收缩,表生植被生长不会受水分缺失胁迫的影响。
10月气温再次降低,月平均气温仅8.4℃,植被进入休眠期,根系耗水能力随之消失,土壤水分进一步得到补充。
2.3 影响因素
研究表明,地裂缝的侧向裸露面、土壤类型及结构、气温、埋深等是土壤水分散失的主要影响要素。
(1)地裂缝侧向裸露面
土壤水分蒸发的主要驱动要素即是外界气象因子,土壤蒸发的水分是自土-气界面散失,界面的大小及形态直接影响土壤水分散失的强度,地裂缝侧向裸露面的存在增加了土-气界面的面积,且自上而下垂直地表的存在形态,加速了靠近地裂缝裸露面土壤水分的散失,这对于土壤水分保持极为不利。
(2)土壤类型
粉土属于细颗粒物质,一方面风积沙+粉土上粗下细的土壤结构,使粉土大部分降水入渗水量拦截保存在近地表的风积沙中,更有利于表生植被生长;另一方面细颗粒物质的比表面积增加,蒸发面扩大,其水分也更容易蒸发散失。
(3)土壤结构
土壤水分散失空间上以下部的细颗粒物质区域为主,上部0.1~0.2m属于水分储存区,因此,该结构地裂缝的存在对植被生长所需水分的影响较小。
(4)气温
气温提供土壤水分散失所需的能量,相同条件下,气温同表土蒸发量呈正相关关系,气温越高,土壤水分散失量更大,反之亦然。
(5)埋深
鉴于土壤蒸发的水分是自土-气界面散失,因此,在远离地裂缝裸露面的情况下,外界太阳辐射对土壤水分的蒸发作用随着埋深的增加而逐步减小,在靠近裸露面的情况下,外界太阳辐射对土壤水分的蒸发作用来自地表和地裂缝侧向裸露面,埋深越大、距离裸露面越远,其外界太阳辐射对土壤水分的蒸发作用越弱,反之亦然。
3. 讨论
土壤含水率的时空变化受外界气象要素、土壤岩性、埋深等多种要素制约,6~10月剖面土壤含水率的观测数据,具有如下几个方面的特征。
(1)将同种土壤岩性含水率监测数据按照同一月份取平均值,绘制出不同土壤岩性含水率随时间变化曲线(图 9),从空间上看,随着土壤埋深的增加,土壤含水率呈现由小(根植土)变大(风积沙)、再由大变小(粉土)、再由小变大(黄土)的“S”形特征变化,不同岩性的含水率大小排序是:风积沙>黄土>根植土>粉土;从时间上看,6~10月期间,7月气温最高,土壤含水率最小(图 3)。
(2)利用6~10月的全部观测数据取平均值进行处理,分为13个监测截面,绘制不同监测截面埋深含水率变化曲线(图 10),不同监测截面土壤含水率自上而下呈现出变化一致、且相互平行的变化特征;随着截面向地裂缝靠近,土壤水分呈现出逐步减小的变化特征,表明地裂缝侧面裸露面的蒸发作用对土壤水分散失起主要作用。
(3)将不同月份所有含水率监测数据取平均值,建立土壤含水率同平均气温、降水量的相关关系(图 11),土壤含水率同降水量呈较差的负相关关系,相关系数R2是0.24,表明在宁东煤矿区降水量对土壤含水率变化的影响较小;土壤含水率与平均气温呈很好的负相关关系,相关系数R2是0.97,与实际情况吻合,表明宁东煤矿区气温对土壤含水率变化的影响较大。
4. 结论
(1)土壤水分散失空间上水分优先自裂缝裸露面散失,其次才会自地表散失,且增加的土壤水分散失面加剧了表土水分散失的程度和范围,与裂缝发育的宽度、长度以及走向均有关系;
(2)土壤水分散失时间上呈动态变化,同气温呈正相关关系。其中,7月最大,呈现地表及裂缝侧向裸露面2个蒸发面,地表土壤水分散失加剧,土壤水分散失水平影响范围大于3.6m;8月份次之,随着8月降水量的适当增加,地表的蒸发面逐渐消失,但是裂缝区裸露的蒸发面依然散失水分,呈现1个蒸发面,土壤水分散失水平影响范围为3m,对土壤水分保持造成极不利的影响。
(3)宁东煤矿区气温对土壤含水率变化的影响较大;受土壤岩性的影响,随着土壤埋深增加,土壤含水率呈现“S”形变化特征,不同岩性持水性大小顺序为风积沙>黄土>根植土>粉土。
(4)地裂缝的存在使土壤水分呈现出二维、三维的迁移特征,相关机理研究需要考虑土壤水分的侧向流动。
致谢: 野外工作得到中国地质调查局沈阳地质调查中心邴志波教授级高工、辽宁省地质矿产勘察院卢崇海研究员、辽宁省地质矿产调查院同事的指导,审稿专家提出宝贵修改意见,实验测试过程得到中国地质科学院矿产资源研究所国土资源部成矿作用与资源评价重点实验室工作人员的指导,在此一并致谢。 -
图 1 研究区地质简图(大地构造位置图据参考文献①修改)
K1d-德仁组火山岩;γπK1-早白垩世花岗斑岩;ρηγJ2-中侏罗世似斑状二长花岗岩;γJ1-早侏罗世黑云母花岗岩;ηγJ1-早侏罗世二长花岗岩;γδT2-宝兴岩体;1-研究区范围;2-同位素采样位置;3-韧性剪切带;4-断层;Ⅰ1-准格尔-兴蒙陆缘造山带;Ⅰ2-锡林浩特微板块;Ⅰ3-松嫩-佳木斯微板块;Ⅱ1-晚古生代陆缘造山带;Ⅱ2-中元古代陆缘造山带;Ⅱ3-华北地块
Figure 1. Geological map of the study area
表 1 宝兴岩体主量、微量和稀土元素分析结果
Table 1 The analyses of major, trace and rare earth elements of Baoxing rock mass
样品编号 PM023-4-TY1 PM023-11-TY1 PM023-11-TY2 PM023-16-TY1 PM023-20-TY1 平均值 SiO2 68.24 74.06 52.36 57.02 60.98 62.53 TiO2 0.58 0.16 0.82 1.11 0.65 0.66 Al2O3 15.12 14.5 15.97 17.34 15.7 15.73 Fe2O3 2.69 1.27 3.83 3.3 2 2.62 FeO 1.35 0.58 3.86 5.02 3.6 2.88 MnO 0.18 0.05 0.2 0.18 0.12 0.15 MgO 1.64 0.61 3.13 3.66 3.28 2.46 CaO 0.76 0.44 5.45 2.97 4.43 2.81 Na2O 3.48 3.00 0.76 3.96 3.55 2.95 K2O 3.2 3.95 4.27 1.86 2.51 3.16 P2O3 0.18 0.04 0.15 0.24 0.22 0.17 烧失量 2.55 1.32 8.37 3 2 3.45 合计 99.97 99.98 99.17 99.66 99.04 99.57 K 26553 32776 35431 15434 20827 26204 Rb 107 109 179 62.2 62 103.84 Ba 588 500 432 886 868 654.8 Th 10.1 8.05 6.18 8.86 5.23 7.68 Ta 0.73 0.37 0.47 0.73 0.46 0.55 Nb 11 8.52 6.86 17.5 13.9 11.56 Zr 112 100 116 210 257 159 Hf 3.58 3.71 4.33 7.3 6.79 5.14 Sr 234 208 354 431 675 380.4 La 36.9 27 19.9 38.5 33.3 31.12 Ce 66.1 45 40.9 77.1 67.3 59.28 Pr 8.64 5.24 5.03 10.3 8.41 7.52 Nd 31.8 16.1 17.9 39.8 33.8 27.88 Sm 5.97 2.63 3.80 7.61 6.99 5.4 Eu 1.56 0.72 1.17 1.78 1.67 1.38 Gd 5.37 2.39 3.32 6.74 5.98 4.76 Tb 0.77 0.3 0.58 1.03 0.86 0.71 Dy 4.47 1.6 3.46 5.65 4.85 4.01 Ho 0.86 0.29 0.72 1.27 0.92 0.81 Er 2.32 0.82 1.98 3.41 2.64 2.23 Tm 0.35 0.1 0.32 0.56 0.43 0.35 Yb 1.97 0.71 1.97 3.41 2.54 2.12 Lu 0.32 0.1 0.35 0.53 0.47 0.35 Y 22.5 9.68 20.1 32.3 27.1 22.34 ΣREE 167.4 103 101.4 197.69 170.16 147.93 LREE 150.97 96.69 88.70 175.09 151.47 132.58 HREE 22.90 38.84 15.76 16.84 17.63 15.35 LREE/HREE 9.19 15.34 6.98 7.75 8.10 9.47 (La/Yb)N 12.63 25.64 6.81 7.61 8.84 12.31 δEu 0.83 0.86 0.99 0.74 0.77 0.84 (La/Sm)N 3.89 6.46 3.29 3.18 3 3.63 (Gd/Yb)N 2.20 2.72 1.36 1.59 1.90 1.81 Sm/Nd 0.19 0.16 0.21 0.19 0.21 0.19 A/CNK 1.43 1.45 1.01 1.25 0.95 1.22 σ43 1.74 1.55 2.09 2.27 1.97 1.92 注:主量元素含量单位为%,微量和稀土元素单位为10-6 表 2 宝兴岩体花岗闪长岩LA-ICP-MS锆石U-Th-Pb同位素分析结果
Table 2 LA-ICP-MS zircon U-Th-Pb isotope analytical data of Baoxing granodiorite rock
点号 Pb/10-6 232Th/10-6 238U/10-6 Th/U 同位素比值 同位素年龄值/Ma 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 208Pb/232Th 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 208Pb/232Th 1σ 1 20.8 330 450 0.73 0.0508 0.002 0.2599 0.011 0.0370 0.001 0.0115 0.001 232 89 235 9 234 3 230 10 2 26.2 459 554 0.83 0.0533 0.003 0.2726 0.016 0.0372 0.001 0.0112 0.001 339 140 245 13 235 4 224 10 3 35.0 520 779 0.67 0.0513 0.002 0.2613 0.008 0.0372 0.001 0.0109 0.001 254 68 236 7 235 3 218 8 4 30.0 819 544 1.50 0.0519 0.002 0.2653 0.009 0.0372 0.001 0.0104 0.001 283 74 239 7 235 3 208 8 5 19.7 336 431 0.78 0.0525 0.002 0.2694 0.010 0.0371 0.001 0.0100 0.001 306 81 242 8 235 3 201 8 6 21.8 224 491 0.46 0.0517 0.002 0.2670 0.009 0.0376 0.001 0.0105 0.001 272 78 240 7 238 3 212 9 7 48.3 1366 879 1.55 0.0514 0.003 0.2649 0.018 0.0373 0.001 0.0104 0.001 257 182 239 14 236 3 209 8 8 14.5 235 314 0.75 0.0532 0.002 0.2718 0.012 0.0374 0.001 0.0108 0.001 339 102 244 10 237 4 217 9 9 33.0 767 658 1.17 0.0518 0.002 0.2657 0.009 0.0374 0.001 0.0104 0.001 276 74 239 7 236 3 210 8 10 22.9 477 467 1.02 0.0541 0.002 0.2720 0.011 0.0366 0.001 0.0103 0.001 376 103 244 9 232 3 208 9 11 24.7 733 528 1.39 0.0530 0.003 0.2683 0.013 0.0369 0.001 0.0058 0.001 328 109 241 10 233 3 116 5 12 21.8 293 493 0.59 0.0521 0.002 0.2654 0.011 0.0372 0.001 0.0106 0.001 287 87 239 9 236 4 212 10 13 32.2 482 689 0.70 0.0528 0.002 0.2698 0.009 0.0372 0.001 0.0104 0.001 317 74 243 7 235 3 208 9 14 29.8 352 644 0.55 0.0519 0.002 0.2671 0.009 0.0377 0.001 0.0106 0.001 280 78 240 7 238 4 213 10 15 30.5 599 607 0.99 0.0515 0.002 0.2619 0.008 0.0370 0.001 0.0105 0.001 261 69 236 7 234 4 210 9 16 21.4 362 433 0.84 0.0526 0.003 0.2669 0.016 0.0369 0.001 0.0103 0.001 309 137 240 13 234 3 208 11 17 30.6 593 613 0.97 0.0515 0.003 0.2570 0.013 0.0363 0.001 0.0099 0.001 261 94 232 11 230 4 200 8 18 22.7 404 447 0.90 0.0522 0.002 0.2655 0.008 0.0372 0.001 0.0104 0.001 295 74 239 7 236 3 210 8 19 25.6 342 544 0.63 0.0513 0.002 0.2617 0.009 0.0371 0.001 0.0105 0.001 254 69 236 7 235 3 212 8 20 28.2 505 551 0.92 0.0528 0.002 0.2687 0.010 0.0370 0.001 0.0107 0.001 320 116 242 8 234 3 216 8 21 37.5 614 758 0.81 0.0508 0.001 0.2567 0.007 0.0366 0.001 0.0105 0.001 232 66 232 6 232 3 211 8 22 33.1 589 621 0.95 0.0513 0.002 0.2627 0.009 0.0374 0.001 0.0107 0.001 254 80 237 7 236 3 215 8 23 26.9 416 526 0.79 0.0518 0.003 0.2692 0.015 0.0374 0.001 0.0104 0.001 280 120 242 12 237 4 209 9 24 16.1 281 304 0.92 0.0529 0.003 0.2700 0.016 0.0371 0.001 0.0110 0.001 324 135 243 13 235 4 221 11 25 7.5 62 157 0.40 0.0515 0.005 0.2619 0.026 0.0370 0.001 0.0125 0.001 265 225 236 21 234 5 252 18 -
张拴宏, 赵越, 刘建民, 等.华北地块北缘晚古生代-早中生代岩浆活动期次、特征及构造背景[J].岩石矿物学杂志, 2010, 29(6):824-842. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yskw201006017&dbname=CJFD&dbcode=CJFQ 肖庆辉, 王涛, 邓晋福, 等.中国典型造山带花岗岩与大陆地壳生长研究[M].北京:地质出版社, 2009:20-30. 高秉璋, 洪大卫.花岗岩类区1:5万区域地质填图方法指南[M].武汉:中国地质大学出版社, 1991:25-40. 吴福元, 李献华.花岗岩成因研究的若干问题[J].岩石学报, 2007, 23(6):1217-1238. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200706000&dbname=CJFD&dbcode=CJFQ 赵建华.关于岩石微量元素构造环境判别图解使用的有关问题[J].大地构造与成矿学, 2007, 31(1):92-103. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dgyk200701012&dbname=CJFD&dbcode=CJFQ 马昌前, 杨坤光.花岗岩类岩浆动力学[M].武汉:中国地质大学出版社, 1994:56-62. Wiedenbeck M, Alle P, Corfu F, et al. Three natural zircon standards for U-Th-Pb, Lu-Hf, trace element and REE analyses[J]. Geostandards and Geoanalytical Research, 1995, 19(1):1-23. doi: 10.1111/ggr.1995.19.issue-1
Ludwig K R. ISOPLOT3. 00:AGeochronological Toolkit for Microsoft Excel[M]. Berkeley Geochronology Center, California, Berkeley, 2003:39.
李惠民, 李怀坤, 陈志宏, 等.基性岩斜锆石U-Pb同位素定年3种方法之比较[J].地球化学, 2007, 26(2):128-135. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dgyk200701012&dbname=CJFD&dbcode=CJFQ 李艳广, 汪双双, 刘民武, 等.斜锆石LA-ICP-MS U-Pb定年方法及应用[J].地质学报, 2015, 89(12):2400-2418. doi: 10.3969/j.issn.0001-5717.2015.12.015 Wu F Y, Yang Y H, Xie L W, et al. Hf isotopic compositions of the standard zircons and baddeleyites used in U-Pb geochronology[J]. Chem. Geol., 2006, 234:105-126. doi: 10.1016/j.chemgeo.2006.05.003
Chu N C, Taylor R N, Chavagnac V, et al. Hf isotope ratio analysis using multi-collector inductively coupled plasma mass spectrometry:an evaluation of isobaric interference corrections[J]. J. Anal. At. Spectrom., 2002, 17, 1567-1574. doi: 10.1039/b206707b
Maniar P D, Piccoli P M.从微量元素方面来对花岗岩构造背景进行判别[J]. Mineralium Deposita, 1990, 25(1):1-18. http://www.oalib.com/paper/4891814 郭春丽, 吴福元.花岗岩定位机制研究进展[J].世界地质, 2002, 21(4):313-320. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sjdz200204000&dbname=CJFD&dbcode=CJFQ 吴福元.壳-幔物质交换的岩浆岩石学研究[J].地学前缘, 1998, 5(3):94-103. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dxqy803.012&dbname=CJFD&dbcode=CJFQ 万天丰.中国大地构造学纲要[M].北京:地质出版社, 2004. 吉林省地质调查院. 辽源市幅(K51C002004)1: 25 万区域地质矿产调查地质图及成果报告. 2004. 辽宁省地质矿产调查院. 辽宁 1: 5万房木镇等四幅区域地质矿产调查地质图及成果报告. 2015. 辽宁省地质勘察院. 辽宁省区域地质志报告. 2013. 辽宁省地质矿产调查院. 建平县幅(K50C003004)1: 25 万区域地质调查地质图及成果报告. 2003. 辽宁省地质勘察院. 辽宁 1: 5 万阜新市(K51E010007)、新丘(K51E010008)幅区域地质调查地质图及成果报告. 2010.