生态学报  2016, Vol. 36 Issue (22): 7453-7461

文章信息

吴未, 陈明, 范诗薇, 欧名豪
WU Wei, CHEN Ming, FAN Shiwei, OU Minghao.
基于空间扩张互侵过程的土地生态安全动态评价——以(中国)苏锡常地区为例
A dynamic approach to land ecological security assessment: a case study of Su-Xi-Chang area, China
生态学报[J]. 2016, 36(22): 7453-7461
Acta Ecologica Sinica[J]. 2016, 36(22): 7453-7461
http://dx.doi.org/10.5846/stxb201601040021

文章历史

收稿日期: 2016-01-04
修订日期: 2016-07-15
基于空间扩张互侵过程的土地生态安全动态评价——以(中国)苏锡常地区为例
吴未1,2, 陈明1, 范诗薇1, 欧名豪1,2     
1. 南京农业大学, 土地管理学院, 南京 210095;
2. 农村土地资源利用与整治国家地方联合工程研究中心, 南京 210095
摘要: 从景观单元水平过程出发,通过采用生态用地和建设用地两种不同源空间扩张互侵过程的方法,提出了以分区形式的土地生态安全动态评价方法。以快速城市化苏锡常地区为研究区域,分别以生态系统服务价值高的水域和林地及建设用地现状为扩张源,生态系统服务价值和生态风险值为阻力赋值依据,在模拟不同生态安全格局水平和城市发展模式的基础上测算互侵结果并进行区划。结果表明:1)不同城市发展模式下建设用地空间扩张随生态用地从底线到满意到理想水平变化呈空间集聚收敛趋势;生态用地呈似圈层状空间集聚趋势。2)城市发展模式对土地生态安全水平具有显著影响;高等级源优先发展模式下土地资源利用更集约。3)将研究区划分为生态核心区(面积占比35.27%)、生态安全区(29.07%)、缓冲区(7.76%)及建设区(27.90%),提出相应土地用途管制措施。探讨了新方法不足和未来研究方向。
关键词: 土地生态安全     动态评价     空间扩张互侵     生态系统服务     生态风险     快速城市化地区    
A dynamic approach to land ecological security assessment: a case study of Su-Xi-Chang area, China
WU Wei1,2, CHEN Ming1, FAN Shiwei1, OU Minghao1,2     
1. College of Land Management, Nanjing Agricultural University, Nanjing 210095, China;
2. National & Joint Local Engineering, Research Center for Rural Land Resources Use and Consolidation, Nanjing 210095, China
Abstract: Land ecological security (LES) assessment is a main focus of LES research. The aim of this study was to develop a dynamic approach to achieve an objective assessment. Su-Xi-Chang area was selected as the study area. Based on the source-sink landscape theory, ecological land and construction land were treated as the source and sink, respectively. The minimum cumulative resistance model was applied to simulate the dynamic changes of these landscapes. The resistance values of the ecological land were obtained from ecosystem service values of different land-use types, whereas those of the construction land were derived from ecological risk values. The spatial expansion and simultaneous invasion of both construction land and ecological land were simulated with ArcGIS software. In the case of ecological land expansion, three scenarios were simulated, i.e., low-, middle-, and high-level LES. In the case of construction land expansion, three urban development models were designed, i.e., metropolitan development in priority (MD), small and medium cities development in priority (SD), and all urban sources development in equilibrium (ED). In the case of invasion, nine scenarios were formed from the spatial overlapping of these two kinds of expansions and their respective scenarios. The study area was zoned into five levels of LES according to the invasion results, i.e., suitable construction (very low), ecologically fragile (low), buffer (middle), ecological security (high), and ecological core (very high) zones. The results of this study showed that:(a) Along with the scenarios of LES ranging from low-to high-level, the expansions of construction land in different development models displayed a significant convergence trend of spatial agglomeration; the modeling of ecological land, however, indicated a trend of circular sprawling from the core to the outward. In the low-level scenario, the percentage of suitable construction zoned land was the highest in the MD model, at 27.90%, and the lowest in the SD model, at 21.59%. In the high-level scenario, the percentage of suitable construction zoned land was the lowest in the MD model, at 9.42%, and the highest in the SD model, at 13.10%. (b) In the same ecological security scenario, the different urban development models displayed different impacts on construction land expansion. In the low-level scenario, the combined percentage of suitable construction zoned land and ecologically fragile zoned land was the lowest in the MD model, at 38.24%, and the highest in the SD model, at 46.03%. In the high-level scenario, the corresponding value was also the lowest in the MD model, at 14.41%, and the highest in the ED model, at 21.20%. (c) Urban development models have significant influence on the levels of LES. Among the three urban development models in this study, the MD model benefitted intensive land use and found the land to be suitable for rapid urbanization, improving the area's LES from the current level to a higher level. A LES zoning plan was recommended, as well as measures of land-use control. In the plan, the area percentages of the suitable construction, buffer, ecological security, and ecological core zones were 27.90%, 7.76%, 29.07%, and 35.27%, respectively. The developed approach was found to be effective in dynamically assessing LES and beneficial in the differential management of land resources. The dynamic approach investigated in this study highlights the methodology in the optimization of ecological land patterning. Finally, issues on how to improve this developed method and future research directions have also been discussed.
Key words: land ecological security     dynamic assessment     spatial expansion and invasion     ecosystem service     ecological risk     rapidly urbanizing areas    

