陶健，劉好寶，辛玉華，劉光亮，王程棟，徐宜民，梁洪波，王樹聲

1 中國農(nóng)業(yè)科學院煙草研究所，青島 266101；2 海南省雪茄研究所，海口 571100

古巴Pinar del Rio省優(yōu)質(zhì)雪茄煙種植區(qū)主要生態(tài)因子特征研究

陶健1，劉好寶1，辛玉華2，劉光亮1，王程棟1，徐宜民1，梁洪波1，王樹聲1

1 中國農(nóng)業(yè)科學院煙草研究所，青島 266101；2 海南省雪茄研究所，?？?571100

為提供中國雪茄煙種植區(qū)域生態(tài)適宜性評價的基礎(chǔ)數(shù)據(jù)支持，研究采用NCDC-GSOD氣象數(shù)據(jù)和ISRIC-GSIF土壤信息數(shù)據(jù)，解析古巴Pinar del Rio省優(yōu)質(zhì)雪茄煙產(chǎn)區(qū)主要生態(tài)因子的動態(tài)特征與數(shù)值規(guī)律。結(jié)果表明：研究區(qū)大田期氣溫呈上升趨勢，月平均氣溫（21.21～24.41℃）、月平均最高氣溫（27.25～30.25℃）、月平均最低氣溫（15.91～19.21℃）。其中，1月份氣溫的變異性較大，月平均最低氣溫在整個生育期內(nèi)的溫差最大（6.84℃）。大田期（次年1—4月）內(nèi)的月降水總量均穩(wěn)定在90 mm以下（55.97～87.02 mm），移栽后至旺長階段1—2月份，平均56.74 mm。東、中部地區(qū)濕潤老成土和潮濕老成土為主，西部地區(qū)濕潤氧化土所占比例最高；東、中部地區(qū)土壤pH值（5.93和6.00）高于西部（5.62）；有機質(zhì)含量（18.39 g/kg和19.40 g/kg）則低于西部（22.43 g/kg）；東、中部地區(qū)的土壤質(zhì)地以沙粒為主（平均值分別為44.18%和43.89%），而西部地區(qū)則以沙粒和粘粒為主（平均值分別為42.45%和33.38%）。

氣象因子；土壤性狀；時間動態(tài)；數(shù)值特征

雪茄煙自清代傳入中國以來，消費群體及銷量不斷上升。據(jù)統(tǒng)計，2015年上半年，在中國卷煙生產(chǎn)與銷售量分別下降3.11%和2.23%的同時，中國雪茄煙產(chǎn)量增長37.20%，銷量增長82.00%，銷售額增長55.40%[1]，雪茄煙已成為中國煙草行業(yè)新的增長點。然而，當前中國的優(yōu)質(zhì)雪茄煙原料嚴重匱乏，絕大部分依賴于進口，給中國優(yōu)質(zhì)雪茄煙的產(chǎn)品開發(fā)和品牌發(fā)展帶來了極大的困難。

氣象、土壤等生態(tài)因子是決定煙葉風格特色的重要因素，就雪茄煙而言，不同雪茄結(jié)構(gòu)部位（茄衣、茄套和茄芯）的雪茄煙原料對區(qū)域生態(tài)因子的需求各異且較為苛刻[2]，即使在優(yōu)質(zhì)產(chǎn)區(qū)，也僅限于小部分煙區(qū)生產(chǎn)優(yōu)質(zhì)雪茄煙[3]。以茄衣為例，優(yōu)質(zhì)茄衣需要在高空氣溫度、高空氣濕度、多云少風的氣象條件下栽培[4]，在不同生育階段對氣象條件的需求各異[3,5]。在大田生長期，日平均溫度在20℃左右，在煙葉采收期和調(diào)制期，日平均溫度需超過22 ℃，空氣相對濕度則需要在70%以上[5]。土壤理化特性是決定優(yōu)質(zhì)雪茄煙特色形成的關(guān)鍵因素，土壤類型、質(zhì)地、有機質(zhì)含量、養(yǎng)分元素含量等因子是優(yōu)質(zhì)煙葉產(chǎn)區(qū)的關(guān)鍵生態(tài)要素[6-7]，也是雪茄煙品牌的重要生態(tài)標識[3]。土壤類型、質(zhì)地、容重等物理指標制約土壤肥力水平，影響植物根系的養(yǎng)分利用效率[8-10]。土壤養(yǎng)分元素含量則是煙葉產(chǎn)量與品質(zhì)的決定因素，土壤有機質(zhì)與土壤肥力呈正相關(guān)關(guān)系，增加其含量更能夠有效增強煙株對土壤病害的自我防控能力[11]。世界上的優(yōu)質(zhì)雪茄煙產(chǎn)地均分布在南北回歸線之間，特別是在23°N附近地區(qū)，該地區(qū)具有優(yōu)質(zhì)雪茄煙生產(chǎn)的獨特生態(tài)條件，其中尤其以古巴共和國（以下簡稱古巴）最為著名。古巴最重要的優(yōu)質(zhì)雪茄煙種植區(qū)為Pinar del Rio省的Vuelta Abajo地區(qū)。通過分析古巴Vuelta Abajo地區(qū)的氣象、土壤等主要生態(tài)因子，解析優(yōu)質(zhì)雪茄煙種植區(qū)關(guān)鍵生態(tài)因子的數(shù)值規(guī)律及其時空動態(tài)特征，把握優(yōu)質(zhì)雪茄煙生長發(fā)育對氣象、土壤等因子的生態(tài)需求，能夠有針對性地在中國選取優(yōu)質(zhì)雪茄煙種植區(qū)、對比中國雪茄煙種植區(qū)相應因子的差異，為發(fā)展相應的生產(chǎn)技術(shù)配套體系提供科學依據(jù)和數(shù)據(jù)支持。

