土壤水分是地球表层多个过程的重要纽带, 直接影响着陆地生态系统的结构、过程和功能^[1]。土壤水分不足是干旱、半干旱地区植被恢复和生态系统服务维持的重要制约因素, 不同时空尺度植被与土壤水分的相互作用关系一直是生态恢复领域的核心科学问题和重点研究内容。尤其在我国黄土高原地区, 土壤水分是植被生长所需的直接水分来源, 土壤水分不足已成为这一区域生态恢复与植被重建的关键制约因素^{[2, 3]}。长期以来, 黄土高原人工植被恢复有效缓解了这一地区严重的水土流失状况, 尤其是大规模退耕还林还草工程的实施, 有效提升了包括土壤保持、固碳等在内的多种生态系统服务功能^[4]。全国生态环境十年变化(2000—2010年)遥感调查表明, 黄土高原是我国植被覆盖增加最为明显的区域^[5], 人工植被已成为这一地区主要的植被类型。然而, 当前的植被恢复和布局方式也导致了土壤水分的过度消耗, 人工植被对土壤水分的长期多度消耗超出了降水的补偿, 导致干燥化现象的出现^{[3, 6-7]}。土壤干燥化是土壤水分含量某一段时间或长期低于自然植被的水文现象。土壤出现干燥化以后难以为地上植被提供充足的水分来源, 不利于人工植被的持续稳定生长^[8], 大范围的土壤干燥化势必影响植被恢复的可持续性和生态系统的稳定性, 对区域生态安全造成威胁。本文针对黄土高原土壤干燥化现象形成的气候背景, 着重从土壤干燥化的植被驱动过程、土壤干燥化时空分异的环境因素以及土壤水分的可持续利用等几个方面, 讨论黄土高原大规模植被建设与土壤干燥化研究面临的问题与挑战。

1 土壤干燥化的气候背景

土壤水分含量是水分输入(降雨)和输出(蒸发、蒸腾和渗漏)综合作用的结果。黄土高原大部分地区年降水量为400—600 mm, 而潜在蒸发量为800—1000 mm^[9], 降水量少而潜在蒸发量大的干旱、半干旱气候是这一地区土壤水分不足的气候背景。另一方面, 黄土高原降水主要集中在夏季, 且多以暴雨形式呈现, 雨强大、历时短, 不利于降水入渗, 使得土壤水分有效补给不足^[10]。从区域分异上来看, 黄土高原从东南到西北降水逐渐减少, 使得土壤水分含量沿着这一降水梯度也逐渐减少^[11], 和降水量与蒸散量的比值相对应^[9], 这也是土壤干燥化在黄土高原区域尺度空间分异的重要影响因素。近几十年来的气象观测资料表明, 黄土高原地区气候存在较大的年际波动, 耦合气象和土壤水分监测的研究也证明了这一地区生长季降水量的减少使得土壤储水呈减少趋势^[12], 这在一定程度上加剧了土壤水分的不足。

2 植被建设对土壤干燥化的影响

植被对土壤水分的作用主要表现在两个方面：一是通过水分利用来影响土壤水分, 尤其是植被生长阶段大量蒸腾耗水对不同深度土壤水分施加影响, 造成土壤持续干化, 二是通过植株形态和群落结构影响降水分配从而影响土壤水分的收入过程^[13-15]。植被在土壤干燥化的形成过程中起到关键作用, 其中人工植被大量蒸腾耗水是土壤干燥化的直接原因。对于黄土高原来说, 人工植被对土壤水分的消耗尤为严重, 土壤水分含量及其季节动态均主要受植被类型的影响^[16]。此外, 植被的生长年限、群落结构、冠层盖度、根系密度等都对土壤干燥化的形成和分布有重要影响^{[17, 18]}。

2.1 植被类型和生长年限对土壤干燥化的影响

