黃梅 周運超 白云星

摘 要：為研究喀斯特系統(tǒng)土壤DOC淋失在不同土壤厚度下對氮沉降的響應過程。以喀斯特石灰土為供試體，采用土柱淋溶法，試驗設計3種不同土壤厚度（10、25、50 cm）的淋溶管，淋溶管為 PVC 材質，直徑 20 cm，選擇3種不同氮沉降LN（191.08 g/m2）、MN（254.78 g/m2）、HN（318.47 g/m2），土表施予多年生黑麥草為研究對象。結果表明：氮沉降對土壤DOC淋失有一定的促進作用;土壤厚度對土壤DOC淋失有顯著影響，土壤厚度越淺，氮沉降越大，DOC淋失越大，反之則越小。該研究可為解決喀斯特系統(tǒng)土壤碳循環(huán)問題與土壤質量評價提供一定參考和依據(jù)。

關鍵詞：氮沉降;土壤厚度;喀斯特;土壤DOC;淋失

中圖分類號：S152;Q948 ? ?文獻標識碼：A ? 文章編號：1006-8023（2019）05-0043-07

Abstract：In order to study how the response of soil DOC leaching to nitrogen deposition under different soil thickness in karst system， using karst lime soil as the test sample， three different soil thickness （10 cm， 25 cm， 50 cm） leaching tubes were designed by soil column leaching method. The leaching tube was made of PVC material with a diameter of 20cm. Three different nitrogen depositions， LN （191.08 g/m2）， MN （254.78 g/m2）， and HN （318.47 g/m2） were selected， and perennial ryegrass was applied to the soil surface. The results showed that nitrogen deposition had a certain effect on soil DOC leaching. Soil thickness had a significant effect on soil DOC leaching. The shallower the soil thickness， the larger the nitrogen deposition， the larger the DOC leaching loss， and vice versa. This study can provide some reference and basis for solving the soil carbon cycle problem and soil quality evaluation of karst system.

Keywords：Nitrogen deposition; soil thickness; karst; soil DOC; leaching

0 引言

隨著社會發(fā)展，工業(yè)生產(chǎn)進程的加快，施肥和化學染料燃燒等人類活動導致在全球范圍內大氣氮沉降的現(xiàn)象日益增加。目前，人為氮沉降的污染主要集中在北美、歐美及亞洲部分地區(qū)[1]，我國已經(jīng)成為僅次于歐洲和美國的第三大氮沉降區(qū)域[2]，且存在氮沉降持續(xù)升高的巨大風險[3]，對生態(tài)系統(tǒng)平衡產(chǎn)生極大的負面影響。氮素的輸入是影響陸地生態(tài)系統(tǒng)碳源/匯的四大機制之一[4]，而土壤作為氮沉降的受體，外氮源的輸入會使土壤有機碳礦化[5]，直接或間接影響土壤碳的輸入輸出[6]，并對土壤碳循環(huán)的動態(tài)發(fā)展不利[7-8]。

土壤溶解性有機碳（DOC）有一定溶解性，且在土壤中移動較快、不穩(wěn)定、易氧化、易分解和易礦化[9]，是土壤有機碳的組成部分之一。作為生態(tài)系統(tǒng)碳循環(huán)的基本組成部分，DOC淋失是土壤有機碳損失的重要途徑[9]，DOC淋失變化對土壤碳循環(huán)有著重要的影響。碳、氮作為陸地生態(tài)系統(tǒng)最重要的兩大生源要素[10]，已有部分學者研究了氮沉降與土壤DOC的關系，Deforest等[11]，Gundersen等[12]，許凱等[13]，魏春蘭等[14]都研究了氮沉降對森林土壤中DOC的變化，也有部分學者研究了不同的施肥方式[15]和生物炭添加[16]和地表徑流[17]等與土壤DOC淋失量的影響關系，但鮮見氮沉降對土壤DOC淋失量的影響。因石漠化問題的加劇，使得喀斯特生態(tài)系統(tǒng)極為敏感和脆弱[18]，喀斯特地區(qū)地形破碎、土壤分布不均勻、土層淺薄、土地利用類型復雜多樣[18]，且水熱狀況、植被的立地條件以及土壤的發(fā)育程度等都與非喀斯特地區(qū)不同[19]?？λ固厥攸c使得氮沉降變化下，研究土壤碳循環(huán)與土壤質量評價等問題帶來了的巨大挑戰(zhàn)。故從喀斯特系統(tǒng)背景出發(fā)，研究喀斯特土壤厚度下土壤DOC淋失對氮沉降的響應，對喀斯特土壤碳庫儲量的估算與土壤質量評價具有重要意義。為此本文以喀斯特石灰土為供試體，采用土柱淋溶法，試圖明晰喀斯特系統(tǒng)土壤DOC淋失在不同土壤厚度下如何響應氮沉降的變化，以期為喀斯特系統(tǒng)土壤碳循環(huán)問題與土壤質量評價提供一定參考和依據(jù)。

