王淑彬，徐慧芳，，宋同清 ，黃國勤，彭晚霞，杜 虎

(1．江西農(nóng)業(yè)大學作物生理生態(tài)與遺傳育種江西省/教育部重點實驗室，南昌 330045;2．中國科學院亞熱帶農(nóng)業(yè)生態(tài)研究所亞熱帶農(nóng)業(yè)生態(tài)過程重點實驗室，長沙 410125;

3．中國科學院環(huán)江喀斯特生態(tài)系統(tǒng)觀測研究站，環(huán)江 547100)

森林土壤是一個無比巨大的碳匯，也是一個天然的肥料庫，能有效地協(xié)調(diào)植物生長所需的水肥氣熱等條件，是影響植物生存最為重要的因素之一［1］。土壤是一個復雜的自然綜合體，內(nèi)部性質(zhì)存在著強烈的空間異質(zhì)性，即使在幾厘米的空間距離上其性質(zhì)也存在強烈變異［2］。因此，土壤屬性較大尺度的變化，對森林格局具有決定性作用［3］。土壤養(yǎng)分的空間分布特征為精確農(nóng)業(yè)及生態(tài)建模提供依據(jù)，幾乎所有的分布式水文模型［4］、氣候變化模型［5］都需要土壤養(yǎng)分的空間信息。土壤肥力是土壤最重要的生態(tài)功能之一，其空間布局與組成結(jié)構(gòu)直接影響著土壤生產(chǎn)力的高低、生態(tài)系統(tǒng)恢復的途徑與方向。在脆弱生境中，土壤N、P、K的空間分布特征與灌叢植被的分布區(qū)域具有高度相關(guān)性［6-7］。因此，開展森林土壤養(yǎng)分空間變異的研究對探索森植物群落結(jié)構(gòu)及其生態(tài)系統(tǒng)能量流動具有指導意義［8-9］。

廣西壯族自治區(qū)位于我國西南部，森林資源非常豐富，是中國南方重要林區(qū)之一。屬于低緯度地區(qū)，地處中、南亞熱帶季風氣候區(qū)，林種分布地域差異明顯，以北回歸線為界，南部為具有熱帶特點的森林，北部為亞熱帶常綠闊葉林。近年來，森林面積有了大幅度的增長，通過引入國外樹種，森林類型更加豐富，但缺乏廣西區(qū)森林土壤養(yǎng)分空間分布規(guī)律的科學認識，優(yōu)勢林種選擇的不合理導致森林質(zhì)量下降、土壤資源浪費，限制了廣西林業(yè)的發(fā)展。地統(tǒng)計學是在傳統(tǒng)統(tǒng)計學基礎(chǔ)上發(fā)展起來的空間分析方法，不僅能夠有效地揭示屬性變量在空間上的分布、變異和相關(guān)特征，而且可以將空間格局與生態(tài)過程聯(lián)系起來，有效解釋空間格局對生態(tài)過程與功能的影響［10-12］。國內(nèi)外學者廣泛應用地統(tǒng)計學研究了不同尺度上土壤屬性的空間變異特征［13-16］，喀斯特峰叢洼地森林土壤水分、有機質(zhì)、礦物質(zhì)、養(yǎng)分的研究也不少［17-21］，但關(guān)于區(qū)域尺度森林土壤肥力空間變異的報道甚少。本文以廣西區(qū)森林主要土壤養(yǎng)分為研究對象，用地統(tǒng)計學的空間特征和空間比較定量化方法，分析廣西區(qū)森林土壤養(yǎng)分的空間異質(zhì)性和分布格局，初步探討了其生態(tài)學過程和機制，為提高廣西區(qū)脆弱生態(tài)系統(tǒng)土壤肥力、促進該區(qū)域植被迅速恢復與生態(tài)重建提供理論依據(jù)和參考。

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

1．1 自然概況

