《HydroResearch》：Connectivity over proximity: a directional analysis of aquifer–river coupling across three baseflow-dominated catchments
编辑推荐：
摘要
含水层与河流的耦合控制着基流的产生、流域记忆及干旱可预测性，然而其在一个含水层内部的空间强度与方向却鲜有被解析。本研究对英国三个基流主导型流域系统——两个白垩（Chalk）含水层系统与一个鲕粒灰岩（oolitic limestone）系统——中的17
摘要
含水层与河流的耦合控制着基流的产生、流域记忆及干旱可预测性，然而其在一个含水层内部的空间强度与方向却鲜有被解析。本研究对英国三个基流主导型流域系统——两个白垩（Chalk）含水层系统与一个鲕粒灰岩（oolitic limestone）系统——中的17口共置观测井，采用相关分析、超前-滞后分析（lead–lag analysis）、谱相干性（spectral coherence）、传递熵（transfer entropy）及部分信息分解（partial information decomposition）等方法，对去季节化异常（deseasonalised anomalies）进行了无模型（model-free）的定向耦合分析。研究发现，耦合在所有系统中均为正相关，但与井至流量站或河道距离之间均无一致关系：West Beck、Dour和Coln三个流域的月尺度中位数相关系数分别为0.87、0.57和0.40，而耦合强度与距离的秩相关在不同系统间符号相反。每个系统中最具信息量的观测井均非距离最近者，该井兼具强耦合与净地下水→河流信息传递的特征；在以降雨为条件进行控制后，连通良好的观测井仍保留该信息。因此，相对于几何邻近性，连通性（connectivity）能更好地识别能够反映河流行为的观测井，并为监测网络设计提供支撑。
**连通性主导而非邻近性主导：跨三个基流主导型流域的含水层-河流耦合定向分析——论文解读**
**一、研究背景与科学问题**
在透水性流域中，河流在近似意义上可视为地下水位的表层展现。无雨期间，当坡面流和浅层地下快速径流消退后，河川径流几乎完全由地下水排泄维持，水文过程线的缓退段与含水层的蓄水量实为同一耦合系统的两个侧面。这种耦合的强度与时序决定降雨信号被存储、平滑和缓慢释放的程度，进而决定流域在干旱期间维持径流的能力（Sophocleous, 2002）。由于平滑程度取决于岩石、土壤和地表堆积物的渗透性与储水性，耦合关系编码了流域继承自地质条件的“记忆”，并由此决定低流量与水文干旱的可预测性。
尽管基流指数（baseflow index）等集水区尺度指标已能较好地刻画流域地质条件与径流补给特征的总体关系，且标准化地下水水位指数和地下水响应时间等指标可量化含水层的记忆效应，但现有研究仍存在两个关键局限：其一，地下水对径流的贡献多通过基流分割从径流过程线本身反推，而非直接利用独立观测的井水位数据加以验证；其二，大多数分析在时间上是对称的——互相关、相干性或共享谱峰仅能刻画两个序列共变的强度，而无法区分孰先孰后。对于监测网络的设计与解释而言，方向性恰恰是核心问题：一口水位波动领先于基流的钻孔携带着对河流具有时间领先性的信息，而一口仅是滞后复现河流信号的钻孔，无论相关度多高，均不具信息增益。基于上述认识，本研究提出两个核心假设：H1（几何假说）——在基流主导型系统中，井-河耦合强度不受几何邻近性（距流量站或河道距离）组织；H2（连通性假说）——耦合强度特别是方向性受水力连通性组织，而地下水异常的去相关时间（记忆）可作为连通性的可观测代理指标。
**二、研究内容与主要技术方法**
研究人员选取英国三个基流主导型流域系统作为研究对象：Yorkshire Wolds白垩含水层上的West Beck（NRFA 26009站）、East Kent白垩含水层上的Dour（NRFA 40033站）以及Cotswolds倾坡鲕粒灰岩含水层上的Coln（NRFA 39020站），共17口共置观测井。数据来源为英国国家河流档案库（National River Flow Archive）的逐日流量记录和观测井水位档案，时间跨度覆盖1963–2024年不等。
主要关键技术方法包括：（1）去季节化处理——月尺度数据去除气候态月均值，日尺度数据采用STL（Seasonal-Trend decomposition by Loess，局部加权回归季节趋势分解）方法提取残差异常；（2）耦合强度量化——采用去季节化序列的同期Pearson相关；（3）超前-滞后分析——计算日尺度互相关函数在±120天滞后范围内的峰值位置与显著性（以1000次循环移位替代数据构建零分布）；（4）谱分析——采用Welch周期图法估计交叉谱相干性与相位，窗长为2048天、50%重叠；（5）传递熵（transfer entropy）——基于等概率三态离散化（并验证二态与四态稳健性），以循环移位替代数据检验显著性，净传递熵为两方向差值；（6）部分信息分解（partial information decomposition, PID）——采用I_min测度将所有井对关于河流的联合信息分解为独特、冗余和协同三部分；（7）基流分割——采用Lyne–Hollick递归数字滤波（三遍，α=0.925）和Eckhardt双参数滤波（α=0.98, BFImax=0.80）两种独立方法。
**三、研究结果**
**3.1 耦合强度与距离的有限作用**
三个系统均为强基流主导型流域，Lyne–Hollick基流指数介于0.92–0.95之间，Eckhardt基流指数约0.80。West Beck九口井的月尺度耦合中位数为0.87，Dour四口井为0.57，Coln四口井为0.40；日尺度耦合排序一致。最关键的发现在于：月耦合强度与距流量站距离的Spearman秩相关在三个系统中均不显著且符号不一致（West Beck ρ=?0.60，Dour ρ=+0.80，Coln ρ=+0.40）。Dour系统提供了最鲜明的例证——距流量站最近的三个井（1.1–2.2 km）反而是耦合最弱者（月r介于0.51–0.62），而最远的Wolverton井（3.5 km）耦合最强（r=0.90）。在West Beck日尺度上，耦合强度与距离无关（ρ=?0.40，p=0.51），但与地下水记忆强相关（ρ=?0.90，p=0.037），即记忆越长（信号越被阻尼）的井耦合越弱。
**3.2 耦合的时标、方向与信息结构**
月尺度上几乎所有井的最大相关滞后均为零，说明含水层与河流之间的交换在一个月内基本完成。日尺度互相关函数则揭示：连通良好的井（West Beck的Elmswell Hall和Gameslack、Dour的Wolverton）在零滞后±1天内达到峰值，与河流同步；而弱耦合井（Dalton Estates滞后19天、Willy How Bottom滞后32天）在显著负滞后处达到峰值，即河流异常先于地下水异常。弱耦合井与基流的相关性显著高于与总流量的相关性（Dalton Estates从0.56升至0.76，Willy How Bottom从0.32升至0.46），表明它们主要追踪水文过程线中的慢速含水层补给组分。谱分析进一步显示：连通井（Gameslack）在整个解析时标范围内与河流高度相干且相位接近零；弱耦合井（Willy How Bottom）相干性低且相位随频率变化，呈现阻尼、延迟和仅部分共享的响应特征。
传递熵方向性分析是区分连通井与弱耦合井最具诊断力的指标：West Beck中Gameslack净传递熵为正且显著（+0.012 bit，地下水→河流），Elmswell Hall为双向均显著，而Dalton Estates和Willy How Bottom显著为负（?0.009和?0.005 bit，河流→地下水）。该符号在二、三、四态离散化及一至七天历史长度变化下均保持稳定。将传递熵按源滞后解析发现：Gameslack在零滞后即向河流传递信息，而弱耦合井在约两周滞后处接收河流信号。条件传递熵分析（以流域降雨为条件）表明：五个最强耦合井在条件化后仍保留104–128%的无条件定向信息（p≈0.003），而两个河流主导的West Beck井不显著（p≈0.22–0.24），排除了纯降雨共同驱动的解释。部分信息分解揭示：Gameslack与Willy How Bottom配对时贡献约80%的关于河流的独特信息，与Dalton Estates配对时约50%；弱耦合井几乎不提供独特信息，其贡献几乎完全与连通井冗余；两个连通井之间约70%的信息冗余。全网络分解确认：每个系统中持有最大独特信息项的井均为较远者（West Beck的Gameslack、Dour的Wolverton、Coln的Ampney Crucis），而近流量站井主要贡献冗余信息。
**3.3 跨系统综合**
三个系统的中位月耦合强度从West Beck（0.87）经Dour（0.57）降至Coln（0.40），表面上呈现岩性排序，但两个白垩系统之间的差异（0.87 vs 0.57）与白垩和灰岩之间的差异（0.57 vs 0.40）几乎相当，表明岩性不能作为耦合强度的干净分类变量。真正跨系统稳定再现的是一种内部结构：每个系统均存在少数井兼具高日耦合与显著正的净传递熵（领先河流），而大多数井弱耦合且携带小或河流指向的通量。两个最清晰的领先井——West Beck的Gameslack与Dour的Wolverton——虽相距两百公里、归属不同流域，但特征近乎一致：日耦合系数分别为0.76和0.72，净传递熵分别为+0.012和+0.014 bit，且两者均为各自系统中较远的井而非最近井。
**四、讨论与结论**
讨论部分指出，上述观测结果支持以水力连通性为核心的组织机制：河流优先整合含水层中水力连通性较好部分的响应；连通良好的井几乎与河流同步响应并提供关于河流的独特信息，而连通性较差的井响应缓慢，主要反映同一补给信号的延迟传播。这一解释不依赖单一诊断指标，且与白垩含水层的双重孔隙-双重渗透性结构以及侏罗纪鲕粒灰岩的异质水文地质特征一致。但需强调，研究采用的连通性代理指标（去相关时间）仅从钻孔记录本身导出，不依赖河流观测，属地下水序列的内在动力学属性，避免了循环论证。然而由于缺乏逐井的导水系数、储水系数、示踪试验等独立水力验证，该推断仍属间接证据。Dour系统因受地下水开采和城市化改造影响，其弱耦合井行为与West Beck存在差异，但总体上H1在三个系统中均获得一致支持，H2得到方向性与部分信息分解结果的支持，而记忆作为连通性代理的适用性存在系统依赖性。
研究对监测网络设计具有直接启示：连通性应被评估而非从邻近性推断。建议将超前-滞后分析与信息论测度相结合，作为监测站点评估的补充判据。本研究尚未评估预测能力，不主张任何干旱预警能力。最高优先级的后续研究方向为：对已识别的领先井开展样本外低流量事件预测技能评估（包括事件判别如ROC曲线、概率验证如Brier分数及有效提前时间），同时以独立水力数据验证连通性代理，并将分析扩展到更多岩性类型。
