党怡雯， 刘景刚， 张家琛， 陈亚恒

(1.河北农业大学国土资源学院，河北 保定 071000； 2.天津市滨海新区土地发展中心，天津 300457)

土地是人类生存、生产、建设不可或缺的要素和载体，在全球气候与环境变化的研究领域中，土地利用/土地覆盖变化 (LUCC) 是其中的重点内容。它直接反映人类活动作用于地球陆地表面生态系统产生的结果，是人类社会经济活动行为与自然生态过程交互和链接的纽带[1-2]。2005年全球土地计划(Global land Project，GLP)将土地利用/土地覆盖动态监测与模拟列为核心内容之一[3]。地形作为重要的地理环境要素，对区域内光照、降水等气候资源以及养分等土壤资源的再分配起着主导作用，直接影响着生态系统的物质流动和能量循环，制约着土地利用的方式和程度，进而影响区域土地利用的空间格局[4-5]，也因此成为了山区贫困化的主要自然因素之一[6]。

近年来，探究不同地形梯度下土地利用方式、土地景观格局的转变已经趋于成熟。学者们将目光聚焦于不同研究视角：城市、山区、流域、森林等[6-11]以及不同研究尺度：省、市、县层面[7,12-13]。在此基础上，也有众多学者运用景观生态学原理，以空间粒度为切入点，对景观格局进行深入探究[14-15]。研究结果表明景观格局的演变与热环境效应也有密不可分的联系[16-18]。但目前鲜有人针对地形因素，将景观格局的演变与热环境效应结合起来进行系统分析。

河北省阜平县地形起伏较大，地形梯度明显，地势由东南向西北逐渐升高，自然生物分布也呈现一定的垂直地带性。区内生态环境较为脆弱，水土流失严重，裸地面积广。县内景观类型多样，其热量资源的分布也受到地形、地势的影响，深刻影响着土地利用的方式和格局。为了加快推进乡村振兴战略的实施，需要利用好阜平县的地形特色优势，掌握资源脉络，因地制宜进行生产、旅游建设。本研究以河北省阜平县为研究区，探究阜平县景观格局和热环境效应的地形梯度时空变化特征以及二者的相关性，揭示地形因子对景观格局和热环境效应的影响，为该区域推进乡村振兴战略、生态文明建设，景观格局的优化以及农业生产、旅游业的发展提供一定的理论依据。

1 研究方法与数据来源

1.1 研究区概况

阜平县位于河北省中西部，保定地区西部，地处太行山中北部东麓，属暖温带半干旱和半湿润地带。下辖4镇9乡，209个行政村，1 208个自然村，处于东经113°45′～114°32′、北纬38°39′～39°08′，总面积达2 494 km2。全境地形复杂多样，山峦起伏，由于海拔高度的差异热量资源呈现阶梯状分布，西部深山区气候寒冷，温度较低，东部浅山区温度较高。境内流域面积200 km2以上的中小河流6条，大都是季节性河流。年平均气温12.7 ℃，无霜期约191 d，地区年降水量为440～650 mm。境内景观类型多样，矿产资源丰富，动植物资源有800余种。

1.2 数据来源及处理

本研究所用数据均来源于地理空间数据云，包含分辨率为30 m的DEM数字高程数据和Landsat TM和Landsat OLI_TIRS的2009年7月2日、2014年7月25日、2019年7月23日3期遥感影像数据，成像时间均为7月，云量均低10%，成像质量较好，且无极端异常天气。通过监督分类与人机交互式解译生成，同时，进行野外实地考察验证进行修正，将土地利用类型划分为耕地、园地、林地、草地、建设用地、水体、其他土地，共计7类。

1.3 地形位指数和分布指数

单一分析地形因子不足以科学体现地形因素对土地利用方式和景观格局的影响程度。因此，引入地形位指数，综合高程与坡度双因子对地形特征进行综合描述，其计算公式如下[19]：

T=ln[(E/E0+1)×(S/S0+1)]

(公式1)

式中，T为地形位指数，E、S分别表示研究区内任意点的高程值和坡度值，E0、S0分别表示计算单元内的平均高程和平均坡度。高程越高、坡度越大的区域地形位指数越大，反之则越小。

