塔里木河流域源流区降水量变化的小波分析.docVIP

塔里木河流域源流区降水量变化的小波分析 　　摘要：降水量（P）是估算水资源各分量的一个重要物理量。选取了塔里木河流域源流区5个站点近50年的降水量序列资料，进行了Kendall秩次相关检验和小波分析，揭示了塔里木河流域源流区降水变化的多时间尺度的复杂结构，从而为未来流域源流区降水量及相关水系径流量变化的预测提供依据。结果表明：近40年来流域源流区降水量呈现增加趋势，各站降水量的变化趋势显著性程度有所不同。在低频变化层次上，沙里桂兰克、卡拉贝利、乌鲁瓦提3站的大尺度周期震荡比较明显。在高频变化层次上，各站都存在5～10年不等的震荡主周期，但局部化特征有所差异。 　　关键词：塔里木河流域；降水量；Kendall秩次检验；小波分析 　　降水的形成过程是一个涉及到水文、气象、地形、植被等因素的复杂过程，既受到确定性因素的作用又受到随机因素的影响，因而其变化非常复杂，对其未来的精确描述也十分困难。传统的统计判别方法很难判别其演化过程，而被誉为“数学显微镜”的小波分析由于其在时域和频域同时具有良好的局部化性质，能同时在时间和频率2个尺度上展现数据序列的变化特征，特别适合做非线性现象的检测与分析，尤其是在多尺度的气候诊断中的应用，已取得一些很好的成果。如尤卫红[1]用Morlet小波和Marr小波分析了云南近百年气温的变化；邓子旺等[2]用DB16正交小波对南方涛动与北大西洋涛动进行了分析；戴新刚等利用Morlet小波对华北降水频谱变化及其可能机制进行了分析[3]。笔者采用Morlet小波对塔里木河流域源流区降水量的演变规律进行了初步的分析，探讨了源流区近50年降水时间序列多时间尺度的变化特点，这对于该流域相关水系未来径流量的预测以及合理利用流域水资源具有重要参考价值。 　　1 研究资料与方法 　　1.1 研究区域概况 　　塔里木河流域是中国最大的内陆河流域，主要包括塔里木河干流和源流区域。源流区河流主要由叶尔羌河、阿克苏河、和田河等汇合而成，其中阿克苏河由昆马立克河与托什干河两源汇合而成，见图1。塔里木河流域总面积102万km2（国内流域面积99.6万km2）。流域内气候干燥、干旱少雨、昼夜温差大，是典型的大陆性气候，其生态系统是一个由高山冰川―高山冷湿草甸―中山湿润森林―低山半干旱灌草―平原干旱荒漠绿洲构成的脆弱生态系统，也正是由于高山冰川和山区降水形成了塔里木河水系，而在径流组成中，冰川融水占48.2%，雨、雪混合补给占27.4%。因此，山区气温和降水是影响塔里木河流域源流水量大小最重要的环境要素[4-5]。 　　1.2 研究资料 　　由于新疆特殊的地理环境，塔里木河流域水文气象站点分布稀少，笔者采用塔里木河流域源流区5个与水系发生联系的气象站近40多年来的年降水量资料来代表源流区降水量变化趋势。5个站分别是卡群、协合拉、沙里桂兰克、乌鲁瓦提、卡拉贝利，对应水系见图1所示。 　　1.3 研究方法 　　1.3.1 小波分析 　　小波变换是把小波母函数（t）作位移b、在伸缩尺度a下与所要分析的函数f（t）（笔者以年降水量）作内积，对于函数f（t） ∈L2（R），连续小波变换的公式[6-8]： 　　2 结果与分析 　　图2a～e给出了5个站降水量时间序列的线形函数趋势、5年滑动平均特征，图3a～e为各站多年降水距平值图。图4a～e显示了塔里木河源流区各站降水量在不同时间尺度上的周期震荡。信号的强弱通过小波系数的大小来表示。从小波变换图上可看出不同时间层次上降水量的变化情况，图中实线表示降水增加，虚线表示降水的减少；各图上半部分为低频，等值线相对稀疏，对应较长尺度周期的震荡，各图下半部分是高频，等值线相对密集，对应较短尺度周期的震荡。实际降水量的变化可以看成各个层次变化叠加的结果，它的大小取决于各时间层次上的变化情况及其对时间降水量的方差贡献。 　　卡群站的5年降水量趋势线倾向率为1.9070 mm，多年的年平均降水量为69.95 mm，年际变化比较剧烈（CV=0.49）。见表1。经5年滑动平均后卡群站降水量的年代际变化得以明显表现（见图2a）。20世纪60年代中期-70年代中期有2次明显的增加。20世纪70年代末-80年代末有2次明显的减少，期间略有回升。20世纪90年代回升趋势比较明显。降水距平值（见图3a）：20世纪70年代初-80年代中期以负距平值为主，最大降水正距平值为71.25 mm，出现在1974年，最大降水负距平值为-61.05 mm，出现在1985年。由图4a可见，其年际尺度和年代际尺度特征十分明显。60年代中期以后主要存在5～8年周期，且周期由长变短。在2000年附近，系数线为虚线，表明降水量减少，这与距平值图反映的情况一致。 　　卡拉贝利站的5年降水倾向率为5.4155 mm，年际变