1978—2017 年西藏耕地面积变化态势及驱动力分析

郭健斌， 曾维莲， 袁庆娟， 刘天平

（西藏农牧学院 西藏农村经济发展研究中心，西藏 林芝 860000）

土地是人类赖以生存和发展最重要的物质基础，耕地是各类土地利用类型中最重要的一部分。随着社会经济的不断发展和一定区域内人口的增加，耕地生产压力不断增强，耕地面积也随之变化。由于耕地不仅具有生产功能，还具有生态服务功能和景观文化功能[1]，因此影响耕地面积变化的因子错综复杂，归纳起来有社会经济和自然因素两方面，短期内人类的社会经济活动对耕地面积变化起决定性作用[2]。学者们从不同尺度对耕地变化进行了研究，如Goldewijk 等认为：近3 个世纪以来，全球耕地面积增加了5.5 倍，新增加的耕地主要由草地和林地转变而来[3]。张丽娟等利用全球1982—2011 年土地利用/覆被的矢量数据，分析了20 世纪80 年代以来全球耕地变化的基本特征及空间格局，认为自20 世纪80 年代以来，全球耕地面积增加了528.768 万km2，北美洲、南美洲、大洋洲耕地面积呈显著增加趋势[4]。赵晓丽等研究表明：1987—2000 年中国耕地总面积略有增加，2000—2010 年中国耕地总面积逐步减少，近30 年间，耕地总量相对稳定[5]。曹银贵等分析全国1996—2005 年耕地变化的情况，认为全国耕地面积持续减少，主要是生态退耕所致[2]。更多学者关注我国省域和市县域尺度的耕地面积变化及驱动力研究，发现农业科技进步[6]、社会经济发展的强度和规模[7]、城市化水平[8]、居民储蓄存款[9]等是影响耕地面积变化的主要因素。

也有学者对西藏耕地面积变化做了研究，如杨春艳等以遥感资料为数据源，对1980—2010 年西藏耕地面积时空变化特征进行研究，结果表明近30 年西藏耕地面积总体呈增长趋势[10]。王建林等采用1981—2001 年的统计资料，研究了近20 年来西藏耕地与人口相关性，发现西藏耕地面积在波动增加，但增加速度缓慢，西藏耕地面积与乡村人口数之间呈正相关关系[11]。本文采用1978—2017 年的长时间序列数据，综合分析影响西藏耕地面积变化的社会经济因素，对于认识西藏耕地数量变化的主导因素具有重要意义。

一、 数据来源及研究方法

（一） 研究区域概况

西藏位于青藏高原西南部，地处北纬26°50′～36°53′，东经78°25′～99°06′之间，平均海拔在4 000 m 以上，素有“世界屋脊”之称。全区面积122.84 万km2，约占全国总面积的八分之一。全区辖6 个地级市、1 个地区，8 个市辖区，66 个县。2017 年末总人口337.15 万人，其中城镇人口104.15 万人，农村人口233 万人。2017 年西藏全区生产总值1 477.63 亿元，比2016 年增长9.1%。

西藏高原地貌大致可分为喜马拉雅山区、藏南谷地、藏北高原和藏东南高山峡谷区，耕地集中分布在藏南河谷及河谷盆地中，即“一江两河”（雅鲁藏布江及其支流拉萨河、年楚河流域）地区（图1）。这一地区是西藏高原最主要的农业区，享有西藏粮仓的美誉。

（二） 数据来源

文中使用的西藏耕地资源空间分布数据来源于国家科技基础条件平台——国家地球系统科学数据共享平台（http://www.geodata.cn）；耕地面积、人口数量、地区生产总值、粮食作物播种面积等数据全部来源于《西藏统计年鉴》（2018—1993），其中1978—1992 年数据来源于《西藏统计年鉴》（1993—1996）。

（三） 研究方法

