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

    ENSO 事件對(duì)山東蘋果生產(chǎn)的影響Ⅱ: 不同ENSO 年型下農(nóng)業(yè)氣象災(zāi)害變化及對(duì)蘋果產(chǎn)量的影響*

    2023-06-28 15:58:02劉布春孫彥坤張曉男劉珊珊朱永昶賀金娜
    關(guān)鍵詞:旱澇拉尼厄爾尼諾

    崔 成,劉 園,劉布春,孫彥坤,楊 凡,3,張曉男,4,劉珊珊,朱永昶,賀金娜

    (1.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所/作物高效用水與抗災(zāi)減損國(guó)家工程實(shí)驗(yàn)室/農(nóng)業(yè)部農(nóng)業(yè)環(huán)境重點(diǎn)實(shí)驗(yàn)室 北京 100081;2.東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院 哈爾濱 150036;3.新疆農(nóng)業(yè)大學(xué)食品科學(xué)與藥學(xué)學(xué)院 烏魯木齊 830000;4.張家口機(jī)場(chǎng)有限公司航務(wù)管理部 張家口 075000)

    蘋果(Malus pumila)是中國(guó)第一大水果,在我國(guó)林果產(chǎn)業(yè)中有著重要的地位,主要分布在環(huán)渤海、黃土高原、黃河故道和秦嶺北麓、新疆以及西南冷涼高地等地[1]。山東作為我國(guó)蘋果傳統(tǒng)環(huán)渤海優(yōu)勢(shì)主產(chǎn)區(qū)的主產(chǎn)省之一,具有獨(dú)特的地理位置和優(yōu)越的氣候資源。由于受臨海季風(fēng)、極端氣候的影響,農(nóng)業(yè)氣象災(zāi)害呈頻發(fā)、重發(fā)態(tài)勢(shì),給蘋果安全生產(chǎn)帶來(lái)了挑戰(zhàn)[2]。

    厄爾尼諾-南方濤動(dòng)事件(ENSO 事件)是指中東太平洋地區(qū)海溫異常變化與太平洋和印度洋地區(qū)東西方向海平面氣壓變化的蹺蹺板現(xiàn)象,是目前人類所能觀測(cè)到最強(qiáng)的全球尺度的海洋和大氣耦合系統(tǒng)[3]。20 世紀(jì)60 年代以來(lái)國(guó)內(nèi)外專家研究發(fā)現(xiàn),ENSO 事件常引發(fā)中國(guó)不同地區(qū)旱澇、高低溫等氣象災(zāi)害的發(fā)生,進(jìn)而影響農(nóng)業(yè)氣候資源的時(shí)空格局,對(duì)農(nóng)業(yè)生產(chǎn)造成一定的損失[4-8]。氣候變化背景下,環(huán)渤海地區(qū)旱澇災(zāi)害頻發(fā),ENSO 事件與降水強(qiáng)度、旱澇指數(shù)呈顯著負(fù)相關(guān),與旱澇等級(jí)呈顯著正相關(guān)[4-6];厄爾尼諾年降水減少且有明顯滯后性,拉尼娜年與之相反。ENSO 事件對(duì)山東氣溫影響呈空間差異性,總體而言,厄爾尼諾年極端高溫日數(shù)增加,極端低溫日數(shù)減少,不同時(shí)段氣溫與ENSO 事件存在不同共振周期[9];但苗正偉等[10]認(rèn)為京津冀地區(qū)極端溫度指數(shù)與ENSO 事件無(wú)顯著相關(guān)性。與此同時(shí),旱澇災(zāi)害的發(fā)生對(duì)環(huán)渤海地區(qū)蘋果生產(chǎn)影響較大[11]。降水量異常變化影響了蘋果產(chǎn)量和品質(zhì)[12-13]。程雪等[14]、楊建瑩等[15]分析了北方蘋果種植區(qū)干旱時(shí)空分布特征,指明了蘋果干旱的觸發(fā)閾值,并在區(qū)域進(jìn)行了驗(yàn)證。馬麗君[16]、程婷婷等[17]分析了夏季干旱、冬春連旱、低溫等農(nóng)業(yè)氣象災(zāi)害,發(fā)現(xiàn)這是導(dǎo)致蘋果坐果率低的主要原因,并提出相應(yīng)的田間管理措施建議。但分析大尺度環(huán)流異常的ENSO 事件對(duì)山東農(nóng)業(yè)氣象災(zāi)害、蘋果產(chǎn)量的影響鮮見(jiàn)報(bào)道。鑒于旱澇災(zāi)害、極端溫度災(zāi)害對(duì)蘋果生產(chǎn)影響較大,本文以此為切入點(diǎn),探究在大尺度環(huán)流異常ENSO 事件影響下的山東旱澇、極端溫度等農(nóng)業(yè)氣象災(zāi)害的區(qū)域演變特征,探討其對(duì)蘋果產(chǎn)量的影響,以期為保障山東蘋果優(yōu)勢(shì)產(chǎn)區(qū)的健康發(fā)展,規(guī)避農(nóng)業(yè)氣象災(zāi)害提供科學(xué)參考。

    1 材料與方法

    1.1 研究區(qū)域

    山東省位于我國(guó)東部環(huán)渤海沿海地區(qū),屬于溫帶季風(fēng)氣候,氣候資源季節(jié)分配不均,雨熱同期(圖1)。山東省也是我國(guó)重要的農(nóng)業(yè)生產(chǎn)基地之一[18-19]。

    圖1 研究區(qū)域及氣象站點(diǎn)分布圖Fig.1 Location of the study area and the distribution of meteorological stations

    1.2 研究數(shù)據(jù)

    本文選取山東省18 個(gè)國(guó)家標(biāo)準(zhǔn)氣象站(圖1) 1991-2019 年逐日氣象數(shù)據(jù),包括最高、最低和平均溫度,降水量,日照時(shí)數(shù)和風(fēng)速等;山東省17 個(gè)地市1991-2019 年蘋果種植面積、產(chǎn)量等統(tǒng)計(jì)數(shù)據(jù)。同時(shí),選取1991-2019 年逐月南方濤動(dòng)指數(shù)(Southern Oscillation Index,簡(jiǎn)稱SOI)數(shù)據(jù)[20-21],全球有且唯一的SOI 數(shù)據(jù)主要用于劃分大尺度環(huán)流異常ENSO 事件年型[19],統(tǒng)計(jì)結(jié)果如表1 所示。

    表1 1991-2019 年厄爾尼諾-南方濤動(dòng)(ENSO)事件不同年型的統(tǒng)計(jì)[20-21]Table 1 Classification of different El Ni?o-Southern Oscillation (ENSO) events from 1991 to 2019

    1.3 研究方法

    1.3.1 山東蘋果可生長(zhǎng)期及生長(zhǎng)的基本氣候條件

    根據(jù)查閱相關(guān)文獻(xiàn)[19]和實(shí)地走訪,收集整理了山東各地不同品種蘋果實(shí)際生長(zhǎng)發(fā)育日期。本文以可覆蓋蘋果所有品種的關(guān)鍵生育期為目標(biāo),劃定每年3 月1 日-10 月31 日為山東蘋果的可生長(zhǎng)期(潛在生長(zhǎng)期),即在此時(shí)間范圍內(nèi)的氣候條件可滿足蘋果正常生長(zhǎng)發(fā)育的最低要求;同時(shí),定義3 月1 日-5 月31 日為開(kāi)花期,6 月1 日-8 月31 日為果實(shí)膨大期。

    1.3.2 旱澇災(zāi)害的判定

    參考國(guó)家標(biāo)準(zhǔn)《氣象干旱等級(jí)》(GB/T 20481-2017)及相關(guān)文獻(xiàn)[22-24],本文選取降水距平百分率作為山東旱澇災(zāi)害發(fā)生的判定指標(biāo)。降水距平百分率是描述某時(shí)段內(nèi)某地降水量與同期降水量多年平均狀態(tài)的偏離程度,用百分率(%)表示,降水量越大降水距平偏離氣候平均值越大,發(fā)生干旱或雨澇災(zāi)害的風(fēng)險(xiǎn)越大。計(jì)算公式如下:

    式中:Pa為降水量距平(%),P為某時(shí)段降水量(mm),為同期降水量的多年平均值(mm)。以降水距平百分率為指標(biāo)的旱澇等級(jí)劃分,具體如表2 所示[23-24]。

    表2 根據(jù)降水距平百分率(Pa)劃分旱澇等級(jí)Table 2 Drought and flood grades according to the percentage of precipitation anomaly (Pa)

    1.3.3 旱澇站次比

    旱澇站次比是研究區(qū)域內(nèi)同一時(shí)段發(fā)生干旱或雨澇的站點(diǎn)與區(qū)域內(nèi)全部站點(diǎn)的比值,用于表達(dá)干旱或雨澇災(zāi)害的影響區(qū)域大小。計(jì)算公式如下:

    式中:Pj為旱澇站次比(%),M為研究區(qū)域內(nèi)站點(diǎn)總個(gè)數(shù),m為發(fā)生干旱或雨澇事件的站點(diǎn)數(shù)。以旱澇站次比為指標(biāo)的旱澇等級(jí)劃分[22-25],具體如表3 所示。

    表3 以旱澇站次比(Pj)劃分區(qū)域?yàn)?zāi)害影響范圍Table 3 Influence area of regional disasters by the stations ratio of drought or flooding (Pj)

    1.3.4 旱澇頻率

    研究區(qū)域內(nèi)某時(shí)段干旱或雨澇發(fā)生的頻率用Pi表示,計(jì)算公式如下:

    式中:N為研究時(shí)間段,n為發(fā)生干旱或雨澇的次數(shù)[2]。

    1.3.5 極端溫度

    山東省蘋果極端溫度指標(biāo)如表4 所示。每年3-5 月為蘋果開(kāi)花階段,常受低溫、倒春寒等災(zāi)害的影響,開(kāi)花受到低溫影響,影響蘋果坐果率和果實(shí)品質(zhì)[26-27]。低溫災(zāi)害一般分為低溫冷害和低溫凍害。6-8 月為果實(shí)膨大期,高溫?zé)岷?dǎo)致果實(shí)發(fā)生日灼,影響果實(shí)品質(zhì)和產(chǎn)量[25]。

    表4 山東蘋果極端溫度指標(biāo)Table 4 Extreme temperature indexes of apple in Shandong

    1.3.6 產(chǎn)量分離

    為研究蘋果產(chǎn)量與氣象因子之間的關(guān)系,剔除生產(chǎn)力發(fā)展下的趨勢(shì)產(chǎn)量,得到蘋果的氣候產(chǎn)量,具體公式如下:

    式中:Y為實(shí)際產(chǎn)量(kg·hm-2);Yi為趨勢(shì)產(chǎn)量(kg·hm-2);Y0為氣象產(chǎn)量(kg·hm-2);a為隨機(jī)誤差產(chǎn)量,忽略不計(jì)。本文選取5 年滑動(dòng)平均法計(jì)算蘋果產(chǎn)量的趨勢(shì)產(chǎn)量[28]。

    此外,利用蘋果減產(chǎn)率(Yw)來(lái)判斷氣象條件是否有利于蘋果產(chǎn)量。當(dāng)Y0>0 時(shí),蘋果減產(chǎn)率為0;當(dāng)Y0<0 時(shí),即氣候條件易導(dǎo)致蘋果產(chǎn)量減產(chǎn),取絕對(duì)值后記為蘋果減產(chǎn)率。具體計(jì)算如下[29]:

    2 結(jié)果與分析

    2.1 不同ENSO 年型蘋果果實(shí)膨大期(6-8 月)旱澇災(zāi)害特征

    2.1.1 旱澇頻次/頻率

    6-8 月是蘋果果實(shí)膨大期,降水量多寡影響蘋果果實(shí)著色及糖分積累[30]。1991-2019 年,不同ENSO年型下山東干旱、雨澇具有時(shí)空差異性(圖2)。厄爾尼諾年,山東共發(fā)生干旱站次共78 次,雨澇26 次;拉尼娜年,發(fā)生干旱站次有32 次,雨澇29 次;中性年,干旱站次僅發(fā)生6 次,雨澇站次發(fā)生了60 次。厄爾尼諾年干旱發(fā)生頻次高于雨澇發(fā)生頻次,中性年雨澇發(fā)生頻次高于干旱。從空間來(lái)看,厄爾尼諾年膠東半島、魯南地區(qū)干旱發(fā)生站次分別為26 次和16 次,雨澇發(fā)生4 次和2 次;拉尼娜年干旱站次發(fā)生高值區(qū)分布在膠東半島(16 次)和魯北地區(qū)(7 次);中性年,海陽(yáng)、煙臺(tái)、莒縣等地雨澇站次均為5 次??傮w而言,厄爾尼諾年更易發(fā)生干旱,中性年更易發(fā)生雨澇災(zāi)害。在地區(qū)分布上,魯西、魯中地區(qū)容易發(fā)生干旱災(zāi)害,膠東半島、魯南地區(qū)容易發(fā)生雨澇災(zāi)害。

    圖2 不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型山東6-8 月降水距平百分率及旱澇等級(jí)(1991-2019 年)Fig.2 Percentage of precipitation anomalies and grade of droughts and floods from June to August under different El Ni?o-Southern Oscillation (ENSO) years in Shandong (1991-2019)

    不同ENSO 年型下山東旱澇災(zāi)害發(fā)生頻率如圖3所示。1991-2019 年,山東地區(qū)干旱和旱澇發(fā)生頻率基本持平;厄爾尼諾年干旱發(fā)生頻率約35%,高于雨澇發(fā)生頻率(10%);中性年雨澇發(fā)生頻率集中在35%~45%,最高達(dá)55%,干旱頻率為0~10%,高值約在20%;拉尼娜年干旱發(fā)生頻率為15%~45%,雨澇發(fā)生頻率在0~45%??傮w而言,厄爾尼諾年干旱發(fā)生頻率高于拉尼娜年和中性年,中性年雨澇災(zāi)害發(fā)生頻率高于厄爾尼諾年和拉尼娜年。

    圖3 不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型山東6-8 月旱澇頻率分布Fig.3 Frequency of drought and flood from June to August under different El Ni?o-Southern Oscillation (ENSO)years in Shandong

    2.1.2 旱澇站次比

    1991-2019 年山東旱澇站次比的變化如圖4 所示,過(guò)去近30 年,山東干旱發(fā)生站次比多于雨澇發(fā)生站次比。山東全區(qū)共發(fā)生1 次全域性干旱、2 次區(qū)域性干旱及3 次區(qū)域性雨澇。厄爾尼諾年山東發(fā)過(guò)1 次全區(qū)干旱、2 次區(qū)域性干旱,中性年發(fā)生3 次區(qū)域性雨澇。厄爾尼諾年局部性雨澇、拉尼娜年局地性雨澇高于局地性干旱的發(fā)生。

    圖4 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下山東6-8 月旱澇災(zāi)害站次比Fig.4 Stations ratio of drought and flood from June to August under different El Ni?o-Southern Oscillation (ENSO) years in Shandong (1991-2019)

    2.2 極端溫度

    2.2.1 低溫冷害

    過(guò)去近30 年,山東發(fā)生低溫冷害的平均年日數(shù)來(lái)看(圖5a-c),魯東、魯北地區(qū)屬低溫冷害高發(fā)區(qū),魯西地區(qū)屬低溫冷害發(fā)生低值區(qū)。厄爾尼諾年低溫冷害平均年日數(shù)高值區(qū)在膠東半島和魯北地區(qū),約5~7 d·a-1,其中濰坊低溫冷害的平均年日數(shù)達(dá)7~9 d·a-1;魯西、魯南煙臺(tái)、威海地區(qū)低溫冷害發(fā)生平均年日數(shù)為3~5 d·a-1。拉尼娜年,低溫冷害低值區(qū)在濟(jì)南、朝陽(yáng)、菏澤、日照、青島等地區(qū),平均為1~3 d·a-1。中性年,低溫冷害高值區(qū)在魯北、魯中及山東半島地區(qū),平均為5~7 d·a-1,德州、濰坊、煙臺(tái)低溫冷害發(fā)生的年均日數(shù)最高,平均為7~9 d·a-1。中性年低溫冷害的高值區(qū)范圍相對(duì)最大。

    圖5 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下低溫冷害(-2 ℃<極端低溫≤0 ℃)發(fā)生日數(shù)、> 3 d 低溫冷害頻率空間分布Fig.5 Spatial distribution of occurrence days of cold injury (-2 ℃ < extreme low temperature≤ 0 ℃) and frequency of cold injury lasting for more than 3 days under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019

    在蘋果開(kāi)花期發(fā)生持續(xù)3 d 以上的低溫冷害,影響蘋果花粉活力和坐果率[27]。圖5d-f 是>3 d 低溫冷害發(fā)生的頻率區(qū)域分布,>3 d 低溫冷害的高頻率區(qū)在魯中至膠東半島地區(qū),頻率約在60%~100%;魯西、魯東南地區(qū)為低頻率區(qū)。厄爾尼諾年魯西、魯東南及煙臺(tái)地區(qū)>3 d 低溫冷害的發(fā)生頻率在30%~60%。拉尼娜年,>3 d 低溫冷害低值區(qū)(10%~30%)分布在濟(jì)南、日照、青島地區(qū)。中性年>3 d 低溫冷害發(fā)生頻率低值區(qū)(30%~60%)分布在菏澤、濟(jì)南、日照等地,高頻率分布范圍最廣,其中膠東半島地區(qū)低溫冷害發(fā)生頻率為90%~100%。綜上所述,厄爾尼諾年下,發(fā)生低溫冷害的日數(shù)、頻率較低,受低溫冷害的影響較低;而中性年,低溫冷害發(fā)生持續(xù)日數(shù)較長(zhǎng),且發(fā)生頻率的影響范圍最廣,因此中性年受低溫災(zāi)害影響較大。

    2.2.2 低溫凍害

    1991-2019 年,山東蘋果開(kāi)花期低溫凍害平均年日數(shù)的空間分布如圖6a-c 所示。魯東、膠東半島地區(qū)為低溫凍害高值區(qū),魯西地區(qū)為低值區(qū)。相較于低溫冷害發(fā)生情況,低溫凍害影響區(qū)域偏高。厄爾尼諾年魯東及膠東半島地區(qū)低溫凍害年日數(shù)5~7 d·a-1,濟(jì)南、菏澤、日照、青島地區(qū)平均為1~3 d·a-1。拉尼娜年,魯西、魯東南地區(qū)及東營(yíng)地區(qū)為低值區(qū)。中性年,魯東至膠東半島地區(qū)為低溫凍害發(fā)生的高值區(qū)。3 種ENSO 年型下,濰坊地區(qū)低溫凍害年發(fā)生日數(shù)最高,為7~9 d·a-1。

    圖6 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下低溫凍害(極端低溫≤-2 ℃)發(fā)生日數(shù)、>3 d 低溫凍害頻率空間分布Fig.6 Spatial distribution of low temperature freezing injury (extreme low temperature ≤ -2 ℃) days and frequency of low temperature freezing injury lasting for more than 3 days under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019

    1991-2019 年,>3 d 低溫凍害發(fā)生頻率的空間分布如圖6d-f 所示。厄爾尼諾年魯東、膠東半島地區(qū)>3 d 低溫凍害發(fā)生頻率較高,約在60%~90%,其中濰坊地區(qū)>3 d 低溫凍害頻率最高,為90%~100%,厄爾尼諾年下,近乎年年發(fā)生。拉尼娜年、中性年,在菏澤、濟(jì)南、東營(yíng)等地低溫凍害發(fā)生頻率較低,約在10%~30%。綜上而言,厄爾尼諾年低溫凍害在魯東、魯北等地區(qū)分布范圍最廣,發(fā)生頻率也較高。由于濰坊地區(qū)以東、以北為渤海、沿海平原,以南、以西多為山地,冷空氣易進(jìn)難出,使得低溫災(zāi)害發(fā)生重于其他地區(qū)。

    2.2.3 高溫?zé)岷?/p>

    6-8 月是蘋果果實(shí)膨大期,高溫?zé)岷τ绊懱O果果實(shí)品質(zhì),也易引發(fā)病害發(fā)生[30]。在1991-2019 年山東發(fā)生高溫?zé)岷Φ钠骄耆諗?shù)區(qū)域分布如圖7a-c所示。不同ENSO 年型山東發(fā)生高溫?zé)岷δ耆諗?shù)呈東向西增加的空間分布。厄爾尼諾年,高溫?zé)岷ζ骄耆諗?shù)高值區(qū)分布在魯中、魯西等地,其中魯西地區(qū)發(fā)生天數(shù)最高,約在11~15 d·a-1;拉尼娜年、中性年在魯東、膠東半島地區(qū)平均年日數(shù)較低,約在0~5 d·a-1。

    圖7 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下高溫?zé)岷?極端高溫≥35 ℃)發(fā)生日數(shù)空間、> 3 d 高溫?zé)岷︻l率分布Fig.7 Spatial distribution of occurrence days of high temperature (extereme high temperature ≥35 ℃) and frequency of high temperature lasting for more than 3 days under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019

    1991-2019 年山東發(fā)生高溫?zé)岷︻l率的空間分布如圖7d-f 所示,厄爾尼諾年,魯東、膠東半島地區(qū)>3 d 發(fā)生頻率僅有0~30%;在魯北、魯中等地>3 d發(fā)生頻率較高,約在60%~90%,其中魯西地區(qū)頻率最高,在90%~100%。拉尼娜年,僅魯西、魯北地區(qū)>3 d發(fā)生頻率較高,約60%~90%。中性年,魯西、魯北地區(qū)>3 d 發(fā)生頻率約在90%~100%。魯東、膠東半島地區(qū)臨近海洋,受海洋氣候的調(diào)節(jié)作用,高溫災(zāi)害有所減弱。

    2.3 ENSO 事件與農(nóng)業(yè)氣象災(zāi)害的關(guān)系

    2.3.1 旱澇災(zāi)害與ENSO 的相關(guān)性

    為了探明區(qū)域旱澇關(guān)系與大尺度環(huán)流的關(guān)系,進(jìn)行了山東地區(qū)蘋果旱澇災(zāi)害指數(shù)與不同ENSO 年型的南方濤動(dòng)指數(shù)(SOI)值相關(guān)關(guān)系的分析,結(jié)果如圖8 所示。3-10 月蘋果可生長(zhǎng)期內(nèi),厄爾尼諾年SOI 與山東大部分地區(qū)干旱呈正相關(guān),拉尼娜年SOI與山東大部分地區(qū)干旱呈負(fù)相關(guān),與膠東半島和魯東南地區(qū)干旱呈正相關(guān)。中性年SOI 與魯東南、膠東半島北部干旱呈負(fù)相關(guān),與魯中地區(qū)干旱呈顯著正相關(guān)(通過(guò)顯著性90%檢驗(yàn))。果實(shí)膨大期(6-8 月),厄爾尼諾年SOI 與魯中、魯西等地干旱呈正相關(guān)關(guān)系;拉尼娜年SOI 與山東大部分地區(qū)呈負(fù)相關(guān),中性年SOI 與膠東半島干旱呈負(fù)相關(guān);其中,拉尼娜年、中性年SOI 與魯東地區(qū)干旱呈顯著負(fù)相關(guān)關(guān)系(通過(guò)顯著性90%檢驗(yàn))。

    圖8 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下旱澇災(zāi)害(A: 3-10 月干旱;B: 3-10 月雨澇;C: 6-8月干旱;D: 6-8 月雨澇)與ENSO 事件南方濤動(dòng)指數(shù)相關(guān)性的空間分布Fig.8 Spatial distribution of correlation between drought and flood disasters (A: drought from March to October;B: flood from March to October;C: drought from June to August;D: flood from June to August) and El Ni?o-Southern Oscillation (ENSO)events in different ENSO years from 1991 to 2019

    對(duì)于雨澇災(zāi)害而言,3-10 月蘋果可生長(zhǎng)期,厄爾尼諾年SOI 數(shù)與山東大部分地區(qū)雨澇呈正相關(guān),與膠東半島雨澇呈負(fù)相關(guān);拉尼娜年、中性年SOI與雨澇呈負(fù)相關(guān)。果實(shí)膨大期(6-8 月),厄爾尼諾年SOI 與雨澇呈負(fù)相關(guān)關(guān)系,與魯東、魯西部分地區(qū)干旱呈顯著負(fù)相關(guān)(通過(guò)90%顯著性檢驗(yàn))。拉尼娜年SOI 與魯中、魯西地區(qū)雨澇呈正相關(guān);與魯東、膠東半島地區(qū)雨澇呈負(fù)相關(guān)。中性年SOI 與山東大部分地區(qū)雨澇呈正相關(guān),但與魯東、魯東南地區(qū)雨澇呈負(fù)相關(guān)。

    2.3.2 極端溫度災(zāi)害與ENSO 事件的相關(guān)性

    1991-2019 年,山東低溫災(zāi)害發(fā)生天數(shù)與不同ENSO 年型的SOI 相關(guān)關(guān)系的空間分布如圖9 所示。厄爾尼諾年SOI 與山東大部分地區(qū)低溫冷害、凍害呈負(fù)相關(guān)關(guān)系,與魯中地區(qū)相關(guān)性通過(guò)90%顯著性檢驗(yàn)。拉尼娜年、中性年SOI 與山東低溫冷害、凍害呈正相關(guān)關(guān)系;拉尼娜年SOI 與魯西地區(qū)、中性年SOI 與魯中地區(qū)低溫冷害、低溫凍害呈負(fù)相關(guān)。總體來(lái)說(shuō),厄爾尼諾年低溫災(zāi)害受SOI 影響較小;拉尼娜年、中性年隨著SOI 增強(qiáng),低溫冷害、凍害發(fā)生日數(shù)有所增加。

    圖9 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下低溫冷害(A)和低溫凍害(B)與ENSO 事件相關(guān)性的空間分布Fig.9 Spatial distribution of correlation between chilling damage (A),low temperature freezing injury (B) and El Ni?o-Southern Oscillation (ENSO) events in different ENSO years from 1991 to 2019

    1991-2019 年,山東蘋果果實(shí)膨大期高溫災(zāi)害發(fā)生天數(shù)與不同ENSO 年型SOI 相關(guān)關(guān)系的空間分布如圖10 所示,山東大部分地區(qū)高溫?zé)岷εcENSO 事件SOI 呈負(fù)相關(guān)。厄爾尼諾年SOI 與膠東半島,拉尼娜年SOI 與沂源、日照、威海等地,中性年SOI 與龍口、濰坊、日照等地高溫?zé)岷Τ收嚓P(guān),但相關(guān)關(guān)系未通過(guò)90%顯著性檢驗(yàn)??傮w而言,隨著南方SOI 變化,對(duì)山東高溫?zé)岷τ绊懹邢蕖?/p>

    圖10 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下高溫?zé)岷εcENSO 事件相關(guān)性的空間分布Fig.10 Spatial distribution of correlation between high temperature heat damage and El Ni?o-Southern Oscillation (ENSO) events under different ENSO years during 1991-2019

    2.4 農(nóng)業(yè)氣象災(zāi)害與蘋果減產(chǎn)率的關(guān)系

    2.4.1 旱澇災(zāi)害與蘋果減產(chǎn)率的關(guān)系

    不同ENSO 年型,山東蘋果減產(chǎn)率與旱澇災(zāi)害相關(guān)關(guān)系的空間分布如圖11 所示。3-10 月蘋果可生長(zhǎng)期,厄爾尼諾年,山東魯東北、膠東半島干旱與蘋果減產(chǎn)率呈正相關(guān),隨著干旱程度增強(qiáng),蘋果減產(chǎn)率也隨之增加;其余地區(qū)干旱與蘋果減產(chǎn)率呈負(fù)相關(guān),其中濟(jì)寧、濰坊干旱與蘋果減產(chǎn)率呈顯著負(fù)相關(guān)(通過(guò)90%顯著性檢驗(yàn))。拉尼娜年德州、濰坊地區(qū),中性年德州、濱州、菏澤、青島地區(qū)蘋果減產(chǎn)率與干旱呈顯著正相關(guān)(通過(guò)顯著性90%)。拉尼娜年魯西南、魯東北及中性年濟(jì)南、臨沂地區(qū)干旱災(zāi)害與蘋果減產(chǎn)率呈負(fù)相關(guān),干旱災(zāi)害與蘋果減產(chǎn)率的相關(guān)系數(shù)較小。對(duì)3-10 月蘋果可生長(zhǎng)期雨澇災(zāi)害而言,山東大部分地區(qū)厄爾尼諾年膠東半島雨澇災(zāi)害與蘋果減產(chǎn)率呈正相關(guān)關(guān)系,其余地區(qū)蘋果減產(chǎn)率受雨澇影響較低;中性年山東大部分地區(qū)雨澇與蘋果減產(chǎn)率呈正相關(guān),降水增多導(dǎo)致蘋果減產(chǎn)率增加。

    圖11 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下旱澇災(zāi)害(降水距平百分率)與蘋果減產(chǎn)率相關(guān)性的空間分布(A: 3-10 月干旱;B: 3-10 月雨澇;C: 6-8 月干旱;D: 6-8 月雨澇)Fig.11 Spatial distribution of correlation between drought and flood disasters (percentage of precipitation anomaly) and apple yield reduction rate under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019 (A: drought from March to October;B: flood from March to October;C: drought from June to August;D: flood from June to August)

    6-8 月蘋果果實(shí)膨大期,厄爾尼諾年菏澤、魯東北及膠東半島地區(qū)干旱與蘋果減產(chǎn)率呈正相關(guān),尤其東營(yíng)、淄博、青島地區(qū)相關(guān)關(guān)系通過(guò)顯著性90%檢驗(yàn);濰坊、臨沂等地干旱與蘋果減產(chǎn)率呈顯著負(fù)相關(guān)(通過(guò)顯著性90%檢驗(yàn))。拉尼娜年膠東半島、淄博、日照、菏澤地區(qū)干旱災(zāi)害與蘋果減產(chǎn)率呈顯著負(fù)相關(guān)(通過(guò)顯著性90%檢驗(yàn)),濰坊、臨沂地區(qū)干旱與蘋果減產(chǎn)率呈顯著正相關(guān)(通過(guò)顯著性90%檢驗(yàn))。中性年魯西、魯東北、煙臺(tái)、日照等地干旱與蘋果減產(chǎn)率呈顯著負(fù)相關(guān)(通過(guò)顯著性90%檢驗(yàn)),濱州、臨沂地區(qū)干旱與蘋果減產(chǎn)率呈顯著正相關(guān)(通過(guò)90%顯著性檢驗(yàn))。厄爾尼諾年、拉尼娜年雨澇災(zāi)害與蘋果減產(chǎn)率呈負(fù)相關(guān),中性年絕大部分地區(qū)雨澇災(zāi)害與蘋果減產(chǎn)率呈正相關(guān)。這與不同ENSO 年型下氣候條件差異有著密切的關(guān)聯(lián)。

    2.4.2 極端氣溫與蘋果減產(chǎn)率的關(guān)系

    1991-2019 年,不同ENSO 年型低溫災(zāi)害發(fā)生天數(shù)與蘋果減產(chǎn)率關(guān)系的空間分布如圖12 所示。由圖可知,厄爾尼諾年魯北、威海等地低溫冷害與蘋果減產(chǎn)率呈正相關(guān),隨低溫冷害日數(shù)的增加,蘋果減產(chǎn)減收幅度增大;其中德州、威海地區(qū)低溫冷害與蘋果減產(chǎn)率呈顯著性正相關(guān)(通過(guò)顯著性90%檢驗(yàn))。拉尼娜年魯東、魯東北及中性年山東半島地區(qū)低溫冷害與蘋果減產(chǎn)率呈負(fù)相關(guān),低溫冷害對(duì)蘋果減產(chǎn)率影響較小。厄爾尼諾年山東大部分地區(qū)低溫凍害與蘋果減產(chǎn)率呈正相關(guān),僅膠東半島、德州等地呈負(fù)相關(guān),煙臺(tái)地區(qū)低溫凍害與蘋果減產(chǎn)率呈顯著負(fù)相關(guān)(通過(guò)顯著性90%檢驗(yàn))。拉尼娜年僅魯西、魯南及膠東半島部分地區(qū)低溫凍害與蘋果減產(chǎn)率呈正相關(guān)關(guān)系。中性年僅魯北及魯南地區(qū)低溫凍害對(duì)蘋果減產(chǎn)影響較低。

    圖12 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下低溫冷害(A)和低溫凍害(B)與蘋果減產(chǎn)率相關(guān)性的空間分布Fig.12 Spatial distribution of correlation between chilling damage (A),low termperature freezing injury (B) and apple yield reduction rate under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019

    不同ENSO 年型下高溫?zé)岷Πl(fā)生天數(shù)與蘋果減產(chǎn)率相關(guān)性的空間分布如圖13 所示,厄爾尼諾年魯東及威海地區(qū)高溫?zé)岷εc蘋果減產(chǎn)率呈負(fù)相關(guān);其余地區(qū)呈正相關(guān)關(guān)系,即隨高溫?zé)岷θ諗?shù)的增加,蘋果減產(chǎn)減收越大。拉尼娜年膠東半島、魯西地區(qū)高溫?zé)岷εc蘋果減產(chǎn)率呈負(fù)相關(guān)。中性年僅魯東南及濟(jì)南地區(qū)高溫?zé)岷εc蘋果減產(chǎn)率呈正相關(guān),濟(jì)南地區(qū)正相關(guān)關(guān)系通過(guò)顯著性90%檢驗(yàn)。

    圖13 1991-2019 年不同厄爾尼諾-南方濤動(dòng)(ENSO)事件年型下高溫?zé)岷εc蘋果減產(chǎn)率相關(guān)性的空間分布Fig.13 Spatial distribution of correlation between high temperature and apple yield reduction rate under different El Ni?o-Southern Oscillation (ENSO) years during 1991-2019

    3 討論

    降水距平百分率指數(shù)能較好地反映區(qū)域旱澇的時(shí)空變化特征。1991-2019 年,不同ENSO 年型的山東旱澇災(zāi)害變化差異明顯。由于厄爾尼諾年氣候資源整體趨于暖干性,熱量資源較為豐富[19],厄爾尼諾年易發(fā)生干旱;中性年降水資源充沛,易發(fā)生雨澇,是影響中性年雨澇災(zāi)害頻發(fā)的重要因素之一[19];雖然拉尼娜年降水量較為豐富,但干旱發(fā)生頻次仍高于雨澇。從空間分布來(lái)看,魯南地區(qū)雨澇災(zāi)害發(fā)生的頻次高于干旱,尤其在拉尼娜年和中性年;而魯中、魯西地區(qū)干旱發(fā)生頻次高于雨澇頻次,主要是因?yàn)閰^(qū)域地理位置的差異性,導(dǎo)致不同區(qū)域的熱量資源和降水資源的差異。在不同ENSO 年型下,氣候資源分配也具有差異性。魯南地區(qū)臨近海洋,受季風(fēng)氣候影響較為顯著,加之拉尼娜年、中性年的降水資源偏多,導(dǎo)致該區(qū)域雨澇過(guò)多。而魯中、魯西地處內(nèi)陸,該區(qū)域熱量資源豐富、氣溫偏高,在厄爾尼諾年下降水偏少,干旱發(fā)生較為頻繁[19]。因此,魯西、魯中地區(qū)為干旱災(zāi)害的易發(fā)、重發(fā)地區(qū),尤其是厄爾尼諾年;魯南等地則成為雨澇災(zāi)害的高發(fā)區(qū),在中性年,魯南地區(qū)雨澇災(zāi)害尤為嚴(yán)重。從影響范圍來(lái)看,區(qū)域性干旱多發(fā)生厄爾尼諾年,中性年、拉尼娜年多發(fā)生局部或者地區(qū)性雨澇災(zāi)害。

    1991-2019 年,山東3-10 月蘋果可生長(zhǎng)期的干旱與厄爾尼諾年呈正相關(guān),但與拉尼娜年呈負(fù)相關(guān);山東大部分地區(qū)干旱與中性年呈正相關(guān)。厄爾尼諾年與雨澇呈正相關(guān),但拉尼娜年、中性年與雨澇呈負(fù)相關(guān)。6-8 月果實(shí)膨大期,魯中、魯西干旱與厄爾尼諾年呈正相關(guān);山東大部分地區(qū)干旱與拉尼娜年呈負(fù)相關(guān);其中拉尼娜年、中性年下,魯東地區(qū)與干旱呈顯著負(fù)相關(guān)(通過(guò)顯著性90%檢驗(yàn))。厄爾尼諾年與6-8 月雨澇災(zāi)害呈負(fù)相關(guān),但魯中、魯西地區(qū)雨澇與拉尼娜年呈正相關(guān);山東大部分地區(qū)雨澇與中性年呈正相關(guān)。厄爾尼諾年多發(fā)易發(fā)干旱災(zāi)害,而中性年卻頻發(fā)雨澇災(zāi)害,這與不同ENSO 年型下氣候條件的差異有關(guān)。厄爾尼諾年熱量資源充沛,但降水資源相較于拉尼娜年、中性年較少。由于地理位置的差異性,導(dǎo)致各地區(qū)氣候資源存在差異,這使得魯西、魯中等地成為了干旱的高發(fā)區(qū),且因?yàn)檫h(yuǎn)離海洋,氣候調(diào)節(jié)能力有限,使該地受大尺度氣候影響比沿海地區(qū)更為敏感??傮w而言,1991-2019 年,受ENSO 事件影響,熱量資源相對(duì)優(yōu)越的厄爾尼諾年,干旱災(zāi)害發(fā)生的程度、影響范圍高于其他年型,魯西、魯中等地區(qū)尤為突出;降水資源充沛的中性年,魯南地區(qū)雨澇災(zāi)害發(fā)生頻率和范圍較大。蘋果可生長(zhǎng)期內(nèi),厄爾尼諾年應(yīng)多加防范區(qū)域性干旱災(zāi)害,中性年謹(jǐn)防雨澇災(zāi)害。徐澤華等[2]認(rèn)為ENSO 事件下山東降水減少、溫度升高,致使山東氣候趨于暖干化,加劇了山東干旱災(zāi)害的發(fā)生;魯西及魯西北等地干旱較為嚴(yán)重;通過(guò)標(biāo)準(zhǔn)化蒸散指數(shù)分析山東旱澇災(zāi)害,發(fā)現(xiàn)厄爾尼諾年指數(shù)趨于干旱,且程度大于拉尼娜年。這一結(jié)果與本文研究一致。由于地理位置差異,使魯東、魯西兩地旱澇災(zāi)害發(fā)生情況有所差異,不同ENSO 年型下氣候分配不均勻,使得在相對(duì)溫暖干燥的厄爾尼諾年,魯西地區(qū)干旱加劇,而雨量較多的拉尼娜年、中性年,則導(dǎo)致魯東地區(qū)雨澇災(zāi)害頻發(fā)。

    1991-2019 年,不同ENSO 年型的山東蘋果開(kāi)花期低溫、果實(shí)膨大期的高溫災(zāi)害也存在時(shí)空變化差異。中性年低溫冷害、凍害發(fā)生的日數(shù)及發(fā)生頻率高于厄爾尼諾年和拉尼娜年,且影響范圍大。中性年低溫冷害在膠東半島地區(qū)發(fā)生頻率為90%~100%,發(fā)生日數(shù)在7~9 d·a-1,>3 d 的低溫凍害發(fā)生頻率在60%~90%。這與中性年降水資源過(guò)于豐富,且熱量資源相較于其他年份偏低有關(guān)。膠東半島低溫災(zāi)害的高發(fā)率,也與其地理位置有關(guān)。高溫?zé)岷t恰恰與低溫災(zāi)害相反,少有發(fā)生低溫災(zāi)害的魯西、魯中地區(qū)成為了高溫?zé)岷Φ母甙l(fā)區(qū),尤其是厄爾尼諾年影響范圍最大,持續(xù)天數(shù)也較多(約8~15 d·a-1),厄爾尼諾年魯西、魯中地區(qū)熱量資源優(yōu)越,也是其成為高溫?zé)岷Ω甙l(fā)區(qū)的重要原因。

    1991-2019 年,低溫冷害、凍害與厄爾尼諾年呈負(fù)相關(guān),與拉尼娜年、中性年呈正相關(guān);高溫?zé)岷εcENSO 事件均呈負(fù)相關(guān)。厄爾尼諾年,低溫冷害、低溫凍害與SOI 值呈負(fù)相關(guān),低溫災(zāi)害受ENSO 事件影響較小。拉尼娜年、中性年,低溫災(zāi)害與SOI值呈正相關(guān),低溫冷害、凍害發(fā)生天數(shù)有所增加。在魯中、膠東半島地區(qū)極端溫度變化受ENSO 事件影響較大,發(fā)生天數(shù)有所增加。董旭光等[31]、劉玄等[32]研究發(fā)現(xiàn)山東極端氣溫分布空間不均勻,魯東南、膠東半島等沿海地區(qū)與內(nèi)陸地區(qū)的極端熱氣溫指數(shù)變化差異明顯;姜德娟等[33]發(fā)現(xiàn)山東極端溫度發(fā)生的頻率空間差異較突出,膠東半島地區(qū)為氣候變化敏感區(qū),極端溫度事件發(fā)生更為頻繁且變化幅度較大。這與本文結(jié)果一致。在膠東半島、魯東南等地,高溫?zé)岷Πl(fā)生的天數(shù)、頻率遠(yuǎn)低于位于內(nèi)陸的魯西、魯中地區(qū)。這可能由于山東地理位置等因素的影響。因此魯西等地,尤為厄爾尼諾年,謹(jǐn)防高溫?zé)岷?lái)的危害。山東局部地區(qū)的ENSO 事件表現(xiàn)與相鄰地區(qū)變化有所差異,由于數(shù)據(jù)收集的有限性,可能是地區(qū)上在地理位置、地勢(shì)等因素的存在差異所導(dǎo)致。后續(xù)會(huì)進(jìn)一步收集相關(guān)數(shù)據(jù)研究。

    氣候變化背景下不同ENSO 年型旱澇災(zāi)害對(duì)蘋果減產(chǎn)率造成一定的影響。在3-10 月干旱災(zāi)害中,厄爾尼諾年,膠東半島地區(qū)干旱頻發(fā)導(dǎo)致蘋果減產(chǎn)率上升;中性年,山東魯北、魯西等地區(qū)旱災(zāi)與蘋果減產(chǎn)率呈正相關(guān),即隨著干旱程度的加深,蘋果減產(chǎn)減收的程度也增大。3-10 月雨澇,膠東半島在厄爾尼諾年與蘋果減產(chǎn)率呈正相關(guān);而中性年,山東大部分地區(qū)與蘋果減產(chǎn)率呈正相關(guān)。6-8 月干旱,厄爾尼諾年膠東半島與蘋果減產(chǎn)率呈正相關(guān),而拉尼娜年、中性年僅在魯中部分地區(qū)與蘋果減產(chǎn)率呈正相關(guān)。6-8 月雨澇,厄爾尼諾年、拉尼娜年山東大部分地區(qū)與蘋果減產(chǎn)率呈負(fù)相關(guān);但是中性年卻與蘋果減產(chǎn)率呈正相關(guān)。膠東半島地區(qū)由于緊鄰海洋,同時(shí)其蘋果種植面積最廣,蘋果減產(chǎn)減收更容易受到旱澇災(zāi)害的影響。魯西、魯中地區(qū)深居內(nèi)陸,熱量資源充沛(尤在厄爾尼諾年),旱災(zāi)的加劇,在一定程度上導(dǎo)致蘋果產(chǎn)量受災(zāi)損失加重。中性年降雨過(guò)多下,雨澇災(zāi)害頻發(fā)、重發(fā),使得蘋果減產(chǎn)減收加劇。因此,厄爾尼諾年等熱量資源充沛的年份,做好蘋果預(yù)旱防旱的措施;中性年降水較多的年份,謹(jǐn)防雨澇災(zāi)害對(duì)蘋果生產(chǎn)的影響。

    不同ENSO 年型下,極端溫度對(duì)蘋果產(chǎn)量影響也具有差異性,厄爾尼諾年,魯西、魯中地區(qū)的高溫?zé)岷εc蘋果減產(chǎn)率呈正相關(guān);魯中、魯西地區(qū)為高溫?zé)岷Ω甙l(fā)區(qū),高溫?zé)岷Πl(fā)生導(dǎo)致蘋果減產(chǎn)率上升。拉尼娜年和中性年,魯東、膠東半島地區(qū)極端低溫災(zāi)害與蘋果減產(chǎn)率呈正相關(guān),極端低溫災(zāi)害的發(fā)生,阻礙了蘋果產(chǎn)量的增產(chǎn)增收;煙臺(tái)地區(qū)與蘋果減產(chǎn)率呈負(fù)相關(guān),這可能由于煙臺(tái)地區(qū)種植技術(shù)、管理措施等因素降低了氣候條件對(duì)蘋果生產(chǎn)的影響。厄爾尼諾年,高溫?zé)岷μ鞌?shù)的增多導(dǎo)致蘋果減產(chǎn)風(fēng)險(xiǎn)增高。拉尼娜、中性年,高溫對(duì)山東蘋果減產(chǎn)影響較小。因此,厄爾尼諾、中性年應(yīng)謹(jǐn)防高溫?zé)岷?、低溫凍害?duì)蘋果產(chǎn)量的影響,拉尼娜年,謹(jǐn)防低溫冷害對(duì)蘋果減產(chǎn)的危害。

    但是,本文是基于站點(diǎn)氣象數(shù)據(jù)對(duì)山東蘋果旱澇災(zāi)害開(kāi)展研究,受數(shù)據(jù)年限及數(shù)據(jù)精細(xì)化的限制,并未充分考慮到蘋果需水量及蒸散量等因素,也未考慮蘋果品種、人工灌溉等因素的影響。在今后研究中,會(huì)收集相關(guān)數(shù)據(jù)進(jìn)一步開(kāi)展更深入的研究。

    4 結(jié)論

    本文通過(guò)山東1991-2019 年逐日氣象數(shù)據(jù)及山東蘋果地級(jí)市蘋果生產(chǎn)數(shù)據(jù),探究了不同ENSO 年型旱澇、極端溫度災(zāi)害的時(shí)空演變規(guī)律,分析不同ENSO 事件農(nóng)業(yè)氣象災(zāi)害對(duì)山東蘋果減產(chǎn)率的影響。過(guò)去29 年,厄爾尼諾年干旱發(fā)生頻率、影響程度和范圍較大;中性年雨澇災(zāi)害發(fā)生頻繁。從地區(qū)來(lái)看,魯西地區(qū)易發(fā)生干旱、魯南易發(fā)生雨澇災(zāi)害。膠東半島、魯東等地易發(fā)生低溫冷害、凍害;魯西、魯中及魯北等地高溫?zé)岷ο鄬?duì)嚴(yán)重。極端低溫災(zāi)害與厄爾尼諾南方濤動(dòng)指數(shù)呈負(fù)相關(guān),與拉尼娜年、中性年的關(guān)系反之。極端高溫與中性年南方濤動(dòng)指數(shù)呈負(fù)相關(guān)。對(duì)于蘋果產(chǎn)量,厄爾尼諾年干旱、高溫導(dǎo)致蘋果減產(chǎn)量增加,尤其在魯西地區(qū)。拉尼娜年、中性年低溫冷害使得蘋果減產(chǎn)率增加。總體而言,厄爾尼諾年蘋果栽培需預(yù)防干旱和高溫,中性年需預(yù)防雨澇、蘋果花期低溫和果實(shí)膨大期的高溫災(zāi)害,拉尼娜年需謹(jǐn)防低溫災(zāi)害的發(fā)生。

    猜你喜歡
    旱澇拉尼厄爾尼諾
    拉尼娜來(lái)了,這個(gè)“小女孩”不簡(jiǎn)單
    四川省1960-2017年旱澇時(shí)空變化特征分析
    厄爾尼諾現(xiàn)象橫行中國(guó)
    梅拉尼婭:世界上最受欺凌的人是我
    最小的一個(gè)
    厄爾尼諾現(xiàn)象及其影響
    氣候變化背景下西北干旱區(qū)旱澇的變化規(guī)律
    三步圖解厄爾尼諾與拉尼娜
    北方農(nóng)牧交錯(cuò)帶氣候變化與旱澇響應(yīng)特征
    寧夏旱澇災(zāi)害對(duì)氣候變化的響應(yīng)
    老女人水多毛片| 免费观看性生交大片5| 禁无遮挡网站| 久久久久久久久久久丰满| 不卡视频在线观看欧美| 2022亚洲国产成人精品| 欧美性猛交╳xxx乱大交人| 亚洲国产精品成人久久小说| 国产色爽女视频免费观看| 大香蕉久久网| 街头女战士在线观看网站| 男女边吃奶边做爰视频| 一级黄片播放器| 国产探花在线观看一区二区| 久久热精品热| 国产一区二区亚洲精品在线观看| 在线免费十八禁| 国产色婷婷99| av网站免费在线观看视频| av免费在线看不卡| 男人和女人高潮做爰伦理| 国产在线一区二区三区精| 黄色配什么色好看| 3wmmmm亚洲av在线观看| 日韩大片免费观看网站| av国产精品久久久久影院| av国产精品久久久久影院| 一本一本综合久久| 美女被艹到高潮喷水动态| 大陆偷拍与自拍| 最新中文字幕久久久久| 久久97久久精品| 韩国高清视频一区二区三区| 国产探花极品一区二区| 国产精品无大码| 激情 狠狠 欧美| 亚洲人与动物交配视频| 中文资源天堂在线| 黄片无遮挡物在线观看| 亚洲精品久久午夜乱码| 亚洲国产成人一精品久久久| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 黑人高潮一二区| 97热精品久久久久久| 中文字幕免费在线视频6| 人妻一区二区av| 午夜福利在线在线| 国产av国产精品国产| 韩国av在线不卡| 国产淫语在线视频| 亚洲av中文av极速乱| 丰满乱子伦码专区| 又爽又黄无遮挡网站| 欧美激情国产日韩精品一区| 欧美bdsm另类| 久久午夜福利片| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看 | 三级国产精品欧美在线观看| 岛国毛片在线播放| 九九爱精品视频在线观看| 蜜桃亚洲精品一区二区三区| 啦啦啦啦在线视频资源| 成人二区视频| 三级国产精品片| 国产欧美日韩一区二区三区在线 | 久久久欧美国产精品| 国产成人精品福利久久| 国产有黄有色有爽视频| 国产精品成人在线| 国产成人免费观看mmmm| 欧美日韩视频高清一区二区三区二| 日韩免费高清中文字幕av| 最近的中文字幕免费完整| 99久久精品一区二区三区| 欧美精品一区二区大全| 久久人人爽av亚洲精品天堂 | 一个人看视频在线观看www免费| 久久国产乱子免费精品| 一区二区三区四区激情视频| 国产乱人偷精品视频| 熟女电影av网| 久久久久久久午夜电影| 亚洲av不卡在线观看| 国产精品一区二区三区四区免费观看| 成人国产av品久久久| 精品国产一区二区三区久久久樱花 | 大又大粗又爽又黄少妇毛片口| 国产高潮美女av| 国产亚洲91精品色在线| av在线播放精品| 能在线免费看毛片的网站| 亚洲国产精品成人久久小说| 99热6这里只有精品| 国产爽快片一区二区三区| 亚洲天堂国产精品一区在线| 国产伦精品一区二区三区四那| 白带黄色成豆腐渣| 天天躁夜夜躁狠狠久久av| 草草在线视频免费看| 一级毛片电影观看| 亚洲va在线va天堂va国产| 超碰97精品在线观看| 午夜精品一区二区三区免费看| 女人久久www免费人成看片| 亚洲激情五月婷婷啪啪| 成人国产av品久久久| 人妻制服诱惑在线中文字幕| 精品久久久久久电影网| 麻豆成人午夜福利视频| 一个人看视频在线观看www免费| 又大又黄又爽视频免费| 色婷婷久久久亚洲欧美| 亚洲性久久影院| 精华霜和精华液先用哪个| 色播亚洲综合网| 久久精品国产亚洲av天美| 欧美日韩一区二区视频在线观看视频在线 | 一二三四中文在线观看免费高清| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 嫩草影院精品99| 男人狂女人下面高潮的视频| 一级毛片黄色毛片免费观看视频| 亚洲av电影在线观看一区二区三区 | 成年女人看的毛片在线观看| 久久久久久久午夜电影| 亚洲图色成人| 亚洲精品亚洲一区二区| 久久97久久精品| 中文在线观看免费www的网站| 亚洲欧美一区二区三区黑人 | 永久网站在线| 日日撸夜夜添| 午夜福利视频精品| 亚洲三级黄色毛片| 欧美日韩视频精品一区| 五月天丁香电影| 欧美老熟妇乱子伦牲交| 99久久精品热视频| 久久久久精品性色| 国产成人精品福利久久| 亚洲国产欧美人成| 亚洲人成网站在线播| 精品熟女少妇av免费看| 91久久精品国产一区二区三区| 欧美丝袜亚洲另类| 网址你懂的国产日韩在线| 最近最新中文字幕免费大全7| 黄色视频在线播放观看不卡| 免费大片黄手机在线观看| 国产成人91sexporn| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 亚洲色图av天堂| 我的女老师完整版在线观看| 久久人人爽人人片av| 人妻制服诱惑在线中文字幕| 一级a做视频免费观看| 成年免费大片在线观看| 欧美三级亚洲精品| 嫩草影院新地址| 中文字幕人妻熟人妻熟丝袜美| 我的老师免费观看完整版| 日韩强制内射视频| 精品人妻一区二区三区麻豆| 婷婷色麻豆天堂久久| 亚洲国产精品成人综合色| av在线老鸭窝| 男女啪啪激烈高潮av片| 性插视频无遮挡在线免费观看| 天堂网av新在线| 国内少妇人妻偷人精品xxx网站| 国产精品女同一区二区软件| 老师上课跳d突然被开到最大视频| 两个人的视频大全免费| 69av精品久久久久久| 99热这里只有是精品50| 国产极品天堂在线| 久久午夜福利片| 亚洲av中文字字幕乱码综合| 国产精品一及| 欧美丝袜亚洲另类| 欧美最新免费一区二区三区| 久久久国产一区二区| 国产精品伦人一区二区| 国产高潮美女av| 亚洲色图综合在线观看| 久久精品熟女亚洲av麻豆精品| 国产精品福利在线免费观看| 国产在视频线精品| 久久精品国产亚洲av涩爱| av免费在线看不卡| 亚洲精品456在线播放app| 国产精品久久久久久久电影| av国产精品久久久久影院| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久一区二区三区| av卡一久久| 免费人成在线观看视频色| 亚洲成人精品中文字幕电影| 中文乱码字字幕精品一区二区三区| 国产毛片在线视频| xxx大片免费视频| 青春草国产在线视频| 九九在线视频观看精品| 久久精品久久精品一区二区三区| 亚洲av中文字字幕乱码综合| 日韩三级伦理在线观看| 国产精品伦人一区二区| 亚洲av中文av极速乱| 亚洲国产最新在线播放| 天天一区二区日本电影三级| 男女国产视频网站| 欧美一级a爱片免费观看看| 亚洲天堂av无毛| 日韩亚洲欧美综合| 一个人看的www免费观看视频| 亚洲内射少妇av| 99热这里只有精品一区| 亚洲成人精品中文字幕电影| 狂野欧美白嫩少妇大欣赏| 老司机影院毛片| 亚洲av电影在线观看一区二区三区 | 精品人妻视频免费看| 最近手机中文字幕大全| 黄色一级大片看看| 我要看日韩黄色一级片| 成人欧美大片| 免费播放大片免费观看视频在线观看| 99热6这里只有精品| 九九在线视频观看精品| 国产免费视频播放在线视频| 特级一级黄色大片| 欧美国产精品一级二级三级 | 极品教师在线视频| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 国产黄片美女视频| 青春草亚洲视频在线观看| 久久久a久久爽久久v久久| 国产精品国产三级专区第一集| 久久午夜福利片| 色5月婷婷丁香| av又黄又爽大尺度在线免费看| 又爽又黄无遮挡网站| 国产女主播在线喷水免费视频网站| 国产黄片视频在线免费观看| 可以在线观看毛片的网站| 男女下面进入的视频免费午夜| 下体分泌物呈黄色| 精品熟女少妇av免费看| 亚洲aⅴ乱码一区二区在线播放| 成人午夜精彩视频在线观看| 午夜激情久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 91aial.com中文字幕在线观看| 亚洲四区av| 看免费成人av毛片| 午夜亚洲福利在线播放| 免费av观看视频| 日韩大片免费观看网站| 亚洲美女搞黄在线观看| 久久国产乱子免费精品| 成年免费大片在线观看| 久久久久久九九精品二区国产| 国产精品久久久久久av不卡| 免费看不卡的av| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线 | 高清在线视频一区二区三区| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| freevideosex欧美| 日韩 亚洲 欧美在线| 人人妻人人澡人人爽人人夜夜| 激情 狠狠 欧美| 欧美97在线视频| 国产精品av视频在线免费观看| 99热6这里只有精品| 身体一侧抽搐| 一区二区三区免费毛片| 五月开心婷婷网| 国产精品精品国产色婷婷| 亚洲人成网站高清观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产高潮美女av| 一级毛片我不卡| 在线观看一区二区三区激情| 只有这里有精品99| av女优亚洲男人天堂| kizo精华| 精品国产乱码久久久久久小说| 国产乱来视频区| av在线蜜桃| 18禁在线无遮挡免费观看视频| av在线亚洲专区| 色哟哟·www| 国产成人精品一,二区| 欧美bdsm另类| 国产欧美日韩一区二区三区在线 | 最近最新中文字幕大全电影3| 久久精品久久久久久噜噜老黄| 色视频www国产| 久久热精品热| 亚洲人与动物交配视频| 国产乱人视频| 狂野欧美白嫩少妇大欣赏| 国内精品美女久久久久久| 毛片一级片免费看久久久久| 少妇高潮的动态图| 国产真实伦视频高清在线观看| 五月天丁香电影| 三级男女做爰猛烈吃奶摸视频| 久久6这里有精品| 久久亚洲国产成人精品v| 夫妻午夜视频| 国产伦精品一区二区三区四那| 在线观看三级黄色| 中文字幕人妻熟人妻熟丝袜美| 免费播放大片免费观看视频在线观看| 99热这里只有精品一区| 少妇人妻 视频| 国产精品久久久久久精品古装| 天天躁夜夜躁狠狠久久av| 久久国内精品自在自线图片| 国产老妇女一区| 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看 | 禁无遮挡网站| 身体一侧抽搐| 在线免费十八禁| 国内精品美女久久久久久| 春色校园在线视频观看| 老司机影院毛片| 天天躁日日操中文字幕| 国产成人freesex在线| 舔av片在线| 日韩欧美一区视频在线观看 | 三级经典国产精品| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 国产精品一区二区三区四区免费观看| 亚洲国产欧美人成| 国产高清三级在线| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类| 制服丝袜香蕉在线| 午夜亚洲福利在线播放| 成人综合一区亚洲| 亚洲欧美成人综合另类久久久| 久久久久网色| 黄色日韩在线| 听说在线观看完整版免费高清| 青春草视频在线免费观看| 草草在线视频免费看| 成年女人看的毛片在线观看| 日本午夜av视频| 简卡轻食公司| av.在线天堂| 免费大片18禁| 免费黄频网站在线观看国产| 亚洲成人中文字幕在线播放| 国产精品99久久99久久久不卡 | 美女脱内裤让男人舔精品视频| 亚洲av电影在线观看一区二区三区 | 不卡视频在线观看欧美| 一级毛片电影观看| 黄色欧美视频在线观看| 97人妻精品一区二区三区麻豆| 国产精品99久久99久久久不卡 | av天堂中文字幕网| 日本猛色少妇xxxxx猛交久久| 国产精品蜜桃在线观看| 1000部很黄的大片| 偷拍熟女少妇极品色| 国产91av在线免费观看| 日韩人妻高清精品专区| 免费观看无遮挡的男女| 国产高清三级在线| 亚洲欧美精品专区久久| 欧美xxxx黑人xx丫x性爽| 亚洲av中文字字幕乱码综合| 亚洲欧美成人综合另类久久久| 最近2019中文字幕mv第一页| 色综合色国产| 久久久久久国产a免费观看| 国产欧美日韩一区二区三区在线 | 有码 亚洲区| 国产精品蜜桃在线观看| 男的添女的下面高潮视频| 国精品久久久久久国模美| 久热这里只有精品99| 精品久久久久久电影网| 国产久久久一区二区三区| 一区二区av电影网| 国产一区有黄有色的免费视频| 丰满人妻一区二区三区视频av| av在线老鸭窝| 亚洲国产欧美人成| 国产成人一区二区在线| 亚洲国产精品成人综合色| 波野结衣二区三区在线| 大香蕉97超碰在线| 日韩av不卡免费在线播放| 午夜激情福利司机影院| 97人妻精品一区二区三区麻豆| 中文字幕制服av| 在线观看一区二区三区| 亚洲国产精品国产精品| 五月玫瑰六月丁香| 国产精品国产av在线观看| 99九九线精品视频在线观看视频| 搞女人的毛片| 色播亚洲综合网| 欧美日韩综合久久久久久| 美女cb高潮喷水在线观看| 我的老师免费观看完整版| 亚洲欧美一区二区三区国产| 日韩欧美精品v在线| 欧美极品一区二区三区四区| 99热网站在线观看| 免费少妇av软件| 国产一区二区亚洲精品在线观看| 亚洲性久久影院| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 一区二区三区乱码不卡18| 日本-黄色视频高清免费观看| 在线 av 中文字幕| 日日啪夜夜撸| 18禁动态无遮挡网站| 少妇人妻久久综合中文| 夜夜爽夜夜爽视频| 欧美性猛交╳xxx乱大交人| 国产一区二区三区综合在线观看 | 观看免费一级毛片| 一级毛片电影观看| 天天一区二区日本电影三级| 精品久久久精品久久久| av.在线天堂| 亚洲精品国产av成人精品| 国产精品一区二区三区四区免费观看| 亚洲精品一二三| 亚州av有码| 高清毛片免费看| 国产精品无大码| 亚洲精品第二区| 制服丝袜香蕉在线| 亚洲精品国产成人久久av| 欧美zozozo另类| 伦精品一区二区三区| 国产av不卡久久| 国产大屁股一区二区在线视频| 久久精品久久久久久噜噜老黄| 99热这里只有是精品在线观看| 成人毛片60女人毛片免费| 欧美日韩国产mv在线观看视频 | 国产又色又爽无遮挡免| 国产亚洲91精品色在线| 国产日韩欧美亚洲二区| 美女主播在线视频| 国产69精品久久久久777片| 久久精品人妻少妇| 毛片女人毛片| 亚洲电影在线观看av| 亚洲国产最新在线播放| 看十八女毛片水多多多| 一区二区三区免费毛片| 在线观看av片永久免费下载| 18禁动态无遮挡网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 哪个播放器可以免费观看大片| 亚洲国产精品国产精品| 国国产精品蜜臀av免费| 人妻制服诱惑在线中文字幕| 亚洲最大成人手机在线| av播播在线观看一区| 直男gayav资源| 大片电影免费在线观看免费| 亚洲在久久综合| 性色avwww在线观看| 久久精品人妻少妇| 熟妇人妻不卡中文字幕| 九色成人免费人妻av| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 日韩一本色道免费dvd| 国产v大片淫在线免费观看| 日韩大片免费观看网站| 熟女av电影| 国产免费一区二区三区四区乱码| 天天躁夜夜躁狠狠久久av| 极品教师在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 成人一区二区视频在线观看| 久久99热这里只频精品6学生| 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 免费播放大片免费观看视频在线观看| 男人添女人高潮全过程视频| 最近中文字幕高清免费大全6| 肉色欧美久久久久久久蜜桃 | 精品一区二区免费观看| 亚洲人成网站在线观看播放| 一二三四中文在线观看免费高清| 日韩欧美 国产精品| 极品教师在线视频| 又粗又硬又长又爽又黄的视频| 亚洲,欧美,日韩| 国模一区二区三区四区视频| 男女边摸边吃奶| 中文字幕制服av| 又粗又硬又长又爽又黄的视频| 特级一级黄色大片| 国产亚洲精品久久久com| 国产精品成人在线| 又黄又爽又刺激的免费视频.| 搡老乐熟女国产| 三级男女做爰猛烈吃奶摸视频| 欧美xxxx性猛交bbbb| 少妇人妻 视频| a级一级毛片免费在线观看| 亚洲国产精品成人综合色| 亚洲av欧美aⅴ国产| 美女视频免费永久观看网站| 爱豆传媒免费全集在线观看| 99re6热这里在线精品视频| 国产亚洲av嫩草精品影院| 国产午夜精品久久久久久一区二区三区| 老师上课跳d突然被开到最大视频| 搡老乐熟女国产| 99精国产麻豆久久婷婷| 建设人人有责人人尽责人人享有的 | 国产永久视频网站| 国产欧美日韩一区二区三区在线 | 99久久九九国产精品国产免费| 国产精品久久久久久精品电影小说 | 日韩制服骚丝袜av| 又爽又黄a免费视频| 美女国产视频在线观看| 国产精品国产av在线观看| 好男人视频免费观看在线| 精品熟女少妇av免费看| 亚洲成人久久爱视频| 亚洲欧美精品专区久久| 精品久久久噜噜| 国产av国产精品国产| 成人毛片a级毛片在线播放| 特大巨黑吊av在线直播| 色哟哟·www| 波野结衣二区三区在线| 国产精品国产三级国产专区5o| 色吧在线观看| 在线精品无人区一区二区三 | 中文字幕久久专区| 国精品久久久久久国模美| 午夜福利高清视频| 国产乱人偷精品视频| 美女xxoo啪啪120秒动态图| 免费观看在线日韩| 中文字幕亚洲精品专区| 久久久精品94久久精品| 99视频精品全部免费 在线| 美女高潮的动态| 中文精品一卡2卡3卡4更新| 岛国毛片在线播放| 97精品久久久久久久久久精品| 亚洲av成人精品一二三区| 97超碰精品成人国产| 国产亚洲av嫩草精品影院| 婷婷色综合大香蕉| 精品久久久久久久人妻蜜臀av| 天天躁日日操中文字幕| 免费高清在线观看视频在线观看| 国产大屁股一区二区在线视频| 日韩成人av中文字幕在线观看| 国产高潮美女av| 国产乱人偷精品视频| 亚洲av一区综合| 午夜福利网站1000一区二区三区| 亚洲美女搞黄在线观看| 99九九线精品视频在线观看视频| 日本黄大片高清| 九九爱精品视频在线观看| 青春草亚洲视频在线观看| 日本黄大片高清| 九九爱精品视频在线观看| 99九九线精品视频在线观看视频| 大话2 男鬼变身卡| 亚洲精品乱码久久久v下载方式| 欧美性猛交╳xxx乱大交人| 亚洲欧美一区二区三区黑人 | 18禁动态无遮挡网站| a级毛片免费高清观看在线播放| 国产老妇女一区| 国产v大片淫在线免费观看| 欧美另类一区| 国产熟女欧美一区二区| 久久久久久久大尺度免费视频| 激情五月婷婷亚洲| 欧美三级亚洲精品| 蜜桃久久精品国产亚洲av| 18禁在线播放成人免费| 边亲边吃奶的免费视频| av黄色大香蕉| 欧美日韩视频高清一区二区三区二| 边亲边吃奶的免费视频| 极品少妇高潮喷水抽搐| 亚洲色图综合在线观看| 日日啪夜夜撸| 国产探花在线观看一区二区| 欧美成人一区二区免费高清观看| 日本色播在线视频|