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

    海南島橡膠林生態(tài)系統(tǒng)碳通量及其影響因子研究

    2014-04-29 00:44:03吳志祥等
    熱帶作物學(xué)報 2014年11期
    關(guān)鍵詞:影響因子海南島

    吳志祥等

    摘 要 以位于海南島西北部儋州地區(qū)的橡膠林為研究對象,采用渦度相關(guān)技術(shù)進行長期CO2通量觀測,獲得海南島橡膠林生態(tài)系統(tǒng)碳通量的動態(tài)特征,并分析其對各影響因子的響應(yīng)情況。結(jié)果表明:(1)橡膠林生態(tài)系統(tǒng)凈碳交換量(NEE)日動態(tài)均為“U”型曲線,白天均為碳吸收(NEE為負(fù)值),夜間均為碳排放(NEE為正值);全年各月除2月外,NEE月總值均為負(fù)值,吸收碳,雨季4~10月各月NEE值均較大,旱季1~3月、11~12月,各月NEE值均較??;橡膠林年固碳能力強大,超過11.0 t C/(hm2·a)。(2)橡膠林生態(tài)系統(tǒng)NEE與光合有效輻射(PAR)、冠層內(nèi)大氣平均溫度(Ta)、地下5 cm土壤溫度(T_-5 cm)和飽和水汽壓差(VPD)及地下5 cm土壤含水量(VWC_-5 cm)等均負(fù)相關(guān),但影響強度不同,其中PAR影響最為顯著;另外NEE還與橡膠林本身林分狀況及周圍大氣環(huán)境CO2濃度等相關(guān)。

    關(guān)鍵詞 碳通量;影響因子;橡膠林生態(tài)系統(tǒng);海南島

    中圖分類號 P422.4;Q948;S718 文獻標(biāo)識碼 A

    Abstract The rubber plantation in the northwest of Danzhou,Hainan area was used as the study object with the eddy covariance technique for long-term CO2 flux observations, Hainan rubber forest ecosystem carbon flux dynamic characteristics were obtained, and the responses on the environmental factors were analyzed in the research. The results showed as follows:(1)The daily dynamic of net ecosystem carbon exchange(NEE)of rubber plantation ecosystem was an“U”shaped curve; and rubber trees absorbed carbon during the daytime(NEE negative), released carbon nighttime(NEE positive). All the monthly NEE values throughout the year except February were negative(rubber trees absorb carbon); and the NEE values from April to October during the rainy season were large, from January to March, November to December during the dry season little. Rubber plantation ecosystem had the strong carbon sequestration capability, and the yearly NEE values were more than 11.0 t C/(hm2·a).(2)The NEE of rubber plantation ecosystem and all the environmental factors such as photosynthetically active radiation(PAR), the average atmosphere temperature within the canopy(Ta), -5 cm soil temperature(T_-5 cm), vapor pressure deficit(VPD)and -5 cm soil volume water content(VWC_-5 cm)were negatively correlated with different intensity, and PAR was the most significantly related factor. NEE was also associated with the rubber plantation itself stand conditions and ambient atmospheric CO2 concentration. The results could provide a reference for tropical plantation ecosystems carbon estimates and offer basic data reduction for the government trade negotiations.

    Key words Carbon fluxes; Factors; Rubber plantation ecosystem; Hainan Island

    doi 10.3969/j.issn.1000-2561.2014.11.001

    巴西橡膠樹(Hevea brasiliensis Müll. Arg.)人工林(本文簡稱為橡膠林)是我國熱帶地區(qū)最重要的森林生態(tài)系統(tǒng),在我國海南、云南、廣東等地的種植面積已超過1×106 hm2。多年種植橡膠樹的生產(chǎn)實踐和研究證明,橡膠人工林是世界上熱帶地區(qū)開發(fā)、在旱地建立的最好的生態(tài)系統(tǒng)之一,在我國熱區(qū)及東南亞種植橡膠樹,是應(yīng)用生態(tài)學(xué)原理開發(fā)熱區(qū)獲得成功的一個范例[1]。在熱帶次生林、稀樹草原和草地灌叢等基礎(chǔ)上開發(fā)種植的橡膠林(林下間種綠肥覆蓋),其結(jié)構(gòu)與熱帶雨林、熱帶季風(fēng)性雨林相近似[2-4]。但另一方面,因近10年來天然橡膠干膠價格上漲,橡膠樹種植面積擴大很快,尤其在云南西雙版納更是擴大更快,導(dǎo)致一系列負(fù)面影響,諸如生物多樣性降低、原始次生林和水源林的破壞、土壤碳儲量的減低等爭議,在生態(tài)學(xué)界對熱區(qū)種植橡膠樹產(chǎn)生反感、抵觸甚至反對的聲音[5-9]。

    橡膠林生態(tài)方面研究,前人進行了大量的工作。尤其近年來,許多研究者從橡膠林生態(tài)系統(tǒng)服務(wù)功能[2]、土壤碳固存[3-4]、橡膠林水分平衡[5,7]、膠林生物多樣性[8]及原始次生林會遭到破壞[9]等許多方面進行了相關(guān)研究,但是許多研究還是從短期數(shù)據(jù)獲得的結(jié)論,或者還處于初步階段,橡膠林的生態(tài)問題仍未透徹研究。為了對橡膠林生態(tài)系統(tǒng)進行深入研究,2009年底建成農(nóng)業(yè)部儋州熱帶作物科學(xué)觀測實驗站(簡稱儋州實驗站)50 m通量觀測塔,采用渦度相關(guān)方法對橡膠林進行長期定位觀測,以期能為橡膠林的生態(tài)功能研究提供支持。本研究應(yīng)用儋州實驗站2010~2013年的CO2通量定位觀測數(shù)據(jù),估算凈生態(tài)系統(tǒng)CO2交換量(Net Ecosystem Exchange, NEE)并分析其季節(jié)特征,從固碳的角度量化分析橡膠林的生態(tài)效益,科學(xué)評價橡膠林的碳匯效益,為政府決策和碳匯貿(mào)易提供基礎(chǔ)數(shù)據(jù);同時可為橡膠林經(jīng)營管理提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)自然概況

    研究地區(qū)位于海南島西北部,試驗地位于農(nóng)業(yè)部儋州熱帶作物科學(xué)觀測實驗站內(nèi)(19°32′47″N,109°28′30″E),距儋州市區(qū)約15 km,地形相對平坦,高差不超過20 m。樣地所在地為第二代膠園,2001年種植,2009年正式開始產(chǎn)膠,橡膠樹栽培品系為熱研7-33-97。膠林群落林冠平均高度13.0 m,結(jié)構(gòu)單一,分層明顯,上層為橡膠林喬木層,高度為11~14 m,下層為林下草本層,高度為0.4 m左右。試驗區(qū)氣候為熱帶海島季風(fēng)氣候,一年分明顯干濕兩季,每年5月至10月是濕季,11月至第二年4月是干季。年均溫21.5~28.5 ℃,歷年平均降雨量為1 607 mm,主要分布在7、8、9三個月,占全年降雨量的72%以上(以上為中國熱帶農(nóng)業(yè)科學(xué)院氣象站1954~2010年數(shù)據(jù))。研究年份(2010年)氣溫降水月動態(tài)如圖1所示。土壤為花崗巖風(fēng)化而成的磚紅壤,多為砂質(zhì)粘壤土。

    1.2 通量觀測方法

    橡膠林中50 m高觀測鐵塔安裝的儀器較多,與本研究相關(guān)的包括2部分:開路渦度相協(xié)方差系統(tǒng)(Open Path Eddy Covariance, OPEC)和常規(guī)氣象觀測系統(tǒng)(Routine Meteorological System, RMET)。安裝情況如下:OPEC安裝高度25 m,系統(tǒng)由三維超聲風(fēng)速儀CSAT-3(Campbell Scientific Inc, USA)、開路紅外CO2/H2O氣體分析儀Li-7500(Li-Cor, USA)和數(shù)據(jù)采集器CR3000(Campbell Scientific Inc, USA)組成,CSAT3和LI-7500采樣頻率均為10 Hz。RMET由分別安裝在鐵塔上7層和地下3層梯度觀測系統(tǒng)收集風(fēng)速、溫濕度、土壤熱通量和土壤水分等相關(guān)數(shù)據(jù)組成,有2個CR3000數(shù)據(jù)采集器收集、存儲數(shù)據(jù)。原始數(shù)據(jù)記錄和實時計算由Campbell公司提供的系統(tǒng)軟件Loggernet(Campbell)在線完成。鐵塔其他儀器詳細(xì)安裝詳見文獻[10]。

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

    為使通量觀測能真實反映橡膠林與大氣間的交換過程,首先要對觀測數(shù)據(jù)進行嚴(yán)格的篩選,方法是:對每個觀測日0 : 00~23 : 30每隔0.5 h共48個CO2通量觀測值用Excel作圖, 對數(shù)據(jù)的連續(xù)性和質(zhì)量進行判斷,剔除一些異?;蚱x較大的數(shù)據(jù)(“去除野點”,標(biāo)準(zhǔn)采用如±3δ);然后采用平均晝夜變化法(Mean Diurnal Variation, MDV)[11-12]的數(shù)據(jù)處理方法對異常數(shù)據(jù)進行替換。用MDV方法插補缺失或被剔除的異常數(shù)據(jù)時,異?;蛉笔?shù)據(jù)用鄰近7~10 d的相同時段觀測的正常數(shù)據(jù)平均值來替換,如此進行可得到一套完整可靠的觀測數(shù)據(jù)。

    2 結(jié)果與分析

    2.1 橡膠林生態(tài)系統(tǒng)碳通量變化特征

    2.1.1 NEE日動態(tài)特征 圖2給出橡膠林生態(tài)系統(tǒng)不同季節(jié)NEE平均日動態(tài),不同季節(jié)平均日動態(tài)為當(dāng)季各月每天半小時觀測值的平均值。整體來看,橡膠林所有季節(jié)NEE日動態(tài)均為“U”型曲線;NEE白天均為碳吸收(即NEE為負(fù)值),夜間均為碳排放(即NEE為正值,本文討論NEE值均指其絕對值)。橡膠林雨季(4~10月)碳吸收時段為7 : 30~19 : 00,年初旱季(1~3月)碳吸收時段為8 : 30~18 : 30,年末旱季(11~12月)碳吸收時段為8 : 00~18 : 30。NEE雨季峰值出現(xiàn)在12 : 30,橡膠林最強碳吸收為-0.856 0 mg CO2/(m2·s);NEE年初旱季峰值出現(xiàn)在14 : 00,橡膠林本季最強碳吸收為

    -0.400 5 mg CO2/(m2·s);年末旱季峰值出現(xiàn)在11 : 00,本季最強碳吸收為-0.559 3 mg CO2/(m2·s)。

    從2010年全年各月的半小時日平均動態(tài)(圖略)來看,橡膠林生態(tài)系統(tǒng)NEE最大值出現(xiàn)在8月份的13 : 00,最強碳吸收為-1.051 9 mg CO2/(m2·s)。全年NEE月均最大值為6月份,碳吸收為-0.225 8 mg CO2/(m2·s);全年僅2月份NEE為正值,碳排放為0.021 0 mg CO2/(m2·s),其余月份均為碳吸收。

    2.1.2 NEE月動態(tài)特征 海南島橡膠林生態(tài)系統(tǒng)2010年全年凈生態(tài)系統(tǒng)碳交換量各月總值如圖3所示。全年除2月份為碳排放(NEE為正值),其余月份均為碳吸收(NEE為負(fù)值)。雨季4~10月各月NEE值均較大,分別為-91.26、-144.94、-159.62、-149.91、-158.20、-148.98、-76.702 g C/(m2·mon),這與雨季光合有效輻射較強、氣溫較高、降水較多,橡膠林生長旺盛有關(guān)。旱季1~3月、11~12月,各月NEE值均較小,分別為-35.0、13.85、-40.57、-82.50、-59.60 g C/(m2·mon),這與當(dāng)季光合有效輻射較弱、氣溫降低、降水減少,橡膠林落葉、生長減弱等有關(guān)。但仔細(xì)研究會發(fā)現(xiàn),2010年當(dāng)年10月NEE比11月稍低,這可能與當(dāng)年10月份多陰雨天氣、PAR較弱造成光合作用減弱而導(dǎo)致。

    2010年度橡膠林生態(tài)系統(tǒng)碳吸收最強月份為6月,達到-159.62 g C/(m2·mon),碳排放最強(吸收最弱)月份為2月,為13.85 g C/(m2·mon)。2013年度(2012-09-01~2013-08-31)橡膠林生態(tài)系統(tǒng)最強碳吸收月份為7月,達到-168.56 g C/(m2·mon),碳排放最強月份為2月,為18.56 g C/(m2·mon)。

    2.1.3 NEE年變化 把全年各月NEE值加總,就可得全年NEE值,即全年總碳吸收值。2010年海南島橡膠林生態(tài)系統(tǒng)年總NEE為-1 133.45 g C/(m2·a)[相當(dāng)于-11.33 t C/(hm2·a)],2012年9月至2013年8月,NEE為-1 087.58 g C/(m2·a)[相當(dāng)于-10.88 t C/(hm2·a)],平均為1 110.52 g C/(m2·a)[相當(dāng)于-11.10 t C/(hm2·a)]。

    橡膠林為熱帶典型人工林,研究地點林齡為10齡,生長是最為旺盛時候,因此它吸收二氧化碳能力十分強大。2010年橡膠林生態(tài)系統(tǒng)NEE年總量為-1 133.45 g C/(m2·a),高于亞熱帶的楊樹人工林-579.0 g C/(m2·a)[24],高于亞熱帶和熱帶的紙漿林(桉樹林)[25-26],也高于位置比較接近的海南島尖峰嶺熱帶山地雨林平均值(-236±42)g C/(m2·a)[27]。

    2.2 橡膠林凈生態(tài)系統(tǒng)交換量影響因素分析

    2.2.1 光合有效輻射對橡膠林生態(tài)系統(tǒng)NEE的影響

    圖4給出海南島橡膠林生態(tài)系統(tǒng)半小時尺度上白天NEE對光合有效輻射(Photosynthetically Active Radiation, PAR)的響應(yīng)情況。白天NEE與PAR呈負(fù)相關(guān),相關(guān)系數(shù)達0.656 4。在半小時尺度上,NEE的最主要影響因子是PAR,這與葉片光合作用相一致,短時間內(nèi)PAR影響生態(tài)系統(tǒng)光合作用。也就是說,橡膠林生態(tài)系統(tǒng)白天NEE對PAR的響應(yīng)實際就是橡膠樹葉片本身光合作用對PAR的響應(yīng)(只不過存在眾多葉片的加和作用而已)。

    采用式(3)Michaelis-Menten模型的方程,對海南島橡膠林生態(tài)系統(tǒng)2010年全年各月白天NEE對PAR進行擬合,R2除1月、2月份外,其余月份在0.187~0.481間,并且達到極顯著水平(表1)。因1、2月回歸擬合較差,不在后續(xù)討論之列。

    雨季橡膠林生長旺盛,從4~10月,初始光能利用率α和生態(tài)系統(tǒng)呼吸速率Rd呈遞減趨勢,α從0.005 136 mg CO2·μmol·photon-1遞減到0.000 275 mg CO2·μmol·photon-1,Rd從0.160 734 mg CO2/(m2·s)遞減到0.107 810 mg CO2/(m2·s);最大光合速率Pmax從5月份的4.596 569 mg CO2/(m2·s)遞減到10月份的1.646 265 mg CO2/(m2·s)(也有個別月份變化趨勢不一致)。

    利用Michaelis-Menten模型擬合的橡膠林生態(tài)系統(tǒng)初始光能利用率α、最大光合作用速率Pmax、生態(tài)系統(tǒng)呼吸速率Rd最大值分別是0.005 136 mg CO2·μmol·photon-1(4月)、4.596 569 mg CO2/(m2·s)(5月)和0.190 734 mg CO2/(m2·s)(4月),大于長白山溫帶針闊混交林光合作用相關(guān)參數(shù)[0.004 1 mg CO2·μmol·photon-1、1.40 mg CO2/(m2·s)、0.34 mg CO2/(m2·s)][28],更大于黑龍江寒溫帶針葉林光合作用相關(guān)參數(shù)[0.003 2 mg CO2·μmol·photon-1、0.855 4 mg CO2/(m2·s)、0.273 4 mg CO2/(m2·s)][29],說明熱帶橡膠林生態(tài)系統(tǒng)比溫帶針闊混交林生態(tài)系統(tǒng)光合與呼吸作用活躍,遠(yuǎn)比寒溫帶針葉林生態(tài)系統(tǒng)活躍。

    不僅如此,橡膠林生態(tài)系統(tǒng)光合作用相關(guān)參數(shù)(初始光能利用率α、最大光合作用速率Pmax和生態(tài)系統(tǒng)呼吸速率Rd)最大值出現(xiàn)的月份較早,分別為4月、5月和4月,而溫帶針闊混交林均出現(xiàn)在6月,到北方寒溫帶針葉林則分別出現(xiàn)在6月、7月和6月。也就是,越往低緯度,光合作用參數(shù)越早達到最大值;越往北方,光合作用參數(shù)越晚達到最大值。這實際上也證明了植物光合作用與環(huán)境因素關(guān)系密切,尤其緯度因素影響的PAR,并與氣溫和降水的協(xié)同影響所致。

    2.2.2 溫度因子對橡膠林生態(tài)系統(tǒng)NEE的影響

    圖5給出橡膠林NEE對冠層內(nèi)大氣平均溫度(Ta)的響應(yīng)情況,白天NEE與大氣平均溫度呈負(fù)相關(guān),相關(guān)系數(shù)R為0.470 9,相對PAR較小。對于海南島橡膠林生態(tài)系統(tǒng)而言,因地處熱帶,終年氣溫較高,因此系統(tǒng)NEE受溫度影響不明顯。

    圖6給出橡膠林生態(tài)系統(tǒng)白天NEE對地表以下5 cm土壤溫度(T_-5 cm)的響應(yīng)情況,白天NEE與土壤溫度呈負(fù)相關(guān),相關(guān)系數(shù)R為0.355 4,相對PAR和Ta較小。對于海南島橡膠林生態(tài)系統(tǒng)而言,因地處熱帶,常年氣溫較高,土壤溫度也較高,生態(tài)系統(tǒng)呼吸作用終年旺盛,因此系統(tǒng)NEE受土壤溫度影響不明顯。

    2.2.3 水分因子對橡膠林生態(tài)系統(tǒng)NEE的影響

    橡膠樹屬于熱帶森林樹種,其生長發(fā)育需大量水分。海南島橡膠林生態(tài)系統(tǒng)地處熱帶海島季風(fēng)氣候,全年降水充沛,但有明顯的旱雨季節(jié),所以水分因素一直是影響橡膠林生態(tài)系統(tǒng)的重要因素,盡管可能水分在短期或一定時期內(nèi)與NEE的相關(guān)性不顯著。本節(jié)探討水分對橡膠林生態(tài)系統(tǒng)NEE的影響主要從3個因子入手:飽和水汽壓差(Vapor Pressure Deficit, VPD)、地表以下5 cm土壤含水量(VWC_-5 cm)和月降水量(Pr)。

    圖7給出橡膠林生態(tài)系統(tǒng)白天NEE對冠層飽和水汽壓差的響應(yīng)情況,白天NEE與VPD呈負(fù)相關(guān),相關(guān)系數(shù)R為0.418 9,相對PAR和Ta較小,但相對T_-5 cm較大。對于海南島橡膠林生態(tài)系統(tǒng)而言,因一年分為旱雨兩季,大氣中水汽多少(在當(dāng)時氣溫下,飽和與否)會影響橡膠林生態(tài)系統(tǒng)的光合作用和呼吸作用,進而影響NEE。但總的來說,橡膠林生態(tài)系統(tǒng)水汽虧缺還是不大,并沒有成為其限制因子(至少在研究的2010年和2012~2013年如此)。實際上,水汽對NEE的作用有兩方面,過少或過多均影響整個系統(tǒng)碳吸收或排放作用。本研究中在不同月份可能VPD相差不大(因為VPD是大氣溫度和大氣濕度的復(fù)合變量,2個月大氣溫度和大氣濕度不盡相同),但橡膠林生態(tài)系統(tǒng)NEE可能相差甚大,說明VPD不是影響NEE的關(guān)鍵因素。

    圖8給出橡膠林生態(tài)系統(tǒng)白天NEE對地表以下5 cm土壤含水量(VWC_-5 cm)的響應(yīng)情況,白天NEE與土壤含水量呈負(fù)相關(guān),相關(guān)系數(shù)R為0.370 4,相對PAR、Ta和VPD均較小。對于海南島橡膠林生態(tài)系統(tǒng)而言,一年分旱雨季節(jié),土壤含水量也呈明顯季節(jié)變化。土壤含水量不僅會影響生態(tài)系統(tǒng)光合作用,更會影響系統(tǒng)呼吸作用,因此系統(tǒng)NEE受土壤含水量的影響。但是,研究年份全年降水較為充沛,土壤含水量沒有成為NEE的限制因子。

    比較分析圖7和圖8,水分并不是海南島橡膠林生態(tài)系統(tǒng)NEE的限制因子。對于森林生態(tài)系統(tǒng)而言,大多數(shù)研究均表明,水分并不是其NEE的限制因子[26,30],尤其在生長季節(jié)(或者在海南島橡膠林生態(tài)系統(tǒng)稱為雨季)更是如此。

    上述半小時尺度分析結(jié)果表明,VPD和地表5 cm以下土壤含水量并不是NEE的限制因子,二者雖有相關(guān)性,但并不顯著。這可能與較小的時間尺度(瞬間尺度)有關(guān)。如果把時間間隔加大到月尺度(實際也就相當(dāng)于把整年數(shù)據(jù)劃分12個數(shù)據(jù)子集),可能效果不同。

    圖9給出橡膠林生態(tài)系統(tǒng)各月NEE對各月降水量的響應(yīng)情況,各月NEE與當(dāng)月降水量呈負(fù)相關(guān),相關(guān)系數(shù)R為0.882 4,相關(guān)性極為顯著。對于海南島橡膠林生態(tài)系統(tǒng)而言,屬于熱帶海島季風(fēng)氣候,全年降水在1 700 mm左右。一年分旱雨兩季,旱季并不是橡膠林旺盛生長季節(jié),其光合作用等系統(tǒng)生化過程受到制約;雨季(此時PAR強烈、氣溫也高)則是其生長旺季,凈生態(tài)系統(tǒng)碳交換量增加迅速。當(dāng)然,這既可說是雨季降水的促進,實際也是當(dāng)季PAR、氣溫等的促進作用,或者說是各個環(huán)境因子的共同作用結(jié)果。

    2.2.4 葉面積指數(shù)和CO2濃度對橡膠林生態(tài)系統(tǒng)NEE的影響 對陸地生態(tài)系統(tǒng)而言,系統(tǒng)NEE的影響因素除環(huán)境因素外,植物本身的生理生態(tài)學(xué)特性或過程對其NEE也具有重要影響。林分葉面積指數(shù)(Leaf Area Index, LAI)是影響生態(tài)系統(tǒng)NEE的重要指標(biāo)。

    對橡膠林生態(tài)系統(tǒng)而言,橡膠林林分LAI的增加對系統(tǒng)NEE的影響包括2個方面:一是通過影響橡膠林本身光合作用,在季節(jié)尺度上制約生態(tài)系統(tǒng)光合生產(chǎn)力(Gross Ecosystem Production, GEP),從而影響了NEE;二是橡膠林LAI的增加,本身就是增加了森林的葉量,積累生物量,當(dāng)然導(dǎo)致二氧化碳的固存。因此,林分LAI的增加,會增加生態(tài)系統(tǒng)的NEE,二者呈負(fù)相關(guān)關(guān)系。

    圖10給出了橡膠林生態(tài)系統(tǒng)2010年各月LAI均值與各月NEE總量的對應(yīng)變化趨勢。橡膠林生態(tài)系統(tǒng)2月份為落葉期,LAI達到當(dāng)年最小值,接近1;而此時NEE達到最大值,為正。隨著3月份抽葉,LAI急劇升高,到4、5月后,LAI升幅變緩;而此段時間NEE也類似變化,NEE急劇下降,至5月份后降幅變緩。然后6~10月,LAI維持高位;NEE也多維持當(dāng)年低值附近,然而在10月份,因葉片光合能力降低,再加上2010年當(dāng)月多陰雨天氣,此月NEE降幅明顯。此后,從11、12月到次年1月,隨著林分LAI降低,NEE逐漸升高。周而復(fù)始。在橡膠林整年生長過程中,橡膠林林分LAI較好地解釋了生態(tài)系統(tǒng)NEE的季節(jié)變化特征。

    圖11回歸分析了海南島橡膠林生態(tài)系統(tǒng)2010年全年各月NEE對林分LAI的響應(yīng)情況。從月尺度看,橡膠林生態(tài)系統(tǒng)NEE與LAI呈顯著的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)為0.796 2。

    陸地生態(tài)系統(tǒng)植被生長環(huán)境CO2濃度也會影響生態(tài)系統(tǒng)NEE。實際上空氣CO2濃度與橡膠林生態(tài)系統(tǒng)NEE影響過程比較復(fù)雜,它對橡膠林光合作用和呼吸作用均會產(chǎn)生重大影響。本研究選取不同旱雨季節(jié)半小時數(shù)據(jù)平均,得到CO2濃度與NEE日動態(tài)關(guān)系(圖12)。無論旱季或雨季,NEE白天均為負(fù)值,此時空氣CO2濃度也相對較低;NEE夜間均為正值,此時空氣CO2濃度相對較高;在晝夜交替時段,空氣CO2濃度則是處于相對濃度較高值。這樣的變化趨勢比較容易解釋橡膠林白天作為碳匯,固定較多的CO2導(dǎo)致白天濃度降低;夜間呼吸釋放CO2,導(dǎo)致CO2濃度會略有升高;晝夜交替時段,因氣溫升高,植被和土壤呼吸均快速升高,但植被光合作用固定CO2還較小,因此大氣CO2處于較高峰。從變化趨勢來看,旱季和雨季沒有顯著性差異。

    2.2.5 橡膠林生態(tài)系統(tǒng)NEE影響因素綜合分析

    為綜合分析各個環(huán)境因子對橡膠林生態(tài)系統(tǒng)NEE的影響,以NEE為因變量,以各個環(huán)境光合有效輻射PAR、冠層內(nèi)大氣平均溫度Ta、飽和水汽壓差VPD、地表以下5 cm和20 cm土壤溫濕度共7個因子為自變量,進行多因子逐步回歸分析,可建立環(huán)境因子與NEE在不同水平上的回歸方程。本研究選取99%置信水平的方程,回歸方程剩下3個因子PAR、VPD和T_-5。回歸方程獲得NEE與上述3個因子的模型如下:

    NEE=-0.001PAR+0.096VPD-0.10T_-5-0.243(5)

    上式?jīng)Q定系數(shù)R2=0.357,略低,但相關(guān)性極為顯著?;貧w模型顯示,橡膠林生態(tài)系統(tǒng)NEE與光合有效輻射PAR和地表以下5 cm土壤溫度呈負(fù)相關(guān),與飽和水汽壓差VPD呈正相關(guān),結(jié)果與前面分析基本一致。但是,前面分析中的冠層內(nèi)大氣平均溫度Ta因子,在回歸方程沒有出現(xiàn),分析其主要原因,主要是Ta與T_-5相關(guān)性很大(相關(guān)系數(shù)達0.88)、與PAR相關(guān)性也較大(相關(guān)系數(shù)達到0.60),因此在逐步回歸時,其效應(yīng)被PAR和T_-5替代了,在回歸方程中沒有出現(xiàn)。

    3 討論與結(jié)論

    3.1 橡膠林生態(tài)系統(tǒng)凈碳交換量特征

    因受PAR、氣溫及橡膠樹本身等的影響,橡膠林所有季節(jié)NEE日動態(tài)均為“U”型曲線;NEE白天均為碳吸收(即NEE為負(fù)值),夜間均為碳排放(即NEE為正值)。橡膠林雨季(4~10月)碳吸收時段為7 : 30~19 : 00,年初旱季(1~3月)碳吸收時段為8 : 30~18 : 30,年末旱季(11~12月)碳吸收時段為8 : 00~18 : 30。

    海南島橡膠林生態(tài)系統(tǒng)全年除2月份為碳排放外,其余月份均為碳吸收。雨季4~10月各月NEE值均較大,旱季1~3月、11~12月,各月NEE值均較小。因為旱雨季節(jié),會有光合有效輻射強弱、氣溫降水多少,橡膠林生長旺盛與否的差異。2010年度橡膠林生態(tài)系統(tǒng)碳吸收最強月份為6月,達到-159.62 g C/(m2·mon),碳排放最強(吸收最弱)的月份為2月,為13.85 g C/(m2·mon)。

    橡膠林為典型熱帶人工林,研究樣地樹齡較輕,固碳能力強大。2010年橡膠林生態(tài)系統(tǒng)NEE年總量為-1 133.45 g C/(m2·a)[相當(dāng)于-11.33 t C/(hm2·a)],2012年9月至2013年8月,NEE為-1 087.58 g C/(m2·a)[相當(dāng)于-10.88 t C/(hm2·a)],平均為1 110.52 g C/(m2·a)[相當(dāng)于-11.10 t C/(hm2·a)]。碳匯能力非常強大,高于亞熱帶的楊樹人工林,高于亞熱帶和熱帶的紙漿林(桉樹林),也高于位置比較接近的海南島尖峰嶺熱帶山地雨林。

    3.2 橡膠林生態(tài)系統(tǒng)NEE對環(huán)境因子的響應(yīng)

    半小時尺度上,NEE與PAR呈負(fù)相關(guān),短時間內(nèi)PAR影響生態(tài)系統(tǒng)光合作用相當(dāng)于在眾多葉片的光合作用的加和。利用Michaelis-Menten模型方程擬合,雨季橡膠林生長旺盛,從4~10月,初始光能利用率α和生態(tài)系統(tǒng)呼吸速率Rd呈遞減趨勢,橡膠林生態(tài)系統(tǒng)初始光能利用率α、最大光合作用速率Pmax、生態(tài)系統(tǒng)呼吸速率Rd最大值分別是0.005 136 mg CO2·μmol·photon-1(4月)、4.596 569 mg CO2/(m2·s)(5月)和0.190 734 mg CO2/(m2·s)(4月),大于可比較的溫帶針闊混交林和寒溫帶針葉林生態(tài)系統(tǒng),熱帶橡膠林生態(tài)系統(tǒng)光合與呼吸作用活躍。并且橡膠林生態(tài)系統(tǒng)光合作用相關(guān)參數(shù)最大值出現(xiàn)的月份要早于其它森林生態(tài)系統(tǒng)出現(xiàn)的時間,與熱帶緯度位置較低導(dǎo)致的PAR狀況及相應(yīng)的氣溫和降水的協(xié)同影響相關(guān)。

    半小時尺度上,NEE與冠層內(nèi)大氣平均溫度、地表以下5 cm土壤氣溫(T_-5 cm)也呈負(fù)相關(guān)關(guān)系,但其相關(guān)系數(shù)小于其與PAR的相關(guān)系數(shù)。熱帶橡膠林生態(tài)系統(tǒng)地處熱帶,終年大氣溫度和土壤溫度均相對較高,變化不明顯,故生態(tài)系統(tǒng)光合作用和呼吸作用受溫度因素影響均相對較弱。

    水分因素一直是影響橡膠林生態(tài)系統(tǒng)的重要因素。在短期或一定時期內(nèi)與NEE的相關(guān)性不顯著,研究表明NEE與VPD和地表以下5 cm土壤含水量(VWC_-5 cm)均呈極弱的負(fù)相關(guān);但在月尺度上,各月NEE與當(dāng)月降水量呈較強的負(fù)相關(guān)。短時間尺度上,比如橡膠樹生長季節(jié)的雨季,水分并不是NEE的限制因子;但在長時間尺度(如月尺度),雨季降水豐沛(協(xié)同PAR和氣溫等),凈生態(tài)系統(tǒng)光合能力得到加強,對NEE影響明顯。但雨季降水過多,晴天減少,也會導(dǎo)致PAR減弱、氣溫促進作用減弱。

    另外,橡膠林本身林分狀況也會影響生態(tài)系統(tǒng)NEE,本身葉面積指數(shù)與NEE也呈負(fù)相關(guān)關(guān)系,隨著橡膠林本身物候變化,LAI變化影響當(dāng)年NEE變化。

    空氣CO2濃度對橡膠林NEE影響復(fù)雜。白天CO2濃度較低,NEE為負(fù)(吸收CO2);夜間CO2濃度較高,NEE為正(釋放CO2);晝夜交替時段,NEE變化較為劇烈(此時CO2濃度相對較高)。上述白天CO2濃度低于夜間CO2濃度,與查同剛[31]研究楊樹人工林白天環(huán)境CO2濃度高于夜間CO2濃度不同(不能完全解釋NEE與CO2濃度的變化關(guān)系),橡膠林生態(tài)系統(tǒng)NEE與空氣CO2濃度的關(guān)系較好地解釋了生態(tài)系統(tǒng)本身的光合作用與呼吸作用固定與釋放CO2。因此可說橡膠林環(huán)境空氣CO2濃度與生態(tài)系統(tǒng)NEE是互為因果。

    致 謝 感謝中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所王介民研究員和北京天諾基業(yè)甄曉杰工程師在儀器維護、EdiRe數(shù)據(jù)處理程序完善等方面提供的無私幫助。

    參考文獻

    [1] 吳志祥, 謝貴水, 陶忠良. 中國天然橡膠可持續(xù)發(fā)展思路探討[J]. 世界農(nóng)業(yè), 2009(2): 47-49.

    [2] 蔣菊生, 王如松. 橡膠林固定CO2和釋放O2的服務(wù)功能及其價值估計[J]. 生態(tài)學(xué)報, 2002, 22(9): 1 545-1 551.

    [3] 張 敏, 鄒曉明. 熱帶季節(jié)雨林與人工橡膠林土壤碳氮比較[J]. 應(yīng)用生態(tài)學(xué)報, 2009, 20(5): 1 013-1 019.

    [4] Cheng C, Wang R, Jiang J. Variation of soil fertility and carbon sequestration by planting Hevea brasiliensis in Hainan[J]. Journal of Environmental Sciences, 2007, 19(3): 348-352.

    [5] Guardiola Claramonte M, Troch P A, Ziegler A D, et al. Hydrologic effects of the expansion of rubber(Hevea brasiliensis)in a tropical catchment[J]. Ecohydrology, 2010, 3(3): 306-314.

    [6] Li H, Ma Y, Liu W, et al. Soil Changes Induced by Rubber and Tea Plantation Establishment: Comparison with Tropical Rain Forest Soil in Xishuangbanna, SW China[J]. Environmental Management, 2012, 50(5): 837-848.

    [7] 劉 珉, 吳志祥, 楊 川, 等. 海南島橡膠林水分循環(huán)過程特征研究[J]. 中國農(nóng)學(xué)通報, 2012, 28(22): 27-33.

    [8] Yi Z F, Cannon C H, Chen J, et al. Developing indicators of economic value and biodiversity loss for rubber plantations in Xishuangbanna, southwest China: A case study from Menglun township[J]. Ecological Indicators, 2014, 36: 788-797.

    [9] de Blécourt M, Brumme R, Xu J, et al. Soil Carbon Stocks Decrease following Conversion of Secondary Forests to Rubber (Hevea brasiliensis)Plantations[J]. PloS one, 2013, 8(7): e69357.

    [10] 吳志祥, 謝貴水, 楊 川, 等. 橡膠林生態(tài)系統(tǒng)干季微氣候特征和通量的初步觀測[J]. 熱帶作物學(xué)報, 2010, 31(12): 2 081-2 090.

    [11] Falge E, Baldocchi D D, Olson R, et al. Gap filling strategies for long term energy flux data sets, a short communication[J]. Agricultural and Forest Meteorology, 2001a, 107(1): 71-77.

    [12] Falge E R, Baldocchi D, Olson R, et al. Gap filling strategies for defensible annual sums of net ecosystem exchange[J]. Agricultural and Forest Meteorology, 2001b, 107(1): 43-69.

    [13] Baldocchi D D, Hicks B B, Meyers T P. Measuring biosphere-atmosphere exchanges of biological related gases with micrometeorological methods[J]. Ecology, 1988, 69: 1 331-1 340.

    [14] Grelle A, Lindroth A. Eddy-correlation system for long-term monitoring of fluxes of heat, water vapour and CO2[J]. Global Change Biology, 1996, 2: 297-308.

    [15] Mcmillen R T. An eddy correlation technique with extended applicability to non-simple terrain[J]. Boundary-Layer Meteorology, 1988, 43: 231-245.

    [16] Aubinet M, Grelle A, Ibrom A, et al. Estimates of the annual net carbon and water exchange of forests: the EUROFLUX methodology[J]. Advances in Ecological Research, 2000, 30: 113-175.

    [17] Webb E K, Pearman G I, Leuning R. Correction of flux measurements for density effects due to heat and water vapour transfer[J]. Quarterly Journal of the Royal Meteorological Society, 1980, 106: 85-100.

    [18] Pilegaard K, Hummelsho/j P, Jensen N O, et al. Two years of continuous CO2 eddy-flux measurements over a Danish beech forest[J]. Agricultural and Forest Meteorology, 2001, 107(1): 29-41.

    [19] Carrara A, Kowalski A S, Neirynck J, et al. Net ecosystem CO2 exchange of mixed forest in Belgium over 5 years[J]. Agricultural and Forest Meteorology, 2003, 119(3): 209-227.

    [20] Michaelis L, Menten M L. Die kinetik der invertinwirkung[J]. Biochem. z, 1913, 49(333-369): 352.

    [21] Lasslop G, Reichstein M, Papale D, et al. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: critical issues and global evaluation[J]. Global Change Biology, 2010, 16(1): 187-208.

    [22] Hoff vant J H. Lectures on theoretical and physical chemistry[M]. Part 1. Chemical dynamics. London: Edward Arnold, 1898: 224-229.

    [23] Ojanen P, Minkkinen K, Lohila A, et al. Chamber measured soil respiration: A useful tool for estimating the carbon balance of peatland forest soils?[J]. Forest Ecology and Management, 2012, 277: 132-140.

    [24] 魏 遠(yuǎn), 張旭東, 江澤平, 等. 湖南岳陽地區(qū)楊樹人工林生態(tài)系統(tǒng)凈碳交換季節(jié)動態(tài)研究[J]. 林業(yè)科學(xué)研究, 2010, 23(5): 656-665.

    [25] Rodrigues A, Pita A, Mateus J. Turbulent fluxes of carbon dioxide and water vapour over an eucalyptus forest in Portugal[J]. Silva Lusitana, 2005, 13(2): 169-180.

    [26] Cabral O M, Gash J H, Rocha H R, et al. Fluxes of CO2 above a plantation of Eucalyptus in southeast Brazil[J]. Agricultural and Forest Meteorology, 2011, 151(1): 49-59.

    [27] 陳德祥. 尖峰嶺熱帶山地雨林碳交換的動態(tài)特征和影響因素研究[D]. 北京: 中國林業(yè)科學(xué)研究院, 2010.

    [28] 周麗艷, 賈丙瑞, 曾 偉, 等. 原始興安落葉松林生長季凈生態(tài)系統(tǒng)CO2交換及其光響應(yīng)特征[J]. 生態(tài)學(xué)報, 2010a, 30(24): 6 919-6 926.

    [29] Zhang J, Han S, Yu G. Seasonal variation in carbon dioxide exchange over a 200-year-old Chinese broad-leaved Korean pine mixed forest[J]. Agricultural and Forest Meteorology, 2006, 137(3): 150-165.

    [30] 關(guān)德新, 吳家兵, 金昌杰, 等. 長白山紅松針闊混交林CO2通量的日變化與季節(jié)變化[J]. 林業(yè)科學(xué), 2006, 42(10): 123-128.

    [31] 查同剛. 北京大興楊樹人工林生態(tài)系統(tǒng)碳平衡的研究[D]. 北京: 北京林業(yè)大學(xué), 2007.

    猜你喜歡
    影響因子海南島
    海南島上石風(fēng)流
    寶藏(2017年4期)2017-05-17 03:33:57
    基于個性化的協(xié)同過濾圖書推薦算法研究
    農(nóng)作物病蟲害氣象環(huán)境影響因子分析
    手機閱讀平臺用戶體驗影響因子分析
    目的論視角下旅游英語的語言特色對漢譯旅游文本的啟示
    科技視界(2016年23期)2016-11-04 13:29:09
    云霧物理生長過程及其影響因子的虛擬仿真實驗
    考試周刊(2016年50期)2016-07-12 23:09:38
    “影響因子”是用來賺大錢的
    讀書(2016年5期)2016-05-21 22:36:29
    是海南省還是海南島?
    關(guān)于海南島戰(zhàn)役作戰(zhàn)方針的考察
    軍事歷史(1994年6期)1994-08-15 08:56:38
    回憶海南島渡海登陸戰(zhàn)役
    軍事歷史(1987年2期)1987-08-20 06:10:48
    中文在线观看免费www的网站| 三级国产精品欧美在线观看| 久久这里只有精品中国| 日韩欧美一区二区三区在线观看| 99国产精品一区二区蜜桃av| 美女高潮的动态| 黄色丝袜av网址大全| 深夜精品福利| 好男人电影高清在线观看| 亚洲色图av天堂| 亚洲无线观看免费| 18+在线观看网站| 免费av观看视频| 91狼人影院| 国产毛片a区久久久久| 三级男女做爰猛烈吃奶摸视频| 成人特级黄色片久久久久久久| 老熟妇仑乱视频hdxx| 又粗又爽又猛毛片免费看| 美女大奶头视频| 亚洲综合色惰| 深夜精品福利| av在线天堂中文字幕| 精品久久久久久久末码| 可以在线观看的亚洲视频| 久99久视频精品免费| 国产亚洲欧美在线一区二区| 一区二区三区激情视频| 国产精品国产高清国产av| 黄色配什么色好看| 女同久久另类99精品国产91| 999久久久精品免费观看国产| 中文字幕av成人在线电影| 搡老岳熟女国产| 久久婷婷人人爽人人干人人爱| 亚洲精品色激情综合| a级毛片a级免费在线| 丝袜美腿在线中文| 日韩亚洲欧美综合| 久久99热6这里只有精品| 免费av观看视频| 97碰自拍视频| 观看美女的网站| 99热这里只有是精品在线观看 | 国产精品影院久久| 观看美女的网站| 国产成人欧美在线观看| 搞女人的毛片| 国产一区二区在线观看日韩| 免费搜索国产男女视频| 精品久久久久久久久久久久久| 免费无遮挡裸体视频| 午夜福利欧美成人| 两个人视频免费观看高清| 欧美中文日本在线观看视频| 国产白丝娇喘喷水9色精品| 少妇丰满av| 久久6这里有精品| 国产单亲对白刺激| 国产白丝娇喘喷水9色精品| 久久亚洲精品不卡| 亚洲精品久久国产高清桃花| 国产在线精品亚洲第一网站| 久久精品人妻少妇| 亚洲无线在线观看| 极品教师在线视频| 又紧又爽又黄一区二区| 天堂动漫精品| 非洲黑人性xxxx精品又粗又长| 真人做人爱边吃奶动态| 久久人人精品亚洲av| 99国产综合亚洲精品| 成人av在线播放网站| 免费在线观看日本一区| 免费看a级黄色片| 国产一区二区三区在线臀色熟女| 国产v大片淫在线免费观看| 国产探花极品一区二区| 久久精品国产自在天天线| 欧美在线一区亚洲| 亚洲欧美日韩高清在线视频| 国产精品影院久久| 免费搜索国产男女视频| 高潮久久久久久久久久久不卡| 精品国产亚洲在线| or卡值多少钱| 老司机午夜福利在线观看视频| 搞女人的毛片| 免费电影在线观看免费观看| 亚洲人与动物交配视频| 国产亚洲精品av在线| ponron亚洲| 国产真实伦视频高清在线观看 | 啪啪无遮挡十八禁网站| 成年女人毛片免费观看观看9| 国产探花极品一区二区| 在线观看一区二区三区| 免费人成视频x8x8入口观看| 国产又黄又爽又无遮挡在线| 中文字幕av成人在线电影| 舔av片在线| 免费在线观看亚洲国产| 欧美激情国产日韩精品一区| av欧美777| 18美女黄网站色大片免费观看| 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 久久人人爽人人爽人人片va | 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 亚洲欧美日韩高清专用| 亚洲一区二区三区不卡视频| 亚洲国产精品999在线| 欧美在线一区亚洲| 日本 av在线| 亚洲成人久久性| 男女做爰动态图高潮gif福利片| 啪啪无遮挡十八禁网站| 成人欧美大片| 男人舔奶头视频| 国产一区二区三区在线臀色熟女| 亚洲第一电影网av| 日韩欧美 国产精品| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| av天堂中文字幕网| 亚洲自偷自拍三级| 午夜激情福利司机影院| 亚洲国产精品合色在线| 美女高潮喷水抽搐中文字幕| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 国内揄拍国产精品人妻在线| 亚洲欧美日韩高清专用| 成人永久免费在线观看视频| 一进一出好大好爽视频| 国内精品久久久久久久电影| 内射极品少妇av片p| 在线观看一区二区三区| 狂野欧美白嫩少妇大欣赏| 在线观看舔阴道视频| 色噜噜av男人的天堂激情| 美女xxoo啪啪120秒动态图 | 99国产综合亚洲精品| 一区二区三区四区激情视频 | 久久九九热精品免费| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品久久男人天堂| 国产亚洲精品av在线| 12—13女人毛片做爰片一| 久久婷婷人人爽人人干人人爱| 日韩欧美免费精品| 三级毛片av免费| 中文在线观看免费www的网站| 日韩精品中文字幕看吧| 最好的美女福利视频网| 国产高清视频在线观看网站| 岛国在线免费视频观看| 成年人黄色毛片网站| 国产精品久久久久久久久免 | 亚洲精品一区av在线观看| 免费高清视频大片| 国产精品久久久久久久久免 | 99国产精品一区二区蜜桃av| 天堂av国产一区二区熟女人妻| 久久人人精品亚洲av| 亚洲五月天丁香| 怎么达到女性高潮| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 一本综合久久免费| 久久人人爽人人爽人人片va | 91字幕亚洲| 天美传媒精品一区二区| 免费人成在线观看视频色| 亚洲真实伦在线观看| 中亚洲国语对白在线视频| 十八禁国产超污无遮挡网站| 丁香欧美五月| 简卡轻食公司| 又紧又爽又黄一区二区| 久久久久久久午夜电影| 亚洲第一电影网av| 婷婷精品国产亚洲av在线| 欧美黑人巨大hd| 国产色婷婷99| 亚洲精品成人久久久久久| 午夜免费激情av| 亚洲国产日韩欧美精品在线观看| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色| 亚洲精品一卡2卡三卡4卡5卡| 五月玫瑰六月丁香| 91麻豆av在线| 神马国产精品三级电影在线观看| 禁无遮挡网站| 看十八女毛片水多多多| 18禁在线播放成人免费| 久久久久久久亚洲中文字幕 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美成人性av电影在线观看| 啦啦啦韩国在线观看视频| 国产午夜精品久久久久久一区二区三区 | 亚洲在线观看片| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 老司机深夜福利视频在线观看| 久久这里只有精品中国| 欧美丝袜亚洲另类 | 国产精品一及| 两人在一起打扑克的视频| 欧美区成人在线视频| 极品教师在线免费播放| 高潮久久久久久久久久久不卡| 一级a爱片免费观看的视频| 在现免费观看毛片| 两性午夜刺激爽爽歪歪视频在线观看| 黄片小视频在线播放| www.www免费av| 国产色爽女视频免费观看| 伦理电影大哥的女人| 男人舔女人下体高潮全视频| 一区二区三区免费毛片| 黄色丝袜av网址大全| 欧美激情久久久久久爽电影| 一区福利在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成人啪精品午夜网站| 看十八女毛片水多多多| 97人妻精品一区二区三区麻豆| 淫秽高清视频在线观看| 亚洲成av人片免费观看| 国产精品一区二区免费欧美| av天堂中文字幕网| 精品一区二区免费观看| 男插女下体视频免费在线播放| av黄色大香蕉| 变态另类成人亚洲欧美熟女| 色综合亚洲欧美另类图片| 嫩草影院入口| 色尼玛亚洲综合影院| 午夜两性在线视频| 97超视频在线观看视频| 人妻制服诱惑在线中文字幕| 成人高潮视频无遮挡免费网站| 婷婷色综合大香蕉| 老司机午夜十八禁免费视频| 免费在线观看影片大全网站| 嫁个100分男人电影在线观看| 老司机深夜福利视频在线观看| 免费高清视频大片| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 99riav亚洲国产免费| 欧美日韩国产亚洲二区| 女生性感内裤真人,穿戴方法视频| 校园春色视频在线观看| 嫩草影视91久久| 天天躁日日操中文字幕| 日韩国内少妇激情av| 国产一区二区在线av高清观看| 高清在线国产一区| av国产免费在线观看| 可以在线观看的亚洲视频| 在线观看免费视频日本深夜| 丝袜美腿在线中文| 天天一区二区日本电影三级| 最近视频中文字幕2019在线8| 中文字幕人妻熟人妻熟丝袜美| 久久人人精品亚洲av| 别揉我奶头~嗯~啊~动态视频| 白带黄色成豆腐渣| 69人妻影院| 久久久久久久久中文| 欧美黄色片欧美黄色片| 亚洲不卡免费看| 村上凉子中文字幕在线| 在线看三级毛片| 1000部很黄的大片| 桃红色精品国产亚洲av| 国产精品1区2区在线观看.| 国产精品女同一区二区软件 | 亚洲成人免费电影在线观看| 亚洲精华国产精华精| 欧美乱色亚洲激情| 琪琪午夜伦伦电影理论片6080| 老女人水多毛片| av国产免费在线观看| 老司机午夜十八禁免费视频| 亚洲精品亚洲一区二区| 精品熟女少妇八av免费久了| 1000部很黄的大片| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 久久亚洲真实| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区久久| 亚洲av成人精品一区久久| 国内揄拍国产精品人妻在线| 校园春色视频在线观看| 少妇丰满av| 首页视频小说图片口味搜索| 午夜影院日韩av| 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 嫩草影院入口| .国产精品久久| 性色av乱码一区二区三区2| 在线国产一区二区在线| 国产精品爽爽va在线观看网站| 99热精品在线国产| 99在线人妻在线中文字幕| 久久人妻av系列| 国产伦精品一区二区三区视频9| 两人在一起打扑克的视频| 91午夜精品亚洲一区二区三区 | 国产av麻豆久久久久久久| 99热只有精品国产| 日本一本二区三区精品| 色吧在线观看| 真实男女啪啪啪动态图| 色吧在线观看| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 久久人妻av系列| 亚洲综合色惰| 国产又黄又爽又无遮挡在线| 99热只有精品国产| 淫妇啪啪啪对白视频| 国产精品,欧美在线| 亚洲精品乱码久久久v下载方式| 久久国产乱子免费精品| 欧美3d第一页| 2021天堂中文幕一二区在线观| 午夜激情欧美在线| 婷婷精品国产亚洲av在线| 欧美xxxx性猛交bbbb| 日本黄大片高清| 午夜福利视频1000在线观看| 亚洲第一区二区三区不卡| 久久久久九九精品影院| 中文亚洲av片在线观看爽| 欧美黑人巨大hd| 国产精品一区二区三区四区免费观看 | 日本五十路高清| 精品熟女少妇八av免费久了| 成人国产一区最新在线观看| 无人区码免费观看不卡| 日韩精品青青久久久久久| 久99久视频精品免费| 每晚都被弄得嗷嗷叫到高潮| 一级作爱视频免费观看| 国产乱人视频| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 午夜精品久久久久久毛片777| av在线蜜桃| 国产精品嫩草影院av在线观看 | 免费看a级黄色片| 久久久久性生活片| 自拍偷自拍亚洲精品老妇| 久久久久性生活片| 亚洲国产精品999在线| 欧美3d第一页| 国产不卡一卡二| 精品99又大又爽又粗少妇毛片 | 男女视频在线观看网站免费| 国产精品98久久久久久宅男小说| 十八禁人妻一区二区| 全区人妻精品视频| 亚洲最大成人手机在线| 91久久精品电影网| 欧美xxxx性猛交bbbb| 成人国产一区最新在线观看| 亚洲成a人片在线一区二区| 国产精品影院久久| 免费看日本二区| 精华霜和精华液先用哪个| 久久久国产成人精品二区| 久久国产乱子免费精品| 国产成人啪精品午夜网站| 亚洲美女视频黄频| 国产在视频线在精品| 色吧在线观看| 啦啦啦韩国在线观看视频| 床上黄色一级片| 亚洲av中文字字幕乱码综合| 俄罗斯特黄特色一大片| 黄色配什么色好看| 久久久久国内视频| 国内精品美女久久久久久| 国产淫片久久久久久久久 | 美女cb高潮喷水在线观看| 国产极品精品免费视频能看的| 亚洲av免费高清在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲av嫩草精品影院| 久久精品91蜜桃| 午夜精品在线福利| 午夜福利高清视频| 免费在线观看成人毛片| 亚洲成人中文字幕在线播放| 日韩有码中文字幕| 国产午夜福利久久久久久| 天美传媒精品一区二区| 国产高清三级在线| 波多野结衣高清无吗| 日日干狠狠操夜夜爽| 韩国av一区二区三区四区| 国产三级黄色录像| 欧美最黄视频在线播放免费| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 国产极品精品免费视频能看的| 亚洲七黄色美女视频| 国产欧美日韩精品一区二区| av欧美777| 很黄的视频免费| 日韩欧美免费精品| 看免费av毛片| www.色视频.com| 成人鲁丝片一二三区免费| 两个人的视频大全免费| 国产精品久久久久久人妻精品电影| 淫妇啪啪啪对白视频| 男女做爰动态图高潮gif福利片| 99久久精品国产亚洲精品| 我的老师免费观看完整版| 国产高清有码在线观看视频| 国产在线男女| 亚州av有码| 国产色婷婷99| 久久久久久久午夜电影| 日韩国内少妇激情av| 性插视频无遮挡在线免费观看| 亚洲中文字幕日韩| 狂野欧美白嫩少妇大欣赏| 亚洲成人精品中文字幕电影| 欧美国产日韩亚洲一区| АⅤ资源中文在线天堂| 在线播放无遮挡| 两个人视频免费观看高清| 三级男女做爰猛烈吃奶摸视频| 亚洲成av人片免费观看| av天堂在线播放| 国产成人欧美在线观看| 欧美午夜高清在线| 美女被艹到高潮喷水动态| 欧美老熟妇乱子伦牲交| 久久综合国产亚洲精品| 老司机影院毛片| 好男人在线观看高清免费视频| 中文字幕av成人在线电影| 又爽又黄无遮挡网站| 少妇人妻 视频| 国产 一区 欧美 日韩| 国产精品人妻久久久影院| 亚洲丝袜综合中文字幕| 精品少妇久久久久久888优播| 国产精品.久久久| 国产亚洲av嫩草精品影院| 亚洲激情五月婷婷啪啪| 亚洲av不卡在线观看| 最近最新中文字幕免费大全7| 一级毛片aaaaaa免费看小| 一级毛片久久久久久久久女| 在线a可以看的网站| 亚洲国产精品成人久久小说| 精品人妻视频免费看| 高清在线视频一区二区三区| 97在线视频观看| 日韩国内少妇激情av| 性色avwww在线观看| 天天躁夜夜躁狠狠久久av| 美女视频免费永久观看网站| 99热网站在线观看| 赤兔流量卡办理| 午夜免费鲁丝| 99视频精品全部免费 在线| 国产 精品1| 国产 一区精品| 乱码一卡2卡4卡精品| 久久久国产一区二区| 亚洲性久久影院| 一区二区三区免费毛片| 久久精品人妻少妇| 免费黄网站久久成人精品| 在线观看三级黄色| 特级一级黄色大片| 美女cb高潮喷水在线观看| 国产黄色免费在线视频| 成人鲁丝片一二三区免费| 国产精品久久久久久精品电影| 可以在线观看毛片的网站| 免费大片18禁| 人体艺术视频欧美日本| 久久精品国产亚洲网站| www.色视频.com| 波多野结衣巨乳人妻| 国产黄片视频在线免费观看| 欧美三级亚洲精品| 成年免费大片在线观看| 最近最新中文字幕免费大全7| 免费看不卡的av| 欧美成人a在线观看| 国产爱豆传媒在线观看| 国产v大片淫在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 国产精品人妻久久久久久| 亚洲,一卡二卡三卡| av.在线天堂| 色视频www国产| 亚洲av男天堂| 中国美白少妇内射xxxbb| 国模一区二区三区四区视频| freevideosex欧美| 久久精品国产a三级三级三级| 夫妻午夜视频| 国产白丝娇喘喷水9色精品| 亚洲国产色片| 最新中文字幕久久久久| 亚洲三级黄色毛片| 91在线精品国自产拍蜜月| 大话2 男鬼变身卡| 欧美成人精品欧美一级黄| 精品午夜福利在线看| 亚洲性久久影院| 春色校园在线视频观看| 精品人妻视频免费看| 成人美女网站在线观看视频| 超碰av人人做人人爽久久| 免费观看av网站的网址| 一级av片app| 自拍欧美九色日韩亚洲蝌蚪91 | 男女无遮挡免费网站观看| av在线播放精品| 91久久精品电影网| 内地一区二区视频在线| 精品一区在线观看国产| 亚洲经典国产精华液单| 最近2019中文字幕mv第一页| 菩萨蛮人人尽说江南好唐韦庄| 人妻 亚洲 视频| 天天躁日日操中文字幕| 2022亚洲国产成人精品| 老女人水多毛片| 欧美日韩视频高清一区二区三区二| av.在线天堂| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 搞女人的毛片| 亚洲国产色片| 亚洲怡红院男人天堂| 欧美高清成人免费视频www| 一级毛片我不卡| 亚洲精品,欧美精品| 国产高清三级在线| 免费观看性生交大片5| 日韩欧美 国产精品| 亚洲美女搞黄在线观看| 97在线视频观看| 久久久久久久久久人人人人人人| 久久99热6这里只有精品| 日韩三级伦理在线观看| 久久鲁丝午夜福利片| 亚洲国产色片| 九九久久精品国产亚洲av麻豆| 小蜜桃在线观看免费完整版高清| 最后的刺客免费高清国语| 亚洲成人精品中文字幕电影| 男人添女人高潮全过程视频| 一个人观看的视频www高清免费观看| 男人舔奶头视频| 最近最新中文字幕大全电影3| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| av线在线观看网站| 亚洲av电影在线观看一区二区三区 | 在线观看美女被高潮喷水网站| 国产精品一区二区性色av| 免费黄网站久久成人精品| 搡老乐熟女国产| 久久女婷五月综合色啪小说 | 国内少妇人妻偷人精品xxx网站| 久久精品国产a三级三级三级| 欧美性感艳星| 亚洲欧美精品自产自拍| 亚洲人成网站在线播| av一本久久久久| 大陆偷拍与自拍| 天堂俺去俺来也www色官网| 国产亚洲91精品色在线| 亚洲成色77777| 麻豆乱淫一区二区| 亚洲一级一片aⅴ在线观看| 久久久久国产网址| 亚洲欧美日韩另类电影网站 | 青春草亚洲视频在线观看| 寂寞人妻少妇视频99o| 内地一区二区视频在线| 深爱激情五月婷婷| 大话2 男鬼变身卡| 丝瓜视频免费看黄片| 超碰av人人做人人爽久久| 91久久精品电影网| 中文在线观看免费www的网站| av国产精品久久久久影院| 亚洲一区二区三区欧美精品 | 免费大片18禁| 久久久久久九九精品二区国产| 亚洲国产精品国产精品| 久久99热这里只频精品6学生| 亚洲精品自拍成人|