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

    馬尾松林土壤呼吸組分對不同營林措施的響應(yīng)

    2016-11-28 07:23:32肖文發(fā)曾立雄黃志霖譚本旺
    生態(tài)學(xué)報(bào) 2016年17期
    關(guān)鍵詞:森林土壤礦質(zhì)土壤濕度

    雷 蕾,肖文發(fā),曾立雄,*,黃志霖,譚本旺

    1 中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所, 國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100091 2 南京林業(yè)大學(xué)南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,南京 210037 3 湖北省秭歸縣林業(yè)局,宜昌 443631

    ?

    馬尾松林土壤呼吸組分對不同營林措施的響應(yīng)

    雷 蕾1,2,肖文發(fā)1,2,曾立雄1,2,*,黃志霖1,2,譚本旺3

    1 中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所, 國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100091 2 南京林業(yè)大學(xué)南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,南京 210037 3 湖北省秭歸縣林業(yè)局,宜昌 443631

    針對不同營林措施(對照、除灌、采伐1(15%)、采伐2(70%)后的三峽庫區(qū)馬尾松飛播林,采用LI- 8100對其土壤呼吸組分的呼吸速率和土壤溫度、濕度進(jìn)行為期1年的連續(xù)觀測分析表明,不同營林措施對土壤呼吸組分的影響不同。1)觀測期內(nèi),各營林措施下凋落物層呼吸速率差異并不顯著,對照、除灌、采伐1、采伐2的根呼吸速率均值分別為:1.00、 0.83、0.86、1.11μmolCO2m-2s-1;采伐處理下礦質(zhì)土壤呼吸顯著高于對照和除灌(P<0.05); 2)與對照相比,營林措施并未顯著改變凋落物呼吸對于土壤總呼吸的貢獻(xiàn)率(18.78%—23.70%),但降低了根呼吸的貢獻(xiàn)率,其中以采伐1最為顯著(P<0.05);除灌的礦質(zhì)土壤呼吸貢獻(xiàn)率(37.00%)與對照(38.32%)相近,而采伐1(45.63%)和采伐2(43.07%)均顯著增加了礦質(zhì)土壤呼吸的貢獻(xiàn)率,礦質(zhì)土壤呼吸的變化是造成采伐措施下土壤呼吸變化的主要土壤呼吸組分;3)營林后僅采伐2措施下土壤溫濕度顯著高于對照,土壤溫濕度雙因子模型較單因子模型能更好的解釋土壤呼吸組分變化,但僅能解釋其部分變化(4.6%—59.3%),仍需對營林后其他相關(guān)因子進(jìn)行深入的綜合研究。

    營林措施;土壤呼吸組分;凋落物層呼吸;根系呼吸;礦質(zhì)土壤呼吸

    土壤碳庫的微小的變化都會對全球碳循環(huán)和氣候產(chǎn)生重要的影響[1- 3]。土壤呼吸作為森林土壤碳庫主要輸出過程之一,占森林生態(tài)系統(tǒng)碳排放的30%—80%[4],是森林參與全球碳循環(huán)的關(guān)鍵部分。多樣的森林類型和不同營林措施的干擾造成了森林土壤碳庫維持機(jī)制、碳固定過程的研究較大的不確定性[5]。森林管理措施如采伐、火燒、林床清理等,可能通過改變地上植被組成、土壤環(huán)境因子、有機(jī)物質(zhì)分解速率等進(jìn)而干擾著土壤呼吸過程[6-7]。近些年來許多學(xué)者針對采伐措施對土壤總呼吸速率影響進(jìn)行了大量研究[7- 9],土壤呼吸速率呈現(xiàn)增加[10-11]、減少[12]或不變[7, 13]的結(jié)果,仍需對特定區(qū)域不同森林經(jīng)營方式下森林土壤呼吸進(jìn)行研究,以定量預(yù)測營林措施對于碳循環(huán)的影響進(jìn)而制定出合理的森林經(jīng)營方案。土壤呼吸主要由根呼吸、礦質(zhì)土壤呼吸以及凋落物層呼吸三部分組成,森林土壤呼吸組分與環(huán)境要素關(guān)系密切[14-15],估算營林措施后森林土壤呼吸組分的變化是評估營林干擾對森林土壤碳庫過程的首要內(nèi)容[9]。

    馬尾松(Pinusmassoniana)是中國南方主要造林樹種,是三峽庫區(qū)分布面積最大的森林類型[16]。目前關(guān)于本區(qū)域馬尾松林人為管理后土壤呼吸變化尚未見報(bào)道。本研究以三峽庫區(qū)人為管理后馬尾松飛播林為研究對象,對不同營林措施(對照(未采伐)、除灌(清除灌叢)、采伐1(強(qiáng)度15%)、采伐2(強(qiáng)度70%))后土壤呼吸組分以及5cm深處土壤溫濕度進(jìn)行連續(xù)觀測,以探討不同營林措施后馬尾松林土壤呼吸組分的變化情況以及土壤溫濕度對其變化的影響作用,旨在為本區(qū)域科學(xué)經(jīng)營馬尾松林提供理論支撐以促進(jìn)馬尾松固碳效應(yīng)。

    1 研究區(qū)概況

    研究地點(diǎn)位于湖北省秭歸縣三峽庫區(qū)九嶺頭林場,地理位置30°59′N, 110°47′E,海拔156—2030.6m,年均氣溫16.9℃,年降水量1000—1250mm,多集中于4—9月,屬于亞熱帶大陸性季風(fēng)氣候。土壤類型以黃壤、黃棕壤為主[17]。調(diào)查樣地內(nèi)馬尾松林為20世紀(jì)70年代飛播造林,樣地內(nèi)林木分布均勻,馬尾松為主要優(yōu)勢種,蓋度為80%,伴生有少量光皮樺(Betulaluminifera)、漆樹(Toxicodendronvernicifluum)、杉木(Cunninghamialanceolata)。灌木主要有火棘(Pyracanthafortuneana)、胡枝子(Lespedezabicolor)、木姜子(Litseapungens)等,草本主要有狗脊(Woodwardiajaponica)、苔草(Carextristachya)、三脈紫菀(Asterageratoides)、中日金星蕨(Parathelypterisnipponica)等。

    樣地設(shè)置于三峽庫區(qū)九嶺頭林場經(jīng)營后的馬尾松飛播林內(nèi),樣地內(nèi)采用典型采樣法,于2013年9月按照不同營林措施,除灌、采伐1、采伐2,分別設(shè)立3塊20m×20m固定樣地,除灌(SM)措施清除樣地內(nèi)所有灌木,并對灌木清理中產(chǎn)生的剩余物(灌叢枝葉等)進(jìn)行仔細(xì)清除;采伐1(H1)措施清除樣地內(nèi)徑階在4cm以上的非馬尾松(主要為光皮樺、漆樹、杉木,包括部分高大灌木,如木姜子、火棘等),強(qiáng)度為15%(按胸高斷面積計(jì)算);采伐2(H2)措施清除樣地內(nèi)平均胸徑(17.9cm)以上的馬尾松,強(qiáng)度為70%。同時(shí)在固定樣地相鄰處,土壤條件基本相同的林地,設(shè)置對照固定樣地3塊(20m×20m)。伐木作業(yè)采用人力油鋸伐木,僅對采伐的樹干進(jìn)行了移除,不對采伐產(chǎn)生的剩余物如枝葉等進(jìn)行清除。各種經(jīng)營措施并未對活地被層草本以及地表枯落物進(jìn)行任何處理,上述營林措施于2013年10月中旬完成(表1)。

    括號內(nèi)數(shù)值為標(biāo)準(zhǔn)誤

    2 研究方法

    2.1 土壤呼吸及土壤溫濕度測定

    在設(shè)定的固定樣地內(nèi)選擇條件基本一致的地點(diǎn)隨機(jī)設(shè)置3種觀測小區(qū)(1m×1m):對照小區(qū)(保留凋落物層+保留根系);去凋小區(qū)(去除凋落物層+保留根系);斷根小區(qū)(保留凋落物層+切斷根系,在小區(qū)四周垂直挖深至1m左右,直至看不到根系,切斷根系,并在四周埋入厚塑料板以阻止外圍根系向樣方內(nèi)生長)。每種小區(qū)重復(fù)3次。每個小區(qū)內(nèi)各布設(shè)1個直徑 20 cm、深度 8 cm 的聚氯乙烯(polyvinylchloride,PVC,5cm埋入地下[18])土壤環(huán)置放于土壤表層內(nèi)。觀測小區(qū)內(nèi)部均無灌木和喬木,草本、枯落物等分布均勻。土壤呼吸由 Li- 8100 便攜式土壤呼吸儀測得。同時(shí)用Li- 8100自帶的土壤溫濕度探頭測定5cm土壤的溫度以及體積含水量(本文中土壤濕度均為體積含水量),每個土壤呼吸圈附近測定 3 次取平均值。測定時(shí)間為2013年11月—2014年10月,避開雨天測定,生長季(3—9月)每月于中旬和月底各測定1次,共測定2次,非生長季每月于月底測定1次;其中,2014年1月、2月由于下雪封路未能進(jìn)行測定,4月生長季由于連續(xù)降雨,僅測定1次。

    2.2 數(shù)據(jù)統(tǒng)計(jì)分析

    月土壤呼吸組分的呼吸速率是每月測量值的平均值,土壤呼吸組分的呼吸速率均值是觀測期內(nèi)(2013年11月—2014年10月)測量的平均值。

    土壤呼吸組分計(jì)算公式為[19]:

    RL=R1-R2;RR=R1-R3;RM=R3-RL

    (1)

    式中,RL為凋落物層呼吸速率(μmolCO2m-2s-1);RR為根系呼吸速率(μmolCO2m-2s-1);RM為礦質(zhì)土壤呼吸速率(μmolCO2m-2s-1);R1為對照小區(qū)土壤呼吸速率(μmolCO2m-2s-1)即總呼吸Rtot;R2為去凋小區(qū)土壤呼吸速率(μmolCO2m-2s-1);R3為斷根小區(qū)土壤呼吸速率(μmolCO2m-2s-1)。斷根小區(qū)于首次測定前3個月設(shè)置,相關(guān)研究認(rèn)為挖根后3個月可忽略死根分解的干擾[15]。

    各土壤呼吸組分對土壤呼吸速率貢獻(xiàn)率(%)通過式(1)中計(jì)算的各組分值除以R1乘以100計(jì)算獲得。凋落物層、根系的土壤溫度T1(℃)和濕度W1(%)均為對照小區(qū)土壤的溫度和濕度;礦質(zhì)土壤溫度T2(℃)和濕度W2(%)均為斷根小區(qū)土壤的溫度和濕度。

    經(jīng)驗(yàn)指數(shù)方程用于描述土壤呼吸與土壤表層溫度的關(guān)系(方程(2)),

    Rs=aebT

    (2)

    式中,Rs為測定的土壤呼吸速率(μmolCO2m-2s-1),T為表層5cm土壤溫度(℃),a是溫度為0℃的土壤呼吸,b為溫度反應(yīng)系數(shù)。

    土壤呼吸速率與土壤濕度(W)之間的關(guān)系通過二次回歸方程(3)計(jì)算:

    Rs=b1+b2×W+b3×W2

    (3)

    式中,W為土壤體積含水量(%),b1、b2、b3為方程擬合參數(shù)。

    為檢驗(yàn)土壤濕度與土壤溫度對于土壤呼吸速率的疊加效應(yīng),將土壤溫度、土壤含水量及其乘積作為獨(dú)立變量引入到對數(shù)變換后的方程(2)中,如下式所示:

    lnRs=c1+c2×T+c3×W+c4×T×W

    (4)

    式中,T、W如上所述,c1、c2、c3、c4為方程擬合參數(shù)。

    所有的統(tǒng)計(jì)分析都在SPSS19軟件中進(jìn)行,用Repeated measured ANOVA檢驗(yàn)土壤呼吸組分季節(jié)變化的顯著性,用SigmaPlot 11.0軟件作圖。

    3 結(jié)果與分析

    3.1 呼吸組分時(shí)間變異性

    從觀測期內(nèi)土壤呼吸速率月均值來看,與對照樣地相比,除灌措施顯著降低了森林土壤總呼吸Rtot,降低了17.74%,采伐2措施則對森林土壤總呼吸Rtot具有顯著的促進(jìn)作用,提高了17.34%,采伐1措施對森林土壤總呼吸Rtot沒有產(chǎn)生明顯的影響。從不同土壤呼吸組分來看,除灌措施顯著降低了礦質(zhì)土壤呼吸RM,采伐1和采伐2措施均顯著提高了礦質(zhì)土壤呼吸RM;除灌和采伐1措施都顯著降低了森林土壤的根呼吸RR,采伐2措施對森林土壤的根呼吸RR無顯著影響;而所有營林措施對凋落物呼吸RL影響均不顯著,均值在0.47—0.53之間(表2)。

    表2 不同營林措施土壤呼吸組分平均值多重比較

    同列a、b、c表示顯著性差異組別(α=0.005)

    由圖1不同呼吸組分土壤呼吸速率月變化可知,不同營林措施下土壤總呼吸Rtot、RR、RM月均值在夏季生長旺盛期7—8月期間達(dá)到最高,而RL中除灌措施(SM)下土壤呼吸速率均值在6月達(dá)到最高,其余營林措施下RL也均在7—8月間達(dá)到最高(圖1)。相較于其他呼吸組分,凋落物呼吸RL月際變化較小,對照(CK)、除灌(SM)、采伐1(H1)、采伐2(H2)月均最低值范圍為0.20—0.31μmolCO2m-2s-1,最高值范圍為0.94—1.08μmolCO2m-2s-1。根呼吸RR和礦質(zhì)土壤呼吸RM月際變化較大,其中H2的RR與RM在觀測期內(nèi)變化波動最大,RR、RM月均值差值分別為:1.65、1.75μmolCO2m-2s-1,CK次之,其RR、RM月均值差值分別為1.17、1.16μmolCO2m-2s-1。

    3.2 不同營林措施下土壤呼吸組分對土壤呼吸貢獻(xiàn)率

    不同營林措施下根呼吸RR和礦質(zhì)土壤呼吸RM是土壤總呼吸的主要組成部分,兩者對于總呼吸的貢獻(xiàn)率均分別在35.75%—40.44%和37.00%—45.63%之間,而凋落物呼吸貢獻(xiàn)率則在18.62%—23.70%之間。

    不同營林措施下SM的RL貢獻(xiàn)率最大(23.70%),RM貢獻(xiàn)率最小(37.00%);H1的RL和RR貢獻(xiàn)率在這幾種處理中均最小分別為18.78%、35.81%,而其RM貢獻(xiàn)率最大為45.42%;對照CK的RR貢獻(xiàn)率最大,40.44%。不同營林措施下土壤呼吸組分的貢獻(xiàn)率差異性也呈現(xiàn)出不同的結(jié)果,采伐措施(H1和H2)RL的貢獻(xiàn)率顯著低于除灌措施(SM)(P<0.01),而與對照(CK)差異并不顯著;對照樣地的RR貢獻(xiàn)率顯著高于采伐1(P<0.01)而與其它措施較相近;礦質(zhì)土壤呼吸RM貢獻(xiàn)率隨著干擾的強(qiáng)度而顯著分化為兩組,對照和除灌顯著低于采伐措施(P<0.01)(表3)。

    表3 不同營林措施下土壤呼吸組分對土壤呼吸貢獻(xiàn)率

    3.3 土壤溫度與土壤呼吸組分

    圖2 不同營林措施下土壤溫度和土壤濕度變化 Fig.2 Soil temperature and soil moisture under different forest management誤差項(xiàng)為標(biāo)準(zhǔn)誤Error bars show SE

    不同營林措施下土壤溫度均在7月達(dá)到最高值(圖2),不同措施下僅H2的土壤溫度月均值顯著高于CK(表4)。土壤呼吸組分均與土壤溫度呈現(xiàn)顯著地正相關(guān)關(guān)系(P<0.05),但通過本研究選取的土壤呼吸與溫度關(guān)系方程擬合可知,溫度僅能很好的解釋部分呼吸組分的變化,如除灌和采伐2措施下土壤溫度僅能解釋土壤呼吸變化的4.2%和4.5%,對照和采伐1略高,分別為21.4%和16.5%。而根呼吸RR與土壤溫度擬合程度高于RL,R2從高到低依次為H2(0.498)>CK(0.352)>SM(0.339)>H1(0.203)。不同營林措施下土壤溫度與礦質(zhì)土壤呼吸RM關(guān)系呈現(xiàn)顯著地差異性,采伐措施下土壤溫度與RM擬合模型的R2顯著高于對照和除灌,土壤溫度能更好的解釋采伐措施下RM的變化規(guī)律(表5)。

    表4 不同營林措施下土壤溫度與濕度平均值多重比較

    3.4 土壤濕度與土壤呼吸組分

    土壤濕度與土壤呼吸組分的關(guān)系較于土壤溫度更加復(fù)雜,土壤濕度對于土壤各呼吸組分并未表現(xiàn)出一致的相關(guān)性。除對照RM(R2=0.001,P=0.966)外,土壤濕度能更好的解釋其他處理措施下RM的變化,其中以采伐1為最優(yōu)(R2=0.235,P<0.01)。而凋落物呼吸RL則相反,除灌、采伐1、采伐2下土壤濕度的變化并不能反映RL的變化,而對照樣地的RL與土壤濕度存在顯著的相關(guān)關(guān)系,但模型模擬效果并不理想(R2=0.09,P<0.01)。不同營林措施下,除灌和采伐2的土壤濕度分別能解釋RR變化的19.8%和14.8%(P<0.01),對照和采伐1的土壤濕度則并未與RR呈現(xiàn)出一定的相關(guān)性(表6)??傊?,土壤濕度僅能從較小的程度影響不同呼吸組分的變化。

    表5 不同營林措施下土壤呼吸組分與土壤溫度(T)關(guān)系模型參數(shù)

    表6 不同營林措施下土壤呼吸組分與土壤濕度(W)關(guān)系模型參數(shù)

    3.5 土壤溫度和土壤濕度對土壤呼吸組分協(xié)同作用

    對比表5、表6、表7可知,土壤溫度與土壤濕度雙因子的協(xié)同作用能更好的闡明土壤呼吸各組分的變異規(guī)律。與對照相比,其他營林措施的土壤溫度和土壤濕度對凋落物呼吸RL的協(xié)同作用減弱。不同于RL,不同營林措施下土壤溫度、土壤濕度與根呼吸RR和礦質(zhì)土壤呼吸RM的關(guān)系并未表現(xiàn)出一致性,采伐2措施提高了土壤溫濕度對于RR和RM的協(xié)同作用,溫濕度能共同解釋RR和RM的50.2%和59.3%,分別高出對照13.4%和35%,而除灌措施則相反。較之對照,采伐1降低了土壤溫濕度對于RR的影響,但提高了土壤溫濕度對于RM的復(fù)合作用。

    4 討論

    4.1 不同措施下土壤呼吸總量及其組分變化

    營林措施通常會造成地上碳儲量和碳通量的顯著改變,但是其對于土壤的碳通量影響目前并未達(dá)到統(tǒng)一的結(jié)論[20]。本研究中,營林措施實(shí)施后至第1個生長季結(jié)束時(shí),3種營林措施對森林土壤總呼吸的影響差異明顯(圖1,表2)。除灌措施的森林土壤總呼吸Rtot降低了17.74%,采伐2措施的森林土壤總呼吸Rtot提高了17.34%,采伐1措施對森林土壤總呼吸Rtot沒有產(chǎn)生明顯的影響。進(jìn)一步分析表明,不同營林措施下根呼吸RR和礦質(zhì)土壤呼吸RM是土壤總呼吸的主要組成部分,兩者對于總呼吸的貢獻(xiàn)率均分別在35.75%—40.44%和37.00%—45.63%之間,而凋落物呼吸貢獻(xiàn)率則在18.62%—23.70%之間。其中,除灌措施主要是降低了礦質(zhì)土壤呼吸RM和根呼吸RR,采伐2措施則主要是提高了礦質(zhì)土壤呼吸RM,而對根呼吸沒有產(chǎn)生明顯影響,采伐1措施雖然也顯著提高了礦質(zhì)土壤呼吸RM,但又同時(shí)顯著降低了根呼吸RR。因此,不同營林措施后森林土壤總呼吸Rtot差異顯著,主要是由于不同營林措施對礦質(zhì)土壤呼吸RM和根呼吸RR的影響存在差異造成的。

    表7 不同營林措施土壤溫度(T)、土壤濕度(W)與土壤呼吸組分的復(fù)合模型參數(shù)

    Table 7 Composite model parameters derived for the combined influence of soil temperature (T) and moisture (W) on the soil respiration components from different forest management

    營林措施Forestmanagement參數(shù)Parameter凋落物層呼吸速率RLLitterlayerRespiration根系呼吸速率RRRootrespiration礦質(zhì)土壤呼吸速率RMMineralsoilrespiration數(shù)值Value標(biāo)準(zhǔn)誤SER2P數(shù)值Value標(biāo)準(zhǔn)誤SER2P數(shù)值Value標(biāo)準(zhǔn)誤SER2P對照(CK)c1-0.6910.6760.271<0.01-0.6470.4490.368<0.01-1.7090.5120.243<0.01Controlc2-0.0260.0410.0340.0270.1290.035c3-0.1300.055-0.0560.0360.0120.032c40.0090.0030.0040.002-0.0020.002除灌(SM)c1-1.3570.3910.0460.119-1.8440.3970.358<0.01-0.8650.4710.200<0.01Shrubs-cuttingc20.0200.0260.0780.0260.0110.032c30.0050.0360.0230.036-0.0440.034c40.0000.0020.0000.0020.0030.002采伐1(H1)c1-3.2340.7220.185<0.01-1.4270.4290.211<0.01-2.0560.4250.553<0.01Harvest1c20.1240.0420.0680.0250.1010.029c30.0890.0540.0350.0320.0390.026c4-0.0040.003-0.0020.002-0.0010.002采伐2(H2)c1-0.8121.0110.107<0.01-1.2730.6260.502-1.5190.7130.593<0.01Harvest2c2-0.0380.0470.0740.0290.0870.030c3-0.0380.047-0.0110.029-0.0050.032c40.0040.0020.0000.0010.0000.001

    不同營林措施并未造成凋落物層呼吸速率RL的顯著差異,這與相關(guān)研究不同。孟春等研究表明擇伐后4年針闊混交林生長季RL呈逐年增加的趨勢,且與采伐強(qiáng)度呈正相關(guān)[21];劉可等發(fā)現(xiàn)不同擇伐強(qiáng)度的油松幼齡林人工林生長季RL顯著高于對照,但RL值與采伐強(qiáng)度無相關(guān)關(guān)系[22]。上述研究均認(rèn)凋落物的分解與土壤溫度、濕度顯著相關(guān)[23],采伐造成了土壤溫度的顯著上升和濕度的顯著下降,從而導(dǎo)致凋落物層呼吸速率的顯著變化。而本研究中采伐1和除灌干擾強(qiáng)度小,并未造成土壤溫度和濕度的顯著變異性,而干擾程度較大的采伐2,土壤在溫度升高的同時(shí),濕度也顯著上升,RL無明顯差異可能是土壤溫度和水分綜合作用的結(jié)果(表4),此外,研究結(jié)果的不同也可能源于不同的森林類型下凋落物的質(zhì)量的差異,針葉林分分解速率低于針闊混交林[24],較短的觀測期時(shí)間(1年)內(nèi),可能RL的差異性并不顯著。

    營林措施對于根呼吸速率RR的影響可能是正效應(yīng)也可能是負(fù)效應(yīng)[20],采伐后地上植被的減少會直接降低地下的活根生物量進(jìn)而造成RR的減少,但地上植被減少后,充足的養(yǎng)分會促進(jìn)剩余植被的生長,導(dǎo)致活根生物量的增加進(jìn)而造成RR的增加[6]。本研究中H2與CK的RR無顯著差異,可能是上述正負(fù)效應(yīng)綜合作用的結(jié)果。除灌和采伐1措施主要是清除了灌木和林下樹種,這些植被的生物量占整個林分的比重較低,清除后對剩余植被養(yǎng)分利用的影響有限,但這些植被細(xì)根比重較大,而且分布較淺,清除后可能會顯著降低林分表層土壤根系生物量,尤其是灌木,這是否是除灌和采伐1措施下根呼吸RR顯著低于CK原因,仍需進(jìn)一步進(jìn)行調(diào)查證實(shí)。

    不同的營林措施造成了礦質(zhì)土壤呼吸RM顯著的差異性(p<0.005),RM從大到小依次為H2(1.27±0.07)μmolCO2m-2s-1>H1(1.17±0.07)μmolCO2m-2s-1>CK(0.95±0.06)μmolCO2m-2s-1>SM(0.75±0.05)μmolCO2m-2s-1。RM的變化與土壤微生物的活動和和底物的供應(yīng)有關(guān)。采伐后,遺留在林地上的伐根死亡并開始分解,為礦質(zhì)土壤的呼吸增加了底物;同時(shí),地上活立木的減少,改變了林內(nèi)微環(huán)境,為土壤微生物活動創(chuàng)造了適宜的條件,使得礦質(zhì)土壤呼吸增加[25-26]。不同于采伐措施,除灌措施對剩余物進(jìn)行了仔細(xì)的清除,地表直接暴露在空氣中,缺乏灌木對于陽光和雨水的調(diào)節(jié),但監(jiān)測結(jié)果表明,在除灌的第1年內(nèi)地表年平均溫濕度與CK沒有顯著差異(表4),造成除灌后RM顯著降低的原因可能是土壤微生物生物量或群落結(jié)構(gòu)發(fā)生了變化,這些需要下一步深入研究進(jìn)行證實(shí)。不同營林措施下,RM對土壤呼吸的貢獻(xiàn)率較高,且采伐措施下RM的貢獻(xiàn)率顯著高于對照和除灌,RM的差異是造成采伐與對照土壤呼吸速率差異的主要原因。

    4.2 土壤呼吸組分與土壤溫濕度

    多數(shù)研究表明,土壤呼吸與土壤溫度呈現(xiàn)顯著地相關(guān)關(guān)系[10-11, 27-28],其季節(jié)變化與氣溫、物候、植物的生理活動、微生物活性等密切相關(guān)[29]。本研究中,土壤呼吸各組分均與土壤溫度有著顯著地正相關(guān)關(guān)系。不同營林措施下土壤溫度月均值變化呈單峰趨勢于7月達(dá)到最大值,RM、RR的月均值與土壤溫度變化呈現(xiàn)相同的趨勢,均在生長旺盛期7—8月達(dá)到最大值,與相關(guān)研究一致[22]。而RL則呈現(xiàn)出一定的波動性,在生長旺盛期6—9月間達(dá)到最大值。營林措施從不同程度改變了凋落物層對于土壤呼吸底物的C供應(yīng),使得土壤溫度對不同營林措施RL的變化解釋度降低(表5)。與RL相比,土壤溫度能更好的解釋RR的變化,是造成H2的總呼吸高于其他措施最主要的原因。

    在觀測期內(nèi),土壤濕度并未表現(xiàn)出明顯的季節(jié)變化,不同營林措施下,H2造成了土壤濕度均值的顯著增加,H1和SM的土壤濕度均值略高于CK(表3)。同時(shí),土壤濕度僅能解釋部分土壤呼吸組分的變化(表6)。土壤濕度會影響土壤中溶解性有機(jī)質(zhì)的變化、植物和微生物能量的分配,并且對植物和微生物能量的分配產(chǎn)生作用[30],因此對土壤呼吸的影響更為復(fù)雜。本研究結(jié)果與相關(guān)研究一致。如太原天龍山自然保護(hù)區(qū)油松林4年土壤呼吸觀測研究表明,土壤濕度僅能解釋土壤呼吸的2%—25%[31];與本研究區(qū)域較接近的神農(nóng)架不同海拔梯度上4種森林類型的土壤呼吸研究表明土壤呼吸與土壤濕度并無顯著地相關(guān)關(guān)系[32]。

    相較于單一因子,土壤溫度與土壤濕度的雙因子模型能更好的說明不同營林措施下土壤呼吸組分的變化[7, 10],如土壤溫濕度共同解釋了東北小興安嶺地區(qū)擇伐后生長季針闊混交林土壤組分呼吸調(diào)變化的55%—99%[21],山西太岳山擇伐后幼齡林人工油松林土壤呼吸組分變化的46.4%—99.3%[22],但本研究中剔除解釋率最低的RL外,土壤溫濕度的僅能解釋土壤呼吸組分變化的20.0%—59.3%(表7),鑒于目前有限的研究,并未能找到造成該本文結(jié)果與其他研究結(jié)論差異的具體原因,仍需進(jìn)行進(jìn)一步的調(diào)查分析。但本結(jié)果在一定程度上反映出土壤溫濕度并不是營林后土壤呼吸組分變化的主要原因,應(yīng)對其他因子,如光照、降雨、微生物組成、活性、植被生長狀況等進(jìn)行全面的調(diào)查,同時(shí)加強(qiáng)對營林后土壤呼吸組分的長期觀測,短期內(nèi)僅能觀測到不同呼吸組分的初步動態(tài)變化,長期連續(xù)的觀測才能找到營林后土壤呼吸組分變化的機(jī)制。

    5 結(jié)論

    3種營林措施對森林土壤呼吸組分的影響差異明顯。不同營林措施下RL無顯著變化;除灌措施顯著降低了林分的根呼吸RR和礦質(zhì)土壤呼吸RM;采伐1措施雖然也顯著提高了礦質(zhì)土壤呼吸RM,但又同時(shí)顯著降低了根呼吸RR;措施采伐2措施則主要是提高了礦質(zhì)土壤呼吸RM,而對跟呼吸沒有產(chǎn)生明顯影響。

    同時(shí)本研究對不同經(jīng)營措施后土壤溫度和濕度因子進(jìn)行了探討,但土壤溫濕度僅能解釋土壤呼吸主要組分RR和RM的部分變化(20%—59.3%),仍需對其它相關(guān)因子如降雨、光照、微生物組成結(jié)構(gòu)等進(jìn)行深入的研究,同時(shí)對營林后土壤呼吸速率及其組成進(jìn)行長期連續(xù)的監(jiān)測以確定營林干擾對不同土壤組分的影響機(jī)制。

    總之,在營林措施實(shí)施后至第1個生長季結(jié)束時(shí),除灌措施使森林土壤總呼吸Rtot降低了17.74%,采伐2措施使森林土壤總呼吸Rtot提高了17.34%,而采伐1措施對森林土壤總呼吸Rtot沒有產(chǎn)生明顯的影響。針對三峽庫區(qū)大面的人工飛播林,經(jīng)過30—40a的自然生長后,定期對林下植被進(jìn)行清理將有利于減少降低林分土壤呼吸,而對部分成熟林的采伐利用,在短期內(nèi)會增加林分的土壤呼吸。

    致謝:感謝湖北秭歸三峽庫區(qū)森林生態(tài)系統(tǒng)國家定位觀測研究站對本研究的支持,感謝紐芬蘭紀(jì)念大學(xué)欒軍偉博士對本文寫作的幫助。

    [1] Pan Y D, Birdsey R A, Fang J Y, Houghton R, Kauppi P E, Kurz W A, Phillips O L, Shvidenko A, Lewis S L, Canadell J G, Ciais P, Jackson R B, Pacala S W, McGuire A D, Piao S L, Rautiainen A, Sitch S, Hayes D. A large and persistent carbon sink in the world′s forests. Science, 2011, 333(6045): 988- 993.

    [2] Williams C A, Vanderhoof M K, Khomik M, Ghimire B. Post-clearcut dynamics of carbon, water and energy exchanges in a midlatitude temperate, deciduous broadleaf forest environment. Global Change Biology, 2013, 20(3): 992- 1007.

    [3] Dixon R K,, Brown S, Houghton R ASolomom A M, Trexler M C, Wisniewski J. Carbon pools and flux of global forest ecosystems. Science, 1994, 263(5144): 185- 190.

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

    [5] Jandl R, Lindner M, Vesterdal L, Bauwens B, Baritz R, Hagedorn F, Johnson D W, Minkkinen K, Byrne K A. How strongly can forest management influence soil carbon sequestration?. Geoderma, 2007, 137(3/4): 253- 268.

    [6] Ma S Y, Chen J Q, North M, Erickson H E, Bresee M, Le Moine J. Short-term effects of experimental burning and thinning on soil respiration in an old-growth, mixed-conifer forest. Environmental Management, 2004, 33(S1): S148-S159.

    [7] Tang J W, Qi Y, Xu M, Misson L, Goldstein A H. Forest thinning and soil respiration in a ponderosa pine plantation in the Sierra Nevada. Tree Physiology, 2005, 25(1): 57- 66.

    [8] Saunders M, Tobin B, Black K, Gioria M, Nieuwenhuis M, Osborne B A. Thinning effects on the net ecosystem carbon exchange of a Sitka spruce forest are temperature-dependent. Agricultural and Forest Meteorology, 2012, 157: 1- 10.

    [9] Ryu S R, Concilio A, Chen J Q, North M, Ma S Y. Prescribed burning and mechanical thinning effects on belowground conditions and soil respiration in a mixed-conifer forest, California. Forest Ecology and Management, 2009, 257(4): 1324- 1332.

    [10] 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.

    [11] Concilio A, Ma S Y, Li Q L, LeMoine J, Chen J Q, North M, Moorhead D, Jensen R. Soil respiration response to prescribed burning and thinning in mixed-conifer and hardwood forests. Canadian Journal of Forest Research, 2005, 35(7): 1581- 1591.

    [12] Laporte M F, Duchesne L C, Morrison I K. Effect of clearcutting, selection cutting, shelterwood cutting and microsites on soil surface CO2efflux in a tolerant hardwood ecosystem of northern Ontario. Forest Ecology and Management, 2003, 174(1/3): 565- 575.

    [13] Fleming R L, Laporte M F, Hogan G D, Hazlett P W. Effects of harvesting and soil disturbance on soil CO2efflux from a jack pine forest. Canadian Journal of Forest Research, 2006, 36(3): 589- 600.

    [14] Bowden R D, Nadelhoffer K J, Boone R D, Melillo J M, Garrison J B. Contributions of aboveground litter, belowground litter, and root respiration to total soil respiration in a temperate mixed hardwood forest. Canadian Journal of Forest Research, 1993, 23(7): 1402- 1407.

    [15] Lee N Y, Koo J W, Noh N J, Kim J, Son Y. Autotrophic and heterotrophic respiration in needle fir andQuercus-dominated stands in a cool-temperate forest, central Korea. Journal of Plant Research, 2010, 123(4): 485- 495.

    [16] 肖文發(fā), 雷靜品. 三峽庫區(qū)森林植被恢復(fù)與可持續(xù)經(jīng)營研究. 長江流域資源與環(huán)境, 2004, 13(2): 138- 144.

    [17] 葛曉改, 肖文發(fā), 曾立雄, 黃志霖, 黃玲玲, 譚本旺. 三峽庫區(qū)不同林齡馬尾松土壤養(yǎng)分與酶活性的關(guān)系. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(2): 445- 451.

    [18] Luan J W, Liu S R, Zhu X L, Wang J X, Liu K. Roles of biotic and abiotic variables in determining spatial variation of soil respiration in secondary oak and planted pine forests. Soil Biology and Biochemistry, 2012, 44(1): 143- 150.

    [19] 褚金翔, 張小全. 川西亞高山林區(qū)三種土地利用方式下土壤呼吸動態(tài)及組分區(qū)分. 生態(tài)學(xué)報(bào), 2006, 26(6): 1693- 1700.

    [20] Peng Y Y, Thomas S C, Tian D. Forest management and soil respiration: Implications for carbon sequestration. Environmental Reviews, 2008, 16(NA): 93- 111.

    [21] 孟春, 王立海, 沈微. 擇伐對生長季針闊混交林土壤分室呼吸的影響. 林業(yè)科學(xué), 2008, 44(8): 23- 28.

    [22] 劉可, 韓海榮, 康峰峰, 程小琴, 宋婭麗, 周彬, 李勇. 山西太岳山油松人工林生長季土壤呼吸對擇伐強(qiáng)度的響應(yīng). 生態(tài)學(xué)雜志, 2013, 32(12): 3173- 3181.

    [23] 葛曉改, 曾立雄, 黃志霖, 肖文發(fā), 譚本旺, 周本智. 土壤溫度和水分含量對三峽庫區(qū)馬尾松林凋落物葉分解的影響. 林業(yè)科學(xué), 2013, 49(9): 153- 157.

    [24] 李志安, 鄒碧, 丁永禎, 曹裕松. 森林凋落物分解重要影響因子及其研究進(jìn)展. 生態(tài)學(xué)雜志, 2004, 23(6): 77- 83.

    [25] Epron D, Ngao J, Granier A. Interannual variation of soil respiration in a beech forest ecosystem over a six-year study. Annuals of Forest Science, 2004, 61(6): 499- 505.

    [26] 孟春, 王立海, 沈微. 擇伐對小興安嶺地區(qū)針闊混交林土壤呼吸溫度敏感性的影響. 林業(yè)科學(xué), 2011, 47(3): 102- 106.

    [27] 劉穎, 韓士杰, 胡艷玲, 戴冠華. 土壤溫度和濕度對長白松林土壤呼吸速率的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2005, 16(9): 1581- 1585.

    [28] 朱凡, 王光軍, 田大倫, 閆文德, 項(xiàng)文化, 梁小翠. 馬尾松人工林根呼吸的季節(jié)變化及影響因子. 林業(yè)科學(xué), 2010, 46(7): 36- 41.

    [29] Ngao J, Epron D, Delpierre N, Bréda N, Granier A, Longdoz B. Spatial variability of soil CO2efflux linked to soil parameters and ecosystem characteristics in a temperate beech forest. Agricultural and Forest Meteorology, 2012, 154- 155: 136- 146.

    [30] 陳全勝, 李凌浩, 韓興國, 閻志丹. 水分對土壤呼吸的影響及機(jī)理. 生態(tài)學(xué)報(bào), 2003, 23(5): 972- 978.

    [31] 嚴(yán)俊霞, 秦作棟, 張義輝, 李洪建. 土壤溫度和水分對油松林土壤呼吸的影響. 生態(tài)學(xué)報(bào), 2009, 29(12): 6366- 6376.

    [32] 羅璐, 申國珍, 謝宗強(qiáng), 周利光. 神農(nóng)架海拔梯度上4種典型森林的土壤呼吸組分及其對溫度的敏感性. 植物生態(tài)學(xué)報(bào), 2011, 35(7): 722- 730.

    Responses of soil respiration and its components to forest management in Pinus massoniana stands

    LEI Lei1,2, XIAO Wenfa1,2, ZENG Lixiong1,2,*, HUANG Zhilin1,2, TAN Benwang3

    1 Research Institute of Forest Ecology, Environment and Protection, Chinese Academy of Forestry, State Forestry Administration Key Laboratory of Forest Ecology and Environment, Beijing 100091, China 2Co-InnovationCenterforSustainableForestryinSouthernChina,NanjingForestryUniversity,Nanjing210037,China3ForestryBureauofZiguiCounty,Yichang443631,China

    Forest soils play a critical role in the carbon cycle and carbon sequestration at both global and local scales, and forest management practices (e.g., harvesting, burning, and thinning) influence soil carbon processes by altering organic matter quality and quantity, key microclimatic conditions, microbial communities, and other factors. However, the effects of forest management on soil carbon effluxes in different ecosystems are still largely unknown, despite being critical to estimate global carbon fluxes. Quantifying the responses of soil respiration and its components to forest management is vital to accurately evaluate forest carbon balance. Thus, an aerially seededPinusmassonianaforest was chosen in the Three Gorges reservoir area to evaluate the effects on soil respiration of different forest management practices (i.e., control; shrub-cutting: harvesting all shrubs and removing all harvest residues. Harvest strategy 1%—15% harvest intensity and removing main harvest residues without leaves and small branches; Harvest strategy 2%—70% harvest intensity and the same harvest residue management as that in harvest 1). All experimental treatments were located in similar habitats and consisted of three 20 m × 20 m plots. The treatments were conducted in October 2013. A combination of trenching and litter removing methods were employed in order to partition soil respiration into components of litter layer respiration, root respiration, and mineral soil respiration. The soil temperature, soil moisture, and rate of soil respiration and its components were observed continuously for one year (from November 2013 to October 2014) using a Li- 8100 system. Management did not affect litter layer respiration within the measuring period. Root respiration in the control, shrub-cutting, harvest strategy 1, and harvest strategy 2 treatments were 1.00, 0.83, 0.86, and 1.11 μmol CO2m-2s-1, respectively. The mineral soil respiration of harvested stands was significantly higher than that of control and shrub-cutting stands (P< 0.05). The proportion of litter layer respiration to total respiration was not significantly influenced by forest management (18.78%—23.70%), but the contribution of root respiration to total respiration was reduced, especially in the harvest strategy 1 treatment (P< 0.05). The contribution of mineral respiration to total respiration significantly increased in harvested stands, which was mainly attributed to the decrease in total soil respiration. Management effects on soil temperature and soil moisture were only observed in the stands of harvest strategy 2. A two-factor model that included soil temperature and moisture better explained the variations in soil respiration (4.6%—59.3%) than that by a model using temperature (4.2%—59.1%) or moisture (0.3%—23.5%) alone. Other factors that influence soil respiration and its components need to be further elucidated.

    forest management; soil respiration components; litter layer respiration; root respiration; mineral soil respiration

    中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(CAFRIFEEP201101);林業(yè)公益性行業(yè)科研專項(xiàng)(201104008)

    2015- 03- 02;

    日期:2015- 12- 14

    10.5846/stxb201503020398

    *通訊作者Corresponding author.E-mail: zlxcaf@163.com

    雷蕾,肖文發(fā),曾立雄,黃志霖,譚本旺.馬尾松林土壤呼吸組分對不同營林措施的響應(yīng).生態(tài)學(xué)報(bào),2016,36(17):5360- 5370.

    Lei L, Xiao W F, Zeng L X, Huang Z L, Tan B W.Responses of soil respiration and its components to forest management inPinusmassonianastands.Acta Ecologica Sinica,2016,36(17):5360- 5370.

    猜你喜歡
    森林土壤礦質(zhì)土壤濕度
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    基于51單片機(jī)控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    不同產(chǎn)地三葉青中27種礦質(zhì)元素的綜合評價(jià)
    中成藥(2018年11期)2018-11-24 02:57:08
    四川盆地土壤濕度時(shí)空分布及影響因子分析
    中國不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    兩種亞熱帶森林土壤甲烷氧化活性的垂直分布特征
    基于礦質(zhì)元素的絞股藍(lán)產(chǎn)地鑒別研究
    龍灣自然保護(hù)區(qū)森林土壤動物群落多樣性及功能類群
    中亞熱帶地區(qū)森林土壤真菌生物量及影響因子研究
    縉云山森林土壤微生物數(shù)量與群落特征
    久久婷婷青草| 国产精品国产三级国产av玫瑰| 2022亚洲国产成人精品| 美女主播在线视频| 中文字幕人妻熟女乱码| 日本色播在线视频| 亚洲精品国产av成人精品| av不卡在线播放| 久久韩国三级中文字幕| 久久青草综合色| av国产精品久久久久影院| 日日撸夜夜添| 春色校园在线视频观看| 九九在线视频观看精品| 日本黄色日本黄色录像| 免费看光身美女| 大香蕉久久网| 国产男女超爽视频在线观看| 女人久久www免费人成看片| 日韩熟女老妇一区二区性免费视频| 国产白丝娇喘喷水9色精品| 精品国产一区二区三区久久久樱花| 满18在线观看网站| 久久精品国产自在天天线| 中文字幕最新亚洲高清| 精品卡一卡二卡四卡免费| 97人妻天天添夜夜摸| 在线 av 中文字幕| 午夜久久久在线观看| 亚洲精品色激情综合| 在线观看国产h片| 午夜久久久在线观看| 国产1区2区3区精品| 侵犯人妻中文字幕一二三四区| 欧美成人午夜精品| 9热在线视频观看99| 国产精品99久久99久久久不卡 | 国产亚洲av片在线观看秒播厂| 一级毛片我不卡| 男女高潮啪啪啪动态图| 国产亚洲欧美精品永久| 看非洲黑人一级黄片| 亚洲精品一二三| 日本欧美视频一区| 久久精品久久久久久久性| 飞空精品影院首页| 18禁国产床啪视频网站| 自线自在国产av| 亚洲五月色婷婷综合| 丰满少妇做爰视频| 日韩免费高清中文字幕av| 国产成人精品婷婷| 哪个播放器可以免费观看大片| 午夜视频国产福利| 亚洲av综合色区一区| 天美传媒精品一区二区| 色吧在线观看| 免费观看av网站的网址| 日本黄大片高清| 日日啪夜夜爽| av不卡在线播放| 99久国产av精品国产电影| 国产在视频线精品| 欧美成人午夜精品| 日韩 亚洲 欧美在线| 麻豆精品久久久久久蜜桃| 久久久久网色| 一二三四在线观看免费中文在 | 亚洲欧美成人精品一区二区| 亚洲成av片中文字幕在线观看 | 免费少妇av软件| 久久久欧美国产精品| 五月开心婷婷网| 亚洲欧美一区二区三区黑人 | 精品久久久久久电影网| 一级片'在线观看视频| 婷婷色综合www| 国产日韩欧美视频二区| 97在线人人人人妻| 久久av网站| 国产精品一区二区在线观看99| 日韩av不卡免费在线播放| 精品一区二区三区四区五区乱码 | 亚洲情色 制服丝袜| 美女国产高潮福利片在线看| 飞空精品影院首页| 日本欧美视频一区| 伦理电影大哥的女人| 在线天堂最新版资源| av在线观看视频网站免费| 看非洲黑人一级黄片| 日本猛色少妇xxxxx猛交久久| 少妇被粗大猛烈的视频| 爱豆传媒免费全集在线观看| 国产一区二区三区综合在线观看 | 夫妻性生交免费视频一级片| 一本大道久久a久久精品| 毛片一级片免费看久久久久| 成年人午夜在线观看视频| 黄色 视频免费看| 午夜福利乱码中文字幕| 国产精品成人在线| 国精品久久久久久国模美| 三上悠亚av全集在线观看| 久久女婷五月综合色啪小说| 亚洲av.av天堂| 国产精品久久久久成人av| 91成人精品电影| 男女免费视频国产| av天堂久久9| 亚洲成人一二三区av| 深夜精品福利| 亚洲av男天堂| 免费黄色在线免费观看| 黄色配什么色好看| av又黄又爽大尺度在线免费看| 久久久久精品性色| 不卡视频在线观看欧美| 中国三级夫妇交换| 亚洲欧美中文字幕日韩二区| 免费看光身美女| 亚洲国产av新网站| 97人妻天天添夜夜摸| 中国美白少妇内射xxxbb| 日本av免费视频播放| 一级,二级,三级黄色视频| 国产成人精品婷婷| 亚洲av综合色区一区| 五月开心婷婷网| 十分钟在线观看高清视频www| av不卡在线播放| 欧美变态另类bdsm刘玥| 69精品国产乱码久久久| 午夜老司机福利剧场| 国产精品久久久av美女十八| 免费人妻精品一区二区三区视频| 久久午夜福利片| 久久久久久久亚洲中文字幕| 国产黄色免费在线视频| 午夜91福利影院| 美女内射精品一级片tv| 国产午夜精品一二区理论片| 亚洲伊人久久精品综合| 51国产日韩欧美| 高清视频免费观看一区二区| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 不卡视频在线观看欧美| 又粗又硬又长又爽又黄的视频| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 五月伊人婷婷丁香| 国产精品久久久久久精品电影小说| 宅男免费午夜| 99久久中文字幕三级久久日本| 9热在线视频观看99| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码 | h视频一区二区三区| 在线观看免费日韩欧美大片| 丝袜美足系列| 18禁国产床啪视频网站| 精品国产一区二区久久| 国产老妇伦熟女老妇高清| 18禁国产床啪视频网站| 一本色道久久久久久精品综合| 国产精品99久久99久久久不卡 | 久久这里有精品视频免费| 黄片无遮挡物在线观看| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 少妇人妻久久综合中文| 韩国av在线不卡| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 久久久久久久大尺度免费视频| 色网站视频免费| 九色成人免费人妻av| 热re99久久精品国产66热6| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 欧美成人午夜免费资源| 亚洲精品久久久久久婷婷小说| 啦啦啦视频在线资源免费观看| 欧美精品高潮呻吟av久久| 日产精品乱码卡一卡2卡三| 高清不卡的av网站| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 熟女人妻精品中文字幕| 亚洲经典国产精华液单| 男女国产视频网站| 少妇人妻精品综合一区二区| 女人久久www免费人成看片| 少妇的逼好多水| 蜜桃国产av成人99| 日韩人妻精品一区2区三区| 日本av手机在线免费观看| 在线观看人妻少妇| 久久久久久久久久久久大奶| 人人妻人人添人人爽欧美一区卜| 在线观看免费视频网站a站| 婷婷色综合大香蕉| 韩国av在线不卡| 国产男女超爽视频在线观看| 久久久久久久精品精品| 精品一区二区三卡| 寂寞人妻少妇视频99o| 久久热在线av| 免费人成在线观看视频色| av免费观看日本| 国内精品宾馆在线| 精品人妻熟女毛片av久久网站| 亚洲国产欧美在线一区| 在线观看www视频免费| 亚洲av电影在线进入| 久久久精品94久久精品| 巨乳人妻的诱惑在线观看| 一区在线观看完整版| 老司机影院成人| 精品人妻熟女毛片av久久网站| xxx大片免费视频| 国产精品蜜桃在线观看| 国产精品 国内视频| 免费观看无遮挡的男女| 熟女av电影| 女人精品久久久久毛片| 久久久久久久国产电影| 欧美日韩成人在线一区二区| 国产精品一区www在线观看| 国产免费现黄频在线看| 最近中文字幕高清免费大全6| 久久国产精品男人的天堂亚洲 | 老女人水多毛片| 视频在线观看一区二区三区| 男人爽女人下面视频在线观看| 九草在线视频观看| 亚洲国产精品一区三区| 一区二区日韩欧美中文字幕 | 国产精品久久久久久精品古装| 国产成人精品久久久久久| 日本av手机在线免费观看| 男女无遮挡免费网站观看| 亚洲av成人精品一二三区| 99热6这里只有精品| 女人久久www免费人成看片| 91久久精品国产一区二区三区| 日日啪夜夜爽| 99热国产这里只有精品6| 99视频精品全部免费 在线| 有码 亚洲区| 久久精品国产亚洲av天美| 最新的欧美精品一区二区| 午夜老司机福利剧场| 国产探花极品一区二区| 91久久精品国产一区二区三区| 午夜福利视频在线观看免费| 国产精品秋霞免费鲁丝片| 欧美精品亚洲一区二区| 亚洲av日韩在线播放| 91精品国产国语对白视频| 一区二区日韩欧美中文字幕 | 18禁动态无遮挡网站| 亚洲国产精品一区二区三区在线| 91精品国产国语对白视频| 国产精品人妻久久久久久| 国产一区亚洲一区在线观看| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 午夜激情av网站| 99精国产麻豆久久婷婷| 91精品伊人久久大香线蕉| 日韩成人伦理影院| 欧美精品av麻豆av| 午夜免费观看性视频| 中国三级夫妇交换| 成年av动漫网址| 亚洲精品,欧美精品| 亚洲精品日本国产第一区| 人人妻人人澡人人爽人人夜夜| 精品第一国产精品| 亚洲国产看品久久| 日本vs欧美在线观看视频| 亚洲精品日本国产第一区| 夫妻性生交免费视频一级片| av卡一久久| 久久精品国产亚洲av涩爱| 国产国拍精品亚洲av在线观看| 成年av动漫网址| 街头女战士在线观看网站| 久久久久网色| 成人毛片60女人毛片免费| 午夜91福利影院| 免费大片黄手机在线观看| 亚洲av电影在线观看一区二区三区| 亚洲欧美色中文字幕在线| 精品国产一区二区久久| 黑人巨大精品欧美一区二区蜜桃 | 国产亚洲av片在线观看秒播厂| 最黄视频免费看| 韩国av在线不卡| av网站免费在线观看视频| 美国免费a级毛片| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 久久狼人影院| 高清av免费在线| av黄色大香蕉| av免费在线看不卡| 国产日韩欧美亚洲二区| 国产成人91sexporn| 久久久久久久大尺度免费视频| 只有这里有精品99| 黄网站色视频无遮挡免费观看| 亚洲av欧美aⅴ国产| 高清欧美精品videossex| 亚洲经典国产精华液单| 大片免费播放器 马上看| 不卡视频在线观看欧美| 久久青草综合色| 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 日韩欧美一区视频在线观看| 国产亚洲欧美精品永久| 在线观看免费视频网站a站| 五月玫瑰六月丁香| 日本欧美视频一区| 一本色道久久久久久精品综合| 伊人久久国产一区二区| 国产成人91sexporn| 内地一区二区视频在线| 黄片播放在线免费| 日韩成人伦理影院| 综合色丁香网| 十八禁网站网址无遮挡| 91久久精品国产一区二区三区| 午夜福利乱码中文字幕| 一级毛片电影观看| 欧美97在线视频| 成人国产av品久久久| 美女内射精品一级片tv| 国产男女超爽视频在线观看| 国语对白做爰xxxⅹ性视频网站| 久久久国产精品麻豆| 看免费av毛片| 女性生殖器流出的白浆| 久久午夜福利片| 亚洲经典国产精华液单| 一级毛片我不卡| 黑丝袜美女国产一区| 免费大片18禁| 在线天堂中文资源库| 最近中文字幕高清免费大全6| 深夜精品福利| 亚洲色图综合在线观看| 国产成人精品福利久久| 男男h啪啪无遮挡| 亚洲,欧美精品.| 亚洲国产精品一区三区| 精品少妇黑人巨大在线播放| 看免费av毛片| 久久久国产精品麻豆| 亚洲成人手机| 一级毛片黄色毛片免费观看视频| 美国免费a级毛片| 国产有黄有色有爽视频| 欧美国产精品va在线观看不卡| 一级毛片黄色毛片免费观看视频| 国产精品熟女久久久久浪| 国产成人精品久久久久久| 日本欧美国产在线视频| 在线看a的网站| 老司机亚洲免费影院| 亚洲第一区二区三区不卡| 夜夜骑夜夜射夜夜干| 精品国产露脸久久av麻豆| 热re99久久国产66热| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| 成年人免费黄色播放视频| 中国国产av一级| 宅男免费午夜| 欧美性感艳星| 国产在线视频一区二区| 啦啦啦视频在线资源免费观看| 精品酒店卫生间| 国产综合精华液| 亚洲人成网站在线观看播放| 黑人欧美特级aaaaaa片| 99国产精品免费福利视频| 久久精品久久精品一区二区三区| 久久鲁丝午夜福利片| 亚洲欧美日韩卡通动漫| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 你懂的网址亚洲精品在线观看| 中文字幕精品免费在线观看视频 | 亚洲精品成人av观看孕妇| 免费高清在线观看视频在线观看| 蜜桃国产av成人99| 一区二区av电影网| videos熟女内射| 国产国语露脸激情在线看| 亚洲美女搞黄在线观看| 亚洲成国产人片在线观看| 在线观看一区二区三区激情| 有码 亚洲区| 国产又色又爽无遮挡免| 视频中文字幕在线观看| 色94色欧美一区二区| 性色avwww在线观看| 日本免费在线观看一区| 欧美日韩av久久| 男女啪啪激烈高潮av片| 啦啦啦啦在线视频资源| 欧美bdsm另类| 国产日韩欧美亚洲二区| 日韩av免费高清视频| 午夜免费男女啪啪视频观看| 看免费av毛片| 国产av精品麻豆| 欧美国产精品va在线观看不卡| 亚洲av欧美aⅴ国产| 欧美成人午夜精品| 五月伊人婷婷丁香| videos熟女内射| 亚洲av福利一区| 欧美变态另类bdsm刘玥| 看免费av毛片| 精品国产露脸久久av麻豆| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 香蕉精品网在线| 搡女人真爽免费视频火全软件| 新久久久久国产一级毛片| 五月玫瑰六月丁香| 18禁国产床啪视频网站| 久久国产亚洲av麻豆专区| 欧美最新免费一区二区三区| 色婷婷久久久亚洲欧美| 大香蕉久久成人网| 久久人人97超碰香蕉20202| 成人国产av品久久久| 国产日韩欧美视频二区| 最新中文字幕久久久久| 国产一区有黄有色的免费视频| 边亲边吃奶的免费视频| 欧美3d第一页| 日本vs欧美在线观看视频| 国产免费又黄又爽又色| 国产亚洲最大av| 精品少妇黑人巨大在线播放| 欧美日韩国产mv在线观看视频| 亚洲国产精品一区二区三区在线| a级毛色黄片| 中文字幕人妻熟女乱码| videos熟女内射| 有码 亚洲区| 欧美人与善性xxx| 美女内射精品一级片tv| 中文天堂在线官网| 久热久热在线精品观看| 久久久久久伊人网av| 亚洲丝袜综合中文字幕| 综合色丁香网| 国产av码专区亚洲av| av电影中文网址| 国产永久视频网站| 亚洲国产成人一精品久久久| 日日撸夜夜添| 赤兔流量卡办理| 欧美日韩视频高清一区二区三区二| 一区二区日韩欧美中文字幕 | 亚洲精品久久成人aⅴ小说| 97在线视频观看| 色5月婷婷丁香| 少妇猛男粗大的猛烈进出视频| av.在线天堂| 久久人人爽人人片av| av福利片在线| 欧美精品高潮呻吟av久久| 亚洲精品国产色婷婷电影| 欧美激情极品国产一区二区三区 | av在线观看视频网站免费| 狠狠精品人妻久久久久久综合| av免费在线看不卡| 日韩在线高清观看一区二区三区| 嫩草影院入口| 日韩一本色道免费dvd| 女的被弄到高潮叫床怎么办| 成人国语在线视频| 国产又色又爽无遮挡免| 精品少妇内射三级| 人成视频在线观看免费观看| 国产69精品久久久久777片| 国产成人91sexporn| 午夜视频国产福利| 卡戴珊不雅视频在线播放| 精品熟女少妇av免费看| 人妻系列 视频| 日韩精品有码人妻一区| 国产淫语在线视频| 最近中文字幕高清免费大全6| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美日韩另类电影网站| 最新中文字幕久久久久| 日韩一区二区视频免费看| 久久99热这里只频精品6学生| 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 亚洲欧洲国产日韩| 亚洲精品久久成人aⅴ小说| av在线播放精品| 最近最新中文字幕大全免费视频 | 日本wwww免费看| 久久这里有精品视频免费| 亚洲欧美日韩卡通动漫| 交换朋友夫妻互换小说| 18禁观看日本| 多毛熟女@视频| 女性生殖器流出的白浆| 亚洲欧美一区二区三区国产| 人人妻人人澡人人看| 极品少妇高潮喷水抽搐| 国产欧美另类精品又又久久亚洲欧美| 成年女人在线观看亚洲视频| 亚洲,一卡二卡三卡| 中文字幕人妻丝袜制服| 男人爽女人下面视频在线观看| av网站免费在线观看视频| 国产精品久久久久久久久免| 国产亚洲欧美精品永久| 黄色一级大片看看| 欧美日韩av久久| 又大又黄又爽视频免费| 赤兔流量卡办理| 王馨瑶露胸无遮挡在线观看| 汤姆久久久久久久影院中文字幕| 亚洲精品一区蜜桃| 精品久久国产蜜桃| 日韩熟女老妇一区二区性免费视频| 成年人午夜在线观看视频| 少妇被粗大猛烈的视频| 免费人成在线观看视频色| 免费观看无遮挡的男女| 精品福利永久在线观看| 韩国精品一区二区三区 | 黄色配什么色好看| 下体分泌物呈黄色| 亚洲欧洲国产日韩| 蜜桃在线观看..| 美国免费a级毛片| 91aial.com中文字幕在线观看| 丝袜喷水一区| 日韩成人av中文字幕在线观看| 日本黄大片高清| 中文字幕制服av| 欧美国产精品一级二级三级| 久久久久久久久久人人人人人人| 99热这里只有是精品在线观看| 国产精品国产三级国产av玫瑰| 亚洲,欧美精品.| av不卡在线播放| av福利片在线| 少妇被粗大的猛进出69影院 | 在线精品无人区一区二区三| 五月天丁香电影| 26uuu在线亚洲综合色| 成人午夜精彩视频在线观看| 热re99久久国产66热| 国产极品天堂在线| 午夜福利视频在线观看免费| 日本爱情动作片www.在线观看| 交换朋友夫妻互换小说| 乱人伦中国视频| 国产精品一区二区在线不卡| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久久久免| 少妇的逼好多水| videosex国产| 99热这里只有是精品在线观看| 国产一区亚洲一区在线观看| 少妇高潮的动态图| 侵犯人妻中文字幕一二三四区| 最近最新中文字幕免费大全7| 国产探花极品一区二区| 水蜜桃什么品种好| 桃花免费在线播放| 免费看光身美女| 水蜜桃什么品种好| 99久久中文字幕三级久久日本| 久久久久精品人妻al黑| 国产黄色视频一区二区在线观看| 精品一区二区三区四区五区乱码 | 国产精品人妻久久久影院| 成年人免费黄色播放视频| 51国产日韩欧美| 18禁裸乳无遮挡动漫免费视频| 精品人妻熟女毛片av久久网站| 国产 精品1| 日本vs欧美在线观看视频| 91精品三级在线观看| 国产精品99久久99久久久不卡 | av福利片在线| 国产精品人妻久久久久久| 2021少妇久久久久久久久久久| 久久久久精品久久久久真实原创| 欧美老熟妇乱子伦牲交| 亚洲精品日本国产第一区| 国产精品秋霞免费鲁丝片| 亚洲,欧美,日韩| 熟女人妻精品中文字幕| 国产精品国产三级国产av玫瑰| 国产成人精品福利久久| 黄网站色视频无遮挡免费观看|