文章信息
- 王旭熙, 彭立, 苏春江, 徐定德, 陈田田
- WANG Xuxi, PENG Li, SU Chunjiang, XU Dingde, CHEN Tiantian.
- 基于景观生态安全格局的低丘缓坡土地资源开发利用——以四川省泸县为例
- Development and utilization of low-slope hilly land resources based on a landscape security pattern theory: a case study in Luxian County, Sichuan Province
- 生态学报[J]. 2016, 36(12): 3646-3654
- Acta Ecologica Sinica[J]. 2016, 36(12): 3646-3654
- http://dx.doi.org/10.5846/stxb201506291340
-
文章历史
- 收稿日期: 2015-06-29
- 修订日期: 2016-01-14
2. 中国科学院大学, 北京 100049
2. University of Chinese Academy of Sciences, Beijing 100049, China
随着中国城镇化的快速推进,使得建设用地需求数量迅速上升,经济社会的发展与土地资源利用管理、生态环境保护之间的矛盾日益尖锐,已经成为了制约经济发展的重要因素。合理开发低丘缓坡土地资源能为建设用地拓展空间,有效保护耕地和节约集约用地,是当前解决用地矛盾的重要途径。因此,低丘缓坡土地资源的开发具有非常重要的现实意义。目前,关于低丘缓坡土地资源的研究已日益受到关注,其研究主要集中在低丘缓坡土地资源规划以及相关政策研究、低丘缓坡土地资源利用现状、土地资源开发适宜性等方面[1-5]。由于低丘缓坡土地资源大多具有生态敏感性和生态脆弱性等特征,因而在其开发前应该优先识别与保护对维持区域生态安全具有至关重要作用的生态屏障,构建生态安全格局,从而促进区域可持续发展。
生态安全格局被认为是实现区域生态安全的基本保障和重要途径[6]。近10年,生态安全格局的研究受到我国学者高度关注,对区域生态安全格局的概念、理论基础和构建方法等方面进行了研究[7-10],马克明等[11]提出了区域生态安全格局的概念、理论基础、设计原则与方法,适应了生态系统恢复和生物多样性保护的发展需求;俞孔坚等[6]基于景观安全格局理论,运用GIS和空间分析技术,构建了北京市生态安全格局网络;蒙吉军等[12]借助多目标化模型和空间分析技术,构建了鄂尔多斯市土地利用生态安全格局。其中,景观安全格局的构建是研究的重点和难点,如何针对景观安全格局、功能和过程之间的耦合关系,在空间上构建维护整个生态环境安全不可替代的景观组分,划分不同类型的生态安全功能区,是保障生态系统服务功能,缓解当前土地资源开发与生态环境保护之间矛盾的关键问题。
泸县位于四川南部丘陵区,属于川东平行岭谷地貌区,平原面积比重小,人口和产业多集中在山前小平原,拥有良好的低丘缓坡资源开发条件。随着泸县城镇化、工业化和农业现代化同步加速推进,耕地逐年减少,根据相关数据表明泸县从1996年到2010年间,耕地面积净减少了4502.41hm2,减少2.94%,年均减少450hm2,县人均耕地0.045hm2,低于四川省平均水平(0.074hm2/人),而对新增建设用地需求量越来越大。根据近几年泸县用地需求实际情况来看,平均每年需求约400hm2,则至2020年新增建设用地总需求为3200hm2,而《四川省泸县土地利用总体规划(2006—2020年)》中规定2020年泸县尚能供给的新增建设用地规模共计2556.65hm2,不能满足需求,亟需合理开发利用低丘缓坡土地资源,缓解用地矛盾。因此,本研究以泸县为研究区域,利用GIS空间分析技术和景观安全格局理论方法,提出不同生态安全水平的低丘缓坡土地资源用地规划,为研究区低丘缓坡土地资源的合理开发提供理论依据,促进区域土地资源持续性利用。
1 研究区概况泸县位于四川南部丘陵区,属于川东平行岭谷地貌区,界于105°10′50″—105°45′30″E,28°54′40″—29°20′00″N之间,海拔215—757m,行政管辖19个镇,全县总面积1532km2,多丘陵山区,全县雨量充沛,四季分明,无霜期长,雨热同季,冬暖春早,年均气温17.8℃,年日照时数1434.2h,年降水量1142mm,属南亚热带兼季风型气候。2010年年末全县总人口108.57万人,城镇化率为30.4%。
对各类土地利用现状分析表明,泸县土地利用中存在的问题一是建设用地需求量增加,人地矛盾日益突出。泸县从2005年到2012年间,建设用地面积净增加了2749.72hm2,增加13.20%,年均增加549.94hm2。至2012年人均建设用地0.019hm2,约低于四川省平均水平(0.020hm2/人)。二是耕地后备资源不足。适宜开发的耕地后备资源不足,可开发的面积有限,开发难度较大,并且基本上只能开垦为旱地,水田后备资源稀少。
2 数据和研究方法 2.1 数据来源与处理本研究中土地利用土地利用类型分布数据的获取以卫星影像为信息源,2007年版《土地利用现状分类》(国家标准),结合泸县土地利用景观类型结构特点,将土地利用类型划分为:耕地、园地、林地、牧草地、居民点及工矿用地、交通运输用地、水域与未利用地等8种类型,最终建立研究区土地利用属性数据库。坡度和坡向则利用ArcGIS 10.0软件的分析功能由DEM生成。
2.2 低丘缓坡土地资源范围界定在地貌学上,丘陵通常被定义为海拔高度在500m以下,相对起伏在200m以下坡度较缓、连绵不断的低矮隆起高地[13]。而目前国内对于对低丘缓坡土地资源概念表述研究较少,且现有的研究中对于低丘缓坡土地资源概念表述具有比较明显的地域性,不同研究区有不同的界定。陈军等[14]指出低丘缓坡是指海拔在250m以下,坡度小于25°的平缓坡地,它通常集中分布在盆地底部河谷平原与周围山地之间的过度地带。刘卫东等[2]指出低丘缓坡土地是坡度在25°以下,面积大于2hm2,且集中连片的缓坡度,主要包括荒草地、低效林地、废弃园地等多种土地后备开发资源;王志清等[3]认为广义低丘缓坡是指海拔300m以下(包括丘陵和山地),坡度5°—14°的山坡地,狭义低丘缓坡是指海拔低于300m的丘陵区坡度,5°—14°的山坡地,两者分别扣除其中已有的建设用地、划入基本农田保护范围的土地、生态公益林用地、湿地水域面积、优质农产品的原产地保护面积和坡度25°以上的面积后,余下的面积即为广义或狭义可供建设用地面积。从上述表述中看出对于低丘缓坡土地资源的定义主要是从海拔、坡度、面积、开发利用实际等4个方面界定的。但实际上,由于研究尺度的差异,关于低丘缓坡资源的界定和大的地貌背景是密切相关的,从以上几个方面对不同的地貌区低丘缓坡进行界定也存在一定的局限性。如某一地块所处的海拔高度属于宏观条件,而对于该地块土地资源的开发利用属于微观行动,两者之间没有绝对的联系,并且中国东西海拔差异较大,对海拔的规定不适宜青藏高原、云贵高原等高海拔地区,用海拔作为低丘缓坡的界定要素存在一定的不合理性。从坡度的范围划分来看,各研究区域差异较大,其中以25°作为低丘缓坡土地资源坡度限定的上限是比较合理的,也是《中华人民共和国水土保持法》等法律法规中所规定的禁垦坡度,而对于低丘缓坡土地资源坡度下限没有统一的合理界定。
本研究依据上述相关原则和泸县低丘缓坡土地资源开发相关规定,选取海拔500m以下,坡度为6°—25°的区域为低丘缓坡(图 2)。然后,本研究从土地开发的实际情况、开发成本和经济效率等角度出发,对提取的区域进行修正,主要考虑地块大小和坡度两个因素,面积较小限制土地开发规模,降低经济效益,同时造成空间布局分散,资源不能集中分配;地形坡度大的地区,开发难度大,而且危险隐患造成的后期防护和治理成本高。最后,结合泸县实际情况(属于川南丘陵区),研究区大部分耕地都分布在2°—10°,在上面选取的低丘缓坡土地资源基础上,选取坡度为10°—25°,地块面积大于2hm2的区域作为研究区低丘缓坡的开发后备资源(图 2)。利用ArcGIS 10.0软件分析完成。
2.3 研究区景观生态安全格局分析 2.3.1 源地的确定“源”是事物或事物向外扩散的起点和基地,具有内部同质性和向四周扩张和吸引的能力[15],是指在格局与过程研究中,能促进生态过程发展的景观类型。本研究中泸县景观安全格局的功能划分是以保护生态用地为目的,尤其是对大片连块的林地和水源的保护。因此,泸县进行生态安全源地的确定,一方面借鉴相关的研究结果,另一方面,结合泸县森林资源分布特征以及泸县内河流属于长江水系,筛选出面积>20hm2的空间上具有连续性的核心林地和水域作为生态保护源地。
2.3.2 阻力面的建立生态用地间的联系是对空间水平方向的竞争性控制和覆盖过程,它是通过克服阻力来实现的,其阻力面反映了生态用地空间连通的趋势[15]。本研究采用理论地理学中表面模型——最小累计阻力模型来建立阻力面,这一模型考虑源、距离和景观介面特征3个因素,基本公式如下:
这一公式根据Knaapen等人的模型和GIS中常用的费用距离修改而来[7, 16]。其中f是一个未知的正函数,反映空间任一点的最小阻力与其到所有源的距离和景观基面特征的正相关关系[7, 12]。Dij代表从空间中某景观单元i到源j的实地距离,Ri是某一景观i的阻力值。尽管函数f一般是未知的。但(Dij×Ri)之累计值可以被认为是从源到空间某一点的某一路径的相对易达性的衡量。其中从所有源到某点阻力的最小值被用来衡量该点的易达性。模型计算主要利用ArcGIS 10.0软件的空间分析模块完成。
2.3.3 不同安全水平的景观组分识别不同安全水平的阻力面的范围和边界的确定是通过阻力面中阻力等值线突变处来确定的[17],本研究根据不同阻力等级的发展阈值,将研究区的景观生态阻力面分析图进行重分类,划分出不同安全水平的景观生态区。然后,将图层与提取的低丘缓坡土地资源斑块图层叠加,得出低丘缓坡土地资源景观安全格局,了解低丘缓坡土地资源开发现状与未来发展方向。
3 结果与分析 3.1 泸县低丘缓坡土地资源利用ArcGIS空间分析工具,得出海拔500m以下,坡度6°—25°的泸县低丘缓坡土地资源总量为872.72km2,占土地总面积的56.91%,根据低丘缓坡土地资源调查结果(表 1),耕地是泸县的主要土地资源类型,面积为552.20km2,占低丘缓坡总量的63.27%;其次是林地、居民点及工矿用地、水域、园地、牧草地、交通用地、未利用地,其面积分别占16.01%、12.24%、3.74%、3.31%、0.45%、0.09%。
土地利用类型Land use types | 修正前Before correction | 修正后After correction | ||
面积Area/km2 | 比例Proportion/% | 面积Area/km2 | 比例Proportion/% | |
耕地 Cropland | 552.20 | 63.27 | 176.54 | 53.22 |
园地 Garden plot | 28.91 | 3.31 | 14.78 | 4.46 |
林地 Forest | 139.69 | 16.01 | 91.63 | 27.62 |
牧草地 Grass land | 7.71 | 0.88 | 5.40 | 1.63 |
居民点及工矿用地Settlements and mining sites | 106.86 | 12.24 | 33.03 | 9.96 |
交通用地 Transportation area | 3.97 | 0.45 | 1.35 | 0.41 |
水域 Water | 32.63 | 3.74 | 8.61 | 2.60 |
未利用地 Unused land | 0.75 | 0.09 | 0.36 | 0.11 |
合计 Total | 872.72 | 56.91 | 331.71 | 21.65 |
从土地开发的实际情况、开发成本和经济效率等角度出发,主要考虑了地块大小和坡度两个因素,对提取的低丘缓坡区域进行修正,最终得到泸县低丘缓坡土地资源总量为331.71km2,占土地总面积的21.65%,土地资源类型以耕地为主,占低丘缓坡总量的53.22%,林地、居民点及工矿用地、园地、水域、牧草地、交通用地、未利用地的面积分别占27.62%、9.96%、4.46%、2.6%、0.41%、0.11%。
3.2 泸县景观生态安全格局 3.2.1 保护源地的确定本研究结合泸县森林资源分布特征以及泸县内河流属于长江水系,筛选出面积较大空间上具有连续性的核心林地和水域作为生态保护源地。用ArcGIS 10.0软件提取相应图层,保护“源”的总面积为77.18km2,占研究区总面积的5.04%。
3.2.2 阻力面的建立坡度和土地覆盖类型是生态保护源地向外扩散中所遇到的阻力主要来源[18]。相关研究表明[19-20],随着坡度的增加,发生水土流失的机率增加,景观保护阻力变大;土地覆被类型与源地中景观特征越接近,其对生态流的阻力就越小;而受人类干扰较大的覆被类型,其对源地内生态流的阻力就越大。因此,本文以坡度和土地覆被类型作为阻力因子。(1)坡度因子分级。根据研究区坡度大小与水土流失强度的关系,将坡度分为4级:平缓坡0°—15°;斜坡15°—25°;陡坡25°—35°;极陡坡>35°;相对阻力值大小顺序为:平缓坡 <斜坡<陡坡<极陡坡。(2)土地覆盖因子分级:本文在参照同类相关研究所采用的分级标准的基础上,结合专家经验判断各景观类型的相对阻力值大小为:林地<水域和未利用地<牧草地<耕地和园地<交通用地、居民点及工矿用地。各因子相对阻力分值的赋予通过专家打分获得(表 2)。
阻力因子Resistance factor | 权重Weight | 阻力等级Resistance level | 相对阻力系数Resistance value |
坡度Slope | 0.4 | 0°—15° | 10 |
15°—25° | 30 | ||
25°—35° | 80 | ||
>35° | 100 | ||
土地覆盖类型Land cover types | 0.6 | 林地 | 1 |
水域 未利用地 | 10 | ||
牧草地 | 30 | ||
耕地 园地 | 50 | ||
居民点及工矿用地交通用地 | 100 |
利用ArcGIS 10.0软件空间分析中的Raster Calculator工具,进行不同要素阻力值得叠加分析,得到研究区源地扩张的生态环境阻力值(图 3)。依据图 3阻力表面,利用ArcGIS 10.0空间分析中Cost distance模型,得到最小累计阻力分布图(图 4)。从图 4看出,离源地越近,阻力值越大;反之越小。
3.2.3 泸县不同安全水平的景观组分识别
根据不同阻力等级的发展阈值,利用ArcGIS 10.0空间分析将研究区的景观生态阻力面分析图进行重分类,划分出不同安全水平的景观生态区(图 5)。
(1) 低安全水平区。该区由源地组成,是维护区域生态系统平衡的核心地带,对区域水源涵养、生物多样性保护、水土保持等起重要作用。对于该区域应加强建设和完善生态廊道、生态节点,把这些区域划分为生态保护的重点区域,禁止开发。
(2) 较低安全水平区。它即是景观安全格局中保护源地恢复或扩展的缓冲区地带,又是核心源地之间联系的重要通道,对于维护源地的生态稳定性具有重要作用。该区应划分为限制开发区域,加强生态环境保护,构建生态源地保护屏障。
(3) 中等安全水平。该区域处于生态缓冲带和人类活动区的过渡地带,对人类的开发和干扰敏感度相对较低,具有一定的抗干扰能力。该区域划分为适度开发区域,在不影响源地景观稳定性的情况下,适度开发土地资源,优化土地资源配置。
(4) 高安全水平区。该区对于保护“源”的扩展相对阻力较高,对于人类的开发和干扰敏感性不高,抗外界干扰能力强。把这些区域划为优先开发区域,应协调土地利用与生态建设,形成合理的景观布局,营造人类宜居环境。
3.3 泸县低丘缓坡土地资源开发类型及用地规划将研究区景观安全格局与界定的低丘缓坡土地资源现状进行叠加,根据景观安全水平,把研究区低丘缓坡土地资源开发类型划分优先开发型、适度开发型、限制开发型、禁止开发型4种类型(图 6),各种类型的低丘缓坡土地资源面积及其比例见表 3所示。从表 3得出,泸县范围内处于高安全水平状态的低丘缓坡土地资源面积总量最大,其面积共计108.72km2,占低丘缓坡土地资源总量的32.78%,开发类型为优先开发型。位于中等安全水平的低丘缓坡土地资源面积总量次之,其面积为97.03km2,占低丘缓坡土地资源总量的29.25%,开发类型为适度开发型。位于较低安全水平的低丘缓坡土地资源面积总量为68.36km2,占低丘缓坡土地资源总量的20.60%,开发类型为限制开发型。位于低安全水平的低丘缓坡土地资源面积总量最低,其面积为57.60km2,占低丘缓坡土地资源总量的17.36%。并依据泸县土地利用、土地整治相关规划,提出各类型用地规划(图 7)。
土地资源开发类型Land resources development type | 低丘缓坡土地资源Land resource development type of low-slope hilly region | 景观安全水平Levels of Landscape security | |
面积Area/km2 | 比例Proportion/% | ||
优先开发型Priority development | 108.72 | 32.78 | 高安全水平 |
适度开发型Moderate development | 97.03 | 29.25 | 中等安全水平 |
限制开发型Limited development | 68.36 | 20.60 | 较低安全水平 |
禁止开发型Prohibited development | 57.60 | 17.36 | 低安全水平 |
(1) 优先开发型。该类型低丘缓坡土地资源主要以耕地和林地为主,未利用地较少。相对于其它3种类型低丘缓坡土地资源,居民点及工矿用地最多。因此,该类型低丘缓坡土地资源的开发利用,除了合理开发利用未利用地外,需按照泸县土地利用、整治相关规划对耕地、居民点及工矿用地进行整治。根据《四川省泸县土地利用总体规划(2006—2020年)》中针对泸县农村居民点布局分散且占地较多的现状提出的关于在福集镇、太伏镇、海潮镇、立石镇、毗邻镇等各镇实施土地综合整治工程、土地复垦工程等相关措施,并结合泸县低丘缓坡的空间分布,建议在低丘缓坡开发利用时划定出相应的土地整治区(图 7),提高低丘缓坡土地资源开发潜力。依据《泸县国民经济和社会发展第十二个五年计划纲要》提出的提升城西工业园区,打造临港经济园区的目标,建议因地制宜地对区域内低丘缓坡土地资源进行开发,进而缓解建设用地压力。
(2) 适度开发型。该类型低丘缓坡土地资源主要以农耕地为主,城镇化的快速发展势必会引起大量优质农地被占用,为了解决城镇扩张与耕地红线保护之间的矛盾,对于此类型的低丘缓坡土地资源利用,应该对零星耕地、田坎、田间道路等进行调整,提高低丘缓坡土地资源开发潜力,缓解建设用地压力,并通过不断完善农田人工廊道系统、形成人工建设区和农业耕作区交错分布的景观布局,提高农业生产的效益。按照因地制宜、合理开发的原则,稳步实施土地开发工程,加强农田基础设施建设;抓好废弃地和塌陷地复耕工程。
(3) 限制开发型和禁止开发型。该类型低丘缓坡土地资源主要分布在泸县薄刀岭、黄瓜岭、龙贯山山脉,处于地质条件脆弱或水源涵养区,对于区域生态环境保护具有重要作用。因此,该区域应限制低丘缓坡土地资源开发(特别是禁止开发型土地资源),保存原有植被,并对已经被人类改造利用的关键土地资源采取相应的生态恢复措施,构筑泸县生态屏障。另外,限制开发型低丘缓坡土地资源主要以农耕地为主,较低水平生态安全格局下进行农业活动,将存在地表水污染、水土流失等隐患,应适当缩减耕地,推进耕地退草、退林等生态退耕措施,以及建立土地复垦区。
4 结论随着中国城镇化的快速推进,建设用地与生态环境保护用地之间的矛盾日益尖锐,低丘缓坡土地资源开发对于缓解这一矛盾起着重要作用。本研究从土地开发的实际情况、开发成本和经济效率等角度出发,主要考虑地块大小和坡度两个因素,界定了研究区低丘缓坡土地资源的范围,共提取低丘缓坡土地资源面积达331.71km2。然后以生态用地保护为目的,选取土地覆盖类型与地形坡度为阻力因子,应用最小累计阻力模型建立泸县景观生态安全格局,这为研究区生态安全策略提供了重要参考。
本研究在上述分析的基础上,将研究区景观安全格局与界定的低丘缓坡土地资源现状进行叠加分析,对不同景观安全水平的低丘缓坡土地资源进行分类,将低丘缓坡土地资源类型分为优先开发型、适度开发型、限制开发型、禁止开发型4种类型,其面积分别为108.72、97.03、68.36、57.60km2,并依据泸县土地利用、土地整治相关规划,提出各类型用地规划,这为低丘缓坡土地资源合理开发提供了重要科学依据。
本研究对低丘缓坡土地资源开发安全格局的研究,主要考虑了土地资源的自然属性,而没有考虑低丘缓坡土地资源的社会经济属性,需进一步考虑土地社会经济属性如土地资源人口承载力、土地的生产潜力、土地的经济效益等,优化低丘缓坡土地资源开发格局。另外,目前有关阻力系数的确定没有一个明确的方法,在已有的相关研究中,阻力系数的确定很多都采用专家经验打分法[7, 15, 19],本研究中景观生态安全格局阻力系数及其权重的确定也是参考已有的研究成果基础上采用专家打分的方式得到。但这种赋值方式较为主观,直接影响了评价结果的客观性,因此,在今后的研究过程中还需考虑引入一些数学方法改进阻力赋值方式。
[1] | 吕杰, 袁希平, 甘淑.低丘缓坡土地资源开发利用战略分析研究.中国农学通报,2013,29(35): 225–229. |
[2] | 刘卫东, 严伟.经济发达地区的低丘缓坡土地资源合理开发利用以——浙江省永康市为例.国土资源科技管理,2007,24(3): 1–5. |
[3] | 王志清, 林飞, 苗国丽, 陈侃, 魏洋, 陈开超.低丘缓坡林地开发利用研究——以浙江省宁波市为例.林业经济,2011(12): 74–77. |
[4] | 刘焱序, 彭建, 韩忆楠, 魏海, 杜悦悦.基于OWA的低丘缓坡建设开发适宜性评价——以云南大理白族自治州为例.生态学报,2014,34(12): 3188–3197. |
[5] | 李可, 邓文胜, 刘海.基于GIS的郧县低丘缓坡土地开发建设适宜性评价.湖北大学学报:自然科学版,2014,36(3): 283–288. |
[6] | 俞孔坚, 王思思, 李迪华, 李春波.北京市生态安全格局及城市增长预景.生态学报,2009,29(3): 1189–1204. |
[7] | 李晖, 易娜, 姚文璟, 王思琪, 李志英, 杨树华.基于景观安全格局的香格里拉县生态用地规划.生态学报,2011,31(20): 5928–5936. |
[8] | 黎晓亚, 马克明, 傅伯杰, 牛树奎.区域生态安全格局:设计原则与方法.生态学报,2004,24(5): 1055–1062. |
[9] | 刘吉平, 吕宪国, 杨青, 王海霞.三江平原东北部湿地生态安全格局设计.生态学报,2009,29(3): 1083–1090. |
[10] | 刘洋, 蒙吉军, 朱利凯.区域生态安全格局研究进展.生态学报,2010,30(24): 6980–6989. |
[11] | 马克明, 傅伯杰, 黎晓亚, 关文彬.区域生态安全格局:概念与理论基础.生态学报,2004,24(4): 761–768. |
[12] | 蒙吉军, 赵春红, 刘明达.基于土地利用变化的区域生态安全评价——以鄂尔多斯市为例.自然资源学报,2011,26(4): 578–590. |
[13] | 魏海, 秦博, 彭建, 金鑫.基于GRNN模型与邻域计算的低丘缓坡综合开发适宜性评价——以乌蒙山集中连片特殊困难片区为例.地理研究,2014,33(5): 831–841. |
[14] | 陈军, 洪新岩.金华市低丘缓坡地资源合理开发利用的探讨.浙江师范大学学报:自然科学版,1989,12(2): 134–138. |
[15] | 周锐, 王新军, 苏海龙, 娄翼来.平顶山新区生态用地的识别与安全格局构建.生态学报,2015,35(6): 2003–2012. |
[16] | Knaapen J P, Scheffer M, Harms B.Estimating habitat isolation in landscape planning.Landscape and Urban Planning,1992,23(1): 1–16. |
[17] | 叶玉瑶, 苏泳娴, 张虹鸥, 吴旗韬, 刘凯.生态阻力面模型构建及其在城市扩展模拟中的应用.地理学报,2014,69(4): 485–496. |
[18] | 赵筱青, 王海波, 杨树华, 徐晓雅.基于GIS支持下的土地资源空间格局生态优化.生态学报,2009,29(9): 4892–4901. |
[19] | 刘孝富, 舒俭民, 张林波.最小累积阻力模型在城市土地生态适宜性评价中的应用——以厦门为例.生态学报,2010,30(2): 421–428. |
[20] | 文博, 刘友兆, 夏敏.基于景观安全格局的农村居民点用地布局优化.农业工程学报,2014,30(8): 181–191. |