• <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
    国产一区二区激情短视频 | 成人影院久久| a 毛片基地| 熟妇人妻不卡中文字幕| 高清视频免费观看一区二区| 91aial.com中文字幕在线观看| 少妇的丰满在线观看| 亚洲,一卡二卡三卡| 国产乱来视频区| 成人综合一区亚洲| 男的添女的下面高潮视频| 中国三级夫妇交换| 2021少妇久久久久久久久久久| 熟妇人妻不卡中文字幕| 国产av精品麻豆| 亚洲欧洲国产日韩| 99九九在线精品视频| 丰满饥渴人妻一区二区三| 亚洲成人一二三区av| 国产永久视频网站| 欧美亚洲 丝袜 人妻 在线| 婷婷色麻豆天堂久久| 免费观看av网站的网址| 精品人妻熟女毛片av久久网站| 亚洲精品久久成人aⅴ小说| 免费日韩欧美在线观看| 久久ye,这里只有精品| 亚洲av福利一区| 国产欧美日韩综合在线一区二区| 成人毛片60女人毛片免费| 免费观看av网站的网址| 成人亚洲欧美一区二区av| 黄色配什么色好看| 日韩,欧美,国产一区二区三区| 女性生殖器流出的白浆| 亚洲经典国产精华液单| 人妻一区二区av| 在线精品无人区一区二区三| 国产在线免费精品| 欧美激情极品国产一区二区三区 | 亚洲欧美清纯卡通| 精品一区二区三卡| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 2021少妇久久久久久久久久久| 在线观看人妻少妇| 天天操日日干夜夜撸| xxx大片免费视频| av电影中文网址| videos熟女内射| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 在现免费观看毛片| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 男女国产视频网站| 亚洲av综合色区一区| 人妻一区二区av| 国产精品无大码| 国产综合精华液| 妹子高潮喷水视频| 国产精品.久久久| 看非洲黑人一级黄片| 久久久久久久精品精品| 亚洲图色成人| 久久人妻熟女aⅴ| 爱豆传媒免费全集在线观看| 欧美日韩av久久| 精品视频人人做人人爽| freevideosex欧美| 免费观看性生交大片5| 999精品在线视频| 国产探花极品一区二区| 国产国语露脸激情在线看| 国产免费视频播放在线视频| 国产高清三级在线| 高清黄色对白视频在线免费看| 自拍欧美九色日韩亚洲蝌蚪91| 色婷婷久久久亚洲欧美| kizo精华| av福利片在线| 春色校园在线视频观看| 一区二区日韩欧美中文字幕 | 咕卡用的链子| 青春草视频在线免费观看| 久久99一区二区三区| 中文精品一卡2卡3卡4更新| 亚洲精品日本国产第一区| 日韩精品免费视频一区二区三区 | 丰满少妇做爰视频| 久久婷婷青草| 制服诱惑二区| 9191精品国产免费久久| 亚洲三级黄色毛片| 超色免费av| 国产精品女同一区二区软件| 精品亚洲乱码少妇综合久久| 最近最新中文字幕大全免费视频 | 午夜福利,免费看| av在线观看视频网站免费| 综合色丁香网| 亚洲精品一区蜜桃| 人成视频在线观看免费观看| 中国国产av一级| 国产一级毛片在线| 国产av精品麻豆| 久久久久视频综合| 色婷婷久久久亚洲欧美| 97超碰精品成人国产| av播播在线观看一区| 国产 精品1| 午夜福利视频在线观看免费| a级毛片在线看网站| 欧美精品av麻豆av| 狠狠精品人妻久久久久久综合| 色94色欧美一区二区| 亚洲成国产人片在线观看| 99久国产av精品国产电影| 亚洲av在线观看美女高潮| 夫妻性生交免费视频一级片| 亚洲内射少妇av| 少妇人妻精品综合一区二区| 久久久国产精品麻豆| 久久精品久久久久久噜噜老黄| av国产久精品久网站免费入址| tube8黄色片| 精品福利永久在线观看| 免费大片18禁| 老女人水多毛片| 日本黄大片高清| 美女主播在线视频| 久久鲁丝午夜福利片| 人妻系列 视频| 久久久久精品性色| 熟女人妻精品中文字幕| 久久久久精品性色| 建设人人有责人人尽责人人享有的| 国产在线一区二区三区精| 久久久国产欧美日韩av| 一级爰片在线观看| 国产成人一区二区在线| 老司机影院毛片| 国产精品国产三级专区第一集| 亚洲欧洲国产日韩| 国产精品欧美亚洲77777| xxxhd国产人妻xxx| 国产精品久久久久久av不卡| a 毛片基地| 国产亚洲av片在线观看秒播厂| 人人妻人人澡人人爽人人夜夜| 日韩av在线免费看完整版不卡| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜爱| 国产亚洲精品久久久com| 亚洲五月色婷婷综合| 自线自在国产av| 亚洲人成网站在线观看播放| 人体艺术视频欧美日本| 成人亚洲精品一区在线观看| 欧美激情极品国产一区二区三区 | 国产一区二区在线观看日韩| 国产一区二区激情短视频 | videossex国产| 亚洲欧美一区二区三区国产| 亚洲图色成人| 日韩欧美一区视频在线观看| 成人亚洲欧美一区二区av| 母亲3免费完整高清在线观看 | 久久久亚洲精品成人影院| 久久久久久久精品精品| 成年动漫av网址| 久久99蜜桃精品久久| 亚洲成人手机| kizo精华| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 国产精品一区二区在线观看99| xxx大片免费视频| 国产亚洲一区二区精品| 日本-黄色视频高清免费观看| 国产片内射在线| 国产精品久久久久久久久免| 男的添女的下面高潮视频| 日本爱情动作片www.在线观看| 久久久国产欧美日韩av| 婷婷色综合大香蕉| 97人妻天天添夜夜摸| 亚洲成av片中文字幕在线观看 | √禁漫天堂资源中文www| 超碰97精品在线观看| 黄色一级大片看看| av一本久久久久| 韩国精品一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 久久久亚洲精品成人影院| 啦啦啦视频在线资源免费观看| 少妇精品久久久久久久| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 精品国产一区二区三区久久久樱花| 欧美国产精品一级二级三级| 91成人精品电影| 亚洲精品久久久久久婷婷小说| 在线观看免费高清a一片| 成年女人在线观看亚洲视频| av片东京热男人的天堂| 一二三四在线观看免费中文在 | 天堂8中文在线网| av黄色大香蕉| 欧美日韩视频精品一区| 中文乱码字字幕精品一区二区三区| 日本黄色日本黄色录像| 亚洲国产欧美在线一区| 自线自在国产av| 日韩一本色道免费dvd| 卡戴珊不雅视频在线播放| 精品人妻在线不人妻| 视频区图区小说| 久久久久网色| 亚洲av电影在线进入| 亚洲欧洲日产国产| 涩涩av久久男人的天堂| 免费在线观看黄色视频的| 国产精品99久久99久久久不卡 | 深夜精品福利| 多毛熟女@视频| 成年av动漫网址| 满18在线观看网站| 免费观看a级毛片全部| 一本—道久久a久久精品蜜桃钙片| 日韩三级伦理在线观看| 国产黄色视频一区二区在线观看| 精品亚洲成国产av| 王馨瑶露胸无遮挡在线观看| 一级黄片播放器| 巨乳人妻的诱惑在线观看| 国产免费福利视频在线观看| videossex国产| 国产伦理片在线播放av一区| 亚洲国产欧美在线一区| 麻豆精品久久久久久蜜桃| 精品国产露脸久久av麻豆| 亚洲国产精品一区三区| 久久久精品免费免费高清| 日本vs欧美在线观看视频| 天堂俺去俺来也www色官网| 国产熟女欧美一区二区| 午夜福利视频精品| 久久久欧美国产精品| 成人国产av品久久久| 亚洲国产精品专区欧美| 在线观看免费日韩欧美大片| 国产色爽女视频免费观看| 丁香六月天网| 欧美人与善性xxx| 久久久久久久久久久久大奶| 日本vs欧美在线观看视频| 美女福利国产在线| 国产淫语在线视频| av又黄又爽大尺度在线免费看| 午夜精品国产一区二区电影| 国产精品成人在线| 韩国av在线不卡| 亚洲第一区二区三区不卡| 亚洲欧洲精品一区二区精品久久久 | 日韩精品有码人妻一区| 人妻人人澡人人爽人人| 夜夜爽夜夜爽视频| 伊人亚洲综合成人网| 日本猛色少妇xxxxx猛交久久| 性色avwww在线观看| 久久狼人影院| 日韩av在线免费看完整版不卡| 两个人看的免费小视频| 18禁动态无遮挡网站| 国产高清不卡午夜福利| 国产男人的电影天堂91| 一级黄片播放器| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 高清欧美精品videossex| 亚洲丝袜综合中文字幕| 亚洲av中文av极速乱| 日本av手机在线免费观看| 另类精品久久| 日本91视频免费播放| 久久精品国产亚洲av涩爱| 欧美xxⅹ黑人| 香蕉国产在线看| 国产激情久久老熟女| 国产成人免费无遮挡视频| 精品人妻一区二区三区麻豆| av网站免费在线观看视频| 中文字幕最新亚洲高清| 久久久久久人人人人人| 99国产综合亚洲精品| 亚洲欧美成人综合另类久久久| 观看美女的网站| 啦啦啦视频在线资源免费观看| 国产精品久久久久久精品古装| 国产成人精品福利久久| 爱豆传媒免费全集在线观看| 在线免费观看不下载黄p国产| 免费看光身美女| 七月丁香在线播放| 日韩 亚洲 欧美在线| 国产成人一区二区在线| 最近最新中文字幕大全免费视频 | 亚洲精品日韩在线中文字幕| 69精品国产乱码久久久| 久久久久人妻精品一区果冻| 在线免费观看不下载黄p国产| 蜜桃在线观看..| 亚洲精品,欧美精品| 男男h啪啪无遮挡| 国产又爽黄色视频| 男女啪啪激烈高潮av片| 免费观看性生交大片5| 丁香六月天网| 亚洲av中文av极速乱| 嫩草影院入口| 精品国产一区二区三区四区第35| 大香蕉久久成人网| 国产免费一区二区三区四区乱码| 欧美97在线视频| 国产精品一国产av| 国产亚洲精品久久久com| 亚洲成av片中文字幕在线观看 | 午夜免费观看性视频| 伊人亚洲综合成人网| 亚洲国产欧美在线一区| 国内精品宾馆在线| 九九爱精品视频在线观看| 久久精品人人爽人人爽视色| 18在线观看网站| 女人精品久久久久毛片| 我的女老师完整版在线观看| 国产免费视频播放在线视频| 国产高清三级在线| 天堂中文最新版在线下载| 日韩不卡一区二区三区视频在线| 久久久久人妻精品一区果冻| 免费日韩欧美在线观看| 国产亚洲精品第一综合不卡 | 亚洲av电影在线观看一区二区三区| 热re99久久精品国产66热6| 午夜免费观看性视频| 亚洲一区二区三区欧美精品| 日本猛色少妇xxxxx猛交久久| 欧美人与性动交α欧美软件 | 免费在线观看黄色视频的| 五月玫瑰六月丁香| 最后的刺客免费高清国语| 国产乱人偷精品视频| 亚洲欧美精品自产自拍| 欧美国产精品一级二级三级| 国产欧美日韩一区二区三区在线| 国产成人aa在线观看| 国产探花极品一区二区| 亚洲内射少妇av| 国产无遮挡羞羞视频在线观看| 国产熟女欧美一区二区| 哪个播放器可以免费观看大片| 久久97久久精品| 亚洲少妇的诱惑av| 街头女战士在线观看网站| 久久国内精品自在自线图片| 欧美人与善性xxx| 又粗又硬又长又爽又黄的视频| 欧美 日韩 精品 国产| 日本黄大片高清| 观看美女的网站| 亚洲图色成人| 中文字幕另类日韩欧美亚洲嫩草| 男女高潮啪啪啪动态图| 香蕉精品网在线| 插逼视频在线观看| 国产视频首页在线观看| 亚洲图色成人| 22中文网久久字幕| 欧美精品亚洲一区二区| 国产永久视频网站| 亚洲人与动物交配视频| 91午夜精品亚洲一区二区三区| 在线 av 中文字幕| 在线看a的网站| 春色校园在线视频观看| 日本免费在线观看一区| 春色校园在线视频观看| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 久久久久久久久久人人人人人人| 99久久综合免费| 大香蕉久久网| 亚洲精品av麻豆狂野| 国产成人一区二区在线| 在线观看美女被高潮喷水网站| 亚洲精品一二三| 97在线视频观看| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 在线观看人妻少妇| 爱豆传媒免费全集在线观看| 午夜影院在线不卡| 国产成人精品婷婷| 成年女人在线观看亚洲视频| 国产亚洲午夜精品一区二区久久| 久久久精品免费免费高清| 国产激情久久老熟女| 精品一区二区三卡| 国产男女内射视频| 成人二区视频| 免费av不卡在线播放| 国产麻豆69| 啦啦啦中文免费视频观看日本| 丰满饥渴人妻一区二区三| 欧美老熟妇乱子伦牲交| 汤姆久久久久久久影院中文字幕| 国产精品久久久久久av不卡| 美女福利国产在线| 亚洲av综合色区一区| 九九在线视频观看精品| 欧美另类一区| 欧美成人午夜免费资源| 国产黄色免费在线视频| 毛片一级片免费看久久久久| 亚洲人成网站在线观看播放| av免费在线看不卡| 天天躁夜夜躁狠狠躁躁| 男女啪啪激烈高潮av片| 丰满乱子伦码专区| 久久青草综合色| 在线观看www视频免费| 亚洲国产精品999| 极品人妻少妇av视频| 欧美日韩av久久| 波野结衣二区三区在线| 一级毛片电影观看| 黑人巨大精品欧美一区二区蜜桃 | 国产在视频线精品| 亚洲国产精品一区二区三区在线| 五月天丁香电影| 满18在线观看网站| 九色亚洲精品在线播放| 亚洲av国产av综合av卡| xxxhd国产人妻xxx| 男人舔女人的私密视频| xxxhd国产人妻xxx| 三级国产精品片| 免费看光身美女| 亚洲综合色网址| 黄网站色视频无遮挡免费观看| 欧美激情国产日韩精品一区| 免费高清在线观看视频在线观看| 精品少妇久久久久久888优播| 国产日韩欧美在线精品| 成人亚洲欧美一区二区av| 国产熟女午夜一区二区三区| 午夜福利影视在线免费观看| 亚洲精品久久成人aⅴ小说| 国产永久视频网站| 国国产精品蜜臀av免费| 久久人人97超碰香蕉20202| 搡女人真爽免费视频火全软件| 国产福利在线免费观看视频| 一级片'在线观看视频| 中文字幕制服av| freevideosex欧美| 又黄又粗又硬又大视频| 久久精品国产亚洲av涩爱| 午夜久久久在线观看| 在线观看人妻少妇| 视频区图区小说| 国产精品久久久久久久久免| 在线观看免费日韩欧美大片| 久久久久国产精品人妻一区二区| 在线看a的网站| 中文字幕精品免费在线观看视频 | 一级a做视频免费观看| 日本91视频免费播放| 成人国产av品久久久| 最近2019中文字幕mv第一页| 婷婷色综合大香蕉| 亚洲欧美色中文字幕在线| 五月开心婷婷网| 2021少妇久久久久久久久久久| 婷婷色麻豆天堂久久| 狠狠精品人妻久久久久久综合| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 建设人人有责人人尽责人人享有的| 国产国语露脸激情在线看| 国产又爽黄色视频| 两个人看的免费小视频| 精品国产一区二区三区久久久樱花| 女性被躁到高潮视频| 国产成人精品久久久久久| 国产亚洲精品久久久com| 国产熟女欧美一区二区| 中文字幕av电影在线播放| 考比视频在线观看| 一本色道久久久久久精品综合| 成年av动漫网址| 国产精品一国产av| 久久人人爽人人片av| freevideosex欧美| 久久99热6这里只有精品| 这个男人来自地球电影免费观看 | 久久国产精品大桥未久av| 国精品久久久久久国模美| 欧美人与善性xxx| 欧美性感艳星| 老司机影院毛片| 日韩 亚洲 欧美在线| 日韩伦理黄色片| 美女中出高潮动态图| xxxhd国产人妻xxx| 久久99一区二区三区| 婷婷成人精品国产| 大码成人一级视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品aⅴ在线观看| 亚洲av综合色区一区| 国产精品国产三级国产av玫瑰| 国产日韩欧美在线精品| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 国产乱来视频区| 久久久久久伊人网av| kizo精华| 亚洲精品国产av成人精品| 亚洲精品,欧美精品| 如何舔出高潮| 久久这里只有精品19| 亚洲少妇的诱惑av| 亚洲天堂av无毛| 在线免费观看不下载黄p国产| 精品福利永久在线观看| 久久久久网色| 中文字幕人妻丝袜制服| av电影中文网址| 亚洲第一区二区三区不卡| 久久精品国产综合久久久 | 午夜免费鲁丝| 国产熟女欧美一区二区| 激情视频va一区二区三区| 18禁裸乳无遮挡动漫免费视频| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 亚洲在久久综合| 国产精品女同一区二区软件| 久久久欧美国产精品| 欧美性感艳星| 免费大片黄手机在线观看| 人人妻人人添人人爽欧美一区卜| 欧美+日韩+精品| 丰满乱子伦码专区| 国产成人一区二区在线| 成人18禁高潮啪啪吃奶动态图| av视频免费观看在线观看| 欧美国产精品一级二级三级| 精品一区二区三卡| 国产精品三级大全| 亚洲丝袜综合中文字幕| 一级a做视频免费观看| 国产av一区二区精品久久| 少妇猛男粗大的猛烈进出视频| 国产永久视频网站| 亚洲精品一区蜜桃| 久久热在线av| 亚洲av日韩在线播放| 纯流量卡能插随身wifi吗| 国产男人的电影天堂91| 国产精品.久久久| 99热6这里只有精品| 色吧在线观看| 国产在线视频一区二区| 欧美日韩av久久| 极品人妻少妇av视频| 国产毛片在线视频| 国产老妇伦熟女老妇高清| 免费播放大片免费观看视频在线观看| 国产av一区二区精品久久| 午夜91福利影院| 欧美精品一区二区免费开放| 美女中出高潮动态图| 大片免费播放器 马上看| 中文天堂在线官网| 久久精品久久久久久久性| 国产国语露脸激情在线看| 老熟女久久久| av片东京热男人的天堂| 欧美97在线视频| 久久综合国产亚洲精品| 黄色一级大片看看| 9191精品国产免费久久| 国产精品.久久久| 精品国产国语对白av| 国产av精品麻豆| 日韩欧美精品免费久久| 国产男女超爽视频在线观看| 99久久精品国产国产毛片| 狠狠婷婷综合久久久久久88av| 在现免费观看毛片| 视频在线观看一区二区三区| 999精品在线视频| 日韩不卡一区二区三区视频在线| 亚洲精品国产av蜜桃| 菩萨蛮人人尽说江南好唐韦庄| av在线app专区| av又黄又爽大尺度在线免费看| 18禁国产床啪视频网站| 久久午夜福利片| 亚洲国产精品一区二区三区在线| 国产成人aa在线观看|