廣西壯族自治區(qū)位于我國西南部，東經(jīng)104°26'—112°04'之間，北緯 20°54'—26°24'之間，屬熱帶、亞熱帶地區(qū)，氣候溫暖，日照充足，降水豐富，年降水量范圍為1080—2760 mm，降水季節(jié)分布不均，有明顯的干濕季，4—9月為雨季，總降水量可達到全年降水量的70%以上，10月至次年3月是干季。廣西地勢西北較高，由西北向東南傾斜。河流廣布，流向多與地質(zhì)構(gòu)造一致，四周多被山地、高原環(huán)繞，呈盆地狀。盆地邊緣多缺口，桂東北、桂東、桂南沿江一帶有大片谷地。廣西區(qū)內(nèi)土壤類型多樣，有18個土類，34個亞類，109個土屬，327個土種，主要有磚紅壤、赤紅壤、紅壤、石灰?guī)r土、黃壤、黃棕壤、紫色土、水稻土等［22］。

廣西地處亞熱帶熱帶地區(qū)，溫度從南到北由于緯度的升高而降低，直接影響著森林植被的分布，呈現(xiàn)規(guī)律性的更替顯現(xiàn)。廣西南部以北熱帶季節(jié)性雨林為主;中部與南部分別以棒科植物的常綠闊葉林和殼斗科植物為代表的常綠葉林。廣西區(qū)森林劃分為7個片區(qū)，天然林以常綠闊葉林為主，亞熱帶落葉闊葉林、亞熱帶針葉闊葉混交林、亞熱帶針葉林分布面積也較大。其中，亞熱帶針葉闊葉混交林僅在百色地區(qū)有分布。集中連片的天然闊葉林分布在九萬大山、大瑤山、海洋山、西大明山、貓兒山、富川西嶺、大明山、花坪林區(qū)、姑婆山等。廣西的人工林以松、杉、桉等用材林和油桐、油茶、八角、肉桂、栲膠等經(jīng)濟林為主。樣地的森林種類在各個地區(qū)屬于代表性林種，其中桉樹主要分布在桂東南地區(qū)，石山林主要分布于桂西北，杉木主要樣地主要分布在桂東北地區(qū)，松類樣地主要分布在桂南。

1．2 土樣采集方法與分析

在廣西區(qū)境內(nèi)設(shè)立共115個樣點，每個樣點有3個重復樣地(20 m×50 m)，共345個樣地，其中每個樣地之間直線距離＞100m，每個樣地圍成10個10m×10 m的小樣方，在每個樣地1、3、8小樣方用土鉆取樣法，去除表面枯枝落葉，取表層0—10 cm土壤，混合后土壤代表該樣點土樣，并詳細記錄該樣點的植被類型、巖石裸露率等環(huán)境信息，坡度、坡向、土層厚度等地理信息，GPS定位，記錄樣地四點經(jīng)緯度及海拔。

土壤樣品指標室內(nèi)分析主要包括全氮(TN)、全磷(TP)、全鉀(TK)、堿解氮(AN)、速效磷(AP)、速效鉀(AK)。其中TN采用半微量開氏法測定;TP采用NaOH熔融鉬銻抗顯色紫外分光光度法測定;TK采用NaOH熔融原子吸收法測定;AN采用堿解擴散法測定;AP采用0．05 mol/L NaHCO3提取鉬銻抗顯色紫外分光光度法測定;AK采用NH4Ac浸提原子吸收法測定。

1．3 數(shù)據(jù)處理

所用統(tǒng)計計算及模型的擬合在SPSS和專業(yè)地統(tǒng)計軟件GS+中完成。地統(tǒng)計學有關(guān)方法及原理限于篇幅不作詳細介紹，具體參見文獻［23-25］。

2 結(jié)果與分析

2．1 經(jīng)典統(tǒng)計描述

