吳炳孫 馮家東 王晶晶 王桂花

（1中國熱帶農(nóng)業(yè)科學院橡膠研究所農(nóng)業(yè)部橡膠樹生物學與遺傳資源利用重點實驗室/省部共建國家重點實驗室培育基地-海南省熱帶作物栽培生理學重點實驗室/農(nóng)業(yè)部儋州熱帶作物科學觀測實驗站 海南儋州 571737；2海南大學林學院 海南?？?570228；3中國熱帶農(nóng)業(yè)科學院橡膠研究所 海南?？?571101）

海南省天然橡膠的種植面積超過53萬hm2［1］，其林下有巨大的發(fā)展空間。在膠價低迷常態(tài)化下，發(fā)展林下經(jīng)濟是一種有效提高膠園綜合效益的途徑，其中膠園林下間作能同時充分利用膠園的水肥光熱資源，還能提高膠園經(jīng)濟效益，已成為提高橡膠生產(chǎn)效益的一項重要措施［2］。

土壤水分是土壤各種生物化學反應過程的重要影響因子［3］，也是橡膠林吸收養(yǎng)分的必需載體，對橡膠林生長有重要影響。土壤水分受到氣候、植被、地形等自然因素及施肥、耕作、灌溉等人為因素共同影響［4］，在不同空間位置上具有明顯的差異，具有空間變異性。橡膠間作復合種植對土壤理化性質(zhì)產(chǎn)生一定的影響［5-8］，尤其對土壤水分利用方面，研究表明膠園間作能顯著提高土壤含水量，還能提高橡膠樹水分利用效率［9］，但膠園間作對土壤水分空間變異的影響未見報道，為此以膠園間作益智為研究對象，橡膠單作為對照，采用地統(tǒng)計方法分析空間異同性，探明空間分布特性及合理采樣數(shù)量，為膠園間作土壤水分評價奠定基礎，而且也為科學開展橡膠林的土壤水分和培育管理提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究樣區(qū)位于海南省西部中國熱帶農(nóng)業(yè)科學院試驗農(nóng)場三隊（19°32′47″N，109°28′30″E），海拔高度為114 m，地形平坦，土壤為花崗巖磚紅壤，土層厚度約100 cm，多為沙質(zhì)黏壤土。屬典型熱帶海島季風氣候，年平均氣溫23.5～24.1℃，年日照時間2 450.0 h，無霜期365d，受熱帶季風影響旱雨季分化明顯，一年當中5～10月份是雨季，11月至翌年4月為旱季，年均降水量約1 500～1 900 mm，雨季降雨量占全年降雨量80%以上，干季降雨量少。選取益智-橡膠復合種植（以下簡稱益-膠）與橡膠單作為研究樣地，研究樣地的橡膠林品系均為熱研7-33-97，2005年種植，2012年開割，橡膠樹種植規(guī)格均3.0 m×7.0 m，地勢平坦、海拔高度一致，樣地相距500 m內(nèi)，其中益智-橡膠復合種植樣地的益智（Alpinia oxyphylla Miq.）間作距離橡膠樹約2 m左右，種植規(guī)格為1.5 m×1.5 m，冠高1.2 m左右；橡膠單作樣地林下為草本層，草本層冠高約為0.2 m左右。

1.2 土壤水分采樣方法

在益智-橡膠復合種植樣地和橡膠單作樣地分別設置30 m×70 m面積為取樣區(qū)，即橡膠樹11行×11株的范圍區(qū)域。取樣位置位于橡膠樹行間，在每行間的南北方向以離橡膠樹0.5 m為起始線，沿著南北方向間隔2 m，東西方向間隔3 m劃分網(wǎng)格（圖1），每行間共劃分30個網(wǎng)格，每個樣地共布設300個取樣網(wǎng)格。于連續(xù)幾天晴天、無降雨的2020年8月30日在每個網(wǎng)格的中心用TDR300土壤水分測定儀隨機測3個點0～10 cm深度土壤含水量，將3個點的平均數(shù)據(jù)代表該網(wǎng)格樣點的樣本值。

1.3 空間變異分析

半方差函數(shù)是地統(tǒng)計學應用最廣泛的空間格局描述的基本工具［10-12］，其計算公式如下：

式（1）中，γ(h)為半方差函數(shù)值；N(h)是間距為向量h的點對總數(shù)；Z(xi)為區(qū)域化變量Z在xi處的實測值；Z(xi+h)是與xi距離為向量h處樣點的值。

圖1 三行橡膠樹間的土壤含水量測定樣點示意圖

基臺值（sill，C0+C）、塊金值（nugget，C0）和變程（range）是半方差函數(shù)中重要的參數(shù)，其中基臺值表示變量在空間上的總變異性大小，塊金值則反映隨機因素引起的空間變異，變程表示空間變異的范圍。而塊金值與基臺值的比值［C0/（C0+C）］的大小，又稱塊金效應，可確定結(jié)構(gòu)性因素（如成土母質(zhì)、氣候、地形等自然因素）與隨機性因素（如耕作措施、施肥及種植制度等人為因素）在空間變異的主導性。一般認為，該比值＜0.25時，具有強烈的空間相關(guān)性，即要素空間變異主要受結(jié)構(gòu)性因素影響；比值為0.25～0.75時，具有中等程度的空間自相關(guān)，即要素的空間變異受到結(jié)構(gòu)性因素與結(jié)構(gòu)性因素共同影響；比值＞0.75時，空間自相關(guān)性很弱，即要素的空間變異主要受隨機性因素影響。

