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

    西南地區(qū)水稻洪澇災(zāi)害風(fēng)險(xiǎn)評估與區(qū)劃*

    2016-11-12 06:16:38楊建瑩霍治國張桂香汪天穎
    中國農(nóng)業(yè)氣象 2016年5期
    關(guān)鍵詞:孕災(zāi)抗災(zāi)洪澇

    楊建瑩,霍治國,2**,吳 立 ,張桂香,汪天穎

    ?

    西南地區(qū)水稻洪澇災(zāi)害風(fēng)險(xiǎn)評估與區(qū)劃*

    楊建瑩1,霍治國1,2**,吳 立1,張桂香1,汪天穎1

    (1. 中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,北京 100081;2. 南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,南京 210044)

    為全面評估水稻洪澇的綜合風(fēng)險(xiǎn),基于自然災(zāi)害系統(tǒng)理論和農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)評估方法,利用西南地區(qū)(重慶、四川、貴州和云南)193個(gè)氣象站1961-2012年逐日降水資料、396個(gè)縣(市)1981-2012年水稻產(chǎn)量、面積資料和17個(gè)農(nóng)氣站點(diǎn)水稻生育期數(shù)據(jù),以及西南地區(qū)數(shù)字高程(DEM)數(shù)據(jù),構(gòu)建區(qū)域水稻洪澇災(zāi)害致災(zāi)因子危險(xiǎn)性、承災(zāi)體暴露性、孕災(zāi)環(huán)境敏感性和區(qū)域抗災(zāi)能力指數(shù),以及綜合風(fēng)險(xiǎn)評價(jià)模型,對西南地區(qū)水稻洪澇進(jìn)行風(fēng)險(xiǎn)分析與區(qū)劃。結(jié)果表明:(1)水稻不同生育階段洪澇等級風(fēng)險(xiǎn)概率分布存在明顯的地區(qū)差異,洪澇危險(xiǎn)性表現(xiàn)為移栽分蘗期>拔節(jié)孕穗期>抽穗成熟期;全生育期高、次高危險(xiǎn)區(qū)主要分布于云南南部和東北部、貴州南部,以及四川的成都、眉山和德陽地區(qū)。(2)基于不同時(shí)間序列的水稻相對暴露率明顯波動,水稻生產(chǎn)承災(zāi)體高、次高暴露區(qū)主要集中在四川東北部和重慶地區(qū);孕災(zāi)環(huán)境高、次高敏感區(qū)主要位于云南北部、四川南部和貴州東南部地區(qū);水稻洪澇低抗災(zāi)能力區(qū)主要位于貴州。(3)西南地區(qū)水稻洪澇綜合風(fēng)險(xiǎn)呈由中部向四周遞增的趨勢,高、次高風(fēng)險(xiǎn)區(qū)主要位于貴州南部、云南南部和四川東北部地區(qū),低風(fēng)險(xiǎn)區(qū)位于重慶南部和云南北部地區(qū)。

    洪澇災(zāi)害;農(nóng)業(yè)氣象;風(fēng)險(xiǎn)分析;區(qū)劃;一季稻

    中國西南地區(qū)地形、地貌復(fù)雜,受季風(fēng)和青藏高原環(huán)流系統(tǒng)的影響,雨季降水過于集中,洪澇災(zāi)害發(fā)生頻率高、強(qiáng)度大[1]。水稻是西南地區(qū)主要糧食作物,常年種植面積近467萬hm2。由于水稻的生長期恰逢雨季,經(jīng)常會出現(xiàn)連續(xù)降雨或暴雨形成洪澇災(zāi)害,造成水稻受淹或沖毀,導(dǎo)致減產(chǎn)甚至絕收[2]。如1998年8月上旬,云南瀾滄地區(qū)發(fā)生強(qiáng)降水過程,水稻大面積被淹,減產(chǎn)20%~30%。2015年7月21-24日,四川盆地累計(jì)雨量100~250mm,局部地區(qū)雨量超過300mm,造成水稻大面積成災(zāi)。因此,開展西南地區(qū)水稻洪澇風(fēng)險(xiǎn)研究,有助于提升區(qū)域?yàn)?zāi)害風(fēng)險(xiǎn)管理和決策水平,減輕洪澇災(zāi)害給水稻生產(chǎn)造成的損失。

    洪澇災(zāi)害風(fēng)險(xiǎn)研究,主要集中在某一地區(qū)某一時(shí)間內(nèi)洪澇發(fā)生的可能性、活動程度、破壞損失及其對經(jīng)濟(jì)、社會和自然環(huán)境系統(tǒng)造成的影響和危害[3-4],多從致災(zāi)因子、承災(zāi)體、孕災(zāi)環(huán)境、防災(zāi)能力等單因素、集成多因素出發(fā),通過構(gòu)建指標(biāo)體系,進(jìn)行區(qū)域風(fēng)險(xiǎn)的定性和定量評估[5-7];研究對象涉及城市[8]、人口[9]、經(jīng)濟(jì)[10-11]、農(nóng)業(yè)[12]等領(lǐng)域。近幾年,針對作物的洪澇風(fēng)險(xiǎn)研究逐漸發(fā)展起來,如盛紹學(xué)等[13]基于災(zāi)害風(fēng)險(xiǎn)分析理論,針對江淮地區(qū)各縣小麥從澇漬脆弱性、自然氣候風(fēng)險(xiǎn)、災(zāi)損風(fēng)險(xiǎn)和抗災(zāi)能力4個(gè)方面對江淮地區(qū)小麥澇漬災(zāi)害風(fēng)險(xiǎn)進(jìn)行空間區(qū)域劃分;張蕾等[14]綜合苗期濕澇危險(xiǎn)性、孕災(zāi)環(huán)境敏感性、承災(zāi)體易損性和防災(zāi)減災(zāi)能力進(jìn)行瓜菜苗期濕澇綜合風(fēng)險(xiǎn)分析和區(qū)劃。目前,作物洪澇研究主要集中在冬小麥[15-16]、玉米[17]、蘋果[18]、油菜[19]等旱生作物上,針對水生作物的洪澇風(fēng)險(xiǎn)研究則甚少。

    水稻屬沼澤類作物,抗?jié)衬芰?qiáng),但若遭遇持續(xù)性降雨或強(qiáng)降雨,仍會造成大面積減產(chǎn)甚至絕收。關(guān)于水稻洪澇,眾多學(xué)者從水稻受澇特征[20]、抗?jié)衬芰21-22]以及水稻洪澇致災(zāi)指標(biāo)[23]等方面進(jìn)行了研究。但是,由于農(nóng)業(yè)生產(chǎn)系統(tǒng)的復(fù)雜性和相關(guān)數(shù)據(jù)獲得的局限性,結(jié)合農(nóng)業(yè)生產(chǎn)過程、針對作物生長特點(diǎn)的水稻洪澇風(fēng)險(xiǎn)研究還未見報(bào)道。

    因此,基于自然災(zāi)害系統(tǒng)理論[24]和農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)評估方法[25],本文以西南地區(qū)一季稻洪澇災(zāi)害為研究對象,針對水稻生產(chǎn)中容易遭受的洪澇災(zāi)害,構(gòu)建致災(zāi)因子危險(xiǎn)性、承災(zāi)體暴露性、孕災(zāi)環(huán)境敏感性(脆弱性)、抗災(zāi)能力指標(biāo),以及綜合風(fēng)險(xiǎn)評估模型,評估區(qū)域水稻洪澇風(fēng)險(xiǎn),以期為因地制宜開展水稻洪澇災(zāi)害風(fēng)險(xiǎn)管理、優(yōu)化生產(chǎn)布局提供科技支撐。

    1 資料與方法

    1.1 資料來源

    西南地區(qū)水稻種植區(qū)域界定依據(jù)國家氣象中心2012年編著的《農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)手冊》(非出版物),包括四川、重慶、云南、貴州四?。ㄊ校┑?96個(gè)縣(市),區(qū)內(nèi)193個(gè)氣象站點(diǎn)1961-2012年的逐日降水量資料來源于國家氣象信息中心,水稻生育期資料來源于西南地區(qū)17個(gè)農(nóng)業(yè)氣象觀測站的觀測數(shù)據(jù),包括1991-2012年水稻播種、移栽、拔節(jié)、抽穗和成熟日期;縣級尺度的農(nóng)業(yè)生產(chǎn)數(shù)據(jù)來源于統(tǒng)計(jì)數(shù)據(jù),包括1981-2012年西南地區(qū)396個(gè)縣(市)主要作物的播種面積、產(chǎn)量等;DEM數(shù)據(jù)來源于美國地質(zhì)勘探局(USGS)生產(chǎn)的全球30″數(shù)字高程模型(GTOPO 30),其空間分辨率為30″(約1km)。區(qū)域范圍及各站點(diǎn)分布見圖1。

    根據(jù)水稻發(fā)育期資料,將研究區(qū)水稻發(fā)育期分成3個(gè)階段,即移栽分蘗期、拔節(jié)孕穗期和抽穗成熟期。水稻在不同?。ㄊ校┑陌l(fā)育時(shí)期差別較大,這主要與氣候特點(diǎn)和水稻的品種熟型有關(guān)。本文中各站點(diǎn)的生育期數(shù)據(jù)采用分省統(tǒng)計(jì)的水稻各發(fā)育階段的開始和結(jié)束日期,4?。ㄊ校┧景l(fā)育期情況見表1。

    表1 西南地區(qū)各?。ㄊ校┧旧?/p>

    Note: E-,M-,L- is the first,middle and last ten-day of a month, respectively.

    1.2 構(gòu)建評估指標(biāo)體系

    1.2.1 致災(zāi)因子危險(xiǎn)性指數(shù)

    (1)水稻洪澇等級評價(jià)指標(biāo)

    基于前期研究成果[23],根據(jù)水稻不同生育階段降水量和降水持續(xù)時(shí)間分別評定移栽分蘗期、拔節(jié)孕穗期和抽穗成熟期的洪澇等級,標(biāo)準(zhǔn)見表2。其中,輕度洪澇的受災(zāi)表現(xiàn)為稻田受淹、水稻受災(zāi)、部分成災(zāi);中度洪澇的受災(zāi)表現(xiàn)為淹沒、水稻成災(zāi)、部分水稻絕收;重度洪澇的受災(zāi)表現(xiàn)為水稻絕收、無收、滅產(chǎn)。

    表2 水稻各生育階段洪澇等級評價(jià)指標(biāo)

    注:R為過程雨量(mm)。

    Note: R is rainfall amount of a rainfall event.

    (2)水稻洪澇危險(xiǎn)性指數(shù)

    水稻洪澇致災(zāi)因子危險(xiǎn)性由洪澇發(fā)生程度(強(qiáng)度)和頻次(概率)決定。一般洪澇強(qiáng)度越大,頻次越高,洪澇災(zāi)害所造成的水稻損失越嚴(yán)重,洪澇災(zāi)害的風(fēng)險(xiǎn)也越大。某地水稻整個(gè)生育期洪澇致災(zāi)風(fēng)險(xiǎn)指數(shù)(M)[12]計(jì)算式為

    式中,Qij為第j個(gè)生育階段第i個(gè)等級洪澇災(zāi)害的等級即輕度、中度和重度,分別取1、2和3,Pij為第j個(gè)生育階段第i個(gè)等級洪澇災(zāi)害發(fā)生的概率。n為洪澇的等級,n=3;m為水稻的3個(gè)發(fā)育時(shí)期,即移栽分蘗期、拔節(jié)孕穗期和抽穗成熟期,m=3。洪澇發(fā)生概率為依據(jù)洪澇等級評價(jià)指標(biāo)計(jì)算得到的1961-2012年各等級洪澇的風(fēng)險(xiǎn)概率,計(jì)算方法采用信息擴(kuò)散法[14],通過線性分配得到控制區(qū)間的原始信息分布,建立概率分布模型,計(jì)算得到不同等級洪澇風(fēng)險(xiǎn)概率。

    1.2.2 水稻暴露性指數(shù)

    水稻空間分布以及種植密度等均可反映水稻生產(chǎn)對洪澇災(zāi)害的暴露性。水稻暴露在洪澇災(zāi)害中的密度越大,受災(zāi)的風(fēng)險(xiǎn)也就越大。水稻相對暴露率(AR)在一定程度上可以進(jìn)行區(qū)域單元間水稻暴露性的比較。近年來,西南地區(qū)水稻種植面積整體呈下降趨勢,但區(qū)域間水稻種植密度波動存在明顯差異。通常情況下,基于長時(shí)間序列資料(1981-2012年)計(jì)算得到的水稻相對暴露率,對近年來水稻暴露風(fēng)險(xiǎn)的表達(dá)較弱;若以近5a(2008-2012年)和近3a(2010-2012年)水稻相對暴露率作為評估水稻暴露風(fēng)險(xiǎn)的依據(jù),所得結(jié)果則忽略了區(qū)域水稻長期種植的穩(wěn)定性。為更好地反映水稻暴露特征在時(shí)間序列上的區(qū)域穩(wěn)定性和變異特征,在水稻暴露性指數(shù)構(gòu)建過程中,分別對1981-2012年(32a)、1993-2012年(近20a)、2003-2012年(近10a)、2008-2012年(近5a)和2010-2012年(近3a)水稻相對暴露率進(jìn)行計(jì)算,以各時(shí)段水稻相對暴露率的加權(quán)得到水稻暴露性指數(shù)。某地水稻暴露性指數(shù)的計(jì)算式為

    式中,E為水稻暴露性指數(shù);ARi為某地每個(gè)計(jì)算時(shí)段的水稻相對暴露率;wi為每個(gè)時(shí)段相應(yīng)的權(quán)重,利用層次分析法[26]綜合確定。

    每個(gè)時(shí)段相對暴露率[27]的計(jì)算式為

    式中,ai為某縣第i年水稻播種面積(hm2),zi為該縣幅員面積(hm2),Ai為該年區(qū)域水稻播種面積(hm2),Zi為區(qū)域幅員面積,n為計(jì)算時(shí)段的年數(shù)。

    1.2.3 孕災(zāi)環(huán)境敏感性指數(shù)

    西南地區(qū)地形地貌極為復(fù)雜,境內(nèi)高原、山地、丘陵、平原、河谷均有分布,還有廣泛分布的喀斯特地貌。從洪澇形成的背景和機(jī)理分析,地形因子是洪澇災(zāi)害形成的主要孕災(zāi)因子。在地勢低洼平坦、地形閉塞的地區(qū),積水難以迅速宣泄,易于暴雨洪澇災(zāi)害的形成[28]。本文中西南地區(qū)地形因子的影響通過高程及高程標(biāo)準(zhǔn)差的不同組合來體現(xiàn)。其中,高程標(biāo)準(zhǔn)差表征地形的起伏程度,利用DEM數(shù)據(jù),采用1000m×1000m的網(wǎng)格大小,計(jì)算數(shù)據(jù)中某一格點(diǎn)的高程與周圍8個(gè)格點(diǎn)的高程標(biāo)準(zhǔn)差。依據(jù)柵格單元占區(qū)域總柵格數(shù)的百分比,將高程及高程標(biāo)準(zhǔn)差劃分為5個(gè)等級。通過專家打分法對高程及高程標(biāo)準(zhǔn)差進(jìn)行賦值[29](表3)。

    表3 高程和高程標(biāo)準(zhǔn)差組合的賦值標(biāo)準(zhǔn)

    1.2.4 抗災(zāi)能力指數(shù)

    水稻抗災(zāi)能力,從以下兩個(gè)方面進(jìn)行評估:(1)水稻產(chǎn)量水平。單產(chǎn)在一定程度上可以代表西南地區(qū)水稻生產(chǎn)的水平,亦可代表這一區(qū)域的農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平[30];(2)減產(chǎn)率的波動程度。水稻減產(chǎn)率的波動大小能體現(xiàn)水稻對外來干擾的適應(yīng)能力,波動大的區(qū)域水稻抗外來干擾能力較脆弱,波動小的區(qū)域水稻抗外來干擾能力較強(qiáng)。

    (1)水稻產(chǎn)量水平采用相對單產(chǎn)系數(shù)表示,即

    式中,Yr為某縣相對產(chǎn)量系數(shù),yi為第i年該縣水稻單產(chǎn)值,yi為第i年西南地區(qū)水稻單產(chǎn)平均值,n為總年數(shù)。

    (2)減產(chǎn)率的波動程度采用減產(chǎn)率變異系數(shù)表示,即

    式中,Ri為第i年某縣水稻減產(chǎn)率,R為該縣1981-2012年水稻減產(chǎn)率的平均值。將水稻減產(chǎn)率定義為當(dāng)年水稻實(shí)際產(chǎn)量與趨勢產(chǎn)量的比值。當(dāng)減產(chǎn)率為正值時(shí),則表示增產(chǎn);若減產(chǎn)率為負(fù)值,則表示減產(chǎn),計(jì)算方法參照《小麥干旱災(zāi)害等級》標(biāo)準(zhǔn)[31]。

    構(gòu)建的抗災(zāi)能力指數(shù)(F)可表示為

    (6)

    1.2.5 水稻洪澇風(fēng)險(xiǎn)指數(shù)

    為了消除不同因子量綱的影響,構(gòu)建水稻洪澇災(zāi)害綜合風(fēng)險(xiǎn)指數(shù)之前,根據(jù)數(shù)據(jù)極差標(biāo)準(zhǔn)化方法,將洪澇發(fā)生風(fēng)險(xiǎn)指數(shù)(M)和水稻承災(zāi)暴露性指數(shù)(E)、敏感性指數(shù)(S)和抗災(zāi)能力指數(shù)(F)極差標(biāo)準(zhǔn)化。為了標(biāo)準(zhǔn)化后數(shù)據(jù)非零,改進(jìn)的計(jì)算式[32]為

    式中,Xs為極差標(biāo)準(zhǔn)化后的數(shù)據(jù);Xj為原始數(shù)據(jù);Xmin為原始數(shù)據(jù)序列中的最小值;Xmax原始數(shù)據(jù)序列中的最大值。

    由于在構(gòu)成災(zāi)害風(fēng)險(xiǎn)的4個(gè)要素中,災(zāi)害的危險(xiǎn)性、承載體暴露性、孕災(zāi)環(huán)境敏感性和災(zāi)害風(fēng)險(xiǎn)生成的作用方向一致,而抗災(zāi)能力與災(zāi)害風(fēng)險(xiǎn)生成的作用方向相反,即某區(qū)域的抗災(zāi)能力越強(qiáng),災(zāi)害危險(xiǎn)性、承載體暴露性和孕災(zāi)環(huán)境敏感性對災(zāi)害生成的作用力就會受到一定抑制,從而減少災(zāi)害的風(fēng)險(xiǎn)度[27]。因此,在實(shí)際應(yīng)用時(shí),應(yīng)考慮到抗災(zāi)能力的反向作用力。由標(biāo)準(zhǔn)化后的洪澇災(zāi)害危險(xiǎn)性指數(shù)(Mjs)、承災(zāi)體暴露性指數(shù)(Ejs)、孕災(zāi)環(huán)境敏感性指數(shù)(Sjs)和抗災(zāi)能力(Fjs)共同構(gòu)建的西南地區(qū)水稻洪澇災(zāi)害綜合風(fēng)險(xiǎn)指數(shù)(FI)為

    式中,Mjs為標(biāo)準(zhǔn)化的j站洪澇風(fēng)險(xiǎn)指數(shù);Ejs為標(biāo)準(zhǔn)化的j站水稻暴露性風(fēng)險(xiǎn)指數(shù);Sjs為標(biāo)準(zhǔn)化的j區(qū)敏感性指數(shù);Fjs為標(biāo)準(zhǔn)化的j區(qū)抗災(zāi)能力指數(shù)。

    2 結(jié)果與分析

    2.1 水稻各生育階段洪澇致災(zāi)因子危險(xiǎn)性評估

    2.1.1 不同等級洪澇概率分布

    頸椎按摩操主要是通過主動活動鍛煉,達(dá)到疏通血脈和調(diào)暢氣機(jī)的目的,配合其他治療重新恢復(fù)頸椎的活動調(diào)節(jié)功能,從而達(dá)到消除緩解頸椎病的臨床癥狀,在治療康復(fù)頸椎病上有其應(yīng)用的理論依據(jù)。頸椎按摩操同時(shí)具有護(hù)士易于指導(dǎo)督促患者正確練習(xí)的特點(diǎn),且易學(xué)易練,是值得推廣的康復(fù)鍛煉方法。

    基于水稻洪澇等級評價(jià)指標(biāo),采用信息擴(kuò)散方法,計(jì)算1961-2012年西南地區(qū)193個(gè)站點(diǎn)水稻輕、中、重度洪澇的風(fēng)險(xiǎn)概率,采用IDW插值法[33]繪制水稻洪澇危險(xiǎn)性空間分布圖(圖2)。由圖2可見,移栽分蘗期洪澇危險(xiǎn)性風(fēng)險(xiǎn)概率呈東高西低的分布,高發(fā)區(qū)域主要位于貴州南部和云南南部地區(qū),該地區(qū)1961-2012年輕、中和重度洪澇風(fēng)險(xiǎn)概率分別在0.16、0.1和0.08以上,其中云南南部的金平、龍陵、綠春、羅平、西盟,貴州的都勻、興仁和六枝地區(qū),輕、中和重度洪澇風(fēng)險(xiǎn)概率分別達(dá)到0.35、0.2和0.15以上,是水稻移栽分蘗期洪澇的高發(fā)區(qū)域。水稻拔節(jié)孕穗期洪澇風(fēng)險(xiǎn)概率較移栽分蘗期降低,高值區(qū)域主要位于云南南部地區(qū),輕、中和重度洪澇風(fēng)險(xiǎn)概率分別達(dá)到0.16、0.12和0.12以上。由于水稻拔節(jié)孕穗期正值汛期,四川中部和東部地區(qū),洪澇風(fēng)險(xiǎn)概率較移栽分蘗期有所增加,尤其是峨眉、名山、洪雅、安縣、夾江以及綿竹等地,輕、中和重度洪澇風(fēng)險(xiǎn)概率分別在0.08、0.04和0.02以上。抽穗成熟期洪澇風(fēng)險(xiǎn)概率空間特征與拔節(jié)孕穗期基本一致,風(fēng)險(xiǎn)概率高值區(qū)域主要位于云南南部和四川中部地區(qū),其中云南南部地區(qū)水稻重度洪澇的風(fēng)險(xiǎn)概率高于輕度和中度洪澇,達(dá)到0.16以上??傮w來說,水稻洪澇風(fēng)險(xiǎn)概率表現(xiàn)為移栽分蘗期>拔節(jié)孕穗期>抽穗成熟期。移栽分蘗期洪澇和拔節(jié)孕穗期以輕度洪澇為主,而抽穗成熟期易發(fā)生水稻重度洪澇。

    2.1.2 洪澇危險(xiǎn)性評估與分布

    利用各生育階段洪澇發(fā)生等級和概率,根據(jù)危險(xiǎn)性評價(jià)模型,計(jì)算得到193個(gè)氣象站點(diǎn)水稻洪澇危險(xiǎn)性指數(shù)。通過IDW插值方法得到水稻洪澇危險(xiǎn)性指數(shù)空間分布,采用百分位法,將水稻洪澇危險(xiǎn)性指數(shù)劃分成4個(gè)等級,對應(yīng)得到水稻洪澇低、中、次高和高危險(xiǎn)性評價(jià)指標(biāo)的分級臨界值(表4)及水稻洪澇危險(xiǎn)性分布(圖3)。由圖3可見,移栽分蘗期的高風(fēng)險(xiǎn)區(qū)主要位于云南南部和貴州西南部地區(qū),危險(xiǎn)性指數(shù)≥1.25;次高危險(xiǎn)區(qū)包括云南東北部和貴州南部,危險(xiǎn)性指數(shù)介于0.27~1.25。云南南部和四川東部地區(qū)水稻拔節(jié)孕穗期處于洪澇高或次高危險(xiǎn)區(qū),水稻洪澇危險(xiǎn)性指數(shù)≥0.53;抽穗成熟期高、次高危險(xiǎn)性包括云南南部和四川東北部地區(qū),危險(xiǎn)性指數(shù)均在0.55以上。

    表4 西南地區(qū)水稻洪澇危險(xiǎn)性評價(jià)指標(biāo)(M)的分級界限值

    西南地區(qū)受季風(fēng)影響顯著,汛期降水主要受西南季風(fēng)和東亞季風(fēng)的交替影響,雨季一般為5-10月,該時(shí)段降水集中,大部分地區(qū)的雨量占年降水總量的85%以上[34]。水稻生長初期降水量從南向北減少,云南南部地區(qū)降水量明顯較其它區(qū)域偏高,容易造成水稻生長前期發(fā)生洪澇危害。隨著水稻生長,四川東北部逐漸進(jìn)入主汛期,該地區(qū)同期平均降水量超過500mm[35],極易造成水稻生長后期的洪澇災(zāi)害。從整個(gè)水稻生長階段來看,洪澇災(zāi)害的高、次高危險(xiǎn)區(qū)主要位于云南南部地區(qū),全生育期危險(xiǎn)性指數(shù)≥1.52;中危險(xiǎn)區(qū)主要包括云南東北部、貴州南部,以及四川的成都、眉山和德陽地區(qū),危險(xiǎn)性指數(shù)在0.41~1.52。

    2.2 水稻生產(chǎn)承災(zāi)體暴露性評估

    暴露性反映了水稻種植面積在幅員面積中所占的比重,比重越高意味著暴露于洪澇危險(xiǎn)中的風(fēng)險(xiǎn)越高。根據(jù)前文對水稻相對暴露率的定義,計(jì)算得到西南地區(qū)396個(gè)縣(市)的水稻相對暴露率,并賦值到這396縣(市)的地理中心點(diǎn),通過IDW插值方法得到區(qū)域1981-2012年(32a)、1993-2012年(近20a)、2003-2012年(近10a)、2008-2012年(近5a)和2010-2012年(近3a)的水稻相對暴露率空間分布(圖4)。由圖4可見,四川北部以及重慶地區(qū)水稻種植比重較大,大部分地區(qū)水稻相對暴露率在1.6以上;云南以及貴州南部水稻相對暴露率不足0.8。貴州地區(qū)相對暴露率在時(shí)序上存在明顯波動,水稻相對暴露率>1.6的區(qū)域在2003-2012年(近10a)、2008-2012年(近5a)和2010-2012年(近3a),明顯高于1993-2012年(近20a)和1981-2012年全時(shí)間序列;云南南部地區(qū)水稻種植在2008-2012年(近5a)和2010-2012年(近3a)減少,相對暴露率<0.8的區(qū)域有所擴(kuò)大,>0.8的區(qū)域明顯減少。

    采用層次分析方法,得到各時(shí)間段水稻相對暴露率的權(quán)重。1981-2012年(32a)、1993-2012年(近20a)、2003-2012年(近10a)、2008-2012年(近5a)和2010-2012年(近3a)的水稻暴露率權(quán)重分別為0.29、0.19、0.11、0.20和0.21。計(jì)算得到西南地區(qū)水稻暴露性指數(shù),采用百分位法進(jìn)行水稻暴露性低、中、次高和高4個(gè)等級評價(jià)指標(biāo)的界限值率定(表5)并制圖(圖5)。結(jié)果表明:1961-2012年西南地區(qū)水稻高、次高暴露區(qū)主要集中在四川東北部和重慶地區(qū),該區(qū)域水稻種植密度高于西南地區(qū)水稻種植密度平均水平,暴露性指數(shù)≥0.72,意味著該區(qū)域暴露于洪澇危險(xiǎn)中的風(fēng)險(xiǎn)也較高;中度暴露區(qū)主要分布于云南西南部和東部地區(qū),該區(qū)域暴露性指數(shù)在0.45~0.72,水稻種植密度與西南地區(qū)水稻種植密度平均水平相當(dāng),水稻暴露于洪澇危險(xiǎn)中的風(fēng)險(xiǎn)也介于中間;低暴露區(qū)包括云南北部、四川南部和貴州西南部地區(qū)(<0.45),該區(qū)域水稻種植密度低于西南地區(qū)水稻種植密度平均水平,即暴露于洪澇危險(xiǎn)中的風(fēng)險(xiǎn)相對較低。

    表5 西南地區(qū)水稻暴露性評價(jià)指標(biāo)(E)的分級界限值

    2.3 水稻生產(chǎn)孕災(zāi)環(huán)境敏感性評估

    利用1000m×1000m DEM數(shù)據(jù),通過柵格計(jì)算與賦值,采用百分位分級方法,得到西南地區(qū)孕災(zāi)環(huán)境敏感性評價(jià)指標(biāo)(S)的分級界限值(表6)和孕災(zāi)環(huán)境敏感性分布(圖6)。由圖6可見,洪澇孕災(zāi)環(huán)境敏感性的空間分布表現(xiàn)為中間高四周低,孕災(zāi)敏感性由四周向中部地區(qū)逐漸增加。西南地區(qū)中部地區(qū),主要包括云南北部、四川南部和貴州東南部地區(qū),是孕災(zāi)環(huán)境高敏感區(qū)域和次高敏感區(qū),敏感性指數(shù)≥0.80,該區(qū)域地勢低洼,高差小,當(dāng)發(fā)生強(qiáng)降水過程時(shí),水澇不易排出;中度敏感區(qū)主要分布于云南南部、貴州中部地區(qū),敏感性指數(shù)在0.70~0.80,該區(qū)域地勢和高差均處中等水平,孕災(zāi)環(huán)境敏感性與整個(gè)區(qū)域敏感性相當(dāng);低敏感性區(qū)域集中在四川東北部、重慶和貴州的山地丘陵區(qū),S均在0.70以下,該區(qū)域海拔較高且地面起伏程度較大,易于水澇的排出。

    表6 西南地區(qū)孕災(zāi)環(huán)境敏感性評價(jià)指標(biāo)(S)的分級界限值

    2.4 水稻生產(chǎn)抗災(zāi)能力評估

    西南地區(qū)水稻抗災(zāi)能力通過當(dāng)?shù)厮井a(chǎn)量高低水平以及水稻減產(chǎn)率的波動程度綜合表現(xiàn)。根據(jù)前文對相對單產(chǎn)系數(shù)(Yr)和減產(chǎn)率變異系數(shù)(Ys)的定義,分別計(jì)算得到西南地區(qū)396個(gè)縣(市)的Yr和Ys并賦值于縣(市)的地理中心點(diǎn),采用IDW插值方法,得到西南地區(qū)Yr和Ys的空間分布(圖7)。從圖7a可知,西南地區(qū)水稻產(chǎn)量水平的高值區(qū)域分布于四川以及云南北部地區(qū),相對單產(chǎn)系數(shù)Yr均在1以上;低產(chǎn)區(qū)位于貴州西部以及云南南部,相對單產(chǎn)系數(shù)<0.7。從水稻減產(chǎn)率的波動程度(圖7b)看,貴州以及云南東南部的紅河、文山地區(qū),水稻減產(chǎn)率變異較大,減產(chǎn)率變異系數(shù)>1;四川中北部,以及云南地區(qū)水稻減產(chǎn)率的波動較小,變異系數(shù)<0.25。綜合來看(表7和圖8),西南地區(qū)低抗災(zāi)能力區(qū)主要位于貴州,抗災(zāi)能力指數(shù)F<0.95,該區(qū)域水稻產(chǎn)量相對較低,且減產(chǎn)率波動較大,水稻生產(chǎn)力整體水平低;高、次高抗災(zāi)能力區(qū)主要位于四川和云南北部地區(qū),抗災(zāi)能力指數(shù)≥2.74,該區(qū)域水稻產(chǎn)量高于西南地區(qū)水稻產(chǎn)量平均水平,且產(chǎn)量穩(wěn)定,意味著當(dāng)?shù)厮酒贩N的適應(yīng)性強(qiáng),水稻生產(chǎn)力整體水平較高;其它區(qū)域?qū)僦卸瓤篂?zāi)能力區(qū),抗災(zāi)能力指數(shù)介于0.95~2.74,水稻生產(chǎn)水平居中。

    表7 西南地區(qū)抗災(zāi)能力評價(jià)指標(biāo)(F)的分級界限值

    2.5 西南地區(qū)水稻生育期洪澇災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃

    中度風(fēng)險(xiǎn)區(qū)主要分布于貴州南部和四川東北部地區(qū),綜合風(fēng)險(xiǎn)指數(shù)介于0.08~0.14。該區(qū)域水稻全生育期大部處于中度洪澇危險(xiǎn)區(qū),部分區(qū)域處于次低洪澇危險(xiǎn)區(qū),水稻種植面積比例、產(chǎn)量波動大小和抗災(zāi)能力大部高于或與西南地區(qū)平均水平相當(dāng)。

    低風(fēng)險(xiǎn)區(qū)主要位于重慶、貴州和云南北部地區(qū),該區(qū)域水稻洪澇綜合風(fēng)險(xiǎn)指數(shù)<0.08。盡管貴州和云南北部地區(qū)地勢低洼且平坦,孕災(zāi)環(huán)境敏感性高,易于洪澇的形成,但該區(qū)域水稻生產(chǎn)遭受洪澇危害的機(jī)率較低,且水稻種植較少,因此是水稻洪澇災(zāi)害的低風(fēng)險(xiǎn)區(qū)域。重慶地區(qū)水稻種植廣泛,但該區(qū)域多山地丘陵,易于水澇及時(shí)排出,盡管抗災(zāi)能力較低,但該區(qū)域確為西南地區(qū)水稻生產(chǎn)遭受洪澇危害的低風(fēng)險(xiǎn)區(qū)。

    3 結(jié)論與討論

    3.1 討論

    基于自然災(zāi)害系統(tǒng)理論和農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)評估方法,以西南地區(qū)一季稻洪澇災(zāi)害為研究對象,綜合考慮水稻洪澇風(fēng)險(xiǎn)形成的危險(xiǎn)性、暴露性、敏感性、抗災(zāi)能力4個(gè)影響因子,構(gòu)建具有區(qū)域特征的指標(biāo)體系評估模型。王春乙等[27,30]認(rèn)為,危險(xiǎn)性評價(jià)是風(fēng)險(xiǎn)評估的主體,災(zāi)害危險(xiǎn)性分布直接影響災(zāi)害風(fēng)險(xiǎn)的大小和空間分布。對于水生作物,短時(shí)大暴雨是洪澇災(zāi)害的主要誘因。針對水稻不同生育階段的抗?jié)程卣?,本文在前期研究的基礎(chǔ)上,選擇能表征水稻洪澇短時(shí)暴雨指標(biāo)進(jìn)行危險(xiǎn)性評估。危險(xiǎn)性評價(jià)結(jié)果與吳立等[12]以整個(gè)農(nóng)業(yè)為承災(zāi)體的洪澇災(zāi)害分析結(jié)果基本相同,但局部地區(qū)存在差異。這與承災(zāi)體種類、洪澇指標(biāo)的選取密切相關(guān)。水稻抗?jié)衬芰?qiáng),采用的降水強(qiáng)度指標(biāo),略高于以整個(gè)農(nóng)業(yè)為承災(zāi)體的降水強(qiáng)度指標(biāo)。最終得到的西南地區(qū)水稻洪澇綜合風(fēng)險(xiǎn)結(jié)果,與洪澇危險(xiǎn)性空間分布圖大體一致。

    表8 西南地區(qū)水稻洪澇綜合風(fēng)險(xiǎn)評價(jià)指標(biāo)(FI)的分級界限值

    考慮水稻暴露特征在時(shí)間序列上的區(qū)域穩(wěn)定性和變異特征,在水稻暴露性指數(shù)構(gòu)建過程中,引入權(quán)重系數(shù),構(gòu)建基于不同時(shí)間序列水稻相對暴露率的暴露性指數(shù),有助于減小水稻暴露性評估中由于時(shí)序數(shù)據(jù)較長或較短對評估結(jié)果造成的影響。針對西南地區(qū)地形地貌復(fù)雜的特點(diǎn),在孕災(zāi)環(huán)境敏感性指數(shù)構(gòu)建過程中,選取地理高程和地表起伏程度指標(biāo),從降水產(chǎn)流與匯流這一主導(dǎo)因素反映地形地貌對降水集聚的影響。單產(chǎn)水平和減產(chǎn)率的波動情況從不同角度反映區(qū)域抗災(zāi)能力,但都存在不全面性,如某縣的水稻單產(chǎn)水平高,但水稻生產(chǎn)波動性大,并不能說明該縣的抗災(zāi)能力強(qiáng)。因此,本文采用了單產(chǎn)水平與減產(chǎn)率波動兩個(gè)指標(biāo)來構(gòu)建抗災(zāi)能力指數(shù)。承災(zāi)體暴露性、孕災(zāi)環(huán)境敏感性和抗災(zāi)能力間接影響洪澇風(fēng)險(xiǎn)的形成。例如,從全生育期水稻洪澇危險(xiǎn)性來看,貴州大部分地區(qū)處于中危險(xiǎn)區(qū),孕災(zāi)敏感性和水稻暴露性均處中等水平,但是,該區(qū)域受陡坡、瘠薄等問題的影響嚴(yán)重,土地質(zhì)量和土壤的肥沃程度相對低下[36],水稻單產(chǎn)水平低。加之該區(qū)域農(nóng)業(yè)生產(chǎn)設(shè)施不完善,抗災(zāi)能力弱,當(dāng)遭遇強(qiáng)降水天氣過程,受洪澇的影響較大。通過綜合風(fēng)險(xiǎn)分析,最終被歸為高、次高洪澇風(fēng)險(xiǎn)區(qū)。為驗(yàn)證水稻洪澇綜合風(fēng)險(xiǎn)模型的合理性,將本文水稻洪澇綜合風(fēng)險(xiǎn)評估結(jié)果與《中國氣象災(zāi)害大典(四川、重慶、貴州和云南)》[31,37-39]中記載的歷史水稻洪澇情況進(jìn)行對比分析。350條關(guān)于水稻洪澇的記錄中,有321條記錄的水稻洪澇發(fā)生地區(qū)處于中度以上風(fēng)險(xiǎn)區(qū)。收集得到的58條關(guān)于水稻絕收的記錄,有32條記錄地點(diǎn)位于云南南部的臨滄、普洱和西雙版納地區(qū),20條記錄位于四川東北部的成都、德陽地區(qū)。從歷史災(zāi)情的回代可以發(fā)現(xiàn),80%以上的水稻洪澇災(zāi)害發(fā)生在水稻洪澇災(zāi)害的次高、高風(fēng)險(xiǎn)區(qū)。歷史文字記錄資料雖然一定程度上可以反應(yīng)區(qū)域水稻洪澇發(fā)生的風(fēng)險(xiǎn)狀況,但是由于這些歷史記載對水稻洪澇發(fā)生的范圍、具體時(shí)段和損失等的定量化描述并不全面,因此,在水稻洪澇風(fēng)險(xiǎn)模型的精確、定量驗(yàn)證方面還有一定局限性。

    西南地區(qū)土質(zhì)變異性較大,不同的土壤類型飽水性、滲水性不同,同樣的降水過程,偏黏質(zhì)土壤透水性差,持水性強(qiáng),更易積水;土壤結(jié)構(gòu)松軟的地區(qū)遇暴雨或強(qiáng)降水天氣容易滯水[2]。另外,洪澇發(fā)生與河網(wǎng)密度和距離水體的遠(yuǎn)近有關(guān),河網(wǎng)越密集距離水體越近的地方,發(fā)生洪澇災(zāi)害的危險(xiǎn)程度越高[40]。因此,可以在孕災(zāi)環(huán)境敏感性中量化土壤質(zhì)地、河網(wǎng)密度等因素對洪澇形成的影響。本文暴露性、抗災(zāi)能力評價(jià)模型構(gòu)建的基礎(chǔ)資料是西南地區(qū)縣域尺度水稻農(nóng)業(yè)統(tǒng)計(jì)資料,雖然農(nóng)業(yè)統(tǒng)計(jì)資料最直接、最客觀反映了西南不同地區(qū)水稻生產(chǎn)狀況,但是,縣域內(nèi)部各要素的分布差異還無法辨識??篂?zāi)能力指標(biāo)主要考慮水稻本身的情況,指標(biāo)選擇相對單一,另外,作物產(chǎn)量受氣候適宜性的影響,采用水稻產(chǎn)量作為反應(yīng)區(qū)域抗災(zāi)能力的一個(gè)指標(biāo),無法剔除氣候適宜對水稻生產(chǎn)的正向作用,還需進(jìn)一步將社會經(jīng)濟(jì)因素等考慮進(jìn)去,例如,農(nóng)村勞動力、農(nóng)民人均純收入、農(nóng)業(yè)機(jī)械總動力等。洪澇風(fēng)險(xiǎn)的形成是一個(gè)非常復(fù)雜的過程,每個(gè)影響因子從不同角度、不同程度影響著洪澇風(fēng)險(xiǎn)的形成,以當(dāng)前有限的評估資料和技術(shù)手段全面準(zhǔn)確分析評價(jià)水稻洪澇災(zāi)害風(fēng)險(xiǎn)仍有一定難度。隨著資料的完善,不斷優(yōu)化模型和相應(yīng)因子、指標(biāo),可使評估結(jié)果對水稻洪澇防災(zāi)減災(zāi)更具針對性和指導(dǎo)意義。

    3.2 結(jié)論

    西南地區(qū)水稻移栽分蘗期洪澇危險(xiǎn)性最高,其次為拔節(jié)孕穗期和抽穗成熟期,全生育期洪澇危險(xiǎn)性空間分布則表現(xiàn)為從中間向四周增加。高危險(xiǎn)性區(qū)域主要包括云南南部和東北部、貴州南部,以及四川的成都、眉山和德陽地區(qū)。西南地區(qū)水稻高、次高暴露區(qū)主要集中在四川東北部和重慶地區(qū),該區(qū)域水稻種植密度高于西南地區(qū)平均水平。孕災(zāi)敏感性由四周向中部地區(qū)逐漸增加。云南北部地區(qū)由于地勢低洼且較平坦,是孕災(zāi)環(huán)境高敏感區(qū)域。低抗災(zāi)能力區(qū)主要位于貴州,該區(qū)域水稻生產(chǎn)水平低于西南地區(qū)。

    總體來說,西南地區(qū)水稻洪澇高、次高風(fēng)險(xiǎn)區(qū)主要位于貴州、四川東北部和云南南部地區(qū);中度風(fēng)險(xiǎn)區(qū)主要分布于貴州南部和四川東北部地區(qū);低風(fēng)險(xiǎn)區(qū)主要位于重慶和云南北部地區(qū)。

    參考文獻(xiàn):References

    [1]Zhou F.Trends in total precipitation and frequency of daily precipitation extremes over China[J].Journal of Climate,2005, 18(7):1096-1108.

    [2]鄭大瑋,李茂松,霍治國.農(nóng)業(yè)災(zāi)害與減災(zāi)對策[M].北京:中國農(nóng)業(yè)大學(xué)出版社,2013.

    Zhen D W,Li M S,Huo Z G.Agricultural disasters and reduction countermeasures[M].Beijing:China Agricultural University Press,2013.(in Chinese)

    [3]William J P,Arthur A A.Natural hazard risk assessment and public policy anticipating:unexpected[M].Springer-Verlag New York Inc,1982.

    [4]Jena P P,Chatterjee C,Pradhan G,et al.Are recent frequent high floods in Mahanadi basin in eastern India due to increase in extreme rainfalls[J].Journal of Hydrology,2014,517(1): 847-862.

    [5]Maskrey A.Disaster mitigation:a community based approach[M].Oxford:Oxfam,1989.

    [6]史培軍.三論災(zāi)害研究的理論與實(shí)踐[J].自然災(zāi)害學(xué)報(bào), 2002,11(3):1-9.

    Shi P J.Theory on disaster science and disaster dynamics[J]. Journal of Natural Disasters,2002,11(3):1-9.(in Chinese)

    [7]張繼權(quán),李寧.主要?dú)庀鬄?zāi)害風(fēng)險(xiǎn)評價(jià)與管理的數(shù)量化方法及其應(yīng)用[M].北京:北京師范大學(xué)出版社,2007:27-244.

    Zhang J Q,Li N.Quantitative methods and applications of risk assessment and management on main meteorological disasters[M].Beijing:Beijing Normal University Press, 2007:27-244.(in Chinese)

    [8]Okada N,Tatano H,Hagihara Y,et al.Integrated research on methodological development of urban diagnosis for disaster risk and its applications[J].Disaster Prevention Research Institute Annuals,Kyoto University,2004,47(C):1-8.

    [9]Koks E E,Jongman B,Husby T G,et al.Combining hazard,exposure and social vulnerability to provide lessons for flood risk management[J].Environmental Science & Policy, 2015, 47:42-52.

    [10]Haddad E A,Teixeira E.Economic impacts of natural disasters in megacities:the case of floods in S?o Paulo, Brazil[J].Habitat International,2013,45:106-113.

    [11]Tezuka S,Takiguchi H,Kazama S,et al.Estimation of the effects of climate change on flood-triggered economic losses in Japan[J].International Journal of Disaster Risk Reduction, 2014,9(1):1619-1649.

    [12]吳立,霍治國,張蕾,等.西南地區(qū)農(nóng)業(yè)洪澇等級指標(biāo)構(gòu)建及時(shí)空分布特征分析[J].應(yīng)用生態(tài)學(xué)報(bào),2015, 26(8): 2473- 2481.

    Wu L,Huo Z G,Zhang L,et al.Level indicators construction and temporal-spatial distribution features of agricultural flood in the southwest of China[J].Chinese Journal of Applied Ecology,2015,26(8):2473-2481.(in Chinese)

    [13]盛紹學(xué),霍治國,石磊.江淮地區(qū)小麥澇漬災(zāi)害風(fēng)險(xiǎn)評估與區(qū)劃[J].生態(tài)學(xué)雜志,2010,29(5):985-990.

    Sheng S X,Huo Z G,Shi L.Risk assessment and regionalization of waterlogging disaster for wheat in Jianghuai region[J].Chinese Journal of Ecology, 2010, 29(5):985-990.(in Chinese)

    [14]張蕾,霍治國,黃大鵬,等.10-11月海南省瓜菜苗期濕澇風(fēng)險(xiǎn)評估與區(qū)劃[J].應(yīng)用氣象學(xué)報(bào),2015,26(4):432-441.

    Zhang L,Huo Z G,Huang D P,et al.Assessment and distribution of waterlogging damage risks for melons and vegetables in Hainan province from October to November[J].Journal of Applied Meteorological Science, 2015,26(4):432-441.(in Chinese)

    [15]王越,江志紅,張強(qiáng),等.用Palmer濕潤指數(shù)作西北地區(qū)東部冬小麥旱澇評估[J].應(yīng)用氣象學(xué)報(bào),2008,19(3):342-349.

    Wang Y,Jiang Z H,Zhang Q,et al.Evaluating aridity and wetness of the wheat with Palmer moisture anomaly index in the East of Northwest China[J].Journal of Applied Meteorological Science,2008,19(3):342-349.(in Chinese)

    [16]吳洪顏,高蘋,徐為根,等.江蘇省冬小麥濕漬害的風(fēng)險(xiǎn)區(qū)劃[J].生態(tài)學(xué)報(bào),2012,32(6):1871-1879.

    Wu H Y,Gao P,Xu W G,et al.Risk division on winter wheat suffering from spring wet damages in Jiangsu province[J].Acta Ecologica Sinica,2012,32(6):1871-1879.(in Chinese)

    [17]陸魁東,黃晚華,方麗,等.氣象災(zāi)害指標(biāo)在湖南春玉米種植區(qū)劃中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2007,18(4):548-554.

    Lu K D,Huang W H,Fang L,et al.The climatic zoning of spring maize in Hunan based on meteorological disaster indexs[J].Journal of Applied Meteorological Science,2007, 18(4):548-554.(in Chinese)

    [18]李美榮,李星敏,李艷莉,等.基于連陰雨災(zāi)害指數(shù)的陜西省蘋果生長風(fēng)險(xiǎn)分析[J].干旱氣象,2012,29(1):106-109.

    Li M R,Li X M,Li Y L,et al.Analysis of apple growth risk based on continuous rainy days disaster index in Shanxi province[J].Journal of Arid Meteorology,2012, 29(1):106- 109. (in Chinese)

    [19]吳利紅,婁偉平,柳苗,等.油菜花期降水適宜度變化趨勢及風(fēng)險(xiǎn)評估[J].中國農(nóng)業(yè)科學(xué),2011,44(3):620-626.

    Wu L H, Lou W P,Liu M,et al.Changes of precipitation suitability and assessment of yield loss risk of cole at floweing period[J].Scientia Agricultura Sinica,2011, 44(3): 620-626.(in Chinese)

    [20]Gouranga K,Narayan S,Ashwani K.Deep-water rice production as influenced by time and depth of flooding on the east coast of India[J].Archives of Agronomy and Soil Science,2012,58(6):573-592.

    [21]寧金花,霍治國,陸魁東,等.不同生育期淹澇脅迫對雜交稻形態(tài)特征和產(chǎn)量的影響[J].中國農(nóng)業(yè)氣象,2013,34 (6):678-684.

    Ning J H,Huo Z G,Lu K D,et al.Effects of water logging on morphological characteristics and yield of hybrid rice during growth stages[J].Chinese Journal of Agrometeorology, 2013, 34(6):678-684.(in Chinese)

    [22]吳啟俠,朱建強(qiáng),晏軍,等.澇脅迫對雜交中稻形態(tài)和產(chǎn)量的影響[J].中國農(nóng)業(yè)氣象,2016,37(2):188-198.

    Wu Q X,Zhu J Q,Yan J,et al.Morphology of middle-season hybrid rice in Hubei province and its yield under different waterlogging stresses[J].Chinese Journal of Agrometeor- ology, 2016,37(2):188-198.(in Chinese)

    [23]楊建瑩,霍治國,吳立,等.西南地區(qū)水稻洪澇等級評價(jià)指標(biāo)構(gòu)建及風(fēng)險(xiǎn)分析[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(16):135-144.

    Yang J Y,Huo Z G,Wu L,et al.Evaluation level construction and analysis of the risk on rice flood in Southwest China[J].Transactions of the CSAE,2015,31(16):135-144.(in Chinese)

    [24]王平,史培軍.中國農(nóng)業(yè)自然災(zāi)害綜合區(qū)劃方案[J].自然災(zāi)害學(xué)報(bào),2000,9(4):16-18.

    Wang P,Shi P J.Comprehensive regionalization of agricultural natural disaster in China[J].Journal of Natural Disasters,2000, 9(4):16-18.(in Chinese)

    [25]霍治國,李世奎,王素艷,等.主要農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)評估技術(shù)及其應(yīng)用研究[J].自然資源學(xué)報(bào),2003,6(6):692-703.

    Huo Z G,Li S K,Wang S Y,et al.Study on the risk evaluation technologies of main agrometeorological disasters and their application[J].Journal of Natural Resources,2003, 6(6): 692-703.(in Chinese)

    [26]戴娟,潘益農(nóng),劉青,等.改進(jìn)的AHP在縣域尺度暴雨洪澇風(fēng)險(xiǎn)評價(jià)的應(yīng)用[J].氣象科學(xué),2014,4(4):428-434.

    Dai J,Pan Y N,Liu Q,et al.Application of improved AHP in county-scale rainstorm and flood risk assessment[J].Journal of the Meteorological Sciences,2014,4(4):428-434.(in Chinese)

    [27]王春乙,蔡菁菁,張繼權(quán).基于自然災(zāi)害風(fēng)險(xiǎn)理論的東北地區(qū)玉米干旱、冷害風(fēng)險(xiǎn)評價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2015, 31(6):238-245.

    Wang C Y,Cai J J,Zhang J Q.Risk assessment of drought and chilling injury of maize in Northeast China[J].Transactions of the CSAE,2015,31(6):238-245.(in Chinese)

    [28]Getirana A C V,Bonnet M P,Filho O C R,et al.Improving hydrological information acquisition from DEM processing in floodplains[J].Hydrological Processes,2009,23(3): 502- 514.

    [29]殷潔,裴志遠(yuǎn),陳曦?zé)?等.基于GIS 的武陵山區(qū)洪水災(zāi)害風(fēng)險(xiǎn)評估[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(24):110-117.

    Yin J,Pei Z Y,Chen X W,et al.GIS-based flood disaster risk assessment in Wuling Mountain Region[J].Transactions of the CSAE,2013,29(24):110-117.(in Chinese)

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

    Jia J Y,He N,Han L Y,et al.Analysis on drought risk of maize in Southwest China based on natural disaster risk theory and ArcGIS[J].Transactions of the CSAE,2015,31(4):152-159.(in Chinese)

    [31]溫克剛,詹兆渝.中國氣象災(zāi)害大典(四川卷) [M].北京:氣象出版社,2006.

    Wen K G, Zan Z Y.Editorial board of China meteorological disasters ceremony(Sichuan)[M].Beijing:China Meteorological Press,2006.(in Chinese)

    [32]單琨,劉布春,劉園,等.基于自然災(zāi)害系統(tǒng)理論的遼寧省玉米干旱風(fēng)險(xiǎn)分析[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(8):186-194.

    Shan K,Liu B C,Liu Y,et al.Analysis on drought risk of maize based on natural disaster system theory in Liaoning province[J].Transactions of the CSAE,2012,28(8): 186-194. (in Chinese)

    [33]李慧芳,殷淑燕.海河流域不同等級降水強(qiáng)度和雨日的時(shí)空變化特征[J].中國農(nóng)業(yè)氣象,2014,35(6):603-610.

    Li H F,Yin S Y.Temporal and spatial variation of precipitation intensity and days with different grades in Haihe River Basin[J].Chinese Journal of Agrometeorology, 2014,35(6): 603-610.(in Chinese)

    [34]Qin N,Wang J,Yang G,et al.Spatial and temporal variations of extreme precipitation and temperature events for the Southwest China in 1960-2009[J].Geoenvironmental Disasters,2015,2(1):1-14.

    [35]晏紅明,李清泉,孫丞虎,等.中國西南區(qū)域雨季開始和結(jié)束日期劃分標(biāo)準(zhǔn)的研究[J].大氣科學(xué),2013,37(5):1111-1128.

    Yan H M,Li Q Q,Sun C H,et al.Criterion for determining the onset and end of the rainy season in southwest China[J].Chinese Journal of Atmospheric Sciences,2013, 37(5):1111-1128.(in Chinese)

    [36]劉金華,楊軍昌.影響貴州農(nóng)業(yè)可持續(xù)發(fā)展的因素分析與對策思考[J].貴州農(nóng)業(yè)科學(xué),2005,1(1):101-104.

    Liu J H,Yang J C.Factors analysis and countermeasure on the sustainable development of agriculture in Guizhou[J]. Guizhou Agricultural Sciences,2005,1(1):101-104.(in Chinese)

    [37]溫克剛,劉建華.中國氣象災(zāi)害大典(云南卷)[M].北京:氣象出版社,2006.

    Wen K G,Liu J H.Editorial board of China meteorological disasters ceremony(Yunnan)[M]. Beijing:China Meteorological Press,2006.(in Chinese)

    [38]溫克剛,馬力.中國氣象災(zāi)害大典(重慶卷) [M].北京:氣象出版社,2008.

    Wen K G,Ma L.Editorial board of China meteorological disasters ceremony(Chongqing)[M].Beijing:China Meteorological Press,2008.(in Chinese)

    [39]溫克剛,羅寧.中國氣象災(zāi)害大典(貴州卷) [M].北京:氣象出版社,2006.

    Wen K G,Luo N.Editorial board of China meteorological disasters ceremony(Guizhou)[M].Beijing:China Meteorological Press,2006.(in Chinese)

    [40]莫建飛,陸甲,李艷蘭,等.基于GIS的廣西洪澇災(zāi)害孕災(zāi)環(huán)境敏感性評估[J].災(zāi)害學(xué),2010,25(4):33-37.

    Mo J F,Lu J,Li Y L,et al.GIS-based sensitivity assessment of environment of developing flood hazards in Guangxi province[J].Journal of Catastrophology,2010,25(4):33-37.(in Chinese)

    Flood Risk Assessment and Zoning for Rice in Southwest China

    YANG Jian-ying1, HUO Zhi-guo1,2, WU Li1, ZHANG Gui-xiang1, WANG Tian-ying1

    (1.State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China; 2.Collaborative Innovation Center of Meteorological Disaster Forecast, Early-Warning and Assessment, Nanjing University of Information Science & Technology, Nanjing 210044)

    As the main crop in Southwest China (Sichuan, Chongqing, Guizhou and Yunnan province), rice is widely accepted to be seriously threatened by the flood disaster. Agro-flood risk analysis is helpful for improving the ability of regional disaster management and reducing potential flood risk. In this study, the rice flood risk was analyzed and then zoned using the integrated rice flood risk assessment model including hazard, exposure, sensitivity, and flood prevention and mitigation ability, and based on a data set of the daily rainfall data in 193 meteorological stations from 1961 to 2012, rice production data in 396 counties including plant and yield data from 1981 to 2012, the phenophase data of rice from 17 agrometeorological stations and geographic data in Southwest China. The results indicated that flood hazard risk probabilities variated with phenophase and hazard level. High flood hazard risk happened in transplanting-tillering stage, followed by jointing-booting and tasselling-maturity stages. High-risk areas of flood hazard in whole rice growth stage was detected to be located in south and northeast areas of Yunnan, southern Guizhou, and Chengdu, Meishan and Deyang in Sichuan province. High and sub high exposure zones were mainly located in northeast part of Sichuan province and Chongqing. Subsequently, high and sub high sensitive zones were mainly found in north Yunnan province, south Sichuan province and southeast Guizhou province. Low flood prevention and mitigation zones were located in parts of Guizhou. Finally, the high and sub high zones for integrated rice flood risk were detected in northeast Sichuan province, south Guizhou province and south part of Yunnan province, while south Chongqing and north Yunnan were recognized with low rice flood risk.

    Flood disaster; Agrometeorology; Risk analysis; Zoning; Single cropping rice

    10.3969/j.issn.1000-6362.2016.05.009

    2016-02-04**

    。E-mail:huozhigg@camscma.cn

    基本科研業(yè)務(wù)費(fèi)專項(xiàng)基金(2015Y003);“十二五”國家科技支撐計(jì)劃課題(2012BAD20B02)

    楊建瑩(1985-),女,助理研究員,主要從事農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)預(yù)測與評估研究。E-mail:yangjy@camscma.cn

    猜你喜歡
    孕災(zāi)抗災(zāi)洪澇
    典型山區(qū)斜坡地質(zhì)災(zāi)害孕災(zāi)地質(zhì)條件分析
    尋甸縣塊河流域地質(zhì)災(zāi)害孕災(zāi)地質(zhì)條件分區(qū)評價(jià)
    四川省北川縣崩滑災(zāi)害孕災(zāi)地質(zhì)條件分析
    洪澇造成孟加拉損失25.4萬噸大米
    胡春華強(qiáng)調(diào):切實(shí)抓好春耕生產(chǎn)和抗災(zāi)奪豐收工作
    今年水稻抗災(zāi)奪豐收技術(shù)意見
    王文濤主持召開省政府常務(wù)會議 全力減災(zāi)救災(zāi)確保今年抗災(zāi)奪豐收
    衡陽暴雨洪澇特征及對后期農(nóng)業(yè)干旱影響
    基于孕災(zāi)敏感性分析的高植被覆蓋區(qū)滑坡地質(zhì)災(zāi)害遙感解譯
    打好農(nóng)業(yè)防汛抗災(zāi)主動仗
    種子科技(2016年7期)2016-12-29 08:49:28
    国产熟女午夜一区二区三区 | 亚洲激情五月婷婷啪啪| 免费看光身美女| a级片在线免费高清观看视频| 在线观看人妻少妇| 国产无遮挡羞羞视频在线观看| 我要看黄色一级片免费的| 日日爽夜夜爽网站| 国产精品熟女久久久久浪| 伊人久久国产一区二区| 成年人免费黄色播放视频| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级| 2021少妇久久久久久久久久久| 亚洲精品乱码久久久v下载方式| 在线精品无人区一区二区三| 欧美日韩视频精品一区| 久久综合国产亚洲精品| 狂野欧美激情性bbbbbb| 精品一区二区三卡| 精品熟女少妇av免费看| 国产成人91sexporn| 天堂8中文在线网| 久久人人爽人人片av| 日韩精品免费视频一区二区三区 | 中文天堂在线官网| 国产一区亚洲一区在线观看| 五月天丁香电影| 91成人精品电影| 亚洲av二区三区四区| 久久精品国产自在天天线| 观看美女的网站| 久久人人爽人人片av| 亚洲人成网站在线观看播放| 欧美日韩成人在线一区二区| 一级二级三级毛片免费看| 亚洲欧美日韩另类电影网站| 久久亚洲国产成人精品v| 亚洲精品乱码久久久久久按摩| 久久人人爽人人爽人人片va| 在线 av 中文字幕| 久久国产精品男人的天堂亚洲 | 亚州av有码| 日韩免费高清中文字幕av| 亚洲精品一区蜜桃| 国产成人精品在线电影| 精品国产一区二区久久| 国产成人免费无遮挡视频| 久热久热在线精品观看| 不卡视频在线观看欧美| 免费观看a级毛片全部| 美女xxoo啪啪120秒动态图| 一级黄片播放器| 亚洲人成网站在线观看播放| 制服诱惑二区| 毛片一级片免费看久久久久| 国产亚洲av片在线观看秒播厂| 亚洲人成77777在线视频| 人人妻人人澡人人看| 中文字幕av电影在线播放| 国产精品久久久久久久久免| 色94色欧美一区二区| 国产极品粉嫩免费观看在线 | 人妻一区二区av| 亚洲精品,欧美精品| 久久综合国产亚洲精品| 久热久热在线精品观看| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 能在线免费看毛片的网站| 高清在线视频一区二区三区| 高清黄色对白视频在线免费看| 五月玫瑰六月丁香| 国产精品一区www在线观看| 久久久精品免费免费高清| 母亲3免费完整高清在线观看 | 亚洲精品自拍成人| 亚洲精品国产色婷婷电影| 99热6这里只有精品| 久久久欧美国产精品| 午夜日本视频在线| 在线观看国产h片| 国产精品一二三区在线看| 国产精品不卡视频一区二区| 日韩中字成人| 免费观看的影片在线观看| 国产精品麻豆人妻色哟哟久久| 国产av精品麻豆| 七月丁香在线播放| 99久久精品国产国产毛片| 亚洲av不卡在线观看| 热99久久久久精品小说推荐| 97超碰精品成人国产| 免费看不卡的av| 岛国毛片在线播放| 成人免费观看视频高清| 少妇高潮的动态图| 嘟嘟电影网在线观看| 一本久久精品| 男男h啪啪无遮挡| 一本大道久久a久久精品| 国内精品宾馆在线| 免费黄色在线免费观看| 超碰97精品在线观看| 久久久久久久国产电影| 一区二区日韩欧美中文字幕 | 国产精品99久久久久久久久| 久久久久国产网址| 99视频精品全部免费 在线| 欧美三级亚洲精品| 男女免费视频国产| 国产精品 国内视频| 日本色播在线视频| 精品久久国产蜜桃| 欧美性感艳星| 精品少妇内射三级| 老司机影院成人| 最近的中文字幕免费完整| 99热这里只有是精品在线观看| 亚洲国产成人一精品久久久| 亚洲国产毛片av蜜桃av| 80岁老熟妇乱子伦牲交| 十八禁网站网址无遮挡| 免费av不卡在线播放| 人成视频在线观看免费观看| 哪个播放器可以免费观看大片| 日韩三级伦理在线观看| 美女福利国产在线| 亚洲高清免费不卡视频| 狂野欧美激情性bbbbbb| 天堂俺去俺来也www色官网| 在线观看www视频免费| 国产成人aa在线观看| 内地一区二区视频在线| 久久狼人影院| 国产不卡av网站在线观看| 性高湖久久久久久久久免费观看| 成年美女黄网站色视频大全免费 | 99九九在线精品视频| 亚洲四区av| 国产深夜福利视频在线观看| 在线天堂最新版资源| 久久久国产一区二区| 一边摸一边做爽爽视频免费| 精品国产一区二区三区久久久樱花| 久久人人爽人人爽人人片va| 中文字幕亚洲精品专区| 肉色欧美久久久久久久蜜桃| 久久影院123| 久久久久久久国产电影| 欧美xxⅹ黑人| 国产高清国产精品国产三级| 如何舔出高潮| 大话2 男鬼变身卡| 黑人欧美特级aaaaaa片| 欧美成人午夜免费资源| 91久久精品电影网| 天堂8中文在线网| 乱人伦中国视频| 国产免费一级a男人的天堂| 日韩强制内射视频| 纯流量卡能插随身wifi吗| 全区人妻精品视频| 欧美少妇被猛烈插入视频| 婷婷色av中文字幕| 精品国产露脸久久av麻豆| 九色成人免费人妻av| 国产乱来视频区| 人人妻人人澡人人爽人人夜夜| 亚洲av福利一区| 麻豆乱淫一区二区| 亚洲成人av在线免费| 夫妻午夜视频| 国产在线一区二区三区精| 成人毛片60女人毛片免费| 在线观看国产h片| 九九爱精品视频在线观看| 成人毛片60女人毛片免费| 综合色丁香网| 在线观看免费视频网站a站| 成人手机av| 18禁在线无遮挡免费观看视频| 日本欧美视频一区| 久久99精品国语久久久| 嫩草影院入口| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产毛片在线视频| 一区二区av电影网| 最近的中文字幕免费完整| 女人久久www免费人成看片| 久久久国产一区二区| 欧美亚洲日本最大视频资源| 久久国产精品大桥未久av| 国产精品免费大片| 晚上一个人看的免费电影| 国产精品久久久久久久电影| 成人漫画全彩无遮挡| 亚洲国产精品国产精品| 视频在线观看一区二区三区| 久久国内精品自在自线图片| 欧美人与善性xxx| 国产成人免费观看mmmm| 天堂俺去俺来也www色官网| 最近中文字幕2019免费版| 91久久精品国产一区二区三区| 自线自在国产av| 男的添女的下面高潮视频| 日日爽夜夜爽网站| 精品人妻在线不人妻| 欧美日韩精品成人综合77777| 免费观看的影片在线观看| 久久av网站| 一本色道久久久久久精品综合| a级毛色黄片| 国产av码专区亚洲av| 亚洲国产精品一区三区| 国产高清三级在线| 制服诱惑二区| 成年人午夜在线观看视频| 亚洲国产日韩一区二区| 亚洲综合色网址| 中文字幕亚洲精品专区| 久久免费观看电影| 大片免费播放器 马上看| 热99久久久久精品小说推荐| 欧美精品国产亚洲| 在线亚洲精品国产二区图片欧美 | 中文字幕亚洲精品专区| 国产精品 国内视频| 在线播放无遮挡| 极品人妻少妇av视频| 一边摸一边做爽爽视频免费| 纵有疾风起免费观看全集完整版| 18禁在线无遮挡免费观看视频| 国产精品国产av在线观看| 亚洲国产精品专区欧美| 在线观看www视频免费| 精品熟女少妇av免费看| 高清不卡的av网站| 日韩一区二区视频免费看| 99热国产这里只有精品6| 亚洲人与动物交配视频| 国产av码专区亚洲av| 综合色丁香网| 亚洲精品久久成人aⅴ小说 | 狠狠精品人妻久久久久久综合| 母亲3免费完整高清在线观看 | 最近最新中文字幕免费大全7| 天堂俺去俺来也www色官网| 永久网站在线| 国模一区二区三区四区视频| 女人久久www免费人成看片| 国产视频首页在线观看| 久久久a久久爽久久v久久| 亚洲天堂av无毛| 中文乱码字字幕精品一区二区三区| 伦理电影大哥的女人| 精品一区在线观看国产| tube8黄色片| 一区在线观看完整版| 精品一区二区三区视频在线| 免费观看在线日韩| 中文字幕av电影在线播放| 国产精品国产av在线观看| 最近最新中文字幕免费大全7| 婷婷成人精品国产| 91午夜精品亚洲一区二区三区| 一区二区三区免费毛片| 久久久久久久久久成人| √禁漫天堂资源中文www| 极品少妇高潮喷水抽搐| 亚洲色图 男人天堂 中文字幕 | 中文字幕免费在线视频6| 最新的欧美精品一区二区| 免费观看的影片在线观看| 国产成人精品福利久久| 在线观看三级黄色| 久久精品国产亚洲av天美| 黄片无遮挡物在线观看| 欧美日韩av久久| 美女脱内裤让男人舔精品视频| 久久精品久久精品一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品熟女少妇av免费看| 亚洲成人手机| 国产白丝娇喘喷水9色精品| 亚洲情色 制服丝袜| 精品一品国产午夜福利视频| av又黄又爽大尺度在线免费看| 亚洲精品乱码久久久v下载方式| 黄片无遮挡物在线观看| 在线免费观看不下载黄p国产| 成年美女黄网站色视频大全免费 | 天天躁夜夜躁狠狠久久av| 婷婷色综合www| 欧美日本中文国产一区发布| 美女xxoo啪啪120秒动态图| 秋霞伦理黄片| 有码 亚洲区| 少妇精品久久久久久久| 成人国产麻豆网| 婷婷成人精品国产| 久久99精品国语久久久| 亚洲人与动物交配视频| 亚洲精品aⅴ在线观看| videossex国产| 黄色一级大片看看| 亚洲精品久久午夜乱码| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线不卡| 两个人免费观看高清视频| 成人国语在线视频| 十八禁高潮呻吟视频| 精品久久久久久久久亚洲| 国产精品一二三区在线看| 久久免费观看电影| 久久久久精品性色| 日本黄大片高清| 日韩欧美精品免费久久| 亚洲经典国产精华液单| 日日爽夜夜爽网站| 伊人久久精品亚洲午夜| 国产成人av激情在线播放 | 国产成人免费观看mmmm| 久久精品国产a三级三级三级| 美女主播在线视频| 免费看av在线观看网站| 久久99精品国语久久久| 麻豆乱淫一区二区| 久久久午夜欧美精品| 精品人妻熟女av久视频| 亚洲精品aⅴ在线观看| 大香蕉97超碰在线| 国产av码专区亚洲av| 精品卡一卡二卡四卡免费| 伦理电影大哥的女人| 人人妻人人添人人爽欧美一区卜| 色5月婷婷丁香| 色哟哟·www| 美女中出高潮动态图| 一本色道久久久久久精品综合| 亚洲国产精品国产精品| www.色视频.com| 久久久久久久国产电影| 看非洲黑人一级黄片| 丝袜在线中文字幕| 国产精品一国产av| 在线观看www视频免费| 国产极品天堂在线| 韩国av在线不卡| 亚洲精品亚洲一区二区| 最近2019中文字幕mv第一页| 久久免费观看电影| 少妇 在线观看| 天天躁夜夜躁狠狠久久av| 亚洲欧美中文字幕日韩二区| 日韩一本色道免费dvd| 精品人妻熟女av久视频| 亚洲av欧美aⅴ国产| 国产老妇伦熟女老妇高清| 欧美日韩一区二区视频在线观看视频在线| 精品久久久精品久久久| 国产有黄有色有爽视频| 国语对白做爰xxxⅹ性视频网站| .国产精品久久| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 亚洲成人手机| 熟女电影av网| 国产精品嫩草影院av在线观看| 欧美xxⅹ黑人| 女性被躁到高潮视频| 肉色欧美久久久久久久蜜桃| 国产老妇伦熟女老妇高清| 大话2 男鬼变身卡| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 建设人人有责人人尽责人人享有的| 91午夜精品亚洲一区二区三区| 最近的中文字幕免费完整| 亚洲av中文av极速乱| 99热全是精品| 欧美日韩在线观看h| 欧美3d第一页| 国产极品粉嫩免费观看在线 | 国产午夜精品久久久久久一区二区三区| 久久精品国产亚洲网站| 91国产中文字幕| 国产精品久久久久久精品古装| 成人黄色视频免费在线看| 如何舔出高潮| 国产一区二区在线观看av| 久久精品国产鲁丝片午夜精品| 日本vs欧美在线观看视频| 久久久国产欧美日韩av| 国产亚洲欧美精品永久| 久久毛片免费看一区二区三区| av天堂久久9| 97超碰精品成人国产| 黑人欧美特级aaaaaa片| 免费看光身美女| 黄色毛片三级朝国网站| 欧美日韩综合久久久久久| 最近手机中文字幕大全| 边亲边吃奶的免费视频| 久久精品国产亚洲网站| 一区在线观看完整版| 欧美 日韩 精品 国产| 3wmmmm亚洲av在线观看| 乱人伦中国视频| 色吧在线观看| av专区在线播放| 成人影院久久| 亚洲久久久国产精品| 亚洲精品久久午夜乱码| 国产精品一区二区在线观看99| 插阴视频在线观看视频| 人成视频在线观看免费观看| 国产视频内射| 国产国语露脸激情在线看| 高清在线视频一区二区三区| 亚洲综合色网址| 成人午夜精彩视频在线观看| 成人漫画全彩无遮挡| 日本黄色日本黄色录像| 午夜91福利影院| 午夜福利在线观看免费完整高清在| 久久久久精品性色| 久久99蜜桃精品久久| 久热久热在线精品观看| 亚洲国产精品专区欧美| 免费人成在线观看视频色| 欧美另类一区| 国产一区二区三区综合在线观看 | 十八禁高潮呻吟视频| 各种免费的搞黄视频| 久久婷婷青草| 18禁在线播放成人免费| 久久久久久久久久久免费av| 高清欧美精品videossex| 男的添女的下面高潮视频| 日韩一区二区三区影片| 免费播放大片免费观看视频在线观看| 久久热精品热| 亚洲欧美日韩另类电影网站| 日韩中字成人| 国产一区有黄有色的免费视频| 国产亚洲午夜精品一区二区久久| 在线精品无人区一区二区三| videossex国产| 满18在线观看网站| 国产精品国产三级国产专区5o| 日本午夜av视频| 欧美成人午夜免费资源| 精品久久蜜臀av无| 久久亚洲国产成人精品v| 日日爽夜夜爽网站| 永久网站在线| 国产亚洲一区二区精品| 九九久久精品国产亚洲av麻豆| 免费少妇av软件| 日韩强制内射视频| 国产精品.久久久| 春色校园在线视频观看| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 国产精品国产三级国产av玫瑰| 一区二区日韩欧美中文字幕 | 国产在视频线精品| 亚洲激情五月婷婷啪啪| 狠狠精品人妻久久久久久综合| 最近最新中文字幕免费大全7| 久久这里有精品视频免费| 成人漫画全彩无遮挡| 久久久久久久久久久免费av| 精品一区二区三区视频在线| 如日韩欧美国产精品一区二区三区 | 亚洲av在线观看美女高潮| 777米奇影视久久| a级片在线免费高清观看视频| 久久 成人 亚洲| 永久网站在线| 欧美丝袜亚洲另类| 我的女老师完整版在线观看| 精品人妻一区二区三区麻豆| 亚洲av成人精品一二三区| 超碰97精品在线观看| 亚洲精品日韩av片在线观看| 自线自在国产av| 国产成人精品婷婷| 99国产精品免费福利视频| 美女福利国产在线| 久久久久久久久久久久大奶| 晚上一个人看的免费电影| 一本色道久久久久久精品综合| 久久毛片免费看一区二区三区| 高清不卡的av网站| 午夜福利网站1000一区二区三区| 91精品国产九色| 爱豆传媒免费全集在线观看| 午夜福利,免费看| 成人无遮挡网站| 亚洲精品乱码久久久久久按摩| 亚洲国产成人一精品久久久| 久久人人爽人人爽人人片va| 国产在线视频一区二区| 亚洲精品国产av成人精品| 91午夜精品亚洲一区二区三区| 久久久久久久大尺度免费视频| 日韩中文字幕视频在线看片| 尾随美女入室| 最近手机中文字幕大全| 亚洲精品,欧美精品| 高清在线视频一区二区三区| 日韩强制内射视频| 亚洲成人av在线免费| 久久精品夜色国产| 视频区图区小说| 一级毛片aaaaaa免费看小| 日本爱情动作片www.在线观看| 久久人妻熟女aⅴ| 妹子高潮喷水视频| 黑人高潮一二区| 伦精品一区二区三区| av电影中文网址| 蜜臀久久99精品久久宅男| 中文字幕制服av| 国产成人精品婷婷| 九草在线视频观看| 午夜激情av网站| 久久久久精品久久久久真实原创| 久久这里有精品视频免费| 久久 成人 亚洲| av网站免费在线观看视频| 18禁动态无遮挡网站| 大片电影免费在线观看免费| 九草在线视频观看| 国产精品99久久99久久久不卡 | 男人添女人高潮全过程视频| 国精品久久久久久国模美| 亚洲精品aⅴ在线观看| 亚洲国产成人一精品久久久| 色哟哟·www| 国产白丝娇喘喷水9色精品| 国产精品99久久99久久久不卡 | 秋霞伦理黄片| 精品少妇内射三级| 2022亚洲国产成人精品| videos熟女内射| 久久午夜综合久久蜜桃| 黑丝袜美女国产一区| 99热全是精品| 国产精品国产三级国产av玫瑰| 国产成人一区二区在线| 在现免费观看毛片| 精品国产乱码久久久久久小说| 秋霞伦理黄片| 久久午夜福利片| 免费高清在线观看视频在线观看| 丝袜喷水一区| 久久久久久久精品精品| 国产免费一级a男人的天堂| 午夜精品国产一区二区电影| 一级毛片电影观看| 99视频精品全部免费 在线| 热re99久久精品国产66热6| 亚洲国产欧美在线一区| 黄色一级大片看看| 狂野欧美激情性bbbbbb| 99热网站在线观看| 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 成人黄色视频免费在线看| 日韩电影二区| 18禁在线无遮挡免费观看视频| 国产精品国产三级专区第一集| 又黄又爽又刺激的免费视频.| xxxhd国产人妻xxx| 91久久精品电影网| 精品亚洲乱码少妇综合久久| 免费大片18禁| 天天操日日干夜夜撸| 欧美亚洲日本最大视频资源| 国产精品免费大片| 成人亚洲欧美一区二区av| 欧美3d第一页| 日本猛色少妇xxxxx猛交久久| 伦理电影大哥的女人| 韩国av在线不卡| 国产精品人妻久久久影院| 少妇高潮的动态图| 国产乱来视频区| 亚洲av福利一区| 人人妻人人澡人人爽人人夜夜| 建设人人有责人人尽责人人享有的| 亚洲国产欧美在线一区| 久久久久久久久久久免费av| 久久99热6这里只有精品| 爱豆传媒免费全集在线观看| 国产在线免费精品| 少妇熟女欧美另类| 精品久久蜜臀av无| 久久久国产精品麻豆| 国产亚洲一区二区精品| 国产国拍精品亚洲av在线观看| 夫妻性生交免费视频一级片| 热re99久久精品国产66热6| 国产精品免费大片| 九九在线视频观看精品| 晚上一个人看的免费电影| 国产视频首页在线观看| 国产免费福利视频在线观看| 久久精品国产亚洲av天美| 天堂8中文在线网| 天天躁夜夜躁狠狠久久av|