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

    模擬氮沉降對(duì)落葉松人工林土壤呼吸的影響

    2016-11-10 01:03:47郝龍飛王慶成付嬌嬌朱凱月
    植物研究 2016年4期
    關(guān)鍵詞:中氮異養(yǎng)落葉松

    張 嬌 郝龍飛 王慶成 付嬌嬌 朱凱月

    (東北林業(yè)大學(xué)林學(xué)院,哈爾濱 150040)

    模擬氮沉降對(duì)落葉松人工林土壤呼吸的影響

    張 嬌 郝龍飛 王慶成*付嬌嬌 朱凱月

    (東北林業(yè)大學(xué)林學(xué)院,哈爾濱 150040)

    在東北林業(yè)大學(xué)帽兒山實(shí)驗(yàn)林場(chǎng)26年生落葉松人工林中,連續(xù)2年(2013~2014年)施加NH4NO3模擬氮沉降試驗(yàn)((對(duì)照(CK,0 g·m-2·a-1N)、低氮(N1,5 g·m-2·a-1N)、中氮(N2,10 g·m-2·a-1N)、高氮(N3,15 g·m-2·a-1N)),研究不同氮沉降水平對(duì)土壤呼吸的影響。結(jié)果表明:(1)2013年模擬氮沉降處理均促進(jìn)年平均土壤呼吸速率(P<0.05);(2)2014年中氮和高氮處理抑制年平均土壤呼吸和異養(yǎng)呼吸速率(P<0.05),低氮處理促進(jìn)年均土壤呼吸速率(P<0.05),對(duì)異養(yǎng)呼吸速率影響不顯著(P>0.05);(3)土壤微生物生物量碳在低氮處理下顯著提高(P<0.05),在中氮和高氮處理下與對(duì)照間差異不顯著(P>0.05);(4)土壤呼吸速率與5和10 cm土壤溫度呈指數(shù)正相關(guān)關(guān)系(P<0.01),相比對(duì)照,各土層土壤呼吸溫度敏感系數(shù)(Q10)均在低氮處理下增加,在中氮和高氮處理下則降低。不同水平的模擬氮沉降改變了土壤呼吸速率及其溫度敏感性,表明短期內(nèi)低水平氮沉降可加快土壤碳排放過(guò)程,相對(duì)較高水平氮沉降則減緩?fù)寥捞寂欧胚^(guò)程。

    落葉松(Larixolgensis);人工林;氮沉降;土壤呼吸;根呼吸;異養(yǎng)呼吸

    森林生態(tài)系統(tǒng)是陸地上最大的碳庫(kù),約占全球植物碳儲(chǔ)量和全球土壤碳儲(chǔ)量的86%和73%,其中土壤碳庫(kù)是其重要的組成部分,超過(guò)植被碳庫(kù)的2倍以上[1~2]。土壤呼吸是土壤與大氣碳交換的主要過(guò)程,也是全球土壤碳循環(huán)重要環(huán)節(jié)。全球每年通過(guò)土壤呼吸作用釋放的C為75 Pg[3],即使土壤呼吸速率發(fā)生很小的變化,也會(huì)在很大程度上影響大氣CO2濃度和土壤碳吸存的變化。以往研究表明氮素是影響土壤呼吸的主要因素之一[4]。近幾十年來(lái),受化肥生產(chǎn)、化石燃料燃燒以及人類生產(chǎn)活動(dòng)的影響,大氣中含氮化合物激增,導(dǎo)致大氣氮沉降明顯增加[5]。大氣氮沉降進(jìn)入森林生態(tài)系統(tǒng),影響其土壤呼吸過(guò)程和土壤碳儲(chǔ)量[6~7]已引起國(guó)內(nèi)外學(xué)者的廣泛關(guān)注,并進(jìn)行了一系列模擬氮沉降對(duì)林地土壤呼吸影響的試驗(yàn)。研究結(jié)果不盡相同,部分研究表明氮輸入抑制土壤呼吸[8~9],然而也有研究表明氮輸入促進(jìn)土壤呼吸[10~12]或無(wú)明顯影響[13]。

    我國(guó)學(xué)者針對(duì)氮沉降主要區(qū)域的土壤呼吸變化進(jìn)行了一系列研究,提高了我們對(duì)氮沉降與土壤呼吸響應(yīng)的認(rèn)識(shí)。涂利華等[14~15]對(duì)華西雨屏區(qū)的研究發(fā)現(xiàn)氮沉降促進(jìn)苦竹和撐綠雜交竹林分土壤呼吸;孫素琪等[16]和黃玉梓等[17]在縉云山常綠闊葉林和杉木人工林的研究發(fā)現(xiàn)氮沉降抑制土壤呼吸;李凱等[18]和李偉斌等[19]的研究結(jié)果發(fā)現(xiàn)氮沉降存在閾值效應(yīng),低水平氮促進(jìn)土壤呼吸,高水平氮抑制土壤呼吸;莫江明等[20]在鼎湖山保護(hù)區(qū)的研究發(fā)現(xiàn)短期(90 d)內(nèi)的氮沉降顯著促進(jìn)季風(fēng)林的土壤呼吸,且這種促進(jìn)作用隨氮水平提高而增大。目前為止,我國(guó)氮沉降研究主要集中在中部和東南沿海地區(qū),而東北地區(qū)研究相對(duì)缺乏[21]。在東北林業(yè)大學(xué)帽兒山的研究結(jié)果表明施氮肥降低了落葉松和水曲柳人工林的土壤呼吸速率,但研究?jī)H設(shè)置了一個(gè)施肥處理(10 g·m-2·a-1N),無(wú)法了解該區(qū)不同氮沉降水平下林下土壤呼吸的變化規(guī)律。

    落葉松(Larixolgensis)是我國(guó)東北地區(qū)重要的造林樹(shù)種,其人工林面積居該區(qū)人工林面積之首。我國(guó)作為全球三大氮沉降區(qū)之一,2012年年均氮沉降量達(dá)到12.9 kg·hm2·a-1N[22]。近年來(lái),東北地區(qū)氮沉降量逐漸增加,濕氮沉降量為21.5 kg·hm2·a-1N,已達(dá)到國(guó)內(nèi)重度氮沉降的二級(jí)水平(20.01~25.00 kg·hm2·a-1N)[21]。因此研究東北地區(qū)落葉松人工林土壤呼吸對(duì)氮沉降的響應(yīng)具有重要意義。本文通過(guò)模擬不同氮沉降水平,研究落葉松人工林土壤呼吸速率變化趨勢(shì),為大氣氮沉降持續(xù)增加背景下東北地區(qū)人工林碳動(dòng)態(tài)變化提供參考。

    1 材料與方法

    1.1 研究地區(qū)概況

    研究地位于東北林業(yè)大學(xué)帽兒山實(shí)驗(yàn)林場(chǎng)尖砬溝森林培育實(shí)驗(yàn)站(E127°30′~127°34′,N45°21′~45°25′)。低山丘陵地貌,平均海拔300 m。大陸性季風(fēng)氣候。年降水量600~800 mm,年蒸發(fā)量1 093 mm,年平均濕度70%;年平均溫度2.7℃,1月平均溫度-19.6℃,7月平均溫度20.9℃,全年≥10℃的積溫2 638℃;年平均總?cè)照諘r(shí)數(shù)1 856.8 h;無(wú)霜期約120~140 d。土壤為地帶性暗棕色森林土,平均土層厚度在40 cm左右[23],有機(jī)質(zhì)含量在A層最高,自A1層向下明顯降低;質(zhì)地為壤土;土壤整體呈微酸性,pH值為5.5~6.0[24]。

    試驗(yàn)林為落葉松人工林,栽植于1987年。于2013年5月對(duì)其進(jìn)行調(diào)查,林分密度為975株·hm-2,平均樹(shù)高19.0 m,平均胸徑15.7 cm;0~20 cm土壤pH為5.57,總碳、總氮和全磷分別為53.8 g·kg-1,5.3 g·kg-1和1.21 mg·g-1,碳氮比為10.17。

    1.2 試驗(yàn)設(shè)計(jì)與研究方法

    2013年5月,在落葉松人工林中設(shè)置3塊標(biāo)準(zhǔn)地(20 m×20 m),樣地間間隔10 m以上。參照國(guó)內(nèi)外模擬氮沉降研究[10,25~26]以及該地區(qū)氮沉降現(xiàn)狀(2.15 g·m-2·a-1N),在每個(gè)樣地內(nèi)均勻設(shè)置4個(gè)5 m×5 m的模擬氮沉降試驗(yàn)地,呈田字排列,各試驗(yàn)地間隔5 m。試驗(yàn)地處在林下自然環(huán)境中,均接受正常的大氣氮沉降量,同時(shí)額外添加氮素模擬不同水平大氣氮沉降,添加量分別為對(duì)照(CK,0 g·m-2·a-1N)、低氮(N1,5 g·m-2·a-1N)、中氮(N2,10 g·m-2·a-1N)和高氮(N3,15 g·m-2·a-1N),每個(gè)處理3次重復(fù)。于2013和2014年生長(zhǎng)季內(nèi)(5~9月)每月初進(jìn)行1次模擬氮沉降施肥處理。采用2 L清水溶解氮肥(NH4NO3),用手持式噴霧器在對(duì)應(yīng)氮沉降試驗(yàn)地內(nèi)均勻噴灑,對(duì)照噴灑清水。

    1.2.1 土壤呼吸測(cè)定

    2013年5月初,在每個(gè)氮沉降樣地地表隨機(jī)布設(shè)5個(gè)內(nèi)徑和高均為10 cm的PVC環(huán),共60個(gè),將PVC環(huán)一端削尖,壓入土壤中,露出地面2 cm左右,同時(shí)保持連接環(huán)在整個(gè)測(cè)定期間位置不變[27],用于定期測(cè)定土壤呼吸速率。測(cè)定土壤呼吸速率的同時(shí),用數(shù)字瞬時(shí)溫度計(jì)(11025,Delta TRAK,Inc.,USA)測(cè)定土壤(5 cm,10 cm)溫度,用土壤水分速測(cè)儀(TDR-300,Spectrum Technologies,Inc.,USA)測(cè)定土壤(5 cm,10 cm)濕度。

    分別于2013年6~10月和2014年5~10月,每月初采用動(dòng)態(tài)氣室法使用便攜式LI-6400土壤碳通量自動(dòng)測(cè)量系統(tǒng)(Li-Cor Inc.,Lincoln,NE,USA)測(cè)定土壤呼吸速率[28],測(cè)定時(shí)間選定在14~17時(shí),各氮沉降處理固定測(cè)定15個(gè)點(diǎn),每個(gè)點(diǎn)3~5次重復(fù),取平均值為該點(diǎn)呼吸值。土壤呼吸溫度敏感系數(shù)Q10值計(jì)算公式:Q10=e10b,其中b是土壤呼吸與溫度單因素指數(shù)曲線模型RS=aebT中的溫度反應(yīng)常數(shù)(RS為土壤呼吸速率,a為溫度為0℃時(shí)的土壤呼吸速率,b為溫度反應(yīng)常數(shù))。

    1.2.2 根呼吸和異養(yǎng)呼吸測(cè)定

    2014年8月采用離體PVC管氣室法[29],用LI-6400土壤碳通量自動(dòng)測(cè)量系統(tǒng)對(duì)根系呼吸進(jìn)行測(cè)定。每個(gè)處理3次重復(fù),在試驗(yàn)地內(nèi)用鐵鍬挖取0~20 cm的土壤,取出林木根系,迅速去除根系土,將直徑小于2 mm的細(xì)根(直徑大于2 mm的根系活性較低[30],土壤呼吸速率較低,本研究中忽略不計(jì))放置在封住底部的PVC連接環(huán)中,連接土壤碳通量自動(dòng)測(cè)量系統(tǒng)測(cè)定根系呼吸速率。將測(cè)定呼吸的根系帶回實(shí)驗(yàn)室烘干并測(cè)定質(zhì)量,得出單位質(zhì)量根系呼吸速率,結(jié)合樣地內(nèi)根生物量密度可計(jì)算出單位面積單位時(shí)間內(nèi)根系的呼吸速率,計(jì)算公式如下:

    VR=VTR/MR*S*DR

    (1)

    式中:VTR為測(cè)定的根系總呼吸速率(μmol·m-2·s-1);MR為用于測(cè)定根系呼吸速率的根系質(zhì)量(g);DR為氮沉降試驗(yàn)地內(nèi)的根生物量密度(g·m-2);S為PVC呼吸環(huán)面積(0.008 1 m2)。

    用土壤總呼吸速率和根系自養(yǎng)呼吸速率差值估算微生物和動(dòng)物的異養(yǎng)呼吸速率(設(shè)置標(biāo)準(zhǔn)地時(shí)已清除標(biāo)準(zhǔn)地內(nèi)灌木和草本,故排出其產(chǎn)生根呼吸的可能性)。異養(yǎng)呼吸速率計(jì)算公式:

    VH=VT-VR

    (2)

    式中:VH為異養(yǎng)呼吸速率(μmol·m-2·s-1),VT為土壤總呼吸速率(μmol·m-2·s-1),VR為根呼吸速率(μmol·m-2·s-1)。

    1.2.3根生物量密度、氮濃度和土壤微生物生物量碳測(cè)定

    2014年8月,測(cè)定細(xì)根呼吸同時(shí),在氮沉降試驗(yàn)地中測(cè)定林分的根生物量密度。在試驗(yàn)地中隨機(jī)選取5個(gè)點(diǎn),共60個(gè),用內(nèi)徑53 mm的土鉆采集0~20 cm土樣,裝袋標(biāo)號(hào),帶回野外實(shí)驗(yàn)站進(jìn)行流水沖洗,除去泥土和雜物,保留落葉松直徑小于2 mm的根系,并根據(jù)根系外形、顏色和彈性特征剔除死根,于75℃下烘干至恒重,測(cè)定生物量,從而計(jì)算各處理下根生物量密度[23]。稱重后將根系粉碎過(guò)60目土壤篩,通過(guò)元素分析儀(Vario Macro,Elementer,German)測(cè)定根系氮濃度。

    隨后在各氮沉降試驗(yàn)地內(nèi)隨機(jī)選取5個(gè)樣點(diǎn),共60個(gè),采集0~20 cm土壤,挑除植物根系和石礫,裝袋標(biāo)號(hào),低溫保存帶回實(shí)驗(yàn)室用于測(cè)定土壤微生物生物量碳,選用氯仿熏蒸方法(0.5 mol·L-1K2SO4)提取(做熏蒸和未熏蒸2種處理)[31],通過(guò)碳氮分析儀(multi N/C 2100,Analytikjena,German)進(jìn)行測(cè)定。

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

    用Excel(Microsoft Office 2007)對(duì)原始數(shù)據(jù)進(jìn)行整合計(jì)算和正則化處理,運(yùn)用SPSS(SPSS PASW Statistics,v18.0)軟件檢驗(yàn)數(shù)據(jù)方差齊性,通過(guò)單因素方差分析檢驗(yàn)不同氮沉降處理下土壤呼吸速率月動(dòng)態(tài)變化的差異,采用最小顯著性差異法(LSD)檢驗(yàn)差異顯著性(α=0.05)。用SigmaPlot(Systat Sigma Plot,v10.0)軟件作圖。

    圖1 連續(xù)2年模擬氮沉降對(duì)26年生落葉松人工林年均土壤呼吸速率的影響 CK、N1、N2、N3分別表示對(duì)照(CK:0 g·m-2·a-1 N)、低氮(N1:5 g·m-2·a-1 N)、中氮(N2:10 g·m-2·a-1 N)和高氮(N3:15 g·m-2·a-1 N),不同字母表示各處理間差異顯著(P<0.05),下同。Fig.1 Average annual soil respiration rate under simulated N deposition of 2 years in a 26-year-old L.olgensis plantation forest CK,N1,N2,N3 indicate control,low N,intermediate N and high N,respectively;different letters indicate significant difference among treatments(P<0.05),the same as below.

    2 結(jié)果與分析

    2.1 氮沉降對(duì)落葉松人工林土壤呼吸的影響

    2013年,落葉松人工林各氮沉降處理下的年平均土壤呼吸速率均顯著高于對(duì)照(P<0.05),各處理間差異不顯著(P>0.05)(圖1)。各處理下土壤呼吸速率均隨時(shí)間發(fā)生規(guī)律性變化,6~8月土壤呼吸速率逐漸升高,氮沉降處理與對(duì)照間的差異逐漸增大;8月土壤呼吸速率均達(dá)到最高,氮沉降處理與對(duì)照間的差異也達(dá)到最大(P>0.05);9~10月土壤呼吸速率迅速降低,各處理間差異縮小(圖2)。

    圖2 連續(xù)2年模擬氮沉降對(duì)26年生落葉松人工林土壤呼吸速率月動(dòng)態(tài)的影響Fig.2 Monthly soil respiration rate under simulated N deposition of 2 years in a 26-year-old L.olgensis plantation forest

    圖3 2014年8月模擬氮沉降對(duì)26年生落葉松人工林根系和土壤異養(yǎng)呼吸速率的影響Fig.3 Root and heterotrophic respiration rate under simulated N deposition in a 26-year-old L.olgensis plantation forest in August 2014

    2014年,隨模擬氮沉降量增加,年平均土壤呼吸速率在低氮處理下顯著升高(P<0.05),在中氮和高氮處理下顯著降低(P<0.05)(圖1)。5~8月土壤呼吸速率逐漸升高,到8月達(dá)到最大值,而后逐漸降低,呈明顯單峰型(圖2)。5月各處理間的土壤呼吸速率差異不顯著(P>0.05)。6~9月中氮和高氮處理的土壤呼吸速率均顯著低于低氮處理和對(duì)照(P<0.05),各月份內(nèi)N2和N3分別比CK降低24.55%、21.3%、25.67%、18.88%(P<0.05)和28.2%、28.84%、32.63%、42.97%(P<0.05);低氮處理的土壤呼吸速率均高于對(duì)照(P>0.05),各月份內(nèi)N1分別比CK提高12.63%、12.11%、13.01%、5.86%。10月,土壤呼吸速率迅速降低,各處理間差異縮小,高氮處理顯著低于低氮處理(P<0.05)(圖2)。

    2.2氮沉降對(duì)落葉松人工林根呼吸和異養(yǎng)呼吸的影響

    落葉松人工林根系呼吸速率在各處理間差異不顯著(P>0.05),與對(duì)照相比,氮沉降處理下的根呼吸速率均有提高的趨勢(shì),提高程度以中氮處理較大(圖3)。根系生物量密度和氮濃度隨氮沉降增加均有所改變,二者在低氮處理下與對(duì)照差異不顯著(P>0.05);中氮和高氮處理下根密度顯著降低(P<0.05),根系氮濃度顯著升高(P<0.05)(表1)。

    表12014年8月模擬氮沉降對(duì)26年生落葉松人工林根生物量密度和根系氮濃度的影響

    Table1RootbiomassdensityandNconcentrationundersimulatedNdepositionina26-year-oldL.olgensisplantationforestinAugust2014

    處理Treatment密度Biomassdensity(g·m-2)氮濃度Nconcentration(g·kg-1)CK50.79±3.39a10.99±0.35bN153.10±3.67a11.13±0.18bN234.54±3.39b12.44±0.22aN334.96±1.81b12.84±0.40a

    注:表中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)誤差(n=15),不同字母表示各處理間差異顯著(P<0.05)。

    Note:Data for mean±standard(n=15);different letters indicate significant difference among treatments(P<0.05).

    氮沉降處理對(duì)土壤異養(yǎng)呼吸速率的影響較大,在中氮和高氮處理下顯著降低(P<0.05),在低氮處理下與對(duì)照間差異不顯著(P=0.32),但土壤異養(yǎng)呼吸速率有升高的趨勢(shì)(圖3)。各處理下土壤微生物生物量碳變化也存在差異,隨氮沉降量增加呈先升高后降低趨勢(shì),低氮處理顯著提高土壤微生物生物量碳(P<0.05),中氮和高氮處理與對(duì)照無(wú)顯著差異(P>0.05)(圖4)。

    圖4 2014年8月模擬氮沉降對(duì)26年生落葉松人工林土壤微生物生物量碳的影響Fig.4 Soil microbial biomass carbon under simulated N deposition in a 26-year-old L.olgensis plantation forest in August 2014

    2.3氮沉降處理下落葉松人工林土壤呼吸與土壤溫濕度的關(guān)系

    各氮沉降處理下土壤呼吸速率均與5和10 cm土壤溫度呈指數(shù)正相關(guān)關(guān)系(P<0.01)(表2),且變化趨勢(shì)與土壤溫度季節(jié)性變化一致,均在8月達(dá)到最高值(圖2)。不同處理下土壤呼吸的溫度敏感性有所差異,相比對(duì)照,兩個(gè)土層的Q10值均在低氮處理下增加,在中氮和高氮處理下降低,各處理下土壤10 cm處的Q10值均高于5 cm處(表2)。

    表2模擬氮沉降下26年生落葉松人工林土壤呼吸速率與土壤溫度的指數(shù)模型和Q10值

    Table2ExponentialmodelsofsoilrespirationandsoiltemperatureandQ10undersimulatedNdepositionina26-year-oldL.olgensisplantationforest

    處理Treatment土壤深度Soildepth(cm)abR2PQ10CK50.56540.07140.6902<0.00012.04100.56250.07640.7078<0.00012.15N150.48800.07660.7497<0.00012.15100.40110.08950.6902<0.00012.45N250.53330.06730.6922<0.00011.96100.48710.07340.7224<0.00012.08N350.61960.06540.6548<0.00011.92100.59210.07220.6753<0.00012.06

    注:指數(shù)模型:RS=aebTRS.土壤呼吸速率;a,b.回歸系數(shù);e.常數(shù)(2.718281828);T.土壤溫度

    Note:Exponential model:RS=aebTRS. Soil respiration rate;a,b. Regression coefficient; e. Constant;T. Soil temperature

    不同氮沉降處理下土壤呼吸速率與各土層土壤濕度無(wú)顯著相關(guān)(P>0.05)。

    3 討論

    3.1 氮沉降對(duì)落葉松人工林土壤呼吸的影響

    模擬氮沉降初期(2013年),各月份土壤呼吸速率在氮沉降處理下與對(duì)照間差異不顯著(P>0.05)(圖2),而年平均土壤呼吸速率顯著高于對(duì)照(P<0.05)(圖1),說(shuō)明短期的氮沉降能在一定程度上促進(jìn)土壤呼吸。以往研究發(fā)現(xiàn)短期氮添加提高興安落葉松人工林[32]以及季風(fēng)林[20]的土壤呼吸速率,與本研究結(jié)果一致。氮沉降處理促進(jìn)土壤呼吸可能與森林生態(tài)系統(tǒng)土壤的氮素水平和狀態(tài)有關(guān)[33],可能原因是施氮提高土壤有效氮水平加強(qiáng)微生物對(duì)氮素的固定、提高細(xì)根對(duì)氮素的吸收利用率以及加速土壤氮礦化,進(jìn)而提高土壤呼吸強(qiáng)度。

    氮沉降一年后(2014年),低氮處理仍促進(jìn)土壤呼吸(P<0.05),而中氮和高氮處理則顯著抑制土壤呼吸(P<0.05),并且抑制效應(yīng)隨氮沉降水平提高而增加(圖2),說(shuō)明氮沉降對(duì)土壤呼吸的影響存在閾值效應(yīng)[19,34],與以往研究結(jié)果一致[8,10~11,19]。氮沉降處理顯著影響土壤呼吸速率,可能與土壤呼吸組分的變化有關(guān)。土壤呼吸作為土壤碳排放的主要過(guò)程之一,包括自養(yǎng)呼吸(根呼吸)和異養(yǎng)呼吸(微生物呼吸)兩部分,二者對(duì)氮沉降的總響應(yīng)決定土壤呼吸的變化。本研究中,氮沉降對(duì)根呼吸速率無(wú)顯著影響(P>0.05)(圖3),說(shuō)明氮沉降處理下根呼吸的變化對(duì)土壤呼吸變化的貢獻(xiàn)不大。以往研究表明根呼吸速率受根生物量以及氮濃度共同影響[35],并隨細(xì)根生物量密度[35]和氮濃度[36~37]的增加而提高。本研究中氮和高氮處理根密度顯著低于對(duì)照和低氮處理(P<0.05),而氮濃度顯著高于對(duì)照和低氮處理(P<0.05)(表1)。因此,氮沉降對(duì)根呼吸速率無(wú)顯著影響可能與細(xì)根生物量密度降低和氮濃度升高的共同作用有關(guān)。本研究中,土壤異養(yǎng)呼吸速率在中氮和高氮處理下受到抑制,在低氮處理下與對(duì)照差異不顯著,但有升高的趨勢(shì),且變化趨勢(shì)同土壤呼吸一致(圖2~3),說(shuō)明異養(yǎng)呼吸的變化對(duì)土壤呼吸變化的貢獻(xiàn)較大。異養(yǎng)呼吸主要是土壤微生物活動(dòng)產(chǎn)生的,陸地生態(tài)系統(tǒng)中,土壤微生物大都生存在碳豐富而氮不足的環(huán)境中[37]。本研究土壤碳氮比為10,而土壤細(xì)菌的碳氮比一般為6[38],相比微生物對(duì)碳和氮的需求,土壤氮素存在限制,進(jìn)而抑制微生物的活性[39~40]并降低其生物量[41]。低氮處理提高土壤中的氮含量,緩解氮素對(duì)微生物的限制,促使微生物活性提高(圖3)[39~40],同時(shí)微生物生物量碳顯著增加(圖4)[41],二者的共同作用促使異養(yǎng)呼吸速率提高,從而間接提高土壤呼吸速率。然而生態(tài)系統(tǒng)的氮輸入存在閾值效應(yīng)[34],隨氮沉降量的增加,氮素輸入達(dá)到閾值從而影響土壤微生物和林分根系生長(zhǎng),中氮和高氮處理的微生物生物量碳與對(duì)照相間無(wú)顯著差異(圖4),表明高水平氮沉降并未改變微生物生物量,然而中氮和高氮處理下根生物量密度顯著降低(表1),導(dǎo)致根際有機(jī)質(zhì)分解和釋放減少,進(jìn)一步抑制土壤微生物活性[10],降低異養(yǎng)呼吸速率,進(jìn)而導(dǎo)致土壤呼吸速率降低。李偉斌等在長(zhǎng)白山闊葉紅松林中模擬氮沉降試驗(yàn),研究發(fā)現(xiàn)土壤呼吸速率在低氮(2.3 g·m-2·a-1N)和中氮(4.6 g·m-2·a-1N)處理下均有所提高,在高氮(6.9 g·m-2·a-1N)處理下降低,也認(rèn)為氮沉降對(duì)土壤呼吸存在閾值效應(yīng)[19]。

    隨氮沉降時(shí)間延長(zhǎng),綜合兩個(gè)生長(zhǎng)季內(nèi)土壤呼吸速率的變化情況,低氮處理均顯著促進(jìn)土壤呼吸速率,而中氮和高氮處理則由促進(jìn)轉(zhuǎn)為抑制效應(yīng),與孫濤[32]模擬氮沉降前2年的研究結(jié)果一致,表明氮沉降可能存在累積效應(yīng),長(zhǎng)時(shí)間持續(xù)高濃度的氮沉降處理會(huì)抑制土壤呼吸作用。

    3.2氮沉降處理下溫度和濕度對(duì)落葉松人工林土壤呼吸的影響

    土壤溫度變化也是影響土壤呼吸速率的主要因素[42]。本研究不同模擬氮沉降處理下土壤呼吸速率均與土壤溫度(5,10 cm)顯著相關(guān)(P<0.01),可以解釋土壤呼吸速率變化的絕大部分[14]。然而不同氮沉降處理下的土壤呼吸的溫度敏感性(Q10)存在差異,說(shuō)明不同水平的氮沉降處理改變了土壤呼吸的溫度敏感性。隨土壤溫度升高,低氮處理土壤微生物和土壤酶活性相應(yīng)增強(qiáng),微生物對(duì)氮的利用率提高,促進(jìn)土壤有機(jī)質(zhì)分解[43],從而表現(xiàn)出異養(yǎng)呼吸強(qiáng)度提高,間接提高土壤呼吸速率;而中氮和高氮處理的閾值效應(yīng)導(dǎo)致土壤微生物和酶活性相應(yīng)降低,從而表現(xiàn)出異養(yǎng)呼吸強(qiáng)度降低,間接降低土壤呼吸速率(圖2)。本研究10 cm處Q10值均高于5 cm,表明10 cm土壤溫度變化可以更好的解釋土壤呼吸速率的變化趨勢(shì)。

    一般地,土壤濕度也是影響土壤呼吸速率的主要因素之一。當(dāng)土壤過(guò)于干旱或者土壤濕度過(guò)大甚至發(fā)生洪澇災(zāi)害時(shí),都會(huì)造成土壤呼吸速率急劇降低[44]。以往研究發(fā)現(xiàn),當(dāng)土壤濕度低于15%[45]或超過(guò)60%[44]時(shí),水分成為限制植物根系生長(zhǎng)和微生物活動(dòng)的主要條件,土壤呼吸速率與土壤濕度顯著相關(guān)。本研究土壤呼吸與土壤濕度間無(wú)顯著相關(guān)性,原因可能是土壤濕度處于15%~55%,并未限制林分根系的生長(zhǎng)以及土壤微生物的活動(dòng),說(shuō)明當(dāng)土壤濕度處于合適范圍時(shí),氮沉降處理不能改變土壤呼吸和土壤濕度間的關(guān)系。

    研究結(jié)果表明:在中齡落葉松人工林及暗棕壤條件下:(1)短期低水平氮沉降提高土壤微生物的生物量及其活性,而高水平氮沉降則未表現(xiàn)出顯著影響;(2)短期氮沉降或連續(xù)低水平氮沉降加快土壤碳排放過(guò)程,促進(jìn)土壤和大氣的碳循環(huán)過(guò)程,高水平氮沉降則減緩?fù)寥捞寂欧胚^(guò)程,有利于短期內(nèi)土壤碳的固持。

    致謝感謝東北林業(yè)大學(xué)王昌亮、李雯和張書(shū)娜碩士在野外和實(shí)驗(yàn)室內(nèi)工作中給予的幫助。

    1.Dixon R K,Solomon A M,Brown S,et al.Carbon pools and flux of global forest ecosystems[J].Science,1994,263(5144):185-190.

    2.Post W M,Emanuel W R,Zinke P J,et al.Soil carbon pools and world life zones[J].Nature,1982,298(5870):156-159.

    3.Schlesinger W H,Andrews J A.Soil respiration and the global carbon cycle[J].Biogeochemistry,2000,48(1):7-20.

    4.張東秋,石培禮,張憲洲.土壤呼吸主要影響因素的研究進(jìn)展[J].地球科學(xué)進(jìn)展,2005,20(7):778-785.

    Zhang Dongqiu,Shi Peili,Zhang Xianzhou.Some advance in the main factors controlling soil respiration[J].Advances in Earth Science,2005,20(7):778-785.

    5.Wright R F,Rasmussen L.Introduction to the NITREX and EXMAN projects[J].Forest Ecology and Management,1998,101(1-3):1-7.

    6.Pretzsch H.Analysis and modeling of spatial stand structures.Methodological considerations based on mixed beech-larch stands in Lower Saxony[J].Forest Ecology and Management,1997,97(3):237-253.

    7.Lewandowski A U,Gadow K V.Ein heuristischer Ansatz zur reproduktion von Waldbest?nden[J].Allgemeine Forest und Jagdzeitung,1997,168(9):170-174.

    8.Franklin O,H?berg P,Ekblad A,et al.Pine forest floor carbon accumulation in response to N and P K additions:Bomb14C modelling and respiration studies[J].Ecosystems,2003,6(7):644-658.

    9.Maier C A,Kress L W.Soil CO2evolution and root respiration in 11 year-old loblolly pine(Pinustaeda) plantations as affected by moisture and nutrient availability[J].Canadian Journal of Forest Research,2000,30(3):347-359.

    10.Bowden R D,Davidson E,Savage K,et al.Chronic nitrogen additions reduce total soil respiration and microbial respiration in temperate forest soils at the Harvard Forest[J].Forest Ecology and Management,2004,196(1):43-56.

    11.Pregitzer K S,Zak D R,Maziasz J,et al.Interactive effects of atmospheric CO2and soil-N availability on fine roots ofPopulustremuloides[J].Ecological Applications,2000,10(1):18-33.

    12.涂利華,胡庭興,張健,等.模擬氮沉降對(duì)華西雨屏區(qū)慈竹林土壤活性有機(jī)碳庫(kù)和根生物量的影響[J].生態(tài)學(xué)報(bào),2010,30(9):2286-2294.

    Tu L H,Hu T X,Zhang J,et al.Effects of simulated nitrogen deposition on soil active carbon pool and root biomass inNeosinocalamusaffinisplantation,Rainy Area of West China[J].Acta Ecological Sinica,2010,30(9):2286-2294.

    13.Hanson P J,Edwards N T,Garten C T,et al.Separating root and soil microbial contributions to soil respiration:A review of methods and observations[J].Biogeochemistry,2000,48(1):115-146.

    14.涂利華,胡庭興,黃立華,等.華西雨屏區(qū)苦竹林土壤呼吸對(duì)模擬氮沉降的響應(yīng)[J].植物生態(tài)學(xué)報(bào),2009,33(4):728-738.

    Tu L H,Hu T X,Huang L H,et al.Response of soil respiration to simulated nitrogen deposition inPleioblastusamarusForest,Rainy Area of West China[J].Chinese Journal of Plant Ecology,2009,33(4):728-738.

    15.涂利華,戴洪忠,胡庭興,等.模擬氮沉降對(duì)華西雨屏區(qū)撐綠雜交竹林土壤呼吸的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(4):829-836.

    Tu L H,Dai H Z,Hu T X,et al.Effects of simulated nitrogen deposition on soil respiration in aBambusapervariabilis×Dendrocalamopsiplantation in Rainy Area of West China[J].Chinese Journal of Applied Ecology,2011,22(4):829-836.

    16.孫素琪,王玉杰,王云琦,等.縉云山常綠闊葉林土壤呼吸對(duì)模擬氮沉降的響應(yīng)[J].林業(yè)科學(xué),2014,50(1):1-8.

    Sun S Q,Wang Y J,Wang Y Q,et al.Responses of soil respiration to simulated nitrogen deposition in an evergreen broad-leaved forest in Jinyun Mountain[J].Scientia Silvae Sinicae,2014,50(1):1-8.

    17.黃玉梓,樊后保,李燕燕,等.氮沉降對(duì)杉木人工林土壤呼吸與土壤纖維素酶活性的影響[J].福建林學(xué)院學(xué)報(bào),2009,29(2):97-102.

    Huang Y Z,Fan H B,Li Y Y,et al.Impacts of nitrogen deposition on the respiration rate and soil cellulose enzyme activity in Chinese fir plantation[J].Journal of Fujian College of Forestry,2009,29(2):97-102.

    18.李凱,江洪,由美娜,等.模擬氮沉降對(duì)石櫟和苦櫧幼苗土壤呼吸的影響[J].生態(tài)學(xué)報(bào),2011,31(1):82-89.

    Li K,Jiang H,You M N,et al.Effect of simulated nitrogen deposition on the soil respiration ofLithocarpusglabraandCastanopsissclerophylla[J].Acta Ecologica Sinica,2011,31(1):82-89.

    19.李偉斌,金昌杰,井艷麗,等.長(zhǎng)白山闊葉紅松林土壤呼吸對(duì)氮沉降增加的響應(yīng)[J].東北林業(yè)大學(xué)學(xué)報(bào),2014,42(12):89-93.

    Li W B,Jing C J,Jing Y L,et al.Response of soil respiration to enhanced nitrogen deposition in broadleaved Korean Pine forest in Changbai Mountains[J].Journal of Northeast Forestry University,2014,42(12):89-93.

    20.莫江明,方云霆,徐國(guó)良,等.鼎湖山苗圃和主要森林土壤CO2排放和CH4吸收對(duì)模擬N沉降的短期響應(yīng)[J].生態(tài)學(xué)報(bào),2005,25(4):682-690.

    Mo J M,Fang Y T,Xu G L,et al.The short-term responses of soil CO2emmission and CH4uptake to simulated N deposition in nursery and forests of Dinghushan in subtropical China[J].Acta Ecologica Sinica,2005,25(4):682-690.

    21.Du E Z,Liu X J.High rates of wet nitrogen deposition in China:a synthesis[M].//Sutton M A,Mason K E,Sheppard L J,et al.Nitrogen deposition,critical loads and biodiversity.Netherlands:Springer,2014:49-56.

    22.Lü C Q,Tian H Q.Spatial and temporal patterns of nitrogen deposition in China:Synthesis of observational data[M].Journal of Geophysical Research:Atmospheres,2007,112 (D22):D22S05.

    23.賈淑霞,王政權(quán),梅莉,等.施肥對(duì)落葉松和水曲柳人工林土壤呼吸的影響[J].植物生態(tài)學(xué)報(bào),2007,31(3):372-379.

    Jia S X,Wang Z Q,Mei L,et al.Effect of nitrogen fertilization on soil resppitation in Larix Gmelinii and Fraxinus Mandshurica plantation in China[J].Chinese Journal of Plant Ecology,2007,31(3):372-379.

    24.陳喜全,歐陽(yáng)華.帽兒山實(shí)驗(yàn)林場(chǎng)的土壤[J].東北林業(yè)大學(xué)學(xué)報(bào),1986,16(S1):79-84.

    Chen X H,OuYang H.Soil in Mao’er shan Experimental Forest Farm[J].Journal of Northeast Forestry University,1986,14:79-84.

    25.莫江明,薛璟花,方運(yùn)霆.鼎湖山主要森林植物凋落物分解及其對(duì)N沉降的響應(yīng)[J].生態(tài)學(xué)報(bào),2004,24(7):1413-1420.

    Mo J M,Xue J H,Fang Y T.Litter decomposition and its responses to simulated N deposition for the major plants of Dinghushan forests in subtropical China[J].Acta Ecologica Sinica,2004,24(7):1413-1420.

    26.Magill A H,Aber J D,Berntson G M,et al.Long-term nitrogen additions and nitrogen saturation in two temperate forests[J].Ecosystems,2000,3(3):238-253.

    27.Bond-lamberty B,Wang C K,Gower S T.Contribution of root respiration to soil surface CO2flux in a boreal black spruce fire chronosequence[J].Tree Physiology,2004,24(12):1387-1395.

    28.Livingston G P,Hutchinson G L.Enclosure-based measurement of trace gas exchange:applications and sources of Error[M].//Matson P A,Harris R C.Biogenic Trace Gases:Measuring Emissions from Soil and Water.Oxford:Blackwell Scientific Publications,1995:14-51.

    29.易志剛.離體根呼吸研究及其在原位根呼吸測(cè)量中的應(yīng)用[J].福建林學(xué)院學(xué)報(bào),2008,28(3):212-215.

    Yi Z G.Study on excised root respiration and its application to in situ measurement of root respiration[J].Journal of Fujian College of Forestry,2008,28(3):212-215.

    30.George K,Norby R J,Hamilton J G,et al.Fine-root respiration in a loblolly pine and sweetgum forest growing in elevated CO2[J].New Phytologist,2003,160(3):511-522.

    31.Jenkinson D S,Powlson D S.The effects of biocidal treatments on metabolism in soil-V:A method for measuring soil biomass[J].Soil Biology and Biochemistry,1976,8(3):209-213.

    32.孫濤.模擬氮沉降對(duì)東北地區(qū)興安落葉松人工林生態(tài)系統(tǒng)呼吸主要組分影響研究[D].哈爾濱:東北林業(yè)大學(xué),2014.

    Sun T.Effects of simulated nitrogen deposition on the main components of ecosystem respiration inLarixgmeliniiplantation in Northeastern China[D].Northeast Forestry University,2014.

    33.Zhou L Y,Zhou X H,Zhang B C,et al.Different responses of soil respiration and its components to nitrogen addition among biomes:A meta-analysis[J].Global Change Biology,2014,20(7):2332-2343.

    34.Aber J D,Mcdowell W,Nadelhoffer K J,et al.Nitrogen saturation in Temperate forest ecosystems[J].Bioscience,1998,48(11):921-934.

    35.涂利華,胡庭興,張健,等.模擬氮沉降對(duì)華西雨屏區(qū)苦竹林細(xì)根特性和土壤呼吸的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2010,21(10):2472-2478.

    Tu L H,Hu T X,Zhang J,et al.Effects of simulated nitrogen deposition on the fine root characteristics and soil respiration in aPleioblastusamarusplantation in Rainy Area of West China[J].Chinese Journal of Applied Ecology,2010,21(10):2472-2478.

    36.Chapin F S III,Matson P A,Mooney H A.Principles of terrestrial ecosystem ecology[M].New York:Springer-Verlag,2002:136-137.

    37.Burton A J,Pregitzer K,Ruess R,et al.Root respiration in North American forests:effects of nitrogen concentration and temperature across biomes[J].Oecologia,2002,131(4):559-568.

    38.Vitousek P M,H?ttenschwiler S,Olander L,et al.Nitrogen and nature[J].AMBIO,2002,31(2):97-101.

    39.Melillo J M,Aber J D,Muratore J F.Nitrogen and lignin control of hardwood leaf litter decomposition dynamics[J].Ecology,1982,63(3):621-626.

    40.Singh J S,Gupta S R.Plant decomposition and soil respiration in terrestrial ecosystems[J].The Botanical Review,1997,43(4):449-528.

    41.袁穎紅,樊后保,劉文飛,等.模擬氮沉降對(duì)杉木人工林土壤可溶性有機(jī)碳和微生物量碳的影響[J].水土保持學(xué)報(bào),2012,26(2):138-143.

    Yuan Y H,Fan H B,Liu W F,et al.Effects simulated nitrogen deposition deposition on soil dissolved organic carbon and microbial biomass carbon in a Chinese Fir plantation[J].Journal of Soil and Water Conservation,2012,26(2):138-143.

    42.Chapin F S III,Matson P A,Mooney H A.Principles of terrestrial ecosystem ecology[J].New York:Springer-Verlag,2002:159-166.

    43.楊萬(wàn)勤,王開(kāi)運(yùn).森林土壤酶的研究進(jìn)展[J].林業(yè)科學(xué),2004,40(2):152-159.

    Yang W Q,Wang K Y.Advances in forest soil enzymology[J].Scientia Silvae Sinicae,2004,40(2):152-159.

    44.Linn D M,Doran J W.Effect of water-filled pore space on carbon dioxide and nitrous oxide production in tilled and nontilled soils[J].Soil Science Society of America Journal,1984,48(6):1267-1272.

    45.鄧琦,周?chē)?guó)逸,劉菊秀,等.CO2濃度倍增、高氮沉降和高降雨對(duì)南亞熱帶人工模擬森林生態(tài)系統(tǒng)土壤呼吸的影響[J].植物生態(tài)學(xué)報(bào),2009,33(6):1023-1033.

    Deng Q,Zhou G Y,Liu X J,et al.Effects of CO2enrichment,high nitrogen deposition and high precipitation on a model forest ecosystem southern China[J].Chinese Journal of Plant Ecology,2009,33(6):1023-1033.

    《植物研究》榮獲“2015年中國(guó)高校技術(shù)類優(yōu)秀期刊”獎(jiǎng)

    National “twelve five” scientific and technological project.

    introduction:ZHANG Jiao(1990—),female,Master,Forestry college of Northeast Forestry University.

    date:2016-01-05

    EffectofSimulatedNDepositiononSoilRespirationinaLarixolgensisPlantationForest

    ZHANG Jiao HAO Long-Fei WANG Qing-Cheng*FU Jiao-Jiao ZHU Kai-Yue

    (Northeast Forestry University,Harbin 150040)

    A two-year simulated N deposition experiment was conducted to explore the effects of different levels of nitrogen(N)deposition on soil respiration, with four different N supply levels, including control(CK, 0 g·m-2·a-1N), low(N1, 5 g·m-2·a-1N), intermediate(N2, 10 g·m-2·a-1N)and high level(N3, 15 g·m-2·a-1N), in a 26-year-oldLarixolgensisplantation forest in Mao’ershan experimental forest farm of northeast forestry university. Average annual soil respiration rate under different treatments was significantly promoted in the first year(P<0.05); Average annual soil respiration rate and heterotrophic respiration rate were significantly decreased under N2 and N3 treatments(P<0.05), respectively, and average annual soil respiration rate was significantly increased(P<0.05), heterotrophic respiration rate had no difference compared with control(P>0.05)under N1 treatment, in the second year. Soil microbial biomass carbon was significantly increased under N1 treatment(P<0.05), while N2 and N3 had no difference compared with control(P>0.05); Significantly positive exponential relationship was detected between soil respiration and soil temperature in 5 and 10 cm(P<0.01). The temperature sensitivity of soil respiration(Q10)in both soil layer was increased under N1 treatment and decreased under N2 and N3 treatments compared with control, respectively. Soil respiration rate and its temperature sensitivity were altered under different N deposition levels, suggesting that the progress of soil carbon emission was accelerated under low levels of N deposition and slowed down under relatively high levels of N deposition in a short period.

    Larixolgensis;plantation forest;N deposition;soil respiration;root respiration;heterotrophic respiration

    國(guó)家“十二五”科技支撐項(xiàng)目(2011BAD37B02)

    張嬌(1990—),女,碩士研究生,主要從事森林生態(tài)學(xué)研究。

    * 通信作者:E-mail:wqcnefu@163.com

    2016-01-05

    * Corresponding author:E-mail:wqcnefu@163.com

    S718.55;S714.5

    A

    10.7525/j.issn.1673-5102.2016.04.016

    猜你喜歡
    中氮異養(yǎng)落葉松
    山西落葉松雜交良種逾10萬(wàn)畝
    異養(yǎng)同化降解氯代烴的研究現(xiàn)狀、微生物代謝特性及展望
    長(zhǎng)白落葉松離體再生體系的建立
    神奇的落葉松提取物
    清水江水體和沉積物中氮、磷的分布及變化趨勢(shì)
    蒸餾滴定法測(cè)定高溫合金中氮的含量
    西藏濕地生態(tài)系統(tǒng)中氮循環(huán)微生物數(shù)量和多樣性研究
    西藏科技(2015年4期)2015-09-26 12:12:59
    馬尾松種子生長(zhǎng)與幼苗成長(zhǎng)過(guò)程的養(yǎng)分動(dòng)態(tài)研究
    長(zhǎng)期不同施肥對(duì)土壤各粒級(jí)組分中氮含量及分配比例的影響
    異養(yǎng)小球藻的篩選鑒定及其特性研究
    禁无遮挡网站| 国产老妇伦熟女老妇高清| 欧美又色又爽又黄视频| 欧美三级亚洲精品| 人妻夜夜爽99麻豆av| 边亲边吃奶的免费视频| 亚洲自拍偷在线| 少妇丰满av| 99久久精品热视频| 国产一区亚洲一区在线观看| 国产男人的电影天堂91| 99在线人妻在线中文字幕| 久久久久久久久久成人| 成年版毛片免费区| 日韩成人伦理影院| 又爽又黄a免费视频| 晚上一个人看的免费电影| 女的被弄到高潮叫床怎么办| 秋霞在线观看毛片| 亚洲欧美一区二区三区国产| 亚洲国产精品成人综合色| 欧美一区二区亚洲| 人人妻人人澡欧美一区二区| 亚洲av中文av极速乱| 又粗又爽又猛毛片免费看| 波野结衣二区三区在线| 一个人看的www免费观看视频| 白带黄色成豆腐渣| 在线观看av片永久免费下载| 欧美精品国产亚洲| 美女cb高潮喷水在线观看| 亚洲一区高清亚洲精品| 男女那种视频在线观看| 日本黄色视频三级网站网址| 三级国产精品欧美在线观看| av在线蜜桃| av天堂中文字幕网| 日韩中字成人| 色噜噜av男人的天堂激情| 高清av免费在线| 亚洲av熟女| 青春草国产在线视频| 日本一本二区三区精品| 插逼视频在线观看| 精品一区二区免费观看| 一区二区三区四区激情视频| 又爽又黄无遮挡网站| 亚洲图色成人| 成年av动漫网址| 一本久久精品| 成人av在线播放网站| 男人狂女人下面高潮的视频| 91狼人影院| 麻豆久久精品国产亚洲av| 国产精品av视频在线免费观看| or卡值多少钱| 久久人人爽人人片av| 乱码一卡2卡4卡精品| 伊人久久精品亚洲午夜| 精华霜和精华液先用哪个| 久久久久久国产a免费观看| 午夜福利在线观看吧| 色视频www国产| 日韩成人伦理影院| 亚洲欧洲国产日韩| 毛片女人毛片| 亚洲精华国产精华液的使用体验| 97人妻精品一区二区三区麻豆| 18禁在线播放成人免费| 一卡2卡三卡四卡精品乱码亚洲| 国产免费男女视频| 久久精品国产自在天天线| av视频在线观看入口| 一卡2卡三卡四卡精品乱码亚洲| 99久久九九国产精品国产免费| 婷婷色综合大香蕉| 色吧在线观看| 欧美成人一区二区免费高清观看| 国产单亲对白刺激| 毛片女人毛片| 精品国内亚洲2022精品成人| 国产精品综合久久久久久久免费| 国产成人91sexporn| 最近中文字幕高清免费大全6| 乱系列少妇在线播放| 成年免费大片在线观看| 能在线免费观看的黄片| 人体艺术视频欧美日本| 联通29元200g的流量卡| 男女那种视频在线观看| av在线老鸭窝| 午夜精品一区二区三区免费看| 精华霜和精华液先用哪个| 精品人妻偷拍中文字幕| 久久热精品热| 女的被弄到高潮叫床怎么办| 免费不卡的大黄色大毛片视频在线观看 | 免费黄色在线免费观看| 日本av手机在线免费观看| 性插视频无遮挡在线免费观看| 亚洲av电影在线观看一区二区三区 | 亚洲怡红院男人天堂| 亚洲欧美精品专区久久| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 久久精品夜色国产| 国产精品永久免费网站| 午夜视频国产福利| 网址你懂的国产日韩在线| 国产精品久久久久久久电影| 欧美又色又爽又黄视频| 我的女老师完整版在线观看| 日本与韩国留学比较| 99久久中文字幕三级久久日本| 国产男人的电影天堂91| 噜噜噜噜噜久久久久久91| 一二三四中文在线观看免费高清| 国产精华一区二区三区| 91精品国产九色| 免费看av在线观看网站| 大话2 男鬼变身卡| 伊人久久精品亚洲午夜| 亚洲精华国产精华液的使用体验| 黄色日韩在线| 伦理电影大哥的女人| 国产精品乱码一区二三区的特点| 18+在线观看网站| 成人亚洲欧美一区二区av| 中文字幕av成人在线电影| 国产三级在线视频| 国产老妇女一区| 老司机影院成人| 男人和女人高潮做爰伦理| 性插视频无遮挡在线免费观看| 欧美zozozo另类| 99久久精品国产国产毛片| 麻豆av噜噜一区二区三区| 内地一区二区视频在线| 天堂网av新在线| 高清视频免费观看一区二区 | 欧美极品一区二区三区四区| 久久韩国三级中文字幕| av卡一久久| 91在线精品国自产拍蜜月| 一区二区三区四区激情视频| 亚洲国产精品成人久久小说| 99热这里只有精品一区| 久久久精品大字幕| 精品酒店卫生间| 国产乱人视频| 日本午夜av视频| 亚洲综合色惰| 国产成人福利小说| 欧美精品国产亚洲| 精品免费久久久久久久清纯| 亚洲四区av| 蜜桃久久精品国产亚洲av| 欧美日韩国产亚洲二区| 国产免费又黄又爽又色| 熟妇人妻久久中文字幕3abv| 99久国产av精品| 免费看美女性在线毛片视频| 天天躁日日操中文字幕| 国产成人aa在线观看| 国产精品蜜桃在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久亚洲国产成人精品v| 熟女电影av网| 久久99精品国语久久久| 九九久久精品国产亚洲av麻豆| 七月丁香在线播放| 深爱激情五月婷婷| 成人av在线播放网站| 十八禁国产超污无遮挡网站| 熟妇人妻久久中文字幕3abv| 91午夜精品亚洲一区二区三区| 久久婷婷人人爽人人干人人爱| 校园人妻丝袜中文字幕| 久久久久性生活片| 亚洲三级黄色毛片| 久久久欧美国产精品| 国产伦在线观看视频一区| 亚洲国产精品合色在线| 我的老师免费观看完整版| 色吧在线观看| 久久综合国产亚洲精品| 国产精品久久久久久av不卡| 国产人妻一区二区三区在| 亚洲精品乱码久久久v下载方式| 国产成人freesex在线| 国产高潮美女av| 高清日韩中文字幕在线| 欧美性猛交╳xxx乱大交人| 狂野欧美白嫩少妇大欣赏| 直男gayav资源| 伦理电影大哥的女人| 26uuu在线亚洲综合色| 天美传媒精品一区二区| 男人狂女人下面高潮的视频| 亚洲无线观看免费| 又粗又爽又猛毛片免费看| 麻豆国产97在线/欧美| 在现免费观看毛片| 国产精品伦人一区二区| 欧美高清性xxxxhd video| 精品国内亚洲2022精品成人| 伦精品一区二区三区| 91久久精品电影网| 亚洲av成人av| 少妇裸体淫交视频免费看高清| 麻豆成人午夜福利视频| 久久人妻av系列| av天堂中文字幕网| 久久99热这里只有精品18| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| 国产精品人妻久久久久久| 日韩一本色道免费dvd| 成年女人永久免费观看视频| 老司机影院成人| 国产日韩欧美在线精品| 精品熟女少妇av免费看| 精品国产三级普通话版| 一区二区三区四区激情视频| 如何舔出高潮| 亚洲在久久综合| 午夜日本视频在线| 亚洲精品成人久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 久久久亚洲精品成人影院| 国内揄拍国产精品人妻在线| 最近中文字幕高清免费大全6| 久久精品夜色国产| 美女被艹到高潮喷水动态| 一个人看的www免费观看视频| 日韩中字成人| 最近中文字幕高清免费大全6| 国产乱人偷精品视频| 国产熟女欧美一区二区| av免费在线看不卡| 最后的刺客免费高清国语| 日韩欧美国产在线观看| 插逼视频在线观看| 精品一区二区免费观看| 一本久久精品| 久久人妻av系列| 亚洲成人av在线免费| av又黄又爽大尺度在线免费看 | av卡一久久| 欧美成人免费av一区二区三区| 内射极品少妇av片p| 日本免费一区二区三区高清不卡| 亚洲av免费在线观看| 午夜福利在线观看吧| 欧美成人免费av一区二区三区| 亚洲最大成人中文| 成人三级黄色视频| 女人十人毛片免费观看3o分钟| 欧美3d第一页| 国产探花极品一区二区| www日本黄色视频网| 亚洲va在线va天堂va国产| 亚洲国产成人一精品久久久| 国产又色又爽无遮挡免| 精品一区二区三区人妻视频| 久久亚洲精品不卡| 欧美zozozo另类| 国产一级毛片在线| 国产成人午夜福利电影在线观看| 两个人的视频大全免费| 麻豆国产97在线/欧美| 亚洲国产色片| 日韩一区二区视频免费看| av播播在线观看一区| 国产高清有码在线观看视频| 亚洲av不卡在线观看| 国产精品.久久久| 日韩av在线免费看完整版不卡| 国产精品一及| 午夜福利在线观看吧| 老女人水多毛片| 床上黄色一级片| 毛片女人毛片| 国产欧美日韩精品一区二区| 日韩欧美精品v在线| 国产人妻一区二区三区在| 伦精品一区二区三区| 亚洲欧洲国产日韩| 久久精品久久精品一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲丝袜综合中文字幕| 国产伦一二天堂av在线观看| 日本五十路高清| 色综合站精品国产| 日韩一区二区三区影片| 国产精品不卡视频一区二区| 亚洲乱码一区二区免费版| 国产高清有码在线观看视频| 亚洲精品一区蜜桃| 最近中文字幕2019免费版| 在线观看一区二区三区| 精品久久久久久久末码| 中文在线观看免费www的网站| 国产乱人视频| 国产一区二区亚洲精品在线观看| 啦啦啦啦在线视频资源| 成人亚洲欧美一区二区av| 国产精品三级大全| 如何舔出高潮| 亚洲一区高清亚洲精品| 久久久久久久久久成人| 少妇丰满av| 国产视频内射| 午夜福利网站1000一区二区三区| 有码 亚洲区| 欧美xxxx性猛交bbbb| av天堂中文字幕网| 精品久久国产蜜桃| 中文字幕制服av| 熟妇人妻久久中文字幕3abv| 一边亲一边摸免费视频| 久久久国产成人精品二区| 干丝袜人妻中文字幕| 麻豆精品久久久久久蜜桃| 国产av码专区亚洲av| 最近中文字幕高清免费大全6| 成人综合一区亚洲| 日本色播在线视频| 国产在视频线精品| 91精品一卡2卡3卡4卡| 亚洲自偷自拍三级| 乱系列少妇在线播放| 美女脱内裤让男人舔精品视频| 噜噜噜噜噜久久久久久91| 国产综合懂色| 精品人妻视频免费看| 日本欧美国产在线视频| 特级一级黄色大片| 国产成人精品婷婷| 极品教师在线视频| 久久热精品热| 亚洲国产最新在线播放| 国产成人午夜福利电影在线观看| 七月丁香在线播放| 69人妻影院| 欧美日本视频| 我的老师免费观看完整版| 97热精品久久久久久| 国产免费男女视频| 国产成人福利小说| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 六月丁香七月| 亚洲精华国产精华液的使用体验| 免费av毛片视频| 成年免费大片在线观看| 免费观看在线日韩| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 国产成人freesex在线| 国产av在哪里看| 国产精品日韩av在线免费观看| 日韩大片免费观看网站 | 最近中文字幕高清免费大全6| 欧美成人免费av一区二区三区| 夫妻性生交免费视频一级片| 国产精品一区二区三区四区久久| 精品欧美国产一区二区三| 97超碰精品成人国产| eeuss影院久久| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 日韩亚洲欧美综合| 亚洲在久久综合| 亚洲精品久久久久久婷婷小说 | 亚洲欧美日韩高清专用| 国产在线男女| 亚洲怡红院男人天堂| 日韩 亚洲 欧美在线| 九草在线视频观看| 99久久中文字幕三级久久日本| 在现免费观看毛片| 久久这里有精品视频免费| 日韩欧美 国产精品| 久久久精品94久久精品| 国产精品美女特级片免费视频播放器| 国产在线男女| 日韩制服骚丝袜av| 日韩 亚洲 欧美在线| 少妇熟女aⅴ在线视频| 国产成年人精品一区二区| 免费在线观看成人毛片| av福利片在线观看| 中文字幕熟女人妻在线| 国产乱人偷精品视频| 村上凉子中文字幕在线| 国产人妻一区二区三区在| 国产爱豆传媒在线观看| 韩国高清视频一区二区三区| 黄片wwwwww| 三级国产精品片| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件| 久久精品国产鲁丝片午夜精品| 一区二区三区高清视频在线| 九九热线精品视视频播放| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 欧美精品国产亚洲| 久久精品夜色国产| 在线a可以看的网站| 亚洲av电影在线观看一区二区三区 | 爱豆传媒免费全集在线观看| 又爽又黄无遮挡网站| 国产女主播在线喷水免费视频网站 | 欧美激情国产日韩精品一区| av免费观看日本| 99久久精品国产国产毛片| 国产伦在线观看视频一区| 少妇的逼好多水| 老司机福利观看| 波多野结衣高清无吗| 久久久久久国产a免费观看| 国产精品麻豆人妻色哟哟久久 | 最近手机中文字幕大全| 两个人的视频大全免费| 99热这里只有是精品在线观看| av播播在线观看一区| 欧美xxxx性猛交bbbb| 亚洲成人av在线免费| 观看美女的网站| 乱码一卡2卡4卡精品| 国产 一区精品| 亚洲欧美清纯卡通| 国产在线男女| 亚洲伊人久久精品综合 | 国产高清国产精品国产三级 | 女的被弄到高潮叫床怎么办| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 国产亚洲精品av在线| 亚洲精品日韩av片在线观看| 精品欧美国产一区二区三| 国产免费男女视频| 美女黄网站色视频| 热99re8久久精品国产| 三级男女做爰猛烈吃奶摸视频| 一个人免费在线观看电影| 久久人妻av系列| 少妇丰满av| 久久99热6这里只有精品| av在线蜜桃| 春色校园在线视频观看| 欧美精品一区二区大全| 超碰97精品在线观看| 91午夜精品亚洲一区二区三区| 1024手机看黄色片| 亚洲精华国产精华液的使用体验| 免费av毛片视频| 国产亚洲一区二区精品| 99视频精品全部免费 在线| 少妇裸体淫交视频免费看高清| 婷婷六月久久综合丁香| 美女内射精品一级片tv| 欧美人与善性xxx| 成年女人永久免费观看视频| kizo精华| 久久6这里有精品| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区国产| 人人妻人人看人人澡| 国产免费福利视频在线观看| 日韩一区二区三区影片| 亚洲欧洲国产日韩| 中文字幕免费在线视频6| av国产免费在线观看| 桃色一区二区三区在线观看| 国产免费福利视频在线观看| av视频在线观看入口| 夜夜爽夜夜爽视频| 精品久久久久久电影网 | 亚洲中文字幕日韩| 日韩av在线大香蕉| 91精品一卡2卡3卡4卡| 中文精品一卡2卡3卡4更新| 免费黄网站久久成人精品| 两个人的视频大全免费| 99久久精品国产国产毛片| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 国产精品国产三级国产av玫瑰| 小蜜桃在线观看免费完整版高清| 国产免费视频播放在线视频 | 日本午夜av视频| 免费看日本二区| 欧美日本亚洲视频在线播放| 国产精华一区二区三区| www.av在线官网国产| 亚洲av电影在线观看一区二区三区 | .国产精品久久| 久久精品久久精品一区二区三区| 看黄色毛片网站| 精品久久国产蜜桃| 国产精品电影一区二区三区| 亚洲精品色激情综合| 成人亚洲精品av一区二区| 亚洲精品乱码久久久久久按摩| 欧美xxxx性猛交bbbb| 男的添女的下面高潮视频| 久久草成人影院| 亚洲精品乱久久久久久| 99久久精品一区二区三区| 五月玫瑰六月丁香| 99热这里只有是精品50| av在线观看视频网站免费| 老司机影院成人| 精品国产露脸久久av麻豆 | 一区二区三区高清视频在线| 国产精华一区二区三区| 九九在线视频观看精品| 午夜老司机福利剧场| 女人久久www免费人成看片 | 国产精品人妻久久久久久| 乱人视频在线观看| 人体艺术视频欧美日本| a级毛色黄片| 久久久久久久久大av| kizo精华| 午夜老司机福利剧场| 精品熟女少妇av免费看| 大香蕉97超碰在线| 国产女主播在线喷水免费视频网站 | 精华霜和精华液先用哪个| 久久99热6这里只有精品| 永久免费av网站大全| 国模一区二区三区四区视频| 国产精品国产高清国产av| 欧美激情在线99| 亚洲国产精品合色在线| 最近中文字幕高清免费大全6| 亚洲欧美日韩东京热| 日韩av在线大香蕉| 欧美三级亚洲精品| 男人狂女人下面高潮的视频| 亚洲av福利一区| 国内精品一区二区在线观看| 人妻少妇偷人精品九色| 国内揄拍国产精品人妻在线| 最近最新中文字幕大全电影3| 欧美精品国产亚洲| a级一级毛片免费在线观看| 最近视频中文字幕2019在线8| 精品酒店卫生间| 国产精品熟女久久久久浪| 国产精品国产三级专区第一集| 国产视频内射| 久久国产乱子免费精品| 人妻系列 视频| 久久欧美精品欧美久久欧美| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 精品一区二区三区视频在线| 一区二区三区高清视频在线| 国产真实乱freesex| 国内精品美女久久久久久| 最近2019中文字幕mv第一页| 亚洲国产高清在线一区二区三| 美女cb高潮喷水在线观看| 欧美三级亚洲精品| 国产精品久久久久久精品电影小说 | 午夜a级毛片| 久久久欧美国产精品| 性色avwww在线观看| 国产真实伦视频高清在线观看| 啦啦啦观看免费观看视频高清| 色播亚洲综合网| 日日啪夜夜撸| 老女人水多毛片| 免费观看性生交大片5| 天天躁夜夜躁狠狠久久av| 又爽又黄a免费视频| 少妇裸体淫交视频免费看高清| 1000部很黄的大片| 别揉我奶头 嗯啊视频| 美女xxoo啪啪120秒动态图| 寂寞人妻少妇视频99o| 久久久久久久国产电影| 国产男人的电影天堂91| 久久久欧美国产精品| 国产高清有码在线观看视频| 亚洲成人av在线免费| 亚洲精品影视一区二区三区av| 九九在线视频观看精品| 好男人在线观看高清免费视频| 国产女主播在线喷水免费视频网站 | 亚洲欧美精品专区久久| 国产精品久久久久久久久免| 3wmmmm亚洲av在线观看| 内地一区二区视频在线| av免费在线看不卡| 一级毛片aaaaaa免费看小| 国产精品女同一区二区软件| 1000部很黄的大片| 国产在线一区二区三区精 | 久久这里有精品视频免费| 能在线免费看毛片的网站| 少妇人妻一区二区三区视频| 亚洲av中文字字幕乱码综合| 午夜免费激情av| 久久精品人妻少妇| 如何舔出高潮| 免费不卡的大黄色大毛片视频在线观看 |