土壤水分特征曲线的滞后现象研究综述

土壤水分特征曲线的滞后现象研究综述土壤水分特征曲线是描述土壤水势与土壤含水率之间定量关系的核心曲线，是土壤水动力学研究的基础参数，在农田灌溉管理、水文模型模拟、污染迁移预测等领域具有不可替代的应用价值。在实际观测中，土壤在脱水（由湿变干）和吸水（由干变湿）过程中，相同水势下对应的含水率往往存在显著差异，这种非唯一对应关系即被称为滞后现象。滞后现象的存在会直接导致土壤水分运动模拟的结果偏差，据现有研究统计，忽略滞后效应的入渗模型预测误差最高可达35%，蒸发模型误差甚至超过40%，因此对滞后现象的形成机制、定量表达及影响效应的系统研究，是提升土壤水文过程模拟精度的关键前提。滞后现象的形成机制研究滞后现象的本质是土壤孔隙结构、水-固-气界面相互作用及流体性质共同作用的结果，现有研究对其形成机制的解释主要围绕孔隙几何特性、接触角效应、空气截留及固体基质膨胀收缩四大核心理论展开。墨水瓶效应是最早被广泛认可的滞后形成机制，该理论由Haines于1930年首次提出，认为土壤孔隙由大小不等的孔隙通过狭窄的孔隙喉道连接，呈现典型的“颈粗肚大”墨水瓶结构。在脱水过程中，土壤水从大孔隙开始排出，只有当水势降低到足以克服喉道的毛管阻力时，内部较大孔隙的水分才能流出，因此相同水势下脱水过程的含水率更低；而在吸水过程中，水分首先填充小孔隙，只有当水势升高到能够进入喉道时，大孔隙才会被湿润，导致相同水势下吸水过程含水率更高。后续通过同步辐射显微CT对土壤孔隙结构的三维观测结果证实，土壤中直径小于50μm的孔隙中，墨水瓶结构占比可达42%~67%，是砂质土壤滞后现象的主要贡献因素。接触角滞后效应是解释滞后现象的另一重要理论，指固体表面的前进接触角（液体润湿固体表面时的接触角）始终大于后退接触角（液体脱离固体表面时的接触角）。土壤颗粒表面通常存在有机质、铁铝氧化物等包覆层，且表面粗糙度不均，导致脱水过程中水土界面的后退接触角通常在15°~40°之间，而吸水过程中的前进接触角可达30°~70°。根据毛管上升公式，接触角越大，相同孔隙半径对应的毛管吸力越小，因此在相同含水率下，脱水过程的吸力值显著高于吸水过程。最新研究表明，当土壤有机质含量超过2%时，接触角滞后对整体滞后效应的贡献度可超过墨水瓶效应，成为主导因素。空气截留效应主要发生在土壤快速吸水过程中，当水分入渗速率超过孔隙中空气的排出速率时，部分空气会被封闭在土壤孔隙中无法排出，占据了原本可被水分填充的孔隙空间，导致吸水过程的饱和含水率比脱水过程的初始饱和含水率低5%~15%。这部分截留空气的存在不仅直接降低了相同水势下的含水率，还会改变孔隙的连通性，进一步放大滞后效应。通过氙气示踪法测定发现，在降雨入渗强度超过12mm/h的砂壤土中，空气截留量可达总孔隙体积的8%~12%，是短历时降雨条件下滞后效应不可忽视的组成部分。对于黏粒含量超过30%的黏性土壤，固体基质的膨胀收缩效应是滞后现象的核心成因。黏性土壤中的蒙脱石、伊利石等黏土矿物具有层状结构，吸水时晶层间距扩大，导致土壤体积膨胀，总孔隙度增加；脱水时晶层收缩，孔隙度降低。这种孔隙结构的动态变化使得土壤在干湿循环过程中的孔隙分布本身存在差异，进一步加剧了水势与含水率关系的非唯一性。研究数据显示，膨胀性黏土在饱和含水率下的体积比风干状态下高40%~80%，由此产生的滞后效应占总滞后量的60%以上。滞后现象的定量模型研究为了将滞后效应纳入土壤水分运动模拟，国内外学者已经开发了数十种滞后定量模型，大致可分为经验模型、域模型和物理基础模型三大类，不同模型的适用场景和精度存在显著差异。经验模型是最早发展的滞后表达方法，其核心思路是通过数学函数直接拟合脱水和吸水两条主曲线，再通过插值方法计算主曲线之间的扫描曲线。最具代表性的是Scott等人于1983年提出的三次样条插值模型，该模型将主曲线离散为多个节点，通过三次多项式拟合节点间的曲线段，扫描曲线则通过相邻主曲线节点的线性或非线性插值生成。这类模型的优点是计算简单、能够较好地拟合实测数据，缺点是缺乏物理基础，需要大量实测数据确定模型参数，且无法预测未覆盖干湿范围的扫描曲线。目前经验模型主要应用于农田灌溉制度优化等对模拟精度要求适中、实测数据较充足的场景。域模型（DomainModel）是目前应用最广泛的滞后模型，其核心假设是土壤由大量独立的孔隙域组成，每个孔隙域具有固定的脱水和吸水特征水势，整个土壤的水分特征曲线是所有孔隙域的统计叠加。其中最具代表性的是Philip提出的独立域模型和Mualem提出的相似域模型。独立域模型假设每个孔隙域的干湿过程互不影响，通过孔隙域的水势分布函数计算整体滞后曲线；相似域模型则假设所有孔隙域的滞后因子（脱水与吸水吸力的比值）相同，仅需要主脱水曲线和一个滞后因子即可预测任意扫描曲线。Mualem模型仅需3个参数即可描述滞后效应，且具有明确的物理意义，因此被广泛集成到HYDRUS、SWAP等主流水文模型中。最新研究通过引入孔隙连通性因子对相似域模型进行改进，使得模拟精度较原模型提升了15%~20%。物理基础模型是近年来滞后模型研究的前沿方向，这类模型直接从孔隙尺度的物理过程出发，基于孔隙结构分布、接触角变化、空气截留等机制构建滞后关系。代表性的如基于孔隙网络模型的滞后模拟，通过CT扫描获取土壤三维孔隙结构，构建真实孔隙网络，再模拟单孔隙的干湿过程，进而升尺度得到宏观尺度的水分特征曲线。这类模型的优点是能够揭示滞后现象的微观机制，不需要实测滞后数据即可进行预测，缺点是计算量大，对孔隙结构数据要求高，目前主要应用于基础理论研究，尚未实现大规模的田间应用。随着X射线层析成像技术的普及和计算能力的提升，物理基础模型的应用门槛正在逐步降低，预计未来10年内将成为滞后模拟的主流方法。滞后现象的影响因素研究土壤滞后效应的强弱受到多种因素的共同影响，现有研究表明土壤质地、有机质含量、孔隙结构、干湿循环历史、温度和溶质浓度是最主要的影响因子，不同因子的作用机制和贡献度存在显著差异。土壤质地是决定滞后效应强弱的核心因素，通常表现为砂质土壤的滞后效应弱于壤质土壤，黏质土壤的滞后效应最强。对我国不同地区12种典型土壤的测定结果显示，砂粒含量超过70%的砂质土壤，滞后环的面积（表征滞后效应强弱的指标，单位为cm³/cm³·kPa）通常在0.2~0.5cm³/cm³·kPa之间；壤质土壤的滞后环面积为0.5~1.2cm³/cm³·kPa；而黏粒含量超过40%的黏质土壤，滞后环面积可达1.5~3.0cm³/cm³·kPa。这种差异主要源于不同质地土壤的孔隙分布和矿物组成差异：砂质土壤以大孔隙为主，墨水瓶结构占比低，且黏土矿物含量少，膨胀收缩效应弱；黏质土壤则以小孔隙为主，黏土矿物含量高，多重机制共同作用导致滞后效应显著。土壤有机质含量对滞后效应有显著的增强作用，一方面有机质会增加土壤颗粒表面的疏水性，增大接触角滞后的幅度；另一方面有机质能够促进土壤团聚体的形成，增加土壤中的孔隙异质性，提升墨水瓶效应的贡献。田间试验数据表明，当土壤有机质含量从1%提升至5%时，滞后环面积可增加30%~60%。对于长期施用有机肥的农田土壤，有机质在表层的积累会导致表层土壤的滞后效应显著高于底层，在进行农田土壤水分模拟时需要分层考虑有机质的影响。干湿循环的幅度和次数是影响滞后效应的重要动态因素。首次干湿循环的滞后效应通常最强，随着循环次数的增加，滞后环面积会逐渐减小，一般经过3~5次循环后趋于稳定，稳定后的滞后环面积比首次循环小20%~40%。这种现象主要与反复干湿过程中土壤孔隙结构的调整、空气截留量的减少以及接触角的稳定有关。同时，干湿循环的幅度也会影响滞后效应，当干湿循环仅发生在高水势范围（>-10kPa）时，滞后效应主要由大孔隙的墨水瓶效应和空气截留主导；当循环延伸到低水势范围（<-100kPa）时，小孔隙的接触角滞后和膨胀收缩效应的贡献逐渐增加。温度对滞后效应的影响长期存在争议，早期研究认为温度变化主要通过改变水的表面张力和黏度影响水分特征曲线，对滞后效应的影响可以忽略。但近年的高精度试验结果表明，当温度变化超过15℃时，滞后效应会发生显著变化：温度升高时，水的表面张力降低，接触角滞后幅度减小，同时黏土矿物的膨胀性受到抑制，整体滞后环面积会减小10%~25%。这种温度效应在黏性土壤中表现更为显著，在进行季节性冻土区或变温环境下的土壤水分模拟时，必须考虑温度对滞后效应的修正。滞后现象的应用研究进展滞后现象的研究成果已经在多个领域得到应用，有效提升了相关过程模拟和管理决策的精度。在农田水文过程模拟中，纳入滞后效应的模型能够更准确地预测降雨入渗、土壤蒸发和根系吸水过程。华北平原冬小麦田间试验结果表明，考虑滞后效应的HYDRUS模型对0~100cm土层含水率的模拟纳什效率系数从0.72提升至0.89，对农田蒸散量的预测误差从18%降低至9%，为精准灌溉决策提供了更可靠的依据。在节水灌溉场景中，基于滞后效应优化的交替沟灌制度，能够利用湿润锋推进过程中的滞后特性，减少深层渗漏，提升水分利用效率10%~15%。在水文水资源领域，滞后效应的考虑显著提升了流域产汇流模型的模拟精度。对于南方丘陵区的小流域研究发现，忽略滞后效应的流域水文模型会高估暴雨期的地表径流量20%~30%，低估壤中流的贡献，而纳入滞后效应后，洪水过程线的峰值预测误差从25%降低至8%，峰现时间的预测误差也从2小时缩小至0.5小时。在干旱区地下水补给研究中，滞后效应的存在使得土壤水分的深层渗漏过程存在显著的时间滞后，补给量的计算需要考虑数月甚至数年的干湿循环历史，传统的瞬时平衡法会低估地下水年补给量15%~40%。在土壤污染修复领域，滞后效应会显著影响重金属和有机污染物的迁移转化过程。非饱和带中的污染物迁移与土壤水分运动密切相关，滞后效应导致的含水率分布偏差，会直接影响污染物的对流弥散过程模拟结果。研究表明，在石油污染场地的修复模拟中，忽略滞后效应会导致污染物羽流的扩散范围预测偏差超过30%，修复药剂的注入量设计误差可达40%，考虑滞后效应后能够更准确地圈定污染范围，优化修复方案，降低修复成本。在气候变化影响评估中，滞后效应的存在会改变土壤水分对降雨和蒸发的响应关系，进而影响陆面过程模型对陆地生态系统碳、水、能量循环的模拟精度。最新的CMIP6模式对比结果显示，未考虑滞后效应的陆面模型会高估半干旱区土壤水分的年际波动幅度20%~30%，低估生态系统的干旱抵抗力，导致对植被净初级生产力的预测偏差超过15%。将滞后过程纳入陆面