高麗楠　張宏

摘要：應(yīng)用經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)相集合的方法，研究青藏高原東緣典型的高寒草甸、高寒沼澤樣地中小尺度下土壤活性鐵的空間變異特征。結(jié)果表明，研究區(qū)土壤活性鐵變異系數(shù)隨著樣地尺度的增加而增大，變異系數(shù)為142%～48.8%，均屬中等程度變異。半方差分析結(jié)果，高寒草甸3個(gè)尺度的樣地從小到大M1、M2、M3的最佳擬合模型依次為高斯、指數(shù)、球狀模型；高寒沼澤S1、S2、S3依次為指數(shù)、球狀、指數(shù)模型。高寒草甸和高寒沼澤的土壤活性鐵的空間變異性的塊金值/基臺(tái)值比值[C0/（C0+C1）]均介于50%～68%，具有中等強(qiáng)度的空間相關(guān)性，結(jié)構(gòu)性因素和隨機(jī)性因素共同對(duì)變異起作用?？死锝鸩逯祱D顯示高寒草甸土壤活性鐵明顯呈帶狀和斑塊狀分布，空間連續(xù)性較好，南部活性鐵含量明顯高于北部；高寒沼澤土壤活性鐵呈現(xiàn)斑狀分布，插值圖比較破碎。在相同尺度下高寒沼澤群落土壤活性鐵的分形維數(shù)值要小于高寒草甸群落，高寒沼澤群落土壤活性鐵呈現(xiàn)更強(qiáng)的空間異質(zhì)性。

關(guān)鍵詞：土壤活性鐵；地統(tǒng)計(jì)學(xué)；空間變異

中圖分類號(hào)： S153.6+1文獻(xiàn)標(biāo)志碼： A

文章編號(hào)：1002-1302（2017）15-0239-04

地殼和土壤中鐵元素的豐度僅次于氧、硅、鋁，不僅影響土壤結(jié)構(gòu)，還可直接以化學(xué)途徑或間接通過微生物/生物途徑干預(yù)大量元素循環(huán)，從而影響土壤功能[1]。特別是對(duì)濕地土壤的氧化還原環(huán)境來說，土壤鐵的變化具有非常重要的指示性意義[2]。同時(shí)，鐵可以直接影響植物的生長發(fā)育，土壤對(duì)植物營養(yǎng)元素的供給情況可以通過土壤中鐵的含量及分布來反映[3]。其中，活性鐵是在一定條件下能夠被活化的氧化鐵，在調(diào)節(jié)Fe3+和 Fe2+含量上具有重要作用，是土壤對(duì)植物供鐵的重要指標(biāo)[1]。所以在植物營養(yǎng)學(xué)和環(huán)境化學(xué)等領(lǐng)域中，土壤活性鐵均起著重要的作用。

青藏高原是地球上最高的高原，平均海拔在4 000 m以上，地域廣闊，面積為2 572.4×103 km2，占我國陸地總面積的26.8%[4]，被譽(yù)為地球的“第三極”。青藏高原具有高海拔、強(qiáng)太陽輻射、空氣稀薄、低氣溫等地理氣候特征，使得高原植被和土壤對(duì)全球變化極為敏感，因此，被認(rèn)為是研究全球變化的“敏感區(qū)”。青藏高原土壤元素的自然含量及其特征的研究，是該區(qū)開發(fā)和規(guī)劃中必不可少的資料，對(duì)于了解該區(qū)的環(huán)境地球化學(xué)與元素化學(xué)地理上具有重要意義。土壤環(huán)境的重要屬性之一是土壤鐵的空間異質(zhì)性[5]。土壤鐵的空間異質(zhì)性的研究具有相當(dāng)重要的意義，不僅可以理解土壤的形成過程、結(jié)構(gòu)和功能，而且可以了解植物和土壤之間的關(guān)系，如土壤的更新過程，植物的空間分布格局，對(duì)土壤養(yǎng)分和水分變化對(duì)根系產(chǎn)生的影響也具有非常重要意義。自20世紀(jì)90年代以來，土壤鐵的空間異質(zhì)性成為生態(tài)學(xué)研究的重點(diǎn)問題，相關(guān)學(xué)者開始關(guān)注土壤鐵含量的空間異質(zhì)性特征[6-7]。本研究在中小尺度下通過高密度網(wǎng)格嵌套采樣，從不同群落類型、樣地大小、采樣密度等角度探討土壤活性鐵的空間異質(zhì)性變化，揭示青藏高原東緣高寒地區(qū)不同草地類型土壤活性鐵的空間分布及其變異特征，有助于認(rèn)識(shí)該地區(qū)土壤活性鐵在不同空間范圍上的變異特征，最終為大尺度的土壤活性鐵的分布、估算模型模擬和真實(shí)信息獲取提供實(shí)踐基礎(chǔ)，對(duì)評(píng)價(jià)青藏高原在全球環(huán)境變化下生物地球化學(xué)循環(huán)的響應(yīng)有重要的意義。

