生态风险是指生态系统面临未来损失的可能性及其严重程度，即在一定范围内，具有不确定性的事故或灾害对生态系统及其组分可能产生的不利作用，包括生态系统结构和功能的损害，从而危及生态系统的安全和健康^[1-3]。生态风险评价起于20世纪80年代初，在90年代基本形成，是目前环境科学领域的研究热点，其主要集中在重金属或洪涝干旱等自然灾害作为风险源的生态风险评价中^[4-7]；而水土流失是我国最大的生态环境问题，对我国的生态安全有很大的影响，国内对于水土流失的研究，最常见的是通过计算土壤流失量来评价一个地区水土流失的现状，最常用的方法是利用美国的通用土壤流失方程。经过大量研究证实，对于中国特有的地质地貌，这种方法所计算出的流失量，并不能很准确的描述区域水土流失状况，而且此方程中，对于水土保持工程因子的计算存在很大困难，这一因子恰好决定着当前开展水土保持工程的程度和规模，以便以最合理的人力物力投入，得到最好的水土流失预防效果。另有学者通过测定其他一些指标，以求间接代替以上因素。例如一些学者把污染物作为主要的风险源^[8-11]，而实际上对于一个完整的生态系统来说，由于环境的差异，相同污染物对不同地区的风险是不可能相同的，这样评价会造成评价结果的片面性；除此之外，一些极端的自然灾害也在很大程度上影响着区域生态风险预测的可靠性；对于山区的生态风险评价，很多指标不能很好的被观测记录，这就导致对于这些不发达地区的生态风险评价无法完成，生态风险无法评估。

基于早期的生态风险评估理论，笔者克服了传统的单因子评价、综合指数评价和生态风险指数评价的不足，综合考虑在生态风险预测中遇到的种种问题，寻找符合北京地区，特别是北京山区水土流失特点的指标，构建一种新的生态风险评估方法和评估体系，为生态风险评估提供参考依据。

北京地区E 39°57′，N 116°26′主要分为山区和平原区。生态风险可能发生的地区多为北京山区，其主要分布于北京西部和北部的房山、延庆、密云、怀柔、昌平、平谷和门头沟等7个远郊区县，总面积约为1.03万km²，占北京市面积的61.29%。北京山区分为2大山系，分布为西南端的太行山系与东北段的燕山山系。北京山区属于典型的北温带半湿润大陆性季风气候，夏季高温多雨，冬季寒冷干燥，春秋短促。对于山区而言，由于其跨度不超过200 km，山幅宽约为50 km，并没有跨越不同的经纬度，同时都处于内陆，气候在水平方向上的差异并不明显，影响气候的主要因素为地形。北京山区的年平均气温为11 ℃，年降雨量在450~660 mm之间，为华北地区降雨最多的地区之一，但是降雨年度分布并不平均，全年有80%的降雨集中在夏季6—8月。特殊的降雨分布特征和地质地貌条件，使得北京山区的水土流失问题极为特殊，导致山区生态环境恶化，生产力下降^[14]。

风险是一种可能造成的损失值，必须处于一定区域和给定时间段内，由某种自然灾害而造成，其结果是人民生命财产和经济活动遭受损失，其风险表达式^[12-13]为

式中：R为风险事件的本身；H为风险发生的可能性；V为风险造成的损失。

当存在一个既定的风险事件k时, 其风险的定义表达为

而对于相互独立的一组风险事件, 风险的表达式为

在评价一个区域的生态风险时，往往由于指标无法获取，而得不到较为精确的评价结果，为了更为有效的对研究区域进行生态风险评价，需要对现有的评价模型进行修改，突出水土流失等灾害对生态风险评价结果的影响；因此，对式 (1) 进行修改

式中：R为生态风险的量化值；I为地区生态重要程度指数量化值；H为水土流失灾害以及人类活动发生的概率；T为区域生态易损性指数量化值。

通过对前人研究成果的归纳与总结，首先对北京地区生态造成影响的指标进行选择，形成初步筛选结果；进而采用专家打分法，同时，考虑指标的可获取性，对初选指标进行再次筛选与分类，得到指标体系如表 1所示。

表 1(Table 1)

