流域土地利用／覆被變化對水文過程的影響研究進展

史曉亮，李 穎，嚴(yán)登華，趙 凱

（1.中國科學(xué)院 東北地理與農(nóng)業(yè)生態(tài)研究所，長春130102；2.中國水利水電科學(xué)研究院 水資源研究所，北京100038；3.中國科學(xué)院大學(xué)，北京100049）

近年來，由于全球變化和人類活動的影響，地球上的水循環(huán)正在發(fā)生變化，世界不同國家或地區(qū)正在發(fā)生嚴(yán)重的水資源問題和危機。水問題已經(jīng)成為限制國家和區(qū)域可持續(xù)發(fā)展的關(guān)鍵性因子，水科學(xué)問題也成為國際地球科學(xué)發(fā)展中的一個重要方面［1］。

隨著水問題的日益突出，國際研究計劃以及世界各國的研究重點也在不斷變化。國際有關(guān)組織實施了一系列國際水科學(xué)計劃，如I HP、WCPR、GWSP等，其目的是從全球、區(qū)域和流域不同尺度探討環(huán)境變化（全球變化和人類活動影響）下的水循環(huán)及與之關(guān)聯(lián)的資源與環(huán)境問題。目前所面臨的水科學(xué)問題與人類非理性的開發(fā)活動密切相關(guān)，有研究表明，這些影響可能超過了氣候變化所造成的影響［2］。而土地利用／覆被變化（LUCC）代表了一種人為的“系統(tǒng)干擾”，最直接體現(xiàn)和反映了人類活動的影響水平，是直接或間接影響水文過程的第二個主要邊界條件［3］。在L UCC計劃結(jié)束之后，IGBP和I HDP聯(lián)合提出全球土地計劃（GLP），而從L UCC到GLP，土地利用／覆被變化分析一直被作為研究氣候變化和人類干擾之間復(fù)雜的相互作用所產(chǎn)生的生態(tài)環(huán)境效應(yīng)的突破點［4］。水文循環(huán)的生物圈（BAHC）作為IGBP的核心計劃之一，更加注重陸面生態(tài)—水文過程與空間格局的變化規(guī)律和受人類影響的關(guān)鍵問題。因此，變化環(huán)境下的水循環(huán)研究成為21世紀(jì)水科學(xué)研究的熱點，在目前乃至未來幾十年，LUCC和水文過程的相互作用研究將持續(xù)是全球生態(tài)與環(huán)境研究的熱點和前沿問題［5－6］。

1 流域土地利用／覆被變化對水文過程的影響特性

人類土地利用方式的改變而產(chǎn)生的土地覆被變化，深刻地改變了地表的性質(zhì)如粗糙度、反射率、植被冠層葉面積指數(shù)和影響水文通量的其它物理性質(zhì)，從而影響冠層截留、入滲和蒸散發(fā)等水文過程；另一方面，地表及其覆蓋的植被決定著太陽輻射在地表的分配，而土地覆被發(fā)生變化時則引起能量的重新分配，從而導(dǎo)致潛熱通量的變化（圖1）。因此，流域土地利用／覆被通過影響地表水熱通量，從而改變流域產(chǎn)匯流過程，改變水循環(huán)的空間格局，加速水循環(huán)要素時空分異的復(fù)雜性和不確定性，最終導(dǎo)致水資源供需關(guān)系發(fā)生變化。因此，土地利用／覆被變化是自然—人為作用共同導(dǎo)致的下墊面變化，也是影響水文過程的重要因素，對流域生態(tài)和社會經(jīng)濟發(fā)展等多方面具有顯著影響。

圖1 流域土地利用／覆被變化與地表水熱平衡關(guān)系

2 流域土地利用／覆被變化對水文過程影響的研究回顧

隨著對LUCC驅(qū)動力及其生態(tài)環(huán)境效應(yīng)研究的不斷深入，土地利用／覆被變化對水文過程的影響研究得到了很大的發(fā)展。根據(jù)國內(nèi)外研究中的研究方法和內(nèi)容的不同，其主要發(fā)展過程可以概括為3個階段（圖2）。相關(guān)研究最初主要集中于小流域森林變化的水文效應(yīng)，最早開始于1900年瑞士Emmental山區(qū)的兩個小流域?qū)Ρ仍囼?。在此之后，國?nèi)外許多學(xué)者采用試驗流域或者水文特征參數(shù)時間序列法開展了大量小流域森林變化、以及中小尺度土地利用／覆被變化對地表徑流的影響研究。1995年隨著LUCC研究計劃的制定與實施，LUCC的水文效應(yīng)研究受到更大程度的關(guān)注。21世紀(jì)以來，隨著基于DEM的分布式水文模型的逐步發(fā)展，結(jié)合遙感技術(shù)在獲取大尺度流域水文過程關(guān)鍵變量、水文模型的輸入和率定等方面的優(yōu)勢，大尺度流域土地利用／覆被變化對各水文參量的影響成為主要研究內(nèi)容；并且將土地利用變化模型與分布式水文模型相結(jié)合預(yù)測未來土地利用／覆被格局對水文過程的潛在影響成為近年來的研究重點。

