《Sci》：Concealed Groundwater Day Zero: An Uncertainty–Scarcity–Vulnerability Framework for Assessing Groundwater Unsustainability and Hydro-Social Risk in Bangladesh
编辑推荐：
地下水枯竭是孟加拉国面临的重要环境与社会经济挑战。本研究利用孟加拉国水资源开发委员会（BWDB）400口监测井1990—2019年的观测数据，评估地下水枯竭的时空动态及水—社会脆弱性。概念基础为可靠性—恢复性—脆弱性（RRV）框架；在操作层面，该框架被调整为三
地下水枯竭是孟加拉国面临的重要环境与社会经济挑战。本研究利用孟加拉国水资源开发委员会（BWDB）400口监测井1990—2019年的观测数据，评估地下水枯竭的时空动态及水—社会脆弱性。概念基础为可靠性—恢复性—脆弱性（RRV）框架；在操作层面，该框架被调整为三个归一化的地下水压力分量：不确定性（U，阈值超越频率，因而作为不可靠性的代理指标）、稀缺性（S，干旱事件后的持续性或不恢复程度）与脆弱性（V，地下水水位亏缺的严重程度）。这三个分量被组合为地下水日零指数（DZI），同时人口与灌溉图层被单独用于刻画暴露度。当标准化地下水水位低于长期均值的?1个标准差时，即界定为干旱事件。结果表明地下水枯竭广泛存在，在受压的西北部与中部地区报告的速率约为0.15—0.21 m yr?1。达卡具有严重的DZI（约0.91；>0.75），而Barind Tract的大部分区域处于高至严重状态（约0.72—0.78）。全国分量均值（U = 0.46，S = 0.52，V = 0.57）表明总体以中至高度压力为主。DZI与灌溉/城市强度（r = 0.74，p < 0.01）、人口密度（r = 0.76，p < 0.01）及灌溉强度（r = 0.81，p < 0.01）之间的空间关联提供了趋同证据，表明该指数能够识别具有显著人为地下水压力的区域。DZI旨在作为一种诊断性筛查指标，而非对地下水失效日历日期的预测。其主要价值在于将地下水压力的频率、持续性与严重程度整合于一体，以支持空间优先级排序，并随新的监测与气候信息的获取而定期更新。
研究背景方面，地下水是孟加拉国饮用水供应与旱季农业特别是Boro稻作的核心水源，灌溉、市政与工业的集约开采已导致西北部与中部地区地下水持续下降。气候变率通过改变补给与灌溉需求进一步调节这一压力：气温升高增加作物需水，季风降水的时序与强度决定季节性补给机会。因此该国地下水压力是开采、补给变率、土地利用变化与人口增长交互作用的结果，而非单一驱动因素所致。已有研究借助监测井、GRACE（Gravity Recovery and Climate Experiment）、GLDAS、InSAR（Interferometric Synthetic Aperture Radar）、数值模型、GIS与遥感等手段，记录了水位下降、储量损失、补给过程、地下水潜力与多灾种脆弱性。这些方法回答的是不同问题：趋势分析量化下降方向与幅度，潜力或脆弱性制图刻画空间易感性，卫星产品量化区域储量变化，而地下水干旱研究通常将干旱定义为水位或储量低于阈值并持续一定时长。现有方法学缺口并非缺少地下水风险研究，而是缺少一个一致的国家尺度框架，能够利用长期井观测描述阈值压力发生的频率、恢复失败的持续性以及亏缺的严重程度，并将该物理压力与人口和灌溉暴露相联系。
针对这一较窄的缺口，研究人员将成熟的可靠性—恢复性—脆弱性（RRV）概念适配于地下水水位时间序列。USV操作化表述应被理解为对RRV的适配而非独立竞争理论：U代表压力频率（可靠性的补集），S代表持续性/不恢复（恢复性的倒数），V代表亏缺严重程度。三者整合为日零指数（DZI），保留了趋势系数无法捕捉的失效频率、持续性与量级信息。与常规地下水潜力或易感性图不同，DZI基于观测地下水水位的时间行为；与基于需求—可利用量比值的常规水压力指数不同，它诊断含水层系统自身的响应。“日零”概念在水安全文献中被用于描述严重稀缺下常规供水预期失效的情形，在本研究中该术语仅作为地下水系统功能渐进丧失的隐喻，DZI不代表日历失效日期或水文地质临界点阈值，而是代表相对地下水压力与干预优先级。研究设定了四个目标：利用国家监测网量化长期与年代际地下水枯竭格局；构建RRV衍生的USV表述与由水位干旱频率、持续性和严重程度组成的复合DZI；识别地下水压力空间热点并考察其与灌溉和人口压力的对应关系；将由此形成的水—社会格局转化为区域特定管理启示，同时明确识别框架的数据与方法学局限。其新颖性在于从同一国家井数据集整合地下水压力的三个时间属性——频率、持续性与严重程度，并与暴露图层比较；该框架旨在补充而非替代已有的枯竭趋势、地下水潜力图、水压力指标与过程型地下水模型。
研究人员开展的研究与结论：基于400口井1990—2019年月度观测，结果显示全国范围内地下水普遍加深，中部城市带与西北灌溉带持续下降最为显著。以达卡DH005、DH007、DH010为例，水位由早期约2—8 mbgl（meters below ground level）变为2018年以8—14 mbgl为主，旱季深度更大、季节变幅更宽；拉杰沙希RJ031、RJ081、RJ144由约4.5—7.2 mbgl加深至10.8—14.6 mbgl，三口井平均变化约6.5 ± 1.2 m，季节变幅由约1.5—2.0 m增至3—4 m，RJ144最大深度超过15 mbgl。年代际分布显示达卡三口井典型深度由1995—1999年约2—6 mbgl移至2010—2019年约7—12 mbgl，拉杰沙希三口井由约4—7 mbgl移至10—15 mbgl，月变率由约1.5—2.0 m增至3.5—4.5 m。三分量空间格局上，不确定性高值（>0.75）集中于达卡、Naogaon、Joypurhat与Sylhet，合计约占国土27%；稀缺性超过0.70的区域主导中部与东北部，包括达卡、Moulvibazar、Habiganj与Rangamati，合计约占31%；脆弱性超过0.80的区域为达卡、Sunamganj与拉杰沙希。复合DZI范围为约0.15（低）至0.93（严重），约22%国土属高值类（>0.50—0.75）、8%属严重类（>0.75）；达卡DZI约0.91，西北部多为约0.65—0.78。暴露图层显示西北与中北部灌溉覆盖指数超过88—96%，高密度人口（>6000 persons km?2）集中于达卡、Gazipur、Narayanganj、Chattogram及Mymensingh部分地区；约28%人口即近4700万人处于高至严重风险区，达卡人口加权DZI为0.87；约34%灌溉土地处于高至严重风险区。全国21个县属高或严重脆弱类，约占国土35%与人口40%，形成达卡与Chattogram的城市压力、拉杰沙希—Naogaon—Dinajpur的灌溉驱动枯竭、以及Habiganj—Sylhet—Sunamganj的气候脆弱区三大风险簇。其重要意义在于DZI将频率、持续性与严重程度整合为可用于国家筛查与优先级排序的诊断工具，并强调其应作为分级评估的第一阶段而非场地特定管理的替代品。
关键技术方法（不超过250字）：研究采用孟加拉国水资源开发委员会（BWDB）全国监测网400口井1990—2019年月度地下水埋深数据，经统一为米地下水埋深（mbgl）、去重、单位与时序核查、异常值标记等质量控制，缺失约6—10%，短内部缺口用三次样条插值且不外推。以标准化水位低于?1个标准差界定干旱事件，按RRV框架计算不确定性（U）、稀缺性（S）、脆弱性（V）三分量并等权合成日零指数（DZI）。空间插值采用反距离权重法（IDW）与普通克里金法（OK）并经交叉验证；分量经最小—最大归一化。暴露分析使用BARC行政边界、BBS 2022人口普查并与WorldPop交叉核对的人口密度、以及FAO、MODIS MCD12Q1与国家土地利用数据库的灌溉/耕地数据，统一至WGS 84与1 km×1 km网格叠加，并以Moran's I评估空间自相关。
研究结果：5.1 地下水枯竭——通过400井网与六口代表性井的水位分布对比，得出全国普遍加深、达卡与拉杰沙希显著下移且季节变幅扩大的结论。5.2 年代际地下水分布动态——通过1995—1999、2000—2009、2010—2019三期月度水位概率分布，得出典型深度持续加深、旱季分布展宽的结论。5.3 不确定性、稀缺性与脆弱性空间评估——通过三分量归一化制图，得出高U集中于达卡、Naogaon、Joypurhat、Sylhet，高S主导中部与东北部，高V集中于达卡、Sunamganj与拉杰沙希的结论。5.4 复合地下水日零指数——通过U、S、V等权合成，得出DZI约0.15—0.93、约22%高值、8%严重、达卡约0.91的结论。5.5 灌溉覆盖与人口密度空间格局——通过灌溉与人口图层制图，得出西北与中北部灌溉覆盖极高、城市群人口密度极高、二者叠加形成复合压力的结论。5.6 人口与灌溉风险制图——通过DZI与暴露图层叠加，得出约28%人口与约34%灌溉土地处于高至严重风险区，且DZI与灌溉/城市强度（r = 0.74）、人口密度（r = 0.76）、灌溉强度（r = 0.81）显著相关的结论。5.7 地下水资源综合总风险——通过水文不可持续性与社会经济暴露综合，得出达卡综合风险>0.85、拉杰沙希0.74、Naogaon 0.71、Habiganj与Sunamganj 0.65—0.69，以及21县高或严重、三大风险簇的结论。
讨论部分总结：全国格局与既有野外、卫星与水文地质评估一致，DZI不替代这些估计，而是区分压力发生的频率、恢复的持续性与亏缺的严重程度，因为相同长期趋势可能来自反复季节压力与长期多季干旱的不同组合。水文气候与开采在各区域交互方式不同：Barind Tract旱季灌溉需求与较低降水、厚层黏土质表层及有限补给机会相叠加，达卡则是高开采与城市地表下渗减少相叠加；年降水量本身是不完整的地下水响应预测因子，季节时序、事件强度、前期水位、地表渗透性与抽水共同决定降水能否形成有效补给。DZI与主要比较需求或取水与可再生水量的常规水压力指数（WSI）不同，前者反映地下水系统对累积水文气候与开采压力的响应，二者互补但仅在空间单元、时段与水量平衡定义一致时才宜直接比较。研究未识别1990—2019年下降的单一因果因素，因全国一致的抽水、人均消费、井结构与补给时间序列在400个监测点不可得；与人口密度和灌溉强度的相关性应解释为关联而非因果归因。水质是DZI计算之外的重要后果，水位下降可改变垂向水力梯度并可能增加咸水运移、砷再分配或抽取更深层水的风险，处理需求具有场地特异性。2019年后状况未纳入指数，近期气候变率可能双向改变DZI，框架应随新监测数据定期重算而非视为静态风险图。与趋势制图、GIS—AHP（Analytic Hierarchy Process）易感性评估、遥感储量分析与WSI类指数相比，DZI保留长期时间压力动态并分离频率、持续性与严重程度，但经验性?1个标准差阈值、样本依赖的归一化、未校准的等权、空间插值与有限的暴露数据约束其解释，故最宜与过程型地下水与水量平衡评估并用。可移植性要求长期质控地下水记录、一致时间分辨率、适当压力阈值、暴露数据与可比归一化，阈值与恢复定义应按当地水文地质条件验证。空间关系提供趋同而非预测性验证，未来需与趋势分类、WSI、卫星储量异常与观测到的服务失效正式比较。管理启示按热点类型区分：城市中心优先需求管理、漏损削减、计量与开采控制、补给区保护、可行时增加处理地表水利用，以及在水文地质与水质适宜性评估后实施管理含水层补给；灌溉依赖区优先高效灌溉、在本地可行处减少最耗水旱季生产的作物多样化、地表水与地下水联合运用、基于地下水状况的调度以及补给通道的保护与修复。实施应依托现有主管机构，BWDB可用于定期地下水状况更新，农业与灌溉部门可据高DZI灌溉区优先安排效率与作物管理干预，城市水务企业可据持续高DZI区指导开采控制与水源多元化；业务化采用前应以最新监测数据重算并与本地抽水、井性能与需水记录评估。结果支持分级决策：先由DZI识别优先地点，再由详细水文地质、需求、水质与社会经济调查确定适当干预，从而避免将复合指数当作场地特定地下水管理的替代。
研究结论翻译：本研究构建了RRV衍生的地下水日零指数（DZI），用于孟加拉国地下水压力的国家筛查。基于400井网月度观测，框架将阈值压力频率、持续性/不恢复与亏缺严重程度概括为U、S、V分量，并将由此形成的空间格局与人口和灌溉暴露相联系。最高DZI值出现在达卡与西北灌溉带，与城市和农业地下水压力强烈的独立证据一致。DZI应被解释为地下水系统性能下降的诊断指标，而非含水层枯竭日期的预测。主要局限为1990—2019年监测时段、全国一致的抽水与人均用水数据不完整、井结构与水质信息有限、经验性?1个标准差阈值、样本依赖的归一化与等权分量。未来发展应包括：以2019年后观测更新分析；对全部监测网按井分类下降速率；检验替代阈值与分量权重；将DZI与WSI、趋势型、卫星储量与过程型地下水指标直接比较；以观测到的井失效、抽水限制与服务中断进行验证；整合气候调整的补给、开采、水质与水—能源—粮食反馈。管理上该指数最宜作为分级评估的第一阶段：持续高DZI的城市区需开采控制、漏损削减、水源多元化与补给保护，高DZI农业区需灌溉效率提升、联合运用、作物多样化与基于地下水状况的抽水管理；业务化使用应与BWDB及相关水、农业与城市服务机构共同开展，并随新监测数据定期更新。
