摘 要：為了研究桂陽(yáng)煙區(qū)烤煙生產(chǎn)氣候適宜性，以1991—2020年桂陽(yáng)的逐日氣象觀測(cè)數(shù)據(jù)為基礎(chǔ)資料，采用趨勢(shì)分析、適宜性評(píng)價(jià)等方法，分析了桂陽(yáng)煙區(qū)氣候資源變化，烤煙種植的主要?dú)夂蛞蜃蛹翱緹煔夂蜻m宜性的變化。結(jié)果表明：桂陽(yáng)煙區(qū)氣候因子年際間差異較大，大田各生育期均溫顯著升高，其中伸根期均溫上升速率達(dá)0.128℃/a；大田期降雨量和大田期日照時(shí)數(shù)年際變化劇烈，但未通過顯著性檢驗(yàn)（P＞0.05）；氣候適宜性的分析結(jié)果顯示，桂陽(yáng)烤煙生產(chǎn)主要限制因子是伸根期均溫和大田期可用天數(shù)，均有向好趨勢(shì)；成熟期均溫和旺長(zhǎng)期降雨量適宜度較高，但呈變差的趨勢(shì)；桂陽(yáng)煙區(qū)烤煙氣候適宜性指數(shù)有所增加，但未通過顯著性檢驗(yàn)（P＞0.05）。綜上，近30 a桂陽(yáng)煙區(qū)氣候資源發(fā)生了較大變化，生產(chǎn)上必須趨利避害合理利用氣候的變化，合理布置好烤煙生產(chǎn)物候期。

關(guān)鍵詞：烤煙；氣候變化；氣候適宜性；桂陽(yáng)

中圖分類號(hào)：S572.06 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1006-060X（2023）06-0055-06

Abstract：In order to explore the suitable climate for tobacco production in Guiyang tobacco-growing area， using" the daily meteorological observation data of Guiyang from 1991 to 2020 as reference， the trend analysis and suitability evaluation methods were adopted to analyze the changes of climate resources， the main climate factors affecting tobacco planting and climate suitability for tobacco in Guiyang tobacco-growing area. The results showed that the annual difference of climatic factors in Guiyang tobacco-growing area was big， and the average temperature in each growth period of field increased significantly， among which the average temperature in root-extending period increased at a rate of 0.128 ℃/a; the rainfall and sunshine duration in the field period changed sharply year by year， but failed to pass the 95% significance test （Pgt;0.05）. Climate suitability analysis showed that the main factors affecting tobacco production in Guiyang were the average temperature in root-extending period and the available days in the field period， which had a positive trend; the average temperature in maturity period and the rainfall in vigorous growing period had high suitability， but showing a trend of deterioration; the climate feasibility index （CFI） of flue-cured tobacco in Guiyang increased， but failed to pass the 95% significance test （Pgt;0.05）. To sum up， the climate resources in Guiyang tobacco-growing area have changed greatly in recent 30 years. In production， we must make good use of the climate change and suitably arrange the phenological period of tobacco production.

Key words： tobacco; climate change; climate suitability; Guiyang

