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

    應(yīng)用高頻觀測探討不同森林經(jīng)營方式下礦質(zhì)土壤呼吸的晝夜動(dòng)態(tài)特征

    2017-02-08 03:05:47劉小飛周嘉聰蘇瑞蘭林偉盛熊德成陳岳民楊玉盛
    生態(tài)學(xué)報(bào) 2017年1期
    關(guān)鍵詞:經(jīng)營方式礦質(zhì)土壤溫度

    鄭 永, 劉小飛,*, 鄭 蔚, 周嘉聰, 蘇瑞蘭, 林偉盛, 熊德成, 胥 超, 陳岳民, 楊玉盛

    1 福建師范大學(xué)地理科學(xué)學(xué)院, 福州 350007 2 福建師范大學(xué)濕潤亞熱帶山地生態(tài)國家重點(diǎn)實(shí)驗(yàn)室培育基地, 福州 350007 3 三明市三元區(qū)林業(yè)局, 三明 365000

    應(yīng)用高頻觀測探討不同森林經(jīng)營方式下礦質(zhì)土壤呼吸的晝夜動(dòng)態(tài)特征

    鄭 永1,2, 劉小飛1,2,*, 鄭 蔚1,2, 周嘉聰1,2, 蘇瑞蘭3, 林偉盛1,2, 熊德成1,2, 胥 超1,2, 陳岳民1,2, 楊玉盛1,2

    1 福建師范大學(xué)地理科學(xué)學(xué)院, 福州 350007 2 福建師范大學(xué)濕潤亞熱帶山地生態(tài)國家重點(diǎn)實(shí)驗(yàn)室培育基地, 福州 350007 3 三明市三元區(qū)林業(yè)局, 三明 365000

    礦質(zhì)土壤呼吸是森林生態(tài)系統(tǒng)土壤碳庫損失的重要途徑之一,也是森林生態(tài)系統(tǒng)碳(C)平衡估算中的關(guān)鍵因子。了解礦質(zhì)土壤呼吸在不同時(shí)間尺度上的變化,對理解森林生態(tài)系統(tǒng)C循環(huán)應(yīng)對全球變化的響應(yīng)至關(guān)重要,而高頻觀測是探討礦質(zhì)土壤呼吸在不同時(shí)間尺度變化的重要手段之一。通過高頻自動(dòng)觀測系統(tǒng)與Li- 8100土壤CO2通量測量系統(tǒng),對福建省三明市陳大鎮(zhèn)國有林場的米櫧(Castanopsiscarlesii)次生林在不同森林經(jīng)營方式下(CK對照,RR皆伐,RB火燒)的礦質(zhì)土壤呼吸與土壤溫度和含水量的晝夜動(dòng)態(tài)進(jìn)行分析,并比較2種采樣策略下礦質(zhì)土壤呼吸的年、日均通量差異。結(jié)果表明:1)不同森林經(jīng)營方式的礦質(zhì)土壤呼吸與土壤溫度和土壤含水量均存在著明顯的季節(jié)動(dòng)態(tài),礦質(zhì)土壤呼吸速率年均值表現(xiàn)為CK(2.18 μmol m-2s-1)>RB(1.93 μmol m-2s-1)>RR(1.89 μmol m-2s-1)。 2)在不同森林經(jīng)營方式下,采用手動(dòng)觀測的礦質(zhì)土壤呼吸年平均日通量顯著低于高頻觀測結(jié)果,而采用高頻觀測09:00—11:00時(shí)間段內(nèi)觀測數(shù)據(jù)計(jì)算日通量與高頻自動(dòng)觀測系統(tǒng)全天(24h)結(jié)果無顯著差異;3)不同森林經(jīng)營方式下的林地,土壤水熱條件的變化是影響礦質(zhì)土壤呼吸的重要因素之一。雙因子模型擬合結(jié)果表明,土壤溫度和含水量共同解釋了CK、RR和RB礦質(zhì)土壤呼吸速率的年變化的96.8%,62.8%,95.4%,擬合結(jié)果明顯優(yōu)于以溫度為單因子的指數(shù)模型。因此,未來氣候變化背景下,為準(zhǔn)確評估和預(yù)測不同森林經(jīng)營方式對土壤與大氣間碳通量交換的影響,采用高頻自動(dòng)觀測技術(shù)觀測礦質(zhì)土壤呼吸,將有利于提高碳通量估算精度。

    不同森林經(jīng)營方式;高頻自動(dòng)觀測;礦質(zhì)土壤呼吸;晝夜動(dòng)態(tài);碳通量

    土壤是陸地生態(tài)系統(tǒng)中最大的碳庫,據(jù)估計(jì),全球土壤有機(jī)碳庫約1500 Pg (1 Pg = 1015g),約是大氣碳庫的2倍,生物碳庫的3倍[1]。土壤呼吸是土壤碳庫和大氣碳庫交換的主要方式之一,其觀測手段主要為手動(dòng)觀測和高頻自動(dòng)觀測,由于土壤呼吸的時(shí)空異質(zhì)性,現(xiàn)今對土壤呼吸通量的評估仍存在著許多不確定性[2]。如:Savage等[3]在美國哈佛森林發(fā)現(xiàn)手動(dòng)觀測低估土壤呼吸日通量的13%,而Wang等[4]在加拿大農(nóng)田生態(tài)系統(tǒng)發(fā)現(xiàn)手動(dòng)觀測值相較于高頻觀測高估了土壤呼吸日通量的17%,手動(dòng)觀測的局限性在于不能進(jìn)行短時(shí)間內(nèi)頻繁的觀測,也無法準(zhǔn)確評估土壤呼吸對降雨、寒潮等特殊事件的響應(yīng),而高頻自動(dòng)觀測方法能準(zhǔn)確、迅速捕捉土壤呼吸對各種事件的響應(yīng)過程,對評估未來氣候變化下不同生態(tài)系統(tǒng)土壤呼吸的變化,進(jìn)而估算全球碳通量有著重要的意義[5-7]。

    土壤異養(yǎng)呼吸是森林生態(tài)系統(tǒng)土壤碳庫損失的重要途徑,也是森林生態(tài)系統(tǒng)碳平衡主要分量之一,與凈初級生產(chǎn)力共同決定森林碳匯的大小[8-9],熱帶和亞熱帶森林生態(tài)系統(tǒng)中異養(yǎng)呼吸占土壤總呼吸的60%—76%[10],而礦質(zhì)土壤呼吸占異養(yǎng)呼吸的30%—83%[11]。礦質(zhì)土壤呼吸主要來源于土壤微生物分解有機(jī)質(zhì)釋放的CO2,而土壤水分和溫度是影響微生物活性的關(guān)鍵因子,其表現(xiàn)在土壤溫度升高提高微生物活性,加快有機(jī)質(zhì)的分解速度,土壤含水量的增加則會(huì)有利于土壤中溶解性有機(jī)質(zhì)的擴(kuò)散,提高微生物可利用養(yǎng)分,從而促進(jìn)礦質(zhì)土壤呼吸[12-13]。

    有關(guān)森林土壤呼吸及其影響因素的研究已有大量報(bào)道[14-16],但主要集中在溫帶森林[17-18]。煉山是亞熱帶地區(qū)的主要造林方式,現(xiàn)今有關(guān)營林活動(dòng)(特別是皆伐和火燒等)對亞熱帶森林礦質(zhì)土壤呼吸的晝夜動(dòng)態(tài)的研究則鮮有報(bào)道。因而,本研究通過對亞熱帶森林不同經(jīng)營方式下礦質(zhì)土壤呼吸進(jìn)行手動(dòng)和高頻自動(dòng)觀測,探討不同經(jīng)營方式的中亞熱帶米櫧林礦質(zhì)土壤呼吸之間的晝夜動(dòng)態(tài)特征,以及手動(dòng)觀測與高頻自動(dòng)觀測對礦質(zhì)土壤呼吸年通量估算的差異性,提高未來氣候變化情景下不同森林經(jīng)營方式下土壤與大氣間CO2交換量的估算精度,進(jìn)一步改進(jìn)土壤呼吸經(jīng)驗(yàn)?zāi)P汀?/p>

    1 研究區(qū)概況

    試驗(yàn)地位于福建三明森林生態(tài)系統(tǒng)與全球變化野外觀測研究站-陳大觀測點(diǎn)(26°19′N,117°36′E),平均海拔330 m,平均坡度33 °,年均氣溫19.1 ℃,年均降水量1749 mm,屬中亞熱帶季風(fēng)性氣候。本區(qū)內(nèi)土壤厚度超過1 m,由黑云母花崗巖發(fā)育的紅壤。米櫧次生林前身為米櫧天然林,1978年天然林皆伐后天然更新演替形成,林分密度為2650株/hm2,平均樹高19.7 m,平均胸徑13.5 cm。主要樹種為米櫧(Castanopsiscarlesii)、閩粵栲(CastanopsisfissaRehd)、木荷(Schimasuperba)等,以米櫧為優(yōu)勢樹種。林下植被主要有毛冬青(IlexpubescensHook)、狗骨柴(Tricalysiadubia)、矩圓葉鼠刺(Iteachinensis)、沿海紫金牛(Ardisiapunctate)、狗脊蕨(Woodwardiajaponica)等。試驗(yàn)地土壤理化性質(zhì)見表1。

    表1 試驗(yàn)地土壤理化性質(zhì)

    不同小寫字母表示不同處理之間存在顯著差異(P<0.05),圖中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差(n=3);CK: Control; RR: Residual retention after clear cutting; RB: Residual burning after clear cutting

    2 研究方法

    2.1 樣地設(shè)置與礦質(zhì)土壤呼吸測定

    2011年11月在米櫧次生林(總面積為17.1 hm2)設(shè)置9塊20 m× 20 m標(biāo)準(zhǔn)地,分布于4個(gè)山坡,山坡間有30 m左右的坡溝相隔。樣地于2011年12月、2012 年3月相繼完成皆伐和火燒處理,各處理設(shè)置3個(gè)重復(fù),并按隨機(jī)區(qū)組排列,其處理分別為:(1)將原有次生林作為對照處理(CK);(2)在采伐跡地上采取人促更新,即皆伐(Residue Retention,RR);(3)對采伐跡地上火燒后人工造林,即火燒(Retention Burning,RB)。采用壕溝法去根: 2012 年4月,在各標(biāo)準(zhǔn)樣地上、中、下各布設(shè)1 m× 1 m的正方形區(qū)域,在其四周用鐵鍬挖壕溝,垂直挖深0.6—0.8 m (看不到根系),切斷根(不移走)后插入細(xì)密的尼龍網(wǎng)(100目)以阻止根向小區(qū)內(nèi)生長[19],貼地面剪除小區(qū)內(nèi)的植被,所有挖壕溝小區(qū)定期清理地表植被。

    礦質(zhì)土壤呼吸測定采用Li- 8100A (美國,基因公司) 土壤CO2通量全自動(dòng)測量系統(tǒng)測定,礦質(zhì)土壤呼吸測定頻率為每半個(gè)小時(shí)測定1次,全天24 h不間斷測定。在測定礦質(zhì)土壤呼吸速率的同時(shí),應(yīng)用Li- 8100A配套的土壤溫濕度傳感器測定土壤5 cm深處的溫度和體積含水量。觀測時(shí)間為2013年8月到2014年7月,并分別在2013年8月、10月、12月和2014年3月,選擇連續(xù)3天晴朗天氣,其觀測的數(shù)據(jù)分別代表2013年夏季、秋季、冬季和2014年春季礦質(zhì)土壤呼吸。

    于相同的觀測時(shí)間段內(nèi),使用Li- 8100開路式土壤碳通量測量系統(tǒng)在09:00—12:00手動(dòng)測定礦質(zhì)土壤呼吸[20],每個(gè)測定點(diǎn)測定時(shí)間為2 min,觀測時(shí)間在每個(gè)月月中(15日左右)和月末(30日左右)。同時(shí),采用時(shí)域反射儀(Model TDR300, Specturm Company, Aurora, USA)測量附近0—12 cm范圍內(nèi)的土壤體積含水量,每個(gè)測定點(diǎn)周圍測量3次,取平均值,采用長桿電子溫度探針(SK- 250WP, Sato Keir-Yoki, Kanda, Japan)測定附近土壤5 cm深處的土壤溫度。

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

    土壤呼吸日通量R(g C/m2)計(jì)算公式為:

    R=Rs×3600×24×12/1000000

    (1)

    運(yùn)用以下模型分析土壤呼吸速率與土壤溫度、土壤含水量的關(guān)系:

    Rs=aW+b

    (1)

    Rs=aebt

    (2)

    Rs=aebtWc

    (3)

    式中,Rs表示土壤呼吸速率,t表示土壤5 cm深處溫度,W表示0—12 cm土壤含水量,a、b、c為待定參數(shù)。

    土壤呼吸的溫度敏感性(Q10) 計(jì)算公式為:

    Q10=e10b

    (4)

    式中,b為(2)式中計(jì)算的常量。

    運(yùn)用SPSS 19.0進(jìn)行數(shù)據(jù)分析,并用Origin9.0進(jìn)行圖形繪制,采用重復(fù)測量方差(repeated-measures ANONA)檢驗(yàn)不同處理之間礦質(zhì)土壤呼吸速率、溫度、含水量的差異性,采用單因素方差(one-way ANONA)分析米櫧次生林不同森林經(jīng)營方式日平均碳排放通量的差異。

    3 結(jié)果與分析

    3.1 不同季節(jié)礦質(zhì)土壤呼吸的晝夜動(dòng)態(tài)特征

    各處理年均礦質(zhì)土壤呼吸表現(xiàn)出CK(2.18 μmol m-2s-1)>RB(1.93 μmol m-2s-1)>RR(1.89 μmol m-2s-1)的趨勢(圖1),但不同處理的礦質(zhì)土壤呼吸在不同季節(jié)的變化趨勢并不一致。CK和RB礦質(zhì)土壤呼吸速率表現(xiàn)為夏季高,冬季低;RR礦質(zhì)土壤呼吸速率則秋季最高,春季最低。春夏季,CK和RB的礦質(zhì)土壤呼吸顯著高于RR(P<0.05);秋冬季RR的礦質(zhì)土壤呼吸顯著高于CK和RB (P<0.05)。RR的礦質(zhì)土壤呼吸季節(jié)變化則與土壤溫度出現(xiàn)了季節(jié)滯后效應(yīng),CK和RB的土壤含水量季節(jié)變化均呈現(xiàn)春季高,冬季低的趨勢,RR的土壤含水量春季最高,夏季最低(圖1)。

    圖1 不同森林經(jīng)營方式土壤CO2速率、土壤溫度和含水量的季節(jié)動(dòng)態(tài)Fig.1 Different forest management of mineral soil CO2 rate, soil temperature and water content of seasonal dynamics不同大寫字母表示同一處理不同季節(jié)之間存在顯著差異(P<0.05),不同小寫字母表示同一季節(jié)不同處理之間存在顯著差異(P<0.05);圖中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差(n=3)

    礦質(zhì)土壤呼吸的晝夜變化較大,不同季節(jié)的晝夜動(dòng)態(tài)也不一致(圖2)。CK夏季礦質(zhì)土壤呼吸速率最高值出現(xiàn)在06:00和22:00左右,最低值出現(xiàn)在18:00左右,而秋、冬和春季節(jié)最高值出現(xiàn)于17:00左右,而最低值出現(xiàn)在08:00左右,其日變幅分別達(dá)到為32%,40%,49%和48%。RR夏、秋、冬和春季節(jié)的礦質(zhì)土壤呼吸速率最高值出現(xiàn)在于06:00和晚上23:00左右,最低值出現(xiàn)在18:00左右,其日變幅分別達(dá)到為54%,26%,77%和66%。RB夏、秋、冬和春季節(jié)的礦質(zhì)土壤呼吸速率最高值出現(xiàn)在13:00左右,最低值出現(xiàn)在06:00左右,其日變幅分別達(dá)到為45%,36%,74%和89%。

    圖2 不同森林經(jīng)營方式的礦質(zhì)土壤呼吸速率、土壤溫度和含水量的各季節(jié)的晝夜動(dòng)態(tài)變化Fig.2 Different forest management of mineral soil respiration rate, and soil temperature and moisture content of dynamic changes in different season

    3.2 高頻觀測與手動(dòng)觀測CO2通量比較

    如表2所示,與高頻觀測結(jié)果相比,CK、RR和RB全年手動(dòng)觀測土壤呼吸年通量分別低估約38%、15%和26%,而09:00—11:00時(shí)間段內(nèi)采用高頻自動(dòng)觀測系統(tǒng)觀測值計(jì)算年通量與高頻自動(dòng)觀測結(jié)果無顯著差異。此外,如圖3所示,CK和 RB手動(dòng)觀測的礦質(zhì)土壤呼吸速率均低于高頻自動(dòng)觀測,并存在顯著差異(P<0.05),RR手動(dòng)觀測的礦質(zhì)土壤呼吸速率除秋季顯著低于高頻自動(dòng)觀測,在其他季節(jié)與其均無顯著性差異(P>0.05)。CK、RR和RB的09:00—11:00 時(shí)間段內(nèi)觀測的土壤呼吸速率在不同季節(jié)中與高頻觀測的日平均呼吸速率均無顯著差異(P>0.05)(表2)。

    3.3 不同森林經(jīng)營方式下礦質(zhì)土壤呼吸與溫度和水分的關(guān)系

    通過對全年數(shù)據(jù)分析發(fā)現(xiàn),不同處理礦質(zhì)土壤呼吸與土壤溫度均呈極顯著正相關(guān)(P<0.01),其中CK(R2=0.90)和RB(R2=0.90)的土壤呼吸與溫度的相關(guān)性高于RR(R2=0.22) (表3),但在不同季節(jié)表現(xiàn)略有不同。CK、RR和RB處理在春季礦質(zhì)土壤呼吸與土壤溫度的相關(guān)性最高(R2分別為0.85,0.54和0.79)。CK處理的年均Q10為1.72,而RR處理顯著降低Q10值(1.27),但RB處理則提高了Q10值(1.84)。與土壤溫度相比,不同森林經(jīng)營方式的礦質(zhì)土壤呼吸與土壤含水量線性模型的相關(guān)性水平較差,僅在RR處理出現(xiàn)統(tǒng)計(jì)學(xué)上的相關(guān)性(表3),CK與RB處理均與土壤含水量無顯著相關(guān)性(P>0.05)。此外,采用溫度和含水量雙因素模型分別共同解釋CK、RR和RB的礦質(zhì)土壤呼吸季節(jié)變化的54.6%—86.3%、25.3%—59.4%和32.5%—75.9%,年變化的96.8%,62.8%和95.4%,擬合結(jié)果明顯優(yōu)于以溫度或含水量為單因子的經(jīng)驗(yàn)?zāi)P?且整體上顯示CK和RB的礦質(zhì)土壤呼吸與土壤溫度和含水量的相關(guān)性更強(qiáng)(CK>RB),說明本研究地土壤溫度與含水量在較大程度上共同影響著礦質(zhì)土壤呼吸。

    表 2 不同森林經(jīng)營方式不同季節(jié)礦質(zhì)土壤呼吸不同觀測方式的日通量差異

    圖3 不同森林經(jīng)營方式礦質(zhì)土壤呼吸速率在不同觀測方式的顯著性Fig.3 The significant of in different forest management of mineral soil respiration rate in different observation methods不同小寫字母表示同一處理同一季節(jié)不同觀測方式之間存在顯著差異(P<0.05),圖中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差(n=3)

    4 討論

    4.1 礦質(zhì)土壤呼吸晝夜動(dòng)態(tài)變化

    本研究結(jié)果發(fā)現(xiàn)皆伐(RR)和火燒(RB)一年后礦質(zhì)土壤呼吸低于對照(CK)(圖1),可能與煉山措施和水土流失造成土壤有機(jī)碳含量下降有關(guān)(表1)。各季節(jié)CK和RB的礦質(zhì)土壤呼吸晝夜呈單峰型變化,與Zhao等[21]研究結(jié)果相同。RB土壤呼吸速率日峰值出現(xiàn)于13:00左右,最小值出現(xiàn)于6點(diǎn)左右,與土壤溫度晝夜變化一致,國內(nèi)外研究也得出類似結(jié)果[22-23],表明RB的礦質(zhì)土壤呼吸在晝夜尺度上的變化主要受土壤溫度晝夜變化的驅(qū)動(dòng)[24]。CK的礦質(zhì)土壤呼吸在春、秋和冬季的礦質(zhì)土壤呼吸峰值出現(xiàn)時(shí)間滯后于土壤溫度(圖 2c,b,e,f,g和h),可能與CK處理林冠覆蓋較密,達(dá)到土壤微生物的最適溫度和濕度需更長的時(shí)間有關(guān)。與其不同的是,CK夏季的礦質(zhì)土壤呼吸和RR整年的礦質(zhì)土壤呼吸速率日峰值與土壤溫度的相關(guān)性并不顯著,反而隨著土壤溫度的上升其礦質(zhì)土壤呼吸速率出現(xiàn)下降趨勢 (圖 2),主要原因是在溫度不是主要限制因子時(shí),較低土壤含水量對礦質(zhì)土壤呼吸起限制作用[25-26]。本研究中除夏季的CK處理外,RR的土壤含水量均較高于CK和RR,過高的土壤含水量導(dǎo)致土壤氧氣含量下降,不利于土壤微生物活性,從而造成礦質(zhì)土壤呼吸與土壤溫度的滯后現(xiàn)象[27-28]。

    表3 不同森林經(jīng)營方式礦質(zhì)土壤呼吸與土壤溫度(T)和含水量(M)的關(guān)系

    *表示達(dá)到顯著水平(0.05),**表示達(dá)到極顯著性水平(0.01)

    4.2 手動(dòng)觀測與高頻自動(dòng)觀測

    手動(dòng)觀測和高頻自動(dòng)觀測是目前土壤呼吸的主要兩種觀測方法。本研究發(fā)現(xiàn),采用手動(dòng)觀測的結(jié)果對礦質(zhì)土壤呼吸年通量進(jìn)行估算,CK、RR和RB分別比高頻自動(dòng)觀測結(jié)果低估38%,15%和26%,且各處理之間顯差異著(P<0.05),其主要原因是手動(dòng)觀測受觀測頻率限制,導(dǎo)致手動(dòng)觀測所估算的礦質(zhì)土壤呼吸平均日通量均小于高頻自動(dòng)觀測結(jié)果[29]。CK、RB和RR的礦質(zhì)土壤呼吸速率日平均值分別為04:00—06:00,07:00—09:00和08:00—10:00之間(圖2),與森林,農(nóng)業(yè)和沼澤生態(tài)系統(tǒng)土壤呼吸速率的日平均值的結(jié)果并不一致,其認(rèn)為土壤呼吸日平均值一般出現(xiàn)于09:00—13:00之間[30,31,32]。因此,若采用09:00—13:00之間的手動(dòng)觀測數(shù)據(jù)代替日平均值,將大大低估礦質(zhì)土壤呼吸日通量。但是,高頻自動(dòng)觀測09:00—11:00的礦質(zhì)土壤呼吸速率平均值與高頻自動(dòng)觀測日均值無顯著差異(P>0.05),主要是因?yàn)楦哳l采樣可以減小兩者之間的誤差。

    4.3 礦質(zhì)土壤呼吸與土壤溫濕度關(guān)系

    本研究CK和RB的礦質(zhì)土壤呼吸與土壤溫度均呈指數(shù)相關(guān)(表3,R2=0.90),均在夏季達(dá)到最大值,與其他亞熱帶森林樣地土壤呼吸研究結(jié)果一致[33-34]。然而,RR的礦質(zhì)土壤呼吸與土壤溫度的相關(guān)性較低(表3,R2=0.22),這可能由于RR處理更新初期,植物生長迅速,土壤微生物對水分和養(yǎng)分的需求大于溫度,而我們研究結(jié)果也表明RR的礦質(zhì)土壤呼吸速率與土壤含水量相關(guān)性(表3,R2=0.47)高于溫度。與RR處理不同的是CK和RB的礦質(zhì)土壤呼吸年均與土壤含水量均無顯著線性相關(guān)(表3,P>0.05),林內(nèi)土壤含水量處于相對穩(wěn)定狀態(tài)[35],CK處理的土壤含水量日變幅均小于RR和RB(圖2)。

    相較于單一因子,土壤溫度與土壤濕度的雙因子模型能更好解釋土壤呼吸的變化[36]。本研究期內(nèi)土壤溫度和含水量共同解釋了CK和RB礦質(zhì)土壤呼吸97%和95%的變化,但RR土壤含水量和土壤溫度卻只解釋了土壤呼吸速率63%的變化,這表明皆伐地除溫度和水分外,可能還受其它因素的影響,如皆伐地仍存在未分解完的采伐剩余物和新近死亡根系,其所排放的CO2并入礦質(zhì)土壤呼吸。Karhu等[37]認(rèn)為,土壤含水量有效性可以影響土壤呼吸的溫度敏感性,土壤溫度敏感性與土壤含水量呈負(fù)相關(guān)性。本研究中皆伐后土壤含水量的增加和火燒后土壤含水量的降低(圖1),因此使礦質(zhì)土壤呼吸的Q10值呈現(xiàn)RB(1.84)>CK(1.72)>RR(1.27)的趨勢。

    5 結(jié)論

    皆伐和火燒的營林方式會(huì)顯著降低礦質(zhì)土壤呼吸速率,并以火燒影響最大。采用手動(dòng)觀測的礦質(zhì)土壤呼吸年平均日通量均低于高頻自動(dòng)觀測的結(jié)果,而采用高頻觀測09:00—11:00 時(shí)間段內(nèi)觀測數(shù)據(jù)計(jì)算日通量與高頻觀測結(jié)果無顯著差異。土壤溫度和濕度可以共同解釋次生林和火燒處理的礦質(zhì)土壤呼吸大部分變化(97%和95%),卻僅能解釋皆伐處理的礦質(zhì)土壤呼吸部分變化(63%)。

    因而,采用高頻自動(dòng)高頻觀測方式對土壤呼吸進(jìn)行觀測,不僅可以準(zhǔn)確評估土壤呼吸與溫度和水分的相關(guān)性,而且還可以對日均、月均和年均土壤呼吸通量進(jìn)行精確的評估,并且能通過高頻自動(dòng)觀測技術(shù)有可能可以量化不同生態(tài)系統(tǒng)的土壤呼吸的動(dòng)態(tài)變化,改進(jìn)其經(jīng)驗(yàn)?zāi)P汀?/p>

    [1] Tarnocai C, Canadell J G, Schuur E A G, Kuhry P, Mazhitova G, Zimov S. Soil organic carbon pools in the northern circumpolar permafrost region. Global Biogeochemical Cycles, 2009, 23(2): GB2023.

    [2] Liang N S, Inoue G, Fujinuma Y. A multichannel automated chamber system for continuous measurement of forest soil CO2efflux. Tree Physiology, 2003, 23(12): 825- 832.

    [3] Savage K E, Davidson E A. A comparison of manual and automated systems for soil CO2flux measurements: trade-offs between spatial and temporal resolution. Journal of Experimental Botany, 2003, 54(384): 891- 899.

    [4] Wang H, McConkey B, Curtin D, Cutforth H. Estimation of daily soil CO2flux using a single-time-point measurement. Canadian Journal of Soil Science, 2010, 90(3): 517- 522.

    [5] Liu Q, Edwards N T, Post W M, Gu L, Ledford J, Lenhart S. Temperature-independent diel variation in soil respiration observed from a temperate deciduous forest. Global Change Biology, 2006, 12(11): 2136- 2145.

    [6] Carbone M S, Vargas R. Automated soil respiration measurements: new information, opportunities and challenges. The New Phytologist, 2008, 177(2): 295- 297.

    [7] Decina S M, Hutyra L R, Gately C K, Getson J M, Reinmann A B, Gianotti A G S, Templer P H. Soil respiration contributes substantially to urban carbon fluxes in the greater Boston area. Environmental Pollution, 2016, 212: 433- 439.

    [8] Bond-Lamberty B, Thomson A. Temperature-associated increases in the global soil respiration record. Nature, 2010, 464(7288): 579- 582.

    [9] Savage K, Davidson E A, Tang J. Diel patterns of autotrophic and heterotrophic respiration among phenological stages. Global Change Biology, 2013, 19(4): 1151- 1159.

    [10] Sayer E J, Tanner E V J. 2010. A new approach to trenching experiments for measuring root-rhizosphere respiration in a lowland tropical forest. Soil Biology and Biochemistry, 42(2): 347- 352.

    [11] Cusack D F, Torn M S, McDowell W H, SILVER W L. The response of heterotrophic activity and carbon cycling to nitrogen additions and warming in two tropical soils. Global Change Biolog, 2010, 16(9): 2555- 2572.

    [12] Craine J M, Gelderman T M. Soil moisture controls on temperature sensitivity of soil organic carbon decomposition for a mesic grassland. Soil Biology and Biochemistry, 2011, 43(2): 455- 457.

    [13] Wang G B, Zhou Y, Xu X, Ruan H H, Wang J S. Temperature sensitivity of soil organic carbon mineralization along an elevation gradient in the Wuyi Mountains, China. PLoS One, 2013, 8(1): e53914.

    [14] 楊玉盛, 董彬, 謝錦升, 陳光水, 高人, 李靈, 王小國, 郭劍芬. 森林土壤呼吸及其對全球變化的響應(yīng). 生態(tài)學(xué)報(bào), 2004, 24(3): 583- 591.

    [15] Davidson E A, Samanta S, Caramori S S, Savage K. The dual Arrhenius and Michaelis-Menten kinetics model for decomposition of soil organic matter at hourly to seasonal time scales. Global Change Biology, 2012, 18(1): 371- 384.

    [16] Page K L, Dalal R C, Raison R J. The impact of harvesting native forests on vegetation and soil C stocks, and soil CO2, N2O and CH4fluxes. Australian Journal of Botany, 2011, 59(7): 654- 669.

    [17] Elliott K J, Vose J M. Initial effects of prescribed fire on quality of soil solution and streamwater in the southern Appalachian mountains. Southern Journal of Applied Forestry, 2005, 29(1): 5- 15.

    [18] Kulmala L, Aaltonen H, Berninger F, Kieloahoc A J, Levulad J, B?ckb J, Harib P, Kolarib P, Korhonenb J F J, Kulmalac M, Nikinmaab E, Pihlatiec M, Vesalac T, Pumpanen J. Changes in biogeochemistry and carbon fluxes in a boreal forest after the clear-cutting and partial burning of slash. Agricultural and Forest Meteorology, 2014, 188: 33- 44.

    [19] Kuzyakov Y. Sources of CO2efflux from soil and review of partitioning methods. Soil Biology and Biochemistry, 2006, 38(3): 425- 448.

    [20] Sheng H, Yang Y S, Yang Z J, Chen G S, Xie J S, Guo J F, Zou S Q. The dynamic response of soil respiration to land-use changes in subtropical China. Global Change Biology, 2010, 16(3): 1107- 1121.

    [21] Zhao J, Qi Y C, Dong Y S. Diurnal and seasonal dynamics of soil respiration in desert shrubland ofArtemisiaOrdosicaon Ordos plateau of inner Mongolia, China. Journal of Forestry Research, 2007, 18(3): 231- 235.

    [22] Wang M, Liu X T, Zhang J T, Li X J, Wang G D, Li X Y, Lu X R. Diurnal and seasonal dynamics of soil respiration at temperateLeymuschinensismeadow steppes in western Songnen plain, China. Chinese Geographical Science, 2014, 24(3): 287- 296.

    [23] Shi W Y, Zhang J G, Yan M J, Guan J H, Du S. Diurnal and seasonal dynamics of soil respiration in a Platycladus orientalis forest stand on the semiarid Loess Plateau, China. Journal of Earth Environment, 2012, 3(6): 1144- 1148.

    [24] Savage K, Phillips R, Davidson E. High temporal frequency measurements of greenhouse gas emissions from soils. Biogeosciences Discussions, 2013, 10(11): 18277- 18308.

    [25] Vargas R, Baldocchi D D, Allen M F, Bahn M, Black T A, Collins S L, Yuste J C, Hirano T, Jassal R S, Pumpanen J, Tang J W. Looking deeper into the soil: biophysical controls and seasonal lags of soil CO2production and efflux. Ecological Applications, 2010, 20(6): 1569- 1582.

    [26] Phillips C L, Nickerson N, Risk D, Bond B J. Interpreting diel hysteresis between soil respiration and temperature. Global Change Biology, 2011, 17(1): 515- 527.

    [27] Vargas R, Allen M F. Diel patterns of soil respiration in a tropical forest after Hurricane Wilma. Journal of Geophysical Research: Biogeosciences, 2008, 113(G3): G03021.

    [28] Barron-Gafford G A, Scott R L, Jenerette G D, Huxman T E. The relative controls of temperature, soil moisture, and plant functional group on soil CO2efflux at diel, seasonal, and annual scales. Journal of Geophysical Research: Biogeosciences, 2011, 116(G1): G01023.

    [29] Savage K, Davidson E A, Richardson A D. A conceptual and practical approach to data quality and analysis procedures for high-frequency soil respiration measurements. Functional Ecology, 2008, 22(6): 1000- 1007.

    [30] 蔣延玲, 周廣勝, 趙敏, 王旭, 曹銘昌. 長白山闊葉紅松林生態(tài)系統(tǒng)土壤呼吸作用研究. 植物生態(tài)學(xué)報(bào), 2005, 29(3): 411- 414.

    [31] Shi X H, Zhang X P, Yang X M, McLaughlin N B, Liang A Z, Fan R Q. An appropriate time-window for measuring soil CO2efflux: a case study on a Black soil in north-east China. Acta Agriculturae Scandinavica, Section B-Soil & Plant Science, 2012, 62(5): 449- 454.

    [32] Wang Y D, Wang H M, Wang Z L, Zhang W J, Guo C C, Wen X F, Liu Y F. Optimizing manual sampling schedule for estimating annual soil CO2efflux in a young exotic pine plantation in subtropical China. European Journal of Soil Biology, 2012, 52: 41- 47.

    [33] Tang X L, Liu S G, Zhou G Y, Zhang D Q, Zhou C Y. Soil-atmospheric exchange of CO2, CH4, and N2O in three subtropical forest ecosystems in southern China. Global Change Biology, 2006, 12(3): 546- 560.

    [34] Wang W, Peng S S, Wang T, Fang J Y. Winter soil CO2efflux and its contribution to annual soil respiration in different ecosystems of a forest-steppe ecotone, north China. Soil Biology and Biochemistry, 2010, 42(3): 451- 458.

    [35] Heinemeyer A, Di Bene C, Lloyd A R, Tortorella D, Baxter R, Huntley B, Gelsomino A, Ineson P. Soil respiration: implications of the plant-soil continuum and respiration chamber collar-insertion depth on measurement and modelling of soil CO2efflux rates in three ecosystems. European Journal of Soil Science, 2011, 62(1): 82- 94.

    [36] Olajuyigbe S, Tobin B, Saunders M, Nieuwenhuis M. Forest thinning and soil respiration in a Sitka spruce forest in Ireland. Agricultural and Forest Meteorology, 2012, 157: 86- 95.

    [37] Karhu K, Fritze H, Tuomi M, Vanhala P, Spetz P, Kitunen V, Liski J. Temperature sensitivity of organic matter decomposition in two boreal forest soil profiles. Soil Biology and Biochemistry, 2010, 42(1): 72- 82.

    High-frequency analysis of the diel patterns of mineral soil respiration under different forest management

    ZHENG Yong1,2, LIU Xiaofei1,2,*, ZHENG Wei1,2,ZHOU Jiacong1,2, SU Ruilan3, LIN Weisheng1,2, XIONG Decheng1,2, XUchao1,2, CHEN Yuehmin1,2, YANG Yusheng1,2

    1SchoolofGeographicalScience,FujianNormalUniversity,Fuzhou350007,China2StateKeyLaboratoryofSubtropicalMountainEcology(FoundedbyMinistryofScienceandTechnologyandFujianProvince),FujianNormalUniversity,Fuzhou350007,China3SanyuanDistrictForestryBureauofSanming,Sanming365000,China

    Mineral soil respiration is an important route for the loss soil carbon (C) in forest ecosystems, and it is a key factor for estimating the forest ecosystem carbon balance. Making sense of mineral soil respiration variation on temporal and spatial scales is a critical component of understanding forest ecosystem C cycle responses to climate change. Thus, high-frequency measurements have become one of the primary tools used to measure mineral soil respiration on both temporal and spatial scales. The objective of this study was to compare a manual system (Li- 8100A) with an automated system for measuring mineral soil respiration in a subtropical forest. We carried out the study in a secondaryCastanopsiscarlesiiforest with different forest management (Control, CK; Residual retention after clear cutting, RR; Residual burning after clear cutting, RB) in Sanming, Fujian. The results obtained were as follows. (1) The mineral soil respiration flux, soil temperature, and soil moisture content under different forest managements all showed obvious seasonal dynamics, and mineral soil respiration rates were in the order of CK (2.18 μmol m-2s-1) > RB (1.93 μmol m-2s-1) > RR (1.89 μmol m-2s-1). (2) The annual flux determined using high-frequency measurement was significantly higher than that determined using the manual system, whereas there was no difference between the soil respiration rate calculated between 09:00 and 11:00 (high-frequency data) for different forest managements. (3) Soil hydrothermal condition is one of the important factors affecting mineral soil respiration under different forest management. A two-factor model show that included soil temperature and moisture explained the annual variation of mineral soil respiration rate as 96.8%, 62.8% and 95.4% in CK, RR, and RB treatments, respectively. The double factors model demonstrated that combined with soil temperature and moisture was explained the annual variation of mineral soil respiration rate 96.8%, 62.8% and 95.4% in CK, RR and RB, respectively. The results obtained using this model were better than those obtained using a single-factor model. Therefore, in order to accurately assess and predict the effects of different forest management on the C flux exchange between the soil and atmosphere under future climate change, high-frequency measurements should be considered.

    different forest management; high frequency automatic observation; mineral soil respiration; diurnal patterns; carbon flux

    國家自然科學(xué)基金重點(diǎn)資助項(xiàng)目(31130013);國家自然科學(xué)基金資助項(xiàng)目(31500407);福建省教育廳B類科技資助項(xiàng)目(JB14025)

    2016- 08- 19;

    2016- 11- 04

    10.5846/stxb201608191698

    *通訊作者Corresponding author.E-mail: xfliu@fjnu.edu.cn

    鄭永, 劉小飛, 鄭蔚, 周嘉聰, 蘇瑞蘭, 林偉盛, 熊德成, 胥超, 陳岳民, 楊玉盛.應(yīng)用高頻觀測探討不同森林經(jīng)營方式下礦質(zhì)土壤呼吸的晝夜動(dòng)態(tài)特征.生態(tài)學(xué)報(bào),2017,37(1):93- 101.

    Zheng Y, Liu X F, Zheng W,Zhou J C, Su R L, Lin W S, Xiong D C, Xu C, Chen Y M, Yang Y S.High-frequency analysis of the diel patterns of mineral soil respiration under different forest management.Acta Ecologica Sinica,2017,37(1):93- 101.

    猜你喜歡
    經(jīng)營方式礦質(zhì)土壤溫度
    賽加城址先民的生產(chǎn)經(jīng)營方式
    不同產(chǎn)地三葉青中27種礦質(zhì)元素的綜合評價(jià)
    中成藥(2018年11期)2018-11-24 02:57:08
    遼東山區(qū)3種人工林土壤呼吸對土壤溫度和土壤水分的響應(yīng)
    新聞傳媒企業(yè)加強(qiáng)財(cái)務(wù)管理創(chuàng)新研究
    創(chuàng)新是生存與發(fā)展的原動(dòng)力
    中國廣播(2016年10期)2016-11-18 13:11:11
    “互聯(lián)網(wǎng)+”農(nóng)業(yè)時(shí)代新型農(nóng)業(yè)經(jīng)營方式構(gòu)建研究
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    基于礦質(zhì)元素的絞股藍(lán)產(chǎn)地鑒別研究
    稻草還田對油菜生長、土壤溫度及濕度的影響
    嫩草影院新地址| 啪啪无遮挡十八禁网站| 欧美日韩福利视频一区二区| 十八禁国产超污无遮挡网站| 欧美三级亚洲精品| 草草在线视频免费看| 亚洲av电影在线进入| 嫁个100分男人电影在线观看| 校园春色视频在线观看| 黄色日韩在线| 国产不卡一卡二| 乱人视频在线观看| 国产中年淑女户外野战色| 精品久久久久久久久久免费视频| 在线a可以看的网站| 男人舔女人下体高潮全视频| 欧美成人a在线观看| 女人被狂操c到高潮| 啦啦啦韩国在线观看视频| 亚洲成人免费电影在线观看| 精品乱码久久久久久99久播| 亚洲精品久久国产高清桃花| 国产精品久久电影中文字幕| 国产午夜精品论理片| 97碰自拍视频| 成年女人永久免费观看视频| aaaaa片日本免费| 天堂av国产一区二区熟女人妻| 久久精品91蜜桃| 欧美bdsm另类| 亚洲av熟女| 男人狂女人下面高潮的视频| 成年版毛片免费区| 国产精品一及| 给我免费播放毛片高清在线观看| 成人欧美大片| 久久欧美精品欧美久久欧美| 欧美日韩乱码在线| 一卡2卡三卡四卡精品乱码亚洲| 在线十欧美十亚洲十日本专区| 亚洲自拍偷在线| 舔av片在线| 久久久国产成人免费| 久久久久久久亚洲中文字幕 | 国产主播在线观看一区二区| 一区二区三区免费毛片| 亚洲第一电影网av| 丁香六月欧美| 91九色精品人成在线观看| 国产熟女xx| 禁无遮挡网站| 久久性视频一级片| 亚洲在线观看片| 国产又黄又爽又无遮挡在线| 中文资源天堂在线| 国产一级毛片七仙女欲春2| 亚洲欧美日韩无卡精品| 国产一区二区在线av高清观看| 欧美成人性av电影在线观看| 好男人在线观看高清免费视频| 国产精品久久久久久久久免 | 少妇熟女aⅴ在线视频| 国产一级毛片七仙女欲春2| 精品久久久久久,| 国产伦精品一区二区三区视频9| 久久久久久大精品| 毛片女人毛片| 亚洲男人的天堂狠狠| 精品一区二区三区av网在线观看| 天天一区二区日本电影三级| 国产免费av片在线观看野外av| 亚洲五月天丁香| av黄色大香蕉| 全区人妻精品视频| 成人特级av手机在线观看| bbb黄色大片| 国产在线男女| 国产精品久久久久久亚洲av鲁大| 91久久精品国产一区二区成人| 精品99又大又爽又粗少妇毛片 | 丁香欧美五月| 嫁个100分男人电影在线观看| 国产精品一区二区免费欧美| 国产三级中文精品| 欧美成人性av电影在线观看| 99热这里只有是精品50| 欧美又色又爽又黄视频| 我要搜黄色片| 天堂影院成人在线观看| 免费在线观看影片大全网站| 美女xxoo啪啪120秒动态图 | 精品午夜福利在线看| 久久久色成人| 搞女人的毛片| 男女下面进入的视频免费午夜| 天天躁日日操中文字幕| 精品久久久久久久久久免费视频| 又粗又爽又猛毛片免费看| 黄色视频,在线免费观看| 成人鲁丝片一二三区免费| 人妻制服诱惑在线中文字幕| 精品一区二区免费观看| 国产真实乱freesex| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久亚洲 | 免费黄网站久久成人精品 | 成人性生交大片免费视频hd| 成人永久免费在线观看视频| eeuss影院久久| 日本 欧美在线| 国产综合懂色| 中文字幕av成人在线电影| 黄色视频,在线免费观看| 亚洲欧美精品综合久久99| 久久久久国产精品人妻aⅴ院| 国产免费av片在线观看野外av| 久久性视频一级片| 亚洲一区二区三区不卡视频| 色综合亚洲欧美另类图片| 午夜精品在线福利| 亚洲午夜理论影院| 亚洲成人中文字幕在线播放| 1024手机看黄色片| 两性午夜刺激爽爽歪歪视频在线观看| 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 亚洲三级黄色毛片| 国产野战对白在线观看| 18美女黄网站色大片免费观看| 精品一区二区三区av网在线观看| 成人三级黄色视频| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片 | 欧美在线一区亚洲| 国内精品久久久久精免费| 中文字幕av成人在线电影| 国产精品女同一区二区软件 | 亚洲色图av天堂| 丰满人妻一区二区三区视频av| 国产欧美日韩一区二区精品| 亚洲,欧美,日韩| 最新中文字幕久久久久| 一本一本综合久久| 少妇熟女aⅴ在线视频| 国产私拍福利视频在线观看| 国产成人av教育| 欧美乱色亚洲激情| 免费电影在线观看免费观看| 成年免费大片在线观看| 国产亚洲精品久久久久久毛片| 精品久久久久久久久av| 99久国产av精品| 亚洲人成网站在线播放欧美日韩| 18禁在线播放成人免费| 亚洲人成网站高清观看| eeuss影院久久| 男插女下体视频免费在线播放| 亚洲av中文字字幕乱码综合| 精品一区二区免费观看| 国产日本99.免费观看| 亚洲av五月六月丁香网| 成人特级av手机在线观看| 99久国产av精品| 亚洲成人免费电影在线观看| 国内毛片毛片毛片毛片毛片| 观看免费一级毛片| 美女免费视频网站| 麻豆成人午夜福利视频| h日本视频在线播放| 国产精品一区二区三区四区久久| 婷婷丁香在线五月| 国产一级毛片七仙女欲春2| 极品教师在线免费播放| 又黄又爽又刺激的免费视频.| 亚洲av成人不卡在线观看播放网| www.www免费av| 丁香六月欧美| 国产91精品成人一区二区三区| 精品熟女少妇八av免费久了| 小说图片视频综合网站| 深爱激情五月婷婷| 夜夜夜夜夜久久久久| 高潮久久久久久久久久久不卡| 中文字幕av在线有码专区| 在线观看66精品国产| 在线天堂最新版资源| 国产精品,欧美在线| 午夜福利在线观看吧| 国产高清视频在线观看网站| 欧美成人a在线观看| 人妻久久中文字幕网| 亚洲,欧美,日韩| 欧美日本亚洲视频在线播放| 免费看美女性在线毛片视频| 国产成人aa在线观看| 免费人成在线观看视频色| 精品一区二区三区视频在线观看免费| 久久婷婷人人爽人人干人人爱| 一区二区三区四区激情视频 | 日日摸夜夜添夜夜添小说| 深夜a级毛片| 午夜福利视频1000在线观看| 99久久九九国产精品国产免费| 男女床上黄色一级片免费看| 中文字幕免费在线视频6| 亚洲成人久久性| 一个人免费在线观看电影| 九色国产91popny在线| 亚洲精品在线观看二区| 日韩高清综合在线| 免费av观看视频| 久久欧美精品欧美久久欧美| 1000部很黄的大片| 欧美高清成人免费视频www| 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| 日韩欧美精品免费久久 | 黄片小视频在线播放| 男人和女人高潮做爰伦理| 日本一本二区三区精品| 免费看a级黄色片| 嫩草影院入口| 成人一区二区视频在线观看| 日韩欧美精品v在线| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| 久久久久久久久中文| 免费一级毛片在线播放高清视频| 色av中文字幕| 可以在线观看毛片的网站| 老熟妇乱子伦视频在线观看| 老熟妇仑乱视频hdxx| 日韩成人在线观看一区二区三区| 中国美女看黄片| 麻豆一二三区av精品| 少妇丰满av| 毛片女人毛片| 伦理电影大哥的女人| 一二三四社区在线视频社区8| 国产精品亚洲av一区麻豆| 中文亚洲av片在线观看爽| 午夜福利免费观看在线| 午夜福利视频1000在线观看| 一级黄片播放器| 国产精品美女特级片免费视频播放器| 蜜桃亚洲精品一区二区三区| 午夜激情福利司机影院| 51午夜福利影视在线观看| 99国产极品粉嫩在线观看| 国产精品久久久久久亚洲av鲁大| av国产免费在线观看| 久久精品国产清高在天天线| 51午夜福利影视在线观看| 在线看三级毛片| 真人做人爱边吃奶动态| 又粗又爽又猛毛片免费看| 床上黄色一级片| 国产真实乱freesex| 亚洲久久久久久中文字幕| 成人美女网站在线观看视频| 精品一区二区三区视频在线观看免费| 亚洲自偷自拍三级| 成人av在线播放网站| 精品人妻一区二区三区麻豆 | 午夜日韩欧美国产| av天堂在线播放| 色综合站精品国产| 一区二区三区免费毛片| 欧美激情在线99| 亚洲精品456在线播放app | 麻豆久久精品国产亚洲av| 国产欧美日韩一区二区三| 天美传媒精品一区二区| 真人一进一出gif抽搐免费| 亚洲三级黄色毛片| 91久久精品电影网| 亚洲激情在线av| 性插视频无遮挡在线免费观看| 亚洲av不卡在线观看| 欧美日韩瑟瑟在线播放| 少妇人妻一区二区三区视频| 日本a在线网址| 99热这里只有精品一区| 波野结衣二区三区在线| 热99在线观看视频| 久久6这里有精品| 日本黄色片子视频| 在线播放无遮挡| 成人国产一区最新在线观看| 亚洲无线在线观看| 99在线视频只有这里精品首页| 亚洲av不卡在线观看| 国内精品一区二区在线观看| 成人特级黄色片久久久久久久| 欧美最新免费一区二区三区 | 亚州av有码| 国产高清有码在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 日本免费a在线| 国产成人av教育| 午夜日韩欧美国产| 午夜福利18| 噜噜噜噜噜久久久久久91| 成人av在线播放网站| 99热这里只有是精品50| 日韩免费av在线播放| 亚洲精品在线观看二区| 男女下面进入的视频免费午夜| 亚洲av免费在线观看| 欧美性猛交黑人性爽| 欧美在线一区亚洲| 在现免费观看毛片| 日韩 亚洲 欧美在线| 国产日本99.免费观看| www.熟女人妻精品国产| 搞女人的毛片| 成年女人看的毛片在线观看| 又黄又爽又免费观看的视频| 一个人看的www免费观看视频| 三级男女做爰猛烈吃奶摸视频| 日韩中字成人| 嫩草影院新地址| 欧美中文日本在线观看视频| 久久精品国产亚洲av涩爱 | 精品久久久久久久久亚洲 | 99久久99久久久精品蜜桃| 波多野结衣高清作品| 精品久久久久久久久久久久久| 日韩成人在线观看一区二区三区| 国产麻豆成人av免费视频| 蜜桃久久精品国产亚洲av| 国产精品亚洲av一区麻豆| 日本 av在线| 九九热线精品视视频播放| 18禁在线播放成人免费| 深爱激情五月婷婷| 97热精品久久久久久| 色综合欧美亚洲国产小说| 国产亚洲精品综合一区在线观看| 一进一出抽搐gif免费好疼| 亚洲激情在线av| 99热这里只有精品一区| av欧美777| 天天躁日日操中文字幕| 亚洲av电影不卡..在线观看| 久久精品国产亚洲av涩爱 | 亚洲人成网站在线播| 又黄又爽又免费观看的视频| 深夜a级毛片| 国产精品一及| 亚洲五月天丁香| 两性午夜刺激爽爽歪歪视频在线观看| 日本免费a在线| 91九色精品人成在线观看| 欧美性猛交黑人性爽| 免费人成在线观看视频色| www.999成人在线观看| 中文字幕人妻熟人妻熟丝袜美| 少妇丰满av| 特大巨黑吊av在线直播| 婷婷色综合大香蕉| 特大巨黑吊av在线直播| av欧美777| 精品99又大又爽又粗少妇毛片 | 欧美日韩乱码在线| 国产成+人综合+亚洲专区| 亚洲无线在线观看| 国产亚洲精品av在线| 久久久久久大精品| 亚洲无线观看免费| 18禁在线播放成人免费| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av熟女| 国产精品不卡视频一区二区 | 欧美乱色亚洲激情| 亚洲av一区综合| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 国产亚洲精品综合一区在线观看| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 国产色爽女视频免费观看| 精品乱码久久久久久99久播| 日本一本二区三区精品| 国产欧美日韩精品亚洲av| 18+在线观看网站| www.999成人在线观看| 看免费av毛片| 琪琪午夜伦伦电影理论片6080| 99久久精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 人人妻人人看人人澡| 午夜福利免费观看在线| 久久精品国产亚洲av涩爱 | 亚洲国产高清在线一区二区三| 日本一二三区视频观看| 三级国产精品欧美在线观看| 精品人妻熟女av久视频| 日本与韩国留学比较| 麻豆av噜噜一区二区三区| 亚洲美女黄片视频| 亚洲欧美精品综合久久99| 日韩中字成人| 直男gayav资源| 亚洲成人久久性| 狂野欧美白嫩少妇大欣赏| 成人国产一区最新在线观看| 国产单亲对白刺激| 国产伦精品一区二区三区视频9| 午夜福利成人在线免费观看| 一区福利在线观看| 中文字幕人妻熟人妻熟丝袜美| 88av欧美| 看片在线看免费视频| 九九热线精品视视频播放| 亚洲无线观看免费| 国产精品久久久久久人妻精品电影| 在线十欧美十亚洲十日本专区| 久久6这里有精品| 变态另类成人亚洲欧美熟女| 最近在线观看免费完整版| 日本免费一区二区三区高清不卡| 狂野欧美白嫩少妇大欣赏| 精品人妻1区二区| 99久国产av精品| 久久婷婷人人爽人人干人人爱| 日韩欧美国产在线观看| 一二三四社区在线视频社区8| 真人做人爱边吃奶动态| 天美传媒精品一区二区| 欧美最新免费一区二区三区 | 看十八女毛片水多多多| 毛片女人毛片| 欧美又色又爽又黄视频| 如何舔出高潮| 丰满的人妻完整版| 亚洲在线观看片| 色哟哟·www| 小蜜桃在线观看免费完整版高清| 亚洲av第一区精品v没综合| 伦理电影大哥的女人| 桃红色精品国产亚洲av| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 最新中文字幕久久久久| 国产色爽女视频免费观看| 宅男免费午夜| 免费观看精品视频网站| 偷拍熟女少妇极品色| 少妇高潮的动态图| 亚洲欧美清纯卡通| 日本免费一区二区三区高清不卡| 老司机午夜福利在线观看视频| 99久久精品热视频| 99精品在免费线老司机午夜| 欧美日韩乱码在线| 国产精品亚洲一级av第二区| bbb黄色大片| 丁香欧美五月| 精品久久久久久久久久久久久| 天堂√8在线中文| 成人精品一区二区免费| 欧美zozozo另类| 小说图片视频综合网站| 日韩 亚洲 欧美在线| 亚洲国产精品sss在线观看| 99热这里只有精品一区| 长腿黑丝高跟| 欧美+日韩+精品| 最近最新免费中文字幕在线| 三级国产精品欧美在线观看| 日本撒尿小便嘘嘘汇集6| 97人妻精品一区二区三区麻豆| 久久精品国产亚洲av天美| 久久久久久久久久黄片| 9191精品国产免费久久| 中亚洲国语对白在线视频| 国产免费一级a男人的天堂| 桃色一区二区三区在线观看| 国产精品久久久久久精品电影| 精品国内亚洲2022精品成人| 国产伦人伦偷精品视频| 在线观看免费视频日本深夜| 成人性生交大片免费视频hd| 五月伊人婷婷丁香| 九色国产91popny在线| 亚洲国产欧美人成| 国产精品国产高清国产av| 免费观看的影片在线观看| 12—13女人毛片做爰片一| 欧美黄色淫秽网站| 他把我摸到了高潮在线观看| 久久久久久久久中文| 欧美在线一区亚洲| 91久久精品电影网| 精品日产1卡2卡| 久久婷婷人人爽人人干人人爱| 久久精品国产亚洲av天美| 色噜噜av男人的天堂激情| 久久久国产成人精品二区| 久久香蕉精品热| 亚洲欧美日韩卡通动漫| 人人妻人人看人人澡| 国产亚洲欧美在线一区二区| 人妻丰满熟妇av一区二区三区| 欧美在线黄色| 国产成人欧美在线观看| 国产野战对白在线观看| 欧美中文日本在线观看视频| 91av网一区二区| 久久久国产成人精品二区| 一级a爱片免费观看的视频| 美女 人体艺术 gogo| 又粗又爽又猛毛片免费看| 久久久久九九精品影院| 成年女人看的毛片在线观看| 男人和女人高潮做爰伦理| 人妻夜夜爽99麻豆av| 免费看日本二区| 亚洲综合色惰| 国产亚洲精品av在线| 日本a在线网址| 亚洲一区二区三区不卡视频| 九九热线精品视视频播放| 小蜜桃在线观看免费完整版高清| 国产一区二区在线av高清观看| 怎么达到女性高潮| 久久中文看片网| 十八禁国产超污无遮挡网站| 国产亚洲精品综合一区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 又爽又黄无遮挡网站| 最好的美女福利视频网| 国产伦人伦偷精品视频| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 波多野结衣巨乳人妻| 村上凉子中文字幕在线| 久久午夜亚洲精品久久| 亚洲avbb在线观看| 欧美区成人在线视频| 亚洲五月婷婷丁香| 成人特级黄色片久久久久久久| 成人无遮挡网站| 麻豆成人av在线观看| 51午夜福利影视在线观看| 赤兔流量卡办理| 亚洲国产欧美人成| 色哟哟·www| 中文字幕久久专区| 久久精品影院6| 久久久精品大字幕| 精品久久国产蜜桃| 国产伦在线观看视频一区| 亚洲黑人精品在线| 动漫黄色视频在线观看| 国产中年淑女户外野战色| 91av网一区二区| 在线播放国产精品三级| 丰满人妻一区二区三区视频av| 日韩欧美国产在线观看| 国产精品av视频在线免费观看| 99热这里只有精品一区| 久久99热这里只有精品18| 日韩欧美在线二视频| 日本精品一区二区三区蜜桃| 国产色婷婷99| 亚洲狠狠婷婷综合久久图片| 久久久久国产精品人妻aⅴ院| 97热精品久久久久久| 99国产精品一区二区三区| 天天躁日日操中文字幕| 亚洲精品在线美女| 永久网站在线| 精品熟女少妇八av免费久了| 国产视频一区二区在线看| 精品久久久久久久末码| 一区二区三区激情视频| 国产精品美女特级片免费视频播放器| 亚洲 国产 在线| 51国产日韩欧美| 国产精品久久久久久亚洲av鲁大| 老司机午夜十八禁免费视频| 嫩草影视91久久| 99热6这里只有精品| 欧美bdsm另类| 老女人水多毛片| 丰满乱子伦码专区| 日韩精品中文字幕看吧| 久久久久九九精品影院| 日韩人妻高清精品专区| 亚洲五月婷婷丁香| 成年女人看的毛片在线观看| 麻豆国产av国片精品| 老熟妇仑乱视频hdxx| 国产美女午夜福利| 日韩人妻高清精品专区| 亚洲欧美精品综合久久99| 欧美黄色片欧美黄色片| 国产精品野战在线观看| 亚洲av电影在线进入| 国产白丝娇喘喷水9色精品| 一区二区三区免费毛片| 国产精品一区二区性色av| 亚洲一区高清亚洲精品| 亚洲精品在线观看二区| 午夜福利成人在线免费观看| 日韩欧美三级三区| 国产熟女xx| 性色avwww在线观看|