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

    中國西北干旱區(qū)蒸散發(fā)時(shí)空動態(tài)特征

    2017-06-22 14:04:27鄧興耀劉志輝姚俊強(qiáng)
    生態(tài)學(xué)報(bào) 2017年9期
    關(guān)鍵詞:干旱區(qū)西北趨勢

    鄧興耀,劉 洋,劉志輝,姚俊強(qiáng)

    1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046 2 新疆大學(xué)干旱生態(tài)環(huán)境研究所,烏魯木齊 830046 3 新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046 4 中國氣象局烏魯木齊沙漠氣象研究所,烏魯木齊 830002

    ?

    中國西北干旱區(qū)蒸散發(fā)時(shí)空動態(tài)特征

    鄧興耀1,2,3,劉 洋1,2,3,劉志輝2,3,*,姚俊強(qiáng)4

    1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046 2 新疆大學(xué)干旱生態(tài)環(huán)境研究所,烏魯木齊 830046 3 新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046 4 中國氣象局烏魯木齊沙漠氣象研究所,烏魯木齊 830002

    利用MODIS ET數(shù)據(jù)集中2000—2014年的地表實(shí)際蒸散發(fā)量產(chǎn)品,運(yùn)用變異系數(shù)、Theil-Sen median趨勢分析與Mann-Kendall檢驗(yàn)和Hurst指數(shù)法,研究了中國西北干旱區(qū)蒸散發(fā)的空間格局、不同維度的空間異質(zhì)性和時(shí)間變化特征及未來趨勢預(yù)測。結(jié)果表明:(1)2000—2014年全區(qū)蒸散發(fā)量總體較小,蒸散發(fā)量小于200 mm的區(qū)域占總面積的38.329%。在空間上ET自山區(qū)向兩側(cè)平原減少,不同土地覆蓋的ET大小為:林地>農(nóng)用地>草地>稀疏植被。受降水和土地覆蓋的綜合影響,ET的高值區(qū)(>400 mm)主要在降水豐富的山區(qū)林地和草地,而低值區(qū)(<200 mm)主要在降水較少的平原稀疏植被區(qū)和草地。(2)近15年全區(qū)蒸散發(fā)變異程度不明顯,以相對較低的波動變化為主。各亞區(qū)內(nèi)波動較低區(qū)域的比例為:北疆>天山>祁連山>內(nèi)蒙西部>河西走廊>南疆。(3)15年間全區(qū)年均蒸散發(fā)量呈波動變化,總體有微弱的減小趨勢,變化率為-0.9348 mm/a?;谙裨叨鹊姆治鲆脖砻魅珔^(qū)ET以減小的變化趨勢為主,但各亞區(qū)的減小程度各異:天山>內(nèi)蒙西部>河西走廊>北疆>祁連山,僅南疆有增加趨勢。(4)全區(qū)ET的Hurst指數(shù)均值為0.689,Hurst指數(shù)大于0.5的范圍所占比例為80.033%,未來全區(qū)蒸散發(fā)的變化趨勢以持續(xù)性減小為主。其中22.003%區(qū)域的變化趨勢無法確定。未來各亞區(qū)ET的減少趨勢為:內(nèi)蒙西部>天山>河西走廊>北疆>祁連山>南疆。

    MODIS;蒸散發(fā);西北干旱區(qū);趨勢分析;氣候變化

    蒸散發(fā)(Evapotranspiration, ET)包括土壤、水面蒸發(fā)和植被蒸騰,連接著地表水分、能量和碳的循環(huán)過程,是構(gòu)成氣候系統(tǒng)的中心環(huán)節(jié)[1]。研究陸面過程中的蒸散發(fā),對天氣預(yù)報(bào)、旱澇監(jiān)測、水資源和農(nóng)業(yè)管理以及全球變化等領(lǐng)域有重要意義[2]。

    基于地面觀測的資料可獲得長時(shí)間序列的蒸散發(fā)信息,但站點(diǎn)觀測值并不能提供蒸散發(fā)的空間分布特征,尤其是在觀測站點(diǎn)稀疏的西北干旱區(qū)。結(jié)合遙感技術(shù)可以反映蒸散發(fā)的空間異質(zhì)性,滿足全球和區(qū)域尺度的研究。Jung等[3]發(fā)布了1982—2008年全球陸地蒸散發(fā)年尺度數(shù)據(jù)集;歐洲氣象衛(wèi)星應(yīng)用組織發(fā)布了覆蓋歐洲、非洲和南美洲東部的LSA-SAF MSG ET數(shù)據(jù)集[4];美國蒙大拿大學(xué)森林學(xué)院工作組制作了2000年至今的全球MODIS ET數(shù)據(jù)集[5]。

    MODIS ET數(shù)據(jù)集憑借較高模擬精度和時(shí)空分辨率,已成功應(yīng)用于全球和區(qū)域蒸散發(fā)的動態(tài)監(jiān)測。賀添等[6]利用中國陸地生態(tài)系統(tǒng)通量觀測研究網(wǎng)絡(luò)數(shù)據(jù)和水文數(shù)據(jù),驗(yàn)證了MODIS ET數(shù)據(jù)的精度。Kim等[7]利用亞洲通量網(wǎng)17個(gè)站點(diǎn)的渦度相關(guān)儀數(shù)據(jù)對MODIS ET產(chǎn)品進(jìn)行驗(yàn)證,并分析了該數(shù)據(jù)集在不同氣候和下墊面條件下的適用性;Liu等[8]利用中國海河平原的大孔徑閃爍儀觀測數(shù)據(jù),驗(yàn)證2008—2010年MODIS ET產(chǎn)品的精度,并分析了3種下墊面蒸散發(fā)量的季節(jié)和年變化規(guī)律;Jang等[9]利用MODIS ET產(chǎn)品和MODIS-KLDAS ET驅(qū)動數(shù)據(jù)監(jiān)測2006—2008年東北亞地區(qū)日蒸散發(fā)量,MODIS ET產(chǎn)品的精度通過了區(qū)域通量塔觀測的驗(yàn)證。

    中國西北干旱區(qū)是亞洲中部干旱區(qū)的重要組成部分之一,區(qū)內(nèi)地形復(fù)雜,山地盆地相間分布,沙漠與綠洲共存,各自然要素時(shí)空分布極不均勻,是生態(tài)環(huán)境脆弱地區(qū),也是對全球變化響應(yīng)最敏感地區(qū)之一[10]。全球氣候變化改變了陸面水循環(huán)要素和蒸散發(fā)量[3]。水資源是干旱區(qū)農(nóng)業(yè)和社會經(jīng)濟(jì)發(fā)展的制約因素。在此背景下,全球變化和人類活動共同驅(qū)動下的干旱區(qū)蒸散發(fā)時(shí)空分布及其變化受到諸多學(xué)者的關(guān)注。李寶富[11]、王海波[12]分別以西北干旱區(qū)典型內(nèi)陸河塔里木河流域、黑河流域?yàn)檠芯繀^(qū),發(fā)展了干旱區(qū)遙感估算蒸散發(fā)量的模型。馬金龍[13]、張鑫[14]利用渦度相關(guān)系統(tǒng)和波文比系統(tǒng)監(jiān)測綠洲農(nóng)田蒸散發(fā),為干旱區(qū)作物生長階段不同灌溉時(shí)期和灌溉量的確定以及田間水分管理提供科學(xué)依據(jù)。劉春雨[15]研究了氣候變化背景下的西北干旱區(qū)蒸散發(fā)響應(yīng)。李稚等[16]利用1958—2010年的氣象站觀測數(shù)據(jù),得出西北干旱區(qū)蒸發(fā)皿蒸發(fā)量在過去50年內(nèi)以1993年為轉(zhuǎn)折點(diǎn),由下降趨勢逆轉(zhuǎn)為顯著上升的趨勢的結(jié)論。但是由于不同的干濕背景,干旱區(qū)蒸發(fā)皿蒸發(fā)和實(shí)際蒸散發(fā)之間的關(guān)系尚存在較大爭論[17- 18]。

    可見,目前針對西北干旱區(qū)蒸散發(fā)的研究成果較多,但這些研究選用的站點(diǎn)資料有限,或以整體分析為主,缺乏不同時(shí)空維度的變異性研究,難以反映研究區(qū)特殊的“山地-綠洲-荒漠”系統(tǒng)的異質(zhì)性,對全區(qū)未來蒸散發(fā)變化趨勢的定量研究還相對比較薄弱。研究西北干旱區(qū)蒸散發(fā)的時(shí)空變化規(guī)律和未來的變化趨勢,對提高應(yīng)對氣候變化能力有重要意義,為促進(jìn)區(qū)域經(jīng)濟(jì)社會可持續(xù)發(fā)展提供參考。

    1 研究區(qū)概況

    中國西北干旱區(qū)位于亞洲中部,深居大陸腹地,為多年平均降水量小于200 mm的極端干旱區(qū)和干旱區(qū),屬于典型的大陸性氣候[19],自然條件惡劣。該地區(qū)東以賀蘭山為界,南至昆侖山—阿爾金山—祁連山,北側(cè)和西側(cè)直抵國界,介于73—107°E和35—50°N之間,包括新疆維吾爾自治區(qū)全境、甘肅河西走廊、祁連山區(qū)、內(nèi)蒙古阿拉善高原及黃河寧夏段以西的寧夏自治區(qū)部分,區(qū)內(nèi)有天山、阿爾泰山、昆侖山、祁連山等一系列高大山系,包圍著準(zhǔn)噶爾盆地、塔里木盆地等內(nèi)陸盆地和河西走廊,分布著大片沙漠和戈壁,構(gòu)成以山地-綠洲-荒漠三大生態(tài)系統(tǒng)為基本特征的特殊自然單元。根據(jù)前人的研究成果[19]和區(qū)域自然地理差異,將西北干旱區(qū)分為北疆(新疆北部)、南疆(新疆南部)、天山、祁連山、河西走廊和內(nèi)蒙古西部等6個(gè)亞區(qū)(圖1)。

    圖1 研究區(qū)概況Fig.1 Sketch map of arid region of Northwest China

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

    2.1 數(shù)據(jù)來源

    MODIS ET數(shù)據(jù)集包括地表實(shí)際蒸散發(fā)量(ET)、潛熱通量(LE)、潛在蒸散發(fā)量(PET)和潛在潛熱通量(PLE)產(chǎn)品,空間分辨率為1 km×1 km,時(shí)間分辨率有8天、1月和1年。其算法是Mu等[20]在Penman-Monteith 公式基礎(chǔ)上改進(jìn)的,具體反演流程如圖2。該數(shù)據(jù)集考慮了土壤表面蒸發(fā)、冠層截流水分蒸發(fā)和植物蒸騰,較好地反映了荒漠和綠洲下墊面的非均勻性,適用于干旱區(qū)地表蒸散發(fā)的研究。本文選用地表實(shí)際蒸散發(fā)數(shù)據(jù)的年合成產(chǎn)品(MOD16A3),時(shí)間序列從2000年1月至2014年12月,在美國蒙大拿大學(xué)森林學(xué)院工作組網(wǎng)站下載(http://www.ntsg.umt.edu/project/mod16),使用MRT(Modis Reprojection Tool)工具對數(shù)據(jù)進(jìn)行重投影等預(yù)處理。

    DEM數(shù)據(jù)為SRTM3(Shuttle Radar Topography Mission),空間分辨率為90 m×90 m,來自中國科學(xué)院數(shù)據(jù)云(http://www.csdb.cn/)。

    土地覆蓋數(shù)據(jù)為2001年和2013年的MODIS MCD12Q1 產(chǎn)品(下載地址: https://lpdaac.usgs.gov/),空間分辨率為500 m,該產(chǎn)品采用國際地圈生物圈計(jì)劃全球植被分類方案,將地表類型分為17種,本文將其合并為水體、林地、草地、農(nóng)用地、城市和建筑區(qū)、冰雪與稀疏植被7大類。

    降水資料來自英國East Anglia大學(xué)氣候研究中心(Climatic Research Unit, CRU)發(fā)布的全球陸地表面月平均氣候數(shù)據(jù)集(http://www.uea.ac.uk/),空間分辨率0.5°×0.5°,選取的時(shí)間序列為2000年1月至2014年12月。該資料不包含衛(wèi)星觀測,不使用模式同化,僅用數(shù)學(xué)方法對數(shù)據(jù)進(jìn)行整合和插值[21]。CRU氣候資料盡管包含插值帶來的誤差,但經(jīng)對比,CRU降水除青藏高原西部外與中國臺站的觀測具有很好的一致性[22],故可用于本文的研究。

    徑流數(shù)據(jù)為水文站的觀測,首先考慮河流上游的控制水文站,以盡可能避免選取有大型水庫及大規(guī)模灌溉用地的流域。同時(shí),盡量選擇資料系列較長的流域。按此標(biāo)準(zhǔn),選取研究區(qū)5個(gè)典型流域,包括天山北坡的博爾塔拉河與呼圖壁河,天山南坡的阿克蘇河與開都河,祁連山水系的黑河。降水?dāng)?shù)據(jù)為各流域內(nèi)及周邊24個(gè)氣象站的實(shí)測,所有站點(diǎn)均在圖1中標(biāo)注。

    圖2 MOD16 ET反演算法流程圖Fig.2 Flow chart of MOD16 ET inversion algorithm

    2.2 研究方法

    2.2.1 流域水量平衡法

    采用流域水量平衡法驗(yàn)證MOD16A3產(chǎn)品在西北干旱區(qū)的模擬精度[6]。

    2.2.2 變異系數(shù)法

    變異系數(shù)是描述隨機(jī)變量分散程度的統(tǒng)計(jì)量,用來分析蒸散發(fā)空間格局與空間分異規(guī)律[23]。

    (1)

    2.2.3 Theil-Sen median趨勢分析與Mann-Kendall檢驗(yàn)

    Theil-Sen median趨勢分析與Mann-Kendall檢驗(yàn)相結(jié)合,用以判斷長時(shí)間序列數(shù)據(jù)的趨勢[24]。其中,Theil-Sen median趨勢分析是一種穩(wěn)健的非參數(shù)統(tǒng)計(jì)的趨勢計(jì)算方法,計(jì)算公式為:

    (2)

    式中,ETi和ETj為樣本數(shù)據(jù)值(2000≤i0時(shí),反映了這一時(shí)間序列的蒸散發(fā)量呈增強(qiáng)趨勢;反之,則為衰減趨勢。

    Mann-Kendall屬于非參數(shù)檢驗(yàn)方法,與其他參數(shù)檢驗(yàn)方法相比,不需要樣本遵從一定的分布,也不受少數(shù)異常值干擾,更適合順序變量。Mann-Kendall檢驗(yàn)已經(jīng)在水文、氣象的時(shí)間序列分析中得到成功應(yīng)用[25-26],用于判斷時(shí)間序列數(shù)據(jù)是否具有上升或下降的趨勢。

    在用Mann-Kendall法進(jìn)行ET趨勢檢驗(yàn)時(shí),將某時(shí)序的ET值看作一組獨(dú)立分布的樣本數(shù)據(jù),以參數(shù)Zc作為像元ET衰減指標(biāo),計(jì)算公式如下:

    (3)

    (4)

    (5)

    (6)

    式中,sign為符號函數(shù)。本文在置信水平α=0.05上判斷ET變化趨勢的顯著性[24]。當(dāng)|Zc|>1.96時(shí)表示時(shí)間序列置信水平α<0.05,|Zc|<1.96表示置信水平α>0.05。

    2.2.4 Hurst指數(shù)

    Hurst指數(shù)用于定量描述時(shí)間序列數(shù)據(jù)的可持續(xù)性。在水文學(xué)、經(jīng)濟(jì)學(xué)、氣候?qū)W等領(lǐng)域有著廣泛的應(yīng)用[27]。考慮ET時(shí)間序列{ET(t)} (t=1,2,3,4,…,n) 對于任意正整數(shù)t≥1,定義該時(shí)間序列的均值序列:

    (7)

    計(jì)算累積離差為:

    (8)

    (9)

    (10)

    計(jì)算Hurst指數(shù):

    (11)

    式中,c為常數(shù)。對公式(11)兩邊同時(shí)取對數(shù)即得到Hurst經(jīng)驗(yàn)公式?;跁r(shí)間序列并利用Hurst經(jīng)驗(yàn)公式得到一簇H值進(jìn)行最小二乘法擬合,得出的直線斜率即為修正后的Hurst指數(shù)(H),它揭示了時(shí)間序列的分形特征。

    不同的H對應(yīng)不同的時(shí)間序列趨勢變化:當(dāng)H=0.5時(shí),表明時(shí)間序列是完全獨(dú)立的,沒有相關(guān)性或只是短程相關(guān);00.5時(shí),意味著未來的變化狀況與過去趨勢一致,整個(gè)過程具有持續(xù)性,H越大,持續(xù)性越強(qiáng)。

    3 結(jié)果與分析

    3.1 MOD16A3蒸散發(fā)產(chǎn)品精度評估

    分析表1,水量平衡法計(jì)算的多年蒸散發(fā)量均值與MOD16A3數(shù)據(jù)的年均蒸散發(fā)量較為吻合,二者的平均絕對誤差為38.33 mm,平均相對誤差為12.29%,均方根誤差為39.88 mm。MODIS ET產(chǎn)品的年蒸散發(fā)值總體偏高,但其精度基本滿足區(qū)域尺度的研究,可以用于研究西北干旱區(qū)蒸散發(fā)的時(shí)空動態(tài)特征。

    表1 西北干旱區(qū)典型流域MOD16A3數(shù)據(jù)精度評價(jià)

    3.2 蒸散發(fā)量的空間分布特征

    將2000—2014年的MODIS 蒸散發(fā)量和CRU降水?dāng)?shù)據(jù)逐像元逐年平均(圖3),用以分析ET的空間分布特征。

    圖3 2000—2014年西北干旱區(qū)年均ET和降水量空間分布Fig.3 Spatial distribution of average annual evapotranspiration and precipitation in arid region of Northwest China from 2000 to 2014

    將年均ET值分為8級,進(jìn)行像元統(tǒng)計(jì)分析(表2)。全區(qū)低蒸散發(fā)區(qū)域(ET<200mm)占總面積的38.329%,各亞區(qū)內(nèi)低蒸散發(fā)量區(qū)面積百分比為:內(nèi)蒙古西部(90.197%)>河西走廊(54.379%)>北疆(53.018%)>南疆(50.557%)>天山(10.149%)>祁連山(2.686%)。全區(qū)高蒸散發(fā)區(qū)域(ET>400mm)占總面積的21.834%。

    表2 西北干旱區(qū)年均ET分級

    近15年西北干旱區(qū)ET值總體較小,空間分布上表現(xiàn)為自山區(qū)向兩側(cè)平原減少的特點(diǎn)。這種差異是因?yàn)楦珊祬^(qū)的實(shí)際蒸散發(fā)主要受水分狀況(降水量)影響,降水直接影響地表土壤含水率大小,從而影響蒸散發(fā)量大小[28]。西北干旱區(qū)雖然海拔較高、太陽輻射多,蒸散發(fā)的能量充足,但是地處內(nèi)陸,降水量稀少、土壤濕度低,使得蒸散發(fā)量低值區(qū)面積廣大。分析近15年降水量(圖3)可以看出,區(qū)域東南部及祁連山區(qū)受西南暖濕氣流和東亞季風(fēng)影響,降水量大;西部為西風(fēng)環(huán)流的通道,帶來大西洋的濕潤氣流,在山地迎風(fēng)坡形成豐富的降水,而平原地區(qū)降水量較小。

    土地覆蓋亦影響ET的空間格局。結(jié)合圖4可以看出:ET的高值區(qū)主要分布在天山、阿爾泰山和祁連山等山區(qū)的林地和草地。低值區(qū)主要在南疆塔里木盆地邊緣、北疆準(zhǔn)噶爾盆地、河西走廊的草地和稀疏植被區(qū)。同時(shí),綠洲區(qū)農(nóng)用地的ET高于綠洲邊緣植被稀疏區(qū)。為進(jìn)一步分析不同土地覆蓋類型的蒸散發(fā)量特征,統(tǒng)計(jì)研究區(qū)典型的4種土地覆蓋的ET平均值(圖5),2001年和2013年各土地覆蓋類型的ET平均值為:林地>農(nóng)用地>草地>稀疏植被,這是因?yàn)椴煌恋馗采w的動力和熱力性質(zhì)存在差異,導(dǎo)致地氣相互作用中能量的重新分配[29]。綠洲區(qū)農(nóng)用地因人工種植和灌溉,其植被覆蓋度和土壤水分都高于綠洲邊緣的稀疏植被區(qū),使得蒸散發(fā)表現(xiàn)出顯著的空間異質(zhì)性。

    由于MODIS ET數(shù)據(jù)集對于無植被覆蓋的裸土、沙漠、戈壁等區(qū)域的蒸散發(fā)量不進(jìn)行計(jì)算,故在本研究中將塔克拉瑪干沙漠、古爾班通古特沙漠、巴丹吉林沙漠和戈壁區(qū)域的ET值設(shè)置為NoData,且不計(jì)入面積統(tǒng)計(jì)范圍。

    圖4 2001年和2013年西北干旱區(qū)土地覆蓋的空間變化Fig.4 Spatial distribution of land cover types in arid region of Northwest China in 2001 and 2013

    圖5 2001年和2013年西北干旱區(qū)4種土地覆蓋的ET平均值Fig.5 The average evapotranspiration of four kinds of land cover types in arid region of Northwest China in 2001 and 2013

    3.3 蒸散發(fā)量的區(qū)域分異特征

    運(yùn)用變異系數(shù)法分析西北干旱區(qū)ET的空間格局變異性(圖6),將研究區(qū)像元尺度的變異系數(shù)CVET分為5級(表3)。可以看出,全區(qū)不同變異程度的面積比例為:相對較低的波動變化(40.905%)>中等波動變化(38.549%)>相對較高的波動變化(11.061%)>高波動變化(6.515%)>低波動變化(2.970%),即全區(qū)蒸散發(fā)的變異程度不明顯。各亞區(qū)內(nèi)波動較低區(qū)的面積百分比依次為:北疆(60.482%)>天山(50.667%)>祁連山(25.54%)>內(nèi)蒙古西部(16.238%)>河西走廊(13.619%)>南疆(9.694%)。

    從圖6可以看出,南疆昆侖山北坡蒸散發(fā)呈高波動變化(CVET>0.20),因?yàn)樵搮^(qū)域?qū)儆诘湫偷母吆0胃珊瞪鷳B(tài)脆弱區(qū)[30],分析2000—2014年西北干旱區(qū)和昆侖山區(qū)降水量的變化(圖7)發(fā)現(xiàn),昆侖山區(qū)降水量的變化幅度明顯大于全區(qū)。年降水量變化較大,使得蒸散發(fā)表現(xiàn)出明顯的不穩(wěn)定性;北疆額敏河流域、博爾塔拉河流域、奎屯河流域、瑪納斯河流域、呼圖壁河流域,天山伊犁河谷,南疆塔里木河流域、喀什噶爾河流域有零星的斑點(diǎn)狀高波動變化區(qū)域,這是因?yàn)榫G洲面積擴(kuò)大、種植結(jié)構(gòu)調(diào)整和種植品種變化等人為因素引起綠洲NDVI的波動[31],導(dǎo)致區(qū)域蒸散發(fā)的高波動變化.

    表3 西北干旱區(qū)ET值的變異系數(shù)統(tǒng)計(jì)

    圖6 2000—2014年西北干旱區(qū)ET值變異程度Fig.6 Spatial distribution of variation coefficient of evapotranspiration in arid region of Northwest China from 2000 to 2014

    圖7 2000—2014年西北干旱區(qū)和昆侖山區(qū)年際降水量變化Fig.7 Inter-annual variation of precipitation in arid region of Northwest China and Kunlun Mountain from 2000 to 2014

    3.4 蒸散發(fā)量的時(shí)間變化特征

    為說明西北干旱區(qū)蒸散發(fā)量隨時(shí)間變化的特點(diǎn),將2000—2014年MODIS ET的區(qū)域均值作一元線性回歸分析(圖8)。15年間,全區(qū)蒸散發(fā)量大致穩(wěn)定分布在225—285 mm之間,其中最小值出現(xiàn)在2008年(224.69 mm),最大值出現(xiàn)在2003年(282.13 mm)。各亞區(qū)中,祁連山和天山蒸散發(fā)量相對較高,最大值分別出現(xiàn)在2003年(414.92 mm)和2002年(387.15 mm),最小值分別在2000年(295.12 mm)和2008年(305.96 mm);內(nèi)蒙古西部亞區(qū)蒸散發(fā)量最小,穩(wěn)定在117.71 mm(2013年)至165.90 mm(2003年)之間。

    近15年全區(qū)年均蒸散發(fā)量呈波動變化,總體有微弱的減小趨勢,變化率為-0.9348 mm/a.各亞區(qū)的ET在時(shí)間序列上呈不同程度的減小趨勢,變化率為:天山 (-2.9116 mm/a) >內(nèi)蒙西部 (-1.8266 mm/a) >河西走廊 (-0.573 mm/a) >北疆 (-0.533 mm/a) >祁連山 (-0.2559 mm/a),僅南疆有增加趨勢,變化率為1.4216 mm/a.

    結(jié)合Theil-Sen median趨勢分析與Mann-Kendall檢驗(yàn)的結(jié)果,得到2000—2014年全區(qū)像元尺度的ET變化趨勢(表4、圖9)。在Theil-Sen median趨勢分析中,根據(jù)β值的計(jì)算結(jié)果,分為增強(qiáng)趨勢(β>0)和衰減趨勢(β<0)兩類;在Mann-Kendall檢驗(yàn)中選取顯著性檢驗(yàn)的置信水平為0.05,將結(jié)果劃分為顯著變化(Zc>1.96或Zc<-1.96)和變化不顯著(-1.96≤Zc≤1.96)。

    從表4可以看出:2000—2014年,全區(qū)ET各變化趨勢的面積比例依次為:輕微減小(53.518%)>輕微增加(25.475%)>顯著減小(15.761%)>顯著增加(5.246%),即西北干旱區(qū)蒸散發(fā)變化趨勢以減小為主,各亞區(qū)ET減小程度為:內(nèi)蒙西部(94.363%)>天山(83.022%)>北疆(76.367%)>河西走廊(70.039%)>祁連山(44.556%)>南疆(35.056%)。

    圖8 2000—2014年西北干旱區(qū)年際ET變化Fig.8 Inter-annual variation of evapotranspiration in arid region of Northwest China from 2000 to 2014

    表4 西北干旱區(qū)ET值的變化趨勢統(tǒng)計(jì)

    從變化趨勢分布(圖9)可以看出:北疆額敏河流域、博爾塔拉河流域、奎屯河流域、瑪納斯河流域、呼圖壁河流域,天山伊犁河谷中部,南疆塔里木河流域、喀什噶爾河流域蒸散發(fā)有顯著增加趨勢。這是因?yàn)橐陨狭饔蚓懈呙芏鹊木G洲農(nóng)業(yè)發(fā)展,土地覆蓋的變化影響地表蒸散過程的信息鏈。結(jié)合研究區(qū)代表年份的土地覆蓋類型空間分布(圖4)和轉(zhuǎn)移矩陣(表5),分析區(qū)域蒸散發(fā)量對土地覆蓋變化的響應(yīng)。結(jié)合表6可以看出,2001—2013年,綠洲邊緣的草地和稀疏植被區(qū)演變?yōu)檗r(nóng)用地,導(dǎo)致草地凈減少16569.63 km2,稀疏植被區(qū)凈減少6366.32 km2,而農(nóng)用地凈增加22935.95 km2。稀疏植被演變?yōu)檗r(nóng)用地的區(qū)域,平均蒸散發(fā)量增加了82.41 mm,草地轉(zhuǎn)化為農(nóng)用地的區(qū)域,平均蒸散發(fā)量亦增加了62.77 mm。這是因?yàn)槿斯す喔鹊母仄渫寥罎穸群椭脖桓采w度都高于同氣象條件下的草地和稀疏植被區(qū),故伴隨著農(nóng)用地?cái)U(kuò)張的進(jìn)程,以上流域的蒸散發(fā)都有顯著增加趨勢。

    同時(shí),天山山地和平原過渡帶的蒸散發(fā)為顯著減小趨勢。這是因?yàn)榈乇韺?shí)際蒸散發(fā)亦受植被覆蓋變化的影響,據(jù)相關(guān)研究[32-33],1982—2012年和1982—2013年天山山地和平原過渡帶植被NDVI呈顯著下降趨勢。

    圖9 2000—2014年西北干旱區(qū)ET變化趨勢Fig.9 Spatial distribution of the coefficient of variation of evapotranspiration in arid region of Northwest China from 2000 to 2014

    表5 2001—2013年西北干旱區(qū)土地覆蓋類型轉(zhuǎn)移矩陣/km2

    表6 2001—2013典型土地覆蓋演變區(qū)域蒸散發(fā)量的變化/mm

    3.5 蒸散發(fā)量的預(yù)測

    分析2000—2014年全區(qū)ET的Hurst指數(shù)空間分布圖(圖10),研究ET變化趨勢的復(fù)雜度及對未來變化趨勢進(jìn)行預(yù)測。全區(qū)ET的Hurst指數(shù)均值為0.689,Hurst指數(shù)小于0.5的范圍所占比例為19.67%,大于0.5的范圍所占比例為80.033%,表明西北干旱區(qū)蒸散發(fā)的正向持續(xù)性較強(qiáng),即未來的變化狀況與過去一致。

    圖10 2000—2014年西北干旱區(qū)ET的Hurst指數(shù)分布Fig.10 Spatial distribution of Hurst Index of evapotranspiration in arid region of Northwest China from 2000 to 2014

    為揭示ET未來的變化趨勢及其持續(xù)性,將2000—2014年ET的變化趨勢結(jié)果與Hurst指數(shù)結(jié)果進(jìn)行疊加分析,得到變化趨勢與持續(xù)性的耦合結(jié)果(圖11、表7),全區(qū)各預(yù)測類型的面積比重依次為:持續(xù)性輕微減小(37.922%)>無法確定(22.003%)>持續(xù)性輕微增加(20.783%)>持續(xù)性顯著減小(14.282%)>持續(xù)性顯著增加(5.01%),即西北干旱區(qū)未來ET變化趨勢以持續(xù)性減小為主.各亞區(qū)內(nèi)ET未來的變化趨勢減小程度為:內(nèi)蒙西部(87.265%)>天山(69.888%)>河西走廊(57.542%)>北疆(50.895%)>祁連山(40.483%)>南疆(26.122%)。全區(qū)22.003%的區(qū)域未來變化趨勢無法確定,主要分布在阿爾泰山南坡、準(zhǔn)噶爾盆地北部,該區(qū)域未來蒸散發(fā)的變化狀況需要持續(xù)關(guān)注。

    圖11 西北干旱區(qū)ET的預(yù)測Fig.11 The predictions for the future of evapotranspiration in arid region of Northwest China

    表7 西北干旱區(qū)ET預(yù)測的類型統(tǒng)計(jì)

    4 結(jié)論與討論

    4.1 結(jié)論

    (1)基于水量平衡法評估了MOD16A3產(chǎn)品在西北干旱區(qū)的模擬精度,平均絕對誤差為38.33 mm,平均相對誤差為12.29%,均方根誤差為39.88 mm。MODIS ET產(chǎn)品的模擬值總體偏高,但其精度基本滿足區(qū)域尺度的研究。

    (2)空間格局上,近15年西北干旱區(qū)ET值總體較小,ET小于200 mm的區(qū)域占總面積的38.329%。受降水量影響,蒸散發(fā)量在空間上表現(xiàn)為自山區(qū)向兩側(cè)平原減少的特點(diǎn);受土地覆蓋影響,ET的高值區(qū)(>400 mm)主要在天山、阿爾泰山和祁連山等山區(qū)的林地和草地,ET的低值區(qū)(<200 mm)主要在南疆塔里木盆地邊緣、北疆準(zhǔn)噶爾盆地、河西走廊的草地和稀疏植被區(qū),不同土地覆蓋的蒸散發(fā)量差異顯著:林地>農(nóng)用地>草地>稀疏植被。

    (3)區(qū)域分異特征方面,15年間西北干旱區(qū)蒸散發(fā)變異程度不明顯,各變異程度的比例為:相對較低的波動變化(40.905%)>中等波動變化(38.549%)>相對較高的波動變化(11.061%)>高波動變化(6.515%)>低波動變化(2.970%).各亞區(qū)內(nèi)波動較低區(qū)的比例為:北疆(60.482%)>天山(50.667%)>祁連山(25.54%)>內(nèi)蒙古西部(16.238%)>河西走廊(13.619%)>南疆(9.694%).值得一提的是,昆侖山北坡蒸散發(fā)的高波動變化是年際降水變率較高所致,額敏河流域等的高波動變化則與綠洲面積擴(kuò)大、種植結(jié)構(gòu)調(diào)整和種植品種變化等人為因素有關(guān)。

    (4)在時(shí)間序列上,近15年西北干旱區(qū)年均蒸散發(fā)量呈波動變化,總體有微弱的減小趨勢,變化率為-0.9348mm/a。各亞區(qū)的減小趨勢程度各異:天山 (-2.9116 mm/a) >內(nèi)蒙西部 (-1.8266 mm/a) >河西走廊 (-0.573 mm/a) >北疆 (-0.533 mm/a) >祁連山 (-0.2559 mm/a),僅南疆有增加趨勢,變化率為1.4216 mm/a。

    基于像元尺度的分析也表明全區(qū)ET以減小的變化趨勢為主,各變化比例為:輕微減小(53.518%)>輕微增加(25.475%)>顯著減小(15.761%)>顯著增加(5.246%)。其中,北疆瑪納斯河流域、天山伊犁河谷中部、南疆塔里木河流域等蒸散發(fā)有顯著增加趨勢,與土地覆蓋演變有關(guān)。天山山地和平原過渡帶的蒸散發(fā)為顯著減小趨勢,是受植被覆蓋變化的影響。

    (5)對未來趨勢的預(yù)測方面,西北干旱區(qū)ET的Hurst指數(shù)均值為0.689,Hurst指數(shù)大于0.5的范圍所占比例為80.033%,表明全區(qū)蒸散發(fā)未來的變化狀況與過去一致。具體來講,各預(yù)測類型的面積比重為:持續(xù)性輕微減小(37.922%)>無法確定(22.003%)>持續(xù)性輕微增加(20.783%)>持續(xù)性顯著減小(14.282%)>持續(xù)性顯著增加(5.01%),即未來全區(qū)的ET變化趨勢以持續(xù)性減小為主。其中22.003%區(qū)域的變化趨勢無法確定,主要分布在阿爾泰山南坡、準(zhǔn)噶爾盆地北部,該區(qū)域未來蒸散發(fā)的變化狀況需要持續(xù)關(guān)注。各亞區(qū)的空間異質(zhì)性分析中,減小趨勢為:內(nèi)蒙西部(87.265%)>天山(69.888%)>河西走廊(57.542%)>北疆(50.895%)>祁連山(40.483%)>南疆(26.122%)。

    4.2 討論

    (1)MODIS ET數(shù)據(jù)集的反演算法考慮了土壤表面蒸發(fā)、冠層截流水分蒸發(fā)和植物蒸騰,較好地反映了荒漠和綠洲下墊面的非均勻性。通過水量平衡法的驗(yàn)證,該產(chǎn)品的模擬值總體偏高,但其精度基本滿足區(qū)域尺度的研究。因此,高空間分辨率的MODIS ET數(shù)據(jù)集可以被用于揭示區(qū)域蒸散發(fā)的時(shí)空動態(tài)特征,尤其在觀測站點(diǎn)稀疏的西北干旱區(qū)。

    (2)西北干旱區(qū)蒸散發(fā)的空間格局受降水和土地覆蓋的綜合影響。地處內(nèi)陸,降水稀少,決定了實(shí)際蒸散發(fā)低值區(qū)面積廣大。祁連山和天山等高大山脈迎風(fēng)坡有豐富的降水,而周邊盆地和走廊降水稀少,使蒸散發(fā)量表現(xiàn)出山區(qū)大于平原的特點(diǎn);同時(shí),由于動力和熱力性質(zhì)差異,不同土地覆蓋的蒸散發(fā)量差異顯著:林地>農(nóng)用地>草地>稀疏植被。王海波等[12]在西北干旱區(qū)典型內(nèi)陸河黑河流域,模擬了高寒草地和干旱區(qū)農(nóng)田的蒸散發(fā),得出生長季農(nóng)田蒸散發(fā)量大于草地的結(jié)論。田靜等[34]利用NOAH陸面過程模擬了中國陸地蒸散發(fā),認(rèn)為年降雨量是決定蒸散發(fā)量大小的主要因素,西北地區(qū)不同土地利用類型的蒸散發(fā)量為:林地>耕地>草地>未利用地,與本研究的結(jié)論相似。因此,本研究的結(jié)論對認(rèn)識西北干旱區(qū)不同空間維度的蒸散發(fā)格局增添了新的證據(jù)。

    (3)受人類活動和氣候變化的共同影響,2000—2014年全區(qū)蒸散發(fā)以相對較低的波動變化與中等波動變化為主。西北干旱區(qū)內(nèi)陸河流域發(fā)展了典型的綠洲農(nóng)業(yè),土地覆蓋演變、綠洲面積擴(kuò)大、種植結(jié)構(gòu)調(diào)整和種植品種變化均會引起蒸散發(fā)的波動;除人類活動,氣候變化對蒸散發(fā)的分異亦產(chǎn)生深刻影響,尤其是生態(tài)環(huán)境脆弱的高寒山區(qū)。氣候變化背景下的地表蒸散發(fā)響應(yīng)[15,17- 18]和中國陸地蒸散發(fā)的時(shí)間序列變化[34-35]已有較多探討。與以上研究相比,本研究有以下優(yōu)勢:根據(jù)自然地理差異,將西北干旱區(qū)分為6個(gè)亞區(qū),分析不同尺度的蒸散發(fā)時(shí)間變化和分異特征,揭示了空間異質(zhì)性,得出蒸散發(fā)波動變化在6個(gè)亞區(qū)內(nèi)程度不一的結(jié)論,北疆和天山以相對較低的波動變化為主,南疆、祁連山、河西走廊和內(nèi)蒙古西部以中等波動變化占主導(dǎo)。西北干旱區(qū)“山地—綠洲—荒漠”系統(tǒng)自然要素的分異特征十分鮮明,在全球干旱區(qū)具有很強(qiáng)的代表性。本研究的結(jié)論,對認(rèn)識干旱區(qū)內(nèi)部不同干濕背景下蒸散發(fā)的異質(zhì)性提供了參考。

    (4)2000年以來西北干旱區(qū)大部分區(qū)域蒸散發(fā)有微弱的減小趨勢,變化率為-0.9348 mm/a,這與區(qū)域降水變化有一定聯(lián)系。據(jù)陳亞寧等[36]的研究,最近10年,西北干旱區(qū)降水量的增加幅度降低,并且約有45%臺站的降水量較20世紀(jì)90年代表現(xiàn)為減少趨勢。降水量降低使土壤水分減少,導(dǎo)致西北干旱區(qū)實(shí)際蒸散發(fā)量減少。本研究與Nature[3]上發(fā)表的由于水分供給不足導(dǎo)致的全球陸面蒸散發(fā)呈下降趨勢的結(jié)論是一致的。與楊秀芹等[35]利用GLEAM遙感模型計(jì)算中國1980—2011年地表蒸散發(fā),得到的西北西部地表蒸散發(fā)量存在下降趨勢有較好的一致性。本研究對西北干旱區(qū)蒸散發(fā)的變化趨勢增加了定量的分析,為認(rèn)識全球變化和人類活動共同驅(qū)動下的干旱區(qū)蒸散發(fā)變化研究提供了一定參考。

    (5)未來,全區(qū)蒸散發(fā)持續(xù)減小的面積比重近一半,各亞區(qū)的減小趨勢依次為:內(nèi)蒙古西部>天山>河西走廊>北疆>祁連山>南疆。這種持續(xù)減小的變化趨勢,對綠洲、沙塵暴和湖泊等干旱區(qū)特殊的生態(tài)環(huán)境要素造成影響[10]。蒸散發(fā)量的減小,對高山草原和綠洲灌溉農(nóng)業(yè)等有有利的影響,同時(shí)也能夠緩解荒漠綠洲過渡帶的退化。隨著區(qū)域氣候系統(tǒng)的變化和人類活動干預(yù),蒸散發(fā)減小的趨勢是否會持續(xù)發(fā)展,干旱區(qū)水文循環(huán)、生態(tài)系統(tǒng)和地表過程將如何響應(yīng),尤其是對荒漠生態(tài)環(huán)境的保護(hù)和修復(fù)將產(chǎn)生怎樣的影響,需要全面的、長序列的數(shù)據(jù)進(jìn)行深入研究。

    [1] Xiong Y J, Zhao S H, Tian F, Qiu G Y. An evapotranspiration product for arid regions based on the three-temperature model and thermal remote sensing. Journal of Hydrology, 2015, 530: 392- 404.

    [2] Zhang K X, Pan S M, Zhang W, Xu Y H, Cao L G, Hao Y P, Wang Y. Influence of climate change on reference evapotranspiration and aridity index and their temporal-spatial variations in the Yellow River Basin, China, from 1961 to 2012. Quaternary International, 2015, 380- 381: 75- 82.

    [3] Jung M, Reichstein M, Ciais P, Seneviratne S I, Sheffield J, Goulden M L, Bonan G, Cescatti A, Chen J Q, de Jeu R, Dolman A J, Eugster W, Gerten D, Gianelle D, Gobron N, Heinke J, Kimball J, Law B E, Montagnani L, Mu Q Z, Mueller B, Oleson K, Papale D, Richardson A D, Roupsard O, Running S, Tomelleri E, Viovy N, Weber U, Williams C, Wood E, Zaehle S, Zhang K. Recent decline in the global land evapotranspiration trend due to limited moisture supply. Nature, 2010, 467(7318): 951- 954.

    [4] Ghilain N, Arboleda A, Gellens-Meulenberghs F. Evapotranspiration modelling at large scale using near-real time MSG SEVIRI derived data. Hydrology and Earth System Sciences, 2011, 15(3): 771- 786.

    [5] Mu Q Z, Heinsch F A, Zhao M S, Running S W. Development of a global evapotranspiration algorithm based on MODIS and global meteorology data. Remote Sensing of Environment, 2007, 111(4): 519- 536.

    [6] 賀添, 邵全琴. 基于MOD16產(chǎn)品的我國2001—2010年蒸散發(fā)時(shí)空格局變化分析. 地球信息科學(xué)學(xué)報(bào), 2014, 16(6): 979- 988.

    [7] Kim H W, Hwang K, Mu Q Z, Lee S O, Choi M. Validation of MODIS 16 global terrestrial evapotranspiration products in various climates and land cover types in Asia. KSCE Journal of Civil Engineering, 2012, 16(2): 229- 238.

    [8] Liu S M, Xu Z W, Zhu Z L, Jia Z Z, Zhu M J. Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai River Basin, China. Journal of Hydrology, 2013, 487: 24- 38.

    [9] Jang K, Kang S, Lim Y J, Jeong S, Kim J, Kimball J S, Hong S Y. Monitoring daily evapotranspiration in Northeast Asia using MODIS and a regional land data assimilation system. Journal of Geophysical Research: Atmospheres, 2013, 118(23): 12927- 12940.

    [10] 姚俊強(qiáng), 楊青, 陳亞寧, 胡文峰, 劉志輝, 趙玲. 西北干旱區(qū)氣候變化及其對生態(tài)環(huán)境影響. 生態(tài)學(xué)雜志, 2013, 32(5): 1283- 1291.

    [11] 李寶富, 陳亞寧, 李衛(wèi)紅, 曹志超. 基于遙感和SEBAL模型的塔里木河干流區(qū)蒸散發(fā)估算. 地理學(xué)報(bào), 2011, 66(9): 1230- 1238.

    [12] 王海波, 馬明國. 基于遙感和Penman-Monteith模型的內(nèi)陸河流域不同生態(tài)系統(tǒng)蒸散發(fā)估算. 生態(tài)學(xué)報(bào), 2014, 34(19): 5617- 5626.

    [13] 馬金龍, 劉麗娟, 李小玉, 王進(jìn), 楊會巾. 干旱區(qū)綠洲膜下滴灌棉田蒸散過程. 生態(tài)學(xué)雜志, 2015, 34(4): 974- 981.

    [14] 張鑫, 佟玲, 李思恩, 康紹忠. 干旱區(qū)兩種微氣象學(xué)法測定農(nóng)田蒸散發(fā)的比較研究. 水利學(xué)報(bào), 2011, 42(12): 1470- 1478.

    [15] 劉春雨, 董曉峰, 劉英英. 西北干旱區(qū)遙感ET與潛在ET對氣候變化的響應(yīng)——以甘南草原區(qū)域?yàn)槔? 蘭州大學(xué)學(xué)報(bào): 自科科學(xué)版, 2014, 50(2): 194- 199, 207- 207.

    [16] Li Z, Chen Y N, Shen Y J, Liu Y B, Zhang S H. Analysis of changing pan evaporation in the arid region of Northwest China. Water Resources Research, 2013, 49(4): 2205- 2212.

    [17] Penman H L. Natural evaporation from open water, bare soil and grass. Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 1948, 193(1032): 120- 145.

    [18] Bouchet R J. Evapotranspiration réelle et potentielle, signification climatique. General Assembly of Berkeley, Red Book, 1963, 62: 134- 142.

    [19] 姚俊強(qiáng), 楊青, 劉志輝, 李誠志. 中國西北干旱區(qū)降水時(shí)空分布特征. 生態(tài)學(xué)報(bào), 2015, 35(17): 5846- 5855.

    [20] Mu Q Z, Zhao M S, Running S W. Improvements to a MODIS global terrestrial evapotranspiration algorithm. Remote Sensing of Environment, 2011, 115(8): 1781- 1800.

    [21] 聞新宇, 王紹武, 朱錦紅, Viner D. 英國CRU高分辨率格點(diǎn)資料揭示的20世紀(jì)中國氣候變化. 大氣科學(xué), 2006, 30(5): 894- 904.

    [22] 馬柱國, 符淙斌. 中國干旱和半干旱帶的10年際演變特征. 地球物理學(xué)報(bào), 2005, 48(3): 519- 525.

    [23] Milich L, Weiss E. GAC NDVI interannual coefficient of variation (CoV) images: ground truth sampling of the Sahel along north-south transects. International Journal of Remote Sensing, 2000, 21(2): 235- 260.

    [24] Jiang W G, Yuan L H, Wang W J, Cao R, Zhang Y F, Shen W M. Spatio-temporal analysis of vegetation variation in the Yellow River Basin. Ecological Indicators, 2015, 51: 117- 126.

    [25] Yue S, Pilon P, Cavadias G. Power of the Mann-Kendall and Spearman′s rho tests for detecting monotonic trends in hydrological series. Journal of Hydrology, 2002, 259(1/4): 254- 271.

    [26] Fensholt R, Langanke T, Rasmussen K, Reenberg A, Prince S D, Tucker C, Scholes R J, Le Q B, Bondeau A, Eastman R, Epstein H, Gaughan A E, Hellden U, Mbow C, Olsson L, Paruelo J, Schweitzer C, Seaquist J, Wessels K. Greenness in semi-arid areas across the globe 1981—2007—an earth observing satellite based analysis of trends and drivers. Remote Sensing of Environment, 2012, 121: 144- 158.

    [27] Jiapaer G, Liang S L, Yi Q X, Liu J P. Vegetation dynamics and responses to recent climate change in Xinjiang using leaf area index as an indicator. Ecological Indicators, 2015, 58: 64- 76.

    [28] Roderick M L, Hobbins M T, Farquhar G D. Pan evaporation trends and the terrestrial water balance. II. Energy balance and interpretation. Geography Compass, 2009, 3(2): 761- 780.

    [29] 何延波, Su Z, Jia L, 王石立. 遙感數(shù)據(jù)支持下不同地表覆蓋的區(qū)域蒸散. 應(yīng)用生態(tài)學(xué)報(bào), 2007, 18(2): 288- 296.

    [30] 朱軍濤, 李向義, 張希明, 曾凡江, 林麗莎, 楊尚功, 桂東偉, 王輝, 劉波. 昆侖山北坡前山帶塔里木沙拐棗對不同海拔生境的生理生態(tài)響應(yīng). 生態(tài)學(xué)報(bào), 2010, 30(3): 602- 609.

    [31] 郭鈮, 朱燕君, 王介民, 鄧朝平. 近22年來西北不同類型植被NDVI變化與氣候因子的關(guān)系. 植物生態(tài)學(xué)報(bào), 2008, 32(2): 319- 327.

    [32] 劉憲鋒, 朱秀芳, 潘耀忠, 李宜展, 趙安周. 1982—2012年中國植被覆蓋時(shí)空變化特征. 生態(tài)學(xué)報(bào), 2015, 35(16): 5331- 5342.

    [33] 許玉鳳, 楊井, 陳亞寧, 楊雅雪. 近32年來新疆地區(qū)植被覆蓋的時(shí)空變化. 草業(yè)科學(xué), 2015, 32(5): 702- 709.

    [34] 田靜, 蘇紅波, 陳少輝, 孫曉敏, 陳慶美. 近20年來中國內(nèi)陸地表蒸散的時(shí)空變化. 資源科學(xué), 2012, 34(7): 1277- 1286.

    [35] 楊秀芹, 王國杰, 潘欣, 張余慶. 基于GLEAM遙感模型的中國1980—2011年地表蒸散發(fā)時(shí)空變化. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(21): 132- 141.

    [36] 陳亞寧, 李稚, 范煜婷, 王懷軍, 方功煥. 西北干旱區(qū)氣候變化對水文水資源影響研究進(jìn)展. 地理學(xué)報(bào), 2014, 69(9): 1295- 1304.

    Temporal-spatial dynamic change characteristics of evapotranspiration in arid region of Northwest China

    DENG Xingyao1, 2, 3, LIU Yang1, 2, 3, LIU Zhihui2, 3,*, YAO Junqiang4

    1CollegeofResourceandEnvironmentScience,XinjiangUniversity,Urumqi830046,China2InstituteofAridEcologyandEnvironment,XinjiangUniversity,Urumqi830046,China3KeyLaboratoryofOasisEcologyofEducationMinistry,XinjiangUniversity,Urumqi830046,China4InstituteofDesertMeteorology,ChinaMeteorlogicalAdministration,Urumqi830002,China

    Using MODIS ET data of actual surface evapotranspiration products that were concentrated between 2000 to 2014, we used variable coefficients, the Theil-Sen median trend analysis, Mann-Kendall test, and Hurst index, to investigate the spatial pattern of evapotranspiration, spatial heterogeneity of different dimensionalities, characteristics of time variation, and future trends considering an arid region in Northwest China. The results showed the following: 1) Evapotranspiration over the entire region was very low from 2000 to 2014, and the area of evapotranspiration, less than 200 mm, accounts for 38.329% of the total area. Evapotranspiration decreases from mountainous areas to the plains. Evapotranspiration data were obtained for different land covers, which are forest land, cropland, grassland, and sparse vegetation, arranged from high to low. Mountain forests and grasslands that receive high rainfall have higher values of evapotranspiration (more than 400 mm) than sparse vegetation plains and grasslands that experience low rainfall (less than 200 mm). 2) The degree of variation in evapotranspiration over the past 15 years for the entire region is not obvious, but shows a slight fluctuation. The proportion of areas with slight fluctuations in each subregion ranged from large to small in the order of Northern Xinjiang, Tianshan Mountain, Qilian Mountain, Western Inner Mongolia, Hexi Corridor, and Southern Xinjiang. 3) The degree of variation in evapotranspiration over the past 15 years for the entire region shows changes in the fluctuations with a weakly decreasing trend, at a change rate of-0.9348 mm/a. An analysis based on the pixel scale also shows a mainly decreasing trend. The degree of decreases in each subregion ranged from large to small in the order of Tianshan Mountain, Western Inner Mongolia, Hexi Corridor, Northern Xinjiang, and Qilian Mountain. Only Southern Xinjiang showed an increasing tendency. 4) The Hurst index average of evapotranspiration for the entire region is 0.689. The area of the Hurst index greater than 0.5 accounts for 80.033% of the total area. The trend of changes in evapotranspiration in the future for the entire region is mainly towards a persistent decrease. However, the trend of changes for 22.003% of the area cannot be determined. The future trend of change in the degree of decreases for each subregion from large to small follows the order of Western Inner Mongolia, Tianshan Mountain, Hexi Corridor, Northern Xinjiang, Qilian Mountain, and Southern Xinjiang.

    MODIS; evapotranspiration; the arid region of Northwest China; trend analysis; climate change

    國家科技支撐計(jì)劃項(xiàng)目課題(2012BAC23B01);水利部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201301103);新疆維吾爾自治區(qū)研究生科研創(chuàng)新項(xiàng)目(XJGRI2015019);國家國際科技合作項(xiàng)目(2010DFA92720-12)

    2016- 01- 27; 網(wǎng)絡(luò)出版日期:2016- 12- 19

    10.5846/stxb201601270190

    *通訊作者Corresponding author.E-mail: lzh@xju.edu.cn

    鄧興耀,劉洋,劉志輝,姚俊強(qiáng).中國西北干旱區(qū)蒸散發(fā)時(shí)空動態(tài)特征.生態(tài)學(xué)報(bào),2017,37(9):2994- 3008.

    Deng X Y, Liu Y, Liu Z H, Yao J Q.Temporal-spatial dynamic change characteristics of evapotranspiration in arid region of Northwest China.Acta Ecologica Sinica,2017,37(9):2994- 3008.

    猜你喜歡
    干旱區(qū)西北趨勢
    黑龍江半干旱區(qū)飼用谷子高產(chǎn)栽培技術(shù)
    趨勢
    一座西北小城
    干旱區(qū)生態(tài)修復(fù)的實(shí)踐——以古爾班通古特沙漠為例
    科學(xué)(2020年6期)2020-02-06 08:59:54
    西北不惑
    黔西北行吟
    散文詩(2017年15期)2018-01-19 03:07:58
    初秋唇妝趨勢
    Coco薇(2017年9期)2017-09-07 21:23:49
    SPINEXPO?2017春夏流行趨勢
    氣候變化背景下西北干旱區(qū)旱澇的變化規(guī)律
    干旱區(qū)影響水面蒸發(fā)的氣象因素多元回歸分析
    久久久久久久大尺度免费视频| 久久久午夜欧美精品| 在线观看免费日韩欧美大片 | .国产精品久久| 这个男人来自地球电影免费观看 | 哪个播放器可以免费观看大片| 久久久久视频综合| 精品少妇久久久久久888优播| 国产精品99久久99久久久不卡 | 高清毛片免费看| 在线观看一区二区三区激情| 亚洲av不卡在线观看| 国内精品宾馆在线| 精品亚洲乱码少妇综合久久| 国产毛片在线视频| 日本与韩国留学比较| 色哟哟·www| 久久国内精品自在自线图片| 高清毛片免费看| 亚洲国产日韩一区二区| 自拍偷自拍亚洲精品老妇| 国产国拍精品亚洲av在线观看| 97热精品久久久久久| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费| kizo精华| 免费人妻精品一区二区三区视频| 国产极品天堂在线| 丰满乱子伦码专区| 免费在线观看成人毛片| 欧美老熟妇乱子伦牲交| 欧美日韩视频精品一区| av国产精品久久久久影院| 美女高潮的动态| 亚洲精华国产精华液的使用体验| 久久久亚洲精品成人影院| 国产国拍精品亚洲av在线观看| 超碰av人人做人人爽久久| 嫩草影院新地址| 久久ye,这里只有精品| 久久久久视频综合| 日韩亚洲欧美综合| av专区在线播放| 成人毛片60女人毛片免费| 亚州av有码| 夜夜骑夜夜射夜夜干| 极品少妇高潮喷水抽搐| 天堂中文最新版在线下载| 视频区图区小说| av免费在线看不卡| 国内少妇人妻偷人精品xxx网站| 日本黄大片高清| 午夜福利在线观看免费完整高清在| 永久免费av网站大全| 欧美精品一区二区免费开放| 久久99热这里只有精品18| 亚洲自偷自拍三级| 国产精品精品国产色婷婷| 欧美xxxx性猛交bbbb| 国产av一区二区精品久久 | 午夜免费鲁丝| 亚洲av综合色区一区| 麻豆成人av视频| 国产永久视频网站| av又黄又爽大尺度在线免费看| 美女福利国产在线 | 亚洲激情五月婷婷啪啪| 天堂8中文在线网| 超碰97精品在线观看| 中文资源天堂在线| 日本午夜av视频| 黄色怎么调成土黄色| 成人二区视频| 亚洲在久久综合| 欧美日韩视频精品一区| 青青草视频在线视频观看| 卡戴珊不雅视频在线播放| 久久人人爽av亚洲精品天堂 | 伊人久久精品亚洲午夜| 日本免费在线观看一区| 插逼视频在线观看| av天堂中文字幕网| 女人久久www免费人成看片| 亚洲欧美日韩东京热| 永久免费av网站大全| 日韩中字成人| 亚洲国产欧美人成| 欧美zozozo另类| 久久久久久久国产电影| 亚洲av日韩在线播放| 少妇的逼好多水| 乱系列少妇在线播放| 人妻 亚洲 视频| 亚洲av中文av极速乱| 中国国产av一级| 韩国高清视频一区二区三区| 国产探花极品一区二区| 男女下面进入的视频免费午夜| 91精品伊人久久大香线蕉| 精品亚洲乱码少妇综合久久| av.在线天堂| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 少妇的逼水好多| 只有这里有精品99| 美女高潮的动态| 亚洲精品一区蜜桃| 免费黄频网站在线观看国产| 久久久精品免费免费高清| 中国美白少妇内射xxxbb| 国产黄片美女视频| 婷婷色麻豆天堂久久| 亚洲怡红院男人天堂| 久久久精品94久久精品| 我要看黄色一级片免费的| 一级爰片在线观看| 久久99蜜桃精品久久| 久久青草综合色| 久久精品国产亚洲av天美| 人体艺术视频欧美日本| 91久久精品国产一区二区三区| 日韩一区二区三区影片| 黄色日韩在线| 国产欧美另类精品又又久久亚洲欧美| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 午夜福利在线在线| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 深夜a级毛片| 久久国产亚洲av麻豆专区| 又爽又黄a免费视频| 欧美xxⅹ黑人| 免费观看性生交大片5| 国产乱人偷精品视频| 日本黄大片高清| 伦理电影免费视频| 久久精品国产a三级三级三级| 你懂的网址亚洲精品在线观看| 亚洲内射少妇av| a级毛色黄片| 国产成人精品福利久久| 一边亲一边摸免费视频| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花 | 亚洲av成人精品一二三区| 免费观看的影片在线观看| h视频一区二区三区| 亚洲美女黄色视频免费看| 王馨瑶露胸无遮挡在线观看| h视频一区二区三区| 91aial.com中文字幕在线观看| 久久韩国三级中文字幕| 国产av一区二区精品久久 | 久久精品夜色国产| 黄色日韩在线| 在线观看三级黄色| 欧美日韩亚洲高清精品| 久久热精品热| 在线观看免费高清a一片| 国产成人精品一,二区| 亚洲精品日韩在线中文字幕| 亚洲,一卡二卡三卡| 中文字幕人妻熟人妻熟丝袜美| 国产女主播在线喷水免费视频网站| 久久ye,这里只有精品| 久久久久久人妻| 日本黄大片高清| 18禁在线播放成人免费| 成人综合一区亚洲| 在线观看一区二区三区| 97精品久久久久久久久久精品| 日韩国内少妇激情av| 成人国产麻豆网| 亚洲aⅴ乱码一区二区在线播放| av女优亚洲男人天堂| 午夜激情福利司机影院| 久久精品国产a三级三级三级| 国产精品99久久久久久久久| 纵有疾风起免费观看全集完整版| 久久久久久久久大av| 久久6这里有精品| 国产精品久久久久久精品古装| 亚洲人成网站在线观看播放| 边亲边吃奶的免费视频| 国产精品人妻久久久影院| 国语对白做爰xxxⅹ性视频网站| 亚洲美女搞黄在线观看| 中文资源天堂在线| 男女边摸边吃奶| 一边亲一边摸免费视频| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线| 久久影院123| 欧美成人午夜免费资源| freevideosex欧美| 日韩免费高清中文字幕av| 欧美日本视频| 我的女老师完整版在线观看| 午夜激情福利司机影院| 成年美女黄网站色视频大全免费 | 少妇的逼好多水| 国产精品蜜桃在线观看| 欧美精品人与动牲交sv欧美| 91久久精品国产一区二区成人| 免费少妇av软件| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡 | 一个人看视频在线观看www免费| 哪个播放器可以免费观看大片| 国产精品爽爽va在线观看网站| 亚洲美女黄色视频免费看| 天堂8中文在线网| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久成人av| 秋霞伦理黄片| 国产免费福利视频在线观看| 一级毛片电影观看| 少妇 在线观看| 亚洲,欧美,日韩| 精品国产乱码久久久久久小说| 国产免费一区二区三区四区乱码| 国产在线一区二区三区精| 国产精品一区二区在线观看99| 久久久久精品久久久久真实原创| 嫩草影院新地址| 欧美 日韩 精品 国产| 欧美+日韩+精品| 成人黄色视频免费在线看| 免费大片18禁| 日韩一区二区三区影片| 中文资源天堂在线| 国产精品偷伦视频观看了| 深爱激情五月婷婷| 51国产日韩欧美| 精品久久久噜噜| 人妻夜夜爽99麻豆av| 国产白丝娇喘喷水9色精品| 亚洲欧美精品专区久久| 久久ye,这里只有精品| 欧美精品亚洲一区二区| 成人综合一区亚洲| 国产精品久久久久久久电影| 亚洲伊人久久精品综合| a 毛片基地| www.av在线官网国产| 一级a做视频免费观看| 欧美高清性xxxxhd video| 亚洲激情五月婷婷啪啪| 久久久色成人| 老司机影院毛片| 在线观看国产h片| 丝袜喷水一区| 久久久久久伊人网av| 亚洲av综合色区一区| 国产免费一级a男人的天堂| 国产成人a∨麻豆精品| 欧美97在线视频| 久久精品久久久久久噜噜老黄| 午夜激情福利司机影院| 亚洲av欧美aⅴ国产| 国产免费一区二区三区四区乱码| 欧美激情国产日韩精品一区| 亚洲国产欧美人成| 国产色爽女视频免费观看| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看| 熟女电影av网| 天美传媒精品一区二区| 国产又色又爽无遮挡免| 精品久久久久久久久av| 精品一区二区三卡| 日韩中文字幕视频在线看片 | 国产精品福利在线免费观看| 亚洲国产色片| 精品熟女少妇av免费看| 美女内射精品一级片tv| 伦理电影大哥的女人| 国产精品爽爽va在线观看网站| 亚洲第一av免费看| 丰满人妻一区二区三区视频av| 在线亚洲精品国产二区图片欧美 | 人妻夜夜爽99麻豆av| 六月丁香七月| 五月开心婷婷网| 搡老乐熟女国产| 一级爰片在线观看| 日韩三级伦理在线观看| 国产深夜福利视频在线观看| 亚洲成人中文字幕在线播放| 久久久久性生活片| 国产高潮美女av| 97热精品久久久久久| 人妻制服诱惑在线中文字幕| 99九九线精品视频在线观看视频| 久久这里有精品视频免费| av免费在线看不卡| 老师上课跳d突然被开到最大视频| 欧美日韩视频精品一区| 干丝袜人妻中文字幕| 午夜福利高清视频| 乱系列少妇在线播放| 天堂中文最新版在线下载| 亚洲高清免费不卡视频| 亚洲av二区三区四区| 校园人妻丝袜中文字幕| 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| 美女中出高潮动态图| 一级毛片 在线播放| 亚洲av不卡在线观看| 街头女战士在线观看网站| 97精品久久久久久久久久精品| 国产免费福利视频在线观看| 国产在视频线精品| 国产av码专区亚洲av| 久久人人爽av亚洲精品天堂 | 一区二区三区精品91| 男女国产视频网站| 美女高潮的动态| 国产精品一区二区三区四区免费观看| 成人国产麻豆网| 精品一区二区三区视频在线| 80岁老熟妇乱子伦牲交| 欧美精品一区二区免费开放| 亚洲婷婷狠狠爱综合网| 亚洲第一av免费看| 黄片无遮挡物在线观看| 啦啦啦中文免费视频观看日本| 国产日韩欧美亚洲二区| 成人毛片a级毛片在线播放| 免费大片18禁| 新久久久久国产一级毛片| 22中文网久久字幕| 免费观看性生交大片5| 欧美高清成人免费视频www| 大码成人一级视频| 男女无遮挡免费网站观看| 亚洲精品色激情综合| 精品人妻一区二区三区麻豆| 少妇人妻 视频| 欧美成人精品欧美一级黄| 亚洲美女黄色视频免费看| 天堂俺去俺来也www色官网| 亚洲一区二区三区欧美精品| 简卡轻食公司| 亚洲精品自拍成人| 久久99热6这里只有精品| 午夜福利在线在线| 国产无遮挡羞羞视频在线观看| 国产免费一级a男人的天堂| 91精品国产国语对白视频| 中国三级夫妇交换| 亚洲精品成人av观看孕妇| 国产成人freesex在线| 99久久中文字幕三级久久日本| 99热6这里只有精品| 99精国产麻豆久久婷婷| 99九九线精品视频在线观看视频| 成年女人在线观看亚洲视频| 最近最新中文字幕免费大全7| 欧美高清性xxxxhd video| 久久久午夜欧美精品| 下体分泌物呈黄色| 五月开心婷婷网| 青青草视频在线视频观看| 亚洲国产毛片av蜜桃av| 深夜a级毛片| 亚洲欧美日韩另类电影网站 | 黑人高潮一二区| 欧美成人午夜免费资源| 在线观看免费日韩欧美大片 | 国产成人91sexporn| 性色av一级| 国产免费福利视频在线观看| 卡戴珊不雅视频在线播放| 久久久欧美国产精品| 成人国产av品久久久| 尾随美女入室| 欧美精品国产亚洲| 亚洲不卡免费看| 国产精品成人在线| 美女cb高潮喷水在线观看| 亚洲在久久综合| 一级毛片电影观看| 少妇精品久久久久久久| 婷婷色综合大香蕉| 最近2019中文字幕mv第一页| 大香蕉久久网| 欧美另类一区| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| 乱系列少妇在线播放| 亚洲欧美精品专区久久| 亚洲精品国产av成人精品| 日韩国内少妇激情av| 97在线视频观看| 国产色爽女视频免费观看| 精品久久久久久久久亚洲| 99热6这里只有精品| 国产免费福利视频在线观看| 亚洲婷婷狠狠爱综合网| 精品酒店卫生间| 日本av手机在线免费观看| 亚洲人成网站在线观看播放| 日韩av免费高清视频| 老师上课跳d突然被开到最大视频| 草草在线视频免费看| 一级av片app| 亚洲一区二区三区欧美精品| 国内精品宾馆在线| 麻豆成人av视频| 国产一区亚洲一区在线观看| 99热这里只有是精品在线观看| 18+在线观看网站| 久久久久久伊人网av| 亚洲精品第二区| 我的老师免费观看完整版| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 日本一二三区视频观看| 天天躁夜夜躁狠狠久久av| 青春草国产在线视频| 国产 一区精品| 中文字幕制服av| 99久久精品热视频| 国产免费又黄又爽又色| 在线观看国产h片| 亚洲最大成人中文| 午夜老司机福利剧场| 精品熟女少妇av免费看| 少妇 在线观看| 国产精品无大码| 狠狠精品人妻久久久久久综合| 精品一区二区免费观看| 成人一区二区视频在线观看| 青春草亚洲视频在线观看| 亚洲最大成人中文| 最近手机中文字幕大全| 国产爱豆传媒在线观看| 免费看av在线观看网站| 草草在线视频免费看| 亚洲国产毛片av蜜桃av| 国产精品爽爽va在线观看网站| 亚洲在久久综合| 六月丁香七月| 制服丝袜香蕉在线| 三级国产精品片| 日本午夜av视频| 精品人妻一区二区三区麻豆| 色哟哟·www| 麻豆成人av视频| 欧美成人a在线观看| 国产精品久久久久久久电影| 狂野欧美白嫩少妇大欣赏| 国产精品偷伦视频观看了| 精品99又大又爽又粗少妇毛片| 日日撸夜夜添| 另类亚洲欧美激情| 国产高清有码在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 欧美xxxx性猛交bbbb| 亚洲天堂av无毛| 国产精品久久久久久久久免| 亚洲av成人精品一二三区| 一级爰片在线观看| 成人漫画全彩无遮挡| 网址你懂的国产日韩在线| 春色校园在线视频观看| 国产极品天堂在线| 亚洲图色成人| 亚洲欧洲日产国产| 久久久久国产网址| 亚洲四区av| 国产亚洲5aaaaa淫片| 亚洲人成网站在线播| 久久久久视频综合| 尾随美女入室| 王馨瑶露胸无遮挡在线观看| 久久热精品热| 欧美极品一区二区三区四区| 中文精品一卡2卡3卡4更新| 亚洲无线观看免费| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 美女高潮的动态| www.色视频.com| 久久这里有精品视频免费| 久久99热这里只频精品6学生| 国产av精品麻豆| a 毛片基地| 亚洲天堂av无毛| 最近最新中文字幕大全电影3| 有码 亚洲区| 成人高潮视频无遮挡免费网站| 我要看黄色一级片免费的| 久久久午夜欧美精品| 欧美成人午夜免费资源| 夫妻午夜视频| 黄色配什么色好看| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区| 久久久久久久久久久丰满| 夫妻午夜视频| 性高湖久久久久久久久免费观看| 下体分泌物呈黄色| 少妇被粗大猛烈的视频| 黄色欧美视频在线观看| 久久久久久久精品精品| 制服丝袜香蕉在线| 国产一级毛片在线| 各种免费的搞黄视频| 国产伦理片在线播放av一区| 丰满少妇做爰视频| 国产精品国产三级专区第一集| 亚洲精品乱码久久久v下载方式| 国产男女内射视频| 国产欧美亚洲国产| 麻豆成人午夜福利视频| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 久久ye,这里只有精品| 人人妻人人澡人人爽人人夜夜| 在线免费十八禁| 成人黄色视频免费在线看| 国产成人a区在线观看| 国产精品女同一区二区软件| 欧美日韩一区二区视频在线观看视频在线| 成年免费大片在线观看| 日韩大片免费观看网站| 国产精品不卡视频一区二区| 插阴视频在线观看视频| 黄片无遮挡物在线观看| 免费少妇av软件| 嫩草影院新地址| 高清黄色对白视频在线免费看 | 麻豆成人av视频| 99久久精品国产国产毛片| 欧美日韩在线观看h| 久久久a久久爽久久v久久| 超碰av人人做人人爽久久| 亚洲电影在线观看av| 亚洲av综合色区一区| 亚洲自偷自拍三级| 人妻制服诱惑在线中文字幕| 国产精品一区www在线观看| 日韩精品有码人妻一区| 免费人成在线观看视频色| 国产成人免费无遮挡视频| 国产爱豆传媒在线观看| 久久精品国产a三级三级三级| 伦理电影免费视频| 国产av码专区亚洲av| 国产亚洲一区二区精品| 日韩中文字幕视频在线看片 | 国模一区二区三区四区视频| 久久99精品国语久久久| 亚洲av成人精品一区久久| 久久99热这里只频精品6学生| 高清黄色对白视频在线免费看 | av天堂中文字幕网| 天堂俺去俺来也www色官网| 97超视频在线观看视频| 亚洲国产日韩一区二区| 国产成人aa在线观看| www.色视频.com| 国产成人一区二区在线| 在现免费观看毛片| 男女下面进入的视频免费午夜| 80岁老熟妇乱子伦牲交| 一级毛片我不卡| 少妇人妻一区二区三区视频| av在线app专区| 又粗又硬又长又爽又黄的视频| 国产在线男女| av一本久久久久| 日韩成人伦理影院| av国产精品久久久久影院| 精品熟女少妇av免费看| 在线亚洲精品国产二区图片欧美 | xxx大片免费视频| 免费看光身美女| 亚洲国产精品成人久久小说| 国产精品秋霞免费鲁丝片| 久久亚洲国产成人精品v| 亚洲av免费高清在线观看| 少妇被粗大猛烈的视频| 国产又色又爽无遮挡免| 春色校园在线视频观看| 久久精品国产自在天天线| 久久久欧美国产精品| 久久国内精品自在自线图片| 亚洲成人中文字幕在线播放| av视频免费观看在线观看| 亚洲精品国产av成人精品| 欧美日韩精品成人综合77777| 欧美变态另类bdsm刘玥| 18禁裸乳无遮挡免费网站照片| 少妇人妻一区二区三区视频| 99久久精品一区二区三区| 最后的刺客免费高清国语| 日本色播在线视频| 欧美变态另类bdsm刘玥| 高清欧美精品videossex| 日本午夜av视频| 久久久欧美国产精品| 大香蕉97超碰在线| 少妇人妻 视频| 国产伦在线观看视频一区| 久久久久精品性色| 高清日韩中文字幕在线| 亚洲最大成人中文| 秋霞在线观看毛片| 久久人人爽人人爽人人片va| 蜜桃在线观看..|