氣候是影響農(nóng)業(yè)生產(chǎn)的主要因素，高質(zhì)量的氣候條件對(duì)農(nóng)林業(yè)生產(chǎn)力的形成起著關(guān)鍵作用[1-2]。工業(yè)革命以來，全球性的氣候變化極其巨大[3-4]，而氣候的變化極易造成氣象災(zāi)害頻發(fā)[5-6]，病蟲害加劇[7-8]。氣候變化引起作物氣候適宜性變化的研究已有大量報(bào)道。越昆等[9]研究發(fā)現(xiàn)，40 a來內(nèi)蒙古中東部大部分玉米產(chǎn)區(qū)適宜氣溫增長(zhǎng)在0.005℃/a以上，致使綜合適宜度發(fā)展向好，但降水適宜度呈下降趨勢(shì)，并指出合理的農(nóng)業(yè)水資源利用是內(nèi)蒙古玉米穩(wěn)產(chǎn)增產(chǎn)的關(guān)鍵。胡春麗等[10]研究發(fā)現(xiàn)，近60 a來遼寧水稻各生育期溫度適宜度呈顯著上升趨勢(shì)，日照適宜度呈顯著下降趨勢(shì)，水稻全生育期平均氣溫每升高1℃，適宜度指數(shù)增加0.05。陳愛莉等[11]對(duì)楊梅在當(dāng)前和未來氣候下的適宜生長(zhǎng)區(qū)進(jìn)行了模擬，認(rèn)為在未來氣候背景下，廣東、云南、廣西和海南等地區(qū)的楊梅適生區(qū)面積會(huì)顯著銳減，而秦嶺—淮河一線適生區(qū)面積將出現(xiàn)大幅增加，總體上楊梅的適宜分布區(qū)有北擴(kuò)趨勢(shì)。煙草作為一種適應(yīng)能力很強(qiáng)的作物，在北緯60°至南緯40°的廣闊區(qū)域內(nèi)都可以完成其生命周期，但煙葉產(chǎn)量和品質(zhì)的形成對(duì)氣候條件的要求很高[12]。因此，有學(xué)者對(duì)氣候變化條件下烤煙種植的適宜性進(jìn)行了研究，胡穎梅等[13]的研究顯示，陜西各產(chǎn)區(qū)烤煙生育期間降雨量年際波動(dòng)較大，是烤煙種植的主要限制因素；多數(shù)產(chǎn)區(qū)氣溫適宜度呈顯著升高趨勢(shì)，寶雞和咸陽(yáng)產(chǎn)區(qū)的烤煙氣候適宜性顯著升高。胡雪瓊等[14]利用氣候模擬數(shù)據(jù)和氣象臺(tái)站實(shí)測(cè)數(shù)據(jù)， 擴(kuò)的趨勢(shì)，未來云南烤煙可種植區(qū)域?qū)⒊手饾u增加的趨勢(shì)，并認(rèn)為這主要與未來云南煙區(qū)的熱量分布變化有關(guān)。

桂陽(yáng)種煙歷史悠久，該煙區(qū)屬中亞熱帶向北亞熱帶過渡的大陸性季風(fēng)濕潤(rùn)氣候，四季分明，寒冷期短，炎熱期長(zhǎng)。桂陽(yáng)煙葉是廣受國(guó)內(nèi)卷煙企業(yè)的青睞，是南京、金圣、利群等多個(gè)重點(diǎn)骨干卷煙品牌的主體原料。開展氣候變化背景下桂陽(yáng)烤煙種植氣候適宜性的研究，對(duì)合理利用當(dāng)?shù)貧夂蛸Y源、挖掘烤煙生產(chǎn)氣候潛力具有重要的指導(dǎo)意義。筆者利用1991—2020年桂陽(yáng)氣象站的逐日地面氣象資料，采用趨勢(shì)分析、適宜性評(píng)價(jià)等方法，分析了桂陽(yáng)煙區(qū)氣候變化，不同時(shí)期烤煙種植的主要?dú)夂蛳拗埔蜃蛹翱緹煔夂蜻m宜性的變化，旨在明確氣候變化對(duì)桂陽(yáng)烤煙種植的影響，進(jìn)而趨利避害，順應(yīng)氣候的變化，合理優(yōu)化桂陽(yáng)烤煙的生產(chǎn)節(jié)令。

1 材料與方法

1.1 數(shù)據(jù)來源

