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

    新疆精河流域?qū)嶋H蒸散發(fā)時(shí)空變化特征

    2017-07-31 21:02:11代鵬超牛蘇娟毋兆鵬施云霞
    關(guān)鍵詞:精河流域研究

    代鵬超,牛蘇娟,毋兆鵬①,施云霞

    (1.新疆師范大學(xué)地理科學(xué)與旅游學(xué)院,新疆 烏魯木齊 830054;2.新疆干旱區(qū)湖泊環(huán)境與資源實(shí)驗(yàn)室,新疆 烏魯木齊 830054)

    新疆精河流域?qū)嶋H蒸散發(fā)時(shí)空變化特征

    代鵬超1,2,牛蘇娟1,毋兆鵬1,2①,施云霞1

    (1.新疆師范大學(xué)地理科學(xué)與旅游學(xué)院,新疆 烏魯木齊 830054;2.新疆干旱區(qū)湖泊環(huán)境與資源實(shí)驗(yàn)室,新疆 烏魯木齊 830054)

    精河流域是新疆天山北坡經(jīng)濟(jì)帶的重要組成部分。利用近60 a氣象數(shù)據(jù)及Landsat 5 TM影像數(shù)據(jù),采用Penman-Monteith公式和SEBAL模型、morlet小波分析和M-K突變檢驗(yàn),研究了精河流域?qū)嶋H蒸散量時(shí)空格局、變化特征及周期性。結(jié)果表明:(1)全流域?qū)嶋H蒸散量的時(shí)空分布變化受到氣象要素及地表下墊面的影響,近60 a蒸散量整體呈顯著波動減小趨勢,實(shí)際蒸散量變化速率在季節(jié)上表現(xiàn)為夏季>春季>秋季>冬季。(2)實(shí)際蒸散量于1982年突變后年平均蒸散量減少150.654 mm(17.361%),且蒸散量變化存在以29 a為主的3個(gè)振蕩周期,未來15 a蒸散量將呈先小幅上升后下降的趨勢。(3)蒸散量的空間分布變化與地表土地利用類型具有顯著相關(guān)關(guān)系,相關(guān)性表現(xiàn)為水域>耕地>林地>草地>未利用地。

    精河流域;實(shí)際蒸散發(fā);Penman-Monteith公式;SEBAL模型

    蒸散發(fā)是陸面水循環(huán)中重要的水文過程[1],也是聯(lián)系植物氣孔行為、碳交換和水分利用的關(guān)鍵生態(tài)過程[2],更是評價(jià)區(qū)域水熱平衡的重要參數(shù)和評價(jià)因子。作為大氣和土壤之間的“橋梁”,蒸散量是唯一同時(shí)出現(xiàn)在地表水量和能量平衡方程中的水文參量,是陸面過程數(shù)值模擬研究中不可缺少的邊界條件[3]。蒸散量估算是農(nóng)業(yè)、水文、氣象和土壤等學(xué)科的重要研究內(nèi)容,準(zhǔn)確估算蒸散量對生態(tài)系統(tǒng)管理、水土資源評價(jià)與利用、農(nóng)業(yè)生產(chǎn)都具有重要意義[4]。

    利用傳統(tǒng)的水文氣象學(xué)方法進(jìn)行蒸散量觀測時(shí),大多數(shù)只能在“點(diǎn)”或景觀尺度上進(jìn)行,很難應(yīng)用到大的區(qū)域尺度[5]。20世紀(jì)60年代遙感技術(shù)在地理研究中的應(yīng)用為定量估算更大區(qū)域的蒸散量提供了新思路。將具有堅(jiān)實(shí)理論基礎(chǔ)和明確物理概念的、傳統(tǒng)點(diǎn)尺度上的模型如Penman-Monteith公式等,與遙感定量反演技術(shù)結(jié)合形成互補(bǔ)相關(guān)模型,則既可以擴(kuò)展傳統(tǒng)“點(diǎn)”方法的空間尺度,又可提升遙感反演結(jié)果的時(shí)間分辨率,近年來已成為研究熱點(diǎn)[6]。

    精河流域地處我國西北干旱區(qū),是天山北坡經(jīng)濟(jì)帶的重要組成部分,同時(shí)也是我國生態(tài)環(huán)境問題最突出的內(nèi)陸河流域[7]。近年來,該地區(qū)大力發(fā)展工業(yè)、種植經(jīng)濟(jì)作物以提高人們生活水平,但由于水資源的過度開發(fā)利用,加上近幾年徑流性水資源的減少,導(dǎo)致流域內(nèi)生態(tài)逐年退化,地下水水位持續(xù)下降,天然植被嚴(yán)重退化,進(jìn)而引發(fā)了土地荒漠化、鹽漬化等一系列生態(tài)問題。該研究通過準(zhǔn)確估算精河流域?qū)嶋H蒸散量,分析未來水分耗散趨勢,為合理規(guī)劃精河流域農(nóng)業(yè)生產(chǎn)、高效利用和適度分配水資源提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    精河流域作為新疆第2大生態(tài)退化區(qū),位于新疆維吾爾自治區(qū)西北部準(zhǔn)噶爾盆地西南邊緣,天山支脈婆羅科努山北麓,介于44°00′21″~45°00′56″ N,81°07′52″~83°05′48″ E之間,總面積2 572.30 km2。經(jīng)濟(jì)以農(nóng)為主,多年形成了“一紅(枸杞)二白(棉花、鹽業(yè))一綠(瓜果蔬菜)和突破發(fā)展奶牛業(yè)”的經(jīng)濟(jì)發(fā)展戰(zhàn)略格局,棉花為當(dāng)?shù)刂еa(chǎn)業(yè)。

    由于地處亞歐大陸腹地,遠(yuǎn)離海洋,加之高大山脈的阻擋,導(dǎo)致該區(qū)夏季降水稀少,冬季干燥寒冷,大陸性氣候特征顯著。年均氣溫7.3 ℃,平原區(qū)多年平均降水量91 mm。該區(qū)是北疆沿天山一帶風(fēng)沙天氣最多的區(qū)域,大風(fēng)天氣多,且持續(xù)時(shí)間較長。干燥的氣候和強(qiáng)勁、活躍的大風(fēng)造成了強(qiáng)烈的蒸散發(fā)。

    1.2 數(shù)據(jù)材料與處理

    研究所用氣象資料來自1953—2015年精河流域地面氣象觀測站日值數(shù)據(jù)。個(gè)別缺測日數(shù)據(jù)采用相關(guān)分析進(jìn)行內(nèi)補(bǔ)。年資料統(tǒng)計(jì)按照自然年算法,即每年1—12月之和。季節(jié)的劃分采用氣象季節(jié),即以上一年12月和當(dāng)年1—2月為冬季,3—5月為春季,6—8月為夏季,9—11月為秋季。

    研究所用空間數(shù)據(jù)為Landsat 5 TM影像數(shù)據(jù),時(shí)間包括1998年9月25日、2007年9月18日和2011年9月13日3期。所有數(shù)據(jù)在ENVI 5.1軟件平臺中預(yù)處理后進(jìn)行空間建模。

    1.3 時(shí)間尺度蒸散量反演

    為了較好地在長時(shí)間尺度上獲得連續(xù)、穩(wěn)定的蒸散量數(shù)據(jù),首先利用聯(lián)合國糧農(nóng)組織推薦的Penman-Monteith公式[8]計(jì)算精河流域參考作物蒸散量(ET0,TE,0)。

    (1)

    式(1)中,TE,0為參考作物蒸散量,mm·d-1;Rn為冠層表面凈輻射,MJ·m-2·d-1;G為土壤熱通量,MJ·m-2·d-1;Δ為飽和水氣壓曲線在氣溫為T時(shí)的斜率,kPa·℃-1;T為日平均氣溫,℃;γ為干濕表常數(shù),kPa·℃-1;U2為距地面2 m處的風(fēng)速,m·s-1;ea為飽和水汽壓,kPa;ed為實(shí)際水汽壓,kPa。

    在ET0計(jì)算基礎(chǔ)上,借助表征作物在不同生長時(shí)期中需水量的生物學(xué)特性的作物系數(shù)Kc,計(jì)算實(shí)際蒸散量ET(TE)。

    TE=TE,0KC。

    (2)

    由于研究區(qū)支柱產(chǎn)業(yè)為棉花種植,因此,根據(jù)新疆當(dāng)?shù)孛藁ǖ纳L周期并考慮不同生育階段,逐月選取相應(yīng)的作物系數(shù)Kc值進(jìn)行反演(表1)。而在冬季非作物生長期,地表被積雪覆蓋,地表實(shí)際蒸散量近似于潛在蒸散量,即TE=TE,0。

    表1 新疆棉花不同生育階段作物系數(shù)取值

    Table 1Kcof cropsrelative to development stages

    生長階段生長周期/d生長時(shí)期作物系數(shù)Kc出苗期303月(春季)0.35苗期504—5月(春季)0.76蕾期606月(夏季)0.97花鈴期557—8月(夏季)1.18成熟期609—10月(秋季)0.60

    數(shù)據(jù)來源于GB/T 32136—2015《農(nóng)業(yè)干旱等級》。

    1.4 空間尺度蒸散量反演

    采用由荷蘭BASTINNASSEN開發(fā)的Surface Energy Balance Algorithm for Land(SEBAL)模型[9],主要借助地表能量的平衡方程計(jì)算區(qū)域的日蒸散量。該模型自1998年起用于實(shí)際研究,現(xiàn)已成為研究蒸散量空間分布特征的重要方法。其基本方程見式(3),凈輻射、感熱通量和土壤熱通量的計(jì)算見式(4)~(7),模型運(yùn)行流程見圖1。

    潛熱通量的計(jì)算公式為

    Eλ=Rn-G-H。

    (3)

    式(3)中,Rn為凈輻射,W·m-2;G為土壤熱通量,W·m-2;Eλ為潛熱通量,W·m-2;H為感熱通量,W·m-2。

    Rn=(1-α)Rswd+εRlwd-εσTs4,

    (4)

    (5)

    (6)

    (7)

    式(4)~(7)中,α為地表反照率;Rswd和Rlwd分別為太陽下行輻射和下行長波輻射,W·m-2;ε為地表比輻射率;σ為Stefan-Bolzmann常數(shù),取值5.67×10-8;Ts為地表輻射溫度,K;INDV為歸一化植被指數(shù);ρ為空氣密度,kg·m-3;cp為空氣定壓比熱容,J·kg-1·K-1;TZ1和TZ2分別為地面Z1與Z2高度處的空氣溫度;ra為空氣動力學(xué)阻抗,s·m-1;u*為摩擦速度,m·s-1;k為Karman常數(shù),取值0.41;Z1和Z2分別為距離地面高度,m。

    圖1 SEBAL模型流程圖Fig.1 Flowchart of the SEBAL model

    2 結(jié)果與分析

    2.1 精河流域?qū)嶋H蒸散量時(shí)間變化特征分析

    2.1.1 實(shí)際蒸散量年際及季節(jié)變化

    新疆精河流域的年實(shí)際蒸散量在各年代的變化如圖2所示。研究區(qū)1953—2015年的近60 a間,年ET值在492.521~1 098.199 mm之間,多年平均值為785.92 mm,1967年出現(xiàn)最大值,1994年出現(xiàn)最小值。近60 a總體呈波動下降趨勢,變化幅度為-38.58 mm·(10 a)-1。

    實(shí)際蒸散量在不同季節(jié)存在較大差異。近60 a中,平均實(shí)際蒸散量以夏季最高,為454.955 mm;春季和秋季次之,分別為209.894和82.403 mm;冬季最小,為38.671 mm,但各季節(jié)實(shí)際蒸散量與年實(shí)際蒸散量變化趨勢基本一致,均呈現(xiàn)較顯著下降趨勢(表2)。其中,實(shí)際蒸散量變化速率最大值出現(xiàn)在夏季,為-21.41 mm·(10 a)-1,最小值出現(xiàn)在冬季,為-1.089 1 mm·(10 a)-1,相差近20倍。這主要由于精河流域的冬季氣候寒冷干燥,使得該季節(jié)對年實(shí)際蒸散量的變化影響甚小。

    圖2 近60 a精河流域年際、季節(jié)實(shí)際蒸散量Fig.2 Variation of the annual and seasonal ETs in the Jinghe Watershed over the past 60 years

    表2 精河流域?qū)嶋H蒸散量的年際、季節(jié)線性變化趨勢統(tǒng)計(jì)

    Table 2 Linear variation trends of annual and seasonal ETs in the Jinghe Watershed

    時(shí)間趨勢系數(shù)趨勢/〔mm·(10a)-1〕擬合趨勢方程全年-3.8580-38.580Y=-3.858X+8441春季-0.9636-9.636Y=-0.9636X+2122夏季-2.1410-21.410Y=-2.141X+4702秋季-0.6451-6.451Y=-0.6451X+1362冬季-0.1089-1.089Y=-0.1089X+254.8

    X為年份;Y為實(shí)際蒸散量。

    2.1.2 實(shí)際蒸散量突變特征分析

    分析年代際距平百分比的變化發(fā)現(xiàn),近60 a精河流域?qū)嶋H蒸散發(fā)表現(xiàn)為3個(gè)階段:20世紀(jì)50年代,研究區(qū)實(shí)際蒸散量的累積距平為負(fù),實(shí)際蒸散量略低于多年平均值;20世紀(jì)60—80年代,研究區(qū)實(shí)際蒸散量的累積距平為正,年代最大實(shí)際蒸散量超出多年平均值約20%,但在經(jīng)歷了20世紀(jì)60年代的極值后迅速下降;20世紀(jì)90年代至今,研究區(qū)實(shí)際蒸散量的累積距平為負(fù)值,且逐步增大,說明實(shí)際蒸散量在緩步提升,但依然小于多年平均值。從年代變化幅度看,蒸散量累積距平經(jīng)歷正、負(fù)值明顯變化的年代為20世紀(jì)60和90年代(圖3)。

    1950s指1953—1959年;2010s指2010—2015年。

    Mann-Kendall(M-K)突變檢驗(yàn)是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn),序列大于0表明參數(shù)呈上升趨勢,反之呈下降趨勢,當(dāng)其超過0.05顯著性水平檢驗(yàn)臨界線(±1.96)時(shí),則上升或下降趨勢顯著。對1953—2015年新疆精河流域年實(shí)際蒸散量序列進(jìn)行M-K檢驗(yàn)(圖4),結(jié)果表明,1953—2015年正序列曲線(UF)總體為持續(xù)下降趨勢,尤其1978年后小于0,即蒸散量1978年后是下降趨勢;1961—1974年UF曲線突破了顯著性水平α=0.05的置信線1.96,蒸散量呈顯著增大趨勢;在1992年突破了置信線-1.96,蒸散量呈顯著性減小趨勢;UF和反序列曲線(UB)于1982年在臨界線區(qū)間內(nèi)有1個(gè)明顯的交點(diǎn),表明1982年發(fā)生突變,突變后精河流域平均年蒸散量較突變前減少150.654 mm,降低幅度為17.361%。

    圖4 精河流域蒸散量Mann-Kendall突變檢驗(yàn)Fig.4 Mann-Kendall tests of the evapotranspiration in the Jinghe Watershed for sudden changes

    對60 a各季節(jié)實(shí)際蒸散量做突變檢驗(yàn),發(fā)現(xiàn)春、夏和冬3季突變點(diǎn)與年際突變點(diǎn)一致,即1982年后3個(gè)季節(jié)的蒸散量均明顯小于之前對應(yīng)各季,蒸散量平均減少幅度為15%左右。秋季蒸散量由高轉(zhuǎn)低的突變則出現(xiàn)在1980年,減少幅度為25%。

    2.1.3 實(shí)際蒸散量變化周期分析

    根據(jù)1953—2012年精河流域?qū)嶋H蒸散量數(shù)據(jù),繪制Morlet小波變換系數(shù)的實(shí)部等值線圖(圖5)。結(jié)果表明,研究區(qū)蒸散量存在6~8、17~20和28~31 a的3類尺度周期變化規(guī)律。其中,6~8 a的周期明顯存在于1953—1970年,17~20 a的周期則僅明顯存在于1957—1973年,而28~31 a的周期明顯始終貫穿于近60 a,但該周期正負(fù)值中心明顯存在由30向26 a下降的趨勢。

    圖5 精河流域蒸散量Morlet小波系數(shù)實(shí)部等值線圖Fig.5 Real part of Morlet wavelet coefficients of the ET in the Jinghe watershed

    進(jìn)一步對其進(jìn)行小波方差計(jì)算,可發(fā)現(xiàn)研究區(qū)實(shí)際蒸散量變化過程中存在的顯著主周期為29 a(圖6)。同時(shí),基于實(shí)際蒸散量在29 a處的小波系數(shù)值〔f(t)〕與年份t建立回歸方程:f(t)=334sin (0.334 5t-52.86),R2=0.93,P<0.05,并選擇2016—2030年為預(yù)測期對流域蒸散量變化趨勢進(jìn)行預(yù)測。

    結(jié)果表明,預(yù)測期2016—2030年整體蒸散量呈波動變化。其中,2016—2020年蒸散量將持續(xù)增加,2020年將發(fā)生突變,轉(zhuǎn)為下降期,蒸散量逐年減少。至2029年,蒸散量將再次進(jìn)入上升周期。

    2.2 精河流域蒸散量空間尺度變化分析

    基于SEBAL模型,選擇1998年9月25日、2007年9月18日和2011年9月13日夏末秋初3期Landsat 5 TM影像數(shù)據(jù),對研究區(qū)蒸散量空間分布特征進(jìn)行分析。由于研究區(qū)沒有實(shí)測日蒸散量數(shù)據(jù),只有1個(gè)國家級氣象觀測站點(diǎn)——精河氣象站,因此,參考劉志武等[10]、王霞[11]和姜紅[12]的方法,對日蒸散量反演精度進(jìn)行驗(yàn)證(表3),反演精度符合要求。

    圖6 精河流域?qū)嶋H蒸散發(fā)量小波方差和小波系數(shù)Fig.6 Wavelet variances and coefficients of the ET in the Jinghe Watershed

    表3 SEBAL模型精度驗(yàn)證結(jié)果

    Table 3 Accuracy evaluation of the SEBAL model with Penman-Monteith formula

    影像時(shí)間FAO-PM蒸散量/mmSEBAL蒸散量/mm相對誤差/%實(shí)測地表溫度/℃遙感反演地表溫度/℃相對誤差/%1998-09-252.6232.381-9.2224.025.2995.412007-09-182.1762.110-3.0320.523.10912.732011-09-131.8421.9194.1820.922.8069.12

    由圖7可見,研究區(qū)極高值出現(xiàn)在艾比湖湖區(qū),這主要是由于湖面水分充足,且湖水比熱容較大,湖區(qū)溫度上升和下降的幅度較小。1998—2007年,湖區(qū)的日蒸散量高值范圍有縮減,2007—2011年湖區(qū)的日蒸散量高值范圍有一定程度恢復(fù),主要是由于區(qū)域氣候暖濕化,精河上游實(shí)施節(jié)水灌溉,2007年艾比湖濕地自然保護(hù)區(qū)晉升為國家級自然保護(hù)區(qū),生態(tài)環(huán)境得到恢復(fù),湖區(qū)面積增加,日蒸散量高值范圍也增加。研究區(qū)中部主要是綠洲農(nóng)業(yè)種植區(qū),有較好的地表作物植被覆蓋,使得日蒸散量分布較均勻,在近20 a保持在一個(gè)相對穩(wěn)定的水平,其蒸散量僅次于湖面日蒸散量,平均值在6.5 mm左右。由于耕地面積的擴(kuò)大,農(nóng)田蒸散區(qū)也隨之?dāng)U大。在研究區(qū)東北部和東南部有小部分山地,因山區(qū)有冰雪和冰川,造成該區(qū)日蒸散量相對較高。低值區(qū)域出現(xiàn)在研究區(qū)北部、西北部和西南部,范圍為0.53~1.27 mm,主要原因是這些區(qū)域植被覆蓋度較低,植物體蒸發(fā)、蒸騰量少。

    在不同土地利用類型下地表土壤的濕度和地表溫度狀況有所不同,會對區(qū)域日蒸散量產(chǎn)生顯著影響,因此,通過參考中國科學(xué)院土地利用覆蓋分類體系,解譯獲取了研究區(qū)1998、2007和2011年土地利用類型圖(圖8)。圖8表明,研究區(qū)的水域(包括艾比湖地面)反射率較低,瞬間獲取的太陽輻射能較多,故實(shí)際蒸散發(fā)在整個(gè)研究區(qū)最大,日蒸散量在7.3~9.32 mm之間。研究區(qū)耕地區(qū)域相對集中,主要種植作物為棉花。1998—2011年精河縣耕地面積不斷增大,耕地邊緣地帶的草地逐漸轉(zhuǎn)變?yōu)楦?使得日蒸散量在5.31~7.55 mm之間的區(qū)域不斷擴(kuò)大。研究區(qū)內(nèi)林地以干旱區(qū)的落葉闊葉林為主,零星分布在耕地的邊緣和河流附近水分比較充足的地帶,由于林地土壤涵養(yǎng)的水分比其他未利用地多,樹木發(fā)達(dá)的根系汲水加之葉面蒸騰,使該區(qū)域蒸散量也較高,日蒸散量在3.95~5.62 mm之間。居民用地在研究區(qū)內(nèi)的面積不大,呈零星分布,較低的地表自然植被覆蓋度以及人工不透水面的建造在較大程度上限制了蒸散發(fā),日蒸散量主要在1.35~3.32 mm之間,呈較低水平。研究區(qū)內(nèi)的未利用地多為裸露巖石、戈壁、荒漠和沙地,由于干旱區(qū)未利用地區(qū)域地下水位較深,地表含水量較少,導(dǎo)致該類型土地蒸散量很小,僅處于0.02~5.69 mm之間。西北部是艾比湖濕地轉(zhuǎn)移出的鹽堿地,含水量高于其他未利用地,其日蒸散量在未利用地區(qū)域中相對較高。

    2.3 精河流域?qū)嶋H蒸散量影響因素分析

    2.3.1 自然因素

    分別對平均降水量、氣溫、相對濕度、日照時(shí)數(shù)、風(fēng)速、太陽輻射、日較差與實(shí)際蒸散量進(jìn)行相關(guān)性分析,其相關(guān)系數(shù)分別為-0.337(P<0.01)、0.692(P<0.05)、-0.533(P<0.01)、0.385(P<0.05)、0.954(P<0.01)、0.652(P<0.01)和0.500(P<0.01)。精河西北部存在著名大風(fēng)口——阿拉山口,使得風(fēng)速成為實(shí)際蒸散發(fā)生的最主要驅(qū)動力。同時(shí),太陽輻射和平均氣溫影響著近地面的熱交換,也成為影響區(qū)域?qū)嶋H蒸散量不可忽視的氣象因素,這也與干旱區(qū)蒸散量相關(guān)研究結(jié)論保持一致。

    圖7 精河流域日蒸散量空間分布圖Fig.7 Spatial distribution of daily evapotranspiration in the Jinghe watershed

    圖8 精河流域土地利用分類圖Fig.8 Classification of land use types in the Jinghe Watershed

    2.3.2 人文因素

    人文因素的作用尤其是土地利用類型對地表實(shí)際蒸散量的影響顯著而強(qiáng)烈。自20世紀(jì)50年代以來,精河流域大規(guī)模開墾后農(nóng)田灌溉面積和灌溉引水量大幅增加是造成精河流域蒸散量緩慢上升的主要人為因素。1968—1994年是精河流域蒸散量變化波動較顯著的時(shí)段。這一時(shí)期,由于開荒面積不斷增長和水資源不合理利用,導(dǎo)致土地鹽堿化和植被退化加劇,20世紀(jì)80年代后大量開墾地撂荒,精河流域蒸散量1968年出現(xiàn)最高值后迅速下降,在1978年出現(xiàn)小幅波動上升后,進(jìn)一步急速下降,蒸散量變化的波動性強(qiáng)烈。而蒸散量在1994年到達(dá)歷史最低點(diǎn)627.31 mm后,開始緩步上升,這與精河入湖徑流量的增加和當(dāng)?shù)乜茖W(xué)引水灌溉、提高水資源利用率有緊密關(guān)系。20世紀(jì)90年代后期,在新疆整體降水減少和氣溫升高的雙重作用下,研究區(qū)徑流總量減少,蒸散量總體偏低。精河流域近60 a的年際實(shí)際蒸散量變化趨勢與我國西部大部分地區(qū)的變化趨勢[13]基本一致。

    3 結(jié)論

    精河流域近60 a的年際實(shí)際蒸散量整體呈現(xiàn)38.58 mm·(10 a)-1的下降趨勢,且以29 a為主周期波動變化,預(yù)測表明,2016—2020年蒸散量將持續(xù)增加,2020年起將突變轉(zhuǎn)為下降期,2029年蒸散量將再次進(jìn)入上升周期。在季節(jié)變化特征上,研究區(qū)實(shí)際蒸散量整體表現(xiàn)為夏季>春季>秋季>冬季的特征。

    研究區(qū)實(shí)際蒸散量受各種因素的綜合影響,而不同因素之間又相互制約。自然因素層面,平均風(fēng)速是影響精河流域?qū)嶋H蒸散量的最重要因子,但太陽輻射和平均氣溫的影響不可忽視。而近60 a來精河流域大規(guī)模開墾導(dǎo)致的農(nóng)田灌溉面積和灌溉引水量大幅增加,成為影響研究區(qū)實(shí)際蒸散量的主要人文因素,不同土地利用類型蒸散量整體表現(xiàn)為水域>耕地>林地>草地>未利用地。

    綜上,對于研究區(qū)而言,農(nóng)業(yè)蒸散耗水是水分消耗的主要部分,恢復(fù)生態(tài)、節(jié)水灌溉、降低農(nóng)田蒸散量、減少作物需水量和灌溉量以及降低地表干燥度是改善研究區(qū)脆弱農(nóng)業(yè)生產(chǎn)及生態(tài)環(huán)境的重點(diǎn),如何平衡生態(tài)環(huán)境與農(nóng)業(yè)經(jīng)濟(jì)的發(fā)展將是未來發(fā)展的關(guān)鍵。

    [1] MCCABE M F,WOOD E F.Scale Influences on the Remote Estimation of Evapotranspiration Using Multiple Satellite Sensors[J].Remote Sensing of Environment,2006,105(4):271-285.

    [2] YAIR A,RAZ-YASSIF N.Hydrological Processes in a Small Arid Catchment:Scale Effects of Rainfall and Slope Length[J].Geomorphology,2004,6(1/2):155-169.

    [3] PAUWEILS V R N,SAMSON R.Comparison of Different Methods to Measure and Model Actual Evapotranspiration Rates for a Wet Sloping Grassland[J].Agricultural Water Management,2006,82(1/2):1-24.

    [4] 魯向暉,白樺,穆興民,等.江西省潛在蒸發(fā)量變化規(guī)律及“蒸發(fā)悖論”成因分析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2016,32(4):552-557.[LU Xiang-hui,BAI Hua,MU Xing-min,etal.Trends of Potential Evapotranspiration in Jiangxi Province and Formation of Evaporation Paradox[J].Journal of Ecology and Rural Environment,2016,32(4):552-557.]

    [5] 張淑蘭,于澎濤,王彥輝,等.涇河上游流域?qū)嶋H蒸散量及其各組分的估算[J].地理學(xué)報(bào),2011,66(3):385-395.[ZHANG Shu-lan,YU Peng-tao,WANG Yan-hui,etal.Estimation of Actual Evapotranspiration and Its Component in the Upstream of Jinghe Basin[J].Acta Geographica Sinica,2011,66(3):385-395.]

    [6] 賈志軍,吳文心,趙海亮.三江平原大豆田蒸散量模擬[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2014,30(1):50-56.[JIA Zhi-jun,WU Wen-xin,ZHAO Hai-liang.Modeling of Evapotranspiration of Soybean Field in the Sanjiang Plain[J].Journal of Ecology and Rural Environment,2014,30(1):50-56.]

    [7] 施云霞,王范霞,毋兆鵬.基于CLUE-S模型的精河流域綠洲土地利用空間格局多情景模擬[J].國土資源遙感,2016,28(2):154-160.[SHI Yun-xia,WANG Fan-xia,WU Zhao-peng.Multi-Simulation of Spatial Distribution of Land Use Based on CLUE-S in Jinghe Watershed[J].Remote Sensing for Land & Resources,2016,28(2):154-160.]

    [8] 楊振常,李玉霖,崔奪,等.半干旱典型沙區(qū)1951—2005年主要?dú)庀笠丶皾撛谡羯⒘康淖兓厔菅芯縖J].中國沙漠,2012,32(5):1384-1392.[YANG Zhen-chang,LI Yu-lin,CUI Duo,etal.Changes of Main Climatic Parameters and Potential Evapotransipiration in Typical Semi-Arid Sandy Lands of Northern China During 1951-2005[J].Journal of Desert Research,2012,32(5):1384-1392.]

    [9] BASTIAANSSEN W G M,MENENTI M,FEDDES R A,etal.A Remote Sensing Surface Energy Balance Algorithm for Land (SEBAL):1.Formulation[J].Journal of Hydrology,1998,212/213:198-212.

    [10]劉志武,雷志棟,黨安榮,等.遙感技術(shù)和SEBAL模型在干旱區(qū)騰發(fā)量估算中的應(yīng)用[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,44(3):421-424.[LIU Zhi-wei,LEI Zhi-dong,DANG An-rong,etal.Remote Sensing and the SEBAL Model for Estimating Evapotranspirationin Arid Regions[J].Journal of Tsinghua University(Science and Technology),2004,44(3):421-424.]

    [11]王霞.基于SEBAL模型的極端干旱氣候區(qū)域蒸散發(fā)的遙感估算[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2012.[WANG Xia.Study on Regional Evapotranspiration Estimation Using Remote Sensing Technology of Extreme Arid Climate Area Based on SEBAL[D].Urumqi:Xinjiang Agricultural University,2012.]

    [12]姜紅.基于MODIS影像的新疆奇臺縣區(qū)域蒸散發(fā)量的研究[D].烏魯木齊:新疆大學(xué),2007.[JIANG Hong.Research of Evapotranspiration in Qitai County,Xinjiang Based on MODIS Data[D].Urumqi:Xinjiang University,2007.]

    [13]李琴,陳曦,劉英,等.干旱區(qū)區(qū)域蒸散發(fā)量遙感反演研究[J].干旱區(qū)資源與環(huán)境,2012,26(8):108-112.[LI Qin,CHEN Xi,LIU Ying,etal.Regional Evapotranspiration Retrieval in Arid Areas[J].Journal of Arid Land Resources and Environment,2012,26(8):108-112.]

    (責(zé)任編輯: 許 素)

    Temporal and Spatial Characteristics of Actual Evapotranspiration in Jinghe Watershed, Xinjiang.

    DAIPeng-chao1,2,NIUSu-juan1,WUZhao-peng1,2,SHIYun-xia1

    (1.College of Geography and Tourism Science, Xinjiang Normal University, Urumqi 830054, China;2.Xinjiang Laboratory of Lake Environment and Resources in Arid Zone, Urumqi 830054, China)

    The Jinghe watershed is an important part of the North Slope Economic Development Belt of the Tianshan Mountains in Xinjiang. The meteorological data of the period from 1953 to 2015 and the Landsat 5 TM images of the years 1998, 2007 and 2011 of the region were analyzed for spatio-temporal variation of the actual evapotranspiration, its characteristics and periodicity with the help of the Penman-Monteith formula, SEBAL model, Morlet wavelet analysis model and M-K mutation test model. Results show: (1) Temporal and spatial distribution of evapotranspiration is influenced by meteorological factors and underlying surface in the Jinghe watershed. The evapotranspiration of the past 60 years as a whole exhibits a significant declining trend with fluctuation in the year 1967 and 1994. In terms of temporal evapotranspiration variation rate, the four seasons follows an order of summer > spring > fall > winter. (2) The mean annual evapotranspiration dropped about 150.654 mm (17.361%) after the sudden change in 1982. The variation of the actual evapotranspiration consisted of three oscillation periods, each lasted for around 29 years, which indicates that in the 15 years to follow the mean annual evapotranspiration tends to ascend first by a small margin and then descend. And (3) spatial distribution of the evapotranspiration is significantly related to the type of land use, and it declines with the type of surface land use follows an order of waters>cultivated land>forest land>grass land>unused land.

    Jinghe Watershed; actual evapotranspiration; Penman-Monteith formula; SEBAL model

    2017-02-05

    國家自然科學(xué)基金(41261055);新疆維吾爾自治區(qū)科技基金(201533109);新疆師范大學(xué)地理學(xué)博士點(diǎn)開放基金(XJNU-DL-201602);新疆師范大學(xué)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練項(xiàng)目(201710762004);新疆師范大學(xué)“十三五”校級重點(diǎn)學(xué)科地理學(xué)

    P426.2

    A

    1673-4831(2017)07-0600-07

    10.11934/j.issn.1673-4831.2017.07.004

    代鵬超(1990—),男,新疆烏魯木齊人,碩士生,研究方向?yàn)橘Y源環(huán)境與地理信息系統(tǒng)。E-mail: daipengchaoxj@sina.com

    ① 通信作者E-mail: wuzhaopengxj@sina.com

    猜你喜歡
    精河流域研究
    壓油溝小流域
    FMS與YBT相關(guān)性的實(shí)證研究
    遼代千人邑研究述論
    堡子溝流域綜合治理
    羅堰小流域
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    水利規(guī)劃與設(shè)計(jì)(2018年1期)2018-01-31 01:53:37
    博爾塔拉河·精河及賽里木湖水質(zhì)比較研究
    地溫、氣溫、氣壓對精河臺水平擺觀測的影響
    中國地震(2015年2期)2015-09-04 01:34:46
    美女国产视频在线观看| 午夜福利在线观看吧| 久久久色成人| 国产一区有黄有色的免费视频 | 国产精品国产三级专区第一集| 精品一区二区免费观看| 国产精品久久久久久av不卡| 内地一区二区视频在线| 色综合色国产| 黄片无遮挡物在线观看| 观看美女的网站| 亚洲精品第二区| 欧美潮喷喷水| 三级国产精品欧美在线观看| 午夜激情久久久久久久| 亚洲av免费高清在线观看| 国产在线一区二区三区精| 高清视频免费观看一区二区 | 国产精品久久久久久久久免| 欧美激情国产日韩精品一区| 国产黄频视频在线观看| 日本爱情动作片www.在线观看| 午夜激情欧美在线| 91精品一卡2卡3卡4卡| 国内少妇人妻偷人精品xxx网站| 亚洲18禁久久av| 成人无遮挡网站| 欧美极品一区二区三区四区| 成人亚洲精品一区在线观看 | 岛国毛片在线播放| www.色视频.com| 男女视频在线观看网站免费| 日日啪夜夜爽| 一个人看的www免费观看视频| 日本色播在线视频| 嫩草影院入口| 精品人妻一区二区三区麻豆| 精品不卡国产一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 午夜免费激情av| 精品久久久久久久久久久久久| 少妇高潮的动态图| 草草在线视频免费看| 久久6这里有精品| 国产单亲对白刺激| 白带黄色成豆腐渣| 精品久久久久久电影网| 久久久久九九精品影院| 日本三级黄在线观看| 国产精品一区二区三区四区免费观看| 欧美bdsm另类| 国产免费视频播放在线视频 | 久久人人爽人人爽人人片va| 国产麻豆成人av免费视频| 成人一区二区视频在线观看| 1000部很黄的大片| 色综合站精品国产| 麻豆av噜噜一区二区三区| 97超碰精品成人国产| 久久热精品热| 少妇裸体淫交视频免费看高清| 非洲黑人性xxxx精品又粗又长| 赤兔流量卡办理| 国内少妇人妻偷人精品xxx网站| 色哟哟·www| av女优亚洲男人天堂| 亚洲人成网站高清观看| 国产成人午夜福利电影在线观看| 在线免费十八禁| 亚洲成人久久爱视频| 欧美高清性xxxxhd video| 亚洲性久久影院| 大话2 男鬼变身卡| 亚洲av男天堂| 亚洲欧美日韩卡通动漫| 男人和女人高潮做爰伦理| 直男gayav资源| 久久久久久久亚洲中文字幕| 乱码一卡2卡4卡精品| 99久久精品一区二区三区| 亚洲成色77777| 久久精品人妻少妇| 在线免费观看的www视频| 国产精品精品国产色婷婷| 青春草亚洲视频在线观看| 亚洲成人久久爱视频| 十八禁国产超污无遮挡网站| 免费看光身美女| 久久精品人妻少妇| 亚洲激情五月婷婷啪啪| 免费看a级黄色片| 精品一区二区免费观看| 成人特级av手机在线观看| 国产黄色小视频在线观看| 国语对白做爰xxxⅹ性视频网站| 久久久色成人| 成年人午夜在线观看视频 | 免费少妇av软件| 黄片无遮挡物在线观看| 99九九线精品视频在线观看视频| 亚洲av电影不卡..在线观看| 亚洲成色77777| 丰满人妻一区二区三区视频av| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 日韩av不卡免费在线播放| 少妇的逼好多水| 一级黄片播放器| 看黄色毛片网站| 1000部很黄的大片| 久久久久久伊人网av| 国产成人精品婷婷| 两个人视频免费观看高清| 免费黄色在线免费观看| 男女那种视频在线观看| 岛国毛片在线播放| 人人妻人人澡欧美一区二区| 777米奇影视久久| 91精品伊人久久大香线蕉| 免费在线观看成人毛片| 免费少妇av软件| 欧美最新免费一区二区三区| 亚洲乱码一区二区免费版| 久久久国产一区二区| 日韩欧美精品免费久久| 如何舔出高潮| 人人妻人人看人人澡| 亚洲欧美一区二区三区黑人 | 丰满人妻一区二区三区视频av| 嫩草影院新地址| 亚洲精品成人av观看孕妇| 欧美97在线视频| 久久久久久久久久人人人人人人| 亚洲国产精品sss在线观看| 一级毛片电影观看| 国产亚洲av嫩草精品影院| 高清毛片免费看| 亚洲性久久影院| 日韩av在线免费看完整版不卡| 97超视频在线观看视频| 欧美激情在线99| 青春草国产在线视频| 日韩一区二区视频免费看| 永久网站在线| 亚洲综合色惰| 91久久精品国产一区二区成人| 国内精品宾馆在线| 狂野欧美白嫩少妇大欣赏| 欧美成人精品欧美一级黄| 成年女人在线观看亚洲视频 | 国产免费福利视频在线观看| 青春草亚洲视频在线观看| 狠狠精品人妻久久久久久综合| 99视频精品全部免费 在线| 亚洲内射少妇av| 日本黄色片子视频| 国内精品美女久久久久久| 亚洲精品国产成人久久av| 亚洲av不卡在线观看| 亚洲国产高清在线一区二区三| 亚洲欧洲日产国产| 亚洲精品乱码久久久v下载方式| 欧美日韩一区二区视频在线观看视频在线 | 人妻夜夜爽99麻豆av| 国产精品爽爽va在线观看网站| 特大巨黑吊av在线直播| 中文天堂在线官网| 欧美激情在线99| 亚洲在线自拍视频| av免费观看日本| 国语对白做爰xxxⅹ性视频网站| 可以在线观看毛片的网站| 街头女战士在线观看网站| 亚洲婷婷狠狠爱综合网| 尾随美女入室| 国产精品久久久久久精品电影| 久久久色成人| 日韩电影二区| 欧美xxⅹ黑人| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 搡老乐熟女国产| 亚洲精品影视一区二区三区av| 精品欧美国产一区二区三| 能在线免费观看的黄片| 国产黄色视频一区二区在线观看| 亚洲av不卡在线观看| 三级经典国产精品| 午夜激情久久久久久久| 91精品国产九色| 国产精品1区2区在线观看.| 亚洲国产欧美在线一区| 成人鲁丝片一二三区免费| h日本视频在线播放| 国产成人91sexporn| 在现免费观看毛片| 99久久精品一区二区三区| av播播在线观看一区| 日韩av在线大香蕉| 成年女人在线观看亚洲视频 | 久久精品国产自在天天线| 精品久久久精品久久久| 一级毛片aaaaaa免费看小| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 亚洲av日韩在线播放| 久久精品综合一区二区三区| 如何舔出高潮| 中文天堂在线官网| 久久韩国三级中文字幕| 日韩一区二区视频免费看| 国产永久视频网站| 国产 亚洲一区二区三区 | 日本免费在线观看一区| 寂寞人妻少妇视频99o| 亚洲精华国产精华液的使用体验| 亚洲无线观看免费| 亚洲色图av天堂| av在线观看视频网站免费| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲av天美| 日本-黄色视频高清免费观看| 日韩国内少妇激情av| 亚洲av成人av| 日韩欧美国产在线观看| 国产亚洲5aaaaa淫片| 能在线免费看毛片的网站| 久久久久免费精品人妻一区二区| 国产av不卡久久| 91精品国产九色| 国产精品综合久久久久久久免费| 国模一区二区三区四区视频| 岛国毛片在线播放| 国产精品爽爽va在线观看网站| 乱码一卡2卡4卡精品| 中文欧美无线码| 亚洲av免费高清在线观看| 日韩在线高清观看一区二区三区| xxx大片免费视频| 99久国产av精品| 国产女主播在线喷水免费视频网站 | av在线老鸭窝| 国产 一区精品| 男人和女人高潮做爰伦理| 国产毛片a区久久久久| 有码 亚洲区| 欧美极品一区二区三区四区| 国产一区二区三区av在线| 成年人午夜在线观看视频 | 国产综合精华液| 大又大粗又爽又黄少妇毛片口| 国产一区二区在线观看日韩| 色视频www国产| 国产av码专区亚洲av| 国产综合精华液| 啦啦啦韩国在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 久久人人爽人人片av| 国产av码专区亚洲av| 免费大片18禁| 免费观看性生交大片5| 国产成人a区在线观看| 国产不卡一卡二| 欧美xxxx性猛交bbbb| 色网站视频免费| 天美传媒精品一区二区| a级一级毛片免费在线观看| 99久久中文字幕三级久久日本| 一级黄片播放器| 肉色欧美久久久久久久蜜桃 | 国产精品伦人一区二区| 2022亚洲国产成人精品| 国产精品一区www在线观看| 最近中文字幕2019免费版| 超碰av人人做人人爽久久| 啦啦啦韩国在线观看视频| 99热网站在线观看| 亚洲精品一区蜜桃| 亚洲美女视频黄频| 晚上一个人看的免费电影| 久久久久久九九精品二区国产| 性色avwww在线观看| 免费大片黄手机在线观看| 狂野欧美激情性xxxx在线观看| 天堂网av新在线| 又大又黄又爽视频免费| 久久精品国产亚洲网站| 91午夜精品亚洲一区二区三区| 亚洲最大成人手机在线| 亚洲欧美日韩东京热| 欧美精品国产亚洲| 2021少妇久久久久久久久久久| 亚洲精品久久午夜乱码| 国产乱来视频区| 久久亚洲国产成人精品v| 2021天堂中文幕一二区在线观| 床上黄色一级片| 成人鲁丝片一二三区免费| 国产一区二区三区综合在线观看 | 日本av手机在线免费观看| 亚洲美女搞黄在线观看| 亚洲成人av在线免费| 美女xxoo啪啪120秒动态图| 久久精品熟女亚洲av麻豆精品 | 女的被弄到高潮叫床怎么办| 97人妻精品一区二区三区麻豆| 麻豆乱淫一区二区| 国产永久视频网站| 最近2019中文字幕mv第一页| 久久久亚洲精品成人影院| 99久国产av精品国产电影| 黄色日韩在线| 欧美成人a在线观看| 少妇熟女aⅴ在线视频| 乱人视频在线观看| 亚洲无线观看免费| 综合色丁香网| 亚洲,欧美,日韩| 最近视频中文字幕2019在线8| 美女被艹到高潮喷水动态| 五月玫瑰六月丁香| 免费观看精品视频网站| 国产男人的电影天堂91| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 久久久久久久久久久丰满| 美女主播在线视频| 国产精品蜜桃在线观看| 黄片wwwwww| 亚洲成人一二三区av| 国产亚洲5aaaaa淫片| 免费黄色在线免费观看| 2022亚洲国产成人精品| 青春草国产在线视频| 亚洲精品亚洲一区二区| 水蜜桃什么品种好| 精品人妻偷拍中文字幕| 狂野欧美白嫩少妇大欣赏| 成人毛片a级毛片在线播放| 美女cb高潮喷水在线观看| 在线a可以看的网站| 久久韩国三级中文字幕| 国产高清不卡午夜福利| 卡戴珊不雅视频在线播放| 97超碰精品成人国产| 一级毛片 在线播放| 久久99热6这里只有精品| 成人亚洲精品一区在线观看 | 久热久热在线精品观看| 欧美三级亚洲精品| or卡值多少钱| 蜜桃亚洲精品一区二区三区| 免费av不卡在线播放| 国产成人福利小说| 欧美成人精品欧美一级黄| 亚洲av日韩在线播放| 搡老乐熟女国产| 国产精品久久久久久久久免| 免费看光身美女| 久久久久久久久大av| 免费观看的影片在线观看| 尾随美女入室| 在线免费观看的www视频| 99热这里只有是精品在线观看| 亚洲成色77777| 成人特级av手机在线观看| 亚洲不卡免费看| 男的添女的下面高潮视频| 少妇裸体淫交视频免费看高清| 午夜久久久久精精品| 午夜福利在线在线| 一区二区三区乱码不卡18| 午夜免费观看性视频| 亚洲欧美一区二区三区黑人 | 国产老妇伦熟女老妇高清| 中文字幕av在线有码专区| 国产av在哪里看| 欧美激情在线99| 秋霞伦理黄片| 免费av毛片视频| 国产69精品久久久久777片| 国产有黄有色有爽视频| 搡老乐熟女国产| 成人综合一区亚洲| 日韩中字成人| 婷婷色综合www| 日韩,欧美,国产一区二区三区| 色哟哟·www| 高清在线视频一区二区三区| 蜜臀久久99精品久久宅男| 国产精品久久久久久av不卡| 观看免费一级毛片| av免费观看日本| 乱系列少妇在线播放| 精品久久久久久久人妻蜜臀av| 黄片无遮挡物在线观看| 一级毛片久久久久久久久女| 亚洲av电影不卡..在线观看| 国产女主播在线喷水免费视频网站 | 老司机影院成人| 免费看不卡的av| 听说在线观看完整版免费高清| 国产午夜精品论理片| 观看美女的网站| 国内精品一区二区在线观看| 亚洲高清免费不卡视频| 亚洲成人av在线免费| 高清欧美精品videossex| 国产精品一区二区性色av| 大片免费播放器 马上看| 人妻一区二区av| 欧美精品国产亚洲| 亚洲不卡免费看| 男女边摸边吃奶| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 国产淫语在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产单亲对白刺激| 日日啪夜夜撸| 老司机影院成人| 欧美三级亚洲精品| 一本久久精品| 成人亚洲精品一区在线观看 | 别揉我奶头 嗯啊视频| 亚洲精品,欧美精品| 免费观看的影片在线观看| 99视频精品全部免费 在线| 亚洲av一区综合| 亚洲丝袜综合中文字幕| 国内揄拍国产精品人妻在线| 国内精品宾馆在线| 免费大片黄手机在线观看| 高清av免费在线| 尤物成人国产欧美一区二区三区| 91aial.com中文字幕在线观看| 成人av在线播放网站| 国产男女超爽视频在线观看| 国产精品精品国产色婷婷| 国产淫语在线视频| 一级毛片 在线播放| 女的被弄到高潮叫床怎么办| 国产午夜精品久久久久久一区二区三区| 久久久久精品性色| 欧美不卡视频在线免费观看| 日韩强制内射视频| 国产伦精品一区二区三区四那| 亚洲av免费在线观看| 晚上一个人看的免费电影| 91av网一区二区| 国产精品国产三级国产专区5o| 男的添女的下面高潮视频| 亚洲va在线va天堂va国产| 人妻一区二区av| 亚洲欧洲日产国产| av在线老鸭窝| 美女黄网站色视频| 国产v大片淫在线免费观看| 日韩一本色道免费dvd| 美女cb高潮喷水在线观看| 免费无遮挡裸体视频| 国产 一区 欧美 日韩| 国产成年人精品一区二区| 国产午夜精品一二区理论片| 亚洲欧美日韩卡通动漫| 国产一级毛片七仙女欲春2| 国产成人aa在线观看| 2021天堂中文幕一二区在线观| 中文字幕免费在线视频6| 久久久久久久久久黄片| 国产亚洲精品av在线| 色播亚洲综合网| 天堂中文最新版在线下载 | 欧美日韩亚洲高清精品| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品日韩av片在线观看| 乱码一卡2卡4卡精品| 国产一区二区亚洲精品在线观看| 免费大片黄手机在线观看| 禁无遮挡网站| 亚洲精品日韩在线中文字幕| 国产91av在线免费观看| 久久久国产一区二区| 国产伦一二天堂av在线观看| 国产精品久久视频播放| 91久久精品电影网| 欧美成人精品欧美一级黄| 美女黄网站色视频| 欧美日本视频| 能在线免费看毛片的网站| 中文字幕av成人在线电影| 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 亚洲国产日韩欧美精品在线观看| 18禁在线无遮挡免费观看视频| 国产免费福利视频在线观看| 成人鲁丝片一二三区免费| 久久久久久伊人网av| 欧美变态另类bdsm刘玥| 日日摸夜夜添夜夜爱| 肉色欧美久久久久久久蜜桃 | 国产精品99久久久久久久久| 欧美日韩国产mv在线观看视频 | 国产免费福利视频在线观看| 亚洲一区高清亚洲精品| 色吧在线观看| 亚洲无线观看免费| 亚洲欧美清纯卡通| 久久99蜜桃精品久久| 国产久久久一区二区三区| 麻豆国产97在线/欧美| 搞女人的毛片| 午夜精品在线福利| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 成人毛片a级毛片在线播放| 老司机影院成人| 三级国产精品欧美在线观看| 在现免费观看毛片| 国内精品美女久久久久久| 午夜福利视频1000在线观看| 岛国毛片在线播放| 欧美激情在线99| 女人久久www免费人成看片| 一区二区三区四区激情视频| 日本熟妇午夜| 搡老妇女老女人老熟妇| 看免费成人av毛片| 欧美zozozo另类| 欧美最新免费一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲经典国产精华液单| 大片免费播放器 马上看| 波野结衣二区三区在线| 国产免费视频播放在线视频 | 丝袜美腿在线中文| 亚洲最大成人中文| 美女cb高潮喷水在线观看| 22中文网久久字幕| 亚洲久久久久久中文字幕| 亚洲国产欧美人成| 国产精品1区2区在线观看.| 最近的中文字幕免费完整| 中文欧美无线码| 日本三级黄在线观看| 爱豆传媒免费全集在线观看| 午夜精品一区二区三区免费看| 亚洲av福利一区| 成人一区二区视频在线观看| a级毛片免费高清观看在线播放| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频 | 国语对白做爰xxxⅹ性视频网站| 一级毛片黄色毛片免费观看视频| 午夜精品在线福利| 中文天堂在线官网| 三级毛片av免费| 婷婷色麻豆天堂久久| 91精品国产九色| 日韩不卡一区二区三区视频在线| 高清在线视频一区二区三区| 中文精品一卡2卡3卡4更新| 午夜视频国产福利| 中国国产av一级| 色5月婷婷丁香| 久久97久久精品| h日本视频在线播放| 毛片女人毛片| 男的添女的下面高潮视频| 欧美xxxx性猛交bbbb| 国产真实伦视频高清在线观看| 免费黄网站久久成人精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最后的刺客免费高清国语| 欧美区成人在线视频| 大香蕉97超碰在线| 日产精品乱码卡一卡2卡三| 亚洲人成网站在线播| 美女国产视频在线观看| 狠狠精品人妻久久久久久综合| 色综合亚洲欧美另类图片| 成年女人在线观看亚洲视频 | 国产不卡一卡二| 日韩欧美国产在线观看| 少妇的逼好多水| 久久精品国产亚洲av涩爱| 国产淫片久久久久久久久| 久久久久精品久久久久真实原创| 国产精品人妻久久久久久| 毛片女人毛片| 日韩精品有码人妻一区| 国产精品不卡视频一区二区| 国产精品嫩草影院av在线观看| 少妇丰满av| 国产亚洲精品av在线| 五月天丁香电影| 日本wwww免费看| 91狼人影院| 简卡轻食公司| 高清日韩中文字幕在线| 99re6热这里在线精品视频| 乱系列少妇在线播放| 一级毛片 在线播放| 成人高潮视频无遮挡免费网站| 最近的中文字幕免费完整| 免费不卡的大黄色大毛片视频在线观看 | 国产精品三级大全| 欧美最新免费一区二区三区| 又爽又黄无遮挡网站| 国产一区二区亚洲精品在线观看| 精品国产露脸久久av麻豆 | 只有这里有精品99| 联通29元200g的流量卡| 欧美日韩国产mv在线观看视频 |