表1  一级景观面积及斑块数量变化总览图
 2009 2016 变化分析
景观类型 面积/hm² 百分比/% 小班数/个 面积/hm² 百分比/% 小班数/个 面积变化/hm² 小班数变化/个
乔木林 9577.01 12.58 437 12446.94 16.36 1600 2869.77 1163
疏林地 3240.89 4.26 151 2299.41 3.02 227 -941.51 76
灌木林地 233.74 0.31 8 467.31 0.61 88 233.57 80
未成林地 1924.46 2.52 89 3972.15 5.22 666 2047.7 577
荒山沙地 30746.07 40.37 97 20110.29 26.43 224 -10635.78 127
牧草地 0 0 0 4378.03 5.75 14 4378.02 14
采伐迹地 0 0 0 818.36 1.08 322 818.36 322
耕地居民点 30359.08 39.86 377 31597.16 41.53 468 1229.74 91
总计 76156.0 100 1159 76081.26 100 3610 79.1 2450
从表1可知，景观总面积减少了79.1hm2，斑块个数却增加了2450个，这反映出了2016年区域景观破碎化程度相较2009年增加很多，空间异质性增大，景观格局的稳定性有所降低。在7年的时间里，荒山荒地面积大幅度减少。该区域景观格局仍然呈现出耕地景观为主导，乔木林地、疏林地、灌木林地、未成林造林地、荒山沙地、牧草地、采伐迹地等景观格局呈补丁状散布之中的景观格局现状。
景观转移的轨迹及特征：表2表明荒山沙地和疏林地为主要的面积输出景观，而且其他景观皆为面积转入，荒山沙地、耕地向其中转移最多为别为2358.41hm2、1736.46hm2； 295.54hm2的疏林地由于合理的经营管理变为乔木林，但也有483.34hm2的天然疏林地由于林地退化不得不规划成完全受人类活动控制的耕地；牧草地对研究区域的自然因素的适应性最好，4378.03hm2的还草工程牧草地景观全部成活并且涨势良好。

