蘆芽山亞高山草甸、云杉林土壤有機碳、全氮含量的小尺度空間異質(zhì)性

武小鋼，郭晉平，田旭平，楊秀云

(山西農(nóng)業(yè)大學(xué)林學(xué)院，太谷 030801)

生境的異質(zhì)性，尤其是土壤要素在空間上呈現(xiàn)復(fù)雜的鑲嵌性，與氣候以及陸地植被和生物發(fā)生復(fù)雜的相互作用，從而使得分析土壤的空間分布格局成為異質(zhì)性研究的一個重要領(lǐng)域［1-3］。不同植被類型下土壤由于承接其凋落物和根系分泌物類型的不同及氣候因子等的差異，因而形成的土壤碳、氮庫狀況存在差異［4-6］。以往的研究，多集中在陸地不同生態(tài)系統(tǒng)土壤碳、氮空間分布特征，及不同干擾和管理條件下土壤碳氮特征［1，5，7-8］。山地區(qū)域海拔高度的變化為研究生態(tài)系統(tǒng)過程的空間異質(zhì)性提供了條件，沿海拔梯度土壤碳氮變化特征的研究案例不斷有報道［9-11］。海拔梯度上土壤屬性空間異質(zhì)性研究多在較大尺度上進行［10-12］。空間異質(zhì)性是一個依賴于尺度的生態(tài)學(xué)概念，生態(tài)系統(tǒng)特性在不同尺度域上有著不同的變化速率，這種多尺度格局反映了生態(tài)系統(tǒng)的等級特征，指示著控制不同尺度格局的不同的生態(tài)學(xué)過程［13］。在等級關(guān)聯(lián)的生態(tài)系統(tǒng)中，小尺度上的空間異質(zhì)性研究可以為大尺度上的生態(tài)學(xué)格局與過程提供機制方面的解釋［14］。

因此，本研究利用地統(tǒng)計學(xué)的理論和方法，分析比較了山西蘆芽山不同海拔處分布的亞高山草甸和云杉林群落土壤有機碳和全氮的小尺度空間異質(zhì)性特征，旨在了解暖溫帶中部山區(qū)地帶性植被土壤碳氮的等級結(jié)構(gòu)特征，同時可以為不同尺度土壤的采樣設(shè)置提供理論依據(jù)，并希望有助于理解海拔對植被群落結(jié)構(gòu)和土壤碳、氮循環(huán)過程的影響，進而為不同植被類型的土壤碳匯管理技術(shù)研究提供基礎(chǔ)數(shù)據(jù)，為亞高山草甸和云杉林的合理利用和保護提供科學(xué)理論依據(jù)。

1 研究區(qū)概況

研究地點位于蘆芽山國家級自然保護區(qū)(38°36'—39°02'E，111°46'—112°54'N)，位于山西省呂梁山脈北端，黃土丘陵區(qū)的東部邊緣，山體由東北向西南斜向延伸，地勢高峻，最高峰荷葉坪海拔2772 m，是管涔山主峰。植被具有明顯的垂直地帶性，從高海拔到低海拔依次分布著亞高山草甸帶、云杉林帶、針闊葉混交林帶，灌草叢及農(nóng)墾帶［14］。

該區(qū)年均溫4.3—6.7 ℃，1月均溫－19.6 ℃，極端最低溫－36.6 ℃，7月均溫19.9 ℃，極端最高溫34.2℃，氣溫年較差和日較差大;年均降雨量453.9 mm，分布不均，6—9月降雨量約占全年的70%，降水年際變化大，歷年最大降水量為711.0 mm，最少降水量為252.9 mm;年蒸發(fā)量1800 mm，年均相對濕度50—55%;無霜期130—170 d。

2 研究方法

2.1 樣地基本情況

依據(jù)海拔高度和植被類型選擇不同海拔的亞高山草甸(A:2756.3 m;B:2542.3 m)2塊樣地和云杉林(C:2656.8 m;B:2387.2 m)2塊樣地。樣地基本情況見表1，針葉林每木檢尺測定樣地內(nèi)樹木的胸徑及位置，立木斷面積分布見圖1。

