氣候及下墊面變化對(duì)嘉陵江流域徑流與輸沙的影響

胡云華，馮精金，王銘烽，田風(fēng)霞，賀秀斌?

(1.四川省第三測繪工程院，610500，成都；2.華中農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，430070，武漢; 3.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所，610041，成都)

?

氣候及下墊面變化對(duì)嘉陵江流域徑流與輸沙的影響

胡云華1，馮精金2，王銘烽3，田風(fēng)霞3，賀秀斌3?

(1.四川省第三測繪工程院，610500，成都；2.華中農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，430070，武漢; 3.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所，610041，成都)

為定量評(píng)價(jià)氣候、下墊面2種驅(qū)動(dòng)因子的變化對(duì)嘉陵江流域產(chǎn)流輸沙變化的貢獻(xiàn)率，采用SWAT模型，對(duì)嘉陵江流域1975年以來的徑流和輸沙過程進(jìn)行模擬，通過模型控制變量的方法，分析嘉陵江流域徑流和泥沙變化的驅(qū)動(dòng)因子。結(jié)果表明：1988—2010年的23年間，由于下墊面變化，導(dǎo)致嘉陵江流域徑流總量比1956—2010年的平均徑流總量減少了1 718.54億m3，輸沙總量減少了14.62億t，下墊面變化對(duì)徑流量減少的貢獻(xiàn)率為76.19%，對(duì)輸沙量減少的貢獻(xiàn)率高達(dá)91.8%；因此，說明下墊面變化是造成嘉陵江流域徑流量和輸沙量變化的主要原因。

SWAT模型； 水沙變化； 氣候變化； 嘉陵江

嘉陵江流域作為長江上游地區(qū)重要的徑流和泥沙來源之一，具有水土流失面廣量大、類型多樣、受人類活動(dòng)影響大等多種特點(diǎn)。根據(jù)多年觀測資料發(fā)現(xiàn)[1-2]，近幾十年來，嘉陵江流域輸沙量明顯減少，1988—2010年的23年間，嘉陵江流域年徑流量同1956—2010年多年平均徑流量相比，累計(jì)減少2 255.16億m3，占多年平均徑流量累積量的15.08%；輸沙量同1956—2010年多年平均輸沙量相比，總共減少15.85億t，占多年平均輸沙量23年累積量的62.65%。嘉陵江流域作為三峽庫尾直接的徑流和泥沙入口，流域水沙變化將對(duì)三峽庫區(qū)水沙運(yùn)動(dòng)以及泥沙沖淤變化等產(chǎn)生重要影響，引發(fā)了社會(huì)的廣泛關(guān)注[3-5]。

張信寶等[6]、范建容等[7]和許全喜等[8]曾利用水文統(tǒng)計(jì)數(shù)據(jù)、遙感和GIS技術(shù)、典型調(diào)查與分析等技術(shù)方法，分析嘉陵江流域水沙變化趨勢及原因；但由于調(diào)查和分析的局限性，研究得出的結(jié)論各有不同，主要原因在于流域出口水文站監(jiān)測數(shù)據(jù)的統(tǒng)計(jì)分析，無法很好的區(qū)分水沙來源，而且，同期嘉陵江流域的降雨量也有所減少，很難定量評(píng)價(jià)由于氣候變化導(dǎo)致的流域水沙變化[9]。

SWAT模型是由美國農(nóng)業(yè)部農(nóng)業(yè)研究所于20世紀(jì)90年代中后期開發(fā)的具有物理機(jī)制、以日為步長運(yùn)行的長時(shí)段流域分布式水文模型[10]。很多學(xué)者將SWAT模型應(yīng)用于國內(nèi)小流域徑流、泥沙模擬，都取得較好的效果[11-13]，使用SWAT模型參數(shù)變量控制的方法，可以對(duì)不同驅(qū)動(dòng)因素、對(duì)流域徑流和泥沙變化的貢獻(xiàn)進(jìn)行定量估算；因此，本研究通過SWAT模型，對(duì)嘉陵江流域徑流和泥沙變化過程進(jìn)行模擬分析，以分別量化近幾十年來，氣候和下墊面變化對(duì)嘉陵江流域徑流和輸沙變化的貢獻(xiàn)率。

