新疆伊犁地區(qū)草地土壤全碳含量空間格局分析

周李磊，朱華忠，鐘華平，楊 華

(1.中國科學院地理科學與資源研究所資源與環(huán)境信息系統(tǒng)國家重點實驗室，北京 100101；2.重慶師范大學地理與旅游學院三峽庫區(qū)地表過程與環(huán)境遙感重慶市重點實驗室，重慶 400047;3.中國科學院地理科學與資源研究所陸地表層格局與模擬院重點實驗室，北京 100101)

?

新疆伊犁地區(qū)草地土壤全碳含量空間格局分析

周李磊1,2，朱華忠1，鐘華平3，楊 華2

(1.中國科學院地理科學與資源研究所資源與環(huán)境信息系統(tǒng)國家重點實驗室，北京 100101；2.重慶師范大學地理與旅游學院三峽庫區(qū)地表過程與環(huán)境遙感重慶市重點實驗室，重慶 400047;3.中國科學院地理科學與資源研究所陸地表層格局與模擬院重點實驗室，北京 100101)

為了研究伊犁地區(qū)草地土壤碳庫的空間分布規(guī)律，本研究基于146個草地樣地的土壤全碳含量實測數據，結合遙感及氣象數據，并利用地理學相關分析方法，構建土壤全碳含量空間化模型，分析了伊犁地區(qū)草地不同土壤深度全碳含量空間格局特征，并對影響土壤全碳含量分布的主要因素進行探討。結果表明,1)伊犁地區(qū)草地不同土壤深度(0-10、10-20、20-30、30-40 cm和4個土層均值)的全碳含量與海拔、年均氣溫、≥10 ℃年積溫和8月份NDVI極顯著相關(P<0.01)，與年均降水和濕潤度顯著相關(P<0.05)，以這6個生態(tài)要素建立土壤全碳含量空間化模型，反演得到空間分辨率1 km的不同草地深度土壤全碳含量空間格局分布數據，檢驗結果表明空間化結果精度較高。2)伊犁地區(qū)草地0-40 cm土層土壤全碳平均含量為4.251 9%，并隨著土層深度的增加，草地土壤全碳含量呈下降趨勢。3)伊犁地區(qū)空間化結果反映不同深度土壤全碳含量空間格局存在差異，土壤全碳含量最小值出現在伊犁河沿岸兩側及特克斯河與伊犁河交匯處人類活動密集的區(qū)域，最大值均出現在南、北天山兩側，并沿著伊犁河向兩側高海拔區(qū)域逐漸增加。4)伊犁地區(qū)草地不同深度土壤全碳含量空間分布規(guī)律與草地類型分布規(guī)律一致；人為活動干擾(過度放牧、盲目開墾)對草地土壤全碳含量影響較大。

土壤全碳；空間格局；地統(tǒng)計學；草地；伊犁地區(qū)

隨著全球氣候變暖日益加劇，碳循環(huán)作為當前全球氣候變化研究的核心問題[1-2]，備受人們的關注，而草地生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的重要組成部分，在全球碳收支平衡中起重要作用[3-4]。中國天然草地面積4×108hm2，占國土面積的41.7%，主要分布在東北平原、內蒙古高原、黃土高原、青藏高原以及新疆山地[4]，是中國最大的陸地生態(tài)系統(tǒng)[5]。草地生態(tài)系統(tǒng)不具有明顯的地上碳庫，其90%碳庫儲存在土壤中[6]，其微弱變化均會對大氣CO2濃度產生巨大影響，進而影響全球碳平衡格局[7]。因此，分析草地土壤碳素含量的空間變化過程及其影響機理對了解全球碳循環(huán)有著重要意義。

目前，國內關于草地生態(tài)系統(tǒng)土壤碳庫(包括有機碳庫和無機碳庫)的研究主要基于不同草地類型[8-10]、不同利用方式[11-13]、人為干擾強度[14-15]以及氣候變化[16]對土壤碳庫的影響。不同植被類型對應的土壤碳含量差異明顯，對藏北高寒草地植被碳密度與碳貯量的研究[17]發(fā)現，不同草地間植被碳密度和碳貯量差異顯著；對黃土高原草地土壤有機碳含量研究[18]發(fā)現，有機碳含量最高的是高寒草甸，其空間變異最大，最小的是荒漠草原，其變異最小。隨著人口的增加，草畜矛盾突出，草地資源遭到破壞，人為干擾活動嚴重影響草地生態(tài)系統(tǒng)的碳儲量，在科爾沁草原上[19]研究發(fā)現，圍封草地土壤有機碳含量隨著土層的加深呈遞減的變化，放牧地土壤有機碳含量不同土層間無顯著差異；在錫林郭勒地區(qū)研究[20]發(fā)現，放牧強度在短期內對典型草原植被土壤有機碳含量沒有顯著影響；同時，溫度、降水等氣候因素的變化對草地土壤碳庫的積累產生重要影響，楊紅飛等[21]對氣候變化對草地系統(tǒng)土壤碳含量的影響進行了綜述，指出溫度和降水對草地植物土壤碳庫影響較大。伊犁地區(qū)位于新疆維吾爾自治區(qū)的西北部，是新疆最濕潤、降水最豐富、植被土壤發(fā)育最好的區(qū)域[22]，作為新疆的重點牧區(qū)，伊犁地區(qū)典型草地土壤無機碳含量隨著海拔的升高而降低，從土壤剖面垂直變化來看，除高寒草甸隨著土層深度的加深無明顯變化外，其它各類型草地無機碳含量均隨著土層加深而增大[23]；而短期放牧對昭蘇草甸草原土壤有機碳含量影響較小[24]。這些研究成果從氣候因子、植被因子和人為干擾因子等方面，對土壤碳含量進行了較為全面的研究，但都是基于經典統(tǒng)計方法，不能準確地描述土壤碳素積累的空間異質性；同時，新疆伊犁地區(qū)屬于半干旱區(qū)，土壤無機碳庫含量較高，單一的研究土壤有機碳或者土壤無機碳都不能很好地反映伊犁地區(qū)的碳庫分布情況。

基于此，以新疆伊犁地區(qū)草地為研究對象，利用地理學相關分析方法，構建草地土壤全碳含量空間化模型，分析土壤碳庫空間分布規(guī)律，以期為進一步開展伊犁地區(qū)草地土壤的碳循環(huán)研究奠定基礎。

