武夷山市茶園土壤有機質的空間分布特征

(福建省生態產業綠色技術重點實驗室，武夷學院生態與資源工程學院，福建 武夷山 354300) (武夷學院土木工程與建筑學院，福建 武夷山 354300) 福建省生態產業綠色技術重點實驗室，武夷學院生態與資源工程學院，福建武夷山354300

土壤是茶樹賴以生存的物質基礎，對茶樹高質、優產、高效至關重要[1]，其中土壤有機質(soil organic matter，SOM)是土壤肥力的重要組成部分，是反映茶園土壤熟化程度和肥力水平高低的重要指標[2]，SOM含量的高低直接影響茶樹生長發育以及茶葉的產量和品質。我國優質名茶產地SOM含量都在20g/kg以上[3]。土壤在形成和演化過程中受多種時空尺度環境因素(如氣候、土地利用、海拔、地貌地形等)的控制，因此SOM具有復雜的空間變異性和非平穩性[4,5]。隨著GIS技術的發展，將其與地統計學相結合的研究方法已經成為熱點[6]，趙汝東等[7]、陳琳等[8]、鄧小華等[9]、李枝樺等[10]對農田SOM的空間分異特征進行了研究，對農業生產起到了一定的指導作用。關于茶園SOM的空間分布已有大量研究[3，11～13]，但目前尚無有關武夷山市茶園SOM空間分布特征研究的報道。

武夷山市作為我國重要的茶產地，生產的武夷巖茶被贊譽為“茶中之王”，是中國十大名茶之一。武夷巖茶是地理標志保護產品，保護范圍限于《地理標志產品保護規定》批準的范圍，即福建省武夷山市所轄行政區域范圍[14]。武夷巖茶獨具神韻,素負盛名,蜚聲中外,主要是源于其優異的品質[15]。武夷山茶業在中國的茶葉市場中具有不可或缺的地位，同時也是閩北地區重要的經濟支柱產業。為此，以武夷巖茶不同種植區土壤為研究對象，探討表層土壤(0～20cm) SOM的空間分布特征，一方面可充實武夷巖茶土壤研究的相關內容，另一方面可為武夷巖茶的農業生產提供技術指導。

1 材料與方法

1.1 樣品采集與處理

圖1 采樣點的分布

通過實地踏查，選擇具有代表性且成片性較好的茶園作為研究區域。2018年10月14日至11月4日(冬茶采收基本結束，茶園未施肥之前)在福建省武夷山市境內星村鎮、五夫鎮、興田鎮、吳屯鄉、嵐谷鄉、洋莊鄉、武夷學院茶學科教基地、武夷山風景區(簡稱景區)等茶園種植區，根據茶園的種植面積大小布設采樣點數量，每個采樣點s形布設取樣點5個，混合成一個樣品。采樣深度0～20cm，共采集茶園土樣266份(圖1)。同時記錄樣點的植被情況、土壤類型、地形狀況、周圍環境等信息以及GPS地理坐標，茶園的經營管理措施通過咨詢茶農獲取。土壤樣品帶回實驗室在揀去石礫、植物根系和大于2mm的碎屑后自然風干，過0.149mm篩后儲藏備用， SOM含量采用重鉻酸鉀容量法測定。

1.2 數據處理

采用IBM SPSS Statistics 22.0對數據進行統計分析，運用GS+9.0進行半方差函數分析，利用ArcGIS 10.5繪制空間分布圖。

2 結果與分析

2.1 武夷山市茶園SOM含量的描述性統計分析

從SOM含量的描述性統計(見表1)可以看出，武夷山市不同種植區茶園SOM的平均含量大小表現為興田鎮(50.98g/kg)>武夷學院(43.45g/kg)>洋莊鄉(39.44g/kg)>星村鎮(37.55g/kg)>嵐谷鄉(35.08g/kg)>武夷山風景區(32.91g/kg)>吳屯鄉(32.57g/kg)>五夫鎮(28.64g/kg)。依據全國第2次土壤普查的土壤有機質分級標準[16]，武夷山市不同茶園SOM含量在五級及以上每個等級基本均有分布，在三級及以上分布的比例較大。興田鎮和武夷學院茶園SOM含量一、二級偏多，占到樣本總數85%以上，而五夫鎮的一、二級僅占到樣本總數36%，導致平均值較低，其他茶園的分布相對較為分散。從平均值來看，興田鎮和武夷學院茶園SOM含量處于一級(極高)水平，五夫鎮的處于三級(中等偏上)水平，其余茶園的處于二級(高)水平，而整個武夷山市的茶園土壤SOM含量(37.24g/kg)處于二級(高)水平。武夷山市的不同茶園土壤0～20cm表層SOM含量約92%達到國家綠色食品產地環境質量標準(NY/T 391—2013)中園地土壤的肥力Ⅰ級標準(>20g/kg)，約98%達高產茶園的標準(>15g/kg)。