1　研究區(qū)概況

嘉陵江發(fā)源于陜西秦嶺南麓(E 102°27′38″～109°1′16.9″，N 29°18′3″～34°32′26″)，干流全長1 119 km，流經(jīng)陜西、甘肅、四川3省，流域面積約16萬km2。嘉陵江流域可以分為渠江流域、涪江流域和嘉陵江干流3大子流域，其基本情況見表1。

表1　渠江、涪江和嘉陵江干流3個(gè)子流域基本情況

2　研究方法

2.1數(shù)據(jù)處理與方法

2.1.1流域邊界的提取和水系及子流域的劃分流域邊界提取和水系及子流域的劃分?jǐn)?shù)據(jù)來源于SRTM的90 m分辨率的DEM數(shù)據(jù)，利用ArcSWAT 2012自帶的Watershed Delineation工具自動(dòng)進(jìn)行，為了便于進(jìn)行模型的校準(zhǔn)和驗(yàn)證，選取七里沱、靜邊和東林水文站分別作為流域內(nèi)的校準(zhǔn)和驗(yàn)證水文站，選取羅渡溪、武勝和小河壩水文站分別作為流域出口的校準(zhǔn)和驗(yàn)證水文站。共將渠江流域劃分為47個(gè)子流域，嘉陵江干流劃分為92個(gè)子流域，涪江劃分為26個(gè)子流域。

2.1.2土地利用數(shù)據(jù)的準(zhǔn)備本研究使用的流域1986年土地利用類型圖，來源于中國1∶10萬土地利用數(shù)據(jù)，源自中國科學(xué)院“八五”重大應(yīng)用項(xiàng)目“全國資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究”。整個(gè)嘉陵江流域共包含6種一級(jí)類，22種二級(jí)類土地利用類型，為減小模型的復(fù)雜度，SWAT模型建議用戶土地利用類型不超過10種，本研究將按土地利用一級(jí)分類，建立模型土地利用數(shù)據(jù)圖層，流域共包括6大地類：耕地、林地、草地、水域、居民地和荒地。

2.1.3土壤數(shù)據(jù)的處理SWAT模型中土壤數(shù)據(jù)是主要的輸入?yún)?shù)之一，本研究使用的土壤物理屬性，來源于糧農(nóng)組織FAO提供的中國1∶100萬土壤類型圖，該數(shù)據(jù)粒徑分類采用USDA 簡化的美制標(biāo)準(zhǔn)，與SWAT 模型要求的土壤粒徑級(jí)配標(biāo)準(zhǔn)相符。而SOL_BD(土壤濕密度)、SOL_AWC(有效持水量)、SOL_K(飽和導(dǎo)水率)3個(gè)變量由SPAW軟件計(jì)算得到。

2.1.4坡度數(shù)據(jù)的處理坡度數(shù)據(jù)以SRTM 90 m分辨率的DEM數(shù)據(jù)計(jì)算而來，在SWAT模型中使用Multiple Slope坡度分類，模型支持的是百分比坡度，為了減小模型的計(jì)算量，參照模型使用手冊，執(zhí)行二級(jí)坡度分類，設(shè)置“slope class 1”的上限為40%，對(duì)應(yīng)實(shí)際坡度約25°左右，“slope class 2”的上限為9 999%，對(duì)應(yīng)實(shí)際坡度90°。

2.1.5水文響應(yīng)單元(HRU)的劃分水文響應(yīng)單元是流域內(nèi)具有相同水文特性的最小水文單元，稱為HRU，是SWAT模型計(jì)算的最小單元。根據(jù)之前輸入的土地利用類型數(shù)據(jù)、土壤數(shù)據(jù)和坡度數(shù)據(jù)，在劃分的子流域的前提下，根據(jù)土壤、坡度和植被等因素，將整個(gè)嘉陵江流域共劃分出945個(gè)水文響應(yīng)單元。

