宋雄剛 王鴻斌 李國宏

(國家林業(yè)局森林保護(hù)學(xué)重點(diǎn)實(shí)驗(yàn)室(中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所))，北京，100091)

于志軍 陳國發(fā) 孔祥波 張真

(國家林業(yè)局森林病蟲害防治總站) (國家林業(yè)局森林保護(hù)學(xué)重點(diǎn)實(shí)驗(yàn)室(中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所))

責(zé)任編輯:程 紅。

農(nóng)林病蟲害是影響農(nóng)林生產(chǎn)及其生態(tài)系統(tǒng)穩(wěn)定的 重 大 自 然 災(zāi) 害 之 一［1－2］。松 毛 蟲(Dendrolimus spp．)是我國松屬樹種主要的食葉害蟲，大發(fā)生時(shí)會(huì)取食全部針葉，并引起人畜中毒;油松毛蟲(D．tabulaeformis)是松毛蟲中危害嚴(yán)重者之一，主要在北方發(fā)生［3］5，［4］。昆蟲生活史的各個(gè)階段對(duì)溫度都較敏感［5］，種群數(shù)量的變化及暴發(fā)規(guī)律除受到自身遺傳特性的影響外，還受到環(huán)境因素的重大影響［3］19，［6］166－192，［7－8］;氣候因素是影響松毛蟲種群發(fā)生規(guī)律的主要環(huán)境因素，不同的氣候因子在松毛蟲發(fā)育的不同時(shí)期可單獨(dú)或交互的起到重要作用［3］22－23，［8］。在氣象因素和松毛蟲的發(fā)生具有顯著相關(guān)性的基礎(chǔ)上，根據(jù)多種統(tǒng)計(jì)學(xué)或統(tǒng)計(jì)力學(xué)原理可建立各種監(jiān)測(cè)預(yù)報(bào)模型。張玉書等運(yùn)用氣象期距及有效積溫法對(duì)阜新縣油松毛蟲的始見期預(yù)測(cè)，其預(yù)測(cè)結(jié)果比實(shí)際上樹期晚一天［9］;周廣學(xué)等運(yùn)用相關(guān)分析和主成分分析法建立朝陽市油松毛蟲越冬死亡率和發(fā)生量的預(yù)測(cè)模型，結(jié)果顯示對(duì)死亡率影響最大的是春季溫度、對(duì)發(fā)生量影響最大的是前一年10月份降水和當(dāng)年1月中旬平均溫度［10］;陳繪畫等運(yùn)用逐步回歸法和人工神經(jīng)網(wǎng)絡(luò)模型建立了仙居縣災(zāi)害預(yù)測(cè)模型，其準(zhǔn)確性能達(dá)到100%［11］;Liang et al．用一般線性模型(GLM)研究了山松大小蠹(Dendroctonus ponderosae)在不同氣候情形下死亡率分布［12］。上述均以氣象因子的月或旬值為基礎(chǔ)展開研究，而對(duì)氣候因子動(dòng)態(tài)變化及氣候因子在昆蟲不同發(fā)育期的累積效應(yīng)考慮較少;費(fèi)海澤等考慮了上述影響后用人工神經(jīng)網(wǎng)絡(luò)法建立廣西馬尾松毛蟲(Dendrolimus punctatus)不同發(fā)生程度的預(yù)測(cè)模型且準(zhǔn)確性較高［8］。氣候因子變化對(duì)昆蟲的影響還體現(xiàn)在災(zāi)害發(fā)生期的提前或分布區(qū)的擴(kuò)大上［13－16］。上述研究大多以縣級(jí)或具有相同立地條件的區(qū)域?yàn)檠芯繉?duì)象，在小區(qū)域或相同的某一棲息環(huán)境中，可得到較準(zhǔn)確的結(jié)果;而在包括多種立地、林型的更大尺度上或更寬區(qū)域，則其預(yù)測(cè)結(jié)果代表性下降。故本研究以遼寧、北京、河北和山西等主要發(fā)生區(qū)為代表，以油松毛蟲2002—2011年的災(zāi)害監(jiān)測(cè)數(shù)據(jù)(包括監(jiān)測(cè)面積和不同發(fā)生程度)和相應(yīng)的氣候數(shù)據(jù)為基礎(chǔ)，結(jié)合氣候因子的動(dòng)態(tài)變化和油松毛蟲的不同生活階段，運(yùn)用典型相關(guān)分析法和主成分分析法，分析北方區(qū)域尺度上影響油松毛蟲發(fā)生的主要物候因子。這些因子在氣候變化背景下分析預(yù)測(cè)油松毛蟲的災(zāi)情變化有重要意義。

