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

    黔中地區(qū)近60年潛在蒸散量時(shí)空變化特征及主導(dǎo)因素識(shí)別

    2021-10-11 05:14:36李洪廣梁任剛楊大方
    水土保持研究 2021年6期
    關(guān)鍵詞:日照時(shí)數(shù)站點(diǎn)風(fēng)速

    肖 楊,周 旭,羅 雪,李洪廣,梁任剛,楊大方

    (貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院,貴陽 550025)

    蒸散量作為區(qū)域能量平衡和水分循環(huán)的重要控制因素,是評(píng)估氣候變化、水資源配置、作物需水量、干旱預(yù)測(cè)與預(yù)警、生態(tài)環(huán)境保護(hù)等領(lǐng)域的關(guān)鍵指標(biāo)[1-2]。在全球氣候變化的背景下,不同區(qū)域潛在蒸散量變化存在差異性,其影響潛在蒸散量變化原因不同,逐漸成為區(qū)域水文氣候變化研究的重點(diǎn)與熱點(diǎn)[3-4]。全球氣候變暖已成為不爭(zhēng)的事實(shí),氣溫上升會(huì)加速水文循環(huán),導(dǎo)致潛在蒸散量增加[5-6]。但是,有些地區(qū)潛在蒸散量呈減少趨勢(shì),被稱為“蒸發(fā)悖論”現(xiàn)象[7-8]。如前蘇聯(lián)[9]、美國(guó)[10]、新西蘭[11]、泰國(guó)等[12]。然而,在其他地區(qū)潛在蒸散量呈增加趨勢(shì),如伊朗[13]、北美東北部[14]、非洲尼羅河沿岸國(guó)家等[15]。在中國(guó)長(zhǎng)江上游[16]、東北北部[17]的年潛在蒸散量顯著增加,在淮河流域[18]、長(zhǎng)江中下游地區(qū)[19]和珠江流域[20]的年潛在蒸散量顯著減少。同時(shí),上述研究確定了影響潛在蒸散量變化的主要?dú)夂蛞蛩?,由于區(qū)域差異性,研究結(jié)果存在異質(zhì)性。一些學(xué)者將潛在蒸散量的變化歸因于日照時(shí)數(shù)或太陽輻射[10-11,13],另外學(xué)者則認(rèn)為風(fēng)速是導(dǎo)致下降趨勢(shì)的主要因素[18-19]。此外,氣候因素之間的相互影響也是促成潛在蒸散量變化的重要原因。因此,在區(qū)域尺度上研究潛在蒸散量的時(shí)空變化趨勢(shì)及其驅(qū)動(dòng)因素,對(duì)厘清區(qū)域與氣候變化的響應(yīng)具有重要意義[21]。

    黔中地區(qū)水資源豐富,但時(shí)空分布不均,且該區(qū)域?yàn)榭λ固厣絽^(qū),自然基底的異質(zhì)性和生態(tài)環(huán)境的脆弱性,加上二元三維水文地質(zhì)結(jié)構(gòu)影響[22],水資源不易儲(chǔ)存于地表,導(dǎo)致水資源稀缺問題突出。受氣候變暖影響,加之農(nóng)業(yè)基礎(chǔ)條件較為薄弱,出現(xiàn)水旱災(zāi)害較多[23]。近年來,工業(yè)經(jīng)濟(jì)發(fā)展快速,對(duì)水資源的需求量較大,水量平衡各項(xiàng)的變化對(duì)地區(qū)工農(nóng)業(yè)發(fā)展及水文過程產(chǎn)生巨大影響,進(jìn)而對(duì)地區(qū)水安全和糧食安全帶來嚴(yán)峻挑戰(zhàn)。潛在蒸散量的研究可分析黔中地區(qū)可獲得水資源量及水分消耗過程,且潛在蒸散量的變化可為各種作物需水量的計(jì)算提供科學(xué)依據(jù),進(jìn)而制定合理的灌溉制度,達(dá)到節(jié)水目的?;诖?,本文利用1960—2019年黔中地區(qū)日氣象數(shù)據(jù),采用FAO推薦的Penman-Monteith方法,計(jì)算黔中地區(qū)日潛在蒸散量,基于Theil-Sen median趨勢(shì)分析與Mann-Kendall檢驗(yàn)、小波分析等方法分析年均及四季潛在蒸散量時(shí)空變化特征與周期變化特征,通過相關(guān)性分析研究黔中地區(qū)潛在蒸散量和氣象因素的變化趨勢(shì)并識(shí)別影響潛在蒸散量變化的主導(dǎo)因素??蔀樵摰貐^(qū)水資源優(yōu)化配置、農(nóng)業(yè)需水預(yù)測(cè)、制定合理灌溉方案和農(nóng)作物種植結(jié)構(gòu)調(diào)整提供科學(xué)參考。

    1 研究區(qū)概況與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    黔中地區(qū)位于貴州省中部地區(qū),包括貴陽、遵義、安順、畢節(jié)、黔東南和黔南6個(gè)市(州),涉及貴安新區(qū)在內(nèi)的33個(gè)縣(市、區(qū)),總面積53 800 km2,占貴州省面積的31%[24]。地勢(shì)自西向東傾斜,屬亞熱帶高原季風(fēng)性氣候,大部分地區(qū)為喀斯特丘陵地貌,年均溫約為15.04℃,年降水量約為1 163.79 mm,年日照時(shí)數(shù)1 175.62 h。土壤類型以黃壤、石灰土、水稻土、沼澤土等為主[25]。區(qū)域森林覆蓋率高,以常綠闊葉林為主,自然資源豐富,是適宜居住的綠色生態(tài)聚居地[24]。

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

    本文采用的氣象數(shù)據(jù)來自中國(guó)氣象數(shù)據(jù)網(wǎng)(http:∥data.com.cn),根據(jù)數(shù)據(jù)的連續(xù)性和時(shí)序性盡可能長(zhǎng)的標(biāo)準(zhǔn),選用黔中地區(qū)11個(gè)氣象站點(diǎn)(圖1),時(shí)間序列為1960—2019年,數(shù)據(jù)包括日平均氣溫、日平均最高氣溫、日平均最低氣溫、日平均相對(duì)濕度、日平均風(fēng)速、日平均日照時(shí)數(shù)等。計(jì)算各站點(diǎn)的日潛在蒸散量,然后進(jìn)行月、季、年潛在蒸散量的統(tǒng)計(jì),季節(jié)的劃分為3—5月為春季,6—8月為夏季,9—11月為秋季,12—2月為冬季[26]。氣象資料空間插值法在通過多種插值法進(jìn)行準(zhǔn)確度分析后采用反距離權(quán)重法[26]。

    圖1 黔中地區(qū)地理位置及氣象站點(diǎn)位置

    2 研究方法

    2.1 Penman-Monteith模型

    計(jì)算潛在蒸散量的方法很多,應(yīng)用最廣泛的是聯(lián)合國(guó)糧農(nóng)組織(FAO)1998年推薦的Penman-Monteith模型[5-6,10],該模型基于能量平衡和水汽擴(kuò)散理論,考慮作物的生理特性和空氣動(dòng)力學(xué)參數(shù)的變化,具有堅(jiān)實(shí)的理論基礎(chǔ),物理意義明確,在各種氣候和時(shí)間步長(zhǎng)計(jì)算中證明效果良好,無需任何局部校準(zhǔn)[8,10]。

    2.2 線性回歸分析

    假設(shè)公式y(tǒng)=ax+b中,y表示特定時(shí)間序列下氣候要素的變化趨勢(shì)[27],x為年序,a表示一元線性方程的斜率,即為氣象要素的線性變化趨勢(shì),若a為正時(shí),則表示氣候要素呈增加趨勢(shì),a為負(fù)時(shí),表示氣候要素呈減小趨勢(shì)。

    2.3 Theil-Sen Median趨勢(shì)分析和Mann-Kendall檢驗(yàn)

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

    (1)

    式中:10表示“上升或增加趨勢(shì)”,β<0表示“下降或減少趨勢(shì)”,β=0表示“趨勢(shì)不變”。

    Mann-Kendall檢驗(yàn)是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn),不需要樣本服從一定的分布,同時(shí)也不受少數(shù)異常值的干擾,被廣泛應(yīng)用于分析具有水文和氣象的時(shí)間序列分析中,用于判斷時(shí)間序列數(shù)據(jù)是否具有上升或下降的趨勢(shì),計(jì)算公式詳見文獻(xiàn)[26]。

    當(dāng)Mann-Kendall進(jìn)一步用于檢驗(yàn)序列突變時(shí),檢驗(yàn)統(tǒng)計(jì)量與Z值有所差異,通過構(gòu)造一秩序列再定義統(tǒng)計(jì)變量。通過分析統(tǒng)計(jì)序列UBk,UFk可以進(jìn)一步分析x的趨勢(shì)變化,判斷明確突變時(shí)間。若UFk>0,則表明序列呈上升趨勢(shì),若UFk<0,則表明序列呈下降趨勢(shì),超過臨界值曲線,則表示上升或下降趨勢(shì)顯著[26]。若UBk,UFk兩曲線出現(xiàn)交點(diǎn),且處于臨界值曲線間,則交點(diǎn)對(duì)應(yīng)的時(shí)間點(diǎn)就是突變時(shí)間。以上公式通過IDL語言和Matlab程序代碼實(shí)現(xiàn)。

    2.4 小波分析

    Morlet小波分析方法能夠清晰地揭示時(shí)間序列中的多種變化,充分反映不同時(shí)間尺度下水文氣象數(shù)據(jù)的變化趨勢(shì)[28]。與多時(shí)間尺度分析方法如移動(dòng)濾波分析和傅里葉分析相比,小波分析在時(shí)域和頻域上都具有更好的局部化功能,并且可以對(duì)時(shí)間序列進(jìn)行定位分析。本文采用此方法對(duì)研究區(qū)進(jìn)行周期分析。

    3 結(jié)果與分析

    3.1 潛在蒸散量時(shí)空變化

    3.1.1 潛在蒸散量時(shí)間變化

    (1)潛在蒸散量年際變化。1960—2019年黔中地區(qū)多年平均潛在蒸散量為1 004.12 mm,最高值出現(xiàn)在1963年,最低值在2012年,波動(dòng)范圍為915.22~1 139.15 mm,極差為223.93 mm。近60 a來年ET0的波動(dòng)幅度較大,總體呈現(xiàn)波動(dòng)下降趨勢(shì)(圖2),線性擬合表明其遞減傾向率約為-6.00 mm/10 a。

    圖2 黔中地區(qū)潛在蒸散量年際變化

    (2)潛在蒸散量年代際變化。年代際尺度上年均和季節(jié)潛在蒸散量變化見表1,在10 a尺度上,黔中地區(qū)潛在蒸散量呈現(xiàn)出“減—減—減—增—減”的波動(dòng)變化狀態(tài)趨勢(shì),其中1970s—1990s一直呈減小趨勢(shì);在2000s潛在蒸散量出現(xiàn)增加趨勢(shì),但在2010s又出現(xiàn)減少趨勢(shì)(圖2),2000s潛在蒸散量的變化幅度為2.8%。從各年代的各不同季節(jié)上看,秋季和冬季1970s—1980s呈減小趨勢(shì),1990s—2000s呈增加趨勢(shì),2010s呈減少趨勢(shì)。2000s年代4個(gè)季節(jié)都呈增加趨勢(shì),其中春季上升幅度最大,變化幅度為4.9%。因此,在2000s年代潛在蒸散量呈現(xiàn)明顯的增加趨勢(shì)。春季和夏季1970s—1990s呈減小趨勢(shì),2000s年代呈增加趨勢(shì)。1970s年代春季減小幅度最大,變化幅度為-8.2%;1990s年代夏季減小幅度較大,變化幅度為-4.0%;1980s和2010s年代秋季減小幅度較大,變化幅度為-4.4%;1980s年代冬季減少幅度最大,變化幅度為-8.0%。從變化幅度上看,四季均呈負(fù)值,春季減少趨勢(shì)最大。從年平均潛在蒸散量來看,整體上呈減小趨勢(shì),變化幅度為-0.5%,比上一年代表現(xiàn)下降趨勢(shì)的年代是1970s,1980s和1990s,其中80年代的變化幅度最大為-2.5%??傮w得出,年均潛在蒸散量下降主要原因是由春季和冬季潛在蒸散量的下降導(dǎo)致的。

    表1 黔中地區(qū)年及季節(jié)潛在蒸散量的年代際變化

    (3)潛在蒸散量季節(jié)變化。圖3是黔中地區(qū)潛在蒸散量的季節(jié)變化,四季整體上呈現(xiàn)減小趨勢(shì),下降速率分別為-2.10 mm/10 a,-3.01 mm/10 a,-0.17 mm/10 a,-0.62 mm/10 a。夏季平均潛在蒸散量最大為373.07 mm,約占全年平均的37.16%,其次是春季和秋季,分別為288.17 mm和218.34 mm,占全年平均ET0的28.75%和21.75%,冬季平均ET0為124.01 mm,是四季中最小值,約占年平均ET0的12.35%。夏春兩季潛在蒸散量占全年的65.91%,對(duì)全年的潛在蒸散量貢獻(xiàn)最大。4個(gè)季節(jié)的R2都較小,說明下降趨勢(shì)受到干擾嚴(yán)重,相比較其他季節(jié),夏季的下降趨勢(shì)略好,但都不夠顯著,季節(jié)變化在年周期上更多地是一種振蕩。

    圖3 黔中地區(qū)潛在蒸散量季節(jié)變化

    3.1.2 潛在蒸散量空間變化

    (1)年際潛在蒸散量空間變化。圖4A顯示了黔中地區(qū)ET0的空間分布情況,除南部地區(qū)的貴陽站、惠水站和凱里站和北部地區(qū)的仁懷站ET0值較高外,西部地區(qū)畢節(jié)站ET0處于最低值,總體表現(xiàn)為自南向北遞減趨勢(shì)。南部的貴陽、惠水、凱里ET0在961.20~1 223.10 mm,西部的畢節(jié)ET0為930.50~1 047.00 mm。

    圖4B顯示了黔中地區(qū)變化趨勢(shì)的空間分布格局,地區(qū)中部(東南—西北走向)呈下降趨勢(shì)、地區(qū)的南北部分呈現(xiàn)增加趨勢(shì)。地區(qū)東北部的匯川與西南部的安順、惠水呈上升趨勢(shì),其余的站點(diǎn)呈下降趨勢(shì)。ET0上升的站點(diǎn)僅有3個(gè),其中一個(gè)站點(diǎn)達(dá)到90%的顯著性水平,其余2個(gè)站點(diǎn)未達(dá)到顯著性水平,但呈小幅度上升趨勢(shì)。其余的8個(gè)站點(diǎn)呈下降趨勢(shì),6個(gè)站點(diǎn)通過置信度90%的顯著性檢驗(yàn),2個(gè)站點(diǎn)未通過顯著性檢驗(yàn),呈不顯著減少變化趨勢(shì)。貴陽站是下降趨勢(shì)最大的站點(diǎn)(-16.20 mm/10 a),匯川站是上升趨勢(shì)最大的站點(diǎn)(6.12 mm/10 a)。

    圖4 黔中地區(qū)潛在蒸散量與M-K趨勢(shì)Z值變化趨勢(shì)空間分布

    (2)不同季節(jié)潛在蒸散量空間變化。春季ET0空間上呈由南向北遞減趨勢(shì),匯川站與畢節(jié)站明顯低于其他地區(qū),處于研究區(qū)最低值狀態(tài)(圖5)。其變化趨勢(shì)呈現(xiàn)東北部和西南部增加、由西北部和東南部向中部遞減趨勢(shì),3個(gè)站點(diǎn)(匯川、安順和惠水)通過90%的置信度檢驗(yàn),呈不顯著增加趨勢(shì),3個(gè)站點(diǎn)(貴陽、息烽和黔西)呈顯著下降趨勢(shì),其余站點(diǎn)均為不顯著下降趨勢(shì)。夏季ET0空間上除北部仁懷ET0高于其他地區(qū),呈由東北部向西南部遞減趨勢(shì)。夏季和秋季ET0變化趨勢(shì)基本呈現(xiàn)東北部和西南部增加、由東南部向西北部遞減趨勢(shì)。夏季ET0有3個(gè)站點(diǎn)為不顯著上升趨勢(shì),中部4站點(diǎn)通過90%的置信度檢驗(yàn),呈顯著下降趨勢(shì)。秋季ET0空間上呈由東南部向西北部遞減趨勢(shì)。秋季ET0有1個(gè)站點(diǎn)(匯川)通過90%的置信度檢驗(yàn),呈顯著上升趨勢(shì),其余8個(gè)站點(diǎn)未通過90%顯著性檢驗(yàn)。冬季ET0空間上呈由南部向北部遞減趨勢(shì)。冬季ET0有1個(gè)站點(diǎn)(匯川)通過90%的置信度檢驗(yàn),呈不顯著上升趨勢(shì),另1個(gè)站點(diǎn)(貴陽)通過90%顯著性檢驗(yàn),呈顯著減少趨勢(shì)。

    圖5 黔中地區(qū)四季潛在蒸散量與M-K趨勢(shì)Z值變化趨勢(shì)空間分布

    3.2 潛在蒸散量突變分析

    運(yùn)用Mann-Kendall檢驗(yàn)法分析黔中地區(qū)近60 a潛在蒸散量的突變特征,從圖6可以看出,UF與UB的交點(diǎn)位于1969年前后,這意味著黔中地區(qū)ET0突變的時(shí)間約為1969年,在1969年以前有1963年和1966年位于零分界線以上,說明該時(shí)段內(nèi)ET0處于上升趨勢(shì),但在1969年以后,黔中地區(qū)UF統(tǒng)計(jì)值小于0,表明1969年以后ET0整體處于下降趨勢(shì),在1984年ET0下降趨勢(shì)達(dá)到0.05顯著性水平。對(duì)1969年前后年潛在蒸散量的計(jì)算,發(fā)現(xiàn)1970—2019年平均潛在蒸散量(996.80 mm)比1960—1969年(1 040.71 mm)下降了43.91 mm。

    圖6 黔中地區(qū)年潛在蒸散量Mann-Kendall檢驗(yàn)

    3.3 潛在蒸散量周期性變化

    為了分析黔中地區(qū)ET0的變化周期,利用Morlet小波分析對(duì)1960—2019年的年ET0進(jìn)行變換,得到小波系數(shù)實(shí)部等值線圖(圖7A)。圖中實(shí)線表示為取正值的小波系數(shù)等值線,即為潛在蒸散量偏高期,虛線表示小波系數(shù)取負(fù)值的等值線,即為潛在蒸散量偏低期。研究區(qū)ET0存在3 a,9 a,21 a左右的周期變化,其中3 a,21 a周期變化貫穿于1960—2019年,而9 a周期變化主要存在于1960—1986年。從小波系數(shù)符號(hào)正負(fù)相位交替變化過程分析,表明ET0在于1960—2019年經(jīng)歷了高—低—高—低—高—低的循環(huán)交替過程,1960—1970年、1982—1993年、2002—2011年為ET0偏高期,其余時(shí)間段為ET0偏低區(qū),研究結(jié)果與ET0的年際變化結(jié)果相吻合。

    通過小波方差來鑒定對(duì)ET0變化規(guī)律具有影響力的時(shí)間尺度(圖7B),小波方差存在兩個(gè)明顯的波峰,分別是3 a和21 a的時(shí)間尺度上,且在21 a的震蕩周期最強(qiáng)烈,取得最大值,為強(qiáng)顯著周期,第二周期為3 a。

    圖7 黔中地區(qū)潛在蒸散量小波分析與小波方差

    3.4 潛在蒸散量變化的氣象因素識(shí)別

    通過分析可得出潛在蒸散量與各氣象因素的相關(guān)性程度(表2),季節(jié)潛在蒸散量、年潛在蒸散量與各氣象因素的相關(guān)性基本一致。從年尺度上分析,平均氣溫、最高氣溫、最低氣溫、日照時(shí)數(shù)、風(fēng)速和水汽壓與ET0呈正相關(guān)關(guān)系,這些氣象因素的增加會(huì)導(dǎo)致ET0的增加,且平均氣溫、最高氣溫、日照時(shí)數(shù)、風(fēng)速和水汽壓與ET0的相關(guān)性通過了0.01水平檢驗(yàn);其中相關(guān)性系數(shù)最大的是日照時(shí)數(shù),然后是風(fēng)速。相對(duì)濕度和降水量與ET0呈負(fù)相關(guān)關(guān)系,這些氣象因素的增加會(huì)導(dǎo)致ET0的減少;這兩個(gè)氣象因素與ET0都在0.01水平上顯著相關(guān)。從氣象因素的年變化趨勢(shì)來看(圖8),風(fēng)速、相對(duì)濕度、降水量和日照時(shí)數(shù)均呈下降趨勢(shì),其氣候變化率分別為-0.03 (m/s)/10 a,-0.15%/10 a,-5.97 mm/10 a,-46.44 h/10 a,其中日照時(shí)數(shù)的變化對(duì)ET0的下降趨勢(shì)影響最大,其次是風(fēng)速。雖然平均氣溫呈上升趨勢(shì),變化率為0.09℃/10 a,其在時(shí)段內(nèi)變化相比其他氣象因素而言不顯著。因此,在年尺度上,日照時(shí)數(shù)和風(fēng)速是影響ET0發(fā)生變化的主導(dǎo)氣象因素。

    圖8 黔中地區(qū)氣象因素變化趨勢(shì)

    表2 黔中地區(qū)年及季節(jié)潛在蒸散量與氣象因素相關(guān)性分析

    季節(jié)尺度上具有與年特征相似性,平均氣溫、最高氣溫、日照時(shí)數(shù)和風(fēng)速在四季上都與ET0呈現(xiàn)正相關(guān)關(guān)系,相反的,相對(duì)濕度、降水量都和ET0呈負(fù)相關(guān)關(guān)系。四季中,日照時(shí)數(shù)與ET0的相關(guān)性最高,均達(dá)到了0.01的正顯著性水平。春季,潛在蒸散量呈下降趨勢(shì)。平均氣溫、最高氣溫、日照時(shí)數(shù)和風(fēng)速與ET0呈顯著正相關(guān)關(guān)系,而最低氣溫、相對(duì)濕度、水汽壓和降水量呈負(fù)相關(guān)關(guān)系。其中呈正相關(guān)性最大的日照時(shí)數(shù)和風(fēng)速呈下降趨勢(shì),變化率分別為-11.60 h/10 a,-0.07 (m/s)/10 a,雖然氣溫呈上升趨勢(shì),降水量處于下降趨勢(shì),為-6.56 mm/10 a,也沒有改變ET0的下降。因此,導(dǎo)致春季ET0下降的主要原因是日照時(shí)數(shù)、風(fēng)速,其次才是降水量。夏季,日照時(shí)數(shù)、風(fēng)速與ET0呈正相關(guān)關(guān)系,相對(duì)濕度、降水量與ET0呈負(fù)相關(guān)關(guān)系。日照時(shí)數(shù)變化率為-20.47 h/10 a,其變化率較大,日照時(shí)數(shù)與風(fēng)速的顯著減少是夏季ET0下降的原因。秋季,平均氣溫、最高氣溫、日照時(shí)數(shù)與風(fēng)速與ET0呈正相關(guān)關(guān)系,其中,平均氣溫、最高氣溫呈上升趨勢(shì),其余呈下降趨勢(shì),日照時(shí)數(shù)與降水量下降趨勢(shì)較大,分別為-7.55 h/10 a,-5.45 mm/10 a,這兩個(gè)因素是影響ET0下降的主要因子。冬季,日照時(shí)數(shù)和風(fēng)速與ET0呈正相關(guān)關(guān)系,均達(dá)到顯著性水平,且都是呈下降趨勢(shì)。相對(duì)濕度、水汽壓和降水量與ET0呈負(fù)相關(guān)關(guān)系,冬季氣溫較低,且日照時(shí)數(shù)較少,整個(gè)地區(qū)蒸發(fā)微弱,因此,這是導(dǎo)致ET0減少的最大原因。

    綜上所述,在年尺度和季節(jié)尺度上各個(gè)氣象因素與ET0的變化分析均存在差異,是氣象因素共同作用的結(jié)果。黔中地區(qū)年ET0的減小的主要?dú)庀笠蛩厥侨照諘r(shí)數(shù)的減少,其次為風(fēng)速。

    4 討 論

    潛在蒸散量作為一項(xiàng)重要的農(nóng)業(yè)與水文變量,對(duì)區(qū)域灌溉排水工程設(shè)計(jì)、流域規(guī)劃具有指導(dǎo)意義[4]??諝鈩?dòng)力學(xué)、輻射、下墊面條件等對(duì)潛在蒸散量產(chǎn)生不同程度影響[29],據(jù)已有研究表明[17],潛在蒸散量受氣溫、降水、風(fēng)速、相對(duì)濕度、日照時(shí)數(shù)等綜合氣象因素直接影響。氣溫的升高與降水量的減少可能會(huì)導(dǎo)致不同程度的干旱,近年來西南地區(qū)的干旱較為頻繁[26],特別是2009年的西南大旱,對(duì)整個(gè)黔中地區(qū)的農(nóng)業(yè)、人口等經(jīng)濟(jì)造成巨大損失。全球變暖和厄爾尼諾現(xiàn)象的加劇導(dǎo)致海面溫度異常,降水少,雨季提前結(jié)束,最終導(dǎo)致氣溫高,降水更少[30]。過去50 a的全球氣溫變化率為0.13℃/10 a,黔中地區(qū)的氣溫增長(zhǎng)率低于全球氣溫上升速率,為0.09℃/10 a。黔中地區(qū)的日照時(shí)數(shù)、風(fēng)速和相對(duì)濕度減少,這與中國(guó)大部分地區(qū)變化相一致[8,31],但是區(qū)域空間下存在異質(zhì)性。氣候特征的變化將改變或已改變黔中地區(qū)的水文循環(huán),這將導(dǎo)致更頻繁的干旱和洪水,水資源的合理配置對(duì)整個(gè)地區(qū)的農(nóng)業(yè)發(fā)展、生態(tài)環(huán)境保護(hù)和經(jīng)濟(jì)協(xié)調(diào)發(fā)展等產(chǎn)生積極的意義。

    黔中地區(qū)近60 a的平均氣溫顯著上升,但是潛在蒸散量顯著減少,說明黔中地區(qū)存在“蒸發(fā)悖論”,這與學(xué)者在研究西南地區(qū)潛在蒸散量的變化分析結(jié)果相同[26]。黔中地區(qū)氣溫升高和濕度降低可能會(huì)導(dǎo)致ET0的增加,但其影響ET0變化的主要?dú)庀笠蛩厥侨照諘r(shí)數(shù)和風(fēng)速,最終導(dǎo)致黔中地區(qū)的“蒸發(fā)悖論”。但是,除氣象因素的影響外,潛在蒸散量在不同下墊面環(huán)境條件下的“蒸發(fā)悖論”是否存在差異,其研究結(jié)果將影響主導(dǎo)因素的識(shí)別,是下一步開展研究的重點(diǎn),研究結(jié)果可為黔中地區(qū)的灌溉管理、水文過程和抗旱減災(zāi)提供科學(xué)依據(jù),

    5 結(jié) 論

    (1)1960—2019年黔中地區(qū)潛在蒸散量呈顯著下降趨勢(shì),其遞減傾向率約為-6.00 mm/10 a。四季潛在蒸散量整體上呈減小趨勢(shì),夏春兩季潛在蒸散量占全年的65.91%,是影響年潛在蒸散量下降的主要原因。

    (2)黔中地區(qū)潛在蒸散量空間分布特征呈自南向北遞減趨勢(shì),其變化趨勢(shì)在地區(qū)中部(東南—西北走向)呈下降趨勢(shì)、在地區(qū)的南北部分呈現(xiàn)增加趨勢(shì)。

    (3)黔中地區(qū)潛在蒸散量在1969年左右發(fā)生一次明顯的減少突變,1970—2019年平均潛在蒸散量(996.80 mm)比1960—1969年(1 040.71 mm)下降了43.91 mm。從小波分析結(jié)果可得,研究區(qū)潛在蒸散量存在3 a,9 a,21 a的左右的周期變化,其中3 a,21 a周期變化貫穿于1960—2019年,而9 a周期變化主要在1960—1986年,且在21 a的震蕩周期最強(qiáng)烈為主周期。

    (4)黔中地區(qū)平均氣溫、最高、最低氣溫呈顯著上升趨勢(shì),而平均風(fēng)速、相對(duì)濕度、降水量與日照時(shí)數(shù)呈顯著下降趨勢(shì),研究區(qū)存在“蒸發(fā)悖論”現(xiàn)象。其中日照時(shí)數(shù)是影響潛在蒸散量變化的主要?dú)庀笠蛩?,其次是風(fēng)速,且各氣象因素在季節(jié)上對(duì)潛在蒸散量影響存在差異性。

    猜你喜歡
    日照時(shí)數(shù)站點(diǎn)風(fēng)速
    福州市近70年日照變化趨勢(shì)分析
    西昌近60年日照時(shí)數(shù)的變化特征分析
    1961~2020年曲麻萊縣日照時(shí)數(shù)變化特征
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
    基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
    基于Web站點(diǎn)的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    2017~2018年冬季西北地區(qū)某站點(diǎn)流感流行特征分析
    1980年~2017年大冶市日照時(shí)數(shù)變化特征分析
    首屆歐洲自行車共享站點(diǎn)協(xié)商會(huì)召開
    怕被人認(rèn)出
    2021天堂中文幕一二区在线观| 中文字幕av成人在线电影| www.熟女人妻精品国产| 亚洲欧美激情综合另类| 日本免费a在线| 草草在线视频免费看| 又粗又爽又猛毛片免费看| 又紧又爽又黄一区二区| 国产亚洲精品综合一区在线观看| 男女那种视频在线观看| 中出人妻视频一区二区| 能在线免费观看的黄片| 一本久久中文字幕| 免费大片18禁| 琪琪午夜伦伦电影理论片6080| 无遮挡黄片免费观看| 亚洲成人免费电影在线观看| av天堂中文字幕网| 午夜日韩欧美国产| 天堂网av新在线| 别揉我奶头 嗯啊视频| 久久久久性生活片| av女优亚洲男人天堂| 18禁裸乳无遮挡免费网站照片| 欧美一区二区亚洲| 中文字幕精品亚洲无线码一区| 日韩人妻高清精品专区| 又紧又爽又黄一区二区| 国产乱人伦免费视频| 内射极品少妇av片p| 亚洲三级黄色毛片| 国内精品久久久久精免费| 精品人妻一区二区三区麻豆 | 亚洲第一欧美日韩一区二区三区| 国产美女午夜福利| 高清毛片免费观看视频网站| 久久久久久久精品吃奶| 国产精品久久视频播放| 亚洲精品影视一区二区三区av| 色综合站精品国产| 久久精品影院6| 久久精品国产亚洲av香蕉五月| 热99在线观看视频| 男人舔女人下体高潮全视频| 精品人妻1区二区| 好男人电影高清在线观看| 成人av一区二区三区在线看| 悠悠久久av| 我的老师免费观看完整版| 久久精品91蜜桃| 日本免费一区二区三区高清不卡| 亚洲一区二区三区色噜噜| 国产一级毛片七仙女欲春2| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一区二区三区四区免费观看 | 亚洲美女搞黄在线观看 | 97碰自拍视频| 中文字幕av成人在线电影| 久久久久久九九精品二区国产| 精品人妻偷拍中文字幕| 999久久久精品免费观看国产| 欧美三级亚洲精品| 欧美最新免费一区二区三区 | 色在线成人网| 国产精品久久久久久亚洲av鲁大| 大型黄色视频在线免费观看| 波多野结衣高清无吗| 成人av一区二区三区在线看| 国产爱豆传媒在线观看| 精品人妻视频免费看| 国产男靠女视频免费网站| 一区二区三区激情视频| 日本黄大片高清| 身体一侧抽搐| 身体一侧抽搐| 国产精品亚洲一级av第二区| 精品一区二区三区人妻视频| 亚洲精品在线美女| 国产精品1区2区在线观看.| 网址你懂的国产日韩在线| 国产精品98久久久久久宅男小说| 成人高潮视频无遮挡免费网站| 搞女人的毛片| 成人国产综合亚洲| 欧美激情久久久久久爽电影| 黄色一级大片看看| 国产精品一区二区三区四区免费观看 | 亚洲欧美日韩东京热| 国产伦一二天堂av在线观看| 国产伦一二天堂av在线观看| 少妇高潮的动态图| 黄色日韩在线| 欧美日韩综合久久久久久 | 热99re8久久精品国产| 99热6这里只有精品| 国产亚洲欧美98| 又爽又黄无遮挡网站| 在线看三级毛片| 99久国产av精品| 中文字幕免费在线视频6| 特大巨黑吊av在线直播| 精品日产1卡2卡| av黄色大香蕉| 欧美潮喷喷水| 国产亚洲精品综合一区在线观看| 成人永久免费在线观看视频| 成人国产综合亚洲| 在线a可以看的网站| 永久网站在线| 亚洲精品一区av在线观看| 亚洲精品成人久久久久久| 精品国产三级普通话版| АⅤ资源中文在线天堂| 午夜两性在线视频| 久久人妻av系列| 特大巨黑吊av在线直播| 国产伦精品一区二区三区视频9| 久久精品国产99精品国产亚洲性色| 欧美日韩福利视频一区二区| 非洲黑人性xxxx精品又粗又长| 亚洲成a人片在线一区二区| 婷婷丁香在线五月| 老女人水多毛片| 欧美一级a爱片免费观看看| 日日摸夜夜添夜夜添av毛片 | 欧美日韩国产亚洲二区| 亚洲无线观看免费| 在线观看免费视频日本深夜| 一个人看的www免费观看视频| 精品一区二区三区人妻视频| 嫩草影视91久久| 无遮挡黄片免费观看| 九九久久精品国产亚洲av麻豆| 少妇的逼好多水| 在线观看美女被高潮喷水网站 | 欧美成人免费av一区二区三区| 性欧美人与动物交配| 在线播放国产精品三级| 国产精品野战在线观看| 国产探花极品一区二区| 女同久久另类99精品国产91| 亚洲精品亚洲一区二区| 亚洲国产精品999在线| 午夜福利欧美成人| 亚州av有码| 中文字幕熟女人妻在线| 天天躁日日操中文字幕| 男人和女人高潮做爰伦理| 国产美女午夜福利| 特大巨黑吊av在线直播| 怎么达到女性高潮| 国产免费男女视频| 在线观看免费视频日本深夜| 国产高潮美女av| 91av网一区二区| 精品乱码久久久久久99久播| 黄色视频,在线免费观看| 日韩欧美精品免费久久 | 国产精品久久久久久久电影| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美精品v在线| 欧美黑人欧美精品刺激| 九九热线精品视视频播放| 日本a在线网址| 18禁裸乳无遮挡免费网站照片| 51国产日韩欧美| 一卡2卡三卡四卡精品乱码亚洲| 亚洲不卡免费看| 国产成人a区在线观看| 中文字幕高清在线视频| 又紧又爽又黄一区二区| 精品久久久久久,| 国产精品久久久久久精品电影| 亚洲自偷自拍三级| 久久久久久国产a免费观看| 舔av片在线| a在线观看视频网站| 非洲黑人性xxxx精品又粗又长| 色尼玛亚洲综合影院| 搡老熟女国产l中国老女人| 欧美日韩综合久久久久久 | 十八禁国产超污无遮挡网站| 可以在线观看毛片的网站| 亚洲,欧美,日韩| av福利片在线观看| 成人毛片a级毛片在线播放| 国产三级在线视频| 精品午夜福利在线看| 欧美性感艳星| 久久性视频一级片| 欧美极品一区二区三区四区| 国产 一区 欧美 日韩| 国产激情偷乱视频一区二区| 亚洲,欧美,日韩| 亚洲中文字幕一区二区三区有码在线看| 亚洲乱码一区二区免费版| 日韩欧美三级三区| 亚洲av一区综合| 国产亚洲精品久久久com| 两人在一起打扑克的视频| 桃色一区二区三区在线观看| 国产综合懂色| 亚洲va日本ⅴa欧美va伊人久久| 国产高清视频在线播放一区| 国产高清激情床上av| 九色成人免费人妻av| 丝袜美腿在线中文| 好男人在线观看高清免费视频| 91麻豆av在线| 国产精品电影一区二区三区| 99国产精品一区二区三区| 91午夜精品亚洲一区二区三区 | 亚洲乱码一区二区免费版| 成熟少妇高潮喷水视频| 黄色丝袜av网址大全| 91九色精品人成在线观看| 久久国产精品影院| 在线观看av片永久免费下载| 男人舔女人下体高潮全视频| 长腿黑丝高跟| 一级毛片久久久久久久久女| 国产精华一区二区三区| 免费黄网站久久成人精品 | 午夜福利18| 宅男免费午夜| 69人妻影院| 一进一出抽搐动态| 国产精品永久免费网站| 亚洲三级黄色毛片| 亚洲国产精品合色在线| 亚洲一区二区三区色噜噜| 日本与韩国留学比较| 国产伦精品一区二区三区视频9| 国产精品女同一区二区软件 | 日本一二三区视频观看| 校园春色视频在线观看| 色综合婷婷激情| 在线观看66精品国产| 91av网一区二区| 国产三级黄色录像| 亚洲精品久久国产高清桃花| 很黄的视频免费| 久久久精品大字幕| 黄色配什么色好看| 午夜福利18| 亚洲内射少妇av| 97热精品久久久久久| 亚洲一区高清亚洲精品| 一本精品99久久精品77| 99国产精品一区二区三区| 色哟哟哟哟哟哟| 精品国产三级普通话版| 久久久久久久午夜电影| 亚洲精品影视一区二区三区av| 成人鲁丝片一二三区免费| 午夜免费男女啪啪视频观看 | 国产欧美日韩一区二区精品| 丰满人妻熟妇乱又伦精品不卡| 内射极品少妇av片p| 91字幕亚洲| 美女黄网站色视频| 51午夜福利影视在线观看| bbb黄色大片| 亚洲天堂国产精品一区在线| 亚洲最大成人中文| 日韩欧美国产一区二区入口| 日本精品一区二区三区蜜桃| 91久久精品国产一区二区成人| 别揉我奶头~嗯~啊~动态视频| 亚洲久久久久久中文字幕| 亚洲无线观看免费| 成人国产一区最新在线观看| 色综合亚洲欧美另类图片| 制服丝袜大香蕉在线| 成人性生交大片免费视频hd| 午夜福利成人在线免费观看| 免费av观看视频| 身体一侧抽搐| 99久久精品热视频| 在线观看午夜福利视频| 日韩欧美国产在线观看| 久久久久免费精品人妻一区二区| 国产免费一级a男人的天堂| 亚洲激情在线av| 高清毛片免费观看视频网站| 国产成人啪精品午夜网站| 国产精品99久久久久久久久| 成人av一区二区三区在线看| 欧美黄色片欧美黄色片| 国模一区二区三区四区视频| 亚洲男人的天堂狠狠| 韩国av一区二区三区四区| 一区二区三区高清视频在线| 国产成年人精品一区二区| 午夜福利在线观看吧| 夜夜夜夜夜久久久久| 亚洲午夜理论影院| ponron亚洲| 日韩欧美一区二区三区在线观看| av视频在线观看入口| 人妻久久中文字幕网| 伊人久久精品亚洲午夜| 国产野战对白在线观看| 亚洲狠狠婷婷综合久久图片| 好男人在线观看高清免费视频| 国内精品久久久久精免费| 国产精品1区2区在线观看.| 中文字幕人成人乱码亚洲影| 性色av乱码一区二区三区2| 一进一出抽搐动态| av在线老鸭窝| 在线天堂最新版资源| 午夜a级毛片| 日韩av在线大香蕉| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 搡老岳熟女国产| 露出奶头的视频| 亚洲成a人片在线一区二区| 18+在线观看网站| 最好的美女福利视频网| 亚洲精品粉嫩美女一区| 午夜老司机福利剧场| 尤物成人国产欧美一区二区三区| 日本成人三级电影网站| 亚洲18禁久久av| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 一进一出抽搐gif免费好疼| 亚洲国产欧美人成| 国产精品久久电影中文字幕| 真实男女啪啪啪动态图| 亚洲国产高清在线一区二区三| 国产精品一区二区三区四区久久| 1024手机看黄色片| 看免费av毛片| 亚洲电影在线观看av| 露出奶头的视频| 国产v大片淫在线免费观看| 国产熟女xx| 国产精品影院久久| 亚洲无线在线观看| 亚洲成人久久爱视频| 中文字幕高清在线视频| 午夜免费成人在线视频| 中文字幕熟女人妻在线| 最新在线观看一区二区三区| 婷婷亚洲欧美| 丰满人妻一区二区三区视频av| 好男人在线观看高清免费视频| 亚洲三级黄色毛片| av在线观看视频网站免费| 夜夜爽天天搞| 男女那种视频在线观看| 我的老师免费观看完整版| 国产成人福利小说| 99热这里只有是精品在线观看 | 听说在线观看完整版免费高清| 日本黄色片子视频| 免费看光身美女| 日韩欧美一区二区三区在线观看| 一个人观看的视频www高清免费观看| 乱码一卡2卡4卡精品| 在线天堂最新版资源| 天堂√8在线中文| 午夜激情欧美在线| 亚洲va日本ⅴa欧美va伊人久久| 国产免费av片在线观看野外av| 欧美xxxx黑人xx丫x性爽| 精品一区二区三区av网在线观看| 久久久久久九九精品二区国产| 最近中文字幕高清免费大全6 | 精品久久久久久久末码| 国产在视频线在精品| 91久久精品电影网| 最近中文字幕高清免费大全6 | 久久天躁狠狠躁夜夜2o2o| 啦啦啦观看免费观看视频高清| 亚洲av成人不卡在线观看播放网| 深爱激情五月婷婷| 黄色丝袜av网址大全| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清| 色综合欧美亚洲国产小说| 在线天堂最新版资源| 长腿黑丝高跟| 男人狂女人下面高潮的视频| 国产老妇女一区| a级毛片免费高清观看在线播放| 极品教师在线视频| 中文亚洲av片在线观看爽| 淫秽高清视频在线观看| 人妻久久中文字幕网| 日本黄色视频三级网站网址| 亚洲最大成人av| 国内久久婷婷六月综合欲色啪| ponron亚洲| 国产精品女同一区二区软件 | 人妻久久中文字幕网| 亚洲精品一卡2卡三卡4卡5卡| 精品午夜福利在线看| x7x7x7水蜜桃| 黄色配什么色好看| 亚洲国产欧洲综合997久久,| 麻豆av噜噜一区二区三区| 老司机深夜福利视频在线观看| av欧美777| 成年人黄色毛片网站| 国产高潮美女av| 成年女人看的毛片在线观看| 国产精品一及| 一级黄片播放器| .国产精品久久| 欧美性猛交黑人性爽| 免费看日本二区| 一区二区三区激情视频| a级毛片免费高清观看在线播放| 久久草成人影院| 国产视频一区二区在线看| 噜噜噜噜噜久久久久久91| 男人狂女人下面高潮的视频| 国产亚洲精品久久久久久毛片| 国产国拍精品亚洲av在线观看| 两人在一起打扑克的视频| 在线天堂最新版资源| 综合色av麻豆| 尤物成人国产欧美一区二区三区| 亚洲 国产 在线| 性色av乱码一区二区三区2| 国产探花在线观看一区二区| 婷婷色综合大香蕉| 欧美色视频一区免费| 伊人久久精品亚洲午夜| 精品乱码久久久久久99久播| 最新在线观看一区二区三区| 永久网站在线| 久久久久久久精品吃奶| 欧美一级a爱片免费观看看| 搡老岳熟女国产| 久久午夜亚洲精品久久| 搞女人的毛片| 亚洲av免费在线观看| 国产精品亚洲av一区麻豆| 亚洲国产高清在线一区二区三| 有码 亚洲区| 黄色丝袜av网址大全| 国产高潮美女av| 免费搜索国产男女视频| 成人精品一区二区免费| 国产不卡一卡二| x7x7x7水蜜桃| 一级黄色大片毛片| 久久久久久国产a免费观看| 婷婷丁香在线五月| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩精品一区二区| 草草在线视频免费看| 亚洲成a人片在线一区二区| 久久精品91蜜桃| 国产黄色小视频在线观看| 色综合婷婷激情| 哪里可以看免费的av片| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 久久6这里有精品| 天堂av国产一区二区熟女人妻| 九九久久精品国产亚洲av麻豆| 老司机午夜十八禁免费视频| 久久精品91蜜桃| 欧美色视频一区免费| 欧美激情在线99| 怎么达到女性高潮| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 小说图片视频综合网站| av在线天堂中文字幕| 国产色婷婷99| 色综合亚洲欧美另类图片| 国产日本99.免费观看| 国产午夜福利久久久久久| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久免费视频| 18禁黄网站禁片免费观看直播| 日韩 亚洲 欧美在线| 日本黄色视频三级网站网址| 日日干狠狠操夜夜爽| 真人做人爱边吃奶动态| avwww免费| 很黄的视频免费| 悠悠久久av| 亚洲精品影视一区二区三区av| 国产午夜精品论理片| 亚洲国产欧洲综合997久久,| 国产单亲对白刺激| 韩国av一区二区三区四区| 变态另类成人亚洲欧美熟女| 欧美日本视频| 91在线观看av| 亚洲中文日韩欧美视频| 国产精品久久电影中文字幕| www.www免费av| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 99久久99久久久精品蜜桃| 久久久久久久久久黄片| 国产精品永久免费网站| 中文字幕人成人乱码亚洲影| 两人在一起打扑克的视频| 亚洲 欧美 日韩 在线 免费| 99国产精品一区二区蜜桃av| 91九色精品人成在线观看| 精品久久国产蜜桃| 精品一区二区三区av网在线观看| 国产精品一及| 天堂影院成人在线观看| 国产高清激情床上av| 老司机福利观看| 久久人人精品亚洲av| 午夜激情欧美在线| 12—13女人毛片做爰片一| 禁无遮挡网站| 九九在线视频观看精品| 一区二区三区四区激情视频 | av专区在线播放| 午夜a级毛片| 成年女人毛片免费观看观看9| 欧美最黄视频在线播放免费| 人妻夜夜爽99麻豆av| 十八禁人妻一区二区| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 黄色一级大片看看| 制服丝袜大香蕉在线| 9191精品国产免费久久| 露出奶头的视频| 国产av麻豆久久久久久久| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 97人妻精品一区二区三区麻豆| 老司机午夜十八禁免费视频| 亚洲 国产 在线| 99在线视频只有这里精品首页| 熟女电影av网| 18禁裸乳无遮挡免费网站照片| 看十八女毛片水多多多| 欧美性感艳星| 成人一区二区视频在线观看| 亚洲av成人av| 久久久久久久亚洲中文字幕 | 蜜桃亚洲精品一区二区三区| 91av网一区二区| 久久草成人影院| 搡老岳熟女国产| 男人狂女人下面高潮的视频| 婷婷精品国产亚洲av在线| 国产人妻一区二区三区在| 日韩欧美 国产精品| 国产淫片久久久久久久久 | 免费av毛片视频| 亚洲精品乱码久久久v下载方式| 色噜噜av男人的天堂激情| 91麻豆av在线| 3wmmmm亚洲av在线观看| 99热这里只有精品一区| 精品国产三级普通话版| 免费黄网站久久成人精品 | 首页视频小说图片口味搜索| 无遮挡黄片免费观看| 如何舔出高潮| 精华霜和精华液先用哪个| 成年女人毛片免费观看观看9| 国产野战对白在线观看| 女人被狂操c到高潮| av欧美777| 久久久久精品国产欧美久久久| 国产免费av片在线观看野外av| 一本精品99久久精品77| 在线免费观看的www视频| 国产在线精品亚洲第一网站| 亚洲av成人av| 国产av麻豆久久久久久久| 亚洲av日韩精品久久久久久密| 99久久99久久久精品蜜桃| 欧美性感艳星| 天天一区二区日本电影三级| 日韩av在线大香蕉| 三级男女做爰猛烈吃奶摸视频| 成人无遮挡网站| 日韩中文字幕欧美一区二区| 亚洲av免费高清在线观看| 成人毛片a级毛片在线播放| 国产中年淑女户外野战色| 人人妻,人人澡人人爽秒播| 欧美丝袜亚洲另类 | 12—13女人毛片做爰片一| 精品人妻熟女av久视频| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | 成人三级黄色视频| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 老女人水多毛片| 十八禁网站免费在线| 日日夜夜操网爽| or卡值多少钱| 日本熟妇午夜| 久久天躁狠狠躁夜夜2o2o| 桃色一区二区三区在线观看| 黄色女人牲交| 国产精品电影一区二区三区| 一个人免费在线观看电影| 国产精品98久久久久久宅男小说| 欧美成狂野欧美在线观看| 人人妻,人人澡人人爽秒播|