1 材料與方法

1.1 研究區(qū)概況

伊犁地區(qū)(42°14′16″-44°53′30″ N，80°09′42″-84°56′50″ E)，位于歐亞大陸腹地，北與新疆境內的博爾塔拉蒙古自治州接壤，南與阿克蘇地區(qū)的拜城、溫宿縣隔山相望，東北與塔城地區(qū)的烏蘇縣相連，東南與巴音郭勒蒙古自治州的和靜縣相鄰，西與哈薩克共和國交界。區(qū)轄1市7縣和1個自治縣，分別為伊寧市、伊寧縣、霍城縣、尼勒克縣、鞏留縣、新源縣、特克斯縣、昭蘇縣和察布查爾錫伯自治縣，總面積55 848.23 km2。伊犁地區(qū)東、南、北三面高山環(huán)繞，地勢東高西低，東窄西寬，呈楔形向西敞開，來自西部里海的濕潤水汽和巴爾喀什湖的暖流沿伊犁河長驅直入谷地深處[25]，形成了溫帶大陸性半干旱氣候區(qū)。獨特的地形地貌使得伊犁地區(qū)降水充沛，境內有伊犁河、喀什河和特克斯河3條主要河流，年均降水量為200～800 mm，年均氣溫為2.9～9.1 ℃，年均日照時數為2 700～3 000 h[26]，是新疆最濕潤的地區(qū)，素有“中亞濕島”之稱。新疆伊犁地區(qū)自然條件優(yōu)越，農、牧業(yè)發(fā)展優(yōu)勢明顯，農畜產品豐富，是新疆的糧倉和全國著名的牧場。伊犁地區(qū)天然草地類型主要有溫性草甸草原類、溫性草原類、溫性荒漠草原類、溫性荒漠類、低地草甸類、山地草甸類和高寒草甸類，天然草地面積約為3.71×106hm2，占新疆草地資源的6.38%，其產草量以及載畜能力均處于全疆最好水平[27-28]。

1.2 數據來源及處理

1.2.1 樣品采集與測定 綜合考慮伊犁地區(qū)草地類型的分布特點、草地利用方式和利用強度等，在全區(qū)范圍內設置了146個草地樣地(圖1)，于2013年7月-9月對全部樣地進行實地GPS定位，同步記錄經緯度坐標、海拔和利用方式，依據代表性原則，選擇代表整個樣地植被、地形及土壤等特征的地段，設置3個樣方，在完成草地樣方調查的基礎上，清理樣方土壤表面的植物殘留物和雜質，用土鉆在0-40 cm土壤內每隔10 cm取土，即0-10、10-20、20-30、30-40 cm，帶回實驗室自然晾干后，去除植物根系及礫石后，用研缽磨碎過149 μm篩。全碳(TC)含量由中科院植物研究所采用干燒法，運用元素分析儀vario MACRO cube測得。

圖1 伊犁地區(qū)位置及樣點分布Fig.1 The location of Ili area and distribution of sample sites