氣象資料由湖南省氣象局提供，包括1991—2020年逐日平均氣溫、降雨量、日照時(shí)數(shù)和相對(duì)濕度等。根據(jù)當(dāng)?shù)乜緹熒a(chǎn)節(jié)令，先進(jìn)行烤煙生育期劃分：大田期（3月20日—7月15日）、伸根期（3月20日—4月18日）、旺長(zhǎng)期（4月19日～5月18日）和成熟期（5月19日—7月15日）。計(jì)算出相應(yīng)時(shí)段的烤煙氣候因子，用于分析的指標(biāo)有大田期均溫、大田期日照時(shí)數(shù)、大田期降雨量、≥10 ℃活動(dòng)積溫、伸根期均溫、成熟期均溫、旺長(zhǎng)期降雨量和成熟期降雨量等。

1.2 數(shù)據(jù)分析方法

1.2.1 趨勢(shì)分析 氣象要素的隨時(shí)間變化規(guī)律可用時(shí)間序列和自然序列的一元線性回歸方程來反映[15-16]。回歸方程如下。

y（t）=ati+b （i=1， 2， 3…， n）

式中：y（t） 為氣象要素的擬合值，t為時(shí)間序列；回歸系數(shù)a和常量b用最小二乘法求算?；貧w系數(shù)a的絕對(duì)值反映氣候要素的變化程度，a為負(fù)表示在氣象要素分析期內(nèi)呈下降趨勢(shì)，反之為上升趨勢(shì)；以a×10 作為氣候傾向率。用方程的決定系數(shù)（R2）和P值來檢驗(yàn)顯著性。

1.2.2 氣候適宜性分析 參考李強(qiáng)等[17]建立的烤煙氣候適宜性評(píng)價(jià)模型，采用隸屬函數(shù)計(jì)算單個(gè)氣候因素的隸屬度值表征該指標(biāo)的適宜度，再利用加乘法則計(jì)算烤煙綜合氣候適宜性指數(shù)（Climate Feasibility Index，CFI），計(jì)算公式如下。

式中，Nij和Wij分別為第i個(gè)評(píng)價(jià)單元的第j個(gè)評(píng)價(jià)指標(biāo)的隸屬度值和其權(quán)重。

2 結(jié)果與分析

2.1 桂陽(yáng)烤煙氣候變化特征

2.1.1 生育期均溫變化 利用桂陽(yáng)縣1991—2020年逐日氣溫?cái)?shù)據(jù)繪制30 a年限的烤煙伸根期、旺長(zhǎng)期、成熟期均溫演變趨勢(shì)圖（圖1～圖3）。伸根期均溫為15.94℃，最高為19.28 ℃（2014年），最低為10.48 ℃（1996年），變異系數(shù)為10.93%，為中等程度變異；在30 a中僅有2 a的伸根期均溫低于13 ℃，均出現(xiàn)在1991—2000年間（圖1）。由歷年伸根期均溫與時(shí)間序列擬合的趨勢(shì)線可知，伸根期均溫呈線性升高趨勢(shì)，并通過了顯著性檢驗(yàn)（P＜0.01），線性擬合傾向率達(dá)0.128 ℃/a。旺長(zhǎng)期均溫為21.36℃，最高為24.15 ℃（2005年），最低為17.90 ℃（1991年），變異系數(shù)為5.79%，為弱變異；在30 a中僅有1991年的旺長(zhǎng)期均溫低于18 ℃（圖2）。由歷年旺長(zhǎng)期均溫與時(shí)間序列擬合的趨勢(shì)線可知，旺長(zhǎng)期均溫呈線性升高趨勢(shì)，并通過了顯著性檢驗(yàn)（P＜0.05），線性擬合傾向率為0.046 ℃/a。成熟期均溫為26.06 ℃，最高為27.87 ℃（2020年），最低為24.22 ℃（1992年），變異系數(shù)為10.93%，為中等程度變異（圖3）；由歷年成熟期均溫與時(shí)間序列擬合的趨勢(shì)線可知，成熟期均溫呈線性升高趨勢(shì)，并通過了顯著性檢驗(yàn)（P＜0.05），線性擬合傾向率為0.035 ℃/a。這表明桂陽(yáng)煙區(qū)烤煙各生育期均溫均呈上升趨勢(shì)，其中伸根期均溫升高幅度最大，旺長(zhǎng)期均溫和成熟期均溫升高幅度較小。