2.2 取樣方法

樣點布設(shè)依據(jù)地統(tǒng)計學(xué)理論和空間格局分析的小支撐、多樣點的取樣設(shè)計原則進行［15］。首先將樣地(30 m×30 m)等距離的劃分為100個3 m×3 m的小樣方。在大樣方內(nèi)選取45個樣點進行取樣。然后在樣地對角線上的兩個3 m×3 m樣方內(nèi)，分別設(shè)立100個小樣方(0.3 m×0.3 m的間隔距離)，從中各選取37個小樣方鉆取土樣(圖2)。土壤取樣用土鉆法進行，取樣時，先除去表層的枯枝落葉，然后在每個取樣點(共計119個)鉆取0—10 cm表層土壤樣品裝入塑料袋內(nèi)帶回實驗室分析測定。

表1 研究樣基本情況表Table1 General situation of research sites

圖1 樣地立木斷面積分布圖Fig.1 The map of basal area distribution in plots

2.3 分析測試方法

土壤全氮的測定用半微量凱氏定氮法，土壤有機碳的測定采用重鉻酸鉀氧化-外加熱法進行［16］。秤取風(fēng)干好的土壤樣品進行分析測定，每個土樣做3個重復(fù)，求其平均值作為土壤樣品全氮及有機碳含量值。

2.4 數(shù)據(jù)分析

2.4.1 經(jīng)典統(tǒng)計分析

用SPSSfor windows 18.0統(tǒng)計軟件進行土壤有機碳含量、全氮含量的平均數(shù)、標準差、變異系數(shù)分析。

2.4.2 異常值的識別、處理和原始數(shù)據(jù)的正態(tài)檢驗和轉(zhuǎn)換

進行特異值的判斷和處理。采用域法識別特異值，即樣本平均值(?)加減3倍的標準差(s)，在區(qū)間(?±3s)以外的數(shù)據(jù)為特異值，而后分別用正常的最大值或最小值來代替。

采用柯爾莫哥洛夫-斯米諾夫(Kolonogorov-Semirnov(K－S))正態(tài)性檢驗方法檢驗所測數(shù)據(jù)的正態(tài)分布，符合正態(tài)分布的數(shù)據(jù)直接進行地統(tǒng)計學(xué)分析(P(K－S)＞0.05);對于不符合正態(tài)分布的數(shù)據(jù)，要經(jīng)過對數(shù)轉(zhuǎn)換或方根轉(zhuǎn)換后再進行地統(tǒng)計學(xué)分析。

2.4.3 土壤有機碳和全氮含量的半方差函數(shù)模型分析地統(tǒng)計學(xué)分析用GS+Win5.0軟件進行。半方差函數(shù)用r(h)來表示，為區(qū)域化變量Z(xi)和Z(xi+h)增量平方的數(shù)學(xué)期望，即區(qū)域化變量的方差［17］。其通式為:

式中，r(h)為變異;h為步長，即為減少各樣點組合對的空間距離個數(shù)而對其進行分類的樣點空間間隔距離;N(h)為距離為h的點對的數(shù)量;Z(xi)和Z(xi+h)分別為變量Z在空間位置xi和xi+h的取值。

地統(tǒng)計學(xué)中的變異與經(jīng)典統(tǒng)計學(xué)中方差的根本差異在于變異考慮了空間尺度，即公式中的h。把植被的某一特征作為依賴變量，則變異能夠反映出統(tǒng)計意義上該變量在各個不同尺度的異質(zhì)性。

最優(yōu)模型的選擇，首先考慮決定系數(shù)(R2)和殘差平方和(RSS);殘差平方和RSS是對回歸模型進行顯著性檢驗的重要參數(shù)，其取值愈小，說明實際觀測值與回歸線靠近，擬合曲線與實際配合愈好。

