典型岩溶流域土地利用变化的水文响应.docxVIP

典型岩溶流域土地利用变化的水文响应 摘要：为了研究土地利用方式变化在典型岩溶流域产生的水文响应，构建了倒天河流域的SWAT模型，利用该模型探讨了不同土地利用情景下流域流量的变化规律和子流域产水量的变化规律及其调蓄能力分析，结果表明：参数率定后的SWAT模型在模拟倒天河流域月流量过程中表现良好，模型适用于该研究区；草地转为耕地会导致流域流量增加，调蓄能力减弱；林地向耕地、草地向耕地转化易改变子流域产水量，草地向林地的转化不易改变子流域产水量；林地向耕地、草地向耕地、草地向林地的土地利用转化均会导致子流域调蓄能力变弱，草地向耕地的转化对调蓄能力的改变尤为显著。研究成果对减轻岩溶区旱涝灾害，防治石漠化有重大意义。 中国西南岩溶地区干旱、洪涝灾害交替频发，生态环境脆弱大部分学者对土地利用变化引起的流域径流变化的研究是基于模型展开的。ZHANG目前，大多数研究都聚焦在径流量变化上，研究土地利用变化引起流域调蓄能力变化还不够深入，且对典型岩溶流域土地利用变化的水文响应研究较少。因此，本文基于SWAT模型从流域整体和子流域角度进行分析，评估典型岩溶流域土地利用变化引起流域流量、产水量、调蓄能力的变化，以期为决策者提供减轻岩溶旱涝、防治石漠化的技术支持。1 研究区概况研究区（105°5E～105°16E,27°15N～27°22N）位于贵州省毕节市七星关区（图1），属倒天河流域上游部分，区内主要分布村镇，倒天河下游位于毕节城区内。倒天河发源于毕节市黄泥大娄山脉南麓，属乌江水系支流白甫河（落脚河）上游。该区属亚热带季风湿润气候，冬无严寒，夏无酷暑，季风气候比较明显，降雨量较为充沛。1951-2017年多年平均降水量为1 134 mm，多年平均气温为12.6℃，多年平均相对湿度为81.5%。研究区地势西南高，东北低，海拔1 552～2 199 m，流域面积约105 km2 研究区数据收集与处理2.1 地理空间数据收集数据主要包括数字高程模型（DEM数据），土地利用/覆盖数据（Land use数据），土壤数据（Soil数据），详细信息见表1。2.2 气象水文数据气象数据主要包括研究区多年日值降水量、日值高低温、日值相对湿度、日值太阳辐射量、日值风速，水文数据主要包括多年逐日平均径流量。气象站、水文站分布见图1，详细信息介绍见表2。其中，毕节站不在研究区内，但数据缺失时亦可使用，但效果不佳。SWAT模型需要保证输入气象数据的起始日期和结束日期相同，其中缺测值可以用-99代替，因此，综合气象和水文数据的时间序列，本文采用1990—2006年的气象数据和流量数据来运行SWAT模型。2.3 搭建属性数据库数据属性库主要包括Land Use属性库、Soil属性库、坡度。(1)Land Use属性库。原始Land Use数据自带二级分类数据，按模型要求重分类后添加对应属性，见表3。(2)Soil属性库。模型需要手动添加Soil属性数据，根据原始数据自带的属性表，利用SPAW软件的SWC模块计算新的Soil属性数据，搭建了适用的土壤库。(3）坡度。研究区坡度按照自然间断点分级法（Jenks）分为3级，即0～8.6,8.6～18.2,18.2～72。3 模型原理及搭建3.1 SWAT模型基本原理SWAT模型由美国农业部（USDA）农业研究机构的工程师J G Arnold开发式中：?SW、P、Q、ET、DP、QR分别表示土壤含水量、降水、地表径流、实际蒸散发、深层下渗和浅层回归流。基于Arc GIS平台的Arc SWAT按照输入的DEM数据自动提取河网，生成子流域，再按照不同的Land Use数据、Soil数据、坡地将子流域划分为HRU（水文响应单元），本文划分阈值为5%、10%、10%，通过输入降水量、温度、太阳辐射量、相对湿度、风速等信息，计算每个HRU的产汇流情况。3.2 模型的校准与率定通过SWAT-CUP软件式中：Q表示流量，m目标函数的值越接近1，表示模拟值与实际值差别越小，模拟效果越好。目标函数置信区间见表4。3.3 模拟的结果与评价为模拟真实土地利用情况，1990-1997年使用1995年土地利用，1998-2002年使用2000年土地利用，2003-2006年使用2005年土地利用。以月步长运转模型，将1991-1992年作为预热期，1993-1997年作为校准期，1998-2002年作为验证期，2003-2006年作为验证期。(1）流域划分结果。模型划定的研究区面积为104.970km(2）参数率定结果。根据大量文献及学者建议，初步选取12个相关参数进行率定，在1 000次模拟后，选择较敏感的8个参数作为选定参数，利用SWAT-C