• <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)記
    亚洲成人国产一区在线观看| 亚洲在线自拍视频| 久久香蕉激情| 丝袜人妻中文字幕| 中文字幕人妻丝袜制服| 中文字幕人妻丝袜一区二区| 一个人免费在线观看的高清视频| 满18在线观看网站| 亚洲九九香蕉| 成人亚洲精品一区在线观看| 黑丝袜美女国产一区| 亚洲第一青青草原| 欧美黄色片欧美黄色片| 色94色欧美一区二区| 黑丝袜美女国产一区| 电影成人av| 夜夜爽天天搞| 日韩熟女老妇一区二区性免费视频| 精品卡一卡二卡四卡免费| 色播在线永久视频| 国产激情欧美一区二区| 久久久久国产一级毛片高清牌| 中国美女看黄片| 亚洲专区国产一区二区| 人人妻人人澡人人看| 国产色视频综合| 亚洲男人天堂网一区| 成人永久免费在线观看视频| 欧美中文综合在线视频| 别揉我奶头~嗯~啊~动态视频| 国产精品偷伦视频观看了| 久久婷婷成人综合色麻豆| 人妻丰满熟妇av一区二区三区 | 精品国产一区二区三区四区第35| 中文字幕精品免费在线观看视频| 一级片'在线观看视频| 国产91精品成人一区二区三区| 变态另类成人亚洲欧美熟女 | 午夜精品在线福利| 久久精品熟女亚洲av麻豆精品| 脱女人内裤的视频| 成年人黄色毛片网站| 19禁男女啪啪无遮挡网站| 亚洲国产欧美日韩在线播放| 一级毛片高清免费大全| x7x7x7水蜜桃| 亚洲男人天堂网一区| 悠悠久久av| 午夜91福利影院| 免费在线观看黄色视频的| 午夜福利一区二区在线看| 色婷婷久久久亚洲欧美| 国产97色在线日韩免费| 最新在线观看一区二区三区| 久久久久久久久久久久大奶| 99国产精品免费福利视频| 两个人免费观看高清视频| 一级a爱视频在线免费观看| 黄色毛片三级朝国网站| 亚洲av成人一区二区三| 嫁个100分男人电影在线观看| 欧美日韩亚洲综合一区二区三区_| 自线自在国产av| 日日夜夜操网爽| 亚洲人成电影免费在线| 成年人免费黄色播放视频| 国产成人系列免费观看| 中文字幕av电影在线播放| 丝袜人妻中文字幕| 99在线人妻在线中文字幕 | 亚洲中文av在线| 国产亚洲精品久久久久久毛片 | 久久人妻av系列| 最近最新免费中文字幕在线| 成年人免费黄色播放视频| 久久人人97超碰香蕉20202| 啦啦啦视频在线资源免费观看| 无人区码免费观看不卡| 丝袜美腿诱惑在线| 女同久久另类99精品国产91| 搡老乐熟女国产| 亚洲免费av在线视频| 1024视频免费在线观看| 午夜福利,免费看| 高潮久久久久久久久久久不卡| 麻豆av在线久日| 男人舔女人的私密视频| 男女下面插进去视频免费观看| 亚洲三区欧美一区| x7x7x7水蜜桃| 高清黄色对白视频在线免费看| 宅男免费午夜| 丁香欧美五月| 看片在线看免费视频| 我的亚洲天堂| 日韩一卡2卡3卡4卡2021年| 美女高潮到喷水免费观看| 免费少妇av软件| 国产欧美日韩一区二区三区在线| 亚洲 国产 在线| 欧美日韩瑟瑟在线播放| 亚洲中文日韩欧美视频| 美女扒开内裤让男人捅视频| 看免费av毛片| 国产高清国产精品国产三级| 成人国语在线视频| 国产欧美日韩一区二区精品| 国产精品秋霞免费鲁丝片| 一进一出抽搐动态| 18禁黄网站禁片午夜丰满| 黄色片一级片一级黄色片| 中文字幕人妻熟女乱码| 在线观看免费午夜福利视频| 亚洲人成伊人成综合网2020| 久久久国产一区二区| 女人精品久久久久毛片| 在线av久久热| 大码成人一级视频| 熟女少妇亚洲综合色aaa.| 久久香蕉激情| 久久久久久亚洲精品国产蜜桃av| 国产成人影院久久av| 国产成人影院久久av| 色精品久久人妻99蜜桃| 亚洲av熟女| 18禁国产床啪视频网站| 久久久国产欧美日韩av| 飞空精品影院首页| 黄网站色视频无遮挡免费观看| 美女国产高潮福利片在线看| 午夜91福利影院| 日日爽夜夜爽网站| 国产高清videossex| 欧美亚洲 丝袜 人妻 在线| 久久精品人人爽人人爽视色| 一区福利在线观看| 黑人巨大精品欧美一区二区蜜桃| 在线十欧美十亚洲十日本专区| 亚洲国产中文字幕在线视频| 久久久久视频综合| 国产成人免费观看mmmm| 国产精品久久视频播放| 国产色视频综合| 亚洲人成伊人成综合网2020| tube8黄色片| 国产精品久久久久久精品古装| 欧美 亚洲 国产 日韩一| 我的亚洲天堂| 亚洲avbb在线观看| e午夜精品久久久久久久| 免费在线观看黄色视频的| 99久久精品国产亚洲精品| 一级片免费观看大全| 成人手机av| 久久香蕉激情| www.精华液| 国精品久久久久久国模美| 啦啦啦视频在线资源免费观看| 国产成+人综合+亚洲专区| 王馨瑶露胸无遮挡在线观看| 国产乱人伦免费视频| 黑人巨大精品欧美一区二区mp4| 美女视频免费永久观看网站| 国产精品综合久久久久久久免费 | 又紧又爽又黄一区二区| 热re99久久精品国产66热6| 国产深夜福利视频在线观看| 国产一卡二卡三卡精品| av有码第一页| 如日韩欧美国产精品一区二区三区| 麻豆av在线久日| 久久香蕉激情| www.熟女人妻精品国产| √禁漫天堂资源中文www| 免费在线观看日本一区| 欧美性长视频在线观看| 在线观看免费视频日本深夜| ponron亚洲| 乱人伦中国视频| 巨乳人妻的诱惑在线观看| 久久精品熟女亚洲av麻豆精品| 男女床上黄色一级片免费看| 欧美日本中文国产一区发布| 99re6热这里在线精品视频| 国产精品自产拍在线观看55亚洲 | 国产成人啪精品午夜网站| 法律面前人人平等表现在哪些方面| 满18在线观看网站| 捣出白浆h1v1| 日韩人妻精品一区2区三区| 嫁个100分男人电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁裸乳无遮挡动漫免费视频| 桃红色精品国产亚洲av| 在线观看免费午夜福利视频| 精品欧美一区二区三区在线| 成年女人毛片免费观看观看9 | 超碰成人久久| 亚洲美女黄片视频| 中文字幕人妻丝袜一区二区| 精品乱码久久久久久99久播| 国产精品久久久av美女十八| 久久中文看片网| 一区二区日韩欧美中文字幕| 啦啦啦免费观看视频1| 啦啦啦免费观看视频1| 免费在线观看完整版高清| 搡老岳熟女国产| 999久久久精品免费观看国产| 999久久久国产精品视频| 久久精品国产a三级三级三级| 天堂俺去俺来也www色官网| 纯流量卡能插随身wifi吗| 巨乳人妻的诱惑在线观看| 国产精品久久久av美女十八| 美女午夜性视频免费| 亚洲第一青青草原| 日韩欧美免费精品| 国产精品久久久久久人妻精品电影| 麻豆av在线久日| 国产精品秋霞免费鲁丝片| 中出人妻视频一区二区| 最新的欧美精品一区二区| 18禁美女被吸乳视频| 最近最新免费中文字幕在线| 一级作爱视频免费观看| 成年女人毛片免费观看观看9 | 久久久久国产一级毛片高清牌| 最近最新中文字幕大全电影3 | 精品人妻在线不人妻| 国产欧美亚洲国产| 国产精品久久视频播放| 久久久精品国产亚洲av高清涩受| 午夜久久久在线观看| 精品一区二区三区av网在线观看| 亚洲国产精品合色在线| 欧美人与性动交α欧美软件| 男人舔女人的私密视频| 19禁男女啪啪无遮挡网站| 欧美日韩福利视频一区二区| 日日摸夜夜添夜夜添小说| 自线自在国产av| 人人妻人人澡人人爽人人夜夜| 香蕉国产在线看| 一级,二级,三级黄色视频| 国产一区二区三区在线臀色熟女 | 精品一区二区三区av网在线观看| 中文字幕制服av| 伊人久久大香线蕉亚洲五| 国产在线观看jvid| 一级黄色大片毛片| 男女之事视频高清在线观看| 99久久精品国产亚洲精品| 757午夜福利合集在线观看| 午夜日韩欧美国产| 老司机靠b影院| 丝袜在线中文字幕| 久久久久久久久免费视频了| 午夜日韩欧美国产| 老司机靠b影院| 十八禁网站免费在线| 亚洲熟妇熟女久久| 岛国毛片在线播放| 亚洲熟女精品中文字幕| 老司机福利观看| 亚洲色图综合在线观看| 亚洲成人手机| 首页视频小说图片口味搜索| 日韩精品免费视频一区二区三区| av在线播放免费不卡| 久久狼人影院| 久久久国产精品麻豆| 亚洲av欧美aⅴ国产| 精品免费久久久久久久清纯 | 亚洲熟女精品中文字幕| 99香蕉大伊视频| 别揉我奶头~嗯~啊~动态视频| 女人精品久久久久毛片| 亚洲在线自拍视频| 日韩免费高清中文字幕av| 曰老女人黄片| 国产欧美日韩一区二区三区在线| 国产真人三级小视频在线观看| 日韩欧美国产一区二区入口| 久久久久久久午夜电影 | 中文欧美无线码| 多毛熟女@视频| 成人影院久久| 天堂中文最新版在线下载| 久久精品成人免费网站| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添小说| 黑人巨大精品欧美一区二区mp4| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看| 亚洲精品久久午夜乱码| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| 久久久久精品国产欧美久久久| 国产精品香港三级国产av潘金莲| 在线永久观看黄色视频| 久久精品国产亚洲av高清一级| 法律面前人人平等表现在哪些方面| 亚洲午夜精品一区,二区,三区| 亚洲少妇的诱惑av| 欧美激情极品国产一区二区三区| 亚洲国产欧美一区二区综合| 悠悠久久av| 亚洲国产精品sss在线观看 | 757午夜福利合集在线观看| 国产亚洲精品一区二区www | 国产精品一区二区在线观看99| 男女免费视频国产| 很黄的视频免费| 午夜日韩欧美国产| 曰老女人黄片| av福利片在线| 成人三级做爰电影| 亚洲情色 制服丝袜| 国产精品欧美亚洲77777| 国产视频一区二区在线看| 18禁国产床啪视频网站| 欧美在线一区亚洲| 欧美精品人与动牲交sv欧美| 夜夜夜夜夜久久久久| 欧美最黄视频在线播放免费 | 国产高清激情床上av| 亚洲一区二区三区欧美精品| 99精品在免费线老司机午夜| 日韩制服丝袜自拍偷拍| 亚洲熟妇熟女久久| 黄色片一级片一级黄色片| 久久国产精品大桥未久av| 午夜福利一区二区在线看| 在线十欧美十亚洲十日本专区| 丁香欧美五月| 人人妻人人爽人人添夜夜欢视频| 精品视频人人做人人爽| 中文字幕人妻丝袜制服| 91大片在线观看| 国产亚洲精品第一综合不卡| 亚洲精品自拍成人| 精品欧美一区二区三区在线| 曰老女人黄片| www.999成人在线观看| 久99久视频精品免费| 午夜激情av网站| 欧美午夜高清在线| 美女 人体艺术 gogo| 欧美性长视频在线观看| 女人高潮潮喷娇喘18禁视频| 日韩人妻精品一区2区三区| 91精品国产国语对白视频| 水蜜桃什么品种好| 国产真人三级小视频在线观看| 欧美日韩亚洲综合一区二区三区_| 满18在线观看网站| 老司机福利观看| 黄色成人免费大全| 午夜福利欧美成人| 欧美精品高潮呻吟av久久| 999久久久国产精品视频| 日本精品一区二区三区蜜桃| 午夜精品在线福利| 久久久久久久午夜电影 | 色94色欧美一区二区| 99热国产这里只有精品6| 91大片在线观看| 国产97色在线日韩免费| 亚洲综合色网址| 免费在线观看影片大全网站| 成人亚洲精品一区在线观看| 亚洲欧洲精品一区二区精品久久久| 女人被躁到高潮嗷嗷叫费观| 99国产精品免费福利视频| 日本vs欧美在线观看视频| 午夜两性在线视频| 男女下面插进去视频免费观看| 在线观看免费高清a一片| 国产av一区二区精品久久| 在线av久久热| 欧美日韩亚洲高清精品| 黄频高清免费视频| 免费观看精品视频网站| 亚洲精品在线观看二区| av超薄肉色丝袜交足视频| 丰满迷人的少妇在线观看| 国产又爽黄色视频| 99香蕉大伊视频| 美女国产高潮福利片在线看| 午夜福利欧美成人| 老鸭窝网址在线观看| 精品久久久久久电影网| 正在播放国产对白刺激| 女性被躁到高潮视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美一区二区三区黑人| 国产成人av激情在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 大香蕉久久成人网| 91大片在线观看| 日韩免费av在线播放| 欧美精品一区二区免费开放| 91在线观看av| 亚洲久久久国产精品| 日本vs欧美在线观看视频| 一二三四在线观看免费中文在| 啦啦啦 在线观看视频| 亚洲第一欧美日韩一区二区三区| 女同久久另类99精品国产91| 在线观看免费高清a一片| 国产成人系列免费观看| 亚洲色图 男人天堂 中文字幕| 国产精品99久久99久久久不卡| 国产男女超爽视频在线观看| 淫妇啪啪啪对白视频| 久久亚洲精品不卡| 99热网站在线观看| 婷婷成人精品国产| 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看免费视频日本深夜| 宅男免费午夜| 99热网站在线观看| 国产亚洲精品久久久久久毛片 | 国产精品电影一区二区三区 | 久久婷婷成人综合色麻豆| 999久久久国产精品视频| 一个人免费在线观看的高清视频| 亚洲精品国产一区二区精华液| 一二三四在线观看免费中文在| 在线天堂中文资源库| 欧美黑人精品巨大| 国产色视频综合| 狠狠婷婷综合久久久久久88av| 高清在线国产一区| 亚洲成人免费av在线播放| 在线观看日韩欧美| 国产精品99久久99久久久不卡| 成年人免费黄色播放视频| 免费av中文字幕在线| 国产不卡一卡二| 最新在线观看一区二区三区| 涩涩av久久男人的天堂| 午夜老司机福利片| 成人av一区二区三区在线看| 91国产中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产精品自产拍在线观看55亚洲 | 亚洲成人免费av在线播放| 国产麻豆69| 激情视频va一区二区三区| 满18在线观看网站| 一级毛片精品| 在线免费观看的www视频| 老司机福利观看| 好男人电影高清在线观看| 又大又爽又粗| 国产在视频线精品| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区 | 高清视频免费观看一区二区| 亚洲精品国产一区二区精华液| 纯流量卡能插随身wifi吗| 久久精品成人免费网站| 91精品国产国语对白视频| 国产精品二区激情视频| 免费av中文字幕在线| 久久狼人影院| 国产亚洲欧美98| 18禁观看日本| 两性夫妻黄色片| 欧美成人午夜精品| 90打野战视频偷拍视频| 怎么达到女性高潮| 国产精品免费一区二区三区在线 | 精品国产亚洲在线| 交换朋友夫妻互换小说| 欧美人与性动交α欧美软件| 国产精品av久久久久免费| 欧美日韩国产mv在线观看视频| 国产三级黄色录像| 夜夜躁狠狠躁天天躁| 99精品欧美一区二区三区四区| 国产aⅴ精品一区二区三区波| 国产不卡一卡二| 国产激情久久老熟女| 欧美精品高潮呻吟av久久| 女人精品久久久久毛片| 久久婷婷成人综合色麻豆| 男女下面插进去视频免费观看| 亚洲欧美激情综合另类| 亚洲精品在线美女| av有码第一页| 国产激情欧美一区二区| 久久中文看片网| 久久热在线av| 久久精品aⅴ一区二区三区四区| 久久久久国产一级毛片高清牌| 动漫黄色视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲精品国产区一区二| 成人黄色视频免费在线看| 中文亚洲av片在线观看爽 | 日韩一卡2卡3卡4卡2021年| 天天影视国产精品| av网站在线播放免费| 精品一区二区三区视频在线观看免费 | 18禁观看日本| 国产成人av教育| 乱人伦中国视频| 婷婷丁香在线五月| 一本综合久久免费| av中文乱码字幕在线| 18在线观看网站| 最新美女视频免费是黄的| 午夜两性在线视频| 亚洲精华国产精华精| 黄色成人免费大全| a在线观看视频网站| 国产视频一区二区在线看| 久久香蕉激情| 99久久人妻综合| 两人在一起打扑克的视频| 亚洲欧美色中文字幕在线| 欧美中文综合在线视频| 在线永久观看黄色视频| 人人澡人人妻人| 精品亚洲成国产av| 成熟少妇高潮喷水视频| 久热这里只有精品99| 欧美成人免费av一区二区三区 | 亚洲精品中文字幕一二三四区| 欧美 亚洲 国产 日韩一| videosex国产| 国产成+人综合+亚洲专区| 免费高清在线观看日韩| 欧美精品啪啪一区二区三区| 久久久久久久久久久久大奶| 激情在线观看视频在线高清 | 欧美日韩瑟瑟在线播放| 欧美在线一区亚洲| 少妇 在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲全国av大片| aaaaa片日本免费| 在线观看一区二区三区激情| 精品国产一区二区久久| 欧美av亚洲av综合av国产av| 午夜久久久在线观看| 国产蜜桃级精品一区二区三区 | 国产免费av片在线观看野外av| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 国产99久久九九免费精品| 午夜福利在线免费观看网站| 成人精品一区二区免费| 国产有黄有色有爽视频| 午夜激情av网站| 国产精品成人在线| av线在线观看网站| 亚洲第一欧美日韩一区二区三区| 久久天堂一区二区三区四区| 日日夜夜操网爽| 欧美激情高清一区二区三区| 久久精品亚洲精品国产色婷小说| 高清黄色对白视频在线免费看| 女人爽到高潮嗷嗷叫在线视频| 午夜亚洲福利在线播放| 国产乱人伦免费视频| 一区二区三区精品91| 亚洲五月色婷婷综合| 一进一出抽搐动态| 久久人人97超碰香蕉20202| 亚洲成a人片在线一区二区| 国产精品久久久av美女十八| 免费av中文字幕在线| 久久精品aⅴ一区二区三区四区| 久久午夜综合久久蜜桃| 亚洲第一青青草原| 久久午夜综合久久蜜桃| 19禁男女啪啪无遮挡网站| 香蕉久久夜色| 亚洲五月天丁香| 国产一区二区三区视频了| 国产成人一区二区三区免费视频网站| 国产亚洲精品第一综合不卡| 黄色丝袜av网址大全| 啦啦啦 在线观看视频| 成人永久免费在线观看视频| 色精品久久人妻99蜜桃| 热re99久久精品国产66热6| 日本撒尿小便嘘嘘汇集6| 国产日韩欧美亚洲二区| 欧美国产精品一级二级三级| 亚洲av美国av| av福利片在线| videos熟女内射| 亚洲精品国产一区二区精华液| 亚洲精品美女久久久久99蜜臀| 亚洲成国产人片在线观看| 精品亚洲成国产av| 久久久久国产一级毛片高清牌| 精品视频人人做人人爽| 人人妻人人添人人爽欧美一区卜| 天堂√8在线中文| 一本大道久久a久久精品| 精品久久久久久久毛片微露脸| 香蕉久久夜色| 99国产精品一区二区蜜桃av | 狠狠狠狠99中文字幕| 黄色成人免费大全| av福利片在线| 国产淫语在线视频| 窝窝影院91人妻|