基于GIS的长江口海域生态系统脆弱性综合评价

ＡＣＴＡＥＣＯＬＯＧＩＣＡＳＩＮＩＣＡ

生态学报

Ｖｏｌ．３９，Ｎｏ．１１Ｊｕｎ．，２０１９

ＤＯＩ：１０．５８４６／ｓｔｘｂ２０１８０７０９１４９１

何彦龙，袁一鸣，王腾，张昊飞，陈耀辉．基于ＧＩＳ的长江口海域生态系统脆弱性综合评价．生态学报，２０１９，３９（１１）：３９１８⁃３９２５．ＥｃｏｌｏｇｉｃａＳｉｎｉｃａ，２０１９，３９（１１）：３９１８⁃３９２５．

ＨｅＹＬ，ＹｕａｎＹＭ，ＷａｎｇＴ，ＺｈａｎｇＨＦ，ＣｈｅｎＹＨ．ＩｎｔｅｇｒａｔｅｄａｓｓｅｓｓｍｅｎｔｏｆｍａｒｉｎｅｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙｕｓｉｎｇＧＩＳ．Ａｃｔａ

基于ＧＩＳ的长江口海域生态系统脆弱性综合评价

何彦龙１，２，∗，袁一鸣１，２，王　腾１，２，张昊飞１，２，陈耀辉２，３

１国家海洋局东海环境监测中心，上海　２０１２０６

２国家海洋局海洋赤潮灾害立体监测技术与应用重点实验室，上海　２０１２０６３上海海洋大学，海洋生态与环境学院，上海　２０１３０６

摘要：气候变化、富营养化、生境破碎等是全球普遍面临的生态问题，科学评估生态系统外部压力及其弹性力，对生态系统管理和生态修复具有重要的指导意义。使用空间主成分分析（ＳＰＣＡ）和层次分析法（ＡＨＰ）构建评价指标体系，结合地理信息系统软件，对长江口海域生态环境脆弱性进行综合评价，并根据生态环境脆弱性指数（ＥＶＩ）值，将研究区生态环境脆弱性分为５级：微度脆弱（＜０．５）、轻度脆弱（０．５—０．８）、中度脆弱（０．８—１．０）、重度脆弱（１．０—１．２）、极度脆弱（１．２—１．５）。结果表明：空间尺度上，长江口口门内生态环境脆弱度最高，生态环境脆弱度从口门内向口门外呈显著的降低趋势，近五年，口门内极度脆弱区空间分布南移；评估区域内，约２０００ｋｍ２的极度脆弱区发生了转变，极度脆弱区、重度脆弱区面积占比分别下降了７％和５％，长江口海域生态环境脆弱性明显好转。总体上，近年来大量陆源污染物输入以及生态系统结构变化，是导致长江口生态环境脆弱度较高的重要因素。

关键词：生态系统；综合评价；脆弱性；地理信息系统；长江口

ＩｎｔｅｇｒａｔｅｄａｓｓｅｓｓｍｅｎｔｏｆｍａｒｉｎｅｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙｕｓｉｎｇＧＩＳ

１ＥａｓｔＣｈｉｎａＳｅａＥｎｖｉｒｏｎｍｅｎｔａｌＭｏｎｉｔｏｒｉｎｇＣｅｎｔｅｒ，ＳｔａｔｅＯｃｅａｎｉｃＡｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｐｅｏｐｌｅ′ｓＲｅｐｕｂｌｉｃｏｆＣｈｉｎａ，Ｓｈａｎｇｈａｉ２０１２０６，ＣｈｉｎａＣｈｉｎａ，Ｓｈａｎｇｈａｉ２０１２０６，Ｃｈｉｎａ

ＨＥＹａｎｌｏｎｇ１，２，∗，ＹＵＡＮＹｉｍｉｎｇ１，２，ＷＡＮＧＴｅｎｇ１，２，ＺＨＡＮＧＨａｏｆｅｉ１，２，ＣＨＥＮＹａｏｈｕｉ２，３

２ＫｅｙＬａｂｏｒａｔｏｒｙｏｆＩｎｔｅｇｒａｔｅｄＭｏｎｉｔｏｒｉｎｇａｎｄＡｐｐｌｉｅｄＴｅｃｈｎｏｌｏｇｙｆｏｒＭａｒｉｎｅＨａｒｍｆｕｌＡｌｇａｌＢｌｏｏｍｓ，ＳｔａｔｅＯｃｅａｎｉｃＡｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｐｅｏｐｌｅ′ｓＲｅｐｕｂｌｉｃｏｆ３ＣｏｌｌｅｇｅｏｆＭａｒｉｎｅＥｃｏｌｏｇｙａｎｄＥｎｖｉｒｏｎｍｅｎｔ，ＳｈａｎｇｈａｉＯｃｅａｎＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ２０１３０６，Ｃｈｉｎａ

