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

    基于SPI的云南省多尺度干旱時(shí)空演變特征識別

    2023-05-09 11:53:52陳正發(fā)李靖相彪段青松李淑芳王樹仿
    灌溉排水學(xué)報(bào) 2023年4期
    關(guān)鍵詞:時(shí)空云南省站點(diǎn)

    陳正發(fā),李靖,相彪,段青松,李淑芳,王樹仿

    基于的云南省多尺度干旱時(shí)空演變特征識別

    陳正發(fā)1,2,李靖1,2,相彪1*,段青松1,2,李淑芳1,王樹仿1,2

    (1.云南農(nóng)業(yè)大學(xué) 水利學(xué)院,昆明 650201;2.云南省農(nóng)業(yè)節(jié)水工程技術(shù)研究中心,昆明 650201)

    【目的】識別云南省多尺度干旱的時(shí)空演變特征。【方法】基于云南省36個(gè)氣象站31 a的氣象數(shù)據(jù),基于標(biāo)準(zhǔn)化降水指數(shù)(),結(jié)合GIS空間分析、非參數(shù)Mann-kendall趨勢檢驗(yàn),識別云南省多尺度氣象干旱時(shí)空分布及演變特征。【結(jié)果】①在年尺度、季節(jié)尺度均呈波動(dòng)變化趨勢,除春季略有增長外,夏季、秋季、冬季和年尺度的均表現(xiàn)為減小趨勢,秋季變幅最大,冬季變幅最小。②年尺度干旱頻率分布在32.26%~50.00%之間,主導(dǎo)的干旱等級為輕微干旱和極端干旱,干旱頻率高值區(qū)主要分布在文山、紅河、昆明、楚雄、德宏。③春、夏、秋、冬4個(gè)季節(jié)的干旱頻率分別為39.86%、39.83%、38.43%、41.33%,盡管季節(jié)間干旱頻率差異較小,但不同季節(jié)主導(dǎo)性干旱等級差異較大。除春季外,各季節(jié)的主導(dǎo)性干旱等級均包含極端干旱。④年尺度上,除5個(gè)站點(diǎn)干旱呈增強(qiáng)趨勢外,其余站點(diǎn)均無顯著變化趨勢;春季大部分站點(diǎn)干旱演變呈減小趨勢,而夏季、秋季、冬季則分別有11、6、2個(gè)站點(diǎn)呈顯著增強(qiáng)趨勢?!窘Y(jié)論】云南省年尺度和季節(jié)尺度干旱頻率分布及演變趨勢呈明顯的時(shí)空分異特征,區(qū)域干旱以季節(jié)性干旱為主,其中春旱、冬旱的干旱頻率和等級總體較高。

    時(shí)空演變;氣象干旱;季節(jié)性干旱;標(biāo)準(zhǔn)化降水指數(shù);云南省

    0 引言

    【研究意義】干旱作為一種具有復(fù)雜性和一定規(guī)律性的氣候?yàn)?zāi)害,對農(nóng)業(yè)、水資源、生態(tài)環(huán)境均會造成負(fù)面影響[1]。全球每年因干旱造成的經(jīng)濟(jì)損失高達(dá)60億~80億美元,遠(yuǎn)高于其他氣象災(zāi)害。伴隨著全球氣候變暖,干旱表現(xiàn)出進(jìn)一步加重的趨勢[2]。干旱成因復(fù)雜,目前針對干旱的概念還沒有統(tǒng)一的定義。Palmer[3]認(rèn)為干旱是一個(gè)持續(xù)、異常的水分虧缺事件。世界氣象組織將干旱定義為在較大范圍內(nèi)相對長期平均水平而言的降水減少,從而導(dǎo)致自然生態(tài)系統(tǒng)和雨養(yǎng)農(nóng)業(yè)生產(chǎn)力下降的過程[4]。美國氣象學(xué)會將干旱劃分為氣象干旱、農(nóng)業(yè)干旱、水文干旱、社會經(jīng)濟(jì)干旱4種類型[5]。定量識別區(qū)域干旱時(shí)空分布特征及其演變規(guī)律,對區(qū)域干旱防控和制定科學(xué)的水土資源管理措施具有重要的指導(dǎo)意義。

    【研究進(jìn)展】目前,國內(nèi)外較為成熟的干旱評估方法主要有降水距平百分比[6]、Palmer指標(biāo)[7]、指數(shù)[8]、指數(shù)[9]、標(biāo)準(zhǔn)化降水指數(shù)()[10]、標(biāo)準(zhǔn)化降水蒸發(fā)指數(shù)()[11]等。近年來,部分學(xué)者對灰色模糊聚類理論[12]、模糊物元模型[13]等方法進(jìn)行了研究,取得了一定研究進(jìn)展。是McKee等[14]在評估美國科羅拉干旱狀況時(shí)提出的一種氣象干旱評估方法,該方法通過概率密度函數(shù)求解累積概率,再將累積概率進(jìn)行標(biāo)準(zhǔn)化而得到,消除了降水的時(shí)空分布差異,在各區(qū)域和多時(shí)間尺度上均能有效地反映旱澇狀況。Merabti等[15]研究表明,基于的干旱評估結(jié)果與具有較好的一致性。Xu等[16]研究表明,對于干旱區(qū)而言,比更適合作為干旱評估指標(biāo)。袁文平[17]等研究表明,計(jì)算簡便,資料容易獲取,且計(jì)算結(jié)果與指數(shù)具有極好的一致性。

    云南省作為我國西南季風(fēng)和東南季風(fēng)的交匯區(qū),降水時(shí)空分布不均導(dǎo)致的季節(jié)性干旱災(zāi)害頻發(fā),干旱導(dǎo)致的農(nóng)業(yè)減產(chǎn)已成為區(qū)域最大的致災(zāi)因素。此外,云南省近50年的氣溫變化與全球、北半球和中國的變化趨勢基本一致,全省降水日數(shù)和雨季降水量減少,但暴雨、大暴雨事件的發(fā)生頻率卻有所上升[18]。降水分配不均和降水量減小將進(jìn)一步加劇區(qū)域氣象干旱程度。近年來,眾多學(xué)者對西南地區(qū)干旱時(shí)空變化及成因進(jìn)行了研究。Xu等[16]采用對中國干旱時(shí)空分布進(jìn)行了研究,得出云南等中國西南地區(qū)的干旱事件有加劇的趨勢。Kim等[19]基于和相結(jié)合的研究表明,云南省約29.4%的地區(qū)屬于干旱易發(fā)區(qū)。【切入點(diǎn)】截至目前,針對云南省干旱時(shí)空演變規(guī)律的研究多側(cè)重于干旱分布特征描述,對多時(shí)空尺度干旱分布、干旱強(qiáng)度變化、連續(xù)性干旱、干旱年際變化及變異特征等方面的研究相對較少?!緮M解決的關(guān)鍵問題】為此,本研究選取勻分布于云南省的36個(gè)氣象觀測站近31年的降水觀測數(shù)據(jù),基于指數(shù)、GIS空間分析、非參數(shù)Mann-Kendall趨勢檢驗(yàn),識別云南省年尺度干旱和季節(jié)性干旱的時(shí)空分布規(guī)律及演變特征,為區(qū)域干旱災(zāi)害防控和農(nóng)業(yè)水土資源管理提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    云南省地處中國西南邊陲,處于東亞季風(fēng)和南亞季風(fēng)交匯區(qū)域。該省地形由高原山地向喀斯特地貌逐漸演變,84%的區(qū)域面積為山區(qū),丘陵區(qū)僅占該省總面積的10%,壩區(qū)面積僅為該省總面積的6%。該省降水充沛,河流眾多,但在時(shí)空分布上嚴(yán)重不均,導(dǎo)致干旱災(zāi)害頻發(fā),春旱和秋旱是云南省最為頻發(fā)的干旱災(zāi)害。研究區(qū)高程及降水觀測站空間分布如圖1所示。

    圖1 云南省高程及降水觀測站空間分布

    1.2 研究方法

    1.2.1 干旱表征方法

    本研究采用McKee等[14]提出的指數(shù)識別云南省不同尺度干旱的時(shí)空演變特征,計(jì)算方法詳見文獻(xiàn)[10]。依據(jù)劃分干旱強(qiáng)度等級,干旱等級劃分標(biāo)準(zhǔn)詳見表1。具有多時(shí)間尺度特征,3月(季節(jié)尺度)、12月(年尺度)時(shí)間尺度的應(yīng)用最多,且季節(jié)尺度的3和年尺度的12更能反映干旱年內(nèi)和年際變化[20]。因此,本研究采用年尺度12、季節(jié)尺度3來評估云南省年尺度、季節(jié)尺度的干旱時(shí)空分布及演變特征。

    表1 SPI干旱等級劃分標(biāo)準(zhǔn)

    1.2.2 干旱頻率

    P表示干旱頻率,其數(shù)值等于干旱發(fā)生頻數(shù)與統(tǒng)計(jì)序列樣本數(shù)之比,計(jì)算式為:

    式中:為不同站點(diǎn)的代號;為對應(yīng)站點(diǎn)發(fā)生干旱的頻數(shù);為統(tǒng)計(jì)序列樣本數(shù)。

    1.2.3 干旱范圍

    P表示干旱范圍,其計(jì)算式為:

    式中:為不同年份的代號;為發(fā)生干旱的站點(diǎn)數(shù);為研究區(qū)內(nèi)的總站點(diǎn)數(shù)。

    根據(jù)P數(shù)值劃分干旱范圍等級。當(dāng)P≥50%時(shí),為全域性干旱;當(dāng)33%≤P<50%時(shí),為區(qū)域性干旱;當(dāng)25%≤P<33%時(shí),為部分區(qū)域性干旱;當(dāng)10%≤P<25%時(shí),為局域性干旱;當(dāng)P<10%時(shí)可認(rèn)為區(qū)域內(nèi)無明顯干旱發(fā)生。

    1.2.4 干旱演變特征分析

    非參數(shù)Mann-kendall(M-K)趨勢檢驗(yàn)是由世界氣象組織(WMO)推薦并使用的一種統(tǒng)計(jì)檢驗(yàn)方法,其優(yōu)點(diǎn)是不需要樣本遵從一定的分布,也不受少數(shù)異常值的干擾。本研究采用非參數(shù)M-K趨勢檢驗(yàn),探究36個(gè)站點(diǎn)的演變特征,分析干旱演變趨勢。對于時(shí)間序列(1,2,…,X)樣本,定義檢驗(yàn)統(tǒng)計(jì)量:

    式中:XX分別為第、時(shí)間序列所對應(yīng)的觀測值,且<;當(dāng)X-X小于、等于或大于0時(shí),sign(X-X)分別為-1、0或1;為樣本數(shù)。

    Man-Kendall統(tǒng)計(jì)量的計(jì)算式為:

    對于給定的置信水平,若||>1-α/2,表明該時(shí)間序列在置信水平上存在顯著的增加或減小趨勢。其中,為正值表示增加趨勢,為負(fù)值表示減少趨勢,||≥1.28、1.64、2.32時(shí),分別表示通過了置信度90%、95%、99%的顯著性檢驗(yàn)。本研究對36個(gè)氣象站觀測值計(jì)算得到的12、3系列進(jìn)行趨勢檢驗(yàn)。如果滿足||>1-α/2且>0,表明該站點(diǎn)干旱演變呈減小趨勢;滿足||>1-α/2且<0,表明該站點(diǎn)干旱演變呈增加趨勢;||≤1-α/2的站點(diǎn)代表干旱無明顯的增大或減小趨勢。

    1.2.5 數(shù)據(jù)分析

    計(jì)算和非參數(shù)Mann-kendall趨勢檢驗(yàn)基于Matlab 2012軟件編程實(shí)現(xiàn)。空間分析基于ArcGIS10.2軟件的地統(tǒng)計(jì)分析模塊實(shí)現(xiàn)。常用的空間數(shù)據(jù)插值方法有反距離權(quán)重插值法(IDW)、多項(xiàng)式插值法(Spline)和克里金(Kriging)插值法3種,本研究采用交叉驗(yàn)證法評價(jià)這3種方法的插值效果。經(jīng)對比分析,在3種插值方法中,IDW的插值精度相對較高,因此本研究采用IDW進(jìn)行空間數(shù)據(jù)插值分析??臻g數(shù)據(jù)插值分析時(shí),采用自然間斷法進(jìn)行數(shù)據(jù)分類。根據(jù)云南省氣候特點(diǎn),季節(jié)劃分標(biāo)準(zhǔn)為:3、4、5月為春季,6、7、8月為夏季,9、10、11月為秋季,12、1、2月為冬季。依據(jù)云南省綜合農(nóng)業(yè)區(qū)劃報(bào)告,將云南省進(jìn)一步劃分為滇中區(qū)、滇東南區(qū)、滇西區(qū)、滇西南區(qū)、滇西區(qū)、滇東北區(qū)、滇西北區(qū)共7個(gè)分區(qū)。

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

    本研究采用的基礎(chǔ)氣象數(shù)據(jù)為均勻分布于云南省的36個(gè)國家基本氣象站1982—2012年逐月降水觀測資料,來源于中國氣象數(shù)據(jù)網(wǎng)。36個(gè)國家基本降水站空間分布見圖1。

    2 結(jié)果與分析

    2.1 SPI12、SPI3變化過程

    圖2為1982—2012年云南省36個(gè)降水站在年尺度、季節(jié)尺度上的平均變化。云南省近30年不同時(shí)間尺度均呈波動(dòng)變化趨勢,但不同時(shí)段的變化趨勢存在差異。除春季略有增長外,夏季、秋季、冬季和年尺度均表現(xiàn)為小幅度減小趨勢。春季變化范圍為-1.92~2.07,對應(yīng)的干旱等級范圍為嚴(yán)重干旱~無旱;夏季變化范圍為-2.22~1.01,對應(yīng)的干旱等級范圍為極端干旱~無旱;秋季變化范圍為-2.67~1.37,對應(yīng)的干旱等級范圍為極端干旱~無旱;冬季變化范圍為-1.78~2.23,對應(yīng)的干旱等級范圍為嚴(yán)重干旱~無旱??傮w來看,秋季變化幅度最大,冬季變化幅度最小。

    2.2 年尺度干旱頻率分布

    圖3為云南省年尺度干旱頻率空間分布。由圖3(a)可知,干旱頻率的高值區(qū)分布在滇東南區(qū)的文山、紅河,滇中區(qū)的昆明、楚雄,以及滇西區(qū)的德宏等地區(qū),干旱頻率分布在41%~50%之間,這些區(qū)域發(fā)生干旱的頻率總體較高。干旱頻率的低值區(qū)主要分布在滇東北區(qū)的昭通,滇西區(qū)的保山以及滇西南區(qū)的普洱等地區(qū),干旱頻率分布在32%~37%之間。其余地區(qū)的干旱頻率介于二者之間,干旱頻率分布在37%~41%之間。

    由圖3(b)可知,輕微干旱的頻率分布在6.5%~21%之間,干旱頻率平均值為11.82%,發(fā)生干旱頻率總體較小。高值區(qū)位于滇中區(qū)的楚雄、昆明、玉溪地區(qū),以及德宏、保山、普洱、怒江的部分地區(qū),其余地區(qū)發(fā)生輕微干旱的頻率相對較小。由圖3(c)年尺度中等干旱頻率分布可知,高值區(qū)主要分布在怒江、麗江、臨滄的部分地區(qū),其他地區(qū)的頻率均小于10%。圖3(d)為年尺度嚴(yán)重干旱頻率分布,高值區(qū)主要分布在西雙版納、麗江、紅河、曲靖的部分區(qū)域,其他地區(qū)的頻率均小于9.5%。圖3(e)為年尺度極端干旱頻率分布,高值區(qū)主要分布在滇西北區(qū)的迪慶,以及滇東南區(qū)的文山、紅河的部分地區(qū),其他地區(qū)的頻率小于14%??傮w來看,年尺度上4個(gè)干旱強(qiáng)度等級的發(fā)生頻率為:輕微干旱(11.82%)>極端干旱(11.14%)>中等干旱(8.95%)>嚴(yán)重干旱(7.17%),表明云南省發(fā)生輕微干旱和極端干旱的頻率總體較高。

    圖3 年尺度干旱頻率空間分布

    2.3 年尺度連續(xù)干旱特征

    為識別連續(xù)干旱特征,本研究統(tǒng)計(jì)了36個(gè)站點(diǎn)1982—2012年連續(xù)干旱年的出現(xiàn)次數(shù)。統(tǒng)計(jì)結(jié)果表明,在36個(gè)站點(diǎn)中僅有貢山、孟定2個(gè)站點(diǎn)未出現(xiàn)連續(xù)干旱,其余34個(gè)站點(diǎn)均發(fā)生過至少1次的連續(xù)2 a以上的干旱。從發(fā)生頻次來看,出現(xiàn)2 a連續(xù)干旱的頻次最高,合計(jì)有29個(gè)站點(diǎn)共發(fā)生54次連續(xù)2 a干旱,占比高達(dá)57.45%。其次為連續(xù)3 a干旱,合計(jì)有17個(gè)站點(diǎn)共發(fā)生23次連續(xù)3 a干旱,占比高達(dá)24.47%。連續(xù)4 a干旱的發(fā)生次數(shù)為15次,在連續(xù)干旱年中的占比為15.96%。干旱持續(xù)時(shí)間最長為6 a,分別是昭通地區(qū)的2000—2005年干旱和瀘西地區(qū)的2007—2012年連續(xù)干旱??傮w來看,昆明、楚雄、大理等區(qū)域發(fā)生連續(xù)干旱的頻率較高,這些區(qū)域也是云南省干旱易發(fā)多發(fā)區(qū),應(yīng)更加注重連續(xù)干旱的防控。

    2.4 年尺度干旱范圍變化

    圖4為年尺度干旱范圍變化過程曲線。云南省干旱范圍變化呈顯著波動(dòng)趨勢。干旱范圍最大的是1987—1989、1992、2003、2005、2009—2012年。除2009年外,干旱范圍均超過了60%,屬于全域性干旱。1982、1984、1993、1996、1998、2006年干旱范圍分布在33%~50%之間,屬區(qū)域性干旱。1994、2004年的干旱范圍分布在25%~33%之間,屬于部分區(qū)域性干旱。其余年份的干旱范圍均在25%以下,屬于局域性干旱或無明顯干旱。總體來看,干旱范圍越大,干旱延續(xù)時(shí)間也越長。

    不同等級干旱范圍曲線也隨著時(shí)間變化呈波動(dòng)趨勢。在同一時(shí)點(diǎn),干旱范圍曲線下方第一條干旱等級范圍曲線為主導(dǎo)該時(shí)點(diǎn)的干旱等級。從干旱范圍變化來看,大部分區(qū)域性干旱由極端干旱主導(dǎo),干旱過程表現(xiàn)為范圍廣、等級高,干旱導(dǎo)致的災(zāi)害損失也較為嚴(yán)重,這也是云南省大范圍干旱發(fā)生、發(fā)展過程中的顯著特征。

    圖4 年尺度干旱范圍變化曲線

    2.5 季節(jié)性干旱頻率分布

    季節(jié)性干旱與區(qū)域耕地利用、種植制度、城鄉(xiāng)供水保障等活動(dòng)相關(guān)。云南省春、夏、秋、冬4個(gè)季節(jié)對應(yīng)的干旱頻率分別為39.86%、39.83%、38.43%、41.33%,4個(gè)季節(jié)的干旱頻率差異較小。表2為季節(jié)性干旱等級頻率統(tǒng)計(jì)結(jié)果。盡管4個(gè)季節(jié)的干旱頻率差異較小,但不同季節(jié)干旱等級對應(yīng)的頻率差異較大,春季、夏季、冬季發(fā)生輕微干旱的頻率最大,而秋季發(fā)生極端干旱的頻率最大。秋季干旱頻率最高的干旱等級為極端干旱,表明在干旱頻率差異不大的情況下,秋季較高等級的干旱強(qiáng)度對干旱災(zāi)變過程影響也最大。

    圖5為不同季節(jié)干旱頻率分布。春季干旱頻率分布呈自滇西北向東南遞減的趨勢;干旱頻率的高值區(qū)主要分布在怒江、大理及麗江的部分地區(qū),干旱頻率分布在40.6%~50.0%之間,這些地區(qū)也是春旱防范的重點(diǎn)區(qū)域;干旱頻率的低值區(qū)主要分布在滇東南區(qū)的文山、紅河等地,干旱頻率分布在32.7%~39.4%之間。夏季干旱頻率分布相對較均勻,除文山等部分地區(qū)干旱頻率較高外,其余大部分地區(qū)的干旱頻率均在39.5%以下。秋季干旱頻率分布也相對均勻,高值區(qū)主要分布在怒江地區(qū),其余地區(qū)的干旱頻率均在39.3%以下。冬季干旱頻率空間分布變異性較大,高值區(qū)分布在滇中區(qū)的昆明、玉溪,滇東南區(qū)的文山以及南部區(qū)的西雙版納,干旱頻率分布在42.8%~50.0%之間;低值區(qū)分布在紅河、臨滄、昭通等部分地區(qū),干旱頻率分布在30.4%~41.4%之間。

    圖5 不同季節(jié)干旱頻率空間分布

    2.6 季節(jié)性干旱范圍變化

    圖6為季節(jié)性干旱范圍變化曲線。不同季節(jié)干旱范圍變化呈顯著的波動(dòng)變化趨勢。由圖6(a)可知,多數(shù)年份的春季干旱處于區(qū)域干旱以上范圍,干旱范圍最高值發(fā)生在1989年,干旱范圍接近100%;1982、1987、1995、2003、2005、2012年的春季干旱范圍也較大,均超過了60%;不同年份主導(dǎo)春季干旱的強(qiáng)度等級也存在差異,1997年前為輕微干旱和極端干旱,1997年后為輕微干旱和中等干旱。由圖6(b)可知,多數(shù)年份的夏季干旱也處于區(qū)域干旱以上范圍,干旱范圍最高值發(fā)生在1992年,干旱范圍接近100%;1986—1987、1992、2003、2006、2011年的夏季干旱范圍也較大,均超過了60%;除2006年外,干旱范圍超過50%的夏季全域性干旱均由極端干旱主導(dǎo)。從圖6(c)可以看出,多數(shù)年份的秋季干旱處于區(qū)域干旱以上范圍,干旱范圍最高值發(fā)生在2006年,干旱范圍超過80%;1996、2003—2004、2005、2009年的秋季干旱范圍也較大,均超過了60%;1997年前干旱范圍無明顯變化,1997年后干旱范圍有明顯增大趨勢,且主導(dǎo)的干旱等級為極端干旱。從圖6(d)可以看出,冬季干旱范圍波動(dòng)較大,干旱范圍均在20%以上;干旱范圍最高值發(fā)生在1986年,干旱范圍也接近100%;1984—1986、1996、2001、2004、2006、2009、2012年的秋季干旱范圍也較大,均超過了60%。

    上述結(jié)果表明,云南省不同季節(jié)干旱范圍變化存在顯著差異,干旱范圍高值的年份也各不相同,同一年內(nèi)4個(gè)季節(jié)的干旱范圍也存在差異。整體來看,干旱等級和干旱范圍間存在著密切聯(lián)系,總體表現(xiàn)為干旱等級越高,干旱范圍也越大。除春季外,主導(dǎo)各季節(jié)干旱的等級均包含極端干旱,這也是導(dǎo)致云南夏、秋、冬3個(gè)季節(jié)干旱范圍廣、干旱災(zāi)害影響程度大的根本原因。

    圖6 季節(jié)性干旱范圍變化過程曲線

    2.7 干旱的年際演變趨勢檢驗(yàn)

    表3為36個(gè)站點(diǎn)非參數(shù)Mann-Kendall 趨勢檢驗(yàn)結(jié)果。從年尺度來看,不同站點(diǎn)趨勢檢驗(yàn)的值主要分布在||≤1-α/2范圍內(nèi);滿足||>1-α/2,且>0的站點(diǎn)數(shù)為0;滿足||>1-α/2,且<0的站點(diǎn)也較少,在90%、95%、99%的3個(gè)置信區(qū)間的站點(diǎn)數(shù)分別為5、2、2。因此,從年尺度干旱演變趨勢檢驗(yàn)結(jié)果來看,除廣南、景洪、瀘西、沾益、昭通5個(gè)站點(diǎn)在90%的置信區(qū)間內(nèi)有減小趨勢外,云南省大部分地區(qū)值沒有明顯的增長或遞增趨勢。

    從季節(jié)尺度來看,不同季節(jié)的干旱演變趨勢存在差異。春季有18個(gè)站點(diǎn)干旱演變有減小趨勢,其中12個(gè)站點(diǎn)有顯著減小趨勢,3個(gè)站點(diǎn)有極顯著減小趨勢,只有1個(gè)站點(diǎn)干旱演變有增強(qiáng)趨勢。夏季大部分地區(qū)降水充沛,大部分站點(diǎn)干旱演變無明顯的增加或減小趨勢,但仍有保山、德欽、江城、景洪、臨滄、瀘水、瀘西、瑞麗、普洱、宜良、昭通11個(gè)站點(diǎn)表現(xiàn)為增強(qiáng)趨勢,其中江城、景洪、臨滄、瀘西、普洱、昭通6個(gè)站點(diǎn)呈顯著增強(qiáng)趨勢,3個(gè)站點(diǎn)呈極顯著增強(qiáng)趨勢。秋季合計(jì)有3個(gè)站點(diǎn)的干旱演變表現(xiàn)為減小趨勢,有6個(gè)站點(diǎn)表現(xiàn)為增強(qiáng)趨勢,其他站點(diǎn)干旱無明顯的增大或減小趨勢。冬季除沾益站點(diǎn)的干旱演變表現(xiàn)為減小趨勢,景洪、勐臘表現(xiàn)為增強(qiáng)趨勢外,其余大部分站點(diǎn)干旱均無明顯的增大或減小趨勢。不同季節(jié)干旱演變趨勢的差異性特征表明,在干旱防控過程中,可針對不同季節(jié)干旱演變的特點(diǎn)實(shí)施差異化的應(yīng)對措施,以實(shí)現(xiàn)干旱防控的精準(zhǔn)化要求。

    表3 非參數(shù)Mann-Kendall趨勢檢驗(yàn)結(jié)果

    3 討論

    干旱評估方法是近年來國內(nèi)外干旱研究領(lǐng)域的熱點(diǎn)問題。目前,國內(nèi)外較成熟的干旱評估方法主要有降水距平百分比、Palmer指標(biāo)、指數(shù)、指數(shù)、指數(shù)、、等。其中,因具有穩(wěn)定的計(jì)算特性,消除了降水時(shí)空分布差異,同時(shí)能有效地反映旱澇狀況等特點(diǎn),在世界各地得到了廣泛應(yīng)用[17]。通過反演得到的干旱頻率、干旱歷時(shí)、干旱強(qiáng)度等信息與云南省發(fā)生的干旱記錄具有較好的一致性,這也說明在云南省干旱評估中具有優(yōu)秀的適宜性。

    同一地區(qū)月(1)、季度(3)、半年(6)、年(12)時(shí)間尺度表現(xiàn)出不同的變化特性。選擇的時(shí)間尺度過長,干旱評估結(jié)果可能存在滯后效應(yīng);而時(shí)間尺度過短,則使評估結(jié)果存在較大的波動(dòng)性,從而失去對干旱的診斷意義。從國內(nèi)外研究來看,3、12月時(shí)間尺度的應(yīng)用較多,且3和12更能反映干旱年際變化[15]。1和3均能夠較好地反映短期氣象的旱澇特征,與農(nóng)業(yè)干旱間存在較為密切的關(guān)系;6和12則可以更精準(zhǔn)地反映長期的旱澇變化特征[20]。由于時(shí)間尺度不同,同一地區(qū)年尺度、季節(jié)尺度上的干旱演變存在差異。因此,干旱評估中應(yīng)根據(jù)實(shí)際需要合理選擇時(shí)間尺度,以提高干旱評估的有效性。

    云南省近50年降水日數(shù)和雨季降水量逐漸減少,但暴雨、大暴雨等極端降水頻率上升,由此對區(qū)域水資源時(shí)空分布和旱澇狀況造成了顯著影響[21]。也有研究結(jié)果表明,近46年云南省年平均降水量趨于減少,特別是夏季降水量明顯減小[22]。降水變化在區(qū)域干旱的演變中扮演著重要角色,隨著降水特征(時(shí)空分布、暴雨特征等)的變化,干旱特征也隨之發(fā)生變化,最終對農(nóng)業(yè)生產(chǎn)、生態(tài)環(huán)境和城鄉(xiāng)供水安全造成重大影響。任菊章等[23]研究表明,在全球變暖背景下,云南省雨季有變干的趨勢,干旱略有加強(qiáng)趨勢。從本文研究來看,云南省不同地區(qū)的干旱演變趨勢存在差異;全年、夏季、秋季、冬季時(shí)間尺度的干旱演變呈小幅度增強(qiáng)趨勢,這與任菊章等[23]研究結(jié)果基本一致。未來應(yīng)重點(diǎn)關(guān)注氣候變化驅(qū)動(dòng)下干旱演變規(guī)律及驅(qū)動(dòng)過程,為區(qū)域干旱災(zāi)害精準(zhǔn)防控提供科學(xué)依據(jù)。

    4 結(jié)論

    1)云南省在不同時(shí)間尺度上均呈波動(dòng)變化趨勢,除春季有小幅度增長趨勢外,夏季、秋季、冬季和全年尺度的均表現(xiàn)為小幅減小趨勢。

    2)年尺度干旱頻率分布在32.26%~50.0%之間,多年平均干旱頻率為39.08%,輕微干旱和極端干旱為年尺度干旱的主導(dǎo)性等級。干旱頻率高值區(qū)主要分布在文山、紅河、昆明、楚雄、德宏等地區(qū),低值區(qū)主要分布在昭通、曲靖、保山、普洱等地區(qū)。

    3)春、夏、秋、冬4個(gè)季節(jié)對應(yīng)的干旱頻率分別為39.86%、39.83%、38.43%、41.33%,其中春、夏、冬3個(gè)季節(jié)頻率最高的干旱等級為輕微干旱,而秋季發(fā)生頻率最高的干旱等級為極端干旱。

    4)年尺度除有5個(gè)站點(diǎn)的干旱演變表現(xiàn)為增強(qiáng)趨勢外,大部分地區(qū)年尺度干旱無明顯的增強(qiáng)或減小趨勢。

    [1] MA M W, SONG S B, REN L L, et al. Multivariate drought characteristics using trivariate Gaussian and Student t copulas[J]. Hydrological Processes, 2013, 27(8): 1 175-1 190.

    [2] WILHITE Donald A. Drought as a natural hazard: concepts and definitions[C]. London, 2000.

    [3] PALMER W C. Meteorological drought US [C]. New York, 1965.

    [4] World Meteorologic Organization. International meteoroljogical vocabulary (2nd)[M]. WMO, 1966.

    [5] American Meteorological Society. Meteorological drought-Policy statement[J]. Bulletin of American Meteorological Society, 1997, 78: 847-849.

    [6] 韋開, 王全九, 周蓓蓓, 等. 基于降水距平百分率的陜西省干旱時(shí)空分布特征[J]. 水土保持學(xué)報(bào), 2017, 31(1): 318-322.

    WEI Kai, WANG Quanjiu, ZHOU Beibei, et al. Analysis of drought characteristics in Shaanxi Province based on precipitation anomaly percentage[J]. Journal of Soil and Water Conservation, 2017, 31(1): 318-322.

    [7] 劉衛(wèi)國, 王曼, 丁俊祥, 等. 帕默爾干旱指數(shù)在天山北坡典型綠洲干旱特征分析中的適用性[J]. 中國沙漠, 2013, 33(1): 249-257.

    LIU Weiguo, WANG Man, DING Junxiang, et al. A test on the palmer drought index with the data collected in the oasis in the northern slope of Tianshan Mountain[J]. Journal of Desert Research, 2013, 33(1): 249-257.

    [8] 包云軒, 孟翠麗, 申雙和, 等. 基于CI指數(shù)的江蘇省近50年干旱的時(shí)空分布規(guī)律[J]. 地理學(xué)報(bào), 2011, 66(5): 599-608.

    BAO Yunxuan, MENG Cuili, SHEN Shuanghe, et al. Temporal and spatial patterns of droughts for recent 50 years in Jiangsu based on meteorological drought composite index[J]. Acta Geographica Sinica, 2011, 66(5): 599-608.

    [9] 王勁松, 郭江勇, 傾繼祖. 一種K干旱指數(shù)在西北地區(qū)春旱分析中的應(yīng)用[J]. 自然資源學(xué)報(bào), 2007, 22(5): 709-717.

    WANG Jinsong, GUO Jiangyong, QING Jizu. Application of a kind of K drought index in the spring drought analysis in northwest China[J]. Journal of Natural Resources, 2007, 22(5): 709-717.

    [10] 李斌, 解建倉, 胡彥華, 等. 基于標(biāo)準(zhǔn)化降水指數(shù)的陜西省干旱時(shí)空變化特征分析[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2017, 33(17): 113-119.

    LI Bin, XIE Jiancang, HU Yanhua, et al. Analysis on spatiotemporal variability characteristics of drought in Shaanxi Province using standardized precipitation index[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(17): 113-119.

    [11] 陳淑紅, 劉宇, 劉招, 等. 基于SPEI指數(shù)涇河流域近50年干旱時(shí)空演變特征[J]. 水土保持學(xué)報(bào), 2022, 36(2): 181-188, 196.

    CHEN Shuhong, LIU Yu, LIU Zhao, et al. Temporal and spatial characteristics of drought in the Jinghe River Basin in the past 50 years based on the SPEI index[J]. Journal of Soil and Water Conservation, 2022, 36(2): 181-188, 196.

    [12] 陳南祥, 李士輝, 楊莉. 灰色模糊聚類理論在干旱評價(jià)中的應(yīng)用[J].華北水利水電學(xué)院學(xué)報(bào), 2007, 28(2): 4-7.

    CHEN Nanxiang, LI Shihui, YANG Li. Application about the grey-fuzzy cluster theory in the arid evaluation[J]. Journal of North China Institute of Water Conservancy and Hydroelectric Power, 2007, 28(2): 4-7.

    [13] 武翡翡, 王烜. 基于熵權(quán)的模糊物元模型在白洋淀干旱評價(jià)中的應(yīng)用[C]//Proceedings of 2011 AASRI Conference on Applied Information Technology (AASRI-AIT 2011 V3). Kota Kinabalu, Malaysia, 2011: 325-328.

    [14] MCKEE T B, DOESKEN N J, KLEIST J. The relationship of drought frequency and duration to time scales. Proceeding of 8th Conference on Applied Climatology[C]. American Meteorological Society, Boston: Massa chusetts, 1993, 179-184.

    [15] MERABTI A, MARTINS D S, MEDDI M, et al. Spatial and time variability of drought based on SPI and RDI with various time scales[J]. Water Resources Management, 2018, 32(3): 1 087-1 100.

    [16] XU K, YANG D W, YANG H B, et al. Spatio-temporal variation of drought in China during 1961–2012: A climatic perspective[J]. Journal of Hydrology, 2015, 526: 253-264.

    [17] 袁文平, 周廣勝. 標(biāo)準(zhǔn)化降水指標(biāo)與Z指數(shù)在我國應(yīng)用的對比分析[J].植物生態(tài)學(xué)報(bào), 2004, 28(4): 523-529.

    YUAN Wenping, ZHOU Guangsheng. Comparison between standardized precipitation index and z-index in China[J]. Acta Phytoecologica Sinica, 2004, 28(4): 523-529.

    [18] 程建剛, 解明恩. 近50年云南區(qū)域氣候變化特征分析[J]. 地理科學(xué)進(jìn)展, 2008, 27(5): 19-26.

    CHENG Jiangang, XIE Ming’en. The analysis of regional climate change features over Yunnan in recent 50 years[J]. Progress in Geography, 2008, 27(5): 19-26.

    [19] KIM K, WANG M C, RANJITKAR S, et al. Using leaf area index (LAI) to assess vegetation response to drought in Yunnan Province of China[J]. Journal of Mountain Science, 2017, 14(9): 1 863-1 872.

    [20] 張利利, 周俊菊, 張恒瑋, 等. 基于的石羊河流域氣候干濕變化及干旱事件的時(shí)空格局特征研究[J]. 生態(tài)學(xué)報(bào), 2017, 37(3): 996-1 007.

    ZHANG Lili, ZHOU Junju, ZHANG Hengwei, et al. Temporal and spatial patterns of climate drought-wet and drought event based on standard precipitation index in Shiyang River Basin[J]. Acta Ecologica Sinica, 2017,37(3): 996-1 007.

    [21] 吳志杰. 城市化對云南高原中部區(qū)域氣候變化的影響研究[D]. 昆明: 云南大學(xué), 2015.

    WU Zhijie. Quantitative estimates of urbanization effects on climate change in central Yunnan Plateau[D]. Kunming: Yunnan University, 2015.

    [22] LI Y G, HE D M, HU J M, et al. Variability of extreme precipitation over Yunnan Province, China 1960—2012[J]. International Journal of Climatology, 2015, 35(2): 245-258.

    [23] 任菊章, 黃中艷, 鄭建萌. 基于相對濕潤度指數(shù)的云南干旱氣候變化特征[J]. 中國農(nóng)業(yè)氣象, 2014, 35(5): 567-574.

    REN Juzhang, HUANG Zhongyan, ZHENG Jianmeng. Analysis on drought climate change in Yunnan based on relative moisture index[J]. Chinese Journal of Agrometeorology, 2014, 35(5): 567-574.

    Spatiotemporal Variation of Drought in Yunnan Province Calculated Using the Standardized Precipitation Index

    CHEN Zhengfa1,2, LI Jing1,2, XIANG Biao1*, DUAN Qingsong1,2, LI Shufang1, WANG Shufang1,2

    (1. College of Water Conservancy, Yunnan Agricultural University, Kunming 650201, China;2. Research Center of Agricultural Water Saving Engineering and Technology in Yunnan Province, Kunming 650201, China)

    【Objective】Drought is the most important abiotic stress facing agricultural management and economic development in many regions and countries. The aim of this paper is to identify the spatiotemporal variation of drought in Yunnan Province, using the standardized precipitation index ().【Method】Meteorological data measured over the past 31 years from 36 meteorological stations across the province were used in the analysis. Multi-scale variations of the droughts in both time and space were calculated using, GIS, nonparametric Mann-Kendall trend test and other methods.【Result】①Both annual and seasonalfluctuated temporally, but we identified a small annual increase in spring, small annual decreases in summer, autumn and winter. The autumnvaried mostly, and the winterleast. The annualvaried from -1.95 to 1.44. ②The annual drought frequency varied from 32.26% to 50.00%, with the severity of most droughts varying from mild to extreme range. Spatially, drought occurred more frequently in Wenshan, Honghe, Kunming, Chuxiong and Dehong than in other regions in the province. ③Seasonally, drought frequency in spring, summer, autumn and winter was 39.86%, 39.83%, 38.43%, and 41.33%, respectively. Although drought frequency in different seasons was close, their severity differed considerably. ④Annually, except for five meteorological stations which found increased drought frequency, drought frequency in other stations remained statistically the same. Seasonally, most stations showed decreased drought frequency in spring, eleven stations showed increased drought frequency in summer, and six and two stations showed increased autumn droughts and winter droughts, respectively.【Conclusion】Annual and seasonal drought frequency is spatially heterogeneous across the province. Droughts in the province are seasonal, with droughts occurring in spring and winter more frequently and severely.

    space-time evolution; meteorological drought; seasonal drought; standardized precipitation index; Yunnan Province

    1672 - 3317(2023)04 - 0092 - 08

    P237

    A

    10.13522/j.cnki.ggps.2022250

    陳正發(fā), 李靖, 相彪, 等. 基于的云南省多尺度干旱時(shí)空演變特征識別[J]. 灌溉排水學(xué)報(bào), 2023, 42(4): 92-99.

    CHEN Zhengfa, LI Jing, XIANG Biao, et al. Spatiotemporal Variation of Drought in Yunnan Province Calculated Using the Standardized Precipitation Index[J]. Journal of Irrigation and Drainage, 2023, 42(4): 92-99.

    2022-05-06

    云南省基礎(chǔ)研究專項(xiàng)面上項(xiàng)目(202201AT070272);云南省教育廳科研項(xiàng)目(2022J0303);云南農(nóng)業(yè)大學(xué)引進(jìn)人才科研啟動(dòng)項(xiàng)目(KY2022-29)

    陳正發(fā)(1985-),男。講師,博士,主要從事水土生態(tài)工程研究。E-mail: chenzhengfa2013@126.com

    相彪(1979-),男。正高級工程師,博士,主要從事農(nóng)業(yè)水工建筑研究。E-mail: 103931114@qq.com

    責(zé)任編輯:韓 洋

    猜你喜歡
    時(shí)空云南省站點(diǎn)
    云南省安寧市老年大學(xué)之歌
    跨越時(shí)空的相遇
    民國初期云南省議會述論
    鏡中的時(shí)空穿梭
    基于Web站點(diǎn)的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    2017~2018年冬季西北地區(qū)某站點(diǎn)流感流行特征分析
    玩一次時(shí)空大“穿越”
    云南省瑞麗市老年大學(xué)之歌
    首屆歐洲自行車共享站點(diǎn)協(xié)商會召開
    中國自行車(2017年1期)2017-04-16 02:53:52
    怕被人認(rèn)出
    故事會(2016年21期)2016-11-10 21:15:15
    神马国产精品三级电影在线观看 | 亚洲中文日韩欧美视频| 亚洲欧美精品综合久久99| 2021天堂中文幕一二区在线观 | 在线观看午夜福利视频| 午夜福利18| 成人亚洲精品av一区二区| 国内毛片毛片毛片毛片毛片| 最好的美女福利视频网| 国产伦一二天堂av在线观看| 精品无人区乱码1区二区| 丝袜人妻中文字幕| 99在线视频只有这里精品首页| 黄片大片在线免费观看| 久久亚洲真实| 给我免费播放毛片高清在线观看| 欧美中文日本在线观看视频| 国产精品久久电影中文字幕| 精品一区二区三区av网在线观看| 国产亚洲精品av在线| 9191精品国产免费久久| 女同久久另类99精品国产91| 国产野战对白在线观看| 国产爱豆传媒在线观看 | 免费在线观看影片大全网站| 国产av在哪里看| 午夜a级毛片| 日韩欧美一区二区三区在线观看| 亚洲熟妇中文字幕五十中出| 欧美成人性av电影在线观看| 精品欧美国产一区二区三| 精品一区二区三区四区五区乱码| 一卡2卡三卡四卡精品乱码亚洲| 成人亚洲精品av一区二区| 91大片在线观看| 怎么达到女性高潮| 久久热在线av| 中文资源天堂在线| 黑人巨大精品欧美一区二区mp4| 999精品在线视频| 日韩免费av在线播放| 免费看a级黄色片| 亚洲av成人一区二区三| 91大片在线观看| 午夜日韩欧美国产| 又紧又爽又黄一区二区| 精品人妻1区二区| 老熟妇仑乱视频hdxx| 最新美女视频免费是黄的| 好男人在线观看高清免费视频 | 天天躁夜夜躁狠狠躁躁| 99精品在免费线老司机午夜| 在线观看午夜福利视频| 国产成年人精品一区二区| 在线观看www视频免费| 可以在线观看毛片的网站| x7x7x7水蜜桃| 男女床上黄色一级片免费看| 韩国精品一区二区三区| 国产爱豆传媒在线观看 | 12—13女人毛片做爰片一| 亚洲 国产 在线| 国产又色又爽无遮挡免费看| 一二三四社区在线视频社区8| av在线天堂中文字幕| 男人舔奶头视频| 欧美大码av| 日韩大码丰满熟妇| 视频在线观看一区二区三区| 久久九九热精品免费| 亚洲成a人片在线一区二区| 日韩一卡2卡3卡4卡2021年| 特大巨黑吊av在线直播 | 在线观看舔阴道视频| 国产亚洲精品综合一区在线观看 | 老司机靠b影院| 国产精品久久久av美女十八| 97碰自拍视频| 亚洲国产看品久久| 国产在线观看jvid| 国产精品免费一区二区三区在线| 国产亚洲欧美在线一区二区| 亚洲精华国产精华精| av福利片在线| 婷婷丁香在线五月| 亚洲av电影不卡..在线观看| 国产激情偷乱视频一区二区| 日本成人三级电影网站| 午夜福利欧美成人| 十八禁网站免费在线| 少妇裸体淫交视频免费看高清 | 色在线成人网| 国产一区二区在线av高清观看| 大型av网站在线播放| 制服诱惑二区| 欧美激情高清一区二区三区| 亚洲男人的天堂狠狠| 麻豆成人午夜福利视频| 91av网站免费观看| 天堂√8在线中文| 人人妻,人人澡人人爽秒播| 色综合亚洲欧美另类图片| 久久久久久大精品| 免费看十八禁软件| 成人亚洲精品av一区二区| 中文字幕高清在线视频| 国产高清视频在线播放一区| 久久精品国产亚洲av香蕉五月| 色综合站精品国产| 久久性视频一级片| 欧美成人免费av一区二区三区| 黄片播放在线免费| 香蕉久久夜色| 国产亚洲av高清不卡| 一区二区三区高清视频在线| 亚洲,欧美精品.| 在线天堂中文资源库| 国产成人精品久久二区二区91| 99久久无色码亚洲精品果冻| 国产成人av教育| 国产高清videossex| 脱女人内裤的视频| 亚洲一区二区三区色噜噜| 亚洲第一av免费看| 亚洲人成77777在线视频| 午夜免费激情av| www国产在线视频色| 国产精品野战在线观看| 一区二区三区激情视频| 国产精品1区2区在线观看.| 成人永久免费在线观看视频| 亚洲欧美精品综合一区二区三区| 国产三级在线视频| 亚洲无线在线观看| 欧美成人一区二区免费高清观看 | www.自偷自拍.com| 90打野战视频偷拍视频| 一区二区三区高清视频在线| 侵犯人妻中文字幕一二三四区| 天天一区二区日本电影三级| 日日爽夜夜爽网站| 国产激情久久老熟女| tocl精华| 在线看三级毛片| 啦啦啦观看免费观看视频高清| 777久久人妻少妇嫩草av网站| 91成人精品电影| 成人亚洲精品一区在线观看| 麻豆成人午夜福利视频| 一级毛片女人18水好多| 国产1区2区3区精品| 精品少妇一区二区三区视频日本电影| 亚洲第一av免费看| 搡老妇女老女人老熟妇| 午夜亚洲福利在线播放| 欧美色视频一区免费| 亚洲国产精品久久男人天堂| 熟妇人妻久久中文字幕3abv| 国产高清视频在线播放一区| 欧美激情 高清一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 91国产中文字幕| www.精华液| 搞女人的毛片| 国产爱豆传媒在线观看 | 久久国产精品男人的天堂亚洲| 久久香蕉国产精品| 亚洲成av人片免费观看| 一级作爱视频免费观看| 丝袜人妻中文字幕| tocl精华| 国产区一区二久久| 欧美一级毛片孕妇| 精品久久蜜臀av无| 很黄的视频免费| 免费在线观看成人毛片| 国产精品亚洲av一区麻豆| 精品久久蜜臀av无| 国产激情偷乱视频一区二区| 国产精品综合久久久久久久免费| 日本 av在线| 国产在线观看jvid| 国内毛片毛片毛片毛片毛片| 欧美成人一区二区免费高清观看 | 国产精品一区二区三区四区久久 | 免费在线观看亚洲国产| 看片在线看免费视频| 国产一区二区三区视频了| 18禁国产床啪视频网站| 啦啦啦 在线观看视频| 亚洲五月婷婷丁香| 中文字幕久久专区| 精品国产超薄肉色丝袜足j| 一区福利在线观看| 一个人观看的视频www高清免费观看 | 老司机靠b影院| 2021天堂中文幕一二区在线观 | 九色国产91popny在线| 天天一区二区日本电影三级| 久久狼人影院| 香蕉丝袜av| 国产成人av教育| 日日爽夜夜爽网站| ponron亚洲| 国产精品一区二区三区四区久久 | 午夜福利视频1000在线观看| 午夜免费鲁丝| 此物有八面人人有两片| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产中文字幕在线视频| 热re99久久国产66热| 美女免费视频网站| 免费人成视频x8x8入口观看| 亚洲人成网站高清观看| 午夜福利18| 国产一区二区激情短视频| 给我免费播放毛片高清在线观看| 欧美乱码精品一区二区三区| 香蕉av资源在线| 97超级碰碰碰精品色视频在线观看| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 亚洲国产精品成人综合色| 日本撒尿小便嘘嘘汇集6| 欧美不卡视频在线免费观看 | 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 久久亚洲真实| 黑人欧美特级aaaaaa片| 久9热在线精品视频| 亚洲人成网站高清观看| 巨乳人妻的诱惑在线观看| 一夜夜www| xxx96com| 亚洲国产欧美日韩在线播放| 久久99热这里只有精品18| 热99re8久久精品国产| 亚洲中文日韩欧美视频| 天堂动漫精品| 亚洲国产欧洲综合997久久, | 777久久人妻少妇嫩草av网站| 国产av一区二区精品久久| 一个人观看的视频www高清免费观看 | 亚洲欧美精品综合久久99| 国内久久婷婷六月综合欲色啪| 99国产精品一区二区三区| 亚洲精品在线观看二区| 91成年电影在线观看| 99在线视频只有这里精品首页| 中文字幕精品免费在线观看视频| 久久精品aⅴ一区二区三区四区| 长腿黑丝高跟| 女人被狂操c到高潮| 国产精品乱码一区二三区的特点| 成人三级黄色视频| 一级作爱视频免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产亚洲av麻豆专区| 999精品在线视频| 久久亚洲精品不卡| 母亲3免费完整高清在线观看| 99在线视频只有这里精品首页| www国产在线视频色| 免费在线观看视频国产中文字幕亚洲| 在线永久观看黄色视频| 欧美日韩一级在线毛片| 在线观看舔阴道视频| 丁香欧美五月| 99在线视频只有这里精品首页| 丁香六月欧美| 欧美在线黄色| 国产麻豆成人av免费视频| 热re99久久国产66热| 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 中亚洲国语对白在线视频| 中文字幕久久专区| 久久国产乱子伦精品免费另类| 午夜亚洲福利在线播放| av片东京热男人的天堂| 亚洲电影在线观看av| 宅男免费午夜| 免费在线观看黄色视频的| 91九色精品人成在线观看| 婷婷六月久久综合丁香| 中文字幕人妻熟女乱码| 国产亚洲精品综合一区在线观看 | 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 免费看日本二区| 欧美日韩瑟瑟在线播放| 亚洲成av片中文字幕在线观看| 99精品久久久久人妻精品| 亚洲中文日韩欧美视频| 午夜福利欧美成人| 欧美色视频一区免费| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观 | 亚洲精品av麻豆狂野| 好男人电影高清在线观看| 久久中文字幕人妻熟女| 亚洲国产中文字幕在线视频| 久久婷婷成人综合色麻豆| a级毛片a级免费在线| 别揉我奶头~嗯~啊~动态视频| 欧美成人性av电影在线观看| 亚洲成人久久性| 欧美三级亚洲精品| 男人舔奶头视频| 777久久人妻少妇嫩草av网站| 在线视频色国产色| 久久精品国产亚洲av香蕉五月| 黄片小视频在线播放| 亚洲精品美女久久av网站| 91av网站免费观看| 狠狠狠狠99中文字幕| 中亚洲国语对白在线视频| 午夜精品在线福利| 首页视频小说图片口味搜索| 女人被狂操c到高潮| 最新美女视频免费是黄的| 免费无遮挡裸体视频| 波多野结衣高清无吗| 欧美日韩中文字幕国产精品一区二区三区| 无人区码免费观看不卡| 性欧美人与动物交配| 精品国产国语对白av| 久久国产精品男人的天堂亚洲| АⅤ资源中文在线天堂| 国产区一区二久久| 色综合欧美亚洲国产小说| 成人三级做爰电影| 热re99久久国产66热| 女人爽到高潮嗷嗷叫在线视频| 成人av一区二区三区在线看| av在线天堂中文字幕| 精品日产1卡2卡| 中文资源天堂在线| www日本黄色视频网| 中国美女看黄片| 啦啦啦观看免费观看视频高清| 久久青草综合色| 成人欧美大片| 亚洲国产毛片av蜜桃av| 一区二区三区精品91| 啦啦啦观看免费观看视频高清| 免费看a级黄色片| 一夜夜www| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 最新美女视频免费是黄的| 欧美+亚洲+日韩+国产| 美女高潮到喷水免费观看| 国产精品亚洲av一区麻豆| 午夜日韩欧美国产| 欧美一区二区精品小视频在线| 亚洲 欧美一区二区三区| 麻豆成人av在线观看| 久久婷婷成人综合色麻豆| 18禁国产床啪视频网站| 日本 欧美在线| 亚洲欧美激情综合另类| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 一个人免费在线观看的高清视频| 久久天堂一区二区三区四区| 国内精品久久久久久久电影| 久久久久久久久中文| 草草在线视频免费看| 一级作爱视频免费观看| 国产精品98久久久久久宅男小说| 亚洲精品美女久久av网站| 久久久久久亚洲精品国产蜜桃av| 最近最新中文字幕大全电影3 | 亚洲av成人不卡在线观看播放网| 亚洲人成77777在线视频| 成人国语在线视频| 成人精品一区二区免费| 丝袜在线中文字幕| 精华霜和精华液先用哪个| 精品不卡国产一区二区三区| 1024香蕉在线观看| 久久精品国产亚洲av香蕉五月| 欧美黑人巨大hd| 久久久久精品国产欧美久久久| 欧美激情久久久久久爽电影| 国产三级在线视频| 少妇被粗大的猛进出69影院| 久久香蕉国产精品| 制服诱惑二区| 色av中文字幕| 欧美性长视频在线观看| 亚洲熟女毛片儿| 午夜福利欧美成人| 亚洲国产精品成人综合色| 我的亚洲天堂| bbb黄色大片| www国产在线视频色| 国产午夜福利久久久久久| 日本a在线网址| 国产精品野战在线观看| 国产亚洲精品第一综合不卡| 成人手机av| 国产亚洲精品一区二区www| 在线观看www视频免费| 淫秽高清视频在线观看| 久久午夜综合久久蜜桃| 亚洲专区中文字幕在线| 男女视频在线观看网站免费 | 欧美av亚洲av综合av国产av| 一级作爱视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产三级黄色录像| 夜夜夜夜夜久久久久| 女同久久另类99精品国产91| 色综合站精品国产| 可以在线观看的亚洲视频| 色婷婷久久久亚洲欧美| 国产亚洲欧美在线一区二区| 亚洲人成网站高清观看| 国产麻豆成人av免费视频| 高潮久久久久久久久久久不卡| 男女午夜视频在线观看| 久久久国产成人免费| 亚洲欧美精品综合一区二区三区| 国产一区二区激情短视频| 亚洲国产中文字幕在线视频| 少妇裸体淫交视频免费看高清 | 巨乳人妻的诱惑在线观看| 麻豆国产av国片精品| 欧美日韩精品网址| 国产精品久久久人人做人人爽| 香蕉国产在线看| 国产成人精品久久二区二区免费| 国产精品一区二区三区四区久久 | 婷婷六月久久综合丁香| 在线观看日韩欧美| 亚洲色图av天堂| 亚洲无线在线观看| 亚洲人成电影免费在线| 脱女人内裤的视频| 一a级毛片在线观看| 国产极品粉嫩免费观看在线| 精品国产超薄肉色丝袜足j| 老司机福利观看| 欧美 亚洲 国产 日韩一| 亚洲精品一区av在线观看| 一本一本综合久久| 精品少妇一区二区三区视频日本电影| 亚洲七黄色美女视频| 久久久久久久久中文| 精品乱码久久久久久99久播| 国产亚洲欧美精品永久| 我的亚洲天堂| 久热这里只有精品99| bbb黄色大片| 国产伦在线观看视频一区| 午夜影院日韩av| 18禁美女被吸乳视频| 首页视频小说图片口味搜索| 久久久久久九九精品二区国产 | 亚洲免费av在线视频| 手机成人av网站| 免费看十八禁软件| 精品高清国产在线一区| 欧美性猛交黑人性爽| 国产久久久一区二区三区| 亚洲,欧美精品.| 无人区码免费观看不卡| 久久性视频一级片| 成人免费观看视频高清| 国产熟女xx| 黄色女人牲交| 无人区码免费观看不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 777久久人妻少妇嫩草av网站| 国产午夜精品久久久久久| 一本精品99久久精品77| av电影中文网址| 亚洲人成网站高清观看| 岛国视频午夜一区免费看| www.熟女人妻精品国产| 免费高清视频大片| 满18在线观看网站| 久久精品亚洲精品国产色婷小说| 国产一区在线观看成人免费| 制服诱惑二区| 日本成人三级电影网站| 两个人视频免费观看高清| 国产亚洲精品av在线| 午夜免费观看网址| 国产成人av教育| 久久久久久九九精品二区国产 | www日本黄色视频网| 无人区码免费观看不卡| 亚洲三区欧美一区| 999久久久国产精品视频| 亚洲真实伦在线观看| 国内毛片毛片毛片毛片毛片| 亚洲精品av麻豆狂野| h日本视频在线播放| 免费大片18禁| 亚洲中文日韩欧美视频| 乱人视频在线观看| 国产在线男女| 少妇熟女欧美另类| 精品少妇黑人巨大在线播放 | 在线国产一区二区在线| 亚洲精品国产成人久久av| 乱码一卡2卡4卡精品| 草草在线视频免费看| 国产精品久久久久久av不卡| 免费av观看视频| 久久精品综合一区二区三区| 国产精品亚洲美女久久久| 俄罗斯特黄特色一大片| 美女xxoo啪啪120秒动态图| 午夜视频国产福利| 一个人观看的视频www高清免费观看| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| ponron亚洲| 中出人妻视频一区二区| 最新中文字幕久久久久| 青春草视频在线免费观看| 最新中文字幕久久久久| 成年女人毛片免费观看观看9| 最好的美女福利视频网| 国产精品久久电影中文字幕| 亚洲美女视频黄频| 免费人成在线观看视频色| 床上黄色一级片| 国产探花极品一区二区| 国产高清有码在线观看视频| 麻豆久久精品国产亚洲av| 国产淫片久久久久久久久| 深夜a级毛片| 国产v大片淫在线免费观看| 国产高清视频在线观看网站| 啦啦啦观看免费观看视频高清| 最近中文字幕高清免费大全6| 国产成人freesex在线 | 久久亚洲国产成人精品v| 色综合站精品国产| 欧美激情久久久久久爽电影| 国产淫片久久久久久久久| 91久久精品电影网| 午夜福利成人在线免费观看| 久久这里只有精品中国| 天堂av国产一区二区熟女人妻| 成人特级av手机在线观看| 亚洲色图av天堂| 久久午夜福利片| 久久精品91蜜桃| 国产人妻一区二区三区在| 成人国产麻豆网| 欧美高清性xxxxhd video| 亚洲熟妇熟女久久| 成人欧美大片| 国产精品久久视频播放| 日本五十路高清| 麻豆成人午夜福利视频| 日本 av在线| 国产精品野战在线观看| 亚洲婷婷狠狠爱综合网| 久久久久久久久久成人| 亚洲最大成人手机在线| eeuss影院久久| 桃色一区二区三区在线观看| 高清日韩中文字幕在线| 国产免费男女视频| 亚洲av第一区精品v没综合| 成人av一区二区三区在线看| 高清毛片免费看| 久久精品夜夜夜夜夜久久蜜豆| 波多野结衣巨乳人妻| 一级毛片我不卡| 在线观看免费视频日本深夜| 性欧美人与动物交配| 亚州av有码| 99精品在免费线老司机午夜| 亚洲精品影视一区二区三区av| 国产探花极品一区二区| 欧美日本亚洲视频在线播放| 久久人人爽人人片av| 在线免费观看的www视频| 亚洲婷婷狠狠爱综合网| 亚洲乱码一区二区免费版| 真实男女啪啪啪动态图| 深夜精品福利| 老熟妇乱子伦视频在线观看| а√天堂www在线а√下载| 1000部很黄的大片| 欧美另类亚洲清纯唯美| 一个人看的www免费观看视频| АⅤ资源中文在线天堂| 国产精品一区二区免费欧美| 免费观看在线日韩| 国产乱人视频| 欧美xxxx黑人xx丫x性爽| 国产大屁股一区二区在线视频| 国产一区二区在线观看日韩| 看黄色毛片网站| 久久国内精品自在自线图片| 国产av麻豆久久久久久久| 午夜老司机福利剧场| 日本色播在线视频| 国产av麻豆久久久久久久| 久久综合国产亚洲精品| 一进一出抽搐gif免费好疼| 尤物成人国产欧美一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| 大型黄色视频在线免费观看|