基于GIS的普定縣茶樹種植氣候適宜性區劃

楊 益，陳貞宏，吳 龍，朱曦嶸，張嬌艷

(1.蘭州大學，甘肅 蘭州 730000;2.貴州省普定縣氣象局，貴州 普定 562100;3.貴州省安順市氣象局，貴州 安順 561000;4.貴州省冊亨縣氣象局，貴州 冊亨 552200;5.貴州省氣候中心，貴州 貴陽 550002)

1 引言

茶樹是一種多年生長的常綠植物，對氣候條件有一定的要求，主要分布在亞熱帶和熱帶地區，但是茶樹經過人工栽培后，它的適應范圍已遠遠超過原生長地條件。普定朵貝茶是貴州省著名的歷史名茶之一，近年來，普定縣大力實施農業產業化工程，以茶產業為主導，以規模化建設、產業化經營、品牌化建設打造朵貝品牌。為了充分挖掘氣候資源潛力，對普定縣茶葉種植氣候進行區劃，有著重要的現實意義。國內學者在茶樹氣候適宜性方面有些研究［1－3］，但在茶樹區劃方面研究甚少，目前對茶樹氣候適宜性研究［4－5］的多以數理統計圖表為主或是以行政區劃為基本評價單元，精細化程度不夠。筆者以普定縣為研究區域，基于GIS技術［6］及氣候資源的小網格分析方法［7］，綜合分析氣候條件，對該縣茶葉適宜種植區進行氣候區劃細化研究，旨在因地制宜地發展茶葉生產，將為該縣茶葉種植尤其是名茶種植的合理布局與規劃提供科學依據。

2 資料來源

2.1 研究區概況

普定縣位于貴州省中部偏西，東西長51.4 km，南北寬40 km，總面積1 085.38 km2，轄7個鄉，4個鎮，327個行政村，屬于亞熱帶季風濕潤氣候，全年氣候溫和，冬無嚴寒，夏無酷暑，春干秋涼，無霜期長，雨量充沛，日照少，輻射能量低。

2.2 氣候資料

氣象數據為安順市6個氣象站1961—2010年氣溫、降雨和日照數據。經過統計得出近50 a年平均氣溫、極端最低氣溫、≥10℃活動積溫、年平均降雨量、年平均日照時數資料;地理信息數據主要應用安順市普定縣100 m×100 m的DEM數據。

3 研究方法

3.1 茶樹生長發育的氣候條件

茶樹是一種喜溫、喜濕、喜漫射光的多年生長的亞熱帶常綠植物［8］，具有不耐高溫、嚴寒、干燥天氣的特性。因此需要在一定的氣候條件下才能正常生長，否則將受到不同程度的影響甚至危害，光照、溫度、水分是影響茶葉生長和品質的最主要氣象因子。貴州茶葉生長氣候條件一般要求年平均氣溫13℃以上，≥10℃的年積溫4 000℃以上，年降水量1 300 mm以上，冬季氣溫＜－13℃會出現嚴重凍害，夏季出現35℃以上高溫，并持續幾天時，生長會受到抑制，茶葉品質也將受到嚴重影響。因此，培育優質名茶不僅要有好的品種，適宜的土壤，而且要有適宜的氣候條件。貴州降水豐富，氣溫偏低，多云霧［9］，多陰天，少日照，太陽輻射強度小，散射光多，有利于茶葉鮮葉氨基酸含量的增加，提高茶葉的品質。

普定縣年平均氣溫15.1℃，年平均日照時數1 175.3 h，無霜期301 d，年平均降水1 348.4 mm，屬全省3大降雨中心地區之一。以上氣候條件對茶葉的生長發育比較有利。

3.2 區劃因子的選擇及適宜性指標的確定