1 研究區(qū)與數(shù)據(jù)來源

1.1 研究區(qū)簡介

古巴Pinar del Rio省位于古巴西部，是哈伯納斯公司（Habanos）雪茄煙原料的主要產(chǎn)區(qū)，也是古巴唯一種植所有種類原料的地區(qū)。該省份內(nèi)的Vuelta Abajo煙區(qū)為古巴乃至全球最著名的優(yōu)質(zhì)雪茄煙種植區(qū)，可以細分為Vuelta Abajo、San Luis、San Juan y Martínez三個片區(qū)（圖1）。整個Vuelta Abajo煙區(qū)僅有不到四分之一區(qū)域的生態(tài)條件達到哈伯納斯公司的雪茄煙最佳產(chǎn)地（古巴Vegas Finas de Primera種植園）標準[3]，其中Vuelta Abajo和San Luis片區(qū)為古巴優(yōu)質(zhì)茄衣的著名產(chǎn)區(qū)。研究基于Vuelta Abajo片區(qū)內(nèi)的氣象站點數(shù)據(jù)分析氣象因子特征，并根據(jù)Vuelta Abajo、San Luis、San Juan y Martínez 三個片區(qū)的邊界數(shù)據(jù)提取土壤信息數(shù)據(jù)分析各片區(qū)的土壤理化特性。

圖1 研究區(qū)海拔及優(yōu)質(zhì)雪茄煙種植區(qū)格局Fig. 1 Patterns of elevation and premium cigar planting areas

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

1.2.1 氣象數(shù)據(jù)

研究采用的氣象數(shù)據(jù)來源于美國國家氣候中心的全球地表日值數(shù)據(jù)集（NCDC-GSOD），選取了古巴Pinar del Rio省Vuelta Abajo片區(qū)內(nèi)氣象站點（圖1）的日值觀測數(shù)據(jù)。然后，綜合考慮數(shù)據(jù)的可用性及氣象特征分析的代表性，選取了1980—2014年的平均氣溫、最高氣溫、最低氣溫和降水總量共4項指標的日值數(shù)據(jù)處理為月平均氣溫、月平均最高氣溫、月平均最低氣溫和月降水總量。對缺失值超過10%的月統(tǒng)計數(shù)據(jù)進行剔除后，以月為時間尺度對4項氣象指標進行分析。研究區(qū)雪茄煙生育期確定為9月至次年4月，其中9—12月的氣象數(shù)據(jù)為1980—2013年平均值，次年1—4月氣象數(shù)據(jù)為1981—2014年平均值。研究區(qū)雪茄煙生育期的農(nóng)事安排時間一般為：9—10月開始育苗，45～60d后開始移栽，移栽45～70d后煙葉成熟（一般為2—3月）。

1.2.2 土壤數(shù)據(jù)

采用ISRIC-GSIF土壤數(shù)據(jù)，結(jié)合古巴Pinar del Rio省Vuelta Abajo煙區(qū)內(nèi)三大片區(qū)的邊界數(shù)據(jù)，提取土壤理化信息。該數(shù)據(jù)基于全球共計約100000個土壤剖面數(shù)據(jù)，采用線性-克里金方法結(jié)合樣條函數(shù)，對全球土壤理化性狀數(shù)據(jù)進行制圖。該數(shù)據(jù)的空間分辨率為1 km×1 km，已被廣泛應用于全球生態(tài)學研究中[12-13]。綜合考慮土壤對雪茄煙種植的影響和對人類施肥與農(nóng)事操作響應的敏感性，研究中選取土壤亞綱類型（USDA分類標準）、土壤pH、土壤有機質(zhì)含量、土壤質(zhì)地（包括沙粒含量、粉粒含量和粘粒含量）共6項重要而較為穩(wěn)定的指標進行分析，其中土壤pH、有機質(zhì)含量和質(zhì)地均為土壤深度為0～22.5 cm平均值。