采用分布指数来描述研究区各景观类型在地形梯度上的空间分布情况，确定优势分布区间，其计算公式如下[20]：

P=(Sie/Si)/(Se/S)

(公式2)

式中，P为分布指数，Sie表示在e级地形内第i种景观类型的分布面积，Si表示第i种景观类型总面积，Se表示研究区内e地形区间的总面积，S表示研究区总面积。当P>1时，表示该类型景观在e级地形上呈优势分布，P值越大，表示分布优势越明显[21]。

1.4 土地利用图谱分析法

土地利用图谱法是揭示土地利用时空变化过程的重要方法[19,22-23]。借助ArcGIS10.2软件进行空间叠加，得到变化信息图谱，可直观反映各个景观类型的演变过程，其计算公式：

Y=(G1)×10n-1+(G2)×10n-2+…+(Gn)×10n-n

(公式3)

式中，Y表示变化图谱，n为土地利用变化时期数目，Gn表示第n期数据中的景观类型。

1.5 景观格局指数

景观格局指数反映景观格局的空间结构信息。剔除指标的冗余度，本研究选取景观水平上的斑块密度(PD)、最大斑块指数(LPI)、景观形状指数(LSI)、蔓延度(CONTAG)、边缘密度(ED)、香农多样性指数(SHDI)共6个指标，应用Fragstats4.2软件进行计算，可较为全面地反映景观格局的特征。

1.6 地表温度反演

借鉴覃志豪等[24]提出的单窗算法对Landsat TM 5和Landsat OLI 8的热红外波段进行地表温度反演。对于Landsat TM 5(波段6)数据，计算公式如下：

T6=1 260.560 0/ln[1+60.776 0/(0.123 8+0.005 632 56DNTM6)]

(公式4)

式中，T6为亮度温度值，单位为K；DNTM6为TM6的像元DN值，取值范围为[0,255]。

LST6={-67.355 351(1-C6-D6)+[0.458 606(1-C6-D6)+C6+D6)]T6-D6Ta}/C6

(公式5)

Ta=17.976 90+0.917 15T0

(公式6)

C6=ε6τ6

(公式7)

D6=(1-ε6)[1+(1-ε6)τ6]

(公式8)

式中，LST6为地表温度，单位为K，地表温度取值范围通常在0～70 ℃；Ta表示大气平均作用温度；C6、D6为中间变量；T0为地面附近气温，单位为K；τ6表示大气透射率；ε6表示地表辐射率。

对于Landsat OLI 8(波段10)数据，计算公式如下：

Lλ=MLDN+AL

(公式9)

式中，Lλ为波段λ的辐射亮度值，ML为波段λ的增益值，AL为波段λ的偏移值。

Tk=K2/ln(K1/Lλ+1)

(公式10)

式中，Tk表示亮温，是运算过程中一个重要参数；K1和K2为定标常数，Lλ为上述辐射定标后的影像波段。其中K1=774.89 W/(m2·sr·μm)，K2=1 321.08 K。

LST10={-67.355 351(1-C10-D10)+[0.458 606(1-C10-D10)+C10+D10)]Tk-D10Ta}/C10

(公式11)

C10=ε10τ10

(公式12)

D10=(1-ε10)[1+(1-ε10)τ10]

(公式13)

式中，C10、D10为中间变量，ε10表示地表辐射率，τ10表示大气透射率。

t=LST-273.15

(公式14)

式中，LST为开氏温度，单位为K；t为摄氏温度，单位为℃。

2 结果与分析

2.1 阜平县景观类型空间分布的地形梯度效应

2.1.1 不同等级地形位指数的分布 借助ArcGIS10.2软件提取研究区的高程和坡度，进行叠加分析，得到阜平县地形位指数图(图1)。参考路昌等[9]、孙丕苓等[19]的研究结果，结合阜平县实际地貌特征，采用等间隔的方法进行重分类[9,25-27]，将地形位指数分为10级，界定低、中低、中高和高地形位指数4个不同梯度[20]。研究区内地形梯度以中低等级为主，低等级和中高等级次之，高等级面积仅占0.05%。在地形位指数等级分布上，集中分布在3～6级，其中4级地形区面积占比18.21%，5级地形区和3级地形区次之，分别占17.66%和17.28%(表1)。

