何興潼，羅再潛，賴(lài)自力，雷 曉，李曉勇，黃忠華，王 翔*

（1.瀘州市氣象局，四川瀘州 646000；2.瀘州市敘永縣氣象局，四川瀘州 646000；3.四川省煙草公司瀘州市公司，四川瀘州 646000）

烤煙（煙葉），一年生草本，茄科，高1 m 左右，原產(chǎn)南美洲。由于其適應(yīng)性強(qiáng)，不論山區(qū)、平原或丘陵地區(qū)均能種植。自引入我國(guó)后，很多省市均有種植[1]。近年來(lái)，我國(guó)烤煙種植業(yè)重心區(qū)位變化趨向西南煙區(qū)。

瀘州市位于長(zhǎng)江上游、四川盆地南部，地處川滇黔三省交界處，是四川盆地與云貴高原過(guò)渡地帶，分布著山地、丘陵等不同地貌，地勢(shì)呈北低南高，海拔落差可達(dá)1 500 m。瀘州市下轄三區(qū)四縣，屬亞熱帶濕潤(rùn)季風(fēng)氣候，雨熱同季，無(wú)霜期長(zhǎng)，光照適宜，非常適宜烤煙種植，北部地區(qū)（江陽(yáng)區(qū)、龍馬潭區(qū)、納溪區(qū)、瀘縣、合江縣）以丘陵地形為主，南部地區(qū)（敘永縣、古藺縣）以高山地形為主。

瀘州從20 世紀(jì)90 年代初開(kāi)始種植烤煙。隨著產(chǎn)業(yè)不斷擴(kuò)大，烤煙逐漸成為瀘州市重要經(jīng)濟(jì)作物，對(duì)瀘州市發(fā)展有重要支撐作用。其中古藺縣、敘永縣南部等適宜種植烤煙的地區(qū)，常年種植面積約1.33 萬(wàn)hm2，年均收購(gòu)煙葉2 500 萬(wàn)kg，已成為古藺、敘永兩縣財(cái)政增效、煙農(nóng)增收的主要支柱產(chǎn)業(yè)之一。對(duì)助力貧困山區(qū)人民脫貧致富、鄉(xiāng)村振興及促進(jìn)地區(qū)社會(huì)穩(wěn)定和發(fā)展發(fā)揮重要作用[2]。但每年烤煙生育期由高溫干旱、陰雨寡照等不利氣象條件造成的損失約為20%，已成為制約瀘州烤煙產(chǎn)業(yè)發(fā)展的重要因素之一。

瀘州市現(xiàn)共有198個(gè)氣象自動(dòng)觀測(cè)站，其中敘永、古藺主要烤煙區(qū)內(nèi)共有21個(gè)區(qū)域自動(dòng)氣象觀測(cè)站，自2008 年逐步開(kāi)始觀測(cè)，觀測(cè)要素多為雨量、氣溫。由于受敘永、古藺山區(qū)地形影響，氣候要素時(shí)空差異大，因此利用現(xiàn)有區(qū)域自動(dòng)氣象觀測(cè)資料開(kāi)展研究時(shí)，存在氣象資料時(shí)間序列長(zhǎng)短不一、與烤煙種植區(qū)實(shí)況氣象數(shù)據(jù)可能存在差異等問(wèn)題。

隨著社會(huì)對(duì)精細(xì)化氣象產(chǎn)品需求的不斷提升，中國(guó)氣象局公共服務(wù)中心以LAPS 系統(tǒng)[3]為基礎(chǔ)進(jìn)行本地化，通過(guò)調(diào)整權(quán)重等參數(shù)，使得分析場(chǎng)與觀測(cè)實(shí)況非常接近[4]，最終研發(fā)出一套基本覆蓋全國(guó)、基于實(shí)際地形、接近實(shí)際觀測(cè)的全國(guó)3 km×3 km 地面多要素客觀分析產(chǎn)品[5]，在氣象服務(wù)中發(fā)揮了顯著的社會(huì)效益。