1.2.2 遙感數據 研究區(qū)內中分辨率成像光譜儀(moderate-resolution imaging spectroradiometer，MODIS)產品數據，主要是指歸一化植被指數(normalized difference vegetation index，NDVI)，數據來源于美國地質調查局(http://ladsweb.nascom.nasa.gov)的成品數據產品(MOD 13Q1)，數據采集的時間為2013年8月11日，用MODIS Reprojection Tool(MRT)拼接轉投影，空間分辨率1 km。

1.2.3 氣候數據 包括1957-2012年伊犁地區(qū)的年均氣溫、年均降雨量、≥10 ℃年積溫、濕潤度(伊萬諾夫濕潤度[28])，根據氣象站點實測數據，運用ANUSPLIN專業(yè)氣象插值軟件，空間插值得到，空間分辨率1 km。

1.2.4 其它數據 數字高程模型(DEM)從國家地球系統(tǒng)科學數據共享平臺(www.geodata.cn)下載，空間分辨率30 m；伊犁地區(qū)草地類型數據為20世紀80年代全國草地調查1∶100萬草地類型圖矢量數據。

1.3 分析方法

1.3.1 生態(tài)要素的選取 從影響土壤碳積累的氣候、地形和植被因素考慮，選取年均氣溫、年均降水量、≥10 ℃年積溫、濕潤度、海拔和NDVI這6個生態(tài)要素，來評價草地土壤全碳的空間分布情況。

1.3.2 數據分析 將土壤全碳含量數據按0-10、10-20、20-30、30-40 cm分層統(tǒng)計，0-40 cm取0-10、10-20、20-30和30-40 cm 4層土壤全碳含量的均值。按照樣地編號，將能被3整除的取出，每層共48個樣地數據作為檢驗數據，剩下的98個樣地數據作為分析統(tǒng)計基礎數據；刪除明顯異常點，最終每層得到41個檢驗數據和93個分析數據，每層共134個有效數據；利用ArcGIS 10.3軟件對樣點數據及生態(tài)要素數據做數據提取及空間化表達，用Origin 9.1進行土壤全碳與各生態(tài)要素間的相關回歸分析及制圖，建立回歸方程，探討不同土壤深度全碳含量與各生態(tài)要素間的相互關系。

1.3.3 土壤全碳含量空間化模型 確定各生態(tài)要素對土壤全碳的權重，構建了多元數據的土壤全碳含量空間化模型：

(1)

(2)

(3)

式中，P表示土壤全碳，wi表示某種生態(tài)要素對土壤全碳的影響權重，yi表示某種生態(tài)要素的與土壤全碳的回歸方程，即生態(tài)要素與實測土壤全碳回歸所得到的土壤全碳預測值，i=1、2……6，表示6個生態(tài)要素；Re表示比值；R2表示反演結果與實測結果的復相關系數，R2越大，表示反演結果與實測結果擬合度越高；K表示離散程度，即實際測定值和反演值回歸模型的斜率，K越接近1說明反演結果越和實測結果一致。

1.3.4 空間化結果精度評價方法 空間化結果精度通過41組檢驗數據與空間化模型計算結果進一步計算均方根誤差(rootmeansquareerror，RMSE)、偏離度(E)和相關系數(R)來檢驗[29]。

(4)

(5)

(6)

2 結果與分析

2.1 草地土壤全碳含量與生態(tài)要素回歸分析

分別建立草地不同土壤深度全碳含量與海拔、年均氣溫、年均降水量、≥10 ℃年積溫、濕潤度和NDVI的回歸方程。由于篇幅所限，以4個土層土壤全碳含量均值(0-40cm)為例，回歸方程如圖2所示。其中，土壤全碳含量與海拔、年均氣溫、≥10 ℃年積溫和濕潤度呈拋物線狀分布趨勢，與年均降水和NDVI呈逐步遞增趨勢。

回歸分析結果顯示(表1)，草地土壤全碳含量(0-40cm)與海拔、年均氣溫、≥10 ℃年積溫、8月份NDVI呈極顯著相關關系(P<0.01)，草地土壤全碳含量(0-40cm)與年均降水、濕潤度呈顯著相關關系(P<0.05)；草地土壤全碳含量(0-40cm)與NDVI和海拔的協相關系數最大，R2分別為0.446和0.461，草地土壤全碳含量(0-40cm)與年均降水的協相關系數最小，R2為0.226，但草地土壤全碳含量(0-40cm)與生態(tài)要素的回歸分析全部滿足99%可信度的F檢驗，說明草地土壤全碳含量(0-40cm)與生態(tài)要素的相關關系成立且非常顯著，線性回歸擬合的方程是有意義的。

2.2 生態(tài)評價因子權重的確定

根據0-10、10-20、20-30、30-40和0-40cm土壤全碳含量與6個生態(tài)要素的回歸方程，分別進行單要素草地全碳含量空間反演，提取參與分析的93個樣點所對應的反演結果，與樣點實測數據進行線性擬合，計算擬合方程斜率k和復相關系數R2，按照公式(2)和公式(3)計算各生態(tài)要素的權重(表2)。

圖2 草地土壤全碳含量均值與生態(tài)要素的回歸分析Fig.2 Regression analysis of the soil total carbon with different ecological factors表1 土壤全碳含量均值與生態(tài)要素的回歸分析結果Table 1 Test of regression analysis between soil total content and ecological factors

生態(tài)要素Ecologicalfactor PR2FSig.海拔Altitude0.303**0.46138.4900.000年均氣溫Averageannualtemperature-0.271**0.39729.4500.000年均降水量Averageannualrainfall0.236*0.22612.0040.000≥10℃年積溫≥10℃Annualaccumulatedtemperature-0.382**0.41531.8870.000濕潤度Humidity0.227*0.36423.7720.0008月份NDVI指數NDVIindexinAugust0.630**0.44673.2030.000

注：**表示在0.01水平達到極顯著相關，*表示在0.05水平達到顯著相關。

Note：* and ** indicate significant correlation at 0.05 and 0.01 level, respectively.

表2 不同草地土壤深度生態(tài)要素權重Table 2 Ecological factors weight at different soil depth

2.3 草地土壤全碳含量空間化

按照土壤全碳含量空間化模型，在ArcGIS中加權反演伊犁地區(qū)草地土壤全碳含量空間格局分布，用伊犁地區(qū)草地類型矢量數據掩膜裁剪出草地區(qū)域，排除非草地區(qū)域(包括伊犁河兩側的大片耕地、天山和博羅科努山山脊的裸巖以及居住用地等非草地區(qū)域)對空間化結果的影響，最終得到伊犁地區(qū)草地不同土壤深度全碳含量的空間格局分布(圖3)。

2.4 草地土壤全碳含量評價結果檢驗

在0-10、10-20、20-30、30-40 cm這4個土壤層之間，30-40 cm土層的均方根誤差相對較小，為0.75%，壤層的偏離度最大，為26.66%，平均預測精度達到73.34%；0-10，10-20 cm土層的復相關系數全部大于0.60，30-40 cm土層復相關系數最小，只有0.42(表3)。整體上(0-40 cm)，空間化結果與實測數據間的均方根誤差為0.96%，整體偏離度為23.33%，平均預測精度達到76.67%，復相關系數0.59，整體上空間化結果精度較高。

表3 精度檢驗Table 3 Testing of the accuracy

2.5 草地土壤全碳含量空間格局分析

關于這四首鋼琴小曲的音樂分析，我們曾在《中國鋼琴音樂的拓荒者—新發(fā)現的早期鋼琴家李樹化的照片史料》一文③中有簡單的藝術分析，可供讀者參考。

伊犁地區(qū)全碳含量較低區(qū)域出現在霍城縣南部、伊寧市南部，察布查爾錫伯自治縣北部、鞏留縣北部以及新源縣西部的低海拔區(qū)域，主要分布著溫性荒漠類草地，草地質量較差；而與這些草地相鄰的伊犁河沿岸兩側區(qū)域水土肥沃，分布著大片耕地，是伊犁地區(qū)主要的糧食產地，也是主要城市的聚集地，人類活動對草地的干擾較大；全碳含量較高區(qū)域，主要分布在南北天山兩側及中部高山區(qū)，該區(qū)域主要是山地草甸類草地，草地質量較好，人類活動干擾較少。

其中，0-10、10-20、20-30和30-40 cm土層的全碳含量空間分布格局顯示。不同土壤深度草地土壤全碳含量空間格局分布差異明顯，土壤全碳含量最大值都出現在南、北天山兩側，最小值出現在伊犁河沿岸兩側及特克斯河與伊犁河交匯處人類活動密集的區(qū)域，土壤全碳含量沿著河流向兩側高山逐漸增加；且隨著土壤層深度的增加，土壤全碳含量呈現出減少的趨勢，草地土壤全碳含量大于5%的區(qū)域在0-10 cm土層全碳空間分布最廣，10-20 cm土層相對減少，20-30 cm土層近一步減少，直到30-40 cm土層不存在土壤全碳含量大于5%的區(qū)域。

伊犁地區(qū)0-40 cm草地土壤全碳含量平均水平為4.251 9%，其中，0-10 cm土層草地土壤全碳含量平均水平為5.101 6%，10-20 cm土層草地土壤全碳含量平均水平為4.538 7%，20-30 cm土層草地土壤全碳含量平均水平為4.069 3%，30-40 cm土層草地土壤全碳含量平均水平為3.535 1%(表4)。隨著土層深度的增加，草地土壤全碳含量呈下降趨勢，與空間格局分布規(guī)律一致。

不同深度不同草地類型的土壤全碳含量不同，0-40 cm土層草地土壤全碳含量最高的是溫性山地草甸,為5.188 4%，其次是溫性草甸草原,為4.882 3%，最小的是溫性荒漠類,為1.899 6%。0-10、10-20、20-30和30-40 cm土層，土壤全碳含量最高的是溫性山地草甸，土壤全碳含量分別為6.279 6%、5.547 9%、4.868 7%和4.207 4%；其次是溫性草甸草原，分別為5.835 2%、5.210 4%、4.614 1%和4.032 5%；最小的是溫性荒漠類，分別為2.155 5%、1.950 9%、1.970 4%和1.917 7%。但不同草地類型土壤全碳含量在各深度土壤層里表現出相同的規(guī)律，即溫性山地草甸>溫性草甸草原>高寒草甸>改良草地>溫性草原>溫性荒漠草原>溫性荒漠。其中，溫性山地草甸、溫性草甸草原、高寒草甸和改良草地這4類草地類型的土壤全碳含量在0-10、10-20和20-30 cm土層中，大于各層平均值；溫性草原、溫性荒漠草原和溫性荒漠這3類小于各層平均值；在30-40 cm土層中，溫性山地草甸、溫性草甸草原、改良草地和溫性草原這4類草地類型的土壤全碳含量大于全區(qū)平均值；高寒草甸、溫性荒漠草原和溫性荒漠這3類小于全區(qū)平均值。

圖3 不同土壤深度草地土壤全碳含量空間分布Fig.3 Spatial distribution of soil total carbon content of different soil depth表4 不同草地類型土壤全碳含量(%)統(tǒng)計結果Table 4 Soil total carbon content(%)of different depths of different type grassland

草地類型Grasslandtype 0-10cm10-20cm20-30cm30-40cm0-40cm溫性草甸草原類Temperatemeadow-steppetype5.83525.21044.61414.03254.8823溫性草原類Temperatesteppetype4.79124.37833.97893.62744.1349溫性荒漠草原類Temperatedesert-steppetype4.40154.04693.72513.41883.8268溫性荒漠類Temperatedeserttype2.15551.95091.97041.91771.8996溫性山地草甸類Temperatemontanemeadowtype6.27965.54794.86874.20745.1884高寒草甸類Alpinemeadowtype5.52114.90354.39483.39054.4681改良草地Improvedgrassland5.23734.63434.13603.76634.4062均值Mean5.10164.53874.06933.53514.2519

3 討論

3.1 生態(tài)評價因子與草地土壤全碳含量的關系

土壤性質的空間異質性是地形、氣候、植被、土壤母質以及外界干擾等因素共同作用的結果[29-32]。伊犁地區(qū)草地土壤全碳含量與海拔、年均氣溫、≥10 ℃年積溫和8月份NDVI呈極顯著相關關系(P<0.01),與年均降水和濕潤度呈顯著相關關系(P<0.05);隨著海拔、年均溫度、≥10 ℃年積溫和濕潤度的增加，土壤全碳含量也呈現出先增加后減少的趨勢；隨著降水和NDVI的增加，土壤全碳含量呈現出正向增加的趨勢。這是由于隨著海拔的增高，氣溫降低，主要是溫性山地草甸類草地，草地植被覆蓋度和生物量高，地上生物量及地下根系生物量充足，土壤有機碳分解比較慢，根系和微生物活動弱[23]，沉淀形成的碳酸鈣少，積累的有機質多[18,30]，進而導致土壤全碳含量增加；但海拔超過一定高度，氣溫進一步降低，主要是高寒草甸類草地，草地植被長勢較差，地上生物量及地下根系生物量作為土壤全碳含量的來源相比其它區(qū)域減少[13]，導致土壤全碳含量降低；相對應的低海拔河谷地區(qū)氣溫高，主要是溫性荒漠類草地，草地質量較差，土壤有機碳分解較快，不利于土壤有機碳的積累，區(qū)域全碳含量以無機碳為主，相比其它區(qū)域全碳含量減少。

3.2 不同草地類型對土壤全碳含量空間格局分布的影響

楊玉海等[22]對伊犁地區(qū)土壤有機碳含量的研究表明，高山草甸>草甸草原>荒漠草原>沼澤類，韓文娟等[23]對伊犁地區(qū)土壤無機碳含量的研究表明，溫性荒漠>溫性荒漠草原>溫性草原>溫性草甸草原>溫性山地草甸>高寒草甸。本研究表明，伊犁地區(qū)不同草地土壤深度(0-10、10-20、20-30、30-40 cm和4個土層均值)全碳含量表現為溫性山地草甸>溫性草甸草原>高寒草甸>改良草地>溫性草原>溫性荒漠草原>溫性荒漠。這種不同的結果可能與溫度的高低，植被凋落物的多少及微生物活動[35]的強弱有關。

從土壤全碳含量空間分布(圖3)及分區(qū)統(tǒng)計表(表4)可知，伊犁地區(qū)溫性荒漠分布區(qū)土壤全碳含量低，而溫性山地草甸分布區(qū)土壤全碳含量高，不同草地類型所對應的土壤全碳含量差異較大，草地土壤全碳含量空間分布格局符合草地類型分布規(guī)律。

3.3 人類活動干擾對土壤全碳含量空間格局分布的影響

隨著人口的增加，社會經濟的發(fā)展，過度放牧[24]、草地開墾[36]等人為活動已嚴重干擾了天然草地資源，草地退化、水土流失[37]和荒漠化問題[38-39]逐漸突出。伊犁河河谷兩側，包括霍城縣南部、察布查爾錫伯自治縣北部、鞏留縣北部、新源縣西部、伊寧市北部和伊寧市南部，土壤全碳含量較其它區(qū)域低，其原因可能是因為伊犁河兩側居民點較多，耕地分布比較集中，受傳統(tǒng)觀念及人口增長過快等因素的影響，天然草地被大量開墾，據統(tǒng)計1950-2002年伊犁河谷已累積開墾29.18萬hm2[40]，人為盲目地就近開墾及過度的放牧破壞了草地資源[41]，使得草地退化，水土流失加劇，土壤有機碳大量流失造成土壤碳素減少。而高海拔區(qū)域土壤及氣候不適合農作物的種植，且居民定居點較為松散及稀少，人類活動對草地干擾強度小，故土壤全碳含量較河谷平原地區(qū)高。

3.4 草地土壤全碳含量空間化方法

草地土壤全碳含量的實測數據往往是分散的點數據，實現草地土壤全碳含量的空間化表達，采用最多的是克里格插值法，但是空間插值的結果，在插值點數據較少或者分布不均勻的情況下，往往存在局部插值結果偏大或者偏小的情況。本研究通過分析土壤全碳含量與生態(tài)因子間的相關性，構建了土壤全碳含量空間化模型，實現了土壤全碳含量的空間化表達。與直接用克里格插值法對實測土壤有機碳含量做空間插值[42]的結果相比較，不存在局部偏大偏小的異常值；有研究表明[43]，基于土地類型修正的克里格插值法的有機碳空間化插值的結果優(yōu)于普通克里格插值法、反距離權重法和徑向基函數法(RBF)，而本研究同時考慮了6個生態(tài)因子，空間化結果更優(yōu)；與傳統(tǒng)插值方法相比，空間分布更加連續(xù)自然，不存在“牛眼”現象。本研究采用單因子回歸分析而沒有采用多因子回歸，是因為土壤全碳含量與生態(tài)因子間存在的不是簡單線性關系，多因子回歸是線性擬合，會降低空間化結果的精度。

由于缺乏0―40 cm土壤全碳含量的實測數據，本研究以0-10、10-20、20-30和30-40 cm這4層土壤全碳含量的均值代替，不同土層(0-10、10-20、20-30、30-40和0-40 cm)全碳含量與生態(tài)因子的回歸分析趨勢相同，因此，本研究沒有將所有數據分析全部列出來，而是以0-40 cm土層為代表進行了分析討論，結果較好地反映了伊犁地區(qū)草地土壤全碳含量的空間格局分布規(guī)律，但如何賦予不同土層(0-10、10-20、20-30和30-40 cm)不同的權重來求取0-40 cm土層土壤全碳的含量，提高數據的精度，在今后的研究中有待進一步探討。

4 結論

本研究利用2013年伊犁地區(qū)146個樣地土壤全碳含量實測數據，結合1957-2012年伊犁地區(qū)的氣候數據(年均氣溫、年均降水量、≥10 ℃年積溫和濕潤度氣候數據)、DEM數據以及2013年8月份的NDVI數據，采用地理學相關分析方法，構建了伊犁地區(qū)草地土壤全碳含量空間化模型，分析伊犁地區(qū)草地土壤全碳含量空間分布特征，得到以下結論：

1)伊犁地區(qū)不同草地土壤深度(0-10、10-20、20-30、30-40 cm和4個土層均值)全碳含量與海拔、年均氣溫、≥10 ℃年積溫、8月份NDVI呈極顯著相關關系(P<0.01)，與年均降水量和濕潤度呈顯著相關關系(P<0.05)。