生态安全是社会经济环境可持续发展不可缺少的基础[1], 也是学界研究热点问题[2-3], 研究主题包括土地生态安全问题[4]。土地生态安全评价作为土地生态安全核心内容具有重要研究意义和价值[5]

国内外土地生态安全评价研究多采用PSR概念框架模型与3S技术等空间分析手段相结合的方法[4], 通过系统分解综合考量生态系统与社会经济系统间相互关系, 构建出适宜评价指标体系并赋权, 对评价单元逐项评价后经空间叠加、模型计算等步骤实现[6-12]。但是, 主观性较强的指标选取和赋权会直接影响评价结果[11, 13];行政区作为评价基本单元[9-11]不能反映出单元内部差异[5];生态因子叠加(景观单元垂直过程)不能体现出景观水平动态变化[14-17];时间序列趋势外推预测方法因较少考虑突发性情景也会影响评价结果[18-19]。2000年我国城市化率为36.2%, 2011年末增至51.3%[20]。快速城市化下建设用地急剧扩张、农用地破碎化明显、土地生态安全脆弱化加速, 影响了区域持续发展[21-22]。有必要探究新的土地生态安全评价方法[23-24]拓展土地生态安全研究。

本文以快速城市化苏锡常地区为研究区域, 从景观单元水平过程出发, 采用不同源空间扩张互侵过程模拟和分区的方法, 进行土地生态安全动态评价, 以期拓展土地生态安全评价方法论的同时为土地用途空间管制实施提供理论依据。

1 研究区域 1.1 研究区概况

苏锡常地区(119°08′-121°15′E, 30°46′-32°04′N)位于江苏省南部太湖之滨, 是长江三角洲腹心地带, 属长江冲积平原, 区内地势平坦(图 1)。

图 1 研究区位置示意图 Fig. 1 Location of the study area

地区总面积1.75万km2, 水域面积占32.5%;地区以江苏省约17%的国土面积和人口, 实现了约40%的地方国民生产总值和地方财政收入。随着城市化快速发展, 苏锡常地区土地利用结构和空间格局发生显著变化, 2000-2010年地区农用地比重从56.69%降至44.41%, 建设用地比重从14.71%增至27.82%, 已处于生态脆弱区[25]

1.2 数据来源与处理

数据主要包括中国科学院国际科学数据服务平台2010年苏锡常地区TM遥感数据、2010年苏锡常地区行政区划图以及苏锡常3市《统计年鉴(2011)》。遥感数据分辨率为30m, 在ENVI遥感软件支持下完成了几何校正、图像配准等处理, 并经地区2010年土地利用现状图校对。依据全国土地资源分类系统和研究需要, 将土地利用类型划分为林地、草地、水域、耕地、建设用地及未利用地6类。

