Habitat Quality of Farming-pastoral Ecotone in Bairin Right Banner, Inner Mongolia Based on Land Use Change and InVEST Model From 2005 to 2015
-
摘要:
土地利用变化与生态效应的关联研究是全球环境变化研究领域的热点。土地利用变化下的生境质量变化研究为生态保护及土地管理政策的制定提供了重要依据。运用GIS技术和InVEST模型,基于2005-2015年遥感数据对内蒙古巴林右旗生境质量变化进行分析。结果表明:巴林右旗土地利用类型以草地、林地和耕地为主,3者面积占全旗土地总面积的85%以上。2005-2015年,湿地、建设用地和林地面积分别增加335.01、298.29和258.67 km2,增长率分别为91.80%、301.85%和12.07%,耕地、草地和未利用地面积分别减少533.98、279.28和78.71 km2,减少率分别为28.34%、5.63%和15.40%。10 a间,草地、林地、旱地和沙地转换面积较大,其中,旱地和沙地转换为重要生态用地(林地、草地、湿地)的面积为948.70和303.53 km2。巴林右旗生境质量指数表现出明显的空间分异特征,其中,2005年较低等级生境质量区域面积占比最大,而2010、2015年则均以高等级生境质量区域面积占比为最大。10 a间,高、较高和低等级生境质量区域面积均有所增加,高等级生境质量区域面积增长和较低等级生境质量区域面积减少最显著,变化率分别为39.65%和37.55%。10 a间,研究区平均生境质量指数由0.457提高到0.528,总体表现为上升趋势,表明该区域生境质量总体趋于好转,退耕还林还草工程的实施对该地区生态恢复起到了积极的促进作用。
Abstract:The ecological effects of land-use changes have emerged as a hot topic within research on global environmental changes. The findings of studies on changes in habitat quality resulting from land-use changes provide important basis for policy-making regarding ecological protection and land management. In this study, GIS technology was used in conjunction with the InVEST model to analyze remote sensing data compiled for the period 2005-2015 to assess changes in the quality of habitat in Bairin Right Banner in Inner Mongolia. The results of the analysis indicated that grassland, woodland, and farmland were the dominant types of land use in this region, accounting for more than 85% of the total land area. During the period 2005-2015, areas of wetland, construction land, and woodland increased by 335.01, 298.29, and 258.67 km2, evidencing corresponding growth rates of 91.80%, 301.85%, and 12.07%, respectively. Conversely, the areas of farmland, grassland, and unused land decreased by 533.98, 279.28, and 78.71 km2, with associated reduction rates of 28.34%, 5.63%, and 15.40%, respectively. In the decade, large areas of grassland, woodland, dryland, and sandy land had undergone conversion. Of these areas, 948.70 and 303.53 km2 of dryland and sandy land, respectively, had been converted into land-use types of ecological importance (forestland, grassland, and wetland). The habitat quality index applied in the study show distinct spatial differentiation characteristics. Whereas areas of land of lower habitat quality were predominant in 2005, areas of land of higher habitat quality were predominant in 2010 and 2015. During this 10-year period, land areas with habitat quality of high, higher, and low grades in the index increased, with marked increases in areas with high-grade habitat quality (a change rate of 39.65%) and decreases in areas of lower-grade habitat quality (a change rate of 37.55%). The average habitat quality index increased from 0.457 to 0.528, indicating an improvement in the overall habitat quality in this region. Thus, the implementation of the project of returning farmland to forest and grassland has evidently been effective, contributing to ecological recovery in this region.
-
Keywords:
- land use change /
- habitat quality /
- farming-pastoral ecotone /
- ecological restoration /
- InVEST model
-
土地利用变化研究是全球生态环境领域研究的重要内容[1-5]。土地利用除具有经济、社会、文化和生产性功能外, 还具有显著的生态效应。土地利用及其结构的变化能够深刻影响生境斑块之间的物质流、能量流循环过程, 进而改变区域生境分布格局和功能[6-9]。通过研究土地利用变化背景下生境质量的变化状况, 可以及时获取土地利用变化产生的生态响应, 从而为生态保护及土地管理政策的制定或调整提供参考[10-11]。
农牧交错带是以草地和农田大面积交错出现并以典型景观为特征的自然群落与人工群落相互镶嵌的生态复合体, 是典型的生态交错区和生态脆弱区[12]。由于生态交错带自身所具有的边缘效应, 农牧交错带对气候变化和人类干扰表现出极端敏感的特征, 成为外界干扰信号的放大器, 农牧交错带土地利用变化对研究区域生态变化状况具有重要的预警作用[13]。近年来, 随着政府和学术界对我国农牧交错带生态环境退化问题的广泛关注, 一系列生态保护和生态恢复措施相继开展, 如禁牧、休牧、轮牧等工程, 对该地区生态环境的改善具有重要的促进作用。目前, 农牧交错带生态恢复成效评估研究逐渐引起关注[14-16], 但已有的成效评估多是对植被数量特征[17]、植被指数[18]和生态环境背景状况[19]等单个指标的具体评价, 而基于土地利用变化的综合评估研究相对较少。为此, 笔者利用GIS技术和InVEST模型, 对位于北方典型农牧交错带的巴林右旗土地利用变化及其生态质量响应状况进行分析, 通过对其不同时期土地利用状况的研究, 获取生态恢复措施下的土地利用变化状况, 在此基础上, 分析生境质量响应状况, 掌握其生态恢复成效, 以期为农牧交错地区生态恢复建设与管理提供支持。
1. 研究区概况
研究区位于我国北方典型农牧交错区——内蒙古自治区赤峰市巴林右旗, 处于大兴安岭与燕山交接地带, 分布范围在43°12′~44°28′ N、118°11′~120°05′ E之间, 面积为9 959.91 km2。该区属温带半干旱大陆性季风气候区, 四季分明, 冬季寒冷漫长, 夏季短促且降雨集中, 春秋两季气温日变化剧烈, 年平均气温约为4.9 ℃, 年平均降水量约为358 mm, 年平均无霜期为125 d, 光照资源丰富, 年日照时数约为3 000~3 200 h[20]。该区自北向南依次分布有山地、丘陵、平川、沙地、滩地和水域等地貌类型。土壤类型包括黑土、灰色森林土、棕壤、栗钙土、草甸土和风沙土6个类型, 土壤肥力较差, 有机质含量少, 区域差异较大[21]。
2. 数据来源与研究方法
2.1 数据来源及处理
土地利用数据来源于中国科学院地理科学与资源研究所提供的分辨率为30 m的栅格数据集(http://www.resdc.cn/data.aspx?DATAID=99), 该数据集是目前我国精度最高的土地利用遥感监测数据产品, 已经在国家土地资源调查、水文、生态研究中发挥着重要作用。参照GB/T 21010—2007《土地利用现状分类标准》, 根据巴林右旗土地资源特征及分析要求, 对研究区土地类型进行归集处理, 各土地利用类型见表 1。
表 1 土地利用遥感解译分类系统Table 1. Land use classification system一级地类 二级地类 耕地 水田、旱地 林地 有林地、灌木林、疏林地、其他林地 草地 高覆盖度草地、中覆盖度草地、低覆盖度草地 湿地 河渠、湖泊、水库/坑塘、滩地、沼泽地 建设用地 城镇用地、农村居民点用地、工交建设用地 未利用地 沙地、盐碱地、裸地 2.2 土地利用变化计算
受不同驱动因子的作用, 研究区内不同土地利用类型的变化幅度和速度均存在差异, 采用净变化、总变化、趋势指数和净变化速率等指标计算土地利用变化动态。其中, 净变化指某一土地利用类型在一定时期内转入面积与转出面积之差; 总变化指转入面积与转出面积之和; 趋势指数描述的是某一土地利用类型面积变化的正负, 反映该土地利用类型扩张或萎缩状况; 净变化速率表征的是某一土地利用类型面积随时间变化的快慢。具体计算公式[22]为
$$ \begin{aligned} C_{\mathrm{n}} &=\left[\left(U_{\mathrm{b}}-U_{\mathrm{a}}\right) / U_{\mathrm{a}}\right] \times 100 \% \\ &=\left[\left(\Delta U_{\mathrm{in}}-\Delta U_{\mathrm{out}}\right) / U_{\mathrm{a}}\right] \times 100 \%, \end{aligned} $$ (1) $$ C_{\mathrm{t}}=\left[\left(\Delta U_{\mathrm{in}}+\Delta U_{\mathrm{out}}\right) / U_{\mathrm{a}}\right] \times 100 \%, $$ (2) $$ \begin{array}{r} S_{\mathrm{p}}=C_{\mathrm{n}} / C_{\mathrm{t}}=\left(\Delta U_{\mathrm{in}}-\Delta U_{\mathrm{out}}\right) /\left(\Delta U_{\mathrm{in}}+\right. \\ \left.\Delta U_{\mathrm{out}}\right)\left(\Delta U_{\mathrm{in}}+\Delta U_{\mathrm{out}} \neq 0 k-1 \leqslant S_{\mathrm{p}} \leqslant 1\right), \end{array} $$ (3) $$ \begin{aligned} S_{\mathrm{r}} &=(\sqrt[T]{\frac{U_{\mathrm{b}}}{U_{\mathrm{a}}}-1}) \times 100 \% \\ &=[\sqrt[T]{\frac{U_{\mathrm{a}}+\left(\Delta U_{\mathrm{in}}-\Delta U_{\mathrm{out}}\right)}{U_{\mathrm{a}}}}-1] \times 100 \%。 \end{aligned} $$ (4) 式(1)~(4)中, Cn、Ct和Sr分别为某一土地利用类型面积的净变化、总变化和净变化速率, %; Sp为趋势指数; Ua、Ub分别为研究初期和末期该土地利用类型面积, km2; ΔUout为研究时段内某土地利用类型转变为其他土地利用类型的面积之和, km2; ΔUin为同时期其他土地利用类型转变为某土地利用类型的面积之和, km2; T为研究时段, a。
利用ArcGIS 10.2软件进行相关运算, 得到2005—2015年巴林右旗土地利用转移矩阵。
2.3 基于InVEST模型的生境质量评估
生境质量指生态系统能够提供给单个物种和种群生存条件的能力, 在模型中是一个从低值到高值的连续变量。基于InVEST模型, 利用土地利用信息和生物多样性威胁因子计算生境质量指数, 可以反映土地利用动态变化对生态系统服务功能变化的影响。模型以栅格为基本评价单元, 每个栅格被分配一种土地利用类型, 综合考虑威胁因子的影响距离、空间权重等因素, 计算每个栅格的生境退化指数, 进而得到生境质量指数。土地利用类型j中栅格x的生境退化指数计算公式[9, 23]为
$$ D_{x j}=\sum\limits_{r=1}^{R} \sum\limits_{y=1}^{Y_{r}}\left(w_{r} / \sum\limits_{r=1}^{R} w_{r}\right) r_{y} i_{r x y} \beta_{x} S_{j r}, $$ (5) $$ i_{r x y}=1-\left(d_{x y} / d_{r, \max }\right)(线性衰退), $$ (6) $$ i_{r x y}=\exp \left[-\left(\frac{2.99}{d_{r, \max }}\right) d_{x y}\right](指数衰退)。 $$ (7) 式(5)~(7)中, Dxj为生境退化程度; R为胁迫因子个数; wr为胁迫因子r的权重; Yr为胁迫因子层在土地利用类型图中的栅格个数; ry为土地利用类型图中栅格y中胁迫因子个数; Sjr为j类型土地覆被对胁迫因子r的敏感性; irxy为栅格y中的威胁因子r对栅格x的影响; dxy为栅格x(生境)与栅格y(威胁因子)之间的距离; dr, max为威胁因子r的影响范围; βx为法律保护程度, 该研究不考虑法律保护程度因子, 将βx设为1。
$$ Q_{x j}=H_{j}\left(1-\frac{D_{x j}^{2}}{D_{x j}^{2}+k^{2}}\right)。 $$ (8) 式(8)中, Qxj为土地利用类型图j中栅格x的生境质量; Hj为土地利用类型图j的生境属性; k为半饱和常数, 一般设置为生境退化程度最大值的1/2[24]。依据InVEST模型使用说明, z为归一化常量, 通常其取值为2.5, 生境质量得分值域为0~1(生境适宜度最高时生境质量得分为1, 最低时为0)[23-24]。
通过InVEST模型运算及ArcGIS 10.2软件叠加分析等过程, 得到巴林右旗2005—2015年生境质量分布及变化状况。
2.4 模型参数的确定
InVEST模型中需要根据研究区状况进行调整的参数主要包括威胁因子的最大影响距离及权重、各土地利用类型的生境质量及其对威胁因子的敏感程度。
生态威胁因子对生境质量影响的大小通过空间距离变量来测算。参考已有相关研究[20, 25-27]及专家咨询结果, 确定农村居民点用地、城镇用地、工交建设用地、耕地和盐碱地5种土地利用类型作为巴林右旗生境质量变化的生态威胁因子, 对研究区生态威胁因子参数的设置见表 2。
表 2 威胁因子主要调整参数Table 2. Adjustment parameters of threat factors威胁因子 最大影响距离/km 权重 空间衰退类型 农村居民点用地 5 0.8 指数 城镇用地 8 1.0 指数 工交建设用地 4 0.9 指数 耕地 3 0.7 线性 盐碱地 1 0.1 线性 不同土地覆盖类型对于生态威胁因子具有不同的敏感度。敏感度取值范围为0~1, 0表示对某威胁因子无敏感度, 1表示对某威胁因子具有较高敏感度。参照模型自带案例和相关研究成果[28-31], 各土地利用类型对生态威胁因子敏感度的设置见表 3。
表 3 各土地利用类型对生态威胁因子的敏感度Table 3. Sensitivity of the land use types to habitat threat factor土地利用类型 生境质量得分 对生态威胁因子敏感度 居民点 城镇 建设 耕地 盐碱地 水田 0.5 0.3 0.4 0.4 0.3 0.1 旱地 0.5 0.3 0.4 0.4 0.3 0.1 有林地 0.9 0.7 0.8 0.8 0.6 0.2 灌木林 0.8 0.6 0.7 0.7 0.5 0.1 疏林地 0.3 0.4 0.5 0.55 0.4 0.1 其他林地 0.1 0.2 0.2 0.45 0.2 0.1 高覆盖 0.8 0.6 0.65 0.75 0.4 0.2 中覆盖 0.5 0.5 0.55 0.6 0.3 0.2 低覆盖 0.3 0.4 0.5 0.55 0.2 0.2 河渠 0.8 0.7 0.8 0.9 0.65 0.1 湖泊 0.9 0.7 0.8 0.9 0.65 0.1 水库/坑塘 0.8 0.7 0.8 0.9 0.65 0.1 滩地 0.5 0.2 0.3 0.3 0.2 0.1 沼泽地 0.9 0.6 0.7 0.7 0.65 0.2 城镇 0 0 0 0 0 0 居民点 0 0 0 0 0 0 建设 0 0 0 0 0 0 沙地 0 0 0 0 0 0 盐碱地 0.2 0.2 0.3 0.3 0.1 0 裸地 0 0 0 0 0 0 高、中、低覆盖分别指高、中、低覆盖度草地,居民点指农村居民点用地,建设指工交建设用地,城镇指城镇用地。 在ArcGIS 10.2软件空间分析功能支持下, 对生境质量栅格图进行运算, 得到研究区2005—2015年生境质量变化栅格图。对输出栅格图进行重分类, 划分为低、较低、中等、较高和高共5个等级, 制定生境质量分级标准(表 4)。在生境质量等级划分基础上, 得到研究区生境质量的空间分布及其变化情况。
表 4 巴林右旗生境质量等级划分Table 4. Habitat quality hierarchy in Bairin Right Banner等级 取值范围 等级特征描述 低 0~0.15 生境质量低, 生物多样性低 较低 >0.15~0.35 生境质量较低, 生物多样性较低 中等 >0.35~0.60 生境质量中等, 生物多样性中等 较高 >0.60~0.75 生境质量较高, 生物多样性较高 高 >0.75 生境质量高, 生物多样性高 3. 结果与分析
3.1 巴林右旗土地利用变化
2005—2015年, 巴林右旗土地利用现状及变化见表 5~6。由表 5~6可知, 巴林右旗土地利用类型以草地、林地和耕地为主, 2015年巴林右旗土地利用类型以草地、林地和耕地为主, 面积分别为4 679.51、2 401.00和1 350.03 km2, 分别占全旗土地总面积的46.98%、24.11%和13.55%。
表 5 2005—2015年巴林右旗土地利用现状Table 5. Land use types in Bairin Right Banner in the period of 2005-2015km2 地类 2005年 2010年 2015年 水田 15.72 54.83 7.28 旱地 1 868.29 1 436.14 1 342.75 有林地 782.79 1 020.09 924.23 灌木林 810.26 903.54 1 365.92 疏林地 518.06 165.38 72.60 其他林地 31.22 59.69 38.25 高覆盖度草地 1 281.46 3 215.62 2 143.56 中覆盖度草地 2 363.98 1 811.28 1 777.36 低覆盖度草地 1 313.35 452.56 758.59 河渠 203.05 204.72 172.79 湖泊 36.48 39.02 34.41 水库/坑塘 7.93 7.21 27.20 滩地 61.57 53.14 277.17 沼泽地 55.89 171.38 188.36 城镇用地 7.69 9.03 8.21 农村居民点用地 66.12 88.86 139.56 工交建设用地 25.01 80.13 249.34 沙地 506.22 154.56 382.16 盐碱地 1.44 23.69 38.28 裸地 3.38 9.04 11.89 总计 9 959.91 9 959.91 9 959.91 表 6 巴林右旗土地利用变化指数Table 6. Index of land use change in Bairin Right Banner地类 2005—2010年 2010—2015年 2005—2015年 A/% B/% C D/% A/% B/% C D/% A/% B/% C D/% 水田 248.79 450.00 0.56 28.47 -86.72 108.05 -0.80 -33.09 -53.69 146.94 -0.36 -14.04 旱地 -23.13 91.83 -0.25 -5.12 -6.50 76.45 -0.09 -1.35 -28.13 96.78 -0.29 -6.40 有林地 30.31 136.83 0.22 5.46 -9.40 107.02 -0.09 -2.05 18.07 141.76 0.12 3.28 灌木林 11.51 174.24 0.07 2.18 51.17 193.67 0.27 8.70 68.58 229.44 0.30 11.05 疏林地 -68.08 122.65 -0.56 -20.43 -56.10 143.29 -0.39 -15.37 -85.99 113.30 -0.76 -32.65 其他林地 91.19 285.85 0.32 13.84 -35.92 113.95 -0.31 -8.41 22.52 208.00 0.11 4.12 高覆盖 150.93 274.45 0.55 20.22 -33.34 102.36 -0.33 -7.82 67.27 218.44 0.31 10.85 中覆盖 -23.38 125.80 -0.19 -5.19 -1.87 140.26 -0.01 -0.35 -24.81 129.21 -0.19 -5.52 低覆盖 -65.54 108.73 -0.60 -19.19 67.62 220.28 0.31 10.93 -42.24 131.32 -0.32 -10.36 河渠 0.82 100.93 0.00 0.01 -15.60 101.33 -0.16 -3.46 -14.90 114.20 -0.14 -3.44 湖泊 6.96 85.95 0.07 1.17 -11.81 121.00 -0.09 -2.26 -5.67 142.63 -0.04 -1.11 水库/坑塘 -9.08 189.61 -0.04 -1.46 277.25 461.39 0.59 29.98 243.00 435.74 0.56 28.08 滩地 -13.69 151.84 -0.09 -2.99 421.58 605.29 0.75 40.85 350.17 539.61 0.70 36.65 沼泽地 206.64 352.61 0.59 25.19 9.91 154.53 0.06 1.75 237.02 404.76 0.58 27.34 城镇用地 17.43 49.81 0.35 3.26 -9.08 83.33 -0.11 -1.87 6.76 115.19 0.06 1.33 居民点 34.39 145.43 0.24 6.13 57.06 132.04 0.43 9.40 111.07 215.75 0.51 16.09 建设用地 220.39 295.33 0.75 26.22 211.17 366.87 0.58 25.51 896.96 1 056.12 0.85 58.44 沙地 -69.47 94.46 -0.73 -21.10 147.26 231.55 0.64 19.94 -24.51 115.58 -0.21 -5.44 盐碱地 1 545.14 1 675.00 0.88 73.56 61.59 217.63 0.29 10.33 2 558.33 2 667.57 0.93 91.48 裸地 167.46 364.62 0.46 21.81 31.53 146.13 0.19 5.05 251.78 442.98 0.55 27.95 A为净变化, B为总变化, C为趋势指数, D为净变化速率。高、中、低覆盖分别指高、中、低覆盖度草地, 居民点指农村居民点用地, 建设用地指工交建设用地。 2005—2010年, 巴林右旗土地利用/覆盖类型中, 草地、湿地、建设用地和林地面积有所增加, 耕地和未利用地面积有所减少。5 a间, 以草地面积增长最为显著, 净增加520.67 km2, 增长率为10.50%, 其中以高覆盖度草地扩张最为突出, 净变化率为150.93%, 其后依次为湿地和建设用地, 面积分别净增加110.55和79.20 km2, 增长率分别为30.29%和80.15%, 其中以沼泽地和工交建设用地扩张较为突出, 净变化率分别为206.64%和220.39%。5 a间, 耕地和未利用地面积分别减少393.04和323.75 km2, 减少率分别为20.86%和63.35%, 其中以旱地和沙地萎缩较为突出, 净变化率分别为23.13%和69.47%。
2010—2015年, 林地、未利用地、建设用地和湿地面积有所增加, 草地和耕地面积有所减少。5 a间, 以林地面积增长最为显著, 净增加252.3 km2, 增长率为11.74%, 其中以灌木林扩张最为突出, 净变化率为51.17%;其后依次为未利用地和建设用地, 面积分别净增加245.04和219.09 km2, 增长率分别为130.83%和123.07%, 其中以沙地和工交建设用地扩张较为突出, 净变化率分别为147.26%和211.17%。5 a间, 草地和耕地面积分别减少799.95和140.94 km2, 减少率分别为14.60%和9.45%, 其中以高覆盖度草地和水田萎缩较为突出, 净变化率分别为-33.34%和-86.72%。
2005—2015年, 湿地、林地和建设用地面积有所增加, 耕地、草地和未利用地面积有所减少。10 a间, 以湿地面积增长最为显著, 净增加335.01 km2, 增长率为91.80%, 其中以滩地扩张最为突出, 净变化率为350.17%;其后依次为建设用地和林地, 面积分别净增加298.29和258.67 km2, 增长率分别为301.85%和12.07%, 其中以工交建设用地和灌木林扩张较为突出, 净变化率分别为896.96%和68.58%。10 a间, 耕地、草地和未利用地面积分别减少533.98、279.28和78.71 km2, 减少率分别为28.34%、5.63%和15.40%, 其中, 以水田、低覆盖度草地和沙地萎缩较为突出, 净变化率分别为-53.69%、-42.24%和-24.51%。10 a间, 林地、湿地和建设用地面积持续增加, 耕地面积则持续减少。其中, 林地、沼泽地、农村居民点用地、工交建设用地、盐碱地和裸地均表现为持续扩张趋势, 旱地、疏林地和中覆盖度草地则表现为持续萎缩趋势。
由表 7可知, 草地、林地、旱地和沙地转换面积较大, 旱地和沙地转换为林地、草地和湿地这类生态用地总面积分别为948.70和303.53 km2, 说明退耕还林还草工程的实施对于该地区生态恢复起到积极的促进作用。由表 7可知, 中覆盖度草地转换为高覆盖度草地最为明显, 转换面积为660.42 km2, 其后依次为中覆盖度草地转换为灌木林和旱地转换为高覆盖度草地, 转换面积分别为327.12和322.47 km2, 转换明显的各地类基本表现为趋好方向发展; 但是也存在较明显的中覆盖度草地转换为旱地的情况, 表明部分地区仍然存在人为活动开发耕地的现象。
表 7 2005—2015年巴林右旗土地利用类型转移矩阵Table 7. Land use types transfer in 2005-2015 in Bairin Right Bannerkm2 类型 水田 旱地 有林地 灌木林 疏林地 其他林地 高覆盖度草地 中覆盖度草地 低覆盖度草地 河渠 湖泊 水库/坑塘 滩地 沼泽地 城镇用地 农村居民点用地 工交建设用地 沙地 盐碱地 裸地 总计 水田 0 2.69 1.23 0.11 0 0 0.47 0.01 0.05 0.06 0 0 0.09 0.16 0 0.62 10.19 0 0 0 15.68 旱地 1.27 0 78.31 193.75 3.99 16.75 322.47 232.29 84.79 14.15 2.20 1.14 38.43 32.56 3.66 45.23 80.06 12.18 2.17 0.62 1 166.11 有林地 0.53 47.03 0 157.50 0.26 1.54 116.20 92.86 27.66 4.14 4.62 1.91 8.28 13.28 0.07 2.07 4.79 2.19 0.29 0.01 485.23 灌木林 4.16 62.06 111.80 0 4.77 0.11 186.50 149.10 55.48 6.45 1.52 0.92 14.74 12.57 0.02 6.47 14.35 16.71 1.05 1.02 649.80 疏林地 0.07 18.38 141.40 100.80 0 0.51 106.60 87.69 28.86 1.03 0.02 0.21 7.40 5.96 0 2.35 5.36 9.30 0.01 0 515.95 其他林地 0 7.11 1.50 3.45 1.20 0 6.82 4.20 2.36 0 0 0 0.05 0.59 0 0.32 1.35 0.03 0 0.03 29.01 高覆盖度草地 0.01 159.40 134.80 188.80 2.82 6.79 0 262.97 83.30 9.43 4.58 1.61 35.45 17.32 0.14 9.37 39.87 6.85 1.36 0.27 965.17 中覆盖度草地 0.37 211.15 93.34 327.12 15.51 7.56 660.42 0 227.25 19.69 2.10 1.28 52.23 22.19 0 21.48 63.05 72.78 14.92 5.37 1 817.87 低覆盖度草地 0.15 86.09 37.95 126.80 27.26 2.21 316.06 300.83 0 15.29 1.83 1.11 40.69 40.77 0.14 12.24 11.43 103.20 10.20 3.44 1 137.70 河渠 0.21 14.30 6.23 11.93 0.38 0.21 14.95 8.48 5.09 0 2.58 0.01 32.30 15.77 0.16 1.21 4.56 3.23 3.59 0.02 125.21 湖泊 0.33 1.53 0.04 0.16 0 0 1.95 1.68 0.79 0.17 0 17.03 0.28 1.48 0 0 0.04 0.90 0.59 0 26.97 水库/坑塘 0 1.58 0.57 0.62 0.10 0 0.83 0.80 0.66 0.16 0.71 0 0.25 0.75 0 0 0.36 0.04 0.11 0 7.54 滩地 0.03 3.62 1.50 3.75 0.18 0.01 3.42 3.27 2.73 16.86 1.24 0.17 0 7.18 0.02 0.31 0.78 1.31 0.79 0.01 47.18 城镇用地 0 0.44 0.07 0.07 0 0 0.12 0 0 0 0 0 0 0 0 3.05 0.46 0 0 0 4.21 农村居民点用地 0 7.04 2.07 4.57 0.03 0.13 5.45 4.59 3.11 0.62 0.05 0.05 1.57 1.71 0.50 0 2.47 0.50 0.06 0.02 34.54 工交建设用地 0.03 3.86 0.80 3.54 0.08 0.06 5.56 2.50 0.93 0.28 0.06 0 0.89 0.37 0 0.45 0 0.26 0.01 0 19.68 沙地 0 6.17 7.43 73.82 13.57 0.07 71.59 77.18 56.62 2.57 0.68 0.34 31.35 4.81 0 1.53 2.52 0 2.48 0.84 353.57 盐碱地 0 0.73 0 0 0 0 0.04 0.01 0.22 0 0 0 0 0 0 0.43 0 0 0 0 1.43 沼泽地 0.20 6.56 3.02 9.96 0.09 0.06 6.52 4.79 4.32 3.27 2.76 1.06 0.80 0 0 0.18 2.44 0.72 0.34 0.08 47.17 裸地 0 0.74 0.08 0.26 0 0 0.22 0.56 0.83 0.14 0.04 0 0.13 0 0.03 0.24 0.06 0.01 0.08 0 3.42 3.2 巴林右旗生境质量时空演变
2005—2015年巴林右旗生境质量分布及变化状况如图 1、表 8所示, 由此可知, 2005、2010和2015年研究区生境质量指数表现出明显的空间分异特征, 高、中等和较低等级生境质量区域面积所占比例较大, 2005、2010和2015年土地覆盖面积分别占总面积的87.24%、80.73%和78.70%。2005年, 低、较低等级生境质量区域面积占比很大, 占总面积的43.46%, 而2010、2015年占比明显降低, 分别为31.54%和36.52%, 高、较高等级生境质量区域面积占比明显提高, 由2005年的28.16%提高到46.24%(2010年)和39.92%(2015年), 表明研究区生境质量呈明显好转趋势。
表 8 2005—2015年巴林右旗不同等级生境质量区域面积变化Table 8. Area change of habitat quality index in 2005-2015 in Bairin Right Banner生境质量等级 生境质量指数值 2005—2010年 2010—2015年 2005—2015年 面积变化/km2 变化率/% 面积变化/km2 变化率/% 面积变化/km2 变化率/% 低 0~0.15 232.18 30.66 355.37 35.92 587.55 77.60 较低 >0.15~0.35 -1 519.55 -41.39 141.05 6.55 -1 378.50 -37.55 中等 >0.35~0.60 -512.83 -18.81 133.01 6.01 -379.82 -13.93 较高 >0.60~0.75 416.27 80.96 -153.91 -16.54 262.36 51.03 高 >0.75 1 383.94 60.40 -475.54 -12.94 908.40 39.65 巴林右旗生境质量空间分布的总体特征表现为北部高、南部低, 生境质量高值区主要集中在研究区北部的索博日嘎镇、查干沐沦镇和幸福之路苏木等地, 低值区集中在南部的查干诺尔镇和西拉沐沦镇。分析原因, 可能是由于北部地区主要为牧业, 生产方式对生境质量影响相对较小, 同时, 赛罕乌拉国家级自然保护区也位于其中, 对该地区生境质量提升有较好的促进作用; 而南部地区主要为农业种植业, 生境质量更容易受到人类生产生活方式的影响。
通过计算可知, 10 a间, 高、较高和低等级生境质量区域面积均有所增加;较低等级生境质量区域面积减少最显著, 净减少1 378.50 km2, 减少率为37.55%;高等级生境质量区域面积增长最显著, 净增长908.40 km2, 增长率为39.65%。巴林右旗2005、2010和2015年平均生境质量指数分别为0.457、0.577和0.528, 10 a间, 平均生境质量总体表现为上升趋势, 表明研究时段内该区域生境质量总体趋于好转。
由图 1可知, 10 a间, 巴林右旗生境质量明显降低的区域主要集中于北部的索博日嘎镇和南部的西拉沐沦镇, 生境质量指数明显升高区域则集中于查干诺尔镇、查干沐沦镇和大板镇。
4. 讨论
长期以来, 受历史自然条件和人为因素的影响, 西部地区生态条件恶劣, 生态环境质量低下。为此, 我国实施了天然林保护工程、退耕还林还草工程、京津风沙源治理专项工程和草原三牧(禁牧、休牧、轮牧)工程等一系列生态保护与恢复工程。随着各项工程及后续工程的阶段推进, 生态恢复成效评估也提上日程。土地利用变化是生态恢复最直观的变化结果, 由此产生的生境质量变化反映了生态恢复的目的和效果。基于土地利用变化的生境质量评估为生态恢复成效评估提供了途径。
退耕还林还草、围封禁牧和季节性休牧是内蒙古农牧交错区实施的主要生态工程类型[32]。自2005年以来, 巴林右旗基本建立了全旗的禁牧、休牧制度, 对自然资源生态功能的恢复起到了促进和示范作用[33]。2005—2015年的10 a间, 由于实施了退耕还林还草等生态工程措施, 巴林右旗旱地向草地和林地的转换效果显著, 耕地面积持续减少, 湿地、林地面积显著增加, 生态环境质量总体表现为好转趋势, 生态恢复取得了一定成效。
InVEST模型弥补了传统生态系统服务评估方法的不足, 为生态系统服务评估提供了有效途径。近年来, 作为InVEST模型功能模块之一的生境质量评估得到了较多应用。然而, 在生境质量评估中, 该模型仍存在一定不足, 如模型通过累加胁迫因子对生境质量的影响来获取研究区生境质量, 但胁迫因子的简单累加并不完全等同于各胁迫因子的综合影响作用, 所以模型在使用过程中仍需进一步改进和完善[24]。此外, 农牧交错带属于生态交错带地区, 受外界因素影响会产生显著的边缘效应, 与一般区域生态系统相比, 相同威胁因子可能会产生加倍效应, 如何进一步科学合理地确定生态交错带生态系统对生境威胁因子的敏感度, 将是需要进一步研究的重要内容。
5. 结论
(1) 2015年, 巴林右旗土地利用类型以草地、林地和耕地为主, 面积分别为4 679.51、2 401.00和1 350.03 km2, 分别占全旗土地总面积的46.98%、24.11%和13.55%。
(2) 2005—2015年, 湿地、林地和建设用地面积持续增加, 其中以湿地面积增长最为显著, 增长率为91.80%, 耕地、草地和未利用地面积有所减少, 特别是耕地面积持续减少。旱地向草地和林地的转换较显著, 说明退耕还草、退耕还林等工程的实施对于该地区草地生态系统的恢复起到了一定促进作用。
(3) 2005—2015年, 生境质量指数表现出明显的空间分异特征, 2005年, 较低等级生境质量区域面积占比最大, 而2010、2015年则均以高等级生境质量区域面积占比为最大; 生境质量空间分布的总体特征表现为北部高、南部低, 10 a间, 巴林右旗生境质量明显降低的区域主要集中于北部的索博日嘎镇和南部的西拉沐沦镇, 生境质量指数明显升高区域则集中于查干诺尔镇、查干沐沦镇和大板镇。
(4) 10 a间, 高、较高和低等级生境质量指数面积均有所增加, 其中, 较低等级生境质量指数面积减少最显著, 减少率为37.55%;高等级生境质量指数面积增长最显著, 增长率为39.65%。10 a间, 平均生境质量指数由0.457提高到0.528, 总体表现为上升趋势, 表明该区域生境质量总体趋于好转。
-
表 1 土地利用遥感解译分类系统
Table 1 Land use classification system
一级地类 二级地类 耕地 水田、旱地 林地 有林地、灌木林、疏林地、其他林地 草地 高覆盖度草地、中覆盖度草地、低覆盖度草地 湿地 河渠、湖泊、水库/坑塘、滩地、沼泽地 建设用地 城镇用地、农村居民点用地、工交建设用地 未利用地 沙地、盐碱地、裸地 表 2 威胁因子主要调整参数
Table 2 Adjustment parameters of threat factors
威胁因子 最大影响距离/km 权重 空间衰退类型 农村居民点用地 5 0.8 指数 城镇用地 8 1.0 指数 工交建设用地 4 0.9 指数 耕地 3 0.7 线性 盐碱地 1 0.1 线性 表 3 各土地利用类型对生态威胁因子的敏感度
Table 3 Sensitivity of the land use types to habitat threat factor
土地利用类型 生境质量得分 对生态威胁因子敏感度 居民点 城镇 建设 耕地 盐碱地 水田 0.5 0.3 0.4 0.4 0.3 0.1 旱地 0.5 0.3 0.4 0.4 0.3 0.1 有林地 0.9 0.7 0.8 0.8 0.6 0.2 灌木林 0.8 0.6 0.7 0.7 0.5 0.1 疏林地 0.3 0.4 0.5 0.55 0.4 0.1 其他林地 0.1 0.2 0.2 0.45 0.2 0.1 高覆盖 0.8 0.6 0.65 0.75 0.4 0.2 中覆盖 0.5 0.5 0.55 0.6 0.3 0.2 低覆盖 0.3 0.4 0.5 0.55 0.2 0.2 河渠 0.8 0.7 0.8 0.9 0.65 0.1 湖泊 0.9 0.7 0.8 0.9 0.65 0.1 水库/坑塘 0.8 0.7 0.8 0.9 0.65 0.1 滩地 0.5 0.2 0.3 0.3 0.2 0.1 沼泽地 0.9 0.6 0.7 0.7 0.65 0.2 城镇 0 0 0 0 0 0 居民点 0 0 0 0 0 0 建设 0 0 0 0 0 0 沙地 0 0 0 0 0 0 盐碱地 0.2 0.2 0.3 0.3 0.1 0 裸地 0 0 0 0 0 0 高、中、低覆盖分别指高、中、低覆盖度草地,居民点指农村居民点用地,建设指工交建设用地,城镇指城镇用地。 表 4 巴林右旗生境质量等级划分
Table 4 Habitat quality hierarchy in Bairin Right Banner
等级 取值范围 等级特征描述 低 0~0.15 生境质量低, 生物多样性低 较低 >0.15~0.35 生境质量较低, 生物多样性较低 中等 >0.35~0.60 生境质量中等, 生物多样性中等 较高 >0.60~0.75 生境质量较高, 生物多样性较高 高 >0.75 生境质量高, 生物多样性高 表 5 2005—2015年巴林右旗土地利用现状
Table 5 Land use types in Bairin Right Banner in the period of 2005-2015
km2 地类 2005年 2010年 2015年 水田 15.72 54.83 7.28 旱地 1 868.29 1 436.14 1 342.75 有林地 782.79 1 020.09 924.23 灌木林 810.26 903.54 1 365.92 疏林地 518.06 165.38 72.60 其他林地 31.22 59.69 38.25 高覆盖度草地 1 281.46 3 215.62 2 143.56 中覆盖度草地 2 363.98 1 811.28 1 777.36 低覆盖度草地 1 313.35 452.56 758.59 河渠 203.05 204.72 172.79 湖泊 36.48 39.02 34.41 水库/坑塘 7.93 7.21 27.20 滩地 61.57 53.14 277.17 沼泽地 55.89 171.38 188.36 城镇用地 7.69 9.03 8.21 农村居民点用地 66.12 88.86 139.56 工交建设用地 25.01 80.13 249.34 沙地 506.22 154.56 382.16 盐碱地 1.44 23.69 38.28 裸地 3.38 9.04 11.89 总计 9 959.91 9 959.91 9 959.91 表 6 巴林右旗土地利用变化指数
Table 6 Index of land use change in Bairin Right Banner
地类 2005—2010年 2010—2015年 2005—2015年 A/% B/% C D/% A/% B/% C D/% A/% B/% C D/% 水田 248.79 450.00 0.56 28.47 -86.72 108.05 -0.80 -33.09 -53.69 146.94 -0.36 -14.04 旱地 -23.13 91.83 -0.25 -5.12 -6.50 76.45 -0.09 -1.35 -28.13 96.78 -0.29 -6.40 有林地 30.31 136.83 0.22 5.46 -9.40 107.02 -0.09 -2.05 18.07 141.76 0.12 3.28 灌木林 11.51 174.24 0.07 2.18 51.17 193.67 0.27 8.70 68.58 229.44 0.30 11.05 疏林地 -68.08 122.65 -0.56 -20.43 -56.10 143.29 -0.39 -15.37 -85.99 113.30 -0.76 -32.65 其他林地 91.19 285.85 0.32 13.84 -35.92 113.95 -0.31 -8.41 22.52 208.00 0.11 4.12 高覆盖 150.93 274.45 0.55 20.22 -33.34 102.36 -0.33 -7.82 67.27 218.44 0.31 10.85 中覆盖 -23.38 125.80 -0.19 -5.19 -1.87 140.26 -0.01 -0.35 -24.81 129.21 -0.19 -5.52 低覆盖 -65.54 108.73 -0.60 -19.19 67.62 220.28 0.31 10.93 -42.24 131.32 -0.32 -10.36 河渠 0.82 100.93 0.00 0.01 -15.60 101.33 -0.16 -3.46 -14.90 114.20 -0.14 -3.44 湖泊 6.96 85.95 0.07 1.17 -11.81 121.00 -0.09 -2.26 -5.67 142.63 -0.04 -1.11 水库/坑塘 -9.08 189.61 -0.04 -1.46 277.25 461.39 0.59 29.98 243.00 435.74 0.56 28.08 滩地 -13.69 151.84 -0.09 -2.99 421.58 605.29 0.75 40.85 350.17 539.61 0.70 36.65 沼泽地 206.64 352.61 0.59 25.19 9.91 154.53 0.06 1.75 237.02 404.76 0.58 27.34 城镇用地 17.43 49.81 0.35 3.26 -9.08 83.33 -0.11 -1.87 6.76 115.19 0.06 1.33 居民点 34.39 145.43 0.24 6.13 57.06 132.04 0.43 9.40 111.07 215.75 0.51 16.09 建设用地 220.39 295.33 0.75 26.22 211.17 366.87 0.58 25.51 896.96 1 056.12 0.85 58.44 沙地 -69.47 94.46 -0.73 -21.10 147.26 231.55 0.64 19.94 -24.51 115.58 -0.21 -5.44 盐碱地 1 545.14 1 675.00 0.88 73.56 61.59 217.63 0.29 10.33 2 558.33 2 667.57 0.93 91.48 裸地 167.46 364.62 0.46 21.81 31.53 146.13 0.19 5.05 251.78 442.98 0.55 27.95 A为净变化, B为总变化, C为趋势指数, D为净变化速率。高、中、低覆盖分别指高、中、低覆盖度草地, 居民点指农村居民点用地, 建设用地指工交建设用地。 表 7 2005—2015年巴林右旗土地利用类型转移矩阵
Table 7 Land use types transfer in 2005-2015 in Bairin Right Banner
km2 类型 水田 旱地 有林地 灌木林 疏林地 其他林地 高覆盖度草地 中覆盖度草地 低覆盖度草地 河渠 湖泊 水库/坑塘 滩地 沼泽地 城镇用地 农村居民点用地 工交建设用地 沙地 盐碱地 裸地 总计 水田 0 2.69 1.23 0.11 0 0 0.47 0.01 0.05 0.06 0 0 0.09 0.16 0 0.62 10.19 0 0 0 15.68 旱地 1.27 0 78.31 193.75 3.99 16.75 322.47 232.29 84.79 14.15 2.20 1.14 38.43 32.56 3.66 45.23 80.06 12.18 2.17 0.62 1 166.11 有林地 0.53 47.03 0 157.50 0.26 1.54 116.20 92.86 27.66 4.14 4.62 1.91 8.28 13.28 0.07 2.07 4.79 2.19 0.29 0.01 485.23 灌木林 4.16 62.06 111.80 0 4.77 0.11 186.50 149.10 55.48 6.45 1.52 0.92 14.74 12.57 0.02 6.47 14.35 16.71 1.05 1.02 649.80 疏林地 0.07 18.38 141.40 100.80 0 0.51 106.60 87.69 28.86 1.03 0.02 0.21 7.40 5.96 0 2.35 5.36 9.30 0.01 0 515.95 其他林地 0 7.11 1.50 3.45 1.20 0 6.82 4.20 2.36 0 0 0 0.05 0.59 0 0.32 1.35 0.03 0 0.03 29.01 高覆盖度草地 0.01 159.40 134.80 188.80 2.82 6.79 0 262.97 83.30 9.43 4.58 1.61 35.45 17.32 0.14 9.37 39.87 6.85 1.36 0.27 965.17 中覆盖度草地 0.37 211.15 93.34 327.12 15.51 7.56 660.42 0 227.25 19.69 2.10 1.28 52.23 22.19 0 21.48 63.05 72.78 14.92 5.37 1 817.87 低覆盖度草地 0.15 86.09 37.95 126.80 27.26 2.21 316.06 300.83 0 15.29 1.83 1.11 40.69 40.77 0.14 12.24 11.43 103.20 10.20 3.44 1 137.70 河渠 0.21 14.30 6.23 11.93 0.38 0.21 14.95 8.48 5.09 0 2.58 0.01 32.30 15.77 0.16 1.21 4.56 3.23 3.59 0.02 125.21 湖泊 0.33 1.53 0.04 0.16 0 0 1.95 1.68 0.79 0.17 0 17.03 0.28 1.48 0 0 0.04 0.90 0.59 0 26.97 水库/坑塘 0 1.58 0.57 0.62 0.10 0 0.83 0.80 0.66 0.16 0.71 0 0.25 0.75 0 0 0.36 0.04 0.11 0 7.54 滩地 0.03 3.62 1.50 3.75 0.18 0.01 3.42 3.27 2.73 16.86 1.24 0.17 0 7.18 0.02 0.31 0.78 1.31 0.79 0.01 47.18 城镇用地 0 0.44 0.07 0.07 0 0 0.12 0 0 0 0 0 0 0 0 3.05 0.46 0 0 0 4.21 农村居民点用地 0 7.04 2.07 4.57 0.03 0.13 5.45 4.59 3.11 0.62 0.05 0.05 1.57 1.71 0.50 0 2.47 0.50 0.06 0.02 34.54 工交建设用地 0.03 3.86 0.80 3.54 0.08 0.06 5.56 2.50 0.93 0.28 0.06 0 0.89 0.37 0 0.45 0 0.26 0.01 0 19.68 沙地 0 6.17 7.43 73.82 13.57 0.07 71.59 77.18 56.62 2.57 0.68 0.34 31.35 4.81 0 1.53 2.52 0 2.48 0.84 353.57 盐碱地 0 0.73 0 0 0 0 0.04 0.01 0.22 0 0 0 0 0 0 0.43 0 0 0 0 1.43 沼泽地 0.20 6.56 3.02 9.96 0.09 0.06 6.52 4.79 4.32 3.27 2.76 1.06 0.80 0 0 0.18 2.44 0.72 0.34 0.08 47.17 裸地 0 0.74 0.08 0.26 0 0 0.22 0.56 0.83 0.14 0.04 0 0.13 0 0.03 0.24 0.06 0.01 0.08 0 3.42 表 8 2005—2015年巴林右旗不同等级生境质量区域面积变化
Table 8 Area change of habitat quality index in 2005-2015 in Bairin Right Banner
生境质量等级 生境质量指数值 2005—2010年 2010—2015年 2005—2015年 面积变化/km2 变化率/% 面积变化/km2 变化率/% 面积变化/km2 变化率/% 低 0~0.15 232.18 30.66 355.37 35.92 587.55 77.60 较低 >0.15~0.35 -1 519.55 -41.39 141.05 6.55 -1 378.50 -37.55 中等 >0.35~0.60 -512.83 -18.81 133.01 6.01 -379.82 -13.93 较高 >0.60~0.75 416.27 80.96 -153.91 -16.54 262.36 51.03 高 >0.75 1 383.94 60.40 -475.54 -12.94 908.40 39.65 -
[1] GIBBES C, HAVLICK D G, ROBB J R.Land Use and Land Cover in a Transitioning Militarized Landscape[J].Journal of Land Use Science, 2017, 12(2/3):182-196. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/1747423X.2017.1313325
[2] DE ROSA M, VESTERGAARD O M, STAUNSTRUP J K, et al.Identifying Land Use and Land-use Changes (LULUC):A Global LULUC Matrix[J].Environmental Science & Technology, 2017, 51(14):7954-7962. http://cn.bing.com/academic/profile?id=f3604576b5a0f3593e054f2ff11a57e1&encoded=0&v=paper_preview&mkt=zh-cn
[3] COX A B, TAYLOR N T, REBEIN M A, et al.Land Use Changes From Unconventional Gas Development in Public Lands of the Fayetteville Shale[J].Natural Areas Journal, 2017, 37(2):233-239. DOI: 10.3375/043.037.0212
[4] MALKAMÄKI A, TOPPINEN A, KANNINEN M.Impacts of Land Use and Land Use Changes on the Resilience of Beekeeping in Uruguay[J].Forest Policy and Economics, 2016, 70:113-123. DOI: 10.1016/j.forpol.2016.06.002
[5] GUILLAUME T, MARANGUIT D, MURTILAKSONO K, et al.Sensitivity and Resistance of Soil Fertility Indicators to Land-use Changes:New Concept and Examples From Conversion of Indonesian Rainforest to Plantations[J].Ecological Indicators, 2016, 67:49-57. DOI: 10.1016/j.ecolind.2016.02.039
[6] BOLLIGER J, BTTIG M, GALLATI J, et al.Landscape Multifunctionality:A Powerful Concept to Identify Effects of Environmental Change[J].Regional Environmental Change, 2011, 11(1):203-206. DOI: 10.1007/s10113-010-0185-6
[7] LOPES L F G, DOS SANTOS-BENTO J M R, AREDE-CORREIA-CRISTOVO A F, et al.Exploring the Effect of Land Use on Ecosystem Services:The Distributive Issues[J].Land Use Policy, 2015, 45:141-149. DOI: 10.1016/j.landusepol.2014.12.008
[8] FLESKENS L, DUARTE F, EICHER I.A Conceptual Framework for the Assessment of Multiple Functions of Agro-ecosystems:A Case Study of Trás-os-Montes Olive Groves[J].Journal of Rural Studies, 2009, 25(1):141-155. DOI: 10.1016/j.jrurstud.2008.08.003
[9] 吴健生, 曹祺文, 石淑芹, 等.基于土地利用变化的京津冀生境质量时空演变[J].应用生态学报, 2015, 26(11):3457-3466. http://d.old.wanfangdata.com.cn/Periodical/yystxb201511027 WU Jian-sheng, CAO Qi-wen, SHI Shu-qin, et al.Spatio-temporal Variability of Habitat Quality in Beijing-Tianjin-Hebei Area Based on Land Use Change[J].Chinese Journal of Applied Ecology, 2015, 26(11):3457-3466. http://d.old.wanfangdata.com.cn/Periodical/yystxb201511027
[10] BHAGABATI N K, RICKETTS T, SULISTYAWAN T B S, et al.Ecosystem Services Reinforce Sumatran Tiger Conservation in Land Use Plans[J].Biological Conservation, 2014, 169:147-156. DOI: 10.1016/j.biocon.2013.11.010
[11] 贾有余, 李辉, 王钊齐, 等.苏州市吴中区生态保护红线区生态环境质量定量评价[J].生态与农村环境学报, 2017, 33(6):525-532. DOI: 10.11934/j.issn.1673-4831.2017.06.006 JIA You-yu, LI Hui, WANG Zhao-qi, et al.Quantitative Evaluation of Eco-environment Quality of the Terrestrial Ecological Red Line Area in Wuzhong District, Suzhou City[J].Journal of Ecology and Rural Environment, 2017, 33(6):525-532. DOI: 10.11934/j.issn.1673-4831.2017.06.006
[12] 包玉斌, 刘康, 李婷, 等.基于InVEST模型的土地利用变化对生境的影响:以陕西省黄河湿地自然保护区为例[J].干旱区研究, 2015, 32(3):622-629. http://d.old.wanfangdata.com.cn/Periodical/ghqyj201503029 BAO Yu-bin, LIU Kang, LI Ting, et al.Effects of Land Use Change on Habitat Based on InVEST Model:Taking Yellow River Wetland Nature Reserve in Shaanxi Province as an Example[J].Arid Zone Research, 2015, 32(3):622-629. http://d.old.wanfangdata.com.cn/Periodical/ghqyj201503029
[13] 刘洪来, 王艺萌, 窦潇, 等.农牧交错带研究进展[J].生态学报, 2009, 29(8):4420-4425. DOI: 10.3321/j.issn:1000-0933.2009.08.048 LIU Hong-lai, WANG Yi-meng, DOU Xiao, et al.Progress and Perspective of Agro-pasturage Ecotone[J].Acta Ecologica Sinica, 2009, 29(8):4420-4425. DOI: 10.3321/j.issn:1000-0933.2009.08.048
[14] 刘洪来, 张卫华, 王堃, 等.华北农牧交错带农田-草地界面土壤水分影响域分析[J].应用生态学报, 2009, 20(3):659-664. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200903028 LIU Hong-lai, ZHANG Wei-hua, WANG Kun, et al.Edge Influence of Soil Moisture at Farmland-grassland Boundary in Agriculture-pasturage Ecotone of Northern China[J].Chinese Journal of Applied Ecology, 2009, 20(3):659-664. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200903028
[15] WANG X H, LU C H, FANG J F, et al.Implications for Development of Grain-for-green Policy Based on Cropland Suitability Evaluation in Desertification-affected North China[J].Land Use Policy, 2007, 24(2):417-424. DOI: 10.1016/j.landusepol.2006.05.005
[16] 白光峰, 白光瑞, 李晓红, 等.东乌珠穆沁旗实施京津风沙源工程分析评价[J].内蒙古林业调查设计, 2003, 26(增刊1):67-68. http://d.old.wanfangdata.com.cn/Periodical/nmglydcsj2003z1033 [17] 李美君, 杜庆, 张克斌, 等.北方农牧交错带草地植被数量波动特征:以宁夏盐池县为例[J].东北林业大学学报, 2016, 44(1):48-51. DOI: 10.3969/j.issn.1000-5382.2016.01.009 LI Mei-jun, DU Qing, ZHANG Ke-bin, et al.Grassland Vegetation Quantitative Fluctuation in the Northern Farming Pastoral Zone:A Case Study of Yanchi County[J].Journal of Northeast Forestry University, 2016, 44(1):48-51. DOI: 10.3969/j.issn.1000-5382.2016.01.009
[18] 罗海江, 白海玲, 方修琦, 等.农牧交错带近十五年生态环境变化评价:以鄂尔多斯地区为例[J].干旱区地理, 2007, 30(4):474-481. http://d.old.wanfangdata.com.cn/Periodical/ghqdl200704002 LUO Hai-jiang, BAI Hai-ling, FANG Xiu-qi, et al.Ecosystem Change and Assessment of Ecotone Between Agriculture and Animal Husbandry Regions in North China in the Last Fifteen Years:A Case Study of Ordos Prefecture, Inner Mongolia[J].Arid Land Geography, 2007, 30(4):474-481. http://d.old.wanfangdata.com.cn/Periodical/ghqdl200704002
[19] 何立环, 董贵华, 王伟民, 等.中国北方农牧交错带2000-2010年生态环境状况分析[J].中国环境监测, 2014, 30(5):63-68. DOI: 10.3969/j.issn.1002-6002.2014.05.014 HE Li-huan, DONG Gui-hua, WANG Wei-min, et al.Ecosystem Status and Change Assessment of Agro-pastoral Ecotone of North China in 2000-2010[J].Environmental Monitoring in China, 2014, 30(5):63-68. DOI: 10.3969/j.issn.1002-6002.2014.05.014
[20] 岳书平, 张树文.内蒙古东部地区草地退化的时空分布特征:以巴林右旗为例[J].资源科学, 2007, 29(6):154-161. DOI: 10.3321/j.issn:1007-7588.2007.06.022 YUE Shu-ping, ZHANG Shu-wen.Temporal and Spatial Distribution of Grassland Degradation in Eastern Inner Mongolia:A Case Study in Balinyouqi[J].Resources Science, 2007, 29(6):154-161. DOI: 10.3321/j.issn:1007-7588.2007.06.022
[21] 贾晓燕.不同土地利用类型水源涵养服务特征及其时空动态[D].呼和浩特: 内蒙古大学, 2013. JIA Xiao-yan.Water Conservation Services of Different Land Uses[D].Hohhot: Inner Mogolia University, 2013.
[22] 闫国振, 张征, 梁康, 等.鄂尔多斯高原泊江海子流域土地利用动态变化及驱动因素分析[J].中国生态农业学报, 2017, 25(11):1693-1706. http://d.old.wanfangdata.com.cn/Periodical/stnyyj201711013 YAN Guo-zhen, ZHANG Zheng, LIANG Kang, et al.Characteristics and Driving Factors of Land Use Change in the Bojiang Lake Basin in Ordos Plateau, China[J].Chinese Journal of Eco-agriculture, 2017, 25(11):1693-1706. http://d.old.wanfangdata.com.cn/Periodical/stnyyj201711013
[23] SHARP R, KRAMER R C, WOOD S, et al.InVEST 3.3.0 User's Guide[M/OL].(2015)[2019-04-09].http://data.naturalcapitalproject.org/nightly-build/invest-users-guide/html/.
[24] 钟莉娜, 王军.基于InVEST模型评估土地整治对生境质量的影响[J].农业工程学报, 2017, 33(1):250-255. http://d.old.wanfangdata.com.cn/Periodical/nygcxb201701034 ZHONG Li-na, WANG Jun.Evaluation on Effect of Land Consolidation on Habitat Quality Based on InVEST Model[J].Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(1):250-255. http://d.old.wanfangdata.com.cn/Periodical/nygcxb201701034
[25] 王丽妍.生态空间质量约束下的城镇用地增长情景模拟[D].南京: 南京大学, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10284-1016122077.htm WANG Li-yan.Scenario Simulation for Urban Land Growth Under the Constraint of Ecological Space Quality[D].Nanjing: Nanjing University, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10284-1016122077.htm
[26] 贾树海, 李金典, 谢东祺, 等.内蒙古巴林右旗耕地及耕地后备资源生态安全研究[J].土壤通报, 2016, 47(4):777-782. http://d.old.wanfangdata.com.cn/Periodical/trtb201604002 JIA Shu-hai, LI Jin-dian, XIE Dong-qi, et al.Study on the Ecology Security Research of Cultivated Land and Its Reserve Resources in Balinyouqi of Inner Mongolia Autonomous Region[J].Chinese Journal of Soil Science, 2016, 47(4):777-782. http://d.old.wanfangdata.com.cn/Periodical/trtb201604002
[27] 白智慧.巴林右旗土地利用规划与生态评价研究[D].呼和浩特: 内蒙古师范大学, 2011. http://cdmd.cnki.com.cn/article/cdmd-10135-1011149569.htm BAI Zhi-hui.Bairin Youqi Land Use Planning of Ecological Evaluation Research[D].Hohhot: Inner Mogolia Normal University, 2011. http://cdmd.cnki.com.cn/article/cdmd-10135-1011149569.htm
[28] 刘军会, 高吉喜, 马苏, 等.内蒙古生态环境敏感性综合评价[J].中国环境科学, 2015, 35(2):591-598. http://d.old.wanfangdata.com.cn/Periodical/zghjkx201502035 LIU Jun-hui, GAO Ji-xi, MA Su, et al.Comprehensive Evaluation of Eco-environmental Sensitivity in Inner Mongolia, China[J].China Environmental Science, 2015, 35(2):591-598. http://d.old.wanfangdata.com.cn/Periodical/zghjkx201502035
[29] 张晶, 陈新明, 李琴琴, 等.宜城市生态环境敏感性评价及空间分布[J].中国土地科学, 2011, 25(6):76-83. DOI: 10.3969/j.issn.1001-8158.2011.06.014 ZHANG Jing, CHEN Xin-ming, LI Qin-qin, et al.Eco-environmental Sensitivity Assessment and Its Spatial Distribution in Yicheng City[J].China Land Science, 2011, 25(6):76-83. DOI: 10.3969/j.issn.1001-8158.2011.06.014
[30] 金宇.内蒙古科尔沁沙地生态系统受威胁状况评估[D].南京: 南京信息工程大学, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10300-1015578919.htm JIN Yu.Ecosystems Threatened Status Evaluation of Horqin Sandy Land in Inner Mongolia[D].Nanjing: Nanjing University of Information Science & Technology, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10300-1015578919.htm
[31] 哈斯巴根, 宝音, 李百岁.呼和浩特市土地资源人口承载力的系统研究[J].干旱区资源与环境, 2008, 22(3):26-32. DOI: 10.3969/j.issn.1003-7578.2008.03.005 Hasibagen, BAO Yin, LI Bai-sui.The Systematical Study on the Land Carrying Capacity of Huhhot[J].Journal of Arid Land Researches and Environment, 2008, 22(3):26-32. DOI: 10.3969/j.issn.1003-7578.2008.03.005
[32] 胡云锋, 刘纪远, 齐永青, 等.内蒙古农牧交错带生态工程成效实证调查和分析[J].地理研究, 2010, 29(8):1452-1460. http://d.old.wanfangdata.com.cn/Periodical/dlyj201008010 HU Yun-feng, LIU Ji-yuan, QI Yong-qing, et al.Positivist Analysis on the Effects of Ecological Projects in the Farming-pastoral Transition Belt of Inner Mongolia Autonomous Region[J].Geographical Research, 2010, 29(8):1452-1460. http://d.old.wanfangdata.com.cn/Periodical/dlyj201008010
[33] 国润才, 傅恒, 刘利胜.内蒙古地区退耕还林(草)工程效益浅析[J].内蒙古林业科技, 2004, 30(2):33-34, 52. DOI: 10.3969/j.issn.1007-4066.2004.02.009 GUO Run-cai, FU Heng, LIU Li-sheng.Analysis on Benefits of Program for Returning Cropland to Forestry in Inner Mongolia[J].Inner Mongolia Forestry Science & Technology, 2004, 30(2):33-34, 52. DOI: 10.3969/j.issn.1007-4066.2004.02.009