黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应.pdf
《黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应.pdf》由会员分享,可在线阅读,更多相关《黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应.pdf(11页珍藏版)》请在咨信网上搜索。
1、黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应汪滨,张志强(1.北京林业大学水土保持学院,水土保持与荒漠化防治教育部重点实验室,北京100083;2.山西吉县森林生态系统国家野外科学观测研究站,临汾042200)摘要:黄土高原土壤侵蚀严重,为此中国从 1999 年起实施了大规模的退耕还林工程。为了分析退耕还林土地利用变化对土壤侵蚀的影响,该研究以黄土高原清水河流域为研究区域,将 20002020 年流域退耕还林工程的实施依据主要措施的不同划分为 4 个阶段,应用 RUSLE(reviseduniversalsoillossequation)模型分析土壤侵蚀强度的变化特征,采用情景模拟方法
2、提出一个区分土地利用变化和降雨变化对土壤侵蚀影响程度的算法,判别土地利用变化对土壤侵蚀的影响程度,将土地利用变化分解为土地利用转换和改造 2 种形式,在剔除降雨变化影响的基础上分析土地利用变化对土壤侵蚀的影响过程。结果表明:1)2000、2005、2011、2014 和 2020 年流域平均侵蚀模数分别为 36.21、41.02、24.93、23.72 和 8.24t/(hm2a),土壤侵蚀强度明显下降;土地利用变化和降雨变化对土壤侵蚀的阶段平均影响程度分别为75.23%和 24.77%,土地利用变化在流域土壤侵蚀的变化中起了主导作用。2)流域土地利用转换区侵蚀强度的变化直接受转换过程中地类类
3、别变更及所实施主要相关措施差异的影响,改造区侵蚀强度的变化直接受改造过程中所实施主要相关措施的影响。剔除降雨变化的影响后:改造区的阶段平均起始侵蚀模数较转换区高 43.47%,其水土流失综合治理的难度总体上大于转换区;转换区侵蚀模数的阶段平均下降量较改造区高 50.80%,改造区侵蚀量的阶段合计减少量占流域阶段合计减少量的 71.16%,土地利用转换在降低其实施地区土壤侵蚀强度方面发挥了重要作用,而土地利用改造因实施面积较大在减少流域土壤侵蚀总量方面发挥了重要作用。3)剔除降雨变化的影响后,草地侵蚀量的阶段合计变化量占流域阶段合计变化量的 70.51%,且草地阶段合计变化量中改造区占 67.4
4、1%,其变化特别是其改造对流域土壤侵蚀的影响最大。该研究在分析土地利用变化对土壤侵蚀的影响程度和过程方面作了一些尝试,研究结果可为黄土高原退耕还林成果巩固及高质量发展有效措施的制定提供科学依据。关键词:土壤;侵蚀;土地利用;土地利用改造;土地利用转换;影响程度;退耕还林工程;RUSLE 模型doi:10.11975/j.issn.1002-6819.202302075中图分类号:S157文献标志码:A文章编号:1002-6819(2023)-12-0060-11汪滨,张志强.黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应J.农业工程学报,2023,39(12):60-70.doi:10.
5、11975/j.issn.1002-6819.202302075http:/www.tcsae.orgWANGBin,ZHANGZhiqiang.ResponseofsoilerosiontolandusechangedrivenbyGrainforGreenProjectinatypicalwatershedoftheLoessPlateauJ.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2023,39(12):60-70.(inChinesewithEnglishabstra
6、ct)doi:10.11975/j.issn.1002-6819.202302075http:/www.tcsae.org0引言由于土壤抗蚀能力较弱以及人类长期的不合理土地开发利用,黄土高原属于世界上土壤侵蚀最严重的地区1-2,为了从根本上有效控制土壤侵蚀,中国于 1999年起在黄土高原实施了大规模的退耕还林工程。该工程以改变不合理的土地利用、遏制土壤侵蚀、减少入黄泥沙并兼顾生态、经济和社会效益为主要目的,以坡耕地退耕还林还草、荒山造林、封山育林和梯田建设等为主要措施3-4,实施后土地利用结构明显优化,取得了显著的综合效益尤其是水土保持效益4-6。退耕还林工程历时较长,不同阶段实施了不同的主要
7、措施,直接驱动土地利用持续发生变化,深入研究这一变化对土壤侵蚀的影响,对于促进退耕还林成果巩固和推动退耕还林高质量发展具有重要意义。土壤侵蚀的发生是降雨、土地利用、地形、土壤等各种因素相互影响与制约的综合结果7,其中降雨和土地利用是 2 个主要影响因素8-9。黄土高原有关试验与观测研究表明,降雨量、降雨强度和土地利用都对土壤侵蚀产生了明显影响10-12。运用模型情景模拟(如 RUSLE(reviseduniversalsoillossequation)13、SWAT(soilandwaterassessmenttool)14、MMF(Morgan,MorganandFinney)15)和数理统
8、计(如双累积曲线16、线性回归17、弹性系数18)等方法,黄土高原已开展了降雨和土地利用对退耕还林期间土壤侵蚀变化影响程度方面的研究,得出的共同结论是土地利用变化在一定程度上减轻了土壤侵蚀,但在影响程度方面存在着 2 种完全不同的结论,一些研究观点认为土地利用变化的影响程度相对较大16,18-20,另一些研究的结论则相反13,17,21。退耕还林工程的具体实施对象是相关的土地利用类型(以下简称地类),实施的结果是这些地类中部分地区(或地块)发生了地类的转换(即类别变更),如坡耕地通过退耕还林转至林地、草地通过工程整地转至耕地14,22,与此同时这些地类内未发生转换的地区(或地块)也在发生变收稿
9、日期:2023-02-14修订日期:2023-04-13基金项目:国家科技支撑计划项目(2015BAD07B030303)作者简介:汪滨,博士生,研究方向为流域治理。Email:通信作者:张志强,教授,博士生导师,研究方向为森林水文、土壤侵蚀与流域管理。Email:Z第39卷第12期农 业 工 程 学 报Vol.39No.12602023年6月TransactionsoftheChineseSocietyofAgriculturalEngineeringJune2023化,如坡耕地实施了梯改、林草地实施了封禁,且从已有的研究来看未发生转换的面积较大14,22-23。土地利用是人类对土地自然属性
10、的利用方式和状况24,地类转换地区所发生的变化和未转换地区所发生的变化通过改变侵蚀环境,如耕地转草地后植物由作物转变为草丛、林草地的封禁提高了植被覆盖度,都会对土壤侵蚀产生直接影响,故同时研究这 2 种变化对土壤侵蚀的影响尤其是剔除降雨变化后的影响,可完整地揭示土地利用变化对土壤侵蚀的影响过程。目前黄土高原已开展了大量有关这方面的研究,其中以对地类或区域土壤侵蚀整体变化影响的研究为主25-27,也有一些研究侧重转换对侵蚀的影响28-30,但前者没有区分这 2 种变化,后者没有关注未转换地区的变化,且均叠加了降雨变化的影响,不能直接反映土地利用变化对土壤侵蚀的影响,关于剔除降雨变化后这 2 种变
11、化影响方面的研究少有报道。实施退耕还林工程后,一些地类对土壤侵蚀的影响最大,黄土高原已有的研究表明影响最大的地类主要有旱地和中低覆盖草地31、林地和草地32等,有的地区建设用地也发挥了重要作用33,这些研究主要依据降雨变化下地类整体或其转换地区侵蚀变化的重要性,关于剔除降雨变化后地类转换地区和未转换地区侵蚀变化重要性方面的研究也少有报道。本研究以山西省吉县清水河流域为对象,应用 RUSLE模型分析流域 20002020 年土壤侵蚀强度的变化特征,通过情景模拟判别土地利用变化和降雨变化对土壤侵蚀的影响程度,将土地利用变化分解为土地利用转换和土地利用改造 2 种形式,并在剔除降雨变化影响的基础上,
12、分析退耕还林工程实施过程中这 2 种变化形式对土壤侵蚀的影响及其在影响中的相互关系和各自作用,识别影响流域土壤侵蚀变化的主要地类,以完整地反映土地利用变化对土壤侵蚀的影响程度和过程,为黄土高原退耕还林成果巩固及高质量发展有效措施的制定提供科学依据。1研究流域及其退耕还林工程实施概况清水河发源于山西省吉县高天山西麓,经州川河向西南汇入黄河,是黄河中游河口镇至龙门区间左岸的一级支流。吉县水文站控制的清水河流域地理位置为1103647E110560E,36218N361623N(图1)。#!11055E11050E11045E11040E11035E11055E11050E11045E11040E1
13、1035E3615N3615N3610N3610N365N365N人祖山高天山吉县郭家垛前杨家峪段家堡山头麦原贾家塬柏坡底川庄雨院沟清水河024 km1图例Legend水文站Hydrologic station雨量站Rainfall station!山峰Mountain peak#河流RiverWatershed boundary流域边界115E110E105E100E115E110E105E40N35N黄土高原黄河N图 1清水河流域区位及水文站和雨量站分布图Fig.1LocationoftheQingshuihewatershedanddistributionsofhydrologicand
14、rainfallstations跨吉县、乡宁两县,总土地面积 436km2,其中吉县境内占 94.67%。为黄土残塬沟壑区典型流域,土壤类型以褐土为主。属暖温带半湿润大陆性季风气候,多年平均气温 10.2,多年平均降水量 541.5mm,年内降水不均,根据吉县水文站 20002020 年的观测结果,汛期510 月降雨量占全年降水量的 86.58%。清水河流域从 2000 年起属于黄土高原退耕还林工程的重点实施地区,经历了不同的发展时期:2000 年下半年开始实施,2005 年完成退耕和造林任务;从 2008 年起进入退耕还林成果巩固期,流域所在的吉县在继续实施荒山造林和封山育林的同时,大力开展
15、坡改梯与基本农田建设;2010 年起吉县发挥地方优势,将发展苹果产业作为成果巩固的重要举措;2014 年下半年起开始实施新一轮工程。为了完整地反映流域退耕还林工程的实施过程,本研究以 2000、2005、2010、2014和 2020年 5个标志年份为节点,将流域已进行的退耕还林实施时期(20002020 年)分为20002005、20052010、20102014 和 20142020 年 4 个时段,并依据工程实施期间主要措施的不同,将其对应的实施阶段分别分为第 1 阶段(20002005 年)、第 2 阶段(20062009 年)、第 3阶段(20102014 年)和第 4 阶段(201
16、52020 年)。2数据与方法2.1数据来源与处理2.1.1遥感影像与土地利用分类本研究选取 2000 年 5 月 21 日 Landsat7、2005 年 6月 12 日 Landsat5、2011 年 6 月 29 日 Landsat5、2014 年5 月 20 日 Landsat8 和 2020 年 6 月 5 日 Landsat8 影像分别作为 5 个标志年份的代表影像(因 2009 年 G22 临吉高速公路的开工建设对 2010 年沿线地区土地利用的影响较大,故用 2011 年影像替代,相应的时段分别调整为20052011 和 20112014 年,对应的实施阶段分别调整为第 2 阶
17、段(20062010 年)和第 3 阶段(20112014 年),这些影像均源自中国科学院计算机网络信息中心地理空间数据云(http:/ 8 类。在 ERDAS9.2 中经监督分类并借助 GoogleEarth 和 91卫图的高空间分辨率影像,进行人工目视判读与修正34后得到 5 期土地利用空间分布图。利用野外实地采点数据、地块利用历史实地访问调查资料等对 5 期分类结果进行精度检验,总体解译精度均达到 90%以上。2.1.2降雨数据与侵蚀性降雨标准流域土壤侵蚀主要发生在汛期,流域内及周边共 10个站点的汛期逐日降雨数据来源于中华人民共和国水文年鉴黄河流域水文资料。依据谢云等35基于黄土高原试
18、验小区实测数据的研究结果,将 12mm/d 作为侵蚀性降雨标准。第12期汪滨等:黄土高原典型流域土壤侵蚀对退耕还林土地利用变化的响应612.1.3土壤类型的分布和属性数据流域土壤类型的分布数据来源于山西省土壤图集,属性数据来源于山西土种志。2.1.4DEM 和 NDVI 数据DEM 数据为空间分辨率 30m 的 ASTERGDEM 影像,源自中国科学院计算机网络信息中心地理空间数据云(http:/ 数据中 2000、2005、2011 和 2014 年数据为 510 月共 24 期分辨率 500m 的MODISNDVI 月合成影像,同样源自地理空间数据云(http:/ 2020 年分辨率 50
19、0m 的月合成影像,2020 年数据选取美国航空航天局(NASA)地球观测数据与信息系统(EOSDIS)(https:/search.earthdata.nasa.gov/)4 月 22 日至 10 月 30日共 12 期分辨率 500m 的 MODISNDVI16 日合成影像。2.2研究方法2.2.1RUSLE 模型的应用黄土高原沟蚀问题突出36,其中浅沟侵蚀分布较广37,模拟试验结果表明浅沟侵蚀量占坡面侵蚀总量的26.6%59.7%38。流域属于残塬沟壑区,沟蚀严重,参考黄土高原地区的有关研究39-41,其侵蚀强度的估算采用增加了浅沟侵蚀影响因子 G 的改进 RUSLE 模型42:A=RK
20、LS GCP(1)式中 A 为年平均土壤侵蚀模数,t/(hm2a);R 为降雨侵蚀力因子,MJmm/(hm2ha);K 为土壤可蚀性因子,thm2h/(hm2MJmm);L 和 S 分别为坡长和坡度因子,无量纲;G 为浅沟侵蚀影响因子,无量纲;C 为覆盖管理因子,无量纲;P 为水土保持措施因子,无量纲。本研究根据数据的可获取情况,选用了各因子的适当算法。1)降雨侵蚀力因子 R。降雨侵蚀力反映了降雨引起土壤侵蚀的潜在能力43,采用章文波等44基于日降雨资料提出的算法:Ri=kj=1(Pj)(2)=21.5867.1891(3)=0.8363+18.177/Pd+24.455/Py(4)式中 Ri
21、为第 i 个半月的降雨侵蚀力值,MJmm/(hm2h);Pj为该半月内第 j 日的侵蚀性降雨量,mm;k 为该半月内的侵蚀性降雨日数,d;和 为参数;Pd和 Py分别为研究时段内侵蚀性日降雨量和侵蚀性年降雨量的多年平均值,mm。采用泰森多边形法为各雨量站所控制的区域赋值。2)土壤可蚀性因子 K。土壤可蚀性反映了土壤性质对侵蚀的敏感程度45,采用 EPIC 模型46中的算法以及张科利等47基于中国土壤属性提出的因子值修正算法:KEPIC=0.13170.2+0.3exp0.0256Sa(1Si100)(SiCl+Si)0.310.25CoCo+exp(3.722.95Co)10.7SnSn+ex
22、p(5.51+22.9Sn)(5)K=0.01383+0.51575KEPIC(6)式中 Sa为砂粒质量分数,%;Si为粉粒质量分数,%;Cl为黏粒质量分数,%;Co为有机碳质量分数,%;Sn=1Sa/100;K 为基于中国土壤属性对 KEPIC的修正值。3)坡长因子 L 和坡度因子 S。坡长和坡度是诱发土壤侵蚀的重要地形因素,L 因子和 10以下坡度的 S 因子直接应用 RUSLE 用户手册48中的算法,10以上的 S 因子采用 LIU 等49基于实测数据提出的算法:L=(/22.13)m(7)m=/(1+)(8)=(sin/0.0896)/3(sin)0.8+0.56(9)S=10.8si
23、n+0.03 516.8sin0.55 15115(11)5)覆盖管理因子 C。植被是土壤侵蚀的重要抑制因素,采用蔡崇法等50基于实测数据提出的算法:C=1f=00.65080.3436lg f0 78.3%(12)式中 f 为植被覆盖度,%。通过采用像元二分模型51对NDVI 数据进行反演获取。6)水土保持措施因子 P。水土保持措施也是土壤侵蚀的重要抑制因素。采用谢怡凡等52关于耕地、建设用地和水域的赋值方法,耕地(不包括梯田)按坡度分级赋值,坡度5、510、1015、1520、2025和25分别赋值 0.100、0.221、0.305、0.575、0.705 和 0.800,建设用地和水域
24、按不发生侵蚀赋值 0。采用秦伟等53关于梯田、园地、林地和草地的赋值方法,梯田赋值 0.084,园地参照鱼鳞坑造林地一并赋值 0.187,未实施水土保持措施的林地和草地赋值 1。2.2.2土地利用变化形式的分解及数据提取在人类活动驱动下土地利用与土地覆被变化关系的分析中,TURNER等54分解并定义了不同土地利用下土地覆被的 2 种变化形式:一是土地覆被的转换,即一种覆被类型完全转至另一种覆被类型,如森林皆伐改种牧草后转至草地,二是土地覆被的改造,即土地覆被属性的改变,如草地过牧后成为退化草地,并且无论是转换还是改造,均需对土地覆被进行维护,如梯田维修、森林抚育等。本研究应用这一土地覆被变化形
25、式分解原理,将土地利用变化分解为土地利用转换和土地利用改造,两者分别指同一空间位置上土地利用类别发生变更和未发生变更的变化。转换区和改造区的面积数据提取62农业工程学报(http:/www.tcsae.org)2023年基于土地利用转移矩阵。2.2.3土地利用和降雨变化对土壤侵蚀影响程度的判别为了判别土地利用变化和降雨变化各自对土壤侵蚀的影响程度,本研究采用 RUSLE 模型情景模拟方法,提出了依据侵蚀模数升降变化量差值的算法:首先以土地利用与降雨共同变化实际情景下研究时段初的侵蚀模数为参照,分别计算研究时段末较时段初土地利用与降雨共同变化实际情景下侵蚀模数的升降变化量(A)、仅土地利用变化模
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 黄土高原 典型 流域 土壤侵蚀 退耕还林 土地利用 变化 响应
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。