• <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)*
    立冬
    精品人妻偷拍中文字幕| 九色国产91popny在线| 三级毛片av免费| 一个人免费在线观看电影| 亚洲精品粉嫩美女一区| 国产蜜桃级精品一区二区三区| 中文字幕人妻丝袜一区二区| 天堂动漫精品| 国产一区二区亚洲精品在线观看| 久久久精品欧美日韩精品| 少妇裸体淫交视频免费看高清| 日本a在线网址| 欧美日韩中文字幕国产精品一区二区三区| or卡值多少钱| 最近视频中文字幕2019在线8| 欧美一区二区国产精品久久精品| 亚洲av成人av| 久久精品91蜜桃| 变态另类成人亚洲欧美熟女| 国产极品精品免费视频能看的| 亚洲av成人不卡在线观看播放网| 国产不卡一卡二| 国产欧美日韩一区二区三| 波多野结衣高清作品| 波多野结衣高清作品| 久久久久久大精品| 色综合站精品国产| 久久久久久久精品吃奶| 99在线视频只有这里精品首页| 99久久精品国产亚洲精品| 亚洲精品456在线播放app | 欧美又色又爽又黄视频| 美女大奶头视频| 美女黄网站色视频| 90打野战视频偷拍视频| 免费观看人在逋| 欧美黑人巨大hd| 嫩草影院入口| 亚洲五月婷婷丁香| 无人区码免费观看不卡| 亚洲一区高清亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 久久草成人影院| 欧美另类亚洲清纯唯美| 校园春色视频在线观看| 色尼玛亚洲综合影院| 午夜视频国产福利| 亚洲欧美日韩东京热| av女优亚洲男人天堂| 久久性视频一级片| 成年人黄色毛片网站| 欧美3d第一页| 欧美国产日韩亚洲一区| 免费av观看视频| 在线十欧美十亚洲十日本专区| 悠悠久久av| 日本与韩国留学比较| 97超视频在线观看视频| 久久久精品欧美日韩精品| 亚洲aⅴ乱码一区二区在线播放| 好男人在线观看高清免费视频| 色视频www国产| 亚洲天堂国产精品一区在线| 九色国产91popny在线| 精品久久久久久久人妻蜜臀av| 在线十欧美十亚洲十日本专区| 麻豆久久精品国产亚洲av| 啦啦啦免费观看视频1| 久久精品综合一区二区三区| 最近最新中文字幕大全免费视频| 男插女下体视频免费在线播放| 女警被强在线播放| 成人亚洲精品av一区二区| 国产色爽女视频免费观看| 熟妇人妻久久中文字幕3abv| 日韩中文字幕欧美一区二区| 每晚都被弄得嗷嗷叫到高潮| 村上凉子中文字幕在线| 男人舔女人下体高潮全视频| 美女黄网站色视频| 露出奶头的视频| 亚洲黑人精品在线| 俺也久久电影网| 无人区码免费观看不卡| 法律面前人人平等表现在哪些方面| 淫妇啪啪啪对白视频| 午夜免费男女啪啪视频观看 | 欧美一区二区精品小视频在线| 国产真人三级小视频在线观看| 女生性感内裤真人,穿戴方法视频| 日日干狠狠操夜夜爽| 免费观看的影片在线观看| 久久精品国产99精品国产亚洲性色| 精品久久久久久成人av| 国产 一区 欧美 日韩| 亚洲午夜理论影院| 欧美不卡视频在线免费观看| 久久久久亚洲av毛片大全| 91麻豆精品激情在线观看国产| 欧美性猛交黑人性爽| 日韩欧美精品免费久久 | 成人亚洲精品av一区二区| 亚洲av美国av| 欧美色欧美亚洲另类二区| 午夜福利视频1000在线观看| 一区福利在线观看| 亚洲成人久久性| 久久久久精品国产欧美久久久| 天堂√8在线中文| 午夜a级毛片| 88av欧美| 精品国产美女av久久久久小说| 午夜福利在线在线| 国产精品一区二区三区四区免费观看 | 人妻久久中文字幕网| 国产淫片久久久久久久久 | 欧美一区二区精品小视频在线| 欧美区成人在线视频| 99久久九九国产精品国产免费| 中文字幕高清在线视频| 日本在线视频免费播放| 国产精品美女特级片免费视频播放器| 狠狠狠狠99中文字幕| 亚洲av成人av| 免费看a级黄色片| 一进一出抽搐gif免费好疼| 女同久久另类99精品国产91| 国产精品一区二区三区四区免费观看 | 久9热在线精品视频| 嫁个100分男人电影在线观看| 欧美在线一区亚洲| 51午夜福利影视在线观看| 听说在线观看完整版免费高清| 19禁男女啪啪无遮挡网站| 美女黄网站色视频| 日韩欧美一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 波多野结衣高清无吗| 亚洲色图av天堂| 少妇高潮的动态图| 岛国在线观看网站| 少妇高潮的动态图| 激情在线观看视频在线高清| 内射极品少妇av片p| 久久精品亚洲精品国产色婷小说| 91在线观看av| 国产亚洲精品一区二区www| 成人永久免费在线观看视频| 又黄又粗又硬又大视频| 亚洲欧美日韩东京热| 色综合婷婷激情| 一本综合久久免费| 一进一出抽搐gif免费好疼| 啦啦啦免费观看视频1| 变态另类丝袜制服| 久久精品影院6| 久久久久九九精品影院| 精品午夜福利视频在线观看一区| 亚洲男人的天堂狠狠| 国产精品电影一区二区三区| 成人鲁丝片一二三区免费| 老司机福利观看| 91麻豆av在线| 国产精品野战在线观看| 欧美一区二区精品小视频在线| 久久久久久久午夜电影| 免费观看人在逋| 美女cb高潮喷水在线观看| 夜夜夜夜夜久久久久| 国内精品久久久久精免费| 日日摸夜夜添夜夜添小说| 国产亚洲精品av在线| 国产一级毛片七仙女欲春2| 国产男靠女视频免费网站| 少妇丰满av| 免费高清视频大片| 精品久久久久久久毛片微露脸| 成人无遮挡网站| 熟女少妇亚洲综合色aaa.| 日韩欧美国产在线观看| 久久精品国产综合久久久| 99视频精品全部免费 在线| 欧美日韩综合久久久久久 | 精品久久久久久久毛片微露脸| 午夜老司机福利剧场| 国产亚洲欧美98| 在线播放无遮挡| 久久精品91蜜桃| 午夜福利在线在线| 成人av一区二区三区在线看| 精品人妻偷拍中文字幕| 亚洲专区中文字幕在线| 亚洲激情在线av| 国产一区二区三区在线臀色熟女| 欧美激情在线99| 熟女人妻精品中文字幕| 久久久久精品国产欧美久久久| 精品福利观看| 九九久久精品国产亚洲av麻豆| 欧美成人性av电影在线观看| 国产美女午夜福利| 亚洲精品456在线播放app | 岛国视频午夜一区免费看| 制服丝袜大香蕉在线| 亚洲最大成人手机在线| 精品99又大又爽又粗少妇毛片 | 久久99热这里只有精品18| 国产精品,欧美在线| 亚洲av日韩精品久久久久久密| 天天添夜夜摸| 51午夜福利影视在线观看| 黄色片一级片一级黄色片| 久久久久免费精品人妻一区二区| 波多野结衣高清无吗| 美女黄网站色视频| 天天添夜夜摸| 日韩有码中文字幕| 久久久精品大字幕| 日本黄大片高清| 久久久国产精品麻豆| 在线观看日韩欧美| 精品人妻一区二区三区麻豆 | 高清毛片免费观看视频网站| 久久6这里有精品| 狠狠狠狠99中文字幕| 免费观看人在逋| 69人妻影院| 黄色片一级片一级黄色片| 国产97色在线日韩免费| 两人在一起打扑克的视频| 天堂av国产一区二区熟女人妻| 三级男女做爰猛烈吃奶摸视频| 久久亚洲精品不卡| 长腿黑丝高跟| 国产真实乱freesex| 日本三级黄在线观看| 成人国产综合亚洲| 人妻夜夜爽99麻豆av| 美女高潮喷水抽搐中文字幕| 99国产精品一区二区蜜桃av| svipshipincom国产片| 亚洲av五月六月丁香网| 变态另类成人亚洲欧美熟女| 国产亚洲精品久久久久久毛片| 国产精品98久久久久久宅男小说| 欧美av亚洲av综合av国产av| 亚洲美女黄片视频| 婷婷六月久久综合丁香| 亚洲熟妇熟女久久| 琪琪午夜伦伦电影理论片6080| av在线天堂中文字幕| 亚洲av熟女| 操出白浆在线播放| 小蜜桃在线观看免费完整版高清| 一个人看视频在线观看www免费 | 免费人成视频x8x8入口观看| 美女cb高潮喷水在线观看| 国产单亲对白刺激| 成人一区二区视频在线观看| 国产精品98久久久久久宅男小说| 亚洲真实伦在线观看| 国产一区二区在线观看日韩 | 亚洲av免费高清在线观看| 特大巨黑吊av在线直播| 免费在线观看成人毛片| 国产97色在线日韩免费| 久久久精品大字幕| 一个人看的www免费观看视频| 久久精品人妻少妇| 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 麻豆国产av国片精品| 老司机午夜福利在线观看视频| 真实男女啪啪啪动态图| 91九色精品人成在线观看| 亚洲av五月六月丁香网| 91字幕亚洲| 国产不卡一卡二| 日韩欧美三级三区| 极品教师在线免费播放| 精品熟女少妇八av免费久了| 一本精品99久久精品77| 日本熟妇午夜| 午夜福利在线在线| 男女床上黄色一级片免费看| 国产一区二区三区在线臀色熟女| 久久性视频一级片| 亚洲男人的天堂狠狠| 99久久99久久久精品蜜桃| 手机成人av网站| 久久人人精品亚洲av| 欧美日本亚洲视频在线播放| 有码 亚洲区| 高清毛片免费观看视频网站| 精品久久久久久成人av| 天堂网av新在线| 亚洲中文字幕日韩| 日韩精品青青久久久久久| 亚洲av美国av| 中国美女看黄片| 校园春色视频在线观看| 欧美乱码精品一区二区三区| 精品国产亚洲在线| 国产野战对白在线观看| 变态另类丝袜制服| 天美传媒精品一区二区| 精品熟女少妇八av免费久了| 日本一本二区三区精品| 99久久九九国产精品国产免费| 亚洲自拍偷在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲五月婷婷丁香| 99久久精品国产亚洲精品| 国产又黄又爽又无遮挡在线| 国产探花极品一区二区| 午夜两性在线视频| 亚洲精品成人久久久久久| 亚洲av电影在线进入| 麻豆国产av国片精品| 两个人看的免费小视频| 欧美日本视频| 亚洲av成人不卡在线观看播放网| АⅤ资源中文在线天堂| 国产精品 国内视频| 在线观看免费午夜福利视频| 人妻丰满熟妇av一区二区三区| 日韩中文字幕欧美一区二区| 搡老岳熟女国产| 搡老妇女老女人老熟妇| 欧美黄色淫秽网站| 欧美日韩国产亚洲二区| 午夜福利在线在线| 久久久精品欧美日韩精品| 国产中年淑女户外野战色| 国产精品野战在线观看| 久久久久久人人人人人| 女人高潮潮喷娇喘18禁视频| 国产精品三级大全| 亚洲最大成人手机在线| 国产在视频线在精品| 黄色日韩在线| 欧美av亚洲av综合av国产av| 黄色视频,在线免费观看| aaaaa片日本免费| 久久久精品大字幕| h日本视频在线播放| 亚洲美女黄片视频| 国内毛片毛片毛片毛片毛片| 欧美3d第一页| 日本一本二区三区精品| 又紧又爽又黄一区二区| 亚洲美女黄片视频| 亚洲欧美日韩东京热| 欧美成人性av电影在线观看| 亚洲av免费高清在线观看| 男女那种视频在线观看| 亚洲精品在线美女| 十八禁网站免费在线| 毛片女人毛片| 欧美中文综合在线视频| 韩国av一区二区三区四区| 村上凉子中文字幕在线| 日本a在线网址| 日本黄色片子视频| 国产精品亚洲美女久久久| 国产午夜精品久久久久久一区二区三区 | 人妻久久中文字幕网| 成人永久免费在线观看视频| 丁香欧美五月| 啪啪无遮挡十八禁网站| 国产在视频线在精品| 欧美日韩黄片免| 国产三级黄色录像| 在线看三级毛片| 人人妻人人看人人澡| 好男人电影高清在线观看| 亚洲狠狠婷婷综合久久图片| www国产在线视频色| 老熟妇乱子伦视频在线观看| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 成人精品一区二区免费| 日本与韩国留学比较| 99国产精品一区二区蜜桃av| 免费看日本二区| 天堂网av新在线| 成人国产一区最新在线观看| 18禁在线播放成人免费| 亚洲人成网站在线播| 女同久久另类99精品国产91| 中文字幕人妻丝袜一区二区| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| avwww免费| 99久久无色码亚洲精品果冻| 中文字幕人妻熟人妻熟丝袜美 | 久久精品91无色码中文字幕| 亚洲av成人精品一区久久| 国产一区二区三区视频了| 久久久久久久午夜电影| 亚洲人成伊人成综合网2020| a级一级毛片免费在线观看| 中文字幕av成人在线电影| av在线天堂中文字幕| 亚洲国产中文字幕在线视频| 久久精品91蜜桃| 淫秽高清视频在线观看| 91麻豆av在线| 夜夜看夜夜爽夜夜摸| 久久久国产精品麻豆| 欧美成人免费av一区二区三区| 女同久久另类99精品国产91| 久久香蕉国产精品| 老司机午夜福利在线观看视频| 日本成人三级电影网站| 亚洲18禁久久av| 久久久久性生活片| h日本视频在线播放| 久久精品影院6| 性欧美人与动物交配| 国产视频一区二区在线看| 久久久久久国产a免费观看| 少妇的丰满在线观看| 国产毛片a区久久久久| 国产高清视频在线观看网站| 99久久精品国产亚洲精品| 波野结衣二区三区在线 | 制服丝袜大香蕉在线| 亚洲第一电影网av| 国产精品综合久久久久久久免费| 我要搜黄色片| 九色国产91popny在线| 亚洲av成人不卡在线观看播放网| 69av精品久久久久久| 一级毛片女人18水好多| 免费看a级黄色片| 一边摸一边抽搐一进一小说| 国产精品电影一区二区三区| 国产野战对白在线观看| 久久天躁狠狠躁夜夜2o2o| 免费在线观看日本一区| 亚洲美女黄片视频| 欧美成人一区二区免费高清观看| 脱女人内裤的视频| 久久国产乱子伦精品免费另类| 在线观看午夜福利视频| 精品久久久久久久人妻蜜臀av| 在线观看一区二区三区| 午夜精品久久久久久毛片777| 免费看光身美女| 亚洲国产欧洲综合997久久,| 欧美中文日本在线观看视频| 最新美女视频免费是黄的| 亚洲成人久久性| 国产高清有码在线观看视频| 99久久精品一区二区三区| 免费av毛片视频| 欧美极品一区二区三区四区| 久久国产乱子伦精品免费另类| 非洲黑人性xxxx精品又粗又长| 国产乱人视频| 毛片女人毛片| 女同久久另类99精品国产91| 国产高清有码在线观看视频| 欧美一级a爱片免费观看看| 国产不卡一卡二| av片东京热男人的天堂| 全区人妻精品视频| 俄罗斯特黄特色一大片| 国产精品野战在线观看| 一进一出好大好爽视频| 日本黄色片子视频| 午夜免费成人在线视频| 在线观看免费午夜福利视频| 欧美乱色亚洲激情| 国产精品av视频在线免费观看| 麻豆国产av国片精品| 国产熟女xx| 他把我摸到了高潮在线观看| 成人性生交大片免费视频hd| 亚洲va日本ⅴa欧美va伊人久久| 久久99热这里只有精品18| 精品日产1卡2卡| 午夜免费成人在线视频| 国产成人aa在线观看| 国产亚洲精品久久久久久毛片| 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 久久国产乱子伦精品免费另类| 欧美黑人巨大hd| 99热精品在线国产| 91av网一区二区| 青草久久国产| 国产午夜精品论理片| 高清日韩中文字幕在线| 免费在线观看成人毛片| 在线a可以看的网站| 丰满乱子伦码专区| 精品99又大又爽又粗少妇毛片 | 久久久国产精品麻豆| 欧美中文综合在线视频| 日本免费一区二区三区高清不卡| av视频在线观看入口| 看免费av毛片| 国产精品av视频在线免费观看| 在线看三级毛片| 每晚都被弄得嗷嗷叫到高潮| 中国美女看黄片| 久久久久国产精品人妻aⅴ院| 天天添夜夜摸| 国产av麻豆久久久久久久| 宅男免费午夜| 法律面前人人平等表现在哪些方面| 国产野战对白在线观看| 亚洲成人久久爱视频| 动漫黄色视频在线观看| 欧美成人a在线观看| 伊人久久精品亚洲午夜| 嫩草影院精品99| av专区在线播放| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 亚洲精品成人久久久久久| 在线看三级毛片| 好男人电影高清在线观看| 1024手机看黄色片| 日本黄色视频三级网站网址| 中亚洲国语对白在线视频| www日本在线高清视频| 香蕉丝袜av| 色综合婷婷激情| 国产av不卡久久| 国产成人av激情在线播放| 国产91精品成人一区二区三区| 亚洲久久久久久中文字幕| av国产免费在线观看| 国产成人欧美在线观看| 国产男靠女视频免费网站| 少妇人妻精品综合一区二区 | 夜夜爽天天搞| 欧美精品啪啪一区二区三区| 我要搜黄色片| 搡老熟女国产l中国老女人| 露出奶头的视频| 色综合婷婷激情| 色综合站精品国产| 国产伦精品一区二区三区四那| 哪里可以看免费的av片| 国产 一区 欧美 日韩| 91麻豆精品激情在线观看国产| 国内久久婷婷六月综合欲色啪| 久久亚洲精品不卡| 精品久久久久久成人av| 最新中文字幕久久久久| 精品日产1卡2卡| 色噜噜av男人的天堂激情| 成人亚洲精品av一区二区| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 欧美bdsm另类| 亚洲男人的天堂狠狠| 一区福利在线观看| 内地一区二区视频在线| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 亚洲美女视频黄频| 高清在线国产一区| www.色视频.com| 中文字幕人妻熟人妻熟丝袜美 | 法律面前人人平等表现在哪些方面| 少妇丰满av| 亚洲不卡免费看| 99久久精品热视频| 级片在线观看| 免费电影在线观看免费观看| xxx96com| 日本 av在线| 亚洲第一电影网av| 国产精品,欧美在线| 久久久久久久亚洲中文字幕 | 九色成人免费人妻av| 观看免费一级毛片| 丰满乱子伦码专区| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 亚洲精品一卡2卡三卡4卡5卡| 一区二区三区免费毛片| 熟女电影av网| 国产精华一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲人与动物交配视频| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| 美女大奶头视频| 欧美极品一区二区三区四区| 男女做爰动态图高潮gif福利片| 午夜久久久久精精品| 天堂动漫精品| 一级毛片高清免费大全| 香蕉丝袜av| 成年版毛片免费区| a在线观看视频网站| av在线蜜桃| 在线免费观看不下载黄p国产 | av片东京热男人的天堂| 18禁黄网站禁片免费观看直播| 国产久久久一区二区三区| 亚洲,欧美精品.| 国产激情欧美一区二区| 天天添夜夜摸| 成人国产综合亚洲| 久久久久九九精品影院| 国产高清视频在线播放一区|