喀斯特常綠落葉闊葉混交林土壤磷鉀養分空間異質性

王 華,陳 莉,3,*,宋 敏,宋同清,曾馥平,彭晚霞,杜 虎,蘇 樑,3

1 湖南農業大學生物科學技術學院，長沙 410128 2 中國科學院亞熱帶農業生態研究所亞熱帶農業生態過程重點實驗室，長沙 410125 3 中國科學院環江喀斯特生態系統觀測研究站，環江 547100 4 湖南農業大學農學院,長沙 410128

土壤是一種形態和演化過程都十分復雜的自然綜合體,受氣候、生物、母質、地形、成土時間等成土因素的影響,具有復雜性和時空變異性[1]。即使在土壤類型和質地相同的區域內,同一時刻土壤特性在空間上也有明顯的差異[2]。土壤養分是土壤提供的植物生長發育所必需的營養元素,是土壤肥力的物質基礎,也是土地評價和管理的重要指標之一。由于受到自然因素和人為因素的共同作用,土壤養分也具有高度的空間變異性[3- 6]。土壤P、K不僅是植物生長發育所必需的營養元素,也是影響區域水體生態環境的重要屬性,因此,研究其時空變異對喀斯特生態系統的恢復和重建等具有重要的理論與現實指導意義,以實現喀斯特地區的生態、經濟、社會的和諧和可持續發展[7]。

一些學者[8- 9]利用土壤學家Fisher所創立的傳統統計方法對土壤性狀的空間分布進行了定量描述,但由于此方法要求樣本的取值相互獨立,給土壤空間變異的研究帶來了很大的局限性。20世紀60年代Matheron創立了地統計學理論[10],為土壤空間變異性研究提供了有效的方法,掀起了國外學者對土壤空間變異性研究的熱潮[11- 13]。目前,土壤空間變異性研究已成為土壤科學的重要研究內容和熱點[14- 16]。如：Corss和Schlesinger[17]認為小尺度上的土壤異質性分布與群落優勢種分布一致；Farlely等[18]在研究闊葉林地時發現,在20 cm小尺度上,土壤的養分含量差異非常大,同時也有研究表明[19- 21]土層深度對森林土壤養分含量的影響差異顯著,且隨著土層的加深土壤有機質和速效鉀的含量降低；宋同清等[22]認為中小尺度上微生境土壤理化特性等因素可能是推動植物群落結構變化的主力軍。有關土壤養分空間分布的研究頗多,但國內外的研究多數集中在小尺度、單因素[23- 25],而在較大尺度上的喀斯特土壤養分空間分布研究較少[26- 38]。喀斯特地區人口密度高、增長速度快、可利用土地少、生產方式落后等不利因素使得該區域生態系統承受壓力大,水土流失嚴重,基巖裸露率高,最終導致愈來愈嚴重的石漠化現象。本文以廣西木論國家級自然保護區典型森林(常綠落葉闊葉混交林)土壤為研究對象,應用經典統計和地統計法學方法分析了該森林土壤磷鉀養分TP、TK、AP、AK的空間格局、相關性與變異性等特征；同時,通過Pearson相關性分析該區域內坡度、坡向、土層深度等環境因子對土壤養分分布的影響,旨在為喀斯特地區的土壤養分庫的評價、生態保護和恢復研究提供理論支持。

1 研究區概況與研究方法

1.1 研究區概況

研究區位于廣西環江毛南族自治縣西北部木論國家級自然保護區25°07′01″—25°12′22″N,107°54′01″—108°05′51″E,南緣云貴高原,北連貴州茂蘭國家級自然保護區,總面積8969 hm2,海拔400—1000 m。該區域為中亞熱帶季風氣候,年平均氣溫15—18.7℃,最冷月(1月)平均氣溫3.4—8.7℃,最熱月(7月)平均氣溫23—26.7℃,≥10℃,年積溫4700—6300℃；年降雨量1530—1820 mm,集中于4—8月,年平均相對濕度為80%—90%,無霜期235—290 d,氣候溫暖,雨量充沛。土壤主要為石灰土和零星分布的硅質土,河流及水時空分布不均,地域性來水差異大。研究區植被為中亞熱帶石灰巖常綠落葉闊葉混交林,屬隱域性喀斯特森林植被頂極群落類型,是目前世界上同緯度地區殘存下來的僅有的、原生性強、相對穩定的喀斯特森林生態系統,也是喀斯特區原生性森林分布面積最大的區域。植被成層現象比較明顯。

