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

    草地農(nóng)業(yè)生態(tài)系統(tǒng)的碳平衡分析方法

    2017-02-08 05:57:04婁珊寧陳先江侯扶江
    生態(tài)學報 2017年2期
    關(guān)鍵詞:草業(yè)家畜牧場

    婁珊寧,陳先江,侯扶江

    蘭州大學草地農(nóng)業(yè)科技學院, 草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點實驗室, 蘭州 730020

    草地農(nóng)業(yè)生態(tài)系統(tǒng)的碳平衡分析方法

    婁珊寧,陳先江,侯扶江*

    蘭州大學草地農(nóng)業(yè)科技學院, 草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點實驗室, 蘭州 730020

    根據(jù)草地農(nóng)業(yè)生態(tài)系統(tǒng)的結(jié)構(gòu),它的碳平衡為4個生產(chǎn)層的碳平衡之和,也是3個界面的碳平衡之和,而某一生產(chǎn)層或者某一界面的碳平衡則是其固定、輸入、排放和輸出的碳之和。草地農(nóng)業(yè)生態(tài)系統(tǒng)4個生產(chǎn)層的碳平衡分析方法定量重要生產(chǎn)環(huán)節(jié)的碳匯與碳源過程,便于草業(yè)生產(chǎn)改進碳匯管理;草地農(nóng)業(yè)生態(tài)系統(tǒng)3個界面的碳平衡分析方法顯示碳源和碳匯的發(fā)生機理,及其空間和數(shù)量關(guān)系,便于調(diào)控草業(yè)生產(chǎn)組分以增匯減排;但是,這兩個方法不易區(qū)分碳的來源和去向,難以明確其利用效率。草地農(nóng)業(yè)生態(tài)系統(tǒng)碳平衡分析的輸入/輸出法定量地指示碳的來源和去向,以及碳效率,計算簡單,但是較為概括,不利于牧場尺度的草業(yè)碳匯管理。以中國祁連山甘肅馬鹿牧場和澳大利亞塔斯瑪尼亞奶牛牧場為例,用3種方法分析了兩個牧場的碳平衡,結(jié)果表明,放牧管理的草業(yè)系統(tǒng)的主要碳源是休閑旅游、產(chǎn)品加工流通環(huán)節(jié)產(chǎn)生的溫室氣體,主要碳匯是草地和土壤中貯存的碳,好的草地管理可以增匯減排。

    草地農(nóng)業(yè)生態(tài)系統(tǒng);生產(chǎn)層;界面;碳平衡;放牧;溫室氣體

    草地是全球面積最大的陸地生態(tài)系統(tǒng),全世界天然草地有機碳儲量7610—10730億t[1-2],占陸地碳儲量的37.1%—52.3%[3]。在一定時間尺度上,草地為碳中性,年際間波動較小[4]。過去20年,中國草地生物量和土壤有機碳庫沒有顯著變化[5-6],處于中性碳匯狀態(tài)[7];期間,它輸出的旅游產(chǎn)品、文化產(chǎn)品和草畜產(chǎn)品卻逐年增多[8-9];如果把這部分碳損失計算在內(nèi),草地本質(zhì)上是一個碳匯。草地碳庫受管理的影響較大[10-11],碳源或碳匯取決于放牧管理水平[12-13],一般放牧系統(tǒng)的溫室氣體(Greenhouse gases,GHG)排放低于同等畜群規(guī)模的舍飼系統(tǒng)[14]。放牧系統(tǒng)能夠數(shù)十年(如英國和新西蘭的栽培草地)乃至上千年(如天然草地)持續(xù)輸出畜產(chǎn)品,草畜耦合可以提升這種能力[15]。

    草業(yè)碳匯管理以碳測定技術(shù)為基礎(chǔ)。以往較為關(guān)注土壤和家畜GHG排放測定,多以某一組分為對象,而非生態(tài)系統(tǒng)整體[16]。土壤GHG的測定有間接測定法和直接測定法[17]。間接法通過測定土壤ATP含量等相關(guān)因子,建立其與土壤GHG排放的數(shù)量預(yù)測模型[18],適用于大尺度或跨區(qū)域的土壤呼吸預(yù)測,以及直接測定土壤呼吸困難或成本較高的領(lǐng)域。直接測定法常用靜態(tài)氣室法(Static Chamber Method)、動態(tài)氣室法和微氣象技術(shù)、室內(nèi)測定法等,測定容器內(nèi)GHG濃度的變化來計算其通量[17],與成分綜合法、底物誘導(dǎo)呼吸法、同位素稀釋法、13C自然豐度法、脈沖標記法等結(jié)合能測定土壤GHG的來源,進而定量根呼吸和土壤微生物呼吸[19]。家畜GHG排放主要是CH4和N2O,放牧家畜呼出的CO2被認為是對草地光合作用的平衡而不計測,CH4主要通過腸道發(fā)酵排放,N2O主要由排泄物排放[20]。家畜GHG排放的測定包括呼吸代謝室(箱)、呼吸面罩和呼吸面具、微氣象技術(shù)、隧道技術(shù)、激光技術(shù)、模型法、體外模擬法,其中排泄物GHG測定基本同土壤。呼吸室法將家畜置于艙室,根據(jù)室內(nèi)GHG變化,計算家畜GHG排放量,有密閉式、開閉式和開放式3種[20]。呼吸面罩常與六氟化硫(SF6)等指示物結(jié)合,在家畜鼻孔處采集氣體樣品,不限制家畜活動,適用于放牧試驗。微氣象技術(shù)、激光技術(shù)、隧道技術(shù)可測定放牧畜群的GHG的排放量[20-22]。大中尺度家畜GHG排放多根據(jù)IPCC的參數(shù)和統(tǒng)計數(shù)據(jù)計算[23]。

    草地農(nóng)業(yè)生態(tài)系統(tǒng)(草業(yè)系統(tǒng))以草地資源為基礎(chǔ),在生物因素、環(huán)境因素和人類社會管理因素的共同作用下發(fā)生與發(fā)展,其中土-草-畜-人是碳流通的主干,它具有前植物(Pre-plant production level, PPP)、植物(Plant production level, PP)、動物(Animal production level, AP)和后生物(Post-biotic production level, PBP)4個生產(chǎn)層,以及草叢-地境(Interface between herbage and site, HIS; interface A, IA)、草地-家畜(Interface between grassland and livestock, IGL; interface B, IB)、草畜系統(tǒng)-社會經(jīng)營管理(Interface between grassland-livestock system and social and economic management, IGLSE; interface C, IC)3個界面[24]。它的碳循環(huán)從CO2同化開始,到碳排放和人類享用各類產(chǎn)品結(jié)束,其碳平衡機制取決于草地利用方式及管理水平,系統(tǒng)的核心驅(qū)動力是草-畜互作,動物生產(chǎn)是源、匯轉(zhuǎn)換的重要閥門[10]。目前,草地碳循環(huán)的研究較多,特別是草地植被、土壤有機碳,對家畜生產(chǎn)及其與草原碳匯的關(guān)系、草畜產(chǎn)品加工、流通等過程碳平衡的報道較少。

    當前,草業(yè)系統(tǒng)碳匯的研究結(jié)果差異較大[25-28],生產(chǎn)管理對碳動態(tài)的影響及其機制尚不明確,限制了人類對草業(yè)碳匯過程的理解和管理,主要因為估算方法不統(tǒng)一及其導(dǎo)致的統(tǒng)計數(shù)據(jù)差異[25-26],根本原因是缺乏以草業(yè)系統(tǒng)結(jié)構(gòu)和關(guān)鍵生態(tài)過程為基礎(chǔ)的評估方法。為此,本研究擬從草業(yè)系統(tǒng)4個生產(chǎn)層(結(jié)構(gòu))、3個界面(關(guān)鍵生態(tài)過程)和系統(tǒng)整體,重視家畜和草畜互作,建立簡單、適用、準確的草業(yè)系統(tǒng)碳平衡分析方法。

    文中,草地農(nóng)業(yè)生態(tài)系統(tǒng)的碳平衡(Carbon balance of grassland agro-ecosystem, CBGAE)指一定時間內(nèi)碳貯量的變化,主要指碳,兼顧氮,因為N2O是草業(yè)系統(tǒng)重要的GHG。計量單位為所含的碳量或相當?shù)奶籍斄俊?/p>

    1 草地農(nóng)業(yè)生態(tài)系統(tǒng)4個生產(chǎn)層的碳流動

    草業(yè)系統(tǒng)的碳沿著前植物、植物、動物、后生物4個生產(chǎn)層運動(圖1)。

    圖1 草地農(nóng)業(yè)生態(tài)系統(tǒng)4個生產(chǎn)層的碳流Fig.1 Carbon flow of grassland agro-ecosystem based on four production levels實線表示碳在4個生產(chǎn)層內(nèi)流動或者由生產(chǎn)層流動至外界,虛線表示外界將碳輸送至4個生產(chǎn)層中

    1.1 前植物生產(chǎn)層的碳平衡

    前植物生產(chǎn)層以水土保持、休閑娛樂、狩獵旅游等景觀生產(chǎn)為主[29]。塵降、水土保持、生物固氮等活動積累碳氮,是碳匯過程。草地旅游、狩獵等景觀、休閑產(chǎn)品是碳源過程,包括生產(chǎn)層管理的能源消耗、草地踐踏等引起的碳排放(圖1)。該生產(chǎn)層的GHG排放研究雖少,但是草業(yè)觀光旅游的碳排放已受到關(guān)注[30]。塵降是碳匯,但短期內(nèi)貢獻較小。

    前植物生產(chǎn)層的碳平衡CBPPP(carbon balance of PPP)=Cinput-Coutput。Cinput為碳輸入,包括碳氮塵降和土壤沉積、生物固氮及施肥、補播等碳輸入。Coutput為碳輸出,主要包括人和草地以及機械、灌溉等管理活動的GHG排放。CBPPP>0,生產(chǎn)層為碳匯;CBPPP=0,生產(chǎn)層為碳中性,處于生態(tài)臨界點,面臨退化風險;CBPPP<0,生產(chǎn)層為碳源,處于退化階段[10]。提高單位碳排放的經(jīng)濟效益是該生產(chǎn)層增碳減排的主要方向。

    1.2 植物生產(chǎn)層的碳平衡

    植物生產(chǎn)層是傳統(tǒng)意義上天然草地和栽培草地的牧草生產(chǎn),是草業(yè)系統(tǒng)的主要碳庫,包括植物碳庫和土壤碳庫。碳匯過程主要是光合作用和生物固氮的碳氮同化,以及施肥、補播等碳輸入。碳源主要是土、草和生產(chǎn)層管理中人力、機械、能量消耗等GHG排放,以及牧草輸出等(圖1)。

    植物生產(chǎn)層的碳平衡CBPP(carbon balance of PP)=Cinput-Coutput。一般,CBPP>0,尤其是天然草地。CBPP=0,表明天然草地面臨退化風險,栽培草地需要水、肥、能量等投入以維持運轉(zhuǎn)因而屬正常。CBPP<0,草地碳貯量正在減少,處于退化狀態(tài)[13,31]。提高牧草產(chǎn)量和質(zhì)量、減少土壤碳排放和管理的投入是該生產(chǎn)層的主要增匯減排途徑。

    1.3 動物生產(chǎn)層的碳平衡

    動物生產(chǎn)層是傳統(tǒng)意義上的食草家畜飼養(yǎng),是草業(yè)系統(tǒng)的主要碳源之一。植物生產(chǎn)層的植物碳,絕大部分通過家畜腸道發(fā)酵、排泄物揮發(fā)和分解而排放到環(huán)境(圖1)。畜群管理中人力、機械、疫病防控等排放GHG,是另一主要碳源。家畜生產(chǎn)因購入飼草料而輸入碳。家畜排泄物除了淋溶、GHG排放,大部分碳氮返還到草地,成為土壤碳。極少部分牧草碳固定在畜產(chǎn)品中,而且家畜骨骼、皮毛、角蹄中的碳易于長期保存,也是碳庫。畜肉血中的碳很快排放,排放量可根據(jù)屠宰率計算。

    動物生產(chǎn)層的碳平衡CBAP(carbon balance of AP)=Cinput-Coutput。Cinput主要是牧草。家畜對牧草轉(zhuǎn)化效率遠小于1,故CBAP<0。提高牧草尤其是N的轉(zhuǎn)化效率、促進排泄物循環(huán)利用是該生產(chǎn)層增碳減排的主要途徑。

    1.4 后生物生產(chǎn)層的碳平衡

    后生物生產(chǎn)層主要是草畜產(chǎn)品加工、貯運以及整個草業(yè)系統(tǒng)管理,能夠大幅提升草業(yè)生產(chǎn)的社會效益和經(jīng)濟效益[32]。同時,需要投入大量能量和物質(zhì)。碳源包括產(chǎn)品加工、貯運等,及其基礎(chǔ)設(shè)施和設(shè)備制造的碳排放。牧草和家畜作為加工原料輸入該生產(chǎn)層,部分草畜產(chǎn)品可長久保存,為碳匯。

    后生物生產(chǎn)層的碳平衡CBPBP(carbon balance of PBP)=Cinput-Coutput。因生產(chǎn)層運轉(zhuǎn)所排放的碳遠多于輸入的碳,CBAP<0。草業(yè)系統(tǒng)如果以經(jīng)濟效益為主要目標,兼顧食物生產(chǎn),可以通過該生產(chǎn)層提升單位碳排放的經(jīng)濟效益,或可總體上減少GHG排放;如果以獲取食物為主要目標,兼顧經(jīng)濟效益,應(yīng)該控制草畜產(chǎn)品加工、貯運的碳排放;可見,該生產(chǎn)層的管理有助于草業(yè)系統(tǒng)整體的增匯減排。

    1.5 系統(tǒng)碳平衡與案例分析

    草地農(nóng)業(yè)生態(tài)系統(tǒng)的碳平衡是4個生產(chǎn)層碳平衡之和,CBGAE=CBPPP+CBPP+CBAP+CBPBP。在某一生產(chǎn)周期,CBGAE>0,草業(yè)系統(tǒng)為碳匯,反之則為碳源。它的碳效益為碳固定量與經(jīng)濟效益的比,值越大,說明生態(tài)系統(tǒng)的物質(zhì)和能量轉(zhuǎn)化效率越高,環(huán)境負債率越低,生態(tài)和經(jīng)濟可持續(xù)性越強。結(jié)合各生產(chǎn)層的碳平衡分析,可以針對性地改進碳管理。

    CBGAE的測算可以基于某一農(nóng)牧場或某一企業(yè)、某一生產(chǎn)區(qū)域、生態(tài)區(qū)域或行政區(qū)域,故以我國祁連山中段的甘肅馬鹿牧場和澳大利亞塔斯瑪尼亞西北部的奶牛牧場為例(表1)[15,31-39]。

    祁連山牧場的甘肅馬鹿(CervuselaphuskansuensisPocock)冬季和春秋季主要在高寒草原放牧,夏季主要在高寒灌叢草甸放牧[40-41],并且發(fā)展了高山牧場旅游、中草藥和馬鹿茸、血產(chǎn)品加工[33]。塔斯瑪尼亞牧場的奶牛常年在多年生黑麥草(LoliumperenneL.)栽培草地放牧,發(fā)展了農(nóng)業(yè)觀光和奶酪加工[38]。這兩個系統(tǒng)均有相對完整的4個生產(chǎn)層,體現(xiàn)了世界農(nóng)業(yè)系統(tǒng)多功能性增強的演化趨勢[15]。

    表1 祁連山牧場馬鹿牧場(DRQ)和塔斯馬尼亞奶牛牧場(DRT)4個生產(chǎn)層的碳平衡

    Table 1 Carbon balance analysis of four production levels from the deer ranch in Qilian Mountain (DRQ) and dairy ranch in Tasmania (DRT)

    生產(chǎn)層Productionlevel主要碳源(1)Maincarbonsource/(tC/a)主要碳匯Maincarbonsink/(tC/a)平移碳平衡Carbonbalance/(tC/a)項目ItemsDRQDRT項目ItemsDRQDRTDRQDRT前植物生產(chǎn)層旅游人數(shù)/(人次/a)50006000碳塵降,水土保持/(t/a)-36000.0-103809.8Pre-plantproductionlevel淋溶/(gm-1a-1)00.7人均碳排放/(t/a)7.217.3植物生產(chǎn)層Plantproductionlevel草地面積/hm210201400SOC積累/(t/m2)0.1220.1221243298.4167580.0動物生產(chǎn)層Animalproductionlevel土壤GHG排放/(kghm-2a-1)10802300日采食量/(kg頭-1d-1)3.1520-15.8-173.3家畜數(shù)量/頭10002100消化率/%5770糞便發(fā)酵系數(shù)/(kg頭-1d-1)0.2229年排糞量/(kg頭-1d-1)5265840后生物生產(chǎn)層腸道發(fā)酵系數(shù)/(kg頭-1d-1)2081屠宰率/%50.7249.00-1005.2-2475.8Post-biological燃煤量/(t/a)108productionlevel燃油的消耗量/(t/a)1015勞動力/(人/a)1010CBGAE(2)1206277.461121.2

    DRQ:祁連山牧場馬鹿牧場 Deer ranch in Qilian Mountain;DRT:塔斯馬尼亞奶牛牧場Dairy ranch in Tasmania;CBGAE:草地農(nóng)業(yè)生態(tài)系統(tǒng)碳平衡 Carbon balance of grassland ago-ecosystem

    祁連山牧場CBGAE>0,為碳匯。碳輸出中,前植物、動物和后生物3個生產(chǎn)層分別占97.2%、0.04%和2.7%(表1)。動物生產(chǎn)層中,43%的C隨排泄物返還草地,57%排放到環(huán)境中。后生物生產(chǎn)層的碳排放主要來自茸、血加工。馬鹿的消化率低導(dǎo)致更多的GHG排放,除了食性外[33],主要因為冬季牧場牧草粗老而且利用時間長,夏季和秋季放牧也不在牧草幼嫩期[40]。

    塔斯馬尼亞牧場CBGAE>0,也為碳匯。植物生產(chǎn)層因為施肥、灌溉等,碳排放多于祁連山牧場(表1)。由于劃區(qū)輪牧,奶??偸窃谀敛葸m宜期采食,因而消化率高。兩個牧場后生物生產(chǎn)層的碳源占比低于3%。

    2 草地農(nóng)業(yè)生態(tài)系統(tǒng)3個界面的碳平衡分析

    草業(yè)系統(tǒng)的3個界面是碳固定或排放的活躍區(qū)域和通道(圖2)。

    圖2 草地農(nóng)業(yè)生態(tài)系統(tǒng)3個界面的碳流Fig.2 Carbon flow of grassland agro-ecosystem based on three interfaces

    界面A過程形成草地[42],它是草類植物與非生物環(huán)境之間碳交換的區(qū)域(圖2)。其碳平衡(Carbon balance of IHS,CBIHS,CBIA)CBIA=Cinput-Coutput。Cinput主要是光合作用、牧草碳氮吸收、大氣碳氮塵降、生物固氮等碳輸入。Coutput包括牧草和土壤GHG排放,GHG主要是CO2,CH4和N2O較少。一般,CBIA>0,該界面是碳匯過程;但是,退化草地的碳等物質(zhì)及其所含能量入不敷出[43],CBIA≤0。

    界面B過程形成草畜系統(tǒng),它是草地與家畜相互作用的區(qū)域(圖2)。該界面的碳平衡(Carbon balance of IGL,CBIGL,CBIB)CBIB=Cinput-Coutput。Coutput主要是家畜腸道發(fā)酵、排泄物的揮發(fā)和分解,GHG以CH4和N2O為主。Cinput主要是凋落物和家畜排泄物返還到草地,毛、皮、骨骼的碳保存。一般,家畜采食的牧草超過20%、甚至50%作為排泄物返還草地,其中一部分揮發(fā)、淋溶,或分解中損失,故CBIB<0或CBIB≈0,界面B是碳源。

    界面C過程進化完整地草地農(nóng)業(yè)生態(tài)系統(tǒng)[42],它的碳平衡(Carbon balance of IGLSE,CBIGLSS,CBIC)CBIC=Cinput-Coutput。Coutput包括草地管理、畜群管理,以及草畜產(chǎn)品加工和貯運、旅游和狩獵等過程的碳排放。Cinput主要是固定在草畜產(chǎn)品的碳,以及草業(yè)系統(tǒng)管理、保護和恢復(fù)所增加的碳。其碳平衡較為復(fù)雜:管理得當,CBIC>0;管理不當,碳排放加速,CBIC<0。

    草業(yè)系統(tǒng)的碳平衡是3個界面碳平衡的和,CBGAE=CBIA+CBIB+CBIC,可以據(jù)此有針對性地加強某一界面的碳管理。仍以祁連山馬鹿牧場和塔斯瑪尼亞奶牛牧場為例(表2)。

    表2 祁連山牧場和塔斯馬尼亞牧場3個界面的碳平衡

    兩個牧場系統(tǒng)的CBGAE>0,均為碳匯。

    草業(yè)系統(tǒng)的碳排放主要發(fā)生在前植物生產(chǎn)層或界面C,碳匯主要在植物生產(chǎn)層或界面A,在兩個牧場中均超過90%(表1和表2)。因此,碳匯管理應(yīng)減少前植物和后生物生產(chǎn)層的碳排放,增加植物和動物生產(chǎn)層的碳匯;或是增加界面A和B的碳匯,減少界面C的碳排放。

    3 草地農(nóng)業(yè)生態(tài)系統(tǒng)碳平衡分析的輸入/輸出法

    草業(yè)系統(tǒng)某一生產(chǎn)層或界面的碳平衡可以根據(jù)四個部分確定。輸入(Carbon input,CI),從草業(yè)系統(tǒng)外部輸入的化肥、機械、種子、畜力、廄肥等物質(zhì),光合作用吸收、同化CO2,微生物活動固定的CH4和N,塵降和水土保持積累碳氮。排放(Carbon emission,CE),人類、家畜消費草畜產(chǎn)品,包括食物、能源、畜力、日用品、廄肥等,在一定時間內(nèi)還原為GHG。固定部分(Carbon fixation,CF),牧草、家畜、排泄物等,在一定時間內(nèi)以可貯存的形態(tài)存在于草業(yè)系統(tǒng)內(nèi)。輸出(Carbon output,CO),輸出的草畜產(chǎn)品、種子、畜力、廄肥等。

    草業(yè)系統(tǒng)的碳平衡為CBGAE=CI+CF-CE-CO。如果CBGAE>0,即CI+CF>CE+CO,草業(yè)系統(tǒng)為碳匯,反之為碳源。某一生產(chǎn)層或某一界面的碳平衡同理計算,它們的初始碳量與碳平衡之和為當前碳量。

    表3 動物生產(chǎn)層碳平衡分析的輸入/輸出法舉例

    祁連山和塔斯瑪尼亞兩個草業(yè)系統(tǒng)的動物生產(chǎn)層均為碳源(表3),與生產(chǎn)層法和界面法的分析結(jié)果相同(表1,表2)。

    4 討論

    4.1 3種碳平衡計算方法的比較

    3種草業(yè)系統(tǒng)碳平衡分析方法各有側(cè)重(表4):4個生產(chǎn)層的方法關(guān)注能量與物質(zhì)流動規(guī)律,對應(yīng)著相應(yīng)的生產(chǎn)環(huán)節(jié);界面的方法強調(diào)關(guān)鍵生態(tài)過程,簡明地體現(xiàn)了碳平衡的作用機制;碳輸入/輸出法側(cè)重碳相對于生態(tài)系統(tǒng)的運動,忽略碳產(chǎn)生的環(huán)節(jié)和機制。目前碳平衡分析常用方法還有生活史分析(Life cycle analysis)[30,44-46]和元分析(Meta-analysis)[11,47-49];前者對碳流程自上而下梳理,與生產(chǎn)層法和界面法有相似之處,但是它忽略生產(chǎn)環(huán)節(jié)或系統(tǒng)組分之間的聯(lián)系和相互作用;后者對已有文獻的實測指標再一次的統(tǒng)計分析,相對定量地、概括地反映大時空尺度的碳動態(tài),但是對碳平衡的機理難以作具體的、有針對性的明晰。

    表4 草地農(nóng)業(yè)生態(tài)系統(tǒng)碳匯計算方法的比較

    草業(yè)系統(tǒng)碳匯,既有光合作用、碳塵降等自然過程,也有施廄肥等人工碳。當前的碳平衡分析方法難以區(qū)別碳的這兩種來源,也難以定量輸出碳的來源,不利于提高碳的管理效率,是今后研究需要改進的方面。

    4.2 影響碳平衡分析的因素

    4.2.1 大氣碳氮塵降

    每年有36億t碳塵降,近10億t不知去向[50],90%返還到土壤或水體,對濕潤地區(qū)生態(tài)系統(tǒng)碳平衡有重要影響。文中案例分析未考慮碳氮塵降,其生態(tài)效應(yīng)卻是當前研究熱點之一,短期內(nèi)氮塵降提高植物生產(chǎn)力,增加土壤微生物對有機質(zhì)的分解[51],直接影響植物碳庫和土壤碳庫[51];長期則降低土壤微生物量[52]。

    4.2.2 草地碳氮淋溶

    草業(yè)系統(tǒng)的可溶性碳氮會淋溶到地下水系統(tǒng)[53-54]。碳氮淋溶的影響因素主要是降水、灌溉和土地利用方式[38,53,55-57]。天然草地因主要位于干旱和極端干旱地區(qū),只存在微弱的淋溶過程[58];因而保持大量無機碳,其碳量超過土壤有機碳[2]。草地灌溉量越多,土壤無機碳的淋溶量越大[58]。放牧會增強土壤的淋溶作用,一方面因為放牧降低植被蓋度,加快土表枯落物分解,減少了降水截留[55];另一方面放牧促進植物地上部分產(chǎn)物向根系分配運輸,加快土壤有機質(zhì)分解,促進分解產(chǎn)物向下部轉(zhuǎn)移和聚集[59]。此外,土壤中空氣CO2分壓和土壤溶液的pH值也控制土壤無機碳濃度和通量[38],影響草業(yè)系統(tǒng)的碳平衡[60]。

    4.2.3 家畜排泄物的GHG排放

    集約化草業(yè)系統(tǒng)中家畜排泄物的GHG排放不容忽視。IPCC的GHG排放目錄中,家畜糞便的CH4排放是重要指標[35]。放牧家畜的排泄物多返還草地,少部分歸牧后存留于畜圈[10],這部分碳發(fā)酵后施肥到耕地中;排泄物的GHG排放多發(fā)生于存儲、堆漚過程,進入土壤的碳還影響土壤微生物呼吸和碳儲量。但是,放牧家畜或者廄肥有多少排泄物碳返還到土壤、多少成為土壤碳,尚不清楚。

    4.2.4 后生物生產(chǎn)層的碳足跡

    國際碳足跡核算多用在能源行業(yè)[61],尚沒有草業(yè)系統(tǒng)通用的碳足跡計算標準,而且后生物生產(chǎn)層的碳足跡仍是研究空白,用生命周期評價法(Life cycle assessment, LCA)可對其各個組分梳理,明確其耗能排放數(shù)據(jù)[62]。我國天然草地面積大,栽培草地方興未艾,草畜產(chǎn)品加工是草業(yè)發(fā)展的關(guān)鍵,卻與增碳減排相矛盾,草業(yè)系統(tǒng)發(fā)展后生物生產(chǎn)層以提高系統(tǒng)整體碳排放的經(jīng)濟效率是兼顧生態(tài)與生產(chǎn)的重要途徑,實施清潔發(fā)展機制(Clean Development Mechanism, CDM)對于草地農(nóng)業(yè)生態(tài)系統(tǒng)后生物生產(chǎn)層的健康發(fā)展有借鑒意義[63-64]。

    [1] 蔣延玲, 周廣勝. 興安落葉松林碳平衡和全球變化影響研究. 應(yīng)用生態(tài)學報, 2001, 12(4): 481-484.

    [2] 鐘華平, 樊江文, 于貴瑞, 韓彬, 胡中民, 岳燕珍, 梁飚. 草地生態(tài)系統(tǒng)碳循環(huán)研究進展. 草地學報, 2005, 13(Z1): 67-73.

    [3] 方精云, 陳安平. 中國森林植被碳庫的動態(tài)變化及其意義. 植物學報, 2001, 43(9): 967-973.

    [4] Hou F J. Adaptation of mixed crop-livestock systems in Asia//Fuhrer J, Gregory P, eds. Climate Change Impact and Adaptation in Agricultural Systems: Soil Ecosystem Management in Sustainable Agriculture. UK: CAB International, 2014: 155-166.

    [5] Yang Y H, Fang J Y, Smith P, Tang Y N, Chen A P, Ji C J, Hu H F, Rao S, Tan K, He J S. Changes in topsoil carbon stock in the Tibetan grasslands between the 1980s and 2004. Global Change Biology, 2009, 15(11): 2723-2729.

    [6] Yang Y H, Fang J Y, Ma W H, Smith P, Mohammat A, Wang S P, Wang W. Soil carbon stock and its changes in northern China′s grasslands from 1980s to 2000s. Global Change Biology, 2010, 16(11): 3036-3047.

    [7] 方精云, 楊元合, 馬文紅, 安尼瓦爾·買買提, 沈?;? 中國草地生態(tài)系統(tǒng)碳庫及其變化. 中國科學: 生命科學, 2010, 40(7): 566-576.

    [8] 王春梅, 楊林杰, 常生華, 侯扶江. 國內(nèi)外主要畜產(chǎn)品與飼料價格分析. 草業(yè)學報, 2014, 23(1): 300-311.

    [9] 張重麗, 于應(yīng)文, 付強國, 孫義, 侯扶江. 我國畜產(chǎn)品生產(chǎn)與貿(mào)易. 草業(yè)學報, 2015, 24(7): 163-171.

    [10] 侯扶江, 楊中藝. 放牧對草地的作用. 生態(tài)學報, 2006, 26(1): 244-264.

    [11] McSherry M E, Ritchie M E. Effects of grazing on grassland soil carbon: a global review. Global Change Biology, 2013, 19(5): 1347-1357.

    [12] Frank A B, Liebig M A, Hanson J D. Soil carbon dioxide fluxes in northern semiarid grasslands. Soil Biology and Biochemistry, 2002, 34(9): 1235-1241.

    [13] Hou F J, Nan Z B, Xie Y Z, Li X L, Lin H L, Ren J Z. Integrated crop-livestock production systems in China. The Rangeland Journal, 2008, 30(2): 221-231.

    [14] Casey J W, Holden N M. Quantification of GHG emissions from sucker-beef production in Ireland. Agricultural Systems, 2006, 90(1/3): 79-98.

    [15] 侯扶江, 南志標, 任繼周. 作物-家畜綜合生產(chǎn)系統(tǒng). 草業(yè)學報, 2009, 18(5): 211-234.

    [16] Chen J B, Hou F J, Chen X J, Wan X L, Millner J. Stocking rate and grazing season modify soil respiration on the Loess Plateau, China. Rangeland Ecology & Management. 2015, 68(1): 48-53.

    [17] 崔驍勇, 陳佐忠, 陳四清. 草地土壤呼吸研究進展. 生態(tài)學報, 2001, 21(2): 315-325.

    [18] Schlentner R E, Van Cleve K. Relationships between CO2 evolution from soil, substrate temperature, and substrate moisture in four mature forest types in interior Alaska. Canadian Journal of Forest Research, 1985, 15(1): 97-106.

    [19] Panikov N S, Paleeva M V, Dedysh S N, Dorofeev A G. Kinetic methods of determining the biomass and activity of different groups of soil microorganisms. Pochvovedenie, 1991, (8): 109-120.

    [20] 趙一廣, 刁其玉, 鄧凱東, 劉潔, 姜成鋼, 屠焰. 反芻動物甲烷排放的測定及調(diào)控技術(shù)研究進展. 動物營養(yǎng)學報, 2011, 23(5): 726-734.

    [21] 應(yīng)洪倉, 黃丹丹, 汪開英. 畜牧業(yè)溫室氣體檢測方法與技術(shù). 中國畜牧雜志, 2011, 47(10): 56-59, 63-63.

    [22] 王成杰, 汪詩平, 周禾. 放牧家畜甲烷氣體排放量測定方法研究進展. 草業(yè)學報, 2006, 15(1): 113-116.

    [23] Canadell P, Raupach M. Global carbon report: emissions will hit new heights in 2014. ECOS, 2014.

    [24] 任繼周, 李向林, 侯扶江. 草地農(nóng)業(yè)生態(tài)學研究進展與趨勢. 應(yīng)用生態(tài)學報, 2002, 13(8): 1017-1021.

    [25] Fan J W, Zhong H P, Harris W, Yu G R, Wang S Q, Hu Z M, Yue Y Z. Carbon storage in the grasslands of China based on field measurements of above-and below-ground biomass. Climatic Change, 2008, 86(3/4): 375-396.

    [26] Fang J Y, Guo Z D, Piao S L, Chen A P. Terrestrial vegetation carbon sinks in China, 1981—2000. Science in China Series D: Earth Sciences, 2007, 50(9): 1341-1350.

    [27] Ni J. Carbon storage in grasslands of China. Journal of Arid Environments, 2002, 50(2): 205-218.

    [28] 任繼周, 林慧龍. 草地土壤有機碳儲量模擬技術(shù)研究. 草業(yè)學報, 2013, 22(6): 280-294.

    [29] 趙雪. 草地旅游在草地生態(tài)系統(tǒng)中的作用及其持續(xù)發(fā)展. 中國草地, 2000, (5): 68-73.

    [30] 唐承財, 鐘林生, 成升魁. 旅游業(yè)碳排放研究進展. 地理科學進展, 2012, 31(4): 451-460.

    [31] 侯扶江, 安玉峰. 祁連山高寒牧區(qū)甘肅馬鹿產(chǎn)茸量的分析. 中國農(nóng)業(yè)科學, 2002, 35(10): 1269-1274.

    [32] 侯扶江, 李廣, 楊逢剛. 甘肅馬鹿夏冬季在祁連山高山草地的放牧行為. 生態(tài)學報. 2003, 23(9): 1807-1815.

    [33] 任繼周, 林惠龍, 侯向陽. 發(fā)展草地農(nóng)業(yè) 確保中國食物安全. 中國農(nóng)業(yè)科學, 2007, 40(3): 614-621.

    [34] 向志民, 何敏. 廄肥資源質(zhì)量分析與評價. 陜西農(nóng)業(yè)科學, 2000, (5): 14-16.

    [35] IPCC. Climate change 2013: the physical science basis. Cambridge: Cambridge University Press. 2013-09-30[2014-07-17]. http://www.ipcc.ch/report/ar5/wg1/#.Uq-tD7KBRR1.

    [36] 錢文熙, 敖維平, 玉蘇普. 塔里木馬鹿采食量與消化率研究. 中國草食動物科學, 2014, 34(2): 31-32, 38-38.

    [37] 任繼周, 朱興運, 王欽, 賈篤敬, 徐長林, 汪璽, 戈棠. 高山草原草地有效生產(chǎn)能力(P4)以前諸轉(zhuǎn)化階生產(chǎn)能力動態(tài)的研究. 生態(tài)學雜志, 1982, (2): 1-8.

    [38] Kindler R, Siemens J, Kaiser K, Walmsley D C, Bernhofer C, Buchmann N, Cellier P, Eugster W, Gleixner G, Grünwald T, Heim A, Ibrom A, Jones S K, Jones M, Klumpp K, Kutsch W, Larsen K S, Lehuger S, Loubet B, Mckenzie R, Moors E, Osborne B, Pilegaard K, Rebmann C, Saunders M, Schmidt M W I, Schrumpf M, Seyfferth J, Skiba U, Soussana J F, Sutton M A, Tefs C, Vowinckel B, Zeeman M J, Kaupenjohann M. Dissolved carbon leaching from soil is a crucial component of the net ecosystem carbon balance. Global Change Biology, 2011, 17(2): 1167-1185.

    [39] Parfitt R L, Percival H J, Dahlgren R A, Hill L F. Soil and solution chemistry under pasture and radiata pine in New Zealand. Plant and Soil, 1997, 191(2): 279-290.

    [40] Yuan H, Hou F J. Grazing intensity and soil depth effects on soil properties in alpine meadow pastures of Qilian Mountain in northwest China. Acta Agriculturae Scandinavica, 2015, 65(3): 222-232.

    [41] Zhang Y, Chen X J, Cheng Y X, Chang S H, Hou F J. Effects of stocking rates on functional group diversity and forage quality in rangeland of Qilian Mountain, China. Journal of Environmental Biology, 2015, 36(Special issue): 713-719.

    [42] Day G S. The product life cycle: analysis and applications issues. Journal of Marketing, 1981, 45(4): 60-67.

    [43] 侯扶江, 南志標, 肖金玉, 常生華. 重牧退化草地的植被、土壤及其耦合特征. 應(yīng)用生態(tài)學報, 2002, 13(8): 915-922.

    [44] Helin T, Sokka L, Soimakallio S, Pingoud K, Pajula T. Approaches for inclusion of forest carbon cycle in life cycle assessment: a review. GCB Bioenergy, 2013, 5(5): 475-486.

    [45] Perez-Garcia J, Lippke B, Comnick J, Manriquez C. An assessment of carbon pools, storage, and wood products market substitution using life-cycle analysis results. Wood and Fiber Science, 2005, 37(Special Issue): 140-148.

    [46] Weber C L, Clavin C. Life cycle carbon footprint of shale gas: Review of evidence and implications. Environmental Science & Technology, 2012, 46(11): 5688-5695.

    [47] Maillard é, Angers D A. Animal manure application and soil organic carbon stocks: A meta-analysis. Global Change Biology, 2014, 20(2): 666-679.

    [48] Bárcena T G, Kir L P, Vesterdal L, Stefánsdóttir H M, Gundersen P, Sigurdsson B D. Soil carbon stock change following afforestation in Northern Europe: a meta-analysis. Global Change Biology, 2014, 20(8): 2393-2405.

    [49] Mueller K E, Hobbie S E, Tilman D, Reich P B. Effects of plant diversity, N fertilization, and elevated carbon dioxide on grassland soil N cycling in a long-term experiment. Global Change Biology, 2013, 19(4): 1249-1261.

    [50] 劉再華. 大氣CO2兩個重要的匯. 科學通報, 2000, 45(21): 2348-2351.

    [51] 肖勝生, 董云社, 齊玉春, 彭琴, 何亞婷, 楊智杰. 草地生態(tài)系統(tǒng)土壤有機碳庫對人為干擾和全球變化的響應(yīng)研究進展. 地球科學進展, 2009, 24(10): 1138-1148.

    [52] Gloser J, Tesarova M. Litter, soil, and root respiration measurement: An improved compartmental analysis method. Pedobiologia, 1978, 18: 76-81.

    [53] Ciais P, Borges A V, Abril G, Meybeck M, Folberth G, Hauglustaine D, Janssens I A. The impact of lateral carbon fluxes on the European carbon balance. Biogeosciences, 2008, 5(5): 1259-1271.

    [54] Janssens I A, Freibauer A, Ciais P, Smith P, Nabuurs G J, Folberth G, Schlamadinger B, Hutjes R W A, Ceulemans R, Schulze E D, Valentini R, Johannes Dolman A. Europe′s terrestrial biosphere absorbs 7 to 12% of European anthropogenic CO2 emissions. Science, 2003, 300(5625): 1538-1542.

    [55] 陸晴, 王玉剛, 李彥, 唐立松. 干旱區(qū)不同土壤和作物灌溉量的無機碳淋溶特征實驗研究. 干旱區(qū)地理, 2013, 36(3): 450-456.

    [56] 鄒誠, 徐福利, 閆亞丹. 黃土高原丘陵溝壑區(qū)不同土地利用模式對土壤氮素淋溶的影響. 水土保持研究, 2009, 16(3): 114-116.

    [57] Schulze E D, Luyssaert S, Ciais P, Freibauer A, Janssens I A. Importance of methane and nitrous oxide for Europe′s terrestrial greenhouse-gas balance. Nature Geoscience, 2009, 2(12): 842-850.

    [58] 樊恒文, 賈曉紅, 張景光, 馬鳳云, 李新榮. 干旱區(qū)土地退化與荒漠化對土壤碳循環(huán)的影響. 中國沙漠, 2002, 22(6): 525-533.

    [59] 耿元波, 羅光強, 袁國富, 李明峰, 孟維奇, 董云社. 農(nóng)墾及放牧對溫帶半干旱草原土壤碳素的影響. 農(nóng)業(yè)環(huán)境科學學報, 2008, 27(6): 2518-2523.

    [60] 曹裕松, 李志安, 江遠清, 丁思統(tǒng). 陸地生態(tài)系統(tǒng)土壤呼吸研究進展. 江西農(nóng)業(yè)大學學報, 2004, 26(1): 138-143.

    [61] 劉亮. 基于LCA模型的工程項目碳足跡評估實證研究. 贛州: 江西理工大學, 2011.

    [62] Huang Y, Sun W J, Zhang W, Yu Y Q. Changes in soil organic carbon of terrestrial ecosystems in China: A mini-review. Science China Life Sciences, 2010, 53(7): 766-775.

    [63] 胡秀蓮, 李俊峰. 關(guān)于建立我國清潔發(fā)展機制(CDM)項目運行管理機制的幾點建議(一). 中國能源, 2001, (8): 6-9.

    [64] 孫慧麗. 清潔發(fā)展機制下我國碳排放市場研究. 南京: 南京林業(yè)大學, 2012.

    Carbon balance analysis methods of grassland agro-ecosystems

    LOU Shanning, CHEN Xianjiang, HOU Fujiang*

    StateKeyLaboratoryofGrasslandAgro-ecosystems,CollegeofPastoralAgricultureScienceandTechnology,LanzhouUniversity,Lanzhou730020,China

    According to the structures of grassland agro-ecosystems, their carbon balance represent the sum of the carbon balances of the four production levels, or the sum from its three interfaces. The carbon balance per production level or interface is the sum of the carbon that is taken up and lost through the different processes in that level or interface. The carbon balance analysis method based on four production levels can quantify which of the processes in the production are either a carbon sink or source. For example, carbon and nitrogen accumulation such as Carbon and Nitrogen sink, soil and water conservation are carbon sink process in Pre-plant production level. The grassland tourism hunting and so on, are the process of carbon source in Pre-plant production level. Photosynthesis and carbon and nitrogen assimilation of biological nitrogen fixation are the process of carbon sink in Plant production level. Greenhouse gases emission is the process of carbon source in Plant production level. The process of carbon sink and source in Animal production level are mostly Carbon fixation in animal product and greenhouse gases emission of ruminant. Carbon fixation in forage and animal products are the process of carbon sink in Post-biological production level. Processing and transport activities are the process of carbon source in Post-biological production level. This can subsequently facilitate improvements in the management of carbon sequestration in grassland production. Furthermore, the carbon balance analysis method based on three interfaces can determine the production mechanisms of carbon sinks and sources, as well as their spatial and quantitative relationships. This can contribute to regulating and controlling the carbon emissions from grassland agro-ecosystems. However, both methods cannot accurately quantify a grassland system′s paths of carbon uptake and output, or its utilization efficiency. The carbon balance analysis by the balance between carbon input and output can quantify the pathway of carbon movement. This method is simple and more succinct; however, it is not applicable to farmland scale management of carbon sequestration. Considering a ranch in the Qilian Mountains, China, and a dairy farm in Tasmania, Australia, as examples, three methods were used for analyzing their carbon balances. The results showed that carbon emissions from tourism, product processing, and marketing accounted for the major portions of the total emissions from the grazing management system. Their main carbon sinks were carbon stored in the forage and the soil. Overall, optimal grassland management would benefit from adding carbon sinks and reducing carbon emissions.

    grassland agro-ecosystems; production level; interface; carbon balance; grazing; greenhouse gases

    國家重點基礎(chǔ)研究發(fā)展計劃資助項目(2014CB138706);國家自然科學基金資助項目(31172249);長江學者和創(chuàng)新團隊發(fā)展計劃資助項目(IRT13019);國家科技支撐計劃課題資助項目(2011BAD17B0203)

    2015-07-17;

    日期:2016-06-13

    10.5846/stxb201507171503

    * 通訊作者Corresponding author.E-mail: cyhoufj@lzu.edu.cn

    婁珊寧,陳先江,侯扶江.草地農(nóng)業(yè)生態(tài)系統(tǒng)的碳平衡分析方法.生態(tài)學報,2017,37(2):557-565.

    Lou S N, Chen X J, Hou F J.Carbon balance analysis methods of grassland agro-ecosystems.Acta Ecologica Sinica,2017,37(2):557-565.

    猜你喜歡
    草業(yè)家畜牧場
    《草食家畜》征稿簡則
    草食家畜(2024年1期)2024-02-06 09:59:00
    青海草業(yè)
    冷瀉冷治 熱瀉熱治 中醫(yī)治療家畜腹瀉有妙招
    中醫(yī)治療家畜腹瀉有妙招
    海上牧場
    海峽姐妹(2020年6期)2020-07-25 01:26:28
    2017 年第1 期《草業(yè)科學》審稿專家
    叮當牧場
    2015年3月草業(yè)科學大事記
    2015年6月草業(yè)科學大事記
    Gift Horse
    99精品在免费线老司机午夜| 叶爱在线成人免费视频播放| 啦啦啦免费观看视频1| 国产乱人伦免费视频| 亚洲欧洲精品一区二区精品久久久| 欧美色视频一区免费| 成年女人永久免费观看视频| 人人妻人人澡欧美一区二区| 少妇熟女aⅴ在线视频| xxx96com| av女优亚洲男人天堂 | 亚洲七黄色美女视频| 国产精品av视频在线免费观看| 亚洲精品在线美女| 色综合婷婷激情| 国产高清videossex| 久久亚洲真实| 最近视频中文字幕2019在线8| 亚洲国产欧美人成| 男女下面进入的视频免费午夜| 女同久久另类99精品国产91| 最新美女视频免费是黄的| 午夜精品在线福利| 精品久久久久久久毛片微露脸| 最好的美女福利视频网| 无限看片的www在线观看| 人妻久久中文字幕网| 欧洲精品卡2卡3卡4卡5卡区| 亚洲九九香蕉| 亚洲精品在线观看二区| 精品福利观看| 国产欧美日韩精品一区二区| 国产日本99.免费观看| 91av网站免费观看| 国产精品免费一区二区三区在线| 精品电影一区二区在线| 亚洲专区字幕在线| 成人亚洲精品av一区二区| 1024手机看黄色片| 久久中文字幕一级| 日本 av在线| 身体一侧抽搐| 国产成+人综合+亚洲专区| 法律面前人人平等表现在哪些方面| 成人无遮挡网站| 亚洲av电影在线进入| 精华霜和精华液先用哪个| 国产一区二区三区视频了| 免费无遮挡裸体视频| 成人一区二区视频在线观看| 国产单亲对白刺激| 桃红色精品国产亚洲av| 偷拍熟女少妇极品色| 国产精品久久电影中文字幕| 久久午夜亚洲精品久久| 最近最新中文字幕大全免费视频| 草草在线视频免费看| 亚洲国产高清在线一区二区三| 欧美高清成人免费视频www| 久久久久国产一级毛片高清牌| 国模一区二区三区四区视频 | 亚洲国产色片| 亚洲欧美一区二区三区黑人| 久久这里只有精品中国| 国产99白浆流出| 一个人看视频在线观看www免费 | 禁无遮挡网站| 国内毛片毛片毛片毛片毛片| 久久香蕉国产精品| 18禁国产床啪视频网站| 午夜a级毛片| 亚洲美女黄片视频| 免费在线观看影片大全网站| 超碰成人久久| 一级毛片精品| 久久午夜综合久久蜜桃| 少妇的丰满在线观看| 黄色片一级片一级黄色片| tocl精华| 高清毛片免费观看视频网站| 欧美日韩乱码在线| 听说在线观看完整版免费高清| 亚洲熟妇熟女久久| 他把我摸到了高潮在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产av在哪里看| 999久久久精品免费观看国产| 色老头精品视频在线观看| netflix在线观看网站| 99国产精品99久久久久| 国产综合懂色| 老司机福利观看| 亚洲成a人片在线一区二区| 欧美日韩一级在线毛片| 亚洲精品久久国产高清桃花| 成人午夜高清在线视频| 亚洲aⅴ乱码一区二区在线播放| 宅男免费午夜| 老汉色∧v一级毛片| 亚洲专区国产一区二区| 欧美一级a爱片免费观看看| 精品久久蜜臀av无| 午夜激情福利司机影院| 亚洲黑人精品在线| 国产一区二区激情短视频| 18禁黄网站禁片免费观看直播| 老汉色∧v一级毛片| 熟女电影av网| 黄色成人免费大全| 亚洲最大成人中文| 不卡av一区二区三区| 香蕉久久夜色| 亚洲午夜理论影院| 亚洲av中文字字幕乱码综合| 怎么达到女性高潮| 久久久久国产一级毛片高清牌| 精品国产乱子伦一区二区三区| 欧美一级a爱片免费观看看| 一卡2卡三卡四卡精品乱码亚洲| 在线看三级毛片| 热99re8久久精品国产| 18禁裸乳无遮挡免费网站照片| 91久久精品国产一区二区成人 | 日本黄色视频三级网站网址| 无限看片的www在线观看| 最近最新免费中文字幕在线| 夜夜躁狠狠躁天天躁| 99国产精品一区二区蜜桃av| 午夜福利18| 一进一出抽搐gif免费好疼| 国内少妇人妻偷人精品xxx网站 | 国产精品久久久久久人妻精品电影| 日本免费a在线| 亚洲欧美日韩卡通动漫| 久久中文字幕一级| 岛国在线免费视频观看| 国产精品 国内视频| 中亚洲国语对白在线视频| 这个男人来自地球电影免费观看| 五月玫瑰六月丁香| 久久久成人免费电影| 少妇丰满av| 精品国产三级普通话版| 久久久成人免费电影| 精品乱码久久久久久99久播| 91麻豆精品激情在线观看国产| 欧美黄色片欧美黄色片| 亚洲精品美女久久久久99蜜臀| 亚洲熟妇熟女久久| 叶爱在线成人免费视频播放| 男女之事视频高清在线观看| 免费人成视频x8x8入口观看| 久久久久国内视频| 久久九九热精品免费| 1024手机看黄色片| 日本撒尿小便嘘嘘汇集6| 精品日产1卡2卡| 国产激情久久老熟女| 精品久久久久久久久久久久久| 观看免费一级毛片| 黑人欧美特级aaaaaa片| 国产成人av激情在线播放| 日韩大尺度精品在线看网址| 一级毛片精品| 日韩欧美在线乱码| 国产精品99久久久久久久久| 不卡av一区二区三区| 老熟妇乱子伦视频在线观看| 国产毛片a区久久久久| 国产免费男女视频| 亚洲专区字幕在线| 国语自产精品视频在线第100页| 婷婷精品国产亚洲av| 又黄又爽又免费观看的视频| 成人高潮视频无遮挡免费网站| 老司机在亚洲福利影院| 午夜福利在线观看免费完整高清在 | 一a级毛片在线观看| 久久久国产欧美日韩av| 热99在线观看视频| 色综合亚洲欧美另类图片| 国产三级黄色录像| 最近最新中文字幕大全免费视频| 啪啪无遮挡十八禁网站| 亚洲欧美日韩高清专用| 国产日本99.免费观看| 99久国产av精品| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人精品一区久久| 国产欧美日韩精品一区二区| 一本综合久久免费| 久久久久国产精品人妻aⅴ院| 亚洲 欧美 日韩 在线 免费| 18禁观看日本| 美女大奶头视频| 国产蜜桃级精品一区二区三区| 国产精品,欧美在线| 岛国视频午夜一区免费看| 嫩草影院精品99| 51午夜福利影视在线观看| 亚洲国产精品成人综合色| 亚洲国产色片| 女人被狂操c到高潮| 午夜福利成人在线免费观看| 久久国产精品人妻蜜桃| 中文字幕av在线有码专区| 亚洲国产精品合色在线| 麻豆成人av在线观看| 午夜两性在线视频| 色吧在线观看| 国产成人精品久久二区二区91| 国产精品久久久久久人妻精品电影| 国产亚洲av嫩草精品影院| 精品一区二区三区四区五区乱码| 一卡2卡三卡四卡精品乱码亚洲| 夜夜躁狠狠躁天天躁| 香蕉久久夜色| 国产亚洲精品久久久久久毛片| 日本与韩国留学比较| 久久精品国产清高在天天线| 嫩草影院精品99| 国产精品久久视频播放| 成年人黄色毛片网站| 成人鲁丝片一二三区免费| 国产在线精品亚洲第一网站| 少妇的逼水好多| av福利片在线观看| 精品电影一区二区在线| 亚洲国产看品久久| 欧美黄色片欧美黄色片| 午夜亚洲福利在线播放| 亚洲色图 男人天堂 中文字幕| 国产69精品久久久久777片 | 国产高清视频在线播放一区| 精品久久久久久,| 日本精品一区二区三区蜜桃| 91av网站免费观看| 国产亚洲精品久久久com| 国产成人欧美在线观看| 亚洲人成网站高清观看| 久久精品亚洲精品国产色婷小说| 九九在线视频观看精品| 国产亚洲精品久久久久久毛片| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线观看免费| 亚洲熟女毛片儿| 亚洲 国产 在线| 少妇丰满av| 深夜精品福利| 偷拍熟女少妇极品色| 欧美xxxx黑人xx丫x性爽| 国产亚洲av嫩草精品影院| 成年免费大片在线观看| 最近最新中文字幕大全电影3| 亚洲国产欧美网| 搞女人的毛片| 动漫黄色视频在线观看| 18禁黄网站禁片免费观看直播| 久久天堂一区二区三区四区| 精品久久久久久久毛片微露脸| 国产成人影院久久av| 精品熟女少妇八av免费久了| 国产精品 欧美亚洲| 韩国av一区二区三区四区| 亚洲熟女毛片儿| 五月伊人婷婷丁香| 亚洲精品国产精品久久久不卡| 91在线精品国自产拍蜜月 | 麻豆av在线久日| 精品熟女少妇八av免费久了| 成人三级做爰电影| 中文字幕av在线有码专区| 免费无遮挡裸体视频| 欧美又色又爽又黄视频| 热99在线观看视频| 中文字幕熟女人妻在线| 午夜福利在线观看吧| 久久这里只有精品19| 最近最新中文字幕大全电影3| 日韩国内少妇激情av| 精品国产亚洲在线| 午夜福利欧美成人| 久久精品综合一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产97色在线日韩免费| 亚洲欧美日韩高清在线视频| 日韩精品中文字幕看吧| 亚洲国产高清在线一区二区三| 中文字幕人成人乱码亚洲影| 欧美一区二区国产精品久久精品| 又大又爽又粗| 亚洲成a人片在线一区二区| 久久人妻av系列| 黄色日韩在线| 精品久久久久久,| 国产一区二区激情短视频| 免费搜索国产男女视频| 日本撒尿小便嘘嘘汇集6| 久久久久久久久免费视频了| 好看av亚洲va欧美ⅴa在| 日本黄色视频三级网站网址| 热99在线观看视频| 国产黄色小视频在线观看| 国产麻豆成人av免费视频| 黄色丝袜av网址大全| 特级一级黄色大片| 97碰自拍视频| 久久精品影院6| 国产精品久久视频播放| 精品日产1卡2卡| 国产精品99久久久久久久久| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 狠狠狠狠99中文字幕| 操出白浆在线播放| 精品久久久久久久久久久久久| 欧美不卡视频在线免费观看| 波多野结衣高清无吗| 97超视频在线观看视频| 99在线视频只有这里精品首页| 亚洲国产欧美网| 国产视频一区二区在线看| 午夜福利欧美成人| 久久久久久大精品| 给我免费播放毛片高清在线观看| 久久伊人香网站| 白带黄色成豆腐渣| 在线观看一区二区三区| 一个人免费在线观看的高清视频| 精品一区二区三区视频在线 | 美女高潮的动态| 韩国av一区二区三区四区| 婷婷精品国产亚洲av在线| 床上黄色一级片| tocl精华| 国产成人av激情在线播放| 在线a可以看的网站| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 桃色一区二区三区在线观看| 9191精品国产免费久久| www.自偷自拍.com| 国产av不卡久久| 美女cb高潮喷水在线观看 | 精品人妻1区二区| 亚洲精品美女久久久久99蜜臀| 又大又爽又粗| 最近最新中文字幕大全电影3| 欧美日韩黄片免| 在线免费观看的www视频| 黄片大片在线免费观看| 久久午夜综合久久蜜桃| av天堂中文字幕网| 亚洲国产看品久久| 亚洲精品中文字幕一二三四区| 又爽又黄无遮挡网站| 99久国产av精品| 日本熟妇午夜| av国产免费在线观看| 岛国在线免费视频观看| 精品久久久久久成人av| 欧美激情久久久久久爽电影| 美女大奶头视频| 18禁美女被吸乳视频| 欧美中文日本在线观看视频| 91在线观看av| 日本五十路高清| 91av网站免费观看| 老司机午夜福利在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 美女大奶头视频| 制服丝袜大香蕉在线| www国产在线视频色| 精品一区二区三区av网在线观看| 91麻豆av在线| 国产欧美日韩一区二区精品| 91字幕亚洲| 精品日产1卡2卡| 日本精品一区二区三区蜜桃| 国产精品精品国产色婷婷| 成人鲁丝片一二三区免费| 亚洲午夜理论影院| av中文乱码字幕在线| 精品国产乱码久久久久久男人| 少妇的丰满在线观看| 国产精品99久久99久久久不卡| 国内揄拍国产精品人妻在线| 桃红色精品国产亚洲av| 窝窝影院91人妻| 亚洲精品一区av在线观看| 亚洲七黄色美女视频| 国产人伦9x9x在线观看| 麻豆久久精品国产亚洲av| 亚洲av日韩精品久久久久久密| av片东京热男人的天堂| 91九色精品人成在线观看| 少妇丰满av| 男人的好看免费观看在线视频| 亚洲欧美日韩高清在线视频| 午夜福利在线在线| 国产成人精品久久二区二区91| 国产成人啪精品午夜网站| 久久精品91蜜桃| 欧美性猛交黑人性爽| 欧美成人免费av一区二区三区| 熟女电影av网| 久久久久久国产a免费观看| 欧美中文日本在线观看视频| 一二三四在线观看免费中文在| 岛国视频午夜一区免费看| 很黄的视频免费| 曰老女人黄片| 亚洲精品美女久久久久99蜜臀| xxxwww97欧美| 国产精品99久久99久久久不卡| 十八禁网站免费在线| 两人在一起打扑克的视频| 九九在线视频观看精品| 九色成人免费人妻av| 桃色一区二区三区在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品日韩av在线免费观看| 久久草成人影院| 亚洲精品在线美女| 精品国产乱码久久久久久男人| 最近最新中文字幕大全电影3| 国产男靠女视频免费网站| 嫩草影院入口| 麻豆一二三区av精品| 丁香欧美五月| 美女cb高潮喷水在线观看 | 日韩三级视频一区二区三区| 91麻豆av在线| 国产99白浆流出| 亚洲av电影在线进入| 国产成人啪精品午夜网站| 国产爱豆传媒在线观看| 精品免费久久久久久久清纯| 国产aⅴ精品一区二区三区波| 此物有八面人人有两片| 欧美在线一区亚洲| 丁香六月欧美| 午夜a级毛片| 美女 人体艺术 gogo| 手机成人av网站| 久久久久久国产a免费观看| 欧美中文综合在线视频| 日韩人妻高清精品专区| a级毛片在线看网站| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 久久久国产成人免费| 免费在线观看视频国产中文字幕亚洲| 美女高潮喷水抽搐中文字幕| 国产精品美女特级片免费视频播放器 | 国产精品爽爽va在线观看网站| 亚洲av电影不卡..在线观看| 日本黄大片高清| 国产精品精品国产色婷婷| 久久午夜综合久久蜜桃| 日本五十路高清| 国产成人影院久久av| 亚洲av片天天在线观看| 热99在线观看视频| 在线观看免费午夜福利视频| 手机成人av网站| 亚洲国产欧美人成| 精品久久久久久成人av| 久久久久久人人人人人| 国产高清有码在线观看视频| 中文字幕久久专区| 18禁黄网站禁片免费观看直播| 精品欧美国产一区二区三| 国产伦在线观看视频一区| 日韩欧美一区二区三区在线观看| 在线免费观看不下载黄p国产 | 手机成人av网站| 亚洲精品美女久久久久99蜜臀| 亚洲成a人片在线一区二区| tocl精华| 久久国产精品影院| 久久亚洲真实| 亚洲最大成人中文| 免费观看人在逋| 亚洲专区字幕在线| 国产精品女同一区二区软件 | 日韩免费av在线播放| 精品一区二区三区视频在线 | 啦啦啦免费观看视频1| 久久久久国产一级毛片高清牌| 午夜福利18| 亚洲一区高清亚洲精品| 免费看a级黄色片| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| 国产成人av教育| 91字幕亚洲| 黄色片一级片一级黄色片| 999久久久国产精品视频| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 精品欧美国产一区二区三| 国产高潮美女av| 亚洲人与动物交配视频| 日韩免费av在线播放| av片东京热男人的天堂| 欧美午夜高清在线| netflix在线观看网站| 精品久久久久久久末码| 久久久久久九九精品二区国产| 中出人妻视频一区二区| 久久久久久久久久黄片| 韩国av一区二区三区四区| 国模一区二区三区四区视频 | 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线av高清观看| 欧美绝顶高潮抽搐喷水| 亚洲五月婷婷丁香| 无遮挡黄片免费观看| 九色国产91popny在线| 国产精品爽爽va在线观看网站| 一个人看的www免费观看视频| 操出白浆在线播放| 久久精品aⅴ一区二区三区四区| 日日摸夜夜添夜夜添小说| 欧美激情在线99| 午夜视频精品福利| 99热这里只有精品一区 | 真实男女啪啪啪动态图| 在线观看美女被高潮喷水网站 | 久久人妻av系列| 欧美性猛交╳xxx乱大交人| 国产精品亚洲一级av第二区| e午夜精品久久久久久久| 久久中文看片网| 一本综合久久免费| 香蕉国产在线看| 91av网一区二区| 午夜精品一区二区三区免费看| 18禁美女被吸乳视频| 我的老师免费观看完整版| 色av中文字幕| 国产精品日韩av在线免费观看| 欧美日韩精品网址| 日韩有码中文字幕| 九九久久精品国产亚洲av麻豆 | 五月玫瑰六月丁香| 在线十欧美十亚洲十日本专区| 亚洲色图 男人天堂 中文字幕| 一级a爱片免费观看的视频| 脱女人内裤的视频| 久久亚洲真实| 亚洲精品在线观看二区| 久久国产精品人妻蜜桃| 亚洲av第一区精品v没综合| 我的老师免费观看完整版| 一卡2卡三卡四卡精品乱码亚洲| 国产精品九九99| 国产精品98久久久久久宅男小说| 很黄的视频免费| 久久中文看片网| 免费在线观看视频国产中文字幕亚洲| 久9热在线精品视频| 午夜福利欧美成人| 中文字幕av在线有码专区| 这个男人来自地球电影免费观看| 亚洲成人精品中文字幕电影| 啦啦啦观看免费观看视频高清| 日韩av在线大香蕉| 日韩欧美国产一区二区入口| 最新美女视频免费是黄的| 男女之事视频高清在线观看| 老熟妇仑乱视频hdxx| 成人18禁在线播放| 国内毛片毛片毛片毛片毛片| 午夜福利在线观看免费完整高清在 | 亚洲性夜色夜夜综合| 成年女人看的毛片在线观看| 色综合欧美亚洲国产小说| 琪琪午夜伦伦电影理论片6080| 精品不卡国产一区二区三区| 欧美3d第一页| 国产黄a三级三级三级人| h日本视频在线播放| av片东京热男人的天堂| 亚洲熟妇中文字幕五十中出| 亚洲国产精品合色在线| 精品无人区乱码1区二区| 岛国视频午夜一区免费看| 国产毛片a区久久久久| 美女 人体艺术 gogo| 亚洲av美国av| 欧美不卡视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧洲综合997久久,| 网址你懂的国产日韩在线| e午夜精品久久久久久久| 久久精品影院6| 久久国产精品影院| 亚洲第一欧美日韩一区二区三区| 琪琪午夜伦伦电影理论片6080| 中文亚洲av片在线观看爽| 精品日产1卡2卡| 噜噜噜噜噜久久久久久91| 中文亚洲av片在线观看爽| 老司机福利观看| 嫩草影院精品99| 1000部很黄的大片| av黄色大香蕉| 91在线精品国自产拍蜜月 | 久久久精品欧美日韩精品| 国内精品美女久久久久久| 成人18禁在线播放| 一夜夜www| 国产精品久久久久久亚洲av鲁大|