賴晶晶，張世熔，劉艷婷，李 婷，徐小遜，姚 蘋，蒲玉琳

(1 四川農(nóng)業(yè)大學(xué)環(huán)境學(xué)院，成都 611130；2 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130)

?

基于多元分析的岷江下游土壤有機(jī)碳密度空間分布及影響因素研究①

賴晶晶1,2，張世熔1*，劉艷婷1，李 婷2，徐小遜1，姚 蘋1,2，蒲玉琳2

(1 四川農(nóng)業(yè)大學(xué)環(huán)境學(xué)院，成都 611130；2 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130)

摘 要：采用傳統(tǒng)統(tǒng)計(jì)、灰色關(guān)聯(lián)度、逐步回歸和通徑分析相結(jié)合的多元分析方法對(duì)岷江下游1 138個(gè)樣點(diǎn)的表層(0～20 cm)土壤有機(jī)碳密度及其影響因素進(jìn)行研究。結(jié)果表明，該區(qū)表層土壤有機(jī)碳密度在0.89～6.49 kg/m2之間，均值為3.24 kg/m2。其在空間上呈不規(guī)則的斑塊狀分布，具中等空間相關(guān)性，表現(xiàn)出中部淺丘區(qū)高而西北和東部高丘區(qū)低的分布趨勢(shì)。土地利用類型、土壤酸堿度(pH)、成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地均是影響土壤有機(jī)碳密度的主要因素。其中，土地利用類型、pH、地貌類型、坡度和土壤質(zhì)地與有機(jī)碳密度呈極顯著相關(guān)?；疑P(guān)聯(lián)度的分析表明與土壤有機(jī)碳密度關(guān)聯(lián)最為密切的是地貌類型，其次是pH和成土母質(zhì)，逐步回歸和通徑分析的結(jié)果均表明地貌類型是最主要因素，成土母質(zhì)和pH次之。綜合比較以上方法，得出地貌類型、成土母質(zhì)和pH為影響岷江下游研究區(qū)表層土壤有機(jī)碳密度的主要因素。

關(guān)鍵詞：有機(jī)碳密度；影響因素；灰色關(guān)聯(lián)度；逐步回歸；通徑分析

土壤碳庫是陸地生態(tài)系統(tǒng)中最大的碳庫，土壤有機(jī)碳是土壤碳庫的重要組成部分[1]。土壤有機(jī)碳不僅直接影響土壤肥力和作物產(chǎn)量，其與大氣中碳的周轉(zhuǎn)對(duì)全球碳平衡格局及氣候變化也有重要影響[2]。通過研究表層土壤有機(jī)碳密度的空間分布和影響因素可以更有效地為土壤有機(jī)碳的評(píng)價(jià)提供科學(xué)依據(jù)[3-4]，這對(duì)環(huán)境監(jiān)測(cè)、生產(chǎn)管理和預(yù)測(cè)陸地生態(tài)系統(tǒng)與氣候變化之間的反饋效應(yīng)都有一定的參考價(jià)值[5]。

近年來，灰色關(guān)聯(lián)度法因其計(jì)算量小且適用樣本量少或無規(guī)律的樣本數(shù)據(jù)，已在土壤質(zhì)量評(píng)價(jià)方面得以應(yīng)用[6-7]；逐步回歸分析是通過明確各因素影響研究主體的相對(duì)重要性篩選出主要因素，構(gòu)造“最優(yōu)”回歸方程的一種方法[8]；而通徑分析是回歸分析的拓展，可用于分析多個(gè)自變量與因變量的線性關(guān)系，判斷回歸模型是否存在多重共線性及其多重共線性的來源和強(qiáng)弱，明確各影響因素間的相互重要性[9]。采用傳統(tǒng)統(tǒng)計(jì)、灰色關(guān)聯(lián)度、逐步回歸和通徑分析相結(jié)合的多元分析方法研究土壤有機(jī)碳密度的空間分布和影響因素，將更全面地揭示其空間變異規(guī)律。

目前，已有研究采用傳統(tǒng)統(tǒng)計(jì)和地統(tǒng)計(jì)學(xué)相結(jié)合的方法對(duì)土壤有機(jī)碳密度的影響因素進(jìn)行探討[10-12]，但在流域尺度下將傳統(tǒng)統(tǒng)計(jì)、灰色關(guān)聯(lián)度、逐步回歸和通徑分析相結(jié)合分析土壤有機(jī)碳密度影響因素的研究還鮮見報(bào)道。岷江下游地區(qū)以丘陵地貌為主，人口稠密，是四川經(jīng)濟(jì)發(fā)展速度較快的區(qū)域之一。本文采用多元分析方法對(duì)該流域表層土壤有機(jī)碳密度的影響因素進(jìn)行分析，揭示其空間變異規(guī)律，以期為該區(qū)的土壤生產(chǎn)管理和環(huán)境監(jiān)測(cè)提供基礎(chǔ)數(shù)據(jù)和決策依據(jù)，也為綜合分析表層土壤有機(jī)碳密度影響因素提供參考。