1材料與方法

1.1試驗(yàn)地概況

研究區(qū)位于青藏高原東南邊緣的四川省阿壩藏族羌族自治州紅原縣（101°51′～103°23′E，31°51′～33°19′N），平均海拔3 600 m。該區(qū)氣候?yàn)榇箨懶愿咴瓪夂?，年均氣?.6～11 ℃，最暖為7月，平均溫度為10.8 ℃，最冷為1月，平均溫度為-10.6 ℃；區(qū)內(nèi)年均雨量781.75 mm，降水主要集中在5—9月；日照時(shí)間長，太陽輻射強(qiáng)，太陽輻射年總量為 6 194 MJ/m2。

1.2樣品采集與分析

本研究在紅原縣龍日種畜場(chǎng)的草場(chǎng)設(shè)置樣地，分別選取淺丘山地草甸（102°22.91′E，32°27.21′N，海拔3 563 m）、丘間低地沼澤（102°22.10′E，32°26.73′N，海拔3 560 m）2塊樣地。前者以矮生嵩草（Kobresia humilis）、翻白委陵菜（Potentilla discolor）、鵝絨委陵菜（Potentilla anserine）、乳白香青（Anaphalis lactea）、米口袋（Gueldenstaedtia verna）、火絨草（Leontopodium lentopodioides）、剪股穎（Agrostis matsumurae）等為群落優(yōu)勢(shì)種；后者以烏拉草（Carex meyeriana）、金蓮花（Trollius chinensis）、小白花地榆（Sanguiorba parviflora）等為群落優(yōu)勢(shì)種。同時(shí)草場(chǎng)處于圍欄禁牧的時(shí)期。

在植被生長盛期的2015年7月中旬，分別對(duì)高寒草甸樣地和高寒沼澤樣地進(jìn)行土壤取樣。在2塊樣地內(nèi)，均采用大樣地套小樣地方法設(shè)置10 m×5 m 、5 m×2.5 m、2 m×1 m樣地尺度。在3種樣地尺度下分別按0.5 m×0.5 m、0.25 m×0.25 m、0.1 m×0.1 m面積劃分成200個(gè)小樣方，用網(wǎng)格法布設(shè)土壤活性鐵采樣樣點(diǎn)。采用蛇形順序在各樣地內(nèi)取樣，取樣深度均為0～20 cm，將編號(hào)裝入布袋的土樣帶回實(shí)驗(yàn)室，經(jīng)風(fēng)干、研磨、過篩，制備裝瓶備用。

土壤活性鐵的測(cè)定采用原子吸收光度計(jì)法測(cè)定。土壤中活性氧化物采用酸性草酸-草酸銨溶液提取，在不損壞土壤礦質(zhì)部分時(shí)，使氫氧化鐵等凝膠進(jìn)入提取液。然后用鹽酸羥胺為還原劑，將三價(jià)鐵還原為二價(jià)鐵，二價(jià)鐵與鄰啡啰林反應(yīng)生成橙紅色的絡(luò)合物進(jìn)行比色。利用UVWIN5紫外-可見光分光光度計(jì)進(jìn)行測(cè)量。

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

采用地統(tǒng)計(jì)學(xué)方法進(jìn)行分析。地統(tǒng)計(jì)學(xué)中，區(qū)域化變量在一定尺度上的相關(guān)程度和空間變異可以通過半方差函數(shù)的塊金值、基臺(tái)值等參數(shù)來反映。相關(guān)統(tǒng)計(jì)學(xué)的方法和原理可見參考文獻(xiàn)[8]。描述性統(tǒng)計(jì)利用SPSS 12.0 完成，變異函數(shù)計(jì)算和理論模型的擬合均利用ArcGIS 10.1完成。endprint

