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

    基于標(biāo)準(zhǔn)化降水指數(shù)的陜西省干旱時空變化特征分析

    2017-11-01 23:03:55解建倉胡彥華姜仁貴
    農(nóng)業(yè)工程學(xué)報 2017年17期
    關(guān)鍵詞:趨勢標(biāo)準(zhǔn)化

    李 斌,解建倉,胡彥華,姜仁貴

    ?

    基于標(biāo)準(zhǔn)化降水指數(shù)的陜西省干旱時空變化特征分析

    李 斌1,2,解建倉1,胡彥華2,姜仁貴1※

    (1. 西安理工大學(xué)水利水電學(xué)院,西安 710048; 2. 陜西省水利廳,西安 710004)

    干旱作為當(dāng)今全球范圍內(nèi)危害嚴(yán)重且頻發(fā)的一種自然災(zāi)害之一,對人類和自然系統(tǒng)都產(chǎn)生較大影響。該文基于陜西省長歷時的典型測站逐日降水資料,采用標(biāo)準(zhǔn)化降水指數(shù)(standardized precipitation index,SPI)、干旱頻率、干旱站次比和干旱強(qiáng)度等干旱指標(biāo),揭示陜西省干旱時間與空間尺度上的變化特征,為區(qū)域水資源合理調(diào)配提供參考。結(jié)果表明:1)研究區(qū)域整體上干旱范圍有擴(kuò)大的趨勢,研究時段內(nèi)干旱較為頻繁且持續(xù)時間較長,其中,區(qū)域內(nèi)各個站點干旱發(fā)生頻率均值為70.0%,輕旱和中旱發(fā)生頻率較高,重旱和特旱發(fā)生頻率相對較低;2)干旱站次比呈現(xiàn)出季節(jié)性特點,其中,春季、冬季干旱站次比呈增加趨勢,夏季、秋季干旱站次比呈減少趨勢;3)干旱強(qiáng)度變化不明顯,其中,春季、冬季的干旱強(qiáng)度呈增強(qiáng)趨勢,夏季、秋季干旱強(qiáng)度呈減弱趨勢;干旱站次比和干旱強(qiáng)度變化趨勢基本一致;4)干旱呈現(xiàn)出較明顯的區(qū)域變化特征,其中,既有覆蓋全省范圍的全域性干旱,也有局部地區(qū)干旱,并且干旱頻率、干旱站次比和干旱強(qiáng)度空間上分布上也存在差異。

    氣象;干旱;降水;標(biāo)準(zhǔn)化降水指數(shù);陜西省

    0 引 言

    干旱作為全球范圍內(nèi)影響范圍最大、災(zāi)害損失最為嚴(yán)重的自然災(zāi)害之一,對區(qū)域的農(nóng)業(yè)、經(jīng)濟(jì)、生態(tài)和環(huán)境造成較為嚴(yán)重的負(fù)面影響[1-2]。中國是1個干旱頻發(fā)的國家,尤其是全球氣候變化,引起區(qū)域降雨減少和蒸發(fā)增加,進(jìn)而使干旱呈現(xiàn)出發(fā)生頻率增加、強(qiáng)度加大和范圍不斷擴(kuò)大的變化趨勢,受到學(xué)者的廣泛關(guān)注,是研究熱點問題[3-6]。陜西省地處中國內(nèi)陸腹地,屬典型的大陸性季風(fēng)氣候,素有“十年九旱”之說,全省每年都有不同程度的旱災(zāi)發(fā)生。特殊的地理位置和復(fù)雜的大氣環(huán)流特征,使得旱災(zāi)成為陜西省主要的自然災(zāi)害之一。根據(jù)國家防汛抗旱總指揮部和水利部聯(lián)合發(fā)布的《2016年中國水旱災(zāi)害公報》統(tǒng)計數(shù)據(jù)表明:2016年陜西省作物因旱受災(zāi)面積高達(dá)2.40×105hm2,成災(zāi)面積達(dá)到1.11× 105hm2。干旱給城鄉(xiāng)居民生活和工農(nóng)業(yè)生產(chǎn)造成不同程度的影響,嚴(yán)重制約著全省經(jīng)濟(jì)社會的健康持續(xù)發(fā)展。

    干旱影響因素較多,為了更好對干旱狀況進(jìn)行定量化描述,學(xué)者提出許多用于分析干旱發(fā)生頻率、干旱強(qiáng)度和歷時的干旱指標(biāo)[7-10]。目前在評價氣象干旱時,通常選取標(biāo)準(zhǔn)化降水指數(shù)(standardized precipitation index,SPI)、土壤相對濕度指數(shù)、Palmer干旱指數(shù)、綜合氣象干旱指數(shù)、溫度干旱植被指數(shù)、農(nóng)業(yè)干旱參考指數(shù)等[11-14]。基于站點降水量的SPI以其簡單實用和時間序列長而受到眾多青睞[15-16]。張巧鳳等[17]采用SPI分析錫林郭勒盟干旱特征,發(fā)現(xiàn)該區(qū)域主要為輕旱和中旱,重旱和特旱頻率較低;高西寧等[18]采用SPI分析遼寧省近54 a干旱時空規(guī)律,發(fā)現(xiàn)春季干旱和秋季干旱發(fā)生的范圍較大,持續(xù)時間較長,夏季干旱發(fā)生的范圍較??;劉梅等[19]發(fā)現(xiàn)20世紀(jì)90年代以后渭河流域干旱事件的頻率、強(qiáng)度和烈度均有增加趨勢;孫智輝等[20]發(fā)現(xiàn)陜西省黃土高原地區(qū)年度和夏、秋、冬季干旱強(qiáng)度都有不同程度降低,春季干旱有增強(qiáng)趨勢。

    SPI通過概率密度函數(shù)求解累積概率,進(jìn)一步將累積概率進(jìn)行標(biāo)準(zhǔn)化處理,消除了降水的時空分布差異,具有穩(wěn)定的計算特性且計算相對簡單,在國內(nèi)外有廣泛的應(yīng)用[14-20],但在陜西省應(yīng)用相對較少。本文選取SPI作為干旱指標(biāo),基于近40 a日值降水?dāng)?shù)據(jù),從氣象干旱角度,分析陜西省干旱發(fā)生頻率和強(qiáng)度的時空變化特征,揭示區(qū)域干旱時空演變規(guī)律,可為區(qū)域抗旱減災(zāi)提供參考。

    1 資料與方法

    1.1 研究區(qū)域與數(shù)據(jù)來源

    陜西省位于105°29′-111°15′ E,31°42′-39°35′ N之間,地處中國東南濕潤地區(qū)到西北干旱地區(qū)的過渡帶,地勢特點是南北高、中間低,由西向東傾斜,屬于典型的大陸性季風(fēng)氣候區(qū)。全省分為陜北高原、關(guān)中平原和陜南秦巴山地3個自然區(qū)域。陜北黃土高原區(qū)由深厚的黃土覆蓋,海拔在900~1 500 m之間,占全省土地總面積約45%;關(guān)中平原區(qū)由渭河干流以及兩岸支流日久沖積、沉積形成,海拔在300~800 m之間,占全省土地總面積約19%;陜南秦巴山地處于秦嶺和大巴山系,海拔200~1 200 m,占全省總面積約36%[21]。

    本文所采用的數(shù)據(jù)為1971-2013年陜西省內(nèi)的34個氣象站點逐日降水資料,數(shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.cma.gov.cn)。采用SPI計算年尺度和季節(jié)尺度的干旱指數(shù),分析陜西省年度和季節(jié)干旱變化特征,統(tǒng)計其干旱頻率、干旱站次比和干旱強(qiáng)度,分析研究區(qū)干旱的時空變化特征。研究區(qū)域及氣象站空間分布如圖1所示。

    圖1 研究區(qū)域及氣象站空間分布

    1.2 研究方法

    1.2.1 標(biāo)準(zhǔn)化降水指數(shù)(SPI)

    由于不同時間、不同地區(qū)降水量變化幅度很大,直接用降水量很難在不同時空尺度上相互比較,而且降水分布是一種偏態(tài)分布,不是正態(tài)分布,所以在降水分析中,采用Γ分布概率來描述降水量的變化,然后再經(jīng)正態(tài)標(biāo)準(zhǔn)化求得SPI值[22-25]。SPI的計算步驟為:

    假設(shè)某時段降水量服從Γ分布,則其概率密度計算公式為:

    式中為形狀參數(shù),為尺度參數(shù),采用極大似然估計(maximum likelihood estimation,MLE)可知

    式(2)中計算式如下

    式中為降水?dāng)?shù)據(jù)的時間序列長度。

    確定概率密度函數(shù)中的參數(shù)后,對于某一年的降水量0可求出多年平均降水量小于該年降水量事件的概率為

    降水量為0時的概率為

    (=0)=/(6)

    式中為樣本數(shù),為總樣本數(shù)。對Γ分布概率從一般正態(tài)分布轉(zhuǎn)為標(biāo)準(zhǔn)正態(tài)分布,近似求解可得

    根據(jù)干旱等級標(biāo)準(zhǔn)(GB/T20481-2006氣象干旱等級),?1.0

    1.2.2 干旱評估指標(biāo)

    除了SPI[26],本文還采用干旱發(fā)生頻率、干旱發(fā)生站次比和干旱強(qiáng)度等干旱指標(biāo)分析陜西省干旱變化特征[27]。

    1)干旱頻率(P

    P是用來評價研究區(qū)某站一段時期內(nèi)發(fā)生干旱的頻繁程度,計算公式為[17,26]

    P=(S/)×100% (8)

    式中為計算總年數(shù),=43;S為站出現(xiàn)干旱的年數(shù),文中為1到34。根據(jù)不同干旱等級的發(fā)生年數(shù)計算各自的發(fā)生頻率。

    2)干旱站次比(Q

    Q為研究區(qū)內(nèi)發(fā)生干旱站數(shù)占全部站數(shù)的比例,用來評價干旱影響范圍,計算公式為[17]

    θ=(Z/) (9)

    式中為總站數(shù),Z為年發(fā)生干旱的測站數(shù)量,文中為1到43。干旱的影響范圍定義:當(dāng)絕大部分測站數(shù)據(jù)指標(biāo)均顯示為干旱時,即θ≥70%時,定義為全域性干旱;當(dāng)50%≤θ<70%,定義為區(qū)域性干旱;當(dāng)30%≤θ<50%,定義為部分地區(qū)干旱;當(dāng)10%≤θ<30%,定義為局部地區(qū)干旱;當(dāng)θ<10%時,定義為無明顯干旱。

    3)干旱強(qiáng)度()

    S用來評價干旱嚴(yán)重程度,計算公式為[17,27]

    式中SPI為年發(fā)生干旱的站的SPI值,為發(fā)生干旱的站數(shù),當(dāng)0.5≤S<1時為輕旱、1≤S<1.5為中旱,S≥1.5為重旱,S越大,干旱越嚴(yán)重。

    4)變化趨勢率

    變化趨勢率通常為測站歷年氣候要素時間序列擬合直線的斜率,降水的變化趨勢率單位為mm/10 a。一般來講,降水的氣候趨勢用一次直線方程和二次曲線方程就能滿足擬合要求。

    2 結(jié)果與分析

    2.1 干旱年際變化特征

    SPI 12表示12個月尺度的SPI值,采用SPI 12可以對陜西省干旱的年際變化特征進(jìn)行分析,如圖2所示,可以看出陜西省在研究時段內(nèi)干旱變化的總體情況。其中,1971-1973年這3 a之中每年都有干旱發(fā)生,1993年之后的12 a內(nèi)發(fā)生干旱頻率也相對較高(7 a發(fā)生干旱);1971年3月-1973年12月這33個月;1993年9月-1996年12月這39個月每個月都發(fā)生干旱;干旱持續(xù)時間最長的時間出現(xiàn)在1999年2月,直至2003年7月結(jié)束。

    干旱發(fā)生的主要年份有1971-1973年、1977-1978年、1986-1987年、1991-1992年、1993-1997年、1999-2002年、2006-2009年。從干旱發(fā)生的等級來看,最嚴(yán)重的年份為1997、1998年,其中1997年11月-1998年12月的SPI 12值達(dá)到了?1.5,等級為重旱。綜上可知,從全省范圍上看,陜西省在研究時段內(nèi)干旱發(fā)生較頻繁,持續(xù)時間也較長。

    注:SPI 12表示12個月的標(biāo)準(zhǔn)化降水指數(shù)。

    2.2 干旱年度變化特征

    2.2.1 干旱發(fā)生頻率(P)

    基于SPI 12可計算得到陜西省不同區(qū)域干旱發(fā)生頻率,如表1所示。各個站點發(fā)生干旱的頻率約為70.0%,其中,發(fā)生輕旱、中旱、重旱和特旱的頻率分布為23.7%、23.3%、14.2%和8.8%。從不同等級干旱發(fā)生的頻率可知,隨著干旱強(qiáng)度增加,發(fā)生頻率呈現(xiàn)遞減的趨勢。

    表1 陜西省及省內(nèi)各區(qū)域的干旱發(fā)生頻率

    由表1可知,陜西省各個地區(qū)干旱發(fā)生頻率接近,陜北平均為70%,關(guān)中平均為72.1%,陜南平均為68%;不同區(qū)域發(fā)生不同類型的干旱的最大頻率則存在較大差異。

    2.2.2 干旱站次比(Q)

    采用SPI 12時間序列,根據(jù)式(8)可計算得到陜西省干旱站次比,結(jié)果表明:陜西省在研究時段內(nèi)每年均出現(xiàn)了不同等級的干旱;1971-2013年,共發(fā)生全域性干旱25次,比例最大;區(qū)域性的干旱事件相對較少,占到陜西省干旱的一半左右(12次),部分地區(qū)和局部地區(qū)的較少,分別為4次和2次。其中1992年、2000年干旱站次比為100%,1973年、1986-1987年、1997-1999年和2001年干旱站次比達(dá)到90%。20世紀(jì)90年代干旱嚴(yán)重,1991-2000年除1993年發(fā)生區(qū)域性干旱外,其余年份均發(fā)生全域性干旱。

    由圖3可知,陜西省各年代中發(fā)生不同等級干旱的年均站次比存在差異,而上世紀(jì)90年代年均干旱站次比最高,達(dá)到了70%。1971-2013年,干旱站次比的線性傾向率為0.36%/10 a,表明從長期來看陜西省干旱站次比呈增加趨勢。干旱范圍有擴(kuò)大的趨勢;輕旱、中旱站次比呈微弱減少趨勢;特旱、重旱站次比呈增加趨勢。

    2.2.3干旱發(fā)生強(qiáng)度()

    基于各站SPI 12序列計算陜西省干旱強(qiáng)度S,如圖4所示,結(jié)果表明:陜西省在研究時段內(nèi)年均干旱S為0.87,其中1998年的S達(dá)到1.79,為重度干旱。70年代6 a為輕度干旱,4 a為中旱;80年代6 a為輕度干旱,4 a為中旱;90年代6 a為輕度干旱,2 a為中旱,2 a為重旱;2000年以后的13 a中,9 a為輕度干旱,4 a為中旱。43 a中,有27 a為輕度干旱,14 a為中旱,2 a為重旱,說明陜西省的干旱主要為輕旱和中旱。陜西省干旱強(qiáng)度變化趨勢不明顯。20世紀(jì)90年代平均干旱強(qiáng)度最大。

    圖3 陜西省各年代年均干旱站次比

    圖4 陜西省干旱強(qiáng)度變化

    2.3 季節(jié)性干旱特征分析

    1971-2013年陜西省發(fā)生的季節(jié)性干旱如表2所示,采用3個月的SPI值SPI 3來分析陜西省季節(jié)性的干旱。結(jié)果表明:43 a間,陜西省共發(fā)生春旱25次,18 a無明顯春旱;發(fā)生夏旱22次,21 a無明顯夏旱;發(fā)生秋旱24次,19 a無明顯秋旱;發(fā)生冬旱27次,16 a無明顯秋旱。

    表2 陜西省1971-2013年季節(jié)干旱次數(shù)統(tǒng)計

    為了進(jìn)一步揭示陜西省季節(jié)性干旱變化特征,將1971-2013年劃分為5個時段,計算不同年代不同季節(jié)和全年的干旱站次比和干旱強(qiáng)度情況,計算結(jié)果如表3所示。

    1)春旱。由表2可知,1971-2013年發(fā)生25次春旱,全域性春旱8次,其中,20世紀(jì)80和90年代分別發(fā)生2次全域性春旱;2000年后共發(fā)生4次全域性春旱,其中1984、2000、2004和2005年春旱最嚴(yán)重,干旱站次比達(dá)到79%,其次為1995年,春旱站次比為77%。春季干旱站次比的線性傾向率為4.317/10 a,表明陜西省春季干旱范圍有擴(kuò)大趨勢。春旱強(qiáng)度在0.75~1.71之間波動變化,與干旱站次比變化曲線較為類似。

    表3 各年代陜西省不同季節(jié)干旱站次比和和干旱強(qiáng)度

    2)夏旱。相比春旱,夏旱發(fā)生頻率也較為頻繁,主要集中在1972-1977、1979、1981-1982、1985-1986、1991、1994-1997、1999、2001、2004-2008年,其中,1974、1976、1982、1999和2004年為夏季中旱年,1995和1997年為夏季重旱年。夏旱站次比在2%~79%之間變化,占總年數(shù)的4.7%。共發(fā)生4次全域性夏旱,分別發(fā)生在上世紀(jì)1976、1982年和90年代的1995和1997年,80年代和21世紀(jì)以來未發(fā)生全域性夏旱。夏旱強(qiáng)度在0.73~1.89之間波動變化,上世紀(jì)90年代最重,其中1997和1995年干旱強(qiáng)度分別為1.89和1.55。夏旱強(qiáng)度的線性傾向率為?0.01/10 a,表明陜西省夏季干旱強(qiáng)度呈減弱趨勢。

    3)秋旱。1971-2013年發(fā)生秋旱24次,集中在1971-1972、1976-1977、1979、1986-1991、1993-1995、1997-2000、2006-2007、2010和2013年,其中,1972、1977、1986、1994、1999和2002年為秋季中旱年,1991、1997和1998年為秋季重旱年。發(fā)生4次全域性秋旱,其中,1991年秋旱站次比最高,達(dá)到79%,1997和1998年秋旱站次比為77%,1986 年全域性秋旱站次比為72%。秋旱站次比在0~79%之間變化,占總年數(shù)的9.3%。90年代陜西省秋旱站次比最高,年平均秋旱站次比達(dá)到73.5%。秋旱站次比的線性傾向率為?1.439%/10 a,表明陜西省秋旱范圍呈微弱減少趨勢。

    4)冬旱。1971-2013年發(fā)生冬旱次數(shù)高達(dá)27次,集中在1973、1976、1978-1988、1990-1991、1995、1997-2001、2005-2006、2008、2010-2013年。發(fā)生8次全域性冬旱,其中,1982和1998年冬旱站次比最高,達(dá)到79%,1985年冬旱站次比為77%,1979和1995年全域性冬旱站次比為74%。上世紀(jì)80年代冬旱站次比最高,年平均達(dá)到69.7%。冬旱站次比的線性傾向率為2.891%/10 a,干旱強(qiáng)度在0.66~2.13之間波動變化,其中,1998、1979和2013年冬旱最重,干旱強(qiáng)度分別為2.13、1.37和1.37。上述結(jié)果表明陜西省冬旱范圍有擴(kuò)大趨勢。

    姜仁貴等[14]研究表明陜西省總體呈現(xiàn)出干旱化趨勢增強(qiáng)的特征。孫智輝等[20]采用SPI分析陜西黃土高原地區(qū)1971-2010年干旱月、季度和年度干旱變化特征。結(jié)果表明:不同時段和不同區(qū)域干旱變化趨勢存在一定差異。周丹等[28]研究結(jié)果表明陜西省干旱發(fā)生頻率呈明顯的增長趨勢;全省以20世紀(jì)90年代干旱最為嚴(yán)重,2000年以來的干旱次之。王曉峰等[29]基于1960-2015年9個測站月降雨數(shù)據(jù),采用Copula函數(shù)和游程理論分析陜北黃土高原地區(qū)干旱事件發(fā)生的強(qiáng)度和歷時,結(jié)果表明研究區(qū)內(nèi)不同測站發(fā)生干旱的強(qiáng)度、歷時及其聯(lián)合重現(xiàn)期變化趨勢存在差異。

    3 結(jié) 論

    1)1971-2013年期間,陜西省干旱范圍總體上有擴(kuò)大的趨勢;其中,春季和冬季干旱范圍呈現(xiàn)擴(kuò)大趨勢;夏季和秋季干旱范圍有減小趨勢。

    2)研究區(qū)各站點發(fā)生干旱頻率均值為70.0%,輕旱、中旱、重旱和特旱的頻率分布為23.7%、23.3%、14.2%和8.8%;輕旱和中旱頻率較高,而重旱和特旱頻率較低。干旱發(fā)生頻率隨著干旱強(qiáng)度增加呈現(xiàn)遞減趨勢。

    3)研究區(qū)干旱站次比在17.6%~100%之間變化;輕旱、中旱站次比呈微弱減少趨勢,特旱、重旱站次比呈增加趨勢。陜西整體干旱強(qiáng)度變化趨勢不明顯;20世紀(jì)90年代平均干旱強(qiáng)度最大。干旱強(qiáng)度和干旱站次比變化趨勢基本一致。

    [1] Schwalm C R, Anderegg W R L, Michalak A M, et al. Global patterns of drought recovery[J]. Nature, 2017, 548(7666): 202-205.

    [2] Sheffield J, Wood E F, Roderick M L. Little change in global drought over the past 60 years[J]. Nature, 2012, 491(7424): 435-438.

    [3] 賈建英,賀楠,韓蘭英,等. 基于自然災(zāi)害風(fēng)險歷來和ArcGIS的西南地區(qū)玉米干旱風(fēng)險分析[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(4):152-159.

    Jia Jianying, He Nan, Han Lanying, et al. Analysis on drought risk of maize in Southwest China based on natural disaster risk theory and ArcGIS[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(4): 152-159. (in Chinese with English abstract)

    [4] Craine J M, Ocheltree T W, Nippert J B, et al. Global diversity of drought tolerance and grassland climate-change resilience[J]. Nature Climate Change, 2013, 3(1): 63-67.

    [5] 張建平,劉宗元,王靖,等. 西南地區(qū)綜合干旱監(jiān)測模型構(gòu)建與驗證[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(5):102-107.

    Zhang Jianping, Liu Zongyuan, Wang Jing, et al. Construction and validation of comprehensive drought monitoring model in Southwest China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(5): 102-107. (in Chinese with English abstract)

    [6] Liu Zhiping, Wang Yunqiang, Shao Mingan, et al. Spatiotemporal analysis of multiscalar drought characteristics across the Loess Plateau of China[J]. Journal of Hydrology, 2016, 534: 281-299.

    [7] Vicente-Serrano S M, Van der Schrier G, Beguería S, et al. Contribution of precipitation and reference evapotranspiration to drought indices under different climates[J]. Journal of Hydrology, 2015, 526: 42-54.

    [8] Touma D, Ashfaq M, Nayak MA, et al. A multi-model and multi-index evaluation of drought characteristics in the 21st century[J]. Journal of Hydrology, 2015, 526: 196-207.

    [9] 王兆禮,黃澤勤,李軍,等. 基于SPEI和NDVI的中國流域尺度氣象干旱及植被分布時空演變[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(14):177-186.

    Wang Zhaoli, Huang Zeqin, Li Jun, et al. Assessing impacts of meteorological drought on vegetation at catchment scale in China based on SPEI and NDVI[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(14): 177-186. (in Chinese with English abstract)

    [10] 趙煥,徐宗學(xué),趙捷. 基于CWSI及干旱稀缺程度的農(nóng)業(yè)干旱指數(shù)構(gòu)建及應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(9):116-125.

    Zhao Huan, Xu Zongxue, Zhao Jie. Development and application of agricultural drought index based on CWSI and drought event rarity[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(9): 116-125. (in Chinese with English abstract)

    [11] 劉宗元,張建平,羅紅霞,等. 基于農(nóng)業(yè)干旱參考指數(shù)的西南地區(qū)玉米干旱時空變化分析[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(2):106-115.

    Liu Zongyuan, Zhang Jianping, Luo Hongxia, et al. Temporal and spatial distribution of maize drought in Southwest of China based on agricultural reference index for drought[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(2): 106-115. (in Chinese with English abstract)

    [12] 鮑艷松,嚴(yán)婧,閔錦忠,等. 基于溫度植被干旱指數(shù)的江蘇淮北地區(qū)農(nóng)業(yè)旱情監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(7):163-172.

    Bao Yansong, Yan Jing, Min Jinzhong, et al. Agricultural drought monitoring in north Jiangsu by using temperature vegetation dryness index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(7): 163-172. (in Chinese with English abstract)

    [13] 中華人民共和國國家標(biāo)準(zhǔn),氣象干旱等級:GB/T20481-2006[S]. 北京:中國標(biāo)準(zhǔn)出版社,2006.

    [14] Jiang Rengui, Xie Jiancang, He Hailong, et al. Use of four drought indices for evaluating drought characteristics under climate change in Shaanxi, China: 1951-2012[J]. Natural Hazards, 2015, 75(3): 2885-2903.

    [15] 胡毅鴻,李景保. 1951-2015 年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(7):107-115.

    Hu Yihong, Li Jingbao. Analysis on evolution of drought-flood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(7): 107-115. (in Chinese with English abstract)

    [16] 王媛媛,張勃. 基于標(biāo)準(zhǔn)化降水指數(shù)的近40 a隴東地區(qū)旱澇時空特征[J]. 自然資源學(xué)報,2012,27(12):2135-2144.

    Wang Yuanyuan, Zhang Bo. Anlysis of drought-flood spatial-temporal characteristics based on standard precipitation index in east region of Gansu in recent 40 years[J]. Journal of Natural Resources, 2012, 27(12): 2135-2144. (in Chinese with English abstract)

    [17] 張巧鳳,劉桂香,于紅博,等. 基于標(biāo)準(zhǔn)化降水指數(shù)的錫林郭勒盟干旱特征分析[J]. 自然災(zāi)害學(xué)報,2015,24(5):119-128.

    Zhang Qiaofeng, Liu Guixiang, Yu Hongbo, et al. Anlysis of drought characteristics in Xilingol League based on standardized precipitation index[J]. Journal of Natural Disasters, 2015, 24(5): 119-128. (in Chinese with English abstract)

    [18] 高西寧,徐慶喆,叢俊霞,等. 基于標(biāo)準(zhǔn)化降水指數(shù)的遼寧省近54年干旱時空規(guī)律分析[J]. 生態(tài)環(huán)境學(xué)報,2015,24(11):1851-1857.

    Gao Xining, Xu Qingzhe, Cong Junxia, et al. Temporal and spatial patterns of droughts based on standard precipitation index (SPI) in Liaoning Province in recent 54 a[J]. Ecology and Environmental Sciences, 2015, 24(11): 1851-1857. (in Chinese with English abstract)

    [19] 劉梅,魏加華,王峰,等. 渭河流域降水時空變化與干旱特征分析[J]. 南水北調(diào)與水利科技,2015,13(2):193-197.

    Liu Mei, Wei Jiahua, Wang Feng, et al. Spatial-temporal variations of precipitation and drought characteristics in Wei River Basin[J]. South-to-North Water Transfers and Water Science & Technology, 2015, 13(2): 193-197. (in Chinese with English abstract)

    [20] 孫智輝,王治亮,曹雪梅,等. 基于標(biāo)準(zhǔn)化降水指數(shù)的陜西省黃土高原地區(qū)1971-2010年干旱變化特征[J]. 中國沙漠,2013,33(5):1560-1567.

    Sun Zhihui, Wang Zhiliang, Cao Xuemei, et al. Characteristics of drought change in the loess plateau area of shaanxi based on the standardized precipitation index during 1971-2010[J]. Journal of Desert Research, 2013, 33(5): 1560-1567. (in Chinese with English abstract)

    [21] 任志遠(yuǎn),黃青,李晶. 陜西省生態(tài)安全及空間差異定量分析[J]. 地理學(xué)報, 2005, 60(4): 597-606.

    Ren Zhiyuan, Huang Qing, Li Jing. Quantitative analysis of dynamic change and spatial difference of the ecological safety: The case of Shaanxi province[J]. Acta Geograhica Sinica, 2005, 60(4): 597-606.

    [22] 翟祿新,馮起. 基于SPI的西北地區(qū)氣候干濕變化[J]. 自然資源學(xué)報,2011,26(5):847-857.

    Zhai Luxin, Feng Qi. Dryness/wetness climate variation based on standardized precipitation index in Northwest China[J]. Journal of Natural Resources, 2011, 26(5): 847-857. (in Chinese with English abstract)

    [23] 楊艷昭,張偉科,封志明,等. 干旱條件下南方紅壤丘陵地區(qū)水分平衡[J]. 農(nóng)業(yè)工程學(xué)報,2013,29(12):110-119.

    Yang Yanzhao, Zhang Weike, Feng Zhiming, et al. Water balance in red soil hilly regions of Southern China under drought climate situations[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(12): 110-119. (in Chinese with English abstract)

    [24] Angelidis P, Maris F, Kotsovinos N, et al. Computation of drought index SPI with alternative distribution functions[J]. Water Resources Management, 2012, 26(9): 2453-2473.

    [25] Moreira E, Coelho C A, Paulo A A, et al. SPI-based drought category prediction using loglinear models[J]. Journal of Hydrology, 2008, 354(1): 116-130.

    [26] 袁云,李棟梁,安迪. 基于標(biāo)準(zhǔn)化降水指數(shù)的中國冬季干旱分區(qū)及氣候特征[J]. 中國沙漠,2010,30(4):917-925.

    Yuan Yun, Li Dongliang, An Di. Winter aridity division in china based on standardized precipitation index and circulation characteristics[J]. Journal of Desert Research, 2010, 30(4): 917-925. (in Chinese with English abstract)

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

    Huang Wanhua, Yang Xiaoguang, Li Maosong, et al. Evolution characteristics of seasonal drought in the south of China during the past 58 years based on standardized precipitation index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(7): 50-59. (in Chinese with English abstract)

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

    Zhou Dan, Zhang Bo, Ren Peigui, et al. Analysis of drought characteristics of Shaanxi Province in recent 50 years based on standardized precipitation evapotranspiration index[J]. Journal of Natural Resources, 2014, 29(4): 677-688. (in Chinese with English abstract)

    [29] 王曉峰,張園,馮曉明,等. 基于游程理論和Copula函數(shù)的干旱特征分析及應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(10):206-214.

    Wang Xiaofeng, Zhang Yuan, Feng Xiaoming, et al. Analysis and application drought characteristics based on run theory and Copula function[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(10): 206-214. (in Chinese with English abstract)

    Analysis on spatiotemporal variability characteristics of drought in Shaanxi Province using standardized precipitation index

    Li Bin1,2, Xie Jiancang1, Hu Yanhua2, Jiang Rengui1※

    (1.,,710048,; 2.,710004,)

    Drought has become one of the most severe and frequent natural hazards, with major impacts on both human and natural systems. The large spatial coverage and long duration characteristics make it one of the most widespread and costliest hazards. With the rapid development of the economy and society, the acceleration of the urbanization process, and the sharp increase of the urban population, the water consumption of the various industries has increased rapidly, which puts forward higher requirements for the limited water resources. The increase of urban living and industrial water consumption occupies the agricultural water consumption, which makes the agricultural drought and water shortage more severe. Taking Shaanxi Province in the Northwest China, located in the 105°29′E-111°15′E and 31°42′N-39°35′N as the study area, the paper uses the standardized precipitation index (SPI), the drought frequency, ratio of number of stations with drought to total number of stations and drought intensity to analyze the spatial and temporal change patterns of drought at different time scales in Shaanxi Province based on the long-term daily observed precipitation data, provided by the China National Climate Center (CNCC). The study on the regional spatial and temporal variability of agriculture drought helps to promote the efficient use of the regional water resources and improve the regional drought resistance. The results showed that: 1) Generally, the drought in Shaanxi Province had an expanding tendency for the period of 1971-2013. The drought over the studied time period occurred frequently and the duration of drought was relatively long. The average drought frequency for the 34 meteorological stations was 70.0%. The frequencies of light drought and moderate drought were 23.7% and 23.3%, respectively, and the frequencies of severe drought and specially severe drought were relatively low, with the values of 14.2% and 8.8%, respectively. 2) The ratio of number of stations with drought to total number of stations had a seasonal change pattern. It had increasing trends in spring and winter, while had decreasing trends in summer and autumn. 3) The change pattern for drought intensity was not obvious. The drought intensity in spring and winter increased, while decreased in summer and autumn. We can conclude that the change patterns for ratio of number of stations with drought to total number of stations and drought intensity were similar. 4) Different regions had different drought change patterns, including provincial drought, regional drought and single station drought. The spatial distributions of drought frequency, ratio of number of stations with drought to total number of stations anddrought intensity also had difference.

    meteorology; drought; precipitation; standardized precipitation index; Shaanxi Province

    10.11975/j.issn.1002-6819.2017.17.015

    P429

    A

    1002-6819(2017)-17-0113-07

    2017-05-25

    2017-08-29

    國家重點研發(fā)計劃資助(2016YFC0401409);國家自然科學(xué)基金項目(51509201,51679188);中國博士后科學(xué)基金資助項目(2016M590964);陜西省高??茀f(xié)青年人才托舉計劃項目(20160217)。

    李 斌,陜西西安人,博士生,高級工程師,研究方向為水利信息化。西安 西安理工大學(xué)水利水電學(xué)院,710048。Email:157506520@qq.com

    姜仁貴,江西玉山人,博士,副教授,主要從事城市防災(zāi)減災(zāi)與應(yīng)急管理研究。西安 西安理工大學(xué)水利水電學(xué)院,710048。Email:jrengui@163.com

    猜你喜歡
    趨勢標(biāo)準(zhǔn)化
    標(biāo)準(zhǔn)化簡述
    趨勢
    企業(yè)標(biāo)準(zhǔn)化管理信息系統(tǒng)
    初秋唇妝趨勢
    Coco薇(2017年9期)2017-09-07 21:23:49
    標(biāo)準(zhǔn)化是綜合交通運輸?shù)谋U稀庾x《交通運輸標(biāo)準(zhǔn)化體系》
    中國公路(2017年9期)2017-07-25 13:26:38
    SPINEXPO?2017春夏流行趨勢
    “去編”大趨勢
    以標(biāo)準(zhǔn)化引領(lǐng)科技創(chuàng)新
    論汽車維修診斷標(biāo)準(zhǔn)化(上)
    趨勢
    汽車科技(2015年1期)2015-02-28 12:14:44
    性色av一级| av女优亚洲男人天堂| 亚洲 欧美一区二区三区| 国产亚洲最大av| 极品人妻少妇av视频| 国产精品成人在线| 狠狠婷婷综合久久久久久88av| av有码第一页| 国产一区二区激情短视频 | 啦啦啦在线观看免费高清www| 亚洲精品乱久久久久久| 少妇人妻久久综合中文| 久久性视频一级片| 亚洲一区二区三区欧美精品| 亚洲成人av在线免费| 欧美最新免费一区二区三区| 精品酒店卫生间| 高清视频免费观看一区二区| 久久狼人影院| 天天躁夜夜躁狠狠躁躁| 一级a爱视频在线免费观看| 无限看片的www在线观看| 午夜日韩欧美国产| 日韩欧美精品免费久久| 亚洲,一卡二卡三卡| 中国国产av一级| 成人三级做爰电影| 国产一区亚洲一区在线观看| 视频区图区小说| 国产精品国产三级专区第一集| 亚洲一码二码三码区别大吗| 国产av码专区亚洲av| 国产成人精品在线电影| 国产男女超爽视频在线观看| videosex国产| 亚洲成人av在线免费| 色94色欧美一区二区| 捣出白浆h1v1| 久久亚洲国产成人精品v| 91精品三级在线观看| 亚洲成人av在线免费| 叶爱在线成人免费视频播放| av.在线天堂| 看十八女毛片水多多多| 91精品三级在线观看| 高清视频免费观看一区二区| 80岁老熟妇乱子伦牲交| 国产免费又黄又爽又色| 波多野结衣av一区二区av| 成人午夜精彩视频在线观看| 一本一本久久a久久精品综合妖精| 亚洲成av片中文字幕在线观看| 80岁老熟妇乱子伦牲交| tube8黄色片| 黄色 视频免费看| av线在线观看网站| 伊人久久国产一区二区| 国产色婷婷99| 国产精品久久久av美女十八| 在线观看一区二区三区激情| 三上悠亚av全集在线观看| 亚洲熟女毛片儿| 婷婷色av中文字幕| 中文字幕制服av| 超碰97精品在线观看| 丰满饥渴人妻一区二区三| 午夜激情av网站| 亚洲自偷自拍图片 自拍| 2021少妇久久久久久久久久久| 日韩欧美精品免费久久| 嫩草影院入口| 高清不卡的av网站| 啦啦啦中文免费视频观看日本| 校园人妻丝袜中文字幕| 欧美日韩视频高清一区二区三区二| av又黄又爽大尺度在线免费看| 亚洲国产av影院在线观看| 99精品久久久久人妻精品| 国产精品成人在线| 十八禁人妻一区二区| 国产精品 国内视频| 亚洲精品aⅴ在线观看| av网站免费在线观看视频| 亚洲一区二区三区欧美精品| 中文字幕人妻丝袜制服| 久久毛片免费看一区二区三区| 久热这里只有精品99| 国产深夜福利视频在线观看| 国产亚洲av高清不卡| 午夜福利视频精品| 一二三四中文在线观看免费高清| 99久久综合免费| av不卡在线播放| 下体分泌物呈黄色| 一二三四在线观看免费中文在| 黄色视频在线播放观看不卡| 国产精品久久久久久精品电影小说| 晚上一个人看的免费电影| 亚洲欧美精品综合一区二区三区| 欧美乱码精品一区二区三区| 男男h啪啪无遮挡| 中文字幕精品免费在线观看视频| 最近手机中文字幕大全| 中文字幕精品免费在线观看视频| 国产精品久久久久成人av| 女人爽到高潮嗷嗷叫在线视频| 国产免费一区二区三区四区乱码| 欧美国产精品va在线观看不卡| 日日爽夜夜爽网站| 搡老乐熟女国产| 日本av手机在线免费观看| 丰满迷人的少妇在线观看| 青春草亚洲视频在线观看| 国产精品 欧美亚洲| 午夜激情久久久久久久| 欧美久久黑人一区二区| av天堂久久9| av在线播放精品| 好男人视频免费观看在线| 一级毛片我不卡| 99久国产av精品国产电影| 欧美日韩精品网址| 下体分泌物呈黄色| 黄色视频在线播放观看不卡| 老司机影院成人| 国产精品熟女久久久久浪| 精品少妇一区二区三区视频日本电影 | 国产毛片在线视频| 少妇被粗大的猛进出69影院| a级片在线免费高清观看视频| 1024香蕉在线观看| 日本黄色日本黄色录像| 亚洲av日韩精品久久久久久密 | 水蜜桃什么品种好| 成年人午夜在线观看视频| www日本在线高清视频| 777久久人妻少妇嫩草av网站| 九色亚洲精品在线播放| 九色亚洲精品在线播放| 亚洲熟女毛片儿| 午夜精品国产一区二区电影| 精品亚洲成国产av| 久久天堂一区二区三区四区| 成年人免费黄色播放视频| 国产精品av久久久久免费| 亚洲色图综合在线观看| 9热在线视频观看99| 久久久久久久大尺度免费视频| 丁香六月天网| 久久久欧美国产精品| 国产伦理片在线播放av一区| 一级毛片黄色毛片免费观看视频| 精品少妇黑人巨大在线播放| 国产在线视频一区二区| 看免费成人av毛片| 日韩,欧美,国产一区二区三区| 99re6热这里在线精品视频| 90打野战视频偷拍视频| 欧美激情极品国产一区二区三区| 成人免费观看视频高清| 亚洲精品久久久久久婷婷小说| 亚洲一级一片aⅴ在线观看| 啦啦啦在线免费观看视频4| 青春草国产在线视频| 成年动漫av网址| 免费在线观看完整版高清| 19禁男女啪啪无遮挡网站| 90打野战视频偷拍视频| 你懂的网址亚洲精品在线观看| 男女下面插进去视频免费观看| av在线观看视频网站免费| 视频区图区小说| 精品久久蜜臀av无| 亚洲 欧美一区二区三区| av卡一久久| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 成年av动漫网址| 波多野结衣av一区二区av| 热99久久久久精品小说推荐| 国产成人精品在线电影| 久久久久久久久久久久大奶| 亚洲熟女精品中文字幕| 成人国语在线视频| 亚洲欧洲精品一区二区精品久久久 | 日韩人妻精品一区2区三区| 免费高清在线观看视频在线观看| 又大又黄又爽视频免费| av又黄又爽大尺度在线免费看| 国产高清国产精品国产三级| 亚洲国产欧美网| 少妇 在线观看| 999久久久国产精品视频| 美女午夜性视频免费| 亚洲男人天堂网一区| 天堂中文最新版在线下载| 丝袜人妻中文字幕| 在线观看免费日韩欧美大片| 亚洲色图 男人天堂 中文字幕| e午夜精品久久久久久久| 天美传媒精品一区二区| 成年人午夜在线观看视频| 亚洲精品国产区一区二| 国产深夜福利视频在线观看| 男女无遮挡免费网站观看| 亚洲精品一二三| 国产一区二区三区综合在线观看| 欧美xxⅹ黑人| 别揉我奶头~嗯~啊~动态视频 | 日本91视频免费播放| 亚洲第一av免费看| 国产在线视频一区二区| 亚洲精品aⅴ在线观看| 日本wwww免费看| 成人影院久久| 女性被躁到高潮视频| 老熟女久久久| 看十八女毛片水多多多| 只有这里有精品99| 亚洲精品乱久久久久久| 中文字幕高清在线视频| 久热爱精品视频在线9| 国产精品麻豆人妻色哟哟久久| 精品少妇黑人巨大在线播放| 天天操日日干夜夜撸| 精品久久久久久电影网| 校园人妻丝袜中文字幕| 麻豆乱淫一区二区| 国产精品香港三级国产av潘金莲 | www.av在线官网国产| 亚洲人成电影观看| 亚洲成国产人片在线观看| 日韩成人av中文字幕在线观看| 久久人人爽人人片av| 精品视频人人做人人爽| 肉色欧美久久久久久久蜜桃| 久久这里只有精品19| 久久久久久久精品精品| 一边摸一边做爽爽视频免费| 久久狼人影院| 久久人妻熟女aⅴ| 人人妻人人澡人人看| 日日摸夜夜添夜夜爱| 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 久久99一区二区三区| 久久这里只有精品19| 啦啦啦中文免费视频观看日本| 亚洲三区欧美一区| 国产色婷婷99| 老司机影院毛片| 免费观看a级毛片全部| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 精品亚洲乱码少妇综合久久| 亚洲av日韩精品久久久久久密 | 热re99久久精品国产66热6| 秋霞伦理黄片| 国产成人精品无人区| av片东京热男人的天堂| 欧美日韩国产mv在线观看视频| 免费看av在线观看网站| 精品一品国产午夜福利视频| 中文字幕另类日韩欧美亚洲嫩草| e午夜精品久久久久久久| 我要看黄色一级片免费的| 欧美黑人精品巨大| 交换朋友夫妻互换小说| 大片电影免费在线观看免费| 欧美日韩国产mv在线观看视频| 国产乱人偷精品视频| 久久国产精品大桥未久av| 色吧在线观看| av在线播放精品| 亚洲精品中文字幕在线视频| 热99久久久久精品小说推荐| 国产成人精品无人区| 国产成人欧美在线观看 | 午夜福利在线免费观看网站| 波多野结衣一区麻豆| 久久热在线av| 一边摸一边做爽爽视频免费| 高清视频免费观看一区二区| 日本色播在线视频| 丁香六月欧美| 18禁动态无遮挡网站| 久热这里只有精品99| 人成视频在线观看免费观看| 久久久久国产一级毛片高清牌| 成年人午夜在线观看视频| 午夜福利一区二区在线看| 视频在线观看一区二区三区| 女人精品久久久久毛片| 在线观看免费高清a一片| 丝袜在线中文字幕| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人看| 日韩成人av中文字幕在线观看| 在线免费观看不下载黄p国产| 精品一品国产午夜福利视频| 欧美xxⅹ黑人| 狂野欧美激情性xxxx| av卡一久久| 国产欧美日韩一区二区三区在线| 人妻一区二区av| 久久99精品国语久久久| 亚洲av电影在线观看一区二区三区| 丝袜脚勾引网站| 一级爰片在线观看| 久久国产精品大桥未久av| 日本色播在线视频| 国产精品女同一区二区软件| 精品一区二区三区av网在线观看 | 婷婷色综合大香蕉| 欧美另类一区| 亚洲伊人色综图| av线在线观看网站| 99国产精品免费福利视频| 哪个播放器可以免费观看大片| 亚洲成人手机| 丝袜喷水一区| 欧美精品人与动牲交sv欧美| 亚洲精品国产色婷婷电影| 啦啦啦啦在线视频资源| 中文字幕另类日韩欧美亚洲嫩草| 丰满迷人的少妇在线观看| tube8黄色片| 一级毛片黄色毛片免费观看视频| 日韩一卡2卡3卡4卡2021年| 中文字幕高清在线视频| 韩国精品一区二区三区| 不卡视频在线观看欧美| 亚洲熟女精品中文字幕| 午夜免费鲁丝| 麻豆精品久久久久久蜜桃| 国产午夜精品一二区理论片| 亚洲激情五月婷婷啪啪| xxx大片免费视频| 日本欧美国产在线视频| 国产精品国产av在线观看| 免费观看a级毛片全部| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 电影成人av| 欧美久久黑人一区二区| 天天躁夜夜躁狠狠躁躁| 丝袜美腿诱惑在线| 午夜激情久久久久久久| 一二三四在线观看免费中文在| 国产乱人偷精品视频| 亚洲久久久国产精品| 久久久久久久久久久久大奶| 性少妇av在线| 最近中文字幕高清免费大全6| 国产片特级美女逼逼视频| 一级爰片在线观看| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 男女午夜视频在线观看| 人成视频在线观看免费观看| 多毛熟女@视频| 久久人人爽人人片av| 51午夜福利影视在线观看| 成人漫画全彩无遮挡| 亚洲av在线观看美女高潮| 久久精品国产亚洲av高清一级| 久久人人爽av亚洲精品天堂| 狠狠精品人妻久久久久久综合| 亚洲成人av在线免费| 最新在线观看一区二区三区 | 日韩人妻精品一区2区三区| 多毛熟女@视频| 少妇人妻精品综合一区二区| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 色播在线永久视频| 赤兔流量卡办理| 黄片无遮挡物在线观看| 热re99久久国产66热| 日本av手机在线免费观看| 侵犯人妻中文字幕一二三四区| 久久鲁丝午夜福利片| 最黄视频免费看| 国产精品久久久久成人av| 免费看不卡的av| 精品久久蜜臀av无| 女的被弄到高潮叫床怎么办| 国产福利在线免费观看视频| 久久99热这里只频精品6学生| 搡老岳熟女国产| 国产 一区精品| 亚洲精品日本国产第一区| 老汉色av国产亚洲站长工具| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 久久女婷五月综合色啪小说| 多毛熟女@视频| 纯流量卡能插随身wifi吗| 亚洲av日韩在线播放| 我的亚洲天堂| 搡老岳熟女国产| 哪个播放器可以免费观看大片| 热99国产精品久久久久久7| 亚洲精品美女久久久久99蜜臀 | 日韩制服骚丝袜av| 亚洲精品成人av观看孕妇| 午夜福利在线免费观看网站| 久久ye,这里只有精品| 国产野战对白在线观看| 亚洲精品,欧美精品| 亚洲精品国产av成人精品| 亚洲久久久国产精品| 各种免费的搞黄视频| 亚洲国产欧美日韩在线播放| 精品国产乱码久久久久久男人| 激情五月婷婷亚洲| 国产野战对白在线观看| 日韩一本色道免费dvd| 七月丁香在线播放| 韩国av在线不卡| 亚洲欧美一区二区三区黑人| 人人妻人人添人人爽欧美一区卜| h视频一区二区三区| av网站在线播放免费| 一级爰片在线观看| 最近最新中文字幕免费大全7| 久久人人爽人人片av| 久久久国产精品麻豆| 黑人猛操日本美女一级片| 天天躁夜夜躁狠狠躁躁| 视频在线观看一区二区三区| 一区二区三区乱码不卡18| 不卡视频在线观看欧美| 色综合欧美亚洲国产小说| 中文字幕制服av| 一级,二级,三级黄色视频| av国产久精品久网站免费入址| 久久婷婷青草| 久久久精品国产亚洲av高清涩受| 精品少妇一区二区三区视频日本电影 | 中文字幕人妻丝袜制服| 久久久久精品久久久久真实原创| 五月开心婷婷网| 欧美激情极品国产一区二区三区| 精品一区二区三卡| 十八禁人妻一区二区| 午夜免费鲁丝| 国产色婷婷99| 国产极品粉嫩免费观看在线| av在线观看视频网站免费| 国产亚洲欧美精品永久| 超碰成人久久| 波多野结衣av一区二区av| a级毛片在线看网站| 日本爱情动作片www.在线观看| 纯流量卡能插随身wifi吗| 一本一本久久a久久精品综合妖精| 精品国产乱码久久久久久男人| 又大又黄又爽视频免费| 欧美日韩视频高清一区二区三区二| a 毛片基地| 国产精品久久久av美女十八| 夫妻性生交免费视频一级片| 最新的欧美精品一区二区| av又黄又爽大尺度在线免费看| 国产免费一区二区三区四区乱码| 中文字幕亚洲精品专区| 不卡视频在线观看欧美| 亚洲美女黄色视频免费看| 久久热在线av| 亚洲熟女毛片儿| 久久久精品免费免费高清| 人成视频在线观看免费观看| 中文字幕精品免费在线观看视频| 久久国产精品大桥未久av| 亚洲国产最新在线播放| 亚洲av电影在线进入| 美女福利国产在线| 不卡视频在线观看欧美| 久久精品久久精品一区二区三区| 亚洲欧美成人精品一区二区| 久久久久人妻精品一区果冻| 亚洲精华国产精华液的使用体验| 各种免费的搞黄视频| 国产野战对白在线观看| 99久久精品国产亚洲精品| 亚洲图色成人| 久久久久精品久久久久真实原创| 国产精品久久久久久精品古装| 一区二区三区精品91| 建设人人有责人人尽责人人享有的| 亚洲精品日韩在线中文字幕| 丝袜喷水一区| 久久ye,这里只有精品| 久久狼人影院| 女人精品久久久久毛片| 午夜激情av网站| 久久女婷五月综合色啪小说| 中文字幕av电影在线播放| 亚洲精品一区蜜桃| 国产极品粉嫩免费观看在线| 久久久久久久国产电影| 在线观看免费日韩欧美大片| 国产日韩欧美视频二区| 婷婷色麻豆天堂久久| 国产精品99久久99久久久不卡 | 久久99热这里只频精品6学生| 亚洲精品美女久久久久99蜜臀 | 欧美乱码精品一区二区三区| 国产精品av久久久久免费| 午夜福利影视在线免费观看| 成年美女黄网站色视频大全免费| 99久国产av精品国产电影| 亚洲av日韩在线播放| 久久精品国产综合久久久| 久久久国产精品麻豆| 亚洲国产欧美网| 黑丝袜美女国产一区| av在线播放精品| 国产精品一区二区在线观看99| av国产久精品久网站免费入址| 久久国产精品大桥未久av| 2021少妇久久久久久久久久久| 男女下面插进去视频免费观看| 午夜激情av网站| 只有这里有精品99| 黄色 视频免费看| 日韩不卡一区二区三区视频在线| 国产在线一区二区三区精| 亚洲av成人精品一二三区| 在线观看人妻少妇| 一边亲一边摸免费视频| 午夜免费观看性视频| 搡老乐熟女国产| 亚洲成人国产一区在线观看 | 天天躁日日躁夜夜躁夜夜| 深夜精品福利| 国产1区2区3区精品| 欧美日本中文国产一区发布| 黄色毛片三级朝国网站| 丝袜美腿诱惑在线| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 成人午夜精彩视频在线观看| av福利片在线| 亚洲美女搞黄在线观看| 亚洲av男天堂| 亚洲国产看品久久| 日韩不卡一区二区三区视频在线| 国产日韩欧美亚洲二区| 777米奇影视久久| 少妇被粗大的猛进出69影院| 色视频在线一区二区三区| 精品酒店卫生间| 嫩草影院入口| 精品国产乱码久久久久久男人| 五月天丁香电影| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区精品视频观看| 久久免费观看电影| 丝袜脚勾引网站| 国产黄色免费在线视频| 另类精品久久| 日本午夜av视频| 国产日韩欧美视频二区| 亚洲精品视频女| 亚洲熟女毛片儿| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 久久精品亚洲av国产电影网| 久久人人爽人人片av| 国产精品成人在线| 久久精品久久久久久久性| 日韩av不卡免费在线播放| 国产黄色视频一区二区在线观看| a级毛片在线看网站| 亚洲三区欧美一区| 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 人人妻人人澡人人爽人人夜夜| av在线播放精品| 久久久国产精品麻豆| 嫩草影院入口| 欧美最新免费一区二区三区| 在现免费观看毛片| 天天操日日干夜夜撸| 一本久久精品| 亚洲av日韩在线播放| 少妇人妻精品综合一区二区| 亚洲熟女精品中文字幕| 如何舔出高潮| 极品人妻少妇av视频| 一级爰片在线观看| 国产激情久久老熟女| 尾随美女入室| 亚洲国产欧美网| 欧美精品人与动牲交sv欧美| 久久久久国产精品人妻一区二区| 亚洲精品久久久久久婷婷小说| 超碰97精品在线观看| 人人妻人人澡人人看| 国产精品久久久av美女十八| 精品一区在线观看国产| 久久av网站| 久久韩国三级中文字幕| 国产成人精品在线电影| 午夜91福利影院| 精品免费久久久久久久清纯 | 国产一区亚洲一区在线观看| 久久久欧美国产精品| 青春草亚洲视频在线观看| 女性生殖器流出的白浆| 日本欧美国产在线视频| 中文欧美无线码| 老熟女久久久|