1 材料與方法

1.1 研究區(qū)概況

岷江下游流域(103°40′～104°03′E，30°06′～30°27′N)位于四川盆地與川西北丘狀高原山地過渡地帶前緣，地處古隆中新陷雁行褶皺帶內(nèi)。該區(qū)地勢(shì)西北高、東南低，岷江由北向南縱貫區(qū)境，東西兩側(cè)多為圓狀淺丘，中部構(gòu)成沖積平原。其地貌兩側(cè)是丘陵，中部為平原。區(qū)內(nèi)東西兩側(cè)土壤由侏羅系和白堊系頁巖、砂巖坡殘積物發(fā)育而成，中部土壤主要由第四系河流沖積物和老沖積物發(fā)育而成。土壤類型主要有水稻土、潮土、紫色土和黃壤。土地利用類型主要為耕地、林地和園地，其中耕地以水田和旱地為主。區(qū)內(nèi)作物主要是水稻、玉米、小麥、油菜，為岷江下游的典型農(nóng)業(yè)區(qū)。區(qū)域?qū)賮啛釒駶櫦撅L(fēng)氣候，冬無嚴(yán)寒，夏無酷暑，年均溫16.8℃，無霜期長，年平均降水量955.2 mm，年日照時(shí)數(shù)1 106.8 h。研究區(qū)域面積為705 km2，境內(nèi)有一條鐵路，兩條高速公路(成雅、成樂)和一條省級(jí)公路。

1.2 研究方法

1.2.1 土壤樣品采集及分析 結(jié)合研究區(qū)衛(wèi)星影像圖、地形圖、成土母質(zhì)分布圖、第二次土壤普查的土壤類型分布圖及其67個(gè)原樣點(diǎn)分布圖和剖面記載表等資料，在考慮樣點(diǎn)代表性和空間分布均勻性的基礎(chǔ)上，按網(wǎng)格法和分層抽樣法進(jìn)行區(qū)域尺度下的樣點(diǎn)設(shè)計(jì)，于2007年調(diào)查樣點(diǎn)1 138個(gè)(圖1)。采樣同時(shí)記錄每個(gè)采樣點(diǎn)的地理坐標(biāo)、海拔、土壤類型、土地利用方式等信息。每個(gè)樣點(diǎn)土樣由同一地塊的3個(gè)表層土壤(0～20 cm)樣品混合組成。土壤樣品經(jīng)自然風(fēng)干，去掉植物殘?bào)w和石塊，研磨過2 mm篩，采用重鉻酸鉀容量法－外加熱法測(cè)定其有機(jī)碳含量。

圖1　樣點(diǎn)分布圖Fig.1 Distribution map of sampling points

1.2.2 有機(jī)碳密度計(jì)算 某一土層i的有機(jī)碳密度(SOCDi，kg/m2)的計(jì)算公式為[13]：

式中：δ% 為直徑大于2 mm的石礫含量(體積百分比)，Ci為土壤有機(jī)碳含量(g/kg)，Bi為土壤體積質(zhì)量(g/cm3)，Li為土層厚度(cm)。

1.2.3 數(shù)據(jù)統(tǒng)計(jì)分析 運(yùn)用SPSS19.0軟件對(duì)研究區(qū)1 138個(gè)表層土壤樣本的SOCD進(jìn)行描述性統(tǒng)計(jì)，采用K-S檢驗(yàn)分析土壤有機(jī)碳密度分布的正態(tài)性，通過方差分析、相關(guān)回歸分析、逐步回歸分析和通徑分析，研究岷江下游流域土壤有機(jī)碳密度實(shí)測(cè)值與土地利用類型、土壤酸堿度、成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地等因素及各因素之間的相關(guān)性。其中平均數(shù)的比較采用最小顯著差數(shù)法(LSD)。土地利用類型、成土母質(zhì)、地貌類型和土壤質(zhì)地為定性分類變量，研究中采用啞變量進(jìn)行賦值。在ArcGIS 10.0平臺(tái)上繪制樣點(diǎn)分布圖，利用DEM計(jì)算獲取坡度數(shù)據(jù)。在ArcGIS 10.0平臺(tái)上，制作樣點(diǎn)分布圖；利用DEM計(jì)算坡度；進(jìn)行趨勢(shì)效應(yīng)分析，獲得最佳變異函數(shù)模型，通過普通克里格插值獲取土壤有機(jī)碳密度空間分布圖。