2結(jié)果與分析

2.1土壤活性鐵經(jīng)典統(tǒng)計(jì)

從表1 可以看出，研究區(qū)高寒草甸土壤活性鐵含量為704～779 g/kg，高寒沼澤土壤活性鐵含量為13.08～15.71 g/kg，2種群落下土壤活性鐵含量差異明顯。隨機(jī)變量離散程度的大小通過變異系數(shù)（CV）來表示，即土壤活性鐵空間變異性的大小。CV<10%，表現(xiàn)空間上為弱變異性；10%100%，說明空間上屬于強(qiáng)變異性[9]。在不同植被覆蓋、不同采樣尺度下，土壤活性鐵變異系數(shù)在14.2%～48.8%，均屬中等變異，且變異系數(shù)隨樣地尺度的增加而增大。進(jìn)行地統(tǒng)計(jì)分析的前提是檢驗(yàn)數(shù)據(jù)是否符合正態(tài)分布，因原始數(shù)據(jù)存在正偏態(tài)效應(yīng)，所以對(duì)數(shù)據(jù)進(jìn)行自然對(duì)數(shù)轉(zhuǎn)換，其結(jié)果符合正態(tài)分布要求，滿足地統(tǒng)計(jì)學(xué)分析要求。

2.2不同采樣尺度下土壤活性鐵的空間結(jié)構(gòu)特征

半方差函數(shù)參數(shù)既能反映出分析對(duì)象空間相關(guān)的類型又能說明空間分布的結(jié)構(gòu)，同時(shí)能夠量化空間相關(guān)的范圍，其中，塊金值（C0） 反映隨機(jī)性變異，受到試驗(yàn)誤差和采樣尺度及系統(tǒng)屬性本身變異特征所控制；偏基臺(tái)值（C1） 反映由結(jié)構(gòu)性因素（地形、氣候等）引起的變異；基臺(tái)值（C0+C1） 表示系統(tǒng)內(nèi)部總的變異，包括隨機(jī)性變異和結(jié)構(gòu)性變異，值的大小表示總變異程度的高低，基臺(tái)值則表示系統(tǒng)總的空間異質(zhì)性，其值越大則空間異質(zhì)性越高[10]。隨機(jī)因素引起的變異與總變異的關(guān)系通過塊基比[C0 /（C0+C1）]能夠得到反映，當(dāng)塊基比較高時(shí)，表明引起空間變異的主要作用為隨機(jī)部分；當(dāng)比值較低時(shí)，則表明引起空間變異的主導(dǎo)因素為空間自相關(guān)部分[11-12]。在高寒草甸樣地，采樣尺度土壤活性鐵含量半方差函數(shù)最佳擬合模型M1符合高斯模型，M2符合指數(shù)模型，而M3符合球狀模型。從表2 可以看出，塊金值、基臺(tái)值和塊基比均隨采樣尺度增大先減小后增大，M1、M2、M3樣地的塊基比分別為63%、32%、35%。對(duì)于高寒草甸來說，M1樣地土壤活性鐵結(jié)構(gòu)性因子影響占的比重小，為37%，而隨機(jī)性因子的影響則較大；M2、M3樣地的結(jié)構(gòu)性因子影響分別占68%、65%，以結(jié)構(gòu)性因子的影響為主。高寒草甸采樣尺度M1表現(xiàn)出各向同性，而M2、M3表現(xiàn)出各向異性，各向異性比均為6.51。

在高寒沼澤S2樣地，土壤活性鐵半方差函數(shù)最佳擬合模型符合球狀模型，S1、S3符合指數(shù)模型。塊金值、基臺(tái)值和塊基比均隨著采樣尺度的增大而增大，S1、S2、S3樣地的塊基比分別為17%、50%、68%。高寒沼澤樣地采樣尺度S1、S2均表現(xiàn)出各向異性，各向異性比分別為2.06、2.32，但采樣尺度S3表現(xiàn)出各向同性。高寒沼澤S1樣地土壤活性鐵呈現(xiàn)較好的空間相關(guān)性，結(jié)構(gòu)性因子影響的空間異質(zhì)性明顯高于S1、S2樣地尺度。高寒沼澤土壤活性鐵隨機(jī)部分的空間異質(zhì)性隨樣地尺度的增大而逐漸升高，逐漸掩蓋小尺度的結(jié)構(gòu)特征。