一般情況下，茶樹生長要求年平均氣溫≥13℃、年內≥10℃的活動積溫在4 000℃/a以上和年降水量達1 300 mm以上，以及空氣相對濕度在70%以上。本文根據金志鳳等人的研究結果［10－12］，結合貴州茶樹生長發育的氣候條件綜合分析，確定以年平均氣溫(T年)、歷史極端最低氣溫(Tmin)、≥10℃活動積溫(∑T≥10)、年平均降雨量(R年)、年平均日照時數(S年)5個氣象因子為主要氣候區劃指標，其權重分別為 0.3、0.2、0.3、0.1、0.1，即以年平均氣溫、≥10℃活動積溫為主導因子對普定縣茶樹適宜性種植進行氣候區劃，建立茶樹種植適宜性區劃指標(表1)。

表1 茶葉適宜種植的評價要素及指標

3.3 茶樹適宜種植區的空間分析模型的建立

3.3.1 地理信息數據的編輯與處理 應用普定縣100 m×100 m的DEM數據，在地理信息系統Arc-GIS9.3平臺下，進行數據的導入、投影、配準等圖層處理，最終生成普定全縣100 m×100 m的海拔高度、經度和緯度的柵格圖層。

3.3.2 區劃指標小網格推算模型的建立 氣候資源由于緯度、海陸分布、地形地貌以及下墊面性質的影響存在著明顯的空間差異。普定縣已建成兩要素區域自動觀測站，然而區域自動站建站時間不長，并不能代表氣候資料，若僅靠普定縣1個站的氣候資料只能反映該站附近區域內的氣候概況，不能代表整個縣實際情況。為解決這一難題，通過建立整個安順市(包括6個縣市)山地氣候的空間分析模型，以此來推算氣候區劃指標在無測站地區的分布狀況。區劃因子與地理因子的關系模型可表示為［13］:

式中Y為氣候區劃因子(如年平均氣溫、歷年極端最低氣溫等);φ、λ、h分別表示經度、緯度、海拔高度。

應用安順市6個縣市1961—2010年50 a氣象觀測資料，統計分析得出年平均氣溫、歷年極端最低氣溫、≥10℃活動積溫、年平均降雨量、年平均日照時數，根據多元回歸分析方法［14］，求出5個區劃指標的氣候學方程(見表2)。

表2 茶樹種植指標各要素的空間模型

從表2看出，各方程的復相關系數在0.820 4～0.994 6之間，從回歸效果看，降雨量推算不如溫度的好，但都通過了α=0.05的顯著性檢驗［14］，表明可以用推算的方程進行區劃指標的小網格推算。

3.3.3 種植氣候區劃專題圖的制作 在地理信息系統ArcGIS9.3平臺中，根據生成的經度、緯度和海拔高度柵格圖層，利用表2中的空間推算模型，分別推算出年平均氣溫、歷年極端最低氣溫、≥10℃活動積溫、年平均降雨量、年平均日照時數柵格圖層數據，根據表1茶樹種植適宜性區劃指標，將上述5個柵格圖層數據，按照最適宜、適宜、次適宜、不適宜的標準同化(此處同化采用ArcGIS9.3中的重分類Reclassify方法而得)成一定的數值(如最適宜區范圍內的數值統一同化為20，適宜區為15，次適宜區為10，不適宜區為5)，最后將5個同化后的區劃指標權重系數分別為 0.3、0.2、0.3、0.1、0.1，利用ArcGIS9.3的柵格疊加工具(Raster Calculator)將5項指標進行權重疊加，再按照總分值的大小根據ArcGIS9.3自然斷點分級法［15］來劃分普定縣茶樹種植的最適宜、適宜、次適宜和不適宜氣候區，結果見圖1所示。

圖1 普定縣茶樹種植氣候適宜區劃

4 結果分析與區劃評述

由圖1可知，普定縣茶樹種植的氣候適宜性自東北向西南呈明顯遞增趨勢，充分反映了氣候條件及地形地貌因素的影響，在其他條件相同的情況下，氣溫隨海拔高度升高而降低，因此在普定縣北部地勢較高鄉鎮不適宜種植茶樹。具體分析如下:

