第42章 从入渗到状态记忆:降水怎样成为径流
42.1 1933年的入渗问题
雨落到土地之后,并不立即成为河流中的水。它可能停在叶面,进入土壤,充满洼地,沿地表移动,也可能在更晚的时候经地下路径抵达河道。如果只记下雨量与出口流量,这些路径便被折叠在一个比例里。比例有时有效,却很容易把不同的发生过程当成同一种转换。水文学的重要问题之一,就是怎样打开这个比例,识别雨水进入不同路径的条件。
霍顿1933年的《入渗在水文循环中的作用》把入渗放在这一转换的中心。[C02] 进入土壤的水由此成为分配后续路径的关键环节:地表供水有多少能够被接纳,剩余部分又以何种方式积留与移动?本章以该论文提出的问题为历史起点,以下采用现代水文学概念分析入渗、储留与产流;这些教学推演不等同于原论文的公式、实验程序或术语原貌。如此分层,既保留1933年的问题位置,也能说明后来形成的分析尺度。
理解这一转向,可以从能力与实际通量的区别开始。入渗能力表达在相应条件下地表允许水进入土壤的速率,实际入渗还受到供水约束。地表没有足够水,即使土壤能够接纳更多,实际通量也不会凭空增加。因此,雨强小于入渗能力时,不能把能力直接当成进入土壤的水量;雨强超过相应能力时,超过的部分也要先考虑洼地储留与其他路径,才能讨论出流。
这一差别使研究从“降了多少雨”转向“何时以多大强度降雨,地表当时处于什么状态”。总雨量的量纲是长度,雨强与入渗速率的量纲则是长度除以时间。前者可以比较一次事件的总体规模,后者才能检验特定时刻是否出现供水超过入渗的条件。如果不先区分量纲,研究者容易把一个总量阈值当成过程阈值,解释也就失去明确的作用位置。
土壤表面还不是一个恒定孔洞集合。含水程度、表层结皮、植被根系与可连通孔隙,会影响进入路径;积水压力与水分再分配又可能改变实际入渗。雨水到来之前的土地因此不能用“同一块土”四个字完全概括。同一地点具有连续身份,却可能处于不同状态。水文学将状态纳入解释,并不是向外部添加一个抽象的“记忆”,而是承认先前过程留下的物理条件仍在参与眼前过程。
这里也须防止另一种简化:入渗能力被超过不是所有地表径流的唯一形成机制。某些地方由于土壤已经接近饱和,或者地下水抬升与地形汇聚而出现地表出水,即使局地降雨强度没有超过一块较干土壤的潜在接纳能力,也可能有径流。不同机制可能同时存在,并随区域、季节与事件阶段改变其份额。把一个机制讲清楚,不能以删除其他机制为代价。
以发生学语言说,研究者面对的S不再只是某场雨后的总出流,而是具有前期条件与路径约束的响应事件。D包括雨强、土壤水分和河流过程的顺序差异。E的关键不在“土壤、雨、河流都有关”,而在土壤含水改变可用储留空间,储留与入渗改变地表积水,积水和连通路径再改变到达河道的时序。每一项的解释身份由其他项共同规定。
这种拆分也改变了调查安排。总量比例只要求比较雨量与出口出流;要区分入渗受限、储留填满与路径连通,就需要同时记录雨强、前期含水和产流位置。变量的增加有明确分工:它们分别约束输入速率、可用容量和空间连接。思想史所要追踪的,正是这一问题转换怎样使原先未被记录的状态进入观测,而非把所有新增变量当作同一种资料积累。
还可以把地表产流与河道出口出流进一步分开。地表某处出现积水,不意味着水立即抵达测站;它可能先填满小洼地,随后连接到浅沟,再进入河道。连通性具有阈值:孤立积水区尚未连成路径时,局地有水而出口响应有限;降雨继续后,原来分离的区域连通,出口可能迅速增加。于是,响应的非线性不必全由某个土壤参数突然变化造成,也可能来自空间连接的形成。
这使“贡献面积”成为随事件阶段变化的对象。地形流域面积可以保持稳定,而真正向河流提供快速径流的区域并不固定。将二者混同,会使研究者误以为整片土地始终以同样比例响应降雨。更细的解释需要辨认坡脚、近河区与高处土壤在不同状态下的作用,但是否采用动态贡献区,仍要由观测和任务决定。
入渗、贡献面积与汇流时延共同打开了雨量—出流比例。比较两场雨时,首先需要保持总量相同,再分别改变输入速度、前期状态或空间连接,才能看清究竟是哪一项使响应分化。
42.2 相同雨量的不同过程
以下用教学情境隔离总量与速率。甲、乙两场雨都提供三十毫米降水,甲均匀分布在六小时内,乙集中在半小时内,相应平均雨强分别为每小时五毫米与六十毫米。若土壤入渗和地表储留能够接纳甲的输入,却无法及时接纳乙的输入,两场雨便可能产生不同的积水与径流。三十毫米仅用于保持事件总量一致,不代表某地预警阈值;能否形成快速出流,还取决于积水是否接入河网。
现在再改变前期条件。一次较小降雨落在已经湿润的土壤上,可能比一次较大降雨落在干燥土壤上产生更大的即时出流份额。然而这一判断也不是处处适用的经验律。干燥土壤的表层状况、疏水性、裂隙与结皮可能使实际响应偏离简单想象。正确的问题不是“湿土一定比干土产流多吗”,而是哪些可测状态在什么地表条件下改变了输入到输出的分配。
洪峰还取决于汇流同步性。若两个支流分别产生同样的出流总量,但峰值到达干流的时刻相差很大,下游峰值可能被摊开;如果到达时刻重合,峰值则可能增加。雨量空间分布与河网传播时间因此共同参与洪峰形成。只比较全流域平均雨量,会把“下在哪儿”这一差异消去。对年度水量而言可以容忍的平均化,对短时洪峰而言可能恰好删掉最有决定力的条件。
这给事件比较提出了严格要求。首先,雨量必须与所研究的汇水面积相匹配。其次,降雨开始、结束与出流观察窗口需要明确。再次,前期状态不能只由一个方便的统计指标代表而不检验其适用性。最后,要说明出口流量已经经过何种水位—流量转换,洪峰期间是否超出了实测率定范围。否则,“同量不同流”可能部分源于资料误差,而不是新机制。
一个有解释力的事件表应同时保留降雨过程、前期状态与输出过程。降雨过程记录强度、持续时间和空间覆盖;状态记录土壤含水、积雪、水库或其他相关储量;输出过程则记录峰值、到达时间、总量与退水形状。不能把这些指标全部揉成一个总分,因为不同机制在不同指标上留下的痕迹不同。一个模型可能正确预测总量,却错误预测峰值;另一个模型可能抓住峰值,却误判退水所反映的储留时间。
若要判断加入状态变量是否有用,应先提出它应改善哪一部分响应。把前期湿度加入模型,可能改善径流总量;引入河网传播可能改善峰值时序。若声称前期湿度同时解释全部变化,就需要更强证据。模型里增加变量总会增加调整空间,只有在没有参与调参的完整事件上保持改善,才较有理由认为抓住了可迁移关系。
同样重要的是“没有改善”的结果。如果土壤湿度观测只代表一小片地表,其信号可能与真正参与洪水形成的区域不一致。此时无改善并不自动证明土壤状态无关,而可能说明测量对象与作用对象失配。反过来,统计改善也可能只是湿度变量代理了未明确记录的季节性。需要通过季节、土地覆盖和降雨类型的分组检查,把机制意义与预测便利分开。
同量事件的推演得到的是一种可区分关系:总雨量保持不变时,改变输入历时或传播时延,峰值仍可变化。真实流域中的主导原因则要靠完整过程线判断。如果短历时输入增强而峰值时序基本保持,产流分配值得先查;如果总出流变化较小而峰形显著展宽,汇流时延就是重要候选。两者也可能同时改变,因而需要把总量、峰值和时序作为相互约束的结果。
用卷积思想还能说明汇流怎样改变洪峰。在一个线性近似的教学模型中,出口响应可表示为此前各时刻有效输入,经由不同传播时延加权后的叠加。权重函数若较窄,输入集中到达,峰值较尖;若较宽,水被分散到更长时间,峰值较低而持续较久。即使输入总量完全相同,传播时间分布也能改变过程线。这个推演隔离了汇流作用,并不假定真实流域必然线性。
如果权重函数随前期状态变化,问题又进一步复杂。较湿状态下部分路径可能更快连接,原本适用的传播关系便改变了。此时仅调整一个产流比例,可能正确改变总量却仍误判峰值。把产流与汇流的误差分开,能够避免一个参数代替两个不同机制。
退水段提供另一种信息。降雨停止后仍有出流,说明储存的水继续释放。在无新输入的线性储库教学近似中,设Q=kW,W为可排储量,Q为体积流出率,k的单位为时间的倒数;由dW/dt=-kW,得到W与Q按指数衰减。释放时间尺度为1/k。退水斜率因此能够约束有效释放时间;不过,多储库叠加或随状态变化的k也会改变退水形状,仅凭局部拟合尚不能唯一辨认内部储库数目。
降雨与流量的比较还涉及面积加权。三个雨量站的算术平均,不一定代表流域平均降水,因为它们覆盖的地形与面积不同。若高处降水显著而站点集中于谷地,平均值可能系统偏低。对同量事件的筛选应保留面雨量估计方法,否则所谓相同总雨量可能只是站点平均相同。
蒸发损失在不同事件时距中的地位也不同。对很短的强降雨过程,它可能不是峰值形成的主导项;对长时段湿润与干燥演变,它却影响前期状态。一个过程在事件方程中可以近似忽略,并不意味着它在准备该事件的状态形成中无关。区分直接作用与通过前期状态进入的作用,有助于解释为什么同一变量在不同模型中时而显眼、时而隐去。
事件结束也不是纯粹日历判断。若前一场雨的退水尚未结束,后一场雨已经开始,两场事件的水在出口叠加。把它们人为切开后分别计算径流系数,可能把储量延续误当作某场雨的异常效率。完整过程分析应保留事件边界的物理理由,而不是只按固定间隔划分。
42.3 历史记忆怎样被误读
水文学中的记忆至少有三种含义。第一种是物理状态的延续:先前降雨使土壤、水库或地下含水系统保留了水。第二种是统计表达:模型用前几个时刻的变量预测后一个时刻。第三种是资料处理造成的伪记忆:训练和检验之间共享了本应隔离的信息。这三者可以同时出现,却绝不能混为一谈。
物理记忆具有可追问的载体。问“流域记得上一场雨”,实际应问什么储库改变了、改变维持多久、通过哪条路径影响本次出流。储库的时间尺度不同,记忆也就不同。浅层土壤可能对近期降雨敏感,深部地下系统可能保持更长的响应;河道中的传播与水库调度又产生其他时间尺度。一个滞后项能够在数学上描述迟到,不意味着已经识别迟到的物理位置。
统计记忆有其正当价值。即使不能逐项测量内部储量,出口流量的历史序列也可能含有可预测信息。但历史系数必须在研究范围内解释。某个模型使用过去三天流量,并不证明自然系统的记忆恰好为三天;三天可能只是资料分辨率与模型结构的折中。若换一个时间步长,系数随之改变,也不能直接说自然机制发生变化。参数身份依赖模型形式,需要与可观测状态建立桥梁。
时间泄漏则属于证据组织错误。比如把同一场洪水的相邻时刻随机分配到训练和检验集合,模型已经从训练片段看过该事件的大部分结构,检验结果却被解释成对一场全新洪水的预测。这种方法可以回答局部插值问题,却不能直接回答新事件泛化问题。若研究目标是未见事件,就应按完整事件或完整时段划分,并保证标准化、参数选择和缺测填补不提前使用检验期信息。
更隐蔽的泄漏发生在特征构建中。所谓“前期降雨指标”如果包含预测时刻之后的雨量,名称再合理也失去预报资格。某个经过居中平滑的流量序列,可能同时利用过去与未来时刻;在事后重建中这可以有用,在实时预警中却不可直接使用。需要分别说明信息何时真实可得、何时经过处理进入数据库,以及模型在何时作出判断。
验证设计因此应先于结果比较。如果研究者在看过一组结果后不断改变事件边界、剔除困难样本或选择最有利时段,最后漂亮的性能不能承担原本的证据责任。预先写清分组原则并保留失败事件,既是统计纪律,也是思想史研究的重要提醒:被保留下来的成功模型往往有清晰故事,被排除的失败条件却可能揭示其真正范围。
误差之间还会相关。一个流域连续数天的流量误差可能来自同一场雨量低估,不能把每小时误差都当作相互独立的失败或成功。对有效样本量的判断应考虑事件关联。否则,几场长期过程会被拆成数千个貌似独立的证据,置信程度随记录频率虚增。数据更密不一定意味着独立信息更多,这是水文学与古气候研究共同面对的问题。
从SDE的差异序列D来看,物理演变、仪器采样、数据库处理和模型读取各有自己的次序。它们之间需要明确映射:一小时平均值归属于哪个时段,缺测补值用了哪些邻近资料,预测时刻以前究竟有哪些信息可用。E则指储量与通量维持响应的具体约束,例如同一含水状态通过排水改变退水,又被新降雨改写。将这两类关系分清,才不会把处理程序引入的连续性当作流域自身的记忆。
状态重建还存在“等效状态”问题。两个流域时刻可能具有相同出口流量,却有不同土壤与地下储水分布。它们在没有新降雨时的短期退水相似,遇到下一场雨却可能分化。出口流量因此可以携带状态信息,却未必完整标识状态。模型若只依赖一个滞后流量,容易把这些不同内部配置压成同一输入条件。
这种压缩是否可接受,取决于预测时距与事件类型。若短时预测已经足够准确,使用有效状态可以是合理封装;若某类新事件反复失败,就应检查压缩是否丢失关键差异。研究进步不是永远拒绝黑箱,而是知道何时黑箱的已验证范围被越过。打开内部结构的理由应来自具体失败模式,而非对复杂机制的抽象偏好。
前期降雨指数也是一种状态近似。它通常让较近降雨获得较大权重,较远降雨逐渐衰减。这样的指数能够概括湿润程度,却不直接测量土壤水分;相同指数在蒸散、土层深度或排水条件不同的地方,可能对应不同真实状态。把它当作方便变量可以,把它当作普适物理量则需要额外论证。
42.4 何时需要重建对象
模型失误之后,最节制的解释通常是先检查测量、参数与已有过程。若某段异常可以由雨量空间代表性不足解释,就没有必要立即宣称流域机制改变。若补入一个物理上已有的前期湿度状态即可改善独立检验,首先应肯定成熟水文学的解释能力,而不是把改善归功于新的哲学框架。
但也有情况需要重建对象。例如同一出口的响应在土地开发、排水工程或水库运行改变后呈现系统差异。要把差异归因于土地利用,至少需要变化的空间与时间记录,并排除同时发生的气候、测站与管理变化。仅凭某时期径流增加,不能判定不透水面积是唯一原因。对象重建要求新的连通路径确实进入了作用链,而不只是参数数值变大。
重建也未必意味着推翻整个流域模型。研究者可以把原先合并的地表分成不同响应区,可以把人工调度从自然出流中显式分离,也可以为变化前后设置不同结构。应根据证据决定改动层级。把所有失败都解释为“模型不够复杂”,会使任何反证被更多自由度吸收;把所有变化都塞进一个固定参数,也可能使模型失去物理意义。
一个可检验的重建应留下旧、新模型的共同比较面。它们应使用相同的训练信息和相同的独立事件,报告总量、峰值、时序等不同指标,而不只挑选新模型最有利的一项。新增结构还应给出独特预期,例如某类短历时降雨响应改变明显而另一类事件改变较小。若新模型只提高总体拟合,却不能产生区别性预期,机制解释仍然有限。
洪峰异常与响应区重划之间可以建立一条待检验的研究链:某类事件反复失配,促使研究者区分原先合并的区域;分区提出新的雨量或土壤观测需求;新观测再决定哪些区域可以合并、哪些必须区分。这里的关键是后一步能否修正前一步。如果新增测站并未改变分区判断,便只能确认原有结构;若要把这条研究链写成历史经过,还需模型版本、调查安排与事件记录相互印证。
水文对象的成熟,表现为能够把失败分层处理。有的失败要求纠正数据,有的要求调整参数,有的要求增补状态,有的才要求重建连通结构。这样的区分让理论既能够修正,又不至于每遇困难就失去身份。下一章将把状态与通量的区别带入海洋:温度和盐度可以帮助辨认水团,但要说明水怎样流动,还需要另一类约束。
本章引用的来源条目
C02 Horton, R. E. (1933). The Rôle of infiltration in the hydrologic cycle. Eos, 14(1), 446–460. DOI:10.1029/TR014i001p00446 原始来源入口