2 研究方法 2.1 研究思路

PSR概念框架模型具有鲜明的生态意义[9, 12], 但对景观单元水平过程重视不足。从水平过程出发, 土地生态安全与土地利用类型变化关系密切。后者指生态用地被侵占为建设用地及将建设用地恢复为生态用地的过程, 受城镇化政策与发展模式、耕地保护与土地管理模式、生态保护等因素影响显著。依据源汇景观理论[26], 将建设用地、生态用地作为源或汇可表述景观类型间水平生态过程[27]。建设用地空间扩张入侵意味着土地生态安全水平下降, 生态风险/威胁增加;生态用地空间扩张入侵意味着土地生态安全水平提升, 生态系统服务功能和价值增强/加。

建设用地和生态用地源空间扩张互侵是一种时空同步的竞争性控制过程, 通过克服阻力实现并反映出土地利用类型及生态安全变化趋势。该过程模拟就是土地生态安全动态评价。为便于测算可将扩张与互侵分开。空间扩张可采用最小累积阻力模型实现[28-29]。不同地类的生态系统服务、区域生态风险/威胁客观存在[30-31]可作为阻力赋值依据, 避免指标选取和赋权主观的问题。

2.2 生态用地源空间扩张模拟

依据相关成果[27]选取区内单位面积生态系统服务价值高和土地利用结构占比大的水域和林地为主要生态用地源。通过实现区内生态系统服务总价值最大化时的土地利用格局达到提升区域土地生态安全水平目的。从用地结构占比出发, 参照生态安全格局划分的3个安全水平(底线/低、满意/中、理想/高)及每个水平生态用地占比最小值(47.32%、70.45%和85.11%)为约束条件[32]模拟生态用地扩张。扩张时, 单位面积生态系统服务价值越高受到的阻力就越小, 反之越大[14, 33-34]。生态系统服务处于动态变化中[30]需当量校正[35]。依据相关成果[27]对研究区2010年不同地类单位面积生态系统服务价值进行极差标准化处理, 得到0-100之间的阻力值(表 1)。以上计算通过ArcGIS10.0实现。

表 1 不同土地利用类型单位面积生态系统服务价值及阻力值/(元hm-2 a-1) Table 1 Ecosystem service value and its corresponding resistance value of land-use type of Su-Xi-Chang area
类型
Land-use type
单位面积生态
系统服务价值
Ecosystem
service value
阻力值
Resistance value
水域Water 1208.30 0
林地Forest 569.69 53
草地Grassland 188.77 84
耕地Farmland 180.16 85
未利用地Unutilized land 11.51 99
建设用地Construction land 0 100
2.3 建设用地源空间扩张模拟

不同建设用地在受到相同阻力时扩张能力不同。这种扩张差异可通过建设用地等级差异反映[28]。参照相关成果[28]从区位、人均GDP、地均GDP、人口密度等指标考量将研究区建设用地源(3个地级市下辖26个行政区/县级市及开发区)由高至低划分为4个等级(表 2)。

表 2 苏锡常地区建设用地源等级划分 Table 2 Expansion source grading of construction land of Su-Xi-Chang area
等级(数量)/个
Grade and number
范围District
一级Level 1(9) 崇安区、苏州工业园区、无锡新区、虎丘区、昆山市、南长区、天宁区、北塘区、张家港市
二级Level 2(5) 戚墅堰区、江阴市、姑苏区、钟楼区、惠山区
三级Level 3(5) 常熟市、新北区、太仓市、锡山区、武进区
四级Level 4(7) 吴江区、滨湖区、相城区、吴中区、宜兴市、金坛市、溧阳市

城市发展通常包括大城市(高等级源)优先发展、大中小城市均衡发展以及中小城市(低等级源)优先发展3种模式。不同等级建设用地源受到相同阻力的扩张差异修正公式为[28]