1 研究區(qū)域概況

研究區(qū)由遼寧、北京、河北及山西4個(gè)省行政區(qū)組成。該區(qū)域?yàn)橛退擅x在我國北方的主要發(fā)生區(qū)。我國近12年來，油松毛蟲的發(fā)生面積共158．69 hm2，其中2002年最為嚴(yán)重。從2009年開始油松毛蟲發(fā)生面積圍繞10 hm2波動(dòng)(圖1)(數(shù)據(jù)來源于森防總站的歷年監(jiān)測(cè)數(shù)據(jù)和總站的內(nèi)部資料《林業(yè)有害資料監(jiān)測(cè)預(yù)警年報(bào)》)，到2014年遼寧凌源大面積發(fā)生，松林鎮(zhèn)成災(zāi)，至2015年凌源持續(xù)大發(fā)生，單株蟲口密度達(dá)到2 000頭(數(shù)據(jù)由遼寧凌源森防站提供)，所以油松毛蟲的防治依然面臨嚴(yán)峻形勢(shì)。油松毛蟲的災(zāi)害性暴發(fā)主要出現(xiàn)于我國北方的9個(gè)省市，且以遼寧省發(fā)生面積最多，達(dá)到33%，天津最少，接近0(表1)。研究區(qū)的氣候因子、地形地勢(shì)、寄主情況能反應(yīng)油松毛蟲在我國北方分布與發(fā)生的一般特征(圖2)。

圖1 2002—2013年全國油松毛蟲發(fā)生面積變化趨勢(shì)

圖2 研究區(qū)地理地勢(shì)、松屬分布示意圖

表1 2002—2013年各省市油松毛蟲平均災(zāi)害面積及發(fā)生比例

2 數(shù)據(jù)資料來源