Ａｂｓｔｒａｃｔ：Ｃｌｉｍａｔｅｃｈａｎｇｅ，ｅｕｔｒｏｐｈｉｃａｔｉｏｎ，ａｎｄｈａｂｉｔａｔｆｒａｇｍｅｎｔａｔｉｏｎａｒｅｕｂｉｑｕｉｔｏｕｓｅｃｏｌｏｇｉｃａｌｐｒｏｂｌｅｍｓｉｎｔｈｅｗｏｒｌｄ．Ｉｔｉｓｉｍｐｏｒｔａｎｔｔｏｃｏｎｄｕｃｔｓｃｉｅｎｔｉｆｉｃａｓｓｅｓｓｍｅｎｔｏｆｔｈｅｅｘｔｅｒｎａｌｐｒｅｓｓｕｒｅａｎｄｔｈｅｒｅｓｉｌｉｅｎｃｅｏｆｅｃｏｓｙｓｔｅｍｓｆｏｒｔｈｅｅｃｏｌｏｇｉｃａｌｍａｎａｇｅｍｅｎｔａｎｄｒｅｓｔｏｒａｔｉｏｎ．ＩｎｔｈｅｃａｓｅｏｆｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙ，ａｓｔｈｅｅｎｄｏｆｔｈｅｂａｓｉｎ，ｈｕｇｅａｍｏｕｎｔｓｏｆｐｏｌｌｕｔａｎｔｓｗｅｒｅｃａｒｒｉｅｄｉｎｔｏｔｈｅａｒｅａ，ｗｈｉｃｈｒｅｓｕｌｔｅｄｉｎｔｈｅｄｅｇｅｎｅｒａｔｉｏｎｏｆｔｈｅｅｃｏｓｙｓｔｅｍｄｕｒｉｎｇｔｈｅｐａｓｔｄｅｃａｄｅｓ．Ｈｏｗｅｖｅｒ，ｌｉｔｔｌｅｉｓｋｎｏｗｎａｂｏｕｔｔｈｅｔｅｎｄｅｎｃｙｏｆｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｔｈｉｓａｒｅａｕｎｄｅｒｔｈｅｅｆｆｅｃｔｓｏｆｈｕｍａｎａｃｔｉｖｉｔｉｅｓａｎｄｃｌｉｍａｔｅｃｈａｎｇｅ．Ｃｏｎｓｅｑｕｅｎｔｌｙ，ｔｈｅｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙｗａｓａｓｓｅｓｓｅｄｕｓｉｎｇｇｅｏｇｒａｐｈｉｃｉｎｆｏｒｍａｔｉｏｎｓｙｓｔｅｍ（ＧＩＳ）ｓｏｆｔｗａｒｅｃｏｍｂｉｎｅｄｗｉｔｈｓｐａｔｉａｌｐｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔａｎａｌｙｓｉｓ（ＳＰＣＡ）．Ａｎａｎａｌｙｔｉｃａｌｈｉｅｒａｒｃｈｙｐｒｏｃｅｓｓ（ＡＨＰ）ｗａｓｕｓｅｄｔｏｄｅｔｅｒｍｉｎｅｔｈｅｅｖａｌｕａｔｉｏｎｆａｃｔｏｒｗｅｉｇｈｔｉｎｇｓ．Ａｃｃｏｒｄｉｎｇｔｏｔｈｅｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｄｅｘ（ＥＶＩ）ｖａｌｕｅｓ，ｔｈｅｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｉｅｓｏｆｔｈｅｓｔｕｄｙａｒｅａｗｅｒｅｃｌａｓｓｉｆｉｅｄｉｎｔｏｆｉｖｅｌｅｖｅｌｓ：ｐｏｔｅｎｔｉａｌ（＜０．５），ｌｉｇｈｔ（０．５—０．８），ｍｏｄｅｒａｔｅ（０．８—１．０），ｈｅａｖｙ（１．０—１．２）ａｎｄｖｅｒｙｈｅａｖｙ（１．２—１．５）．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅｉｎｎｅｒｍｏｕｔｈａｒｅａｏｆｔｈｅ

基金项目：海洋公益性行业科研专项经费项目（２０１５０５００８）收稿日期：２０１８⁃０７⁃０９；　　网络出版日期：２０１９⁃０３⁃２１∗通讯作者Ｃｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒ．Ｅ⁃ｍａｉｌ：ｙａｎｌｈｅｌ＠１６３．ｃｏｍ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１１期　　　何彦龙　等：基于ＧＩＳ的长江口海域生态系统脆弱性综合评价　３９１９

ＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙｗａｓｔｈｅｍｏｓｔｖｕｌｎｅｒａｂｌｅ．ＴｈｅｒｅｗａｓａｓｉｇｎｉｆｉｃａｎｔｄｅｃｒｅａｓｅｔｅｎｄｅｎｃｙｉｎＥＶＩｆｒｏｍｔｈｅｉｎｓｉｄｅｔｏｔｈｅｏｕｔｓｉｄｅｏｆｔｈｅｍｏｕｔｈ．Ｔｈｅｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｖｅｒｙｈｅａｖｙｖｕｌｎｅｒａｂｌｅａｒｅａｓｗｉｔｈｉｎｔｈｅｍｏｕｔｈｈａｓｓｈｉｆｔｅｄｓｏｕｔｈｗａｒｄｉｎｖｕｌｎｅｒａｂｌｅａｒｅａｄｅｃｒｅａｓｅｄｂｙ７％ａｎｄｓｅｖｅｒｅｌｙｖｕｌｎｅｒａｂｌｅａｒｅａｄｅｃｒｅａｓｅｄｂｙ５％．Ｏｎｅｔｈｅｏｔｈｅｒｈａｎｄ，ｔｈｅｅｃｏｌｏｇｉｃａｌｔｈｅｌａｓｔｆｉｖｅｙｅａｒｓ．Ｔｈｅｅｘｔｒｅｍｅｌｙｖｕｌｎｅｒａｂｌｅａｒｅａｏｆａｂｏｕｔ２０００ｋｍ２ｈａｓｃｈａｎｇｅｄｉｎｔｈｅａｓｓｅｓｓｍｅｎｔａｒｅａｗｉｔｈｅｘｔｒｅｍｅｌｙｅｎｖｉｒｏｎｍｅｎｔｓｔａｔｕｓｏｕｔｏｆｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙｈａｓｏｂｖｉｏｕｓｌｙｉｍｐｒｏｖｅｄｆｒｏｍ２０１３—２０１７．Ｉｎｇｅｎｅｒａｌ，ｔｈｅｉｎｐｕｔｏｆｌａｒｇｅａｍｏｕｎｔｏｆｌａｎｄ⁃ｂａｓｅｄｐｏｌｌｕｔａｎｔｓａｎｄｔｈｅｉｎｓｔａｂｉｌｉｔｙｏｆｔｈｅｅｃｏｓｙｓｔｅｍｃａｕｓｅｄｂｙｈｕｍａｎａｃｔｉｖｉｔｙｗｅｒｅｔｈｅｍａｉｎｆａｃｔｏｒｓｆｏｒｔｈｅｈｉｇｈｖｕｌｎｅｒａｂｉｌｉｔｙｏｆｔｈｅＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙ．

ＫｅｙＷｏｒｄｓ：ｅｃｏｓｙｓｔｅｍ；ｉｎｔｅｇｒａｔｅｄａｓｓｅｓｓｍｅｎｔ；ｖｕｌｎｅｒａｂｉｌｉｔｙ；ＧＩＳ；ＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙ

