原文作者：朱友娟

　　摘要：降水量（p）是估算水资源各分量的一个重要物理量。选取了塔里木河流域源流区5个站点近50年的降水量序列资料，进行了kendall秩次相关检验和小波分析，揭示了塔里木河流域源流区降水变化的多时间尺度的复杂结构，从而为未来流域源流区降水量及相关水系径流量变化的预测提供依据。结果表明：近40年来流域源流区降水量呈现增加趋势，各站降水量的变化趋势显著性程度有所不同。在低频变化层次上，沙里桂兰克、卡拉贝利、乌鲁瓦提3站的大尺度周期震荡比较明显。在高频变化层次上，各站都存在5～10年不等的震荡主周期，但局部化特征有所差异。
 　　关键词：塔里木河流域；降水量；kendall秩次检验；小波分析
　　降水的形成过程是一个涉及到水文、气象、地形、植被等因素的复杂过程，既受到确定性因素的作用又受到随机因素的影响，因而其变化非常复杂，对其未来的精确描述也十分困难。传统的统计判别方法很难判别其演化过程，而被誉为“数学显微镜”的小波分析由于其在时域和频域同时具有良好的局部化性质，能同时在时间和频率2个尺度上展现数据序列的变化特征，特别适合做非线性现象的检测与分析，尤其是在多尺度的气候诊断中的应用，已取得一些很好的成果。如尤卫红[1]用morlet小波和marr小波分析了云南近百年气温的变化；邓子旺等[2]用db16正交小波对南方涛动与北大西洋涛动进行了分析；戴新刚等利用morlet小波对华北降水频谱变化及其可能机制进行了分析[3]。笔者采用morlet小波对塔里木河流域源流区降水量的演变规律进行了初步的分析，探讨了源流区近50年降水时间序列多时间尺度的变化特点，这对于该流域相关水系未来径流量的预测以及合理利用流域水资源具有重要参考价值。WwW.11665.com
 　　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，年际变化较剧烈（cv=0.38），见表1。从5年滑动平均曲线上可清楚看到1965-1970年降水有一次明显的增加（见图2b）。20世纪50年代中期-80年代中期，降水距平以负距平值为主，20世纪90年代初-2000年降水量持续增加，见图3b。从其小波变换图的高频部分上可看出，年际时间尺度在时间域上分布不均匀，降水周期由较短向较长周期发展。
 　　沙里桂兰克站的5年降水倾向率为7.7645 mm，增加趋势最为明显，见表1。其降水年际变化相对稳定。从距平值图2c上可看出从20世纪50-80年代初期，降水距平以负距平值为主，即偏干旱年出现的频率大于偏湿润年。其中20世纪70-80年代降水偏少，80年代末期降水逐渐增加，一直到90年代末，以正距平值为主。从沙里桂兰克1957-2000年降水量变化趋势图上可见（图3c），总体上降水量呈上升趋势。由小波变换图（图4c）可见其降水量在20世纪60-90年代初时，5～8年时间尺度上的周期震荡最明显。
 　　乌鲁瓦提站的5年降水倾向率为0.197 mm，增加趋势最不明显，年际变化最为稳定（cv=0.35），见表1。20世纪50年代中期-60年代末，降水距平以正距平值为主（图3d），20世纪70年代初-80年代中后期，降水距平以负距平值为主。20世纪80年代以前，主要以5年左右的周期为主，从大的周期尺度来看，整个时间序列来看，期间经历了由多—少—多3个循环交替（图4d）。

　　协合拉站的5年降水倾向率为3.5610 mm，其滑动平均曲线波动不是很明显（见图2e），20世纪70年代以前以负距平值为主（图3e）。其小波变换图上以5年周期为主（图4e），20世纪70年代中期以前还存在一个15年左右的长周期。
 　　3 初步结论
　　3.1 源流区乌鲁瓦提、卡拉贝利、沙里桂兰克站降水年际变率较为接近。从各站多年降水量5年滑动平均曲线的变化过程大致可以看出：塔里木河流域