變異系數(CV)反映樣本數據的離散程度，一般情況下，CV≤10%為較弱空間變異性，10%

表1 武夷山市茶園SOM含量描述性統計結果

2.2 武夷山市茶園SOM含量與其他地區的差異性分析

武夷山市不同種植區茶園SOM含量及變異系數(見表1)與江浙地區不同鄉鎮[3]的較為接近，土壤有機質都屬于中等變異。

武夷山市茶園SOM含量的平均值為37.24g/kg，浙江省松陽縣、江西省南昌縣、湖北省秭歸縣和云南省景洪市的平均值為17.17～31.51g/kg[18～21]，而陜西省商南縣的平均值僅為5.79g/kg[1]，相比而言，武夷山市茶園SOM含量相對較高。武夷山市茶園SOM含量的變異系數為37.69%，與浙江省松陽縣、江西省南昌縣、湖北省秭歸縣和云南省景洪市都屬于中等變異，而陜西省商南縣的變異系數為9.93%，屬于弱變異性。

2.3 半方差函數的擬合

半方差函數的計算一般要求數據符合正態或近似正態分布[22～25]。茶園SOM含量的數據K-S檢驗(P<0.05，2-tailed)結果顯示，將武夷山市8個種植區數據對數轉換后均符合正態分布。通過GS+9.0對武夷山市不同種植區茶園的SOM含量數據進行半方差函數分析，擬合參數詳見表2。武夷山市不同種植區茶園SOM含量的擬合模型不同：武夷山風景區、吳屯鎮和興田鎮為球形模型，嵐谷鄉為線性模型，五夫鎮、星村鎮和洋莊鄉為高斯模型，武夷學院為指數模型。

表2 半方差函數擬合參數

半方差函數的塊金值主要反映實驗誤差和小于實際取樣尺度引起的變異等隨機部分引起的空間變異程度，基臺值主要反映包括隨機變異和結構變異在內的總體的變異程度[26]，而基底效應(塊基比)，即塊金值與基臺值之比，表示由隨機部分引起的空間變異性占總體變異的比例[27]。通常塊基比小于0.25時，結構因子引起的空間變異占比較高，介于0.25～0.75時，表示空間變異受結構因子的影響為中等；大于0.75時，則是隨機因子引起的空間變異所占比例較高[28,29]。

由表4可知，洋莊鄉的塊基比為0.0345，遠小于0.25，具有強烈的空間相關性，在一定程度上反映洋莊鄉的茶園SOM含量受土壤類型、成土母質、地形條件、氣候等環境作用的結構性因素引起的空間自相關變異性占到總變異的 95%以上；嵐谷鄉的塊基比為1，大于0.75，說明嵐谷鄉的茶園SOM空間相關性較弱，且空間變異的來源受采樣、檢測、施肥、耕作措施、種植制度等隨機因素引起的空間變異較大。武夷山市其他6個種植區的塊基比為0.25～0.75，具有中等的空間相關性。

變程是半方差達到基臺值時相應的距離，代表SOM的自相關范圍，SOM在變程范圍內具有空間自相關性,在變程范圍外則沒有。武夷山市8個研究區的茶園SOM變程由大到小依次為星村鎮>武夷山風景區>武夷學院>五夫鎮>吳屯鄉>興田鎮>嵐谷鄉>洋莊鄉(見表2)，說明影響不同種植區茶園SOM含量的過程作用在不同的尺度上[30]。

2.4 武夷山市茶園SOM空間分布特征

圖2 武夷山市茶園SOM含量的空間分布

武夷山市茶園SOM的空間分布[31～33]如圖2所示，SOM含量極高(大于40g/kg)的地方主要集中在南部的興田鎮和中部的武夷學院，西北部SOM含量中等，東部SOM含量相對較低。經調查走訪發現，地處東部的嵐谷鄉、吳屯鄉和五夫鎮的茶產業在當地的農業生產中所占的比重較小，在茶產業中的投入相對較少，致使東部鄉鎮茶園SOM含量較低。