文中數(shù)據(jù)采用樣本均值加減3倍標準差識別特異值，在此區(qū)間外的數(shù)據(jù)均定為特異值，分別用正常的最大和最小值代替［26］，后續(xù)計算均采用處理后的原始數(shù)據(jù)。由廣西森林土壤養(yǎng)分含量的統(tǒng)計特征(表1)可看出，TN、TP、TK的變化范圍分別為0．27—6．85、0．1—1．46、0．41—30．75 g/kg，AN、AP、AK 分別為 34．67—561．19、0．93—78．76、7．23—417 mg/kg。不同土壤養(yǎng)分變異幅度不同，除AP外變異系數(shù)均介于10%—75%，呈中等變異，這可能與養(yǎng)分元素在土壤中的化學行為及肥料施用狀況、林間管理措施等有關(guān)。速效成分受隨機因素影響較大，所以變異程度較大，土壤全量變異系數(shù)較小，說明其在土壤中含量較穩(wěn)定。土壤各養(yǎng)分的偏態(tài)數(shù)(Skewness)均大于零，表明均呈正偏態(tài)分布，其中TK的偏度最小，AP的偏度最大。與標準正態(tài)分布相比，峰值系數(shù)(Kurtosis)大于3時樣本數(shù)據(jù)為高狹峰，低于3時為低闊峰。除AP外，其他養(yǎng)分均為低闊峰，TP最接近正態(tài)分布，其他養(yǎng)分均不服從正態(tài)分布，進行對數(shù)轉(zhuǎn)換后均符合正態(tài)分布。描述性統(tǒng)計分析能夠反映樣本全體的基本信息，而不能定量地刻畫土壤養(yǎng)分的隨機性、結(jié)構(gòu)性、獨立性和相關(guān)性［27］。因此，需要采用地統(tǒng)計學分析土壤養(yǎng)分的空間變異特征。

表1 廣西區(qū)森林土壤養(yǎng)分描述性統(tǒng)計特征Table 1 Descriptive statistics ofmain soil nutrients in Guangxi forests

2．2 空間變異特征

由圖1可知，廣西森林土壤AN試驗半變異函數(shù)的最佳擬合模型為球狀模型，其他養(yǎng)分均為指數(shù)模型，R2除TK 外在0．653—0．941之間，RSS除 TK 外均較小，說明各變量最佳模型擬合度較高，能很好地反映各土壤養(yǎng)分的空間結(jié)構(gòu)特征，各半變異函數(shù)模型的結(jié)構(gòu)參數(shù)見表2。TN、AN的變程相近且均較大，分別為 5．304°和 5．667°，TK 變程最小，為 0．186°，TP、AP、AK 的變程介于其間。TN、TP、AP 的塊金值/基臺值較大，介于28%—49．9%之間，為中等程度的空間自相關(guān)。TK、AN、AK的塊金值/基臺值相對較小，分別為 10．4%、23．7%、21．8%，都小于25%，表現(xiàn)為強烈的空間自相關(guān)性。各土壤養(yǎng)分的變異系數(shù)與塊金值/基臺值并不對應，AN變異程度雖然較大，但塊金值/基臺值較小，說明人類的隨機干擾對AN空間變異的貢獻較小，土壤氮素的空間異質(zhì)性主要來源于結(jié)構(gòu)因素。除AN、AK外各土壤養(yǎng)分的半變異函數(shù)曲線在超過一定滯后距后不再增加，而是圍繞基臺值呈周期性上下波動的特征，即孔穴效應［25，28］，這種現(xiàn)象主要是由區(qū)域化變量周期性變化引起的，說明研究區(qū)土壤養(yǎng)分空間異質(zhì)性具有周期性變化的特征。

表2 土壤養(yǎng)分全向半變異函數(shù)理論模型及其結(jié)構(gòu)參數(shù)Table 2 The best-fitted sem ivartiogram model of soil nutrients and the structure parameters

圖1 廣西森林土壤養(yǎng)分半變異函數(shù)圖Fig．1 Isotropic sem ivariogram of forest soil nutrients in Guangxi