(1)

式(1)中, Kj为源j所属等级扩张修正系数, 系数越小扩张能力越强。3种模式不同等级建设用地源修正系数比为:高等级源优先发展(0.7:0.8:0.9:1.0)、均衡发展(1.0:1.0:1.0:1.0)、低等级源优先发展(1.0:0.9:0.8:0.7)。按上述比例对不同等级建设用地源进行空间扩张模拟。

建设用地扩张与生态风险正相关:建设用地密度越高, 扩张能力越强, 生态风险越大[36]、到建设用地距离越小, 被侵占可能性越大, 生态风险越高[37-38]、生态风险越大, 阻力值越小;反之亦然。快速城市化地区生态风险测算方法很多[39-42], 考虑到建设用地和生态用地空间扩张的可比性, 选用生态损失度指数法[36]测算生态风险。主要包括:测算不同地类生态损失度指数;按2 km×2 km正方形网格将研究区等间距系统空间化采样得到4591个样方;依据生态损失度指数测算样方中不同地类生态损失度指数加权和及样方生态风险值;以样方生态风险值为建设用地扩张阻力赋值依据, 设定0-100的阻力值。

生态损失度指数计算表达式为:

(2)
(3)
(4)
(5)
(6)

式(2)中, RiGiDi分别为第i类生态系统生态损失度指数、干扰度指数和脆弱度指数。Ci、Fi、Bi分别为第i类生态系统破碎度、分离度和分维数, 干扰度指数为上述3指标归一化后的加权和;α、β、γ为对应权重。ni为景观类型i的斑块数;Ai为景观类型i的面积。Ii为景观类型i的距离指数;A为景观总面积。pijaij分别为景观类型i中第j个斑块的周长和面积。

依据相关成果[36], 对α、β、γ赋值0.6、0.3、0.1;将景观类型脆弱性分为6级归一化处理后得到对应脆弱度指数和损失度指数:未利用地(0.29, 0.401)、水域(0.24, 0.148)、草地(0.19, 0.122)、林地(0.14, 0.086)、耕地(0.10, 0.061)、建设用地(0.05, 0.032)。

样方生态风险值计算表达式为:

(7)

式(7)中, ERIk为第k个样方的景观生态风险值;Aki为第k个样方中第i类生态系统类型的面积;TAki为第k个样方的总面积。

生态风险等级划分尚没有统一标准, 故采用等间距法[36]将介于0-0.15的生态风险值由低至高分为5级:低(0, 0.03]、较低(0.03, 0.06]、中(0.06, 0.09]、较高(0.09, 0.12]和高(0.12, 0.15], 分别对应阻力值100、75、50、25和0。

以上计算通过ArcGIS10.0、FRAGSTATS4.2及Excel软件实现。

2.4 互侵结果计算及评价

建设用地和生态用地空间互侵过程反映了区域土地生态安全水平变化过程。结果由公式(8)计算得到:

(8)

式(8)中, 当MCR差值> 0且数值趋大时, 说明评价单元土地生态安全水平下降并趋于恶化;当MCR差值=0时, 说明评价单元土地生态安全水平处于临界状态;当MCR差值 < 0且数值趋小时, 说明评价单元土地生态安全水平上升并持续好转。采用自然断点法将评价单元MCR差值划分为5个等级:适宜建设区(较低水平区)、生态脆弱区(低水平区)、缓冲区(临界区)、生态安全区(较高水平区)和生态核心区(高水平区)。

3 结果与分析 3.1 不同源空间扩张结果

图 2(a, b, c)分别是建设用地高等级源优先发展、均衡发展和低等级源优先发展3种模式空间扩张结果。差异主要是对水域侵占和对非建设用地破碎化影响程度不同:高等级源多集中在环太湖北-东北方向;高等级源优先发展中, 建设用地从高等级源向四周辐射、相向集中连片、对周边水域产生明显侵占(图 2a);低等级源优先发展中, 低等级源建设用地比高等级源建设用地向四周扩张态势更明显、对水域侵占范围更广(图 2c);均衡发展中, 各等级源建设用地对周边非建设用地侵占程度较均衡(图 2b)。