數(shù)據(jù)資料包括氣象數(shù)據(jù)和蟲情數(shù)據(jù)。其中氣象數(shù)據(jù)由中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc．nmic．cn/home．do)和中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所(http://www．isi－mip．org)提供，包括日、月、年的溫度、相對(duì)濕度、降雨、日照和風(fēng)速等;2002—2011年的油松毛蟲縣市級(jí)蟲情數(shù)據(jù)(包括監(jiān)測(cè)面積和各種程度的發(fā)生面積)由中國森防總站提供。

3 研究方法

3．1 衍生物候因子

氣候因子影響油松毛蟲生活史的各個(gè)期。油松毛蟲在研究區(qū)1年1代，且跨年度完成發(fā)育，即油松毛蟲的災(zāi)害情況和前一年、當(dāng)年的環(huán)境因素相關(guān)，為了更科學(xué)的分析氣候條件對(duì)油松毛蟲世代生長發(fā)育和其暴發(fā)的影響，結(jié)合研究區(qū)油松毛蟲的各個(gè)生長發(fā)育期對(duì)自然年界進(jìn)行如下劃分:首先將油松毛蟲的生活史劃分為4個(gè)生長發(fā)育期，即羽化—產(chǎn)卵—孵化—幼齡幼蟲期、越冬期、上樹期和老熟幼蟲食葉—化蛹期;然后將這4個(gè)期松毛蟲活動(dòng)特征和氣候結(jié)合并分別對(duì)應(yīng)為其生長發(fā)育的Ⅰ、Ⅱ、Ⅲ和Ⅳ期(表2)。結(jié)合油松毛蟲生活史不同時(shí)期對(duì)不同氣象因子不同的敏感性和相關(guān)的研究文獻(xiàn)［3］5－25、［5］、［7］，衍生物候因子。

表2 油松毛蟲生長期對(duì)應(yīng)時(shí)間

3．2 蟲情災(zāi)害數(shù)據(jù)資料

在2002—2011年的縣級(jí)數(shù)據(jù)庫中，選擇研究區(qū)中具有完整數(shù)據(jù)記錄的縣市級(jí)災(zāi)害發(fā)生點(diǎn)77個(gè)，并分別計(jì)算輕、中和重度發(fā)生程度的面積比例(p)。p為各個(gè)發(fā)生程度的發(fā)生面積和監(jiān)測(cè)面積的比值。計(jì)算公式如下:

式中:Pij為i點(diǎn)j災(zāi)害程度比例;i為發(fā)生點(diǎn);j為不同的發(fā)生程度，其中1代表輕度發(fā)生、2代表中度發(fā)生、3代表重度發(fā)生;Sij為i點(diǎn)j發(fā)生程度災(zāi)害面積;Si為i點(diǎn)的監(jiān)測(cè)面積。

3．3 物候因子篩選

典型相關(guān)分析法是利用綜合變量對(duì)之間的相關(guān)關(guān)系來反應(yīng)2組變量間的整體相關(guān)性的多元降維統(tǒng)計(jì)方法［17］，可用來篩選與油松毛蟲災(zāi)害發(fā)生情況相關(guān)的衍生物候因子。考慮到單一統(tǒng)計(jì)方法在相關(guān)性分析中存在的缺陷，同時(shí)運(yùn)用主成分分析法對(duì)與油松毛蟲災(zāi)害發(fā)生的相關(guān)氣候因子作分析，經(jīng)過2種方法的結(jié)合，得到主要物候因子。其一般過程為:首先對(duì)所有的衍生因子作正態(tài)性檢測(cè);其次用典型相關(guān)分析篩選衍生因子，并對(duì)顯著性進(jìn)行檢驗(yàn);最后用主成分分析法得出相關(guān)結(jié)果，并與典型相關(guān)分析結(jié)果做比較分析。用SPSS 19．0和Correl函數(shù)進(jìn)行數(shù)據(jù)分析。用典型相關(guān)分析篩選物候衍生因子的原理過程詳見文獻(xiàn)［17］。

4 結(jié)果與分析

4．1 衍生的物候因子

根據(jù)相關(guān)研究和本研究設(shè)定的條件，共衍生出可能和油松毛蟲災(zāi)害發(fā)生相關(guān)的物候因子73個(gè)。這些因子有溫度、濕度、降水、日照時(shí)間和風(fēng)速等(表3)，包括了油松毛蟲在北方區(qū)域中完成生活史不同階段的基礎(chǔ)物候因子和極端物候條件，如Ⅰ期平均溫度(X6)指前一年的8、9、10月份的平均溫度，此期間是油松毛蟲幼齡幼蟲進(jìn)食期，適宜的溫度是其能否完成各個(gè)期的基礎(chǔ)條件，而期間的濕度、風(fēng)速和最大日降雨(X43、X59、X62等)等條件在很大程度上影響了成活率、擴(kuò)散范圍［3］19－25，可直接造成其種群分布和數(shù)量的波動(dòng)，從而影響其災(zāi)害的暴發(fā)。

表3 區(qū)域尺度油松毛蟲預(yù)報(bào)衍生因子

4．2 衍生因子典型相關(guān)性

用SPSS 19．0對(duì)衍生物候數(shù)據(jù)和油松毛蟲災(zāi)害發(fā)生比率數(shù)據(jù)矩陣塊分析計(jì)算，得到物候因子和油松毛蟲災(zāi)害情況的均值和標(biāo)準(zhǔn)差，同時(shí)分析得出典型相關(guān)系數(shù)及其檢驗(yàn)結(jié)果(表4)。由表4可見，第1、2典型相關(guān)系數(shù)分別為1．000和0．988，其統(tǒng)計(jì)檢驗(yàn)達(dá)到極顯著水平(P1＜0．001、P2＜0．001)。第3典型相關(guān)系數(shù)雖然也較大但其統(tǒng)計(jì)檢驗(yàn)無法達(dá)到顯著水平，故剔除。因此，選前2對(duì)典型變量來分析物候因子與油松毛蟲災(zāi)害發(fā)生之間的相關(guān)關(guān)系。

表4 預(yù)報(bào)因子對(duì)油松毛蟲災(zāi)害發(fā)生影響的典型相關(guān)分析

由于物候因子和蟲情數(shù)據(jù)的量綱不同，為了降低分析的誤差，故在相關(guān)分析中對(duì)變量進(jìn)行標(biāo)準(zhǔn)化處理。標(biāo)準(zhǔn)化處理后的第1組典型相關(guān)變量為:

在第1組典型相關(guān)變量中，U1主要受到溫度因素(X6、X16、X26、X28、X30)、濕度因素(X2、X17、X22、X29、X33、X35、X37)和風(fēng)速因子(X9)的影響。由偏相關(guān)系數(shù)分析可得Ⅲ期平均溫度(X16)、Ⅰ期平均溫度(X6)及構(gòu)成Ⅲ期各月份平均溫度(X26、X28、X30)為其主要影響因子。其中，X16、X26、X28、X30呈正相關(guān)，X6呈負(fù)相關(guān)。即在一定范圍內(nèi)油松毛蟲上樹活動(dòng)期溫度越高，災(zāi)害發(fā)生的概率就越大;反之，油松毛蟲災(zāi)害發(fā)生的概率就越小。V1主要受到輕度發(fā)生比率y1的影響。在上樹期溫度因子的作用下，油松毛蟲災(zāi)害的發(fā)生主要為輕度發(fā)生。

標(biāo)準(zhǔn)化處理后的第2組典型相關(guān)變量為:

在第2組典型相關(guān)變量中，U2主要受到Ⅳ期平均濕度X22及Ⅳ期各月份平均濕度(X33、X35、X37)的影響，且呈負(fù)相關(guān)性。即Ⅳ期的濕度越高，油松毛蟲災(zāi)害發(fā)生的概率就越低，Ⅳ期的濕度越低則其災(zāi)害發(fā)生的概率就會(huì)越高。因此，將Ⅳ期平均濕度稱為影響油松毛蟲災(zāi)害發(fā)生的第2因子。V2則主要決定于輕度發(fā)生比例和中度發(fā)生比例。此外，典型相關(guān)表達(dá)式中也反映出Ⅰ期風(fēng)速對(duì)油松毛蟲災(zāi)害發(fā)生的影響較大且呈現(xiàn)正相關(guān)。

綜合第1、第2典型相關(guān)變量，Ⅲ期平均溫度、2月份平均溫度、3月份平均溫度、4月份平均溫度、Ⅰ期平均溫度、Ⅰ期平均風(fēng)速構(gòu)成影響油松毛蟲災(zāi)害發(fā)生的主要物候因子。

4．3 衍生因子的主成分分析

典型相關(guān)分析可得出和油松毛蟲災(zāi)害發(fā)生最相關(guān)的物候因子，但同時(shí)也有可能漏掉一些相關(guān)因子而使得篩選的相關(guān)因子不全面，故對(duì)衍生因子運(yùn)用主成分分析法對(duì)其進(jìn)行補(bǔ)充。主成分分析法在保證信息可靠性的前提下，可通過少數(shù)因子反映全部因子的大部分信息(＞85%)。而當(dāng)分析過程中設(shè)定抽取主成分和衍生因子相等時(shí)，可保證篩選物候因子的全面性，故而可補(bǔ)充典型相關(guān)分析的不足。對(duì)上述73個(gè)因子運(yùn)用SPSS 19．0做主成分分析，當(dāng)選取前11個(gè)主成分時(shí)，累計(jì)貢獻(xiàn)率達(dá)到86．46%(表5)，故可用前11個(gè)成分來代替分析原來的73個(gè)因子。

表5 前11個(gè)主成分方差

將前11個(gè)主成分的載荷矩陣?yán)肧PSS 19．0的Transfrom計(jì)算得到11個(gè)主成分的特征向量。結(jié)果表明，第1主成分中變量系數(shù)絕對(duì)值在0．30～0．45的變量共33個(gè)，是對(duì)溫度和風(fēng)速的綜合。溫度與第1主成分呈正相關(guān)，相關(guān)順序由大到小依次為X1、X6、X16、X73、X60、X21、X11;風(fēng)速和第1主成分呈負(fù)相關(guān)，相關(guān)順序由大到小依次為X4、X9。第2主成分集中反應(yīng)了相對(duì)濕度的影響，其變量系數(shù)為0．30～0．37，由大到小依次為X2、X44、X12、X22、X17。第3主成分到第11主成分中未集中反應(yīng)某種因子的相關(guān)性，但反應(yīng)了多種因子的綜合影響，其偏相關(guān)系數(shù)都較小(＜0．1)。故在分析影響油松毛蟲災(zāi)害發(fā)生的相關(guān)因子時(shí)，以第1和第2主成分為主。