2 結(jié)果與分析

2.1 氣象因子的時間動態(tài)特征

2.1.1 空氣溫度

由圖2可知，研究區(qū)空氣溫度在育苗期（9—12月）內(nèi)呈下降趨勢，移栽期（次年1月）的氣溫最低，次年2—4月的氣溫逐漸升高。月平均氣溫（圖2a）、月平均最高氣溫（圖2b）、月平均最低氣溫（圖2c）在全生育期內(nèi)的平均值分別為23.41℃、29.19℃、18.62℃，其中育苗期平均值分別為24.21℃、29.80℃、20.00℃，大田期平均值分別為22.61℃、28.58℃、17.25℃，成熟期平均值分別為22.41℃、28.41℃、16.94℃，在整個生育期內(nèi)的溫差分別為5.05℃（21.21～26.26℃）、4.66℃（27.25～31.91℃）、6.84℃（15.91～22.75℃）。綜合對比箱形圖和各月份的標準偏差可知，月平均氣溫（圖2a）在大田期內(nèi)1月、3月和育苗期內(nèi)11月的變異性較大，標準偏差分別為1.26℃、1.17℃、1.15℃，育苗期初期9—10月及末期12月平均氣溫的變異性最??；月平均最高氣溫（圖2b）在大田期（次年1—4月）內(nèi)的變異性較大，尤其是次年4月和1月份，標準偏差分別為1.28℃、1.07℃；月平均最低氣溫（圖2c）在11月—次年4月內(nèi)的變異性較大，其中次年1月平均最低氣溫的變異性最大，標準偏差為1.71℃。

圖2 研究區(qū)氣溫月平均值變化趨勢Fig. 2 Trends of monthly average temperatures

2.1.2 降水總量

研究區(qū)月降水總量在全生育期內(nèi)的總量為737.93 mm，其中育苗期總量為455.11 mm，大田期總量為282.82 mm，成熟期總量為142.98 mm。由圖3可知，研究區(qū)育苗期初期（9—10月）的月降水總量最高，均在150 mm以上，末期（12月）降水最低（1980—2013年平均值：31.51 mm），大田期內(nèi)（次年1—4月）的月降水總量穩(wěn)定在90 mm以下。育苗期9—10月降水總量的變異性最大，標準偏差分別為121.60 mm、73.55 mm，11月到次年2月降水總量的變異性較小。

圖3 研究區(qū)月降水總量變化趨勢Fig. 3 Trend of monthly total precipitation

2.2 土壤因子的空間格局

2.2.1 土壤類型

由圖4a可知，研究區(qū)東部、中部地區(qū)以潮濕與濕潤土壤為主，而西部則以半濕潤土壤為主。東、中部地區(qū)位于奧爾加諾斯山南部，為地勢平緩的山麓平原，空氣濕度高，河流密集而短促，河谷農(nóng)業(yè)發(fā)達，土綱以老成土為主；西部為奧爾加諾斯山西南末端，海拔較高而空氣濕度較低，土綱以氧化土和軟土為主。

圖4 研究區(qū)土壤分類格局（a）及各亞綱面積比例（b）Fig.4 Pattern of soil class (a) and area percentage of each suborder (b)

由圖4b可知，研究區(qū)土壤類型主要以濕潤老成土和潮濕老成土為主，兩者面積占研究區(qū)總面積的51.22%（濕潤老成土占30.20%，潮濕老成土占21.02%）。對比三個片區(qū)來看，Vuelta Abajo片區(qū)中的濕潤老成土所占比例最高（46.06%），San Luis片區(qū)則以潮濕老成土（占37.56%）和濕潤老成土（占29.85%）為主，San Juan y Martínez片區(qū)中的濕潤氧化土所占比例最高（30.58%）。

2.2.2 土壤pH

由圖5a可知，研究區(qū)土壤pH呈西北部低而東部、中部高的空間格局，中部土壤pH最高。由圖5b可知，研究區(qū)土壤pH的平均值為5.89（范圍：5.30～6.30），其中pH在5.80～6.10間的土壤面積占研究區(qū)總面積的63.49%。對比三個片區(qū)來看，Vuelta Abajo和San Luis片區(qū)的土壤pH平均值分別為5.93和6.00，主要在5.80～6.10范圍內(nèi)（Vuelta Abajo占84.62%，San Luis占 76.69%），San Juan y Martínez片區(qū)的土壤pH平均值為5.62，主要在5.40～5.60范圍內(nèi)（占63.03%）。