2．3 空間分布特征

土壤養(yǎng)分的空間分布格局受海拔、人為干擾、植被、地形等人為與自然因子的綜合影響，不同因子的主導作用導致不同養(yǎng)分的空間格局。廣西森林土壤TN含量比較豐富，高于0．5 g/kg，空間變異規(guī)律明顯，呈北部地區(qū)含量較高、向南部和西北部逐漸降低的空間格局。TN含量高的森林土壤集中在西北部地區(qū)的河池、柳州、桂林、百色、賀州和崇左，森林類型主要為石山灌木及松杉類。AN的空間分布格局與TN極為相似(圖2)。

圖2 廣西森林土壤養(yǎng)分K riging等值線圖Fig．2 K riging contouring maps of forest soil nutrients in Guangxi

廣西森林土壤TP含量并不豐富，尤其廣西南部森林土壤TP極為貧乏;有兩個高峰區(qū)——桂柳和河池地區(qū)，含量大于1 g/kg，由這兩個區(qū)向四周逐漸減少，在廣西西南部百色與崇左出現(xiàn)了次高含量區(qū)域。AP的空間變異比較復雜，分布格局并沒有明顯規(guī)律，AP高值區(qū)分布在廣西區(qū)東北部與西南部邊緣，中等值集中在廣西西北、東北、西南3個邊角及位于南寧與百色交接處，含量最低的區(qū)域分布在柳州與梧州，這可能與河流的分布有關(guān)。

廣西森林土壤TK的空間分布比較均勻、平衡，高值分布在東南邊緣地區(qū)，次高值主要分布在東南部和東北部，中等值TK含量(10—15、5—10 g/kg)，占總面積的86．2%，說明廣西區(qū)森林土壤TK含量并不豐富。AK的空間變異格局比較有規(guī)律，總體的分布格局是北高南低，高值區(qū)主要分布在廣西東北部，次高值主要分布在東北部、西北部及中北部邊緣地區(qū)，中等值主要分布在廣西中部和北部大片區(qū)域，低含量區(qū)域主要分布在廣西南部。

3 結(jié)論和討論

3．1 廣西森林土壤養(yǎng)分的總體特征

廣西森林分為7個片區(qū)，不同片區(qū)養(yǎng)分含量狀況不同，變異程度也不相同?？傮w而言，北部森林養(yǎng)分含量高于南部森林片區(qū)(數(shù)據(jù)未顯示)，其中廣西南部ⅠA2a、ⅡA1a兩個片區(qū)各種養(yǎng)分含量均很低，造成片區(qū)養(yǎng)分不平衡的原因有多種，這可能與森林類型有關(guān)。南部片區(qū)森林類型主要為桉樹林，桉樹是廣西區(qū)主要人工經(jīng)濟林樹種，生長快，養(yǎng)分消耗大，導致土壤養(yǎng)分急劇流失，北部森林以石山灌木、松、杉類為主，生長速率較慢，養(yǎng)分循環(huán)也較好，受人為干擾小，土壤養(yǎng)分含量維持較高水平。土壤養(yǎng)分差異除受森林類型影響外，還與地形、氣候、土壤類型有密切關(guān)系。

廣西區(qū)森林土壤養(yǎng)分基本屬于中等變異，速效養(yǎng)分變異程度大于全量養(yǎng)分，其中AP、AK變異程度最大。進一步的地統(tǒng)計學分析，發(fā)現(xiàn)廣西森林土壤養(yǎng)分除TK外擬合效果均很好，變程大小為AK＞AN＞TN＞TP＞AP＞TK，TN、TP、AP 的塊金值/基臺值較大，表現(xiàn)為中等程度的空間自相關(guān)，說明在當前觀測尺度上，隨機因素對這些養(yǎng)分影響較大，與林間管理施肥、砍伐及實驗誤差有關(guān)。TK、AN、AK表現(xiàn)為強烈的空間自相關(guān)。說明隨機擾動對這些養(yǎng)分的影響相對較小，主要受研究區(qū)結(jié)構(gòu)性影響。