C0 /（C0+C1）<25%，表明變量空間相關(guān)性較強(qiáng)；C0 /（C0+C1） 介于25%～75%，則說明變量具有中等的空間相關(guān)性；C0 /（C0+C1）>75%，則表明變量的空間相關(guān)性很弱[13]。本研究區(qū)內(nèi)的高寒沼澤和高寒草甸土壤活性鐵在不同采樣尺度下塊基比值均較大，幾乎在25%和75%范圍之間，則表示中等強(qiáng)度的空間變異性，其空間變異是由隨機(jī)性因素和結(jié)構(gòu)性因素共同作用。僅有高寒沼澤樣地S1的塊基比值小于25%，主要是由結(jié)構(gòu)性因素（地形、氣候等）引起的，且空間相關(guān)性較為強(qiáng)烈。

2.3不同采樣尺度下土壤活性鐵的分形維數(shù)

通過半方差函數(shù)的雙對(duì)數(shù)曲線的斜率可以求得分形維數(shù)，其代表空間異質(zhì)性程度，如果空間異質(zhì)性程度越強(qiáng)，則分形維數(shù)的數(shù)值越小，反之則增大，當(dāng)分形維數(shù)值趨于2 時(shí)，表明區(qū)域化變量基本上是空間不相關(guān)的，若等于2，可以說整個(gè)樣本是空間獨(dú)立的[14]。從表2 可以看出，高寒草甸群落下土壤活性鐵的分形維數(shù)值較高，3個(gè)樣地均大于1.9，其中高寒草甸群落采樣尺度M1的分形維數(shù)值最大，達(dá)到2。隨著采樣尺度的增大，土壤活性鐵的分形維數(shù)值在高寒草甸群落下表現(xiàn)為減小趨勢(shì)， 而在高寒沼澤的群落下表現(xiàn)為逐漸增大的變化趨勢(shì)；在相同尺度下土壤活性鐵的分形維數(shù)值，高寒沼澤群落均小于高寒草甸群落，說明高寒沼澤群落土壤活性鐵比高寒草甸群落具有更強(qiáng)的空間異質(zhì)性。

2.4 不同尺度下土壤活性鐵的空間插值

利用ArcGIS 軟件中的地統(tǒng)計(jì)模塊對(duì)不同采樣尺度土壤活性鐵分布的Kriging 插值結(jié)果見圖1、圖2。高寒草甸樣地M1的土壤活性鐵呈現(xiàn)斑狀和帶狀分布，高值區(qū)主要分布在南部，低值區(qū)主要分布在北部，具有較好的空間連續(xù)性，緩和的梯度變化。樣地M2的土壤活性鐵分布呈現(xiàn)出條帶狀的梯度變化，均一性非常好，其中低值區(qū)與樣地內(nèi)的裸露斑塊的位置和延伸方向一致，高值區(qū)主要分布在南部；樣地 M3尺度下土壤活性鐵呈現(xiàn)出帶狀和斑塊狀鑲嵌的分布，且具有明顯的由北部向南部逐漸增大的趨勢(shì)，空間散布較均勻（圖1）。高寒沼澤樣地S1尺度下土壤活性鐵呈現(xiàn)比較分散的斑狀分布，說明其具有不太好的空間連續(xù)性以及較高的破碎化程度，在小尺度上形成強(qiáng)烈的變異；土壤活性鐵在樣地S2中呈現(xiàn)出弧狀和斑狀分布，分布總體呈現(xiàn)從西北向東南逐漸減少的趨勢(shì)，高值區(qū)主要分布在西北部，低值區(qū)主要分布在東南部；S3樣地的土壤活性鐵分布完全細(xì)化且梯度不明顯，反映較高的隨機(jī)性，但有明顯的聚鐵中心（圖2）。高寒沼澤群落土壤活性鐵空間分布空間連續(xù)性不好，隨機(jī)性較大，具有顯著的聚鐵中心，這也說明其變異函數(shù)各向異性比相對(duì)較低。

3討論