圖5 研究區(qū)土壤pH格局（a）及面積比例（b）Fig. 5 Pattern (a) and area percentage (b) of soil pH

2.2.3 土壤有機質(zhì)含量

圖6 研究區(qū)土壤有機質(zhì)含量格局（a）及面積比例（b）Fig. 6 Pattern (a) and area percentage (b) of soil organic matter content

由圖6a可知，研究區(qū)土壤有機質(zhì)含量呈東北部、中部低而西部高的空間格局，西部奧爾加諾斯山的西南末端土壤有機質(zhì)含量較高。由圖6b可知，研究區(qū)土壤有機質(zhì)含量的平均值為20.40 g/kg（范圍：10.67～31.33 g/kg），其中有機質(zhì)含量在14.00～26.00 g/kg間的土壤面積占研究區(qū)總面積的92.02%。對比三個片區(qū)來看，Vuelta Abajo和San Luis片區(qū)的土壤有機質(zhì)含量平均值分別為18.39 g/kg和19.40 g/kg，主要在15.00～21.00 g/kg范圍內(nèi)（Vuelta Abajo占63.84%，San Luis占 63.97%），San Juan y Martínez片區(qū)的土壤有機質(zhì)含量平均值為22.43 g/kg，主要在20.00～26.00 g/kg范圍內(nèi)（占84.41%）。

2.2.4 土壤質(zhì)地

研究區(qū)土壤沙粒含量呈南部高而北部、西北部低的空間格局（圖7a），土壤粉粒含量呈北部高而西北部低的空間格局（圖7c），土壤粘粒含量呈北部、西北部高而南部低的空間格局（圖7e）。北部、西北部地區(qū)分別為奧爾加諾斯山的山麓和西南末端，土壤中沙粒含量較低，而粘粒含量較高；南部則為河谷平原區(qū)域，土壤中沙粒含量高，而土壤粉粒和粘粒含量低。

研究區(qū)土壤沙粒含量的平均值為43.52%（范圍：32.67～50.67%）,土壤粉粒含量的平均值為27.09%（范圍：18.00～33.33%）,土壤粘粒含量的平均值為29.40%（范圍：24.00～41.00%）。其中，土壤沙粒含量主要在40.00～47.00%間（占76.73%），土壤粉粒含量主要在26.00～31.00%間（占73.00%），土壤粘粒含量主要在25.00～28.00%間（占62.92%）。

對比來看，Vuelta Abajo和San Luis片區(qū)的土壤質(zhì)地相近：兩者的土壤沙粒含量（圖7b）平均值分別為44.18%和43.89%，主要在42.00～47.00%范 圍 內(nèi)（Vuelta Abajo占70.31%，San Luis占69.85%）；土壤粉粒含量（圖7d）平均值分別為28.44%和28.55%，主要在26.00～31.00%范圍內(nèi)（Vuelta Abajo占 88.74%，San Luis占 91.18%）；土壤粘粒含量（圖7f）平均值分別為27.42%和27.52%，主要在25.00～28.00%范圍內(nèi)（Vuelta Abajo占78.67%，San Luis占80.39%）。San Juan y Martínez片區(qū)的土壤沙粒含量平均值為42.45%，主要在40.00～43.00%范圍內(nèi)（占54.55%）；粉粒含量平均值為24.18%，主要在19.00～22.00%范圍內(nèi)（占50.62%）；粘粒含量平均值為33.38%，主要在35.00～39.00%范圍內(nèi)（占55.59%）。

圖7 研究區(qū)土壤質(zhì)地格局及其面積比例Fig. 7 Patterns and area percentages of soil texture factors

3 討論

3.1 生態(tài)因子特征