圖2 流域土地利用／覆被變化對水文過程影響的研究回顧

3 流域土地利用／覆被變化對水文過程的影響研究的主要進展

3.1 流域土地利用／覆被變化對徑流的影響

Calder［21］認(rèn)為，在區(qū)域尺度上影響水文過程的主要土地利用變化是植被變化（如造林和毀林）、農(nóng)業(yè)開發(fā)活動以及城市化等。由于徑流能夠反映整個流域的生態(tài)狀況，也能用于預(yù)測未來土地利用／覆被潛在變化對水文水資源的影響［22］。因此目前土地利用／覆被變化對水文過程的影響研究多集中在對徑流的影響方面。

3.1.1 植被變化對徑流的影響 在土地利用／覆被變化的水文響應(yīng)研究中，森林的水文效應(yīng)倍受關(guān)注。上世紀(jì)80年代左右，劉昌明［9］針對森林對年徑流量的影響做了初步分析，黃秉維［23］也強調(diào)不要過分夸大森林的作用，應(yīng)正確認(rèn)識森林效應(yīng)。針對“森林植被的存在能否減少流域年徑流量和削減洪峰”等問題，國內(nèi)外出現(xiàn)了不同甚至相反的觀點。

一般認(rèn)為，在面積較小的集水區(qū)和流域，森林植被的存在，植被蓋度的增加會減少徑流量，并認(rèn)為其原因在于流域面積較小時，森林蒸騰大量水分起到主要作用［24］。但也有研究在大尺度流域得到類似的結(jié)論［25］，Jackson等［26］基于南半球600多個觀測點的數(shù)據(jù)研究認(rèn)為，為固碳而栽種樹木后，10 a之內(nèi)徑流量將大幅度減少，甚至導(dǎo)致溪流干涸。Cao等［14］、陳瑩等［18］、李佳等［19］、王盛萍等［27］也得出了類似的結(jié)論。其它許多研究在不同的流域，利用不同的研究方法同樣得到了植被減少將導(dǎo)致徑流量增加的結(jié)論［28－30］，這也從反面證明了森林覆蓋率增加會減少年徑流量的觀點。

但也有研究認(rèn)為，森林植被的存在會增加年徑流量［31］。郝芳華等［13］在黃河下游支流洛河上游盧氏水文站以上流域的研究表明，森林的存在增加了徑流量。王根緒等［32］的研究等表明，1967年以來，馬營河流域上游林草地大規(guī)模轉(zhuǎn)為耕地，使流域年均徑流量減少了28.12%。

一般認(rèn)為森林通過喬、灌、草以及枯枝落葉層的截持含蓄、大量蒸騰、土壤滲透和延緩融雪等過程，調(diào)節(jié)洪水發(fā)生過程，削減洪峰流量。但這一作用又與降雨強度與歷時、森林類型、林分結(jié)構(gòu)、土壤前期含水量等多種因素有關(guān)，隨著降水強度的增加和時間的延長，森林的削洪作用也將逐漸減弱［33］，連續(xù)性的大雨超過400 mm，森林與洪水的關(guān)系不復(fù)存在［34］。因此，森林的削洪作用有一定限度。

由此可見，由于區(qū)域氣候和地理特征的差異性，以及研究尺度與方法的不同，特定流域上森林植被變化對水文過程的影響較為復(fù)雜，難以取得一致性的結(jié)論。某一自然條件和不同尺度上得出的結(jié)論不具可比性，也不能簡單地進行外推。

3.1.2 城市化對徑流的影響 城市化作為最典型的土地利用變化過程之一，隨著其進程的不斷加快，城市周邊其它土地利用類型逐漸被開發(fā)為城鎮(zhèn)用地，導(dǎo)致城市地區(qū)不透水層面積增加，下墊面的滲透性、滯水性和熱力狀況隨之發(fā)生變化，使得截留、填洼、下滲和蒸發(fā)量減少，從而導(dǎo)致地表徑流和洪峰流量增加，徑流系數(shù)增大，顯著影響區(qū)域水循環(huán)過程［35］。葛怡等［36］、鄭璟等［37］利用SCS模型研究發(fā)現(xiàn)，上海和深圳地區(qū)城市化的發(fā)展會增加徑流系數(shù)，徑流量趨于增大。UNESCO的研究報告認(rèn)為城市化過程引起的水文效應(yīng)中最突出的表現(xiàn)為下滲減少和洪峰流量增大［8］。

