旧金山湾区(California, USA)大气河流风暴捕获的极端地下水文条件评估及其在浅层滑坡早期预警中的应用

研究人员评估了2022—2023年冬季影响加利福尼亚州旧金山湾区(San Francisco Bay Area, SFBA)的一系列大气河流(Atmospheric River, AR)风暴期间所记录的极端次表面水文条件，并将其应用于浅层滑坡(shallow landslide)早期预警(early warning)。通过分布于东湾(East Bay)、马林县(Marin County)及旧金山半岛(San Francisco Peninsula)的三个USGS BALT监测站点，获取降雨、体积土壤含水量(换算为有效饱和度Se)及基岩界面上方孔隙水压力(pore water pressure / piezometric head)，并结合土工反分析(geotechnical back-analysis)与时空标记滑坡编目，探讨导致广泛浅层滑坡发生的因素。研究发现，当有效饱和度Se> ~65%时孔隙水压力呈线性增长，广泛滑坡触发需Se= 100%且持续数小时并伴有高正孔隙水压力(>40 cm水头)；仅依据AR强度或静态峰值孔隙水压力不足以判定区域滑坡规模，需结合前期土壤饱和度和饱和持续时间。研究成果支持将原位土壤湿度监测整合入基于降雨强度—历时(intensity-duration, I-D)阈值的滑坡预警体系，并为考虑三维坡面流场中持续孔隙水压迁移的动态视角提供依据。

论文解读：旧金山湾区大气河流风暴极端地下水文条件评估及其在浅层滑坡早期预警中的应用

本文发表于《Landslides》，研究对象为2022年12月26日至2023年1月16日袭击美国加州旧金山湾区(San Francisco Bay Area, SFBA)的八场连续大气河流(Atmospheric River, AR)风暴所引发的广泛浅层滑坡(shallow landslide/debris slide-flow)。既往研究显示仅靠降雨强度—历时(Intensity-Duration, I-D)阈值进行滑坡预警存在局限，Keefer等(1987)早先提出需引入前期土壤饱和度(antecedent soil moisture)、孔隙水压力(pore water pressure)及地质—水文控制因素，但该类建议在加州长期未系统实施。2022—2023年AR风暴序列虽AR等级不一，却在第二场中等强度AR(AR2—3)后引发东湾(BALT1)及半岛(BALT3)区域多达数百处浅层滑坡，而更强AR4首场风暴因前期土体未饱和几乎无滑坡，凸显单纯AR等级或雨量阈值之不足。本文依托USGS预先布设之BALT监测网(代表区内典型易滑凹形—平面残积土坡面)，结合土工参数测试、无限斜坡稳定性分析(infinite slope stability analysis)及地面—遥感滑坡编目，阐明极端饱和条件与孔隙水压动态对区域浅层滑坡触发之控制作用，并讨论其对改进浅层滑坡早期预警(early warning system, EWS)之意义。

主要关键技术方法简述：

研究人员利用旧金山湾区三个BALT监测站(BALT1东湾、BALT2马林县、BALT3旧金山半岛)近实时体积土壤含水量传感器(Meter Group ECH₂O EC-5)及敞开式渗压计(vented piezometer, Greenspan PS7000/Stevens SDX)获取次表面水文数据；以瞬态释放与吸渗法(Transient Release and Imbibition Method, TRIM)测定土—水特征曲线(soil-water characteristic curve, SWCC)参数(θ_s、θ_r、k_s)；以ASTM D3080固结排水直剪试验获Mohr-Coulomb有效强度参数(c′、φ′)；依实验室θ_s、θ_r按Mualem公式算有效饱和度S_e＝(θ－θ_r)/(θ_s－θ_r)并取浅深传感器均值；以最小8 h无雨界定独立风暴(Inter-Event Time, MIT＝8 h)并比对GFS模式积分水汽输送(Integrated Vapor Transport, IVT)定AR等级(Ralph AR Scale)；以地面踏勘、航拍、卫星影像(GIS lidar坡度筛选≥35°易滑坡面)建立各子区(10—20 km²)滑坡编目；采用无限斜坡安全系数公式F_s＝[c′＋(γHcos²β－μ_w)tanφ′]/(γHsinβcosβ)，以实测孔压头h_p计算μ_w＝γ_w·h_p，取BALT1临界坡角β_crit(F_s＝1.0)横向比较各站稳定性。

