郭曉雷 申雙和 張磊

摘要：利用寧夏枸杞種植區(qū)15個(gè)氣象站1961—2017年最低氣溫觀測(cè)資料，結(jié)合寧夏枸杞春霜凍致災(zāi)指標(biāo)，采用線性傾向估計(jì)、Moflet小波分析、正交經(jīng)驗(yàn)函數(shù)（EOF）分析等方法，研究寧夏枸杞種植區(qū)春霜凍發(fā)生的時(shí)空分布特征。結(jié)果表明，（1）寧夏枸杞種植區(qū)春霜凍發(fā)生次數(shù)呈現(xiàn)南北多、中間少的分布特征，晚霜凍日也呈現(xiàn)南部偏晚、中北部偏早的空間分布特征，霜凍發(fā)生次數(shù)最多和春霜凍日結(jié)束最晚都集中在南部地區(qū)，尤以興仁最為嚴(yán)重;（2）1961—2017年寧夏枸杞終霜凍日以2.1 d/10年的速率提前，且階段性變化特征明顯，20世紀(jì)70年代春霜凍都在5月上旬后結(jié)束，而到20世紀(jì)90年代末之后，大部分年份春霜凍都在5月上旬前結(jié)束;57年尺度上枸杞春霜凍發(fā)生次數(shù)以 2.1次/10年的速率減少，在1992年之前發(fā)生次數(shù)偏多，之后發(fā)生明顯突變，且1992年前存在明顯的7年、準(zhǔn)4年的周期變化，之后周期性變化不明顯。（3）寧夏枸杞春霜凍發(fā)生次數(shù)在大部分年份都呈現(xiàn)偏多或偏少的一致性特征，也存在明顯的南北差異分布，但總體一致性是寧夏枸杞春霜凍發(fā)生次數(shù)變化的主要特征。

寧夏枸杞馳名中外，是唯一被載入《中華人民共和國藥典》（2010年版）的品種，因其較高的營養(yǎng)價(jià)值和藥用保健功效受到廣大消費(fèi)者青睞。近年來，隨著政府的大力扶持，寧夏已形成了以中寧為核心、清水河流域和賀蘭山東麓為兩翼的枸杞產(chǎn)業(yè)帶，枸杞產(chǎn)業(yè)已經(jīng)成為寧夏農(nóng)業(yè)特色優(yōu)勢(shì)產(chǎn)業(yè)的重要組成部分。霜凍是威脅寧夏枸杞產(chǎn)業(yè)發(fā)展的主要?dú)庀鬄?zāi)害之一，且寧夏是霜凍多發(fā)區(qū)，正常年份春霜凍重于秋霜凍[1-2]，由于春季正值枸杞嫩梢期至花果期，春季霜凍災(zāi)害的發(fā)生直接影響枸杞頭茬果和夏果的產(chǎn)量及品質(zhì)[3]，而枸杞頭茬果和夏果的營養(yǎng)成分均高于秋果[4-5]，經(jīng)濟(jì)價(jià)值更高，春季霜凍造成的損失尤為嚴(yán)重。因此，摸清寧夏枸杞種植區(qū)春季霜凍發(fā)生規(guī)律，可為寧夏枸杞春季霜凍預(yù)報(bào)提供基礎(chǔ)理論依據(jù)，從而提高防御枸杞春霜凍災(zāi)害能力，保障寧夏枸杞產(chǎn)業(yè)健康發(fā)展。

多位研究者以氣象學(xué)上的霜凍致災(zāi)指標(biāo)（地面最低溫度 ≤ 0 ℃，最低氣溫 ≤ 0 ℃或最低氣溫 ≤ 2 ℃）來分析區(qū)域霜凍變化特征[6-13]，但對(duì)具體特色作物而言，因作物的種類、品種不同，其對(duì)霜凍災(zāi)害的耐受性也不同， 普適性的霜凍致災(zāi)指標(biāo)還不能準(zhǔn)確反映其受凍情況，鑒于此，部分學(xué)者利用茶樹[14-16]、蘋果花期[17]的凍害指標(biāo)分析得出了區(qū)域內(nèi)特色作物的霜凍變化規(guī)律。目前針對(duì)枸杞霜凍氣候特征方面的研究還鮮有報(bào)道，本試驗(yàn)在前人研究的基礎(chǔ)上，擬利用枸杞霜凍的致災(zāi)指標(biāo)，統(tǒng)計(jì)分析寧夏枸杞種植區(qū)霜凍災(zāi)害發(fā)生的時(shí)空分布特征，以期為優(yōu)化寧夏枸杞種植產(chǎn)業(yè)布局、提升枸杞種植區(qū)霜凍災(zāi)害的預(yù)測(cè)服務(wù)能力提供參考依據(jù)。

