郭健斌，张 英，张志伟，侯 磊，曾维莲
(1.西藏农牧学院 资源与环境学院，西藏 林芝 860000；
2.西藏农牧学院 西藏乡村振兴研究院，西藏 林芝 860000；
3.西藏农牧学院 植物科学学院，西藏 林芝 860000)
土地利用变化研究是近年来全球变化研究的热点[1]，主要包括土地覆被变化和土地利用程度变化2 个方面。但当前学术界普遍关注土地覆被变化的空间形态、类型与数量、演化机理和驱动力研究[2-4]，而对土地利用程度的关注度不高。土地利用程度是一个地区土地开发数量和程度的度量指标，反映了土地自然属性和人类活动对土地利用的综合影响[5]。探究区域土地利用程度的变化，有利于人类可持续利用土地资源。国外关于土地利用程度的研究侧重于农用地的投入产出分析[6]；
国内有学者选用1 个[5]或多个指标[7]对土地利用程度进行量化研究，为了便于区域之间的横向比较，多数学者采用[8]或改进[9]楚玉山等[10]提出的土地利用程度综合指数量化分析方法对中国不同自然地理区域[11-12]和行政区域[8,13]的土地利用程度进行综合评价。随着遥感和地理信息系统技术的快速发展，多源高精度影像数据在土地利用程度变化的研究中也得到普遍运用[13-15]。空间研究中经常使用的方法是空间自相关分析[16-19]。由于土地利用变化是一类具有高度空间异质性的地理现象，利用空间自相关技术分析其空间聚集规律的同时，还可以将空间数据可视化。空间自相关性使用全域和局域2 种指标，Moran’sⅠ指数在判断区域土地利用程度的空间聚集上具有可靠性[20]和不易受偏离正态分布影响的特点[19]，得到广泛应用[11-12,21]。
在中国，土地利用程度的研究多聚焦于黄土高原[9,14,22]、西南喀斯特山区[12]和西北内陆干旱区[23-24]等生态脆弱地区以及中东部人口密集地区[8,11,21]。青藏高原作为中国典型的生态环境脆弱地区和国家重要生态安全屏障，目前仅见三江源地区土地利用程度变化的时空分析[25]和《西藏自治区土地利用》对西藏各县土地利用程度的对比分析[10]。林芝市位于西藏东南部的高山峡谷地区，垂直落差大，土质疏松，局部暴雨多，滑坡和泥石流等自然灾害频繁。随着人口和人类社会经济活动的增加，不但会改变林芝市的社会经济发展模式，也会威胁脆弱的生态环境。因此，客观和准确地认识林芝市土地利用强度及其空间分异对高山峡谷地区的可持续发展具有重要意义。本研究以藏东南林芝市为研究区域，以2010 年和2020 年遥感影像为基础，建立土地利用强度及其变化模型，采用Moran’sⅠ指数对林芝市土地利用程度进行空间异质性分析，探究区域土地利用程度的内部差异和空间相关性，可为当地国土空间规划和生态环境保护决策提供理论基础。
1.1 研究区概况
林芝市是西藏自治区下辖地级市，位于E92°09′～98°47′，N26°52′～30°40′ (图1)。东与昌都市和云南省迪庆藏族自治州毗邻，西与拉萨市和山南市交界，北与那曲市相连，南与缅甸和印度接壤。东西长646.7 km，南北宽353.2 km，平均海拔3 100 m，总面积11.7 万km2(包括印度非法控制区)[26]。林芝市地处青藏高原南缘，地势自西北向东南倾斜，雅鲁藏布江和怒江等水系切割高原，形成典型的高山峡谷地貌，海拔落差大。市域气候分为热带北缘和亚热带山地湿润气候区以及高原温带半湿润气候区，立体气候明显[27]。
图1 研究区位置Fig.1 Location of the study area
林芝史称贡布，是藏民族发祥地之一。和平解放后，1960—1963 年曾设林芝专区，随后撤销，1986 年林芝地区行政公署恢复成立，2015 年3 月国务院批复同意撤销林芝地区，设立林芝市。2020 年全市下辖1 区6 县，包括54 个乡(镇)、2 个街道办，共504 个村(居)，总人口23.89 万，是以藏族为主体，汉、门巴和珞巴等35 个民族和僜人多民族聚居区。境内分布有中国最大的原始林区，耕地散布于河谷地带，草地位于林线上下，垂直地带性明显。2020 年全市实现地区生产总值191.34 亿元，可比增长7.9%，三次产业的比例为6.4∶38.7∶54.9[28]。
1.2 数据来源
土地利用数据来源于中国科学院计算机网络信息中心创立的地学大数据平台——地理空间数据云(https://www.gscloud.cn/)的2010 年Landsat-7 ETM+和2020 年Landsat-8OLI 2 期各11 景卫星影像(空间分辨率为30 m×30 m)，为确保云量小于10%，影像日期选择当年10—12 月或次年1—3 月。利用ENVI 5.1 对2 期遥感数据进行辐射定标、大气校正和裁剪等预处理，对于部分ETM+数据条带缺失的问题，利用ENVI 软件中的条带修复补丁，并基于单文件三角插值的条带填充进行修复。结合研究区实际情况，依据《土地利用现状分类标准》(GB/T 21010—2017)，选取典型林地、水域、耕地、建设用地、草地和未利用地为训练样本，对研究区典型地物进行最大似然法监督分类。同时借助Google Earth 高分辨率影像和部分实地考察数据，每一地类随机选取100 个样本进行精度评价，2010 年和2020 年2 期土地利用类型识别精度分别为90%和92%，kappa系数分别为0.88 和0.90，满足研究需求。
1.3 研究方法
1.3.1 土地利用程度指数
土地利用程度指数是以土地开发利用程度来反映人类活动对某一区域土地利用变化的影响水平[14]。本研究参考庄大方等[29]数量化土地利用程度分析方法，将土地利用程度按照土地自然综合体在社会因素影响下的自然平衡保持状态分为4 级，并分级赋予指数，从而给出土地利用程度的定量表达。其计算公式为
式中：LA指某一区域土地利用综合指数，为100～400 的连续函数，值越大表明该区域土地利用程度越高；
Ai为 第i级土地利用程度的分级指数，参照已有研究[8,10-11,29]中的分级标准并结合研究区实际，将未利用地赋值为1，林地、草地和水域赋值为2，耕地赋值为3，建设用地赋值为4；
Ci为第i级土地利用程度分级面积百分比，即第i级土地利用程度分级类型与区域总面积的比值；
n为土地利用程度分级数，本研究取值为4。
由于全局空间自相关大小对采样方式和采样尺度较敏感[30]，介于研究区域林草地和未利用地占比大，垂直地带性明显，综合考虑评价单元的实效性和工作量，最终采用10 km×10 km 的网格对研究区2 期土地利用类型图进行采样，共计获得1 314 个采样网格进行土地利用程度指数计算及土地利用程度变化分析；
再利用ArcGIS 10.6，在自然断点基础上结合当地实际，依据土地利用程度指数将林芝市土地利用程度划分为强、较强、中等、较弱和弱5 级(表1)。
表1 土地利用程度分级表Tab.1 Grades of land use intensity classification
1.3.2 土地利用程度变化
本研究利用2 个时期土地利用程度指数的变化来衡量某一区域土地利用程度的变化趋势[23]，其表达式为
式中：CLAb-a为 该区域研究期初到期末土地利用程度变化指数；
LAb和 LAa分别为该区域研究期末和期初土地利用程度指数。
1.3.3 空间自相关分析
空间自相关是指某一变量在一定区域内空间上的相关性，土地利用程度(变化)作为连续空间的度量指标也具有空间相关性。本研究选用全域和局域Moran’sⅠ指数衡量土地利用程度(变化)的空间聚集程度。
全域Moran’sⅠ统计量用于描述土地利用程度(变化)在整个研究区域的空间依赖程度，可反映土地利用程度(变化)是否在空间上聚类，其表达式为
式中：n为样本数，Xi和Xj分 别为区域内第i和第j处的土地利用程度(变化)；
X为区域内所有样本的平均值；
Wij为空间权重矩阵，本研究采用基于距离标准的二进制空间权重矩阵，当Wij=1时，区 域i与j相 邻，当Wij=0 时，区 域i与j不相邻。I值的显著性用以下公式检验[31]：
式中：E(I) 和v(I)分别为Moran’sⅠ的数学期望和方差。
全域Moran’sⅠ指数范围在[-1,1]，在既定的显著性水平下，I＞0 时，表示土地利用程度(变化)之间存在显著正相关，值越大相关性越强，空间上聚合；
I＜0 时，表示土地利用程度(变化)之间存在显著负相关，绝对值越大相关性越强，空间上离散；
I=0 时，表示土地利用程度(变化)不存在空间自相关，随机分布。
虽然全域空间自相关可以反映土地利用程度是否在空间集聚，但它不能进一步显示具体的空间集聚特征[19]。局域空间自相关则可以进一步展示土地利用程度空间集聚或离散的位置和范围，本研究选用的局部空间关联指标(local indicators spatial association，LISA)为局域Moran’sⅠ指数。
局域Moran’sⅠ指数用于描述土地利用程度在局域空间是否存在空间自相关，能更直观地反映研究区土地利用程度的空间集聚状况，其表达式为
式中变量的含义与全域Moran’sⅠ统计量公式相同。局域Moran’sⅠ指数表示土地利用程度(变化)在不同空间位置上的集聚模式，当I＞0 时，表示以高—高值或低—低值的模式集聚；
当I＜0时，表示以高—低值或低—高值的模式集聚。
2.1 土地利用结构及利用程度的时空特征
2.1.1 土地利用结构变化
2010—2020 年，林芝市的草地、林地、未利用地和水域4 类用地稍有减少，其中，林地的占比从2010 年的53.49%降至2020 年的53.43%，未利用地的占比从2010 年的16.34%降至2020年的16.32%，草地的占比从26.20%降至26.17%，水域的占比从2.86%降至2.85%。2010—2020年，林芝市的耕地占比从2010 年的1.10%增加至2020 年的1.17%；
建设用地明显增加，从2010的25.05 km2增加至2020 年的66.57 km2，增长了2.66 倍，占比从0.02%增至0.06% (表2)。
表2 2010—2020 年林芝市各用地类型面积及其占比Tab.2 Area and proportion of land use types in Nyingchi City from 2010 to 2020
2.1.2 土地利用程度时间变化
由表3 可知：2010—2020 年，弱等级网格单元面积的占比保持在0.51%；
较弱等级网格单元面积的占比从6.50%增加至6.55%；
中等等级网格单元面积的占比从26.02%降至25.82%；
较强等级网格单元面积的占比从 48.48% 降至 45.61%；
强等级网格单元面积的占比从 18.48% 增至 21.51%。整体而言，土地利用程度呈现出由中等和较强等级向强等级转变的态势。中等和较强等级土地利用程度网格单元面积从85 486.70 km2降至81 959.90 km2；
除弱和较弱等级土地利用程度网格单元面积增加1.07 和58.09 km2外，其余3 467.60 km2(占比98.33%)全部转变为强等级土地利用程度。
表3 2010—2020 年林芝市土地利用程度等级面积及其占比Tab.3 Area and proportion of land use intensity in Nyingchi City from 2010 to 2020
2.1.3 土地利用程度空间特征
林芝市土地利用程度空间分布(图2)表明：较强和强等级土地利用程度网格单元主要分布在墨脱县、察隅县南部的低海拔区域和人类活动密集的雅鲁藏布江及其支流尼洋河宽河谷流域，弱和较弱等级土地利用程度网格单元主要分布在高海拔的市域东北部、西北部和中部东喜马拉雅南迦巴瓦区域。2010 年林芝市土地利用程度指数平均值为181.27，2020 年为181.37。2010—2020 年，强等级土地利用程度网格单元在尼洋河流域和墨脱县有明显增加。
图2 林芝市土地利用程度等级Fig.2 Grades of land use intensity in Nyingchi City
2.2 土地利用程度指数空间自相关分析
2.2.1 全域空间自相关
由图3 可知：2010 年和2020 年的Moran’sⅠ指数均为正数，数值分别为0.588 和0.587，大于0、接近1；
z值分别为31.21 和26.32，均大于临界值(2.58)，并且通过了 α=0.01 的显著性检验，表明林芝市土地利用程度存在较强的正空间相关性，即土地利用程度空间分布上不是随机的，而是表现为高值或低值的集聚。可见，林芝市土地利用程度的空间联系特征为：土地利用程度高的网格单元相对趋向于与土地利用程度高的网格单元相邻，土地利用程度低的网格单元相对趋向于与土地利用程度低的网格单元相邻。2020年Moran’sⅠ指数较2010 年仅减少了0.001，说明2010—2020 年土地利用程度的空间分布呈现集聚的趋势基本稳定。
图3 2010 年(左)和2020 年(右)林芝市土地利用程度散点图Fig.3 Scatterplots of land use intensity in Nyingchi City in 2010 (left) and 2020 (right)
2.2.2 局域空间自相关
由图4 可知：低—低值集聚区主要分布在市域东北部和西北部高海拔地区；
高—高值集聚区分布在市域低海拔的南部和东南部部分地区、尼洋河流域和雅鲁藏布江宽河谷地区。2020 年尼洋河流域高—高值集聚区较2010 年明显增加。
图4 2010 年(左)和2020(右)年林芝市土地利用程度LISA 集聚图Fig.4 LISA agglomeration graphs of land use intensity in Nyingchi City in 2010 (left) and 2020 (right)
2.3 土地利用程度变化空间自相关分析
2010—2020 年，林芝市土地利用程度变化全域Moran’sⅠ指数为0.223，z值为22.21，大于临界值，通过了 α=0.01 的显著性检验，表明其土地利用程度变化量存在较强的空间正相关性，即土地利用程度变化剧烈的网格单元相对趋向于与土地利用程度变化剧烈的网格单元相邻，土地利用程度变化微弱的网格单元相对趋向于与土地利用程度变化微弱的网格单元相邻(图5)。
图5 2010—2020 年林芝市土地利用程度变化散点图Fig.5 Scatterplots of land use intensity change in Nyingchi City from 2010 to 2020
由局域空间自相关LISA 集聚地图(图6)可知：2010—2020 年林芝市土地利用程度变化高—高值集聚区集中分布在尼洋河流域的巴宜区的八一镇、百巴镇和更张门巴民族乡以及工布江达县的巴河镇、错高乡和工布江达镇，在巴宜区的鲁朗镇、波密县的扎木镇、墨脱县的墨脱镇和格当乡以及察隅县的下察隅镇和上察隅镇也有分布；
低—高值集聚区分布在上述高—高值集聚范围的周边；
低—低值集聚区分布在墨脱县的墨脱镇和背崩乡、米林县的丹娘乡和羌纳乡以及巴宜区的米瑞乡、林芝镇和布久乡，在高海拔区域如工布江达县朱拉乡和察隅县古玉乡等也有分布；
高—低值在波密县松宗镇和市域边缘有少量分布。
图6 2010—2020 年林芝市土地利用程度变化LISA 集聚图Fig.6 LISA agglomeration graph of land use intensity change in Nyingchi City from 2010 to 2020
2010—2020 年，林芝市人口从19.51 万增加到23.89 万，固定资产投资从65.17 亿元增加到125.12 亿元[28]，人口和资金集中在海拔低、地势平坦的河谷地带，导致这一区域建设用地和耕地增加，草地、未利用地和林地减少。建设用地的增长反映了林芝市非农经济活动的快速增长；
耕地的增减变化说明林芝市人地关系渐趋紧张，城市规模扩张态势明显。在全球变暖背景下，青藏高原内流区得益于冰川融雪和降水增加，湖泊水域扩张明显[32-33]；
而位于藏东南外流区的林芝市冰川消融[33]，水域面积减少。
杨玉婷等[23]研究表明张掖市甘州区2005 年的土地利用程度指数为255.25，而林芝市2010 年和2020 年的土地利用程度指数分别为181.27 和181.37，远低于西北干旱内陆河流域，土地利用程度低。林芝市土地利用程度空间分布明显受地形条件和社会经济发展的影响，高海拔地区难利用土地占比高，土地利用程度低；
低海拔区域林草地占比高，土地利用程度较高；
平坦的河谷地带耕地占比高，村庄和城镇分布广泛，土地利用程度最高。土地利用程度综合指数能够诠释林芝市土地利用程度的空间分布与自然环境间的关系，但尚不能准确厘定其作用机制，需要进一步研究。2010—2020 年，林芝市土地利用程度由中等和较强等级向强等级转变，说明随着当地人口数量和社会经济活动的不断增加，除人类难以涉足的区域外，林芝市土地利用程度都在不断增强。
林芝市域西北部和东北部为高海拔区域，是土地利用程度低—低聚集区域；
而南部和东南部为低海拔区域，是土地利用程度高—高集聚区域；
以巴宜市区为中心的尼洋河和雅鲁藏布江宽河谷区域是土地利用程度高—高集聚区域。这些LISA 集聚特征与实际情况较吻合，能够准确地反映林芝市自然地理环境条件和社会经济发展的特点。2010—2020 年，土地利用程度变化(增加)在低海拔、宽河谷和人类活动密集地区呈高—高集聚的特征，其外围呈低—高集聚的过渡特征，说明随着林芝市城镇化水平的不断提高，农村地区土地利用程度在不断地缓释，而城镇用地利用程度不断增强，且有持续向外扩张的趋势。在可利用土地资源极其稀缺的藏东南高山峡谷地区，未来土地利用应该以挖掘现有建设用地潜力，提高城镇用地的集约性，引导其有序扩张。
本研究利用Landsat 影像解译土地利用数据，计算土地利用程度及其变化，采用空间自相关法分析土地利用程度及其变化的空间特性，结果表明：林芝市人类社会经济活动密集的低海拔河谷平地土地利用程度高；
人类社会经济活动稀少的高海拔区域土地利用程度低。2010—2020年该区域土地利用程度不断增强，低海拔河谷平地的强等级土地利用程度网格单元增加更为明显。市政府驻地尼洋河流域是土地利用程度变化的高—高值集聚区，即土地利用程度变化的热点区域，今后国土空间规划和生态环境保护应着重关注这一区域的土地集约利用。
——林芝地理教学(2014年7期)2014-01-13幸福的程度环球时报(2012-03-21)2012-03-21
