Zircon U-Pb age, Hf isotopic composition of the Bulangshan Ordovician granite in the Menghai area, southwestern Yunnan Province, and its tectonic significance
-
摘要:
对滇西南勐海布朗山花岗岩进行LA-ICP-MS锆石U-Pb年龄测试，获得该花岗岩的206Pb/238U年龄加权平均值为458.5±3.0Ma（n=21，MSWD=2.3），属中-晚奥陶世的产物。样品锆石的εHf（t）变化范围为-2.4~0.9之间，平均值为-0.9；亏损地壳模式年龄tDMC变化范围为1.41~1.58Ga，加权平均值为1.49Ga。岩石地球化学特征表明，该花岗岩具有高SiO2（75.79%~77.56%）、富碱（K2O+Na2O=7.39%~8.42%）、中-高钾（K2O/Na2O=1.23~1.95）、低MgO（0.14%~0.27%）、低CaO（0.05%~0.64%）的特征。铝过饱和指数A/CNK值介于1.02~1.27之间，岩石属高钾钙碱性强过铝S型花岗岩。岩石稀土元素总体表现为由轻稀土元素富集型向平坦的"海鸥"型过渡，具强烈负Eu异常，δEu值为0.03~0.14；在原始地幔标准化微量元素蛛网图上，明显富集Rb、Th、Nd、Ta等元素，具有明显的Ba、Sr、P、Ti亏损。结合区域地质资料分析，该岩体的形成可能与原特提斯洋俯冲，以及古特提斯洋扩张成盆时，该区处于强烈挤压状态，持续的俯冲使澜沧陆壳残片地壳不断加厚并导致区域重力均衡隆升，引起深部地壳物质在加热后抬升减压发生部分熔融有关。
Abstract:The LA-ICP-MS U-Pb dating of the zircon from the Bulangshan granite in the Menghai area of Yunnan Province was carried out. The weighted average age of the granite 206Pb/238U was 458.5±3.0Ma (n=21, MSWD=2.3), suggesting that the Bulangshan granite was not formed in Late Yanshanian period as originally thought, but was produced in Middle-Late Ordovician. The zircon εHf(t) of the samples range from -2.4 to 0.9 with an average of -0.9. The tDMC of the crustal model vary from 1.41Ga to 1.58Ga with a weighted average of 1.49Ga. It is inferred that the original rock was probably the second cloud (quartz) schist. The geochemical characteristics of the rock indicate that the granite has high SiO2 (75.79%~77.56%), high alkali (K2O+Na2O=7.39%~8.42%), medium-high potassium (K2O/Na2O=1.23~1.95), and low MgO (0.14%~0.27%) and CaO(0.05%~0.64%). Aluminum supersaturated index A/CNK values range from 1.02 to 1.27, and the rock is a high-potassium-calcium-alkali-strong over-aluminum S-type granite. The total amount of rare earth elements in the rocks appears to be transition from LREE enrichment to a flat"gull"type, with a strong negative Eu anomaly and δEu values ranging from 0.03 to 0.14. There exists significant enrichment of elements such as Rb, Th, Nd, and Ta and significant loss of Ba, Sr, P, and Ti in the original mantle standard trace element spider pattern. Combined with regional geological data, the authors hold that the formation of the rock mass may be related to the following factors:the subduction of the original Tethys Ocean and the Paleo-Tethys Ocean expanded into basins, the area was under the condition of intense compression, and the continued intracontinental subduction and collision would continuously increase the regional crust and led to regional equilibrium gravity uplift, which caused deep crust materials to migrate upward after heating and decompression and partial melting.
-
Keywords:
- Bulangshan granite /
- Ordovician /
- original Tethys /
- zircon U-Pb age /
- Hf isotope /
- Menghai
-
生态系统服务指生态系统所形成和维持的人类赖以生存的自然环境条件与效用，是人类直接或间接从生态系统得到的所有收益[1]。黄河流域是中国重要的生态屏障，拥有多个国家重点生态功能区，生态系统服务功能十分重要。黄河流域在近40年的经济快速发展的同时，面临水资源短缺、环境污染、自然生态系统面积减少、生物多样性降低、生境破碎、生态灾难频发等系列资源环境问题[2]。从“十八大”报告提出的“把生态文明建设放在突出地位”，到2019年习近平总书记提出的实现黄河流域生态保护与高质量发展，都是要求在发展经济的同时加强生态保护，在发展与保护之间找到平衡点，实现“双赢”。对黄河流域生态系统服务价值变化进行评估，可为流域生态地质调查、生态环境保护、国土空间规划、环境经济核算、生态补偿等决策提供重要依据，也是实现流域高质量发展的关键。
生态系统服务评价开始于20世纪50年代。1951年，美国水资源委员会发表的《流域经济分析的实践建议》从经济学的角度对流域进行定量分析，成为生态系统服务经济价值研究的起始点[3]。之后，通过损害成本法、收益转移、条件价值法、享乐价格和交通成本法、当量因子法、能值分析法等多种方法将生态系统服务的经济价值和一部分社会文化价值币值化[4]。国外已经对密西西比河流域[5]、亚马逊流域[6-7]、湄公河流域[8]进行了生态系统服务的评估，为系统了解流域生态环境问题及产品和服务供给能力提供了基础。中国也针对黄河流域的区域生态系统服务价值估算开展了大量理论与技术方法研究，如牛叔文[9]估算了黄河上游玛曲地区的生态系统服务价值，宋伟伟等[10]计算了兰州地表水生态服务价值，张楠等[11]定量计算了陕西省安塞县1999—2010年的退耕林生态系统服务价值，丁辉等[12]估算了黄河上游甘南段生态系统服务价值，杨荣等[13]对黄河包头湿地生态系统服务价值进行了货币化估算等。但大多数研究集中于流域局部区域或个别生态系统类型中，缺乏对黄河流域整体性、包含各生态系统的生态系统服务价值及其变化的研究。
本次研究结合国内外经验，根据黄河流域实际情况，以黄河流域整体为研究对象，并考虑公众和决策者对生态服务的理解状况，将生态服务划分为供给服务：食物生产、原材料生产、水资源供给，调节服务：气体调节、气候调节、水文调节、净化环境，支持服务：土壤保持、维持养分循环、生物多样性，文化服务：美化景观，四大类11项服务，将生态系统类型划为森林、草地、农田、湿地、水体和荒漠6类。定量分析流域生态系统服务功能的空间变化和价值变化，揭示黄河流域生态系统服务功能的强弱及生态系统的稳定性状况，有助于管理者对黄河流域生态系统所提供服务的类型、大小与相对组合进行权衡管理，可以为流域国土空间结构的优化调整和自然资源综合管理提供决策依据。
1. 研究区概况
黄河干流河道全长5464 km，流经9省71个市(包括州、盟)，本次研究以流经市域范围作为流域范围，流域总面积198.46×104 km2。流域内土地、矿产、水资源等自然资源非常丰富，具有极大的发展潜力，在中国经济社会发展和生态环境安全方面具有十分重要的意义。从黄河流域2015年生态系统服务价值分布图(图 1)可知，黄河流域11种生态系统服务价值总量在2015年达到5.03亿元，占当年流域GDP总量的43%；流域生态系统服务以水文、气候、土壤保持调节服务为主，三者占总生态服务价值比重分别为25.7%、20.0%、18.0%；流域生态系统服务价值空间分布上呈“南高北低、上下游低、中游高”的空间特征，中游地区大于下游地区，下游地区大于上游地区。其中，生态系统总服务价值最高的区域为甘肃省南部、阿坝藏族羌族自治州和陕西省，每公顷超过3万元，部分地区大于8万元。
![]()
图 1 2015年黄河流经市域生态系统服务价值分布图Figure 1. Distribution of ecosystem service value of the city through which the Yellow River was flowing in 20152. 数据来源与方法
研究数据来自中国科学院资源环境科学数据中心提供的2000年、2005年、2010年、2015年4期中国陆地生态系统服务价值空间分布数据集、基于DEM提取的中国流域、河网数据集等[14-15]。其中，陆地生态系统服务价值空间分布数据是以中国陆地生态系统类型遥感分类结果为基础，根据谢高地等[16-17]的生态服务价值当量因子法，借助全国净初级生产力、土壤保持空间分布数据、降水量，分别对生态系统内各服务价值当量因子价值进行一定的调整，计算得到全国范围内食物生产、水资源供给、水文调节等11种生态服务价值，单种生态服务价值求和得到总服务价值。研究还使用了北京大学城市与环境学院地理数据共享服务平台(http://geodata.pku.edu.cn)提供的2015年中国地市行政边界数据。
借助ArcGIS空间分析功能，按黄河流经地市边界范围重新提取生成黄河流域边界数据，从4期中国陆地生态系统服务价值数据中提取生成黄河流域生态系统服务价值数据。运用Google Earth Engine平台提供的统计分析函数工具，采用线性回归[18-19]、滑动平均[20-21]和Sen’s斜率估计[22]等趋势分析方法，评估黄河流域四大类共11种生态系统服务价值的变化速率，最终得到1 km×1 km分辨率流域每公顷十一种生态服务价值，以及总生态服务价值年变化速率的空间分布图。
3. 单项服务价值变化评估
(1) 食物生产和原材料生产价值变化
2000—2015年，流域食物生产价值和原料生产价值均呈递增趋势(图 2)，玉树—阿坝一带山区和下游引黄灌区减少明显。其中，食物生产价值由1027亿元上涨到1727亿元，上涨了68%；原料生产价值由990亿元上涨到1554亿元，上涨了57%。
![]()
图 2 2000—2015年原料和食物生产价值变化Figure 2. Changes in the value of raw materials and food production during 2000—2015宁夏河套平原、黄河三角洲地区农田面积呈增加趋势，且宁蒙河套平原、汾渭盆地农田生产潜力增长趋势较大，食物生产价值增加。下游引黄灌区的河南和山东两省，农田总面积显著下降，食物生产价值呈减少趋势。其余地区变化速率较低。
(2) 气体和气候调节价值变化
黄河流域的不同生态系统通过其系统内部水分的蒸发、蒸腾、光合作用等一定的调温和减弱温室效应功能，其中森林和湿地是气体和气候调节功能的主要承担者。2000—2015年，黄河流域气体调节价值总体上呈递增的趋势(图 3)，由2974亿元上涨到4690亿元，增长幅度为58%，其中在2005—2010年间增速放缓。2000年黄河流域气候调节价值为7223亿元，到2015年为止为10889亿元，增长近51%。气候调节价值变化和气体调节变化趋势大致相同，增长速率高的地区主要为黄土丘陵沟壑区、榆林地区、陇中黄土高原区、汾河盆地等地区，主要是由于该地区植被(森林和草地)逐渐恢复。
![]()
图 3 2000—2015年气体和气候调节价值变化Figure 3. Change in the value of gas and climate regulation during 2000—2015上游西部玉树地区作为全球气候变化的敏感区，受温室气体排放加剧影响比流域其他地区更显著，增温幅度最大，调节价值变小。上游非保护区域过度放牧、中游地区工农业用水激增，导致河流下泄量日趋减少，使下游河床、湖泊干涸，地下水位下降，湿地、草地生态系统退化明显，导致气体和气候价值减少。
(3) 水资源供给价值变化
提供淡水是流域生态系统的服务功能之一。流域资源型缺水禀赋条件长期存在，不考虑冰川融水的生态系统服务流的变化，受整体气候暖干化影响，流域水资源供给价值整体减少(图 4)。在2000—2005的5年间由818亿元上升22.9%，达到1005亿元，但是在2005年后持续下降，下降20.6%，降至798亿元。
![]()
图 4 2000—2015年水资源供给和水文调节价值变化Figure 4. Change of water resources supply and hydrological regulation value from 2000 to 2015黄河上游青铜峡以上流域受气候暖干化影响，径流量和水资源供给价值呈减少趋势，主要是汛期降水减少，非汛期降水增加，非汛期产流系数低于汛期；大范围长历时有利于产流的降水过程减少，降水日数减少；降水强度有所增加，局地暴雨增多，以及增温影响。局部海西地区因为气候暖湿化影响，降水量增多，径流量增加。中部宁夏平原、河套平原、黄土高原区、泾河流域、渭河流域，以及下游黄河三角洲等地，由于植被恢复，水土流失得到有效控制，降雨量增加，径流量呈增加趋势，水资源供给价值呈增加趋势。
(4) 水文调节价值变化
生态水文调节可以理解为生态系统通过水库、湖泊、沼泽等对降水截留、过滤、吸收作用等改变降雨径流的时空分配，在一定程度上起到削峰补枯、缓和地表径流、增加地下径流等的作用。黄河流域水文调节价值在2000—2015年间呈先增后减的趋势(图 4)，2005年为转折点。其中，在2000—2005年间，黄河流域水文调节价值由1.07万亿增长到1.29万亿元，增长了20%。但是2005年后，水文调节价值不断降低，降到1.06万亿元，下滑18%。
黄河流域水少沙多，水沙关系不协调。除中游鄂尔多斯高原、汾渭盆地，上游海西、海北及阿拉善盟西部暖湿化地区，下游黄河三角洲地区增加外，大部分地区减少。流域水文调节价值的减少除与大尺度的气候变暖有关外，区域尺度的植被覆盖变化和人类用水量(如农业灌溉)的增加也是不可忽视的原因。在黄河上游，气候变化的影响占主导地位；在中下游，人类活动的影响甚至与气候变化作用旗鼓相当。
(5) 生物多样性调节价值变化
生物多样性对于保持水土、调节气候、维持生态平衡、稳定环境具有关键性作用。黄河流域生物多样性调节价值在2000—2015年呈递增趋势(图 5)，由2000年的3012亿元上涨到4775亿元，增幅59%，其中在2005—2010年间变化较平缓。
![]()
图 5 2000—2015年生物多样性调节和净化环境价值变化Figure 5. Changes in the environmental value of biodiversity regulation and purification from 2000 to 2015宁蒙河套平原、黄土高原、西北暖湿化地区等大部分地区生物多样性价值呈增长趋势。而上游生态环境脆弱区、下游城市化扩张严重地区，存在河流空间过度管控、生态空间被挤占、水资源过度开发等问题，导致一些重要生态空间与干支流重要廊道的生态功能遭到破坏。
(6) 净化环境价值变化
净化服务是生态系统重要的生态服务之一，通过生态系统的生态过程，在一系列物理、化学和生物作用下，将废弃物降解和净化。2000—2015年黄河流域的净化环境价值呈递增趋势(图 5)，由2396亿元增加到4099亿元，增加了1703亿元，占2000年价值总量的71%。
山区生态脆弱区、关中-天水经济区、下游中段等地区存在农业面源污染、河流底泥污染问题，部分河段纳污超载问题严重，净化环境价值呈减少态势。流域大部分地区净化环境价值呈增加趋势，流域内生态环境正在得到积极有效的恢复。
(7) 土壤保持价值变化
2000—2015年呈现“增-减-增”的变化趋势(图 6)，2000年初总价值为6810亿元，到2005年为止增加到8916亿元；之后到2010年减少到8088亿元。2010年后，土壤保持价值缓慢上升到8498亿元，整体上呈增加的趋势，增长了19%。土壤保持价值减少地区主要位于晋陕峡谷段和汾河地堑地区。黄河晋陕峡谷段黄土侵蚀严重，是主要的产沙段，尤其是西部黄土高原，黄土高原沟谷发育与黄河阶地相对应[23]，下游阶地级数多，黄土侵蚀破碎严重，沟谷长[24]；上游阶地级数少，黄土较为完整，沟谷较短。汾河地堑地区第四纪早中期发育的河流也对盆地边缘的黄土造成强烈侵蚀，黄土多支离破碎，很少残留塬面，沟谷多为老年U型沟谷[25-28]。流域整体上土壤保持价值呈增加趋势，和流域内荒漠的减少、植被覆盖的增加有密不可分的关系[27]。从2000年开始流域荒漠化扩展得到遏制，流域中游黄土高原蓄水保土能力显著增强，实现了“人进沙退”的治沙奇迹[28]。
![]()
图 6 2000—2015年土壤保持和维持养分价值变化Figure 6. Changes in the value of soil conservation and maintenance nutrients from 2000 to 2015(8) 维持养分循环价值变化
养分元素的循环利用是黄河流域内生态系统主要的功能之一。2000—2015年，流域内维持养分循环价值呈递增趋势(图 6)，2000年价值总量为324亿元，15年内上涨了57%，达到509亿元。下游沿黄农业生产基地、关中—天水经济区、玉树—果洛—阿坝一带山区等地维持养分循环价值呈减少趋势。宁夏平原、内蒙河套平原、汾渭盆地、甘肃南部，以及延安黄土丘陵沟壑区等地区维持养分循环价值显著增加[29]。
(9) 景观美学价值变化
景观美学是景观服务的重要组成部分。2000—2015年，黄河流域景观美学价值总体上呈递增趋势(图 7)，由1353亿元跃为2206亿元，增长幅度达63%。中上游和西北部地区随着植被恢复、荒漠化遏制、水多沙少等情况的改善，景观美学价值不断增加。但玉树—阿坝—商洛一带林木分布转向由单一树种或生态特性相似的树种组成的单纯林为主趋势，景观美学价值降低。近30年流域内湖泊湿地、沼泽湿地、河口滩涂湿地等重要自然湿地分别减少25%、21%和40%，而景观单一的人工湿地面积增加60%以上。黄河下游滩区居住有190万居民，开发建设活动密集，存在三角洲侵蚀衰退、河流廊道破碎化严重现象[30-32]。
![]()
图 7 2000—2015年美学景观价值变化Figure 7. The change of aesthetic landscape value from 2000 to 20154. 总服务价值变化评估
(1) 总生态系统服务价值除水文调节和水资源供给服务呈下降趋势外，总体呈上升趋势。流域总服务价值呈上升趋势。除水文调节和水资源供给服务先增后减，呈下降趋势外，其他9项服务均呈增加趋势。通过对流域食物生产、原材料生产、水资源供给等11项生态服务价值进行求和，得到流域总生态服务价值。2000年流域总生态服务价值为3.77万亿元，2005年上涨22%，为4.59万亿元，2010年下降到4.48万亿元，到2015年又上涨至5.03亿元，相比2000年增长了33.4%。
(2) 总生态系统服务价值呈现上游南部山区明显减少、中游黄土高原区显著增加、下游轻微降低的空间格局特点(图 8)。
![]()
图 8 2000—2015年总生态服务价值变化速率空间分布Figure 8. Spatial distribution of the change rate of total ecological service value from 2000 to 2015① 上游属于青藏高原和西北部干旱风沙区，物质变化和能量转换缓慢，抗干扰水平和自我平衡水平较差，一旦遭到破坏将很难逆转，生态环境脆弱性高，受全球气候变化影响，玉树—果洛—阿坝一带山区以及青海湖服务价值呈明显下降趋势；②中游黄土高原区生态服务价值增加尤为明显，增幅达48%，其南部改善面积较大而北部地区改善速度缓慢，价值增加主要是因为近年生态保护与修复取得积极有效的成果，森林植被覆盖率大幅度提高，荒漠化趋势得到遏制，生物多样性逐步提高，以国家公园为主题的自然保护地体系逐步完善；③中游太行山区和下游引黄灌区地区等城镇化率高的地区资源开发、工程建设等人类活动易于破坏原有的自然环境，导致其生态功能受损，服务价值降低。
(3) 价值总量增长了33.4%，但增速低于GDP增速，生态环境与经济发展处于低度协调状态。流域生态系统服务价值增速低于经济发展增速。2000年、2005年、2010年、2015年GDP分别为14555亿元、52249亿元、66089亿元、116236亿元，生态服务价值占GDP比例由315%、85%、76%，下降至43%，生态系统服务价值和GDP变化率比值为0.048，生态环境与经济发展处于低度协调状态。
(4) 服务供给与下游人口聚居的服务需求之间存在空间不匹配问题。服务供给与需求之间存在空间不匹配问题。如黄河流域生态系统服务分布呈南高北低、上下游低、中游高的空间特征，而上下游由于生态退化其生态系统服务价值在持续下降。下游河南、山东人口密度高、数量大，生态系统服务供给与人口聚居的服务需求之间存在空间不匹配(图 8)。
5. 基于生态系统服务价值评估的生态地质调查工作建议
为恢复、维持进而提高黄河流域生态系统服务功能，针对黄河流域生态系统服务价值变化，从生态地质工作角度出发，提出如下建议。
(1) 动态评估流域生态系统承载力。根据生态系统服务价值变化格局，动态评估黄河流域生态系统承载力，分析生态系统服务在区域之间的流动，积极调整和改善国土空间利用类型结构、地质调查工作部署，实施生态恢复措施，为保护生态系统服务和增强流域可持续性提供有用的信息。
(2) 聚焦关键过程、关键区域、关键要素，开展生态地质调查工作，提升生态系统服务价值，优化流域国土空间结构。其中，在食物生产和原材料生产价值减少明显区，开展土地资源及其利用问题调查研究，开展清洁能源和关键矿产绿色勘查开发；在水资源供给价值减少区，开展“气候-水-人类活动”相互作用机理研究，开展水资源系统行星边界研究，优化配置流域水资源；在气体和气候调节价值下降区，开展气候变化应对研究，分区实施植被恢复措施；在水文调节价值下降区，开展气候变化对水源涵养功能退化影响机制研究，开展重大水利工程环境效应调查研究；在生物多样性调节价值减少明显区，在典型地区建立生态修复示范，开展湿地生态系统修复；在净化环境价值呈减少态势区，开展水体富营养化、土壤污染问题和矿山地质环境问题调查研究；在土壤保持价值明显减少区，开展黄土斜坡水文-应力响应机制研究，查明淤地坝溃决风险；在维持养分循环价值减少区，开展流域生态系统承载力研究；在景观美学价值下降明显区，开展生态林业建设，开展地质遗迹调查、旅游廊道和国家公园建设。
6. 结论
(1) 黄河流域生态系统服务价值以水文调节、气候调节、土壤保持调节服务为主，空间分布上呈“南高北低、上下游低、中游高”的空间特征。
(2) 2000—2015年黄河流域生态系统11种类型服务中，水文调节和水资源供给服务以2005年为转折点，呈先增后减整体小幅下降的趋势，分别减少1%、2.3%；其他9项服务呈波动或递增上升趋势，除土壤保持价值小幅增长19%外，其他均大幅上涨50%以上，净化环境价值涨幅高达71%。
(3) 相比2000年，2015年黄河流域生态系统总服务价值增长了33.4%，呈现“上游玉树—阿坝一带山区明显减少、中游黄土高原区显著增加、下游轻微降低”的空间格局特点；服务价值增速远低于GDP增速，生态环境与经济发展处于低度协调状态；下游生态系统服务与需求之间存在空间不匹配问题。
(4) 从生态地质调查工作角度出发，建议动态评估黄河流域生态系统承载能力，聚焦关键过程、关键区域、关键要素，开展黄河流域生态地质调查工作，提升生态系统服务价值，优化国土空间结构。
致谢: 云南省地质调查院李静教授级高工认真审阅了本文并提出了详细的修改意见，LA-ICPMS锆石U-Pb同位素测试得到武汉上普分析测试有限公司陈红芳老师的支持和帮助，审稿专家提出了宝贵的修改意见，在此一并表示衷心的感谢。 -
![]()
Qh—第四系；N1s—三营组；Em—勐腊组；J2h1—花开左组一段; J2h2—花开左组二段；Pz1dx.—大新山岩组；Pt2h.—惠民岩组；Pt2ml.—曼来 岩组；Pt2m.—勐井山岩组；Pt2n—南木林岩组；Pt1D.—大勐龙岩群；ηοE—古近纪石英二长岩；δηοE—古近纪石英二长闪长岩；βE—古近 纪玄武岩；γO—奥陶纪花岗岩；γκT—三叠纪花岗细晶岩；γT—三叠纪花岗岩；ηT—三叠纪二长岩；ητT—三叠纪花岗细晶岩；ηγT—三 叠纪二长花岗岩；ηγcT—三叠纪中粗粒黑云二长花岗岩；γδP—二叠纪花岗闪长岩；δμC-P—石炭纪-二叠纪闪长玢岩；νC-P—石 炭纪-二叠纪辉长岩；βμC-P—石炭纪-二叠纪辉绿岩；δC—石炭纪闪长岩；νC—石炭纪辉长岩；βμC—石炭纪辉绿岩
Figure 1. Geological map of the study area
![]()
图 2 布朗山白云钾长花岗岩野外及镜下特征
Pl—斜长石；Qtz—石英；Or—正长石；Ms—白云母；Ep—绿帘石；Zrn—锆石
Figure 2. Photograph and microphotograph of muscovite potash feldspar in Bulangshan
![]()
图 3 布朗山花岗岩锆石阴极发光（CL）图像及打点位置
（实线圆圈代表U-Pb分析点，虚线圆圈代表Hf分析点）
Figure 3. CL images and analytical positions of zircons in Bulangshan granites
![]()
图 4 布朗山花岗岩锆石U-Pb谐和图和年龄加权平均值
Figure 4. U-Pb concordia and weighted average diagrams for zircons of Bulangshan granites
![]()
图 7 布朗山花岗岩B-A图解[32]
指数A反映花岗岩类的含铝特征，B反映花岗岩浆的演化程度，单位为每100g物质中每种元素的g-原子数×10-3）
Figure 7. B-A diagram of Bulangshan granites
表 1 布朗山花岗岩LA-ICP-MS锆石U-Th-Pb同位素数据
Table 1 LA-ICP-MS zircon U-Th-Pb data of Bulangshan granites
点号 Pb Th U Th/U 同位素比值 同位素年龄/Ma 10-6 207Pb/206Pb 1σ 207Pb/235U 1σ Pb/U206238 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 01 43.2 327 484 0.68 0.0583 0.0015 0.6059 0.0159 0.0748 0.0007 543 53.69 481 10.08 465 4.50 02 24.4 315 251 1.25 0.0567 0.0016 0.5728 0.0166 0.0729 0.0007 480 64.81 460 10.69 454 4.10 03 21.9 186 248 0.75 0.0596 0.0019 0.6004 0.0180 0.0729 0.0007 591 68.51 477 11.46 454 4.08 04 17.2 143 195 0.73 0.0570 0.0017 0.5819 0.0184 0.0736 0.0009 500 66.66 466 11.81 458 5.45 05 89.7 654 1023 0.64 0.0556 0.0013 0.5605 0.0130 0.0728 0.0006 435 49.07 452 8.49 453 3.80 06 21.1 201 231 0.87 0.0552 0.0018 0.5636 0.0183 0.0739 0.0007 420 78.69 454 11.86 459 4.14 07 14.0 129 155 0.83 0.0550 0.0019 0.5525 0.0194 0.0726 0.0008 413 77.77 447 12.66 452 4.62 08 23.4 237 254 0.93 0.0593 0.0017 0.5969 0.0173 0.0727 0.0007 589 61.1 475 11.02 452 4.10 09 19.5 135 206 0.66 0.0561 0.0017 0.6109 0.0179 0.0790 0.0009 457 66.66 484 11.27 490 5.11 10 21.9 182 269 0.68 0.0578 0.0019 0.5505 0.0190 0.0689 0.0010 524 70.36 445 12.42 429 6.28 11 30.4 279 328 0.85 0.0587 0.0016 0.6062 0.0164 0.0747 0.0007 567 59.25 481 10.39 464 4.11 12 31.57 240 355 0.68 0.0566 0.0015 0.5871 0.0159 0.0748 0.0008 476 89.80 469 10.19 465 4.73 13 88.4 840 981 0.86 0.0566 0.0012 0.5683 0.0117 0.0726 0.0006 476 46.29 457 7.595 452 3.50 14 20.6 228 215 1.06 0.0584 0.0017 0.5956 0.0169 0.0740 0.0008 543 66.65 474 10.78 460 4.74 15 45.8 378 502 0.75 0.0576 0.0016 0.6072 0.0169 0.0764 0.0009 522 62.96 482 10.70 475 5.62 16 20.8 163 230 0.71 0.0580 0.0019 0.5993 0.0199 0.0747 0.0007 528 72.21 477 12.67 464 4.10 17 15.46 131 182 0.72 0.0601 0.0023 0.5804 0.0209 0.0700 0.0007 609 81.47 465 13.41 436 4.44 18 25.09 269 267 1.00 0.0591 0.0023 0.6162 0.0243 0.0751 0.0010 572 83.32 487 15.30 467 5.71 19 35.61 249 394 0.63 0.0555 0.0018 0.5956 0.0191 0.0774 0.0008 432 70.36 474 12.18 480 4.79 20 40.87 234 438 0.53 0.0571 0.0016 0.6525 0.0194 0.0825 0.0010 494 67.59 510 11.94 511 5.93 21 33.49 255 373 0.68 0.0585 0.0018 0.6087 0.0182 0.0752 0.0007 550 64.80 483 11.48 467 4.39 22 20.70 186 233 0.80 0.0565 0.0018 0.5655 0.0171 0.0725 0.0006 472 70.36 455 11.12 451 3.86 23 18.65 190 198 0.96 0.0575 0.0019 0.5925 0.0192 0.0745 0.0008 509 76.84 472 12.23 463 4.56 24 16.50 152 176 0.86 0.0634 0.0021 0.6567 0.0217 0.0748 0.0007 724 70.36 513 13.32 465 4.07 25 39.9 467 358 1.30 0.0638 0.0019 0.7196 0.0210 0.0819 0.0010 744 62.19 550 12.40 507 6.21 26 44.9 415 442 0.94 0.0571 0.0016 0.6293 0.0171 0.0795 0.0006 494 59.26 496 10.64 493 3.84 27 21.40 172 245 0.70 0.0578 0.0019 0.5804 0.0184 0.0728 0.0007 520 70.36 465 11.82 453 4.09 28 30.9 139 174 0.80 0.1889 0.0147 3.3627 0.4494 0.1031 0.0041 2733 129.5 1496 104.6 633 23.77 29 67.7 387 733 0.53 0.0564 0.0012 0.6348 0.0143 0.0811 0.0008 478 50.92 499 8.86 503 4.79 30 20.16 187 221 0.85 0.0590 0.0019 0.5965 0.0195 0.0730 0.0007 565 70.36 475 12.39 454 4.48 
下载: 导出CSV
表 2 布朗山花岗岩全岩主量、微量和稀土元素测试数据
Table 2 Whole-rock major, trace and rare earth element analyses of the Bulangshan granites
样品号 MBbls-32-1 MBbls-32-1-1 MBbls-32-1-2 MBbls-32-1-3 MBxng-36-1 MBxng-36-1-1 MBxng-36-1-2 样品号 MBbls-32-1 MBbls-32-1-1 MBbls-32-1-2 MBbls-32-1-3 MBxng-36-1 MBxng-36-1-1 MBxng-36-1-2 SiO2 77.54 76.76 77.56 76.61 77.02 77.25 75.79 Sc 1.82 1.83 1.56 2.78 2.02 1.89 2.23 TiO2 0.04 0.04 0.04 0.04 0.09 0.08 0.09 Nb 60.20 45.10 47.60 60.20 28.90 26.80 40.00 Al2O3 12.33 12.60 12.11 12.62 12.17 12.05 12.41 Zr 135.00 55.90 64.80 59.80 29.50 35.30 38.20 Fe2O3 0.84 0.92 0.71 0.89 0.96 1.02 1.04 Hf 7.14 3.23 3.44 3.38 1.54 1.59 2.09 FeO 0.70 0.65 0.59 0.68 0.78 0.83 0.89 U 10.00 4.30 6.67 4.46 5.47 4.11 6.98 MnO 0.02 0.01 0.02 0.02 0.03 0.03 0.03 Th 38.70 37.30 32.90 41.9 25.50 28.20 21.30 MgO 0.22 0.27 0.25 0.14 0.26 0.25 0.22 La 14.60 20.10 32.10 15.50 20.60 31.20 29.10 CaO 0.10 0.05 0.11 0.07 0.01 0.11 0.64 Ce 21.90 26.80 45.20 12.50 35.20 59.70 55.70 Na2O 2.54 3.64 2.56 3.31 3.57 3.39 3.65 Pr 4.28 6.41 9.25 5.27 4.69 6.84 6.50 K2O 4.95 4.47 4.83 4.58 4.85 4.58 4.67 Nd 16.9 27.00 38.10 22.90 16.70 24.90 24.20 P2O5 0.01 0.01 0.01 0.01 0.02 0.02 0.02 Sm 4.52 8.46 10.6 7.39 3.43 5.50 5.85 烧失量 0.91 0.72 1.28 1.18 0.44 0.7 0.94 Eu 0.06 0.08 0.13 0.07 0.15 0.23 0.23 总计 100.20 100.14 100.07 100.15 100.29 100.31 100.39 Gd 4.66 8.29 10.10 7.52 2.99 5.31 5.56 A/NK 1.29 1.16 1.28 1.21 1.09 1.14 1.12 Tb 1.13 1.90 2.13 1.84 0.60 1.10 1.25 A/CNK 1.27 1.15 1.26 1.20 1.08 1.12 1.02 Dy 7.79 11.60 13.10 11.60 3.41 6.43 7.69 Na2O+K2O 7.49 8.11 7.39 7.89 8.42 7.97 8.32 Ho 1.66 2.39 2.72 2.47 0.67 1.33 1.67 K2O/ Na2O 1.95 1.23 1.89 1.38 1.36 1.35 1.28 Er 4.77 6.53 7.59 6.86 1.84 3.78 4.82 CaO/ Na2O 0.04 0.01 0.04 0.02 0.03 0.03 0.18 Tm 0.88 1.18 1.35 1.23 0.33 0.66 0.92 Cu 2.72 1.89 1.32 2.02 1.87 2.42 8.93 Yb 5.71 7.18 8.29 7.93 2.10 4.03 6.06 Zn 41.20 40.90 35.50 48.20 31.10 36.00 26.50 Lu 0.80 0.99 1.15 1.03 0.28 0.53 0.84 Cr 2.07 4.10 7.32 5.47 2.79 2.58 122.00 Y 42.10 60.70 73.50 62.80 15.80 36.60 48.30 Ni 0.002 0.002 0.30 0.05 0.17 0.69 48.80 ΣREE 131.76 189.61 255.31 166.91 108.80 188.14 198.68 Co 0.26 0.16 0.16 0.23 0.63 0.66 2.18 δEu 0.04 0.03 0.04 0.03 0.14 0.13 0.12 Rb 322 264 293 339 202 197 223.00 (La/Yb)N 1.72 1.89 2.61 1.32 6.61 5.22 3.24 Sr 7.72 4.38 8.13 3.25 9.97 7.37 13.00 (Gd/Yb)N 0.66 0.93 0.98 0.77 1.15 1.06 0.74 Ba 33.70 25.30 51.30 26.40 164.00 158.00 173.00 (La /Sm)N 2.03 1.49 1.90 1.32 3.78 3.57 3.13 V 4.23 3.37 2.80 3.15 4.61 4.73 5.10 注：主量元素含量单位为%，微量和稀土元素含量单位为10-6；A/NK=摩尔Al2O3/(Na2O+K2O)，A/CNK=摩尔Al2O3/(CaO+Na2O+K2O)；δEu=2EuN/(SmN+GdN)，其中N为球粒陨石标准化值
下载: 导出CSV
表 3 布朗山花岗岩锆石Hf同位素数据
Table 3 Zircon Hf isotopic compositions of the Bulangshan granites
测点号 年龄/Ma 176Hf/177Hf 2σ 176Yb/177Hf 2σ 176Lu/177Hf 2σ (176Hf/177Hf)t εHf(1t) tDM/Ma tDMC/Ma fLu/Hf 01 465 0.282460 0.000008 0.036197 0.000196 0.001428 0.000009 0.282448 -1.4 1132 1519 -0.96 02 454 0.282476 0.000008 0.029577 0.000224 0.001135 0.000005 0.282466 -0.7 1101 1478 -0.97 03 454 0.282429 0.000009 0.020115 0.000126 0.000786 0.000003 0.282422 -2.3 1157 1576 -0.98 04 458 0.282495 0.000008 0.019454 0.000128 0.000769 0.000004 0.282488 0.1 1064 1429 -0.98 05 453 0.282460 0.000009 0.031819 0.000139 0.001241 0.000006 0.282449 -1.3 1127 1516 -0.96 06 459 0.282427 0.000008 0.02309 0.000168 0.000880 0.000005 0.282419 -2.4 1163 1583 -0.97 07 452 0.282468 0.000009 0.02575 0.000232 0.000994 0.000010 0.282459 -1.0 1109 1494 -0.97 08 452 0.282457 0.000008 0.019423 0.000847 0.000737 0.000010 0.282451 -1.3 1116 1513 -0.98 11 464 0.282462 0.000008 0.021789 0.000618 0.000848 0.000004 0.282455 -1.1 1112 1504 -0.97 12 465 0.282504 0.000010 0.021653 0.000115 0.000875 0.000003 0.282496 0.3 1055 1412 -0.97 13 452 0.282472 0.000007 0.019145 0.000831 0.000767 0.000003 0.282466 -0.7 1096 1480 -0.98 14 460 0.282469 0.000008 0.023286 0.000913 0.000914 0.000003 0.282461 -0.9 1104 1489 -0.97 15 475 0.282526 0.000008 0.039758 0.000499 0.001592 0.000017 0.282512 0.9 1044 1377 -0.95 16 464 0.282468 0.000009 0.021784 0.000239 0.000867 0.000008 0.282461 -0.9 1104 1490 -0.97 18 467 0.282478 0.000008 0.023974 0.000177 0.000915 0.000007 0.282470 -0.6 1093 1471 -0.97 21 467 0.282470 0.000008 0.030781 0.000822 0.001181 0.000035 0.282460 -0.9 1111 1492 -0.96 22 451 0.282461 0.000008 0.030479 0.000491 0.001167 0.000020 0.282451 -1.3 1124 1513 -0.96 23 463 0.282492 0.000011 0.03716 0.000416 0.001412 0.000018 0.282480 -0.2 1087 1448 -0.96 24 465 0.282437 0.000007 0.023064 0.000447 0.000872 0.000002 0.282430 -2.0 1148 1560 -0.97 27 465 0.282495 0.000009 0.021947 0.000872 0.000877 0.000003 0.282488 0.0 1067 1431 -0.97 注：εHf（t）=10000×{[（176Hf/177Hf）s -（176Lu/177Hf）s×（eλt-1）]/[ （176Hf/177Hf）CHUR,0-（176Lu/177Hf）CHUR×（eλt-1）]-1},tDM=1/λ×ln{1+ [（176Hf/177Hf）s-（176Hf/177Hf）DM]/[（176Lu/177Hf）s-（176Lu/177Hf）DM]}, tDMC=tDM-(tDM-t)×[(fcc-fs)/(fcc-fDM)], fLu/Hf=（176Lu/177Hf）s/（176Lu/ 177Hf）CHUR-1.其中：λ=1.867×10-11/a;（176Hf/177Hf）s 和（176Hf/177Hf）CHUR=0.0332,（176Hf/177Hf）CHUR,0=0.282772；（176Lu/177Hf）DM= 0.0384，（176Hf/177Hf）DM=0.28325；（176Lu/177Hf）平均地壳=0.015；fcc=[（176Lu/177Hf）平均地壳/（176Lu/177Hf）CHUR]-1；fs=fLu/Hf；fDM=[（176Lu/177Hf）DM/ （176Lu/177Hf）CHUR]-1; t为锆石结晶年龄
下载: 导出CSV
-
莫宣学, 沈上越, 朱勤文, 等.三江中南段火山岩-蛇绿岩与成矿[M].北京:地质出版社, 1998:123-132. Sone M, Metcalfe I. Parallel Tethyan sutures in mainland Southeast Asia:New insights for Paleo-Tethys closure and implications for the Indosinian orogeny[J]. Comptes Rendus eosciences, 2008, 340(2/3):166-179. http://cn.bing.com/academic/profile?id=5b3488b9c160f6376e1fe1baa10c93f1&encoded=0&v=paper_preview&mkt=zh-cn
彭头平, 王岳军, 范蔚名, 等.澜沧江南段早中生代酸性火成岩SHRIMP锆石U-Pb定年及构造意义[J].中国科学(D辑), 2006, 36(2):123-132. http://d.old.wanfangdata.com.cn/Periodical/zgkx-cd200602002 廖世勇, 尹福光, 王冬兵, 等.滇西"三江"地区临沧花岗岩基中三叠世碱长花岗岩的发现及其意义[J].岩石矿物学杂志, 2013, 33(1):1-12. doi: 10.3969/j.issn.1000-6524.2013.01.001 刘德利, 刘继顺, 张彩华, 等.滇西南澜沧江结合带北段云县花岗岩的地质特征及形成环境[J].岩石矿物学杂志, 2008, 27(1):23-31. doi: 10.3969/j.issn.1000-6524.2008.01.003 秦元季.滇西临沧花岗岩基的基本特征和构造侵位机制[D].中国科学院地质与地球物理研究所博士学位论文, 1991. 李兴林.临沧复式花岗岩基的基本特征及形成环境的研究[J].云南地质, 1996, 15(1):1-18. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600616589 陈吉琛.滇西花岗岩类形成的构造环境及岩石特征[J].云南地质, 1989, 8(3/4):205-212. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000003601933 刘昌实, 朱金初, 徐夕生.滇西临沧复式岩基特征研究[J].云南地质, 1989, 8(3/4):189-204. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000003601932 吕伯西, 王增, 张德能, 等.三江地区花岗岩类及其成矿专属性[M].成都:四川科技出版社, 1993:1-133. 王舫, 刘福来, 刘平华, 等.澜沧江南段临沧花岗岩的锆石U-Pb年龄及构造意义[J].岩石学报, 2014, 30(10):3034-3050. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201410018 Hennig D, Lehmanneh B D, Frei D, et al. Early Permian seafloor to continental arc magmatism in the eastern Paleo-Tethys:U-Pb age and Nd-Sr isotope data from the southern Lancangjiang zone, Yunnan China[J]. Lithos, 2009, 113(3):408-422. http://cn.bing.com/academic/profile?id=f50d8c6d10296b28611653136451672d&encoded=0&v=paper_preview&mkt=zh-cn
孔会磊, 董国臣, 莫宣学, 等.滇西三江地区临沧花岗岩的岩石成因——地球化学、锆石U-Pb年代学及Hf同位素约束[J].岩石学报, 2012, 28(5):1438-1452. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201205008 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
Zong K Q, Klemd R, Yuan Y, et al. The assembly of Rodinia:The correlation of early Neoproterozoic (ca. 900Ma) high-grade metamorphism and continental arc formation in the southern Beishan Orogen, southern Central Asian Orogenic Belt (CAOB)[J]. Precambrian Research, 2017, 290:32-48. doi: 10.1016/j.precamres.2016.12.010
Hu Z C, Zhang W, Liu Y S, et al."Wave" signal smoothing and mercury removing device for laser ablation quadrupole and multiple collector ICP-MS analysis:application to lead isotope analysis[J]. Analytical Chemistry, 2015, 87:1152-1157. doi: 10.1021/ac503749k
Liu Y S, Gao S, Hu Z C, et al. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North 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
Lin J, Liu Y S, Yang Y H, et al. Calibration and correction of LA-ICP-MS and LA-MC-ICP-MS analyses for element contens and isotopic ratios[J]. Solid Earth Science, 2016, 1:5-27. doi: 10.1016/j.sesci.2016.04.002
Ludwig K R. Isoplot 3.00: A Geochronological Toolkit for Microsoft Excel[M]. Berkeley Geochronology Center, California, Berkeley, 2003.
耿建珍, 李怀坤, 张健, 等.锆石Hf同位素组成的LA-ICP-MS测定[J].地质通报, 2011, 30(10):1508-1513. doi: 10.3969/j.issn.1671-2552.2011.10.004 Taylor S R, Mclennan S M. The continental crust:Its composition and evolution[M]. Oxford:Blackell Scientific Publications, 1985:54-372.
刘振声, 王洁民.青藏高原南部花岗岩地质地球化学特征[M].成都:四川科技出版社, 1994:1-133. Alther R, Holl A, Hegner E. High-potassium, calc-alkaline plutonism in the European Variscides:northern Vosges (France) and northern Schwarzwald (Germany)[J]. Lithos, 2000, 50:51-73. doi: 10.1016/S0024-4937(99)00052-3
Zhou X M, Li W X. Origin of Late Mesozoic igneous rocks in South eastern China:Implications for lithosphere subduction and underplating of mafic magmas[J]. Tectonophysics, 2000, 326(3/4):269-287. http://www.sciencedirect.com/science/article/pii/S0040195100001207/pdf?md5=498f70655f1dbfbb3a2e4e7738d90458&pid=1-s2.0-S0040195100001207-main.pdf&_valck=1
吴福元, 李献华, 郑永飞, 等. Lu-Hf同位素体系及其岩石学应用[J].岩石学报, 2007, 23(2):185-220. http://d.old.wanfangdata.com.cn/Periodical/ysxb98200702001 Amelin Y, Lee D C, Halliday A N. Early-middle Archaean crustal evolution deduced from Lu-Hf and U-Pb isotopic studies of single zircon grains[J]. Geochimical et Cosmochimica Acta, 2000, (64):4205-4225. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9483bc569d4dfe9bebecf8aee168fcb9
Vervoort J D, Pachelt P J, Gehrels G E, et al. Constraints on early Earth differentiation from hafnium and neodymium isotopes[J]. Nature, 1996, (379):624-627. http://www.nature.com/nature/journal/v379/n6566/abs/379624a0.html
Sylvester P J. Post-collisional strongly peraluminous granites[J]. Lithos, 1998, (45):29-44. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0223027258/
Patino D A E, Johnson A D. Phase equilibrium and melting productivity in the politic system:Implication for the origin of peraluminous granitoids and aluminous guanulites[J]. Contributions to Mineralogy and Petrology, 1991, (107):202-218. doi: 10.1007/BF00310707
Vielzeuf D, Montel J M. Partial melting of metagreywackes. Part Ⅰ:Fluid-absent experiments and phase relationships[J]. Contributions to Mineralogy and Petrology, 1994, 117:375-393. doi: 10.1007/BF00307272
Peter D K, Roland M. Lu-Hf and Sm-Nd isotope systems in zircon[J]. Reviews in Mineralogy and Geochemistry, 2003, 53(1):327-341. doi: 10.2113/0530327
Barbarin B. Genesis of the two main types of peraluminous granitoids[J]. Geology, 1996, 24(3):295-298. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1ec164c8ba3726a963883690198e1b39
云南省地质调查院. 1: 25万景洪、勐腊幅区域地质矿产调查报告.2013. 云南省地质调查院. 1: 25万澜沧、勐海幅区域地质矿产调查报告.2013. 云南省地质局.1: 20万勐海幅区域地质报告. 1980.
计量
- 文章访问数: 2769
- HTML全文浏览量: 336
- PDF下载量: 2405
