羅由林，李啟權(quán)，2 *，王昌全，辛志遠(yuǎn)，王永豪，宋安軍，萬傲然，唐 科（ 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 60；2 四川農(nóng)業(yè)大學(xué)資源與地理信息技術(shù)研究所，成都 60； 重慶大學(xué)法學(xué)院，重慶 400044）

近30年川中丘陵區(qū)不同土地利用方式土壤碳氮磷生態(tài)化學(xué)計量特征變化①

羅由林1，李啟權(quán)1，2 *，王昌全1，辛志遠(yuǎn)1，王永豪1，宋安軍1，萬傲然1，唐 科3
（1 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130；2 四川農(nóng)業(yè)大學(xué)資源與地理信息技術(shù)研究所，成都 611130；3 重慶大學(xué)法學(xué)院，重慶 400044）

基于1981年第二次土壤普查數(shù)據(jù)和2012年實地采樣數(shù)據(jù)，分析了紫色丘陵區(qū)1981—2012年不同土地利用方式下表層土壤碳氮磷含量及生態(tài)化學(xué)計量特征變化。結(jié)果表明：研究區(qū)近 30年來土壤碳、氮含量增幅分別為109.98%、27.27%，磷含量基本穩(wěn)定。1981年土壤C︰N︰P比為7.28︰1︰1，2012年C︰N︰P比為16.41︰1.37︰1。各土地利用方式土壤碳、氮含量均有不同程度提升，尤以林草地、園地提升最為明顯；水田、旱地、園地土壤磷含量基本穩(wěn)定而林草地下降明顯；土壤C/N、C/P、N/P均有提升。30年來表層土壤碳儲量明顯增加，氮儲量基本穩(wěn)定而磷儲量有所下降。

土地利用方式；碳；氮；磷；生態(tài)化學(xué)計量特征

生態(tài)化學(xué)計量學(xué)是一門研究生態(tài)系統(tǒng)能量平衡和多重化學(xué)元素（主要是C、N、P）平衡等對生態(tài)交互作用的科學(xué)［1-2］。作為一種新興生態(tài)學(xué)研究領(lǐng)域，已成為生態(tài)學(xué)與生物化學(xué)、土壤化學(xué)研究領(lǐng)域的新方向，為研究土壤-植物相互作用與C、N、P循環(huán)提供了新思路［3-5］，在種群動態(tài)、養(yǎng)分循環(huán)、限制性元素判斷、養(yǎng)分利用效率乃至全球C、N、P生物地球化學(xué)循環(huán)等領(lǐng)域得到廣泛應(yīng)用［6-8］，并取得許多研究成果［9-11］。其中，C與N、P等養(yǎng)分元素的循環(huán)過程具有相互耦合關(guān)系，強(qiáng)烈的影響著生態(tài)系統(tǒng)的C、N循環(huán)過程［12］，成為全球變化的研究熱點(diǎn)之一。土壤作為陸地生態(tài)系統(tǒng)最大的有機(jī)碳庫，是植物、微生物生態(tài)系統(tǒng)C、N、P等元素的重要來源［10］。土壤系統(tǒng)中的食物鏈正日益被當(dāng)作養(yǎng)分循環(huán)的重要調(diào)節(jié)者，土壤C、N、P作為重要營養(yǎng)元素和生態(tài)環(huán)境因子，其含量變化會影響土壤中微生物數(shù)量、有機(jī)物分解速率及養(yǎng)分的長期積累［13-14］。因而研究土壤C、N、P的生態(tài)化學(xué)計量在養(yǎng)分平衡管理及C、N循環(huán)變化等方面具有重大意義。