1.2 樣品采集與分析

圖1 木論樣地地形圖Fig.1 The contour map of the 25 hm2 Mulun forest plot

通過勘測,于2014年3月在研究區內設置一塊500 m×500 m的大樣地,作為長期監測樣地(圖1)。參照CTFS標準采用規則網格法(20 m×20 m)進行樣點布置,用全站儀和GPS基站相結合的測量方法設置和劃分樣地,參考木論自然保護區的森林狀況和植被、地形分布狀況,共設置625個小樣方,樣方4個角用水泥樁作永久標記。在每個20 m×20 m樣方內,按照其中心點半徑2 m范圍內五點取樣法,用直徑5 cm的土鉆取表層(0—10 cm)土壤5個混合成一個樣品,代表該樣點樣本。同時記錄各樣方坡向、坡位、海拔、土層厚度等相關立地因子。將所有混合土樣帶回實驗室風干、過篩,以備分析測定全磷(TP)、全鉀(TK)、速效磷(AP)、速效鉀(AK)。坡向視為坡度方向或山體所面對的羅盤方向,以度為單位按逆時針方向進行測量。巖石裸露率以樣點周圍2米范圍內巖石出露的面積估算。凹凸度為樣方的海拔減去該樣方相鄰的8個樣方海拔的平均值。全磷(TP)的測定采用氫氧化鈉熔融—鉬銻抗顯色—紫外分光光度法；全鉀(TK)的測定采用氫氧化鈉熔融—火焰原子吸收分光光度法；速效磷(AP)的測定采用0.5 mol/L碳酸氫鈉浸提—鉬銻抗顯色—紫外分光光度法；速效鉀(AK)的測定采用乙酸銨浸提—火焰原子吸收分光光度法[29]。

1.3 數據處理

使用SPSS 18.0(統計分析與檢驗、Pearson相關性)、Minitab 16(Box-Cox法)、GS+(地統計分析)、ArcGIS 9.3(克里格插值圖)等軟件,采用經典統計學和地統計學方法對樣本數據進行分析處理。

1.3.1 半方差函數分析

半方差函數是地統計學研究中特有的基本工具,用來描述屬性變量在空間上的差異程度,它是距離的函數,也是應用最廣泛的空間格局描述工具,其公式如下:

(1)

式中,γ(h)為半方差函數值；h為兩樣本點的空間距離；N(h)為間距等于h的點對數；Z(xi)為區域化變量Z在位置xi的實測值；Z(xi+h)為與xi距離為h處樣點的值。一般認為半方差函數只有在最大間隔的l/2內才有意義。塊金值 (C0)、基臺值(C0+C)和變程(A)為半方差分析的重要參數,塊金值與基臺值的比值用來表示空間變異程度,一般認為塊基比C0/(C0+C)<25%時,空間變量為強烈的空間自相關,在25%—75%之間時,為中等空間自相關,>75%為弱空間自相關。本研究中所有半方差函數的有效滯后距均設為最大滯后距的1/2,采用指數模型進行擬合。

1.3.2 克里格插值

克里格插值是建立在半方差函數理論的基礎上,在一定區域內對變量的取值進行無偏最優估計的一種方法,也是最簡單、應用最廣泛的一種插值方法。在普通克里格插值中,對于任意待估點(塊段)x0的實際值Z(x0),其估計值 Z#(x0)通過該待估點(塊段)周圍的n個有效樣點Z(xi)的線性組合得到。

(2)