3.1.3 水土保持措施對徑流的影響 水土保持措施主要包括以植樹造林、種草為主的植被措施、坡面工程措施（包括梯田、反坡梯田、水平階及魚鱗坑）、耕作措施和溝道工程措施等［38］，這些措施通過改變流域下墊面條件，從而影響區(qū)域的水文過程。一般來說，植樹造林和種草等植被措施具有減少年徑流量、增加枯水徑流和削減洪水的作用［34］，黃河水利委員會綏德水土保持試驗站的實測資料表明，當(dāng)林地蓋度為30%，50%，70%時，地表徑流將分別減少53%，86%，94%；而在同樣的蓋度條件下，草地分別減少地表徑流45%，75%，89%［39］。此外，林草地的空間分布格局也會影響水文效應(yīng)，草地—林地—坡耕地和林地—草地—坡耕地（由坡頂至坡底）的分布結(jié)構(gòu)，由于在徑流的形成區(qū)進行了控制，可以相對有效地減少徑流量，增強水保效果［40］。各種坡面工程措施具有類似的水保作用，其中以梯田最為典型。梯田通過改變局地的下墊面條件，不僅具有減少入河徑流總量的作用，同時也可以減緩河道泥沙的淤積，在一定的降雨條件下減水減沙作用突出。山西水保所1957—1966年在王家溝流域的觀測結(jié)果表明，水平梯田10 a平均減水減蝕效益達到92.4%［41］。總的來說，實施植樹種草和梯田等水土保持措施后，下墊面條件發(fā)生了變化，增加入滲、蒸散（發(fā)）量，對流域年徑流量會有不同程度的減少［42］。但不同地理位置、不同氣象條件也可能使不同地區(qū)水土保持的水文作用有所差異［34］。

3.1.4 農(nóng)業(yè)開發(fā)活動對徑流的影響 隨著社會的不斷發(fā)展和人口的急劇增長，人類對于糧食的需求日益增加，導(dǎo)致農(nóng)業(yè)用地不斷被開墾，并且土地利用的開發(fā)強度持續(xù)增加，但是耕作措施會使土壤結(jié)皮和壓實，從而使得土壤入滲速率減小，土壤容重增加［43］。Costa等［25］、Githui等［15］的研究表明，農(nóng)業(yè)用地的增加，降低了入滲和蒸發(fā)，從而增加了年徑流量。而在我國東北黑土區(qū)，由自然土壤—植被系統(tǒng)轉(zhuǎn)化為人工農(nóng)田生態(tài)系統(tǒng)，黑土儲水量將減少17.6～30.8%［44］。總體來說，農(nóng)業(yè)開發(fā)活動具有增加徑流量和洪峰流量的作用。

3.2 研究方法

目前研究流域土地利用／覆被變化對水文過程影響的方法可歸結(jié)為三類：試驗流域法、水文特征參數(shù)時間序列法和流域水文模型模擬法。

土地利用／覆被變化對水文過程的影響研究，早期大都采用試驗流域法。該方法有利于揭示植被—土壤—大氣相互作用的機理，比較容易獲得試驗結(jié)果。但這一方法僅適用于較小的流域，也不可能找到地理和氣象條件完全一致的兩個小流域，研究結(jié)果難以外推到其它流域，可對比性較差［45］。而水文特征參數(shù)時間序列法通過選擇較長時間段上能夠反映L UCC水文響應(yīng)的特征參數(shù)，從特征參數(shù)的變化趨勢上評價LUCC的水文效應(yīng)［46］。由于表征水文響應(yīng)的特征參數(shù)計算簡單，因此，該方法成為研究LUCC對水文過程影響的一種簡單有效的方法，并得到了廣泛的應(yīng)用［27，47－49］，但此方法僅適于下墊面條件比較均勻的流域，并且水文特征參數(shù)的變化受氣候和LUCC等多種因素的綜合影響，統(tǒng)計時間尺度也會影響評價結(jié)果的準(zhǔn)確性。

