Temporal and spatial dynamic monitoring of land desertification and landscape pattern evolution characteristics in the Xiaoluan River Basin

围场县(40°35′~42°40′N，116°32′~118°14′E)隶属河北省承德市，东、西、北三面分别与内蒙古赤峰市、克什克腾旗、多伦县毗邻，西南分别与丰宁满族自治县、隆化县接壤。小滦河流域位于围场县西部，流域面积约1 655.0 km²。作为滦河上游的主要支流，流域气候类型属于季风型高原山地气候，年均气温为3.3 ℃，年均降水量为445.0 mm，地势东北高，西南低，包括坝上草原区和冀北山地区，海拔为900~2 000 m。植被资源丰富，森林多为天然次生林和人工林。土壤类型有灰色森林土、潮土、风沙土、沼泽土等。由于大量沙丘存在于地表，抗蚀力弱，极易产生侵蚀，导致该区域土地沙化严重^[21−22]。自20世纪60年代开始，流域内大面积实施人工造林，生态环境有所改善。随着周边多个景点和公路的开发，小滦河流域成为围场县生态旅游观光热点区域，致使人类活动干扰强度较大，耕地、草场均存在不同程度退化等生态问题。

土地利用类型来源于中国科学院资源环境数据中心(http://www.resdc.cn)；地形数据为数字高程模型(DEM)，来源于地理空间数据云(http://www.gscloud.cn)，空间分辨率均为30 m。基于研究区边界完成掩膜、拼接和裁剪，得到各期荒漠化差值指数，采取ArcGIS 10.6、R 4.1.0和Origin 2021对数据进行统计分析。

Landsat数据源于谷歌地球引擎(Google Earth Engine, GEE)提供的3个时段Landsat TM5 (2000、2010年)及Landsat-8 OLI (2020年)数据集为基础，筛选云量小于10%的生长期影像，2000年共调用8景Landsat TM影像，2010年共调用15景Landsat TM影像，2020年共调用20景Landsat OLI影像，采用均值合成法计算各时期NDVI、Albedo，并构建Albedo-NDVI特征空间。

地表反照率是指地面反射的太阳辐射与到达地面的太阳辐射之比，与土壤颜色、土壤湿度、太阳高度角和天气状况等有关^[23]。采用LIANG^[24]的Landsat数据反演模型对研究区进行估算，公式为：

式(1)中：A为地表反照率；$ {\rho }_{\mathrm{T}\mathrm{M}1}\mathrm{、}{\rho }_{\mathrm{T}\mathrm{M}3}\mathrm{、}{\rho }_{\mathrm{T}\mathrm{M}4}{\mathrm{、}\rho }_{\mathrm{T}\mathrm{M}5}\mathrm{、}{\rho }_{\mathrm{T}\mathrm{M}7} $为不同波段的权重值。

采用近红外波段和红光波段两者的反射率计算NDVI，计算公式为：

式(2)中：N_DVI表示归一化植被指数；$ {\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}} $、$ {\rho }_{\mathrm{R}\mathrm{E}\mathrm{D}} $分别表示近红外波段和红光波段所对应的反射率值。由于 N_DVI与Albedo的量纲与范围不同，分别计算两者的最大值和最小值，并进行归一化处理。

研究表明：Albedo-NDVI特征空间分布图一般呈现梯形或者三角形^[25]，可将高植被覆盖、部分覆盖和全裸露很好地区分。

式(3)~(5)中：D_DI为荒漠化差值指数；$ A $为地表反照率；N_DVI为归一化植被指数；$ k、b $为方程拟合系数；$ \alpha $值由拟合曲线的斜率($ k $)决定。研究表明^[26]：基于Jenk最优化法的分界点使各分组方差之和最小，将荒漠化差值指数(D_DI)分为5个区间：极重度沙化(D_DI≤1.110)、重度沙化(1.110＜D_DI≤1.453)、中度沙化(1.453＜D_DI≤1.942)、轻度沙化(1.942＜D_DI≤2.472)和非沙化土地(D_DI＞2.472)。

土地沙化动态度($ K $)可定量描述某区域某研究时段内各类沙化土地面积的变化速度，对沙化区域性差异及预测未来趋势有一定优势。单一沙化程度动态度公式如下^[27]：

