• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    馬尾松純林闊葉化改造對(duì)土壤碳氮固持的短期效應(yīng)

    2024-12-31 00:00:00王浩東陳夢(mèng)袁叢軍何爽丁訪(fǎng)軍楊瑞

    摘 要:【目的】為了解馬尾松純林補(bǔ)植不同闊葉樹(shù)種對(duì)土壤碳庫(kù)的影響,篩選不同類(lèi)型的闊葉化改造土壤高效固碳模式?!痉椒ā吭谫F州省獨(dú)山縣國(guó)有林場(chǎng)內(nèi)選擇林分結(jié)構(gòu)相似且具有代表性的馬尾松Pinus massoniana人工純林,經(jīng)擇伐后補(bǔ)植香樟Cinnamomum camphora、楠木Phoebe zhennan、南酸棗Choerospondias axillaris、鵝掌楸Liriodendron chinense闊葉樹(shù)種8~14 a,以未補(bǔ)植闊葉樹(shù)的馬尾松純林為對(duì)照,分別設(shè)置3塊樣地,共15塊,采集0~20、20~40和40~60 cm的土壤,測(cè)定土壤理化性質(zhì)和活性有機(jī)碳組分?!窘Y(jié)果】補(bǔ)植不同闊葉樹(shù)種對(duì)土壤碳庫(kù)影響不同,對(duì)0~20 cm影響較為顯著,對(duì)更深層土壤影響較小,0~60 cm的土壤碳儲(chǔ)量變化范圍約1.39~12.77 kg·m-2,其中,馬尾松純林和馬尾松+南酸棗林的土壤碳儲(chǔ)量較高,分別達(dá)到212.35和203.51 kg·m-2,馬尾松+香樟林土壤碳儲(chǔ)量最低,約為100.78 kg·m-2;補(bǔ)植闊葉樹(shù)種后土壤pH值顯著降低,而MBC顯著增加;在0~20 cm土層,南酸棗+馬尾松林土壤有機(jī)碳和全氮顯著提升;補(bǔ)植闊葉樹(shù)種后短期內(nèi)土壤碳儲(chǔ)量降低,但土壤氮儲(chǔ)量無(wú)顯著變化;補(bǔ)植楠木后土壤碳庫(kù)活度和碳庫(kù)活度指數(shù)提升,補(bǔ)植南酸棗后土壤碳庫(kù)穩(wěn)定性顯著提升?!窘Y(jié)論】?jī)?yōu)先選擇葉片中初始木質(zhì)素/氮值高的落葉闊葉樹(shù)種能有效提升馬尾松純林土壤碳氮固持能力。馬尾松純林闊葉化改造短期內(nèi)對(duì)土壤氮儲(chǔ)量影響較小,土壤碳庫(kù)一定程度下降后,可能需要15 a以上的恢復(fù)期。

    關(guān)鍵詞:馬尾松人工純林;碳氮固存;補(bǔ)植闊葉樹(shù);短期效應(yīng)

    中圖分類(lèi)號(hào):S714.2 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1673-923X(2024)10-0126-12

    基金項(xiàng)目:貴州省科技計(jì)劃項(xiàng)目(黔科合服企〔2020〕4010);2023年貴州天然林保護(hù)管理補(bǔ)助資金項(xiàng)目“天然林資源保護(hù)與修復(fù)效益監(jiān)測(cè)”;2023年貴州森林資源管理補(bǔ)助資金項(xiàng)目“貴州人工商品純林樹(shù)種結(jié)構(gòu)優(yōu)化調(diào)整對(duì)森林碳匯效益提升監(jiān)測(cè)研究”;貴州雷公山森林生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站項(xiàng)目資助。

    The short-term effects of converting pure Pinus massoniana forests into mixed broadleaved forests on soil carbon and nitrogen sequestration

    WANG Haodong1,2,3,4, CHEN Meng1,2, YUAN Congjun1,2, HE Shuang1,2, DING Fangjun1,2, YANG Rui4

    (1. Key Laboratory of National Forestry and Grassland Administration on Biodiversity Conservation in Karst Mountainous Areas of Southwestern China, Guizhou Academy of Forestry, Guiyang 550005, Guizhou, China; 2. Guizhou Leigongshan Observation and Research Station for Forest Ecosystem, Leishan 557100, Guizhou, China; 3. Longli Forestry Farm of Guizhou Prorince, Longli 551200, Guizhou, China; 4. College of Forestry, Guizhou University, Guizhou 550025, Guiyang, China)

    Abstract:【Objective】In order to understand the impact of replanting different broad-leaved tree species on soil carbon storage in Pinus massoniana forests, a study was conducted to identify an effective carbon sequestration model for the broadleaf transformation of P. massoniana plantations.【Method】Within the state-owned forest farm in Dushan county, Guizhou province, a representative artificially planted pure forest of P. massoniana with similar forest structure was selected. After selective logging, Cinnamomum camphora, Phoebe zhennan, Choerospondias axillaris, and Liriodendron chinense of broad-leaved tree species were replanted in the same forest plot for 8-14 a. Three plots were designated for each condition, with a total of 15 plots, and soil samples of 0-20, 20-40, and 40-60 cm were collected for determination of soil physiochemical properties and active organic carbon components, using the pure forest of P. massoniana without replanting the broad-leaved trees as the control.【Result】Replanting different broad-leaved tree species had varying effects on soil carbon pools, with significant effects on the 0-20 cm soil layer and smaller effects on deeper soil layers. The range of soil carbon storage changes in the 0-60 cm layer was approximately 1.39-12.77 kg·m-2, with the highest soil carbon storage observed in pure forests of P. massoniana and P. massoniana + C. axillaris, reaching 212.35 and 203.51 kg·m-2, respectively, while the P. massoniana and C. camphora forest had the lowest soil carbon storage, approximately 100.78 kg·m-2. Replanting broad-leaved trees significantly reduced soil pH but increased microbial biomass carbon (MBC). In the 0-20 cm soil layer, the soil organic carbon and total nitrogen significantly increased in the P. massoniana + C. axillaris forest. Although soil carbon storage decreased in the short term after replanting broad-leaved trees, soil nitrogen storage did not significantly change. The replanting of P. zhennan significantly increased soil carbon pool activity and carbon pool activity index, while the replanting of C. axillaris significantly improved soil carbon pool stability.【Conclusion】Prioritizing the selection of deciduous broad-leaved tree species with high initial lignin /N values in the leaf effectively enhanced the soil carbon and nitrogen retention capacity of pure P. massoniana forests. The transformation of pure P. massoniana forests into broad-leaved forests had a minor short-term effect on soil nitrogen storage. After a certain degree of decline in soil carbon stocks, it might require a recovery period of fifteen years or more.

    Keywords: Pinus massoniana artificial pure forest; carbon and nitrogen sequestration; replanting broad-leaved trees; short-term effects

    近年來(lái),氣候變暖已成為人們必須重視的全球性問(wèn)題[1]??刂拼髿鉁厥覛怏w二氧化碳濃度最有效的途徑是人為減排和森林碳匯。森林固碳可分為地上和地下兩部分,地下土壤是陸地生態(tài)系統(tǒng)最大的碳庫(kù),在全球碳循環(huán)中起著至關(guān)重要的作用[2-3]。與植物碳庫(kù)相比,土壤碳庫(kù)能保存更久[4],對(duì)穩(wěn)定大氣二氧化碳濃度至關(guān)重要。即使土壤碳庫(kù)發(fā)生微小變化也會(huì)引起大氣二氧化碳含量的顯著變化[5],而氣候變暖反過(guò)來(lái)會(huì)加快土壤碳庫(kù)的分解[6],因此,土壤碳庫(kù)穩(wěn)定性備受關(guān)注[7]。土壤有機(jī)碳主要來(lái)源于凋落物和根系分泌物,因而土壤碳庫(kù)穩(wěn)定性與地上植被類(lèi)型密切相關(guān)。目前的研究按照周轉(zhuǎn)的快慢、轉(zhuǎn)化控制因素、分組方法的不同等將土壤有機(jī)碳分為不同的獨(dú)立組分,如活性有機(jī)碳、慢性有機(jī)碳、惰性有機(jī)碳、顆粒有機(jī)碳、可溶性有機(jī)碳、輕組有機(jī)碳、重組有機(jī)碳、黏粉粒有機(jī)碳、團(tuán)聚體保護(hù)有機(jī)碳等[8],共同參與土壤腐殖質(zhì)和團(tuán)聚體結(jié)構(gòu)的形成、土壤供水供肥過(guò)程[9-10]。土壤碳循環(huán)與氮存在一定關(guān)聯(lián),因?yàn)樵诮^大多數(shù)生態(tài)系統(tǒng)中氮與凈初級(jí)產(chǎn)量直接關(guān)聯(lián)[11],土壤氮?jiǎng)討B(tài)對(duì)土壤碳固存有較強(qiáng)影響,研究土壤碳庫(kù)時(shí)分析土壤氮固持能力不可或缺[12],二者共同參與地球化學(xué)循環(huán),對(duì)環(huán)境變化響應(yīng)迅速,備受關(guān)注[13]。

    針葉林連栽經(jīng)營(yíng)常引起土壤肥力下降和土壤酸化等問(wèn)題[14],為更好地解決這些問(wèn)題,需采取有目的性的森林經(jīng)營(yíng)措施。研究表明,改變森林類(lèi)型和林分結(jié)構(gòu)等森林經(jīng)營(yíng)措施影響土壤固持碳氮的作用[15-16]。通過(guò)調(diào)整樹(shù)種結(jié)構(gòu)將針葉林變?yōu)獒橀熁旖涣?,可加快凋落物的分解,也進(jìn)一步豐富森林物種多樣性,有利于土壤碳固存。當(dāng)前這一結(jié)論已得到大量研究的證實(shí)[17-19]。但是,林分結(jié)構(gòu)調(diào)整后短期內(nèi)對(duì)土壤碳氮固持的影響以及不同針闊混交林土壤碳儲(chǔ)量是否存在差異等問(wèn)題仍不十分清楚,還需要深入探索。

    馬尾松(Pinus massoniana,PM)作為中國(guó)亞熱帶主要造林樹(shù)種,在全國(guó)人工林中占據(jù)重要地位。目前的研究主要集中在馬尾松純林土壤理化性質(zhì)[20]和碳儲(chǔ)量[21]。盧立華等[19]研究發(fā)現(xiàn),與純林相比,馬尾松混交林土壤碳儲(chǔ)量顯著提升。徐芷君等[22]研究發(fā)現(xiàn),與木荷混交種植的馬尾松混交林同樣提升了土壤碳氮固持能力。相較于闊葉樹(shù)種,馬尾松針葉木質(zhì)素和樹(shù)脂含量高,比葉面積小[23],使得微生物分解速率較慢,分解后產(chǎn)生的酸性成分,抑制微生物的捕食和繁殖等活動(dòng)[24],導(dǎo)致馬尾松針葉分解緩慢。但是,較多研究?jī)H有一種馬尾松混交林,關(guān)于馬尾松與不同闊葉樹(shù)種混交林的對(duì)比研究較少,缺乏足夠的對(duì)比研究,可能會(huì)限制我們對(duì)森林生態(tài)系統(tǒng)土壤碳氮固持動(dòng)態(tài)變化的認(rèn)識(shí),還影響森林經(jīng)營(yíng)管理決策的科學(xué)性。因此,比較研究馬尾松闊葉化改造,有利于維持馬尾松林生態(tài)系統(tǒng)穩(wěn)定并提高固碳能力。

    基于此,本研究在樹(shù)種結(jié)構(gòu)調(diào)整措施的基礎(chǔ)上,在馬尾松林下補(bǔ)植香樟(Cinnamomum camphora,CC)、楠木(Phoebe zhennan,PZ)、南酸棗(Choerospondias axillaris,CA)和鵝掌楸(Liriodendron chinense,LC)。由于前期對(duì)馬尾松純林進(jìn)行了間伐,導(dǎo)致在試驗(yàn)初期人工林內(nèi)物種多樣性大幅下降,作出以下假設(shè):短期內(nèi)森林的碳儲(chǔ)量無(wú)法恢復(fù)到最初的狀態(tài),且這一過(guò)程中不同混交林土壤碳儲(chǔ)量差異顯著,這通常與生態(tài)系統(tǒng)植物多樣性和土壤養(yǎng)分有關(guān)。本研究通過(guò)對(duì)比分析不同樹(shù)種搭配模式下,土壤碳氮固持能力的變化趨勢(shì)以及差異性,評(píng)估短期內(nèi)馬尾松人工林闊葉化改造后土壤碳氮固持效應(yīng)。研究馬尾松純林闊葉化改造后土壤碳氮固持能力,可為針葉林改造調(diào)整提供理論依據(jù),有助于篩選馬尾松純林闊葉化改造土壤高效固碳模式。

    1 研究區(qū)概況與方法

    1.1 研究區(qū)概況

    研究區(qū)位于貴州省黔南布衣族苗族自治州東南部獨(dú)山縣國(guó)有林場(chǎng)(圖1),獨(dú)山縣位于貴州省,地處貴州最南端,素有“貴州南大門(mén)”“西南門(mén)戶(hù)”之稱(chēng),是“中國(guó)花燈藝術(shù)之鄉(xiāng)”。全縣平均海拔850~1 100 m,屬中亞熱帶溫潤(rùn)季風(fēng)性氣候,四季分明,冬無(wú)嚴(yán)寒、夏無(wú)酷暑,年平均氣溫15 ℃,年均降水量1 430 mm,無(wú)霜期297 d。獨(dú)山縣國(guó)有林場(chǎng)占地面積約18 588.00 hm2,其中:林地15 839.14 hm2,占總面積的85.21%,非林地2 748.86 hm2,占總面積的 14.79%。森林面積 13 564.32 hm2,森林覆蓋率 72.96%,森林總蓄積865 465.48 m3,是貴州最大的國(guó)有林場(chǎng)。林場(chǎng)內(nèi)出露的巖石以砂頁(yè)巖和碳酸巖為主,土壤以硅鋁質(zhì)和鐵鋁質(zhì)黃壤為主,土壤厚度一般在中層以上,肥力中等,多呈酸性或微酸性。林場(chǎng)栽植的樹(shù)種以馬尾松、濕地松Pinus elliottii、香樟和鵝掌楸等為主。

    1.2 樣地設(shè)置與調(diào)查

    本研究土壤樣品采自4種不同的馬尾松-闊葉混交林(馬尾松+香樟、馬尾松+楠木、馬尾松+南酸棗和馬尾松+鵝掌楸)和未改造的馬尾松純林。選擇林分結(jié)構(gòu)相近且具有代表性的馬尾松純林,通過(guò)林分結(jié)構(gòu)調(diào)整,將其改造成與香樟、楠木、南酸棗和鵝掌楸混交的森林。每種類(lèi)型設(shè)置3塊樣地,樣地面積為667 m2,共15塊樣地。調(diào)查固定樣地基本信息(包括海拔、坡度等),對(duì)樣地內(nèi)所有胸徑大于5 cm的喬木和補(bǔ)植的闊葉樹(shù)種進(jìn)行每木檢尺,結(jié)果如表1所示。在每個(gè)樣地的三個(gè)角分別布設(shè)1 m2的草本樣方,于草本樣方內(nèi)收集全部的凋落物,稱(chēng)量樣方內(nèi)凋落物的鮮質(zhì)量,并帶回實(shí)驗(yàn)室備用。馬尾松純林在文中縮寫(xiě)為PP,馬尾松+香樟混交林縮寫(xiě)為PCC,馬尾松+楠木混交林縮寫(xiě)為PPZ,馬尾松+南酸棗混交林縮寫(xiě)為PCA,馬尾松+鵝掌楸混交林縮寫(xiě)為PLC。

    1.3 土壤樣品采集與指標(biāo)測(cè)定

    于2023年5月29日采集土壤樣品,每個(gè)樣地挖開(kāi)一個(gè)土壤剖面,依次從0~20、20~40和40~60 cm采集不同土層的土壤,并用環(huán)刀挖取土壤用于測(cè)定土壤容重,挖出的土壤去除其中包含的根和石塊等雜物,完全風(fēng)干后,分別研磨過(guò)2 mm和0.149 mm篩并裝入密封袋中保存?zhèn)溆茫糜谕寥览砘再|(zhì)的測(cè)定,共45個(gè)樣品。其中,土壤有機(jī)碳含量采用重鉻酸鉀氧化-分光光度法測(cè)定,土壤全氮采用凱氏法測(cè)定,土壤全磷采用酸溶劑熱處理法測(cè)定,土壤pH值采用電位法測(cè)定[25],易氧化有機(jī)碳(ROC)采用高錳酸鉀氧化法測(cè)定[26],顆粒有機(jī)碳(POC)采用濕篩法-重鉻酸鉀外加熱法測(cè)定[27],可溶性有機(jī)碳(DOC)采用蒸餾水浸提法測(cè)定[28],微生物量碳(MBC)采用氯仿熏蒸法測(cè)定[29]。凋落物帶回后立即放入80 ℃恒溫烘箱中烘干,稱(chēng)其干質(zhì)量,通過(guò)計(jì)算得到每公頃的凋落物質(zhì)量。

    1.4 統(tǒng)計(jì)分析

    數(shù)據(jù)統(tǒng)計(jì)與分析在Excel 2016和SPSS 26.0軟件中進(jìn)行,采用單因素方差分析檢驗(yàn)馬尾松不同改造模式下的差異顯著性(P<0.05)。在R 4.1.3中運(yùn)用“GGally”程序包繪制相關(guān)性矩陣圖,運(yùn)用“Vegan”包計(jì)算Shannon指數(shù)和Rarefied SR指數(shù),RDA分析在Canoco 5.0軟件中繪制,其余圖片運(yùn)用Origin 2021軟件繪制。

    2 結(jié)果與分析

    2.1 不同改造模式群落植物多樣性與凋落物量特征

    由表2可知,馬尾松+香樟混交林林下植物多樣性顯著高于馬尾松純林(P<0.05),馬尾松純林和馬尾松+鵝掌楸混交林凋落物量顯著高于其他處理(P<0.05)。

    2.2 不同改造模式土壤理化性質(zhì)

    如表3所示,不同土層馬尾松純林土壤pH值均為最高,且補(bǔ)植闊葉樹(shù)種土壤pH值與馬尾松純林之間差異顯著(P<0.05)。0~20 cm土層中,馬尾松+南酸棗林土壤容重顯著低于馬尾松純林,馬尾松+南酸棗林土壤有機(jī)碳含量最高,馬尾松+南酸棗林和馬尾松+鵝掌楸林土壤全氮顯著高于馬尾松純林(P<0.05);而20~40 cm土壤中,僅馬尾松+南酸棗林和馬尾松+鵝掌楸林土壤有機(jī)碳含量差異顯著(P<0.05),在40~60 cm土壤中,不同改造模式下土壤有機(jī)碳含量均無(wú)顯著差異(P>0.05)。

    2.3 不同改造模式土壤碳氮固存特征

    結(jié)果表明,土壤碳固存受到不同混交模式的影響,氮固存受到的影響較?。▓D2)。在0~60 cm土層中,馬尾松純林碳儲(chǔ)量最高,且顯著高于馬尾松+楠木混交林(P<0.05),其次是馬尾松+南酸棗混交林,其余處理均較低,馬尾松純林與其余處理無(wú)顯著差異(P>0.05)。對(duì)于0~20 cm土層,除馬尾松+南酸棗混交林外,馬尾松純林與其他預(yù)處理差異達(dá)顯著水平(P<0.05)。對(duì)于20~40和40~60 cm土層,各處理間差異較小(P>0.05);土壤氮固存在各處理間均無(wú)顯著差異(P>0.05)。

    2.4 不同改造模式土壤碳組分分布特征

    結(jié)果顯示(圖3),在0~20、20~40和40~60 cm土層中,馬尾松純林混交林土壤微生物量碳顯著低于其余處理(P<0.05),顆粒有機(jī)碳在各處理間無(wú)顯著差異(P>0.05),在0~20和40~60 cm土層中,馬尾松+鵝掌楸混交林顯著高于馬尾松+南酸棗林,在20~40 cm土層,各處理無(wú)顯著差異(P>0.05),不同混交林的土壤易氧化有機(jī)碳僅表層的馬尾松+南酸棗混交林和馬尾松+香樟混交林差異顯著(P<0.05),其余處理差異不顯著(P>0.05)。

    2.5 不同改造模式土壤碳庫(kù)管理指數(shù)特征

    如表4所示,混入闊葉樹(shù)種后,0~20 cm僅混入楠木的碳庫(kù)活度和碳庫(kù)活度指數(shù)得到提升,碳庫(kù)指數(shù)僅混入南酸棗后提高,增長(zhǎng)了256%,混入闊葉樹(shù)種后0~40 cm的碳庫(kù)管理指數(shù)只有楠木和南酸棗顯著提升(P<0.05),其余指標(biāo)在更深土層無(wú)顯著差異(P>0.05)。

    2.6 土壤理化性質(zhì)與土壤碳氮儲(chǔ)存的關(guān)系

    土壤碳氮儲(chǔ)量和碳組分相關(guān)性分析如圖4所示。相關(guān)性表現(xiàn)為,碳儲(chǔ)量與ROC和POC極顯著正相關(guān)(P<0.001),POC和ROC存在極顯著的兩兩正相關(guān)關(guān)系(P<0.001),DOC與MBC顯著正相關(guān)(P<0.05)。

    冗余分析結(jié)果(圖5)顯示,環(huán)境因子共解釋了土壤碳氮儲(chǔ)量變異的93.36%,前兩軸分別解釋了93.19%和0.17%(圖5a)。SOC顯著影響土壤氮儲(chǔ)量(P<0.01),TN顯著影響土壤碳儲(chǔ)量(P<0.01),土壤碳、氮儲(chǔ)量存在正相關(guān)關(guān)系,但并未達(dá)到顯著水平,植物多樣性抑制了土壤碳儲(chǔ)量的增長(zhǎng),土壤全磷含量和調(diào)落物量顯著促進(jìn)了土壤碳氮固存(P<0.05)(圖5b)。

    3 討 論

    3.1 不同改造模式對(duì)土壤碳庫(kù)垂直分布的影響

    植物生長(zhǎng)發(fā)育受到土壤有機(jī)碳的轉(zhuǎn)化與礦質(zhì)元素釋放的影響[30],植物根系活動(dòng)影響土壤環(huán)境[31],二者相互作用、相互影響,共同維持生態(tài)系統(tǒng)穩(wěn)定。由于有機(jī)碳的來(lái)源和組成以及環(huán)境條件的差異,底土中的碳可能比表土中的碳更穩(wěn)定[32],僅僅關(guān)注表層土壤無(wú)法揭示整個(gè)土壤剖面的土壤碳特征。研究表明,土壤養(yǎng)分含量隨土層的深入逐漸減少[33-34]。本研究中,不同混交模式下,土壤SOC、TN和TP隨土層的增加逐漸減少,或隨土層的增加先增加再減少。造成這一現(xiàn)象的原因可能是由于植物根系主要集中在土壤表層,根系活動(dòng)改變了土壤結(jié)構(gòu),改善了土壤環(huán)境[35]。同時(shí),土壤容重隨土層的增加逐漸變大,深層土壤由于受更強(qiáng)的擠壓作用不利于植物根系和土壤微生物的呼吸。加上有機(jī)質(zhì)和地下生物量的減少,導(dǎo)致土壤碳、氮儲(chǔ)量和碳組分呈土層加深而遞減的趨勢(shì)[36]。結(jié)合本研究結(jié)果來(lái)看,補(bǔ)植不同樹(shù)種后土壤碳氮儲(chǔ)量隨著土層由淺至深的變化趨勢(shì)不完全一致,暗示了不同種類(lèi)的闊葉樹(shù)對(duì)土壤碳轉(zhuǎn)化和循環(huán)的影響是不同的,這可能與不同植物的根系垂直分布格局和光合產(chǎn)物的分配差異有關(guān)[37-38]。本研究?jī)H研究了馬尾松純林闊葉改造后短期內(nèi)的土壤碳氮儲(chǔ)量變化,并未研究植物根系的分布特征以及樹(shù)種對(duì)資源的吸收、利用的能力。未來(lái)的研究可以著重于探討這些因素,并通過(guò)綜合分析以揭示更深層次的機(jī)理。

    3.2 不同改造模式對(duì)土壤碳氮固持的影響

    MBC與整個(gè)土壤生態(tài)系統(tǒng)養(yǎng)分周轉(zhuǎn)密切相關(guān),其大小主要與土壤微生物的生存環(huán)境有關(guān)[39]。本研究中,補(bǔ)植闊葉樹(shù)種使土壤表層MBC得到顯著提升,說(shuō)明補(bǔ)植闊葉樹(shù)種后微生物的生存環(huán)境得到優(yōu)化。ROC是主要的活性有機(jī)碳,反映土壤碳庫(kù)的活躍度[40]。結(jié)果顯示,ROC和SOC成顯著正相關(guān),表征土壤ROC對(duì)土壤碳庫(kù)變化的強(qiáng)敏感性[41]。研究結(jié)果表明,DOC作為土壤最具動(dòng)態(tài)特征的碳組分,與MBC常有較強(qiáng)的相關(guān)性[42],這與本研究結(jié)果一致,二者對(duì)環(huán)境和氣候變化的敏感性較強(qiáng),都可以作為評(píng)價(jià)土壤微生物分解利用土壤速效養(yǎng)分的重要指標(biāo)。

    林分結(jié)構(gòu)調(diào)整明顯改變森林群落的物種組成和群落結(jié)構(gòu),植物多樣性將受到影響[43]。Shannon指數(shù)和Rarefied SR反映的是在相同樣本數(shù)量或大小的情況下,物種出現(xiàn)的豐富程度,通常與生物量存在正相關(guān)關(guān)系[44],且對(duì)土壤碳固存具有促進(jìn)作用[18,45]。本研究中,植物多樣性抑制了馬尾松林和混交林的土壤碳固存,這與前人研究結(jié)果不一致[46]。原因可能是前期的森林經(jīng)營(yíng)措施大大降低了森林群落總體生物量和物種多樣性,短期內(nèi)難以恢復(fù)[47],而馬尾松林林分密度的降低將嚴(yán)重影響土壤碳儲(chǔ)量[48],導(dǎo)致土壤碳儲(chǔ)量甚至低于未開(kāi)展森林經(jīng)營(yíng)的林分。這完全符合前期做出的假設(shè),但這并不意味著采取的森林經(jīng)營(yíng)措施是不合適的。從結(jié)果中可看出,只有混入香樟的林分土壤碳儲(chǔ)量顯著低于馬尾松純林,其余林分土壤碳儲(chǔ)量與馬尾松純林差異不顯著,且年限較久的南酸棗+馬尾松林土壤碳儲(chǔ)量即將超過(guò)馬尾松純林。由表1可知,南酸棗雖年限較長(zhǎng),但其增長(zhǎng)速度明顯較其他樹(shù)種更快,即在有限的生長(zhǎng)周期內(nèi),生長(zhǎng)越快的闊葉樹(shù)種越能提升人工針葉林的土壤碳匯能力,且林分結(jié)構(gòu)調(diào)整后可能需要15 a以上才能恢復(fù)至原來(lái)的狀態(tài)。香樟生長(zhǎng)周期較其他樹(shù)種慢,相較于同期補(bǔ)植的其他闊葉樹(shù)種,其體型較小,而南酸棗生長(zhǎng)迅速,一方面更多容易分解的凋落物為土壤微生物提供了更多的分解底物[49],另一方面,合理的林分郁閉度避免了太陽(yáng)長(zhǎng)期照射抑制土壤微生物活力[50],這在一定程度上加快了土壤碳固存。郭婧等[51]對(duì)比馬尾松+石櫟針闊混交林、南酸棗落葉闊葉林、石櫟+青岡常綠闊葉林3種次生林和杉木人工純林的凋落物量與周轉(zhuǎn)期,發(fā)現(xiàn)次生林的凋落物量顯著高于杉木人工林,且落葉樹(shù)種南酸棗的凋落物分解效率相較于常綠闊葉林更高,周轉(zhuǎn)期更短。說(shuō)明樹(shù)種特性所決定的凋落物量及其分解速率影響了森林生態(tài)系統(tǒng)的養(yǎng)分循環(huán)功能。森林生態(tài)系統(tǒng)養(yǎng)分歸還的能力與土壤碳氮固持能力密切關(guān)聯(lián),不同樹(shù)種搭配不僅增加了森林生態(tài)系統(tǒng)物種多樣性,提升了林分的穩(wěn)定性與抗逆能力。從凋落物的分解角度來(lái)看,某些凋落物混合分解后通常都能加快分解速率,這些過(guò)程可促進(jìn)土壤對(duì)碳和氮的固定[52]。綜合來(lái)看,補(bǔ)植闊葉樹(shù)種有利于土壤對(duì)碳氮的固持,但樹(shù)種的配置模式是其中的關(guān)鍵,通常認(rèn)為選擇初始木質(zhì)素/N值較高的物種分解效率更高[53],選擇補(bǔ)植此類(lèi)樹(shù)種更契合經(jīng)營(yíng)目的的闊葉樹(shù)種可有效提升馬尾松純林土壤碳固存。陳歆宇等[54]研究發(fā)現(xiàn)混交比例和林分密度對(duì)馬尾松混交林的生產(chǎn)功能有較大的影響。因此,未來(lái)需要在確定混交樹(shù)種后,探究馬尾松混交林最佳混交比例和最佳林分密度。

    3.3 不同改造模式對(duì)土壤碳庫(kù)管理指數(shù)的影響

    土壤碳庫(kù)管理指數(shù)可以表征土地利用或管理方式對(duì)土壤碳庫(kù)質(zhì)量提升的能力[55]。本研究中,對(duì)于0~20 cm土層,僅補(bǔ)植楠木的碳庫(kù)活度和碳庫(kù)活度指數(shù)得到提升,說(shuō)明馬尾松純林補(bǔ)植闊葉樹(shù)種初期,補(bǔ)植除楠木外的其余樹(shù)種均會(huì)降低土壤有機(jī)碳的分解速率。深層土壤(20~40 cm)的碳庫(kù)活度和碳庫(kù)活度指數(shù)較接近,可能是由于該區(qū)域受植物根系影響較小,土壤環(huán)境差異較小[35]。碳庫(kù)管理指數(shù)是土壤碳儲(chǔ)量與土壤總碳中不穩(wěn)定碳活度指標(biāo)的最終平衡結(jié)果,直接反映土壤碳庫(kù)穩(wěn)定性。本研究?jī)H碳庫(kù)活度指數(shù)僅0~20 cm土層在補(bǔ)植南酸棗后提高,說(shuō)明補(bǔ)植南酸棗后土壤質(zhì)量得到了提升。土壤碳庫(kù)管理指數(shù)可作為評(píng)估最佳補(bǔ)闊樹(shù)種的指標(biāo)之一。

    4 結(jié) 論

    本研究基于馬尾松人工純林,通過(guò)林分結(jié)構(gòu)調(diào)整混入不同闊葉樹(shù)種,以期找出馬尾松高效固碳模式。結(jié)果顯示,補(bǔ)植不同闊葉樹(shù)種后對(duì)馬尾松林土壤碳庫(kù)影響不一致,0~20 cm土壤受影響較強(qiáng),20~60 cm受到影響較弱;補(bǔ)植闊葉樹(shù)種后土壤pH值顯著降低;對(duì)于0~20 cm土壤,南酸棗+馬尾松混交林土壤有機(jī)碳和全氮顯著提升;總體來(lái)看,經(jīng)林分結(jié)構(gòu)調(diào)整后的馬尾松林分,補(bǔ)植闊葉樹(shù)種在短期內(nèi)無(wú)法提升土壤碳儲(chǔ)量,但未來(lái)具有較好的提升潛力;評(píng)估補(bǔ)植闊葉樹(shù)種對(duì)土壤碳儲(chǔ)量和碳庫(kù)穩(wěn)定性的影響,發(fā)現(xiàn)南酸棗在目前研究的幾個(gè)闊葉樹(shù)種中表現(xiàn)最佳,馬尾松闊葉化改造后的土壤碳庫(kù)可能需要15 a以上的恢復(fù)期,在馬尾松純林中補(bǔ)植初始木質(zhì)素/N值較高的闊葉樹(shù)種更有利于土壤碳氮的固持。

    參考文獻(xiàn):

    [1] NOORMETS A, EPRON D, DOMEC J, et al. Effects of forest management on productivity and carbon sequestration: a review and hypothesis[J]. Forest Ecology and Management,2015,355: 124-140.

    [2] LEI W Y, PAN Q, TENG P J, et al. How does soil organic matter stabilize with soil and environmental variables along a black soil belt in northeast China? An explanation using FTIR spectroscopy data[J]. Catena,2023,228:107152.

    [3] LANGE M, EISENHAUER N, SIERRA C A, et al. Plant diversity increases soil microbial activity and soil carbon storage[J]. Nature Communications,2015,6(1):6707.

    [4] RATTAN L. Sequestering carbon and increasing productivity by conservation agriculture[J]. Journal of Soil and Water Conservation,2015,70(3):55A.

    [5] ZHOU X H, XU X, ZHOU G Y, et al. Temperature sensitivity of soil organic carbon decomposition increased with mean carbon residence time: field incubation and data assimilation[J]. Global Change Biology,2018,24(2):810-822.

    [6] YOU M Y, HE P, DAI S S, et al. Priming effect of stable C pool in soil and its temperature sensitivity[J]. Geoderma,2021, 401:115216.

    [7] PENG Y Y, THOMAS S C, TIAN D L. Forest management and soil respiration: Implications for carbon sequestration[J]. Environmental Reviews,2008,16:93-111.

    [8] 趙鑫,宇萬(wàn)太,李建東.不同經(jīng)營(yíng)管理?xiàng)l件下土壤有機(jī)碳及其組分研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2006,17(11):2203-2209. ZHAO X, YU W T, LI J D. Research advances in soil organic carbon and its fractions under different management patterns[J]. Chinese Journal of Applied Ecology,2006,17(11):2203-2209.

    [9] 張維理, KOLBE H,張認(rèn)連.土壤有機(jī)碳作用及轉(zhuǎn)化機(jī)制研究進(jìn)展[J].中國(guó)農(nóng)業(yè)科學(xué),2020,53(2):317-331. ZHANG W L, KOLBE H, ZHANG R L. Research progress of SOC functions and transformation mechanisms[J]. Scientia Agricultura Sinica,2020,53(2):317-331.

    [10] LAVALLEE J M, SOONG J L, COTRUFO M F. Conceptualizing soil organic matter into particulate and mineral-associated forms to address global change in the 21st century[J]. Global Change Biology,2020,26(1):261-273.

    [11] LEBAUER D S, TRESEDER K K. Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed[J]. Ecology,2008,89(2):371-379.

    [12] JELINSKI N A, KUCHARIK C J. Land-use effects on soil carbon and nitrogen on a US Midwestern floodplain[J]. Soil Science Society of America Journal,2009,73(1):217-225.

    [13] LAN J C, HU N, FU W L. Soil carbon-nitrogen coupled accumulation following the natural vegetation restoration of abandoned farmlands in a karst rocky desertification region[J]. Ecological Engineering,2020,158:106033.

    [14] DONG X D, GAO P, ZHOU R, et al. Changing characteristics and influencing factors of the soil microbial community during litter decomposition in a mixed Quercus acutissima Carruth. and Robinia pseudoacacia L. forest in northern China[J]. Catena, 2021,196:104811.

    [15] XIONG Y M, XIA H P, LI Z A, et al. Impacts of litter and understory removal on soil properties in a subtropical Acacia mangium plantation in China[J]. Plant and Soil,2008,304:179-188.

    [16] ZHAO J, WANG X L, SHAO Y H, et al. Effects of vegetation removal on soil properties and decomposer organisms[J]. Soil Biology and Biochemistry,2011,43(5):954-960.

    [17] BAI Y X, ZHOU Y C, CHEN X L, et al. Tree species composition alters the decomposition of mixed litter and the associated microbial community composition and function in subtropical plantations in China[J]. Forest Ecology and Management,2023, 529:120743.

    [18] CHEN X L, TAYLOR A R, REICH P B, et al. Tree diversity increases decadal forest soil carbon and nitrogen accrual[J]. Nature,2023,618:94-101.

    [19] 盧立華,郭文福,蔡道雄,等.馬尾松與紅椎純林及混交林生態(tài)系統(tǒng)碳儲(chǔ)量研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2019,39(7): 78-84. LU L H, GUO W F, CAI D X, et al. Study on carbon storage of monoculture and mixed plantation of Pinus massoniana and Castanopsis hystrix[J]. Journal of Central South University of Forestry Technology,2019,39(7):78-84.

    [20] 陳蓉,王韋韋,曹麗榮,等.馬尾松和杉木人工林細(xì)根碳氮磷化學(xué)計(jì)量特征隨土層深度的變化[J].生態(tài)學(xué)報(bào),2023,43(9): 3709-3718. CHEN R, WANG W W, CAO L R, et al. Variation of carbon, nitrogen and phosphorus stoichiometric characteristics of fine roots in masson pine and Chinese fir plantations with soil depth[J]. Acta Ecologica Sinica,2023,43(9):3709-3718.

    [21] 陶玉華,馮金朝,馬麟英,等.廣西羅城馬尾松、杉木、桉樹(shù)人工林碳儲(chǔ)量及其動(dòng)態(tài)變化[J].生態(tài)環(huán)境學(xué)報(bào),2011,20(11): 1608-1613. TAO Y H, FENG J C, MA L Y, et al. Carbon storage and distribution of massion pine, Chinese fir and eucalyptus plantations at Liuzhou, Guangxi[J]. Ecology and Environmental Sciences,2011,20(11):1608-1613.

    [22] 徐芷君,劉苑秋,方向民,等.亞熱帶2種針葉林土壤碳氮磷儲(chǔ)量及化學(xué)計(jì)量比對(duì)混交的響應(yīng)[J].水土保持學(xué)報(bào),2019, 33(1):165-170. XU Z J, LIU W Q, FANG X M, et al. The responses of soil carbon, nitrogen and phosphorus storage and their stoichiometry in two coniferous forests to mixed effect in subtropical area[J]. Journal of Soil and Water Conservation,2019,33(1):165-170.

    [23] 劉濤,孫守琴,邱陽(yáng).川西亞高山生態(tài)系統(tǒng)三種典型植物凋落物分解動(dòng)態(tài)特征[J].山地學(xué)報(bào),2017,35(5):663-668. LIU T, SUN S Q, QIU Y. Dynamics and differences in the decomposition of litters from three dominating plants in subalpine ecosystems in western Sichuan, China[J]. Mountain Research, 2017,35(5):663-668.

    [24] ZHANG Y, LI X, ZHANG D, et al. Characteristics of fungal community structure during the decomposition of mixed foliage litter from Pinus massoniana and broadleaved tree species in southwestern China[J]. Journal of Plant Ecology,2020,13(5): 574-588.

    [25] 鮑士旦.土壤農(nóng)化分析[M].北京:中國(guó)農(nóng)業(yè)出版社,2000. BAO S D. Soil and agrochemistry analytical methods[M]. Beijing: Chinese Agriculture Press,2000.

    [26] LEFROY R D B, BLAIR G J, STRONG W M. Changes in soil organic matter with cropping as measured by organic carbon fractions and 13C natural isotope abundance[J]. Plant and Soil,1993,155(1):399-402.

    [27] CAMBARDELLA C A, ELLIOTT E T. Particulate soil organic- matter changes across a grassland cultivation sequence[J]. Soil Science Society of America Journal,1992,56(3):777-783.

    [28] 呂國(guó)紅,周廣勝,周莉,等.土壤溶解性有機(jī)碳測(cè)定方法與應(yīng)用[J].氣象與環(huán)境學(xué)報(bào),2006,22(2):51-55. LYU G H, ZHOU G S, ZHOU L, et al. Methods of soil dissolved organic carbon measurement and their applications[J]. Journal of Meteorology and Environment,2006,22(2):51-55.

    [29] JENKINSON D S, POWLSON D S. The effects of biocidal treatments on metabolism in soil-V: a method for measuring soil biomass[J]. Soil Biology and Biochemistry,1976,8(3):209-213.

    [30] ZHANG C, LIU G B, XUE S, et al. Soil bacterial community dynamics reflect changes in plant community and soil properties during the secondary succession of abandoned farmland in the loess plateau[J]. Soil Biology and Biochemistry,2016,97:40-49.

    [31] WANG G L, LIU G B, XU M X. Above-and belowground dynamics of plant community succession following abandonment of farmland on the loess plateau, China[J]. Plant and Soil,2009,316: 227-239.

    [32] RUMPEL C, K?GEL-KNABNER I. Deep soil organic matter: a key but poorly understood component of terrestrial C cycle[J]. Plant and Soil,2011,338:143-158.

    [33] 孫雙紅,陳立新,李少博,等.闊葉紅松林不同演替階段土壤酶活性與養(yǎng)分特征及其相關(guān)性[J].北京林業(yè)大學(xué)學(xué)報(bào), 2016,38(2):20-28. SUN S H, CHEN L X, LI S B, et al. Characteristics of soil enzyme activity and nutrient content and their correlations at different succession stages of broadleaf-Korean pine forest[J]. Journal of Beijing Forestry University,2016,38(2):20-28.

    [34] PAN J W, GUO Q Q, LI H E, et al. Dynamics of soil nutrients, microbial community structure, enzymatic activity, and their relationships along a chronosequence of Pinus massoniana plantations[J]. Forests,2021,12(3):376.

    [35] 宋日,吳春勝,郭繼勛.東北草原植物殘?bào)w腐解動(dòng)態(tài)研究(簡(jiǎn)報(bào))[J].草業(yè)學(xué)報(bào),2002,11(2):105-108. SONG R, WU C S, GUO J X. Decomposition dynamics of plant residues in natural meadow in northeast China[J]. Acta Prataculturae Sinica,2002,11(2):105-108.

    [36] HU P, LIU S, YE Y, et al. Soil carbon and nitrogen accumulation following agricultural abandonment in a subtropical karst region[J]. Applied Soil Ecology,2018,132:169-178.

    [37] 趙巧巧,趙筱青,黃佩,等.云南亞熱帶地區(qū)主要林地類(lèi)型土壤碳含量變化及影響因素研究[J].林業(yè)科學(xué)研究,2024,37(1): 63-72. ZHAO Q Q, ZHAO X Q, HUANG P, et al. Soil carbon changes and its influencing factors in major forest types in the subtropical area of Yunnan province[J]. Forest Research,2024,37(1):63-72.

    [38] 唐雪婭,徐明,文春玉,等.黔中地區(qū)馬尾松針闊混交林細(xì)根性狀空間變異特征[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2023,43(11): 142-150. TANG X Y, XU M, WEN C Y, et al. Spatial variation of root functional traits of Pinus massoniana mixed forest in central Guizhou[J]. Journal of Central South University of Forestry Technology,2023,43(11):142-150.

    [39] DEVI N B, YADAVA P S. Seasonal dynamics in soil microbial biomass C, N and P in a mixed-oak forest ecosystem of Manipur, northeast India[J]. Applied Soil Ecology,2006,31(3):220-227.

    [40] 沈宏,曹志洪,胡正義.土壤活性有機(jī)碳的表征及其生態(tài)效應(yīng)[J].生態(tài)學(xué)雜志,1999,18(3):33-39. SHEN H, CAO Z H, HU Z Y. Characteristics and ecological effects of the active organic carbon in soil[J]. Chinese Journal of Ecology,1999,18(3):33-39.

    [41] 王義祥,翁伯琦,邢世和,等.果園土壤有機(jī)碳及其影響因素的研究進(jìn)展[J].福建農(nóng)業(yè)學(xué)報(bào),2011,26(6):1113-1122. WANG Y X, WENG B Q, XING S H, et al. Advance in soil organic carbon stock and the impact factors on orchard ecosystem research[J]. Fujian Journal of Agricultural Sciences,2011,26(6): 1113-1122.

    [42] CRESSEY E L, DUNGAIT J, JONES D L, et al. Soil microbial populations in deep floodplain soils are adapted to infrequent but regular carbon substrate addition[J]. Soil Biology and Biochemistry,2018,122:60-70.

    [43] 馮廣,李俊清,臧潤(rùn)國(guó),等.皆伐與刀耕火種后常綠-落葉闊葉混交林的動(dòng)態(tài)恢復(fù)機(jī)制[J].北京林業(yè)大學(xué)學(xué)報(bào),2019,41(10): 1-10. FENG G, LI J Q, ZANG R G, et al. Dynamics and mechanisms of natural restoration of evergreen-deciduous broadleaved mixed forest following clear cutting and slash-and-burn[J]. Journal of Beijing Forestry University,2019,41(10):1-10.

    [44] 王麗紅,辛穎,鄒夢(mèng)玲,等.大興安嶺火燒跡地植被恢復(fù)中植物多樣性與生物量分配格局[J].北京林業(yè)大學(xué)學(xué)報(bào), 2015,37(12):41-47. WANG L H, XIN Y, ZOU M L, et al. Plants diversity and biomass distribution of vegetation restoration in burned area of Great Xing’an mountains[J]. Journal of Beijing Forestry University,2015,37(12):41-47.

    [45] LANGE M, EISENHAUER N, SIERRA C A, et al. Plant diversity increases soil microbial activity and soil carbon storage[J]. Nature Communications,2015,6(1):6707.

    [46] CONG W F, VAN RUIJVEN J, MOMMER L, et al. Plant species richness promotes soil carbon and nitrogen stocks in grasslands without legumes[J]. Journal of Ecology,2014,102(5):1163-1170.

    [47] 范春楠,劉強(qiáng),鄭金萍,等.采伐強(qiáng)度對(duì)闊葉紅松林生態(tài)系統(tǒng)碳密度恢復(fù)的影響[J].北京林業(yè)大學(xué)學(xué)報(bào),2022,44(10):33-42. FAN C N, LIU Q, ZHENG J P, et al. Effects of logging intensity on restoration of carbon density in broadleaved Korean pine forest ecosystem[J]. Journal of Beijing Forestry University,2022, 44(10):33-42.

    [48] 李玉鳳,李妹珍,馬姜明,等.林分密度對(duì)馬尾松人工林土壤碳儲(chǔ)量及其分配特征的影響[J].廣西林業(yè)科學(xué),2021,50(1): 54-59. LI Y F, LI M Z, MA J P, et al. Effects of stand density on soil carbon storage and distribution characteristics of Pinus massoniana plantations[J]. Guangxi Forestry Science,2021,50(1): 54-59.

    [49] BALAMI S, VA?UTOVá M, KO?NAR J, et al. Soil fungal communities in abandoned agricultural land has not yet moved towards the seminatural forest[J]. Forest Ecology and Management,2021,491:119181.

    [50] LEI W, PAN Q, TENG P, et al. How does soil organic matter stabilize with soil and environmental variables along a black soil belt in northeast China? An explanation using FTIR spectroscopy data[J]. Catena,2023,228:107152.

    [51] 郭婧,喻林華,方晰,等.中亞熱帶4種森林凋落物量、組成、動(dòng)態(tài)及其周轉(zhuǎn)期[J].生態(tài)學(xué)報(bào),2015,35(14):4668-4677. GUO J, YU L H, FANG X, et al. Litter production and turnover in four types of subtropical forests in China[J]. Acta Ecologica Sinica, 2015,35(14):4668-4677.

    [52] ZHAO H, YANG R, YUAN C, et al. Chemical stoichiometry and enzyme activity changes during mixed decomposition of Camellia sinensis pruning residues and companion tree species litter[J]. Agronomy,2023,13(7):1717.

    [53] 劉莎茜,楊瑞,侯春蘭,等.貴州山區(qū)生態(tài)茶園不同凋落物木質(zhì)素、纖維素分解特征[J].茶葉科學(xué),2021,41(5):654-668. LIU S Q, YANG R, HOU C L, et al. Decomposition characteristics of lignin and cellulose in different litters of ecological tea gardens in mountainous areas of Guizhou[J]. Journal of Tea Science,2021,41(5):654-668.

    [54] 陳歆宇,譚偉,楊深鈞,等.不同類(lèi)型馬尾松混交林結(jié)構(gòu)與生產(chǎn)功能的耦合關(guān)系[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2023,43(4): 53-61. CHEN X Y, TAN W, YANG S J, et al. The coupling relationship between structure and production function of different types of Pinus massoniana mixed forest[J]. Journal of Central South University of Forestry Technology,2023,43(4):53-61.

    [55] MEENA V S, MONDAL T, PANDEY B M, et al. Land use changes: strategies to improve soil carbon and nitrogen storage pattern in the mid-Himalaya ecosystem, India[J]. Geoderma, 2018,321:69-78.

    [本文編校:吳 彬]

    亚洲成a人片在线一区二区| 中文亚洲av片在线观看爽| 99视频精品全部免费 在线| 99在线人妻在线中文字幕| 19禁男女啪啪无遮挡网站| 叶爱在线成人免费视频播放| 女人十人毛片免费观看3o分钟| 国产亚洲av嫩草精品影院| 99久久综合精品五月天人人| 日本在线视频免费播放| 三级毛片av免费| 久久精品国产亚洲av香蕉五月| 美女大奶头视频| 亚洲专区国产一区二区| 日韩中文字幕欧美一区二区| 97碰自拍视频| 亚洲av成人不卡在线观看播放网| 搡女人真爽免费视频火全软件 | 日韩有码中文字幕| 免费观看精品视频网站| 精品电影一区二区在线| 精品久久久久久久毛片微露脸| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 亚洲无线观看免费| 人人妻,人人澡人人爽秒播| 久久香蕉国产精品| 少妇人妻精品综合一区二区 | 日韩国内少妇激情av| 人人妻,人人澡人人爽秒播| 久久精品综合一区二区三区| 久久久国产成人免费| 最近最新中文字幕大全电影3| 波野结衣二区三区在线 | 久久久久性生活片| av专区在线播放| 99久久精品国产亚洲精品| 可以在线观看毛片的网站| 国产午夜福利久久久久久| 18禁国产床啪视频网站| 久久久久九九精品影院| 亚洲在线观看片| 精品国产亚洲在线| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 午夜免费男女啪啪视频观看 | 老司机午夜十八禁免费视频| 欧美日韩国产亚洲二区| www日本黄色视频网| 精品免费久久久久久久清纯| 欧美国产日韩亚洲一区| 国产一级毛片七仙女欲春2| 99热只有精品国产| 女警被强在线播放| 午夜视频国产福利| 欧美日本亚洲视频在线播放| 久久6这里有精品| 中出人妻视频一区二区| 国产精品98久久久久久宅男小说| 免费电影在线观看免费观看| 可以在线观看的亚洲视频| 成人一区二区视频在线观看| 日本五十路高清| 90打野战视频偷拍视频| 免费在线观看亚洲国产| 婷婷亚洲欧美| 久久久久久久久久黄片| 午夜精品在线福利| 国产av麻豆久久久久久久| 国产高潮美女av| АⅤ资源中文在线天堂| bbb黄色大片| 真人做人爱边吃奶动态| 精华霜和精华液先用哪个| 成人午夜高清在线视频| 国产在视频线在精品| 小说图片视频综合网站| 变态另类丝袜制服| 国产69精品久久久久777片| 成人永久免费在线观看视频| 天天一区二区日本电影三级| 久久婷婷人人爽人人干人人爱| 国产精品 国内视频| 国产久久久一区二区三区| 国产探花极品一区二区| 欧美日韩国产亚洲二区| 麻豆国产av国片精品| 99久久无色码亚洲精品果冻| 九九热线精品视视频播放| 夜夜爽天天搞| 欧美xxxx黑人xx丫x性爽| 毛片女人毛片| 熟女人妻精品中文字幕| 成年人黄色毛片网站| 麻豆国产av国片精品| 99久久精品国产亚洲精品| svipshipincom国产片| 国产成人av激情在线播放| 欧美日韩国产亚洲二区| 国产精品亚洲一级av第二区| 日本免费一区二区三区高清不卡| 欧美bdsm另类| 丰满乱子伦码专区| 中国美女看黄片| 欧美一级毛片孕妇| 内射极品少妇av片p| 亚洲专区中文字幕在线| 一进一出好大好爽视频| 欧美黑人欧美精品刺激| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 婷婷六月久久综合丁香| 精品乱码久久久久久99久播| 制服丝袜大香蕉在线| 女人十人毛片免费观看3o分钟| 婷婷精品国产亚洲av在线| 五月伊人婷婷丁香| 免费av毛片视频| 欧美性感艳星| 午夜精品久久久久久毛片777| 亚洲不卡免费看| 成人永久免费在线观看视频| 亚洲 国产 在线| 欧美+亚洲+日韩+国产| 欧美日本亚洲视频在线播放| 可以在线观看毛片的网站| 日本一本二区三区精品| 无限看片的www在线观看| 国产高潮美女av| 美女黄网站色视频| 日本a在线网址| e午夜精品久久久久久久| 国产 一区 欧美 日韩| 不卡一级毛片| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| 国产极品精品免费视频能看的| 久久香蕉国产精品| 欧美国产日韩亚洲一区| 国产精品影院久久| 最后的刺客免费高清国语| 免费人成在线观看视频色| 婷婷精品国产亚洲av在线| 日日干狠狠操夜夜爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本熟妇午夜| 十八禁网站免费在线| 午夜激情欧美在线| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 亚洲精品乱码久久久v下载方式 | 欧美3d第一页| 国产精品久久久久久久电影 | 少妇高潮的动态图| 哪里可以看免费的av片| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 黄色视频,在线免费观看| 蜜桃久久精品国产亚洲av| 精品国产亚洲在线| av视频在线观看入口| 国产 一区 欧美 日韩| 好男人电影高清在线观看| 18禁美女被吸乳视频| av黄色大香蕉| 亚洲国产精品999在线| 欧美xxxx黑人xx丫x性爽| 国产欧美日韩精品一区二区| 男人舔女人下体高潮全视频| 午夜激情福利司机影院| h日本视频在线播放| 亚洲内射少妇av| 亚洲精品456在线播放app | 欧美激情久久久久久爽电影| 国产一区二区亚洲精品在线观看| 亚洲男人的天堂狠狠| 不卡一级毛片| 日韩av在线大香蕉| 午夜视频国产福利| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 我的老师免费观看完整版| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 久久久久久九九精品二区国产| 国产一区二区三区在线臀色熟女| 在线播放国产精品三级| 久久国产乱子伦精品免费另类| 在线观看午夜福利视频| 搞女人的毛片| 中文字幕人妻丝袜一区二区| x7x7x7水蜜桃| 亚洲电影在线观看av| 欧美性猛交╳xxx乱大交人| 日本三级黄在线观看| 久久久久久九九精品二区国产| 亚洲七黄色美女视频| 久久久精品大字幕| 日本 av在线| 一级作爱视频免费观看| 久久精品国产亚洲av涩爱 | 久久香蕉精品热| av中文乱码字幕在线| 亚洲熟妇熟女久久| 国产真实乱freesex| 久久精品国产综合久久久| 国产精品国产高清国产av| 亚洲18禁久久av| 夜夜夜夜夜久久久久| 国产老妇女一区| 中文字幕熟女人妻在线| 美女 人体艺术 gogo| 老司机在亚洲福利影院| 久久精品国产清高在天天线| 露出奶头的视频| 两人在一起打扑克的视频| 国产成人欧美在线观看| 午夜免费男女啪啪视频观看 | 亚洲,欧美精品.| 一二三四社区在线视频社区8| 日韩有码中文字幕| 免费人成在线观看视频色| 午夜两性在线视频| 日韩 欧美 亚洲 中文字幕| 欧美性猛交╳xxx乱大交人| 久久久久性生活片| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 18禁美女被吸乳视频| 亚洲人成网站在线播| 手机成人av网站| 亚洲激情在线av| 老司机午夜十八禁免费视频| 亚洲电影在线观看av| 精品免费久久久久久久清纯| 亚洲成av人片在线播放无| 宅男免费午夜| 精品久久久久久成人av| 精品久久久久久,| 欧美色视频一区免费| 婷婷亚洲欧美| 香蕉丝袜av| 精品一区二区三区av网在线观看| 久久99热这里只有精品18| 国内精品久久久久久久电影| 久久久久久九九精品二区国产| 精品久久久久久久末码| 成人国产综合亚洲| 午夜亚洲福利在线播放| 成人亚洲精品av一区二区| 亚洲中文字幕日韩| 老汉色av国产亚洲站长工具| 国内精品久久久久久久电影| 99热精品在线国产| 中出人妻视频一区二区| 久久精品国产综合久久久| 亚洲第一电影网av| 噜噜噜噜噜久久久久久91| 欧美中文日本在线观看视频| 免费无遮挡裸体视频| 日韩 欧美 亚洲 中文字幕| 免费搜索国产男女视频| 日韩欧美在线二视频| 在线a可以看的网站| 一二三四社区在线视频社区8| 国产亚洲欧美98| 成年版毛片免费区| 美女黄网站色视频| 搡老妇女老女人老熟妇| 校园春色视频在线观看| 动漫黄色视频在线观看| 99在线人妻在线中文字幕| 少妇人妻精品综合一区二区 | 18禁在线播放成人免费| 色综合欧美亚洲国产小说| 国产视频一区二区在线看| 男女之事视频高清在线观看| 国产亚洲精品综合一区在线观看| 欧美一区二区国产精品久久精品| 日日夜夜操网爽| 一进一出好大好爽视频| 色尼玛亚洲综合影院| 亚洲精品日韩av片在线观看 | 精品国产超薄肉色丝袜足j| 亚洲av一区综合| 亚洲欧美精品综合久久99| 动漫黄色视频在线观看| 午夜福利在线观看吧| 内射极品少妇av片p| 麻豆国产97在线/欧美| 一级黄片播放器| 激情在线观看视频在线高清| 亚洲性夜色夜夜综合| 欧美丝袜亚洲另类 | 久久久久久久亚洲中文字幕 | 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av| 国产精品爽爽va在线观看网站| 精品一区二区三区av网在线观看| 免费电影在线观看免费观看| 久久久久精品国产欧美久久久| 美女免费视频网站| 亚洲第一欧美日韩一区二区三区| 人人妻人人澡欧美一区二区| 午夜激情福利司机影院| 97碰自拍视频| 精品无人区乱码1区二区| 日韩欧美 国产精品| 亚洲精华国产精华精| 伊人久久精品亚洲午夜| 国产黄色小视频在线观看| 国产精品野战在线观看| 国产私拍福利视频在线观看| 精品久久久久久久久久久久久| eeuss影院久久| 夜夜夜夜夜久久久久| 亚洲电影在线观看av| 国产黄色小视频在线观看| 99热这里只有是精品50| 天堂网av新在线| 欧美色视频一区免费| 国产精品98久久久久久宅男小说| АⅤ资源中文在线天堂| 狂野欧美白嫩少妇大欣赏| 亚洲av五月六月丁香网| 最新中文字幕久久久久| 亚洲av免费高清在线观看| 欧美+日韩+精品| 变态另类丝袜制服| 两人在一起打扑克的视频| 小说图片视频综合网站| 免费看美女性在线毛片视频| 窝窝影院91人妻| 丰满的人妻完整版| 夜夜夜夜夜久久久久| 黄色日韩在线| 亚洲av电影在线进入| 日韩欧美国产在线观看| 伊人久久大香线蕉亚洲五| 亚洲七黄色美女视频| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩无卡精品| 欧美日韩一级在线毛片| 老熟妇乱子伦视频在线观看| 午夜福利免费观看在线| 蜜桃久久精品国产亚洲av| 国产毛片a区久久久久| 国产三级黄色录像| 高清在线国产一区| 在线观看av片永久免费下载| 欧美一区二区国产精品久久精品| 老司机午夜十八禁免费视频| 特级一级黄色大片| 中文字幕人妻丝袜一区二区| 亚洲精品456在线播放app | 午夜福利视频1000在线观看| 日本熟妇午夜| 午夜福利免费观看在线| 制服人妻中文乱码| 日韩欧美精品v在线| 好男人在线观看高清免费视频| 成年版毛片免费区| 成年女人永久免费观看视频| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 91久久精品国产一区二区成人 | 国产极品精品免费视频能看的| 噜噜噜噜噜久久久久久91| 99久久成人亚洲精品观看| 久久精品亚洲精品国产色婷小说| 国产精品一区二区免费欧美| 999久久久精品免费观看国产| 最近视频中文字幕2019在线8| 特大巨黑吊av在线直播| 黄色日韩在线| 亚洲精品美女久久久久99蜜臀| 黄色丝袜av网址大全| 午夜a级毛片| 精品一区二区三区视频在线 | 精品免费久久久久久久清纯| 色综合婷婷激情| 在线免费观看的www视频| 国产久久久一区二区三区| 丁香欧美五月| 18禁黄网站禁片午夜丰满| 午夜福利在线观看免费完整高清在 | 99久久精品一区二区三区| 五月伊人婷婷丁香| 在线观看日韩欧美| 老司机在亚洲福利影院| 精品久久久久久久毛片微露脸| 18禁美女被吸乳视频| 久9热在线精品视频| 午夜日韩欧美国产| 一区福利在线观看| 无人区码免费观看不卡| 麻豆一二三区av精品| 日韩精品中文字幕看吧| 九九热线精品视视频播放| 亚洲人成网站高清观看| 搡老岳熟女国产| 99精品久久久久人妻精品| 亚洲av二区三区四区| 听说在线观看完整版免费高清| 久久伊人香网站| 欧美av亚洲av综合av国产av| www.www免费av| 天天一区二区日本电影三级| 欧美一区二区国产精品久久精品| 99热这里只有精品一区| 久久久久久人人人人人| 国产精品三级大全| 亚洲av免费在线观看| 黄色片一级片一级黄色片| 国产伦精品一区二区三区四那| 1000部很黄的大片| 最近最新免费中文字幕在线| 精品久久久久久久人妻蜜臀av| 国产成年人精品一区二区| 亚洲不卡免费看| 国产三级中文精品| 给我免费播放毛片高清在线观看| 网址你懂的国产日韩在线| 久久久国产精品麻豆| 久久国产精品人妻蜜桃| 国产三级在线视频| 日韩亚洲欧美综合| 久久国产精品影院| 午夜亚洲福利在线播放| a级一级毛片免费在线观看| 亚洲激情在线av| 亚洲av电影不卡..在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲国产色片| 国产极品精品免费视频能看的| 母亲3免费完整高清在线观看| 婷婷亚洲欧美| 国产亚洲精品久久久com| 嫩草影视91久久| 中文在线观看免费www的网站| 亚洲欧美日韩东京热| 国产毛片a区久久久久| 亚洲精品国产精品久久久不卡| 国产在视频线在精品| 久久精品国产自在天天线| 非洲黑人性xxxx精品又粗又长| 99视频精品全部免费 在线| 亚洲国产精品成人综合色| 免费av观看视频| 成人特级av手机在线观看| 99热这里只有是精品50| 美女免费视频网站| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 一级黄色大片毛片| 亚洲精品国产精品久久久不卡| 中文字幕久久专区| 免费高清视频大片| 少妇的丰满在线观看| 人人妻,人人澡人人爽秒播| 国产麻豆成人av免费视频| 欧美一区二区精品小视频在线| 男插女下体视频免费在线播放| 国产黄片美女视频| 亚洲精品国产精品久久久不卡| 啦啦啦韩国在线观看视频| 长腿黑丝高跟| 美女被艹到高潮喷水动态| 琪琪午夜伦伦电影理论片6080| 夜夜躁狠狠躁天天躁| 亚洲av美国av| 99热这里只有精品一区| 伊人久久大香线蕉亚洲五| 国产成人aa在线观看| 在线视频色国产色| 一本综合久久免费| 久久久久国内视频| 伊人久久精品亚洲午夜| 欧美日韩乱码在线| 成人精品一区二区免费| 成年女人永久免费观看视频| 看免费av毛片| 美女黄网站色视频| 在线播放国产精品三级| 熟妇人妻久久中文字幕3abv| 欧美zozozo另类| 久久香蕉精品热| 1024手机看黄色片| 午夜福利18| 丰满人妻熟妇乱又伦精品不卡| 少妇的逼好多水| 国产伦在线观看视频一区| 中国美女看黄片| 91字幕亚洲| 色av中文字幕| 国产av不卡久久| 国产免费男女视频| 国产精品 国内视频| 日韩欧美 国产精品| 亚洲精品456在线播放app | 国产极品精品免费视频能看的| 18禁黄网站禁片午夜丰满| 美女免费视频网站| 亚洲在线自拍视频| 免费大片18禁| 熟女电影av网| 99热6这里只有精品| 色噜噜av男人的天堂激情| 国内精品久久久久精免费| 免费观看人在逋| 女同久久另类99精品国产91| 18禁裸乳无遮挡免费网站照片| 亚洲中文字幕日韩| 琪琪午夜伦伦电影理论片6080| 国产精华一区二区三区| 国产精品国产高清国产av| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 两个人视频免费观看高清| 亚洲国产日韩欧美精品在线观看 | 两个人看的免费小视频| 深夜精品福利| 淫妇啪啪啪对白视频| 美女黄网站色视频| 久久久精品大字幕| 久久久久久国产a免费观看| 成人精品一区二区免费| 亚洲av日韩精品久久久久久密| 国产精品电影一区二区三区| 床上黄色一级片| 看片在线看免费视频| 欧美一区二区国产精品久久精品| 中文字幕av在线有码专区| 丰满乱子伦码专区| 中文资源天堂在线| а√天堂www在线а√下载| 色老头精品视频在线观看| 久久精品91无色码中文字幕| 伊人久久大香线蕉亚洲五| 欧美国产日韩亚洲一区| 天天添夜夜摸| 亚洲电影在线观看av| 性欧美人与动物交配| 欧美黑人欧美精品刺激| 天天一区二区日本电影三级| 午夜精品在线福利| 午夜免费成人在线视频| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| www.www免费av| 又爽又黄无遮挡网站| 2021天堂中文幕一二区在线观| 午夜激情福利司机影院| 亚洲av成人不卡在线观看播放网| 熟女人妻精品中文字幕| 亚洲专区中文字幕在线| 51国产日韩欧美| 欧美又色又爽又黄视频| e午夜精品久久久久久久| 丁香欧美五月| 成人18禁在线播放| 亚洲熟妇中文字幕五十中出| 天堂av国产一区二区熟女人妻| 天天添夜夜摸| 成人高潮视频无遮挡免费网站| 色精品久久人妻99蜜桃| 亚洲国产高清在线一区二区三| 欧美中文日本在线观看视频| 亚洲熟妇熟女久久| 国产伦一二天堂av在线观看| 欧美zozozo另类| 午夜两性在线视频| 久9热在线精品视频| 国产高潮美女av| 女警被强在线播放| 波多野结衣高清无吗| www国产在线视频色| 搡老岳熟女国产| 丰满人妻熟妇乱又伦精品不卡| 99精品欧美一区二区三区四区| 午夜视频国产福利| 免费观看人在逋| 黄色日韩在线| 欧美色视频一区免费| 亚洲av成人不卡在线观看播放网| 免费看日本二区| 日本 av在线| 麻豆一二三区av精品| 久久久久国内视频| 亚洲精品日韩av片在线观看 | 国产精品一区二区三区四区免费观看 | 丰满的人妻完整版| 亚洲天堂国产精品一区在线| 久久精品影院6| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 免费av不卡在线播放| 久久性视频一级片| 亚洲精品在线观看二区| 亚洲av不卡在线观看| 88av欧美| 国产亚洲精品av在线| 天天躁日日操中文字幕| 国产一级毛片七仙女欲春2| 美女高潮的动态| 他把我摸到了高潮在线观看| 免费在线观看日本一区| 性色av乱码一区二区三区2| 国产精品99久久99久久久不卡| 香蕉丝袜av| 国产美女午夜福利| 啦啦啦免费观看视频1| 欧美大码av| 噜噜噜噜噜久久久久久91| 中文资源天堂在线|