人工植被建设改变了植被自然演替的方向, 从而改变了土壤水分的演化过程。在黄土高原水土流失治理的初期, 为了快速达到水土保持和改善生态环境的效果, 植被恢复一般选择抗逆性和抗旱性强、速生高产的外来种和栽植种^[19]，并且种植密度较大。这类植被及其种植方式一般耗水性强, 生长阶段需要大量水分供给, 而黄土高原天然降水不足, 植被为维持生长必然要通过根系不断下延吸收深层土壤水分, 而深层土壤水分难以得到降水的有效补充, 从而使得深层土壤干化难以恢复, 对植被恢复的可持续性与生态功能的发挥形成了制约。黄土高原不同地区的研究报道显示, 人工林草植被可消耗5m、10m、甚至20m以下土层的水分^{[20, 21]}。植被类型是影响土壤水分的重要因素, 人工植被生产力过高、种植密度过大是形成土壤干燥化的首要原因^{[3, 22-24]}。研究表明, 乔木、灌木、人工牧草等不同植被类型对土壤干燥化的影响差异很大, 同一植被类型也因密度、生物量以及生长年限不同, 对土壤干燥化的影响程度也不同^{[25, 26]}。乔木、灌木和人工牧草植被地上部分叶片数量和生物量远高于荒草和农田, 强烈的蒸腾作用使得相同条件下其土壤水分含量明显偏低, 这种不同植被类型之间土壤水分的差异已得到广泛的证实^[27-29]。

一般而言, 人工植被随着生长年限越长, 地上生物量累积越多, 其对土壤水分的消耗也就越高, 在降水不能有效补充土壤水分的时候, 土壤干燥化随着植被生长年限的增加而加剧。例如, 李玉山^[30]早在1983年就发现多年生人工牧草地和茂密林地所形成的土壤水分亏缺, 不论在深度和干燥化程度上, 都远远超过一年生作物地。在黄土高原不同地区, 分别针对刺槐、油松、山杏、柠条、苜蓿、沙打旺、苹果等不同植被类型时间序列的土壤水分对比分析均发现, 人工乔灌草植被生长年限越长, 土壤水分含量就越低^{[26, 31-33]}。但并不是随着植被生长年限延长, 土壤水分含量呈一直减少的状态, 而是在一定年限以后, 土壤水分呈一个较为稳定的波动状态, 并且随着人工植被的生长衰退, 其对土壤水分的过度消耗也会相应减少, 土壤水分会有一定的恢复。

2.2 植被群落特征与土壤干燥化的关系

相比自然植被, 人工植被群落往往结构单一、密度过大、生物量较高, 会消耗更多的土壤水分, 加剧土壤干燥化程度^[18]。例如, 李军等^[34]模拟不同种植密度刺槐林土壤干燥化效应发现, 种植密度越高土壤干燥化速度越快。另外, 黄土高原南部半湿润区农田在中低产年代未发现土壤干燥化现象, 而在高产以后出现土壤水分不足也是一个很好的证明^[9]。王力等^{[35, 36]}对比了子午岭天然林地和人工林地的结构特征和耗水特性, 发现天然林地在结构上形成了典型的“乔—灌—草”复层稳定的空间层片结构, 自我调节能力较强, 只在浅层形成了轻度水分亏缺, 且在雨季后能及时恢复, 并未影响天然植被的发育与演替, 而人工林地则形成较为严重的土壤干燥化。总体而言, 当前黄土高原大规模人工植被恢复中, 人工植被物种单一, 群落结构简单, 容易对土壤水分过度消耗, 反过来土壤水分不足又影响了人工植被的生长发育, 并且单一的群落结构对自然灾害、病虫害的抵御能力也相对较弱。

土壤水分及其空间变异与叶片面积、冠层盖度、叶面积指数的空间异质性密切相关^{[37, 38]}。群落结构、叶片形态等对不同尺度的环境变化具有极强的敏感性^[39], 植被系统在水分限制条件下可以通过调节植被密度、物种组成、叶片面积、叶片数目等来降低对水分需求的压力^[40], 因而地表植被的覆盖度、叶面积、形态结构等功能性状及相应的水分利用过程对土壤干燥化有一个明显的反映^{[41, 42]}, 这使得利用植被动态信息识别土壤干燥化过程成为可能。黄土高原浅层土壤水分时空变异方面的研究积累丰富^[43], 但坡面、小流域和区域尺度上植被动态对土壤干燥化的作用还比较缺乏定量和系统的研究, 这种研究上的不足还难以形成对植被合理空间布局的科学指导。尤其是当前对土壤干燥化时空格局的研究, 还主要停留在野外采样—数据分析比较的阶段, 由于土壤干燥化尤其是深层土壤干燥化野外采样费时费力, 在一定程度上制约了相关研究的发展。虽然已经有学者开展了不同土地利用/植被格局下土壤干燥化特征的研究^{[16, 28, 44-45]}, 但仍多限于立地或坡面尺度。基于空间上的植被群落功能性状、覆盖度等植被动态信息对土壤干燥化时空格局进行定量识别, 通过植被格局优化来调整不同尺度土壤水分状况, 以维持土壤水分的持续利用, 亟待进一步的深入。