本研究區(qū)內(nèi)高寒草甸表層土壤活性鐵含量為7.04～779 g/kg，高寒沼澤土壤活性鐵含量為13.08～15.71 g/kg，2 種群落下土壤活性鐵含量差異明顯，高寒沼澤土壤活性鐵的平均含量高于高寒草甸，土壤活性鐵的含量和分布與研究區(qū)的群落結(jié)構(gòu)、植被蓋度、微地貌、土壤含水量、土壤pH值和土壤有機(jī)質(zhì)有一定關(guān)系。高寒沼澤群落具有土壤含水量高、冠層相對(duì)濕度大、溫度低等特征，抑制了土壤微生物的活動(dòng)，降低了有機(jī)質(zhì)分解的速度，土壤有機(jī)質(zhì)明顯高于高寒草甸[15]，有機(jī)質(zhì)與土壤活性鐵含量呈正相關(guān)關(guān)系[16]。同時(shí)，高寒沼澤群落的土壤含水量高，增強(qiáng)了土壤還原性，土壤氧化還原電位下降到一定水平后，鐵活性增強(qiáng)[17]。因此，高寒沼澤樣地土壤活性鐵含量高于高寒草甸。

本研究區(qū)內(nèi)高寒草甸土壤活性鐵屬于中等強(qiáng)度的空間變異性，其空間變異是由隨機(jī)性因素和結(jié)構(gòu)性因素共同作用引起的。但是高寒草甸群落土壤活性鐵在不同取樣尺度下均呈現(xiàn)一定的結(jié)構(gòu)性，5 m×2.5 m、10 m×5 m樣地總體來看塊基比較小，最大僅為35%，表現(xiàn)為強(qiáng)烈的空間自相關(guān)性或中等強(qiáng)度的空間自相關(guān)性，說明主要由地形、氣候等結(jié)構(gòu)性因素對(duì)高寒草甸土壤活性鐵空間變異起主導(dǎo)作用，而隨機(jī)性因素對(duì)其影響相對(duì)較小[18]。隨著采樣尺度的增大半方差函數(shù)的變程呈增大趨勢(shì)，而塊金值、基臺(tái)值和塊基比均隨采樣尺度增大先減小后增大，M1表現(xiàn)出各向同性，而M2、M3表現(xiàn)出各向異性，在M2中土壤活性鐵的方向變異與樣地內(nèi)裸露斑塊的出現(xiàn)，斑塊大小和分布方向基本一致。高寒草甸M1樣地的土壤活性鐵的分形維數(shù)為2，在空間自相關(guān)距離范圍內(nèi)同質(zhì)性較強(qiáng)，同質(zhì)性明顯[19]，趨于均質(zhì)，這與前面的半方差函數(shù)分析結(jié)果一致。隨樣點(diǎn)間距增大，樣地M2、M3表現(xiàn)出一定的空間自相關(guān)?？梢?，高寒草甸的樣地尺度不能設(shè)置太小，否則不能表現(xiàn)出土壤活性鐵的空間變異特征。另外，土壤水分和養(yǎng)分的運(yùn)移與地形變化密切相關(guān)，地勢(shì)低洼處，土壤水分和養(yǎng)分資源顯著聚集[20]。高寒草甸地勢(shì)有一定的起伏，屬于北高南低，土壤含水量南部區(qū)域可能較高，土壤還原性增強(qiáng)，故土壤活性鐵的含量高值都出現(xiàn)在研究區(qū)的南部，而且呈現(xiàn)由北向南增加的趨勢(shì)，印證了隨著采樣尺度的增大，隨機(jī)因素逐漸減弱，地形起伏等結(jié)構(gòu)因素逐漸加強(qiáng)。