表 1 指标体系筛选结果 Table 1 Results of screened indexes 指标类别 Classification 指标名称 Index 生态危险性指标 (H) 生态重要性指标 (I) 生态易损性指标 (T) 气象因素 Meteorological elements 降水、风速风向 Precipitation, wind speed, and direction 多年平均降雨量 Average annual rainfall 地质和土壤因素 Geological and soilelements 岩性、土壤类型 Lithology and soil type 地层岩性、土壤类型 Lithology and soil type 地形因素 Topographical elements 坡度、相对高差、坡长、坡向、流域面积、沟床纵比降、流域拓扑关系 Slope gradient, relative elevation, slope length, slope direction, drainage area, gully bed longitudinal slope, and basin topology 坡度、坡向Slope gradient and slope direction 生境敏感性指数 Habitat sensitivity index 植被因素 Vegetation elements 植被类型、植被覆盖度 Vegetation type and vegetation coverage 植被类型、植被覆盖度Vegetation type and vegetation coverage 土地利用因素 Land-use elements 土地利用类型 Land-use type 风景区分布Distribution of scenic area 人为因素 Human elements 人口密度分布、主要人类活动分布Distribution of population density and human activity 人类活动 Human activity 受威胁人数、潜在经济损失Threatened population and potential economic loss 注：H为生态危险性指标，I为生态重要性指标，T为生态易损性指标。以下类同。Note: H refers to ecological hazard index, I to ecological importance, and T to ecological vulnerability. The same below.

表 1 指标体系筛选结果 Table 1 Results of screened indexes

通过查阅大量的历史资料，建立指标分级体系 (表 2)。利用Arc GIS 9.3中的重分类功能，对以上12个指标进行分级，进而得到不同指标的栅格数据图。

表 2(Table 2)

表 2 各指标分级标准以及分类 Table 2 Grading and classification of each index 评价指标 Evaluation index 等级划分Grading and classification of each indicator 极低Very low 低Low 中等Middle 高High 极高Very high 生态危险性 指标 (H) 多年平均降雨量Average annual rainfall 均分Average 地层岩性Lithology 花岗、闪长岩Granite and diorite 辉长、玄武岩Gabbro and basalts. 石灰、页岩Limestone and sand shale, 白云岩Dolomite 砂岩Malmstone 土壤类型Soil type K=0, 0.001 K=0.010 K=0.074 K=0.147, 0.177 K=0.211, 0.219, 0.222, 0.232, 0.235 坡度Slope gradient/(°) ＜5 5~8 8~15 15~25 >25 坡向Slope direction 阴坡Cloudy 半阴Semi-cloudy 平地Level ground 半阳Semi-sunny 阳坡Sunny 植被类型Vegetation type 阔叶林Broad-leaved forest 针叶林Coniferous forest 灌木林Shrub forest 草地Glass land 农田Farmland 植被覆盖度Vegetation coverage/% 75 75~55 55~35 35~10 ＜10 人类活动Human activity 天然草地、林地，未开发利用地Natural grassland, forest land, and undeveloped land 耕地、园地、水域及水利设施用地Cultivated land, garden land, water area, and land for water conservancy facilities 住宅用地、商服用地、公共管理与公共服务用地等Residential land, commercial land, and land for public management and public service 工矿仓储用地Land for industrial storage 交通运输用地Land for transportation 生态重要性 指标 (I) 生境敏感性指数Habitat sensitivity index 其他Others 农田Farmland 草地Glass land 水体Water area 森林Forest 风景区等级Scenic 非景区Non-scenic 区级、县级景区 District and county 市级景区City 国家级景区 National 生态易损性 指标 (T) 受威胁人数Threatened population 均分Average 潜在经济损失Potential economic loss 其他Others 水体Water area 农业区 (耕地、林地、园地、草地等) Agricultural area (farmland, forestland, garden, grassland, etc.) 房屋住宅、商服用地、公共管理与公共服务用地Housing, commercial land, and land for public management and public service 工矿区、重要基础设施 (交通干线、水利工程、电力工程等) Industrial and mining areas, important infrastructure (traffic trunk, hydraulic engineering, power engineering, etc.) 等级赋值Value 1 2 3 4 5

表 2 各指标分级标准以及分类 Table 2 Grading and classification of each index

1) 多年平均降雨量。通过搜集北京地区现有的气象站点降雨资料，利用Arc GIS中的空间插值分析模块，得到多年平均降雨量等值线，并按照降雨量的大小对其进行分级。