研究结果

Antecedent conditions and storm meteorological sequence（前期条件与风暴气象序列）：

八场风暴总雨量10—150 mm，AR等级AR2—AR4不等，但滑坡响应不与AR等级成正比。首场AR4风暴S1前有效饱和度S_e仅35—65%，未触发广泛滑坡；紧接S2(AR2—3)前东湾及半岛前期S_e已升至高位，S2致BALT1、BALT3达S_e＝100%并维持5—16 h，伴最高孔隙水压力(东湾>40 cm，半岛>60 cm)，同期发生最多滑坡。表明前期土体饱和度为是否发灾之首要门槛。

Subsurface hydrologic response（次表面水文响应）：

BALT1与BALT3于S_e≈65%起孔隙水压力随S_e线性增加(R²＝0.86)，S_e>90%时出现显著正超静孔压(>40 cm)；BALT2始终未超S_e≈85%且无记录到正孔压(较高比产水率specific yield S_y＝θ_s－θ_r＝0.33、较深土层、可能透水基岩)。广泛滑坡仅发生于S_e＝100%且持续饱和(东湾5 h、半岛16 h)叠加持续高孔压时段，同等峰值孔压若无持续完全饱和则滑坡较少。

Characteristics and spatial distribution of landslides（滑坡特征与空间分布）：

滑坡均为浅层碎屑滑动—碎屑流(shallow debris slide-flow)，独立源区无重复同源失败。BALT1子区(12.4 km²)录得222处(18 landslides/km²，归一化陡坡≥35°后更高)，BALT3约3 landslides/km²，BALT2＜1 landslides/km²。东湾与半岛大量滑坡集中于S2，后续风暴极少新增，马林县全序列少灾。

Analysis of triggering conditions（触发条件分析）：

无限斜坡稳定性计算示BALT1典型1.2 m深35°坡于实测S2孔压下F_s≈1.0可解释广泛失败；BALT2同35°坡F_s≈1.7(无正孔压、高抗剪强度)故稳定；BALT3虽孔压更大(最大60 cm，接近土柱高)但因土体粘聚力c′≈三倍BALT1(粉砂—含黏土 vs 粉砂)，35°坡F_s≈1.3，仅最陡最弱坡段失效，故滑坡密度低于BALT1。差异主因为土层岩性(土体强度)及水力传导度(k_sBALT3≈3.3×10⁻⁴m/s vs BALT1≈2.6×10⁻⁸m/s)导致饱和持时不同。

讨论与结论(D原文Conclusions浓缩翻译)：

研究人员发现孔隙水压力之于浅层滑坡触发非静态单一指标，等效幅值孔压在不同风暴未产生同等滑坡规模，广泛触发需降雨、土体饱和及孔压三者同时达峰且饱和具持续性(S_e＝100%并维持)。S_e＞~65%时土体孔隙连通建立，近基底失效面正孔压线性发育；完全饱和消除基质吸力(suction stress)并使孔压于基底不平整(fill-and-spill bedrock depressions)更广泛传递，扩大潜在破坏面范围从而增加区域滑坡概率。AR等级本身不足以发布滑坡预警——本例最强AR4因低前期S_e未致灾，次强AR2—3在高S_e背景下引发广泛灾害；传统前期雨量阈值(250—760 mm)在本序列未满足仍发生大灾，而原位土壤有效饱和度监测可更准确捕捉前期条件。降雨I-D阈值结合前期S_e之水文—气象阈值(hydro-meteorological threshold, 初值S_evs 24 h最大雨量)可区分高孔压与低孔压风暴，有助于EWS分级预警。尽管地质—水文异质性限制通用阈值推广，持续直接监测次表面饱和与孔压(或经校准水文模型预测饱和历时)是改进浅层滑坡区域早期预警之核心方向。本研究证实广泛浅层滑坡触发属极端水文—土质耦合事件，最可能发生於易滑(弱)坡面当降雨、有效饱和度及孔隙水压力共同达极值且饱和状态得以维持之时。