桂西北喀斯特地區(qū)石灰土養(yǎng)分空間變異特征

俞月鳳,何鐵光,*,宋同清,李麗娟,韋彩會(huì),蒙炎成,唐紅琴,李忠義,李婷婷,胡 芳

1 廣西壯族自治區(qū)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所, 南寧 530007 2 中國(guó)科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所, 長(zhǎng)沙 410125 3 廣西桂江林業(yè)調(diào)查規(guī)劃設(shè)計(jì)有限公司,南寧 530022

土壤是在母質(zhì)、氣候、生物、地形及人為因素共同作用下形成的,這些因素在時(shí)空上具有不一致性,使土壤物理、化學(xué)及生物學(xué)性質(zhì)也都表現(xiàn)出普遍的空間變異性[1- 3]。這種空間變異性在不同區(qū)域、不同尺度以及不同時(shí)間段上表現(xiàn)出普遍性和復(fù)雜性[4]。土壤養(yǎng)分作為土壤資源中最基本的屬性之一,其豐缺程度及存在形態(tài)直接決定著土壤的肥力狀況,影響著植被的生長(zhǎng)與演替[5]。土壤養(yǎng)分的空間分布格局是土壤空間異質(zhì)性的具體表現(xiàn)[6- 7]。國(guó)內(nèi)外大量學(xué)者對(duì)不同尺度和不同區(qū)域的土壤養(yǎng)分進(jìn)行了空間異質(zhì)性研究[8- 11]?？梢娡寥鲤B(yǎng)分異質(zhì)性深刻影響到管理者對(duì)土壤利用相應(yīng)策略和政策的制定,同時(shí)也是精準(zhǔn)農(nóng)業(yè)和精確施肥以及生態(tài)建模的信息基礎(chǔ)[12- 13]。

桂西北喀斯特地區(qū)是中國(guó)貧困與環(huán)境退化問題最為突出的地區(qū)之一,面臨環(huán)境退化、人口密集、經(jīng)濟(jì)與社會(huì)落后等多重難題。在喀斯特環(huán)境中,石灰?guī)r風(fēng)化作用和土壤形成過程是以土壤為媒介的生物地球化學(xué)過程[14]。喀斯特生態(tài)系統(tǒng)是以碳酸鹽為物質(zhì)基礎(chǔ),具有土壤淺薄、土被不連續(xù)及土壤pH值高等特征的脆弱生態(tài)系統(tǒng)[15]。正是這種特殊的形成基質(zhì)和生態(tài)環(huán)境條件,使得喀斯特地區(qū)土壤具有高度的空間異質(zhì)性[16-17]。土壤養(yǎng)分作為生態(tài)系統(tǒng)結(jié)構(gòu)與功能的重要因子,對(duì)喀斯特生態(tài)系統(tǒng)的穩(wěn)定起著決定作用[18]。在生態(tài)脆弱地區(qū),土壤養(yǎng)分的空間分布特征直接影響該地區(qū)土壤生產(chǎn)力的高低和生態(tài)恢復(fù)的途徑和方向,因而探明土壤養(yǎng)分空間變異是合理指導(dǎo)喀斯特地區(qū)農(nóng)業(yè)生產(chǎn)和生態(tài)重建的關(guān)鍵[17]。由于空間變異尺度效應(yīng)的存在,研究尺度的選擇成為土壤養(yǎng)分空間變異研究的關(guān)鍵一步[4]。就喀斯特地區(qū)土壤養(yǎng)分空間異質(zhì)性而言,小尺度的研究較多,主要集中在森林、坡地及不同土地利用方式下的土壤[19-23]。小尺度區(qū)域土壤養(yǎng)分空間變異的研究能夠直接指導(dǎo)生產(chǎn)實(shí)踐,但大尺度區(qū)域的研究也必不可少,其能為決策者提供更廣闊的視角和更完整的認(rèn)識(shí)[4]。然而,桂西北喀斯特區(qū)域尺度上,石灰土養(yǎng)分空間異質(zhì)性有待進(jìn)一步研究,從而為桂西北喀斯特地區(qū)石漠化綜合治理中植被恢復(fù)、生態(tài)重建等宏觀決策提供理論依據(jù)。因此,本研究從大尺度區(qū)域出發(fā),將整個(gè)桂西北喀斯特地區(qū)作為研究區(qū)域,以大面積野外調(diào)查采樣和室內(nèi)分析測(cè)定為基礎(chǔ),采用經(jīng)典統(tǒng)計(jì)學(xué)、地統(tǒng)計(jì)學(xué)等方法,探索喀斯特區(qū)域內(nèi)石灰土養(yǎng)分的含量水平、空間變異特性、空間分布特征,通過克里格插值繪制整個(gè)桂西北區(qū)域各養(yǎng)分指標(biāo)的空間分布圖。

