曹意 丁蘇雅 覃云斌 何昕諾 馬姜明

摘?要：為探究毛竹林下種植茶樹對(duì)土壤有機(jī)碳儲(chǔ)量與碳組分的影響，該研究以毛竹純林、竹茶混交林和常綠闊葉林為研究對(duì)象，采集這3種林分類型的表層（0～10 cm）土壤，測(cè)定土壤有機(jī)碳（SOC）、碳組分、生物與非生物因素指標(biāo)。結(jié)果表明：（1）竹茶混交林林下植物多樣性相較于毛竹純林顯著降低，但其土壤有機(jī)碳密度（22.54±2.09） t·hm-2、碳組分與毛竹純林無顯著差異（P>0.05）。竹茶混交林的礦物結(jié)合態(tài)有機(jī)碳（MOC）為（20.13±1.83） g·kg-1，占總有機(jī)碳的92.66%。常綠闊葉林土壤有機(jī)碳密度比竹茶混交林和毛竹純林高土壤有機(jī)碳密度分別高41.15%和41.00%（P<0.05）。（2）3種林分類型土壤微生物量碳（MBC）含量范圍為0.58～3.08 g·kg-1，土壤16S rRNA豐度范圍為2.18×1010 ～5.65×1010copies·g-1，固碳基因cbbL豐度范圍為0.37×108～ 1.10 ×108 copies·g-1，土壤微生物碳利用效率范圍為0.03～0.28；3種林分類型之間微生物相關(guān)指標(biāo)不存在顯著差異（P>0.05）。（3）3種林分類型SOC與土壤pH、砂粒含量和地上凋落物生物量呈顯著負(fù)相關(guān)，與土壤黏粒含量、粉粒含量、總氮、C∶N、總磷和銨態(tài)氮含量呈顯著正相關(guān)（P<0.05）。（4）就不同碳組分而言，顆粒有機(jī)碳（POC）和MOC均與土壤pH、砂粒含量和根系生物量呈顯著負(fù)相關(guān)，與土壤含水量、黏粒含量、粉粒含量、總氮、C∶N、總磷和銨態(tài)氮含量呈顯著正相關(guān)（P<0.05）。綜上表明，竹茶混交改造會(huì)造成原生毛竹純林林下植被多樣性下降，但并未造成土壤碳儲(chǔ)量下降；而相較于常綠闊葉林，毛竹經(jīng)營(yíng)措施需要改進(jìn)，以提升其碳匯效益。

關(guān)鍵詞： 顆粒有機(jī)碳， 礦物結(jié)合有機(jī)碳， 土壤有機(jī)碳密度， 毛竹林， 竹茶混交林

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

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-3142（2023）09-1668-10

收稿日期：2022-07-10

基金項(xiàng)目：廣西科技重大專項(xiàng)課題（桂科AA20161002-1）； 廣西重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（桂科AB21220057）； 廣西科技基地和人才專項(xiàng)（桂科AD21220163）。

第一作者： 曹意（1996-），碩士研究生，主要從事恢復(fù)生態(tài)學(xué)研究，（E-mail）cywiyii@163.com。

*通信作者：覃云斌，博士，主要從事森林生態(tài)學(xué)研究，（E-mail）shuibaoqinyunbin@163.com。

Effects of bamboo-tea mixed model on surface soil

organic carbon storage and components

CAO Yi1，2，3， DING Suya1，2，3， QIN Yunbin 1，2，3*， HE Xinnuo1，2，3， MA Jiangming1，2，3

（ 1. Key Laboratory of Ecology and Environmental Protection of Rare and Endangered Animals and Plants， Ministry of Education （Guangxi Normal

University）， Guilin 541004， Guangxi， China; 2. Guangxi Key Laboratory for Conservation and Sustainable Utilization of Landscape

Resources in Lijiang River Basin， Guilin 541006， Guangxi， China; 3. Institute of Sustainable Development

and Innovation， Guangxi Normal University， Guilin 541006， Guangxi， China ）

Abstract：In order to explore the effects of planting tea trees under Phyllostachys edulis forest on soil organic carbon storage and carbon components， we took pure bamboo forest， bamboo-tea mixed forest and evergreen broad-leaved forest as the study objects， collected the surface soil （0-10 cm） of these three forest types， and measured soil organic carbon（SOC）， carbon components， biotic and abiotic factors. The results were as follows： （1） The diversity of understory plants in bamboo-tea mixed forest was significantly lower than that in pure bamboo forest， but there was no significant differences in soil organic carbondensity［ （22.54 ± 2.09） t·hm-2 ］ and carbon composition between bamboo-tea mixed forest and pure bamboo forest （P>0.05）. Mineral-associated organic carbon （MOC） was （20.13 ± 1.83） g · kg-1， accounting for 92.66% of total organic carbon. The soil organic carbon density of evergreen broad-leaved forest was 41.15% and 41.00% higher than that of bamboo-tea mixed forest and pure bamboo forest （P>0.05）. （2） Soil microbial biomass carbon（MBC） content， 16S rRNA gene abundance， cbbL gene abundance and microbial carbon use efficiency of these forest types were （0.58 - 3.08） g·kg-1， （2.18×1010 - 5.65×1010） copies·g-1， （0.37×108-1.10×108） copies·g-1， 0.03 - 0.28， respectively. But there were significant differences about these microbial indicators between three forest types （P>0.05）. （3） SOC of the three forest types was significantly negatively correlated with soil pH， gravel content and aboveground litter biomass， and significantly positively correlated with soil clay content， silt content， total nitrogen， C∶N， total phosphorus and ammonium nitrogen （P<0.05）. （4） In terms of different carbon fractions， particulate organic carbon （POC） and MOC were significantly negatively correlated with soil pH， gravel content and root biomass， and significantly positively correlated with soil water content， clay content， silt content， total nitrogen， C∶N， total phosphorus and ammonium nitrogen （P<0.05）. In conclusion， the mixed transformation of bamboo and tea will cause the decrease of understory vegetation diversity in the original pure bamboo forest， but it will not cause the decrease of soil carbon storage; compared with evergreen broad-leaved forest， it is necessary to improve the management measures of bamboo to improve its carbon sink efficiency.