分析理論模型參數(shù)，基臺值(Sill，C0+C)表示變量的最大變異程度，它的值越大表示變量的異質(zhì)性程度越高。而塊金值C0是空間距離為零時的變異值，表示隨機部分的空間變異性，較大的塊金值表明較小的尺度上某種生態(tài)學(xué)過程不容忽視。空間結(jié)構(gòu)比(Spatially structure variance)C/(C0+C)可度量空間自相關(guān)的變異所占的比例。塊金值(Nugget，C0)與基臺值之比C0/(C0+C)可用于估計隨機因素在所研究的空間異質(zhì)性中的相對重要性。

變程表示研究變量空間變異中空間自相關(guān)變異的尺度范圍，在變程內(nèi)，空間越靠近的點之間其相關(guān)性越大，距離大于變程的點之間不具備自相關(guān)性。

分維數(shù)可對不同變量之間的空間自相關(guān)強度進行比較，D值越大表示格局變異中隨機因素引起的異質(zhì)性的比重越大，D值越小，格局變異的空間依賴性越強。

3 結(jié)果與分析

3.1 土壤有機碳含量的空間變異性

3.1.1 土壤有機碳含量的描述性統(tǒng)計

對不同海拔高度的亞高山草甸和云杉林土壤表層(0—10 cm)有機碳含量的描述統(tǒng)計結(jié)果表明(表2)，樣地A(海拔2756.3 m)亞高山草甸土壤有機碳含量均值最高為49.84 g/kg，樣地B(海拔2542.3 m)的草甸土壤有機碳含量最低(38.33 g/kg);不同海拔云杉林樣地C(海拔2656.8 m)和樣地D(海拔2387.2 m)云杉純林土壤有機碳含量分別為47.06 g/kg和40.67 g/kg。從有機碳含量的均值比較分析，同一類型植被下，較高海拔的土壤有機碳含量高于較低海拔的土壤。從有機碳的變異分析，樣地A土壤有機碳的波動范圍23.76—67.31 g/kg，最大值是最小值的2.83倍;變異系數(shù)為15.81%。樣地B有機碳的波動范圍21.01—60.84 g/kg，最大值是最小值的2.89倍;變異系數(shù)(16.88%)大于樣地A。樣地C土壤有機碳的波動范圍21.34—83.81 g/kg，最大值是最小值的3.93倍;變異系數(shù)為28.62%。樣地D波動范圍17.03—92.39 g/kg，變異系數(shù)(32.06%)大于樣地C。以上結(jié)果表明，土壤有機碳的變異表現(xiàn)為云杉林有機碳的變異明顯大于亞高山草甸，同一植被類型，較低海拔土壤有機碳變異高于較高海拔的變異。

圖2 空間取樣設(shè)計Fig.2 Spatial sampling design

表2 土壤有機碳含量的描述性統(tǒng)計結(jié)果Table2 Statistics of SOC in the research area

3.1.2 土壤有機碳含量的半方差函數(shù)分析

土壤有機碳含量變異函數(shù)理論模型擬合結(jié)果及參數(shù)見表3。亞高山草甸和云杉林土壤有機碳含量的各向同性半方差理論模型為球狀模型。同一植被類型不同海拔土壤有機碳半方差函數(shù)的基臺值相比較，亞高山草甸樣地 B(2542.3m，C0+C=0.107)＞樣地 A(2756.3m，C0+C=0.046);云杉林樣地 D(2387.2m，C0+C=0.102)＞樣地C(2656.8m，C0+C=0.089)，總體上表現(xiàn)為較低海拔的樣地土壤有機碳含量的基臺值較高，進一步論證了上述變異系數(shù)所反映的同一植被類型下較低海拔土壤有機碳含量的異質(zhì)性較高的現(xiàn)象。

亞高山草甸樣地A土壤有機碳的空間結(jié)構(gòu)比(C/C0+C)為0.524，表現(xiàn)為中等強度的空間自相關(guān)變異特征。樣地B及云杉純林樣地C和樣地D土壤有機碳含量表現(xiàn)為強烈的空間自相關(guān)性。亞高山草甸樣地A和樣地B空間自相關(guān)的范圍為61.0 m。云杉純林樣地C和樣地D表現(xiàn)為小尺度的空間自相關(guān)變異(變程分別為 6.87，6.55 m)。