研究區(qū)大田期氣溫均呈上升趨勢（月平均氣溫：21.21～24.41℃、月平均最高氣溫：27.25～30.25℃、月平均最低氣溫：15.91～19.21℃），大田期降水則基本穩(wěn)定在90 mm以下（55.97～87.02 mm）。盡管Pinar del Rio省位于古巴最濕潤地區(qū)，年降水量在1600 mm以上，但大田期內(nèi)平均僅70.71 mm，尤其是移栽后至旺長階段1—2月份，平均56.74 mm，研究表明在該階段內(nèi)月降水量在50.80 mm左右符合雪茄煙生長最佳條件[3]。對土壤亞綱類型、pH、有機質(zhì)含量、質(zhì)地（包括土壤沙粒、粉粒、粘粒含量）的空間格局和數(shù)值規(guī)律進行分析發(fā)現(xiàn)東部、中部的Vuelta Abajo、San Luis片區(qū)與西部的San Juan y Martínez片區(qū)間的土壤理化特性有明顯差異。Vuelta Abajo和San Luis片區(qū)屬于著名的優(yōu)質(zhì)茄衣產(chǎn)區(qū)，而San Juan y Martínez片區(qū)則以種植優(yōu)質(zhì)茄芯為主。不同雪茄煙原料的生態(tài)條件需求和種植管理方式存在明顯的差異，因此未來研究需要對優(yōu)質(zhì)茄衣產(chǎn)區(qū)（如：Vuelta Abajo）和優(yōu)質(zhì)茄芯產(chǎn)區(qū)（如：Semi Vuelta）的關(guān)鍵氣象、土壤因子進行對比，解析不同雪茄煙原料的生態(tài)條件差異，揭示不同雪茄煙原料特色風格的生態(tài)成因。

中國雪茄煙傳統(tǒng)種植區(qū)以四川省什邡市以及浙江省桐鄉(xiāng)市最為著名，兩者的地理分布距23°N較遠，因此種植時間與古巴不同。以什邡煙區(qū)為例，一般12月上旬育苗，3月上中旬移栽，5月中下旬開始采收晾制，大田期4—5月日平均氣溫20℃作用，低于古巴雪茄煙區(qū)，降水總量150.90 mm，高于古巴雪茄煙區(qū)[14,5]。什邡煙區(qū)土壤為沙壤土，土壤pH大多在5.5～6.5之間，有機質(zhì)含量為22.8～52.5 g/kg[5]，有機質(zhì)含量高于古巴雪茄煙區(qū)。此外，中國海南省、雷州半島及臺灣省南部地區(qū)屬于23°N以南地區(qū)，并與古巴同屬于海島和半島的海洋性氣候，土壤理化特性類型也極為相近。因此，下一步將基于本研究結(jié)果，依據(jù)Vuelta Abajo煙區(qū)數(shù)據(jù)，對上述地區(qū)的氣象、土壤因子進行對比評價。

3.2 特征分析的意義與不確定性

解析古巴Pinar del Rio省Vuelta Abajo煙區(qū)氣溫、降水因子的月尺度動態(tài)和土壤主要理化指標的數(shù)值特征，能夠為中國雪茄煙種植區(qū)的選址和種植、管理措施的制定提供數(shù)據(jù)基礎(chǔ)?；谘芯拷Y(jié)果，結(jié)合中國傳統(tǒng)雪茄煙種植區(qū)和中國熱帶農(nóng)業(yè)區(qū)的氣象和土壤理化信息，采用數(shù)據(jù)統(tǒng)計分析或物種分布模型，選擇與Vuelta Abajo煙區(qū)生態(tài)因子最為相似的區(qū)域作為中國雪茄煙種植的重點發(fā)展區(qū)域；通過捕捉與Vuelta Abajo煙區(qū)大田期氣象因子動態(tài)的相關(guān)程度最高的時段，調(diào)整中國雪茄煙種植區(qū)的種植時間節(jié)律；通過對比中國雪茄煙種植區(qū)生態(tài)因子與Vuelta Abajo煙區(qū)的差異及相應時間點，改進管理措施以提升本地生態(tài)條件。

研究中，氣象數(shù)據(jù)和土壤信息數(shù)據(jù)均來源于全球數(shù)據(jù)集，因此在研究區(qū)內(nèi)應用時受區(qū)域內(nèi)數(shù)據(jù)量及數(shù)據(jù)來源站點或采樣點空間分布的限制，數(shù)據(jù)精度有一定的不確定性。在氣象指標選取時，NCDC-GSOD數(shù)據(jù)中的指標包括氣溫、降水、風速、露點溫度、海平面高度、大氣壓力，缺乏雪茄煙種植的重要指標空氣濕度，因此研究中未對該指標進行分析。土壤指標的選取則主要考慮到對雪茄煙種植的影響和對施肥與農(nóng)事操作響應的敏感性。未來研究需要結(jié)合多源氣象數(shù)據(jù)確保數(shù)據(jù)精度和指標完善性，如NCEP Reanalysis，并獲取完整的土壤理化指標數(shù)據(jù)，對雪茄煙種植的關(guān)鍵生態(tài)因子進行全面分析，如：太陽輻射、空氣濕度、土壤容重等。此外，研究中以月為時間尺度對雪茄煙生育期氣象因子進行動態(tài)分析，然而就煙草作物種植而言，月尺度的時間分辨率較低，未來研究需要進一步細化至旬或按照研究區(qū)雪茄煙種植的不同生育期進行時間劃分，深入細化優(yōu)質(zhì)雪茄煙葉生長發(fā)育的氣象要素特征及不同原料產(chǎn)地間的差異。

