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

    不同溫度條件下亞熱帶森林土壤碳礦化對氮磷添加的響應(yīng)

    2018-01-04 05:57:35劉玉槐嚴員英張艷杰趙玉皓陳伏生葛體達魯順保
    生態(tài)學報 2017年23期
    關(guān)鍵詞:礦化排放量顯著性

    劉玉槐, 嚴員英, 張艷杰, 嚴 月, 趙玉皓, 徐 燕, 陳伏生, 葛體達, 魯順保,*

    1 江西師范大學生命科學學院,南昌 330022 2 江西省林業(yè)科學院濕地中心,南昌 330032 3 江西農(nóng)業(yè)大學林學院,南昌 330045 4 中國科學院亞熱帶農(nóng)業(yè)生態(tài)研究所區(qū)域農(nóng)業(yè)生態(tài)研究中心, 長沙 410125

    不同溫度條件下亞熱帶森林土壤碳礦化對氮磷添加的響應(yīng)

    劉玉槐1,4, 嚴員英2, 張艷杰1, 嚴 月1, 趙玉皓1, 徐 燕1, 陳伏生3, 葛體達4, 魯順保1,*

    1 江西師范大學生命科學學院,南昌 330022 2 江西省林業(yè)科學院濕地中心,南昌 330032 3 江西農(nóng)業(yè)大學林學院,南昌 330045 4 中國科學院亞熱帶農(nóng)業(yè)生態(tài)研究所區(qū)域農(nóng)業(yè)生態(tài)研究中心, 長沙 410125

    在氮沉降增加的背景下,土壤可利用性氮(N)素增加是多數(shù)陸地生態(tài)系統(tǒng)正經(jīng)歷的重要生態(tài)學過程。氮沉降增加對土壤呼吸會產(chǎn)生擾動,進而影響到森林土壤碳循環(huán)。以江西省泰和縣石溪杉木林為研究對象,開展N沉降模擬試驗,采用堿液吸收法分析4種不同N沉降處理在不同溫度(5、15、25、35℃)下對土壤碳排放量的影響。N沉降對土壤呼吸具有抑制作用,隨著溫度的升高,各處理對土壤呼吸碳排放呈顯著性增加趨勢(P< 0.05)。不同溫度條件下添加磷(P)處理對CO2累積排放量與添加高N處理變化規(guī)律類似。在5℃條件下,N會抑制土壤碳礦化,且隨N濃度升高抑制作用越強,P處理能促進土壤碳排放;15℃和 35℃條件下,N和P處理對土壤碳排放量均有抑制作用。在較高溫度(25℃和35℃)培養(yǎng)下,土壤全碳對CO2累積排放量均有顯著性影響。在25℃條件下,土壤DOC含量對土壤碳排放量的影響顯著(P< 0.05),其他溫度無影響。在5℃時,不同處理下土壤含水量對土壤碳排放量影響顯著。對N沉降與土壤Q10值進行方差分析表明不同處理對Q10值無顯著性差異。通過單庫模式方程Cm=Co(1-exp-kt)對土壤潛在碳排放進行模擬得出5、15℃和35℃對土壤潛在排放量有顯著影響。

    氮沉降;杉木;碳礦化;溫度

    土壤呼吸是陸地生態(tài)系統(tǒng)土壤和大氣之間CO2交換的重要組成部分。全球氣候變化引起氣溫升高加速了土壤有機碳在微生物和動物的礦化作用,其產(chǎn)物通過土壤呼吸作用以CO2形式釋放反饋回大氣圈,進一步加劇全球變暖,因此,土壤呼吸對環(huán)境變化反應(yīng)最為敏感[1-2],這也是目前全球變化研究的熱點問題,各國學者對此研究異?;钴S,他們正嘗試通過量化各種生態(tài)系統(tǒng)碳通量與環(huán)境因子之間的關(guān)系模型來預(yù)測未來全球變化對陸地生態(tài)系統(tǒng)的影響[3]。

    N、P是植物生長的必需元素,也是植物從土壤中吸收最多的元素。在土壤中99%以上的氮素是通過微生物礦化作用將其轉(zhuǎn)化為植物可直接吸收的有效態(tài)氮[4]。當前全球氣溫升高、降水格局改變以及大氣氮沉降增加等因素正在強烈改變著陸地生態(tài)系統(tǒng)的碳匯強度及其分布格局[5]。這些動態(tài)的環(huán)境因子迫使植被改變其小氣候進而影響土壤呼吸。大量的研究表明土壤水分和溫度敏感性對土壤碳礦化具有重要影響[6-7],施氮肥對土壤呼吸的影響也存在較大爭議,如Gallardo等[8]認為施氮肥能夠增加土壤呼吸,Kane等[9]的研究表明施氮肥對土壤呼吸速率沒有顯著的影響,Haynes等[10]對針葉林地施肥結(jié)果卻顯著降低土壤的呼吸速率。

    近年來,大氣氮沉降加重[11],并成為全球化趨勢[12],我國現(xiàn)已是世界第三大氮沉降集中區(qū)[13]。全球化背景下大氣氮沉降量的增加已成為許多森林的新生態(tài)因子[14],長期的高氮輸入可提高土壤碳礦化速率,導(dǎo)致土壤酸化。土壤呼吸對氣候變化(氮沉降、溫度和水分)又非常敏感,且氮沉降、溫度及水分對其土壤碳礦化影響機制尚未得到一致結(jié)論,氮沉降和溫度變化是亞熱帶地區(qū)面臨的重要生態(tài)環(huán)境問題,它們對森林地土壤碳礦化的影響如何?土壤呼吸對二者間的響應(yīng)是否存在差異?均具有重要科學意義。本文以杉木人工林為研究對象,擬通過模擬氣候變化試驗,研究該地區(qū)土壤碳礦化對氮沉降以及溫度敏感性的響應(yīng),研究結(jié)果將有利于了解氮沉降如何引起土壤碳排放的變化,為降低土壤CO2排放提供參考。

    1 研究地區(qū)與研究方法

    1.1 試驗地概況

    試驗地位于江西省泰和縣中國科學院千煙洲生態(tài)試驗站(26°42′N,115°04′E),該地區(qū)屬典型的中亞熱帶濕潤季風氣候,夏季炎熱干燥,冬季寒冷潮濕。年均相對濕度89%,年均溫度在18℃,6月最高溫達29.7℃,1月最低溫可達6.5℃,年均日照時數(shù)1400 h,年總輻射量約為4349 MJ/m2,年平均降水量約為1500 mm,降雨集中在4—6月,夏季降水量占全年降水量的50%以上。該區(qū)域為典型紅壤丘陵地貌,土壤母質(zhì)是第四紀紅粘土。研究區(qū)域大氣沉降向林地輸入的氮為82.8 kg hm-2a-1),干沉降為主,占67%,大氣氨沉降是大氣氮沉降主要貢獻者,占大氣干沉降82.0%,占大氣氮沉降總量的56%[15]。

    1.2 實驗設(shè)計及樣品采集

    試驗樣地起始時間為1990年,1991年4月開始種植植物。施肥自2012年春季起,每季度各施肥1次,施肥方式是將氮肥和磷肥拌少量細沙進行均勻撒施。土壤樣品采集時間為2015年1月,在立地條件一致的杉木人工林中,選取12個20 m×20 m的樣方,將這12個樣方隨機各選出3個樣方作為重復(fù),每個重復(fù)分別進行CK、N1、N2和P處理,4種處理分別為:CK為不施肥只施沙;N1:氮肥(NH4NO3)量為142.857 kg hm-2a-1(含氮量為50 kg N hm-2a-1);N2:285.714 kg hm-2a-1(含氮量為100 kg N hm-2a-1);P:磷肥(Na2HPO4)量為229.043 kg hm-2a-1(含磷量為50 kg P hm-2a-1),相鄰樣方之間留出一個間距為10 m的緩沖帶以減少不同施肥帶來的影響。在每個樣方中以S型選取5個點,分別用直徑7.5 cm的土鉆在0—10 cm取土壤樣品,然后將從同一樣方的土壤混合,每份土壤樣品約1000 g,裝入無菌自封袋帶回實驗室。將新鮮土壤樣品過2 mm篩去除土壤動物、小石頭、可見的植物殘體和根系等雜質(zhì),再將每份土壤樣品平均分為2份,一份用于測定土壤含水量、pH、全氮和全碳的含量,另一份放入4℃冰箱保存,用于測定土壤可溶性碳氮及生物活性。土壤基本理化性質(zhì)見表1。

    表1 土壤的基本理化性質(zhì)

    CK為不施肥只施沙;N1:含氮量為50 kg N hm-2a-1;N2:含氮量為100 kg N hm-2a-1;P:含磷量為50 kg P hm-2a-1。同列不同字母a、b、c表示不同處理間土壤理化性質(zhì)差異顯著(P< 0.05)

    1.3 實驗方法

    方法參照Chen等[16],采用元素分析儀(Isoprime-EuroEA3000,意大利)測定土壤全碳、全氮含量,SOC、SON含量采用2 mol/L KCl浸提,然后用TOCN 分析儀測定SOC 和可溶性總氮(total soluble N,TSN),pH值用1∶25水溶液、pH計測定,土壤碳礦化測試方法:先稱取新鮮土壤40.0 g,將其含水量調(diào)到田間持水量的60%,放入4個恒溫培養(yǎng)箱(溫度分別設(shè)為5、15、25、35℃)里進行42 d培養(yǎng)。土壤碳礦化測定方法采用堿液吸收法:2個70 mL的帶蓋塑料瓶置于1 L的容器內(nèi),其中一個為CO2吸收裝置,其內(nèi)加入用于吸收CO2的0.1 mol/L NaOH,另一個裝有40.0 g鮮土,在22℃條件下進行培養(yǎng),同時做5個對照試驗(僅放置堿液)。培養(yǎng)第1、3、7、14、21、28、35、42 d后將盛有0.1 mol/L NaOH吸收瓶取出后,先加入1 mol/L BaCl2使得土壤排放的CO2以碳酸鹽的形式沉淀,加入酚酞指示劑后攪拌,用0.05 mol/L HCl進行滴定剩余的NaOH。對土壤碳礦化進行模擬方程

    Cm=Co(1-exp-kt)

    式中,Cm為瞬時碳排放量,Co是土壤碳排放量極大值,k為常數(shù),t為培養(yǎng)時間。做潛在最高碳排放量[17],結(jié)合溫度敏感性Q10分析,即溫度每升高10℃土壤碳礦化速率所增加的倍數(shù)。

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

    實驗數(shù)據(jù)采用SPSS 17.0軟件進行分析,對不同處理間的差異顯著性進行雙因素方差分析,并用LSD法進行差異顯著性檢驗,并用SigmaPlot 12.0軟件進行數(shù)據(jù)分析和制圖,顯著性差異水平α=0.05。

    2 結(jié)果與分析

    2.1 氮沉降對土壤碳礦化的影響

    2.1.1 氮沉降對不同溫度條件下土壤累積碳排放量的影響

    從圖1可知,在5℃時,土壤碳礦化累積排放量隨氮濃度的升高呈下降趨勢;在15℃時,氮沉降對土壤碳礦化累積排放量具有一定的抑制作用,但P肥對土壤累積碳排放量具有一定的促進作用;在25℃時,氮沉降對土壤呼吸累積碳排放量均低于對照組,且低氮處理抑制作用更強;在35℃條件下,兩者氮沉降處理對土壤呼吸累積碳排放量均表現(xiàn)為一定的抑制作用。在同一溫度下,各處理間均無顯著性差異(P> 0.05)。隨著溫度的升高,對照CK和N、P添加處理下土壤碳礦化累積碳排放量顯著增加(P< 0.05,圖2)。

    圖1 NP添加對杉木林在相同溫度下土壤碳礦化量的影響Fig.1 Effects of substrate addition on CO2-C emission under Cunninghamia lanceolata in the same temperatureCK為不施肥只施沙;N1:含氮量為50 kg N hm-2 a-1;N2:含氮量為100 kg N hm-2 a-1;ns:差異不顯著( P > 0.05)。圖中數(shù)據(jù)為平均值 ± 標準誤

    從圖2可以看出,溫度升高顯著提升了林地土壤碳礦化量,不同處理間,不同溫度對土壤碳排放量均有一定的影響。在CK處理中,5℃和15℃間無顯著性差異,但均顯著低于25℃和35℃對土壤碳礦化的影響(P< 0.05),35℃對土壤碳排放量影響顯著高于25℃。N1處理中,于35℃時CO2累計排放量顯著高于其他溫度土壤碳排放量(P< 0.05),但5、15、25℃間無顯著差異。N2處理中,隨溫度升高CO2土壤累計排放量逐漸增加,其中5℃與15℃,15℃與25℃之間差異不顯著,但5℃顯著低于25℃、后者顯著低于35℃(P< 0.05)。添加P處理中,不同溫度對CO2累計排放量與N2變化規(guī)律類似,即35℃處理下碳排放量最大,其次為25℃和15℃,最小排放量為5℃處理。

    圖2 底物添加在不同溫度下對杉木林土壤碳礦化量的影響Fig.2 Effects of substrate addition on CO2-C evolved under Cunninghamia lanceolata in different temperature不同小寫字母表示同一處理不同溫度間差異顯著(P < 0.05),圖中數(shù)據(jù)為平均值 ± 標準差

    2.1.2 氮沉降對不同溫度條件下土壤碳礦化速率的影響

    由圖3可知,在同一溫度下,不同處理土壤礦化速率大致趨勢是先急速下降然后趨于平緩。在5℃時,由于溫度較低,各處理土壤礦化率偏低,N2、P二者呈先緩慢上升再趨于平緩的單峰曲線,CK、N1培養(yǎng)3 d內(nèi)呈急速下降,并且N1在前3天土壤礦化速率高于CK,P處理最低,在第7天后變化緩慢;在15℃時,對照CK在第7天略有升高外,其他處理都是先下降再趨于平緩;在25℃時,P處理后的第3天后土壤礦化速率開始下降,其他處理都是先下降再趨于穩(wěn)定,N1處理土壤礦化速率始終為最低;在35℃時,各處理都是先下降后趨于平穩(wěn),NP添加對土壤礦化速率影響的規(guī)律為:P > N1> N2。從圖3還可以看出,溫度對土壤碳礦化速率的影響顯著。

    2.2 環(huán)境因子對土壤累積碳排放量的影響

    2.2.1 土壤全碳對土壤累積碳排放量的影響

    從圖4可以看出,對土壤在不同溫度條件下,做土壤全碳和累積碳排放量進行分析得出,在5℃和15℃條件下,土壤全碳對土壤碳排放量無顯著性差異,在25℃條件下,土壤全碳與土壤碳排放量呈類似開口朝上的拋物線的變化趨勢,即隨土壤全碳含量的增加,土壤碳排放量先減小后增加。在35℃條件下,CO2累積排放量隨土壤全碳含顯著增加(P< 0.05),對其全碳和碳排放量進行回歸分析,得回歸方程為:y=-10.1595+143.201x,R2=0.38。

    圖3 不同底物添加對杉木林土壤礦化速率的影響Fig.3 Effects of different substrate addition on rate of Cunninghamia lanceolata soil respiration

    圖4 全碳對不同氮添加處理杉木林土壤碳礦化的影響Fig.4 Effects of total carbon for Cunninghamia lanceolata soil carbon minerlization under different nitrogen addition

    2.2.2 土壤可溶性碳氮對土壤累積碳排放量的影響

    從圖5可以看出,在5、15、35℃條件下,土壤DOC對土壤碳排放量影響無顯著性差異(P> 0.05)。在25℃條件下,土壤碳排放量隨DOC含量先增后減,呈類似開口朝下的拋物線變化規(guī)律,且差異達顯著性(P< 0.05)。在同一溫度條件下,土壤DON對土壤累積碳排放量均無顯著性差異(P> 0.05)。

    圖5 可溶性有機碳對不同氮添加處理杉木林土壤碳礦化的影響Fig.5 Effects of Dissolved Organic for Cunninghamia lanceolata soil carbon minerlization under different nitrogen addition

    2.2.3 土壤含水量對土壤累積碳排放量的影響

    從圖6可以看出,在5℃時,不同處理下杉木林土壤碳排放量隨土壤含水量增加而增加,且差異達顯著性水平(P< 0.05)。在其他溫度條件下并無顯著性差異,在同一含水量的條件下,土壤碳排放量隨溫度的升高而增加,說明土壤含水量并不是影響土壤碳排放量的調(diào)控因子。

    圖6 土壤含水量對不同氮添加處理杉木林土壤碳礦化的影響Fig.6 Effect of moisture for Cunninghamia lanceolata soil carbon minerlization under different nitrogen addition

    2.3 土壤礦化潛在排放量及其對溫度敏感性的響應(yīng)

    在沒有外界干擾的情況下長時間培養(yǎng)后土壤累積碳排放量將會達到一個峰值。對杉木林地土壤進行氮添加處理,在42 d的室內(nèi)培養(yǎng)實驗中,運用方程Cm=Co(1-exp-kt)對土壤潛在碳排放進行模擬,得出在不同溫度不同處理下杉木林土壤最大碳排放量(表2),相關(guān)系數(shù)R2基本在0.9以上,說明擬合效果很好。在5℃條件下,N會抑制土壤碳礦化,且隨N濃度升高抑制作用越強,P處理能促進土壤碳排放;15℃和35℃條件下,NP處理對土壤碳排放量均有抑制作用;在25℃條件下,不同處理對杉木林土壤潛在排放量影響不顯著。

    土壤礦化與溫度的關(guān)系還可以從溫度敏感性角度來研究(Q10),因為二氧化碳產(chǎn)生的溫度敏感性可作為有機質(zhì)降解溫度敏感性的主要指標。本研究中Q10培養(yǎng)42d的結(jié)果見表3,其范圍為1.5149—1.8354,從表3可以看出,N和P處理對土壤Q10值均有下降趨勢,其中N濃度升高,Q10值下降,對不同處理的Q10值做方差分析得出它們無顯著差異。

    3 討論

    土壤礦化是一個非常復(fù)雜的生物學過程,受土壤溫度、土壤含水量、土壤C、N含量等非生物因子以及植被類型、根系等生物因子綜合影響[18]。溫度是影響土壤礦化一個重要的因子,土壤礦化和溫度之間具有顯著的相關(guān)關(guān)系[19-20]。本研究結(jié)果表明,溫度升高能顯著提升林地土壤碳礦化量。在不同氮沉降處理下,溫度對土壤碳排放量均存在不同程度的影響。在對照處理中,25℃和35℃對土壤碳排放量影響顯著,且顯著高于5℃和15℃間對土壤碳礦化的影響(P< 0.05),可能原因是隨著溫度的升高,土壤中的微生物活性增加促進了土壤礦化,增溫在短時間內(nèi)可刺激土壤礦化產(chǎn)生大量 CO2[21-22]。在森林生態(tài)系統(tǒng)中,土壤碳礦化受溫度的影響較為復(fù)雜,主要是溫度對土壤生物活性的影響,往往同時還受土壤外源N添加的影響[20]。在N1處理中,35℃時CO2累計排放量顯著高于其他溫度對土壤碳排放量的影響(P< 0.05),且其他溫度間無顯著差異。在N2處理中,隨溫度升高CO2土壤累計排放量逐漸增加,且5℃與25、35℃間的差異顯著(P< 0.05),外源氮的輸入可能導(dǎo)致根系氮濃度的增加,進而影響土壤礦化作用[23]。在同一溫度下,與對照組相比,土壤碳排放量均隨氮濃度的升高而減少,說明氮沉降對土壤碳呼吸具有抑制作用,在25℃和35℃條件下,氮抑制呼吸作用,高溫度條件下抑制作用更強[24-25],主要原因是高氮處理降低土壤pH值,導(dǎo)致土壤微生物活性降低,分解有機碳的速率降低,從而抑制了土壤礦化[26],這種現(xiàn)象在針葉林和闊葉林中的影響更為明顯[27]。也可能是由于氮素的增加降低土壤C/N進而改變土壤微生物的呼吸進而影響土壤礦化[28-30],也有學者認為氮添加對土壤礦化的無影響或促進作用[31-32]。在不同溫度條件下,土壤有機碳礦化速率呈日趨下降,培養(yǎng)21 d后礦化速率趨于平緩,其原因可能是土壤可利用性碳、氮基質(zhì)的量是影響其礦化的主要因素,因此,在礦化培養(yǎng)實驗中,前期土壤易礦化的有機碳數(shù)量較多,已被微生物分解,培養(yǎng)到中后期土壤碳礦化主要受到可利用碳含量的限制,微生物開始難礦化有機碳,使得分解速率下降,CO2釋放量開始降低[33]。較高溫度培養(yǎng)過的樣品 CO2釋放速率顯著低于較低溫度,很多長期培養(yǎng)試驗中,都觀測到基質(zhì)大量損耗后,試驗后期土壤礦化速率的降低[34],Yoshitake等發(fā)現(xiàn)添加氮在開始會刺激土壤礦化,但不會持續(xù)下去[26],這也與本研究中土壤礦化率在后期降低并趨于穩(wěn)定相一致[35]。隨著溫度的升高,CK和N、P添加在短時間內(nèi)對土壤礦化碳排放呈顯著性增加趨勢(P< 0.05)。添加P處理中,不同溫度對CO2累計排放量與N2變化規(guī)律類似,可能的原因是P對土壤微生物比較敏感抑制土壤礦化[36-37]。

    表2不同溫度下各處理土壤潛在碳礦化最大排放量及參數(shù)

    Table2ComparisonofcalculatedpotentiallymineralizableorganicCpoolsCoandfirstorderrateconstantsKofdifferenttreatmentsatdifferenttemperaturesinsoil

    溫度/℃Temperature處理Treatments潛在碳排放量Potentiallycarbonmineralizable參數(shù)Constant(K)相關(guān)系數(shù)Relationcoefficient(R2)5CK138.663±22.858ab0.032±0.005a0.941±0.018aN1134.848±18.117ab0.026±0.007ab0.921±0.040aN277.133±1.532b0.037±0.020a0.960±0.014aP206.517±104.477a0.016±0.009b0.925±0.023a15CK315.946±21.335a0.054±0.034b0.887±0.007aN185.397±20.987bc0.102±0.079a0.954±0.015aN274.193±18.464c0.062±0.040ab0.840±0.040aP111.495±12.358b0.098±0.085a0.680±0.139b25CK176.355±22.848a0.062±0.004a0.955±0.017aN1174.745±26.953a0.041±0.010b0.987±0.005aN2187.938±21.846a0.041±0.001b0.987±0.002aP186.869±6.041a0.053±0.007a0.983±0.012a35CK488.887±35.656a0.025±0.002b0.991±0.003aN1348.949±64.454ab0.033±0.002b0.972±0.011aN2315.799±48.418b0.041±0.007a0.981±0.004aP477.196±130.299a0.031±0.009a0.974±0.021a

    同一列不同小寫字母表示同一溫度不同處理間差異顯著(P< 0.05),圖中數(shù)據(jù)為平均值 ± 標準差

    表3 第42天的Q10值

    有關(guān)土壤礦化速率與溫度的關(guān)系研究比較多,他們之間擬合的曲線方程較為復(fù)雜。本研究從溫度敏感性角度來分析(Q10) 土壤礦化速率受溫度的影響,Q10培養(yǎng)范圍為1.5149—1.8354,變異比較小,說明氮沉降的影響表明氮添加處理對杉木人工林土壤有機碳礦化對溫度變化的敏感性不是特別明顯,這與王紅等研究的結(jié)果不一致[38],其原因可能是土壤礦化是一個十分復(fù)雜的生物學過程,Q10值除受溫度的影響外,還受土壤水分、根生物量及其他因子的綜合結(jié)果[28],因此,氮添加對土壤礦化溫度敏感性影響的機理需要進一步細化研究。總之,未來氣候變化對亞熱帶森林生態(tài)系統(tǒng)碳循環(huán)的影響及其機理會更為復(fù)雜,這將是今后的主要研究方向之一。

    在森林生態(tài)系統(tǒng)中,土壤礦化是土壤碳庫主要輸出途徑,因此,土壤碳庫含量的高低對土壤碳礦化作用產(chǎn)生重要影響。從土壤全碳與土壤碳礦化的關(guān)系來看,在5℃和15℃條件下,土壤全碳對土壤碳排放量無顯著性差異,在25℃和35℃條件下,土壤全碳對土壤碳排放量的影響均達差異顯著(P< 0.05),說明土壤全碳對土壤碳礦化影響同時也受溫度調(diào)控。土壤碳庫以及土壤氮素作為土壤有機碳礦化的底物,其含量必將直接對土壤有機碳礦化作用產(chǎn)生影響[39-40]。但生態(tài)系統(tǒng)中非生物因子和生物因子間相互關(guān)系的復(fù)雜性,使得土壤有機碳礦化研究仍存在不確定性,如Giardina等認為溫度的升高不會對土壤有機碳礦化產(chǎn)生影響[41],本研究結(jié)果表明在25℃條件下,土壤碳排放量隨DOC含量先增后減(P< 0.05),而其他溫度條件下,土壤DOC對土壤碳排放量并無顯著性影響。且在同一溫度條件下,土壤DON對土壤累積碳排放量均無顯著性差異。土壤礦化除受土壤碳庫影響外,可能還與土壤含水量密切相關(guān),由于土壤含水量對呼吸速度影響較為復(fù)雜,往往同時取決于溫度的配置狀況。在一定溫度條件下,土壤含水量適度提高不僅可以增強微生物的活性,也可以增加土壤中微生物的數(shù)量,從而加快土壤碳的礦化速度。本研究結(jié)果表明,在5℃時,杉木林土壤碳礦化量隨土壤含水量的增加而顯著增加(P< 0.05),在其他溫度條件下并無顯著性差異。在同一含水量的條件下,土壤碳排放量隨溫度的升高而增加,這與結(jié)論與Hartley研究的結(jié)果一致[35]。溫度對土壤有機碳礦化作用的影響可能會因土壤有機碳庫以及土壤水分等因子的限制作用而減弱,并且對于土壤水分以及土壤碳、氮含量等因子對土壤有機碳礦化的影響也存在較大的爭議[40]。因此,土壤碳礦化對溫度的響應(yīng)是土壤性質(zhì)、生物活性以及可利用性碳、氮綜合作用的結(jié)果[33]。

    目前,室內(nèi)培養(yǎng)方法廣泛應(yīng)用于土壤潛在碳排放的研究,由于土壤礦化在沒有根系或氣候因子等干擾的情況下,土壤礦化在長時間培養(yǎng)后終將達到平衡。本研究采用單庫模式一級動力學方程Cm=Co(1-exp-kt)對土壤礦化潛在碳排放進行模擬,得出相關(guān)系數(shù)均在0.9以上。在25℃條件下,不同處理對杉木林土壤潛在排放量影響不顯著。在其他溫度條件下,N添加均會抑制土壤潛在碳排放量,且隨N濃度升高抑制作用越強(P< 0.05)。

    4 結(jié)論

    氮沉降對土壤礦化碳的累積排放量具有抑制作用。隨著溫度的升高,CK、N和P添加對土壤礦化碳排放呈顯著性增加趨勢(P< 0.05)。不同處理間,不同溫度對土壤碳排放量均有一定的影響。不同溫度條件下添加P處理中對CO2累計排放量與N2變化規(guī)律類似。在5℃條件下,N會抑制土壤碳礦化,且隨N濃度升高抑制作用越強,P處理能促進土壤碳排放;15℃和 35℃條件下,NP處理對土壤碳排放量均有抑制作用;25℃時不同處理對土壤潛在排放量影響不顯著。

    在較高溫度(25℃和35℃)培養(yǎng)下,土壤全碳對CO2累積排放量均有顯著性影響,較低的溫度(5℃和15℃)培養(yǎng)下,土壤全碳對土壤碳排放量無顯著性差異。在25℃條件下,土壤DOC含量對土壤碳排放量的影響顯著(P< 0.05),其他溫度無影響。在同一溫度條件下,土壤DON對土壤累積碳排放量均無顯著性差異。在5℃時,不同處理下土壤含水量對土壤碳排放量影響顯著,在其他溫度條件下并無顯著性差異。

    氮沉降和P處理對土壤Q10值與對照相比均有下降趨勢,隨N濃度升高Q10值下降,對不同處理的Q10值做方差分析得出它們無顯著差異。

    [1] 姜艷, 王兵, 汪如玉, 楊清培. 亞熱帶林分土壤呼吸及其與土壤溫濕度關(guān)系的模型模擬. 應(yīng)用生態(tài)學報, 2010, 21(7): 1641- 1648.

    [2] Zerva A,Mencuccini M. Short-term effects ofclearfelling on soil CO2, CH4, and N2O fluxes in a Sitka spruce plantation. Soil Biology and Biochemistry, 2005, 37(11): 2025- 2036.

    [3] Li S G, Asanuma J,Eugster W,Kotani A, Liu J J,Urano T, Oikawa T,Davaa G,Oyunbaatar D, Sugita M. Net ecosystem carbon dioxide exchange over grazed steppe in central Mongolia. Global Change Biology, 2005, 11(11): 1941- 1955.

    [4] 欒軍偉, 向成華, 駱宗詩, 宮淵波. 森林土壤呼吸研究進展. 應(yīng)用生態(tài)學報, 2006, 17(12): 2451- 2456.

    [5] Piao S L, Fang J Y,Ciais P, Peylin P, Huang Y, Sitch S, Wang T. The carbon balance of terrestrial ecosystems in China. Nature, 2009, 458(7241): 1009- 1013.

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

    [7] Xu M, Qi Y. Soil-surface CO2efflux and its spatial and temporal variations in a young ponderosa pine plantation in northern California. Global Change Biology, 2001, 7(6): 667- 677.

    [8] Gallardo A, Schlesinger W H. Factors limiting microbial biomass in the mineral soil and forest floor of a warm-temperate forest. Soil Biology and Biochemistry, 1994, 26(10): 1409- 1415.

    [9] Kane E S,Pregitzer K S, Burton A J. Soil respiration along environmental gradients in Olympic National Park. Ecosystems, 2003, 6(4): 326- 335.

    [10] Haynes B E, Gower S T. Belowground carbon allocation in unfertilized and fertilized red pine plantations in northern Wisconsin. Tree Physiology, 1995, 15(5): 317- 325.

    [11] Vitousek P M, Aber J D,Howarth R W, Likens E L, Matson P A, Schindler D W, Schlesinger W H, Tilman D G. Human alteration of the global nitrogen cycle: sources and consequences. Ecological Applications, 1997, 7(3): 737- 750.

    [12] Matson P, Lohse K A, Hall S J. The globalization of nitrogen deposition: consequences for terrestrial ecosystems.Ambio: A Journal of the Human Environment, 2001, 31(2): 113- 119.

    [13] Mo J M, Brown S,Xue J H, Fang Y T, Li Z. Response of litter decomposition to simulated N deposition in disturbed, rehabilitated and mature forests in subtropical China. Plant and Soil, 2006, 282(1/2): 135- 151.

    [14] Kazda M. Indications of unbalanced nitrogen nutrition of Norway spruce stands. Plant and Soil, 1990, 128(1): 97- 101.

    [15] 樊建凌, 胡正義, 莊舜堯, 周靜, 王體健, 劉翠英. 林地大氣氮沉降的觀測研究. 中國環(huán)境科學, 2007, 27(1): 7- 9.

    [16] Chen C R, Xu Z H, Mathers N J. Soil carbon pools in adjacent natural and plantation forests of subtropical Australia. Soil Science Society of America Journal, 2004, 68(1): 282- 291.

    [17] Smith J L,Schobel R R, McNeal B L, Campbell G S. Potential errors in the first-order model for estimating soil nitrogen mineralization potentials. Soil Science Society of America Journal, 1980, 44(5): 996- 1000.

    [18] 李玉寧, 王關(guān)玉, 李偉. 土壤呼吸作用和全球碳循環(huán). 地學前緣, 2002, 9(2): 351- 357.

    [19] 嚴俊霞, 湯億, 李洪建. 城市綠地土壤呼吸與土壤溫度土壤水分的關(guān)系研究. 干旱區(qū)地理, 2009, 32(4): 604- 609.

    [20] Menyailo O V,Huwe B. Denitrification and C, N mineralization as function of temperature and moisture potential in organic and mineral horizons of an acid spruce forest soil. Journal of Plant Nutrition and Soil Science, 1999, 162(5): 527- 531.

    [21] 楊慶朋, 徐明, 劉洪升, 王勁松, 劉麗香, 遲永剛, 鄭云普. 土壤呼吸溫度敏感性的影響因素和不確定性. 生態(tài)學報, 2011, 31(8): 2301- 2311.

    [22] Mikan C J, Schimel J P, Doyle A P. Temperature controls of microbial respiration in arctic tundra soils above and below freezing. Soil Biology and Biochemistry, 2002, 34(11): 1785- 1795.

    [23] Eissenstat D M, Wells C E,Yanai R D,Whitbeck J L. Building roots in a changing environment: implications for root longevity. New Phytologist, 2000, 147(1): 33- 42.

    [24] 吳浩浩, 徐星凱, 李團勝, 段存濤. 碳氮添加對不同濕潤程度的溫帶森林土壤氧化亞氮排放影響. 科學通報, 2014, 59(14): 1337- 1347.

    [25] Zhao S C,Qiu S J, Cao C Y, Zheng C L, Zhou W, He P. Responses of soil properties, microbial community and crop yields to various rates of nitrogen fertilization in a wheat-maize cropping system in north-central China. Agriculture, Ecosystems & Environment, 2014, 194: 29- 37.

    [26] Yoshitake S, Sasaki A, Uchida M,Funatsu Y, Nakatsubo T. Carbon and nitrogen limitation to microbial respiration and biomass in an acidicsolfatara field. European Journal of Soil Biology, 2007, 43(1): 1- 13.

    [27] 王效科,白艷瑩,歐陽志云,苗鴻. 全球碳循環(huán)中的失匯及其形成原因. 生態(tài)學報, 2001, 22(1): 94- 103.

    [28] Janssens I A,Dieleman W,Luyssaert S,Subke J A, Reichstein M,Ceulemans R,Ciais P, Dolman A J, Grace J,Matteucci G,Papale D,Piao S L, Schulze E D, Tang J, Law B E. Reduction of forest soil respiration in response to nitrogen deposition. Nature Geoscience, 2010, 3(5): 315- 322.

    [29] Traore S,Thiombiano L,Millogo J R,Guinko S. Carbon and nitrogen enhancement incambisols andvertisols byAcaciaspp. in eastern Burkina Faso: relation to soil respiration and Microbial biomass. Applied Soil Ecology, 2007, 35(3): 660- 669.

    [30] Persson T,Wirén A. Microbial activity in forest soils in relation to acid/base and carbon/nitrogen status. Communications of the Norwegian Forest Research Institute, 1989, 81(suppl): S15- 21.

    [31] Cheshire M V, Chapman S J. Influence of the N and P status of Plant material and of added N and P on the mineralization of C from 14C-labelled ryegrass in soil. Biology and Fertility of Soils, 1996, 21(3): 166- 170.

    [32] Bradford M A,Fierer N, Reynolds J F. Soil carbon stocks in experimentalmesocosms are dependent on the rate of labile carbon, nitrogen and Phosphorus inputs to soils. Functional Ecology, 2008, 22(6): 964- 974.

    [33] 沈征濤, 施斌, 王寶軍, 姜洪濤. 土壤有機質(zhì)轉(zhuǎn)化及CO2釋放的溫度效應(yīng)研究進展. 生態(tài)學報, 2013, 33(10): 3011- 3019.

    [34] Hartley I P, Hopkins D W, Sommerkorn M, Wookey P A. The response of organic mattermineralisation to nutrient and substrate additions in sub-arctic soils. Soil Biology and Biochemistry, 2010, 42(1): 92- 100.

    [35] Jarvis P, Linder S. Botany: constraints to growth of boreal forests. Nature, 2000, 405(6789): 904- 905.

    [36] Keith H, Jacobsen K L, Raison R J. Effects of soil phosphorus availability, temperature and moisture on soil respiration in Eucalyptuspauciflora forest. Plant and Soil, 1997, 190(1): 127- 141.

    [37] 王紅, 范志平, 鄧東周, 陳德龍, 孫學凱, 高俊剛, 曾德慧. 不同環(huán)境因子對樟子松人工林土壤有機碳礦化的影響. 生態(tài)學雜志, 2008, 27(9): 1469- 1475.

    [38] Davidson E A,Janssens I A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature, 2006, 440(7081): 165- 173.

    [39] Hopkins D W, Sparrow A D,Elberling B,Gregorich E G,Novis P M, Greenfield L G,Tilston E L. Carbon, nitrogen and temperature controls on microbial activity in soils from an Antarctic dry valley. Soil Biology and Biochemistry, 2006, 38(10): 3130- 3140.

    [40] 楊繼松, 劉景雙, 孫麗娜. 溫度、水分對濕地土壤有機碳礦化的影響. 生態(tài)學雜志, 2008, 27(1): 38- 42.

    [41] Giardina C P, Ryan M G. Evidence that decomposition rates of organic carbon in mineral soil do not vary with temperature. Nature, 2000, 404(6780): 858- 861.

    EffectsofnitrogenandphosphorusadditionunderdifferenttemperaturesonthesoilcarbonmineralizationinaCunninghamialanceolataplantationinthesubtropics

    LIU Yuhuai1,4, YAN Yuanying2, ZHANG Yanjie1, YAN Yue1, ZHAO Yuhao1, XU Yan1, CHEN Fusheng3, GE Tida4, LU Shunbao1,*

    1CollegeofLifeSciences,JiangxiNormalUniversity,Nanchang330022,China2CentreforWetlandJiangxiForestryAcademy,Nanchang330032,China3CollegeofForestry,JiangxiAgriculturalUniversity,Nanchang330045,China4KeyLaboratoryofAgro-ecologicalProcessesinSubtropicalRegion,InstituteofSubtropicalAgriculture,ChineseAcademyofSciences,Changsha410125,China

    Increase in the soil nitrogen (N) availability under increasing N deposition is an important ecological process that most terrestrial ecosystems experience. Increasing N deposition could affect soil respiration and consequently, affect the soil carbon (C) and N cycling. In this study, an N deposition simulation experiment was carried out in theCunninghamialanceolataforest in Shixi, Taihe County, Jiangxi Province, China. Soil carbon emissions were analyzed, using an alkali absorption method, under four treatments (no N and Phosphorus (P) addition (CK), N1 addition, N2 addition, and P addition) at four different temperatures (5℃, 15℃, 25℃, and 35℃). Soil respiration was inhibited by N deposition. As the temperature increased, soil respiration increased substantially (P< 0.05). Different treatments had different effects on soil carbon emissions. The change regulation of CO2accumulation emission and high N change regulation were similar for P addition under different temperature conditions. At 5℃, N addition inhibited soil carbon mineralization and this effect was stronger with increasing N concentration, whereas P addition promoted the soil carbon emission with increasing N concentration. However, at 15℃ and 35℃, N and P addition had an inhibitory effect on the soil carbon emission. At higher temperatures (25℃ and 35℃), total soil carbon significantly affected the CO2accumulation; however, at lower temperatures (5℃ and 15℃) total soil carbon had the opposite effect on CO2accumulation. At 5℃, soil water content in different treatments had a significant effect on soil carbon emission, whereas the other three temperatures (15℃, 25℃, and 35℃) had no significant effect. At 25℃, soil carbon emission was significantly affected by soil DOC (P< 0.05), whereas no effect was observed at the other three soil temperatures (5℃, 15℃, and 35℃). At 25℃, soil DON had no significant effect on soil carbon accumulation emission. At 5℃, soil water content in different treatments had a significant effect on soil carbon emission but no significant effect at the other three temperatures (15℃, 25℃, and 35℃). The Q10value under N and P addition was lower than that under the control condition (no N and P addition) and theQ10value decreased with increasing N concentration. However, the difference in N deposition had no significant effect on theQ10values. The potential carbon emission from soil respiration was simulated by the single reservoir model,Cm=Co(1-exp-kt). The results indicated that the four treatments had a significant effect on potential soil carbon emission at 5℃, 15℃, and 35℃, but no significant effect at 25℃.

    N deposition;Cunninghamialanceolataplantation; carbon mineralization; temperature

    國家自然科學基金(31360136, 31560168, 31660072);江西省自然基金(20151BAB204007, 20161BAB204175)

    2016- 10- 27; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 08- 14

    *通訊作者Corresponding author.E-mail: luxunbao8012@126.com

    10.5846/stxb201610272195

    劉玉槐, 嚴員英, 張艷杰, 嚴月, 趙玉皓, 徐燕, 陳伏生, 葛體達, 魯順保.不同溫度條件下亞熱帶森林土壤碳礦化對氮磷添加的響應(yīng).生態(tài)學報,2017,37(23):7994- 8004.

    Liu Y H, Yan Y Y, Zhang Y J, Yan Y, Zhao Y H, Xu Y, Chen F S, Ge T D, Lu S B.Effects of nitrogen and phosphorus addition under different temperatures on the soil carbon mineralization in aCunninghamialanceolataplantation in the subtropics.Acta Ecologica Sinica,2017,37(23):7994- 8004.

    猜你喜歡
    礦化排放量顯著性
    礦化劑對硅酸鹽水泥煅燒的促進作用
    大麥蟲對聚苯乙烯塑料的生物降解和礦化作用
    天然氣輸配系統(tǒng)甲烷排放量化方法
    煤氣與熱力(2021年6期)2021-07-28 07:21:40
    黑龍江省碳排放量影響因素研究
    基于顯著性權(quán)重融合的圖像拼接算法
    電子制作(2019年24期)2019-02-23 13:22:26
    基于視覺顯著性的視頻差錯掩蓋算法
    一種基于顯著性邊緣的運動模糊圖像復(fù)原方法
    論商標固有顯著性的認定
    全國機動車污染物排放量
    ——《2013年中國機動車污染防治年報》(第Ⅱ部分)
    江蘇省火力發(fā)電機組二氧化碳排放量估算
    天堂网av新在线| 中文字幕精品亚洲无线码一区| 九九在线视频观看精品| 淫妇啪啪啪对白视频| 亚洲美女搞黄在线观看 | 国产主播在线观看一区二区| 国语自产精品视频在线第100页| 精品午夜福利视频在线观看一区| 欧美丝袜亚洲另类 | 欧美高清成人免费视频www| 人妻久久中文字幕网| 国产亚洲av嫩草精品影院| 偷拍熟女少妇极品色| 精品免费久久久久久久清纯| 丰满的人妻完整版| 两个人的视频大全免费| 国产精品电影一区二区三区| 级片在线观看| 51国产日韩欧美| 久99久视频精品免费| 日韩高清综合在线| 99热网站在线观看| 成人高潮视频无遮挡免费网站| 欧美日韩国产亚洲二区| 在线观看舔阴道视频| 成人精品一区二区免费| 99久久成人亚洲精品观看| 婷婷丁香在线五月| 欧美中文日本在线观看视频| 51国产日韩欧美| 精品一区二区三区视频在线| 久久久久九九精品影院| 日本a在线网址| 少妇裸体淫交视频免费看高清| 极品教师在线免费播放| 国产亚洲91精品色在线| 久久久久久伊人网av| 搡女人真爽免费视频火全软件 | 男女那种视频在线观看| 国产大屁股一区二区在线视频| 男人狂女人下面高潮的视频| 午夜激情欧美在线| 国产美女午夜福利| x7x7x7水蜜桃| 久久久久性生活片| 亚洲av免费高清在线观看| 可以在线观看的亚洲视频| 乱系列少妇在线播放| 免费观看的影片在线观看| 嫩草影院入口| 久久午夜福利片| 日本欧美国产在线视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品色激情综合| 亚洲精品色激情综合| АⅤ资源中文在线天堂| 天堂√8在线中文| 亚洲精品国产成人久久av| 国产成人福利小说| 国产精品自产拍在线观看55亚洲| 丰满人妻一区二区三区视频av| 哪里可以看免费的av片| 亚洲内射少妇av| 欧美一区二区国产精品久久精品| 最近视频中文字幕2019在线8| 国产精品自产拍在线观看55亚洲| 亚洲欧美日韩卡通动漫| 国产黄a三级三级三级人| 精品人妻视频免费看| 日韩欧美一区二区三区在线观看| 九九在线视频观看精品| 看片在线看免费视频| 99久久精品一区二区三区| 成人精品一区二区免费| 老熟妇乱子伦视频在线观看| 国产一区二区三区视频了| 1000部很黄的大片| av中文乱码字幕在线| 婷婷色综合大香蕉| 亚洲一区二区三区色噜噜| 亚洲一区二区三区色噜噜| 我的老师免费观看完整版| 欧美一区二区精品小视频在线| 人妻夜夜爽99麻豆av| 国产精品国产高清国产av| 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情国产日韩精品一区| 亚洲精华国产精华液的使用体验 | 午夜亚洲福利在线播放| 国产精品自产拍在线观看55亚洲| 日韩欧美三级三区| 九九久久精品国产亚洲av麻豆| 日韩欧美精品v在线| 婷婷精品国产亚洲av在线| 日韩欧美精品v在线| 无遮挡黄片免费观看| 色在线成人网| 麻豆久久精品国产亚洲av| 亚洲第一区二区三区不卡| 精品一区二区三区视频在线观看免费| 免费av不卡在线播放| 搞女人的毛片| 日日夜夜操网爽| 联通29元200g的流量卡| 少妇的逼水好多| 成人欧美大片| 搡老熟女国产l中国老女人| 日韩精品中文字幕看吧| 欧美+日韩+精品| 色5月婷婷丁香| 波多野结衣巨乳人妻| 91狼人影院| 小说图片视频综合网站| 免费无遮挡裸体视频| 女同久久另类99精品国产91| 麻豆久久精品国产亚洲av| 动漫黄色视频在线观看| 久久热精品热| 国产精品无大码| 久久久久久久久大av| 深夜精品福利| 露出奶头的视频| 黄片wwwwww| 国产精品1区2区在线观看.| 色哟哟哟哟哟哟| 亚洲精品亚洲一区二区| 99国产极品粉嫩在线观看| 黄片wwwwww| 一进一出好大好爽视频| 国产精品人妻久久久影院| 免费不卡的大黄色大毛片视频在线观看 | 日韩,欧美,国产一区二区三区 | 精品久久久久久久人妻蜜臀av| 久久久久久久久久成人| 亚洲色图av天堂| 欧美区成人在线视频| 免费无遮挡裸体视频| 欧美日韩综合久久久久久 | eeuss影院久久| 熟女电影av网| 99热精品在线国产| 精品99又大又爽又粗少妇毛片 | 九色国产91popny在线| 久久久久性生活片| 99久久无色码亚洲精品果冻| 成年女人永久免费观看视频| 一a级毛片在线观看| 91久久精品电影网| 在线天堂最新版资源| 亚洲va日本ⅴa欧美va伊人久久| 赤兔流量卡办理| 久久久久精品国产欧美久久久| 一区二区三区免费毛片| 在线观看美女被高潮喷水网站| 97热精品久久久久久| 精品无人区乱码1区二区| 日本三级黄在线观看| 啦啦啦观看免费观看视频高清| 嫩草影院新地址| 一区二区三区四区激情视频 | 91午夜精品亚洲一区二区三区 | 三级毛片av免费| 亚洲男人的天堂狠狠| 国产精品野战在线观看| 日本熟妇午夜| 日日啪夜夜撸| 亚洲综合色惰| 噜噜噜噜噜久久久久久91| 色播亚洲综合网| 99热这里只有是精品在线观看| 国内久久婷婷六月综合欲色啪| 亚洲中文日韩欧美视频| 亚洲经典国产精华液单| 又粗又爽又猛毛片免费看| 黄色女人牲交| 俄罗斯特黄特色一大片| 午夜精品久久久久久毛片777| 日本三级黄在线观看| 亚洲性久久影院| 日本-黄色视频高清免费观看| av中文乱码字幕在线| 国产精品,欧美在线| 欧美性猛交╳xxx乱大交人| a在线观看视频网站| 男人和女人高潮做爰伦理| 可以在线观看的亚洲视频| 噜噜噜噜噜久久久久久91| 欧美xxxx黑人xx丫x性爽| 亚洲国产色片| 美女 人体艺术 gogo| 99热6这里只有精品| 九色国产91popny在线| 伦精品一区二区三区| 人妻夜夜爽99麻豆av| av国产免费在线观看| 99九九线精品视频在线观看视频| 18禁黄网站禁片午夜丰满| 又黄又爽又刺激的免费视频.| 欧美区成人在线视频| 国产高清激情床上av| 一夜夜www| 精品欧美国产一区二区三| 精品久久久久久成人av| 99热精品在线国产| 小说图片视频综合网站| 久久精品夜夜夜夜夜久久蜜豆| 美女黄网站色视频| 亚洲va在线va天堂va国产| 蜜桃久久精品国产亚洲av| 欧美另类亚洲清纯唯美| 国产综合懂色| 欧美日韩综合久久久久久 | www.色视频.com| 69av精品久久久久久| 很黄的视频免费| 亚洲天堂国产精品一区在线| 国产精品一区二区免费欧美| 亚洲av五月六月丁香网| 免费看a级黄色片| 精品乱码久久久久久99久播| 亚洲五月天丁香| 直男gayav资源| 亚洲,欧美,日韩| 久久久午夜欧美精品| 亚洲熟妇中文字幕五十中出| 免费人成视频x8x8入口观看| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 日韩欧美三级三区| 一区二区三区免费毛片| 性欧美人与动物交配| 美女被艹到高潮喷水动态| 日韩精品中文字幕看吧| 日本与韩国留学比较| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 国产精品久久电影中文字幕| 人人妻人人看人人澡| 日韩欧美 国产精品| 1000部很黄的大片| 亚洲最大成人中文| 久久热精品热| 免费高清视频大片| 国产老妇女一区| 精品久久久噜噜| 老司机深夜福利视频在线观看| 1000部很黄的大片| 亚洲精品粉嫩美女一区| 亚洲精华国产精华精| xxxwww97欧美| 亚洲午夜理论影院| 不卡一级毛片| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6 | 国产毛片a区久久久久| av在线蜜桃| 亚洲av.av天堂| 国产精品爽爽va在线观看网站| 欧美xxxx性猛交bbbb| 亚洲图色成人| 欧美最新免费一区二区三区| 成人av一区二区三区在线看| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 成人一区二区视频在线观看| 精品福利观看| 欧美bdsm另类| 色播亚洲综合网| 国产精品久久电影中文字幕| .国产精品久久| 日日撸夜夜添| 99久久精品国产国产毛片| 国产精品久久电影中文字幕| 一个人看的www免费观看视频| 色噜噜av男人的天堂激情| 亚洲黑人精品在线| 高清在线国产一区| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 日韩在线高清观看一区二区三区 | 欧美精品国产亚洲| 桃红色精品国产亚洲av| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 黄色女人牲交| 久久国内精品自在自线图片| 九色国产91popny在线| 亚洲人成网站在线播放欧美日韩| 五月玫瑰六月丁香| 一本久久中文字幕| 国产在视频线在精品| 我的老师免费观看完整版| 亚洲国产欧美人成| 熟女人妻精品中文字幕| 国产精品一区二区免费欧美| 国内精品美女久久久久久| 香蕉av资源在线| 久久久久久久午夜电影| 久99久视频精品免费| 国产精品不卡视频一区二区| 少妇人妻一区二区三区视频| 国产黄色小视频在线观看| 国产单亲对白刺激| 色综合色国产| 国产不卡一卡二| 欧美在线一区亚洲| 午夜久久久久精精品| 一进一出抽搐gif免费好疼| 狂野欧美激情性xxxx在线观看| 如何舔出高潮| 小说图片视频综合网站| 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 国产精品电影一区二区三区| 男人舔奶头视频| а√天堂www在线а√下载| 欧美国产日韩亚洲一区| 免费观看在线日韩| 五月玫瑰六月丁香| 18禁黄网站禁片午夜丰满| 韩国av一区二区三区四区| 少妇人妻一区二区三区视频| 久久午夜亚洲精品久久| 色尼玛亚洲综合影院| 日韩中文字幕欧美一区二区| 国产av不卡久久| 一级黄片播放器| 白带黄色成豆腐渣| 国产精品国产三级国产av玫瑰| 最近在线观看免费完整版| 日日摸夜夜添夜夜添小说| 国产午夜精品久久久久久一区二区三区 | 老女人水多毛片| 欧美日韩黄片免| 国产男人的电影天堂91| 一级黄片播放器| 老司机深夜福利视频在线观看| 成人国产麻豆网| 午夜老司机福利剧场| 亚洲av成人精品一区久久| 在线观看av片永久免费下载| 噜噜噜噜噜久久久久久91| 蜜桃亚洲精品一区二区三区| 亚洲经典国产精华液单| 国产男人的电影天堂91| 国产人妻一区二区三区在| 国产伦在线观看视频一区| 男女下面进入的视频免费午夜| 毛片女人毛片| 日韩一区二区视频免费看| 日本与韩国留学比较| 搡老妇女老女人老熟妇| 在线播放国产精品三级| 国产高清视频在线播放一区| 免费大片18禁| 高清在线国产一区| 国产免费一级a男人的天堂| 乱人视频在线观看| 欧美成人免费av一区二区三区| 亚洲,欧美,日韩| 香蕉av资源在线| 国产麻豆成人av免费视频| 黄色配什么色好看| 日韩欧美一区二区三区在线观看| 国产人妻一区二区三区在| 黄色日韩在线| 美女高潮喷水抽搐中文字幕| 亚洲人成伊人成综合网2020| 国产乱人伦免费视频| 88av欧美| 亚洲成a人片在线一区二区| 欧美激情久久久久久爽电影| 日韩人妻高清精品专区| 超碰av人人做人人爽久久| 午夜久久久久精精品| 亚洲中文字幕一区二区三区有码在线看| 91久久精品国产一区二区三区| 99热这里只有精品一区| 成人av一区二区三区在线看| 日韩欧美免费精品| 色尼玛亚洲综合影院| 成人毛片a级毛片在线播放| 嫩草影院新地址| 成人综合一区亚洲| 3wmmmm亚洲av在线观看| 国产精品嫩草影院av在线观看 | av天堂中文字幕网| 国产高清激情床上av| 午夜免费激情av| 国产伦精品一区二区三区四那| 国产一区二区在线观看日韩| 99精品久久久久人妻精品| h日本视频在线播放| 日韩强制内射视频| 日韩中文字幕欧美一区二区| 噜噜噜噜噜久久久久久91| 高清在线国产一区| 久久久午夜欧美精品| 久久久久九九精品影院| 欧美日本视频| 精品国产三级普通话版| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av涩爱 | 人妻丰满熟妇av一区二区三区| 五月伊人婷婷丁香| 国产亚洲av嫩草精品影院| 亚洲无线观看免费| 一进一出好大好爽视频| 色视频www国产| 三级男女做爰猛烈吃奶摸视频| 一级毛片久久久久久久久女| 九色国产91popny在线| 亚洲国产精品sss在线观看| 在现免费观看毛片| 黄色配什么色好看| 亚洲国产欧美人成| 日本三级黄在线观看| 亚洲精品一区av在线观看| 亚洲四区av| 伦理电影大哥的女人| 国产一区二区三区在线臀色熟女| 日韩欧美国产在线观看| 人人妻,人人澡人人爽秒播| 一区二区三区免费毛片| 日日摸夜夜添夜夜添av毛片 | 天堂影院成人在线观看| 国产真实伦视频高清在线观看 | 午夜福利在线观看免费完整高清在 | 成年免费大片在线观看| 亚洲色图av天堂| 亚洲av美国av| 精品无人区乱码1区二区| 午夜福利在线观看免费完整高清在 | 乱码一卡2卡4卡精品| 午夜影院日韩av| 真人一进一出gif抽搐免费| 成人国产综合亚洲| 久久精品综合一区二区三区| 中文在线观看免费www的网站| 亚洲精品成人久久久久久| 三级毛片av免费| 国产免费一级a男人的天堂| 日韩av在线大香蕉| 国产精品人妻久久久久久| 欧美黑人欧美精品刺激| 如何舔出高潮| 国产av在哪里看| 麻豆国产av国片精品| 亚洲欧美日韩东京热| 我要看日韩黄色一级片| 免费观看精品视频网站| 一区福利在线观看| 免费在线观看成人毛片| 九九在线视频观看精品| 欧美zozozo另类| 在线国产一区二区在线| 久久人人爽人人爽人人片va| 精品日产1卡2卡| 亚洲真实伦在线观看| 精品一区二区三区视频在线| 国产成人aa在线观看| 99热这里只有精品一区| 国产精品无大码| 在线看三级毛片| 国产乱人视频| 男插女下体视频免费在线播放| 99久久九九国产精品国产免费| 欧美日韩中文字幕国产精品一区二区三区| 最近最新中文字幕大全电影3| 国产不卡一卡二| 非洲黑人性xxxx精品又粗又长| 动漫黄色视频在线观看| 性插视频无遮挡在线免费观看| 免费人成在线观看视频色| 国产免费av片在线观看野外av| 老司机福利观看| 久久人人精品亚洲av| 一级黄片播放器| 22中文网久久字幕| 丝袜美腿在线中文| 熟妇人妻久久中文字幕3abv| 国产一区二区三区在线臀色熟女| 窝窝影院91人妻| 变态另类丝袜制服| 亚洲在线自拍视频| 他把我摸到了高潮在线观看| 日本三级黄在线观看| 日韩精品有码人妻一区| 老司机福利观看| 97人妻精品一区二区三区麻豆| 久久久久性生活片| 午夜福利在线在线| 亚洲欧美清纯卡通| 中文亚洲av片在线观看爽| 一区二区三区免费毛片| 日本成人三级电影网站| 日本免费a在线| 人人妻人人看人人澡| 九九热线精品视视频播放| 91麻豆av在线| 日本与韩国留学比较| 亚洲精品456在线播放app | 舔av片在线| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看| 国产精品一区www在线观看 | 色精品久久人妻99蜜桃| 在线a可以看的网站| 日本 欧美在线| 婷婷亚洲欧美| 精品免费久久久久久久清纯| 欧美+日韩+精品| 在线观看美女被高潮喷水网站| 国产精品亚洲美女久久久| 国内毛片毛片毛片毛片毛片| 成人永久免费在线观看视频| 欧美日韩乱码在线| 亚洲va日本ⅴa欧美va伊人久久| 制服丝袜大香蕉在线| 欧美+日韩+精品| 免费观看人在逋| 国产探花极品一区二区| 欧美性猛交╳xxx乱大交人| 欧美日本视频| 国产精品精品国产色婷婷| 精品久久久久久久人妻蜜臀av| 女生性感内裤真人,穿戴方法视频| 国产精品一区二区免费欧美| 亚洲av免费在线观看| 日韩欧美国产在线观看| 久久精品夜夜夜夜夜久久蜜豆| 免费高清视频大片| 国产精品永久免费网站| 久久99热这里只有精品18| 国产精品永久免费网站| 国产成年人精品一区二区| 淫妇啪啪啪对白视频| 精品久久国产蜜桃| 久久久久久久久中文| 一级毛片久久久久久久久女| 欧美+日韩+精品| 久久久精品欧美日韩精品| 久久久久免费精品人妻一区二区| 欧美日韩乱码在线| 少妇熟女aⅴ在线视频| 一夜夜www| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频 | 国产真实伦视频高清在线观看 | 国产不卡一卡二| 日韩欧美在线二视频| 日韩精品有码人妻一区| 亚洲男人的天堂狠狠| 久久99热这里只有精品18| 老师上课跳d突然被开到最大视频| 在线免费观看的www视频| 黄色丝袜av网址大全| 亚洲精品一区av在线观看| 男人狂女人下面高潮的视频| 少妇被粗大猛烈的视频| 欧美高清性xxxxhd video| 他把我摸到了高潮在线观看| 国产色爽女视频免费观看| 亚洲熟妇中文字幕五十中出| 成人精品一区二区免费| 日本精品一区二区三区蜜桃| 真实男女啪啪啪动态图| 日本精品一区二区三区蜜桃| a级一级毛片免费在线观看| 中文字幕熟女人妻在线| 一个人看视频在线观看www免费| 岛国在线免费视频观看| 国产精品98久久久久久宅男小说| 国产亚洲av嫩草精品影院| 91久久精品电影网| 精品一区二区三区视频在线| 亚洲av中文av极速乱 | 国产精品一区二区性色av| 精品一区二区三区视频在线| 97碰自拍视频| 亚洲国产色片| 日本五十路高清| 极品教师在线视频| ponron亚洲| 亚洲人与动物交配视频| 精品人妻视频免费看| 亚洲精品乱码久久久v下载方式| 国产高清激情床上av| 精品人妻一区二区三区麻豆 | av在线蜜桃| 非洲黑人性xxxx精品又粗又长| 亚洲成人精品中文字幕电影| 欧美日韩亚洲国产一区二区在线观看| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆 | 91久久精品电影网| 久久99热6这里只有精品| 久久精品国产自在天天线| 国产 一区 欧美 日韩| 精品久久久久久久久av| 日日摸夜夜添夜夜添av毛片 | 天天一区二区日本电影三级| 久久久久久久精品吃奶| 日韩中文字幕欧美一区二区| 乱系列少妇在线播放| 2021天堂中文幕一二区在线观| 啦啦啦啦在线视频资源| 丰满的人妻完整版| 国内精品久久久久精免费| 色综合站精品国产| 一个人看视频在线观看www免费| 成人无遮挡网站|