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

    等碳量添加生物炭和秸稈對(duì)煙田土壤呼吸及凈碳收支的影響

    2023-03-23 07:22:50符云鵬李耀鑫柳淵博金佳威蔣偉峰劉詠艷
    作物學(xué)報(bào) 2023年5期
    關(guān)鍵詞:煙田土壤溫度土壤水分

    符云鵬 劉 天 李耀鑫 柳淵博 王 靜 金佳威 蔣偉峰 劉詠艷 楊 洋 云 菲

    等碳量添加生物炭和秸稈對(duì)煙田土壤呼吸及凈碳收支的影響

    符云鵬 劉 天 李耀鑫 柳淵博 王 靜 金佳威 蔣偉峰 劉詠艷 楊 洋 云 菲*

    河南農(nóng)業(yè)大學(xué)煙草學(xué)院 / 煙草行業(yè)煙草栽培重點(diǎn)實(shí)驗(yàn)室, 河南鄭州 450002

    研究生物炭和秸稈添加對(duì)煙田生態(tài)系統(tǒng)土壤呼吸及碳收支的影響, 闡明生物炭、秸稈以及二者配施處理煙田的固碳效應(yīng)。2020—2021年, 設(shè)置常規(guī)施肥(CK)、常規(guī)施肥+2.25 t hm–2生物炭-C (T1)、常規(guī)施肥+2.25 t hm–2秸稈-C (T2)、常規(guī)施肥+1.125 t hm–2生物炭-C+1.125 t hm–2秸稈-C (T3) 4個(gè)處理, 對(duì)不同組分土壤呼吸及土壤主要環(huán)境因子、土壤碳增量、作物凈初級(jí)生產(chǎn)力固碳量以及進(jìn)行農(nóng)業(yè)生產(chǎn)造成的碳排放進(jìn)行了測(cè)定。結(jié)果表明, 添加秸稈、秸稈與生物炭配施處理烤煙生育期內(nèi)土壤呼吸累計(jì)排放碳量顯著高于對(duì)照(<0.05), 提升幅度分別為21.40%~35.45%、5.90%~9.89%; 添加生物炭、秸稈處理土壤自養(yǎng)呼吸占比分別較對(duì)照提高6.35%~7.34%、3.21%~5.97%, 生物炭與秸稈配施處理僅2021年較對(duì)照提升3.91%。添加生物炭、秸稈提高了生育期內(nèi)土壤溫度和水分, 單施生物炭處理土壤水分顯著提高了1.93%~7.07%。土壤主要環(huán)境因子中, 土壤溫度對(duì)土壤呼吸速率影響最大, 二者呈極顯著正相關(guān); 添加秸稈、生物炭與秸稈配施處理土壤呼吸速率與土壤濕度呈顯著正相關(guān)。各處理生態(tài)系統(tǒng)固碳量均為正值, 表示為“碳匯”, 其中添加生物炭不僅可增加凈初級(jí)生產(chǎn)力(NPP)固碳量以及土壤固碳量, 還可降低煙草生長(zhǎng)季內(nèi)土壤呼吸累積碳排放量, “碳匯”能力最強(qiáng)。因此, 施用生物炭是降低煙田土壤呼吸碳排放, 增強(qiáng)煙田生態(tài)系統(tǒng)“碳匯”能力的最佳途徑。

    生物炭; 秸稈; 土壤呼吸; 碳收支

    溫室氣體排放導(dǎo)致全球氣候變暖問題日益嚴(yán)峻,各種極端天氣頻頻出現(xiàn)。中國(guó)CO2排放總量位居世界第一[1], 我國(guó)實(shí)現(xiàn)2030年“碳達(dá)峰”目標(biāo)與2060年“碳中和”愿景挑戰(zhàn)巨大。土壤呼吸是陸地生態(tài)系統(tǒng)到大氣的第二大碳通量, 陸地土壤CO2年平均排放通量約為80.4 PgC[2]。全球土壤儲(chǔ)存了至少2倍于大氣的碳[3-4], 土壤呼吸向大氣中排放的CO2通量正逐年增加[5]。土壤呼吸由自養(yǎng)呼吸(植物根際及根際微生物)和異養(yǎng)呼吸(非根際相關(guān)的土壤微生物和土壤動(dòng)物)組成。1987—1999年間, 全球土壤呼吸速率以每年0.161 PgC的速率增加[6]。1990—2014年間, 異養(yǎng)呼吸(heterotrophic respiration Rh)與總呼吸(soil respiration, Rs)的比值由0.54增加到0.63, 異養(yǎng)呼吸的增加可能是響應(yīng)環(huán)境變暖帶來的環(huán)境變化[7]。農(nóng)田耕地土壤退化問題日益嚴(yán)峻, 在提高土壤肥力的同時(shí)增強(qiáng)土壤碳庫(kù)貯存已被公認(rèn)為全球范圍內(nèi)應(yīng)對(duì)減緩?fù)寥劳嘶皻夂蜃兣呗訹8-9], 對(duì)于占全球土壤近33%的有機(jī)碳損失退化土壤來言尤其重要[10]。農(nóng)田土壤應(yīng)用秸稈已被公認(rèn)為是一種環(huán)境友好的做法,可增加有機(jī)碳的輸入, 提高礦質(zhì)元素的有效性進(jìn)而阻止土壤退化, 提高土壤質(zhì)量及作物產(chǎn)量, 同時(shí)還可減少化肥的施用[11-12]。有研究表明, 秸稈還田會(huì)增加土壤CO2的排放, 加劇全球變暖問題[12]?;诮斩挼葟U棄生物質(zhì)制備而成的生物炭被認(rèn)為是另一種增加土壤碳儲(chǔ)量的有效方式, 并可起到增加土壤養(yǎng)分含量及土壤碳庫(kù)穩(wěn)定性的作用[13]。研究表明, 生物炭在土壤中周轉(zhuǎn)率極低(0.0046% d-1)[14], 因此可將生物炭施于土壤實(shí)現(xiàn)碳封存進(jìn)而達(dá)到緩解溫室效應(yīng)的目的。農(nóng)田添加生物炭可通過改變土壤物理結(jié)構(gòu)、水分含量、陽(yáng)離子交換能力、酶活性及土壤微生物群落結(jié)構(gòu)進(jìn)而減少土壤呼吸CO2排放[15]。秸稈與生物炭的交互作用研究較少。二者配合施用對(duì)農(nóng)田生態(tài)系統(tǒng)凈碳收支的影響尚不清楚。煙草作為一種壟作作物, 農(nóng)事操作較為復(fù)雜, 涉及到施肥起壟、中耕培土、除草等, 對(duì)土壤造成擾動(dòng)較多, 勢(shì)必影響到煙田土壤呼吸。目前, 關(guān)于農(nóng)田生態(tài)系統(tǒng)碳收支的研究主要集中在糧食作物(玉米、小麥、水稻), 而有關(guān)煙田生態(tài)系統(tǒng)碳收支研究相對(duì)較少。煙草作為一種重要的經(jīng)濟(jì)作物, 在我國(guó)種植面積較大。秸稈還田是煙區(qū)常用的土壤改良措施, 生物炭作為新興的土壤改良劑近幾年在植煙土壤改良中的應(yīng)用日漸增多, 然而關(guān)于生物炭和秸稈添加對(duì)煙田生態(tài)系統(tǒng)碳收支影響的研究較少。因此, 研究不同碳源輸入下煙田土壤主要環(huán)境因子特征及土壤呼吸組分變化規(guī)律、煙田生態(tài)系統(tǒng)凈碳收支, 對(duì)正確評(píng)價(jià)煙田生態(tài)系統(tǒng)碳源/匯效應(yīng)具有重要的理論意義及實(shí)用價(jià)值。

    為避免土壤環(huán)境不同對(duì)準(zhǔn)確區(qū)分土壤呼吸組分造成的影響, 本試驗(yàn)采用壟溝法對(duì)土壤呼吸組分進(jìn)行了分離。通過田間原位連續(xù)監(jiān)測(cè), 探究不同碳源輸入對(duì)煙田土壤環(huán)境及不同呼吸組分的影響, 對(duì)生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力、土壤固碳量、土壤呼吸作用排放碳量、物料投入碳量及農(nóng)業(yè)投入碳排放量進(jìn)行收集測(cè)定, 定量評(píng)價(jià)生物炭、秸稈及二者混施處理生態(tài)系統(tǒng)的固碳潛力。

    1 材料與方法

    1.1 研究地點(diǎn)及試驗(yàn)材料

    定位試驗(yàn)開始于2017年, 按照“等碳量”的原則連年添加生物炭和秸稈。本試驗(yàn)于2020年4月至9月、2021年4月至9月在河南農(nóng)業(yè)大學(xué)許昌校區(qū)現(xiàn)代煙草科教園區(qū)(東經(jīng)113°85′, 北緯34°14′)進(jìn)行。該區(qū)域?qū)儆诘湫偷呐瘻貛啙駶?rùn)季風(fēng)氣候, 年平均氣溫15.6℃, 多年平均降雨量約為700 mm。土壤類型為褐土, 質(zhì)地為中壤。土壤基礎(chǔ)養(yǎng)分含量(質(zhì)量分?jǐn)?shù))為: 有機(jī)質(zhì)19.09 g kg-1、堿解氮74.7 mg kg-1、速效磷8.7 mg kg-1、速效鉀114.5 mg kg-1和水溶性氯24.8 mg kg-1, pH 7.61。試驗(yàn)所用花生殼生物炭由河南省生物炭技術(shù)工程實(shí)驗(yàn)室提供, 在400~450℃下限氧熱解30 min制備而成, 其基本性質(zhì)為總碳(total carbon, TC) 47.28%、總氮(total nitrogen, TN) 1.25%、總硫(total sulfur, TS) 0.48%、總磷(total phosphorus, TP) 1.31%、總鉀(total potassium, TK) 1.77%、碳氮比(carbon nitrogen ratio, C/N) 37.82、酸堿度(pH) 8.5。秸稈為腐熟小麥秸稈, 基本性質(zhì)為TC 44.82%、TN 1.21%、TS 0.80%、TP 0.72%、TK 0.96%。

    1.2 試驗(yàn)設(shè)計(jì)

    試驗(yàn)設(shè)4個(gè)處理, 隨機(jī)區(qū)組排列, 重復(fù)3次, 小區(qū)面積0.02 hm2。各處理分別為: CK: 常規(guī)施肥; T1: 常規(guī)施肥+2.25 t hm-2生物炭-C; T2: 常規(guī)施肥+2.25 t hm-2秸稈-C; T3: 常規(guī)施肥+1.125 t hm–2生物炭-C+ 1.125 t hm–2秸稈-C。生物炭、秸稈均于烤煙移栽前2周施入土壤, 與表層0~20 cm土層充分混合后起壟,壟距1.2 m。各處理氮用量為30 kg hm–2, 氮磷鉀比例1∶2∶4?;视蔁煵輰S脧?fù)合肥(15% N、15% P2O5、15% K2O), 硫酸鉀(52% K2O)和過磷酸鈣(12% P2O5)提供, 追肥為硝酸鉀(13.5% N、44.5% K2O), 基肥、追肥比例為4∶1。種植品種為中煙100, 行距1.2 m, 株距0.55 m, 密度為15,000株 hm-2, 2020年5月3日移栽, 9月27日采收結(jié)束, 2021年5月6日移栽, 9月17日采收結(jié)束。

    1.3 測(cè)定指標(biāo)與方法

    1.3.1 土壤呼吸及土壤環(huán)境因子 土壤呼吸使用Li-8100土壤碳通量測(cè)量?jī)x(LI-COR, Lincoln, 美國(guó))測(cè)定, 每次測(cè)量前24 h清除地表雜草。于測(cè)定日的08:30—11:30之間進(jìn)行土壤呼吸測(cè)定, 代表本日土壤呼吸平均值, 為剔除土壤呼吸日變化對(duì)測(cè)量結(jié)果的影響, 每次測(cè)量時(shí)均按照不同的順序?qū)τ^測(cè)點(diǎn)進(jìn)行測(cè)定, 每5 d測(cè)量1次, 雨天除外。每個(gè)小區(qū)共放置3個(gè)PVC環(huán), 其中壟間和壟體各安裝1個(gè)PVC環(huán)。安裝好的PVC環(huán)高出地面3 cm, PVC環(huán)在生育期內(nèi)不再移動(dòng), 直至測(cè)量結(jié)束。使用壕溝法將土壤總呼吸(total soil respiration, Rs)劃分為土壤自養(yǎng)呼吸(soil autotrophic respiration, Ra)和土壤異養(yǎng)呼吸(soil heterotrophic respiration, Rh)。在每個(gè)處理的壟間設(shè)置一個(gè)矩形區(qū)域, 在矩形區(qū)域邊界向內(nèi)5 cm左右挖一個(gè)60 cm深的溝, 保留長(zhǎng)寬約60 cm×40 cm的未擾動(dòng)區(qū)域。沿溝槽內(nèi)壁放置孔隙為0.038 mm的尼龍網(wǎng), 然后將土壤回填。在未干擾區(qū)域中放置一個(gè)PVC環(huán)測(cè)量土壤呼吸, 此測(cè)定值為異養(yǎng)呼吸(Rh), 將壟體呼吸與壟間呼吸的平均值視為田間土壤總呼吸(Rs), 土壤總呼吸(Rs)與異養(yǎng)呼吸(Rh)的差值為自養(yǎng)呼吸(Ra)。土壤溫度、濕度采用WET三參儀(WET-2, 英國(guó))與土壤呼吸一同測(cè)定, 測(cè)定深度為壟體土壤5 cm深處。

    土壤呼吸CO2日排放量:

    R=R×3600×24×12×10–5(1)

    整個(gè)生長(zhǎng)季CO2累計(jì)排放量:

    式中,R為當(dāng)日土壤呼吸CO2日排放量(kg hm–2),s為土壤呼吸速率(μmol m–2s–1), 12為CO2-C的摩爾質(zhì)量(g mol–1), 3600、24、10–5均為轉(zhuǎn)換系數(shù);r為整個(gè)生長(zhǎng)季CO2累計(jì)排放量(kg hm–2),RR1分別為第次和第+1次測(cè)定時(shí)CO2日排放量(kg hm–2),為相鄰2次土壤呼吸測(cè)定間隔天數(shù),first為第1次測(cè)定的土壤呼吸CO2排放量(kg hm–2)。

    1.3.2 煙草生物量的測(cè)定 烤煙采收期, 各小區(qū)選擇3株長(zhǎng)勢(shì)均勻、代表性強(qiáng)的煙株, 采集整株植物樣品, 分為根、莖、葉3個(gè)部分, 對(duì)生育期內(nèi)進(jìn)行的農(nóng)事操作如打頂抹杈、打底腳葉去除的植株生物量以及煙株凋落物進(jìn)行收集, 于105℃殺青30 min后在65℃烘干至恒重, 冷卻后測(cè)定各器官干物質(zhì)質(zhì)量。

    1.3.3 土壤有機(jī)碳含量 移栽前及采收后采用“S”形取樣法采集0~20 cm耕層土壤樣品。自然風(fēng)干后過0.15 mm篩, 采用德國(guó)Elementar公司生產(chǎn)的有機(jī)碳分析儀(VarioMACRO TOC)測(cè)定土壤有機(jī)碳含量。

    1.3.4 農(nóng)田生態(tài)系統(tǒng)碳收支估算 農(nóng)田生態(tài)系統(tǒng)凈碳收支估算參考田冬[16]方法進(jìn)行。采用全生命周期評(píng)價(jià)系統(tǒng)凈固碳效應(yīng), 農(nóng)田生態(tài)系統(tǒng)的源匯效應(yīng)取決于土壤碳儲(chǔ)量增量、作物凈初級(jí)生產(chǎn)力固定碳量、土壤呼吸排放碳量、物料還田碳量、農(nóng)業(yè)生產(chǎn)物資消耗及人工管理排放碳平衡進(jìn)行估算。

    (1) 土壤碳儲(chǔ)量增量

    ΔSOC=(SOCi-SOC0)×BD×H×10–1。

    ΔSOC為土壤碳增量(t hm–2), SOCi為第次測(cè)定的土壤有機(jī)碳含量(g kg–1), SOC0試驗(yàn)開始前基礎(chǔ)土樣土壤有機(jī)碳含量(g kg–1), BD為土壤容重(g cm–3), H為土層深度。

    (2) 系統(tǒng)碳吸收(Ca)

    Ca=CNPP=Ccrop+Clitter

    CNPP為系統(tǒng)內(nèi)作物的凈初級(jí)生產(chǎn)力固定碳, Ccrop為作物固碳(包括作物地上部的莖稈、葉、煙杈、花及地下部的根系所包含的碳), Clitter為田間凋落物殘留碳。

    (3) 系統(tǒng)碳排放(Ce)

    Ce=Cr+Ceh

    Cr為系統(tǒng)土壤呼吸碳排放, Ceh為生產(chǎn)活動(dòng)所產(chǎn)生的碳排放。

    (4) 系統(tǒng)凈碳匯(Cs)

    系統(tǒng)凈碳匯(Cs)通過土壤固碳增量、吸收與排放碳平衡計(jì)算:

    Cs=Ci+CNPP+ΔSOC-Cr-Ceh

    Ci為物料碳還田量, CNPP為凈初級(jí)生產(chǎn)力, ΔSOC為土壤固碳量, Cr為土壤呼吸碳排放量, Ceh為采用碳足跡法, 從作物播種到收獲期農(nóng)田生產(chǎn)活動(dòng)中的農(nóng)業(yè)投入碳排放量, 主要包括農(nóng)藥、化肥、耕地起壟、灌溉、田間管理如打頂抹杈及農(nóng)機(jī)使用過程中產(chǎn)生的碳排放(表1)。表1中煙苗碳排放系數(shù)無(wú)參考文獻(xiàn), 本試驗(yàn)采用德國(guó)Elementar公司生產(chǎn)的全自動(dòng)CNS元素分析儀(VarioMACRO cube)測(cè)定移栽所用煙苗碳積累量。

    1.3.5 相關(guān)性 采用經(jīng)驗(yàn)指數(shù)方程描述土壤呼吸與表層土壤溫度的關(guān)系:=e,為土壤呼吸速率(μmol m–2s–1),為地表下5 cm處土壤溫度,、為擬合參數(shù)。

    10: 溫度敏感系數(shù), 表示溫度每增加10℃, 土壤呼吸增加的倍數(shù),10=e10。

    土壤濕度與土壤呼吸速率之間的關(guān)系采用二次回歸方程進(jìn)行擬合:=++2, 式中,為土壤呼吸速率,為地表下5 cm處土壤體積含水量,、、均為方程擬合參數(shù)。

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

    采用Microsoft Excel 2010、SPSS 22.0進(jìn)行數(shù)據(jù)處理及分析, 用Duncan’s新復(fù)極差法進(jìn)行差異顯著性檢驗(yàn), 使用Origin 2018和Microsoft Excel 2010制作圖表。

    2 結(jié)果與分析

    2.1 添加生物炭和秸稈對(duì)土壤溫度和濕度的影響

    由圖1可以看出, 2020年與2021年烤煙大田生育期壟體土壤溫度、濕度變化趨勢(shì)基本一致。生育期內(nèi), 各處理土壤溫度均呈先上升后下降的趨勢(shì), 最高溫度出現(xiàn)在8月初。2020年土壤溫度變化范圍在21.92~39.72℃之間, 添加生物炭、秸稈及生物炭和秸稈配施處理平均土壤溫度分別較CK提高0.55%、0.41%、0.78%, 處理間無(wú)顯著差異。2021年土壤溫度整體低于2020年, 變化幅度范圍為25.02~31.77℃, 各處理平均土壤溫度分別較CK提高0.63%、0.52%、0.58%, 處理間差異不顯著。整個(gè)生育期內(nèi), 土壤水分變化波動(dòng)較為明顯, 主要受農(nóng)田灌溉及降雨的影響; 2020年T1、T3處理生育期內(nèi)土壤水分均值分別較CK顯著提高7.07%、2.03%, 2021年T1、T2、T3處理土壤水分均值較CK分別提高1.93%、0.01%、3.38%, T1、T3處理土壤水分含量均值與CK差異顯著(<0.05)。

    2.2 添加生物炭和秸稈對(duì)煙草生育期內(nèi)土壤呼吸速率及土壤CO2-C累積排放量的影響

    由圖2可知, 在煙草生育期內(nèi), 各處理土壤總呼吸表現(xiàn)出明顯的季節(jié)性, 2020年與2021年土壤呼吸速率變化趨勢(shì)基本一致, 總體呈先增加后降低, 與土壤溫度變化趨勢(shì)相似, 8月初達(dá)到峰值。2020年各處理土壤呼吸速率變化范圍分別為1.66~4.76、1.93~4.01、2.32~5.94、2.00~4.61 μmol m–2s–1, 土壤呼吸平均速率分別為2.87、2.74、3.82、3.01 μmol m–2s–1, 處理間差異均達(dá)顯著水平(<0.05), 其中T1處理平均土壤呼吸速率較CK顯著降低4.53%, T2、T3處理分別較CK顯著增加33.40%、4.88%。2021年各處理土壤呼吸速率變化規(guī)律相同, 土壤呼吸速率變化范圍分別為1.77~6.22、2.12~5.54、2.79~7.49、2.55~6.84 μmol m–2s–1, 土壤呼吸平均速率分別為3.96、3.75、4.82、4.36 μmol m–2s–1, 均高于2020年相對(duì)應(yīng)處理, 且處理間差異顯著(<0.05), 其中T1處理較CK顯著降低5.30%, T2、T3處理分別較CK顯著增加21.72%、10.10%??梢? 添加生物炭可顯著降低煙草生育期內(nèi)土壤呼吸平均速率。

    表1 各項(xiàng)農(nóng)業(yè)物資的碳排放系數(shù)

    圖1 生物炭和秸稈添加下土壤溫度、濕度的動(dòng)態(tài)變化

    CK: 對(duì)照; T1: 常規(guī)施肥+2.25 t hm–2生物炭-C; T2: 常規(guī)施肥+2.25 t hm–2秸稈-C; T3: 常規(guī)施肥+1.125 t hm–2生物炭-C+1.125 t hm–2秸稈-C。

    CK: the control; T1: the conventional fertilization + 2.25 t hm–2biochar-C; T2: the conventional fertilization + 2.25 t hm–2straw-C; T3: the conventional fertilization + 1.125 t hm–2biochar-C + 1.125 t hm–2straw-C.

    2020年與2021年土壤異養(yǎng)呼吸平均速率表現(xiàn)為T2>T3>CK>T1, T2處理顯著高于其他處理(<0.05), 2年T1處理土壤異養(yǎng)呼吸平均速率分別較CK顯著降低9.28%、10.66%。2020年各處理土壤自養(yǎng)呼吸平均速率表現(xiàn)為T2>T1>T3>CK, T2與其他處理差異顯著(<0.05); CK、T1、T2、T3處理土壤自養(yǎng)呼吸占總呼吸比例分別為43.99%、46.87%、45.40%、42.45%, 其中T1、T2處理較CK提高了自養(yǎng)呼吸在總呼吸的占比。2021年各處理土壤自養(yǎng)呼吸平均速率表現(xiàn)為T2>T3>T1>CK, 與CK相比, 各處理均提高了自養(yǎng)呼吸占比。

    由表2可知, 2020年煙草生育期內(nèi)各處理土壤呼吸作用CO2-C累積排放量在3833.65~5395.87kg hm–2之間, 2021年土壤呼吸作用CO2-C累積排放量高于2020年, 在5356.42~6916.35 kg hm–2之間, 2年均表現(xiàn)為T2>T3>CK>T1, 且T2、T3處理顯著高于CK。

    2.3 添加生物炭和秸稈對(duì)煙草生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力的影響

    由圖3可知, 煙株各器官固碳量均表現(xiàn)為葉>莖> 根>杈>花>凋落物; 2020年CK、T1、T2、T3處理葉片固碳量占總固碳量的比例分別為41.47%、41.03%、42.62%、40.75%, 2021年各處理葉片固碳量占總固碳量的比例分別為41.63%、40.91%、41.71%、41.18%。2020年, 各處理葉片固碳量表現(xiàn)為T2>T3>T1>CK, 處理間差異顯著; 莖和根系固碳量表現(xiàn)為T3>T2>T1>CK, 其中T3處理顯著高于其他處理, T2、T1處理顯著高于CK; 其他器官固碳量表現(xiàn)為T3>T1>T2>CK, 處理間差異顯著。2021年, 根系固碳量表現(xiàn)為T3、T1顯著高于CK, 各處理葉、莖和其他器官固碳量均顯著高于CK。

    圖2 生物炭和秸稈添加下煙田土壤呼吸速率的動(dòng)態(tài)變化

    處理同圖1。Treatments are the same as those given in Fig. 1.

    表2 生物炭和秸稈添加下土壤呼吸CO2-C累計(jì)排放量

    同行不同小寫字母表示處理間差異顯著(< 0.05)。處理同圖1。

    Different lowercase letters in the same row indicate significant differences between treatments at the 0.05 probability level. Treatments are the same as those given in Fig. 1.

    圖3 添加生物炭和秸稈對(duì)煙株各部位固碳量的影響

    同顏色橫柱不同小寫字母表示處理間差異顯著(< 0.05)。處理同圖1。

    Different lowercase letters in the same color horizontal column indicate significant differences at the 0.05 probability level between treatments. Treatments are the same as those given in Fig. 1.

    由圖4可知, 2020年與2021年各處理煙株固碳量均表現(xiàn)為T3>T2>T1>CK, 2020年T1、T2、T3處理煙株固碳量分別較CK顯著增加14.32%、16.59%、16.63%, 2021年各處理分別較對(duì)照顯著增加15.29%、16.73%、18.66%。

    圖4 添加生物炭和秸稈對(duì)煙田生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力的影響

    不同小寫字母表示同一年份處理間差異顯著(< 0.05)。處理同圖1。

    Different lowercase letters indicate significant differences between treatments in same year at the 0.05 probability level. Treatments are the same as those given in Fig. 1.

    2.4 土壤呼吸對(duì)土壤溫度、濕度的響應(yīng)

    土壤溫度、濕度是影響土壤呼吸的兩大主要因素。由表3可知, 各處理土壤溫度與土壤呼吸呈指數(shù)相關(guān)。2020年, T1和T3處理Q10值分別較CK降低11.45%、13.86%, T2處理較CK增加1.20%。2021年, 各處理Q10值表現(xiàn)為CK>T3>T2>T1, 與2020年相比, 除T1外, 各處理Q10值均有所增加。整體來看, 單施生物炭可有效降低土壤呼吸敏感系數(shù)。2020年, 土壤濕度與土壤呼吸擬合度較低,2變化范圍在0.032~0.168之間, 其中T3處理擬合度最高; T2、T3處理土壤濕度與土壤呼吸速率為顯著正相關(guān)(<0.05)。2021年, 各處理土壤呼吸與土壤濕度擬合度均大幅度增高, 這可能與當(dāng)年降雨量較大有關(guān)。

    2.5 煙田生態(tài)系統(tǒng)碳收支

    由表4可知, CK表現(xiàn)為 “碳源”, 各處理均呈“碳匯”。處理間凈固碳量表現(xiàn)為T1>T3>T2。2021年各處理“碳匯”效應(yīng)與2020年相比較弱, 這與2021年土壤呼吸作用CO2-C排放量大幅度增加有關(guān)。相較于CK, 各處理增加了碳投入量, 提高了凈初級(jí)生產(chǎn)力固碳量及土壤固碳量。綜合來看, 單施生物炭的T1處理“碳匯”能力最強(qiáng)。添加不同碳源的3個(gè)處理凈初級(jí)生產(chǎn)力固碳量無(wú)顯著差異, 但均顯著高于CK (<0.05)。2020年T1處理土壤固碳量較CK、T2、T3分別提高11.18、0.24、0.05倍, 土壤呼吸CO2-C排放量分別較其他處理降低3.76%、28.95%、9.12%。2021年T1處理土壤固碳量較CK、T2、T3分別提高10.88、0.20、0.04倍, 土壤呼吸CO2-C排放量較其他處理降低5.98%、22.55%、14.44%。因此, 欲提高煙田生態(tài)系統(tǒng)固碳量, 應(yīng)從減少土壤呼吸碳排放、增加凈初級(jí)生產(chǎn)力固碳量及土壤固碳量入手。

    表3 土壤呼吸速率與土壤環(huán)境因子的擬合方程

    處理同圖1。Treatments are the same as those given in Fig. 1.

    表4 生物炭和秸稈添加對(duì)煙田生態(tài)系統(tǒng)碳收支的影響

    處理同圖1。Treatments are the same as those given in Fig. 1.

    3 討論

    3.1 不同碳源輸入對(duì)土壤呼吸特征及其組分的影響

    本研究發(fā)現(xiàn), 烤煙生長(zhǎng)前期土壤呼吸速率較低, 8月初土壤呼吸CO2排放達(dá)到峰值, 隨后下降(圖2-a), 這與何甜甜等[22]、李亞森等[23]在豫中煙田生態(tài)系統(tǒng)研究烤煙生育期內(nèi)土壤呼吸速率的動(dòng)態(tài)變化結(jié)果一致。本研究結(jié)果表明, 添加秸稈顯著增加了煙草生育期內(nèi)土壤呼吸平均速率, 而添加生物炭顯著降低了土壤呼吸平均速率, 生物炭與秸稈配施較單施秸稈降低了土壤呼吸速率。添加生物炭降低土壤呼吸平均速率可能是因?yàn)? (1) 少量存在于生物炭表面易被微生物分解利用的碳源消耗殆盡后, 剩下穩(wěn)定性極強(qiáng)的以芳香環(huán)不規(guī)則疊層存在的碳不易被分解[24]; (2) 生物炭促進(jìn)土壤團(tuán)聚體生成及穩(wěn)定有利于農(nóng)田土壤碳固持[25]; (3) 生物炭較秸稈可顯著增加一些礦化能力弱、腐殖化能力強(qiáng)的細(xì)菌豐度[26], 表現(xiàn)為土壤呼吸速率降低。劉杏認(rèn)等[27]在冬小麥-夏玉米輪作體系中研究認(rèn)為添加玉米秸稈生物炭可提高輪作周期內(nèi)CO2排放通量, 這可能歸因于生物炭用量、還田年限的差異。許多研究表明高量生物炭添加下可大幅度提高土壤呼吸速率[23,28], 本研究中生物炭添加量較小(4.76 t hm–2), 表現(xiàn)為土壤呼吸速率降低。此外, 生物炭施用年限同樣會(huì)對(duì)其作用效果產(chǎn)生影響, Singh等[29]認(rèn)為由桉木、桉葉、家禽糞便、牛糞制備而成的生物炭添加于沙姜土壤的前二三年起正激發(fā)效應(yīng), 促進(jìn)本土有機(jī)碳礦化, 之后一直到第5年正激發(fā)效應(yīng)逐漸減弱, 對(duì)本土有機(jī)碳礦化作用由促進(jìn)轉(zhuǎn)為抑制。本研究為定位試驗(yàn), 已連續(xù)添加5年生物炭, 導(dǎo)致了土壤呼吸速率的降低。

    明確土壤呼吸各組分對(duì)土壤總呼吸的貢獻(xiàn), 對(duì)定量評(píng)價(jià)生態(tài)系統(tǒng)碳收支具有較大的理論及實(shí)踐價(jià)值[30]。自養(yǎng)呼吸在土壤總呼吸中占比在不同生態(tài)系統(tǒng)中差異較大。本研究中, 2020年烤煙生育期內(nèi)各處理自養(yǎng)呼吸占比在42.45%~46.87%之間, 2021年各處理自養(yǎng)呼吸占比在43.44%~46.63%之間。田冬等[31]研究表明, 油菜季自養(yǎng)呼吸在土壤總呼吸中占比為17.31%~22.9%, 玉米季為26.49%~53.65%。Kuzyakov等[32]采用14C標(biāo)定示蹤法研究自養(yǎng)呼吸對(duì)土壤總呼吸的貢獻(xiàn)率, 結(jié)果顯示自養(yǎng)呼吸貢獻(xiàn)率在15%~60%之間, 均值為41%。魯齊飛等[33]在豫中煙田生態(tài)系統(tǒng)研究各土壤呼吸組分在土壤總呼吸中占比, 其中自養(yǎng)呼吸占比為40.26%~43.43%。本研究中自養(yǎng)呼吸占比相對(duì)較高, 這可能與土壤呼吸組分的區(qū)分測(cè)定方法有關(guān), 本試驗(yàn)采用壕溝法, 對(duì)土壤微生態(tài)影響相對(duì)較小, 裸地處無(wú)植被覆蓋, 土壤溫度較高, 會(huì)導(dǎo)致異養(yǎng)呼吸測(cè)定值偏高。單施生物炭處理和單施秸稈處理較對(duì)照提高了土壤自養(yǎng)呼吸占比, 主要是添加生物炭、秸稈可提高土壤有機(jī)碳含量, 改善土壤物理特性, 促進(jìn)根系發(fā)育, 進(jìn)而提高了土壤自養(yǎng)呼吸占比。此外, 添加生物炭對(duì)原土有機(jī)碳礦化的抑制作用導(dǎo)致異養(yǎng)呼吸速率降低[34]。

    3.2 不同碳源輸入下土壤環(huán)境因素對(duì)土壤呼吸的影響

    土壤水分和溫度是影響土壤CO2排放的兩大土壤環(huán)境因子[35]。本研究中, 土壤溫度與土壤呼吸呈極顯著正相關(guān), 是對(duì)土壤呼吸速率影響最大的土壤環(huán)境因子, 對(duì)土壤呼吸變化的貢獻(xiàn)率達(dá)到30.7%~ 63.15%, 趙亞麗等[36]將玉米、小麥秸稈添加于粉質(zhì)黏壤土取得研究結(jié)果與本項(xiàng)研究結(jié)果一致。土壤CO2排放溫度敏感系數(shù)(Q10)是預(yù)測(cè)土壤碳排放的重要指標(biāo), 本研究中添加生物炭處理較對(duì)照降低Q10。有研究發(fā)現(xiàn)生物炭施于農(nóng)田短期內(nèi)增加Q10, 長(zhǎng)期則降低Q10[37], 本試驗(yàn)連續(xù)5年施用生物炭可能導(dǎo)致Q10降低。添加秸稈刺激土壤中微生物活化, 加速土壤呼吸速率, Q10隨著底物增加而增加[38]。

    土壤水分對(duì)土壤呼吸的影響具有復(fù)雜性和不確定性, 同時(shí)也有可能被其他因子所掩蓋[39]。有研究表明, 土壤CO2排放量與水分存在顯著的負(fù)相關(guān)關(guān)系[40]。本研究中, 2020年單施秸稈處理、秸稈與生物炭配施處理土壤水分與呼吸速率之間呈顯著正相關(guān)。2021年添加生物炭、秸稈處理土壤呼吸速率與土壤水分呈顯著正相關(guān), 土壤水分含量增大可大幅度提高生育期內(nèi)煙田土壤呼吸速率。這可能與該年度7月至8月份降雨較多導(dǎo)致土壤含水量較高有關(guān)。有研究表明[41]土壤水分含量大于12% (質(zhì)量含水率)時(shí), 土壤呼吸速率隨土壤水分含量增加而增大, 土壤水分含量增大時(shí), 土壤微生物對(duì)有機(jī)碳的分解能力增強(qiáng)。Singh等[42]研究發(fā)現(xiàn), 土壤呼吸速率隨土壤水分含量提高而增大, 土壤EOC濃度隨之降低。風(fēng)干土壤中EOC濃度始終較高, 表示由于水分限制, 有機(jī)碳的微生物礦化作用較低。Keiluweit等[43]認(rèn)為土壤水分含量過高時(shí)導(dǎo)致O2的擴(kuò)散減少, 從而限制好氧微生物呼吸。本項(xiàng)試驗(yàn)中抑制土壤呼吸作用的土壤水分含量上限有待進(jìn)一步研究。

    3.3 不同碳源輸入對(duì)煙田生態(tài)系統(tǒng)碳收支的影響

    煙田生態(tài)系統(tǒng)碳收支取決于碳輸入及碳支出。Liu等[44]在中國(guó)8個(gè)作物主要產(chǎn)區(qū)進(jìn)行田間試驗(yàn)發(fā)現(xiàn), 生物炭添加不會(huì)通過增加農(nóng)業(yè)土壤中CO2排放碳損失來抵消凈碳積累。本研究中, 添加生物炭后土壤碳累計(jì)排放量較對(duì)照降低, 與Liu等[44]研究結(jié)果一致。單施秸稈處理及秸稈與生物炭配施處理顯著增加了烤煙生育期內(nèi)土壤呼吸累計(jì)排放碳量。添加不同碳源處理均較對(duì)照顯著提高了煙田生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力, 其中生物炭與秸稈配施處理較對(duì)照提升效果最大; 各處理煙株葉片固碳量最高, 在40.75%~42.62%之間。生育期內(nèi)各處理土壤固碳量均高于對(duì)照, 這與張鵬等[45]、Tang等[27]研究結(jié)果一致; 碳輸入項(xiàng)中煙株固碳量均顯著高于土壤固碳量。煙草田間生產(chǎn)管理措施較多, 包括施肥、起壟、中耕培土、打頂抹杈等, 農(nóng)事操作投入造成的碳排放較多, 在286.38~288.73 kg hm-2之間。

    綜合比較農(nóng)田生態(tài)系統(tǒng)碳收支各項(xiàng)指標(biāo), 各處理均增加了凈初級(jí)生產(chǎn)力(NPP)固碳量, 同時(shí)增加了土壤固碳量; 各項(xiàng)碳源物料添加均能提高煙田生態(tài)系統(tǒng)凈固碳量, 表現(xiàn)為“碳匯”。劉青麗等[46]認(rèn)為煙田生態(tài)系統(tǒng)在無(wú)外源碳輸入情況下呈弱“碳匯”, 與本研究結(jié)果不同, 這是由于土壤呼吸測(cè)定方法及碳收支計(jì)算方式不同所致。本項(xiàng)試驗(yàn)采用壟溝法測(cè)定土壤呼吸, 劉青麗等[46]采用靜態(tài)箱-氣象色譜法測(cè)定, 將植株罩入采樣箱, 煙株進(jìn)行光合作用同時(shí)吸收部分土壤呼吸排放CO2, 土壤呼吸速率偏低。此外,本項(xiàng)研究采用全生命周期評(píng)價(jià)煙田生態(tài)系統(tǒng)凈固碳效應(yīng), 對(duì)進(jìn)行各項(xiàng)農(nóng)事操作造成的碳排放進(jìn)行收集, 因而整個(gè)系統(tǒng)碳排放相對(duì)較高。本項(xiàng)研究中, 添加生物炭處理系統(tǒng)凈固碳量最高, “碳匯”能力最強(qiáng)。Tang等[27]認(rèn)為添加生物炭可有效提高作物生物量和土壤碳貯存, 施用生物炭是增強(qiáng)煙田“碳匯”能力的有效途徑, 與本項(xiàng)研究結(jié)果一致。生物炭添加不僅能提高土壤固碳量及凈初級(jí)生產(chǎn)力, 還可降低土壤呼吸作用造成的碳損失, 單施秸稈、生物炭與秸稈混施處理土壤呼吸作用造成的CO2排放顯著高于單施生物炭處理。因此, 增加煙田生態(tài)系統(tǒng)碳匯量首先應(yīng)減少土壤呼吸作用碳排放, 其次增加土壤固碳量和作物凈初級(jí)生產(chǎn)力。

    4 結(jié)論

    添加生物炭可降低土壤呼吸速率及生育期內(nèi)CO2總排放量, 單施秸稈、生物炭與秸稈配施顯著促進(jìn)土壤呼吸CO2-C排放; 單施生物炭、單施秸稈處理較對(duì)照可提高土壤自養(yǎng)呼吸貢獻(xiàn)率。添加生物炭、秸稈對(duì)土壤溫度影響不大, 單施秸稈及生物炭與秸稈配施可顯著提高土壤水分含量。土壤主要環(huán)境因子中, 土壤溫度對(duì)土壤呼吸速率影響最大。添加生物炭、秸稈及二者配施均提高了土壤固碳量、凈初級(jí)生產(chǎn)力固碳量。添加生物炭可降低生育期內(nèi)土壤呼吸累計(jì)CO2-C排放量。施用生物炭是提高煙田生態(tài)系統(tǒng)“碳匯”能力的最佳途徑。

    [1] 2000–2020年全球主要地區(qū)和國(guó)家CO2排放量. 世界石油工業(yè), 2021, 28(4): 81–82.

    CO2emissions of major regions and countries in the world from 2000 to 2020., 2021, 28(4): 81–82 (in Chinese).

    [2] Raich J W, Potter C S, Bhagawati D. Interannual variability in global soil respiration, 1980–1994., 2002, 8: 800–812.

    [3] Kochy M, Hiederer R, Freibauer A. Global distribution of soil organic carbon—Part 1: Masses and frequency distributions of SOC stocks for the tropics, permafrost regions, wetlands, and the world., 2015, 1: 351–365.

    [4] Scharlemann J P W, Tanner E V J, Hiederer R, Kapos V. Global soil carbon: understanding and managing the largest terrestrial carbon pool., 2014, 5: 81–91.

    [5] Bond-Lamberty B, Thomson A. Temperature-associated increases in the global soil respiration record., 2010, 464: 579–582.

    [6] Lei J S, Guo X, Zeng Y F, Zhou J Z, Gao Q, Yang Y F. Temporal changes in global soil respiration since 1987., 2021, 12: 1801.

    [7] Bond-Lamberty B, Bailey V L, Chen M, Gough C M, Vargas RGlobally rising soil heterotrophic respiration over recent decades., 2018, 560:80–83.

    [8] Amelung W, Bossio D, de Vries W, Kogel-knabner I, Lehmann J, Amundson R, Bol R, Collins C, Lal R, Leifeld J, Minasny B, Pan G, Paustian K, Rumpel C, Sanderman J, van Groenigen J W, Mooney S, van Wesemael B, Wander M, Chabb A. Towards a global-scale soil climate mitigation strategy., 2020, 11: 5427.

    [9] Schmidt M W I, Torn M S, Abiven S, Dittmar T, Guggenberger G, Janssens I A, Kleber M, K?gel-Knabner I, Lehmann J, Manning D A C, Nannipieri P, Rasse D P, Weiner S, Trumbore S E. Persistence of soil organic matter as an ecosystem property., 2011, 478: 49–56.

    [10] FAO. Recarbonization of Global Soils: a Dynamic Response to Offset Global Emissions. 2019.

    [11] Lin Y X, Ye G P, Kuzyakov Y, Liu D Y, Fan J B, Ding W X. Long-term manure application increases soil organic matter and aggregation, and alters microbial community structure and keystone taxa., 2019, 51: 187–196.

    [12] Bhattacharyya P, Roy K S, Neogi S, Adhya T K, Rao K S, Manna M C. Effects of rice straw and nitrogen fertilization on greenhouse gas emissions and carbon storage in tropical flooded soil planted with rice., 2012, 32: 119–130.

    [13] Singh B P, Cowie A L. Long-term influence of biochar on native organic carbon mineralisation in a low-carbon clayey soil., 2014, 4: 3687.

    [14] Wang J Y, Xiong Z Q, Kuzyakov Y. Biochar stability in soil: meta-analysis of decomposition and priming effects., 2016, 8: 512–523.

    [15] Lehmann J, Rillig M C, Thies J, Masiello C A, Hockaday W C, Crowley D. Biochar effects on soil biota: a review., 2011, 43: 1812–1836.

    [16] 田冬. 農(nóng)田生態(tài)系統(tǒng)碳收支對(duì)秸稈與生物炭還田的響應(yīng). 西南大學(xué)碩士學(xué)位論文, 重慶, 2017.

    Tian D. Response of Carbon Budget of Farmland Ecosystem to Straw and Biochar Returning. MS Thesis of Southwest University, Chongqing, China, 2017 (in Chinese with English abstract).

    [17] 逯非, 王效科, 韓冰, 歐陽(yáng)志云, 段曉男, 鄭華. 中國(guó)農(nóng)田施用化學(xué)氮肥的固碳潛力及其有效性評(píng)價(jià). 應(yīng)用生態(tài)學(xué)報(bào), 2008, 19: 2239–2250.

    Lu F, Wang X K, Han B, Ouyang Z Y, Duan X N, Zheng H. Carbon sequestration potential and effectiveness evaluation of chemical nitrogen fertilizer in farmland of China., 2008, 19: 2239–2250 (in Chinese with English abstract).

    [18] Lal R. Agricultural activities and the global carbon cycle., 2004, 70: 103–116.

    [19] Lal R. Soil carbon sequestration impacts on global climate change and food security., 2004, 304: 1623–1627.

    [20] West T O, Marland G. A synthesis of carbon sequestration, carbon emissions, and net carbon fiux in agriculture: comparing tillage practices in the United States., 2002, 91: 217–232.

    [21] 劉巽浩, 徐文修, 李增嘉, 褚慶全, 楊曉琳, 陳阜. 農(nóng)田生態(tài)系統(tǒng)碳足跡法:誤區(qū)、改進(jìn)與應(yīng)用: 兼析中國(guó)集約農(nóng)作碳效率(續(xù)). 中國(guó)農(nóng)業(yè)資源與區(qū)劃, 2014, 35(1): 1–7.

    Liu X H, Xu W X, Li Z J, Chu Q Q, Yang X L, Chen F. Farmland ecosystem carbon footprint method: misunderstanding, improvement and application: also an analysis of carbon efficiency of intensive farming in China (continued), 2014, 35(1): 1–7 (in Chinese).

    [22] 何甜甜, 王靜, 符云鵬, 符新妍, 劉天, 李亞坤, 李建華. 等碳量添加秸稈和生物炭對(duì)土壤呼吸及微生物生物量碳氮的影響. 環(huán)境科學(xué), 2021, 42: 450–458.

    He T T, Wang J, Fu Y P, Fu X Y, Liu T, Li Y K, Li J H. Effects of equal carbon addition of straw and biochar on soil respiration and microbial biomass carbon and nitrogen., 2021, 42: 450–458 (in Chinese with English abstract).

    [23] 李亞森, 丁松爽, 殷全玉, 李佳軼, 周迪, 劉國(guó)順. 多年施用生物炭對(duì)河南烤煙種植區(qū)土壤呼吸的影響. 環(huán)境科學(xué), 2019, 40: 915–923.

    Li Y S, Ding S S, Yin Q Y, Li J Y, Zhou D, Liu G S. Effects of multi-year application of biochar on soil respiration in Henan flue-cured tobacco planting area., 2019, 40: 915–923 (in Chinese with English abstract).

    [24] 徐敏, 伍鈞, 張小洪, 楊剛. 生物炭施用的固碳減排潛力及農(nóng)田效應(yīng). 生態(tài)學(xué)報(bào), 2018, 38: 393–404.

    Xu M, Wu J, Zhang X H, Yang G. Carbon sequestration and emission reduction potential of biochar application and farmland effect., 2018, 38: 393–404 (in Chinese with English abstract).

    [25] Huang R, Lan M L, Liu J, Gao M. Soil aggregate and organic carbon distribution at dry land soil and paddy soil: the role of different straws returning.,2017, 24: 27942–27952.

    [26] Tian X P, Wang L, Hou Y H, Wang H, Tsang V F, Wu J H. Responses of soil microbial community structure and activity to incorporation of straws and straw biochars and their effects on soil respiration and soil organic carbon turnover., 2019, 29: 492–503.

    [27] Tang Y, Gao W C, Cai K, Chen Y, Li C B, Lee X Q, Cheng H G, Zhang Q H, Cheng J Z. Effects of biochar amendment on soil carbon dioxide emission and carbon budget in the karst region of southwest China., 2021, 55: 114895.

    [28] 劉杏認(rèn), 張星, 張晴雯, 李貴春, 張慶忠. 施用生物炭和秸稈還田對(duì)華北農(nóng)田CO2、N2O排放的影響. 生態(tài)學(xué)報(bào), 2017, 37: 6700–6711.

    Liu X R, Zhang X, Zhang Q W, Li G C, Zhang Q Z. Effects of biochar application and straw returning on CO2and N2O emissions from farmland in North China, 2017, 37: 6700–6711 (in Chinese with English abstract).

    [29] Singh B P, Cowie A L. Long-term influence of biochar on native organic carbon mineralisation in a low-carbon clayey soil., 2014, 4: 3687.

    [30] Pangle R E, Seiler J. Influence of seedling roots, environmental factors and soil characteristics on soil CO2efflux rates in a 2-year-old loblolly pine (L.) plantation in the Virginia Piedmont., 2002, 116: S85–S96.

    [31] 田冬, 高明, 黃容, 呂盛, 徐暢. 油菜/玉米輪作農(nóng)田土壤呼吸和異養(yǎng)呼吸對(duì)秸稈與生物炭還田的響應(yīng). 環(huán)境科學(xué), 2017, 38: 2988–2999.

    Tian D, Gao M, Huang R, Lyu S, Xu C. Response of soil respiration and heterotrophic respiration to straw and biochar returning in rape/maize rotation farmland., 2017, 38: 2988–2999 (in Chinese with English abstract).

    [32] Kuzyakov Y, Ehrensberger H, Stahr K. Carbon partitioning and below-ground translocation by Lolium perenne., 2001, 33: 61–74.

    [33] 魯琪飛, 葉協(xié)鋒, 韓金, 潘昊東, 張明杰, 王晶, 楊佳豪, 姚鵬偉, 李雪利. 添加有機(jī)物料對(duì)豫中煙田土壤呼吸的影響. 環(huán)境科學(xué), 2022, 43: 3825–3834.

    Lu Q F, Ye X F, Han J, Pan H D, Zhang M J, Wang J, Yang J H, Yao P W, Li X L. Effects of adding organic materials on soil respiration in tobacco fields in central Henan., 2022, 43: 3825–3834 (in Chinese with English abstract).

    [34] Jones D L, Murphy D V, Khalid M, Ahmad W, Edwards-Jones G, DeLuca T H. Short-term biochar-induced increase in soil CO2release is both biotically and abiotically mediated., 2011, 43: 1723–1731.

    [35] Davidson E A, Janssens I A, Luo Y Q. On the variability of respiration in terrestrial ecosystems: moving beyond Q10., 2006, 12: 154–164.

    [36] 趙亞麗, 薛志偉, 郭海斌, 穆心愿, 李潮海. 耕作方式與秸稈還田對(duì)土壤呼吸的影響及機(jī)理. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(19): 155–165.

    Zhao Y L, Xue Z W, Guo H B, Guo H B, Li C H. Effects of tillage methods and straw returning on soil respiration and its mechanism., 2014, 30(19): 155–165 (in Chinese with English abstract).

    [37] Wang X J, Chen G H, Wang S Y, Zhang L Y, Zhang R D. Temperature sensitivity of different soil carbon pools under biochar addition., 2019, 26: 4130–4140.

    [38] 李艾蒙. 外源玉米秸稈碳對(duì)土壤有機(jī)碳激發(fā)效應(yīng)和溫度敏感性的影響. 沈陽(yáng)農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 遼寧沈陽(yáng), 2020.

    Li A M. Effects of Exogenous Corn Straw Carbon on Soil Organic Carbon Excitation Effect and Temperature Sensitivity. MS Thesis of Shenyang Agricultural University, Shenyang, Liaoning, China, 2020 (in Chinese with English abstract).

    [39] 王新源, 李玉霖, 趙學(xué)勇, 毛偉, 崔奪, 曲浩, 連杰, 羅永清. 干旱半干旱區(qū)不同環(huán)境因素對(duì)土壤呼吸影響研究進(jìn)展. 生態(tài)學(xué)報(bào), 2012, 32: 4890–4901.

    Wang X Y, Li Y L, Zhao X Y, Mao W, Cui D, Qu H, Lian J, Luo Y Q. Research progress on the effects of different environmental factors on soil respiration in arid and semi-arid areas., 2012, 32: 4890–4901 (in Chinese with English abstract).

    [40] Lu N, Liu X R, Du Z L, Wang Y D, Zhang Q Z. Effect of biochar on soil respiration in the maize growing season after 5 years of consecutive application., 2014, 38: 505–512.

    [41] 張丁辰, 蔡典雄, 代快, 馮宗會(huì), 張曉明, 王小彬. 旱作農(nóng)田不同耕作土壤呼吸及其對(duì)水熱因子的響應(yīng). 生態(tài)學(xué)報(bào), 2013, 33: 1916–1925.

    Zhang D C, Cai D X, Dai K, Feng Z H, Zhang X M, Wang X B. Soil respiration under different tillage and its response to hydrothermal factors in dryland farmland., 2013, 33: 1916–1925 (in Chinese with English abstract).

    [42] Singh S, Jagadamma S, Liang J Y, Kivlin S N, Mayes M A. Differential organic carbon mineralization responses to soil moisture in three different soil orders under mixed forested system., 2021, 9: 682450.

    [43] Keiluweit M, Nico P S, Kleber M, Fendorf S. Are oxygen limitations under recognized regulators of organic carbon turnover in upland soils?, 2016, 127: 157–171.

    [44] Liu X Y, Zheng J F, Zhang D X, Cheng K, Zhou H M, Zhang A F, Li L Q, Joseph S, Smith P, Crowley D, Kuzyakov Y, Pan G X. Biochar has no effect on soil respiration across Chinese agricultural soils., 2016, 45: 259–265.

    [45] 張鵬, 李涵, 賈志寬, 王維, 路文濤, 張惠, 楊寶平. 秸稈還田對(duì)寧南旱區(qū)土壤有機(jī)碳含量及土壤碳礦化的影響. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2011, 30: 2518–2525.

    Zhang P, Li H, Jia Z K, Wang W, Lu W T, Zhang H, Yang B P. Effect of straw on soil organic carbon content in dryland., 2011, 30: 2518–2525 (in Chinese with English abstract).

    [46] 劉青麗, 蔣雨洲, 鄒焱, 張?jiān)瀑F, 張恒, 石俊雄, 李志宏. 煙田生態(tài)系統(tǒng)碳收支研究. 作物學(xué)報(bào), 2020, 46: 1258–1265.

    Liu Q L, Jiang Y Z, Zou Y, Zhang Y G, Zhang H, Shi J X, Li Z H. Study on carbon budget of tobacco field ecosystem., 2020, 46: 1258–1265 (in Chinese with English abstract).

    Effects of adding biochar and straw with equal carbon content on soil respiration and net carbon budget in tobacco field

    FU Yun-Peng, LIU Tian, LI Yao-Xin, LIU Yuan-Bo, WANG Jing, JIN Jia-Wei, JIANG Wei-Feng, LIU Yong-Yan, YANG Yang, and YUN Fei*

    College of Tobacco Science, Henan Agricultural University / Key Laboratory of Tobacco Cultivation of China Tobacco Industry, Zhengzhou 450002, Henan, China

    To clarify the carbon sequestration effect of biochar, straw and their combined application in tobacco field, the effects of biochar and straw addition on soil respiration and carbon budget of tobacco ecosystem were explored. From 2020 to 2021, four treatments including conventional fertilization (CK), conventional fertilization +2.25 t hm–2biochar-C (T1), conventional fertilization +2.25 t hm–2straw-C (T2), and conventional fertilization +1.125 t hm–2biochar-C + 1.125 t hm–2straw-C (T3) were conducted to measure soil respiration and main environmental factors of soil, soil carbon increment, carbon sequestration of net primary productivity of crops and carbon emissions caused by agricultural production. The results showed that the cumulative carbon emissions from soil respiration during the growth period of flue-cured tobacco treated with straw, straw and biochar were significantly higher than that of the control (< 0.05), and the increase ranges were 21.40%–35.45% and 5.90%–9.89%, respectively. Compared with the control, the proportion of soil autotrophic respiration increased by 6.35%–7.34% and 3.21%–5.97% in the treatment of adding biochar and straw, respectively. The combined application of biochar and straw increased by 3.91% in 2021 compared with the control. The addition of biochar and straw increased the soil temperature and moisture during the growth period, and soil water content increased significantly by 1.93% to 7.07% under biochar treatment alone. Among the main environmental factors in soil, soil temperature had the greatest impact on soil respiration rate. Soil respiration rate was significantly positively correlated with soil temperature. There was a significant positive correlation between soil respiration rate and soil moisture in the treatment of adding straw, biochar, and straw. The carbon sequestration of ecosystem treated by all exogenous additives was positive, which was expressed as “carbon sink”. The addition of biochar could not only increase the carbon sequestration of net primary productivity (NPP) and soil carbon sequestration, but also reduce the cumulative carbon emissions of soil respiration in tobacco growing season, and the “carbon sink” ability was the strongest. Therefore, the application of biochar was the best way to reduce the carbon emission of soil respiration, which could enhance the “carbon sink” capacity of tobacco ecosystem.

    biochar; straw; soil respiration; carbon budget

    10.3724/SP.J.1006.2023.24054

    本研究由中國(guó)煙草總公司河南省公司科技項(xiàng)目(2020410000270020), 河南省煙草公司濟(jì)源市公司項(xiàng)目(2020410881240045)和河南省自然科學(xué)基金青年科學(xué)基金項(xiàng)目(212300410160)資助。

    This study was supported by the Science and Technology Project of China National Tobacco Corporation Henan Province (2020410000270020), the Jiyuan Company Project of Henan Tobacco Company (2020410881240045), and the Youth Science Fund Project of Henan Natural Science Foundation (212300410160).

    云菲, E-mail: yunfeifei55@126.com

    E-mail: yunpengfu@henau.edu.cn

    2022-03-10;

    2022-10-10;

    2022-10-24.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20221021.1638.006.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    煙田土壤溫度土壤水分
    煙田施肥起埂機(jī)的設(shè)計(jì)分析
    遼東山區(qū)3種人工林土壤呼吸對(duì)土壤溫度和土壤水分的響應(yīng)
    西藏高原土壤水分遙感監(jiān)測(cè)方法研究
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    四川省煙田空心蓮子草危害調(diào)查及化學(xué)防除
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    不同覆蓋措施對(duì)棗園土壤水分和溫度的影響
    基于模糊綜合評(píng)價(jià)法的煙田水利設(shè)施經(jīng)濟(jì)評(píng)價(jià)
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    土壤水分的遙感監(jiān)測(cè)方法概述
    久久久精品94久久精品| 少妇的逼水好多| 九九爱精品视频在线观看| 免费大片18禁| 亚洲美女视频黄频| freevideosex欧美| 九九在线视频观看精品| 三级经典国产精品| 成年人免费黄色播放视频 | av免费观看日本| 日韩精品有码人妻一区| 欧美日韩视频精品一区| 国产一区二区在线观看日韩| 国产极品粉嫩免费观看在线 | 日韩欧美 国产精品| 纵有疾风起免费观看全集完整版| 久久热精品热| 日韩精品有码人妻一区| 日韩伦理黄色片| 亚洲精品乱码久久久久久按摩| 日韩一本色道免费dvd| 免费观看av网站的网址| 永久网站在线| 国产av精品麻豆| 一个人看视频在线观看www免费| a级一级毛片免费在线观看| 国产一区有黄有色的免费视频| 成年女人在线观看亚洲视频| 日韩一区二区视频免费看| 国产精品国产av在线观看| 亚洲在久久综合| .国产精品久久| 午夜福利在线观看免费完整高清在| 国产女主播在线喷水免费视频网站| 国产免费一级a男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| 国产淫语在线视频| 亚洲av中文av极速乱| 久久久国产一区二区| 婷婷色综合www| 久久精品国产a三级三级三级| 黄色配什么色好看| 少妇人妻 视频| 青青草视频在线视频观看| 久久精品国产亚洲网站| 熟女av电影| 男的添女的下面高潮视频| 国内少妇人妻偷人精品xxx网站| 日韩三级伦理在线观看| 国产亚洲欧美精品永久| 春色校园在线视频观看| 国产亚洲av片在线观看秒播厂| 久久女婷五月综合色啪小说| 日日撸夜夜添| 中文字幕精品免费在线观看视频 | 色视频www国产| 久热这里只有精品99| 亚洲精品aⅴ在线观看| av国产精品久久久久影院| 五月天丁香电影| 国产爽快片一区二区三区| 乱系列少妇在线播放| 男男h啪啪无遮挡| 国产一区二区在线观看日韩| 91久久精品电影网| 亚洲国产精品一区二区三区在线| 亚洲av国产av综合av卡| 日韩三级伦理在线观看| 曰老女人黄片| 国产淫语在线视频| 最新的欧美精品一区二区| 一级毛片电影观看| 一级黄片播放器| 欧美97在线视频| 国产免费又黄又爽又色| 国产成人精品一,二区| 亚洲成人一二三区av| 亚洲成人一二三区av| 亚洲美女搞黄在线观看| 女性生殖器流出的白浆| 草草在线视频免费看| 日本av手机在线免费观看| 亚洲欧美日韩另类电影网站| 久久久久国产网址| 精品一区二区免费观看| 久久久久精品性色| 国产毛片在线视频| 少妇人妻 视频| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 成年人免费黄色播放视频 | 狂野欧美激情性bbbbbb| 亚洲,欧美,日韩| 日本与韩国留学比较| 偷拍熟女少妇极品色| 51国产日韩欧美| 久热这里只有精品99| 亚洲av欧美aⅴ国产| 亚洲内射少妇av| 丝袜脚勾引网站| 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 嘟嘟电影网在线观看| 免费看光身美女| 欧美日韩精品成人综合77777| 久久综合国产亚洲精品| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 免费观看无遮挡的男女| 大又大粗又爽又黄少妇毛片口| 午夜福利网站1000一区二区三区| 青青草视频在线视频观看| 美女视频免费永久观看网站| 肉色欧美久久久久久久蜜桃| 爱豆传媒免费全集在线观看| 欧美丝袜亚洲另类| 国产av码专区亚洲av| 亚洲av福利一区| 国产成人aa在线观看| 在线播放无遮挡| av播播在线观看一区| 嫩草影院新地址| 免费黄网站久久成人精品| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 一本大道久久a久久精品| 成人综合一区亚洲| 国产熟女午夜一区二区三区 | 中文字幕免费在线视频6| 九九在线视频观看精品| 国产精品一区二区在线不卡| 女的被弄到高潮叫床怎么办| 日本黄大片高清| 一级a做视频免费观看| 国产精品一区www在线观看| 国产精品一二三区在线看| 免费在线观看成人毛片| 久久精品国产a三级三级三级| 日韩电影二区| 亚洲人成网站在线观看播放| 内射极品少妇av片p| av福利片在线观看| 亚洲激情五月婷婷啪啪| 免费观看a级毛片全部| 看非洲黑人一级黄片| 久久精品国产a三级三级三级| 高清视频免费观看一区二区| 国产一区二区三区av在线| 黄色怎么调成土黄色| 青春草国产在线视频| 国产探花极品一区二区| 久久精品国产自在天天线| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品成人久久小说| 日本欧美视频一区| 免费观看av网站的网址| 日韩不卡一区二区三区视频在线| 久久影院123| 黄色毛片三级朝国网站 | 国产探花极品一区二区| 91午夜精品亚洲一区二区三区| 国产伦精品一区二区三区视频9| 亚洲精品国产色婷婷电影| 啦啦啦中文免费视频观看日本| av不卡在线播放| 亚洲成人一二三区av| 日韩在线高清观看一区二区三区| 我要看日韩黄色一级片| xxx大片免费视频| 在线观看美女被高潮喷水网站| 热99国产精品久久久久久7| 一级黄片播放器| 国产深夜福利视频在线观看| 中国美白少妇内射xxxbb| 99视频精品全部免费 在线| 深夜a级毛片| 久久久国产一区二区| 91精品国产国语对白视频| 欧美日韩亚洲高清精品| 国产精品99久久久久久久久| 亚洲精品,欧美精品| 多毛熟女@视频| 最近手机中文字幕大全| 日日啪夜夜爽| 精品熟女少妇av免费看| 精品人妻熟女毛片av久久网站| 黄色视频在线播放观看不卡| 午夜精品国产一区二区电影| 亚洲精品一区蜜桃| 97超碰精品成人国产| 男女免费视频国产| 亚洲av.av天堂| 久久久久久久精品精品| 在线 av 中文字幕| 亚洲美女搞黄在线观看| 久久狼人影院| 五月开心婷婷网| 男的添女的下面高潮视频| 极品人妻少妇av视频| 激情五月婷婷亚洲| 久久青草综合色| 啦啦啦啦在线视频资源| 国产在线男女| 精品熟女少妇av免费看| 18禁裸乳无遮挡动漫免费视频| 久久韩国三级中文字幕| 免费看av在线观看网站| 精品一品国产午夜福利视频| 国产熟女欧美一区二区| 久久久久精品久久久久真实原创| 一边亲一边摸免费视频| 91久久精品国产一区二区成人| 欧美3d第一页| 18+在线观看网站| 日本黄大片高清| 在线观看三级黄色| 亚洲欧美清纯卡通| 精品视频人人做人人爽| 老司机影院毛片| 777米奇影视久久| 欧美丝袜亚洲另类| 久久免费观看电影| 欧美 亚洲 国产 日韩一| 麻豆成人av视频| 精品午夜福利在线看| 国产精品久久久久久精品古装| 欧美三级亚洲精品| 午夜福利视频精品| 一个人看视频在线观看www免费| 在线观看三级黄色| 精品人妻一区二区三区麻豆| 日韩中字成人| a级毛色黄片| 国产精品久久久久久久电影| 国产免费视频播放在线视频| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 日日啪夜夜爽| 日韩精品有码人妻一区| 国产精品伦人一区二区| 婷婷色麻豆天堂久久| 国产在视频线精品| 国产白丝娇喘喷水9色精品| 91久久精品国产一区二区三区| 亚洲av在线观看美女高潮| 亚洲av中文av极速乱| 日韩免费高清中文字幕av| 内射极品少妇av片p| 免费久久久久久久精品成人欧美视频 | 国产成人freesex在线| 欧美精品国产亚洲| 久久久国产欧美日韩av| 精品国产乱码久久久久久小说| 国产美女午夜福利| 精品国产一区二区三区久久久樱花| 成人漫画全彩无遮挡| 国产免费又黄又爽又色| 国产成人免费观看mmmm| 亚洲人成网站在线播| 这个男人来自地球电影免费观看 | 一级爰片在线观看| 精品一区二区免费观看| 99热这里只有精品一区| 观看av在线不卡| 国产成人aa在线观看| 国产熟女欧美一区二区| 18禁在线无遮挡免费观看视频| 亚洲欧美成人精品一区二区| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久久丰满| 午夜av观看不卡| 国产 精品1| 黑丝袜美女国产一区| 久久综合国产亚洲精品| 熟女av电影| 大码成人一级视频| 观看免费一级毛片| 精品久久久久久电影网| 亚洲精品国产av蜜桃| 久久久午夜欧美精品| 黄色毛片三级朝国网站 | 黄片无遮挡物在线观看| 久久ye,这里只有精品| 久久人人爽人人片av| 国产在线免费精品| 日本与韩国留学比较| av一本久久久久| 国产精品一区www在线观看| 欧美变态另类bdsm刘玥| 自拍偷自拍亚洲精品老妇| 日日爽夜夜爽网站| 精品少妇黑人巨大在线播放| 日韩强制内射视频| 久久毛片免费看一区二区三区| 黑人高潮一二区| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区 | 三级国产精品欧美在线观看| 免费观看的影片在线观看| 久久久久久久久久久免费av| 国产成人免费无遮挡视频| 在线播放无遮挡| 精品久久久久久久久av| 成年女人在线观看亚洲视频| 国产成人精品婷婷| 免费少妇av软件| 99热国产这里只有精品6| 青春草国产在线视频| 国内少妇人妻偷人精品xxx网站| 春色校园在线视频观看| av.在线天堂| 精品少妇久久久久久888优播| 丁香六月天网| 女性生殖器流出的白浆| 一本大道久久a久久精品| 日韩在线高清观看一区二区三区| 亚洲精品日韩av片在线观看| videossex国产| 亚洲精品一二三| 亚洲精品乱久久久久久| 成人二区视频| 国产免费一级a男人的天堂| 欧美性感艳星| 黄色毛片三级朝国网站 | 日本黄大片高清| 777米奇影视久久| 色视频www国产| 成年美女黄网站色视频大全免费 | 久久影院123| 国产无遮挡羞羞视频在线观看| 国产欧美亚洲国产| 国产精品一区二区在线不卡| 国产高清三级在线| 亚洲精品自拍成人| 天美传媒精品一区二区| 亚洲国产精品国产精品| 街头女战士在线观看网站| 乱系列少妇在线播放| 性高湖久久久久久久久免费观看| 黑人猛操日本美女一级片| 国产淫片久久久久久久久| 少妇高潮的动态图| av又黄又爽大尺度在线免费看| 男人狂女人下面高潮的视频| 天堂中文最新版在线下载| 国产视频首页在线观看| 纯流量卡能插随身wifi吗| 精品酒店卫生间| 亚洲一级一片aⅴ在线观看| 成人二区视频| 肉色欧美久久久久久久蜜桃| 少妇猛男粗大的猛烈进出视频| 午夜免费男女啪啪视频观看| 亚洲av男天堂| 中文字幕制服av| 一级,二级,三级黄色视频| 卡戴珊不雅视频在线播放| 赤兔流量卡办理| 在线观看三级黄色| 亚洲国产成人一精品久久久| 亚洲电影在线观看av| 久久久久人妻精品一区果冻| 一级毛片aaaaaa免费看小| 亚洲国产精品专区欧美| 欧美 日韩 精品 国产| 22中文网久久字幕| 丰满迷人的少妇在线观看| 亚洲精品国产av蜜桃| 寂寞人妻少妇视频99o| 视频中文字幕在线观看| 国产精品麻豆人妻色哟哟久久| 国产精品成人在线| 老司机亚洲免费影院| 亚洲中文av在线| 91久久精品国产一区二区三区| 久热这里只有精品99| 国产精品欧美亚洲77777| 成人午夜精彩视频在线观看| 99九九线精品视频在线观看视频| 日本黄色日本黄色录像| av卡一久久| 久久久久久久久久人人人人人人| 国产亚洲精品久久久com| 你懂的网址亚洲精品在线观看| 国语对白做爰xxxⅹ性视频网站| 国产熟女欧美一区二区| 男人爽女人下面视频在线观看| 人妻夜夜爽99麻豆av| 丝袜脚勾引网站| a级毛片在线看网站| 在线免费观看不下载黄p国产| a级毛片在线看网站| 久久人人爽人人片av| 在线观看免费日韩欧美大片 | 国产精品不卡视频一区二区| 观看美女的网站| 国产在线免费精品| 亚洲综合精品二区| 男的添女的下面高潮视频| 少妇 在线观看| 亚洲真实伦在线观看| 另类精品久久| 国产伦精品一区二区三区视频9| 精品人妻一区二区三区麻豆| 久久久国产一区二区| 国产av一区二区精品久久| 噜噜噜噜噜久久久久久91| 亚洲性久久影院| 久久99热这里只频精品6学生| 国产成人freesex在线| 在线观看免费高清a一片| 99精国产麻豆久久婷婷| 日韩亚洲欧美综合| 男女边摸边吃奶| 日韩一区二区三区影片| 国产一区二区在线观看日韩| 精品亚洲成a人片在线观看| 成人国产av品久久久| 黄色配什么色好看| 精品久久久精品久久久| 亚洲av欧美aⅴ国产| 少妇熟女欧美另类| 精华霜和精华液先用哪个| 久久精品国产亚洲av天美| 一本—道久久a久久精品蜜桃钙片| 新久久久久国产一级毛片| av又黄又爽大尺度在线免费看| av专区在线播放| 丰满迷人的少妇在线观看| 国产欧美亚洲国产| 夜夜骑夜夜射夜夜干| 伊人亚洲综合成人网| 国产视频内射| 一区二区av电影网| 成人综合一区亚洲| 亚洲av电影在线观看一区二区三区| 精品人妻熟女毛片av久久网站| 毛片一级片免费看久久久久| 久久99一区二区三区| 亚洲精品第二区| 街头女战士在线观看网站| 成年av动漫网址| 免费不卡的大黄色大毛片视频在线观看| 一级毛片aaaaaa免费看小| 在线观看国产h片| 黄色一级大片看看| 最近手机中文字幕大全| 亚洲精品自拍成人| 中文字幕制服av| 我的女老师完整版在线观看| 成年女人在线观看亚洲视频| 欧美日韩视频高清一区二区三区二| 性高湖久久久久久久久免费观看| 在现免费观看毛片| 国产精品福利在线免费观看| 内射极品少妇av片p| 欧美精品亚洲一区二区| 日韩熟女老妇一区二区性免费视频| 久久久精品94久久精品| 日韩电影二区| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产a三级三级三级| 亚洲av国产av综合av卡| 国产精品熟女久久久久浪| 久久国产精品大桥未久av | 久久精品国产a三级三级三级| 精品久久久久久久久亚洲| 国产视频内射| 汤姆久久久久久久影院中文字幕| 乱系列少妇在线播放| 高清黄色对白视频在线免费看 | 插逼视频在线观看| 好男人视频免费观看在线| 国产白丝娇喘喷水9色精品| 国产av国产精品国产| 精品少妇黑人巨大在线播放| 欧美日韩视频高清一区二区三区二| 亚洲精品一二三| 亚洲精品国产色婷婷电影| h视频一区二区三区| 国产亚洲欧美精品永久| 亚洲三级黄色毛片| 亚洲一级一片aⅴ在线观看| 日韩欧美一区视频在线观看 | av.在线天堂| 久久精品国产亚洲网站| 精品久久久久久电影网| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩卡通动漫| 男女免费视频国产| 国产精品伦人一区二区| 久久久久视频综合| 少妇的逼水好多| 免费观看无遮挡的男女| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| 爱豆传媒免费全集在线观看| 久久亚洲国产成人精品v| 精品卡一卡二卡四卡免费| 全区人妻精品视频| 国产精品国产三级国产专区5o| 亚洲欧美成人综合另类久久久| 亚洲三级黄色毛片| 啦啦啦中文免费视频观看日本| 免费观看av网站的网址| 热re99久久精品国产66热6| 日日撸夜夜添| 在线播放无遮挡| 久久午夜福利片| 男女免费视频国产| 国产精品久久久久久久久免| 国产男人的电影天堂91| 精品人妻熟女毛片av久久网站| 高清欧美精品videossex| 亚洲精品久久久久久婷婷小说| 老熟女久久久| 日韩强制内射视频| 99精国产麻豆久久婷婷| 人妻系列 视频| 色婷婷av一区二区三区视频| 日韩精品免费视频一区二区三区 | 亚洲精品aⅴ在线观看| 欧美成人精品欧美一级黄| 在线观看三级黄色| 久久鲁丝午夜福利片| 国产精品嫩草影院av在线观看| 少妇人妻一区二区三区视频| 国产免费又黄又爽又色| 亚洲色图综合在线观看| 久久99一区二区三区| 如何舔出高潮| 国产男女内射视频| 女性生殖器流出的白浆| 中文乱码字字幕精品一区二区三区| 国产 一区精品| 中文在线观看免费www的网站| 成人18禁高潮啪啪吃奶动态图 | 婷婷色麻豆天堂久久| 高清在线视频一区二区三区| 青春草亚洲视频在线观看| 九九在线视频观看精品| 啦啦啦啦在线视频资源| 爱豆传媒免费全集在线观看| 一个人看视频在线观看www免费| 亚洲av免费高清在线观看| 日本色播在线视频| 精品视频人人做人人爽| 美女福利国产在线| 国产成人精品无人区| 在线亚洲精品国产二区图片欧美 | 国产精品99久久99久久久不卡 | 久久午夜福利片| 日产精品乱码卡一卡2卡三| 极品少妇高潮喷水抽搐| 国产综合精华液| 亚洲国产精品一区二区三区在线| 国产精品伦人一区二区| 高清av免费在线| 成人免费观看视频高清| 久久人人爽人人片av| 丰满人妻一区二区三区视频av| 内射极品少妇av片p| 亚洲精品日本国产第一区| 26uuu在线亚洲综合色| 国产av一区二区精品久久| av国产久精品久网站免费入址| 亚洲真实伦在线观看| 一级爰片在线观看| 中文字幕av电影在线播放| 久久热精品热| 国产日韩欧美亚洲二区| 人人妻人人添人人爽欧美一区卜| 成人毛片60女人毛片免费| 午夜精品国产一区二区电影| 男女国产视频网站| 黑丝袜美女国产一区| 热re99久久精品国产66热6| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久精品性色| 国产一区二区在线观看av| 午夜福利,免费看| 中文乱码字字幕精品一区二区三区| 最黄视频免费看| 人妻 亚洲 视频| 人妻人人澡人人爽人人| 亚洲天堂av无毛| 日韩欧美一区视频在线观看 | 日韩人妻高清精品专区| av黄色大香蕉| 久久99热这里只频精品6学生| 22中文网久久字幕| 乱码一卡2卡4卡精品| 亚洲精品久久午夜乱码| 伊人亚洲综合成人网| av线在线观看网站| 久久影院123| 天堂中文最新版在线下载| 亚洲人成网站在线播| 欧美国产精品一级二级三级 | 午夜免费鲁丝| 黄色日韩在线| 一本一本综合久久| 日韩一本色道免费dvd| 久久精品国产亚洲网站| 亚洲精品日韩在线中文字幕| 国产69精品久久久久777片| 国产极品粉嫩免费观看在线 | 乱码一卡2卡4卡精品| av在线播放精品| 2022亚洲国产成人精品| 成人综合一区亚洲| 王馨瑶露胸无遮挡在线观看|