3.3 遮蔭效應

優(yōu)質(zhì)茄衣產(chǎn)區(qū)（如古巴、尼加拉瓜、多米尼加等）一般采用遮光棚（Cheese-cloth/muslin）對雪茄煙田進行遮光處理，實現(xiàn)優(yōu)質(zhì)茄衣的最佳生長環(huán)境[15-16]，而茄套和茄芯一般在日光條件下生長，不需要進行遮光處理。采用光度測量法對遮光棚效果的研究發(fā)現(xiàn)，棚內(nèi)的光照總量比棚外減少約1/3，但散射光基本不變，棚內(nèi)散射光所占光照總量的比率（0.33～0.72）遠高于棚外（0.27～0.52）。遮光棚對日平均氣溫的影響較小，遮光棚內(nèi)氣溫比棚外低約0.08℃（60d平均值），但是遮光棚內(nèi)外氣溫日較差的差異較大。氣溫日較差的差異易導致作物葉片呼吸和蒸騰作用的日夜間差異[17-19]，進而影響煙葉同化物質(zhì)積累[20-21]。遮光棚對日間空氣相對濕度的影響則遠遠高于氣溫，并且作物葉片蒸騰作用的日夜間差異能夠增強遮光棚內(nèi)外空氣濕度的差異[18,22]。因此，未來研究需要進一步就遮光棚對氣象條件的調(diào)控作用及其生態(tài)效應進行深入解析。

4 結(jié)論

（1）研究采用NCDC-GSOD氣象數(shù)據(jù)解析古巴Vuelta Abajo煙區(qū)氣溫、降水因子的時間動態(tài)特征，為中國雪茄煙種植區(qū)域評價提供基礎(chǔ)氣象信息。研究區(qū)平均氣溫、最高氣溫、最低氣溫三者在整個生育期內(nèi)的時間動態(tài)特征基本一致，大田期氣溫均呈一致的上升趨勢。其中1月份氣溫的變異性較大，最低氣溫在整個生育期內(nèi)的溫差最大（6.84℃），然后為平均氣溫（5.05℃）、最高氣溫（4.66℃）。雪茄煙生育期內(nèi)的降水量主要集中在育苗期初期（9—10月），并且該時段內(nèi)月降水總量的變異性較大，末期（12月）降水最低，大田期（次年1—4月）內(nèi)的月降水總量均穩(wěn)定在90 mm以下（55.97～87.02 mm），尤其是移栽后至旺長階段1—2月份，平均56.74 mm。

（2）研究區(qū)主要為地勢平緩的山麓平原，河谷農(nóng)業(yè)發(fā)達，土壤類型以濕潤老成土和潮濕老成土為主，不同雪茄煙原料種植區(qū)域的土壤理化性狀呈明顯的差異。東、中部地區(qū)（優(yōu)質(zhì)茄衣產(chǎn)區(qū)）與西部地區(qū)（優(yōu)質(zhì)茄芯產(chǎn)區(qū)）的土壤理化性狀有明顯的空間差異：東、中部地區(qū)濕潤老成土和潮濕老成土所占比例最高，西部地區(qū)濕潤氧化土所占比例最高；東、中部地區(qū)的土壤pH值（5.93和6.00）高于西部（5.62），有機質(zhì)含量（18.39 g/kg和19.40 g/kg）低于西部（22.43 g/kg），東、中部地區(qū)的土壤質(zhì)地中沙粒含量最高（平均值分別為44.18%和43.89%），土壤粉粒和粘粒含量相近（粉粒含量平均值分別為28.44%和28.55%，粘粒含量平均值分別為27.42%和27.52%），西部地區(qū)的土壤沙粒含量最高（平均值為42.45%），其次為土壤粘粒含量（平均值為33.38%），土壤粉粒含量最低（平均值為24.18%）。

（3）研究能為中國雪茄煙種植區(qū)生態(tài)適宜性評價提供數(shù)據(jù)基礎(chǔ)：基于研究結(jié)果，結(jié)合中國植煙區(qū)域生態(tài)數(shù)據(jù)，能夠準確評估中國種植雪茄煙的生態(tài)適宜區(qū)，并精確捕捉適宜的種植時間節(jié)律，量化中國雪茄煙種植區(qū)的生態(tài)差異及相應時間點，在相應時間點改進管理措施以彌補差異。未來將就上述內(nèi)容開展生態(tài)適宜性評價、移栽期及種植節(jié)律時間提取和生態(tài)資源匹配評估三方面的相關(guān)研究。