1.2.4 灰色關(guān)聯(lián)度分析 灰色關(guān)聯(lián)度分析根據(jù)序列曲線幾何形狀的相似程度來判斷其聯(lián)系是否緊密，通過定量各序列之間的聯(lián)系，確定相關(guān)程度最大的因素。序列作為系統(tǒng)行為的映射量，間接地表征系統(tǒng)行為。本研究采用區(qū)域土壤有機(jī)碳密度實(shí)測(cè)值作為參考序列Xo，將相關(guān)影響因素的賦值作為比較序列Xi，計(jì)算公式如下[14]：

1.2.5 逐步回歸分析 土壤有機(jī)碳密度受多個(gè)因素的共同影響，但眾多影響因素并非都能進(jìn)入線性回歸模型，逐步回歸分析將每個(gè)影響因素不斷進(jìn)入回歸方程，再次判斷是否存在刻意剔除出方程的影響因素，從而篩選出對(duì)有機(jī)碳密度影響最大的因素。多元線性回歸的數(shù)學(xué)模型是：

式中：a為模型中的未知參數(shù)，B為隨機(jī)誤差。

1.2.6 通徑分析 通徑分析通過分解自變量與因變量之間的表面直接相關(guān)性，研究自變量對(duì)因變量的直接和間接重要性，是一種用于分析多個(gè)自變量與因變量之間的線性關(guān)系的多元統(tǒng)計(jì)分析方法，是回歸分析的拓展[15]。

式中：Riy為自變量Xi與Y之間的簡單相關(guān)系數(shù)，Piy為Xi與Y之間的直接通徑系數(shù)，即Xi與Y標(biāo)準(zhǔn)化后的偏相關(guān)系數(shù)，表示Xi對(duì)Y的直接影響效應(yīng)；Iiy為間接通徑系數(shù)，表示Xi通過其他變量對(duì)因變量Y的總間接影響效應(yīng)。

2 結(jié)果與討論

2.1 土壤有機(jī)碳密度的描述性統(tǒng)計(jì)特征

運(yùn)用SPSS19.0軟件對(duì)土壤樣本的有機(jī)碳密度值進(jìn)行描述性統(tǒng)計(jì)，經(jīng)單一樣本K-S檢驗(yàn)，結(jié)果表明研究區(qū)1 138個(gè)樣點(diǎn)的表層土壤有機(jī)碳密度呈正態(tài)分布(P ＞ 0.05)，從偏度值和峰度值可知其分布呈右偏態(tài)且較為平緩，土壤有機(jī)碳密度值在0.89～6.49 kg/m2，均值為3.24 kg/m2(表1)。根據(jù)20世紀(jì)80年代第二次土壤普查的有機(jī)質(zhì)含量分級(jí)標(biāo)準(zhǔn)，該區(qū)表層土壤有機(jī)碳含量屬于中上水平(20.1～30.0 g/kg)。

表1　土壤有機(jī)碳密度統(tǒng)計(jì)特征Table 1 Descriptive statistical characteristics of soil organic carbon density

2.2 空間分布特征

利用ArcGIS10.0軟件的地統(tǒng)計(jì)分析功能對(duì)土壤有機(jī)碳密度進(jìn)行趨勢(shì)效應(yīng)分析，結(jié)果表明研究區(qū)土壤有機(jī)碳密度分布在南-北方向和東-西方向皆呈二階趨勢(shì)效應(yīng)。通過變異函數(shù)的套合模型模擬，得到的最優(yōu)模型為指數(shù)模型(表2)。結(jié)果表明，該區(qū)土壤有機(jī)碳密度的分布長軸方向約為北-南向，具有顯著的異向性特征，隨機(jī)變異所占比例為47.96%，土壤有機(jī)碳密度具有中等空間相關(guān)性，空間變異以隨機(jī)變異為主。采用指數(shù)模型擬合結(jié)果的參數(shù)進(jìn)行普通克里格插值，得到土壤有機(jī)碳密度的空間分布圖(圖2)。結(jié)果表明，土壤有機(jī)碳密度總體上呈不規(guī)則的斑塊狀分布，高值區(qū)(＞ 4.00 kg/m2)主要位于中部低丘區(qū)，低值區(qū)(＜2.50 kg/m2)位于西北和東部高丘區(qū)域，這表明其受結(jié)構(gòu)性因素和隨機(jī)性因素的雙重影響。

表2　土壤有機(jī)碳密度半方差函數(shù)模型及相關(guān)參數(shù)Table 2 Theory model and parameters of sermi-variogram for soil organic carbon density

圖2　土壤有機(jī)碳密度空間分布圖Fig.2 Spatial distribution of topsoil organic carbon density

2.3 影響因素分析