表3 土壤有機碳含量變異函數(shù)理論模型參數(shù)Table3 Parameters of semivariogram for SOC

用分維數(shù)D可對不同變量之間的空間自相關(guān)強度進行比較，本研究結(jié)果表明，土壤有機碳含量空間變異分維數(shù)值大小依次為樣地A＞樣地D＞樣地C＞樣地B。進一步表明樣地A空間變異中隨機因素引起的異質(zhì)性的比重大，而其他樣地格局變異的空間依賴性強。

3.2 土壤全氮含量的空間異質(zhì)性

3.2.1 土壤全氮含量的描述性統(tǒng)計

對不同海拔高度的亞高山草甸和云杉林土壤表層(0—10 cm)全氮含量的描述統(tǒng)計結(jié)果表明(表4)，草甸土壤全氮含量高于云杉純林;同一植被類型中，高海拔的樣地高于低海拔樣地。樣地A(海拔2756.3 m)土壤全氮含量均值為4.85 g/kg，全氮含量的波動范圍為2.89—10.10 g/kg，最大值是最小值3.49倍;土壤全氮含量的變異系數(shù)為24.33%。樣地B(海拔2542.3 m)土壤全氮含量均值為3.76 g/kg，全氮含量的波動范圍為2.48—5.09 g/kg，最大值是最小值2.05倍，土壤全氮含量的變異表現(xiàn)為低海拔(Cv=13.83%，樣地B)小于高海拔。樣地C(海拔2656.8 m)云杉純林土壤全氮量的均值為3.01 g/kg，全氮含量的變化范圍1.96—4.77 g/kg，最大值為最小值的2.43倍;變異系數(shù)為33.22%。樣地D(海拔2387.2 m)云杉純林土壤全氮量的均值為2.90 g/kg，全氮含量的變化范圍2.09—4.44 g/kg。云杉純林土壤全氮在不同海拔的變異與亞高山草甸相似，即較低海拔土壤全氮含量的變異(Cv=17.94%)小于較高海拔。分析表明森林土壤全氮含量的變異大于草甸土壤。

表4 土壤全氮含量的描述性統(tǒng)計結(jié)果Table 4 Statistics of TN in the research area

3.2.2 土壤全氮含量的地統(tǒng)計學(xué)分析

小尺度下所研究的4塊樣地土壤全氮含量的空間變異均表現(xiàn)為球狀模型的變化趨勢(表5)。同一植被類型不同海拔土壤有機碳半方差函數(shù)的基臺值相比較，亞高山草甸樣地A(2756.3m，C0+C=6.02)和樣地B(2542.3m，C0+C=1.041)遠大于云杉林樣地C和D，表明亞高山草甸土壤全氮含量的異質(zhì)性遠遠高于云杉林。而高海拔處的樣地A土壤全氮含量的空間異質(zhì)性也明顯高于低海拔處的樣地B。

所研究4個樣地的空間結(jié)構(gòu)比(C/C0+C)值均大于75%，說明土壤全氮含量在4個樣地種均表現(xiàn)出強的空間自相關(guān)性，其空間變異主要由結(jié)構(gòu)性因素造成的。亞高山草甸植被的空間自相關(guān)性明顯大于寒溫性針林樣地。亞高山草甸樣地A和樣地B空間自相關(guān)的范圍分別為47.65 m和61.0 m。云杉林樣地C和樣地D表現(xiàn)為小尺度的空間變異(變程分別為7.87，8.67 m)。D值大小依次為樣地C＞樣地D＞樣地B＞樣地A。亞高山草甸樣地A土壤全氮含量格局變異的空間依賴性較強。

表5 土壤全氮含量變異函數(shù)理論模型參數(shù)Table 5 Parameters of semivariogram for TN

4 討論