2)伊犁地區(qū)草地土壤全碳含量空間化結果，直觀反映了伊犁地區(qū)草地土壤全碳含量空間分布格局。土壤全碳含量最小值出現在伊犁河沿岸兩側及特克斯河與伊犁河交匯處人類活動密集的區(qū)域，最大值都出現在南、北天山兩側，沿著伊犁河向兩側高海拔區(qū)域逐漸增加。

3)不同土壤深度，其草地土壤全碳含量差異顯著。伊犁地區(qū)0-40 cm土層(4個土層均值)草地土壤全碳含量為4.251 9%，其中，0-10、10-20、20-30和30-40 cm土層草地土壤全碳含量平均水平分別為5.101 6%、4.538 7%、4.069 3%和3.53 1%，隨著土層深度的增加，草地土壤全碳含量呈下降趨勢。

4)伊犁地區(qū)草地土壤全碳含量空間格局分布受地形和植被影響較大，不同深度草地土壤全碳含量空間分布規(guī)律與草地類型分布一致；人為活動干擾(過度放牧、盲目開墾)對草地土壤全碳含量影響較大。

致謝：感謝課題組鐘橫、張祎、劉清和周星杰等同學對野外采樣的幫助，特此感謝！

References：

[1] 陳曉鵬,尚占環(huán).中國草地生態(tài)系統(tǒng)碳循環(huán)研究進展.中國草地學報,2011,33(4):99-110.