土壤有機(jī)碳密度的變化實(shí)質(zhì)是其碳輸入和輸出動(dòng)態(tài)平衡的結(jié)果，它受結(jié)構(gòu)性因素和隨機(jī)性因素的共同影響，其影響因素隨空間尺度的不同而改變[16]。岷江下游地區(qū)社會(huì)經(jīng)濟(jì)和農(nóng)業(yè)生產(chǎn)發(fā)展迅速，土地利用方式和土壤酸堿度是自然條件和人類活動(dòng)的綜合反映，直接影響土壤有機(jī)碳密度的分布。區(qū)內(nèi)地貌以丘陵為主，樣點(diǎn)在低丘、中丘和高丘皆有分布，成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地都與該區(qū)土壤有機(jī)碳密度有一定關(guān)系。本文根據(jù)研究區(qū)的實(shí)際情況主要探討土地利用類型、土壤酸堿度、成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地對(duì)土壤有機(jī)碳密度的影響。

2.3.1 土地利用類型 土地利用類型主要通過控制有機(jī)碳的輸入來源及其積累和分解速率影響土壤有機(jī)碳密度[17]。研究區(qū)4種不同土地利用類型的土壤其有機(jī)碳密度存在極顯著差異(P＜0.05)，主要表現(xiàn)為水田的土壤有機(jī)碳密度(4.68 kg/m2)顯著高于菜地、果園和旱地，但菜地、果園和旱地之間均無顯著差異(P＞0.05)。這是由于該區(qū)菜地、果園和旱地因在輪作方式下土壤透水通氣性好，長期處于好氧狀態(tài)下土壤有機(jī)質(zhì)易于分解[13]，因此盡管其有機(jī)肥施用量較高，但土壤有機(jī)碳不易累積；水田因長期處于淹水的嫌氣狀態(tài)，土壤有機(jī)質(zhì)分解緩慢，因而有機(jī)碳易于累積使其含量較高。菜地、果園和旱地在區(qū)內(nèi)分布較集中，其土壤理化性質(zhì)、地表植被覆蓋情況、人為干擾程度相似，土壤有機(jī)碳密度無顯著差異[18]。

2.3.2 土壤酸堿度 土壤酸堿度(pH)是土壤的重要化學(xué)性質(zhì)，對(duì)土壤微生物的活性和有機(jī)質(zhì)的形成及分解均有著重要作用，與土壤有機(jī)碳密度關(guān)系密切[19]。研究區(qū)土壤pH范圍在4.34～9.01，現(xiàn)將它們按 ＜5.0、5.0～6.5、6.5～7.5、7.5～8.5、＞ 8.5分為強(qiáng)酸性、酸性、中性、堿性、強(qiáng)堿性土壤。結(jié)果表明，pH與土壤有機(jī)碳密度呈極顯著相關(guān)(P＜0.001)，酸性和中性土壤之間的有機(jī)碳密度差異不顯著，但均顯著高于其他的pH范圍(表3)。

表3　不同pH、成土母質(zhì)和土壤質(zhì)地下土壤有機(jī)碳密度值與方差分析Table 3 Mean and variance analyses of soil organic carbon densities in different pH values,parent materials and textures

2.3.3 成土母質(zhì) 成土母質(zhì)作為土壤形成的物質(zhì)基礎(chǔ)直接影響著土壤的成土過程，是決定土壤結(jié)構(gòu)和有機(jī)碳最初含量與分布特征的主要因素，與土壤的通氣性、保肥保水性及耕作難易性關(guān)系密切，不同成土母質(zhì)的土壤其有機(jī)碳密度存在差異[20]。該區(qū)土壤主要由第四系成土母質(zhì)發(fā)育而成，統(tǒng)計(jì)結(jié)果表明(表3)，老沖積物和河流沖積物上發(fā)育的土壤其有機(jī)碳密度值差異不明顯，但均顯著高于砂質(zhì)泥巖和泥質(zhì)砂巖上發(fā)育的土壤(P ＜ 0.05)。產(chǎn)生這種差異的原因主要是該區(qū)砂質(zhì)泥巖和泥質(zhì)砂巖上發(fā)育的土壤保肥保水能力較差，有機(jī)質(zhì)易分解，土壤有機(jī)碳儲(chǔ)量少且難積累，而更新統(tǒng)老沖積物和全新統(tǒng)近代河流沖積物上發(fā)育的土壤養(yǎng)分含量豐富且質(zhì)地黏重，土壤黏粒及其吸附的有機(jī)質(zhì)所形成的團(tuán)粒結(jié)構(gòu)有利于有機(jī)碳的長期積累[21]，因而土壤有機(jī)碳密度較大。在相同母質(zhì)下，不同質(zhì)地間表層土壤有機(jī)碳密度也存在差異。例如，同為河流沖積物上發(fā)育的土壤，黏土表層土壤有機(jī)碳密度為3.63 kg/m2，顯著高于重壤(3.31 kg/m2)、中壤(3.11 kg/m2)、輕壤(3.04 kg/m2)和砂壤(3.02 kg/m2)，而中壤、輕壤和砂壤間無顯著差異(P＜0.05)。顯然，在同類母質(zhì)下，質(zhì)地亦可能影響土壤有機(jī)碳密度。