3．2 廣西森林土壤的空間變異和格局

廣西區(qū)森林土壤養(yǎng)分均表現(xiàn)出高度的空間異質(zhì)性，不同養(yǎng)分的空間變異特征和分布格局不同。廣西區(qū)森林土壤AN和TN含量水平都很高，表明廣西區(qū)氮素含量比較豐富。TN和AN空間變異特征相似度很高，呈廣西北部地區(qū)含量比較豐富、向南部和西北部地區(qū)逐漸降低的空間格局。TN含量與有機質(zhì)含量具有很強的相關(guān)性，有機質(zhì)也是TN的主要來源，廣西北部主要是石山林和松類，表層枯枝落葉量多，經(jīng)過長期的腐殖化作用，富集了大量的腐殖質(zhì)，有機質(zhì)含量豐富，而且廣西北部較南部寒冷，更加有利于氮素的積累［29-30］;南方主要為桉樹人工林，枯枝落葉比較少，有機質(zhì)含量低，同時桉樹的快速生長急劇地消耗土壤養(yǎng)分，從而影響了土壤氮素的含量［31］。廣西區(qū)氮素在植被類型、土壤類型、地形、土壤母質(zhì)、氣候及人類活動等因素綜合作用下形成了獨特的空間格局。

廣西森林土壤磷庫含量較小，TP的空間變異趨勢是從桂柳和河池兩個含量最高區(qū)向四周逐漸降低的趨勢，AP的空間變異比較復雜，規(guī)律不明顯。氣候的因素對土壤TP含量影響很大，在熱帶、亞熱帶地區(qū)氣候高溫多雨加快了土壤分化速度與磷元素的淋溶［32-33］。Miller等［34］認為在熱帶季風區(qū)隨著降水量的增加TP降低，相反，溫度低降水較少的地區(qū)TP含量比較高，與本研究結(jié)論一致，廣西北部氣溫較南部較低，降水較少，且北部喀斯特地貌廣布，地下河較多，降水快速滲透到地下河，地表水蓄積量少，不能充分參加與TP的理化反應。AP含量最低區(qū)主要分布在柳州與梧州兩個區(qū)域，這與河流的分布有關(guān)，這兩個區(qū)域均位于柳江流經(jīng)區(qū)域，河流匯聚細小河流，加速了速效磷的流失。不同土壤類型對AP的含量也有重要影響，不同地球化學類型上土壤有效磷密度的差異，反映了土壤中的物理、化學過程對土壤有效磷含量的控制［35-36］。AP空間格局分布復雜，與隨機因素有直接關(guān)系，不同林地的施肥管理措施不同，長期使用磷肥會提高土壤中AP的含量。

廣西森林土壤K含量并不豐富，TK和AK均屬中等水平，TK東南部含量稍高，AK呈北高南低。土壤TK、AK的變異性和其它養(yǎng)分一樣都受結(jié)構(gòu)性因素(土壤母質(zhì)、氣候、土壤類型、地形等)和隨機因素如不同的管理施肥措施等共同作用的結(jié)果。研究區(qū)AK與TK空間變異有一定的差異，說明兩者分別有自身形成的主導因素，土壤中新鮮植物殘體、根系、微生物中也含有一定量的鉀，但并不組成穩(wěn)定的含鉀有機物，因而對植物的有效性高。有機體一旦死亡，其所含鉀可被水淋洗或浸提出，也是土壤速效鉀的組成部分、補充液態(tài)鉀和交換性鉀的重要來源。由TK、AK空間格局可反映出廣西區(qū)鉀素總體成中等水平，在農(nóng)業(yè)活動中要針對鉀素的含量格局，科學地使用鉀肥來提高速效鉀的含量。

