張立勇， 高鵬?， 王成軍， 劉勝濤， 李肖

(1.山東農(nóng)業(yè)大學(xué)林學(xué)院,國家林業(yè)局泰山森林生態(tài)站, 271018, 山東泰安; 2.山東省臨朐縣水利局, 262600, 山東臨朐)

魯中南藥鄉(xiāng)小流域林地土壤有機碳空間分布特征

張立勇1， 高鵬1?， 王成軍2， 劉勝濤1， 李肖1

(1.山東農(nóng)業(yè)大學(xué)林學(xué)院,國家林業(yè)局泰山森林生態(tài)站, 271018, 山東泰安; 2.山東省臨朐縣水利局, 262600, 山東臨朐)

土壤有機碳庫是全球碳循環(huán)和氣候變化研究的核心內(nèi)容，以魯中南山區(qū)藥鄉(xiāng)小流域林地土壤為研究對象，應(yīng)用地統(tǒng)計學(xué)和GIS技術(shù)，研究林地土壤有機碳的空間異質(zhì)性分布特征。結(jié)果表明：1)土壤有機碳和全氮、全磷質(zhì)量分數(shù)隨土壤層深度增加而減小，土壤有機碳質(zhì)量分數(shù)在0～10 和10～30 cm土層之間差異極顯著；全氮和全磷質(zhì)量分數(shù)在0～10和10～30 cm土層之間差異顯著。2)0～10 cm土層的有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈極顯著正相關(guān)，與全磷質(zhì)量分數(shù)呈顯著正相關(guān)；在10～30 cm土層有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈顯著正相關(guān)。3)土壤有機碳質(zhì)量分數(shù)的半方差函數(shù)模型屬于高斯模型，0～10 cm土層有機碳質(zhì)量分數(shù)具有強烈的架構(gòu)化區(qū)域模式；10～30 cm土層有機碳質(zhì)量分數(shù)表現(xiàn)為中等強度的空間自相關(guān)性。4)土壤有機碳質(zhì)量分數(shù)是流域的東部和西北部較高、中南部較低；有機碳質(zhì)量分數(shù)分布由東部和西北部向中南部遞減，基本與藥鄉(xiāng)小流域DEM趨勢一致。5)0～10 cm土層的土壤有機碳質(zhì)量分數(shù)分形維數(shù)較小，土壤有機碳空間格局較簡單；10～30 cm土層的有機碳質(zhì)量分數(shù)分形維數(shù)較大，土壤有機碳空間格局較為復(fù)雜。不同方位的土壤有機碳質(zhì)量分數(shù)分形維數(shù)表明，東-西方向是藥鄉(xiāng)小流域林地土壤有機碳分布的優(yōu)勢格局。

土壤有機碳； 空間異質(zhì)性； 地統(tǒng)計學(xué)； GIS； 魯中南山區(qū)

土壤有機碳作為陸地生態(tài)系統(tǒng)中重要的生態(tài)因子，是評價土壤肥力和土壤質(zhì)量的重要指標[1-2]，尤其是林地土壤有機碳庫是全球碳循環(huán)的重要組成部分，是全球碳循環(huán)和氣候變化研究的核心內(nèi)容，它不僅為植被生長提供碳源，而且在很大程度上影響著土壤結(jié)構(gòu)、土壤抗侵蝕的物理穩(wěn)定性及土壤生物多樣性等[3-5]?？臻g異質(zhì)性是土壤有機碳的一個基本特征，近年來，應(yīng)用地統(tǒng)計學(xué)和GIS技術(shù)對小流域土壤養(yǎng)分空間變異特征的研究備受眾多學(xué)者關(guān)注[6-9]。張偉等[10]對西南喀斯特地區(qū)峰叢洼地土壤有機碳的空間分布進行了研究，并對其進行了預(yù)測分析；薛志婧等[11]在小流域尺度上研究了黃土丘陵區(qū)小流域土壤有機碳空間異質(zhì)性，并得出不同土地利用類型下的土壤有機碳存在顯著性差異；耿廣坡等[12-13]利用回歸克里格插值法對山東省馬蹄峪小流域土壤有機碳質(zhì)量分數(shù)空間分布特征進行了研究，并得出流域內(nèi)海拔高的區(qū)域土壤有機碳質(zhì)量分數(shù)相對較高的結(jié)論；劉玲等[14]對吉林省金倉林場土壤有機碳空間變異進行了研究，發(fā)現(xiàn)空間變異因土壤深度而異，表層土壤有機碳的空間變異性較大?？傊?，針對小流域土壤養(yǎng)分空間變異特征的研究主要集中在我國南部和西部，而且多為不同土地利用方式下土壤有機碳分布特征方面的研究；但是關(guān)于魯中南山區(qū)小流域尺度下土壤有機碳的分布，尤其是基于地統(tǒng)計學(xué)的小流域林地土壤空間變異與分布規(guī)律的研究成果鮮見報道。筆者以魯中南山區(qū)典型區(qū)域的藥鄉(xiāng)小流域林地土壤為研究對象，應(yīng)用地統(tǒng)計學(xué)和GIS技術(shù)，研究分析藥鄉(xiāng)小流域林地土壤有機碳的空間異質(zhì)性分布特征，并探討林地土壤有機碳的空間分布異質(zhì)性與土壤全氮、全磷和海拔的相關(guān)性，可為魯中南山區(qū)小流域林地土壤養(yǎng)分空間變異規(guī)律研究提供科學(xué)依據(jù)，同時對于小尺度的土壤有機碳空間變異研究，精確估算區(qū)域碳儲量和制訂區(qū)域水土保持與生態(tài)恢復(fù)措施具有重要意義。

