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

    基于MOD16的陜西省蒸散量時空分布特征

    2014-02-09 07:58:46范建忠李登科高茂盛
    生態(tài)環(huán)境學(xué)報 2014年9期
    關(guān)鍵詞:陜西省土地利用變化

    范建忠,李登科,高茂盛

    陜西省農(nóng)業(yè)遙感信息中心,陜西 西安 710014

    基于MOD16的陜西省蒸散量時空分布特征

    范建忠,李登科,高茂盛

    陜西省農(nóng)業(yè)遙感信息中心,陜西 西安 710014

    基于MOD16遙感數(shù)據(jù)集,在ERDAS IMAGINE 2013遙感圖像處理系統(tǒng)的支持下,通過空間建模,計算蒸散多年年平均值和月平均值,并生成圖像;結(jié)合陜西省矢量邊界圖、土地利用矢量圖,統(tǒng)計不同時間尺度統(tǒng)計行政區(qū)域和不同土地利用類型的蒸散值。在ARCGIS 10系統(tǒng)中,制作陜西省2000——2013年年、月平均蒸散分布圖。利用線性回歸進行蒸散時間趨勢分析,采用相關(guān)系數(shù)的統(tǒng)計檢驗方法進行顯著性趨勢檢驗。進而研究了陜西省2000——2013年蒸散量的空間分布特征和時間變化規(guī)律,分析了不同類型下蒸散量的差異性變化特征。結(jié)果表明:(1)全省年蒸散量在波動中緩慢上升,波動范圍為448.0~533.3 mm·a-1,年平均值493.3 mm·a-1。各月蒸散量的年際變化具有季節(jié)分異特征,秋末至仲春的月蒸散具有減少的趨勢,春末至仲秋的蒸散具有增加的趨勢。年內(nèi)蒸散量呈單峰型分布,季節(jié)性變化特征明顯,蒸散主要集中在5—9月份,最高、最低值分別出現(xiàn)在8月和11月。(2)多年平均蒸散空間格局呈現(xiàn)北低南高的分布規(guī)律,高植被覆蓋區(qū)蒸散量較大。蒸散變化趨勢不明顯的面積占77.2%,蒸散顯著、極顯著增加的像元主要分布在陜北地區(qū)、關(guān)中地區(qū)西部和陜南丘陵淺山區(qū),蒸散顯著和極顯著減少的像元主要分布在關(guān)中城市群。(3)土地利用特點影響著陜西省蒸散量的分布狀況,蒸散強度大小按類型排序依次為森林>草地>農(nóng)田>荒漠。研究結(jié)果對于陜西有限水資源的合理利用以及水資源短缺問題的解決、旱澇監(jiān)測和預(yù)警等研究具有重要意義。

    MOD16;蒸散量;時空分布;土地利用類型;陜西省

    全球暖化會影響大氣中的水汽含量和大氣環(huán)流,進而使降水、蒸散等水循環(huán)系統(tǒng)發(fā)生明顯變化(尹云鶴等,2010)。蒸散(Evapotranspiration,ET)是植被及地面向大氣輸送的水汽總通量,既包括從土壤和植物表面的水分蒸發(fā),也包括植被的蒸騰,它是土壤-植被-大氣系統(tǒng)中能量、水分交換的主要途徑,對ET過程的研究既是了解能量平衡和水分循環(huán)的重要環(huán)節(jié),也是深入認識陸面過程的基礎(chǔ)(王萬同等,2012;Liu等,2012)。水旱災(zāi)害是影響陜西農(nóng)業(yè)生產(chǎn)的主要災(zāi)種,旱澇災(zāi)害的發(fā)生意味著水量平衡關(guān)系發(fā)生了變化(Mahmood和Hubbard,2013)。由于蒸散量在水量平衡關(guān)系中占據(jù)重要比重,因此研究陜西省不同覆蓋下的蒸散及其空間分布特征,不僅可以加深對陸面過程的認識,而且對于陜西有限水資源的合理利用以及水資源短缺問題的解決、旱澇監(jiān)測和預(yù)警等研究具有重要意義。

    傳統(tǒng)計算方法都是以點的觀測資料為基礎(chǔ)開展蒸散研究的(YIN,2013;陽伏林等,2013;云文麗等,2013)。蒸散的時空分布與區(qū)域氣象條件、土壤水分狀況、植被覆蓋等因素緊密相關(guān)。較大尺度的陸面蒸散特征差異, 使得傳統(tǒng)的局地尺度研究方法很難直接應(yīng)用到區(qū)域尺度(郭曉寅和程國棟,2004;Batra等,2006)。而遙感技術(shù)具有快速、經(jīng)濟、宏觀等特點,尤其是它的可見光、近紅外和熱紅外波段數(shù)據(jù)可提供大范圍的特征參數(shù)和熱信息,這使蒸散研究從站點走向區(qū)域、從定性走向定量半定量成為現(xiàn)實(陳云浩等,2002;孫志剛等,2004;張曉濤等,2006;鄭有飛等,2011;何延波等,2007)。結(jié)合遙感信息和地面氣象要素的蒸散模型的發(fā)展,為區(qū)域蒸散的可靠估算提供了科學(xué)基礎(chǔ)(曾麗紅等,2010;劉朝順等,2010;付剛等,2010;莫興國等,2011)。

    目前,MODIS作為新一代資源衛(wèi)星傳感器,其數(shù)據(jù)和產(chǎn)品已經(jīng)越來越多地應(yīng)用于能量平衡的監(jiān)測過程中(杜嘉等,2010;徐永明等,2011;金曉媚等2013)。2011年美國NASA研究團隊在MODIS遙感數(shù)據(jù)蒸散反演算法上取得了重要成果,并通過NASA地球觀測系統(tǒng)發(fā)布了全球MODIS陸地蒸散產(chǎn)品數(shù)據(jù)(MOD16)(Mu等,2007,2011;吳桂平等2013),該產(chǎn)品不僅提供了蒸散量的特征參數(shù),還具有高時間分辨率以及免費獲取等特點,因此利用MOD16產(chǎn)品來反映陜西省蒸散量的時空分布及其與土地利用類型的關(guān)系具有一定的優(yōu)勢。

    近年來,一些學(xué)者對陜西的蒸散進行了研究(陳丹等,2006;楊文峰等,2013;張俊俊等,2013)。但這些研究僅涉及短時間、局部區(qū)域的模擬方法或變化特征研究,還未見對整個區(qū)域蒸散的分異特征進行系統(tǒng)研究。

    本文擬以1 km空間分辨率的MOD16產(chǎn)品數(shù)據(jù)為基礎(chǔ),在探討利用該產(chǎn)品進行陜西省蒸散量分析的可行性基礎(chǔ)之上,根據(jù)流域土地利用數(shù)據(jù)和地形數(shù)據(jù),通過GIS空間分析技術(shù),對流域不同土地利用類型下蒸散量進行統(tǒng)計分析,從而深入分析陜西省在年、月時間尺度下,蒸散的數(shù)量特征和空間格局特征,更進一步地探尋這種空間分布特征與不同土地利用類型之間的關(guān)系,以期揭示陜西省蒸散量的變化規(guī)律,為干旱災(zāi)害的成因分析和生態(tài)需水量研究提供科學(xué)依據(jù),同時為加強水資源管理與水分高效利用提供參考。

    1 研究地區(qū)和研究方法

    1.1 研究區(qū)概況

    陜西省(31°42′~39°35′ N,105°29′~111°15′ E)地處我國西北地區(qū),北部跨黃土高原中部,總面積2.1×105km2,平均海拔1127 m。地勢的總特點是南北高,中部低。北山和秦嶺把陜西分為三大自然區(qū)域,北部是陜北高原,中部是關(guān)中平原,南部是秦巴山區(qū),分別占全省總面積的45%、19%和36%。陜西縱貫三個氣候帶,南北氣候差異較大。陜南屬亞熱帶氣候,關(guān)中及陜北大部屬暖溫帶氣候,陜北北部長城沿線屬中溫帶氣候。年降水量310~1274 mm,降水南多北少,陜南為濕潤區(qū),關(guān)中為半濕潤區(qū),陜北為半干旱區(qū)。自然環(huán)境復(fù)雜,生態(tài)條件多樣,植物資源種類繁多,植被分布具有明顯的地帶性。從北到南依次分布溫帶草原、森林草原、暖溫帶落葉闊葉林和北亞熱帶常綠闊葉林,全省森林覆蓋率41.4%。自然災(zāi)害發(fā)生頻率較高,且區(qū)域性、季節(jié)性、伴生性特征突出。在自然災(zāi)害中,干旱最嚴重。

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

    蒸散產(chǎn)品(MOD16)下載地址為http://www.ntsg.umt.edu/project/mod16。根據(jù)MOD16產(chǎn)品數(shù)據(jù)軌道號的排列規(guī)律及陜西省所在地理位置,選擇的衛(wèi)星軌道號為h26v05和h27v05,其中涵蓋了陜西省2000—2013年共14年的數(shù)據(jù)。MOD16產(chǎn)品包括全球植被覆蓋區(qū)域的8天、月、年時間尺度的蒸散量(ET)、潛熱通量(LE)、潛在蒸散量(PET)、潛在潛熱通量(PLE),空間分辨率為1 km。

    陜西省土地利用數(shù)據(jù)來自地球系統(tǒng)科學(xué)共享網(wǎng)(http://www.geodata.cn/),該數(shù)據(jù)是由中國科學(xué)院地理科學(xué)與資源研究所、中國科學(xué)院遙感應(yīng)用研究所等8個單位共同完成的2005年全國32省、自治區(qū)1∶25萬的土地覆被矢量數(shù)據(jù),其投影方式為Albers正軸等面積雙標(biāo)準(zhǔn)緯線圓錐投影,坐標(biāo)系為Krasovsky坐標(biāo)系。該土地覆被遙感分類系統(tǒng)是針對中國土地覆被實際情況,在聯(lián)合國糧農(nóng)組織提出的LCCS(Land Cover Classification System)土地覆被分類體系和中國科學(xué)院環(huán)境數(shù)據(jù)庫土地利用分類體系的基礎(chǔ)上,從遙感制圖角度和陸地生態(tài)系統(tǒng)觀點出發(fā)而建立的一種新土地覆被遙感分類體系,包括6個一級分類和25個二級分類(白燕和廖順寶,2010)。

    1.3 數(shù)據(jù)處理

    原始的MODIS產(chǎn)品是采用分級數(shù)據(jù)格式(HDF,Hierarchical DataFormat)、正弦曲線投影存儲的,因此首先需要利用NASA提供的MRT軟件,將MOD16-ET產(chǎn)品的HDF文件轉(zhuǎn)換為WGS-1984經(jīng)緯度坐標(biāo)系統(tǒng)下的GeoTiff格式文件,并進行投影轉(zhuǎn)換、軌道鑲嵌和重采樣等操作。

    在ERDAS IMAGINE 2013遙感圖像處理系統(tǒng)的支持下,通過空間建模,計算蒸散多年年平均值和月平均值,并生成圖像;利用陜西省的矢量邊界圖、土地利用矢量圖,對不同時間尺度序列統(tǒng)計行政區(qū)域和各種土地利用類型的蒸散平均值,進而分析研究區(qū)域及其不同土地利用類型蒸散量的年際和月際變化特征。MOD16產(chǎn)品覆蓋范圍為有植被區(qū)域,不包括水體、城市等無植被覆蓋區(qū)域,因此在統(tǒng)計各土地利用類型蒸散量時不包括水體和城鎮(zhèn),進行行政區(qū)域內(nèi)蒸散量時不考慮水體、城市等無植被覆蓋區(qū)域的蒸散量。在ARCGIS 10系統(tǒng)中,制作陜西省2000—2013年年、月平均蒸散分布圖。

    利用線性傾向估計進行蒸散時間趨勢分析,采用相關(guān)系數(shù)的統(tǒng)計檢驗方法進行顯著性趨勢檢驗。隨時間變化,蒸散常表現(xiàn)為序列整體的上升或下降趨勢、空間分布格局變化以及在某時刻出現(xiàn)的轉(zhuǎn)折或突變。這些變量可以看作是時間的一元線性回歸,線性傾向值用最小二乘法估計:

    式中:B為線性傾向值,x為蒸散量,t為年份,n=14。當(dāng)B>0時,隨時間t的增加,x呈上升趨勢;當(dāng)B<0時,隨時間t的增加,x呈下降趨勢。B值大小反映了蒸散量上升或下降的速率,即表示上升或下降的傾向程度。如果回歸方程的相關(guān)系數(shù)通過信度為0.05、0.01的顯著性水平(P<0.05、P<0.01),則蒸散減小或增加趨勢分別達到顯著和極顯著水平。

    圖1 陜西省2000—2013年蒸散量的年際變化Fig. 1 Annual variation of evapotranspiration in Shaanxi Province during 2000─2013

    表1 陜西省蒸散量線性回歸方程Table 1 Equations of linear regress of evapotranspiration in Shaanxi Province

    2 結(jié)果與分析

    2.1 蒸散的時間變化

    2000—2013年陜西省年蒸散量波動范圍為448.0~533.3 mm·a-1,多年平均蒸散量為493.3 mm·a-1(圖1)。年蒸散量明顯超出多年平均值的年份有4年(2002、2003、2004、2012年),其中2003年蒸散量最高,超出多年平均值40.0 mm·a-1,相對變化率達到了8.1%;明顯低于多年平均值的年份有5年(2000、2001、2009、2010、2011年),2001年蒸散量最低,比多年平均值少45.3 mm·a-1,相對變化率達到了-9.2%。

    陜西省各月蒸散量的年際變化具有季節(jié)分異特征(表1)。12月至翌年4月的蒸散的線性傾向值為負值,說明秋末至仲春的蒸散年際變化具有減少的趨勢,1月和12月蒸散減少趨勢達到顯著水平(P<0.05);5月至10月的蒸散的線性傾向值為正值,說明春末至仲秋的蒸散年際變化具有增加的趨勢,8月蒸散增加趨勢達到顯著水平(P<0.05)。

    陜西省蒸散量年內(nèi)分布呈現(xiàn)先增大后減小的單峰型分布趨勢,蒸散量主要集中在5—9月份,其中11月蒸散發(fā)量最小,僅為23.0 mm,8月蒸散發(fā)量最大,達到75.8 mm(圖2)。究其原因主要在于:1、2、3、4、10、11、12月份陜西省的氣溫較低,不利于蒸散發(fā);4、5、6月份氣溫迅速回升,因此蒸散量也隨之升高;7、8月份氣溫達到最高值,并且在此期間降雨量大,供水充分,日照充足,提供了有利于蒸散發(fā)的充分條件,同時在該時間段內(nèi)8月份的蒸散量達到了最高值;9、10月份氣溫迅速降低,又向著不利于蒸散發(fā)的條件轉(zhuǎn)變,因此該階段的蒸散量又陡然回落。

    2.2 蒸散的空間變化

    陜西省多年蒸散平均值具有較強的空間分異性規(guī)律,呈現(xiàn)出北低南高的變化趨勢(圖3),該特征與覆蓋的地帶性變化大體一致。在植被覆蓋度高的地方,蒸散量較高。具體表現(xiàn)為,延安市以北蒸散在200~400 mm·a-1之間,延安市以南至秦嶺北坡蒸散大部在400~600 mm·a-1之間,秦巴山區(qū)蒸散大部在600~800 mm·a-1之間,蒸散>mm·a-1的像元零散地分布在秦巴山地。陜西省平均年蒸散量在200~800 mm之間的像元占全省94.9%,其中200~400、400~600、600~800 mm的像元分別占總像元的34.8%、28.1%和32.0%。

    受到太陽輻射、溫度、風(fēng)速等因素的影響,陜西省內(nèi)不同季節(jié)蒸散量的空間分布差異明顯(圖4)。1月份,陜西省蒸散量絕大部分區(qū)域在20~40 mm之間,全省分布比較均勻,全省平均為25.6 mm。4月份,陜西省蒸散量平均為28.0 mm;與1月份相比,4月份蒸散量陜北地區(qū)有了一定程度的降低,關(guān)中喝陜南有一定程度增加,其余地區(qū)變化較?。魂儽钡貐^(qū)大部在1~20 mm之間;陜北地區(qū)南部的子午嶺、黃龍山林區(qū),渭北地區(qū)以及秦嶺北坡在20~40 mm之間;關(guān)中地區(qū)大部在20~100 mm之間,局部地區(qū)〉100 mm;陜南的漢中、安康地區(qū)在40~60 mm之間,局部地區(qū)>60 mm。7月份,由于氣溫較高,降雨量大,光照充足,因此月平均蒸散量相對于4月份總體上得到了較大升高,不同地形地貌、不同植被的蒸散量差異較大,層次分明,全省月平均蒸散量達到74.1 mm;榆林市以北的長城沿線風(fēng)沙區(qū)蒸散量在1~20 mm之間,榆林市已南、延安市北部的黃土高原丘陵溝壑區(qū)蒸散量在20~40 mm之間,延安市中、南部以及關(guān)中地區(qū)北蒸散量40~100 mm(林區(qū)除外),子午嶺、黃龍山林區(qū),寶雞以北的關(guān)山林區(qū)以及秦巴山地蒸散量較高,大部分區(qū)域在100 mm以上。自9月份開始,受氣溫、土地覆蓋類型、以及風(fēng)速等因素的影響,陜西省月蒸散量又開始降低,到10月份全省月平均蒸散量僅為31.0 mm。全省大部月平均蒸散量在20~40 mm之間,長城沿線風(fēng)沙區(qū)在1~20 mm之間,秦巴山地大部在40~60 mm之間。

    圖2 陜西省2000—2013年平均蒸散月變化Fig. 2 The land surface ET of different months in Shaanxi Province during 2000—2013

    圖3 陜西省2000—2013年平均蒸散分布圖Fig.3 Averages of evapotranspiration during 2000—2013

    對于2000—2013年陜西省蒸散圖,計算每個像元蒸散與時間的相關(guān)系數(shù),并進行顯著性檢驗,制作成植被覆蓋度變化趨勢圖(圖5)。從圖5可見,蒸散顯著、極顯著增加的像元主要分布在陜北地區(qū)、關(guān)中地區(qū)西部和陜南丘陵淺山區(qū),這些地方是1998年以來陜西省實施退耕還林工程的重點區(qū)域,植被覆蓋有了顯著增加(李登科等,2010),蒸散也隨之顯著增加;蒸散顯著和極顯著減少的像元主要分布在關(guān)中城市群,說明隨著城鎮(zhèn)化的發(fā)展,建筑面積不斷擴大,植被覆蓋面積減少,蒸散也隨之顯著減少。研究期間,陜西省蒸散極顯著(P<0.01)、顯著(P<0.05)增加的面積分別占總面積的8.7%和7.7%,變化趨勢不明顯的面積占77.2%,減少極顯著(P<0.01)、顯著(P<0.05)的面積分別占3.0%和3.4%。

    圖4 陜西省1月、4月、7月、10月份蒸散量空間分布Fig.4 Spatial distribution of evapotranspiration in January, April, July and October over Shaanxi Province

    2.3 不同土地利用類型的蒸散量分布特征

    不同土地利用類型由于其本身的生理生態(tài)特性以及所處區(qū)域的降水差異,平均蒸散量及頻率分布特征均表現(xiàn)出不同的變化特點。為了進一步分析各種土地利用類型的蒸散特征,利用ERDAS IMAGINE的空間統(tǒng)計功能,分別提取了各種土地利用類型年均蒸散量統(tǒng)計特征(圖6)。從圖6可以發(fā)現(xiàn),陜西省年均蒸散量與土地利用類型密切相關(guān),不同土地利用類型的平均蒸散量存在著較大差異,森林的年蒸散量均值最高,達到576.0 mm·a-1,大大超過了其他土地利用類型的蒸散量,表明森林一般比較濕潤,土壤水分供應(yīng)充沛,因此年均蒸散量相對較高;荒漠的年ET均值最低,僅為224.3 mm·a-1。除聚落、水體外,各種土地覆蓋類型的年均蒸散量按照“林地>草地>農(nóng)田>荒漠”的順序遞減。

    根據(jù)MOD16月產(chǎn)品數(shù)據(jù),結(jié)合陜西省土地利用類型,可以計算得到各種不同土地覆蓋類型上的多年月平均蒸散量:森林為48.1 mm,草地為41.6 mm,農(nóng)田為36.4 mm,荒漠為18.6 mm。圖7顯示了2000—2013年陜西省各種土地利用類型每月平均蒸散量變化特征??傮w上來講,不同土地利用類型的月平均蒸散量大小順序與年蒸散量平均值比較一致,但是具體到不同土地利用類型上,其全年各月平均蒸散量的分布特征也不盡相同。森林、草地、農(nóng)田、荒漠最大、最小月蒸散量出現(xiàn)的月份不同,森林最大月蒸散量出現(xiàn)在7月,草地、農(nóng)田、荒漠則為8月;森林最小月蒸散量出現(xiàn)在12月,草地、農(nóng)田則為11月,而荒漠最小月蒸散量出現(xiàn)在4月。森林、草地、農(nóng)田蒸散量的季節(jié)變化呈現(xiàn)單峰型,而荒漠蒸散量的季節(jié)變化呈現(xiàn)一高一低的波浪型變化。

    圖5 陜西省蒸散量變化趨勢分布Fig.5 Spatial distribution of evapotranspiration change trend in Shaanxi Province

    圖6 不同土地利用類型的年平均蒸散量Fig. 6 Annual mean of evapotranspiration over different types of landuse

    3 結(jié)論

    本文基于MOD16蒸散產(chǎn)品,在ERDAS IMAGINE遙感圖像處理系統(tǒng)和ARCGIS地理信息系統(tǒng)的支持下,研究了陜西省蒸散時空分布特征,得到以下結(jié)論:

    1)2000—2013年,陜西省蒸散量年際變化呈現(xiàn)在波動中緩慢增加趨勢,波動范圍為448.0~533.3 mm·a-1,多年平均蒸散量為493.3 mm·a-1。各月蒸散量的年際變化具有季節(jié)分異特征,秋末至仲春的月蒸散具有減少的趨勢,春末至仲秋的蒸散具有增加的趨勢。陜西省蒸散量年內(nèi)分布呈現(xiàn)先增大后減小的單峰型分布趨勢,蒸散量主要集中在5—9月份,其中11月蒸散發(fā)量最小,8月蒸散發(fā)量最大。

    2)陜西省多年蒸散平均值具有較強的空間分異性規(guī)律,呈現(xiàn)出北低南高的變化趨勢,即陜北北部年蒸散在200~400 mm之間,陜北南部和關(guān)中地區(qū)蒸散在400~600 mm之間,秦巴山區(qū)蒸散大部在600~800 mm·a-1之間,局部蒸散>800 mm·a-1。該特征與覆蓋的地帶性變化大體一致。

    3)14年來,陜西省蒸散變化趨勢不明顯的面積占77.2%,蒸散極顯著(P<0.01)、顯著(P<0.05)增加的面積占總面積的16.4%,極顯著(P<0.01)、顯著(P<0.05)減少的面積占6.4%。蒸散顯著、極顯著增加的像元主要分布在陜北地區(qū)、關(guān)中地區(qū)西部和陜南丘陵淺山區(qū),蒸散顯著、極顯著減少的像元主要分布在關(guān)中城市群。

    4)陜西省不同土地利用類型平均蒸散量表現(xiàn)出不同的變化特點。森林的年蒸散均值最高,草地、農(nóng)田次之,荒漠最低;各土地利用類型蒸散量的季節(jié)變化大小順序與年尺度的基本一致,但是由于氣象條件、農(nóng)作物特性以及人類活動等綜合作用的影響,不同土地利用類型蒸散量的分布特征不盡相同。

    圖7 不同土地利用類型的月蒸散量分布Fig. 7 monthly mean of evapotranspiration over landuse types

    BATRA N, ISLAM S, VENTURINI V et al. 2006. Estimation and comparison of evapotranspiration from MODIS and AVHRR sensors for clear sky days over the Southern Great Plains[J]. Remote Sensing of Environment, 103: 1-15.

    LIU J, CHEN J M, CHILAR J. 2003. Mapping evapotranspiration based on remote sensing: An application to Canada's landmass. Water Resources Research, 280:1189-1200.

    Mahmood R, Hubbard K G. 2003. Simulating sensitivity of soil moisture and evapotranspiration under heterogeneous soils and land uses[J]. Journal of Hydrology, 280: 72-90.

    MU Q Z, HEINSCH F A, ZHAO M S, et al. 2007.Development of a global evapotranspiration algorithm based on MODIS and global meteorology data[J]. Remote sensing of Environment, 111: 519-536.

    MU Q Z, ZHAO M S, RUNNING S W. 2011. Improvements to a MODIS global terrestrial evapotranspiration algorithm[J]. Remote sensing of Environment, 115: 1781-1800.

    YIN Y H, WU S H, Z D S, et al. 2013. Modeled effects of climate change on actual evapotranspiration on the different eco-geographical regions in the Tibetan Plateau[J]. Journal of Geographical Sciences, 23(2): 195-207.

    白燕, 廖順寶. 2010. 矢量數(shù)據(jù)屬性信息無損柵格化的實現(xiàn)方法——以全國1∶ 25萬土地覆被數(shù)據(jù)為例[J]. 地球信息科學(xué)學(xué)報, 12(3): 385-391.

    曾麗紅, 宋開山, 張柏, 等. 2010.松嫩平原不同覆蓋蒸散特征的遙感研究[J].農(nóng)業(yè)工程學(xué)報, 26(9): 233-242.

    陳丹, 莫興國, 林忠輝, 等. 2006. 基于MODIS 數(shù)據(jù)的無定河流域蒸散模擬[J]. 地理研究, 25(4): 617-623.

    陳云浩, 李曉兵, 史培軍. 2002. 非均勻陸面條件下區(qū)域蒸散量計算的遙感模型[J]. 氣象學(xué)報, 60(4): 508-512.

    杜嘉, 張柏, 宋開山, 等. 2010. 基于MODIS產(chǎn)品和SEBAL模型的三江平原日蒸散量估算[J]. 中國農(nóng)業(yè)氣象, 31(1): 104-110.

    付剛, 沈振西, 張憲洲, 等. 2010. 基于MODIS影像的藏北高寒草甸的蒸散模擬[J]. 草業(yè)學(xué)報, 19(5): 103-112.

    郭曉寅, 程國棟. 2004. 遙感技術(shù)應(yīng)用于面蒸散發(fā)的研究進展[J]. 地球科學(xué)進展, 19(1): 107-114.

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

    金曉媚, 郭任宏, 夏薇. 2013. 基于MODIS數(shù)據(jù)的柴達木盆地區(qū)域蒸散量的變化特征[J]. 水文地質(zhì)工程地質(zhì), 40(11): 8-13.

    李登科, 范建忠, 王娟. 2010.陜西省植被覆蓋度變化特征及其成因[J].應(yīng)用生態(tài)學(xué)報, 21(11): 2896-2903.

    劉朝順, 施潤和, 高煒, 等. 2010.利用區(qū)域遙感ET分析山東省水分盈虧的研究[J]. 自然資源學(xué)報, 25(11): 1938-1948.

    莫興國, 劉蘇峽, 林忠輝, 等. 2011. 華北平原蒸散和GPP格局及其對氣候波動的響應(yīng)[J]. 地理學(xué)報, 66(5): 589-598.

    孫志剛, 王勤學(xué), 歐陽竹, 等. 2004.MODIS水汽通量估算方法在華北平原農(nóng)田的適應(yīng)性驗證[J]. 地理學(xué)報, 59(1): 49-55.

    王萬同, 趙慶良, 杜佳. 2012. 基于遙感技術(shù)的區(qū)域蒸散研究進展[J]. 國土資源遙感, (1): 1-7.

    吳桂平, 劉元波, 趙曉松, 等. 2013. 基于MOD16產(chǎn)品的鄱陽湖流域蒸散量時空分布特征[J]. 地理研究, 32(4): 617-627.

    徐永明, 趙巧華, 巴雅爾, 等. 2011. 基于MODIS數(shù)據(jù)的博斯騰湖流域蒸散時空變化[J]. 地理科學(xué), 32(11): 1353-1357

    陽伏林, 張強, 王潤元, 等. 2013.黃土高原半干旱區(qū)農(nóng)田生態(tài)系統(tǒng)蒸散與作物系數(shù)特征[J]. 應(yīng)用生態(tài)學(xué)報, 24(5): 1209-1214.

    楊文峰, 李星敏, 盧玲. 2013. 基于能量平衡的蒸散遙感估算模型的應(yīng)用研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 41(2): 46-52.

    尹云鶴, 吳紹洪, 戴爾阜. 2010. 1971─2008年我國潛在蒸散時空演變的歸因[J]. 科學(xué)通報, 55(22): 2226-2234.

    云文麗, 侯瓊, 王永利, 等. 2013.內(nèi)蒙古典型草原作物系數(shù)和實際蒸散量的時空分布特征[J]. 自然資源學(xué)報, 28(2): 300-311.

    張俊俊, 劉文兆, 韓曉陽. 2013. 陜西長武塬區(qū)冬小麥的蒸散特征[J]. 江蘇農(nóng)業(yè)科學(xué), 41(12): 60-63.

    張曉濤, 康紹忠, 王鵬新, 等. 2006. 估算區(qū)域蒸發(fā)蒸騰量的遙感模型對比分析[J]. 農(nóng)業(yè)工程學(xué)報, 22(7):6-3.

    鄭有飛, 陳鵬, 吳榮軍, 等. 2011. 蒸散的遙感估算模型及其在農(nóng)業(yè)干旱監(jiān)測中的應(yīng)用[J]. 生態(tài)學(xué)雜志, 30(4): 837-844.

    Spatio-temporal Variations of Evapotranspiration in Shaanxi Province Using MOD16 Products

    FAN Jianzhong, LI Dengke, GAO Maosheng
    Shaanxi Remote Sensing Information Center for Agricultuer, Xi’an 710014, China

    Based on MOD16 dataset, annual and monthly average evapotranspiration were calculated by the tool of spatial model maker in ERDAS IMAGINE 2013, evapotranspiration of the administrative regions and landuse types for different time scales was statistically computed with the vector maps of Shaanxi Province boundary and landuse. Annual and monthly average evapotranspiration images were drew by ARCGIS 10. The time trends of evapotranspiration were analyzed by using linear regression and the significance was tested by correlation coefficient statistical test. Then the spatio-temporal variation characteristics of evapotranspiration from 2000 to 2013 in Shaanxi Province and the differences for evapotranspiration under different landuse types were analyzed. The results show that the mean annual evapotranspiration of Shaanxi Province increased slowly in the fluctuations with values ranging from 448.0 mm·a-1to 533.3 mm·a-1, and with a mean value of 493.3 mm·a-1. The annual variation of the mean monthly evapotranspiration presented seasonal difference, the mean evapotranspiration in the period from December to April of next year showed decreasing trends, the mean evapotranspiration in the period from May to October showed increasing trends. The monthly variation of evapotranspiration within the year shows a unimodal pattern, and its peak value occurs in August, its bottom value occurs in November, it is mainly concentrated in the period from May to September. The spatial distribution of mean annual evapotranspiration from 2000 to 2013 was generally high over the south of the Shaanxi Province, but low in the north, especially in the high vegetation coverage areas. The areas with insignificant variation trends were accounted for 77.2% of the total, the pixels with significantly increased evapotranspiration were mainly distributed in north Shaanxi, Mid-Shaanxi, and hilly-mountainous areas in south Shaanxi, the pixels with significantly decreased evapotranspiration were mainly distributed in Mid-Shaanxi urban agglomeration. The spatial distribution of evapotranspiration is related to land use types, which have effect on the evapotranspiration distribution in the Shaanxi Province, and the order of evapotranspiration intensity by land cover types is woodland > grassland >farmland >desert. The research results have important significance to utilization of the limited water resources, to resolution of the shortage of water resources, and to research of monitor and early warning for drought and flood disaster.

    MOD16; evapotranspiration; spatio-temporal variation; landuse type; Shaanxi province

    P426.2

    A

    1674-5906(2014)09-1536-08

    范建忠,李登科,高茂盛. 基于MOD16的陜西省蒸散量時空分布特征[J]. 生態(tài)環(huán)境學(xué)報, 2014, 23(9): 1536-1543.

    FAN Jianzhong, LI Dengke, GAO Maosheng. Spatio-temporal Variations of Evapotranspiration in Shaanxi Province Using MOD16 Products [J]. Ecology and Environmental Sciences, 2014, 23(9): 1536-1543.

    陜西省科學(xué)研究發(fā)展計劃項目(2013K02-20)

    范建忠(1961年生),男,高級工程師,從事氣候和遙感應(yīng)用研究。

    *通信作者:李登科,男,正研級高級工程師,從事農(nóng)業(yè)氣象和遙感應(yīng)用研究。E-mail:ldk81711@sohu.com

    2014-06-20

    猜你喜歡
    陜西省土地利用變化
    陜西省自強中等專業(yè)學(xué)校簡介
    從9到3的變化
    陜西省抓黨建促脫貧攻堅的實踐與思考
    聚焦兩會
    陜西畫報(2018年1期)2018-11-17 19:33:14
    這五年的變化
    特別文摘(2018年3期)2018-08-08 11:19:42
    土地利用生態(tài)系統(tǒng)服務(wù)研究進展及啟示
    陜西省閱讀文化節(jié)
    濱??h土地利用挖潛方向在哪里
    鳥的變化系列
    詩選刊(2015年6期)2015-10-26 09:47:11
    “讓路”讓出的變化
    中國火炬(2012年2期)2012-07-24 14:17:54
    av网站在线播放免费| svipshipincom国产片| √禁漫天堂资源中文www| 亚洲欧美日韩高清在线视频| 久久国产精品影院| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av在线 | 51午夜福利影视在线观看| 嫁个100分男人电影在线观看| 国产精华一区二区三区| 国产男女内射视频| 亚洲av片天天在线观看| 免费在线观看影片大全网站| 精品国产一区二区三区四区第35| 国产精品av久久久久免费| 在线看a的网站| 欧美日韩视频精品一区| 久久久国产一区二区| av国产精品久久久久影院| 老熟妇乱子伦视频在线观看| 国产亚洲精品第一综合不卡| a级片在线免费高清观看视频| 俄罗斯特黄特色一大片| 久久久久久久精品吃奶| 国产精品 欧美亚洲| 两个人看的免费小视频| 黑人欧美特级aaaaaa片| 久久影院123| 免费人成视频x8x8入口观看| 91成年电影在线观看| 国产一区二区三区综合在线观看| 男女床上黄色一级片免费看| 宅男免费午夜| 午夜福利一区二区在线看| 国产成人系列免费观看| 男女午夜视频在线观看| 老熟女久久久| 丝瓜视频免费看黄片| 久久久久久免费高清国产稀缺| 999久久久精品免费观看国产| 国产成人欧美| 亚洲一区中文字幕在线| 国产视频一区二区在线看| 欧美久久黑人一区二区| 校园春色视频在线观看| 一进一出抽搐动态| 日本wwww免费看| 别揉我奶头~嗯~啊~动态视频| 最新的欧美精品一区二区| 精品国产乱码久久久久久男人| 欧美日韩中文字幕国产精品一区二区三区 | 精品少妇久久久久久888优播| 正在播放国产对白刺激| 欧美精品av麻豆av| 国产aⅴ精品一区二区三区波| 国产在线精品亚洲第一网站| 国产亚洲一区二区精品| 国产精品免费一区二区三区在线 | 亚洲avbb在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久人人97超碰香蕉20202| 人人妻人人爽人人添夜夜欢视频| 久久性视频一级片| 欧洲精品卡2卡3卡4卡5卡区| 9热在线视频观看99| 涩涩av久久男人的天堂| 黄色视频,在线免费观看| 后天国语完整版免费观看| www.精华液| av电影中文网址| 欧美日韩av久久| 十八禁高潮呻吟视频| 成年人免费黄色播放视频| 大香蕉久久成人网| 最近最新中文字幕大全电影3 | 国产精品1区2区在线观看. | 狠狠狠狠99中文字幕| 亚洲欧美激情综合另类| 搡老熟女国产l中国老女人| 老熟妇仑乱视频hdxx| 天天影视国产精品| 一夜夜www| 很黄的视频免费| 国产人伦9x9x在线观看| 一级a爱视频在线免费观看| 国产精品久久久av美女十八| 日韩视频一区二区在线观看| 99香蕉大伊视频| 人人澡人人妻人| 久久香蕉精品热| 亚洲国产精品sss在线观看 | 老熟妇仑乱视频hdxx| 别揉我奶头~嗯~啊~动态视频| 在线观看免费视频日本深夜| 久久国产精品影院| 亚洲熟女毛片儿| 一本综合久久免费| 国产精品98久久久久久宅男小说| 一级黄色大片毛片| 国产亚洲一区二区精品| 久久婷婷成人综合色麻豆| 国产精品影院久久| 亚洲第一av免费看| 99re在线观看精品视频| 欧美老熟妇乱子伦牲交| 一二三四社区在线视频社区8| 成年人免费黄色播放视频| 在线免费观看的www视频| 在线视频色国产色| 18禁观看日本| av天堂在线播放| 91成人精品电影| 免费看a级黄色片| 国产精品影院久久| 国产成人系列免费观看| 精品人妻1区二区| 亚洲国产毛片av蜜桃av| 美女午夜性视频免费| 国产欧美日韩一区二区三| 高清视频免费观看一区二区| 国产精品九九99| 91麻豆精品激情在线观看国产 | 久久久久精品人妻al黑| 欧美在线黄色| 精品一区二区三卡| 亚洲第一青青草原| 国产亚洲欧美在线一区二区| 日本a在线网址| 亚洲精品自拍成人| 在线观看免费午夜福利视频| 久久精品国产99精品国产亚洲性色 | 亚洲欧美精品综合一区二区三区| 中文字幕制服av| 婷婷成人精品国产| 天天躁夜夜躁狠狠躁躁| 在线天堂中文资源库| 亚洲第一欧美日韩一区二区三区| 在线观看免费高清a一片| 91字幕亚洲| 亚洲熟妇熟女久久| e午夜精品久久久久久久| 成人精品一区二区免费| 精品亚洲成国产av| 中文字幕人妻丝袜制服| 国产97色在线日韩免费| 欧美精品亚洲一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 大型av网站在线播放| 亚洲国产中文字幕在线视频| 亚洲专区国产一区二区| 80岁老熟妇乱子伦牲交| 高清毛片免费观看视频网站 | 欧美成狂野欧美在线观看| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 在线永久观看黄色视频| a级毛片在线看网站| 国产精品 欧美亚洲| 国产精品99久久99久久久不卡| 国产高清激情床上av| 久久亚洲真实| 国产高清videossex| av免费在线观看网站| 热99久久久久精品小说推荐| 国产欧美日韩一区二区精品| 一级毛片精品| 午夜免费观看网址| 亚洲第一欧美日韩一区二区三区| 午夜视频精品福利| 国产精品久久久久久人妻精品电影| 两个人看的免费小视频| 在线观看一区二区三区激情| 国产欧美日韩一区二区精品| 欧美 亚洲 国产 日韩一| 日韩欧美三级三区| 91麻豆av在线| 日韩免费av在线播放| 午夜亚洲福利在线播放| 久久影院123| 动漫黄色视频在线观看| 欧美黄色片欧美黄色片| 黄色片一级片一级黄色片| 免费一级毛片在线播放高清视频 | 亚洲国产中文字幕在线视频| 亚洲五月色婷婷综合| 国产免费av片在线观看野外av| 国产精品久久视频播放| 精品一区二区三区av网在线观看| 不卡av一区二区三区| 一区二区三区国产精品乱码| 欧美激情高清一区二区三区| 一夜夜www| av超薄肉色丝袜交足视频| 午夜福利在线观看吧| 69精品国产乱码久久久| 两人在一起打扑克的视频| 久久久久久久久免费视频了| 啦啦啦在线免费观看视频4| 久久精品国产综合久久久| 久久香蕉精品热| 亚洲伊人色综图| 精品亚洲成a人片在线观看| 国产男靠女视频免费网站| 亚洲精品成人av观看孕妇| 成人av一区二区三区在线看| 日韩成人在线观看一区二区三区| 另类亚洲欧美激情| 丰满的人妻完整版| 日韩欧美一区视频在线观看| 国产又色又爽无遮挡免费看| 国产精品久久久av美女十八| www.999成人在线观看| 欧美激情高清一区二区三区| 99久久国产精品久久久| 少妇的丰满在线观看| 免费日韩欧美在线观看| 国产91精品成人一区二区三区| 一边摸一边抽搐一进一小说 | 国产高清视频在线播放一区| 91成年电影在线观看| 日本vs欧美在线观看视频| 国产精品一区二区免费欧美| 极品少妇高潮喷水抽搐| 久久精品国产亚洲av高清一级| 色婷婷av一区二区三区视频| 欧美久久黑人一区二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲精华国产精华精| 后天国语完整版免费观看| 大片电影免费在线观看免费| 99国产精品99久久久久| 成人特级黄色片久久久久久久| 免费黄频网站在线观看国产| 国产欧美日韩精品亚洲av| 亚洲av第一区精品v没综合| 美女 人体艺术 gogo| 黑人欧美特级aaaaaa片| 99久久人妻综合| 男女免费视频国产| 欧美久久黑人一区二区| av天堂在线播放| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 国产亚洲精品第一综合不卡| 欧美精品一区二区免费开放| 欧美乱色亚洲激情| 日韩欧美三级三区| 女人精品久久久久毛片| 亚洲三区欧美一区| 岛国毛片在线播放| 老司机福利观看| 麻豆成人av在线观看| 在线天堂中文资源库| 欧美黄色淫秽网站| 美女国产高潮福利片在线看| 成年人免费黄色播放视频| 免费看十八禁软件| 精品国内亚洲2022精品成人 | 亚洲精品在线美女| 一二三四社区在线视频社区8| 女人精品久久久久毛片| 欧美成人午夜精品| 成人免费观看视频高清| 自拍欧美九色日韩亚洲蝌蚪91| 侵犯人妻中文字幕一二三四区| 欧美老熟妇乱子伦牲交| 亚洲一区高清亚洲精品| 一区二区三区激情视频| 免费高清在线观看日韩| www.精华液| 国产一区二区三区视频了| 露出奶头的视频| 国产精品 欧美亚洲| 又黄又爽又免费观看的视频| 亚洲国产毛片av蜜桃av| 久久午夜亚洲精品久久| 91九色精品人成在线观看| 久久亚洲精品不卡| 桃红色精品国产亚洲av| 免费日韩欧美在线观看| 亚洲av日韩精品久久久久久密| 亚洲avbb在线观看| 麻豆成人av在线观看| 大陆偷拍与自拍| 黄色a级毛片大全视频| 精品人妻在线不人妻| 亚洲成国产人片在线观看| 国产精品综合久久久久久久免费 | 80岁老熟妇乱子伦牲交| 午夜久久久在线观看| 精品久久久久久,| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区| 九色亚洲精品在线播放| 亚洲av欧美aⅴ国产| 日韩免费高清中文字幕av| 亚洲中文av在线| 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 亚洲免费av在线视频| 日韩免费av在线播放| 成人手机av| 欧美精品亚洲一区二区| 亚洲国产欧美一区二区综合| 国产精品99久久99久久久不卡| 黄频高清免费视频| 超碰成人久久| 亚洲欧美一区二区三区黑人| 在线十欧美十亚洲十日本专区| 人人妻人人澡人人爽人人夜夜| 午夜免费鲁丝| 国产精品二区激情视频| 亚洲欧美日韩高清在线视频| 亚洲av美国av| 午夜免费鲁丝| 天堂√8在线中文| 国产人伦9x9x在线观看| 1024香蕉在线观看| 正在播放国产对白刺激| 久久中文看片网| 国产片内射在线| 久久狼人影院| 免费久久久久久久精品成人欧美视频| 成年动漫av网址| 99久久综合精品五月天人人| 亚洲国产精品sss在线观看 | 亚洲在线自拍视频| 亚洲久久久国产精品| 国内久久婷婷六月综合欲色啪| av在线播放免费不卡| 国产精品亚洲av一区麻豆| 十八禁高潮呻吟视频| 国产精华一区二区三区| 不卡av一区二区三区| 久久久久久久午夜电影 | 久久国产亚洲av麻豆专区| 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美三级三区| 午夜免费鲁丝| 亚洲精品国产精品久久久不卡| 久久久久久久久免费视频了| www.999成人在线观看| 黑人欧美特级aaaaaa片| 国产欧美日韩精品亚洲av| 亚洲在线自拍视频| 欧美日本中文国产一区发布| av国产精品久久久久影院| 91精品三级在线观看| 国产成人av激情在线播放| 在线免费观看的www视频| 久久午夜综合久久蜜桃| 亚洲久久久国产精品| 久久久久久久久免费视频了| 99在线人妻在线中文字幕 | 免费女性裸体啪啪无遮挡网站| 午夜91福利影院| 黄色视频,在线免费观看| netflix在线观看网站| 91成年电影在线观看| 男男h啪啪无遮挡| 黄色视频不卡| 精品高清国产在线一区| 丝袜美腿诱惑在线| www日本在线高清视频| 高清黄色对白视频在线免费看| 中文字幕精品免费在线观看视频| 国产成人精品久久二区二区91| 老司机深夜福利视频在线观看| 欧美日韩乱码在线| 欧美性长视频在线观看| 亚洲一码二码三码区别大吗| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利在线观看吧| 在线观看66精品国产| 一级片免费观看大全| 欧美不卡视频在线免费观看 | 亚洲综合色网址| 久久狼人影院| 成人18禁在线播放| 日韩欧美一区视频在线观看| 中文字幕最新亚洲高清| 午夜影院日韩av| 午夜精品久久久久久毛片777| 99久久99久久久精品蜜桃| 丁香欧美五月| a在线观看视频网站| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看| 亚洲成国产人片在线观看| 欧美日韩黄片免| 大型av网站在线播放| 亚洲成人手机| 日韩欧美三级三区| 亚洲第一欧美日韩一区二区三区| 国产精品综合久久久久久久免费 | 水蜜桃什么品种好| 精品国产乱码久久久久久男人| 日日夜夜操网爽| 国产高清激情床上av| www.熟女人妻精品国产| 人妻一区二区av| 脱女人内裤的视频| 国产精品自产拍在线观看55亚洲 | 十八禁人妻一区二区| 动漫黄色视频在线观看| 视频区图区小说| 亚洲少妇的诱惑av| av网站在线播放免费| 久久精品国产亚洲av香蕉五月 | 久久久久久久午夜电影 | 免费黄频网站在线观看国产| 国产激情久久老熟女| 黄色女人牲交| 最近最新中文字幕大全电影3 | 成人18禁高潮啪啪吃奶动态图| 黄片小视频在线播放| 人妻丰满熟妇av一区二区三区 | 最近最新中文字幕大全电影3 | 国产日韩一区二区三区精品不卡| 久久久久精品人妻al黑| 久久精品国产综合久久久| 午夜亚洲福利在线播放| 岛国在线观看网站| 中国美女看黄片| 中文亚洲av片在线观看爽 | 免费看a级黄色片| 可以免费在线观看a视频的电影网站| 免费少妇av软件| 777米奇影视久久| 亚洲第一av免费看| 精品午夜福利视频在线观看一区| 国产主播在线观看一区二区| 亚洲在线自拍视频| 老司机影院毛片| 国产在视频线精品| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 别揉我奶头~嗯~啊~动态视频| 老熟女久久久| 高清av免费在线| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 国内毛片毛片毛片毛片毛片| 一级a爱视频在线免费观看| 两性夫妻黄色片| 韩国精品一区二区三区| 色综合欧美亚洲国产小说| 欧美激情极品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 一边摸一边做爽爽视频免费| bbb黄色大片| 亚洲国产毛片av蜜桃av| av天堂久久9| 国产一区二区激情短视频| 日韩免费高清中文字幕av| 久久亚洲真实| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 久热这里只有精品99| 久久热在线av| 黑丝袜美女国产一区| 激情在线观看视频在线高清 | av天堂在线播放| 国产亚洲欧美精品永久| 99热只有精品国产| 色播在线永久视频| cao死你这个sao货| 精品国产亚洲在线| 99久久精品国产亚洲精品| 99re在线观看精品视频| 亚洲在线自拍视频| 亚洲九九香蕉| 国产亚洲av高清不卡| 久99久视频精品免费| 在线视频色国产色| 久久人妻av系列| 91成年电影在线观看| a级毛片黄视频| 超色免费av| 亚洲av成人不卡在线观看播放网| 午夜久久久在线观看| 亚洲熟妇中文字幕五十中出 | 国产一区有黄有色的免费视频| 亚洲精品国产色婷婷电影| 国产午夜精品久久久久久| 老司机深夜福利视频在线观看| 18禁裸乳无遮挡免费网站照片 | 久久性视频一级片| 男女午夜视频在线观看| 中文字幕最新亚洲高清| 精品久久久久久久毛片微露脸| www日本在线高清视频| 十八禁人妻一区二区| 黑人巨大精品欧美一区二区蜜桃| www.精华液| 成年人免费黄色播放视频| 国产免费男女视频| 99riav亚洲国产免费| 大型黄色视频在线免费观看| 亚洲国产中文字幕在线视频| 国产精品九九99| 久久久久国产精品人妻aⅴ院 | 国产亚洲精品第一综合不卡| 欧美大码av| 欧美日韩中文字幕国产精品一区二区三区 | 天天操日日干夜夜撸| 精品一区二区三区视频在线观看免费 | 久久精品91无色码中文字幕| 久久久久久人人人人人| 欧美另类亚洲清纯唯美| 亚洲精品国产精品久久久不卡| 免费日韩欧美在线观看| 麻豆国产av国片精品| 国产在线观看jvid| 国产成人av激情在线播放| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲一区二区精品| 国产精品亚洲av一区麻豆| 黄片小视频在线播放| 搡老熟女国产l中国老女人| 18在线观看网站| 国产在线一区二区三区精| 成人18禁在线播放| 亚洲人成77777在线视频| 一级a爱视频在线免费观看| 国产成人av教育| 精品人妻在线不人妻| 成人精品一区二区免费| 久久精品91无色码中文字幕| 99国产极品粉嫩在线观看| 成人三级做爰电影| 欧美激情极品国产一区二区三区| 亚洲色图av天堂| bbb黄色大片| 欧美黑人精品巨大| 啦啦啦视频在线资源免费观看| 天天操日日干夜夜撸| 精品电影一区二区在线| 国产精品欧美亚洲77777| 久久久久久久午夜电影 | 啦啦啦免费观看视频1| 怎么达到女性高潮| 天天操日日干夜夜撸| 亚洲 国产 在线| 99久久综合精品五月天人人| 亚洲熟妇熟女久久| 精品福利观看| x7x7x7水蜜桃| 宅男免费午夜| 在线观看舔阴道视频| 十八禁人妻一区二区| 久久久国产欧美日韩av| 精品少妇久久久久久888优播| 久久99一区二区三区| 中文字幕制服av| svipshipincom国产片| 免费高清在线观看日韩| 国产亚洲一区二区精品| 精品午夜福利视频在线观看一区| 91精品三级在线观看| a级毛片黄视频| 热99re8久久精品国产| 午夜免费观看网址| 老司机深夜福利视频在线观看| 一区在线观看完整版| 亚洲国产毛片av蜜桃av| 夜夜躁狠狠躁天天躁| 狠狠狠狠99中文字幕| 中文亚洲av片在线观看爽 | 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 飞空精品影院首页| bbb黄色大片| 日韩成人在线观看一区二区三区| 色综合欧美亚洲国产小说| 一区二区日韩欧美中文字幕| 久久午夜综合久久蜜桃| 国产欧美亚洲国产| 日韩免费高清中文字幕av| 欧美精品人与动牲交sv欧美| 欧美成人午夜精品| 操出白浆在线播放| 午夜福利视频在线观看免费| 日韩欧美在线二视频 | 欧美久久黑人一区二区| 国产亚洲欧美在线一区二区| 老司机深夜福利视频在线观看| 亚洲一区二区三区欧美精品| 成人av一区二区三区在线看| 国产区一区二久久| 亚洲精品久久成人aⅴ小说| 欧美大码av| 大香蕉久久网| 高清欧美精品videossex| 麻豆av在线久日| 身体一侧抽搐| 免费不卡黄色视频| 精品一品国产午夜福利视频| 男男h啪啪无遮挡| 成人手机av| 亚洲精品久久成人aⅴ小说| 欧美日韩国产mv在线观看视频| tocl精华| 在线永久观看黄色视频| 国产精品久久久久成人av| 国产国语露脸激情在线看| 在线免费观看的www视频| 一边摸一边抽搐一进一出视频| 美女 人体艺术 gogo| 日本精品一区二区三区蜜桃| 色老头精品视频在线观看| 操美女的视频在线观看|