图1 阜平县地形梯度分布图Fig.1 Topographic gradient distribution map of Fuping County

表1 阜平县地形位指数分级及其面积占比

2.1.2 不同景观类型在地形梯度上的分布特征 景观类型的空间分布呈现出地带差异性，不同地形位等级上各景观类型优势分布区间有所不同(图2)。时间序列上，各景观类型的优势分布区间不变，且总体曲线形态基本相似。空间序列上，耕地、草地、建设用地、园地优势区段分布在1～3级地形区，水体优势分布区段分布在1～2级地形区，分布指数跨度大，一级地形区优势明显，说明其分布严格受地形位的限制，水体的优势分布区域也体现出研究区河流和水库的分布。其他土地优势分布区段分布在4～6级地形区，林地优势分布区段分布在6～10级地形区。阜平县全境森林覆盖率较高，林地对各种地形具有较强的适应性，高程越高、坡度越大的地区，林地的分布优势越为明显。

图2 各景观类型在地形梯度上的分布指数Fig.2 Distribution index of each landscape type on topographic gradient

2.2 景观格局时空变化的地形梯度效应

2.2.1 景观格局变化图谱 景观格局变化图谱可深入剖析连续时间序列内景观类型的动态变化过程，明确地类转换的先后次序，从侧面反映人类活动的作用结果及程度。2009-2019年2个等时间间隔内，耕地、建设用地均以转入为主，其余土地利用类型以转出为主。基于土地利用转移矩阵(表2、表3)，对景观格局图谱进行分析。在ArcGIS10.2中进行空间叠加处理，得到阜平县连续2个时间间隔内的土地景观格局演变轨迹(图3)。在景观类型图谱变化模式中，稳定型面积占比达到98.489%，其余类型面积占比极小(表4)。稳定型中最大变化图谱类型为其他土地-其他土地-其他土地，其次为建设用地-建设用地-建设用地。其他土地为研究区内的主导景观类型，占据重要地位。阜平县境内存在大量裸地，占其他土地景观类型的绝大部分，未来应该进一步合理开发利用。前期变化型中，有其他土地-耕地-耕地、草地-建设用地-建设用地、耕地-建设用地-建设用地、其他土地-建设用地-建设用地，其中其他土地-耕地-耕地占比57.73%，说明在该时段内阜平县的耕地需求量增加。后期变化型中，以其他土地-其他土地-耕地和草地-草地-耕地为主，其中，其他土地-其他土地-耕地所占该类型比重较大，达57.42%，凸显了耕地的重要战略地位，进一步体现了粮食安全的重要性，阜平县其他土地仍有较大利用空间。反复变化型中主要以水体-建设用地-水体为主，占比95.68%，水源保护对区域的可持续发展起着重要的作用，反映出当地政府管控建设用地无序随意扩增政策的推进。持续变化型面积占比仅0.001%，主要由耕地-其他土地-建设用地和林地-其他土地-建设用地构成。

表2 2009-2014年土地利用转移矩阵

表3 2014-2019年土地利用转移矩阵

2.2.2 景观格局信息图谱的地形梯度 将景观类型变化图谱与地形位指数进行空间叠加处理，景观格局信息图谱的地形梯度效应如图4所示。稳定型变化图谱对于地形等级高的地区具有较好的适应性，优势区间为4～10级区，其分布指数在地形区间内始终保持相对稳定。低地形梯度区间内，前期、后期变化型图谱分布优势明显。反复变化型随地形位指数增加呈现急速下降的趋势，人类活动主要集中在地势平坦，地形起伏小的1级区，优势度明显。持续变化型随地形位增加呈减小趋势，尤其在1～2级地形梯度间呈现急速下降状态，优势分布区间在1～2级区。

图4 各类型图谱在地形梯度上的分布指数Fig.4 The distribution index of each type of map on the topographic gradient

2.3 景观格局的地形梯度效应