3 土壤和地形因素对土壤干燥化空间分异的影响

黄土高原地下水埋藏深, 土壤蒸发和植物蒸腾无法利用地下水, 在降雨不能满足植被消耗的情况下, 植被只能消耗土壤储水, 而这部分水分难以及时补充, 人工植被若过度消耗就容易造成土壤水分的持续亏缺, 而深厚的黄土覆盖是这一地区形成深层土壤干燥化的重要因素。另一方面, 黄土的质地均一, 毛管孔隙发达, 具有极强的蒸发能力, 黄土特殊的土壤属性同样影响了土壤干燥化的形成与分布。在黄土高原, 黄土粒度组成存在由东南向西北逐渐变粗的现象, 从而导致在这一方向上土壤持水性能逐渐降低, 土壤蒸发性能逐渐增强^[46]。不少学者的研究表明土壤质地、容重等的差异也是黄土高原区域尺度土壤干燥化空间分异的重要影响因素, 而小流域尺度的研究也同样发现, 土壤类型和质地的差异影响了土壤干燥化的空间分异^{[26, 47-48]}。

干旱、半干旱地区土壤水分的时空分异与地形密切相关^{[49, 50]}, 尤其在地形破碎、沟壑纵横的黄土丘陵区, 地形是土壤水分空间变异的重要影响因素^{[51, 52]}。地形因素中的坡度、坡向、坡位、坡形等都能通过影响水文路径、太阳辐射等作用于不同深度土壤水分^[53]。多个研究表明, 植被一致的情况下, 随着坡度增大, 坡面径流的流速和流量增加, 入渗减少, 土壤含水量降低, 土壤干燥化的程度就越严重^{[31, 54]}。阳坡受太阳辐射相对较多, 土壤温度高, 蒸发蒸腾作用强烈, 因而土壤干燥化程度也大于阴坡^{[29, 31, 55]}。下坡位因径流的汇集, 土壤水分含量一般高于上坡位, 坡脚、沟底、川道水分条件相对较好, 土壤干燥化程度小于坡顶^{[19, 31]}。也有研究认为地形对土壤水分的作用仅限于浅层土壤^[54], 不一定能有效作用于深层土壤, 深层土壤水分还主要受植被类型的影响^{[56, 57]}。同样也有研究发现, 若不考虑植被因素, 土壤水分主要受地形和土壤影响, 其空间格局呈现出较强的时间稳定性；但若考虑植被因素, 土壤水分的空间格局并不稳定^[58]。比如, 植被生长状况等作为土壤干燥化时空分异的一个重要影响因素, 有时甚至能改变地形条件对土壤水分的影响^{[27, 59]}。

4 土壤水分的恢复与可持续利用

黄土高原的降水入渗深度除南部半湿润区在大丰水年可达到5m以下外, 一般都在3m以内, 且没有深层渗漏^{[9, 60]}, 黄土高原深层土壤水分的恢复较为困难。王进鑫等^[61]通过定位监测研究发现, 自然条件下单纯依赖天然降水难以使人工林地土壤干燥化状况得到有效改善。也有研究发现, 3年生沙打旺草地形成的土壤干燥化至少在沙打旺铲除3年以后方能恢复到最初状态, 而乔木形成的土壤干化在砍伐几十年后才能恢复如初^[62]。干旱、半干旱气候背景下土壤水分的恢复虽然较为困难, 但合理的植被配置、整地措施以及土地管理对土壤水分的恢复也有显著作用。例如, 王晗生等^[23]认为撂荒休闲、间伐或加大树木空间的人工植被管理措施, 同样可以缓和郁闭林分中的缺水状况, 恢复深层土壤水分。也有研究发现种植一年生农作物, 土壤水分能得到一定的恢复, 且恢复的深度和程度随年限的增加而增大。为维持植被生态系统的稳定和土壤水分的持续利用, 当前多数研究均建议黄土高原植被建设应以自然恢复为主, 遵循植被演替规律和耗水规律^{[46, 63]}, 人工植被恢复需要因地制宜, 依据水分承载力调控群落密度和生产力^{[22, 64]}, 并改变单一的人工植被种植方式, 丰富群落结构^[17], 使土壤水分得以维持和改善。梯田、水平阶、水平沟、鱼鳞坑等坡面整地措施在一定条件下可以有效拦蓄降水、增加入渗, 具有改善土壤水分状况、防止和削弱土壤干燥化的作用；另一方面, 坡面整地措施通过对水土进行拦蓄, 还有利于提高植被成活率, 促进植被生长和恢复^{[65, 66]}。采取合理的措施与方法, 干燥化的土壤能够得到一定程度的恢复, 维持水分的可持续利用与生态系统的稳定。黄土高原地区大规模人工植被恢复已有几十年的时间, 而生态系统演化到与环境相适应的状态往往在60年甚至百年以上。在探讨人工植被恢复与土壤水分的可持续利用上, 也同样需要考虑在更长时间尺度上, 人工植被生态系统是否最终会自我调节到与环境相适应的状态, 土壤水分在更长的时间尺度上是否会同人工植被的演替达到相互适应的状态。