基于此，本文利用烤煙種植區(qū)內(nèi)現(xiàn)有21個(gè)氣象觀測(cè)站資料，以及LAPS 多源實(shí)況格點(diǎn)資料，以期在全球變暖的大背景下[6]，分析得到瀘州烤煙區(qū)的氣候適宜種植區(qū)劃，做到趨利避害，提高煙葉產(chǎn)量，保障烤煙產(chǎn)業(yè)可持續(xù)發(fā)展。

1 資料與方法

1.1 研究區(qū)概述及數(shù)據(jù)來(lái)源

古藺、敘永兩縣地處川南邊緣，與云南、貴州接壤，境內(nèi)以山地為主，屬盆周山區(qū)，經(jīng)濟(jì)發(fā)展欠發(fā)達(dá)[7]。敘永南部、古藺作為瀘州主要烤煙種植區(qū)，主要烤煙種植品種為云煙87，蜜甜香型，截至2020年共有20 個(gè)鄉(xiāng)鎮(zhèn)、街道（敘永9 個(gè)，古藺11 個(gè)），總種植面積1.26 萬(wàn)hm2（古藺0.63 萬(wàn)hm2、敘永0.63 萬(wàn)hm2），海拔為538～1 319 m。

本研究所用數(shù)據(jù)分為2部分：1）敘永縣、古藺縣內(nèi)21個(gè)氣象觀測(cè)站2009—2020年日降水、平均2 m 氣溫資料（本文平均氣溫均代指2 m 平均氣溫）來(lái)源于瀘州市氣象局；2）中國(guó)氣象局LAPS（3 km×3 km）2009—2020 年日降水、平均2 m 氣溫多源實(shí)況融合格點(diǎn)資料（下文簡(jiǎn)稱(chēng)LAPS 格點(diǎn)數(shù)據(jù)）來(lái)源于中國(guó)氣象局，本文研究區(qū)分布如圖1所示。

圖1 瀘州烤煙種植區(qū)分布圖

1.2 研究方法

依據(jù)瀘州烤煙種植生育期研究和調(diào)查，將敘永、古藺烤煙種植區(qū)大田生育期（4—9 月）大致劃分為：團(tuán)棵期（4月上旬至5月中旬）、旺長(zhǎng)期（5月下旬至6月下旬）、成熟期（7上旬至9月下旬）。

將各烤煙種植區(qū)的氣溫、降水觀測(cè)日資料和LAPS 日資料均處理為旬資料，選取均方根誤差（RMSE）、平均誤差（ME）對(duì)LAPS 資料進(jìn)行差異檢驗(yàn)評(píng)估，方法如下：

均方根誤差（RMSE）：反映了該區(qū)域內(nèi)誤差幅度的平均狀況，表示隨機(jī)誤差，計(jì)算公式為：

平均誤差（ME）：表征當(dāng)站點(diǎn)觀測(cè)平均值高于格點(diǎn)融合產(chǎn)品的時(shí)，平均誤差大于0，反之小于0。計(jì)算公式為[8]：

式（1）、式（2）中，Xi為旬平均氣溫、旬累計(jì)降水的觀測(cè)值，Yi為旬平均氣溫、旬累計(jì)降水的LAPS格點(diǎn)值。

活動(dòng)積溫（K1）以及有效積溫（K2）計(jì)算公式分別如下：

其中Ti為第i 日平均溫度，T0定為0 ℃，T1定為10 ℃，n為生育期天數(shù)。

2 結(jié)果與分析

2.1 觀測(cè)資料與LAPS差異性檢驗(yàn)

本節(jié)選取烤煙生長(zhǎng)發(fā)育受外界氣溫、降水等因素影響最大且最為關(guān)鍵的大田生育期進(jìn)行差異檢驗(yàn)，其中部分區(qū)域站受建站時(shí)間限制，起始觀測(cè)時(shí)間段不盡相同，故將各站點(diǎn)的LAPS 數(shù)據(jù)的時(shí)間序列長(zhǎng)度與觀測(cè)時(shí)間序列長(zhǎng)度保持一致，將烤煙種植區(qū)內(nèi)各站點(diǎn)日觀測(cè)資料與LAPS 格點(diǎn)融合日數(shù)據(jù)均處理為旬資料，再應(yīng)用式（1）、式（2）進(jìn)行誤差、均方根誤差計(jì)算，結(jié)果見(jiàn)表1。