1970年Onstad等［7］最早嘗試?yán)盟哪Ｐ脱芯客恋乩玫乃捻憫?yīng)，LUCC的水文效應(yīng)研究也由傳統(tǒng)的統(tǒng)計分析方法轉(zhuǎn)向具有一定物理基礎(chǔ)的水文模型方法。至今已有許多水文模型用于LUCC的水文效應(yīng)研究，包括集總式水文模型，如 HBV、SCS、HSPF、CHARM等，半分布式／分布式水文模型，如TOPMODEL、VIC、SHE、SWAT等。由于基于物理基礎(chǔ)的分布式模型能夠考慮降水等因子和下墊面的空間異質(zhì)性，并且可以靈活地設(shè)置土地利用變化情景，從而有效地模擬不同土地利用情景下的水文響應(yīng)［34，50］，因而成為土地利用／覆被變化水文效應(yīng)研究的重要工具［30，50］。

目前，分布式水文模型主要有兩種建模思路和計算方法。第一種是應(yīng)用數(shù)值分析來建立相鄰網(wǎng)格單元之間的時空關(guān)系，如各種版本的SHE模型、I HDM模型等。該類模型均以1969年Freeze和Harlan發(fā)表的“一個具有物理基礎(chǔ)數(shù)值模擬的水文響應(yīng)模型的藍圖”作為描述徑流過程的基本框架［51］，以質(zhì)量、能量和動量方程描述自然系統(tǒng)，并考慮各變量參數(shù)的空間差異性，參數(shù)原則上具有明確的物理意義，可通過實測資料估計，但其結(jié)構(gòu)復(fù)雜，計算繁瑣，當(dāng)前很難適用于較大的流域［52］。第二種建模思路是在每一個網(wǎng)格單元（或子流域）上應(yīng)用傳統(tǒng)的概念性模型來推求凈雨，再進行匯流演算，最后求得出口斷面流量，如SWAT模型等。該類模型結(jié)構(gòu)與計算過程比較簡單，適用于較大的流域。此外還有一類模型基于DEM推求地形空間變化信息，利用地形信息模擬水文響應(yīng)，最后用統(tǒng)計方法求算流域總出口斷面流量。如1979年Beven和Kir kby［53］提出的以變源產(chǎn)流為基礎(chǔ)的TOPMODEL模型。這類模型的優(yōu)勢在于可以反映流域下墊面中影響水文響應(yīng)的關(guān)鍵因子（地形等）的空間變異性，但卻不需要像全分布式模型那樣引入大量參數(shù)，減小了參數(shù)率定帶來的不確定性［54］。

關(guān)于模型數(shù)據(jù)和參數(shù)問題，由于水文循環(huán)過程的高度非線性和復(fù)雜性，導(dǎo)致分布式水文模型需要大量的觀測數(shù)據(jù)或水文過程參數(shù)，不僅包括氣象、DEM、土壤、植被等自然環(huán)境要素數(shù)據(jù)，而且需要復(fù)雜的人工取用水、水庫徑流調(diào)節(jié)、農(nóng)業(yè)灌溉等人類活動數(shù)據(jù)，從而使得某些區(qū)域／流域建模所需的與水文過程尺度相匹配的觀測數(shù)據(jù)很難滿足分布式水文模擬的需求。Seyfried和 Wilcox［55］就曾認(rèn)為，物理基礎(chǔ)模型的缺點之一是較集總式模型需要更多的輸入資料，從而導(dǎo)致參數(shù)化和檢驗的工作量大大加大。

理論上，分布式水文模型的參數(shù)可以通過實測或者其它方式獲取。但在實際建模過程中，大量參數(shù)并不能由實測物理量確定，這就需要利用實測流量等數(shù)據(jù)率定網(wǎng)格或者子流域尺度上的有效參數(shù)值，由之得出合理的流域尺度上的模擬結(jié)果［54］。但是參數(shù)率定過程復(fù)雜，并且水文模擬和預(yù)測的參數(shù)估計都是點估計，結(jié)合擬合程度檢驗?zāi)Ｐ?，并不能夠降低不確定性的區(qū)間預(yù)測［56］。因此近年來水文模擬中的不確定性問題越來越受到重視，主要研究方法有SUFI-2（sequential uncertainty fitting）、GLUE（Generalized Likelihood Uncertainty Esti mation）、Ba RE（Bayesian Recursive Esti mation）、ParaSol（parameter sol ution）等。在參數(shù)優(yōu)選方面，啟發(fā)式算法和進化算法，如基因算法（GA）、SCE-UA等有效的全局參數(shù)率定方法已經(jīng)取代了局部尋優(yōu)的方法，而今后的發(fā)展趨勢是采用多目標(biāo)全局優(yōu)化技術(shù)［56］，從而降低參數(shù)估計的不確定性。