Chen X P,Shang Z H.Progress of carbon cycle research in China grassland ecosystem.Chinese Journal of Grassland,2011,33(4):99-110.(in Chinese)

[2] 劉淑麗,林麗,郭小偉,李婧,歐陽經政,杜巖功,張法偉,李以康,曹廣民.青海省高寒草地土壤無機碳儲量空間分異特征.生態(tài)學報,2014,34(20):5953-5961.

Liu S L,Lin L,Guo X W,Li J,Ouyang J Z,Du Y G,Zhang F W,Li Y K,Cao G M.The variation feature of soil inorganic carbon storage in alpine grassland in Qinghai Province.Acta Ecologica Sinica,2014,34(20):5953-5961.(in Chinese)

[3] Scurlock J M O,Johnson K,Olson R J.Estimating net primary productivity from grassland biomass dynamics measurements.Global Change Biology,2002,8(8):736-753.

[4] 方精云,楊元合,馬文紅,安尼瓦爾·買買提,沈?；?中國草地生態(tài)系統(tǒng)碳庫及其變化.中國科學:生命科學,2010,40(7):566-576.

Fang J Y,Yang Y H,Ma W H,Mohammat A,Shen H H.Ecosystem carbon stocks and their changes in China’s grasslands.Science China：Life Sciences,2010,40(7):566-576.(in Chinese)

[5] 李學斌,樊瑞霞,劉學東.中國草地生態(tài)系統(tǒng)碳儲量及碳過程研究進展.生態(tài)環(huán)境學報,2014,23(11):1845-1851.

Li X B,Fan R X,Liu X D.Advance in studies on carbon storage and carbon process in grassland ecosystem of China.Ecology and Environmental Sciences,2014,23(11):1845-1851.(in Chinese)

[6] 肖勝生,董云社,齊玉春,彭琴,何亞婷,楊智杰.草地生態(tài)系統(tǒng)土壤有機碳庫對人為干擾和全球變化的響應研究進展.地球科學進展,2009,24(10):1138-1148.