1 材料與方法

本研究所用惠農(nóng)、平羅、陶樂、賀蘭、銀川、永寧、靈武、吳忠、青銅峽、同心、中衛(wèi)、中寧、海原、興仁、固原15個(gè)站點(diǎn)的1961年到2017年逐日最低氣溫觀測(cè)資料來源于寧夏氣象信息中心，所選站點(diǎn)均為寧夏枸杞種植區(qū)內(nèi)具有長序列氣象資料的站點(diǎn)，站點(diǎn)分布見圖1;筆者所在項(xiàng)目組通過調(diào)查多年春季受災(zāi)情況，結(jié)合人工霜箱模擬和野外霜凍模擬試驗(yàn)，擬定日最低氣溫Tmin ≤-1 ℃作為寧夏枸杞春季霜凍的致災(zāi)指標(biāo)，本研究也將采用這一指標(biāo)開展統(tǒng)計(jì)分析;由于寧夏枸杞正常年份在4月下旬老眼枝開始顯蕾，所以春霜凍統(tǒng)計(jì)時(shí)段選定為4月21日至5月31日;將枸杞春季最后1次發(fā)生霜凍的日期定義為晚霜凍日。

本試驗(yàn)采用反距離權(quán)重法進(jìn)行空間插值分析寧夏枸杞種植區(qū)春霜凍發(fā)生次數(shù)的空間分布特征[18]，運(yùn)用ArcGIS 9.3制作了春霜凍發(fā)生次數(shù)的空間分布圖;運(yùn)用線性傾向估計(jì)分析枸杞春霜凍發(fā)生的年際變化特征[19];運(yùn)用Morlet小波變換分析逐年春霜凍發(fā)生總次數(shù)周期性變化[20-22];對(duì)1961—2017年寧夏逐年春霜凍發(fā)生總次數(shù)資料進(jìn)行標(biāo)準(zhǔn)化處理，采用經(jīng)驗(yàn)正交函數(shù)（EOF）綜合分析寧夏枸杞種植區(qū)春霜凍發(fā)生的時(shí)空變化特征[23-24]。

2 結(jié)果與分析

2.1 寧夏枸杞春霜凍發(fā)生的空間分布特征

2.1.1 枸杞春霜凍發(fā)生次數(shù)的空間分析 圖2反映了寧夏枸杞種植區(qū)1961—2017年春霜凍發(fā)生次數(shù)的空間分布情況。可以看出，寧夏枸杞種植區(qū)春霜凍具有南部地區(qū)發(fā)生次數(shù)多、中間區(qū)域發(fā)生次數(shù)少的空間分布特征，具體情況如下：中衛(wèi)、同心以南的區(qū)域春霜凍發(fā)生次數(shù)總體偏多，包括興仁、海原、固原3地發(fā)生春霜凍的次數(shù)占總數(shù)的46%，其中發(fā)生次數(shù)的最多的區(qū)域?yàn)榕d仁，57年來共發(fā)生了113次，其次是固原，發(fā)生了60次;而北部的陶樂地區(qū)也是春霜凍發(fā)生相對(duì)較多的區(qū)域，57年來也發(fā)生了47次，其他地區(qū)均未超過40次;春霜凍發(fā)生較少的區(qū)域主要位于永寧、青銅峽、中寧等地，分析期內(nèi)共發(fā)生了36次，其中永寧、中寧發(fā)生次數(shù)最少為11次。

2.1.2 枸杞晚霜凍日的空間分析 表1給出了1961—2017年寧夏枸杞種植區(qū)各站平均和最晚晚霜凍日的具體時(shí)間，可以看出，各站平均晚霜凍日結(jié)束時(shí)間除固原外均在4月下旬的中期到末期，平均晚霜凍日最晚的是固原，為5月2日，最早的是平羅和青銅峽，為4月24日;但從各站晚霜凍結(jié)束日期來看，最晚的區(qū)域集中在南部地區(qū)的興仁、海原、固原等地，都在5月中旬以后，結(jié)束最晚的是興仁，為5月28日;而北部地區(qū)的平羅、陶樂等地的最晚晚霜凍日也相對(duì)較晚，集中在5月中旬左右，其他地區(qū)除靈武外，最晚晚霜凍日都在5月上旬，最早的為永寧、中寧、同心，結(jié)束日期為5月4日。

