赤峰市位于内蒙古自治区东南部，地理坐标北纬(41˚7′10ʺ~45˚24′15ʺ)，东经(16˚07′43ʺ~120˚58′52ʺ)。西、北与锡林郭勒盟临近，东、东南与通辽市和朝阳市相连，西南连通承德市(见图1)。辖区东西宽375 km，南北长457.5 km，总面积90275 km²。下设三区七旗两县，赤峰市行政区内所含21个街道办事处，132个苏木、乡镇 [ 8 ] 。

图1. 赤峰市行政图

通过阅读文献，搜集相关资料，本文所利用的土地利用数据来自中国科学院资源与环境科学数据中心，分辨率为30米，根据研究区实际和研究需要将土地利用类型归并为耕地、林地、草地、水域、建设用地、未利用地6个一级地类。

土地利用转移矩阵是一个指标，表明土地利用的结构特征和数量方向的变化，可以充分反映土地在特定时期的变化方向 [ 9 ] [ 10 ] ，其表达式为：

S i j = [ S 11 S 12 ⋯ S 1 n S 21 S 22 ⋯ S 2 n ⋮ ⋮ ⋱ ⋮ S n 1 S n 2 ⋯ S n n ]

式中，n：土地利用类型数，i：研究时段内初期土地利用类型，j：研究时段内末期土地利用类型，S_ij：土地利用类型从i转变为j的转移面积。利用1990~2020年四期土地利用数据，运用ArcGIS软件空间叠加功能，获得研究区土地利用转移数据。

土地利用重心迁移模型可从空间上阐述各土地类型的变化轨迹，直观地反映土地利用变化的方向和速度 [ 11 ] ，重心坐标计算公式为：

X i = ∑ j = 1 n ( C i j × X i j ) ∑ j = 1 n C i j

Y i = ∑ j = 1 n ( C i j × Y i j ) ∑ j = 1 n C i j

式中， X i 、 Y i 分别为第i种地类重心经纬度， C i j 为第i种地类中的第j个斑块面积， X i j 、 Y i j 分别为第i种地类中的第j个斑块的经纬度。利用4期土地利用数据和重心迁移模型 [ 12 ] ，可直观地反映研究区40年间土地利用变化方向和趋势。

土地利用程度变化模型可以根据城市土地利用程度进行量化分析 [ 13 ] [ 14 ] 。计算公式为：

Δ L b - a = L b − L a = [ ∑ i = 1 n A i C i b − ∑ i = 1 n A i C i a ] × 100

式中： Δ L b - a 为土地利用程度变化指数； L a 和 L b 时是分别a点和b时点的土地利用程度指数； A i 为第i级土地利用程度分级指数； C i a 和 C i b 分别为a时点和b时点的第i级的土地利用程度面积比。若 Δ L b - a > 0 ，则表明该区域土地利用处于发展时期，反之，则处于衰退或调整时期。

参照景观生态学观点和《地理格网国家标准GB/T 12409》相关理论要求，借鉴已有学者的研究，根据研究区实际情况和各土地利用类型分布情况，采用等间距系统抽样法将赤峰市土地利用数据划分为8 km × 8 km空间网格 [ 15 ] [ 16 ] ，共得到1519个风险单元。

本文参考相关研究并选取以下景观指数，构建生态风险指数模型计算生态风险指数 [ 15 ] [ 17 ] (见表1)。

通过计算景观干扰度指数、景观脆弱度指数和景观损失度指数来建立景观生态风险指数，公式为：

E R I k = ∑ i = 1 n A k i Z i A k

E R I k 为第k个风险单元的景观生态风险指数；n为景观类型数量； A k i 为第k个风险单元景观类型i的面积； A k 为第k个风险单元的总面积； Z i 为景观损失度指数 [ 18 ] 。

表1. 景观格局指数构建模型表

从表2可以看出1990~2000年间，赤峰市土地利用类型发生了明显的变化。一方面草地、水域、未利用地总体减少，分别减少了1974.44 km²、89.37 km²、28.21 km²、430.66 km²；耕地、林地和建设用地面积增加，增加面积分别为2444.38 km²、29.40 km²、20.91 km²。2000~2010年间耕地、林地、建设用地面积持续增加，分别为1134.79 km²、7137.53 km²、225.66 km²。相比于1990~2000年林地面积扩大最多，建设用地相对增加了较多面积，在这一时期，社会经济的快速发展和人口的快速增长对建筑用地产生了强烈的需求。除此水域和未利用地也增加了，增加面积分别为63.35 km²和585.86 km²。而草地面积持续较少状态，减少面积为9146.83 km²。2010~2020年间除了草地面积其余各类土地面积全部减少，耕地、林地、水域、建设用地、未利用地依次减少为376.52 km²、7120.74 km²、91.04 km²、592.81 km²。而草地面积增加了8269.97 km²。