借助Fragstats4.4软件完成对阜平全县景观格局指数的测算，由图5可知，时间序列上，景观格局指数最终都保持稳定不变的状态。斑块密度(PD)随时间推移保持小幅上升的趋势,最大斑块指数(LPI)随时间推移保持小幅下降状态,景观形状指数(LSI)随时间推移保持小幅上升的趋势，蔓延度(CONTAG)随时间推移保持下降状态，边缘密度(ED)随时间推移保持上升趋势，香农多样性指数(SHDI)总体保持相对稳定且呈微幅上涨的趋势。

空间序列上，PD值在1～7级地形位区段呈现下降状态，7～8级地形位区段保持相对稳定，8～10级地形位区段呈现上升趋势。PD值在地形梯度内变化跨度较大，景观格局破碎化现象集中发生在低、高等级地形位区域，这些地区景观破碎度较高，低等级地形区间主要是受到人类活动的强烈干扰所致，高等级地形区间主要受复杂地貌的影响。LPI值在1～2级地形区急剧下降，在2～7级地形区保持相对稳定，在7～10级地形区呈现增长的趋势，表明优势景观主要分布在自然质量较好的低等级地形区和人类活动破坏少、保持原始自然风貌的高等级地形区。LSI值在1～4级地形区保持上升趋势，在4～10级地形区呈现下降状态，峰值出现在第4级地形区上，表明景观形状指数高的区域分布在海拔低、坡度较小的低和中低等级地形区。中高和高等级地形区的景观形状指数较小，表明在该地区分布形状较为规则，集中连片。CONTAG值随地形位梯度的升高呈现下降-上升-相对稳定-上升-下降的趋势，表明随地形梯度的升高，景观格局的完整性越高。ED值在1～2级地形区处于上升状态，2～10级地形区处于下降状态，且2～5级地形区下降速度明显大于5～10级地形区，随时间推移，各时点ED峰值均出现在1级地形区上，反映出低和中低地形梯度区域的景观类型易受到人类活动的影响，被分割的程度高，布局较为分散。SHDI值的高低反映其景观异质性，其值越高，表明各景观类型在景观水平上分布越均衡。SHDI值随着地形位指数的增加呈现下降-上升-下降的趋势，景观类型表现更为多样化，1～2级地形区下降缓慢，2～5级地形区急速下降，5～7级地形区缓慢上升，7～10级地形区稳步下降。表明景观格局受地形因子和人类活动的双重影响。

图5 景观格局指数在地形梯度上的变化图Fig.5 Changes of landscape pattern index on terrain gradient

2.4 热环境分布的地形梯度效应

以夏季7月为例，运用反演手段，得到阜平县2009年7月、2014年7月、2019年7月全域地表温度分布状况(图6)，并与地形位指数分级区间进行空间叠置，得到地表温度在不同等级地形区间上的空间分布情况。剔除3期地表温度反演结果的0.01%的最高和最低异常值，经空间分析计算，得到各级地形区地表温度的平均值(表5)。时间尺度上，夏季1～7级地形区的平均地表温度随时间的推移均呈现先上升后下降的趋势；8～10级地形区的平均地表温度呈现上升态势。生产活动中产生的气体排放量，在人类活动较为频繁的低、中低、中高等级地形，2014年后产生的热环境效应较2014年前逐渐减弱。空间尺度上，2009年夏季地表温度平均值随地形梯度的提升呈现下降-上升-下降的走势，临界区间出现在2级地形区；2014年、2019年夏季地表温度平均值均随地形梯度的提升呈现上升-下降的走势，临界区间也出现在2级地形区。平均地表温度的分布状况在一定程度反映了地形因子和人类活动对地表温度的影响程度。

图6 2009-2019年7月阜平县全域地表温度分布图Fig.6 The global surface temperature distribution from July 2009 to July 2019 in Fuping County

表5 不同梯度地形平均地表温度变化

2.5 景观格局与地表温度相关性的地形梯度效应

