《Water》：Hydrogeological Response of a Karst Aquifer to Extreme Recharge: Insights from the 2026 Grazalema–Líbar Hydro-Seismic Crisis (Southern Spain)
编辑推荐：
极端降水可深刻改变岩溶含水层中的地下水动力学特征，尽管控制其对异常补给响应过程的认识仍不充分。特别是，含水层充填过程中水力连通性与承压性的演化可能强烈影响地下水压力传播及相关地质过程。本研究调查了2026年1–2月一次异常补给事件期间，西班牙南部部分承压的格拉
极端降水可深刻改变岩溶含水层中的地下水动力学特征，尽管控制其对异常补给响应过程的认识仍不充分。特别是，含水层充填过程中水力连通性与承压性的演化可能强烈影响地下水压力传播及相关地质过程。本研究调查了2026年1–2月一次异常补给事件期间，西班牙南部部分承压的格拉萨莱马–利巴尔（Grazalema–Líbar）岩溶含水层的水文地质响应，该事件期间累计降水量超过2800 mm。研究人员整合了气象观测、泉流量记录、水文地质信息及地震数据，以分析地下水补给、水力行为及地震活动。含水层响应迅速，表现为地下水水位大幅上升、泉流量显著增加、先前不连通区域间水力联系逐步建立，以及承压条件扩展。这些变化提高了水力扩散率，并促进了压力通过岩溶管道网络的快速传递。最强的地震响应仅在此水力重组发展之后出现，地震集中于承压区段并逐步向浅部迁移。结果表明，水力连通性的演化（而非补给本身）控制着压力传递与地震激活。该研究为部分承压岩溶含水层在日益频繁的极端降雨事件中的水文-地震响应提供了概念框架。
极端降雨触发岩溶含水层水文-地震响应的机制：2026年格拉萨莱马–利巴尔水文-地震危机研究解读
一、研究背景与科学问题
极端降水事件正日益成为级联地质灾害的重要驱动因素，其影响远不止于洪水和边坡失稳。强烈降雨可改变地下水力条件，产生瞬态应力扰动，进而可能影响地震活动。水文诱发地震活动已在水库蓄水、季节性地下水波动、人工补给、废水注入及地热作业等情境中被记录。然而，关于强降雨后局部地震活动增强（即降雨触发地震）的案例研究仍然较少，部分原因在于此类降雨本身亦属极端事件。值得注意的是，已报道的案例均发生在岩溶地区，这提示岩溶化作用可能在将水文应力传递至地壳的过程中发挥关键作用。
尽管已有研究报道了降雨与地震活动之间的时间相关性，但自然补给过程触发地震的机制仍未完全阐明。大多数研究将孔隙压力扩散视为主导触发机制，类似于水库诱发地震。然而，高度岩溶化的含水层可能表现出其他地质介质中不常见的额外过程，包括管道快速水力传导、高水位期间的瞬态承压、水力扩散率的大幅波动以及局部超压发育。这些机制的相对重要性仍不确定，且鲜有研究利用整合的气象、水文地质和地震学数据集加以探讨。
格拉萨莱马–利巴尔（Grazalema–Líbar）岩溶系统位于西班牙南部，为研究上述过程提供了独特的天然实验室。2026年1–2月，该地区经历了伊比利亚半岛有记录以来最异常的降雨事件之一，累计降水量超过2800 mm。该事件产生了异常的地下水补给、前所未有的泉流量增加、大幅地下水位上升以及含水层内承压条件的逐步扩展。与此同时，在不到两周的时间内发生了超过350次地震的强烈地震序列。同一岩溶系统内广泛的无压、半承压和承压区段共存，为考察水力重组如何影响压力传播及相关地震活动提供了独特机遇。
二、研究设计与方法概述
研究人员综合运用气象学、水文地质学和地震学观测数据，重建了2026年1–2月异常补给事件期间格拉萨莱马–利巴尔岩溶含水层的时空演化过程。气象数据来自西班牙国家气象局（AEMET）；泉流量记录和地下水水位数据来自安达卢西亚自治区水资源总局；地震数据来自西班牙国家地震台网（IGN）。研究方法包括：表征极端降雨事件的气象强迫特征；利用泉流量记录和水文地质观测评估地下水响应；基于泉衰退系数估算水力扩散率；分析地震活动的时间和空间演化；将所有观测整合为描述压力传播与地震激活水文地质过程的概念模型。其中水力扩散率的估算采用Sanz de Ojeda等提出的方法，将泉的衰退系数与含水层特征水流长度相关联。
三、主要研究结果
3.1 水文地质响应
含水层的水文地质响应迅速而强烈。钻孔水位快速上升，多处溢流泉（trop-pleins）出现，表明含水层大部分区域的水位上升超过300 m。主要泉点（包括加托泉、贝纳奥汉泉、希梅拉·德·利巴尔泉和查尔科·德尔·莫罗泉）的流量均出现显著增加。地下水水位上升伴随先前仅部分连通的区段逐步建立水力联系。特别重要的是，随着水位上升，利巴尔含水层中先前无压的区段转变为半承压或承压状态，从而促进了压力在更远距离上的传递。
3.2 有效水力扩散率的演化
利用泉衰退系数和代表性水流路径长度估算的水力扩散率显示，孤立的格拉萨莱马含水层与水力整合的利巴尔含水层之间存在近两个数量级的差异。格拉萨莱马含水层的扩散率约为0.007 m2/s，以无压地下水条件为主；利巴尔含水层无压区段的扩散率介于0.75至1.2 m2/s之间，反映了其发达管道网络带来的更高水力连通性。研究推断，极端补给事件期间，地下水位大幅上升增加了含水层饱和厚度，逐步连接了先前不连通的管道网络，并扩展了低渗透地层下方的承压和半承压区段，从而增强了压力传递效率。
3.3 地震响应
2026年2月4日至16日期间记录了超过350次地震，而此前2025年1月至2026年1月全年仅检测到35次地震。地震活动在最大补给后迅速增加并持续约两周。地震序列可分为两个截然不同的阶段：第一阶段（2月4–8日）包含浅源和深源地震，震源深度可达约35 km，古腾堡–里希特b值接近1（0.96 ± 0.07），符合区域构造地震特征；第二阶段（2月8–16日）以浅源地震为主，约90%的事件发生在1 km深度以内，b值升至1.20 ± 0.12，反映了低震级事件比例的相对增加。地震密度图显示，地震活动与利巴尔含水层存在强烈的空间关联，最高密度出现在瓜迪亚罗河谷下方的承压和半承压区段，以及贝纳奥汉、希梅拉·德·利巴尔和查尔科·德尔·莫罗泉附近的主要地下水排泄区。
四、讨论与结论
研究提供了降雨触发地震活动的综合证据。降雨本身并非直接触发机制，地震活动是含水层对异常补给渐进水文地质响应的结果。主降雨事件与强烈地震活动开始之间数天的延迟表明，地下水入渗、水位上升和压力传播需要时间才能改变既有断层上的应力条件。地震活动集中于利巴尔含水层的水文地质域内，而非均匀分布于整个研究区，且震源深度逐步向浅部迁移，b值增加表明从以构造响应为主向受流体压力扰动影响增强的转变。
地震活动集中于利巴尔含水层的原因在于其水文地质特征：该含水层规模更大、水力整合程度更高，且包含由发达管道网络连接的广泛承压和半承压区段。这些特征有利于补给期间水头的积累和压力扰动在区域距离上的有效传递。主要泉点附近的地震集中进一步支持了这一解释，因为这些排泄区是区域地下水流汇聚之处，异常补给期间水力梯度特别陡峭，瞬态孔隙压力变化预计最大。
岩溶含水层对降雨诱发地震特别敏感，归因于三个水文地质特征的组合：补给期间异常大的地下水位上升（高补给率和低有效孔隙度导致水位上升幅度比其他岩性高4至17倍）、通过发达管道网络的快速水文响应，以及通过高度连通管道系统的高效压力传递。
有效水力扩散率在极端补给事件中并非恒定属性，而是随含水层充填过程动态演化。地下水位上升增加了饱和厚度，扩展了承压和半承压区段，增强了管道网络的水力连通性，共同大幅提高了压力传递效率。这可能存在一个水力阈值，超过该阈值后压力传递显著增强，断层再激活概率大幅提高。
研究提出了一种概念模型，将极端补给与地震激活联系起来：持续强降雨在高度岩溶化的补给区产生异常大的地下水补给；由于有效储水能力低，地下水位迅速上升，逐步增强水力连通性；承压和半承压条件的扩展提高了压力传递效率；压力扰动通过管道网络快速传播至周围裂隙岩体；在接近力学失稳的断层上，瞬态孔隙压力增加降低了有效正应力，最终触发破裂。
该研究对自然灾害评估具有重要启示。2026年2月的水文-地震危机表明，极端降水可通过连接大气强迫、地下水补给和断层再激活的过程，在高度岩溶化地区产生次生地质灾害。研究强调，部分承压岩溶含水层对极端补给的响应无法用静态水文地质模型充分描述，地下水系统可能经历瞬态水力重组。将降水、地下水水位、泉流量和地震活动的综合监测纳入多灾种评估框架，对于 tectonically active 岩溶地区至关重要。
研究结论指出，仅凭异常降雨不足以引起与含水层相关的地震活动，含水层还必须具备允许大量补给的水文地质特征，从而产生快速且显著的水位上升。岩溶含水层因其高补给率和低有效孔隙度，水位上升幅度比其他岩性高4至17倍，特别容易产生地震扰动。并非所有接受相似补给的岩溶含水层都发生了相关地震活动，地震仅发生在具有较高水力扩散率的含水层中，如利巴尔含水层。透水系数和水力扩散率并非恒定参数，而是在异常降雨期间随含水层逐步充填而变化。高水位条件下利巴尔山脉周围承压面积增加约57 km2，改善了与88 km2无压补给区的水力联系，增强了压力传递。在最大饱和和承压条件下，主导机制是通过岩溶化区域内的管道直接传递水压力。浅源震中在承压含水层内分布更密集，与含水层及其半承压区域的边界密切吻合，该区域是岩溶含水层主要排泄区所在，水压力对含水层顶板的液压作用可能在促进上覆盖层最薄弱、最破碎区段的破坏中发挥关键作用。