表1 瀘州烤煙種植點(diǎn)實(shí)況資料與LAPS檢驗(yàn)結(jié)果（2009—2020年）

從各種植烤煙點(diǎn)大田生育期的旬累計(jì)降水、旬平均氣溫的觀測(cè)資料與LAPS 資料的誤差和均方根誤差結(jié)果來(lái)看，總體相差不大，旬累計(jì)降水的LAPS 資料較觀測(cè)資料平均偏高0.3 mm，其中敘永水潦站點(diǎn)的誤差最大，偏高2.1 mm，其余站點(diǎn)誤差均小于1.0 mm；旬降水的均方根誤差總體平均為2.2 mm，敘永水潦站點(diǎn)較大為5.2 mm，其余站點(diǎn)為2.0 mm 左右，反映出水潦站點(diǎn)部分時(shí)段的LAPS 降水?dāng)?shù)據(jù)較觀測(cè)數(shù)據(jù)偏差幅度明顯高于其余站點(diǎn)。

旬平均氣溫的LAPS 數(shù)據(jù)較觀測(cè)數(shù)據(jù)平均偏低0.3 ℃，其中敘永合樂(lè)、古藺觀文站點(diǎn)的誤差較大，分別偏低3.3 ℃、1.9 ℃，其余站點(diǎn)均小于1.0 ℃；其均方根誤差結(jié)果也顯示總體平均為1.3 ℃，偏差幅度較大的有合樂(lè)、后山、丹桂，其均方根誤差可達(dá)2 ℃以上，其余站點(diǎn)1 ℃左右。

LAPS 格點(diǎn)融合數(shù)據(jù)與觀測(cè)數(shù)據(jù)的差異性檢驗(yàn)結(jié)果表明，旬累計(jì)降水、旬平均氣溫總體偏差不大，基本與觀測(cè)資料較為符合，其平均誤差分別為-0.3 mm、0.3 ℃，僅水潦、營(yíng)山、合樂(lè)、觀文等站點(diǎn)誤差相對(duì)較大；此外，旬累計(jì)降水與旬平均氣溫的平均誤差絕對(duì)值雖相近，但從均方根誤差的結(jié)果可看出，旬累計(jì)降水的總體誤差和偏差幅度卻大于旬平均氣溫，反映了LAPS 格點(diǎn)融合數(shù)據(jù)較觀測(cè)數(shù)據(jù)在瀘州烤煙種植區(qū)內(nèi)的旬平均氣溫比旬累計(jì)降水更接近實(shí)際且偏差幅度更小。

2.2 氣象要素年際變化

利用LAPS 格點(diǎn)數(shù)據(jù)將種植烤煙區(qū)作為整體，對(duì)2009—2020 年大田生育期的累計(jì)降水、平均氣溫、活動(dòng)積溫、有效積溫的年際變化進(jìn)行分析，結(jié)果見(jiàn)表2。

表2 大田烤煙生育期各氣象要素年變化（2009—2020年）及上等煙比例

從各氣象要素的年際平均可知，大田生育期的累計(jì)降水、平均氣溫、活動(dòng)積溫、有效積溫的年平均值分別為661.3 mm、21.1 ℃、3 787.6 ℃、1 981.7 ℃，滿(mǎn)足烤煙生育發(fā)育所需降水及光熱條件，較為適宜烤煙種植。

