城生态系统健康评价

1、6结论10参考文献10致谢11城市生态系统健康评价引言近年来城市化进程日益加快，在城市经济持续、快速发展的同时，不可避免地会产生一系列的生态环境问题。如何有效地协调人口、资源、环境、经济（PRED）,在资源约束条件下寻求最优的生态、经济、社会三者的综合效益，成为人们普遍关注的问题。本文通过建立生态系统健康评价的压力-状态-响应概念框架的指标体系及评价模型，旨在确定各项指标值及其权重，从而定量评估各项指标对城市生态系统健康的影响程度，并进一步绘制压力、状态、响应指标分布图，形象的表明烟台市生态系统健康程度等级。为很好的协调城市化进程与生态环境的关系问题提供了依据，同时对于指导城市规划、监控城市生

2、态环境的可持续发展都具有一定的现实意义和实际价值。1城市生态系统健康的内涵当前对于生态系统健康的内涵，比较认可的是由Rapport等总结为：“以符合适宜的目标为标准来定义一个生态系统的状态、条件或表现1”。就目前对于生态系统健康的理解而言，其应该包含两方面的含义：在满足人类社会的合理要求的同时又能满足生态系统本身的自我恢复与维持。对于城市生态系统健康的理解，具不仅意味着为人类提供生态服务的由自然环境和人工环境组成的生态系统的健康和完整，也意味着城市居住着（包括人群和其他生物）的健康和社会健康2'3。前者是后者的基础，后者是前者的目标。在此基础上对于生态系统健康概念的理解应包括以下几个方

3、面：其一，城市生态系统结构合理，城市生态系统内生产生活和周围环境之间的物质和能量交换形成良性循环。其二，系统内功能高效，物质、能量、信息高效利用。其三，人与自然和谐，人类社会和自然环境高度和谐，自然、技术、人文充分融合。其四，废水、废气、废弃物被严格控制在相应的环境承载力范围内，城市生物的健康和成长不受不良影响4。2研究区域概况烟台市地处山东半岛中部，位于东经119°34'121°57',北纬36016'38023'。烟台依山傍海，气候宜人，冬无严寒，更无酷暑，东连威海，西接潍坊，西南与青岛毗邻，北濒渤海、黄海，与辽东半岛对峙，并与大连隔海相

4、望，共同形成拱卫首都北京的海上门户。全市土地面积13745.95平方公里，其中市区面积2643.60平方公里。烟台市辖芝景、福山、牟平、莱山区、经济技术开发区及长岛县，代管龙口、莱阳、莱州、蓬莱、招远、栖霞、海阳7个县级市。3研究方法3.1 数据来源及处理3.1.1 影像数据遥感影像数据主要是选择美国landsatTM影像，分辨率28.5mX28.5m,时相为1986年、1995年和2000年三个时段的影像，成像时间主要是选取植被信息比较好的7、8、9三个月。3.1.2 其他数据主要包括烟台地区高程、水系、居民点、水土流失、道路等1:5万矢量数据。烟台统计年鉴以及经济、社会方面的其它文字资料和

