一个现象值得注意:做遥感蒸散发的人很少关心根区水力再分配,做土壤水文的人很少用到冠层辐射传输,做大气边界层的人几乎不看土壤大孔隙流。三个群体各自发了大量论文,但他们研究的其实是同一杯水从天空到根尖的旅程。这篇文章想讨论的,不是某个模型或某颗卫星,而是一个更根本的问题:我们的学科分工方式,是否正在成为理解农业水循环的最大障碍?
每一本水文学教材都会提到SPAC(Soil-Plant-Atmosphere Continuum,土壤-植物-大气连续体)。这个由Philip在1966年正式提出的概念,至今已近六十年。它的核心思想极其简洁:水从土壤经植物到大气的运移,是由一条连续的水势梯度驱动的,可以用类似于欧姆定律的电阻网络来描述。
问题在于:我们嘴上说着「连续体」,手上做的却全是「切片」。
看看我们的主流研究范式就知道了。遥感反演蒸散发时,多数方案处理植被的方式是一个「大叶模型」或经验性的作物系数Kc——冠层被压缩成一个没有垂直结构的平面。SWAT、VIC等水文模型模拟入渗和壤中流时,植物根系的存在通常被简化为一条「根区深度」参数——根系不主动吸水,它只是一个「水桶的深度」。而陆面过程模型(如Noah-MP、CLM)虽然在植被生理上做得更精细,却往往把土壤水力特性简化为van Genuchten曲线的几个参数,对非均质性视而不见。
每个模型都声称自己「考虑了SPAC过程」,但实际上它们做的是:把连续体拆成三截,各自简化到可以独立求解,再用几个通量项松散地「粘」回去。这不是耦合,这是缝合。Philip如果看到今天的模型,他会说:你们用最精密的数学,求解了一个最粗糙的假设。
这不是在苛责任何一个模型团队。每一步简化在它的历史语境中都是合理的——当年没有高分辨率遥感,没有足够的计算能力,没有通量塔和宇宙射线中子探测器。但今天,我们的观测能力已经发生了量级跃迁,而模型的底层架构却还停留在二十年前的范式中。
— · —
二、遥感的「三明治困境」:看得见顶和底,看不见中间
让我们把遥感能提供什么、不能提供什么做一个冷静的盘点。
大气端:我们有了再分析资料(ERA5)和多源融合降水产品(IMERG、CHIRPS),辐射通量可以通过CERES卫星获得。大气强迫的时空分辨率在过去十年提升了一个数量级。
地表端:Sentinel-2给了我们10米分辨率、5天重访的光学影像;SAR可以穿透云雨;无人机LiDAR可以重建冠层三维结构。我们对「地表长什么样」的知识,前所未有地丰富。
然而,中间发生了什么——几乎是一个黑洞。
从冠层截留降水开始,到水分穿过叶面边界层蒸散,到穿透雨沿树干流和穿冠雨到达土面,到入渗锋推进,到根系从不同深度按水势梯度吸水,到多余的水分沿优先流路径快速下渗或形成壤中流——这整条链条的中间过程,遥感几乎全部缺席。
一个具体的例子
SMAP卫星的土壤水分产品,标称深度是地表以下5厘米。而华北平原冬小麦的有效根区深度在返青期可达60-100厘米,抽穗灌浆期甚至超过150厘米。也就是说,遥感「看到」的那一薄层水分,与真正驱动作物蒸腾和产量形成的根区水分,之间隔着一个巨大的推断鸿沟。用5厘米的信息去约束150厘米的过程,这不是数据同化,这是信息的过度杠杆化。
SIF(日光诱导叶绿素荧光)的出现曾让很多人兴奋——终于有一个信号直接关联光合过程了。但SIF能告诉你的是「碳泵」的速率,而不是「水泵」的阻力。从SIF到气孔导度到蒸腾,中间还隔着叶片水势、气孔响应函数、边界层阻抗等一系列非线性环节。把SIF直接回归蒸散发,在统计上可行,在物理上可疑。
— · —
如果我们真的要把SPAC当回事,那么水文模拟与遥感的耦合方式需要彻底重构。不是「遥感提供输入 → 模型运行 → 遥感验证输出」这种流水线式的合作,而是让遥感信息在模型的每一个时间步内、每一个垂直层上都参与约束。
这意味着什么?至少意味着三个方向的突破:
① 从「状态同化」走向「过程同化」
目前的数据同化(EnKF、4D-Var等)本质上是在修正「状态变量」——土壤含水量、叶面积指数、雪水当量。但真正需要被约束的往往是「过程参数」——根系水力导度、土壤大孔隙比例、冠层消光系数。状态可以随时间漂移回去,而过程参数的修正才能改变模型的行为模式。这要求我们发展出能够联合同化多源遥感观测来约束过程参数的框架,而不是把每种观测分别对应一个状态变量。
② 从「像元思维」走向「功能单元思维」
30米像元内可能包含三种作物、两种土壤质地和一条田间沟渠。用这个像元的平均NDVI去驱动一个假设均质的水文响应单元,误差不是来自传感器,而是来自我们对「空间单元」的定义方式。更合理的路径可能是:用高分辨率遥感先识别「水文功能单元」——按照产流机制、蒸散发特征和地下水补给关系划分的单元,而非行政边界或固定网格。让水文过程本身来定义空间离散化的方式。
③ 从「单向驱动」走向「双向反馈」
在真实的SPAC中,植被不是被动的「水分消耗者」。根系会主动向干燥区域生长(向水性),气孔在午间会因为水力安全阈值而主动关闭(等水策略 vs 非等水策略),作物在长期干旱下会调整根冠比。这些生物反馈机制使得植被既是水循环的「被驱动者」,也是「调节者」。当前多数模型将LAI作为外部遥感输入来强制驱动蒸散发计算,实际上切断了「水分胁迫 → 叶面积缩减 → 蒸散发下降 → 土壤水分恢复」这条负反馈回路。正确的做法是让LAI成为模型的预报变量,而遥感LAI作为观测算子参与同化。
— · —
最后不得不谈谈房间里的大象:深度学习。
Kratzert等人用LSTM在美国数百个流域上取得了超越传统水文模型的径流模拟精度,这个结果是真实的、可复现的。但有两个被广泛忽略的前提:
第一,LSTM学到的是「输入-输出映射」,不是「过程机理」。它可以预测明天的径流量,但无法告诉你这些水是从地表来的还是从基流来的,更无法告诉你改变灌溉制度后径流会怎样响应。而农业水资源管理恰恰需要的是后者——可干预、可推演的因果知识。
第二,深度学习的成功高度依赖训练数据的代表性。美国有CAMELS数据集,有数十年的高质量径流、气象和属性数据。而中国西北的内陆河流域、南亚的季风灌区、非洲的雨养农业区——这些农业水问题最严峻的地方,恰恰是数据最贫瘠的地方。在数据荒漠中,物理模型的先验知识是不可替代的脚手架。
我的判断是:未来十年,真正有突破力的工作不会来自「纯物理模型」阵营,也不会来自「纯数据驱动」阵营,而是来自那些能够在物理模型的骨架上嵌入可学习组件的混合架构。用神经网络去参数化那些我们写不出解析表达式的子过程(如根系吸水的空间分配、优先流的触发条件、作物气孔对VPD的非线性响应),同时保持质量守恒、能量守恒等硬约束不被违反。differentiable modeling这个方向已经有了一些开创性的工作,但在农业水文领域的应用才刚刚起步。
— · —
回到开头的问题。农业水文模拟与遥感的真正瓶颈,表面上是分辨率不够、参数不准、算力不足,本质上是我们还没有建立起一个真正的「植被-土壤-大气」一体化建模范式。
我们有了越来越多的卫星、越来越高的分辨率、越来越强的GPU,但我们对水分穿过一个玉米冠层到达根区这件事的理解,可能并不比三十年前进步多少。不是因为没人研究,而是因为研究冠层的、研究根区的、研究边界层的,长期在不同的期刊、不同的会议、不同的学科评价体系中各自运转。
如果这篇文章只能留下一句话,希望是这句:水不在乎你是搞遥感的还是搞水文的,它只沿着势能梯度走自己的路。我们的模型和观测体系,也该学会用水的方式思考。
— END —