再?gòu)母鳉庀笠刈儺愊禂?shù)（標(biāo)準(zhǔn)差/平均值）計(jì)算結(jié)果可知，累計(jì)降水最大為0.201，其次為有效積溫（0.045），活動(dòng)積溫和平均氣溫最小，均為0.023，反映出降水量的年變化最為明顯，特別是2010、2011年的累計(jì)降水較平均值偏低220.9 mm、262.9 mm，而2014年、2020年則偏高180.9 mm、177.1 mm，年際變化幅度為26.8%～39.7%；而2013 年的平均氣溫、活動(dòng)積溫、有效積溫與其均值相差（年變化幅度）最大，分別為0.9 ℃（4.1%）、147.0 ℃（3.9%）、156.8 ℃（7.9%），也反映出烤煙種植區(qū)的累計(jì)降水年波動(dòng)幅度明顯高于平均氣溫等要素，光熱條件較降水條件年際變化穩(wěn)定。

此外，各氣象要素在2009—2014 年波動(dòng)變化明顯，但從2015年開(kāi)始，氣溫和降水均呈現(xiàn)出緩慢上升趨勢(shì)，也呈現(xiàn)出瀘州烤煙種植區(qū)在全球變暖的大背景下，具有降水、氣溫均上升的特點(diǎn)。而相關(guān)烤煙品質(zhì)數(shù)據(jù)顯示，2011 年、2013 年的上等煙比例分別為32.6%、33.6%較歷年平均（47.2%）偏低15%左右，降水偏少、氣溫偏高均不利于烤煙品質(zhì)的提升，如海拔較低的山腳河谷地帶的烤煙種植區(qū)，但對(duì)海拔相對(duì)較高的高山頂種植區(qū)有改善品質(zhì)的作用。

2.3 適宜烤煙種植氣候區(qū)劃

基于LAPS 數(shù)據(jù)對(duì)各烤煙種植區(qū)的2009—2020 年大田生育期進(jìn)行旬平均氣溫、旬累計(jì)降水、旬活動(dòng)積溫、旬有效積溫的分類(lèi)計(jì)算。并按照?qǐng)F(tuán)棵期、旺長(zhǎng)期、成熟期、大田生育期分類(lèi)統(tǒng)計(jì)，具體結(jié)果見(jiàn)表3。

表3 烤煙種植區(qū)不同生育期下各氣象要素對(duì)比（2009—2020年）

烤煙作為喜溫植物，從所需的熱量條件來(lái)看，晝夜溫差大，有利于煙草的生長(zhǎng)發(fā)育，但對(duì)煙葉的品質(zhì)來(lái)說(shuō)，成熟期溫差小反而有利。移栽后的團(tuán)棵期最適宜平均氣溫為16～18 ℃，如果氣溫低于13 ℃持續(xù)7 d以上，將導(dǎo)致早花現(xiàn)象；旺長(zhǎng)期最適宜平均氣溫為19～22 ℃，成熟期最適宜平均氣溫為22～25 ℃且至少持續(xù)30 d 以上?？緹煷筇锷诨顒?dòng)積溫要大于3 000 ℃，適宜區(qū)、最適宜區(qū)分別要求活動(dòng)積溫在3 300 ℃、3 800 ℃以上，而有效積溫為1 800～2 300 ℃。如果活動(dòng)積溫達(dá)不到烤煙正常生長(zhǎng)發(fā)育所需時(shí)，烤煙為滿(mǎn)足自己所需要的熱量總和，會(huì)自發(fā)地將生育期延長(zhǎng)，進(jìn)而直接影響煙草的產(chǎn)量與品質(zhì)[10]。

從烤煙生長(zhǎng)發(fā)育所需的水分條件來(lái)看，烤煙對(duì)水分的需求特點(diǎn)一般為“少—多—少”特點(diǎn)，即前期移栽至團(tuán)棵期，煙苗較小，生長(zhǎng)發(fā)育相對(duì)慢，需水量較少，月降雨量為60～80 mm；到了中期的旺長(zhǎng)期，烤煙各組織根系發(fā)育完善，生長(zhǎng)旺盛，也是烤煙需水最大的時(shí)間段，月降雨量以120～170 mm 為優(yōu)；最后的成熟期，月降雨量為100～130 mm 即可滿(mǎn)足烤煙的成熟落黃，抑制其部分代謝作用，從而協(xié)調(diào)煙葉的化學(xué)成分配比，對(duì)烤煙品質(zhì)提升有很大幫助[9]。