式(6)中：$ K $为研究时段内某一沙化程度动态度；$ {U}_{\mathrm{a}}、{U}_{\mathrm{b}} $分别表示研究初期与末期某种沙化程度的面积；$ T $为研究时段。$ K $＜0表明该类沙化土地处于缩减状态，$ K $＞0表明该类沙化土地处于扩张状态。当$ T $定义为年时，$ K $值即区域某程度沙化土地的年变化率。

利用回归分析，研究小滦河流域沙化动态随时间变化的趋势^[28]，计算公式如下：

式(7)中：S为荒漠化差值指数(D_DI)变化趋势的斜率；n为监测年数；D_DIi为第$ i $年对应的荒漠化差值指数。根据荒漠化差值指数变化趋势划分为严重退化(S≤−0.20)、轻微退化(−0.20＜S≤−0.10)、基本不变(−0.10＜S≤0.10)、轻微改善(0.10＜S≤0.25)、明显改善(S＞0.25)。

转移矩阵可对一定时期单元内不同沙化程度的转移状态做定量分析，反映其变化的方向和特征^[29]，更好地揭示沙化土地的时空演化过程。数学表达形式为：

式(8)中：${S}_{hj}$表示初期对应h类沙化土地转变为末期$ j $类沙化土地的面积；$ m $表示沙化土地相应程度的数目。

景观格局指数以斑块作为基本单元用以分析各景观类别的结构组成及空间特征^[30]，根据景观生态学理论，结合小滦河流域实际情况，选取面积周长维数、蔓延度指数、凝聚度指数、分离度指数、香农多样性指数、香农均匀度指数、聚集度指数等7类景观格局指数，运用Fragstats 4.2软件对小滦河流域3期土地沙化景观格局特征进行分析。

研究表明^[25]：高反照率线(AC)与沙化过程在Albedo-NDVI特征空间中的轨迹较为接近，选取垂直于高反照率线的直线将不同沙化土地进行区分，荒漠化差值指数中斜率$ \alpha =-1/k $，即荒漠化差值指数表现为Albedo-NDVI特征空间垂直于高反照率线的各分割直线位置，代表不同沙化土地在特征空间中的地表水热组合及动态变化差异(图1)。基于3期数据构建Albedo-NDVI特征方程，衍生出荒漠化差值指数，可用于土地沙化程度分级^[31]。基于ArcGIS和RStudio，随机选取研究区内1 256个样本点进行回归分析(以2000年为例)。从图2可见：Albedo与NDVI具有明显的线性负相关，R²为0.515 1。

Figure 1.  Albedo-NDVI feature space

Figure 2.  Albedo-NDVI feature space scatter plot and linear equation

基于Google Earth Pro历史影像数据，随机选取300个样本点对3期沙化土地进行精度验证，结合2021年8月野外调查数据进行目视解译，以此判断验证点沙化程度类别，生成2000、2010和2020年沙化土地分类精度评价混淆矩阵(表1)。由表1可见：总体精度为79.33%~84.67%，表明基于Albedo-NDVI特征空间方法在沙化程度区分中结果良好，分类精度满足实验要求。

由图3可知：小滦河流域极重度、重度沙化土地集中分布在御道口镇、御道口牧场等。中度沙化土地分布范围较广，涉及中、上游大部分地区。而下游以轻度沙化和非沙化为主，包括老窝铺乡、西龙头乡、南山嘴乡；小滦河流域中、下游沙化情况持续好转，而上游地区在2010年后沙化土地有所增加。

Figure 3.  Spatial and temporal distribution map of land desertification in Xiaoluanhe River Basin based on desertification difference index

从表2可以看出：2000—2010年小滦河流域沙化土地面积变化幅度明显，极重度沙化土地面积减少了165.31 km²，重度沙化土地面积减少了171.21 km²，动态度分别为−6.88%、−4.18%；轻度沙化、非沙化土地面积分别增加了186.47、219.18 km²，其中非沙化土地动态度达6.98%。从空间分布看出，下游地区大范围沙化土地转变为非沙化土地，中游地区由极重度和重度沙化连片分布转变为斑块状零散分布(图3)。2010—2020年，小滦河流域极重度、重度沙化土地面积进一步减小，其中极重度、重度沙化土地变化速率相对较大，动态度分别为−3.40%、−2.44%，不同沙化程度土地动态度均有缩减。极重度、重度沙化土地面积持续缩小，轻度沙化土地面积持续增加，但变化幅度较2000—2010年均有所减少。表明小滦河流域沙化状况整体有所好转，极重度、重度沙化面积和占比有明显减少，上游地区的沙化趋势不容忽视。