运用SPSS 26软件计算研究期内各梯度地形内地表温度的平均值和景观格局指数的相关性系数。结果(表6)显示，斑块密度(PD)和最大斑块指数(LPI)与地表温度的相关性较低。景观形状指数(LSI)与地表温度的相关性随时间推移呈现先上升后下降的走势，2009年和2019年在0.05的置信水平上呈现显著正相关性，2014年在0.01的置信水平上呈现极显著正相关性。表明景观类型不规则性越高，热环境效应就越显著。蔓延度(CONTAG)与地表温度的相关性保持上升状态，2009年和2014年在0.05的置信水平上呈现显著负相关性，2019年在0.01的置信水平上呈现极显著负相关性，反映景观格局的完整度越高，其对应的地表温度越低。边缘密度(ED)与地表温度的相关性呈现随时间推移先上升后稳定的趋势，在0.01的置信水平上呈现极显著正相关性，表明景观格局布局越分散，其对应的地表温度值越高。香农多样性指数(SHDI)与地表温度的相关性随时间推移不断加强，2009、2014、2019年均在0.01的置信水平上呈现极显著正相关性。表明景观组合类型越丰富，整体上其对应的地表温度越高。

表6 地表温度与景观格局指数的相关系数

3 结论与讨论

本研究以阜平县2009年7月、2014年7月 、2019年7月3期遥感、DEM影像为基础，借助地形分析、信息图谱分析、景观格局指数和地表温度反演等方法，探究地形因子对阜平县景观格局与热环境效应时空演变的影响程度，对其进行综合分析。结果表明：(1)地形因子对各景观类型优势分布区间产生较大的影响，呈现出明显的时空分异特征。研究区内耕地、草地、建设用地、园地、水体主要集中分布在自然环境条件良好的1～3级或1～2级的低地形梯度区域，为人类生存提供生活空间以及从事各项生产建设活动的生产空间。其他土地优势区分布在4～6级中低地形梯度区域，林地优势区分布在6～10级中高地形梯度区域，林地对各种地形具有较强的适应性，林地为优势景观类型，分布范围逐渐向高地形位指数区域拓展，且地形位指数越高区域的分布指数越大。(2)景观类型的演变上，稳定型占绝对主导分布，其主要分布于中低、中高和高等级地形梯度区域。其余演变类型全部分布于低等级地形梯度区域。稳定型变化图谱中主要类型以其他土地-其他土地-其他土地以及林地-林地-林地为主。(3)景观格局水平上，低、高等级地形位指数区域的景观格局破碎化现象显著，优势景观也主要分布其中。随地形梯度的升高，景观格局的完整性越高。低和中低地形梯度区域的景观类型受到人类活动的影响，被分割的程度高，布局较为分散。中高和高等级地形区域的景观类型分布形状较为规则，集中连片。在低、中低地形梯度内景观类型表现更为多样化。(4)地表温度在时间序列上呈现先上升后下降的趋势，在空间序列上随地形位指数的升高呈现阶梯状分布，且随地形位指数的升高，地表温度呈现下降趋势。地表温度与景观格局存在一定空间相关性，在空间分布上地表温度与PD、LPI相关性较低。地表温度2009年和2019年与LSI在0.05的置信水平上呈现显著正相关性，2014年在0.01的置信水平上呈现极显著正相关性；地表温度2009年和2014年与CONTAG在0.05的置信水平上呈现显著负相关性，2019年在0.01的置信水平上呈现极显著负相关性。地表温度与ED在2009、2014、2019年在0.01的置信水平上均呈现极显著正相关性。香农多样性指数(SHDI)与地表温度2009、2014、2019年均在0.01的置信水平上呈现极显著正相关性。

景观格局的形成和演变受到人类活动和地形因子的共同作用[28-30]，而热环境效应也深受景观格局变化的影响。阜平县地处太行山区，土地利用方式受地形因子的制约，2020年2月刚刚“脱贫摘帽”，因此未来要夯实脱贫成果，进一步扩大发展空间和加大发展力度。要利用好地形梯度跨度大这一地理特性，遵循景观生态学原理和可持续发展原则，因地制宜制定发展战略。土地利用结构的形成是一个长期的过程，转变并不能依靠一朝一夕，在掌握生态发展要义的前提下，要深入挖掘土地潜力，有针对性地开发、配置土地资源，摸清资源禀赋，在保障资源环境承载力和稳定生态的前提下，充分利用好当地的“红色资源”和“绿色资源”培育特色产业，进一步带动农业、林果业和旅游业的发展，进而实现乡村振兴。