《Water》：Response Characteristics of Karst Water Level to Precipitation and Hydrological Circulation Patterns in the Jinan Spring Basin, Northern China
编辑推荐：
岩溶（karst）地下水系统以高度非均质流网为特征，导致水位对降水变率的响应十分复杂。济南泉域（Jinan Spring Basin）是中国北方最具代表性的岩溶泉系统之一，受降水、地下水开采和水文地质条件变化影响，泉流量发生显著变化。然而，降水信号控制岩溶含水
岩溶（karst）地下水系统以高度非均质流网为特征，导致水位对降水变率的响应十分复杂。济南泉域（Jinan Spring Basin）是中国北方最具代表性的岩溶泉系统之一，受降水、地下水开采和水文地质条件变化影响，泉流量发生显著变化。然而，降水信号控制岩溶含水层水位波动的时间尺度，以及信号在岩溶含水层中传输的机制仍不清楚。研究人员采用30个气象站2016—2018年日降水数据与趵突泉（Baotu Spring）水位记录，运用全局小波谱（Global Wavelet Spectrum, GWS）和小波变换相干性（Wavelet Transform Coherence, WTC）方法，量化主导降水周期、尺度依赖的降水-地下水关系及由相位推导的水位响应滞后。研究识别出17.37、29.22和330.57天三个显著降水周期。地下水响应呈现明确时间尺度依赖性：短周期信号响应快但局部化，中、长周期信号相干性更强且更持久。显著相干面积占比（Percentage of Significant Coherence Area, PASC）从0–17.37天尺度的23.50%增至29.22–330.57天尺度的64.28%，表明累计降水而非单次降雨事件主导水位变化。响应滞后空间分布显示：17.37天快速响应主要位于南部补给区，反映优势通道补给；29.22天响应滞后自南向北渐进增加，指示区域流与储水效应；330.57天响应在大断层附近局部增强，暗示深部循环受构造控制。该研究揭示济南泉域降水信号经多条具不同时序特征的水文循环路径传输，为岩溶含水层结构与功能提供新认识，并为泉水资源的可持续管理提供科学支撑。
研究背景方面，岩溶（karst）含水层由相互连通的管道（conduit）、裂隙（fracture）与基质（matrix）域构成，具有强空间非均质性，降水信号可经优势路径快速传输，也可因扩散入渗与储水过程被延迟和衰减。济南被称为“泉城”，济南泉域是中国北方著名岩溶系统，趵突泉等泉群具有生态、文化与社会经济价值。已有研究多使用相关分析、回归模型与滞后估计等传统统计方法，通常假定单一主导响应时间尺度，难以刻画岩溶地下水循环的多尺度本质。降水-泉流量关系还受地下水开采、断层与地层边界等地质因素扰动，因此降水信号在济南岩溶含水层内的时间传输机制仍未被充分揭示。研究人员以济南泉域为对象，分析2016—2018年30站日降水与趵突泉日水位，旨在区分尺度依赖响应模式、将其空间变异联系水文地质背景，并为护泉与水资源管理提供依据；该论文发表于《Water》。
关键技术方法上，研究人员使用济南降雨水文信息系统提供的30个气象站（P1–P30）日降水数据，以及济南市水务局提供的趵突泉监测站日水位数据，样本覆盖南部山区补给区、中部径流带与北部山前排泄区。方法包括全局小波谱（GWS）识别降水主导周期，小波变换相干性（WTC）在时频域量化降水与水位关联及相位关系，以AR(1)红噪声零假设与1000次蒙特卡洛替代试验做95%显著性检验，用显著相干面积占比（PASC）和平均小波相干性（AWC）比较尺度间关系，并将选定周期相位角转为响应滞后，采用普通点克里金插值制作滞后空间分布图，边缘效应区（Cone of Influence, COI）被排除。
研究结果部分，4.1 主导降水周期性：研究人员对降水序列作GWS，识别出17.37天、29.22天与330.57天三个超95%显著阈值的周期；330.57天接近年降水量循环，但因三年记录仅约三个完整周期，解释需谨慎。4.2 尺度依赖的降水-地下水相干性：WTC显示30站整体平均AWC为0.66、平均PASC为55.05%；0–17.37天带AWC仅0.20、PASC为23.50%，个别降雨事件与水位相干弱且不连续；17.37–29.22天带AWC升至0.25、PASC至30.00%；29.22–330.57天带AWC达0.56、PASC达64.28%，P1与P30等站PASC超70.39%，说明累计降水信号与水位关系最强最持久。4.3 水位响应时间滞后的空间分布：由WTC相位导出滞后，17.37天尺度滞后1.76–3.04天，刘长山—郎猫山—万灵山带滞后较长，南部与东南如P12、P25较短；29.22天尺度滞后约5.42–6.65天，自南部补给区向北部排泄区递增，P9为5.43天、P30为6.44天；330.57天尺度滞后77.18–86.88天，千佛山断层（Qianfoshan Fault）附近P25与P17出现局部短滞后异常，反映断层可作优先水力通道。
讨论部分总结，研究人员指出岩溶含水层具低通滤波效应：短时标降雨波动被衰减，长时标补给信号被保留，传统方法仅给单一平均滞后，而小波方法能分辨17.37天与29.22天等被月、年尺度分析忽略的组分；2016—2018年日开采与人工回补未作为独立协变量，滞后应理解为自然与管理共同条件下的降水-水位响应滞后。短周期滞后受裸露碳酸盐岩与岩体侵入控制，中周期滞后受北倾单斜与区域流路径控制，长周期局部异常得到水化学与同位素证据间接支持：张夏组（Zhangxia Formation）深部循环水年龄10–15年，浅层水小于6年，断层可沟通深浅单元。概念模型分为浅层循环（17.37天，管道与大裂隙流）、中间循环（29.22天，裂隙网与基质储水）、深层循环（330.57天，断层控制水力窗与季节补给）。
结论部分翻译为：理解降水信号如何在不均质岩溶含水层中传播，对揭示地下水循环机制与提升泉系统可持续管理至关重要。研究人员以中国北方典型岩溶系统济南泉域为对象，集成GWS、WTC与基于相位的响应时间分析，系统评价降水-地下水作用的时空特征。（1）GWS识别出17.37、29.22与330.57天三个主导周期且均超95%置信水平；330.57天近似年降水量循环，但三年记录仅约三个完整周期，解释需谨慎。（2）WTC表明水位波动与累计降水信号耦合更强；PASC从短尺度23.50%增至长尺度64.28%，说明济南岩溶含水层通过储水与迟滞水力传输衰减高频降水、保留长期补给信号。（3）17.37天响应快（1.76–3.04天）受局部补给与优势流路控制；29.22天滞后南小北大，反映区域流与储水；330.57天滞后77.18–86.88天，千佛山断层等构造致局部异常。（4）降水-地下水关系可解释为快速浅层岩溶流、裂隙控制中间流与长期区域循环组成的层级系统，为不均质岩溶含水层信号传播提供新视角。