在衍生物候因子的基礎(chǔ)上，綜合典型相關(guān)分析和主成分分析結(jié)果，選擇溫度因子中的年平均溫度、Ⅰ期平均溫度、Ⅳ期平均溫度、Ⅱ期平均溫度、Ⅱ期平均最高溫、9月份最高溫，風(fēng)速因子中的年平均風(fēng)速、Ⅰ期平均風(fēng)速，濕度因子中的年平均濕度、9月份平均濕度、Ⅱ期平均濕度、Ⅳ期平均濕度、Ⅲ期平均濕度為我國北方油松毛蟲災(zāi)害發(fā)生預(yù)測(cè)預(yù)報(bào)上的主要因子。

4．4 篩選物候因子的Correl函數(shù)

由典型相關(guān)分析和主成分分析篩選出和油松毛蟲災(zāi)害發(fā)生關(guān)系最大的13個(gè)物候因子，并用Correl函數(shù)對(duì)這些因子與不同程度的災(zāi)害發(fā)生作相關(guān)性分析。

表6 篩選的物候因子與油松毛蟲不同發(fā)生比例間的相關(guān)系數(shù)(r)

篩選物候因子和不同災(zāi)害程度的發(fā)生相關(guān)性變化范圍較大(－0．74≤r≤0．93)。表6顯示，相關(guān)性最強(qiáng)的是Ⅱ期平均最高溫度X73和重度發(fā)生(r=0．93)，并呈正相關(guān);相關(guān)性最弱的是Ⅱ期平均濕度X12和輕度發(fā)生(r=－0．01)，并呈負(fù)相關(guān);同時(shí)，年平均溫度(X1)和輕度發(fā)生比例、Ⅱ期平均溫度(X11)和重度發(fā)生比例、Ⅲ期平均濕度(X17)和重度發(fā)生比例、9月份最高溫(X60)和重度發(fā)生比例都有較強(qiáng)的相關(guān)，其相關(guān)系數(shù)都在0．75附近。

