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

    基于Biome-BGC模型的西雙版納橡膠林碳收支模擬

    2017-10-13 12:32:57孫燕瓷馬友鑫曹坤芳沈金祥張一平梅岑岑劉文俊
    生態(tài)學(xué)報(bào) 2017年17期
    關(guān)鍵詞:橡膠林碳循環(huán)西雙版納

    孫燕瓷,馬友鑫, 曹坤芳,沈金祥,張一平,梅岑岑,劉文俊

    1 綜合保護(hù)中心,西雙版納熱帶植物園, 勐侖 666303 2 云南經(jīng)濟(jì)管理學(xué)院, 昆明 650106 3 亞熱帶農(nóng)業(yè)資源保護(hù)與利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣西大學(xué)林學(xué)院, 南寧 530004 4 云南國(guó)土資源職業(yè)學(xué)院, 昆明 652501 5 深圳市深港產(chǎn)學(xué)研環(huán)保工程技術(shù)股份有限公司, 深圳 518057

    基于Biome-BGC模型的西雙版納橡膠林碳收支模擬

    孫燕瓷2,馬友鑫1,*, 曹坤芳3,沈金祥4,張一平1,梅岑岑5,劉文俊1

    1 綜合保護(hù)中心,西雙版納熱帶植物園, 勐侖 666303 2 云南經(jīng)濟(jì)管理學(xué)院, 昆明 650106 3 亞熱帶農(nóng)業(yè)資源保護(hù)與利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣西大學(xué)林學(xué)院, 南寧 530004 4 云南國(guó)土資源職業(yè)學(xué)院, 昆明 652501 5 深圳市深港產(chǎn)學(xué)研環(huán)保工程技術(shù)股份有限公司, 深圳 518057

    以西雙版納橡膠適宜種植區(qū)(海拔550—600m)的橡膠林 (Heveabrasiliensis)為研究對(duì)象,應(yīng)用參數(shù)同化后的Biome-BGC模型模擬了1959—2012年橡膠林的碳循環(huán)。結(jié)果表明,(1)與渦度相關(guān)監(jiān)測(cè)結(jié)果相比,橡膠林年總初級(jí)生產(chǎn)力(Gross Primary Productivity,GPP)、年總呼吸(Total Respiration, Rt)的模擬精度分別為98.37%和90%。由于對(duì)年GPP的過低估計(jì)和對(duì)年Rt的過高估計(jì),年凈生態(tài)系統(tǒng)交換量(Net Ecosystem Exchange, NEE)的模擬值比實(shí)測(cè)值低157.35 g C m-2a-1。但若考慮干膠碳(139 g C m-2a-1),模擬值與實(shí)測(cè)值十分接近;(2)橡膠林在模擬進(jìn)行的前8年里因異養(yǎng)呼吸較高,以碳排放為主, NEE平均約357 g C m-2a-1;之后轉(zhuǎn)為以碳固定為主,NEE平均約-146 g C m-2a-1;(3)橡膠林在40年的更新周期中可固定碳1835 g C m-2,是一個(gè)弱的碳匯。但與熱帶雨林相同周期固碳6720 g C m-2相比,仍為碳源。以上結(jié)果為深入了解橡膠種植對(duì)區(qū)域碳循環(huán)的影響提供了科學(xué)依據(jù),建議當(dāng)?shù)卣环矫嬉杏?jì)劃的對(duì)老膠林進(jìn)行更新,以維持當(dāng)前橡膠林生態(tài)系統(tǒng)中的碳平衡;另一方面要注重對(duì)熱帶雨林的保護(hù),從而實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)和生態(tài)環(huán)境保護(hù)的協(xié)調(diào)發(fā)展。

    橡膠林;凈生態(tài)系統(tǒng)交換量;Biome-BGC模型;西雙版納

    KeyWords: rubber plantation; net ecosystem exchange; Biome-BGC model; Xishuangbanna

    自1975年以來,大氣累積CO2濃度增加了40%,21世紀(jì)末期及以后全球氣候系統(tǒng)變暖與累積CO2排放密切相關(guān)[1]。在熱帶地區(qū),橡膠(Heveabrasiliensis)種植是影響區(qū)域碳循環(huán)的一個(gè)重要方面[2- 8],自《京都協(xié)議》簽訂生效以來,一個(gè)國(guó)家或地區(qū)的社會(huì)經(jīng)濟(jì)發(fā)展與碳排放緊密相連,因此關(guān)于橡膠種植過程中碳收支的動(dòng)態(tài)變化及其對(duì)區(qū)域碳循環(huán)的準(zhǔn)確評(píng)估,已成為全球變化研究的熱點(diǎn)之一。

    生理生態(tài)過程模型由于考慮了植物的生理生態(tài)特性,多將土壤—植被—大氣連續(xù)體(SPAC)作為一個(gè)系統(tǒng),從植物的生理過程反映生態(tài)系統(tǒng)與大氣之間的物質(zhì)和能量流動(dòng)[9],因此被廣泛用于生態(tài)系統(tǒng)碳循環(huán)的研究中。比較典型的生理生態(tài)過程模型有Forest-BGC模型[10]、Biome-BGC模型[11]、CENTURY模型[12]、BGC-ES模型[13]等。其中Biome-BGC模型因用C++語言編寫,易于編譯和理解,且源代碼開放,成為被廣泛使用的一種生理生態(tài)過程模型。該模型中設(shè)置了植物光合作用、呼吸作用、蒸騰作用及凋落物分解等主要生理生態(tài)過程的循環(huán)模式與計(jì)算方法[14],目前已被用于農(nóng)作物生態(tài)系統(tǒng)[15]、人為干擾森林生態(tài)系統(tǒng)[16- 17]、濕地紅松生態(tài)系統(tǒng)[18]、草地生態(tài)系統(tǒng)[19]等的碳循環(huán)的模擬。利用Biome-BGC模型對(duì)熱帶人工經(jīng)濟(jì)林生態(tài)系統(tǒng)碳循環(huán)的模擬尚不多見。

    西雙版納地處熱帶北緣,受經(jīng)濟(jì)利益和人口快速增長(zhǎng)的共同驅(qū)動(dòng),自20世紀(jì)50年代以來橡膠種植業(yè)(特別是民營(yíng)橡膠)迅猛發(fā)展,橡膠種植面積已超過了該區(qū)總面積的20%[20],導(dǎo)致該區(qū)約139567hm2的熱帶雨林[21]和364116hm2的熱帶自然林[22]消失,極大地改變了區(qū)域的碳收支格局[3]。橡膠林的最佳固碳周期是40年[23],Song[24]利用渦度監(jiān)測(cè)技術(shù)指出,西雙版納一片33年生的橡膠林表現(xiàn)為碳匯。但目前還沒有關(guān)于其它林齡(或發(fā)展階段)橡膠林的研究結(jié)果,而這方面的數(shù)據(jù)是當(dāng)前有效應(yīng)對(duì)全球氣候變化和國(guó)際氣候談判所迫切需要的。

    本研究以西雙版納海拔550—600m橡膠適宜種植區(qū)[25]的橡膠林為研究對(duì)象,利用Biome-BGC模型模擬1959—2012年橡膠林的碳循環(huán),著重分析40年中橡膠林碳收支的時(shí)間變化規(guī)律及凈生態(tài)系統(tǒng)交換量(Net Ecosystem Exchange, NEE),以期為研究其它海拔橡膠林的碳收支及綜合評(píng)價(jià)橡膠種植對(duì)區(qū)域碳循環(huán)的影響提供科學(xué)依據(jù)。

    1 研究區(qū)域與研究方法

    1.1 研究區(qū)概況

    表1為5種參數(shù)化方案以及集合平均模擬的24 h、6 h和3 h累積降水的評(píng)估結(jié)果。對(duì)于24 h累積降水,距平相關(guān)系數(shù)和均方根誤差的值均變優(yōu)。相對(duì)于最優(yōu)的單方案結(jié)果,集成后的距平相關(guān)系數(shù)提高了2個(gè)百分點(diǎn),均方根誤差降低了0.2 mm。從TS評(píng)分的結(jié)果可以看出,集成后對(duì)于中雨和暴雨的TS評(píng)分值改進(jìn)較為明顯,優(yōu)于最優(yōu)的單方案,而小雨和大雨的TS值不如最優(yōu)的單方案,但也僅次之。該結(jié)論與陳茂欽等(2012)對(duì)江淮和華南兩次暴雨過程進(jìn)行多微物理方案的集成預(yù)報(bào)試驗(yàn)的結(jié)論相一致,均得出多方案集成后結(jié)果較為優(yōu)異且穩(wěn)定的結(jié)論。

    西雙版納(21°08′—22°36′N,99°56′—101°50′E)位于云南省南端,與緬甸、老撾接壤,擁有土地總面積19120 km2。屬熱帶北緣季風(fēng)氣候,年降雨量為1162.1mm,其中約85%的降雨集中在雨季(5—10月);干季降水較少(11月至次年4月),但有霧,有霧日占70—80%,一般在22:30開始起霧,直到第二天11:00才逐漸消散。年均溫為22.3℃,最熱月出現(xiàn)在5—6月25.8℃,最冷月為1月16.5℃,年均日照時(shí)數(shù)為2152.9h。主要植被類型有熱帶季節(jié)雨林、熱帶山地雨林、熱帶季節(jié)性濕潤(rùn)林、熱帶季雨林和熱帶山地常綠闊葉林[26]。境內(nèi)地勢(shì)高差起伏較大,海拔475—2429 m,橡膠在各個(gè)海拔均有種植,目前已突破種植上限1200m[3,21],至2014年橡膠種植面積達(dá)35.29×104hm2[27]。

    1.2 研究方法

    1.2.1 Biome-BGC模型描述

    Biome-BGC模型是從第一代的Forest-BGC模型發(fā)展而來的,與Forest-BGC模型只適于針葉林生態(tài)系統(tǒng)碳氮水循環(huán)的模擬相比,Biome-BGC模型更具有普適性,可模擬常綠針葉林、常綠闊葉林、落葉針葉林、落葉闊葉林、C3草本植物、C4草本植物和灌木林共7種植被類型的碳、氮、水的循環(huán)過程與交互影響。Biome-BGC模型用C語言編寫,可從美國(guó)蒙拿大大學(xué)森林學(xué)院陸地動(dòng)態(tài)數(shù)值模擬團(tuán)隊(duì)(Numerical Terradynamic Simulation Group, NTSG)的網(wǎng)站下載Biome-BGC4.2版。模型的設(shè)計(jì)遵循物質(zhì)與能量守恒定律,即進(jìn)入系統(tǒng)的物質(zhì)和能量等于留在系統(tǒng)中的物質(zhì)和能量加上離開系統(tǒng)的物質(zhì)和能量[28],具體模擬原理參見White[14]。

    1.2.2 Biome-BGC模型參數(shù)

    Biome-BGC模型的輸入?yún)?shù)包括樣區(qū)參數(shù)、以日為步長(zhǎng)的氣象資料和生態(tài)生理常數(shù)。本研究中橡膠林的樣區(qū)參數(shù)如緯度、海拔、反照率均來自中國(guó)科學(xué)院西雙版納熱帶森林生態(tài)系統(tǒng)研究站,其他參數(shù)來自文獻(xiàn)(表1);1980—2012年的日氣象數(shù)據(jù)從中國(guó)科學(xué)院西雙版納熱帶森林生態(tài)系統(tǒng)研究站申請(qǐng)獲得,1980年以前的從國(guó)家氣象信息中心下載。日氣象數(shù)據(jù)的預(yù)處理方法為,先將潤(rùn)年第366天的數(shù)據(jù)刪除[14],再用線性內(nèi)插法對(duì)缺失值進(jìn)行插補(bǔ)后,調(diào)入MT-CLIM程序生成Biome-BGC所需要的文件;經(jīng)過多次試驗(yàn),生理生態(tài)參數(shù)用系統(tǒng)自帶的參數(shù)模擬效果最佳(表1)。

    表1 橡膠林的Biome-BGC模型樣區(qū)參數(shù)和生理生態(tài)參數(shù)

    DIM:Dimensionless; LAI: leaf area index;*數(shù)據(jù)來源于中國(guó)科學(xué)院西雙版納熱帶雨林生態(tài)系統(tǒng)研究站;Data came from Xishuangbanna Station for Tropical Rainforest Ecosystem Studies;**數(shù)據(jù)來源于White[14]

    1.2.3 Biome-BGC模型同化

    Biome-BGC模型的模擬分兩步實(shí)現(xiàn),第一步是獲得模擬的初始狀態(tài),第二步是正式運(yùn)行模擬程序。初始狀態(tài)的獲得有兩種途徑,一是直接輸入長(zhǎng)期定點(diǎn)觀察數(shù)據(jù)最初的觀測(cè)值,二是運(yùn)行模型自帶的Spin-up程序。在Spin-up過程中,將模擬起點(diǎn)時(shí)的碳氮存量設(shè)為極低的值(葉片的碳存量為0.001 kg C/m2,其它庫(kù)的碳氮存量均為0 kg C/m2),反復(fù)模擬數(shù)千年,直到連續(xù)兩年土壤碳含量的差異小于0.0005 kg C/m2,默認(rèn)系統(tǒng)達(dá)到穩(wěn)定狀態(tài)[28],Spin-up過程結(jié)束[35]??紤]到橡膠種植前對(duì)已有植被的清除性砍伐和燒山,若用達(dá)到穩(wěn)定狀態(tài)時(shí)各個(gè)庫(kù)中的碳氮存量,會(huì)對(duì)橡膠林的固碳過程造成錯(cuò)誤的估計(jì)。因此,參考Thornton等的方法[36],將莖的1%和粗根移入粗死木質(zhì)殘?bào)w庫(kù),隨后的燒山使粗死木質(zhì)殘?bào)w庫(kù)損失50%,葉片、地上凋落物和細(xì)根全部被燒掉,橡膠種植使葉片存有10 g C/m2。氮肥是橡膠林生態(tài)系統(tǒng)中氮素的來源之一,平均每年約3.846 g N/m2[34]。

    Biome-BGC模型首先進(jìn)行物候計(jì)算。模型中規(guī)定,在北半球1月至7月為生長(zhǎng)期,11月至次年1月初為落葉期。在完全落葉期內(nèi),系統(tǒng)的總初級(jí)生產(chǎn)力(GPP)為零[37]。盡管橡膠林在西雙版納出現(xiàn)1—2月集中落葉的現(xiàn)象,但落葉和長(zhǎng)新葉幾乎同時(shí)發(fā)生,尤其在海拔550—600m的區(qū)域,橡膠林的日GPP和最低葉面積指數(shù)并未出現(xiàn)為零的情況[24,38]。綜合考慮,將橡膠林視為一種常綠闊葉林進(jìn)行模擬更為合理。同時(shí)考慮霧水的影響,霧水對(duì)降水的貢獻(xiàn)為全年降水的1.1%[39]。

    1.2.4 模型檢驗(yàn)參數(shù)

    本研究模擬了西雙版納海拔550—600m的橡膠林1959—2012年的碳循環(huán)。模型檢驗(yàn)參數(shù)為同一區(qū)域內(nèi)33年生橡膠林的渦度相關(guān)監(jiān)測(cè)數(shù)據(jù),觀測(cè)鐵塔裝在中國(guó)科學(xué)院西雙版納熱帶森林生態(tài)系統(tǒng)研究站內(nèi)單一橡膠林人工群落(21°56′N,101°15′E,海拔580 m)中[24],觀測(cè)時(shí)間為2010年7月—2011年6月。

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

    統(tǒng)計(jì)分析用SPSS 16.0,模擬值與實(shí)測(cè)值的相關(guān)性采用線性回歸分析,差異顯著性水平設(shè)為а=0.05。用SigmaPlot 12.5軟件制作圖表。

    2 結(jié)果與分析

    2.1 模型檢驗(yàn)

    用Biome-BGC模型模擬了西雙版納1959—2012年橡膠林的碳循環(huán),并從模擬結(jié)果中提取了與驗(yàn)證參數(shù)相對(duì)應(yīng)林齡橡膠林的GPP、總呼吸(Rt)和NEE,然后參照Hidy[19]的方法,從提取年份至2011年逐年求平均,以消除氣象數(shù)據(jù)年際間波動(dòng)的影響。與原模型相比,同化后的Biome-BGC模型對(duì)橡膠林碳循環(huán)的模擬效果相對(duì)較好。日GPP模擬值與實(shí)測(cè)值的相關(guān)系數(shù)由0.58提高到0.66(表2),誤差源于對(duì)干季中期至雨季初期(1月中旬至次年5月)的模擬結(jié)果偏低,對(duì)雨季末期(9—10月)的模擬偏低(圖1);模型同化后對(duì)日Rt的模擬效果最好,模擬值與實(shí)測(cè)值的相關(guān)性達(dá)0.92(表2),模擬值比實(shí)測(cè)值稍高(圖1);受GPP模擬值偏低和Rt模擬值偏高的共同影響,模型同化后對(duì)日NEE(NEE為負(fù)值表示碳吸收,NEE為正值表示碳釋放)的模擬效果相對(duì)較差,模擬值與實(shí)測(cè)值的相關(guān)性僅0.39,模擬值低于實(shí)測(cè)值(表2,圖1)。從年積累量來看(表3),模型同化后年GPP的相對(duì)誤差由-21.76%降低至-1.62%,模擬值比實(shí)測(cè)值低41.61 g C m-2a-1,總體精度為98.37%;年Rt的相對(duì)誤差由-11.74%降低到9.04%,模擬值比實(shí)測(cè)值高205.82 g C m-2a-1,整體精度約90%;年NEE模擬值比實(shí)測(cè)值低157.35 g C m-2a-1,相對(duì)誤差為63.7%,這可能與Biome-BGC模型未能模擬干膠產(chǎn)量有關(guān)。33a林齡橡膠林的干膠碳儲(chǔ)量為139 g C m-2a-1[40],因此若考慮干膠固碳量,模型同化后的模擬結(jié)果與實(shí)測(cè)值十分接近。表明同化后的Biome-BGC模型可用于西雙版納橡膠林生態(tài)系統(tǒng)碳循環(huán)的模擬。

    表2 橡膠林日GPP、Rt、NEE實(shí)測(cè)值與模擬值的相關(guān)性分析

    Y:模擬值,simulated valued;X:實(shí)測(cè)值,measured valued

    表3 橡膠林年GPP、Rt、NEE實(shí)測(cè)值與模擬值的誤差分析

    MV:實(shí)測(cè)值,measured value;SV:模擬值,simulated value;RD: 相對(duì)差異,relative difference;RE:相對(duì)誤差;relative error

    圖1 Biome-BGC同化后模擬的橡膠林日GPP(a)、Rt(b)和NEE(c)與實(shí)測(cè)值比較Fig.1 Comparison of simulated daily GPP (a), Rt (b) and NEE (c) in rubber plantation using the assimilated Biome-BGC model with the observed values

    2.2 霧水對(duì)橡膠林碳循環(huán)模擬的影響

    氣孔是植物水氣傳輸?shù)闹匾T戶,它與植物的光合作用、呼吸作用等重要生理生態(tài)過程密切相關(guān)。在Biome-BGC模型中,植物可利用水分的多少直接影響氣孔導(dǎo)度。在西雙版納地區(qū),橡膠林中有長(zhǎng)達(dá)188d的霧日,霧日橡膠林冠層對(duì)霧水的截留量占全年降水量的1.1%,極大的彌補(bǔ)了干季降雨量的不足[39]。為了摸清霧水對(duì)橡膠林碳循環(huán)的影響,本研究模擬了降水中不加入霧水和加入霧水兩種情況。結(jié)果表明,在降水中加入霧水后,GPP和Rt同時(shí)增加,NEE沒有明顯的變化(圖2)。但在降水中加入霧水后,NEE異常值大為減少,NEE在1969年以后即保持以碳吸收為主,表明霧水對(duì)橡膠林碳收支的模擬具有促進(jìn)作用,在后續(xù)分析中均使用降水中加入霧水的模擬結(jié)果。

    圖2 霧水對(duì)橡膠林年GPP、Rt和NEE的影響分析Fig.2 The effects of fog on annual GPP,Rt and NEE in rubber plantation

    2.3 橡膠林碳收支的模擬

    從圖3可以看出,橡膠林自養(yǎng)呼吸隨模擬年限的增加而逐漸增加,在模擬進(jìn)行5a時(shí)自養(yǎng)呼吸即由最初的32 g C m-2a-1增加到2 108 g C m-2a-1,之后維持在2 100—2 600 g C m-2a-1之間。異養(yǎng)呼吸的變化規(guī)律與自養(yǎng)呼吸正好相反,在模擬的前4a里較高,平均約1 276 g C m-2a-1。4a后異養(yǎng)呼吸迅速降低并在第10年左右達(dá)到穩(wěn)定狀態(tài),平均約496 g C m-2a-1,為前4a的2/5。橡膠林總呼吸隨時(shí)間的變化與自養(yǎng)呼吸的變化基本一致,前4a由1959年的1 492 g C m-2a-1迅速增加到1962年的2 374 g C m-2a-1, 4a后總呼吸保持在2 776 g C m-2a-1左右。

    圖3 1959—2012年橡膠林年Rt、異養(yǎng)呼吸(Rh)、自養(yǎng)呼吸(Ra)的模擬Fig.3 Simulation of annual Rt, heterotrophic respiration (Rh), autotrophic respiration (Ra) in rubber plantations during 1959—2012

    橡膠林年GPP和年NPP均在模擬第3—4年快速增加,之后達(dá)到穩(wěn)定狀態(tài),分別保持在-266—-970 g C m-2a-1和-2413—-3226 g C m-2a-1之間(圖4)。年NEE在模擬進(jìn)行的前8年里平均357 g C m-2a-1,且絕大多數(shù)年份大于0,表明該時(shí)期橡膠林以碳排放為主(共釋放碳2855 g C/m2);8年以后除極少數(shù)年份外,NEE均小于0,平均約-146g C m-2a-1,表現(xiàn)為碳匯。由前面的分析可知,年NEE模擬值比觀測(cè)值低63.7%,若將年NEE模擬值增加63.7%,橡膠林仍然是在模擬進(jìn)行8a左右才由碳源轉(zhuǎn)為碳匯(圖4)。由此我們可以初步判斷,橡膠種植的前8a里為碳源,8a以后由碳源轉(zhuǎn)為碳匯。若以40a為橡膠林的最佳固碳周期[23],在適宜種植區(qū)橡膠林一生中可固定碳1835 g C/m2,表現(xiàn)為碳匯。

    通過對(duì)呼吸作用(Ra和Rh)和固碳作用(GPP、NPP)的綜合分析發(fā)現(xiàn),橡膠林在模擬初期之所以表現(xiàn)為碳源,主要是由較高的異養(yǎng)呼吸和較弱的固碳能力(GPP、NPP均較低)引起。

    圖4 1959—2012年橡膠林年NEE、GPP和NPP的模擬Fig.4 Simulation of annual NEE, GPP and NPP in rubber plantations during 1959—2012

    3 討論

    本研究根據(jù)橡膠種植的實(shí)際情況,考慮了橡膠種植對(duì)已有植被的清除性砍伐和燒山行為、施肥等因素,對(duì)模型進(jìn)行了參數(shù)同化,同化后的Biome-BGC對(duì)橡膠林碳收支的模擬效果獲得了較大的改善。對(duì)年GPP的模擬精度由78.34%提高到98.37%,模擬值與實(shí)測(cè)值的相關(guān)性為0.66。從GPP的日變化來看,對(duì)日GPP的模擬值在干季中期至雨季初期(1月中旬至次年5月)偏低,在雨季末期(9—10月)偏高。這可能與橡膠林自身的水分利用策略有關(guān),即橡膠林在干季轉(zhuǎn)向深層土壤吸收[41- 42]。在雨季,橡膠林主要利用表層0—30 cm土壤中的水分,對(duì)70 cm以下土層中的水分利用極少[43];而在干季后期,因70 cm以上土層的水分被大量消耗,橡膠林向深層土壤吸水深度可達(dá)2米以下[41]。同時(shí),雨季來臨第一次降雨過后,橡膠樹就能迅速的吸收表層土壤中的水分,使得橡膠林能在干季生長(zhǎng)發(fā)芽。但在Biome-BGC模型中,一方面植物有效吸水根系深度固定為一層;另一方面氣孔是植物-土壤-大氣系統(tǒng)中水氣交換的唯一門戶,土壤水勢(shì)直接影響氣孔導(dǎo)度。因此,當(dāng)根系吸水深度被固定為1米時(shí),在雨季模型極可能高估了橡膠林對(duì)深層土壤水的利用率,從而對(duì)GPP造成高估;在干季則相反,模型中低估了深層土壤中可利用水分的量并限制了光合作用的進(jìn)行,對(duì)GPP的模擬偏低。此外,在土壤含水量較低的情況下,模型中土壤水勢(shì)對(duì)水分的響應(yīng)是突變型的,而實(shí)際上土壤水的滲透有一定的過程,即滯后性。Hidy等[19]利用Biome-BGC模型模擬草地生態(tài)系統(tǒng)的碳收支時(shí),也發(fā)現(xiàn)了這一問題,對(duì)模型進(jìn)行了改進(jìn),并取得了較好的模擬效果,但Hidy等人并未將改進(jìn)后的源代碼公開,今后的工作中我們將會(huì)在這方面做出更多的努力。

    模型同化后對(duì)Rt的模擬效果較好,模擬值與渦度相關(guān)監(jiān)測(cè)值的相關(guān)性達(dá)0.91,總體模擬精度為91%。但由于對(duì)GPP的模擬結(jié)果偏低和對(duì)Rt的模擬結(jié)果偏高,導(dǎo)致模型對(duì)NEE的模擬誤差較大,模擬值比實(shí)測(cè)值低157.35 g C m-2a-1,相對(duì)誤差為63.7%。NEE模擬值偏低可能源于模型未能模擬干膠產(chǎn)量,若考慮干膠產(chǎn)量(33年生橡膠林約139 g C m-2a-1)[40]的影響, NEE模擬值與實(shí)測(cè)值十分接近,表明同化后的模型可以用來進(jìn)行橡膠林碳收支的模擬。

    水分是影響B(tài)iome-BGC模型模擬結(jié)果的一個(gè)重要方面,鑒于西雙版納具有長(zhǎng)達(dá)188 d的霧日和霧水對(duì)降水的貢獻(xiàn)[39],在模擬橡膠林的碳循環(huán)時(shí)還考慮了霧水的作用。分析發(fā)現(xiàn),雖然霧水對(duì)NEE沒有明顯的影響,但卻有效減少了NEE年際間的異常值,使橡膠林在模擬8a后即表現(xiàn)為相對(duì)穩(wěn)定的碳匯。表明霧水對(duì)橡膠林碳循環(huán)的模擬具有一定的促進(jìn)作用,但不是主要影響因子。

    由NEE的變化規(guī)律可知,橡膠林在模擬的前8a因異養(yǎng)呼吸過高以碳排放為主,表現(xiàn)為碳源;8a后異養(yǎng)呼吸逐漸降低,至第9年達(dá)到相對(duì)穩(wěn)定的狀態(tài),平均為500 g C m-2a-1,橡膠林由碳源轉(zhuǎn)為相對(duì)穩(wěn)定的碳匯。這一結(jié)果與Amiro等[44]對(duì)受干擾的林齡系列的研究結(jié)論基本相符,其主要結(jié)論有:1)在演替初期由于光合作用弱于呼吸作用,生態(tài)系統(tǒng)表現(xiàn)為碳源;2)從干擾期開始至干擾后10a左右,生態(tài)系統(tǒng)即可恢復(fù)為碳匯,之后碳匯維持不變;3)干擾系列由碳源到碳匯的轉(zhuǎn)變是由異養(yǎng)呼吸與光合作用的相對(duì)變化引起。這進(jìn)一步證明同化后的Biome-BGC模型適宜于橡膠林碳收支的模擬。若以40a更新周期計(jì)算,在橡膠適宜種植區(qū)橡膠林一生中可固定碳1 835 g C/m2,是一個(gè)弱的碳匯。但與西雙版納熱帶雨林相同周期固碳6 720 g C/m2相比[45],建立于熱帶雨林上的橡膠林在整個(gè)生命周期中都表現(xiàn)為碳源。

    4 結(jié)論

    在橡膠適宜種植區(qū),橡膠林在種植初期的前8a里以碳排放為主,主要由異養(yǎng)呼吸過高引起,8年以后轉(zhuǎn)為相對(duì)穩(wěn)定的碳匯;若以40a為橡膠林的最佳固碳周期,橡膠林一生中可固定碳1835 g C/m2,是一個(gè)弱的碳匯。但與熱帶雨林相比,極可能在整個(gè)生命周期中都表現(xiàn)為碳匯。因此在橡膠種植過程中,一是要對(duì)老膠林進(jìn)行有計(jì)劃的更新,務(wù)必考慮橡膠林的林齡結(jié)構(gòu),盡量避免盲目更新破壞橡膠林生態(tài)系統(tǒng)現(xiàn)有的碳平衡;二是必須重視對(duì)熱帶雨林的保護(hù),不宜再破壞熱帶雨林建立新的橡膠林。從而最大限度的降低橡膠種植對(duì)區(qū)域性碳收支的影響,實(shí)現(xiàn)經(jīng)濟(jì)和生態(tài)環(huán)境保護(hù)的協(xié)調(diào)發(fā)展。

    致謝:感謝中國(guó)科學(xué)院西雙版納熱帶森林生態(tài)系統(tǒng)研究站,及中國(guó)科學(xué)院西雙版納熱帶植物園唐建維研究員、沙麗清研究員、范澤鑫副研究員、宋清海副研究員和中國(guó)科學(xué)院南京土壤研究所的遞超普等對(duì)本工作所給予的大力支持。

    [1] Stocker T F, 秦大河, Plattner G K, Tignor M M B, Allen S K, Boschung J, Nauels A, Xia Y, Bex V, Midgley P M. IPCC 2013: 決策者摘要. 政府間氣候變化專門委員會(huì)第五次評(píng)估報(bào)告第一工作組報(bào)告氣候變化2013: 自然科學(xué)基礎(chǔ). 英國(guó)劍橋, 美國(guó)紐約: 劍橋大學(xué)出版社, 2013: 9- 10.

    [2] Yang J C, Huang J H, Tang J W, Pan Q M, Han X G. Carbon sequestration in rubber tree plantation established on former arable lands in Xishuangbanna, SW China. Acta Phytoecologica Sinica, 2005, 29(2): 296- 303.

    [3] Li H M, Ma Y X, Aide T M, Liu W J. Past, present and future land-use in Xishuangbanna, China and the implications for carbon dynamics. Forest Ecology and Management, 2008, 255(1): 16- 24.

    [4] Ziegler A D, Fox J M, Xu J C. The rubber juggernaut. Science, 2009, 324(5930): 1024- 1025.

    [5] de Blécourt M, Brumme R, Xu J C, Corre M D, Veldkamp E. Soil carbon stocks decrease following conversion of secondary forests to rubber (Heveabrasiliensis) plantations. PLoS One, 2013, 8(7): e69357.

    [6] de Blécourt M, H?nsel V M, Brumme R, Corre M D, Veldkamp E. Soil redistribution by terracing alleviates soil organic carbon losses caused by forest conversion to rubber plantation. Forest Ecology and Management, 2014, 313: 26- 33.

    [7] Blagodatsky S, Xu J C, Cadisch G. Carbon balance of rubber (Heveabrasiliensis) plantations: a review of uncertainties at plot, landscape and production level. Agriculture, Ecosystems & Environment, 2016, 221: 8- 19.

    [8] Yang X Q, Blagodatsky S, Lippe M, Liu F, Hammond J, Xu J C, Cadisch G. Land-use change impact on time-averaged carbon balances: rubber expansion and reforestation in a biosphere reserve, South-West China. Forest Ecology and Management, 2016, 372: 149- 163.

    [9] 馮險(xiǎn)峰, 劉高煥, 陳述彭, 周文佐. 陸地生態(tài)系統(tǒng)凈第一性生產(chǎn)力過程模型研究綜述. 自然資源學(xué)報(bào), 2004, 19(3): 369- 378.

    [10] Running S W, Gower S T. FOREST-BGC, a general model of forest ecosystem processes for regional applications. II. Dynamic carbon allocation and nitrogen budgets. Tree Physiology, 1991, 9(1/2): 147- 160.

    [11] Running S W, Hunt E R Jr. Generalization of a forest ecosystem process model for other biomes, BIOME-BGC, and an application for global-scale models//Ehleringer J R, Field C B, eds. Scaling Physiological Processes: Leaf to Globe. San Diego: Academic Press, 1993: 141- 158.

    [12] Pennock D J, Frick A H. The role of field studies in landscape-scale applications of process models: an example of soil redistribution and soil organic carbon modeling using CENTURY. Soil and Tillage Research, 2001, 58(3/4): 183- 191.

    [13] Ooba M, Wang Q X, Murakami S, Kohata K. Biogeochemical model (BGC-ES) and its basin-level application for evaluating ecosystem services under forest management practices. Ecological Modelling, 2010, 221(16): 1979- 1994.

    [14] White M A, Thornton P E, Running S W, Nemani R R. Parameterization and sensitivity analysis of the BIOME-BGC terrestrial ecosystem model: net primary production controls. Earth Interactions, 2000, 4(3): 1- 85.

    [15] Wang W L, Ichii K, Hashimoto H, Michaelis A R, Thornton P E, Law B E, Nemani R R. A hierarchical analysis of terrestrial ecosystem model Biome-BGC: equilibrium analysis and model calibration. Ecological Modelling, 2005, 220(17): 2009- 2023.

    [16] Cienciala E, Tatarinov F A. Application of BIOME-BGC model to managed forests: 2. comparison with long-term observations of stand production for major tree species. Forest Ecology and Management, 2006, 237(1/3): 252- 266.

    [17] Peckham S D, Gower S T, Perry C H, Wilson B T, Stueve K M. Modeling harvest and biomass removal effects on the forest carbon balance of the Midwest, USA. Environmental Science & Policy, 2013, 25: 22- 35.

    [18] 曾慧卿, 劉琪璟, 馮宗煒, 王效科, 馬澤清. 基于BIOME-BGC模型的紅壤丘陵區(qū)濕地松(Pinuselliottii)人工林GPP和NPP. 生態(tài)學(xué)報(bào), 2008, 28(11): 5314- 5321.

    [19] Hidy D, Barcza Z, Haszpra L, Churkina G, Pintér K, Nagy Z. Development of the Biome-BGC model for simulation of managed herbaceous ecosystems. Ecological Modelling, 2012, 226: 99- 119.

    [20] Xu J C, Grumbine R E, Becksch?fer P. Landscape transformation through the use of ecological and socioeconomic indicators in Xishuangbanna, Southwest China, Mekong Region. Ecological Indicators, 2014, 36: 749- 756.

    [21] Li H M, Aide T M, Ma Y X, Liu W J, Cao M. Demand for rubber is causing the loss of high diversity rain forest in SW China. Biodiversity and Conservation, 2007, 16(6): 1731- 1745.

    [22] Huang X H, Yuan H, Yu F K, Li X Y, Liang Q B, Yao P, Shao H B. Spatial-temporal succession of the vegetation in Xishuangbanna, China during 1976- 2010: a case study based on RS technology and implications for eco-restoration. Ecological Engineering, 2014, 70: 255- 262.

    [23] Nizami S M, Zhang Y P, Sha L Q, Zhao W, Zhang X. Managing carbon sinks in rubber (Heveabrasilensis) plantation by changing rotation length in SW China. PLoS One, 2014, 9(12): e115234.

    [24] Song Q H, Tan Z H, Zhang Y P, Sha L Q, Deng X B, Deng Y, Zhou W J, Zhao J F, Zhao J B, Zhang X, Zhao W, Yu G R, Sun X M, Liang N S, Yang L Y. Do the rubber plantations in tropical China act as large carbon sinks?. iForest, 2014, 7: 42- 47.

    [25] 宋清海, 張一平. 西雙版納地區(qū)人工橡膠林生物量、固碳現(xiàn)狀及潛力. 生態(tài)學(xué)雜志, 2010, 29(10): 1887- 1891.

    [26] 朱華. 論滇南西雙版納的森林植被分類. 云南植物研究, 2007, 29(4): 377- 387.

    [27] 梅岑岑, 馬友鑫, 李紅梅, 劉文俊, 孫燕瓷, 湯柔馨. 西雙版納土地利用景觀格局分析. 云南大學(xué)學(xué)報(bào): 自然科學(xué)版(S2), 2015, 36(S1): 363- 369.

    [28] Thornton P E, Rosenbloom N A. Ecosystem model spin-up: estimating steady state conditions in a coupled terrestrial carbon and nitrogen cycle model. Ecological Modelling, 2005, 189(1/2): 25- 48.

    [29] 唐炎林, 鄧曉保, 李玉武, 張順賓. 西雙版納不同林分土壤機(jī)械組成及其肥力比較. 中南林業(yè)科技大學(xué)學(xué)報(bào): 自然科學(xué)版, 2007, 27(1): 70- 75.

    [30] Li H M, Ma Y X, Liu W J, Liu W J. Soil changes induced by rubber and tea plantation establishment: comparison with tropical rain forest soil in Xishuangbanna, SW China. Environmental Management, 2012, 50(5): 837- 848.

    [31] 沙麗清, 鄭征, 馮志立, 劉玉洪, 劉文杰, 孟盈, 李明銳. 西雙版納熱帶季節(jié)雨林生態(tài)系統(tǒng)氮的生物地球化學(xué)循環(huán)研究. 植物生態(tài)學(xué)報(bào), 2002, 26(6): 689- 694.

    [32] Han Y G, Fan Y T, Yang P L, Wang X X, Wang Y J, Tian J X, Xu L, Wang C Z. Net anthropogenic nitrogen inputs (NANI) index application in Mainland China. Geoderma, 2014, 231: 87- 94.

    [33] Boring L R, Swank W T, Waide J B, Henderson G S. Sources, fates, and impacts of nitrogen inputs to terrestrial ecosystems: review and synthesis. Biogeochemistry, 1988, 6(2): 119- 159.

    [34] Petsri S, Chidthaisong A, Pumijumnong N, Wachrinrat C. Greenhouse gas emissions and carbon stock changes in rubber tree plantations in Thailand from 1990 to 2004. Journal of Cleaner Production, 2013, 52: 61- 70.

    [35] Pietsch S A, Hasenauer H, Thornton P E. BGC-model parameters for tree species growing in central European forests. Forest Ecology and Management, 2005, 211(3): 264- 295.

    [36] Thornton P E, Law B E, Gholz H L, Clark K L, Falge E, Ellsworth D S, Goldstein A H, Monson R K, Hollinger D, Falk M, Chen J, Sparks J P. Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests. Agricultural and Forest Meteorology, 2002, 113(1/4): 185- 222.

    [37] White M A, Thornton P E, Running S W. A continental phenology model for monitoring vegetation responses to interannual climatic variability. Global Biogeochemical Cycles, 1997, 11(2): 217- 234.

    [38] 賈開心. 西雙版納三葉橡膠林生長(zhǎng)隨海拔高度變化研究[D]. 景洪: 中國(guó)科學(xué)院研究生院(西雙版納熱帶植物園), 2006.

    [39] 劉文杰, 張一平, 劉玉洪, 李紅梅, 段文平. 熱帶季節(jié)雨林和人工橡膠林林冠截留霧水的比較研究. 生態(tài)學(xué)報(bào), 2003, 23(11): 2379- 2386.

    [40] 孫燕瓷, 馬友鑫, 曹坤芳, 李紅梅, 何建芬, 劉昌平, 劉洪平, 李斌. 西雙版納長(zhǎng)期橡膠種植對(duì)土壤固碳和氮的影響研究. 土壤通報(bào), 2015, 46(2): 412- 419.

    [41] Guardiola-claramonte M, Troch P A, Ziegler A D, Giambelluca T W, Vogler J B, Nullet M A. Local hydrologic effects of introducing non-native vegetation in a tropical catchment. Ecohydrology, 2008, 1(1): 13- 22.

    [42] Tan Z H, Zhang Y P, Song Q H, Liu W J, Deng X B, Tang J W, Deng Y, Zhou W J, Yang L Y, Yu G R, Sun X M, Liang N S. Rubber plantations act as water pumps in tropical China. Geophysical Research Letters, 2011, 38: L24406.

    [43] Liu W J, Li J T, Lu H J, Wang P Y, Luo Q P, Liu W Y, Li H M. Vertical patterns of soil water acquisition by non-native rubber trees (Heveabrasiliensis) in Xishuangbanna, southwest China. Ecohydrology, 2014, 7(4): 1234- 1244.

    [44] Amiro B D, Barr A G, Barr J G, Black T A, Bracho R, Brown M, Chen J, Clark K L, Davis K J, Desai A R, Dore S, Engel V, Fuentes J D, Goldstein A H, Goulden M L, Kolb T E, Lavigne M B, Law B E, Margolis H A, Martin T, McCaughey J H, Misson L, Montes-Helu M, Noormets A, Randerson J T, Starr G, Xiao J. Ecosystem carbon dioxide fluxes after disturbance in forests of North America. Journal of Geophysical Research, 2010, 115(G4): G00K02.

    [45] Zhang Y P, Tan Z H, Song Q H, Yu G R, Sun X M. Respiration controls the unexpected seasonal pattern of carbon flux in an Asian tropical rain forest. Atmospheric Environment, 2010, 44(32): 3886- 3893.

    SimulationofcarbonbudgetinrubberplantationsinXishuangbannabasedontheBiome-BGCmodel

    SUN Yanci2, MA Youxin1,*, CAO Kunfang3, SHEN Jinxiang4, ZHANG Yiping1, MEI Cencen5, LIU Wenjun1

    1KeyLaboratoryofTropicalForestEcology,XishuangbannaTropicalBotanicalGarden,ChineseAcademyofSciences,Menglun666303,China2YunnanCollegeofBusinessManagement,Kunming650106,China3StateKeyLaboratoryforConservationandUtilizationofSubtropicalAgro-bioresources,andCollegeofForestry,GuangxiUniversity,Nanning530004,China4DepartmentofDigitalLandandLandManagement,YunnanLandandResourcesVocationalCollege,Kunming652501,China5Shenzhen-HongkongInstitutionofIndustry,Education&ResearchEnvironmentalTechniqueCenter,Shenzhen518057,China

    中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA05050206)

    2016- 06- 08; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 22

    10.5846/stxb201606081104

    *通訊作者Corresponding author.E-mail: may@xtbg.ac.cn

    孫燕瓷,馬友鑫, 曹坤芳,沈金祥,張一平,梅岑岑,劉文俊.基于Biome-BGC模型的西雙版納橡膠林碳收支模擬.生態(tài)學(xué)報(bào),2017,37(17):5732- 5741.

    Sun Y C, Ma Y X, Cao K F, Shen J X, Zhang Y P, Mei C C, Liu W J.Simulation of carbon budget in rubber plantations in Xishuangbanna based on the Biome-BGC model.Acta Ecologica Sinica,2017,37(17):5732- 5741.

    猜你喜歡
    橡膠林碳循環(huán)西雙版納
    食物網(wǎng)與碳循環(huán)
    海南橡膠林生態(tài)系統(tǒng)凈碳交換物候特征
    橡膠林
    西雙版納 一個(gè)與冬天背道而馳的地方
    西雙版納橡膠林土壤有機(jī)碳分布特征研究
    南京城市系統(tǒng)碳循環(huán)與碳平衡分析
    多措并舉構(gòu)筑綠色低碳循環(huán)發(fā)展的制造業(yè)體系
    YUNNAN
    西雙版納散記
    西雙版納叢林歷險(xiǎn)記
    欧美成人精品欧美一级黄| 黄色日韩在线| 此物有八面人人有两片| 天堂√8在线中文| 18禁在线无遮挡免费观看视频 | 人人妻人人澡欧美一区二区| 久久久久久久午夜电影| 欧美色视频一区免费| 久久久久精品国产欧美久久久| 亚洲国产精品合色在线| 国内揄拍国产精品人妻在线| 色播亚洲综合网| 在线看三级毛片| 欧美最黄视频在线播放免费| 在线观看66精品国产| 成人特级黄色片久久久久久久| 三级男女做爰猛烈吃奶摸视频| h日本视频在线播放| 伦理电影大哥的女人| 国产视频一区二区在线看| 国产高清不卡午夜福利| 国产精品久久久久久久电影| 日韩欧美精品v在线| 少妇丰满av| 一进一出好大好爽视频| 色吧在线观看| 十八禁网站免费在线| 女生性感内裤真人,穿戴方法视频| 成人特级av手机在线观看| 麻豆国产97在线/欧美| 男女之事视频高清在线观看| 久久精品国产99精品国产亚洲性色| 亚洲精华国产精华液的使用体验 | 精品久久久久久久久av| 秋霞在线观看毛片| 亚洲中文日韩欧美视频| 午夜影院日韩av| 淫秽高清视频在线观看| 欧美中文日本在线观看视频| 99久久九九国产精品国产免费| 成人三级黄色视频| 嫩草影院入口| 十八禁国产超污无遮挡网站| 国产免费一级a男人的天堂| 国产成年人精品一区二区| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久久丰满| 亚洲精品一卡2卡三卡4卡5卡| 一区福利在线观看| 99riav亚洲国产免费| 又爽又黄a免费视频| av福利片在线观看| 99久久精品热视频| 麻豆国产av国片精品| 国产成人a区在线观看| 免费观看在线日韩| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 国产成人a区在线观看| 白带黄色成豆腐渣| av在线蜜桃| 少妇被粗大猛烈的视频| av在线播放精品| 欧美成人免费av一区二区三区| 黄色配什么色好看| 精品乱码久久久久久99久播| 最近2019中文字幕mv第一页| 一个人看视频在线观看www免费| 村上凉子中文字幕在线| 卡戴珊不雅视频在线播放| 国产91av在线免费观看| 免费av不卡在线播放| 日韩欧美一区二区三区在线观看| 一本一本综合久久| 18+在线观看网站| 神马国产精品三级电影在线观看| 校园人妻丝袜中文字幕| 亚洲精品乱码久久久v下载方式| 你懂的网址亚洲精品在线观看 | 91久久精品电影网| 少妇的逼好多水| 1024手机看黄色片| 搡老熟女国产l中国老女人| 自拍偷自拍亚洲精品老妇| 成人漫画全彩无遮挡| a级毛色黄片| 熟妇人妻久久中文字幕3abv| 在线观看一区二区三区| 亚洲最大成人手机在线| 久久久久久久久中文| 99国产极品粉嫩在线观看| 日韩欧美免费精品| 久久草成人影院| 日韩,欧美,国产一区二区三区 | 又粗又爽又猛毛片免费看| 九九热线精品视视频播放| 久久人人爽人人片av| 波多野结衣巨乳人妻| 亚洲激情五月婷婷啪啪| 亚洲av一区综合| 亚洲欧美精品综合久久99| 大香蕉久久网| 国产成年人精品一区二区| 国产真实乱freesex| 午夜免费男女啪啪视频观看 | 久久精品国产亚洲av天美| av天堂中文字幕网| 久久精品影院6| 最近的中文字幕免费完整| 国产亚洲精品综合一区在线观看| 在线观看一区二区三区| av女优亚洲男人天堂| 国内精品久久久久精免费| 亚洲中文字幕一区二区三区有码在线看| 天堂√8在线中文| 一级黄片播放器| 中国美白少妇内射xxxbb| 美女大奶头视频| 国产精品久久久久久av不卡| 一级毛片aaaaaa免费看小| 男人舔奶头视频| 成人漫画全彩无遮挡| 色av中文字幕| 99视频精品全部免费 在线| 亚洲av成人av| 亚洲一区高清亚洲精品| 欧美三级亚洲精品| 男女边吃奶边做爰视频| 又爽又黄无遮挡网站| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 特大巨黑吊av在线直播| 欧美区成人在线视频| 在线免费观看的www视频| 亚洲国产精品成人久久小说 | 免费看a级黄色片| 亚洲欧美日韩卡通动漫| 直男gayav资源| 看免费成人av毛片| 九九在线视频观看精品| 男人和女人高潮做爰伦理| 久久精品国产99精品国产亚洲性色| 国产探花在线观看一区二区| 国产精品亚洲一级av第二区| 久久久成人免费电影| 日本爱情动作片www.在线观看 | 欧美最黄视频在线播放免费| 日韩,欧美,国产一区二区三区 | 亚洲人成网站在线播| 日韩大尺度精品在线看网址| 免费无遮挡裸体视频| 一进一出抽搐gif免费好疼| 欧美日韩精品成人综合77777| 国产一区二区亚洲精品在线观看| 亚洲精品国产av成人精品 | 国产午夜精品久久久久久一区二区三区 | 免费看a级黄色片| 久久久精品94久久精品| 国模一区二区三区四区视频| 蜜桃久久精品国产亚洲av| 色综合色国产| 国产黄a三级三级三级人| 色视频www国产| 三级经典国产精品| 我要看日韩黄色一级片| 真实男女啪啪啪动态图| 午夜福利18| 天天躁日日操中文字幕| 亚洲人成网站在线播| 国产精品国产高清国产av| 精品乱码久久久久久99久播| 欧美成人免费av一区二区三区| 欧美三级亚洲精品| 狠狠狠狠99中文字幕| 男女那种视频在线观看| 色噜噜av男人的天堂激情| 久久国内精品自在自线图片| 欧美绝顶高潮抽搐喷水| 深夜精品福利| 欧美丝袜亚洲另类| 夜夜看夜夜爽夜夜摸| 男女边吃奶边做爰视频| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 一级毛片久久久久久久久女| 国产又黄又爽又无遮挡在线| 精品人妻视频免费看| 最后的刺客免费高清国语| 国产在线男女| 色吧在线观看| 99热这里只有是精品在线观看| 欧美最新免费一区二区三区| 亚洲av一区综合| 老师上课跳d突然被开到最大视频| 日本免费一区二区三区高清不卡| 一区二区三区免费毛片| 午夜视频国产福利| 精品国产三级普通话版| 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 九九热线精品视视频播放| 别揉我奶头 嗯啊视频| 国内少妇人妻偷人精品xxx网站| 女人十人毛片免费观看3o分钟| 天堂网av新在线| 国产精品无大码| 男插女下体视频免费在线播放| 国产单亲对白刺激| 99久久九九国产精品国产免费| 哪里可以看免费的av片| 欧美极品一区二区三区四区| 亚洲中文字幕日韩| 欧美zozozo另类| av专区在线播放| 在线观看一区二区三区| 麻豆av噜噜一区二区三区| 99视频精品全部免费 在线| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品| 欧美+亚洲+日韩+国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲三级黄色毛片| 亚洲av.av天堂| 麻豆乱淫一区二区| 不卡一级毛片| 午夜老司机福利剧场| 国产亚洲精品综合一区在线观看| 天堂√8在线中文| 99热这里只有精品一区| 国产精品无大码| 好男人在线观看高清免费视频| 亚洲最大成人中文| 美女xxoo啪啪120秒动态图| 国产综合懂色| 男女下面进入的视频免费午夜| 国产成人aa在线观看| 婷婷精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 亚洲国产精品国产精品| 美女大奶头视频| 最近在线观看免费完整版| 国产成人freesex在线 | 免费无遮挡裸体视频| 亚洲精品国产成人久久av| 久久久久精品国产欧美久久久| 国产精品福利在线免费观看| 久久久色成人| 女同久久另类99精品国产91| 99久久无色码亚洲精品果冻| 中文资源天堂在线| 精品久久久久久久久av| 亚洲在线观看片| 日韩中字成人| 最近2019中文字幕mv第一页| 麻豆久久精品国产亚洲av| 亚洲av第一区精品v没综合| 日韩一区二区视频免费看| 搡女人真爽免费视频火全软件 | 日韩欧美 国产精品| 91久久精品国产一区二区成人| 久久久欧美国产精品| 美女黄网站色视频| 日本三级黄在线观看| 蜜臀久久99精品久久宅男| 日本一二三区视频观看| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 十八禁国产超污无遮挡网站| 99久久精品热视频| 亚洲美女搞黄在线观看 | 级片在线观看| 久久婷婷人人爽人人干人人爱| 亚洲成a人片在线一区二区| 国产午夜福利久久久久久| 久久韩国三级中文字幕| 黄色一级大片看看| 狠狠狠狠99中文字幕| 在线国产一区二区在线| 国产一区二区激情短视频| 国产91av在线免费观看| 村上凉子中文字幕在线| 桃色一区二区三区在线观看| 欧美+亚洲+日韩+国产| 99久国产av精品国产电影| 成年版毛片免费区| 天美传媒精品一区二区| 欧美日本视频| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 国产不卡一卡二| 国产视频内射| av在线观看视频网站免费| 国产真实乱freesex| 国产亚洲精品久久久久久毛片| 日本a在线网址| 国产精品人妻久久久久久| 久久人人爽人人爽人人片va| 国产色婷婷99| 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 久久精品国产清高在天天线| 久久久久精品国产欧美久久久| 我的老师免费观看完整版| 乱人视频在线观看| 欧美一级a爱片免费观看看| 国内揄拍国产精品人妻在线| 成年av动漫网址| 村上凉子中文字幕在线| 美女高潮的动态| 亚洲av中文av极速乱| 别揉我奶头~嗯~啊~动态视频| 在线观看66精品国产| 综合色丁香网| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品精品国产色婷婷| 国产精品电影一区二区三区| 黑人高潮一二区| 直男gayav资源| 国产不卡一卡二| 又爽又黄无遮挡网站| 国产高清视频在线播放一区| 久久久久久久午夜电影| 欧美极品一区二区三区四区| 国产 一区精品| 久久人人爽人人片av| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区精品| av在线蜜桃| 给我免费播放毛片高清在线观看| 日本黄色视频三级网站网址| 欧美成人一区二区免费高清观看| 一个人看的www免费观看视频| 岛国在线免费视频观看| 伦理电影大哥的女人| 日韩成人av中文字幕在线观看 | 在线观看午夜福利视频| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 日本黄色视频三级网站网址| 一个人免费在线观看电影| 淫妇啪啪啪对白视频| 丰满的人妻完整版| 九九爱精品视频在线观看| 久久精品国产自在天天线| 欧美中文日本在线观看视频| www日本黄色视频网| 蜜桃亚洲精品一区二区三区| 精品久久久久久久久亚洲| 蜜桃亚洲精品一区二区三区| 亚洲在线观看片| 亚洲激情五月婷婷啪啪| 美女被艹到高潮喷水动态| 中文字幕人妻熟人妻熟丝袜美| 一级黄片播放器| 99久国产av精品| 久久久久久久久久成人| 国产精品乱码一区二三区的特点| 久久久精品大字幕| 久久久久久久亚洲中文字幕| 国产探花在线观看一区二区| 国产伦精品一区二区三区四那| 日日摸夜夜添夜夜添av毛片| 国产一级毛片七仙女欲春2| 夜夜爽天天搞| 黄色一级大片看看| 欧美激情在线99| 小蜜桃在线观看免费完整版高清| 99在线人妻在线中文字幕| 国产综合懂色| 久久精品国产亚洲av香蕉五月| 晚上一个人看的免费电影| 国产精品久久久久久亚洲av鲁大| 三级毛片av免费| 一区二区三区高清视频在线| 亚洲欧美中文字幕日韩二区| 欧美高清性xxxxhd video| 精品久久久噜噜| 免费看光身美女| 亚洲精品日韩av片在线观看| 舔av片在线| 日韩 亚洲 欧美在线| 日韩大尺度精品在线看网址| 日本黄色视频三级网站网址| 香蕉av资源在线| 欧美成人a在线观看| 99热这里只有精品一区| 国内精品一区二区在线观看| 国产男靠女视频免费网站| 淫秽高清视频在线观看| 在线免费十八禁| 国产在线精品亚洲第一网站| 亚洲人成网站在线播放欧美日韩| 国产精品av视频在线免费观看| 欧美激情在线99| 久久久久久久久久久丰满| 亚洲国产精品合色在线| 日韩欧美精品v在线| 欧美性猛交黑人性爽| 久久午夜亚洲精品久久| 1024手机看黄色片| 真实男女啪啪啪动态图| 你懂的网址亚洲精品在线观看 | 真实男女啪啪啪动态图| 国产视频一区二区在线看| www.色视频.com| av卡一久久| 最近最新中文字幕大全电影3| 亚洲成人av在线免费| 九九久久精品国产亚洲av麻豆| 69人妻影院| 干丝袜人妻中文字幕| 久久久欧美国产精品| 人妻夜夜爽99麻豆av| www.色视频.com| 最近中文字幕高清免费大全6| 欧美色欧美亚洲另类二区| 日韩国内少妇激情av| 国产伦精品一区二区三区视频9| 午夜福利在线观看吧| 久久久精品大字幕| 一进一出抽搐gif免费好疼| 有码 亚洲区| 最近的中文字幕免费完整| 波多野结衣巨乳人妻| 亚洲精品日韩在线中文字幕 | 亚洲无线在线观看| 日本黄色片子视频| 少妇熟女aⅴ在线视频| 久久精品夜夜夜夜夜久久蜜豆| 一区二区三区免费毛片| 成人美女网站在线观看视频| ponron亚洲| 嫩草影院精品99| 极品教师在线视频| 国产毛片a区久久久久| 少妇的逼好多水| 欧美3d第一页| 亚洲精品一卡2卡三卡4卡5卡| 国产一区亚洲一区在线观看| 亚洲av免费高清在线观看| 亚洲丝袜综合中文字幕| 亚洲乱码一区二区免费版| 天天躁日日操中文字幕| 日本一本二区三区精品| 男女啪啪激烈高潮av片| 深爱激情五月婷婷| 成年av动漫网址| 欧美日韩一区二区视频在线观看视频在线 | 成人美女网站在线观看视频| 欧美区成人在线视频| 国产精品一及| 可以在线观看的亚洲视频| 少妇的逼水好多| 蜜臀久久99精品久久宅男| 日本a在线网址| 校园春色视频在线观看| 久久精品91蜜桃| 男女边吃奶边做爰视频| 色av中文字幕| 亚洲中文字幕日韩| 成人特级av手机在线观看| 又爽又黄a免费视频| 日韩欧美精品v在线| 少妇熟女欧美另类| 国产乱人视频| 亚洲精品国产av成人精品 | 日韩欧美免费精品| 免费观看精品视频网站| 久久婷婷人人爽人人干人人爱| 亚洲一区高清亚洲精品| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站| 大香蕉久久网| 91在线精品国自产拍蜜月| 日韩欧美在线乱码| 亚洲熟妇中文字幕五十中出| 国产乱人偷精品视频| www日本黄色视频网| 久久精品国产亚洲网站| 国产女主播在线喷水免费视频网站 | 亚洲欧美日韩卡通动漫| 国产亚洲精品综合一区在线观看| 波多野结衣高清作品| 亚洲国产精品国产精品| 国产一区二区三区av在线 | 午夜福利成人在线免费观看| 99久国产av精品国产电影| 中文亚洲av片在线观看爽| 91久久精品国产一区二区成人| 亚洲成av人片在线播放无| 长腿黑丝高跟| 精品久久久久久久久久免费视频| 丰满的人妻完整版| 午夜精品国产一区二区电影 | 国产 一区 欧美 日韩| 一级a爱片免费观看的视频| aaaaa片日本免费| 男女啪啪激烈高潮av片| 欧美xxxx性猛交bbbb| 12—13女人毛片做爰片一| 毛片女人毛片| 中文亚洲av片在线观看爽| 成人无遮挡网站| 欧美日韩综合久久久久久| 国产淫片久久久久久久久| 国产真实伦视频高清在线观看| 亚洲国产精品久久男人天堂| 99在线人妻在线中文字幕| 国产探花极品一区二区| 一本久久中文字幕| 久久久午夜欧美精品| 精品99又大又爽又粗少妇毛片| 亚洲av免费在线观看| 国产亚洲精品综合一区在线观看| 欧美xxxx黑人xx丫x性爽| 99久久久亚洲精品蜜臀av| 国产精品久久久久久久久免| 99热只有精品国产| 中文字幕av在线有码专区| 一个人观看的视频www高清免费观看| 国产女主播在线喷水免费视频网站 | 一进一出抽搐gif免费好疼| 如何舔出高潮| 久久精品影院6| 国产精品乱码一区二三区的特点| 久久午夜亚洲精品久久| 美女免费视频网站| 99热精品在线国产| 最近中文字幕高清免费大全6| 不卡一级毛片| 男人狂女人下面高潮的视频| 十八禁网站免费在线| 国产伦在线观看视频一区| 最后的刺客免费高清国语| av天堂在线播放| 免费黄网站久久成人精品| 国产精品日韩av在线免费观看| 日韩精品中文字幕看吧| 99在线视频只有这里精品首页| eeuss影院久久| 精品久久久噜噜| 91在线精品国自产拍蜜月| 国产免费一级a男人的天堂| 女人被狂操c到高潮| 日本精品一区二区三区蜜桃| 国产精品99久久久久久久久| 禁无遮挡网站| 欧美日韩在线观看h| 精品久久久久久久久久免费视频| 日本爱情动作片www.在线观看 | 久久99热6这里只有精品| 国产精品,欧美在线| 91在线观看av| 18禁裸乳无遮挡免费网站照片| 欧美又色又爽又黄视频| 国产aⅴ精品一区二区三区波| 噜噜噜噜噜久久久久久91| 波野结衣二区三区在线| 非洲黑人性xxxx精品又粗又长| 亚洲美女视频黄频| 日韩在线高清观看一区二区三区| 男插女下体视频免费在线播放| 人人妻人人澡欧美一区二区| 国产又黄又爽又无遮挡在线| 深爱激情五月婷婷| 身体一侧抽搐| 九九爱精品视频在线观看| 91久久精品国产一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 麻豆乱淫一区二区| 日韩三级伦理在线观看| 一个人看视频在线观看www免费| 久久久国产成人免费| 国产不卡一卡二| 人人妻人人看人人澡| 午夜免费激情av| 色在线成人网| 国产探花在线观看一区二区| av卡一久久| 亚洲真实伦在线观看| 精品久久久久久久久久免费视频| 欧美激情在线99| 亚洲精品日韩在线中文字幕 | 我的老师免费观看完整版| 如何舔出高潮| 久久久精品94久久精品| 成人亚洲欧美一区二区av| 麻豆精品久久久久久蜜桃| 国产精品国产三级国产av玫瑰| 天天一区二区日本电影三级| 亚洲欧美中文字幕日韩二区| www.色视频.com| 搞女人的毛片| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 成人永久免费在线观看视频| 国产真实乱freesex| 草草在线视频免费看| 亚洲精华国产精华液的使用体验 | 青春草视频在线免费观看| 免费搜索国产男女视频| 波多野结衣高清无吗| 亚洲天堂国产精品一区在线| 少妇裸体淫交视频免费看高清| 欧美一区二区精品小视频在线| 欧美国产日韩亚洲一区| 观看美女的网站| 中文字幕av成人在线电影| 国产 一区 欧美 日韩| 久久6这里有精品| 蜜臀久久99精品久久宅男| 亚洲三级黄色毛片| 夜夜爽天天搞| 日韩欧美精品v在线| 精品久久久久久久人妻蜜臀av|