式中,λi為權重系數。在滿足無偏條件下,要使預測值的估計方差最小,即:

σ2=Var[Z(x0)-Z#(x0)]=min

(3)

2 結果與分析

2.1 土壤磷鉀養分經典統計描述

采用樣本均值±3倍標準差來識別特異值,在此區間外的數據均定為特異值,分別用最大和最小值代替,后續計算均采用處理后的原始數據。由表1可知,土壤中TP、TK、AP、AK的含量分別為(1.60±0.76) g/kg、(5.42±2.74) g/kg、(5.74±3.63) mg/kg、(5.20±2.96) mg/kg。從變異系數看,土壤磷鉀養分的變異系數范圍為47.66%—63.22%,均屬中等變異水平。其中變異系數最大的指標為AP,為63.22%。通過對偏度、峰度的觀測以及采用 K-S 法對4種土壤磷鉀養分數據進行非參數檢驗,在5%的檢驗水平下,結果均不符合正態分布。用Box-Cox法對不符合的樣本進行轉換,轉換后數據再進行單樣本Kolmogorov-Smirnov檢驗,結果表明均符合正態分布。

表1 土壤磷鉀養分含量的描述性統計分析及正態分布檢驗

*表示經Box-cox轉換后為正態分布

2.2 土壤磷鉀養分空間自相關分析

Moran′sI的數值大小反應了土壤養分指標的空間自相關性大小。由圖2可以看出,研究區土壤磷鉀養分呈現出一定的結構性,TP與AP之間具有相似的空間結構。AP的Moran′sI系數在-0.124—0.497之間,I系數達0.497,表明相對于其他養分指標,AP的對空間的依賴程度較大；TP、TK、AK的Moran′sI在-0.043—0.313之間；各養分空間自相關的排序大小為：AP>TP>AK>TK。所有養分指標的Moran′sI系數均隨著滯后距的增大而減小,在220 m左右出現拐點,Moran′sI下降為負相關并趨于平緩,說明在該研究區220 m范圍內森林原生性較好。

圖2 研究區土壤磷鉀養分空間相關圖Fig.2 Spatial correlation diagram of soil phosphorus and potassium nutrients in the study area

2.3 土壤磷鉀養分半方差函數

土壤特性空間異質性的經典統計分析中將土壤特性看作隨機的、彼此相互獨立的。然而實際中,在一定的空間范圍內,土壤性狀存在一定的空間依賴性。由表2和圖3可看出,研究區土壤磷鉀養分TP、TK、AP、AK變異函數值的最佳擬合模型均為指數模型,決定系數R2均很高,在0.671—0.995之間,能很好地反映土壤磷鉀養分的空間結構特征。研究區各土壤磷鉀養分指標的塊金值/基臺值C0/(C0+C)依次為TK>AK>AP>TP,其中TP、AP塊金值/基臺值C0/(C0+C) 為50.4%、64.4%,空間自相關中等,主要受結構因子影響；TK、AK塊金值/基臺值分別為88.4%、86.6%,為弱空間自相關。變程反映土壤性狀的有效空間自相關距離,研究區TP、AP的變程分別為336.00、373.50 m,表明空間連續性較好；TK、AK的變程在33.30—64.50 m之間,空間依賴性較強。

表2 研究區土壤磷鉀養分半方差函數模型類型及參數

圖3 研究區土壤磷鉀養分半方差函數Fig.3 Semivariogram of soil phosphorus and potassium nutrientsin the study area

2.4 土壤磷鉀養分空間分布格局

為了能更直觀地反映研究區土壤參數的空間分布特征,根據半方差函數模型Arc GIS 9.2 軟件進行Kriging 插值分析,分別繪制全磷(TP)、全鉀(TK)、速效磷(AP)、速效鉀(AK)4種土壤養分的空間分布圖(圖4)。從圖中可看出,研究區土壤養分TP表現為坡下(含洼地)含量高,坡上含量較低,在偏西北方向出現最大值；AK表現為中坡含量高于洼地含量；AP呈斑塊狀分布,規律性不強,樣地正北方向出現最大值；TK斑塊破碎化程度較高。

圖4 土壤磷鉀養分空間分布圖-克里格插值圖Fig.4 Spatial distribution diagram of soil phosphorus and potassium nutrients—Kriging maps

2.5 土壤磷鉀養分空間變異的影響因素

經過Pearson相關性分析計算4個土壤磷鉀養分指標的相關性矩陣(表3),可以看出研究區4種養分指標TP、TK、AP、AK均與海拔呈負相關,其中TP、TK、AP與海拔呈極顯著負相關；各指標與坡向和地面凹凸度均為負相關,其中坡向與TP呈極顯著負相關關系,地面凹凸度與TP、TK呈極顯著負相關；土層深度與AK呈極顯著負相關；巖石裸露率與AP呈顯著負相關,與TP、TK、AK呈正相關關系；坡度與TP、AP呈極顯著負相關。說明土壤磷鉀養分主要受海拔、坡度、地面凹凸度的影響。海拔較低、坡度較緩的點,土壤磷鉀養分的含量也較高。

表3 土壤磷鉀養分與環境因子Pearson相關性分析

**.在0.01水平(雙側)上顯著相關,*.在0.05水平(雙側)上顯著相關

3 討論

研究區位于木論自然保護區緩沖帶內,為一典型的峰叢洼地單元,森林保存完好、人為干擾少、植被多樣性高、地形地貌復雜多變、微生境豐富[7]。經典描述性統計分析發現,研究區土壤各磷鉀養分指標差異較大,這可能與喀斯特地區土壤的高異質性及營養元素的淋溶特征密切相關。喀斯特常綠落葉闊葉林土壤磷鉀養分變異系數在47.66%—63.22%之間,均為中等變異；TP、TK、AP、AK養分含量的變化范圍與劉璐等[30]的研究結果相符,即原生性較好的喀斯特森林的養分水平皆呈現出較高的水平,同時也與自然保護區內的生態結構功能和保護區的頂極群落優勢密切相關。研究區常綠落葉闊葉混交林土壤因子的基本特征與其他亞熱帶森林有所不同。與浙江常綠闊葉林土壤相比[31],本研究區土壤TP含量均值高于浙江常綠闊葉林水平；與鄂西南亞熱帶常綠落葉闊葉混交林相比[32],研究區土壤TP、AP較高,AK較低。

基于前期對喀斯特土壤主要養分和水分的空間異質性研究[33- 36],進一步對喀斯特常綠落葉闊葉混交林土壤磷鉀養分含量進行分析發現,研究區土壤磷鉀養分呈現出一定的結構性,所有養分指標的Moran′sI系數均隨著滯后距的增大而減小。研究區TP、TK、AP、AK變異函數值的最佳擬合模型均為指數模型,對土壤空間變異具有極高的擬合度,表明研究區土壤具有良好的空間結構。土壤性狀的空間異質性是結構因子和隨機因子共同作用的結果。影響土壤空間變異的結構因素為母質、氣候、地形等自然地理要素,隨機因素則為取樣設計、測定誤差與干擾等人為原因。喀斯特常綠落葉闊葉林土壤TP、AP呈空間自相關中等,說明其變異是由結構性變異(環境因子等)和隨機因素共同引起的；TK、AK為弱空間自相關,說明隨機因素對這兩項土壤養分指標空間分布的的影響較大, 可能與干擾、試驗誤差有關。研究區TP、AP的變程較高,空間連續性較好；TK、AK空間依賴性較強。與湖北[32]及浙江常綠落葉闊葉混交林[37]的同類研究相比具有相同點與不同點：相同點在于均發現土壤AP的變異系數相對其他因子較高,這表明相比其他養分指標,土壤AP易受其他因素影響；不同點在于本文發現土壤AP的半方差變程最大,而湖北與浙江同類森林中土壤AP較其他因子的半方差擬合的變程相對較小,這可能與喀斯特生態系統具有獨特的二元水文結構及復雜的微生境有關,生態系統的結構與功能對土壤空間分布格局產生重要影響。

空間異質性導致空間格局的存在,Kriging 等值線圖能更深刻、全面和直觀地反映土壤磷鉀的空間分布。研究區土壤TP表現為坡下(含洼地)含量高,坡上含量較低,這是由于土壤中的P主要來源于基巖的風化,并通過植物的表聚作用在表層土壤積累,因而凋落物P是表層土壤P的直接來源[38-39]；潘復靜等[40]研究表明,在峰叢洼地,凋落物P沿坡面具有明顯的空間分異規律,隨坡位的升高而降低,這也是本研究土壤TP坡上含量較低的主要原因。AP、TK呈斑塊破碎化,這可能與土壤的淋溶特征有關。由于土壤養分受多種因子的綜合作用,不同養分的破碎化程度不一樣,還需進一步分析地形因子在土壤養分變化中的作用,以便對土壤養分空間格局機理進行合理的解釋。

環境因子對喀斯特森林土壤磷鉀養分空間分布特征和變異規律起著非常重要的作用[41]。研究區TP、TK、AP、AK均與海拔、坡向和地面凹凸度呈負相關關系；AK與土層深度極顯著負相關；TP、AK與坡度呈極顯著負相關關系；巖石裸露率與AP負相關。海拔、坡度和地面凹凸度對土壤磷鉀養分的影響最大,說明海拔越高、坡度越陡的點,土壤磷鉀養分的含量較低。這可能是由于淋溶作用強烈、植物的表聚作用以及凋落物相關。

以貴州為中心連帶成片的中國西南喀斯特地區(世界三大巖溶地區之一)受地球內動力、強烈的地質運動、高溫多雨且分布不均、碳酸鹽巖溶蝕性強、水文二維結構明顯以及其適生植物具有嗜鈣性、耐旱性和石生性等限制特點的綜合影響,生境和植被具有高度異質性,生態系統的穩定性和抗干擾性差,退化容易恢復難[42]。土壤特性空間變異研究可為土壤質量的恢復與改善,喀斯特脆弱生態系統植被的迅速恢復和生態重建提供理論依據。

4 結論

(1)喀斯特常綠落葉闊葉混交林土壤各磷鉀養分含量順序依次為AP>TK>AK>TP,變異程度均為中等變異。

(2)各磷鉀養分變異函數值的最佳擬合模型均為指數模型,決定系數均很高；TP、AP空間自相關中等,TK、AK為弱空間自相關。

(3)土壤TP表現為坡下(含洼地)含量高,坡上含量較低；AK表現為中坡含量高于洼地含量；AP、TK呈斑塊破碎化分布。

(4)海拔、坡度和地面凹凸度是影響各土壤磷鉀養分分布的主要因子。

[1] 王煒明. 基于GIS的地統計學方法在土壤科學中的應用. 中國農學通報, 2007, 23(5): 404- 408.

[2] 高艷霞, 徐東瑞. 石家莊市邊緣帶土壤養分空間分布特征研究. 土壤通報, 2009, 40(5): 1063- 1068.

[3] 徐劍波, 宋立生, 彭磊, 張橋. 土壤養分空間估測方法研究綜述. 生態環境學報, 2011, 20(8/9): 1379- 1386.

[4] 陸琦, 馬克明, 張潔瑜, 盧濤, 倪紅偉. 三江平原退化濕地和農田土壤養分的比較研究. 生態與農村環境學報, 2007, 23(2): 23- 28.

[5] Du H, Wang K L, Peng W X, Zeng F P, Song T Q, Zhang H, Lu S Y. Spatial heterogeneity of soil mineral oxide components in depression between karst hills, Southwest China. Chinese Geographical Science, 2014, 24(2): 163- 179.

[6] 劉淑娟, 張偉, 王克林, 陳洪松, 舒世燕, 譚衛寧. 桂西北喀斯特峰叢洼地表層土壤養分時空分異特征. 生態學報, 2011, 31(11): 3036- 3043.

[7] 宋同清, 彭晚霞, 杜虎, 王克林, 曾馥平. 中國西南喀斯特石漠化時空演變特征、發生機制與調控對策. 生態學報, 2014, 34(18): 5328- 5341.

[8] 陳效民, 吳華山, 孫靜紅. 太湖地區農田土壤中銨態氮和硝態氮的時空變異. 環境科學, 2006, 27(6): 1217- 1222.

[9] 張春霞, 郝明德, 王旭剛, 魏孝榮. 黃土高原溝壑區小流域土壤養分分布特征. 水土保持研究, 2003, 10(1): 78- 80.

[10] 孫洪泉. 地質統計學及其應用. 徐州: 中國礦業大學出版社, 1990.

[11] Li J W, Richter D D, Mendoza A, Heine P. Effects of land-use history on soil spatial heterogeneity of macro- and trace elements in the Southern Piedmont USA. Geoderma, 2010, 156(1/2): 60- 73.

[12] Loecke T D, Robertson G P. Soil resource heterogeneity in terms of litter aggregation promotes nitrous oxide fluxes and slows decomposition. Soil Biology and Biochemistry, 2009, 41(2): 228- 235.

[13] 龐夙, 李廷軒, 王永東, 余海英, 吳德勇. 土壤速效氮、磷、鉀含量空間變異特征及其影響因子. 植物營養與肥料學報, 2009, 15(1): 114- 120.

[14] 胡忠良, 潘根興, 李戀卿, 杜有新, 王新洲. 貴州喀斯特山區不同植被下土壤C、N、P含量和空間異質性. 生態學報, 2009, 29(8): 4187- 4195.

[15] 邢永強, 李金榮, 李金玲, 常秋玲, 賀傳閱. 土壤中硝態氮的空間變異研究. 灌溉排水學報, 2008, 27(3): 106- 108.

[16] 徐海燕, 聶宜民, 趙文武, 陳永智. 壽光市土壤養分特性及空間變異性分析. 農業現代化研究, 2009, 30(1): 118- 121.

[17] Cross A F, Schlesinger W H. Plant regulation of soil nutrient distribution in the northern Chihuahuan Desert. Plant Ecology, 1999, 145(1): 11- 25.

[18] Farley R A, Fitter A H. The responses of seven co-occurring woodland herbaceous perennials to localized nutrient-rich patches. Journal of Ecology, 1999, 87(5): 849- 859.

[19] 張珍明, 賀紅早, 張玉武, 劉盈盈, 何云松, 黃麗華. 雷公山自然保護區不同植被類型土壤的肥力及碳含量. 西南農業學報, 2014, 27(3): 1202- 1206.

[20] 張曉琳, 王帥, 王旭, 文關四, 劉海偉, 李佳靈. 海南吊羅山自然保護區土壤有機碳貯量研究. 熱帶作物學報, 2014, 35(2): 362- 368.

[21] 黃承標, 吳仁宏, 何斌, 曹雄亮, 黎家春, 黃必慶, 張傳, 李保平. 三匹虎自然保護區森林土壤理化性質的研究. 西部林業科學, 2009, 38(3): 16- 21.

[22] 宋同清, 彭晚霞, 曾馥平, 王克林, 覃文更, 譚衛寧, 劉璐, 杜虎, 鹿士楊. 木論喀斯特峰叢洼地森林群落空間格局及環境解釋. 植物生態學報, 2010, 34(3): 298- 308.

[23] 楊弘, 李忠, 裴鐵璠, 王安志, 金昌杰, 朱教君. 長白山北坡闊葉紅松林和暗針葉林的土壤水分物理性質. 應用生態學報, 2007, 18(2): 272- 276.

[24] 周莉, 代力民, 谷會巖, 于大炮. 長白山闊葉紅松林采伐跡地土壤養分含量動態研究. 應用生態學報, 2004, 15(10): 1771- 1775.

[25] 張偉, 陳洪松, 王克林, 張繼光, 侯婭. 典型喀斯特峰叢洼地坡面土壤養分空間變異性研究. 農業工程學報, 2008, 24(1): 68- 73.

[26] 張偉, 陳洪松, 王克林, 蘇以榮, 張繼光, 易愛軍. 喀斯特峰叢洼地土壤養分空間分異特征及影響因子分析. 中國農業科學, 2006, 39(9): 1828- 1835.

[27] 謝忠凱, 肖桂秀, 楊振玲, 宿艷萍, 潘永明. 長白山區新林地人參栽培土壤養分動態變化研究. 人參研究, 2006, 18(1): 10- 12.

[28] 楊小波, 胡榮桂. 熱帶濱海沙灘上森林植被的組成成分與土壤性質的研究. 生態學雜志, 2000, 19(4): 6- 11.

[29] 鮑士旦. 土壤農化分析. 北京: 中國農業出版社, 2000.

[30] 劉璐, 曾馥平, 宋同清, 彭晚霞, 王克林, 覃文更, 譚衛寧. 喀斯特木論自然保護區土壤養分的空間變異特征. 應用生態學報, 2010, 21(7): 1667- 1673.

[31] 張娜, 王希華, 鄭澤梅, 馬遵平, 楊慶松, 方曉峰, 謝玉彬. 浙江天童常綠闊葉林土壤的空間異質性及其與地形的關系. 應用生態學報, 2012, 23(9): 2361- 2369.

[32] 馮廣, 姚蘭, 艾訓儒, 黃繼紅, 路興慧, 丁易, 臧潤國. 鄂西南亞熱帶常綠落葉闊葉混交林土壤空間異質性. 生態學報, 2016, 36(23): 7635- 7643.

[33] 范夫靜, 宋同清, 黃國勤, 曾馥平, 彭晚霞, 杜虎, 鹿士楊, 時偉偉, 譚秋錦. 西南峽谷型喀斯特坡地土壤養分的空間變異特征. 應用生態學報, 2014, 25(1): 92- 98.

[34] 杜虎, 宋同清, 彭晚霞, 王克林, 劉璐, 鹿士楊, 曾馥平. 木論喀斯特自然保護區表層土壤礦物質的空間異質性. 農業工程學報, 2011, 27(6): 79- 84.

[35] 宋同清, 彭晚霞, 曾馥平, 歐陽資文, 吳海勇. 喀斯特木論自然保護區旱季土壤水分的空間異質性. 應用生態學報, 2009, 20(1): 98- 104.

[36] 張偉, 陳洪松, 王克林, 侯婭, 張繼光. 桂西北喀斯特洼地土壤有機碳和速效磷的空間變異. 生態學報, 2007, 27(12): 5168- 5175.

[37] 杜華強, 湯孟平, 崔瑞蕊. 天目山常綠闊葉林土壤養分的空間異質性. 浙江農林大學學報, 2011, 28(4): 562- 568.

[38] 劉叢強. 生物地球化學過程與地表物質循環——西南喀斯特土壤-植被系統生源要素循環. 北京: 科學出版社, 2009.

[39] 張偉, 王克林, 陳洪松, 張繼光. 典型喀斯特峰叢洼地土壤有機碳含量空間預測研究. 土壤學報, 2012, 49(3): 601- 606.

[40] 潘復靜, 張偉, 王克林, 何尋陽, 梁士楚, 韋國富. 典型喀斯特峰叢洼地植被群落凋落物C:N:P生態化學計量特征. 生態學報, 2011, 31(2): 335- 343.

[41] 劉淑娟, 張偉, 王克林, 陳洪松, 韋國富. 桂西北喀斯特峰叢洼地土壤物理性質的時空分異及成因. 應用生態學報, 2010, 21(9): 2249- 2256.

[42] 彭晚霞, 王克林, 宋同清, 曾馥平, 王久榮. 喀斯特脆弱生態系統復合退化控制與重建模式. 生態學報, 2008, 28(2): 811- 820.