綜上分析:與油松毛蟲輕度災(zāi)害發(fā)生呈正相關(guān)的主要物候因子由強(qiáng)到弱依次為X1、X17、X60、X21、X4、X9。與中度災(zāi)害發(fā)生呈正相關(guān)的主要物候因子由強(qiáng)到弱依次為X11、X6、X9。與重度災(zāi)害發(fā)生呈正相關(guān)的主要物候因子由強(qiáng)到弱依次為X73、X60、X11、X6、X1、X44，呈負(fù)相關(guān)的物候因子由強(qiáng)到弱依次為X17、X2，即越冬期高溫和上樹期低溫對(duì)油松毛蟲重度災(zāi)害的發(fā)生有重要影響。

5 結(jié)論與討論

研究表明，在眾多的物候因子中，溫度、風(fēng)速、濕度因子是影響我國北方油松毛蟲災(zāi)害發(fā)生的主要物候因子。這些因子包括基礎(chǔ)物候因子和極端物候因子，基礎(chǔ)物候因子有年平均溫度、Ⅰ期平均溫度、Ⅱ期平均溫度、Ⅳ期平均溫度、年均風(fēng)速、Ⅰ期平均風(fēng)速、年平均濕度、9月份平均濕度、Ⅱ期平均濕度、Ⅲ期平均濕度、Ⅳ期平均濕度，極端物候因子Ⅱ期平均最高溫、9月份最高溫。這些物候因子對(duì)油松毛蟲不同程度災(zāi)害發(fā)生響應(yīng)也不同。輕度發(fā)生與年平均溫度和上樹期的平均相對(duì)濕度相關(guān)性較強(qiáng)。中度發(fā)生與羽化—產(chǎn)卵—孵化—幼齡幼蟲期的平均溫度和越冬期平均溫度相關(guān)性較強(qiáng)。重度發(fā)生與越冬期最高溫、9月份最高溫和上樹期平均相對(duì)濕度相關(guān)性較強(qiáng)，即暖冬和早春干旱均會(huì)提高油松毛蟲重度災(zāi)害發(fā)生的可能性。本研究結(jié)果與前人眾多的物候因子方面的研究［10，13－14，18－20］各有異同。就松毛蟲而言，薛賢清［18］通過對(duì)25個(gè)縣、市馬尾松毛蟲發(fā)生情況及氣象因子分析，得出和馬尾松毛蟲災(zāi)害發(fā)生的主要?dú)庀笠蜃佑蓮?qiáng)到弱為平均氣溫、降水量、相對(duì)濕度;而謝偉忠［19］則認(rèn)為3、8月份的平均氣溫，4、5、6月份的平均濕度及其溫濕度系數(shù)是影響馬尾松毛蟲災(zāi)害發(fā)生的主要因子，5—7月份連續(xù)的暴雨可減少松毛蟲的危害。劉志明等［20］研究結(jié)果表明，相對(duì)濕度是影響落葉松毛蟲災(zāi)害發(fā)生的首要因子。對(duì)其他蟲種類而言，F(xiàn)lower et al．的研究證明前一年的低濕對(duì)云杉卷葉蛾的發(fā)生有促進(jìn)作用，而高濕則很少導(dǎo)致災(zāi)害的發(fā)生［13］;Bentz et al．則發(fā)現(xiàn)溫度的升高對(duì)紅脂大小蠹和山松大小蠹行為造成顯著影響［14］;Liang et al．報(bào)道冬季高溫、夏季少雨可提高山松大小蠹的死亡率［12］?？梢?，氣候條件尤其是溫濕度的變化會(huì)對(duì)不同林業(yè)昆蟲造成不同的影響。綜上所述，蟲種不同，氣候條件對(duì)其的作用亦不盡相同，但對(duì)于松毛蟲災(zāi)害的發(fā)生都受到溫濕度的顯著影響。

