• <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)凈碳交換特征及其驅(qū)動因子

    2015-06-26 11:00:43劉廷璽雷慧閩張圣微段利民
    草業(yè)學(xué)報(bào) 2015年11期
    關(guān)鍵詞:草甸葉面積通量

    王 婧,劉廷璽*,雷慧閩,張圣微,段利民

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古呼和浩特010018;2.清華大學(xué)水利水電工程系,北京100084)

    科爾沁草甸生態(tài)系統(tǒng)凈碳交換特征及其驅(qū)動因子

    王 婧1,劉廷璽1*,雷慧閩2*,張圣微1,段利民1

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古呼和浩特010018;2.清華大學(xué)水利水電工程系,北京100084)

    以2013年10月-2014年9月連續(xù)觀測的CO2通量數(shù)據(jù)為基礎(chǔ),分析了科爾沁草甸生態(tài)系統(tǒng)凈碳交換量(NEE)的時間變化特征及其驅(qū)動因素。結(jié)果表明,NEE日變化季節(jié)差異明顯,生長季變化幅度大,凈CO2日吸收速率7月>8月>9月>6月>5月;生長季內(nèi),NEE主要受控于葉面積指數(shù)和光合有效輻射。NEE與光合有效輻射(PAR)之間的關(guān)系可用直角雙曲線方程來描述,擬合得到的表觀初始光能利用率α為0.0015μmol CO2/μmol PAR,最大光合速率Pmax為0.65μmol CO2/(m2·s)。葉面積指數(shù)(LAI)對NEE的影響可由分段函數(shù)表示,當(dāng)LAI>3.08時,表現(xiàn)為漸進(jìn)飽和型,且LAI越大NEE對PAR的響應(yīng)越明顯;當(dāng)高飽和水汽壓差(VPD)在1.5~2.0 kPa時,光合作用開始降低,NEE明顯受到VPD值的抑制;短暫強(qiáng)降雨(累計(jì)降雨量>40 mm/d)對晝間NEE有一定的抑制,而持續(xù)低強(qiáng)度降雨(降雨時長>15 h)對夜間NEE存在激發(fā)作用;夜間NEE隨土壤溫度呈指數(shù)增長,溫度敏感系數(shù)(Q10)為2.63。

    凈生態(tài)系統(tǒng)交換量;驅(qū)動因子;渦度相關(guān)系統(tǒng);科爾沁沙地

    陸地生態(tài)系統(tǒng)的碳循環(huán)及其對環(huán)境因子的響應(yīng)研究是理解生態(tài)系統(tǒng)碳平衡的關(guān)鍵。草地生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的重要組成部分,幾乎覆蓋了陸地面積的20%,其凈初級生產(chǎn)力約占全球陸地生物區(qū)凈初級生產(chǎn)力的1/3。美歐科學(xué)家于1996年初步建立起全球CO2通量網(wǎng),利用渦度技術(shù)在全球范圍內(nèi)開展了不同緯度、不同類型生態(tài)系統(tǒng)CO2通量的全球性研究,但其中70%是針對不同森林植被類型的,對草原生態(tài)系統(tǒng)的研究主要集中在低海拔的平原地區(qū)[1-2],且研究較少。Dugas等[2]對德高草原碳收支狀況研究發(fā)現(xiàn)葉面積指數(shù)(LAI)、光合有效輻射(PAR)及土壤含水量影響碳通量大小,Shaver等[3]應(yīng)用Pan-Arctic模型對生態(tài)系統(tǒng)凈生態(tài)系統(tǒng)碳交換(以下簡稱“NEE”)研究表明,NEE的變異約75%可由LAI、氣溫(Ta)、PAR函數(shù)的單一回歸模型解釋。我國自20世紀(jì)90年代開始應(yīng)用渦度相關(guān)法開展了陸面過程中碳循環(huán)通量的研究,盡管起步稍晚,但也取得了部分有價值的研究成果。伏玉玲等[4]研究了中國草地生態(tài)系統(tǒng)碳收支及其對環(huán)境的響應(yīng),郝彥賓等[5]研究了內(nèi)蒙古羊草(Leymus chinensis)草原碳交換的季節(jié)變異規(guī)律,發(fā)現(xiàn)水分和PAR是影響碳交換的兩個重要生態(tài)因子。Zhao等[6]研究了青藏高寒草甸凈碳交換特征,表明葉面積指數(shù)(GLAI)、土壤含水量(θ)、Ta是3個對NEE季節(jié)變異影響最大的指標(biāo)。

    科爾沁沙地草甸區(qū)作為我國北方典型的半干旱草原區(qū),具有極其脆弱的草原生態(tài)資源。與其他下墊面相比,對氣候變化的響應(yīng)更為敏感。近年來,有關(guān)該地區(qū)的碳通量研究主要局限于短期、不連續(xù)的觀測數(shù)據(jù)[7]及碳排放的觀測[8],針對NEE的研究比較缺乏。而研究半干旱區(qū)草甸生態(tài)系統(tǒng)NEE的日、季變化及其與氣象因子的關(guān)系,對于了解碳的源匯過程、開展草原碳循環(huán)的模擬和氣候變化情景下的預(yù)測研究具有十分重要的意義。本文依據(jù)渦度相關(guān)系統(tǒng)2013年10月-2014年9月的連續(xù)觀測數(shù)據(jù),對科爾沁草甸生態(tài)系統(tǒng)的碳收支及其過程機(jī)制進(jìn)行定量評價,分析各氣象要素對碳交換的影響,揭示草甸草原NEE的日、季變化規(guī)律,解析主要環(huán)境因子對碳源匯的影響及其差異,為評估該地區(qū)碳收支和明晰全球碳循環(huán)提供更為完善的試驗(yàn)數(shù)據(jù)和理論基礎(chǔ)。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗(yàn)研究區(qū)隸屬內(nèi)蒙古通遼市科爾沁左翼后旗,面積約55 km2。該區(qū)多年平均氣溫約6.6℃,1月平均氣溫最低,為-13.3℃,7月平均氣溫最高,為23.8℃;多年平均降水量389 mm,降水主要集中在6-8月,占年總降水量的69.3%;多年平均蒸發(fā)量1412 mm,蒸發(fā)主要集中在4-9月。此間蒸發(fā)量占全年的17%。試驗(yàn)區(qū)為沙丘-草甸相間分布地區(qū),草甸面積約占整個研究區(qū)面積的26.6%,本文選取草甸試驗(yàn)站(122°39′15″E,43°20′14″N)作為研究對象,該試驗(yàn)站地面相對高程約188.5 m,南臨小型湖泊,地勢平坦。天然植被主要有羊草、蘆葦(Phragmites australis)等。6-8月份植被生長旺季,植被平均高度為145 cm,平均葉面積指數(shù)為2.47。

    1.2 觀測儀器與方法

    草甸試驗(yàn)點(diǎn)主要安裝有渦度相關(guān)系統(tǒng)傳感器、微氣象觀測系統(tǒng)。分別在5,10,20,40,80,160 cm處埋有6層土壤溫度、土壤含水量觀測傳感器。每30 d進(jìn)行一次數(shù)據(jù)采集與處理,渦度原始數(shù)據(jù)與氣象數(shù)據(jù)同步下載,本研究主要以2013年10月-2014年9月連續(xù)觀測的CO2通量數(shù)據(jù)及微氣象數(shù)據(jù)為基礎(chǔ)。

    1.2.1 渦度相關(guān)系統(tǒng)和微氣象系統(tǒng)觀測 渦度相關(guān)系統(tǒng)傳感器安裝在距地面2.6 m處,主要用于觀測湍流通量。該系統(tǒng)主要包括三維超聲風(fēng)速儀(CSAT-3,Campbell Scientific,USA)及紅外氣體分析儀(LI-7500,Li-COR Inc,Nebraska,USA),分別用于測定三維風(fēng)速與大氣CO2和H2O密度,由數(shù)據(jù)采集器(CR3000,Campbell Scientific,USA)記錄并存儲系統(tǒng)同步計(jì)算的二氧化碳和水熱通量,采樣頻率為10 Hz,數(shù)據(jù)記錄時間間隔為30 min。

    微氣象觀測系統(tǒng)主要包括安裝在渦度觀測塔上的四分量凈輻射儀(NR-LITE,Campbell Scientific,USA)、光合有效輻射儀(LI190SB,Campbell Scientific,USA),以及3層(1.5,2.4,10.0 m)溫濕度測量儀器(HMP45C,Campbell Scientific,US)、土壤熱通量板(HFP01,Campbell Scientific,USA)、120 cm口徑水面蒸發(fā)量(255,Nova Lynx,US)、降水量(TE525MM,Campbell Scientific,US)等。這些監(jiān)測要素由數(shù)據(jù)采集器(CR1000,Campbell Scientific,USA)按10 min時間間隔計(jì)算其平均值并存儲。

    1.2.2 葉面積指數(shù)觀測 對草甸群落葉面積指數(shù)LAI的動態(tài)測定也與碳通量觀測同步進(jìn)行,利用植物冠層分析儀(LAI-2200,Li-COR Nebraska,USA)通過測定冠層上方(簡稱“A”值)、冠層下方(簡稱“B”值)光照輻射的衰減來計(jì)算葉面積指數(shù)(葉面積指數(shù)=葉面積/地面積)。草甸植被生長初期的4-5月份,每隔15~20 d測量一次,6-9月份每20~30 d測量一次。選取3個1 m×1 m的樣方,每個樣方測定1個A值,4個B值,采用90度遮蓋帽,取3組葉面積指數(shù)的平均值作為代表值。

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

    1.3.1 CO2通量的計(jì)算公式 CO2通量即凈生態(tài)系統(tǒng)碳交換量,由CO2(μmol/m2·s)/H2O(mmol/m2·s)的濃度與垂直風(fēng)速的原始數(shù)據(jù)經(jīng)過協(xié)方差計(jì)算來確定。式中,w'為垂直湍流速度脈動,c'為二氧化碳濃度。NEE由渦度相關(guān)系統(tǒng)測定的CO2通量直接獲得,NEE為正表示生態(tài)系統(tǒng)向大氣釋放CO2,反之則表示從大氣吸收CO2。

    1.3.2 通量數(shù)據(jù)的質(zhì)量控制與插補(bǔ) 對渦動相關(guān)法觀測到的10 Hz原始湍流數(shù)據(jù)運(yùn)用Logger Net軟件進(jìn)行格式轉(zhuǎn)換(TOB3-TOB1),然后采用美國LI-COR公司研發(fā)的Eddy Pro軟件[9]進(jìn)行系列校正計(jì)算,包括異常值及野值點(diǎn)剔除、坐標(biāo)旋轉(zhuǎn)修正、頻率損失修正、超聲虛溫修正、空氣密度效應(yīng)修正(WPL修正)等基本處理,得到采樣間隔為30 min的通量數(shù)據(jù)。

    對30 min數(shù)據(jù),若野點(diǎn)數(shù)大于100,則剔除該時次數(shù)據(jù);或根據(jù)LI-7500紅外氣體分析儀診斷異常標(biāo)志值(diag≠0)[10],以及相關(guān)監(jiān)測數(shù)據(jù)超出CO2濃度閾值(100~1000 mg/m3),對夜間觀測,用給定的一個臨界閾值uc*加以檢驗(yàn)。當(dāng)u*<uc*時,則舍棄該通量數(shù)據(jù)不進(jìn)行分析[11]。uc*值一般取0.1~0.2 m/s,本文取0.1 m/s。

    對原始通量數(shù)據(jù)進(jìn)行質(zhì)量控制后,剔除的低質(zhì)量數(shù)據(jù)會影響數(shù)據(jù)的連續(xù)性。因此,本研究對夜間和晝間缺失的NEE數(shù)據(jù),采用分別的方法進(jìn)行插補(bǔ)。針對夜間缺失的NEE數(shù)據(jù),利用5 cm深度處的土壤溫度(Ts)采用Vants Hoff擬合方程加以插補(bǔ)[12]。

    式中,Reco為夜間生態(tài)系統(tǒng)呼吸速率(由于夜間沒有光合作用,Reco即為NEE;Reco,ref為生態(tài)系統(tǒng)10℃時的參考呼吸速率(本文取0.05 mg/m2·s);Q10為呼吸隨溫度的變化速率(本文取2.63),即5至40℃之間,溫度每增加10℃單位時間土壤呼吸增加的倍數(shù),Raich和Schlesinger[13]通過收集和分析全球各種生態(tài)系統(tǒng)的土壤呼吸數(shù)據(jù),發(fā)現(xiàn)Q10在1.3~3.3之間波動,Ts為5 cm深度處土壤實(shí)測溫度。

    針對晝間缺失的NEE數(shù)據(jù),則利用Michaelis-Menten函數(shù)進(jìn)行插補(bǔ)[13]。

    式中,α為生態(tài)系統(tǒng)表觀光量子效率;pmax為生態(tài)系統(tǒng)最大光合速率,即PAR→∞時NEE的漸進(jìn)值;Reco,d為暗呼吸速率,即PAR→0時的NEE值。

    2 結(jié)果與分析

    2.1 凈二氧化碳通量日內(nèi)變化特征

    草甸NEE的日變化基本呈現(xiàn)單峰型(圖1),且呈現(xiàn)出上午碳吸收較強(qiáng),下午碳吸收較弱的不對稱現(xiàn)象,在日出(5:00-6:00時)、日落(18:00-19:00時)前后碳通量波動較大,此時的呼吸作用與光合作用達(dá)到一個平衡狀態(tài)。夜間NEE為正,整個生態(tài)系統(tǒng)表現(xiàn)為碳排放,白天NEE均為負(fù)值,即隨著光合作用的加強(qiáng),整個生態(tài)系統(tǒng)表現(xiàn)為凈碳吸收,5-9月日均碳吸收速率分別為0.02,0.12,0.33,0.29,0.24 mg/(m2·s),呈現(xiàn)為7月>8月>9月>6月>5月,且碳吸收出現(xiàn)在凌晨的時間依次滯后0.25,0.25,0.50,0.50 h,碳排放出現(xiàn)在傍晚的時間依次提前0.5,0.5,1.0,1.5 h。同時,日出、日落前后碳吸收量與排放量達(dá)到相對平衡的狀態(tài)。

    圖1 生長季草甸生態(tài)系統(tǒng)碳通量月平均日變化Fig.1 Diurnal changes of CO2fluxes during growing season over the meadow grassland ecosystem

    2.2 凈二氧化碳通量及其影響因子的季節(jié)變化特征

    為了解NEE與環(huán)境因子的季節(jié)變化特征,本文選取2013年10月-2014年9月生態(tài)系統(tǒng)NEE及其主要環(huán)境因子(飽和水汽壓差VPD、氣溫Ta、土壤溫度Ts、光合有效輻射PAR)的月變化加以分析(圖2)。氣溫與5 cm深度處的土壤溫度最大值均出現(xiàn)在7月中旬,分別為29.3和26.6℃,草甸生長季節(jié)的氣溫、晝夜溫差、土壤溫度的日平均值分別為19.3,10.0,16.0℃;VPD的日平均值為1.08 kPa,最高、最低值分別出現(xiàn)在7月9日和10月18日;PAR的月平均日變化呈現(xiàn)“單峰型”,最大值出現(xiàn)在7月5日,為1874μmol/(m2·s);8月降水量最高,為245.9 mm,7月次之,為104.6 mm。生長季平均葉面積指數(shù)LAI為2.57。4-6月LAI迅速增加,6-8月增速放緩,8月中下旬達(dá)到最大值3.34,9月微弱下降,但降幅較小。

    NEE在全年表現(xiàn)出“兩峰一谷”的特征(圖2d),其中“峰值”出現(xiàn)在1,10月,均值分別為0.13和0.08 mg/(m2·s),“谷值”出現(xiàn)在7月,均值-0.33 mg/(m2·s),較徐麗君等[14]研究的貝加爾針茅草(Stipa baicalensis)原碳吸收谷值偏小。4月草甸NEE均值開始由正轉(zhuǎn)負(fù),生長季的5-9月,生態(tài)系統(tǒng)主要表現(xiàn)為碳匯,凈CO2吸收總量1.02 mg/(m2·s),5,6,7,8,9月分別貢獻(xiàn)了11.2%,27.0%,30.5%,22.5%,8.4%。但不同緯度的草甸生態(tài)系統(tǒng)NEE正負(fù)轉(zhuǎn)換的時間有較大差異,如Shi等[15]研究嵩草草甸生態(tài)系統(tǒng)7-9月為碳吸收,其余為碳排放。而本研究區(qū)從10月出現(xiàn)碳排放,10月至翌年3月,凈CO2日排放均值0.11 mg/(m2·s)。同時生態(tài)系統(tǒng)碳交換因受到環(huán)境因素影響而產(chǎn)生差異(表1),如生長季NEE與PAR、LAI相關(guān)性顯著(P<0.01),非生長季NEE主要受控于溫度因素(R2=0.568),生態(tài)系統(tǒng)主要表現(xiàn)為碳源。研究時段,草甸作為碳匯從大氣中吸收218.52 g CO2/m2。

    2.3 葉面積指數(shù)對NEE的影響

    LAI是反映植物長勢的重要特征參數(shù),它對植物光能利用、干物質(zhì)積累等都有明顯影響。生長初期(LAI<0.97),由于植被稀疏LAI較小,對散射光的吸收較少,草甸CO2吸收量隨LAI的增加而線性增長;生長旺期(0.97<LAI<3.08)植被生長稠密,LAI增加冠層有利于植物吸收來自各個方向的散射光,CO2吸收量隨LAI的增加而緩慢增加,之后漸進(jìn)飽和狀態(tài);生長末期(LAI>3.08)部分葉片開始枯黃、凋落,LAI有所下降,碳吸收減少且隨LAI的增加而趨于平穩(wěn)(圖3),這與Suyker等[16]對玉米(Zea may)生殖生長研究得出的結(jié)論相似。同時由于生長末期群體底層存在大量枯葉,利用冠層分析儀進(jìn)行測量時,儀器無法自動消除這部分枯葉的影響,導(dǎo)致LAI測定結(jié)果偏大[17]。

    圖2 凈碳交換量NEE與環(huán)境因子的季節(jié)變化Fig.2 Seasonal variations of environmental factors with net carbon exchange(NEE)

    表1 凈碳交換NEE與氣象因子標(biāo)準(zhǔn)化回歸分析Table 1 Normalized regression analysis between net ecosystem CO2exchange(NEE)and meteorology factors

    生長季LAI與PAR的偏相關(guān)分析(表2)表明,當(dāng)控制PAR時,LAI與NEE相關(guān)性不顯著(R2=0.38,P>0.05),而考慮PAR影響時,LAI與NEE相關(guān)系數(shù)提高了40.6%,二者顯著相關(guān)(P<0.01),LAI與PAR存在明顯的交互作用,共同影響著生態(tài)系統(tǒng)凈碳交換過程。

    2.4 光合有效輻射對NEE的影響

    生長季的NEE與PAR變化關(guān)系較為分散,通常用一個直角雙曲線模型表示(公式3)。其中表觀初始光能效率α,利用低光強(qiáng)下NEE-PAR的線性擬合關(guān)系公式反推得到[18](NEE=0.0015PAR+0.264,R2=0.924),即α值為0.0015μmol CO2/μmol PAR。由NEE-PAR的擬合曲線(圖4a)可以估計(jì)光補(bǔ)償點(diǎn)的范圍變化在80~150μmol/(m2·s)之間。

    表2 LAI與PAR偏相關(guān)分析Table 2 Partial correlation analysis between LAI and PAR

    圖4 凈碳交換NEE與光合有效輻射PAR的關(guān)系Fig.4 Relationships between net ecosystem CO2exchange(NEE)and photosynthetically active radiation(PAR)in growing season

    草甸生態(tài)系統(tǒng)NEE對光強(qiáng)PAR的響應(yīng)分為3個不同階段(圖4b):PAR<300μmol/(m2·s)時,NEE吸收值隨PAR的升高呈線性較快速增加,這時光是唯一的環(huán)境限制因素;300<PAR<1200μmol/(m2·s)時,NEE吸收值隨PAR的升高呈曲線式較緩慢增長,此時溫度、外界CO2濃度等環(huán)境因素成為共同的限制因素;PAR>1200μmol/(m2·s)時,NEE吸收量不再隨著PAR的升高而增長,即達(dá)到光合作用的光飽和階段。這與許大全等[19]的研究結(jié)果相一致。

    2.5 水分因素對NEE的影響

    水汽壓差VPD間接地代表空氣的水分狀況,有研究指出高光強(qiáng)和高水汽飽和差是植物光合作用降低的兩個條件[20]。將生長季VPD觀測數(shù)據(jù)劃分為3個區(qū)間:低區(qū)(0.5~1.0 k Pa)、中區(qū)(1.0~1.5 k Pa)、高區(qū)(1.5~2.0 kPa),VPD由低區(qū)到高區(qū)曲線斜率依次增大-減小-增大(圖5)。VPD在低區(qū)時,隨著VPD值的升高,碳吸收速率增加,即更有利于生態(tài)系統(tǒng)的碳蓄積;VPD在中區(qū)時,NEE吸收速率隨VPD的增加繼續(xù)升高,但增加幅度越來越小,二者變化趨于飽和;VPD處于高區(qū)時,NEE吸收速率隨VPD增加開始降低,不利于生態(tài)系統(tǒng)的碳蓄積,原因可能是高水汽壓導(dǎo)致氣孔部分關(guān)閉,氣孔阻力增大,光合作用降低。進(jìn)而影響草甸生態(tài)系統(tǒng)的生產(chǎn)力[21]。

    降水是影響生態(tài)系統(tǒng)碳交換的重要水分因子,Lawerence等[22]在對加拿大北部溫帶C3草地生態(tài)系統(tǒng)連續(xù)觀測結(jié)果表明,降水充足的年份CO2日吸收峰值約為干旱條件下的3倍。同時降雨格局的差異也會影響生態(tài)系統(tǒng)的碳交換能力,以2014年8月16日(短時強(qiáng)降雨)、22日(短時弱降雨)、28日(持續(xù)降雨)為例,探討光照強(qiáng)度、降雨對NEE的共同影響(圖6)。

    圖5 凈碳交換NEE與飽和水汽壓差VPD的關(guān)系Fig.5 Relationships between net ecosystem CO2exchange(NEE)and vapor pressure deficit(VPD)

    圖6 降雨對NEE日變化的影響Fig.6 Diurnal variation of CO2exchange(NEE)under different rainfall conditions

    8月16日,降雨開始前3 h內(nèi),草甸呼吸速率由0.35 mg/(m2·s)迅速降到0.02 mg/(m2·s),降幅達(dá)94%。受凌晨降雨影響,白天PAR僅在0~600μmol/(m2·s)范圍波動,而表層土壤含水量出現(xiàn)15.8%的小幅遞增,CO2日吸收峰值為0.79 mg/(m2·s),出現(xiàn)在雨后14:00,夜間NEE與雨前碳排放速率大體持平。8月22日,白天NEE與PAR呈明顯“單峰”變化,PAR極值達(dá)1347μmol/(m2·s),正午12:30左右,到達(dá)碳吸收峰值0.93 mg/(m2·s),而夜間21:00-22:00發(fā)生了5.8 mm的微量降雨,但未對碳排放造成明顯影響。8月28日,降雨歷時15.5 h,凈輻射(Rn)與PAR均受其影響,數(shù)值僅在0~50 W/m2和0~300μmol/(m2·s)范圍內(nèi)波動,表層土壤含水量在降雨開始3 h內(nèi)緩慢增加,而后穩(wěn)定。午間13:00左右隨雨量增加,生態(tài)系統(tǒng)碳匯功能降低(-0.48 mg/m2·s)。夜間受降雨影響,生態(tài)系統(tǒng)呼吸速率明顯上升。可以看出,降雨削弱了到達(dá)地面的輻射強(qiáng)度,卻增加了表層土壤水含量,進(jìn)而影響生態(tài)系統(tǒng)碳交換。短暫強(qiáng)降雨過程使白天生態(tài)系統(tǒng)碳匯功能降低,而連續(xù)低強(qiáng)度降雨雖然使生態(tài)系統(tǒng)處于低光強(qiáng)(0<PAR<300μmol/m2·s)的環(huán)境,但提高了草甸生態(tài)系統(tǒng)夜間的碳交換能力,增加了地溫對夜間呼吸的控制[23]。

    3 討論

    生長季NEE日變化呈明顯“倒鐘型”,全天出現(xiàn)一個吸收高峰。日出后的6:00-6:30,草甸植被便開始吸收CO2,8:00-10:30之間達(dá)到碳吸收高峰,午間的11:00-13:00由于光飽和作用,碳吸收量減小,出現(xiàn)群落“午休”現(xiàn)象,這與張法偉等[24]研究的青海湖高寒草甸光合午休時間較為一致。

    12月至翌年3月草甸植被尚未復(fù)蘇,4月中下旬植被才開始返青,草甸系統(tǒng)碳吸收微弱,土壤呼吸作用顯著。曾有研究指出,非生長季土壤呼吸是生長季生態(tài)系統(tǒng)釋放CO2的一個重要組成部分,顯著影響著系統(tǒng)的碳收支[25],陳驥等[26]在對青海湖北岸高寒草甸碳收支的研究中發(fā)現(xiàn),非生長季土壤呼吸速率是生長季的土壤呼吸速率的0.30~0.33倍。而生長季5-7月生態(tài)系統(tǒng)白天的碳吸收強(qiáng)度主要受控于光合有效輻射,二者表現(xiàn)出良好的直角雙曲線關(guān)系(NEE=+0.304,R2=0.73),而夜間則通過擬合出Vants Hoff方程(公式2)參數(shù),來估算草甸生態(tài)系統(tǒng)的呼吸速率,本文利用已獲試驗(yàn)數(shù)據(jù)經(jīng)擬合確定Q10為2.63,從而建立起碳通量與土壤溫度的指數(shù)擬合關(guān)系(Reco=0.05eIn(2.63)(Ts-10)/10,R2=0.67)。

    4 結(jié)論

    1)科爾沁草甸生態(tài)系統(tǒng)碳通量在全年呈現(xiàn)“兩峰一谷”的變化規(guī)律,草甸生態(tài)系統(tǒng)碳吸收總量218.52 g CO2/m2。非生長季主要表現(xiàn)為碳源,凈CO2日排放速率均值0.11 mg/(m2·s),生長季主要表現(xiàn)為碳匯,凈CO2吸收速率均值-0.17 mg/(m2·s),較北美大草原的-0.88 mg/(m2·s)值偏?。?7],這可能與海拔及葉面積指數(shù)的差異有關(guān)[17]。NEE日變化呈明顯的“單峰型”,且日出、日落前后達(dá)到碳收支平衡狀態(tài)。7月12日,10月9日分別出現(xiàn)CO2日吸收(-0.89 mg/m2·s)和日釋放(0.56 mg/m2·s)的峰值。

    2)葉面積指數(shù)對草甸草原生長季NEE的變化影響明顯,LAI<0.97時,草甸CO2吸收量隨LAI的增加而線性增長,當(dāng)0.97<LAI<3.08時,CO2吸收量又隨LAI的增加而漸進(jìn)飽和,PAR與LAI的交互作用顯著(P<0.01),葉面積指數(shù)越大,凈碳交換NEE對PAR的響應(yīng)也越明顯。

    3)光合有效輻射也是影響NEE的重要環(huán)境因子。低光強(qiáng)下,生長季的NEE-PAR呈顯著線性相關(guān)(R2=0.924),生長季表觀量子效率均值為0.0015μmol CO2/μmol PAR,而在PAR>1200μmol/(m2·s)時,NEE對光強(qiáng)的響應(yīng)表現(xiàn)為漸進(jìn)飽和型;同時草甸生態(tài)系統(tǒng)的碳匯功能受到高VPD值的抑制;短暫強(qiáng)降雨增加了表層土壤水分含量,但削弱了到達(dá)地面的輻射強(qiáng)度,降低了晝間土壤呼吸速率,但對夜間NEE無明顯影響。而持續(xù)降雨使表層土壤含水量維持在較高水平,雖引起晝間碳交換能力的降低,但激發(fā)了夜間土壤呼吸速率的增加。

    [1] Rank A B,Dugas W A.Carbon dioxide fluxes over a northern,semiarid,mixed grass prairie.Agricultural and Forest Meteorology,2001,108:317-326.

    [2] Dugas W A,Heuer M L,Msyeux H S.Carbon dioxide fluxes over Bermudagrass,native prairie,and sorghum.Agricultural and Forest Meteorology,1999,93:121-139.

    [3] Shaver G R,Rastetter E B,Salmon V,et al.Pan-Arctic modeling of net ecosystem exchange of CO2.Philosophical Transactions of the Royal Society B-biolonical Sciences,2013,368:1-13.

    [4] Fu Y L,Yu G R,Wang Y F,et al.Impact of water stress on photosynthesis and respiration over Leymus chinensis ecosystem in Inner Mongolia.Science in China(Ser.D)Earth Sciences,2006,36(Supp,I):183-193.

    [5] Hao Y B,Wang Y F,Sun X M,et al.Seasonal variation in carbon exchange and its ecological analysis over Leymus chinensis steppe in Inner Mongolia.Science in China(Ser.D)Earth Sciences,2006,36:174-182.

    [6] Zhao Y P,Zhang X Z,Shi P L,et al.Impact of drought stress on net CO2exchange above an alpine grassland ecosystem in the Gentral Tibetan Plateau.Journal of Resources and Ecology,2013,4(4):327-336.

    [7] Li Y Q,Zhao H L,Zhao X Y,et al.Responses of soil respiration to air temperature from June to September in different sanddunes,Horqin sandy land.Journal of Desert Research,2008,28(2):249-254.

    [8] Chen N N,Guan D X,Jin C J,et al.Characteristics of soil respiration on Horqin meadow.Chinese Journal of Grassland,2011,33(5):82-87.

    [9] LI-COR Biosciences,Eddy Pro 4.0:Help and User’s Guide.Lincoln[EB/OL].[2013-12-02].http://envsupport.licor.com/help/EddyPro4/Default.htm.

    [10] Burba G G,McDermitt D K.Addressing the influence of instrument surface heat exchange on the measurements of CO2flux from open-path gas analyzers.Global Change Biology,2008,14:1854-1876.

    [11] Li Y J,Xu Z Z,Wang Y L,et al.Latent and sensible heat fluxes and energy balance in a maize agroecosystem.Chinese Journal of Plant Ecology,2007,31(6):1132-1144.

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

    [13] Raich J W,Schlesinger W H.The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate.Tellus,1992,44:81-99.

    [14] Xu L J,Tang H J,Yang G X,et al.Variation of net ecosystem carbon flux and its impact factors on Stipa baicalensis steppe in the growing season.Acta Prataculturae Sinica,2011,20(6):287-292.

    [15] Shi P L,Sun X M,Xu L L,et al.Net ecosystem CO2exchange and controlling factors in a Steppe-Kobresia meadow on Tibetan Plateau.Science in China(Ser.D),2006,36(Supp,I):194-203.

    [16] Suyker A E,Verma S B,Burba G G,et al.Gross primary production and ecosystem respiration of irrigated maize and irrigated soybean during a growing season.Agricultural and Forest Meteorology,2005,131(3/4):180-190.

    [17] Chen J L,Wang J,Wang Q,et al.LAI and light transmittance measurement of isolated Vitex shrub using LAI-2200 PCA.Chinese Journal of Agrometeorology,2014,35(4):373-379.

    [18] Xu L L,Zhang X Z,Shi P L,et al.Photoionization efficiency in response to temperature under low pressure of alpine meadow on Tibetan Platean.Science in China(Ser.D),2006,36(Supp,II):204-209.

    [19] Xu D Q,Li D Y,Shen Y G,et al.On midday depression of photosynthesis of wheat leaf under field conditions.Acta Phytophysiologia Sinica,1984,10:269-276.

    [20] Bunce J A.Afternoon inhibition of photosynthesis in maize.Field Crop Research,1990,24:261-271.

    [21] Fu Y L,Yu G R,Sun X M,et al.Depression of net ecosystem CO2exchange in semi-arid Leymus chinensis steppe and alpine shrub.Agricultural and Forest Meteorology,2006,137:234-244.

    [22] Lawerence B,F(xiàn)lanagan,Linda A,et al.Seasonal and interannual variation in carbon dioxide exchange and carbon balance in a northern temperate grassland.Global Change Biology,2002,8:599-615.

    [23] Zhang F W,Li Y N,Zhao X Q,et al.Effect of one precipitation process on CO2flux and themal transportation in alpine meadow of Qinghai-Tibetan plateau.Chinese Journal of Ecology,2008,27(10):1685-1691.

    [24] Zhang F W,Li Y N,Cao G M,et al.CO2fluxes and their driving factors over alpine meadow grassland ecosystems in the northern shore of Qinghai Lake,China.Chinese Journal of Plant Ecology,2012,36(3):187-198.

    [25] Wang W,Wang T,Peng S S,et al.Review of winter CO2efflux from soils:a key process of CO2exchange between soil and atmosphere.Chinese Journal of Plant Ecology,2007,31(3):394-402.

    [26] Chen J,Cao J J,Wei Y L,et al.Effect of grazing exclusion on soil respiration during the dormant in alpine meadow grassland ecosystems on the northern shore of Qinghai Lake,China.Acta Prataculturae Sinica,2014,23(6):78-86.

    [27] Verma S B,Kim J,Clement R J,et al.Water vapor and carbon dioxide exchange at a centrally located prairie site during FIFE.Journal of Geophysics Research,1992,97:18629-18639.

    參考文獻(xiàn):

    [4] 付玉玲,于貴瑞,王艷芬.水分脅迫對內(nèi)蒙古羊草草原光合和呼吸作用的影響.中國科學(xué)D輯,地球科學(xué),2006,36(增刊I):183-193.

    [5] 郝彥賓,王艷芬,孫曉敏,等.內(nèi)蒙古羊草草原碳交換季節(jié)變異及其生態(tài)學(xué)解析.中國科學(xué),2006,36:174-182.

    [7] 李玉強(qiáng),趙哈林,趙學(xué)勇,等.科爾沁沙地夏秋(6-9月)季不同類型沙丘土壤呼吸對氣溫變化的響應(yīng).中國沙漠,2008,28(2):249-254.

    [8] 陳妮娜,關(guān)德新,金昌杰,等.科爾沁草甸草地土壤呼吸特征.中國草地學(xué)報(bào),2011,33(5):82-87.

    [11] 李祎君,許振柱,王云龍,等.玉米農(nóng)田水熱通量動態(tài)與能量閉合分析.植物生態(tài)學(xué)報(bào),2007,31(6):1132-1144.

    [14] 徐麗君,唐華俊,楊桂霞,等.貝加爾針茅草原生態(tài)系統(tǒng)生長季碳通量及其影響因素分析.草業(yè)學(xué)報(bào),2011,20(6):287-292.

    [15] 石培禮,孫曉敏,徐玲玲,等.西藏高原草原化嵩草草甸生態(tài)系統(tǒng)CO2凈交換及其影響因子.中國科學(xué)D輯,地球科學(xué),2006,36(增刊I):194-203.

    [17] 陳景玲,王靜,王謙,等.植物冠層分析儀測定荊條孤立冠層葉面積指數(shù)和透光率的分析.中國農(nóng)業(yè)氣象,2014,35(4):373-379.

    [18] 徐玲玲,張憲洲,石培禮,等.青藏高原低氣壓條件下高寒草甸群落光量子效率對溫度的響應(yīng).中國科學(xué)D輯,地球科學(xué),2006,36(增刊I):204-209.

    [19] 許大全,李德耀,沈允剛,等.田間小麥葉片光合作用“午睡”現(xiàn)象的研究.植物生理學(xué)報(bào),1984,10:269-276.

    [23] 張法偉,李英年,趙新全,等.一次降水過程對青藏高原高寒草甸CO2通量和熱量輸送的影響.生態(tài)學(xué)雜志,2008,27(10):1685-1691.

    [24] 張法偉,李英年,曹廣民,等.青海湖北岸高寒草甸草原生態(tài)系統(tǒng)CO2通量特征及其驅(qū)動因子.植物生態(tài)學(xué)報(bào),2012,36(3):187-198.

    [25] 王娓,汪濤,彭書時,等.冬季土壤呼吸:不可忽視的地氣CO2交換過程.植物生態(tài)學(xué)報(bào),2007,31(3):394-402.

    [26] 陳驥,曹軍驥,魏永林,等.青海湖北岸高寒草甸草原非生長季土壤呼吸對溫度和濕度的響應(yīng).草業(yè)學(xué)報(bào),2014,23(6):78-86.

    Net ecosystem carbon exchange and its main driving factors in meadow grassland ecosystems in the Horqin Sand Lands

    WANG Jing1,LIU Ting-Xi1*,LEI Hui-Min2*,ZHANG Sheng-Wei1,DUAN Li-Min1
    1.Water Resources and Hydropower Survey and Design Institude,Inner Mongolia Agricultural University,Hohhot 010018,China;2.Department of Hydraulic Engineering,Tsinghua University,Beijing 100084,China

    The objectives of this study were to investigate the main drivers of ecosystem carbon dynamics in meadow grasslands in the Horqin Sand Lands,and to understand the relationship between net ecosystem CO2exchange(NEE)and key meteorological factors in this region.We continuously measured(NEE)and environmental factors in meadow grasslands in the Horqin Sand Lands from October 2013 to September 2014 using the eddy covariance method.The results showed that the diurnal variation in NEE differed significantly among seasons,especially among vegetation growing seasons.The highest seasonal diurnal absorption rate of CO2was in July,followed by August,September,June,and then May.The NEE was controlled by photosynthetically active radiation(PAR)and the leaf area index(LAI)during vegetation growing seasons.The relationship between NEE and PAR was expressed by a rectangular hyperbolic function.The simulated apparent initial average light-use efficiency was 0.0015μmol CO2/μmol PAR,and the maximum photosynthetic rate(Pmax)was 0.65μmol CO2/(m2·s).The NEE was influenced by LAI,which was expressed by piecewise functions.When the LAI was greater than 3.08,NEE-PAR showed a gradual saturation-type relationship.The higher the LAI,the stronger the response of net carbon exchange to the amount of PAR.A light saturation response occurred when 1.5<VPD<2 k Pa,and NEE was significantly inhibited under high VPD.The day-time net CO2exchange was inhibited by heavy rainfall(cumulative rainfall>40 mm/day),while sustained low-intensity rainfall(rainfall length>15 h)stimulated soil respiration at night.The night-time respiration rate of the ecosystem(Reco)increased exponentially as the soil temperature(Ts)increased,and the temperature-sensitive coefficient Q10was 2.63.

    net ecosystem exchange(NEE);meteorological factors;eddy covariance;Horqin sand

    10.11686/cyxb2014530 http://cyxb.lzu.edu.cn

    王婧,劉廷璽,雷慧閩,張圣微,段利民.科爾沁草甸生態(tài)系統(tǒng)凈碳交換特征及其驅(qū)動因子.草業(yè)學(xué)報(bào),2015,24(11):10-19.

    WANG Jing,LIU Ting-Xi,LEI Hui-Min,ZHANG Sheng-Wei,DUAN Li-Min.Net ecosystem carbon exchange and its main driving factors in meadow grassland ecosystems in the Horqin Sand Lands.Acta Prataculturae Sinica,2015,24(11):10-19.

    2014-12-19;改回日期:2015-03-10

    國家自然科學(xué)基金重點(diǎn)資助項(xiàng)目(51139002),教育部創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃(IRT13069)和內(nèi)蒙古農(nóng)業(yè)大學(xué)寒旱區(qū)水資源利用創(chuàng)新團(tuán)隊(duì)(NDTD2010-6)資助。

    王婧(1989-),女,內(nèi)蒙古阿拉善人,在讀碩士。E-mail:lebishu1234@163.com

    *通訊作者Corresponding author.E-mail:txliu1966@163.com,leihm@m(xù)ail.tsinghua.edu.cn

    猜你喜歡
    草甸葉面積通量
    高山草甸
    The world's narrowest river
    山地草甸
    冬小麥田N2O通量研究
    作物葉面積測量的研究進(jìn)展
    馬奶子葡萄葉面積評估模型的建立
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    “光合作用與細(xì)胞呼吸”重點(diǎn)分析和突破
    苧麻葉面積測定方法比較研究
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    99热6这里只有精品| 精品福利观看| 午夜福利在线在线| 国产欧美日韩精品亚洲av| 国产伦人伦偷精品视频| 亚洲成人久久爱视频| 欧洲精品卡2卡3卡4卡5卡区| 特级一级黄色大片| 波多野结衣高清无吗| 日韩欧美一区二区三区在线观看| 我要搜黄色片| 免费在线观看成人毛片| 99久久99久久久精品蜜桃| 身体一侧抽搐| 国产国拍精品亚洲av在线观看| 久久午夜福利片| 国产v大片淫在线免费观看| 午夜日韩欧美国产| 亚洲专区中文字幕在线| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 国产大屁股一区二区在线视频| 国产白丝娇喘喷水9色精品| 丁香欧美五月| 毛片一级片免费看久久久久 | 精品午夜福利视频在线观看一区| 黄色配什么色好看| 最新中文字幕久久久久| 最后的刺客免费高清国语| 欧美极品一区二区三区四区| 欧美最新免费一区二区三区 | 精品熟女少妇八av免费久了| 国产在线精品亚洲第一网站| 精品一区二区三区人妻视频| 天堂动漫精品| 午夜激情欧美在线| 国产亚洲精品av在线| 久久精品人妻少妇| 观看美女的网站| 97超级碰碰碰精品色视频在线观看| 免费一级毛片在线播放高清视频| 人人妻,人人澡人人爽秒播| 欧美日韩中文字幕国产精品一区二区三区| xxxwww97欧美| 色5月婷婷丁香| 99国产极品粉嫩在线观看| 色综合站精品国产| 欧美色欧美亚洲另类二区| 看黄色毛片网站| 欧美+亚洲+日韩+国产| 国产av不卡久久| 久9热在线精品视频| 天堂网av新在线| 我的老师免费观看完整版| 十八禁网站免费在线| 日本一二三区视频观看| 窝窝影院91人妻| 一进一出抽搐gif免费好疼| 国产高清视频在线观看网站| 51国产日韩欧美| 成人精品一区二区免费| 国产男靠女视频免费网站| 最近视频中文字幕2019在线8| 国产成人av教育| 小蜜桃在线观看免费完整版高清| 亚洲成人精品中文字幕电影| 欧美成人a在线观看| 直男gayav资源| 色5月婷婷丁香| 国内精品美女久久久久久| av在线蜜桃| 亚洲精品成人久久久久久| 性色av乱码一区二区三区2| 成人鲁丝片一二三区免费| 成人永久免费在线观看视频| 草草在线视频免费看| 欧美xxxx黑人xx丫x性爽| 综合色av麻豆| 看黄色毛片网站| 亚洲,欧美,日韩| 精品不卡国产一区二区三区| 亚洲专区中文字幕在线| 三级国产精品欧美在线观看| 亚洲真实伦在线观看| 国产主播在线观看一区二区| 在线免费观看的www视频| or卡值多少钱| 国产精品久久视频播放| 久久性视频一级片| 蜜桃亚洲精品一区二区三区| 欧美一级a爱片免费观看看| 国产午夜精品论理片| 中文亚洲av片在线观看爽| 亚洲真实伦在线观看| 久久久成人免费电影| 免费看美女性在线毛片视频| xxxwww97欧美| 小说图片视频综合网站| 在线免费观看的www视频| 日韩亚洲欧美综合| 人人妻人人看人人澡| 精品人妻1区二区| 免费大片18禁| www.熟女人妻精品国产| 国产高清视频在线播放一区| 我要搜黄色片| 直男gayav资源| 亚洲第一区二区三区不卡| 哪里可以看免费的av片| 日本撒尿小便嘘嘘汇集6| 国产探花极品一区二区| 国产私拍福利视频在线观看| 老司机午夜十八禁免费视频| 亚洲专区国产一区二区| 国产一区二区亚洲精品在线观看| 久久久国产成人精品二区| 欧美精品啪啪一区二区三区| 亚洲av免费在线观看| 99精品久久久久人妻精品| 18+在线观看网站| 色在线成人网| 日韩av在线大香蕉| 脱女人内裤的视频| 精品人妻视频免费看| 国产精品久久久久久人妻精品电影| 亚洲不卡免费看| 熟女电影av网| 国产黄a三级三级三级人| 中文字幕精品亚洲无线码一区| 搡老熟女国产l中国老女人| 久久精品国产自在天天线| 久久久久国内视频| 欧美高清成人免费视频www| 成人三级黄色视频| 日本一二三区视频观看| 国产成+人综合+亚洲专区| 色噜噜av男人的天堂激情| 久久久久九九精品影院| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆| 免费观看的影片在线观看| 可以在线观看毛片的网站| av天堂在线播放| 成人午夜高清在线视频| 99热精品在线国产| 一区福利在线观看| bbb黄色大片| 久久久久久久久久成人| 久久精品国产自在天天线| 亚洲自偷自拍三级| 老熟妇乱子伦视频在线观看| 欧美日韩综合久久久久久 | 少妇被粗大猛烈的视频| 国产激情偷乱视频一区二区| 亚洲成人久久爱视频| 国产久久久一区二区三区| 在线看三级毛片| 男女那种视频在线观看| 好男人电影高清在线观看| 少妇裸体淫交视频免费看高清| 亚洲不卡免费看| 12—13女人毛片做爰片一| 18禁在线播放成人免费| 亚洲人成网站高清观看| 丁香欧美五月| 熟女人妻精品中文字幕| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 国产精品伦人一区二区| 亚洲成人久久性| 欧美性猛交黑人性爽| a在线观看视频网站| 精品人妻熟女av久视频| 久久久久久久精品吃奶| 午夜福利欧美成人| av在线蜜桃| 亚洲中文日韩欧美视频| 观看免费一级毛片| 久久精品国产亚洲av天美| 看片在线看免费视频| 自拍偷自拍亚洲精品老妇| 久久天躁狠狠躁夜夜2o2o| 啦啦啦韩国在线观看视频| 亚洲精品日韩av片在线观看| 18禁黄网站禁片午夜丰满| 成人国产综合亚洲| 露出奶头的视频| 日韩欧美在线二视频| 国产高潮美女av| 国产精品精品国产色婷婷| 男插女下体视频免费在线播放| 日韩精品青青久久久久久| 久久久成人免费电影| 久久久久九九精品影院| 久久草成人影院| 国产大屁股一区二区在线视频| 成人欧美大片| 国产视频一区二区在线看| 亚洲专区中文字幕在线| 全区人妻精品视频| 天天一区二区日本电影三级| 亚洲中文日韩欧美视频| 小蜜桃在线观看免费完整版高清| 男女那种视频在线观看| 欧美国产日韩亚洲一区| 国产在视频线在精品| 日本一本二区三区精品| 欧美区成人在线视频| 色哟哟哟哟哟哟| 高清在线国产一区| 国产伦一二天堂av在线观看| 在线观看av片永久免费下载| 亚洲精品色激情综合| 97热精品久久久久久| 又爽又黄a免费视频| 又紧又爽又黄一区二区| 草草在线视频免费看| 欧美黑人巨大hd| 亚洲国产精品久久男人天堂| 白带黄色成豆腐渣| 男人舔女人下体高潮全视频| 嫩草影院新地址| 免费看日本二区| 日韩国内少妇激情av| www日本黄色视频网| 91在线精品国自产拍蜜月| 日本 欧美在线| 亚洲精品久久国产高清桃花| 精品久久久久久久久久久久久| 亚洲专区国产一区二区| 久久精品国产清高在天天线| 超碰av人人做人人爽久久| 免费在线观看成人毛片| 直男gayav资源| 久久6这里有精品| 亚洲精品久久国产高清桃花| 一个人观看的视频www高清免费观看| 久久久久久国产a免费观看| 噜噜噜噜噜久久久久久91| 亚洲狠狠婷婷综合久久图片| 精品久久国产蜜桃| а√天堂www在线а√下载| 一边摸一边抽搐一进一小说| 校园春色视频在线观看| 亚洲成人免费电影在线观看| 欧美成人性av电影在线观看| 国产熟女xx| 亚洲av免费在线观看| 精品午夜福利视频在线观看一区| 日本与韩国留学比较| 久久九九热精品免费| 黄色一级大片看看| 午夜亚洲福利在线播放| 91九色精品人成在线观看| 国产高清视频在线观看网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲狠狠婷婷综合久久图片| 国产av麻豆久久久久久久| 九九久久精品国产亚洲av麻豆| 桃色一区二区三区在线观看| 亚洲精品亚洲一区二区| 成人av在线播放网站| 99热这里只有是精品在线观看 | 禁无遮挡网站| 变态另类成人亚洲欧美熟女| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 亚洲av一区综合| 国产又黄又爽又无遮挡在线| 全区人妻精品视频| 一级a爱片免费观看的视频| 哪里可以看免费的av片| 欧美黑人欧美精品刺激| 亚洲熟妇中文字幕五十中出| 亚洲无线在线观看| 美女cb高潮喷水在线观看| 美女cb高潮喷水在线观看| 欧美乱妇无乱码| 人妻久久中文字幕网| 亚洲中文日韩欧美视频| 好男人电影高清在线观看| 欧美日韩乱码在线| 波多野结衣巨乳人妻| 国产成年人精品一区二区| 亚洲不卡免费看| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片| 国产高清视频在线播放一区| 精品欧美国产一区二区三| av福利片在线观看| 色精品久久人妻99蜜桃| 国产精品98久久久久久宅男小说| 一个人观看的视频www高清免费观看| 亚洲一区二区三区不卡视频| 18禁在线播放成人免费| 国产熟女xx| 女生性感内裤真人,穿戴方法视频| 国产又黄又爽又无遮挡在线| 国产精品,欧美在线| 特大巨黑吊av在线直播| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| 午夜影院日韩av| 少妇人妻精品综合一区二区 | 亚洲五月天丁香| 赤兔流量卡办理| 悠悠久久av| 久久人人爽人人爽人人片va | 亚洲七黄色美女视频| 久久性视频一级片| 国产中年淑女户外野战色| 国产成人啪精品午夜网站| 亚洲av成人精品一区久久| 波野结衣二区三区在线| 国产真实伦视频高清在线观看 | 亚洲av第一区精品v没综合| 精品99又大又爽又粗少妇毛片 | 成人国产一区最新在线观看| 最近中文字幕高清免费大全6 | 在线看三级毛片| 亚洲经典国产精华液单 | 欧美日韩亚洲国产一区二区在线观看| 夜夜看夜夜爽夜夜摸| 91麻豆精品激情在线观看国产| 国产在视频线在精品| 国产视频一区二区在线看| 我的女老师完整版在线观看| 欧美乱妇无乱码| 亚洲精品456在线播放app | 午夜激情欧美在线| 日本在线视频免费播放| 欧美高清成人免费视频www| 久久亚洲精品不卡| 国产欧美日韩精品亚洲av| 国产一区二区三区在线臀色熟女| 在线观看美女被高潮喷水网站 | 看黄色毛片网站| 午夜福利欧美成人| 在现免费观看毛片| 日日夜夜操网爽| 人人妻,人人澡人人爽秒播| 一级毛片久久久久久久久女| 丰满人妻熟妇乱又伦精品不卡| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 国产成人啪精品午夜网站| 成人毛片a级毛片在线播放| 性色av乱码一区二区三区2| 久久伊人香网站| 欧美一区二区精品小视频在线| 丰满乱子伦码专区| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| av视频在线观看入口| 欧美日韩亚洲国产一区二区在线观看| 变态另类成人亚洲欧美熟女| av在线蜜桃| 国产高潮美女av| 久久国产乱子伦精品免费另类| 国产伦精品一区二区三区视频9| 国产一区二区亚洲精品在线观看| 嫩草影院新地址| av女优亚洲男人天堂| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 精品久久久久久久久久免费视频| 美女黄网站色视频| 欧美一级a爱片免费观看看| 高清日韩中文字幕在线| 国产成人a区在线观看| 亚洲18禁久久av| 国产色婷婷99| 女生性感内裤真人,穿戴方法视频| 成人欧美大片| 噜噜噜噜噜久久久久久91| 99久久99久久久精品蜜桃| 我要看日韩黄色一级片| 一个人看视频在线观看www免费| 免费看美女性在线毛片视频| 两性午夜刺激爽爽歪歪视频在线观看| 两个人视频免费观看高清| 少妇熟女aⅴ在线视频| 一区二区三区四区激情视频 | 日韩欧美精品v在线| 国产伦人伦偷精品视频| 一本久久中文字幕| 美女 人体艺术 gogo| 国产aⅴ精品一区二区三区波| 午夜免费男女啪啪视频观看 | 99久久精品国产亚洲精品| 国产亚洲精品综合一区在线观看| 97超视频在线观看视频| 精品久久久久久久久久免费视频| 亚洲成人中文字幕在线播放| 欧美高清性xxxxhd video| 国产免费av片在线观看野外av| 一本精品99久久精品77| 桃色一区二区三区在线观看| 久久久久亚洲av毛片大全| 亚洲成a人片在线一区二区| 亚洲国产精品sss在线观看| 亚洲av不卡在线观看| 国产精品av视频在线免费观看| 99久久成人亚洲精品观看| 国产91精品成人一区二区三区| 国产精品一区二区免费欧美| 床上黄色一级片| 国产激情偷乱视频一区二区| 亚洲av一区综合| 欧美激情久久久久久爽电影| 97热精品久久久久久| 亚洲精品在线美女| 18禁黄网站禁片午夜丰满| 精品人妻1区二区| 老司机午夜福利在线观看视频| 在线观看一区二区三区| 九色国产91popny在线| 久久精品影院6| 中文字幕高清在线视频| 免费电影在线观看免费观看| 很黄的视频免费| 又紧又爽又黄一区二区| 国产主播在线观看一区二区| 精品一区二区三区视频在线| 99热这里只有是精品50| 两个人的视频大全免费| 久久人人爽人人爽人人片va | 99久久成人亚洲精品观看| 欧美区成人在线视频| 免费看美女性在线毛片视频| 级片在线观看| 婷婷色综合大香蕉| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品色激情综合| 免费一级毛片在线播放高清视频| 免费观看人在逋| 国产黄片美女视频| 最近在线观看免费完整版| 日韩精品青青久久久久久| 国产单亲对白刺激| 他把我摸到了高潮在线观看| 又爽又黄a免费视频| 久久人人精品亚洲av| 日韩亚洲欧美综合| 波野结衣二区三区在线| 搡老妇女老女人老熟妇| 窝窝影院91人妻| 我的女老师完整版在线观看| 在线看三级毛片| 亚洲av五月六月丁香网| 一进一出好大好爽视频| 国产亚洲av嫩草精品影院| 成人一区二区视频在线观看| 成人国产一区最新在线观看| 欧美一区二区亚洲| 久久精品国产亚洲av天美| 啪啪无遮挡十八禁网站| 美女被艹到高潮喷水动态| 久久久久久久久久黄片| av专区在线播放| 欧美日韩乱码在线| 久久香蕉精品热| 亚洲人成电影免费在线| 一级a爱片免费观看的视频| 精品久久久久久久久久久久久| 长腿黑丝高跟| 欧美绝顶高潮抽搐喷水| 国产成人影院久久av| 日韩欧美 国产精品| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产| 美女被艹到高潮喷水动态| 欧美日韩乱码在线| 日韩中文字幕欧美一区二区| 亚洲人成伊人成综合网2020| 亚洲美女视频黄频| 国产av麻豆久久久久久久| 一个人免费在线观看电影| 小说图片视频综合网站| 在线免费观看不下载黄p国产 | 我要搜黄色片| 夜夜看夜夜爽夜夜摸| 中文字幕精品亚洲无线码一区| 3wmmmm亚洲av在线观看| 亚洲欧美日韩高清在线视频| 观看美女的网站| 欧美日韩综合久久久久久 | 欧美另类亚洲清纯唯美| 长腿黑丝高跟| 亚洲欧美日韩东京热| 国产亚洲av嫩草精品影院| 亚洲 国产 在线| 精品无人区乱码1区二区| 国产高潮美女av| 99riav亚洲国产免费| 婷婷精品国产亚洲av| 成年女人永久免费观看视频| 亚洲久久久久久中文字幕| 永久网站在线| 一个人免费在线观看电影| 国产成人a区在线观看| 最新中文字幕久久久久| 三级毛片av免费| 亚洲国产高清在线一区二区三| 天美传媒精品一区二区| 日韩欧美在线乱码| 亚洲自偷自拍三级| 午夜精品一区二区三区免费看| 亚洲熟妇中文字幕五十中出| 熟妇人妻久久中文字幕3abv| 嫁个100分男人电影在线观看| 久久久久亚洲av毛片大全| 久久精品国产亚洲av天美| 亚洲av成人av| 国产精品乱码一区二三区的特点| 日韩亚洲欧美综合| 五月伊人婷婷丁香| 免费人成在线观看视频色| 狠狠狠狠99中文字幕| 亚洲精品在线美女| 久久精品国产99精品国产亚洲性色| 91麻豆精品激情在线观看国产| 18+在线观看网站| 国产精品久久视频播放| 国产精品影院久久| 中文在线观看免费www的网站| 露出奶头的视频| 在线a可以看的网站| 一级黄片播放器| 他把我摸到了高潮在线观看| 欧美日本视频| 午夜老司机福利剧场| 欧美又色又爽又黄视频| 国产大屁股一区二区在线视频| 国产精品精品国产色婷婷| 国产亚洲精品久久久久久毛片| 99热这里只有是精品在线观看 | 天堂av国产一区二区熟女人妻| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 黄色丝袜av网址大全| 淫秽高清视频在线观看| 国产精品亚洲一级av第二区| 美女高潮的动态| 国产精品,欧美在线| 国产精品亚洲美女久久久| 亚洲最大成人手机在线| www日本黄色视频网| 欧美潮喷喷水| 色精品久久人妻99蜜桃| 亚洲真实伦在线观看| 国产色爽女视频免费观看| 观看美女的网站| 色哟哟哟哟哟哟| 午夜福利成人在线免费观看| 有码 亚洲区| 黄色日韩在线| 国内精品久久久久久久电影| 黄片小视频在线播放| 内射极品少妇av片p| 久9热在线精品视频| 国内揄拍国产精品人妻在线| 午夜日韩欧美国产| 日韩欧美精品免费久久 | 中文在线观看免费www的网站| 欧美性感艳星| 精品一区二区三区人妻视频| 亚洲avbb在线观看| 九九热线精品视视频播放| 小说图片视频综合网站| 首页视频小说图片口味搜索| 长腿黑丝高跟| 亚洲经典国产精华液单 | 色综合亚洲欧美另类图片| 搡老熟女国产l中国老女人| 成年人黄色毛片网站| 色噜噜av男人的天堂激情| 看黄色毛片网站| 亚洲五月婷婷丁香| 一级a爱片免费观看的视频| 国产精品野战在线观看| 欧美xxxx黑人xx丫x性爽| 午夜福利18| 啦啦啦韩国在线观看视频| 久久天躁狠狠躁夜夜2o2o| 久久人人精品亚洲av| 脱女人内裤的视频| 精品久久久久久久末码| 村上凉子中文字幕在线| 天堂网av新在线| 欧美成人a在线观看| 国产野战对白在线观看| 国内精品久久久久精免费| 99国产精品一区二区蜜桃av| 国产精品影院久久| 天堂网av新在线| 亚洲五月天丁香| 精品久久久久久,| 窝窝影院91人妻| 99riav亚洲国产免费| 精品久久久久久久久亚洲 | av天堂在线播放| 国内精品久久久久精免费| 色精品久久人妻99蜜桃| 国产真实伦视频高清在线观看 | 久久欧美精品欧美久久欧美| 亚洲av成人精品一区久久| 久久亚洲真实| 亚洲成人精品中文字幕电影| 丁香六月欧美| 美女cb高潮喷水在线观看| 午夜a级毛片| 91九色精品人成在线观看| 麻豆成人午夜福利视频| 久久99热这里只有精品18|