总体来看，赤峰主要以草地、耕地和林地为主，1990~2020年间草地面积均占总面积百分之40%以上，耕地面积每时段占总面积的18%以上，林地均占14%以上，三类土地占总面积的三分之二以上土地，研究区建设用地占比最少，其次为水域用地。1990~2010年间耕地面积持续增加状态，2010~2020年间变化较少，而建设用地无太大扩张反而较稳定发展，说明近年来，主要是通过严格的土地管理规划和审批，积极规范新建建设用地和结构控制，实施耕地保护、环境保护等举措得以改善(见表2)。

表2. 1990~2020年研究区各土地利用面积与占比表

表3. 研究区土地利用转移变化

根据表3可知，1990~2000年间，草地转出面积最大，转为耕地面积为2284.09 km²，其次为未利用地、耕地和林地，未利用地转出和转入面积分别为492.99 km²和62.33 km²，耕地转出面积为182.08 km²，而转入面积2626.43 km²，林地转出面积为157.01 km²，转入面积为186.26 km²，草地面积大幅度减少，转出面积和转入面积分别为2440.70 km²和466.22 km²，水域用地和建设用地波动较少，转出面积和转入面积分别为115.22 km²和25.84 km²，建设用地变动较为稳定。2000~2010年间，各类土地转移较大，耕地转出和转入面积分别为4966.67 km²和6101.37 km²，林地转出和转入面积分别2464.48 km²和9610.92 km²，草地面积依然减少，转出面积为16023.41 km²，其中转变为林地面积最大，草地转入面积为6876.48 km²，水域用地转出和转入面积较均匀，建设用地转出面积和转入面积分别为626.15 km²和851.78 km²，未利用地转入面积比转出面积大585.83 km²。2010~2020年间，林地转出面积最大，转为草地面积的为8018.68 km²，林地转入面积为2354.30 km²，其中草地转入面积最大，为15375.77 km²，耕地、水域和建设用地转出和转入面积相差不大，未利用地转出面积比转入面积大592.76 km²(见表3)。

从图2可以看出赤峰市四期的土地利用类型重心迁移轨迹，很好地反映了赤峰市6大地类的时空演变过程。1990~2020年赤峰市各类土地利用重心迁移具有差异性，根据图可知，赤峰市耕地，林地，草地，水域和建设用地土地利用重心均在翁牛特旗，而未利用地土地利用重心在巴林右旗(见图2)。

图2. 1990~2020年赤峰市土地利用类型重心迁移图

耕地重心迁移特征

1990~2020年耕地重心向南方向迁移了2.07 km，1990~2000年向西偏南方向迁移了5.86 km；2000~2010年耕地重心向东偏南方向迁移了15.87 km，该时期是重心迁移距离最大且最快的时期，2010~2020年耕地重心向西偏北方向迁移了9.76 km。

林地重心迁移特征

1990~2020年林地重心向南偏东方向迁移了1.58 km。其中1990~2000年迁移方向为西偏南；2000~2010年林地重心迁移方向为南偏东，迁移距离为37.19 km，该时期是林地重心迁移距离最大且迁移最快的时期；2010~2020年林地重心迁移方向为北偏西，迁移距离为34.18 km。

草地重心迁移特征

1990~2020年草地重心从西南向东北方向迁移，迁移距离为11.45 km。其中1990~2000年向西南方向迁移了0.67 km，该时期是重心迁移距离最小、重心迁移最慢的时期；2000~2010年草地重心向西偏南迁移了25.77 km；2010~2020年草地重心迁移方向为南偏西，迁移距离为0.32 km。

水域重心迁移特征

1990~2020年水域重心向北偏东迁移了6.14 km。其中2010~2020年是水域重心迁移距离最大且最快的时期，迁移距离为25.55 km，迁移方向为北偏东；1990~2000年是水域重心迁移距离最小、重心迁移最慢的时期，迁移距离为2.05 km，迁移方向为东南；2000~2010年水域重心向南偏西迁移了19.53 km。

建设用地重心迁移特征

1900~2020年建设用地重心向北迁移了1.09 km。其中1990~2000年向西迁移了0.09 km，该时期是重心迁移最慢的时期；2000~2010年建设用地重心迁移方向为南偏西，迁移距离16.82 km；2010~2020年建设用地重心向北偏东迁移了17.82 km，该时期是建设用地重心迁移最快的时期。

未利用地重心迁移特征

1990~2020年未利用地重心向东北迁移了3.11 km。其中1990~2000年向东北方向迁移了3.15 km，2010年未利用地重心与2020年的重心几乎重合，2000~2010年未利用地重心向西偏南迁移了23.42 km，该时期是未利用地重心迁移最快的时期；2010~2020年未利用地重心迁移方向为向东偏北。