2.3.4 地貌類型 地貌類型支配著地表和土壤中水熱資源的重新分配，影響土壤生態(tài)系統(tǒng)的物質(zhì)循環(huán)過程，對(duì)土壤有機(jī)碳密度有著深遠(yuǎn)影響[22]。研究區(qū)地貌主要以丘陵為主，樣點(diǎn)在低丘、中丘和高丘皆有分布。分析結(jié)果表明(圖3)，不同地貌條件下土壤有機(jī)碳密度差異極顯著(P ＜ 0.01)，低丘、中丘和高丘的土壤有機(jī)碳密度分別為3.70、3.25、2.38 kg/m2。土壤有機(jī)碳密度隨著海拔的升高而降低，這與不同海拔條件下植被覆蓋、土壤侵蝕、耕作制度、成土物質(zhì)的再分配和景觀單元內(nèi)水分分配差異等有一定關(guān)系。土壤侵蝕使研究區(qū)內(nèi)高丘區(qū)大量肥沃表土流失，土壤肥力和植物量降低導(dǎo)致有機(jī)碳輸入量較低。低丘區(qū)多為農(nóng)用地，在人為管理下有機(jī)肥料投入量大，土壤有機(jī)碳含量較高。

圖3　不同地貌類型下土壤有機(jī)碳密度均值與方差分析Fig.3 Statistical characters of soil organic carbon densities in different terrains

2.3.5 坡度 坡度通過影響微地形的水土平衡和太陽輻射，一定程度上使不同坡度土壤有機(jī)碳密度產(chǎn)生差異[23]。將坡度用Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ級(jí)分別表示0°～5°、5°～10°、10°～15°、15°～20°、20°～25°、＞ 25°。分析結(jié)果表明，在不同的坡度條件下土壤有機(jī)碳密度存在極顯著差異(P＜0.01)。隨著坡度的增加，土壤有機(jī)碳密度呈現(xiàn)出在Ⅰ級(jí)～Ⅴ級(jí)坡度先降低，Ⅴ級(jí)～Ⅵ級(jí)坡度再升高的趨勢(shì)(圖4)，這是由于在0°～25° 坡度級(jí)別內(nèi)，研究區(qū)土壤侵蝕使坡面流失的土壤、水肥集中下坡位，使該位置土壤的有機(jī)物質(zhì)量增大，利于有機(jī)碳的積累，因而土壤有機(jī)碳密度隨著坡度的增加呈下降趨勢(shì)；Ⅴ級(jí)～Ⅵ級(jí)為坡度＞ 25° 區(qū)域，農(nóng)用地較少，多為植被覆蓋率高且人類活動(dòng)少的林地，有機(jī)碳密度相對(duì)較高。多重比較表明，Ⅰ級(jí)(0°～5°)與Ⅱ級(jí)(5°～10°)和Ⅲ級(jí)(10°～15°)的土壤有機(jī)碳密度差異極顯著(P＜0.01)，這和Ⅰ級(jí)坡度區(qū)域多水田有一定關(guān)系，表明土壤有機(jī)碳密度受坡度和土地利用類型的共同影響。

圖4　不同坡度下的土壤有機(jī)碳密度均值Fig.4 Soil organic carbon densities in different slope gradients

2.3.6 土壤質(zhì)地 土壤質(zhì)地對(duì)有機(jī)碳密度的影響與土壤礦物顆粒組合狀況關(guān)系密切[24]。該區(qū)土壤質(zhì)地主要包括砂壤、輕壤、中壤、重壤和黏土。土壤有機(jī)碳密度均值比較為黏土＞重壤＞中壤＞輕壤＞砂壤(表3)，其中有機(jī)碳密度在輕黏土中顯著高于其他4種質(zhì)地類型(P＜0.05)。這是由于黏土中黏粒礦物含量高，土壤孔隙度相對(duì)較小、通氣性差，從而緩解有機(jī)碳的分解，同時(shí)土壤碳與黏粒礦物的相互作用形成如團(tuán)聚體結(jié)構(gòu)的物理保護(hù)，使土壤有機(jī)碳不易被微生物或酶分解[25]，因此能夠較長期蓄積在土壤中；另一方面，土壤中較高的黏粒含量和粉粒含量增強(qiáng)了土壤的保水能力[26]，促進(jìn)植物生長，使更多的碳輸入到土壤中。此外，研究區(qū)內(nèi)黏土多由老沖積物及河流沖積物發(fā)育而來，用作耕地后因施用有機(jī)肥多而土壤有機(jī)碳密度高。

2.4 不同因素對(duì)土壤有機(jī)碳密度的影響程度

土壤有機(jī)碳密度受到眾多因素的共同影響，土地利用、pH、成土母質(zhì)、地貌類型、坡度和質(zhì)地都在不同程度下影響著岷江下游流域的土壤有機(jī)碳密度，但傳統(tǒng)統(tǒng)計(jì)方法難以反映各影響因素的主次性和重要性，因此本文將采用灰色關(guān)聯(lián)度法和逐步回歸分析法進(jìn)一步探討各影響因素對(duì)土壤有機(jī)碳密度的影響程度，運(yùn)用通徑分析法研究各因素間的相互重要性及其對(duì)土壤有機(jī)碳密度的直接和間接重要性。