2.1.2 生育期降雨變化特征 桂陽(yáng)煙區(qū)近30 a各生育期的降雨量演變趨勢(shì)如圖4～圖6所示，伸根期降雨量平均為178.87 mm，降雨量最多的年份2016年達(dá)378.70 mm，最少為24.20 mm（2011年），極差達(dá)354.50 mm，變異系數(shù)為44.14%，表現(xiàn)為中等程度變異（圖4）。旺長(zhǎng)期降雨量平均為187.04 mm，最高的年份為2019年，達(dá)313.20 mm，最低的年份為2020年，為85.40 mm，極差達(dá)227.8 mm，變異系數(shù)為29.16%，表現(xiàn)為中等變異（圖5）。成熟期降雨量平均為369.47 mm，最高的年份2006年達(dá)626.10 mm，最低的年份2011年為213.20 mm，極差達(dá)412.9 mm，變異系數(shù)為32.13 %，表現(xiàn)為中等變異（圖6）。

由各生育期降雨量與時(shí)間序列擬合的趨勢(shì)線可知，整個(gè)分析期內(nèi)伸根期降雨量和成熟期呈小幅下降，旺長(zhǎng)期降雨量呈小幅上升，但均未通過顯著性檢驗(yàn)（P＞0.05）。以上分析表明，桂陽(yáng)煙區(qū)各生育期降雨量年度間變化劇烈，與時(shí)間序列的關(guān)系無明顯規(guī)律。

2.1.3 大田期日照時(shí)數(shù)變化特征 近30 a大田期日照時(shí)數(shù)平均為553.03 h，變異系數(shù)為11.39%，為中等變異，其中日照時(shí)數(shù)最長(zhǎng)的年份為2000年，達(dá)717.40 h；最少的年份為2019年，僅為393.40 h，大田期日照時(shí)數(shù)最長(zhǎng)的年份和最短的年份相差324 h。與最適日照時(shí)數(shù)（＞600 h）相比，30 a中有4 a的日照時(shí)數(shù)多于600 h，分別出現(xiàn)在2000年、2007年、2013年和2014年（圖7）。歷年大田期日照時(shí)數(shù)與時(shí)間序列擬合的趨勢(shì)線未通過顯著性檢驗(yàn)（P＞0.05），表明大田期日照時(shí)數(shù)隨時(shí)間線性變化不明顯。

2.1.4 大田期可用天數(shù) 近30 a桂陽(yáng)大田期可用天數(shù)平均為116.17 d，其中可用天數(shù)最長(zhǎng)的年份2008年達(dá)146 d；最少的年份2007年僅有78 d，可用天數(shù)最多的年份和最少的年份相差68 d，變異系數(shù)為14.97%，為中等變異。與最少天數(shù)相比（90 d）相比，30 a中有2 a的日照時(shí)數(shù)少于90 d，分別出現(xiàn)在2007年和2010年（圖8）。歷年大田期可用天數(shù)與時(shí)間序列擬合的趨勢(shì)線可見，大田期可用天數(shù)呈線性增加趨勢(shì)，但未通過顯著性檢驗(yàn)（P＞0.05），表明大田期可用天數(shù)隨時(shí)間線性變化不明顯。

2.2 桂陽(yáng)煙區(qū)氣候適宜性評(píng)價(jià)