油松毛蟲災(zāi)害發(fā)生與溫度、相對(duì)濕度和風(fēng)速相關(guān)。油松毛蟲在北方以2～4齡幼蟲越冬，越冬期和早春復(fù)蘇期其幼蟲的自我調(diào)節(jié)能力及抗環(huán)境脅迫能力都較差，越冬期高溫可提高越冬代的成活率。高濕環(huán)境利于昆蟲病害的流行，從而在很大程度上對(duì)松毛蟲種群密度的增長起到抑制作用。在高濕環(huán)境下，赤松毛蟲(D．spectabilis)CPV在林間使用良好［21］，思茅松毛蟲(D．kikuchii)CPV病毒感染率可達(dá)到97%以上［22］。所以，經(jīng)過越冬期復(fù)蘇的幼蟲生命力和抵抗力不強(qiáng)極易感染病害，而上樹期的低濕可抑制昆蟲病害的流行提高蟲口基數(shù)。不同的風(fēng)速對(duì)油松毛蟲的災(zāi)害發(fā)生的影響主要體現(xiàn)在幼齡幼蟲期，過大的風(fēng)速如日風(fēng)速＞5 m·s－1或月風(fēng)速＞10 m·s－1則會(huì)造成幼蟲的大量死亡而降低幼蟲成活率［3］24，而適當(dāng)?shù)娘L(fēng)速會(huì)對(duì)松毛蟲的傳播起到促進(jìn)作用［23］。