1 材料與方法

1.1 試驗材料

1.1.1 試驗地概況

試驗地點設在貴州省貴陽市貴州大學南校區(qū)林學院苗圃（E106°27′-106°52′，N26°11′-26°34′），其位于貴州高原中部、貴陽市南郊。該區(qū)海拔1 100 m，屬中國亞特帶濕潤溫和型氣候類型，年平均降水量1 178.3 mm，最高的1977 年為1 539.2 mm，最低的1962年為891.6 mm;年日照時數(shù)在1 125～1 140 h之間，年日平均氣溫為15.3 ℃，最熱月（7月）和最冷月（1月）的平均氣溫分別為23.7、4.6 ℃。該區(qū)供試的土壤類型為喀斯特發(fā)育形成的石灰土，土壤養(yǎng)分含量見表1。

1.1.2 淋溶試驗裝置

采用模擬野外巖溶土壤系統(tǒng)的土柱淋溶法，土柱裝置為PVC材質的淋溶管，高分別為10、25、50 cm，內徑20 cm，在圓柱管下端鋪上2 cm厚的用去離子水洗凈的干燥石英砂（起過濾水樣作用）并在底部管口處鋪上尼龍網(wǎng)，將圓柱底部封緊。按容重為1.2 g/cm3混合均勻的石灰土分別裝入10、25、50 cm的淋溶管內，將各土柱分別裝入10、25、50 cm土層，安裝土柱時，將土柱壁邊緣的土壤壓實，確保無貼壁水流入滲，盡量減少邊緣效應。

1.2 試驗方法

1.2.1 試驗設計

于2016年5月，在貴州大學南校區(qū)苗圃選擇9塊5 m×5 m立地條件及長勢基本一致的試驗地，土柱安裝好之后，將之前稱好的重量為5 g的多年生黑麥草種子均勻撒在距離土壤表面只有2～3 cm的土壤上，表面再用2～3 cm的土壤覆蓋。依據(jù)隨機區(qū)組試驗設計原理進行樣地設置，本試驗共分3個處理水平，氮沉降處理水平的設置主要基于目前自然氮沉降水平及其至本世紀末的可能變化，參考國內已有施氮試驗[20-21]。按氮沉降量由低到高分別標記為低氮（LN，191.08 g/m2）、中氮（MN，254.78 g/m2）、高氮（HN，318.47 g/m2），每個處理重復3次，施用氮肥均為尿素（含N 460 g/kg），共計 27 個處理，具體施用方法是： 將需要施入的氮量換算成每個小區(qū)所需的尿素量，將尿素充分溶于 20 L 水中，在模擬自然降雨前，以噴壺在草地人工來回噴灑。

1.2.2 淋溶試驗

裝填好土柱之后，通過模擬自然降雨淋溶，從2016年5月13日開始，每半個月模擬淋洗一次降雨，每次500 ml，共淋洗12次，累計淋溶量為6L，淋洗液的總量相當于1 178 mm的降雨量，與當?shù)氐钠骄涤炅恳恢?。淋溶液的pH調至6.0，與當?shù)赜晁膒H大致相同。淋洗開始后，用500 ml的燒杯收集淋濾液，直到淋濾液滴完為止，將收集到的淋濾液用真空抽濾泵過 0.45 μm玻璃纖維濾膜后，保存在 4 ℃條件下，待測土壤淋濾液DOC（保存 1～2 d）。