2000—2020年，小滦河流域土地沙化总体得到改善，不同程度沙化土地的变化趋势和幅度各不相同。如图4A所示：在2000—2010年，极重度转为重度沙化土地面积为95.36 km²，占极重度沙化土地面积的39.69%，重度转为中度沙化土地面积为164.59 km²，占重度沙化土地面积的40.18%，中度转为轻度沙化土地面积为201.08 km²，占中度沙化土地面积的48.65%，轻度转为非沙化土地面积为129.10 km²，占轻度沙化土地面积的46.33%，表明前期小滦河流域沙化土地改善幅度较大，以相邻土地沙化发生转移。如图4B所示：2010—2020年，小滦河流域中度沙化转为轻度沙化面积达137.82 km²，轻度转为非沙化土地面积103.88 km²，重度转中度沙化土地面积为100.80 km²，但存在非沙化转向轻度沙化，轻度转向中度沙化土地的现象，转移面积分别为118.06 、81.14 km²。相较2010年前，小滦河流域改善面积有所缩减，且293.81 km²的非沙化、轻度沙化等土地向沙化程度严重的类型转变，流域部分区域的沙化状况恶化加剧。

Figure 4.  Matrix of desertification land change of the Xiaoluanhe River Basin

如图5所示：2000—2020年，流域呈明显改善和轻度改善状态的土地面积分别占总面积的17.18%和24.28%，呈严重退化和轻度退化的土地面积之和不足10%。表明20 a间小滦河流域的沙化状况总体有明显改善，主要分布于北部塞罕坝林场及小滦河流域下游，而以轻度退化为主的土地沙化加剧状况分别集中在御道口牧场、御道口镇等区域及周边。

Figure 5.  Changes in the spatial and temporal pattern of land desertification of the Xiaoluanhe River Basin from 2000 to 2020

如表3所示：20 a间小滦河流域凝聚度指数呈先减后增趋势，总体为增加趋势，说明不同沙化景观程度间同类斑块的自然连通度提高。蔓延度指数、聚集度指数均为连年上升，前者增加了8.215，聚集度指数增加了3.589，且2010年前增幅超过后期，说明小滦河流域沙化景观集聚程度逐渐提升，流域景观中小斑块数量减少，在空间分布上趋向于集中，流域景观破碎化程度持续降低。面积周长维数持续减小，流域沙化景观的几何形状复杂程度及分布情况趋于简单。同时，分离度指数呈微弱的减少趋势，景观沙化土地分散程度减小，香农多样性和均匀度指数均有降低，分别降低了0.182、0.113，在2000—2010年下降趋势较为明显，表明小滦河流域沙化景观的多样性和异质性有所减缓，破碎化程度随景观类型复杂程度的减少而变小。

人为因素驱动下的土地利用类型对区域土地沙化有着改善或者恶化的作用，从而引起土地沙化并发生改变。由表4可见：2000—2010年，耕地的极重度和重度沙化面积共减少了76.33 km²，占同期沙化土地改善面积的15.86%，林地的极重度、重度和中度沙化面积减幅最大，达225.52 km²，未利用土地的沙化面积减少82.08 km²。10 a内草地面积增加最多，达83.25 km²，草地的沙化土地改善面积占草地总面积的20.06%。在2010—2020年，草地为小滦河流域沙化改善的主要地类。其中，耕地的极重度、重度和中度沙化土地面积持续减少，林地非沙化土地面积减少62.26 km²，总面积未有明显改变，而重度、中度和轻度沙化面积均有所增加。这一时期未利用土地面积变化较小，且各类沙化土地均有好转。后10 a间草地面积变幅较小，但其土地沙化好转面积最为突出，其中极重度、重度和中度沙化土地面积减少63.90 km²，占同期沙化土地改善面积的57.70%。

小滦河流域海拔为910~1 868 m，结合沙化土地分布，以200 m为单位划分5个等级，对不同海拔的沙化土地面积进行分析。如图6所示：20 a间小滦河流域各类沙化土地主要分布于海拔1 300~1 500 m。前期土地沙化在海拔1 300~1 500、1 500~1 700和1 100~1 300 m缩减最为明显，其中，海拔1 100~1 300 m以极重度沙化土地减轻为主，1 300~1 500 m以重度沙化减轻为主，1 500~1 700 m以中度沙化面积减少为主，其余海拔范围的沙化好转面积较小。2010—2020年，不同海拔之间极重度、重度沙化土地面积有轻微减少，在海拔1 300~1 500 m，非沙化、轻度沙化土地面积增幅最大。海拔1 500~1 700 m非沙化土地面积减少了57.64 km²，而轻度沙化土地面积增加了30.55 km²，中度沙化土地面积增加了21.49 km²；且分布在海拔1 100~1 300、1 500~1 700 m的土地沙化逆转趋势有所加重，今后需加强治理。