利用SPSS17.0 对影响耕地面积变化的因子做主成分分析，利用降维的思想，把原来的多指标转化为少数几个主成分，使复杂问题简单化，从而获取更加有效的数据信息，进而计算1978—2017 年西藏耕地面积变化驱动力；并采用回归分析法，分析各驱动因素对于耕地变化的影响程度。

二、 结果与分析

（一） 西藏耕地面积变化分析

由图2 可知，西藏耕地面积呈现波动增加趋势，1978—2017 年西藏耕地面积从227 600 hm2增加到242 750 hm2，年增加耕地面积为388.46 hm2。西藏和平解放以后耕地面积持续增加，至1979 年耕地总面积为229 900 hm2；1979—1987年波动减少至221 400 hm2；1987—1999 年缓慢增加至231 160 hm2；1999—2004 又持续减少至222 740 hm2；2004—2010 年 波 动 增 加 至229 530 hm2；2010—2017 年快速增加到历史最高值。

西藏人均耕地面积呈现出持续减少的趋势，1978 年西藏人均耕地面积为0.13 hm2，2017 年已减少为0.07 hm2。

（二） 西藏耕地面积变化驱动因子主成分分析

1. 指标选取与评价

根据主成分分析法的要求[12]，结合诸多学者研究成果，选择西藏1978—2017 年数据作为分析样本，将耕地面积作为因变量Y，将年末总人口（X1）、人口城镇化率（X2）、全社会固定资产投资（X3）、财政收入（X4）、农村居民人均纯收入（X5）、粮食作物（指小麦、青稞为主的谷物、薯类和豆类）单位面积产量（X6）、有效灌溉面积（X7）、每公顷耕地平均化肥施用量（X8）、地区生产总值（X9）、第二产业占比（X10）、第三产业占比（X11）、农业总产值（X12）等12 项涉及社会经济发展的因素作为自变量。

自变量相关系数矩阵显示（表1），影响耕地面积变化的12 个自变量中存在着不同程度的相关性，相关系数大于0.90 的因子共有22 对，表明各自变量包含信息存在重叠，运用主成分分析法是必要的。对自变量进行KMO 检验和Bartlett 球形检验显示：KMO 统计量达到0.881，Bartlett’ s 球形检验显著性水平小于0.01，适合进行主成分分析。

表 1 变量相关系数矩阵

2. 主成分分析

主成分分析结果显示，2 个主成分解释的累计贡献率已经达到94.33%，第3 主成分特征值为0.329（＜1.000）（表2），说明该主成分的解释力度还不如直接引入原变量大，因此，提取F1与F2这2 个主成分就基本能够解释原来所有变量所包含的大部分信息。

由表3 可知，第1 主成分在全社会固定资产投资、财政收入、农村居民人均纯收入、地区生产总值指标上的负荷较大，其特征值5.912，贡献率49.270%。第2 主成分在年末总人口、人口城镇化率、粮食作物单位面积产量、每公顷耕地平均化肥施用量、第三产业占比指标上的负荷较大，其特征值5.407，贡献率45.060%。可将西藏耕地面积变化驱动力因子归纳为经济因子和人口、科技、产业因子。

表 2 主成分特征值和贡献率

表 3 主成分载荷矩阵

1978—2017 年，西藏自治区地区总产值由6.65 亿元增长到1 477.63 亿元，增长了近200 多倍；固定资产投资由1.85 亿元增长到2 252 亿元，增长了1 200 多倍。随着经济的不断增长，人民生活逐渐富裕，对农副产品的需求逐渐增强，间接促进耕地面积增加。人口是耕地面积变化最根本和较活跃的驱动因子，不断增加的人口需要有更多的耕地提供粮食以满足其生存需要。西藏总人口从1978 年的178.82 万人增加到2017 年的337.15 万人，几乎翻了一番。科技进步使耕地生产力不断提升，单位面积粮食产量增加，一定程度上缓解了耕地压力。近40 年来，西藏农药化肥的施用量、农业机械的普及率和农田灌溉条件都有大幅度的增加，粮食作物单位面积产量由1987年的2 505 kg/hm2增加到2017 年的5 711 kg/hm2，极大地缓减了耕地压力。随着社会经济不断发展，西藏自治区第三产业所占比率由1978 年的22%增加到2017 年的51%，产业结构的变化使得大量农业人口从第一产业调整到第三产业，从事农业生产的劳动力减少亦缓减了耕地压力。

