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

    異齡復(fù)層混交對(duì)馬尾松人工林土壤團(tuán)聚體碳組分和轉(zhuǎn)化的影響

    2022-05-13 04:22:57肖納莫雪青譚許脈蘇小艷顏金柳高冠女張文黃雪蔓尤業(yè)明
    廣西植物 2022年4期
    關(guān)鍵詞:土壤酶活性

    肖納 莫雪青 譚許脈 蘇小艷 顏金柳 高冠女 張文 黃雪蔓 尤業(yè)明

    摘要:? 復(fù)層混交種植的近自然化經(jīng)營(yíng)是我國(guó)亞熱帶地區(qū)最常用的營(yíng)林模式之一,但目前其對(duì)土壤有機(jī)碳(SOC)的積累及其穩(wěn)定性的影響仍不明確。為探討該營(yíng)林模式對(duì)土壤碳組分及其轉(zhuǎn)化的影響,該文選取南亞熱帶復(fù)層混交經(jīng)營(yíng)的馬尾松近自然林(CP)作為研究對(duì)象,并以相鄰未改造的馬尾松純林(PP)作為對(duì)照(CK),重點(diǎn)探究異齡復(fù)層混交后原土和不同粒徑土壤團(tuán)聚體的穩(wěn)定性、有機(jī)碳組分、微生物生物量和酶活性等的變化規(guī)律及其之間的關(guān)系。結(jié)果表明:(1)復(fù)層混交后表征土壤團(tuán)聚體穩(wěn)定性的平均重量直徑(MWD)顯著提高;與PP相比,CP原土和不同粒徑土壤團(tuán)聚體的SOC、高度不穩(wěn)定有機(jī)碳(LOC)、中度不穩(wěn)定有機(jī)碳(IOC)、惰性有機(jī)碳(ROC)、土壤碳氮比(C/N)、微生物生物量碳和氮(MBC和MBN)、微生物呼吸(MR)、β-葡萄糖苷酶(BG)、纖維二糖水解酶(CB)、N-乙酰-葡萄糖苷酶(NAG)和酚氧化酶(PO)均顯著提高,而其MBC/MBN和微生物代謝熵(qCO2)卻顯著降低。(2)幾乎所有原土和土壤團(tuán)聚體的SOC、總活性有機(jī)碳(LOC+IOC)、ROC和ROC數(shù)(RCI)均與MWD呈顯著或極顯著正相關(guān)。(3)SOC及其組分的變化主要來(lái)自于土壤理化性質(zhì)、微生物性質(zhì)和植物性質(zhì)三者之間的交互作用,而復(fù)層混交后導(dǎo)致的凋落物碳氮比(C/NLF)和土壤全氮(TN)的顯著變化是驅(qū)動(dòng)SOC及其組分發(fā)生改變的最關(guān)鍵因子。研究認(rèn)為,復(fù)層混交是一種有效促進(jìn)馬尾松人工林SOC及其穩(wěn)定性組分形成的近自然經(jīng)營(yíng)模式,該研究結(jié)果為提高馬尾松人工林的土壤質(zhì)量和穩(wěn)定性提供了科學(xué)依據(jù)。

    關(guān)鍵詞: 馬尾松人工林, 土壤團(tuán)聚體, 有機(jī)碳組分, 土壤酶活性, 近自然化經(jīng)營(yíng)

    中圖分類號(hào):? Q948文獻(xiàn)標(biāo)識(shí)碼:? A文章編號(hào):? 1000-3142(2022)04-0595-13

    Effects of multi-layer and mixed-age forest management?of Pinus massoniana plantations on carbon components?and transformation of soil aggregates

    XIAO Na MO Xueqing TAN Xumai SU Xiaoyan YAN Jinliu GAO Guannü ZHANG WenHUANG Xueman YOU Yeming

    ( 1. Guangxi Key Laboratory of Forest Ecology and Conservation, Forestry College of? Guangxi University, Nanning 530004, China;

    2. Guangxi Youyiguang Forest Ecosystem National Research Station, Pingxiang 532600, Guangxi, China )

    Abstract:? Close-to-nature forest management of multi-layer mixing is one of the most common forest conversion patterns in subtropical regions of China, but its effects on soil organic carbon (SOC) accumulation and stability are still unclear. In this study, the close-to-nature plantations of masson pine(Pinus massoniana) (CP) by multi-layer mixed management in the southern subtropical region were selected as the research object, and the adjacent pure plantations of masson pine (PP) without modification were taken as the control (CK). And then we investigated the stability, organic carbon component, microbial biomass and enzyme activity of soil aggregates with different particle sizes and their relationships. The results were as follows: (1) Mean weight diameter (MWD), which was used to characterize the stability of soil aggregates, increased significantly after near-natural modification. Compared to PP, SOC, labile organic carbon (LOC), intermediate organic carbon (IOC), recalcitrant organic carbon (ROC), soil C/N, microbial biomass carbon and nitrogen (MBC and MBN), microbial respiration (MR), β-glucosidase (BG), cellobiohydrolase (CB), N-acetyl-glucosaminidase (NAG) and phenol oxidase (PO) were significantly increased, but the MBC/MBN and microbial metabolic quotient (qCO2) were significantly decreased. (2) The SOC, total organic carbon(LOC+IOC), ROC and ROC index (RCI) of in soils or most soil aggregates were significantly positively correlated with MWD. (3) The change of SOC content and its components mainly resulted from the interaction among soil physical and chemical, microbial and plant properties. Furthermore, the significant changes of litter carbon to nitrogen ratio (C/NLF) and soil total nitrogen (TN) caused by masson pine with multi-layer mixed forest management were the primary drivers of differences in SOC content and its components. Our findings indicate that that the close-to-natural forest management model of multi-layer mixing is an effective management measure to improve SOC and its stability components, and so as to promote masson pine plantations to sequester more and more stabilized SOC in the subtropical region of China.

    Key words: masson pine (Pinus massoniana) plantations, soil aggregates, organic carbon components, soil enzyme activities, close-to-nature forest management

    增加森林碳匯,是實(shí)現(xiàn)“碳達(dá)峰、碳中和”的一個(gè)重要途徑,而森林土壤碳庫(kù)由于具有巨大的碳固持潛力,其在減緩大氣CO2濃度快速升高和調(diào)控全球生態(tài)系統(tǒng)的碳平衡方面發(fā)揮著極其關(guān)鍵的作用,因此成為目前世界各國(guó)應(yīng)對(duì)全球氣候變化重點(diǎn)探索的領(lǐng)域(Bossio et al.,2020)。我國(guó)人工林面積位居世界第一,約7 954.28萬(wàn)hm2,具有巨大的碳匯功能,但其土壤碳扣押潛力受樹(shù)種、種植模式和管理措施等方面的影響。因此,如何才能增加人工林土壤碳的儲(chǔ)量和穩(wěn)定性,已成為當(dāng)前國(guó)內(nèi)外許多科學(xué)家關(guān)注的焦點(diǎn)和亟需解決的問(wèn)題。

    “近自然森林經(jīng)營(yíng)”是一種最早源于德國(guó)的兼容林業(yè)生產(chǎn)和生態(tài)保護(hù)的森林經(jīng)營(yíng)模式 ,長(zhǎng)期實(shí)踐證明了該模式能在保持森林自然環(huán)境健康穩(wěn)定發(fā)展的前提下達(dá)到不同的收獲目標(biāo),從而使得森林綜合效益達(dá)到最大化(Moradi et al.,2012;牟耀杰等,2019)。針葉純林間伐后套種多種闊葉樹(shù)是一種快速形成針闊混交、異齡、復(fù)層結(jié)構(gòu)人工林群落的重要方法,是近自然改造經(jīng)營(yíng)的一項(xiàng)主要技術(shù)措施(陸元昌和甘敬,2002)。針葉純林經(jīng)復(fù)層混交的近自然化經(jīng)營(yíng)后,其樹(shù)種組成和林分結(jié)構(gòu)發(fā)生了明顯改變,影響了人工林群落的生產(chǎn)力、凋落物的數(shù)量和質(zhì)量、根系及其分泌物以及物種多樣性。前人的研究表明,經(jīng)近自然化改造的針葉純林可達(dá)到優(yōu)化林分結(jié)構(gòu),增加物種多樣性,改變土壤理化性質(zhì),提升人工林生產(chǎn)力等目的(孫冬婧等,2015;明安剛等,2017;Ming et al.,2018)。因此,人工針葉純林經(jīng)過(guò)近自然化改造形成針/闊異齡混交林后,通過(guò)改變林分微環(huán)境、凋落物的數(shù)量和質(zhì)量、根系及其分泌物等,勢(shì)必會(huì)對(duì)土壤的物理結(jié)構(gòu)(尤其是團(tuán)聚體結(jié)構(gòu))、土壤微生物群落組成和功能以及其他土壤理化性質(zhì)產(chǎn)生更為復(fù)雜的影響,最終影響到土壤碳的積累和穩(wěn)定性。

    馬尾松(Pinus massoniana)作為重要的用材林樹(shù)種及荒山造林的先鋒樹(shù)種之一,因其具有耐干旱瘠薄、速生豐產(chǎn)、適應(yīng)性強(qiáng)等特性而在我國(guó)亞熱帶地區(qū)被廣泛種植,并給人們帶來(lái)了巨大的經(jīng)濟(jì)價(jià)值與社會(huì)效益(Fu et al.,2011;楊予靜等,2018)。但是,早期馬尾松的營(yíng)林措施主要關(guān)注地上木材產(chǎn)量,而忽視了地下土壤質(zhì)量的提升,造成馬尾松人工林經(jīng)營(yíng)結(jié)構(gòu)單一,從而導(dǎo)致其生產(chǎn)力低下、土壤的團(tuán)粒結(jié)構(gòu)失去穩(wěn)定性、養(yǎng)分流失嚴(yán)重和固碳潛力受限等(You et al.,2018;梁艷等,2021)。利用引入闊葉樹(shù)種,如:紅錐(Castanopsis hystrix)、香梓楠(Michelia hedyosperma)和格木(Erythrophleum fordii)等對(duì)馬尾松進(jìn)行混交改造,建立系統(tǒng)生產(chǎn)力更高和生態(tài)系統(tǒng)服務(wù)功能更強(qiáng)的針/闊異齡混交林替代原有大面積分布的低質(zhì)低效馬尾松純林,這一做法已成為我國(guó)亞熱帶地區(qū)最常用的營(yíng)林改造模式之一,也是目前林學(xué)、生態(tài)學(xué)研究的熱點(diǎn)問(wèn)題之一(劉世榮等,2018)。然而,馬尾松人工純林通過(guò)引入闊葉樹(shù)種異齡復(fù)層混交后,其凋落物數(shù)量和質(zhì)量、根系生物量和分泌物等均發(fā)生了不同程度的改變,這些改變將對(duì)土壤團(tuán)聚體的碳組分和碳轉(zhuǎn)化速率產(chǎn)生哪些影響,是否有利于提高土壤碳的固持和穩(wěn)定性,與土壤團(tuán)聚體穩(wěn)定性存在何種關(guān)系等。對(duì)這些問(wèn)題仍缺乏深入的探究和理解。因此,本研究選擇南亞熱帶的馬尾松純林間伐后套種闊葉樹(shù)種形成的異齡復(fù)層近自然林為試驗(yàn)對(duì)象,并以相鄰的馬尾松純林(未改造)作為對(duì)照(CK),重點(diǎn)探究異齡復(fù)層混交后土壤團(tuán)聚體穩(wěn)定性的變化特征;原土和不同徑級(jí)土壤團(tuán)聚體碳組分和土壤酶活性的有關(guān)變化;土壤碳組分與土壤酶活性、環(huán)境因子和土壤團(tuán)聚體穩(wěn)定性之間的關(guān)系,以期為提高低效馬尾松人工林的土壤碳固持及穩(wěn)定性等方面提供有效的科學(xué)依據(jù)。

    1材料與方法

    1.1 研究區(qū)和試驗(yàn)地概況

    本試驗(yàn)區(qū)域位于廣西友誼關(guān)森林生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站內(nèi)的中國(guó)林科院熱帶林業(yè)實(shí)驗(yàn)中心的人工林試驗(yàn)示范基地(106°51′—106°53′ E,22°02′—22°04′ N)。該區(qū)域地處亞熱帶,具有亞熱帶季風(fēng)氣候,夏季高溫多雨,冬季溫暖少雨,土壤類型以花崗巖風(fēng)化后形成的紅壤為主,年均溫21 ℃,年均降雨量1 400 mm。其典型的人工針葉林以馬尾松和杉木(Cunninghamia lanceolata)為主,人工闊葉林主要有紅錐、格木、大葉櫟(Quercus griffithii)、日本厚樸(Magnolia obotava)和鐵力木(Mesua ferrea)等(莫雪青等,2022)。

    本研究選擇的馬尾松近自然林(close-to-nature plantations of masson pine, CP)的改造過(guò)程如下:選擇立地條件、經(jīng)營(yíng)歷史和生長(zhǎng)狀況基本一致的馬尾松純林(1993種植,初植密度為 2 500株·hm-2),分別在2000年、2004年和2007年對(duì)其進(jìn)行間伐,最終保留的密度為 450株·hm-2。于2008年春季,對(duì)間伐后的馬尾松林均勻補(bǔ)植紅錐和香梓楠1年生幼苗(實(shí)生苗),補(bǔ)植的紅錐和香梓楠密度均為 375株·hm-2,最終經(jīng)過(guò)異齡復(fù)層混交后的近自然林的密度為 1 200株·hm-2(馬尾松∶紅錐∶香梓楠=6∶5∶5),同時(shí)選擇鄰近未進(jìn)行改造的馬尾松純林(pure plantations of masson pine, PP)(林分密度與近自然林總密度一致,1 200株·hm-2)作為對(duì)照(CK)。

    根據(jù)隨機(jī)區(qū)組設(shè)計(jì)原則,于2020年7月在試驗(yàn)區(qū)隨機(jī)選擇了4個(gè)試驗(yàn)區(qū)組,其中每個(gè)區(qū)組均包含CP和PP 2種林分類型,區(qū)組內(nèi)的CP和PP各設(shè)置一個(gè)20 m × 20 m的樣方,整個(gè)實(shí)驗(yàn)共設(shè)置了8個(gè)400 m2的樣方,并采用傳統(tǒng)的樣地調(diào)查法對(duì)每個(gè)實(shí)驗(yàn)樣方進(jìn)行調(diào)查。其中,在每個(gè)樣方中隨機(jī)布設(shè)6個(gè)的尼龍網(wǎng)收集框(1 m × 1 m)用于監(jiān)測(cè)其年凋落物的量(litter fallmass, LF,kg·hm-2·yr-1),每個(gè)樣方的細(xì)根生物量(fine root mass, FR)采用土鉆法獲得,所有的植物樣品烘干(65 ℃)至恒重,經(jīng)粉碎后測(cè)定其碳和氮的含量(陳永康等,2021)。樣地的基本信息詳見(jiàn)表1。

    1.2 土壤樣品的采集和處理

    于2020年8月分別對(duì)每個(gè)實(shí)驗(yàn)樣方進(jìn)行土壤樣品采集??紤]混交后的馬尾松人工林樹(shù)木具有隨機(jī)分布的特點(diǎn),本研究采用系統(tǒng)樣點(diǎn)布設(shè)法來(lái)確定土壤采樣點(diǎn),即把每個(gè)樣方(20 m × 20 m)分成16個(gè)5 m × 5 m大小相等的網(wǎng)格,網(wǎng)格線的交點(diǎn)即為采樣點(diǎn)(9個(gè))。在采用點(diǎn)去除表面凋落物和雜質(zhì)后在進(jìn)行原狀土(0~10 cm土層)的采集,其規(guī)格(長(zhǎng)×寬×高)為15 cm × 10 cm × 10 cm。為防止原狀土結(jié)構(gòu)遭到破壞,特將采集好的土壤塊放入硬質(zhì)PVC塑料盒后置于保溫箱中低溫保存。把所有土壤帶回實(shí)驗(yàn)室后將同一樣方的9個(gè)采樣點(diǎn)土塊混合均勻, 剔除雜質(zhì)后將其分成兩份。一份土壤過(guò)2 mm土篩作為原土,另一過(guò)2 mm和0.25 mm的套篩,得到粒級(jí)為>2 mm、0.25~2 mm和<0.25 mm的3種不同粒徑土壤團(tuán)聚體。將處理后的土壤樣品保存在-20 ℃的冰箱中,以用于測(cè)定酶活性、有機(jī)碳組分、微生物呼吸和其他土壤理化性質(zhì)等指標(biāo)。

    1.3 土壤樣品的理化性質(zhì)分析

    土壤容重(bulk density, BD)指標(biāo)采用環(huán)刀法測(cè)定;土壤含水量(soil water content, SWC,%)指標(biāo)是將新鮮土壤于105 ℃烘干至恒重后測(cè)定;土壤pH值是利用pH計(jì)測(cè)定(土∶水=1∶2.5,w∶v);土壤總氮(total nitrogen, TN)依照凱氏定氮法進(jìn)行測(cè)定;土壤銨態(tài)氮(ammonia nitrogen, NH4+-N)和硝態(tài)氮(nitrate nitrogen, NO3--N)含量在KCl(2 mol·L-1)溶液浸提后,將浸提液置于連續(xù)流動(dòng)分析儀(SEAL Auto Analyzer 3)上分析測(cè)定(鮑士旦,2000)。

    土壤碳化學(xué)組分采用酸解法進(jìn)行測(cè)定(Guo et al.,2021):稱取風(fēng)干的土壤樣品(0.50 g)置于100 mL離心管中,①加入20 mL 1.25(2.5N)H2SO4溶液,在105 ℃下水解0.5 h,離心(5 000 r·min-1,10 min)取其上清液,稀釋后測(cè)定的上清液碳含量即為高度不穩(wěn)定有機(jī)碳(labile organic carbon, LOC);②加入2 mL 6.5(13N)H2SO4溶液,振蕩(25 ℃搖床)充分水解(12 h)后加入24 mL超純水,使其溶液的濃度降至0.5(1N),在105 ℃下水解3 h后離心(5 000 r·min-1,10 min)取上清液,稀釋后測(cè)定的上清液碳含量即為中度不穩(wěn)定有機(jī)碳(intermediate organic carbon, IOC),高度不穩(wěn)定的碳與中度不穩(wěn)定的碳之和統(tǒng)稱總活性碳(總活性碳=LOC+ IOC);③剩余土壤殘?jiān)娓桑?0 ℃)至恒重后稱重,測(cè)定其惰性有機(jī)碳含量(recalcitrant organic carbon, ROC)。

    1.4 土壤酶活性、微生物生物量和呼吸的測(cè)定

    采用96孔微孔板法測(cè)定6種土壤酶活性,包括4種水解酶,分別為參與土壤C循環(huán)的β-葡萄糖苷酶(β-glucosidase, BG)和纖維二糖水解酶(cellobiohydrolase, CB)、參與N循環(huán)的N-乙酰-葡萄糖苷酶(N-acetyl-glucosaminidase, NAG)和亮氨酸氨基肽酶(leucine-amino-peptidase, LAP);2種氧化酶,分別為酚氧化酶(phenol oxidase, PO)和過(guò)氧化物酶(peroxidase, PER)。 其中,水解酶活性的測(cè)定是基于甲基傘形酮(MUB)為底物的熒光法,置于酶標(biāo)儀(Thermo Multiskan Spectrum)在波長(zhǎng)365~450 nm下進(jìn)行測(cè)定,每個(gè)樣品8個(gè)重復(fù)(Steinauer et al.,2015;You et al.,2020)。氧化酶 以l -3,4-二羥基苯丙氨酸(L-DOPA)作為底物的比色法進(jìn)行測(cè)定,其中,L-DOPA加0.3% H2O2用于測(cè)定過(guò)氧化物酶活性,兩種氧化酶的活性均置于酶標(biāo)儀在波長(zhǎng)460 nm下進(jìn)行測(cè)定,每個(gè)樣品8個(gè)重復(fù)(Saiya-Cork et al.,2002)。 水解酶和氧化酶的活性單位均統(tǒng)一為nmol·h-1·g-1soil。

    土壤微生物生物量碳(microbial biomass carbon, MBC)和微生物量氮(microbial biomass nitrogen, MBN)采用氯仿熏蒸浸提法進(jìn)行測(cè)定(Vance et al.,1987),而可溶性有機(jī)碳(dissolved organic carbon, DOC)和可溶性有機(jī)氮(dissolved nitrogen, DON)為未經(jīng)過(guò)氯仿熏蒸處理的新鮮土壤用硫酸鉀(K2SO4,0.5 mol·L-1)浸提出的有機(jī)碳和氮含量(You et al.,2014)。土壤微生物呼吸(microbial respiration, MR)采用室內(nèi)培養(yǎng)法進(jìn)行測(cè)定,其分析過(guò)程主要參照(Jin et al.,2010)描述的測(cè)定方法進(jìn)行。微生物代謝熵(microbial metabolic quotient, qCO2)為MR和MBC的比值。

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

    平均重量直徑(mean weight diameter, MWD)根據(jù)Castro提出的公式計(jì)算獲得(Castro et al.,2002),主要用于土壤團(tuán)聚體穩(wěn)定性的評(píng)價(jià),其計(jì)算公式如下:

    式中:n表示團(tuán)聚體個(gè)數(shù);Xi為第i 粒徑團(tuán)聚體直徑的平均值(mm); Wi為第i 粒徑團(tuán)聚體的重量百分含量(%)。

    利用SPSS 19.0的獨(dú)立樣本t檢驗(yàn)(independent sample t test)對(duì)PP和CP之間土壤理化性質(zhì)、土壤酶活性、有機(jī)碳組分和微生物呼吸等進(jìn)行顯著性檢驗(yàn),顯著性水平設(shè)置為P<0.05;采用Pearson 相關(guān)性分析土壤有機(jī)碳(soil organic carbon, SOC)及其組分與MWD之間的相關(guān)性,顯著性水平設(shè)置為P<0.05。

    基于野外監(jiān)測(cè)和室內(nèi)分析,一共得到22個(gè)生物和非生物因子納入影響SOC及其組分的多元統(tǒng)計(jì)分析,包括土壤理化性質(zhì): BD、MWD、TN、NH4+-N、NO3--N、pH、DOC、DON;微生物性質(zhì): MBC、MBN、MBC/MBN、BG、CB、NAG、LAP、PO、PER、MR、qCO2;植物性質(zhì): FR、LF、C/NLF

    采用“Vegan”數(shù)據(jù)包的“Varpart”程序?qū)τ绊慡OC及其組分的環(huán)境因子進(jìn)行方差分解分析,其分析過(guò)程是在R 3.01上運(yùn)行完成的。不同類型環(huán)境因子的純效應(yīng)及其之間的相互效應(yīng)采用韋恩圖進(jìn)行顯示,圖中的數(shù)字代表各類環(huán)境因子對(duì)SOC及其組分的解釋率,其中解釋率為負(fù)值時(shí),在圖中不作顯示,但其解釋率等同于“0” (Sattler et al.,2010)。最后運(yùn)用冗余分析(RDA)的蒙特卡洛檢驗(yàn)(Monte Carlo,499 iterations)對(duì)基于野外監(jiān)測(cè)和室內(nèi)分析得到的22個(gè)環(huán)境因子進(jìn)行預(yù)選,并把具有顯著性影響(P<0.05)的環(huán)境因子作為最終的環(huán)境變量納入RDA模型分析,并通過(guò)模型的排序確定影響SOC及其組分的最主要環(huán)境因子(P<0.05)及其解釋率。RDA程序是在多元統(tǒng)計(jì)分析軟件Canoco for Windows 4.5上運(yùn)行完成,所有數(shù)據(jù)在分析之前均經(jīng)過(guò)對(duì)數(shù)轉(zhuǎn)換(You et al.,2014)。運(yùn)用軟件Sigma Plot 14.0完成所有的作圖。

    2結(jié)果與分析

    2.1 復(fù)層混交對(duì)馬尾松人工林土壤團(tuán)聚體粒徑分布和理化性質(zhì)的影響

    研究結(jié)果表明,相對(duì)于PP, CP的不同粒徑土壤團(tuán)聚體重量百分比發(fā)生了不同程度的變化,其中,CP的>2 mm團(tuán)聚體重量百分比顯著高于PP(P<0.05),而0.25~2 mm和<0.25 mm團(tuán)聚體在兩種林分中無(wú)顯著差異(P>0.05)。本研究中,CP土壤團(tuán)聚體的MWD顯著高于PP,表明馬尾松人工純林經(jīng)復(fù)層混交后其土壤團(tuán)聚體的穩(wěn)定性顯著提高(表2)。

    相比于PP,CP的SOC、C/N、MBC和MBN在原土和各個(gè)徑級(jí)的土壤團(tuán)聚體中均顯著提高(P<0.05),CP的TN和NH4+-N在原土和大團(tuán)聚體(>2 mm)中顯著提高(P<0.05);而CP的MBC/MBN、DOC和DON在原土和大多數(shù)徑級(jí)土壤團(tuán)聚體中卻顯著降低(P<0.05),NO3--N和pH值在原土和各個(gè)徑級(jí)的土壤團(tuán)聚體中均無(wú)顯著變化(P>0.05)(表3)。

    2.2 復(fù)層混交對(duì)馬尾松人工林原土和土壤團(tuán)聚體碳組分的影響

    CP的LOC、IOC、ROC含量在原土及所有粒徑團(tuán)聚體均與PP存在顯著或極顯著差異。其中, LOC在原土、 >2 mm、 0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體分別增加了79.59%、75.22%、63.22%和53.95%(圖1:A);IOC在原土、>2 mm、0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體分別增加了18.29%、12.91%、16.18%和24.62%(圖1:B);ROC在原土、>2 mm、0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體分別增加了56.33%、72.67%、48.55%和101.67%(圖1:C),但CP的RCI與PP只在原土、>2 mm和<0.25 mm粒徑存在顯著和極顯著差異,并分別增加了10.17%、13.42% 和16.71%,而在0.25~2 mm粒徑土壤團(tuán)聚體無(wú)顯著差異(P>0.05)(圖1:D)。

    MWD是一種表示土壤物理結(jié)構(gòu)穩(wěn)定性的指標(biāo)。如圖2所示,SOC在原土、0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體中與MWD呈顯著線性正相關(guān)(P<0.05),而在>2 mm粒徑土壤團(tuán)聚體中與MWD呈極顯著線性正相關(guān)(P<0.01)(圖2:A);總活性有機(jī)碳(LOC+IOC)在原土、0.25~2 mm和<0.25 mm與MWD呈顯著或極顯著線性正相關(guān),但在>2 mm粒徑土壤團(tuán)聚體中與MWD的關(guān)系不顯著(P>0.05)(圖2:B);ROC在原土、>2 mm和0.25~2 mm粒徑土壤團(tuán)聚體中與MWD呈顯著線性正相關(guān)(P<0.05),但在<0.25 mm粒徑土壤團(tuán)聚體中與MWD的關(guān)系不顯著(P>0.05)(圖2:C);RCI只在>2 mm粒徑土壤團(tuán)聚體中與MWD呈顯著線性正相關(guān)(P<0.05),而在原土、0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體中均與MWD的關(guān)系不顯著(P>0.05)(圖2:D)。

    2.3 復(fù)層混交對(duì)馬尾松人工林原土和土壤團(tuán)聚體酶活性和微生物呼吸的影響

    數(shù)據(jù)分析結(jié)果表明,CP的BG、CB、NAG、PO在原土及所有粒徑團(tuán)聚體均與PP存在顯著或極顯著差異。其中,BG在原土、>2 mm、0.25~2 mm和<0.25 mm徑級(jí)團(tuán)聚體分別增加了237.34%、264.66%、152.42%和227.29%(圖3:A); CB 在原土、>2 mm、0.25~2 mm和<0.25 mm粒徑的團(tuán)聚體中分別增加了109.67%、124.89%、154.57%和82.88%(圖3:B); NAG在原土、>2 mm、0.25~2 mm及<0.25 mm團(tuán)聚體中分別極顯著增長(zhǎng)了71.63%、74.36%、57.81%和81.26%(圖3:C); PO在原土、>2 mm、0.25~2 mm及<0.25 mm團(tuán)聚體分別增加了40.59%、44.67%、78.12%和34.96%(圖3:E);CP的LAP僅在原土和0.25~2 mm土壤小團(tuán)聚體與PP存在極顯著差異,分別增加了33.47%和56.16%,在>2 mm及<0.25 mm團(tuán)聚體中均無(wú)顯著差異(圖3:D),但PER不管在原土還是團(tuán)聚體中均無(wú)顯著變化(圖3:F)。

    CP的MR在原土及>2 mm、0.25~2 mm和<0.25 mm粒徑的團(tuán)聚體中與PP存在極顯著差異,并分別提升了14.24%、17.90%、17.33%和20.66%(圖4:A),但CP的qCO2不僅在原土、0.25~2 mm和<0.25 mm粒徑土壤團(tuán)聚體中與PP存在顯著或極顯著差異,且分別降低了23.49%、34.51%和19.47%,而在>2 mm粒徑的大團(tuán)聚體中不存在顯著差異(圖4:B)。

    2.4 土壤碳組分的影響因素

    如圖5所示,本研究所調(diào)查的3類環(huán)境因子(土壤理化性質(zhì)、微生物性質(zhì)和植物性質(zhì))的總疊加效應(yīng)解釋了SOC及其組分變異的69.5%。其中,影響SOC及其組分的因素主要來(lái)自于土壤理化性質(zhì)、微生物性質(zhì)和植物性質(zhì)三類因子之間的交互作用,解釋了其變化的55.0%(圖5:g);其次是土壤理化性質(zhì)和微生物性質(zhì)兩類因子之間的交互作用(12.6%)(圖5:d);而微生物性質(zhì)和植物性質(zhì)兩類因子之間的交換作用對(duì)SOC及其組分變化的解釋率較?。?.4%)(圖5:e)。植物性質(zhì)的獨(dú)立效應(yīng)解釋了7.8%的SOC及其組分變化(圖5:c),而單獨(dú)的土壤理化性質(zhì)、微生物性質(zhì)以及土壤理化性質(zhì)和微生物性質(zhì)兩類因子之間的交互作用沒(méi)有對(duì)SOC及其組分變化產(chǎn)生影響(圖5:a, b, f)。

    在所調(diào)查的土壤理化性質(zhì)、微生物性質(zhì)和植物性質(zhì)因子3類環(huán)境因子中,先通過(guò)RDA程序?qū)ζ鋯我画h(huán)境因子進(jìn)行蒙特卡洛檢驗(yàn)預(yù)選后得到17個(gè)環(huán)境因子與SOC及其組分顯著相關(guān)并被選入作為模型分析的環(huán)境變量,再通過(guò)RDA排序確定C/NLF和TN是驅(qū)動(dòng)馬尾松人工林SOC及其組分發(fā)生變化的最主要環(huán)境因子(P<0.05)(圖6)。

    3討論

    對(duì)森林SOC的起源及其組分的探索,將有助于確保森林SOC能夠成為減緩全球大氣CO2快速上升的一個(gè)可靠的碳匯。土壤中的SOC主要來(lái)源于植物的光合作用,地表凋落物、根系及根系分泌物等的合成轉(zhuǎn)化,因此,地上和地下有機(jī)碳輸入的數(shù)量和質(zhì)量將影響SOC的儲(chǔ)量和穩(wěn)定性(Huang et al.,2017)。Wang 等(2018)研究發(fā)現(xiàn)馬尾松與紅錐混交種植可加速頑固針葉凋落物的分解,減少不穩(wěn)定闊葉凋落物的碳排放,從而有效促進(jìn)碳固存;而盧立華等(2019)的研究結(jié)果表明馬尾松與紅錐混交種植可顯著提升生態(tài)系統(tǒng)碳儲(chǔ)量(喬木層、地被層和土壤層的碳儲(chǔ)量之和);此外,還有研究表明,馬尾松經(jīng)復(fù)層混交套種兩種鄉(xiāng)土闊葉樹(shù)種后均能促進(jìn)林分空間結(jié)構(gòu)的優(yōu)化,可明顯提升林分的生物量和植被的固碳能力(劉志龍等,2017;明安剛等,2017)。在本研究中,馬尾松純林經(jīng)間伐后套種鄉(xiāng)土闊葉樹(shù)種的復(fù)層混交后,SOC在原土和各個(gè)粒徑土壤團(tuán)聚體中均顯著提高,這與Schütz等(2011)的研究結(jié)果基本一致,經(jīng)復(fù)層混交后林分的凋落物產(chǎn)量和細(xì)根生物量均顯著提高,為SOC的積累提供了豐富的有機(jī)物輸入源,有效提高SOC含量。根據(jù)SOC的穩(wěn)定性差異及其在土壤中周轉(zhuǎn)速率的快慢,可將其劃分為活性碳組分(LOC和IOC)和惰性碳組分(ROC)兩大類。在森林生態(tài)系統(tǒng)中,活性碳對(duì)維持其土壤肥力和碳平衡具有重要的積極作用,而ROC有利于土壤SOC的積累和長(zhǎng)期固持,在減緩全球氣候變化中扮演重要角色(Jandl et al.,2007;Kooch et al.,2019)。鮑勇等(2018)的研究發(fā)現(xiàn),SOC的組分特別是活性組分主要來(lái)源于新近的葉凋落物、根系分泌物、土壤腐殖質(zhì)和微生物殘?bào)w等;而根據(jù)Fan等(2004)和Simpson等(2007)的研究表明,不同植物凋落物的理化性質(zhì)存在較大差異將影響其分解速率和進(jìn)入土壤的化學(xué)組分,最終影響SOC及其組分。此外,越來(lái)越多的研究表明MBC不僅是形成SOC的最重要前體,還是形成SOC中復(fù)雜穩(wěn)定性組分的重要組成部分,輸入土壤的有機(jī)物經(jīng)微生物轉(zhuǎn)化后將可能形成更加復(fù)雜的聚合物(Simpson et al.,2007;Miltner et al.,2012)。本研究結(jié)果印證了前人的研究觀點(diǎn),復(fù)層混交的近自然營(yíng)林模式不但可改善林分結(jié)構(gòu),提高土壤肥力和生產(chǎn)力,而且可有效提高原土和不同徑級(jí)土壤團(tuán)聚體的MBC、SOC含量及其穩(wěn)定性組分。

    MWD是評(píng)價(jià)團(tuán)聚體穩(wěn)定性高低的重要指標(biāo),而團(tuán)聚體的穩(wěn)定性又直接關(guān)系到土壤碳含量及其組分穩(wěn)定性(Dungait et al.,2012)。有研究表明,團(tuán)聚體儲(chǔ)存著超過(guò)90%的土壤碳,其組成和穩(wěn)定性與土壤碳固持和養(yǎng)分的有效性具有密切關(guān)系(Six et al.,2000;錢婧等,2018)。本研究發(fā)現(xiàn),原土和大部分土壤團(tuán)聚體中的SOC及其組分均與MWD呈顯著或極顯著正相關(guān)關(guān)系,而近自然化改造也顯著促進(jìn)了MWD,這可能與近自然化經(jīng)營(yíng)提高C源輸入在促進(jìn)土壤團(tuán)聚體形成的同時(shí)也增加了土壤SOC及其組分的積累有關(guān)。這與顧偉等(2007)的研究相似,認(rèn)為隨著土壤團(tuán)聚體穩(wěn)定性的提高,可以顯著影響土壤中的不同類型碳組分含量,從而改變?cè)械奶脊坛稚踔撂佳h(huán),并最終有效地改善土壤質(zhì)量。

    土壤胞外酶是森林土壤的重要組成部分,其大部分由土壤微生物代謝產(chǎn)生,只有少數(shù)一部分由土壤動(dòng)物、根系分泌及其殘?bào)w的分解釋放。然而,一些胞外酶由特定的微生物類群分泌用于降解某些特定的化合物,不同的酶活性在參與土壤生物地球化學(xué)循環(huán)過(guò)程具有不同作用(Allison et al.,2010)。因此,林分結(jié)構(gòu)、物種組成、凋落物數(shù)量和質(zhì)量、根系及其分泌物等能夠影響土壤微生物物種多樣性和群落組成的環(huán)境要素均可能改變土壤酶活性,影響土壤碳的循環(huán)過(guò)程(You et al.,2014; Lucas-Borja & Delgado-Baquerizo,2019)。本研究發(fā)現(xiàn),馬尾松純林經(jīng)復(fù)層混交后原土及各粒徑團(tuán)聚體中參與土壤碳和氮循環(huán)有關(guān)的酶活性(BG、 CB、NAG和 PO)均有顯著提高,這可能與改造后近自然林的微生物生物量(采用MBC和MBN表征)顯著提高有關(guān),Brockett等(2012)研究也發(fā)現(xiàn)土壤酶活性的提高與土壤微生物生物量增加密切相關(guān)。但是,PER的活性在未改造的純林和改造后的近自然林之間并未發(fā)生明顯的變化,這與Sinsabaugh 等(2008)的研究結(jié)論相似,發(fā)現(xiàn)PER的活性與土壤SOC無(wú)顯著的相關(guān)性,而是隨著 pH 的增加而呈現(xiàn)增加的趨勢(shì),而本研究的近自然林SOC雖有顯著提高,但其pH值與未改造的馬尾松純林并無(wú)明顯差異。此外,本研究中不同的土壤酶活性對(duì)團(tuán)聚體粒徑大小的響應(yīng)無(wú)特定規(guī)律,這也暗示了土壤團(tuán)聚體粒徑大小對(duì)土壤酶活性的影響過(guò)程和機(jī)制仍較為復(fù)雜,尚需開(kāi)展更加長(zhǎng)期和科學(xué)的研究。

    本研究通過(guò)方差分解分析得到所調(diào)查的3類環(huán)境因子(土壤理化性質(zhì)、微生物性質(zhì)和植物特性)的總疊加效應(yīng)解釋了SOC及其組分變化的69.5%,這表明了所調(diào)查的環(huán)境因子能夠解釋土壤SOC及其組分的大部分變化,但也暗示一些未調(diào)查的因子(如:凋落物的木質(zhì)素含量、土壤動(dòng)物種類和土壤微生物多樣性等)仍可能是驅(qū)動(dòng)土壤SOC及其組分發(fā)生變化的重要因素(Prescott et al.,2010;Zhang et al.,2013)。本研究也發(fā)現(xiàn)了土壤SOC及其組分的大部分變化主要來(lái)自于土壤理化性質(zhì)、微生物性質(zhì)和植物特性三者之間交互作用的解釋,而三者中的單獨(dú)效應(yīng)或兩者之間的交互作用對(duì)土壤SOC及其組分變化的解釋率較低甚至沒(méi)有,這與Schmidt等(2011)的觀點(diǎn)一致,認(rèn)為土壤SOC的固持及其組分的形成主要取決于生物和非生物環(huán)境復(fù)雜的綜合作用。而通過(guò)RDA進(jìn)一步分析的結(jié)果揭示了馬尾松近自然化改造后導(dǎo)致的凋落物碳氮比(C/NLF)和土壤TN的顯著變化是驅(qū)動(dòng)SOC及其組分發(fā)生改變的最主要因子。這與前人的研究結(jié)論一致,高的凋落物質(zhì)量(低的C/N)和土壤氮含量更有利于提高土壤碳的輸入、土壤微生物生物量和土壤腐殖質(zhì)的形成,更有利于土壤碳的積累和穩(wěn)定性組分的形成,促使森林固持更多更穩(wěn)定的SOC(Prescott et al.,2010;Huang et al.,2017;Chen et al.,2018)。

    4結(jié)論

    馬尾松人工純林經(jīng)異齡復(fù)層混交改造后改變了林分結(jié)構(gòu),顯著增加了凋落物產(chǎn)量和細(xì)根生物量,并顯著提高了土壤團(tuán)聚體的穩(wěn)定性、MBC、MBN、MR、參與土壤碳和氮循環(huán)有關(guān)的酶活性(BG、 CB、NAG和 PO)、SOC及其活性和惰性組分,但顯著降低了MBC/MBN和qCO2;而幾乎所有原土和土壤團(tuán)聚體的SOC、總活性有機(jī)碳、ROC和RCI均與表征土壤團(tuán)聚體穩(wěn)定性的MWD呈顯著或極顯著正相關(guān)關(guān)系。本研究揭示了近自然化經(jīng)營(yíng)后導(dǎo)致的林分C/NLF和土壤TN的顯著變化是驅(qū)動(dòng)SOC及其組分發(fā)生改變的最關(guān)鍵因子,近自然化改造是一種有效提高低質(zhì)低效馬尾松人工林土壤團(tuán)聚體穩(wěn)定性,促使森林固持更多更穩(wěn)定的SOC的經(jīng)營(yíng)管理措施。

    參考文獻(xiàn):

    ALLISON SD, WEINTRAUB MN, GARTNER TB, et al.,2010. Evolutionary-economic principles as regulators of soil enzyme production and ecosystem function [M]//Soil enzymology. Heidelberg: Springer: 229-243.

    BAO SD, 2000. Soil agrochemical analysis [M]. Beijing: China Agriculture Press: 1-120.? [鮑士旦, 2000. 土壤農(nóng)化分析 [M]. 北京: 中國(guó)農(nóng)業(yè)出版社: 1-120.]

    BAO Y, GAO Y, ZENG XM, et al., 2018. Relationships between carbon and nitrogen contents and enzyme activities in soil of three typical subtropical forests in China [J]. Chin J Plant Ecol, 42(4): 508-516.? [鮑勇, 高穎, 曾曉敏, 等, 2018. 中亞熱帶3種典型森林土壤碳氮含量和酶活性的關(guān)系 [J]. 植物生態(tài)學(xué)報(bào), 42(4): 508-516.]

    BOSSIO DA, COOK-PATTON SC, ELLIS PW, et al., 2020. The role of soil carbon in natural climate solutions [J]. Nat Sustain, 3(5): 391-398.

    BROCKETT BFT, PRESCOTT CE, GRAYSTON SJ, 2012. Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada [J]. Soil Boil Biochem, 44(1): 9-20.

    CASTRO FILHO C, LOURENO A, GUIMARES MF et al., 2002. Aggregate stability under different soil management systems in a red latosol in the state of Parana, Brazil [J]. Soil Till Res, 65(1): 45-51.

    CHEN H, LI DJ, ZHAO J, et al., 2018. Effects of nitrogen addition on activities of soil nitrogen acquisition enzymes: A meta-analysis [J]. Agric Ecosyst Environ: 126-131.

    CHEN YK, TAN XM, LI M, et al., 2021. The effects of the mixture of the valuable nitrogen-fixing tree species Dalbergia odorifera and second-generation Eucalyptus urophylla on the structure and function of soil microbial community in subtropical China [J]. Guihaia, 41(9): 1476-1485. [陳永康, 譚許脈, 李萌, 等, 2021. 珍貴固氮樹(shù)種降香黃檀與二代巨尾桉混交種植對(duì)土壤微生物群落結(jié)構(gòu)和功能的影響 [J]. 廣西植物, 41(9): 1476-1485.]

    DUNGAIT JAJ, HOPKINS DW, GREGORY AS, et al., 2012. Soil organic matter turnover is governed by accessibility not recalcitrance [J]. Global Change Biol, 18(6): 1781-1796.

    FAN TWM, LANE AN, CHEKMENEV E, et al., 2004. Synthesis and physico-chemical properties of peptides in soil humic substances [J]. J Pept Res, 63(3): 253-264.

    FU LY, ZENG WS, TANG SZ, XIA ZS, et al., 2011. Using linear mixed model and dummy variable model approaches to construct compatible single-tree biomass equations at different scales-A case study for masson pine in Southern China [J]. J For Sci, 41(7): 1547-1554.

    GU W, LI ZA, ZOU B, et al., 2007. Carbon sequestration in soils of rehabilitated plantations on severely eroded lands in tropical China [J]. J Trop Subtrop Bot, 15(5): 369-376.? [顧偉, 李志安, 鄒碧, 等, 2007. 華南熱帶人工林土壤有機(jī)碳含量及其穩(wěn)定性特征 [J]. 熱帶亞熱帶植物學(xué)報(bào), 15(5): 369-376.]

    GUO X, LUO Z, SUN OJ, 2021. Long-term litter type treatments alter soil carbon composition but not microbial carbon utilization in a mixed pine-oak forest [J]. Biogeochemistry, 152(2): 327-343.

    HUANG XM, LIU SR, YOU YM, et al., 2017. Microbial community and associated enzymes activity influence soil carbon chemical composition in Eucalyptus urophylla plantation with mixing N2-fixing species in subtropical China [J]. Plant Soil, 414(1-2): 199-212.

    JANDL R, LINDNER M, VESTERDAL L, et al., 2007. How strongly can forest management influence soil carbon sequestration?? [J]. Geoderma, 137(3-4): 253-268.

    JIN H, SUN OJ, LIU J, 2010. Changes in soil microbial biomass and community structure with addition of contrasting types of plant litter in a semiarid grassland ecosystem [J]. J Plant Ecol, 3(3): 209-217.

    KOOCH Y, SANJI R, TABARI M, 2019. The effect of vegetation change in C and N contents in litter and soil organic fractions of a Northern Iran temperate forest [J]. Catena, 178: 32-39.

    LIANG Y, MING AG, HE YJ, et al., 2021. Structure and function of soil bacterial communities in the monoculture and mixed plantation of Pinus massoniana and Castanopsis hystrix in southern subtropical China [J]. Chin J Appl Ecol, 32(3): 878-886.? [梁艷, 明安剛, 何友均, 等, 2021. 南亞熱帶馬尾松-紅椎混交林及其純林土壤細(xì)菌群落結(jié)構(gòu)與功能 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 32(3): 878-886.]

    LIU SR, YANG YJ, WANG H, 2018. Development strategy and management countermeasures of planted forests in China: transforming from timber-centered single objective management towards multi-purpose management for enhancing quality and benefits of ecosystem services [J]. Acta Ecol Sin, 38(1): 1-10.? [劉世榮, 楊予靜, 王暉, 2018. 中國(guó)人工林經(jīng)營(yíng)發(fā)展戰(zhàn)略與對(duì)策: 從追求木材產(chǎn)量的單一目標(biāo)經(jīng)營(yíng)轉(zhuǎn)向提升生態(tài)系統(tǒng)服務(wù)質(zhì)量和效益的多目標(biāo)經(jīng)營(yíng) [J]. 生態(tài)學(xué)報(bào), 38(1): 1-10.]

    LIU ZL, MING AG, JIA HY, et al., 2017. Effects of close-to-nature transformation on structure characteristics of Pinus massoniana and Cunninghamia lanceolata plantations [J]. J Nanjing For Univ (Nat Sci Ed), 41(4): 101-107.? [劉志龍, 明安剛, 賈宏炎, 等, 2017. 近自然化改造對(duì)桂南馬尾松和杉木人工林結(jié)構(gòu)特征的影響 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 41(4): 101-107.]

    LU LH, GUO WF, CAI DX, et al., 2019. Study on carbon storage of monoculture and mixed plantation of Pinus massoniana and Castanopsis hystrix [J]. J Cent S Univ For Technol, 39(7): 78-84.? [盧立華, 郭文福, 蔡道雄, 等, 2019. 馬尾松與紅椎純林及混交林生態(tài)系統(tǒng)碳儲(chǔ)量研究 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 39(7): 78-84.]

    LU YC, GAN J, 2002. The development trend of forest management in the 21st century [J]. World For Res, 15(1): 1-11.? [陸元昌, 甘敬, 2002. 21世紀(jì)的森林經(jīng)理發(fā)展動(dòng)態(tài) [J]. 世界林業(yè)研究, 15(1): 1-11.]

    LUCAS-BORJA ME, DELGADO-BAQUERIZO M, 2019. Plant diversity and soil stoichiometry regulates the changes in multifunctionality during pine temperate forest secondary succession [J]. Sci Total Environ, 697: 134204.

    MO XQ, XIAO N, TAN XM, et al., 2022. Effects of nitrogen-fixing tree species on soil aggregate-associated enzyme activities and ecoenzymatic stoichiometry in Eucalyptus plantation [J/OL]. Guihaia, https://kns.cnki.net/kcms/detail/45.1134.q.20211231.1215.002.html.? [莫雪青, 肖納, 譚許脈, 等, 2022. 固氮樹(shù)種對(duì)桉樹(shù)人工林土壤團(tuán)聚體酶活性及其化學(xué)計(jì)量比的影響 [J/OL]. 廣西植物, https://kns.cnki.net/kcms/detail/45.1134.q.20211231.1215.002.html.]

    MILTNER A, BOMBACH P, SCHMIDT-BR?CKEN B, et al., 2012. SOM genesis: microbial biomass as a significant source [J]. Biogeochemistry, 111(1): 41-55.

    MING AG, LIU SR, LI H, et al., 2017. Effects of close-to-nature transformation on biomass and its allocation in Pinus massoniana and Cunninghamia lanceolata plantations [J]. Acta Ecol Sin, 37(23): 7833-7842.? [明安剛, 劉世榮, 李華, 等, 2017. 近自然化改造對(duì)馬尾松和杉木人工林生物量及其分配的影響 [J]. 生態(tài)學(xué)報(bào), 37(23): 7833-7842.]

    MING AG, YANG YJ, LIU SR, et al., 2018. Effects of near natural forest management on soil greenhouse gas flux in Pinus massoniana (Lamb.) and Cunninghamia lanceolata (Lamb.) Hook. Plantation [J]. Forests, 9(5): 229.

    MORADI M, MOHADJER MR, SEFIDI K, et al., 2012. Over-mature beech trees (Fagus orientalis Lipsky) and close-to-nature forestry in northern Iran [J]. J For Res, 23(2): 289-294.

    MOU YJ, SHAO JA, GUO Y et al., 2019. The outlook of researches on forest based health and wellness environment construction under the concept of near-nature management [J]. For Econ, 41(8): 49-55.? [牟耀杰, 邵景安, 郭躍, 等, 2019. 近自然經(jīng)營(yíng)理念下森林康養(yǎng)環(huán)境營(yíng)造研究展望 [J]. 林業(yè)經(jīng)濟(jì), 41(8): 49-55.]

    PRESCOTT CE, 2010. Litter decomposition: what controls it and how can we alter it to sequester more carbon in forest soils?? [J]. Biogeochemistry, 101(1): 133-149.

    QIAN J, ZHANG LP, WANG WY, 2018. The relationship between soil aggregates and eroded sediments from sloping vegetated red soils of South China [J]. Acta Ecol Sin, 38(5): 1590-1599.? [錢婧, 張麗萍, 王文艷, 2018. 紅壤坡面土壤團(tuán)聚體特性與侵蝕泥沙的相關(guān)性 [J]. 生態(tài)學(xué)報(bào), 38(5): 1590-1599.]

    SAIYA-CORK KR, SINSABAUGH RL, ZAK DR, 2002. The effects of long-term nitrogen deposition on extracellular enzyme activity in an Acer saccharum forest soil [J]. Soil Boil Biochem, 34(9): 1309-1315.

    SATTLER T, BORCARD D, ARLETTAZ R, et al., 2010. Spider, bee, and bird communities in cities are shaped by environmental control and high stochasticity [J]. Ecology, 91(11): 3343-3353.

    SCHMIDT MWI, TORN MS, ABIVEN S, et al., 2011. Persistence of soil organic matter as an ecosystem property [J]. Nature, 478(7367): 49-56.

    SCHTZ JP, 2011. Development of close to nature forestry and the role of ProSilva Europe [J]. Zbornik Gozdarstva in Lesarstva, 94: 39-42.

    SIMPSON AJ, SIMPSON MJ, SMITH E, et al., 2007. Microbially derived inputs to soil organic matter: are current estimates too low?? [J]. Environ Sci Technol, 41(23): 8070-8076.

    SINSABAUGH RL, LAUBER CL, WEINTRAUB MN, et al., 2008. Stoichiometry of soil enzyme activity at global scale [J]. Ecol Lett, 11(11): 1252-1264.

    SIX J, ELLIOTT ET, PAUSTIAN K, 2000. Soil macroaggregate turnover and microaggregate formation: a mechanism for C sequestration under no-tillage agriculture [J]. Soil Boil Biochem, 32(14): 2099-2103.

    STEINAUER K, TILMAN D, WRAGG PD, et al., 2015. Plant diversity effects on soil microbial functions and enzymes are stronger than warming in a grassland experiment [J]. Ecology, 96(1): 99-112.

    SUN DJ, WEN YG, LUO YH, et al., 2015. Effect of close-to-nature management on species diversity in a Cunninghamia lanceolata plantation [J]. For Res, 28(2): 202-208.? [孫冬婧, 溫遠(yuǎn)光, 羅應(yīng)華, 等, 2015. 近自然化改造對(duì)杉木人工林物種多樣性的影響 [J]. 林業(yè)科學(xué)研究, 28(2): 202-208.]

    VANCE ED, BROOKES PC, JENKINSON DS, 1987. An extraction method for measuring soil microbial biomass C [J]. Soil Boil Biochem, 19(6): 703-707.

    WANG H, LIU SR, WANG JX,et al., 2018. Mixed-species plantation with Pinus massoniana and Castanopsis hystrix accelerates C loss in recalcitrant coniferous litter but slows C loss in labile broadleaf litter in southern China [J]. For Ecol Manag, 422: 207-213.

    YANG YJ, LIU SR, CHEN L, et al., 2018. Responses of soil and aggregate-associated organic carbon mineralization to throughfall reduction in a Pinus massoniana plantation [J]. Chin J Appl Ecol, 29(6): 1779-1786.? [楊予靜, 劉世榮, 陳琳, 等, 2018. 馬尾松人工林土壤和團(tuán)聚體有機(jī)碳礦化對(duì)穿透雨減少的響應(yīng) [J]. 應(yīng)用生態(tài)學(xué)報(bào), 29(6): 1779-1786.]

    YOU YM, HUANG XM, ZHU HG, et al., 2018. Positive interactions between Pinus massoniana and Castanopsis hystrix species in the uneven-aged mixed plantations can produce more ecosystem carbon in subtropical China [J]. For Ecol Manag, 410: 193-200.

    YOU YM, WANG J, HUANG XM, et al., 2014. Relating microbial community structure to functioning in forest soil organic carbon transformation and turnover [J]. Ecol Evol, 4(5): 633-647.

    YOU YM, XU HC, WU XP, et al., 2020. Native broadleaf tree species stimulate topsoil nutrient transformation by changing microbial community composition and physiological function, but not biomass in subtropical plantations with low P status [J]. For Ecol Manag, 477: 118491.

    ZHANG W, HENDRIX PF, DAME LE, et al., 2013. Earthworms facilitate carbon sequestration through unequal amplification of carbon stabilization compared with mineralization [J]. Nat Comm, 4(1): 1-9.

    (責(zé)任編輯李莉)

    猜你喜歡
    土壤酶活性
    不同施肥處理對(duì)植煙土壤關(guān)鍵酶活性及根系活力的影響
    秸稈還田及腐熟劑對(duì)土壤微生物特性和酶活性的影響
    海岸帶森林生態(tài)系統(tǒng)中土壤酶活性研究進(jìn)展
    大棚蘿卜菜田土壤理化性質(zhì)及相關(guān)酶活性特征
    新疆泥火山不同生境土壤微生物數(shù)量、酶活性與理化因子的關(guān)系
    不同施肥年限對(duì)贛南臍橙果園土壤酶活性及微生物種群的影響
    伊犁河谷不同土地利用方式對(duì)土壤養(yǎng)分與酶活性的影響
    重金屬污染對(duì)土壤微生物及土壤酶活性影響的研究進(jìn)展
    遼西半干旱區(qū)果糧間作對(duì)土壤微生物和酶的影響
    精河沙區(qū)土壤酶分布特征及其對(duì)土壤理化性狀的響應(yīng)
    婷婷丁香在线五月| 国产精品精品国产色婷婷| 亚洲一区二区三区色噜噜| 老司机福利观看| 天天躁日日操中文字幕| www日本黄色视频网| 老司机在亚洲福利影院| h日本视频在线播放| 亚洲在线自拍视频| 国产视频内射| 天美传媒精品一区二区| 国产国拍精品亚洲av在线观看 | 精品福利观看| 淫妇啪啪啪对白视频| 久久久久国内视频| 久久99热这里只有精品18| 欧美日韩综合久久久久久 | 精品日产1卡2卡| 偷拍熟女少妇极品色| 国产不卡一卡二| 国产精品98久久久久久宅男小说| 亚洲av五月六月丁香网| 精品国产超薄肉色丝袜足j| 国产亚洲精品一区二区www| 免费观看精品视频网站| 精品久久久久久成人av| 在线观看免费午夜福利视频| 18禁在线播放成人免费| 午夜福利欧美成人| 一区二区三区高清视频在线| 窝窝影院91人妻| 亚洲欧美日韩高清专用| 女生性感内裤真人,穿戴方法视频| 亚洲中文字幕一区二区三区有码在线看| 真人做人爱边吃奶动态| 狂野欧美激情性xxxx| 欧美成人一区二区免费高清观看| 亚洲精品美女久久久久99蜜臀| 亚洲 欧美 日韩 在线 免费| 成人精品一区二区免费| 欧美成狂野欧美在线观看| av国产免费在线观看| 欧美乱色亚洲激情| 国产黄色小视频在线观看| 国产老妇女一区| 欧美成人免费av一区二区三区| 禁无遮挡网站| 乱人视频在线观看| 亚洲成人精品中文字幕电影| 亚洲专区国产一区二区| 黄色丝袜av网址大全| 欧美乱妇无乱码| 日本免费a在线| 热99在线观看视频| 国产乱人视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲中文日韩欧美视频| 人妻丰满熟妇av一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐动态| 亚洲精品久久国产高清桃花| 69人妻影院| 亚洲国产欧美人成| 久久久久亚洲av毛片大全| 国产老妇女一区| av福利片在线观看| 韩国av一区二区三区四区| 亚洲在线自拍视频| 九色国产91popny在线| 欧洲精品卡2卡3卡4卡5卡区| 可以在线观看的亚洲视频| 国产精品久久视频播放| 69av精品久久久久久| 麻豆国产97在线/欧美| 18禁国产床啪视频网站| 村上凉子中文字幕在线| 成人av一区二区三区在线看| 好看av亚洲va欧美ⅴa在| 国产精品 欧美亚洲| 亚洲欧美激情综合另类| 一本久久中文字幕| 久久精品91无色码中文字幕| 国产精品久久久久久精品电影| 国产精品久久久久久久电影 | 国产成+人综合+亚洲专区| 中文字幕人妻熟人妻熟丝袜美 | 又紧又爽又黄一区二区| 脱女人内裤的视频| 亚洲av电影不卡..在线观看| 亚洲欧美日韩高清专用| 亚洲av第一区精品v没综合| 国产精品久久久人人做人人爽| 一进一出抽搐gif免费好疼| 国产毛片a区久久久久| 日本三级黄在线观看| av片东京热男人的天堂| a在线观看视频网站| 亚洲五月天丁香| 天天添夜夜摸| 女生性感内裤真人,穿戴方法视频| av国产免费在线观看| 免费看a级黄色片| 亚洲一区二区三区色噜噜| 国产毛片a区久久久久| 无限看片的www在线观看| 亚洲av成人精品一区久久| 欧美3d第一页| 国产久久久一区二区三区| 亚洲国产欧美网| 国产精品久久久久久久久免 | 国产高潮美女av| 伊人久久精品亚洲午夜| 可以在线观看毛片的网站| 欧美国产日韩亚洲一区| xxxwww97欧美| 1024手机看黄色片| 国产精品久久视频播放| 观看美女的网站| www.999成人在线观看| 国产三级在线视频| 亚洲欧美日韩无卡精品| 国产精品久久久久久精品电影| 亚洲av免费在线观看| 99久久精品热视频| 亚洲av美国av| 麻豆国产97在线/欧美| 亚洲狠狠婷婷综合久久图片| 免费搜索国产男女视频| 欧美在线一区亚洲| 搡老妇女老女人老熟妇| 天天添夜夜摸| 麻豆国产97在线/欧美| 最新中文字幕久久久久| 日韩欧美 国产精品| 手机成人av网站| 午夜视频国产福利| svipshipincom国产片| 国产中年淑女户外野战色| 国产亚洲av嫩草精品影院| 最近最新中文字幕大全免费视频| 国产一区二区亚洲精品在线观看| 精品国产美女av久久久久小说| 国产午夜精品久久久久久一区二区三区 | 国产爱豆传媒在线观看| 国内精品一区二区在线观看| 免费人成在线观看视频色| 法律面前人人平等表现在哪些方面| 国产v大片淫在线免费观看| 婷婷丁香在线五月| 久久草成人影院| av女优亚洲男人天堂| 十八禁人妻一区二区| 亚洲精华国产精华精| 少妇裸体淫交视频免费看高清| 淫秽高清视频在线观看| 久久精品影院6| 欧美黄色片欧美黄色片| 亚洲美女视频黄频| 色噜噜av男人的天堂激情| 亚洲精品日韩av片在线观看 | 一a级毛片在线观看| 国产野战对白在线观看| 99国产极品粉嫩在线观看| 身体一侧抽搐| 每晚都被弄得嗷嗷叫到高潮| 精品福利观看| 久久久久国产精品人妻aⅴ院| 日韩大尺度精品在线看网址| 级片在线观看| 国产精品 欧美亚洲| 亚洲欧美精品综合久久99| 国产精品久久久久久久久免 | 男女视频在线观看网站免费| 午夜福利在线观看吧| 一本一本综合久久| 真人一进一出gif抽搐免费| 一本一本综合久久| 欧美丝袜亚洲另类 | 天堂av国产一区二区熟女人妻| 91九色精品人成在线观看| 内射极品少妇av片p| 欧美日韩精品网址| 亚洲人成网站在线播放欧美日韩| 精品人妻偷拍中文字幕| 欧美日本亚洲视频在线播放| 丰满的人妻完整版| 亚洲 国产 在线| 成人特级黄色片久久久久久久| 亚洲成人免费电影在线观看| 欧美一级a爱片免费观看看| 精品久久久久久久人妻蜜臀av| 99热6这里只有精品| 亚洲av不卡在线观看| av黄色大香蕉| 怎么达到女性高潮| 国产亚洲精品综合一区在线观看| 18禁裸乳无遮挡免费网站照片| 国产伦人伦偷精品视频| 欧美激情在线99| 麻豆国产av国片精品| 欧美色欧美亚洲另类二区| 久久草成人影院| 国产精品一区二区免费欧美| 久久久久免费精品人妻一区二区| 网址你懂的国产日韩在线| 成人特级黄色片久久久久久久| 伊人久久大香线蕉亚洲五| 久久国产乱子伦精品免费另类| 国产精品久久久久久精品电影| 十八禁网站免费在线| 偷拍熟女少妇极品色| 国产高清有码在线观看视频| 高清毛片免费观看视频网站| 亚洲欧美精品综合久久99| 一本精品99久久精品77| 中文字幕高清在线视频| av在线蜜桃| 五月玫瑰六月丁香| 51午夜福利影视在线观看| 51国产日韩欧美| 国产欧美日韩精品一区二区| 日本黄色视频三级网站网址| eeuss影院久久| 亚洲av成人av| 国产伦人伦偷精品视频| 99riav亚洲国产免费| 黄色片一级片一级黄色片| 特大巨黑吊av在线直播| 亚洲在线自拍视频| 老司机深夜福利视频在线观看| 日本免费一区二区三区高清不卡| 亚洲精品成人久久久久久| 久久久国产精品麻豆| 女警被强在线播放| 中文在线观看免费www的网站| 免费无遮挡裸体视频| 又粗又爽又猛毛片免费看| 国产一区二区在线av高清观看| 禁无遮挡网站| 女人被狂操c到高潮| 国产在视频线在精品| 人人妻人人看人人澡| 激情在线观看视频在线高清| 嫩草影视91久久| 露出奶头的视频| bbb黄色大片| 久久久久久久亚洲中文字幕 | 亚洲在线观看片| 九九在线视频观看精品| 国产精品av视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产在线观看| 欧美精品啪啪一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲国产欧美人成| 亚洲第一欧美日韩一区二区三区| 窝窝影院91人妻| 99久久九九国产精品国产免费| 欧美午夜高清在线| 国产精品爽爽va在线观看网站| 18禁黄网站禁片午夜丰满| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲| 嫩草影院入口| 人妻丰满熟妇av一区二区三区| 午夜福利在线在线| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 网址你懂的国产日韩在线| 九九久久精品国产亚洲av麻豆| 男女视频在线观看网站免费| 法律面前人人平等表现在哪些方面| 少妇熟女aⅴ在线视频| 日本五十路高清| 亚洲性夜色夜夜综合| aaaaa片日本免费| 香蕉久久夜色| 亚洲美女黄片视频| 欧美性感艳星| 两人在一起打扑克的视频| 国产熟女xx| 国产av在哪里看| av在线天堂中文字幕| www国产在线视频色| 国内精品一区二区在线观看| 久久这里只有精品中国| 欧美xxxx黑人xx丫x性爽| 久久人妻av系列| 国产极品精品免费视频能看的| 波多野结衣高清作品| xxxwww97欧美| 少妇熟女aⅴ在线视频| 综合色av麻豆| 美女被艹到高潮喷水动态| 亚洲激情在线av| 欧美成人免费av一区二区三区| 日韩欧美国产一区二区入口| 久久午夜亚洲精品久久| 一级黄色大片毛片| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看| 亚洲av一区综合| 99久久成人亚洲精品观看| 99久久久亚洲精品蜜臀av| 国产一区二区三区在线臀色熟女| 女人被狂操c到高潮| 激情在线观看视频在线高清| 宅男免费午夜| 亚洲成人久久爱视频| 99久久精品一区二区三区| 国产高清视频在线播放一区| av在线天堂中文字幕| 久久精品国产99精品国产亚洲性色| 亚洲va日本ⅴa欧美va伊人久久| 欧美av亚洲av综合av国产av| 欧美黑人巨大hd| 中文字幕人妻熟人妻熟丝袜美 | 两个人视频免费观看高清| 国产精品亚洲一级av第二区| 日本黄色片子视频| 成年版毛片免费区| 天堂√8在线中文| 午夜久久久久精精品| 国产不卡一卡二| 性欧美人与动物交配| 亚洲av成人不卡在线观看播放网| 黄色片一级片一级黄色片| 日本黄大片高清| 日韩高清综合在线| 中文亚洲av片在线观看爽| 天天躁日日操中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 性色av乱码一区二区三区2| 一级毛片高清免费大全| 亚洲无线在线观看| 精品久久久久久久人妻蜜臀av| 亚洲av成人av| 免费人成在线观看视频色| 欧美黄色片欧美黄色片| 男女那种视频在线观看| 久久久成人免费电影| 精品99又大又爽又粗少妇毛片 | 国产高潮美女av| 毛片女人毛片| av在线蜜桃| 给我免费播放毛片高清在线观看| 精品一区二区三区视频在线观看免费| 禁无遮挡网站| 国产一级毛片七仙女欲春2| 18美女黄网站色大片免费观看| 老司机在亚洲福利影院| 中亚洲国语对白在线视频| 欧美激情久久久久久爽电影| 国产精品久久久久久人妻精品电影| 欧美+日韩+精品| 中文字幕熟女人妻在线| 好看av亚洲va欧美ⅴa在| 中文字幕熟女人妻在线| 好看av亚洲va欧美ⅴa在| 舔av片在线| 看黄色毛片网站| 丰满乱子伦码专区| 91在线观看av| 观看美女的网站| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 日韩欧美在线乱码| 香蕉丝袜av| 久久精品亚洲精品国产色婷小说| 观看美女的网站| 特级一级黄色大片| 亚洲欧美日韩卡通动漫| 手机成人av网站| 亚洲国产精品成人综合色| 国产精品一区二区三区四区免费观看 | 超碰av人人做人人爽久久 | 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| bbb黄色大片| 在线a可以看的网站| 亚洲无线在线观看| 亚洲天堂国产精品一区在线| av天堂在线播放| 亚洲精品久久国产高清桃花| 精品久久久久久久人妻蜜臀av| 婷婷六月久久综合丁香| 热99re8久久精品国产| 亚洲国产精品合色在线| 亚洲内射少妇av| 午夜久久久久精精品| 午夜福利在线在线| 最近最新免费中文字幕在线| 成人精品一区二区免费| 3wmmmm亚洲av在线观看| tocl精华| 久久久久免费精品人妻一区二区| 欧美午夜高清在线| 日韩有码中文字幕| av天堂中文字幕网| 丰满乱子伦码专区| 国产av一区在线观看免费| 日本与韩国留学比较| 亚洲av免费高清在线观看| 国产毛片a区久久久久| 国产亚洲欧美在线一区二区| 亚洲不卡免费看| 久久欧美精品欧美久久欧美| 午夜免费激情av| 少妇熟女aⅴ在线视频| 九九热线精品视视频播放| 日本一本二区三区精品| av国产免费在线观看| bbb黄色大片| 成人特级黄色片久久久久久久| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 亚洲av美国av| 亚洲国产中文字幕在线视频| 人人妻人人澡欧美一区二区| 欧美日韩一级在线毛片| 999久久久精品免费观看国产| 亚洲精品影视一区二区三区av| 国产精品99久久久久久久久| 久久亚洲真实| 亚洲 欧美 日韩 在线 免费| 国产美女午夜福利| 久久国产乱子伦精品免费另类| 亚洲欧美精品综合久久99| 久久亚洲精品不卡| 亚洲av成人精品一区久久| 亚洲内射少妇av| 色综合婷婷激情| 在线观看美女被高潮喷水网站 | 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 国产午夜福利久久久久久| 精品熟女少妇八av免费久了| 亚洲乱码一区二区免费版| 欧美成人免费av一区二区三区| 久久精品国产亚洲av香蕉五月| 久久人妻av系列| 亚洲内射少妇av| 亚洲午夜理论影院| 91久久精品国产一区二区成人 | 亚洲欧美日韩卡通动漫| 国产老妇女一区| 亚洲五月婷婷丁香| 久久精品国产综合久久久| 美女大奶头视频| 亚洲人与动物交配视频| 欧美日韩黄片免| 成人无遮挡网站| 国产免费男女视频| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 国产一区二区亚洲精品在线观看| 日韩大尺度精品在线看网址| 亚洲av免费在线观看| 国内久久婷婷六月综合欲色啪| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 成人av一区二区三区在线看| 久久久色成人| 欧美日韩黄片免| 18禁裸乳无遮挡免费网站照片| 手机成人av网站| 18禁美女被吸乳视频| 91九色精品人成在线观看| 成人特级av手机在线观看| 免费在线观看日本一区| 久久精品影院6| 国产精品一区二区免费欧美| 99热只有精品国产| 亚洲人成网站高清观看| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 18禁美女被吸乳视频| 久久香蕉精品热| 黄色片一级片一级黄色片| 国产真人三级小视频在线观看| 99久久成人亚洲精品观看| 床上黄色一级片| 又黄又粗又硬又大视频| 国产成人影院久久av| 嫩草影视91久久| 欧美日韩亚洲国产一区二区在线观看| 白带黄色成豆腐渣| 麻豆成人午夜福利视频| 亚洲精品日韩av片在线观看 | 亚洲精品成人久久久久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产一区二区三区在线臀色熟女| 老熟妇乱子伦视频在线观看| 国产欧美日韩一区二区精品| 欧美av亚洲av综合av国产av| 男女那种视频在线观看| 深爱激情五月婷婷| 九九在线视频观看精品| 老司机午夜十八禁免费视频| 色综合亚洲欧美另类图片| 国产伦精品一区二区三区四那| 亚洲av第一区精品v没综合| 夜夜看夜夜爽夜夜摸| 久久久久性生活片| 欧美午夜高清在线| 亚洲午夜理论影院| 成人无遮挡网站| 欧美+亚洲+日韩+国产| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出| 免费看光身美女| 亚洲成av人片免费观看| 亚洲国产精品sss在线观看| 久久国产精品人妻蜜桃| 天天添夜夜摸| 熟女电影av网| 国产亚洲欧美98| 国产私拍福利视频在线观看| 韩国av一区二区三区四区| 色精品久久人妻99蜜桃| 好男人电影高清在线观看| 内地一区二区视频在线| 久久久久久久午夜电影| 婷婷精品国产亚洲av| 国产精品美女特级片免费视频播放器| 最新中文字幕久久久久| 国产激情偷乱视频一区二区| 老鸭窝网址在线观看| 亚洲专区中文字幕在线| 日本在线视频免费播放| 精品久久久久久久毛片微露脸| 免费av不卡在线播放| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月| 真人做人爱边吃奶动态| 中亚洲国语对白在线视频| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 90打野战视频偷拍视频| 小说图片视频综合网站| 97超级碰碰碰精品色视频在线观看| 国产免费av片在线观看野外av| 天堂av国产一区二区熟女人妻| 99久久无色码亚洲精品果冻| 伊人久久大香线蕉亚洲五| 欧美国产日韩亚洲一区| 中文字幕熟女人妻在线| 亚洲成人免费电影在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品一及| 国产伦一二天堂av在线观看| 欧美乱妇无乱码| www日本黄色视频网| 国产精品电影一区二区三区| 老鸭窝网址在线观看| 波多野结衣高清无吗| 久久久久精品国产欧美久久久| 久久久精品欧美日韩精品| 国产亚洲精品久久久久久毛片| 日本撒尿小便嘘嘘汇集6| 国产精品免费一区二区三区在线| 国产精品亚洲美女久久久| 亚洲av免费高清在线观看| 啦啦啦免费观看视频1| 国产欧美日韩精品一区二区| 亚洲精华国产精华精| 女人高潮潮喷娇喘18禁视频| 国产在线精品亚洲第一网站| 性色av乱码一区二区三区2| 每晚都被弄得嗷嗷叫到高潮| 国产久久久一区二区三区| 国产一区二区三区在线臀色熟女| 全区人妻精品视频| 国产精品久久久久久久久免 | 18禁美女被吸乳视频| 日韩欧美 国产精品| 国产久久久一区二区三区| 亚洲狠狠婷婷综合久久图片| 精品国产超薄肉色丝袜足j| 人妻久久中文字幕网| 亚洲无线在线观看| 亚洲国产高清在线一区二区三| 亚洲第一欧美日韩一区二区三区| 国产精品三级大全| 小说图片视频综合网站| 欧美中文日本在线观看视频| 国产欧美日韩精品亚洲av| 亚洲av二区三区四区| 夜夜夜夜夜久久久久| АⅤ资源中文在线天堂| 黄色成人免费大全| 国产97色在线日韩免费| 久久精品影院6| 变态另类成人亚洲欧美熟女| 欧美黄色淫秽网站| 亚洲精品粉嫩美女一区| 麻豆成人av在线观看| 国产成人欧美在线观看| 国产精品三级大全| 九九在线视频观看精品| 欧美一级毛片孕妇| 淫秽高清视频在线观看| 老熟妇乱子伦视频在线观看| 成年人黄色毛片网站| 日韩欧美国产一区二区入口| 男女那种视频在线观看| 免费人成视频x8x8入口观看| 中文字幕av成人在线电影| 九色国产91popny在线| 久久久久久久精品吃奶| 国产av在哪里看| 观看免费一级毛片| 国产精品爽爽va在线观看网站| 国产精品一区二区三区四区免费观看 | 久久精品国产99精品国产亚洲性色|