2.2 寧夏枸杞春霜凍發(fā)生的時(shí)間變化

2.2.1 枸杞春霜凍發(fā)生的年際變化 由圖3可見，1961—2017年寧夏枸杞種植區(qū)逐年春霜凍發(fā)生次數(shù)呈減少趨勢(shì)，氣候傾向率為2.1次/10年，并且年際變化差異大，發(fā)生次數(shù)最多的年份為1976年，高達(dá)48次，而57年中有17年未發(fā)生春霜凍災(zāi)害;特別是從20世紀(jì)90年代初開始，春霜凍的發(fā)生次數(shù)明顯減少，除2008年超過了10次外，其他年份均未超過10次，21世紀(jì)10年代開始更是連續(xù)3年未發(fā)生枸杞春霜凍。

從逐年晚霜凍日年際變化來看（圖4），在57年尺度上，晚霜凍日呈提前趨勢(shì)，氣候傾向率為 2.1 d/10年，晚霜凍日最晚日期在5月28日，分別發(fā)生在1961年和1975年，最早日期為4月23日，發(fā)生在1988年，晚霜凍日較晚的時(shí)期主要集中在20世紀(jì)70年代到80年代初期，都在5月上旬后結(jié)束， 而從20世紀(jì)90年代末開始，晚霜凍日明顯提前，基本都在5月上旬前結(jié)束，只有3年是在5月上旬后結(jié)束，最晚結(jié)束的發(fā)生在2015年，為5月22日。

2.2.2 枸杞春霜凍發(fā)生次數(shù)的周期變化及突變分析 對(duì)寧夏枸杞種植區(qū)1961—2017年春霜凍逐年發(fā)生次數(shù)的小波分析圖結(jié)合小波方差圖進(jìn)行分析（圖5），結(jié)果顯示，在90年代之前，枸杞春霜凍發(fā)生次數(shù)主要存在7年的周期性變化，其次有4年的周期性變化，還存在較弱的15、23年的周期變化，而90年代之后的周期性變化規(guī)律不明顯。

計(jì)算57年來寧夏春霜凍發(fā)生次數(shù)的累計(jì)距平可知，枸杞春霜凍發(fā)生次數(shù)有2個(gè)突變點(diǎn)，即1972、1992年。從1972年開始春霜凍發(fā)生次數(shù)逐年增加，而1992年之后又逐年減少。利用滑動(dòng)t檢驗(yàn)方法進(jìn)行突變分析，為減少人為主觀性，本研究選擇2～20年，共18個(gè)時(shí)間長度，得到霜凍發(fā)生次數(shù)的突變點(diǎn)為1992年，該突變點(diǎn)通過了0.001的置信度檢驗(yàn)。

2.3 寧夏枸杞種植區(qū)春霜凍發(fā)生次數(shù)EOF分析

為了將時(shí)間和空間結(jié)合起來對(duì)寧夏枸杞種植區(qū)逐年春霜

凍發(fā)生次數(shù)的變化特征進(jìn)行綜合分析，對(duì)1961—2017年寧夏枸杞種植區(qū)各站逐年春霜凍發(fā)生次數(shù)資料進(jìn)行標(biāo)準(zhǔn)化處理后，再進(jìn)行EOF分解，所得方差貢獻(xiàn)率見表2?？梢钥闯?，前3個(gè)模態(tài)方差貢獻(xiàn)率分別為63.90%、8.60%、6.49%，累積方差貢獻(xiàn)率達(dá)到了78.99%，且3個(gè)模態(tài)均通過了North檢驗(yàn)，所以本研究取前3個(gè)模態(tài)來反映寧夏枸杞春霜凍發(fā)生次數(shù)的時(shí)空分布型。

圖6-a、圖6-b給出了第1模態(tài)的空間分布及時(shí)間序列，可以看出57年來春霜凍發(fā)生呈現(xiàn)出一致性變化，即整個(gè)枸杞種植區(qū)春霜凍發(fā)生一致偏多或者一致偏少，其中吳忠、青銅峽及以北大部分地區(qū)的載荷向量值普遍較大，說明該區(qū)域春霜凍發(fā)生次數(shù)的年際變化幅度大。從對(duì)應(yīng)的時(shí)間序列可以看出，1962、1968、1976年寧夏枸杞春霜凍發(fā)生次數(shù)偏多明顯，而到1992年之后呈現(xiàn)一致偏少的狀態(tài)，這與前文中突變檢驗(yàn)結(jié)果一致。