1.4 數(shù)據(jù)處理與制圖

土壤水分含量的描述性統(tǒng)采用SPSS（17.0）計算，半方差函數(shù)擬合與分析在GS+9.0完成，空間制圖由ArcGIS10.6實現(xiàn)。

2 結(jié)果與分析

2.1 橡膠園單/間作模式土壤水分統(tǒng)計特征

各膠園樣地土壤含水量的描述性統(tǒng)計結(jié)果如表1所示。橡膠間作模式的土壤含水量顯著高于橡膠單作（p＜0.05），其中，益-膠復合樣地土壤含水量均值為21.92%，而橡膠單作樣地為19.96%。益-膠復合樣地和橡膠單作樣地土壤含水量的變異系數(shù)分別為12.28%、9.55%，均小于25%，屬于弱變異性。經(jīng)樣本Kolmogorov-Smirnov（K-S）正態(tài)檢驗和對數(shù)正態(tài)檢驗結(jié)果表明，各樣地土壤含水量呈對數(shù)正態(tài)分布。

表1 橡膠單作與間作模式的土壤水分的描述性統(tǒng)計

2.2 橡膠園單/間作模式土壤水分含量的空間異質(zhì)性分析

運用GS+9.0對各樣地土壤水分含量的半方差函數(shù)進行反復模擬，獲取最優(yōu)半方差擬合函數(shù)參數(shù)（見表2）和擬合曲線（見圖2），其中橡膠單作樣地的土壤水分空間變異用指數(shù)模型擬合較優(yōu)，而益-膠復合樣地則為球形模型擬合較優(yōu)。益-膠復合樣地的塊金值大于橡膠單作樣地，較大的塊金值可能與間作益智人為干擾有關(guān)，反映出在小于2 m研究尺度上存在作物的株間差異和土壤擾動等影響。另外，塊金值與基臺值之比（塊金效應）表示隨機部分引起的空間異質(zhì)性占系統(tǒng)總變異的比例。益-膠復合樣地的塊金值與基臺值之比為28.39%，表明土壤水分具有中等程度的空間相關(guān)性，說明結(jié)構(gòu)性（如地形、降水、氣溫等自然因素）和隨機性因素（如施肥、耕作模式等人為因素）共同引起的空間變異，而橡膠單作樣地的塊金值與基臺值之比為12.63%，說明土壤水分具有強烈的空間相關(guān)性，這主要是由自然因素引起空間變異，并且顯著小于益-膠復合樣地，表明橡膠園間作作物會削弱自然因素引起的土壤水分空間變異性。

表2 橡膠單作與間作模式的土壤水分半方差擬合函數(shù)參數(shù)

圖2 不同膠園土壤水分半變異函數(shù)圖

2.3 橡膠園單/間作模式土壤水分空間分布特性

根據(jù)各膠園半方差模型擬合的參數(shù)繪制出橡膠單作與間作模式的土壤水分含量空間分布圖（圖3～4），從圖3、圖4中可以看出各膠園土壤含水量表現(xiàn)出一定的空間異質(zhì)性，益-膠復合膠園土壤水分空間異質(zhì)性較橡膠單作膠園土壤水分強度高，且益-膠復合膠園土壤含水量也比橡膠單作膠園高。橡膠單作膠園土壤水分空間分布呈大圖斑分布，碎片化圖斑較少，空間聚群少，行間與非行間土壤含水量差異不顯著（F=0.011，p＞0.05），整體空間變異小，這可能與橡膠單作地表覆蓋為草本層，草本層密集且單層覆蓋，影響單一使土壤含水量空間變化不大。而益-膠復合膠園土壤水分空間分布呈多環(huán)圖斑分布，一般由內(nèi)環(huán)土壤含水量低逐漸向外環(huán)土壤含水量高擴散，且存在以內(nèi)環(huán)為中心多環(huán)擴散形式的多個空間聚群區(qū)域，呈明顯的階梯狀且連續(xù)性較好，這可能與間作益智使橡膠園呈多層立體冠層相關(guān)，相比橡膠單作橡膠園，增加了益智冠層，影響橡膠園內(nèi)降水分配與地表蒸發(fā)，造成空間上存相關(guān)性，導致斑塊變化相似，連續(xù)性好，這也是使變程較大的原因。

2.4 橡膠園單/間作模式土壤水分合理采樣數(shù)量

土壤存在空間變異性，通過增加土壤采樣數(shù)量可以提高估算土壤含水量精度，但需消耗大量的勞力與財力，而樣本數(shù)太小，則又會造成估算不精確。因此科學地確定采樣數(shù)量對實際應用中具有重要的意義。根據(jù)區(qū)域隨機抽樣理論，一般采用Cochran提出針對隨機采樣而構(gòu)建的采樣數(shù)量計算公式［13］。