表4. 赤峰市及其各辖区土地利用程度变化量和变化率

根据表4来看，1990~2020年期间赤峰市各辖旗县区土地利用程度变化量和变化率均大于0，这表明，赤峰市的土地利用率总体上升，普遍提高。其中阿鲁科尔沁旗、敖汉旗、巴林左旗、翁牛特旗、元宝山区和红山区土地利用变化量大于平均变化量，这表明，这些地区的经济增长速度更快，土地利用正在发生巨大变化。而克什克腾旗土地利用变化量小于1，处于衰退期，其余区域变化量均大于1，但整体小于平均变化量，说明处于发展期，发展状况低于上述六个区域。

分时期来看，1990~2000年赤峰市全部区域土地利用程度变化量均是正数，最为突出的区域集中在非中心城区敖汉旗，土地利用程度变化量达到8.95，变化率为3.69%，其次是阿鲁科尔沁旗、巴林左旗、宁城县、翁牛特旗、元宝山区土地利用程度变化量均大于该时期平均程度；而巴林右旗、克什克腾旗、林西县和红山区的土地利用程度变化较不明显，变化率均不到1%。

2000~2010年期间，赤峰市敖汉旗和克什克腾旗土地利用程度变化量和变化率小于0，说明该阶段两个区域土地利用处于调整期或衰退期；而元宝山区土地利用程度的涨幅相对较大，处于快速发展期，变化量在该阶段中最大，变化率为2.94%。相比于1990~2000年红山区和林西县土地利用程度的涨幅也较大，此外喀喇沁旗、翁牛特旗、宁城县和松山区土地利用程度发展较慢，变化量均低于平均变化量。其余区县变化量均大于2.69，说明处于快速发展期。

2010~2020年期间，赤峰市整体处于缓慢发展，平均土地利用程度变化两位−1.38。其中克什克腾旗相比于前两段时期快速发展，是土地利用程度变化较为剧烈的地区，增幅明显，远高于平均变化量，变化率为1.84%。此外阿鲁科尔沁旗、敖汉旗、翁牛特旗变化量为正数，发展较快。其余区县土地利用程度变化量和变化率小于0，该阶段这些区域土地利用处于调整期或衰退期(见表4)。

(1) 景观格局指数分析

表5. 赤峰市四期各用地类型景观格局指数

从景观破碎度指数来看，各地类在研究时段内的破碎度均小于0.05，说明研究区整体景观破碎化程度较低，建设用地破碎化程度在整体景观中最高，但呈现不稳定趋势。1990~2000年破碎度指数变化最为明显的用地类型为耕地，表明在此时段内研究区耕地破碎化程度大，其他地类破碎度相对稳定。2000~2010年草地和建设用地破碎化程度增加，耕地、林地、水域和未利用地的破碎度则表现为降低的态势。2010~2020年研究区破碎度趋势与上一阶段相反，耕地、林地、水域和未利用地的破碎度增大(见表5)。

从景观分离度指数来看，各景观类型的景观分离度指数差异较大，建设用地受外界影响较大，面积小，且分散分布，分离度指数最大。1990~2000年草地、水域和未利用地受人类活动影响，分离度指数增幅较大，表明分布分散，而耕地、林地、建设用地的分离度指数下降，表明斑块间距缩小。2000~2010年草地和建设用地分离度上升，说明离散程度增大，其他地类则呈现下降趋势，2010~2020年与上一阶段相反，耕地、林地、水域和未利用地分离度指数明显增加，草地和建设用地分离度指数则减少(见表5)。

从景观优势度指数来看，各地类的优势度指数变化不稳定。1990~2000年草地优势度下降，其他地类优势度均有所增加。2000~2010年建设用地优势度指数增加，其他地类均下降。2010~2020年耕地、林地、草地和未利用地优势度上升，水域和建设用地指数则下降。景观优势度指数显示了它对景观结构的重要性，随着景观优势度指数的增加，景观稳定的影响也在增加(见表5)。

景观损失度是景观干扰度与景观脆弱度的直接表现，脆弱度指数从大到小为：未利用地 > 水域 >耕地 > 草地 > 林地 > 建设用地。研究时段内各地类的景观干扰度指数相差不大，建设用地在每个时期内干扰度最大，主要该地类面积和斑块数目均较少，分布较为分散，且受到外界干扰的程度减大。水域和未利用地破碎度和脆弱度都较高，受外界干扰最明显，损失度较大，耕地、林地、草地是研究区占比面积最大的三个地类，受人类活动和自然环境影响较大，损失度不稳定(见表5)。

(2) 景观生态风险值变化特征