源流区的站点的降水量都存在逐步增多的趋势；各站都曾发生过几段多雨期和少雨期，且各段的持续时间并不固定。根据kendall秩次相关检验结果来看，除沙里桂兰克外，其余4站的降水量增加趋势在统计上都不明显。周期分析表明，各站的降水具有不同的长短周期震荡。小波分析很清楚地反映了各站降水量不同时间尺度震荡的交替作用。在高频变化层次上，各站都存在5～10年左右的震荡周期，其中卡群、卡拉贝利的短周期震荡比较明显。沙里桂兰克、乌鲁瓦提两站降水量时间序列的高频部分在时间域上分布不均匀，具有明显的局部化特征。在低频变化层次上，沙里桂兰克、卡拉贝利、乌鲁瓦提三站的大尺度周期震荡比较明显。虽然各站的降水都有增加趋势，但在2000年附近，卡拉贝利、沙里桂兰克、协合拉站的降水小波系数在低频部分为虚线，这预示着21世纪的最初几年，这几个站点可能处于一个降水偏少的时期。
 　　3.2 塔里木河流域自20世纪六七十年代以来，灌溉面积不断扩大，引水量不断增大，平原区的水蒸发量增大，空气湿度增大，当出现有利的降水环流时，首先会形成更多的降水量，其次在降水环流条件弱，原来不能产生降水的情况下，可能出现弱的降水天气。因此，从人类活动的角度来看，这有可能是引起源流区降水量增加的一个原因。
 　　3.3 20世纪80年代中期以来，西北地区的气候转型问题已引起众多学者的关注[13-16]，一些学者认为新疆的气候开始由暖干向暖湿转型，这种转变首先出现在受西风环流带影响的地区，并且以天山西部最为强劲。还有一些学者认为流域径流量的增加只能说明塔里木河正处于一个丰水期。笔者对降水量增加趋势的显著性检验还不能证明源流区降水增加的假设。因此，“转型学说”是否适用于塔里木河流域源流区还有待于进一步深入探究。
 　　3.4 一般而言，不同的周期变化是由不同的机制产生的，降水量多时间层次结构的变化也可以说明：不同的机制在不同的时期内作用强度不同，就形成了降水量变化的非周期性和复杂性。因此，塔里木河流域降水量特别是源流区降水量的演变规律及其产生机制的研究还有待于进一步深入探究。
 　　参考文献
　　[1] 尤卫红.气候变化的多尺度诊断分析和预测的多种技术方法研究[m].北京：气象出版社，1998：9-42.
　　[2] 邓自旺，施能，王永波.南方涛动（so）与北大西洋涛动（nao）百年变化的正交小波分析[j].热带气象学报，2001，17（1）：57-64.
　　[3] 戴新刚，汪萍，张培群，等.华北降水频谱变化及其可能机制分析[j].自然科学进展，2003，13（11）：1182-1188.
　　[4] 徐海量，叶茂，宋郁东，等.塔里木河流域水资源变化的特点与趋势[j].地理学报，2005，60（3）：487-494.
　　[5] 陈亚宁，崔旺诚，李卫红，等.塔里木河的水资源利用与生态保护[j].地理学报，2003，58（2）：215-222.
　　[6] 飞思科技产品研发中心.小波分析理论与matlab7实现[m].北京：电子工业出版社，2005，29-37.
　　[7] 丁晶，邓育仁.随机水文学[m].成都：成都科技大学出版社，1988.
　　[8] jury m r， melice j l. analysis of durban rainfall and nile river flow 1871-1999[j].theoretical and applied climatology，2000，67：161-169.
 　　[9] torrence c. a practical guide to wavelet analysis[j]. bulletin of the american meterological society，1997，79（1）：61-78.
 　　[10] 宋郁东，樊自立，雷志栋，等.中国塔里木河流域水资源与生态问题研究[m].乌鲁木齐：新疆人民出版社，2000：39.
　　[11] kendall m， gibbons j d. rank correlation methods[m]. london： a charles griffin book，1990：10-35.
 　　[12] 张家宝，史玉光.新疆气候变化及短期气候预测研究[m].北京：气象出版社，2002.
　　[13] 施雅风，沈永平，胡汝骥.西北气候由暖干向暖湿转型的信号、影响和前景初步讨论[j].冰川动土，2002，24（3）：219-226.
　　[14] 胡汝骥，姜逢清，王亚俊，等.新疆气候由暖干向暖湿转变的信号及影响[j].干旱区地理，2002，25（3）：194-201.
　　[15] 胡汝骥，王东亮，冯国华，等.怎能让塔里木河下游干涸的悲剧重现中游地区：2002年塔里木河流域3条源流区间耗水评析[j].干旱区研究，2004，21（3）：199-203.
 　　[16] 胡汝骥，马虹，樊自立，等.新疆水资源对气候变化的响应[j].自然资源学报，2002，17（1）：22-26.