滇池流域降水与径流的因果关系分析.docxVIP

滇池流域降水与径流的因果关系分析 全球气候变化及其对生态环境、植物、农业和水资源的影响已成为国内外的研究热点。姚凤梅等利用DSSAT中的CERES-rice作物模式分析发现, 在A2和B2情景下, 昆明和海口 (海南省) 的水稻产量下降幅度大于南方其他站点。宁夏引黄灌区1961年以来的调查统计表明, 气候变暖对春小麦单产的贡献率为-2.6%, 不利于春小麦生产。吴绍洪等利用参考作物腾发量和干湿指数分析了1971-2000年中国陆地表层的干湿状况。在气候变化科学中, Mann-Kendall法主要用于检验水文、气象等要素的长期变化趋势。Granger因果关系检验法过去多用于计量经济学中, 现已逐渐在地理、水资源方面得到应用。 西南地区由于云贵高原、青藏高原及中南半岛的低纬度不同地形地貌共同作用, 加上靠近太平洋和印度洋两个热带气团发源地, 气候变化异于中国北方等其他地区, 如池再香等发现, 1961年以来黔东南地区降水增加、气温下降。马振锋等认为西南地区气候变化与全球变暖存在非同步性。孙永亮等的研究发现青海湖流域的气温升高、降水增加, 但主要入湖河流径流量却减少。滇池流域是滇中经济核心区, 流域水资源开发利用远远超过水环境承载能力, 滇池水质长期处于Ⅴ类～劣Ⅴ类, 研究滇池流域的气候变化及其对水资源利用的影响, 对于全社会进一步明晰滇池高原湖泊生态修复治理的方向和重点具有极其重要意义。 1 材料和方法 1.1 农业统计年2000年 (1) 滇池流域内昆明站1953-2000年的逐日平均最高温度、平均最低温度、平均温度、相对湿度、日照时数、风速、降水量、蒸发量等资料。 (2) 2000年度滇池流域内各县 (区) 的农业统计年报。 (3) 滇池的1961-2006年逐月平均水位、出湖水量。 (4) 滇池出口及滇池下游螳螂川的海口、蔡家村两个水文站1953-2005年的实测逐月径流量系列。 1.2 h方程计算方程 参考作物腾发量ET0采用Penman Montieth方程计算。水稻灌溉定额为生育期内各次灌水定额之和。农业综合灌溉定额指灌区内同一时期各种作物灌水定额按种植面积的加权平均值, 并考虑大、小春作物轮作和复种的影响。 1.3 文气象资料趋势检验图1 Mann-Kendall非参数秩次相关检验法 (简称M-K检验) 主要用于水文气象资料趋势检验, 如水质、流量、气温和降雨等。给定显著性水平α, 在正态分布表查临界值Mα/2, 若|M|Mα/2, 趋势不显著;若|M|Mα/2, 即趋势显著。且当M为正时表示上升或增加的趋势, 反之则为下降或减少趋势。 1.4 时间序列平稳性检验adf 时间序列之间的Granger因果关系定义为:设有两个时间序列变量X、Y, 若采用X的过去值情况下, 得到Y的预测值比在不采用X的过去值情况下所得到的Y预测值更好, 则说X是Y的原因, 或X会引起Y。反之, 则X对于Y在统计上独立, 在时间上不相关。设有两个平稳的时间序列变量X、Y, 检验它们间有无Granger因果关系的模型如下 Yt=λ+∑i=1mαXt?i+∑j=1nBjYt?j+μt(1)Yt=λ+∑i=1mαXt-i+∑j=1nBjYt-j+μt(1) 式中:λ为截距;m、n为滞后值;t为时间;μt为误差项。变量X、Y为含一个单位根的不平稳的I (1) 时间序列, 检验它们间有无Granger因果关系可用如下的第一差分滞后向量自回归 (VAR) 模型: ΔYt=λ+∑i=1mαΔXt?i+∑j=1nBjΔYt?j+μt(2)ΔYt=λ+∑i=1mαΔXt-i+∑j=1nBjΔYt-j+μt(2) 式中, 所有项都是平稳的, 故消除了时间序列X、Y的不平稳问题。究竟采用式 (1) 或者式 (2) 分析时间序列间是否存在Granger因果关系, 取决于时间序列是否含有单位根, 即时间序列的平稳性。时间序列平稳性检验通常采用ADF (Augmented Dickey-Fuller) 检验技术, 回归方程为: ΔYt=α+βT+?Yt?1+∑j=1pψjΔYt?j+εt(3)ΔYt=α+βΤ+?Yt-1+∑j=1pψjΔYt-j+εt(3) 式中:Yt是待检验之时间序列;α是常数项;T为时间趋势;p是时间滞后值;εt是随机误差项。 2 结果分析 2.1 年平均温度的变化趋势 对昆明站的月平均温度、最高温度、最低温度, 相对湿度、日照时数、风速等6个气象要素, 按1-12月份及年值分别采用Mann-kendall法进行趋势变化检验, 结果如表1。变化趋势是年平均温度、最高温度和最低温度都是普遍升高, 只有4、5、7、12月的最高温度呈下降趋势;但从1900年以来的变化趋势看, 昆明站的年平均温度呈现70年左右的“增-减-增”的周期性波动, 即1900