由于LUCC水文效應(yīng)的各種研究方法均存在局限性，近年來出現(xiàn)了多種綜合方法，如水文模型與統(tǒng)計學(xué)方法相結(jié)合的方法、模型耦合法以及模型對比法等。

3.3 模擬尺度

Bloschl和Sivapalan［57］將水文尺度劃分為過程尺度、觀測尺度和模擬尺度。隨著水文模型在LUCC與水文過程的相互作用研究中的廣泛應(yīng)用，研究并確定合適的模擬尺度就顯得愈加重要。理論上，模擬尺度應(yīng)盡量與過程尺度吻合，但是受到模擬水平的限制，實際很難達到，需要通過尺度轉(zhuǎn)換彌補這種缺陷［58］。

分布式水文模擬方法是尺度轉(zhuǎn)換的途徑之一［59］。分布式水文模型通過將流域離散化，尋找一種臨界單元考慮其內(nèi)部的地形、土壤和降雨空間分布的一致性，在這些離散單元中建立陸面過程與水文響應(yīng)的關(guān)系。但是通常這些單元的分辨率大于輸入數(shù)據(jù)的空間分辨率，使得每個計算單元中的參數(shù)值會存在一定程度的集總，而參數(shù)的集總對模擬結(jié)果有一定影響［60］。

Fitz Hugh等［60］研究認(rèn)為，徑流模擬結(jié)果對子流域劃分?jǐn)?shù)目的敏感性較小。郝芳華等［61］、張雪松等［62］、Jha 等［63］、Muleta 等［64］得 到 了 類 似 的 結(jié) 論。但是上述研究大多是在年和月尺度得出的結(jié)果，而Ku mar等［65］認(rèn)為子流域劃分?jǐn)?shù)量會顯著影響日水文過程模擬結(jié)果。因此可以認(rèn)為，在年和月尺度上，不同流域離散化水平下的徑流模擬精度無明顯差別，然而日徑流模擬結(jié)果受到子流域劃分層次的影響。

盡管許多研究認(rèn)為，子流域劃分水平對徑流的影響較小，但同時也存在合理的子流域劃分水平，當(dāng)超出此水平時可能會使模擬精度降低［66］，或產(chǎn)生過多面積較小且狹長的子流域，導(dǎo)致模型對水文過程的不真實模擬［67］。Ro manowicz等［66］認(rèn)為SWAT模型中合理的最小集水面積為模型推薦值的60%，也有研究認(rèn)為應(yīng)該為流域面積的4%［68］。

在不同的研究尺度下，分布式水文模擬尺度均會對模擬結(jié)果產(chǎn)生不同程度的影響，從而導(dǎo)致土地利用／覆被變化對水文過程影響結(jié)果的差異。需要針對具體的區(qū)域特征確定合理的子流域劃分水平，減小分布式水文模擬的不確定性，從而準(zhǔn)確評價土地利用／覆被變化對水文過程的影響。

4 存在不足及研究展望

土地利用／覆被變化對水文過程的影響研究已取得較大進展，但仍存在許多不足有待解決：① 由于有地域特性差異，以及研究方法的不同而導(dǎo)致流域LUCC對水文過程的影響存在差異，研究結(jié)果不具有可比性。需要在相同的研究尺度和方法下開展大量典型區(qū)域的實例研究，正確評價土地利用／覆被變化對水文過程的影響；② 目前主要通過水文站實測流量來率定和驗證分布式水文模型的有效性，但參數(shù)等效性問題在分布式水文模型中仍然沒有得到很好地解決。如何結(jié)合遙感技術(shù)為分布式水文模型提供可信度較高的水文過程關(guān)鍵變量的空間反演結(jié)果還需進一步研究；③需根據(jù)研究區(qū)域的特點和研究目的，選取合適的研究方法，而多種方法的集成和耦合仍有待深入；④ 土地利用／覆被變化對水文過程的影響研究較多，但水文過程的生態(tài)效應(yīng)研究相對較少。

在今后土地利用／覆被變化對水文過程的影響研究中，有必要從空間格局上研究LUCC對流域水文過程的影響機制，而不是僅僅關(guān)注土地利用類型變化的影響結(jié)果。并且在全球洪澇和干旱等極端天氣頻發(fā)的背景下，有必要加強土地利用變化對極值水文事件的影響，以及生態(tài)環(huán)境效應(yīng)方面的研究；在研究方法上，基于物理機制的分布式水文模型與GIS和RS技術(shù)更加緊密地結(jié)合應(yīng)用，提高對遙感數(shù)據(jù)的利用以及與GIS功能的耦合程度將會是重要的研究方向。

