• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    土地利用變化下的水文響應(yīng)及其可檢測(cè)性分析

    2018-08-15 03:33:24,,
    關(guān)鍵詞:涪江水文林地

    , ,

    (1.長(zhǎng)江水利委員會(huì)水文局 水文情報(bào)預(yù)測(cè)中心, 武漢 430010; 2.長(zhǎng)江水利委員會(huì)水文局 長(zhǎng)江口水文水資源勘測(cè)局,上海 200136; 3.河海大學(xué) 水文水資源學(xué)院,南京 210098)

    1 研究背景

    土地覆被/土地利用作為水循環(huán)過程中重要的下墊面條件之一,對(duì)于全球和區(qū)域的水資源保護(hù)與生態(tài)可持續(xù)發(fā)展十分關(guān)鍵。土地利用變化不僅改變了自然景觀面貌,而且影響景觀中的物質(zhì)循環(huán)和能量分配,在客觀記錄地球表面特征空間格局活動(dòng)的同時(shí),也深刻地影響著水文循環(huán)的各個(gè)環(huán)節(jié)[1]。

    近年來,國(guó)內(nèi)外學(xué)者在土地利用變化下的水文響應(yīng)方面開展了大量研究工作,并得出了一些結(jié)論:森林植被的蒸散發(fā)量遠(yuǎn)遠(yuǎn)大于低矮植被[2],甚至接近水面蒸發(fā)[3];森林植被減少或者森林砍伐都會(huì)導(dǎo)致蒸散發(fā)量的明顯減少[4-6]。然而,一些研究卻表明盡管流域內(nèi)土地利用發(fā)生較為劇烈的變化,但是在流域上并沒有明顯地檢測(cè)到這些土地利用變化引起的水文要素變化[7-10]。因此,有必要進(jìn)一步開展土地利用變化引起的水文要素變化的可檢測(cè)性問題研究。本文以SWAT(Soil and Water Assessment Tool)模型為研究工具,采用情景分析方法,結(jié)合土地利用空間轉(zhuǎn)移轉(zhuǎn)化圖,來討論土地利用變化引起的水文要素變化的可檢測(cè)性問題。

    2 研究區(qū)概況

    涪江發(fā)源于四川省松潘縣境內(nèi)岷山雪寶頂北麓,是嘉陵江右岸最大支流。本文研究區(qū)在涪江中上游地區(qū)(圖1),介于平武縣和綿陽市區(qū)之間,平武為上游控制斷面,涪江橋?yàn)榱饔虺隹诳刂茢嗝妫饔蛎娣e7 570.7 km2,由于流域出口在涪江橋水文站,故稱之為涪江橋流域。

    圖1 涪江橋流域地理位置Fig.1 Location of Fujiangqiao catchment

    由圖1知,流域緯度范圍在31°29′N—32°33′N之間,經(jīng)度范圍在103°42′E—105°00′E之間。流域平均海拔1 683 m,西北部為高山地區(qū),最大海拔達(dá)4 753 m,東南是四川盆地,最低海拔為452 m。

    涪江橋流域?qū)賮啛釒駶?rùn)性季風(fēng)氣候區(qū),多年平武縣平均氣溫為14.8 ℃,綿陽市區(qū)為 16.5 ℃,多年平均降雨為822 mm。

    3 數(shù)據(jù)與方法

    3.1 數(shù)據(jù)源

    本文利用的數(shù)據(jù)包括水文氣象、土地利用、土壤和數(shù)字高程數(shù)據(jù)等。

    氣象數(shù)據(jù)為“中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)”(http:∥cdc.cma.gov.cn )提供的中國(guó)地面氣候資料日值數(shù)據(jù)集,采用的數(shù)據(jù)集要素包括:日最高氣溫、日最低氣溫、平均風(fēng)速。涪江橋流域內(nèi)共有2個(gè)氣象站點(diǎn),分別為平武站和綿陽站(見圖1)。

    水文數(shù)據(jù)為日降雨和日流量數(shù)據(jù),來自于水文年鑒《長(zhǎng)江流域水文資料第9冊(cè)嘉陵江區(qū)》。收集了研究區(qū)及附近25個(gè)雨量站的日降雨數(shù)據(jù)和4個(gè)水文站(平武、甘溪、麥地灣和涪江橋,圖1)的日流量數(shù)據(jù)。數(shù)據(jù)時(shí)段為1980—1987年。

    本文采用的土地利用數(shù)據(jù)是由中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心( http:∥www.resdc.cn )提供的中國(guó)1∶100 000土地利用土地利用/覆被現(xiàn)狀圖,共3期,分別為1985年、1995年和2000年土地利用土地利用/覆被現(xiàn)狀圖,如圖2所示。

    圖2 涪江橋流域土地利用空間分布Fig.2 Land use patterns of Fujiangqiao catchment

    數(shù)字高程模型(DEM)資料來源于SRTM[11],全稱是Shuttle Radar Topography Mission。SRTM3數(shù)據(jù)覆蓋全球,空間分辨率為3″(約90 m),平面基準(zhǔn)為WGS84。

    土壤數(shù)據(jù)為中國(guó)科學(xué)院南京土壤研究所編制的《中國(guó)1∶100萬數(shù)字化土壤屬性圖》[12],數(shù)據(jù)格式為柵格格式。該數(shù)據(jù)包括砂粒含量、粉粒含量、黏粒含量、重度、有機(jī)質(zhì)含量以及孔隙度等土壤屬性。

    3.2 研究方法

    本文采用SWAT模型模擬流域的水文過程。SWAT模型由Jeff Arnold博士為美國(guó)農(nóng)業(yè)部農(nóng)業(yè)科學(xué)研究所(USDA-ARS)開發(fā)的流域尺度的水文模型。多年來,該模型在美國(guó)、歐洲、中國(guó)、印度、伊朗以及韓國(guó)等地區(qū)得到了廣泛應(yīng)用[13]。大量學(xué)者將該模型應(yīng)用于氣候變化和土地利用變化的水文響應(yīng)研究中[14-16],并取得了很好的效果。SWAT模型由3個(gè)部分組成:子流域水文循環(huán)過程、河道徑流演算、水庫水量平衡過程。模型結(jié)構(gòu)及基本原理可參考文獻(xiàn)[17]。

    分別采用1985年和1995年2個(gè)時(shí)期土地利用數(shù)據(jù),保持其他水文模型參數(shù)不變,進(jìn)行水文過程模擬。為了消除氣候變化對(duì)水文過程的影響,模擬期間均采用了1980—1987年的水文氣象驅(qū)動(dòng)數(shù)據(jù)。2種情景下水文要素的差值被認(rèn)為是由土地利用變化引起的。將土地利用變化引起的水文要素變化量與土地利用轉(zhuǎn)移轉(zhuǎn)化圖相結(jié)合,來討論土地利用引起的水文要素變化的可檢測(cè)性問題。

    4 結(jié)果與討論

    4.1 土地利用變化分析

    涪江橋流域土地利用面積為7 570.7 km2,以耕地、林地和草地為主,3類土地利用占流域總面積的99%以上。其中,林地面積最大,約占流域面積的55%~64%;其次為耕地,約占流域面積的20%~21%;草地面積約占流域面積的14%~24%。

    涪江流域土地利用面積變化情況見表1。1985—1995年間,林地大面積減少,減少面積為620.4 km2;草地面積大幅度增加了682.6 km2;耕地面積略有減少;其他土地利用面積有所增加,但數(shù)量不大。從變化幅度上看,草地、水域和城鎮(zhèn)面積增幅較大,其中草地面積增幅高達(dá)61.0%。水域和城鎮(zhèn)的面積變化不大,但變化幅度很大,主要由于其面積基數(shù)小。1995—2000年間,草地面積減少最多,達(dá)394.9 km2,減小幅度為21.9%;林地面積增加速度最快,增加了7.6%,增加面積為320.7 km2。1985—2000年間,耕地、水域面積基本維持不變,林地、草地轉(zhuǎn)移轉(zhuǎn)化劇烈。最終林地面積減少了299.7 km2,草地面積增加了287.7 km2。期間城鎮(zhèn)面積持續(xù)擴(kuò)張,并主要發(fā)生在1995—2000年間。

    表1 涪江橋流域土地利用面積變化情況Table 1 Change of land use area in Fujiangqiao catchment

    采用土地利用轉(zhuǎn)移矩陣來描述和衡量土地利用類型之間的轉(zhuǎn)換關(guān)系,其優(yōu)點(diǎn)在于能夠全面具體地刻畫區(qū)域土地利用變化的結(jié)構(gòu)特征以及各土地類型變化的方向。該方法來源于系統(tǒng)分析中對(duì)系統(tǒng)狀態(tài)與狀態(tài)轉(zhuǎn)移的定量描述,為國(guó)內(nèi)外學(xué)者所常用[18]。涪江橋流域土地利用轉(zhuǎn)移矩陣見表2—表4。

    表2 1985—1995年涪江橋流域土地利用轉(zhuǎn)移矩陣Table 2 Conversion matrix among land use types from 1985 to 1995 in Fujiangqiao catchment km2

    表3 1995—2000 年涪江橋流域土地利用轉(zhuǎn)移矩陣Table 3 Conversion matrix among land use types from 1995 to 2000 in Fujiangqiao catchment km2

    表4 1985—2000年涪江橋流域土地利用轉(zhuǎn)移矩陣Table 4 Conversion matrix among land use types from 1985 to 2000 in Fujiangqiao catchment km2

    由表2—表4可以看出,各年代土地利用類型間的主要轉(zhuǎn)換關(guān)系為耕地、林地、草地間的相互轉(zhuǎn)化(1985—1995年,累計(jì)轉(zhuǎn)入、轉(zhuǎn)出1 753.4 km2;1995—2000年,累計(jì)轉(zhuǎn)入、轉(zhuǎn)出1 853.4 km2; 1985—2000年, 累計(jì)轉(zhuǎn)入、轉(zhuǎn)出612.8 km2。

    以1985—1995年為例,期間有412.6 km2耕地轉(zhuǎn)化為其他土地利用類型,主要轉(zhuǎn)化為林地和草地,面積分別為244.4 km2和165.2 km2;有1 082.1 km2林地轉(zhuǎn)化為其他土地利用類型,主要轉(zhuǎn)化為耕地和草地,轉(zhuǎn)化面積為306.0 km2和776.1 km2;有258.7 km2草地轉(zhuǎn)化為其他土地利用類型,主要轉(zhuǎn)化為林地,轉(zhuǎn)化面積為217.3 km2。

    此外,1985—2000年間城鎮(zhèn)面積不斷增長(zhǎng),其增加面積主要由耕地改造而來。城鎮(zhèn)用地土地利用化程度最高,且最為穩(wěn)定,不易轉(zhuǎn)化為其他類型的土地利用。

    流域內(nèi)土地利用變化是一個(gè)動(dòng)態(tài)過程,在某種土地利用大面積轉(zhuǎn)化為其他土地利用類型的同時(shí)(轉(zhuǎn)出),也伴隨著其他土地利用類型轉(zhuǎn)化為該種土地利用(轉(zhuǎn)入)。這使得最終發(fā)生變化的土地利用數(shù)量遠(yuǎn)遠(yuǎn)小于實(shí)際發(fā)生變化的土地利用數(shù)量,如1985—1995年耕地類型為例,研究時(shí)段始末,耕地面積減少65.2 km2,占流域面積的4.1%,而耕地實(shí)際變化過程并非如此簡(jiǎn)單,期間有412.6 km2的耕地轉(zhuǎn)化為其他類型土地利用,有374.4 km2其他類型土地利用轉(zhuǎn)化為耕地,耕地相關(guān)的土地利用面積為787.0 km2(轉(zhuǎn)入轉(zhuǎn)出之和),變化面積占流域面積的10.4%。這種轉(zhuǎn)入轉(zhuǎn)出的動(dòng)態(tài)過程容易使人們忽略土地利用的真實(shí)變化,也勢(shì)必會(huì)對(duì)分析土地利用變化下的流域水文響應(yīng)產(chǎn)生影響。

    4.2 SWAT模型率定與驗(yàn)證

    采用相對(duì)誤差(RE)和納什效率系數(shù)(NSE)分別對(duì)模擬與實(shí)測(cè)過程的水量和擬合程度進(jìn)行控制。研究區(qū)降雨徑流系列數(shù)據(jù)為1980—1987年日值數(shù)據(jù)。以1980年作為模型預(yù)熱期,初始化模型各狀態(tài),1981—1985年作為模型率定期,1986—1987年作為模型驗(yàn)證期。采用人工試錯(cuò)法,分別率定SWAT模型參數(shù)。率定與驗(yàn)證結(jié)果見表5。由表5可知,相對(duì)誤差在±5%左右,月納什效率系數(shù)>0.9,日效率系數(shù)>0.85,該模型可用于土地利用變化下的水文響應(yīng)及其可檢測(cè)性研究。

    涪江橋流域各控制站點(diǎn)月流量過程如圖3所示,涪江橋流域各控制站點(diǎn)日流量過程如圖4所示。由圖3和圖4可見,實(shí)測(cè)流量和模擬流量很接近。

    表5 涪江橋流域各控制站點(diǎn)SWAT模型率定與驗(yàn)證結(jié)果Table 5 Calibration and validation results of SWAT model for control stations in Fujiangqiao catchment

    注:由于平武以上流域存在積雪融水,因此并未參與模擬計(jì)算,而是以該站實(shí)測(cè)流量數(shù)據(jù)作為直接徑流,接入模型

    圖3 涪江橋流域各控制站點(diǎn)月流量過程Fig.3 Monthly hydrographs at control stations of Fujiangqiao catchment for calibration and validation periods

    圖4 涪江橋流域各控制站點(diǎn)日流量過程Fig.4 Daily hydrographs at control stations of Fujiangqiao catchment for calibration and validation periods

    4.3 不同時(shí)期土地利用情景下的水文模擬

    涪江橋各子流域不同土地利用情景下水文要素統(tǒng)計(jì)見表6。根據(jù)表6,略加計(jì)算,1995年的水文要素相對(duì)于1985年的變化幅度均在1%以內(nèi)。

    在各產(chǎn)匯流分區(qū)和全流域水平,均不能明顯地檢測(cè)到兩土地利用情景下模擬所得水文要素的差別。通常,表6中所示的水文要素變化的量級(jí),我們會(huì)認(rèn)為其在模型和數(shù)據(jù)的不確定性之內(nèi),可以忽略不計(jì)。1985—1995年間,流域內(nèi)23.4%面積上的土地利用發(fā)生了變化,但這些變化只引起了1%以內(nèi)的水文要素變化。因此有必要進(jìn)一步探究產(chǎn)生這種現(xiàn)象的原因所在。

    表6 涪江橋各子流域不同土地利用情景下水文要素統(tǒng)計(jì)Table 6 Hydrological regimes in subcatchments under different land use scenarios

    4.4 水文響應(yīng)檢測(cè)性分析

    以蒸散發(fā)為例,將蒸散發(fā)變化量反映在土地利用轉(zhuǎn)移轉(zhuǎn)化圖中(圖5),以便更深入地了解土地利用變化與流域水文要素變化的關(guān)系。

    圖5 1985—1995年涪江橋流域各子流域土地利用和蒸散發(fā)量變化Fig.5 Changes of land use and evapotranspiration in subcatchments from 1985 to 1995

    由圖5可知,流域范圍內(nèi)蒸散發(fā)量的變化范圍在-15.4~17.6 mm之間,占流域平均蒸散發(fā)量的-3.9%~4.4%。蒸散發(fā)量的變化與土地利用變化空間分布一致。該圖基本反映出了不同土地利用的水文效應(yīng)差異,即林地蒸散發(fā)量的能力大于草地,草地蒸散發(fā)量的能力大于耕地。但圖5中同樣反映出一些問題,即劇烈的土地利用變化并沒有引起水文要素的響應(yīng),如子流域4和26等。

    表7列出了土地利用變化最劇烈的5個(gè)子流域的蒸散發(fā)量變化情況。

    表7 1985—1995年土地利用變化最劇烈的子流域Table 7 Subcatchments with the most dramatic land use change from 1985 to 1995

    注:正值為增加,負(fù)值為減少

    由表7可知,區(qū)內(nèi)主要土地利用變化為林地向草地的轉(zhuǎn)化,夾雜著一些林地向耕地、耕地向草地的變化??傮w上看,林地面積的減少(不管是轉(zhuǎn)化為耕地還是草地),使得流域的多年平均蒸散發(fā)量減少,以子流域11,13,15最為明顯。盡管子流域4中有大部分林地減少,但其蒸散發(fā)量的變化不明顯。這主要是由于該子流域內(nèi)伴隨著部分耕地面積轉(zhuǎn)化為草地。耕地蒸散發(fā)的能力在多年平均水平小于草地,耕地轉(zhuǎn)化為草地的變化使得子流域蒸散發(fā)增加,而林地轉(zhuǎn)化為草地使得子流域蒸散發(fā)量減小,2種效應(yīng)疊加在一起,子流域整體的蒸散發(fā)量變化就不明顯了。至于子流域26,大面積的林地轉(zhuǎn)為蒸散發(fā)能力較小的耕地和草地,卻只引起了微小的蒸散發(fā)量變化,可能是由于不同土地利用在豐枯年的水文效應(yīng)差異和土壤特性差異引起的[19]。

    表8給出了蒸散發(fā)量變化最劇烈的5個(gè)子流域的土地利用變化情況。

    表8 1985—1995年蒸散發(fā)量變化最劇烈的子流域

    由表7和表8可見,相同的子流域只有子流域11,說明水文要素變化最劇烈的子流域并非一定是土地利用變化最劇烈的子流域。值得注意的是子流域39,盡管林地面積減少,但蒸散發(fā)量增加,并且十分顯著,其原因可能是林地減少轉(zhuǎn)化為草地,耕地減少轉(zhuǎn)化為林地,林地和耕地之間蒸散發(fā)的差異要大于林地和草地之間的差異,從而使得總體蒸散發(fā)量增加。

    由上面分析可知?jiǎng)×业耐恋乩米兓]有引起顯著的蒸散發(fā)量變化。導(dǎo)致該現(xiàn)象的原因主要有以下3個(gè)方面:

    (1)不同土地利用的水文效應(yīng)相似性。當(dāng)子流域內(nèi)土地利用變化主要是兩種具有相似水文效應(yīng)的土地利用間的相互轉(zhuǎn)化時(shí),這些變化并沒有對(duì)水文過程產(chǎn)生太大影響。如草地和耕地的轉(zhuǎn)化。張濤[19]從多年平均尺度和水文響應(yīng)單元(HRU)尺度分析得出兩種土地利用水文效應(yīng)在某些情況下差別不大。

    (2)不同土地利用類型的水文效應(yīng)差異性。在同一子流域內(nèi),發(fā)生具有相反水文效應(yīng)的土地利用變化。如在子流域39,林地減少為草地,會(huì)使得子流域的蒸發(fā)量減少,而耕地向草地的轉(zhuǎn)變又會(huì)使得子流域的蒸散發(fā)量增加,兩種效應(yīng)相互中和,增加了水文要素變化的可檢測(cè)性難度。

    (3)統(tǒng)計(jì)中和效應(yīng)和均化效應(yīng)。統(tǒng)計(jì)中和效應(yīng)和均化效應(yīng)是統(tǒng)計(jì)水文要素變化時(shí)的尺度問題。通常,統(tǒng)計(jì)水文要素變化時(shí)是以全流域或者控制站以上控制面積為統(tǒng)計(jì)單位。這里我們將控制站的控制面積稱作產(chǎn)匯流分區(qū),如前文中所提及的甘溪、麥地灣和涪江橋,這些產(chǎn)匯流分區(qū)由多個(gè)子流域構(gòu)成。當(dāng)流域或產(chǎn)匯流分區(qū)內(nèi)統(tǒng)計(jì)的水文要素正向變化和負(fù)向變化同時(shí)存在時(shí),所得到的水文要素的統(tǒng)計(jì)結(jié)果將會(huì)因?yàn)榻y(tǒng)計(jì)中和效應(yīng)而減少。在進(jìn)行水文要素變化統(tǒng)計(jì)計(jì)算時(shí),當(dāng)產(chǎn)匯流分區(qū)內(nèi)存在多個(gè)未變化子流域或變化量很小的子流域時(shí),分區(qū)的水文要素統(tǒng)計(jì)結(jié)果會(huì)被平均化,從而影響水文要素變化的檢測(cè)。

    上述各種效應(yīng)組合到一起,最終導(dǎo)致土地利用變化引起的水文要素變化難以檢測(cè)。解決這一問題的可行途徑之一是根據(jù)土地利用動(dòng)態(tài)變化過程來研究其水文響應(yīng)。

    4.5 水文動(dòng)態(tài)響應(yīng)度

    為了明確水文響應(yīng)動(dòng)態(tài)變化,本節(jié)提出水文動(dòng)態(tài)響應(yīng)度指標(biāo)來評(píng)判土地利用引起的水文要素變化的相對(duì)劇烈程度。

    前面分析了影響水文響應(yīng)可檢測(cè)性的的3個(gè)因素。其中第1和第2個(gè)因素是土地利用的固有屬性,發(fā)生在子流域內(nèi)部,是客觀存在的,不能改變,只能通過不斷減小子流域面積(如SWAT 模型)或網(wǎng)格大小(如VIC 模型)來減弱這2種因素的影響。最理想的程度就是每個(gè)子流域面積或網(wǎng)格大小與土地利用數(shù)據(jù)分辨率的大小相同,這樣就不存在子流域或網(wǎng)格內(nèi)的相互中和。然而,由于水文模型自身限制,這種分辨率通常不易實(shí)現(xiàn)。統(tǒng)計(jì)均化效應(yīng)是統(tǒng)計(jì)單元內(nèi)水文響應(yīng)動(dòng)態(tài)變化的一部分,不能被消除。因此本文嘗試消除統(tǒng)計(jì)中和效應(yīng)對(duì)水文可檢測(cè)性的影響。

    本文定義水文動(dòng)態(tài)響應(yīng)度為

    (1)

    式中:D為水文動(dòng)態(tài)響應(yīng)度,表征水文響應(yīng)變化的相對(duì)劇烈程度,D=0時(shí),統(tǒng)計(jì)單元內(nèi)(相當(dāng)于前文提到的產(chǎn)匯流分區(qū))水文要素如降雨量、蒸散發(fā)量、徑流深等,無變化,D值越大,水文水文要素變化程度越劇烈;Δi為統(tǒng)計(jì)單元中第i個(gè)計(jì)算單元內(nèi)(相當(dāng)于前文提到的子流域)某水文要素變化量;Si為第i個(gè)計(jì)算單元的面積(km2);S為統(tǒng)計(jì)單元面積,即所有子流域面積之和(km2)。

    式(1)中將水文要素變化量通過絕對(duì)值處理后,消除了統(tǒng)計(jì)單元內(nèi)的中和效應(yīng),能更好地反映統(tǒng)計(jì)單元間水文響應(yīng)的劇烈程度。

    應(yīng)用式(1)計(jì)算所得甘溪、麥地灣以及涪江橋產(chǎn)匯流分區(qū)蒸散發(fā)的水文動(dòng)態(tài)響應(yīng)度分別為2.48,4.04,4.84。由此可知各產(chǎn)匯流分區(qū)中,涪江橋產(chǎn)匯流分區(qū)內(nèi)蒸散發(fā)變化劇烈程度最大,其次是麥地灣,甘溪最小。該指標(biāo)可以在一定程度上幫助研究者了解研究區(qū)特定水文要素變化的相對(duì)劇烈程度。

    5 結(jié) 論

    本文以SWAT模型為工具,采用情景分析方法,結(jié)合土地利用轉(zhuǎn)移轉(zhuǎn)化圖,分析土地利用變化引起的水文要素變化的可檢測(cè)性問題。得到以下結(jié)論:

    (1)通過涪江橋流域3個(gè)水文控制站的徑流數(shù)據(jù)對(duì)SWAT模型進(jìn)行率定與驗(yàn)證,結(jié)果表明該模型模擬結(jié)果較好,可用于土地利用變化下的水文響應(yīng)及其可檢測(cè)性研究。

    (2)2個(gè)時(shí)期土地利用情景下的水文模擬結(jié)果表明涪江橋流域發(fā)生較為劇烈的土地利用變化,但并未檢測(cè)到足夠顯著的水文要素變化。

    (3)通過將1985年和1995年2個(gè)時(shí)期土地利用情景的SWAT模型模擬結(jié)果與土地利用空間轉(zhuǎn)化圖進(jìn)行耦合研究,將影響可檢測(cè)性的原因歸納為:①不同土地利用類型水文效應(yīng)的相似性;②不同土地利用類型水文效應(yīng)差異性;③統(tǒng)計(jì)中和效應(yīng)和統(tǒng)計(jì)均化效應(yīng)。

    (4)為消除統(tǒng)計(jì)中和效應(yīng),提出了水文動(dòng)態(tài)響應(yīng)度指標(biāo)評(píng)價(jià)水文要素變化的相對(duì)劇烈程度,該指標(biāo)可以在一定程度上幫助研究者了解研究區(qū)特定水文要素變化的相對(duì)劇烈程度。

    猜你喜歡
    涪江水文林地
    2022年《中國(guó)水文年報(bào)》發(fā)布
    有驚無險(xiǎn)
    田間趣事
    水文
    水文水資源管理
    李培基《絢彩涪江岸》
    慶祝中華人民共和國(guó)成立70 周年禮贊母親河涪江綿陽作家“涪江溯源”文學(xué)采風(fēng)創(chuàng)作活動(dòng)成功舉辦
    丹東市林地分類研究
    水文
    淺談林地保護(hù)及恢復(fù)措施
    国产又黄又爽又无遮挡在线| 亚洲av中文av极速乱 | 午夜激情福利司机影院| 韩国av在线不卡| 18禁黄网站禁片免费观看直播| 国产精品av视频在线免费观看| 久久精品综合一区二区三区| 精品久久久久久久久av| 淫秽高清视频在线观看| 最近中文字幕高清免费大全6 | 成人二区视频| 天堂av国产一区二区熟女人妻| 日韩精品中文字幕看吧| 午夜福利高清视频| 久久精品国产自在天天线| 亚州av有码| 欧美日韩乱码在线| 久久天躁狠狠躁夜夜2o2o| 中文字幕精品亚洲无线码一区| 性欧美人与动物交配| 国产精品久久电影中文字幕| 国产精品伦人一区二区| 蜜桃久久精品国产亚洲av| 久久久久久久久久久丰满 | 亚洲美女黄片视频| 国产一级毛片七仙女欲春2| 日本五十路高清| 成人国产一区最新在线观看| 午夜福利在线观看免费完整高清在 | 人人妻人人澡欧美一区二区| 亚洲成人免费电影在线观看| 日韩国内少妇激情av| 国产伦精品一区二区三区四那| 亚洲天堂国产精品一区在线| 久久久久精品国产欧美久久久| 高清在线国产一区| 在线播放国产精品三级| 波多野结衣巨乳人妻| 国产午夜精品论理片| 欧美区成人在线视频| 国产高清不卡午夜福利| 亚洲中文字幕一区二区三区有码在线看| 成人av一区二区三区在线看| 国产三级中文精品| 少妇人妻一区二区三区视频| 国产主播在线观看一区二区| 欧美黑人欧美精品刺激| or卡值多少钱| 日本 欧美在线| 一本久久中文字幕| 精品无人区乱码1区二区| 国产精品久久电影中文字幕| 人妻丰满熟妇av一区二区三区| 尾随美女入室| 亚洲精品一卡2卡三卡4卡5卡| 国产真实伦视频高清在线观看 | 成人国产综合亚洲| 高清日韩中文字幕在线| 午夜精品久久久久久毛片777| 人妻丰满熟妇av一区二区三区| 不卡一级毛片| 日韩强制内射视频| 人妻久久中文字幕网| 天天躁日日操中文字幕| 中文字幕精品亚洲无线码一区| 午夜日韩欧美国产| 亚洲成人精品中文字幕电影| 免费观看人在逋| 精品无人区乱码1区二区| netflix在线观看网站| 久久欧美精品欧美久久欧美| 日本免费a在线| 国产单亲对白刺激| 日本成人三级电影网站| 日本色播在线视频| 俄罗斯特黄特色一大片| bbb黄色大片| 91麻豆av在线| 波野结衣二区三区在线| 有码 亚洲区| 在现免费观看毛片| 亚洲欧美日韩高清在线视频| 国产精品综合久久久久久久免费| 成年女人永久免费观看视频| 日韩 亚洲 欧美在线| 麻豆一二三区av精品| 嫩草影视91久久| 黄色视频,在线免费观看| 色综合色国产| 久久久久久九九精品二区国产| 亚洲精品粉嫩美女一区| 美女高潮的动态| 欧美+日韩+精品| 亚洲成人中文字幕在线播放| 午夜免费激情av| 精品一区二区三区av网在线观看| 国产女主播在线喷水免费视频网站 | 国产91精品成人一区二区三区| 自拍偷自拍亚洲精品老妇| 欧美又色又爽又黄视频| 3wmmmm亚洲av在线观看| 精华霜和精华液先用哪个| 午夜福利欧美成人| 国产精品国产三级国产av玫瑰| videossex国产| 如何舔出高潮| 欧美高清成人免费视频www| av.在线天堂| 欧美日韩精品成人综合77777| 国产精品福利在线免费观看| 成年女人看的毛片在线观看| 三级国产精品欧美在线观看| 国产精品乱码一区二三区的特点| 亚洲专区国产一区二区| 男女啪啪激烈高潮av片| 成人国产麻豆网| a级毛片免费高清观看在线播放| 观看美女的网站| 亚洲精品久久国产高清桃花| 久久6这里有精品| 高清毛片免费观看视频网站| 熟女人妻精品中文字幕| 中文字幕av在线有码专区| 精品人妻1区二区| 亚洲七黄色美女视频| 亚洲图色成人| 国产乱人视频| 欧美成人性av电影在线观看| 老女人水多毛片| 两个人的视频大全免费| 嫩草影院新地址| 直男gayav资源| 三级国产精品欧美在线观看| 乱人视频在线观看| 蜜桃久久精品国产亚洲av| 亚洲国产高清在线一区二区三| 亚洲内射少妇av| 国产亚洲欧美98| 免费无遮挡裸体视频| 在线免费观看的www视频| 人妻少妇偷人精品九色| 国产白丝娇喘喷水9色精品| 国产精品人妻久久久影院| 午夜激情福利司机影院| 欧美日韩综合久久久久久 | 午夜老司机福利剧场| 国产色婷婷99| 国产精品福利在线免费观看| 99久久精品国产国产毛片| 天堂av国产一区二区熟女人妻| 精品人妻视频免费看| 男女视频在线观看网站免费| 国产色婷婷99| 国产精品福利在线免费观看| 亚洲av免费高清在线观看| 午夜免费男女啪啪视频观看 | 女人十人毛片免费观看3o分钟| 天堂影院成人在线观看| 久久人人精品亚洲av| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av| 色在线成人网| 日本一二三区视频观看| 久久精品国产亚洲av涩爱 | 国产一区二区三区av在线 | 一区二区三区四区激情视频 | 免费黄网站久久成人精品| 亚洲中文字幕一区二区三区有码在线看| 久久久精品欧美日韩精品| 亚洲人成网站在线播| 麻豆久久精品国产亚洲av| 美女xxoo啪啪120秒动态图| 亚洲久久久久久中文字幕| 人妻夜夜爽99麻豆av| 美女xxoo啪啪120秒动态图| 别揉我奶头 嗯啊视频| 一进一出抽搐gif免费好疼| 欧美性感艳星| 亚洲av成人精品一区久久| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 色哟哟·www| 国产美女午夜福利| 麻豆成人av在线观看| 久久久久久久午夜电影| 国产大屁股一区二区在线视频| 日本免费一区二区三区高清不卡| 亚洲五月天丁香| 国产av麻豆久久久久久久| av国产免费在线观看| 国产精华一区二区三区| 欧美激情久久久久久爽电影| 久久久久久大精品| 九色国产91popny在线| 国产精品久久久久久亚洲av鲁大| 日本一二三区视频观看| av在线亚洲专区| 久久久久久伊人网av| 美女 人体艺术 gogo| 日韩精品中文字幕看吧| 黄色视频,在线免费观看| 级片在线观看| 黄片wwwwww| 国产精品嫩草影院av在线观看 | 国产一区二区在线av高清观看| 国产v大片淫在线免费观看| 日本黄大片高清| 天美传媒精品一区二区| 日本一本二区三区精品| av国产免费在线观看| 精品一区二区三区视频在线| 日本成人三级电影网站| 日日干狠狠操夜夜爽| 亚洲三级黄色毛片| 久久久久久伊人网av| 国产亚洲欧美98| 成年女人看的毛片在线观看| 欧美成人a在线观看| 美女被艹到高潮喷水动态| 最新中文字幕久久久久| 麻豆av噜噜一区二区三区| 欧美日韩综合久久久久久 | 天堂av国产一区二区熟女人妻| 日韩精品有码人妻一区| 美女高潮的动态| 国产精品一区二区性色av| 狠狠狠狠99中文字幕| 免费人成在线观看视频色| 男人舔女人下体高潮全视频| 一进一出抽搐动态| 亚洲电影在线观看av| 日本a在线网址| 免费在线观看日本一区| 他把我摸到了高潮在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久久午夜欧美精品| 两性午夜刺激爽爽歪歪视频在线观看| 三级国产精品欧美在线观看| 国产伦人伦偷精品视频| 亚洲精华国产精华液的使用体验 | 日本色播在线视频| 人妻少妇偷人精品九色| 欧美日韩黄片免| 两个人视频免费观看高清| av专区在线播放| 国产精品av视频在线免费观看| 在线看三级毛片| 欧美日韩综合久久久久久 | 真人做人爱边吃奶动态| 亚洲av美国av| 黄色丝袜av网址大全| 久久人人精品亚洲av| 99在线视频只有这里精品首页| 成年版毛片免费区| 国产精品自产拍在线观看55亚洲| 免费人成在线观看视频色| 九九久久精品国产亚洲av麻豆| 99热6这里只有精品| 亚洲中文字幕日韩| 亚洲精品影视一区二区三区av| 欧美日本视频| 精品午夜福利在线看| 欧美高清性xxxxhd video| 午夜福利在线在线| 亚洲人与动物交配视频| 人妻夜夜爽99麻豆av| 亚洲五月天丁香| 国产精品亚洲美女久久久| 久久精品影院6| 九九爱精品视频在线观看| 国产精品福利在线免费观看| 韩国av在线不卡| 身体一侧抽搐| 国产精华一区二区三区| 国产视频一区二区在线看| 欧美成人性av电影在线观看| 18+在线观看网站| 丰满人妻一区二区三区视频av| 极品教师在线视频| 精品久久久久久,| bbb黄色大片| 日本黄色视频三级网站网址| 久久久久国内视频| 欧美精品啪啪一区二区三区| 一夜夜www| 色综合婷婷激情| 免费高清视频大片| 一本精品99久久精品77| 男女做爰动态图高潮gif福利片| 免费看av在线观看网站| 最好的美女福利视频网| 免费在线观看成人毛片| 亚洲专区中文字幕在线| av女优亚洲男人天堂| 久久精品国产亚洲av涩爱 | 一区二区三区激情视频| 91午夜精品亚洲一区二区三区 | 尤物成人国产欧美一区二区三区| .国产精品久久| 久久久久久久久久黄片| 欧美xxxx性猛交bbbb| 91狼人影院| 搡女人真爽免费视频火全软件 | 日本精品一区二区三区蜜桃| 真人一进一出gif抽搐免费| 国产在线精品亚洲第一网站| 精品一区二区三区人妻视频| 久久久久久大精品| 国产亚洲欧美98| 亚洲色图av天堂| 久久这里只有精品中国| 欧美色视频一区免费| 搡老熟女国产l中国老女人| 久久99热这里只有精品18| 在线免费十八禁| 免费av观看视频| 最新中文字幕久久久久| 国产精品嫩草影院av在线观看 | 美女免费视频网站| 一边摸一边抽搐一进一小说| 成人欧美大片| 18禁黄网站禁片免费观看直播| 美女cb高潮喷水在线观看| 亚洲国产日韩欧美精品在线观看| 日韩欧美一区二区三区在线观看| 一个人看的www免费观看视频| 久久亚洲精品不卡| 午夜福利成人在线免费观看| 18禁黄网站禁片免费观看直播| 最新在线观看一区二区三区| 免费看a级黄色片| 午夜激情福利司机影院| 午夜免费成人在线视频| 最新在线观看一区二区三区| 深夜a级毛片| aaaaa片日本免费| 97碰自拍视频| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 熟女人妻精品中文字幕| 男人舔奶头视频| 国产真实乱freesex| 欧美三级亚洲精品| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 桃红色精品国产亚洲av| 亚洲电影在线观看av| a在线观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 少妇猛男粗大的猛烈进出视频 | 在线播放无遮挡| 91久久精品国产一区二区三区| a级一级毛片免费在线观看| 国产成年人精品一区二区| 最近最新中文字幕大全电影3| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 99久久精品热视频| 又爽又黄a免费视频| 日韩中字成人| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 欧美+亚洲+日韩+国产| 两个人的视频大全免费| 亚洲内射少妇av| 亚洲一区高清亚洲精品| 两个人的视频大全免费| 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 国产精品福利在线免费观看| 日韩亚洲欧美综合| 伊人久久精品亚洲午夜| 波多野结衣高清作品| 精品无人区乱码1区二区| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影| av黄色大香蕉| 最好的美女福利视频网| 亚洲黑人精品在线| 精品国内亚洲2022精品成人| 精品人妻视频免费看| 国产淫片久久久久久久久| 国产亚洲av嫩草精品影院| 少妇高潮的动态图| 亚洲精品久久国产高清桃花| 不卡一级毛片| 无人区码免费观看不卡| 国产午夜精品久久久久久一区二区三区 | 国产一区二区在线av高清观看| 日日摸夜夜添夜夜添小说| 欧美一区二区国产精品久久精品| 欧美日本视频| 女人十人毛片免费观看3o分钟| 少妇高潮的动态图| 午夜免费激情av| 永久网站在线| 午夜激情欧美在线| 亚洲精品日韩av片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av免费在线观看| 国产在线精品亚洲第一网站| 成年女人永久免费观看视频| 欧美日韩国产亚洲二区| 丝袜美腿在线中文| 全区人妻精品视频| 国产人妻一区二区三区在| 亚洲 国产 在线| 两人在一起打扑克的视频| 亚洲国产精品sss在线观看| 国产v大片淫在线免费观看| 男女下面进入的视频免费午夜| 直男gayav资源| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 极品教师在线免费播放| 黄色配什么色好看| 嫁个100分男人电影在线观看| 日本黄大片高清| bbb黄色大片| 午夜影院日韩av| 老司机福利观看| 精品一区二区免费观看| 国产探花在线观看一区二区| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添av毛片 | 亚洲中文字幕日韩| 亚洲人与动物交配视频| 97人妻精品一区二区三区麻豆| 精品久久久久久成人av| 国产精品人妻久久久影院| 国产伦精品一区二区三区四那| 乱人视频在线观看| 国产乱人伦免费视频| 99riav亚洲国产免费| 免费在线观看影片大全网站| aaaaa片日本免费| 亚洲精品成人久久久久久| 内地一区二区视频在线| 国产男人的电影天堂91| 亚洲中文日韩欧美视频| 男女之事视频高清在线观看| 欧美日韩综合久久久久久 | 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 久久精品国产亚洲av天美| 亚洲自偷自拍三级| 免费av观看视频| 一区二区三区四区激情视频 | 午夜久久久久精精品| 亚洲综合色惰| av专区在线播放| 偷拍熟女少妇极品色| 伦理电影大哥的女人| 嫁个100分男人电影在线观看| 国产精品综合久久久久久久免费| 欧美中文日本在线观看视频| 日日撸夜夜添| 看黄色毛片网站| 日日干狠狠操夜夜爽| 国产熟女欧美一区二区| 午夜激情福利司机影院| 最新中文字幕久久久久| 大又大粗又爽又黄少妇毛片口| 亚洲最大成人av| 国内精品一区二区在线观看| 午夜福利视频1000在线观看| 午夜福利在线在线| 内地一区二区视频在线| 哪里可以看免费的av片| 国产午夜精品论理片| 亚洲一区高清亚洲精品| 国产伦精品一区二区三区视频9| 亚洲国产精品久久男人天堂| 变态另类丝袜制服| 亚洲,欧美,日韩| 欧美成人一区二区免费高清观看| 在线观看av片永久免费下载| 男人狂女人下面高潮的视频| 在线免费观看的www视频| 看黄色毛片网站| 久久久久久久精品吃奶| 国产精品一区二区性色av| 桃色一区二区三区在线观看| 久久久成人免费电影| 窝窝影院91人妻| 日本爱情动作片www.在线观看 | 亚洲18禁久久av| 久久久久免费精品人妻一区二区| 一边摸一边抽搐一进一小说| 男人舔奶头视频| 蜜桃久久精品国产亚洲av| 亚洲av成人av| 亚洲最大成人中文| 亚洲欧美日韩高清专用| 少妇人妻一区二区三区视频| 男女那种视频在线观看| 偷拍熟女少妇极品色| 麻豆成人av在线观看| 51国产日韩欧美| 午夜福利欧美成人| 免费人成在线观看视频色| 天美传媒精品一区二区| 999久久久精品免费观看国产| 欧美色欧美亚洲另类二区| 国产高清激情床上av| 免费人成在线观看视频色| 嫩草影院新地址| 久久精品91蜜桃| 亚洲欧美激情综合另类| 日日撸夜夜添| 中文资源天堂在线| 嫁个100分男人电影在线观看| 久久欧美精品欧美久久欧美| 在线播放国产精品三级| 亚洲最大成人av| 三级男女做爰猛烈吃奶摸视频| bbb黄色大片| 国产一区二区三区av在线 | 亚洲国产欧美人成| 啦啦啦韩国在线观看视频| 长腿黑丝高跟| 久久国产精品人妻蜜桃| 免费电影在线观看免费观看| 亚洲人成网站高清观看| 欧美日韩乱码在线| 身体一侧抽搐| 色av中文字幕| 亚洲熟妇熟女久久| 人妻制服诱惑在线中文字幕| 久9热在线精品视频| 久久国产乱子免费精品| 亚洲三级黄色毛片| 69av精品久久久久久| ponron亚洲| 九九爱精品视频在线观看| 久久久久久久午夜电影| 在线观看66精品国产| 一级黄片播放器| 国产精品久久久久久精品电影| 日韩欧美在线乱码| a在线观看视频网站| 午夜影院日韩av| 国产在线精品亚洲第一网站| 亚洲av第一区精品v没综合| av在线蜜桃| 特大巨黑吊av在线直播| АⅤ资源中文在线天堂| 久久久久国产精品人妻aⅴ院| 成人性生交大片免费视频hd| 久久这里只有精品中国| 亚洲精华国产精华液的使用体验 | 两个人视频免费观看高清| 悠悠久久av| 精品久久久久久成人av| 亚洲精品国产成人久久av| 在线a可以看的网站| 桃色一区二区三区在线观看| 五月玫瑰六月丁香| 黄色配什么色好看| 51国产日韩欧美| 久久久久久久久大av| 国产伦在线观看视频一区| 成年人黄色毛片网站| 亚洲国产精品sss在线观看| 亚洲三级黄色毛片| 久久亚洲真实| 国产av在哪里看| 精品久久久久久久久久免费视频| 成人午夜高清在线视频| 又粗又爽又猛毛片免费看| 尤物成人国产欧美一区二区三区| 深爱激情五月婷婷| 1000部很黄的大片| 国产久久久一区二区三区| 欧美性猛交黑人性爽| 在线观看免费视频日本深夜| 美女高潮的动态| 国产av不卡久久| 免费av不卡在线播放| 最新在线观看一区二区三区| АⅤ资源中文在线天堂| .国产精品久久| 国产69精品久久久久777片| 免费高清视频大片| 婷婷精品国产亚洲av| 日本黄色片子视频| 99在线人妻在线中文字幕| 白带黄色成豆腐渣| 国内精品久久久久久久电影| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 久99久视频精品免费| 91麻豆精品激情在线观看国产| 内射极品少妇av片p| 可以在线观看的亚洲视频| 一进一出好大好爽视频| 日本熟妇午夜| 乱系列少妇在线播放| 久久久久久伊人网av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成人久久爱视频| 免费电影在线观看免费观看| 又黄又爽又免费观看的视频| 久久精品国产亚洲网站| 国产男人的电影天堂91| 亚洲专区国产一区二区| 国产精品伦人一区二区| 成人特级av手机在线观看| 在线天堂最新版资源| 成熟少妇高潮喷水视频| 乱人视频在线观看| 99久久久亚洲精品蜜臀av| 黄色丝袜av网址大全| 日本a在线网址|