2.4.1 灰色關(guān)聯(lián)度分析 為探討影響土壤有機(jī)碳密度的主要因素，剔除次要因素，首先采用灰色關(guān)聯(lián)度分析法進(jìn)行研究。根據(jù)公式(1～5)求得土壤有機(jī)碳密度與各影響因素的灰色關(guān)聯(lián)度，并進(jìn)行排序。從表4可以看出，在所選取的影響因素中，與土壤有機(jī)碳密度關(guān)聯(lián)最為密切的是地貌類型，其次是pH和成土母質(zhì)。地貌類型在研究區(qū)內(nèi)很大程度上影響著人類活動(dòng)、土地利用類型、pH、成土母質(zhì)和土壤質(zhì)地的差異，因其反映的環(huán)境信息量大，對(duì)土壤有機(jī)碳密度產(chǎn)生的空間分布差異解釋能力較強(qiáng)[27]。岷江下游研究區(qū)頻繁的農(nóng)業(yè)活動(dòng)受到成土母質(zhì)的限制較大，同時(shí)人類活動(dòng)對(duì)土壤pH產(chǎn)生了深刻影響，因此pH和成土母質(zhì)也是該區(qū)影響土壤有機(jī)碳密度的主要因素。

表4　土壤有機(jī)碳密度與影響因素的灰色關(guān)聯(lián)度及排序Table 4 Gray correlation degrees and sorting of soil organic carbon densities with influence factors

2.4.2 逐步回歸分析 本文主要選擇了土地利用、pH、成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地6個(gè)影響因子進(jìn)行研究，分別設(shè)為X1、X2、X3、X4、X5和X6，采用逐步回歸分析法建立逐步回歸方程(表5)。為了充分反映回歸模型的可靠性，確定建立逐步回歸方程的顯著性水平需要達(dá)到0.05。按照所設(shè)定的P＜0.05一定進(jìn)入方程，P＞0.1一定不進(jìn)入方程的原則，X2、X3和X4先后進(jìn)入方程，最終的回歸方程中常數(shù)項(xiàng)、X2、X3和X4的偏相關(guān)系數(shù)t的顯著性分別為0.000、0.000、0.003和0.006，偏相關(guān)系數(shù)是X4(-2.742)＞X3(-2.950)＞X2(-7.695)，表明X2(pH)、X3(成土母質(zhì))和X4(地貌類型)都是影響有機(jī)碳密度的重要因素，而又以地貌類型為影響土壤有機(jī)碳密度的最主要因素，成土母質(zhì)和pH次之。

表5　土壤有機(jī)碳密度與影響因素逐步回歸方程Table 5 Stepwise regression equations of soil organic carbon densities with influence factors

2.4.3 通徑分析 對(duì)各影響因素進(jìn)行逐步回歸分析后，再對(duì)pH、成土母質(zhì)和地貌類型與土壤有機(jī)碳密度進(jìn)行通徑分析。由表6可知，X4(地貌類型)是對(duì)有機(jī)碳密度直接影響中作用最大的因素，X3(成土母質(zhì))次之，X2(pH)的直接作用最小。分析各個(gè)間接通徑系數(shù)發(fā)現(xiàn)，X2(pH)通過X4(地貌類型)對(duì)有機(jī)碳密度的間接作用最大，其間接通徑系數(shù)R24×P4Y= -0.007 7。這可能是因?yàn)橥寥纏H受人類活動(dòng)的影響較大，而不同地貌條件下的土壤受到的人為干擾程度差異較明顯，因此土壤pH在不同的地貌類型區(qū)對(duì)土壤有機(jī)碳密度的影響較大[28]。

表6　影響因素對(duì)土壤有機(jī)碳密度值的通徑分析Table 6 Path analysis of relative values to soil organic carbon densities under influence factors

綜合比較灰色關(guān)聯(lián)度分析、逐步回歸分析和通徑分析 3 種方法，得出 6 個(gè)影響因素中以 pH、成土母質(zhì)和地貌類型為主。灰色關(guān)聯(lián)度的分析表明，與土壤有機(jī)碳密度關(guān)聯(lián)最為密切的是地貌類型，其次是pH 和成土母質(zhì)。逐步回歸和通徑分析的結(jié)果均表明地貌類型是影響表層土壤有機(jī)碳密度的最主要因素，成土母質(zhì)和 pH 次之。由于進(jìn)行灰色關(guān)聯(lián)度分析時(shí)僅通過比較參考序列和比較序列的曲線幾何形狀的相似程度定量關(guān)聯(lián)程度，而逐步回歸和通徑分析在考慮數(shù)據(jù)的分布特征后進(jìn)行計(jì)算，因此灰色關(guān)聯(lián)度分析與逐步回歸和通徑分析在對(duì)成土母質(zhì)和 pH 的影響程度大小上有不同的分析結(jié)果[14]，但均以地貌類型、成土母質(zhì)和 pH 為表層土壤有機(jī)碳密度的主要影響因素。