在10、25、50 cm土層采集土壤樣品，2016年5月初開始采集土樣，采集頻率為每個月 1 次，每次樣品采集后， 將相同區(qū)組的同一土壤厚度下土壤樣品混合，每個處理 3個重復，土壤樣品經(jīng)自然風干后，挑去肉眼可見的細根和石塊，過 2 mm篩研磨保存?zhèn)溆谩?/p>

1.3 測定方法

土壤樣品中有機碳含量利用重鉻酸鉀-外加熱法分析;淋濾液樣品先加入磷酸前處理去除無機碳對DOC測量結果的干擾[22]，淋濾液DOC采用Elementar High TOCⅡ型總有機碳/氮分析儀分析。

試驗數(shù)據(jù)采用SPSS 18.0 統(tǒng)計分析軟件和 Excel 軟件進行處理。

2 結果與分析

2.1 不同土壤厚度對土壤有機碳與DOC淋失的影響

通過模擬自然降雨，氮沉降影響下3種土壤厚度的土壤有機碳變化情況（表1），研究表明3種氮沉降處理土壤有機碳的含量隨著氮沉降量的增加而增加，LN、MN和HN處理的土壤有機碳含量大小順序表現(xiàn)為LN10～25 cm>25～50 cm，以0～10 cm土壤厚度的土壤有機碳含量最大。不同氮沉降處理土壤有機碳含量之間差異極顯著（P<0.01）。

土壤DOC是土壤中碳元素的一種重要存在形式，經(jīng)過分析不同土壤厚度處理土壤有機碳含量與DOC淋失量變化之間的關系，根據(jù)圖1可以看出，兩者之間的變化表現(xiàn)出顯著的正相關關系（P<0.01）。經(jīng)線性回歸擬合得到，0～10 cm土層厚度下為：y=16.998x-31.185（R2=0.732 2）;10～25 cm土層厚度下為：y=13.077x-35.199（R2=0.753 6）;25～50 cm土層厚度下為：y=6.976 9x-6.532 3（R2=0.707 7）。根據(jù)表1，土壤有機碳含量變化是：0～10 cm>10～2 5 cm>25～50 cm，同時，經(jīng)擬合方程可知，土壤DOC淋失變化也是：0～10 cm>10～25 cm>25～50 cm。這表明，土壤有機碳含量的變化在一定程度上制約著土壤DOC 淋失的變化。同時，比較土壤有機碳含量與土壤DOC淋失量的土層厚度變異可知，不同氮沉降處理土壤DOC淋失幅度隨土層厚度的變異均遠大于土壤有機碳，與土壤有機碳的變化相比，土壤DOC淋失量的變化更為敏感和迅速。

2.2 十二次降雨土壤DOC淋失過程

土柱淋溶系統(tǒng)中淋濾液DOC的遷移狀況響應于降雨次數(shù)的動態(tài)。根據(jù)圖2可知，十二次降雨的動態(tài)變化，在10 cm土層HN、MN、LN處理土壤DOC的淋失量隨降雨次數(shù)增多，逐漸減少，但累計淋失量逐漸增大。25 cm土層土壤DOC淋失量隨降雨次數(shù)增多而減小，后再第7次時出現(xiàn)峰值，然后又減小，但累計淋失量逐漸增大（圖3）。50 cm土層土壤DO淋失量隨降雨次數(shù)增多緩慢減小，也在第7次時出現(xiàn)峰值，隨后逐漸趨于平緩，累計淋失量逐漸增大。3種土壤厚度淋濾液DOC的淋出量高低順序為：10 cm>25 cm>50 cm，說明土壤厚度增加土壤DOC淋失量減少。

2.3 不同土壤厚度下氮沉降對土壤DOC淋失影響

人工模擬氮沉降對淋濾液DOC的淋失量影響（圖3），由圖3可知，在10 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，其中以HN氮沉降淋濾液DOC的淋失量最大為204.85 mg，HN、MN、LN處理之間差異不顯著;在25 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，其中以HN氮沉降淋濾液DOC的淋失量最大為175.35 mg，HN氮沉降淋濾液DOC的淋失量最小為76.71 mg，HN與LN處理之間差異極顯著，MN與HN、LN差異顯著;在50 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，其中以HN氮沉降淋濾液DOC的淋失量最大為153.02 mg，HN與MN、LN處理之間差異顯著，LN與MN差異不顯著。由此可以看出不同土壤厚度下，土壤DOC淋失響應氮沉降處理變化不同。表明增施氮肥給土壤提供了豐富的有機物，提高了土壤DOC的含量，通過模擬降雨試驗，土壤DOC逐漸由土體向水體遷移，淋濾液DOC的淋出量增多。

