土地利用和覆被变化(land use and cover change, LUCC)是全球环境变化的重要组成部分和驱动因素之一[1], 研究土地利用变化规律和发展趋势, 对于土地资源的合理开发与利用, 促进区域协调可持续发展具有重要意义[2-3]. 近年来, 国内外学者已对土地利用变化时空动态预测模型进行了深入研究, 采用的预测模型主要有CLUE-S模型、SLEUTH、多智能体系统模型(multi-agent system, MAS)、元胞自动机模型(cellular automata, CA)等[4-9]. 其中, CA-Markov模型耦合了元胞自动机和马尔可夫链的各自优势, 可以实现对长时间复杂土地利用类型时空演变过程的可视化模拟和预测, 并具有较高的精度. 国内不少学者应用CA-Markov模型在新疆克州、昌化江流域、云南德宏州、三峡库区、黄河流域、锡林河流域、喀什市、长沙市等地区开展土地利用变化分析与预测研究[10-16].
大凌河是辽宁西部地区最大的河流, 位于大凌河流域的锦州、阜新、朝阳三市是东接辽沈, 西联京津唐地区的重要城市. 大凌河流域内地形复杂, 石漠化、水土流失等问题严重, 其生态系统稳定直接关系着辽西的经济发展[17-18]. 一些学者在大凌河流域开展的土地利用变化研究有力揭示了流域生态环境变化特征[19-20], 只是研究数据较陈旧, 不能反映流域新变化. 鉴于此, 本研究基于自然资源部新发布的土地覆盖数据, 从土地利用动态度、土地利用程度和土地利用转移矩阵等不同角度, 分析研究区过去20年来来土地利用变化规律, 并基于CA-Markov模型对研究区2030年土地利用格局进行预测, 以期为大凌河流域生态环境保护与区域经济和社会发展提供科学依据.
1 研究区简介研究区大凌河流域覆盖的地理范围为: 118°53'-121°52'E、40°28'-42°38'N. 河流全长435 km, 上游分南、北两支, 于喀左县大城子东南汇合后流经朝阳、北票、凌海、义县等地, 最终汇入渤海(图 1). 流域总面积2.33×104 km2, 其中2.08×104 km2位于辽宁省境内. 地貌以山地丘陵为主, 少量平原区. 气候类型属于中温带气候, 四季冷暖干湿分明, 温度变化较大, 多年平均气温8.3℃, 平均相对湿度53%, 日照时数2 800 h, 年均降雨量465 mm, 年蒸发量1 974.4 mm, 年均径流量17.91×108 m3 [19].
本研究主要基于自然资源部发布的30 m全球地表覆盖数据(GlobeLand30). 该数据从2014年开始发布, 目前已完成2020年版更新. GlobeLand30数据研制所使用的分类影像主要有美国陆地资源卫星(Landsat)、中国环境减灾卫星(HJ-1)和高分一号(GF-1)等多光谱影像, 优选植被生长季少云影像[21]. GlobeLand30数据共包括10个一级类型, 分别是: 耕地、林地、草地、灌木地、湿地、水体、苔原、人造地表(建设用地)、裸地、冰川和永久积雪, 研究区内共有其中7类数据(如表 1所示). GlobeLand30数据精度评价由同济大学和中国科学院空天信息创新研究院等牵头完成, 平均总体精度为84.61%, 平均Kappa系数为0.80.
本次研究共收集了6幅数据, 数据列表为: N50_40_2000LC030、N51_40_2000LC030、N50_40_2010LC030、N51_40_2010LC030、N50_40_2020LC030、N51_40_2020LC030. 数据处理与分析采用ArcGIS和IDRISI软件完成.
3 研究方法 3.1 土地利用动态度土地利用动态主要用于表征一段时间内土地利用类型的变化情况, 即土地利用的变化速度, 分为单一土地利用动态度和综合土地利用动态度[22-24]. 计算公式如下:
$ L=\frac{U_{\mathrm{b}}-U_{\mathrm{a}}}{U_{\mathrm{a}}} \times \frac{1}{T} \times 100 \% $ | (1) |
$ L C=\frac{\sum\limits_{i=1}^{n} \Delta L U_{i-j}}{2 \sum\limits_{i=1}^{n} L U_{i}} \times \frac{1}{T} \times 100\% $ | (2) |
式中, L为研究时段内单一土地利用类型的动态度; Ua和Ub分别表示该类土地利用类型研究始、末的面积; LC为综合土地利用动态度, LUi为第i类土地类型在研究初期的面积; ΔLUi-j为研究时段内第i类土地利用类型转变为非i类(j类)土地利用类型面积的绝对值; T为研究时段, 当该时段设定为年时, L为该研究区某种土地利用类型年变化率, LC为研究区土地利用年变化率[23].
3.2 土地利用程度土地利用程度是在某个时间内由自然因素和人为因素双重作用的结果, 它反映了土地利用的广度和深度. 根据刘纪远等[1]提出的土地利用程度的综合分析方法, 按照土地自然综合体在社会因素影响下的自然平衡状态将土地利用程度分为4个级别, 依次赋予分级指数(表 2), 在此基础上建立土地利用程度综合指数及变化模型[25].
1) 土地利用程度综合指数
土地利用程度综合指数计算公式如下:
$ L_{i}=\sum\limits_{i=1}^{n} A_{i} C_{i} \quad L \in[100,400] $ | (3) |
式中, Li表示土地利用程度综合指数; Ai表示第i级的土地利用程度分级指数; Ci表示第i级的土地利用类型面积百分比; n为土地利用程度分级数. Li数值位于100~400区间内, 数值大小反映土地利用程度的高低.
2) 土地利用程度变化
土地利用程度的变化是指特定区域在一定的时段内的综合变化, 通常由土地利用程度的变化量表征, 计算公式如下:
$ \Delta L_{b-a}=L_{b^{-}} L_{a}=100 \times\left(\sum\limits_{i=1}^{n} A_{i} C_{i b^{-}} \sum\limits_{i=1}^{n} A_{i} C_{i a}\right) $ | (4) |
式中, La和Lb分别为a、b两个时间的研究区土地利用程度综合指数; Ai为第i级的土地利用程度分级指数; Cia和Cib分别为研究区a时间和b时间第i级土地类型利用程度面积百分比. 如果ΔLb-a>0, 表示研究区土地利用处于发展时期; 如果ΔLb-a < 0, 表示该区域的土地利用处于衰退期.
3.3 土地利用转移矩阵土地利用变化不仅是面积数量上的增减, 而且还存在类型之间的转换关系. 土地利用转移矩阵能够定量描述各种土地利用类型之间的变化, 有助于了解期初各地类的流失去向和期末各地类的来源和构成. 通过对任意两期土地利用数据进行空间叠加, 获取各研究时段土地利用类型转移矩阵, 其表达式为:
$ S_{i j}=\left[\begin{array}{cccc} S_{11} & S_{12} & \cdots & S_{1 n} \\ S_{21} & S_{22} & \cdots & S_{2 n} \\ \vdots & \vdots & \vdots & \vdots \\ S_{n 1} & S_{n 2} & \cdots & S_{n n} \end{array}\right] $ | (5) |
式中, S为面积; n为土地利用类型数量; i、j分别为研究初期与末期的土地利用类型.
3.4 CA-Markov模型CA-Markov模型结合了元胞自动机模型和马尔可夫链的优势, 常用于土地利用时空格局的模拟预测, 具有高精度、可视化等优点[26]. Markov链描述的是土地利用从一个时刻(t)到另一个时刻(t+1)的变化趋势, 并以此为规则预测将来的变化趋势[26-27]. 公式如下:
$ S_{t+1}=P_{i j} \times S_{t} $ | (6) |
式中, St、St+1分别为t、t+1时刻的系统状态, Pij为状态转移概率.
CA模型同时具备时间和空间预测功能, 模型变量包括状态、领域空间和时间等离散变量, 公式如下:
$ S_{t+1}=f\left(S_{t}, N\right) $ | (7) |
式中, S为元胞有限、离散的状态集合; t和t+1表示不同时刻; N为元胞的邻域; f为局部空间元胞转换规则[28].
4 结果分析 4.1 土地利用空间格局研究区土地利用以耕地和草地为主. 其中耕地占总面积的50%左右, 主要集中分布在大凌河中、下游, 少量位于上游丘陵地区谷地; 草地占比超过30%, 主要分布在上游丘陵区; 林地占比12%左右; 建设用地约占比5%;湿地和水体面积较小, 约占总面积的1%;裸地面积可忽略不计(图 2, 表 3), 表明研究区土地开发程度处于较高水平.
1) 土地利用结构变化
2000、2010、2020年研究区土地利用结构有较大变化. 2000-2020年期间, 研究区内耕地持续减少, 其中前10年减少了132.85 km2, 后10年减少了1 360.13 km2, 减少速度在加快. 近20年间, 耕地共减少了1 492.97 km2, 面积占比由2000年的53.23%减少为2020年的46.02%. 与此同时, 林地、草地和建设用地面积不断增加, 其中, 草地增加了616.17 km2, 建设用地增加492.23 km2, 林地共增加346.97 km2 (表 3).
2) 土地利用动态度
从土地利用动态度计算结果(表 4)来看, 2000-2010年间, 研究区综合土地利用动态度为43.53%, 单地类中动态度最大的是湿地, 平均每年增加86.32%, 动态度最小的是耕地, 平均每年减少0.12%. 2010-2020年间, 研究区综合土地利用动态度为10.36%, 单地类中动态度最大的是水体, 平均每年增加5.00%, 动态度最小的是湿地, 平均每年减少7.22%. 总体来看, 近20年间, 研究区综合土地利用动态度为15.60%, 处于不断减少的状态. 单地类中动态度最大的是湿地, 平均每年增加8.38%, 动态度最小的是耕地, 平均每年减少0.68%.
3) 土地利用程度
根据式(3)(4)计算了研究区2000、2010和2020年土地利用综合指数, 从数据(表 5)来看, 研究区平均土地利用程度综合指数为261.5, 土地利用程度处于中等偏高水平. 由2000年的262.53变为2020年259.91, 土地利用变化量在持续降低, 表明大凌河流域近20年土地处于衰退期, 且后半期(2010-2020年)衰退速度大于前半期(2000-2010).
4) 土地利用变化矩阵
通过建立土地利用转移矩阵(图 3, 表 6), 可知2000-2010年间, 耕地减少了132.85 km2, 主要转为草地、林地和建设用地. 其中转为草地和林地的耕地面积分别为130.74和54.96 km2, 主要分布于大凌河下游义县-锦州一带; 共有33.34 km2耕地转为建设用地, 主要分布于朝阳市和锦州市等地. 2010-2020年间, 耕地减少面积达1 360.13 km2, 其中有723.68 km2耕地转为草地, 主要发生在大凌河上游和下游; 共有403.15 km2耕地转为林地, 主要位于大凌河上游地区; 另有473.15 km2耕地转为建设用地, 集中分布在阜新市一带. 与此同时, 有少量的草地和林地转化为耕地, 面积分别为143.18和77.74 km2 (如图 3, 表 7).
本研究利用IDRISI软件进行土地利用的模拟预测分析. 第一步在Markov模块设置相对误差为0.15, 分别计算研究区2000-2010年、2010-2020年的土地利用转移概率矩阵; 第二步在CA-Markov模块以2000-2010年的土地利用转移概率矩阵为条件, 设置30 m×30 m元胞、5×5邻域结构、元胞自动机循环次数为10等参数, 模拟预测2020年研究区土地利用分布图; 第三步在CROSSTAB模块对2020年预测值和实际值进行交叉验证, 得到Kappa系数为0.8835, 满足精度要求. 然后重复第二步, 以2020年实际土地利用数据为基础, 模拟预测研究区2030年土地利用格局[29](如图 4a所示).
根据预测结果, 建立2020-2030年土地利用转移矩阵(图 4b, 表 8). 结果表明2030年研究区耕地将进一步减少1 699.92 km2, 其中有630.42 km2耕地转为草地, 主要位于大凌河上游和下游; 共有571.65 km2耕地转为建设用地, 主要发生在朝阳市和义县-锦州一带; 另有454.19 km2转为林地, 主要发生在大凌河中上游地区.
通过分析和预测大凌河流域土地利用格局, 得出如下结论:
1) 2000-2020年大凌河流域土地利用以耕地和草地为主, 二者占比超过80%. 近20年间, 研究区内林地、草地和建设用地面积不断增加, 耕地面积持续减少, 且减少的速度在加快, 前后两个10年分别减少了132.85 km2、1 360.13 km2, 耕地面积占比由2000年的53.23%, 减少为2020年的46.02%. 林地和草地增加主要位于大凌河上游和下游, 建设用地增加主要位于朝阳市和锦州市周边. 这些变化表明区域内"退耕还林还草"工程实施效果明显, 建设用地的扩张规模基本可控.
2) 从土地利用动态度分析来看, 近20年来, 大凌河流域综合土地利用动态度从2000-2010年期间的43.53%, 下降到2010-2020年期间为10.36%, 处于不断减少的状态. 单地类来看, 湿地动态度最大, 平均每年增加8.38%. 表明区域土地利用方式趋于稳定, 湿地等生态系统正在恢复.
3) 2000-2020年间, 大凌河流域土地利用程度综合指数由2000年的262.53变为2020年259.91, 表明区域土地处于衰退期, 土地开发利用程度在逐步降低, 自然生态系统正在恢复之中.
4) 基于2010-2020年土地利用变化概率矩阵, 利用CA-Markov模型预测了大凌河流域在自然条件下2030年土地利用格局, 模拟精度Kappa系数为0.8835, 说明结果可靠. 统计表明, 2030年, 大凌河流域耕地面积将进一步减少1 699.92 km2, 主要变为草地、建设用地和林地. 当然, 土地利用变化预测是一个系统工程, 除了历史土地变化规律外, 还受到人口、经济、环境和土地调控政策等多因素的影响, 今后研究中需综合考虑其他因素, 进一步优化预测模型.
[1] |
刘纪远, 宁佳, 匡文慧, 等. 2010-2015年中国土地利用变化的时空格局与新特征[J]. 地理学报, 2018, 73(5): 789-802. Liu J Y, Ning J, Kuang W H, et al. Spatio-temporal patterns and characteristics of land-use change in China during 2010-2015[J]. Acta Geographica Sinica, 2018, 73(5): 789-802. |
[2] |
侯西勇, 常斌, 于信芳. 基于CA-Markov的河西走廊土地利用变化研究[J]. 农业工程学报, 2004, 20(5): 286-291. Hou X Y, Chang B, Yu X F. Land use change in Hexi Corridor based on CA-Markov methods[J]. Transactions of the Chinese Society of Agricultural Engineering, 2004, 20(5): 286-291. |
[3] |
陈龙泉, 郑海金. 基于Markov-CA的土地利用/土地覆盖变化动态模型研究[J]. 测绘信息与工程, 2004, 29(1): 36-38. Chen L Q, Zheng H J. Dynamic models of land use/land cover changes based on Markov-CA[J]. Journal of Geomatics, 2004, 29(1): 36-38. |
[4] |
汤洁, 汪雪格, 李昭阳, 等. 基于CA-Markov模型的吉林省西部土地利用景观格局变化趋势预测[J]. 吉林大学学报(地球科学版), 2010, 40(2): 405-411. Tang J, Wang X G, Li Z Y, et al. The tendency forecast on land use landscape pattern change in western Jilin Province based on CA-Markov model[J]. Journal of Jilin University (Earth Science Edition), 2010, 40(2): 405-411. |
[5] |
马梓策, 于红博, 张巧凤, 等. 基于CA-Markov模型的锡林河流域土地利用变化及预测分析[J]. 西北林学院学报, 2020, 35(2): 169-177. Ma Z C, Yu H B, Zhang Q F, et al. Land use change and forecast analysis in the Xilin River Basin based on CA-Markov model[J]. Journal of Northwest Forestry University, 2020, 35(2): 169-177. |
[6] |
张永民, 赵士洞, Verburg P H. CLUE-S模型及其在奈曼旗土地利用时空动态变化模拟中的应用[J]. 自然资源学报, 2003, 18(3): 310-318. Zhang Y M, Zhao S D, Verburg P H. CLUE-S and its application for simulating temporal and spatial change of land use in Naiman Banner[J]. Journal of Natural Resources, 2003, 18(3): 310-318. DOI:10.3321/j.issn:1000-3037.2003.03.008 |
[7] |
涂小松, 濮励杰, 吴骏, 等. 基于SLEUTH模型的无锡市区土地利用变化情景模拟[J]. 长江流域资源与环境, 2008, 17(6): 860-865. Tu X S, Pu L J, Wu J, et al. Using SLUETH model to simulate land use change scenarios in Wuxi district[J]. Resources and Environment in the Yangtze Basin, 2008, 17(6): 860-865. |
[8] |
陈海, 王涛, 梁小英, 等. 基于MAS的农户土地利用模型构建与模拟——以陕西省米脂县孟岔村为例[J]. 地理学报, 2009, 64(12): 1448-1456. Chen H, Wang T, Liang X Y, et al. Simulation and application of household's LUCC based on a multi-agent system: A case study for Mengcha Village of Mizhi County of Shaanxi Province[J]. Acta Geographica Sinica, 2009, 64(12): 1448-1456. |
[9] |
Muller M R, Middleton J. A Markov model of land-use change dynamics in the Niagara Region, Ontario, Canada[J]. Landscape Ecology, 1994, 9(2): 151-157. |
[10] |
Subedi P, Subedi K, Thapa B. Application of a hybrid cellular automaton-Markov (CA-Markov) model in land-use change prediction: A case study of Saddle Creek Drainage Basin, Florida[J]. Applied Ecology and Environmental Sciences, 2013, 1(6): 126-132. DOI:10.12691/aees-1-6-5 |
[11] |
Halmy M W A, Gessler P E, Hicke J A, et al. Land use/land cover change detection and prediction in the north-western coastal desert of Egypt using Markov-CA[J]. Applied Geography, 2015, 63: 101-112. |
[12] |
于涛, 沈浩, 仲嘉亮. 基于CA-Markov模型的新疆克州土地利用动态模拟研究[J]. 新疆环境保护, 2008, 30(1): 11-14. Yu T, Shen H, Zhong J L. Dynamic analogue research of land utilization of Kizilsu Kirgiz Autonomous Prefecture Xinjiang base on CA-Markov model[J]. Environmental Protection of Xinjiang, 2008, 30(1): 11-14. |
[13] |
杨国清, 刘耀林, 吴志峰. 基于CA-Markov模型的土地利用格局变化研究[J]. 武汉大学学报(信息科学版), 2007, 32(5): 414-418. Yang G Q, Liu Y L, Wu Z F. Analysis and simulation of land-use temporal and spatial pattern based on CA-Markov model[J]. Geomatics and Information Science of Wuhan University, 2007, 32(5): 414-418. |
[14] |
施成超, 周英, 刘滢, 等. 基于CA-Markov模型的云南省德宏州土地利用变化和预测研究[J]. 西部林业科学, 2022, 51(1): 102-109. Shi C C, Zhou Y, Liu Y, et al. Modeling and prediction of land use change in Dehong Prefecture of Yunnan Province based on CA-Markov model[J]. Journal of West China Forestry Science, 2022, 51(1): 102-109. |
[15] |
张晓娟, 周启刚, 王兆林, 等. 基于MCE-CA-Markov的三峡库区土地利用演变模拟及预测[J]. 农业工程学报, 2017, 33(19): 268-277. Zhang X J, Zhou Q G, Wang Z L, et al. Simulation and prediction of land use change in Three Gorges reservoir area based on MCE-CA-Markov[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(19): 268-277. |
[16] |
黎云云, 畅建霞, 王义民, 等. 基于CA-Markov模型的黄河流域土地利用模拟预测研究[J]. 西北农林科技大学学报(自然科学版), 2020, 48(12): 107-116. Li Y Y, Chang J X, Wang Y M, et al. Land use simulation and prediction in the Yellow River Basin based on CA-Markov model[J]. Journal of Northwest A&F University (Natural Science Edition), 2020, 48(12): 107-116. |
[17] |
王禹浪, 孙军, 王文轶. 大、小凌河流域的古代文明与历史文化[J]. 黑龙江民族丛刊, 2008(1): 83-88. Wang Y L, Sun J, Wang W Y. Ancient civilization and culture along the Big and Small Ling River[J]. Heilongjiang National Series, 2008(1): 83-88. |
[18] |
邸志强, 金洪涛, 苗英, 等. 辽西蓄水构造及找水方向[J]. 地质与资源, 2007, 16(2): 112-115. Di Z Q, Jin H T, Miao Y, et al. Water storage structure and water search target in western Liaoning Province[J]. Geology and Resources, 2007, 16(2): 112-115. |
[19] |
裴亮, 陈晨, 戴激光, 等. 基于马尔科夫模型的大凌河流域土地利用/覆被变化趋势研究[J]. 土壤通报, 2017, 48(3): 525-531. Pei L, Chen C, Dai J G, et al. Research on forecast trend of land use and land cover change in Daling River Basin based on Markov model[J]. Chinese Journal of Soil Science, 2017, 48(3): 525-531. |
[20] |
魏晶, 王涌翔, 吴钢, 等. 辽西大凌河流域土地利用变化及驱动力分析[J]. 生态环境, 2006, 15(3): 559-563. Wei J, Wang Y X, Wu G, et al. Driving forces analysis of land utilization in Dalinghe River watershed of western Liaoning Province during 1987-2002[J]. Ecology and Environment, 2006, 15(3): 559-563. |
[21] |
陈军, 陈晋, 廖安平, 等. 全球地表覆盖遥感制图[M]. 北京: 科学出版社, 2016. Chen J, Chen J, Liao A P, et al. Remote sensing mapping of global land cover[M]. Beijing: Science Press, 2016. |
[22] |
王秀兰, 包玉海. 土地利用动态变化研究方法探讨[J]. 地理科学进展, 1999, 18(1): 81-87. Wang X L, Bao Y H. Study on the methods of land use dynamic change research[J]. Progress in Geography, 1999, 18(1): 81-87. |
[23] |
朱会义, 李秀彬. 关于区域土地利用变化指数模型方法的讨论[J]. 地理学报, 2003, 58(5): 643-650. Zhu H Y, Li X B. Discussion on the index method of regional land use change[J]. Acta Geographica Sinica, 2003, 58(5): 643-650. |
[24] |
杨佳佳, 张一鹤, 冯雨林, 等. 松嫩平原东部土地利用时空动态变化分析[J]. 地质与资源, 2020, 29(6): 627-634, 602. Yang J J, Zhang Y H, Feng Y L, et al. Spatiotemporal dynamic change analysis of land use in eastern Songnen Plain[J]. Geology and Resources, 2020, 29(6): 627-634, 602. |
[25] |
张安定, 李德一, 王大鹏, 等. 山东半岛北部海岸带土地利用变化与驱动力——以龙口市为例[J]. 经济地理, 2007, 27(6): 1007-1010. Zhang A D, Li D Y, Wang D P, et al. Analysis of land use dynamic changes and its driving forces in the north of Shandong Peninsula: Taking Longkou as an example[J]. Economic Geography, 2007, 27(6): 1007-1010. |
[26] |
张俊平, 李净. 基于CA-Markov模型的甘州区土地利用变化预测分析[J]. 中国农学通报, 2017, 33(4): 105-110. Zhang J P, Li J. Prediction and analysis of land use change based on CA-Markov model in Ganzhou District[J]. Chinese Agricultural Science Bulletin, 2017, 33(4): 105-110. |
[27] |
刘淑燕, 余新晓, 李庆云, 等. 基于CA-Markov模型的黄土丘陵区土地利用变化[J]. 农业工程学报, 2010, 26(11): 297-303. Liu S Y, Yu X X, Li Q Y, et al. Land use change in loess hilly region based on CA-Markov model[J]. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(11): 297-303. |
[28] |
何丹, 周璟, 高伟, 等. 基于CA-Markov模型的滇池流域土地利用变化动态模拟研究[J]. 北京大学学报(自然科学版), 2014, 50(6): 1095-1105. He D, Zhou J, Gao W, et al. An integrated CA-Markov model for dynamic simulation of land use change in Lake Dianchi watershed[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2014, 50(6): 1095-1105. |
[29] |
席梅竹, 赵中秋, 吴攀升, 等. 基于改进CA-Markov模型的滹沱河流域山区段土地利用变化模拟及预测[J]. 西北林学院学报, 2021, 36(4): 150-158. Xi M Z, Zhao Z Q, Wu P S, et al. Changes and predictions of land use in mountain section of the Hutuo River Basin based on improved CA-Markov model[J]. Journal of Northwest Forestry University, 2021, 36(4): 150-158. |