Key words： particulate organic carbon，mineral-associated organic carbon，soil orgnaic carbon density，Phyllostachys edulis forest，bamboo-tea mixed forest

隨著“碳中和”目標(biāo)的提出，未來40年中國(guó)會(huì)實(shí)施更有力的政策減少碳排放、增加碳匯。但是，合理政策的制定依賴于對(duì)增匯減排的科學(xué)認(rèn)知。土壤是陸地最大碳儲(chǔ)存庫(kù)，通過增加土壤碳儲(chǔ)量被認(rèn)為是實(shí)現(xiàn)“碳中和”目標(biāo)的有效途徑（楊元合等，2022）。近自然造林以及植被恢復(fù)的土壤碳匯效應(yīng)已經(jīng)被廣泛認(rèn)知（Don et al.， 2011; Lu et al.， 2018），如何提升經(jīng)濟(jì)林的土壤碳匯效益已成為當(dāng)前關(guān)注的熱點(diǎn)（Yang et al.， 2021b; Villa et al.， 2022）。

毛竹（Phyllostachys edulis）是我國(guó)價(jià)值較高的竹類品種，廣泛分布于南方山地丘陵地區(qū)，約占全國(guó)竹林面積的74%，其主要經(jīng)濟(jì)價(jià)值在于生產(chǎn)竹材與竹筍（楊傳寶等，2020）。毛竹在經(jīng)濟(jì)林種植中占據(jù)重要地位，以及快速生長(zhǎng)特性，使其碳匯效益受到了廣泛關(guān)注。毛竹地上部分的年均固碳量較高，為（8.13±2.15）Mg·hm-2（Yen & Lee， 2011）。目前，主要圍繞毛竹入侵與經(jīng)營(yíng)管理措施兩個(gè)方面開展毛竹土壤碳庫(kù)研究。相較于常綠闊葉林，毛竹的種植與入侵往往會(huì)造成對(duì)土壤碳儲(chǔ)量的下降（Wang et al.， 2019；祁雪連等，2021），同時(shí)引起顆粒有機(jī)碳（particulate organic carbon， POC）與礦物結(jié)合態(tài)有機(jī)碳（mineral-associated organic carbon， MOC）含量與占比的變化（Yang et al.， 2021b）。不同的管理措施（如間伐、林下植被清除）和混交模式（如竹-杉木混交）會(huì)對(duì)毛竹林土壤碳庫(kù)產(chǎn)生不同的影響（漆良華等，2013；李光敏等，2019；Yang et al.， 2021a）。Yang等（2021a）在浙江安吉的研究表明，粗放經(jīng)營(yíng)的毛竹林比無經(jīng)營(yíng)和集約經(jīng)營(yíng)的毛竹林具有更高的團(tuán)聚體穩(wěn)定性和有機(jī)碳儲(chǔ)量。生物與非生物因素對(duì)毛竹土壤碳庫(kù)影響因研究地區(qū)與經(jīng)營(yíng)管理措施不同也存在差異性（張厚喜等，2019；Yang et al.， 2021a）。

目前，毛竹單一的經(jīng)濟(jì)收入模式難以滿足農(nóng)民收入水平提高的熱切需求。同時(shí)，近年來受竹材加工業(yè)低迷的影響，毛竹種植戶的積極性也備受打擊。因此，發(fā)展林下經(jīng)濟(jì)成為了增加農(nóng)戶收入、推動(dòng)毛竹產(chǎn)業(yè)可持續(xù)發(fā)展的重要舉措（蔡春菊等，2018；Zhang et al.， 2019）。在關(guān)注毛竹林林下種植經(jīng)濟(jì)效益的同時(shí)，其生態(tài)效益也是研究的重點(diǎn)，涉及土壤養(yǎng)分、微生物、水土保持等方面（蔡澤宇等，2017；王勤等，2020；曹小青等，2022）。王勤等（2020）對(duì)不同毛竹-多花黃精復(fù)合模式的研究顯示，低毛竹立竹密度的竹-藥復(fù)合模式土壤有機(jī)碳（soil organic carbon， SOC）含量要比高密度毛竹立竹模式高13%。曹小青等（2022）的研究結(jié)果顯示，隨著套種年限的增加，毛竹-白及復(fù)合模式的土壤微生物多樣性逐漸提高。但是目前，毛竹林下復(fù)合經(jīng)營(yíng)的生態(tài)效益評(píng)估研究有限，并且多是關(guān)注毛竹與中草藥復(fù)合經(jīng)營(yíng)帶來的變化，缺乏對(duì)毛竹與茶樹復(fù)合經(jīng)營(yíng)的碳匯/源作用研究。同時(shí)，竹茶種植模式同樣是提升毛竹經(jīng)濟(jì)效益的重要林下經(jīng)濟(jì)舉措，特別是具有特色茶資源的竹產(chǎn)區(qū)。因此，開展竹茶混交模式的土壤碳庫(kù)變化研究對(duì)增進(jìn)該改造模式的生態(tài)效益認(rèn)知、促進(jìn)毛竹林改造的經(jīng)濟(jì)與生態(tài)效益協(xié)同提升建設(shè)具有重要意義。