2) 地层岩性。按照不同地层岩性抗侵蚀能力的强弱划分等级，岩层的抗侵蚀性越高，其存在的生态风险越小。

3) 土壤类型。不同土壤类型存在不同的抗侵蚀能力，采用K因子数值，对土壤类型进行分级。

4) 坡度和坡向。坡度主要用于水力侵蚀类型的面蚀分级，依据“水土保持行业标准”，坡度分为＜5°、5°~8°、8°~15°、15°~25°、25°和＞25°等5个等级。坡向划分一般只说阳坡和阴坡，阴阳坡划分为阴坡 (0°~45°)、半阴坡 (45°~135°)、阳坡 (135°~225°)、半阳坡 (225°~315°)、阴坡 (315°~360°)。

5) 植被类型和植被覆盖度。植被类型和植被覆盖度的变化，能够直接影响区域水土流失分布规律，根据2005年北京土地利用现状图，对植被类型进行分类和提取。利用ENVI 4.8中的植被覆盖度计算模块，根据2010年北京地区ETM+影像，计算植被覆盖度。

6) 人类活动。人类活动的分级是根据土地利用类型划分的，不同的土地利用类型，表征不同的人类活动干扰程度。

7) 生境敏感性指数。生物多样性维持中包括体系中的授粉、生物控制、栖息地和基因资源4项功能。根据中国陆地生态系统单位面积生态服务价值当量表 (每km²)，将生物多样性保护作为评价区域生境敏感性指数的指标特征，对于不同类别的土地利用类型赋以不同值，来代替区域生境敏感性指数。

8) 风景区分布。按照风景区的景区等级，对风景区分布进行等级划分。

9) 受威胁人数。人是风险的受体之一，受威胁人数越多，易损性越大。各区县的统计数据采用插值得到。根据《北京市2010年第六次全国人口普查主要数据公报》^[15]中的各区县人口普查数据，计算得到各区县的人口密度数据。

10) 潜在经济损失。根据土地利用类型进行分级，不同的土地利用情况，遭受水土流失时的损失不同，易损性也存在差异。

利用层次分析法，对12个指标进行权重的分配，得到的分配结果见表 3。利用Arc GIS的重分类模块，对3个指标组和最后的风险等级进行划分，划分的依据为自然分类，等级由低到高为安全 (Ⅰ)、较安全 (Ⅱ)、轻度危险 (Ⅲ)、中度危险 (Ⅳ) 和高度危险 (Ⅴ)。

表 3(Table 3)

表 3 风险源指标权重划分 Table 3 Weight definition of risk source index 生态危险性指标 (H) 生态重要性指标 (I) 生态易损性指标 (T) 指标名称 Index 权重 Weight 指标名称 Index 权重 Weight 指标名称 Index 权重 Weight 指标名称 Index 权重 Weight 多年平均降雨量Average annual rainfall 0.261 6 植被覆盖度Vegetation coverage 0.175 9 生境敏感性指数Habitat sensitivity index 0.750 0 受威胁人数Threatened population 0.666 7 土壤类型Soil type 0.033 1 人类活动Human activity 0.110 2 风景区类型Scenic type 0.250 0 潜在经济损失Potential economic loss 0.333 3 坡度Slope gradient 0.269 0 植被类型Vegetation type 0.053 8 坡向Slope direction 0.071 1 地层岩性Lithology 0.025 2

表 3 风险源指标权重划分 Table 3 Weight definition of risk source index

数据来源于北京地区1989—2009年各气象站年均降雨数据、北京市DEM和ETM影像、北京市各区县人口密度统计表、北京市09年土地利用分类图、土壤属性图、地层岩性图和北京市景区分布图。

处理软件为ArcGIS9.3和ENVI4.8。所有数据都采用WGS84坐标系，并且经过投影矫正和几何校正。植被覆盖度F=(NDVI-NDVI_S)/(NDVI_V-NDVI_S)。式中：NDVI为归一化植被指数；NDVI_S为完全是裸土或无植被覆盖区域的归一化植被指数；NDVI_V为完全被植被所覆盖区域的归一化植被指数。生境敏感性指数按照不同土地利用类型的生物多样性服务价值，采用下式计算：

式中：HIS为生境敏感性指数；l为土地利用；n为生物多样性服务当量；m为保护级别修正值。其中，自然保护区赋值为1.75，公园为1.5，风景旅游区为1.25，其他地方为1。

