塔里木河流域降水时空变化分析.docxVIP

塔里木河流域降水时空变化分析 塔里木河流域总面积为107.54km2，包括五个生态系统：山地、绿洲、水域、森林（灌木和草）和沙漠。这是一个封闭的内部和水循环和平衡的区域。流域内人口占新疆人口44.00%;耕地占新疆耕地面积39.03%;国内生产总值占全新疆国内生产总值的26.58%[l], 是新疆主要经济区。旱涝是塔里木河流域的重大灾害问题。近年来, 许多研究得出新疆大部地区降水增多, 塔里木河流域气候变湿, 这种变化在该流域内是否达到突变水平, 是否引起旱涝等级、趋势变化和改变, 本文采用Z指数变换法和塔里木河流域1961～2001年资料, 分段计算41年来4～9月降水的旱涝等级, 分析旱涝时空变化及特点, 为这一地区较长时间旱涝预测和水资源调配奠定基础。 1 单站旱涝指标确定 文中灾害资料取自文献和各年新疆年鉴, 降水资料选用新疆气象局整编的塔里木河流域13个气象站降水数据。考虑各站资料完整性和可比较性, 统一选用无缺测的196l～2001年4～9月月降水记录。对塔里木河流域降水均一性检测, 仅若羌站1966年出现不均一断点, 但可能是虚假的, 查看若羌台站历史沿革, 1961～2001年观测降水方法和仪器并未改变, 而判断若羌降水不均一的参考站与之距离 (最大相距900 km) 、海拔 (最大相差1 097.3 m) 和降水相差很大, 可能是误判。用方差比X2检验若羌1966年前后资料变率差异显著性表明, 在信度0.05水平下, 1966年前后资料没显著差异, 能认为若羌降水序列是均一的。为表征旱涝变化的区域特征, 以阿合奇、乌什、阿克苏站代表塔里木河的主要源流区;以阿拉尔、拜城、沙雅、库车、新和站代表上游区, 以轮台、库尔勒、尉犁站代表中游区;以铁干里克、若羌气象站代表下游区进行分段统计。 本文用文献Z指数作为单站的旱涝指标划分早涝级别。假设某时段降水量服从Person-Ⅲ型分布, 对降水量X进行正态化处理, 可以把概率密度函数Person-Ⅲ型分布转换为以Z为变量的标准正态分布, 其转换公式为: zi=6Cs(Cs2xi+1)1/3?6Cs+Cs6zi=6Cs(Cs2xi+1)1/3-6Cs+Cs6 其中,Xi为降水标准化变量,Cs为降水偏态系数, 二者都可通过降水资料序列求得。根据z变量正态分布曲线, 划分7个等级其相应z指数界限值作为单站旱涝指标。据文献计算得塔里木河源流区、上游、中游和下游降水空间变异系数较小, 分别是0.62、0.43、0.36和0.48, 认为各段区域内降水分布均匀, 区域旱涝指数由下式确定:z=1m∑i=1mZiz=1m∑i=1mΖi 式中,Zi为区域内各站Z指标,m为区域内总站数 (表1) 。 文中降水量旱涝等级突变检验使用Mann-Kendall方法。取95%信度水平, 当序列c1和c2超过Y=±1.96信度线, 且c1和c2曲线相交于信度线之间, 则交点是突变点开始;对旱涝周期变化分析使用最大熵谱分析, 具有频谱短且光滑、分辨率高等优势。对旱涝趋势分析使用气候趋势系数: r=∑i=1n(xi?xˉ)(i?tˉ)/∑i=1n(xi?xˉ)2∑i=1n(i?tˉ)2???????????????????√r=∑i=1n(xi-xˉ)(i-tˉ)/∑i=1n(xi-xˉ)2∑i=1n(i-tˉ)2 其中,xˉxˉ是平均值,tˉ=(n+1)/2tˉ=(n+1)/2。 2 结果和分析 2.1 降水突变的趋势 近40年来, 新疆平均降水量表现出上升趋势, 20世纪80年代中期以来增多尤其明显, 降水量增加成为新疆当代气候变化的重要特点。 图1a可以看出, 1990年出现降水增多突变点, 1990年以前降水增多趋势和减少趋势的时期基本交替出现, 且大致相等;1990～1999年降水呈持续增多趋势, 1998年以来, 增多达到显著水平;上游 (图1b) 降水增多突变点出现在1975年, 1975年以前, 除1968～1970期间降水减少外, 其余时段属降水增多期, 但增多不显著, 1975年以来降水增多波动性增长, 1982年以来达到显著水平;中游 (图1c) 降水增多突变点出现在1964和1973年, 1964年突变后增长的趋势时间不长, 1966年增多达到显著水平后, 随即在1967年回落, 低于显著水平。1973年突变之后, 降水增多波动性增大, 1991～1998年达到显著水平, 之后低于显著水平。下游 (图1d) 降水变化与源流区、上游和中游的相同, 降水基本呈增多的趋势, 突变时间是1969年, 1974年以后达到显著水平, 之后波动性下降, 其中1981、1993年又达到增多显著水平。 对流域内各站4～9月降水量突变检验表明:源流区三个站在20世纪80年代中期降水量减少趋势由增大转变