高寒沼澤群落土壤活性鐵具有中等強(qiáng)度的空間變異性，表明由隨機(jī)性因素和結(jié)構(gòu)性因素共同引起的空間變異。但是高寒沼澤S2、S3樣地土壤活性鐵的塊基比值均較大，分別為50%、68%，表明在該樣地區(qū)域內(nèi)土壤活性鐵的空間變異主要是由微域內(nèi)變異和人為活動(dòng)等隨機(jī)性因素引起，而結(jié)構(gòu)性因素對(duì)其影響相對(duì)較小。在相同尺度下的土壤活性鐵的分形維數(shù)值，高寒沼澤均小于高寒草甸群落，說明高寒沼澤土壤活性鐵呈現(xiàn)出更強(qiáng)的空間異質(zhì)性。土壤水分含量較高的沼澤草甸，在某種程度上加劇了水分脅迫，養(yǎng)分資源的空間異質(zhì)性降低，促使物種競(jìng)爭逐漸增大，同時(shí)高原的溫度低以及水熱條件的不一致[21]，由于極端環(huán)境限制了其他種的存活和侵入，促使物種豐富度和均勻度較低，因而使群落多樣性較低[22]。高寒沼澤物種多樣性的降低通常會(huì)釋放出更加多樣化的生態(tài)位空間，存留物種會(huì)具有更廣的適應(yīng)性，也就是說物種多樣性低的群落種間性狀差異將更大[23]，其中植物根際環(huán)境以及元素吸收的顯著差異促使不同種植被下土壤活性鐵空間分布差異顯著[24]，土壤活性鐵的空間異質(zhì)性增強(qiáng)。另外，高寒沼澤樣地分布著大量的點(diǎn)狀草丘，面積不大，呈橢圓形或圓形，散布在沼澤水分多變化的地段，地表有薄層常年積水，在沼澤邊緣最多，是由烏拉草、木里薹草等形成的。水位從丘間洼地到丘上逐漸降低甚至消失，土壤含水量的降低促使通氣性逐漸增強(qiáng)，土壤由還原環(huán)境逐漸變?yōu)檠趸h(huán)境，丘上的土壤活性鐵含量低于丘間洼地，這是點(diǎn)狀草丘促使土壤活性鐵產(chǎn)生水平分異的一個(gè)主要原因[25]。所以，高寒沼澤微地貌的起伏對(duì)土壤活性鐵的空間分布產(chǎn)生了顯著的影響。另外，高寒沼澤的S1樣地的土壤活性鐵在小尺度空間上具有強(qiáng)烈的自相關(guān)性，受到了系統(tǒng)本身變異影響，主要由結(jié)構(gòu)性的因素引起，這是由沼澤特有的點(diǎn)狀草丘的微地貌表現(xiàn)出來的分布差異。但是半方差函數(shù)的塊金值、塊基比、基臺(tái)值均隨著采樣尺度的增大而增大，逐漸掩蓋小尺度的結(jié)構(gòu)特征，而變?yōu)橹饕呻S機(jī)性因素引起，雖然高寒沼澤處于圍欄禁牧期，但曾經(jīng)人為放牧以及放牧的強(qiáng)度會(huì)顯著影響土壤活性鐵含量的變化，管理粗放，其插值圖也相對(duì)零碎。

綜上所述，高寒沼澤群落土壤活性鐵含量明顯高于高寒草甸群落，同采樣尺度下土壤活性鐵變異系數(shù)變化范圍為14.2%～488%，均屬中等變異，且變異系數(shù)隨樣地尺度的增加而增大。本研究區(qū)內(nèi)的高寒沼澤和高寒草甸土壤活性鐵在不同采樣尺度下塊基比值均較大，介于25%和75%之間，屬于中等強(qiáng)度的空間變異性，其空間變異是由隨機(jī)性因素和結(jié)構(gòu)性因素共同作用引起的。另外，高寒沼澤群落土壤活性鐵比高寒草甸群落具有更強(qiáng)的空間異質(zhì)性。高寒草甸的空間變異，隨著采樣尺度的增大，地形起伏等結(jié)構(gòu)性因素的影響逐漸加強(qiáng)，而植被分布、人為活動(dòng)等隨機(jī)因素的影響逐漸減弱；但高寒沼澤正好相反。

參考文獻(xiàn)：

[1]袁可能.植物營養(yǎng)元素的土壤化學(xué)[M].北京：科學(xué)出版社，1983：344-351.

[2]Zhang M，Karathanasis A D. Characterization of iron-manganese concretions in Kentucky Alfisols with perched water tables[J]. Clays and Clay Minerals，1997，45（3）：428-439.

[3]文勇立，李輝，李學(xué)偉，等. 川西北草原土壤及冷暖季牧草微量元素含量比較[J]. 生態(tài)學(xué)報(bào)，2007，27（7）：2837-2846.

[4]張鐿鋰，李炳元，鄭度. 論青藏高原范圍與面積[J]. 地理研究，2002，21（1）：1-8.