5 黄土高原土壤干燥化研究的问题与展望

对黄土高原土壤干燥化形成过程的研究, 从最初的不同植被类型土壤水分的比较, 到土壤水分的动态监测以揭示其形成过程, 也有学者开始通过水分循环过程、蒸散发过程的观测, 以及同位素示踪技术的应用^[67], 对这一过程进行分析, 从降水—土壤水的转化过程、土壤水分的消耗过程等多个方面来解析土壤干燥化过程。例如, Wang等^[68]通过次降雨过程中湿润锋和入渗量的对比分析来研究植被类型对土壤水分收支过程的作用, Fan等^[27]和Chen等^[69]则结合降雨特征, 通过剖面土壤水分动态揭示了土壤干燥化的形成过程。这些研究有效推动了土壤干燥化驱动机制的研究, 但向坡面、流域或区域尺度的扩展还需要更多的案例探讨。

对土壤干燥化动态的研究表明, 浅层土壤水分受降雨、气象条件、植被耗水等影响季节波动较大, 而较深层次的土壤水分则相对较为稳定^{[18, 69]}, 当前对小流域土壤干燥化状况及其动态的研究多以年为单位进行对比分析, 还较为缺乏空间格局及其季节动态的探讨。不同尺度土壤干燥化时空分异规律的解析, 比较多的集中于立地和坡面尺度, 侧重于不同植被类型、群落结构特征、种植年限、生物量等因素影响下的土壤水分的比较。植被随着降雨、温度等的季节变化而存在较大的季节波动, 目前还较为缺乏耦合植被群落动态特征、植被配置结构、植被空间布局对土壤干燥化时空动态影响的研究, 以明确植被动态及其空间异质性对土壤干燥化时空分异的影响。此外, 明确人工植被与土壤水分之间的互馈过程是黄土高原人工植被生态系统可持续性维持的重要基础, 而植被与土壤水分之间的互馈机制尚需要更多长期和深入的研究。

土壤干燥化是多重因素综合作用的结果, 后续研究需要进一步揭示多重环境因素对土壤干燥化时空动态的直接和间接影响并明确其定量关系, 构建植被与土壤水分相互作用的空间模型, 开展变化环境下不同尺度人工植被恢复对土壤水分以及生态系统影响的模拟预测研究。科学的植被配置是黄土高原水土保持和生态安全的重要基础, 对于大规模人工植被恢复区域, 需要基于植被与水分的相互作用关系, 探讨如何合理地配置植被组成和空间布局, 进行科学的植被建设以维持土壤水资源的可持续利用和生态系统服务功能的提升。小流域是黄土高原水土流失综合治理的基本单元, 小流域内不同植被类型的比例如何影响流域整体的土壤水分状况？植被的空间格局及其覆盖动态是否对土壤干燥化过程产生影响？黄土高原小流域的植被覆盖较为复杂, 从流域整体的角度研究复杂植被覆盖格局下土壤干燥化的时空分异规律, 探讨基于植被格局优化的土壤干燥化调控机制, 对于维护植被生态系统稳定和土壤安全具有重要的科学意义和应用价值。