[1]王文華. 煙草行業(yè)2015年上半年經(jīng)濟運行穩(wěn)步提升[EB/OL].東方煙草網(wǎng), (2015-7-15), [2016-5-17]. http://www.eastobacco.com/pub/web/zxbk/dyzxzx/gjjwj/201507/t20150715_374569.html.WANG Wenhua. A steadily promotion of tobacco industry economy in the first half Year of 2015[EB/OL].East tobacco，(2015-7-15), [2016-5-17]. http://www.eastobacco. com/pub/web/zxbk/dyzxzx/gjjwj/201507/t20150715_374569.html.

[2]秦艷青, 李愛軍, 范靜苑, 等. 優(yōu)質(zhì)雪茄茄衣生產(chǎn)技術(shù)探討[J]. 江西農(nóng)業(yè)學報, 2012, 07: 101-103.QIN Yanqing, LI Aijun, FAN Jingyuan, et al. Discussion on production technology of high-quality cigar wrapper[J].Acta Agriculturae Jiangxi, 2012, 07: 101-103.

[3]Espino M. Cuban cigar tobacco. Why Cuban cigars are the world’s best[M]. TFH Publications, 1996: 23-36.

[4]王浩雅, 左興俊, 孫福山, 等. 雪茄煙外包葉的研究進展[J]. 中國煙草科學, 2009, 05: 71-76.WANG Haoya, ZUO Xingjun, SUN Fushan, et al. Advance in cigar wrapper tobacco[J].Chinese Tobacco Science, 2009,05: 71-76.

[5]李愛軍, 秦艷青, 范靜苑, 等. 四川省發(fā)展雪茄煙葉生產(chǎn)的優(yōu)勢分析及對策[J]. 安徽農(nóng)業(yè)科學, 2013, (07): 3188-3189, 3228.LI Aijun, QIN Yanqing, FAN Jingyuan, et al. Advantage analysis and counter measure research on cigar tobacco production in Sichuan Province[J]. Journal of Anhui Agricultural Sciences, 2013, (07): 3188-3189, 3228.

[6]Davis D L, Nielsen M T. Tobacco: production, chemistry and technology[M]. Blackwell Science Ltd, 1999: 10-30.

[7]Delibacak S, OngunA R, Ekren S. Influence of soil properties on yield and quality of tobacco plant in Akhisar region of Turkey[J]. Eurasian Journal of Soil Science, 2014,3(4): 286–292.

[8]郝葳, 田孝華. 優(yōu)質(zhì)煙區(qū)土壤物理性狀分析與研究[J].煙草科技, 1996(05): 34-35.HAO Wei, TIAN Xiaohua. Research on soil physical properties in premium cigar planting area[J]. Tobacco Science & Technology, 1996(05): 34-35.

[9]張海林, 秦耀東, 朱文珊. 耕作措施對土壤物理性狀的影響[J]. 土壤, 2003, 02: 140-144.ZHANG Hailin, QIN Yaodong, ZHU Wenshan. Effects of tillage on soil physical properties[J]. Soil, 2003, 02: 140-144.

[10]Olivet Y E, Sanchez-Giron V, Hernanz J L. Reduced tillage for tobacco (Nicotiana tabacum L.) production in East Cuba[J]. Soil Physical Properties and Crop Yield, 2014,12(3): 12.

[11]Yulianti T, Hidayah N. Soil amendments with organic matter for the control of hollow stalk (Pectobacterium carotovorum subsp. carotovorum) of Besuki cigar tobacco[J]. Agrivita,2011,33(2): 199.

[12]Arrouays D, Grundy M G,Hartemink A E, et al. Chapter Three-Global soil map: toward a fi ne-resolution global grid of soil properties, in Advances in Agronomy[M]. Academic Press, 2014: 93-134.

[13]Hengl T, de Jesus J M, MacMillan R A, et al. Soil grids 1km-global soil information based on automated mapping[J]. PLoS ONE, 2014,9(8): e105992.

[14]李愛軍, 秦艷青, 代惠娟, 等. 國產(chǎn)雪茄煙葉科學發(fā)展芻議[J]. 中國煙草學報, 2012, (01): 112-114.LI Aijun, QIN Yanqing, DAI Huijuan, et al. On scientific development of China’s cigar leaf[J]. Acta Tabacaria Sinica,2012, (01): 112-114.