[5]田玉強(qiáng)，歐陽華，徐興良，等. 青藏高原土壤有機(jī)碳儲(chǔ)量與密度分布[J]. 土壤學(xué)報(bào)，2008，45（5）：933-942.

[6]陳義強(qiáng)，劉國順，習(xí)紅昂. 微尺度下煙田鐵的空間變異性及其與煙葉鐵的相關(guān)分析[J]. 生態(tài)學(xué)報(bào)，2009，29（3）：1448-1458.

[7]馬媛，師慶東，楊建軍，等. 干旱區(qū)典型流域土壤微量元素的空間變異特征研究[J]. 干旱區(qū)地理，2006，29（5）：682-687.

[8]Goovaerts P. Geostatistical tools for characterizing the spatial variability of microbiological and physic-chemical soil properties[J]. Biology and Fertility of Soils，1998，27（4）：315-334.endprint

[9]Goovaerts P. Geostatistical modelling of uncertainty in soil science[J]. Geoderma，2001，103（1/2）：3-26.

[10]劉曉光，張少良，張興義，等. 典型黑土坡面速效鉀空間分布特征研究[J]. 土壤通報(bào)，2015，46（1）：133-141.

[11]張文敏，姜小三，吳明，等. 杭州灣南岸土壤有機(jī)碳空間異質(zhì)性研究[J]. 土壤學(xué)報(bào)，2014，51（5）：1087-1095.

[12]張忠華，胡剛. 南亞熱帶次生林土壤pH值與含水量的空間異質(zhì)性[J]. 生態(tài)科學(xué)，2014，33（1）：148-153.

[13]王正權(quán). 地統(tǒng)計(jì)學(xué)及其在生態(tài)學(xué)中的應(yīng)用[M].北京：科學(xué)出版社，1999：71-99.

[14]呂貽忠，李保國，崔燕. 不同植被群落下土壤有機(jī)質(zhì)和速效磷的小尺度空間變異[J]. 中國農(nóng)業(yè)科學(xué)，2006，39（8）：1581-1588.

[15]王長庭，龍瑞軍，丁路明. 青藏高原高寒嵩草草甸基本特征的研究[J]. 草業(yè)科學(xué)，2004，21（8）：16-19.

[16]陳家坊，何群，邵宗臣. 土壤中氧化鐵的活化過程的探討[J]. 土壤學(xué)報(bào)，1983，20（4）：387-393.

[17]唐羅忠，生原喜久雄，戶田浩人，等. 濕地林土壤的Fe2+、Eh及pH值的變化[J]. 生態(tài)學(xué)報(bào)，2005，25（1）：103-107.

[18]Cambardella C A，Moarmman T B，Novak J M，et al. Field scale variability of soil properties in Central lowa soils[J]. Soil Science Society of America Journal，1994，58（5）：1501-1511.

[19]劉丹，姜勇，梁文舉，等. 沈陽市郊耕地土壤交換性Ca含量的空間異質(zhì)性特征[J]. 生態(tài)學(xué)雜志，2003，22（4）：6-9.

[20]Li J R，Okin G S，Alvarez L，et al. Quantitative effects of vegetation cover on wind erosion and soil nutrient loss in a desert grassland of[LM]

southern New Mexico，USA[J]. Biogeochemistry，2007，85（3）：317-332.

[21]李英年，鮑新奎，曹廣民. 青藏高原正常有機(jī)土與草氈寒凍雛形土地溫觀測(cè)的比較研究[J]. 土壤學(xué)報(bào)，2001，38（2）：145-152.

[22]王長庭，龍瑞軍，王啟基，等. 高寒草甸不同草地群落物種多樣性與生產(chǎn)力關(guān)系研究[J]. 生態(tài)學(xué)雜志，2005，24（5）：483-487.

[23]Stubbs W J，Bastow W J. Evidence for limiting similarity in a sand dune community[J]. Journal of Ecology，2004，92（4）：557-567.

[24]吳沿友，郝建朝，李萍萍，等. 北固山濕地土壤金屬元素空間分布與變化[J]. 江蘇大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，26（4）：340-344.

[25]劉吉平，楊青，呂憲國. 三江平原典型環(huán)型濕地土壤營養(yǎng)元素的空間分異規(guī)律[J]. 水土保持學(xué)報(bào)，2005，19（2）：76-79.endprint