3 結(jié)論

1)本研究區(qū)位于人口稠密的岷江下游，1 138個(gè)表層(0～20 cm)土壤樣本中有機(jī)碳密度呈正態(tài)分布，其值為3.24 kg/m2，空間上呈不規(guī)則的斑塊狀分布，具有中等空間相關(guān)性，表現(xiàn)出中部淺丘區(qū)高而西北和東部高丘區(qū)低的分布趨勢(shì)。

2)該區(qū)土壤有機(jī)碳密度受到結(jié)構(gòu)性因素和隨機(jī)性因素的共同影響，土地利用類型、土壤酸堿度、成土母質(zhì)、地貌類型、坡度和土壤質(zhì)地與其關(guān)系密切，其中土地利用類型、土壤酸堿度、地貌類型、坡度和土壤質(zhì)地與有機(jī)碳密度呈極顯著相關(guān)，成土母質(zhì)對(duì)表層土壤有機(jī)碳密度的影響達(dá)顯著水平。

3)灰色關(guān)聯(lián)度的分析表明，與土壤有機(jī)碳密度關(guān)聯(lián)最為密切的是地貌類型，其次是pH和成土母質(zhì)。逐步回歸和通徑分析的結(jié)果均表明地貌類型是影響土壤有機(jī)碳密度的最主要因素，成土母質(zhì)和pH次之。綜合比較以上3種方法，得出土地利用類型、pH、成土母質(zhì)、地貌類型、坡度和質(zhì)地中以地貌類型、成土母質(zhì)和pH為影響岷江下游研究區(qū)表層土壤有機(jī)碳密度的主要因素。

參考文獻(xiàn)：

[1]Perruchoud D,Walthert L,Zimmermann S,et al.Contemporary carbon stocks of mineral forest soils in the Swiss Alps[J].Biogeochemistry,2000,50(2):111-136

[2]Puget P,Lal R.Soil organic carbon and nitrogen in a Mollisol in central Ohio as affected by tillage and land use[J].Soil and Tillage Research,2005,80(1-2):201-213

[3]Korsaeth A.Relations between nitrogen leaching and food productivity in organic and conventional cropping systems in a long-term field study[J].Agriculture,Ecosytems and Environment,2008,127(3):177-188

[4]Bajes N H.Total carbon and nitrogen in soils of the world[J].European Journal of Soil Science,2014,65(1):10-21

[5]Jobbagy E G,Jackson R B.The vertical distribution of soil organic carbon and its relation to climate and vegetation[J].Ecological Applications,2000,10(2):423-436

[6]Gau H S,Hsieh C Y,Liu C W.Application of grey correlation method to evaluate potential groundwater recharge sites[J].Stochastic Environmental Research and Risk Assessment,2006,20(6):407-421

[7]張同娟,楊勁松,劉廣明,等.基于灰色關(guān)聯(lián)度法河口地區(qū)土壤鹽分影響因子分析[J].土壤通報(bào),2010,41(4):793-796

[8]Shacham M,Brauner N.Application of stepwise regression for dynamic parameter estimation[J].Computers and Chemical Engineering,2014,69(1 016):26-38

[9]劉廣深,徐冬梅,許中堅(jiān),等.用通徑分析法研究土壤水解酶活性與土壤性質(zhì)的關(guān)系[J].土壤學(xué)報(bào),2009,40(5):756-762

[10]Zhang C S,McGrath D.Geostatistical and GIS analyses on soil organic carbon concentrations in grassland of southeastern Ireland from two different periods[J].Geoderma,2004,119(3):261-275

[11]解憲麗,孫波,周慧珍,等.不同植被下中國土壤有機(jī)碳的儲(chǔ)量與影響因子[J].土壤學(xué)報(bào),2004,41(5):687-699

[12]Selma Y K.Effects of afforestation on soil organic carbon and other soil properties[J].Catena,2014,123(10):62-69

[13]Scott N A,Tare K R,Giltrap D J.Monitoring land-use effects on soil carbon in New Zealand:Quantifying baseline soil carbon stocks[J].Environmental Pollution,2002,116(1):167-186

[14]劉思峰,黨耀國,方志耕,等.灰色系統(tǒng)理論及其應(yīng)用[M].5版.北京:科學(xué)出版社,2010:62-63

[15]劉廣深,徐冬梅,許中堅(jiān),等.用通徑分析法研究土壤水解酶活性與土壤性質(zhì)的關(guān)系[J].土壤學(xué)報(bào),2003,40(5):756-762