2.1.6氣象觀測數(shù)據(jù)的處理SWAT要求用戶必須輸入模擬期內(nèi)的日觀測降雨和氣溫?cái)?shù)據(jù)，由于實(shí)際觀測過程中人為失誤或者系統(tǒng)誤差，實(shí)際觀測數(shù)據(jù)中有大量的缺測值；因此，在模型建立過程中，需要利用WXGEN天氣發(fā)生器模擬缺測值。WXGEN天氣發(fā)生器進(jìn)行天氣模擬時(shí)，需要輸入至少30年的天氣統(tǒng)計(jì)參數(shù)，主要的統(tǒng)計(jì)參數(shù)有月平均最高氣溫、月平均最低氣溫、最高氣溫標(biāo)準(zhǔn)偏差、月平均降雨量、降雨量標(biāo)準(zhǔn)偏差、月內(nèi)干日時(shí)間(d)、露點(diǎn)溫度和月平均太陽輻射量等。筆者研究獲取了嘉陵江流域及其周邊24個(gè)氣象站(圖1)，1975—2010年36年的日氣象觀測數(shù)據(jù)，用作模型計(jì)算及天氣發(fā)生器數(shù)據(jù)庫的構(gòu)建。

2.1.7水文和泥沙數(shù)據(jù)的處理模型校準(zhǔn)期1975—1982年和驗(yàn)證期1983—1987年的徑流和泥沙數(shù)據(jù)來源于1975—1987年的“長江流域水文資料：嘉陵江區(qū)”，整理得到渠江流域內(nèi)的東林、七里沱、渡邊3個(gè)子流域和渠江流域下游的羅渡溪、嘉陵江干流下游的武圣、涪江下游的小河壩6個(gè)水文觀測站(圖1)，1975—1987年逐月的徑流和泥沙觀測數(shù)據(jù)。

在模型的模擬期，收集到羅渡溪、小河壩、武圣站1988—2000年的月徑流數(shù)據(jù)。羅渡溪水文站控制流域面積3萬6 936.71 km2，小河壩水文站控制流域面積2萬5 634.03 km2，武勝水文站控制流域面積7萬9 079.98 km2，總控制區(qū)面積14.17萬km2，占整個(gè)嘉陵江流域面積的91.55%。北碚水文站的控制流域面積14萬9 165 km2，占整個(gè)嘉陵江流域面積的96.4%。本研究中，用北碚水文站的觀測數(shù)據(jù)代表整個(gè)嘉陵江的徑流和泥沙狀況。北碚站的徑流和泥沙主要來源于渠江、涪江和嘉陵江干流流域，由于沒有收集到北碚站的1975—2000年的觀測數(shù)據(jù)，將武勝、羅渡溪和小河壩3個(gè)站的數(shù)據(jù)相加，以替代北碚水文站的月徑流和泥沙數(shù)據(jù)，和模擬結(jié)果進(jìn)行對(duì)比，以判斷擬合效果。

北碚站2001—2010年的月徑流和泥沙資料來源于2001—2010年的“長江泥沙公報(bào)”，由于沒有在北碚水文站整個(gè)控制區(qū)內(nèi)建立模型，將武勝、羅渡溪和小河壩3個(gè)站2001—2010的模擬數(shù)據(jù)相加，同北碚水文站的月徑流和泥沙觀測數(shù)據(jù)對(duì)比，以判斷模型擬合效果。

圖1　研究區(qū)氣象站和水文站分布圖Fig.1　Spatial distribution of meteorological and hydrologic stations in the studied area

2.2SWAT運(yùn)行、校準(zhǔn)和驗(yàn)證

