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

    東北地區(qū)森林凋落葉分解速率與氣候、林型、林分光照的關(guān)系

    2016-07-21 01:25:31王襄平張新平周海城方精云
    生態(tài)學(xué)報(bào) 2016年8期
    關(guān)鍵詞:林型氣候

    吳 鵬, 王襄平,*, 張新平, 朱 彪, 周海城, 方精云

    1 北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083 2 北京大學(xué)生態(tài)學(xué)系, 北京大學(xué)地表過程與分析模擬教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100871 3 吉林省長白山自然保護(hù)管理中心, 安圖 133613

    ?

    東北地區(qū)森林凋落葉分解速率與氣候、林型、林分光照的關(guān)系

    吳鵬1, 王襄平1,*, 張新平2, 朱彪2, 周海城3, 方精云2

    1 北京林業(yè)大學(xué)森林資源與生態(tài)系統(tǒng)過程北京市重點(diǎn)實(shí)驗(yàn)室, 北京100083 2 北京大學(xué)生態(tài)學(xué)系, 北京大學(xué)地表過程與分析模擬教育部重點(diǎn)實(shí)驗(yàn)室, 北京100871 3 吉林省長白山自然保護(hù)管理中心, 安圖133613

    摘要:在東北長白山、張廣才嶺、小興安嶺、大興安嶺的主要森林類型中設(shè)置26塊樣地,進(jìn)行為期3a(2004—2006年)凋落葉分解實(shí)驗(yàn),以研究氣候、林型、林冠透光率對凋落葉分解速率的相對影響大小。結(jié)果表明,不同林型凋落葉分解速率依次為:落葉闊葉林>針闊葉混交林>落葉針葉林>常綠針葉林>岳樺林。 對分解速率影響因素的分析表明,氣候因子(熱量和水分)對分解速率有較強(qiáng)的解釋力,分別解釋了分解常數(shù)k和分解95%所需時(shí)間(t95%)的55.5%和65.0%的變異。但是,氣候?qū)Ψ纸馑俾实挠绊懺诤艽蟪潭壬鲜峭ㄟ^與林型、林冠透光率的協(xié)同作用而實(shí)現(xiàn)的,其獨(dú)立解釋力并不大(<9%)。氣候的變化導(dǎo)致林型(物種組成)的變化、進(jìn)而影響分解速率,這一因素解釋了分解參數(shù)變異的46.8%(k) 和56.8%(t95%)。與此同時(shí),氣候和林型的變化還導(dǎo)致林冠透光率的變化,隨著熱量水平的上升林冠透光率下降、間接提高分解速率。這一因素分別解釋了k值和t95%變異的23.9%和22.3%。研究結(jié)果表明,氣候?qū)Φ蚵淙~分解的影響主要是通過對物種組成、林冠結(jié)構(gòu)(影響透光率)等生物因素的間接作用實(shí)現(xiàn)的。忽視這些生物因素、簡單研究氣候和分解速率的關(guān)系可能難以正確預(yù)測未來氣候變化對凋落物分解的影響。

    關(guān)鍵詞:凋落葉分解速率; 氣候; 林型; 林冠透光率; 中國東北地區(qū)

    凋落物分解是森林碳循環(huán)和全球碳平衡的一個(gè)重要環(huán)節(jié)[1- 3]。由凋落物分解返回大氣中的碳含量是全球碳平衡的重要組成部分[4]。據(jù)估計(jì),凋落物(包括根凋落物)分解產(chǎn)生的碳大約占到全球陸地生態(tài)系統(tǒng)年際碳流通總量的70%,達(dá)到68 Pg C/a(Pg=1015g)[5]。

    我國東北地區(qū)森林覆蓋率高,碳儲(chǔ)量大,僅黑龍江、吉林、遼寧三省的森林碳儲(chǔ)量就占全國森林總碳儲(chǔ)量的24%—31%[6],是東亞地區(qū)的一個(gè)重要碳匯[7]。而高緯度、寒冷地區(qū)凋落物分解由于受熱量制約,對氣候變暖更為敏感[3, 8]。據(jù)估計(jì),到21世紀(jì)后期,中國東北地區(qū)氣溫將可能較目前變暖3.0℃或以上,降水將可能增加[9]。在氣候變暖的背景下,研究我國東北地區(qū)凋落物分解與氣候的關(guān)系具有重要意義。

    研究表明,凋落物分解速率受多種因素的作用,包括氣候和局域環(huán)境條件,凋落物質(zhì)量,植被類型及土壤生物等[10- 13]。關(guān)于東北森林凋落物分解已經(jīng)開展了不少研究,如王鳳友等對涼水地區(qū)紅松和幾種闊葉樹凋落葉的分解和養(yǎng)分釋放的研究[14],沈海龍對帽兒山樟子松人工林針闊葉凋落物分解動(dòng)態(tài)的研究[15],張東來對帽兒山6種主要樹種凋落物分解速率的研究[16],以及郭忠玲等對長白山主要樹種和林型凋落物分解的比較研究等[17]。這些研究雖然揭示了某一地點(diǎn)不同林型或樹種的分解速率差異,但地理跨度小,未能探討大尺度上凋落物分解的影響因素。在整個(gè)東北地區(qū)的尺度上,宏觀氣候因素、森林類型、局域環(huán)境條件如何共同作用影響凋落葉的分解?這些因素的相對重要性如何?相關(guān)問題還需要定量研究,以深入理解森林碳循環(huán)過程,并可為生物地球化學(xué)模型提供重要參數(shù)。

    大尺度的研究表明氣候和植被類型對凋落物分解有重要影響[13, 18- 19]。凋落物分解速率在植被或森林類型之間的差異,一方面是因?yàn)椴煌闹脖环植荚诓煌臍夂驐l件下,氣候的差異對分解速率有著直接的影響[20];另一方面是由于氣候通過其它因素的間接作用。比如,不同地區(qū)的植被物種組成不同、因而植物功能性狀不同[21- 23],后者直接影響凋落物質(zhì)量,是凋落物分解的重要影響因素[13, 24]。前述東北地區(qū)的研究也表明,相同氣候下的不同林型由于樹種組成不同,分解速率有顯著差異。然而,這兩種因素對凋落物分解的相對影響大小,卻很少有研究涉及。氣候變化顯然會(huì)導(dǎo)致物種組成的變化、甚至森林類型的變化[25- 26],如果后一種機(jī)制對凋落物分解的影響占主導(dǎo)作用,則意味著不考慮物種組成變化的因素,是無法準(zhǔn)確預(yù)測氣候變化對森林碳循環(huán)的影響的。

    此外,小尺度的研究常發(fā)現(xiàn)林內(nèi)光照對凋落物分解的重要影響。如林冠空隙可通過改變林內(nèi)微環(huán)境影響土壤生物活動(dòng),進(jìn)而加速或阻礙凋落物分解[27- 28],對于針葉林其作用更是不可低估[29]。在人為干擾不大的情況下,林內(nèi)光照是和氣候條件導(dǎo)致的林型、物種組成差異有關(guān)的。如大興安嶺的明亮針葉林的林內(nèi)光照條件要好于東北東部山地的針闊混交林,有研究表明這種差異對林下植物多樣性是有明顯影響的[30]。因此光照條件也有可能和氣候、森林類型等因素協(xié)同作用于凋落物的分解。但迄今為止,大尺度上的凋落物分解研究尚未考慮過這一因素。

    本文對東北主要山地(大小興安嶺、張廣才嶺、長白山)主要森林類型的凋落葉分解進(jìn)行了系統(tǒng)研究。目的在于:(1)研究氣候?qū)|北地區(qū)主要森林類型的凋落葉分解速率的影響;(2)探討森林類型、林冠透光率如何與氣候因子共同作用,影響凋落葉的分解速率。

    1研究地點(diǎn)及方法

    1.1研究地點(diǎn)及樣地設(shè)置

    研究選取東北林區(qū)的4個(gè)地點(diǎn),分別是長白山(14塊樣地),張廣才嶺的帽兒山(6塊樣地),小興安嶺的涼水(4塊樣地)和大興安嶺的根河(2塊樣地)共設(shè)置26塊固定觀測樣地,樣地面積為1 000 m2(20 m×50 m)。各研究地點(diǎn)的概況見張新平等[31]。長白山的樣地沿海拔梯度設(shè)置(海拔從650m至1940m),包含了各森林垂直帶。其余海拔梯度不大的地點(diǎn)則選擇當(dāng)?shù)氐湫偷纳诸愋驮O(shè)立樣地。樣地設(shè)置時(shí)選擇沒有明顯干擾痕跡的郁閉林分。本研究的樣地包括了東北地區(qū)5種主要森林類型,在緯度上跨越了溫帶和寒溫帶,在海拔梯度上包括了從低山帶到亞高山帶的各種群落類型,具有較好的代表性。林型包括落葉闊葉林、針闊葉混交林、落葉針葉林、常綠針葉林、岳樺林。岳樺林本屬于落葉闊葉林,但由于其生長在高山林線附近,與低海拔的落葉闊葉林不論在氣候條件、物種組成、群落結(jié)構(gòu)等方面均有很大的不同(表1),因此在分析時(shí)將其單獨(dú)作為一種林型予以考慮。

    表1 各林型氣候及主要樹種組成概況

    DB: Deciduous broadleaf forest; DNM: Deciduous broadleaf and needle-leaf mixed forest; DN: Deciduous needle-leaf forest; EN: Evergreen needle-leaf forest; BE:Betulaermaniiforest, WI, warmth index; MAP, mean annual precipitation.

    1.2凋落葉分解測定

    用網(wǎng)袋法進(jìn)行分解實(shí)驗(yàn),網(wǎng)袋大小為15 cm×20 cm,網(wǎng)孔為1 mm×1 mm[32]。于2003年9月底在4個(gè)地點(diǎn)的每塊樣地內(nèi)收集新鮮凋落葉、自然風(fēng)干。每個(gè)樣地取風(fēng)干葉樣300 g于65 ℃烘至恒重,測定含水率。每個(gè)分解袋內(nèi)裝入30 g風(fēng)干葉樣品放回原樣地的土壤表層。每個(gè)樣地放置60袋,于每年4—11月的月初取2袋。研究持續(xù)了3a,到2006年11月為止。分解袋取回后去除雜質(zhì)、土壤,在65℃下烘至恒重后稱重。

    凋落葉分解過程采用修正的[33]指數(shù)衰減函數(shù)進(jìn)行擬合[34- 35]:

    Xt/X0=a×e-kt

    (1)

    式中,Xt為分解時(shí)間X0時(shí)的凋落葉干重,X0為凋落葉初始干重,Xt/X0即為凋落葉剩余比例,k是分解常數(shù),以年為時(shí)間單位,a為修正系數(shù)。通過擬合出的函數(shù)計(jì)算得出凋落葉分解95%所需時(shí)間(t95%),以更全面地反映凋落葉分解速率[33, 36]。

    1.3林冠空隙度的測定

    林冠空隙會(huì)影響到光照、熱量、水分等在林冠層下的再分配,是影響林內(nèi)微環(huán)境的重要因素,因此進(jìn)行了林冠透光狀況的測定,以研究它對凋落物分解速率的可能影響。在每塊樣地內(nèi)均勻布設(shè)5個(gè)點(diǎn),通過WinScanopy林冠分析系統(tǒng)拍攝數(shù)碼照片。拍攝時(shí)間選在天氣晴朗的8:00—10:00和14:00—16:00,以消除太陽直射光斑的影響。對拍攝的180°角的林冠層數(shù)碼照片,應(yīng)用WinScanopy For Canopy Analysis軟件進(jìn)行分析(加拿大Regent Instruments公司),得到林冠空隙度、林冠開度等光照參數(shù),以反映林冠層的透光率。本研究使用的指標(biāo)為林冠空隙度(取5個(gè)點(diǎn)的平均值),林冠開度與其生物學(xué)意義、解釋力接近,因而在最終數(shù)據(jù)分析中沒有使用。

    1.4氣候指標(biāo)的計(jì)算

    采用如下方法推算各樣地的氣溫、降水等氣候指標(biāo):

    T (P) =a + bLatitude+ cLongitude+ dAltitude

    (2)

    式中,T(或P) 為1—12 月的月平均溫度(月平均降水),a,b,c,d為常數(shù)[37]。采用Wang等利用國家氣象臺(tái)站數(shù)據(jù)建立的模型估算各樣地的氣候指標(biāo)[38]。對這些模型的驗(yàn)證結(jié)果表明,年均溫預(yù)測的精度為96%,年降水的預(yù)測精度為79%,足以滿足本研究的需要[38]。

    利用上述方法得到所有樣地的1—12月月平均溫度和降水,并計(jì)算如下氣候指標(biāo):

    1)溫暖指數(shù)(WI)WI=∑(ti-5),其中ti為5℃以上的月均溫。

    2)年降水量(MAP)一年各月降水量之和。

    1.5數(shù)據(jù)分析

    采用回歸分析和方差分析探討三類因素對分解速率的影響。1)氣候因子:包括溫暖指數(shù)、年降水量,分別反映熱量和水分的有效性[38];2)林分光照狀況,以林冠空隙度反映林冠層透光率;3)林型,即落葉闊葉林、針闊葉混交林、落葉針葉林、常綠針葉林、岳樺林等5種林型。數(shù)據(jù)分析之前均先進(jìn)行正態(tài)性檢驗(yàn)和方差齊性檢驗(yàn)。各林型之間凋落葉分解參數(shù)的比較采用多重比較Duncan檢驗(yàn)。

    為了分析氣候、林冠透光率和林型之間的協(xié)同作用和獨(dú)立作用,采用變異分離的方法,將總變異分解為如下組分:(1)3種因素各自的獨(dú)立作用(a, b, c);(2)兩因素之間的協(xié)同作用(d, e, f);(3)3種因素之間的協(xié)同作用(g);(4)未解釋的變異。各組分的計(jì)算方法見Heikkinen等[39]。采用Wang等提出的隨機(jī)化檢驗(yàn)方法進(jìn)行各組分顯著性的檢驗(yàn)[40]。

    數(shù)據(jù)分析使用R 3.1.1(R Development Core Team 2014)軟件進(jìn)行。

    2結(jié)果與討論

    2.1各林型的凋落葉分解速率

    比較不同林型凋落葉分解參數(shù)可以看出(圖1),各林型凋落葉分解常數(shù)k的平均值依次為:落葉闊葉林(0.338)、針闊葉混交林(0.293)、落葉針葉林(0.207)、常綠針葉林(0.146)、岳樺林(0.097),并且落葉闊葉林、針闊葉混交林與常綠針葉林、岳樺林具有顯著性差異。這與長白山其他研究中[17, 41]得出的相同(似)林型的k值相比偏小,可能主要是由于這些研究的時(shí)間跨度一般為1—2a,而本研究為3a。凋落葉分解速率前期較大,后期則逐漸遞減,因此研究時(shí)間較短時(shí)得到的k值會(huì)偏大。不少研究都指出過長時(shí)間凋落物分解實(shí)驗(yàn)的重要性[42- 43]。

    對于分解95%所需時(shí)間,落葉闊葉林、針闊葉混交林的平均值(分別為9.2、9.9a)均顯著小于常綠針葉林、岳樺林(分別為19.7、24.6a),落葉針葉林介于中間(14.5a)。

    圖1 不同林型凋落葉分解速率的差異Fig.1 Comparison of litter decomposition rate between forest types decomposition coefficient (k)上標(biāo)字母有相同者在P<0.05水平上無顯著差異;BE、DB、DN、DNM、EN,森林林型,分別為岳樺林(Betula ermanii forest,)、落葉闊葉林(Deciduous broadleaf forest)、落葉針葉林(Deciduous needle-leaf forest)、針闊葉混交林(Deciduous broadleaf and needle-leaf mixed forest)、常綠針葉林(Evergreen needle-leaf forest)

    2.2凋落葉分解與各影響因素的關(guān)系

    2.2.1分解參數(shù)與海拔、緯度的關(guān)系

    為考察分解速率與海拔的關(guān)系,選取長白山的樣地進(jìn)行分析(表2, 圖2)。結(jié)果表明,分解常數(shù)k、分解95%所需時(shí)間與海拔之間的回歸關(guān)系均極為顯著(P<0.001)。隨海拔的升高,k值下降,分解95%所需時(shí)間延長。由于各研究地點(diǎn)樣地的海拔差異較大(長白山樣地在650—1940 m之間,其余4個(gè)地點(diǎn)的樣地在300—840 m之間),沒有單獨(dú)分析分解參數(shù)和緯度的相關(guān)性。不過,對分解參數(shù)與緯度、海拔的多元回歸分析表明,在回歸模型中緯度的作用并不顯著(k:P=0.92;t95%:P=0.83),而海拔則在P=0.001水平顯著。這一關(guān)于緯度的統(tǒng)計(jì)結(jié)果可能與本研究的不同地點(diǎn)海拔差異較大、緯度梯度還不夠大有關(guān),在全國和全球大尺度上的研究都發(fā)現(xiàn),分解速率隨緯度上升顯著下降[13, 19]。

    表2長白山各林型樣地海拔分布

    Table 2Altitudinal ranges for plots of each forest type on Mt. Changbai

    林型Foresttype海拔Altitude/m樣地?cái)?shù)No.plots落葉闊葉林DB9501針闊葉混交林DNM650—9703落葉針葉林DN1420—14402常綠針葉林EN1010—16605岳樺林BE1850—19403

    2.2.2分解參數(shù)與氣候因子的關(guān)系

    海拔、緯度等地理因子與分解速率的關(guān)系,實(shí)質(zhì)是通過對熱量、水分的分配而起作用。對氣候因子WI、MAP與分解參數(shù)關(guān)系的分析表明(圖3),各相關(guān)關(guān)系均極為顯著,分解常數(shù)k與熱量(WI)呈正相關(guān),與水分(MAP)呈負(fù)相關(guān);凋落葉分解95%所需時(shí)間則反之。東北森林區(qū)降水較充沛,熱量是生產(chǎn)力的主要限制因子,而水分不構(gòu)成限制[30-31]。因此分解速率與降水的負(fù)相關(guān)關(guān)系可能主要是由于MAP和WI的負(fù)相關(guān)關(guān)系所導(dǎo)致(R2=0.51;P<0.001)。一般認(rèn)為,在熱帶生態(tài)系統(tǒng),降水量增加對凋落物分解有明顯的促進(jìn)作用,但在一些溫帶生態(tài)系統(tǒng)中,較高的降水量反而可能使凋落物分解減慢[44- 45]。一方面,過多的雨水導(dǎo)致土壤嫌氣條件,不利于土壤生物的活動(dòng);另一方面,降水可能改變凋落物基質(zhì)質(zhì)量,比如降低磷濃度、增加木質(zhì)素濃度等[46],從而對分解產(chǎn)生負(fù)效應(yīng)。本文中分解速率與降水的負(fù)相關(guān)也可能與這種因素有關(guān),但這種因素的貢獻(xiàn)有多大還需進(jìn)一步研究。

    2.2.3分解參數(shù)與林分光照的關(guān)系

    對林冠空隙度與分解參數(shù)關(guān)系的分析表明(圖4),林冠空隙度與分解常數(shù)k呈顯著負(fù)相關(guān),與分解95%所需時(shí)間相反。即林冠空隙越大,凋落葉分解速率減小。Zhang和Zak的研究發(fā)現(xiàn)林隙較小處的凋落物比林隙大的分解得更快,并且認(rèn)為較大的林冠空隙會(huì)通過改變微環(huán)境條件而顯著減小微生物活性和分解速率[47]。一些其他研究也得出了類似的結(jié)果[28, 48]。圖4與這些研究結(jié)果是一致的。

    圖2 分解常數(shù)k、分解95%所需時(shí)間(t95%)隨海拔的變化Fig.2 Changes of leaf litter decomposition coefficient (k) and the time for 95% litter decomposed (t95%) with altitude

    圖3 分解常數(shù)k、分解95%所需時(shí)間(t95%)與氣候因子溫暖指數(shù)(WI)和年降水量(MAP)的關(guān)系Fig.3 Relationship of decomposition coefficient k and the time for 95% litter decomposed (t95%) with climate indices WI, warmth index; MAP, Mean annual precipitation

    圖4 分解常數(shù)k、分解95%所需時(shí)間(t95%)與林冠空隙度的關(guān)系Fig.4 The relationship between decomposition coefficient (k) and the time for 95% litter decomposed (t95%) with gap fraction

    2.3氣候、林型、林冠透光率對凋落葉分解的相對影響

    前述結(jié)果表明,氣候、林冠透光率、林型等因素對凋落葉分解都有顯著的影響。為了研究這些因素如何共同作用于分解速率,采用變異分離來反映各因素的協(xié)同作用和獨(dú)立作用。圖5表明,氣候因子(WI和MAP)共計(jì)解釋了k值的55.5%的變異(a+d+e+g),而林型的總解釋力為66.0% (b+d+f+g),林冠空隙度的總解釋力為24.2% (c+e+f+g)。分解95%所需時(shí)間的結(jié)果與k值相似。這一結(jié)果與以往的研究一致,即氣候和林型是凋落葉分解的重要影響因素[13, 19]。

    不過,圖5進(jìn)一步揭示,氣候、林型、林冠透光率對分解速率的解釋力主要屬于三種因素之間的協(xié)同作用,各因素的獨(dú)立解釋力(a、b、c)并不大。隨機(jī)化檢驗(yàn)的結(jié)果表明,在各方差組分中,僅d和g的作用是顯著的。其中氣候和林型的協(xié)同解釋力(d+g)高達(dá)46.8% (k) —56.8% (t95%)。這說明氣候?qū)Φ蚵淙~分解的作用,在很大程度上是由于氣候變化導(dǎo)致林型(物種組成)的變化、進(jìn)而影響分解速率(見引言)。至少在本研究中,氣候?qū)Ψ纸馑俾实闹苯佑绊懖⒉淮?<9%)。

    圖5 氣候、林型、林冠空隙度對分解參數(shù)影響的變異分解圖Fig.5 Variance partitioning for the effects of climate, forest type and gap fraction on decomposition constant (k) and the time for 95% litter decomposed(t95%), in terms of the percentage of sum of squares explained圖中給出了各組分占總變異的百分比;a,b,c表示各因素之間的獨(dú)立解釋力;d,e,f,g表示各因素之間的協(xié)同解釋力

    雖然在很多小尺度研究中已經(jīng)證實(shí)林內(nèi)光照、林冠空隙對凋落物分解有重要的影響[28, 48],但由于大尺度的凋落物分解研究沒有測定過林分光照狀況,宏觀氣候和林內(nèi)光照對分解速率的相對作用大小還沒有進(jìn)行過研究。事實(shí)上,一些研究早就指出過,氣候條件的變化會(huì)顯著影響林內(nèi)的光環(huán)境[49]:在高緯度地區(qū)的針葉林,樹冠呈錐形、樹木間間距較大、林冠垂直層次簡單;而在南方的闊葉林中,樹冠趨向扁平、林冠垂直層次復(fù)雜且林分密度較大。這種氣候條件導(dǎo)致的林冠結(jié)構(gòu)差異,會(huì)顯著影響林內(nèi)光照,并對森林更新、多樣性產(chǎn)生重要影響[30, 49]。因此,在人為干擾不大的情況下(如本文的郁閉林分樣地),林內(nèi)光照條件應(yīng)當(dāng)是隨氣候條件而變化的。本研究的測定結(jié)果也證實(shí)了這一點(diǎn):隨著熱量(WI)的增加,林冠透光率顯著下降(圖6)。在東北地區(qū),從北到南隨著熱量的增加,森林類型由垂直結(jié)構(gòu)簡單、透光率高的興安落葉松林逐漸過渡到垂直層次復(fù)雜的針闊混交林[30, 50]。圖6的結(jié)果與這一事實(shí)是一致的。

    圖6 林冠透光率和溫暖指數(shù)的關(guān)系 Fig.6 The relationship between gap fraction and warmth index (WI)

    變異分離的結(jié)果表明(圖5), 氣候、林型和林冠透光率的協(xié)同作用(g)是最重要的變異組分之一(而且作用顯著),分別解釋了k和t95%變異的23.9%和22.3%。因此,氣候通過林分光照等生物因素對凋落物分解的間接作用是不可忽視的,即熱量的增加除了對分解速率的直接影響外,還可通過森林類型(物種組成和林冠結(jié)構(gòu))的改變導(dǎo)致的林冠透光率下降,從而間接地促進(jìn)凋落葉的分解(圖4)[48, 51]。

    3結(jié)論與討論

    本研究結(jié)果表明,各林型凋落物分解速率大小依次為:落葉闊葉林、針闊葉混交林、落葉針葉林、常綠針葉林、岳樺林。從總體上看,闊葉>針葉,落葉>常綠,這和大多數(shù)研究得出的結(jié)論相同[52]。一方面,相比于針葉凋落物,闊葉凋落物具有更高的養(yǎng)分含量和更低的木質(zhì)素、多酚含量,因此可能分解得更快[53]。另一方面,闊葉林或混交林凋落物形成的枯枝落葉層相比于針葉林的更有利于凋落物分解。因?yàn)橐话阏J(rèn)為闊葉形成的細(xì)腐殖質(zhì)枯枝落葉層(mull forest floor)所含營養(yǎng)元素更為豐富,大型土壤動(dòng)物物種豐富度高,可加速凋落物的機(jī)械粉碎,刺激微生物活性,從而加速分解進(jìn)程[54]。

    凋落物分解是多種因素共同作用的結(jié)果,并且不同因素在凋落物分解的不同階段所起的相對作用大小是不同的。因此,弄清楚這些因素是如何相互影響并且共同作用于凋落物分解過程就尤其重要。

    氣候、林型、林冠透光率對分解速率的影響結(jié)果表明,三者之間的協(xié)同作用對于凋落物分解起著控制性作用。氣候中溫度、水分在不同地區(qū)對于凋落物分解起著明顯不同的作用。在干旱或半干旱地區(qū),隨海拔升高,分解速率反而加快,這是由于在這些地區(qū)海拔較高的地方,盡管溫度降低,但濕度也會(huì)增大,總的作用還是加快分解速率[55]。而對于水分不構(gòu)成限制性因素的東北地區(qū),分解速率則隨海拔的升高而減小。本研究結(jié)果表明,氣候?qū)Ψ纸馑俾实挠绊懺诤艽蟪潭壬鲜峭ㄟ^如下機(jī)制實(shí)現(xiàn)的:1) 氣候的變化導(dǎo)致林型的變化,即通過影響物種組成導(dǎo)致凋落葉質(zhì)量的變化[24],進(jìn)而影響分解速率;2) 氣候和林型的變化還導(dǎo)致林冠透光率的變化,間接作用于分解速率。林型對分解速率的解釋力為66.0%,可見不同林型間凋落物質(zhì)量的差異對分解速率具有較強(qiáng)的解釋力,這可以通過不同樹種所形成的獨(dú)特周圍環(huán)境條件的影響來起作用,并且植物物種特性和植物-土壤的長期互相作用深刻地影響著凋落物的分解[56]。林冠透光率相比于氣候、林型的解釋力較小,但由于其受氣候、林型雙重影響,對于凋落物的分解作用依然不可忽視,不少研究發(fā)現(xiàn)林冠空隙大小對于凋落物分解具有重要作用[28,51]。因此,需要深入研究氣候和生物因素(如物種組成、林冠結(jié)構(gòu))如何交互作用影響凋落物分解,才能更好地理解氣候變化對森林碳循環(huán)的影響。

    參考文獻(xiàn)(References):

    [1]Aerts R. Climate, leaf litter chemistry and leaf litter decomposition in terrestrial ecosystems: a triangular relationship. Oikos, 1997, 79(3): 439- 449.

    [2]Couteaux M M, Bottner P, Berg B. Litter decomposition, climate and liter quality. Trends in Ecology & Evolution, 1995, 10(2): 63- 66.

    [3]Robinson C H. Controls on decomposition and soil nitrogen availability at high latitudes. Plant and Soil, 2002, 242(1): 65- 81.

    [4]Shaver G R, Billings W D, Chapin F S, Giblin A E, Nadelhoffer K J, Oechel W C, Rastetter E B. Global change and the carbon balance of arctic ecosystems. BioScience, 1992, 42(6): 433- 441.

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

    [6]Fang J Y, Chen A P, Peng C H, Zhao S Q, Ci L J. Changes in forest biomass carbon storage in China between 1949 and 1998. Science, 2001, 292(5525): 2320- 2322.

    [7]Tan K, Piao S L, Peng C H, Fang J Y. Satellite- based estimation of biomass carbon stocks for northeast China′s forests between 1982 and 1999. Forest Ecology and Management, 2007, 240(1/3): 114- 121.

    [8]Hobbie S E, Nadelhoffer K J, H?gberg P. A synthesis: the role of nutrients as constraints on carbon balances in boreal and arctic regions. Plant and Soil, 2002, 242(1): 163- 170.

    [9]趙宗慈, 羅勇. 21世紀(jì)中國東北地區(qū)氣候變化預(yù)估. 氣象與環(huán)境學(xué)報(bào), 2007, 23(3): 1- 4.

    [10]Aber J D, Melillo J M, Nadelhoffer K J, Pastor J, Boone R D. Factors controlling nitrogen cycling and nitrogen saturation in northern temperate forest ecosystems. Ecological Applications, 1991, 1(3): 303- 315.

    [11]Wardle D, Lavelle P. Linkages between soil biota, plant litter quality and decomposition // Cadisch G, Giller K E, eds. Driven by Nature: Plant Litter Quality and Decomposition. Wallingford: CAB International, 1997: 107- 124.

    [12]Austin A T, Vivanco L. Plant litter decomposition in a semi- arid ecosystem controlled by photodegradation. Nature, 2006, 442(7102): 555- 558.

    [13]Zhang D Q, Hui D F, Luo Y Q, Zhou G Y. Rates of litter decomposition in terrestrial ecosystems: global patterns and controlling factors. Journal of Plant Ecology, 2008, 1(2): 85- 93.

    [14]王鳳友, 王業(yè)蘧. 紅松和五種闊葉樹枯葉的分解及養(yǎng)分釋放 // 周曉峰. 森林生態(tài)系統(tǒng)定位研究.第一集. 哈爾濱: 東北林業(yè)大學(xué)出版社, 1991: 250- 259.

    [15]沈海龍, 丁寶永, 沈國舫, 陳愛民. 樟子松人工林下針闊葉凋落物分解動(dòng)態(tài). 林業(yè)科學(xué), 1996, 32(5): 393- 402.

    [16]張東來, 毛子軍, 朱勝英, 周彪. 黑龍江省帽兒山林區(qū)6種主要林分類型凋落物研究. 植物研究, 2008, 28(1): 104- 108.

    [17]郭忠玲, 鄭金萍, 馬元丹, 李慶康, 于貴瑞, 韓士杰, 范春楠, 劉萬德. 長白山各植被帶主要樹種凋落物分解速率及模型模擬的試驗(yàn)研究. 生態(tài)學(xué)報(bào), 2006, 26(4): 1037- 1046.

    [18]王健健, 王永吉, 來利明, 趙學(xué)春, 王飛, 申國珍, 賴江山, 魯洪斌, 趙春強(qiáng), 鄭元潤. 我國中東部不同氣候帶成熟林凋落物生產(chǎn)和分解及其與環(huán)境因子的關(guān)系. 生態(tài)學(xué)報(bào), 2013, 33(15): 4818- 4825.

    [19]唐仕姍, 楊萬勤, 殷睿, 熊莉, 王海鵬, 王濱, 張艷, 彭艷君, 陳青松, 徐振鋒. 中國森林生態(tài)系統(tǒng)凋落葉分解速率的分布特征及其控制因子. 植物生態(tài)學(xué)報(bào), 2014, 38(6): 529- 539

    [20]Liski J, Nissinen A, Erhard M, Taskinen O. Climatic effects on litter decomposition from arctic tundra to tropical rainforest. Global Change Biology, 2003, 9(4): 575- 584.

    [21]Liu C, Wang X P, Wu X, Dai S, He J S, Yin W L. Relative effects of phylogeny, biological characters and environments on leaf traits in shrub biomes across central Inner Mongolia, China. Journal of Plant Ecology, 2013, 6(3): 220- 231.

    [22]Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H, Diemer M. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821- 827.

    [23]He J S, Wang X P, Schmid B, Flynn D F B, Li X F, Reich P B, Fang J Y. Taxonomic identity, phylogeny, climate and soil fertility as drivers of leaf traits across Chinese grassland biomes. Journal of Plant Research, 2010, 123(4): 551- 561.

    [24]Santiago L S. Extending the leaf economics spectrum to decomposition: evidence from a tropical forest. Ecology, 2007, 88(5): 1126- 1131.

    [25]Chapin III F S, Zavaleta E S, Eviner V T, Naylor R L, Vitousek P M, Reynolds H L, Hooper D U, Lavorel S, Sala O E, Hobbie S E. Consequences of changing biodiversity. Nature, 2000, 405(6783): 234- 242.

    [26]Reitalu T, Sepp? H, Sugita S, Kangur M, Koff T, Avel E, Kihno K, Vassiljev J, Renssen H, Hammarlund D. Long- term drivers of forest composition in a boreonemoral region: the relative importance of climate and human impact. Journal of Biogeography, 2013, 40(8): 1524- 1534.

    [27]Prescott C E. The influence of the forest canopy on nutrient cycling. Tree Physiology, 2002, 22(15/16): 1193- 1200.

    [28]Prescott C E, Hope G D, Blevins L L. Effect of gap size on litter decomposition and soil nitrate concentrations in a high- elevation spruce fir forest. Canadian Journal of Forest Research, 2003, 33(11): 2210- 2220.

    [29]McCarthy J. Gap dynamics of forest trees: a review with particular attention to boreal forests. Environmental Reviews, 2001, 9(1): 1- 59.

    [30]Wang X P, Fang J Y, Sanders N J, White P S, Tang Z Y. Relative importance of climate vs local factors in shaping the regional patterns of forest plant richness across northeast China. Ecography, 2009, 32(1): 133- 142.

    [31]張新平, 王襄平, 朱彪, 宗占江, 彭長輝, 方精云. 我國東北主要森林類型的凋落物產(chǎn)量及其影響因素. 植物生態(tài)學(xué)報(bào), 2008, 32(5): 1031- 1040.

    [32]Gartner T B, Cardon Z G. Decomposition dynamics in mixed-species leaf litter. Oikos, 2004, 104(2): 230- 246.

    [33]Olson J S. Energy storage and the balance of producers and decomposers in ecological systems. Ecology, 1963, 44(2): 322- 331.

    [34]張德強(qiáng), 葉萬輝, 余清發(fā), 孔國輝, 張佑倡. 鼎湖山演替系列中代表性森林凋落物研究. 生態(tài)學(xué)報(bào), 2000, 20 (6): 938- 944.

    [35]黃建輝, 陳靈芝, 韓興國. 森林生態(tài)系統(tǒng)中凋落物分解的研究進(jìn)展 // 李承森. 植物科學(xué)進(jìn)展Vol. 1. 北京: 高等教育出版社, 1998: 218- 316.

    [36]Ndagurwa H G T, Dube J S, Mlambo D. Decomposition and nutrient release patterns of mistletoe litters in a semi-arid savanna, southwest Zimbabwe. Austral Ecology, 2015, 40(2): 178-185.

    [37]方精云. 地理要素對我國溫度分布影響的數(shù)量評價(jià). 生態(tài)學(xué)報(bào), 1992, 12(2): 97- 104.

    [38]Wang X P, Tang Z Y, Fang J Y. Climatic control on forests and tree species distribution in the forest region of Northeast China. Journal of Integrative Plant Biology, 2006, 48(7): 778- 789.

    [39]Heikkinen R K, Luoto M, Kuussaari M, P?yry J. New insights into butterfly-environment relationships using partitioning methods. Proceedings of the Royal Society B: Biological Sciences, 2005, 272(1577): 2203- 2210.

    [40]Wang X P, Tang Z Y, Shen Z H, Zheng C Y, Luo J C, Fang J Y. Relative influence of regional species richness vs local climate on local species richness in China′s forests. Ecography, 2012, 35(12): 1176- 1184.

    [41]劉穎, 武耀祥, 韓士杰, 林鹿. 長白山四種森林類型凋落物分解動(dòng)態(tài). 生態(tài)學(xué)雜志, 2009, 28(3): 400- 404.

    [42]Trofymow J A, Moore T R, Titus B, Prescott C, Morrison I, Siltanen M, Smith S, Fyles J, Wein R, Camiré C, Duschene L, Kozak L, Kranabetter M, Visser S. Rates of litter decomposition over 6 years in Canadian forests: influence of litter quality and climate. Canadian Journal of Forest Research, 2002, 32(5): 789- 804.

    [43]Gholz H L, Wedin D A, Smitherman S M, Harmon M E, Parton W J. Long-term dynamics of pine and hardwood litter in contrasting environments: toward a global model of decomposition. Global Change Biology, 2000, 6(7): 751- 765.

    [44]Vitousek P M, Turner D R, Parton W J, Sanford R L. Litter decomposition on the Mauna Loa environmental matrix, Hawai′i: patterns, mechanisms, and models. Ecology, 1994, 75(2): 418- 429.

    [45]Latter P M, Howson G, Howard D M, Scott W A. Long- term study of litter decomposition on a Pennine peat bog: which regression?. Oecologia, 1997, 113(1): 94- 103.

    [46]Austin A T, Vitousek P M. Nutrient dynamics on a precipitation gradient in Hawai′i. Oecologia, 1998, 113(4): 519- 529.

    [47]Zhang Q S, Zak J C. Effects of gap size on litter decomposition and microbial activity in a subtropical forest. Ecology, 1995, 76(7): 2196- 2204.

    [48]Ritter E. Litter decomposition and nitrogen mineralization in newly formed gaps in a Danish beech (Fagussylvatica) forest. Soil Biology and Biochemistry, 2005, 37(7): 1237- 1247.

    [49]Terborgh J. The vertical component of plant species diversity in temperate and tropical forests. The American Naturalist, 1985, 126(6): 760- 776.

    [50]周以良. 中國東北植被地理. 北京: 科學(xué)出版社, 1997: 289- 299.

    [51]Sariyildiz T. Effects of gap-size classes on long-term litter decomposition rates of beech, oak and chestnut species at high elevations in Northeast Turkey. Ecosystems, 2008, 11(6): 841- 853.

    [52]Prescott C E, Zabek L M, Staley C L, Kabzems R. Decomposition of broadleaf and needle litter in forests of British Columbia: influences of litter type, forest type, and litter mixtures. Canadian Journal of Forest Research, 2000, 30(11): 1742- 1750.

    [53]Perry D A, Choquette C, Schroeder P. Nitrogen dynamics in conifer-dominated forests with and without hardwoods. Canadian Journal of Forest Research, 1987, 17(11): 1434- 1441.

    [54]Anderson J M. Spatiotemporal effects of invertebrates on soil processes. Biology and Fertility of Soils, 1988, 6(3): 216- 227.

    [55]Murphy K L, Klopatek J M, Klopatek C C. The effects of litter quality and climate on decomposition along an elevational gradient. Ecological Applications, 1998, 8(4): 1061- 1071.

    [56]Vivanco L, Austin A T. Tree species identity alters forest litter decomposition through long- term plant and soil interactions in Patagonia, Argentina. Journal of Ecology, 2008, 96(4): 727- 736.

    Effects of climate, forest type and light availability on litter decomposition rate in forests of Northeast China

    WU Peng1, WANG Xiangping1, *, ZHANG Xinping2, ZHU Biao2, ZHOU Haicheng3, FANG Jingyun2

    1KeyLaboratoryforForestResources&EcosystemProcessesofBeijing,BeijingForestryUniversity,Beijing100083,China2DepartmentofEcologyandKeyLaboratoryforEarthSurfaceProcessesoftheMinistryofEducation,PekingUniversity,Beijing100871,China3CenterforNaturalprotectionandmanagementoftheMt.Changbaireserve,Antu133613,China

    Abstract:Decomposition of plant litter is a key process in nutrient and carbon cycling in terrestrial ecosystems. Disentangling the relative effects of different factors, including climate, local environmental factors, and species identities (which determine litter quality), on leaf litter decomposition rate provides insight into forest carbon cycles. In this study, we conducted a three-year decomposition experiment of leaf litter in 26 plots across Northeast China (Changbai, Zhuangguangcai, Xiaoxing′an, and Daxing′an Mountains) to examine the relative effects of climate, forest type, and gap fraction on leaf litter decomposition rate. The leaf litter decomposition rate differed among forest types in the following order: deciduous broadleaf forest > deciduous broadleaf and needle-leaf mixed forest > deciduous needle-leaf forest > evergreen needle-leaf forest > Betula ermanii forest. A multivariate analysis showed that climate factors (energy and water availability) were strong predictors of the leaf decomposition rate, which explained 55.5% and 65.0% of the variation in decomposition coefficient (k) and the time required for decomposition of 95% of the litter (t95%), respectively. However, the independent effect of climate on decomposition was weak, and most of the variation was attributed to the joint effect of climate variables with forest type and light availability (<9%). Climate drives changes in forest types (species composition), which in turn affects the leaf decomposition rate. This joint effect explained 46.8% (k) and 56.8% (t95%) of the total variation in leaf decomposition rate. Moreover, changes in climate and forest types can affect under-canopy light availability (measured using the gap fraction), which declined with increasing temperature and indirectly improved leaf decomposition rate. The joint effects of climate, forest type, and light availability explained 23.9% and 22.3% of the variation in k and t95%, respectively. Our results indicate that climate affects litter decomposition mainly through its indirect effects on biological factors such as species composition and canopy structure, which influence litter quality and light availability. It is necessary to understand the interactions of climate with these biological factors for predicting changes in litter decomposition in response to future climate change.

    Key Words:leaf litter decomposition rate; climate; forest type; gap fraction; Northeast China

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(31370620, 31321061, 31270656); 國家留學(xué)基金項(xiàng)目(2011811457)

    收稿日期:2014- 10- 10; 網(wǎng)絡(luò)出版日期:2015- 05- 20

    *通訊作者

    Corresponding author.E-mail: wangxiangping@bjfu.edu.cn

    DOI:10.5846/stxb201410101991

    吳鵬, 王襄平, 張新平, 朱彪, 周海城, 方精云.東北地區(qū)森林凋落葉分解速率與氣候、林型、林分光照的關(guān)系.生態(tài)學(xué)報(bào),2016,36(8):2223- 2232.

    Wu P, Wang X P, Zhang X P, Zhu B, Zhou H C, Fang J Y.Effects of climate, forest type and light availability on litter decomposition rate in forests of Northeast China.Acta Ecologica Sinica,2016,36(8):2223- 2232.

    猜你喜歡
    林型氣候
    不同林分類型葉片穩(wěn)定碳、氮同位素的變化特征
    山東半島昆崳山地區(qū)主要森林類型可燃物垂直分布及影響因子
    不同紅樹林群落結(jié)構(gòu)與植被碳分布
    咦?氣候會(huì)變化?
    瞧,氣候大不同
    氣候變暖會(huì)怎樣?
    都是氣候變暖惹的禍
    寬甸地區(qū)不同林型枯落物及土壤蓄水功能研究
    青海不同林分土壤微生物群落結(jié)構(gòu)(PLFA)*
    立冬
    97超视频在线观看视频| 99久国产av精品国产电影| 丰满饥渴人妻一区二区三| 久久久a久久爽久久v久久| 日韩熟女老妇一区二区性免费视频| 美女大奶头黄色视频| 一本色道久久久久久精品综合| 中文字幕av电影在线播放| 狂野欧美白嫩少妇大欣赏| av黄色大香蕉| 大陆偷拍与自拍| 91国产中文字幕| 一级a做视频免费观看| 婷婷色综合www| 国产深夜福利视频在线观看| 免费大片黄手机在线观看| 日韩人妻高清精品专区| 黑人猛操日本美女一级片| av不卡在线播放| 精品亚洲乱码少妇综合久久| 高清午夜精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片我不卡| 狠狠精品人妻久久久久久综合| 少妇的逼水好多| 交换朋友夫妻互换小说| 一级毛片电影观看| 黄色视频在线播放观看不卡| 亚洲精品一二三| 亚洲成人一二三区av| 岛国毛片在线播放| av福利片在线| 王馨瑶露胸无遮挡在线观看| 老女人水多毛片| 国产成人精品在线电影| 亚洲av成人精品一二三区| 一区二区三区四区激情视频| 亚洲av不卡在线观看| 欧美激情国产日韩精品一区| 久久精品久久久久久噜噜老黄| 最近中文字幕高清免费大全6| 啦啦啦啦在线视频资源| 国产欧美亚洲国产| 精品久久久久久电影网| 一本一本综合久久| h视频一区二区三区| 久久久久久久久久人人人人人人| 久久久久久久久大av| 好男人视频免费观看在线| 欧美成人午夜免费资源| 一级毛片黄色毛片免费观看视频| 国产男女内射视频| 建设人人有责人人尽责人人享有的| 国产又色又爽无遮挡免| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月| 嘟嘟电影网在线观看| 十八禁网站网址无遮挡| 免费播放大片免费观看视频在线观看| 久久 成人 亚洲| 中文字幕制服av| 精品酒店卫生间| 视频区图区小说| 丝袜美足系列| 国产亚洲精品第一综合不卡 | 人体艺术视频欧美日本| 午夜激情久久久久久久| 国产视频内射| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产一区有黄有色的免费视频| 成人亚洲欧美一区二区av| 精品久久久久久电影网| 大香蕉97超碰在线| 久久毛片免费看一区二区三区| 日韩人妻高清精品专区| av在线播放精品| 狠狠精品人妻久久久久久综合| 久久精品人人爽人人爽视色| 五月玫瑰六月丁香| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 日韩中字成人| 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 一级爰片在线观看| 2022亚洲国产成人精品| 久久99热这里只频精品6学生| 亚洲国产成人一精品久久久| 久久国产精品大桥未久av| 亚洲av男天堂| 日韩 亚洲 欧美在线| 热99久久久久精品小说推荐| av免费观看日本| 91久久精品国产一区二区成人| 蜜桃久久精品国产亚洲av| 国产成人精品在线电影| 国产精品免费大片| 91精品国产九色| 女人久久www免费人成看片| 日韩不卡一区二区三区视频在线| 老司机亚洲免费影院| 国产在线一区二区三区精| 丁香六月天网| 日韩欧美精品免费久久| 春色校园在线视频观看| 国产精品久久久久久精品古装| 建设人人有责人人尽责人人享有的| 免费黄网站久久成人精品| 在线观看人妻少妇| 美女中出高潮动态图| 午夜福利,免费看| 在线天堂最新版资源| 国产午夜精品久久久久久一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 人人妻人人添人人爽欧美一区卜| av女优亚洲男人天堂| 久久久久久久久久久丰满| 精品人妻熟女毛片av久久网站| a级毛片在线看网站| 最近最新中文字幕免费大全7| av在线播放精品| 一边摸一边做爽爽视频免费| 一级爰片在线观看| 国产男女超爽视频在线观看| 国产午夜精品久久久久久一区二区三区| 国产亚洲欧美精品永久| 亚洲国产精品一区三区| av一本久久久久| 99九九在线精品视频| 精品少妇内射三级| a级毛片黄视频| 少妇精品久久久久久久| 亚洲成人av在线免费| 少妇丰满av| 国产成人av激情在线播放 | 亚洲欧美一区二区三区国产| 午夜福利视频在线观看免费| 免费看光身美女| 午夜福利在线观看免费完整高清在| 日韩一本色道免费dvd| av福利片在线| 久久人人爽人人爽人人片va| 免费大片18禁| 一级毛片我不卡| 下体分泌物呈黄色| 高清毛片免费看| 麻豆成人av视频| 亚洲激情五月婷婷啪啪| 各种免费的搞黄视频| 91精品伊人久久大香线蕉| 啦啦啦啦在线视频资源| 国产又色又爽无遮挡免| 纵有疾风起免费观看全集完整版| 亚洲国产欧美日韩在线播放| 国产又色又爽无遮挡免| 国产精品 国内视频| 久久人人爽人人爽人人片va| 2021少妇久久久久久久久久久| 精品99又大又爽又粗少妇毛片| 精品一品国产午夜福利视频| 欧美日本中文国产一区发布| 狠狠婷婷综合久久久久久88av| 男的添女的下面高潮视频| av免费在线看不卡| 一区二区三区四区激情视频| 丝袜脚勾引网站| 久久久久久久亚洲中文字幕| a级片在线免费高清观看视频| 久久鲁丝午夜福利片| 色94色欧美一区二区| 国产男女超爽视频在线观看| 国产精品无大码| 91国产中文字幕| 中国三级夫妇交换| 男男h啪啪无遮挡| 美女内射精品一级片tv| 亚洲婷婷狠狠爱综合网| 青春草国产在线视频| 国产在线一区二区三区精| 欧美激情极品国产一区二区三区 | 国产永久视频网站| 2018国产大陆天天弄谢| 最近手机中文字幕大全| 纯流量卡能插随身wifi吗| 久久久欧美国产精品| 日本黄大片高清| 亚洲精品日韩av片在线观看| 一区二区av电影网| 精品久久国产蜜桃| 国产成人免费无遮挡视频| 亚洲欧美中文字幕日韩二区| 午夜福利视频精品| 国产精品久久久久久精品古装| 国产不卡av网站在线观看| 久久久久久久久久久久大奶| 精品久久久久久久久av| 老司机亚洲免费影院| 九九爱精品视频在线观看| 日韩视频在线欧美| 伦理电影大哥的女人| 午夜福利在线观看免费完整高清在| 超色免费av| 亚洲av.av天堂| 狂野欧美激情性xxxx在线观看| 在现免费观看毛片| 欧美精品高潮呻吟av久久| 人体艺术视频欧美日本| 日日撸夜夜添| 婷婷色综合www| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 99国产精品免费福利视频| 人妻夜夜爽99麻豆av| 欧美亚洲 丝袜 人妻 在线| 黑人高潮一二区| 精品国产国语对白av| 男人添女人高潮全过程视频| 如日韩欧美国产精品一区二区三区 | 在线观看国产h片| 一本大道久久a久久精品| 少妇 在线观看| av在线app专区| 色94色欧美一区二区| 人人妻人人澡人人看| 精品一区在线观看国产| 亚洲四区av| 最近的中文字幕免费完整| 精品人妻在线不人妻| 亚洲一区二区三区欧美精品| 国产片内射在线| 国产69精品久久久久777片| 久热这里只有精品99| 肉色欧美久久久久久久蜜桃| 91aial.com中文字幕在线观看| 桃花免费在线播放| av天堂久久9| 男人爽女人下面视频在线观看| 午夜福利视频在线观看免费| 女人久久www免费人成看片| 亚洲欧美成人综合另类久久久| 91精品国产国语对白视频| 在线亚洲精品国产二区图片欧美 | 久久影院123| 色婷婷久久久亚洲欧美| 欧美成人精品欧美一级黄| 午夜免费男女啪啪视频观看| 成年美女黄网站色视频大全免费 | 久热久热在线精品观看| 国产在线免费精品| 亚洲欧美成人综合另类久久久| 一区二区三区乱码不卡18| 国产一区有黄有色的免费视频| 69精品国产乱码久久久| 亚洲av欧美aⅴ国产| 亚洲国产av新网站| 一二三四中文在线观看免费高清| 99热全是精品| 久久女婷五月综合色啪小说| 亚洲伊人久久精品综合| 激情五月婷婷亚洲| 最近手机中文字幕大全| 亚洲精品一区蜜桃| 九草在线视频观看| 91国产中文字幕| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 午夜免费男女啪啪视频观看| 两个人的视频大全免费| 日韩不卡一区二区三区视频在线| videossex国产| 久久久精品免费免费高清| 久久精品久久久久久噜噜老黄| 2021少妇久久久久久久久久久| 在线播放无遮挡| 日本wwww免费看| 99热6这里只有精品| 国产免费现黄频在线看| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久久久免| 夜夜看夜夜爽夜夜摸| 亚洲精品国产色婷婷电影| 纯流量卡能插随身wifi吗| 欧美激情国产日韩精品一区| 日本免费在线观看一区| 亚洲一级一片aⅴ在线观看| 青青草视频在线视频观看| 亚洲综合精品二区| 中文字幕制服av| 国产精品.久久久| 国产成人一区二区在线| 哪个播放器可以免费观看大片| 大片免费播放器 马上看| 亚洲国产毛片av蜜桃av| 亚洲第一区二区三区不卡| 日韩伦理黄色片| 18禁在线无遮挡免费观看视频| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 97超视频在线观看视频| 人妻人人澡人人爽人人| 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 哪个播放器可以免费观看大片| 黄色毛片三级朝国网站| 亚洲国产毛片av蜜桃av| 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 97在线人人人人妻| 男女啪啪激烈高潮av片| xxxhd国产人妻xxx| 中文字幕免费在线视频6| 日韩熟女老妇一区二区性免费视频| 久久精品国产自在天天线| 国产成人免费无遮挡视频| 高清毛片免费看| 人妻 亚洲 视频| 亚洲av日韩在线播放| 亚洲精品色激情综合| 欧美xxxx性猛交bbbb| 免费观看在线日韩| 国产成人免费无遮挡视频| 久久久国产欧美日韩av| 亚洲欧美精品自产自拍| 大片电影免费在线观看免费| 亚洲精品日本国产第一区| 成人综合一区亚洲| 国产女主播在线喷水免费视频网站| 欧美bdsm另类| 久久青草综合色| 欧美成人午夜免费资源| 国产av精品麻豆| 国产亚洲最大av| 涩涩av久久男人的天堂| 日本-黄色视频高清免费观看| 菩萨蛮人人尽说江南好唐韦庄| 国产成人免费无遮挡视频| 99久久人妻综合| 国产av精品麻豆| 国产精品人妻久久久影院| 亚洲人与动物交配视频| 两个人的视频大全免费| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 99国产精品免费福利视频| 国产熟女午夜一区二区三区 | 久久99热6这里只有精品| 成人综合一区亚洲| a级毛片在线看网站| 制服人妻中文乱码| 久热久热在线精品观看| 免费人妻精品一区二区三区视频| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠久久av| 男人爽女人下面视频在线观看| 国产成人freesex在线| 桃花免费在线播放| 久久久久网色| 久久久精品免费免费高清| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 亚洲成人一二三区av| 精品国产一区二区三区久久久樱花| 国产熟女午夜一区二区三区 | 精品一区二区三卡| 国产精品一区二区在线观看99| 亚洲一级一片aⅴ在线观看| av女优亚洲男人天堂| 成年人免费黄色播放视频| 日本爱情动作片www.在线观看| 国产成人一区二区在线| 男女国产视频网站| 欧美成人午夜免费资源| 国产男女超爽视频在线观看| 一区在线观看完整版| 在线观看免费高清a一片| 亚洲精品久久成人aⅴ小说 | 国产国语露脸激情在线看| 在现免费观看毛片| 久久鲁丝午夜福利片| 丰满乱子伦码专区| av黄色大香蕉| 蜜臀久久99精品久久宅男| 国产黄色视频一区二区在线观看| 最近中文字幕2019免费版| 激情五月婷婷亚洲| 一级黄片播放器| 99国产精品免费福利视频| 成人国产麻豆网| 丰满少妇做爰视频| 天天影视国产精品| av在线观看视频网站免费| 午夜福利视频精品| 国产精品国产av在线观看| 国产精品欧美亚洲77777| 久热这里只有精品99| 伦理电影免费视频| 日日摸夜夜添夜夜添av毛片| 观看美女的网站| 免费看光身美女| 91精品国产国语对白视频| 国产成人精品在线电影| 在线观看美女被高潮喷水网站| av在线老鸭窝| 18在线观看网站| 夜夜骑夜夜射夜夜干| 秋霞伦理黄片| 最黄视频免费看| 性色avwww在线观看| 成人午夜精彩视频在线观看| av免费观看日本| 美女国产高潮福利片在线看| 成人免费观看视频高清| 亚洲精品av麻豆狂野| 久久精品人人爽人人爽视色| 人妻人人澡人人爽人人| 成人毛片60女人毛片免费| 精品人妻熟女av久视频| 亚洲精品久久久久久婷婷小说| 性色av一级| 精品久久久久久久久亚洲| 香蕉精品网在线| 日韩人妻高清精品专区| 一二三四中文在线观看免费高清| 日本91视频免费播放| 日韩人妻高清精品专区| 日本爱情动作片www.在线观看| 老熟女久久久| 欧美xxⅹ黑人| 久久久久人妻精品一区果冻| 狠狠精品人妻久久久久久综合| 国产免费一级a男人的天堂| av电影中文网址| 国产av码专区亚洲av| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站| 狂野欧美白嫩少妇大欣赏| 欧美精品一区二区免费开放| 国产爽快片一区二区三区| 如日韩欧美国产精品一区二区三区 | 国产精品秋霞免费鲁丝片| 老司机影院成人| 尾随美女入室| 少妇的逼好多水| 夜夜爽夜夜爽视频| 精品一区在线观看国产| 人妻 亚洲 视频| 久久人人爽人人爽人人片va| 国产精品国产三级国产专区5o| 日韩亚洲欧美综合| 精品人妻熟女毛片av久久网站| 久久影院123| 国产在线免费精品| 大香蕉97超碰在线| 成人毛片60女人毛片免费| 亚洲高清免费不卡视频| 伊人久久精品亚洲午夜| 又粗又硬又长又爽又黄的视频| 纯流量卡能插随身wifi吗| 99精国产麻豆久久婷婷| 亚洲欧洲日产国产| 亚洲国产精品国产精品| 国产免费一级a男人的天堂| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 99热国产这里只有精品6| 亚洲av不卡在线观看| 欧美激情 高清一区二区三区| 亚洲伊人久久精品综合| 国产精品一二三区在线看| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 乱码一卡2卡4卡精品| 天天躁夜夜躁狠狠久久av| 精品卡一卡二卡四卡免费| 久久亚洲国产成人精品v| 精品久久久久久久久av| 新久久久久国产一级毛片| 成人黄色视频免费在线看| 亚洲成人一二三区av| 国产伦精品一区二区三区视频9| 欧美一级a爱片免费观看看| 久久影院123| 色94色欧美一区二区| 成年人午夜在线观看视频| 亚洲激情五月婷婷啪啪| 简卡轻食公司| 成人毛片a级毛片在线播放| 午夜老司机福利剧场| a级片在线免费高清观看视频| 久久久久视频综合| 99热这里只有精品一区| 一区二区三区精品91| 日韩免费高清中文字幕av| 免费看光身美女| 97在线视频观看| av在线app专区| 各种免费的搞黄视频| 97超碰精品成人国产| 毛片一级片免费看久久久久| 午夜福利网站1000一区二区三区| 欧美亚洲日本最大视频资源| 亚洲精品色激情综合| 日日摸夜夜添夜夜添av毛片| 国产高清国产精品国产三级| 国产av一区二区精品久久| 中文精品一卡2卡3卡4更新| 国产日韩欧美亚洲二区| 日本色播在线视频| 丝袜脚勾引网站| 日本免费在线观看一区| 亚洲高清免费不卡视频| 欧美人与善性xxx| 久久久久人妻精品一区果冻| 青春草视频在线免费观看| 亚洲av福利一区| 新久久久久国产一级毛片| 丰满乱子伦码专区| 免费av不卡在线播放| 夫妻午夜视频| videos熟女内射| 视频区图区小说| 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 色婷婷久久久亚洲欧美| 国产精品一区二区三区四区免费观看| 夫妻午夜视频| 极品少妇高潮喷水抽搐| 日韩av在线免费看完整版不卡| 久久 成人 亚洲| 91精品一卡2卡3卡4卡| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区三区av在线| 水蜜桃什么品种好| 久久人妻熟女aⅴ| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 中文字幕最新亚洲高清| 成年女人在线观看亚洲视频| 一级片'在线观看视频| 一区二区日韩欧美中文字幕 | 亚洲国产av新网站| 国产精品久久久久久久久免| 汤姆久久久久久久影院中文字幕| a 毛片基地| 欧美精品一区二区大全| 女性被躁到高潮视频| 国产成人精品久久久久久| 曰老女人黄片| 国产黄色免费在线视频| 久久99热6这里只有精品| 亚洲精品,欧美精品| 国产日韩欧美视频二区| 我的老师免费观看完整版| 日韩亚洲欧美综合| 国产免费福利视频在线观看| 午夜老司机福利剧场| .国产精品久久| 啦啦啦视频在线资源免费观看| 久久久久久久久久成人| 久久综合国产亚洲精品| 亚洲久久久国产精品| 国产精品久久久久久av不卡| 日韩,欧美,国产一区二区三区| 毛片一级片免费看久久久久| 久久女婷五月综合色啪小说| 91精品国产国语对白视频| 久久鲁丝午夜福利片| av在线播放精品| 岛国毛片在线播放| av卡一久久| 久久女婷五月综合色啪小说| 午夜91福利影院| 永久免费av网站大全| 国产爽快片一区二区三区| 91久久精品国产一区二区三区| 欧美日韩综合久久久久久| 妹子高潮喷水视频| 免费播放大片免费观看视频在线观看| videos熟女内射| 啦啦啦在线观看免费高清www| 中文字幕久久专区| 在线观看三级黄色| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| 亚洲精华国产精华液的使用体验| 丰满迷人的少妇在线观看| 久久久久久久久久成人| 91精品国产国语对白视频| 国产男女超爽视频在线观看| 亚洲av免费高清在线观看| 欧美精品一区二区大全| 国产精品国产三级国产专区5o| 国产熟女午夜一区二区三区 | 成人无遮挡网站| 在线看a的网站| 下体分泌物呈黄色| 91国产中文字幕| 狠狠精品人妻久久久久久综合| 精品一品国产午夜福利视频| 亚洲欧洲精品一区二区精品久久久 | 91精品国产国语对白视频| 久久久久久久久久久久大奶| 久久久久久久久久久丰满| 综合色丁香网| 天天躁夜夜躁狠狠久久av| 国产色爽女视频免费观看| 人人澡人人妻人| 22中文网久久字幕| 国产精品久久久久久久久免| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久| 多毛熟女@视频| 麻豆成人av视频|