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

    土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學(xué)計量特征的影響

    2015-01-20 02:47:42周正虎王傳寬張全智
    生態(tài)學(xué)報 2015年20期
    關(guān)鍵詞:氮磷土地利用計量

    周正虎,王傳寬,張全智

    東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

    土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學(xué)計量特征的影響

    周正虎,王傳寬*,張全智

    東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

    土地利用方式的改變打破森林生態(tài)系統(tǒng)原有的碳氮磷平衡,從而顯著地影響森林生態(tài)系統(tǒng)的生物地球化學(xué)循環(huán)過程。以地段相鄰、林齡相同(10年生)、原始植被一致但土地利用方式不同(無土壤翻動的天然次生林[NS]、間作大豆而土壤翻動中等的人工林[MS]、間作人參而土壤翻動嚴重的次生林[SS])的溫帶幼齡林為對象,探索土地利用變化對土壤碳、氮、磷含量及相互關(guān)系的影響。結(jié)果顯示:(1)土地利用方式顯著改變表層和深層土壤碳含量,各土壤層次碳含量均呈現(xiàn)NS>MS>SS;而氮含量僅在0—20 cm具有顯著性差異(P<0.05);不同土地利用類型之間磷含量無顯著差異(P>0.05);表明碳氮磷對土地利用變化敏感程度不同。(2)SS土壤碳氮比(C/N)和碳磷比(C/P)低于NS和MS,而NS和MS之間C/N和C/P因土壤層次而異。不同土壤層次氮磷比(N/P)均隨土壤翻動強度的增加而顯著減小(NS>MS>SS,P<0.05),且隨土層加深而降低;表明N/P相對于C/N和C/P可能對土地利用變化具有更優(yōu)生態(tài)指示功能。(3)土地利用變化顯著改變土壤碳-氮、碳-磷、氮-磷的耦合關(guān)系。土壤碳-氮(C-N)之間存在極顯著(P<0.001)的線性關(guān)系,其中3種土地利用方式的土壤C-N關(guān)系的斜率差異不顯著(P=0.458,共同斜率為11.1),但截距差異顯著(P<0.001)。結(jié)合本地區(qū)和全球文獻數(shù)據(jù)分析指出,森林土壤碳氮關(guān)系既有大尺度上的普適性,又有小尺度上對土地利用方式響應(yīng)的局域分異性。

    土地利用變化; 生物地球化學(xué); 化學(xué)計量; 碳氮耦合

    人為活動造成土地利用方式的改變勢必會打破陸地生態(tài)系統(tǒng)原有的碳氮磷平衡[1-3]。土壤翻動不但破壞了土壤團聚體,而且也加速了土壤中動物植物微生物殘體、排泄物、分泌物以及土壤腐殖質(zhì)的破碎速率,使其更易被微生物分解、化學(xué)氧化、物理釋放進入大氣和水域生態(tài)系統(tǒng)[4-5],從而使土壤成為大氣碳庫的重要來源[6]。與此同時,土地利用變化也會造成氮磷元素的丟失[7-9],但丟失程度取決于干擾的強度[1, 10]。例如,Guo等[1]綜合分析全球數(shù)據(jù)表明:天然林轉(zhuǎn)變?yōu)槿斯ち趾娃r(nóng)田過程中土壤碳分別減少了13%和42%,而農(nóng)田轉(zhuǎn)變?yōu)樘烊淮紊趾腿斯ち謺r土壤碳分別增加了18%和53%;Xu等[10]研究表明,在天然林到第一代人工林再到第二代人工林轉(zhuǎn)變過程中土壤碳氮磷含量均顯著減小。

    生態(tài)化學(xué)計量受國內(nèi)外學(xué)者的廣泛關(guān)注,特別是研究碳、氮、磷化學(xué)計量特征可以反映土壤內(nèi)部碳氮磷循環(huán),具有重要的生態(tài)指示作用,因而對認識生態(tài)系統(tǒng)碳匯潛力及其對氣候變化的響應(yīng)具有重要意義[11-15]。Tian[16]、Cleveland等[17]分別從區(qū)域尺度和全球尺度上探討土壤C/N、N/P和C/P化學(xué)計量比時得出:雖然在大尺度上不同生態(tài)系統(tǒng)之間土壤碳、氮、磷含量具有顯著的變異性,但C/N的空間變異較?。欢捎诹自睾吞嫉貋碓吹牟町愋?,土壤C/P和N/P隨氣候、土壤類型、土壤深度和風(fēng)化階段而呈現(xiàn)顯著的空間變異性。而且,天然生態(tài)系統(tǒng)土壤C/N在不同層次間基本保持恒定[16, 18]。Yang等[18]研究也發(fā)現(xiàn),天然高寒草原土壤碳-氮之間具有顯著的共斜率線性關(guān)系,且這種共斜率關(guān)系也存在于不同天然生態(tài)系統(tǒng)中[17]及其不同演替階段[19]。

    東北地區(qū)歷經(jīng)半個多世紀的采伐、墾荒等人為干擾,林地面積顯著減小,原始林基本殆盡,森林年齡結(jié)構(gòu)趨于單一化和低齡化。隨著退耕還林、天然林保護等林業(yè)生態(tài)工程的實施,出現(xiàn)大面積的天然次生林和人工林[20-21]。然而,這種土地利用方式轉(zhuǎn)變對土壤碳氮磷含量、生態(tài)化學(xué)計量特征及其相互關(guān)系的影響還不清楚。本研究以帽兒山生態(tài)站的相鄰地段、相同林齡(均為10年生)、原始植被一致但土地利用方式不同的3種溫帶幼齡林為對象,最大程度地消除土壤本身異質(zhì)性,以便探索土壤碳、氮、磷含量及其生態(tài)化學(xué)計量關(guān)系對土地利用變化的響應(yīng)。

    1 研究方法

    1.1 研究地概況與試驗設(shè)計

    研究地設(shè)在帽兒山森林生態(tài)站(45°20′N,127°30′E)。平均海拔400 m,平均坡度10°—15°, 地帶性土壤為暗棕色森林土。氣候?qū)俅箨懶约撅L(fēng)氣候,夏季溫暖濕潤,冬季寒冷干燥,平均年降水量約629 mm,50%的降水集中在6—8月,平均年蒸發(fā)量約864 mm,年均氣溫3.1 ℃,1月份平均氣溫-18.5 ℃,7月份平均氣溫 22.0 ℃。無霜期約為120—140 d[22]?,F(xiàn)有植被是原地帶性植被闊葉-紅松林經(jīng)干擾后形成的天然次生林和人工林。

    本研究選擇帽兒山生態(tài)站相鄰地段、相同林齡(土壤擾動后自然更新10a,林分現(xiàn)已郁閉)但土地利用方式不同的3種溫帶森林類型,分別為:無土壤翻動的天然次生林 [NS]、間作大豆而土壤翻動中等的人工林[MS]、間作人參而土壤翻動嚴重的次生林[SS]。NS為皆伐后天然更新次生林,土壤無任何擾動,其中優(yōu)勢樹種為白樺(BetulaplatyphyllaSuk.),伴生有春榆(Ulmusjaponica(Rehc.) Sarg.)、水曲柳(FraxinusmandshuricaRupr.)、大青楊(PopulousussuriensisKom.)、黃菠蘿(PhellodendronamurenseRupr.)和五角槭(AcermonoMaxim.)。MS為皆伐后人工栽植水曲柳(株行距1.5 m × 1.5 m),同時翻土深度15 cm左右,間作大豆2a之后混入天然更新樹種白樺、春榆等。SS為皆伐后翻土深度40 cm左右,作壟種植人參10a后天然更新的次生林,其樹種組成與NS相似。于2013年9月在每林型中隨機設(shè)立3塊20 m × 30 m的固定樣地,并對樣地進行植被調(diào)查。

    1.2 樣品采集與分析

    2013年10月,在每塊樣地內(nèi)隨機選取5個樣點,采用內(nèi)徑為50 mm的土鉆分層取樣,取樣深度分別為:0—10 cm、10—20 cm、20—30 cm、30—40 cm、40—50 cm、50—60 cm。土壤樣品風(fēng)干后全部過2 mm篩,除去所有根系和石礫,將土樣粉碎后作為分析土樣。稱取65 ℃烘干粉碎土樣約0.100 g,采用multi N/C 3000分析儀和HT 1500 Solids Module分析儀(Analytik Jena AG, Germany)燃燒法測定土壤全碳含量。另稱取烘干的粉碎土樣約0.200 g,加入5 mL濃硫酸,3 mL過氧化氫420 ℃消煮90 min(消煮完全至白色),待冷卻后定容至100 mL,靜置12 h后吸取上清液利用連續(xù)流動分析儀(BRAN+LUEBBE-AA3, Germany)測定氮、磷含量。

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

    以樣地為統(tǒng)計分析試驗單位,采用單因素方差分析(one-way ANOVA)和最小顯著差異法(LSD)比較不同林型、土壤層次的碳、氮、磷含量及其生態(tài)化學(xué)計量特征;回歸分析法建立不同土地利用方式下土壤碳、氮、磷含量之間的最優(yōu)耦合模型(以每樣方各土層平均土壤碳、氮、磷含量為數(shù)據(jù)分析單位),并采用SMATR檢驗耦合模型的斜率同質(zhì)性[23]。直線回歸系數(shù)差異性檢驗由SMATR Version 2.0軟件完成,其它統(tǒng)計分析均由SPSS19.0軟件完成。

    2 結(jié)果

    2.1 不同土地利用方式下土壤碳、氮、磷含量的垂直分布

    圖1 不同土地利用類型下土壤碳、氮、磷含量的垂直分布Fig.1 Vertical distribution of soil carbon (C), nitrogen (N) and phosphorus (P) contents for different land use types誤差線為標準誤(n=3);誤差線上的不同字母表示同一土壤層次不同土地利用間顯著性差異分組(α=0.05);NS: 無土壤翻動的天然次生林Naturally regenerated stands; MS: 間作大豆而土壤翻動中等的人工林Moderately disturbed stands;SS: 間作人參而土壤翻動嚴重的次生林Severely disturbed stands

    土地利用方式顯著影響土壤碳含量(圖1)。各層次土壤碳含量均呈現(xiàn)NS最高、MS居中、SS最低。NS各土壤層次碳含量平均比MS高19.4%(變化范圍為7.4%—32.4%),比SS高77.4%(變化范圍為69.5%—83.1%)。其中,除0—10 cm外,其它各層次碳含量在NS和MS之間差異不顯著(P>0.05),但兩者均顯著高于SS(P<0.05)。NS(F5,12=19.52,P<0.001)、MS(F5,12=38.71,P<0.001)和SS(F5,12=31.55,P<0.001)土壤碳含量均隨土層加深而顯著降低(圖1A),但同一林型20 cm以下不同土層之間碳含量均無顯著性差異(P>0.05)。

    與碳含量相同,各土壤層次氮含量均呈現(xiàn)NS>MS>SS,且具有相同垂直分布格局(圖1)。但是,除了0—10 cm(F2,6=7.18,P=0.026)、10—20 cm(F2,6=5.77,P=0.040)外,其它各土層氮含量在3種土地利用類型之間差異均不顯著(P>0.05)。NS各土壤層次氮含量平均比MS高38.0%(變化范圍為29.5% — 47.3%),比SS高56.8%(變化范圍為51.8% — 60.8%)。

    然而,與碳、氮含量不同,不同土地利用類型之間各土層的磷含量均無顯著差異(P>0.05)。NS和MS土壤磷含量隨著土層的加深而遞減,而SS的土壤磷含量呈現(xiàn)先遞減、20 cm以下保持基本穩(wěn)定的趨勢(圖1)。

    2.2 土地利用變化對土壤碳、氮、磷生態(tài)化學(xué)計量特征的影響

    土地利用變化顯著影響土壤C/N, C/P和N/P(圖2)。其中,除0—10 cm外,其它各層C/N在不同土地利用方式之間具有顯著差異(P<0.05),并呈現(xiàn)MS>NS>SS。NS不同土層之間C/N基本保持恒定(F5,12=0.923,P=0.499),平均值為15.5;而MS的C/N隨著土壤深度的增加而顯著增加(F5,12=7.64,P=0.002),變化范圍為13.4—28.6;SS的C/N隨著土壤深度的增加而減小,但不同土層之間無顯著差異(F5,12=1.744,P=0.199),平均C/N為7.7。NS和MS各土層C/P均顯著高于SS(P<0.05),而NS和MS之間僅在0—10 cm層具有顯著差異(P=0. 006)。土壤翻動顯著降低了土壤N/P,不同土地利用方式下N/P呈現(xiàn):NS>MS>SS,且隨土層加深而顯著減小(P<0.05)。

    圖2 不同土地利用類型下土壤C/N、C/P和N/P的垂直格局Fig.2 Vertical patterns of C/N,C/P and N/P for different land use types誤差線為標準誤(n=3);誤差線上的不同字母表示同一土壤層次不同土地利用間顯著性差異分組(α=0.05)

    土壤碳氮磷含量的耦合模型因元素和土地利用方式而異(圖3)。碳-氮(C-N)之間具有極顯著(P<0.001)的線性關(guān)系。進一步的SMATR斜率同質(zhì)性檢驗表明,不同土地利用方式之間C-N化學(xué)計量關(guān)系斜率無顯著性差異(t=2.048,P=0.458),共同斜率為11.1,NS和MS之間截距無顯著差異(P>0.05),但兩者均顯著大于SS(P<0.001)(表1)。土壤碳磷(C-P)、氮磷(N-P)相互關(guān)系因土地利用方式而異,NS和MS呈現(xiàn)極顯著(P<0.001)的指數(shù)關(guān)系,而SS卻為直線關(guān)系(圖3)。

    3 討論

    3.1 土地利用方式對土壤碳、氮、磷含量及其垂直分布的影響

    本研究中相鄰地段、相同林齡的不同土地利用方式下3種溫帶幼齡林,因土壤翻動歷史、程度和類型不同,其土壤碳、氮含量產(chǎn)生了顯著的差異。無論是表層還是深層,土壤翻動均促使土壤碳、氮含量的丟失(圖1)。其中,MS平均碳丟失率為19.4%,與Guo等[1]報道的全球平均值(13%)相當。這種中等程度土壤擾動導(dǎo)致的碳丟失主要發(fā)生在土壤表層(0—10 cm),而對深層土壤的碳含量影響不顯著(圖1),這是因為間作大豆2a過程中只翻動了0—15 cm,深層土壤結(jié)構(gòu)沒有受到破壞,但對表層土壤的翻動過程中無疑增加了土壤的通透性,一定程度上促進了深層土壤碳的釋放,從而呈現(xiàn)出輕微的碳丟失,這與以往的研究相符[24-25]。然而,翻動強度大(作壟深達40 cm左右)、時間長(10a)的SS平均碳丟失率高達77.4%,遠高于Guo等[1]報道的全球森林向農(nóng)田轉(zhuǎn)變過程中42%的碳丟失率,也高于Ellert等[8]報道的長期耕作導(dǎo)致34%的碳丟失率。這種重度土壤翻動顯著促使了表層和深層土壤的碳丟失,這是因為種植人參過程中作壟可達40 cm處,一方面土壤結(jié)構(gòu)遭到嚴重破壞,另一方面也加大了深層土壤和外界環(huán)境的接觸機會與面積,從而促進了微生物對土壤有機質(zhì)的分解導(dǎo)致的碳釋放以及隨地表徑流的丟失。此外,Don等[26]整合分析也表明土地利用變化過程中土地利用變化的類型、土壤水熱狀況以及土壤屬性(如土壤類型,粘土含量等)只解釋了55%的土壤碳變化,而將近一半的變異很可能是由土壤翻動強度的差異性所致。

    圖3 不同土地利用類型下土壤碳、氮、磷相互關(guān)系Fig.3 Relationships between soil carbon (C), nitrogen (N) and phosphorus (P) contents for different land use types誤差線為標準誤(n=18)

    表1 不同土地利用方式下土壤碳-氮關(guān)系斜率、截距同質(zhì)性檢驗

    將本研究結(jié)果與本地區(qū)的頂極群落和成熟林相比:NS的碳含量相對于天然次生成熟林(蒙古櫟林、硬闊葉林、楊樺林[27])的表層土壤(0—10 cm)碳丟失率為6.6%;而3種天然次生成熟林表層土壤平均碳含量比本地區(qū)闊葉紅松原始林[3]低47.9%。與碳含量相同,各土壤層次氮含量也均呈現(xiàn)NS>MS>SS(圖1)。NS各土壤層次氮含量平均比MS高38.0%,比SS高56.8%,可見土壤氮含量對土地利用方式的響應(yīng)與碳含量相似。土壤翻動造成碳氮含量的同時丟失已得到眾多學(xué)者的證實[7-9],這是因為土壤有機質(zhì)在分解過程中每釋放25份碳素就會消耗1份氮素[28]。綜上所述,在溫帶森林從原始林-天然次生林-人工林-農(nóng)田轉(zhuǎn)變過程中,土壤碳氮含量不斷降低。因此,采用適當?shù)耐烁€林和森林更新措施可以逆轉(zhuǎn)土壤的碳氮丟失,從而充分發(fā)揮森林土壤的碳氮固持潛能[26, 29]。

    土壤碳和氮主要來源于植物地上地下凋落物的分解,其含量受植被、氣候、人為活動等影響較大[28]。因此,土壤對土壤碳氮的垂直分布的影響會因土地利用而發(fā)生變化,但其主導(dǎo)能力隨著土壤深度的增加而降低;而深層土壤碳氮含量主要決定于土壤本身質(zhì)地,如粘土含量、土壤類型等[18, 30-32]。然而,長期的重度土壤翻動也會導(dǎo)致深層碳氮含量減小(圖1)。這表明量化在退耕還林過程中長期耕作農(nóng)田的深層土壤碳氮匯功能也是必要的[26, 33]。與土壤碳氮不同,土壤磷主要來源于巖石風(fēng)化,其含量大小更大程度上取決于立地的地質(zhì)特征,而受人為影響較小[3, 34]。本研究中0—30 cm土層磷含量呈現(xiàn)NS>MS>SS,這可能是土壤翻動促使土壤有機磷的分解,并隨地表徑流而流失的緣故[32]。

    3.2 土壤碳、氮、磷生態(tài)化學(xué)計量特征對土地利用變化的響應(yīng)

    Sistla等[35]指出,從局部到整體研究生態(tài)化學(xué)計量的可塑性和適應(yīng)性是理解生態(tài)系統(tǒng)碳和養(yǎng)分循環(huán)對全球變化響應(yīng)的有效途徑。與Tian[16]、Cleveland等[17]大尺度研究結(jié)果不同,本研究中NS和MS土壤C/N隨土層加深而呈現(xiàn)出增加的趨勢。Hamer等[36]也報道了C/N的這種垂直分布格局。進一步分析劉爽等[27]報道的當?shù)?種溫帶成熟林結(jié)果也顯示,10—20 cm土壤C/N(20.1)高于0—10 cm(16.7),而C/N與土壤有機質(zhì)分解速率成反比[14],這說明本地區(qū)深層土壤具有較高的碳固存潛力[37],且因本地區(qū)深層土壤溫度較低、土壤微生物分解速率較慢、粘粒含量較高而相對穩(wěn)定。NS和MS之間的差異性可能是因為間作大豆2a中對土壤的翻動促進了表層有機質(zhì)向深層的輸入而補充了翻動導(dǎo)致的有機質(zhì)輸出,而在這些動態(tài)變化過程中碳對外界干擾較氮更敏感[38];另一方面本研究中NS和MS具有不同的植被組成,而植被特征一定程度上會影響土壤C/N的垂直分布。例如,最近王建林等[39]發(fā)現(xiàn),不同高寒草地類型下土壤C/N具有5種不同的剖面分布型。此外,也可能是由于深層土壤碳氮的小尺度空間異質(zhì)性所導(dǎo)致[40]。種植人參過程中的重度土壤翻動導(dǎo)致了C/N的顯著減小,經(jīng)過10a的恢復(fù),SS僅表層土壤C/N和NS與MS無顯著差異,這說明盡管本地區(qū)森林深層土壤具有穩(wěn)定而又較高的碳含量[37],但是在土地利用變化過程中一旦發(fā)生碳丟失,便需要相當長的時間才能恢復(fù)。

    同樣,土地利用變化也顯著影響土壤C/P和N/P。3種不同土地利用方式下的C/P的剖面分布與C/N類似。相對于C/N、C/P,N/P垂直分布,與Tian等[16]研究結(jié)果相符,且不同土壤層次之間均呈現(xiàn)NS>MS>SS。這說明在消除了土壤本身異質(zhì)性的小尺度上,土壤N/P可能比C/N和C/P與土地利用變化具有更高的同步性(隨土壤翻動強度的增加而顯著減小,圖2)。林麗等[13]也發(fā)現(xiàn),土壤N/P對草地退化演替的敏感性較高。因此,N/P比C/N和C/P對外界環(huán)境變化可能具有更優(yōu)的生態(tài)指示功能。表層土壤C/N、C/P和N/P均呈現(xiàn)SS < MS < NS,說明在森林恢復(fù)過程中會受到氮和磷限制。然而,盡管Koerselman[41]、Tessier[42]等先后通過綜述植被N/P來分析其對養(yǎng)分限制的指示作用,但由于受土壤本身碳氮磷絕對含量的影響,目前有關(guān)土壤碳氮磷化學(xué)計量對土壤養(yǎng)分限制的指示作用還無定論[14]。

    本研究結(jié)果表明,3種林型土壤碳-氮-磷之間存在極顯著的耦聯(lián)關(guān)系(圖3)。Cleveland等[17]在綜述陸地生態(tài)系統(tǒng)土壤和微生物的“Redfield ratio[43]”時、Yang等[19]在綜述全球森林演替過程中土壤碳氮動態(tài)變化以及Tian等[16]在綜合我國土壤碳氮磷化學(xué)計量特征時均發(fā)現(xiàn)土壤碳氮之間存在著顯著的耦合關(guān)系。進一步的回歸系數(shù)差異性分析顯示,土壤碳-氮之間均存在顯著的共斜率線性關(guān)系,但土地利用變化顯著影響其截距(P<0.001,表2)。這表明土地利用變化并沒有顯著改變土壤碳-氮關(guān)系,而是通過改變其直線關(guān)系的截距來影響土壤碳氮的丟失以及C/N,這與最近研究相符[34, 44-45]。

    圖4 普適性土壤碳-氮關(guān)系 Fig.4 Generalized relationships between soil carbon (C) and nitrogen (N)NS、MS、SS分別表示本研究中不同土地利用方式幼齡林;MF表示本地區(qū)的成熟林(mature forest)土壤碳-氮關(guān)系,土壤碳氮數(shù)據(jù)來自劉爽等[27];Global表示全球森林生態(tài)系統(tǒng)土壤碳-氮關(guān)系,數(shù)據(jù)來源于Cleveland等[17];方框中數(shù)據(jù)表示異常值

    為了進一步驗證森林土壤碳氮化學(xué)計量關(guān)系的普適性,將本研究中的幼齡林和同一地區(qū)的5種典型森林類型的成熟林(即紅松人工林、落葉松人工林、硬闊葉林、蒙古櫟和楊樺林)[27]以及Cleveland等[17]的全球不同森林生態(tài)系統(tǒng)土壤碳氮數(shù)據(jù)作綜合分析。結(jié)果發(fā)現(xiàn),不同尺度的土壤碳-氮之間線性關(guān)系的斜率均無顯著性差異(t=2.336,P=0.749),其共同斜率為11.2(95%置信區(qū)間為[10.3, 11.8],P<0.001,圖4)。全球尺度上碳-氮之間線性關(guān)系的截距與0無顯著性差異(95%置信區(qū)間為[-4.5, 15.9])。這些結(jié)果反映了土壤C/N在大尺度上恒定性,這種森林土壤碳氮關(guān)系的普適性在陸地生態(tài)系統(tǒng)生物地球化學(xué)循環(huán)研究中具有重要意義[18]。然而,大尺度上的異質(zhì)性又掩蓋了C-N關(guān)系截距的差異性,這從另一方面說明小尺度上研究碳-氮之間線性關(guān)系截距的差異性也具有重要的生態(tài)指示作用。

    與碳-氮關(guān)系不同,在最大程度消除土壤本身異質(zhì)性的情況下,土壤碳-磷、氮-磷之間的相互關(guān)系因土地利用方式而異。NS和MS呈現(xiàn)極顯著的指數(shù)關(guān)系,而SS卻為直線關(guān)系(圖3)。Chen等[3]發(fā)現(xiàn)土壤磷含量對土地利用變化的響應(yīng)與碳和氮并不相同,土壤全磷含量變化的不確定性同樣在落演替過程中發(fā)現(xiàn)[13]。這說明磷元素和碳氮元素不僅在來源上不同,且對外界干擾的響應(yīng)機制可能也不同[16],從而導(dǎo)致碳-磷、氮-磷關(guān)系的局域性。

    4 結(jié)論

    (1)土壤翻動會導(dǎo)致整個土壤剖面碳含量的顯著減小,而氮含量僅在0—20 cm顯著變化,而土壤磷含量保持恒定,說明土壤碳氮磷含量對土地利用變化響應(yīng)機制可能不同。NS各土壤層次碳含量平均比MS和SS分別高19.4%和77.4%,氮含量平均比MS和SS分別高38%和56.8%,表明土壤翻動強度越大、持續(xù)時間越長,造成土壤碳氮丟失越嚴重。

    (2)相比與NS和MS,各土層SS的C/N、C/P均最低,而NS和MS之間大小因土壤層次而異。3種土地利用方式下C/N和C/P具有不同的垂直分布。與C/N和C/P不同,各土層N/P與土地利用變化同步(NS>MS>SS),且均為從高到低剖面分布。表明土壤碳氮磷學(xué)計量比具有可塑性與適應(yīng)性,氮磷化學(xué)計量比作為土地利用變化的敏感指標優(yōu)于碳氮和碳磷計量比。

    (3)不同土地利用類型的C-P和N-P具有不同的耦合關(guān)系。其中C-N關(guān)系的斜率無顯著性差異,而其截距的差異顯著。這反映了森林土壤碳氮關(guān)系在大尺度上的普適性和小尺度上對土地利用方式的局域差異性。

    致謝:帽兒山森林生態(tài)站提供野外基礎(chǔ)支持。

    [1] Guo L B, Gifford R M. Soil carbon stocks and land use change:a meta analysis. Global Change Biology, 2002, 8(4):345-360.

    [2] Post W M, Kwon K C. Soil carbon sequestration and land-use change:processes and potential. Global Change Biology, 2000, 6(3):317-327.

    [3] Chen X W, Li B L. Change in soil carbon and nutrient storage after human disturbance of a primary Korean pine forest in Northeast China. Forest Ecology and Management, 2003, 186(1/3):197-206.

    [4] Baker J M, Ochsner T E, Venterea R T, Griffis T J. Tillage and soil carbon sequestration—What do we really know? Agriculture, Ecosystems & Environment, 2007, 118(1/4):1-5.

    [5] Cambardella C A, Elliott E T. Carbon and nitrogen distribution in aggregates from cultivated and native grassland soils. Soil Science Society of America Journal, 1993, 57(4):1071-1076.

    [6] Lal R. World cropland soils as a source or sink for atmospheric carbon. Advances in Agronomy, 2001, 71:145-191.

    [7] Ding F, Hu Y L, Li L J, Li A, Shi S W, Lian P Y, Zeng D H. Changes in soil organic carbon and total nitrogen stocks after conversion of meadow to cropland in Northeast China. Plant and Soil, 2013, 373(1/2):659-672.

    [8] Ellert B H, Gregorich E G. Storage of carbon, nitrogen and phosphorus in cultivated and adjacent forested soils of Ontario. Soil Science, 1996, 161(9):587-603.

    [9] Gál A, Vyn T J, Michéli E, Kladivko E J, McFee W W. Soil carbon and nitrogen accumulation with long-term no-till versus moldboard plowing overestimated with tilled-zone sampling depths. Soil and Tillage Research, 2007, 96(1/2):42-51.

    [10] Xu Z H, Ward S, Chen C R, Blumfield T, Prasolova N, Liu J X. Soil carbon and nutrient pools, microbial properties and gross nitrogen transformations in adjacent natural forest and hoop pine plantations of subtropical Australia. Journal of Soils and Sediments, 2008, 8(2):99-105.

    [11] 程濱, 趙永軍, 張文廣, 安樹青. 生態(tài)化學(xué)計量學(xué)研究進展. 生態(tài)學(xué)報, 2010, 30(6):1628-1637.

    [12] 李從娟, 雷加強, 徐新文, 唐清亮, 高培, 王永東. 塔克拉瑪干沙漠腹地人工植被及土壤CNP的化學(xué)計量特征. 生態(tài)學(xué)報, 2013, 33(18):5760-5767.

    [13] 林麗, 李以康, 張法偉, 杜巖功, 郭小偉, 李婧, 劉淑麗, 曹廣民. 高寒矮嵩草群落退化演替系列氮、磷生態(tài)化學(xué)計量學(xué)特征. 生態(tài)學(xué)報, 2013, 33(17):5245-5251.

    [14] 王紹強, 于貴瑞. 生態(tài)系統(tǒng)碳氮磷元素的生態(tài)化學(xué)計量學(xué)特征. 生態(tài)學(xué)報, 2008, 28(8):3937-3947.

    [15] 王維奇, 徐玲琳, 曾從盛, 仝川, 張林海. 河口濕地植物活體-枯落物-土壤的碳氮磷生態(tài)化學(xué)計量特征. 生態(tài)學(xué)報, 2012, 31(23):7119-7124.

    [16] Tian H Q, Chen G S, Zhang C, Melillo J M, Hall C A S. Pattern and variation of C:N:P ratios in China′s soils:a synthesis of observational data. Biogeochemistry, 2010, 98(1/3):139-151.

    [17] Cleveland C C, Liptzin D. C:N:P stoichiometry in soil:is there a “Redfield ratio” for the microbial biomass? Biogeochemistry, 2007, 85(3):235-252.

    [18] Yang Y H, Fang J Y, Guo D L, Ji C J, Ma W H. Vertical patterns of soil carbon, nitrogen and carbon:nitrogen stoichiometry in Tibetan grasslands. Biogeosciences Discussions, 2010, 7(1):1-24.

    [19] Yang Y H, Luo Y Q. Carbon:nitrogen stoichiometry in forest ecosystems during stand development. Global Ecology and Biogeography, 2011, 20(2):354-361.

    [20] Chen X W, Li B L, Lin Z S. The acceleration of succession for the restoration of the mixed-broadleaved Korean pine forests in Northeast China. Forest Ecology and Management, 2003, 177(1/3):503-514.

    [21] Deng X Z, Jiang Q O, Zhan J Y, He S J, Lin Y Z. Simulation on the dynamics of forest area changes in Northeast China. Journal of Geographical Sciences, 2010, 20(4):495-509.

    [22] Wang C K, Han Y, Chen J Q, Wang X C, Zhang Q Z, Bond-Lamberty B. Seasonality of soil CO2efflux in a temperate forest:Biophysical effects of snowpack and spring freeze-thaw cycles. Agricultural and Forest Meteorology, 2013, 177:83-92.

    [23] Warton D I, Wright I J, Falster D S, Westoby M. Bivariate line-fitting methods for allometry. Biological Reviews, 2006, 81(2):259-291.

    [24] Mishra U, Ussiri D A N, Lal R. Tillage effects on soil organic carbon storage and dynamics in Corn Belt of Ohio USA. Soil and Tillage Research, 2010, 107(2):88-96.

    [25] Tiessen H, Stewart J W B, Bettany J R. Cultivation effects on the amounts and concentration of carbon, nitrogen, and phosphorus in grassland soils. Agronomy Journal, 1982, 74(5):831-835.

    [26] Don A, Schumacher J, Freibauer A. Impact of tropical land-use change on soil organic carbon stocks-a meta-analysis. Global Change Biology, 2011, 17(4):1658-1670.

    [27] 劉爽, 王傳寬. 五種溫帶森林土壤微生物生物量碳氮的時空格局. 生態(tài)學(xué)報, 2010, 30(12):3135-3143.

    [28] Chapin III F S, Matson P A, Vitousek P M. Principles of Terrestrial Ecosystem Ecology. 2nd ed. New York:Springer, 2011.

    [29] Meyfroidt P, Rudel T K, Lambin E F. Forest transitions, trade, and the global displacement of land use. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(49):20917-20922.

    [30] Jobbágy E G, Jackson R B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecological Applications, 2000, 10(2):423-436.

    [31] Jobbágy E G, Jackson R B. The distribution of soil nutrients with depth:global patterns and the imprint of plants. Biogeochemistry, 2001, 53(1):51-77.

    [32] Meersmans J, Van Wesemael B, De Ridder F, Fallas Dotti M, De Baets S, Van Molle M. Changes in organic carbon distribution with depth in agricultural soils in northern Belgium, 1960-2006. Global Change Biology, 2009, 15(11):2739-2750.

    [33] Rumpel C, K?gel-Knabner I. Deep soil organic matter—a key but poorly understood component of terrestrial C cycle. Plant and Soil, 2011, 338(1/2):143-158.

    [34] Wang W J, Qiu L, Zu Y G, Su D X, An J, Wang H Y, Zheng G Y, Sun W, Chen X Q. Changes in soil organic carbon, nitrogen, pH and bulk density with the development of larch (Larixgmelinii) plantations in China. Global Change Biology, 2011, 17(8):2657-2676.

    [35] Sistla S A, Schimel J P. Stoichiometric flexibility as a regulator of carbon and nutrient cycling in terrestrial ecosystems under change. New Phytologist, 2012, 196(1):68-78.

    [36] Hamer U, Potthast K, Burneo J I, Makeschin F. Nutrient stocks and phosphorus fractions in mountain soils of Southern Ecuador after conversion of forest to pasture. Biogeochemistry, 2013, 112(1/3):495-510.

    [37] 楊金艷, 王傳寬. 東北東部森林生態(tài)系統(tǒng)土壤碳貯量和碳通量. 生態(tài)學(xué)報, 2005, 25(11):2875-2882.

    [38] 羅亞勇, 張宇, 張靜輝, 卡召加, 尚倫宇, 王少影. 不同退化階段高寒草甸土壤化學(xué)計量特征. 生態(tài)學(xué)雜志, 2012, 31(2):254-260.

    [39] 王建林, 鐘志明, 王忠紅, 陳寶雄, 余成群, 胡興祥, 沈振西, 大次卓嘎, 張憲洲. 青藏高原高寒草原生態(tài)系統(tǒng)土壤氮磷比的分布特征. 應(yīng)用生態(tài)學(xué)報, 2013, 24(12):3399-3406.

    [40] 武小鋼, 郭晉平, 田旭平, 楊秀云. 蘆芽山亞高山草甸、云杉林土壤有機碳、全氮含量的小尺度空間異質(zhì)性. 生態(tài)學(xué)報, 2013, 33(24):7756-7764.

    [41] Koerselman W, Meuleman A F M. The vegetation N:P ratio:a new tool to detect the nature of nutrient limitation. Journal of Applied Ecology, 1996, 33(6):1441-1450.

    [42] Tessier J T, Raynal D J. Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. Journal of Applied Ecology, 2003, 40(3):523-534.

    [43] Redfield A C. The biological control of chemical factors in the environment. American Scientist, 1958, 46(3):205-221.

    [44] Li H, Li J, He Y L, Li S J, Liang Z S, Peng C H, Polle A, Luo Z B. Changes in carbon, nutrients and stoichiometric relations under different soil depths, plant tissues and ages in black locust plantations. Acta Physiologiae Plantarum, 2013, 35(10):2951-2964.

    [45] Tischer A, Potthast K, Hamer U. Land-use and soil depth affect resource and microbial stoichiometry in a tropical mountain rainforest region of southern Ecuador. Oecologia, 2014, 175(1):375-393.

    The effect of land use change on soil carbon, nitrogen, and phosphorus contents and their stoichiometry in temperate sapling stands in northeastern China

    ZHOU Zhenghu, WANG Chuankuan*, ZHANG Quanzhi

    CenterforEcologicalResearch,NortheastForestryUniversity,Harbin150040,China

    Anthropogenic land use change alters the equilibrium among soil carbon (C), nitrogen (N), and phosphorus (P) in the forest ecosystem, likely affecting the biogeochemical cycles. In this study, conducted in northeastern China, we measured C, N, and P contents along soil profiles in adjacent sapling stands identical in age (10 years) and primary vegetation, but differing in land use type. Land use types included naturally regenerated stands without soil disturbance (NS), moderately disturbed stands (two-year interplanting of soybean after artificial reforestation, MS), and severely disturbed stands (natural reforestation after plantingPanaxginsengfor 10 years, SS). Our objective was to examine the effects of land use change on soil C, N, and P contents and their stoichiometry. Our results showed that:(1) Land use change significantly altered both top-and sub-soil C contents, which descended in the order NS>MS>SS for each soil horizon. However, land use markedly influenced only top-soil N contents and did not significantly influence P content along the entire soil profile (P> 0.05). These results suggest differing sensitivities of elements to land use change. (2) The C/N and C/P ratios among the three land use types were lowest for SS and differed between NS and MS depending on the soil horizon. Vertical patterns of C/N and C/P ratios differed among land use types, but N/P ratio was consistently constrained by land use change. N/P ratios descended in the order NS>MS>SS within each soil horizon and decreased with soil depth (P< 0.05), suggesting that N/P may be a better ecological indicator of response to land use change than C/N and C/P. (3) Land use change significantly affected coupled relationships between C-N, C-P, and N-P. Soil C and N contents were linearly correlated (P< 0.001), regardless of the land use type, with a common slope of 11.1 but significantly different intercepts (P< 0.001). This result, together with the synthesis of published local and global datasets, reflects a large-scale general C-N relationship and a small-scale heterogeneous response of C-N stoichiometry to local land use change.

    land use change; biogeochemistry; stoichiometry; carbon-nitrogen coupling

    林業(yè)公益性行業(yè)專項(201104009-05); 國家“十二五”科技支撐項目(2011BAD37B01); 教育部長江學(xué)者和創(chuàng)新團隊發(fā)展計劃(IRT1054)

    2014-03-29; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2014-12-18

    10.5846/stxb201403290589

    *通訊作者Corresponding author.E-mail: wangck-cf@nefu.edu.cn

    周正虎,王傳寬,張全智.土地利用變化對東北溫帶幼齡林土壤碳氮磷含量及其化學(xué)計量特征的影響.生態(tài)學(xué)報,2015,35(20):6694-6702.

    Zhou Z H, Wang C K, Zhang Q Z.The effect of land use change on soil carbon, nitrogen, and phosphorus contents and their stoichiometry in temperate sapling stands in northeastern China.Acta Ecologica Sinica,2015,35(20):6694-6702.

    猜你喜歡
    氮磷土地利用計量
    玻利維亞擬建新的氮磷鉀肥料工廠
    《化學(xué)分析計量》2020年第6期目次
    關(guān)注日常 計量幸福
    特別健康(2018年4期)2018-07-03 00:38:26
    計量自動化在線損異常中的應(yīng)用
    土地利用生態(tài)系統(tǒng)服務(wù)研究進展及啟示
    濱??h土地利用挖潛方向在哪里
    植物需要氮磷鉀
    基于因子分析的人力資本計量研究
    不同氮磷水平下AM真菌對棉花的侵染率及氮磷含量的影響
    三種填料對模擬河水中氮磷吸附性能對比
    日本黄大片高清| 麻豆成人午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 一本综合久久免费| ponron亚洲| 精品国产美女av久久久久小说| 国产伦人伦偷精品视频| 亚洲成a人片在线一区二区| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 国产成人福利小说| x7x7x7水蜜桃| 亚洲无线观看免费| 亚洲国产精品999在线| 国产欧美日韩精品亚洲av| 午夜老司机福利剧场| 国产单亲对白刺激| 国产精品一区二区三区四区免费观看 | 日本a在线网址| 免费在线观看成人毛片| 搡老熟女国产l中国老女人| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| 真人一进一出gif抽搐免费| 国产综合懂色| 老汉色∧v一级毛片| 久久久久久久久久黄片| 欧美日韩一级在线毛片| 免费av观看视频| 大型黄色视频在线免费观看| 国产精品女同一区二区软件 | 久久性视频一级片| 黑人欧美特级aaaaaa片| 悠悠久久av| 亚洲精品色激情综合| 日本一二三区视频观看| 脱女人内裤的视频| 国产精品影院久久| 极品教师在线免费播放| 国产精品久久久久久久久免 | 久久国产精品影院| 97超视频在线观看视频| 精品熟女少妇八av免费久了| 一卡2卡三卡四卡精品乱码亚洲| 可以在线观看毛片的网站| 亚洲 欧美 日韩 在线 免费| 国产成人福利小说| x7x7x7水蜜桃| 久久久久国产精品人妻aⅴ院| 亚洲成人精品中文字幕电影| 久久久久久久精品吃奶| 精品人妻一区二区三区麻豆 | 亚洲欧美一区二区三区黑人| 成人18禁在线播放| 熟女电影av网| 亚洲激情在线av| 亚洲精品久久国产高清桃花| 国产成人a区在线观看| 日韩中文字幕欧美一区二区| 亚洲精品国产精品久久久不卡| 久9热在线精品视频| 一级毛片女人18水好多| 亚洲精品一区av在线观看| 国产一区二区三区在线臀色熟女| 很黄的视频免费| 99久久精品一区二区三区| 亚洲成av人片在线播放无| 每晚都被弄得嗷嗷叫到高潮| 91av网一区二区| 99在线人妻在线中文字幕| 久久天躁狠狠躁夜夜2o2o| 青草久久国产| 男插女下体视频免费在线播放| 真人一进一出gif抽搐免费| 国产单亲对白刺激| 老司机福利观看| 日日摸夜夜添夜夜添小说| 久久久久国产精品人妻aⅴ院| 免费在线观看影片大全网站| 亚洲精品在线美女| 午夜影院日韩av| 精华霜和精华液先用哪个| 黄片小视频在线播放| 女生性感内裤真人,穿戴方法视频| 我要搜黄色片| 最近在线观看免费完整版| 免费在线观看亚洲国产| 亚洲av不卡在线观看| 久久久精品大字幕| 午夜免费成人在线视频| 亚洲av熟女| 亚洲久久久久久中文字幕| 国产黄片美女视频| 五月玫瑰六月丁香| 首页视频小说图片口味搜索| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久大av| 十八禁人妻一区二区| 国产成人a区在线观看| 一进一出抽搐gif免费好疼| 久久久国产精品麻豆| 两人在一起打扑克的视频| 日韩国内少妇激情av| 国产不卡一卡二| 天天躁日日操中文字幕| 欧美最黄视频在线播放免费| 亚洲国产精品sss在线观看| 最近最新中文字幕大全免费视频| 国产真实乱freesex| 99热这里只有精品一区| 国产激情偷乱视频一区二区| 在线十欧美十亚洲十日本专区| 老汉色∧v一级毛片| 国产高清激情床上av| or卡值多少钱| 18禁裸乳无遮挡免费网站照片| 亚洲欧美激情综合另类| 欧美极品一区二区三区四区| av片东京热男人的天堂| 一卡2卡三卡四卡精品乱码亚洲| 夜夜躁狠狠躁天天躁| 可以在线观看毛片的网站| av福利片在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲中文字幕日韩| 色噜噜av男人的天堂激情| 亚洲精品影视一区二区三区av| 中文字幕人妻熟人妻熟丝袜美 | 69av精品久久久久久| 一本精品99久久精品77| 亚洲精品在线观看二区| 两个人视频免费观看高清| 亚洲不卡免费看| 熟女电影av网| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 99视频精品全部免费 在线| 少妇的逼好多水| 97超级碰碰碰精品色视频在线观看| 成人高潮视频无遮挡免费网站| 久久99热这里只有精品18| 国产亚洲精品久久久com| 亚洲av中文字字幕乱码综合| 久久这里只有精品中国| 久久九九热精品免费| 少妇的逼水好多| 国产真实伦视频高清在线观看 | 欧美日韩乱码在线| 天堂网av新在线| av中文乱码字幕在线| 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美 | 特大巨黑吊av在线直播| 在线观看免费视频日本深夜| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久久久久久久| 国产高清视频在线播放一区| 亚洲黑人精品在线| 亚洲av中文字字幕乱码综合| 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| 中文字幕久久专区| 精华霜和精华液先用哪个| 午夜免费激情av| 国产精品野战在线观看| 亚洲精品成人久久久久久| 免费看十八禁软件| 欧美乱码精品一区二区三区| 高清日韩中文字幕在线| 国产麻豆成人av免费视频| 国产探花极品一区二区| 精品午夜福利视频在线观看一区| 国产视频一区二区在线看| 国产成人福利小说| 免费在线观看亚洲国产| 亚洲中文字幕日韩| 日本a在线网址| 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 男女视频在线观看网站免费| 99久久99久久久精品蜜桃| 一区福利在线观看| 女警被强在线播放| 在线天堂最新版资源| 波多野结衣高清作品| 首页视频小说图片口味搜索| 中文亚洲av片在线观看爽| 欧美一区二区国产精品久久精品| 99久久成人亚洲精品观看| 中文字幕av成人在线电影| 69av精品久久久久久| 看免费av毛片| 午夜福利欧美成人| 最近在线观看免费完整版| 内射极品少妇av片p| 在线免费观看的www视频| 国内精品久久久久久久电影| 日韩中文字幕欧美一区二区| 男女之事视频高清在线观看| 午夜视频国产福利| 美女黄网站色视频| 一卡2卡三卡四卡精品乱码亚洲| 人妻夜夜爽99麻豆av| 一夜夜www| 日本一本二区三区精品| avwww免费| 亚洲成人久久爱视频| 少妇高潮的动态图| 舔av片在线| 国产又黄又爽又无遮挡在线| 午夜福利欧美成人| e午夜精品久久久久久久| 日本成人三级电影网站| 69人妻影院| 精品国内亚洲2022精品成人| 色综合亚洲欧美另类图片| 国产麻豆成人av免费视频| 波野结衣二区三区在线 | 真实男女啪啪啪动态图| 美女免费视频网站| 99久久精品国产亚洲精品| 麻豆成人av在线观看| 国产精品精品国产色婷婷| 亚洲av免费在线观看| 午夜激情福利司机影院| 国内精品美女久久久久久| 狠狠狠狠99中文字幕| 精品国内亚洲2022精品成人| 久久精品国产清高在天天线| 色av中文字幕| 精品人妻1区二区| 久久久久九九精品影院| 少妇高潮的动态图| 制服丝袜大香蕉在线| 夜夜躁狠狠躁天天躁| 亚洲精品成人久久久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品合色在线| 日韩欧美免费精品| 午夜视频国产福利| 哪里可以看免费的av片| 制服人妻中文乱码| 热99在线观看视频| 国产成人aa在线观看| 欧美+日韩+精品| 成人三级黄色视频| 性色av乱码一区二区三区2| 久久久久性生活片| 国产亚洲精品av在线| 亚洲美女黄片视频| 露出奶头的视频| x7x7x7水蜜桃| 18+在线观看网站| 精品人妻偷拍中文字幕| 久久久久久大精品| 美女 人体艺术 gogo| 搡老熟女国产l中国老女人| 日韩成人在线观看一区二区三区| 日本a在线网址| 99热这里只有精品一区| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 久久久久久九九精品二区国产| 国产精品av视频在线免费观看| 在线观看午夜福利视频| 久久久久精品国产欧美久久久| 国产老妇女一区| 97超视频在线观看视频| 麻豆国产av国片精品| 人妻丰满熟妇av一区二区三区| 久久精品91蜜桃| 亚洲精品成人久久久久久| 国产精品 国内视频| 高潮久久久久久久久久久不卡| 国产不卡一卡二| 国内精品久久久久久久电影| 久久精品人妻少妇| 免费观看的影片在线观看| 日本熟妇午夜| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久免费视频| 国产精品免费一区二区三区在线| 精品人妻一区二区三区麻豆 | 1024手机看黄色片| 国产一级毛片七仙女欲春2| 久久亚洲真实| 一本综合久久免费| 偷拍熟女少妇极品色| 天堂av国产一区二区熟女人妻| 亚洲av美国av| 黑人欧美特级aaaaaa片| 99热只有精品国产| 欧美中文综合在线视频| 一a级毛片在线观看| 国产亚洲欧美98| 大型黄色视频在线免费观看| 黄色片一级片一级黄色片| 麻豆成人av在线观看| 欧美zozozo另类| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| xxx96com| 色吧在线观看| 日韩欧美国产在线观看| 久久精品国产清高在天天线| 级片在线观看| 国产精品免费一区二区三区在线| 国产又黄又爽又无遮挡在线| 精品乱码久久久久久99久播| 久久人妻av系列| 久久精品夜夜夜夜夜久久蜜豆| www国产在线视频色| 日本五十路高清| 一进一出抽搐动态| 亚洲专区国产一区二区| 国产一区二区在线av高清观看| e午夜精品久久久久久久| 欧美绝顶高潮抽搐喷水| 亚洲av免费在线观看| 成年免费大片在线观看| 午夜两性在线视频| 国产欧美日韩精品亚洲av| 国产中年淑女户外野战色| 国产成+人综合+亚洲专区| 久久久精品大字幕| 亚洲专区中文字幕在线| 欧美一区二区亚洲| 色av中文字幕| 国产高潮美女av| 少妇裸体淫交视频免费看高清| 欧美国产日韩亚洲一区| 一区二区三区高清视频在线| 国产一区二区三区在线臀色熟女| 亚洲第一欧美日韩一区二区三区| 特级一级黄色大片| 久久精品夜夜夜夜夜久久蜜豆| 精品欧美国产一区二区三| 亚洲国产日韩欧美精品在线观看 | 非洲黑人性xxxx精品又粗又长| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 日韩有码中文字幕| 欧美精品啪啪一区二区三区| 欧美黄色片欧美黄色片| 在线免费观看不下载黄p国产 | 在线观看免费午夜福利视频| 老熟妇仑乱视频hdxx| 久9热在线精品视频| 人妻夜夜爽99麻豆av| 一本一本综合久久| 亚洲成人久久性| 麻豆成人午夜福利视频| 大型黄色视频在线免费观看| 少妇人妻精品综合一区二区 | 草草在线视频免费看| 亚洲av电影不卡..在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩亚洲国产一区二区在线观看| 午夜a级毛片| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产| 日本在线视频免费播放| 18禁国产床啪视频网站| 99热这里只有是精品50| 国产欧美日韩精品亚洲av| 啦啦啦观看免费观看视频高清| 亚洲国产欧洲综合997久久,| 欧美最黄视频在线播放免费| 90打野战视频偷拍视频| 久久精品国产综合久久久| 淫秽高清视频在线观看| 特级一级黄色大片| 啪啪无遮挡十八禁网站| 黄色成人免费大全| 亚洲最大成人手机在线| 久久久久久九九精品二区国产| 国产成人aa在线观看| 嫩草影视91久久| 中文字幕精品亚洲无线码一区| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 91av网一区二区| 人妻久久中文字幕网| 午夜a级毛片| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| 久久久久久久午夜电影| 亚洲无线观看免费| 91久久精品国产一区二区成人 | 天美传媒精品一区二区| 网址你懂的国产日韩在线| 亚洲精品久久国产高清桃花| 亚洲中文字幕一区二区三区有码在线看| 真人一进一出gif抽搐免费| 欧美一区二区精品小视频在线| 黄色片一级片一级黄色片| 欧美又色又爽又黄视频| 99久久综合精品五月天人人| 日韩欧美一区二区三区在线观看| 国产又黄又爽又无遮挡在线| 国产国拍精品亚洲av在线观看 | 丰满人妻一区二区三区视频av | 久久久久久久久大av| xxx96com| 久久久国产精品麻豆| 亚洲无线观看免费| 免费搜索国产男女视频| 国产精品亚洲一级av第二区| 国产精品嫩草影院av在线观看 | netflix在线观看网站| 听说在线观看完整版免费高清| 又黄又粗又硬又大视频| 一个人观看的视频www高清免费观看| 在线a可以看的网站| 午夜福利视频1000在线观看| 久久精品人妻少妇| 最新在线观看一区二区三区| 国产不卡一卡二| 在线播放国产精品三级| 亚洲精品美女久久久久99蜜臀| 日日摸夜夜添夜夜添小说| 天天添夜夜摸| 成年版毛片免费区| 一本一本综合久久| av中文乱码字幕在线| 精品午夜福利视频在线观看一区| 高清日韩中文字幕在线| 99久久精品一区二区三区| av在线天堂中文字幕| 日本熟妇午夜| 黄片大片在线免费观看| 久久久久久九九精品二区国产| 中文字幕人妻丝袜一区二区| 人人妻人人看人人澡| 18禁在线播放成人免费| 激情在线观看视频在线高清| 精品久久久久久久末码| 成人av一区二区三区在线看| 成人鲁丝片一二三区免费| 亚洲片人在线观看| 中文在线观看免费www的网站| 天堂动漫精品| 国产亚洲精品综合一区在线观看| 久久6这里有精品| 一级黄色大片毛片| 熟女人妻精品中文字幕| 亚洲国产精品999在线| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 18+在线观看网站| 日本a在线网址| 麻豆成人av在线观看| 国产淫片久久久久久久久 | 91麻豆精品激情在线观看国产| 免费在线观看日本一区| 国产精品女同一区二区软件 | 内射极品少妇av片p| 最近最新免费中文字幕在线| 色哟哟哟哟哟哟| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 99久久成人亚洲精品观看| 国产国拍精品亚洲av在线观看 | 国产精品日韩av在线免费观看| 欧美乱色亚洲激情| 18美女黄网站色大片免费观看| 亚洲成人久久性| 欧美3d第一页| 久久久久久久久大av| 日日摸夜夜添夜夜添小说| 国产伦人伦偷精品视频| 成人特级黄色片久久久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品sss在线观看| 尤物成人国产欧美一区二区三区| 国产三级黄色录像| 免费看光身美女| 亚洲精品久久国产高清桃花| 嫩草影视91久久| 欧美另类亚洲清纯唯美| 少妇的逼好多水| 老司机深夜福利视频在线观看| 最后的刺客免费高清国语| 在线观看av片永久免费下载| 亚洲av电影不卡..在线观看| 国产欧美日韩精品亚洲av| 一本综合久久免费| 桃红色精品国产亚洲av| 亚洲精品一区av在线观看| 日本三级黄在线观看| 亚洲国产欧洲综合997久久,| 国产免费男女视频| 俺也久久电影网| 神马国产精品三级电影在线观看| 日本在线视频免费播放| 丰满乱子伦码专区| 婷婷精品国产亚洲av| av女优亚洲男人天堂| 欧美av亚洲av综合av国产av| 欧美大码av| 中文字幕熟女人妻在线| 在线a可以看的网站| 无人区码免费观看不卡| 久久婷婷人人爽人人干人人爱| 日韩大尺度精品在线看网址| 亚洲欧美一区二区三区黑人| 精品一区二区三区人妻视频| 国产精品久久久久久久久免 | 日本撒尿小便嘘嘘汇集6| 久久久成人免费电影| 波多野结衣巨乳人妻| 夜夜爽天天搞| 亚洲精品亚洲一区二区| 熟女少妇亚洲综合色aaa.| 99在线人妻在线中文字幕| 成人高潮视频无遮挡免费网站| 男人和女人高潮做爰伦理| 在线视频色国产色| 淫妇啪啪啪对白视频| 国产伦人伦偷精品视频| 亚洲一区二区三区色噜噜| 久久草成人影院| 99热只有精品国产| 一个人看视频在线观看www免费 | 久久国产乱子伦精品免费另类| 级片在线观看| 国产精品精品国产色婷婷| 亚洲人成网站在线播放欧美日韩| 夜夜夜夜夜久久久久| 国产精品免费一区二区三区在线| 国产伦在线观看视频一区| 嫁个100分男人电影在线观看| 国产三级在线视频| 99热这里只有是精品50| 国产精品女同一区二区软件 | 久99久视频精品免费| 亚洲成人久久性| 女警被强在线播放| 国产三级中文精品| aaaaa片日本免费| 欧美在线黄色| 亚洲不卡免费看| www.色视频.com| 观看免费一级毛片| 18禁国产床啪视频网站| 啦啦啦观看免费观看视频高清| 亚洲精品美女久久久久99蜜臀| a级一级毛片免费在线观看| 18禁在线播放成人免费| 欧美激情久久久久久爽电影| 免费人成视频x8x8入口观看| 国产精品1区2区在线观看.| 性色avwww在线观看| av在线蜜桃| 1024手机看黄色片| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 成人国产综合亚洲| 日日干狠狠操夜夜爽| 国产黄片美女视频| 欧美xxxx黑人xx丫x性爽| 亚洲电影在线观看av| 无人区码免费观看不卡| 精品久久久久久,| 国产久久久一区二区三区| 他把我摸到了高潮在线观看| 天堂网av新在线| 中文字幕熟女人妻在线| 成人永久免费在线观看视频| 日本免费a在线| 日本熟妇午夜| 最近最新中文字幕大全免费视频| 内射极品少妇av片p| 亚洲av一区综合| 国产熟女xx| 18禁在线播放成人免费| 久久久久免费精品人妻一区二区| 偷拍熟女少妇极品色| 激情在线观看视频在线高清| 青草久久国产| 久久久久亚洲av毛片大全| www国产在线视频色| 欧美日本亚洲视频在线播放| 操出白浆在线播放| 999久久久精品免费观看国产| 母亲3免费完整高清在线观看| 国产亚洲精品av在线| 19禁男女啪啪无遮挡网站| 久久国产精品人妻蜜桃| 日韩 欧美 亚洲 中文字幕| aaaaa片日本免费| 日本免费一区二区三区高清不卡| 欧美成人一区二区免费高清观看| 欧美精品啪啪一区二区三区| 最近视频中文字幕2019在线8| 一级毛片女人18水好多| 宅男免费午夜| 高清在线国产一区| 五月伊人婷婷丁香| 在线看三级毛片| 久久精品人妻少妇| 别揉我奶头~嗯~啊~动态视频| 日韩中文字幕欧美一区二区| 亚洲精品影视一区二区三区av| 日本熟妇午夜| 日本免费a在线| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久国产精品人妻蜜桃| 欧美午夜高清在线|