2.2.1SWAT運(yùn)行將模型模擬所需參數(shù)輸入到ArcSWAT當(dāng)中，設(shè)置1973—1974年2年作為模型的預(yù)熱期，運(yùn)行模型，分別模擬渠江流域、涪江流域和嘉陵江干流1975—1982年的徑流和泥沙值。此時(shí)的模擬結(jié)果是使用模型默認(rèn)參數(shù)生成的，模型內(nèi)部參數(shù)大多根據(jù)美國的環(huán)境條件設(shè)置。通過對(duì)比渠江流域的模擬結(jié)果可得，渠江流域出口羅渡溪水文站SWAT模擬的1975—1982徑流量擬合優(yōu)度判斷系數(shù)R2為0.87，Nash-Sutcliffe(NS)系數(shù)僅為0.44，模擬的輸沙量擬合優(yōu)度判斷系數(shù)R2為0.36，NS系數(shù)僅為-0.61。模擬效果比較差，為提高模型模擬的精度，需要利用實(shí)際觀測數(shù)據(jù)對(duì)模型參數(shù)進(jìn)行校準(zhǔn)。

2.2.2SWAT-CUP校準(zhǔn)使用SWAT-CUP2012作為模型校準(zhǔn)的軟件，對(duì)徑流進(jìn)行率定時(shí)，選擇對(duì)徑流較為敏感的13個(gè)參數(shù)，對(duì)泥沙進(jìn)行率定時(shí)，選擇較敏感的21個(gè)參數(shù)。在流域出口武勝、羅渡溪、小河壩和流域內(nèi)東林、七里沱、靜邊6個(gè)水文觀測站，采用從上游到下游的順序依次率定，每輪率定設(shè)置運(yùn)行次數(shù)500次，每輪率定完成后，采用推薦的參數(shù)范圍進(jìn)行下一輪率定，直到R2系數(shù)和NS系數(shù)穩(wěn)定為止。經(jīng)過SWAT-CUP參數(shù)的率定，SWAT模擬的結(jié)果和實(shí)際觀測數(shù)據(jù)逐漸接近，各水文站月徑流模擬值和實(shí)際觀測值的R2系數(shù)和NS系數(shù)均在0.8以上(圖2)，月輸沙量的模擬值和和實(shí)際觀測值的R2系數(shù)和NS系數(shù)均在0.75以上(圖3)，證明模擬結(jié)果可以較好地?cái)M合校準(zhǔn)期的徑流和輸沙過程。

圖2　6個(gè)水文站月徑流量校準(zhǔn)期的效果圖Fig.2　Calibration results of monthly runoff in the 6 hydrologic stations

圖3　3個(gè)水文站月輸沙量校準(zhǔn)期的效果圖Fig.3　Calibration results of monthly sediment in the 3 hydrologic stations

2.2.3SWAT-CUP驗(yàn)證經(jīng)過SWAT-CUP率定后，將率定得到的最佳參數(shù)帶入到SWAT模型當(dāng)中，替換原來模型當(dāng)中對(duì)應(yīng)子流域內(nèi)的參數(shù)，重新運(yùn)行模型，生成渠江流域、涪江流域和嘉陵江干流1983—1987年的徑流量和泥沙輸移量，并與實(shí)際觀測數(shù)據(jù)進(jìn)行對(duì)比，計(jì)算模型在驗(yàn)證期徑流模擬的精度，結(jié)果如圖4和圖5所示。在驗(yàn)證期，6個(gè)流域水文站月徑流量模擬結(jié)果和實(shí)際觀測結(jié)果的R2系數(shù)和NS系數(shù)均在0.75以上，3個(gè)流域水文站月輸沙量模擬結(jié)果和實(shí)際觀測結(jié)果的R2系數(shù)和NS系數(shù)均在0.5以上，模擬結(jié)果較好。說明經(jīng)過校準(zhǔn)的嘉陵江流域SWAT模型，可以較好地?cái)M合驗(yàn)證期流域的徑流和輸沙過程，所建立的嘉陵江流域SWAT模型，在模型的校準(zhǔn)期和驗(yàn)證期比較可靠。