［1］ 夏軍，談戈.全球變化與水文科學(xué)新的進展與挑戰(zhàn)［J］.資源科學(xué)，2002，24（3）：1-7.

［2］ Voros marty C J，Green P，Salisbur y J，et al.Global water resources：vulnerability from climate change and population gr owth［J］.Science，2000，289（5477）：284-288.

［3］ Lah mer W，Pfützner B，Becker A.Assess ment of land use and climate change i mpacts on the mesoscale［J］.Physics and Chemistry of the Earth，2001，26（7／8）：565-575.

［4］ 楊久春，張樹文.近50年來呼倫湖水系土地利用／覆被變化及其生態(tài)環(huán)境效應(yīng)［J］.干旱區(qū)資源與環(huán)境，2009，23（2）：41-46.

［5］ DeFries R，Eshleman K N.Land-use change and hydrologic processes：a major f ocus f or the f uture［J］.Hydrological Processes，2004，18（11）：2183-2186.

［6］ 夏軍，左其亭.國際水文科學(xué)研究的新進展［J］.地球科學(xué)進展，2006，21（3）：256-261.

［7］ Onstad C A，Jamieson D G.Modeling the effects of land use modifications on runoff［J］.Water Resources Research，1970，6（5）：1287-1295.

［8］ Unesco.Hydrological effects of urbanization［R］.Studies and Reports in Hydrology 18，Paris，F(xiàn)rance，1974.

［9］ 劉昌明，鐘駿襄.黃土高原森林對年徑流影響的初步分析［J］.地理學(xué)報，1978（2）：12-16.

［10］ Matheussen B，Kirschbaum R L，Good man I A，et al.Effects of land cover change on streamflow in the interior Colu mbia River Basin （USA and Canada）［J］.Hydrological Pr ocesses，2000，14（5）：867-885.

［11］ 史培軍，袁藝，陳晉.深圳市土地利用變化對流域徑流的影響［J］.生態(tài)學(xué)報，2001，21（7）：1041-1050.

［12］ 陳軍鋒，李秀彬.土地覆被變化的水文響應(yīng)模擬研究［J］.應(yīng)用生態(tài)學(xué)報，2004，15（5）：833-836.

［13］ 郝芳華，陳利群，劉昌明，等.土地利用變化對產(chǎn)流和產(chǎn)沙的影響分析［J］.水土保持學(xué)報，2004，18（3）：5-8.

［14］ Cao W Z，Bowden W B，Davie T，et al.Modeling i mpacts of land cover change on critical water resources in the Motueka river catch ment，New Zealand［J］.Water Resources Management，2009，23（1）：137-151.

［15］ Githui F，Mutua F，Bau wens W.Esti mating the i mpacts of land-cover change on r unoff using the soil and water assessment tool（SWAT）：case study of Nzoia catch ment，Kenya［J］.Hydrological Sciences Jour nal，2009，54（5）：899-908.

［16］ Im S，Kim H，Kim C，et al.Assessing the i mpacts of land use changes on watershed hydrology using MIKE SHE［J］.Envir on mental Geology，2009，57（1）：231-239.

［17］ Mao D，Keith A C.Impacts of land-use change on hydrologic responses in the Great Lakes region［J］.Journal of Hydrology，2009，374（1／2）：71-82.

［18］ 陳瑩，許有鵬，尹義星.基于土地利用／覆被情景分析的長期水文效應(yīng)研究：以西苕溪流域為例［J］.自然資源學(xué)報，2009，24（2）：351-359.

［19］ 李佳，張小詠，楊艷昭.基于SWAT模型的長江源土地利用／土地覆被情景變化對徑流影響研究［J］.水土保持研究，2012，19（3）：119-124.

［20］ Wijesekara G N，Gupta A，Valeo C，et al.Assessing the i mpact of f ut ure land-use changes on hydrological processes in the Elbow River watershed in southern Alberta，Canada［J］.Jour nal of Hydrology，2012，412／413：220-232.

［21］ Calder I R.Hydrologic Effects of Land-use Change［C］∥Maid ment D R.Handbook of hydrolog.Mc Graw-Hill，New Yor k，1993：50.

［22］ Bewket W，Ster k G.Dynamics in land cover and its effect on stream flow in the Chemoga watershed，Blue Nile basin，Ethiopia［J］.Hydrological Processes，2005，19（2）：445-458.

［23］ 黃秉維.確切地估計森林的作用［J］.地理知識，1981（1）：1-5.