［1］ Huang C Y．Soil Science．Beijing:China Agriculture Press，2000:221-223．

［2］ Mallarino A P．Spatial variability patterns of phosphorus and potassium in no-tilled soils for two sampling scales．Soil Science Society of America Journal，1996，60(5):1473-1481．

［3］ Han Y Z，Wang Z Q．Spatial heterogeneity and forest regeneration．Chinese Journal of Applied Ecology，2002，13(5):615-619．

［4］ Abbott M B，Refsgaard JC．Distributed Hydrological Modeling．Netherland:Kluwer Academic Publishers，1996．

［5］ Zhao Y C，Shi X Z，Yu D S，Zhao Y F，Sun W X，Wang H J．Differentmethods for prediction of spatial patterns of soil organic carbon density in Hebei Province， China． Acta Pedologica Sinica，2005，42(3):379-385．

［6］ SchlesingerW H，Abrahams A D，Parsons A J，Wainwright J．Nutrient losses in runoff from grassland and shrubland habitats in southern New Mexico: I． rainfall simulation experiments．Biogeochemistry，1999，45(1):21-34

［7］ SchlesingerW H，Raikes JA，Hartley A E，Cross A F．On the spatial pattern of soil nutrients in desert ecosystems．Ecology，1996，77(2):364-374．

［8］ SchlesingerW H，Reynolds JF，Cunningham G L，Huenneke L F，Jarrell W M，Virginia R A，Whitford W G．Biological feedbacks in global desertification．Science，1990，247(4946):1043-1048．

［9］ Tilman D．Constraintsand tradeoffs:Toward a predictive theory of competition and succession．Oikos，1990，58:3-15．

［10］ Guo X D，F(xiàn)u B J，Ma K M，Chen L D，Yang F L．Spatial variability of soil nutrients based on geostatistics combined with GIS—A case study in Zunghua city of Hebai Province．Chinese Journal of Applied Ecology，2000，11(4):557-563．

［11］ Li B H，Zhang J，Yao X L，Ye J，Wang X G，Hao Z Q．Seasonal dynamics and spatial distribution patterns of herbs diversity in broadleaved Korean pine(Pinus koraiensis)mixed forest in Changbai Mountains．Chinese Journal of Applied Ecology，2008，19(3):467-473．

［12］ Ma FY，LiX R，Zhang JG，LiA X．Spatial heterogeneity of soil moisture in Shapotou sand-fixing artificial vegetation area．Chinese Journal of Applied Ecology，2006，17(5):789-795．

［13］ Ku?el S，N y'dl V，Kolá rˇL，Tich y'R．Spatial variability of cadmium，pH，organic matter in soil and its dependence on sampling scales．Water Air and Soil Pollution，1994，78(1/2):51-59．

［14］ Mishra T K，Banerjee S K．Spatial variability of soil pH and organic matter under Shorea robusta in lateritic region．Indian Journal of Forestry，1995，18(2):144-152．

［15］ Huang Y F，Zhou Z Y，Yuan X Y，Zhang H Y．Spatial variability of soil organic matter content in an arid desert area．Acta Ecologica Sinica，2004，24(12):2776-2781．

［16］ Su Y Z，Zhao H L，Cui J Y．Spatial heterogeneity of soil properties in the desertification process of rainfed farmland in Horqin sandy land．Acta Pedologica Sinica，2004，41(2):210-217．

［17］ PengW X，Song TQ，Zeng FP，Wang K L，Liu L，Du H，Lu S Y，Yin Q C．Spatial variation of soil water and the key impact factors in dry season in karst cluster-peak-depression region．Acta Ecologica Sinica，2010，30(24):6787-6797．

［18］ Song TQ，PengW X，Zeng F P，Ouyang ZW，Wu H Y．Spatial heterogeneity of surface soil moisture content in dry season in Mulun National Natural Reserve in Karst area．Chinese Journal of Applied Ecology，2009，20(1):98-104．