圖4　6個(gè)水文站月徑流量驗(yàn)證期的效果圖Fig.4　Verification results of monthly runoff in the 6 hydrologic stations

圖5　3個(gè)水文站月輸沙量驗(yàn)證期的效果圖Fig.5　Verification results of monthly sediment in the 3 hydrologic stations

由于模型建立在1975—1982年的基礎(chǔ)數(shù)據(jù)之上，經(jīng)過校準(zhǔn)后，從驗(yàn)證結(jié)果可以看出，模型內(nèi)部參數(shù)已經(jīng)基本符合當(dāng)時(shí)流域基本情況，在模型其他參數(shù)不變的情況下，僅輸入氣象觀測數(shù)據(jù)，可以預(yù)測在地表下墊面不變的情況下，嘉陵江流域近30年來的水沙變化過程。將模擬結(jié)果和實(shí)際觀測數(shù)據(jù)對(duì)比，可以估算地表下墊面變化所造成的流域徑流和輸沙變化量，將實(shí)際觀測的變化總量和由地表下墊面變化所造成的流域徑流和輸沙變化量對(duì)比，可以估算因氣候變化導(dǎo)致的流域徑流和輸沙變化量。

3　結(jié)果與分析

3.1地表下墊面變化驅(qū)動(dòng)下的嘉陵江流域徑流和泥沙變化分析

從圖6和圖7可知：在SWAT模型的校準(zhǔn)期和驗(yàn)證期，模型可以較好地模擬嘉陵江流域的徑流和輸沙過程；但隨著模擬時(shí)間的推移，擬合結(jié)果的誤差逐步增大。SWAT模擬的效果由模型參數(shù)決定，由于模型的參數(shù)是用1975—1982年的實(shí)測數(shù)據(jù)進(jìn)行校準(zhǔn)的，模型輸入的土地利用數(shù)據(jù)也是80年代的，隨著時(shí)間的推移，模型參數(shù)中只改變氣象參數(shù)。但近30年來，嘉陵江流域的下墊面已經(jīng)發(fā)生很大的變化，而模型下墊面參數(shù)并未修改，導(dǎo)致模擬結(jié)果和實(shí)際觀測結(jié)果誤差越來越大。在模型參數(shù)不改變的情況下，可以認(rèn)為模擬結(jié)果是地表下墊面未發(fā)生改變時(shí)的嘉陵江流域徑流和泥沙數(shù)以過程。將模擬結(jié)果和實(shí)測結(jié)果進(jìn)行對(duì)比分析，可以得出近30年來，嘉陵江流域地表下墊面變化所造成的嘉陵江流域徑流和泥沙變化量。從統(tǒng)計(jì)結(jié)果得出，僅在地表下墊面變化的影響下，1988—2010年的23年來，嘉陵江流域徑流量同多年平均徑流量相比，累計(jì)減少1 718.54億m3，占實(shí)際觀測徑流減少量的76.19%，占多年平均徑流量23年累積量的11.50%。輸沙量同多年平均輸沙量相比，總共減少14.62億t，占實(shí)際觀測輸沙量減少量的91.8%，占多年平均輸沙量23年累積量的57.79%。

圖6　嘉陵江流域?qū)嶋H觀測和模擬的年徑流量變化過程Fig.6　Change process of simulated and observed annual runoff in Jialing River Basin

圖7　嘉陵江流域?qū)嶋H觀測和模擬的年輸沙量變化過程Fig.7　Change process of simulated and observed annual sediment discharge in Jialing River Basin

3.2氣候變化驅(qū)動(dòng)下的嘉陵江流域徑流和泥沙變化分析

圖8描繪出流域年降雨量的變化過程。從統(tǒng)計(jì)結(jié)果可得1988—2010年的23年來，嘉陵江流域年降雨量同多年平均降雨量相比累計(jì)減少895 mm，占嘉陵江流域多年平均降雨量23年累積量的4.29%。