目前，竹茶混交改造對(duì)毛竹林土壤碳庫(kù)儲(chǔ)量及其碳組分的影響研究相對(duì)缺乏。鑒于此，本文以廣西貓兒山地區(qū)毛竹純林、竹茶混交林和天然常綠闊葉林為研究對(duì)象，探究竹茶混交模式對(duì)毛竹林土壤碳儲(chǔ)量的影響，分析POC與MOC含量及其占比變化，揭示影響土壤碳儲(chǔ)量與碳組分變化的關(guān)鍵因素?；谕寥捞純?chǔ)量與碳組分評(píng)估，本研究成果可為提高毛竹林林下竹茶改造提供理論依據(jù)與技術(shù)支撐。

1?材料與方法

1.1 研究區(qū)概況

研究區(qū)位于廣西桂林市華江瑤族鄉(xiāng)，地處110°27′ E、25°48′ N，海拔在300～500 m之間，屬于中亞熱帶山地氣候，年均氣溫16.40 ℃，年均降水量在2 100 mm以上。成土母質(zhì)主要為花崗巖風(fēng)化母質(zhì)，土壤主要以紅壤和黃壤為主。華江瑤族鄉(xiāng)是全國(guó)十大毛竹生產(chǎn)基地之一，擁有毛竹林面積約（1.53×104）hm2且毛竹林面積有逐年擴(kuò)張的趨勢(shì)，該地區(qū)野生藥用植物和野生茶樹資源豐富（宜麗娜等，2012）。在鄉(xiāng)村振興的背景下，當(dāng)?shù)匾劳刑厣栀Y源開展竹茶復(fù)合經(jīng)營(yíng)的毛竹純林改造，以提升毛竹的經(jīng)濟(jì)效益。

1.2 研究方法

1.2.1 樣地設(shè)置?本研究共設(shè)置3種林分樣地類型，分別為毛竹純林（pure bamboo forest， PBF）、竹茶混交林（bamboo-tea mixed forest， BMF）和天然常綠闊葉林（evergreen broad-leaved forest， EBF）。每種林分類型設(shè)置6個(gè)立地條件接近的重復(fù)樣方（20 m × 20 m），并進(jìn)行群落學(xué)調(diào)查。常綠闊葉林與毛竹林分布在相同的海拔位置。3種林分類型的樣地之間至少間隔1 000 m，重復(fù)樣方之間至少間隔500 m。選取的毛竹純林樣地一直為單一種植毛竹，每年都會(huì)進(jìn)行周期性的間伐與林下植被清耕。選取的竹茶混交林樣地的喬木層為單一毛竹，每年會(huì)進(jìn)行毛竹間伐；其林下茶樹為自然生長(zhǎng)的野生茶樹，不存在人為管理，但每年春季會(huì)進(jìn)行采茶，茶樹因有一定經(jīng)濟(jì)價(jià)值得以保留。本研究以該自然存在的竹茶混交林作為毛竹純林改造的參照對(duì)象。3種林分樣地基本信息見表1。

1.2.2 樣地調(diào)查與樣品采集?在2021年11月對(duì)樣地進(jìn)行調(diào)查與采樣。調(diào)查內(nèi)容包括海拔、坡度、坡向、植物組成、經(jīng)營(yíng)措施等。其中，喬灌植物調(diào)查采用每木檢尺，草本植物調(diào)查采用小樣方法。土壤樣品采集則采用土鉆在每個(gè)樣方內(nèi)隨機(jī)采集5個(gè)0～10 cm子土樣，再混合成一個(gè)土樣，作為該樣地的土樣。同時(shí)，在每個(gè)采樣點(diǎn)附近收集土壤測(cè)定土壤容重，用鋁盒收集土壤測(cè)定含水率。在樣方內(nèi)隨機(jī)設(shè)置5個(gè)50 cm × 50 cm的小樣方收集地上凋落物，設(shè)置5個(gè)20 cm × 20 cm的小樣方收集0～10 cm層土壤根系。收集的土樣帶回實(shí)驗(yàn)室后，一部分自然風(fēng)干，一部分保存在4 ℃冰箱用于鮮土實(shí)驗(yàn)測(cè)定，一部分保存在-80 ℃下以便進(jìn)行微生物DNA提取。

1.2.3 土壤理化性質(zhì)與微生物碳代謝相關(guān)指標(biāo)測(cè)定?土壤全碳與總氮含量采用碳氮元素分析儀（Elementar Vario EL Ⅲ，德國(guó)）測(cè)定。土壤容重采用環(huán)刀法測(cè)定。pH采用電位法測(cè)定（水土質(zhì)量比為2.5∶1）?？偭祝╰otal phosphorus， TP）含量采用鉬銻鈧比色法測(cè)定。土壤速效磷（available phosphorus， AP）采用 1/2硫酸消煮法然后用鉬銻抗比色法測(cè)定。土壤銨態(tài)氮（NH4+-N）與硝態(tài)氮（NO3--N）含量，采用2 mol·L-1 KCl溶液浸提及SKALAR SAN++連續(xù)流動(dòng)分析儀測(cè)定。土壤機(jī)械組成測(cè)定采用比重計(jì)法，即50 g風(fēng)干土樣混合分散劑煮沸后，過濾至1 000 mL沉降筒內(nèi)，采用比重計(jì)進(jìn)行測(cè)定。采用國(guó)際制將土壤機(jī)械組成劃分為砂粒、粉粒和黏粒。土壤有機(jī)碳含量采用濃硫酸-重鉻酸鉀高溫外加熱氧化法測(cè)定。地上凋落物與根系生物量采用烘干法測(cè)定。采用六偏磷酸鈉分散法區(qū)分土壤顆粒態(tài)有機(jī)碳和礦物結(jié)合態(tài)有機(jī)碳。具體方法如下：自然風(fēng)干的土壤樣品加入0.5%的六偏磷酸鈉溶液（土液質(zhì)量比1∶3）和玻璃珠，振蕩后充分分散；利用53 μm篩過濾分組為礦物結(jié)合態(tài)有機(jī)質(zhì)（<53 μm）和顆粒態(tài)有機(jī)質(zhì)（>53 μm）；分離后的組分用去離子水洗后60 °C烘干，待測(cè)；結(jié)合土層厚度、容重等土壤理化性質(zhì)來計(jì)算土壤有機(jī)碳密度（soil organic carbon density， SOCD），公式如下。