不同海拔高度處兩個亞高山草甸樣地(a:2756.6 m;B:2542.3 m)和兩個云杉林樣地(C:2656.8 m;D:2387.2 m)土壤有機碳和全氮含量的空間變異表現(xiàn)出很大不同。單因素方差分析結(jié)果表明(表6)，樣地A有機碳和全氮含量顯著高于樣地B(P＜0.001);樣地C和樣地D土壤全氮含量差異不顯著，而有機碳含量差異達極顯著水平(P＜0.001)。大量研究表明，海拔作為環(huán)境因子的綜合體現(xiàn)，通過對植被類型和植被生產(chǎn)力的制約直接影響輸入土壤的有機物質(zhì)量，通過土壤溫度和水分等條件影響微生物對有機質(zhì)的分解和轉(zhuǎn)化［10-11］。本研究中，樣地A土壤含水量顯著高于樣地B，樣地C土壤含水量顯著高于樣地D(表6);不考慮微地形對小氣候的影響，200—300 m的海拔落差會造成約1—2℃的氣溫差。因此，水熱要素的差異是不同海拔高度土壤有機碳和全氮含量空間異質(zhì)性產(chǎn)生的重要原因。

比較小尺度上亞高山草甸樣地和云杉林樣地的空間異質(zhì)性(表3，表5)，可以發(fā)現(xiàn)同一植被類型下較低海拔土壤有機碳含量的異質(zhì)性較高，其中隨機性因素對樣地A有機碳含量的空間變異影響較大。蘆芽山是山西省重要的夏季牧場，樣地A就位于林線以上亞高山草甸集中分布區(qū)(樣地B則位于林緣空地)，牲畜的啃食、踐踏及排泄等行為直接干擾草甸土壤環(huán)境，這成為樣地A土壤有機碳含量空間異質(zhì)性降低及隨機性變異比例增加的重要原因。對于云杉林樣地，通過比較樣地C和D的立木斷面積分布圖(圖1)不難看出，與樣地D相比，樣地C中林分密度大，立木空間分布較均一，較為一致的林下微環(huán)境可能是樣地C土壤有機碳含量空間異質(zhì)性較低的原因。植物群落的組成和群落中植物種群分布格局的改變制約著土壤組成的異質(zhì)化過程，同時決定著土壤養(yǎng)分循環(huán)［18-19］。

有機碳含量空間異質(zhì)性相反，樣地A的全氮含量表現(xiàn)為高異質(zhì)性空間分布，且空間自相關(guān)性強。在生態(tài)系統(tǒng)的物質(zhì)循環(huán)中，碳氮循環(huán)通過生產(chǎn)和分解緊密聯(lián)系在一起。考察4個樣地有機碳和全氮的相關(guān)性(圖3)可知，樣地B、C、D均呈現(xiàn)極顯著的正相關(guān)性，而樣地A相關(guān)性不顯著。進一步分析樣地A土壤含氮量數(shù)據(jù)，119個樣點中有8個高含氮量的(＞8.0)異常值，而平均值為4.85;如果去除這8個樣點數(shù)據(jù)，土壤有機碳和全氮含量相關(guān)性達極顯著水平。這與亞高山草甸和云杉林土壤有機碳含量和全氮含量極顯著正相關(guān)的研究結(jié)論相一致［20］。8個高含氮量樣點的坐標為(3，3)、(3，6)、3，12)、(3，18)、(3，24)、(3，27)、(6，3)、(6，6)，由空間取樣設(shè)計(圖 2)可知，這些樣點在空間上呈連續(xù)帶狀分布，暗示著樣地A存在一個氮源在地形和降雨的共同作用下遷移形成一條富氮帶。

圖3 土壤有機碳和全氮含量相關(guān)性Fig.3 Correlations between SOC and total nitrogen

不同海拔高度，相同植被類型下土壤有機碳和全氮含量的小尺度空間異質(zhì)性具有一定相似性。變異函數(shù)均呈球狀模型(表3，表5)，空間自相關(guān)變異的尺度大小相似。亞高山草甸土壤有機碳和全氮含量表現(xiàn)為較大尺度的空間自相關(guān)，而云杉林則表現(xiàn)為較小尺度的空間自相關(guān)。土壤屬性的空間分布是潛在的局地異質(zhì)性的總和，它們受生物學(xué)和地質(zhì)學(xué)等過程影響，使得區(qū)域化變量在空間分布上存在差異性，因而產(chǎn)生異質(zhì)的土壤環(huán)境［21］。研究結(jié)果反映出相同植被類型下，影響土壤有機質(zhì)和全氮的生態(tài)過程在相同的尺度上起作用;植被類型發(fā)生變化，則生態(tài)過程的尺度依賴性將發(fā)生顯著改變，而人為干擾如放牧將顯著改變碳氮循環(huán)的生態(tài)過程，表現(xiàn)為空間變異的增大和空間自相關(guān)性的下降。