Figure 6.  Distribution of different desertified land areas in Xiaoluanhe River Basin from 2000 to 2020

2000—2020年，小滦河流域各土地利用的土地沙化均有不同程度改善。前期(2000—2010年)，土地利用类型的变化面积从大到小依次为草地、未利用土地、林地、耕地。沙化土地改善面积从大到小依次为林地、草地、未利用土地、耕地，林地最为突出。不同海拔范围的沙化改善程度有所不同。后期(2000—2010年)，流域土地利用间的转变幅度较小，主要为耕地、林地转为未利用土地。不同地类的沙化改善速率有所降低，从大到小依次为草地、耕地、未利用土地、林地，以草地、林地作为土地沙化逆向转变的主要地类，沙化逆转趋势加重。

本研究表明：小滦河流域土地沙化状况总体向改善的方向发展，主要以轻度、非沙化土地为主，部分土地存在退化趋势，该结论与贾双竹^[20]的研究结果相符。从海拔来看，周智勇^[32]指出：中高海拔山地、高原的植被覆盖度高，沙化程度改善明显，以轻度和非沙化土地为主。沙化严重区集中在中低海拔区，对应区域植被覆盖较低，是今后需重点关注及治理的区域。此外，依据围场县多年植被景观格局变化特征，对流域沙化景观格局指数进行分析发现：随着流域土地沙化状况的明显改善，景观中小斑块数量减少，连通性有所提高，景观破碎化程度有所减缓。

联系围场县实际情况可知：自20世纪90年代开始实施农业生态工程，退化牧草地、未利用土地被用于林业建设，使林地成为该区域的主要用地类型。从1998年开始，“三北”防护林工程、京津风沙源治理、坝上生态农业建设项目等相继启动，使围场县作为重点治理区域的沙化情况有所好转^[33]。适应当地立地条件的林地，对沙化土地的改善具有促进作用，森林根系发达、自身储存水份等特点使该区域土壤中水分较为充足，林冠、枝干可有效减缓风速，通过截留降水进行雨水再分配，凋落物层进一步减轻对土壤的冲刷和侵蚀。研究区通过人工、固沙造林和稀疏林地补植等措施，使沙化有明显改善，表明森林在土地沙化、水土流失等方面作用突出。草地作为小滦河流域沙化改善的主要用地类型，应建设草地保护与修复项目示范区，同时在退化草地区域人工补种优质牧草^[34]。这主要因为草地多为连片分布，根系在土壤密集交错，可使沙化程度降低。

后期(2010—2020年)，小滦河流域土地沙化恶化趋势加剧，御道口镇东北部的草地沙化程度存在加重趋势，且中上游存在着林地被转为耕地、牧场南部的林地转变为草地等现象，均导致流域非沙化、轻度沙化等地块恶化。再者，近年来，由于沙化土地改善与当地经济发展并存，尤其是旅游业的发展使流域受人为因素干扰较大，分布于中低海拔区的建设用地、水域面积有轻微增加，进而对流域沙化动态变化有所影响。同时，草场超载放牧、土地退化、土层浅薄导致水土流失，流域局部地区沙化趋势有所加重。今后需结合小滦河流域实际，加强以草地为重点的沙化改善措施，坚持实施退耕还林还草，以中低海拔区作为今后林草建设的重点区域，优化耕地空间格局及治理，将区域生态旅游发展与居民受益紧密联系，同时对农村环境进行综合改善与整治，实现小滦河流域的可持续发展。

2000—2020年，小滦河流域土地沙化总体呈好转趋势。其中，极重度、重度沙化土地面积分别减少约80%和56%，集中分布于御道口镇及牧场周边，轻度沙化和非沙化分布于下游地区，后期沙化土地恶化趋势加重。流域沙化斑块在空间分布上趋向于集中，景观破碎化程度持续降低。不同地类的土地沙化均有改善，林地为前期土地沙化改善的主要地类，后期以草地、耕地为主，耕地的沙化状况持续好转，林地出现逆向沙化现象。低海拔区以极重度沙化土地减轻为主，中海拔区以重度沙化减轻为主，高海拔区以中度沙化减少为主。