①最適宜區。最適宜區域分布最廣，占該縣整個國土面積的46%，主要分布在該縣中部偏南，城關鎮南側、化處、馬官、馬場等地，尤其是化處鎮，盛產貴州歷史名茶朵貝茶，這與區劃結果相吻合;朵貝茶在古代多次作為貢品進貢給皇帝，近年來也曾多次榮獲中國茶葉博覽會銀獎和金獎，因此，該區是茶樹種植的重點發展區域。

②適宜區。適宜區主要集中在城關北側、龍場、白巖、雞場等鄉鎮，占該縣國土面積的22%。

③次適宜區。次適宜區占該縣國土面積的31%，主要分布在該縣中北部的鄉鎮，其中猴場所占比例最大。但在該區發展茶樹種植要謹慎。

④不適宜區。不適宜種植茶樹的區域分布最少，占地面積不到該縣國土面積的1%，主要分布在猴場的中部地區。

從區劃結果來看，普定縣種植茶樹適宜區和最適宜區已占該縣國土面積的68%，覆蓋該縣2/3以上的土地均可種植茶樹，不適宜種植茶樹的區域占地不到國土面積的1%，然而目前普定縣茶葉種植面積僅占耕地面積的4.92%，這說明該縣將來發展茶產業還有很大的空間。除氣候條件外，影響茶葉產量及品質高低的因素還很多，本文建立區劃因子空間分布模型時只考慮了經度、緯度和海拔高度，事實上氣象要素的空間分布還受坡度及坡向的影響，有待于今后進一步的研究和細化。

［1］王銀峰，張繼紅，翟海國.信陽茶區生態環境分析與評價［J］. 地理學報，1993，48(5):428－436.

［2］李湘閣，閔慶文，佘衛東.南京地區茶樹生長氣候適應性研究［J］.南京氣象學院學報，1995，18(4):572－577.

［3］張麗霞，趙淑娟，王桂雪，等.泰安市茶樹種植氣候條件分析［J］.中國農業氣象，2006，27(3):244－248.

［4］金志鳳，封秀燕.基于GIS的浙江省茶樹栽培氣候區劃［J］. 茶葉，2006，32(1):7 －10.

［5］羅京義，晏理華，徐大紅，等.銅仁地區茶樹生長的氣候適應性分析及優質綠茶種植區劃［J］.茶葉科學，2011，31(2):136－142.

［6］鄔倫，劉瑜，馬修軍，等.地理信息系統原理方法和應用［M］.北京:科學出版社，2011:324－326.

［7］馬旭清，國世友，鄒立堯.氣候資源小網格推算方法［J］.黑龍江氣象，2001，3:26 －28.

［8］韓湘玲，曲曼麗.作物生態學［M］.北京:氣象出版社，1991.

［9］胡家敏，谷曉平，徐永靈，等.貴州茶樹種植的氣候資源優勢及發展建議［J］. 貴州氣象，2008，32(3):19－21.

［10］金志鳳，黃敬峰，李波，等.基于GIS及氣候－土壤－地形因子的浙江省茶樹栽培適宜性評價［J］.農業工程學報，2011，27(3):231－235.

［11］朱壽燕，尹先龍，張文明，等.基于GIS的浙江仙居茶葉氣候生態區劃［J］. 地理空間信息，2012，10(2):25－27.

［12］萬素琴，向華，劉仁進，等.五峰縣采花毛尖系列有機茶種植氣候區劃［J］.湖北氣象，2004(2):13－15.

［13］劉永裕，劉梅，米浦強，等.紅花大果油茶種植氣候區劃［J］. 安徽農業科學，2009，37(24):11 818 －11 820.

［14］蘇永秀，李政，孫涵.基于GIS的廣西甘蔗種植氣候區劃［J］.中國農業氣象，2006，27(3):252－255.

［15］李楠，任穎，顧偉宗，等.基于GIS的山東省暴雨洪澇災害風險區劃［J］.中國農學通報，2010，26(20):313－317.