此外，海拔也一定程度上決定了烤煙生長(zhǎng)的品質(zhì)，如果烤煙種植區(qū)海拔低于550 m或者高于1 400 m，由于溫度的和地形的原因，海拔過(guò)低使得平壩地區(qū)高濕，病蟲(chóng)害易發(fā)生，造成煙草品質(zhì)低下，產(chǎn)量較低；而海拔過(guò)高，又面臨人工種植強(qiáng)度過(guò)大的問(wèn)題，不利管理，一旦遭遇旱災(zāi)，補(bǔ)救措施難以開(kāi)展。瀘州地區(qū)適宜種植烤煙的海拔為550～1 400 m，最適宜烤煙種植的海拔主要為700～1 300 m。

綜上所述，依據(jù)烤煙各生育期最適宜生長(zhǎng)的氣象指標(biāo)，并結(jié)合瀘州市烤煙種植區(qū)的種植情況，將現(xiàn)有烤煙種植區(qū)分為優(yōu)質(zhì)區(qū)、適宜區(qū)和次適宜區(qū)，結(jié)果如下。

2.3.1 優(yōu)質(zhì)區(qū)

烤煙的優(yōu)質(zhì)區(qū)有著很?chē)?yán)格的指標(biāo)，依據(jù)各生育期的不同，降水量、溫度及積溫均要滿(mǎn)足并適宜其生長(zhǎng)發(fā)育需求（見(jiàn)上文），同時(shí)不能太高或太低，否則生長(zhǎng)緩慢且影響烤煙品質(zhì)；同時(shí)海拔為700～1 300 m，受種植區(qū)所處山地的影響，優(yōu)質(zhì)烤煙區(qū)分布比較分散。

優(yōu)質(zhì)區(qū)的主要鄉(xiāng)鎮(zhèn)為敘永縣觀興、枧槽及古藺縣雙沙、皇華、石寶、大寨、德耀、馬嘶、箭竹、丹桂、茅溪等中山區(qū)鄉(xiāng)鎮(zhèn)。本區(qū)的烤煙自然生長(zhǎng)條件較好，栽種、收割時(shí)間整齊，栽種歷史較長(zhǎng)，種植和烘烤經(jīng)驗(yàn)也比較充分，為瀘州市優(yōu)質(zhì)烤煙的生產(chǎn)基地和重點(diǎn)發(fā)展地區(qū)，已經(jīng)作出了古藺優(yōu)質(zhì)煙草區(qū)“煙水”工程預(yù)案。

2.3.2 適宜區(qū)

此區(qū)域種植區(qū)的海拔也較為適宜（550～1 400 m），但由于個(gè)別氣象因子無(wú)法達(dá)到某生育期最優(yōu)的標(biāo)準(zhǔn)，如降水量明顯偏多或偏少，溫度、積溫亦這樣，即大多數(shù)氣象要素都達(dá)到優(yōu)化標(biāo)準(zhǔn)，但因?yàn)橐粋€(gè)或少數(shù)幾個(gè)要素達(dá)不到最優(yōu)，從而影響了烤煙最終品質(zhì)和產(chǎn)量。

適宜區(qū)的主要鄉(xiāng)鎮(zhèn)為敘永縣摩尼、分水、水潦、麻城、后山，以及古藺縣觀文、龍山等鄉(xiāng)鎮(zhèn)。該區(qū)雖然受個(gè)別氣象要素的影響，但總體烤煙品質(zhì)和產(chǎn)量還是較高的，農(nóng)戶(hù)在氣候正常年份種植烤煙是創(chuàng)收的，也是一個(gè)極具發(fā)展?jié)摿Φ目緹煼N植區(qū)域，可以通過(guò)人工影響的方式來(lái)達(dá)到優(yōu)質(zhì)優(yōu)產(chǎn)，如先進(jìn)種植技術(shù)、修建排灌工程等。