圖8　嘉陵江流域年降雨量變化過程Fig.8　Change process of annual precipitation in Jialing River Basin

在模型其他參數(shù)不變的情況下，模型模擬的結(jié)果是只在氣象因素(降雨量、最高氣溫、最低氣溫等)影響下的徑流和輸沙量。模型模擬的年平均徑流量和年輸沙量如圖9和圖10所示?？梢钥闯觯耗M的年徑流量和輸沙量變化過程與流域的年降雨量變化過程高度正相關(guān)，降雨的變化會(huì)引發(fā)流域徑流量和輸沙量的變化。僅在氣象條件的影響下，1988—2010年的23年來，嘉陵江流域徑流量同多年平均徑流量相比，累計(jì)減少了537億m3，占實(shí)際觀測徑流減少量的23.81%，占多年平均徑流量23年累積量的3.59%。輸沙量同多年平均輸沙量相比，總共減少1.3億t，占實(shí)際觀測輸沙減少量的8.2%，占多年平均輸沙量23年累積量的5.14%。

圖9　嘉陵江流域模擬的年徑流量變化過程Fig.9　Change process of simulated annual runoff in Jialing River Basin

圖10　嘉陵江流域模擬的年輸沙量變化過程Fig.10　Change process of simulated annual sediment discharge in Jialing River Basin

4　結(jié)論與討論

1) 1988—2010年的23年來，嘉陵江流域年降雨量同多年平均降雨量相比累計(jì)減少895 mm，在地表下墊面保持不變的情況下，單純由氣象條件導(dǎo)致的嘉陵江流域徑流總量比1956—2010年的平均徑流量減少約3.59%，輸沙總量減少約5.14%，氣象條件變化對(duì)嘉陵江流域徑流減少的貢獻(xiàn)率為23.81%，對(duì)輸沙量減少的貢獻(xiàn)率僅為8.2%。由此可見，單純氣象條件變化僅導(dǎo)致嘉陵江流域徑流和輸沙量略微減少。

2) 排除氣象條件造成的徑流和泥沙變化，由于下墊面變化，導(dǎo)致嘉陵江流域徑流總量比1956—2010年的平均徑流總量減少約11.5%，輸沙總量減少約57.79%，下墊面變化對(duì)嘉陵江流域徑流量減少的貢獻(xiàn)率為76.19%，對(duì)輸沙量減少的貢獻(xiàn)率高達(dá)91.8%。由此可見，下墊面變化是造成嘉陵江流域徑流量和輸沙量變化的主要原因；但目前的模型精度，還無法較準(zhǔn)確地區(qū)分各下墊面變化驅(qū)動(dòng)因子對(duì)流域徑流和輸沙過程變化的貢獻(xiàn)率。地表下墊面因素中，水土保持措施、土地利用變化、水庫水電站等大型水利設(shè)施的修建等因素對(duì)流域徑流和輸沙過程的影響，將是需要進(jìn)一步研究的問題。

[1]長江水利委員會(huì)長江流域水土保持監(jiān)測中心站.嘉陵江流域水土保持減沙效益研究報(bào)告[R].武漢:長江水利委員會(huì),2006:1-64.

Yangtze River Scientific Research Institute,Changjiang Water Resources Commission.Benefit research report of soil-water conservation on water and sand reduction in Jialing River basin[R].Wuhan:Changjiang Water Resources Commission,2006:1-64.(in Chinese)

[2]陳顯維.嘉陵江流域水庫群攔沙量估算及攔沙效應(yīng)分析[J].水文,1992,(4):34.

Chen Xianwei.Weight estimation and benefit analysis of sediment trapped by reservoir group in Jialing River basin[J].Hydrology,1992,(4):34.(in Chinese)

[3]許炯心.人類活動(dòng)和降水變化對(duì)嘉陵江流域侵蝕產(chǎn)沙的影響[J].地理科學(xué),2006,26(4):432.

Xu Jiongxin.Effect of human activities and precipitation on sediment yield of the Jialing River[J].Geographical Science,2006,26(4):432.(in Chinese)

