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

    基于SEBS模型估算遼西北地區(qū)蒸散發(fā)及時空特征

    2017-12-18 01:13:25張繼權(quán)曹永強郭恩亮李丹君馬齊云
    水土保持研究 2017年6期
    關(guān)鍵詞:西北地區(qū)通量反演

    王 蕊, 張繼權(quán), 曹永強, 郭恩亮, 李丹君, 馬齊云

    (1.東北師范大學(xué) 環(huán)境學(xué)院自然災(zāi)害研究所, 長春 130117; 2.遼寧師范大學(xué) 城市與環(huán)境學(xué)院, 遼寧 大連 116029)

    基于SEBS模型估算遼西北地區(qū)蒸散發(fā)及時空特征

    王 蕊1, 張繼權(quán)1, 曹永強2, 郭恩亮1, 李丹君1, 馬齊云1

    (1.東北師范大學(xué) 環(huán)境學(xué)院自然災(zāi)害研究所, 長春 130117; 2.遼寧師范大學(xué) 城市與環(huán)境學(xué)院, 遼寧 大連 116029)

    蒸散發(fā)是陸面過程中的重要環(huán)節(jié),聯(lián)系著陸面水循環(huán)和地表能量平衡。采用MODIS數(shù)據(jù),逐月選取2013年12期遙感數(shù)據(jù),并分季度求其平均值作為四季代表,結(jié)合SEBS模型,對遼西北地區(qū)地表蒸發(fā)蒸騰量進行了估算,并為了驗證遙感反演精度,提取不同土地利用類型的蒸散發(fā)量。同時,利用Penman-Monteith公式計算遼西北地區(qū)近50年潛在蒸散發(fā)值,通過一元線性趨勢分析法、小波分析法、IDW反距離加權(quán)法和R/S分析法分析其時空分布特征。結(jié)果表明:遼西北地區(qū)2013年蒸散發(fā)量夏季>春季>秋季>冬季;近50年來潛在蒸散發(fā)值由高到低依次為夏季、春季、秋季和冬季,且具有明顯的區(qū)域性。同時,經(jīng)驗證,遙感手段反演蒸散發(fā)值誤差約為7.5%,表明該方法在研究區(qū)應(yīng)用合理可行。研究結(jié)果可為遼西北地區(qū)水資源分配和管理提供科學(xué)參考。

    蒸散發(fā); 遙感估算;SEBS模型; 遼西北地區(qū)

    蒸散發(fā)(ET)包括地表、植物表面的蒸發(fā)量以及通過植物表面和體內(nèi)的蒸騰量,是水量循環(huán)之間的關(guān)鍵環(huán)節(jié),地表蒸散發(fā)通常是指土壤蒸發(fā)和植物蒸騰量之和[1-2]。它為農(nóng)作物的生長發(fā)育提供重要的水分和能量,是陸表生態(tài)體系和水文循環(huán)的紐帶,控制著重要的水文過程,對研究流域水文循環(huán)和水資源情況起著至關(guān)重要的作用[3]。因此,準(zhǔn)確估算區(qū)域地表蒸散發(fā)值,分析區(qū)域水資源分布情況成為研究熱點和難點。早期有關(guān)蒸散發(fā)的研究多為采用波文比法、渦度相關(guān)法以及土壤水量平衡法等樣地尺度、景觀尺度的測量,不能提供區(qū)域蒸散發(fā)的空間分布特征[4-5]。遙感數(shù)據(jù)能夠提供非均質(zhì)性的地表能量平衡過程和與地表覆蓋情況緊密相關(guān)的參數(shù),準(zhǔn)確有效地反映時空的連續(xù)性和動態(tài)性[6-8]。因此借助遙感手段對地表蒸散發(fā)進行估算是目前最前沿的研究方向,近年來出現(xiàn)了許多遙感反演蒸散發(fā)的模型。其中SEBS模型是基于SEBAL模型基礎(chǔ)上發(fā)展的遙感模型,具有較高的通量估算精度和應(yīng)用價值,是目前國內(nèi)外學(xué)者反演蒸散發(fā)量中應(yīng)用最為廣泛的模型之一。Vinukollu等通過SEBS模型和傳統(tǒng)的Penman-Monteith,Priestley-Taylor采用3種方法估算全球尺度的蒸散發(fā)值[9];Gokmen等基于SEBS模型,在水分脅迫條件下集成土壤水分參數(shù)估算蒸散發(fā)[10];Alkhaier等基于MODIS數(shù)據(jù),采用SEBS模型分析了溫度、地下水位深度、地表顯熱通量和日蒸散發(fā)之間的關(guān)系[11];陸靜等人基于MODIS數(shù)據(jù),使用SEBS模型估算中國南方亞熱帶常綠松樹種植園得蒸散發(fā)值,并驗證了該模型的適用性[12];呂明全等人基于遙感數(shù)據(jù)和氣象數(shù)據(jù)研究了1960—2013年三峽庫區(qū)蒸散發(fā)時空分布[13];李琴等人基于SEBS模型,選取30 a遙感數(shù)據(jù)反演中亞5國及新疆地區(qū)日實際蒸散發(fā)量[14];馬亮和魏光輝以新疆塔里木盆地西緣的阿克蘇地區(qū)為例,基于彭曼模型計算值為標(biāo)準(zhǔn),評價了H-S,P-T和McCloud模型的適用性[15]。

    遼西北地區(qū)是重干旱區(qū),屬暖溫帶,半濕潤、半干旱大陸性季風(fēng)氣候,降水量由東南向西北遞減,西部降水量最少,多年平均水資源量不足以維持原有生態(tài)系統(tǒng)的最低水資源量,屬于資源性缺水區(qū)域[16-17]。在已有的研究中,對遼西北地區(qū)蒸散發(fā)量的研究較少,因此,本文通過MODIS數(shù)據(jù),采用SEBS模型估算研究區(qū)蒸散發(fā)現(xiàn)狀,結(jié)合土地利用類型和彭曼模型估算其模擬精度,并進一步分析研究區(qū)近50 a蒸散發(fā)量分布特征,最后采用變標(biāo)度極差分析法對遼西北地區(qū)未來蒸散發(fā)趨勢加以分析,以期為區(qū)域水資源綜合管理、分配與規(guī)劃提供理論支持與科學(xué)依據(jù),對進一步研究該區(qū)域旱澇災(zāi)害特征并提供防御措施具有重要的研究意義。

    1 研究區(qū)域與研究方法

    1.1 區(qū)域概況及數(shù)據(jù)來源

    遼西北地區(qū)是遼寧省重點干旱區(qū)域,素有“十年九旱”之稱,主要包括朝陽、阜新、錦州、葫蘆島、鐵嶺5縣市以及沈陽的康平、法庫縣[18]。該區(qū)域多年平均降水量為450~700 mm左右,降水多集中在夏季,年際變化很大[19]。本文選取遼西北地區(qū)及其周邊27個氣象臺站50 a的氣象數(shù)據(jù)(圖1),數(shù)據(jù)來源于氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng),包括日照時數(shù)、日最高最低氣溫、平均相對濕度、平均風(fēng)速等,時間序列為1961—2010年。遙感數(shù)據(jù)來源于NASA數(shù)據(jù)中心,為了更準(zhǔn)確的估算2013年不同季節(jié)蒸散量的變化情況,本文選取2013年內(nèi)平均每月一景無云或少云遙感影像,應(yīng)用MODIS投影變換軟件MRT工具對數(shù)據(jù)進行重采樣及投影變換的處理,主要提取1~7,1~19,31,32波段。并用MODIS LDOPE工具對云量的分布進行檢測,進行除云處理工作,并評價各個波段數(shù)據(jù)的質(zhì)量。最終應(yīng)用內(nèi)嵌于ILWIS中的SEBS模型進行數(shù)據(jù)處理、應(yīng)用和分析。

    圖1 遼西北地區(qū)氣象臺站

    1.2 基于SEBS的遙感ET估算模型

    本文基于能量平衡原理及SEBS模型[20],得到能量平衡公式如下:

    λET=Rn-G-H

    (1)

    式中:Rn是地表凈輻射通量(W/m2);λET是潛熱通量(W/m2);λ為水的汽化熱(J/kg);ET為實際蒸散量[kg/(m2·s)];G是土壤熱通量(W/m2);H是顯熱通量(W/m2)。

    地表凈輻射通量(Rn)的公式如下:

    (2)

    式中:α為地表反照率;Rswd為向下的短波輻射;Rlwd為向下的長波輻射;ε為地表發(fā)射率,又稱輻射率;σ為Stefan-Bolzmann常數(shù),T0為遙感反演的地表輻射溫度(K)。其中α、ε、T0能從遙感數(shù)據(jù)的可見光至熱紅外波段得到。

    土壤熱通量(G0)的計算公式如下:

    G0=Rn×[Γc+(1-fc)×(Γs-Γc)]

    (3)

    式中:Rn為凈輻射通量(W/m2);Γc為植被在完全覆蓋的情況下,G0與凈輻射的比值;Γs為裸土條件下G0與凈輻射的比值;fc為植被覆蓋率。

    顯熱通量(H)的計算公式如下:

    (4)

    式中:ρ為空氣密度(kg/m3);CP為比熱容[J/(kg·k)];ΔT為溫度差(K);ra為空氣動力學(xué)阻抗(s/m)。

    通過以上公式,進而求出λET。

    1.3 Penman-Monteith公式

    本文依據(jù)FAO[21]推薦的Penman-Monteith公式計算遼西北地區(qū)的日蒸散發(fā)值,時間序列為1961—2010年。其公式如下:

    (5)

    式中:ET0為參考作物蒸散量(mm/d);Δ為飽和水汽壓曲線斜率(kpa/℃);Rn為地表凈輻射(MJ/(m2·d));G為土壤熱通量(MJ/(m2·d));r為干濕表常數(shù);T為2 m高處日均氣溫(℃);U2為2 m高處風(fēng)速(m/s);es為飽和水汽壓(kPa);ea為實際水汽壓(KPa)。

    2 結(jié)果與分析

    2.1 遼西北地區(qū)遙感反演地表蒸散發(fā)

    本文通過計算瞬時潛熱通量以及衛(wèi)星過境瞬間提供的每個參量,得到蒸發(fā)比;利用蒸發(fā)比,由地表通量和各項因子并通過日照時數(shù)等地表參數(shù)的引入,最終通過轉(zhuǎn)換獲得實際的日蒸散量。結(jié)果見圖2,春季蒸散發(fā)最高值最大,但全區(qū)范圍內(nèi)高值點較少,夏季值整體較高,遼西北地區(qū)2013年區(qū)域平均蒸散發(fā)量大小依次為夏季>春季>秋季>冬季。春季ET分布范圍較廣且均勻,其低值區(qū)大多分布在中北部,主要是因為該地植被覆蓋率低,植被稀疏,因此蒸散量較低;夏季植被覆蓋率最高,ET值明顯增加,其高值區(qū)分布在朝陽、阜新一帶,主要是因為該地區(qū)屬于遼西北重旱區(qū),夏季溫度高,蒸散發(fā)量明顯增加,低值區(qū)出現(xiàn)在錦州、葫蘆島等地,該地區(qū)靠近海洋,溫度相對低,蒸散發(fā)值也相對降低;秋季ET值沒有明顯下降;而冬季氣溫驟降,植被覆蓋率減少,ET值進一步降低,由西部向東部逐漸降低。綜上所述,蒸散發(fā)量季節(jié)變化較為明顯。由圖2可以看出,春季蒸散發(fā)量分布范圍較廣且均勻;夏季表現(xiàn)為中部地區(qū)相對東西部地區(qū)較高;秋季蒸散分布整體均勻,只有西部零星高值區(qū);冬季各地區(qū)出現(xiàn)差異,由西部向東部逐漸降低。遼西北地區(qū)季節(jié)內(nèi)蒸散發(fā)時空差異較大,易發(fā)生春旱、秋旱,城市缺水問題較為嚴(yán)重。因此,遼西北地區(qū)水資源配置應(yīng)考慮到區(qū)域蒸散發(fā)的四季變化情況進行合理優(yōu)化。

    圖2 遼西北地區(qū)不同季節(jié)日蒸散量(ET)空間分布

    2.2 遼西北地區(qū)不同土地利用類型蒸散發(fā)值

    為進一步探究遼西北地區(qū)蒸散發(fā)的空間分布特征,本文統(tǒng)計了該地區(qū)典型土地利用類型的蒸散發(fā)量,從箱狀圖(圖3)可見,雖然蒸散發(fā)量在季節(jié)上整體趨勢相同,但是不同土地類型的蒸散發(fā)量具有明顯的差異性,。遼西北地區(qū)土地主要以耕作為主,該地區(qū)是玉米主要種植區(qū)之一,因此其農(nóng)田的比例相對較大,其蒸散發(fā)值也較高。林地、草地植被覆蓋率次之,因此其蒸散發(fā)量也相對較高,各個土地利用類型蒸散發(fā)值季節(jié)分布大小為夏季>春季>秋季>冬季。各個季節(jié)不同土地利用類型的蒸散發(fā)量變化大小依次為農(nóng)田、草地、林地和水體。春季、夏季各土地利用類型蒸散發(fā)值變化幅度較小,變化不明顯,農(nóng)田蒸散發(fā)值最大為8 mm/d左右,水體蒸散發(fā)值最小為0.5 mm/d左右;秋季農(nóng)田蒸散發(fā)的最低值為3.24 mm/d,最高值為4.68 mm/d,落差較?。欢緶囟鹊?,植被凋落、水體冰凍,因此各個土地利用類型的蒸散發(fā)量最低。在各個土地利用類型中,水體蒸散發(fā)值最小,主要由于遼西北地區(qū)氣候干旱,水資源匱乏,經(jīng)濟發(fā)展水平較低,水體干枯,城市缺水現(xiàn)象嚴(yán)重,造成蒸散發(fā)量極少。

    2.3 遙感ET與傳統(tǒng)P-M模型值的對比

    為驗證SEBS模型估算蒸散發(fā)值的精度,本文根據(jù)2013年不同土地利用分布提取農(nóng)田中各日蒸散發(fā)值,結(jié)合FAO推薦的Penman-Monteith公式和作物參照法計算研究區(qū)與遙感影像所對應(yīng)的實際日蒸散發(fā)值,并求其平均值,與之進行對比,結(jié)果見圖4。SEBS模型的計算結(jié)果與實際測量的結(jié)果相差不大,誤差大約在7.5%左右,根據(jù)以往研究可知,本研究結(jié)果在允許接受的范圍內(nèi)[22]??傮w來說,本文利用ILWIS軟件,采用SEBS模型計算得到的遼西北地區(qū)2013年蒸散發(fā)時空分布上是合理可靠的。

    2.4 遼西北地區(qū)近50 a參考作物蒸散發(fā)值時空分布特征

    2.4.1 時間分析 本文利用均一化模型RHtest對氣象數(shù)據(jù)進行了均一化處理,然后根據(jù)Penman-Monteith公式計算了近50 a遼西北地區(qū)日參考作物蒸散發(fā)值。采用一元線性回歸法對不同季節(jié)全區(qū)蒸散發(fā)平均值進行分析(圖5,圖6),結(jié)果表明遼西北地區(qū)50 a來日均蒸散發(fā)值季節(jié)差異性較大,整體呈下降趨勢,但是各季節(jié)變化趨勢差異性較小。1980年蒸散發(fā)量下降速率最大,蒸散發(fā)值由高到低依次為夏季、春季、秋季和冬季。其中50 a間夏季和春季蒸散發(fā)值波動較大,冬季和秋季蒸散發(fā)值相對平緩。春、夏、秋和冬季全區(qū)平均蒸散發(fā)值分別以0.020 mm/50 a,0.025 mm/50 a,0.019 mm/50 a和0.011 mm/50 a的速率減少,與近50 a蒸散發(fā)季節(jié)變化趨勢一致。春季蒸散發(fā)傾向率區(qū)間為-0.054 8~0.001 6 mm/50 a;夏季蒸散發(fā)傾向率降幅區(qū)間為0.005 1~0.004 1 mm/50 a;秋季蒸散發(fā)傾向率值區(qū)間為-0.039 5~0.012 1 mm/50 a;冬季蒸散發(fā)傾向率值區(qū)間為-0.024 3~0.005 9 mm/50 a,且全區(qū)蒸散發(fā)傾向率均為中部地區(qū)最高,向四周遞減的趨勢。

    圖3 遼西北地區(qū)不同季節(jié)不同土地利用類型蒸散發(fā)值

    圖4 兩種方法下日蒸散發(fā)量結(jié)果對比

    圖5 1961-2010年遼西北地區(qū)近50 a參考作物蒸散發(fā)時間變化趨勢

    2.4.2 空間分析 本文利用反距離加權(quán)法對遼西北地區(qū)參考作物蒸散發(fā)值進行空間插值(圖7),春季蒸散發(fā)范圍為4.24~13.47 mm;夏季范圍為4.81~15.43 mm;秋季范圍為2.70~8.40 mm;冬季范圍為1.08~3.36 mm。遼西北地區(qū)蒸散發(fā)空間分布具有明顯的區(qū)域性,春、夏兩季蒸散發(fā)值由北部向南部逐漸減少,其較高的地區(qū)位于朝陽、阜新和康平縣一帶,該地區(qū)溫度較高且干燥,植被覆蓋率較多,因此其蒸散發(fā)值較高。秋季蒸散發(fā)值高值區(qū)位于遼西北西部和北部地區(qū),包括錦州、朝陽、建平等地。冬季蒸散發(fā)值高值區(qū)分布在遼西北西部地區(qū),由中北部向中南部逐漸降低??傮w而言4個季節(jié)蒸散發(fā)低值區(qū)均位于遼西北中南部黑山縣。且四季平均日蒸散量排序為夏季>春季>秋季>冬季,這是由于春季隨著氣溫的回升,農(nóng)作物耕種,此時蒸散發(fā)值較高,夏季氣溫升高,太陽輻射較強,降水充沛,農(nóng)作物逐漸成熟,其蒸散發(fā)量達到最大,到了9月份,氣溫逐漸降低,其蒸散發(fā)值也隨之減少,冬季寒冷,作物枯萎,需水量大幅度下降,相應(yīng)蒸散發(fā)量也達到最小。其分布規(guī)律與遙感反演的結(jié)果基本一致,且與蒸散發(fā)50 a間總體變化趨勢一致。

    2.4.3 未來變化趨勢 變標(biāo)度極差分析法(簡稱R/S分析法)是一種研究時間序列的分形理論,R/S→∞時,則說明所分析的時間序列存在Hurst現(xiàn)象,H稱為Hurst指數(shù)。當(dāng)0.5

    圖6 1961-2010年遼西北地區(qū)不同季節(jié)參考作物蒸散發(fā)傾向率變化

    圖7 遼西北地區(qū)不同季節(jié)參考作物蒸散值空間分布

    表1 各站點不同時間尺度的Hurst指數(shù)和分維數(shù)

    3 結(jié) 論

    (1) 遼西北地區(qū)下墊面蒸散發(fā)季節(jié)變化明顯,植被覆蓋情況和地表溫度是影響蒸散發(fā)的主要因素。春季ET分布范圍均勻,該地植被覆蓋率低,植被稀疏,因此蒸散量較低;夏季植被覆蓋率最高,ET值明顯增加;秋季植被覆蓋率較高,因此ET值沒有明顯下降;而冬季氣溫驟降,植被覆蓋率減少,ET值進一步降低??傮w來說遼西北地區(qū)ET值較均勻,沒有明顯的高低之分,在所研究站點的數(shù)據(jù)上分析能夠得出在遼西北地區(qū)偏西北地區(qū)蒸散發(fā)值相對較高,中南部地區(qū)相對較低,只有部分低值區(qū)呈零星的分布,其區(qū)域分異并不明顯。

    (2) 蒸散發(fā)量在不同土地類型中具有明顯的差異性,4個季節(jié)整體趨勢相同。蒸散發(fā)量的季節(jié)分布為夏季>春季>秋季>冬季,且各土地利用類型中蒸散發(fā)值由大到小分別為農(nóng)田、草地、林地和水體。

    (3) 本文利用Penman-Monteith公式和作物參照法計算的日蒸散發(fā)量驗證遙感反演的農(nóng)田日蒸散發(fā)量,其誤差值約為7.5%,證明遙感手段反演區(qū)域蒸散發(fā)量值方法合理可行,可為區(qū)域水資源合理優(yōu)化配置提供科學(xué)參考依據(jù)。

    傳統(tǒng)的蒸散估算模型多是基于觀測站單點或田間尺度的計算,采用經(jīng)驗統(tǒng)計模型,對地表觀測數(shù)據(jù)的依賴性過高。本文依據(jù)地表能量平衡方程,采用內(nèi)嵌于ILWIS軟件中的SEBS模型,對各個波段拉伸偏移量、水汽含量、氣溶膠光學(xué)厚度以及臭氧含量等參數(shù)進行校正,結(jié)合顯熱通量、土壤熱通量、凈輻射通量等地表參數(shù),進而計算區(qū)域蒸散發(fā)。結(jié)果表明采用SEBS模型有效消除了地表溫度和氣象數(shù)據(jù)對蒸散發(fā)估算帶來的不確定性,具有較高的精度。同時選用Penman-Manteith公式計算遼西北地區(qū)日潛在蒸散發(fā),對模型模擬的結(jié)果進行驗證,得到合理可靠的結(jié)果,其誤差在前人所研究的允許誤差之內(nèi)[22]。

    本文也存在一定不足,由于數(shù)據(jù)的缺失,站點氣象數(shù)據(jù)不完善,彭曼公式所采用的氣象數(shù)據(jù),只是對點上蒸散發(fā)的估算,不能對大區(qū)域結(jié)果進行驗證。但文中采用此模型對遼西北地區(qū)蒸散發(fā)進行估算及驗證,仍可為該地區(qū)有關(guān)蒸散發(fā)的研究提供科學(xué)依據(jù),日后將重點研究基于長時間序列的遙感數(shù)據(jù),使用SEBS模型的蒸散發(fā)時空特征分析。

    [1] 施云霞.基于SEBAL模型的新疆精河綠洲蒸散發(fā)研究[D].烏魯木齊:新疆師范大學(xué),2016.

    [2] 黨素珍,蔣曉輝,董國濤,等.涇河上游流域?qū)嶋H蒸散發(fā)變化趨勢及成因分析[J].水土保持研究,2016,23(2):143-147.

    [3] 湯秋鴻,黃忠偉,劉星才,等.人類用水活動對大尺度陸地水循環(huán)的影響[J].地球科學(xué)進展,2015(10):1091-1099.

    [4] 李洋,陳亞寧,李衛(wèi)紅,等.塔里木河下游胡楊群落的蒸散發(fā)觀測研究[J].新疆環(huán)境保護,2013,35(1):1-7.

    [5] 李念,孫維君,秦翔,等.祁連山老虎溝地區(qū)高寒草甸蒸散發(fā)估算[J].干旱區(qū)資源與環(huán)境,2016(6):173-178.

    [6] 張榮華,杜君平,孫睿.區(qū)域蒸散發(fā)遙感估算方法及驗證綜述[J].地球科學(xué)進展,2012,27(12):1295-1307.

    [7] 王健美,楊國范,周林滔.基于SEBS模型反演凌河流域尺度地表蒸散發(fā)量[J].灌溉排水學(xué)報,2016,35(2):95-99.

    [8] 白娟,楊勝天,董國濤,等.基于多源遙感數(shù)據(jù)的三江平原日蒸散量估算[J].水土保持研究,2013,20(3):190-195.

    [9] Vinukollu R K, Wood E F, Ferguson C R, et al. Global estimates of evapotranspiration for climate studies using multi-sensor remote sensing data:Evaluation of three process-based approaches[J]. Remote Sensing of Environment, 2011,115(3):801-823.

    [10] Gokmen M, Vekerdy Z, Verhoef A, et al. Integration of soil moisture in SEBS for improving evapotranspiration estimation under water stress conditions[J]. Remote Sensing of Environment, 2012,121(138):261-274.

    [11] Alkhaier F, Su Z, Flerchinger G N. Reconnoitering the effect of shallow groundwater on land surface temperature and surface energy balance using MODIS and SEBS[J]. Hydrology & Earth System Sciences & Discussions, 2012,16(7):1833-1844.

    [12] Lu J, Li Z L, Tang R, et al. Evaluating the SEBS-estimated evaporative fraction from MODIS data for a complex underlying surface[J]. Hydrological Processes, 2013,27(22):3139-3149.

    [13] Lv Ming-quan, Chen Ji-long, Mirza Z A, et al. Spatial distribution and temporal variation of reference evapotranspiration in the Three Gorges Reservoir area during 1960—2013 [J]. International Journal of Climatology, 2016,36(14):4497-4511.

    [14] 李琴,陳曦,包安明,等.基于SEBS模型干旱區(qū)蒸散發(fā)量研究[J].遙感技術(shù)與應(yīng)用,2014,29(2):195-201.

    [15] 馬亮,魏光輝.新疆塔里木盆地西緣參考作物蒸散發(fā)模型的適用性評價[J].干旱區(qū)資源與環(huán)境,2015,29(8):132-137.

    [16] 鄭帥,李德明,郭成久,等.遼西北地區(qū)風(fēng)力侵蝕過程分析[J].水土保持研究,2015,22(4):138-142.

    [17] 張藝萌,張雪松,郭婷婷,等.遼西北地區(qū)氣溫和降水變化對氣候生產(chǎn)潛力的影響[J].中國農(nóng)業(yè)氣象,2015(2):203-211.

    [18] 梁文章,孫玉華,孫寶杰.遼西北地區(qū)水資源量趨勢分析[J].東北水利水電,2011,37(2):37-38.

    [19] 張繼權(quán),嚴(yán)登華,王春乙,等.遼西北地區(qū)農(nóng)業(yè)干旱災(zāi)害風(fēng)險評價與風(fēng)險區(qū)劃研究[J].防災(zāi)減災(zāi)工程學(xué)報,2012,32(3):300-306.

    [20] 曾麗紅,宋開山,張柏,等. SEBAL模型與MODIS產(chǎn)品的松嫩平原蒸散量研究[J].干旱區(qū)資源與環(huán)境,2011,25(1):140-147.

    [21] Djaman K, Irmak S, Kabenge I, et al. Evaluation of FAO-56Penman-Monteith Model with Limited Data and the Valiantzas Models for Estimating Grass-Reference Evapotranspiration in Sahelian Conditions[J]. Journal of Irrigation and Drainage Engineering, 2016,142(11):04016044.

    [22] 曹永強,張亭亭,徐丹,等.海河流域蒸散發(fā)時空演變規(guī)律分析[J].資源科學(xué),2014,36(7):1489-1499.

    EstimationofEvapotranspirationandAnalysisofTemporal-SpaceCharacteristicinNorthwestofLiaoningProvinceBasedonSEBSModel

    WANG Rui1, ZHANG Jiquan1, CAO Yongqiang2, GUO Enliang1, LI Danjun1, MA Qiyun1

    (1.SchoolofEnvironment/NaturalDisasterResearchInstitute,NortheastNormalUniversity,Changchun130024,China;2.SchoolofUrbanPlanningandEnvironmentalScience,LiaoningNormalUniversity,Dalian,Liaoning116029,China)

    Terrestrial evapotranspiration is the key connection of a wide array of Earth system processes, which links water cycles and enemy balances over land. This paper selected 12 sensing data from NASA as the representative of four seasons, combined with the SEBS model to reckon the evapotranspiration of land surface and analyzed the ET of different land use types to verify the precision of remote sensing. To analyze the temporal-spatial characteristic, trend analysis, wavelet analysis, IDW and R/S analysis method based on Penman-Monteith wereselected. The results show that the ET in northwest of Liaoning Province decreased in the order: summer>autumn>spring>winter; the evapotranspriration decreased in the order: summer>spring>autumn>winter; the error of evapotranspiration estimating is about 7.5% in meantime, indicating that the method is reasonable to be applied in the study area. The results can provide a scientific reference for water resources allocation and management in northwest of Liaoning Province.

    evapotranspiration; estimating by remote sensing; SEBS model; northwest of Liaoning Province

    P426.2;TP79

    A

    1005-3409(2017)06-0382-06

    2016-12-16

    2016-12-29

    中央高?;究蒲袠I(yè)務(wù)費專項資金(2412016KJ046);公益性行業(yè)(氣象)科研專項儲備項目(2015001);國家自然科學(xué)基金(41571491);國家自然科學(xué)基金(51279072)

    王蕊(1993—),女,遼寧省錦州市人,碩士研究生,主要從事農(nóng)業(yè)氣象災(zāi)害風(fēng)險評估等方面的研究。E-mail:wangr967@nenu.edu.cn

    張繼權(quán)(1965—),男,吉林省九臺市人,教授,博士學(xué)位,主要從事環(huán)境災(zāi)害風(fēng)險評估與管理等方面研究。E-mail:zhangjq022@nenu.edu.cn

    猜你喜歡
    西北地區(qū)通量反演
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    冬小麥田N2O通量研究
    近60年西北地區(qū)森林資源變化分析與思考
    三集一體機在西北地區(qū)適用性分析
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    化工進展(2015年6期)2015-11-13 00:26:29
    西北地區(qū)首屆課博會在寶雞召開
    疊前同步反演在港中油田的應(yīng)用
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    免费观看精品视频网站| 亚洲欧洲精品一区二区精品久久久| 这个男人来自地球电影免费观看| 人人妻人人澡人人看| 亚洲情色 制服丝袜| 夜夜躁狠狠躁天天躁| 制服人妻中文乱码| 一边摸一边做爽爽视频免费| 18禁国产床啪视频网站| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 成人国产综合亚洲| 两个人免费观看高清视频| 亚洲欧美一区二区三区黑人| www国产在线视频色| 成人精品一区二区免费| 夜夜看夜夜爽夜夜摸| 男女做爰动态图高潮gif福利片 | 一本久久中文字幕| 久久人人爽av亚洲精品天堂| 中文亚洲av片在线观看爽| 亚洲精品久久成人aⅴ小说| 精品免费久久久久久久清纯| 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 男人舔女人下体高潮全视频| 国产成人精品无人区| 在线av久久热| 国产精品国产高清国产av| 天天躁狠狠躁夜夜躁狠狠躁| 两个人视频免费观看高清| 操美女的视频在线观看| 免费在线观看视频国产中文字幕亚洲| 国产又爽黄色视频| 老鸭窝网址在线观看| 精品久久蜜臀av无| 亚洲av五月六月丁香网| 他把我摸到了高潮在线观看| 大香蕉久久成人网| 国产片内射在线| 欧美精品亚洲一区二区| 国产精品,欧美在线| 夜夜爽天天搞| 亚洲天堂国产精品一区在线| 亚洲精品中文字幕在线视频| 琪琪午夜伦伦电影理论片6080| 午夜精品久久久久久毛片777| 久久国产精品人妻蜜桃| 黄色 视频免费看| 亚洲男人天堂网一区| 18禁裸乳无遮挡免费网站照片 | 看免费av毛片| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 99久久综合精品五月天人人| 欧美精品亚洲一区二区| 88av欧美| 欧美日韩福利视频一区二区| 国产熟女xx| 精品人妻1区二区| 大陆偷拍与自拍| 国内精品久久久久久久电影| 欧美 亚洲 国产 日韩一| 国产野战对白在线观看| 亚洲一码二码三码区别大吗| 国产片内射在线| 亚洲欧美精品综合一区二区三区| 亚洲,欧美精品.| 午夜免费激情av| 自拍欧美九色日韩亚洲蝌蚪91| 国产成+人综合+亚洲专区| 涩涩av久久男人的天堂| 亚洲av第一区精品v没综合| 可以免费在线观看a视频的电影网站| 一本综合久久免费| 日韩欧美一区视频在线观看| 青草久久国产| 成人精品一区二区免费| 国产xxxxx性猛交| 亚洲七黄色美女视频| 中出人妻视频一区二区| 久热爱精品视频在线9| 国产蜜桃级精品一区二区三区| 一边摸一边抽搐一进一小说| 国产激情欧美一区二区| 久久人人精品亚洲av| 男女之事视频高清在线观看| 亚洲欧洲精品一区二区精品久久久| 国产亚洲av高清不卡| 精品国产美女av久久久久小说| 视频在线观看一区二区三区| 看免费av毛片| 精品人妻1区二区| 日韩欧美免费精品| 伦理电影免费视频| 大型黄色视频在线免费观看| 狠狠狠狠99中文字幕| 国产亚洲欧美精品永久| 亚洲一区二区三区色噜噜| 天堂√8在线中文| 高清在线国产一区| 亚洲精品在线美女| 亚洲精品国产精品久久久不卡| 日韩有码中文字幕| 超碰成人久久| 国产蜜桃级精品一区二区三区| 在线观看午夜福利视频| 久久 成人 亚洲| 天天躁夜夜躁狠狠躁躁| 一夜夜www| 国产成人欧美在线观看| 少妇 在线观看| 亚洲视频免费观看视频| 国产精品美女特级片免费视频播放器 | 亚洲熟女毛片儿| 99国产精品免费福利视频| 别揉我奶头~嗯~啊~动态视频| 成人欧美大片| 国产xxxxx性猛交| 亚洲精华国产精华精| 国产激情欧美一区二区| 国产真人三级小视频在线观看| 亚洲av熟女| svipshipincom国产片| 亚洲熟女毛片儿| 99在线人妻在线中文字幕| 国产国语露脸激情在线看| av电影中文网址| 精品国产一区二区三区四区第35| 国产av又大| 久久国产亚洲av麻豆专区| 此物有八面人人有两片| av网站免费在线观看视频| 级片在线观看| 午夜福利高清视频| 国产成人精品无人区| 亚洲天堂国产精品一区在线| 国产精品影院久久| www.www免费av| 看免费av毛片| 久久国产精品男人的天堂亚洲| av网站免费在线观看视频| 精品第一国产精品| 18禁观看日本| 免费人成视频x8x8入口观看| 国产99久久九九免费精品| 国产精品久久久久久亚洲av鲁大| 男女下面插进去视频免费观看| 高清在线国产一区| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 看免费av毛片| 久久狼人影院| 久99久视频精品免费| 精品久久久精品久久久| 欧美激情 高清一区二区三区| 精品熟女少妇八av免费久了| 国产精品99久久99久久久不卡| 91精品三级在线观看| 日韩一卡2卡3卡4卡2021年| 欧美精品啪啪一区二区三区| 19禁男女啪啪无遮挡网站| 性欧美人与动物交配| 97人妻天天添夜夜摸| 亚洲久久久国产精品| 精品国产一区二区久久| 久热爱精品视频在线9| 啦啦啦 在线观看视频| 十分钟在线观看高清视频www| 色老头精品视频在线观看| 日韩欧美三级三区| 亚洲第一欧美日韩一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 国内久久婷婷六月综合欲色啪| 免费在线观看影片大全网站| 亚洲一码二码三码区别大吗| 91精品国产国语对白视频| 在线视频色国产色| 国产精华一区二区三区| 欧美绝顶高潮抽搐喷水| 长腿黑丝高跟| 亚洲男人的天堂狠狠| 丝袜在线中文字幕| a级毛片在线看网站| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 一区二区三区精品91| 色播在线永久视频| 欧美老熟妇乱子伦牲交| 黄色片一级片一级黄色片| 不卡av一区二区三区| 国产91精品成人一区二区三区| 在线观看免费午夜福利视频| 欧美精品亚洲一区二区| 丝袜美腿诱惑在线| 啦啦啦观看免费观看视频高清 | 欧美大码av| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 日本三级黄在线观看| 免费在线观看影片大全网站| 国产精品 国内视频| 日日夜夜操网爽| 成人三级黄色视频| 久99久视频精品免费| 又黄又爽又免费观看的视频| 久久久水蜜桃国产精品网| 中文字幕色久视频| 精品电影一区二区在线| 日韩高清综合在线| 亚洲一区中文字幕在线| 天天添夜夜摸| 免费观看精品视频网站| 午夜免费观看网址| 婷婷丁香在线五月| 女人高潮潮喷娇喘18禁视频| 久久久久九九精品影院| 曰老女人黄片| 嫩草影院精品99| 最新美女视频免费是黄的| 亚洲少妇的诱惑av| 国产熟女午夜一区二区三区| 美女国产高潮福利片在线看| 露出奶头的视频| 欧美国产日韩亚洲一区| 亚洲国产精品sss在线观看| 成人手机av| 一本大道久久a久久精品| 亚洲第一青青草原| www.精华液| 色哟哟哟哟哟哟| 国产成人精品在线电影| 老司机午夜十八禁免费视频| 一级毛片女人18水好多| 日韩精品青青久久久久久| 亚洲aⅴ乱码一区二区在线播放 | 少妇的丰满在线观看| 麻豆久久精品国产亚洲av| 免费观看人在逋| 黄色丝袜av网址大全| 久久久久亚洲av毛片大全| 高清毛片免费观看视频网站| www日本在线高清视频| 天堂√8在线中文| 法律面前人人平等表现在哪些方面| 男人舔女人的私密视频| tocl精华| 在线观看免费视频日本深夜| 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区av网在线观看| 色在线成人网| 国产成人av教育| 久久国产精品影院| 日日干狠狠操夜夜爽| 中文字幕av电影在线播放| 色综合亚洲欧美另类图片| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看 | 99国产精品99久久久久| 人人妻人人爽人人添夜夜欢视频| 99久久综合精品五月天人人| 村上凉子中文字幕在线| 午夜老司机福利片| 久久香蕉精品热| 日日干狠狠操夜夜爽| 一级毛片高清免费大全| 亚洲熟妇中文字幕五十中出| 九色国产91popny在线| 男人舔女人下体高潮全视频| 精品久久蜜臀av无| 91成人精品电影| 伦理电影免费视频| 国产在线观看jvid| 亚洲一区高清亚洲精品| 亚洲精品在线观看二区| 18禁黄网站禁片午夜丰满| 一级a爱视频在线免费观看| 成在线人永久免费视频| 真人做人爱边吃奶动态| 久久精品影院6| 色尼玛亚洲综合影院| 村上凉子中文字幕在线| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区视频了| a级毛片在线看网站| 熟妇人妻久久中文字幕3abv| 久久中文字幕一级| 国产1区2区3区精品| 1024香蕉在线观看| 高清在线国产一区| 国产免费av片在线观看野外av| 国产区一区二久久| 老司机靠b影院| 国产成人影院久久av| 欧美日韩瑟瑟在线播放| 丝袜美腿诱惑在线| av电影中文网址| 琪琪午夜伦伦电影理论片6080| 国产在线精品亚洲第一网站| 久久天堂一区二区三区四区| 亚洲第一电影网av| 午夜福利一区二区在线看| 久久久国产成人精品二区| 天天躁夜夜躁狠狠躁躁| 可以在线观看的亚洲视频| 动漫黄色视频在线观看| 国产精品久久久久久人妻精品电影| 91九色精品人成在线观看| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 人人妻人人爽人人添夜夜欢视频| 国产在线精品亚洲第一网站| 岛国视频午夜一区免费看| 国产精品永久免费网站| 国产精品久久久人人做人人爽| 老司机福利观看| 亚洲精品一区av在线观看| 99精品久久久久人妻精品| 成人特级黄色片久久久久久久| 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 欧美成人性av电影在线观看| 中文字幕精品免费在线观看视频| 成在线人永久免费视频| 女人爽到高潮嗷嗷叫在线视频| 女性生殖器流出的白浆| 成人精品一区二区免费| 国产精品爽爽va在线观看网站 | 国产人伦9x9x在线观看| 亚洲国产精品999在线| 很黄的视频免费| 免费无遮挡裸体视频| 国产一区二区三区综合在线观看| 脱女人内裤的视频| 国产三级黄色录像| 免费一级毛片在线播放高清视频 | 中文字幕精品免费在线观看视频| 日本免费a在线| netflix在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产视频一区二区在线看| 99国产精品一区二区蜜桃av| 亚洲一区中文字幕在线| 母亲3免费完整高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国内毛片毛片毛片毛片毛片| 成人永久免费在线观看视频| 国产亚洲精品久久久久5区| 在线观看舔阴道视频| 亚洲无线在线观看| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 久久中文字幕一级| 国产精品乱码一区二三区的特点 | 精品少妇一区二区三区视频日本电影| 欧美成人一区二区免费高清观看 | 久久香蕉激情| 一区二区三区高清视频在线| 欧美日韩一级在线毛片| 亚洲av成人不卡在线观看播放网| 精品一区二区三区四区五区乱码| av天堂在线播放| 国产成人av教育| 嫩草影视91久久| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久亚洲av鲁大| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 国产精品亚洲av一区麻豆| 波多野结衣av一区二区av| 亚洲人成伊人成综合网2020| 亚洲欧美日韩无卡精品| 亚洲第一欧美日韩一区二区三区| avwww免费| 成人亚洲精品av一区二区| 久久狼人影院| 天堂影院成人在线观看| 99热只有精品国产| 免费在线观看日本一区| 久久久久久大精品| 亚洲专区字幕在线| 日本a在线网址| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 久久国产精品影院| tocl精华| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 精品久久久久久,| 老司机午夜福利在线观看视频| 啦啦啦免费观看视频1| 一卡2卡三卡四卡精品乱码亚洲| 国产xxxxx性猛交| 一卡2卡三卡四卡精品乱码亚洲| 日本 av在线| 黄色丝袜av网址大全| 国产一卡二卡三卡精品| 人妻丰满熟妇av一区二区三区| 黄色视频不卡| 91九色精品人成在线观看| 制服丝袜大香蕉在线| 国产精品乱码一区二三区的特点 | 老司机靠b影院| 在线观看免费午夜福利视频| 免费看a级黄色片| 少妇的丰满在线观看| 午夜福利欧美成人| 天堂影院成人在线观看| 国产一区二区三区综合在线观看| 亚洲国产欧美网| 叶爱在线成人免费视频播放| 中文字幕精品免费在线观看视频| 一个人观看的视频www高清免费观看 | 国产精品,欧美在线| 久久久久久久久免费视频了| 国产亚洲欧美在线一区二区| 久久精品国产综合久久久| 久久中文字幕一级| 97人妻天天添夜夜摸| 高清毛片免费观看视频网站| 欧美一区二区精品小视频在线| 无人区码免费观看不卡| 久久久国产欧美日韩av| 亚洲男人的天堂狠狠| 久久天堂一区二区三区四区| 99精品在免费线老司机午夜| 大陆偷拍与自拍| 在线十欧美十亚洲十日本专区| 老鸭窝网址在线观看| 亚洲中文字幕日韩| 国产国语露脸激情在线看| 最好的美女福利视频网| 午夜福利成人在线免费观看| 国产精品免费一区二区三区在线| 黄色毛片三级朝国网站| 久久人妻熟女aⅴ| 亚洲国产精品合色在线| 大香蕉久久成人网| 99久久久亚洲精品蜜臀av| 国产色视频综合| 999久久久国产精品视频| 夜夜夜夜夜久久久久| av在线播放免费不卡| 天堂影院成人在线观看| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 一夜夜www| 一级a爱视频在线免费观看| 香蕉国产在线看| 此物有八面人人有两片| 亚洲av电影不卡..在线观看| 宅男免费午夜| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| 手机成人av网站| 国产亚洲av嫩草精品影院| 少妇粗大呻吟视频| 18禁黄网站禁片午夜丰满| 精品久久久久久,| svipshipincom国产片| 亚洲久久久国产精品| 亚洲国产欧美网| 高清黄色对白视频在线免费看| 国产精品亚洲一级av第二区| 亚洲人成电影免费在线| 搡老岳熟女国产| 亚洲天堂国产精品一区在线| 久久久国产精品麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产看品久久| 精品一区二区三区av网在线观看| 国产av精品麻豆| 欧美精品亚洲一区二区| av福利片在线| 亚洲第一av免费看| 亚洲在线自拍视频| 国产av精品麻豆| 亚洲精品一卡2卡三卡4卡5卡| 亚洲色图 男人天堂 中文字幕| 午夜免费成人在线视频| 多毛熟女@视频| 国产精品二区激情视频| 91在线观看av| 免费看十八禁软件| 波多野结衣高清无吗| 校园春色视频在线观看| 天堂√8在线中文| 国产成+人综合+亚洲专区| 18美女黄网站色大片免费观看| 久久久久久人人人人人| 老司机福利观看| 日本精品一区二区三区蜜桃| 欧美乱码精品一区二区三区| 1024视频免费在线观看| 一进一出抽搐gif免费好疼| 国产97色在线日韩免费| 18禁国产床啪视频网站| netflix在线观看网站| 日本 欧美在线| 少妇粗大呻吟视频| 午夜福利成人在线免费观看| АⅤ资源中文在线天堂| 色婷婷久久久亚洲欧美| 两性夫妻黄色片| 国产亚洲精品第一综合不卡| 国产aⅴ精品一区二区三区波| 国产精品日韩av在线免费观看 | 桃红色精品国产亚洲av| 午夜免费成人在线视频| 久久久久久亚洲精品国产蜜桃av| 欧美av亚洲av综合av国产av| 欧美国产日韩亚洲一区| 午夜福利免费观看在线| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 欧美日韩精品网址| 国产极品粉嫩免费观看在线| av在线播放免费不卡| 久久青草综合色| 51午夜福利影视在线观看| 精品国产亚洲在线| 久久精品国产99精品国产亚洲性色 | 久久久久国内视频| 国产av在哪里看| 国产人伦9x9x在线观看| 国产精品精品国产色婷婷| 亚洲精品一卡2卡三卡4卡5卡| 一进一出抽搐动态| 国产高清激情床上av| 亚洲色图综合在线观看| 亚洲黑人精品在线| 97碰自拍视频| 久久久久精品国产欧美久久久| 久久久国产精品麻豆| 亚洲精品一卡2卡三卡4卡5卡| 黄色片一级片一级黄色片| 国产高清激情床上av| 99国产精品99久久久久| 国产精品久久视频播放| 波多野结衣高清无吗| 国产高清videossex| 嫩草影院精品99| 久久香蕉精品热| 亚洲一区高清亚洲精品| 亚洲精品美女久久久久99蜜臀| 中文字幕精品免费在线观看视频| 国产av又大| 美女午夜性视频免费| 国产精品野战在线观看| 一级,二级,三级黄色视频| www.自偷自拍.com| 中文字幕最新亚洲高清| 免费高清在线观看日韩| 日本黄色视频三级网站网址| 女人被狂操c到高潮| 亚洲精品美女久久久久99蜜臀| 757午夜福利合集在线观看| 久久香蕉激情| 精品人妻在线不人妻| 精品国产国语对白av| 桃红色精品国产亚洲av| 一级a爱片免费观看的视频| 啦啦啦观看免费观看视频高清 | 黄色 视频免费看| 一卡2卡三卡四卡精品乱码亚洲| av有码第一页| 黑人巨大精品欧美一区二区mp4| 国产区一区二久久| 免费高清视频大片| 又黄又爽又免费观看的视频| 一进一出抽搐动态| 国产一卡二卡三卡精品| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成av片中文字幕在线观看| 色婷婷久久久亚洲欧美| 国产麻豆成人av免费视频| 午夜福利一区二区在线看| 在线观看舔阴道视频| 香蕉国产在线看| 国产精品香港三级国产av潘金莲| 国产成人精品久久二区二区免费| 久久久久九九精品影院| 亚洲第一av免费看| 波多野结衣巨乳人妻| 高清毛片免费观看视频网站| 美女免费视频网站| 真人做人爱边吃奶动态| 男人的好看免费观看在线视频 | 亚洲精品在线观看二区| 久久久久久久久免费视频了| cao死你这个sao货| av天堂在线播放| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 最好的美女福利视频网| 久久青草综合色| 1024视频免费在线观看| 亚洲国产精品sss在线观看| 国产精品野战在线观看| 国产成人影院久久av| 国产又爽黄色视频| 国产成人精品久久二区二区免费| 久久性视频一级片| avwww免费| 91九色精品人成在线观看| 久久久国产成人免费| 免费av毛片视频| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色视频不卡| 国产精品亚洲一级av第二区| 91麻豆av在线| 亚洲av日韩精品久久久久久密|