• <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)呼吸各組分間的相互關系

    2013-12-09 06:01:21朱先進于貴瑞王秋鳳高艷妮趙新全韓士杰閆俊華
    生態(tài)學報 2013年21期
    關鍵詞:異養(yǎng)渦度通量

    朱先進,于貴瑞,王秋鳳,高艷妮,趙新全,韓士杰,閆俊華

    (1. 中國科學院地理科學與資源研究所生態(tài)網(wǎng)絡觀測與模擬重點實驗室CERN綜合研究中心, 北京 100101; 2. 中國科學院大學, 北京 100049; 3. 中國科學院西北高原生物研究所, 西寧 810001; 4. 中國科學院沈陽應用生態(tài)研究所, 沈陽 110016; 5. 中國科學院華南植物園, 廣州 510650)

    典型森林和草地生態(tài)系統(tǒng)呼吸各組分間的相互關系

    朱先進1,2,于貴瑞1,*,王秋鳳1,高艷妮1,2,趙新全3,韓士杰4,閆俊華5

    (1. 中國科學院地理科學與資源研究所生態(tài)網(wǎng)絡觀測與模擬重點實驗室CERN綜合研究中心, 北京 100101; 2. 中國科學院大學, 北京 100049; 3. 中國科學院西北高原生物研究所, 西寧 810001; 4. 中國科學院沈陽應用生態(tài)研究所, 沈陽 110016; 5. 中國科學院華南植物園, 廣州 510650)

    生態(tài)系統(tǒng)呼吸是陸地生態(tài)系統(tǒng)碳收支的重要組成部分,分析其組分間的相互關系對理解生態(tài)系統(tǒng)呼吸過程和精確評價生態(tài)系統(tǒng)碳收支具有重要意義,也是當前碳循環(huán)研究工作的一大難點。利用ChinaFLUX的長白山溫帶針闊混交林(CBS),鼎湖山亞熱帶常綠闊葉林(DHS)和海北灌叢草甸(HBGC)3個典型生態(tài)系統(tǒng)的通量觀測數(shù)據(jù),分析了經驗統(tǒng)計方法在中國典型生態(tài)系統(tǒng)中的適用性及敏感性,揭示了生態(tài)系統(tǒng)呼吸各組分的動態(tài)變化特征及相互關系。結果表明:采用呼吸組分拆分方法所獲結果與理論推測及實測數(shù)據(jù)大致相同,拆分結果對凈初級生產力與總初級生產力比值(NPP/GPP)的變化較為敏感,NPP/GPP變化0.1時,自養(yǎng)呼吸在生態(tài)系統(tǒng)呼吸中的比例(Ra/RE)改變0.05。各生態(tài)系統(tǒng)中,生態(tài)系統(tǒng)呼吸及其組分在年內均表現(xiàn)出明顯的單峰型變化特征,在夏季生長旺盛的時節(jié)達到最大值。異養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值(Rh/RE)也具有明顯的季節(jié)變化,但在生態(tài)系統(tǒng)間表現(xiàn)出明顯差異,CBS和HBGC分別表現(xiàn)出先增大后減小和先減小后增大的變化趨勢,DHS則相對穩(wěn)定,在0.5附近波動, Ra/RE的季節(jié)動態(tài)與Rh/RE相反。在年總量上,HBGC主要通過異養(yǎng)呼吸向大氣排放CO2,異養(yǎng)呼吸占生態(tài)系統(tǒng)呼吸的60%,而CBS和DHS的自養(yǎng)呼吸和異養(yǎng)呼吸所占比重大致相似,異養(yǎng)呼吸占生態(tài)系統(tǒng)呼吸的49%。這說明,經驗統(tǒng)計學模型可以用來進行生態(tài)系統(tǒng)呼吸組分的拆分,進而可以為生態(tài)系統(tǒng)碳循環(huán)過程的精細研究提供參考數(shù)據(jù),但今后應加強NPP/GPP的測定,以提高生態(tài)系統(tǒng)呼吸拆分的精度。

    陸地生態(tài)系統(tǒng);生態(tài)系統(tǒng)呼吸;渦度相關;碳通量;自養(yǎng)呼吸;異養(yǎng)呼吸

    生態(tài)系統(tǒng)呼吸(RE)是碳循環(huán)的重要分量[1],與氣候變化間存在著正反饋效應[2- 3],不僅在當前全球碳收支中占有較大的比重,而且在未來還會呈現(xiàn)不斷增大的趨勢[4]。已有研究表明,全球的RE約為103 PgC a-1[5],占全球總初級生產力(GPP)總量(123 PgC a-1[6])的83.74%。

    按照呼吸底物的不同,RE可分為自養(yǎng)呼吸(Ra)和異養(yǎng)呼吸(Rh)。其中,Ra是植物維持自身生命活動所必需的新陳代謝過程,又分為維持呼吸和生長呼吸。維持呼吸與溫度相關,隨著溫度的增加呈指數(shù)增加[7];生長呼吸與溫度沒有直接關系,但與GPP呈一定的比例[8]。Rh是殘存有機質分解并向大氣釋放CO2的過程[9],主要受土壤溫濕度和有機質含量的影響[10]。因而,影響RE組分的因素存在差異,拆分RE的各個組分、分析各組分間的相互關系有助于深入認識生態(tài)系統(tǒng)碳循環(huán)的過程機理,也可為陸地生態(tài)系統(tǒng)碳收支的精確評估提供理論依據(jù)。同時,Ra的變化對大氣CO2濃度具有重要影響[11],科學界迫切需要拆分和評估RE的各組分,為評價Ra對大氣CO2濃度變化的貢獻提供數(shù)據(jù)。

    目前,拆分RE組分的方法有許多種。第1種是將箱式法與渦度相關系統(tǒng)觀測相結合,分別獲取土壤呼吸(Rs)和RE,實現(xiàn)RE組分的拆分[12- 14]。然而,基于箱式法獲取的Rs包含了Ra和Rh,對Rs進行拆分也是一項艱巨的工作。第2種是利用生物量調查法獲取生態(tài)系統(tǒng)凈初級生產力(NPP),結合渦度相關觀測結果估算Ra和Rh[15- 17],但生物量調查和渦度相關觀測的時間尺度不一致,并且生物量調查法耗費大量的人力物力,難以獲得持續(xù)的、高質量觀測數(shù)據(jù)。第3種是基于過程機理模型模擬Ra和Rh,但該方法所需參數(shù)較多,結果也受到模型設計者對RE過程理解程度的限制[10]。基于同位素分餾比也可以實現(xiàn)對RE組分的拆分[18- 19],但所需儀器設備非常昂貴。Schwalm等[20]提出了基于通量觀測數(shù)據(jù)的統(tǒng)計學拆分方法,并對其不確定性進行了分析。結果表明,該方法具有很好的穩(wěn)定性,不同參數(shù)方案間沒有明顯差異,這為利用通量觀測數(shù)據(jù)分析RE組分提供了新的技術途徑。

    中國陸地生態(tài)系統(tǒng)通量觀測網(wǎng)絡(ChinaFLUX)自2002年建立以來,已積累了大量的觀測數(shù)據(jù),構建了學術界普遍認可的通量數(shù)據(jù)處理流程[21]。然而,現(xiàn)有通量數(shù)據(jù)處理僅將生態(tài)系統(tǒng)碳通量分解為GPP和RE,尚沒有拆分RE的組分。因此,本研究利用Schwalm等[20]分析方法,拆分并評價我國典型森林和草地生態(tài)系統(tǒng)的RE組分及其相互關系,進而闡明:(1)不同生態(tài)系統(tǒng)Ra和Rh的動態(tài)特征,(2)Rh/RE的年內動態(tài)變化特征,(3)不同生態(tài)系統(tǒng)間Ra和Rh年總量的差異。該研究可以為生態(tài)系統(tǒng)碳循環(huán)過程研究提供基礎數(shù)據(jù),也為通量數(shù)據(jù)處理過程中RE拆分方法的確定提供依據(jù)。

    1 研究地點與研究方法

    1.1 站點介紹

    本研究選擇ChinaFLUX中3個陸地生態(tài)系統(tǒng)為研究對象,分別為:長白山溫帶針闊混交林(CBS)、鼎湖山亞熱帶常綠闊葉林(DHS)和海北高寒灌叢草甸(HBGC),各生態(tài)系統(tǒng)的基本信息如表1所示。

    表1 研究站點基本信息

    CBS:長白山溫帶針闊混交林Changbai Mountains temperate conifer and broadleaf mixed forest;HBGC:海北灌叢草甸Haibei alpine Shrub;DHS:鼎湖山亞熱帶常綠闊葉林Dinghu Mountains subtropical evergreen broadleaf forest

    3個生態(tài)系統(tǒng)均利用開路式渦度相關系統(tǒng)(OPEC)進行生態(tài)系統(tǒng)碳水通量觀測,原始采樣頻率為10 Hz,通過數(shù)據(jù)采集器計算并儲存30 min的通量數(shù)據(jù)。同時觀測溫度、濕度、降水、輻射等氣象要素并計算30 min平均值儲存[21- 24]。此外,輔以靜態(tài)箱-氣相色譜法測定土壤呼吸,每隔7—10 d測定1次[25- 27]。

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

    1.2.1 通量數(shù)據(jù)處理

    采用ChinaFLUX通用數(shù)據(jù)處理流程對30 min通量數(shù)據(jù)進行處理[21],包括三次坐標旋轉、WPL校正、儲存項計算及降水剔除、閾值剔除和方差剔除等,并根據(jù)Reichstein等[28]的方法確定夜間u*臨界值、實現(xiàn)夜間數(shù)據(jù)剔除、完成數(shù)據(jù)質量控制。利用非線性回歸方法對碳通量數(shù)據(jù)進行插補,并將碳通量拆分為GPP和RE。

    1.2.2 自養(yǎng)呼吸和異養(yǎng)呼吸的拆分

    本研究根據(jù)自養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值(Ra/RE)對RE進行拆分[20]。具體做法為:假設凈生態(tài)系統(tǒng)生產力(NEP)與RE的比值為θ,Ra/RE為β,則:

    (1)

    因而,根據(jù)式(1)可以得到β,即:

    (2)

    已有研究表明,NPP/GPP介于0.47—0.60之間[20],因而,β的范圍可以表示為:

    (3)

    在生態(tài)學研究中,Schwalm等[20]將Ra/RE或Rh/RE的最小值Ω設為0.10,并分析了Ω上下浮動0.05后的結果,他發(fā)現(xiàn)基于不同Ω得到的β沒有明顯差異,因而此處Ω設為0.10。

    當β取最小值(0.10)時,θ=-0.75;當β取最大值(0.90)即Rh/RE取最小值(0.10)時,θ=0.70。因而,θlt;-0.75或者θgt;0.70時,公式(3)不再適用。當θ=-0.75時,βlower為0.10,βupper為0.13;當θ=0.70時,βlower為0.68,βupper為0.90。針對不同的θ值,設定不同的β范圍,具體范圍如表2所示。

    表2 不同θ值下的β范圍

    根據(jù)θ值,在β給定的范圍內隨機選取1000個數(shù)據(jù),取其平均值為β。夜間時,θ=-1,利用公式(3)無法獲取β值,因而將當日白天最后一個有效的β和次日最早一個有效的β進行線性內插的方法獲取夜間β。

    最后,根據(jù)β值和RE計算Ra和Rh,即:

    Ra=β×RE
    Rh=(1-β)×RE

    (4)

    1.3 生態(tài)系統(tǒng)呼吸拆分結果的驗證

    為了驗證該方法的合理性,利用靜態(tài)箱-氣相色譜法觀測的土壤呼吸(Rs)數(shù)據(jù)與本研究獲取的Rh進行對比。由于靜態(tài)箱-氣相色譜法的觀測時間是9:00—11:00,因此,Rh選用9:00—11:30的平均值。

    1.4 數(shù)據(jù)統(tǒng)計

    采用單因素方差分析法評價不同生態(tài)系統(tǒng)間RE組分的差異,采用線性回歸分析拆分結果與靜態(tài)箱-氣相色譜法所獲結果的差異。所有數(shù)據(jù)統(tǒng)計分析過程均在matlab7.7中進行。

    2 結果與分析

    2.1 生態(tài)系統(tǒng)呼吸及其組分的季節(jié)動態(tài)

    生態(tài)系統(tǒng)碳通量具有明顯的季節(jié)變化特征[29- 31]。由于年際間生態(tài)系統(tǒng)的季節(jié)動態(tài)特征相似[32- 34],本文以2003年數(shù)據(jù)為例分析中國典型森林和草地生態(tài)系統(tǒng)碳通量的季節(jié)動態(tài)(圖1)。

    圖1 2003年典型生態(tài)系統(tǒng)氣溫及碳通量的季節(jié)變化特征Fig.1 Seasonal dynamics of air temperature and carbon fluxes in typical ecosystems in 2003

    氣溫(Ta)在3個典型生態(tài)系統(tǒng)均表現(xiàn)出單峰型變化(圖1),最大值出現(xiàn)在7月份前后,但不同生態(tài)系統(tǒng),其年內氣溫的變化幅度存在明顯差異,DHS常年保持較高的氣溫,北方兩個生態(tài)系統(tǒng)(CBS和HBGC)年內氣溫變幅較大。不同生態(tài)系統(tǒng),GPP的季節(jié)動態(tài)有所差異(圖1),北方生態(tài)系統(tǒng)(CBS和HBGC)表現(xiàn)出明顯的單峰型特點,但DHS的GPP在秋冬季節(jié)較高,在6—8月因溫度過高而略低,這與已有結果[23]相一致。NEP的季節(jié)動態(tài)與GPP相似(圖1),但峰值出現(xiàn)時間較GPP略有提前。

    3個生態(tài)系統(tǒng)的RE均表現(xiàn)出明顯的單峰型變化,峰值均出現(xiàn)在7—8月份(圖1),與氣溫的季節(jié)動態(tài)(圖1)相似,這是因為,生態(tài)系統(tǒng)尺度上,氣溫的波動決定了RE的季節(jié)動態(tài)[35]。生態(tài)系統(tǒng)間,最大呼吸速率存在差異,CBS可以達到9.20 gC m-2d-1,DHS和HBGC明顯小于CBS,均不足4 gC m-2d-1。這主要是生產力的差異引起的,生產力與RE存在著密切的耦合關系[35- 37]。

    與RE的動態(tài)特征相似,Ra和Rh也表現(xiàn)出單峰型的季節(jié)動態(tài)(圖1),峰值出現(xiàn)的時間也與RE相似,峰值大小約為RE峰值的一半。

    2.2 生態(tài)系統(tǒng)呼吸組分與生態(tài)系統(tǒng)呼吸比值(Rh/RE和Ra/RE)的季節(jié)動態(tài)

    異養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值(Rh/RE)表現(xiàn)出明顯的季節(jié)波動,但在生態(tài)系統(tǒng)間表現(xiàn)出不同的變化特征(圖2)。

    Rh/RE在CBS表現(xiàn)為年初和年末略低、生長季開始和結束時略有增加的趨勢,在生長旺盛時節(jié)大致維持在0.5左右,DHS全年均在0.5附近波動。與兩個森林生態(tài)系統(tǒng)的變化趨勢不同,Rh/RE在HBGC表現(xiàn)出明顯的單峰型動態(tài),非生長季可以達到0.9,隨著生長季的開始,Rh/RE逐漸降低,并在生長旺季達到0.5左右的最小值,隨著生長季的結束,異養(yǎng)呼吸占據(jù)主導,Rh/RE又開始增大。

    自養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值(Ra/RE)表現(xiàn)出與Rh/RE相反的季節(jié)變化趨勢。

    圖2 2003年三個典型生態(tài)系統(tǒng)的生態(tài)系統(tǒng)呼吸組分在生態(tài)系統(tǒng)呼吸中所占比例(Rh/RE, Ra/RE)的季節(jié)動態(tài)特征Fig.2 Seasonal dynamics of the ratio of autotrophic respiration and heterotrophic respiration to ecosystem respiration (Ra/RE, Rh/RE) in typical ecosystems in 2003

    2.3 生態(tài)系統(tǒng)呼吸及其組分的年總量

    生態(tài)系統(tǒng)呼吸的組分反映了生態(tài)系統(tǒng)消耗碳量的不同途徑。中國3個典型生態(tài)系統(tǒng)的RE及其組分的年總量如表3所示。

    生態(tài)系統(tǒng)間,RE、Ra和Rh均存在顯著差異,均表現(xiàn)為CBS最高,DHS居中,HBGC最低。生態(tài)系統(tǒng)間Rh/RE的差異表現(xiàn)為:兩個森林生態(tài)系統(tǒng)(CBS和DHS)間沒有明顯差別,約為0.49,但顯著小于HBGC的0.60。Ra/RE表現(xiàn)出相反的特征,即HBGC的數(shù)值為0.40,明顯小于CBS和DHS的0.51。

    表3 典型生態(tài)系統(tǒng)的生態(tài)系統(tǒng)呼吸及其組分的年總量

    表中數(shù)據(jù)為2003—2005年的平均值,括號內數(shù)值為2003—2005年的標準差

    3 討論

    3.1 生態(tài)系統(tǒng)呼吸拆分結果的不確定性

    本研究采用的方法是Schwalm等[20]對加拿大森林RE進行研究時提出的,目前尚沒有在其他生態(tài)系統(tǒng)中應用。但相對于其他研究方法,本方法結構簡單、可操作性強并具有一定的生態(tài)學理論基礎。

    理論上來說,當NPP/GPP保持恒定的比值時,Ra/RE隨著NEP/RE的增加而增大,即隨著NEP/RE的增加,RE中的Rh組分降低,從圖2可知,Ra/RE的變化與NEP的變化有相似的規(guī)律,這證實了該方法在理論上是可行的。

    在理論分析的基礎上,用靜態(tài)箱-氣相色譜法觀測的Rs與本研究中拆分得到的Rh作對比(圖3),以驗證該方法的拆分效果。

    從圖3可以看出,在CBS,Rh與Rs的比值為0.34(圖3),這在之前研究結果[38- 39]的范圍內,并與東北林區(qū)其他研究[40]相似。在DHS和HBGC,Rh與Rs間沒有明顯的關系。這是因為,在CBS,利用箱式法觀測的Rs是生態(tài)系統(tǒng)呼吸的一部分(圖3),可以反映RE的動態(tài),但在DHS和HBGC,箱式法觀測到的Rs幾乎與該時段的RE相一致,這可能是由箱式法本身觀測的誤差引起的。

    進一步整理了已發(fā)表文獻的NPP數(shù)據(jù),結合渦度相關觀測的碳通量結果,估算RE及其組分年總量間的相互關系。在CBS,張娜等[41]模擬發(fā)現(xiàn),1995年該生態(tài)系統(tǒng)Ra/RE高達0.82,高于本研究結果(表3),這可能是他們估算的NEP過高及NPP偏小導致的,NEP值高達392 gC m-2a-1[41],比渦度相關觀測結果[34]高54%。在DHS,2003—2004年的NPP為970 gC m-2a-1[42],結合Yu等[34]的結果可以看到,這兩年該生態(tài)系統(tǒng)的Ra/RE為52.33%,與本研究結果(表3)相似。

    雖然Piao等[11]總結已有研究結果[1, 43- 44]發(fā)現(xiàn),Ra與GPP的比值雖然表現(xiàn)出隨著溫度的增加先降低后增大,但具有較大的變異性,基于本研究的方法,CBS、DHS、HBGC的Ra/GPP分別為0.40、0.31、0.35,也在Piao等[11]的范圍之內。

    因而,本研究所用方法具有簡便易行、無需額外測定即可準確拆分RE,這為分析碳收支組分及其影響因素提供了可能,可以在生態(tài)系統(tǒng)碳循環(huán)研究中應用。

    3.2 生態(tài)系統(tǒng)呼吸拆分結果對NPP/GPP的敏感性

    本研究將NPP/GPP的值設為0.47—0.60,但Litton等[43]整理文獻數(shù)據(jù)發(fā)現(xiàn),NPP/GPP變異很大,最小甚至可能低于0.20。因而有必要進一步分析NPP/GPP的變化對RE組分拆分的影響。

    此處設定兩個方案,即將NPP/GPP(θ)的上下限(0.47—0.60)分別下調0.10(方案I)和上調0.10(方案II),相應的RE/NEP(β)的臨界比值發(fā)生改變(表4)。

    表4 不同處理方案下θ值下的β范圍

    基于不同數(shù)據(jù)處理方案獲得了生態(tài)系統(tǒng)呼吸組分相互關系的季節(jié)動態(tài)(圖4),Ra/RE及Rh/RE的季節(jié)動態(tài)與圖2沒有明顯差異,說明不同NPP/GPP的變化不會對Ra/RE及Rh/RE的季節(jié)動態(tài)產生明顯影響。

    圖4 2003年不同NPP/GPP的典型生態(tài)系統(tǒng)的生態(tài)系統(tǒng)呼吸組分在生態(tài)系統(tǒng)呼吸中所占比例(Rh/RE, Ra/RE)的季節(jié)動態(tài)特征Fig.4 Seasonal dynamics of the ratio of autotrophic respiration and heterotrophic respiration to ecosystem respiration (Ra/RE, Rh/RE) in typical ecosystems in 2003 with different NPP/GPP

    不同數(shù)據(jù)處理方案下,生態(tài)系統(tǒng)呼吸組分在生態(tài)系統(tǒng)呼吸中所占的比例發(fā)生一定的波動(表5),NPP/GPP每改變0.1,Ra/RE變化0.05,Ra/RE隨著NPP/GPP的增加而降低。

    表5不同數(shù)據(jù)處理方案生態(tài)系統(tǒng)呼吸組分在生態(tài)系統(tǒng)呼吸中所占比例

    Table5Theratioofautotrophicrespirationandheterotrophicrespirationtoecosystemrespiration(Ra/RE, Rh/RE)underdifferentdatatreatments

    處理Treatment自養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值Ra/RE長白山CBS鼎湖山DHS海北灌叢HBGC異養(yǎng)呼吸與生態(tài)系統(tǒng)呼吸的比值Rh/RE長白山CBS鼎湖山DHS海北灌叢HBGC本研究Thisstudy0.51(0.01)0.51(0.01)0.40(0.00)0.49(0.01)0.49(0.01)0.60(0.01)方案ⅠTreatmentⅠ0.56(0.01)0.55(0.01)0.43(0.01)0.44(0.01)0.45(0.02)0.57(0.01)方案ⅡTreatmentⅡ0.46(0.01)0.45(0.01)0.36(0.00)0.54(0.01)0.55(0.02)0.64(0.00)

    表中數(shù)據(jù)為2003—2005年的平均值,括號內數(shù)值為觀測期間的標準差

    生態(tài)系統(tǒng)間,不管NPP/GPP發(fā)生何種變化,HBGC的Ra/RE均小于CBS和DHS,并且該值始終小于0.5,因而HBGC的生態(tài)系統(tǒng)呼吸以異養(yǎng)呼吸為主。隨著NPP/GPP的變化,CBS和DHS中Ra/RE在0.5上下浮動,表明森林生態(tài)系統(tǒng)中自養(yǎng)和異養(yǎng)呼吸所占的比例大致相似。

    可見,NPP/GPP影響著RE組分拆分,生物量調查獲取的NPP將為生態(tài)系統(tǒng)碳通量組分的拆分提供重要基礎數(shù)據(jù)。

    3.3 生態(tài)系統(tǒng)呼吸各組分間的相互關系

    在季節(jié)動態(tài)上,Ra/RE隨著NEP的增加有增大的趨勢[20],表明隨著植物生長的開始,生態(tài)系統(tǒng)在不斷固碳的基礎上通過自養(yǎng)呼吸向外排放的碳量也在增加。這是因為Ra中的維持呼吸與溫度相關,隨著溫度的增加呈指數(shù)增加[7];生長呼吸與溫度沒有直接關系,但與GPP呈一定的比例關系[8];而Rh受土壤溫濕度和有機質含量的影響。因此,在生長季開始后,隨著GPP和氣溫的增大,Ra增加速度快于Rh,導致Ra/RE開始增加。

    生態(tài)系統(tǒng)之間,Ra年總量與RE年總量的比值也明顯不同,森林生態(tài)系統(tǒng)明顯高于草地生態(tài)系統(tǒng),表明森林生態(tài)系統(tǒng)Ra所占比例高于草地生態(tài)系統(tǒng),與Schwalm等[20]的結果相同。這是因為,森林生態(tài)系統(tǒng)具有較高的NEP,NEP/RE也較高,導致Ra/RE高。比較中國與加拿大森林生態(tài)系統(tǒng)間[20]的Ra/RE的差異可以發(fā)現(xiàn),中國森林生態(tài)系統(tǒng)的Ra/RE明顯高于加拿大森林生態(tài)系統(tǒng),這也是由中國森林生態(tài)系統(tǒng)中NEP/RE高于加拿大所決定的。

    4 結論

    基于ChinaFLUX3個陸地生態(tài)系統(tǒng)(CBS:長白山溫帶混交林,DHS:鼎湖山亞熱帶常綠闊葉林,HBGC:海北高寒灌叢草甸)的渦度相關觀測數(shù)據(jù),利用統(tǒng)計分析方法拆分了RE組分,分析了該方法在中國陸地生態(tài)系統(tǒng)的適用性及對NPP/GPP的敏感性,闡明了RE組分的動態(tài)特征及其相互關系。結果表明:基于渦度相關觀測數(shù)據(jù)及NPP/GPP的比值,可以準確拆分中國典型生態(tài)系統(tǒng)的RE,但該方法對NPP/GPP較為敏感。3個生態(tài)系統(tǒng)中,RE及其組分年內均表現(xiàn)出單峰型的變化,在生長旺盛時節(jié)(7月份)達到最大值。異養(yǎng)呼吸在生態(tài)系統(tǒng)呼吸中所占的比例(Rh/RE)表現(xiàn)出明顯的季節(jié)變化,CBS表現(xiàn)出先增大后減小的季節(jié)特征,DHS則常年維持為0.5左右,HBGC表現(xiàn)出先降低后增大的特點。年總量上,Rh/RE因生態(tài)系統(tǒng)的不同而不同,CBS和DHS兩個森林生態(tài)系統(tǒng)的自養(yǎng)和異養(yǎng)呼吸大致相當,HBGC的碳排放則主要通過異養(yǎng)呼吸來實現(xiàn)。該方法可以拆分RE組分,進而為生態(tài)系統(tǒng)碳循環(huán)的過程研究提供參考數(shù)據(jù),但今后應加強NPP/GPP的測定,以提高生態(tài)系統(tǒng)呼吸拆分的精度。

    Referrences:

    [1] Luyssaert S, Inglima I, Jung M, Richardson A D, Reichsteins M, Papale D, Piao S L, Schulzes E D, Wingate L, Matteucci G, Aragao L, Aubinet M, Beers C, Bernhoffer C, Black K G, Bonal D, Bonnefond J M, Chambers J, Ciais P, Cook B, Davis K J, Dolman A J, Gielen B, Goulden M, Grace J, Granier A, Grelle A, Griffis T, Grunwald T, Guidolotti G, Hanson P J, Harding R, Hollinger D Y, Hutyra L R, Kolar P, Kruijt B, Kutsch W, Lagergren F, Laurila T, Law B E, Le Maire G, Lindroth A, Loustau D, Malhi Y, Mateus J, Migliavacca M, Misson L, Montagnani L, Moncrieff J, Moors E, Munger J W, Nikinmaa E, Ollinger S V, Pita G, Rebmann C, Roupsard O, Saigusa N, Sanz M J, Seufert G, Sierra C, Smith M L, Tang J, Valentini R, Vesala T, Janssens I A. CO2balance of boreal, temperate, and tropical forests derived from a global database. Global Change Biology, 2007, 13(12): 2509- 2537.

    [2] Davidson E A, Janssens I A, Luo Y Q. On the variability of respiration in terrestrial ecosystems: moving beyond Q10. Global Change Biology, 2006, 12(2): 154- 164.

    [3] Houghton R A, Davidson E A, Woodwell G M. Missing sinks, feedbacks, and understanding the role of terrestrial ecosystems in the global carbon balance. Global Biogeochemical Cycles, 1998, 12(1): 25- 34.

    [4] Cox P M, Betts R A, Jones C D, Spall S A, Totterdell I J. Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model. Nature, 2000, 408(6809): 184- 187.

    [5] Yuan W P, Luo Y Q, Li X L, Liu S G, Yu G R, Zhou T, Bahn M, Black A, Desai A R, Cescatti A, Marcolla B, Jacobs C, Chen J Q, Aurela M, Bernhofer C, Gielen B, Bohrer G, Cook D R, Dragoni D, Dunn A L, Gianelle D, Grünwald T, Ibrom A, Leclerc M Y, Lindroth A, Liu H P, Marchesini L B, Montagnani L, Pita G, Rodeghiero M, Rodrigues A, Starr G, Stoy P C. Redefinition and global estimation of basal ecosystem respiration rate. Global Biogeochemical Cycles, 2011, 25: GB4002, doi: 10.1029/2011gb004150.

    [6] Beer C, Reichstein M, Tomelleri E, Ciais P, Jung M, Carvalhais N, Rodenbeck C, Arain M A, Baldocchi D, Bonan G B, Bondeau A, Cescatti A, Lasslop G, Lindroth A, Lomas M, Luyssaert S, Margolis H, Oleson K W, Roupsard O, Veenendaal E, Viovy N, Williams C, Woodward F I, Papale D. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science, 2010, 329(5993): 834- 838.

    [7] Ryan M G. A simple method for estimating gross carbon budgets for vegetation in forest ecosystems. Tree Physiology, 1991, 9(1/2): 255- 266.

    [8] Chen J M, Liu J, Cihlar J, Goulden M L. Daily canopy photosynthesis model through temporal and spatial scaling for remote sensing applications. Ecological Modelling, 1999, 124(2/3): 99- 119.

    [9] Trumbore S. Carbon respired by terrestrial ecosystems-recent progress and challenges. Global Change Biology, 2006, 12(2): 141- 153.

    [10] Yu G R. Scientific Frontier on Human Activities and Ecosystem Changes. Beijing: Higher Education Press, 2009.

    [11] Piao S L, Luyssaert S, Ciais P, Janssens I A, Chen A P, Cao C, Fang J Y, Friedlingstein P, Luo Y Q, Wang S P. Forest annual carbon cost: a global-scale analysis of autotrophic respiration. Ecology, 2010, 91(3): 652- 661.

    [12] Irvine J, Law B E, Martin J G, Vickers D. Interannual variation in soil CO2efflux and the response of root respiration to climate and canopy gas exchange in mature ponderosa pine. Global Change Biology, 2008, 14(12): 2848- 2859.

    [13] Nagy Z, Pintér K, Pavelka M, Darenová E, Balogh J. Carbon fluxes of surfaces vs. ecosystems: advantages of measuring eddy covariance and soil respiration simultaneously in dry grassland ecosystems. Biogeosciences, 2011, 8(9): 2523- 2534.

    [14] Griffis T J, Black T A, Gaumont-Guay D, Drewitt G B, Nesic Z, Barr A G, Morgenstern K, Kljun N. Seasonal variation and partitioning of ecosystem respiration in a southern boreal aspen forest. Agricultural and Forest Meteorology, 2004, 125(3/4): 207- 223.

    [15] Law B E, Thornton P E, Irvine J, Anthoni P M, Van Tuyl S. Carbon storage and fluxes in ponderosa pine forests at different developmental stages. Global Change Biology, 2001, 7(7): 755- 777.

    [16] Luyssaert S, Reichstein M, Schulze E D, Janssens I A, Law B E, Papale D, Dragoni D, Goulden M L, Granier A, Kutsch W L, Linder S, Matteucci G, Moors E, Munger J W, Pilegaard K, Saunders M, Falge E M. Toward a consistency cross-check of eddy covariance flux-based and biometric estimates of ecosystem carbon balance. Global Biogeochemical Cycles, 2009, 23: GB3009, doi: 10.1029/2008gb003377.

    [17] Tan Z H, Zhang Y P, Yu G R, Sha L Q, Tang J W, Deng X B, Song Q H. Carbon balance of a primary tropical seasonal rain forest. Journal of Geophysical Research, 2010, 115: D00H26, doi: 10.1029/2009jd012913.

    [18] Riveros-Iregui D A, Hu J, Burns S P, Bowling D R, Monson R K. An interannual assessment of the relationship between the stable carbon isotopic composition of ecosystem respiration and climate in a high-elevation subalpine forest. Journal of Geophysical Research, 2011, 116: G02005, doi: 10.1029/2010jg001556.

    [19] Schuur E A G, Trumbore S E. Partitioning sources of soil respiration in boreal black spruce forest using radiocarbon. Global Change Biology, 2006, 12(2): 165- 176.

    [20] Schwalm C R, Black T A, Morgenstern K, Humphreys E R. A method for deriving net primary productivity and component respiratory fluxes from tower-based eddy covariance data: a case study using a 17-year data record from a Douglas-fir chronosequence. Global Change Biology, 2007, 13(2): 370- 385.

    [21] Yu G R, Fu Y L, Sun X M, Wen X F, Zhang L M. Recent progress and future directions of ChinaFLUX. Science in China Series D: Earth Sciences, 2006, 49(Suppl II): 1- 23.

    [22] Zhang J H, Yu G R, Han S J, Guan D X, Sun X M. Seasonal and annual variation of CO2flux above a broad-leaved Korean pine mixed forest. Science in China Series D: Earth Sciences, 2006, 49(Suppl II): 63- 73.

    [23] Wang C L, Yu G R, Zhou G Y, Yan J H, Zhang L M, Wang X, Tang X L, Sun X M. CO2flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan, China. Science in China Series D: Earth Sciences, 2006, 49(Suppl II): 127- 138.

    [24] Li Y N, Sun X M, Zhao X Q, Zhao L, Xu S X, Gu S, Zhang G F, Yu G R. Seasonal variations and mechanism for environmental control of NEE of CO2concerning thePotentillafruticosain alpine shrub meadow of Qinghai-Tibet Plateau. Science in China Series D: Earth Sciences, 2006, 49(Suppl II): 174- 185.

    [25] Zheng Z M, Yu G R, Sun X M, Cao G M, Wang Y S, Du M Y, Li J, Li Y N. Comparison of eddy covariance and static chamber/gas chromatogram methods in measuring ecosystem respiration. Chinese Journal of Applied Ecology, 2008, 19(2): 290- 298.

    [26] Lin L S, Han S J, Wang Y S, Gu Z J. Soil CO2flux in several typical forests of Mt. Changbai. Chinese Journal of Ecology, 2004, 23(5): 42- 45.

    [27] Zhang D Q, Sun X M, Zhou G Y, Yan J H, Wang Y S, Liu S Z, Zhou C Y, Liu J X, Tang X L, Li J, Zhang Q M. Seasonal dynamics of soil CO2effluxes with responses to environmental factors in lower subtropical forests of China. Science in China Series D: Earth Sciences, 2006, 49(Suppl II): 139- 149.

    [28] Reichstein M, Falge E, Baldocchi D, Papale D, Aubinet M, Berbigier P, Bernhofer C, Buchmann N, Gilmanov T, Granier A, Grunwald T, Havrankova K, Ilvesniemi H, Janous D, Knohl A, Laurila T, Lohila A, Loustau D, Matteucci G, Meyers T, Miglietta F, Ourcival J M, Pumpanen J, Rambal S, Rotenberg E, Sanz M, Tenhunen J, Seufert G, Vaccari F, Vesala T, Yakir D, Valentini R. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Global Change Biology, 2005, 11(9): 1424- 1439.

    [29] Law B E, Williams M, Anthoni P M, Baldocchi D D, Unsworth M H. Measuring and modelling seasonal variation of carbon dioxide and water vapour exchange of aPinusponderosaforest subject to soil water deficit. Global Change Biology, 2000, 6(6): 613- 630.

    [30] Valentini R, Gamon J A, Field C B. Ecosystem ags-exchange in a California grassland-seasonal patterns and implications for scaling. Ecology, 1995, 76(6): 1940- 1952.

    [31] Guan D X, Wu J B, Zhao X S, Han S J, Yu G R, Sun X M, Jin C J. CO2fluxes over an old, temperate mixed forest in northeastern China. Agricultural and Forest Meteorology, 2006, 137(3/4): 138- 149.

    [32] Barr A G, Black T A, Hogg E H, Griffis T J, Morgenstern K, Kljun N, Theede A, Nesic Z. Climatic controls on the carbon and water balances of a boreal aspen forest, 1994—2003. Global Change Biology, 2007, 13(3): 561- 576.

    [33] Polley H W, Frank A B, Sanabria J, Phillips R L. Interannual variability in carbon dioxide fluxes and flux-climate relationships on grazed and ungrazed northern mixed-grass prairie. Global Change Biology, 2008, 14(7): 1620- 1632.

    [34] Yu G R, Zhang L M, Sun X M, Fu Y L, Wen X F, Wang Q F, Li S G, Ren C Y, Song X, Liu Y F, Han S J, Yan J H. Environmental controls over carbon exchange of three forest ecosystems in eastern China. Global Change Biology, 2008, 14(11): 2555- 2571.

    [35] Wen X F, Wang H M, Wang J L, Yu G R, Sun X M. Ecosystem carbon exchanges of a subtropical evergreen coniferous plantation subjected to seasonal drought, 2003—2007. Biogeosciences, 2010, 7(1): 357- 369.

    [36] Wang X C, Wang C K, Yu G R. Spatio-temporal patterns of forest carbon dioxide exchange based on global eddy covariance measurements. Science in China Series D: Earth Sciences, 2008, 51(8): 1129- 1143.

    [37] Lasslop G, Reichstein M, Detto M, Richardson A D, Baldocchi D D. Comment on Vickers et al.: self-correlation between assimilation and respiration resulting from flux partitioning of eddy-covariance CO2fluxes. Agricultural and Forest Meteorology, 2010, 150(2): 312- 314.

    [38] Subke J A, Inglima I, Francesca Cotrufo M. Trends and methodological impacts in soil CO2efflux partitioning: a metaanalytical review. Global Change Biology, 2006, 12(6): 921- 943.

    [39] Zhan X Y, Yu G R, Zheng Z M, Wang Q F. Carbon emission and spatial pattern of soil respiration of terrestrial ecosystems in China: based on geostatistic estimation of flux measurement. Progress in Geography, 2012, 31(1): 97- 108.

    [40] Wang C K, Yang J Y. Rhizospheric and heterotrophic components of soil respiration in six Chinese temperate forests. Global Change Biology, 2007, 13(1): 123- 131.

    [41] Zhang N, Yu G R, Zhao S D, Yu Z L. Carbon budget of ecosystem in Changbai Mountain natural reserve. Environmental Science, 2003, 24(1): 24- 32.

    [42] Yan J H, Wang Y P, Zhou G Y, Zhang D Q. Estimates of soil respiration and net primary production of three forests at different succession stages in South China. Global Change Biology, 2006, 12(5): 810- 821.

    [43] Litton C M, Raich J W, Ryan M G. Carbon allocation in forest ecosystems. Global Change Biology, 2007, 13(10): 2089- 2109.

    [44] DeLucia E H, Drake J E, Thomas R B, Gonzalez-Meler M. Forest carbon use efficiency: is respiration a constant fraction of gross primary production?. Global Change Biology, 2007, 13(6): 1157- 1167.

    [10] 于貴瑞. 人類活動與生態(tài)系統(tǒng)變化的前沿科學問題. 北京: 高等教育出版社, 2009.

    [25] 鄭澤梅, 于貴瑞, 孫曉敏, 曹廣民, 王躍思, 杜明遠, 李俊, 李英年. 渦度相關法和靜態(tài)箱/氣相色譜法在生態(tài)系統(tǒng)呼吸觀測中的比較. 應用生態(tài)學報, 2008, 19(2): 290- 298.

    [26] 林麗莎, 韓士杰, 王躍思, 谷志靜. 長白山四種林分土壤CO2釋放通量的研究. 生態(tài)學雜志, 2004, 23(5): 42- 45.

    [39] 展小云, 于貴瑞, 鄭澤梅, 王秋鳳. 中國區(qū)域陸地生態(tài)系統(tǒng)土壤呼吸的碳排放及空間格局——基于通量觀測的地學統(tǒng)計評估. 地理科學進展, 2012, 31(1): 97- 108.

    [41] 張娜, 于貴瑞, 趙士洞, 于振良. 長白山自然保護區(qū)生態(tài)系統(tǒng)碳平衡研究. 環(huán)境科學, 2003, 24(1): 24- 32.

    Theinteractionbetweencomponentsofecosystemrespirationintypicalforestandgrasslandecosystems

    ZHU Xianjin1,2, YU Guirui1, *, WANG Qiufeng1, GAO Yanni1,2, ZHAO Xinquan3, HAN Shijie4, YAN Junhua5

    1SynthesisResearchCenterofCERN,KeyLaboratoryofEcosystemNetworkObservationandModeling,InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China2UniversityofChineseAcademyofSciences,Beijing100049,China3NorthwestPlateauInstituteofBiology,ChineseAcademyofSciences,Xining810001,China4InstituteofAppliedEcology,ChineseAcademyofSciences,Shenyang110016,China5SouthChinaBotanicalGarden,ChineseAcademyofSciences,Guangzhou510650,China

    Ecosystem respiration (RE), the release of carbon from an ecosystem into the atmosphere, is a key component of the terrestrial ecosystem carbon budget and plays an important role in the global carbon balance. Analyzing the interaction between RE components is essential to understand RE and to accurately evaluate the ecosystem carbon budget. There are many methods for separating RE into autotrophic respiration (Ra) and heterotrophic respiration (Rh), but each approach has disadvantages. Large RE data have obtained through long-term eddy covariance measurements, while the interaction between Ra and Rh is poorly documented, which inhibits the accurate assessment of global carbon budget. In this study, we used an empirical statistical method to separate RE into its two components and to examine component relationships and seasonal dynamics at three ChinaFLUX sites: 1) Changbaishan temperate mixed forest (CBS); 2) Dinghushan subtropical evergreen broad-leaf forest (DHS); and 3) Haibei shrub meadow (HBGC). The applicability and sensitivity of this method in typical ecosystems of China were also evaluated. The method used in this study was based on the ratio of Ra to RE (Ra/RE). The range of Ra/RE was obtained by calculating two ratios: the ratio of RE to net ecosystem productivity (NEP) (RE/NEP) and that of net primary productivity (NPP) to gross primary productivity (GPP) (NPP/GPP). Within the range of Ra/RE, 1000 Ra/REs were randomly selected and the value of Ra/RE used in this study was set as the mean of the 1000 random Ra/REs. Ra and Rh were then calculated using Ra/RE and RE.

    Our study showed that the RE separating method produced consistent results with those obtained through static chamber/gas chromatographic techniques at the same sites, as well as with biomass surveys and theoretical speculation. The interaction of RE components was sensitive to the variation of NPP/GPP: a ten-percent increase of NPP/GPP led to a five-percent decrease of Ra/RE. In all three ecosystems, RE and its components showed similar seasonal dynamics, with a single-peak pattern achieving its maximum midway through the growing season. The ratio of Rh to RE (Rh/RE) also showed different seasonal dynamic among the three ecosystems. In CBS, Rh/RE increased during the first half of the year, reached its peak during the growing season then decreased. However, Rh/RE in HBGC decreased during the first half of the year and increased again later in growing-season. In DHS, Rh/RE was relatively stable at approximately 0.5. Seasonal dynamics of Ra/RE were opposite to those in Rh/RE. The annual total Rh accounted for 60% of the RE in HBGC, suggesting that a large proportion of emitted carbon was released by Rh in this ecosystem. In CBS and DHS, Rh was only 49% of RE, indicating that the release of carbon through Ra and Rh was nearly the same in these two forest ecosystems. Results indicate that this statistical method, which requires detailed observations of NPP/GPP, can successfully separate RE into Rh and Ra and can provide necessary data for the detailed analysis of the ecosystem carbon cycle.

    terrestrial ecosystem; ecosystem respiration; eddy covariance; carbon fluxes; autotrophic respiration; heterotrophic respiration

    國家重點基礎研究發(fā)展規(guī)劃資助項目(2010CB833504); 國家自然科學基金資助項目(31061140359, 30590380); 中國科學院戰(zhàn)略性先導科技專項資助項目(XDA05050601)

    2012- 07- 13;

    2013- 06- 21

    *通訊作者Corresponding author.E-mail: yugr@igsnrr.ac.cn

    10.5846/stxb201207130988

    朱先進,于貴瑞,王秋鳳,高艷妮,趙新全,韓士杰,閆俊華.典型森林和草地生態(tài)系統(tǒng)呼吸各組分間的相互關系.生態(tài)學報,2013,33(21):6925- 6934.

    Zhu X J, Yu G R, Wang Q F, Gao Y N, Zhao X Q, Han S J, Yan J H.The interaction between components of ecosystem respiration in typical forest and grassland ecosystems.Acta Ecologica Sinica,2013,33(21):6925- 6934.

    猜你喜歡
    異養(yǎng)渦度通量
    冬小麥田N2O通量研究
    異養(yǎng)同化降解氯代烴的研究現(xiàn)狀、微生物代謝特性及展望
    南海相對渦度的時空變化特征分析*
    斜壓渦度的變化與臺風暴雨的關系研究
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    化工進展(2015年6期)2015-11-13 00:26:29
    馬尾松種子生長與幼苗成長過程的養(yǎng)分動態(tài)研究
    一次東北南部特大暴雪的動力機制分析
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    異養(yǎng)小球藻的篩選鑒定及其特性研究
    久久久久视频综合| 国产高清国产精品国产三级| 午夜91福利影院| 国产无遮挡羞羞视频在线观看| 一级爰片在线观看| 国产免费现黄频在线看| 久久这里只有精品19| 青春草亚洲视频在线观看| 欧美国产精品va在线观看不卡| 欧美成人精品欧美一级黄| 国产精品久久久人人做人人爽| 成年美女黄网站色视频大全免费| 成人国产av品久久久| 国产精品国产三级国产专区5o| 一个人免费看片子| 午夜福利免费观看在线| 一区二区av电影网| 王馨瑶露胸无遮挡在线观看| 男人舔女人的私密视频| 人妻人人澡人人爽人人| 亚洲精品国产av蜜桃| 亚洲,欧美精品.| 国产精品 国内视频| 亚洲欧美成人综合另类久久久| 9色porny在线观看| 日韩一区二区三区影片| 欧美黄色片欧美黄色片| 日本色播在线视频| 黑人欧美特级aaaaaa片| 久久久久久免费高清国产稀缺| 久久久欧美国产精品| 高清黄色对白视频在线免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 一二三四中文在线观看免费高清| 日韩制服丝袜自拍偷拍| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 久久精品久久精品一区二区三区| 国产无遮挡羞羞视频在线观看| 免费看不卡的av| 美女大奶头黄色视频| 久久久久国产一级毛片高清牌| 国产免费一区二区三区四区乱码| 女人高潮潮喷娇喘18禁视频| 性高湖久久久久久久久免费观看| 婷婷色综合大香蕉| 丝袜美足系列| 国产亚洲av片在线观看秒播厂| 看十八女毛片水多多多| 一级毛片 在线播放| 色播在线永久视频| 久久精品国产a三级三级三级| 欧美激情极品国产一区二区三区| 婷婷色麻豆天堂久久| 一区二区av电影网| 精品一区在线观看国产| 国产精品免费大片| 国产精品亚洲av一区麻豆 | 国产一级毛片在线| 久久人人97超碰香蕉20202| 晚上一个人看的免费电影| 狠狠精品人妻久久久久久综合| 亚洲国产毛片av蜜桃av| 久久av网站| 日韩熟女老妇一区二区性免费视频| 午夜91福利影院| 一本大道久久a久久精品| 最近2019中文字幕mv第一页| 免费日韩欧美在线观看| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久精品电影小说| 婷婷成人精品国产| 精品亚洲成a人片在线观看| 高清在线视频一区二区三区| 欧美av亚洲av综合av国产av | 操美女的视频在线观看| 男女边吃奶边做爰视频| 一区二区三区四区激情视频| 91精品国产国语对白视频| 久久精品亚洲熟妇少妇任你| 亚洲欧美精品综合一区二区三区| 超碰97精品在线观看| 波多野结衣av一区二区av| 亚洲成人免费av在线播放| 欧美人与善性xxx| 国产免费现黄频在线看| 街头女战士在线观看网站| 国产男人的电影天堂91| 两个人看的免费小视频| 卡戴珊不雅视频在线播放| 国产精品一区二区在线不卡| 国产精品 国内视频| 大片电影免费在线观看免费| 99热网站在线观看| 国产成人av激情在线播放| 亚洲自偷自拍图片 自拍| 王馨瑶露胸无遮挡在线观看| 亚洲成国产人片在线观看| 国产成人91sexporn| 亚洲国产欧美在线一区| 国产成人a∨麻豆精品| av电影中文网址| 在现免费观看毛片| 国产成人免费无遮挡视频| 青春草国产在线视频| 啦啦啦啦在线视频资源| av线在线观看网站| 在现免费观看毛片| 热re99久久国产66热| 青春草国产在线视频| 新久久久久国产一级毛片| 免费少妇av软件| 亚洲男人天堂网一区| 丝袜美腿诱惑在线| 免费高清在线观看视频在线观看| 国产成人欧美在线观看 | 狠狠婷婷综合久久久久久88av| 大话2 男鬼变身卡| 别揉我奶头~嗯~啊~动态视频 | 满18在线观看网站| 精品国产一区二区久久| 男女国产视频网站| 日本vs欧美在线观看视频| 国产欧美日韩综合在线一区二区| 亚洲欧洲日产国产| 色94色欧美一区二区| 欧美黄色片欧美黄色片| 久久精品国产亚洲av涩爱| 国产熟女欧美一区二区| 亚洲av国产av综合av卡| 国产97色在线日韩免费| 久久久精品94久久精品| 日韩免费高清中文字幕av| 国产精品免费视频内射| 亚洲av在线观看美女高潮| 久久天堂一区二区三区四区| videosex国产| 一级,二级,三级黄色视频| 国产97色在线日韩免费| 老司机影院成人| 国产午夜精品一二区理论片| 亚洲国产av影院在线观看| 老司机亚洲免费影院| 精品一品国产午夜福利视频| 看非洲黑人一级黄片| 人妻人人澡人人爽人人| 亚洲,欧美,日韩| 日韩精品有码人妻一区| 国产精品久久久久久精品电影小说| 亚洲第一青青草原| 中国三级夫妇交换| 婷婷成人精品国产| 亚洲精品aⅴ在线观看| 久久午夜综合久久蜜桃| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕| 卡戴珊不雅视频在线播放| 欧美精品高潮呻吟av久久| a级毛片在线看网站| 看免费av毛片| 久久精品熟女亚洲av麻豆精品| 欧美亚洲 丝袜 人妻 在线| 啦啦啦中文免费视频观看日本| 日本91视频免费播放| 欧美精品人与动牲交sv欧美| 黑人欧美特级aaaaaa片| 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 亚洲,一卡二卡三卡| 男女午夜视频在线观看| 一区二区三区精品91| 亚洲精品久久午夜乱码| 激情视频va一区二区三区| 欧美中文综合在线视频| 19禁男女啪啪无遮挡网站| 91精品三级在线观看| a级毛片在线看网站| 最近最新中文字幕免费大全7| 人人澡人人妻人| 少妇 在线观看| 国产成人午夜福利电影在线观看| 国产欧美亚洲国产| 19禁男女啪啪无遮挡网站| 免费人妻精品一区二区三区视频| 18禁裸乳无遮挡动漫免费视频| 一级毛片 在线播放| 啦啦啦在线免费观看视频4| 人人妻人人添人人爽欧美一区卜| 精品久久久精品久久久| 精品亚洲成a人片在线观看| 久久女婷五月综合色啪小说| 中文字幕人妻熟女乱码| 啦啦啦在线免费观看视频4| 九色亚洲精品在线播放| 亚洲一级一片aⅴ在线观看| 精品一区在线观看国产| 精品国产超薄肉色丝袜足j| 熟女av电影| 制服诱惑二区| svipshipincom国产片| 国产乱来视频区| 少妇被粗大的猛进出69影院| 国产免费福利视频在线观看| 久久久久精品久久久久真实原创| 在线观看免费高清a一片| 丝瓜视频免费看黄片| 色网站视频免费| 又粗又硬又长又爽又黄的视频| 老司机影院成人| 美女中出高潮动态图| 中文天堂在线官网| 中文字幕最新亚洲高清| 国语对白做爰xxxⅹ性视频网站| 老汉色∧v一级毛片| 亚洲成av片中文字幕在线观看| 国产精品无大码| 久久人人爽av亚洲精品天堂| 老司机靠b影院| 亚洲中文av在线| 日韩伦理黄色片| svipshipincom国产片| 午夜福利影视在线免费观看| 精品亚洲成国产av| 啦啦啦 在线观看视频| 日韩大码丰满熟妇| 操美女的视频在线观看| 国产免费视频播放在线视频| 麻豆精品久久久久久蜜桃| 欧美黑人精品巨大| 中文字幕人妻丝袜制服| 亚洲av男天堂| 最近手机中文字幕大全| 丝瓜视频免费看黄片| 欧美最新免费一区二区三区| 69精品国产乱码久久久| 国精品久久久久久国模美| 在线观看免费视频网站a站| 视频区图区小说| 日韩 亚洲 欧美在线| 免费高清在线观看日韩| 国精品久久久久久国模美| 一区二区日韩欧美中文字幕| 最新在线观看一区二区三区 | 精品一区二区免费观看| 建设人人有责人人尽责人人享有的| 免费观看人在逋| 黑人欧美特级aaaaaa片| 欧美变态另类bdsm刘玥| 亚洲成国产人片在线观看| 国产成人欧美| 国产黄色免费在线视频| 日本欧美国产在线视频| 国产精品av久久久久免费| 成人漫画全彩无遮挡| av在线老鸭窝| 波多野结衣一区麻豆| 亚洲国产成人一精品久久久| 人成视频在线观看免费观看| 91精品伊人久久大香线蕉| 亚洲人成网站在线观看播放| 久久久久国产一级毛片高清牌| 在线亚洲精品国产二区图片欧美| 免费高清在线观看视频在线观看| 国产黄色视频一区二区在线观看| 久久99一区二区三区| 亚洲成人手机| 日本黄色日本黄色录像| 精品午夜福利在线看| 美女国产高潮福利片在线看| 国产亚洲精品第一综合不卡| 久久精品久久久久久久性| 久久亚洲国产成人精品v| 欧美日韩成人在线一区二区| 成人18禁高潮啪啪吃奶动态图| 欧美 日韩 精品 国产| 久久精品aⅴ一区二区三区四区| 亚洲第一区二区三区不卡| e午夜精品久久久久久久| av卡一久久| 在线观看人妻少妇| 51午夜福利影视在线观看| 大片电影免费在线观看免费| 九色亚洲精品在线播放| 少妇被粗大的猛进出69影院| 女人高潮潮喷娇喘18禁视频| 国产伦理片在线播放av一区| 国产 一区精品| 亚洲国产精品一区三区| 男人爽女人下面视频在线观看| 日本欧美国产在线视频| 午夜福利一区二区在线看| 伦理电影大哥的女人| 日日撸夜夜添| 卡戴珊不雅视频在线播放| 国产精品久久久av美女十八| 一级片免费观看大全| 男女国产视频网站| 黄网站色视频无遮挡免费观看| 久久久欧美国产精品| 青草久久国产| 黄色毛片三级朝国网站| 欧美成人午夜精品| 亚洲av日韩精品久久久久久密 | 亚洲国产精品999| 午夜福利乱码中文字幕| 夫妻性生交免费视频一级片| 亚洲欧洲国产日韩| 啦啦啦在线免费观看视频4| 国产xxxxx性猛交| 国产 一区精品| 男男h啪啪无遮挡| 免费少妇av软件| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 久久精品国产a三级三级三级| 亚洲精品美女久久久久99蜜臀 | 日韩欧美一区视频在线观看| 亚洲精品第二区| 日日啪夜夜爽| 国产高清不卡午夜福利| 欧美最新免费一区二区三区| 中文精品一卡2卡3卡4更新| 欧美 亚洲 国产 日韩一| 亚洲欧美精品综合一区二区三区| 丝瓜视频免费看黄片| 在线观看人妻少妇| 亚洲四区av| 韩国精品一区二区三区| 国产精品二区激情视频| 成人国产麻豆网| 亚洲婷婷狠狠爱综合网| 日韩精品免费视频一区二区三区| 狠狠婷婷综合久久久久久88av| 久久ye,这里只有精品| 麻豆精品久久久久久蜜桃| 久久女婷五月综合色啪小说| 亚洲av福利一区| 国产日韩欧美视频二区| 国产一区二区激情短视频 | 男人舔女人的私密视频| 日本一区二区免费在线视频| 在线亚洲精品国产二区图片欧美| 免费在线观看黄色视频的| 建设人人有责人人尽责人人享有的| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 看十八女毛片水多多多| netflix在线观看网站| 日韩熟女老妇一区二区性免费视频| 校园人妻丝袜中文字幕| 免费在线观看完整版高清| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 黄色怎么调成土黄色| 1024香蕉在线观看| 人人澡人人妻人| 免费在线观看完整版高清| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| 国产视频首页在线观看| 人妻一区二区av| 国产精品欧美亚洲77777| 精品一区二区三卡| 日本猛色少妇xxxxx猛交久久| 国产一区有黄有色的免费视频| 18禁裸乳无遮挡动漫免费视频| 国产野战对白在线观看| 亚洲免费av在线视频| 新久久久久国产一级毛片| 99热网站在线观看| 在线观看人妻少妇| 午夜激情久久久久久久| 成年动漫av网址| 啦啦啦在线观看免费高清www| 亚洲成人免费av在线播放| 最新在线观看一区二区三区 | 赤兔流量卡办理| av网站在线播放免费| 亚洲成人免费av在线播放| 亚洲一级一片aⅴ在线观看| 亚洲成人免费av在线播放| 一本—道久久a久久精品蜜桃钙片| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 欧美97在线视频| 777米奇影视久久| 国产色婷婷99| www.熟女人妻精品国产| 亚洲,欧美精品.| 国产精品一区二区在线不卡| 国产探花极品一区二区| 天美传媒精品一区二区| 成人免费观看视频高清| 亚洲综合精品二区| 亚洲欧美一区二区三区黑人| 丝袜脚勾引网站| 男女边摸边吃奶| 国产精品国产三级国产专区5o| 免费人妻精品一区二区三区视频| 欧美 日韩 精品 国产| 天天躁夜夜躁狠狠久久av| 国产免费一区二区三区四区乱码| 精品国产一区二区三区四区第35| 欧美精品一区二区大全| 亚洲精品一二三| 菩萨蛮人人尽说江南好唐韦庄| 777久久人妻少妇嫩草av网站| 精品视频人人做人人爽| 亚洲av福利一区| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| www日本在线高清视频| av有码第一页| 久久97久久精品| 不卡视频在线观看欧美| 黄片播放在线免费| 午夜免费观看性视频| 久久久精品国产亚洲av高清涩受| 99热全是精品| 亚洲国产精品一区二区三区在线| 91老司机精品| 在线亚洲精品国产二区图片欧美| 美国免费a级毛片| 日韩 欧美 亚洲 中文字幕| 街头女战士在线观看网站| 免费日韩欧美在线观看| 午夜免费鲁丝| 国产亚洲精品第一综合不卡| 亚洲精品久久久久久婷婷小说| 免费日韩欧美在线观看| av有码第一页| 亚洲国产精品一区三区| 男女无遮挡免费网站观看| 日韩制服骚丝袜av| 日韩 亚洲 欧美在线| 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 美女午夜性视频免费| 国产成人91sexporn| 国产精品99久久99久久久不卡 | 欧美少妇被猛烈插入视频| 国产精品香港三级国产av潘金莲 | 久久综合国产亚洲精品| 考比视频在线观看| 男女高潮啪啪啪动态图| 欧美变态另类bdsm刘玥| www.精华液| www.自偷自拍.com| 婷婷色综合大香蕉| 老汉色∧v一级毛片| 香蕉国产在线看| 黄网站色视频无遮挡免费观看| 最新的欧美精品一区二区| 欧美人与性动交α欧美软件| 免费看不卡的av| 国产激情久久老熟女| 搡老岳熟女国产| 777米奇影视久久| 精品国产一区二区三区四区第35| 精品一区二区三区四区五区乱码 | 精品第一国产精品| 韩国高清视频一区二区三区| 国产av精品麻豆| 天天躁日日躁夜夜躁夜夜| 啦啦啦 在线观看视频| 精品一区在线观看国产| 女人高潮潮喷娇喘18禁视频| 久久精品久久久久久噜噜老黄| 超碰成人久久| av有码第一页| 中文字幕色久视频| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜制服| 午夜91福利影院| 一本一本久久a久久精品综合妖精| 成人国产av品久久久| 亚洲av综合色区一区| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 午夜福利在线免费观看网站| 一级片'在线观看视频| 亚洲av成人不卡在线观看播放网 | 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 精品国产超薄肉色丝袜足j| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的| 欧美精品一区二区大全| 热99国产精品久久久久久7| 国产成人精品久久久久久| 美国免费a级毛片| 久久国产精品大桥未久av| 又黄又粗又硬又大视频| 香蕉丝袜av| 老司机深夜福利视频在线观看 | 不卡av一区二区三区| 国产毛片在线视频| 亚洲情色 制服丝袜| 两个人免费观看高清视频| 亚洲免费av在线视频| 精品电影一区二区在线| 日韩高清综合在线| 少妇 在线观看| 色播亚洲综合网| 女同久久另类99精品国产91| 两个人视频免费观看高清| 亚洲色图综合在线观看| 国产私拍福利视频在线观看| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 亚洲精华国产精华精| 美女国产高潮福利片在线看| 国产成人精品久久二区二区91| 国产乱人伦免费视频| 国产高清有码在线观看视频 | 国产精品爽爽va在线观看网站 | 国产成人啪精品午夜网站| 国产又爽黄色视频| 欧美中文综合在线视频| 国产91精品成人一区二区三区| 亚洲精品国产色婷婷电影| 国产精品一区二区三区四区久久 | 亚洲成国产人片在线观看| 国产xxxxx性猛交| 成人18禁在线播放| 久久精品国产亚洲av香蕉五月| 人妻久久中文字幕网| 两个人免费观看高清视频| 高清在线国产一区| 免费在线观看完整版高清| 日日摸夜夜添夜夜添小说| 亚洲人成77777在线视频| 亚洲国产欧美网| 99国产综合亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 欧美在线一区亚洲| 国产精品99久久99久久久不卡| 嫁个100分男人电影在线观看| 久久婷婷成人综合色麻豆| 久久香蕉激情| 91字幕亚洲| 女人被狂操c到高潮| av在线播放免费不卡| 国产av一区在线观看免费| 国产精品国产高清国产av| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 99久久精品国产亚洲精品| 男人舔女人下体高潮全视频| 亚洲成人久久性| 国产av在哪里看| 99久久精品国产亚洲精品| av天堂在线播放| 国产精品亚洲av一区麻豆| 欧美日本亚洲视频在线播放| 欧美中文日本在线观看视频| 亚洲成人精品中文字幕电影| 精品久久久久久成人av| 久久亚洲精品不卡| 免费人成视频x8x8入口观看| 午夜福利影视在线免费观看| 久久久久久国产a免费观看| 色播在线永久视频| 香蕉丝袜av| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 老司机深夜福利视频在线观看| 搡老妇女老女人老熟妇| 久久 成人 亚洲| 亚洲狠狠婷婷综合久久图片| 一级毛片高清免费大全| 精品久久久久久久毛片微露脸| 中文字幕精品免费在线观看视频| 日日干狠狠操夜夜爽| 国产精品国产高清国产av| 国产精品香港三级国产av潘金莲| 色精品久久人妻99蜜桃| 午夜激情av网站| 乱人伦中国视频| 国产亚洲av嫩草精品影院| 女同久久另类99精品国产91| 日韩大码丰满熟妇| 午夜两性在线视频| 在线观看免费日韩欧美大片| 精品一区二区三区视频在线观看免费| 久久午夜亚洲精品久久| 久久久久久免费高清国产稀缺| 麻豆av在线久日| 国产精品乱码一区二三区的特点 | 亚洲无线在线观看| 国产亚洲精品综合一区在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久av网站| 亚洲男人的天堂狠狠| 国产真人三级小视频在线观看| av片东京热男人的天堂| 精品人妻1区二区| 精品欧美一区二区三区在线| 成人三级黄色视频| 精品久久蜜臀av无| 在线观看一区二区三区| xxx96com| 国产三级黄色录像| 搡老熟女国产l中国老女人| 啦啦啦 在线观看视频| 97超级碰碰碰精品色视频在线观看| 国产三级在线视频| 日本vs欧美在线观看视频| 真人一进一出gif抽搐免费| www.精华液| 丰满的人妻完整版| 亚洲情色 制服丝袜| 一级a爱片免费观看的视频|