[16]姜春,吳志峰,錢樂祥,等.廣東山區(qū)土壤有機(jī)碳空間變異的尺度效應(yīng)[J].生態(tài)學(xué)報(bào),2013,33(16):5 118-5 125

[17]張勇,史學(xué)正,于東升,等.滇黔桂地區(qū)土壤有機(jī)碳密度變異的影響因素研究[J].土壤學(xué)報(bào),2009,46(3):526-531

[18]Nishimura S,Yonemura S,Sawamoto T,et al.Effect of land use change from paddy rice cultivation to upland crop cultivation on soil carbon budget of a cropland in Japan[J].Agriculture Ecosystems and Environment,2008,125(1-4):9-20

[19]David M,Zhang C S.Spatial distribution of soil organic carbon concentrations in grassland of Ireland[J].Applied Geochemistry,2003,18(10):1 629-1 639

[20]Tan Z X,Lal R,Smeck N E,et al.Relationships between surface soil organic carbon pool and site variables[J].Geoderma,2004,121(3-4):187-195

[21]唐曉紅,邵景安,高明,等.保護(hù)性耕作對(duì)紫色水稻土團(tuán)聚體組成和有機(jī)碳儲(chǔ)量的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2007,18(5):1 027-1 032

[22]李林海,郜二虎,夢(mèng)夢(mèng),等.黃土高原小流域不同地形下土壤有機(jī)碳分布特征[J].生態(tài)學(xué)報(bào),2013,33(1):179-187

[23]Nelson J D J,Schoenau J J,Malhi S S.Soil organic carbon changes and distribution in cultivated and restored grassland soils in Saskatchewan[J].Nutrient Cycling in Agroecosystems,2008,82(2):137-148

[24]Six J,Conant R T,Paul E A,et al.Stabilization mechanisms of soil organic matter:Implications for C-saturation of soils[J].Plant and Soil,2002,241(2):155-176

[25]Sysweerda S P,Corbin A T,Mokma D L,et al.Agricultural management and soil carbon storage in surface vs.deep layers[J].Soil Science Society of Americal Journal,2010,75(1):92-101

[26]Gruneberg E,Ziche D,Wellbrock N.Organic carbon stocks and sequestration rates of forest soils in Germany[J].Global Change Biology,2014,20(8):2 644-2 662

[27]羅由林,李啟權(quán),王昌全,等.川中丘陵縣域土壤碳氮比空間變異特征及其影響因素[J].應(yīng)用生態(tài)學(xué)報(bào),2015,26(1):177-185

[28]Xiao K C,Xu J M,Tang C X,et al.Differences in carbon and nitrogen mineralization in soils of differing initial pH induced by electrokinesis and receiving crop residue amendments[J].Soil Biology and Biochemistry,2013,67(1 016):70-84

Spatial Distribution and Its Influence Factors of Soil Organic Carbon Density of Lower Reaches of Minjiang River Based on Multivariate Analysis Method

LAI Jingjing1,2,ZHANG Shirong1*,LIU Yanting1,LI Ting2,XU Xiaoxun1,YAO Ping1,2,PU Yulin2
(1 College of Environment Sciences,Sichuan Agricultural University,Chengdu 611130,China; 2 College of Resources,Sichuan Agricultural University,Chengdu 611130,China)

Abstract:Spatial distribution and its influence factors of topsoil(0-20 cm)organic carbon density(SOCD)of the Lower Reaches of Minjiang River were studied by the combination of traditional statistical analysis,gray related analysis,stepwise regression analysis and path analysis.SOCD ranged from 0.89 to 6.49 kg/m2with a mean value of 3.24 kg/m2.Its overall distribution displayed irregular plaque-shaped,with a moderate spatial correlation,and SOCD was higher in the central low hilly than northeastern and eastern high hilly area.Among the influence factors,it showed extremely significant correlations with the land use,soil pH,topography,slope gradient and soil texture.Gray related analysis indicated that the highest correlation degree of SOCD with terrain,and followed by soil pH and soil parent material.Stepwise regression analysis combined with path analysis showed that terrain was the most significant influence factor,followed by soil parent material and pH.Therefore,terrain,soil parent material and pH were the main influence factors of topsoil organic carbon density in the lower reaches of Minjiang River.

Key words:Soil organic carbon density; Influence factor; Gray related analysis; Stepwise regression analysis; Path analysis

作者簡介：賴晶晶(1990—)，女，江西贛州人，主要從事土壤肥力評(píng)價(jià)研究。E-mail:Lcici1031@163.com

* 通訊作者(rsz01@163.com)

基金項(xiàng)目：①國家自然科學(xué)基金項(xiàng)目(41401328)和四川省科技支撐計(jì)劃項(xiàng)目(2014NZ0044)資助。

DOI:10.13758/j.cnki.tr.2016.01.024

中圖分類號(hào)：S159