1 研究區(qū)概況和研究方法

1.1 研究區(qū)概況

桂西北喀斯特區(qū)位于廣西壯族自治區(qū)西北部,屬云貴高原山前區(qū),地處104°28′—109°09′E,22°51′—25°37′N,包括河池和百色2地級(jí)市23縣(市、區(qū))。區(qū)域內(nèi)以山地、丘陵為主,巖溶地貌廣泛發(fā)育,裸露石灰?guī)r面積約3.38 萬km2,占全區(qū)域的48.5%(河池地區(qū)石山面積占本地區(qū)的65.7%,百色地區(qū)石山面積占本地的31.2%)。區(qū)域以典型的峰叢洼地巖溶地貌為主,海拔85.0—2062.0 m。地處亞熱帶季風(fēng)氣候區(qū),光照充足,年均日照時(shí)數(shù)達(dá)1968.6 h,無霜期315.3 d；氣候溫和,年均氣溫20.3℃,≥10℃積溫7130.4℃；降水集中于6—8月份,年均降雨量1178.8 mm,年均暴雨日數(shù)4.6 d,最大降雨強(qiáng)度67.8 mm/h；干濕季明顯,年均蒸發(fā)量1298.7 mm,干旱日數(shù)181.2 d,干旱指數(shù)1.10?？λ固氐貐^(qū)土壤類型主要以石灰土為主,有黑色石灰土、紅色石灰土、黃色石灰土、棕色石灰土4個(gè)亞類。

1.2 樣品采集與分析

在對(duì)未知區(qū)域進(jìn)行初探時(shí),覆蓋整個(gè)研究區(qū)的均勻網(wǎng)格布點(diǎn)法是較為常用的方法[4]。然而,面積廣大的研究區(qū)域以及喀斯特復(fù)雜多變的地形地貌使得在整個(gè)桂西北喀斯特地區(qū)完成規(guī)則的網(wǎng)格布點(diǎn)非常困難?；诳λ固胤鍏餐莸鼐坝^異質(zhì)性強(qiáng)的特點(diǎn),根據(jù)裸露石灰?guī)r分布特征及面積,確定調(diào)查網(wǎng)格面積為100 km2(約為0.1°經(jīng)緯網(wǎng)格),網(wǎng)格數(shù)量為338個(gè),參照《IPCC優(yōu)良做法指南》對(duì)系統(tǒng)隨機(jī)抽樣的建議,調(diào)查網(wǎng)格數(shù)量為24個(gè),占總網(wǎng)格的7%,按每個(gè)網(wǎng)格約布置7—8個(gè)樣點(diǎn)計(jì),整個(gè)桂西北喀斯特區(qū)域設(shè)置175個(gè)樣點(diǎn)。在實(shí)際采樣過程中,根據(jù)石灰土的分布圖及道路系統(tǒng),綜合考慮植被類型及采樣點(diǎn)的可達(dá)性和實(shí)際可操作性,采樣點(diǎn)的選擇原則是要具有隨機(jī)性并且兼顧均勻分布的采樣方法[24]。所有采樣點(diǎn)分為三種植被類型：草灌、灌叢、次生林,研究區(qū)域樣點(diǎn)分布見圖1。