2.2.1 指標(biāo)選擇及隸屬函數(shù)構(gòu)建 參考相關(guān)文獻(xiàn)[12]并結(jié)合桂陽(yáng)氣候特點(diǎn)，選取伸根期均溫、成熟期均溫、旺長(zhǎng)期降雨量、大田期日照時(shí)數(shù)和大田期可用天數(shù)等5項(xiàng)指標(biāo)作為評(píng)價(jià)桂陽(yáng)氣候適宜性的因子?；谀：龜?shù)學(xué)原理，根據(jù)各參評(píng)指標(biāo)對(duì)烤煙生長(zhǎng)及煙葉品質(zhì)影響，同時(shí)結(jié)合桂陽(yáng)烤煙生產(chǎn)實(shí)際，構(gòu)建各參評(píng)指標(biāo)隸屬函數(shù)公式（1）拋物線形、公式（2）S形。在咨詢當(dāng)?shù)貙＜业幕A(chǔ)上采用層次分析法[18]確定各參評(píng)指標(biāo)的臨界值和權(quán)重（表1、表2）。

2.2.2 桂陽(yáng)煙區(qū)氣候適宜性評(píng)價(jià) 利用隸屬函數(shù)模型計(jì)算得到3個(gè)年代各指標(biāo)的隸屬度值（表3），隸屬度取值越大表示該指標(biāo)越利于烤煙生長(zhǎng)和煙葉品質(zhì)，反之則越差。從各指標(biāo)隸屬度值的年代際變化來看，伸根期均溫和大田期可用天數(shù)有向好的趨勢(shì)，成熟期均溫和旺長(zhǎng)期降雨量有變差的趨勢(shì)，大田期日照時(shí)數(shù)隸屬度值在年代間變化較小。桂陽(yáng)歷年氣候適宜性指數(shù)（climate feasibility index，CFI）通過采用加成法則計(jì)算得到（圖9），整個(gè)分析期內(nèi)歷年烤煙氣候適宜性指數(shù)呈緩慢升高的趨勢(shì)，但未通過顯著性檢驗(yàn)（P＞0.05），CFI的升高主要得益于伸根期均溫的升高和大田期可用天數(shù)的增加。

3 小結(jié)與討論

近30 a來，桂陽(yáng)煙區(qū)氣候變化明顯，主要表現(xiàn)在生育期均溫升高、降雨量和日照時(shí)數(shù)年際間變化劇烈但趨勢(shì)不明顯；桂陽(yáng)烤煙生產(chǎn)主要限制因子是伸根期均溫和大田期可用天數(shù)，這兩項(xiàng)指標(biāo)均有向好的趨勢(shì)；成熟期均溫和旺長(zhǎng)期降雨量適宜度目前較高，但有下降趨勢(shì)；桂陽(yáng)煙區(qū)烤煙氣候適宜性指數(shù)（CFI）有所增加。

聯(lián)合國(guó)政府間氣候變化專門委員會(huì)（IPCC）發(fā)布的氣候報(bào)告（AR6）指出，近代以來的氣候變化及其對(duì)地球的影響極其巨大[3]。氣候變化對(duì)農(nóng)業(yè)生產(chǎn)的影響是廣泛而深遠(yuǎn)的[1-8]，氣候的變化必然導(dǎo)致區(qū)域農(nóng)作物適宜性的變化[8-11]。煙草生長(zhǎng)是一種對(duì)氣候適應(yīng)能力很強(qiáng)的作物，但優(yōu)質(zhì)煙葉的生產(chǎn)對(duì)氣候條件的要求卻十分高，要求生育期內(nèi)光、溫、水、氣高效耦合。此研究總體上，近30 a桂陽(yáng)煙區(qū)氣候變化十分明顯。

