• <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)變化
    99热这里只有是精品50| 99久久中文字幕三级久久日本| 日本午夜av视频| av在线播放精品| 国产高清有码在线观看视频| 简卡轻食公司| 免费看日本二区| 国产成人aa在线观看| 91精品伊人久久大香线蕉| 国产 一区 欧美 日韩| 我要看日韩黄色一级片| 成人毛片60女人毛片免费| 中文字幕人妻熟人妻熟丝袜美| 99久国产av精品国产电影| www.色视频.com| 干丝袜人妻中文字幕| 国产av在哪里看| 啦啦啦韩国在线观看视频| 色视频www国产| 99久久无色码亚洲精品果冻| 久久99热这里只频精品6学生 | 啦啦啦啦在线视频资源| a级毛片免费高清观看在线播放| 夜夜看夜夜爽夜夜摸| 国产成人a区在线观看| 久久午夜福利片| 日本免费一区二区三区高清不卡| 老师上课跳d突然被开到最大视频| 久久久精品大字幕| 黄色配什么色好看| 久久草成人影院| 亚州av有码| 国产一区有黄有色的免费视频 | 亚洲国产精品合色在线| 国产精品.久久久| 91午夜精品亚洲一区二区三区| 国产精品无大码| 伦理电影大哥的女人| 黄色日韩在线| videos熟女内射| 亚洲综合精品二区| 精品午夜福利在线看| 长腿黑丝高跟| 中文乱码字字幕精品一区二区三区 | 亚洲真实伦在线观看| 成人亚洲欧美一区二区av| 2021少妇久久久久久久久久久| 天堂√8在线中文| 天天躁夜夜躁狠狠久久av| 久久久久久久亚洲中文字幕| 久久精品国产99精品国产亚洲性色| 成人美女网站在线观看视频| 老司机影院毛片| 26uuu在线亚洲综合色| 九色成人免费人妻av| 成年版毛片免费区| 中文在线观看免费www的网站| 少妇猛男粗大的猛烈进出视频 | 亚洲国产最新在线播放| 国产精品麻豆人妻色哟哟久久 | 色综合色国产| 日韩在线高清观看一区二区三区| 一本久久精品| 国产精品无大码| www.av在线官网国产| 黄色欧美视频在线观看| 国产精品久久久久久精品电影小说 | 人妻系列 视频| 国产三级在线视频| 美女脱内裤让男人舔精品视频| 天堂影院成人在线观看| 国产精品伦人一区二区| 黑人高潮一二区| 伦理电影大哥的女人| 午夜精品一区二区三区免费看| 不卡视频在线观看欧美| 亚洲欧美精品专区久久| 久久国产乱子免费精品| 午夜福利高清视频| 亚洲成av人片在线播放无| 在线免费观看不下载黄p国产| av免费观看日本| 亚洲av二区三区四区| 国产熟女欧美一区二区| 亚洲18禁久久av| 97热精品久久久久久| 啦啦啦啦在线视频资源| 午夜爱爱视频在线播放| 免费电影在线观看免费观看| 男女啪啪激烈高潮av片| 国产精品1区2区在线观看.| 国产精品久久久久久久久免| 1024手机看黄色片| 精品国内亚洲2022精品成人| 22中文网久久字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大香蕉久久网| 直男gayav资源| 中文字幕精品亚洲无线码一区| 国产一区亚洲一区在线观看| 欧美精品一区二区大全| 人人妻人人澡欧美一区二区| 久久久久久国产a免费观看| 日韩亚洲欧美综合| 人人妻人人看人人澡| 欧美+日韩+精品| 国产精品久久久久久久久免| 在线观看美女被高潮喷水网站| 久久久久久大精品| 色综合亚洲欧美另类图片| ponron亚洲| 最近最新中文字幕免费大全7| 日韩高清综合在线| 免费观看精品视频网站| 国产亚洲5aaaaa淫片| 少妇丰满av| 久久久久久九九精品二区国产| 联通29元200g的流量卡| 内地一区二区视频在线| 亚洲人与动物交配视频| 一级毛片我不卡| 高清毛片免费看| 久久鲁丝午夜福利片| 日日啪夜夜撸| 卡戴珊不雅视频在线播放| 久久99热这里只频精品6学生 | 九九在线视频观看精品| 看免费成人av毛片| 18+在线观看网站| videos熟女内射| 亚洲色图av天堂| 午夜a级毛片| 亚洲欧美日韩卡通动漫| 嘟嘟电影网在线观看| 午夜视频国产福利| 久久久久久久久久成人| av专区在线播放| 国产精品电影一区二区三区| 三级毛片av免费| 舔av片在线| av播播在线观看一区| 免费大片18禁| 国产在线一区二区三区精 | 国产av在哪里看| 日本wwww免费看| 成人亚洲欧美一区二区av| 麻豆成人av视频| 少妇熟女aⅴ在线视频| kizo精华| 久久久久久久久久久免费av| 国产一区二区在线观看日韩| 一边摸一边抽搐一进一小说| 精品久久久久久成人av| 国产日韩欧美在线精品| 欧美成人午夜免费资源| 如何舔出高潮| 日本三级黄在线观看| 久久久精品大字幕| 成人午夜精彩视频在线观看| 丝袜喷水一区| 不卡视频在线观看欧美| 亚洲,欧美,日韩| 韩国高清视频一区二区三区| 欧美xxxx性猛交bbbb| 午夜激情欧美在线| 国产大屁股一区二区在线视频| 精品久久久久久久人妻蜜臀av| 中文字幕久久专区| 国产精品电影一区二区三区| 99热精品在线国产| 成人鲁丝片一二三区免费| 久久草成人影院| 亚洲av成人精品一区久久| ponron亚洲| 成人漫画全彩无遮挡| 久久久国产成人精品二区| 午夜久久久久精精品| 欧美另类亚洲清纯唯美| 久久99热6这里只有精品| 人妻系列 视频| 成人性生交大片免费视频hd| 在现免费观看毛片| 夜夜爽夜夜爽视频| 欧美日韩在线观看h| 欧美成人a在线观看| 久久久久久久亚洲中文字幕| 国产成人freesex在线| 大又大粗又爽又黄少妇毛片口| 丝袜喷水一区| 国产免费又黄又爽又色| 99热精品在线国产| 久久久久久久亚洲中文字幕| 欧美日韩综合久久久久久| 亚洲18禁久久av| 一级黄片播放器| 日韩国内少妇激情av| 精品久久久久久电影网 | 性色avwww在线观看| 国产免费视频播放在线视频 | 亚洲av一区综合| 男女边吃奶边做爰视频| 国产私拍福利视频在线观看| 久久久久久国产a免费观看| 免费在线观看成人毛片| 最近中文字幕2019免费版| 床上黄色一级片| 青春草国产在线视频| 亚洲最大成人中文| 菩萨蛮人人尽说江南好唐韦庄 | 国产国拍精品亚洲av在线观看| 成人特级av手机在线观看| 国产精品精品国产色婷婷| 视频中文字幕在线观看| 丝袜美腿在线中文| 六月丁香七月| 大又大粗又爽又黄少妇毛片口| 国产av在哪里看| 亚洲高清免费不卡视频| av专区在线播放| 亚洲激情五月婷婷啪啪| 国产精品日韩av在线免费观看| 欧美精品国产亚洲| 在线观看一区二区三区| 一级毛片电影观看 | 亚洲综合精品二区| 日韩亚洲欧美综合| 欧美人与善性xxx| 国产亚洲5aaaaa淫片| 菩萨蛮人人尽说江南好唐韦庄 | 九九久久精品国产亚洲av麻豆| av又黄又爽大尺度在线免费看 | 亚洲高清免费不卡视频| 乱系列少妇在线播放| 成年av动漫网址| 国产 一区 欧美 日韩| a级一级毛片免费在线观看| kizo精华| 亚洲一级一片aⅴ在线观看| 久99久视频精品免费| 欧美人与善性xxx| av免费在线看不卡| 成人鲁丝片一二三区免费| 日韩欧美在线乱码| 日韩av在线免费看完整版不卡| 精品国产露脸久久av麻豆 | 99热这里只有精品一区| 联通29元200g的流量卡| 日日摸夜夜添夜夜添av毛片| 精华霜和精华液先用哪个| 男人的好看免费观看在线视频| 精品久久久久久久久久久久久| 男女那种视频在线观看| 天堂影院成人在线观看| 联通29元200g的流量卡| 99久久九九国产精品国产免费| 亚洲av.av天堂| 亚洲精品乱码久久久久久按摩| www日本黄色视频网| 亚洲精品乱久久久久久| 18禁在线播放成人免费| 亚洲欧美清纯卡通| 久久精品国产99精品国产亚洲性色| 国产精品一区二区在线观看99 | 丰满乱子伦码专区| 一级黄片播放器| videos熟女内射| 国产日韩欧美在线精品| 高清视频免费观看一区二区 | 国产亚洲精品av在线| 国产精品.久久久| 色综合亚洲欧美另类图片| 99热全是精品| 你懂的网址亚洲精品在线观看 | 夫妻性生交免费视频一级片| 国产精品av视频在线免费观看| 国产精品一区二区三区四区久久| 男女边吃奶边做爰视频| 欧美成人一区二区免费高清观看| 久久久国产成人精品二区| 真实男女啪啪啪动态图| 伦精品一区二区三区| 国产精品综合久久久久久久免费| 国产精品日韩av在线免费观看| 久久精品国产亚洲网站| 精品一区二区三区人妻视频| 好男人视频免费观看在线| 亚洲欧美成人精品一区二区| 色噜噜av男人的天堂激情| 国模一区二区三区四区视频| 久久久久久久国产电影| 老司机福利观看| 欧美日本视频| 亚洲精品亚洲一区二区| 国产欧美另类精品又又久久亚洲欧美| 国产精品爽爽va在线观看网站| 久久久久九九精品影院| 内地一区二区视频在线| 2021天堂中文幕一二区在线观| 亚洲最大成人手机在线| 特级一级黄色大片| 日韩精品有码人妻一区| 一区二区三区免费毛片| 午夜老司机福利剧场| 久久人人爽人人片av| 在线播放无遮挡| 两个人的视频大全免费| av又黄又爽大尺度在线免费看 | 在现免费观看毛片| .国产精品久久| 国产高清国产精品国产三级 | 亚洲av电影在线观看一区二区三区 | 特级一级黄色大片| 欧美一区二区精品小视频在线| 2022亚洲国产成人精品| 亚洲最大成人av| 一本久久精品| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 日本免费在线观看一区| 九九久久精品国产亚洲av麻豆| 国模一区二区三区四区视频| 国产综合懂色| 日韩av不卡免费在线播放| 国产成人免费观看mmmm| 久久精品久久精品一区二区三区| 老司机影院成人| 成人综合一区亚洲| 女人久久www免费人成看片 | 插阴视频在线观看视频| 国产精品久久久久久久久免| 欧美人与善性xxx| 99久久精品热视频| 亚洲精品456在线播放app| 成人毛片a级毛片在线播放| 国模一区二区三区四区视频| 久久久久久久久久久免费av| 日日干狠狠操夜夜爽| 国产成人午夜福利电影在线观看| 国产成人福利小说| 啦啦啦观看免费观看视频高清| 久久99热这里只有精品18| 亚洲在久久综合| 亚洲综合色惰| 日本爱情动作片www.在线观看| 九九久久精品国产亚洲av麻豆| 淫秽高清视频在线观看| 女人被狂操c到高潮| 国产伦精品一区二区三区四那| 男的添女的下面高潮视频| 黑人高潮一二区| 国产亚洲精品久久久com| 国产成人精品一,二区| 国产黄片美女视频| 精品久久久久久电影网 | 级片在线观看| 亚洲国产高清在线一区二区三| 精品久久国产蜜桃| 亚洲三级黄色毛片| 亚洲国产欧洲综合997久久,| 国产大屁股一区二区在线视频| 少妇丰满av| 国产免费又黄又爽又色| 亚洲一级一片aⅴ在线观看| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 久久人人爽人人爽人人片va| 国产真实乱freesex| 99久久精品一区二区三区| 麻豆久久精品国产亚洲av| 热99在线观看视频| 国产欧美日韩精品一区二区| 人人妻人人澡欧美一区二区| 一本久久精品| 国产在线男女| 黄片无遮挡物在线观看| av在线播放精品| 日本黄色片子视频| 国产激情偷乱视频一区二区| 国内精品美女久久久久久| 美女国产视频在线观看| 亚洲va在线va天堂va国产| 国产精品永久免费网站| 精品酒店卫生间| 国产欧美另类精品又又久久亚洲欧美| 一级av片app| 男的添女的下面高潮视频| 久久综合国产亚洲精品| 久久久久久久国产电影| 日韩,欧美,国产一区二区三区 | 在线观看美女被高潮喷水网站| 水蜜桃什么品种好| 国产成人精品一,二区| 永久网站在线| 日韩强制内射视频| 干丝袜人妻中文字幕| 成人亚洲欧美一区二区av| 久久精品综合一区二区三区| 狂野欧美白嫩少妇大欣赏| 又黄又爽又刺激的免费视频.| 三级经典国产精品| 久久精品国产自在天天线| 大香蕉97超碰在线| 九草在线视频观看| 99在线人妻在线中文字幕| 一级av片app| h日本视频在线播放| 精品午夜福利在线看| 一区二区三区免费毛片| 蜜桃久久精品国产亚洲av| 人妻制服诱惑在线中文字幕| 亚洲国产精品sss在线观看| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 免费看光身美女| 丰满人妻一区二区三区视频av| 在线免费观看不下载黄p国产| 看十八女毛片水多多多| 成年女人永久免费观看视频| 秋霞伦理黄片| 国产精品爽爽va在线观看网站| 少妇的逼水好多| 久久久久久久午夜电影| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 亚洲最大成人中文| av卡一久久| 国产成人一区二区在线| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 精品人妻一区二区三区麻豆| av福利片在线观看| 黄片无遮挡物在线观看| 最后的刺客免费高清国语| 成人一区二区视频在线观看| 神马国产精品三级电影在线观看| 在线观看美女被高潮喷水网站| av在线老鸭窝| 午夜a级毛片| 亚洲美女搞黄在线观看| 亚洲综合色惰| 亚洲美女搞黄在线观看| 麻豆久久精品国产亚洲av| 国产精品精品国产色婷婷| 99热这里只有精品一区| 插逼视频在线观看| 色综合站精品国产| 在线免费观看不下载黄p国产| 麻豆成人av视频| 久久久久久久久大av| 欧美色视频一区免费| 中文资源天堂在线| 五月玫瑰六月丁香| 国产精品蜜桃在线观看| 赤兔流量卡办理| 久久久成人免费电影| 亚洲最大成人手机在线| 欧美不卡视频在线免费观看| 久久综合国产亚洲精品| 午夜视频国产福利| 又粗又硬又长又爽又黄的视频| av线在线观看网站| 国产乱人视频| 超碰av人人做人人爽久久| av黄色大香蕉| 日韩欧美精品v在线| 国产又色又爽无遮挡免| 亚洲精品,欧美精品| 久久精品国产99精品国产亚洲性色| 少妇人妻一区二区三区视频| 熟女人妻精品中文字幕| 亚洲电影在线观看av| 噜噜噜噜噜久久久久久91| 国产精品永久免费网站| 日韩精品青青久久久久久| 国产在视频线精品| 亚洲,欧美,日韩| 亚洲精品色激情综合| 一级av片app| 亚洲电影在线观看av| 欧美三级亚洲精品| 精品久久国产蜜桃| 日韩一区二区视频免费看| 亚洲最大成人av| 久久精品国产亚洲av涩爱| 少妇裸体淫交视频免费看高清| 欧美激情久久久久久爽电影| 一本久久精品| 久久热精品热| 亚洲一区高清亚洲精品| 久久这里只有精品中国| 国产成人精品一,二区| 男女边吃奶边做爰视频| 国产精品99久久久久久久久| 91久久精品国产一区二区成人| 欧美日韩综合久久久久久| 嫩草影院入口| 看免费成人av毛片| 亚洲av成人精品一二三区| 亚洲欧美中文字幕日韩二区| 女人久久www免费人成看片 | 久久久色成人| 久久精品国产亚洲av涩爱| 国产高清视频在线观看网站| 中文字幕免费在线视频6| 三级国产精品欧美在线观看| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| 国产成人a∨麻豆精品| 日本免费在线观看一区| 一本一本综合久久| 欧美区成人在线视频| 乱人视频在线观看| 久久精品国产鲁丝片午夜精品| 身体一侧抽搐| 欧美zozozo另类| 麻豆成人av视频| 久久久久久久午夜电影| 亚洲av成人精品一区久久| 亚洲av成人精品一二三区| 国产中年淑女户外野战色| 色网站视频免费| a级毛片免费高清观看在线播放| 观看免费一级毛片| 亚洲欧美日韩无卡精品| 国产亚洲精品久久久com| 天天一区二区日本电影三级| 日本猛色少妇xxxxx猛交久久| 国产大屁股一区二区在线视频| 精品午夜福利在线看| 97超碰精品成人国产| 一二三四中文在线观看免费高清| 亚洲中文字幕日韩| 秋霞在线观看毛片| 亚洲av成人精品一二三区| 黄色欧美视频在线观看| 国产av不卡久久| 别揉我奶头 嗯啊视频| av国产免费在线观看| 免费在线观看成人毛片| 国产精品久久久久久久电影| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 爱豆传媒免费全集在线观看| 婷婷色av中文字幕| 婷婷色麻豆天堂久久 | 九草在线视频观看| 最近中文字幕高清免费大全6| 日韩国内少妇激情av| 国语对白做爰xxxⅹ性视频网站| 最近中文字幕2019免费版| 91精品伊人久久大香线蕉| 亚洲无线观看免费| 亚洲熟妇中文字幕五十中出| 国产精品久久电影中文字幕| 五月伊人婷婷丁香| 麻豆乱淫一区二区| 国产极品天堂在线| 国产精品一区www在线观看| 纵有疾风起免费观看全集完整版 | 男的添女的下面高潮视频| 日韩高清综合在线| 三级男女做爰猛烈吃奶摸视频| 校园人妻丝袜中文字幕| 美女被艹到高潮喷水动态| 亚洲综合色惰| 欧美成人免费av一区二区三区| 国产极品天堂在线| 51国产日韩欧美| 边亲边吃奶的免费视频| 中文乱码字字幕精品一区二区三区 | 精品不卡国产一区二区三区| 欧美xxxx性猛交bbbb| 男女下面进入的视频免费午夜| 免费观看精品视频网站| 中文字幕熟女人妻在线| 日韩在线高清观看一区二区三区| 九九热线精品视视频播放| 黄色欧美视频在线观看| 美女黄网站色视频| 日韩欧美国产在线观看| 成人国产麻豆网| 老司机福利观看| 热99在线观看视频| 特大巨黑吊av在线直播| 黄片无遮挡物在线观看| 国产精品一二三区在线看| 日日啪夜夜撸| 三级国产精品片| 少妇人妻一区二区三区视频| 真实男女啪啪啪动态图| 超碰av人人做人人爽久久| 亚洲中文字幕日韩| 麻豆av噜噜一区二区三区| 国产精品久久久久久精品电影| 久久久久久久久中文| 丝袜美腿在线中文| 夫妻性生交免费视频一级片| av国产久精品久网站免费入址| 国产毛片a区久久久久| 岛国在线免费视频观看| 日韩国内少妇激情av| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放 | 久久久a久久爽久久v久久| 精品久久久久久电影网 | 亚洲欧美成人综合另类久久久 | 亚洲av免费在线观看| 国产精品99久久久久久久久| 亚洲av免费高清在线观看| av免费观看日本| 精品国内亚洲2022精品成人| 亚洲aⅴ乱码一区二区在线播放| 国内少妇人妻偷人精品xxx网站| 精品久久久久久成人av| a级毛片免费高清观看在线播放|