松毛蟲災(zāi)害的暴發(fā)不僅和種群自身遺傳有關(guān)，還取決于其棲息地諸多環(huán)境因素。而眾多環(huán)境因素中，氣象因子的可量化性使其在模擬預(yù)測(cè)中地位突出。森林昆蟲完成生活史的各個(gè)時(shí)期，都要以適宜的氣候環(huán)境作為依托。故本研究在考慮油松毛蟲和氣象因子的相關(guān)性時(shí)，為了保證其生長發(fā)育的各個(gè)時(shí)期和對(duì)應(yīng)的氣象因子的最大同步性，最大程度上結(jié)合了我國北方區(qū)域?qū)嶋H氣候特征和油松毛蟲在該區(qū)不同時(shí)間的發(fā)育情況，將前一年8—10月份定義為Ⅰ期，11月份至次年1月份為Ⅱ期，2—4月份為Ⅲ期，5—7月份為Ⅳ期。同時(shí)，將各月的氣象因子以月為單位分別處理，并結(jié)合油松毛蟲生長發(fā)育不同時(shí)期的不同敏感性［5］，加入其極端氣候變量。從而可更好地篩選影響其災(zāi)害發(fā)生的主要因子，除上述篩選的主要物候因子外，降雨和光照也有一定作用。降雨的影響取決于降雨的強(qiáng)度和降雨的時(shí)間，月總降雨量大于500 mm或日降雨量超過100 mm蟲口密度會(huì)明顯降低［3］22，［6］169，［23］;松毛蟲對(duì)光周期的敏感性會(huì)引起松毛蟲滯育［24］，有利于其越冬，從而保證其蟲口密度并對(duì)災(zāi)害發(fā)生起積極作用;反之，則起到消極作用。但由于缺乏油松毛蟲長期的監(jiān)測(cè)數(shù)據(jù)，只能用2002—2011年相對(duì)完整的77組歷史數(shù)據(jù)進(jìn)行分析，其結(jié)果可能存在一定的偏差，為了更好地了解研究氣象因子對(duì)森林病蟲害的影響，應(yīng)加強(qiáng)監(jiān)測(cè)樣點(diǎn)布置的代表性和監(jiān)測(cè)數(shù)據(jù)量的積累。建立氣候因子與有害生物災(zāi)情之間的關(guān)系模型將是未來林業(yè)有害生物監(jiān)測(cè)預(yù)報(bào)的研究重點(diǎn)內(nèi)容之一。同時(shí)，以溫度和降水變化為主要特征的氣候變化背景下，結(jié)合氣象關(guān)鍵因子和其他的環(huán)境因子對(duì)林業(yè)有害生物災(zāi)害暴發(fā)［12］的評(píng)價(jià)和預(yù)報(bào)也是研究的重要內(nèi)容。

［1］段永春，袁洪剛，張永亮，等．山東茶樹蟲害的種類及主要品種的發(fā)生規(guī)律研究［J］．中國農(nóng)學(xué)通報(bào)，2010，26(11):284－289．

［2］王志偉，張東霞，馬雅麗，等．山西山西省冬小麥主要病蟲害氣象等級(jí)預(yù)報(bào)模型［J］．中國農(nóng)學(xué)通報(bào)，2010，26(11):267－271．

［3］陳昌潔．松毛蟲綜合管理［M］．北京:中國林業(yè)出版社，1990．

［4］曾菊平，戈峰，蘇建偉，等．我國林業(yè)重大害蟲松毛蟲的災(zāi)害研究進(jìn)展［J］．昆蟲知識(shí)，2010，47(3):451－459．

［5］Yamamura K，Yokozawa M，Nishimori M，et al．How to analyze long-term insect population dynamics under climate change:50-year data of three insect pests in paddy fields［J］．Population Ecology，2006，48(1):31－48．