土地利用變化和土地管理實踐措施都會改變生物地球化學(xué)和生態(tài)系統(tǒng)C、N、P生態(tài)化學(xué)計量特征［15-16］。一方面，不同土地利用方式下進(jìn)入土壤的肥料和植物殘體在數(shù)量和性質(zhì)上均有差異［17-19］；另一方面，耕作方式等農(nóng)業(yè)管理措施差異，會通過影響土壤養(yǎng)分的礦化、運(yùn)輸、吸收利用等差異進(jìn)而影響土壤C、N、P等養(yǎng)分循環(huán)［20-23］。隨著氣候變化和生態(tài)環(huán)境的惡化［24-25］，農(nóng)業(yè)生產(chǎn)理念已發(fā)生重大改變，由高投入高產(chǎn)出的思想逐漸轉(zhuǎn)變?yōu)楦又匾曓r(nóng)業(yè)生產(chǎn)與生態(tài)環(huán)境的協(xié)調(diào)，實現(xiàn)經(jīng)濟(jì)效益、社會效益與環(huán)境效益的統(tǒng)一［26-27］。這對農(nóng)業(yè)施肥管理和土地利用方式優(yōu)化提出了更高的要求。因而研究不同土地利用方式對土壤C、N、P含量及其對生態(tài)系統(tǒng)的長期影響對現(xiàn)代農(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境保護(hù)具有迫切而重要的現(xiàn)實意義。

近幾十年來，中國土地利用變化劇烈，生態(tài)環(huán)境狀況面臨嚴(yán)峻形勢，政府對土地利用結(jié)構(gòu)和方式的進(jìn)行了不斷調(diào)整，如在長江中上游的退耕還林工程、京津生態(tài)防護(hù)林建設(shè)、西部大開發(fā)等策略，這些勢必會對土壤C、N、P及生態(tài)系統(tǒng)造成影響［28-30］。四川盆地屬于西部大開發(fā)地區(qū)，也是中國南方典型的紫色盆地，盆地丘陵區(qū)人口密度大，土地墾殖率高，是重要的農(nóng)耕區(qū)和長江水土流失重點(diǎn)防治區(qū)。該區(qū)已有研究多集中于土壤養(yǎng)分空間變異特征方面［31-32］，且多涉及為小流域或小尺度地區(qū)［33-34］，有關(guān)其土壤C、N、P長期變化及其生態(tài)化學(xué)計量特征的資料缺乏；而研究其土壤C、N、P及其生態(tài)化學(xué)計量特征變化趨勢對該區(qū)土壤C、N、P調(diào)控及生態(tài)效應(yīng)改善皆有重要指導(dǎo)作用。本研究選取川中丘陵第一大縣仁壽縣為案例區(qū)，通過整合該區(qū)1981年的214個第二次土壤普查數(shù)據(jù)并結(jié)合2012年實地采集的555個土樣數(shù)據(jù)比較分析近30年來該區(qū)不同土地利用方式下土壤C、N、P變化及其生態(tài)化學(xué)計量特征，以期為紫色丘陵區(qū)養(yǎng)分動態(tài)平衡管理與生態(tài)環(huán)境的協(xié)調(diào)統(tǒng)一提供科學(xué)依據(jù)。

1 材料和方法

1.1 研究區(qū)概況

研究區(qū)位于四川盆地丘陵區(qū)（29°37′～30°20′N， 103°54′～104°29′E），面積2 606 km2，其中農(nóng)業(yè)用地約占82%，人口162萬，人口密度約為600人/km2，轄22個鎮(zhèn)、38個鄉(xiāng)，為川中丘陵區(qū)第一大縣和第一人口大縣。地貌以丘陵為主，海拔位于345～972 m之間，含部分低山區(qū)域（圖1）?？h域處于亞熱帶濕潤性季風(fēng)氣候區(qū)，四季分明，雨熱同期；年均氣溫17.4 ℃，年均降雨1 009.4 mm，年均日照時數(shù)1 196.6 h，無霜期312天。土壤母質(zhì)除有少量第四系老沖洪積物外，主要為侏羅系、白堊系和三疊系出露的紫色巖層；主要發(fā)育土壤類型為紫色土、水稻土、黃壤。土地利用方式包括耕地、園地、林草地等，土地利用程度較高。20世紀(jì) 80年代全區(qū)森林覆蓋率約為 4.70%，到 21世紀(jì)初恢復(fù)達(dá)到 19.96%，是長江上游人口最為密集的典型農(nóng)業(yè)區(qū)和水土流失重點(diǎn)防治區(qū)。

圖1 研究區(qū)位置、采樣點(diǎn)分布圖及土地利用現(xiàn)狀Fig. 1 Basic geographic information of study area， spatial distribution of soil sampling sites， and map of land use types

1.2 數(shù)據(jù)來源

研究區(qū)1981年土壤數(shù)據(jù)來自于第二次全國土壤普查時采集的典型土壤剖面。這些典型土壤剖面數(shù)據(jù)資料記錄在《仁壽土壤》、《四川土壤》、《區(qū)域土壤地理》等資料中，每一個剖面都較詳細(xì)地記錄了該土壤剖面的采樣地點(diǎn)、成土母質(zhì)、土壤類型、地形特征、土地利用方式以及土壤理化性質(zhì)等信息。本研究通過收集整理這些資料，建立了該區(qū)第二次全國土壤普查的典型土壤剖面數(shù)據(jù)庫，并結(jié)合行政區(qū)劃、居民點(diǎn)空間分布數(shù)據(jù)、土壤類型圖、母質(zhì)圖、植被類型圖、30 m分辨率DEM和1︰1萬地形圖等輔助資料對各剖面點(diǎn)位置逐一進(jìn)行了定位［35］，盡可能使各土壤剖面的位置接近原采樣點(diǎn)位置，同時使定位處的環(huán)境條件符合各土壤剖面記錄的環(huán)境信息；共整理獲得214個具有相對準(zhǔn)確位置信息的典型剖面（圖1）。

1.3 樣品采集

2012年數(shù)據(jù)在根據(jù)研究區(qū)的實際情況并考慮樣點(diǎn)代表性與空間分布均勻性的基礎(chǔ)上，結(jié)合研究區(qū)成土母質(zhì)、土壤類型和土地利用方式等資料進(jìn)行采樣點(diǎn)布設(shè)，在全縣共計采集表層（0 ～ 20 cm）土壤樣點(diǎn)558個（圖1）。每一采樣點(diǎn)周圍5 m范圍內(nèi)取4個點(diǎn)，混合后以四分法取樣。采樣的同時記錄每個采樣點(diǎn)的地理坐標(biāo)、海拔、成土母質(zhì)、土地利用方式、土壤類型、種植制度和施肥狀況等地表環(huán)境信息。土壤樣品在室內(nèi)自然風(fēng)干后，去除石礫和根系等雜物，研磨過100目篩待測。

1.4 實驗方法

土壤有機(jī)碳（SOC）含量采用重鉻酸鉀外源加熱法測定，總磷（TP）含量采用鉬銻抗比色法測定［36］；總氮（TN）含量采用Foss Kjeltec8400全自動凱氏定氮儀進(jìn)行測定。

1.5 數(shù)據(jù)處理

利用均值 ± 3 倍標(biāo)準(zhǔn)差的方法剔除了 2012 年的 3 個異常樣點(diǎn)數(shù)據(jù)，用 1981 年的 214 個樣點(diǎn)數(shù)據(jù)和 2012 年剩余的 555 個樣點(diǎn)數(shù)據(jù)進(jìn)行分析。文中圖表制作在 Excel 2010 實現(xiàn)，基本統(tǒng)計特征、方差分析等均在 SPSS 20 完成。

1.6 土壤碳氮磷儲量計算方法

本文土壤C、N、P儲量估算方法如下：

式中：S為土壤有機(jī)碳（全氮、全磷）儲量（單位：t/hm2）；C為有機(jī)碳（全氮、全磷）含量（單位：g/kg）；ρ為土壤容重（單位：g/cm2）；D為土層厚度（單位：cm）；A為面積（單位：km2）。

2 結(jié)果與分析

2.1 土壤碳氮磷含量及其相關(guān)關(guān)系

統(tǒng)計結(jié)果顯示（表1），30年來該區(qū)C含量（即有機(jī)碳）由6.41 g/kg上升到13.46 g/kg，增幅109.98%；N含量（即全氮）由0.88 g/kg上升到1.12 g/kg，增幅27.27%；P含量（即全磷）在1981年為0.89 g/kg，到2012年為0.82 g/kg，基本穩(wěn)定。C、N、P含量的不同變化特征，引起C︰N︰P三者間比例關(guān)系發(fā)生較大變化，1981年C︰N︰P比約為7.28︰1︰1，2012 年C︰N︰P比為16.41︰1.37︰1。其中土壤C/N由6.85上升到 11.93，增幅 75.96%；相關(guān)系數(shù)由 0.71降至0.68。N/P由1.15升到1.69，增幅46.96%。N/P 由8.25上升到21.67，增幅162.67%。全區(qū)土壤C/N、C/P、N/P之間的比例關(guān)系均有提升，而相關(guān)系數(shù)均有不同程度的下降（圖2）。

表1 不同時期土壤C、N、P含量及比值Table1 Descriptive statistics of contents and ratios of soil carbon， nitrogen and phosphorous contents in different years

圖2 1981年和2012年研究區(qū)土壤C、N、P相關(guān)關(guān)系Fig. 2 Correlations of soil carbon， nitrogen and phosphorus in 1981 and 2012

2.2 不同土地利用方式土壤碳氮磷含量變化

不同土地利用方式下的土壤C、N、P呈現(xiàn)出不同的變化特征（表2）。1981年各土地利用方式間C含量差異明顯，表現(xiàn)為水田＞林草地＞園地、旱地；2012年表現(xiàn)為水田、林草地＞園地＞旱地；30年來，各土地利用方式土壤 C含量均大幅增加，水田與林草地之間C含量已無明顯差異，園地與旱地之間C含量已有顯著差異。N含量也有顯著變化，1981年N含量為水田＞林草地、旱地、園地；2012年表現(xiàn)為水田＞林草地＞園地＞旱地，30年來林草地、園地N含量明顯變化，與C呈現(xiàn)出類似的變化趨勢。P含量變化與C、N含量變化不同，1981年僅園地與林草地之間P含量存在顯著差異。2012年P(guān)含量為：園地、旱地＞水田＞林草地。2個時期P含量均以園地最高，林草地最低，30年來旱地、園地、水田P含量基本穩(wěn)定，林草地顯著下降。

總體來看，30年來不同土地利用方式C含量均有顯著增加；水田、園地、林草地 N含量顯著增加而旱地增長不明顯；P含量在林草地明顯下降，在水田、園地、旱地均無顯著變化。

表2 不同土地利用方式土壤C、N、P含量統(tǒng)計特征Table2 Descriptive statistics of soil carbon， nitrogen and phosphorous contents under different land use types

2.3 土壤碳氮磷化學(xué)計量學(xué)特征

土壤C、N、P生態(tài)化學(xué)計量學(xué)能揭示生態(tài)系統(tǒng)中各組分養(yǎng)分比例調(diào)控機(jī)制，認(rèn)識養(yǎng)分比例在生態(tài)系統(tǒng)過程和功能中的作用，利于揭示元素間相互作用與制約規(guī)律。近30年來，不同土地利用方式C、N、P含量的變化顯著影響著土壤 C︰N︰P比值關(guān)系（表3）。1981年各土地利用方式下土壤C/N表現(xiàn)為：水田、林草地＞旱地、園地；水田與林草地、旱地與園地C/N接近，但水田、林草地土壤C/N顯著高于旱地、園地。C/P表現(xiàn)為水田＞林草地、旱地、園地，N/P表現(xiàn)為水田＞林草地、旱地、園地；水田 C/N、C/P均顯著高于林草地、旱地、園地，而旱地、園地、林草地C/P之間無顯著差異。2012年不同土地利用方式C/N、C/P、N/P發(fā)生變化。其中C/N為：林草地＞水田＞園地、旱地。C/P為：林草地、水田＞園地、旱地。N/P為林草地、水田＞園地、旱地。土壤C/N、C/P、N/P在1981年均以水田最高，園地最低，在2012年以林草地最高而旱地最低。

不同土地利用方式土壤 C/N、C/P、N/P變化幅度呈現(xiàn)不同變化特征。30年來水田和林草地 C/N、C/P、N/P均大幅提升。旱地 C/N、C/P大幅提升，N/P基本保持穩(wěn)定。園地C︰N︰P比值均有提升。各種土地利用方式下C/P均大幅增長，其中林草地、園地增長約2倍，水田、旱地增長約1倍。水田C/N、N/P增幅均較?。缓档谻/N增幅均較大而N/P基本保持穩(wěn)定。這與不同土地利用方式土壤C、N、P含量變化幅度的差異息息相關(guān)。

表3 不同土地利用方式下土壤C、N、P關(guān)系Table3 Correlations of soil carbon， nitrogen and phosphorous contents under different land use types

2.4 土壤碳氮磷儲量的變化

據(jù)《仁壽土壤》記載，1981年該區(qū)土壤容重平均值為1.44 g/m3，根據(jù)該區(qū)土壤有機(jī)質(zhì)含量與土壤容重間的對應(yīng)關(guān)系，2012年土壤容重平均值為1.13 g/m3。全區(qū)農(nóng)業(yè)用地面積約2 135 km2；估算得到其土壤C、N、P儲量信息（表4）：1981年該區(qū)農(nóng)業(yè)土壤C儲量約3.94×106t，2012年約6.96×106t，30年間儲量增幅64.74%。N儲量在1981年約5.41×105t，2012年約為5.40×105t，儲量相對穩(wěn)定。P儲量在1981年約為5.41×105t，2012年約為3.96×105t，降幅26.80%。 30年來，表層土壤C儲量上升，N儲量相對穩(wěn)定，P儲量下降，表層土壤C、N、P含量呈現(xiàn)出不平衡變化；C、N、P儲量變化幅度明顯低于C、N、P含量變化幅度。

表4 不同時期土壤C、N、P儲量Table4 Carbon， nitrogen and phosphorous storages in different years

3 討論

3.1 對土壤碳氮磷含量的影響

近30年來，不同土地利用方式下的土壤C、N、P呈現(xiàn)出不同的變化特征。各土地利用方式土壤C含量均大幅增加，并引起不同土地利用方式之間 C含量差異改變；其中林草地與水田之間 C含量已無顯著差異，園地 C含量顯著高于旱地，反映出林草地與園地土壤C含量明顯上升，C儲能力顯著。這主要是自然土壤受人為翻耕影響小，C匯作用較為明顯［37］。林草地、園地N含量明顯變化，與C呈現(xiàn)出類似的變化趨勢，顯示出園地、林草地C、N之間有著較好相關(guān)性。這主要是隨著時間推移，植被生產(chǎn)力的恢復(fù)和提升，大量地上凋落物歸還土壤，有機(jī)碳易于累積增長，而土壤中 N素多以有機(jī)氮的形態(tài)存在于有機(jī)質(zhì)中，林草地C、N之間又具有良好的相關(guān)性，進(jìn)而N含量也會隨C含量增加而增加［38］。

P含量變化與C、N含量變化不同，表現(xiàn)為旱地、園地、水田P含量基本穩(wěn)定，主要是旱地、園地、水田土壤有一定P肥補(bǔ)充且P作為沉積性礦物，較為穩(wěn)定；林草地顯著下降。這主要與該區(qū)林草地多源于退耕還林有關(guān)和植物生長年限有關(guān)，有研究顯示隨著植被生長年限的增加表層土壤P含量會呈現(xiàn)出現(xiàn)上升后下降的趨勢［34，39］。總體來看，30年來不同土地利用方式土壤C含量均有顯著增加；水田、園地、林草地 N含量顯著增加而旱地增長不明顯；P含量在水田、園地、旱地均無顯著變化，林草地 P含量明顯下降。

3.2 對土壤碳氮磷化學(xué)計量的影響

土壤C、N、P三者之間存在著相互平衡與制約的關(guān)系，土壤C︰N︰P化學(xué)計量比值差異能反映并影響這種關(guān)系；土壤 C/N主要通過影響微生物活性來影響有機(jī)質(zhì)分解速率，進(jìn)而影響土壤C、N循環(huán)；N/P可作為N飽和的診斷指標(biāo)，并被用于確定養(yǎng)分的闕值。在這項研究中，研究區(qū)近30年來，不同土地利用方式土壤C︰N︰P比值呈現(xiàn)出不同變化特征（表3）。其中，水田C/P、N/P增加，水田C/N、C/P、N/P相關(guān)系數(shù)（R2）均明顯下降，表明水田在30年來N素施用量明顯增加，而有機(jī)碳增長速度相對減緩，繼續(xù)增施N肥對提升C儲量意義甚微，其土壤N素持有水平基本已經(jīng)飽和，并可能存在施N過量的現(xiàn)象［40］；同時長期以來雖存在著偏施N肥而相對忽視P肥投入的現(xiàn)象，但P素仍舊保持穩(wěn)定水平，這主要與西部紫色丘陵區(qū)土壤P背景值較高有關(guān)。旱地C/N、C/P增加，N/P基本保持穩(wěn)定。旱地有機(jī)碳提升明顯，全N、全P含量基本處于穩(wěn)定狀態(tài)，30年來旱地土壤肥力狀況有所改善。

園地 C/N、C/P、N/P均有不同程度的增加，C︰N︰P比例關(guān)系明顯改善。這與園地C、N含量明顯增加而P含量基本穩(wěn)定有關(guān)。園地土壤C︰N相關(guān)系數(shù)有較大提升，反映出園地土壤C、N關(guān)系呈良性發(fā)展態(tài)勢。林草地 C/N、C/P、N/P均有大幅提升，其中C/P增長高達(dá)2倍，這主要與林草地C、N含量大幅提升而P含量明顯下降和植被凋落物N/P變化有關(guān)；已有長期觀測研究顯示，隨著植物生長，新鮮凋落物和腐殖質(zhì)的 N/P將增加［34］，凋落物進(jìn)入土壤引起林草地土壤N/P相應(yīng)上升。這反映出林草地近30年來固C能力穩(wěn)步提升，植被對土壤P的消耗量高于其生產(chǎn)補(bǔ)給量，土壤P含量制約著整個林草地C︰N︰P關(guān)系。各種土地利用方式下C/P值均大幅增長，其中林草地、園地增幅約2倍，水田、旱地增幅約1倍。園地C/N增幅最大，林草地N/P增幅最高，水田 C/N、N/P增幅均較小；旱地 C/N增幅均較大而N/P基本保持穩(wěn)定。這與各土地利用方式土壤C、N、P含量變化幅度密切相關(guān)，其C、N、P含量變化影響著C︰N︰P三者之間的比值關(guān)系，同時C︰N︰P三者之間的比值關(guān)系通過影響其養(yǎng)分有效性進(jìn)而作用于土壤C、N、P含量變化。

4 結(jié)論

研究區(qū)近30年來C、N含量均有不同程度增長，P含量相對穩(wěn)定。C/N、C/P、N/P也表現(xiàn)出不同程度的上升，尤以N/P增至1981年的3倍，但其相關(guān)系數(shù)總體均呈現(xiàn)出下降趨勢。不同土地利用方式下土壤C、N、P含量及C︰N︰P比值關(guān)系呈現(xiàn)出不同的變化特征。C、N含量均以林草地、園地增幅最為明顯，P含量除林草地明顯下降外，水田、旱地、園地均基本穩(wěn)定。各土地利用方式土壤 C/N、C/P、N/P均有不同程度上升，旱地土壤C、N、P關(guān)系變化相對穩(wěn)定而園地C、N、P之間的關(guān)系顯著改善；林草地土壤C、N儲量明顯提升。近30年來該區(qū)表層土壤C儲量明顯增加，N儲量基本穩(wěn)定，P儲量有所下降。

［1］ Stemer R W， Elser J J. Ecological stoichiometry： The biology of elements from molecules to the biosphere［M］. Princeton： Princeton University Press， 2002

［2］ 王紹強(qiáng)， 余貴瑞. 生態(tài)系統(tǒng)碳氮磷元素的生態(tài)化學(xué)計量學(xué)特征［J］. 生態(tài)學(xué)報， 2008， 28（8）： 3 937-3 947

［3］ Michaels A F. The ratios of life［J］. Science， 2003， 300：906-907

［4］ Reich P B， Oleksyn J. Global patterns of plant leaf N and P in relation to temperature and latitude［J］. Proceedings of the National Academy of Sciences of the United States of America， 2004，101 （30）：11 001-11 006

［5］ Elser J J， Bracken M E， Cleland E E， et al. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater， marine and terrestrial ecosystems［J］. Ecology Letters， 2007， 10：1 135-1 142

［6］ Hobbie S E， Nadelhoffer K， Hgbery P. A synthesis： The role of nutrients as constraints on carbon balance in boreal and arctic regions［J］. Plant and Soil， 2002， 242：163-170

［7］ 王加恩， 鄭潔， 康占軍， 等. 浙北典型稻作區(qū)近 30年表層土壤碳氮含量變化［J］. 生態(tài)環(huán)境學(xué)報， 2010， 19（5）：1 191-1 196

［8］ Yang Y H， Luo Y Q. Carbon： nitrogen stoichiometry in forest ecosystems during stand development［J］. Global Ecology and Biogeography， 2011， 20（2）： 354-361

［9］ 曾全超， 李鑫， 董揚(yáng)紅， 等. 陜北黃土高原土壤性質(zhì)及生態(tài)化學(xué)計量的緯度變化特征［J］. 自然資源學(xué)報， 2015，30（5）： 870-879

［10］ Gao Y， He N P， Yu G R， et al. Long-term effects of different land use types on C， N， and P stoichiometry and storage in subtropical ecosystems： A case study in China［J］. Ecological Engineering， 2014， 67： 171-181

［11］ Deng Q， Cheng X L， Yang Y H， et al. Carbon-nitrogen interactions during afforestation in central China［J］. Soil Biology and Biochemistry， 2014， 69： 119-122

［12］ Cleveland C C， Liptzin D. C︰N︰P stoichiometry in soil：Is there a “Redfield ratios” for the microbial biomass？［J］. Biogeochemistry， 2007， 85： 235-252

［13］ Tian H Q， Chen G S， Zhang C， et al. Pattern and variation of C︰N︰P ratios in Chin's soils： a synthesis of observational date［J］. Biogeochemistry， 2010， 98（1/3）：139-151

［14］ Kirkby C A， Kirkegaard J A， Richardson A E， et al. Stable soil organic matter： A comparison of C︰N︰P︰S ratios in Australian and other world soils［J］. Geoderma， 2011，163（3/4）： 197-208

［15］ Power J S. Changes in soil carbon and nitrogen after contrasting land-use transitions in Northeastern Costa Rica［J］. Ecosystems， 2004， 7（2）： 134-146

［16］ Wei X R， Shao M G， Fu X L， et al. Distribution of soil organic C， N and P in the adjacent land use patterns in the northern Loess Plateau， China［J］. Biogeochemistry， 2009，96（1）：149-162

［17］ 侯曉靜， 楊勁松， 趙曼， 等. 不同施肥措施對濱海鹽漬土有機(jī)碳含量的影響［J］. 土壤， 2014， 46（5）： 780-786

［18］ 楊文， 周腳根， 王美慧， 等. 亞熱帶丘陵小流域土壤碳氮磷生態(tài)計量特征的空間分異性［J］. 土壤學(xué)報， 2015，52（6）： 1 336-1 344

［19］ 青燁， 孫飛達(dá)， 李勇， 等. 若爾蓋高寒退化濕地土壤碳氮磷比及相關(guān)性分析［J］. 草業(yè)學(xué)報， 2015， 24（3）： 38-47

［20］ 李小剛， 馬有華， 王強(qiáng)， 等. 縣域耕地土壤養(yǎng)分的時空變化及區(qū)域施肥研究-以肥東縣為例［J］.土壤， 2014， 46（6）：976-983

［21］ 周睿， 潘賢章， 王昌坤， 等.上海市城郊土壤有機(jī)質(zhì)的時空變異特征及其影響因素［J］. 土壤， 2014， 46（3）： 433-438

［22］ Alvarez C ， Alvarez C R， Costantini A， et al. Carbon and nitrogen sequestration in soils under different management in the semi-arid Pampa （Argentina）［J］. Soil and Tillage Research， 2014， 142： 25-31

［23］ 候曉靜， 楊勁松， 王相平， 等. 不同施肥方式下灘涂圍墾農(nóng)田土壤有機(jī)碳及團(tuán)聚體有機(jī)碳的分布［J］. 土壤學(xué)報，2015， 52（4）： 818-827

［24］ Lutze J L， Gifford R M. Nitrogen accumulation and distribution in Danthonia richardsonii swards in response to CO2and nitrogen supply over four years of growth［J］. Global Change Biology， 2000， 6（1）： 1-12

［25］ Wei R， Han Q T， Bo T， et al. China's crop productivity and soil carbon storage as influenced by multifactor global change［J］. Global Change Biology， 2012， 18： 2 945-2 957

［26］ Gelaw A M， Singh B R， Lal R. Soil organic carbon and total nitrogen stocks under different land uses in a semi-arid watershed in Tigray， Northern Ethiopia. Agriculture， Ecosystems and Environment， 2014， 188：256-263

［27］ Bayer C， Gomes J， Zanatta J A， et al. Soil nitrous oxide emissions as affected by long-term tillage， cropping systems and nitrogen fertilization in Southern Brazil［J］. Soil and Tillage Research， 2015， 146： 213-222

［28］ Zu B X， Jiang G Z， Gang L， et al. Soil organic carbon stocks in China and changes from 1980s to 2000s［J］. Global Change Biology， 2007， 13： 1 989-2 007

［29］ Wang H J， Shi X Z， Yu D S， et al. Factors determining soil nutrient distribution in a small-scaled watershed in the purple soil region of Sichuan Province， China［J］. Soil and Tillage Research， 2009， 105（2）： 300-306

［30］ Zeng X H， Zhang W J， Cao J S， et al. Changes in soil organic carbon， nitrogen， phosphorus， and bulk density after afforestation of the “Beijing-Tianjin Sandstorm Source Control” program in China［J］. Catena， 2014， 118： 186-194

［31］ 羅由林， 李啟權(quán)， 王昌全， 等. 川中丘陵縣域土壤氮素空間分布特征及其影響因素［J］. 環(huán)境科學(xué)， 2015， 36（2）：652-660

［32］ 羅由林， 李啟權(quán)， 王昌全， 等. 四川省仁壽縣土壤有機(jī)碳空間分布特征及其主控因素［J］. 中國生態(tài)農(nóng)業(yè)學(xué)報，2015， 23（1）： 34-42

［33］ 李婷， 張世熔， 劉潯， 等. 沱江流域中游土壤有機(jī)質(zhì)的空間變異特點(diǎn)及其影響因素［J］. 土壤學(xué)報， 2011， 48（4）：863-868

［34］ Wardle D A， Walker R， Bardgett R D. Ecosystem properties and forest decline in contracting long-term chronosequences［J］. Science， 2004， 305： 509-513

［35］ Li Q Q， Yue T X， Wang C Q， et al. Spatially distributed modeling of soil organic matter across China： An application of artificial neural network approach［J］. Catena，2013， 104： 210-218

［36］ 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法［M］. 北京： 中國農(nóng)業(yè)科技出版社， 2000

［37］ 孔濤， 張德勝， 寇涌蘋， 等. 渾河上游典型植被河岸帶土壤有機(jī)碳、全氮和全磷分布特征［J］. 土壤， 2014， 46（5）：793-798

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

［39］ 李煒， 鄭子成， 李廷軒. 不同植茶年限土壤團(tuán)聚體碳氮磷生態(tài)化學(xué)計量學(xué)特征［J］. 應(yīng)用生態(tài)學(xué)報， 2015， 26（1）：9-16

［40］ 羅由林， 李啟權(quán)， 王昌全， 等. 近 30年來川中紫色丘陵區(qū)土壤碳氮時空演變格局及其驅(qū)動因素［J］. 土壤學(xué)報，2016， 53（3）： 582-593

Last 30a Changes of C， N and P Ecological Stoichiometry of Different Land Use Types in Hilly Area of Mid-Sichuan Basin， Southwest China

LUO Youlin1， LI Qiquan1，2 *， WANG Changquan1， XIN Zhiyuan1， WANG Yonghao1，SONG Anjun1， WAN Aoran1， TANG Ke3
（1 College of Resources， Sichuan Agricultural University， Chengdu 611130， China； 2 Institute of Resources and Geography Information Technology， Sichuan Agricultural University， Chengdu 611130， China； 3 Law School， Chongqing University，Chongqing 400044， China）

Based on the 2ndNational Soil Survey Data in 1980s and the measured data of 555 soil samples （0-20 cm） in 2012， the changes of C， N and P ecological stoichiometry of different land use types over 30 years were analyzed in this paper. The results showed that soil carbon and nitrogen increased by 109.98% and 27.27% respectively， while phosphorus content remained stable in the last 30 years. Soil C︰N︰P ratio of 1981 was 7.28︰1︰1 and was 16.41︰1.37︰1 in 2012. The contents of soil carbon and nitrogen all increased in each land use type， especially in forest and grassland and garden； the contents of soil phosphorus retained stable in paddy field， dry land and orchard while decreased significantly in forest and grassland. Soil C/N，C/P and N/P increased in each land use type. In the last 30 years， carbon storage of topsoil increased significantly， nitrogen storage remained stable and phosphorus storage declined.

Land use type； Carbon； Nitrogen； Phosphorous； Ecological stoichiometry

S152.4；S153.6

10.13758/j.cnki.tr.2016.04.016

國家自然科學(xué)基金項目（4120124）和四川省科技支撐計劃項目（四川省-中科院合作項目2012JZ0003）資助。

（liqq@lreis.ac.cn）

羅由林（1990—），男，四川樂山人，博士研究生，主要研究方向為土壤資源環(huán)境。E-mail： S20132717@163.com