1. 引言

黄河源区是指位于青藏高原东北部(95˚50'E~103˚30'E，32˚30'N~30˚00'N之间)的黄河干流唐乃亥以上流域 [1] 。自20世纪末，由于气候变化和人类活动的影响，源区的生态系统和水源补给系统了发生明显变化，如冰川消融、雪线上升等，加快了地表水量减小、冻土层地下水位的下降、湿地退化等现象的发生，黄河源生态系统对气候变化非常敏感 [2] [3] 。径流的影响因子主要是气候变化因子和人类活动。气候变化因素主要包括气温、降水、相对湿度、日照、风速等；人类活动因素主要包括人类对流域下垫面的改变和水资源的开发利用等 [4] 。

国际上从20世纪80年代中期展开了关于气候变化对水文、水循环和水资源的影响的研究 [5] 。1977年，美国国家研究协会(USNA)组织讨论了气候变化与供水之间的相互影响与联系。气候变化对水文影响的不确定性评估引起人们的高度重视。Shen [6] 通过研究全球气候模式(GCM)选择的不确定性的时间变化(主要是长期的)，以预测气候变化对水文的影响，结果表明，20个GCM预测的未来气候变化(降水和温度)和水文响应具有高度的不确定性，并且不确定性随时间显着增加。Nepal [7] 使用12000模型，在DudhKoshi子流域研究了目前并预测未来条件下不同径流成分的水文特征和贡献，模型预测结果表明，至本世纪中叶，流量将增加13%，随后略有下降；本世纪末蒸发量将增加16%。1985年的Vinach会议，我国加进了气候变化对水资源影响研究 [8] 。刘昌明等 [9] [10] 就黄河源区降水和径流的关系以及径流对土地覆盖情况和气候变化的响应展开研究，根据SWAT模型模拟结果，随着植被覆盖度的增加，流域年径流量增加，蒸散量减小。

黄河源区处在特殊地理位置以及脆弱的生态环境中，形成其独特的水文系统和天然的冻土环境 [11] 。本文分析了源区的气温的变化趋势及突变年代，降水的季节分配和变化趋势，以及径流的年内变化和年际变化，分析气候变化对源区的径流影响。对黄河源区径流变化进行系统分析，探讨其变化趋势和原因，揭示气候变化对水生态环境的影响，可以为流域水资源管理提供科学依据。因此，有必要对黄河上游流域径流变化进行系统的分析，探讨其变化趋势和原因，揭示气候变化对水生态环境的影响，为流域水资源管理提供科学依据。所以，开展气候变化对黄河源的水文影响分析对水资源开发利用以及河道生态保护非常重要。

2. 研究区概况与研究方法

2.1. 研究区概况

黄河发源于青藏高原巴颜克拉山北麓的约古宗列曲，本研究中，黄河源区选用唐乃亥水文站以上的部分，流域面积为12.2万km²。黄河源区地势呈现西南较高东北较低的特点，整体海拔较高，高程大致在2656~6350 m范围内，属于高原大陆性气候，太阳辐射比较强、日照时间较长、大气压力较低、昼夜温差大 [12] 。源区属于半湿润地区，但区域内降水分布不均，呈东南多西北少的特点。独特的地理位置以及气候特征，使得该区冰川广布、冻土较为发育，以及星罗棋布的湖泊湿地，发挥了黄河源区涵养水源和补给水源的重要功能 [13] 。黄河源区的冻土地貌处于季节冻土到片状连续多年冻土的过渡区，相比高原腹地，多年冻土具有温度更高、厚度更薄的特点 [14] 。黄河沿、玛多县城等玛多盆地1970年代曾经属于多年冻土，1990年复位钻探研究显示多年冻土大部分已消融，冻土分布边界线也分别向西和向北移动 [15] 。

2.2. 研究方法

径流过程是水循环最重要的环节之一，其变化过程受降雨、蒸发、下渗、下垫面条件以及人类活动等诸多因素的影响。

2.2.1. 特征参数分析

反映经流年内分配不均的特征值较多，主要用径流年内分配不均匀系数C_v。不均匀系数C_v的计算公式为：

$C_v=\frac{\sqrt{\frac{1}{12}{\displaystyle {\sum }_{i=1}^{12}{\left(Q_i-\bar{Q}\right)}^2}}}{\bar{Q}}$

式中：Q_i为月平均流量，m³/s； $\bar{Q}$为年平均流量，m³/s；C_v是表示径流年内分配不均的一个指标。C_v越大，表示径流年内分配越不均匀；反之，则径流年内分配越均匀。