SOCD=SOC×ρ×D×10 （1）

式中： SOC為土壤有機(jī)碳的平均含量（g·kg-1）； D為土層厚度（cm）； ρ為土壤容重（g·cm-3）。

土壤微生物量碳（microbial biomass carbon， MBC）含量采用氯仿熏蒸法測(cè)定（汪亞芳等，2022）。16S rRNA與cbbL基因采用qPCR測(cè)定，引物序列分別為F515/R907和K2f/V2f（Qin et al.， 2020）。采用 H218O 標(biāo)記土壤樣品中的微生物，測(cè)定微生物碳利用效率（carbon use efficiency， CUE），參考Zheng等（2019）方法。

1.2.4 數(shù)據(jù)處理與分析?采用 Excel 2016、SPSS 25和Origin 2021進(jìn)行數(shù)據(jù)整理與分析。采用單因素方差分析（ANOVA）分析不同處理組間變量顯著性差異。采用 Pearson 相關(guān)分析揭示變量間相互關(guān)系，通過多元線性回歸模型逐步法揭示土壤理化因素與微生物指標(biāo)對(duì)土壤有機(jī)碳的影響。所有統(tǒng)計(jì)的顯著性水平設(shè)置為P<0.05。

2?結(jié)果與分析

2.1 不同林分類型土壤理化性質(zhì)差異

3種林分類型中土壤pH，黏粒和砂粒含量，總磷及硝態(tài)氮含量均存在顯著差異（P<0.05）。其中，3種林分類型的土壤均為酸性土壤且竹茶混交林的土壤pH值最高。毛竹純林的土壤含水率最低且土壤的砂粒含量顯著高于其他2個(gè)林分類型（P<0.05），其土壤養(yǎng)分含量（包括總磷和硝態(tài)氮）均低于竹茶混交林和常綠闊葉林。竹茶混交林及毛竹純林的根系生物量、C∶N、總氮和銨態(tài)氮含量無顯著差異，但與常綠闊葉林存在顯著差異（P<0.05）。其中，二者的根系生物量顯著高于常綠闊葉林，而其他指標(biāo)則顯著低于常綠闊葉林（P<0.05）（表2）。

2.2 不同林分類型SOC與SOCD變化特征

竹茶混交林表層SOC為（22.95±1.91）g·kg-1，與毛竹純林SOC不存在顯著差異（P>0.05），比常綠闊葉林（39.74±5.82）g·kg-1低42.25%（P<0.05）。3種林分類型SOCD變化范圍為（19.65～44.09）t·hm-2。其中，竹茶混交林與毛竹純林表層SOCD分別為（22.54±2.09）、（22.60±2.53）t·hm-2，均顯著低于常綠闊葉林（38.31±5.40）t·hm-2（P<0.05）（圖1）。

2.3 不同林分類型有機(jī)碳組分含量變化特征

3種林分類型土壤均以礦物結(jié)合態(tài)有機(jī)碳為主，MOC的范圍為17.98～36.83 g·kg-1，對(duì)SOC的貢獻(xiàn)占90.19%～94.40%，剩余貢獻(xiàn)部分為POC；不同林分類型之間POC和MOC對(duì)SOC的貢獻(xiàn)不存在顯著差異（P>0.05）。由表3相關(guān)性分析顯示，SOC與POC、MOC呈極顯著正相關(guān)，相關(guān)系數(shù)分別為0.95、0.97（P<0.05）。常綠闊葉林土壤POC和MOC分別為（3.14±0.42）g·kg-1和（32.09±3.54）g·kg-1，均分別顯著高于其毛竹純林和竹茶混交林（P<0.05）。竹茶混交林POC和MOC分別為（1.59±0.32）g·kg-1和（20.13±1.83）g·kg-1，與毛竹純林不存在顯著差異（P>0.05）（圖2）。

2.4 不同林分類型微生物碳代謝相關(guān)指標(biāo)變化特征

3種林分類型土壤微生物量碳含量范圍為0.58～3.08 g·kg-1，其與SOC之間相關(guān)性并不明顯（P>0.05）（表3）。不同林分類型之間微生物量碳含量不存在顯著差異（P>0.05），其中竹茶混交林土壤微生物量碳含量為（1.47±0.34）g·kg-1。竹茶混交林土壤16S rRNA豐度平均為3.91×1010copies·g-1；固碳基因cbbL豐度為0.76×108copies·g-1，占總細(xì)菌豐度的0.21%。竹茶混交林與毛竹純林和常綠闊葉林在16S rRNA豐度和cbbL豐度上不存在顯著差異（P>0.05）。3種林分類型土壤微生物碳利用效率范圍為0.03～0.28，3種林分類型之間微生物碳利用效率（carbon use efficiency， CUE）不存在顯著差異（P>0.05）（表4）。