健康的生态系统取决于生态系统的物质输入、能量流动和信息传递能否保持平衡，研究表明生态系统的动态平衡主要由其弹性力维持［１］。近年来，由于经济开发等人类干扰活动，输入到生态系统的物质、能量等已超出生态系统自恢复的能力，导致生态系统的脆弱性显著增加［２］。对区域生态环境脆弱性进行评估，是国内外学者可持续发展研究普遍采用的方法［３］。

至今对脆弱性定义仍然存在一定争议，但脆弱性被普遍用于表征生态系统对外部压力的响应。如ＩＰＣＣ

定义的生态系统脆弱性：在外部压力下系统倾向受到不利影响，包括敏感性或易受危害、缺乏应对以及适应的能力［４］。总体上，生态脆弱性可从３个方面进行表征：１）相对于生态系统本身的外部压力，２）生态系统自身的敏感性，以及３）生态系统自适应能力或弹性力［５⁃６］。国内外对脆弱性评价方法开展了大量研究，如基于模糊理论、灰色评价模型、人工神经网络、层次分析方法等［７⁃１０］。然而，目前的研究方法普遍存在一些不足，首先，评估等级主要依靠专家打分；其次，大多数研究仅能定性地开展脆弱性评价；再次，基于个体的室内实验数据难以支撑生态系统水平上脆弱性评估需要［５⁃６，１１⁃１２］。

河口是海陆直接交汇的区域，长期以来一直是人类活动最为频繁的区域，生态系统往往较为脆弱。近４０

年来我国河口海域生态环境状况发生了显著变化，长江口及毗邻海域生态系统生物量和环境容量明显下降，重要生境如产卵场和索饵场发生变迁等［１３⁃１４］。然而，对于长江口海洋生态环境脆弱性状况仍不十分清楚，尤其在人类活动和全球气候变化加剧的情景下，长江口生态环境脆弱性变化趋势如何，导致其脆弱性的主导因素是什么。因此，本研究利用空间主成分分析方法建立长江口生态环境脆弱性评估指标体系，分析近五年长江口生态环境脆弱性的动态变化趋势，试图为长江口基于生态系统的管理提供依据。１　实验方法１．１　研究区概况

本研究以长江口（范围３１°４８．５４′Ｎ—３０°４５．３６′Ｎ，１２１°６．６６′Ｅ—１２３°Ｅ）近岸海域为研究区域（图１）。长横沙岛以及九段沙将河口划分为“三级分汊、四口分流”的格局［１３］。

近年来，长江口生态环境状况不容乐观，由于处于江口呈扇形分叉分布，在徐六泾以下，崇明岛、长兴岛、

流域末端，长江径流不断将流域内的污染物输入该区域。目前，富营养化已成为该区域的主要环境污染问题，其污染物来源主要有３个方面，流域面源污染、市政排污、沿岸工业排污，研究表明，长江口废水排放总量约为３３５３万ｍ／ｄ

３

［１５］

图１　研究区域与站位分布图Ｆｉｇ．１　Ｔｈｅｓｔｕｄｙａｒｅａａｎｄｓａｍｐｌｉｎｇｓｉｔｅｓ

。《２０１７年中国生态环境状况公报》

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３９２０　生　态　学　报　　　３９卷　

１．２　数据来源及处理

最新结果显示，长江口海域生态环境状况仍然较差，环境污染问题较为突出［１６］。

研究所用数据均为２０１３年、２０１５年和２０１７年８月份长江口海域生态环境监测数据。由于不同指标的量级不同，在评价之前首先对数据进行标准化处理，标准化公式如下：

（１）当评价指标为正指标（２）当评价指标为负指标

ｐｉｊ＝

ｘｉｊ－ｘｍｉｎ

（１）（２）

ｐｉｊ＝

ｘｍａｘ－ｘｍｉｎｘｍａｘ－ｘｉｊ

式中，Ｐｉｊ为各指标标准化结果，Ｘｉｊ为实际数值，Ｘｍａｘ，Ｘｍｉｎ分别代表最大值和最小值。将标准化后的属性数据在１．３　指标体系构建方法

ＡｒｃＧＩＳ中利用普通克里金插值转换为空间栅格数据，得到指标的空间化数据。

ｘｍａｘ－ｘｍｉｎ

本研究依据脆弱性定义以及国内外研究的相关成果从生态暴露度、生态敏感性和生态适应能力３个方面，构建１１个指标作为初选指标，为减少评价指标间的相关性，避免由于指标重复导致评价结果的准确性降低，本研究采用主成分分析方法对评价指标进行筛选，空间主成分分析在ＡｒｃＧＩＳ中实现。空间主成分法分析是基于ＧＩＳ的统计方法，其优点是将众多指标参数通过正交变换去除关联指标，但仍能保持所提供大量原有信息的方法，从一定程度避免了指标选择的主观性［１７］。本研究通过主成分分析，去除对主成分贡献较小的指标石油类浓度、浮游动物密度及鱼卵密度，最终评价指标由８个指标构成（表１）。

表１　长江口海域生态环境脆弱性评估指标

Ｔａｂｌｅ１　ＴｈｅｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｄｅｘｉｎｔｈｅＹａｎｇｔｚｅＲｉｖｅｒｅｓｔｕａｒｙ

目标层Ｔａｒｇｅｔｌａｙｅｒ脆弱性指数Ｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｄｅｘ，ＥＶＩ

第一层Ｌｅｖｅｌ１生态暴露度

第二层Ｌｅｖｅｌ２环境质量

第三层Ｌｅｖｅｌ３化学需氧量（＋）溶解氧浓度（－）无机氮浓度（＋）悬浮物含量（＋）

生态敏感性生态适应能力

　　“＋”表示脆弱性正指标，“－”表示脆弱性负指标

生态系统结构抗干扰能力

浮游植物密度（－）底栖生物密度（－）叶绿素ａ浓度（－）活性磷酸盐浓度（＋）

５个主成分作为研究区域生态环境脆弱性的计算指标进行分析。前三主成分特征值见表２。１．４　指标权重的确定

一般地，在主成分分析过程中，选取方差累积贡献率大于９５％的主成分作为计算指标［１８］。本研究选取前

本研究采用层次分析法（ＡｎａｌｙｔｉｃＨｉｅｒａｒｃｈｙＰｒｏｃｅｓｓ，ＡＨＰ）确定评价指标的权重，首先根据前五个主成分方差贡献率确定各主成分之间的重要性，其次两两比较各主成分之间的关系，确定各主成分相对的得分值，然后在Ｙａａｈｐ软件中计算不同指标的权重，最后进行一致性检验，若检验结果为随机一致性比率（ＣＲ）小于０．１，１．５　脆弱性指数及分类方法