5、图件资料。3.2 评价单元的确定在确定评价单元方面本文采用面状矢量单元作为评价单元，即以行政单元作为计算单元。针对烟台具体情况，文章主要是对芝景区、福山区、莱山区、牟平区、开发区五区作为评价单元，并分别计算各个单元的评价指标值及其权重。3.3 评价指标体系的建立以及评价指标处理本文以“压力-状态-响应”(PSR为评价模型，在综合考虑自然、经济、社会的情况下，从压力、状态、响应这三个层次共选取了人口密度、耕地与建设用地面积比、NDVI综合弹性值、多样性指数、均匀度指数、土壤侵蚀、人均GDR8个单项指标5作为生态系统健康评价指标，各指标如下图1所示。以下简单介绍其中几个指标的计算公式。生态系统健康

6、评价压力指标人类干扰指数人口压力活力人口密|耕地与度建设比NDVI状态指标响应指标恢复力综合弹性值组织多样性指数自然系统土壤侵人类系统GDP均匀度指数3.3.1 NDVI的计算公式图1生态系统健康评价指标借用美国Landsat5TM影像3波段和4波段数据，在ERDAS中可计算出NDVI值。NDVI计算公式如下：NDVI=TM4TM3,对计算出的所有NDVI值，取其均值作为评价指TM4TM33.3.2 综合弹性值的计算公式综合弹性值是通过生态系统的恢复力来体现的，依据不同土地利用类型对生态恢复的贡献度，分别赋以不同的恢复力值6（表1所示）。表1不同土地利用类型的恢复值，地利用类恢复力特征有林地、

7、水面O.81.O指对维持生态系统恢复力有决定意义的地物类型。林地和水域是在维持区域的稳定性和保持区域的调节滩地、沼泽0.6O.8能方面有极其重要作用的生态类型。水田旱地建设用地0.4O6O.2O48O.2这种生态类型对维持生态系统恢复力有重要作用，对人类社会系统提供重要的物质和活动场所，如果利用不好，则容易导致生态恢复力下降。计算综合弹性值根据公式：ECQs=£SiXP,其中ECQs表示综合弹性值，Si表示土地利用类型i所占总面积的比值，R则表示不同土地利用类型恢复力的分值，m羡示土地利用类型的数目。3.3.3 多样性指数的计算公式多样性指数即指土地利用类型的多样化程度，其计算公式如

8、下所示：H=-ZPJnP,其中H表示多样性指数，R表示土地类型i所占总面积的比值，m表示土地利用类型的数3.3.4 均匀度指数的计算公式均匀度反映区域内各土地利用类型在面积上分布的不均匀程度。表达式:E=H/Hmax其中E表示均匀度指数，H表示多样性指数，Hma刎表示最大多样性指数（Hmax=lnm）。3.3.5 土壤侵蚀度指数的计算公式表2土壤侵蚀程度分级表水蚀状态分值轻度水蚀区0.4中度水蚀区0.8强烈度水蚀区1m依据公式M至ABi,其中Ai表示不同强度水蚀区所占区域面积的比例，B1表示对应i1侵蚀区的分值。3.4 评价指标值的标准化处理依据上述公式在计算所有评彳/T指标值时借鉴以下（表3

9、）数据（备注：面积单位为平方公里，GDPI位为万元）。表3数据表指标芝紧区福山区莱山区牟平区开发区人口总数687234257339205947472126174863土地总面积172.02482.83254.741588.45224.26GDP18396231317146123678817052645806401耕地与建设用地面积比例0.75722.30952.92832.28741.6566林地、水面所占面积比例0.08670.09000.02800.11310.0280水田所占面积比例0.00140.00030.00580.00470.0031滩地、沼泽所占面积比例0.00620.0142

10、0.04450.00990.0028旱地所占面积比例0.31610.44080.83850.54800.5182在借鉴以上数据的基础上，采用上述相应指标计算公式，可得出以下指标值数据表4（备注：人口密度单位为人/平方公里，人均GDPI位为元）。表4指标值数据表指标值芝紧区福山区莱山区牟平区开发区人口密度3995.082532.9806808.4596297.2243779.7333耕地与建设用地面积比0.75722.30952.92832.28741.6566NDVI0.4030.38700.39980.40120.4676综合弹性值0.21980.17420.31610.36860.2119

11、多样性指数0.98120.95690.57190.63300.8043均匀度指数1.64031.68202.81422.54262.0010土壤侵蚀0.30800.24720.32100.28260.2674GDP6105952076589663611236941由于各个参评指标数据性质、量纲不同，需要对资料进行无量纲化处理，本文通过极值归一化处理形成无量纲数据。对于所获取的各项指标值进行归一化处理，采用如下计算公式：P?=XiXmin7其中P?表示归一化处理后的各指标值。对各指标归一化处理Xmax-Xmin后其具体值如表5所示。表5归一化处理后的指标值综合指标单因子指标芝旱区福山区莱山区牟平

12、区开发区人口密度10.06380.138300.1305压旨耕地与标建设用地面积00.102310.56230.0593比NDVI100.83330.66670.7687综合弹0.234600.729910.1939状态指性值多样性标指数10.940600.14930.5678均匀度指数00.035510.76860.307响应指土壤侵0.36110.91670.680601标蚀人均GDP0.08290.05310.0760013.5 评价指标权重的确定计算各项评价指标的权重，即各因素因子的权重时采用德尔菲法（Delphi法）。由专家根据经验对各项因素的重要性进行打分；根据专家打分结果，计算各

13、因素的平均重要性和各专家对其重要性的偏离度，如果偏离度超出可以接受的范围，则将结果反馈给专家重新打分，直至偏离度降到某一限值以下；根据最终打分结果确定各因素权重。本文在参考专家打分结果的基础上恰当的选取各项评价指标的权重。其具体指标权重值如表6所小。表6指标值权重综合指标权重单因子指标权重人口密度0.4455压力指标0.2980耕地与建设用地面积比0.5545NDVI0.2885综合弹性值0.1873状态指标05445.多样性指数0.2255均匀度指数0.2987土壤侵蚀0.6428响应指标0.1575人均GDP0.35724生态系统健康综合评价4.1生态系统健康压力评价本文所采用压力指标包括

14、人口密度、耕地与建设用地面积比，在计算压力指标值时2采用如下公式：Si=zPiCi,其中Pi表示归一化处理后的指标值，Ci表示对应各指标i1的权重。4.2生态系统健康状态评价所采用状态指标包括NDVk综合弹性值、多样性指数、均匀度指数，在计算压力指4标值时，也采用如下公式：S2=£PiG,其中Pi表示归一化处理后的指标值，G表示对i1应各指标的权重。4.3 生态系统健康响应评价响应指标包括土壤侵蚀、人均GDP同理采用如上公示。在分别计算压力、状态、响应指标值的基础上可获得表7。表7各区指标值指标值芝旱区福山区莱山区牟平区开发区压出旨标0.44550.08510.04410.31180

15、.0910状态指标0.55790.22270.05100.64290.6606响应指标0.26170.60820.38550.36870.45424.4 生态系统健康评价等级确定根据压力指标、状态指标、响应指标及其各自的权重进行求和，可分别求得芝景区、福山区、莱山区、牟平区、开发区的生态系统健康分值，如下表8。表8总体健康分值莱山区0.47780.24240.10160.50390.4584在确定各区生态系统健康等级时有以下依据：按照分值高低0.8、0.6、0.4、0.2、0.1依次分为五级，即一级很健康、二级健康、三级亚健康、四级不健康、五级病态。生态系统健康分级及健康状态的含义如表9所示。

16、表9生态系统健康状态分级健康分健康状生态系统特征生态结构非常合理、系统活力十分强，受外界压力较小，无一级很健康生态异常出现，生态系统的生态功能很完善，系统极稳定，处于可持续状态生态系统结构合理、格局比较完整，受外界压力较小，基本二级健康无生态异常出现，生态系统生态功能较完善，系统尚稳定，生态系统可持续生态结构比较合理，系统尚稳定，受外界压力比较大，基本三级亚健康上已经接近生态阈值，且敏感地带较多，已有少量的生态异常出现，可发挥基本的生态功能，生态系统基本可维持生态结构开始出现缺陷，系统活力比较低，此时受外界压力四级不健康大，生态异常较多，生态功能已不能满足维持湿地生态系统的需要，生态系统己开始退化生态结构极不合理，自然植被斑块破碎化严重，系统活力极五级病态低，出现大面积的生态异常区，此时生态系统已经受到严重恶化5烟台市生态系统健康状况通过以上的计算分析，可分别绘制得压力、状态、响应指标分布图，并最后求得各区生态系统健康状况的总体健康等级，并进一步绘制烟台市生态系统健康综合分布图。通过分析也可知，属于二级健康的有3个，即芝景区、牟平区和开发区；属于三级健康的有1个，即福山区；属于四级健康有1个，即莱山区。以下图2、3、4、5详细表明各指标分布图及总体健康状态分布图。图2烟台市压力指标分布图图3烟台市状态指标分布图图4烟台市响应指标分布图图5烟台市生态系统健康综