年内分配完全调节系数C_r是另外一种年内分配的指标。其计算公式如下：

$C_r={\displaystyle \underset{i=1}{\overset{12}{\sum }}\Psi \left(t\right)}\left[Q\left(t\right)-\bar{Q}\right]/{\displaystyle \underset{i=1}{\overset{12}{\sum }}Q(t)}$

式中，

$\Psi \left(t\right)=\{\begin{array}{l}0,Q\left(t\right)<\bar{Q}\hfill \\ 1,Q\left(t\right)\ge \bar{Q}\hfill \end{array}$

年内分配完全调节系数C_r越大，则年内水文气象因子分配越不均匀。

2.2.2. 趋势分析

1) 多年滑动平均

多年滑动平均值序列可以近似表示全序列的变化趋势，此方法是由时间序列的平滑值来表示的，对一个样本数量为n的一个长序列x，多年滑动平均值序列可以表示为：

$x_j=\frac{1}{k}{\displaystyle \underset{i=1}{\overset{k}{\sum }}{x}_{i+j-1}}\left(j=1,2,\cdots ,n-k+1\right)$

式中，k值可以取非零正整数。经过多年滑动平均的计算，通过分析滑动平均序列的曲线图来判断整个时间序列的趋势。

2) 累计距平法

累计距平法又称差积曲线法，对于一个时间序列x，某一时刻t的累积距平值表达如下：

$x_t=\underset{i=1}{\overset{t}{{\displaystyle \sum }}}\left(x_i-\bar{x}\right)\left(t=1,2,\cdots ,n\right)$

式中，x均值表示为

$\bar{x}=\frac{1}{n}\underset{i=1}{\overset{n}{{\displaystyle \sum }}}{x}_i$

在绘制累积距平曲线时，需要求出n个时刻的累积距平值，然后进行趋势性分析；如果曲线呈现上升趋势则表示距平值在增加，相反如果曲线呈现下降趋势则表示距平值减小。通过曲线的变化可以判断出序列大致突变的时间以及长期持续变化的情况和演变趋势 [16] 。

3. 研究结果

3.1. 气候因子变化分析

3.1.1. 气温变化分析

1) 气温年际变化

黄河源平均气温时间序列变化图如图1所示，由图可以看出黄河源流域1960~2014年气温整体呈递增趋势。年平均气温变化倾向率为0.36℃/10a，查相关系数检验表得α = 0.05，n = 55时r_a= 0.266，而相关系数r = 0.79 > r_a，表明研究的系列表现为显著的递增趋势。

2) 气温累积距平变化过程

气温累积距平图如图2所示，可以看到在本次研究中，气温累积距平均为负值，累积距平值都呈现出先降低后增加的趋势。从平均气温来看，源区气温累积距平值出现了2个转折点，在1986年发生转折而后平稳波动，在90年代末，发生第二次转折，此后开始呈现上升趋势。年平均最低气温在1986年发生转折，年平均最高气温在90年代末发生转折，而后呈现显著上升趋势。

3.1.2. 降水变化分析

1) 降水年内变化

由于气候的季节性波动，降水呈现季节性变化，使得降水年内分配不均，降水的年内分配特征通常采用各

月、季占全年降水量的百分比等来反映。图3显示黄河源流域降水年内分配差异也较大，多年最大平均月降水量发生在7月，月降水量为99.2 mm，占多年平均年降水量的21.7%。根据黄河源区降水径流特性，研究划分汛期为5月至10月，占年降水量的90.9%；而11月份至次年2月份的降水量最少，仅占全年降水量的3.1%。黄河源流域降水季节分配也极不平衡，夏季降水占全年降水量的58.1%，冬季降水只占全年降水量的只有2.2%。

2) 降水年际变化

黄河源流域平均降水量年时间序列变化图如图4所示，由图可以看出1960~2014年降水量呈现波动递增趋势。在1990年之前序列点系处于趋势线两侧，呈现震荡增加趋势，在1990年之后降水量处于趋势线以下，说明在1990~2006年之间降水量呈现了下降趋势。年平均降水量变化倾向为6.80 mm/10 a，当α = 0.05，n = 55时，r_a= 0.266，而相关系数r = 0.19 < r_a。表明原系列整体存在不显著的递增趋势。