［19］ Ouyang ZW，PengW X，Song TQ，Zeng FP，Wang K L，Guan X，Wu H Y．Spatial heterogeneity of soil organic matter and its response to disturbance in karst peak cluster depressions．Chinese Journal of Applied Ecology，2009，20(6):1329-1336．

［20］ Du H，Song TQ，Wang K L，Liu L，Lu SY，Zeng F P．Spatial heterogeneity of mineral compositions in surface soil in Mulun National Nature Reserve karst areas．Transactions of the Chinese Society of Agricultural Engineering，2011，27(6):79-84．

［21］ Zhang W，Chen H S，Wang K L，Su Y R，Zhang JG，Yi A J．The heterogeneity of soil nutrients and their influencing factors in peak-cluster depression areas of karst region．Scientia Agricultura Sinica，2006，39(9):1828-1835．

［22］ Yu G Z．Random talk about main soils in Guangxi．Land and Resources of Southern China，2007，(3):39-40．

［23］ Li H B，Wang Z Q，Wang Q C．Theory and methodology of spatial heterogeneity quantification．Chinese Journal of Applied Ecology，1998，9(6):651-657．

［24］ Goovaerts P．Geostatistical tools for characterizing the spatial variability ofmicrobiological and physico-chemical soil properties．Biology and Fertility of Soils，1998，27(4):315-324．

［25］ Wang ZQ．Geostatistics and its application in Ecology．Beijing:Science press，1999．

［26］ Liu F C，Shi X Z，Yu D S，Pan X Z．Characteristics of spatial variability of total soil nitrogen in the typical area of Taihu Lake basin．Geographical Research，2004，23(1):63-70．

［27］ Tsegaye T，Hill R L．Intensive tillage effects on spatial variability of soil test，plantgrowth，and nutrientuptakemeasurements．Soil Science，1998，163(2):155-165．

［28］ Qin Y D．Quantitative analysis for soil spatial variability study．Advance in Earth Science，1992，7(1):44-49．

［29］ Trumbore S E，Chadwick O A，Amundson R．Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change．Science，1996，272(5260):393-396．

［30］ Wang L，Ouyang H，Zhou C P，Zhang F，Bai JH，Peng K．Distribution characteristics of soil organic matter and nitrogen on the eastern slope ofMt．Gongga．Acta Geographica Sinica，2004，59(6):1012-1018．

［31］ Zhang J J，Zhang Q，Yang Z P，Chen M C，Li Y，Chen D L．Spatial variability of topsoil organic matter and total nitrogen in Linfen Basin，Shanxi and its influencing factors．Chinese Journal of Soil Science，2010，41(4):839-844．

［32］ Onthong J，OsakiM，Nilnond C，Tadano T．Phosphorus status of some highly weathered soils in Peninsular Thailand and availability in relation to citrate and oxalate application．Soil Science and Plant Nutrition，1999，45(3):627-637．

［33］ Neufeldt H，da Silva J E，Ayarza M A，Zech W．Land-use Effects on phosphorus fractions in Cerrado oxisols．Biology and Fertility of Soils，2000，31(1):30-37．

［34］ Miller A J，Schuur E A G，Chadwick O A．Redox control of phosphorus pools in Hawaiian montane forest soils．Geoderma，2001，102(3/4):219-237．

［35］ Lajtha K，Schlesinger W H．The biogeochemistry of phosphorus cycling and phosphorus availability along a desert soil chronosequence．Ecology，1988，69(1):24-39．

［36］ Walbridge M R， Richardson C J， Swank W T． Vertical distribution of biological and geochemical phosphorus Subcycles in two southern Appalachian forest soils．Biogeochemistry，1991，13:61-85．

參考文獻:

［1］ 黃昌勇．土壤學．北京:中國農(nóng)業(yè)出版社，2000:221-223．