表2  2009-2016年景观类型转移矩阵（单位：hm2）
2009
2016 乔木林 疏林地 疏林地 未成林 采伐迹地 荒山沙地 牧草地 耕地 2009总计
乔木林 7352.88 1.74 59.23 266.15 710.66 1.03 0 1185.32 9577.01
疏林地 295.54 2134.31 7.82 213.76 106.11 0 0 483.34 3240.88
灌木林 6.53 0 169.18 0 0 0 0 58.05 233.76
采伐迹地 0 0 0 0 0 0 0 0 0
未成林 697.13 8.08 24.17 763.84 1.54 0 0 429.69 1924.45
荒山沙地 2358.41 56.59 130.75 2207.99 0 20109.26 4378.03 1505.04 30746.07
牧草地 0 0 0 0 0 0 0 0 0
耕地 1736.46 98.69 76.16 520.4 0 0 0 27927.38 30359.09
2016总计 12446.95 2299.41 467.31 3972.14 818.31 20110.29 4378.03 31588.82 76081.26
3.景观变化驱动力分析
研究区域土质松散且呈沙性，在天然植被的调节下，保持沙地生态系统相对平衡状态，自然环境不会产生剧烈的退化，然而人为干扰过度，沙性土壤潜在的自然因素便会激化与活化，产生沙地[5]，通过分析林政记录材料、实地踏查验证，总结出了以下8种主要的引起景观类型变化的原因，结果见表3
表3  2009到2016年间景观演变驱动力分析
变化原因 变化面积/hm² 百分比/% 变化小班数/个
人工造林或飞播造林 6914.94 42.09 1349
森林采伐 1342.06 8.16 438
经审批使用林地 6.47 0.05 3
退耕还草 4634.37 28.20 108
干旱 2024.21 12.32 412
沙化 1017.65 6.19 149
其他自然因素 266.31 1.62 189
其他调查因素 225.49 1.37 65
总计 16431.50 100 2713
2009-2016年期间景观变化面积总计16431.50hm2，景观斑块变化个数总计2713个。在研究区域的各种灾害类型之中由于干旱和沙化引发的景观类型变化占主导，分别为2024.21hm²和1017.65hm²占变化面积总量的12.32%和6.19%，与这两种相比包括虫灾、白灾、鼠灾、火灾在内的其他自然因素所引起的景观类型变化总计只有266.31hm²占1.62%。
四、景观格局预测
通过对2009年-2016年各个时段景观变化面积数据的叠加分析，研究区域景观类型转移概率表（见表4）。
表4  2009-2016年景观类型转移概率（%）
2009
2016 乔木林 疏林地 疏林地 未成林 采伐迹地 荒山沙地 牧草地 耕地
乔木林 76.78 0.02 0.62 2.78 7.42 0.01 0 12.38
疏林地 9.12 65.86 0.24 6.60 3.27 0 0 14.91
灌木林 2.79 0.00 72.37 0.00 0.00 0 0 24.83
未成林 36.22 0.42 1.26 39.69 0.08 0 0 22.33
采伐迹地 0 0 0 0 0 0 0 0
荒山沙地 7.67 0.18 0.43 7.18 0 65.40 14.24 4.90
牧草地 0 0 0 0 0 0 0 0
耕地 5.72 0.33 0.25 1.71 0 0 0 91.99
已知2009年转化到2016年的转移矩阵A和2016年的状态向量B，假设在2016-2023年这段时间内，各土地利用类型的转移概率保持不变，那么就可以通过A×B，算出2023年的各种土地利用类型的面积向量以及2016－2023年研究区域的土地利用转移矩阵，利用转移概率矩阵，使用马尔科夫模型可以预测未来相同时间间隔内的土地利用变化，并与得到的2023年景观类型面积用初始状态矩阵，与第7年的转移概率相乘，得到第7年末的状态矩阵。在该研究中可以算出第7年末各景观的预测面积（见表5）。
表5  2023年景观格局预测
景观类型 2016年景观面积（hm2） 2023年景观面积（hm2）
乔木林 12446.94 14567.34
疏林地 2299.41 1672.94
灌木林 467.31 635.39
未成林造林地 3972.15 4059.84
采伐迹地 818.31 1820.40
荒山沙地 20110.29 13154.34
牧草地 4378.03 7241.60
耕地 31597.16 32929.42
由表5可知，未来7年胡花硕地区主要土地类型的变化趋势，荒山沙地面积将进一步减少，这是沙地生态结构进一步优化的趋势，将有利于当地生态环境的改善，面积增加最快的是牧草地，由4378.03hm2增加7241.60hm2，这不但有利于区域沙地生态系统优化，还有利于景观空间结构与功能的优化；值得注意的是疏林地的面积继续降低，作为研究区域原有自然景观类型，天然榆树林是疏林地的主要成分，预测结果表明其面积将进一步下降，这将导致生物多样性降低、珍稀动植物资源丧失以及抗干扰（包括人为和自然干扰）能力的降低。倘若保持现有的转移方式不加以干预，研究区域未来景观结构并不令人满意，维持现有内、外条件的前提下，虽然研究区域植被覆盖率整体是上升趋势，但以树种单一的人工林、结构简单的牧草地为主体的景观结构令人担忧，如不加以干预，该区域景观格局将向一个非常危险的方向发展，这种景观格局发展趋势不利于区域生态系统的稳定。
五、总结
对研究区域景观变化进行驱动力分析发现2009-2016年期间景观面积变化16431.50hm2，景观斑块变化2713个，其中自然原因造成的景观变化为3308.17hm2，占20.13%，人为原因造成的景观变化为13123.33hm2，占79.87%。通过预测得出2023年沙地景观格局，如果保持现有的转移方式，未来景观结构并不令人满意，现有的景观格局发展趋势不利于沙地生态系统的恢复，如果按此转移概率继续模拟，景观格局将不能满足科学发展的总方针和生态恢复的目标。

参考文献：
[1]王韶泽、姜乐平.地物波谱特性及实例[J].内江科技，2017(1):136-136
[2]李书娟、曾辉等.景观空间动态模型研究现状和应重点解决的问题[J].应用生态学报，2004.15(4):701-706
[3]杨学军、姜志林.状态空间模型在区域森林景观生态经济研究中的应用[J].系统工程理论与实践，2001(12):111-116
[4]李德成、徐彬彬、石晓日.利用马氏过程模拟和预测土壤侵蚀的动态演变——以安徽省岳西县为例[J].环境遥感，1995.10(2):89-96
[5]王涛、赵哈林、肖洪浪.中国沙漠化研究的进展[J].中国沙漠，1999.19(4):299-311
作者简介：
王梁泓（1994—），毕业于东北林业大学林学系，硕士研究生学历。研究领域：林业信息工程。