5 結(jié)論

(1)亞高山草甸樣地A(海拔2756.3 m)和樣地B(海拔2542.3 m)土壤有機碳含量均值分別為49.84 g/kg和38.33 g/kg);云杉林樣地C(海拔2656.8 m)和樣地D(海拔2387.2 m)云杉純林土壤有機碳含量分別為47.06 g/kg和40.67 g/kg。有機碳空間異質(zhì)性總體上表現(xiàn)為，相同植被類型下較高海拔樣地有機碳含量高，而較低海拔的樣地土壤有機碳含量的異質(zhì)性較高。亞高山草甸樣地A和樣地B空間自相關(guān)的范圍為61.0 m。云杉純林樣地C和樣地D表現(xiàn)為小尺度的空間自相關(guān)變異(變程分別為6.87，6.55 m)。

(2)亞高山草甸樣地A和樣地B土壤全氮含量平均值為4.85 g/kg和3.76 g/kg，云杉樣地C和樣地D分別為3.01 g/kg和2.90 g/kg。高山草甸土壤全氮含量的異質(zhì)性遠遠高于云杉林。土壤全氮含量在4個樣地中均表現(xiàn)出強的空間自相關(guān)性，其空間變異主要由結(jié)構(gòu)性因素造成。亞高山草甸樣地和云杉林樣地全氮含量空間變異的尺度與有機碳含量空間變異表現(xiàn)相似，分別為47.65、61.0 m和7.87、8.67 m。

［1］Su SJ，Liu JF，He Z S，Hong W，Zhang JB.The spatial heterogeneity of soil nutrients in a mid-subtropical castanopsis kawakamii natural forest.Acta Ecologica Sinica，2012，32(18):5673-5682.

［2］Liu Y L，Zeng F P，Song T Q，Peng W X，Wang K L，Qin W G，Tan W N.Spatial heterogeneity of soil nutrients in karst area's mulun national nature reserve.Chinese Journal of Applied Ecology，2010，21(7):1667-1673.

［3］Pickett ST A，Cadenasso M.Landscape ecology:Spatial heterogeneity in ecological systems.Science，1995，269(5222):331-334.

［4］Xu G L，Mo JM，Zhou G Y，F(xiàn)u SL.Preliminary response of soil fauna to simulated n deposition in three typical subtropical forests.Pedosphere，2006，16(5):596-601.

［5］Conant R T，Six J，Paustian K.Land use effects on soil carbon fractions in the southeastern united states.I.Management-intensive versus extensive grazing.Biology and Fertility of Soils，2003，38(6):386-392.

［6］Peikun J.Soil active carbon pool under different types of vegetation.Scientia Silvae Sinicae，2005，41(1):10-13.

［7］Hungate B A，Dukes J S，Shaw M R，Luo Y，F(xiàn)ield C B.Nitrogen and climate change.Science，2003，302(5650):1512-1513.

［8］Bai J H，Deng W，Zhu Y M，Luan Z Q，Zhang Y X.Spatial distribution characteristics and ecological effects of carbon and nitrogen of soil in huolin river catchment wetland.Chinese Journal of Applied Ecology，2003，14(9):1494-1498.

［9］Xu X，Chen Y Q，Wang JS，F(xiàn)ang Y H，Quan W，Ruan H H，Xu Z K.Variations of soil labile organic carbon along an altitude gradient in Wuyi Mountain.Chinese Journal of Applied Ecology，2008，19(3):539-544.

