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

    基于寒害過程的福建芒果種植氣候風(fēng)險區(qū)劃*

    2019-11-15 12:36:24陳彬彬李麗純陳斌源
    中國農(nóng)業(yè)氣象 2019年11期
    關(guān)鍵詞:區(qū)劃芒果氣候

    楊 凱,陳彬彬,陳 惠,李麗純,陳斌源,李 政

    基于寒害過程的福建芒果種植氣候風(fēng)險區(qū)劃*

    楊 凱1,陳彬彬1,陳 惠1,李麗純1,陳斌源2,李 政3

    (1.福建省氣象科學(xué)研究所,福州 350001;2.漳州市氣象局,漳州 363000;3.廣西壯族自治區(qū)氣象減災(zāi)研究所,南寧 530022)

    通過冬季寒害災(zāi)情調(diào)查、果園定位觀測、地理移置和人工氣候箱低溫致災(zāi)試驗4種渠道,收集芒果寒害的災(zāi)情樣本39例,從而確定芒果寒害臨界溫度為2.0℃。在此基礎(chǔ)上,篩選出福建芒果主產(chǎn)縣市2001?2014年冬季28個寒害過程,通過將寒害過程溫度、降水等氣象因子與相對氣象產(chǎn)量進行相關(guān)分析,確定芒果寒害致災(zāi)因子為過程極端最低氣溫、日最低氣溫≤2.0℃有害積寒、≤2.0℃寒害過程持續(xù)日數(shù)、≤2.0℃寒害過程降溫幅度。利用主成分分析對4個致災(zāi)因子進行綜合簡化,得到芒果寒害評價的綜合氣候指標(biāo),結(jié)合相對氣象產(chǎn)量確定了指標(biāo)分級。根據(jù)寒害綜合氣候指標(biāo)等級,應(yīng)用GIS方法完成福建省芒果寒害種植氣候風(fēng)險區(qū)劃,無寒害?輕度寒害的區(qū)域主要分布在沿海地區(qū),是引種芒果優(yōu)質(zhì)品種寒害風(fēng)險最低的區(qū)域。

    芒果;寒害;臨界溫度;綜合氣候指標(biāo);風(fēng)險區(qū)劃

    芒果原產(chǎn)印度,是漆樹科熱帶常綠大喬木[1],其早結(jié)豐產(chǎn),種植效益高,經(jīng)過多年的引種、試種、擴種,已發(fā)展成為中國主要的特色熱帶水果,2015年中國芒果種植面積位居世界第二位,占17%[2]。福建省是芒果的主產(chǎn)地之一,但是受季風(fēng)氣候影響,時常遭受低溫寒害、臺風(fēng)、干旱等氣象災(zāi)害的影響,其中低溫寒害是影響芒果生產(chǎn)最主要的氣象災(zāi)害。近年來,隨著臺灣優(yōu)質(zhì)芒果品種臺農(nóng)1號、金煌芒等在福建的引種擴種,芒果種植大戶和種植基地迅速增加,優(yōu)勢產(chǎn)業(yè)帶逐步形成。福建省地處熱帶果樹種植的北緣地帶,在全球變暖的大背景下,冬季極端低溫事件頻發(fā),導(dǎo)致芒果遭受寒害[3?5],例如2016年1月22?26日的強寒潮過程,造成寧德市芒果損失嚴(yán)重,部分種植園甚至絕收。加之部分果農(nóng)在不適宜區(qū)盲目引種和擴種,進一步加重了災(zāi)害損失,嚴(yán)重制約了芒果產(chǎn)業(yè)的健康發(fā)展。因此,開展福建省芒果低溫寒害風(fēng)險評估和區(qū)劃,對防寒減災(zāi)決策制定、農(nóng)業(yè)生產(chǎn)布局和區(qū)域可持續(xù)發(fā)展具有重要意義,而合理確定芒果低溫寒害等級指標(biāo)是低溫寒害評估和區(qū)劃的前提。

    國內(nèi)對果樹低溫寒害的風(fēng)險評估及區(qū)劃研究已取得一定進展。劉錦鑾等[6?7]對荔枝、香蕉冬季寒害的災(zāi)損性、易災(zāi)性及防寒抗災(zāi)能力等方面進行分析,但并未確定對果樹寒害有較大影響的具體低溫過程;在果樹寒害等級指標(biāo)研究方面,有學(xué)者針對冬季極端最低氣溫單個因子,對果樹遭受不同等級寒害的極端低溫閾值進行了劃分,但未考慮低溫強度、低溫持續(xù)時間等影響芒果受害程度的因子[8?11]。杜堯東等[12?15]對香蕉、荔枝等開展了寒害綜合氣候指標(biāo)與風(fēng)險區(qū)劃研究,但針對近年來引種的新特優(yōu)熱帶果樹的寒害指標(biāo)研究和風(fēng)險區(qū)劃較少,特別是針對芒果寒害風(fēng)險區(qū)劃方面的研究則更少。本研究以芒果寒害過程為研究對象,在冬季生長季尚未結(jié)束時,針對已發(fā)生的寒害過程對果樹的影響程度進行評估,以期為寒害的動態(tài)評估提供理論基礎(chǔ)。隨著氣候變化、品種更替、新的農(nóng)業(yè)生產(chǎn)技術(shù)的應(yīng)用,沿用舊的寒害指標(biāo)易導(dǎo)致寒害風(fēng)險評估及區(qū)劃出現(xiàn)偏差,造成生產(chǎn)損失。為此,本研究通過對芒果低溫寒害過程的氣象因子進行分析,確定芒果寒害臨界溫度和致災(zāi)因子,構(gòu)建客觀、定量的低溫寒害綜合氣候指標(biāo),得到福建省芒果寒害種植氣候風(fēng)險區(qū)劃,以期為芒果低溫寒害的評估和區(qū)劃提供依據(jù),為農(nóng)業(yè)生產(chǎn)部門制定合理的引種擴種、防災(zāi)減災(zāi)決策和措施提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 資料

    福建省長泰縣、南安市、安溪縣、仙游縣、福清縣、寧德市區(qū)6個芒果主產(chǎn)縣市2001?2014年芒果產(chǎn)量數(shù)據(jù),來源于歷年《福建省農(nóng)村統(tǒng)計年鑒》;6個主產(chǎn)縣市地面氣象觀測臺站2001?2014年冬季(12月?翌年2月)逐日平均氣溫、逐日最低氣溫、逐日降水量資料以及福鼎、永泰、羅源、長泰、平和、龍海、云霄7個縣氣象局觀測場2014/2015年冬季和2015/2016年冬季逐日最低氣溫資料由福建省氣象信息中心提供;芒果發(fā)生寒害所對應(yīng)的當(dāng)時當(dāng)?shù)氐臉O端最低氣溫,來自芒果果園內(nèi)自動氣象站或從鄰近的區(qū)域自動氣象站數(shù)據(jù)訂正獲得。

    1.2 研究思路

    收集芒果寒害災(zāi)情樣本32例,通過分析各樣本極端最低溫度下芒果的受害癥狀,確定導(dǎo)致芒果遭受寒害的臨界溫度,即當(dāng)冬季極端最低溫度達到某個閾值時,芒果開始受害,高于此閾值時寒害結(jié)束。

    第二步,根據(jù)芒果的寒害臨界溫度,篩選出6個主產(chǎn)縣2001?2014年冬季可能造成芒果寒害的過程資料,將過程極端最低溫度、過程平均溫度、過程有害積寒、過程持續(xù)時間、過程降溫幅度、過程降水量與當(dāng)年相對氣象產(chǎn)量進行相關(guān)分析,確定芒果寒害的致災(zāi)因子,分析因子間的相關(guān)性,并采用主成分分析方法構(gòu)建芒果寒害綜合氣候指標(biāo)。

    第三步,參照農(nóng)業(yè)上劃分災(zāi)害年型的方法,結(jié)合相對氣象產(chǎn)量減產(chǎn)率對芒果寒害綜合氣候指標(biāo)進行分級,并利用綜合氣候指標(biāo)閾值進行福建芒果寒害種植氣候風(fēng)險區(qū)劃。

    1.3 芒果寒害災(zāi)情的收集

    1.3.1 田間災(zāi)情調(diào)查

    在2014/2015年冬季和2015/2016年冬季冷空氣過程后,對福建芒果園開展寒害災(zāi)情調(diào)查。根據(jù)《芒果寒害災(zāi)情調(diào)查表》,記錄調(diào)查地點芒果遭受寒害的時間、詳細災(zāi)情及災(zāi)后果樹枝、干、葉的變化,獲得有效寒害災(zāi)情樣本6例。

    1.3.2 果園定位觀測

    2014年10月?2016年3月,在福建南安、南靖、三都澳及廣西田林、陽圩、武鳴六地的芒果易災(zāi)果園和平行觀測點,分別布設(shè)1套自動小氣候觀測系統(tǒng)(RR-9140型,北京),實時獲取溫度、濕度、降水、光照數(shù)據(jù),并用于寒害監(jiān)測;分別安裝1臺紅外自動夜視監(jiān)測儀(SG-008型,中國),設(shè)置每日7:00?19:00每隔10min自動拍照記錄芒果果樹、果苗生長發(fā)育狀況,了解其形態(tài)變化。共獲得有效災(zāi)情樣本18例。

    1.3.3 地理移置試驗

    試驗主要在福建省開展。福建省引種的芒果優(yōu)質(zhì)品種主要為臺農(nóng)1號,因此本試驗采用盆栽方式培育臺農(nóng)1號果苗,待果苗植株健壯后進行地理移置。根據(jù)福建省歷史年極端最低氣溫平均值的分布規(guī)律,同時為了獲得芒果遭受寒害后最準(zhǔn)確氣象數(shù)據(jù),分別于2014/2015年和2015/2016年10月底,提前將芒果果苗3~4盆分別安放在福鼎、永泰、羅源、長泰、平和、龍海、云霄7個縣氣象局觀測場附近,以及漳州天寶試驗點,常規(guī)管理果苗,并在出現(xiàn)冷空氣過程時,按照《芒果寒害記錄表》逐日記錄果樹葉片和葉柄變色情況、枯葉和落葉情況、枯枝和主干受凍情況等寒害癥狀并拍照,共獲得有效災(zāi)情樣本9例。

    1.3.4 人工氣候箱低溫致災(zāi)試驗

    近年來,采用大型步入式人工氣候箱開展人工控制試驗,確定果樹寒害等級指標(biāo)的方法已取得一些研究成果[16?17]。為獲取足夠多的災(zāi)情樣本,在日本產(chǎn)ESPEC大型人工氣候箱內(nèi)開展歷時1d的芒果寒害致災(zāi)試驗。試驗在2015/2016年冬季低于20℃的天氣進行,模擬日較差8℃降溫,溫度設(shè)置范圍為?3.0~5.0℃,設(shè)定6個溫度處理為5.0℃、3.0℃、1.0℃、?1.0℃、?2.0℃、?3.0℃,每個溫度處理為3盆果苗樣本。首先將18盆果苗全部放入人工氣候箱,模擬自然降溫至設(shè)定起始溫度5.0℃,每個溫度處理保持1h恒溫,然后取出3盆果苗置于室內(nèi)并做好標(biāo)記,降溫過程保持每小時降溫1℃的速度,直至溫度降至?3.0℃,持續(xù)1h后取出最后1個處理。試驗結(jié)束后,將試驗所有果苗置于室外,1周后進行觀測,記錄果苗形態(tài)變化。共獲得有效災(zāi)情樣本6例。

    1.4 計算方法

    1.4.1 相對氣象產(chǎn)量的計算

    相對氣象產(chǎn)量比氣象產(chǎn)量能更好地描述以氣象要素為主的各種短期變動因子對產(chǎn)量序列的影響,在對作物或果樹進行產(chǎn)量分析時,可將產(chǎn)量分解為趨勢產(chǎn)量(yt)、氣象產(chǎn)量(yw)和隨機產(chǎn)量(Δy),隨機產(chǎn)量一般可忽略不計[18?19],相對氣象產(chǎn)量可表述為氣象產(chǎn)量與趨勢產(chǎn)量的比值,即

    式中,趨勢產(chǎn)量(yt)采用直線滑動平均進行估算,滑動步長為5a。該方法不需要預(yù)先假定趨勢產(chǎn)量模擬的曲線類型,不會造成資料樣本的損失[20]。

    1.4.2 過程有害積寒值的計算

    研究發(fā)現(xiàn)采用低溫持續(xù)日數(shù)和過程有害積寒可較好地表示中弱冷空氣多次補充造成的平流型寒害的累積作用,因此,采用杜堯東等提出的過程有害積寒[12?13,21],作為芒果寒害過程中重要的氣象因子進行分析。

    式中,X過程為寒害過程的有害積寒值(℃·d);N為寒害過程持續(xù)日數(shù)(d);Tc為寒害臨界溫度(℃),Tmin為日極端最低氣溫(℃);Tm為日平均氣溫(℃)。

    2 結(jié)果與分析

    2.1 芒果寒害臨界溫度的確定

    通過冬季寒害災(zāi)情調(diào)查方法,獲取歷史災(zāi)情樣本;結(jié)合2014/2015年、2015/2016年冬季果園定位觀測、地理移置和人工氣候箱低溫致災(zāi)試驗結(jié)果,整理分析試驗觀測數(shù)據(jù)、拍照記錄、觀測人員記錄資料及相關(guān)氣象數(shù)據(jù),得到有明確時間、地點、果樹受害癥狀、過程極端最低溫度等有效芒果寒害樣本共39例,其中33例為自然條件下調(diào)查和試驗觀測樣本(表1),6例為人工氣候箱低溫致災(zāi)試驗樣本(表2)。

    由表1可知,當(dāng)極端最低溫度低于2.0℃時,芒果出現(xiàn)寒害癥狀,隨著溫度進一步降低,寒害趨于嚴(yán)重甚至芒果被凍死,而當(dāng)溫度高于2.0℃時,芒果無受害癥狀;表2中當(dāng)試驗處理溫度為3.0℃持續(xù)時間1h時,芒果無受害癥狀,但當(dāng)溫度降至1.0℃持續(xù)時間1h時,芒果葉片變色,遭受寒害。由此可以確定芒果的寒害臨界溫度為2.0℃,即當(dāng)極端最低氣溫≤2. 0℃時,芒果寒害過程開始,當(dāng)極端最低氣溫>2. 0℃時,寒害過程結(jié)束。

    表1 芒果寒害災(zāi)情調(diào)查和觀測樣本

    Table 1 Disaster samples of mango cold injury by investigation and observation

    (續(xù)表)

    注:A表示樣本來源于果園定位觀測,B表示樣本來源于地理移置試驗,C表示樣本來源于災(zāi)情調(diào)查。

    Note:“A” indicates that sample was from the orchard positioning observation; “B” indicates that sample was from the geographical displacement experiment;“C”indicates that sample was from the disaster investigation.

    表2 人工氣候箱低溫致災(zāi)試驗芒果苗寒害癥狀

    Table 2 Cold injury symptoms of mango seedlings in low temperature disaster experiment of artificial climate chamber

    2.2 寒害過程篩選結(jié)果

    選取福建省長泰縣、南安市、安溪縣、仙游縣、福清市、寧德市區(qū)6個芒果主產(chǎn)縣市,根據(jù)已確定的芒果寒害臨界溫度2.0℃,篩選6個縣市2001?2014年冬季極端最低溫度≤2. 0℃的寒害過程,確定寒害過程開始和結(jié)束的時間,并從中選取當(dāng)年≤2. 0℃持續(xù)時間最長的寒害過程。由此可得,長泰縣符合條件的寒害過程有4個,南安市符合條件的有1個,安溪縣符合條件的有3個,仙游縣符合條件的有8個,福清市符合條件的有2個,寧德符合條件的有10個,共篩選出28個寒害過程。

    2.3 芒果寒害風(fēng)險綜合氣候指標(biāo)的確定

    2.3.1 綜合氣候指標(biāo)計算

    根據(jù)篩選出的28個芒果寒害過程開展寒害致災(zāi)因子的研究,將過程極端最低溫度、過程平均溫度、過程有害積寒、過程持續(xù)時間、過程降溫幅度、過程降水量與當(dāng)年相對氣象產(chǎn)量進行相關(guān)分析,研究發(fā)現(xiàn)過程極端最低氣溫(X1)、日最低氣溫≤2.0℃有害積寒(X2)、≤2.0℃寒害過程持續(xù)日數(shù)(X3)、≤2.0℃寒害過程降溫幅度(X4)與相對氣象產(chǎn)量的相關(guān)性最好,4個氣象因子均通過了0.05水平的顯著性檢驗,故選取這4個因子作為芒果寒害的致災(zāi)因子。

    芒果寒害4個致災(zāi)因子X1、X2、X3、X4序列間的相關(guān)系數(shù)見表3,由表可以看出,致災(zāi)因子兩兩之間的相關(guān)系數(shù)均通過了0.01水平的顯著性檢驗,表明致災(zāi)因子之間并不獨立,而是互有影響。因此,為了避免信息重疊而影響分析效果,利用主成分分析對4個致災(zāi)因子進行綜合簡化。

    表3 致災(zāi)因子間的相關(guān)系數(shù)矩陣

    注:X1為寒害過程的極端最低氣溫(℃),X2為日最低氣溫≤2.0℃有害積寒、X3為≤2.0℃寒害過程持續(xù)日數(shù)(d)、X4為≤2.0℃寒害過程降溫幅度(℃)

    Note:X1is the extreme minimum temperature of cold injury process, X2is the accumulated value of cold injury process for the extreme minimum daily temperature ≤2.0℃, X3is the sustained days of cold injury process for the temperature ≤2.0℃, X4is the cooling range of cold injury process for the temperature ≤2.0℃.

    運用SPSS統(tǒng)計軟件中的主成分分析工具,對芒果寒害4個致災(zāi)因子X1、X2、X3、X4序列進行分析,軟件提取了第一主成分,其累計方差貢獻率達到76.11%,說明這一主成分已能充分說明數(shù)據(jù)間的波動原因,并得到協(xié)方差矩陣的特征根λ=(2.030,1.014,0.331,0.088),以及特征向量A1、A2、A3、A4分別為?0.656、0.737、0.866和0.492。

    由此可構(gòu)建芒果寒害綜合氣候指標(biāo)(HI)為

    HI=?0.656X1+0.737X2+0.866X3+0.492X4(3)

    由式(3)可知,寒害綜合氣候指標(biāo)HI與冬季極端最低氣溫X1呈負相關(guān)關(guān)系,而與其它3個致災(zāi)因子呈正相關(guān)關(guān)系,即≤2.0℃有害積寒越大,≤2.0℃低溫持續(xù)日數(shù)越多,≤2.0℃寒害過程降溫幅度越大,冬季極端最低氣溫越低,則寒害綜合氣候指標(biāo)HI越大,芒果遭受的寒害等級就越高,這與果樹實際發(fā)生的災(zāi)害過程是相符的。

    2.3.2 綜合氣候指標(biāo)分級

    參照農(nóng)業(yè)上劃分災(zāi)害年型的方法,將相對氣象產(chǎn)量減產(chǎn)率(%)分為0~10%、10.1%~20%、20.1%~30%、30.1%~100%共 4個等級,分別對應(yīng)輕、中、重、嚴(yán)重寒害等級[22?23]。將計算篩選出的芒果主產(chǎn)縣市28個寒害過程綜合氣候指標(biāo)(HI)與當(dāng)年芒果相對氣象產(chǎn)量繪制點聚圖(圖1)。

    圖1 28個寒害過程綜合氣候指標(biāo)[HI,由公式(3)計算]與相對氣象產(chǎn)量間關(guān)系

    由圖1可見,芒果寒害綜合氣候指標(biāo)HI與相對氣象產(chǎn)量呈顯著負相關(guān)關(guān)系(R=?0.383,P<0.05)。采用線性回歸方法對兩個因子進行分析,可得到一元線性回歸方程。

    y=?0.142x+10.12 (4)

    式中,y為芒果寒害綜合氣候指標(biāo),x為相對氣象產(chǎn)量。用相對氣象減產(chǎn)率(%)x=?30、?20、?10、0代入方程,得到芒果各級寒害綜合氣候指標(biāo)的閾值,并據(jù)此確定各寒害等級的指標(biāo)區(qū)間(表4)。

    2.3.3 綜合氣候指標(biāo)驗證

    選取寧德市2001?2014年冬季芒果寒害綜合氣候指標(biāo)與相對氣象產(chǎn)量進行對比印證(圖2)。從圖2可知,芒果寒害綜合氣候指標(biāo)與相對氣象產(chǎn)量呈現(xiàn)明顯相反的變化趨勢,芒果產(chǎn)量較低的年份,寒害綜合氣候指標(biāo)數(shù)值就大,反之產(chǎn)量高的年份,則綜合氣候指標(biāo)就小。其中2005、2009、2011年芒果產(chǎn)量較低,相應(yīng)的寒害綜合氣候指標(biāo)較大,這與寧德市2004?2005年冬季、2008?2009年冬季和2010?2011年冬季芒果遭受低溫寒害,從而導(dǎo)致當(dāng)年產(chǎn)量降低的實際情況相一致。計算芒果寒害綜合氣候指標(biāo)與相對氣象產(chǎn)量的相關(guān)系數(shù)為?0.606(P<0.05)。綜上所述,本研究的寒害綜合氣候指標(biāo)可以用來分析福建省芒果寒害的輕重。

    表4 芒果寒害綜合氣候指標(biāo)等級

    Table 4 Grade of comprehensive climatic index for mango cold injury

    2.4 福建芒果寒害種植氣候風(fēng)險區(qū)劃

    建立芒果寒害綜合氣候指標(biāo)與測站經(jīng)度、緯度和海拔高度的多元回歸模型,利用GIS技術(shù)對其進行地理插值推算,將點狀數(shù)據(jù)轉(zhuǎn)化為面狀格網(wǎng)數(shù)據(jù),解決無氣象觀測資料地區(qū)的寒害指數(shù)計算問題,實現(xiàn)寒害綜合氣候指標(biāo)的空間連續(xù)分布。

    由福建芒果寒害種植氣候風(fēng)險區(qū)劃圖可見(圖3),福建省芒果無寒害?輕度寒害的區(qū)域主要分布在沿海地區(qū),是福建省引種芒果災(zāi)害風(fēng)險最低的區(qū)域,適宜芒果的引種擴種,但遇冬季冷空氣過程時也應(yīng)采取灌水、稻草覆蓋、熏煙、高培土、樹干包扎、果實套袋、增施熱性磷鉀肥、噴施果樹防凍劑等防寒措施,以減輕或避免寒害損失;中度寒害區(qū)主要分布在南靖縣、上杭縣、福鼎市、三明市區(qū)、德化縣,是福建省引種芒果災(zāi)害風(fēng)險較低的區(qū)域,由于芒果樹體較大,采用設(shè)施大棚進行種植成本較高,建議這些區(qū)域不作引種;其余地區(qū)為重度寒害?嚴(yán)重寒害區(qū)域,不適宜芒果的引種試種。

    圖3 福建省芒果寒害種植氣候風(fēng)險區(qū)劃

    3 結(jié)論與討論

    3.1 討論

    通過冬季寒害災(zāi)情調(diào)查、果園定位觀測、地理移置和人工氣候箱低溫致災(zāi)試驗4種渠道收集芒果寒害災(zāi)情樣本39例,其中33例為自然條件下調(diào)查和試驗觀測樣本,6例為人工氣候箱低溫致災(zāi)試驗樣本,從災(zāi)情樣本中可以看到當(dāng)極端最低溫度低于2.0℃時,芒果出現(xiàn)寒害癥狀,從而確定芒果寒害臨界溫度為2.0℃,這與陳惠等[24?25]研究的芒果寒害起始受害極端最低氣溫2.0℃相一致,即當(dāng)最低氣溫≤2.0℃時,芒果開始遭受寒害。

    根據(jù)確定的芒果寒害臨界溫度,從福建芒果主產(chǎn)縣市2001?2014年冬季低溫過程中共篩選出28個芒果寒害過程。通過將寒害過程溫度、降水等氣象因子與相對氣象產(chǎn)量進行相關(guān)分析,確定過程極端最低氣溫、過程有害積寒、過程持續(xù)日數(shù)和過程最大降溫幅度為芒果寒害的4個致災(zāi)因子,因子之間相關(guān)性顯著?;?個致災(zāi)因子,采用主成分分析方法構(gòu)建芒果寒害綜合氣候指標(biāo),其第一主成分累計方差貢獻率達到76.11%,物理意義清晰。根據(jù)災(zāi)害年型的劃分方法,最終確定輕、中、重、嚴(yán)重寒害綜合氣候指標(biāo)等級分別為10.12<HI≤11.54、11.54<HI≤12.96、12.96<HI≤14.38、HI>14.38。為了驗證指標(biāo)的代表性和準(zhǔn)確性,對寧德市歷年相對氣象產(chǎn)量和寒害綜合氣候指標(biāo)進行試算,結(jié)果表明兩者有顯著的負相關(guān)關(guān)系,指標(biāo)很好地反映了芒果災(zāi)害年遭受的寒害輕重,說明該指標(biāo)對評價芒果寒害受害程度有實際參考價值。從芒果寒害綜合氣候區(qū)劃結(jié)果來看,無寒害?輕度寒害的區(qū)域主要分布在福建省沿海地區(qū),是引種芒果優(yōu)質(zhì)品種寒害風(fēng)險最低的區(qū)域,適宜引種擴種,這也與芒果的實際種植區(qū)域基本一致,區(qū)劃結(jié)果可為芒果寒害監(jiān)測預(yù)警、農(nóng)業(yè)生產(chǎn)區(qū)劃及產(chǎn)業(yè)帶建設(shè)提供理論依據(jù)。

    應(yīng)當(dāng)指出的是,影響果樹寒害的因子很復(fù)雜,其發(fā)生及程度還受植物學(xué)因子的影響[26?28]。本研究未對芒果遭受寒害后植物生理指標(biāo)進行測定,僅針對氣象因子進行了分析,如何結(jié)合陰雨、大風(fēng)、日照等設(shè)計寒害指標(biāo)試驗,并得到寒害預(yù)警業(yè)務(wù)所需的更精細、更準(zhǔn)確綜合性指標(biāo),還有待今后進一步研究。

    3.2 結(jié)論

    通過對典型寒害過程分析,芒果寒害臨界溫度為2.0℃;利用芒果寒害綜合氣候指標(biāo),確定輕、中、重、嚴(yán)重寒害綜合氣候指標(biāo)等級分別為10.12<HI≤11.54、11.54<HI≤12.96、12.96<HI≤14.38、HI>14.38;從綜合氣候區(qū)劃結(jié)果看,福建省沿海地區(qū)為無寒害或輕度寒害區(qū)域,是引種優(yōu)質(zhì)品種寒害風(fēng)險最低的區(qū)域。

    [1]陳華蕊,陳業(yè)淵,高愛平,等.芒果光合特性的研究[J].熱帶作物學(xué)報,2010,31(5):729-735.

    Chen H R,Chen Y Y,Gao A P,et al.Photosynthetic characteristics of mango[J].Chinese Journal of Tropical Crops,2010,31(5):729-735.(in Chinese)

    [2]任傳帥,葉回春,崔貝,等.基于面向?qū)ο蠓诸惖拿⒐诌b感提取方法研究[J].資源科學(xué),2017,39(8):1584-1591.

    Ren C S,Ye H C,Cui B,et al.Acreage estimation of mango orchards using object-oriented classification and remote sensing[J].Resources Science,2017,39(8):1584-1591.(in Chinese)

    [3]徐宗煥,陳家金,張容焱,等.引種臺灣愛文芒果的低溫災(zāi)害評估及GIS應(yīng)用[J].自然災(zāi)害學(xué)報,2006,15(6):260-264.

    Xu Z H,Chen J J,Zhang R Y,et al.Low temperature disaster evaluation and GIS application to adventitious Planting of Aiwen mango from Taiwan[J].Journal of Natural Disasters, 2006,15(6):260-264.(in Chinese)

    [4]余會康,陳安芳,阮翠冰.閩東氣候變化與杧果生長分析[J].中國農(nóng)業(yè)資源與區(qū)劃,2014,35(4):63-68.

    Yu H K,Chen A F,Ruan C B.Analysis of climate change and the growth of mango in northeastern Fujian province[J]. Chinese Journal of Agricultural Resources and Regional Planning,2014,35(4):63-68.(in Chinese)

    [5]魏秀清,余東,章希娟,等. 6個杧果品種耐寒性的研究[J].果樹學(xué)報,2017,34(5):594-602.

    Wei X Q,Yu D,Zhang X J,et al.A study on cold resistance in six mango cultivars[J].Journal of Fruit Science,2017,34(5): 594-602.(in Chinese)

    [6]劉錦鑾,杜堯東,毛慧琴.華南地區(qū)荔枝寒害風(fēng)險分析與區(qū)劃[J].自然災(zāi)害學(xué)報,2003,12(3):126-130.

    Liu J L,Du Y D,Mao H Q.Risk analysis and zonation of cold damage to litchi in South China[J].Journal of Natural Disasters,2003,12(3):126-130.(in Chinese)

    [7]鄒海平,王春乙,張京紅等.海南島香蕉寒害風(fēng)險區(qū)劃[J].自然災(zāi)害學(xué)報,2013,22(3):130-134.

    Zou H P,Wang C Y,Zhang J H,et al.Risk zoning of cold damage to banana in Hainan Island[J].Journal of Natural Disasters,2013,22(3):130-134.(in Chinese)

    [8]蔡文華,陳惠,潘衛(wèi)華,等.福建龍眼樹凍害指標(biāo)初探[J].中國農(nóng)業(yè)氣象,2009,30(1):109-112.

    Cai W H,Chen H,Pan W H,et al.A preliminary study on freezing injury indicators of longan trees[J].Chinese Journal of Agrometeorology,2009,30(1):109-112.(in Chinese)

    [9]陳惠,王加義,潘衛(wèi)華,等.南亞熱帶主要果樹凍(寒)害低溫指標(biāo)的確定[J].中國農(nóng)業(yè)氣象,2012,33(1):148-155.

    Chen H,Wang J Y,Pan W H,et al.Determination of frozen(cold) injury indicators of main fruit trees in South Subtropical area[J].Chinese Journal of Agrometeorology,2012,33(1): 148-155.(in Chinese)

    [10]陳惠,楊凱,林晶,等.枇杷幼果凍害低溫等級指標(biāo)的確定[J].中國農(nóng)業(yè)氣象,2016,37(1):91-97.

    Chen H,Yang K,Lin J,et al.Determination of frozen injury grade index for young fruit of loquat[J].Chinese Journal of Agrometeorology,2016,37(1):91-97.(in Chinese)

    [11]陳惠,楊凱,李麗純,等.蓮霧寒凍害低溫等級指標(biāo)的確定[J].生態(tài)學(xué)雜志,2017,36(4):1047-1052.

    Chen H,Yang K,Li L C,et al.Determination of cold and freezing injury grade index of wax apple[J].Chinese Journal of Ecology,2017,36(4):1047-1052.(in Chinese)

    [12]杜堯東,李春梅,毛慧琴.廣東省香蕉與荔枝寒害致災(zāi)因子和綜合氣候指標(biāo)研究[J].生態(tài)學(xué)雜志,2006,25(2):225-230.

    Du Y D,Li C M,Mao H Q.Disaster-inducing factors and integrated climatic index for banana and litchi chilling injuries in Guangdong province[J].Chinese Journal of Ecology,2006,25(2):225-230.(in Chinese)

    [13]李娜,霍治國,賀楠,等.華南地區(qū)香蕉、荔枝寒害的氣候風(fēng)險區(qū)劃[J].應(yīng)用生態(tài)學(xué)報,2010,21(5):1244-1251.

    Li N,Huo Z G,He N,et al.Climatic risk zoning for banana and litchi chilling injury in South China[J].Chinese Journal of Applied Ecology,2010,21(5):1244-1251.(in Chinese)

    [14]楊凱,林晶,陳惠,等.福建枇杷低溫害臨界溫度和綜合氣候指標(biāo)[J].中國農(nóng)業(yè)氣象,2013,34(3):468-473.

    Yang K,Lin J,Chen H,et al.Critical temperature and integrated climatic index of low temperature injury for loquat in Fujian province[J].Chinese Journal of Agrometeorology, 2013,34(3):468-473.(in Chinese)

    [15]霍治國,王石立.農(nóng)業(yè)和生物氣象災(zāi)害[M].北京:氣象出版社,2009:51-64.

    Huo Z G,Wang S L.Agricultural and biometeorological disasters[M].Beijing:China Meteorological Press,2009:51-64. (in Chinese)

    [16]谷曉平,袁小康,胡家敏.火龍果幼苗和成齡書寒害指標(biāo)初探[J].中國農(nóng)業(yè)氣象,2014,35(2):214-220.

    Gu X P,Yuan X K,Hu J M.Indices of chilling injury for pitaya seedling and mature plant[J].Chinese Journal of Agrometeorology,2014,35(2):214-220.(in Chinese)

    [17]楊凱,陳惠,李麗純,等.引種臺灣青棗的寒凍害等級指標(biāo)研究[J].自然災(zāi)害學(xué)報,2017,26(4):91-97.

    Yang K,Chen H,Li L C,et al.Study on cold and freezing injury grade index in introduction ofLam[J].Journal of Natural Disasters,2017,26(4):91-97.(in Chinese)

    [18]姜麗霞,王育光,孫孟梅,等.黑龍江省玉米產(chǎn)量預(yù)報模式的研究[J].中國農(nóng)業(yè)氣象,2004,25(1):13-16.

    Jiang L X,Wang Y G,Sun M M,et al.A Study on forecast model for maize climate yield in Heilongjiang province[J]. Chinese Journal of Agrometeorology,2004,25(1):13-16.(in Chinese)

    [19]王馥棠,李郁竹,王石立.農(nóng)業(yè)產(chǎn)量氣象模擬與模型引論[M].北京:科學(xué)出版社,1990:40-41.

    Wang F T, Li Y Z,Wang S L.Introduction of meteorological simulation and model to agricultural yield[M].Beijing: Science Press,1990:40-41.(in Chinese)

    [20]孫衛(wèi)國,程炳巖,楊沈斌,等.區(qū)域氣候變化對華東地區(qū)水稻產(chǎn)量的影響[J].中國農(nóng)業(yè)氣象,2011,32(2):227-234.

    Sun W G,Cheng B Y,Yang S B,et al.Impact of the regional climate variations on rice yield in East China[J].Chinese Journal of Agrometeorology,2011,32(2):227-234.(in Chinese)

    [21]馮穎竹,梁紅,黃璜.廣東冬季寒害指標(biāo)研究[J].自然災(zāi)害學(xué)報,2005,14(1):59-65.

    Feng Y Z,Liang H,Huang H.Study on index for winter’s chilling damage in Guangdong Province[J].Journal of Natural Disasters,2005,14(1):59-65.(in Chinese)

    [22]王明田,張玉芳,馬均,等.四川盆地區(qū)玉米干旱災(zāi)害風(fēng)險評價及區(qū)劃[J].應(yīng)用生態(tài)學(xué)報,2012,23(10):2803-2811.

    Wang M T,Zhang Y F,Ma J,et al.Risk assessment and regionalization of maize drought disasters in Sichuan Basin,Southeast China[J].Chinese Journal of Applied Ecology,2012,23(10):2803-2811.(in Chinese)

    [23]薛昌穎,劉榮花,馬志紅.黃淮海地區(qū)夏玉米干旱等級劃分[J].農(nóng)業(yè)工程學(xué)報,2014,30(16):147-156.

    Xue C Y,Liu R H,Ma Z H.Drought grade classification of summer maize in Huang-Huai-Hai area[J].Transactions of the CSAE,2014,30(16):147-156.(in Chinese)

    [24]陳惠,楊凱,李政,等.3種熱帶特色果樹寒凍害低溫等級指標(biāo)的確定[J].果樹學(xué)報,2018,35(1):82-93.

    Chen H,Yang K,Li Z,et al.Determination of cold and freezing injury grade indexes for three tropical characteristic fruit trees[J].Journal of Fruit Science,2018,35(1):82-93.(in Chinese)

    [25]李政,蘇永秀,王瑩,等.芒果寒(凍)害等級劃分及低溫指標(biāo)確定[J].災(zāi)害學(xué),2017,32(3):18-22.

    Li Z,Su Y X,Wang Y,et al.Determination of cold (freezing) injury classification and indicator of low temperature on mango[J].Journal of Catastrophology,2017,32(3):18-22.(in Chinese)

    [26]馮美利.杧果葉片組織的細胞結(jié)構(gòu)與耐寒性的初步研究[J].熱帶農(nóng)業(yè)科技,2003,26(1):5-7.

    Feng M L.Relations between the cellularity of mango leaf and resistance to coldness[J].Tropical Agricultural Science and Technology,2003,26(1):5-7.(in Chinese)

    [27]唐力生,王華,胡飛,等.低溫脅迫下芒果苗的受害癥狀及生理響應(yīng)[J].生態(tài)學(xué)雜志,2016,35(10):2627-2636.

    Tang L S,Wang H,Hu F,et al.Cold injury symptoms and physiological responses of mango seedlings under low temperature stress[J].Chinese Journal of Ecology,2016, 35(10):2627-2636.(in Chinese)

    [28]胡春霞,王麗,湯杰.低溫對南果梨的生理生化指標(biāo)的影響[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報,2009,40(3):349-352.

    Hu C X,Wang L,Tang J.Effect of low temperature on physiological indicators of Nanguo pear[J].Journal of Shengyang Agricultural University,2009,40(3):349-352.(in Chinese)

    Climatic Risk Regionalization of Mango Planting in Fujian Province Based on Cold Injury Process

    YANG Kai1,CHEN Bin-bin1,CHEN Hui1,LI Li-chun1,CHEN Bin-yuan2,LI Zheng3

    (1. Fujian Institute of Meteorological Science, Fuzhou 350001, China; 2. Zhangzhou Meteorological Bureau, Zhangzhou 363000; 3. Institute of Guangxi Meteorological Disaster Mitigation, Nanning 530022)

    Based on the39 mango cold disaster samples collected through four channels which were disaster investigation of cold injury in winter, orchard positioning observation, geographic displacement experiment and low temperature disaster experiment of artificial climate chamber, the critical temperature of cold injury for mango was determined to be 2.0℃. On this basis, 28 cold injury processes were screened out in winter from 2001 to 2014 in Fujian major mango producing counties and cities.By correlating the meteorological factors such as temperature and precipitation in the process of cold injury with the relative meteorological yield, the disaster-inducing factors for cold injury of mango was determined as extreme minimum temperature of cold injury process, accumulated value of cold injury process for the extreme minimum daily temperature ≤2.0℃, sustained days of cold injury process for the temperature ≤2.0℃, and cooling range of cold injury process for the temperature ≤2.0℃. By using the method of principal component analysis to simplify comprehensively four disaster-inducing factors, the integrated climatic index of cold injury evaluation of mango was obtained, and the index classification was determined by combining the relative meteorological yield.According to the grade of comprehensive climatic index of cold injury, the climatic risk regionalization of mango planting for cold injury was completed by using GIS method in Fujian Province. The region with no or mild cold injury mainly distributes in the coastal areas, which is the region with the lowest risk of cold injury for introducing high-quality mango varieties.

    Mango; Cold injury; Critical temperature; Comprehensive climatic index; Risk regionalization

    10.3969/j.issn.1000-6362.2019.11.006

    楊凱,陳彬彬,陳惠,等.基于寒害過程的福建芒果種植氣候風(fēng)險區(qū)劃[J].中國農(nóng)業(yè)氣象,2019,40(11):723-732

    2019?04?26

    。E-mail:ch775@126.com

    公益性行業(yè)(氣象)科研專項(GYHY201406027);福建省自然基金項目(2018J01047);福建省科技廳農(nóng)業(yè)引導(dǎo)性(重點)項目(2017N0011);福建省氣象局開放式氣象科學(xué)研究基金項目(2015K02)作者簡介:楊凱(1979?),碩士,高級工程師,從事農(nóng)業(yè)氣象和氣候變化研究。E-mail:yangkai_1998@163.com

    猜你喜歡
    區(qū)劃芒果氣候
    我才不要穿
    林芝市雷電風(fēng)險區(qū)劃
    基于“三線一單”的環(huán)境功能區(qū)劃實證研究
    川渝傳統(tǒng)民居區(qū)劃的環(huán)境蘊意
    小洞會“咬”人
    小洞會“咬”人
    瞧,氣候大不同
    氣候變暖會怎樣?
    立冬
    氣候的未來掌握在我們手中
    变态另类丝袜制服| 亚洲成人久久爱视频| 亚洲av电影不卡..在线观看| 小蜜桃在线观看免费完整版高清| 三级国产精品欧美在线观看| 老汉色av国产亚洲站长工具| 久9热在线精品视频| 欧美极品一区二区三区四区| 亚洲av五月六月丁香网| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 嫁个100分男人电影在线观看| 亚洲av电影不卡..在线观看| 波野结衣二区三区在线 | 国产真实乱freesex| 级片在线观看| 在线观看66精品国产| 久久久久久九九精品二区国产| 国产伦精品一区二区三区视频9 | 免费人成在线观看视频色| 一区二区三区激情视频| 亚洲18禁久久av| 日本 欧美在线| 亚洲av不卡在线观看| 熟妇人妻久久中文字幕3abv| 日韩 欧美 亚洲 中文字幕| 18禁在线播放成人免费| 日本撒尿小便嘘嘘汇集6| 国产激情偷乱视频一区二区| 午夜福利在线在线| 黄色视频,在线免费观看| 少妇人妻精品综合一区二区 | 伊人久久精品亚洲午夜| 国产不卡一卡二| 搡老妇女老女人老熟妇| 久久草成人影院| 好看av亚洲va欧美ⅴa在| 亚洲精品亚洲一区二区| 99国产综合亚洲精品| 国产成人av教育| 国产激情偷乱视频一区二区| 日本与韩国留学比较| 噜噜噜噜噜久久久久久91| 狂野欧美白嫩少妇大欣赏| 欧美最黄视频在线播放免费| 欧美又色又爽又黄视频| 观看美女的网站| 亚洲成a人片在线一区二区| 中文字幕av在线有码专区| 久久久久久九九精品二区国产| 国产黄a三级三级三级人| 国产欧美日韩一区二区精品| 国产激情偷乱视频一区二区| www.999成人在线观看| 99久久无色码亚洲精品果冻| 成熟少妇高潮喷水视频| 免费av不卡在线播放| 草草在线视频免费看| 色播亚洲综合网| 亚洲男人的天堂狠狠| 岛国在线观看网站| 欧美午夜高清在线| 高清毛片免费观看视频网站| 久久久久久久亚洲中文字幕 | 五月玫瑰六月丁香| 在线免费观看不下载黄p国产 | 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| АⅤ资源中文在线天堂| 可以在线观看的亚洲视频| 天堂av国产一区二区熟女人妻| 日本免费一区二区三区高清不卡| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 亚洲欧美日韩东京热| 国产野战对白在线观看| 精品免费久久久久久久清纯| svipshipincom国产片| 婷婷精品国产亚洲av| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 国内精品久久久久精免费| 999久久久精品免费观看国产| 99精品欧美一区二区三区四区| 精品欧美国产一区二区三| 人妻丰满熟妇av一区二区三区| 日本黄大片高清| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 欧美乱码精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 老汉色∧v一级毛片| 国产一区二区在线av高清观看| 成人性生交大片免费视频hd| 又黄又粗又硬又大视频| 中文字幕人成人乱码亚洲影| 午夜福利视频1000在线观看| 欧美日本视频| 国产主播在线观看一区二区| 国产欧美日韩一区二区精品| 国产精品亚洲美女久久久| 91av网一区二区| 久久6这里有精品| 亚洲国产精品久久男人天堂| ponron亚洲| 亚洲性夜色夜夜综合| 狂野欧美激情性xxxx| 激情在线观看视频在线高清| 深夜精品福利| 成人av在线播放网站| 人人妻人人看人人澡| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 人人妻,人人澡人人爽秒播| 欧美一区二区亚洲| 高清在线国产一区| 成人特级黄色片久久久久久久| 五月玫瑰六月丁香| 成年免费大片在线观看| 欧美精品啪啪一区二区三区| 亚洲18禁久久av| 午夜激情欧美在线| 国产成+人综合+亚洲专区| 午夜福利欧美成人| 男人的好看免费观看在线视频| 免费在线观看日本一区| 亚洲精品在线观看二区| 日本 欧美在线| 国产成人av激情在线播放| 欧美一区二区精品小视频在线| 乱人视频在线观看| 日韩欧美国产在线观看| 久久久久久九九精品二区国产| 99国产综合亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 麻豆成人av在线观看| 亚洲精品美女久久久久99蜜臀| 午夜激情欧美在线| 国产老妇女一区| 十八禁网站免费在线| 中文字幕精品亚洲无线码一区| 搞女人的毛片| 久久精品国产亚洲av香蕉五月| 成年人黄色毛片网站| 看片在线看免费视频| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频| 天堂av国产一区二区熟女人妻| 欧美av亚洲av综合av国产av| 午夜福利视频1000在线观看| 精品久久久久久久人妻蜜臀av| 色在线成人网| 男人的好看免费观看在线视频| 一级毛片高清免费大全| 村上凉子中文字幕在线| 久久九九热精品免费| 9191精品国产免费久久| 一进一出好大好爽视频| 国产日本99.免费观看| 日韩精品中文字幕看吧| 19禁男女啪啪无遮挡网站| 成年免费大片在线观看| 又紧又爽又黄一区二区| 免费大片18禁| 亚洲在线自拍视频| 少妇丰满av| 国产亚洲av嫩草精品影院| 中文字幕高清在线视频| 国产欧美日韩一区二区三| 午夜精品久久久久久毛片777| 又黄又粗又硬又大视频| 国产黄a三级三级三级人| 老司机午夜福利在线观看视频| 国产精品亚洲一级av第二区| 99国产综合亚洲精品| 国内久久婷婷六月综合欲色啪| 人妻夜夜爽99麻豆av| 日本免费a在线| 国产av一区在线观看免费| 久久久久久人人人人人| 一边摸一边抽搐一进一小说| 国产亚洲精品一区二区www| 黄色女人牲交| 国产美女午夜福利| 国产亚洲精品久久久久久毛片| 成人特级av手机在线观看| 特大巨黑吊av在线直播| 午夜福利在线观看免费完整高清在 | 日韩国内少妇激情av| 免费电影在线观看免费观看| 变态另类成人亚洲欧美熟女| 少妇的逼好多水| 欧美中文日本在线观看视频| 国产老妇女一区| 亚洲av电影在线进入| 日本在线视频免费播放| 午夜影院日韩av| 成人无遮挡网站| 日日摸夜夜添夜夜添小说| 丁香欧美五月| 国产成人aa在线观看| 手机成人av网站| 一二三四社区在线视频社区8| 欧美中文日本在线观看视频| 日韩免费av在线播放| 亚洲黑人精品在线| 好看av亚洲va欧美ⅴa在| 美女黄网站色视频| 成人鲁丝片一二三区免费| 国产亚洲精品久久久久久毛片| 舔av片在线| 午夜久久久久精精品| 国产av在哪里看| 日韩欧美精品免费久久 | av天堂在线播放| 欧美成人一区二区免费高清观看| 又黄又粗又硬又大视频| 十八禁网站免费在线| 色av中文字幕| 亚洲美女黄片视频| 别揉我奶头~嗯~啊~动态视频| 熟妇人妻久久中文字幕3abv| 九九热线精品视视频播放| 尤物成人国产欧美一区二区三区| 亚洲精品日韩av片在线观看 | 成人高潮视频无遮挡免费网站| 老汉色av国产亚洲站长工具| 国产又黄又爽又无遮挡在线| 日本 av在线| 少妇人妻精品综合一区二区 | av中文乱码字幕在线| 午夜老司机福利剧场| 国产伦一二天堂av在线观看| 久久精品综合一区二区三区| 男女那种视频在线观看| 久久这里只有精品中国| 久久草成人影院| 久久久久久人人人人人| 欧美乱妇无乱码| 岛国在线观看网站| 我要搜黄色片| 在线播放无遮挡| 黄片小视频在线播放| 久久久国产成人精品二区| 亚洲中文字幕日韩| 久久亚洲真实| 老司机深夜福利视频在线观看| 国产精品一区二区三区四区久久| 中文字幕人妻熟人妻熟丝袜美 | 观看免费一级毛片| 午夜a级毛片| 免费在线观看亚洲国产| 在线观看av片永久免费下载| 午夜精品在线福利| 麻豆成人av在线观看| 天美传媒精品一区二区| 18+在线观看网站| 91麻豆精品激情在线观看国产| 夜夜夜夜夜久久久久| 国产午夜福利久久久久久| 51午夜福利影视在线观看| 免费人成在线观看视频色| 久久久成人免费电影| 欧美丝袜亚洲另类 | 免费av不卡在线播放| 欧美性感艳星| 91在线观看av| 成年女人永久免费观看视频| 国产探花在线观看一区二区| 日韩欧美一区二区三区在线观看| 校园春色视频在线观看| 亚洲精品日韩av片在线观看 | 国产一区二区三区视频了| 母亲3免费完整高清在线观看| 国产亚洲精品av在线| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 国产精品久久久久久精品电影| 又爽又黄无遮挡网站| ponron亚洲| 禁无遮挡网站| 欧美三级亚洲精品| 亚洲精品久久国产高清桃花| 亚洲色图av天堂| 日韩欧美免费精品| 成人欧美大片| 99在线人妻在线中文字幕| 亚洲天堂国产精品一区在线| 国产蜜桃级精品一区二区三区| 亚洲第一欧美日韩一区二区三区| av中文乱码字幕在线| 成人性生交大片免费视频hd| 国产高清视频在线播放一区| 99riav亚洲国产免费| 激情在线观看视频在线高清| 免费人成视频x8x8入口观看| 欧美不卡视频在线免费观看| 亚洲七黄色美女视频| 久久久久久久久中文| 亚洲熟妇熟女久久| 看免费av毛片| 亚洲激情在线av| 欧美成人免费av一区二区三区| 2021天堂中文幕一二区在线观| 一进一出抽搐动态| 大型黄色视频在线免费观看| 精品熟女少妇八av免费久了| 国产爱豆传媒在线观看| 操出白浆在线播放| 久久久久久人人人人人| 日韩有码中文字幕| 一进一出抽搐gif免费好疼| 男人的好看免费观看在线视频| 国产精品98久久久久久宅男小说| 嫩草影院入口| 国产黄片美女视频| 欧美xxxx黑人xx丫x性爽| 欧美绝顶高潮抽搐喷水| 18禁美女被吸乳视频| 午夜日韩欧美国产| av在线天堂中文字幕| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 日本一本二区三区精品| 天堂动漫精品| 怎么达到女性高潮| 国产成人福利小说| 美女被艹到高潮喷水动态| www.999成人在线观看| 久久人妻av系列| 热99在线观看视频| 成人18禁在线播放| 亚洲激情在线av| 久久久久精品国产欧美久久久| 男人舔女人下体高潮全视频| 亚洲精品乱码久久久v下载方式 | 99久久精品一区二区三区| 长腿黑丝高跟| 99久久精品一区二区三区| 哪里可以看免费的av片| 看片在线看免费视频| 天堂√8在线中文| 天堂网av新在线| 欧美黑人欧美精品刺激| 免费在线观看日本一区| 国产亚洲欧美在线一区二区| 精品久久久久久久人妻蜜臀av| 国产黄a三级三级三级人| 国产探花在线观看一区二区| 国产精品日韩av在线免费观看| 窝窝影院91人妻| 免费高清视频大片| 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| 亚洲av中文字字幕乱码综合| 嫩草影院精品99| 看黄色毛片网站| 久久久久久久午夜电影| 午夜日韩欧美国产| 亚洲狠狠婷婷综合久久图片| 亚洲一区二区三区不卡视频| 亚洲色图av天堂| 久久久久久久精品吃奶| 欧美最新免费一区二区三区 | 成人性生交大片免费视频hd| 日韩欧美国产一区二区入口| 首页视频小说图片口味搜索| 欧美bdsm另类| 少妇的丰满在线观看| 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| www日本黄色视频网| 狂野欧美激情性xxxx| 91av网一区二区| 欧美成人免费av一区二区三区| 国产真人三级小视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲片人在线观看| 天堂影院成人在线观看| 高清日韩中文字幕在线| 中文在线观看免费www的网站| 国产精品,欧美在线| 国产色爽女视频免费观看| 美女高潮喷水抽搐中文字幕| 可以在线观看毛片的网站| 亚洲18禁久久av| 免费观看精品视频网站| 亚洲内射少妇av| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 韩国av一区二区三区四区| 我要搜黄色片| 精品福利观看| 99riav亚洲国产免费| 动漫黄色视频在线观看| 一级作爱视频免费观看| 此物有八面人人有两片| 在线看三级毛片| 久久人人精品亚洲av| 一级黄片播放器| 99国产极品粉嫩在线观看| 色视频www国产| 无遮挡黄片免费观看| 狂野欧美激情性xxxx| avwww免费| 尤物成人国产欧美一区二区三区| 国产高清视频在线观看网站| 午夜视频国产福利| 国产精品女同一区二区软件 | 在线国产一区二区在线| 在线观看美女被高潮喷水网站 | 亚洲 欧美 日韩 在线 免费| 亚洲成人久久爱视频| 午夜免费成人在线视频| 精品午夜福利视频在线观看一区| 久久婷婷人人爽人人干人人爱| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 少妇的逼水好多| 99国产综合亚洲精品| 成人高潮视频无遮挡免费网站| 欧美黑人巨大hd| 国产精品99久久99久久久不卡| 日韩欧美精品免费久久 | 国产黄a三级三级三级人| 午夜日韩欧美国产| 亚洲av熟女| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 日韩人妻高清精品专区| 欧美丝袜亚洲另类 | 午夜影院日韩av| 婷婷丁香在线五月| 国产精品免费一区二区三区在线| 99精品在免费线老司机午夜| 露出奶头的视频| 特级一级黄色大片| 久久久久国内视频| 每晚都被弄得嗷嗷叫到高潮| av片东京热男人的天堂| 禁无遮挡网站| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 男女下面进入的视频免费午夜| 日韩欧美国产在线观看| 国产精品日韩av在线免费观看| 一本综合久久免费| 久久天躁狠狠躁夜夜2o2o| 全区人妻精品视频| 国产成人aa在线观看| 亚洲自拍偷在线| 亚洲人成网站在线播| 欧美黄色片欧美黄色片| 一区二区三区国产精品乱码| 成人午夜高清在线视频| 亚洲电影在线观看av| 黄色片一级片一级黄色片| 国产成人系列免费观看| 观看免费一级毛片| 亚洲狠狠婷婷综合久久图片| 网址你懂的国产日韩在线| 嫁个100分男人电影在线观看| 国内少妇人妻偷人精品xxx网站| 欧美丝袜亚洲另类 | 美女大奶头视频| 久久久久精品国产欧美久久久| АⅤ资源中文在线天堂| 搞女人的毛片| 嫁个100分男人电影在线观看| 91字幕亚洲| 又爽又黄无遮挡网站| 免费av毛片视频| 免费观看精品视频网站| 国产黄色小视频在线观看| 一个人看的www免费观看视频| 午夜精品在线福利| 色在线成人网| 天堂影院成人在线观看| or卡值多少钱| av福利片在线观看| 国产麻豆成人av免费视频| av国产免费在线观看| 午夜久久久久精精品| 欧美一区二区国产精品久久精品| 亚洲最大成人中文| 精品欧美国产一区二区三| 国产三级在线视频| 欧美大码av| 淫秽高清视频在线观看| 婷婷精品国产亚洲av| 在线观看日韩欧美| 97人妻精品一区二区三区麻豆| 日韩欧美国产在线观看| 国产一区二区三区在线臀色熟女| 亚洲精品一卡2卡三卡4卡5卡| 成人av在线播放网站| 亚洲精品日韩av片在线观看 | 久久精品夜夜夜夜夜久久蜜豆| 中文字幕久久专区| 欧美日韩乱码在线| 熟女电影av网| 成年免费大片在线观看| 免费电影在线观看免费观看| 国产熟女xx| 精品久久久久久久末码| 日韩免费av在线播放| 亚洲av不卡在线观看| 波多野结衣高清作品| 国产精品99久久久久久久久| 国产精华一区二区三区| 亚洲乱码一区二区免费版| 中文在线观看免费www的网站| 国产精品爽爽va在线观看网站| 国产黄片美女视频| 三级毛片av免费| 草草在线视频免费看| 天堂av国产一区二区熟女人妻| 99国产精品一区二区蜜桃av| 久久精品91蜜桃| 99久久精品国产亚洲精品| 久久久久久久久大av| 99久久99久久久精品蜜桃| 一区二区三区免费毛片| 国产精品影院久久| 亚洲天堂国产精品一区在线| 国产精品精品国产色婷婷| 国产精品亚洲一级av第二区| 桃色一区二区三区在线观看| 香蕉av资源在线| 人人妻人人看人人澡| 国模一区二区三区四区视频| 嫩草影视91久久| 欧美国产日韩亚洲一区| 舔av片在线| 国产亚洲精品久久久久久毛片| 97超级碰碰碰精品色视频在线观看| 亚洲激情在线av| 亚洲无线观看免费| 嫁个100分男人电影在线观看| 最好的美女福利视频网| 国产 一区 欧美 日韩| 偷拍熟女少妇极品色| 国产精品久久久久久久电影 | 免费av毛片视频| 欧美日本视频| 禁无遮挡网站| 亚洲成人中文字幕在线播放| 国内精品久久久久精免费| 成人午夜高清在线视频| 日本免费一区二区三区高清不卡| svipshipincom国产片| 国产一区二区三区在线臀色熟女| 中文亚洲av片在线观看爽| 97超级碰碰碰精品色视频在线观看| 国产91精品成人一区二区三区| 9191精品国产免费久久| 午夜福利视频1000在线观看| 一级a爱片免费观看的视频| 脱女人内裤的视频| 久久久久久久精品吃奶| 国产三级中文精品| 中文在线观看免费www的网站| 国内毛片毛片毛片毛片毛片| 国产精品av视频在线免费观看| 久久久久久久久中文| a在线观看视频网站| 亚洲自拍偷在线| 国产探花在线观看一区二区| 一个人观看的视频www高清免费观看| 嫩草影视91久久| 国产综合懂色| 日日干狠狠操夜夜爽| ponron亚洲| 国产av一区在线观看免费| 国产精品女同一区二区软件 | 九九在线视频观看精品| 性欧美人与动物交配| 变态另类丝袜制服| 变态另类成人亚洲欧美熟女| 欧美黄色淫秽网站| 国产精品一区二区三区四区久久| 国产精品国产高清国产av| 熟女人妻精品中文字幕| 香蕉av资源在线| 国产探花极品一区二区| 亚洲人成电影免费在线| 婷婷六月久久综合丁香| 亚洲一区二区三区不卡视频| 麻豆成人午夜福利视频| 亚洲精品日韩av片在线观看 | 五月玫瑰六月丁香| 757午夜福利合集在线观看| 午夜激情欧美在线| 亚洲av免费在线观看| 久久久久国产精品人妻aⅴ院| 在线免费观看不下载黄p国产 | 男人舔女人下体高潮全视频| 国产亚洲精品综合一区在线观看| 午夜免费激情av| 国产一区二区在线观看日韩 | 中文字幕人妻丝袜一区二区| 99久久九九国产精品国产免费| 很黄的视频免费| 十八禁网站免费在线| 免费av不卡在线播放| 少妇的丰满在线观看| 色综合站精品国产| 亚洲人成网站在线播放欧美日韩| 成年女人毛片免费观看观看9| 亚洲国产欧美人成| 亚洲电影在线观看av| 欧美中文综合在线视频|