2.5 各有機(jī)碳組分與環(huán)境因子相關(guān)性

由表5可知，3種林分類型SOC與土壤pH、砂粒含量和地上凋落物生物量呈顯著負(fù)相關(guān)， 與土壤黏粒含量、粉粒含量、總氮、C∶N、總磷和銨態(tài)氮含量呈顯著正相關(guān) （P<0.05）。就不同碳組分而言，POC和MOC均與土壤pH、砂粒含量和根系生物量呈顯著負(fù)相關(guān)，與土壤含水率、黏粒含量、粉粒含量、總氮、C∶N、總磷和銨態(tài)氮含量呈顯著正相關(guān)（P<0.05）（表5）。多元線性回歸結(jié)果顯示，土壤總氮含量與C∶N是預(yù)測(cè)土壤碳密度的關(guān)鍵因素，可解釋土壤碳密度變化的97%，公式如下。

y=-22.12+11.66CN+1.81RC∶N （R2=0.97， P<0.01） （2）

式中：y為土壤碳密度（t·hm-2）; CN為總氮含量（g·kg-1）; RC∶N為碳氮質(zhì)量比。

3?討論

3.1 毛竹林改造對(duì)土壤碳儲(chǔ)量和組分的影響

在研究區(qū)內(nèi)開展毛竹純林的竹茶混交模式改造，從長(zhǎng)期來看并不會(huì)引起土壤有機(jī)碳儲(chǔ)量的下降。研究結(jié)果顯示，毛竹純林與竹茶混交林SOCD不存在顯著差異。以往對(duì)喬木樹種與林下植被多樣性研究表明，物種多樣性與土壤碳儲(chǔ)量多呈正相關(guān)（Zhao et al.， 2014；Chen et al.， 2020）。這種差異性可能是由于毛竹純林與竹茶混交林的碳循環(huán)均主要受毛竹的支配。這兩種林分類型喬木層只有毛竹，而毛竹根系發(fā)達(dá)，占據(jù)著表層土壤，影響著土壤有機(jī)碳與其他養(yǎng)分循環(huán)。劉駿等（2013）研究表明，毛竹林細(xì)根生物量是常綠闊葉林的5.86倍。雖然竹茶模式改造使毛竹的株密度下降了約20%，但兩者的表層土壤根系生物量并不存在顯著性差異。竹茶混交改造遺留的地下空間會(huì)迅速被其他毛竹根系占據(jù)。此外在毛竹砍伐期間，種植者會(huì)清除林下植被，這種周期性更新也削弱了毛竹純林內(nèi)植物源SOC的輸入與積累。李光敏等（2019）的研究結(jié)果顯示，毛竹林林下植被清除2年后其SOC與對(duì)照組之間不存在顯著差異。在毛竹的支配性地位與當(dāng)前管理方式的影響下，毛竹純林的竹茶混交模式改造對(duì)表層土壤碳儲(chǔ)量的影響有限。

雖然毛竹林林下改造不會(huì)造成土壤碳儲(chǔ)量的下降，但這種種植模式并不利于SOC的固持。本研究顯示，毛竹純林與竹茶混交林的SOCD分別比常綠闊葉林低41.15%和41.00% （P<0.05）。這說明將常綠闊葉林轉(zhuǎn)化為毛竹純林和竹茶混交林均會(huì)造成土壤碳儲(chǔ)量的流失，這與其他地區(qū)的研究結(jié)果一致（Lin et al.， 2018；Wang et al.， 2019；祁雪連等，2021）。Lin等（2018）通過對(duì)常綠闊葉林轉(zhuǎn)為毛竹林的對(duì)比研究表明，這種土地變化會(huì)使0～40 cm層土壤碳儲(chǔ)量下降12%。毛竹經(jīng)營(yíng)管理、喬木層樹種差異及其影響下的土壤理化性質(zhì)差異是造成毛竹純林、竹茶混交林與常綠闊葉林土壤碳儲(chǔ)量存在差異的主要原因。因此，在竹茶混交改造過程中探索如何提高土壤有機(jī)碳儲(chǔ)量十分必要。

本研究表明，竹茶混交改造不會(huì)引起毛竹林土壤有機(jī)碳組分含量與占比的顯著性變化。3種林分類型土壤均以MOC為主，占總有機(jī)碳的比例高達(dá)90%以上。這說明這些林地土壤有機(jī)碳庫(kù)以穩(wěn)定性碳為主，但也反映出POC較難保留。每年植物會(huì)通過地上凋落物、根系分泌物等形式向土壤釋放碳源，再經(jīng)分解轉(zhuǎn)化為POC；而3種林分類型POC占比極小，大部分新輸入的碳不可能馬上轉(zhuǎn)化為穩(wěn)定性的MOC，極可能是被損失掉。POC的損失可能由該地區(qū)充沛的降雨造成。較強(qiáng)的降水既能將植物凋落物沖走，減少其原地轉(zhuǎn)化為POC的可能，同時(shí)也能將不穩(wěn)定的POC隨徑流帶走，減少其保留。董玉清等（2020）在貓兒山相近海拔的常綠闊葉林研究中，也發(fā)現(xiàn)POC占總有機(jī)碳的比例在1.26%～14.44%之間。因此，有必要在這個(gè)地區(qū)開展碳示蹤研究，以揭示碳的轉(zhuǎn)化歸趨。

3.2 影響土壤有機(jī)碳及其組分變化的主要因素