表明各指标权重计算结果合理［１９］。本研究前５个主成分的权重见表３。

ｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｄｅｘ，ＥＶＩ），计算公式：

根据Ｗａｎｇ等脆弱性生态环境定量评价方法［１０］，对研究海域每个像元计算脆弱性指数（Ｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌ

ＥＶＩ＝

ｗｉｆｉ∑ｉ＝１

ｎ

（３）

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１１期　　　何彦龙　等：基于ＧＩＳ的长江口海域生态系统脆弱性综合评价　３９２１

式中，ＥＶＩ为生态环境脆弱性指数，ｗｉ为评价指标的权重，ｆｉ为评价指标的等级。

表２　长江口海域生态环境脆弱性评估指标前三主成分特征值

Ｔａｂｌｅ２　Ｔｈｅｅｉｇｅｎｖａｌｕｅｓｏｆｔｈｅｆｉｒｓｔｔｈｒｅｅｐｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔｏｆｔｈｅｉｎｄｉｃｅｓ

评估指标

Ａｓｓｅｓｓｍｅｎｔｉｎｄｅｘ

底栖生物Ｂｅｎｔｈｏｎｉｃｏｒｇａｎｉｓｍ叶绿素ａＣｈｌｏｒｏｐｈｙｌｌａ

化学需氧量Ｃｈｅｍｉｃａｌｏｘｙｇｅｎｄｅｍａｎｄ溶解氧Ｄｉｓｓｏｌｖｅｄｏｘｙｇｅｎ无机氮Ｉｎｏｒｇａｎｉｃｎｉｔｒｏｇｅｎ浮游植物Ｐｈｙｔｏｐｌａｎｋｔｏｎ活性磷酸盐Ｒｅａｃｔｉｖｅｐｈｏｓｐｈａｔｅ悬浮物Ｓｕｓｐｅｎｄｅｄｓｕｂｓｔａｎｃｅ

２０１３０．２３

主成分分析结果ＴｈｅｒｅｓｕｌｔｓｏｆｐｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔａｎａｌｙｓｉｓＰＣＡ３０．１５０．１５０．０１０．２５－０．０９－０．０２－０．５７０．７３

ＰＣＡ１０．０５０．１６０．３８０．０４０．５２０．１００．４７０．４２

ＰＣＡ２２０１５

ＰＣＡ３０．０９－０．０８－０．４６０．３７０．２４０．２１０．２１－０．１１

ＰＣＡ１０．１９０．１１０．２２０．１３０．６４０．１３０．５５０．３０

２０１７

ＰＣＡ１０．１１０．２８０．０４０．２３０．７００．２７０．４５０．２８

ＰＣＡ２－０．４６－０．１１－０．３９０．４８－０．４２－０．１４０．２０

－０．３６０．２１０．０６－０．２８－０．１３－０．０４－０．２７－０．１６

－０．３４－０．０７０．３００．３００．０２０．０３－０．４７０．１９

ＰＣＡ２ＰＣＡ３０．２００．０４－０．３２０．８３－０．２００．１６０．１１－０．１３

表３　长江口海域生态环境脆弱性评估指标各主成分权重Ｔａｂｌｅ３　ＴｈｅｗｅｉｇｈｔｏｆＰｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔｏｆｔｈｅｉｎｄｉｃｅｓ

年Ｙｅａｒ２０１３２０１５２０１７

主成分Ｐｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔ

第１主成分ＰＣＡ１

０．７５５２０．８５９８０．７１５２

第２主成分ＰＣＡ２

０．１０９３０．０９１４０．１２０４

第３主成分ＰＣＡ３

０．０９１０．０３８０．０９３６

第４主成分ＰＣＡ４

０．０２６０．００７４０．０５４８

第５主成分ＰＣＡ５

０．０１８５０．００３４０．０１６

ｖｕｌｎｅｒａｂｉｌｉｔｙｇｒａｄｅｓ，ＲＩＥＶＧ），计算公式：

根据Ｌｉｕ等研究方法［６］，本研究计算了区域综合脆弱性等级指数（Ｒｅｇｉｏｎａｌｉｎｔｅｇｒａｔｅｄｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌ

ＲＩＥＶＧ＝

Ｇｉ∑ｉ

ｍ

×

ＡｉＳｉ

（４）

式中，ＲＩＥＶＧ为区域综合脆弱性等级，Ｇｉ为ｉ区域脆弱性评价等级，Ａｉ为ｉ脆弱性等级对应的面积，Ｓｉ为区域总面积。

本研究在ＡｒｃＧＩＳ软件中，利用Ｒｅｃｌａｓｓｉｆｙ函数，采用自然断点法（ＮａｔｕｒａｌＢｒｅａｋｓ（Ｊｅｎｋｓ））将脆弱性等级分为５级。２　结果与分析

２．１　长江口海域生态环境脆弱性分级

２０１３—２０１７年，长江口海域生态环境脆弱性指数（ＥＶＩ）均值在０．５—１．５之间，平均值为１．２，其中０．６、

１．１、１．３是频率分布较多的ＥＶＩ值。脆弱性等级分级为：微度脆弱（＜０．５）、轻度脆弱（０．５—０．８）、中度脆弱（０．８—１．０）、重度脆弱（１．０—１．２）、极度脆弱（１．２—１．５）。不同的脆弱度面积占比不同，２０１３—２０１７年，极度脆弱度区域面积减少，微度和中度脆弱度区域面积增加（见表４）。２．２　长江口海域生态环境脆弱性时空分布

脆弱区，尤其是长江口南支海域脆弱性普遍较高，面积约为２０３９ｋｍ２，占本次评价海域面积的１５％。微度脆弱区和轻度脆弱区均分布在长江口门外部海域。２０１５年，极度脆弱区、重度脆弱区仍主要分布于长江口南支海域，北支重度脆弱性分布区域减少；２０１７年，长江口海域极度脆弱区和重度脆弱区分布面积显著减少，北支海域逐渐转为中度脆弱区，南支海域极度和重度脆弱区普遍南移，主要分布于上海市南汇嘴两侧近岸海域，与操文颖等人的研究结果基本一致［１５］。总体上，长江口海域生态环境脆弱度从口内向口门外逐渐降低，呈显著

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