［24］ 周曉峰，趙惠勛，孫慧珍.正確評價森林水文效應(yīng)［J］.自然資源學(xué)報，2001，16（5）：420-426.

［25］ Costa M H，Botta A，Car dille J A.Effects of lar gescale changes in land cover on the discharge of the Tocantins River，Southeastern Amazonia［J］.Journal of Hydrology，2003，283（1／4）：206-217.

［26］ Jackson R B，Jobbagy E G，Avissar R，et al.Trading water for car bon with biological car bon sequestration［J］.Science，2005，310（5756）：1944-1947.

［27］ 王盛萍，張志強，孫閣，等.黃土高原流域土地利用變化水文動態(tài)響應(yīng)：以甘肅天水呂二溝流域為例［J］.北京林業(yè)大學(xué)學(xué)報，2006，28（1）：48-54.

［28］ Siri war dena L，F(xiàn)inlayson B L，Mc mahon T A.The i mpact of land use change on catch ment hydrology in lar ge catch ments： the comet River， Central Queensland，Australia［J］.Jour nal of Hydrology，2006，326（1／4）：199-214.

［29］ Zhang Y K，Schilling K E.Increasing streamflow and baseflow in Mississippi River since the 1940s：effect of land use change［J］.Jour nal of Hydrology，2006，324（1／4）：412-422.

［30］ Wang S F，Kang S Z，Zhang L，et al.Modeling hydrological response to different land-use and climate change scenarios in the Za mu River basin of north west China［J］.Hydrological Pr ocesses，2008，22（14）：2502-2510.

［31］ 陳軍鋒，李秀彬.森林植被變化對流域水文影響的爭論［J］.自然資源學(xué)報，2001，16（5）：474-480.

［32］ 王根緒，張鈺，劉桂民，等.馬營河流域1967—2000年土地利用變化對河流徑流的影響［J］.中國科學(xué)D輯：地球科學(xué)，2005，35（7）：671-681.

［33］ Olang L O，F(xiàn)urst J.Effects of land cover change on flood peak discharges and r unoff volu mes：model estimates for the Nyando River Basin，Kenya［J］.Hydrological Processes，2011，25（1）：80-89.

［34］ 李麗娟，姜德娟，李九一，等.土地利用／覆被變化的水文效應(yīng)研究進展［J］.自然資源學(xué)報，2007，22（2）：211-224.

［35］ 冉茂玉.論城市化的水文效應(yīng)［J］.四川師范大學(xué)學(xué)報：自然科學(xué)版，2000，23（4）：436-439.

［36］ 葛怡，史培軍，周俊華，等.土地利用變化驅(qū)動下的上海市區(qū)水災(zāi)災(zāi)情模擬［J］.自然災(zāi)害學(xué)報，2003，12（3）：25-30.

［37］ 鄭璟，方偉華，史培軍，等.快速城市化地區(qū)土地利用變化對流域水文過程影響的模擬研究：以深圳市布吉河流域為例［J］.自然資源學(xué)報，2009，24（9）：1560-1572.

［38］ 黃奕龍，傅伯杰，陳利頂.黃土高原水土保持建設(shè)的環(huán)境效應(yīng)［J］.水土保持學(xué)報，2003，17（1）：29-32.

［39］ 景可，鄭粉莉.黃土高原水土保持對地表水資源的影響［J］.水土保持研究，2004，11（4）：11-12，73.

［40］ 傅伯杰，陳利頂，馬克明.黃土丘陵小流域土地利用變化對生態(tài)環(huán)境的影響［J］.地理學(xué)報，1999，54（3）：241-246.

［41］ 姚云峰，王禮先.水平梯田減蝕作用分析［J］.中國水土保持，1992（12）：40-41.

［42］ 何長高，董增川，石景元，等.水土保持的水文效應(yīng)分布式模擬［J］.水科學(xué)進展，2009，20（4）：584-589.

［43］ Ankeny M D，Kaspar T C，Horton R.Characterization of tillage and traffic effects on unconfined infiltration measurements［J］.Soil Science Society of America Journal，1990，54（3）：837-840.

［44］ 韓曉增，王守宇，宋春雨，等.土地利用／覆蓋變化對黑土生態(tài)環(huán)境的影響［J］.地理科學(xué)，2005，25（2）：203-208.

［45］ 姚允龍，呂憲國，王蕾.流域土地利用／覆被變化水文效應(yīng)研究的方法評述［J］.濕地科學(xué)，2009，7（1）：83-88.