由图5可知黄河源汛期降水量变化基本平稳，在80年代末出现下降趋势，21世纪初表现出上升趋势。结论表明黄河源在80年代末到21世纪初，整体比较干旱，在21世纪初出现恢复变湿的趋势。黄河源非汛期降水量整体呈现增加趋势，当α = 0.05，n = 55时，r_a= 0.266，而相关系数r = 0.52 > r_a。表明原系列整体存在显著的递增趋势。

3.2. 径流变化分析

3.2.1. 径流年内变化

径流的年内分配又称径流的年内变化或季节分配。天然河流由于受气候因素及流域调蓄能力有关的下垫面因素的影响，径流量在年内的分配时不均匀的。

1) 径流百分比法

根据实测径流资料(日平均资料)统计出相应年度的年径流量，然后确定计算时段(月或季)，并根据水文站径流资料统计出相应时段内的径流量，可算出不同时段径流量的百分比。

由表1可知，径流年内分布不均，其多年平均年径流量是199.89亿m³，而且绝大部分集中在6~10月份，其径流量能占到年径流总量的70%。由图6可知，6月份平均径流量为23.272亿m³，占全年径流量的11.6%，7月份径流量最大，多年平均径流量为34.559亿m³，占全年径流量的17.3%，8月份平均径流量为28.473亿m³，占全年径流量的14.2%，9月份平均径流量为31.177亿m³，占全年径流量的15.6%，10月份平均径流量为25.391亿m³，占全年径流量的12.7%。研究时期内，流域12月、1月、2月径流量均很少，三个月的径流量仅占年流量的7.4%，其中最小月径流量出现在2月份，仅占年径流量的2.19%。同时由最大正负误差线，可以看出在研究时期内，历史最大月径流发生在9月份，最低月径流发生在12月份。

2) 年内分配特征

由表2可以看出，黄河源流域径流年内分配不均性和集中度普遍比较高，说明径流年内分配相对比较集中，计算集中期显示出现在8月。黄河流域1960~2014年平均月径流量的不均匀系数和完全调节系数分别是0.71和0.3，并且与集中度具有较好的同步性变化规律。不均匀系数值越高，集中度值也越高，径流量具有明显的丰枯季节变化。

3.2.2. 径流年际变化

1) 径流变化过程

通过分析图7可知，唐乃亥站1960年~2014年径流呈现递减趋势，下降斜率为−5.684/10a，相关性系数为0.175。假设序列置信度α = 0.95，n = 55时r_a= 0.266，相关性系数r = 0.175 < r_a，表明原序列存在不显著的递减趋势。这种下降趋势在80年代中期到20世纪初滑动平均值下降的比较明显，下降趋势比较显著。

2) 累积距平法

通过黄河源区年径流量、汛期、非汛期径流量累积距平随时间分布图(图8)可知，黄河源区年径流量与汛期

径流量累积距平曲线变化趋势十分相近。结合径流滑动平均的分析结果，同样得到汛期多年变化趋势与多年平均径流量变化趋势十分接近，即可间接说明黄河源区降水对径流的影响较大。由累积距平图可知，径流系列转折点发生在1990年左右。年径流量和汛期径流量主要有2个阶段，总体表现为1960~1998年为波动性显著增加，说明该阶段丰水年份多于枯水年份；1998~2014年间为显著减少趋势，说明该阶段枯水年份要多于丰水年份。非汛期累积距平分为两个阶段，1984年之前呈现缓慢增加趋势，1984~1990年基本保持在稳定的范围之内，在1990年以后非汛期累积距平出现逐渐下降趋势，但变化幅度明显小于年径流量和汛期径流量。

4. 结论

本文对黄河源区的气候因子变化和径流变化进行了分析，主要结论如下：

1) 在气候因子变化分析中，黄河源区年平均气温在研究系列内呈现显著的增加趋势。由降水年内变化分析可知，黄河源流域降水年内分配差异较大，降水的年内分布以5~10月份最多，而11月份至次年2月份的降水量最少。降水年际变化分析可知，黄河源区降水量整体呈现不显著的增加趋势。

2) 由年内径流分布可知黄河源区多年平均最大月径流发生在7月份。多年平均不均匀系数比较大，集中度比较高，说明黄河源区径流分配有明显的丰枯季节变化。由径流年际变化可知，研究系列中黄河源区径流量呈整体减少趋势，这种趋势在90年代至21世纪中期尤为明显。

基金项目

国家自然科学基金青年基金项目，项目名称为盐渍化农田水–热–盐传输过程与作物生长的特征及模拟(编号413615048)。

参考文献