1 研究區(qū)概況

研究區(qū)位于山東省泰安市城市水源地黃前流域藥鄉(xiāng)小流域，地理坐標為E 117°05′39″～117°09′26″，N 36°17′58″～36°20′30″，總面積4.20 km2。地處北方土石山區(qū)魯中南山地丘陵區(qū)西北部，屬暖溫帶大陸性半濕潤季風氣候，多年平均氣溫為12.6 ℃，≥10 ℃的活動積溫3 821 ℃，無霜期196 d，多年平均降水量758 mm，降水主要集中在6—9月。流域地形東、北、西三面環(huán)山，海拔629～940 m。侵蝕剝蝕巖溶構(gòu)造地貌，溝壑密度6.67 km/km2，主溝比降18.23‰。流域內(nèi)巖石主要為花崗片麻巖風化層。土壤類型主要是棕壤土，呈微酸性，pH值在6.0左右，土層較薄，平均厚度30 cm。流域內(nèi)林地植被均為人工林(表1)，主要樹種為赤松(PinusdensifloraSieb. et Zucc.)、麻櫟(QuercusacutissimaCarruth.)和刺槐(Robiniapseudoacacia)。

2 材料與方法

2.1 土樣采集與數(shù)據(jù)處理

根據(jù)流域地形地貌、植被類型和采樣點的代表性及可操作性，結(jié)合小流域1∶1萬地形圖生成DEM柵格圖(10 m×10 m)，于2013年11月，共設(shè)置樣地77個(圖1)，樣地面積為0.025 hm2。在樣地內(nèi)以S形布設(shè)5個采樣點，按0～10 cm和10～30 cm分別采集各層土壤的混合樣1 000 g左右；同時，用GPS記錄樣地的經(jīng)緯度及高程值，并詳細記錄采樣點周圍的景觀信息?；旌贤翗咏?jīng)風干、磨細、過篩(2 mm、0.25 mm和0.149 mm土壤篩)，用于土壤化學(xué)性質(zhì)的測定。土壤有機碳采用外加熱重鉻酸鉀氧化—容量法，全氮采用半微量凱氏法，全磷采用氫氧化鈉熔融—鋁銻抗比色法[15]。

表1 研究區(qū)典型植被的基本概況

圖1 藥鄉(xiāng)小流域林地土壤采樣位置布設(shè)圖Fig.1 Locations of soil sampling points (0-30 cm) distribution in the Yaoxiang small watershed

2.2 統(tǒng)計分析

1)描述性及相關(guān)性分析：應(yīng)用SPSS 19.0軟件進行統(tǒng)計分析和回歸分析。

2)地統(tǒng)計分析：根據(jù)小流域1∶1萬地形圖生成DEM柵格圖(10 m×10 m)來提取地形因子，坡向以正北方向為0°，并按順時針方向計為0°到360°。利用GS+9.0和Arcgis10.0軟件進行土壤有機碳的半方差函數(shù)分析和空間預(yù)測分析。半方差函數(shù)計算公式[16]為

(1)