Xiao S S,Dong Y S,Qi Y C,Peng Q,He Y T,Yang Z J.Advance in responses of soil organic carbon pool of grassland ecosystem to human effects and global changes.Advances in Earth Science,2009,24(10):1138-1148.(in Chinese)

[7] Wei J,Wu G,Deng H B,Zhao J Z.Spatial pattern of soil carbon and nutrient storage at the alpine tundra ecosystem of Changbai Mountain,China.Journal of Forestry,2004,15(4):249-254.

[8] 程積民,程杰,楊曉梅,劉偉，陳芙蓉.黃土高原草地植被碳密度的空間分布特征.生態(tài)學報,2012,32(1):226-237.

Cheng J M,Cheng J,Yang X M,Liu W,Chen F R.Spatial distribution of carbon density in grassland vegetation of the Loess Plateau of China.Acta Ecologica Sinica,2012,32(1):226-237.(in Chinese)

[9] 華娟,趙世偉,張揚,馬帥.云霧山不同草地群落土壤活性有機碳分布特征.草地學報,2009,17(3):315-320.

Hua J,Zhao S W,Zhang Y,Ma S.Distribution characteristics of soil labile organic carbon of different grassland communities in Yunwu Mountain.Acta Agrestia Sinica,2009,17(3):315-320.(in Chinese)

[10] 任繼周,林慧龍.草地土壤有機碳儲量模擬技術研究.草業(yè)學報,2013，22(6):280-294.

Ren J Z,Lin H L.Study on the simulation methods of grassland soil organic carbon:A review.Acta Prataculturae Sinica,2013，22(6):280-294.(in Chinese)

[11] 代景忠,衛(wèi)智軍,何念鵬,王若夢,溫學華,張云海,趙小寧,于貴瑞.封育對羊草草地土壤碳礦化激發(fā)效應和溫度敏感性的影響.植物生態(tài)學報,2012,36(12):1226-1236.

Dai J Z,Wei Z J,He N P，Wang R M，Wen X H，Zhang Y H，Zhao X N，Yu G R.Effect of grazing enclosure on the priming effect and temperature sensitivity of soil C mineralization inLeymuschinensisgrasslands,Inner Mongolia,China.Chinese Journal of Plant Ecology，2012，36(12):1226-1236.(in Chinese)

[12] 何念鵬,韓興國,于貴瑞.長期封育對不同類型草地碳貯量及其固持速率的影響.生態(tài)學報,2011,31(15):4270-4276.

He N P,Han X G,Yu G R.Carbon and nitrogen sequestration rate in long-term fenced grasslands in Inner Mongolia,China.Acta Ecologica Sinica，2011，31(15):4270-4276.(in Chinese)

[13] 徐沙,龔吉蕊,張梓榆,劉敏,王憶慧,羅親普.不同利用方式下草地優(yōu)勢植物的生態(tài)化學計量特征.草業(yè)學報,2014,28(6):45-53.

Xu S,Gong J R,Zhang Z Y,Liu M,Wang Y H,Luo Q P.The ecological stoichiometry of dominant species in different land uses type of grassland.Acta Prataculturae Sinica,2014,28(6):45-53.(in Chinese)

[14] 任海彥,鄭淑霞,白永飛.放牧對內蒙古錫林河流域草地群落植物莖葉生物量資源分配的影響.植物生態(tài)學報,2009,33(6):1065-1074.

Ren H Y,Zheng S X,Bai Y F.Effects of grazing on foliage biomass allocation of grassland communities in Xilin river basin, Inner Mongolia.Chinese Journal of Plant Ecology,2009,33(6):1065-1074.(in Chinese)

[15] 張凡,祁彪,溫飛,張德罡,吳紅,張力.不同利用程度高寒干旱草地碳儲量的變化特征分析.草業(yè)學報,2011,20(4):11-18.

Zhang F,Qi B,Wen F,Zhang D G,Wu H,Zhang L.Analysis of the change of carbon storage in alpine arid grassland.Acta Prataculturae Sinica,2011,20(4):11-18.(in Chinese)

[16] 范月君,侯向陽,石紅霄,師尚禮.氣候變暖對草地生態(tài)系統(tǒng)碳循環(huán)的影響.草業(yè)學報,2012,21(3):294-302.

Fan Y J,Hou X Y,Shi H X,Shi S L.Effect of carbon cycling in grassland ecosystems on climate warming.Acta Prataculturae Sinica,2012,21(3):294-302.(in Chinese)

[17] 常天軍,王建林,李鵬,成海宏,方華麗.藏北高寒草地植被的碳密度與碳貯量.生態(tài)科學,2007,26(5):437-442.

Chang T J,Wang J L,Li P,Cheng H H,Fang H L.Carbon density and reserve of alpine grassland vegetation in northern Tibet.Ecological Science,2007,26(5):437-442.(in Chinese)

[18] 劉偉,程積民,高陽,程杰,梁萬鵬.黃土高原草地土壤有機碳分布及其影響因素.土壤學報,2012,49(1):68-76.

Liu W,Cheng J M,Gao Y,Cheng J,Liang W P.Distribution of soil organic carbon in grassland on loess plateau and its influencing factors.Acta Pedologica Sinica，2012,49(1):68-76.(in Chinese)

[19] 孔玉華,姚風軍,鵬爽,劉艷,董文軒,白龍.不同利用方式下草地土壤碳積累及匯/源功能轉換特征研究.草業(yè)科學,2010,27(4):40-45.

Kong Y H,Yao F J,Peng S,Liu Y,Dong W X,Bai L.Study on the characteristics of soil carbon accumulation and conversion of carbon sink and source of grassland under different land use types.Pratacultural Science,2010,27(4):40-45.(in Chinese)

[20] 王艷芬,陳佐忠,Tieszen L T.人類活動對錫林郭勒地區(qū)主要草原土壤有機碳分布的影響.植物生態(tài)學報,1998,22(6):66-72.

Wang Y F,Chen Z Z,Tieszen L T.Distribution of soil organic carbon in the major grasslands of Xilingela,Inner Mongolia,China.Acta Phytoecologica Sinica，1998,22(6):66-72.(in Chinese)

[21] 楊紅飛,穆少杰,李建龍.氣候變化對草地生態(tài)系統(tǒng)土壤有機碳儲量的影響.草業(yè)科學,2012,29(3):392-400.

