观潮 · 模式预报与气象 AI 观察

湘河水库流域 46 年降水:7—8 月占全年 60.70%,10、20、30 年一遇最大 1 日降水 45.3、57.7、68.5 mm

2026-10-23

> 湘河水库是国务院 172 项重大节水供水工程之一,流域海拔 4049—6126 m,2020 年之前只有县政府所在地和少数重点区域有气象站,乡镇级长期没有连续观测。作者用南木林国家基本站 1980—2025 年的 46 年长序列,配上上游 5 站 2021—2025 年的短序列和 ERA5 再分析,先把「一个站能不能代表全流域」这件事验了,再算集中度、周期、极端降水重现期和丰枯环流差异。结果是 7—8 月占全年 60.70%、主周期 23 a、10 年一遇最大 1 日降水 45.3 mm。作者明确列出三件没做的事:未按季节和等级拆分集中度、只用了单一 P-Ⅲ 型分布、驱动机制只做了定性分析没量化贡献度。

缺的是哪一块

降水时序特征的分析方法本身已经很成熟——气候倾向率、滑动平均、降水集中度(RCD)与集中期(RCP)、Mann-Kendall 突变检验、小波分析、百分位法配重现期,这套组合在雅鲁藏布江流域、三江源、四川盆地、云贵高原都验证过。作者指出缺的不是方法,是对象:高原新建大型水利枢纽所在的流域,尤其是以农村水利和生态保护为主功能的高原水库流域,拿长期与短期观测互补起来做系统降水研究的很少,湘河水库流域正是这类典型。

数据条件也是硬约束:2020 年前流域内只在县政府所在地和部分重点区域设站,乡(镇)级长期缺连续观测;2020 年起乡镇全面建自动站实现「乡乡有站」,才补上了这块短板。所以本研究的结构是 46 a 长序列(南木林,距水库 20 km)与 5 a 短序列(上游仁堆、普当、甲措、拉布普、热当)互补。

方法怎么搭的

第一步先解决代表性问题。用 Pearson 相关检验上下游同步性:除拉布普站外,上游仁堆、热当、普当、甲措相互间相关系数介于 0.820—0.984(P<0.01);下游南木林站与仁堆站达 0.940(P<0.01),与热当、普当、甲措分别为 0.667、0.708、0.632(P<0.05)。再用 2021—2025 年汛期(5—9 月)逐月数据复核,南木林与甲措 9 月相关 0.975(P<0.01)、与普当 8 月 0.962(P<0.01)、与热当 8 月 0.919(P<0.05)。结论是南木林站可以代表流域气候背景。

集中度用矢量合成原理:把 12 个月映射为单位圆上的等分角度,以月降水量为模长逐月累加,合成矢量模长占年总量的比例即 RCD,方位角即 RCP。面雨量用泰森多边形定权重后加权平均,各站权重为南木林 0.0344、拉布普 0.4782、普当 0.1050、热当 0.0789、仁堆 0.2322、甲措 0.0713。空间插值用反距离权重法(IDW),分辨率 1 km×1 km,交叉验证误差控制在 5% 以内。

极端降水用皮尔逊Ⅲ型(P-Ⅲ)分布:取 46 年汛期逐日资料里的年最大 1 d 降水量做频率分析,重现期 T 与频率 P 满足 T=1/P。环流分型则以 23 a 主周期小波实部序列的上下 25% 分位数为阈值切丰枯年——Q0.25 = −132.22 判强偏枯年,Q0.75 = 102.10 判强偏丰年,中间态不参与合成,最终筛出强偏丰 12 a(1987—1989、2001—2004、2016—2020)、强偏枯 9 a(1994—1996、2009—2012、2024—2025),再做 500 hPa 位势高度距平、风场与水汽通量的合成场对比。

结果数字

年内分配是标准单峰型,1—4 月少,5 月起增,7 月达峰值 145.60 mm,12 月最低:

项 · 数值

6—9 月降水占全年 · 76.04%

7 月 / 8 月占比 · 30.68% / 30.02%(合计 60.70%)

RCD 多年均值 · 0.827(标准差 0.039)

RCP 多年均值 · 221.67°(标准差 8.26°,对应 8 月中下旬)