式中：γ(h)為半方差函數(shù)；h為樣本間隔距離；z(xi)和z(xi+h)為變量在空間位置xi和xi+h的取值；N(h)為取樣間隔為h時的樣本對總數(shù)。

半方差函數(shù)反映了土壤相鄰采樣點之間的空間關(guān)系，并且把空間變異距離增加到一個恒定值，即基臺值(C0+C1)作為空間依賴范圍，當采樣點間距接近這個恒定值則具有空間相關(guān)性，如果超過此距離則空間相關(guān)性消失；因此，應(yīng)用方差函數(shù)模型分析了土壤有機碳的空間相關(guān)性。半方差函數(shù)理論模型包括球形模型、指數(shù)模型、高斯模型、純塊金效應(yīng)模型和線性模型。半方差模型提供變量的空間結(jié)構(gòu)和克里格插值的輸入?yún)?shù)。土壤變量的空間依賴性可由半方差函數(shù)得到塊金值(C0)、基臺值(C0+C1)、結(jié)構(gòu)比(C0/(C0+C1))、變程Range和決定系數(shù)R2等重要指標。半方差函數(shù)γ(h)的值隨距離h的增大，從非零值達到一個相對穩(wěn)定的常數(shù)，這個常數(shù)為基臺值，基臺值是系統(tǒng)屬性中最大的變異；當距離h為0時，γ(h)=C0，表示隨機部分的空間異質(zhì)性。當函數(shù)達到基臺值時的間隔距離成為變程，當距離h大于變程以后，區(qū)域化變量z(xi)的空間相關(guān)性會消失[17]。R2的數(shù)值越大表示土壤有機碳質(zhì)量分數(shù)分布的優(yōu)勢格局越明顯。

有機碳質(zhì)量分數(shù)的分形維數(shù)的計算公式[18]為

D=(4-M)/2。

(2)

式中：D為有機碳質(zhì)量分數(shù)的分形維數(shù)，表示樣本之間的結(jié)構(gòu)性，D值越大，樣本之間差異越小，即均一性較好；M為logγ(h)與logh進行線性回歸后的直線斜率。

3 結(jié)果與分析

3.1 描述性統(tǒng)計

藥鄉(xiāng)小流域林地土壤有機碳、全氮和全磷質(zhì)量分數(shù)描述性統(tǒng)計見表2，可看出：在0～10 cm，土壤有機碳質(zhì)量分數(shù)均值為20.38 g/kg；在10～30 cm，土壤有機碳質(zhì)量分數(shù)均值為9.36 g/kg，隨著土壤深度的增加，土壤有機碳質(zhì)量分數(shù)減小。在0～10 cm，土壤有機碳質(zhì)量分數(shù)變異系數(shù)為17.58%；在10～30 cm，土壤有機碳質(zhì)量分數(shù)變異系數(shù)為14.52%，它們都屬于中等程度變異。在0～10 cm，全氮質(zhì)量分數(shù)變異系數(shù)均超過40%，屬于中等程度變異，全磷質(zhì)量分數(shù)變異系數(shù)均小于10%，屬于弱度變異；在10～30 cm，全氮和全磷質(zhì)量分數(shù)變異系數(shù)均超過10%，屬于中等程度變異，而且隨土壤深度增加，土壤全氮、全磷質(zhì)量分數(shù)減小。

表2 藥鄉(xiāng)小流域林地土壤養(yǎng)分質(zhì)量分數(shù)的描述性統(tǒng)計表

注：表中大寫字母表示不同土層間極顯著差異(P<0.01)，小寫字母表示不同土層間顯著差異(P<0.05)。Note: The capital letters indicate a significant correlation at the 0.01 level, and the lowercase letters indicate a significant correlation at the 0.05 level.

方差分析結(jié)果表明：土壤有機碳質(zhì)量分數(shù)在0～10 cm和10～30 cm土層之間差異極顯著(P<0.01)；全氮和全磷質(zhì)量分數(shù)在0～10 cm和10～30 cm土層之間差異顯著(P<0.05)。這主要是因為，在0～10 cm土壤層有林地植被枯落物覆蓋，枯落物分解后，對其所在的土壤具有很好的改良效果，增加了土壤養(yǎng)分來源，從而導(dǎo)致林地土壤0～10 cm的有機碳、全氮和全磷質(zhì)量分數(shù)比10～30 cm的高。

3.2 相關(guān)性分析