（三） 西藏耕地面积变化驱动力分析

根据上述研究结果，以每个主成分所对应的特征值占所提取主成分总的特征值之和的比例作为权重计算主成分综合评价模型：

式（1）中，F 为西藏耕地面积变化驱动力得分；λi为第i 主成分的特征值占总的特征值之和的比例（i=1，2），Fi为第i 主成分。

对影响西藏耕地面积变化的2 个主成分进行加权汇总，得出西藏1978—2017 年耕地面积变化驱动力，进而绘制出西藏耕地面积变化驱动力得分图。如图3 所示，1978—2017 年西藏耕地面积变化驱动力从-0.9 增加至1.7，整体呈现持续上升的趋势，表明影响西藏耕地面积变化的驱动力逐渐加强。

（四） 西藏耕地变化驱动因子回归分析

回归分析是一种揭示因变量与自变量之间变化规律的分析方法，将西藏1978—2017 年的耕地面积做标准化处理，将其作为回归分析的因变量(Y)，把驱动力因子F1、F2作为回归分析的自变量，得到如下的回归方程：

将主成分F1、F2的12 个自变量因素代入式（2）中，并计算得出原始自变量在原回归模型中的相应参数（表4），从而得到消除了多重线性的回归模型。

消除多重共线性的标准回归模型能够比较全面地反映西藏耕地面积变化情况，其系数大小反映了耕地面积变化对所选驱动因子的敏感程度。从式（3）可以看出，西藏耕地面积数量变化与年末总人口（X1）、人口城镇化率（X2）、全社会固定资产投资（X3）、财政收入（X4）、农村居民人均纯收入（X5）、粮食作物单位面积产量（X6）、有效灌溉面积（X7）、每公顷耕地平均化肥施用量（X8）、地区生产总值（X9）、第二产业占比（X10）、第三产业占比（X11）、农业总产值（X12）的相关系数分别为0.026、0.040、0.192、 0.182、 0.162、 0.044、 0.025、 0.032、0.160、0.093、-0.083 和0.096。可见，西藏耕地面积变化驱动因子中全社会固定资产投资、财政收入、农村居民人均纯收入和地区生产总值等是影响耕地面积变化的主要因素，并呈现显著正相关关系，即经济快速发展对耕地面积变化产生直接且明显的影响。

表 4 主成分的回归分析结果

三、 结论

（1）1978 年—2017 年，西藏耕地面积从227 600 hm2增加到242 750 hm2，呈现波动增加的趋势，年增加耕地面积为388.46 hm2；人均耕地面积持续减少，1978 年西藏人均耕地面积为0.13 hm2，2017 年已减少为0.07 hm2。

（2）运用主成分分析方法对西藏耕地面积变化的社会经济因素进行综合量化分析，得出影响西藏耕地面积变化的两个主成分，根据各因素在主成分中的负荷大小，将西藏耕地面积变化驱动力因子归纳为经济因子和人口、科技、产业因子。

（3）通过建立耕地面积变化驱动力模型，计算得知1978—2017 年西藏耕地面积变化驱动力从-0.9 增加至1.7，整体呈现持续上升的趋势，表明影响西藏耕地面积变化的驱动力逐渐加强。

（4）利用回归分析法对影响西藏耕地面积变化的因素进行分析，发现全社会固定资产投资、财政收入、农村居民人均纯收入、地区生产总值等是影响耕地面积变化的主要因素，并呈现显著正相关关系，即经济快速发展对耕地面积变化产生直接且明显的影响。