［46］ 張蕾娜，李秀彬.用水文特征參數(shù)變化表征人類活動的水文效應(yīng)初探：以云州水庫流域為例［J］.資源科學(xué)，2004，26（2）：62-67.

［47］ 鄧慧平，李秀彬，張明，等.氣候與地表覆被變化對梭磨河流域水文影響的分析［J］.地理科學(xué)，2001，21（6）：493-497.

［48］ 張曉明，余新曉，武思宏，等.黃土丘陵溝壑區(qū)典型流域土地利用／土地覆被變化水文動態(tài)響應(yīng)［J］.生態(tài)學(xué)報，2007，27（2）：414-423.

［49］ 朱麗，秦富倉，姚云峰.華北土石山區(qū)土地利用／土地覆被變化的水文響應(yīng)研究：以張家口云州水庫流域為例［J］.水土保持研究，2009，16（6）：224-228.

［50］ 徐宗學(xué)，程磊.分布式水文模型研究與應(yīng)用進展［J］.水利學(xué)報，2010，41（9）：1009-1017.

［51］ Freeze R A，Harlan R L.Buleprint for a physically－based，digitally-si mulated hydrological response model［J］.Journal of Hydrology，1969，9（3）：237-258.

［52］ 王中根，劉昌明，左其亭，等.基于DEM的分布式水文模型構(gòu)建方法［J］.地理科學(xué)進展，2002，21（5）：430-439.

［53］ Beven K J，Kir kby M J.A physically based variable contributing model of basin hydrology［J］.Hydrological Sciences Bulletin，1979，24（1）：43-69.

［54］ 張金存，芮孝芳.分布式水文模型構(gòu)建理論與方法述評［J］.水科學(xué)進展，2007，18（2）：286-292.

［55］ Seyfried M S，Wilcox B P.Scale and the nature of spatial variability：field examples having i mplications for hydr ologic modeling［J］.Water Resources Research，1995，31（1）：173-184.

［56］ 沈冰，黃領(lǐng)梅，李懷恩.水文模擬研究評述［J］.西安理工大學(xué)學(xué)報，2004，20（4）：351-355.

［57］ Bloschl G，Sivapalan M.Scale issues in hydrological modeling：a review ［J］.Hydrological Processes，1995，9（3／4）：251-290.

［58］ 任立良，劉新仁，郝振純.水文尺度若干問題述評［J］.水科學(xué)進展，1996，7（增刊）：87-99.

［59］ 劉建梅，裴鐵璠.水文尺度轉(zhuǎn)換研究進展［J］.應(yīng)用生態(tài)學(xué)報，2003，14（12）：2305-2310.

［60］ Fitz Hugh T W，Mackay D S.Impacts of input parameter spatial aggregation on an agricultural nonpoint source pollution model［J］.Jour nal of Hydrology，2000，236（1／2）：35-53.

［61］ 郝芳華，張雪松，程紅光，等.分布式水文模型亞流域合理劃分水平芻議［J］.水土保持學(xué)報，2003，17（4）：75-78.

［62］ 張雪松，郝芳華，程紅光，等.亞流域劃分對分布式水文模型模擬結(jié)果的影響［J］.水利學(xué)報，2004，35（7）：119-123.

［63］ Jha M，Gassman P W，Secchi S，et al.Effects of watershed subdivision on SWAT flow，sedi ments and nutrient predictions［J］.Jour nal of American Water Resources Association，2004，40（3）：811-825.

［64］ Muleta M K，Nicklow J W，Bekele E G.Sensitivity of a distributed watershed simulation to spatial scale［J］.Jour nal of Hydrologic Engineering，2007，12（2）：163-172.

［65］ Ku mar S，Mer wade V.Impact of watershed subdivision and soil data resolution on SWAT model calibration and parameter uncertainty［J］.Jour nal of the A-merican Water Resources Association，2009，45（5）：1179-1196.

［66］ Ro mano wicz A A，Vanclooster M，Rounsevell M，et al.Sensitivity of the SWAT model to the soil and land use data para metrisation：a case study in the Thyle catch ment，Belgiu m［J］.Ecological Modeling，2005，187（1）：27-39.

［67］ 胡連伍，王學(xué)軍，羅定貴，等.不同子流域劃分對流域徑流、泥沙、營養(yǎng)物模擬的影響［J］.水科學(xué)進展，2007，18（2）：235-240.

［68］ Arabi M，ovindaraju R S，antush M M，et al.Role of watershed subdivision on modeling the effectiveness of best management practices with SWAT［J］.Jour nal of the American Water Resources Association，2006，42（2）：513-528.