藥鄉(xiāng)小流域林地不同土層土壤有機碳質(zhì)量分數(shù)與土壤性質(zhì)和海拔的相關(guān)性分析見圖2，可以看出：在0～10 cm土層，土壤有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈極顯著正相關(guān)，相關(guān)系數(shù)分別為0.711(P<0.01)和0.563(P<0.01)，與全磷質(zhì)量分數(shù)呈顯著正相關(guān)，相關(guān)系數(shù)為0.499(P<0.05)；在10～30 cm土層，有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈顯著正相關(guān)，相關(guān)系數(shù)分別為0.424(P<0.05)和0.351(P<0.05)，與全磷質(zhì)量分數(shù)相關(guān)性不顯著(P>0.05)。相關(guān)分析表明，藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)與土壤的全氮和全磷質(zhì)量分數(shù)以及流域的海拔間存在密切關(guān)系。

3.3 地統(tǒng)計分析

3.3.1 半方差分析 藥鄉(xiāng)小流域林地土壤有機碳的半方差函數(shù)模型分析結(jié)果見表3和圖3，可看出，該流域林地土壤的0～10 cm和10～30 cm土層有機碳質(zhì)量分數(shù)的半方差函數(shù)模型均屬于高斯模型，2土層半方差函數(shù)模型的參數(shù)值存在一定差異。0～10 cm土層半方差函數(shù)模型的結(jié)構(gòu)比為14.79%(<25%)，變程是706 m，說明具有強烈的架構(gòu)化區(qū)域模式；10～30 cm土層半方差函數(shù)模型的結(jié)構(gòu)比為30.00%(>25%)，變程是721 m，并表現(xiàn)為中等強度的空間自相關(guān)性。

3.3.2 土壤有機碳質(zhì)量分數(shù)空間分布特征及分形維數(shù) 利用普通克里格插值方法得到藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)的空間分布情況見圖4，可以看出：藥鄉(xiāng)小流域林地土壤的有機碳質(zhì)量分數(shù)東部和西北部較高，流域的中南部較低；有機碳質(zhì)量分數(shù)由東部和西北部向中南部遞減，基本與藥鄉(xiāng)小流域DEM趨勢一致，也就是說海拔高的區(qū)域土壤有機碳質(zhì)量分數(shù)相對較高。

圖2 藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)與土壤性質(zhì)和海拔的相關(guān)性分析Fig.2 Correlation analysis of SOC and different soil properties and elevation of forest land in the Yaoxiang small watershed

土層Soillayer/cm模型Model塊金值C0基臺值C0+C塊基比C0/(C0+C)/%變程Range/m0～10高斯模型Gaussianmodel1.770011.970014.7970610～30高斯模型Gaussianmodel0.52501.750030.00721

圖3 藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)的半方差函數(shù)模型Fig.3 Semivariogram model of SOC content of forest land in the Yaoxiang small watershed

表4是藥鄉(xiāng)小流域林地土壤不同土層有機碳質(zhì)量分數(shù)的分形維數(shù)，從全方位來看：0～10 cm土層的有機碳質(zhì)量分數(shù)分形維數(shù)較小(1.593)，空間格局較簡單，空間依賴性較強；10～30 cm土層的有機碳質(zhì)量分數(shù)分形維數(shù)較大(1.762)，均一程度高，土壤有機碳質(zhì)量分數(shù)空間格局較為復(fù)雜。比較不同方位的分形維數(shù)：0～10 cm土層的土壤有機碳質(zhì)量分數(shù)在各方位分形維數(shù)差異不大，其中西北—東南方向上大，在東—西方向分形維數(shù)小，在此方向上的空間異質(zhì)性較大，具有較強的空間依賴性；10～30 cm土層的土壤有機碳質(zhì)量分數(shù)在東北—西南方向分形維數(shù)大，均一性較好，在南—北方向分形維數(shù)小，在該方向上空間異質(zhì)性和空間依賴性較強。同時，從決定系數(shù)來看，0～10和10～30 cm土層都是在東—西方向大(0.866和0.547)，說明東—西方向是該流域林地土壤有機碳質(zhì)量分數(shù)分布的優(yōu)勢格局。

圖4 藥鄉(xiāng)小流域不同土層土壤有機碳質(zhì)量分數(shù)普通克里格插值圖Fig.4 Spatial distributions of different SOC content interpolated by ordinary-kriging in the Yaoxiang small watershed