[4]丁文峰,張平倉,任洪玉.近50年來嘉陵江流域徑流泥沙演變規(guī)律及驅(qū)動(dòng)因素定量分析[J].長江科學(xué)院院報(bào),2008,25(3):23.

Ding Wenfeng,Zhang Pingcang,Ren Hongyu.Quantitative analysis on evolution characteristics and driving factors of annual runoff and sediment transportation changes for Jialing River[J].Journal of Yangtze River Scientific Research Institute,2008,25(3):23.(in Chinese)

[5]劉惠英,高強(qiáng).三峽庫區(qū)龍河流域水沙關(guān)系驅(qū)動(dòng)因子貢獻(xiàn)率分析[J].中國水土保持科學(xué),2015,13(5):1.

Liu Huiying,Gao Qiang.Contribution rate of driving factors on the rainfall-sediment relationship of Longhe River watershed in the Three Gorges Reservoir Region [J].Science of Soil and Water Conservation,2015,13(5):1.(in Chinese)

[6]張信寶,文安邦.長江上游干流和支流河流泥沙近期變化及其原因[J].水利學(xué)報(bào),2002,27(4):56.

Zhang Xinbao,Wen Anbang.Variations of sediment in upper stream of Yangtze River and its tributary[J].Journal of Hydraulic Engineering,2002,27(4):56.(in Chinese)

[7]范建容,鐘祥浩,劉淑珍.嘉陵江中下游典型流域土壤侵蝕與泥沙輸移遙感監(jiān)測[J].中國科學(xué)(E輯:科學(xué)技術(shù)),2003,33(S1):157.

Fan Jianrong,Zhong Xianghao,Liu Shuzhen.Remote sensing of soil erosion and sediment transport in the typical area of middle and lower reach of Jialing River[J].Science in China Series E:Technological Sciences,2003,33(S1):157.(in Chinese)

[8]許全喜,陳松生,熊明,等.嘉陵江流域水沙變化特性及原因分析[J].泥沙研究,2008,(2):1.

Xu Quanxi,Chen Songsheng,Xiong Ming,et al.Analysis on runoff and sediment characteristics and affecting factors of the Jialing River basin[J].Journal of Sediment Research,2008,(2):1.(in Chinese)

[9]楊麗紅,周正朝.最小二乘與分位數(shù)回歸方法在河流水沙研究中的應(yīng)用：以陜北延河流域?yàn)槔齕J].中國水土保持科學(xué),2014,12(5):45.

Yang Lihong,Zhou Zhengchao.Comparison of ordinary least square regression and quantile regression in the study of hydrological elements of Yanhe River Basin[J].Science of Soil and Water Conservation,2014,12(5):45.(in Chinese)

[10] Arnold J G,Srinivasan R,Muttiah R S,et al.Large-area hydrologic modeling and assessment:Part I.Model development[J].Journal of the American Water Resources Association,1998,34(1):73.

[11] 郭軍庭,張志強(qiáng),王盛萍,等.應(yīng)用SWAT模型研究潮河流域土地利用和氣候變化對(duì)徑流的影響[J].生態(tài)學(xué)報(bào),2014,34(6):1559.

Guo Junting,Zhang Zhiqiang,Wang Shengping,et al.Appling SWAT model to explore the impact of changes in land use and climate on the streamflow in a watershed of northern China[J].Acta Ecologica Sinica,2014,34(6):1559.(in Chinese)

[12] 韓博,周麗麗,范昊明,等.蒲河流域氮污染負(fù)荷模擬及時(shí)空分布[J].中國水土保持科學(xué),2013,11(2):90.

Han Bo,Zhou Lili,Fan Haoming,et al.Nitrogen pollution load simulation of the Puhe River basin and the spatial and temporal distribution[J].Science of Soil and Water Conservation,2013,11(2):90.(in Chinese)