土壤樣品采集于2015年11月,用手持GPS記錄采樣中心點(diǎn)坐標(biāo),采用“之”字形取樣法,在以該點(diǎn)為中心的20 m×20 m的范圍內(nèi)采集0—20 cm深度的5個(gè)土壤副樣本,四分法混合為1個(gè)土壤主樣本保留1 kg左右。同時(shí)記錄海拔高度、植被種類、植被蓋度、坡度、坡向等信息。將采集的土壤樣品去除石頭、雜草和植物根系等非土壤成分后置于通風(fēng)處自然風(fēng)干,磨碎過篩,保存于自封袋中待測(cè)。土壤樣品室內(nèi)分析指標(biāo)為有機(jī)質(zhì)SOM(重鉻酸鉀-外加熱法)、全氮TN(半微量開氏法)、全磷TP(NaOH熔融-鉬銻抗顯色-紫外分光光度法)、全鉀TK(NaOH熔融-火焰光度計(jì)法)、堿解氮AN(堿解-擴(kuò)散法)、速效磷AP(0.5 mol/L NaHCO3浸提法)、AK速效鉀(NH4OAc浸提法)和pH值(電極電位法)。

圖1 桂西北喀斯特區(qū)域及采樣點(diǎn)分布圖Fig.1 Map of the karst area in northwestern Guangxi and distribution of soil sampling sites

1.3 數(shù)據(jù)處理

文中數(shù)據(jù)采用域法識(shí)別特異值,即樣本均值加減3倍標(biāo)準(zhǔn)差識(shí)別特異值,在此區(qū)間外的數(shù)據(jù)均定為特異值,分別用正常的最大和最小值代替[25]。土壤養(yǎng)分的描述性統(tǒng)計(jì)分析、K-S檢驗(yàn)(α =0.05)和普通相關(guān)分析在SPSS 18.0軟件中完成；相關(guān)分析中數(shù)值變量之間使用皮爾遜相關(guān)系數(shù),分類變量之間以及與數(shù)值變量之間使用斯皮爾曼秩相關(guān)系數(shù)。地統(tǒng)計(jì)學(xué)分析在GS+軟件中完成,用理論模型對(duì)各土壤養(yǎng)分指標(biāo)的半方差函數(shù)進(jìn)行擬合。Kriging等值線圖在ArcGIS 10.0中完成。

2 結(jié)果與分析

2.1 石灰土養(yǎng)分的描述性統(tǒng)計(jì)分析