［6］侯陶謙．中國松毛蟲［M］．北京:科學(xué)出版社，1989．

［7］韓瑞東，何忠，戈峰，等．馬尾松毛蟲黃河流影響松毛蟲種群動(dòng)態(tài)的因素［J］．昆蟲知識(shí)，2004，41(6):504－511．

［8］費(fèi)海澤，王鴻斌，孔祥波，等．馬尾松毛蟲發(fā)生相關(guān)氣象因子篩選及預(yù)測(cè)［J］．東北林業(yè)大學(xué)學(xué)報(bào)，2014，42(1):136－140．

［9］張玉書，馮銳，陳鵬獅，等．松毛蟲發(fā)生期與氣象條件關(guān)系［J］．中國農(nóng)業(yè)氣象，2004，25(3):26－28．

［10］周廣學(xué)，張國林，梁群，等．氣象條件對(duì)油松毛蟲的影響及其預(yù)測(cè)模型的構(gòu)建［J］．東北林業(yè)大學(xué)學(xué)報(bào)，2012，40(11):131－134．

［11］陳繪畫，朱壽燕，崔相富，等．基于人工神經(jīng)網(wǎng)絡(luò)的馬尾松毛蟲發(fā)生量預(yù)測(cè)模型的研究［J］．林業(yè)科學(xué)研究，2003，16(2):159－165．

［12］Liang L，Hawbaker T J，Chen Y L，et al．Characterizing recent and projecting future potential patterns of mountain pine beetle outbreaks in the southern rocky mountains［J］．Applied Geography，2014，55:165－175．

［13］Flower A，Gavin D G，Heyerdahl E K，et al．Drought-triggered western spruce budworm outbreaks in the Interior Pacific Northwest:A multi-century dendrochronological record［J］．Forest Ecology and Management，2014，24(3):16－27．

［14］Bentz B J，Régnière J，F(xiàn)ettig C J，et al．Climate change and bark beetles of the western United States and Canada:direct and indirect effects［J］．BioScience，2010，60(8):602－613．

［15］Musolin D L．Insects in a warmer world:ecological，physiological and life-history responses of true bugs(Heteroptera)to climate change［J］．Global Change Biology，2007，13(8):1565－1585．

［16］Kroschel J，Sporleder M，Tonnang H E Z，et al．Predicting climate-change-caused changes in global temperature on potato tuber moth(Phthorimaea operculella)distribution and abundance using phenology modeling and GISmapping［J］．Agricultural and Forest Meteorology，2013，17(6):228－241．

［17］田兵．典型相關(guān)分析及其應(yīng)用［J］．陰山學(xué)刊，2013，27(4):11－22．

［18］薛賢清．在馬尾松毛蟲測(cè)報(bào)中應(yīng)用逐步回歸電算方法的研究［J］．林業(yè)科學(xué)，1984，20(1):42－49．

［19］謝偉忠．溫濕度對(duì)廣東陽春市馬尾松毛蟲種群的影響［J］．中國森林病蟲，2005，24(4):17－19．

［20］劉志明，晏明，張旭東，等．氣象因子對(duì)落葉松毛蟲的影響分析［J］．災(zāi)害學(xué)，2002，17(4):34－37．

［21］陳昌潔，王志賢，劉革，等．赤松毛蟲質(zhì)型多角體病毒的引進(jìn)和利用研究［J］．林業(yè)科學(xué)研究，1988，1(1):14－27．

［22］曾述圣，楊桂，楊中學(xué)，等．思茅松毛蟲質(zhì)型多角體病毒的應(yīng)用［J］．森林病蟲通訊，2000，2(1):28－29．

［23］孫思清，邵守英，劉鴻巖．臺(tái)風(fēng)暴雨對(duì)赤松毛蟲蟲口密度的影響［J］．山東林業(yè)科技，1997(6):21－23．

［24］李兆麟，賈鳳友．油松毛蟲的光周期反應(yīng)［J］．昆蟲學(xué)報(bào)，1989，32(4):410－417．