圖6-c、6-d為第2模態(tài)空間分布及對(duì)應(yīng)的時(shí)間序列，可以看出在個(gè)別年份寧夏春霜凍呈現(xiàn)出南北差異分布。其中，1961、1975、1982年寧夏中寧、中衛(wèi)以北地區(qū)發(fā)生次數(shù)偏少明顯，以南地區(qū)偏多明顯，而 1962、1968年正好相反，中寧、中衛(wèi)以北地區(qū)偏多明顯，而以南地區(qū)相對(duì)偏少。

圖6-e、6-f為第3模態(tài)空間分布及對(duì)應(yīng)的時(shí)間序列，可以看出寧夏春霜凍的發(fā)生在極個(gè)別年份呈現(xiàn)出3中心差異分布。其中，1981、1990年惠農(nóng)、吳忠發(fā)生次數(shù)偏多較為明顯，海原發(fā)生次數(shù)相對(duì)偏少;而1962、1974、2009年則相反。

3 結(jié)論與討論

寧夏枸杞春霜凍發(fā)生的區(qū)域特征明顯，春霜凍發(fā)生次數(shù)呈現(xiàn)南北多、中間少的空間分布特征，發(fā)生次數(shù)較多的區(qū)域?yàn)槟喜康呐d仁和固原以及北部的陶樂地區(qū)，而發(fā)生較少的區(qū)域?yàn)橛缹帯⒅袑幍鹊?而從晚霜凍日的空間分布特征來看，平均和最晚晚霜凍日結(jié)束最晚的區(qū)域還是集中在南部的固原和興仁，其他地區(qū)平均晚霜凍日結(jié)束日期主要集中在4月下旬中期到末期，而最晚晚霜凍日結(jié)束最早的主要區(qū)域仍然是用永寧、中寧、同心等地。總體來看，寧夏枸杞春霜凍防御的重點(diǎn)區(qū)域應(yīng)該集中在南部地區(qū)，而被譽(yù)為“中國枸杞之鄉(xiāng)”的中寧等地，春霜凍影響較輕，對(duì)于發(fā)展枸杞產(chǎn)業(yè)比較有利。

從時(shí)間變化來看，57年來枸杞晚霜凍日以2.1 d/10年的速率提前，且階段性變化特征明顯，20世紀(jì)70年代到80年代初期晚霜凍結(jié)束日期在5月上旬以后，而進(jìn)入21世紀(jì)以來，終霜凍日明顯提前，大都在5月上旬前結(jié)束;57年來枸杞春霜凍發(fā)生次數(shù)呈減少趨勢(shì)，氣候傾向率為-2.1次/10年，20世紀(jì)90年代前枸杞霜凍發(fā)生較多，之后明顯減少，1992年發(fā)生了明顯的突變;57年來寧夏枸杞春霜凍發(fā)生次數(shù)在20世紀(jì)90年代之前存在明顯的7年、準(zhǔn)4年的周期變化，而從90年代初期以后這種周期性變化不明顯。

寧夏枸杞春霜凍發(fā)生次數(shù)的第1模態(tài)的分布型為主導(dǎo)分布型，即大多數(shù)年份內(nèi)春霜凍發(fā)生次數(shù)呈現(xiàn)一致性變化特征，一致偏多或一致偏少，且1992年之前呈振蕩變化而之后呈明顯偏少的特征，吳忠、青銅峽及以北大部地區(qū)春霜凍發(fā)生次數(shù)年際變化幅度較大;從第2、第3模態(tài)可以看出，寧夏枸杞春霜凍發(fā)生次數(shù)還呈現(xiàn)南北差異分布，有些年份中衛(wèi)、中寧以南地區(qū)偏多明顯，而以北地區(qū)偏少明顯，有些年份正好相反。

參考文獻(xiàn)：

[1]董永祥，周仲顯.寧夏農(nóng)業(yè)氣候與農(nóng)業(yè)，4[M]. 銀川：寧夏人民出版社，1986：59-63.

[2]桑建人，舒志亮，邱 旺.近44a寧夏霜凍特征[J]. 干旱區(qū)資源與環(huán)境，2008，22（12）：97-102.