2.3.3 次適宜區(qū)

此區(qū)的氣候條件與適宜區(qū)相比，相差并不很大，影響其品質(zhì)及產(chǎn)量的最主要因素是地形，種植區(qū)的海拔低于550 m 或者高于1 400 m。平壩地區(qū)易積水，而高山地區(qū)又干旱，這對(duì)煙草生長(zhǎng)極為不利。但煙草在此區(qū)還是能形成產(chǎn)量，但平壩地區(qū)遇洪澇或高山地區(qū)遇干旱年份，將會(huì)嚴(yán)重影響產(chǎn)量，造成農(nóng)戶(hù)虧本。

此區(qū)包括敘永縣合樂(lè)、赤水和古藺縣金蘭等鄉(xiāng)鎮(zhèn)的平谷和高山頂區(qū)以及優(yōu)質(zhì)區(qū)山頂或山腳。此區(qū)取得豐收的關(guān)鍵為保障烤煙適宜需水，可修建排灌工程。

最后，利用ArcGIS10.0 對(duì)瀘州市適宜烤煙種植氣候區(qū)劃的結(jié)果如圖2所示。

圖2 瀘州市適宜烤煙種植區(qū)劃

3 結(jié)論與討論

本文通過(guò)對(duì)瀘州各烤煙種植區(qū)2009—2020 年的LAPS 氣溫、降水格點(diǎn)值與監(jiān)測(cè)值做了分析對(duì)比，并基于LAPS 格點(diǎn)值對(duì)各烤煙種植區(qū)做了年際變化趨勢(shì)分析和氣候適宜性種植區(qū)劃，主要結(jié)論如下：1）旬累計(jì)降水、旬平均氣溫的LAPS 格點(diǎn)值與觀測(cè)值總體偏差不大，其平均誤差分別為0.3 mm、-0.3 ℃，其中旬累計(jì)降水的總體誤差和偏差幅度略大于旬平均氣溫。2）現(xiàn)有烤煙種植區(qū)大田生育期內(nèi)的降水及光熱條件較為滿(mǎn)足烤煙生育發(fā)育所需，但從2015 年開(kāi)始，氣溫和降水均呈現(xiàn)出上升趨勢(shì)，而降水偏少、氣溫偏高均不利于烤煙品質(zhì)的提升。3）基于LAPS 格點(diǎn)數(shù)據(jù)并結(jié)合瀘州烤煙種植區(qū)情況，將瀘州現(xiàn)有烤煙種植區(qū)做了適宜氣候分區(qū)，優(yōu)質(zhì)區(qū)為敘永縣觀興、枧槽和古藺縣雙沙、皇華、石寶、大寨、德耀、馬嘶、箭竹、丹桂、茅溪等中海拔山區(qū)鄉(xiāng)鎮(zhèn)；適宜區(qū)為敘永縣摩尼、分水、水潦、麻城、后山和古藺縣觀文、龍山等鄉(xiāng)鎮(zhèn)；次適宜區(qū)為敘永縣合樂(lè)、赤水和古藺縣金蘭等鄉(xiāng)鎮(zhèn)。

已有研究表明，全球氣候變暖不僅僅伴隨氣溫的上升，還會(huì)導(dǎo)致降水量的下降[11]，也會(huì)改變烤煙各生育期時(shí)段如烤煙的移栽、收割等，優(yōu)質(zhì)烤煙種植海拔[12]。此外，烤煙品質(zhì)對(duì)氣溫、降水變化較為敏感，從全球變暖的大趨勢(shì)來(lái)看，預(yù)計(jì)未來(lái)可能會(huì)對(duì)海拔較低的山腳河谷地帶的烤煙種植區(qū)產(chǎn)生品質(zhì)下降的不利影響，但對(duì)海拔相對(duì)較高的高山頂種植區(qū)有改善品質(zhì)的作用。