郑粉莉, 张加琼, 刘 刚, 范昊明, 王 彬, 沈海鸥
(1.西北农林科技大学 水土保持研究所 黄土高原土壤侵蚀与旱地农业国家重点实验室,陕西 杨陵 712100; 2.中国科学院 水利部 水土保持研究所, 陕西 杨陵 712100; 3.沈阳农业大学,辽宁 沈阳 110866; 4.北京林业大学, 北京 100083; 5.吉林农业大学, 吉林 长春 130118)
全球黑土以其土质肥沃和粮食高产对人类的生存与发展做出了巨大贡献。但过度开发利用导致全球黑土地均出现了严重土壤侵蚀。1934年发生在美国的“黑风暴”事件,不但导致全美生态灾难,也敲响了全球农田侵蚀灾害的警钟;20世纪50年代前苏联发生的“黑风暴”事件导致前苏联的粮食危机。中国东北黑土区包括黑龙江省、吉林省、辽宁省和内蒙古自治区,其面积约1.03×106km2,是国家最大的商品粮生产基地,商品率占全国的1/3以上,是国家粮食的“压舱石”[1]。 然而,仅有百余年垦殖历史的东北黑土地,由于高强度的开发利用和掠夺式经营,黑土发生严重退化,黑土生产力降低,平均黑土层厚度已由20世纪50年代60~80 cm下降到目前的20~40 cm[1];每年仅因土壤侵蚀造成的全年粮食减产高达14.7%[2],严重威胁中国粮食安全和生态安全。作为土壤侵蚀“癌症”表现的各类沟蚀在东北黑土区也非常严重。据中国水土流失与生态安全(东北组)报告[1],东北现有侵蚀沟2.50×105余条,侵蚀沟占地面积4.83×105hm2。农田侵蚀沟的发展使“黑土地破碎”,阻碍农机作业,限制农业可持续发展。更为严重的是由于缺乏有效的土壤侵蚀调控技术措施,黑土地土壤退化呈继续恶化的发展态势,是中国唯一仍呈退化加重的农田生态系统。如不及时保护黑土地,将极大削弱区域粮食生产能力和区域生态功能,导致中国最大的粮仓出现问题,影响国家粮食安全和中国人的生存问题。
近几十年来,中国的土壤侵蚀机理研究取得了重要进展, 为水土流失综合治理提供了重要的科学依据。然而,过去的研究主要针对单一侵蚀营力,而对多外营力相互作用的土壤复合侵蚀研究相对较少;尤其与其他水蚀区相比,中国东北黑土区的土壤侵蚀机理研究相对薄弱,从而一方面导致现有的适用于中国其他水蚀区的水土保持措施在该区域的适用性受到质疑,另一方面也导致该区土壤侵蚀治理效果不显著。因此,迫切需要揭示东北黑土区复合土壤侵蚀过程,阐明多营力复合侵蚀的驱动机制,为黑土区土壤复合侵蚀防治和黑土资源的可持续利用提供重要理论依据。
近10 a余来,作者团队通过主持国家重点基础研究发展计划(973)项目课题和多项国家自然科学基金项目以及“十三五”国家重点研发计划项目等,在东北黑土区坡面土壤侵蚀主要驱动机制及其关键因子对坡面土壤侵蚀影响,坡面侵蚀—沉积空间分布规律,黑土区土壤可蚀性估算模型,保护性耕作措施对坡面土壤侵蚀阻控,土壤侵蚀强度对作物产量的影响,适用于坡面复合土壤侵蚀的水土保持措施筛选,复合土壤侵蚀防治分区等方面开展了大量研究工作,取得了重要研究进展。本文基于大量野外调查和理论研究以及文献资料,并集成作者团队长期的研究成果,综合分析东北黑土区坡面土壤侵蚀的特点,以期为研究该区复合土壤侵蚀有效防治措施提供重要科学依据。
受特殊的地理环境和人为过度垦殖的影响,东北黑土区坡耕地最显著的侵蚀特点是多营力作用的复合侵蚀在时间上的更替和空间上的叠加。坡耕地晚春季节发生融雪侵蚀与昼融夜冻互相交替的复合侵蚀(图1),随后发生耕作侵蚀和风蚀叠加的复合侵蚀(图2)。
东北黑土区坡耕地多营力复合侵蚀在时间尺度上表现为季节性的更替。东北黑土区坡耕地夏季和秋季会发生严重的水蚀(图3);晚秋季节由于整地又发生耕作侵蚀以及冬季和早春季节土壤冻胀作用引起土壤颗粒位移和临空部分土体位移(图4)。在空间尺度上表现为水蚀(降雨和融雪径流侵蚀)、风蚀、冻融侵蚀并存,且相互叠加和相互作用,加剧了坡面土壤侵蚀的严重程度。
图1 东北黑土区晚春季节的融雪侵蚀现象 (郑粉莉摄于20170427)
图2 东北黑土区晚春季节的耕作侵蚀与扬尘现象 (张加琼摄于20170425)
图3 东北黑土区夏季的坡耕地水蚀(左:范昊明摄于20180625;右:郑粉莉摄于20180729)
图4 东北黑土区晚春季节的土壤冻胀裂隙与 土体位移现象(刘刚摄于20190526)
黑龙江省沟蚀数量由20世纪50年代4 000余条增加到2000年的15 000余条[1]。其中,黑龙江省沟蚀严重的克东县和拜泉县东部10个乡镇,侵蚀沟数量已从1965年的1 712条增加到2005年的2 565条,侵蚀沟面积占总面积的比例已由5.97%增加到17.62%,对应的沟壑密度也由202 m/km2增加到283 m/km2[1]。东北黑土区14县的21个典型小流域沟蚀调查表明,2005年的沟壑密度较之2004年增加了12.2%[1];仅经过2004年一个雨季,调查的7条小流域的沟头溯源侵蚀速率介于0.05~5 m,平均为1.34 m;沟床加深变化于1.6~9.0 cm之间,平均为5.8 cm;侵蚀沟加宽介于4~145 cm,平均为36.8 cm[1]。农田侵蚀沟的发展不但使“黑土地破碎”,阻碍农机作业,也使耕地面积减少。据张兴义等在美丽中国生态文明建设科技工程专项研究中的调查资料,60%以上的侵蚀沟分布在耕地(图5),损毁黑土地5 000 km2以上。作者野外调查也表明,在克山县英民沟流域仅2018年的春季融雪径流使侵蚀沟沟头平均前进了1.9 m,沟深增加了3.8 cm;在2018年的雨季侵蚀沟沟头平均前进了2.6 m,沟深增加了7.8 cm。同时,坡耕地浅沟侵蚀也非常严重。
受长缓坡地形和顺坡垄作的影响,坡面汇流引起的土壤侵蚀更加明显。东北黑土区坡耕地坡长一般为上百米至上千米,大多介于300~500 m。刘宝元等的观测资料表明,对于顺坡垄作种植大豆的坡面径流小区,坡长为300 m的土壤侵蚀量(2 006.5 t/hm2)是坡长为20 m土壤侵蚀量(5.6 t/hm2)的213倍[1]。加之顺坡垄作形成的一系列“垄—沟”系统,加大了坡面汇流作用对坡面侵蚀的贡献。室内模拟降雨试验表明,对于无垄作裸露坡面,在50或100 mm/h降雨强度下,当坡面汇流速率由10 L/min增加到20 L/min,坡面侵蚀量增加了1.2~4.7倍[3];而对于顺坡垄作坡面,同样在50或100 mm/h降雨强度下,当坡面汇流由10 L/min增加到20 L/min时,坡面侵蚀量增加1.4~7.6倍[4];充分说明了在东北黑土区特殊的长缓坡地形条件下,顺坡垄作明显加大了汇流对坡面土壤侵蚀的作用。
图5 东北黑土区坡耕地的沟蚀现象(左:郑粉莉摄于20170426;右:郑粉莉摄于20170622)
雨滴打击和径流侵蚀是引起坡面侵蚀的主要动力。与中国其他水蚀类型区相比,雨滴打击是东北黑土区坡面土壤侵蚀的主要驱动力。黑土坡面雨滴打击作用对坡面侵蚀的贡献可达到72%~96%,反映了雨滴打击对颗粒的分散占主导作用[5-6];而不同区域的类似研究结果表明,黄土区雨滴打击对坡面侵蚀的贡献占35%~61%[7],红壤坡面雨滴打击对坡面侵蚀的贡献为27%~76%[8]。野外原位模拟降雨试验也表明[9],当黑土坡面秸秆覆盖量仅为2 kg/m2,即可减少87%的坡面侵蚀量,进一步证明了雨滴打击作用是黑土区坡面侵蚀的主要驱动力。同时,由于东北黑土层下覆母质多为黄土状亚黏土或湖相沉积物,土壤质地较为黏重,加之长期耕作形成了坚硬的犁底层,容易导致在雨季和融雪过程中形成壤中流。壤中流的形成极大降低土壤颗粒间的黏结力,使土壤可蚀性增加,从而加剧坡面土壤侵蚀。已有研究结果表明,与自由入渗条件下相比,坡面壤中流的形成使坡面侵蚀量增加了68.3%~74.3%[5]。
目前水蚀和风蚀研究取得了重要成果,并研发了各类水蚀和风蚀模型,为世界范围的土壤侵蚀防治做出了重要贡献。有关复合土壤侵蚀研究,目前主要集中在融雪和降雨复合侵蚀、冻融—水力复合侵蚀等方面,同时在冻融作用对土壤可蚀性影响等方面取得了一定的进展;而有关冻融—风力复合侵蚀,水力—风力—冻融复合侵蚀方面的研究还十分薄弱。
冻融作用与风蚀、水蚀、重力侵蚀等共同作用或交替发生,使复合土壤侵蚀过程不同于单一营力侵蚀过程。当前在冻融作用影响水蚀的方式、速率和土壤团聚体稳定性等方面的研究都取得了重要进展,而对冻融—风力复合侵蚀研究还相对薄弱。有关的研究报道多为冻融作用对土壤风蚀速率和土壤风力可蚀性的影响以及冻融与风蚀的相互关系研究[10-14]。而目前冻融—风力复合侵蚀研究大多集中在风蚀较严重的干旱或半干旱区,而有关湿润和半湿润区冻融—风力复合侵蚀研究还鲜见报导,有待进一步加强。
风水复合侵蚀是风力和水力对同一对象的共同或交替作用,从而产生有别于单一营力(风力或水力)的侵蚀、搬运和沉积过程。按照风力和水力的耦合方式,可分为动力耦合和媒介耦合两种。动力耦合是风力作用于雨滴形成“风驱雨”(wind driven rainfall),从而改变雨滴的侵蚀动能。这实质是风力影响下的水蚀, 风力和水力具有时间和空间同步性。媒介耦合是风力和水力对可蚀性下垫面的交错或交互叠加作用[15]。 风力和水力的交错或交互叠加作用,导致一种侵蚀营力对地表物质的侵蚀、搬运和沉积改变另一种侵蚀发生的物质基础和地表条件,从而形成有别于单一营力作用的复合侵蚀[16]。因此,媒介耦合式的风水复合侵蚀在时间上交错或交互叠加发生,在空间上同地或异地发生。这种模式的风水复合侵蚀较“风驱雨”侵蚀更为复杂,是目前研究的难点和重点。目前在风水复合侵蚀的区域划分、时空变化、侵蚀强度、侵蚀能量特征、驱动因素、过程和机理等方面取得了重要进展[10,17-21]。然而有关研究仍然存在以下几方面的不足: ①研究区域多为干旱区河流—沙丘系统的风水复合侵蚀,而较少涉及半干旱和半湿润区风水复合侵蚀。由于半干旱或半湿润区的地形地貌、土壤和地表物质组成、降水(尤其是降雪)特征、植被条件和土壤含水量动态变化等与干旱和半干旱区完全不同,导致风水复合侵蚀更为复杂,且半干旱和半湿润区大多是粮食主产区,更需要明确风水复合侵蚀过程和机理,以便为风水复合侵蚀防治和维持粮食产能提供科学依据。②关于风蚀和水蚀相互作用的研究有待深入。目前的研究更注重两者的相互促进效应,而对二者的相互抑制效应重视不够[15]。风力的强烈分选作用导致地表粗化,从而影响水蚀过程中降水入渗、径流路径和流速,进而影响水蚀速率[22]。另一方面,水蚀形成沟谷加剧风力磨蚀,从而加剧风蚀[23];但细沟发育也可能增大地表粗糙度,又可减小风蚀。因而,需要加强风蚀和水蚀的相互作用研究。③风水复合侵蚀过程和机理研究亟需加强。由于风蚀和水蚀在发生范围、影响区域、侵蚀物质输移方向与维度、风力和水力的搬运能力、侵蚀事件的时间尺度等方面均存在巨大差异[10],而当前风水复合侵蚀机理多为简单的定性描述[18-19],有待进一步加强。
多外营力复合侵蚀除冻融—水力侵蚀、冻融—风力侵蚀、风水复合侵蚀外,还包含冻融—风力—水力侵蚀、冻融—水力—重力侵蚀、风力—水力—重力侵蚀等。目前,针对多营力复合作用的土壤侵蚀研究较多关注两个侵蚀营力相互作用, 而对三个及以上营力相互作用下的侵蚀过程研究极少[10-11, 20]。中国的东北黑土区正是三个及以上侵蚀营力复合作用的侵蚀区域,尤其是嫩江下游的1.04×105km2(包括黑龙江西南部、内蒙古东部和吉林西部),受水力、风力、 冻融、融雪、重力等作用的共同影响,土壤侵蚀剧烈,但有关多营力复合侵蚀研究进展非常缓慢[24-25]。已有研究主要是通过核素示踪、实验模拟、模型估算、GIS/RS 等方法,初步区分了风蚀和水蚀对坡面土壤侵蚀的贡献,明晰了降雨径流和融雪径流对土壤侵蚀的作用,分析了冻融作用对土壤可蚀性相关土壤性质的影响,明确了不同季节黑土坡面侵蚀的主要影响因素[26-32]。但目前关于多营力复合侵蚀的理论还十分薄弱,有待今后进一步加强。
本文基于以上东北黑土区复合土壤侵蚀特点和研究现状的分析,提出黑土区复合土壤侵蚀研究的重点领域。
(1) 多种外营力相互作用的坡面复合侵蚀过程机制。目前大多研究聚焦单一侵蚀外营力作用下的坡面侵蚀特征和过程研究,但有关不同季节坡面侵蚀过程与主控影响因子研究,以及多种外营力互作过程机制等研究相当薄弱。因此,当前迫切需要强化复合土壤侵蚀的理论研究,揭示多种外营力相互作用的复合侵蚀过程机制。
(2) 冻融作用对土壤抗侵蚀能力的影响机制。虽然这方面的研究取得了重要进展,但冻融作用对土壤抗侵蚀能力的影响机理、融雪侵蚀过程土壤抗侵蚀能力的动态变化机制等方面的研究还相当薄弱。当前研究重点应是冻融作用下土壤可蚀性动态变化机制,土壤抗侵蚀能力对冻融作用的响应机制等。
(3) 区分和量化冻融作用、融雪、降雨径流和风力侵蚀对坡面侵蚀和沟蚀的贡献。由于各种侵蚀外营力相互叠加、相互影响,且彼此互作机理及作用贡献尚不清楚。因此区分多种侵蚀营力对坡面片蚀、细沟和浅沟侵蚀过程的主导作用,量化各侵蚀营力对切沟侵蚀的贡献是当前需要强化的研究领域。
(4) 复合侵蚀作用的流域尺度泥沙来源的诊断。由于多营力相互作用和相互叠加,东北黑土区流域泥沙来源更加复杂,需要建立针对这种特殊侵蚀环境的复合指纹示踪方法,以便能有效地诊断流域泥沙来源,为流域水土保持规划提供重要科学支持。
(5) 复合土壤侵蚀预报模型研发。国内外现有成熟的土壤侵蚀模型,皆是针对单一营力的水蚀或风蚀预报,目前尚未有适用于复合侵蚀的预报模型,因此亟待开展研究。
(6) 防治复合土壤侵蚀措施的有效性评价。数十余年来,国家先后启动了一系列东北黑土区水土流失综合治理试点工程。但由于理论研究滞后于生产实践,从而使治理工程实施效果不佳。因此,迫切需要开展复合土壤侵蚀防治措施评价工作,科学回答保护黑土地应采取哪些技术措施等重大技术问题,为国家宏观决策和工程实施提供重要科技支撑。
(7) 东北黑土区复合土壤侵蚀防治分区。黑土区属于多营力叠加的复合土壤侵蚀区,但现阶段中国尚未有关于黑土区复合侵蚀分区的图件,严重影响了黑土侵蚀防治工作的开展。因此亟待开展东北黑土区复合土壤侵蚀防治分区研究。