表4 藥鄉(xiāng)小流域林地土壤各層有機碳質(zhì)量分數(shù)的分形維數(shù)

注：S-N為南-北方向；NE-SW為東北-西南方向；E-W為東-西方向；NW-SE為西北-東南方向；D為分形維數(shù)。Note: S-N is the south- north direction; NE-SW is the northeast-southwest direction; E-W is the East-west direction; NW-SE is the Northwest-Southeast direction.Dis the fractal dimension of SOC content.

4 討論

4.1 不同土層有機碳質(zhì)量分數(shù)的變化特征

藥鄉(xiāng)小流域0～10 cm土層有機碳質(zhì)量分數(shù)高于10～30 cm土層，二者之間差異極顯著；半方差函數(shù)模型分析表明，0～10 cm土層有機碳質(zhì)量分數(shù)具有強烈的架構(gòu)化區(qū)域模式，10～30 cm土層有機碳質(zhì)量分數(shù)表現(xiàn)為中等強度的空間自相關(guān)性，該研究結(jié)果與劉玲[14]和孫文義等[19]的研究結(jié)果一致。這主要與凋落物和植物根系關(guān)系密切，表層土壤覆蓋有大量的森林凋落物，而凋落物是森林生態(tài)系統(tǒng)物質(zhì)循環(huán)過程中的一個重要物質(zhì)庫，是表土碳、氮元素的主要補給者，調(diào)落物的數(shù)量和分解過程對表層土壤碳、氮的累積有直接的影響，同時，凋落物淋洗水所攜帶的有機碳量是土壤有機碳輸入的較大項。藥鄉(xiāng)小流域在0～10 cm土層有較多的植被凋落物覆蓋，凋落物分解后，對其所在的土壤具有很好的改良效果，增加了土壤養(yǎng)分來源，從而導(dǎo)致林地土壤0～10 cm的有機碳、全氮和全磷的質(zhì)量分數(shù)均比10～30 cm土層的高。

4.2 不同海拔有機碳質(zhì)量分數(shù)的異質(zhì)性變化特征

藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)與流域的海拔間存在顯著的相關(guān)性，總體分布特征是東部和西北部較高，流域的中南部較低；有機碳質(zhì)量分數(shù)分布由東部和西北部向中南部遞減，基本與藥鄉(xiāng)小流域DEM走向一致，即海拔高的區(qū)域土壤有機碳質(zhì)量分數(shù)相對較高。這主要是由于該地區(qū)經(jīng)過多年的封山育林，高海拔地區(qū)的森林植被受到人為的干擾程度較輕，林下枯枝落葉物保存的比較完整，從而較好地改良了土壤結(jié)構(gòu)，土壤有機碳質(zhì)量分數(shù)增加，這一研究結(jié)論與T. Hengl等[20]、R. K. Gupta等[21]和李海東等[22]的研究結(jié)論吻合。可見，保護和恢復(fù)好森林植被，盡可能減少人為干擾和破壞，對森林土壤有機碳固存具有重要的意義。

本研究采用普通克里格插值法預(yù)測并分析了小流域尺度上林地土壤有機碳質(zhì)量分數(shù)分布特點，但由于受到地形、坡向、植被等綜合環(huán)境因子的影響，預(yù)測精度存在一定的局限性。同時，本研究僅對藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)進行了分析，未進行有機碳儲量估算，仍還需進一步研究。

5 結(jié)論

1) 藥鄉(xiāng)小流域林地土壤有機碳和全氮、全磷質(zhì)量分數(shù)隨土壤層深度增加而減小，土壤有機碳0～10 cm和10～30 cm土層間差異極顯著(P<0.01)；全氮和全磷質(zhì)量分數(shù)0～10 cm和10～30 cm土層間差異顯著(P<0.05)。

2) 藥鄉(xiāng)小流域林地0～10 cm土層的有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈極顯著正相關(guān)關(guān)系，與全磷質(zhì)量分數(shù)呈顯著正相關(guān)關(guān)系；在10～30 cm土層有機碳質(zhì)量分數(shù)與全氮質(zhì)量分數(shù)、海拔呈顯著正相關(guān)關(guān)系。