式（2）中：n為合理采樣數(shù)量；CV為土壤含水量的變異系數(shù)；tα/2表示顯著性水平為α、自由度為n-1條件下的t分布分位數(shù)；E表示估計均值的允許誤差。

根據(jù)計算公式，設定95%和90%2種置信水平，以及相對各膠園土壤含水量均值的1%、5%和10%等3種允許誤差等級，分別計算了益-膠復合膠園和橡膠單作膠園所需的合理采樣點數(shù)量（表3）。由表3可知，在相同置信水平與允許誤差下，益-膠復合膠園的所需采樣數(shù)量均大于橡膠單作膠園，且隨允許誤差等級的增加其差值而迅速減少，如在95%置信水平時，益-膠復合膠園所需采樣數(shù)量在1%、5%和10%允許誤差等級分別比橡膠單作膠園大231、9和2。在兩置信水平上，各膠園所需采樣點數(shù)量均隨允許誤差的增加而迅速減少，其中誤差由1%增加為5%時，采樣點減少幅度最大，如在95%置信水平，益-膠復合膠園和橡膠單作膠園所需的合理采樣點數(shù)量分別減少563、341個。

圖3 橡膠單作膠園土壤水分含水量空間分布圖

圖4 益智-橡膠復合種植膠園土壤水分含量空間分布圖

3 討論與結(jié)論

本研究通過對益智-橡膠復合種植膠園與橡膠單作膠園土壤水分進行統(tǒng)計學分析及空間異質(zhì)性分析發(fā)現(xiàn)，橡膠間作對膠園土壤水分含量及空間分布有重要影響。相比橡膠單作，橡膠間作益智能顯著提高土壤含水量，0～10 cm土層土壤提高9.82%，這與前人研究結(jié)果一致［6］。橡膠間作膠園土壤含水量顯著高橡膠單作膠園，這可能因為橡膠間作復合種植膠園存在多層植被，本研究中益智冠層高1.2 m，且?guī)缀跬耆采w膠園行間，多層植被覆蓋可以很大程度上減弱表層土壤水的蒸發(fā)［8］，從而增強土壤的保水能力。橡膠單作膠園土壤水分的塊金值/基臺值比為12.63%，具有強烈的空間相關(guān)性，主要由自然因素引起空間變異；而益智-橡膠復合種植膠園土壤水分的塊金值/基臺值比為28.39%，具有中等程度的空間相關(guān)性，受到自然因素與人為因素共同影響，可見橡膠間作益智的塊金值/基臺值比顯著大于橡膠單作，說明橡膠園間作作物會削弱自然因素引起的土壤水分空間變異性。這可能因為土壤特性具有時間和空間的變化特征，其變化特征取決于自然因素（如土壤母質(zhì)、地形、氣候等）和人為因素（如土壤施肥、耕作、灌溉等）的綜合作用［14-15］，本研究中間作益智時常要進行施肥、摘果等農(nóng)事活動，人為頻繁干擾土壤，同時益智根系主要分布在0～20 cm土層，根系生長活動也深刻影響土壤，提高了人為因素引起的土壤水分空間變異性，導致膠園土壤水分空間變異增大。益智-橡膠復合種植膠園的土壤水分空間分布上呈明顯的階梯狀且連續(xù)性較好，且存在以內(nèi)環(huán)為中心多環(huán)擴散形式的多個空間聚群區(qū)域，空間變異大，而橡膠單作膠園呈大圖斑分布，碎片化圖斑較少，空間聚群少，空間變異小。說明間作益智顯著影響土壤水分空間分布，使土壤水分在空間上以某優(yōu)勢植被或某區(qū)域為中心向外擴散變化，這可能與間作益智增加了益智冠層，使橡膠園呈多層立體冠層，影響橡膠園內(nèi)降水分配與地表蒸發(fā)，同時益智生長發(fā)生蒸騰吸水，造成土壤水分空間上存相關(guān)性，導致斑塊變化相似，連續(xù)性好。

表3 估算各膠園土壤水分均值所需要的樣本數(shù)量

合理采樣點數(shù)量的確定是優(yōu)化土壤采樣調(diào)查方案的基礎。本研究在益智-橡膠復合種植膠園和橡膠單作膠園等樣地小尺度上，通過設定95%和90%2個置信水平上估算了1%、5%和10%3種估算誤差等級所需的土壤水分合理采樣點數(shù)量，發(fā)現(xiàn)益智-橡膠復合種植膠園的所有合理采樣點數(shù)量均比橡膠單作膠園大，這是因為益智-橡膠復合種植膠園的土壤水分變異系數(shù)大，要達到一定的精度需通過增加采樣數(shù)量減少誤差。一般取95%置信水平且采樣精度為均值5%條件下為合理采樣數(shù)量，益智-橡膠復合種植膠園和橡膠單作膠園土壤水分的合理采樣數(shù)量分別為23個和14個，但是制定土壤采樣方案時，要結(jié)合土壤水分空間變異性大小，根據(jù)所需的置信水平和允許誤差來共同確定區(qū)域合理采樣點數(shù)量。