2. 西北农林科技大学 杨陵 712100;
3. 黄河水利科学研究院 郑州 450003;
4. 西北大学 西安 710127;
5. 西安科技大学 西安 710054
2. Northwest A & F University, Yangling 712100, China;
3. Yellow River Institute of Hydraulic Research, Zhengzhou 450003, China;
4. Northwest University, Xi'an 710127, China;
5. Xi'an University of Science and Technology, Xi'an 710054, China
黄土高原地区位于我国中北部,其范围包括青海、甘肃、宁夏、山西、内蒙古、陕西和河南7个省(区)的全部或部分地区,总面积约64万km2,占全国陆地总面积的6.67%。黄土高原是世界上水土流失最严重的地区,是我国生态问题集中体现的地区,也是我国水土保持和生态建设的重点地区。
黄土高原地区水土流失影响着黄河的生态安全,长期制约区域经济社会可持续发展。民国时期,李仪祉、张含英等人先后主持的黄委会就十分重视黄河上中游的水土保持,学者们也开展了水土保持研究。建国以后,国家对水土保持工作愈发重视,将其纳入了国民经济建设的轨道。20世纪80年代后,全国的水土保持工作开始了新的发展历程,先后开展了小流域治理、七大江河水土保持重点工程、全国“八大片”治理(黄土高原有无定河、皇甫川、三川河等)、“三北”防护林建设等[1]。后来开展的退耕还林(草)工程、淤地坝建设和坡耕地整治等一系列大型生态工程,对土壤侵蚀控制、生态建设等均起到了良好作用。与此同时随着社会经济发展和城镇化水平提高,人们对环境的干扰也在逐渐减轻。这些都导致了地区生态和环境的明显变化[2-4],如植被覆盖度增加、黄河输沙量减少、居民生活得到明显改善。
为全面掌握黄土高原地区生态环境状况的动态变化,2014年开始,中科院组织开展了黄土高原生态工程生态成效综合评估工作。本项评估工作以黄土高原地区相关野外站的不同尺度监测及生态系统关键过程长期研究为基础,以遥感解译成果数据(表 1)、社会经济年鉴、水文气象观测数据等为基础,以一组指标测度体系(包括治理措施、土地利用、植被覆盖、土壤水分、土壤侵蚀、径流泥沙、粮食和收入等)为对象,利用多种模型(如土壤侵蚀模型RUSLE/CSLE、植物生长模型DSSAT、土壤水分模型模拟和植被承载力模型SWCCV、土壤有机碳模型RothC等)和统计分析方法(包括DPSIR概念框架、TOPSIS熵权法、向量自回归模型和结构方程等),在地块/农户、小流域/典型生态工程/典型样区、行政区/侵蚀和地貌区等多种尺度上,对土地利用和植被覆盖变化、水土保持与生态建设进展、治理过程中河川径流、土壤侵蚀强度和社会经济变化等进行综合分析与评估,形成黄土高原生态工程成效评估报告以及工程后期建设方向与可持续经营管理建议。
通过对黄河上中游地区186个县级行政单元2000-2012年历年统计的水土保持措施数据整理和核查,得到该区逐年的水土保持措施总面积(表 2),目前潼关以上地区保存各类水土保持措施面积约1 474万hm2,其中梯田面积为349万hm2、人工林面积为747万hm2、人工草地面积为258万hm2、坝地18万hm2、封禁治理面积为102万hm2。
2000年以来,黄河上中游地区梯田面积逐年增加,尤其是2000-2006年增加较快,2007年以后逐年增幅趋于稳定;林地面积自2000年以来逐年增加,至2010年增至最大值748.79万hm2,之后林地面积趋于相对稳定;2000年以来,草地面积逐年增加,其中2000-2007年增长缓慢,之后快速增长至2011年,随后增幅放缓;坝地面积呈先增大后减小的趋势,于2007年增至最大值21.44万hm2,后坝地面积逐年减小至2011年,后趋于相对稳定;封禁治理面积自2000年起增长较快,于2005年达到最大值104.07万hm2,后趋于相对稳定。
2.2 土地利用动态变化对20世纪70-80年代、90年代、2000年和2010年4个时段土地利用结构变化进行分析表明,近40年来黄土高原土地利用变化特征如图 1、图 2所示。
![]() |
图 1 黄土高原土地利用结构变化 |
![]() |
图 2 黄土高原土地利用图(1980-2010年) |
(1)耕地先增加后减少。由20世纪70-80年代的17.96万km2增加到90年代的19.67万km2、2000年达20.64万km2,到2010年减少至19.98万km2,占全区面积的比例分别为28.7%、31.47%、33.00%和31.83%。
(2)林地持续增加。由2 0世纪7 0-8 0年代的7.38万km2增加到90年代的7.72万km2、2000年增至9.25万km2,2010年增至12.35万km2,占全区面积的比例分别为11.80%、12.36%、14.80%和19.68%。
(3)草地呈减少趋势。由20世纪70-80年代的31.15万km2减少至90年代的30.22万km2、2000年减至25.99万km2,2010年减至24.61万km2,占全区面积的比例分别为49.70%、48.34%、41.56%和39.21%。
(4)居民地增加,未利用地减少。黄土高原前期耕地增加的原因是“以粮为纲”的政策,后期减少的主要原因是退耕还林(草)工程中对陡坡耕地的退耕;林地和居民地的增加,除了退耕还林(草)工程外,还与持续多年的水土保持生态建设、社会经济发展后人们对自然的干扰减少有关(图 1、图 2)。
黄土高原中部是水土流失治理重点区,也是土地利用变化最显著的地区。该区域典型样区(面积约400 km2)统计表明(图 3),耕地先略为增加,然后减少,如吴起样区耕地由36.70%增至41.44%,然后减少到14.43%;林地增加,如吴起样区林地由1.33%增至21.05%;草地的变化有增加(如绥德)、稳定(如吴起)和减少(如延安、离石)3种情形。草地的增加或稳定(吴起和绥德),与本地区耕地急剧减少有关;草地的减少(延安和离石),与当地耕地变化幅度相对较小、林地明显增加有关。
![]() |
图 3 黄土高原典型样区土地利用变化 |
自20世纪80年代初开始,30多年来植被覆盖逐年明显增加,黄土高原植被覆盖度变化在时间和空间上有所差异,这与该地区气候变化和人为活动有着密切联系。从空间格局特征看,中部丘陵区(丘I,丘II)和土石山区明显增加,NDVI(植被覆盖度指数)均值由0.21增加到0.48,净增加128.6%(图 4)。从时间看,1982-2013年,NDVI均值由0.30上升到0.45,净增加50.0%,夏天变化尤为明显,由0.36增至0.59,增幅63.9%(图 4)。从时间看,NDVI增加主要发生在春、夏两季,占全年总增量60%以上(图 5)。
![]() |
图 4 黄土高原植被覆盖度(NDVI)变化时空图谱(1982-2013年) |
![]() |
图 5 黄土高原植被覆盖度(NDVI)各年度逐月变化 |
自2000年以来,黄土高原植被变化的区域特征为:(1)按省区看,陕西增幅最大,为0.47%;青海增幅最小,为0.16%;山西、宁夏、甘肃、内蒙古分别增长0.38%、0.33%、0.31%、0.23%。(2)从黄河各支流看,延河流域增幅最大,为1.14%;无定河次之,为0.59%;渭河干流增幅最小,仅为0.43%。(3)从黄土高原侵蚀类型区看,水力侵蚀区植被覆盖增长速率较快,风力侵蚀区则较缓。所有这些变化首先与1999年以来退耕还林(草)政策的实施有关,其次也与社会经济发展后的自然恢复和人工干扰减少有关。
2.4 土壤水文生态动态变化黄土高原植被恢复引起土壤干燥化现象受到普遍关注,目前黄土高原土壤干层分布广泛,表现出明显的空间变异和独特的分布格局(图 6)。土壤干层的平均厚度为160 cm,在剖面上起始形成深度平均为270 cm。区域尺度上土壤干层厚度的变异程度属于强变异(CV=110%),土地利用对土壤干层具有极显著的影响。地统计结果表明黄土高原地区土壤干层厚度具有强烈的空间依赖性,而干层起始形成深度呈中等空间依赖性。干层厚度和起始形成深度的空间变程分别为34 km和125 km,空间异质比分别为79%和50%。
在黄土高原西部(即宁夏盐池-陕西定边-宁夏固原-甘肃静宁-甘谷一线以西)和中部地区特别是陕西和山西交界的沿黄河地区,土壤干层厚度较厚;而在黄土高原沿黄灌区(如宁夏、内蒙灌区)、内陆灌区、汾河灌区、南部关中平原等地较薄[5]。黄土高原地区的林地普遍存在下伏土壤干层(125个样地中有102个样地存在土壤干层),干层发育较草地和农地严重(干层起始形成深度为140 cm,干层厚度为304 cm,干层内土壤含水量为7.9%)。
确定了黄土高原中部沿降水梯度带的最适植被盖度指标(即:叶面积指数LAI)。根据黄土高原降水量从东南向西北递减,土壤质地由细变粗,植被类型由乔木林、乔灌混交林向灌木林和灌木草原过渡的分布特点,利用土壤-植被-大气水量转化模型,确定各降水梯度带适宜的植被叶面积指数[6]:人工乔木林适生区(降水量大于550 mm)最适宜叶面积指数范围是2.5 < LAI < 3.5;人工灌木适生区(降水量250-350 mm)最适宜叶面积指数范围是0.8 < LAI < 1.5;人工乔冠过渡区域(降水量350-500 mm)最适宜叶面积指数范围是1.5 < LAI < 2.5;降水量小于250 mm时最适宜叶面积指数范围是LAI < 0.8。
2.5 土壤侵蚀动态变化黄土高原土壤侵蚀以水蚀为主,2000-2010年11年的平均土壤侵蚀强度为15.2 t/(hm2·a)。中度以下侵蚀强度( < 25 t/(hm2·a))占到整个黄土高原的80.5%,主要分布在河谷平原区、山西东部黄土丘陵区、土石山区及地势比较平坦的风蚀地区(图 7)。土壤侵蚀强度比较大的区域呈东北-西南向分布于黄土高原中部的黄土丘陵沟壑区和黄土高塬沟壑区。
2000-2010年黄土高原土壤侵蚀强度整体上呈显著下降的趋势,并具有明显的区域性差异(图 8)。土壤侵蚀强度变化最大的区域分布在黄土丘陵沟壑区和黄土高塬沟壑区,该地区2000和2010年平均侵蚀强度约32.0 t/(hm2·a),其中坡耕地侵蚀强度为38.7 t/(hm2·a),局部地段土壤侵蚀依然较强[7]。遇到强降水年,土壤侵蚀强度依然较大。如延河流域2013年降水量767.21 mm(是多年平均降水量的140%多),土壤侵蚀模数和输沙模数分别为50.45 t/(hm2·a)和41.23 t/(hm2·a)[7],发生强烈土壤侵蚀的危险依然存在。2000-2010年中度以上侵蚀区土壤侵蚀强度以1-3 t/(hm2·a)速度在减少,其主要分布在黄土高原重点水土流失区,包括黄土丘陵沟壑区和高塬沟壑区的陕西榆林、延安地区和山西吕梁、临汾地区等沿吕梁山脉一带和黄土高塬沟壑区的甘肃庆阳、定西、白银和宁夏固原、平凉等六盘山一带。侵蚀强度以大于1 t/(hm2·a)速度增加的区域零星分布在山西太原、阳泉,内蒙古鄂尔多斯东部沿黄土高原边界一带和青海海南地区[8]。
2.6 河流径流和输沙量动态变化自20世纪70年代以来,黄河流域潼关以上干支流水沙量持续减少,2000年以后减少更多,尤其是2008年以来急剧下降。1919-1959年潼关水文站实测年均径流量426.1亿m3、年均输沙量15.92亿t。20世纪60年代以后,黄河水沙总体呈减少态势,尤以80年代中期以后减少显著。20世纪70年代、80年代、90年代水土保持措施减少泥沙的量分别为3.05亿t、4.29亿t和4.21亿t。2000-2012年潼关水文站年均输沙量减少至2.76亿t,较1919-1959年平均值减少82.7%;同期年均径流量减至231.2亿m3,减少了45.8%。
与1919-1959年基准期相比,2000-2012年黄河流域潼关以上年均减少径流量217.7亿m3,其中降雨、下垫面变化、能源开发等对减流量的贡献分别为2.1%、76.7%和21.2%,在下垫面变化各因素中,梯田、林草地(包括封禁)、坝地、经济社会用水、地表水、地下水超采等对减流量的贡献分别为5.0%、9.8%、3.3%、37.8%、11.7%和9.1%。2000-2012年黄河流域潼关以上年均减沙量为13.7亿t,其中降雨、下垫面变化、河道采砂等对减沙量的贡献分别为22.1%、71.5%和6.4%,在下垫面变化各因素中,梯田、林草地、坝地、水库河道等对减沙量的贡献分别为11. 9 %、2 3. 0 %、22.0%和14.6%。水土保持措施对减沙起到了主导作用,水利措施及经济社会用水对减水起到了主导作用,降雨等自然因素对水沙变化的影响已处于相对次要的地位。自20世纪80年代以来,水利水土保持综合治理等人类活动对水沙变化的影响程度逐渐增强,而降雨的作用则相对减弱。例如,与20世纪70年代相比,80年代水利水土保持综合治理等人类活动对径流量、输沙量减少的贡献率均分别在30%、40%以上;90年代人类活动对减水的作用升至40%以上,减沙作用升至65%以上;2000年以来人类活动对水沙变化的作用明显进一步增强,坡面及沟道工程对控制黄河水沙起了主导作用。可以说,近年来坡面和沟道的生物和工程等多种措施共同作用,使常年降水条件下黄河输沙量被控制到了人类活动影响之前的程度。但应当认识到,随着坝库等工程措施拦沙能力的逐渐下降,维持一个可持续的植被生态系统对有效保持土壤和控制黄河输沙量反弹具有更加重要的作用。
2.7 经济结构动态变化黄土高原地区1995-2010年间社会经济发展情况处在连续向好的态势,生态修复工程带来的社会经济效益显著。从经济结构来看,1995-2010年间该区第二产业和第三产业产值增长态势明显,而第一产业表现为先下降后增长(图 9)。其中,第一产业产值1996年最高为588.52亿元,2001年达到最低的441.91亿元,经历了2000-2004年的低谷期之后开始回升,由2004年的482.75亿元上升至2010年的579.69亿元;1995-2010年间,该区第二产业和第三产业产值逐年增加,第二产业由1995年的1 254.43亿元增至2010年的4 670.19亿元,年均增长率达18.15%,第三产业也由1995年的852.25亿元增至2010年的3 108.40亿元,年均增长率为17.65%。
![]() |
图 9 黄土高原地区三次产业产值变化动态 |
从三次产业占比来看,1995-2010年间黄土高原地区第一产业所占比重逐年下降,第二产业逐年上升,第三产业比重先上升后平稳。其中,1995年第一、二、三产业比重分别为21.40%、46.80%和21.80%;2001年后第一产业占比下降速度逐渐趋于缓慢;第二产业比重自2002年后不断增加,增加幅度最大;第三产业在经历了1995-2000年一段时间的快速上升之后增速渐缓;2010年第一、二、三产业比重分别为6.94%、55.87%和37.19%。2000-2003年是黄土高原三种产业结构转变的关键节点,而这一时期也是该区生态修复工程实施的初始期。前期的产值下降可能是该区实施的退耕还林(草)政策导致了耕地面积的下降,而后期的增长则可能是由于技术进步带来的单位面积产值与农产品附加值的增加。
3 水土保持与生态工程建设的对策与建议(1)加大水土流失防治力度,推进以植被生态功能提升为核心的生态工程。在黄土高原维持一个可持续的植被生态系统对有效保持土壤和控制黄河输沙量反弹具有更加重要的作用。遵循地带性规律和因地制宜原则的前提下,持续增加林草植被覆盖率和提高林草措施质量,提升人工植被稳定性,强化水土保持功能;继续实施退耕还林(草)工程,加强水土保持综合治理工程,继续推进以坡改梯为重点的土地整治和小流域综合治理;加大植被保护力度,全面实行生态公益林补偿机制,实施天然林保护,强化封山育林,充分发挥黄土高原地区的自然修复能力。
(2)实施生态工程跟踪监测,为水土保持与生态建设科学决策提供依据。建立健全各级水土流失监测机构,加强落实监测队伍建设,配备监测设施设备,提高监测工作的组织保障能力;建立并完善水土保持与生态建设工程效益监测评价体系,对生态工程实施进展及成效做出客观的监测,科学评估生态工程实施的生态环境效应、社会经济效益、生态服务功能价值以及生态经济可持续性,为工程建设与科学部署后续生态工程提供基础数据。
(3)强化水土保持与生态建设科学研究,发展与丰富现代水土保持学科。在新的社会经济背景下,要进一步明确黄土高原地区土壤保持、粮食生产、草畜业的战略定位和科学布局,辨析区域土壤侵蚀的区域分异特征以及水沙关系的变异规律,确定维持区域水土保持功能的植被结构特征阈值,发展生态系统水土资源承载力定量评价和管理模型,以推进黄土高原地区退化生态系统生态服务功能的可持续修复。
(4)优化产业结构,重视农户转移就业,减轻生态环境承载压力。对生活条件恶劣,地质灾害频繁的地区,实施生态移民,有计划地开展实施异地扶贫搬迁,落实国家精准扶贫的目标;对水土流失地区的富余劳动力,开展专业性技能培训,为农户非农自主经营提供小额贷款、完善农村社会保障体系,提高农民素质与就业能力,为农户创造更多的非农就业机会,排除农户难以向非农产业转移的障碍,使更多农村剩余劳动力逐渐从传统小农生产中转移出来,并达到更高的收入水平。
(5)完善政策机制,强化依法防治,严格执行实施水土保持防治法律法规。建立稳定的投入机制,加大对水土流失防治的资金投入;完善金融扶持和税收优惠等政策,引导企业、个人等社会资金积极投入水土流失防治事业;尽快研究制订《生态补偿条例》,完善水土保持防治法律体系;加大普法教育力度,增强广大群众的生态保护意识和法制观念;强化退耕还林(草)补助、落实水土保持生态补偿、完善生态工程管护,夯实巩固生态治理成果;严格执行《水土保持法》《森林法》《草原法》等法律法规,加大对破坏生态行为的执法力度;统筹处理生态修复与社会经济发展的关系,建立生态修复政策实施长效机制,在确保生态效益的前提下,最大程度发挥经济效益,兼顾长远的社会效益。
致谢 本研究中,土地利用分析采用了中科院地理科学与资源所、遥感与数字地球所提供的部分数据,特此致谢。[1] | 唐克丽, 主编. 中国水土保持. 北京: 科学出版社, 2004. |
[2] | 李敏. 水土保持对黄河输沙量的影响. 中国水土保持科学, 2014, 12(6): 23–29. |
[3] | 刘晓燕, 杨胜天, 李晓宇, 等. 黄河主要来沙区林草植被变化及对产流产沙的影响机制. 中国科学:技术科学, 2015, 45(10): 1052–1059. |
[4] | 姚文艺, 冉大川, 陈江南. 黄河流域近期水沙变化及其趋势预测. 水科学进展, 2013, 24(5): 607–616. |
[5] | Wang Y Q, Shao M A, Zhu Y J, et al. Impacts of land use and plant characteristics on dried soil layers in different climatic regions on the Loess Plateau of China. Agricultural and Forest Meteorology, 2011, 151(4): 437–448. DOI:10.1016/j.agrformet.2010.11.016 |
[6] | Fu W, Huang M B, Gallichand J, et al. Optimization of plant coverage in relation to water balance in the Loess Plateau of China. Geoderma, 2012, s 173-174: 134–144. |
[7] | 杨勤科, 等. 基于GIS的黄土丘陵区年侵蚀产沙经验模型研究. 西安: 西北大学城市与环境学院, 2015. |
[8] | Sun W, Shao Q, Liu J. Soil erosion and its response to the changes of precipitation and vegetation cover on the Loess Plateau. Journal of Geographical Sciences, 2013, 23(6): 1091–1106. DOI:10.1007/s11442-013-1065-z |