3 討論

3.1 土壤DOC淋失對氮沉降的響應

土壤溶解有機碳（DOC）主要來源于林木下的枯落物、腐殖質分解、微生物生物量和根系分泌物等[23-24]。氮沉降進入土壤系統(tǒng)會使其發(fā)生一系列物理和化學反應，導致DOC遷移淋溶。本文研究得出，不同氮沉降處理下，土壤DOC淋失量順序：HN>MN>LN，HN更容易促進DOC淋失，各處理之間差異不顯著。Findlay 等[25]曾報道，氮沉降可能導致土壤可溶性有機碳大量流失;以及Pregitzer[26]等，在美國北部闊葉林開展的一項連續(xù) 8 a的模擬氮沉降結果表明，氮輸入對土壤DOC的生產(chǎn)和淋溶有明顯的促進作用，研究結果與本文一致。原因可能是由于外加氮素滿足了微生物對氮素的需求[27]，此外，氮沉降增加也會使土壤中可溶性酚類物質逐漸積累[11]，從而增加土壤中DOC的含量，通過模擬降雨試驗，土壤DOC逐漸由土體向水體遷移，淋濾液DOC的淋出量增加。土壤活性有機碳占土壤有機碳總量的比例雖小，但它可以靈敏地反映土壤微小變化，可作為土壤有機碳有效性和土壤質量的早期指標[29-30]。因為土壤DOC是土壤有機碳的組成部分之一[31]，故兩者之間關系敏感。由圖4可以看出，土壤有機碳和DOC之間的變化表現(xiàn)出顯著的正相關關系（P<0.01），兩者相比較來說，DOC淋失幅度隨土層厚度的變異均遠大于土壤有機碳，與土壤有機碳的變化相比，土壤DOC淋失量的變化更為敏感和迅速。原因可能是土壤DOC具有穩(wěn)定性差、易氧化、易礦化、在土壤中移動快、容易受植物和微生物影響和活性高等特點[9，29]。

3.2 土壤厚度對DOC淋失的影響

由于特殊的氣候條件和長期的巖溶作用，使得喀斯特地區(qū)成土速度緩慢，土層極薄，巖石裸露率高，土地存在嚴重的不連續(xù)性和非均質性[32-33]。 土壤厚度是巖溶石漠化地區(qū)農(nóng)林業(yè)發(fā)展的主要限制因子，表征土壤的空間大小[34]，土壤厚度是土壤的一個重要基本特性，能直接反映土壤的發(fā)育程度，是野外鑒別土壤肥力的重要指標，也是喀斯特地貌土地生產(chǎn)力的決定因子之一[11]。研究表明，土壤厚度對DOC淋失量有極顯著的影響，DOC淋失量隨土壤厚度增加而減少，且不同土壤厚度之間差異顯著（P<0.05）（圖4）。3種土壤厚度的土壤淋濾液DOC的淋失量高低順序為：10 cm>25 cm>50 cm。因此土壤層次越薄，土壤DOC的淋失量就越大，表明土壤厚度對土壤DOC的淋失量有直接的影響。由于喀斯特地區(qū)土層厚度深淺不一[35]，對于土層厚的土壤，LN處理作用到深層土壤的能力小，HN處理作用到深層土壤能力強，故HN促進土壤DOC淋失。喀斯特地上—地下二元結構，導致土壤異質性較高，因此，不同土壤厚度DOC淋失對同一氮沉降處理的響應也有所不同。

4 結論

3種土壤厚度淋濾液DOC的淋出量高低順序為：10 cm>25 cm>50 cm，說明土壤厚度增加土壤DOC淋失量減少。在10 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，HN、MN、LN處理之間差異不顯著;在25 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，HN與LN處理之間差異極顯著，MN與HN、LN差異顯著;在50 cm土層時不同氮沉降處理土壤DOC淋失量順序：HN>MN>LN，HN與MN、LN處理之間差異顯著，LN與MN差異不顯著。由此可以看出氮沉降對土壤DOC 淋失有促進作用，喀斯特系統(tǒng)不同土壤厚度下，土壤DOC淋失響應氮沉降處理變化不同。