3) 藥鄉(xiāng)小流域林地土壤有機碳質(zhì)量分數(shù)的半方差函數(shù)模型屬于高斯模型，0～10 cm土層有機碳質(zhì)量分數(shù)具有強烈的架構(gòu)化區(qū)域模式，10～30 cm土層有機碳質(zhì)量分數(shù)表現(xiàn)為中等強度的空間自相關(guān)性。

4) 藥鄉(xiāng)小流域林地土壤的有機碳質(zhì)量分數(shù)東部和西北部較高，流域的中南部較低；有機碳質(zhì)量分數(shù)分布由東部和西北部向中南部遞減，基本與藥鄉(xiāng)小流域DEM趨勢一致，即海拔高的區(qū)域土壤有機碳質(zhì)量分數(shù)相對較高。

5) 藥鄉(xiāng)小流域林地土壤0～10 cm土層的有機碳質(zhì)量分數(shù)分形維數(shù)較小，空間格局較簡單，空間依賴性較強；10～30 cm土層的有機碳質(zhì)量分數(shù)分形維數(shù)較大，均一程度高，土壤有機碳質(zhì)量分數(shù)空間格局較為復(fù)雜。不同方位的土壤有機碳質(zhì)量分數(shù)分形維數(shù)表明東-西方向是該流域林地土壤有機碳質(zhì)量分數(shù)分布的優(yōu)勢格局。

[1] 程先富, 史學(xué)正, 于東升, 等. 江西省興國縣土壤全氮和有機質(zhì)的空間變異及其分布格局[J].應(yīng)用與環(huán)境生物學(xué)報,2004, 10(1):64-67

[2] Lal R. Carbon sequestration in dry lands [J]. Annual of Arid Zone, 2000, 39 (1): 1-10

[3] 藍良就, 張野, 黃炎和,等.退化花崗巖植被恢復(fù)對團聚體及其有機碳的影響[J].水土保持學(xué)報,2012,26(1):190-194

[4] 耿玉清, 孫向陽, 亢新剛,等.長白山林區(qū)不同森林類型下的土壤肥力狀況的研究[J].北京林業(yè)大學(xué)學(xué)報,1999,21(6):97-101

[5] 張超, 劉國彬, 薛萐,等. 黃土丘陵區(qū)不同林齡人工刺槐林土壤抗蝕性演變特征[J]. 中國水土保持科學(xué), 2010, 8(2): 1-7

[6] Kucharik C J, Brye K R, Norman J M, et al. Measurements and modeling of carbon and nitrogen cycling in agro-ecosystems of southern Wisconsin: Potential for SOC sequestration during the next 50 years[J]. Ecosystems, 2001, 4: 237-258

[7] 王清奎, 汪思龍, 高洪, 等. 土地利用方式對土壤有機質(zhì)的影響[J]. 生態(tài)學(xué)雜志,2005,24(4): 360-363

[8] Stacey K F, Lark R M, Whitmore A P, et al. Using a process model and regression kriging to improve predictions of nitrous oxide emissions from soil[J].Geoderma,2006,135:107-117

[9] 張世熔, 黃元仿, 李保國, 等. 黃淮海沖積平原區(qū)土壤有機質(zhì)時空變異特征[J]. 生態(tài)學(xué)報,2002,22(12): 2041-2047

[10] 張偉, 王克林, 陳洪松, 等. 典型喀斯特峰叢洼地土壤有機碳含量空間預(yù)測研究[J]. 土壤學(xué)報,2012,49(3): 601-606

[11] 薛志婧, 侯曉瑞, 程曼, 等. 黃土丘陵區(qū)小流域尺度上土壤有機碳空間異質(zhì)性[J]. 水土保持學(xué)報,2011,25(3): 160-163

[12] 耿廣坡, 高鵬, 呂圣橋, 等. 魯中南山區(qū)馬蹄峪小流域土壤有機質(zhì)和全氮空間分布特征[J]. 中國水土保持科學(xué),2011,9(6): 99-105

[13] 耿廣坡.基于回歸克里格模型的小流域土壤養(yǎng)分空間分布特征及影響因素[D]. 山東泰安:山東農(nóng)業(yè)大學(xué),2012:20-25

[14] 劉玲.長白山低山區(qū)典型森林類型土壤有機碳及養(yǎng)分空間異質(zhì)性研究[D]. 北京: 北京林業(yè)大學(xué),2013:40-43

[15] 魯如坤.土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社,1999:16-39