通过对7个危险性指标的计算 (图 1)，颜色深的区域表示生态危险性指数越高，反则生态危险性指数越低。北京生态危险性指数的分布区域主要呈条带状分布，与年均降雨量分布和山系的分布较为吻合。生态危险性指数最低的区域出现在北京北部区域，该区域为北京市密云水库所在的区域，因为主要生态环境为水域保护良好，所以该区域生态危险性指数很低。

图 1(Figure 1)

图 1 生态危险性指标等级划分图 Figure 1 Classification map of ecological hazard index

由图 2可知，生态重要性高的区域主要分布在西部和北部，而重要性低的地区主要分布在北京东南部和西北部的部分地区，而北京主城区处于中间地带。总体上看，体现西北部重要性要强于东南部，即山区要强于平原地区，造成这种现象的主要原因是植被分布。平原地区，由于人类的过渡开发，造成植被覆盖减少，生态功能下降，并且由于东南部人类活动频繁，造成对原有环境的极大影响；而西北山区多植被覆盖，由于部分地区人烟稀少，反而对原有的环境影响较小；同时，多景区的分布，人为的看管与维护，在一定程度上保护了当地的生态环境，尤其是国家级景区的分布区域，其受生态重要性最高，这充分说明该区域对人类乃至自然的作用重大，对北京地区进行生态风险评估时，需要重点突出研究。

图 2(Figure 2)

图 2 生态重要性指标等级划分图 Figure 2 Classification map of ecological importance index

由图 3可知，处于高危害区域为南部地区，并且有向北部慢慢降低的趋势。从经度方向看，西部的危害等级要高于东部。总的来说，这与区域的人口分布和经济发展情况密不可分。对于一个地区的生态风险评估，需要知道区域生态易损性的危险等级，虽然北京主城区有较好的经济发展，但是其潜在的危害也是最高的；而北京山区虽然经济不发达，但是其潜在的危害要比主城区小。因此，在发生灾害的情况下，北京山区的生态易损性要小于主城区，而在大范围内体现为：北部的生态易损性要小于南部，中部的生态易损性要高于边界地区。

图 3(Figure 3)

图 3 生态易损性指标等级划分图 Figure 3 Classification map of ecological vulnerability index

由图 4可知，整个北京地区基本处于轻度危险级别，西北和东南部边缘处于较安全级别，而东北和西南部边缘以及北部地区处于中度危险级别，高度危险级别只分布在东北部边缘地区，安全区域仅有零星分布。总的来说，整个北京地区生态风险不容乐观，生态环境已经受到一定的破坏，生态系统服务功能在某些方面不太健全，生态恢复能力差，生态灾害少有发生。整体上来看，水系附近的生态风险比较大，其次是山区和东南部平原地区，而主城区虽然有较高的受体重要性等级和暴露与危害等级，但是其生态风险等级并不高，说明其发生风险和造成损失的可能较小。而北京山区虽然经济不发达，但是由于其对整个北京地区的生态服务起到不可忽视的作用，无论从自身还是对其他区域的影响来看，都是非常重要的。虽然这部分地区易发生地质灾害，但是总的来说生态风险等级不是很高。东北部多密集水库，其水源涵养功能非常重要，而恰好这部分处于非常高的生态风险等级，所以应重点关注此地区。总之，从长远发展的角度来看，这些地区都应该加强防范措施。

图 4(Figure 4)

图 4 北京地区生态风险等级划分图 Figure 4 Classification map of ecological risk in Beijing

基于北京地区水土流失特点，分析构建符合北京地区生态评价指标体系，评价北京地区生态风险等级。所构建的生态风险评价体系，克服了传统的单因子评价、综合指数评价和生态风险指数评价的不足，为生态风险评估提供新的研究方法。研究从影响生态风险的因素入手，判断哪些区域容易发生灾害以及发生的概率等。利用层次分析法，在GIS技术的强有力的支持下，采用空间分析、重分类和叠加分析、以及栅格分析等定性或定量分析的方法，结合实地调查结果，判断北京山区生态风险的等级，进而采取相应的预防措施。评估结果表明，北京山区大部分处于轻度风险等级，存在高度危险的地区仅在密云少有分布，北部和南部的边缘地区处于轻度危险等级，而中度危险等级分布在西北和东南部分地区。水源地是北京重点保护区域，其次是山区和东南部平原。