【參 考 文 獻】

[1]徐雷.氮沉降對杉木人工林土壤理化性質和酶活性的影響[D].福州：福建農(nóng)林大學，2007.

XU L. Effects of nitrogen deposition on soil physical-chemical properties and enzyme activities of Chinese fir plantation forest[D]. Fuzhou： Fujian Agriculture and Forestry University， 2007.

[2]樊后保.酸雨與森林衰退關系研究綜述[J].森林與環(huán)境學報，2003，23（1）：88-92.

FAN H B. Review of the relationship between acid rain and forest decline[J]. Journal of Forest and Environment， 2003， 23（1）： 88-92.

[3]方華，莫江明.氮沉降對森林凋落物分解的影響[J].生態(tài)學報，2006，26（9）：3127-3136.

FANG H， MO J M. Effects of nitrogen deposition on the decomposition of forest litter[J]. Chinese Journal of Ecology， 2006， 26（9）： 3127-3136.

[4]HOUGHTON R A. Terrestrial carbon sinks - uncertain explanations[J]. Biologist， 2002， 49（4）：155-160.

[5]劉德燕，宋長春.外源氮輸入對土壤有機碳礦化和凋落物分解的影響[J].土壤通報，2008，39（3）：675-680.

LIU D Y， SONG C C. Effects of exogenous nitrogen input on soil organic carbon mineralization and litter decomposition[J]. Chinese Journal of Soil Science， 2008， 39（3）： 675-680.

XING Y， LI X Q， FANG B， et al. Effects of biochar addition on DOC leaching in two types of soils[J]. Earth and Environment， 2015， 43（2）： 133-137.

[17]花可可，朱波，王小國，等，紫色土坡耕地可溶性有機碳徑流遷移特征[J].農(nóng)業(yè)工程學報，2013，29（5）： 81-88.

HUA K K， ZHU B， WANG X G， et al. Runoff migration characteristics of soluble organic carbon in purple soil slope farmland[J]. Transactions of the Chinese Society of Agricultural Engineering， 2013， 29（5）： 81-88.

[18]王世杰，季宏兵，歐陽自遠，等.碳酸鹽巖風化成土作用的初步研究[J].中國科學地球科學，1999，29（5）：441-449.

WANG S J， JI H B， OUYANG Z Y， et al. Preliminary study on the weathering of carbonate rock into soil[J]. Chinese Science Earth Science， 1999， 29（5）： 441-449.

[19]張殿發(fā)，王世杰，李瑞玲.貴州省喀斯特山區(qū)生態(tài)環(huán)境脆弱性研究[J].地理與地理信息科學，2002，18（1）：77-79.

ZHANG D F， WANG S J， LI R L. Study on the vulnerability of ecological environment in karst mountainous areas of Guizhou Province[J]. Journal of Geography and Geo-Information Science， 2002， 18（1）： 77-79.

[20]況野，龍偲，陳中吉.不同厚度土壤呼吸對氮沉降的響應[J].山地農(nóng)業(yè)生物學報，2014，33（3）：55-59.

KUANG Y， LONG S， CHEN Z J， et al. Response of soil respiration to nitrogen deposition in different thicknesses[J]. Journal of Mountain Agriculture Biology， 2014， 33（3）： 55-59.

[21]陶武輝.施肥對黃土高原水蝕風蝕交錯區(qū)植物產(chǎn)量及土壤礦質氮、水溶性有機碳/氮的影響[D].楊凌：西北農(nóng)林科技大學，2009.

TAO W H. Effects of fertilization on plant yield and soil mineral nitrogen and water-soluble organic carbon/nitrogen in water-erosion and wind-induced erosion zone on the Loess Plateau[D]. Yangling： Northwest A&F University， 2009.

[22]李劍，孫友寶，馬曉玲，陳建立，黃濤宏.離子色譜（IC）前處理法測定土壤中的總有機碳含量[J].環(huán)境化學，2014，33（8）：1425-1426.

LI J， SUN Y B， MA X L， et al. Determination of total organic carbon in soil by ion chromatography （IC） pretreatment[J]. Environmental Chemistry， 2014， 33（8）：1425-1426.