南部的興田鎮近幾年對茶產業的投入逐年增加，可能是茶園SOM含量高值區多分布在南部興田鎮的原因之一。武夷學院的茶學科教基地面積相對較小且茶園分布比較集中，園內布設水土流失綜合治理的坡改梯工程，部分茶園套種烏桕、丹桂、花石榴等，套種模式下植物、落葉等殘體積累更多，有利于有機質的積累；同時在梯壁和梯埂種植麥冬、白花葛、紫葉酢醬草等可有效防止水土流失，從而防止有機質缺失，提升了茶園的環境質量，提高了SOM含量。西部的星村鎮是中國烏龍茶和紅茶的發源地，是“福建產茶明星鄉鎮”，有“武夷巖茶第一鎮”之稱，其茶園多屬于大型茶場和大承包戶的茶園，由于大規模種植的大型茶場種植施肥和經營管理比較科學有效且投入較多，茶園SOM含量相對較高。武夷山風景區是武夷巖茶正巖產區，武夷山正巖茶園面積有限，資源珍貴，正巖茶園的產區較為零散，產區內道路崎嶇，交通不便，有的茶園因此減少了有機肥施用數量，結果導致景區內茶園土壤SOM含量相對較低。

圖3 武夷山市茶園SOM含量與海拔的關系

全市茶園SOM含量大于40g/kg的采樣點占34%，30～40g/kg的占33%，20～30g/kg的占25%，小于20g/kg的僅占8%，因此武夷山市茶園SOM含量處于整體偏高的狀況。

2.5 武夷山市茶園SOM含量與海拔高度的相關性分析

地形差異是造成SOM含量差異的主要因素之一，對采樣區域的平均海拔與SOM的平均含量進行相關性分析，如圖3所示。由圖3可知，茶園SOM含量與取樣點的海拔高度的呈一定的負相關，但差異不顯著，這與楊如興等[34]、王軍等[35]、張建杰等[36]的研究結果一致，但與王宗明等[37]、李超等[38]、張偉等[39]的結論不一致，這可能與研究區域茶園的地理氣候、生態環境、土壤特點及培肥措施等不同有關。長期或大量施用化肥會導致土壤和茶葉的鮮葉質量下降并可能進一步造成環境污染。將化肥與有機肥配合施用，在有效增加土壤肥力的同時可以保持土壤質量。武夷山部分茶企從提升茶葉質量和保護環境可持續發展的角度出發，對茶園土壤只施用有機肥。但在高海拔區域，由于交通原因導致施用有機肥的成本可能會大大增加，因此有的茶園減少了有機肥的施用量，這可能是分布在較高海拔區域的茶園SOM含量有下降趨勢的原因之一。

3 結論

1)武夷山市8個種植區茶園SOM的平均含量大小表現為興田鎮(50.98g/kg)>武夷學院(43.45g/kg)>洋莊鄉(39.44g/kg)>星村鎮(37.55g/kg)>嵐谷鄉(35.08g/kg)>武夷山風景區(32.91g/kg)>吳屯鄉(32.57g/kg)>五夫鎮(28.64g/kg)。武夷山市茶園SOM的平均含量(37.24g/kg)處于二級(高)水平(全國第二次土壤普查的土壤有機質分級標準)。其中約92%的樣品達到國家綠色食品產地環境質量標準(NY/T 391—2013)中園地土壤的肥力Ⅰ級標準(>20g/kg)，約98%達高產茶園的標準(>15g/kg)。

2)武夷山市茶園SOM含量的變異系數介于22.75%～48.70%之間，屬于中等空間變異。不同茶園SOM數據對數轉換后均符合正態分布，武夷山風景區、吳屯鎮和興田鎮茶園SOM空間分布為球形模型，嵐谷鄉為線性模型，五夫鎮、星村鎮和洋莊鄉為高斯模型，武夷學院為指數模型。武夷山風景區、武夷學院和星村鎮茶園SOM空間自相關距離較大，分別為 2.172、2.11、 2.58km，其他種植區茶園SOM空間自相關距離相對較小。洋莊鄉的塊基比為0.0345，具有強烈的空間相關性，表明其空間變異受結構性因素影響較大；嵐谷鄉的塊基比為1，空間相關性較弱，說明受隨機因素引起的空間變異較大，其他6個種植區的塊基比為0.25～0.75，具有中等的空間相關性。

3)武夷山市茶園SOM含量極高(大于40g/kg)的地方主要集中在南部的興田鎮和中部的武夷學院，其次是西部，東部相對較低。茶園SOM含量與取樣點的海拔高度呈一定的負相關，但差異不顯著。