• <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)
    中文字幕人妻丝袜制服| 亚洲欧美清纯卡通| 中文字幕亚洲精品专区| 国产伦理片在线播放av一区| 亚洲精品日韩在线中文字幕| 男女高潮啪啪啪动态图| 亚洲精品日本国产第一区| 国产免费视频播放在线视频| 国产高清不卡午夜福利| 人妻 亚洲 视频| 久久毛片免费看一区二区三区| 亚洲人成网站在线观看播放| 亚洲av中文av极速乱| 亚洲国产欧美在线一区| 美女国产高潮福利片在线看| 成人午夜精彩视频在线观看| 国产黄色免费在线视频| 亚洲情色 制服丝袜| 欧美xxxx性猛交bbbb| 精品久久久噜噜| av.在线天堂| 99热网站在线观看| 国产亚洲精品久久久com| 乱码一卡2卡4卡精品| 亚洲成色77777| 国产一区二区三区av在线| 亚洲欧洲精品一区二区精品久久久 | √禁漫天堂资源中文www| 插阴视频在线观看视频| 久久国产精品男人的天堂亚洲 | 久久久久久伊人网av| 久久国产亚洲av麻豆专区| 中文字幕亚洲精品专区| 久久国产精品大桥未久av| 精品国产乱码久久久久久小说| 观看美女的网站| 久久精品久久久久久久性| 国产精品人妻久久久久久| 国产精品 国内视频| 国产亚洲欧美精品永久| 最近的中文字幕免费完整| 国产淫语在线视频| 男女边摸边吃奶| 男女边吃奶边做爰视频| 久久久精品区二区三区| 久久久久久久久久成人| 99热国产这里只有精品6| 久久97久久精品| 久热久热在线精品观看| 黄色一级大片看看| 我的女老师完整版在线观看| 男人爽女人下面视频在线观看| 亚洲欧美色中文字幕在线| 高清午夜精品一区二区三区| 国产成人免费观看mmmm| 亚洲欧美清纯卡通| 考比视频在线观看| 国产精品一国产av| 免费观看性生交大片5| 亚洲精品456在线播放app| 国产精品99久久久久久久久| 99热国产这里只有精品6| 日本爱情动作片www.在线观看| 水蜜桃什么品种好| 涩涩av久久男人的天堂| 满18在线观看网站| 精品久久久久久电影网| 久久韩国三级中文字幕| 老女人水多毛片| 99国产综合亚洲精品| 色婷婷av一区二区三区视频| 国产老妇伦熟女老妇高清| 精品亚洲成a人片在线观看| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 美女国产高潮福利片在线看| 午夜福利,免费看| 成人手机av| 久久精品久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 99精国产麻豆久久婷婷| 一级a做视频免费观看| 亚洲欧美一区二区三区国产| 国产在线视频一区二区| 国产一区二区在线观看av| 美女内射精品一级片tv| 亚洲av成人精品一区久久| 日本黄大片高清| 亚洲内射少妇av| 亚洲欧美一区二区三区国产| 99re6热这里在线精品视频| 国产 一区精品| av在线播放精品| 久久久久久人妻| 成人亚洲欧美一区二区av| 日韩三级伦理在线观看| 97超视频在线观看视频| 狂野欧美激情性bbbbbb| 熟女av电影| 91aial.com中文字幕在线观看| 国产综合精华液| 丝袜在线中文字幕| 黑丝袜美女国产一区| 成人国语在线视频| 少妇熟女欧美另类| 久久女婷五月综合色啪小说| 久久久国产一区二区| 日韩 亚洲 欧美在线| 熟女av电影| 少妇被粗大猛烈的视频| 亚洲四区av| 在线观看免费视频网站a站| 国产女主播在线喷水免费视频网站| 亚洲内射少妇av| 一边摸一边做爽爽视频免费| 亚洲在久久综合| av福利片在线| av女优亚洲男人天堂| 久久久欧美国产精品| 国产精品一国产av| 亚洲四区av| 欧美变态另类bdsm刘玥| 精品国产乱码久久久久久小说| 成人无遮挡网站| 免费少妇av软件| 亚洲精品视频女| 校园人妻丝袜中文字幕| 国产亚洲av片在线观看秒播厂| 我的老师免费观看完整版| 免费av中文字幕在线| 性色av一级| 黄色一级大片看看| 色哟哟·www| 国产精品一二三区在线看| 国产欧美亚洲国产| 啦啦啦在线观看免费高清www| 一本大道久久a久久精品| 99热这里只有精品一区| 熟女av电影| 国产成人精品在线电影| 亚洲av国产av综合av卡| 十分钟在线观看高清视频www| 9色porny在线观看| 99久久综合免费| 国产精品国产三级国产专区5o| 91成人精品电影| 久久精品国产亚洲av涩爱| 纵有疾风起免费观看全集完整版| 看非洲黑人一级黄片| 男女国产视频网站| 性色avwww在线观看| 国产69精品久久久久777片| 在线观看一区二区三区激情| 久久久国产精品麻豆| 成人手机av| 欧美少妇被猛烈插入视频| 国产亚洲精品久久久com| 精品人妻一区二区三区麻豆| 夫妻性生交免费视频一级片| 欧美国产精品一级二级三级| 嘟嘟电影网在线观看| 老司机影院毛片| 久久久国产一区二区| 国产无遮挡羞羞视频在线观看| 久久女婷五月综合色啪小说| 久久影院123| 亚洲人成网站在线播| 国产精品久久久久久久电影| 人人妻人人添人人爽欧美一区卜| 国产日韩欧美在线精品| 欧美日韩精品成人综合77777| 一区二区三区乱码不卡18| 考比视频在线观看| 最近最新中文字幕免费大全7| 水蜜桃什么品种好| 精品国产一区二区三区久久久樱花| 午夜久久久在线观看| 国产一级毛片在线| 边亲边吃奶的免费视频| 午夜福利视频精品| 国产精品不卡视频一区二区| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放| 中文天堂在线官网| 18禁在线无遮挡免费观看视频| 日韩一区二区视频免费看| 欧美日韩在线观看h| 国产免费一区二区三区四区乱码| 色5月婷婷丁香| 美女国产视频在线观看| a级毛色黄片| 久久久久久久久久人人人人人人| 久久影院123| 国产日韩欧美亚洲二区| 黄片无遮挡物在线观看| 国产亚洲精品久久久com| 日日摸夜夜添夜夜添av毛片| 精品国产一区二区久久| 一级黄片播放器| 亚洲经典国产精华液单| 免费大片18禁| 免费人成在线观看视频色| av免费观看日本| 色婷婷久久久亚洲欧美| 少妇高潮的动态图| 欧美性感艳星| 97在线人人人人妻| 哪个播放器可以免费观看大片| 中文字幕精品免费在线观看视频 | 蜜臀久久99精品久久宅男| 18禁在线无遮挡免费观看视频| 久久久欧美国产精品| 极品人妻少妇av视频| 一级黄片播放器| 免费高清在线观看日韩| 午夜免费鲁丝| 国产午夜精品一二区理论片| 97超视频在线观看视频| 超碰97精品在线观看| 久久人妻熟女aⅴ| 天堂8中文在线网| 国产日韩欧美在线精品| 国产伦精品一区二区三区视频9| 九色亚洲精品在线播放| 天堂中文最新版在线下载| 观看美女的网站| 夜夜爽夜夜爽视频| 少妇人妻久久综合中文| 欧美亚洲日本最大视频资源| 日韩av在线免费看完整版不卡| 丝袜美足系列| 高清视频免费观看一区二区| 日本猛色少妇xxxxx猛交久久| 熟妇人妻不卡中文字幕| 日韩av在线免费看完整版不卡| 国模一区二区三区四区视频| 看免费成人av毛片| 亚洲国产精品一区三区| 国产日韩欧美在线精品| 久久久久久久久大av| 热re99久久精品国产66热6| 成人亚洲欧美一区二区av| 亚洲情色 制服丝袜| 久久国产精品男人的天堂亚洲 | 亚洲成人手机| 99国产精品免费福利视频| 亚洲精品第二区| 国产av码专区亚洲av| 日本av手机在线免费观看| 男女无遮挡免费网站观看| 搡女人真爽免费视频火全软件| 久久精品夜色国产| 成人毛片60女人毛片免费| 久久 成人 亚洲| 午夜视频国产福利| 国产男人的电影天堂91| 国产在线一区二区三区精| 在线观看免费日韩欧美大片 | 2022亚洲国产成人精品| 国产精品99久久99久久久不卡 | 国产精品国产三级国产专区5o| 天堂中文最新版在线下载| 久久99热这里只频精品6学生| 久热这里只有精品99| 精品一区二区免费观看| 成人二区视频| 成人手机av| 午夜激情av网站| 99视频精品全部免费 在线| 天堂8中文在线网| 久久婷婷青草| 国产欧美另类精品又又久久亚洲欧美| 18禁观看日本| 欧美bdsm另类| 晚上一个人看的免费电影| 亚洲精品自拍成人| 精品人妻熟女av久视频| 人妻夜夜爽99麻豆av| 国产精品久久久久久av不卡| 黄色一级大片看看| 制服人妻中文乱码| 在线观看免费视频网站a站| 免费播放大片免费观看视频在线观看| 超碰97精品在线观看| 日韩一本色道免费dvd| 99久久精品国产国产毛片| 国产在线免费精品| 久久青草综合色| 热99久久久久精品小说推荐| 一级毛片电影观看| 日韩一区二区视频免费看| 成人国产av品久久久| 欧美日韩视频精品一区| 韩国av在线不卡| 亚洲av日韩在线播放| 成人无遮挡网站| 黄色配什么色好看| 国产黄片视频在线免费观看| 亚洲,欧美,日韩| 男女国产视频网站| av专区在线播放| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| 蜜臀久久99精品久久宅男| 亚洲欧美中文字幕日韩二区| 三级国产精品欧美在线观看| 欧美 亚洲 国产 日韩一| 看十八女毛片水多多多| 视频中文字幕在线观看| 人妻一区二区av| 一级毛片电影观看| 久久毛片免费看一区二区三区| 久久人人爽人人爽人人片va| 国产淫语在线视频| 一级a做视频免费观看| 国产日韩欧美在线精品| 777米奇影视久久| 十八禁高潮呻吟视频| 看非洲黑人一级黄片| 国产无遮挡羞羞视频在线观看| 亚洲精品视频女| 久久久久精品性色| 伊人亚洲综合成人网| 春色校园在线视频观看| 久久人妻熟女aⅴ| 一区二区三区乱码不卡18| 亚洲av二区三区四区| 免费观看性生交大片5| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本黄大片高清| 一本久久精品| 免费观看av网站的网址| 热99国产精品久久久久久7| 久久久久久久久久人人人人人人| 最后的刺客免费高清国语| 看十八女毛片水多多多| 一区二区三区免费毛片| 亚洲成人av在线免费| 热re99久久国产66热| 精品酒店卫生间| 久久久久久久精品精品| 日本黄色日本黄色录像| av国产精品久久久久影院| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久免费观看电影| 一级黄片播放器| 日韩av在线免费看完整版不卡| 国产 精品1| 成人毛片60女人毛片免费| 国产精品熟女久久久久浪| 肉色欧美久久久久久久蜜桃| 婷婷色综合www| 久久国产精品大桥未久av| av卡一久久| 九草在线视频观看| 亚洲成人手机| 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| 少妇高潮的动态图| 国模一区二区三区四区视频| 欧美日韩一区二区视频在线观看视频在线| 免费观看的影片在线观看| 高清午夜精品一区二区三区| 在线 av 中文字幕| 成人毛片a级毛片在线播放| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 下体分泌物呈黄色| 天堂中文最新版在线下载| 一区二区三区四区激情视频| 麻豆乱淫一区二区| 成年av动漫网址| 亚洲经典国产精华液单| 亚洲精品国产色婷婷电影| 中文字幕人妻熟人妻熟丝袜美| av播播在线观看一区| 久久 成人 亚洲| 秋霞伦理黄片| 国产伦精品一区二区三区视频9| 丝袜在线中文字幕| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 全区人妻精品视频| 少妇猛男粗大的猛烈进出视频| 精品一品国产午夜福利视频| 免费人妻精品一区二区三区视频| 免费观看无遮挡的男女| 久热这里只有精品99| 草草在线视频免费看| 久久久久国产精品人妻一区二区| 国产精品一区二区在线观看99| 中国国产av一级| 中文字幕久久专区| 一区二区av电影网| 韩国av在线不卡| 寂寞人妻少妇视频99o| av国产精品久久久久影院| av视频免费观看在线观看| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 一级,二级,三级黄色视频| 亚洲综合精品二区| 欧美日本中文国产一区发布| 九色成人免费人妻av| 亚洲美女搞黄在线观看| 免费人妻精品一区二区三区视频| 国内精品宾馆在线| 看免费成人av毛片| 成人影院久久| 欧美一级a爱片免费观看看| www.色视频.com| 欧美 日韩 精品 国产| 街头女战士在线观看网站| 最近2019中文字幕mv第一页| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 一级毛片aaaaaa免费看小| 夫妻午夜视频| 国产成人一区二区在线| 国产成人a∨麻豆精品| 精品酒店卫生间| 丰满少妇做爰视频| 天堂中文最新版在线下载| 久久久久精品性色| 91国产中文字幕| 纯流量卡能插随身wifi吗| 日韩不卡一区二区三区视频在线| 国产亚洲精品第一综合不卡 | xxxhd国产人妻xxx| 成人亚洲欧美一区二区av| 大香蕉97超碰在线| 日韩一本色道免费dvd| 久久精品久久久久久久性| 看非洲黑人一级黄片| 高清在线视频一区二区三区| 一边亲一边摸免费视频| a级片在线免费高清观看视频| 午夜免费鲁丝| 99视频精品全部免费 在线| 久久久久网色| 大香蕉久久成人网| 久久精品国产a三级三级三级| 国产精品国产三级国产av玫瑰| 在线观看免费视频网站a站| 男女高潮啪啪啪动态图| 久久久久久久久久成人| 热re99久久国产66热| 欧美xxⅹ黑人| 在线免费观看不下载黄p国产| 最近手机中文字幕大全| 日韩成人伦理影院| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人 | 日韩亚洲欧美综合| 交换朋友夫妻互换小说| 国内精品宾馆在线| 高清午夜精品一区二区三区| 91久久精品电影网| 人人妻人人添人人爽欧美一区卜| 成年女人在线观看亚洲视频| 国产欧美另类精品又又久久亚洲欧美| 大香蕉97超碰在线| 日韩熟女老妇一区二区性免费视频| 黑人欧美特级aaaaaa片| 日韩视频在线欧美| 一级毛片 在线播放| 成年av动漫网址| 久久99蜜桃精品久久| 国产精品久久久久久精品电影小说| 欧美日韩综合久久久久久| 高清不卡的av网站| 国产毛片在线视频| 国产片内射在线| 在线观看免费视频网站a站| av免费在线看不卡| 三上悠亚av全集在线观看| 交换朋友夫妻互换小说| av不卡在线播放| 国产av精品麻豆| 男男h啪啪无遮挡| 亚洲在久久综合| 精品一区在线观看国产| 国产免费又黄又爽又色| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻精品综合一区二区| 欧美日韩成人在线一区二区| 欧美成人午夜免费资源| 亚洲内射少妇av| 国产成人freesex在线| 日韩欧美精品免费久久| 18在线观看网站| 成年av动漫网址| 一级片'在线观看视频| 中国美白少妇内射xxxbb| 久久午夜福利片| 美女大奶头黄色视频| 婷婷色av中文字幕| 一本色道久久久久久精品综合| 久久久精品免费免费高清| 国产av国产精品国产| 国产成人精品久久久久久| 天堂俺去俺来也www色官网| 国产伦精品一区二区三区视频9| 成年女人在线观看亚洲视频| 亚洲精品av麻豆狂野| 精品熟女少妇av免费看| 韩国高清视频一区二区三区| 亚洲av.av天堂| 亚洲人成网站在线观看播放| 99久久综合免费| 久久免费观看电影| 亚洲精品乱码久久久v下载方式| 午夜久久久在线观看| 亚洲国产日韩一区二区| 在线 av 中文字幕| 最黄视频免费看| 欧美一级a爱片免费观看看| 观看美女的网站| 日本-黄色视频高清免费观看| 国产69精品久久久久777片| 色吧在线观看| 伦精品一区二区三区| 高清欧美精品videossex| 国产精品偷伦视频观看了| videos熟女内射| 中国国产av一级| 99视频精品全部免费 在线| 人妻系列 视频| 精品久久久久久电影网| 精品人妻在线不人妻| 美女福利国产在线| 母亲3免费完整高清在线观看 | 国产国拍精品亚洲av在线观看| 伦精品一区二区三区| 久久国内精品自在自线图片| 亚洲熟女精品中文字幕| 中文字幕久久专区| 啦啦啦视频在线资源免费观看| 男女免费视频国产| 欧美日韩亚洲高清精品| 啦啦啦在线观看免费高清www| 国产一区二区三区综合在线观看 | 男人添女人高潮全过程视频| 亚洲欧美日韩卡通动漫| 国产永久视频网站| 亚洲av在线观看美女高潮| 久久久精品94久久精品| 免费高清在线观看日韩| 99热这里只有是精品在线观看| 十分钟在线观看高清视频www| 国产毛片在线视频| 一本一本综合久久| 亚洲精品日韩av片在线观看| 精品视频人人做人人爽| 嫩草影院入口| 国产黄色免费在线视频| 午夜福利视频在线观看免费| 久久人人爽av亚洲精品天堂| 欧美三级亚洲精品| 欧美 日韩 精品 国产| 精品久久国产蜜桃| 人人妻人人澡人人爽人人夜夜| 秋霞在线观看毛片| 七月丁香在线播放| 91精品三级在线观看| 男人爽女人下面视频在线观看| 国产极品粉嫩免费观看在线 | 国产一区亚洲一区在线观看| 日韩,欧美,国产一区二区三区| 精品一区二区三区视频在线| 亚洲美女视频黄频| 中国美白少妇内射xxxbb| 2021少妇久久久久久久久久久| 美女主播在线视频| 久久久精品免费免费高清| 一区二区三区免费毛片| 91国产中文字幕| 国产色婷婷99| 丝袜美足系列| 中文欧美无线码| 日韩电影二区| 美女内射精品一级片tv| 久久99蜜桃精品久久| 亚洲三级黄色毛片| 亚洲高清免费不卡视频| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美在线精品| 午夜福利视频精品| 亚洲精品,欧美精品| 精品一区二区三卡| 精品卡一卡二卡四卡免费| 国产精品偷伦视频观看了| 日韩视频在线欧美| 亚洲婷婷狠狠爱综合网| 黄片无遮挡物在线观看| 精品久久久噜噜| 大片免费播放器 马上看| 中文欧美无线码| 亚洲丝袜综合中文字幕| 最后的刺客免费高清国语| 国产精品一国产av| 久久久精品免费免费高清| 日韩免费高清中文字幕av| 少妇被粗大猛烈的视频| 91精品一卡2卡3卡4卡| 国产精品嫩草影院av在线观看| 亚洲av成人精品一二三区| 国产精品不卡视频一区二区| 女性被躁到高潮视频| 精品一区二区免费观看| 一个人看视频在线观看www免费| 最近2019中文字幕mv第一页| 一本色道久久久久久精品综合| 精品酒店卫生间| 建设人人有责人人尽责人人享有的|