［3］ 韓有志，王政權(quán)．森林更新與空間異質(zhì)性．應用生態(tài)學報，2002，13(5):615-619．

［5］ 趙永存，史學正，于東升，趙彥鋒，孫維俠，王洪杰．不同方法預測河北省土壤有機碳密度空間分布特征的研究．土壤學報，2005，42(3):379-385．

［10］ 郭旭東，傅伯杰，馬克明，陳利頂，楊福林．基于GIS和地統(tǒng)計學的土壤養(yǎng)分空間變異特征研究——以河北省遵化市為例．應用生態(tài)學報，2000，11(4):557-563．

［11］ 李步杭，張健，姚曉琳，葉吉，王緒高，郝占慶．長白山闊葉紅松林草本植物多樣性季節(jié)動態(tài)及空間分布格局．應用生態(tài)學報，2008，19(3):467-473．

［12］ 馬風云，李新榮，張景光，李愛霞．沙坡頭人工固沙植被土壤水分空間異質(zhì)性．應用生態(tài)學報，2006，17(5):789-795．

［15］ 黃元仿，周志宇，苑小勇，張紅艷．干旱荒漠區(qū)土壤有機質(zhì)空間變異特征．生態(tài)學報，2004，24(12):2776-2781．

［16］ 蘇永中，趙哈林，崔建垣．農(nóng)田沙漠化演變中土壤性狀特征及其空間變異性分析．土壤學報，2004，41(2):210-217．

［17］ 彭晚霞，宋同清，曾馥平，王克林，劉璐，杜虎，鹿士楊，殷慶倉．喀斯特峰叢洼地旱季土壤水分的空間變化及主要影響因子．生態(tài)學報，2010，30(24):6787-6797．

［18］ 宋同清，彭晚霞，曾馥平，歐陽資文，吳海勇．喀斯特木論自然保護區(qū)旱季土壤水分的空間異質(zhì)性．應用生態(tài)學報，2009，20(1):98-104．

［19］ 歐陽資文，彭晚霞，宋同清，曾馥平，王克林，關(guān)欣，吳海勇．喀斯特峰叢洼地土壤有機質(zhì)的空間變化及其對干擾的響應．應用生態(tài)學報，2009，20(6):1329-1336．

［20］ 杜虎，宋同清，彭晚霞，王克林，劉璐，鹿士楊，曾馥平．木論喀斯特自然保護區(qū)表層土壤礦物質(zhì)的空間異質(zhì)性．農(nóng)業(yè)工程學報，2011，27(6):79-84．

［21］ 張偉，陳洪松，王克林，蘇以榮，張繼光，易愛軍．喀斯特峰叢洼地土壤養(yǎng)分空間分異特征及影響因子分析．中國農(nóng)業(yè)科學，2006，39(9):1828-1835．

［22］ 喻國忠．漫談廣西主要土壤．南方國土資源，2007，(3):39-40．

［23］ 李哈濱，王政權(quán)，王慶成．空間異質(zhì)性定量研究理論與方法．應用生態(tài)學報，1998，9(6):651-657．

［25］ 王政權(quán)．地統(tǒng)計學及在生態(tài)學中的應用．北京:科學出版社，1999．

［26］ 劉付程，史學正，于東升，潘賢章．太湖流域典型地區(qū)土壤全氮的空間變異特征．地理研究，2004，23(1):63-70．

［28］ 秦耀東．土壤空間變異研究中的定量分析．地球科學進展，1992，7(1):44-49．

［30］ 王琳，歐陽華，周才平，張鋒，白軍紅，彭奎．貢嘎山東坡土壤有機質(zhì)及氮素分布特征．地理學報，2004，59(6):1012-1019．

［31］ 張建杰，張強，楊治平，陳明唱，李勇，陳德利．山西臨汾盆地土壤有機質(zhì)和全氮的空間變異特征及其影響因素．土壤通報，2010，41(4):839-844．