[16] 姜勇, 梁文舉, 李琪. 利用與回歸模型相結(jié)合的克里格方法對農(nóng)田土壤有機碳的估值及制圖[J]. 水土保持學(xué)報,2005,19(5): 97-126

[17] Sumfleth K, Duttmann R. Prediction of soil property distribution in paddy soil landscapes using terrain data and satellite information as indicators[J]. Ecological Indicators, 2008,8 (5): 485-501

[18] 龔元石, 廖超子, 李保國. 土壤含水量和容重的空間變異及其分形特征[J]. 土壤學(xué)報,1998,35(1): 10-15

[19] 孫文義, 郭勝利, 宋小燕. 地形和土地利用對黃土丘陵溝壑區(qū)表層土壤有機碳空間分布影響[J]. 自然資源學(xué)報,2010,25(3): 443-453

[20] Hengl T, Heuvelink G B M, Stein A. A generic framework for spatial prediction of soil variables based on regression-kriging[J]. Geoderma, 2004,120: 75-93

[21] Gupta R K,Vijayan D, Prasad T S. New hyperspectral vegetation characterization parameters[J]. Advances in Space Research,2001,28: 201-206

[22] 李海東, 林杰, 張金池, 等. 小流域尺度下土壤有機碳和全氮空間變異特征[J]. 南京林業(yè)大學(xué)學(xué)報: 自然科學(xué)版,2008,32(4): 38-42

(責任編輯：程 云)

Spatial distribution of soil organic carbon in the forestland of the Yaoxiang small watershed in central and southern Shandong Province

Zhang Liyong1, Gao Peng1, Wang Chengjun2, Liu Shengtao1, Li Xiao1

(1.Shandong Agricultural University, College of Forestry; Taishan Mountain Forest Ecosystem Research Station, 271018, Tai′an, Shandong, China; 2.Water Resources Bureau of Linqu County, 262600, Linqu, Shandong, China)

The spatial heterogeneity of soil organic carbon (SOC) levels is important in researches of both global carbon-nitrogen cycle and climate change. By employing geostatistics methods and GIS technology, we studied the spatially heterogeneous distribution of SOC in the Yaoxiang small watershed in the hilly area of central and southern Shandong Province. Results showed that: 1) SOC, total nitrogen and total phosphorus contents showed a trend of declining with the increase of soil depth, and there was an extremely significant difference in SOC content between 0-10 cm and 10-30 cm soil layers, and significant differences in total nitrogen and total phosphorus contents were found between 0-10 cm and 10-30 cm soil layers. 2) SOC content in 0-10 cm soil layer had an extremely significantly positive correlation with total nitrogen content and elevation, and had a significantly positive correlation with total phosphorus content; SOC content in 10-30 cm soil layer had a significantly positive correlation with total nitrogen and elevation. 3) Semivariogram model of SOC content belongs to the Gaussian model, and SOC content in 0-10 cm soil layer had a strong regional pattern, and that in 10-30 cm soil layer showed a moderate spatial correlation. 4) SOC content was higher in the eastern and northwestern parts and lower in the central and southern parts of the small watershed; the distribution of SOC diminished gradually from the eastern and northwestern parts to central and southern parts, which was basically consistent with the trend of DEM in the Yaoxiang small watershed. 5) Fractal dimension of SOC content in 0-10 cm soil layer was smaller than that in 10-30 cm soil layer, and the spatial pattern of SOC was relatively simple; in contrast, the spatial pattern of SOC in 10-30 cm soil layer was more complex. The fractal dimensions of SOC content in different directions showed that the east-west (E-W) direction was the advantageous pattern of SOC distribution in the Yaoxiang small watershed.

SOC; spatial heterogeneity; geostatistics; GIS; hilly area of central and southern Shandong Province

2014-09-08

2015-04-06

張立勇(1988—)，男，碩士研究生。主要研究方向：流域治理與生態(tài)修復(fù)。E-mail:zhangliyong@yeah.net

?通信作者簡介： 高鵬(1967—)，男，教授，博士。主要研究方向：林業(yè)土壤生態(tài)。E-mail：gaopengy@163.com

S158.2； S153.6

A

1672-3007(2015)03-0083-07

項目名稱： 國家林業(yè)公益性行業(yè)專項課題“森林生態(tài)系統(tǒng)水文和氣象生態(tài)要素連清技術(shù)研究”(2014303- 08)；“全國第八次山東省森林資源生態(tài)服務(wù)功能評估研究”(201204101)