土壤有機(jī)碳的轉(zhuǎn)化和積累與生物、非生物因素存在密切聯(lián)系。本研究中土壤有機(jī)碳及其組分與總氮、C∶N、總磷和銨態(tài)氮含量呈極顯著正相關(guān)（P<0.01），并且回歸分析顯示總氮是影響有機(jī)碳及其組分的最關(guān)鍵因素。微生物生長(zhǎng)需要氮源，當(dāng)?shù)慈狈r(shí)其會(huì)選擇降解凋落物與有機(jī)質(zhì)以獲取氮源。這個(gè)過程會(huì)造成有機(jī)碳含量的下降（Frey et al.， 2014；Lu et al.， 2021）。因此，充足的氮源有助于SOC的積累。在一些研究中，同樣發(fā)現(xiàn)土壤總氮和無機(jī)氮含量與有機(jī)碳含量呈正相關(guān)的現(xiàn)象（李萍，2010；Yang et al.， 2021b）。李萍（2010）對(duì)江西大崗山毛竹林的研究同樣發(fā)現(xiàn)土壤總氮是影響土壤有機(jī)碳變化的最大因素且彼此之間呈正相關(guān)。因此，在毛竹的管理中可以適當(dāng)添加氮肥，改變毛竹林土壤有機(jī)碳較低的現(xiàn)象，以提高土壤碳匯能力。

整體而言，3種林分類型SOC和各碳組分含量與土壤pH呈極顯著負(fù)相關(guān)（P<0.01）。這說明pH越低越有利于土壤有機(jī)碳的積累。土壤 pH可以通過影響微生物群落結(jié)構(gòu)與活性而影響有機(jī)碳的分解與轉(zhuǎn)化。pH過高與過低均不利于微生物的生長(zhǎng)與活動(dòng)，從而抑制有機(jī)碳的分解轉(zhuǎn)化（Zhang et al.， 2019）。張厚喜等（2019）對(duì)武夷山不同海拔梯度毛竹林土壤有機(jī)碳特征及影響因素中也發(fā)現(xiàn)相同現(xiàn)象。本研究中SOC及其組分與土壤黏粒含量、粉粒含量呈極顯著正相關(guān)（P<0.01）。這主要是因?yàn)榉哿：宛ち?duì)SOC存在保護(hù)作用，減少碳的礦化分解與流失（Lehmann & Kleber，2015；Kasmerchak et al.，2018）。但是，竹茶混交林比毛竹純林具有更高的黏粒和粉粒含量，而兩者之間的SOCD并無顯著差異。這說明黏粒和粉粒含量可能需要達(dá)到一定閾值才會(huì)影響到有機(jī)碳的積累，也可能是其他關(guān)鍵因素抵消了竹茶混交林黏粒和粉粒對(duì)有機(jī)碳積累的貢獻(xiàn)。因此，毛竹林林下種植茶樹至少可以改善土壤機(jī)械組成，朝有利于有機(jī)碳積累的方向發(fā)展。

土壤微生物是影響土壤碳循環(huán)的關(guān)鍵因素，參與SOC的積累與分解。但是，基于對(duì)土壤微生物量碳含量、固碳基因豐度、碳利用效率指標(biāo)的分析顯示，本研究區(qū)內(nèi)微生物活動(dòng)不是影響土壤有機(jī)碳儲(chǔ)量變化的決定因素。汪亞芳等（2022）研究同樣發(fā)現(xiàn)毛竹林與常綠闊葉林土壤微生物量之間不存在顯著差異。此外，Bai等（2016）的研究還發(fā)現(xiàn)毛竹入侵闊葉林會(huì)造成SOC儲(chǔ)量下降，但微生物量碳得到了提高。因此，還需要開展動(dòng)態(tài)和多季節(jié)性研究來揭示微生物對(duì)土壤碳庫(kù)的影響。

4?結(jié)論

本研究結(jié)果表明，相較于毛竹純林，毛竹林林下種植茶樹會(huì)顯著提高土壤pH、含水率、黏粒、粉粒含量、總磷與硝態(tài)氮含量。毛竹林下茶樹的種植會(huì)造成林下灌木層植物多樣性下降，但不會(huì)造成土壤有機(jī)碳密度的下降，以及有機(jī)碳組分比例的改變。毛竹純林與竹茶混交林的土壤有機(jī)碳儲(chǔ)量顯著低于常綠闊葉林。3種林分的碳組分均以礦物結(jié)合態(tài)有機(jī)碳為主。因此，需要改進(jìn)現(xiàn)有的毛竹純林與竹茶混交林經(jīng)營(yíng)管理模式，在提升經(jīng)濟(jì)效益的同時(shí)，注重碳匯效益的協(xié)同發(fā)展。此外，本研究從土壤碳匯的角度反映生態(tài)效益，既有一定的代表性，也有一定的局限性。未來還需要開展更多元的生態(tài)效益評(píng)估，以揭示竹茶改造模式對(duì)毛竹純林生態(tài)系統(tǒng)功能的影響，從而為竹茶改造模式提供科學(xué)依據(jù)與指導(dǎo)。

參考文獻(xiàn)：

BAI SB， CONANT RT， ZHOU GM， et al.， 2016. Effects of Moso bamboo encroachment into native， broad-leaved forests on soil carbon and nitrogen pools ［J］. Sci Rep， 6（1）： 31480.

CAI CJ， FAN SH， LIU GL， et al.， 2018. Research and development advance of compound management of bamboo forests ［J］. World Bamboo Ratt， 16（5）： 47-52. ［蔡春菊， 范少輝， 劉廣路， 等， 2018. 竹林復(fù)合經(jīng)營(yíng)研究和發(fā)展現(xiàn)狀［J］. 世界竹藤通訊， 16（5）： 47-52.］