從大田期均溫來看，桂陽(yáng)煙區(qū)各生育期均溫均呈上升趨勢(shì)，其中伸根期均溫升高最為明顯上升速率達(dá)0.128℃/a，高于桑植煙區(qū)伸根期均溫上升速率[16]，更是顯著高于中國(guó)近100 a和近50 a的年均氣溫上升速率[7， 19]。生育期期均溫的升高和有效積溫的增加有利于烤煙早生快發(fā)、形態(tài)建成和干物質(zhì)積累。從大田生育期降雨量的變化來看，桂陽(yáng)煙區(qū)各生育期降雨量總體上較為充沛，年度間變化劇烈，無明顯規(guī)律；充沛而不過的降雨量有利于煙株生長(zhǎng)發(fā)育，對(duì)煙葉品質(zhì)形成有利[17]，但個(gè)別年份降雨量過多則不利于烤煙生產(chǎn)。從大田期日照時(shí)數(shù)來看，桂陽(yáng)煙區(qū)日照時(shí)數(shù)平均為553.03 h，接近最適宜值600 h，但年際間差異較大，部分年份日照時(shí)數(shù)嚴(yán)重偏少（＜500 h），顯著低于國(guó)內(nèi)外優(yōu)質(zhì)煙區(qū)，不利于煙葉產(chǎn)量和品質(zhì)的形成。

進(jìn)行煙草氣候適宜性的綜合評(píng)價(jià)，有助于全面掌握煙區(qū)的有利因子和不利因子，對(duì)合理利用氣候資源，趨利避害，提高煙葉產(chǎn)質(zhì)量有重要的宏觀指導(dǎo)意義。此研究在進(jìn)行充分文獻(xiàn)檢索、煙區(qū)生產(chǎn)調(diào)查和咨詢專家意見的基礎(chǔ)上，確定了氣候適宜性評(píng)價(jià)的參評(píng)指標(biāo)、指標(biāo)閾值和權(quán)重，既考慮了烤煙生長(zhǎng)和品質(zhì)形成的一般規(guī)律，也兼顧了桂陽(yáng)煙區(qū)氣候特點(diǎn)。為便于標(biāo)準(zhǔn)統(tǒng)一，此研究劃定的大田期生育期為3月20日至7月15日，與當(dāng)前桂陽(yáng)煙區(qū)烤煙物侯期基本一致。從隸屬度值來看，桂陽(yáng)煙區(qū)主要限制因子是伸根期均溫和大田期可用天數(shù)，主要表現(xiàn)為伸根期均溫偏低，大田期可用天數(shù)偏短，這兩項(xiàng)指標(biāo)均有向好的趨勢(shì)。值得注意的是，雖然成熟期均溫的隸屬度值下降緩慢，但結(jié)合35℃高溫出現(xiàn)日期越來越早這一現(xiàn)象來看（數(shù)據(jù)未列出），桂陽(yáng)煙區(qū)成熟期烤煙高溫逼熟問題或?qū)㈦S著氣候變暖而逐步突顯[20]。旺長(zhǎng)期降雨量隸屬度值也呈下降趨勢(shì)，主要由旺長(zhǎng)期降雨過多并緩慢增加導(dǎo)致。從桂陽(yáng)烤煙氣候適宜性綜合評(píng)價(jià)的結(jié)果來看，整個(gè)分析期內(nèi)，CFI呈明顯增加的趨勢(shì)。從年代變化來看2001—2010年和2011—2020年CFI較1991—2000年分別平均增加了0.04和0.07，這主要得益于伸根期均溫的升高和大田期可用天數(shù)的增加。

參考文獻(xiàn)：

[1] 穆青云，李 俊，何 亮，等. 青藏高原冬小麥生產(chǎn)潛力及其對(duì)氣候變化的響應(yīng)[J]. 干旱區(qū)資源與環(huán)境，2021，35（7）：92-99.

