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

    覆膜方式對小麥-玉米輪作農(nóng)田生態(tài)系統(tǒng)凈碳匯的影響

    2017-06-05 15:00:30劉晶晶張阿鳳鄒小陽陳海心
    關(guān)鍵詞:足跡覆膜通量

    馮 浩 劉晶晶 張阿鳳 鄒小陽 陳海心

    (1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院, 陜西楊凌 712100; 2.中國科學(xué)院水利部水土保持研究所, 陜西楊凌 712100;3.西北農(nóng)林科技大學(xué)中國旱區(qū)節(jié)水農(nóng)業(yè)研究院, 陜西楊凌 712100; 4.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院, 陜西楊凌 712100;5.西北農(nóng)林科技大學(xué)水土保持研究所, 陜西楊凌 712100)

    覆膜方式對小麥-玉米輪作農(nóng)田生態(tài)系統(tǒng)凈碳匯的影響

    馮 浩1,2劉晶晶1,3張阿鳳3,4鄒小陽2陳海心3,5

    (1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院, 陜西楊凌 712100; 2.中國科學(xué)院水利部水土保持研究所, 陜西楊凌 712100;3.西北農(nóng)林科技大學(xué)中國旱區(qū)節(jié)水農(nóng)業(yè)研究院, 陜西楊凌 712100; 4.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院, 陜西楊凌 712100;5.西北農(nóng)林科技大學(xué)水土保持研究所, 陜西楊凌 712100)

    采用靜態(tài)暗箱-氣相色譜法對關(guān)中平原小麥-玉米輪作(2014—2015年)農(nóng)田溫室氣體(CO2、CH4和N2O)排放通量進(jìn)行監(jiān)測,并用凈增溫潛勢(NGWP)和碳足跡2個(gè)指標(biāo)評估不同覆膜方式對關(guān)中平原農(nóng)田溫室效應(yīng)、作物生產(chǎn)碳足跡的分布和構(gòu)成的影響。試驗(yàn)處理設(shè)置為不覆膜(CK)、半膜覆蓋(BM)、壟作覆膜(LM)、全膜覆蓋(QM)。結(jié)果表明,與CK處理相比,BM、LM、QM處理下作物年際總產(chǎn)量分別增加了9.0%、16.5%、26.6%;CO2的年際排放總量分別增加了33.9%、9.3%、31.6%,N2O的年際排放總量分別增加了22.9%、14.3%、47.1%,但CH4的年際吸收總量無顯著性差異;NGWP分別增加了9.0%、16.7%、26.0%; LM、QM處理碳足跡較CK處理減少了33.2%、21.9%,而BM處理碳足跡與CK處理無顯著差異;BM處理單位產(chǎn)量碳足跡較CK處理增加了16.3%,而LM處理較CK處理減少了13.1%,QM處理與CK處理無顯著性差異。綜合考慮不同覆膜方式的經(jīng)濟(jì)效應(yīng)和環(huán)境效應(yīng),壟作覆膜更有利于關(guān)中平原小麥-玉米輪作系統(tǒng)的固碳減排。

    小麥-玉米; 覆膜; 溫室氣體; 凈增溫潛勢; 碳足跡

    引言

    溫室氣體(CO2、CH4和N2O)大量排放是全球氣候變暖及其引起的一系列環(huán)境問題的主要原因。2013年聯(lián)合國環(huán)境規(guī)劃署發(fā)布的《碳排放評估報(bào)告》指出,農(nóng)業(yè)源溫室氣體排放量占全球溫室氣體排放總量的11%,已超過2020年的排放目標(biāo)[1-2]。農(nóng)業(yè)生態(tài)系統(tǒng)不僅是溫室氣體CO2、CH4和N2O的主要排放源,而且還是重要的固碳減排體系,具有很強(qiáng)的固碳潛力[3]。地膜覆蓋是旱區(qū)農(nóng)業(yè)獲得高產(chǎn)穩(wěn)產(chǎn)的重要手段,中國已成為覆膜種植面積最大的國家[4]。覆膜具有增溫保墑、改善土壤性狀和提高水分利用率等諸多優(yōu)點(diǎn)[5-6],不同的覆膜方式對土壤水熱條件的影響存在差異。馬忠明等[7]研究了不同覆膜方式對旱地土壤水分和溫度的影響,發(fā)現(xiàn)全膜覆蓋和半膜壟作覆蓋的0~40 cm土壤含水量較半膜覆蓋分別提高了10.61%和6.75%,全膜覆蓋和半膜壟作覆蓋的0~15 cm土溫較半膜覆蓋分別增加了1.96℃和1.18℃。土壤水分和土壤溫度是溫室氣體排放與碳循環(huán)的關(guān)鍵驅(qū)動(dòng)因子,在一定變化范圍內(nèi)與溫室氣體排放通量具有顯著相關(guān)性[8-9]。陶麗佳等[10]通過馬鈴薯盆栽試驗(yàn)發(fā)現(xiàn),覆膜促使土壤CO2排放量升高10.4%~94.5%,CH4排放量降低5.1%~47.4%。NAN等[11]研究發(fā)現(xiàn),覆膜處理的CO2和N2O濃度更高,而CH4濃度更低;覆膜明顯促進(jìn)土壤CO2和N2O排放潛勢的提高,而阻礙CH4的排放。NISHIMURA等[12]研究了覆膜對菜地N2O氣體排放的影響,發(fā)現(xiàn)覆膜促進(jìn)了N2O的排放。目前,關(guān)于不同覆膜方式對農(nóng)田溫室效應(yīng)的影響研究較少且多集中在對單一氣體的研究方面。

    碳足跡作為新的研究方法得到學(xué)術(shù)界認(rèn)可,并成為氣候變化領(lǐng)域的研究熱點(diǎn)[13]。農(nóng)業(yè)碳足跡能夠系統(tǒng)地評價(jià)耕作、施肥、灌溉和收獲等農(nóng)業(yè)生產(chǎn)活動(dòng)過程中,由人為因素引起直接和間接的碳排放總量,定量測算農(nóng)業(yè)生產(chǎn)活動(dòng)對溫室效應(yīng)的影響[14]。然而,目前有關(guān)中國農(nóng)田管理對大氣溫室效應(yīng)和碳足跡貢獻(xiàn)影響的文獻(xiàn)中,大多為秸稈還田、免耕、有機(jī)肥施用等管理措施的相關(guān)研究結(jié)果,有關(guān)不同覆膜方式影響效應(yīng)的報(bào)道較少。因此,本文以關(guān)中平原小麥-玉米輪作農(nóng)田土壤為研究對象,通過設(shè)置不同覆膜方式處理,研究其對冬小麥-夏玉米生長季內(nèi)土壤溫室氣體排放規(guī)律的影響,用凈增溫潛勢(NGWP)和碳足跡2個(gè)指標(biāo)評估不同覆膜方式下的固碳減排效應(yīng),優(yōu)化田間管理措施,探尋關(guān)中平原最佳覆膜種植方式,以期為保障小麥、玉米高產(chǎn)、優(yōu)質(zhì)及農(nóng)業(yè)固碳減排提供依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)區(qū)概況

    田間試驗(yàn)位于陜西關(guān)中平原中部楊陵區(qū),西北農(nóng)林科技大學(xué)節(jié)水灌溉試驗(yàn)站(34°20′N、108°24′E,海拔高度521 m),屬半干旱半濕潤氣候區(qū)。2014年10月份—2015年10月份總降水量為540.8 mm,年平均溫度為14.3℃,試驗(yàn)站內(nèi)土壤為塿土,0~1 m土層的平均田間持水率為23%~25%,凋萎系數(shù)為8.5%(均為質(zhì)量含水率),播前0~20 cm土壤的有機(jī)碳含量為8.23 g/kg,全氮含量為0.95 mg/kg,土壤硝態(tài)氮含量為5.41 mg/kg,土壤銨態(tài)氮含量為1.35 mg/kg,有效磷(P2O5)含量為20.91 mg/kg,速效鉀(K2O)含量為287 mg/kg,pH值(水土質(zhì)量比1∶1)為8.32,土壤容重為1.43 g/cm3。

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

    試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì),設(shè)不覆膜(CK)、半膜覆蓋(BM)、壟作覆膜(LM)、全膜覆蓋(QM)4個(gè)處理,試驗(yàn)區(qū)布置如圖1所示。各處理重復(fù)3次,共12個(gè)小區(qū),每個(gè)小區(qū)面積為10 m2(2 m×5 m),使用塑料薄膜(聚乙烯,膜厚0.008 mm,每平方米質(zhì)量為6 g)覆蓋,實(shí)際使用半膜寬為0.4 m,壟膜寬為0.5 m(壟高15 cm),全膜寬為2.4 m。試驗(yàn)區(qū)周圍設(shè)0.5 m寬的作物保護(hù)帶。各處理施肥和灌水制度一致,以尿素和磷酸二氫鈣作基肥,其他管理措施與當(dāng)?shù)匾恢?。小?小偃22)于2014年10月15日播種(播種密度為187.5 kg/hm2),2015年6月6日收獲。玉米(秦龍14)于6月11日播種,10月8日收獲。各處理灌水和施肥時(shí)間及其用量如表1所示。

    圖1 試驗(yàn)布置示意圖Fig.1 Sketches of experiment arrangement1.冬小麥 2.滴灌帶 3.采氣箱 4.土壤 5.夏玉米

    日期灌水量/mm施肥量/(kg·hm-2)NP2O52014-10-1501201002014-11-1930002015-03-1530002015-04-093002015-06-10225902015-06-1230002015-08-023000

    1.3 測定項(xiàng)目和方法

    1.3.1 溫室氣體采集與測定

    溫室氣體采用靜態(tài)暗箱-氣相色譜法進(jìn)行采集和測定。采樣箱由箱體和底座組成,箱體用PVC材料制成,覆有泡沫和反光材料,尺寸為50 cm×50 cm×50 cm,箱內(nèi)安裝攪勻氣體的風(fēng)扇和抽氣接口,接口連接三通閥。底座用不銹鋼制成,尺寸為50 cm×50 cm×5 cm,底座內(nèi)無作物,測定裸地土壤溫室氣體的排放。底座上面附有凹槽,采樣之前槽內(nèi)注水密封以隔絕箱外環(huán)境。氣體采集時(shí)間為09:00—11:00之間??巯浜罅⒓从?0 mL注射器采集樣品,然后分別在扣箱后10、20、30 min時(shí)采樣一次。每10 d采集一次,若施肥、灌溉或降水,則增加采樣頻率,一共增加采樣3次,分別為2015年5月23日、7月16日、8月5日。采樣結(jié)束后立即帶回實(shí)驗(yàn)室用氣相色譜儀(Agilent 7890 A型,美國安捷倫科技有限公司)進(jìn)行分析,CH4和CO2使用氫火焰離子檢測器(FID)測定,載氣為氮?dú)?,柱溫?0℃,檢測器溫度為200℃,流速為40 mL/min;燃?xì)馐菤錃?,其流?5 mL/min;助燃?xì)馐强諝?,流?50 mL/min。N2O使用電子捕獲檢測器(ECD),載氣為氬甲烷,其柱溫是65℃,檢測器溫度330℃,流速為30 mL/min。氣體排放速率由4個(gè)連續(xù)樣品濃度的斜率經(jīng)過線性回歸分析得出,排放通量計(jì)算式[14]為

    (1)

    式中F——?dú)怏w的排放通量,mg/(m2·h)H——采樣箱高度,mM——?dú)怏w的摩爾質(zhì)量,g/molP——采樣點(diǎn)氣壓,PaR——普適氣體常數(shù),取8.314 Pa·m3/(mol·K)

    T——采樣時(shí)箱內(nèi)平均氣溫,℃

    用內(nèi)插法估算小麥、玉米生育期CO2、N2O和CH4的排放總量。

    1.3.2 產(chǎn)量、溫度及含水率測定

    小麥?zhǔn)斋@后,在每個(gè)小區(qū)選取1 m2的樣點(diǎn)進(jìn)行單獨(dú)收割,風(fēng)干脫粒后稱取籽粒質(zhì)量并計(jì)算產(chǎn)量,最終結(jié)果以每公頃小麥產(chǎn)量(t/hm2)表示。試驗(yàn)區(qū)自動(dòng)氣象站測定每日最高、最低氣溫和降水量(圖2),用地溫計(jì)測定5 cm深處土壤溫度、氣密閉箱內(nèi)空氣溫度,測定時(shí)間為密閉箱置于底座20 min后。采用TDR土壤水分速測儀測定土壤體積含水率,每個(gè)小區(qū)隨機(jī)測3個(gè)點(diǎn)取平均值,用土壤孔隙度含水率(Water filled pore space,WFPS)表示,計(jì)算公式為

    WFPS=θv/(2.65-γ)×100%

    (2)

    式中WFPS——土壤孔隙度含水率,%θv——土壤體積含水率,%γ——土壤容重,g/cm3

    圖2 小麥-玉米季的降水量、最低氣溫和最高氣溫Fig.2 Daily precipitation and maximum and minimum air temperatures in wheat-maize growing season

    1.4 生態(tài)系統(tǒng)凈交換量和凈增溫潛勢計(jì)算

    作物凈初級生產(chǎn)力(NPP)作為表征植被活力的關(guān)鍵變量,不僅是有機(jī)體能量和物質(zhì)流動(dòng)的基礎(chǔ),關(guān)系到生態(tài)系統(tǒng)對碳固定能力的強(qiáng)弱。生態(tài)系統(tǒng)凈交換是指生態(tài)系統(tǒng)光合同化作用與呼吸作用之間的平衡,表征了陸地生態(tài)系統(tǒng)吸收大氣CO2的能力。生態(tài)系統(tǒng)凈交換量(NEE)的計(jì)算式[15-16]為

    NEE=RH-NPP

    (3)

    其中NPP=0.446Wmax-0.000 67

    式中NPP——凈初級生產(chǎn)力(以碳等量(Ce)計(jì)),kg/hm2

    Wmax——作物收獲后地上部分和地下部分生物量的總和,kg

    RH——土壤CO2排放量(以Ce計(jì)),kg/hm2

    凈增溫潛勢(NGWP)是參考土壤和生物量排放法(Soil and crop-based approach)[2],由于本研究將物資和農(nóng)事活動(dòng)投入計(jì)入碳足跡進(jìn)行評價(jià),所以不再計(jì)入到增溫潛勢,則NGWP計(jì)算式為

    NGWP=GWPNPP+GWPΔSOC-GWPGHGs

    (4)

    式中NGWP——凈增溫潛勢,kg/hm2,以CO2當(dāng)量計(jì)算

    GWPNPP——凈初級生產(chǎn)力的增溫潛勢,kg/hm2,以CO2當(dāng)量計(jì)算

    GWPΔSOC——試驗(yàn)前后有機(jī)碳含量變化的增溫潛勢,由于本試驗(yàn)只有一年可忽略不計(jì)

    GWPGHGs——土壤CO2(秸稈部分全部以CO2形式排放到大氣中[17])、N2O和CH4排放總量的增溫潛勢,kg/hm2,以CO2當(dāng)量計(jì)算

    當(dāng)NGWP>0時(shí),表現(xiàn)為溫室氣體的匯,反之為源。

    1.5 碳足跡評價(jià)方法與計(jì)算

    碳足跡的邊界為開始種植小麥到收獲玉米結(jié)束(2014年10月份—2015年10月份),評價(jià)農(nóng)資生產(chǎn)過程中(肥料和農(nóng)藥)的碳排放、農(nóng)業(yè)活動(dòng)中的能源消耗(施肥、噴灑農(nóng)藥、耕作、收獲、起壟和覆膜)、地膜生產(chǎn)過程中的碳排放、灌溉過程中的能源消耗(電力)、農(nóng)田生態(tài)系統(tǒng)凈交換量(NEE)表征生態(tài)系統(tǒng)的碳吸收及小麥和玉米生育期內(nèi)CH4和N2O排放總量進(jìn)行分析,評價(jià)不同覆膜方式對農(nóng)田生態(tài)系統(tǒng)的固碳能力。

    小麥-玉米輪作生長過程總碳足跡CF的計(jì)算式[18]為

    CF=∑AiEFj

    (5)

    式中Ai——各農(nóng)業(yè)投入的總量

    EFj——相應(yīng)的排放參數(shù),正值表示碳消耗,負(fù)值表示碳吸收

    作物生產(chǎn)單位產(chǎn)量碳足跡CFy(以CO2-Ce計(jì))[19]計(jì)算式為

    (6)

    式中CFy——作物生產(chǎn)單位產(chǎn)量碳足跡,kg/tTY——小麥和玉米的總產(chǎn)量,t/hm2

    1.6 數(shù)據(jù)處理與分析

    試驗(yàn)數(shù)據(jù)均取3次重復(fù)的平均值,采用Microsoft Excel 2013和JMP 10.0(SAS Institute,Cary,NC,USA,2011)統(tǒng)計(jì)軟件進(jìn)行處理和相關(guān)性分析,用Duncan法檢驗(yàn)差異顯著性(α=0.05)。應(yīng)用OriginPro 9.0軟件對數(shù)據(jù)進(jìn)行制圖。

    2 結(jié)果與分析

    2.1 土壤CO2排放通量的季節(jié)變化

    從圖3a和圖4a中可以看出,各處理小麥-玉米生育期CO2排放通量(以CO2-C計(jì))季節(jié)變化趨勢與土壤溫度變化規(guī)律相似,呈現(xiàn)3月份之前較低、4—8月份較高的趨勢。CK處理的CO2排放通量在冬小麥生長季內(nèi)基本低于覆膜處理,降水或灌溉使不同覆膜方式處理CO2排放通量差異較大。進(jìn)入8月份之后伴隨溫度降低和高頻率降水各處理CO2排放通量呈降低趨勢(圖3)。從表2可知小麥季各覆膜處理CO2排放總量與CK處理存在顯著性差異(p<0.05),玉米季QM處理與CK處理存在顯著性差異(p<0.05)。BM、LM、QM處理下,小麥-玉米季CO2排放總量較CK處理分別增加了33.9%、9.3%、31.6%。

    圖3 小麥-玉米生長期不同覆膜方式CO2、N2O和CH4季節(jié)性排放通量Fig.3 Seasonal emission fluxes of CO2, N2O and CH4 under different mulching patterns during wheat-maize growing season

    圖4 不同覆膜方式土壤溫度和土壤孔隙度含水率變化Fig.4 Dynamics of soil temperature and soil WFPS under different mulching patterns

    2.2 土壤N2O、CH4排放和吸收通量的季節(jié)變化

    由圖3b可知,各處理N2O排放通量(以N2O-N計(jì))在整個(gè)小麥-玉米生育期內(nèi)變化趨勢一致,在施肥或灌溉后共出現(xiàn)了4個(gè)小排放峰和1個(gè)大排放峰,由于小麥季溫度較低在施肥或灌溉后只出現(xiàn)小排放峰,且各處理間差異不明顯。玉米種植后灌溉、施肥及高溫使得各處理N2O排放通量出現(xiàn)極大值,各處理差異較明顯,最高的QM處理達(dá)到了72.3 μg/(m2·h)(6月17日)。從小麥季N2O排放總量來看(表2),BM、LM、CK處理間無顯著性差異(p>0.05),QM處理較CK處理顯著增加(p<0.05)。玉米季BM、QM處理N2O排放總量較CK處理均顯著增加(p<0.05)。總的來說,BM、LM、QM處理下小麥-玉米季N2O排放總量較CK處理分別增加了22.9%、14.3%、47.1%。

    從圖3c可知,小麥-玉米季土壤CH4既有吸收特征又有排放特征,基本上呈現(xiàn)土壤CH4的吸收現(xiàn)象。各處理CH4通量(以CH4-C計(jì))在-0.33~0.3 mg/(m2·h)之間變化。灌溉和降水后出現(xiàn)CH4排放現(xiàn)象,總結(jié)發(fā)現(xiàn)CH4排放現(xiàn)象基本出現(xiàn)在WFPS大于60%時(shí)(圖4b)。由表2可知,小麥、玉米季各處理土壤均為CH4的匯,小麥季各處理無顯著性差異(p>0.05),而玉米季LM、QM處理較CK處理顯著降低了CH4的吸收。從年際CH4吸收總量來看,各處理無顯著性差異(p>0.05)。

    2.3 小麥、玉米產(chǎn)量和NGWP的變化

    由表2可知,覆膜增加了作物產(chǎn)量,小麥季LM處理產(chǎn)量較CK處理顯著增加,BM、QM、CK處理小麥產(chǎn)量無顯著性差異(p>0.05),玉米季QM處理產(chǎn)量較CK處理顯著增加,BM、LM、CK處理玉米產(chǎn)量無顯著性差異(p>0.05)(表2)。相比CK處理,BM、LM、QM處理的年際總產(chǎn)量分別增加了9.0%、16.5%、26.6%。本試驗(yàn)NGWP主要取決于籽粒產(chǎn)量,各處理的凈增溫潛勢均表現(xiàn)為溫室氣體的匯,且籽粒產(chǎn)量越多,固碳潛力越大,然而各處理間NGWP差異與產(chǎn)量差異一致。BM處理小麥-玉米季NGWP與CK處理無顯著性差異(p>0.05),LM、QM處理小麥-玉米季NGWP較CK處理顯著增加16.7%、26.0%。

    2.4 碳足跡的變化

    從小麥-玉米季碳足跡分析可知,不同處理均表現(xiàn)為農(nóng)田生態(tài)系統(tǒng)碳匯(表3)。LM、QM處理碳足跡較CK處理減少了33.2%、21.9%,而BM處理碳足跡與CK處理無顯著差異(p>0.05)。從圖5的碳足跡構(gòu)成可知,不同處理CO2凈交換量(NEE)占總碳匯的99.2%~99.8%。BM、LM、QM處理總碳消耗量較CK處理增加了41.7%、51.7%、76.9%。碳消耗量的排放來源主要是肥料和地膜,其次是灌溉和N2O排放,分別占各處理總碳消耗量的33.1%~58.6%和27.1%~40.7%、9.9%~17.5%和6.5%~8.6%。從圖6可知,BM處理單位產(chǎn)量碳足跡較CK處理增加了16.3%,而LM處理較CK處理減少了13.1%,QM處理與CK處理無顯著性差異(p>0.05)。

    表2 不同覆膜方式的產(chǎn)量、溫室氣體排放總量和凈增溫潛勢

    注:同列不同字母表示處理間差異顯著(p<0.05)。

    表3 不同覆膜方式的碳消耗

    3 討論

    3.1 不同覆膜方式對溫室氣體排放的影響

    覆膜可改善土壤水熱狀況,從而增強(qiáng)土壤微生物生長和活性[26]和其他影響土壤CO2通量的土壤性質(zhì),如土壤孔隙度[27]。土壤CO2排放通量趨勢與氣溫、土壤溫度變化趨勢基本一致,且與土壤溫度呈顯著正相關(guān)(表4),而土壤水分狀況對CO2產(chǎn)生與排放的影響機(jī)制主要體現(xiàn)在土壤通氣性和土壤溶解性有機(jī)質(zhì)含量的變化方面, 土壤過濕或過干都會(huì)抑制土壤微生物活動(dòng)。在小麥和玉米生長期,不同覆膜方式處理均可促進(jìn)CO2排放,因?yàn)楦材じ纳仆寥浪疅?,從而增?qiáng)土壤微生物活性,加速土壤有機(jī)質(zhì)的分解[28],且有利于植物根系和土壤微生物呼吸產(chǎn)生CO2,增加了土壤CO2濃度[29-30],與本研究的結(jié)論一致。但有研究認(rèn)為覆膜的物理阻隔大,從而降低了CO2排放通量[10,29,31],本研究得出覆膜增加了CO2排放通量,這可能是由于膜的滲透性導(dǎo)致CO2排放增加。因此,與前人研究CO2排放通量的差異可能原因有兩方面:一是覆膜土壤儲(chǔ)存的高濃度CO2可能通過膜的滲透[12]和水平擴(kuò)散而進(jìn)入大氣[32];二是未種植作物的氣體采集區(qū)周圍分布了少量作物根系,土壤CO2排放量的測定沒有為溫室氣體效應(yīng)提供一個(gè)可靠的估計(jì),而該效應(yīng)與CO2凈排放量(植株根系光合作用和呼吸作用及土壤呼吸作用的凈值)相關(guān)。

    圖5 不同覆膜方式的小麥-玉米季碳足跡構(gòu)成Fig.5 Composition of carbon footprint of wheat-maize growing season under different mulching patterns

    圖6 不同覆膜方式的單位小麥-玉米產(chǎn)量的碳足跡Fig.6 Carbon footprint of per wheat-maize yield under different mulching patterns

    不同管理措施下N2O排放量與氮肥施用、土壤水分、土壤溫度都呈正相關(guān)[31,33],這與本研究結(jié)果一致(表4)。覆膜措施下,土壤的水熱條件明顯得到協(xié)同改善。此時(shí),土壤微生物的活性和數(shù)量迅速增加,使得土壤硝化和反硝化過程加強(qiáng),土壤中N2O產(chǎn)生和排放量也隨之明顯增加。此外,由于地表的覆蓋,西北干旱半干旱地區(qū)土壤表層較為常見的干濕交替和凍融交替過程被有效地抑制,使得以此方式損耗的氮素減少,膜下的氮素得到一定程度地累積,在土壤水熱條件較好時(shí),土壤微生物過程加強(qiáng),氮素作為反硝化過程的底物,該累積過程反倒有可能增加土壤N2O的排放量[34]。NIAHIMURA等[12]研究發(fā)現(xiàn),透過地膜的N2O通量比從水平方向擴(kuò)散的更高,N2O透過地膜排出,其通透性隨周圍環(huán)境溫度的增加而增強(qiáng),這可能是本研究中QM處理在小麥-玉米季N2O排放量高于其他處理的原因。

    表4 不同覆膜方式與土壤溫度(5 cm)和土壤孔隙度含水量(10 cm)的相關(guān)系數(shù)

    注:*表示在p<0.05水平顯著相關(guān);** 表示在p<0.01水平顯著相關(guān)。

    土壤CH4的產(chǎn)生是產(chǎn)甲烷菌與甲烷氧化菌共同作用的結(jié)果,而旱地一般認(rèn)為是CH4的匯[35]。CH4在土壤中的排放和吸收受農(nóng)業(yè)管理措施和氣候季節(jié)性變化規(guī)律的影響。CUELLO等[28]研究結(jié)果表明地膜覆蓋減少了土壤CH4的吸收,因?yàn)榈啬じ采w使土壤中氧氣含量減少,微生物通常是在土壤厭氧條件下分解土壤中有機(jī)物。本試驗(yàn)各處理在冬小麥返青期之前WFPS較低,基本處于20%~45%之間。KHALIL等[36]試驗(yàn)發(fā)現(xiàn),WFPS對CH4氧化速率有顯著影響,WFPS為30%時(shí)土壤中CH4氧化速率最低。然而小麥季覆膜處理土壤對CH4的吸收與CK處理無顯著性差異。相反,本試驗(yàn)玉米季高頻率降水使各處理WFPS增加,土壤處于厭氧狀態(tài)促進(jìn)了甲烷菌分解有機(jī)物,丁維新等[37]研究表明甲烷氧化率與土壤濕度呈負(fù)相關(guān),這與本研究結(jié)果一致(表4)。因此,壟作的集水作用和覆膜的保水作用使得玉米季LM、QM處理顯著減少了土壤CH4的吸收。

    3.2 不同覆膜方式對NGWP的影響

    趙自超等[2]用土壤碳庫排放法、生物量量法、土壤和生物量排放法來評估硝化和脲酶抑制劑對華北冬小麥-夏玉米輪作固碳減排效果,其中土壤和生物量排放法綜合考慮了溫室氣體和農(nóng)田的作物產(chǎn)量。這里采用此方法評估不同覆膜方式之間對農(nóng)田生態(tài)系統(tǒng)的固碳潛力,主要取決于NPP,雖然各覆膜處理的N2O、CH4排放總量較CK處理有所增加,但其本身所占NGWP比重較小,這與劉巽浩等[38]的研究結(jié)果一致。本試驗(yàn)將NPP移出農(nóng)田部分的收獲物(作物籽粒)作為固碳單元,故整個(gè)系統(tǒng)移出農(nóng)田部分越高(產(chǎn)量越高),固碳潛力越大。岳維云等[39]研究不同覆膜方式下冬小麥增產(chǎn)效應(yīng)由大到小表現(xiàn)為:全膜覆土、膜側(cè)(起壟覆膜側(cè)播)、全膜穴播、露地,這與本試驗(yàn)冬小麥產(chǎn)量LM處理高于QM處理的結(jié)果一致。李尚中等[40]探求不同覆膜方式對旱地農(nóng)田玉米產(chǎn)量的影響時(shí)產(chǎn)量結(jié)果由大到小為:全平膜穴播、寬平膜、膜側(cè)(起壟覆膜側(cè)播)、露地,這與本文的試驗(yàn)結(jié)果一致。各處理之間的籽粒產(chǎn)量差異與NGWP差異基本一致。CK、BM、LM和QM處理NGWP分別高達(dá)19 083、20 799、22 269、24 043 kg/hm2,相比趙自超等[2]的研究結(jié)果要高,這是因?yàn)楣喔?、機(jī)械和肥料施用等農(nóng)事活動(dòng)所造成的CO2排放當(dāng)量沒有計(jì)入,而是把這一部分放入碳足跡來評價(jià),且高籽粒產(chǎn)量和N2O排放是LM、QM處理下NGWP高于其他2個(gè)處理的原因所在。

    3.3 不同覆膜方式對碳足跡的影響

    許多研究在評價(jià)農(nóng)田生態(tài)系統(tǒng)碳平衡時(shí),只考慮了農(nóng)業(yè)投入產(chǎn)生的碳足跡,而沒有考慮土壤呼吸CO2排放量、生態(tài)系統(tǒng)凈初級生產(chǎn)力以及農(nóng)田生態(tài)系統(tǒng)CH4和N2O的排放量[18-20,41]。本研究綜合考慮這些方面全面系統(tǒng)地分析不同灌溉量對農(nóng)田生態(tài)系統(tǒng)的碳平衡的影響。碳足跡是碳匯和碳源(總碳消耗量)的差值,本試驗(yàn)各處理NEE占碳匯的99.2%~99.8%。成功等[14]研究結(jié)果表明各處理NEE占碳匯的97.0%~99.7%。可見NEE是影響碳匯的最關(guān)鍵因子,而NEE主要受溫度、降水等氣候因子及農(nóng)田管理措施影響[42]。徐昔保等[43]研究結(jié)果表明太湖流域典型稻麥輪作農(nóng)田年NEE為-7 494~-7 853 kg/hm2。這與本研究NEE值(-5 332.5~-7 823 kg/hm2)結(jié)果相近。本試驗(yàn)中,相比CK處理,BM處理既增加了NPP也增加了土壤呼吸,使得生態(tài)系統(tǒng)碳足跡無顯著差異,LM、QM處理顯著增加了NPP,但土壤呼吸增加比例不大,從而減少了NEE??梢姴煌材し绞皆诓煌潭壬细纳屏送寥浪疅釛l件,活化土壤養(yǎng)分,促進(jìn)作物地上和地下部分的生長(即碳的固定)。從碳足跡的構(gòu)成可知(圖5),碳消耗的排放來源主要是肥料和地膜,主要原因是化肥本身在生產(chǎn)和運(yùn)輸?shù)倪^程中需要消耗大量的化石燃料,另一方面化肥施用量相對較多,造成化肥部分的碳排放量相對最大[41]。覆膜雖增加了農(nóng)田生態(tài)系統(tǒng)的碳匯,但在生產(chǎn)和使用過程中的碳消耗量占總碳消耗量的27.3%~40.9%。因此,高效利用地膜能有效減少碳消耗。LM處理顯著減少了單位產(chǎn)量碳足跡,是因?yàn)檩^CK處理產(chǎn)量增加幅度小于碳足跡減小幅度。

    4 結(jié)論

    (1)覆膜的增溫保墑作用促進(jìn)水分消耗由物理過程向生理過程轉(zhuǎn)化,由無效消耗向有效消耗轉(zhuǎn)化,從而有效增加作物籽粒產(chǎn)量。不同覆膜方式改善土壤理化性質(zhì)的程度不同,對作物的增產(chǎn)效應(yīng)也不同,本研究半膜覆蓋、壟作覆膜、全膜覆蓋條件下作物年總產(chǎn)量較不覆膜處理分別增加了9.0%、16.5%、26.6%。

    (2)覆膜處理小麥-玉米季CO2和N2O排放總量較不覆膜處理分別增加了9.3%~33.9%和14.3%~47.1%。各處理CH4吸收無顯著性差異(p>0.05)。相關(guān)性分析表明,土壤溫度和水分是影響溫室氣體排放的重要因素。

    (3)凈增溫潛勢(NGWP)和碳足跡2個(gè)指標(biāo)評價(jià)時(shí),結(jié)果存在差異。全膜覆蓋處理NGWP較不覆膜處理顯著增加了26.0%,壟作覆膜處理與全膜覆蓋處理無顯著性差異(p>0.05)。而半膜覆蓋處理NGWP與不覆膜處理無顯著性差異(p>0.05)。半膜覆蓋處理單位產(chǎn)量碳足跡較不覆膜處理增加了16.3%,而壟作覆膜處理較不覆膜處理減少了13.1%,QM處理與CK處理無顯著性差異(p>0.05)。因此,綜合2個(gè)指標(biāo)考慮不同覆膜方式的經(jīng)濟(jì)效應(yīng)和環(huán)境效應(yīng),壟作覆膜處理更有利于關(guān)中平原小麥-玉米輪作系統(tǒng)的固碳減排。

    1 UNEP. The global launch of the emissions gap report[R]. United Nations Environment Programme, Nairobi,2013.

    2 趙自超, 韓笑, 石岳峰, 等. 硝化和脲酶抑制劑對華北冬小麥-夏玉米輪作固碳減排效果評價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(6):254-262. ZHAO Zichao, HAN Xiao, SHI Yuefeng,et al. Effect of nitrification and urease inhibitor on carbon sequestration and greenhouse gas emissions in winter wheat and summer maize rotation system in North China[J]. Transactions of the CSAE, 2016,32(6):254-262.(in Chinese)

    3 IPCC. Change in atmospheric constituents and in radioactive forcing[R]. Climate Change. The Physical Basis. Cambridge University Press, Cambridge, United Kingdom, 2007.

    4 邢旭光, 趙文剛, 馬孝義, 等. 覆膜滴灌條件下棉花根層土壤鹽分時(shí)間穩(wěn)定性研究[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015, 46(7):146-153.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20150722&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2015.07.022. XING Xuguang, ZHAO Wen’gang, MA Xiaoyi, et al. Temporal stability of soil salinity in root zone of cotton under drip irrigation with plastic mulch[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 146-153. (in Chinese)

    5 蘇旺, 屈洋, 馮佰利, 等. 溝壟覆膜集水模式提高糜子光合作用和產(chǎn)量[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(13): 137-145. SU Wang, QU Yang, FENG Baili, et al. Photosynthesis characteristics and yield of broomcorn millet under film mulching on ridge-furrow for harvesting rainwater model in semi-arid region of Northern Shaanxi[J]. Transactions of the CSAE, 2014, 30(13): 137-145. (in Chinese)

    6 劉芳, 李天安, 樊小林. 華南地區(qū)覆膜旱種稻田甲烷排放及其與土壤水分和溫度的關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2013, 29(2):110-116. LIU Fang, LI Tian’an, FAN Xiaolin. Methane emission and its relationship with soil temperature and moisture during rice growth in film mulching upland rice field in South China[J]. Transactions of the CSAE, 2013, 29(2): 110-116. (in Chinese)

    7 馬忠明, 白玉龍, 薛亮, 等. 不同覆膜栽培方式對旱地土壤水熱效應(yīng)及西瓜產(chǎn)量的影響[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(3): 514-522. MA Zhongming, BAI Yulong, XUE Liang, et al. Effects of different plastic film mulching methods on soil water and temperature as well as watermelon yield in loess dry land[J]. Scientia Agricultura Sinica, 2015, 48(3): 514-522.(in Chinese)

    8 FLANAGAN L B, JOHNSON B G. Interacting effects of temperature, soil moisture and plant biomass production on ecosystem respiration in a northern temperate grassland[J]. Agricultural and Forest Meteorology, 2005, 130(3-4): 237-253.

    9 NORTON U, MOSIER A R, MORGAN J A, et al. Moisture pulses, trace gas emissions and soil C and N in cheat grass and native grass-dominated sagebrush-steppe in Wyoming, USA[J]. Soil Biology and Biochemistry, 2008, 40(6): 1421-1431.

    10 陶麗佳, 王鳳新, 顧小小. 膜下滴灌對土壤CO2與CH4濃度的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2012, 20(3): 330-336. TAO Lijia, WANG Fengxin, GU Xiaoxiao. Influence of drip irrigation under plastic film mulching on concentration of CO2and CH4in soil[J].Chinese Journal of Eco-Agriculturte, 2012, 20(3): 330-336.(in Chinese)

    11 NAN W G, YUE S C, HUANG H Z, et al. Effects of plastic film mulching on soil greenhouse gases (CO2, CH4and N2O) concentration within soil profiles in maize fields on the Loess Plateau, China[J]. Journal of Integrative Agriculture, 2016, 15(2): 451-464.

    12 NISHIMURA S, KOMADA M, TAKEBE M, et al. Contribution of nitrous oxide emission from soil covered with plastic mulch film in vegetable field[J]. Journal of Agricultural Meteorology, 2014, 70(2): 117-125.

    13 段華平, 張悅, 趙建波, 等. 中國農(nóng)田生態(tài)系統(tǒng)的碳足跡分析[J]. 水土保持學(xué)報(bào), 2011, 25(5): 203-208. DUAN Huaping, ZHANG Yue, ZHAO Jianbo, et al. Carbon footprint analysis of farmland ecosystem in China[J]. Journal of Soil and Water Conservation, 2011, 25(5): 203-208.(in Chinese)

    14 成功, 張阿鳳, 王旭東, 等. 運(yùn)用“碳足跡”的方法評估小麥秸稈及其生物質(zhì)炭添加對農(nóng)田生態(tài)系統(tǒng)凈碳匯的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2016, 35(3): 604-612. CHENG Gong, ZHANG Afeng, WANG Xudong, et al. Assessment of wheat straw and its biochar effects on carbon sink in agricultural ecosystems using “carbon footprint” method[J]. Journal of Agro-Environment Science, 2016, 35(3): 604-612.(in Chinese)

    15 OSAK M, SHIANNO T, TADANO T, et al. Carbon-nitrogen interaction model in field crop production[J].Plant and Soil, 1993,155-156(1):203-206.

    16 RAICH J W, TUFEKCIOGLU A. Vegetation and soil respiration:correlations and controls[J]. Biogeochemistry, 2000, 48(1):71-90.

    17 黃堅(jiān)雄,隋鵬,高旺盛,等. 華北平原玉米II大豆間作農(nóng)田溫室氣體排放及系統(tǒng)凈溫室效應(yīng)評價(jià)[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào),2015,20(4):66-74. HUANG Jianxiong, SUI Peng, GAO Wangsheng, et al. Effect of maize-soybean inter cropping on greenhouse gas emission and the assessment in North China Plain[J]. Journal of China Agricultural University, 2015, 20(4): 66-74. (in Chinese)

    18 LIU Q, LIU B, AMBUS P, et al. Carbon footprint of rice production under biochar amendment—a case study in a Chinese rice cropping system[J]. Global Chang Biology Bioenergy, 2016, 8(1): 148-159.

    19 王占彪, 王猛, 陳阜. 華北平原作物生產(chǎn)碳足跡分析[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(1): 83-92. WANG Zhanbiao, WANG Meng, CHEN Fu. Carbon footprint analysis of crop production in North China Plain[J]. Scientia Agricultura Sinica, 2015, 48(1): 83-92. (in Chinese)

    20 LAL R. Carbon emission from farm operations[J]. Environment International, 2004, 30(7): 981-990.

    21 HILLIER J, HAWES C, SQUIRE G, et al. The carbon footprints of food crop production[J]. International Journal of Agricultural Sustainability, 2009, 7(2): 107-118.

    22 YANG S H. Primary study on effect of C—O balance of afforestatal trees in cities[J]. Urban Environment and Urban Ecology, 1996, 9(1):37-39

    23 殷文, 史倩倩, 郭瑤, 等. 秸稈還田, 一膜兩年用及間作對農(nóng)田碳排放的短期效應(yīng)[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2016, 24(6): 716-724. YIN Wen, SHI Qianqian, GUO Yao, et al. Short-term response of farmland carbon emission to straw return, two-year plastic film mulching and intercropping[J]. Chinese Journal of Eco-Agriculture,2016, 24(6):716-724.(in Chinese)

    24 MOSIER A R, HALVORSON A D, RETLE C A, et al. Net global warming potential and greenhouse gas intensity in irrigated cropping systems in northeastern Colorado[J]. Journal of Environmental Quality, 2006, 35(4): 1584-1598.

    25 Intergovernmental Panel on Climate Change(IPCC). Climate change 2013: the physical science basis[R]. Cambridge:Cambridge University Press, 2014.

    26 SINGH H V. Soil carbon sequestration and rhizospheric microbial population in apricot orchards following plastic film mulching under cold arid region[J]. International Journal of Horticulture, 2013, 3: 35-41.

    27 RETH S, GOCKEDE M, FALGE E. CO2efflux from agricultural soils in Eastern Germany—comparison of a closed chamber system with eddy covariance measurements[J]. Theoretical and Applied Climatology, 2005, 80(2): 105-120.

    28 CUELLO J P, HWANG H Y, GUTIERREZ J, et al. Impact of plastic film mulching on increasing greenhouse gas emissions in temperate upland soil during maize cultivation[J]. Applied Soil Ecology, 2015, 91: 48-57.

    29 LI Z G, ZHANG R H, WANG X J, et al. Carbon dioxide fluxes and concentrations in a cotton field in northwestern China: effects of plastic mulching and drip irrigation[J]. Pedosphere, 2011, 21(2): 178-185.

    30 ZHAO Z M, ZHAO C Y, YILIHAMU Y, et al. Contribution of root respiration to total soil respiration in a cotton field of Northwest China[J]. Pedosphere, 2013, 23(2): 223-228.

    31 OKUDA H, NODA K, SAWAMOTO T, et al. Emission of N2O and CO2and uptake of CH4in soil from a satsuma mandarin orchard under mulching cultivation in central Japan[J]. Journal of the Japanese Society for Horticultural Science, 2007, 76(4): 279-287.

    32 俞永祥, 趙成義, 賈宏濤,等. 覆膜對綠洲棉田土壤CO2通量和CO2濃度的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(1): 155-160. YU Yongxiang, ZHAO Chengyi, JIA Hongtao, et al. Effects of plastic film mulching on soil CO2efflux and CO2concentration in an oasis cotton field[J]. Chinese Journal of Applied Ecology, 2015, 26(1): 155-160. (in Chinese)

    33 XU Y C, SHEN Q R, LI M L, et al. Effect of soil water status and mulching on N2O and CH4emission from lowland rice field in China[J]. Biology and Fertility of Soils, 2004, 39(3): 215-217.

    34 白紅英, 韓建剛, 張一平. 覆蓋種植措施對農(nóng)田土壤中N2O排放的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2003, 22(4): 394-396. BAI Hongying, HAN Jian’gang, ZHANG Yiping. Effects of mulching and planting on N2O discharge flux from soil[J]. Journal of Agro-Enviroment Science, 2003, 22(4): 394-396. (in Chinese)

    35 DALAL R C, ALLEN D E, LIVESLEY S J, et al. Magnitude and biophysical regulators of methane emission and consumption in the Australian agricultural, forest, and submerged landscapes: a review[J]. Plant and Soil, 2008, 309(1-2): 43-76.

    36 KHALIL M I, BAGGS E M. CH4oxidation and N2O emissions at varied soil water-filled pore spaces and headspace CH4concentrations[J]. Soil Biology & Biochemistry,2005,37(10):1785-1794.

    37 丁維新,蔡祖聰.土壤甲烷氧化菌及水分狀況對其活性的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2003, 11(1): 94-97. DING Weixin, CAI Zucong. Mechanism of methane oxidation by methanotrophs and effect of soil moisture content on their activity[J]. Chinese Journal of Eco-Agriculturte, 2003, 11(1): 94-97. (in Chinese)

    38 劉巽浩, 徐文修, 李增嘉, 等. 農(nóng)田生態(tài)系統(tǒng)碳足跡法:誤區(qū)、改進(jìn)與應(yīng)用—兼析中國集約農(nóng)作碳效率(續(xù))[J]. 中國農(nóng)業(yè)資源與區(qū)劃, 2014, 35(1): 1-7. LIU Xunhao, XU Wenxiu, LI Zengjia, et al. Missteps improvement and application of carbon footprint methodology in farmland ecosystems with the case study of analyzing the carbon efficient of China’s intensive farming[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2014, 35(1): 1-7. (in Chinese)

    39 岳維云, 宋建榮, 趙尚文, 等. 不同覆膜方式對地溫及冬小麥地上干物質(zhì)累積規(guī)律的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2014, 32(4): 133-139. YUE Weiyun, SONG Jianrong, ZHAO Shangwen, et al. Influence of film mulching on soil temperature and dry matter accumulation of winter wheat[J]. Agricultural Research in the Arid Areas, 2014, 32(4): 133-139. (in Chinese)

    40 李尚中, 王勇, 樊廷錄, 等. 旱地玉米不同覆膜方式的水溫及增產(chǎn)效應(yīng)[J]. 中國農(nóng)業(yè)科學(xué), 2010, 43(5): 922-931. LI Shangzhong, WANG Yong, FAN Tinglu, et al. Effects of different plastic film mulching modes on soil moisture, temperature and yield of dryland maize[J]. Scientia Agricultura Sinica, 2010, 43(5): 922-931. (in Chinese)

    41 史磊剛, 陳阜, 孔凡磊, 等. 華北平原冬小麥-夏玉米種植模式碳足跡研究[J]. 中國人口資源與環(huán)境, 2011, 21(9): 93-98. SHI Leigang, CHEN Fu, KONG Fanlei, et al. Carbon footprint of winter wheat-summer maize cropping pattern on North China Plain[J].China Population Resources and Environment, 2011, 21(9): 93-98. (in Chinese)

    42 李俊, 于強(qiáng), 孫曉敏, 等. 華北平原農(nóng)田生態(tài)系統(tǒng)碳交換及其環(huán)境調(diào)控機(jī)制[J]. 中國科學(xué):D輯:地球科學(xué), 2006, 36(A01): 210-223. LI Jun, YU Qiang, SUN Xiaomin, et al. Carbon exchange and its environmental regulation mechanism of farmland ecosystem in North China Plain[J].Science in China: Series D:Earth Sciences,2006, 36(A01): 210-223. (in Chinese)

    43 徐昔保, 楊桂山, 孫小祥. 太湖流域典型稻麥輪作農(nóng)田生態(tài)系統(tǒng)碳交換及影響因素[J]. 生態(tài)學(xué)報(bào), 2015, 35(20):6655-6665. XU Xibao, YANG Guishan, SUN Xiaoxiang. Analysis of net ecosystem CO2exchange (NEE) in the rice-wheat rotation agroecosystem of the Lake Taihu Basin,China[J]. Acta Ecologica Sinica, 2015, 35(20): 6655-6665. (in Chinese)

    Effect of Film Mulching Patterns on Carbon Sequestration in Wheat-Maize Rotation System

    FENG Hao1,2LIU Jingjing1,3ZHANG Afeng3,4ZOU Xiaoyang2CHEN Haixin3,5

    (1.CollegeofWaterResourcesandArchitecturalEngineering,NorthwestA&FUniversity,Yangling,Shaanxi712100,China2.InstituteofWaterandSoilConservation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling,Shaanxi712100,China3.ChineseNationalAcademyofWater-savingAgricultureinAridRegion,NorthwestA&FUniversity,Yangling,Shaanxi712100,China4.CollegeofNaturalResourcesandEnvironment,NorthwestA&FUniversity,Yangling,Shaanxi712100,China5.InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

    The greenhouse gases such as methane (CH4), nitrous oxide (N2O), and carbon dioxide (CO2) were increased in atmospheric concentration since 1750, which attracted more and more attention regarding climate change. Agriculture activities contributed to the increase of greenhouse gas concentration in the atmosphere, resulting in high global warming potential, which was estimated to contribute about 47%~58% of the total anthropogenic emissions of CH4and N2O. Film mulching is reported to be a significant agricultural factor which greatly affected greenhouse gas (GHG) emissions and the carbon footprint. The annual GHG emissions (CO2, CH4and N2O) from a wheat-maize rotation system during 2014—2015 were monitored by using the static opaque chamber and gas chromatography technique in Guanzhong Plain of China. Four mulching treatments were no mulching (CK), semi-film mulching (BM), ridge-furrow planting with film mulching over ridge (LM) and whole film mulching (QM). Net global warming potential (NGWP) and carbon footprint were used to evaluate the effect of film mulching on GHG emissions and composition of carbon footprint production. The results showed that compared with CK, the BM, LM and QM treatments increased annual crop yield by 9.0%, 16.5% and 26.6%, respectively. Similarly, the BM, LM and QM treatments increased annual CO2emission by 33.9%, 9.3% and 31.6% and annual N2O emission by 22.9%, 14.3% and 47.1%, respectively. However, annual CH4emissions showed no significant difference compared with CK. In addition, NGWP was increased by 9.0%, 16.7% and 26.0%, respectively. The carbon footprint of LM and QM treatments was 33.2% and 21.9% lower than that of CK, while BM treatment showed no significant difference compared with CK. Compared with carbon footprint per unit crop yield of CK, that of BM treatment was increased by 16.3%, while that of LM treatment was decreased by 13.1%. Meanwhile, QM treatment showed no significant difference. Considering both the economic and environmental effects of different mulching treatments, ridge-furrow planting with film mulching over ridge in the wheat-maize rotation system was recommended for carbon sequestration and greenhouse gas reduction in Guanzhong Plain of China with the purposes of water saving and carbon sequestration.

    wheat-maize; film mulching; greenhouse gas; net global warming potential; carbon footprint

    10.6041/j.issn.1000-1298.2017.04.024

    2016-08-08

    2016-09-20

    國家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2013AA102904)、國家自然科學(xué)基金項(xiàng)目(41301305)和西北農(nóng)林科技大學(xué)基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2014YB062、2452015355)

    馮浩(1970—),男,研究員,博士生導(dǎo)師,主要從事農(nóng)業(yè)水土資源高效利用研究,E-mail: nercwsi@vip.sina.com

    S162.4+1; X826

    A

    1000-1298(2017)04-0180-10

    猜你喜歡
    足跡覆膜通量
    冬小麥田N2O通量研究
    蘋果秋覆膜 樹體營養(yǎng)好
    成長足跡
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    基于SLS覆膜砂的無模鑄型快速制造
    足跡
    花生新品種錦花20覆膜栽培技術(shù)
    雜豆全覆膜綜合配套高產(chǎn)栽培技術(shù)
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    春的足跡
    黄色配什么色好看| 男人舔奶头视频| 日韩精品中文字幕看吧| 久久久久久久精品吃奶| 国产欧美日韩一区二区三| 精品国产三级普通话版| 草草在线视频免费看| 丰满的人妻完整版| 少妇被粗大猛烈的视频| eeuss影院久久| 黄片小视频在线播放| 我的女老师完整版在线观看| 亚洲,欧美,日韩| 白带黄色成豆腐渣| 亚洲专区国产一区二区| 成人高潮视频无遮挡免费网站| 婷婷亚洲欧美| 午夜福利在线在线| 一a级毛片在线观看| 十八禁网站免费在线| 他把我摸到了高潮在线观看| 我的女老师完整版在线观看| 亚洲18禁久久av| 国产精品日韩av在线免费观看| 午夜亚洲福利在线播放| 久久久国产成人免费| 一级作爱视频免费观看| 三级毛片av免费| 在现免费观看毛片| 在线免费观看不下载黄p国产 | 制服丝袜大香蕉在线| 午夜福利在线观看吧| 国产美女午夜福利| 97超级碰碰碰精品色视频在线观看| 亚洲第一电影网av| 日日夜夜操网爽| 亚洲欧美日韩东京热| 亚洲成人免费电影在线观看| 日韩精品青青久久久久久| 丰满的人妻完整版| 久久久精品欧美日韩精品| 欧美成人一区二区免费高清观看| 日韩有码中文字幕| 极品教师在线视频| 午夜日韩欧美国产| 麻豆久久精品国产亚洲av| 日韩成人在线观看一区二区三区| 香蕉av资源在线| 国产一区二区激情短视频| 欧美日韩乱码在线| 精品午夜福利视频在线观看一区| 亚洲色图av天堂| eeuss影院久久| av视频在线观看入口| 婷婷丁香在线五月| 99视频精品全部免费 在线| 国产蜜桃级精品一区二区三区| 久久亚洲真实| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 亚洲自偷自拍三级| 可以在线观看的亚洲视频| av专区在线播放| 欧美一区二区精品小视频在线| 变态另类丝袜制服| 99久久精品国产亚洲精品| 国产熟女xx| 波多野结衣高清无吗| 综合色av麻豆| 亚洲成a人片在线一区二区| 中文资源天堂在线| av在线观看视频网站免费| 极品教师在线免费播放| 亚洲久久久久久中文字幕| 成人国产综合亚洲| 搡老岳熟女国产| 国产国拍精品亚洲av在线观看| 午夜福利免费观看在线| 日韩高清综合在线| 91麻豆av在线| 少妇高潮的动态图| 亚洲18禁久久av| 麻豆成人午夜福利视频| 色5月婷婷丁香| 午夜影院日韩av| 日本免费a在线| 人妻久久中文字幕网| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 精品久久久久久久久av| 欧美性猛交╳xxx乱大交人| 免费观看人在逋| 97热精品久久久久久| aaaaa片日本免费| 内地一区二区视频在线| 天堂av国产一区二区熟女人妻| 亚洲经典国产精华液单 | 真人一进一出gif抽搐免费| 国产成人欧美在线观看| 亚洲三级黄色毛片| 亚洲av五月六月丁香网| 最好的美女福利视频网| 国产高清视频在线观看网站| 免费看a级黄色片| 男人舔女人下体高潮全视频| 亚洲精品色激情综合| 少妇裸体淫交视频免费看高清| 欧美色欧美亚洲另类二区| 成人永久免费在线观看视频| 色噜噜av男人的天堂激情| 欧美成人一区二区免费高清观看| 真人做人爱边吃奶动态| eeuss影院久久| 精品久久久久久久末码| 成人欧美大片| 国产av在哪里看| 嫩草影院新地址| 校园春色视频在线观看| 在线观看66精品国产| 国产成年人精品一区二区| 男人和女人高潮做爰伦理| 麻豆av噜噜一区二区三区| 男女那种视频在线观看| 亚洲第一区二区三区不卡| 国产黄a三级三级三级人| 欧美成狂野欧美在线观看| 蜜桃亚洲精品一区二区三区| 久久国产精品人妻蜜桃| 欧美潮喷喷水| 国产视频一区二区在线看| 精品久久久久久久久亚洲 | 欧美日韩乱码在线| 国产美女午夜福利| 日本一本二区三区精品| 黄色日韩在线| www日本黄色视频网| 国产精品自产拍在线观看55亚洲| av女优亚洲男人天堂| 高潮久久久久久久久久久不卡| 国产乱人伦免费视频| 男人舔女人下体高潮全视频| 青草久久国产| 欧美激情国产日韩精品一区| 亚洲av.av天堂| 麻豆成人午夜福利视频| 黄色一级大片看看| 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av| 波多野结衣高清作品| 很黄的视频免费| 人妻制服诱惑在线中文字幕| 亚洲自偷自拍三级| 91麻豆av在线| 国产爱豆传媒在线观看| 亚洲aⅴ乱码一区二区在线播放| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 制服丝袜大香蕉在线| 99久久成人亚洲精品观看| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 最近最新中文字幕大全电影3| 91久久精品电影网| 嫩草影院入口| 午夜日韩欧美国产| 热99在线观看视频| 亚洲欧美激情综合另类| 在线a可以看的网站| 亚洲熟妇熟女久久| 婷婷丁香在线五月| 亚洲专区中文字幕在线| 精品久久久久久久末码| 国产精品亚洲美女久久久| 99久久精品热视频| 97超级碰碰碰精品色视频在线观看| 看十八女毛片水多多多| a级毛片a级免费在线| 精品日产1卡2卡| 亚洲电影在线观看av| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 亚洲 国产 在线| 91在线精品国自产拍蜜月| 观看免费一级毛片| 欧美日本亚洲视频在线播放| 真人一进一出gif抽搐免费| 国产伦在线观看视频一区| 男女做爰动态图高潮gif福利片| 亚洲最大成人手机在线| 国内精品一区二区在线观看| 亚洲av免费在线观看| 亚洲成a人片在线一区二区| 日本成人三级电影网站| 国产精品女同一区二区软件 | 欧美最新免费一区二区三区 | 久久国产精品人妻蜜桃| 少妇人妻一区二区三区视频| 可以在线观看的亚洲视频| 黄色一级大片看看| 乱人视频在线观看| 亚洲av成人精品一区久久| 午夜福利18| 免费观看的影片在线观看| 亚洲国产日韩欧美精品在线观看| 亚州av有码| 桃色一区二区三区在线观看| 人妻夜夜爽99麻豆av| 少妇裸体淫交视频免费看高清| 久久人妻av系列| 最好的美女福利视频网| 国产人妻一区二区三区在| 免费看光身美女| 国产精品一区二区三区四区久久| 欧美国产日韩亚洲一区| 欧美日本亚洲视频在线播放| 最近在线观看免费完整版| 天天一区二区日本电影三级| 亚洲熟妇熟女久久| 尤物成人国产欧美一区二区三区| 搡老岳熟女国产| 久久久国产成人免费| 国产精品一区二区三区四区免费观看 | 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线观看免费| 国产精品三级大全| 久9热在线精品视频| 亚洲人成网站在线播| 一个人免费在线观看的高清视频| aaaaa片日本免费| 男女那种视频在线观看| 国产精品98久久久久久宅男小说| 国产精品影院久久| a级毛片a级免费在线| 男人狂女人下面高潮的视频| 国产乱人伦免费视频| 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 亚洲精品乱码久久久v下载方式| 国产免费一级a男人的天堂| 黄色视频,在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩高清专用| 九色国产91popny在线| 亚洲中文字幕一区二区三区有码在线看| 热99在线观看视频| 亚洲欧美日韩卡通动漫| 欧洲精品卡2卡3卡4卡5卡区| av天堂中文字幕网| 美女黄网站色视频| 精品欧美国产一区二区三| 国产伦精品一区二区三区四那| 国产精品久久久久久亚洲av鲁大| 偷拍熟女少妇极品色| 老司机午夜福利在线观看视频| 国内精品一区二区在线观看| avwww免费| 99riav亚洲国产免费| 淫妇啪啪啪对白视频| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 日韩欧美国产在线观看| 1000部很黄的大片| 国产精华一区二区三区| 波野结衣二区三区在线| 欧美在线一区亚洲| 内射极品少妇av片p| 九色成人免费人妻av| 欧美高清性xxxxhd video| 日韩欧美国产一区二区入口| 亚洲av免费高清在线观看| 国产毛片a区久久久久| 欧美性猛交黑人性爽| 九色国产91popny在线| 久久精品国产自在天天线| 中文字幕精品亚洲无线码一区| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 国产精品1区2区在线观看.| 制服丝袜大香蕉在线| 亚洲av中文字字幕乱码综合| 99在线人妻在线中文字幕| 级片在线观看| 男人舔女人下体高潮全视频| 精品日产1卡2卡| av在线观看视频网站免费| 久久久国产成人免费| 欧美日韩综合久久久久久 | 亚洲成人久久性| 少妇高潮的动态图| 成人国产一区最新在线观看| 国产欧美日韩精品亚洲av| 欧美成狂野欧美在线观看| 啦啦啦观看免费观看视频高清| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 国产精品国产高清国产av| 欧美性猛交╳xxx乱大交人| 99精品在免费线老司机午夜| 波多野结衣高清作品| 精品人妻一区二区三区麻豆 | 国产av麻豆久久久久久久| 亚洲美女视频黄频| 婷婷丁香在线五月| 欧美在线黄色| 精品午夜福利视频在线观看一区| 亚洲成人中文字幕在线播放| 免费在线观看亚洲国产| av天堂在线播放| 一级作爱视频免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品嫩草影院av在线观看 | 97碰自拍视频| 精品熟女少妇八av免费久了| 国内少妇人妻偷人精品xxx网站| 18禁裸乳无遮挡免费网站照片| 欧美不卡视频在线免费观看| 色播亚洲综合网| 亚洲va日本ⅴa欧美va伊人久久| 丰满乱子伦码专区| 亚洲中文日韩欧美视频| 国产一区二区亚洲精品在线观看| 丁香六月欧美| 国产国拍精品亚洲av在线观看| 国产又黄又爽又无遮挡在线| 亚洲激情在线av| 日本三级黄在线观看| 中文字幕人成人乱码亚洲影| 精品99又大又爽又粗少妇毛片 | 免费人成视频x8x8入口观看| 欧美激情久久久久久爽电影| 欧美日韩乱码在线| 成人永久免费在线观看视频| 观看美女的网站| 亚洲午夜理论影院| 欧美性猛交╳xxx乱大交人| 中文字幕av在线有码专区| 长腿黑丝高跟| 免费一级毛片在线播放高清视频| 色哟哟·www| 国产探花在线观看一区二区| 国产黄a三级三级三级人| 高潮久久久久久久久久久不卡| 日韩欧美三级三区| 天堂网av新在线| 美女黄网站色视频| 欧美一级a爱片免费观看看| 美女黄网站色视频| 亚洲在线观看片| 深夜精品福利| 成人一区二区视频在线观看| 成年女人永久免费观看视频| 看黄色毛片网站| 欧美性猛交黑人性爽| 真人做人爱边吃奶动态| 久久久精品欧美日韩精品| 色综合婷婷激情| 丰满乱子伦码专区| 91九色精品人成在线观看| 亚洲性夜色夜夜综合| 亚洲国产日韩欧美精品在线观看| 日本免费a在线| 熟妇人妻久久中文字幕3abv| 91在线精品国自产拍蜜月| 亚洲欧美日韩东京热| 69人妻影院| 人妻制服诱惑在线中文字幕| 久久亚洲真实| 中文亚洲av片在线观看爽| 亚洲精华国产精华精| 亚洲精品一卡2卡三卡4卡5卡| 男女视频在线观看网站免费| 在线观看66精品国产| 久久午夜福利片| 变态另类成人亚洲欧美熟女| 国产精华一区二区三区| 91久久精品电影网| 免费在线观看亚洲国产| 亚洲国产精品sss在线观看| 亚洲精品456在线播放app | 女人十人毛片免费观看3o分钟| 国内精品一区二区在线观看| 精品久久久久久久末码| 精品国产亚洲在线| 日韩有码中文字幕| 国产男靠女视频免费网站| 最新在线观看一区二区三区| 精品久久久久久久久亚洲 | av福利片在线观看| 在线观看午夜福利视频| 久久精品91蜜桃| 一本一本综合久久| 亚洲精品亚洲一区二区| 国产69精品久久久久777片| 亚洲自拍偷在线| 中出人妻视频一区二区| 久久99热这里只有精品18| 级片在线观看| 国产av不卡久久| 国产精品野战在线观看| 99国产综合亚洲精品| 欧美+日韩+精品| 51午夜福利影视在线观看| 日日摸夜夜添夜夜添av毛片 | 一级av片app| eeuss影院久久| 亚洲美女搞黄在线观看 | 少妇丰满av| 日韩欧美在线二视频| 亚洲av.av天堂| 99热这里只有是精品50| 波多野结衣高清作品| 色哟哟·www| 色综合婷婷激情| 欧美一区二区亚洲| 国产白丝娇喘喷水9色精品| 中亚洲国语对白在线视频| 在线a可以看的网站| 久久性视频一级片| 国产乱人视频| 国产69精品久久久久777片| 亚洲成人中文字幕在线播放| 国产av一区在线观看免费| 丰满人妻一区二区三区视频av| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| 久久久久久久久中文| 亚洲av一区综合| 好看av亚洲va欧美ⅴa在| 床上黄色一级片| 国产 一区 欧美 日韩| 日韩欧美一区二区三区在线观看| 国产成人av教育| 中文字幕av成人在线电影| 男插女下体视频免费在线播放| 亚洲最大成人中文| 精品久久国产蜜桃| 最近最新中文字幕大全电影3| 男人舔女人下体高潮全视频| 最新在线观看一区二区三区| 大型黄色视频在线免费观看| 亚洲成av人片免费观看| 欧美最新免费一区二区三区 | 久久伊人香网站| 美女免费视频网站| 免费在线观看亚洲国产| 在线观看av片永久免费下载| 午夜精品在线福利| 国产老妇女一区| 免费搜索国产男女视频| 欧美乱色亚洲激情| 久久精品综合一区二区三区| 日韩欧美在线二视频| 久久精品久久久久久噜噜老黄 | 美女cb高潮喷水在线观看| 亚洲av二区三区四区| 成人午夜高清在线视频| 老司机深夜福利视频在线观看| 亚洲av电影在线进入| 精品国内亚洲2022精品成人| 男女下面进入的视频免费午夜| 99久久精品国产亚洲精品| 99久久久亚洲精品蜜臀av| 性色avwww在线观看| 精品人妻视频免费看| www.色视频.com| 亚洲精品456在线播放app | 欧美精品啪啪一区二区三区| 免费大片18禁| 国产精品乱码一区二三区的特点| 最后的刺客免费高清国语| 一进一出抽搐gif免费好疼| 亚洲欧美激情综合另类| 亚洲美女视频黄频| 亚洲欧美日韩高清在线视频| 18禁在线播放成人免费| 亚洲经典国产精华液单 | 亚洲欧美日韩无卡精品| 国产av一区在线观看免费| 国产黄a三级三级三级人| 特级一级黄色大片| 嫁个100分男人电影在线观看| 大型黄色视频在线免费观看| 成人亚洲精品av一区二区| 宅男免费午夜| 国产精品不卡视频一区二区 | 12—13女人毛片做爰片一| 色视频www国产| 精品一区二区免费观看| 窝窝影院91人妻| 嫩草影院新地址| 国产精品98久久久久久宅男小说| 精品熟女少妇八av免费久了| 少妇被粗大猛烈的视频| 日韩大尺度精品在线看网址| 男女之事视频高清在线观看| 午夜两性在线视频| 在线看三级毛片| 久久婷婷人人爽人人干人人爱| 精品国内亚洲2022精品成人| 国产高清有码在线观看视频| 亚洲精品在线美女| 欧美色欧美亚洲另类二区| 99久久精品热视频| 中文字幕高清在线视频| 丝袜美腿在线中文| 欧美一区二区精品小视频在线| 如何舔出高潮| 久久精品影院6| 两性午夜刺激爽爽歪歪视频在线观看| 毛片女人毛片| 日本五十路高清| 直男gayav资源| 琪琪午夜伦伦电影理论片6080| 国产毛片a区久久久久| 亚洲欧美日韩高清专用| 18禁黄网站禁片免费观看直播| 国产精品一区二区免费欧美| 免费观看人在逋| 亚洲电影在线观看av| 久久99热6这里只有精品| 亚洲av.av天堂| 在线播放无遮挡| 欧美激情在线99| 免费在线观看影片大全网站| 欧美色视频一区免费| 丁香欧美五月| 精品一区二区三区av网在线观看| 我要搜黄色片| 一区二区三区激情视频| 又爽又黄无遮挡网站| 日韩欧美三级三区| 欧美激情国产日韩精品一区| 毛片女人毛片| 一级a爱片免费观看的视频| 天堂网av新在线| 日本a在线网址| 亚洲一区二区三区色噜噜| 亚洲第一欧美日韩一区二区三区| 十八禁网站免费在线| 757午夜福利合集在线观看| 午夜精品久久久久久毛片777| 婷婷色综合大香蕉| 可以在线观看毛片的网站| 色播亚洲综合网| 中文字幕av成人在线电影| 久久精品国产亚洲av天美| 丰满乱子伦码专区| 一级黄色大片毛片| 精品无人区乱码1区二区| 色综合亚洲欧美另类图片| 在线播放无遮挡| 男插女下体视频免费在线播放| 村上凉子中文字幕在线| 夜夜躁狠狠躁天天躁| 久久6这里有精品| 老熟妇仑乱视频hdxx| 亚洲精品乱码久久久v下载方式| 国产亚洲欧美在线一区二区| 在线观看66精品国产| 两人在一起打扑克的视频| 亚州av有码| 欧美精品啪啪一区二区三区| 熟女人妻精品中文字幕| 国产精品av视频在线免费观看| 亚洲av中文字字幕乱码综合| 97热精品久久久久久| 欧美日本亚洲视频在线播放| 亚洲欧美日韩高清在线视频| 免费人成视频x8x8入口观看| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 免费观看人在逋| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品乱码久久久v下载方式| or卡值多少钱| 动漫黄色视频在线观看| 亚洲av电影不卡..在线观看| 午夜福利免费观看在线| 亚洲国产高清在线一区二区三| 中出人妻视频一区二区| 久久久久久国产a免费观看| 日本成人三级电影网站| 又爽又黄无遮挡网站| 少妇的逼好多水| 波野结衣二区三区在线| 成年免费大片在线观看| 国产国拍精品亚洲av在线观看| 精品福利观看| 欧美色欧美亚洲另类二区| 精品午夜福利在线看| 国产精品一区二区三区四区免费观看 | 宅男免费午夜| 免费看光身美女| 国产麻豆成人av免费视频| 免费电影在线观看免费观看| 757午夜福利合集在线观看| 美女 人体艺术 gogo| 97超视频在线观看视频| 成人高潮视频无遮挡免费网站| 午夜激情欧美在线| 99热只有精品国产| 欧美不卡视频在线免费观看| xxxwww97欧美| 男女床上黄色一级片免费看| 人妻制服诱惑在线中文字幕| 男插女下体视频免费在线播放| 嫩草影视91久久| 欧美极品一区二区三区四区| 日本精品一区二区三区蜜桃| 中文字幕高清在线视频| 免费观看人在逋| 欧美最新免费一区二区三区 | 日本一本二区三区精品| 国产三级黄色录像|