气候变化及人类活动对延河流域泥沙影响定量探析.doc

气候变化及人类活动对延河流域泥沙影响定量探析 　　摘要：对延河4个子流域1982―2010年的产沙量、降雨量采用Mann-Kendall非参数检验法、累积距平法分析其变化趋势及产沙量突变点，分析汛期和全年产沙量与径流量的关系，使用水文分析法分离了人类活动和气候变化对流域产沙量的影响，最终结合景观指数阐述植被与泥沙的关系。结果表明：①该流域产沙比较集中，主要在6―9月份。研究时段内研究区的产沙量显著下降（P 　　目前关于气候变化和人类活动对水文水资源影响的研究大多集中在较大尺度的区域，而对中小尺度流域造成的影响研究较少。有关延河流域的景观格局变化，已有研究者进行了整体研究[15]，但并未细致到以站点为单位的子流域，而各子流域间存在的差异会影响管理决策的制定[16]。本研究基于水文分析，对延河4个子流域的景观格局进行分析，探讨各景观类型对产沙量的影响 1研究区概况 延河是黄河右岸、中游区上段的河口镇至龙门段的一级支流，位于东经108°45′～110°28′、北纬36°23′～37°17′之间，发源于陕西省靖边县东南天赐湾乡周山，由西北向东南流经志丹、安塞、宝塔、延长四县区，在延长县南河沟乡凉水岸附近汇入黄河[17]。其干流总长390.59 km，流域面积5 891 km2，平均坡度为16.4°，平均海拔为1 318 m，河网密度约为4.7 km/km2。流域多年平均降水量为435.64 mm，年平均气温9℃，主要灾害类型为旱灾、霜冻、冰雹及暴雨。延河流域地势西北高，东南低，主要属于黄土丘陵沟壑区，占流域总面积的 94.6%，其中：上游属梁峁丘陵沟壑区，梁多峁少，河床比降大，侵蚀强烈；中部属峁状丘陵沟壑区，梁窄而短，峁小而圆，侵蚀小于上游[18]。流域内设有20多处分布比较均匀的雨量站，水文站有上游杏河流域的杏河站（379 km2）、安塞站（1 334 km2）、延水流域延安站（3 208 km2）、枣园站（719 km2）及下游控制断面甘谷驿站[6 150 km2（1957―1964年）、5 898 km2（1965―1970年）、5 891 km2（1971―1989年）]5个。延河流域水土流失非常严重，是黄河泥沙的主要来源之一，含沙浓度高，泥沙年内分布集中，6―9月产沙量占全年产沙量的96.0%以上 本次研究所选区域包括杏河、安塞、枣园及甘谷驿站，属于延河中上游的范围 2数据来源与处理方法 2.1数据来源 根据流域气象站点分布，选用资料较完整且相邻的20个气象站点1982―2010年的年降雨资料，以保证空间布点的均匀性和流域界附近等值线的正确延伸。对于数据缺失的站点及年份采用反距离权重（IDW）插值的方式进行延展。气象资料来源于中国气象科学数据共享服务网；选用中国科学院提供的4期（1995、2000、2007、2010年）延河流域土地利用数据；流域径流量以及泥沙量选用由黄河及长江水利委员会（水资源公报）提供的各水文观测站1981―2010年逐年实测资料 2.2研究方法 2.2.1Mann-Kendall检验Mann-Kendall非参数趋势检验法可用于检测要素序列变化趋势[19，20]，其优点是统计测试的样本不需要服从一定的分布，也不受少数异常值的干扰，适用于非正态分布的水文气象等数据。在Mann-Kendall检验中，原假设H0为时间序列数据（x1，x2，…，xn），是n个独立的、随机分布的样本；对于所有的k，j≤n，且k≠j，xk和xj的分布式是不相同的，检验的统计变量S计算公式如下： 2.2.2距平及累计距平法在有序数据系列存在趋势性变化的前提下，采用累计距平法检验其跃变时间[21]。距平是表示气候变量偏离正常情况最常用的量，一组数据的某一个数xi与x之间的差就是该数据的距平值，即xi-x，气候变量的一组数据x1，x2，…，xn与其均值的差就构成了距平序列x1-x，x2-x，…，xi-x 累积距平是一种十分常用的、根据曲线直观判断变化趋势的方法。对于某一序列x某一时刻t的累积距平表示为： 2.2.3水文分析法为区分降水和人类活动对流域水沙变化的影响，本研究借鉴了计算水土保持减水减沙效益的一种方法――水文分析法，简称水文法[22]。河川径流量与输沙量受降水和下垫面条件控制，在自然状态下，它们之间具有统计相关关系，即某一特定流域，若下垫面条件保持不变，在一定的降水条件下产生的径流量和输沙量基本是一定的；而如果下垫面条件发生了变化，那么在相同的降雨条件下产生的径流量和输沙量就可能不同。根据上述原理，利用基准期（1980―1989年）的实测水文气象资料，通过多元回归分析，建立累积降雨与累积泥沙关系的回归方程[23]： 3结果与分析 3.1泥沙、水土保持措施与降