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

    渭河流域多尺度干旱時空特征分析

    2017-07-19 10:03:54劉蕊娟王旭紅趙會娟
    關(guān)鍵詞:渭河流域天水中度

    劉蕊娟,王旭紅,趙會娟

    (1.西北大學(xué)城市與環(huán)境學(xué)院, 陜西 西安 710127;2.北京大學(xué)深圳研究生院,城市人居環(huán)境科學(xué)與技術(shù)重點實驗室, 廣東 深圳 518055)

    渭河流域多尺度干旱時空特征分析

    劉蕊娟1,王旭紅1,趙會娟2

    (1.西北大學(xué)城市與環(huán)境學(xué)院, 陜西 西安 710127;2.北京大學(xué)深圳研究生院,城市人居環(huán)境科學(xué)與技術(shù)重點實驗室, 廣東 深圳 518055)

    基于渭河流域23個氣象站1959—2014年的實測資料,對比分析了標(biāo)準(zhǔn)化降水指數(shù)(SPI)和標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)在渭河流域的適用性,得出SPEI的穩(wěn)定性更佳。繼而采用標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)研究該流域干旱發(fā)生頻率與干旱強(qiáng)度,分析了近55年渭河流域不同時間尺度干旱的時空變化特征。結(jié)果表明:尺度越大,SPEI波動的幅度越小。渭河流域SPEI在1991年發(fā)生突變,1991年后干旱加劇,干旱的發(fā)生存在14 a和7 a左右周期變化。年尺度干旱的發(fā)生頻率呈東高西低的分布,中度以上干旱頻率最高的地區(qū)為西吉、武功和天水周邊,頻率分別為14.28%、14.28%、12.5%;季尺度上整體的秋旱相比其它季節(jié)的干旱發(fā)生頻率高,中度以上春、秋、冬三季干旱發(fā)生的頻率天水盆地均最高,頻率分別為21.4%、19.6%和8.9%,中度以上夏季干旱發(fā)生頻率最高的是平?jīng)?,其頻率為16.1%;月尺度上干旱強(qiáng)度的分布總體上是由東向西減小,其中干旱強(qiáng)度最大的是南部佛坪地區(qū),干旱強(qiáng)度為30%,天水和海源地區(qū)近55年發(fā)生連續(xù)干旱過程的次數(shù)最多,可達(dá)22 a,該過程中連續(xù)干旱的月數(shù)大多為3~4個月,而吳旗站存在一次長達(dá)9個月的持續(xù)干旱過程。

    標(biāo)準(zhǔn)化降水指數(shù)(SPI);標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI);多尺度干旱;時空變化特征;渭河流域

    全球氣候變化是人類迄今為止面臨的最主要的全球環(huán)境問題,這一挑戰(zhàn)引起了世界各國的廣泛關(guān)注[1]。20世紀(jì)70年代以來,全球氣候變暖導(dǎo)致的干旱等極端氣候事件的頻率和強(qiáng)度都呈顯著增加趨勢[2]。全球氣象災(zāi)害造成的損失經(jīng)計算約占各類自然災(zāi)害損失的85%,其中干旱造成的損失占?xì)庀鬄?zāi)害的50%以上[3]。由此,干旱已經(jīng)成為各國政府和學(xué)者關(guān)注的焦點問題。

    國際上通常將干旱劃分為氣象、農(nóng)業(yè)、水文和社會經(jīng)濟(jì)干旱,四類干旱的本質(zhì)均是由于降水不足而導(dǎo)致的水分虧缺現(xiàn)象[3]。氣象干旱表現(xiàn)在月尺度上降水量的減少;農(nóng)業(yè)干旱是以季節(jié)為尺度來衡量土壤中的水分虧缺,繼而表現(xiàn)為農(nóng)作物減產(chǎn)或絕收,糧食供給失衡等;水文干旱以半年為尺度,表現(xiàn)為河流水庫的水資源減少、地下水位下降等;以上三類干旱在一定程度上則會引發(fā)社會經(jīng)濟(jì)干旱[4-5]。國內(nèi)外眾多學(xué)者對于干旱的研究逐步深入,從歷史記載到定量分析、從數(shù)值統(tǒng)計到遙感監(jiān)測。針對多時間尺度的干旱評價眾多學(xué)者也提出了多種干旱監(jiān)測指標(biāo),如Palmer干旱指數(shù)(PDSI)、Z指數(shù)、作物缺水指數(shù)(CWSI)、標(biāo)準(zhǔn)化降水量指數(shù)(SPI)、標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)等[6-10]。通過對不同指數(shù)在不同區(qū)域的應(yīng)用分析發(fā)現(xiàn)標(biāo)準(zhǔn)化蒸散指數(shù)(SPEI)在監(jiān)測年尺度水文干旱和季尺度農(nóng)業(yè)干旱方面的能力要優(yōu)于其它干旱指標(biāo)[11-13]。目前SPEI在中國已有一些應(yīng)用,莊少偉等[12]研究了SPEI在中國區(qū)域的應(yīng)用,表明在我國年均降水量大于200 mm的區(qū)域各種尺度SPEI均適用;李偉光等[13-14]分析了近60年中國和華南地區(qū)的干旱時空特征,得出的結(jié)果與其它干旱指標(biāo)一致;任培貴、馬瓊[15-16]等運(yùn)用SPEI指標(biāo)分析了西北地區(qū)和黃土高原地區(qū)的干旱特征,得出的結(jié)論與年鑒資料一致,說明SPEI在西北地區(qū)也具有較好的適用性。

    渭河是黃河最大的支流,該流域大部分位于陜西省境內(nèi),集中了陜西省61%的人口、56%的耕地,72%的灌溉面積、68%的糧食產(chǎn)量和80%以上的國內(nèi)生產(chǎn)總值,在陜西乃至西部經(jīng)濟(jì)社會發(fā)展中都具有十分重要的戰(zhàn)略地位[17-18]。但隨著經(jīng)濟(jì)社會的快速發(fā)展,渭河流域水資源短缺,水質(zhì)污染和水土流失嚴(yán)重[18]。渭河流域生態(tài)環(huán)境的惡化與旱澇災(zāi)害的頻發(fā)成為眾多學(xué)者的關(guān)注點。諸多學(xué)者依據(jù)各種干旱評價指標(biāo)研究了渭河流域的干旱演變[19-24],但對于多尺度干旱的空間分布研究還有待深入,因此本文運(yùn)用不同尺度的SPI和SPEI對比分析其適用性,最終選取SPEI研究該流域的氣象、農(nóng)業(yè)和水文干旱演變特征,以期為該流域防旱抗旱提供指導(dǎo)。

    1 研究區(qū)概況

    渭河流域位于中國西北部(圖1),是黃河的第一大支流,流經(jīng)甘肅、陜西和寧夏三省。自然區(qū)劃上其北部為黃土高原,南部為秦嶺山脈,流經(jīng)隴東黃土高原、天水盆地、寶雞峽谷進(jìn)入關(guān)中平原,流域面積134 800 km2。該流域?qū)儆诖箨懶詺夂颍彩堑湫偷母珊蛋敫珊祬^(qū),地處我國西北地區(qū)東部生態(tài)脆弱帶[17,25],自然災(zāi)害多發(fā),其中旱災(zāi)發(fā)生頻繁,且程度最重,危害最大,成為制約該地區(qū)國民經(jīng)濟(jì),尤其是農(nóng)業(yè)經(jīng)濟(jì)和農(nóng)業(yè)可持續(xù)發(fā)展的主要因素。

    圖1 渭河流域位置與氣象站點分布

    2 數(shù)據(jù)與方法

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

    本文選取了了渭河流域及周邊地區(qū)共23個氣象站點1959—2014年逐月的降水和氣溫資料,主要來自中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.gov.cn/),天水站數(shù)據(jù)來自甘肅省氣象局,西安站、寶雞站2005和2006年缺測數(shù)據(jù)均來自陜西省統(tǒng)計年鑒。主要利用均值替代方法對吳旗站和洛川站個別缺測數(shù)據(jù)進(jìn)行了插補(bǔ)處理,經(jīng)過訂正后的各臺站氣象資料連續(xù)。本文的研究以年、季、月為尺度,季節(jié)的劃分:3—5月為春季、6—8月為夏季、9—11月為秋季、12月至翌年2月為冬季。

    2.2 研究方法

    2.2.1 標(biāo)準(zhǔn)化降水指數(shù)(SPI) 標(biāo)準(zhǔn)化降水指數(shù)(SPI)認(rèn)為降水量服從Γ分布,通過Γ分布概率密度函數(shù)求累計概率,再將累積的概率標(biāo)準(zhǔn)化。計算過程中因為不涉及不同地域的空間參數(shù),因此在多區(qū)域適用。具體計算步驟和劃分等級可參見文獻(xiàn)[9,26-27],SPI的干旱等級劃分如表1所示。

    2.2.2 標(biāo)準(zhǔn)化蒸散指數(shù)(SPEI) 標(biāo)準(zhǔn)化蒸散指數(shù)(SPEI)是以月降水量、月均溫及氣象站點的緯度作為輸入值,采用Vicente-Serrano所述的Thornthwaite方法計算潛在蒸散發(fā)(PET),繼而計算月降水量與潛在蒸散發(fā)量的差值,再通過差值的分布概率密度函數(shù)(log-logistic函數(shù))求累計概率,最后轉(zhuǎn)化為標(biāo)準(zhǔn)正態(tài)分布 來表征該區(qū)域不同時間尺度的旱澇狀況,其詳細(xì)步驟可參見文獻(xiàn)[12-16]。根據(jù)《氣象干旱等級》規(guī)范[28],將SPEI的干旱等級劃分如表1所示。

    表1 SPI和SPEI干旱等級劃分與發(fā)生累積概率

    2.2.3 干旱評估方法 本文選用旱澇發(fā)生頻率和強(qiáng)度指標(biāo)來評估干旱[27,29]

    (1) 干旱頻率Pi:Pi可按不同程度干旱發(fā)生的年數(shù)計算相應(yīng)程序干旱頻率,計算公式為:

    Pi=(n/N)×100%

    (1)

    式中,n為某站點發(fā)生干旱的年數(shù);N為降水資料年數(shù);i為某站點。

    (2) 干旱強(qiáng)度:當(dāng)SPEI 1的數(shù)值連續(xù)3個月在輕度干旱以上,則發(fā)生一次連續(xù)干旱過程。干旱過程中各月的SPEI 1平均值為這次干旱過程的SPEI值,然后再根據(jù)各站點所有干旱過程累計概率求取平均值,則得到干旱過程的強(qiáng)度。平均值越小,干旱強(qiáng)度越強(qiáng)。本文的干旱等級和累計概率按照表1的標(biāo)準(zhǔn)來劃分[16,29-30]。

    3 結(jié)果與分析

    3.1 SPI和SPEI的對比驗證

    根據(jù)《陜西省救災(zāi)年鑒》和陜西省地情網(wǎng)對旱情的統(tǒng)計資料,驗證SPI和SPEI在渭河流域的適用性。根據(jù)記載,陜西省渭北和渭河流域關(guān)中地區(qū)歷史上的干旱災(zāi)害頻發(fā),據(jù)統(tǒng)計40年中有23次旱災(zāi),并且將關(guān)中中部劃分為重伏旱區(qū)。1959—1962年,西安、渭南、咸陽和安康四地旱象嚴(yán)重,整個陜西省受災(zāi)面積達(dá)145.58 hm2。1994—1995年關(guān)中24個縣(市)降水較歷年偏少5—7成。

    因此本文以渭河流域關(guān)中段的西安站為代表站,根據(jù)SPI和SPEI的月和年尺度的數(shù)值來表現(xiàn)干旱等級。由圖2可知SPI和SPEI在月尺度(SPI 1/SPEI 1)與年尺度(SPI 12/SPEI 12)上變化有所不同。月尺度上,雖然SPI 1和SPEI 1的波動都相對較強(qiáng),但在特殊年份的典型月份,如1987年9月和1994年5月以及2002年7月SPI 1的值都小于SPEI,SPI 1所反映的干旱等級高于SPEI 1,雖然兩者都考慮了累積效應(yīng),但是SPI 1的個別值和前后月份的值差距較大,波動明顯,說明其只考慮降水而忽略溫度和蒸散所得出的結(jié)果對干旱的反應(yīng)過強(qiáng)。年尺度上SPI 12和SPEI 12兩者具有很好的一致性,但是以2013年SPI 12和SPEI 12的值為例,統(tǒng)計后年平均的SPEI 12反映的干旱級別要重于SPI 12,2013年平均的SPEI 12反映的干旱等級為中度干旱,而SPI為正常等級。根據(jù)西安市氣象臺2013年3月發(fā)布的干旱黃色預(yù)警記載,西安市農(nóng)田土壤相對濕度達(dá)到中度干旱等級,4月旱情將持續(xù)發(fā)展。而圖2統(tǒng)計的結(jié)果SPEI 12在3月和4月為中度干旱等級,SPI 12在3月和4月為輕度干旱等級。因此,SPEI 1/SPEI 12對干旱的評估更加準(zhǔn)確,并且對于不同尺度的響應(yīng)也更加符合實際。

    圖2 渭河流域西安站1959—2013年SPI和SPEI月尺度(SPI 1/SPEI 1)和年尺度(SPI 12/SPEI 12)比較

    3.2 基于SPEI的渭河流域干旱時間特征

    3.2.1 SPEI表征的不同時間尺度特征 渭河流域不同時間尺度的SPEI波動規(guī)律明顯不同(圖3),月尺度(SPEI 1)的波動最強(qiáng),隨著時間尺度的增大,季尺度(SPEI 3)和年尺度(SPEI 12)的波動起伏減緩。SPEI 1由于每月降水和溫度的變化,能更準(zhǔn)確地反映氣象干旱和濕潤情況;SPEI 3受到季節(jié)降水和溫度的變化,能反映土壤中的水分虧缺;SPEI 12則反映出長期的水文旱情。根據(jù)年尺度(SPEI 12)每年的旱澇指數(shù)可以得出2000年的SPEI 12的值最小,呈輕度干旱等級。以2000年的各尺度SPEI 1/SPEI 3/SPEI 12值為例,通過統(tǒng)計SPEI 3可以得出春季旱情最強(qiáng),其值為-1.329,屬于中度干旱。月尺度SPEI 1則可以看出2000年的7月份的旱情最嚴(yán)重,SPEI 1值為-1.648,這和SPEI 12各月份的值中7月份的SPEI 12值最小(-1.834)相吻合。說明不同尺度SPEI反映不同類型的干旱。對于尺度的選擇要根據(jù)干旱類型的研究來定奪。

    圖3 1959—2014年渭河流域不同時間尺度歷年SPEI值

    3.2.2 年尺度干旱突變檢驗 對渭河流域1959—2014年的SPEI 12進(jìn)行Mann-Kendall突變檢驗,結(jié)果如圖4所示,UF值在1972年之前都是大于0的,表明該SPEI序列呈上升趨勢;UF值在1991年之后都是小于0的,表明該SPEI序列呈下降趨勢,尤其是1997年之后UF曲線超出了置信區(qū)間,說明其下降趨勢顯著。在置信范圍內(nèi)UF曲線和UB曲線相交于1991年,說明其突變始于1991年;突變前SPEI的均值為0.286,突變后SPEI的均值為-0.373,間接說明自1991年后渭河流域的干旱加劇。

    圖4 渭河流域1959—2014年平均SPEI指數(shù)的M-K突變檢驗

    3.2.3 年尺度干旱周期特征 為了進(jìn)一步研究渭河流域干旱的周期特征,本文對渭河流域1959—2014年的SPEI值進(jìn)行了小波分析,選用了水文氣象常用的Morlet連續(xù)復(fù)小波變換,如圖5(左)的小波實部等值線圖,藍(lán)色表示負(fù)值,紅色表示正值,在流域旱澇演變過程中,明顯存在著6—8 a和10—16 a的周期,在10—16 a的時間尺度上出現(xiàn)6次明顯的旱澇交替;6—8 a的尺度上有明顯的3次旱澇交替。又結(jié)合圖5(右)小波方差圖得知在7 a、11 a和14 a時間尺度上震蕩明顯,最終表明14 a、11 a和7 a分別為第一主周期和第二、三周期。根據(jù)小波實部圖可以看出2014年虛線未完全閉合,值為正,屬于濕潤區(qū)間,說明目前正處在相對濕潤的時期,接下來的一個周期內(nèi),則向負(fù)值過渡,是變干旱的趨勢。

    3.3 渭河流域干旱的空間分布特征

    3.3.1 年、季尺度干旱發(fā)生頻率 根據(jù)干旱頻率的計算得出在渭河流域近60年各地區(qū)干旱發(fā)生的頻率(圖6)有所不同。在渭河上游甘肅段天水和平?jīng)龅貐^(qū)干旱頻率較低,自西南向東北頻率增強(qiáng)。支流北洛河上游的吳旗和延安周邊地區(qū)干旱發(fā)生頻率較高,分別為30.4%和32.1%。渭河流域中游關(guān)中段的長武年干旱發(fā)生頻率較高為30.4%,銅川和西安附近干旱頻率相較于長武較低,年干旱的發(fā)生頻率分別為27.3%和29.1%,關(guān)中地區(qū)周邊的華山和佛坪是整個流域年干旱發(fā)生頻率最高的地區(qū),在近56年中佛坪和華山附近分別發(fā)生了21個和19個年頭的干旱,干旱頻率分別可達(dá)37.5%和33.9%。但是總體而言,渭河流域干流關(guān)中段比其它地區(qū)發(fā)生干旱的頻率要高,上游天水盆地附近干旱發(fā)生頻率最低。其支流北洛河和涇河區(qū)域的年干旱發(fā)生頻率要高于干流。

    春旱(圖7A)主要發(fā)生在陜西關(guān)中地區(qū)的西部和中部,并且以關(guān)中地區(qū)為核心呈環(huán)狀的頻率高-低-高的局勢,其中隴東黃土高原的春旱相比陜西黃土高原的干旱發(fā)生頻率較低。春季干旱頻率最大值出現(xiàn)在寶雞附近,達(dá)33.9%,隴東黃土高原西峰和平?jīng)龈浇拇汉殿l率約為23%~26%,天水地區(qū)發(fā)生春旱的頻率最低,約為23%。夏季發(fā)生干旱的頻率(圖7B)在整個流域較春季高,并且大部分地區(qū)夏旱發(fā)生頻率顯著,尤以渭河上游岷縣附近和天水地區(qū)最高,頻率可達(dá)30%~32.1%;除此以外在關(guān)中平原北部的長武、關(guān)中西部的華山和黃土高原的環(huán)縣附近頻率較高,其頻率可達(dá)30.4%;寧夏南部和甘肅華家?guī)X附近地區(qū)夏季干旱發(fā)生的頻率最低。秋旱的發(fā)生的頻率(圖7C)在渭河流域并不高,只有在渭河流域東北部延安地區(qū)頻率能達(dá)到35%;在渭河流域的西北部,甘肅華家?guī)X附近地區(qū)秋旱的頻率也相對較高,可以達(dá)到33.9%;其余大部分地區(qū)秋旱發(fā)生的頻率都分布在27%~31%,關(guān)中地區(qū)秋旱的頻率相對其它地區(qū)較低。冬旱發(fā)生的頻率(圖7D)尤以關(guān)中地區(qū)最為嚴(yán)重,武功地區(qū)的冬旱發(fā)生頻率約為25%~27%,陜西黃土高原的吳旗附近頻率最低為17.9%,其它區(qū)域冬旱發(fā)生的頻率為21%~23%。

    3.3.2 中度以上干旱發(fā)生頻率 年尺度方面中度以上干旱發(fā)生頻率(圖8)較高的地區(qū)分布在渭河流域西部和南部,尤其以寧夏南部和關(guān)中地區(qū)中部的中度以上干旱發(fā)生頻率最高,其中西吉和武功周邊的年尺度中度以上干旱頻率最高可達(dá)14.28%;關(guān)中地區(qū)北部即黃土高原南部地區(qū)年尺度中度以上干旱反而較低,長武和銅川周邊干旱的發(fā)生頻率低至3.5%~5.5%,整個黃土高原地區(qū)除了隴東黃土高原的西峰以外地區(qū),年尺度的干旱發(fā)生頻率都相對較低,大多介于9%~10%,西峰地區(qū)干旱的發(fā)生頻率可達(dá)12.5%。這與總體各個等級干旱的發(fā)生頻率相較而言存在極大的差異,總的4個干旱等級上渭河流域上游的天水盆地是干旱發(fā)生頻率較低的地區(qū)(圖6),但是就中度以上干旱的發(fā)生頻率而言,天水盆地附近卻是中度和極端干旱發(fā)生的典型區(qū)域,說明天水地區(qū)的干旱雖發(fā)生頻率低,但是強(qiáng)度大,等級高,旱情嚴(yán)重。

    春季渭河流域中度以上干旱發(fā)生頻率(圖9a)9%~14.2%,渭河流域上游天水地區(qū)春季中度以上干旱的發(fā)生頻率較高,整個時間段內(nèi)天水盆地周邊中度和極端干旱發(fā)生的頻率高達(dá)21.4%,成為全域內(nèi)春季干旱級別最高、發(fā)生頻率最強(qiáng)的地區(qū)。除此以外整個渭河流域的西峰和洛川附近屬于春季中度以上干旱發(fā)生頻率較高的地區(qū)。整體而言,渭河流域東南部比西部和北部春季中度以上干旱的發(fā)生頻率略微要高。這與整個干旱等級春季干旱發(fā)生頻率(圖7A)相比較也有差異,最明顯的是天水和寶雞地區(qū),天水盆地各個等級干旱的發(fā)生頻率為23.2%,但是僅春季中度以上干旱可達(dá)21.4%,說明天水盆地附近有嚴(yán)重的春旱;關(guān)中地區(qū)西部寶雞附近各個等級干旱的發(fā)生頻率為33.9%,但中度以上的干旱僅為14.5%,說明關(guān)中地區(qū)西部春季雖然有干旱,但是干旱等級相對較低。

    圖5 渭河流域1959—2014年平均SPEI指數(shù)小波實部圖(左)和小波方差(右)

    圖6 渭河流域年尺度干旱發(fā)生頻率分布

    夏季整個渭河流域中度以上干旱的發(fā)生頻率(圖9b)存在明顯的東西差異,西部隴東黃土高原的平?jīng)龊臀骷貐^(qū)中度以上干旱發(fā)生頻率可達(dá)16.1%,黃土高原中部吳旗地區(qū)和關(guān)中地區(qū)西部的寶雞中度以上干旱發(fā)生頻率為14.3%,其它地區(qū)頻率都較低。相比較總的4個等級夏季干旱發(fā)生的頻率(圖7B),長武、天水、西吉和華山存在明顯的差異,長武、華山和天水周邊所有等級的干旱發(fā)生頻率較高,但是中度和極端干旱的發(fā)生頻率低,說明發(fā)生夏旱的頻率高,但等級低。隴東黃土高原的平?jīng)龅貐^(qū)所有等級干旱的發(fā)生頻率達(dá)25%,相對其它地區(qū)較低,但是夏季中度和極端干旱發(fā)生頻率16.1%,是整個流域地區(qū)最高的,說明該地區(qū)發(fā)生夏旱的頻率低,但是等級高,夏旱嚴(yán)重。

    圖7 渭河流域季尺度干旱發(fā)生頻率分布

    圖8 渭河流域年尺度中度以上干旱發(fā)生頻率分布

    圖9 渭河流域年季尺度中度以上干旱發(fā)生頻率分布

    圖10 干旱發(fā)生強(qiáng)度分布

    秋季整個渭河流域中度以上干旱的發(fā)生頻率(圖9c)最高的區(qū)域分布于渭河流域甘肅段,其中天水地區(qū)中度以上干旱的發(fā)生頻率最高,可達(dá)19.6%;流域內(nèi)其它地區(qū)中度以上干旱的發(fā)生頻率為11%~16%,無明顯的差異,相較于各個干旱等級發(fā)生頻率(圖7C),天水盆地附近是屬于干旱發(fā)生等級高但頻率較低的區(qū)域,陜西境內(nèi)黃土高原東部延安地區(qū)雖然不同等級干旱的發(fā)生頻率相對較高,但是中度以上干旱的發(fā)生頻率低,說明該區(qū)域?qū)儆诟珊蛋l(fā)生等級低、頻率高的地域。

    冬季整個渭河流域中度以上干旱發(fā)生頻率(圖9d)比其它季節(jié)都要低,整個區(qū)域內(nèi)天水、華山和吳旗地區(qū)屬于中度以上干旱發(fā)生頻率較高的地區(qū),其它地區(qū)中度以上干旱發(fā)生的頻率都介于3.7%~7.8%。相比較所有干旱等級其冬季干旱的發(fā)生頻率(圖7D),關(guān)中地區(qū)中西部屬于干旱發(fā)生頻率較高的地區(qū),但是該地區(qū)中度和極端干旱的發(fā)生頻率卻很低,說明冬季該地區(qū)雖然容易發(fā)生干旱,但其干旱程度都相對較低。

    3.3.3 基于月尺度的干旱發(fā)生強(qiáng)度 基于各站點SPEI 1數(shù)據(jù)和干旱強(qiáng)度的統(tǒng)計方法繪制出渭河流域干旱發(fā)生強(qiáng)度分布圖(見105頁圖10),最終統(tǒng)計得出在渭河流域干旱強(qiáng)度最大的是南部佛坪地區(qū),干旱強(qiáng)度達(dá)30%。其次干旱強(qiáng)度有明顯的東西差異,東部支流北洛河流域附近延安、吳旗等地干旱強(qiáng)度可達(dá)22%~23%,其它地區(qū)強(qiáng)度低于20%。天水和海源在近55年中發(fā)生連續(xù)干旱過程的次數(shù)最高,可達(dá)22 a。但是這兩個區(qū)域干旱過程連續(xù)的月數(shù)都是主要集中在連續(xù)3個月或者4個月,并無長期的連續(xù),反之吳旗附近存在長達(dá)9個月的一次連續(xù)干旱過程。

    4 結(jié)論與討論

    通過不同尺度SPI和SPEI指數(shù)對渭河流域近56 a的干旱時空分布特征的分析,得出以下結(jié)論:

    (1) SPEI相比SPI對渭河流域干旱的評估更加貼合實際。根據(jù)SPEI計算,發(fā)現(xiàn)尺度越大波動幅度越小,不同尺度的SPEI可以反映不同類型的干旱;渭河流域年平均SPEI在1991年發(fā)生突變,1991年后干旱加??;SPEI存在14 a左右的第一主周期變化,14 a周期內(nèi)還包括11 a和7 a左右的小周期,目前正處于相對濕潤的時期,接下來的一個周期有變干旱的趨勢。

    (2) 渭河流域年尺度干旱的發(fā)生頻率呈東高西低的分布,其中天水盆地附近干旱發(fā)生頻率較低,為21.8%,渭河流域其支流北洛河和涇河附近干旱發(fā)生的頻率較高,約為30.4%;發(fā)生中度以上干旱頻率較高的地區(qū)為西吉、武功和天水周邊,分別為14.28%、14.28%和12.5%,但是天水總的各個干旱等級的發(fā)生頻率卻相比其它地區(qū)低,說明天水干旱的發(fā)生雖然頻率低,但強(qiáng)度大,等級高,旱情嚴(yán)重。

    (3) 渭河流域季尺度秋季干旱發(fā)生頻率相比其它三個季節(jié)要高,為25.5%~37.5%,主要分布于陜西關(guān)中地區(qū)西部和中部;夏旱的高發(fā)區(qū)為岷縣和天水周邊,頻率分別為32.1%和30.4%;秋旱發(fā)生頻率較高區(qū)為北洛河延安周邊,頻率為37.5%;冬旱在關(guān)中平原頻率最高,頻率為27.3%~28.6%。就發(fā)生中度以上季節(jié)干旱而言,春季天水盆地周邊是發(fā)生中度以上春旱最嚴(yán)重的地區(qū),頻率為21.4%;夏季隴東黃土高原的平?jīng)鲋卸纫陨细珊祰?yán)重,頻率為16.1%;秋季天水盆地干旱發(fā)生等級高、頻率低,中度以上干旱發(fā)生頻率為19.6%,延安周邊干旱發(fā)生的頻率高但等級低,中度以上干旱發(fā)生頻率14.3%;冬季天水周邊發(fā)生中度以上干旱頻率較高,相反關(guān)中地區(qū)中西部干旱等級低。

    (4) 渭河流域干旱強(qiáng)度的分布總體上是由東向西減小,其中干旱強(qiáng)度最大的是南部佛坪地區(qū),干旱強(qiáng)度可達(dá)30%;天水和海源地區(qū)發(fā)生連續(xù)干旱過程的次數(shù)最多,可達(dá)22次,但連續(xù)月數(shù)最長的是吳旗地區(qū),統(tǒng)計期內(nèi)存在長達(dá)9個月的一次連續(xù)干旱過程。

    本文在研究干旱發(fā)生頻率、周期、突變情況和中度以上干旱發(fā)生頻率主要是基于年尺度(SPEI 12)而言,干旱強(qiáng)度和持續(xù)干旱月數(shù)是基于月尺度(SPEI 1),季節(jié)干旱發(fā)生頻率是基于季尺度(SPEI 3),最終得出的結(jié)論和前人的研究結(jié)果基本一致[20-24]。本文采用的SPEI雖然考慮了降水、溫度和蒸散發(fā),但大氣環(huán)流和下墊面特征,如風(fēng)速、地形等并未完全考慮,因此對于其成因和演變機(jī)制有待進(jìn)一步的研究。

    [1] 葉篤正.中國的全球變化預(yù)研究[M].北京:地震出版社,1992.

    [2] Easterling D R, Evans J L, Groismman P Y, et al. Observed variability and trends in extreme climate events:A brief review[J]. Bulletin of the American Meteorological Society, 2000,81(3):417-425.

    [3] 劉憲鋒,朱秀芳,潘耀忠,等.農(nóng)業(yè)干旱監(jiān)測研究進(jìn)展與展望[J].地理學(xué)報,2015,70(11):1835-1848.

    [4] Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation: Special Report of the Intergovernmental Panel on Climate Change[M]. Cambridge: Cambridge University Press, 2012.

    [5] 張 強(qiáng),韓蘭英,張立陽,等.論氣候變暖背景下干旱和干旱災(zāi)害風(fēng)險特征與管理策略[J].地球科學(xué)進(jìn)展,2014,29(1):80-91.

    [6] Palmer W C. Meteorological drought[M]. Washington, D C, USA: US Department of Commerce, Weather Bureau, 1965.

    [7] 么枕生,丁裕國.氣候統(tǒng)汁[M].北京:科學(xué)出版社,1990.

    [8] Jackson R D, Kustas W P. A reexamination of the crop water stress index[J]. Irrigation Science, 1988,9(4):309-317.

    [9] McKee T B, Doesken N J, Kleist J, et al. The relationship of drought frequency and duration to time scales[J]. Proceedings of the 8th Conference on Applied Climatology: American Meteorological Society Boston, M A, 1993:179-183.

    [10] Vicente-Serrano S M, Beguería 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(7):1696-1718.

    [11] Vicente-Serrano S M, Beguería S, Lorenzo-Lacruz J, et al. Performance of drought indices for ecological,agricultural,and hydrological applications[J]. Earth Interactions, 2012,16(10):1-27.

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

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

    [14] 李偉光,侯美亭,陳匯林,等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的華南干旱趨勢研究[J].自然災(zāi)害學(xué)報,2012,21(4):84-90.[15] 任培貴,張 勃,張調(diào)風(fēng),等.基于SPEI的中國西北地區(qū)氣象干旱變化趨勢分析[J].水土保持通報,2014,34(1):182-187,192.[16] 馬 瓊,張 勃,王 東,等.1960—2012年甘肅黃土高原干旱時空變化特征分析——基于標(biāo)準(zhǔn)化降水蒸散指數(shù)[J].資源科學(xué),2014,36(9):1834-1841.

    [17] 李景宜.陜西渭河流域主要環(huán)境問題及其治理[J].干旱區(qū)研究,2011,28(6):967-973.

    [18] Li Qi, Song Jinxi, Wei Anlei, et al. Changes in major factors affecting the ecosystem health of the Weihe River in Shaanxi Province, China[J]. Frontiers of Environmental Science & Engineering, 2013,7(6):875-885.

    [19] 畢彩霞,穆興民,趙廣舉,等.1960—2010年渭河流域旱澇時空分布特征[J].水土保持通報,2013,33(2):281-284,290.

    [20] 趙安周,劉憲鋒,朱秀芳,等.基于SWAT模型的渭河流域干旱時空分布[J].地理科學(xué)進(jìn)展,2015,34(9):1156-1166.

    [21] 陳昱潼,暢建霞,黃生志,等.基于PDSI的渭河流域干旱變化特征[J].自然災(zāi)害學(xué)報,2014,23(5):29-37.

    [22] 黃生志,黃 強(qiáng),王義民,等.基于SPI的渭河流域干旱特征演變研究[J].自然災(zāi)害學(xué)報,2015,24(1):15-22.

    [23] 馬明衛(wèi),宋松柏.渭河流域干旱指標(biāo)空間分布研究[J].干旱區(qū)研究,2012,4(29):681-691.

    [24] 田 甜,黃 強(qiáng),郭愛軍,等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的渭河流域干旱演變特征分析[J].水利發(fā)電學(xué)報,2016,35(2):16-27

    [25] 左德鵬,徐宗學(xué),程 磊,等.渭河流域潛在蒸散量時空變化及其突變特征[J].資源科學(xué),2011,33(5):975-982.

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

    [27] 楊曉靜,左德鵬,徐宗學(xué).基于標(biāo)準(zhǔn)化降水指數(shù)的云南省近55年旱澇演變特征[J].資源科學(xué),2014,36(3):473-480.

    [28] 中國氣象科學(xué)研究院,國家氣象中心,中國氣象局預(yù)測減災(zāi)司.氣象干旱等級(GB/T20481-2006)[S].北京:中國標(biāo)準(zhǔn)出版社,2006.

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

    [30] 周 丹,張 勃,任培貴,等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的陜西省近50 a干旱特征分析[J].自然資源學(xué)報,2014,29(4):677-688.

    Analysisofspatial-temporalcharacteristicsofdroughtbymultiplescalesinWeiheRiverBasin

    Liu Ruijuan1, Wang Xuhong1, Zhao Huijuan2

    (1.CollegeofUrbanandEnvironmentScience,NorthwestUniversity,Xi'an,Shaanxi710127,China; 2.KeyLaboratoryforEnvironmentalandUrbanSciences,ShenzhenGraduateSchool,PekingUniversity,Shenzhen,Guangdong518055,China)

    Based on the measured data of 23 climate stations during 1959—2014 in Weihe River Basin, compared and analyzed the applicability of the standardized precipitation index (SPI) and standardized precipitation evapotranspiration index (SPEI) in Weihe River Basin, obtained that stability of SPEI was more better. Afterwards using the SPEI, has researched the drought frequency and intensity, analyzed the spatial and temporal change characteristics of drought in recent 55 years by different time scale in Weihe River Basin. The results showed that: the more time scale, the smaller fluctuation of SPEI. The SPEI was mutated in 1991, hereafter the drought was become severe, and exsisted about 14 years and 7 years changing cycle in Weihe River Basin. The frequency of yearly scale drought was presented the distribution of high in east and low in west, the maximal drought frequency places above moderate drought were around Xiji, Wugong and Tianshui, the frequencies were 14.28%, 14.28% and 14.28% respectively. From the seasonal scale, the autumn drought frequency was higher than other seasons, the highest frequency of above moderate drought in spring, autumn and winter seasons was in Tianshui Basin, the frequency was 21.4%, 19.6% and 8.9% respectively. The highest frequency of above moderate drought in summer season was in Pingliang, its frequency was 16.1%. As a whole the distribution of monthly scale drought intensity was decreased from east to west, among them the biggest drought intensity was at the south Foping region, the drought intensity was 30%. In recent 55 years, the Tianshui and Haiyuan Regions were occured the maximal numbers of continuous drought process, can be reached 22 years. The month numbers as continuous drought process were mostly for 3 to 4 months, but there were one time of 9 months long continuous drought process in Wuqi station.Through the research for all kinds of drought characteristics by multiple time scales in this river basin, in order to provide guidance for drought control and resistance.

    standardized precipitation index (SPI); standardized precipitation evapotranspiration index (SPEI); multiple scale drought; spatial-temporal change characteristics; Weihe River Basin

    1000-7601(2017)03-0099-08doi:10.7606/j.issn.1000-7601.2017.03.16

    2016-05-13

    :2017-04-18

    :陜西省自然科學(xué)基礎(chǔ)研究計劃(2015JM4132)

    S423

    : A

    猜你喜歡
    渭河流域天水中度
    天水嬸與兩岸商貿(mào)
    華人時刊(2022年7期)2022-06-05 07:33:44
    天水地區(qū)的『秦與戎』
    小米手表
    電腦報(2020年12期)2020-06-30 19:56:42
    重返絲綢之路—從天水到青海湖
    美食(2018年10期)2018-10-18 08:10:58
    《天水之鏡像》
    LEEP治療中度宮頸上皮內(nèi)瘤變對HPV感染的影響
    基于RS/GIS 渭河流域植被覆蓋時空變化特征研究
    燕太子回國
    渭河流域香菜夏秋無公害栽培技術(shù)
    渭河流域糧食生產(chǎn)與灌溉農(nóng)業(yè)發(fā)展的相互作用關(guān)系分析
    中國水利(2015年5期)2015-02-28 15:12:36
    高潮久久久久久久久久久不卡| 欧美日韩精品网址| 黄色丝袜av网址大全| 国产欧美日韩一区二区精品| 国产精品久久视频播放| 日本免费a在线| 搡老妇女老女人老熟妇| 日本免费a在线| 亚洲七黄色美女视频| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影| 99国产综合亚洲精品| 91在线精品国自产拍蜜月 | 国内久久婷婷六月综合欲色啪| 国产综合懂色| 精品一区二区三区av网在线观看| 国产欧美日韩精品一区二区| 国产色爽女视频免费观看| 色老头精品视频在线观看| 国产av在哪里看| 麻豆国产97在线/欧美| 高清毛片免费观看视频网站| 真实男女啪啪啪动态图| 精品一区二区三区人妻视频| 亚洲国产欧洲综合997久久,| 男人的好看免费观看在线视频| 蜜桃久久精品国产亚洲av| 天天一区二区日本电影三级| 亚洲乱码一区二区免费版| 国产一区二区激情短视频| 国产亚洲av嫩草精品影院| 法律面前人人平等表现在哪些方面| 久久久久免费精品人妻一区二区| 色在线成人网| 狠狠狠狠99中文字幕| 夜夜夜夜夜久久久久| av专区在线播放| 日本黄大片高清| 国产精品免费一区二区三区在线| 美女cb高潮喷水在线观看| 久久伊人香网站| 又黄又爽又免费观看的视频| 男人舔女人下体高潮全视频| 嫩草影院入口| 日韩欧美在线二视频| 成人国产综合亚洲| 免费av毛片视频| 最近最新免费中文字幕在线| 中文在线观看免费www的网站| 一个人看的www免费观看视频| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 亚洲精品色激情综合| 黄片大片在线免费观看| 精品99又大又爽又粗少妇毛片 | 成年女人毛片免费观看观看9| 色综合婷婷激情| 男人舔奶头视频| 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器| 好看av亚洲va欧美ⅴa在| 亚洲精品在线观看二区| а√天堂www在线а√下载| 成人鲁丝片一二三区免费| 在线播放无遮挡| 成人18禁在线播放| 欧美绝顶高潮抽搐喷水| 久久99热这里只有精品18| 女人十人毛片免费观看3o分钟| 性欧美人与动物交配| 韩国av一区二区三区四区| 此物有八面人人有两片| 免费观看人在逋| 国产成人aa在线观看| 国产精品野战在线观看| 午夜福利18| 深爱激情五月婷婷| 中出人妻视频一区二区| 国产精品99久久99久久久不卡| 久久久成人免费电影| 亚洲精华国产精华精| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 热99在线观看视频| 给我免费播放毛片高清在线观看| 国产精品野战在线观看| 精华霜和精华液先用哪个| 午夜亚洲福利在线播放| 成人18禁在线播放| 亚洲 欧美 日韩 在线 免费| 久久久久久久午夜电影| x7x7x7水蜜桃| 久久精品国产亚洲av涩爱 | av在线天堂中文字幕| 欧美3d第一页| 一级黄色大片毛片| 亚洲精品美女久久久久99蜜臀| 久久久国产成人免费| 可以在线观看毛片的网站| 大型黄色视频在线免费观看| 悠悠久久av| 19禁男女啪啪无遮挡网站| 亚洲人与动物交配视频| 国产精品乱码一区二三区的特点| 日韩免费av在线播放| 男女午夜视频在线观看| 国产高潮美女av| 国产精品98久久久久久宅男小说| 一个人看视频在线观看www免费 | 久久久久久久精品吃奶| 一本一本综合久久| 老司机午夜福利在线观看视频| 国产成人av激情在线播放| 亚洲avbb在线观看| 欧美最黄视频在线播放免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲va日本ⅴa欧美va伊人久久| 免费av观看视频| 日韩欧美在线乱码| 女人高潮潮喷娇喘18禁视频| 一个人观看的视频www高清免费观看| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 欧美黑人欧美精品刺激| 国产成人aa在线观看| 中文资源天堂在线| 黄色女人牲交| 国产中年淑女户外野战色| 91九色精品人成在线观看| 亚洲精品成人久久久久久| 最近最新免费中文字幕在线| 久久久久久久久中文| 2021天堂中文幕一二区在线观| 日本a在线网址| 国产久久久一区二区三区| 亚洲电影在线观看av| 国产一区二区亚洲精品在线观看| 午夜福利成人在线免费观看| 日韩中文字幕欧美一区二区| 淫秽高清视频在线观看| 中出人妻视频一区二区| 12—13女人毛片做爰片一| 亚洲av电影在线进入| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区视频了| or卡值多少钱| 天堂影院成人在线观看| 午夜福利在线在线| 嫁个100分男人电影在线观看| 欧美最黄视频在线播放免费| 亚洲片人在线观看| 日韩亚洲欧美综合| 日韩国内少妇激情av| 国产欧美日韩一区二区三| 精品一区二区三区视频在线 | 男女视频在线观看网站免费| av专区在线播放| 国产午夜精品久久久久久一区二区三区 | 午夜激情福利司机影院| 免费看光身美女| 婷婷六月久久综合丁香| 99久久精品国产亚洲精品| 在线十欧美十亚洲十日本专区| 久久久成人免费电影| 成人av在线播放网站| 欧美乱码精品一区二区三区| 91麻豆精品激情在线观看国产| 最近在线观看免费完整版| 99久久99久久久精品蜜桃| 男人舔女人下体高潮全视频| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 18禁国产床啪视频网站| 少妇的丰满在线观看| 又紧又爽又黄一区二区| 国产精品三级大全| 黄色丝袜av网址大全| 日韩人妻高清精品专区| 免费看十八禁软件| a级一级毛片免费在线观看| 村上凉子中文字幕在线| 最近在线观看免费完整版| 无限看片的www在线观看| 丰满的人妻完整版| 少妇的丰满在线观看| 夜夜看夜夜爽夜夜摸| 久久久久久久午夜电影| 亚洲人成网站在线播| 精品一区二区三区视频在线观看免费| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| 欧美av亚洲av综合av国产av| 最近在线观看免费完整版| 一个人看视频在线观看www免费 | 少妇熟女aⅴ在线视频| 宅男免费午夜| 蜜桃亚洲精品一区二区三区| 精品国产亚洲在线| 亚洲av熟女| 国产探花在线观看一区二区| 色噜噜av男人的天堂激情| 高清毛片免费观看视频网站| 757午夜福利合集在线观看| 欧美黄色片欧美黄色片| 欧美3d第一页| 91久久精品电影网| 久久精品91无色码中文字幕| 国产亚洲精品一区二区www| 欧美午夜高清在线| 中出人妻视频一区二区| 日韩精品青青久久久久久| 天天躁日日操中文字幕| 欧美乱码精品一区二区三区| 九色成人免费人妻av| 成年女人毛片免费观看观看9| 国产99白浆流出| 精品免费久久久久久久清纯| 成人欧美大片| 久久久久久久久大av| 午夜激情欧美在线| 日本黄色视频三级网站网址| 狠狠狠狠99中文字幕| 99精品在免费线老司机午夜| 久久久成人免费电影| 免费在线观看亚洲国产| 国产99白浆流出| 狠狠狠狠99中文字幕| 免费人成在线观看视频色| 3wmmmm亚洲av在线观看| 亚洲欧美激情综合另类| 欧美日韩黄片免| 男女之事视频高清在线观看| 免费一级毛片在线播放高清视频| 久久中文看片网| 女人被狂操c到高潮| 久久国产精品影院| 国产精品,欧美在线| a级毛片a级免费在线| 国产亚洲精品久久久久久毛片| 一区福利在线观看| 亚洲aⅴ乱码一区二区在线播放| 免费看日本二区| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在 | 久久久久久久午夜电影| 午夜福利在线观看吧| 色综合欧美亚洲国产小说| 老司机深夜福利视频在线观看| 亚洲欧美日韩无卡精品| 99国产综合亚洲精品| 欧美丝袜亚洲另类 | 性色avwww在线观看| 久久久久久久亚洲中文字幕 | 欧美成人性av电影在线观看| 亚洲成av人片在线播放无| 国产精品,欧美在线| 成年女人毛片免费观看观看9| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 亚洲中文日韩欧美视频| 日本精品一区二区三区蜜桃| 制服人妻中文乱码| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 亚洲真实伦在线观看| 国产亚洲精品久久久com| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 精品熟女少妇八av免费久了| 99久久成人亚洲精品观看| 成人性生交大片免费视频hd| 免费在线观看日本一区| 日韩欧美国产一区二区入口| 丝袜美腿在线中文| 亚洲在线观看片| 久久精品人妻少妇| 天美传媒精品一区二区| 91麻豆av在线| 一个人看视频在线观看www免费 | 午夜福利视频1000在线观看| 淫妇啪啪啪对白视频| 国产老妇女一区| 美女高潮喷水抽搐中文字幕| 好男人在线观看高清免费视频| 午夜精品一区二区三区免费看| 两个人的视频大全免费| 九九热线精品视视频播放| 免费在线观看亚洲国产| 两性午夜刺激爽爽歪歪视频在线观看| 午夜免费成人在线视频| 欧美在线一区亚洲| 波多野结衣巨乳人妻| 男女床上黄色一级片免费看| 黑人欧美特级aaaaaa片| 久久中文看片网| 亚洲国产中文字幕在线视频| av国产免费在线观看| 少妇人妻精品综合一区二区 | 成人三级黄色视频| 国产麻豆成人av免费视频| 欧美三级亚洲精品| 免费看十八禁软件| 又紧又爽又黄一区二区| 一边摸一边抽搐一进一小说| 又紧又爽又黄一区二区| 免费看十八禁软件| 欧美区成人在线视频| 欧美在线一区亚洲| 免费无遮挡裸体视频| 天美传媒精品一区二区| 免费看光身美女| 欧美性感艳星| 好看av亚洲va欧美ⅴa在| 国产三级黄色录像| 国产伦精品一区二区三区四那| 美女黄网站色视频| 九色国产91popny在线| 久久精品国产亚洲av香蕉五月| 欧美3d第一页| 亚洲精品日韩av片在线观看 | 欧美三级亚洲精品| 久久性视频一级片| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 熟女电影av网| 亚洲成av人片免费观看| 不卡一级毛片| 午夜精品久久久久久毛片777| 亚洲成人免费电影在线观看| 黄片大片在线免费观看| 国产一级毛片七仙女欲春2| 十八禁人妻一区二区| 免费av不卡在线播放| 亚洲av熟女| 香蕉久久夜色| 婷婷精品国产亚洲av在线| 成人性生交大片免费视频hd| 国产乱人视频| 午夜激情福利司机影院| 欧美日本亚洲视频在线播放| 美女高潮的动态| 亚洲精品在线观看二区| 首页视频小说图片口味搜索| 美女免费视频网站| 99热6这里只有精品| 好男人电影高清在线观看| 最近在线观看免费完整版| 亚洲一区二区三区色噜噜| 给我免费播放毛片高清在线观看| 我的老师免费观看完整版| 久久午夜亚洲精品久久| 99久久久亚洲精品蜜臀av| 精品欧美国产一区二区三| 精品不卡国产一区二区三区| 黄片大片在线免费观看| 麻豆久久精品国产亚洲av| 久久九九热精品免费| 国产极品精品免费视频能看的| 日韩高清综合在线| 熟女电影av网| 一区二区三区免费毛片| 在线国产一区二区在线| 国产精品一及| 母亲3免费完整高清在线观看| 国产精品嫩草影院av在线观看 | 国产高清视频在线观看网站| 桃色一区二区三区在线观看| 三级国产精品欧美在线观看| 色精品久久人妻99蜜桃| 两性午夜刺激爽爽歪歪视频在线观看| 丰满的人妻完整版| 色av中文字幕| 午夜福利在线在线| 一个人免费在线观看电影| 麻豆一二三区av精品| 十八禁网站免费在线| 精品人妻一区二区三区麻豆 | 在线免费观看的www视频| www国产在线视频色| 日本免费一区二区三区高清不卡| 一边摸一边抽搐一进一小说| 成人永久免费在线观看视频| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 国产野战对白在线观看| 操出白浆在线播放| 亚洲电影在线观看av| 中亚洲国语对白在线视频| 中文字幕高清在线视频| 亚洲激情在线av| 国产v大片淫在线免费观看| 伊人久久大香线蕉亚洲五| 天天躁日日操中文字幕| 国内精品一区二区在线观看| 亚洲av电影在线进入| 精品久久久久久久久久免费视频| ponron亚洲| 国模一区二区三区四区视频| 香蕉久久夜色| 高清日韩中文字幕在线| 国产淫片久久久久久久久 | 俄罗斯特黄特色一大片| 精华霜和精华液先用哪个| 国产精品三级大全| 99精品久久久久人妻精品| 精品国内亚洲2022精品成人| 国产免费一级a男人的天堂| 亚洲熟妇熟女久久| 亚洲人成电影免费在线| 成人一区二区视频在线观看| 1000部很黄的大片| 亚洲无线在线观看| 国产视频内射| 伊人久久精品亚洲午夜| 国产一区二区激情短视频| 国产在视频线在精品| 国产真实伦视频高清在线观看 | 午夜a级毛片| 少妇人妻精品综合一区二区 | 久久草成人影院| 日韩人妻高清精品专区| 给我免费播放毛片高清在线观看| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全免费视频| 欧美高清成人免费视频www| 男女床上黄色一级片免费看| 亚洲天堂国产精品一区在线| 一区福利在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 观看美女的网站| 国产精品一区二区三区四区免费观看 | a级一级毛片免费在线观看| 舔av片在线| 老司机深夜福利视频在线观看| 免费搜索国产男女视频| 全区人妻精品视频| 亚洲成人久久爱视频| 男女午夜视频在线观看| 国产精品久久电影中文字幕| 久久久久久久久大av| 日韩大尺度精品在线看网址| a级毛片a级免费在线| 久久精品国产亚洲av涩爱 | 亚洲天堂国产精品一区在线| 亚洲精品美女久久久久99蜜臀| 午夜免费男女啪啪视频观看 | 久久香蕉精品热| 欧美成狂野欧美在线观看| 成人国产一区最新在线观看| 欧美一级毛片孕妇| 高清日韩中文字幕在线| 精品久久久久久久人妻蜜臀av| 免费大片18禁| 亚洲专区国产一区二区| 欧美+日韩+精品| 国产精品久久久久久亚洲av鲁大| 美女高潮的动态| 亚洲片人在线观看| 人妻丰满熟妇av一区二区三区| 久久香蕉国产精品| 91在线观看av| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在 | 欧美色视频一区免费| 国产精品影院久久| 一卡2卡三卡四卡精品乱码亚洲| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添小说| 香蕉丝袜av| av福利片在线观看| 麻豆国产av国片精品| 免费看a级黄色片| 亚洲不卡免费看| 亚洲av成人av| 韩国av一区二区三区四区| 欧美在线一区亚洲| 老汉色∧v一级毛片| 很黄的视频免费| 一个人观看的视频www高清免费观看| 久久精品91蜜桃| 成人av在线播放网站| 精品久久久久久,| 成人18禁在线播放| 欧美一区二区国产精品久久精品| 欧美+日韩+精品| 亚洲精品久久国产高清桃花| 伊人久久大香线蕉亚洲五| 91麻豆av在线| 波多野结衣高清作品| 91在线观看av| 国模一区二区三区四区视频| 精品一区二区三区视频在线观看免费| av片东京热男人的天堂| 欧美午夜高清在线| 午夜日韩欧美国产| 最近最新免费中文字幕在线| 日韩人妻高清精品专区| 亚洲人成网站在线播| 国产精品久久久久久精品电影| 亚洲欧美日韩高清在线视频| 日韩欧美免费精品| 夜夜看夜夜爽夜夜摸| aaaaa片日本免费| 五月伊人婷婷丁香| 精品久久久久久久久久免费视频| 国产精品女同一区二区软件 | 一个人看视频在线观看www免费 | 午夜影院日韩av| 精品久久久久久久久久久久久| 亚洲精品粉嫩美女一区| 国产精品国产高清国产av| 真人一进一出gif抽搐免费| 欧美日本亚洲视频在线播放| 国产综合懂色| 欧美性感艳星| 亚洲中文字幕日韩| 国产97色在线日韩免费| 黄片小视频在线播放| 久久精品国产综合久久久| 国产精品 国内视频| 精品一区二区三区人妻视频| 99久国产av精品| 88av欧美| 日韩国内少妇激情av| 欧美乱妇无乱码| 亚洲欧美日韩无卡精品| 午夜精品一区二区三区免费看| 欧美日韩乱码在线| 精品国产超薄肉色丝袜足j| 97超级碰碰碰精品色视频在线观看| 色综合亚洲欧美另类图片| 人人妻人人看人人澡| 久久天躁狠狠躁夜夜2o2o| 国产三级中文精品| 成人18禁在线播放| 中文字幕熟女人妻在线| 在线a可以看的网站| 国内精品久久久久久久电影| 色综合亚洲欧美另类图片| 在线国产一区二区在线| 日韩欧美精品免费久久 | 18禁在线播放成人免费| 日本成人三级电影网站| av片东京热男人的天堂| 97超视频在线观看视频| 老司机深夜福利视频在线观看| 欧美成狂野欧美在线观看| 91av网一区二区| 免费观看精品视频网站| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区精品小视频在线| 亚洲久久久久久中文字幕| 欧美一区二区国产精品久久精品| 欧美最新免费一区二区三区 | 日韩欧美在线二视频| 成人特级av手机在线观看| 国产精品98久久久久久宅男小说| 99热6这里只有精品| 久久久久久久久大av| 在线观看舔阴道视频| 午夜老司机福利剧场| 99国产极品粉嫩在线观看| 日本在线视频免费播放| 亚洲性夜色夜夜综合| 亚洲 欧美 日韩 在线 免费| 中出人妻视频一区二区| 亚洲精品国产精品久久久不卡| 一个人看的www免费观看视频| x7x7x7水蜜桃| 岛国视频午夜一区免费看| 国产精品乱码一区二三区的特点| 在线观看日韩欧美| 国内揄拍国产精品人妻在线| 在线国产一区二区在线| 国产精华一区二区三区| 一级a爱片免费观看的视频| 国内精品久久久久精免费| 91在线精品国自产拍蜜月 | 色综合欧美亚洲国产小说| 国产三级黄色录像| 国产精品亚洲一级av第二区| www国产在线视频色| 国产成人福利小说| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 国产av不卡久久| 免费av不卡在线播放| 最新在线观看一区二区三区| 国产国拍精品亚洲av在线观看 | 国产精品女同一区二区软件 | 亚洲av不卡在线观看| 欧美区成人在线视频| 国产男靠女视频免费网站| 欧美不卡视频在线免费观看| 日本五十路高清| 桃红色精品国产亚洲av| 亚洲人成伊人成综合网2020| 日韩欧美精品v在线| 色综合亚洲欧美另类图片| 51午夜福利影视在线观看| 一区二区三区国产精品乱码| 国产一区二区在线av高清观看| 午夜免费男女啪啪视频观看 | 啪啪无遮挡十八禁网站| 国产精品99久久99久久久不卡| 国产精品 国内视频| 亚洲国产高清在线一区二区三| 午夜福利视频1000在线观看| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 亚洲色图av天堂|