[2] 韓艷剛，蓋學(xué)瑞，邱思玉，等. 大興安嶺興安落葉松徑向生長(zhǎng)對(duì)氣候響應(yīng)的時(shí)空變化[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2021，32（10）：3397-3404.

[3] IPCC. Summary for policymakers[A]. Climate Change 2013 – The Physical Science Basis[M]. Cambridge：Cambridge University Press，2014. 1-30.

[4] 中國(guó)氣象局氣候變化中心. 中國(guó)氣候變化藍(lán)皮書（2020）[M]. 北京：科學(xué)出版社，2020.

[5] 陳 昆，黃淑娥，景元書. 江西省早稻高溫逼熟規(guī)律及其對(duì)產(chǎn)量的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，2016，44（6）：143-148.

[6] SMOLENAARS W J，PAPARRIZOS S，WERNERS S，et al. Flood risk and adaptation strategies for soybean production systems on the flood-prone pampas under climate change[J]. Agronomy，2021，11（6）：1187.

[7] 霍治國(guó)，李茂松，王 麗，等. 氣候變暖對(duì)中國(guó)農(nóng)作物病蟲害的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2012，45（10）：1926-1934.

[8] 劉 丹，姚俊萌，余焰文，等. 氣候變暖背景下柑橘木虱在中國(guó)的適生區(qū)變化[J]. 植物保護(hù)學(xué)報(bào)，2021，48（4）：872-881.

[9] 越 昆，武榮盛，桑 婧，等. 內(nèi)蒙古春玉米氣候適宜度變化特征及評(píng)價(jià)指標(biāo)[J]. 干旱地區(qū)農(nóng)業(yè)研究，2021，39（3）：209-217.

[10] 胡春麗，陳鵬獅，王婉昭，等. 遼寧省水稻氣候適宜度變化特征及其對(duì)區(qū)域氣候變暖的響應(yīng)[J]. 江西農(nóng)業(yè)學(xué)報(bào)，2021，33（7）：13-18.

[11] 陳愛莉，龔 偉，孔 芬，等. 氣候變化對(duì)楊梅種植適宜區(qū)的影響[J]. 揚(yáng)州大學(xué)學(xué)報(bào)（農(nóng)業(yè)與生命科學(xué)版），2021，42（4）：96-102.

[12] 王彥亭，謝劍平，李志宏. 中國(guó)煙草種植區(qū)劃[M]. 北京：科學(xué)出版社，2010.

[13] 胡穎梅，徐 磊，袁 帥，等. 陜西烤煙氣候適宜性時(shí)空變化特征及風(fēng)險(xiǎn)分析[J]. 煙草科技，2021，54（7）：41-50.

[14] 胡雪瓊，徐夢(mèng)瑩，何雨芩，等. 未來氣候變化對(duì)云南烤煙種植氣候適宜性的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2016，27（4）：1241-1247.

[15] 郭瑞鴿，劉壽東，杜筱玲. 江西氣溫穩(wěn)定通過10℃初日變化及其對(duì)雙季早稻物候期的影響[J]. 中國(guó)農(nóng)業(yè)氣象，2011，32（1）：12-16.

[16] 陳 頤，周清明，楊虹琦，等. 湖南桂陽(yáng)和桑植煙區(qū)烤煙生育期氣溫及降水特征[J]. 中國(guó)煙草學(xué)報(bào)，2015，21（4）：48-55.

[17] 李 強(qiáng)，周冀衡，解 燕，等. 曲靖煙區(qū)烤煙氣候適宜性綜合評(píng)價(jià)和主成分聚類分析[J]. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版），2011， 40（5）：455-460.

[18] 彭勇剛，黃肖寒，莫益江，等. 基于層次分析法的農(nóng)業(yè)氣象災(zāi)害風(fēng)險(xiǎn)區(qū)劃指標(biāo)權(quán)重分析[J]. 氣象研究與應(yīng)用，2018，39（1）：70-72.

[19] 任國(guó)玉，郭 軍，徐銘志，等. 近50年中國(guó)地面氣候變化基本特征[J]. 氣象學(xué)報(bào)，2005，63（6）：942-956.

[20] 胡榮海，云南省煙草科學(xué)研究所，中國(guó)煙草育種研究（南方）中心. 云南煙草栽培學(xué)[M]. 北京：科學(xué)出版社，2007.

（責(zé)任編輯：張煥裕）