[13] 郝芳華，陳利群，劉昌明，等.土地利用變化對(duì)產(chǎn)流和產(chǎn)沙的影響分析[J].水土保持學(xué)報(bào),2004,18(3):5.

Hao Fanghua,Chen Liqun,Liu Changming,et al.Impact of land use change on runoff and sediment yield[J].Journal of Soil and Water Conservation,2004,18(3):5.(in Chinese)

Influences of climate and land surface change on runoff and sediment in Jialing River Basin

Hu Yunhua1,Feng Jingjin2,Wang Mingfeng3,Tian Fengxia3,He Xiubin3

(1.The Third Surveying and Mapping Engineering Institute in Sichuan Province,610500,Chengdu,China; 2.College of Resource and Environment, Huazhong Agricultural University,430070,Wuhan,China; 3.Institute of Mountain Hazards and Environment,Chinese Academy of Sciences, 610041,Chengdu,China)

[Background] The Jialing River Basin is one of the important sources of runoff and sediment in the upper reaches of the Yangtze River Basin.Investigation on the mechanism of the influence of climate and underlying surface conditions on runoff and sediment discharge reveal the reason for the change of runoff and sediment discharge in Jialing River Basin in recent decades,therefore,analyzing the influence contribution rates from 2 factors of climate and underlying surface conditions on the changes of runoff and sediment has the important significance in the Jialing River Basin,and it also can be used as certain significant guidance to control the basin sediment variability of soil and water conservation work.[Methods] SWAT model was adopted to simulate the process of runoff and sediment in the Jialing River Basin since 1975,and was calibrated by processing the basic data (including land use,DEM,soil data and meteorological data) and using the sediment data of the hydrologic stations,furthermore the driving factors resulting in the changes of runoff and sediment as well as their contribution rates to the changes of runoff and sediment were analyzed by the method of model control variables.[Results] The study showed that:during the 23 years from 1988 to 2010,1) the total runoff in Jialing River Basin was 3.59% less than the average runoff,and the total amount of sediment decreased by 5.14% and only resulted from the change of climate condition.The contribution rate of climate changes on the decrease of runoff was 23.81% and on the decrease of sediment was 8.2%.2) Compared with the average total amount of runoff and sediment discharge,the total runoff in Jialing River Basin had been a decrease of 171.854 billion cubic meters and the sediment discharge had been a decrease of 1 462 million tons.The contribution rate of underlying surface changes on the decrease of runoff was 76.19% and on the decrease of sediment was 91.8%.[Conclusions] 1) The influence of weather change on runoff and sediment in Jialing River Basin is little.2) The changes of runoff and sediment in Jialing River Basin are mainly caused by the changes of underlying surface conditions.Although this paper shows the underlying surface is the main reason on the change of runoff and sediment in Jialing River Basin,by the model precision it is infeasible to accurately distinguish the contribution rate of varied driving factors of underlying surface to the change processes of runoff and sediment.The underlying surface factors include soil and water conservation measures,land use change,reservoir hydropower station and the construction of large water conservancy facilities,and their influences on the runoff and sediment need further studying.

SWAT model; runoff and sediment change; climate change; Jialing River

2015-05-15

2016-04-25

項(xiàng)目名稱：國家自然科學(xué)基金“三峽庫區(qū)紫色土柑橘園土質(zhì)道路路面侵蝕的研究”(41201273)；中國科學(xué)院資助課題“氣候和人類影響下的坡地徑流與侵蝕過程”(KZZD-EW-TZ-06-03)；長江水利委員會(huì)資助課題“長江流域水土保持防治技術(shù)調(diào)研和總結(jié)”(Y4R1-450)

胡云華(1988—)，男，碩士，助理工程師。主要研究方向：地理信息系統(tǒng)。E-mail:578652097@qq.com

簡介：賀秀斌(1967—)，男，博士，研究員。主要研究方向：土壤侵蝕與水土保持。E-mail:xiubinh@imde.ac.cn

S157.1

A

1672-3007(2016)04-0075-09

10.16843/j.sswc.2016.04.010