石灰土養(yǎng)分的描述性統(tǒng)計(jì)特征如表1所示。桂西北喀斯特地區(qū)石灰土pH值呈微酸性至微堿性性,為5.65—8.12,平均值為7.04,呈中性；石灰土SOM和養(yǎng)分含量均較高,但其變異相差很大。根據(jù)Nielsen標(biāo)準(zhǔn)[26],除了pH呈弱變異外(CV<10%),其他指標(biāo)均呈中度變異(30%

表1 石灰土養(yǎng)分描述性統(tǒng)計(jì)特征

SOM：有機(jī)質(zhì),Soil organic matter；AN：堿解氮,Available nitrogen；AK：速效鉀,Available potassium；AP：速效磷,Available phosphorus；TN：全氮,Total nitrogen；TP：全磷,Total phosphorus；TK：全鉀,Total potassium

2.2 石灰土養(yǎng)分的空間變異特征

采用GS+軟件進(jìn)行各養(yǎng)分指標(biāo)半方差函數(shù)模型擬合,獲得半方差函數(shù)圖(圖2)和結(jié)構(gòu)參數(shù)(表2)。如果變量在采樣尺度上具有空間依賴性,半變異函數(shù)會(huì)隨著滯后距的增加而增大,并且在超過一定距離后逐漸趨近于基臺(tái)值或圍繞基臺(tái)值上下波動(dòng)[27]。由圖2可以看出,AK的最佳擬合模型為線性模型,SOM的最佳擬合模型為球狀模型,AP的最佳擬合模型為高斯模型,其他指標(biāo)的最佳擬合模型均為指數(shù)模型。決定系數(shù)(R2)或殘差平方和(RSS)用于判斷模型擬合的好壞。由表2的擬合參數(shù)可知,AP的高斯模型和SOM的球狀模型擬合效果最好,其他養(yǎng)分的擬合模型也可以較好地反映桂西北喀斯特石灰土養(yǎng)分的空間結(jié)構(gòu)特征。

圖2 石灰土養(yǎng)分的半方差函數(shù)圖Fig.2 Semi-variance grams of limestone soil nutrientsSOM：有機(jī)質(zhì),Soil organic matter；AN：堿解氮,Available nitrogen；AK：速效鉀,Available potassium；AP：速效磷,Available phosphorus；TN：全氮,Total nitrogen；TP：全磷,Total phosphorus；TK：全鉀,Total potassium

各土壤指標(biāo)的塊金值均為正值,說明本研究存在一定的采樣和試驗(yàn)誤差。土壤指標(biāo)均具有一定的基臺(tái)值,說明其受地形、坡位、土壤結(jié)構(gòu)和母質(zhì)等自然因素所控制。塊金值與基臺(tái)值的比值[C0/(C0+C)]反映了土壤屬性的空間依賴性,一般認(rèn)為<25%時(shí),空間變量為強(qiáng)烈的空間自相關(guān),在25%—75% 之間時(shí),為中等空間自相關(guān),>75%為弱空間自相關(guān)[22]。SOM、AP、TN和TP 的C0/(C0+C)值在50%—55%之間,表現(xiàn)為中等空間自相關(guān)性,說明空間差異是由隨機(jī)因素和結(jié)構(gòu)性因素共同作用引起的；pH、AN、AK和TK的C0/(C0+C)值均>75%,表現(xiàn)為弱空間自相關(guān)性,表明空間變異主要受隨機(jī)因素的控制,而結(jié)構(gòu)性因素對(duì)變異的貢獻(xiàn)很小,同時(shí)也說明對(duì)這些指標(biāo)進(jìn)行研究應(yīng)適當(dāng)減小采樣距離。A 為變程,空間變量只有在變程以內(nèi)具有空間自相關(guān)性。石灰土養(yǎng)分指標(biāo)變程大小各異,范圍在59.4—379.5 km之間,TN和SOM變程最大,說明較大的范圍內(nèi)存在空間相關(guān)性,其空間連續(xù)性較好,具有較差的空間漸變性規(guī)律；AK、AP和TP次之,pH、AN和TK的變程最小,為50—70 km之間,說明其空間自相關(guān)的范圍較小。

2.3 桂西北喀斯特地區(qū)石灰土養(yǎng)分的空間分布格局

喀斯特地區(qū)石灰土養(yǎng)分含量的空間插值圖(圖3)所示,在圖中可以按照顏色深淺分為8個(gè)含量等級(jí),能夠直觀清晰地反映出研究區(qū)域石灰土養(yǎng)分的空間分布特征。SOM、AN、AP、TN和TP分布較為相似,高值區(qū)主要分布在西北部,由西向東逐漸遞減；其中SOM和TN斑塊較大,變化緩和,最低值區(qū)均出現(xiàn)在東部偏南地區(qū)；AN在東部地區(qū)也出現(xiàn)了高值區(qū),中值區(qū)分布的斑塊最大；pH值的分布呈中部高,沿著四周降低的趨勢(shì),高值區(qū)斑塊相對(duì)較小,同時(shí)可以看出pH值主要分布在6.92—7.43之間；AK高值區(qū)出現(xiàn)在東北部,斑塊小,連續(xù)性差；TK的分布無規(guī)律性,空間連續(xù)性差,破碎化程度高,在小尺度上形成強(qiáng)烈的變異。

2.4 石灰土養(yǎng)分空間異質(zhì)性的影響因素分析

為了初步探明桂西北喀斯特區(qū)石灰土養(yǎng)分空間變異特征的影響因素,選取了海拔、坡度、巖石出露率和植被類型等指標(biāo),利用相關(guān)分析對(duì)上述各環(huán)境因子與石灰土養(yǎng)分之間的關(guān)系進(jìn)行了分析。表3列出了根據(jù)原始數(shù)據(jù)計(jì)算出來的各變量之間的相關(guān)系數(shù)。區(qū)域尺度上,大多數(shù)變量之間都存在不同程度的線性相關(guān)關(guān)系。例如,SOM與石灰土養(yǎng)分TN、AN、TP、AK、pH之間有極顯著的正相關(guān)關(guān)系,與AP、TK的相關(guān)關(guān)系并不明顯。 海拔與SOM、AN、TN、TP、裸巖率、坡度和植被類型之間有極顯著的正相關(guān)關(guān)系,與TK正相關(guān)顯著。裸巖率與AK、TK相關(guān)關(guān)系不明顯,而與其他石灰土養(yǎng)分指標(biāo)及環(huán)境因子之間均存在極顯著的正相關(guān)關(guān)系。坡度與SOM、AN、TN、AK、植被類型之間變現(xiàn)為顯著或極顯著的正相關(guān)關(guān)系,與AP為極顯著負(fù)相關(guān)關(guān)系。植被類型與SOM、AN、TN、AP之間為顯著正相關(guān)關(guān)系。

表2 石灰土養(yǎng)分的半方差函數(shù)模型參數(shù)

圖3 石灰土養(yǎng)分的空間分布Fig.3 Kriging maps showing spatial distribution of limestone soil nutrients

指標(biāo)IndexSOMANAKAPTNTPTKPHHLYLPDZBSOM1AN0.507??1AK0.251??0.256??1AP0.021-0.0490.335??1TN0.609??0.589??0.304??0.0561TP0.351??0.337??0.369??0.473??0.680??1TK0.024-0.0220.216??0.0060.171?0.206??1PH0.238??0.276??0.347??-0.0210.343??0.342??0.183?1H0.405??0.220??0.0490.0420.292??0.226??0.159?-0.0021LYL0.415??0.254??0.0590.241??0.454??0.211??0.0480.322??0.198??1PD0.337??0.191?0.233??-0.443??0.296??0.0190.0740.1640.238??0.617??1ZB0.436??0.482??0.1080.201?0.495??0.0490.0590.1480.1200.304??0.394??1

*顯著相關(guān)(P<0.05),**極顯著相關(guān)(P<0.01)； H：海拔Altitude；LYL：裸巖率Rock exposed rate；PD：坡度Slope； ZB：植被類型Vegetation types

3 討論

土壤養(yǎng)分空間異質(zhì)性受到氣候(降雨、太陽(yáng)輻射等)、地形(坡度、坡向、海拔等)、土壤物理性質(zhì)(機(jī)械組成、結(jié)構(gòu)、孔隙等)以及植被類型等因素所影響?？λ固氐貐^(qū)地形地貌復(fù)雜,生境的異質(zhì)性程度高,地形的變化導(dǎo)致了環(huán)境因子的空間變異,使得區(qū)域內(nèi)土壤養(yǎng)分具有明顯的空間變異性。地形是重要的土壤成土因子之一,它能調(diào)控著太陽(yáng)輻射和降水的空間分配,影響局部生境的小氣候與土壤厚度和養(yǎng)分的空間差異[28-29]。本研究中,石灰土養(yǎng)分與地形因子的相關(guān)分析結(jié)果表明,石灰土養(yǎng)分的空間分布與地形特征有密切的關(guān)系。其中,SOM、AN、TN與海拔、坡度、裸巖率和植被類型極顯著的正相關(guān)關(guān)系。土壤碳素和氮素在生物地球化學(xué)過程中緊密聯(lián)系在一起。據(jù)有關(guān)研究發(fā)現(xiàn),海拔高度每上升100 m土壤有機(jī)碳濃度會(huì)增加0.75—2.10 mg/g[30]。海拔高度和坡度對(duì)土壤養(yǎng)分的影響通常是間接的,是通過對(duì)其他自然因素的影響而產(chǎn)生作用。石灰土AK和TK分別與坡度和海拔有極顯著和顯著的正相關(guān)關(guān)系,而與其他環(huán)境因子相關(guān)不明顯。自然情況下,土壤鉀素基本全部由母質(zhì)風(fēng)化而來,其分布特征可能與土壤的淋溶特征有關(guān)[21]。石灰土AP與坡度、裸巖率、植被類型等顯著相關(guān),而TP與海拔和裸巖率有極顯著的正相關(guān)關(guān)系,有關(guān)研究發(fā)現(xiàn),植被對(duì)磷素利用能力的差異可能是導(dǎo)致土壤速效磷素資源空間分配異質(zhì)較大的主要原因[18]。在石灰土條件下,難溶性無機(jī)磷極易與堿性土壤中的鈣結(jié)合生成Ca10-P 沉淀,不易淋失[31]。石灰土pH值和裸巖率有極顯著的正相關(guān)關(guān)系,這與張忠華等[19]的研究結(jié)果相似。因此,裸巖率高的地方,土壤pH及鈣含量增加,土壤TP有較明顯的洼積效應(yīng)。在區(qū)域大尺度下,石灰土養(yǎng)分含量隨著海拔高度、坡度和裸巖率的增加而增加,這與其他地區(qū)明顯不同。張偉等[32]的研究結(jié)果表明立地因子較高(裸巖率大、坡度較陡、土被較薄) 處的土壤養(yǎng)分含量較高,本文的研究結(jié)果與其基本一致。這是由于在強(qiáng)烈的巖溶作用下,微地貌形態(tài)的發(fā)育對(duì)土壤養(yǎng)分積累和循環(huán)具有重要的影響,在坡度較陡、坡位較高的地段容易形成石槽、石洞等陰暗生境,有利于有機(jī)質(zhì)的積累,土壤養(yǎng)分含量也相對(duì)較高[32]。植被類型對(duì)石灰土養(yǎng)分分布特征有著重要的影響,旺盛的生物累積過程是土壤養(yǎng)分形成和維持的基礎(chǔ),而土壤養(yǎng)分的空間變異又作用于植物的生長(zhǎng)和空間分布。

在區(qū)域大尺度上,桂西北喀斯特地區(qū)石灰土SOM、TN、TP和AP呈現(xiàn)較為相似的分布規(guī)律,養(yǎng)分含量呈自西向東逐漸減少的趨勢(shì),空間變異由結(jié)構(gòu)因素和隨機(jī)因素共同影響。這種分布特點(diǎn)與區(qū)域尺度上大地形區(qū)的分布相關(guān)：桂西北喀斯特地區(qū)北與貴州接壤,西與云南毗鄰,海拔高度自東向西,自南向北逐漸升高。海拔差異導(dǎo)致水熱條件及其組合在空間上的變化,并常伴隨著溫度、降水、光照、土壤和植被類型等許多因子的改變,海拔梯度的大小充分反映了環(huán)境間的差異。因此,海拔、坡度、裸巖率等地形因子可能是桂西北區(qū)域尺度上石灰土養(yǎng)分空間分布變異的主要結(jié)構(gòu)因素。當(dāng)然,在區(qū)域范圍內(nèi)影響石灰土空間異質(zhì)性的環(huán)境因素還有很多,比如氣溫、降雨、土壤物理和生物特征等。由于本研究區(qū)域內(nèi)收集到的氣象站數(shù)據(jù)較少,未能較準(zhǔn)確提取各個(gè)采樣點(diǎn)的氣象數(shù)據(jù),因此在分析過程沒有引入氣候因素。今后,需更加完善各方面因素及土壤指標(biāo)的研究,揭示這些因素對(duì)石灰土養(yǎng)分空間變異的影響。

4 結(jié)論

桂西北喀斯特區(qū)域尺度下,石灰土養(yǎng)分具有不同的變異特征,變異程度依次為AP>TP>TN>SOM>AK>AN>TK>pH,除pH為弱變異外,其他均是中等強(qiáng)度變異。SOM和TN具有較為相似的變異特征,變程分別為338.8 km和379.5 km；AK的變異尺度相對(duì)較小,為213.0 km；AP和TP變異尺度近乎相等,變程為142.1和179.2；pH、AN和TK的變異程度最小,變程范圍僅在59.4—70.0 km之間。石灰土SOM、AP、TN和TP表現(xiàn)為中等程度的空間自相關(guān),分布有一定規(guī)律性,連續(xù)性較好,這是由結(jié)構(gòu)性和隨機(jī)性因素共同作用的結(jié)果；pH、AN、AK和TK表現(xiàn)為弱空間相關(guān)性,在區(qū)域尺度上的分布沒有呈現(xiàn)一定的規(guī)律性,連續(xù)性差,斑塊破碎化程度較高,說明這些空間變異主要來自隨機(jī)因素,同時(shí)也說明對(duì)這類養(yǎng)分指標(biāo)進(jìn)行研究時(shí)應(yīng)適當(dāng)減小土壤采樣點(diǎn)的間距。桂西北喀斯特地區(qū)石灰土養(yǎng)分存在不同程度的空間變異性,且空間自相關(guān)范圍各不相同,說明石灰土養(yǎng)分的空間異質(zhì)性對(duì)尺度的依賴程度各不相同。因此采樣距離需要根據(jù)各養(yǎng)分指標(biāo)的特征進(jìn)行調(diào)整。相關(guān)分析結(jié)果表明,地形因子(如海拔、坡度和巖石裸露率等)和植被類型是影響桂西北喀斯特石灰土空間分布特征和變異規(guī)律的重要因素。

參考文獻(xiàn)(References):

[1] Jenny H. Factors of Soil Formation——A System of Quantitative Pedology. New York: Dover Publications, 1941.

[2] Burrough P A. Soil variability: a late 20th century view. Soil and Fertilizers, 1993, 56(5): 529- 562.

[3] Heuvelink G B M, Webster R. Modelling soil variation: past, present, and future. Geoderma, 2001, 100(3/4): 269- 301.

[4] 劉志鵬. 黃土高原地區(qū)土壤養(yǎng)分的空間分布及其影響因素[D]. 咸陽(yáng): 中國(guó)科學(xué)院研究生院(教育部水土保持與生態(tài)環(huán)境研究中心), 2013.

[5] 黃昌勇. 土壤學(xué). 北京: 中國(guó)農(nóng)業(yè)出版社, 2000.

[6] Huggett R J. Soil chronosequences, soil development, and soil evolution: a critical review. Catena, 1998, 32(3/4): 155- 172.

[7] 王國(guó)梁, 劉國(guó)彬, 許明祥. 黃土丘陵區(qū)紙坊溝流域植被恢復(fù)的土壤養(yǎng)分效應(yīng). 水土保持通報(bào), 2002, 22(1): 1- 5.

[8] Tesfahunegn G B, Tamene L, Vlek P L G. Catchment-scale spatial variability of soil properties and implications on site-specific soil management in northern Ethiopia. Soil and Tillage Research, 2011, 117(6): 124- 139.

[9] Marchetti A, Piccini C, Francaviglia R, Mabit L. Spatial distribution of soil organic matter using geostatistics: a key indicator to assess soil degradation status in central Italy. Pedosphere, 2012, 22(2): 230- 242.

[10] 黃紹文, 金繼運(yùn), 楊俐蘋, 程明芳. 縣級(jí)區(qū)域糧田土壤養(yǎng)分空間變異與分區(qū)管理技術(shù)研究. 土壤學(xué)報(bào), 2003, 40(1): 79- 88.

[11] 楊艷麗, 史學(xué)正, 于東升, 王洪杰, 徐茂, 王果. 區(qū)域尺度土壤養(yǎng)分空間變異及其影響因素研究. 地理科學(xué), 2008, 28(6): 788- 792.

[12] 陳文輝, 謝高地, 卓慶卿. 農(nóng)田基礎(chǔ)環(huán)境信息空間變異性分析. 生態(tài)學(xué)報(bào), 2004, 24(2): 347- 351.

[13] 趙永存, 史學(xué)正, 于東升, 趙彥鋒, 孫維俠, 王洪杰. 不同方法預(yù)測(cè)河北省土壤有機(jī)碳密度空間分布特征的研究. 土壤學(xué)報(bào), 2005, 42(3): 379- 385.

[14] 潘根興, 曹建華. 表層帶巖溶作用: 以土壤為媒介的地球表層生態(tài)系統(tǒng)過程——以桂林峰叢洼地巖溶系統(tǒng)為例. 中國(guó)巖溶, 1999, 18(4): 287- 296.

[15] 李陽(yáng)兵, 王世杰, 李瑞玲. 巖溶生態(tài)系統(tǒng)的土壤. 生態(tài)環(huán)境, 2004, 13(3): 434- 438.

[16] 王世杰, 盧紅梅, 周運(yùn)超, 謝麗萍, 肖德安. 茂蘭喀斯特原始森林土壤有機(jī)碳的空間變異性與代表性土樣采集方法. 土壤學(xué)報(bào), 2007, 44(3): 475- 483.

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

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

[19] 張忠華, 胡剛, 祝介東, 倪健. 喀斯特森林土壤養(yǎng)分的空間異質(zhì)性及其對(duì)樹種分布的影響. 植物生態(tài)學(xué)報(bào), 2011, 35(10): 1038- 1049.

[20] 張偉, 劉淑娟, 葉瑩瑩, 陳洪松, 王克林, 韋國(guó)富. 典型喀斯特林地土壤養(yǎng)分空間變異的影響因素. 農(nóng)業(yè)工程學(xué)報(bào), 2013, 29(1): 93- 101.

[21] 劉璐, 曾馥平, 宋同清, 彭晚霞, 王克林, 覃文更, 譚衛(wèi)寧. 喀斯特木論自然保護(hù)區(qū)土壤養(yǎng)分的空間變異特征. 應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(7): 1667- 1673.

[22] 范夫靜, 宋同清, 黃國(guó)勤, 曾馥平, 彭晚霞, 杜虎, 鹿士楊, 時(shí)偉偉, 譚秋錦. 西南峽谷型喀斯特坡地土壤養(yǎng)分的空間變異特征. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(1): 92- 98.

[23] 許聯(lián)芳, 王克林, 朱捍華, 侯婭, 張偉. 桂西北喀斯特移民區(qū)土地利用方式對(duì)土壤養(yǎng)分的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2008, 19(5): 1013- 1018.

[24] 李龍, 姚云峰, 秦富倉(cāng), 高玉寒, 張美麗. 黃花甸子流域土壤全氮含量空間分布及其影響因素. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(5): 1306- 1312.

[25] 宋同清, 彭晚霞, 曾馥平, 歐陽(yáng)資文, 吳海勇. 喀斯特木論自然保護(hù)區(qū)旱季土壤水分的空間異質(zhì)性. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(1): 98- 104.

[26] Nielsen D R, Bouma J. Soil Spatial Variability. Wageningen: Pudoc, 1985.

[27] Ettema C H, Wardle D A. Spatial soil ecology. Trends in Ecology & Evolution, 2002, 17(4): 177- 183.

[28] Enoki T, Kawaguchi H, Iwatsubo G. Topographic variations of soil properties and stand structure in aPinusthunbergiiplantation. Ecological Research, 1996, 11(3): 299- 309.

[29] Tateno R, Takeda H. Forest structure and tree species distribution in relation to topography-mediated heterogeneity of soil nitrogen and light at the forest floor. Ecological Research, 2003, 18(5): 559- 571.

[30] Leifeld J, Bassin S, Fuhrer J. Carbon stocks in Swiss agricultural soils predicted by land-use, soil characteristics, and altitude. Agriculture, Ecosystems & Environment, 2005, 105(1/2): 255- 266.

[31] 張素霞. 黃土高原坡地不同土地利用方式下土壤剖面磷素分布及其有效性研究. 楊凌: 西北農(nóng)林科技大學(xué), 2008.

[32] Zhang W, Chen H S, Wang K L, Su Y R, Zhang J G, Yi A J. The heterogeneity and its influencing factors of soil nutrients in Peak-Cluster Depression Areas of Karst Region. Agricultural Sciences in China, 2007, 6(3): 322- 329.