Yang H F,Mu S J,Li J L.Effects of climate change on soil organic carbon storage of grassland ecosystem.Pratacultural Science,2012,29(3):392-400.(in Chinese)

[22] 楊玉海,陳亞寧,李衛(wèi)紅,陳亞鵬.伊犁河谷不同植被帶下土壤有機碳分布.地理學報,2010,65(5):605-612.

Yang Y M,Chen Y N,Li W H,Chen Y P.Soil organic carbon distribution of different vegetation types in the Ili river valley.Acta Geographica Sinica,2010,65(5):605-612.(in Chinese)

[23] 韓文娟,盛建東,谷海斌,劉利利.新疆伊犁地區(qū)典型草地土壤無機碳分布特征.新疆農業(yè)大學學報,2015,38(2):140-145.

Han W J,Sheng J D,Gu H B,Liu L L.Soil inorganic carbon distribution characteristics in typical grassland in Ili area of Xinjiang.Journal of Xinjiang Agricultural University,2015,38(2):140-145.(in Chinese)

[24] 蒲寧寧,孫宗玖,范燕敏,楊合龍.放牧強度對昭蘇草甸草原土壤有機碳及微生物碳的影響.新疆農業(yè)大學學報,2013,36(1):66-70.

Pu N N,Sun Z J,Fan Y M,Yang H L.Influence of grazing intensity on the soil organic carbon and microbial biomass carbon of meadow steppe in Zhaosu area.Journal of Xinjiang Agricultural University,2013,36(1):66-70.(in Chinese)

[25] 章祖同,劉起.中國重點牧區(qū)草地資源及其開發(fā)利用.北京:中國科學技術出版社,1992.

[26] 周聿超.新疆河流水文水資源.烏魯木齊:新疆科技衛(wèi)生出版社,1999.

[27] 劉芳,張紅旗,董光龍.伊犁河谷草地植被NDVI變化及其降水敏感性特征.資源科學,2014,36(8):1724-1731.

Liu F,Zhang H Q,Dong G L.Vegetation dynamics and precipitation sensitivity in Yili valley grassland.Resources Science,2014,36(8):1724-1731.(in Chinese)

[28] 孫慧蘭,陳亞寧,李衛(wèi)紅,黎楓,木巴熱克·阿尤普.新疆伊犁河流域草地類型特征及其生態(tài)服務價值研究.中國沙漠,2011,31(5):1273-1277.

Sun H L,Chen Y N,Li W H,Li F,Mubareke·Ayousi.Study on types and ecological services values of the grassland in the Ili river basin,Xinjiang,China.Journal of Desert Research,2011,31(5):1273-1277.(in Chinese)

[29] Xu S,Wu C,Wang L,Gonsamo A,Shen Y,Niu Z.A new satellite-based monthly precipitation downscaling algorithm with non-stationary relationship between precipitation and land surface characteristics.Remote Sensing of Environment,2015,162:119-140.

[30] 湯愛坤,劉汝海,許廖奇,王金玉,劉一霆.昌邑海洋生態(tài)特別保護區(qū)土壤養(yǎng)分的空間異質性與植物群落的分布.水土保持通報,2011,31(3):88-93.

Tang A K,Liu R H,Xu L Q,Wang J Y,Liu Y T.Spatial heterogeneity of soil nutrients and distribution of plant community in Changyi marine ecological special protection area.Bulletin of Soil and Water Conservation,2011,31(3):88-93.(in Chinese)

[31] 岳躍民,王克林,張偉,陳洪松,王敏.基于典范對應分析的喀斯特峰叢洼地土壤-環(huán)境關系研究.環(huán)境科學,2008,29(5):1400-1405.

Yue Y M,Wang K L,Zhang W,Chen H S,Wang M.Relationships between soil and environment in Peak-Cluster depression areas of Karst region based on canonical correspondence analysis.Environmental Science,2008,29(5):1400-1405.(in Chinese)

[32] 烏云娜,雒文濤,霍光偉,李海山,胡高娃.草原群落退化演替過程中微斑塊土壤碳氮的空間異質動態(tài).生態(tài)學報,2014:34(19):5549-5557.

Wuyunna,Luo W T,Huo G W,Li H S,Hu G W.Spatial heterogeneity dynamics of soil carbon and nitrogen in the mini-patches during degeneration succession of grassland communities.Acta Ecologica Sinica,2014,34(19):5549-5557.(in Chinese)

[33] 羅上華,毛齊正,馬克明,鄔建國.北京城市綠地表層土壤碳氮分布特征.生態(tài)學報,2014,34(20):6011-6019.

Luo S H,Mao Q Z,Ma K M,Wu J G.Spatial distribution of soil carbon and nitrogen in urban greenspace of Beijing.Acta Ecologica Sinica,2014,34(20):6011-6019.(in Chinese)

[34] 王長庭,龍瑞軍,曹廣民,王啟蘭,丁路明,施建軍.三江源地區(qū)主要草地類型土壤碳氮沿海拔變化特征及其影響因素.植物生態(tài)學報,2006,30(3):441-449.

Wang C T,Long R J,Cao G M,Wang Q L,Ding L M,Shi J J.Soil carbon and nitrogen contents along elevation gradients in the source region of Yangtze,Yellow and Lantsang rivers.Journal of Plant Ecology,2006,30(3):441-449.(in Chinese)

[35] 趙娜,莊洋,趙吉.放牧和補播對草地土壤有機碳和微生物量碳的影響.草業(yè)科學,2014,31(3):367-374.

Zhao N,Zhang Y,Zhao J.Effects of grassland managements on soil organic carbon and microbial biomass carbon.Pratacultural Science,2014,31(3):367-374.(in Chinese)

[36] 李凌浩.土地利用變化對草原生態(tài)系統(tǒng)土壤碳貯量的影響.植物生態(tài)學報,1998,22(4):13-15.

Li L H.Effects of land-use change on soil carbon storage in grassland ecosystems.Acta Phytoecologica Sinica,1998,22(4):13-15.(in Chinese)

[37] 袁曉波,尚振艷,牛得草,傅華.黃土高原生態(tài)退化與恢復.草業(yè)科學,2015,32(3):363-371.

Yuan X B,Shang Z Y,Niu D C,Fu H.Advances in ecological degeneration and restoration of Loess Plateau.Pratacultural Science,2015,32(3):363-371．(in Chinese)