CAI ZY， ZHANG JF， SUN SY， et al.， 2017. Influence of Phyllostachy edulis stand improvement on runoff and nitrogen and phosphorus losses in Anji， Zhejiang， China［J］. J Soil Water Conserv， 31（6）： 46-51. ［蔡澤宇， 張建鋒， 孫士詠， 等， 2017. 安吉毛竹林分改造對(duì)地表徑流及氮磷流失的影響［J］. 水土保持學(xué)報(bào)， 31（6）： 46-51.］

CAO XQ， WANG L， SUN MY， et al.， 2022. Soil microbial community diversity in Phyllostachys pubescens-Bletilla striata ecosystems with different intercropping years［J］. SoilFert Sci Chin， （1）： 147-154. ［曹小青， 王亮， 孫孟瑤， 等， 2022. 不同年限毛竹-白及復(fù)合系統(tǒng)土壤微生物群落多樣性特點(diǎn)［J］. 中國(guó)土壤與肥料， （1）： 147-154.］

CHEN XL， CHEN HYH， CHEN C， et al.， 2020. Effects of plant diversity on soil carbon in diverse ecosystems： a global meta-analysis ［J］. Biol Rev， 95： 167-183.

DON A， SCHUMACHER J， FREIBAUER A， 2011. Impact of tropical land-use change on soil organic carbon stocks — a meta-analysis ［J］. Glob Change Biol， 17（4）： 1658-1670.

DONG YQ， GUAN P， LU Y， et al.，2020. Fractionsand contents of soil organic carbon at different elevations in Maoer Mountain ［J］. Chin J Soil Sci， 51（5）： 1142-1151. ［董玉清， 官鵬， 盧瑛， 等， 2020. 貓兒山不同海拔土壤有機(jī)碳組分構(gòu)成及含量特征 ［J］. 土壤通報(bào)， 51（5）： 1142-1151.］

FREY SD， OLLINGER S， NADELHOFFER K， et al.， 2014. Chronic nitrogen additions suppress decomposition and sequester soil carbon in temperate forests ［J］. Biogeochemistry， 121（2）： 305-316.

KASMERCHAK CS， MASON JA， LIANG MY， 2018. Laser diffraction analysis of aggregate stability and disintegration in forest and grassland soils of northern Minnesota， USA ［J］. Geoderma， 338： 430-444.

LEHMANN J， KLEBER M， 2015. The contentious nature of soil organic matter ［J］. Nature， 528： 60-68.

LI GM， CHEN FS， XU ZW， et al.， 2019. Effects of thinning and understory removal on soil labile organic carbon in moso plantation ［J］. Acta Agric Univ Jiangxi， 41（4）： 733-740. ［李光敏， 陳伏生， 徐志文， 等， 2019. 間伐和林下植被剔除對(duì)毛竹林土壤活性有機(jī)碳的影響 ［J］. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào)， 41（4）： 733-740.］

LI P， 2010. Soil organic carbon and its relationship with soil properties in Phyllostachys forests ［D］. Beijing： Beijing Forestry University. ［李萍， 2010. 毛竹林土壤有機(jī)碳變化及其與土壤性質(zhì)的關(guān)系［D］. 北京： 北京林業(yè)大學(xué).］

LIN ZW， LI YF， TANG CX， et al.， 2018. Converting natural evergreen broadleaf forests to intensively managed Moso bamboo plantations affects the pool size and stability of soil organic carbon and enzyme activities ［J］ Biol Fert Soils， 54（4）： 467-480.

LIU J， YANG QP， YU DK， et al.， 2013. Contribution of fine root to soil nutrient heterogeneity at two sides of the bamboo and broadleaved forest interface ［J］. Chin J Plant Ecol， 37（8）： 739-749. ［劉駿， 楊清培， 余定坤， 等， 2013. 細(xì)根對(duì)竹林-闊葉林界面兩側(cè)土壤養(yǎng)分異質(zhì)性形成的貢獻(xiàn) ［J］. 植物生態(tài)學(xué)報(bào)， 37（8）： 739-749.］

LU F， HU HF， SUN WJ， et al.， 2018. Effect of national ecological restoration projects on carbon sequestration in China from 2001 to 2010 ［J］. PNAS， 115： 4039-4044.

LU XK， VITOUSEK PM， MAO QG， et al.， 2021. Nitrogen deposition accelerates soil carbon sequestration in tropical forests ［J］. PNAS， 118（16）： e2020790118.

QI LH， FAN SH， DU MY， et al.， 2013. Vertical distribution and seasonal dynamicsof soil organic carbon in Phyllostachy edulis plantations and P. edulis-Cunninghamia lanceolata mixed forests in the hilly region of Central Hunan， Southern China［J］. Sci Silv Sin， 49（3）： 17-24. ［漆良華， 范少輝， 杜滿義， 等， 2013. 湘中丘陵區(qū)毛竹純林、毛竹-杉木混交林土壤有機(jī)碳垂直分布與季節(jié)動(dòng)態(tài)［J］. 林業(yè)科學(xué)， 49（3）： 17-24.］

QI XL， GE XM， QIAN ZZ， et al.， 2021. Differences of soil properties between natural mixed coniferous and broad-leaved forest and Moso bamboo plantation in Wuyi Mountains ［J］. Ecol Environ Sci， 30（8）： 1599-1606. ［祁雪連， 葛曉敏， 錢壯壯， 等， 2021. 武夷山天然針闊混交林與毛竹人工林土壤性質(zhì)差異 ［J］. 生態(tài)環(huán)境學(xué)報(bào)， 30（8）： 1599-1606.］