２０１３年，长江口海域生态环境脆弱度主要以极度脆弱和重度脆弱为主，空间分布上，长江口门内为极度

３９２２　生　态　学　报　　　３９卷　

的梯度变化（图２）。

表４　长江口海域生态环境脆弱度等级及分布面积占比

Ｔａｂｌｅ４　Ｔｈｅｔｏｔａｌａｒｅａａｎｄｔｈｅｐｒｏｐｏｒｔｉｏｎｏｆｄｉｆｆｅｒｅｎｔｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｌｅｖｅｌｉｎＹａｎｇｔｚｅＲｉｖｅｒｅｓｔｕａｒｙ

脆弱性等级

Ｖｕｌｎｅｒａｂｉｌｉｔｙｌｅｖｅｒ微度Ｐｏｔｅｎｔｉａｌ轻度Ｌｉｇｈｔ中度Ｍｅｄｉｕｍ重度Ｈｅａｖｙ极度Ｖｅｒｙｈｅａｖｙ

ＥＶＩ值ＥＶＩｖａｌｕｅ＜０．５０．５—０．８０．８—１．０１．０—１．２１．２—１．５

栅格数Ｒａｓｔｅｒｎｕｍｂｅｒ２４７９５２６６８０１７３８２１４２６５１２４７８

２０１３

占总面积比Ｐｅｒｃｅｎｔｏｆｔｏｔａｌａｒｅａ／％

１７２０２４２３１５

栅格数Ｒａｓｔｅｒｎｕｍｂｅｒ２７９８６１７８２９８３７９２３０８６２８６４２

２０１５

占总面积比Ｐｅｒｃｅｎｔｏｆｔｏｔａｌａｒｅａ／％

３６２９１０１１１３

栅格数Ｒａｓｔｅｒｎｕｍｂｅｒ２４２６９１６３４２２０１９２１９３４０２５１０７

２０１７

占总面积比Ｐｅｒｃｅｎｔｏｆｔｏｔａｌａｒｅａ／％

３２２４２１１６８

面积Ａｒｅａ／ｋｍ２２３２４２６６０３２４７３０３３２０３９

面积Ａｒｅａ／ｋｍ２４８３０３８９２１３７３１４６５１７３５

面积Ａｒｅａ／ｋｍ２４２３８３１５７２７９５２０６６１０２０

图２　长江口海域生态环境脆弱度空间分布

Ｆｉｇ．２　Ｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎＹａｎｇｔｚｅＲｉｖｅｒＥｓｔｕａｒｙ

ＥＶＩ：生态环境脆弱性指数，Ｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｄｅｘ

２．３　近５年长江口海域生态环境脆弱度变化趋势

对２０１３年至２０１７年近５年长江口海域生态环境脆弱度趋势分析表明，近五年，长江口海域极度和重度脆弱区面积占比显著减少，极度脆弱区、重度脆弱区面积占比分别下降了７％和５％。微度脆弱区面积增加趋势显著，相比于２０１３年，２０１５年和２０１７年分别增加了１９％和１５％。中度和轻度脆弱区面积变化不大（图３）。

对２０１３年至２０１７年长江口海域生态环境区域综合脆弱性等级（Ｒｅｇｉｏｎａｌｉｎｔｅｇｒａｔｅｄｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌ２０１５年及２０１７年，ＲＩＥＶＧ分别为２．９８、２．３５、２．４３。近五年区域脆弱性等级下降了１８％，极度脆弱区和重度脆弱区显著减少对区域总体脆弱性下降贡献较高，表明长江口生态环境状况好转。２．４　近５年长江口海域生态环境脆弱区转化矩阵分析

近五年长江口生态环境脆弱度存在显著变化，表现为极度脆弱区和重度脆弱区面积显著减少，２０１３年到

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

ｖｕｌｎｅｒａｂｉｌｉｔｙｇｒａｄｅｓ，ＲＩＥＶＧ）分析表明（图４），ＲＩＥＶＧ值２０１３年最大，２０１５年和２０１７年显著减小，２０１３年、

　１１期　　　何彦龙　等：基于ＧＩＳ的长江口海域生态系统脆弱性综合评价　３９２３

２０１５年，约６６１ｋｍ２极度脆弱区转变为其他类型的区域，约１５００ｋｍ２重度脆弱区转变为中度、轻度和微度脆弱区（表５）。２０１５年到２０１７年，长江口区域脆弱性进一步改善，约１２００ｋｍ２的极度脆弱区转变为中度和重度脆弱区（表６）。

图３　近五年长江口海域生态环境脆弱度变化趋势

　Ｆｉｇ．３　Ｔｈｅｃｈａｎｇｉｎｇｔｒｅｎｄｏｆｐｅｒｃｅｎｔａｇｅｓｏｆｖｕｌｎｅｒａｂｉｌｉｔｙｇｒａｄｅｓｉｎｔｈｅｌａｓｔｆｉｖｅｙｅａｒｓ

图４　近五年长江口海域生态环境综合脆弱性等级变化趋势Ｆｉｇ．４　

Ｔｈｅｃｈａｎｇｉｎｇｔｒｅｎｄｏｆｒｅｇｉｏｎａｌｉｎｔｅｇｒａｔｅｄｅｃｏ⁃

ｅｎｖｉｒｏｎｍｅｎｔａｌｇｒａｄｅ（ＲＩＥＶＧ）ｖａｌｕｅｓｉｎｔｈｅｐａｓｔｆｉｖｅｙｅａｒｓ

表５　２０１３—２０１５年长江口海域生态环境脆弱区转换矩阵

Ｔａｂｌｅ５　Ｔｈｅｃｏｎｖｅｒｓｉｏｎｍａｔｒｉｘｏｆｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｆｒｏｍ２０１３ｔｏ２０１５

２０１３

微度脆弱区Ｐｏｔｅｎｔｉａｌ轻度脆弱区Ｌｉｇｈｔ中度脆弱区Ｍｅｄｉｕｍ重度脆弱区Ｈｅａｖｙ极度脆弱区Ｖｅｒｙｈｅａｖｙ

２０１５

微度脆弱区Ｐｏｔｅｎｔｉａｌ１８３８．８１１６３．６１４９８．９２３５．５８８．７

轻度脆弱区

