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

    基于改進(jìn)帕默爾干旱指數(shù)的中國氣象干旱時(shí)空演變分析

    2016-03-21 12:37:56王兆禮黃澤勤鐘睿達(dá)陳佳穎邱卓輝華南理工大學(xué)土木與交通學(xué)院廣州510641
    關(guān)鍵詞:時(shí)空演變干旱中國

    王兆禮,李 軍,黃澤勤,鐘睿達(dá),陳佳穎,邱卓輝(華南理工大學(xué)土木與交通學(xué)院,廣州 510641)

    ?

    基于改進(jìn)帕默爾干旱指數(shù)的中國氣象干旱時(shí)空演變分析

    王兆禮,李軍,黃澤勤,鐘睿達(dá),陳佳穎,邱卓輝
    (華南理工大學(xué)土木與交通學(xué)院,廣州 510641)

    摘要:近幾十年來隨著全球變暖的不斷加劇,中國多地頻繁發(fā)生干旱災(zāi)害,經(jīng)濟(jì)損失也越來越大,因而研究干旱的時(shí)空演變對社會經(jīng)濟(jì)發(fā)展具有重要的現(xiàn)實(shí)意義。論文基于改進(jìn)帕默爾干旱指數(shù)(scPDSI)分析了1961—2009年間年和四季干濕變化,并利用旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(rotated empirical orthogonal function,REOF)探討了中國氣象干旱的時(shí)空演變特征。結(jié)果表明:1)近49 a來,整體上中國年和四季均呈顯著變濕趨勢,且均在20世紀(jì)70年代初期發(fā)生了由干變濕的突變;2)年及四季均存在2~8 a的振蕩周期,其中夏季主周期為4.4 a,其余主周期均為6.2 a;3)根據(jù)REOF時(shí)空分解的前8個(gè)空間模態(tài),將中國劃分成8個(gè)干濕特征區(qū)域,其中大興安嶺地區(qū)、東南地區(qū)、西北地區(qū)和青藏高原南部地區(qū)有變濕趨勢,后兩個(gè)區(qū)域變濕趨勢顯著;東北—內(nèi)蒙古高原區(qū)、青藏高原北部地區(qū)、中部地區(qū)和黃淮海平原地區(qū)呈變干趨勢,前3個(gè)區(qū)變干趨勢顯著,黃淮海平原分區(qū)干濕變化不明顯;各分區(qū)普遍具有2~9 a的振蕩周期;4)極渦指數(shù)、印度洋偶極子和太平洋濤動與部分地區(qū)氣象干旱有較好的相關(guān)關(guān)系。研究結(jié)果可為研究區(qū)防旱、抗旱與干旱預(yù)測提供參考依據(jù)。

    關(guān)鍵詞:干旱;氣象;降水;改進(jìn)帕默爾干旱指數(shù);旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù);時(shí)空演變;中國

    王兆禮,李軍,黃澤勤,鐘睿達(dá),陳佳穎,邱卓輝. 基于改進(jìn)帕默爾干旱指數(shù)的中國氣象干旱時(shí)空演變分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(2):161-168.doi:10.11975/j.issn.1002-6819.2016.02.024http://www.tcsae.org

    Wang Zhaoli, Li Jun, Huang Zeqin, Zhong Ruida, Chen Jiayin, Qiu Zhuohui. Spatiotemporal variations analysis of meteorological drought in China based on scPDSI[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(2): 161-168. (in Chinese with English abstract)doi:10.11975/j.issn.1002-6819.2016.02.024http://www.tcsae.org

    Email:wangzhl@scut.edu.cn

    0 引 言

    干旱是一種自然災(zāi)害,對農(nóng)業(yè)生產(chǎn)、水資源供給和環(huán)境生態(tài)系統(tǒng)產(chǎn)生不良影響[1-2]。氣候變化伴隨著人類活動增加(如,森林砍伐、大規(guī)模修建水庫等)造成全球氣溫升高,這加快全球的水文循環(huán),影響干旱發(fā)生頻次,甚至可能導(dǎo)致旱災(zāi)的發(fā)生[3]。干旱發(fā)生在世界上大部分地區(qū),甚至在濕潤地區(qū)[4]。自20世紀(jì)50年代以來,大范圍的干旱在世界各地都有發(fā)生,例如,澳大利亞、非洲、亞洲、歐洲和美洲等地區(qū)[5]。近10年來中國多地發(fā)生嚴(yán)重旱災(zāi)。例如,長江流域多次發(fā)生旱災(zāi),在2006年長江徑流量跌落到近50年來最底點(diǎn)[6]及2009-2010年在中國西南地區(qū)發(fā)生冬春極端嚴(yán)重干旱[7]。旱災(zāi)給社會經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境帶來諸多不利影響,故分析中國的干旱時(shí)空演變具有重要的現(xiàn)實(shí)意義。

    由于干旱的復(fù)雜性,使用不同的干旱指數(shù)、數(shù)據(jù)集或者計(jì)算蒸發(fā)量方法來評估干旱往往得出不同的結(jié)果[8]。例如,Dai認(rèn)為全球干旱呈增加趨勢[9],但Sheffield等研究發(fā)現(xiàn)過去60年里全球干旱狀況幾乎沒有發(fā)生變化[10]。眾學(xué)者也對中國干旱進(jìn)行諸多研究,獲得了若干新的認(rèn)識,但也出現(xiàn)上述類似結(jié)果,難有共識性結(jié)論。一種觀點(diǎn)認(rèn)為中國干旱無明顯變化,例如,Wang等運(yùn)用SPEI 和SPI指數(shù)分析中國干濕變化,并認(rèn)為1961-2012年期間中國干旱沒有呈增加或減少趨勢[11]。Zou等發(fā)現(xiàn)1951 -2003年間,整個(gè)中國的干旱狀況并無增加或下降的趨勢[12]。徐仁貴等基于蒸散量和濕潤指數(shù),研究發(fā)現(xiàn)中國極端干旱區(qū)面積在顯著減小,西北蒸散量呈顯著減小趨勢[13]。翟盤茂等認(rèn)為1951-2003年間中國干旱面積略有增加,但變化趨勢并不顯著[14]。第二種觀點(diǎn)認(rèn)為中國干旱呈增加趨勢,例如,Li等的研究表明,在1951-2005年間中國土壤水分減少趨勢顯著,干旱加劇[15]。Wang等結(jié)合土壤中水分變化,分析了中國1950-2006年中國干旱的嚴(yán)重程度、持續(xù)時(shí)間和頻率,結(jié)果表明由土壤中缺乏水分引起的干旱變得更加頻繁,在東北和中部地區(qū)這種情況更加顯著,這容易引起農(nóng)業(yè)干旱[16]。Yu等[17]表明近50年來,受旱面積每10年以3.72%的速度增長,且自1990年代以來,中國嚴(yán)重和極端干旱災(zāi)害呈增加趨勢。葉敏等[18]分析中國旱澇時(shí)空變化特征,認(rèn)為21世紀(jì)以來,中國干旱發(fā)生變頻繁,且強(qiáng)度也增加。陳方藻等[19]研究發(fā)現(xiàn)近60年來中國干旱的發(fā)生整體呈上升趨勢。

    上述研究成果為中國干旱的研究提供重要的參考依據(jù),但也出現(xiàn)了類似全球研究的兩種矛盾性結(jié)論。那么,中國氣象干旱的變化趨勢到底如何,是否具有突變特性和顯著的周期性還有待進(jìn)一步深入研究。再者,已有的大多研究使用的干旱指數(shù)只牽涉單一的降雨量(如SPI指數(shù)),或所用的干旱指數(shù)采用桑斯維特(Thornthwaite)算法計(jì)算蒸發(fā)量,忽略了輻射和空氣動力控制等因素(如SPEI_Th,PDSI_Th),不能精確評估干旱變化[8]。因而采用精確性更高的scPDSI_PM指數(shù)(使用Penman-Monteith公式計(jì)算蒸散發(fā))分析中國氣象干旱狀況就顯得十分必要了。

    鑒于此,本文基于2013年發(fā)布的高空間分辨率scPDSI_PM數(shù)據(jù)集,采用線性回歸、Mann-Kendall方法、小波分析和旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(rotated empirical orthogonal function,REOF)等方法,研究1961-2009年中國氣象干旱的時(shí)空分布和演變規(guī)律,并探究干旱的大氣環(huán)流影響因素,以期為研究區(qū)防旱、抗旱提供理論依據(jù)和科學(xué)參考。

    1 研究區(qū)概況和研究方法

    1.1研究區(qū)域

    中國幅員遼闊,跨經(jīng)緯度較廣(73°40′~135°2′E,3°52′~53°33′N),位于亞歐大陸東部,太平洋西岸,距印度洋較近,中國氣候受大陸、大洋影響顯著。由于地形復(fù)雜,氣候類型復(fù)雜多樣,地跨多個(gè)溫度帶,可劃分為多個(gè)不同干濕特征區(qū)域,在同一氣候類型中,熱量和干濕程度差異較大。中國是全球氣象災(zāi)害最頻繁、最嚴(yán)重的國家之一,其中干旱又是最嚴(yán)重的災(zāi)害之一。過去干旱普遍發(fā)生在中國北區(qū)地區(qū),且造成嚴(yán)重危害,但近年來南方地區(qū)嚴(yán)重干旱事件也頻繁發(fā)生[20]。在1959、1960、1972、1978、1997、2006與2009-2010年等年份中國都發(fā)生嚴(yán)重旱災(zāi),小干旱災(zāi)害幾乎每年都會發(fā)生,平均每年有千萬公頃的農(nóng)田因旱受災(zāi),占各類災(zāi)害影響面積的60%以上[21]。

    1.2數(shù)據(jù)獲取與處理

    本文的scPDSI_PM數(shù)據(jù)[22](self-calibrating palmer drought severity index based on Penman-Monteith)來源于英國東英格利亞大學(xué)氣候研究中心(climatic research unit,CRU)編輯的氣候數(shù)據(jù)集。該數(shù)據(jù)集基于CRU TS 3.10.01版本氣候數(shù)據(jù),具有較高的精度[22],可從英國氣象數(shù)據(jù)中心下載(http://www.badc.rl.ac.uk/)。本研究所選的中國地區(qū)3 725個(gè)網(wǎng)格點(diǎn)的數(shù)據(jù)資料均來自scPDSI_PM數(shù)據(jù)集,時(shí)間范圍為1961-2009年,空間分辨率為0.5°×0.5°。有關(guān)土壤田間持水量、徑流的推求過程可參考文獻(xiàn)[22]。該數(shù)據(jù)集主要有3個(gè)優(yōu)勢:1)在多樣的氣候條件下,scPDSI有相似的變動范圍,使其更適合分析不同地區(qū)的有效水分。2)使用物理意義更清晰的Penman-Monteith公式(PM)計(jì)算蒸散量,用實(shí)際的植被代替參考作物, 植被數(shù)據(jù)來源于美國地質(zhì)調(diào)查局, 為衛(wèi)星遙感解譯數(shù)據(jù)(1 km分辨率)。3)計(jì)算水分平衡時(shí)考慮了季節(jié)性積雪動態(tài)變化[22]。中國基礎(chǔ)地理信息(ArcGIS格式)由中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供(http://www.resdc.cn)。74項(xiàng)大氣環(huán)流數(shù)據(jù)由國家氣候中心氣候系統(tǒng)診斷預(yù)測室提供(http://ncc.cma.gov.cn);太平洋濤動(Pacific Oscillation, PDO)數(shù)據(jù)來源于Earth System Research Laboratory (ESRL) of NOAA(http://www.esrl.noaa.gov/psd/ data/correlation/pdo.data);印度洋偶極子(Indian Ocean Dipole,IOD)數(shù)據(jù)由日本低緯度氣候預(yù)測研究所提供(http://www.jamstec.go.jp/frsgc/research/d1/iod/e/index.html)。數(shù)據(jù)處理過程為:首先運(yùn)用ArcGIS10.2軟件裁減出中國部分,再運(yùn)用R軟件計(jì)算年均和四季scPDSI值及其MK值;采用Matlab軟件做M-K突變檢測圖及小波分析圖,最后再運(yùn)用ArcGIS10.2軟件獲得REOF空間分區(qū)圖。

    1.3 研究方法

    1.3.1改進(jìn)的帕默爾干旱指數(shù)(scPDSI)

    PDSI廣泛運(yùn)用于干旱評估,由Palmer提出的[23],最初是為測量陸地表面水分供需的累積距平。PDSI考慮氣溫因素,可以有效反映氣候變化對于干旱的影響[23]。此外,其表征在一段時(shí)間內(nèi),該地區(qū)實(shí)際水分供應(yīng)持續(xù)少于當(dāng)?shù)貧夂蜻m宜水分供應(yīng)的水分虧缺[23]。但因PDSI對于不同空間的干旱分析有差異,認(rèn)為不適用于評估不同地區(qū)的干旱[22]。而scPDSI較PDSI有較大改進(jìn),且使用FAO PM公式計(jì)算蒸散發(fā),精確度更高[22]。此外,scPDSI用各自站點(diǎn)的氣候資料得出的值取代氣候權(quán)重因子和持續(xù)因子,故較少受地域約束,具有更大的空間可比性[24],具體PDSI和scPDSI計(jì)算方法詳見參考文獻(xiàn)[25-26]。scPDSI指數(shù)的干旱等級分布,詳見表1。

    表1 scPDSI干濕等級劃分標(biāo)準(zhǔn)Table 1 Classification of scPDSI values

    1.3.2統(tǒng)計(jì)方法

    對于中國氣象干旱的時(shí)間變化主要采取趨勢分析、突變分析和周期分析,空間變化則采用旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(REOF)。本文采用Mannn-Kendall法(M-K法)[27-28]對年和四季scPDSI值進(jìn)行突變檢測。M-K法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法,具有不需要樣本遵從一定的分布,不受少數(shù)異常值干擾,檢測范圍寬廣和定量化程度高的優(yōu)點(diǎn)。小波分析(continuous wavelet transform,CWT)[29]是建立在經(jīng)典的傅里葉分析和窗口傅里葉變換的基礎(chǔ)上發(fā)展起,是國際上十分熱門的一種周期分析方法。它可以很好表征時(shí)間序列的周期振蕩現(xiàn)象,可以提供連續(xù)、平滑的小波能量譜圖。

    旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(REOF)[30]是基于經(jīng)驗(yàn)正交函數(shù)(empirical orthogonal function,EOF),進(jìn)一步分析要素場空間分布的局部區(qū)域特征,可更好表現(xiàn)地域差異,同時(shí)反映不同地域的相關(guān)分布,將氣候要素一致的地區(qū)劃分為同一的區(qū)域。上述方法的計(jì)算步驟詳見參考文獻(xiàn)[31]。

    2 結(jié)果與分析

    2.1中國干濕變化分析及周期分析

    2.1.1年與四季干濕變化趨勢和突變檢測

    為研究1961-2009年中國年、四季的干濕突變情況,采用M-K法對年和四季的scPDSI時(shí)間序列進(jìn)行突變分析(圖1)。由圖1可知,UF曲線自20世紀(jì)70年代呈上升趨勢,表明整體上中國呈變濕趨勢。一直到21世紀(jì)初期,這種變濕趨勢超過0.05臨界線,表明變濕趨勢顯著(P<0.05);但2005年以來,又呈現(xiàn)明顯的變干趨勢,這與近幾年中國多地頻繁發(fā)生干旱災(zāi)害相吻合。UF與UB曲線相交位置可知,1971年為由干變濕的突變點(diǎn)。基于面積平均的中國春季、夏季、秋季、冬季與年系列的變化趨勢基本一致(圖略),均呈顯著的變濕趨勢,1970年代初期為變濕的突變點(diǎn),均通過了檢驗(yàn)(P<0.05)。

    圖1 中國年scPDSI指數(shù)M-K突變檢測圖Fig.1 Mann-Kendall test of annual scPDSI across China

    2.1.2年及四季的干濕變化周期分析

    為研究干濕序列不同時(shí)間尺度的變化特征,使用小波分析其周期變化特征。由2 a小波功率圖可知,在1961-2009年scPDSI序列存在2.6、4.4、6.2、14.7和24.8 a的周期,其中2.6、4.4和6.2 a的周期通過95%的紅噪聲檢驗(yàn),6.2a為該序列的主周期。在連續(xù)小波變換圖中,能量密度的峰值主要集中在1960年代、1980年代到90年代中期、80年代中期到2000年和1950年代到21世紀(jì)初期。其中,1980年代到90年代中期、80年代中期到2000年分別存在5~8 a和2~3 a的振蕩周期,并均通過紅噪聲檢驗(yàn),此外1960年代和1980年代到21世紀(jì)初期存在2~4 a和4~8 a的振蕩周期,但沒有通過紅噪聲檢驗(yàn),可知基于面積平均的中國年干濕變化存在2.0~8.0 a的振蕩周期。

    圖2 中國年干濕變化小波分析圖Fig.2 Wavelet analysis map of average annual dry/wet changes in China

    春、秋和冬季干濕變化周期基本一致(圖略),均存在2.6、6.2、14.7和24.8a的周期,其中2.6和6.2 a的周期經(jīng)過紅噪聲檢測,主周期均為6.2 a。夏季存在3.1、4.4、14.7和24.8 a的周期,其中3.1和4.4 a的周期經(jīng)過紅噪聲檢測,主周期為4.4 a。四季干濕變化基本一致,存在2~8 a 的振蕩周期。

    2.2中國氣象干旱時(shí)空演變特征分析

    2.2.1REOF空間分區(qū)

    為研究中國干濕變化空間分布規(guī)律,對中國49 a的年和四季scPDSI值作旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(REOF)分解,結(jié)果表明:年和四季空間分區(qū)基本一致。限于篇幅,本文僅分析年尺度空間模態(tài),前8個(gè)主分量旋轉(zhuǎn)后的方差貢獻(xiàn)分別為14.3%、9.8%、8.6%、7.2%、6.9%、5.5%、5.1%、4.5%,累積解釋方差貢獻(xiàn)占總方差的61.9%。

    圖3是年scPDSI值前8個(gè) REOF旋轉(zhuǎn)載荷向量(rotated load vector,RLV)的空間分布。圖3a為第1空間模態(tài)(REOF1),RLV1正的高值區(qū)位于中國西北部地區(qū),其中心旋轉(zhuǎn)載荷向量值為0.055。圖3b為第2空間模態(tài)(REOF2),RLV2正的高值區(qū)位于東北—內(nèi)蒙古高原地區(qū),中心旋轉(zhuǎn)載荷向量值為0.053。第3旋轉(zhuǎn)載荷向量RLV3正的高值區(qū)位于黃淮海平原地區(qū),中心旋轉(zhuǎn)載荷向量值為0.061(圖3c)。圖3d為第4空間模態(tài)(REOF4),RLV4負(fù)的高值位于青藏高原北部地區(qū),中心旋轉(zhuǎn)載荷向量值為?0.057。第5空間模態(tài)(REOF5),RLV5正的高值區(qū)位于中國中部地區(qū)(包括四川盆地、漢中盆地和黃土高原等地區(qū)),其中心旋轉(zhuǎn)載荷向量值為0.055(圖3e)。圖2f為第6空間模態(tài)(REOF6),RLV6的高值區(qū)位于青藏高原南部地區(qū),其中心旋轉(zhuǎn)載荷向量值為0.067。第7空間模態(tài)(圖3g),RLV7的高值區(qū)位于大興安嶺地區(qū),其中心旋轉(zhuǎn)載荷向量值為0.066。第8空間模態(tài)(圖3h),RLV8負(fù)的高值區(qū)位于中國東南地區(qū),其中心旋轉(zhuǎn)載荷向量值為?0. 060。

    根據(jù)各模態(tài)在空間上的高載荷區(qū),將中國劃分成8個(gè)干濕特征區(qū)域:西北地區(qū)、東北—內(nèi)蒙古高原地區(qū)、黃淮海平原地區(qū)、青藏高原北部地區(qū)、中部地區(qū)、青藏高原南部地區(qū)、大興安嶺地區(qū)、東南地區(qū)。相鄰區(qū)域基本沒有重疊,REOF分區(qū)結(jié)果合理地反映了中國地形的實(shí)際情況。

    2.2.2 各分區(qū)年和四季干濕狀況時(shí)間演變趨勢分析

    REOF分解得到的旋轉(zhuǎn)主成分(rotated principal components,RPC)序列實(shí)際上反映了各分區(qū)干濕狀況的時(shí)間變化特征,本文分析8個(gè)分區(qū)年和四季干濕狀況的時(shí)間變化M-K值,結(jié)果見表2。限于篇幅,本文具體分析年尺度干濕變化特征。

    由表2可知,西北地區(qū)、青藏高原南部地區(qū)、大興安嶺地區(qū)和東南地區(qū)有變濕趨勢,其前兩個(gè)區(qū)變濕趨勢顯著(P<0.05);而東北—內(nèi)蒙古高原地區(qū)、青藏高原北部地區(qū)、中部地區(qū)和黃淮海平原地區(qū)呈現(xiàn)出變干趨勢,其中,前3個(gè)分區(qū)變干趨勢顯著(P<0.05),黃淮海平原地區(qū)變干趨勢不明顯。此外,西北地區(qū)年和四季M-K值均大于1.96,表明全年變濕趨勢顯著(P<0.05);而中部地區(qū)年和四季M-K值均小于-1.96,全年變干趨勢顯著(P<0.05)。其他分區(qū)干濕變化季節(jié)特征較明顯,存在一定差異性,但大部分地區(qū)春季均呈變干趨勢,夏季均呈現(xiàn)變濕趨勢。

    圖3 前8個(gè)REOF的旋轉(zhuǎn)載荷向量空間分布Fig.3 Spatial distributions of rotated load vector of first eight REOFs

    表2 中國各分區(qū)年和四季干濕變化M-K值Table 2 M-K trend test of average annual and season dry/wet changes in China

    2.2.3 各分區(qū)年干濕變化周期

    為了進(jìn)一步探索中國49 a來不同區(qū)域干濕的年際及年代際變化特征。采用小波分析方法研究前8個(gè)空間模態(tài)的干濕狀況(即旋轉(zhuǎn)主成分rotated principal components,RPC)的周期特征(圖4)。由圖4可知,不同分區(qū)干濕變化的時(shí)間尺度分布不均勻,各分區(qū)具有2~9 a的振蕩周期,且具有明顯的局部化特征,干濕變化具有區(qū)域差異性。西北地區(qū)(圖4a)1993-2005年存在3~5a振蕩周期,1997-2009年存在6~9 a振蕩周期。由圖4b可知,東北-內(nèi)蒙古高原地區(qū)1982-2009年存在3~9 a的振蕩周期。在黃淮海平原地區(qū)1961-1969年存在4~8 a的周期,1987-1994存在4~6年的周期(圖4c)。青藏高原北部地區(qū)(圖4d)1961-1967年存在2~5 a的周期,1981-1988存在6~7 a的周期。在中部地區(qū)(圖4e),3~6 a的周期大約出現(xiàn)在1961-1970年,1970年至2009年存在6~8 a的周期。由圖4f可知,青藏高原南部地區(qū)1973-1984年存在3~4 a的周期,1980-2003年存在8~10 a的周期。大興安嶺地區(qū)(圖4g)僅在1977-1994年存在3~8 a的周期。在東南地區(qū)(圖4h),1968-1977出現(xiàn)2~4 a的振蕩周期,1994-2005年出現(xiàn)4~6 a的振蕩周期。

    圖4 REOF分解的8個(gè)旋轉(zhuǎn)主成分小波分析Fig.4 Continuous wavelet transform of annual rotated PC series

    2.3與大尺度環(huán)流的關(guān)系

    相關(guān)研究表明,中國地處亞洲季風(fēng)區(qū),容易受各種大氣環(huán)流的影響[32-35]。但有關(guān)大尺度大氣環(huán)流指標(biāo)對中國干旱的演變是否有影響的研究還較薄弱,故本文用年scPDSI值及各分區(qū)旋轉(zhuǎn)主成分(RPC)與大氣環(huán)流因子做相關(guān)分析(表3),識別可能影響中國干旱的大氣環(huán)流因子。

    表3 全國和各分區(qū)scPDSI與大氣環(huán)流指數(shù)的關(guān)系Table 3 Correlation between scPDSI and circulation indices of subregion and whole country

    由表3可知,極渦指數(shù)跟年均scPDSI值密切相關(guān),當(dāng)極渦指數(shù)低值時(shí),scPDSI值較大。西北地區(qū)不僅受極渦指數(shù)影響,也受印度洋偶極子(IOD)和太平洋濤動(PDO)影響,與其呈正相關(guān)趨勢。IOD與東北—內(nèi)蒙古高原地區(qū)呈負(fù)相關(guān)關(guān)系。黃淮海平原地區(qū)PDO呈負(fù)相關(guān)關(guān)系。青藏高原北部地區(qū)與極渦指數(shù)呈正相關(guān)關(guān)系,但與IOD呈負(fù)相關(guān)關(guān)系。中部地區(qū)受北美、亞洲極渦面積指數(shù)影響,與其呈正相關(guān)關(guān)系,但與IOD呈負(fù)相關(guān)關(guān)系。極渦指數(shù)與青藏高原南部地區(qū)呈負(fù)相關(guān)關(guān)系,IOD與青藏高原南部地區(qū)呈正相關(guān)關(guān)系。北半球極渦面積指數(shù)與大興安嶺地區(qū)呈負(fù)相關(guān)關(guān)系,亞洲極渦強(qiáng)度指數(shù)與東南地區(qū)呈正相關(guān)關(guān)系。

    以上相關(guān)分析表明,極渦指數(shù)與中國大部分地區(qū)氣象干旱有較好相關(guān)性。前人眾多研究結(jié)果認(rèn)為極渦與中國不少區(qū)域降水、氣溫有一定的聯(lián)系。例如,極渦指數(shù)與后期中國氣溫呈負(fù)相關(guān),在春季,若太平洋區(qū)極渦面積異常大,亞洲區(qū)極渦面積、北美區(qū)極渦面積強(qiáng)度異常弱,長江、黃河之間中游地區(qū)的降水出現(xiàn)負(fù)距平,廣東、福建沿海降水出現(xiàn)正距平[36-37]。印度洋偶極子與中國部分地區(qū)降雨密切相關(guān),如IOD指數(shù)為正位相時(shí),華北及江淮流域夏季干旱少雨,而華南沿海地區(qū)為多雨帶。此外云南、河套地區(qū)北部、東北北部也為正的降水距平區(qū)。當(dāng)印度洋偶極指數(shù)為負(fù)位相,夏季南北地區(qū)多雨而長江流域少雨[38]。北半球冬季IOD正異常峰值后一年春季歐亞中高緯度地區(qū)呈顯著的負(fù)歐亞太平洋型遙相關(guān),長江以北(特別是黃河中上游地區(qū))大部分地區(qū)持續(xù)多雨,反之,中國東部大部分地區(qū)降雨偏少,容易發(fā)生大范圍干旱。太平洋濤動是影響中國降雨的另一重要因素[39]。亞洲—太平洋濤動指數(shù)可很好指示中國東部水變化。當(dāng)PDO指數(shù)偏高時(shí),長江以北和以南、南亞地區(qū)的降水增加,而長江流域的降水減少[40]。

    3 討 論

    論文應(yīng)用scPDSI_PM指數(shù)研究表明,中國年和四季變濕趨勢均顯著且存在2~8 a的振蕩周期,該結(jié)論有別于前人[11-19]的研究結(jié)論。究其原因主要是采用PM公式計(jì)算蒸散發(fā)所導(dǎo)致。已有的研究多采用桑斯維特(Thornthwaite,Th)算法計(jì)算蒸發(fā)量,其在中國的效用是值得懷疑的[41],用Th算法評估干旱的結(jié)果很可能會被夸大[10]。而Penman-Monteith(PM)算法考慮了輻射因素、風(fēng)速等因素,這可以更加精確地評估干旱[10]。Zheng 等[42]研究表明風(fēng)速、太陽輻射和濕度對中國地區(qū)的蒸發(fā)量有較大影響。而由于輻射、水汽壓差以及(或者)風(fēng)速的減小,可能最終導(dǎo)致蒸發(fā)的減弱[10]。Chen等[41]用PM算法計(jì)算中國1956-2005年的蒸發(fā)量,發(fā)現(xiàn)蒸發(fā)量顯著減少。然而,由于氣溫的顯著增加,采用Th公式計(jì)算得到的潛在蒸發(fā)量大都呈增加趨勢[10,17,41]。因而,在中國大部分地區(qū)蒸發(fā)量減少的背景下,scPDSI指數(shù)采用PM算法計(jì)算蒸發(fā)量,這可能使得scPDSI_PM值更加準(zhǔn)確,在此基礎(chǔ)上得出的相關(guān)研究結(jié)論可能更接于實(shí)際情況。后續(xù)的研究可借鑒Nɑture期刊討論全球干旱變化趨勢研究的思路與方法[8-10],深入探討Th與PM公式以及不同數(shù)據(jù)源降水?dāng)?shù)據(jù)對中國氣象干旱變化趨勢分析的影響。

    4 結(jié) 論

    根據(jù)中國1961-2009年間高空間分辨率(0.5°×0.5°)scPDSI_PM(self-calibrating palmer drought severity index based on Penman-Monteith)數(shù)據(jù)集,Mann-Kendall法、小波分析和旋轉(zhuǎn)經(jīng)驗(yàn)正交函數(shù)(rotated empirical orthogonal function,REOF)等方法對中國氣象干旱時(shí)空演變分析,并探討了干旱與大氣環(huán)流因子的關(guān)系,得到如下結(jié)論。

    1)在1961-2009年間,整體上中國年和四季變濕趨勢均顯著(P<0.05),且均于20世紀(jì)70年代初期發(fā)生了突變。周期分析表明,年和四季的周期基本一致,年干濕變化存在2.6、4.4和6.2 a的周期,其中主周期為6.2 a,春、秋、冬季有2.6和6.2 a的周期,主周期為6.2 a,夏季主周期為4.4a;年和四季均存在2~8 a的振蕩周期。

    2)根據(jù)前8個(gè)REOF空間模態(tài),將中國劃分為8個(gè)干濕特征區(qū):西北地區(qū)、東北—內(nèi)蒙古高原地區(qū)、黃淮平原地區(qū)、青藏高原北部地區(qū)、中部地區(qū)、青藏高原南部地區(qū)、大興安嶺地區(qū)、東南地區(qū)。REOF分區(qū)結(jié)果較好表現(xiàn)了地域特征,客觀反映了中國干旱分區(qū)實(shí)際情況。各分區(qū)年干濕時(shí)間演變趨勢分析表明:西北地區(qū)、青藏高原南部地區(qū)、大興安嶺地區(qū)和東南地區(qū)呈變濕趨勢;其中前兩個(gè)區(qū)變濕趨勢顯著(P<0.05);東北—內(nèi)蒙古高原區(qū)、青藏高原北部地區(qū)、中部地區(qū)和黃淮平原區(qū)呈現(xiàn)出變干趨勢,其中前3個(gè)區(qū)年變干趨勢顯著;8個(gè)分區(qū)所處地理位置不同,受各種氣候因素影響存在差異,但干濕變化普遍存在2~9 a的振蕩周期。

    3)極渦指數(shù)與中國大部分地區(qū)氣象干旱有較好相關(guān)性;西北地區(qū)干濕變化可能受極渦、印度洋偶極子(indian ocean dipole,IOD)和太平洋濤動(pacific oscillation,PDO)影響;東北—內(nèi)蒙古高原區(qū)、青藏高原和中部地區(qū)(包括四川盆地、漢中盆地和黃土高原等地區(qū))干濕變化受IOD影響顯著;PDO與黃淮海平原地區(qū)干旱呈負(fù)相關(guān)關(guān)系。

    [參考文獻(xiàn)]

    [1] Ashok K M, Vijay P S. Drought modeling-a review[J]. Journal of Hydrology, 2011, 403: 157-175.

    [2] Steven M Q, Timothy N P. An evaluation of agricultural drought indices for the Canadian prairies[J]. Agricultural and Forest Meteorology, 2003, 118: 49-62.

    [3] Liu Zhiyong, Zhou Ping, Zhang Fangqiu, et al. Spatiotemporal characteristics of dryness/wetness conditions across Qinghai Province, Northwest China[J]. Agricultural and Forest Meteorology, 2013(182/183): 101-108.

    [4] Chen Hupo, Sun Jianqi. Changes in drought characteristics over China using the standardized precipitation evapotranspiration index[J]. Journal of Climate, 2015, 28(13): 5430-5447.

    [5] Ashok K M, Vijay P S. A review of drought concepts[J]. Journal of Hydrology, 2010, 391(1/2): 204-216.

    [6] Yu Meixiu, Li Qiongfang, Lu Guobin, et al. Investigation into the impact of the Gezhouba and the Three Gores Reservoirs on the flow regime of the Yangtze River[J]. Journal of Hydrologic Engineering, 2013, 18(9): 1098-1106.

    [7] Barriopedro D, Gouveia C M, Trigo R M, et al. The 2009/10 drought in China: possible causes and impacts on vegetation[J]. Journal of Hydrometeorology, 2012, 13: 1251-1267.

    [8] Kevin E T, Dai Aiguo, Gerard VDS, et al. Global warming and changes in drought[J]. Nature climate change, 2014, 4: 17-22.

    [9] Dai Aiguo. Increasing drought under global warming in observations and models[J]. Nature climate change, 2013, 3: 52-58.

    [10] Justin S, Eric FW, Michael LR. Little change in global drought over the past 60 years[J]. Nature Climate Change,2012, 491: 435-441.

    [11] Wang Wen, Zhu Ye, Xu Rengui, et al. Drought severity change in China during 1961-2012 indicated by SPI and SPEI[J]. Natural Hazards, 2015, 75: 2437-2451.

    [12] Zou Xukai, Zhai Panmao, Zhang Qiang. Variations in droughts over China:1951-2003[J]. Geophysical Research Letters, 2005, 32(4): 2437-2451.

    [13] 徐仁貴,王文,汪小菊. 1961-2012年中國地表干濕狀況時(shí)空變化特征[J]. 中國農(nóng)村水利水電,2014,8:30-35. Xu Rengui, Wang wen, Wang Xiaoju. Spatio-temporal changes of surface humidity over China during 1961-2012[J]. China Rural Water Conservancy and Hydropower, 2014, 8: 30-35. (in Chinese with English abstract)

    [14] 翟盤茂,鄒旭愷. 1951-2003年中國氣溫和降水變化及其對干旱的影響[J]. 氣候變化研究進(jìn)展,2005,1(1):16-18. Zhai Panmao, Zou Xukai. Changes in temperature and precipitation and their impacts on drought in China during 1951-2003[J]. Advances in Climate Change Research, 2005,1(1): 16-18. (in Chinese with English abstract)

    [15] Li Jinbao, Edward R, Cook RD, et al. Moisture variability across China and Mongolia: 1951-2005[J]. Climate Dynamics, 2009, 32: 1173-1186.

    [16] Wang AiHui, Dennis PL, Justin S. Soil moisture drought in China, 1950-2006[J]. Journal of Climate, 2011, 24: 3257-3271.

    [17] Yu Meixiu, Li Qiongfang, Michael JH, et al. Are droughts becoming more frequent or severe in China based on the standardized precipitation evapotranspiration Index: 1951-2010?[J]. International Journal of Climatology, 2014, 34: 545-558.

    [18] 葉敏,錢忠華,吳永萍. 中國旱澇時(shí)空分布特征分析[J]. 物理學(xué)報(bào),2013,62(13):139-203. Ye Min, Qian Zhonghua, Wu Yongping. Spatiotemporal evolution of the droughts and floods over China[J]. Acta Physica Sinica, 2013, 62(13): 139-203. (in Chinese with English abstract)

    [19] 陳方藻,劉江,李茂松. 60年來中國農(nóng)業(yè)干旱時(shí)空演替規(guī)律研究[J]. 西南師范大學(xué)學(xué)報(bào),2011,36(4):111-114. Chen Fangzao, Liu Jiang,Li Mao-song. Researches on spatial and temporal succession law of agricultural drought in the past 60 years in China[J]. Journal of Southwest China Normal University, 2011, 36(4): 111-114. (in Chinese with English abstract)

    [20] 黃晚華,楊曉光,李茂松,等. 基于標(biāo)準(zhǔn)化降水指數(shù)的中國南方季節(jié)性干旱近58a演變特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(7):50-59. Huang Wanhua, Yang Xiaoguang, Li Maosong, et al. Evolution characteristics of seasonal drought in the south of China during the past 58 years based on standardized precipitation index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010,26(7): 50-59. (in Chinese with English abstract)

    [21] 中國氣象局. 中國災(zāi)害性天氣氣候圖集(1961-2006)[M].北京:氣象出版社,2007:32-40.

    [22] Schrier GVD, Barichivich J, Briffa KR, et al. A scPDSI-based global data set of dry and wet spells for 1901-2009[J]. Journal of Geophysical Research, 2013, 118: 4025-4048.

    [23] Palmer WC. Meteorological drought[R]. US. Weather Bureau Research Paper No. 45, 1965.

    [24] Wan Huaijung, Chen Yaning, Pan Yingping, et al. Spatial and temporal variability of drought in the arid region of China and its relationships to teleconnection indices[J]. Journal of Hydrology, 2015, 523: 283-296.

    [25] Lu Liu, Yang Hong, Christopher NB, et al. Hydroclimatological drought analyses and projections using meteorological and hydrological drought indices: a case study in Blue River Basin, Oklahoma[J]. Water Resources Manage,2012, 26: 2761-2779.

    [26] Nathan W, Steve G, Michael JH. A self-calibrating palmer drought severity index[J]. Journal of Climate, 2004, 17: 2335-2351.

    [27] Kendall, MG, Rank Correlation Methods, 4th ed. Charles Griffin, London, 1975.

    [28] Henry BM. Nonparametric tests against trend[J]. Econometric Society, 1945, 3(3): 245-259.

    [29] Christopher T, Gilbert PC. A practical guide to wavelet analysis[J]. American Meteorological Society, 1998, 79: 61-78.

    [30] Tao Lian, Chen Dake. An evaluation of rotated EOF analysis and its application to tropical Pacific SST variability[J]. Journal of Climate, 2012, 5(15): 5361-5373.

    [31] 魏鳳英. 現(xiàn)代氣候統(tǒng)計(jì)診斷預(yù)測技術(shù)[M]. 北京:氣象出版社,1999.

    [32] Lee H F, Zhang D D. Relationship between NAO and drought disasters in northwestern China in the last millennium[J]. Journal of Arid Environments, 2011,75(11): 1114-1120.

    [33] Li Yong, Lu Riyu, He Jinhai. Several climate factors influencing the winter temperature over china[J]. Chinese Journal of Atmospheric Sciences, 2007, 31(3): 505-514.

    [34] Xiao Mingzhong, Zhang Qiang, Vijay PS, et al. Influences of ENSO, NAO, IOD and PDO on seasonal precipitation regimes in the Yangtze River basin, China[J]. International Journal of Climatology, 2014, 35(12): 3556-3567.

    [35] Gu Sinan, Yang Xiuqun. Variability of the northern circumpolar vortex and its association with climate anomaly in China[J]. Scientia Meteorologica Sinica, 2006, 26(2): 135-142.

    [36] 張恒德,陸維松,高守亭,等. 北極渦活動對我國同期及后期氣溫的影響[J]. 南京氣象學(xué)院學(xué)報(bào),2006,29(4):507-516. Zhang Hengde, Lu Weisong, Gao Shouting, et al. Influence of the north polar vortex activity on the contemporaneous and subsequent air temperature in China[J]. Journal of Nanjing Insitute of Meteorology, 2006, 29(4): 507-516. (in Chinese with English abstract)

    [37] 張恒德,高守亭,張友姝. 北極渦年代際變化及其與我國春季降水的關(guān)系[J]. 氣候與環(huán)境研究,2006,11(5):593-604. Zhang Hengde, Gao Shouting, Zhang Youshu. The interdecadal variation of north polar vortex and itsrelationships with spring precipitation in China[J]. Climatic and Environmental Research, 2006, 11(5): 593-604. (in Chinese with English abstract)

    [38] 肖子牛,晏紅明,李崇銀. 印度洋地區(qū)異常海溫的偶極振蕩與中國降水及溫度的關(guān)系[J]. 熱帶氣象學(xué)報(bào),2002,18(4):335-344. Xiao Ziniu, Yan Hongming, Li Chongyin. The relationship between indian ocean SSTA dipole index and the precipitation and temperature over China[J]. Journal of Tropical Meteorology, 2002, 18(4): 335-344. (in Chinese with English abstract)

    [39] 楊秋明. 南印度洋副熱帶偶極子型海溫異常與全球環(huán)流和我國降水變化的關(guān)系[J]. 海洋學(xué)報(bào),2006,28(3):47-56. Yang Qiuming. Indian Ocean subtropical dipole and variations of global circulations and rainfall in China[J]. Acta Oceanologica Sinica, 2006, 28(3): 47-56. (in Chinese with English abstract)

    [40] 趙平,陳軍明,肖棟,等. 夏季亞洲—太平洋濤動與大氣環(huán)流和季風(fēng)降水[J]. 氣象學(xué)報(bào),2008,66(5):716-729. Zhao Ping, Chen Junmin, Xiao Dong, et al. Summer Asian-Pacific oscillation and its relation-ship with atmospheric circulation and monsoon rainfall[J]. Acta Meteorologica Sinica, 2008, 66(5): 716-729. (in Chinese with English abstract)

    [41] Chen Deliang, Xu Chong-Yu, Ren Guoyu. Comparison of the Thornthwaite method and pan data with the standard Penman-Monteith estimates of reference ET in China[J]. Climate Research, 2008, 28: 123-132.

    ·農(nóng)業(yè)信息與電氣技術(shù)·

    [42] Zheng Chaolei, Wang Quan. Spatiotemporal pattern of the global sensitivity of the reference evapotranspiration to climatic variables in recent five decades over China[J]. Stochastic Environmental Research and Risk Assessment,2015, 29: 1937-1947.

    Spatiotemporal variations analysis of meteorological drought in China based on scPDSI

    Wang Zhaoli, Li Jun, Huang Zeqin, Zhong Ruida, Chen Jiayin, Qiu Zhuohui
    (School of Civil Engineering ɑnd Trɑnsportɑtion, South Chinɑ University of Technology, Guɑngzhou 510641, Chinɑ)

    Abstract:In recent decades, frequent droughts have been detected in some places in China along with the global climate change, which cause great damage to the development of social economy. Thus, it’s high time to investigate the spatio-temporal evolution of droughts. Based on the newly published scPDSI PM (self-calibrating Palmer drought severity index based on Penman-Monteith) dataset with high spatial resolution, the annual and seasonal variations of droughts in China from 1961 to 2009 were researched. Other research methods like linear regression, Mann-Kendall method, wavelet analysis as well as REOF (rotated empirical orthogonal function) were adopted to investigate the spatio-temporal pattern of droughts in China, and find out how droughts changed under the influence of the general circulation of atmosphere, in order to provide the foundation for drought control and drought resistance. The results indicated that: 1) In 1961-2009, China as a whole became moist significantly at annual and seasonal scale with the mutations of turning drought to wetness detected in early 1970s. 2)Period analysis illustrated that the oscillation period of seasonal and annual drought remained generally unanimous, and the dominant period of annual drought was 6.2 years, while drought in spring, autumn and winter generally had the periods of 2.6 and 6.2 years with 6.2 years as the dominant period; otherwise, drought in summer presented a dominant period of 4.4 years. 3)In the light of the spatial modes disassembled from REOF, the entire country was divided into 8 regions with different characteristics of drought or wetness, which were Northwest China, Northeast - Inner Mongolia Plateau area, the Greater Khingan Range area, the northern Qinghai-Tibet Plateau area, the southern Qinghai-Tibet Plateau area, the Middle China (including Sichuan Basin, Hanzhong Basin and Loess Plateau), Huang-Huai-Hai Plain and Southeast China. It was found that the sub-regions of REOF could primely manifest the geographical features of different regions, and hence objectively testify the actual drought condition in China. Among the 8 regions, the Greater Khingan Range area, Southeast China, the southern Qinghai-Tibet Plateau area, and Northwest China were getting moist while the latter 2 areas had a significant trend (P<0.05); Northeast - Inner Mongolia Plateau area, the northern Qinghai-Tibet Plateau area, the Middle China and Huang-Huai-Hai Plain were getting dry, and except Huang-Huai-Hai Plain, the former 3 regions were prominently getting dry. Take the regions as the objects, no matter at annual or seasonal scale, Northwest China was getting moist, while the Middle China was getting arid (P<0.05). Although the geographical positions of the 8 regions resulted in the discrepancy of influence under different climate factors, the 8 regions generally had the oscillation period of 2-9 years. 4) M-K test (Mann-Kendall test) on the 8 regions showed that the seasonal MK values of different regions varied with each other, which indicated that drought-wet degree of different regions had obvious seasonal feature, and yet most regions had a tendency of getting drought in spring and getting moist in summer. 5) Good correlations were found between polar vortex index and drought in most areas. Correlation analysis indicated that droughts in Northwest China may be affected by the polar vortex, Indian Ocean Dipole (IOD) and Pacific Oscillation (PDO). Droughts in Northeast - Inner Mongolia Plateau area, the Qinghai-Tibet Plateau area and the Middle China were significantly influenced by IOD, and drought of Huang-Huai-Hai Plain was negatively correlated with PDO. Unlike the conventional PDSI index, scPDSI integrates solar radiation and air speed into account, which will be of scientific and practical importance to forecast and distinguish different drought conditions.

    Keywords:drought; meteorology; precipitation; scPDSI; REOF; spatio-temporal evolution; China

    基金項(xiàng)目:國家自然科學(xué)基金(51209095,51579105),國家科技支撐計(jì)劃(2012BAC21B0103),中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)基金(2014ZZ0027)。

    作者簡介:王兆禮,男(漢族),江蘇徐州人,博士,副教授,主要從事水文水資源研究。廣州華南理工大學(xué) 土木與交通學(xué)院,510641。

    收稿日期:2015-11-06

    修訂日期:2015-12-23

    中圖分類號:S165

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1002-6819(2016)-02-0161-08

    doi:10.11975/j.issn.1002-6819.2016.02.024

    猜你喜歡
    時(shí)空演變干旱中國
    基于距平的白城地區(qū)干旱時(shí)間分布特征分析
    臨夏地區(qū)干旱特征及干濕氣候區(qū)劃
    主要客源國對中國入境旅游市場的貢獻(xiàn)分析
    肯尼·格雷特,爵士的“中國”調(diào)子
    基于WebGIS的城市空氣質(zhì)量時(shí)空演變動態(tài)可視化系統(tǒng)研究
    夏季高溫干旱時(shí)節(jié)高山蔬菜種植管理策略
    江蘇省物流業(yè)與經(jīng)濟(jì)發(fā)展耦合協(xié)調(diào)性的時(shí)空演變特征
    小說里的世界 虛擬里的真實(shí)
    山東省工業(yè)經(jīng)濟(jì)增長的水資源效應(yīng)及時(shí)空演變分析
    英媒:“中國”成美國網(wǎng)絡(luò)威脅敏感詞
    久久天堂一区二区三区四区| 国产精品一区二区精品视频观看| 亚洲av片天天在线观看| 欧美 日韩 精品 国产| 国产在线观看jvid| 亚洲九九香蕉| bbb黄色大片| 可以免费在线观看a视频的电影网站| 母亲3免费完整高清在线观看| 久久99精品国语久久久| 国产精品一区二区在线不卡| 亚洲午夜精品一区,二区,三区| 中文字幕亚洲精品专区| 宅男免费午夜| 永久免费av网站大全| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩视频高清一区二区三区二| 女人爽到高潮嗷嗷叫在线视频| 另类精品久久| 午夜免费男女啪啪视频观看| 亚洲激情五月婷婷啪啪| 黑人猛操日本美女一级片| 国产在视频线精品| 妹子高潮喷水视频| 亚洲国产精品999| 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 午夜免费鲁丝| 亚洲精品在线美女| 亚洲人成网站在线观看播放| 久久免费观看电影| 国产人伦9x9x在线观看| 久久人人爽人人片av| 巨乳人妻的诱惑在线观看| 18禁黄网站禁片午夜丰满| 狠狠婷婷综合久久久久久88av| 亚洲国产精品999| 免费av中文字幕在线| 亚洲中文av在线| 秋霞在线观看毛片| 久久精品aⅴ一区二区三区四区| 一区在线观看完整版| 国产成人欧美| 美女扒开内裤让男人捅视频| 亚洲欧美日韩高清在线视频 | 日本av手机在线免费观看| 国产片特级美女逼逼视频| 久久精品国产a三级三级三级| 国产成人一区二区在线| 久久精品成人免费网站| 免费av中文字幕在线| 日本欧美国产在线视频| 中文精品一卡2卡3卡4更新| 亚洲国产精品一区二区三区在线| 欧美日韩视频精品一区| 在现免费观看毛片| 亚洲美女黄色视频免费看| 亚洲熟女精品中文字幕| svipshipincom国产片| 一本久久精品| 亚洲成人免费av在线播放| 少妇被粗大的猛进出69影院| 亚洲av欧美aⅴ国产| 成人国语在线视频| 欧美国产精品va在线观看不卡| 韩国精品一区二区三区| videosex国产| 又大又黄又爽视频免费| 久久国产亚洲av麻豆专区| 丰满饥渴人妻一区二区三| 精品人妻在线不人妻| 香蕉丝袜av| 欧美日韩福利视频一区二区| 国产av国产精品国产| 天堂俺去俺来也www色官网| 亚洲精品久久成人aⅴ小说| 久久人妻福利社区极品人妻图片 | 国产精品国产三级国产专区5o| 欧美97在线视频| 亚洲九九香蕉| 亚洲成国产人片在线观看| 热re99久久国产66热| 欧美日韩精品网址| 国产麻豆69| 日日夜夜操网爽| 巨乳人妻的诱惑在线观看| 精品高清国产在线一区| 国产精品国产三级专区第一集| 深夜精品福利| 黑人巨大精品欧美一区二区蜜桃| 又紧又爽又黄一区二区| netflix在线观看网站| 国产成人91sexporn| 少妇精品久久久久久久| 亚洲精品久久成人aⅴ小说| 另类亚洲欧美激情| 汤姆久久久久久久影院中文字幕| 国产男女内射视频| 精品福利观看| 亚洲av美国av| 99精品久久久久人妻精品| 狠狠婷婷综合久久久久久88av| 韩国高清视频一区二区三区| 黄色毛片三级朝国网站| 亚洲av成人不卡在线观看播放网 | 人人妻人人澡人人爽人人夜夜| 美女主播在线视频| 亚洲国产欧美网| 婷婷色av中文字幕| 国产免费现黄频在线看| 一本色道久久久久久精品综合| 大片免费播放器 马上看| 亚洲 国产 在线| www.av在线官网国产| 国产爽快片一区二区三区| 99久久精品国产亚洲精品| 欧美97在线视频| 亚洲精品国产区一区二| xxxhd国产人妻xxx| 少妇人妻 视频| 亚洲色图 男人天堂 中文字幕| 色网站视频免费| 国产日韩欧美亚洲二区| 亚洲精品第二区| 国产福利在线免费观看视频| 自线自在国产av| 国产一区二区 视频在线| 日本色播在线视频| 久久久亚洲精品成人影院| 美女高潮到喷水免费观看| 两个人免费观看高清视频| 精品亚洲乱码少妇综合久久| 国产视频一区二区在线看| 久久精品亚洲熟妇少妇任你| 精品高清国产在线一区| 亚洲人成电影免费在线| av天堂在线播放| 操美女的视频在线观看| 国产高清视频在线播放一区 | 黄片小视频在线播放| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 国产麻豆69| 一级a爱视频在线免费观看| 亚洲情色 制服丝袜| 日本av免费视频播放| 成人亚洲欧美一区二区av| kizo精华| 别揉我奶头~嗯~啊~动态视频 | 真人做人爱边吃奶动态| 极品少妇高潮喷水抽搐| 母亲3免费完整高清在线观看| 18在线观看网站| 丰满饥渴人妻一区二区三| 国产精品一二三区在线看| 少妇被粗大的猛进出69影院| 少妇粗大呻吟视频| 婷婷色综合www| 高清av免费在线| 天天躁夜夜躁狠狠久久av| √禁漫天堂资源中文www| 麻豆乱淫一区二区| 色网站视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产黄色视频一区二区在线观看| www.999成人在线观看| 欧美日韩福利视频一区二区| 在线 av 中文字幕| 亚洲三区欧美一区| 丁香六月欧美| 一级片'在线观看视频| 日韩 亚洲 欧美在线| 久久天堂一区二区三区四区| 日韩一区二区三区影片| 老司机深夜福利视频在线观看 | 精品福利观看| 亚洲国产精品国产精品| 高清av免费在线| 国产一区二区三区综合在线观看| 国产欧美日韩精品亚洲av| 欧美激情高清一区二区三区| 一级片免费观看大全| 一区二区三区精品91| 大香蕉久久成人网| 大码成人一级视频| 9191精品国产免费久久| 一区二区三区激情视频| 亚洲国产看品久久| 亚洲黑人精品在线| 91老司机精品| 久久精品熟女亚洲av麻豆精品| 精品免费久久久久久久清纯 | 精品亚洲成a人片在线观看| 纵有疾风起免费观看全集完整版| 一区二区三区四区激情视频| 国产真人三级小视频在线观看| 日韩大码丰满熟妇| 好男人视频免费观看在线| 精品亚洲成国产av| 成人手机av| 国产1区2区3区精品| 国产亚洲精品第一综合不卡| 日韩,欧美,国产一区二区三区| 国产精品久久久av美女十八| 国产主播在线观看一区二区 | 久9热在线精品视频| 男女之事视频高清在线观看 | 丝袜美腿诱惑在线| 久久久欧美国产精品| 性色av一级| 国产免费视频播放在线视频| 亚洲精品久久久久久婷婷小说| 国产伦人伦偷精品视频| 只有这里有精品99| 精品久久蜜臀av无| 亚洲,一卡二卡三卡| 另类精品久久| 国产欧美日韩精品亚洲av| 青草久久国产| 亚洲欧美日韩高清在线视频 | 日本wwww免费看| 欧美国产精品一级二级三级| 国产成人a∨麻豆精品| 老司机影院成人| 下体分泌物呈黄色| 一边摸一边做爽爽视频免费| 国产一区有黄有色的免费视频| 婷婷色av中文字幕| 男人爽女人下面视频在线观看| 看免费av毛片| 亚洲综合色网址| 91老司机精品| 国产精品一区二区在线不卡| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 黄片小视频在线播放| 香蕉国产在线看| 十八禁网站网址无遮挡| 一本大道久久a久久精品| 欧美日韩av久久| 一本久久精品| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久人人做人人爽| 性色av乱码一区二区三区2| 人人妻人人澡人人看| 午夜免费鲁丝| 免费看av在线观看网站| 亚洲国产精品999| 麻豆国产av国片精品| 一区二区三区四区激情视频| 精品亚洲乱码少妇综合久久| 精品福利永久在线观看| 免费av中文字幕在线| 脱女人内裤的视频| 亚洲人成网站在线观看播放| 两个人免费观看高清视频| 纵有疾风起免费观看全集完整版| 欧美国产精品一级二级三级| 国产精品一二三区在线看| 日本黄色日本黄色录像| 七月丁香在线播放| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 亚洲免费av在线视频| 婷婷色av中文字幕| 久久久亚洲精品成人影院| 欧美乱码精品一区二区三区| 国产在线免费精品| 久久久久久久久免费视频了| 亚洲精品一二三| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| 免费看十八禁软件| 亚洲一区中文字幕在线| 你懂的网址亚洲精品在线观看| 午夜免费鲁丝| 尾随美女入室| 久久久精品94久久精品| 国产成人系列免费观看| 国产精品亚洲av一区麻豆| 美女视频免费永久观看网站| 亚洲美女黄色视频免费看| 老司机在亚洲福利影院| 日本a在线网址| 男男h啪啪无遮挡| 一本一本久久a久久精品综合妖精| 久久久国产一区二区| 亚洲av在线观看美女高潮| 亚洲五月婷婷丁香| 国产日韩欧美视频二区| 亚洲国产欧美一区二区综合| 美女国产高潮福利片在线看| 中文字幕人妻丝袜制服| netflix在线观看网站| 国产一区二区三区综合在线观看| 我的亚洲天堂| 中文字幕人妻丝袜一区二区| 9色porny在线观看| 亚洲国产精品999| 久久久久久人人人人人| 老司机影院毛片| 亚洲综合色网址| 日韩一区二区三区影片| 老鸭窝网址在线观看| 超碰成人久久| 久久中文字幕一级| 91九色精品人成在线观看| 深夜精品福利| 精品一区二区三区av网在线观看 | 精品少妇久久久久久888优播| 午夜福利一区二区在线看| 成人国语在线视频| 日本91视频免费播放| 欧美日韩一级在线毛片| 老司机靠b影院| 麻豆av在线久日| 欧美精品啪啪一区二区三区 | 伊人亚洲综合成人网| 一本久久精品| 成人亚洲欧美一区二区av| 啦啦啦视频在线资源免费观看| 18禁国产床啪视频网站| 国产日韩一区二区三区精品不卡| 久久久欧美国产精品| 美女扒开内裤让男人捅视频| 欧美在线黄色| 日韩中文字幕视频在线看片| 天天影视国产精品| 又粗又硬又长又爽又黄的视频| 国产亚洲av片在线观看秒播厂| 久久精品久久久久久噜噜老黄| 天堂俺去俺来也www色官网| 1024视频免费在线观看| 97精品久久久久久久久久精品| 国产免费福利视频在线观看| 女性被躁到高潮视频| 99国产综合亚洲精品| 亚洲,一卡二卡三卡| 中文字幕人妻丝袜一区二区| 久久久久久久精品精品| 91精品三级在线观看| 国产人伦9x9x在线观看| 国产精品国产av在线观看| 亚洲综合色网址| 午夜影院在线不卡| 午夜激情久久久久久久| 欧美97在线视频| 大码成人一级视频| 日韩精品免费视频一区二区三区| 免费久久久久久久精品成人欧美视频| 一区二区日韩欧美中文字幕| 亚洲欧洲日产国产| 美女中出高潮动态图| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 日本wwww免费看| 蜜桃国产av成人99| 国产福利在线免费观看视频| 久久久精品94久久精品| 91国产中文字幕| 天堂8中文在线网| 嫩草影视91久久| 最近最新中文字幕大全免费视频 | 亚洲国产欧美日韩在线播放| 欧美精品一区二区大全| 欧美日韩国产mv在线观看视频| 精品人妻在线不人妻| 亚洲欧洲国产日韩| 国产日韩欧美视频二区| 国产成人精品久久二区二区91| 免费在线观看日本一区| 人妻 亚洲 视频| 少妇人妻久久综合中文| 欧美日韩亚洲高清精品| 成人亚洲欧美一区二区av| 成人国语在线视频| 日韩 亚洲 欧美在线| 久久鲁丝午夜福利片| 操出白浆在线播放| 高清不卡的av网站| 亚洲成人免费av在线播放| 如日韩欧美国产精品一区二区三区| 亚洲欧美一区二区三区久久| 午夜免费男女啪啪视频观看| av片东京热男人的天堂| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 日本av免费视频播放| 最近手机中文字幕大全| 国产成人欧美| 免费少妇av软件| 久久av网站| 人妻一区二区av| 国产在线一区二区三区精| 久久鲁丝午夜福利片| 久久久久久久久久久久大奶| 高清视频免费观看一区二区| 男女无遮挡免费网站观看| 亚洲精品一区蜜桃| 男女午夜视频在线观看| 日韩 亚洲 欧美在线| av欧美777| 超色免费av| 精品少妇内射三级| 欧美97在线视频| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 精品一品国产午夜福利视频| 精品久久久久久电影网| 亚洲欧美一区二区三区国产| www.999成人在线观看| 成年女人毛片免费观看观看9 | 久久天堂一区二区三区四区| 波野结衣二区三区在线| 日韩大码丰满熟妇| 国产精品三级大全| 一区二区三区四区激情视频| 国产精品久久久久久精品古装| 日韩,欧美,国产一区二区三区| 欧美精品一区二区免费开放| 久久 成人 亚洲| 午夜福利乱码中文字幕| 欧美黄色淫秽网站| 三上悠亚av全集在线观看| 国产成人精品无人区| 脱女人内裤的视频| 一级a爱视频在线免费观看| 精品国产一区二区三区四区第35| 成年人黄色毛片网站| 久久精品亚洲熟妇少妇任你| 久久毛片免费看一区二区三区| 中文字幕色久视频| 久久国产精品影院| 国产精品99久久99久久久不卡| 亚洲,欧美精品.| 一本色道久久久久久精品综合| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 丝袜在线中文字幕| 91国产中文字幕| 国产精品国产三级国产专区5o| 国产成人精品久久二区二区91| 我要看黄色一级片免费的| 午夜福利,免费看| 丰满饥渴人妻一区二区三| 永久免费av网站大全| 一本久久精品| 天天添夜夜摸| 午夜免费鲁丝| 欧美精品高潮呻吟av久久| 欧美人与性动交α欧美精品济南到| 曰老女人黄片| 大型av网站在线播放| 欧美性长视频在线观看| e午夜精品久久久久久久| 男人爽女人下面视频在线观看| 丝袜在线中文字幕| 91麻豆精品激情在线观看国产 | 狠狠精品人妻久久久久久综合| 国产高清不卡午夜福利| 18禁裸乳无遮挡动漫免费视频| 777久久人妻少妇嫩草av网站| 亚洲国产日韩一区二区| 日韩制服骚丝袜av| 日韩免费高清中文字幕av| 国产主播在线观看一区二区 | 丝袜脚勾引网站| 亚洲国产精品999| av片东京热男人的天堂| 777久久人妻少妇嫩草av网站| 亚洲国产精品一区二区三区在线| 美国免费a级毛片| 丝袜美腿诱惑在线| 男男h啪啪无遮挡| 亚洲成色77777| 精品卡一卡二卡四卡免费| 考比视频在线观看| 精品一区在线观看国产| 青春草视频在线免费观看| 久久久国产一区二区| 日韩视频在线欧美| 欧美日韩一级在线毛片| 美女国产高潮福利片在线看| 国产精品av久久久久免费| 国产成人av激情在线播放| 男女无遮挡免费网站观看| 国产色视频综合| 在现免费观看毛片| 国产精品国产av在线观看| 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 一边摸一边抽搐一进一出视频| 三上悠亚av全集在线观看| 亚洲黑人精品在线| 夜夜骑夜夜射夜夜干| 看免费成人av毛片| 久久人人爽av亚洲精品天堂| 免费在线观看完整版高清| 亚洲国产av新网站| 日韩制服丝袜自拍偷拍| 99热国产这里只有精品6| 国产爽快片一区二区三区| 69精品国产乱码久久久| 中文字幕制服av| 少妇人妻 视频| 乱人伦中国视频| 精品久久久久久电影网| 飞空精品影院首页| 后天国语完整版免费观看| 老司机深夜福利视频在线观看 | 欧美变态另类bdsm刘玥| 亚洲伊人久久精品综合| 丰满饥渴人妻一区二区三| 久久天躁狠狠躁夜夜2o2o | 蜜桃在线观看..| 亚洲成国产人片在线观看| 老汉色∧v一级毛片| 满18在线观看网站| 黄色毛片三级朝国网站| 久久天堂一区二区三区四区| 爱豆传媒免费全集在线观看| 久久久久久免费高清国产稀缺| 91国产中文字幕| 啦啦啦啦在线视频资源| 免费在线观看黄色视频的| 亚洲图色成人| 国产一区二区在线观看av| 极品少妇高潮喷水抽搐| xxx大片免费视频| 久久久精品免费免费高清| 久久热在线av| 久久鲁丝午夜福利片| 欧美精品av麻豆av| 尾随美女入室| 黑人欧美特级aaaaaa片| 蜜桃在线观看..| 伦理电影免费视频| 日韩精品免费视频一区二区三区| 人人妻人人澡人人看| 狠狠精品人妻久久久久久综合| 黄色片一级片一级黄色片| 亚洲专区中文字幕在线| 国产成人av教育| 18禁观看日本| 菩萨蛮人人尽说江南好唐韦庄| 狠狠精品人妻久久久久久综合| 香蕉国产在线看| 婷婷成人精品国产| 在线观看免费日韩欧美大片| 成年人午夜在线观看视频| 波野结衣二区三区在线| 亚洲 欧美一区二区三区| 国产一区有黄有色的免费视频| 亚洲欧洲国产日韩| 女人被躁到高潮嗷嗷叫费观| 超色免费av| 妹子高潮喷水视频| 男的添女的下面高潮视频| 国产精品久久久久久精品电影小说| 成年美女黄网站色视频大全免费| 高清不卡的av网站| 欧美激情 高清一区二区三区| 丰满迷人的少妇在线观看| 午夜免费成人在线视频| 操出白浆在线播放| 亚洲图色成人| 国产成人av教育| 18禁观看日本| 亚洲欧美成人综合另类久久久| 日本91视频免费播放| 欧美日韩亚洲国产一区二区在线观看 | 欧美激情极品国产一区二区三区| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜制服| 亚洲一区中文字幕在线| 亚洲av国产av综合av卡| 久久精品久久久久久久性| 亚洲熟女精品中文字幕| 亚洲国产精品国产精品| 制服人妻中文乱码| 天堂中文最新版在线下载| 午夜视频精品福利| 久久久精品94久久精品| 欧美日韩亚洲综合一区二区三区_| 免费看av在线观看网站| 又粗又硬又长又爽又黄的视频| 女人久久www免费人成看片| 国产黄色视频一区二区在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲图色成人| 又粗又硬又长又爽又黄的视频| 久久亚洲国产成人精品v| 亚洲九九香蕉| 久久 成人 亚洲| av有码第一页| 一级毛片黄色毛片免费观看视频| 国产成人精品久久二区二区免费| 久久中文字幕一级| 亚洲欧美日韩另类电影网站| 久久精品久久久久久噜噜老黄| 国产亚洲av片在线观看秒播厂| 熟女少妇亚洲综合色aaa.| 国产91精品成人一区二区三区 | 九草在线视频观看| 在现免费观看毛片| 欧美成狂野欧美在线观看| 99久久人妻综合| 亚洲欧洲日产国产| 国产一区二区激情短视频 | 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| 两个人免费观看高清视频| 啦啦啦视频在线资源免费观看| 久久人妻熟女aⅴ| 久久久久国产精品人妻一区二区|