7 月 / 8 月降水标准差 · 64.58 mm / 52.23 mm

趋势与周期:1980—2025 年汛期降水量增幅 13.35 mm·(10a)⁻¹,UF 统计量大部分年份未突破 ±1.96 的 0.05 显著性线,属不显著微弱上升,且无明显突变点。7 月、8 月增速较大,分别为 9.60 与 5.11 mm·(10a)⁻¹;9 月反而是减少趋势 −3.02 mm·(10a)⁻¹,其占全年比例的变化达 −2.02%·(10a)⁻¹(P<0.01)。周期上存在 5—8 a 与 15—25 a 的显著振荡,其中 23 a 周期的小波方差峰值最突出,2019 年后转入偏少期且负值持续加深。

面雨量分级与极值:

项 · 数值

强降雨阈值 · 日面雨量 > 14.5 mm

中等强度 · 6.2—14.5 mm

弱降雨 · < 6.2 mm

10 a 一遇年最大 1 d 降水 · 45.3 mm

20 a 一遇 · 57.7 mm

30 a 一遇 · 68.5 mm

空间上,2021—2025 年汛期降水量介于 173.1—479.2 mm,高值区在热当、仁堆一带(354.6—414.4 mm),低值区在流域中部与西部(173.1—325.8 mm),拉布普一带多为 208.8—250.5 mm,是流域内的相对低值中心。年际上 2023 年全域偏多(467.7—479.2 mm),2024—2025 年呈以库区为界的「东多西少」格局。

丰枯环流对比差别很清晰:丰水年流域上空受偏弱负位势高度距平(−20—0 gpm)控制,中纬度西风槽、高原低值系统与切变线活跃,水汽通量散度距平 −1.9×10⁻⁵ kg·m⁻²·s⁻¹,处于明显辐合区;枯水年全域被强正距平中心(40—60 gpm)控制,高压脊与反气旋主导,气流下沉辐散,水汽通量散度距平 1.5×10⁻⁵ kg·m⁻²·s⁻¹。

原文图

月降水量与集中度
月降水量与集中度

两个面板都用 1980—2025 年的月尺度序列。(a) 蓝线为各月降水量、黑虚线为标准差、红柱为气候倾向率:7 月降水量峰值 145.60 mm,7 月与 8 月标准差分别为 64.58 mm 与 52.25 mm,9 月的倾向率为负值(正文给的 −3.02 mm·(10a)⁻¹)。(b) 绿线为月降水占全年比例(MPAP),7 月 30.68%、8 月 30.02%,两者合计 60.70%,9 月那一根橙色大柱是 MPAP 倾向率 −2.02%·(10a)⁻¹ 的显著减少信号。

汛期降水空间格局
汛期降水空间格局

2021—2025 年汛期(5—9 月)降水量的逐年空间分布,色标 173.10—513.90 mm,站点标注为仁堆、普当、甲措、拉布普、热当与下游的南木林。2021 年高值集中在水库附近的南部(甲措、热当一带),2022 年全域偏低,2023 年整片转为高值(橙色系),2024 与 2025 年又回到偏枯的浅绿色调。

丰枯年环流合成
丰枯年环流合成

丰水年(左列)与枯水年(右列)的合成场。(a)(b) 为 500 hPa 位势高度距平叠加风场:丰水年流域处于弱负距平区、上游有槽活动,枯水年则被大片正距平(红色,40—60 gpm)覆盖。(c)(d) 为水汽通量散度距平叠加水汽通量距平:丰水年流域处在蓝绿色的辐合区,枯水年转为红色辐散区,对应正文的 −1.9×10⁻⁵ 与 1.5×10⁻⁵ kg·m⁻²·s⁻¹。

作者承认的限制

一句话判断

对一个此前基本没有乡镇观测的高原水库流域,这篇把「代表站能不能用」「降水有多集中」「极端值到什么量级」「丰枯年环流差在哪」四件事一次性补齐了,10—30 年一遇的 45.3—68.5 mm 是可以直接进防洪规划的数字;但要拿去做季节差异或驱动因子贡献的定量判断,作者自己说了现在还做不到。

在观潮中打开(含配图与交互)→