[3]張 磊，段曉鳳，李紅英，等. 寧夏枸杞生長的氣象條件分析及管理措施[J]. 北方果樹，2014（4）：16-19，前插1-前插2.

[4]葉力勤，趙銀寶.枸杞頭茬果質(zhì)量初探[J]. 中藥材，2004，27（1）：3-4.

[5]張自萍，郭 榮，史曉文，等. 不同采摘期枸杞品質(zhì)變化研究[J]. 西北農(nóng)業(yè)學(xué)報(bào)，2007，16（4）：126-130.

[6]葉殿秀，張 勇.1961—2007年我國霜凍變化特征[J]. 應(yīng)用氣象學(xué)報(bào)，2008，19（6）：661-665.

[7]許 艷，王國復(fù)，王盤興.近50a中國霜期的變化特征分析[J]. 氣象科學(xué)，2009，29（4）：4427-4433.

[8]王秀萍，任國玉，趙春雨，等. 近46年大連地區(qū)初、終霜凍事件和無霜凍期變化[J]. 應(yīng)用氣象學(xué)報(bào)，2008，19（6）：673-678.

[9]李 芬，張建新，閆永剛，等. 山西近50年初霜凍的時(shí)空分布及其突變特征[J]. 中國農(nóng)業(yè)氣象，2012，33（3）：448-456.

[10]杜 軍，石 磊，袁 雷. 近50年西藏主要農(nóng)區(qū)霜凍指標(biāo)的變化特征[J]. 中國農(nóng)業(yè)氣象，2013，34（3）：264-271.

[11]張旭暉，居為民，蒯志敏，等. 江蘇春季霜凍氣候變化特征及其未來可能變化趨勢(shì)[J]. 大氣科學(xué)學(xué)報(bào)，2013，36（6）：666-673.

[12]李 娜，張建新，李 芬. 山西省冬小麥種植區(qū)霜凍氣候指標(biāo)的時(shí)空特征分析[J]. 中國農(nóng)業(yè)氣象，2012，32（3）：436-441.

[13]陳少勇，夏 權(quán)，王勁松，等. 西北地區(qū)晚霜凍結(jié)束日的氣候變化特征及其影響因子[J]. 中國農(nóng)業(yè)氣象，2013，34（1）：8-13.

[14]姜燕敏，金志鳳，李仁忠，等. 浙南春茶早春霜凍的時(shí)空分布特征[J]. 氣象科技，2016，44（6）：1066-1070.

[15]許映蓮，張旭暉，李旭群. 蘇南茶區(qū)晚霜凍與茶樹凍害的研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2012，40（8）：236-238.

[16]吳 楊，金志鳳，葉建剛，等. 浙江茶樹春霜凍發(fā)生規(guī)律及其與太平洋海溫的遙相關(guān)分析[J]. 中國農(nóng)業(yè)氣象，2014，35（4）：434-439.

[17]李美榮，朱 琳，杜繼穩(wěn). 陜西蘋果花期霜凍災(zāi)害分析[J]. 果樹學(xué)報(bào)，2008，25（5）：666-670.

[18]張 宏，溫永寧，劉愛利. 地理信息系統(tǒng)算法基礎(chǔ)[M]. 北京：科學(xué)出版社，2006.

[19]普宗朝，張山清，李景林，等. 近48年新疆烏-昌地區(qū)氣候變化[J]. 干旱區(qū)研究，2010，27（3）：422-432.

[20]萬明波，程 智，王 文. 青藏鐵路沿線氣溫和降水的小波分析[J]. 干旱氣象，2006，24（4）：35-39.

[21]劉曉梅，閔錦忠，劉天龍. 新疆葉爾羌河流域溫度與降水序列的小波分析[J]. 中國沙漠，2009，29（3）：566-570.

[22]張維敏，王景紅. 近40年黃土高原丘陵紅棗種植區(qū)氣溫變化特征及其影響分析[J]. 中國農(nóng)學(xué)通報(bào)，2016，32（10）：139-144.

[23]李 芬，張建新. 山西近50年初霜凍的區(qū)域特征[J]. 災(zāi)害學(xué)，2013，28（4）：104-109，181.

[24]魏鳳英. 現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M]. 北京：氣象出版社，1999：106-113.李淏源，邱新法，曾 燕，等. 溫度植被干旱指數(shù)、 干燥度指數(shù) 、 水體分布與干濕特征的關(guān)系[J]. 江蘇農(nóng)業(yè)科學(xué)，2019，47（6）：243-249.