-
摘要:
长江经济带发展战略实施以来，长江上、中、下游地区经济快速发展，与此同时，生态系统不可避免受到扰动，长江经济带生态环境状况面临压力，亟需开展长江经济带生态环境质量评价。本文选取2001年和2020年MODIS和Landsat遥感影像计算湿度、绿度、干度、热度4个评价指标，通过主成分分析方法计算遥感生态指数RSEI，对长江经济带11省（市）生态环境现状及变化进行定量分析评价。结果表明长江经济带整体生态环境等级为较优，2001—2020年保持总体向好态势。空间分布上，生态环境优等区域主要分布于武夷山、罗霄山、雪峰山、武陵山、巫山、大巴山、大别山、西双版纳等山地丘陵一带，生态环境较差及差等区域主要分布于大中小城镇等人类聚集地、川西横断山区。时空变化上，生态环境等级下降区域主要集中在长三角太湖周边、江汉平原、洞庭湖周边、安徽西北部、湖南南部及四川盆地等人类聚集地周边，生态环境等级上升区域主要集中在四川盆地大部、重庆山区及三峡库区、贵州大部、云南东南部、安徽中部及北部、江苏北部等地区。
Abstract:Since the implementation of the Yangtze River Economic Zone strategy, the economy of provinces and cities along the Yangtze River has developed rapidly. At the same time of economic growth, the ecosystem has been strongly disturbed. The ecological environment of the Yangtze River Economic Zone is facing a severe challenge, which seriously restricts the sustainable development of the Yangtze River Economic Zone. The MODIS and Landsat remote sensing images in 2001 and 2020 were selected as the evaluation indexes of humidity, greenness, dryness and heat to calculate the Remote Sensing Ecological Index (RSEI) by means of principal component analysis method, and finally to quantitatively analyze and evaluate the ecological environment status and changes of 11 provinces and cities in the Yangtze River Economic Zone. The results show that the overall ecological environment grade of the Yangtze River economic zone is optimal, and especially it maintained the overall stable and slightly better trend from 2001 to 2020. The optima ecological environment regions are mainly distributed in Wuyi Mountain, Luoxiao Mountain, Xuefeng Mountain, Wuling Mountain, Wushan Mountain, Dabie Mountain, Xishuangbanna and other areas, while the poor and poor ecological environment regions are mainly distributed in large, medium and small towns, Hengduan Mountain Area in Western Sichuan. The eco-environment grade declining areas are mainly concentrated in the Yangtze River Delta, Taihu Lake, Jianghan Plain, Dongting Lake, Northwest Anhui, southern Hunan and Sichuan Basin. While, the eco-environment grade rising areas are mainly concentrated in most of Sichuan Basin, Chongqing mountain area and Three Gorges Reservoir area, most of Guizhou, Southeast Yunnan, central and Northern Anhui and Northern Jiangsu.
-
1. 引言
蛇绿岩是分布在板块缝合带的一套特殊类型的镁铁质-超镁铁质岩石组合，被认为是仰冲就位于陆地的大洋岩石圈板块残片，是板块汇聚最重要的标志之一。由于蛇绿岩的岩石单元可以和现代大洋岩石圈对比，记录了大洋的形成、扩张、消减到碰撞等一系列过程，因此蛇绿岩被认为蕴含了地球内部动力学和板块运动过程的丰富信息，是研究大洋岩石圈演化的最佳载体（Dilek and Furnes, 2011；Rampone and Hofmann, 2012）。
西藏雅鲁藏布构造带是印度板块和欧亚板块的分界线。狭窄的缝合带不连续分布有一系列蛇绿岩体，西段以大规模出露新鲜地幔橄榄岩和非常少的镁铁质岩石为特征，中段以普遍蛇纹石化的地幔橄榄岩伴生相对完整的镁铁质岩石单元为特征，东段以较少出露镁铁质岩石且地幔橄榄岩中赋存大规模铬铁矿床为特征。然而，雅鲁藏布构造带蛇绿岩体的形成环境目前仍然存在较大的争议（Bao et al., 2013）。一些学者认为雅鲁藏布蛇绿岩形于洋中脊慢速扩张中心（Nicolas et al., 1981；Allegre et al., 1984；Girardeau and Mercier, 1988；吴福元等，2014；Liu et al., 2014；Zhang et al., 2016；Zhang et al., 2016, 2017；Liu et al., 2016）。然而，很多学者认为雅鲁藏布蛇绿岩形成于俯冲带环境（张旗等2003；Dubois-Cote et al., 2005；钟立峰等，2006；Bedard et al., 2009；Hebert et al., 2012；Dai et al., 2012, 2013；李文霞等，2012）。本研究展示新的地质和地球化学证据进一步探讨雅鲁藏布构造带中段日喀则蛇绿岩的形成构造环境。
2. 地质背景
青藏高原由不同的地体拼接而成，由南向北包括拉萨地体、羌唐地体、松潘—甘孜地体、昆仑—祁连地体（Yin and Harrison, 2000）。雅鲁藏布构造带位于拉萨地体和印度大陆之间，一系列断续出露的蛇绿岩体东西向展布超过2500 km（图 1a）（Gansser, 1980）。前人对这些蛇绿岩体中镁铁质岩石开展了大量高精度U-Pb锆石年代学工作，定年结果表明这些镁铁质岩石形成时代在118~149 Ma范围内（图 1a）。根据地理分布和岩体特征，雅鲁藏布蛇绿岩带从西向东可以分为3段：西段、中段和东段（图 1a）。
![]()
图 1 雅鲁藏布构造带及蛇绿岩简图Figure 1. Simplified geological map of Yarlung -Zangbo suture zone and ophiolitic massifsa-Distribution of the Yarlung Zangbo ophiolites and their zircon U-Pb ages from mafic rocks (modified from An et al., 2014; Wu Fuyuan et al., 2014); b-Distribution of the Xigaze ophiolitic massifs and sampling locations of this study (modified from Wang Xibin et al., 1987)雅鲁藏布西段蛇绿岩主要以东波、普兰、休古嘎布、仲巴、萨嘎、桑桑等岩体为代表，以出露大量地幔橄榄岩及少量镁铁质岩石为特征（徐德明等，2007；Bedard et al., 2009；刘钊等，2011；Bezard et al., 熊发挥等，2011；刘飞等，2013；2011; Liu et al., 2014; Liu et al., 2015）。东波、普兰、休古嘎布、仲巴岩体的地幔橄榄岩主要以新鲜的方辉橄榄岩为主，含少量二辉橄榄岩和纯橄岩，局部出露小规模豆荚状铬铁矿化。而萨嘎和桑桑岩体的地幔橄榄岩强烈蛇纹石化（Bedard et al., 2009）。这些岩体的地幔橄榄岩中普遍被不同规模的辉长岩和辉绿岩脉侵入，局部出露少量玄武岩。
雅鲁藏布中段蛇绿岩分布于昂仁到仁布地区，从西向东包括昂仁、吉定、路曲、下鲁、冲堆、群让、得几、白朗、大竹曲和仁布等岩体，统称为日喀则蛇绿岩（图 1a）。日喀则蛇绿岩发育最为典型，研究程度也最高（吴福元等，2014）。日喀则蛇绿岩体以广泛出露地幔橄榄岩和相对完整的镁铁质岩石层序为特征（Nicolas et al., 1981; 王希斌等，1987）。蛇绿岩层序从南向北，大致以地幔橄榄岩-辉长岩-辉绿岩-玄武岩-硅质岩的岩性变化。新鲜的地幔橄榄岩主要分布在路曲和大竹曲岩体，主要由方辉橄榄岩和少量二辉橄榄岩和纯橄岩组成，部分橄榄岩中含豆荚状铬铁矿化，而其他岩体的地幔橄榄岩多已不同程度的蛇纹石化。日喀则蛇绿岩的吉定、白朗、大竹曲等岩体发育洋壳，主要有层状辉长岩、块状辉长岩、辉绿岩席、枕状熔岩组成。而昂仁、路曲、下鲁、冲堆、群让、得几、仁布等岩体则缺失辉长岩，仅仅包含相对较薄的辉绿岩席和枕状熔岩。在壳幔过渡带可见一定规模的辉绿岩和辉长岩脉侵入到下覆的地幔橄榄岩中（陈根文等, 2003, 2015；牛晓露等，2006；李文霞等，2012；Liu et al., 2016；Zhang et al., 2016；Zhang et al., 2017）。
雅鲁藏布东段蛇绿岩主要包括泽当、罗布莎等岩体，以大量地幔橄榄岩和发育纯橄岩和铬铁矿床为特征。泽当岩体由地幔橄榄岩及一套镁铁质杂岩组成，但是与日喀则蛇绿岩不同的是，泽当岩体的地幔橄榄岩部分位于北侧，南侧出现少量辉绿岩及硅质岩。罗布莎岩体主要由地幔橄榄岩、辉长岩-辉绿岩、玄武岩组成。辉绿岩和辉长岩规模延伸较大，位于地幔橄榄岩体的北侧，玄武岩不发育，少量玄武岩位于地幔橄榄岩体的北侧（韦栋梁等，2004；钟立峰等，2006；Zhang et al., 2016；Xiong et al., 2016；刘良志等，2018）。
3. 野外地质和岩相学特征
3.1 路曲岩体
路曲岩体位于日喀则市区向南10 km左右（图 1b）。岩体北侧与古近—新近纪日喀则组磨拉石建造和第四纪堆积物呈构造接触，岩体南侧与厚层的三叠纪嘎学群复理石建造呈断层接触。路曲岩体最底部为蛇绿混杂岩，主要为强烈蛇纹石化的地幔橄榄岩和少量异剥钙榴岩脉组成（图 2a）。新鲜的橄榄岩主要位于地幔橄榄岩段中部，岩性为方辉橄榄岩和少量二辉橄榄岩。壳幔过渡带橄榄岩强烈蛇纹石化，其中被少量不同规模的辉绿岩脉侵入（图 2b）。路曲岩体镁铁质岩石主要为辉绿岩和少量枕状玄武岩（图 2f），缺失辉长岩。辉绿岩位于地幔橄榄岩和玄武岩之间，与上覆玄武岩过渡接触，与下覆地幔橄榄岩呈侵入接触关系。玄武岩呈枕状和气孔-杏仁状构造，主要位于岩体北侧，与顶部硅质岩整合接触。
![]()
图 2 日喀则蛇绿岩野外照片a—路曲岩体蛇绿混杂岩；b—路曲岩体辉绿岩脉侵入地幔橄榄岩（蛇纹岩）；c—大竹曲岩体层状辉长岩；d—大竹曲岩体块状辉长岩和辉绿岩；e—大竹曲岩体含斜长石斑晶的辉绿岩；f—路曲岩体枕状玄武岩Figure 2. Field photos of the Xigaze ophiolitea-Ophiolitic mélange of the Luqu massif; b-Diabase dike that intruded into mantle peridotite (serpentinite) of the Luqu massif; c-Layered gabbro of the Dazhuqu massif; d-Massive gabbro and diabase of the Dazhuqu massif; e-Plagioclase phenocrysts in diabase of the Dazhuqu massif; f-Pillow basalts of the Luqu massif3.2 大竹曲岩体
大竹曲岩体位于日喀则蛇绿岩东部，紧邻仁布岩体，距日喀则市区80 km左右（图 1b）。岩体北侧为古近—新近纪日喀则组，南侧为侏罗—白垩系深海沉积物，为断层接触关系。大竹曲岩体的地幔橄榄岩部分主要由方辉橄榄岩和和少量纯橄岩组成。壳幔过渡带的地幔橄榄岩强烈蛇纹石化且被辉绿岩和辉长岩脉侵入。而洋壳序列底部为镁铁质-超镁铁质堆积杂岩，包括橄长岩、辉石岩、层状辉长岩、块状辉长岩，但是难以观察到这些岩石的界限和接触关系（王希斌等，1987）。层状辉长岩具典型的暗示辉长岩和浅色辉长岩交替韵律层结构（图 2c）。暗色辉长岩一般含有60%~70%单斜辉石和30%~40%斜长石，而浅色辉长岩含50%~60%斜长石和40%~50%单斜辉石。斜长石多呈自形—半自形晶，中-粗粒结构，多已蚀变，单斜辉石多呈半自形粒状，并与斜长石组成典型的辉长结构（图 3a）。
![]()
图 3 日喀则蛇绿岩镁铁质岩石显微照片a—辉长岩中半自形单斜辉石和自形-半自形斜长石构成辉长结构；b—辉长岩中粗粒单斜辉石中包裹自形—半自形的斜长石晶体；c—辉绿岩中细粒辉石和自形斜长石组成典型的辉绿结构；d—辉绿岩中自形斜长石斑晶；e—玄武岩中自形斜长石和他形单斜辉石；f—玄武岩中自形斜长石和他形单斜辉石斑晶Figure 3. Petrographic microphotos of the mafic rocks from the Xigaze ophiolitea-Gabbro texture consisting of subhedral clinopyroxene and subhedral to euhedral plagioclase in gabbro; b-Subhedral-euhedral plagioclase enclosed by coarse-grain clinopyroxene in gabbro; c-Diabasic texture consisting of fine grained clinopyroxene and euhedral plagioclase in diabase; d-Euhedral plagioclase phenocrysts in diabase; e-Fine grained plagioclase and clinopyroxene in basalt; f-Euhedral plagioclase and anhedral clinopyroxene phenocrysts in basalt常见粗大的单斜辉石晶体包裹自形—半自形的斜长石晶体（图 3b）。块状辉长岩常显示块状构造，中粒结构，岩石常被辉绿岩顺层侵入切割（图 2d）。块状辉长岩主要组成为半自形的单斜辉石和已蚀变的自形-半自形的斜长石。大竹曲辉绿岩多为顺层侵入的辉绿岩席，且这些辉绿岩席也普遍侵入到下部的地幔橄榄岩中。部分侵入地幔橄榄岩的辉绿岩发生了异剥钙榴岩化。这些辉绿岩一般粒度较细，呈典型的辉绿结构（图 3c），辉绿岩中常出现自形—半自形的斜长石斑晶（图 2e，图 3d）。在辉绿岩席的上部为一定厚度的玄武岩，这些玄武质喷出岩粒度变化较大（图 3e~f），部分玄武岩出现斑晶，以自形—半自形斜长石为主，单斜辉石较少且自形程度较低（图 3f）。
4. 分析方法
4.1 全岩主量和微量元素分析
全岩主量元素分析在澳实分析检测（广州）有限公司完成，分析仪器为荷兰帕纳科公司生产的X荧光光谱仪，型号为PANalytical MagXFast。通过XRF射线光谱仪对制成四硼酸锂玻璃进行分析，分析精度优于5%。分析方法为：首先将已扣除烧失量的准确称量的粉末样品和Li2B4O7-LiBO2助熔剂，混合均匀，转移至铂金坩埚中，在1200℃下使样品完全熔融，再将熔体倒入铂金模具，冷却后制成玻璃再上机测试。全岩氧化亚铁（FeO）含量分析在澳实分析检测（广州）有限公司完成，根据强酸消解、重铬酸钾滴定测定样品中Fe2+含量。
全岩微量元素测试分析在中国科学院地球化学研究所矿床地球化学国家重点实验室完成，测试仪器为Perkin-Elmer Sciex ELAN DRC-e电感耦合等离子体质谱仪（ICP-MS），分析精度优于10%。准确称取50 mg粉末（200目）样品放入Teflon杯中，加入1 mL亚沸HF，并在电热板上蒸干以去除大部分SiO2，再加入1 mL亚沸HF和0.5 mL亚沸HNO3，然后用不锈钢钢套密封后放入烘箱中，温度185℃加热36 h。冷却至室温后在电热板上150℃蒸干去除HF，然后再加入1 mL亚沸HNO3再次蒸干，并重复一次全部除去HF。最后加入2 mL HNO3、5 mL去离子水和500 ng Rh内标溶液，重新放入不锈钢钢套密封放入烘箱，温度140℃加热8 h，冷却至室温后吸取其中400 μL溶液移到10 mL离心管，加入去离子水定容至10 mL，最后上机测试。
4.2 全岩Pb同位素分析
Pb同位素分析在中国地质科学院地质研究所国土资源部同位素地质研究重点实验室完成。化学流程如下：吸取样品溶液0.2 mL转移至AG1-X8阴离子交换树脂，并用1 mol/L HBr溶液淋洗5次，每次1 mL，然后用2mol/L HCl溶液0.5 mL淋洗1次，然后再用6 mol/L HCl溶液1 mL淋洗1次并接Pb，最后继续用6 mol/L HCl溶液0.5 mL淋洗1次并接Pb。Pb同位素比值用多接收器等离子体质谱法（MC-ICPMS）测定，所用仪器为英国Nu Plasma HR，仪器的质量分馏以Tl同位素外标校正。样品中Tl的加入量约为铅含量的1/2。国际标样NBS 981长期测定的统计结果：208Pb/206Pb =2.16736 ± 0.00066 (± 2σ)，207Pb/206Pb =0.91488 ± 0.00028 (± 2σ)，206Pb/204Pb=16.9386 ± 0.0131 (± 2σ)，207Pb/204Pb=15.4968±0.0107 (±2σ)，208Pb/204Pb= 36.7119±0.0331 (± 2σ)。详细的化学分离方法和质谱测定方法见何学贤等（2007）。
5. 日喀则蛇绿岩镁铁质岩石地球化学特征
5.1 全岩主量元素特征
日喀则蛇绿岩路曲和大竹曲岩石样品主微量元素分析结果见表 1。路曲和大竹曲岩体的玄武岩、辉绿岩、辉长岩均有一定程度的蚀变，烧失量分别为2.28%~3.20%、2.32%~3.42 %、3.06%~3.24%。路曲镁铁质岩石的SiO2含量在49.38%~50.62%，Na2O含量为2.73% ~4.35%，K2O含量为0.25% ~1.14%。TAS投图显示路曲样品均落在玄武岩/辉长岩区域（图 4）。大竹曲岩体镁铁质岩石的SiO2含量为48.89%~51.17%，Na2O含量为2.38%~4.71%，K2O含量为0.07%~1.32 %。TAS投图显示大竹曲样品也全部落在玄武岩/辉长岩范围之内（图 4）。
表 1 日喀则蛇绿岩路曲和大竹曲岩体镁铁质岩石全岩主量元素(%)和微量元素(10-6)成分Table 1. The compositions of whole-rock major (%) and trace (10-6) elements of the mafic rocks from the Luqu and Dazhuqu massifs of the Xigaze ophiolite
![]()
图 4 日喀则蛇绿岩路曲和大竹曲镁铁质岩石全岩TAS图解Figure 4. Diagram of TAS for whole rock composition of the mafic rocks of the Luqu and Dazhuqu massifs from the Xigaze ophiolite5.2 全岩微量元素特征
路曲和大竹曲岩石样品的球粒陨石标准化REE配分模式整体上显示LREE轻微亏损和HREE相对平坦，其特征与N-MORB非常相似（图 5a、5c）。在球粒陨石标准化REE配分模式中，路曲所有样品均无明显的Eu异常（图 5a），大竹曲岩体除了辉长岩样品显示轻微Eu正异常外，玄武岩和辉绿岩样品均没有明显的Eu异常（图 5c）。由于Rb、Ba、Cs、U、K、Sr、P等元素容易受到后期海水和热液的影响发生迁移，而高场强元素Nb、Ta、Zr、Hf、Ti和REE元素在后期热液活动过程中相对不容易发生迁移，它们主要受地幔源区和岩浆演化过程控制（Pearce et al., 2014）。因此，本次研究不选取这些容易受后期影响的元素进行标准化作图。在NMORB标准化微量元素配分模式中，所有路曲样品具有相似的微量元素特征且均显示弱Nb-Ta负异常，而Th、Zr、Hf、Ti、REE等元素特征与N-MORB一致（图 5b）。在大竹曲岩体中，大部分样品的微量元素特征保持一致，均显示弱Nb-Ta负异常（图 5d）。大竹曲玄武岩和辉绿岩的配分模式总体上相似于N-MORB，辉长岩相对于N-MORB显示略低的微量元素含量，但配分图形状接近于N-MORB。
![]()
图 5 日喀则蛇绿岩路曲和大竹曲岩石球粒陨石标准化REE配分模式（a，c）及N-MORB标准化微量元素配分模式（b，d）（球粒陨石和N-MORB标准化值据Sun and McDonough, 1989; 用于对比的OIB、E-MORB、N-MORB数值据Sun and McDonough, 1989）Figure 5. Chlondrite-normalized REE patterns (a, c) and N-MORB normalized trace elements patterns (b, d)of Xigaze ophioliteNormalized values, and N-MORB, OIB, E-MORB values after Sun and McDonough (1989)5.3 全岩Pb同位素特征
路曲和大竹曲样品全岩Pb同位素比值结果见表 2。路曲和大竹曲岩体样品的Pb同位素比值变化范围相对较小，其206Pb/204Pb的变化范围为17.54~18.05，207Pb/204Pb变化范围为15.42~15.56，208Pb/204Pb变化范围为37.67~38.17。根据样品中Th、U、Pb含量对Pb同位素比值的时间校正计算（t=140 Ma），校正结果显示(206Pb/204Pb)t的变化范围为17.46~17.94，(207Pb/204Pb)t变化范围为15.42~15.56，(208Pb/204Pb)t变化范围为37.35~38.13。在(208Pb/204Pb)t-(206Pb/204Pb)t和(207Pb/204Pb)t-(206Pb/204Pb)t图解上（图 6a、6b），路曲和大竹曲样品Pb同位素比值接近于雅鲁藏布蛇绿岩带其他岩体镁铁质岩石，且完全落入印度洋MORB区域。这些特征表明这些镁铁质岩石的地幔源区具有印度洋MORB型同位素组成特征。而雅鲁藏布蛇绿岩镁铁质岩石的Pb同位素比值则明显区别于现代大洋沉积物（Ben Othman et al., 1989）。
表 2 日喀则蛇绿岩路曲和大竹曲镁铁质岩石Pb同位素成分Table 2. The compositions of Pb isotope of mafic rocks from the Luqu and Dazhuqu massifs of the Xigaze ophiolite
![]()
图 6 全岩Pb同位素比值图解a-(208Pb/204Pb)t-(206Pb/204Pb)t; b-(207Pb/204Pb)t-(206Pb/204Pb)t；矫正年龄为t=140 Ma。太平洋和北大西洋MORB和印度洋MORB范围据Mahoney et al.1998和Zhang et al. 2005。现代大洋沉积物(GLOSS, 据Ben Othman et al., 1989, 北半球参考线(NHRL)据Hart, 1984。用于对比的雅鲁藏布蛇绿岩的Pb同位素数据源自文献（牛晓露等, 2006; 钟立峰等, 2006; 徐德明等, 2007; 刘飞等, 2013; Liu et al., 2015）Figure 6. Plots of Pb isotope ratiosa-(208Pb/204Pb)tversus (206Pb/204Pb)t; b-(207Pb/204Pb)t versus (206Pb/204Pb)t at 140 Ma. Fields of Pacific and North Atlantic MORB, and India MORB after Mahoney et al. (1998) and Zhang et al. (2005). Field of GLOSS after Ben Othman et al., 1898 and NHRL after Hart (1984). The data of Pb isotope of the Yarlung-Zangbo ophiolite after Niu et al. (2006), Zhong et al. (2006), Xu et al. (2007), Liu et al.(2013, 2015)6. 讨论
6.1 地幔源区性质
由于单斜辉石相对富集HREE，因此单斜辉石分离结晶将会导致残余岩浆亏损HREE而相对富集LREE（Hauri et al., 1994）。然而，路曲和大竹曲岩石样品的球粒陨石标准化REE配分模式均显示轻微LREE亏损（图 5a，5c），说明岩浆就位之前在深部没有经历重要的单斜辉石分离结晶。此外，路曲和大竹曲镁铁质岩石的球粒陨石标准化REE配分模式没有明显Eu负异常，而Eu相容于斜长石（Aigner-Torres et al., 2007），这说明形成这些镁铁质岩石的岩浆在深部没有经历重要的斜长石分离结晶（5a，5c）。这些特征说明经地幔部分熔融产生而聚集的岩浆在上升过程中没有经历明显的结晶分异作用，因此它们的微量元素能够反演地幔源区性质和部分熔融程度。
路曲和大竹曲岩体的玄武岩、辉绿岩、辉长岩均显示相似的球粒陨石标准化REE配分模式和NMORB标准化微量元素配分模式（图 5），说明形成这些岩石的岩浆来自同一地幔源区。此外，这些岩石的REE与N-MORB的REE配分模式非常接近（图 5a，5c），且N-MORB标准化微量元素配分模式也非常平坦且接近于1（图 5b，5d），表明形成这些岩石的岩浆来自亏损地幔源区。路曲和大竹曲岩石样品的Pb同位素均落在印度洋MORB范围之内（图 6），一方面说明这些岩石起源于亏损地幔，另一方面暗示新特提斯地幔源区与印度洋地幔域有相似的地球化学特征（Mahoney et al., 1998；Xu and Castillo, 2004；Zhang et al., 2005；牛晓露等，2006）。尖晶石二辉橄榄岩部分熔融过程中熔融残留物和熔体之间不会发生中稀土（如Sm）和重稀土（如Yb）分异，但轻稀土（如La）和中稀土（如Sm）之间会随着地幔部分熔融程度的增高而发生明显分异（Aldanmaz et al., 2000）。据图 7所示，大竹曲和路曲岩体样品均落在尖晶石二辉橄榄岩熔融趋势线上，其部分熔融程度约在10%。
![]()
图 7 全岩La/Sm-Sm/Yb图解（据Aldanmaz et al., 2000）PM、N-MORB、E-MORB数据源自Sun and McDonough (1989)Figure 7. Diagram of La/Sm versus Sm/Yb (after Aldanmaz et al., 2000)PM, N-MORB, E-MORB data after Sun and McDonough (1989)6.2 形成构造环境
作为雅鲁藏布蛇绿岩带上重要的组成部分，日喀则蛇绿岩的形成环境一直是众多学者关注的焦点（吴福元等，2014）。从野外地质关系来看，路曲岩体洋壳部分缺失堆晶岩，表明路曲岩体洋壳形成过程中可能不存在岩浆房。而路曲辉绿岩大多以岩席的形式顺层侵入，甚至少量辉绿岩脉侵入到下部的地幔橄榄岩中（图 2b）。这些特征与慢速扩张脊特征一致（Purdy and Detrick, 1986）。大竹曲岩体的辉绿岩和玄武岩野外特征与路曲岩体类似，但是发育规模较小的层状辉长岩和块状辉长岩（图 2c，2d），表明大竹曲岩体洋壳形成过程中存在一定规模的岩浆房。层状辉长岩显示浅色辉长岩层和暗色辉长岩层的韵律交替结构（图 2c），表明岩浆房发生一定程度的分异。辉长岩中常见粗粒单斜辉石包裹自形斜长石晶体（图 3b），说明斜长石结晶明显早于单斜辉石。此外，大竹曲辉绿岩和玄武岩中常出现自形—半自形的斜长石斑晶（图 2d，2f），这些特征都说明斜长石较早结晶并被上升岩浆捕获到浅部。实验岩石学研究表明玄武质岩浆中水增加将降低斜长石结晶温度，从而推迟斜长石结晶时间，然而在无水环境下，斜长石则较早从玄武质岩浆中结晶（Sisson and Grove, 1993; Botcharnikov et al., 2008）。因此，可以推断形成大竹曲镁铁质岩石的岩浆是一种起源于大洋中脊扩张中心的无水玄武质岩浆。
传统上判断蛇绿岩形成于俯冲带（SSZ）环境主要是根据蛇绿岩石序列中镁铁质岩石地球化学特征而没有其他直接地质证据（Miyashiro et al., 1973; Pearce et al., 1984, 2014）。然而，从Zr-Zr/Y图解来看（图 8a），日喀则蛇绿岩路曲和大竹曲岩体大部分岩石覆盖于洋中脊玄武岩范围，明显区别于弧前玄武岩，但是不能与弧后盆地玄武岩区别。在(Nb/Yb)N-(Th/Yb)N相关图中（图 8b），路曲和大竹曲岩体岩石分布范围和趋势与洋中脊玄武岩一致，而弧后盆地玄武岩则显示明显更高(Nb/Yb)N和(Th/Yb)N比值。在(Nb/Th)N-(Th/La)N相关图中（图 8c），路曲和大竹曲样品明显区别于弧前和弧后盆地玄武岩，完全落入洋中脊玄武岩区域。这些特征再次表明路曲和大竹曲镁铁质岩体岩石形成于洋中脊环境。此外，路曲和大竹曲岩体岩石的球粒陨石标准化REE配分模式与N-MORB配分模式完全一致（图 5a，5c），也暗示这些岩石形成于洋中脊环境。在NMORB标准化微量元素配分模式图中（图 5b，5d），路曲和大竹曲镁铁质岩石微量元素特征除了相对于Th和LREE显示轻微Nb-Ta负异常，而Th、LREE、Zr、Hf、Ti总体上与N-MORB保持一致。
![]()
图 8 微量元素构造环境判别图a-Zr/Y-Zr（据Pearce, 2000）；b-Nb/Yb-Th/Yb（据Pearce, 2008）；c-Nb/Th-Th/La（据Godard et al., 2006）；弧后盆地、洋中脊扩张中心、弧前玄武岩数据源自于PETDB数据库（http://www.petdb.org）Figure 8. Discriminant diagrams of tectonic settinga-Zr/Y versus Zr (after Pearce, 2000), b-Nb/Yb versus Th/Yb(after Pearce, 2008); c-Nb/Th versus Th/La(after Godard et al., 2006). The data of back-arc basalt, mid-oceanic ridge basalt, fore-arc basalt after PETDB database虽然路曲和大竹曲样品普遍存在轻微Nb-Ta负异常，而Nb-Ta负异常也被众多学者认为是SSZ成因的主要印记（Hebert et al., 2012；李文霞等，2012）。然而，值得注意的是路曲和大竹曲岩石的Pb同位素结果均落入印度洋MORB地幔域，而这些岩石的Nb-Ta负异常可能与印度洋MORB地幔域的地球化学不均一性有关（Moores et al., 2000）。先前的研究表明，青藏高原晚震旦世以来的蛇绿岩镁铁质岩石的同位素特征与现代印度洋MORB的同位素地球化学特征一致，明显区别于太平洋和北大西洋洋中脊玄武岩，表明现代印度洋MORB地幔域继承了特提斯洋地幔原来的空间位置和地球化学特征（Mahoney et al., 1998; Xu and Castillo, 2004; Zhang et al., 2005）。然而，在地质演化期间，特提斯洋岩石圈地幔域曾发生过大陆地壳、古老俯冲带物质混入（Castillo, 1988；Klein et al., 1988；Le Roex et al., 1989；Storey et al., 1989；Rehkamper and Hofmann, 1997）。由于在俯冲过程中Th容易迁移而Nb-Ta相对不易迁移，因此大多数与岛弧有关的岩石显示较高的Th/Nb和Th/Ta比值。此外，大陆地壳也显示较高的Th/Nb比值（Pearce et al., 2008），因此若地幔源区交代古老俯冲带物质或大陆地壳物质，则交代地幔形成的岩浆可能显示Nb-Ta负异常。而路曲和大竹曲镁铁质岩石的Pb同位素结果均落入印度洋MORB范围之内，暗示这些镁铁质岩石的地幔源区可能属于与现代印度洋MORB地幔域相似的异常地幔。因此，路曲和大竹曲镁铁质岩石中轻微的Nb-Ta负异常可能是继承而来。
7. 结论
（1）日喀则蛇绿岩路曲和大竹曲岩体镁铁质岩石的母岩浆起源亏损地幔源区，其地幔源区部分熔融程度在10%左右。这些镁铁质岩石的地幔源区的地球化学特征与印度洋MORB地幔域相似。
（2）路曲和大竹曲岩体镁铁质岩石的野外地质特征、岩相学特征和地球化学特征均支持这些岩石起源于大洋中脊扩张中心。微量元素比值也接近洋中脊玄武岩且明显区别于弧后盆地玄武岩和弧前玄武岩。N-MORB标准化微量元素配分模式显示弱Nb-Ta负异常可能是由于其地幔源区交代了古老俯冲物质。
-
![]()
图 5 2020年长江经济带生态环境遥感评估图
Figure 5. Remote sensing assessment map of the ecological environment of the Yangtze River Economic Zone (2020)
![]()
图 6 2001年长江经济带生态环境遥感评估图
Figure 6. The remote sensing assessment map of the ecological environment of the Yangtze River Economic Zone (2001)
![]()
图 7 长江经济带生态环境遥感评估动态变化图（2001—2020年）
Figure 7. Ecological environment dynamic changes of the Yangtze River Economic Zone evaluated by remote sensing(2001-2020)
![]()
图 8 长江经济带生态环境等级变化面积及占比
Figure 8. The area and proportion of changes in the ecological environment of the Yangtze River Economic Zone
![]()
图 9 三峡库区及周边生态环境变化状况（2001—2020年）
Figure 9. Changes in the Ecological Environment of the Three Gorges Reservoir Area and surrounding areas (2001-2020)
![]()
![]()
图 11 南京市主城区、莫干山核心区卫星影像（a、c）及区域生态环境质量指数分布（b, d）
Figure 11. Satellite image map (a, c) and ecological environmental quality index map (b, d) of the main urban area of Nanjing and the core area of Mogan Mountain
表 3 2020年长江经济带生态环境遥感评估等级表
Table 3 Remote sensing assessment grade table of the ecological environment of the Yangtze River Economic Zone (2020)
下载: 导出CSV
表 4 2020年长江经济带各省市生态环境统计表
Table 4 Ecological environment statistics of each province and city in the Yangtze River Economic Zone (2020)
下载: 导出CSV
表 5 2001—2020年长江经济带生态环境等级变化及面积占比
Table 5 Changes in the ecological environment grade and area proportion of the Yangtze River Economic Zone from 2001 to 2020
下载: 导出CSV
表 6 2001年和2020年长江经济带各省市平均生态环境指数
Table 6 Average ecological environment index of each province and city in the Yangtze River Economic Zone in 2001 and 2020
下载: 导出CSV
-
Baig M H A, Zhang Lifu, Shuai Tong, Tong Qingxi. 2014. Derivation of a tasselled cap transformation based on Landsat 8 at-satellite reflectance[J]. Remote Sensing Letters, 5(5): 423-431. doi: 10.1080/2150704X.2014.915434
China Geological Survey, Ministry of Natural Resources. 2018. Analysis Report on Resource and Environmental Conditions and Major Issues Supporting the Land and Space Planning of the Yangtze River Economic Zone[R]. Beijing: China Geological Survey.
Chen Jiang, Wan Li, Liang Sihai, Jin Xiaomei, Chen Li. 2007. A tentative discussion on the trend of ecological environment change in Qinghai-Tibet Plateau[J]. Acta Geoscientia Sinica, 28(6): 555(in Chinese with English abstract). http://www.oalib.com/paper/1557546
Cheng Hangxin, Li Kuo, Li Min, Cheng Xiaomeng. 2015. Management target value (MTV) and rectification action value (RAV) of trace metals in urban soil in China[J]. Earth Science Frontiers, 22(5): 215-225 (in Chinese with English abstract). http://www.researchgate.net/publication/283129654_Management_target_value_MTV_and_rectification_action_value_RAV_of_trace_metals_in_urban_soil_in_China
Crist E P. 1985. A TM tasseled cap equivalent transformation for reflectance factor data[J]. Remote Sensing of Environment, 17(3): 301-306. doi: 10.1016/0034-4257(85)90102-6
Duan Yiren, Yang Zhongfang, Yang Qiong, Zheng Guodong, Zhuo Xiaoxiong, Chen Biao. 2020. The distribution of soil germanium and its influencing factors in Beibu Gulf of Guangxi[J]. Geology in China, 47(6): 1826-1837(in Chinese with English abstract).
Gong Peng, Li Xuecao, Zhang Wei. 2019.40-year (1978-2017) human settlement changes in China reflected by impervious surfaces from satellite remote sensing[J]. Chinese Science Bulletin, 64(11): 756-763. http://www.cnki.com.cn/Article/CJFDTotal-JXTW201911006.htm
Gorelick N, Hancher M, Dixon M, SimonIl Y, Thau D, Moore R. 2017. Google Earth Engine: Planetary-scale geospatial analysis for everyone[J]. Remote Sensing of Environment, 202. http://www.sciencedirect.com/science/article/pii/S0034425717302900
Goward S N, Xue Yongkang, Czajkowski K P. 2002. Evaluating land surface moisture conditions from the remotely sensed temperature/vegetation index measurements: An exploration with the simplified simple biosphere model[J]. Remote Sensing of Environment, 79(2/3): 225-242. http://www.remotesensing.utoledo.edu/resch/student/digimgclas/articles/evaluate.pdf
Gupta K, Kumar P, Pathan S K, Sharma K P. 2012. Urban Neighborhood Green Index-A measure of green spaces in urban areas[J]. Landscape and urban planning, 105(3): 325-335. doi: 10.1016/j.landurbplan.2012.01.003
Huang C, Wylie B, Yang L, Homer C, Zylstra G. 2002. Derivation of a tasselled cap transformation based on Landsat 7 at-satellite reflectance[J]. International Journal of Remote Sensing, 23(8): 1741-1748. doi: 10.1080/01431160110106113
Huang Jinliang, Pontius Jr R G, Li Qingsheng, Zhang Yujia. 2012. Use of intensity analysis to link patterns with processes of land change from 1986 to 2007 in a coastal watershed of southeast China[J]. Applied Geography, 34(3): 371-384. https://www.sciencedirect.com/science/article/abs/pii/S0143622812000033
Jiang Yuehua, Lin Liangjun, Chen Lide, Ni Huayong, Ge Weiya, Cheng Hangxin, Zhai Gangyi, Wang Guiling, Ban Yizhong, Li Yuan, Lei Mingtang, Tan Chengxuan, Su Jingwen, Zhou Quanping, Zhang Taili, Li Yun, Liu Hongying, Peng Ke, Wang Hanmei. 2017. Research on conditions of resources and environment and major geological problems in the Yangtze River Economic Zone[J]. Geology in China, 44(6): 1045-1061 (in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DIZI201706002.htm
Jiang Yuehua, Zhou Quanping, Chen Lide, Ni Huayong, Lei Mingtang, Cheng Heqin, Shi Bin, Ma Teng, Ge Weiya, Su Jingwen, Li Yun, Tan Jianmin. 2019. Progresses and main achievements of geological environment comprehensive survey project in the Yangtze River Economic Zone[J]. Geological Survey of China, 6(5): 1-20(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-DZDC201905001.htm
José A Sobrino, Juan C. Jiménez-Muoz, Paolini L. 2004. Land surface temperature retrieval from LANDSAT TM 5[J]. Remote Sensing of Environment, 90(4): 434-440. doi: 10.1016/j.rse.2004.02.003
Kerr J T, Ostrovsky M. 2003. From space to species: Ecological applications for remote sensing[J]. Trends in Ecology & Evolution, 18(6): 299-305. http://www.sciencedirect.com/science/article/pii/S0169534703000715
Liu Yingping, Chen Wende, Xu Wei, Hao Pengpei. 2012. Soil quality geochemical evaluation of Chengdu Economic Area[J]. Geoscience, 26(5): 865(in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-XDDZ201205003.htm
NASA. 2012. Landsat 7 Science Data Users Handbook[EB/OL]. http://landsathandbook.gsfc.nasa.gov,2012-09-05.
Nichol J. 2005. Remote sensing of urban heat islands by day and night[J]. Photogrammetric Engineering & Remote Sensing, 71(5): 613-622. http://www.ingentaconnect.com/contentone/asprs/pers/2005/00000071/00000005/art00004
Shabiti M, Yusufu D. 2008. Analysis on the Characteristics of Landuse Space-Time Change in the Economic Belt on the Northern Slope of Tianshan Mountain, Xinjiang, China[J]. Areal Research and Development, 27(2): 103-108(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DYYY200802027.htm
Shawuti R, Abuli A, Li H. 2003. Dynamic monitoring and analysis of ecological environment in Fukang City based on RSEI model[J]. Research of Soil and Water Conservation, 27(1): 289-295, 303(in Chinese with English abstract). https://www.ncbi.nlm.nih.gov/pubmed/29704350
State Environmental Protection Administration. 2006. Environmental Protection Industry Standard of the People's Republic of China (Trial) HJ/T192-2006[S].
Sun Yangui, Zhang Guowei, Wang Dongqing, Zhang Xianting, Li Dongling. 2003. Application of the remote sensing technique in eco-environmental division in Qinghai Province[J]. Geology in China, 30(2): 214-219(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DIZI200302015.htm
Sobrino J A, Jiménez-Muñoz J C, Paolini L. 2004. Land surface temperature retrieval from LANDSAT TM 5[J]. Remote Sensing of Environment, 90(4): 434-440. doi: 10.1016/j.rse.2004.02.003
Teng Mingjun, Zeng Lixiong, Xiao Wenfa, Zhou Zhixiang, Huang Zhilin, Wang Pengcheng, Dian Yuanyong. 2014. Research progress on remote sensing of ecological and environmental changes in the Three Gorges Reservoir area, China[J]. Chinese Journal of Applied Ecology, 25(12): 3683-3693(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-YYSB201412040.htm
Wang Qinhua, Dong Yanxiang, Song Mingyi, Liu Junbao, Huang Chunlei. 2011. Research on the integration method of land quality geochemical evaluation and agricultural land classification: Taking Jiashan County and Cixi City in Zhejiang Province as examples[J]. Shanghai Land and Resources, 32(4): 20-25 (in Chinese)
Wang Xuxi. 2021. Urban ecosystem health comprehensive assessment of cities on Yangtze River economic belt[J]. Journal of Mianyang Teachers' College, 40(2): 105-110(in Chinese with English abstract). https://www.sciencedirect.com/science/article/abs/pii/S0048969720383704
Wu Zhonghai, Zhou Chunjing, Tan Chengxuan, Sun Yujun, Ma Xiaoxue. 2016. The active tectonics and regional crustal stability features in the area of Yangtze River economic belt[J]. Journal of Geomechanics, 22(3): 379-411(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZLX201603001.htm
Xie Xuejin, Cheng Hangxin, Xie Yuanru. 2002. Analytical methods and quality in the compilation of 76 elements geochemical atlas of Sichuan, Yunnan, Guizhou, and Guangxi Provinces of China (1): Similarity of geochemical maps compiled from data generated by different laboratories-examples for Ag, Cs, Ga, and Ge analysis[J]. Geological Bulletin of China, 21(6): 277-284(in Chinese with English abstract). http://www.cnki.com.cn/article/cjfdtotal-zqyd200206000.htm
Xu H. 2010. Analysis of impervious surface and its impact on urban heat environment using the normalized difference impervious surface index (NDISI)[J]. Photogrammetric Engineering & Remote Sensing, 76(5): 557-565. http://www.cabdirect.org/abstracts/20103149181.html
Xu Hanqiu. 2013a. A remote sensing index for assessment of regional ecological changes[J]. China Environmental Science, 33(5): 889-897(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZGHJ201305023.htm
Xu Hanqiu. 2013b. A remote sensing urban ecological index and its application[J]. Acta Ecologica Sinica, 33(24): 7853-7862(in Chinese with English abstract). https://www.hindawi.com/journals/js/2017/1353691/
Yang Huiting, Xu Hanqiu. 2020. Assessing fractional vegetation cover changes and ecological quality of the Wuyi Mountain National Nature Reserve based on remote sensing spatial information[J]. The Journal of Applied Ecology, 31(2): 533-542(in Chinese with English abstract).
Zhao Yan, Luan Wenlou, Guo Haiquan, Cai Kui, Ma Zhongshe, Dun Yanran. 2021. Characteristics, causes and ecological environment health evaluation of Selenium-enriched soil in Gaocheng District of Shijiazhuang City, Hebei Province[J]. Geology in China, 48(3): 764-776(in Chinese with English abstract).
Zhao Yuelong, Zhang Lingjuan. 1999. Study on method of quantitative assessment of fragile environment[J]. Chinese Geographical Science, 9(2): 141-145(in Chinese). doi: 10.1007/BF02791364
陈江, 万力, 梁四海, 陈立. 2007. 青藏高原生态环境变化趋势的初步探索[J]. 地球学报, 28(6): 555-560. doi: 10.3321/j.issn:1006-3021.2007.06.007 成杭新, 李括, 李敏, 成晓梦. 2015. 中国城市土壤微量金属元素的管理目标值和整治行动值[J]. 地学前缘, 22(5): 215-225. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201505021.htm 段轶仁, 杨忠芳, 杨琼, 郑国东, 卓小雄, 陈彪. 2020. 广西北部湾地区土壤锗分布特征及其影响因素及其生态环境评价[J]. 中国地质, 47(6): 1826-1837. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20200618&flag=1 国家环保总局. 2006. 中华人民共和国环境保护行业标准(试行) HJ/T192-2006[S]. 姜月华, 林良俊, 陈立德, 倪化勇, 葛伟亚, 成杭新, 翟刚毅, 王贵玲, 班宜忠, 李媛. 2017. 长江经济带资源环境条件与重大地质问题[J]. 中国地质, 44(6): 1045-1061. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20170601&flag=1 姜月华, 周权平, 陈立德, 倪化勇, 雷明堂, 程和琴, 施斌, 马腾, 葛伟亚, 苏晶文. 2019. 长江经济带地质环境综合调查工程进展与主要成果[J]. 中国地质调查, 6(5): 1-20. https://www.cnki.com.cn/Article/CJFDTOTAL-DZDC201905001.htm 刘应平, 陈文德, 许伟, 彭培好. 2012. 成都经济区土壤质量地球化学评估[J]. 现代地质, (5): 865-872. doi: 10.3969/j.issn.1000-8527.2012.05.004 满苏尔·沙比提, 迪里木拉提·玉苏甫. 2008. 新疆天山北坡经济带土地利用时空变化特征研究[J]. 地域研究与开发, 27(2): 103-108. doi: 10.3969/j.issn.1003-2363.2008.02.026 茹克亚·萨吾提, 阿不都艾尼·阿不里, 李虎, 尼加提·卡斯木, 李晓航. 2020. 基于遥感生态指数模型的阜康市生态环境动态变化监测与评价[J]. 水土保持研究, 27 (1): 289-295, 303. https://www.cnki.com.cn/Article/CJFDTOTAL-STBY202001041.htm 孙延贵. 2003. 青海省生态环境分区的遥感应用研究[J]. 中国地质, 30(2): 214-219. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20030216&flag=1 滕明君, 曾立雄, 肖文发, 周志翔, 黄志霖, 王鹏程, 佃袁勇. 2014. 长江三峡库区生态环境变化遥感研究进展[J]. 应用生态学报, (12): 3683-3693. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201412040.htm 汪庆华, 董岩翔, 宋明义, 刘军保, 黄春雷. 2011. 土地质量地球化学评估与农用地分等成果整合方法研究——以浙江嘉善县和慈溪市为例[J]. 上海国土资源, 32(4): 20-25. doi: 10.3969/j.issn.2095-1329.2011.04.005 王旭熙. 2021. 长江经济带城市生态系统健康评价[J]. 绵阳师范学院学报, 40(2): 105-110. https://www.cnki.com.cn/Article/CJFDTOTAL-MYSF202102020.htm 吴中海, 周春景, 谭成轩, 孙玉军, 马晓雪. 2016. 长江经济带地区活动构造与区域地壳稳定性基本特征[J]. 地质力学学报, 22(3): 379-411, 372. doi: 10.3969/j.issn.1006-6616.2016.03.001 谢学锦, 成杭新, 谢渊如. 2002. 川滇黔桂76种元素地球化学图编制中分析方法与分析质量研究(一)不同实验室产生地球化学图的相似性——以Ag, Cs, Ga, Ge为例[J]. 地质通报, 21(6): 2-9. https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200206000.htm 徐涵秋. 2013a. 区域生态环境变化的遥感评价指数[J]. 中国环境科学, 33(5): 889-897. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGHJ201305023.htm 徐涵秋. 2013b. 城市遥感生态指数的创建及其应用[J]. 生态学报, 33(24): 7853-7862. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201324027.htm 杨绘婷, 徐涵秋. 2020. 基于遥感空间信息的武夷山国家级自然保护区植被覆盖度变化与生态质量评估[J]. 应用生态学报, 31(2): 187-196. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB202002023.htm 赵燕, 栾文楼, 郭海全, 蔡奎, 马忠社, 敦妍冉. 2021. 河北省石家庄市藁城区富硒土壤特征、成因与生态环境健康评价[J]. 中国地质, 48(3): 764-776. http://geochina.cgs.gov.cn/geochina/ch/reader/view_abstract.aspx?file_no=20210307&flag=1 赵跃龙. 1998. 脆弱生态环境定量评价方法的研究[J]. 地理科学进展, 18(1): 73-79. https://www.cnki.com.cn/Article/CJFDTOTAL-DLKJ801.009.htm
计量
- 文章访问数: 3022
- HTML全文浏览量: 848
- PDF下载量: 4289