［10］Zhang P，Zhang T，Chen N L.Vertical distribution patterns of soil organic carbon and total nitrogen and related affecting factors along northern slope of qilian mountains.Chinese Journal of Applied Ecology，2009，20(3):518-524.

［11］Xiang CH，Luan JW，Luo Z S，Gong Y B.Labile soil organic carbon distribution on influenced by vegetation types along an elevation gradient in West Sichuan，China.Acta Ecological Sinica，2010，30(4):1025-1034.

［12］Wang L，Ouyang H，Zhou CP，Zhang F，Song M H，Tian Y Q.Soil organic matter dynamics along a vertical vegetation gradient in the gongga mountain on the tibetan plateau.Journal of Integrative Plant Biology，2005，47(4):411-420.

［13］Cheng Y F，Dong M.Spatial heterogeneity in ecological systems.Acta Ecologica Sinica，2003，23(2):346-352.

［14］O'Neill R V，Deangelis D L.A Hierarchical Concept of Ecosystems.Princeton:Princeton University Press，1986.

［15］Mou P，Mitchell R J，Jones R H.Root distribution of two tree species under a heterogeneous nutrient environment.Journal of Applied Ecology，1997，34(3):645-656.

［16］Liu G S.Soil Physical and Chemical Analysis and Profile Description.Beijing:China Standards Press，1996.

［17］Wang Z Q.Geostatistic and Application in Ecology.Beijing:Science Press，1999.

［18］Wang H T，He X D，Gao Y B，Lu J G，Xue P P，Ma D.Density in artemisia ordosica successional communnity in response cto spatial heterogenity of soil moisture and organic matter.Chinese Journal of Plant Ecology，2007，31(6):1145-1153.

［19］Binkley D，Giardina C.Why do tree species affect soils?The warp and woof of tree-soil interactions.Biogeochemistry，1998，42(1/2):89-106.

［20］Wu X G，Guo JP，Yang X Y，Tian X P.Soil organic carbon storage and profile inventory in the different vegetation types of Luya Mountain.Acta Ecologica Sinica，2011，31(11):3009-3019.

［21］Gallardo A.Spatial variability of soil properties in a floodplain forest in Northwest Spain.Ecosystems，2003，6(6):564-576.

參考文獻:

［1］蘇松錦，劉金福，何中聲，洪偉，張金彪.格氏栲天然林土壤養(yǎng)分空間異質(zhì)性.生態(tài)學(xué)報，2012，32(18):5673-5682.

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

［8］白軍紅，鄧偉，朱顏明，欒兆擎，張玉霞.霍林河流域濕地土壤碳氮空間分布特征及生態(tài)效應(yīng).應(yīng)用生態(tài)學(xué)報，2003，14(9):1494-1498.

［9］徐俠，陳月琴，汪家社，方燕鴻，權(quán)偉，阮宏華，徐自坤.武夷山不同海拔高度土壤活性有機碳變化.應(yīng)用生態(tài)學(xué)報，2008，19(3):539-544.

［10］張鵬，張濤，陳年來.祁連山北麓山體垂直帶土壤碳氮分布特征及影響因素.應(yīng)用生態(tài)學(xué)報，2009，20(3):518-524.

［11］向成華，欒軍偉，駱宗詩，宮淵波.川西沿海拔梯度典型植被類型土壤活性有機碳分布.生態(tài)學(xué)報，2010，30(4):1025-1034.

［13］陳玉福，董鳴.生態(tài)學(xué)系統(tǒng)的空間異質(zhì)性.生態(tài)學(xué)報，2003，23(2):346-352.

［16］劉光崧.土壤理化分析與剖面描述.北京:中國標準出版社，1996.

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

［18］王海濤，何興東，高玉葆，盧建國，薛蘋蘋，馬迪.油蒿演替群落密度對土壤濕度和有機質(zhì)空間異質(zhì)性的響應(yīng).植物生態(tài)學(xué)報，2007，31(6):1145-1153.

［20］武小鋼，郭晉平，楊秀云，田旭平.蘆芽山典型植被土壤有機碳剖面分布特征及碳儲量.生態(tài)學(xué)報，2011，31(11):3009-3019.