Ｌｉｇｈｔ４８０．１１４９２．５１５５２．９２２１．２１４１．１

中度脆弱区Ｍｅｄｉｕｍ

０．００．０１８８．９１０４７．７１３６．８

重度脆弱区Ｈｅａｖｙ

０．００．００．０１１７０．３２９４．４

极度脆弱区Ｖｅｒｙｈｅａｖｙ

０．００．００．０３５２．３１３７８．３

表６　２０１５—２０１７年长江口海域生态环境脆弱区转换矩阵

Ｔａｂｌｅ６　Ｔｈｅｃｏｎｖｅｒｓｉｏｎｍａｔｒｉｘｏｆｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｆｒｏｍ２０１５ｔｏ２０１７

２０１５

微度脆弱区Ｐｏｔｅｎｔｉａｌ轻度脆弱区Ｌｉｇｈｔ中度脆弱区Ｍｅｄｉｕｍ重度脆弱区Ｈｅａｖｙ极度脆弱区Ｖｅｒｙｈｅａｖｙ

２０１７

微度脆弱区Ｐｏｔｅｎｔｉａｌ２２３２．５１９１５．３７２．００．００．０

轻度脆弱区

Ｌｉｇｈｔ１６４７．９１０９０．４２３６．１１８２．６０．０

中度脆弱区Ｍｅｄｉｕｍ９２２．３８８０．３６８８．３２５４．６４４．０

重度脆弱区Ｈｅａｖｙ

０．００．０３１６．５５５３．９１１９５．７

极度脆弱区Ｖｅｒｙｈｅａｖｙ

０．００．０５８．４４７３．０４８８．６

３　讨论

生态系统的脆弱性具有定量和定性两方面的特征，并且是一个复杂的持续变化的过程［２０⁃２２］，近年来，生

态环境脆弱度被广泛应用于局域生态环境受损程度及主要威胁要素的评估［２３⁃２６］。目前脆弱性评价指标体系中，有关外部压力指标的选择更多考虑人口分布密度和经济发展总量［２７⁃２８］，并赋以较高的权重。但是，人类活动对生态系统的影响并非全是负面，往往也有正向的促进作用［１１］。对于海洋生态系统而言，海域空间并无

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３９２４　生　态　学　报　　　３９卷　

直接的人口密度压力，人类活动的干扰更多体现在陆源污染物汇入到海洋生态系统造成的间接影响。因此，对于海洋生态系统脆弱性评价，选择经济发展总量或人口密度等指标可能会放大外部压力，导致脆弱性评价结果受主观性影响较大。同时，大多数陆源污染通过多种传播途径（河流、大气等）最终汇聚到海洋生态系统本研究指标的构建主要从生态系统内部考虑，通过分析主要污染物质迁入对生态环境造成的影响，从海洋生态环境基本状况和生态系统结构入手，选择主要监测指标作为评估指标，一定程度上避免了过度考虑经济指标对评估的不利影响。

评估结果表明，长江口门内生态环境脆弱度最高，重度脆弱区和极度脆弱区均分布于长江口门内南支海域，而造成该区域脆弱性较高的主要因素是氮、磷等污染物浓度高所导致，这与王保栋、张勇等人的研究结果基本一致［２９⁃３０］。表明长江流域污染物输入是导致河口地区生态环境脆弱度高的重要原因。同时，由于近年来长江流域水沙输入的变化，以及全球气候变化下河口区湿地生态系统萎缩［３１］，降低了河口生态系统的生态服务功能。通过近５年数据分析发现，长江口门内生态系统的初级生产力较低，浮游生物物种密度年际间变幅较大，生态系统结构不稳定，表明该海域生态系统弹性力降低，也进一步加剧了该区域生态环境的脆弱性。

从近５年脆弱性变化趋势来看，长江口生态环境脆弱度具有空间变化的特点，脆弱性等级从口门内到口门外呈显著的梯度变化。长江口海域极度脆弱区和重度脆弱区面积显著减少，口门内大部分海域从极度脆弱区降为重度脆弱区，表明近５年长江口生态环境脆弱度得到改善。由于受暖流北上制约，口门内污染物扩散受到一定程度的抑制，长江口悬浮物和污染物分布主要局限于近岸一带（１２３°以西）［３０］，同时黄浦江径流携带的氮磷、以及上海近岸、杭州湾污染物共同作用，导致长江口南汇嘴海域脆弱度增加，近５年长江口南支极度脆弱区移向南汇嘴至杭州湾北岸海域。研究表明，从２０世纪８０年始，长江口从低度富营养化水平逐渐上升到高度富营养化水平，随之导致了长江口生态系统结构发生了一系列变化，加速了河口区域的脆弱度［２９］。因此，导致长江口生态环境脆弱度较高的主导因素是人类活动的干扰所致，加大对长江口海域重要生态系统的保护以及流域污染物减排是改善该区域脆弱度的重要手段。４　结语

总体上，生态脆弱性评估是认识区域生态环境问题的主要方法，通过ＧＩＳ技术结合空间主成分分析方法和层次分析法，可以避免过多依赖于专家打分而可能造成的主观影响。然而，在生态系统水平上，脆弱性评价仅能辨识关键的影响因素，但就外部压力如何影响生态系统结构与过程，影响程度如何等更为深层次的问题，目前脆弱性评估方法仍然难以解决。因此，将来脆弱性评估研究首先应进一步夯实脆弱性评估的理论基础，在此基础上，拓展脆弱性评估的方法；其次，进一步探讨不同尺度脆弱性定量评估的方法，使得脆弱性评估从定性向定量评估发展。另外，本研究仅是对短期内长江口生态系统脆弱性变化趋势的评估，也未考虑气候变化对生态系统造成的影响。因此，结合气候变化，应进一步开展长江口中长期生态脆弱性评估研究，进一步理解大河口区域人类活动与气候变化双重压力下生态系统运维机制及其主导因素。

参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：

［１］　ＨｏｌｌｉｎｇＣＳ．Ｒｅｓｉｌｉｅｎｃｅａｎｄｓｔａｂｉｌｉｔｙｏｆｅｃｏｌｏｇｉｃａｌｓｙｓｔｅｍｓ．ＡｎｎｕａｌＲｅｖｉｅｗｏｆＥｃｏｌｏｇｙａｎｄＳｙｓｔｅｍａｔｉｃｓ，１９７３，４（１）：１⁃２３．