图 2 不同情景建设用地及生态用地空间扩张结果 Fig. 2 Spatial expansion scenarios of construction land and ecological land

图 2(d, e, f)分别是生态用地底线、满意和理想3个安全水平空间扩张结果。总体表现为以底线水平生态用地为核心, 从满意到理想水平呈近圈层状空间扩张趋势。

3.2 互侵及评价结果

图 3(A, B, C)、(D, E, F)和(G, H, I)是建设用地高等级源优先发展、均衡发展和低等级源优先发展模式分别与生态用地底线、满意和理想水平互侵后的结果。表 3是对应土地生态安全分区情况。

图 3 苏锡常地区土地生态安全分区图 Fig. 3 Scenarios of land ecological security zoning of Su-Xi-Chang area

表 3 苏锡常地区土地生态安全分区情况 Table 3 Information of land ecological security zoning scenarios of Su-Xi-Chang area
情景Scenario 生态核心区面积及强度
Ecological core area &
its percentage
生态安全区面积
Ecological area &
its percentage
缓冲区面积
Buffer area & its
percentage
生态脆弱区面积
Ecological fragile area &
its percentage
适宜建设区面积
Suitable building area &
its percentage
/km2 % /km2 % /km2 % /km2 % /km2 %
高等级-底线Scene A 421.47 2.43 4972.04 28.68 5318.85 30.67 1793.14 10.34 4838.71 27.90
高等级-满意Scene B 5393.54 31.10 5750.45 33.16 1634.92 9.42 1302.22 7.51 3263.09 18.82
高等级-理想Scene C 10494.17 60.51 2580.90 14.88 1769.59 10.20 865.02 4.99 1634.57 9.42
均衡-底线Scene D 326.33 1.88 3747.63 21.61 6143.04 35.42 3332.30 19.21 3796.73 21.89
均衡-满意Scene E 4896.81 28.23 5533.40 31.91 2360.28 13.61 1746.59 10.07 2808.98 16.19
均衡-理想Scene F 9443.77 54.45 2404.86 13.86 1820.77 10.49 1779.16 10.26 1897.55 10.94
低等级-底线Scene G 421.26 2.43 3818.38 22.02 5029.83 29.00 4238.71 24.44 3744.75 21.59
低等级-满意Scene H 3461.45 19.96 6718.89 38.74 2536.65 14.62 1775.25 10.24 2760.75 15.92
低等级-理想Scene I 8444.73 48.69 3553.49 20.49 1728.09 9.96 1254.40 7.23 2272.33 13.10

不同城市发展模式建设用地空间扩张均随着生态用地从底线到满意到理想水平变化呈空间集聚收敛趋势:底线水平时, 高等级源优先发展模式适宜建设区面积最高为27.90%、低等级源优先发展模式最低为21.59%;理想水平时, 高等级源优先发展模式最低为9.42%、低等级源优先发展模式最高为13.10%。相同生态安全水平时, 城市发展模式对建设用地空间扩张影响不同:底线水平时, 高等级源优先发展模式适宜建设区和生态脆弱区面积和最低为38.24%、均衡发展模式居中为41.10%、低等级源优先发展模式最高为46.03%。满意水平时, 3种模式适宜建设区和生态脆弱区面积和比较接近分别为26.33%、26.26%和26.16%。理想水平时, 高等级源优先发展模式适宜建设区和生态脆弱区面积和最低为14.41%、低等级源优先发展模式居中为20.33%、均衡发展模式最高为21.20%。生态用地总体以西部溧阳-宜兴丘陵低山地区和长荡湖、滆湖地区、南部环太湖地区及东部阳澄湖、澄湖、淀山湖地区为核心呈似圈层状集聚(表 3)。

