《Journal of Hydrology: Regional Studies》：Climate-induced changes in flow seasonality and peak discharge with socio-hydrologic consequences in the Hunza River Basin, Upper Indus
编辑推荐：
在喀喇昆仑山脉冰川补给的洪扎流域（Hunza River Basin, HRB），传统灌溉渠系（kuls）因洪水受损的频率日益增加，与此同时农户将作物种植日历提前了数周。气候变暖对该流域的威胁究竟是源于水资源短缺还是时间错配危机？为解答这一问题，研究人员将一个
在喀喇昆仑山脉冰川补给的洪扎流域（Hunza River Basin, HRB），传统灌溉渠系（kuls）因洪水受损的频率日益增加，与此同时农户将作物种植日历提前了数周。气候变暖对该流域的威胁究竟是源于水资源短缺还是时间错配危机？为解答这一问题，研究人员将一个经过率定的逐日水文模型（UBCWM）与覆盖20个村庄共491户家庭的分层抽样调查相结合。该模型基于1995–2015年Danyore站逐日径流数据完成率定，并在2016–2023年时段进行独立验证（NSE = 0.85，KGE = 0.89），基线期内冰川融水贡献了年径流量的70%。基于包含SSP1–2.6、SSP2–4.5和SSP5–8.5情景的五成员CMIP6集合模拟至2100年，水文过程线发生了显著的重组织：洪峰提前12天，晚冬流量急剧上升，高流量日（Q80 = 609 m3 s-1）和极端流量日（Q90 = 895 m3 s-1）每年中位数增加9天，其中5月至6月的信号约在2038年从内部变率中显现出来。尽管春季总径流量增加，但86%的家庭报告春季缺水状况恶化，原因是更早、更集中的融水脉冲在主要作物需水期之前到来，且其强度足以破坏负责配水的引水口设施。与此一致的是，75%的农户将灌溉起始时间提前了2周或以上，67%的家庭报告过去5年内遭受过洪水或滑坡破坏。模拟得到的物理性提前与社会调查揭示的运行压力实际上是同一转变的两个侧面：水流更早、更强烈地到来，超出了既有基础设施的设计承受能力。该流域的水安全关键在于径流时序与渠系韧性，而非绝对的水量短缺，因此研究人员建议优先采取引水口防护和渠系加固措施，而非建设水库规模的供水工程。
研究背景与意义
洪扎河流域位于巴基斯坦吉尔吉特-巴尔蒂斯坦地区的上印度河上游，是兴都库什-喀喇昆仑-喜马拉雅（HKH）这一亚洲大陆"水塔"的关键组成部分。该流域约28%的面积被冰川覆盖，冰川融水贡献了基准期（1995–2023）年径流量的70%，使其成为研究气候变暖下水文响应的理想早期指示器。长期以来，喀喇昆仑地区因冰川质量相对稳定甚至略有增加的"喀喇昆仑异常"（Karakoram Anomaly）现象而备受关注，但自2010年代中期以来观测证据表明该异常正在减弱，区域正转向净物质亏损状态。在此背景下，传统灌溉渠系（kuls）面临日益严峻的挑战：一方面农户为适应气候将作物种植日历提前，另一方面渠系引水口频繁遭受融雪洪水和泥石流破坏。现有研究普遍存在物理维度与社会维度割裂、水文模型对冰川径流组分划分简化、缺乏面向基础设施运营的阈值指标等三方面不足。针对上述问题，研究人员采用社会水文耦合的研究范式，将高程分带水文模拟与分层入户调查相结合，旨在量化气候变化导致的径流时序重组特征，并评估其对社区灌溉管理和基础设施韧性的实际影响，研究成果发表于《Journal of Hydrology: Regional Studies》。
关键技术方法
研究人员构建了基于高程分带（hypsometric）的UBC流域水文模型（UBCWM），将流域划分为8个海拔带，显式模拟积雪融水、冰川融水、基流和降雨径流四个产流组分。模型输入包括ASTER GDEM v3数字高程模型、ICIMOD土地覆盖数据、Randolph冰川编目v7.0、PMD和WAPDA站点实测气象数据及ERA5再分析资料。径流基准序列（1995–2023）由WAPDA实测数据与GloFAS-ERA5再分析产品插补拼接而成。未来气候强迫来自5个CMIP6全球环流模式（ACCESS-CM2、BCC-CSM2-MR、CNRM-CM6-1、INM-CM5-0、MPI-ESM1-2-LR）在SSP1–2.6、SSP2–4.5和SSP5–8.5三种排放情景下的逐日输出，经偏差校正后驱动模型。模型率定期为1995–2015年，验证期为2016–2023年，采用纳什效率系数（NSE）、克林-古普塔效率（KGE）、皮尔逊相关系数和百分比偏差（PBIAS）进行评价。同时，研究人员于2023年在沿喀喇昆仑公路分布的20个村庄开展了覆盖491户家庭的分层抽样入户调查，样本按土地持有规模分层以确保小、中、大规模农户的代表性。社会水文耦合分析采用Logistic回归模型，将村庄尺度的模拟暴露指标（峰值提前天数、极端流量日变化、复合暖湿事件频率）与家庭尺度的调查响应（灌溉起始提前、渠系损坏经历、灌溉冲突、维护出资意愿）进行统计关联，标准误经村庄层面聚类稳健调整。
研究结果
3.1 水文模型性能
UBCWM在率定期（1995–2015）和验证期（2016–2023）均表现出满意的性能，率定期NSE为0.74、KGE为0.75，验证期NSE提升至0.85、KGE达0.89，验证期平均模拟径流（305 m3 s-1）与实测值（308 m3 s-1）几乎重合。季节性评估显示夏季（JJA）和春季（MAM）模拟效果最佳，冬季（DJF）基流因模型缺乏冻土模块而被低估约40%，但该偏差不影响基于高流量尾端（Q80、Q90）的灾害指标计算。值得注意的是，验证期性能优于率定期，这与2016–2023年间喀喇昆仑异常活跃年份（2019、2020、2022）导致积雪融水占比升高有关，间接验证了模型对积雪融水过程的良好刻画能力。
3.1.1 年尺度径流组分划分
基准期内冰川融水平均贡献年径流的70%（年际变幅10%–92%），积雪融水占18%，基流占11%，直接降雨径流仅占2%。冰川融水份额最低值集中在2019–2022年积雪丰沛年份，印证了近期气候态的转变。
3.1.2 月尺度径流组分划分
率定期呈现典型冰雪融水型（nival-glacial）格局：12月至次年3月基流主导（67%–98%）；4月冰川融水启动（占60%）；5–6月积雪与冰川融水共同主导（积雪31%–32%，冰川61%–64%）；7–9月冰川融水占绝对优势（80%–87%）；10–11月过渡回基流主导。验证期5–6月积雪融水占比升至39%–45%，冰川融水降至42%–44%，反映了近期气候态下春季融雪提前增强的特征。
3.1.3 季节尺度径流组分划分
夏季半年（4–9月）冰川融水贡献72%的总径流，积雪融水占20%；冬季半年（10–3月）基流占71%，冰川融水仍贡献27%。季度分解显示：12月至次年2月（DJF）基流占93%；3–5月（MAM）冰川融水领先（54%）并伴随显著积雪融水贡献（26%）；6–8月（JJA）冰川融水占72%；9–11月（SON）冰川融水占70%。
3.2 2015–2100年水文变化预测
所有SSP情景下流域均保持单峰冰雪融水型水文情势，但峰值量级与发生时间显著改变。世纪中期（2031–2060）7月平均径流量较基准期增幅达106%–157%（SSP1–2.6为+106%，SSP5–8.5为+157%），晚冬（2–3月）平均流量相对增幅达196%–377%，尽管绝对值仍远小于夏季峰值。年径流序列显示SSP1–2.6情景下径流量在2030年左右达到平台期（较基准期+132%），SSP2–4.5下逐步增至世纪末的+175%，SSP5–8.5下持续攀升至+293%。峰值出现时间在三个SSP情景下均系统性提前：世纪中期中位数提前12天（四分位距9–15天，由基准期DOY 220提前至DOY 208），SSP5–8.5情景下世纪末进一步提前至DOY 195（提前25天，四分位距18–32天）。峰值提前是所有情景下情景间重叠度最高的信号，是世纪中期最具鲁棒性的水文响应特征。
3.3 灾害相关指标
基准期Q80 = 609 m3 s-1，Q90 = 895 m3 s-1，对应高流量日（超过Q80）中位数73天/年，极端流量日（超过Q90）中位数36天/年。世纪中期两类高流量日中位数均增加9天/年（四分位距6–12天）。若基于2031–2060年非平稳基准重新计算阈值，高流量日增量降至7天/年，表明约四分之三的信号源于流量分布的真实上移而非阈值选取的人为效应。5–6月窗口对融雪启动加速最为敏感，该时段超过Q90的天数中位数增加11%（四分位距8%–14%）。5–6月Q90信号的时间显现点（Time of Emergence, ToE）为2038年（四分位距2035–2042），是规划引水口防护和春季灌溉基础设施的最关键操作指标。敏感性分析表明上述结果对偏差校正方法选择（CF vs QDM）、夏季流量±15%扰动以及重雨阈值（P90–P99）设定均具有鲁棒性。
3.4 喀喇昆仑公路沿线调查证据
491户家庭的调查结果与水文信号高度吻合。86%的受访者报告春季缺水状况恶化（56%显著增加，30%增加），这与模拟显示的晚冬至早春流量上升形成鲜明对比——原因在于更早、更集中的融水脉冲在主要作物需水期之前到来，且其强度破坏了kul引水口设施，导致总流量增加并未转化为可用水量。灌溉日历提前现象普遍：45%的农户将灌溉起始提前超过4周，30%提前2–4周。67%的家庭过去5年内经历过洪水或滑坡对基础设施或田地的破坏，其中19%经历3次，6%经历4次及以上。灌溉冲突整体较少（59%无冲突），但发生的冲突集中在5–6月和7–8月，与模拟的5–6月流量异常时段吻合。洪水与冰湖溃决洪水（GLOF）早期预警系统公众认知度低：41%完全不知晓任何预警系统，26%表示毫无信心，仅15%表示中等或高度信任。适应措施以用水轮灌（25%）、非农收入（22%）、改种作物（15%）和渠道衬砌/蓄水（12%）为主，27%未采取任何适应措施。维护出资意愿呈小额分散特征：38%不愿出资，29%愿出资不超过2000卢比，35%愿出资超过2000卢比。人口流动性显著：65%为临时迁移，16%季节性迁移，6%永久迁移。对未来展望方面，33%预期2035–2050年夏季可用水量大幅减少，30%预期有所减少；对春季缺水预期则分化明显（38%预期显著缓解，33%预期中度缓解），反映出社区对长期冰川消退的认知与对近期融雪提前带来晚冬流量上升的预期并存。
3.5 社会水文耦合
Logistic回归分析定量证实了模拟暴露指标与调查响应之间的统计关联。灌溉起始提前（Y1）与模拟峰值提前（E1）呈显著正相关（比值比OR = 1.45，p < 0.01）；kul/引水口损坏经历（Y2）与极端流量日变化（E2）关联最强（OR = 1.62，p < 0.01），同时也受峰值提前和复合暖湿事件的影响；晚季灌溉冲突（Y3）与峰值提前（OR = 1.38，p = 0.02）和极端流量日变化（OR = 1.41，p = 0.015）均显著相关；维护出资意愿（Y4）与极端流量日变化呈弱正相关（OR = 1.27，p = 0.072），表明家庭财务响应还受到模型未捕捉的其他因素驱动。最强的两个关联——峰值提前与灌溉起始提前、极端流量日增加与kul损坏——直接指向引水口防护、渠系加固和更新维护计划的核心适应优先级。
讨论
本研究的核心发现在于：尽管春季总径流量增加，但水安全的主要制约并非绝对水量短缺，而是径流时序错配与基础设施输配能力不足。物理维度的水文过程线重组织与社会维度的运营压力是同一气候驱动转变的两个表征。相比以往研究，本研究的创新在于将气候信号转化为可操作的时间特异性阈值指标（12天峰值提前、2038年时间显现点、9天额外极端流量日），并通过社会水文耦合提供了从流域尺度到家庭尺度的完整证据链。模型验证期在喀喇昆仑异常活跃年份表现优异，增强了结果的可信度。主要不确定性来源于UBCWM采用静态冰川面积和度日因子法，无法动态反映冰川物质平衡耗尽后的"峰值水"（peak water）拐点，因此世纪末径流量应视为上限估计；冬季基流低估因缺乏冻土模块所致，但不影响高流量灾害指标；5个GCM的集合规模处于CORDEX推荐范围下限，但信号在各情景间高度一致，情景间差异主要体现在2050年后。基于上述证据，研究人员建议将传统的基于年度总量的水资源规划转变为基于时序阶段的运营管理：制定晚冬准备期（2–3月）、初夏承压期（5–6月）和主汛期（7–8月）的分阶段操作规程；按新增9天极端流量日和12天峰值提前量设计引水口防护和渠系加固标准；利用社区小额出资意愿建立规则化维护机制，同时争取外部资金配套支持大型工程；扩大社区级洪水与GLOF预警系统覆盖范围以弥补现有空白。
结论
本研究通过耦合UBC流域水文模型（由5个CMIP6模式在三种SSP情景下驱动）与20个村庄491户家庭的分层抽样调查，评估了1995–2100年洪扎河流域径流季节性与洪峰流量的气候变化响应。基准期内流域高度依赖冰冻圈：冰川融水占Danyore站年径流量的70%（年际变幅10%–92%），积雪融水占18%，基流占11%，直接降雨径流仅占2%。模型率定期NSE为0.74，独立验证期NSE达0.85。预测结果显示水文过程线发生一致性重组织：洪峰提前约12天，晚冬流量相对增幅达196%–377%，世纪中期极端流量日每年增加9天，5–6月Q90异常信号约在2038年从内部变率中显现。物理信号与家庭体验相互印证：86%的家庭报告春季缺水恶化，75%提前灌溉起始2周以上，67%过去5年内遭受洪水或滑坡破坏。Logistic回归证实峰值提前与灌溉起始提前（OR = 1.45）、极端流量日增加与kul损坏（OR = 1.62）之间存在统计学显著关联。径流时序提前与基础设施输配瓶颈而非绝对水量短缺构成了该流域水安全的核心风险。基于12天峰值提前量和9天额外极端流量日的低后悔干预措施包括引水口防护、渠系加固、分散式蓄水设施建设以及社区预警系统完善。将基于年度总量的规划转变为基于时序阶段的运营管理，是将本研究结果转化为实际适应的关键路径。