Ｂｉｏｌｏｇｙ，２００５，１９（１）：８６⁃９７．

中，随着海洋潮汐及洋流的变化处于动态变化中［２］。从而表明海洋生态系统的脆弱性具有动态变化的特点。

［２］　ＺａｃｈａｒｉａｓＭ，ＧｒｅｇｒＥＪ．Ｓｅｎｓｉｔｉｖｉｔｙａｎｄｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｍａｒｉｎｅｅｎｖｉｒｏｎｍｅｎｔｓ：ａｎａｐｐｒｏａｃｈｔｏｉｄｅｎｔｉｆｙｉｎｇｖｕｌｎｅｒａｂｌｅｍａｒｉｎｅａｒｅａｓ．Ｃｏｎｓｅｒｖａｔｉｏｎ［３］　ＳｃｈｅｆｆｅｒＭ，ＣａｒｐｅｎｔｅｒＳ，ＦｏｌｅｙＪＡ，ＦｏｌｋｅＣ，ＷａｌｋｅｒＢ．Ｃａｔａｓｔｒｏｐｈｉｃｓｈｉｆｔｓｉｎｅｃｏｓｙｓｔｅｍｓ．Ｎａｔｕｒｅ，２００１，４１３（６８５６）：５９１⁃５９６．

［４］　ＩＰＣＣ，２０１４：Ｓｕｍｍａｒｙｆｏｒｐｏｌｉｃｙｍａｋｅｒｓ．Ｉｎ：ＣｌｉｍａｔｅＣｈａｎｇｅ２０１４：Ｉｍｐａｃｔｓ，Ａｄａｐｔａｔｉｏｎ，ａｎｄＶｕｌｎｅｒａｂｉｌｉｔｙ．ＰａｒｔＡ：ＧｌｏｂａｌａｎｄＳｅｃｔｏｒａｌＡｓｐｅｃｔｓ．

ＣｏｎｔｒｉｂｕｔｉｏｎｏｆＷｏｒｋｉｎｇＧｒｏｕｐＩＩｔｏｔｈｅＦｉｆｔｈＡｓｓｅｓｓｍｅｎｔＲｅｐｏｒｔｏｆｔｈｅＩｎｔｅｒｇｏｖｅｒｎｍｅｎｔａｌＰａｎｅｌｏｎＣｌｉｍａｔｅＣｈａｎｇｅ［Ｆｉｅｌｄ，Ｃ．Ｂ．，Ｖ．Ｒ．Ｂａｒｒｏｓ，Ｄ．Ｊ．Ｄｏｋｋｅｎ，Ｋ．Ｊ．Ｍａｃｈ，Ｍ．Ｄ．Ｍａｓｔｒａｎｄｒｅａ，Ｔ．Ｅ．Ｂｉｌｉｒ，Ｍ．Ｃｈａｔｔｅｒｊｅｅ，Ｋ．Ｌ．Ｅｂｉ，Ｙ．Ｏ．Ｅｓｔｒａｄａ，Ｒ．Ｃ．Ｇｅｎｏｖａ，Ｂ．Ｇｉｒｍａ，Ｅ．Ｓ．Ｋｉｓｓｅｌ，Ａ．Ｎ．Ｌｅｖｙ，２０１４：１⁃３２．

Ｓ．ＭａｃＣｒａｃｋｅｎ，Ｐ．Ｒ．Ｍａｓｔｒａｎｄｒｅａ，ａｎｄＬ．Ｌ．Ｗｈｉｔｅ（ｅｄｓ．）］．ＣａｍｂｒｉｄｇｅＵｎｉｖｅｒｓｉｔｙＰｒｅｓｓ，Ｃａｍｂｒｉｄｇｅ，ＵｎｉｔｅｄＫｉｎｇｄｏｍａｎｄＮｅｗＹｏｒｋ，ＮＹ，ＵＳＡ，

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１１期　　　何彦龙　等：基于ＧＩＳ的长江口海域生态系统脆弱性综合评价　３９２５

［５］　ＤｅＬａｎｇｅＨＪ，ＳａｌａＳ，ＶｉｇｈｉＭ，ＦａｂｅｒＪＨ．Ｅｃｏｌｏｇｉｃａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｒｉｓｋａｓｓｅｓｓｍｅｎｔ—Ａｒｅｖｉｅｗａｎｄｐｅｒｓｐｅｃｔｉｖｅｓ．ＳｃｉｅｎｃｅｏｆｔｈｅＴｏｔａｌ［６］　ＬｉｕＤ，ＣａｏＣＸ，ＤｕｂｏｖｙｋＯ，ＴｉａｎＲ，ＣｈｅｎＷ，ＺｈｕａｎｇＱＦ，ＺｈａｏＹＪ，ＭｅｎｚＧ．Ｕｓｉｎｇｆｕｚｚｙａｎａｌｙｔｉｃｈｉｅｒａｒｃｈｙｐｒｏｃｅｓｓｆｏｒｓｐａｔｉｏ⁃ｔｅｍｐｏｒａｌ

６１２⁃６２５．

Ｅｎｖｉｒｏｎｍｅｎｔ，２０１０，４０８（１８）：３８７１⁃３８７９．

ａｎａｌｙｓｉｓｏｆｅｃｏ⁃ｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｃｈａｎｇｅｄｕｒｉｎｇ１９９０⁃２０１０ｉｎＳａｎｊｉａｎｇｙｕａｎｒｅｇｉｏｎ，Ｃｈｉｎａ．ＥｃｏｌｏｇｉｃａｌＩｎｄｉｃａｔｏｒｓ，２０１７，７３（２）：

［７］　ＥｎｅａＭ，ＳａｌｅｍｉＧ．Ｆｕｚｚｙａｐｐｒｏａｃｈｔｏｔｈｅｅｎｖｉｒｏｎｍｅｎｔａｌｉｍｐａｃｔｅｖａｌｕａｔｉｏｎ．ＥｃｏｌｏｇｉｃａｌＭｏｄｅｌｌｉｎｇ，２００１，１３６（２／３）：１３１⁃１４７．［８］　郝永红，周海潮．区域生态环境质量的灰色评价模型及其应用．环境工程，２００２，２０（４）：６６⁃６８．

