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

    基于SPEI的北江飛來峽流域干旱時(shí)空演變特征分析*

    2015-05-08 02:16:22冼卓雁武傳號(hào)黃國(guó)如
    災(zāi)害學(xué) 2015年3期
    關(guān)鍵詞:北江時(shí)間尺度尺度

    冼卓雁,武傳號(hào),黃國(guó)如, 2

    (1.華南理工大學(xué) 土木與交通學(xué)院,廣東 廣州 510640;2.華南理工大學(xué) 亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣東 廣州,510640)

    ?

    基于SPEI的北江飛來峽流域干旱時(shí)空演變特征分析*

    冼卓雁1,武傳號(hào)1,黃國(guó)如1, 2

    (1.華南理工大學(xué) 土木與交通學(xué)院,廣東 廣州 510640;2.華南理工大學(xué) 亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣東 廣州,510640)

    氣候變暖導(dǎo)致干旱事件頻發(fā),利用1969-2011年北江飛來峽流域24個(gè)雨量站和4個(gè)氣象站月降雨及氣溫?cái)?shù)據(jù),計(jì)算標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI),采用Mann-Kendall趨勢(shì)分析等方法分析了該地區(qū)多時(shí)間尺度的干旱時(shí)空變化特征。研究結(jié)果表明,SPEI能較好地反映出流域內(nèi)干濕變化特征。在過去的43年中,流域內(nèi)干旱事件頻次總體上呈上升趨勢(shì),且極端干旱事件有所增加,從空間分布看,流域上游和下游出口斷面附近呈現(xiàn)出較為顯著的干旱趨勢(shì),隨著時(shí)間尺度增加,干旱趨勢(shì)愈發(fā)顯著,趨勢(shì)穩(wěn)定性也有所增強(qiáng)。

    干旱;標(biāo)準(zhǔn)化降水蒸散指數(shù);飛來峽流域;時(shí)空特征;趨勢(shì)穩(wěn)定性

    干旱是影響農(nóng)作物生長(zhǎng)、阻礙社會(huì)經(jīng)濟(jì)發(fā)展最嚴(yán)重的自然災(zāi)害之一,水文氣象變化、社會(huì)經(jīng)濟(jì)因素差異以及不同地區(qū)對(duì)水的需求不同,使得干旱難以精準(zhǔn)評(píng)定。干旱指數(shù)是評(píng)估干旱影響并區(qū)別不同干旱期的重要參數(shù),反映內(nèi)容包括干旱強(qiáng)度、持續(xù)時(shí)間、嚴(yán)重程度及空間范圍[1]。目前,國(guó)內(nèi)外干旱研究領(lǐng)域使用最為廣泛的為帕默爾干旱指數(shù)PDSI和標(biāo)準(zhǔn)化降水指數(shù)SPI。PDSI由美國(guó)人Palmer[2]于1965年首次提出,該指數(shù)基于土壤水分平衡原理,綜合考慮前期降水和水分供需,半個(gè)多世紀(jì)以來得到了廣泛驗(yàn)證和進(jìn)一步發(fā)展,由于PDSI計(jì)算過程復(fù)雜且需要資料較多,參數(shù)的空間區(qū)域性較強(qiáng),難以大范圍推廣使用[3]。SPI由Mckee等[4]于1993年提出,該指數(shù)緊扣降水量這個(gè)干旱主要影響因子,不涉及具體的干旱機(jī)理,其時(shí)空適應(yīng)性較強(qiáng),多時(shí)間尺度特性可反映不同方面的水資源狀況。但SPI忽略了蒸發(fā)對(duì)干旱的影響,在全球變暖大背景下,SPI局限性愈發(fā)不可忽略。

    Vicente-Serrano等[5]以SPI為基礎(chǔ),于2010年提出了標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI,以降水與蒸發(fā)的差值作為輸入,采用 log-logistic概率分布函數(shù)來描述其變化,再經(jīng)正態(tài)標(biāo)準(zhǔn)化得到SPEI值。SPEI考慮了地表蒸散發(fā)作用,計(jì)算簡(jiǎn)單,具有多時(shí)間尺度分析特性,可用于干旱監(jiān)測(cè)和識(shí)別。Vicente-Serrano等[6]基于SPEI建立了1901-2006年全球干旱趨勢(shì)的數(shù)據(jù)集,并分析了與SPDI監(jiān)測(cè)結(jié)果的相關(guān)關(guān)系。López-Moreno等[7]研究西班牙埃布羅河流域不同時(shí)間尺度的SPEI的響應(yīng)特征,討論了不同蓄水區(qū)與時(shí)間尺度的相關(guān)性。Potop等[8]計(jì)算1901-2010年間捷克共和國(guó)低洼地區(qū)的SPEI和SPI,發(fā)現(xiàn)兩個(gè)指數(shù)對(duì)不同持續(xù)時(shí)間的干旱事件具有不同的敏感度。莊少偉等[9]選取中國(guó)不同等級(jí)降水區(qū)的代表站點(diǎn),將累積水分虧缺量的經(jīng)驗(yàn)概率分布與log-logistic概率分布進(jìn)行對(duì)比,分析SPEI在中國(guó)的適用性。李偉光等[10-11]以SPEI為基礎(chǔ),對(duì)全國(guó)范圍和華南地區(qū)的干旱長(zhǎng)期趨勢(shì)和突變情況進(jìn)行分析。

    飛來峽流域地處北江中上游,北江流域與西江匯合后流經(jīng)佛山、廣州、江門、珠海等城市,為珠三角的經(jīng)濟(jì)發(fā)展和人民生產(chǎn)生活提供了重要保障,流域內(nèi)干旱將會(huì)干擾樞紐的發(fā)電、航運(yùn)、生態(tài)等功能的正常發(fā)揮,勢(shì)必對(duì)下游珠江三角洲的經(jīng)濟(jì)社會(huì)發(fā)展造成一定的負(fù)面影響。在已有的研究成果中,大多是針對(duì)整個(gè)廣東省或華南地區(qū)進(jìn)行研究,或是以傳統(tǒng)干旱指數(shù)為基礎(chǔ)[12-14]。本文采用更加合理的SPEI作為干旱等級(jí)劃分指標(biāo),綜合描述北江飛來峽流域的干旱時(shí)空演變特征,為飛來峽流域的綜合干旱監(jiān)測(cè)、減少旱澇災(zāi)害損失提供參考依據(jù)。

    1 研究區(qū)域及資料情況

    飛來峽流域位于我國(guó)南方,是北江流域重要組成部分,北江是珠江流域第二大水系,是廣東省最重要的河流之一。飛來峽流域的控制集雨面積為34 097 km2,占北江流域面積的73%。流域?qū)儆趤啛釒Ъ撅L(fēng)氣候,夏季以東南風(fēng)和偏南風(fēng)為主,冬季則為北風(fēng)和偏北風(fēng),年均降水量超過1 700 mm,年內(nèi)降水分布不均,主要集中在4-9月。多年平均氣溫約為21 ℃,年內(nèi)氣溫在10~30 ℃之間變化,夏季最高溫超過30 ℃。

    流域出口斷面為飛來峽水利樞紐,該樞紐坐落在北江干流清遠(yuǎn)市飛來峽鎮(zhèn),以防洪為主,兼顧航運(yùn)和發(fā)電。飛來峽流域?qū)Ρ苯蛷V東省的社會(huì)、人文、環(huán)境意義重大,一旦發(fā)生持續(xù)惡劣的干旱事件,除了對(duì)流域內(nèi)生產(chǎn)生活和生態(tài)環(huán)境具有不良影響外,對(duì)飛來峽水利樞紐的正常運(yùn)作和下游珠江三角洲地區(qū)的社會(huì)經(jīng)濟(jì)發(fā)展也會(huì)造成一定影響。

    本文以飛來峽流域?yàn)檠芯繉?duì)象,研究資料主要包括流域內(nèi)24個(gè)雨量站點(diǎn)的月降雨以及佛岡、連縣、南雄和韶關(guān)等4個(gè)氣象站點(diǎn)的月氣溫?cái)?shù)據(jù),時(shí)間序列長(zhǎng)度均為43年(1969-2011年),站點(diǎn)的位置分布如圖1所示。

    圖1 流域內(nèi)雨量站點(diǎn)和氣象站點(diǎn)分布

    2 研究方法

    2.1 標(biāo)準(zhǔn)化降水蒸散指數(shù)

    根據(jù)研究區(qū)域43年數(shù)據(jù)資料計(jì)算不同時(shí)間尺度的SPEI,尺度為1、3、6、12個(gè)月的計(jì)算結(jié)果分別用SPEI1、SPEI3、SPEI6、SPEI12表示,具體計(jì)算過程如下[5]。

    (1)使用改進(jìn)的Thornthwaite方法計(jì)算月潛在蒸發(fā)量PE。

    (1)

    式中:Ti為月平均氣溫;I為年熱量指數(shù);m為與I相關(guān)的參數(shù);K為緯度和月份校正系數(shù)。

    (2)

    (3)

    式中:當(dāng)月水分虧缺量D=P-PE。

    (4)

    式中:參數(shù)α、β和γ分別為尺度參數(shù)、形狀參數(shù)和位置參數(shù),采用線性矩法擬合得到。

    (5)

    式中:c0=2.515 517,c1=0.802 853,c2=0.010 328,d1=1.432 788,d2=0.189 269,d3=0.001 308。

    目前SPEI對(duì)干旱程度的等級(jí)劃分沒有統(tǒng)一標(biāo)準(zhǔn),綜合比較了多種等級(jí)劃分方法[13-16],本文采用的SPEI干濕等級(jí)分類如表1所示。

    2.2 趨勢(shì)分析

    Mann-Kendall法為國(guó)內(nèi)外學(xué)者廣泛應(yīng)用的非參數(shù)檢驗(yàn)方法[17-19],檢驗(yàn)統(tǒng)計(jì)量|Z|值越大表示增加或減少的趨勢(shì)越明顯,如果|Z|大于或等于1.64、1.96和2.57,分別表示序列通過了置信度為90%、95%和99%的顯著性趨勢(shì)檢驗(yàn)。Z為正值表示趨向濕潤(rùn),負(fù)值表示趨向干旱,且定義如下:①|(zhì)Z|≤1.64為無明顯趨勢(shì);② 1.64<|Z|≤1.96為微弱趨勢(shì);③1.96<|Z|≤2.57為顯著趨勢(shì);④|Z|>2.57為非常顯著趨勢(shì)。

    一般情況下,針對(duì)某一固定時(shí)間序列進(jìn)行趨勢(shì)檢驗(yàn),但可能由于所選序列時(shí)段不同導(dǎo)致趨勢(shì)結(jié)果發(fā)生變化,難以判斷時(shí)間序列真實(shí)變化趨勢(shì),為體現(xiàn)時(shí)段不同對(duì)序列趨勢(shì)檢測(cè)結(jié)果的影響,即判別序列變化趨勢(shì)是否穩(wěn)定,本文計(jì)算極端指數(shù)趨勢(shì)變化的穩(wěn)定性[20]:

    (6)

    式中:c為用于確定趨勢(shì)穩(wěn)定性的指標(biāo);M為站點(diǎn)個(gè)數(shù),本文針對(duì)單個(gè)站點(diǎn)進(jìn)行穩(wěn)定性分析,故M=1;K為全序列中滑動(dòng)30個(gè)值的子序列總數(shù),如全序列長(zhǎng)度為516,則K=487;Np為通過置信度為90%的顯著性檢驗(yàn)的子序列個(gè)數(shù)。

    本文對(duì)趨勢(shì)穩(wěn)定性的判別條件定義如下:①c≤25%時(shí),序列為不穩(wěn)定趨勢(shì)變化;②25%65%時(shí),序列為強(qiáng)穩(wěn)定的趨勢(shì)變化。

    3 結(jié)果分析

    3.1 SPEI多尺度序列分析

    SPEI時(shí)間尺度是前期影響時(shí)長(zhǎng)的表達(dá),尺度為1、3、6、12時(shí),表示考慮前期影響期為1個(gè)月、3個(gè)月、6個(gè)月、12個(gè)月。SPEI1和SPEI3分別受每月和季節(jié)降水氣溫變化影響,前期累積影響時(shí)段較小,可反映短時(shí)段內(nèi)土壤水分變化情況,在農(nóng)業(yè)灌溉方面具有一定的優(yōu)勢(shì);相對(duì)長(zhǎng)尺度的SPEI6和SPEI12可反映下層土壤水、地下徑流量和水庫(kù)蓄水情況。時(shí)間尺度越大,受降水氣溫影響時(shí)間越長(zhǎng)。在給定時(shí)間尺度上計(jì)算水分虧缺量累積概率,分析不同尺度SPEI值,可綜合評(píng)價(jià)研究區(qū)域旱澇情況。對(duì)24個(gè)站點(diǎn)各時(shí)間尺度的SPEI進(jìn)行算術(shù)平均計(jì)算,分析飛來峽流域整體干旱情況,結(jié)果見圖2。

    從圖2可以看出,短尺度下,SPEI1和SPEI3對(duì)短期內(nèi)降水和氣溫的變化反應(yīng)較為敏感,波動(dòng)性較大,充分反映出流域短期內(nèi)干濕頻繁交替特征。隨著時(shí)間尺度增加,SPEI6和SPEI12對(duì)降水和氣溫的響應(yīng)減緩,變化波動(dòng)性降低,能夠很好地反映前期降水和氣溫變化對(duì)干濕程度的累積影響,體現(xiàn)出旱澇的長(zhǎng)期變化特征。長(zhǎng)尺度的SPEI值波動(dòng)減少,主要是因?yàn)榍捌谟绊憰r(shí)長(zhǎng)增加后,短時(shí)間內(nèi)的降水氣溫變化不能徹底改變區(qū)域干濕程度。

    據(jù)統(tǒng)計(jì)分析,1990年代以來,華南地區(qū)進(jìn)入干旱災(zāi)害頻發(fā)期,2003年降水量嚴(yán)重偏少,在此影響下2004年汛期少雨,旱情進(jìn)一步嚴(yán)重,發(fā)生建國(guó)以來少有的干旱事件[21-22]。流域SPEI指數(shù)也有相應(yīng)體現(xiàn),2003-2004年的24個(gè)SPEI1值中有13個(gè)為輕度干旱以上級(jí)別,除了2003年7月達(dá)到了嚴(yán)重干旱級(jí)別,2003年2月和2004年9月達(dá)到中度干旱級(jí)別,其余的10個(gè)月均為輕度干旱級(jí)別。

    圖2 1969-2011年間飛來峽流域內(nèi)1、3、6和12個(gè)月時(shí)間尺度的SPEI動(dòng)態(tài)特征

    而SPEI3值則全為負(fù)值,其中有16個(gè)月達(dá)到輕度干旱以上級(jí)別,特大、嚴(yán)重、中度、輕度干旱的月份分別為0個(gè)、1個(gè)、6個(gè)、9個(gè),反映了該時(shí)段內(nèi)全區(qū)域偏干的特征。SPEI6顯示2003年1月-2004年9月受前期降水及氣溫影響,保持了一個(gè)較久的旱期,2004年10月開始轉(zhuǎn)為正常偏濕潤(rùn)狀態(tài),輕度干旱以上級(jí)別的月份增加至20個(gè),特大、嚴(yán)重、中度、輕度干旱的月份個(gè)數(shù)分別為0個(gè)、2個(gè)、7個(gè)、11個(gè)。SPEI12顯示出穩(wěn)定的干旱特征,多個(gè)月份達(dá)到了輕度干旱級(jí)別,其中特大、嚴(yán)重、中度、輕度干旱的月份個(gè)數(shù)分別為0個(gè)、1個(gè)、13個(gè)、3個(gè)。以上分析表明,隨著時(shí)間尺度增加,干旱程度有所提高,充分反映了前期水分虧缺的累積影響。

    干旱形成主要是由于短時(shí)段內(nèi)降水量減少、蒸發(fā)量增加而引起的土壤水分缺乏,或者長(zhǎng)期的補(bǔ)給水分缺乏[23]。本文所采用的SPEI能夠清楚地反映干濕變化不同方面的影響,多尺度綜合應(yīng)用可實(shí)現(xiàn)對(duì)飛來峽流域干旱的綜合監(jiān)測(cè)評(píng)價(jià)。

    3.2 干旱頻次變化

    為分析流域干旱事件在各年代的頻次變化情況,取1969-2011年43年的SPEI值,根據(jù)表1統(tǒng)計(jì)不同年代間飛來峽流域的干旱事件發(fā)生頻次,如表2所示。

    表1 SPEI干旱等級(jí)分類

    表2 飛來峽流域不同年代干旱事件頻次

    飛來峽流域4個(gè)時(shí)間尺度的SPEI均顯示21世紀(jì)初的干旱事件頻次最多,其次為1980年代。年代干旱累計(jì)頻次在22%~40%范圍內(nèi)波動(dòng)上升,增加幅度較大,長(zhǎng)時(shí)間尺度的增幅更大。影響累計(jì)頻次增加的主要因子隨著時(shí)間尺度增加,由輕度干旱向嚴(yán)重干旱發(fā)展,1個(gè)月尺度的累計(jì)頻次變化主要是輕度干旱頻次變化,3個(gè)月尺度的主要是中度干旱和輕度干旱的變化,6個(gè)月和12個(gè)月尺度的主要是中度干旱和嚴(yán)重干旱的變化。

    特大干旱往往給農(nóng)業(yè)生產(chǎn)和社會(huì)生活造成巨大影響,定義輕度干旱以上為干旱事件,統(tǒng)計(jì)特大干旱占所有干旱事件的比例變化情況,如表3所示。

    表3 1969-2011年特大干旱占干旱事件比例變化

    由表3可以看出,4個(gè)尺度的特大干旱比例不超過4%,1個(gè)月和3個(gè)月尺度的特旱比例趨勢(shì)變化不顯著,6個(gè)月和12個(gè)月尺度特旱比例的M-K趨勢(shì)檢驗(yàn)值通過了置信度為90%的顯著性檢驗(yàn),增加趨勢(shì)顯著,也就是說飛來峽流域特大干旱所占比例顯著增加。

    3.3 趨勢(shì)變化空間分布

    分析SPEI趨勢(shì)變化對(duì)于揭示飛來峽流域內(nèi)干濕變化情況具有重要參考價(jià)值,本文基于M-K趨勢(shì)分析法,計(jì)算了43年來24個(gè)站點(diǎn)的SPEI1、SPEI3、SPEI6、SPEI12的趨勢(shì)變化特征,不同尺度的SPEI趨勢(shì)變化的空間分布情況如圖3所示。

    M-K檢驗(yàn)結(jié)果顯示整個(gè)飛來峽流域的SPEI1整體上呈由東北向西南遞減且無濕潤(rùn)趨勢(shì),流域下游多個(gè)站點(diǎn)呈現(xiàn)顯著下降趨勢(shì),靠近流域出口斷面的高道和紅橋雨量站的SPEI1通過了置信度為99%的顯著性檢驗(yàn),呈現(xiàn)非常顯著的干旱趨勢(shì),也就是說飛來峽下游土壤含水量顯著減少。

    SPEI3的趨勢(shì)變化與SPEI1類似,24個(gè)站點(diǎn)均無上升趨勢(shì),流域中部地區(qū)無明顯變化,上游和下游多個(gè)站點(diǎn)的SPEI3顯示出不同程度的下降趨勢(shì),其中下游靠近流域出口斷面的站點(diǎn)表現(xiàn)出較為顯著的干旱趨勢(shì),由此可以看出受前期降水氣溫變化影響,流域上游和下游季尺度的干旱程度有所增加。

    隨著時(shí)間尺度增加,飛來峽流域整體上顯示出更加明顯穩(wěn)定的趨勢(shì)。SPEI6在流域上部和出口斷面附近表現(xiàn)出不同程度的下降趨勢(shì),它主要反映深層土壤水和地下徑流的干濕狀況,其空間趨勢(shì)變化特征體現(xiàn)出流域內(nèi)受補(bǔ)給不足的影響,上游和出口斷面的干旱加劇。

    圖3 飛來峽流域不同尺度SPEI指數(shù)的變化趨勢(shì)空間分布圖

    較之小尺度的SPEI值,SPEI12在上游與出口斷面區(qū)域的下降趨勢(shì)更顯著,包括樂昌、紅橋在內(nèi)的14個(gè)站點(diǎn)的下降趨勢(shì)通過了置信度水平為99%的顯著性檢驗(yàn),這說明在長(zhǎng)時(shí)間尺度上流域上游和下游干旱程度增加趨勢(shì)非常顯著,這對(duì)飛來峽水利樞紐的水庫(kù)蓄水及周邊生態(tài)環(huán)境產(chǎn)生一定的負(fù)面影響。

    3.4 趨勢(shì)穩(wěn)定性分析

    變化趨勢(shì)穩(wěn)定即波動(dòng)性較小,在短期內(nèi)干濕變化不大。根據(jù)定義對(duì)4個(gè)尺度的SPEI進(jìn)行穩(wěn)定等級(jí)劃分,其空間分布如圖4所示。

    圖4 飛來峽流域不同尺度SPEI指數(shù)的趨勢(shì)穩(wěn)定性空間分布圖

    由圖4可以看出,SPEI1和SPEI3反映流域內(nèi)土壤含水量變化趨勢(shì)的穩(wěn)定性不佳,在月尺度和季尺度上出現(xiàn)較大波動(dòng),同時(shí)流域上下游的穩(wěn)定性沒有明顯差別,空間分布較為均勻。SPEI6顯示地下徑流的干濕趨勢(shì)穩(wěn)定性存在空間差異,呈現(xiàn)出由東南向西北穩(wěn)定性遞增的分布特征。受前期水分虧缺影響,SPEI12顯示全流域24個(gè)站點(diǎn)均表現(xiàn)出穩(wěn)定的干濕趨勢(shì),也就是說水庫(kù)蓄水趨勢(shì)變化較為穩(wěn)定。總的來說,隨著時(shí)間尺度增加,飛來峽流域干濕變化趨勢(shì)穩(wěn)定性有所增加,但沒有出現(xiàn)強(qiáng)穩(wěn)定的變化情況。

    4 結(jié)論

    本文利用飛來峽流域內(nèi)1969-2011年24個(gè)雨量站的逐月降水量和4個(gè)氣象站的月平均氣溫,應(yīng)用標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI),綜合分析了飛來峽流域干旱時(shí)空演變特征,得出以下結(jié)論。

    (1)SPEI能較好地表征流域內(nèi)的干旱變化特征,短時(shí)間尺度的SPEI對(duì)短期內(nèi)干濕變化反應(yīng)靈敏,長(zhǎng)時(shí)間尺度的SPEI反映前期累積影響,可體現(xiàn)旱澇的長(zhǎng)期變化特征。

    (2)在時(shí)間上,飛來峽流域的干旱事件頻次波動(dòng)上升,長(zhǎng)時(shí)間尺度干旱事件增加更快,以中度和嚴(yán)重干旱事件的增加為主,隨著全球氣溫的持續(xù)升高,特大干旱事件所占比例總體上顯著增加,說明流域內(nèi)干旱問題突出。

    (3)在空間上,流域呈現(xiàn)上游和出口斷面附近趨向干旱的特征,流域中部地區(qū)無明顯的干濕趨勢(shì)變化,長(zhǎng)時(shí)間尺度的干旱趨勢(shì)非常顯著;短時(shí)間尺度的干濕變化趨勢(shì)穩(wěn)定性不佳,隨著時(shí)間尺度增加,穩(wěn)定程度有所增加,但未出現(xiàn)強(qiáng)穩(wěn)定情況。

    [1] Mishra A K, Singh V P. A review of drought concepts[J]. Journal of Hydrology, 2010, 391: 202-216.

    [2] Palmer W C. Meteorological drought[R].US: Weather Bureau,1965.

    [3] Guttman N B. Comparing the palmer drought index and the standardized precipitation index[J]. Journal of The American Water Resources Association, 1998, 34(1):113-121.

    [4] Mckee T B, Doesken N J, Kleist J. The relationship of drought frequency and duration to time scales[C]//Eighth Conference on Applied Climatology. Anaheim, CA, American meteorological society, 1993:179-184.

    [5] Vicente-Serrano S M, Begueriá S, López-Moreno J I. A multiscalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index[J]. Journal of Climate, 2010(23):1696-1718.

    [6] Vicente-Serrano S M, Begueriá S, López-Moreno J I, et al. A new global 0.5° gridded dataset (1901-2006) of a multiscalar drought index Comparison with current drought index datasets based on the Palmer Drought Severity Index[J]. Notes and Correspondence, 2011(11):1033-1043.

    [7] López-Moreno J I, Vicente-Serranoa S M, Zabalzaa J, et al. Hydrological response to climate variability at different time scales A study in the Ebro basin[J]. Journal of Hydrology, 2013(477):175-188.

    [8] Potop V, Mozny M, Soukup J. Drought evolution at various time scales in the lowland regions and their impact on vegetable crops in the Czech Republic[J]. Agricultural and Forest Meteorology, 2012(156):121-133.

    [9] 莊少偉, 左洪超, 任鵬程, 等. 標(biāo)準(zhǔn)化降水蒸發(fā)指數(shù)在中國(guó)區(qū)域的應(yīng)用[J]. 氣候與環(huán)境研究, 2013, 18(5):617-625.

    [10]李偉光, 易雪, 侯美亭, 等. 基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的中國(guó)干旱趨勢(shì)研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2012, 20(5):643-649.

    [11]李偉光, 侯美亭, 陳匯林, 等. 基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的華南干旱趨勢(shì)研究[J]. 自然災(zāi)害學(xué)報(bào), 2012, 21(4):84-90.

    [12]郭晶. 廣東省干旱動(dòng)態(tài)監(jiān)測(cè)及其影響定量評(píng)估研究[D]. 南京:南京信息工程大學(xué), 2008.

    [13]陳子燊, 黃強(qiáng), 劉曾美. 1962-2007年廣東干濕時(shí)空變化特征分析[J]. 水科學(xué)進(jìn)展, 2013, 24(4):469-476.

    [14]劉占明, 陳子燊, 黃強(qiáng), 等. 7種干旱評(píng)估指標(biāo)在廣東北江流域應(yīng)用中的對(duì)比分析[J]. 資源科學(xué), 2013, 35(5):1007-1015.

    [15]熊光潔, 張博凱, 李崇銀, 等. 基于SPEI的中國(guó)西南地區(qū)1961-2012年干旱變化特征分析[J]. 氣候變化研究進(jìn)展, 2013, 9(3):192-198.

    [16]王勁廷. 濕潤(rùn)指數(shù)在我國(guó)江淮流域的適應(yīng)性評(píng)估[D]. 南京:南京信息工程大學(xué), 2012.

    [17]鄭杰元, 黃國(guó)如, 王質(zhì)軍, 等. 廣州市近年降雨時(shí)空變化規(guī)律分析[J]. 水電能源科學(xué), 2011, 29(3):5-8.

    [18]康淑媛, 張勃, 柳景峰, 等. 基于Mann-Kendall法的張掖市降水量時(shí)空分布規(guī)律分析[J]. 資源科學(xué), 2009, 31(3):501-508.

    [19]Wang L N, Shao Q X, Chen X H, et al. Flood changes during the past 50 years in Wujiang River, South China[J]. Hydrological Processes, 2012(26):3561-3569.

    [20]Lupikasza E. Spatial and temporal variability of extreme precipitation in Poland in the period 1951-2006[J]. International Journal of Climatology,2006,30(7):991-1007.

    [21]李曉娟, 曾沁, 梁健, 等. 華南地區(qū)干旱氣候預(yù)測(cè)研究[J]. 氣象科技, 2007, 35(1):26-30.

    [22]黃晚華, 楊曉光, 李茂松, 等. 基于標(biāo)準(zhǔn)化降水指數(shù)的中國(guó)南方季節(jié)性干旱近58a演變特征[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2010, 26(7):50-59.

    [23]李偉光, 陳匯林, 朱乃海, 等. 標(biāo)準(zhǔn)化降水指標(biāo)在海南島干旱監(jiān)測(cè)中的應(yīng)用分析[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2009, 17(1):178-182.

    Analysis on Spatial and Temporal Characteristics of Drought in theFeilaixia Catchment of Beijiang River Based on SPEI

    Xian Zhuoyan1, Wu Chuanhao1and Huang Guoru1, 2

    (1.SchoolofCivilEngineeringandTransportation,SouthChinaUniversityofTechnology,Guangzhou510640,China; 2.StateKeyLaboratoryofSubtropicalBuildingScience,SouthChinaUniversityofTechnology,Guangzhou510640,China)

    Droughteventsoccurredfrequentlybecauseofclimatewarming.Basedonmonthlyprecipitationandtemperaturedatafrom28stations(24rainfallstationsand4meteorologicalstations)duringtheperiod1969-2011intheFeilaixiacatchmentoftheBeijiangRiverbasin,thestandardizedprecipitationevapotranspirationindex(SPEI)iscalculated.Themulti-scalespatialandtemporalcharacteristicsofdroughtareanalyzedusingtheMann-Kendalltrendtestmethod.TheresultsshowthattheSPEIisasuitableindextodescribethevariationofwetanddryspells.Overall,thefrequencyofdroughtandtheproportionofseveredroughtshowincreasingtrendsinthestudyarea.Asignificanttrendofdroughtisfoundintheupstreamareaanddownstreamareaneartheoutlet.Furthermore,thetrendsofdroughtbecomemoresignificantandmorestablewiththeincreaseoftimescales.

    drought;standardizedprecipitationevapotranspirationindex(SPEI);Feilaixiacatchment;spatialandtemporalcharacteristics;stabilityoftrends

    2014-11-07

    2015-01-12

    水利部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201301093, 201401048);華南理工大學(xué)亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室自主研究課題項(xiàng)目(2014ZC09)

    冼卓雁(1990-),女,廣西桂平人,碩士研究生,主要從事水文水資源研究. E-mail:867325059@qq.com

    黃國(guó)如(1969-),男,江蘇南京人,教授,博士生導(dǎo)師,主要從事水文水資源研究. E-mail: huanggr@scut.edu.cn

    X43;S423

    A

    1000-811X(2015)03-0198-06

    10.3969/j.issn.1000-811X.2015.03.036

    冼卓雁,武傳號(hào),黃國(guó)如. 基于SPEI的北江飛來峽流域干旱時(shí)空演變特征分析[J].災(zāi)害學(xué), 2015,30(3):198-203. [Xian Zhuoyan, Wu Chuanhao and Huang Guoru. Analysis on spatial and temporal characteristics of drought in the Feilaixia catchment of the Beijiang River based on SPEI [J].Journal of Catastrophology, 2015,30(3):198-203.]

    猜你喜歡
    北江時(shí)間尺度尺度
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    北江,向前
    奔騰北江
    泥娃娃
    暢談(2018年11期)2018-08-26 02:10:28
    大連市暴雨多時(shí)間尺度研究分析
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于MODIS北江流域土地覆蓋變化研究及CA-Markov預(yù)測(cè)
    河北遙感(2015年2期)2015-07-18 11:11:14
    国产一区二区在线观看日韩| 久久久久九九精品影院| 听说在线观看完整版免费高清| 麻豆国产97在线/欧美| 久久婷婷人人爽人人干人人爱| 男女下面进入的视频免费午夜| 五月伊人婷婷丁香| 午夜a级毛片| 99热这里只有是精品50| 岛国在线免费视频观看| 狠狠狠狠99中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 免费无遮挡裸体视频| 欧美日韩国产亚洲二区| 国产午夜精品论理片| 免费av观看视频| 91精品一卡2卡3卡4卡| 亚洲最大成人中文| 看十八女毛片水多多多| 免费av不卡在线播放| 免费av不卡在线播放| 日本黄色视频三级网站网址| 亚洲欧美精品专区久久| 看非洲黑人一级黄片| 91麻豆精品激情在线观看国产| avwww免费| 日韩av不卡免费在线播放| 人人妻人人澡人人爽人人夜夜 | 成人亚洲精品av一区二区| 老司机福利观看| 日韩欧美国产在线观看| 中文字幕久久专区| 亚洲在久久综合| 亚洲国产高清在线一区二区三| 国产成人a区在线观看| 色综合色国产| 国产一区二区三区av在线 | 国产精品乱码一区二三区的特点| 欧美bdsm另类| 夜夜爽天天搞| 国产 一区 欧美 日韩| 美女脱内裤让男人舔精品视频 | 麻豆精品久久久久久蜜桃| 国产激情偷乱视频一区二区| 午夜福利在线在线| 在线观看美女被高潮喷水网站| 色视频www国产| 国产午夜精品一二区理论片| 国产精品福利在线免费观看| 欧美日韩乱码在线| 亚洲精品乱码久久久v下载方式| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 亚洲精品影视一区二区三区av| 中出人妻视频一区二区| 黄色配什么色好看| 亚洲欧美成人综合另类久久久 | 老司机福利观看| 黄色配什么色好看| 99久久中文字幕三级久久日本| 性色avwww在线观看| 久久精品国产鲁丝片午夜精品| 久久久久久久久久久免费av| 搡女人真爽免费视频火全软件| 成人午夜高清在线视频| 看黄色毛片网站| 99久久久亚洲精品蜜臀av| 国产高清视频在线观看网站| 国产精品美女特级片免费视频播放器| 中文字幕av在线有码专区| 日本免费一区二区三区高清不卡| 国产精品爽爽va在线观看网站| 狂野欧美白嫩少妇大欣赏| 日韩一区二区视频免费看| 91久久精品电影网| 国产一区二区在线av高清观看| 亚洲国产精品成人久久小说 | av在线播放精品| 国产精品久久久久久久久免| 亚洲经典国产精华液单| 99热网站在线观看| 国产成人影院久久av| 又爽又黄a免费视频| 国产精品1区2区在线观看.| 亚洲欧美成人综合另类久久久 | 观看免费一级毛片| 不卡一级毛片| 淫秽高清视频在线观看| 中文字幕制服av| 日本熟妇午夜| 国产成年人精品一区二区| 男人的好看免费观看在线视频| 亚洲在线观看片| 精品一区二区三区人妻视频| 97在线视频观看| 国产成人aa在线观看| 国产在线男女| eeuss影院久久| 狠狠狠狠99中文字幕| 内地一区二区视频在线| 亚洲成人av在线免费| 久久久成人免费电影| 亚洲高清免费不卡视频| 亚洲欧美日韩高清在线视频| 91aial.com中文字幕在线观看| 天美传媒精品一区二区| 午夜福利成人在线免费观看| 久久久久国产网址| 精品午夜福利在线看| 1024手机看黄色片| 国产精品不卡视频一区二区| 2021天堂中文幕一二区在线观| 日韩强制内射视频| 日本黄色片子视频| 三级经典国产精品| 五月玫瑰六月丁香| 性欧美人与动物交配| 免费看光身美女| 婷婷六月久久综合丁香| 国产老妇伦熟女老妇高清| 欧美最黄视频在线播放免费| 2022亚洲国产成人精品| 好男人视频免费观看在线| 亚洲av成人精品一区久久| 亚洲中文字幕一区二区三区有码在线看| 精品久久国产蜜桃| 美女高潮的动态| 国产精品99久久久久久久久| av免费在线看不卡| 欧美xxxx性猛交bbbb| 欧美人与善性xxx| 97人妻精品一区二区三区麻豆| 青春草亚洲视频在线观看| 久久精品久久久久久久性| 天天一区二区日本电影三级| 亚洲av免费在线观看| 久久久久网色| 大香蕉久久网| 久久人人爽人人爽人人片va| 晚上一个人看的免费电影| 少妇高潮的动态图| 色哟哟·www| 少妇人妻精品综合一区二区 | 国产一区二区激情短视频| 亚洲成人久久爱视频| 97超碰精品成人国产| 日本与韩国留学比较| 亚洲av一区综合| 日本在线视频免费播放| 欧美精品一区二区大全| a级毛片a级免费在线| 婷婷色av中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美清纯卡通| 99精品在免费线老司机午夜| 日韩精品有码人妻一区| 国产真实乱freesex| 久久亚洲国产成人精品v| www日本黄色视频网| 国模一区二区三区四区视频| 亚洲无线观看免费| 美女大奶头视频| 欧美一区二区国产精品久久精品| 中文欧美无线码| 亚洲电影在线观看av| 最好的美女福利视频网| 国产私拍福利视频在线观看| 久久久久久久午夜电影| а√天堂www在线а√下载| 亚洲欧美日韩高清专用| 亚洲第一电影网av| 三级经典国产精品| 久久亚洲精品不卡| av在线播放精品| 国产v大片淫在线免费观看| 人人妻人人看人人澡| 寂寞人妻少妇视频99o| 在线免费观看的www视频| 亚洲五月天丁香| 成熟少妇高潮喷水视频| av在线亚洲专区| av在线观看视频网站免费| 欧美性猛交黑人性爽| 国产黄色视频一区二区在线观看 | 午夜精品一区二区三区免费看| 欧美三级亚洲精品| 久久久久久伊人网av| 亚洲av熟女| 又爽又黄无遮挡网站| 人人妻人人看人人澡| 国产成人精品久久久久久| 精品一区二区免费观看| a级毛片a级免费在线| 久99久视频精品免费| 中文字幕久久专区| 成年女人看的毛片在线观看| 人妻系列 视频| a级一级毛片免费在线观看| 三级国产精品欧美在线观看| 爱豆传媒免费全集在线观看| 亚洲婷婷狠狠爱综合网| 真实男女啪啪啪动态图| 亚洲一区二区三区色噜噜| 亚洲av成人av| 99riav亚洲国产免费| 欧美成人免费av一区二区三区| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 日日啪夜夜撸| 久久久久久久午夜电影| 别揉我奶头 嗯啊视频| 欧美一区二区精品小视频在线| 成人午夜高清在线视频| 天堂av国产一区二区熟女人妻| 日韩中字成人| 午夜激情欧美在线| 身体一侧抽搐| 性插视频无遮挡在线免费观看| 真实男女啪啪啪动态图| 精品熟女少妇av免费看| 哪个播放器可以免费观看大片| 国国产精品蜜臀av免费| 亚洲五月天丁香| 久久精品久久久久久久性| 亚洲性久久影院| 国产真实伦视频高清在线观看| 久久久精品欧美日韩精品| 在线免费观看不下载黄p国产| 欧美变态另类bdsm刘玥| 在线观看av片永久免费下载| a级毛片免费高清观看在线播放| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| 在线观看av片永久免费下载| 一级毛片电影观看 | 亚洲国产欧美在线一区| 亚洲精品影视一区二区三区av| 国产91av在线免费观看| 亚洲无线观看免费| 久久人人精品亚洲av| 欧美一区二区亚洲| 成人鲁丝片一二三区免费| 亚洲最大成人中文| 九色成人免费人妻av| 午夜精品国产一区二区电影 | 欧美三级亚洲精品| 久久久成人免费电影| 午夜老司机福利剧场| 亚洲在线自拍视频| 亚洲综合色惰| 亚洲不卡免费看| 久久99精品国语久久久| 男插女下体视频免费在线播放| 日本熟妇午夜| 免费大片18禁| 成人午夜精彩视频在线观看| 亚洲av免费在线观看| 亚洲av不卡在线观看| 久久久久久久久久久免费av| 国产午夜精品久久久久久一区二区三区| 国国产精品蜜臀av免费| 在线免费十八禁| av黄色大香蕉| 国产熟女欧美一区二区| 亚洲av免费在线观看| 成人av在线播放网站| 黄色欧美视频在线观看| 九九久久精品国产亚洲av麻豆| 非洲黑人性xxxx精品又粗又长| 日本成人三级电影网站| 国产一区二区亚洲精品在线观看| 一级毛片我不卡| 国产v大片淫在线免费观看| 99热这里只有精品一区| 熟妇人妻久久中文字幕3abv| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 亚洲天堂国产精品一区在线| 日本黄大片高清| 久久久精品94久久精品| 色视频www国产| 大香蕉久久网| 日本av手机在线免费观看| 久久久精品大字幕| 日产精品乱码卡一卡2卡三| 免费看美女性在线毛片视频| 久久精品久久久久久久性| 中国美白少妇内射xxxbb| 国产极品天堂在线| 国产精品三级大全| 秋霞在线观看毛片| 国产精品一区二区三区四区免费观看| 国产精品一及| 久久久久网色| 欧洲精品卡2卡3卡4卡5卡区| 黄色日韩在线| 精品午夜福利在线看| 久久国内精品自在自线图片| 久久精品国产鲁丝片午夜精品| 国产老妇女一区| 麻豆国产97在线/欧美| 亚洲成a人片在线一区二区| 亚洲最大成人中文| 高清日韩中文字幕在线| 国产麻豆成人av免费视频| 久久亚洲国产成人精品v| 热99在线观看视频| 国产视频内射| 在线观看av片永久免费下载| 一区二区三区四区激情视频 | 日韩精品青青久久久久久| 嫩草影院精品99| 嘟嘟电影网在线观看| 大香蕉久久网| 深夜精品福利| 久久精品国产亚洲网站| 久久鲁丝午夜福利片| 美女黄网站色视频| 麻豆精品久久久久久蜜桃| 99视频精品全部免费 在线| 一区二区三区四区激情视频 | 日日干狠狠操夜夜爽| 国产av一区在线观看免费| 成人毛片60女人毛片免费| 男女边吃奶边做爰视频| 中文字幕精品亚洲无线码一区| 国产午夜精品论理片| 久久国产乱子免费精品| 午夜免费激情av| 日韩欧美在线乱码| 特级一级黄色大片| 国产精品精品国产色婷婷| 99精品在免费线老司机午夜| 精品人妻视频免费看| 美女内射精品一级片tv| 国产大屁股一区二区在线视频| 免费观看在线日韩| 欧美性感艳星| 免费看av在线观看网站| 99视频精品全部免费 在线| 噜噜噜噜噜久久久久久91| 一卡2卡三卡四卡精品乱码亚洲| 中文亚洲av片在线观看爽| 欧美一区二区亚洲| 老女人水多毛片| 日韩欧美一区二区三区在线观看| 色尼玛亚洲综合影院| 亚洲精品日韩av片在线观看| 精品一区二区三区视频在线| 亚洲性久久影院| eeuss影院久久| 国产一区二区在线观看日韩| 在线观看美女被高潮喷水网站| 观看免费一级毛片| 亚洲欧美精品专区久久| 男插女下体视频免费在线播放| 亚洲成av人片在线播放无| 我的女老师完整版在线观看| 成人鲁丝片一二三区免费| av在线观看视频网站免费| 国产精品一区www在线观看| 国产毛片a区久久久久| 国产成人福利小说| 最近2019中文字幕mv第一页| kizo精华| 亚洲欧美成人综合另类久久久 | 91久久精品国产一区二区成人| 黄色欧美视频在线观看| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 婷婷亚洲欧美| 少妇裸体淫交视频免费看高清| 国产伦一二天堂av在线观看| 日本免费a在线| 免费av观看视频| 高清日韩中文字幕在线| 精品无人区乱码1区二区| 国产精品99久久久久久久久| 99热网站在线观看| 国产人妻一区二区三区在| 亚洲经典国产精华液单| 男女边吃奶边做爰视频| 亚洲图色成人| 欧美+亚洲+日韩+国产| 欧美xxxx黑人xx丫x性爽| 麻豆成人av视频| 国产国拍精品亚洲av在线观看| 99热精品在线国产| a级毛片a级免费在线| 欧美日本视频| 深夜精品福利| 国产精品伦人一区二区| 一级黄色大片毛片| 亚洲自偷自拍三级| 乱系列少妇在线播放| 在线观看66精品国产| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放 | 一本一本综合久久| 亚洲第一区二区三区不卡| 1024手机看黄色片| av天堂中文字幕网| 国产精华一区二区三区| 22中文网久久字幕| 可以在线观看的亚洲视频| 人妻制服诱惑在线中文字幕| АⅤ资源中文在线天堂| 久久精品久久久久久久性| 国产高清不卡午夜福利| 床上黄色一级片| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 最近的中文字幕免费完整| ponron亚洲| 少妇熟女欧美另类| 婷婷色综合大香蕉| 国产大屁股一区二区在线视频| 国产一区亚洲一区在线观看| 精品久久久噜噜| 国国产精品蜜臀av免费| 国产免费男女视频| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久久免费av| 国产精品女同一区二区软件| 麻豆成人午夜福利视频| 日本熟妇午夜| 欧美性猛交黑人性爽| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线 | 3wmmmm亚洲av在线观看| 中文在线观看免费www的网站| 国产高清视频在线观看网站| kizo精华| 91久久精品国产一区二区成人| 精品不卡国产一区二区三区| av女优亚洲男人天堂| 久久欧美精品欧美久久欧美| 欧美成人免费av一区二区三区| 听说在线观看完整版免费高清| 久久这里有精品视频免费| 欧美一区二区国产精品久久精品| 内射极品少妇av片p| 精品免费久久久久久久清纯| 精品不卡国产一区二区三区| 成人av在线播放网站| 日本免费一区二区三区高清不卡| 春色校园在线视频观看| 少妇高潮的动态图| 久久久精品欧美日韩精品| 日本一本二区三区精品| eeuss影院久久| 亚洲一区高清亚洲精品| 五月伊人婷婷丁香| av在线老鸭窝| 亚洲最大成人av| 看免费成人av毛片| 久久国内精品自在自线图片| 欧美又色又爽又黄视频| 国产精品久久电影中文字幕| 亚洲熟妇中文字幕五十中出| av在线亚洲专区| 99久久无色码亚洲精品果冻| 麻豆久久精品国产亚洲av| 99热全是精品| 日韩中字成人| 国产久久久一区二区三区| 日韩中字成人| 女人被狂操c到高潮| 22中文网久久字幕| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 一个人观看的视频www高清免费观看| kizo精华| 久久精品综合一区二区三区| 国产综合懂色| 国产亚洲91精品色在线| 在线播放无遮挡| 一个人观看的视频www高清免费观看| 成年女人看的毛片在线观看| 国产视频内射| 国产成人91sexporn| 成人欧美大片| 岛国在线免费视频观看| 寂寞人妻少妇视频99o| 国产视频内射| 国产三级在线视频| 国产精品美女特级片免费视频播放器| 国产午夜福利久久久久久| 国产成人a∨麻豆精品| h日本视频在线播放| 亚洲av第一区精品v没综合| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 特级一级黄色大片| 国产亚洲5aaaaa淫片| 夜夜夜夜夜久久久久| 精品一区二区免费观看| 赤兔流量卡办理| 亚洲aⅴ乱码一区二区在线播放| 国产极品精品免费视频能看的| 麻豆精品久久久久久蜜桃| 热99在线观看视频| 久久这里有精品视频免费| 国产精品三级大全| 看片在线看免费视频| 亚洲国产欧美人成| 大又大粗又爽又黄少妇毛片口| 免费不卡的大黄色大毛片视频在线观看 | 成人一区二区视频在线观看| 九九热线精品视视频播放| 内射极品少妇av片p| 一边摸一边抽搐一进一小说| 国产精品久久久久久亚洲av鲁大| 亚洲图色成人| 亚洲精品国产成人久久av| 久久久久久久久久久免费av| 国产成人aa在线观看| 在线观看一区二区三区| 国产色爽女视频免费观看| 久久久久久大精品| 搞女人的毛片| 亚洲av熟女| 成年版毛片免费区| 99国产精品一区二区蜜桃av| 欧美激情久久久久久爽电影| 精品免费久久久久久久清纯| 日韩,欧美,国产一区二区三区 | 欧美成人a在线观看| 欧美精品一区二区大全| 一卡2卡三卡四卡精品乱码亚洲| АⅤ资源中文在线天堂| 国内精品一区二区在线观看| 在现免费观看毛片| 亚洲成人中文字幕在线播放| 亚洲三级黄色毛片| 亚洲av第一区精品v没综合| 天堂中文最新版在线下载 | 麻豆成人av视频| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 麻豆国产97在线/欧美| 成人欧美大片| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 中文字幕熟女人妻在线| 日本成人三级电影网站| 久久99热6这里只有精品| 日韩欧美一区二区三区在线观看| 97超碰精品成人国产| 日韩欧美精品免费久久| 2022亚洲国产成人精品| 老女人水多毛片| 观看美女的网站| 我的老师免费观看完整版| 欧美极品一区二区三区四区| 久久久久久久午夜电影| 日本成人三级电影网站| 美女高潮的动态| 精品一区二区免费观看| 最近最新中文字幕大全电影3| 97超视频在线观看视频| 国产精品麻豆人妻色哟哟久久 | 国产高清三级在线| 国产成人精品婷婷| av专区在线播放| 91久久精品国产一区二区成人| 国产乱人视频| 国产在视频线在精品| 亚洲国产欧美人成| 欧美日韩综合久久久久久| 麻豆av噜噜一区二区三区| 国产高清视频在线观看网站| 国产精品一区二区三区四区免费观看| 国产麻豆成人av免费视频| 久久久久久伊人网av| 久久久久久久久久成人| 久久久国产成人免费| 久久亚洲精品不卡| 国产精品麻豆人妻色哟哟久久 | 岛国毛片在线播放| 国产成人a∨麻豆精品| 全区人妻精品视频| 99久久精品国产国产毛片| 99久久久亚洲精品蜜臀av| 久久草成人影院| 精品久久国产蜜桃| 一级黄色大片毛片| 成人二区视频| 成年女人永久免费观看视频| 欧美日韩综合久久久久久| 久久久国产成人精品二区| 精品久久久久久久久av| 国产老妇女一区| 亚洲精品456在线播放app| 国国产精品蜜臀av免费| 高清午夜精品一区二区三区 | 亚洲av免费在线观看| 黄色一级大片看看| 免费观看的影片在线观看| 此物有八面人人有两片| 国产精品伦人一区二区| 成人国产麻豆网| av在线老鸭窝| 国国产精品蜜臀av免费| 国产伦一二天堂av在线观看| 性插视频无遮挡在线免费观看| 99视频精品全部免费 在线| 熟女人妻精品中文字幕| 久久精品国产自在天天线| kizo精华| 精品一区二区三区视频在线| 日韩欧美一区二区三区在线观看| 特级一级黄色大片| 一进一出抽搐动态| 亚州av有码|