[23]KALBITZ K， SOLINGER S， PARK J， et al. Controls on the dynamics of dissolved organic matter in soils： A review[J]. Soil Science， 2000， 165（4）：277-304.

[24]SMOLANDER A， KITUNEN V. Soil microbial activities and characteristics of dissolved organic C and N in relation to tree species[J]. Soil Biology & Biochemistry， 2002， 34（5）： 651-660.

[25]FINDLAY S E G. Increased carbon transport in the Hudson River： Unexpected consequence of nitrogen deposition[J]. Frontiers in Ecology and the Environment， 2005， 3（3）：133-137.

[26]GUGGENBERGER G， ZECH W. Dissloved organic carbon controls in acid forest soils of the Fichtelgebirge （Germany） as revealed by distribution patterns and structural composition analyses[J]. Geoderma， 1993， 59（1-4）：109-129.

[27]莫江明，薛璟花，方運霆.鼎湖山主要森林植物凋落物分解及其對N沉降的響應[J].生態(tài)學報，2004，24（7）：1413-1420.

MO J M， XUE J H， FANG Y T. Litter decomposition of main forest plants in Dinghu Mountain and its response to N deposition[J]. Chinese Journal of Ecology， 2004， 24（7）： 1413-1420.

[28]張偉，王克林，陳洪松，等.典型喀斯特峰叢洼地土壤有機碳含量空間預測研究[J].土壤學報，2012，49（3）：601-606.

ZHANG W， WANG K L， CHEN H S， et al. Spatial prediction of soil organic carbon content in typical karst peak clusters[J]. Chinese Journal of Soil Science， 2012， 49（3）： 601-606.

[29]宇萬太，馬強，趙鑫.不同土地利用類型下土壤活性有機碳庫的變化[J].生態(tài)學雜志，2007，26（12）：2013-2016.

YU W T， MA Q， ZHAO X， et al. Changes of soil active organic carbon pool under different land use types[J]. Chinese Journal of Ecology， 2007， 26（12）： 2013-2016.

[30]張揚，裴亞蒙，任昊曄，等.模擬氮沉降對長白山興安落葉松林土壤有機碳庫的短期影響[J].林業(yè)科技，2019，44（3）：22-25.

ZHANG Y，PEI Y M，REN H Y，et al.Response of Soil Total Organic Carbon and Its Fractions to Short-term Simulated Nitrogen Deposition in a Larix gmelinii Plantation in Changbai Mountains[J].Forestry Science & Technology，2019，44（3）：22-25.

[31]羅梅，田冬，高明.紫色土壤有機碳活性組分對生物炭施用量的響應[J].環(huán)境科學，2018，39（9）：4327-4337.

LUO M， TIAN D， GAO M， et al. Response of organic carbon active components in purple soil to biochar application rate [J]. Environmental Science， 2018， 39（9）： 4327-4337.

[32]曹建華.袁道先.潘根興.巖溶生態(tài)系統(tǒng)中的土壤[J].地球科學進展，2003，18（1）：37-044.

CAO J H， YUAN D X， PAN G X. Soil in karst ecosystem[J]. Advances in Earth Science， 2003， 18（1）： 37-044.

[33]周運超，王世杰，盧紅梅.喀斯特石漠化過程中土壤的空間分布[J].地球與環(huán)境，2010，38（1）：1-7.

ZHOU Y C， WANG S J， LU H M. Spatial distribution of soil during karst rocky desertification[J]. Earth and Environment， 2010， 38（1）： 1-7.

[34]尹亮.西南巖溶高原石漠化地區(qū)小流域土壤厚度空間分異規(guī)律研究[D].長沙：湖南大學，2012.

YIN L. Spatial variation of soil thickness in small watershed in rocky desertification area of southwestern karst plateau[D]. Changsha： Hunan University， 2012.

[35]張珍明，周運超，李會.喀斯特小流域土壤有機碳分布特征及其影響因素[J].地球與環(huán)境，2017，45（1）： 38-45.

ZHANG Z M， ZHOU Y C， LI H， et al. Distribution characteristics of soil organic carbon and its influencing factors in karst watershed[J]. Earth and Environment， 2017， 45（1）： 38-45.