[38] 邢文淵,李大平,王蕾,石玉,肖然,肖繼東.巴里坤草原面積動態(tài)及其氣象因素分析.草業(yè)科學,2014(3):408-414.

Xing W Y,Li D P,Wang L,Shi Y,Xiao R,Xiao J D.Dynamic variation of grassland area and meteorological influence analysis in Barkol.Pratacultural Science,2014(3):408-414.(in Chinese)

[39] 孟凡喬,肖斐雯,孫向平,焦子偉,依力漢木,陳勇,吳文良,郭巖彬.新疆伊犁河流域有機產業(yè)發(fā)展?jié)摿Ψ治?新疆農業(yè)科學,2010,47(10):2037-2044.

Meng F Q,Xiao F W,Sun X P,Jiao Z W,Yilihanmu,Chen Y,Wu W L,Guo Y B.Analysis on development potentiality of organic industry in Ili river valley of Xinjiang.Xinjiang Agricultural Sciences,2010,47(10):2037-2044.(in Chinese)

[40] 范天文,薛發(fā)祥,范守民.伊犁河谷退耕還草的模式及生態(tài)經濟效益探討.新疆環(huán)境保護,2004,26(zk):42-45.

Fan T W,Xue F X,Fan S M.Approaching returning the grain plots to grass style of Yili river valley and ecology economy benefits.Environmental Protection of Xinjiang,2004,26(zk): 42-45.(in Chinese)

[41] 孫文義,邵全琴,劉紀遠,肖桐.三江源典型高寒草地坡面土壤有機碳變化特征及其影響因素.自然資源學報,2011,26(12):2072-2087.

Sun W Y,Shao Q Q,Liu J Y,Xiao T.Alpine slope grasslands and its influencing factors in the “Three-River Headwaters” region.Journal of Natural Resources,2011,26(12):2072-2087.(in Chinese)

[42] 李明明,韓鳳朋,劉恒博,龐國偉,張興昌.1992-2010年黃土高原小流域土壤有機碳時空變異性研究.干旱區(qū)資源與環(huán)境,2014,28(4):134-139.

Li M M,Han F P,Liu H B,Pang G W,Zhang X C.Analysis of temporal and spatial variability of soil organic carbon change in Loess Plateau small catchment:1992-2010.Journal of Arid Land Resources and Environment,2014,28(4):134-139.(in Chinese)

[43] 文雯,周寶同,汪亞峰,黃勇.基于輔助環(huán)境變量的土壤有機碳空間插值——以黃土丘陵區(qū)小流域為例.生態(tài)學報,2013,33(19):6389-6397.

Wen W,Zhou B T,Wang Y F,Huang Y.Soil organic carbon interpolation based on auxiliary environmental covariates:A case study at small watershed scale in Loss Hilly region.Acta Ecologica Sinica,2013,33(19):6389-6397.(in Chinese)

(責任編輯 王芳)

Spatial analysis of the soil total carbon in Ili,Xinjiang Uygur Autonomous Region, China

Zhou Li-lei1,2, Zhu Hua-zhong1, Zhong Hua-ping3, Yang Hua2

(1.State Key Laboratory of Resources and Environment Information System, Institute of Geographical Sciences and Natural Resources Researches, CAS, Beijing 100101,China; 2.Chongqing Key Laboratory of Earth Surface Processes and Environmental Remote Sensing in Three Gorges; Reservoir Area, College of Geographical and Tourism, Chongqing Normal University, Chongqing 400047, China； 3.Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographical Sciences and Natural Resources Research, CAS, Beijing 100101, China)

Currently, the studies of soil total carbon of grassland are mainly focused on traditional statistical analysis based on survey data and spatial analysis of soil total carbon are limited. The present study attempted to apply multi-source data including field survey soil total carbon data, remote sensing data and meteorological data obtained from 146 grassland sampling plots by using geostatistics to study the spatial analysis of soil total carbon data in the Ili region, located at the northwest margin of Xinjiang Uygur Autonomous Region. The results show that there was a significant correlation between the total carbon content from various depths (0-10, 10-20, 20-30, 30-40 cm and the mean of four soil depths) and altitude, average annual temperature, annual precipitation, ≥10 ℃ annual accumulated temperature, humidity and NDVI. A comprehensive evaluation model was established using this six ecological factors, and the spatial distribution pattern of grassland soil total carbon content was inversed with 1 km spatial resolution. The soil total carbon content in 0-40 cm depth was 4.2519% and it decreased with the increase of soil depth. The results indicated that the spatial distribution of soil total carbon content in different depths significantly varied in areas of Ili. The total carbon content in high altitude region was higher than that of the river valley. The highest total carbon content was distributed on northern slopes of the Southern Tianshan mountains and southern slopes of Boluokenu mountains, and the lowest was found in intensive human activity areas, distributed in both sides of Ili river and interchange of the Teks River and the Ili River. The results indicated that the spatial distribution of soil total carbon content was correlated with the types of grassland in Ili, soil total carbon content significantly influenced by human activities, i.e. overgrazing, blind land reclamation.

soil total carbon; spatial pattern; geostatistics; grassland; Ili region

Zhong Hua-ping E-mail:zhonghp@igsnrr.ac.cn

10.11829/j.issn.1001-0629.2015-0745

2015-12-30 接受日期：2016-05-10

科技部國家科技基礎性工作專項(2012FY111900-2、2011FY110400-3)；國家科技基礎條件平臺——地球系統(tǒng)科學數據共享平臺(2005DKA32300)第一作者：周李磊(1989-)，男，江蘇東海人，在讀博士生，主要從事資源環(huán)境遙感與GIS應用研究。E-mail:zhoull_001@126.com

鐘華平(1964-)，男，江西興國人，副研究員，碩士，主要從事草地資源與草地生態(tài)研究。E-mail:zhonghp@igsnrr.ac.cn

S812.2

A

1001-0629(2016)10-1963-12*

周李磊,朱華忠,鐘華平,楊華.新疆伊犁地區(qū)草地土壤全碳含量空間格局分析.草業(yè)科學,2016,33(10)：1963-1974.

Zhou L L,Zhu H Z,Zhong H P,Yang H.Spatial analysis of the soil total carbon in Ili,Xinjiang Uygur Autonomous Region,China.Pratacultural Science，2016,33(10)：1963-1974.