[15]Hasselbring H. The effect of shading on the transpiration and assimilation of the tobacco plant in Cuba[J]. Botanical Gazette, 1914, 57(4): 257-286.

[16]時向東, 汪文杰, 王衛(wèi)武, 等.遮蔭下氮肥用量對雪茄外包皮煙葉光合特性的調(diào)控效應[J]. 植物營養(yǎng)與肥料學報,2007, 02: 299-304.SHI Xiangdong, WANG Wenjie, WANG Weiwu, et al.Response of photosynthetic characteristics in leaves of cigar wrapper tobacco to nitrogen application under shady condition[J]. Plant Nutrition and Fertilizer Science, 2007,02: 299-304.

[17]Gates D M. Transpiration and leaf temperature[J]. Annual Review of Plant Physiology, 1968,19(1): 211-238.

[18]Evans J R, Von Caemmerer S. Temperature response of carbon isotope discrimination and mesophyll conductance in tobacco[J]. Plant, Cell & Environment, 2013, 36(4): 745-756.

[19]HU Jing, YANG Qiuyun, HUANG Wei, et al. Effects of temperature on leaf hydraulic architecture of tobacco plants[J]. Planta, 2014,240(3): 489-496.

[20]Neales T F, Incoll L. The control of leaf photosynthesis rate by the level of assimilate concentration in the leaf: a review of the hypothesis[J]. The Botanical Review, 1968, 34(2):107-125.

[21]Peterson R B, Zelitch I. Relationship between net CO2assimilation and dry weight accumulation in field-grown tobacco[J]. Plant Physiology, 1982, 70(3): 677-685.

[22]Jauregui I, Aroca R, Garnica M, et al. Nitrogen assimilation and transpiration: key processes conditioning responsiveness of wheat to elevated CO2and temperature[J]. Physiologia Plantarum, 2015, 155(3): 338-354.

Research on characteristics of major ecological factors in growing areas for premium cigar tobacco in Pinar del Rio province of Cuba

TAO Jian1，LIU Haobao1，XIN Yuhua2，LIU Guangliang1，WANG Chengdong1，XU Yimin1，LIANG Hongbo1，WANG Shusheng1
1 Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao 266101, China；2 Hainan Cigar Research Institute, Haikou 571100, China

Meteorological records from NCDC-GSOD and soil property dataset obtained from ISRIC-GSIF were used to investigate temporal tendency of meteorological factors and spatial pattern, numerical feature of soil property in premium cigar planting area of Pinar del Rio Province, Cuba. Results showed that: during fi eld growth period after transplantation, temperature showed an increasing trend (monthly average temperature: 21.21～24.41 ℃ , monthly average highest temperature: 27.25～30.25 ℃ , monthly average lowest temperature: 15.91～19.21℃ ). Temperature changed greatly in January with averaged lowest temperature having the largest temperature di ff erence of 6.84 ℃. At the same time, the monthly precipitation remained below 90 mm (from 55.97 mm to 87.02mm), and at average 56.74 mm between January and February. The spatial patterns of soil class and properties showed distinctive diversity between di ff erent planting areas of di ff erent cigar parts, i.e. wrapper planting areas (WPAs) in the eastern (Vuelta Abajo) and central (San Luis) parts of the study area, fi ller planting area (FPA) in the western (San Juan y Martínez) part of the study area. Soil pH values in WPAs (mean values were 5.93 and 6.00 respectively) were greater than that in FPA (mean value: 5.62) and content of organic substance was lower in WPAs(18.39g/kg and 9.40g/kg respectively) than that in FPA (22.43g/kg). Sand was the main type of soil texture in WPAs, with contents of 44.18%(Vuelta Abajo) and 43.89% (San Luis) respectively. In FPA, soil sand and clay were the main types, with contents of 42.45 % and 33.38 %respectively.

meteorological factors; soil property; temporal tendency; numerical feature

陶健，劉好寶，辛玉華，等. 古巴Pinar del Rio省優(yōu)質(zhì)雪茄煙種植區(qū)主要生態(tài)因子特征研究[J]. 中國煙草學報，2016，22(4)

中國農(nóng)業(yè)科學院煙草研究所青年科學基金資助項目（No. 2015A02）

陶健（1983—），博士，助理研究員，主要研究方向為農(nóng)田生態(tài)學，Tel：0532-66715598，Email：taojiancaas@163.com

2016-01-15

:TAO Jian，LIU Haobao，XIN Yuhua, et al. Research on characteristics of major ecological factors in growing areas for premium cigar tobacco in Pinar del Rio province of Cuba [J]. Acta Tabacaria Sinica, 2016,22(4)