以上说明城市发展模式对土地生态安全水平影响显著。土地生态不安全水平时, 高等级源优先发展模式对非建设用地侵占最多;随着土地生态安全水平的提高, 高等级源优先发展模式下土地集约利用的优势逐步显现[43]

4 结论与讨论

从景观单元水平过程出发, 采用不同源空间扩张互侵过程模拟方法提出一套分区形式的土地生态安全动态评价流程, 根据互侵结果绘制出分区推荐方案(图 4):(1)生态底线水平和高等级源优先发展模式下适宜建设区(图 3A红色部分)为推荐适宜建设区(面积占比27.90%), 建议采取有条件限制土地用途管制措施如建设用地增减挂钩、农地占补平衡等调整区内土地利用格局;(2)生态底线水平和3种城市化模式下生态核心区与生态安全区(图 3ADG深绿与浅绿部分)为推荐生态核心区(占比35.27%略大于生态用地34.82%现状), 建议采取严禁变更土地中类使用性质管制措施维持并提高地区土地生态安全水平;(3)生态满意水平和3种城市化模式下生态核心区与生态安全区(图 3BEH深绿与浅绿部分)为推荐生态安全区(剔除推荐生态核心区部分, 占比29.07%), 建议采取禁止变更土地大类使用性质管制措施为提高地区土地生态安全水平创造条件(与推荐生态核心区面积和达到63.89%接近生态安全满意水平最小约束值);(4)剩余部分设为缓冲区(占比7.76%), 建议采取预警机制和风险防控手段防止建设用地扩张突破缓冲区降低区域土地生态安全水平。

图 4 苏锡常地区土地生态安全分区推荐方案 Fig. 4 Recommended land ecological security zoning of Su-Xi-Chang area

采用空间扩张互侵及分区方法进行土地生态安全动态评价, 虽然把生态过程与土地利用格局变化紧密结合, 具有客观性动态性特征, 但采用不同源空间竞争性控制替代时间序列外推的方法, 评价结果以分区方式替代数值方式反映变化差异不显著。

如何考虑自然保护区、森林公园、水源地、基本农田等对建设用地扩张的刚性生态约束?基于生态风险与建设用地扩张正相关的阻力设定是否存在临界阈[31]?如何判定?如何分析不同源地类具体用途[14]?如何量化不同地区生态用地不同生态安全格局水平占地结构比差异[32]?等都影响着评价结果准确性, 亟待深入探讨。

建设用地源划分等级和引入不同城市发展模式的方法和思路, 为城市规划、土地利用规划、生态规划等空间政策制定提供了多种选择, 具有较好普适性, 对生态安全评价及格局优化等研究具有启示意义。

