Page 98 - 2022年第53卷第10期

图 ９的箱型图分别展示了基于 Ｂｕｄｙｋｏ － Ｆｕ框架分解法和水能分配法的不同下垫面流域实际蒸散发
变异归因结果（贡献率）。经 ｔ检验，两类下垫面类型区的归因结果存在较显著差异。对于大部分低 ｎ
值流域，实际蒸散发变化主要归因于下垫面因素，基于水能分配法的计算结果表明，７５％低 ｎ值流域
下垫面对流域实际蒸散发变化的贡献率 Ｃ超过 ５６．８％（基于分解法的为 ６２．９％），５０％低 ｎ值流域的 Ｃ Ｌ
Ｌ
大于 ７５．９％（基于分解法的为 ７６．８％）；在高 ｎ值区内，下垫面因子对流域实际蒸散发的变化贡献略大
于气候因子的贡献，７５％的高 ｎ值流域下垫面对流域实际蒸散发贡献超过 ３２．１％（基于分解法的为 ２９．
６％），有一半高 ｎ值流域的 Ｃ大于 ５５．３％（分解法为 ５３．４％）。
Ｌ

图 ９ 各下垫面类型区内气候和下垫面因子对实际蒸散发变化的贡献率（水能分配法）
注：在 Ｐ＜０．０５的水平上有显著相关性；在 Ｐ＜０．０１的水平有显著相关性。

３．４ 气候及下垫面变化对蒸散发变异的影响方向分析 图 １０描述了长时间序列内各流域水量及能量
分配模式变化，由 １９００—１９５０年和 １９５１—２００８年（ｆ，ｑ）均值的箭头指向及长度揭示了两个时段间各
流域水能状态的变化方向及幅度，以此表示流域水能状态在 ｆ － ｑ空间坐标系中的动态演变过程。

（点与箭头分别表示各流域在 １９００—１９５０和 １９５１—２００８年期间的状态）
图 １０ 各流域水量及能量分配比的多年平均值

对此 １９００—１９５０年和 １９５１—２００８年两时段，大部分湿润流域的水能状态轨迹指向右上方，表明
流域的气候条件变得更湿润，下垫面特征参数有所增大，且实际蒸散发变化主要受下垫面影响。其
中，有一部分下垫面特征参数较低（ｎ＜２）的流域，如乌达河流域的 ｎ值呈明显的降低态势，且在下垫
面的驱动作用下，流域水热状态朝着水量及能量利用率更低（即 ｑ和 ｆ值更低）的方向移动。对于下垫
面特征参数较高的湿润流域而言，除呈现 ｎ值增大的态势（箭头指向右上角）外，一些流域内气候因素

３
— １ ２ ６ —