Geochronology and Trace Element Compositions of Zircon in Granodiorite in North Baibandi Area, Western Bangong-Nujiang Metallogenic Belt and Their Geological Significance
-
摘要:
白板地北部花岗闪长岩体位于班公湖–怒江缝合带北侧改则县境内。野外地质调查表明花岗闪长岩呈岩脉状侵入二叠系龙格组碳酸盐中，并在接触带上发育矽卡岩化蚀变和铜矿化，成矿潜力良好。但目前研究程度较低，尚未开展高精度年代学研究工作。笔者通过对白板地北部花岗闪长岩开展LA-ICP-MS锆石U-Pb年代学和微量元素研究，以精细厘定成岩年龄，揭示成矿温度并分析成矿潜力。结果表明，白板地北部花岗闪长岩锆石均为岩浆成因锆石，17个测点获得的锆石206Pb/208U谐和年龄为(154.8±1.2) Ma（MSWD=1.7），岩体侵位于晚侏罗世。锆石ƩREE值为6.1×10−6～24.04×10−6，平均值11.68×10−6，ƩLREEs和ƩHREE分别为0.41×10−6～9.44×10−6和4.93×10−6～23.55×10−6，具有亏损轻稀土，富集重稀土的特征；δEu值为0.26～0.64，δCe值为0.91～5.03，呈现显著的负Eu异常和正Ce异常特征。锆石Ti元素含量为0.89×10−6～1.43 ×10−6，估算的结晶温度为600.3～799.3 ℃，平均值为697.6℃，一定程度上反映了花岗闪长岩的成岩温度范围。综合成岩年龄和区域构造演化特征，推测白板地北部花岗闪长岩形成于班公湖-怒江洋壳北向俯冲的构造背景，为深入认识区域构造演化提供了新证据。
Abstract:The granodiorite in the north of Baibandi is located in Gaize County, north of the Bangong Nujiang suture zone. The field geological survey shows that granodiorite intruded into the Permian Longge Formation carbonate. Skarn alteration and copper mineralization developed along the contact zone, showing good metallogenic potential. However, due to the absence of high-precision chronology research, the genesis of the granodiorite remains poorly understanding. In this paper, the LA-ICP-MS zircon U-Pb geochronology and trace element studies were carried out on for zircon from the granodiorite in the north of Baibandi area in order to determining the age and analyzing the metallogenic potential. The results show that zircons from the granodiorite in the north of Baibandi area are all magmatic zircons. 17 zircons show weighted average 206Pb/208U ages of (154.8 ± 1.2) Ma (MSWD=1.7), showing that the granodiorite was formed in the Late Jurassic. The ƩREE values zircon of are 6.1×10−6～24.04 ×10−6 (averaging 11.68 ×10−6). The ƩLREEs and ƩHREEs values are 0.41×10−6～9.44×10−6 and 4.93×10−6～23.55×10−6, respectively, showing the enrichment of heavy rare earths. The δEu and δCe values of zircons range from 0.26 to 0.64 and 0.91 to 5.03 respectively, showing significant negative Eu anomalies and positive Ce anomalies. The Ti contents of zircon range from 0.89×10−6 to 1.43×10−6, with estimated crystallization temperature of 600.3 to 799.3 ℃ (averaging 697.6 ℃). Based on these results and the characteristics of regional tectonic evolution, it is inferred that the granodiorite in the north of Baibandi was formed during northward subduction of Bangong-Nujiang oceanic crust. Our results provide new evidence for understanding the formation of the Bangong-Nujiang metallogenic belt.
-
土地是人类生存与发展的重要基础（Zhang et al.，2022），人类依靠土地来获取必须的食物、能源、生存空间和发展，但同时随着人口的持续增长，对自然资源的需求也进一步增长，生态系统会随之受到压力（Song et al.，2018）。土地利用和土地覆盖变化（LUCC）研究计划在1995年被提出，现已成为国际环境变化研究的核心和热点（Lambin et al.，2001；Foley et al.，2005；乔伟峰等，2013；毛盛林等，2022）。土地利用是指人类活动对土地进行开发利用，产生社会效益的过程和活动。土地利用会对土地覆盖产生影响，土地覆盖改变的结果又会作用于人类进行土地利用的过程。（张靖等，2013；Hong et al.，2022）因此，对一个区域内人类土地利用导致的土地覆盖变化情况以及背后驱动因素的研究，可为人类更合理的进行土地利用提供支持（张景华等，2021a）。
土地利用转移矩阵可以分析得到区域各地类的用地结构以及变化方向，在土地利用变化分析中有着广泛的应用。土地利用转移矩阵通过将不同时期各地类发生转移的面积以矩阵的形式列出，可以分析得到各地类的用地结构以及变化方向（刘瑞等，2010）。目前有很多学者使用土地利用转移矩阵对两个时期间的土地利用变化进行分析（吴琳娜等，2014；徐羽等，2017；马彩虹等，2022；Zhou et al.，2022），也有学者在此基础上通过对同一地区多个时段土地利用转移矩阵的时序性分析，运用马尔可夫模型预测未来土地利用的变化趋势（朱会义等，2003）。
大凉山区林业资源丰富，但20世纪大凉山区林地遭到大面积砍伐，严重影响当地生态环境，1998年“天保工程”在凉山开始试点实施，当地生态环境得到逐步恢复（张景华等，2021b）。前人对于大凉山区的研究主要在地质、水文等方面，对当地土地利用变化的研究较少（Liu et al.，2012；邓宾等，2016；祁玉萍等，2018；刘洪等，2020，2022；张腾蛟等，2020；李樋等，2021；欧阳渊等，2021；李琬欣等，2022）。笔者通过遥感手段，解译得到大凉山区1990～2022年土地利用数据，运用土地利用转移矩阵分析当地土地利用变化与生态变化的关系，为大凉山区生态资源管理提供理论依据。
1. 研究区概况
大凉山区位于四川省南部，西跨横断山脉，东达四川盆地，北接川西平原和山地，其地理坐标为E 100°42′～103°53′，N 26°06′～29°27′，在生态地质分区上主要跨越了云贵高原生态地质区和横断山区生态地质区（聂洪峰等，2021），总面积约
30000 km2，是全国最大的彝族聚居区（图1）。地貌类型以山地为主，地势北高南低，东高西低，由EN向WS倾斜，地形地貌明显受地质构造基本格架的控制，多属中深切割的中山地形。山脉河流与构造线相吻合，多呈SN向展布。大凉山区属中亚热带季风气候区，冬季受南支西风环流控制，天气晴朗，日照充足，降雨很少，气候干燥；夏季受暖湿的海洋季风影响，降水丰富，5～9月降水量可以占到全年的90%，雨量常在800～
1200 mm。大凉山区气候虽然属于中亚热带季风气候区，但境内气候复杂多样，既具有南亚热带气候特点，又具有亚热带和寒温带的垂直气候特点。以大小相岭、黄茅埂为界，东北部温凉湿润，西南部炎热干燥，具有南干北湿、西燥东润和低热高凉的特点，垂直差异和水平差异均明显。2. 研究方法
2.1 遥感解译与野外验证
遥感卫星影像是土地利用研究中的重要数据来源之一（韩海辉等，2022）。笔者获取了美国陆地卫星（LANDSAT）的专题绘图仪（Thematic Mapper，TM）、陆地成像仪（Operational Land Imager，OLI）系列数据，数据获取时间等具体信息见表1。再对数据进行几何校正、图像融合、数字图像增强等预处理，得到遥感解译底图。
表 1 遥感数据源Table 1. Remote sensing data source成像日期 卫星 传感器 条带号 行列号 1989.5.11 LANDSAT5 TM 130 41 2000.1.2 LANDSAT5 TM 130 41 2010.3.18 LANDSAT5 TM 130 41 2022.2.15 LANDSAT8 OLI 130 41 由于大凉山区地形地貌复杂，阴影覆盖较多，现有遥感土地利用分类难以满足所需精度，故采用目视解译进行土地利用分类，以一级地类为分类基本单位，将大凉山区土地利用分为耕地、林地、草地、水域湿地、建设用地和未利用地6类，并利用野外调查，建立6种地类遥感解译标志，进行人机交互解译。
得到初步解译结果后，再结合相关地形、地质资料，进行实地调查，采取路线调查和定点观测方法进行工作，采用GPS定点进行观察记录和描述，验证各种遥感解译标志是否准确，重点验证解译各地类边界是否准确，并及时根据验证情况修改遥感解译结果。最后对土地利用解译成果进行随机抽检，准确率达到91.06%。
2.2 土地利用转移矩阵
在各种人为以及自然因素的影响下，在一定时间区域内的土地利用情况会随着时间的变化而发生改变，且不同土地利用类型变化的速度、方向不尽相同。土地利用的变化可以理解为土地利用的空间特征随着时间的变化。除此以外，不同土地利用类型间的互相转化也会受到初始状态的影响。综合来看，土地利用转移矩阵便能够很好的反映出这一变化（张建国等，2018）。
在土地利用转移矩阵（表2）中常用行与列（表1中T1与T2）来表示两个不同研究时期的土地利用类型。在矩阵中的Sij表示的是从T1到T2时期地类i转换为地类j的面积；Sii表示从T1到T2时期保持空间位置、面积不变的地类i的面积（张云亮等，2017）。在矩阵中任一Ai行就表示T1时期的地类i向T2时期的各地类的转移情况，将这一行所有的数据相加也就是T1时期地类i的总面积；任一Aj列就表示T2时期的地类j自T1时期各地类的来源情况，同理，列中所有数据相加就是T2时期地类j的总面积。为了使研究结果更加客观，通常在转移前后使用相同的地类评价标准，这时的土地利用转移矩阵是一个标准的n阶矩阵。
表 2 土地利用转移矩阵Table 2. Land use transfer matrix不同时期土地利用类型 T2 A1 A2 … An T1 A1 S11 S12 … S1n A2 S21 S22 … S2n … … … … … An Sn1 Sn2 … Snn 3. 结果分析
3.1 遥感解译结果
对大凉山区1990年、2000年、2010年和2022年4期的土地利用/覆盖情况开展了遥感解译，以一级类为分类单位，将地物分为耕地、林地、草地、水域湿地、建设用地、未利用地共6种地物类型（刘纪远等，2014），编制了大凉山区4个时期的土地利用/覆盖图（图2），并统计了各地类的总面积情况（表3）。
![]()
图 2 大凉山区1990～2022年土地利用/覆盖图a.1990年土地利用/覆盖图；b.2000年土地利用/覆盖图；c.2010年土地利用/覆盖图；d.2022年土地利用/覆盖图Figure 2. Land use/cover map of Daliangshan area from 1990 to 2022表 3 大凉山区1990～2022年土地利用/覆盖统计表Table 3. Statistical table of land use/cover in Daliangshan area 1990～2022地类 年份 1990年 2000年 2010年 2022年 面积（km2） 比例（%） 面积（km2） 比例（%） 面积（km2） 比例（%） 面积（km2） 比例（%） 耕地 11153.65 36.06 9006.18 29.11 7445.54 24.07 4090.61 13.22 林地 11713.27 37.87 11435.84 36.97 18001.54 58.19 18658.72 60.30 草地 7708.18 24.92 10160.67 32.85 5128.98 16.58 6907.76 22.32 水域湿地 254.81 0.82 174.86 0.57 211.56 0.68 520.75 1.68 建设用地 92.98 0.30 126.17 0.41 146.55 0.47 748.17 2.42 未利用地 11.10 0.04 30.25 0.10 0.15 0.00 8.01 0.03 3.2 大凉山区土地利用变化情况
3.2.1 1990～2000年大凉山区土地利用转移矩阵
从表4中可以总览大凉山区1990～2000年的土地利用总体变化情况，这10年间，耕地转出面积为6 345.55 km2，而转入面积为4 198.09 km2，耕地面积减少趋势较为明显。林地面积转出为4 345.32 km2，转入为4 067.93 km2，面积减少，但变化不大，基本与10年前持平。草地转出面积为3 519.78 km2，但是有5 972.28 km2面积其他土地转为草地，草地面积增加较多。水域湿地转出面积为163.5 km2，转入面积为83.56 km2，面积略微减少。建设用地转出面积为72.13 km2，转入面积为105.31 km2，面积略微增加。
表 4 1990～2000年大凉山区土地利用转移矩阵（km2）Table 4. Land use transfer matrix of Daliangshan area from 1990 to 2000耕地 林地 草地 水域湿地 建设用地 未利用地 合计 耕地 4808.09 2480.85 3725.63 52.36 78.19 8.52 11153.64 林地 2096.85 7367.94 2211.15 14.69 13.11 9.52 11713.26 草地 1934.86 1547.30 4188.39 15.46 10.48 11.68 7708.17 水域湿地 97.69 34.84 27.12 91.31 3.53 0.32 254.81 建设用地 66.11 2.94 2.24 0.64 20.85 0.20 92.98 未利用地 2.58 2.00 6.14 0.41 0 0 11.13 合计 9006.18 11435.87 10160.67 174.87 126.16 30.24 30933.99 3.2.2 2000～2010年大凉山区土地利用转移矩阵
从表5中可以看出大凉山区从2000～2010年的各地类转移情况，这10年间，耕地有5 556.36 km2面积变成其他地类，3 995.72 km2其他类型土地转入，耕地面积继续有较明显减少。林地有2 095.93 km2土地转出，8 661.26 km2土地转入，林地面积有很大幅度的增长。草地有7 255.31 km2土地转出，2 223.62 km2土地转入，草地面积相反有很大幅度的减少。水域湿地有75.06 km2土地转出，111.76 km2土地转入，面积小幅度增长。建设用地有99.03 km2土地转出，119.41 km2土地转入，面积略微增加。
表 5 2000～2010年大凉山区土地利用转移矩阵（km2）Table 5. Land use transfer matrix of Daliangshan area from 2000 to 2010耕地 林地 草地 水域湿地 建设用地 未利用地 合计 耕地 3449.82 4223.71 1196.57 49.02 87.02 0.04 9006.18 林地 1064.06 9339.91 997.40 23.00 11.46 0.01 11435.84 草地 2835.05 4362.10 2905.36 38.70 19.36 0.09 10160.67 水域湿地 23.16 34.50 16.20 99.81 1.17 0.02 174.86 建设用地 61.31 29.86 7.19 0.67 27.14 0.00 126.17 未利用地 12.13 11.09 6.27 0.36 0.40 0.00 30.25 合计 7445.54 18001.18 5128.99 211.56 146.55 0.15 30933.97 3.2.3 2010～2022年大凉山区土地利用转移矩阵
从表6中可以看出2010～2022年大凉山区的土地利用转移情况，在此12年间，耕地面积转出4 268.77 km2，转入913.87 km2，面积持续明显减少。林地面积转出3 592.39 km2，转入4 249.86 km2，面积增长。草地面积转出2 797.05 km2，转入4 575.83 km2，面积增长明显。水域湿地面积转出69.57 km2，转入378.75 km2，面积持续增长。建设用地面积转出50.33 km2，转入651.95 km2，面积增长。
表 6 2010～2022年大凉山区土地利用转移矩阵（km2）Table 6. Land use transfer matrix of Daliangshan area from 2010 to 2022耕地 林地 草地 水域湿地 建设用地 未利用地 合计 耕地 3176.73 1947.80 1698.90 162.45 458.6 1.02 7445.50 林地 504.37 14408.86 2845.71 123.82 114.88 3.61 18001.25 草地 373.44 2263.38 2331.93 88.34 69.02 2.87 5128.98 水域湿地 8.43 27.47 23.96 142.00 9.45 0.26 211.57 建设用地 27.62 11.15 7.19 4.13 96.22 0.24 146.55 未利用地 0.01 0.06 0.07 0.01 0 0.01 0.16 合计 4090.60 18658.72 6907.76 520.75 748.17 8.01 30934.01 4. 讨论
4.1 大凉山区土地利用变化分析
从1990～2022年大凉山区土地利用转移的动态变化过程（图3、图4）来看，发现大凉山区土地利用主要变化量集中在林地、耕地、草地之间，其变化关系反映了当地生态格局的演变（徐润宏等，2023）。而水域湿地虽变化量相对较小，但巨大的生态价值也是衡量当地生态环境的重要指标（谢启姣等，2020）。
![]()
图 3 1990～2022年大凉山区土地利用转移动态过程Figure 3. Dynamic process of land use transfer in Daliang area from 1990 to 20224.1.1 耕地时空分布与演变
2022年大凉山区耕地总面积约4 096 km2（表3），通过实地考察结合遥感影像分析，其主要以旱作农田为主，多分布于各山地半坡部位，其余的水田主要分布于安宁河谷，越西、甘洛、普格等县的河流谷底。
大凉山区除安宁河谷以外，整体经济落后，截止到2022年，研究区域共12市县中有9县为四川省深度贫困县，其主要人类活动为垦殖（耕种），且主要以坡地旱作耕种为主，生产方式落后，带来一系列植被破坏、水土流失等生态问题。但自1998年“天保工程”实施以来，国家号召退耕还林、还草，1990年至2022年大凉山区耕地面积持续减少（图5a）。从1990年到2000年，在1998年“天保工程”实施之前，林地砍伐情况严峻，被砍伐的林地退化为草地、耕地，1998年“天保工程”实施，但时间较短，一时退耕还林、还草的土地难以快速发育为林地，这段时间大凉山区减少的耕地主要向草地转化、其次向林地转化。后随着“天保工程”的进一步实施，从2000～2022年，减少的耕地都主要转化为林地，其次是草地，这说明国家实施的退耕还林、还草政策取得了较好成效，森林与草地的生态系统都得到了较大改善，提高了大凉山区生态系统的稳定性和服务功能。
![]()
图 5 大凉山区耕地（a）、林地（b）、草地（c）和水域湿地（d）1990～2022年变化趋势图Figure 5. Change trend of (a) cultivated Land, (b) forest Land, (c) grassland and (d) wetland in Daliangshan area from 1990 to 20224.1.2 林地时空分布与演变
2022年大凉山区林地面积约为18 658 km2，以常绿针、阔叶林为主，多分布于锦屏山、牦牛山、小相岭、螺髻山、马鞍山、乌科梁子、大凉山等南北走向的主要山脉，有着水土保持（张景华等，2018）、净化空气等重要生态价值，是区内最重要的生态系统。
从1990～2022年，大凉山区林地面积整体变化趋势为持续增长，仅在1990～2000年有微小下降，尤其在2000～2010年涨势明显。大凉山区自古以来森林资源丰富，是中国著名林场，但20世纪中期，大凉山区由于毁林开荒，森林植被破坏较严重，在实地考察中可见许多人工砍伐留下的大型木桩。林地在被砍伐后，较难快速发育回到原本的生态条件，自然主要演变为草地。从表3中可看出，1990～2000年，林地遭到破坏演化为草地的面积达到了约2 211 km2，占到了总林地面积的20 %，但由于期间“天保工程”开始实施，部分耕地、草地也有向林地演化，林地面积只出现了微小下降。2000～2010年，林地得到进一步保护，退耕还林政策的耕地人为转化以及封山育林的草地自然转化面积远远大于林地遭到破坏转出的面积，林地面积大幅增长。2010～2022年，林地面积小幅度增长，主要来源于人为退耕还林的耕地转入，自然条件下的草地、林地相互演化达到平衡。总体来看，1990～2022年，新增林地主要分布于大凉山区周围，尤其是安宁河两岸、邛海盆地增长最多（图6），这里是耕地主要分布区域，说明退耕还林政策效果明显。
![]()
图 6 大凉山区1990～2022年新增林地空间位置分布图Figure 6. Spatial distribution of new forest land in Daliangshan area from 1990 to 20224.1.3 草地时空分布与演变
2022年大凉山区草地面积约6 907 km2，多呈星岛状分布于锦屏山、牦牛山、小相岭、螺髻山、马鞍山、乌科梁子、大凉山等山脉的林地之间，结合数字高程模型进行空间分析后发现，主要分布于海拔较高的中上部，仅在局部，如美姑县南部、金阳县北部有集中连片分布。在实地考察中发现，天然草地存在较少，主要为人工草地和林地砍伐后形成。
从1990～2022年，大凉山区草地面积先增长，后减少，再增长，其特殊的变化规律值得我们注意。1990～2000年，耕地转化为草地的面积要远大于草地转化为耕地的面积，这是此时草地增长的主要原因。2000～2010年，草地面积出现较大幅度下降，这是因为大凉山区的草地主要是人工草地和林地砍伐逆向演替形成，且这个时间段里林地面积快速增加，分析是因为林地得到保护砍伐变少导致草地转入面积大幅度降低，和封山育林草地自然演化为林地所致。2010～2022年，和林地的变化原因相同，此时林、草地相互演化平衡，草地增长来自于退耕还草政策导致的耕地转入。总的来看，虽然草地面积从1990～2022年变化不大，但中间有很大波动，减少的草地主要正向演替为了林地，增加的草地主要得益于退耕还林、还草政策的实施，由耕地转化而来（图7）。
![]()
图 7 大凉山区1990～2022年新增草地空间位置分布图Figure 7. Spatial location of new grasslands in Daliangshan area from 1990 to 20224.1.4 水域湿地时空分布与演变
大凉山区水域湿地面积约520 km2，主要由河流、湖泊、水库坑塘组成。分布上来看，河流主要分布于区内主要江河如金沙江、雅砻江、安宁河等江河干流及主要支流沿岸，湖泊主要分布于邛海、马湖等湖泊，水库坑塘则零星分布。
从1990年到2022年，大凉山区水域湿地面积先出现了小幅度的降低，主要转化为了耕地，可能是少数地区人为破坏所致，后大幅度增加，主要也是来源于耕地的转入。总体涨势明显，水域湿地面积增加了近一倍（图8）。分析原因主要有：① 水利工程设施的修建，导致水面（主要是水库水面）面积有较大幅度增加，如冕宁县大桥水库（图9）。② 湿地保护得到重视，原湿地耕作区域逐步退耕并还原为湿地，区内一些重要湿地得到保护，面积不断扩大，典型如邛海湿地。
![]()
图 8 大凉山区1990～2022年新增水域湿地空间位置分布图Figure 8. Spatial location of wetlands in new waters in Daliangshan area from 1990 to 2022![]()
图 9 冕宁水库1990年（左图）和2022年（右图）影像Figure 9. Image of mianning reservoir in 1990 (left) and 2022 (right)4.2 结果讨论
对土地利用类型变化原因的研究是个系统性的工程，受到自然条件、地理位置、社会经济发展、人口状况、国家或地方政策等诸多因素的影响（Hemmavanh et al.，2010；李璇琼等，2022）。如何去筛选出土地利用类型变化背后主要的驱动因素，利用土地利用转移矩阵计算得到地类的转移面积，结合发生转移地类的时空分布信息，能有效分析出地类发生转移的相关驱动因素。进一步来看，如果能将自然因素、人为因素以及政府政策因素统一标准进行量化分析，便可以确定不同因素对地类变化的影响程度，能对下一步的政策法规起到指导性作用（Wang et al.，2018；赵银兵等，2022）。然而，量化驱动因素的问题尚未被克服，尤其是如何量化评价政府大型政策法规对土地类型变化的影响，需要进一步的研究（Dong et al.，2022）。
政府的政策法规主要会对诸多人为因素产生影响（唐亚明等，2015）。大凉山区在“天保工程”实施之前，虽然林地资源丰富，但经过长期的人工砍伐，已经产生了很大的破坏，由此带来一系列生态问题（Wang et al.，2020）。自1998年国家“天保工程”计划实施以来，大凉山区人为的进行封山育林以及退耕还林、还草，保护林地生态资源，从1990～2022年，林地面积增长幅度达到约58%（图10），效果显著，尤其是2000年后林地资源恢复明显。政府政策的引导主要是对人类活动做出干扰，大凉山区林地资源丰富，但经济发展水平落后，在实地考察中，仍可见许多粗壮树木被砍伐后留下的树桩。1990～2000年，虽然林地面积总体平稳，但林地遭到破坏面积很大，后得益于1998年国家“天保工程”以及退耕还林政策实施，林地资源得到保护。深究原因，早在1984年，国家就颁布了森林法（樊宝敏等，2021），走上依法治林，但当地经济基础薄弱，农民为了维持生计，砍伐树木。在1998年，“天保工程”开始实施，实施天然林保护并人工种植树木，2001年，建立中央森林生态效益补偿基金，1998～2010年，国家在大凉山区投入“天保工程”专项资金17.76亿元，因此2000～2010年，大凉山区林地面积快速增长。2010年后，国家相关资金减少，大凉山区林业发展转为维护已有成果，林地面积增长放缓。
![]()
图 10 1990～2022年大凉山区林地覆盖率变化Figure 10. Change of forest land coverage in Daliangshan area from 1990 to 2022最后，在大凉山区生态发展良好的同时，还存在一些问题。大凉山区林地多以人工种植为主，且多为针叶林，当地原生阔叶林较少（张景华等，2021c），不利于当地生物多样性的保护以及林地生态功能的发挥。针叶林还容易导致火灾，统计显示，仅在西昌市云南松林区每年发生火灾面积2～5 km2。大凉山区海拔高，火灾扑救难度大，2019年木里县“3•30”火灾就导致了大量经济损失与30人遇难。
5. 结论
（1）1990～2022年，大凉山区林地面积大幅度增长，耕地面积持续减少，且两者在时空分布上吻合，说明增长的林地主要由耕地转化而来。草地面积因转入面积与转出面积相对平衡，整体变化不明显，转出的草地主要通过人工封山育林正向演替为林地，转入的草地主要是退耕还林、还草政策导致的耕地面积转入。水域湿地面积得益于水利设施的建设与湿地保护；建设用地伴随着大凉山区经济的发展，两者面积均有较大幅度增长，大凉山区“天保工程”等国家生态工程的实施实现了人与自然的双赢。但实地考察发现仍存在一定问题：大量人工林存在着火灾等安全隐患、部分林地砍伐后逆向演替为草地等，这些仍需相关部门需要多加监管注意。
（2）笔者基于土地利用转移矩阵，对大凉山区土地利用变化展开研究，结合相关政策背景，着重分析了大凉山区土地利用变化情况与当地生态格局演变的关系，以及相关政策对土地利用变化的影响，可以为大凉山区进一步的区域规划管理与政策制定提供支持。但想要进一步研究政策法规等与土地利用变化的量化关系，还需要更多的社会、经济等信息支撑，建立一个完整的评价体系。
-
![]()
图 3 白板地北部花岗闪长岩锆石阴极发光（CL）图像及测点图
Figure 3. The CL images and the laser location of zircon from the granodiorite in the north of Baibandi area
![]()
图 4 花岗闪长岩中锆石稀土元素配分模式图
Figure 4. The chondrite-normalized REE pattern of zircon from the granodiorite
![]()
图 5 岩浆锆石-热液锆石判别图解（Hoskin, 2005）
Figure 5. La-(Sm/La)N diagram for zircon
![]()
图 6 白板地北部花岗闪长岩锆石U-Pb年龄谐和图
Figure 6. The U-Pb Concordia age of the granodiorite in the north of Baibandi area
表 1 白板地矿区外围花岗闪长岩锆石微量元素(10–6)分析结果
Table 1 The trace element (10–6) result of zircon from granodiorite in Baibandi area
测点 Ti Y Nb Ta La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu ∑REE δEu δCe 1 1.20 12.01 0.08 0.04 0.02 0.14 0.01 0.10 0.15 0.04 0.42 0.11 0.97 0.44 2.00 0.48 5.51 1.13 11.53 0.47 2.89 2 1.13 14.14 0.09 0.03 2.26 4.37 0.46 1.93 0.37 0.05 0.38 0.10 1.23 0.49 2.25 0.55 5.50 1.04 20.98 0.39 0.98 3 0.98 12.32 0.09 0.04 0.09 0.27 0.03 0.23 0.17 0.04 0.51 0.13 1.25 0.47 1.88 0.47 3.79 0.65 9.98 0.38 1.23 4 1.14 5.56 0.08 0.05 0.14 0.30 0.04 0.18 0.14 0.03 0.32 0.07 0.63 0.22 1.00 0.28 2.93 0.64 6.92 0.48 0.91 5 1.33 17.87 0.09 0.04 0.01 0.19 0.01 0.09 0.14 0.04 0.46 0.16 1.61 0.63 2.78 0.59 5.59 1.02 13.32 0.41 4.28 6 1.41 9.19 0.07 0.04 0.01 0.18 0.01 0.12 0.16 0.05 0.38 0.10 0.91 0.34 1.59 0.35 3.44 0.70 8.36 0.64 5.03 7 1.43 7.98 0.07 0.03 0.01 0.15 0.02 0.24 0.24 0.06 0.62 0.13 1.01 0.33 1.40 0.41 3.84 0.83 9.29 0.49 2.30 8 0.94 4.69 0.07 0.04 0.01 0.11 0.01 0.10 0.15 0.03 0.29 0.07 0.60 0.19 0.79 0.26 2.81 0.68 6.08 0.41 2.36 9 1.10 31.51 0.10 0.04 0.01 0.15 0.01 0.11 0.17 0.04 0.78 0.22 2.67 1.05 4.99 1.12 10.76 1.96 24.03 0.26 3.28 10 1.25 21.91 0.09 0.04 0.01 0.17 0.01 0.12 0.16 0.05 0.47 0.14 1.81 0.73 3.31 0.82 8.09 1.74 17.63 0.47 4.47 11 1.06 8.81 0.08 0.03 0.01 0.14 0.01 0.13 0.17 0.04 0.36 0.09 0.82 0.31 1.62 0.43 4.79 1.09 10.00 0.45 2.85 12 1.24 8.60 0.08 0.04 0.01 0.12 0.01 0.10 0.16 0.04 0.40 0.10 0.93 0.36 1.45 0.34 3.32 0.66 8.00 0.43 3.51 13 1.24 4.10 0.09 0.04 0.01 0.14 0.01 0.10 0.14 0.04 0.32 0.07 0.59 0.20 0.92 0.26 2.81 0.63 6.25 0.61 2.88 14 0.89 4.32 0.05 0.03 0.27 0.55 0.06 0.32 0.12 0.03 0.28 0.06 0.49 0.20 0.88 0.23 2.34 0.45 6.28 0.41 0.97 15 1.09 11.60 0.09 0.04 0.02 0.16 0.01 0.11 0.16 0.04 0.49 0.12 1.01 0.42 1.87 0.48 4.89 0.93 10.70 0.39 2.29 16 1.16 27.24 0.08 0.03 0.01 0.13 0.01 0.13 0.16 0.04 0.62 0.22 2.53 0.89 4.10 0.89 8.77 1.60 20.09 0.35 3.10 17 1.25 10.00 0.08 0.04 0.01 0.20 0.01 0.14 0.17 0.05 0.41 0.09 0.94 0.32 1.50 0.40 3.99 0.86 9.08 0.58 3.41 
下载: 导出CSV
表 2 白板地矿区外围花岗闪长岩样品锆石U-Pb定年结果
Table 2 The U-Pb results of zircon from granodiorite in Baibandi area
测点号 含量（10−6） 同位素比值及误差 年龄（Ma）及误差 238U 232Th U/Th 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 206Pb/238U 1σ 1 138.54 366.17 2.64 0.0536 0.0022 0.1811 0.0075 0.0246 0.0002 156.64 1.51 2 148.67 397.41 2.67 0.0518 0.0020 0.1711 0.0062 0.0240 0.0002 153.01 1.50 3 130.68 295.38 2.26 0.0504 0.0024 0.1631 0.0075 0.0236 0.0003 150.12 1.80 4 203.13 472.42 2.33 0.0507 0.0018 0.1726 0.0061 0.0247 0.0002 157.14 1.44 5 68.52 144.69 2.11 0.0469 0.0027 0.1525 0.0084 0.0240 0.0004 152.62 2.28 6 64.27 105.97 1.65 0.0605 0.0036 0.1938 0.0109 0.0241 0.0004 153.25 2.42 7 107.75 184.78 1.71 0.0543 0.0028 0.1821 0.0096 0.0246 0.0004 156.69 2.44 8 165.62 400.21 2.42 0.0527 0.0022 0.1736 0.0066 0.0243 0.0003 154.55 1.91 9 107.27 314.87 2.94 0.0459 0.0022 0.1506 0.0071 0.0240 0.0003 152.73 1.86 10 144.11 333.98 2.32 0.0489 0.0021 0.1596 0.0067 0.0240 0.0003 152.82 1.68 11 118.63 294.75 2.48 0.0508 0.0023 0.1714 0.0074 0.0248 0.0003 158.00 1.84 12 80.83 216.32 2.68 0.0516 0.0028 0.1696 0.0090 0.0240 0.0003 152.90 1.99 13 125.69 343.74 2.73 0.0506 0.0020 0.1684 0.0067 0.0241 0.0003 153.82 1.69 14 83.72 257.36 3.07 0.0523 0.0024 0.1737 0.0075 0.0245 0.0003 155.85 1.77 15 88.07 238.58 2.71 0.0532 0.0026 0.1798 0.0084 0.0248 0.0003 157.90 2.10 16 73.90 206.91 2.80 0.0580 0.0027 0.1982 0.0095 0.0249 0.0004 158.38 2.27 17 93.93 157.44 1.68 0.0518 0.0026 0.1756 0.0093 0.0245 0.0004 156.32 2.29
下载: 导出CSV
表 3 锆石Ti含量温度计算结果
Table 3 Ti LA-ICP-MS zircon data and TTi-in-zircon thermometry calculation results
点号 1 2 3 4 5 6 7 8 9 T（℃） 600.3 609.0 607.5 617.0 638.3 665.4 695.3 690.0 698.2 点号 10 11 12 13 14 15 16 17 T（℃） 718.0 723.3 743.4 764.3 752.9 761.3 776.2 799.3
下载: 导出CSV
-
常青松, 朱弟成, 赵志丹, 等. 西藏羌塘南缘热那错早白垩世流纹岩锆石U-Pb年代学和Hf同位素及其意义[J]. 岩石学报, 2011, 27（7）: 2034−2044. CHANG Qingsong, ZHU Dicheng, ZHAO Zhidan, et al. Zircon U-Pb geochronology and Hf isotopes of Early Cretaceous rhyolites from the southern margin of the Qiangtang area, Tibet, and their significance[J]. Acta Petrologica Sinica,2011,27(7):2034−2044.
陈华安, 祝向平, 马东方, 等. 西藏波龙斑岩铜金矿床成矿斑岩年代学、岩石化学特征及其成矿意义[J]. 地质学报, 2013, 87（10）: 1593−1611. CHEN Hua'an, ZHU Xiangping, MA Dongfang, et al. Chronology, petrology, and metallogenic significance of mineralizing porphyries in the Bolong porphyry Cu-Au deposit, Tibet[J]. Acta Geologica Sinica,2013,87(10):1593−1611.
陈澍民, 徐宏根, 吴金虹, 等. 湖南天明金矿区云斜煌斑岩年代学、同位素地球化学及成矿意义[J]. 西北地质, 2023, 56（6）: 285−300. CHEN Shumin, XU Honggen, WU Jinhong, et al. Geochronology, Isotopic Geochemistry of Diorite Porphyrite in Tianming Gold Deposit, Hunan[J]. Northwestern Geology,2023,56(6):285−300.
冯国胜, 廖六根, 陈振华, 等. 西藏西部日土县材玛铁多金属矿地质特征及找矿意义[J]. 地质通报, 2006, 25（S1）: 267−272. FENG Guosheng, LIAO Liugen, CHEN Zhenhua, et al. Geological characteristics and prospecting significance of the Caimar iron-polymetallic deposit in the western part of Tibet, China[J]. Geological Bulletin of China,2006,25(S1):267−272.
胡为正, 廖辉宝, 黄东荣. 西藏日土县材玛铁矿地质特征及找矿方向[J]. 资源调查与环境, 2014, 35（2）: 120−129. HU Weizheng, LIAO Huibao, HUANG Dongrong. Geological characteristics and prospecting direction of the Caimar iron deposit in Ritong County, Tibet[J]. Resources Survey and Environment,2014,35(2):120−129.
代新宇, 周斌, 李新林, 等. 西昆仑奇台达坂北中新世石英二长岩侵入岩年代学、地球化学及其构造意义[J]. 西北地质, 2024, 57（4）: 191−205. DAI Xinyu,ZHOU Bin,LI Xinlin,et al. Geochronology, Geochemistry and Tectonic Significance of Miocene Quartz Monzonite from the Northern of Qitai Mountain in Western Kunlun[J]. Northwestern Geology,2024,57(4):191−205.
李金祥, 李光明, 秦克章, 等. 班公湖带多不杂富金斑岩铜矿床斑岩-火山岩的地球化学特征与时代: 对成矿构造背景的制约[J]. 岩石学报, 2007, 24（3）: 531−543. LI Jinxiang, LI Guangming, QIN Kezhang, et al. Geochemical characteristics and age of porphyry-volcanic rocks in the Duobuza rich-gold porphyry copper deposit in the Bangong Lake area: Implications for metallogenic tectonic settings[J]. Acta Petrologica Sinica,2007,24(3):531−543.
潘桂棠, 莫宣学, 侯增谦, 等. 冈底斯造山带的时空结构及演化[J]. 岩石学报, 2007, 22（3）: 521−533. PAN Guitang, MO Xuanxue, HOU Zengqian, et al. The spatial-temporal structure and evolution of the Gangdese orogenic belt[J]. Acta Petrologica Sinica,2007,22(3):521−533.
李平, 朱涛, 吕鹏瑞, 等. 西天山早寒武世夏特辉长岩: 南天山洋早期俯冲的岩浆记录[J]. 西北地质, 2024, 57（3）: 44−58. LI Ping,ZHU Tao,LÜ Pengrui,et al. Early Cambrian Xiate Gabbro in Western Tianshan: Magmatic Records of Initial Subduction of the South Tianshan Ocean[J]. Northwestern Geology,2024,57(3):44−58.
潘桂棠, 王立全, 朱弟成. 青藏高原区域地质调查中几个重大科学问题的思考[J]. 地质通报, 2004, 23（1）: 12−19. PAN Guitang, WANG Liquan, ZHU Dicheng. Reflections on some major scientific issues in regional geological surveys of the Qinghai-Tibet Plateau[J]. Geological Bulletin of China,2004,23(1):12−19.
曲晓明, 王瑞江, 辛洪波, 等. 西藏西部与班公湖特提斯洋盆俯冲相关的火成岩年代学和地球化学[J]. 地球化学, 2009, 38（6）: 523−535. doi: 10.3321/j.issn:0379-1726.2009.06.002 QU Xiaoming, WANG Ruijiang, XIN Hongbo, et al. Geochronology and geochemistry of igneous rocks associated with the subduction of the Bangong Lake Tethys Ocean Basin in western Tibet[J]. Geochimica,2009,38(6):523−535. doi: 10.3321/j.issn:0379-1726.2009.06.002
曲晓明, 辛洪波. 藏西班公湖斑岩铜矿带的形成时代与成矿构造环境[J]. 地质通报, 2006, 25（7）: 792−799. doi: 10.3969/j.issn.1671-2552.2006.07.004 QU Xiaoming, XIN Hongbo. Formation age and metallogenic tectonic environment of the Bangong Lake porphyry copper belt in western Tibet[J]. Geological Bulletin of China,2006,25(7):792−799. doi: 10.3969/j.issn.1671-2552.2006.07.004
佘宏全, 李进文, 马东方, 等. 西藏多不杂斑岩铜矿床辉钼矿Re-Os和锆石U-Pb SHRIMP测年及地质意义[J]. 矿床地质, 2009, 28（6）: 737−746. doi: 10.3969/j.issn.0258-7106.2009.06.003 SHE Hongquan, LI Jinwen, MA Dongfang, et al. Molybdenite Re-Os and zircon U-Pb SHRIMP dating and their geological significance in the Duobuza porphyry copper deposit, Tibet[J]. Mineral Deposits,2009,28(6):737−746. doi: 10.3969/j.issn.0258-7106.2009.06.003
宋扬, 唐菊兴, 曲晓明, 等. 西藏班公湖-怒江成矿带的广义概念及成矿特点[J]. 矿床地质, 2014, 33（S1）: 815−816. SONG Yang, TANG Juxing, QU Xiaoming, et al. General concept and metallogenic characteristics of the Bangong Lake-Nujiang metallogenic belt, Tibet[J]. Mineral Deposits,2014,33(S1):815−816.
唐菊兴. 青藏高原及邻区重要成矿带矿产资源基地调查与研究进展[J]. 岩石学报, 2019, 35（3）: 617−624. doi: 10.18654/1000-0569/2019.03.01 TANG Juxing. Progress in the investigation and research of mineral resource bases in important metallogenic belts in the Qinghai-Tibet Plateau and adjacent areas[J]. Acta Petrologica Sinica,2019,35(3):617−624. doi: 10.18654/1000-0569/2019.03.01
王立强, 王勇, 旦真王修, 等. 班公湖—怒江成矿带西段主要岩浆热液型矿床成矿特征初探[J]. 地球学报, 2017, 38（5）: 615−626. doi: 10.3975/cagsb.2017.05.03 WANG Liqiang, WANG Yong, DANZHEN Wangxiu, et al. Preliminary study on metallogenic characteristics of major magmatic hydrothermal deposits in the western Bangong Lake-Nujiang metallogenic belt[J]. Acta Geoscientica Sinica,2017,38(5):615−626. doi: 10.3975/cagsb.2017.05.03
王立强, 谢富伟, 王勇. 西藏巴嘎拉东铅锌矿床黑云母花岗岩锆石U-Pb年龄、微量元素组成及地质意义[J]. 岩矿测试, 2016, 35（6）: 650−657. WANG Liqiang, XIE Fuwei, WANG Yong. Zircon U-Pb age, trace element composition, and geological significance of biotite granite in the Bagala East Pb-Zn deposit, Tibet[J]. Rock and Mineral Analysis,2016,35(6):650−657.
王新雨, 王书来, 吴锦荣, 等. 青海省牛苦头铅锌矿床成矿时代研究: 来自成矿岩体年代学和黄铁矿Re–Os地球化学证据[J]. 西北地质, 2023, 56（6）: 71−81. WANG Xinyu, WANG Shulai, WU Jinrong, et al. Mineralization Age and Ore forming–Source of Niukutou Pb–Zn Deposit, Qinghai: Evidence from Geochronology of Ore–forming Rock Bodies and Re–Os Geochemistry of Pyrite[J]. Northwestern Geology,2023,56(6):71−81.
吴元保, 郑永飞. 锆石成因矿物学研究及其对U-Pb年龄解释的制约[J]. 科学通报, 2004, 49（16）: 1589−1604. doi: 10.3321/j.issn:0023-074X.2004.16.002 WU Yuanbao, ZHENG Yongfei. Genetic mineralogy of zircon and its constraints on U-Pb age interpretation[J]. Chinese Science Bulletin,2004,49(16):1589−1604. doi: 10.3321/j.issn:0023-074X.2004.16.002
周金胜, 孟祥金, 臧文栓, 等. 西藏青草山斑岩铜金矿含矿斑岩锆石U-Pb年代学、微量元素地球化学及地质意义[J]. 岩石学报, 2013, 29（11）: 3755−3766. ZHOU Jinsheng, MENG Xiangjin, ZANG Wenshuan, et al. Zircon U-Pb geochronology, trace element geochemistry, and geological significance of mineralized porphyries in the Qingcaoshan porphyry Cu-Au deposit, Tibet[J]. Acta Petrologica Sinica,2013,29(11):3755−3766.
AndersenT, Correction of common lead in U-Pb analyses that do not report Pb-204[J]. Chemical Geology,2002,192:59−79. doi: 10.1016/S0009-2541(02)00195-X
Belousova E A, Griffin W L, Suzanne Y, et al. Igneous zircon: trace element composition as an indicator of source rock type[J]. Contributions to Mineralogy and Petrology,2002,143(5):602−622. doi: 10.1007/s00410-002-0364-7
Cao M J, Qin K Z, Li G M, et al. Tectono-magmatic evolution of Late Jurassic to Early Cretaceous granitoids in the west central Lhasa subterrane, Tibet[J]. Gondwana Research,2016,39:386−400. doi: 10.1016/j.gr.2016.01.006
Fan J J, Li C, Xie X M, et al. Petrology and U–Pb zircon geochronology of bimodal volcanic rocks from the Maierze Group, northern Tibet: Constraints on the timing of closure of the Banggong-Nujiang Ocean[J]. Lithos,2015,227:148−160. doi: 10.1016/j.lithos.2015.03.021
Guynn, J H , Kapp P , Pullen A, et al. Tibetan basement rocks near Amdo reveal “missing” Mesozoic tectonism along the Bangong suture central Tibet[J]. Geology,2006,34:505−508.
Hao L L, Wang Q, Wyman D A, et al. Underplating of basaltic magmas and crustal growth in a continental arc: Evidence from Late Mesozoic intermediate-felsic intrusive rocks in southern Qiangtang, central Tibet[J]. Lithos, 2016, 245: 223-242.
Hoskin P W O, Schaltegger U. The composition of zircon and igneous and metamorphic petrogenesis[J]. Reviews in Mineralogy and Geochemistry,2003,53(1):27−62. doi: 10.2113/0530027
Hoskin P W O. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills[J]. Australia. Geochimica et Cosmochimica Acta,2005,69(3):637−648. doi: 10.1016/j.gca.2004.07.006
Li S M, Zhu D C, Wang Q, et al. Northward subduction of Bangong-Nujiang Tethys: Insight from Late Jurassic intrusive rocks from Bangong Tso in western Tibet[J]. Lithos,2014,205:284−297.
Li Y L, He J, Han Z P, et al. Late Jurassic sodium-rich adakitic intrusive rocks in the southern Qiangtang terrane, central Tibet, and their implications for the Bangong–Nujiang Ocean subduction[J]. Lithos,2016,245:34−46.
Ludwig, K. User's manual for Isoplot 3.00: a geochronological toolkit for Microsoft Excel[J]. Berkeley Geochronology Center Special Publication,2003,4:1−70.
Shi R D, Yang J S, Xu Z Q, et al. The Bangong Lake ophiolite (NW Tibet) and its bearing on the tectonic evolution of the Bangong–Nujiang suture zone[J]. Journal of Asian Earth Sciences,2008,32(5):438−457.
Shi R D. SHRIMP dating of the Bangong Lake SSZ-type ophiolite: Constraints on the closure time of ocean in the Bangong Lake-Nujiang River, northwestern Tibet[J]. Chinese Science Bulletin,2007,52(7):936−941. doi: 10.1007/s11434-007-0134-z
Watson E B, Harrison T M. Zircon thermometer reveals minimum melting conditions on earliest Earth[J]. Science,2005,308:841−844.
Watson E B, Wark D A, Thomas J B. Crystallization thermometers for zircon and rutile[J]. Contributions to Mineralogy and Petrology,2006,151:413−433.
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 Newsletter,1995,19(1):1−23. doi: 10.1111/j.1751-908X.1995.tb00147.x
Wiedenbeck M, Hanchar J M, Peck W H, et al. Further characterisation of the 91500 zircon crystal[J]. Geostandards and Geoanalytical Research,2004,28(1):9−39. doi: 10.1111/j.1751-908X.2004.tb01041.x
Yuan H, Gao S, Liu X, et al. Accurate U-Pb age and trace element determinations of zircon by laser ablation-inductively coupled plasma-mass spectrometry[J]. Geostandards and Geoanalytical Research,2004,28(3):353−370. doi: 10.1111/j.1751-908X.2004.tb00755.x
Zhu D C, Li S M, Cawood P A, et al. Assembly of the Lhasa and Qiangtang terranes in central Tibet by divergent double subduction[J]. Lithos,2016,245:7−17.
计量
- 文章访问数: 65
- HTML全文浏览量: 8
- PDF下载量: 41
