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

    北方蘋果干旱觸發(fā)判識(shí)方法

    2020-12-03 02:24:48楊建瑩霍治國徐建文王培娟鄔定榮毛紅丹
    關(guān)鍵詞:萌芽期發(fā)育階段盛花期

    楊建瑩,霍治國,2,徐建文,王培娟,鄔定榮,毛紅丹,孔 瑞

    北方蘋果干旱觸發(fā)判識(shí)方法

    楊建瑩1,霍治國1,2※,徐建文3,王培娟1,鄔定榮1,毛紅丹1,孔 瑞1

    (1. 中國氣象科學(xué)研究院生態(tài)與農(nóng)業(yè)氣象研究所,北京 100081;2. 南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,南京 210044;3. 大連市氣象服務(wù)中心,大連 116001)

    基于小樣本歷史災(zāi)害數(shù)據(jù)和長(zhǎng)序列氣象、林果生長(zhǎng)數(shù)據(jù)的林果災(zāi)害判識(shí),對(duì)目前歷史災(zāi)害數(shù)據(jù)匱乏的林果等經(jīng)濟(jì)作物氣象災(zāi)害研究具有重要意義。該研究以中國陜西省富士系蘋果干旱災(zāi)害為例,利用氣象資料、蘋果干旱災(zāi)情史料和富士系蘋果發(fā)育期資料,充分考慮蘋果不同發(fā)育階段的水分需求和降水供給情況,以及前期水分盈虧狀況對(duì)當(dāng)前發(fā)育階段蘋果生長(zhǎng)的影響,在水分盈虧指數(shù)計(jì)算的基礎(chǔ)上,構(gòu)建蘋果干旱指數(shù)。通過概率分析、K-Means聚類、歐式距離等方法,厘定陜西省富士系蘋果的干旱觸發(fā)閾值。采用致災(zāi)因子序列對(duì)比分析、預(yù)留樣本驗(yàn)證相結(jié)合的方法,驗(yàn)證蘋果干旱觸發(fā)閾值有效性。結(jié)果表明:1)蘋果干旱觸發(fā)閾值分別為:蘋果果樹萌動(dòng)-萌芽期0.87,萌芽-盛花期0.84,盛花-成熟期0.73;2)基于閾值提取的蘋果干旱年份的干旱指數(shù)序列與歷史災(zāi)害樣本干旱指數(shù)序列具有同一性;預(yù)留獨(dú)立樣本指標(biāo)判識(shí)準(zhǔn)確率為85.58%;典型站點(diǎn)長(zhǎng)時(shí)間序列檢驗(yàn)判識(shí)結(jié)果準(zhǔn)確率為80.95%。研究結(jié)果可為林果災(zāi)害指標(biāo)研究提供技術(shù)支撐。

    干旱;氣象;蘋果;閾值厘定;有效性檢驗(yàn)

    0 引 言

    干旱是中國最主要的農(nóng)業(yè)氣象災(zāi)害,發(fā)生頻率高、致災(zāi)范圍廣、持續(xù)時(shí)間長(zhǎng)、損失影響大。中國常年干旱受災(zāi)面積約占?xì)庀鬄?zāi)害總受災(zāi)面積的60%以上,平均每年農(nóng)業(yè)干旱面積在2 000~3 000萬hm2[1]。林果干旱是農(nóng)業(yè)干旱的一種,是指在林果生長(zhǎng)發(fā)育過程中,因降水不足并得不到適時(shí)適量的灌溉致使供水不能滿足林果的正常需水而造成損傷[2]。近50年來中國干旱面積迅速擴(kuò)大,極端干旱事件頻發(fā)[3-5],加劇了農(nóng)業(yè)干旱的形成,給林果生產(chǎn)帶來了嚴(yán)重影響。

    農(nóng)業(yè)干旱的發(fā)生時(shí)間、發(fā)展過程和影響范圍,可通過構(gòu)建干旱指標(biāo)來綜合判識(shí)和表達(dá)。從氣象、土壤條件、作物特性出發(fā),以降水[6]、蒸發(fā)[7]、濕度和土壤含水率[8]等指標(biāo)作為主要數(shù)據(jù)源的區(qū)域農(nóng)業(yè)干旱研究相對(duì)豐富,干旱判識(shí)指標(biāo)和方法在農(nóng)田作物干旱識(shí)別[9-10]、時(shí)空發(fā)展趨勢(shì)[11-12]、災(zāi)損量化中應(yīng)用較為廣泛。常用的農(nóng)業(yè)干旱判識(shí)指數(shù)包括降水量距平百分率(Percentage of precipitation anomalies,Pa)[13]、土壤相對(duì)濕度(Relative soil moisture,Rsm)[14-15]、連續(xù)無有效降水(降雪、積雪)日數(shù)(continuous days without available precipitation,Dnp)[16]、作物水分虧缺指數(shù)(Crop Water Deficit Index,CWDI)[17-18]等。針對(duì)不同作物,采用的干旱指標(biāo)也各不相同。2015年發(fā)布的國家標(biāo)準(zhǔn)《農(nóng)業(yè)干旱等級(jí)》(GB/T 32136—2015)[19]中,采用作物水分虧缺距平指數(shù)、土壤相對(duì)濕潤度指數(shù)、農(nóng)田與作物干旱形態(tài)指標(biāo)來進(jìn)行農(nóng)業(yè)干旱的界定。其中,作物水分虧缺距平指數(shù)>35時(shí)認(rèn)為發(fā)生農(nóng)業(yè)干旱,涉及的主要農(nóng)作物包括小麥、玉米、水稻、棉花、大豆、馬鈴薯等。

    蘋果是中國北方地區(qū)農(nóng)業(yè)產(chǎn)業(yè)中經(jīng)濟(jì)效益較高的林果。然而,受大陸性季風(fēng)氣候影響,降水時(shí)空分布不均,且80%的果園無灌溉條件,蘋果生產(chǎn)過程水分供需矛盾突出[20],制約光熱和土肥資源潛力的有效發(fā)揮[21],對(duì)果樹光合產(chǎn)物形成和積累產(chǎn)生顯著影響[22-23]。與農(nóng)田作物農(nóng)業(yè)干旱識(shí)別指標(biāo)相比較,蘋果干旱識(shí)別研究具有特殊性。蘋果生長(zhǎng)周期內(nèi)耗水較多,可達(dá)1 000 mm以上,盛花-成熟期階段是果樹關(guān)鍵需水階段,需水量一般500 mm以上[20];根系發(fā)達(dá),對(duì)地表水分?jǐn)r截、保持和利用作用顯著,對(duì)深層土壤水分吸收能力強(qiáng)。與農(nóng)田作物相比,蘋果對(duì)水分和養(yǎng)分的利用效率,以及對(duì)干旱天氣的抵御能力均高于農(nóng)田作物。目前,國內(nèi)外對(duì)蘋果干旱的研究,主要是基于果園控制試驗(yàn),研究干旱條件對(duì)蘋果生理生化[21,24]、產(chǎn)量[25]、品質(zhì)[26]的影響。蘋果干旱判識(shí)指標(biāo)研究還相對(duì)較少。馬延慶等[27]基于干燥度、土壤含水率、果園蒸散量等農(nóng)業(yè)干旱指標(biāo)綜合分析了蘋果干旱發(fā)生的主要時(shí)段和蘋果的需水耗水規(guī)律;王景紅等[28]結(jié)合蘋果不同發(fā)育階段的水分供需特征和氣象干旱強(qiáng)度,構(gòu)建了黃土高原區(qū)富士系蘋果干旱指數(shù),評(píng)價(jià)了陜西省蘋果干旱風(fēng)險(xiǎn);程雪等[20]利用該指標(biāo)研究了中國北方地區(qū)蘋果干旱時(shí)空分布特征。上述指標(biāo)能較好地表征蘋果不同發(fā)育階段的水分供需特征及區(qū)域尺度蘋果干旱風(fēng)險(xiǎn)的差異,但蘋果干旱發(fā)生的臨界氣象條件及觸發(fā)閾值尚未明確,蘋果干旱判識(shí)針對(duì)性有所不足。相關(guān)研究表明,適度干旱脅迫會(huì)促進(jìn)同化物從營養(yǎng)器官向生殖器官轉(zhuǎn)移,產(chǎn)生增產(chǎn)節(jié)水效應(yīng)的同時(shí),增強(qiáng)蘋果品質(zhì)[29];但若干旱天氣持續(xù)加強(qiáng)且無補(bǔ)充灌溉,則發(fā)生蘋果干旱,影響果實(shí)產(chǎn)量和品質(zhì)。明確蘋果干旱觸發(fā)閾值、準(zhǔn)確判識(shí)蘋果干旱的發(fā)生發(fā)展,對(duì)針對(duì)性的開展蘋果防災(zāi)減災(zāi)意義重大。

    歷史災(zāi)害數(shù)據(jù)對(duì)于理解區(qū)域?yàn)?zāi)害系統(tǒng)的相互作用機(jī)制有著重要的作用?;跉v史災(zāi)害樣本重建與再分析的農(nóng)業(yè)氣象指標(biāo)構(gòu)建方法,在農(nóng)業(yè)洪澇[30]、澇漬[31-32]、高溫?zé)岷33]、冷害[34]等災(zāi)害觸發(fā)閾值研究中得到應(yīng)用和證實(shí)。然而,與農(nóng)田作物相比,林果等經(jīng)濟(jì)作物的歷史災(zāi)害記錄相對(duì)匱乏,單純采用災(zāi)情反演與再分析的災(zāi)害研究方法存在樣本不足、精度不高等問題。因此,如何綜合利用小樣本歷史災(zāi)害數(shù)據(jù)、長(zhǎng)序列氣象和林果生長(zhǎng)數(shù)據(jù),實(shí)現(xiàn)區(qū)域林果氣象災(zāi)害的閾值厘定和判識(shí),對(duì)林果等經(jīng)濟(jì)作物氣象災(zāi)害的監(jiān)測(cè)預(yù)警與防控有重要意義。

    鑒于此,本文以陜西富士系蘋果干旱災(zāi)害為例,參考農(nóng)田作物農(nóng)業(yè)干旱指數(shù)構(gòu)建方法,在綜合考慮蘋果不同發(fā)育階段水分供需的基礎(chǔ)上,構(gòu)建能表征蘋果階段水分狀況的蘋果干旱指數(shù);通過災(zāi)情反演、指數(shù)提取、聚類分析和樣本距離等方法,研發(fā)基于小樣本災(zāi)害數(shù)據(jù)、長(zhǎng)序列歷史氣象數(shù)據(jù)和林果生長(zhǎng)數(shù)據(jù)的災(zāi)害判識(shí)方法,厘定陜西富士系蘋果的干旱觸發(fā)閾值,并采用致災(zāi)因子序列對(duì)比分析、預(yù)留樣本驗(yàn)證相結(jié)合的方法,驗(yàn)證蘋果干旱觸發(fā)指標(biāo)有效性。研究結(jié)果可為林果等經(jīng)濟(jì)作物災(zāi)害判識(shí)指標(biāo)構(gòu)建提供范式,為針對(duì)性地開展林果干旱防災(zāi)減災(zāi)提供技術(shù)支撐。

    1 材料與方法

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

    陜西省內(nèi)35個(gè)國家級(jí)基本氣象站1981—2018年逐日氣象資料來源于國家氣象信息中心,包括逐日的降水量、平均風(fēng)速、平均溫度、日最高氣溫和日最低氣溫等,對(duì)其中個(gè)別缺測(cè)值采用該日相鄰兩日該要素的平均值替代[35]。蘋果發(fā)育期包括果樹萌動(dòng)期、萌芽期、盛花期和成熟期,其中蘋果果樹萌動(dòng)期以平均溫度穩(wěn)定通過3 ℃為依據(jù)確定[36]。蘋果萌芽期、盛花期和成熟期數(shù)據(jù)來源于《中國農(nóng)業(yè)氣象資源圖集》[37],通過紙質(zhì)資料掃描、ArcGIS地理校正等方法,對(duì)發(fā)育期資料進(jìn)行數(shù)字化。提取等值線附近氣象站點(diǎn)發(fā)育期日期,利用Kriging插值結(jié)合物候定律法[38]將蘋果發(fā)育期數(shù)據(jù)插值到研究區(qū)域內(nèi)的35個(gè)氣象站點(diǎn)。蘋果干旱災(zāi)情資料來源于《中國氣象災(zāi)害大典》(陜西)[39]、《中國氣象災(zāi)害年鑒》[40]、蘋果災(zāi)情調(diào)查、果樹災(zāi)害專題調(diào)查,以及果業(yè)基地縣縣志、果業(yè)基地縣民政部門的災(zāi)情記錄以及媒體報(bào)道。

    1.2 蘋果干旱指數(shù)構(gòu)建

    參考農(nóng)田作物(如小麥、玉米等)干旱指標(biāo)構(gòu)建方法,充分考慮蘋果不同發(fā)育階段的水分需求和降水供給情況,以及前期水分盈虧狀況對(duì)當(dāng)前發(fā)育階段蘋果生長(zhǎng)的影響,在蘋果不同發(fā)育階段水分盈虧指數(shù)計(jì)算的基礎(chǔ)上,構(gòu)建蘋果干旱指數(shù)。

    1.2.1 水分盈虧指數(shù)及其計(jì)算方法

    借鑒農(nóng)作物干旱識(shí)別指標(biāo),水分盈虧指數(shù)計(jì)算公式如下所示。

    式中,WD為某年某站點(diǎn)蘋果第個(gè)發(fā)育階段內(nèi)的水分盈虧指數(shù);ETP分別為某年某站點(diǎn)蘋果第個(gè)發(fā)育階段需水量和降水量,mm。作物需水量(ET)采用FAO推薦的“參考作物蒸散量乘以作物系數(shù)法”計(jì)算蘋果各發(fā)育階段需水量,具體計(jì)算公式如下

    式中ET為蘋果某一發(fā)育階段的需水量,mm;ET0為對(duì)應(yīng)時(shí)段的參考作物蒸散量,mm,計(jì)算方法采用FAO推薦的Penman-Monteith公式[41];K為蘋果該發(fā)育階段作物系數(shù),參考1998年FAO-56推薦的分段單值平均作物系數(shù)法[41]和相關(guān)文獻(xiàn)[42-43],蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期3個(gè)發(fā)育階的K分別為0.65、0.8和0.95。

    1.2.2 蘋果干旱指數(shù)及其計(jì)算方法

    果樹萌動(dòng)-萌芽期和萌芽-盛花期一般歷時(shí)約30~40 d,當(dāng)前發(fā)育階段之前的降水條件分別對(duì)果樹萌動(dòng)-萌芽期和萌芽-盛花期的干旱形成影響較大;盛花-成熟期歷時(shí)60~100 d,該階段雖然處于一年中的水分充盈期,但如遇干旱年,果樹盛花期前的水分條件依然影響當(dāng)前階段的水分供給。參考國家標(biāo)準(zhǔn)《農(nóng)業(yè)干旱等級(jí)》(GB/T 32136—2015)中作物水分虧缺指數(shù) CWDI(Crop Water Deficit Index)的構(gòu)建方法,以及林果作物與農(nóng)田作物在生理需水、根系水分吸收及土壤深層水分利用上的差異性,在構(gòu)建蘋果干旱指數(shù)時(shí),考慮當(dāng)前階段及前60 d水分盈虧情況。采用權(quán)重遞減的思路[44],即假定當(dāng)前的水分盈虧指數(shù)對(duì)該階段的干旱指數(shù)貢獻(xiàn)最大,隨著時(shí)間的前移,過去60 d的水分盈虧指數(shù)對(duì)當(dāng)前階段的干旱指數(shù)隨時(shí)間距離遞減??紤]到水分盈虧指數(shù)適用于旬以上尺度的干旱監(jiān)測(cè)和評(píng)估,故以10 d為步長(zhǎng),累計(jì)統(tǒng)計(jì)過去60 d(即6個(gè)步長(zhǎng))的水分盈虧指數(shù),并設(shè)置所有因子總權(quán)重為1。WD的權(quán)重系數(shù)為7/28;WD-1即前1~10 d的水分盈虧指數(shù),它的權(quán)重系數(shù)為6/28;WD-2即前11~20 d的水分盈虧指數(shù),它的權(quán)重系數(shù)為5/28;以此類推。

    具體公式如下

    式中DI為第階段的蘋果干旱指數(shù),最大值為1;為計(jì)算水分盈虧指數(shù)的步長(zhǎng)個(gè)數(shù),為1~7。

    1.3 歷史災(zāi)害樣本反演

    根據(jù)災(zāi)情記錄中的蘋果干旱發(fā)生時(shí)間、地點(diǎn),結(jié)合蘋果發(fā)育期資料,分別反演對(duì)應(yīng)年份、對(duì)應(yīng)氣象站點(diǎn)、對(duì)應(yīng)發(fā)育階段(蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期)的蘋果干旱指數(shù),構(gòu)建蘋果3個(gè)發(fā)育階段的災(zāi)害樣本干旱指數(shù)集合。例如災(zāi)情記錄為“1994年6—9月,受干旱影響,陜西禮泉縣全縣果園受災(zāi),果樹整體發(fā)育生長(zhǎng)緩慢,生理性落葉、落果”。根據(jù)蘋果發(fā)育期資料,6—9月蘋果處于盛花-成熟期間,計(jì)算1994年禮泉縣盛花-成熟期蘋果干旱指數(shù)為0.79,因此反演得到該樣本為:1994,禮泉,盛花-成熟期,0.79。同理,反演得到歷史蘋果干旱實(shí)際受災(zāi)樣本序列3組共130個(gè),其中蘋果果樹萌動(dòng)-萌芽期干旱樣本19個(gè),萌芽-盛花期干旱樣本41個(gè)、盛花-成熟期干旱樣本70個(gè)。

    在以往基于歷史災(zāi)害樣本重建與再分析的農(nóng)業(yè)氣象災(zāi)害研究中,通常預(yù)留10%~20%的獨(dú)立樣本進(jìn)行閾值驗(yàn)證,閾值厘定樣本數(shù)量通常為10~50個(gè)[30-34]。為研發(fā)基于小樣本災(zāi)害數(shù)據(jù)的林果干旱閾值判識(shí)方法,本文隨機(jī)選取蘋果各發(fā)育階段災(zāi)害總樣本量的20%用于蘋果干旱觸發(fā)閾值厘定,預(yù)留80%災(zāi)害樣本用于閾值的驗(yàn)證。蘋果干旱樣本數(shù)量詳細(xì)信息如表1所示。

    表1 蘋果干旱樣本數(shù)量詳細(xì)信息

    1.4 蘋果干旱觸發(fā)閾值厘定

    1.4.1 蘋果干旱初始判識(shí)指標(biāo)

    概率密度或累積概率函數(shù)可以表征總體樣本特征,天氣學(xué)中,常被應(yīng)用于定義特定地區(qū)和時(shí)間的極端天氣事件[45]。本文采用累積概率密度表征1981—2018年蘋果各發(fā)育期干旱指數(shù)的總體分布。計(jì)算得到1981—2018年35個(gè)氣象站點(diǎn)蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期干旱指數(shù)樣本各1 330個(gè)。假設(shè)蘋果各發(fā)育階段的1 330個(gè)干旱指數(shù)為連續(xù)型隨機(jī)變量,分別對(duì)干旱指數(shù)總體樣本進(jìn)行累計(jì)概率擬合。以5%為步長(zhǎng),計(jì)算總體蘋果干旱指數(shù)累積概率下的歷史災(zāi)害樣本覆蓋率()。以>0的起始值所對(duì)應(yīng)的總體樣本累積概率反函數(shù)值,作為該發(fā)育階段蘋果干旱的初始判識(shí)指標(biāo),認(rèn)為只有干旱指數(shù)值達(dá)到初始判識(shí)指標(biāo),才有可能發(fā)生蘋果干旱事件。災(zāi)害樣本覆蓋率(i)的計(jì)算方法如下式所示

    式中R為蘋果第個(gè)發(fā)育階段的災(zāi)害樣本覆蓋率;n為該發(fā)育階段達(dá)到1981—2018年總體樣本某一累積概率的歷史災(zāi)害樣本數(shù);N為蘋果第個(gè)發(fā)育階段歷史災(zāi)害樣本數(shù)總數(shù)。

    1.4.2 觸發(fā)閾值厘定方法

    以初始判識(shí)指標(biāo)為基礎(chǔ),將1981—2018年干旱指數(shù)大于初始判識(shí)指標(biāo)的干旱指數(shù)值定義為聚類樣本。分別對(duì)蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期聚類樣本進(jìn)行分類,判識(shí)聚類樣本中心點(diǎn);計(jì)算歷史災(zāi)害樣本干旱指數(shù)與聚類中心點(diǎn)距離,以最小距離對(duì)應(yīng)的聚類中心點(diǎn)作為該階段蘋果干旱發(fā)生的觸發(fā)閾值。

    1)基于K-Means算法的聚類中心點(diǎn)識(shí)別

    聚類分析是依據(jù)研究對(duì)象的個(gè)體特征建立分類的多元統(tǒng)計(jì)分析方法。它能夠?qū)⒁慌鷺颖緮?shù)據(jù)按照性質(zhì)上的親疏程度對(duì)其進(jìn)行分類,在經(jīng)濟(jì)、管理、社會(huì)學(xué)等領(lǐng)域都有廣泛應(yīng)用[46]。K-Means算法是常用的聚類分析方法之一[47],其基本思想是,對(duì)于給定的樣本集,按照樣本之間的距離大小,將樣本集劃分為個(gè)簇。讓簇內(nèi)的點(diǎn)盡量緊密的連在一起,而讓簇間的距離盡量的大。假設(shè)簇劃分為(1,2,...C),則目標(biāo)是最小化平方誤差

    式中為聚類數(shù);為簇C的樣本點(diǎn);u是簇C的均值向量,也稱為聚類中心點(diǎn),表達(dá)式為

    本文采用K-Means算法,對(duì)1981—2018年蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期聚類樣本進(jìn)行分類,判識(shí)聚類樣本中心點(diǎn)。

    2)基于最小距離的干旱觸發(fā)閾值厘定

    歐式距離一定程度上能表征樣本間的相似程度,是分類研究中的常用計(jì)算方法。本文采用歐式距離公式計(jì)算蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期歷史災(zāi)害樣本與聚類樣本中心點(diǎn)距離,判識(shí)歷史災(zāi)害樣本所屬聚類類別,以最小歐式距離對(duì)應(yīng)的聚類中心點(diǎn),作為該階段蘋果干旱發(fā)生的觸發(fā)閾值。

    歐式距離公式[48]計(jì)算方法為

    式中d是第個(gè)災(zāi)害樣本和聚類中心點(diǎn)的歐式距離;是災(zāi)害樣本數(shù)量;Z是第個(gè)樣本的干旱指數(shù)值;Z是聚類中心點(diǎn)對(duì)應(yīng)的干旱指數(shù)值。

    1.5 蘋果干旱閾值合理性檢驗(yàn)

    1.5.1 致災(zāi)因子(干旱指數(shù))序列的統(tǒng)計(jì)學(xué)檢驗(yàn)

    計(jì)算1981—2018年干旱指數(shù),基于蘋果干旱觸發(fā)閾值,將提取得到的干旱年份干旱指數(shù)序列與歷史記錄反演得到的災(zāi)害樣本干旱指數(shù)序列作對(duì)比,計(jì)算干旱指數(shù)分布頻率、平均值和標(biāo)準(zhǔn)差,并對(duì)3個(gè)發(fā)育階段的共3組樣本集合做檢驗(yàn),驗(yàn)證是否來著同一整體,評(píng)估基于閾值提取的蘋果干旱年份的干旱指數(shù)序列與歷史災(zāi)害樣本干旱指數(shù)序列是否具有同一性。

    1.5.2 預(yù)留獨(dú)立樣本檢驗(yàn)

    1)歷史災(zāi)害樣本判識(shí)準(zhǔn)確率檢驗(yàn)

    參考水稻洪澇、高溫?zé)岷Φ燃?jí)閾值檢驗(yàn)方法[30,33],采用預(yù)留驗(yàn)證樣本,依據(jù)蘋果干旱觸發(fā)閾值,計(jì)算蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期干旱判識(shí)準(zhǔn)確率,進(jìn)行不同發(fā)育階段干旱判別結(jié)果檢驗(yàn)。

    2)典型站干旱判識(shí)檢驗(yàn)

    采用單一站點(diǎn)1981—2018年干旱指數(shù)值為典型站災(zāi)害判識(shí)序列,根據(jù)驗(yàn)證樣本中記載的蘋果干旱發(fā)生年份和發(fā)育階段,對(duì)照檢驗(yàn)基于干旱指標(biāo)計(jì)算的干旱判識(shí)結(jié)果,驗(yàn)證其與歷史記錄是否一致。站點(diǎn)選擇依據(jù)以下原則:①該站常年有蘋果種植,是蘋果生產(chǎn)的典型區(qū);②該站歷史蘋果干旱災(zāi)情記錄相對(duì)詳實(shí)。

    2 結(jié)果與分析

    2.1 蘋果干旱觸發(fā)閾值厘定

    2.1.1 初始判識(shí)指標(biāo)的確定

    1981—2018年干旱指數(shù)累積概率及歷史災(zāi)害樣本干旱指數(shù)分布如圖1所示。以1981—2018年干旱指數(shù)序列的5%累積概率為步長(zhǎng),逐步計(jì)算不同累積概率下的災(zāi)害樣本覆蓋率(),其結(jié)果如圖2所示。隨著干旱指數(shù)累積概率的增加,干旱指數(shù)(Drought Index, DI)逐漸增大,災(zāi)害樣本覆蓋率()逐漸增大。以蘋果果樹萌動(dòng)-萌芽期初始判識(shí)指標(biāo)確定過程為例,初始判識(shí)指標(biāo)的確定蘋果果樹萌動(dòng)-萌芽期的4個(gè)歷史災(zāi)害樣本干旱指數(shù)均分布于1981—2018年干旱指數(shù)序列的75%累積概率以上,80%、85%和90%累積概率下的值分別為25%、25%和100%,即當(dāng)蘋果干旱指數(shù)達(dá)到1981—2018年干旱指數(shù)序列的80%和85%反函數(shù)值時(shí),25%的歷史災(zāi)害樣本干旱指數(shù)值小于該指標(biāo),不能被判識(shí);當(dāng)達(dá)到90%反函數(shù)值時(shí),全部歷史災(zāi)害樣本的干旱指數(shù)值在該判識(shí)指標(biāo)以下,均無法通過該指標(biāo)得到判識(shí)。以5%累積概率為步長(zhǎng),判斷>0的起始值為1981—2018年干旱指數(shù)的75%累積概率,全部歷史災(zāi)害樣本干旱指數(shù)值均大于該指標(biāo),均可得到判識(shí)。以75%累積概率反函數(shù)值作為該發(fā)育階段蘋果干旱的初始判識(shí)指標(biāo),判定1981—2018年干旱指數(shù)序列的75%累積概率反函數(shù)值(0.74)作為該發(fā)育階段蘋果干旱的初始判識(shí)指標(biāo)。同理,判定蘋果萌芽-盛花期、盛花-成熟期的干旱初步判識(shí)指標(biāo)為1981—2018年干旱指數(shù)80%、70%累積概率對(duì)應(yīng)干旱指數(shù)值,即0.75、0.50。蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期干旱的初始判識(shí)指標(biāo)分別為0.74、0.75、0.50,只有蘋果干旱指數(shù)達(dá)到0.74、0.75和0.50時(shí),蘋果干旱才可能發(fā)生。

    圖1 不同發(fā)育期1981—2018年干旱指數(shù)(DI)累積概率及歷史災(zāi)害樣本干旱指數(shù)分布

    圖2 不同發(fā)育期1981—2018年干旱指數(shù)累積概率下歷史災(zāi)害樣本覆蓋率(R)變化

    2.1.2 閾值厘定

    以0.74、0.75、0.50為初始判識(shí)指標(biāo),將1981—2018年干旱指數(shù)大于初始判識(shí)指標(biāo)的干旱指數(shù)值定義為聚類樣本,采用K-Mean聚類將蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期聚類樣本進(jìn)行分類。參考K-Mean聚類算法的使用和適用原則,將蘋果不同發(fā)育階段聚類樣本分為5個(gè)聚類。聚類中心點(diǎn)及方差分析如表2所示。對(duì)1981—2018年蘋果果樹萌動(dòng)-萌芽期干旱指數(shù)>0.74的干旱指數(shù)樣本進(jìn)行聚類,得到的聚類中心點(diǎn)分別為:0.81、0.87、0.92、0.96、0.99;分別對(duì)1981—2018年萌芽-盛花期干旱指數(shù)>0.75和盛花-成熟期干旱指數(shù)>0.50的干旱指數(shù)進(jìn)行聚類,得到聚類中心點(diǎn)分別為0.78、0.84、0.90、0.95、0.99和0.65、0.73、0.80、0.86、0.91。組間方差分別為0.198、0.228、0.378,大于組內(nèi)方差(<0.01)(表2)。

    表2 聚類中心點(diǎn)及方差分析

    采用歐式距離公式計(jì)算得到蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期災(zāi)害樣本與聚類樣本中心點(diǎn)距離。采用離差標(biāo)準(zhǔn)化方法,對(duì)不同發(fā)育階段的歐氏距離進(jìn)行線性變換,使歐氏距離值映射到[0,1] 區(qū)間,結(jié)果如圖3所示。蘋果果樹萌動(dòng)-萌芽期災(zāi)害樣本對(duì)應(yīng)聚類中心點(diǎn)0.81、0.87、0.92、0.96和0.99的歐式距離分別為0.38、0.10、0.46、0.77和1.00。以最小歐式距離對(duì)應(yīng)的聚類中心點(diǎn)0.87,作為該階段蘋果干旱的觸發(fā)閾值,即蘋果果樹萌動(dòng)-萌芽期干旱發(fā)生的觸發(fā)閾值為0.87,認(rèn)為該階段干旱指數(shù)達(dá)到0.87時(shí),發(fā)生蘋果干旱災(zāi)害。同理,確定蘋果萌芽-盛花期和盛花-成熟期的干旱觸發(fā)閾值分別為0.84和0.73(圖3)??傮w來說,蘋果干旱觸發(fā)閾值由大到小順序?yàn)樘O果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期。蘋果果樹萌動(dòng)-萌芽期與萌芽-盛花期干旱標(biāo)閾值較為接近,均>0.8。盛花-成熟期時(shí)間跨度大,覆蓋整個(gè)夏季,降雨充沛,相較于果樹萌動(dòng)-萌芽期和萌芽-盛花期,氣象干旱條件略有緩解,干旱觸發(fā)閾值較蘋果果樹萌動(dòng)-萌芽期和萌芽-盛花期偏低。

    圖3 不同發(fā)育期歷史災(zāi)害樣本干旱指數(shù)與各聚類樣本中心點(diǎn)的歐式距離

    2.2 蘋果干旱觸發(fā)閾值檢驗(yàn)

    2.2.1 干旱指數(shù)序列的同一性檢驗(yàn)

    基于蘋果干旱觸發(fā)閾值判識(shí)1981—2018年蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期和盛花-成熟期的干旱年份,提取干旱年份的干旱指數(shù)序列。提取得到的1981—2018年蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期和盛花-成熟期干旱樣本數(shù)分別149、142和110個(gè);預(yù)留的歷史災(zāi)害樣本數(shù)分別為15、33和56個(gè)。對(duì)比蘋果不同發(fā)育階段干旱指數(shù)序列,結(jié)果如圖4所示。基于閾值提取得到的干旱年份干旱指數(shù)序列與歷史災(zāi)害干旱指數(shù)序列的頻率分布規(guī)律基本一致。蘋果果樹萌動(dòng)-萌芽期干旱指數(shù)主要分布于0.87~0.92,提取樣本與災(zāi)害樣本占比分別為72%和80%。萌芽-盛花期干旱指數(shù)分布頻率最高的數(shù)值區(qū)域?yàn)閇0.84 0.86),提取樣本與災(zāi)害樣本干旱指數(shù)在該數(shù)值區(qū)域占比分別為33.10%和33.33%;其次為[0.86 0.88),提取樣本與災(zāi)害樣本干旱指數(shù)占比分別為20.42%和20.22%。盛花-成熟期干旱指數(shù)隨著指數(shù)增加分布頻率逐漸遞減,[0.73 0.76)數(shù)值區(qū)域的干旱指數(shù)分布頻率最高,提取樣本與災(zāi)害樣本干旱指數(shù)在該數(shù)值區(qū)域占比分別為31.00%和32.42%。3個(gè)發(fā)育階段干旱指數(shù)平均值和標(biāo)準(zhǔn)差如表3所示。蘋果果樹萌動(dòng)-萌芽期提取樣本與災(zāi)害樣本干旱指數(shù)平均值分別為0.901和0.899;萌芽-盛花期提取樣本與災(zāi)害樣本干旱指數(shù)平均值分別為0.881和0.877;盛花-成熟期取樣本與災(zāi)害樣本干旱指數(shù)平均值分別為0.794和0.792。3個(gè)發(fā)育階段的3對(duì)樣本做獨(dú)立樣本檢驗(yàn),檢驗(yàn)顯著性分別為0.356、0.264和0.871,均大于0.05,結(jié)果表明,3個(gè)發(fā)育階段的樣本分別來自同一整體,表明基于閾值提取的蘋果干旱年份的干旱指數(shù)序列與歷史災(zāi)害樣本干旱指數(shù)序列具有同一性。

    2.2.2 預(yù)留獨(dú)立樣本檢驗(yàn)

    采用預(yù)留的104個(gè)蘋果歷史干旱樣本,根據(jù)記錄中的受災(zāi)時(shí)間和地點(diǎn),提取蘋果干旱指數(shù),根據(jù)構(gòu)建的蘋果觸發(fā)指標(biāo)閾值,判斷蘋果是否受災(zāi),進(jìn)行預(yù)留樣本的干旱判別準(zhǔn)確率檢驗(yàn)。由圖5可以看出,果樹萌動(dòng)-萌芽期的14個(gè)災(zāi)害樣本中,11個(gè)樣本可通過干旱指數(shù)計(jì)算判識(shí)得到,3個(gè)樣本干旱指數(shù)值置于閾值線以下,蘋果干旱的指標(biāo)判識(shí)結(jié)果與實(shí)際發(fā)生的吻合率為78.57%;萌芽-盛花期的33個(gè)樣本中,27個(gè)樣本干旱指數(shù)顯示蘋果干旱,6個(gè)樣本未能通過指標(biāo)判識(shí)得到,該階段蘋果干旱的指標(biāo)判識(shí)結(jié)果與實(shí)際發(fā)生的吻合率為81.82%;盛花-成熟期的56個(gè)災(zāi)害樣本中,5個(gè)樣本未通過指標(biāo)判識(shí)得到,51個(gè)樣本干旱指數(shù)計(jì)算結(jié)果顯示該階段蘋果發(fā)生干旱,吻合率為91.07%??傮w來說,104個(gè)獨(dú)立樣本中,89個(gè)判識(shí)結(jié)果與實(shí)際一致,指標(biāo)判識(shí)驗(yàn)證的準(zhǔn)確率為85.58%。

    圖4 基于觸發(fā)閾值的干旱指數(shù)與歷史災(zāi)害樣本干旱指數(shù)頻率特征

    表3 基于觸發(fā)閾值的干旱指數(shù)與歷史災(zāi)害樣本干旱指數(shù)序列對(duì)比分析及T檢驗(yàn)

    圖5 預(yù)留災(zāi)害樣本反演與判識(shí)準(zhǔn)確性檢驗(yàn)

    2.2.3 典型站干旱特征及閾值檢驗(yàn)

    對(duì)比各地區(qū)歷史災(zāi)情記錄及預(yù)留災(zāi)害樣本情況,選擇禮泉作為干旱判識(shí)指標(biāo)有效性檢驗(yàn)典型站。計(jì)算禮泉站1981—2018年干旱指數(shù)值,得到禮泉蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期干旱指數(shù)時(shí)序變化特征,如圖6所示。根據(jù)災(zāi)情樣本中記載的禮泉地區(qū)蘋果干旱發(fā)生時(shí)間,提取該站發(fā)生蘋果干旱的年份和發(fā)育階段,對(duì)照檢驗(yàn)干旱觸發(fā)閾值的判識(shí)結(jié)果。歷史災(zāi)情反演得到禮泉地區(qū)蘋果果樹萌動(dòng)-萌芽期干旱年份,分別是1985、1986、1995、1999和2001年。依據(jù)1981—2018年干旱指數(shù)和觸發(fā)閾值,4個(gè)年份禮泉蘋果果樹萌動(dòng)-萌芽期判識(shí)為干旱,分別為1985、1986、1995和2001年,1999年判識(shí)結(jié)果為無旱,與實(shí)際結(jié)果有差異。果樹萌動(dòng)-萌芽期的指標(biāo)判識(shí)吻合率為80%。歷史災(zāi)情反演得到禮泉地區(qū)蘋果萌芽-盛花期干旱樣本5個(gè),分別是1986、1995、1997、2008和2012年。通過干旱指數(shù)計(jì)算,1986、1995、2008和2012年判識(shí)結(jié)果為干旱,1997年判識(shí)結(jié)果為無旱。萌芽-盛花期的判識(shí)吻合率為80%。歷史災(zāi)情記錄中關(guān)于禮泉蘋果盛花-成熟期的干旱記載較完善,災(zāi)情反演得到禮泉地區(qū)蘋果盛花-成熟期干旱樣本11個(gè),分別是1988—1991、1994—1996、2002、2005—2006和2016年。通過干旱指數(shù)計(jì)算,11個(gè)年中的9個(gè)年份判識(shí)結(jié)果為干旱,2個(gè)年份判識(shí)結(jié)果為無旱,判識(shí)吻合率為81.82%??傮w來說,典型站點(diǎn)長(zhǎng)時(shí)間序列檢驗(yàn)樣本共21個(gè),判識(shí)結(jié)果與實(shí)際干旱發(fā)生一致的樣本17個(gè),判識(shí)吻合率為80.95%。結(jié)果表明,構(gòu)建的陜西蘋果干旱判識(shí)指標(biāo),可以用于氣候變化背景下長(zhǎng)時(shí)間序列典型站的蘋果干旱判識(shí)。

    圖6 禮泉站蘋果不同發(fā)育階段1981—2018干旱指數(shù)時(shí)序變化特征

    3 討 論

    本文參考農(nóng)田作物干旱指標(biāo)構(gòu)建方法,考慮蘋果水分供需特征的同時(shí),進(jìn)一步考慮前期水分狀況對(duì)當(dāng)前階段蘋果干旱形成的影響。農(nóng)田作物小麥、玉米等干旱研究中,通??紤]前50 d的水分供需情況來評(píng)價(jià)當(dāng)前階段作物干旱程度[17,19]??紤]到林果作物在土壤底墑利用、抗干旱能力等方面均優(yōu)于農(nóng)田作物,采用當(dāng)前階段和前60 d水分供需狀況構(gòu)建蘋果干旱指數(shù)。與以往單純考慮當(dāng)前階段水分供需的蘋果干旱研究相比,本文構(gòu)建的蘋果干旱指數(shù)更具有理論可行性。

    以往的災(zāi)害觸發(fā)判識(shí)研究中,閾值厘定往往采用人工或半人工方法,主觀因素干擾較大。受歷史災(zāi)情記錄數(shù)量和質(zhì)量的限制,災(zāi)害樣本往往單純用于指標(biāo)的驗(yàn)證[30,49]。本文深入挖掘1981—2018年蘋果3個(gè)發(fā)育階段干旱指數(shù)和歷史災(zāi)害樣本干旱指數(shù)的特點(diǎn)和聯(lián)系,將有限歷史災(zāi)害樣本信息,投射到1981—2018年長(zhǎng)序列樣本的擬合和分類規(guī)律中,研發(fā)適用于林果小樣本災(zāi)害資料的閾值厘定方法。該研究方法保證了小樣本數(shù)據(jù)(歷史災(zāi)害樣本干旱指數(shù))與長(zhǎng)序列數(shù)據(jù)(1981—2018年干旱指數(shù)序列)規(guī)律表達(dá)上的全面型、規(guī)范性和一致性,指標(biāo)閾值厘定方法客觀、科學(xué)。蘋果干旱觸發(fā)閾值果樹萌動(dòng)-萌芽期>萌芽-盛花期>盛花-成熟期,與程雪等[20]的研究結(jié)論基本一致。蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期正處春季少雨時(shí)期,干旱指數(shù)普遍偏高;盛花-成熟期蘋果生長(zhǎng)干旱環(huán)境略有緩解,蘋果干旱觸發(fā)閾值也略有下降。結(jié)合相關(guān)文獻(xiàn)[50]及災(zāi)情記錄,蘋果果樹萌動(dòng)-萌芽期若發(fā)生干旱(干旱指數(shù)大于0.87),則影響蘋果樹體及花芽萌發(fā),導(dǎo)致萌芽延遲或萌芽不整齊,影響新梢生長(zhǎng);萌芽-盛花期干旱如遇干旱天氣(干旱指數(shù)大于0.84),則出現(xiàn)嚴(yán)重落花;盛花-成熟期階段是蘋果關(guān)鍵需水臨界期,花后40 d內(nèi)如果發(fā)生干旱(干旱指數(shù)大于0.73),則影響果肉細(xì)胞分裂。果實(shí)膨大期干旱條件直接影響蘋果產(chǎn)量和品質(zhì)。

    值得注意的是,本文隨機(jī)選取20%的災(zāi)害樣本用于指標(biāo)閾值厘定,盡管樣本資料可以較為客觀地表征歷史災(zāi)害發(fā)生特點(diǎn),但難以避免某些樣本的代表性較差,抽樣比例差異亦使閾值厘定結(jié)果存在偏差。通常情況下,樣本量越大,抽樣誤差就越小,對(duì)規(guī)律的模擬越準(zhǔn)確。但是,考慮到本文的研究目的,少量樣本被用于指標(biāo)閾值的厘定,預(yù)留了大量歷史災(zāi)害樣本用于閾值的檢驗(yàn)與驗(yàn)證。盡管預(yù)留了80%的災(zāi)害樣本用于指標(biāo)閾值的驗(yàn)證,但驗(yàn)證樣本依然無法實(shí)現(xiàn)時(shí)序和空間尺度全覆蓋,尤其是蘋果開花前災(zāi)情數(shù)據(jù)稀缺,難以在空間尺度對(duì)典型年份蘋果干旱分布進(jìn)行區(qū)域性驗(yàn)證。驗(yàn)證樣本的區(qū)域分布特征、樣本年代際分布的豐富程度、典型站點(diǎn)數(shù)量及地理特征等因素一定程度上影響蘋果不同發(fā)育階段干旱判識(shí)的準(zhǔn)確率。隨著災(zāi)情資料的不斷完善,還需逐漸增加指標(biāo)構(gòu)建及檢驗(yàn)樣本,綜合氣象、水文、地理、果園管理等多要素,不斷完善和優(yōu)化干旱閾值和判識(shí)方法,提高蘋果干旱判識(shí)準(zhǔn)確率。

    4 結(jié) 論

    本文充分考慮蘋果不同發(fā)育階段的水分需求和降水供給情況,以及前期水分盈虧狀況對(duì)當(dāng)前發(fā)育階段蘋果生長(zhǎng)的影響,在水分盈虧指數(shù)計(jì)算的基礎(chǔ)上,構(gòu)建蘋果干旱指數(shù)。通過概率分析、K-Means聚類、歐式距離等方法,綜合分析1981—2018年蘋果干旱指數(shù)與歷史災(zāi)害樣本干旱指數(shù)分布與分類特征,研發(fā)基于小樣本歷史災(zāi)害數(shù)據(jù)和長(zhǎng)序列氣象、林果生長(zhǎng)數(shù)據(jù)的林果災(zāi)害判識(shí)方法,厘定適用于陜西省富士系蘋果的干旱觸發(fā)判識(shí)閾值,主要結(jié)論如下:

    1)蘋果干旱觸發(fā)閾值由大到小順序?yàn)樘O果果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期。蘋果果樹萌動(dòng)-萌芽期、萌芽-盛花期正處春季少雨時(shí)期,干旱指數(shù)偏高,干旱觸發(fā)閾值分別為0.87和0.84;盛花-成熟期時(shí)間跨度大,覆蓋整個(gè)夏季,降雨充沛,相較于果樹萌動(dòng)-萌芽期和萌芽-盛花期干旱略有緩解,但該時(shí)期為蘋果需水關(guān)鍵期,干旱觸發(fā)閾值為0.73。

    2)基于閾值提取的蘋果干旱年份的干旱指數(shù)序列與歷史災(zāi)害樣本干旱指數(shù)序列具有同一性,果樹萌動(dòng)-萌芽期、萌芽-盛花期、盛花-成熟期干旱指數(shù)分布頻率最高的數(shù)值區(qū)域?yàn)閇0.87 0.92)、[0.84 0.86)和[0.73 0.76);預(yù)留獨(dú)立樣本指標(biāo)判識(shí)驗(yàn)證的準(zhǔn)確率為85.58%;典型站點(diǎn)(禮泉站)長(zhǎng)時(shí)間序列檢驗(yàn)判識(shí)吻合率為80.95%。

    [1]王春乙. 中國重大農(nóng)業(yè)氣象災(zāi)害研究[M]. 北京:氣象出版社. 2010.

    [2]柏秦鳳,霍治國,王景紅,等. 中國主要果樹氣象災(zāi)害指標(biāo)研究進(jìn)展[J]. 果樹學(xué)報(bào),2019,36(9):1229-1243. Bai Qinfeng, Huo Zhiguo, Wang Jinghong, et al. Progress in research on meteorological disaster indicators of major fruit trees in China[J]. Journal of Fruit Science, 2019, 36(9): 1229-1243. (in Chinese with English abstract)

    [3]孫悅,高繼卿,楊曉光. 西北各省季節(jié)降水變化及其貢獻(xiàn)的差異分析[J]. 中國農(nóng)業(yè)氣象,2019,40(8):489-501. Sun Yue, Gao Jiqing, Yang Xiaoguang. Seasonal variation about the occurrence frequency of different types of precipitation and their relative contribution over Northwest China[J]. Chinese Journal of Agrometeorology, 2019, 40(8): 489-501. (in Chinese with English abstract)

    [4]齊月,陳海燕,房世波,等. 1961—2010 年西北地區(qū)極端氣候事件變化特征[J]. 干旱氣象,2015,33(6):963-969. Qi Yue, Chen Haiyan, Fang Shibo, et al. Variation characteristics of extreme climate events in Northwest China during 1961-2010[J]. Journal of Arid Meteorology, 2015, 33(6): 963-969. (in Chinese with English abstract)

    [5]徐建文,居輝,劉勤,等. 黃淮海地區(qū)干旱變化特征及其對(duì)氣候變化的響應(yīng)[J]. 生態(tài)學(xué)報(bào),2014,34(2):460-470. Xu Jianwen, Ju Hui, Liu Qin, et al. Variation of drought and regional response to climate change in Huang-Huai-Hai Plain[J]. Acta Ecologica Sinica, 2014, 34(2): 460-470. (in Chinese with English abstract)

    [6]鞠笑生,鄒旭愷,張強(qiáng). 氣候旱澇指標(biāo)方法及其分析[J]. 自然災(zāi)害學(xué)報(bào),1998,7(3):51-57. Ju Xiaosheng, Zou Xukai, Zhang Qiang. The method of the climatic drought-flood index and its analysis[J]. Journal of Natural Disasters, 1998, 7(3): 51-57. (in Chinese with English abstract)

    [7]Yang Jianying, Liu Qin, Mei Xurong, et al. Spatiotemporal characteristics of reference evapotranspiration and its sensitivity coefficients to climate factors in Huang-Huai-Hai Plain, China[J]. Journal of Integrative Agriculture, 2013, 12(12): 2280-2291

    [8]劉曉靜,馬東來,張繼權(quán),等. 遼西北地區(qū)玉米干旱的識(shí)別方法[J]. 干旱區(qū)研究,2018,35(3):722-727. Liu Xiaojing, Ma Donglai, Zhang Jiquan, et al. Identification method of maize drought in northwest Liaoning province[J]. Arid Zone Research, 2018, 35(3): 722-727. (in Chinese with English abstract)

    [9]Wu Xia, Wang Peijuan, Huo Zhiguo, et al. Crop drought identification index for winter wheat based on evapotranspiration in the Huang-Huai-Hai Plain, China[J]. Agriculture, Ecosystems and Environment, 2018, 263: 18-30.

    [10]侯陳瑤,朱秀芳,肖名忠,等. 基于游程理論和Copula函數(shù)的遼寧省農(nóng)業(yè)氣象干旱特征研究[J]. 災(zāi)害學(xué),2019,34(2):222-227. Hou Chenyao, Zhu Xiufang, Xiao Mingzhong, et al. Research on agricultural meteorological drought in Liaoning province based on run theory and copula function[J]. Journal of Catastrophology, 2019, 34(2): 222-227. (in Chinese with English abstract)

    [11]高超,李學(xué)文,孫艷偉,等. 淮河流域夏玉米生育階段需水量及農(nóng)業(yè)干旱時(shí)空特征[J]. 作物學(xué)報(bào),2019,45(2):297-309. Gao Chao, Li Xuewen, Sun Yanwei, et al. Spatiotemporal characteristics of water requirement and agricultural drought during summer maize season in Huaihe River Basin[J]. Acta Agronomica Sinica, 2019, 45(2): 297-309. (in Chinese with English abstract)

    [12]李新堯,楊聯(lián)安,聶紅梅,等. 基于植被狀態(tài)指數(shù)的陜西省農(nóng)業(yè)干旱時(shí)空動(dòng)態(tài)[J]. 生態(tài)學(xué)雜志,2018,37(4):1172-1180. Li Xinyao, Yang Lian’an, Nie Hongmei. Assessment of temporal and spatial dynamics of agricultural drought in Shaanxi Province based on vegetation condition index[J]. Chinese Journal of Ecology, 2018, 37(4): 1172-1180. (in Chinese with English abstract)

    [13]陳少勇,郭俊瑞,吳超. 基于降水量距平百分率的中國西南和華南地區(qū)的冬旱特征[J]. 自然災(zāi)害學(xué)報(bào),2015,24(1):23-31. Chen Shaoyong, Guo Junrui, Wu Chao. Characteristics of winter drought in southwest-south China based on precipitation anomaly percentage[J]. Journal of Natural Disasters, 2015, 24(1): 23-31. (in Chinese with English abstract)

    [14]宋艷玲,王建林,田靳峰,等. 氣象干旱指數(shù)在東北春玉米干旱監(jiān)測(cè)中的改進(jìn)[J]. 應(yīng)用氣象學(xué)報(bào),2019,30(1):25-34. Song Yanling, Wang Jianlin, Tian Jinfeng, et al. The spring maize drought index in Northeast China based on meteorological drought index[J]. Journal of Applied Meteorological Science, 2019, 30(1): 25-34. (in Chinese with English abstract)

    [15]安雪麗,武建軍,周洪奎,等. 土壤相對(duì)濕度在東北地區(qū)農(nóng)業(yè)干旱監(jiān)測(cè)中的適用性分析[J]. 地理研究,2017,36(5):837-849. An Xueli, Wu Jianjun, Zhou Hongkui, et al. Assessing the relative soil moisture for agricultural drought monitoring in Northeast China[J]. Geographical Research, 2017, 36(5): 837-849. (in Chinese with English abstract)

    [16]黃晚華,隋月,楊曉光,等. 基于連續(xù)無有效降水日數(shù)指標(biāo)的中國南方作物干旱時(shí)空特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2014,30(4):125-135. Huang Wanhua, Sui Yue, Yang Xiaoguang, et al. Spatio-temporal characteristics of crop drought in southern China based on drought index of continuous days without available precipitation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(4): 125-135. (in Chinese with English abstract)

    [17]李崇瑞,游松財(cái),武永峰,等. 改進(jìn)作物水分虧缺指數(shù)用于東北地區(qū)春玉米干旱災(zāi)變監(jiān)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(21):175-185. Li Chongrui, You Songcai, Wu Yongfeng, et al. Improved crop water deficit index for monitoring drought disaster change process of spring maize in the Northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(21): 175-185. (in Chinese with English abstract)

    [18]張力,陳阜,雷永登. 近60年河北省冬小麥干旱風(fēng)險(xiǎn)時(shí)空規(guī)律[J]. 作物學(xué)報(bào),2019,45(9):1407-1415. Zhang Li, Chen Fu, Lei Yongdeng. Spatial and temporal patterns of drought risk for winter wheat grown in Hebei province in past 60 years[J]. Acta Agronomica Sinica, 2019, 45(9): 1407-1415. (in Chinese with English abstract)

    [19]中華人民共和國國家標(biāo)準(zhǔn). 農(nóng)業(yè)干旱等級(jí):GB/T 32136—2015[S]. 北京:中國標(biāo)準(zhǔn)出版社. 2015.

    [20]程雪,孫爽,張方亮,等. 我國北方地區(qū)蘋果干旱時(shí)空分布特征[J]. 應(yīng)用氣象學(xué)報(bào),2020,31(1):63-73. Cheng Xue, Sun Shuang, Zhang Fangliang, et al. Spatial and temporal distributions of apple drought in Northern China[J]. Journal of Applied Meteorological Science, 2020, 31(1): 63-73. (in Chinese with English abstract)

    [21]Yang Qiliang, Zhang Fucang, Li Fusheng. Effect of different drip irrigation methods and fertilization on growth, physiology and water use of young apple tree[J]. Scientia Horticulture, 2011, 129(1): 119-126.

    [22]郭愛霞,石曉昀,王延秀,等. 干旱脅迫對(duì)3種蘋果砧木葉片光合、葉綠體超微結(jié)構(gòu)和抗氧化系統(tǒng)的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2019,37(1):178-186. Guo Aixia, Shi Xiaoyun, Wang Yanxiu, et al. Effect of drought stress on the photosynthesis, chloroplast ultrastructure and antioxidant system in leaves of three apple rootstocks[J]. Agricultural Research in the Arid Areas, 2019, 37(1): 178-186. (in Chinese with English abstract)

    [23]王冰,丁錫強(qiáng),林曲鳳,等. 2016年煙臺(tái)蘋果與氣象條件分析[J]. 山西果樹,2017(4):18-19. Wang Bing, Ding Xiqiang, Lin Qufeng, et al. Analysis of apple and meteorological conditions of Yantai in 2016[J]. Shanxi Fruits, 2017(4): 18-19. (in Chinese with English abstract)

    [24]Yang Q L, Zhang F C, Li F S. Effect of different drip irrigation methods and fertilization on growth, physiology and water use of young apple tree[J]. Scientia Horticulture, 2011, 129(1): 119-126.

    [25]尤彩虹,景娟,馬遠(yuǎn)衛(wèi). 寶塔區(qū)蘋果生長(zhǎng)主要?dú)庀鬄?zāi)害及應(yīng)對(duì)措施[J]. 農(nóng)業(yè)災(zāi)害研究,2019,9(3):84-85. You Caihong, Jing Juan, Ma Yuanwei. Main meteorological disasters and counter measures of apple growing in Baota District[J]. Journal of Agricultural Catastrophe, 2019, 9(3): 84-85. (in Chinese with English abstract)

    [26]王元基. 干旱對(duì)蘋果品質(zhì)的影響及其與糖代謝的關(guān)系[D]. 楊凌:西北農(nóng)林科技大學(xué),2017. Wang Yuanji. Effect of Drought on Apple Quality and Its Relationship with Sugar Metabolism[D]. Yangling: Northwest A&F University, 2017. (in Chinese with English abstract)

    [27]馬延慶,徐志達(dá),劉長(zhǎng)民,等.陜西渭北旱塬蘋果種植區(qū)域農(nóng)田水分特征分析[J]. 干旱地區(qū)農(nóng)業(yè)研究,2009,27(2):54-59. Ma Yanqing, Xu Zhida, Liu Changmin, et al. Analysis of characteristics of farmland moisture content in apple producing area of Weibei dry plateau in Shaanxi[J]. Agricultural Research in the Arid Areas, 2009, 27(2): 54-59.(in Chinese with English abstract)

    [28]王景紅,柏秦鳳,梁軼,等. 陜西蘋果干旱指數(shù)研究及基于縣域單元的蘋果干旱風(fēng)險(xiǎn)分布[J]. 氣象科技,2014,42(3):516-523. Wang Jinghong, Bai Qinfeng, Liang Yi, et al. Apple drought index and apple drought risk distribution in Shaanxi at county level[J]. Meteorological Science and Technology, 2014, 42(3): 516-523. (in Chinese with English abstract)

    [29]Wang Yuanji, Liu Li, Wang Ying, et al. Effects of soil water stress on fruit yield, quality and their relationship with sugar metabolism in 'Gala' apple[J]. Scientia Horticulturae, 2019, 258: 108753.

    [30]楊建瑩,霍治國,吳立,等. 西南地區(qū)水稻洪澇等級(jí)評(píng)價(jià)指標(biāo)構(gòu)建及風(fēng)險(xiǎn)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(16):135-144. Yang Jianying, Huo Zhiguo, Wu Li, et al. Evaluation level construction and analysis of risk on rice flood in Southwest China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(16): 135-144. (in Chinese with English abstract)

    [31]汪天穎,霍治國,楊建瑩,等. 湖南晚稻洪澇過程等級(jí)指標(biāo)構(gòu)建與演變特征[J]. 應(yīng)用氣象學(xué)報(bào),2019,30(1):35-48. Wang Tianying, Huo Zhiguo, Yang Jianying, et al. Process grade indicator construction and evolution characteristics of late rice flood in Hunan[J]. Journal of Applied Meteorological Science, 2019, 30(1): 35-48. (in Chinese with English abstract)

    [32]楊宏毅,霍治國,楊建瑩,等.江漢和江南西部春玉米澇漬指標(biāo)及風(fēng)險(xiǎn)評(píng)估[J]. 應(yīng)用氣象學(xué)報(bào),2017,28(2):237-246. Yang Hongyi, Huo Zhiguo, Yang Jianying, et al. Indicators and risk of spring corn waterlogging disaster in Jianghan and west region of Jiangnan[J]. Journal of Applied Meteorological Science, 2017, 28(2): 237-246. (in Chinese with English abstract)

    [33]楊建瑩,霍治國,王培娟,等. 江西早稻高溫?zé)岷Φ燃?jí)動(dòng)態(tài)判識(shí)及時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào),2020,31(1):199-207. Yang Jianying, Huo Zhiguo, Wang Peijuan, et al. Dynamic identification of double-early rice heat and its spatiotemporal characteristics in Jiangxi Province, China[J]. Chinese Journal of Applied Ecology, 2020, 31(1): 199-207. (in Chinese with English abstract)

    [34]王培娟,霍治國,楊建瑩,等. 基于熱量指數(shù)的東北春玉米冷害指標(biāo)[J]. 應(yīng)用氣象學(xué)報(bào),2019,30(1):13-24. Wang Peijuan, Huo Zhiguo, Yang Jianying, et al. Indicators of chilling damage for spring maize based on heat index in Northeast China[J]. Journal of Applied Meteorological Science, 2019, 30(1): 13-24. (in Chinese with English abstract)

    [35]吳霞,王培娟,霍治國,等. 1961—2015 年中國潛在蒸散時(shí)空變化特征與成因[J]. 資源科學(xué),2017,39(5):964-977. Wu Xia, Wang Peijuan, Huo Zhiguo, et al. Spatio-temporal distribution characteristics of potential evapotranspiration and impact factors in China from 1961 to 2015[J]. Resources Science, 2017, 39(5): 964-977. (in Chinese with English abstract)

    [36]曾驤,劉金銅. 蘋果與氣象[M]. 北京:氣象出版社,1988.

    [37]梅旭榮. 中國農(nóng)業(yè)氣象資源圖集[M]. 杭州:浙江科學(xué)技術(shù)出版社,2015.

    [38]龔高法,簡(jiǎn)慰民. 我國植物物候期的地理分布[J]. 地理學(xué)報(bào),1983,50(1):33-40. Gong Gaofa, Jian Weimin. Geographical distribution of plant phenology in China[J]. Acta Geographica Sinica, 1983,50(1): 33-40. (in Chinese with English abstract)

    [39]溫克剛,翟佑安. 中國氣象災(zāi)害大典(陜西卷)[M]. 北京:氣象出版社,2005.

    [40]中國氣象局.中國氣象災(zāi)害年鑒[M]. 北京:氣象出版社,2006-2015.

    [41]Allen R G, Pereira L S, Rase D, et al. Crop Evaportranspiration Guide Lines for Computing Crop Water Requirements[M]. RomeP: FAO Irrigation and Drainage Paper, 1998.

    [42]龔道枝,康紹忠,張建華,等. 蘋果樹蒸發(fā)蒸騰量的測(cè)定和計(jì)算[J]. 沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào),2004,35(5):429-431. Gong Daozhi, Kang Shaozhong, Zhang Jianhua, et al. Measuring and estimating evapotranspiration of an apple (. ) orchard[J]. Journal of Shenyang Agricultural University, 2004, 35(5): 429-431. (in Chinese with English abstract)

    [43]左余寶,田昌玉,唐繼偉,等. 魯北地區(qū)主要作物不同生育期需水量和作物系數(shù)的試驗(yàn)研究[J]. 中國農(nóng)業(yè)氣象,2009,30(1):70-73. Zuo Yubao, Tian Changyu, Tang Jiwei, et al. Studies on ETc and Kc of main crops in Northern Shandong province[J]. Chinese Journal of Agrometeorology, 2009, 30(1): 70-73. (in Chinese with English abstract)

    [44]李娜,霍治國,錢錦霞,等. 基于改進(jìn)后相對(duì)濕潤度指數(shù)的山西省氣象干旱時(shí)空特征[J]. 生態(tài)學(xué)雜志,2019,38(7):2249-2257. Li Na, Huo Zhiguo, Qian Jinxia, et al. Spatiotemporal distribution of drought in Shanxi Province based on modified relative moisture index[J]. Chinese Journal of Ecology, 2019, 38(7): 2249-2257. (in Chinese with English abstract)

    [45]Yang Jianying, Huo Zhiguo, Li Xiangxiang, et al. Hot weather event-based characteristics of double-early rice heat risk: A study of Jiangxi province, South China[J]. Ecological Indicators, 2020, 113: 106148.

    [46]劉大海,李寧,晁陽. SPSS 15. 0統(tǒng)計(jì)分析從入門到精通[M]. 北京:清華大學(xué)出版社,2008.

    [47]王冰,伯玥,姜如茵,等. 基于聚類分析的膠東半島蘋果種植氣候區(qū)劃[J]. 江蘇農(nóng)業(yè)科學(xué),2019,47(7):134-138. Wang Bing, Bo Yue, Jiang Ruyin, et al. Application of optimized K-means clustering algorithm in China offshore climate regionalization[J]. Jiangsu Agricultural Sciences, 2019, 47(7): 134-138. (in Chinese with English abstract)

    [48]邱美娟,劉布春,劉園,等. 兩種不同產(chǎn)量歷史豐歉氣象影響指數(shù)確定方法在農(nóng)業(yè)氣象產(chǎn)量預(yù)報(bào)中的對(duì)比研究[J]. 氣象與環(huán)境科學(xué),2019,42(1):41-46. Qiu Meijuan, Liu Buchun, Liu Yuan, et al. Comparative study of two different methods for determining meteorological impact index of historical yield in agrometeorological yield prediction[J]. Meteorological and Environmental Sciences, 2019, 42(1): 41-46. (in Chinese with English abstract)

    [49]王春乙,蔡菁菁,張繼權(quán). 基于自然災(zāi)害風(fēng)險(xiǎn)理論的東北地區(qū)玉米干旱、冷害風(fēng)險(xiǎn)評(píng)價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(6):238-245. Wang Chunyi, Cai Jingjing, Zhang Jiquan. Risk assessment of drought and chilling injury of maize in Northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(6): 238-245. (in Chinese with English abstract)

    [50]伏健民,束懷瑞. 春季干旱對(duì)金冠蘋果不同部位葉片衰老和脫落的影響[J]. 果樹學(xué)報(bào),1993,10(2):65-68. Fu Jianmin, Shu Huairui. The effect of spring drought on the senescence and shedding of different parts of the golden crown apple[J]. Journal of Fruit Science, 1993, 10(2): 65-68. (in Chinese with English abstract)

    Method for trigger identification of apple drought in northern China

    Yang Jianying1, Huo Zhiguo1,2※, Xu Jianwen3, Wang Peijuan1, Wu Dingrong1, Mao Hongdan1, Kong Rui1

    (1.,100081,; 2.,,210044,; 3.116001,)

    Historical disaster documents can greatly contribute to machine learning from disaster experiences, particularly in understanding the interaction mechanism of regional disaster systems. Historical, agricultural, and meteorological disaster characteristics can be explored by the representation and reanalysis of disaster processes and events, integrating meteorological data and agricultural information.However, compared with field crops, the historical disaster records of cash crops, such as forests and fruits, are relatively scarce, which can make the disaster representation and reanalysis with low accuracy. Therefore, it is of great merit to develop a method for identifying forest and fruit disasters, based on the limited historical disaster data, as well as long series meteorological and fruit growth data, for the cash crop meteorological disasters researches which are lack of historical disaster records. In this study, taking the drought disaster of Fuji Apple in Shaanxi Province as an example, the meteorological data, historical drought disaster, and phenological data were integrated to identify the trigger threshold of apple drought. According to the phenological data of apple collected in this study, tree growth stages, that is, tree germinating to flower budding, flower budding to full bloom, and full bloom to mature, were analyzed as target growth stages for apple. Referring to the drought index construction method for crops, such as wheat and corn, pre- and current water demand and precipitation supply were fully considered in the construction of the apple Drought Index (DI). The probability analysis, K-means clustering, and Euclidean distance were used to comprehensively analyze the distribution and classification characteristics of DI between 35 stations from 1981 to 2018 and historical disaster samples. According to the Euclidean distance between the DI in historical disaster sample and the center point of the cluster sample, the drought trigger threshold of Fuji Apple in tree germinating to flower budding, flower budding to full bloom and full bloom to mature were identified by the corresponding minimum Euclidean distance. Afterwards, the trigger threshold was verified by comparing the sequence of disaster-causing factors and reserving samples. The results showed that: 1) The trigger thresholds of apple drought in tree germinating to flower budding, flower budding to full bloom and full bloom to mature were 0.87, 0.84, and 0.73, respectively. 2) The DI sequences that extracted based on the threshold value in tree germinating to flower budding, flower budding to full bloom and full bloom to mature stages were detected the same characters with that in historical disaster samples. The apple drought data that identified by the calculated of DI and trigger threshold were generally consistent with that disaster records in historical documents, with 85.58% of trigger threshold-based results completely consistent with historical records. The identification coincidence rate was 80.95% in the long-time series validation for typical sites. Generally, the trigger threshold of apple drought can provide a sound technical support for apple drought monitoring, early warning, and assessment in Northern China. The agrometeorological disaster trigger identification method based on small samples of historical disaster data can also offer a paradigm in the current research on the meteorological disasters of cash crops with insufficient historical disaster data.

    drought; meteorology; apple; threshold determination; validity test

    楊建瑩,霍治國,徐建文,等. 北方蘋果干旱觸發(fā)判識(shí)方法[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(18):99-108.doi:10.11975/j.issn.1002-6819.2020.18.013 http://www.tcsae.org

    Yang Jianying, Huo Zhiguo, Xu Jianwen, et al. Method for trigger identification of apple drought in northern China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(18): 99-108. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.18.013 http://www.tcsae.org

    2020-06-11

    2020-08-12

    國家重點(diǎn)研發(fā)計(jì)劃(2017YFC1502801;2019YFD1002203);中國氣象科學(xué)研究院科技發(fā)展基金(2019KJ008)

    楊建瑩,副研究員,主要從事農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)預(yù)測(cè)與評(píng)估研究。Email:jyyang@cma.gov.cn

    霍治國,研究員,博士生導(dǎo)師,主要從事農(nóng)業(yè)氣象災(zāi)害、病蟲害預(yù)測(cè)與評(píng)估研究。Email:huozg@ cma.gov.cn

    10.11975/j.issn.1002-6819.2020.18.013

    A

    1002-6819(2020)-18-0099-10

    猜你喜歡
    萌芽期發(fā)育階段盛花期
    麗水市桂花盛花期與氣候條件的關(guān)系研究
    花海泛波春意濃
    小麥生殖發(fā)育階段對(duì)低溫的敏感性鑒定
    花海泛波春意濃
    萌芽期|春秋戰(zhàn)國
    對(duì)森工林區(qū)在商品林基地培養(yǎng)速生楊樹探討
    冠縣春季低溫霜凍與早熟大櫻桃凍害氣候分析
    不同玉米品種萌芽期抗旱篩選
    根域加溫對(duì)塑料大棚內(nèi)葡萄芽萌發(fā)及生理特性的影響
    大花黃牡丹葉片發(fā)育過程中氣孔密度和氣孔指數(shù)的動(dòng)態(tài)變化
    高清在线视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 成人漫画全彩无遮挡| 亚洲精品国产av成人精品| 纯流量卡能插随身wifi吗| 欧美老熟妇乱子伦牲交| freevideosex欧美| 丰满乱子伦码专区| 欧美精品亚洲一区二区| av视频免费观看在线观看| 久久人妻熟女aⅴ| www.色视频.com| 久久精品久久精品一区二区三区| 中文字幕亚洲精品专区| 成人免费观看视频高清| 99久久精品热视频| 性色avwww在线观看| 久久久久久人妻| 色5月婷婷丁香| 亚洲精品乱码久久久久久按摩| 中国三级夫妇交换| 日本欧美国产在线视频| 大又大粗又爽又黄少妇毛片口| 国产av码专区亚洲av| 日韩制服骚丝袜av| 91精品国产国语对白视频| 啦啦啦啦在线视频资源| 精品少妇黑人巨大在线播放| 观看av在线不卡| 午夜激情久久久久久久| 日日摸夜夜添夜夜爱| 国产成人aa在线观看| 一级黄片播放器| 日韩精品免费视频一区二区三区 | 国产在视频线精品| 99热这里只有是精品50| 高清黄色对白视频在线免费看 | 99久久精品一区二区三区| 日本-黄色视频高清免费观看| www.色视频.com| 一区二区三区免费毛片| 高清毛片免费看| 亚洲中文av在线| 最黄视频免费看| 嫩草影院新地址| 免费看av在线观看网站| 国产成人精品福利久久| 99热国产这里只有精品6| 免费观看性生交大片5| 久久精品久久久久久久性| 国产精品不卡视频一区二区| 国产午夜精品久久久久久一区二区三区| 国产高清国产精品国产三级| 国产精品蜜桃在线观看| 老司机影院毛片| 精品久久久久久电影网| 热99国产精品久久久久久7| tube8黄色片| 国产精品久久久久久久电影| 九色成人免费人妻av| 麻豆成人午夜福利视频| 97超碰精品成人国产| 精品人妻偷拍中文字幕| 国产日韩一区二区三区精品不卡 | 纯流量卡能插随身wifi吗| 精品人妻熟女毛片av久久网站| 亚洲精品,欧美精品| 成人影院久久| 国产精品久久久久久久久免| 免费大片黄手机在线观看| 最黄视频免费看| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放 | 一区二区三区乱码不卡18| 国产深夜福利视频在线观看| 国产免费福利视频在线观看| 91久久精品国产一区二区三区| 日本欧美国产在线视频| 久热这里只有精品99| 美女大奶头黄色视频| 一级片'在线观看视频| 国内精品宾馆在线| 午夜视频国产福利| 自拍偷自拍亚洲精品老妇| 在线观看国产h片| 日韩一区二区三区影片| av有码第一页| 国产极品粉嫩免费观看在线 | av女优亚洲男人天堂| 久久国产精品男人的天堂亚洲 | 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看| 免费观看的影片在线观看| 亚洲av成人精品一区久久| 大话2 男鬼变身卡| 国产视频首页在线观看| 中文字幕亚洲精品专区| 免费观看在线日韩| 成人18禁高潮啪啪吃奶动态图 | 久久精品国产a三级三级三级| 久久久国产精品麻豆| 69精品国产乱码久久久| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 看非洲黑人一级黄片| 免费观看无遮挡的男女| 日韩视频在线欧美| 成人毛片60女人毛片免费| 夫妻午夜视频| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 少妇丰满av| 男男h啪啪无遮挡| av黄色大香蕉| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 久久久久久久久久久丰满| 插逼视频在线观看| 男女国产视频网站| 三上悠亚av全集在线观看 | 国产精品99久久99久久久不卡 | 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 自线自在国产av| 夫妻性生交免费视频一级片| 欧美少妇被猛烈插入视频| 亚洲精品成人av观看孕妇| 亚洲精华国产精华液的使用体验| 亚洲精品日韩av片在线观看| 91精品伊人久久大香线蕉| 99久久人妻综合| 人体艺术视频欧美日本| 夫妻午夜视频| 国产一级毛片在线| 欧美三级亚洲精品| 精品一区在线观看国产| 另类精品久久| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 热99国产精品久久久久久7| 国产av一区二区精品久久| 国产精品国产三级专区第一集| 人人澡人人妻人| 狠狠精品人妻久久久久久综合| 夫妻性生交免费视频一级片| 欧美少妇被猛烈插入视频| 热re99久久精品国产66热6| 亚洲国产日韩一区二区| 各种免费的搞黄视频| 亚洲欧美中文字幕日韩二区| 人人妻人人爽人人添夜夜欢视频 | 色视频www国产| 成年人免费黄色播放视频 | 国产精品蜜桃在线观看| 国产成人精品无人区| 一级毛片电影观看| 亚洲精品日本国产第一区| 内地一区二区视频在线| 亚洲国产精品一区二区三区在线| 人人妻人人澡人人看| 少妇猛男粗大的猛烈进出视频| 久久久国产欧美日韩av| 国产午夜精品久久久久久一区二区三区| 黄色毛片三级朝国网站 | 国产高清国产精品国产三级| 精品人妻一区二区三区麻豆| 亚洲av中文av极速乱| 天美传媒精品一区二区| 久久精品国产a三级三级三级| 制服丝袜香蕉在线| 国产黄片视频在线免费观看| 精品视频人人做人人爽| 精华霜和精华液先用哪个| 美女内射精品一级片tv| 蜜桃在线观看..| 91午夜精品亚洲一区二区三区| 大码成人一级视频| 免费不卡的大黄色大毛片视频在线观看| 日日啪夜夜撸| 国产精品福利在线免费观看| 国产精品三级大全| 国产精品嫩草影院av在线观看| 搡老乐熟女国产| 午夜91福利影院| 精品亚洲成a人片在线观看| 免费av不卡在线播放| 最黄视频免费看| 国产亚洲最大av| 春色校园在线视频观看| 18+在线观看网站| 日韩伦理黄色片| 亚洲欧美清纯卡通| 国产欧美亚洲国产| xxx大片免费视频| 天堂俺去俺来也www色官网| 精品酒店卫生间| 大码成人一级视频| 老熟女久久久| 国产精品成人在线| 精品久久久精品久久久| 久久久欧美国产精品| 26uuu在线亚洲综合色| 女性被躁到高潮视频| 欧美激情国产日韩精品一区| 国产亚洲欧美精品永久| 99热这里只有精品一区| 人妻夜夜爽99麻豆av| 丰满迷人的少妇在线观看| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| av国产精品久久久久影院| 精品一区二区三卡| 日韩不卡一区二区三区视频在线| 在线免费观看不下载黄p国产| www.av在线官网国产| 久久久久视频综合| 国产爽快片一区二区三区| 成人毛片60女人毛片免费| 日日爽夜夜爽网站| 午夜福利,免费看| 一本色道久久久久久精品综合| 麻豆乱淫一区二区| 你懂的网址亚洲精品在线观看| 久久久久人妻精品一区果冻| 久久久a久久爽久久v久久| 国产精品偷伦视频观看了| 一级毛片我不卡| 午夜福利视频精品| 亚洲人成网站在线观看播放| 自拍欧美九色日韩亚洲蝌蚪91 | 男女边摸边吃奶| 久久99热这里只频精品6学生| 久久久国产欧美日韩av| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| av.在线天堂| 18禁在线播放成人免费| 成人午夜精彩视频在线观看| av又黄又爽大尺度在线免费看| 久久久久久久大尺度免费视频| 交换朋友夫妻互换小说| 春色校园在线视频观看| 欧美最新免费一区二区三区| 日本av手机在线免费观看| 精品视频人人做人人爽| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 欧美精品一区二区免费开放| 亚洲国产色片| 亚洲第一区二区三区不卡| 久久99热这里只频精品6学生| 精品一品国产午夜福利视频| 国产成人91sexporn| 欧美3d第一页| 日产精品乱码卡一卡2卡三| 午夜福利影视在线免费观看| 熟女av电影| 亚洲av.av天堂| 国产日韩欧美视频二区| 久久久亚洲精品成人影院| 人妻人人澡人人爽人人| 久久女婷五月综合色啪小说| 中文字幕av电影在线播放| 久久毛片免费看一区二区三区| 国产成人91sexporn| 国产精品伦人一区二区| 这个男人来自地球电影免费观看 | 日本欧美视频一区| 久久99蜜桃精品久久| 男人爽女人下面视频在线观看| 女性生殖器流出的白浆| 亚洲av电影在线观看一区二区三区| 99热这里只有是精品50| 最近最新中文字幕免费大全7| 久热久热在线精品观看| 国产精品无大码| 男人爽女人下面视频在线观看| .国产精品久久| 少妇裸体淫交视频免费看高清| 久久婷婷青草| 啦啦啦在线观看免费高清www| 乱人伦中国视频| 少妇人妻一区二区三区视频| 人妻系列 视频| 我的老师免费观看完整版| 久久韩国三级中文字幕| 少妇高潮的动态图| 啦啦啦中文免费视频观看日本| 人妻 亚洲 视频| 中国三级夫妇交换| 亚洲不卡免费看| 欧美精品一区二区免费开放| 色视频www国产| 亚洲精品aⅴ在线观看| 成人无遮挡网站| 精品人妻熟女毛片av久久网站| 精品一区在线观看国产| 精品视频人人做人人爽| 97超视频在线观看视频| 99久久人妻综合| 婷婷色综合大香蕉| 香蕉精品网在线| 这个男人来自地球电影免费观看 | 两个人免费观看高清视频 | 国产又色又爽无遮挡免| 黑人猛操日本美女一级片| 一级爰片在线观看| 亚洲自偷自拍三级| 久久久久久久亚洲中文字幕| 国产一区二区在线观看av| 国产精品伦人一区二区| 日本色播在线视频| 免费看日本二区| 久久久午夜欧美精品| av不卡在线播放| 美女脱内裤让男人舔精品视频| 亚洲,一卡二卡三卡| 亚洲国产精品999| 久久久精品94久久精品| 国产午夜精品一二区理论片| 一级毛片电影观看| 亚洲丝袜综合中文字幕| 精品久久久久久久久av| 国产欧美日韩精品一区二区| 国国产精品蜜臀av免费| 国产色婷婷99| 黑人高潮一二区| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 欧美精品亚洲一区二区| 亚洲精品乱码久久久久久按摩| 亚洲欧美中文字幕日韩二区| 天天躁夜夜躁狠狠久久av| 在现免费观看毛片| 国产精品一区二区性色av| 我要看黄色一级片免费的| 日韩熟女老妇一区二区性免费视频| 91在线精品国自产拍蜜月| 久久亚洲国产成人精品v| 91久久精品国产一区二区三区| 日韩视频在线欧美| 日本猛色少妇xxxxx猛交久久| 一级爰片在线观看| 高清av免费在线| 欧美精品高潮呻吟av久久| www.色视频.com| 亚洲人成网站在线播| 午夜av观看不卡| a级毛色黄片| 日本黄大片高清| 国精品久久久久久国模美| 国产色婷婷99| 欧美日韩精品成人综合77777| 亚洲va在线va天堂va国产| 免费大片黄手机在线观看| 午夜久久久在线观看| 97精品久久久久久久久久精品| 亚洲欧美日韩东京热| 美女xxoo啪啪120秒动态图| kizo精华| 偷拍熟女少妇极品色| 尾随美女入室| 久久青草综合色| 亚洲性久久影院| 亚洲情色 制服丝袜| 亚洲国产精品999| 偷拍熟女少妇极品色| 国产成人一区二区在线| 午夜免费鲁丝| 女人久久www免费人成看片| 一二三四中文在线观看免费高清| 国产精品一区二区三区四区免费观看| 六月丁香七月| 国产有黄有色有爽视频| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 热99国产精品久久久久久7| 男女免费视频国产| 最近2019中文字幕mv第一页| 国产极品天堂在线| 免费看日本二区| 少妇被粗大猛烈的视频| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看| 少妇精品久久久久久久| 国产男女内射视频| kizo精华| 国产av国产精品国产| 黄色一级大片看看| 精品卡一卡二卡四卡免费| .国产精品久久| 一级毛片我不卡| 中文字幕人妻丝袜制服| 国产老妇伦熟女老妇高清| 成人毛片a级毛片在线播放| 简卡轻食公司| 在线观看免费日韩欧美大片 | 99re6热这里在线精品视频| 日本av手机在线免费观看| 免费大片黄手机在线观看| 久久韩国三级中文字幕| 中文字幕人妻丝袜制服| 国产成人aa在线观看| 这个男人来自地球电影免费观看 | 夜夜骑夜夜射夜夜干| 日本-黄色视频高清免费观看| 国产精品嫩草影院av在线观看| 欧美精品一区二区大全| 麻豆乱淫一区二区| av福利片在线| 狠狠精品人妻久久久久久综合| 日本wwww免费看| 伊人久久国产一区二区| 一级片'在线观看视频| 日韩精品免费视频一区二区三区 | 久久99精品国语久久久| 麻豆乱淫一区二区| 久久精品国产自在天天线| 亚洲av在线观看美女高潮| 18禁裸乳无遮挡动漫免费视频| 成人美女网站在线观看视频| 亚洲精品国产色婷婷电影| 亚洲国产精品999| av在线播放精品| 夫妻性生交免费视频一级片| 老司机影院成人| 日韩一本色道免费dvd| 99热这里只有是精品在线观看| 亚洲天堂av无毛| 久久精品国产亚洲av天美| 18禁裸乳无遮挡动漫免费视频| 日韩电影二区| 人人澡人人妻人| 五月开心婷婷网| 欧美日韩国产mv在线观看视频| 制服丝袜香蕉在线| 三级国产精品片| 婷婷色综合大香蕉| 高清在线视频一区二区三区| 国产男女内射视频| 99国产精品免费福利视频| 啦啦啦视频在线资源免费观看| 成人亚洲欧美一区二区av| 一本一本综合久久| 欧美亚洲 丝袜 人妻 在线| 免费看日本二区| 国产一区二区在线观看av| 亚洲第一av免费看| 亚洲精品国产色婷婷电影| 国产成人freesex在线| 精品国产国语对白av| 国产熟女欧美一区二区| 狂野欧美激情性xxxx在线观看| 国产精品不卡视频一区二区| 久久韩国三级中文字幕| 久久99一区二区三区| 在线免费观看不下载黄p国产| 亚洲熟女精品中文字幕| 日本91视频免费播放| 国产精品免费大片| 十八禁网站网址无遮挡 | 亚洲激情五月婷婷啪啪| 丝袜在线中文字幕| 国产精品国产三级国产av玫瑰| 成人美女网站在线观看视频| 欧美性感艳星| 大陆偷拍与自拍| 人妻少妇偷人精品九色| 深夜a级毛片| 伊人久久国产一区二区| 男女边吃奶边做爰视频| 人妻制服诱惑在线中文字幕| 大香蕉97超碰在线| 蜜臀久久99精品久久宅男| 18禁动态无遮挡网站| 七月丁香在线播放| 国产精品国产av在线观看| 久久精品久久精品一区二区三区| av.在线天堂| 少妇熟女欧美另类| 最新中文字幕久久久久| 久久国产精品男人的天堂亚洲 | 精品久久国产蜜桃| 久久久久国产网址| 在现免费观看毛片| 街头女战士在线观看网站| 亚洲不卡免费看| 午夜激情久久久久久久| 一个人看视频在线观看www免费| 午夜老司机福利剧场| 男的添女的下面高潮视频| 免费看日本二区| 最近最新中文字幕免费大全7| 另类亚洲欧美激情| 嫩草影院入口| 日韩av不卡免费在线播放| 欧美 日韩 精品 国产| 爱豆传媒免费全集在线观看| 国产伦在线观看视频一区| 丰满少妇做爰视频| 久久免费观看电影| 久久久久久人妻| 日韩av免费高清视频| 亚洲在久久综合| 亚洲va在线va天堂va国产| 日日啪夜夜撸| 又爽又黄a免费视频| 久久久精品94久久精品| 亚洲av二区三区四区| 乱系列少妇在线播放| 精品久久久久久久久亚洲| 欧美精品高潮呻吟av久久| 一级二级三级毛片免费看| .国产精品久久| 狂野欧美激情性bbbbbb| 亚洲欧洲国产日韩| 精品一区二区三卡| 久久97久久精品| 人妻夜夜爽99麻豆av| 嘟嘟电影网在线观看| 免费不卡的大黄色大毛片视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲综合精品二区| 亚洲美女黄色视频免费看| 国产色爽女视频免费观看| 美女xxoo啪啪120秒动态图| 插逼视频在线观看| 色视频www国产| 一本—道久久a久久精品蜜桃钙片| 91精品一卡2卡3卡4卡| 精品少妇内射三级| 亚洲精品日韩在线中文字幕| 丝瓜视频免费看黄片| 极品人妻少妇av视频| videos熟女内射| 多毛熟女@视频| h日本视频在线播放| 日韩免费高清中文字幕av| 永久网站在线| 黄色视频在线播放观看不卡| 精品国产乱码久久久久久小说| 精品视频人人做人人爽| 卡戴珊不雅视频在线播放| 日韩人妻高清精品专区| 久久精品国产鲁丝片午夜精品| 一级毛片电影观看| 大陆偷拍与自拍| 天堂俺去俺来也www色官网| 亚洲一级一片aⅴ在线观看| 国产伦在线观看视频一区| 欧美 亚洲 国产 日韩一| 日本色播在线视频| 在线观看免费视频网站a站| 91久久精品国产一区二区成人| 欧美精品人与动牲交sv欧美| 久久精品国产亚洲av天美| 国产高清国产精品国产三级| 高清午夜精品一区二区三区| 免费看日本二区| 日韩av免费高清视频| 三级经典国产精品| 久久精品夜色国产| 一级毛片 在线播放| 欧美人与善性xxx| 久久久久网色| 久久精品久久久久久噜噜老黄| av视频免费观看在线观看| 欧美日韩国产mv在线观看视频| 自拍偷自拍亚洲精品老妇| 午夜影院在线不卡| 中文字幕av电影在线播放| 日韩av不卡免费在线播放| av在线app专区| 视频中文字幕在线观看| 交换朋友夫妻互换小说| 久久精品国产鲁丝片午夜精品| 欧美性感艳星| 国产免费一级a男人的天堂| 熟妇人妻不卡中文字幕| 中国国产av一级| 亚洲国产精品一区二区三区在线| 亚洲人成网站在线播| 亚洲精品国产色婷婷电影| 亚洲电影在线观看av| 美女视频免费永久观看网站| 国产亚洲午夜精品一区二区久久| 久久国产精品大桥未久av | 一级毛片 在线播放| 伊人久久精品亚洲午夜| 日韩 亚洲 欧美在线| 日本vs欧美在线观看视频 | 美女cb高潮喷水在线观看| 69精品国产乱码久久久| 赤兔流量卡办理| 免费观看的影片在线观看| 精品国产乱码久久久久久小说| 18+在线观看网站| 亚洲四区av| 欧美bdsm另类| 国产午夜精品一二区理论片| 99久久精品国产国产毛片| 国产欧美日韩综合在线一区二区 | 久久久国产一区二区| 性高湖久久久久久久久免费观看| 亚洲欧美日韩另类电影网站| 五月玫瑰六月丁香| 色94色欧美一区二区| 另类精品久久| 97在线视频观看| 熟女电影av网| 三上悠亚av全集在线观看 | 久久久久久久久久人人人人人人| 国产69精品久久久久777片| 有码 亚洲区| 欧美区成人在线视频|