QIN YB， CHEN ZH， DING BJ， et al.， 2020. Impact of sand mining on the carbon sequestration and nitrogen removal ability of soil in the riparian area of Lijiang River， China ［J］. Environ Poll， 261： 114220.

VILLA MVE， CRISSTIANO PM， DE D， et al.， 2022. Do selective logging and pine plantations in humid subtropical forests affect aboveground primary productivity as well as carbon and nutrients transfer to soil？ ［J］. For Ecol Manage， 503（1） ： 119736.

WANG Q， SUN MY， DI JH， et al.， 2020. Effects of different Phyllostachys pubescens-Polygonatum sibiricum intercropping models on soil physicochemical properties ［J］. Ecol Sci， 39（6）： 54-59. ［王勤， 孫夢(mèng)瑤， 遆建航， 等， 2020. 毛竹-多花黃精復(fù)合經(jīng)營(yíng)模式對(duì)土壤理化特性的影響［J］. 生態(tài)科學(xué)， 39（6）： 54-59.］

WANG YF， LIU ZY， ZHANG BG， et al.， 2022.Effects of invasive Moso bamboo removal on soil microbial biomass and enzyme activities in subtropical forests ［J］. Chin J Appl Ecol， 33（5）： 1233-1239. ［汪亞芳， 劉宗悅， 張寶剛， 等， 2022. 入侵毛竹皆伐對(duì)亞熱帶森林土壤微生物生物量和酶活性的影響 ［J］. 應(yīng)用生態(tài)學(xué)報(bào)， 33（5）： 1233-1239.］

WANG HB， JIN J， YU PY， et al.， 2019. Converting evergreen broad-leaved forests into tea and Moso bamboo plantations affects labile carbon pools and the chemical composition of soil organic carbon ［J］. Sci Total Environ， 711（1）： 135225.

YANG CB， NI HJ， SU WH， et al.， 2020. Effects of management measures on organic carbon， nitrogen and chemical structure of different soil fractions in Phyllostachys edulis plantations ［J］. Chin J Appl Ecol， 31（1） ： 25-34. ［楊傳寶， 倪惠菁， 蘇文會(huì)， 等， 2020. 經(jīng)營(yíng)措施對(duì)毛竹林土壤不同組分有機(jī)碳、氮及化學(xué)結(jié)構(gòu)的影響 ［J］. 應(yīng)用生態(tài)學(xué)報(bào)， 31（1）： 25-34.］

YANG CB， WANG AK， ZHU ZX， et al.， 2021a. Impact of extensive management system on soil properties and carbon sequestration under an age chronosequence of Moso bamboo plantations in subtropical China ［J］. For Ecol Manage， 497（1）： 119535.

YANG CB， ZHANG XP， NI HJ， et al.， 2021b. Soil carbon and associated bacterial community shifts driven by fine root traits along a chronosequence of Moso bamboo （Phyllostachys edulis） plantations in subtropical China ［J］. Sci Total Environ， 752（15）： 142333.

YANG YH， SHI Y， SUN WJ， et al.， 2022. Terrestrial carbon sinks in China and around the world and their contribution to carbon neutrality ［J］. Sci Sin （Vit）， 52（4）： 534-574. ［楊元合， 石岳， 孫文娟， 等， 2022. 中國(guó)及全球陸地生態(tài)系統(tǒng)碳源匯特征及其對(duì)碳中和的貢獻(xiàn) ［J］. 中國(guó)科學(xué)（生命科學(xué)）， 52（4）： 534-574.］

YEN TM， LEE JS， 2011. Comparing aboveground carbon sequestration between Moso bamboo （Phyllostachys heterocycla） and China fir （Cunninghamia lanceolata） forests based on the allometric model ［J］. For Ecol Manage， 261（6） ： 995-1002.

YI LN， BAI KD， LI MX， et al.， 2012. Comparison of physiology and biochemical index of some evergreen and deciduous broad-leaved species in Maoer Mountain， Guangxi， China ［J］. J Guangxi Norm Univ （Nat Sci Ed）， 30（1）： 83-88. ［宜麗娜， 白坤棟， 李明霞， 等， 2012. 貓兒山常綠和落葉闊葉樹抗性生理生化指標(biāo)的比較 ［J］. 廣西師范大學(xué)學(xué)報(bào) （自然科學(xué)版）， 30（1）： 83-88.］

ZHAO FZ， KANG D， HAN XH， et al.， 2014. Soil stoichiometry and carbon storage in long-term afforestation soil affected by understory vegetation diversity ［J］. Ecol Eng， 74： 415-422.

ZHANG HX， LIN C， CHENG H， et al.， 2019. Variation of soil organic carbon content of Moso bamboo forest along altitudinal gradient in Wuyi Mountain in China ［J］. Soils， 51（4）： 821-828. ［張厚喜， 林叢， 程浩， 等， 2019. 武夷山不同海拔梯度毛竹林土壤有機(jī)碳特征及影響因素 ［J］. 土壤， 51（4）： 821-828.］

ZHANG XP， GUO GB， WU ZZ， et al.， 2019. Agroforestry alters the rhizosphere soil bacterial and fungal communities of Moso bamboo plantations in subtropical China ［J］. Appl Soil Ecol， 143： 192-200.

ZHENG Q， HU YT， ZHENG SS， et al.， 2019. Growth explains microbial carbon use efficiency across soils differing in land use and geology ［J］. Soil Biol Biochem， 128： 45-55.

（責(zé)任編輯?李?莉?王登惠）