ＳｃｉｅｎｃｅｏｆｔｈｅＴｏｔａｌＥｎｖｉｒｏｎｍｅｎｔ，２００４，３２７（１／３）：１０５⁃１２２．

［９］　ＰａｒｋＹＳ，ＣｈｏｎＴＳ，ＫｗａｋＩＳ，ＬｅｋＳ．Ｈｉｅｒａｒｃｈｉｃａｌｃｏｍｍｕｎｉｔｙｃｌａｓｓｉｆｉｃａｔｉｏｎａｎｄａｓｓｅｓｓｍｅｎｔｏｆａｑｕａｔｉｃｅｃｏｓｙｓｔｅｍｓｕｓｉｎｇａｒｔｉｆｉｃｉａｌｎｅｕｒａｌｎｅｔｗｏｒｋｓ．［１０］　ＷａｎｇＸＤ，ＺｈｏｎｇＸＨ，ＬｉｕＳＺ，ＬｉｕＧＷ，ＷａｎｇＺＹ，ＬｉＭＨ．ＲｅｇｉｏｎａｌａｓｓｅｓｓｍｅｎｔｏｆｅｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｉｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ：［１１］　廖富强，刘影，叶慕亚，郑林．鄱阳湖典型湿地生态环境脆弱性评价及压力分析．长江流域资源与环境，２００８，１７（１）：１３３⁃１３７．

ＥｃｏｌｏｇｉｃａｌＩｎｄｉｃａｔｏｒｓ，２０１７，７８（７）：４０５⁃４１５．

ｄｅｖｅｌｏｐｍｅｎｔａｎｄａｐｐｌｉｃａｔｉｏｎｏｆａｎｅｗｍｅｔｈｏｄ．ＪｏｕｒｎａｌｏｆＡｒｉｄＥｎｖｉｒｏｎｍｅｎｔｓ，２００８，７２（１０）：１９２９⁃１９３９．

［１２］　ＺｏｕＴＨ，ＹｏｓｈｉｎｏＫ．ＥｎｖｉｒｏｎｍｅｎｔａｌｖｕｌｎｅｒａｂｉｌｉｔｙｅｖａｌｕａｔｉｏｎｕｓｉｎｇａｓｐａｔｉａｌｐｒｉｎｃｉｐａｌｃｏｍｐｏｎｅｎｔｓａｐｐｒｏａｃｈｉｎｔｈｅＤａｘｉｎｇ′ａｎｌｉｎｇｒｅｇｉｏｎ，Ｃｈｉｎａ．［１３］　陈吉余，陈沈良．近２０年来长江口生态环境变化．净水技术，２００２，２１（Ｓ１）：１⁃３．

［１４］　单秀娟，线薇薇，武云飞．长江河口生态系统鱼类浮游生物生态学研究进展．海洋湖沼学报，２００４，４８（４）：８７⁃９３．

［１５］　操文颖．长江口“健康状况”诊断研究／／第三届湖北科技论坛水利分会场学术研讨会论文集．武汉：湖北省水利学会，２００５：１７４⁃１７９．

Ｐ０２０１８０５３１５３４６４５０３２３７２．ｐｄｆ

［１６］　中华人民共和国生态环境部．２０１７年中国生态环境状况公报．２０１８：３３⁃３４．ｈｔｔｐ：／／ｗｗｗ．ｍｅｅ．ｇｏｖ．ｃｎ／ｈｊｚｌ／ｚｇｈｊｚｋｇｂ／ｌｎｚｇｈｊｚｋｇｂ／２０１８０５／［１７］　刘正佳，于兴修，李蕾，黄玫．基于ＳＲＰ概念模型的沂蒙山区生态环境脆弱性评价．应用生态学报，２０１１，２２（８）：２０８４⁃２０９０．［１９］　黄方，刘湘南，张养贞．ＧＩＳ支持下的吉林省西部生态环境脆弱态势评价研究．地理科学，２００３，２３（１）：９５⁃１００．［２０］　沈珍瑶，杨志峰，曹瑜．环境脆弱性研究述评．地质科技情报，２００３，２２（３）：９１⁃９４．

［１８］　徐庆勇，黄玫，陆佩玲，李仁强．基于ＲＳ与ＧＩＳ的长江三角洲生态环境脆弱性综合评价．环境科学研究，２０１１，２４（１）：５８⁃６５．

［２１］　刘康，欧阳志云，王效科，徐卫华，苗鸿．甘肃省生态环境敏感性评价及其空间分布．生态学报，２００３，２３（１２）：２７１１⁃２７１８．［２２］　欧阳志云，王效科，苗鸿．中国生态环境敏感性及其区域差异规律研究．生态学报，２０００，２０（１）：９⁃１２．

［２３］　王效科，欧阳志云，肖寒，苗鸿，傅伯杰．中国水土流失敏感性分布规律及其区划研究．生态学报，２００１，２１（１）：１４⁃１９．［２４］　李东梅，吴晓青，于德永，高正文，吴钢．云南省生态环境敏感性评价．生态学报，２００８，２８（１１）：５２７０⁃５２７８．［２５］　潘竟虎，董晓峰．基于ＧＩＳ的黑河流域生态环境敏感性评价与分区．自然资源学报，２００６，２１（２）：２６７⁃２７３．

［２６］　陈杰，欧阳志云，郑华，徐卫华．淮河流域脆弱生态区生态系统特征及区划．中国人口·资源与环境，２０１０，２０（１０）：１６９⁃１７４．［２７］　肖荣波，欧阳志云，王效科，赵同谦．中国西南地区石漠化敏感性评价及其空间分析．生态学杂志，２００５，２４（５）：５５１⁃５５４．［２８］　黄宝荣，欧阳志云，张慧智，张利华，郑华．海南岛生态环境脆弱性评价．应用生态学报，２００９，２０（３）：６３９⁃６４６．［２９］　王保栋．长江口及邻近海域富营养化状况及其生态效应［Ｄ］．青岛：中国海洋大学，２００６：６２⁃７７．［３０］　张勇．基于ＧＩＳ的长江口及邻近海域环境时空分析［Ｄ］．青岛：中国海洋大学，２００８：４１⁃６３．

［３１］　崔利芳，王宁，葛振鸣，张利权．海平面上升影响下长江口滨海湿地脆弱性评价．应用生态学报，２０１４，２５（２）：５５３⁃５６１．

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