参考文献
[1] 肖笃宁, 陈文波, 郭福良. 论生态安全的基本概念和研究内容. 应用生态学报 , 2002, 13 (3) : 354–358.
[2] 崔胜辉, 洪华生, 黄云凤, 薛雄志. 生态安全研究进展. 生态学报 , 2005, 25 (4) : 861–868.
[3] 杨青生, 乔纪纲, 艾彬. 快速城市化地区景观生态安全时空演化过程分析--以东莞市为例. 生态学报 , 2013, 33 (4) : 1230–1239.
[4] 秦晓楠, 卢小丽, 武春友. 国内生态安全研究知识图谱--基于Citespace的计量分析. 生态学报 , 2014, 34 (13) : 3693–3703.
[5] 余敦, 陈文波. 基于物元模型的鄱阳湖生态经济区土地生态安全评价. 应用生态学报 , 2011, 22 (10) : 2681–2685.
[6] 南颖, 吉喆, 冯恒栋, 张冲冲. 基于遥感和地理信息系统的图们江地区生态安全评价. 生态学报 , 2013, 33 (15) : 4790–4798.
[7] Li Y F, Sun X, Zhu X D, Cao H H. An early warning method of landscape ecological security in rapid urbanizing coastal areas and its application in Xiamen, China. Ecological Modelling , 2010, 221 (19) : 2251–2260. DOI:10.1016/j.ecolmodel.2010.04.016
[8] Lin J Y, Lin T, Cui S H. Quantitative selection model of ecological indicators and its solving method. Ecological Indicators , 2012, 13 (1) : 294–302. DOI:10.1016/j.ecolind.2011.06.024
[9] Tian J Y, Gang G S. Research on regional ecological security assessment. Energy Procedia , 2012, 16 : 1180–1186. DOI:10.1016/j.egypro.2012.01.188
[10] Han B L, Liu H X, Wang R S. Urban ecological security assessment for cities in the Beijing-Tianjin-Hebei metropolitan region based on fuzzy and entropy methods. Ecological Modelling , 2015, 318 : 217–225. DOI:10.1016/j.ecolmodel.2014.12.015
[11] Pei L, Du L M, Yue G J. Ecological security assessment of Beijing based on PSR model. Procedia Environmental Sciences , 2010, 2 : 832–841. DOI:10.1016/j.proenv.2010.10.094
[12] Liu D, Chang Q. Ecological security research progress in China. Acta Ecologica Sinica , 2015, 35 (5) : 111–121. DOI:10.1016/j.chnaes.2015.07.001
[13] 左太安, 苏维词, 马景娜, 肖艳. 三峡重庆库区针对水土流失的土地资源生态安全评价. 水土保持学报 , 2010, 24 (2) : 74–78.
[14] 刘孝富, 舒俭民, 张林波. 最小累积阻力模型在城市土地生态适宜性评价中的应用--以厦门为例. 生态学报 , 2010, 30 (2) : 421–428.
[15] 俞孔坚. 生物保护的景观生态安全格局. 生态学报 , 1999, 19 (1) : 8–15.
[16] 李卫锋, 王仰麟, 蒋依依, 李贵才. 城市地域生态调控的空间途径--以深圳市为例. 生态学报 , 2003, 23 (9) : 1823–1831.
[17] 俞孔坚, 黄刚, 李迪华, 刘海龙. 景观网络的构建与组织--石花洞风景名胜区景观生态规划探讨. 城市规划学刊 , 2005 (3) : 76–81.
[18] 赵宏波, 马延吉. 基于变权-物元分析模型的老工业基地区域生态安全动态预警研究--以吉林省为例. 生态学报 , 2014, 34 (16) : 4720–4733.
[19] 周彬, 钟林生, 陈田, 周睿. 基于变权模型的舟山群岛生态安全预警. 应用生态学报 , 2015, 26 (6) : 1854–1862.
[20] 中国经济增长前沿课题组. 中国经济转型的结构性特征、风险与效率提升路径. 经济研究 , 2013 (10) : 4–17.
[21] 张利, 陈影, 王树涛, 门明新, 许皞. 滨海快速城市化地区土地生态安全评价与预警--以曹妃甸新区为例. 应用生态学报 , 2015, 26 (8) : 2445–2454.
[22] 袁艺, 史培军, 刘颖慧, 谢锋. 快速城市化过程中土地覆盖格局研究--以深圳市为例. 生态学报 , 2003, 23 (9) : 1832–1840.
[23] 陈燕飞, 杜鹏飞, 郑筱津, 林瑾. 基于GIS的南宁市建设用地生态适宜性评价. 清华大学学报(自然科学版) , 2006, 46 (6) : 801–804.
[24] 邬建国. 景观生态学中的十大研究论题. 生态学报 , 2004, 24 (9) : 2074–2076.
[25] 许妍, 高俊峰, 郭建科. 太湖流域生态风险评价. 生态学报 , 2013, 33 (9) : 2896–2906.
[26] 陈利顶, 傅伯杰, 赵文武. "源""汇"景观理论及其生态学意义. 生态学报 , 2006, 26 (5) : 1444–1449.
[27] 刘桂林, 张落成, 张倩. 长三角地区土地利用时空变化对生态系统服务价值的影响. 生态学报 , 2014, 34 (12) : 3311–3319.
[28] Ye Y Y, Su Y X, Zhang H O, Liu K, Wu Q T. Construction of an ecological resistance surface model and its application in urban expansion simulations. Journal of Geographical Sciences , 2015, 25 (2) : 211–224. DOI:10.1007/s11442-015-1163-1
[29] Li F, Ye Y P, Song B W, Wang R S. Evaluation of urban suitable ecological land based on the minimum cumulative resistance model:a case study from Changzhou, China. Ecological Modelling , 2015, 318 : 194–203. DOI:10.1016/j.ecolmodel.2014.09.002
[30] 毛齐正, 黄甘霖, 邬建国. 城市生态系统服务研究综述. 应用生态学报 , 2015, 26 (4) : 1023–1033.
[31] 彭建, 党威雄, 刘焱序, 宗敏丽, 胡晓旭. 景观生态风险评价研究进展与展望. 地理学报 , 2015, 70 (4) : 664–677.
[32] 俞孔坚, 王思思, 李迪华, 李春波. 北京市生态安全格局及城市增长预景. 生态学报 , 2009, 29 (3) : 1189–1204.
[33] 吴蒙, 车越, 杨凯. 基于生态系统服务价值的城市土地空间分区优化研究--以上海市宝山区为例. 资源科学 , 2013, 35 (12) : 2390–2396.
[34] 金妍, 车越, 杨凯. 基于最小累积阻力模型的江南水乡河网分区保护研究. 长江流域资源与科学 , 2013, 22 (1) : 8–14.
[35] Costanza R, d'Arge R, de Groot R, Farber S, Grasso M, Hannon B, Limburg K, Naeem S, O'Neill R V, Paruelo J, Raskin R G, Sutton P, van den Belt M. The value of the world's ecosystem services and natural capital. Nature , 1997, 387 (6630) : 253–260. DOI:10.1038/387253a0
[36] 胡和兵, 刘红玉, 郝敬锋, 安静. 流域景观结构的城市化影响与生态风险评价. 生态学报 , 2011, 31 (12) : 3432–3440.
[37] 石浩朋, 于开芹, 冯永军. 基于景观结构的城乡结合部生态风险分析--以泰安市岱岳区为例. 应用生态学报 , 2013, 24 (3) : 705–712.
[38] 苏海民, 何爱霞. 基于RS和地统计学的福州市土地利用分析. 自然资源学报 , 2010, 25 (1) : 91–99.
[39] Malekmohammadi B, Blouchi L R. Ecological risk assessment of wetland ecosystems using multi criteria decision making and geographic information system. Ecological Indicators , 2014, 41 : 133–144. DOI:10.1016/j.ecolind.2014.01.038
[40] He W, Qin N, Kong X Z, Liu W X, Wu W J, He Q S, Yang C, Jiang Y J, Wang Q M, Yang B, Xu F L. Ecological risk assessment and priority setting for typical toxic pollutants in the water from Beijing-Tianjin-Bohai area using Bayesian matbugs calculator (BMC). Ecological Indicators , 2014, 45 : 209–218. DOI:10.1016/j.ecolind.2014.04.008
[41] Xu F L, Li Y L, Wang Y, He W, Kong X Z, Qin N, Liu W X, Wu W J, Jorgensen S E. Key issues for the development and application of the species sensitivity distribution (SSD) model for ecological risk assessment. Ecological Indicators , 2015, 54 : 227–237. DOI:10.1016/j.ecolind.2015.02.001
[42] McDonald K S, Ryder D S, Tighe M. Developing best-practice bayesian belief networks in ecological risk assessments for freshwater and estuarine ecosystems:a quantitative review. Journal of Environmental Management , 2015, 154 : 190–200. DOI:10.1016/j.jenvman.2015.02.031
[43] 周翔, 陈亮, 象伟宁. 苏锡常地区建设用地扩张过程的定量分析. 应用生态学报 , 2014, 25 (5) : 1422–1430.