伊洛河(33° 34'— 35° 4'N, 111° 7'— 112° 58'E)主要流经河南省境内, 是黄河十大支流之一, 由伊河和洛河组成, 全长约447 km, 流域面积约1.86万k m2[15](图1)。 该流域地处暖温带向北亚热带的过渡区域, 年平均气温约为14.6 ℃, 年平均降水约为600 mm。 伊洛河流域海拔57~2 655 m, 呈西高东低特征, 生态与环境过渡性明显, 其上游为丘陵山区, 以森林、 草地为主, 而下游主要为黄土冲积平原, 是重要的粮食生产基地。 受其地形地貌特征和人类活动影响, 伊洛河流域的生态环境复杂且脆弱。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 研究区域示意图Fig.1 Schematic diagram of the study area

研究数据主要包括遥感影像数据、 基础地理数据、 气象数据和社会经济数据。 其中, 遥感影像数据涉及1990年— 2020年期间伊洛河流域的Landsat系列图像, 数据来源于地理空间数据云(http://www.gscloud.cn/)。 基础地理数据包括伊洛河流域边界, 数据来源于中国科学院资源环境数据中心(https://www.resdc.cn/)。 气象数据、 社会经济数据包括人口, GDP, 第一、 二、 三产业的产值等数据, 主要来源于伊洛河流域县区人民政府官网与历年统计年鉴。

(1)土地利用转移矩阵

土地利用转移矩阵来源于系统分析中对系统状态与状态转移的定量描述^[16]。 通常的土地利用转移矩阵中, 行表示T₁时点土地利用类型, 列表示T₂时点土地利用类型。 P_MN表示土地利用类型N转化为土地利用类型M的面积。

P=P11P12...P1NP21P22...P2N............PN1PN2...PNN(1)

(2) 土地利用变化量模型

土地利用/土地覆被数量变化模型参照现有研究成果^[17], 选取土地利用/土地覆被净变化量指数、 土地利用/土地覆被交换变化量指数、 土地利用/土地覆被总变化量指数系统分析研究区土地利用类型在数量上的动态变化。

Dj=|P+j-Pj+|(2)

Sj=2min(Pj+-Pjj, P+j-Pjj)(3)

Cj=Dj+Sj(4)

式中, D_j为特定时期内某区域土地利用类型净变化量百分比; P_+j为特定时期内期末第j类土地利用类型面积百分比; P_j+为特定时期内期初第i类土地利用类型面积百分比; S_j为特定时期内第j类土地利用/土地覆被类型的交换量百分比; P_jj为特定时期内第i类土地利用类型未发生变化的面积百分比; C_j为特定时期内第j类土地利用/土地覆被类型的总变化量指数。

(3)土地利用变化程度

选取土地利用变化面积比重D来评估研究区域内土地利用变化的剧烈程度, 土地利用变化面积比重D指各土地变化面积之和占区域总面积的比例, 见式(5)^[11]

D=AS×100%(5)

式(5)中: A为区域内各土地变化面积之和; S为研究区面积。

(4)弦图可视化分析

弦图(Chord diagram)是使用圆形坐标系展现多组数据间相互关系的可视化形式^[18]。 为呈现研究区域内各种土地利用类型的转变过程, 利用R语言中的Circlize package将土地利用转移矩阵在弦图中可视化, 构建土地利用类型变化的轨迹模型。

(5)冗余分析

冗余分析(RDA)是响应变量矩阵与解释变量之间多元多重线性回归的拟合值矩阵的PCA分析, 是基于对应分析发展而来的一种排序方法^[17]。 此分析常用于研究群落与环境因子之间关系, 本工作将不同的土地利用类型视为群落, 采用Canoco 5软件对土地利用类型和环境因子进行RDA分析, 对RDA分析后的数据采用土地利用类型、 样本与环境因子变量的三序图呈现。

图2给出了伊洛河流域1990年— 2020年不同土地利用类型的时空分布及各类土地利用类型占比。 由图可知, 1990年— 2020年伊洛河流域的土地利用类型主要为林地、 草地、 耕地、 建设用地、 水域和未利用土地。 在空间上, 林地和草地主要分布在伊洛河上游地区和洛河中游南岸, 整体分布呈现聚集状态。 耕地主要分布在伊洛河中下游两岸地区。 建设用地主要聚集在下游靠近水域的平坦地带, 在中上游平坦地带也有少量聚集。 水域主要为伊洛河及其支流和水库。 未利用土地占比较少, 主要分布在建设用地周围。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 伊洛河流域1990年— 2020年土地利用类型变化图Fig.2 Changes in land use types in the Yiluo River Basin from 1990 to 2020

在数量上, 林地占主导地位, 占比为48.4%~70.7%, 整体呈现先减小后增加的趋势; 其次为耕地, 占比为20.8%~49.1%, 整体呈现先增加后减小的趋势; 而建设用地所占比例较小, 为1.8%~7.0%, 整体呈现快速增长的趋势; 水域所占比例为0.4%~1.1%, 整体呈现减少趋势; 草地和未利用土地所占比例很低, 分别为0.0%~1.6%和0.0%~0.6%。

2.2.1 土地利用变化量

图3给出了伊洛河流域每10年间各种土地利用类型的净变化量指数D_j、 交换变化量指数S_j以及总变化量指数C_j, 3个指数从变化量的角度表现了各种土地利用类型的变化情况。 1990年— 2020年耕地、 林地的总变化量最大, 各贡献了该流域土地变化量的40.45%和41.51%以上, 占绝对的主导地位。 建设用地次之, 其他类型土地变化量较少。 1990年— 2020年各时段林地和耕地总变化量不断增加, 林地从14.52%上升到19.72%, 耕地从14.15%上升到19.73%; 建设用地在1990年— 2020年各时段变化量处于2.03%~4.81%。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 伊洛河流域每10年间各种土地利用类型的变化量指数Fig.3 The change index of various land use types in the Yiluo River Basin every 10 years

2.2.2 土地利用迁移变化

表1给出了1990年— 2020年的土地转移矩阵。 由表可知, 各土地利用类型面积的相互转移变化的总体表现为林地的转入转出面积最大, 耕地次之, 再则是建筑用地, 水域、 草地和未利用土地转入转出面积最小。 林地1990年— 2020年间向其他土地利用类型转变量为876.8 km², 由其他土地利用类型向林地转变量为4 009.9 km², 整体面积增长较大。 减少的林地主要转变为耕地, 增加的林地主要来源于耕地, 部分来源于草地、 建设用地。 耕地1990年— 2020年向其他土地类型转变量为4 859.1 km², 由其他土地利用类型向耕地转变量为663.3 km², 整体面积减少较多, 减少的耕地主要转化为林地和水域, 还有部分转化成建设用地, 耕地面积增加主要由林地、 建设用地转化而来。 建设用地1990年— 2020年向其他土地类型转变量为166.0 km², 由其他土地利用类型向建设用地转变量为1 226.9 km², 整体面积增长较大, 减少的建设用地主要被耕地占用, 还有部分转化为林地和水域。 建设用地面积的增加来源于林地转化以及部分草地的开垦, 还有一部分是来源于耕地。 水域、 草地和未利用土地面积较小, 在此不做讨论。

表1

Table 1

表1(Table 1)

表1 1990年— 2020年土地转移矩阵(km2) Table 1 1990— 2020 land change index (km2) 1990年 总计 林地 耕地 建设用地 水域 草地 未利用土地 2020年 林地 10 272.6 3 755 115.7 56.3 70.2 12.8 14 282.6 耕地 581.1 3 564.6 38.9 36.3 3 4 4 227.9 建设用地 154.9 995.8 200.1 62.4 0.6 13.3 1 427.1 水域 1.8 24.2 6.1 56.2 0 1.5 89.8 草地 118.3 2.8 0 0 5.6 0 126.7 未利用土地 20.7 81.4 5.4 5.6 0.2 3.4 116.7 总计 1 1149.4 8 423.8 366.2 216.8 79.6 35

表1 1990年— 2020年土地转移矩阵(km2) Table 1 1990— 2020 land change index (km2)

图4给出1990年— 2020年每10年伊洛河流域土地利用类型数量的变化轨迹。 从土地利用类型转变的轨迹来看, 1990年— 2000年间土地利用类型之间的相互转化以林地的流出为主, 主要流向为耕地; 2000年— 2010年以耕地的流出为主, 主要流向为林地; 2010年— 2020年依然以耕地的流出为主, 主要流向为林地。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 1990年— 2020年伊洛河流域土地利用变化轨迹Fig.4 Land use change trajectory in Yiluo River Basin from 1990 to 2020

1990年— 2000年、 2000年— 2010年、 2010年— 2020年的土地利用变化面积比重分别为7.11%, 22.34%和23.61%, 土地利用变化总面积分别为3 467.66, 4 528.35与4 785.76 km²。 3个10年的时间段里, 土地利用变化的总面积在不断增加, 说明土地利用变化剧烈程度随时间呈上升趋势, 在研究时间段内2010年— 2020年达到最高。

图5给出了伊洛河流域林地、 耕地以及建设用地的重心变化, 从不同土地利用类型重心在空间上的分布来看, 林地重心始终处于伊洛河流域的上游区域, 伊洛河上游为丘陵地区, 人类活动较少适合树木的生长。 耕地始终处于伊洛河流域的中游区域, 该区域是冲积平原, 海拔和坡度较低, 适宜发展耕作, 而建设用地的重心则始终围绕洛阳市主城区。 1990年— 2020年间林地重心整体是向东北方向偏移, 结合图2可以看出由于洛河北岸退耕还林效果显著, 林地由山区向平原扩张, 并且2010年— 2020年间迁移幅度最大; 耕地的重心总体向南迁移, 由于林地的扩张导致伊洛河流域平原地带西部和东北部耕地减少。 耕地重心呈“ z” 型迁移特征, 尤其以2000年— 2010年的迁移幅度最大; 建设用地的重心在1990年— 2020年间先向东北迁移, 后向西南迁移, 表明洛阳市的城市发展先以洛阳城区发展为主, 1990年— 2000年的迁移幅度最大, 在2000年后, 上游各城镇的社会经济发展逐渐活跃。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 伊洛河流域土地利用类型变化轨迹图Fig.5 Change track map of land use types in the Yiluo River Basin

表2给出伊洛河流域的社会经济和气候数据。 由表2可知, 在社会经济方面, 伊洛河流域的人口在1990年— 2020年间整体呈增长趋势, 近年来趋向稳定。 而GDP在30年间发展迅速, 呈现出指数式增长。 1990年— 2020年间整体看第二产业产值占比最大, 其中第三产业增速最快, 2020年已经超过第二产业产值, 表明伊洛河流域整体经济由中低收入水平向中上收入水平转化, 反映了工业化和城市化发展的程度较高。 伊洛河流域的气温在过去30年间缓慢上升, 由年平均14.0 ℃上升到15.4 ℃, 年平均降水量变动较大, 但均值处在500~800 mm左右, 没有明显的增加或减少趋势。

表2

Table 2

表2(Table 2)

表2 伊洛河流域社会经济及气候参数 Table 2 Socio-economic and climatic parameters of the Yiluo River Basin 年份 GDP /亿元 人口 /万人 第一产业 /亿元 第二产业 /亿元 第三产业 /亿元 降水 /mm 气温 /℃ 退耕还林政策 实施情况 1990 73.0 563.6 12.6 39.7 20.7 754.4 14.0 未实施 1995 245.1 596.4 22.6 143.4 79.1 609.6 14.5 未实施 2000 422.8 626.0 39.3 230.8 152.6 794.0 14.0 未实施 2005 1 112.4 653.4 110.5 648.4 353.5 531.7 14.2 已实施 2010 2 320.3 703.5 187.7 1 396.3 736.5 527.3 14.4 已实施 2015 3 508.8 728.5 236.4 1740.7 1 531.7 652.7 14.8 已实施 2020 5 128.4 720.0 254.1 2 312.2 2 562.1 606.9 15.4 已实施

表2 伊洛河流域社会经济及气候参数 Table 2 Socio-economic and climatic parameters of the Yiluo River Basin

通过RDA分析绘制了样方、 土地利用类型和环境变量关系的RDA三序图, 如图6所示。 图6中轴1、 轴2的特征值分别为0.951和0.048。 人口、 第一、 二、 三产业生产总值和建设用地呈现显著正相关; 退耕还林政策与林地呈现显著正相关, 与耕地呈现显著负相关, 表明社会经济发展与退耕还林政策实施对土地利用类型的转变有显著影响。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 RDA三序图Fig.6 RDA three sequence diagram

由RDA三序图还可看出, 降水与土地利用类型变化相关性很低, 而气温与林地变化则呈正相关。 伊洛河流域气温在过去30年间呈现缓慢上升的趋势, 有研究^{[3, 10]}表示气温对植被覆盖有重要影响, 这些影响可能会引起土地利用类型的变化, 主要体现在植被面积的变化。 所以, 气温对林地变化起到了一定影响。

自中国改革开放以来, 人口规模扩张和社会经济发展迅速。 几乎所有城市都呈现出极快的扩张趋势, 建设用地面积越来越多^[19]。 姜朋辉^[20]等研究表明社会经济因素是导致黑河中游地区土地利用变化的根本原因。 麦丽开· 艾麦提^[21]等研究表明人口增长是叶尔羌河中下游流域耕地、 建设用地面积增加的直接原因。 同时, 城市的发展也对生活用水和生产用水产生了更高的需求, 导致水域面积减少^[20]。 伊洛河流域人口增长和经济发展直接导致了建设用地需求日益增长, 使得建设用地的面积持续扩大, 并且向上中游平坦区域扩张。

除了社会经济因素外, 政策因素往往在短期内就能引起土地利用格局重大调整^{[22, 23]}。 大量研究表明, 农业土地利用政策的实施能够很大程度上影响区域农业生产方式和发展路径, 进而影响土地利用格局的变化^[24]。 随着社会生产力的提高, 政府对生态建设更加重视, 2002年1月10日, 国务院确定全面启动退耕还林工程。 退耕还林政策的实施, 直接导致了林地面积的增长以及耕地面积的减少。 政策和经济发展的综合作用, 使得经济效益较差的耕地逐渐转为林地, 在伊洛河流域的直接体现是中部的大量耕地转为林地。

黄河流域发展要求统筹推进山水林田湖草沙综合治理、 系统治理、 源头治理, 着力保障黄河长治久安, 着力改善黄河流域生态环境, 着力优化水资源配置, 着力促进全流域高质量发展。 伊洛河流域30年来生态向好, 林地增加是生态环境改善的重要指标, 但是由于生活需水和生产生活需水增加, 伊洛河流域水域面积在减少, 未来要研究城镇化/经济发展对水资源的影响, 未来土地利用规划要充分考虑水资源保护和调配。