利用双曲模型和相位依赖模型诊断热带岛屿流域的表观线性退水特征

《Journal of Hydrology X》：Diagnosing apparent linear recession signatures using hyperbolic and phase-dependent models in tropical Island catchments

　　

主退水曲线(Master Recession Curves, MRCs)呈现的表观半对数线性通常被解释为均质线性水库行为的证据。然而，基于物理过程的地下水理论预测，随着流域耗竭过程中储水量、导水系数(transmissivity)和水力连通性(hydraulic connectivity)的演变，退水动态可能呈现非线性。表观MRC线性是代表结构不变的流域属性，还是聚合退水过程的涌现特征，这一问题仍未得到解决，特别是在以陡峭地形、快速径流产生和异质地下水流路径为特征的热带岛屿流域中。本研究利用印度尼西亚东部安汶岛(Ambon Island)五个热带季风流域的217次径流退水事件来回答这一问题。研究人员采用分层多诊断框架(hierarchical multi-diagnostic framework)，整合了全局指数退水建模(global exponential recession modelling)、双曲退水分析(hyperbolic recession analysis)、相位特定退水因子估计(phase-specific recession-factor estimation)、残差诊断(residual diagnostics)、连续分段回归(continuous segmented regression)、信息论模型比较(information-theoretic model comparison)以及退水陡度指数(Recession Steepness Index, RSI=?ln k)的不确定性感知解释。尽管全局指数模型始终取得优异的描述性能(R²=0.979–0.987)，但所有流域均表现出系统性的相位依赖特征(k_earlylate)、时间上有组织的残差曲率，以及分段模型获得显著更强的支持。统计支持的断点集中在约7至10天之间，表明有效退水行为存在可重复的转变。这些发现表明，较高的全局半对数拟合优度可能掩盖内部演化的耗竭动态，且单一退水常数不一定能提供流域退水的结构完整表示。由于地下水观测、导水系数测量、示踪剂证据、蒸散发估算和海岸边界效应未得到独立量化，推断的机制仍是过程一致的假说(process-consistent hypotheses)，而非直接的水文地质验证。总体而言，本研究提供了一个可重复、具有不确定性感知的框架，明确区分了描述性模型性能、结构性退水行为和物理过程解释，从而加强了对热带岛屿流域低流量持续性和储水–出流动力学的解释。

**热带岛屿流域退水动态的多诊断解析：表观线性背后的相位依赖机制**

**研究背景与意义**

退水曲线(master recession curve, MRC)分析是流域水文学中刻画无雨期径流消退过程、评估低流量持续性和储水–出流关系的经典方法。传统研究通常将半对数空间中呈现近似线性的主退水曲线(MRCs)解释为均质线性水库行为的证据，认为径流以时间不变的退水常数按比例随流域蓄量衰减。然而，基于物理过程的地下水理论（尤其是Boussinesq型公式）预测，随着流域耗竭过程中饱和厚度、导水系数和水力连通性的演变，储水–出流关系可能呈现非线性特征，退水速率应随储水状态而变化。因此，表观MRC线性究竟是结构不变的流域属性，还是多个非线性退水过程聚合后产生的涌现特征，仍是流域水文学中悬而未决的核心问题。该问题在热带岛屿流域尤为突出，因为此类流域以陡峭地形、快速径流产生、强季风降雨、深风化火山物质、裂隙基岩和高度动态的储水连通性为特征，其退水轨迹可能整合地下水排水、土壤水耗竭、河道与河岸储水释放、蒸散发及海岸边界影响等多种过程。该研究发表在《Journal of Hydrology X》上。

**研究内容与结论**

研究人员利用印度尼西亚东部安汶岛五个热带季风流域（Wae Batu Merah、Wae Ruhu、Wae Batu Gantung、Wae Batu Gajah和Wae Tomu）2010–2020年的日径流和降水记录，经严格筛选获得217次独立退水事件，构建了分层多诊断分析框架。该框架整合了全局指数退水建模、双曲退水分析、相位特定退水因子估计、残差诊断、连续分段回归、信息论模型比较(AICc)以及退水陡度指数(RSI=?ln k)的不确定性感知解释。

**关键技术方法**

研究人员采用的主要技术方法包括：(1)基于四重筛选标准的退水事件识别（日降雨量<1 mm day^?1、径流单调递减、持续≥10天且无测量异常）；(2)事件归一化、累积和对齐的复合MRC构建方法（采用中位数聚合）；(3)基于Levenberg-Marquardt算法的非线性最小二乘双曲模型拟合；(4)10,000次事件级bootstrap重采样量化参数与断点不确定性；(5)基于hinge函数的连续分段回归结合校正Akaike信息准则(AICc)、Akaike权重进行结构模型比较；(6)剖面似然置信区间估计断点位置不确定性。

**研究结果**

**3.1 双曲退水作为基于过程的框架：超越表观对数线性**

五个流域的复合MRC在半对数空间中呈现光滑、单调且近平行的轨迹，全局指数模型的R²为0.9789–0.9872，表明单一指数模型可作为描述性的一阶近似。然而，双曲模型Q(t)=Q₀(1+αt)^?m的拟合优度更高（R²=0.999977–0.999993），其曲率参数m在所有流域均接近1（0.9902–1.0034），但该结果对应的是双曲衰减形式而非时间不变的指数水库行为。

**3.2 相位依赖退水敏感性**

所有流域均表现出k_earlylate的系统性模式，即早期退水阶段比例耗竭速率较快、后期逐渐减缓。相位敏感度Δk=k_late?k_early范围为0.0201–0.0289，bootstrap 95%置信区间均不包含零，表明相位依赖具有统计显著性。优化断点集中在7.60–10.49天。

**3.3 残差曲率诊断**

全局单斜率模型的残差并非随机分布：四个流域呈现不同程度的U形曲率，Wae Batu Merah以集中后期偏差为主。残差的系统时间组织表明单一退水斜率在结构上不足以代表完整退水轨迹。

**3.4 信息论结构检验**

分段模型的ΔAICc值为?45.55至?65.34，远低于?10的强支持阈值，Akaike权重几乎完全支持分段表示。断点集中在约7–10天，与相位分析结果高度一致。

**3.5 退水陡度指数与低流量弹性**

RSI=?ln(k)=?b=τ^?1，范围0.0420–0.0506 day^?1，对应特征退水时间尺度19.75–23.84天。Wae Batu Gajah具有最低RSI和最长持久性，但其置信区间与Wae Ruhu部分重叠，Wae Tomu与Wae Batu Gantung之间无法进行确定性排序。

**讨论与结论**

讨论部分指出