基于构建的景观生态风险评价模型，计算研究区三个时期内1519个样区的生态风险指数可得，赤峰市四个时期的生态风险指数平均值分别为0.012234、0.012516 、0.012244和0.012267，1990~2000和2010~2020年上升，结果表明，由于土地使用利用结构的变化，研究地区的生态安全恶化。

根据表6分析各土地利用类型的景观生态风险指数可知，1990~2020年各地类的景观生态风险平均值均小于0.025。其中，虽然水域和未利用地景观的损失度较高，但由于面积相对较小，不会对生态环境安全构成威胁。1990~2000年耕地、林地、草地和水域的景观生态风险指数呈增加的趋势，而建设用地和未利用地生态风险指数下降，这表明生态安全正在改善。2000~2010年耕地、林地和未利用地景观生态风险指数下降，林地指数增加，水域则稳定不变。2010~2020年耕地、草地和水域景观生态风险指数增加，林地和建设用地指数下降，而未利用不变(见表6)。

表6. 赤峰市四个时期各用地类型生态风险指数

为对比四个时期景观生态风险的时空演变，通过ArcGIS软件中的克里金插值法得到赤峰市四个时期的生态风险空间分布状况，并通过自然间断点分级法将其划分为五个等级：低生态风险区(ERI ≤ 0.0081)、较低生态风险区(0.0081 < ERI ≤ 0.0109)中生态风险区(0.0109 < ERI ≤ 0.0137)、较高生态风险区(0.0137 < ERI ≤ 0.0165)和高生态风险区(ERI > 0.0165) (见图3)。

图3. 赤峰市1900、2000、2010、2020年土地利用的生态风险空间分布

由图3可知，研究区生态风险等级主要以中风险区和较高风险区为主。从各景观生态风险在各旗县的分布看，整体上赤峰市1990~2020年的景观生态风险空间分布变化并无差异。其中，中风险区分布在研究每一个旗县，整体风险由中心向四周梯度递减，呈“较高–中–较低”的空间分布格局，即研究区整体风险表现为西部及中部地区景观生态风险值较大。低和较低生态风险等级区主要分布在南部区域及边缘区域，与中生态风险等级、较高生态风险等级紧邻并占区域全部风险等级的少部分。在研究时期研究区耕地、草地、建设用地、未利用地分散分布，而耕地和建设用地的破碎度和分离度较大，明显高于区域内其他用地类型，因此该区域整体处于中和较高风险区。其中生态风险区较高风险区分布在阿鲁科尔沁旗中南部、巴林左旗南部、巴林右旗东南部、敖汉旗中部、翁牛特旗及克什克腾旗，高风险区主要分布在克什克腾旗、松山区、翁牛特旗中部和阿鲁科尔沁旗西南部，而较低和低风险区分布在研究区南部、阿鲁科尔沁旗北部，分布比较分散，中风险区几乎分布在研究区每一个区域，分布范围最广，所占面积最大。

(3) 景观生态风险等级变化分析

为进一步解释赤峰市不同等级生态风险区的变化，对1990、2000、2010和2020年生态风险空间分布图进行空间分析，得到赤峰市1990~2020年各生态风险等级的面积、比例及转移面积矩阵(见表7、表8)。

表7. 赤峰市四个时期各生态风险等级面积及比例

据表7可知在四个时期研究区中风险区占比最高，占总面积45%左右，其次为较高风险区，占比为33%左右，低风险和高风险区占比最低。1990~2000年研究区高风险区面积增加，增幅为0.41%，增加面积为353.05 km²，其余等级风险区面积均减少。2000~2010年低风险、较低风险和中风险区面积持续减少，减少面积分别为77.28 km²、143.24 km²和633.17 km²，降幅分别为0.09%、0.16%和0.73%，中风险和高风险区分别增加829.28 km²和24.41 km²，占比分别从原来的33.12%、2.29%增到34.07%和2.32%。2010~2020年低风险、较低风险区面积增加，增幅分别为0.35%和0.29%，除此高风险区面积持续增加，增加面积为23.63 km²，而中风险和较高风险区面积减少，减少面积分别为422.32 km²和159.87 km²。

表8. 赤峰市1990~2020年生态风险等级转移面积矩阵(km²)

据表8可知赤峰市生态风险等级保持不变的区域最大，面积为80277.84 km²，约占总面积的92.22%；生态风险等级降低的区域面积4267.29 km²，生态风险等级升高的区域总面积为2504.93 km²。具体而言，较高生态风险区向中生态风险区转移面积最多，达到2262.67 km²；其次是较低生态风险区转向中生态风险区、中生态风险区向较高生态风险区以及高风险区向较高风险区转移，转移面积分别为959.86 km²和858.26 km²和824.91 km²。由此表明：研究期间内赤峰市低风险和较低风险区面积占比虽小，但景观生态风险等级呈无太大升高趋势。