• <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ù)量與群落特征
    天堂俺去俺来也www色官网| 亚洲国产日韩一区二区| freevideosex欧美| av一本久久久久| 99久久综合免费| 另类精品久久| 蜜桃在线观看..| 国产免费一级a男人的天堂| 最黄视频免费看| 欧美性感艳星| 国产乱来视频区| 国产av一区二区精品久久| 日本-黄色视频高清免费观看| 亚洲怡红院男人天堂| 久久青草综合色| 少妇熟女欧美另类| 亚洲性久久影院| 亚洲精品视频女| 爱豆传媒免费全集在线观看| 美女内射精品一级片tv| 欧美成人午夜免费资源| 大香蕉久久网| 亚洲国产av新网站| 久久午夜福利片| 国产欧美日韩一区二区三区在线 | 一本—道久久a久久精品蜜桃钙片| 在线观看免费高清a一片| 国产无遮挡羞羞视频在线观看| 亚洲精华国产精华液的使用体验| 国产女主播在线喷水免费视频网站| 十分钟在线观看高清视频www| 午夜福利视频精品| 日日啪夜夜爽| 高清午夜精品一区二区三区| 啦啦啦在线观看免费高清www| 亚洲精品国产av成人精品| 日日啪夜夜爽| 亚洲第一区二区三区不卡| 日本vs欧美在线观看视频| 久久 成人 亚洲| 18禁动态无遮挡网站| 精品人妻熟女毛片av久久网站| 大香蕉久久成人网| 国产av国产精品国产| 久热久热在线精品观看| 亚洲av二区三区四区| 日本av免费视频播放| 最近最新中文字幕免费大全7| 免费久久久久久久精品成人欧美视频 | 男人爽女人下面视频在线观看| 午夜福利影视在线免费观看| 婷婷色综合www| 最近中文字幕高清免费大全6| 校园人妻丝袜中文字幕| 国产在视频线精品| 亚洲精品第二区| 自线自在国产av| av专区在线播放| 成年人免费黄色播放视频| 免费观看av网站的网址| 亚洲精品久久久久久婷婷小说| 亚洲av成人精品一二三区| 2018国产大陆天天弄谢| 日日摸夜夜添夜夜爱| 简卡轻食公司| 亚洲国产精品999| 在线观看免费日韩欧美大片 | 久久鲁丝午夜福利片| 五月伊人婷婷丁香| 王馨瑶露胸无遮挡在线观看| 欧美97在线视频| 男的添女的下面高潮视频| 一级片'在线观看视频| 日韩亚洲欧美综合| 免费看不卡的av| 日产精品乱码卡一卡2卡三| 久久久久人妻精品一区果冻| 三级国产精品片| 久久亚洲国产成人精品v| 2022亚洲国产成人精品| 日日摸夜夜添夜夜爱| 亚洲图色成人| 日本黄大片高清| 2018国产大陆天天弄谢| 超碰97精品在线观看| 国产成人freesex在线| av免费观看日本| 高清不卡的av网站| 国产免费又黄又爽又色| 久久久午夜欧美精品| 欧美性感艳星| 十八禁高潮呻吟视频| 国产av码专区亚洲av| 黑人高潮一二区| 国产精品熟女久久久久浪| 男人添女人高潮全过程视频| 一区二区三区精品91| 国产成人av激情在线播放 | a级毛片在线看网站| 国产高清国产精品国产三级| 亚洲美女搞黄在线观看| 满18在线观看网站| 精品一区在线观看国产| 亚洲av不卡在线观看| 日韩熟女老妇一区二区性免费视频| 亚洲国产av新网站| 九色亚洲精品在线播放| 国产在线视频一区二区| 精品久久久久久电影网| 婷婷成人精品国产| 伦理电影大哥的女人| 日韩制服骚丝袜av| 丝袜脚勾引网站| 午夜免费男女啪啪视频观看| 又大又黄又爽视频免费| 好男人视频免费观看在线| 亚洲av成人精品一区久久| 交换朋友夫妻互换小说| 色视频在线一区二区三区| 18禁动态无遮挡网站| 久久国产精品大桥未久av| 国产精品人妻久久久影院| 午夜激情av网站| 久久av网站| 精品国产乱码久久久久久小说| 精品国产国语对白av| tube8黄色片| 亚洲欧美一区二区三区黑人 | 成人毛片a级毛片在线播放| 啦啦啦在线观看免费高清www| 亚洲婷婷狠狠爱综合网| 丝瓜视频免费看黄片| 高清午夜精品一区二区三区| 亚洲情色 制服丝袜| 亚洲国产欧美日韩在线播放| 男女免费视频国产| 国产国拍精品亚洲av在线观看| 久久人人爽人人片av| 欧美人与善性xxx| 成人国产麻豆网| 高清欧美精品videossex| 日韩一区二区视频免费看| 超色免费av| av国产精品久久久久影院| 久久精品国产a三级三级三级| 狠狠婷婷综合久久久久久88av| 有码 亚洲区| 日韩一区二区三区影片| 一级毛片黄色毛片免费观看视频| 美女xxoo啪啪120秒动态图| 18禁观看日本| 久久免费观看电影| 天美传媒精品一区二区| 久久久久久久亚洲中文字幕| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 精品熟女少妇av免费看| 久久ye,这里只有精品| 中文字幕亚洲精品专区| a级毛片在线看网站| 久久久久久久精品精品| 国产色婷婷99| 3wmmmm亚洲av在线观看| 国产精品嫩草影院av在线观看| 蜜桃在线观看..| 国产国语露脸激情在线看| 大香蕉97超碰在线| 亚洲性久久影院| 精品一区二区免费观看| 久久精品国产亚洲网站| 天堂中文最新版在线下载| 亚洲国产精品一区三区| 丝瓜视频免费看黄片| 交换朋友夫妻互换小说| 精品国产国语对白av| 久热这里只有精品99| 美女主播在线视频| 亚洲综合精品二区| 麻豆成人av视频| 国产成人91sexporn| 黄色视频在线播放观看不卡| 免费高清在线观看日韩| 黑人欧美特级aaaaaa片| 美女中出高潮动态图| 欧美日韩视频精品一区| 国产 精品1| 国产一区二区三区综合在线观看 | 亚洲欧美清纯卡通| 国产精品人妻久久久久久| 免费播放大片免费观看视频在线观看| 99热6这里只有精品| 亚洲av二区三区四区| √禁漫天堂资源中文www| 美女视频免费永久观看网站| 22中文网久久字幕| 黄色毛片三级朝国网站| av视频免费观看在线观看| 精品久久久久久久久亚洲| 五月天丁香电影| 99国产精品免费福利视频| 日韩大片免费观看网站| 久久久久久久国产电影| 久久99精品国语久久久| 在现免费观看毛片| 国产高清不卡午夜福利| 18禁在线播放成人免费| 男人添女人高潮全过程视频| 草草在线视频免费看| 国产 精品1| 性高湖久久久久久久久免费观看| 97在线视频观看| 狠狠精品人妻久久久久久综合| 午夜福利网站1000一区二区三区| 久久97久久精品| 日韩精品免费视频一区二区三区 | 国产精品人妻久久久影院| 婷婷色av中文字幕| 永久网站在线| 成人无遮挡网站| 热99久久久久精品小说推荐| av又黄又爽大尺度在线免费看| 蜜桃国产av成人99| 这个男人来自地球电影免费观看 | av网站免费在线观看视频| 亚洲色图综合在线观看| 18+在线观看网站| 美女国产视频在线观看| 亚洲国产成人一精品久久久| 婷婷色麻豆天堂久久| 王馨瑶露胸无遮挡在线观看| 99国产精品免费福利视频| 五月玫瑰六月丁香| 亚洲精品一区蜜桃| 免费观看a级毛片全部| 国内精品宾馆在线| 久久女婷五月综合色啪小说| 在线观看国产h片| 国产亚洲av片在线观看秒播厂| 国产黄频视频在线观看| 国产成人aa在线观看| 亚洲成人av在线免费| 国产日韩一区二区三区精品不卡 | 99热国产这里只有精品6| a 毛片基地| 久久女婷五月综合色啪小说| 亚洲欧美日韩卡通动漫| 国产毛片在线视频| 国国产精品蜜臀av免费| 欧美97在线视频| 日韩不卡一区二区三区视频在线| 亚洲欧美日韩另类电影网站| 国产熟女欧美一区二区| 亚洲高清免费不卡视频| a级片在线免费高清观看视频| 十分钟在线观看高清视频www| 国产一级毛片在线| av网站免费在线观看视频| 一区二区三区免费毛片| 女人久久www免费人成看片| 一区二区三区四区激情视频| 久久精品熟女亚洲av麻豆精品| 日本wwww免费看| 国产视频内射| 91午夜精品亚洲一区二区三区| 亚洲综合精品二区| 久久久久久久久久人人人人人人| 亚洲欧美日韩卡通动漫| 日韩成人伦理影院| 免费av中文字幕在线| av在线app专区| 日本av免费视频播放| 精品人妻熟女av久视频| 三级国产精品欧美在线观看| 中文字幕精品免费在线观看视频 | 视频区图区小说| 国产日韩欧美视频二区| 91午夜精品亚洲一区二区三区| 一本大道久久a久久精品| 国产精品一区二区三区四区免费观看| 国产精品女同一区二区软件| 日韩欧美一区视频在线观看| 大片电影免费在线观看免费| 久久国产精品男人的天堂亚洲 | 中国国产av一级| 精品一区二区三卡| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品456在线播放app| 久久久久久久久久久久大奶| 精品一区二区三区视频在线| 肉色欧美久久久久久久蜜桃| 午夜老司机福利剧场| 不卡视频在线观看欧美| 亚洲情色 制服丝袜| 人妻 亚洲 视频| 免费不卡的大黄色大毛片视频在线观看| 在线亚洲精品国产二区图片欧美 | 97在线人人人人妻| 菩萨蛮人人尽说江南好唐韦庄| 免费看av在线观看网站| 大香蕉久久成人网| 国产精品一二三区在线看| 两个人的视频大全免费| 美女国产高潮福利片在线看| 国产女主播在线喷水免费视频网站| 插阴视频在线观看视频| 国产精品.久久久| 综合色丁香网| 18+在线观看网站| 少妇的逼好多水| 国产精品麻豆人妻色哟哟久久| 亚洲欧美成人综合另类久久久| 欧美亚洲日本最大视频资源| 亚洲天堂av无毛| 少妇熟女欧美另类| 免费高清在线观看视频在线观看| 国产日韩欧美视频二区| av在线观看视频网站免费| 国产成人aa在线观看| 在线看a的网站| av卡一久久| 亚洲国产日韩一区二区| 亚洲色图综合在线观看| 不卡视频在线观看欧美| 亚洲,一卡二卡三卡| 人妻一区二区av| 欧美少妇被猛烈插入视频| 成人综合一区亚洲| 九色成人免费人妻av| 久久久久久久亚洲中文字幕| 在线观看三级黄色| 日韩视频在线欧美| 精品久久蜜臀av无| 日本欧美视频一区| 一边亲一边摸免费视频| 一区二区日韩欧美中文字幕 | 国产片内射在线| 久久精品国产鲁丝片午夜精品| tube8黄色片| 丝袜喷水一区| 你懂的网址亚洲精品在线观看| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线 | 国产成人av激情在线播放 | 久久久久久久久久久免费av| 亚洲av.av天堂| 亚洲精品日本国产第一区| 赤兔流量卡办理| 日本色播在线视频| 午夜老司机福利剧场| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 我要看黄色一级片免费的| 国产精品国产av在线观看| 亚洲av成人精品一二三区| 最近2019中文字幕mv第一页| 亚洲中文av在线| 国产亚洲欧美精品永久| 狠狠精品人妻久久久久久综合| 成人毛片60女人毛片免费| 最近中文字幕高清免费大全6| 国产日韩欧美在线精品| 国产一区二区三区av在线| 色婷婷久久久亚洲欧美| 大香蕉久久网| 亚洲高清免费不卡视频| 精品国产露脸久久av麻豆| 最后的刺客免费高清国语| 99精国产麻豆久久婷婷| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 日韩亚洲欧美综合| 久久免费观看电影| 国产精品一二三区在线看| av在线观看视频网站免费| 在线 av 中文字幕| 大香蕉久久网| 国产成人精品在线电影| 高清黄色对白视频在线免费看| 午夜免费鲁丝| 中文字幕免费在线视频6| 亚洲成色77777| 国产亚洲最大av| a 毛片基地| 久久国产精品大桥未久av| 免费av中文字幕在线| 久久久久久久久久久丰满| 最黄视频免费看| 亚洲精品成人av观看孕妇| 免费观看在线日韩| 精品一区在线观看国产| 极品少妇高潮喷水抽搐| 天天操日日干夜夜撸| 精品国产国语对白av| 久久毛片免费看一区二区三区| 又大又黄又爽视频免费| 国产视频内射| 视频在线观看一区二区三区| 精品一区二区三区视频在线| 国产成人精品在线电影| av播播在线观看一区| 午夜91福利影院| 久久毛片免费看一区二区三区| 三级国产精品片| 亚洲欧美一区二区三区黑人 | 日本黄大片高清| 久久 成人 亚洲| 少妇的逼水好多| 有码 亚洲区| 伊人久久国产一区二区| 国产av码专区亚洲av| 欧美性感艳星| 欧美97在线视频| 美女国产视频在线观看| 99久久人妻综合| 欧美bdsm另类| 亚洲综合精品二区| 国产色婷婷99| 成人国语在线视频| 精品人妻一区二区三区麻豆| 五月伊人婷婷丁香| 午夜激情av网站| 永久网站在线| 免费观看的影片在线观看| 久久久国产精品麻豆| 亚洲一级一片aⅴ在线观看| 99久久人妻综合| 哪个播放器可以免费观看大片| 精品一区二区三区视频在线| 国产白丝娇喘喷水9色精品| xxx大片免费视频| 亚洲精品456在线播放app| 男女免费视频国产| 91精品伊人久久大香线蕉| 久久久久久久久久成人| 国产精品欧美亚洲77777| 亚洲av在线观看美女高潮| 中文欧美无线码| 97精品久久久久久久久久精品| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 国产免费又黄又爽又色| 日韩av在线免费看完整版不卡| 精品人妻熟女av久视频| 国产熟女欧美一区二区| 校园人妻丝袜中文字幕| 国产精品国产av在线观看| 美女国产视频在线观看| 久久99热6这里只有精品| 国产精品久久久久久久电影| 久久 成人 亚洲| 国产精品99久久久久久久久| 中文字幕最新亚洲高清| 精品一品国产午夜福利视频| 97超视频在线观看视频| 亚洲不卡免费看| av免费观看日本| 建设人人有责人人尽责人人享有的| 亚洲精品aⅴ在线观看| 婷婷色综合www| 国产成人精品在线电影| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 欧美97在线视频| 18禁裸乳无遮挡动漫免费视频| 乱人伦中国视频| 又大又黄又爽视频免费| 欧美日韩亚洲高清精品| 一区在线观看完整版| 99九九在线精品视频| 久久精品人人爽人人爽视色| 成年人午夜在线观看视频| 这个男人来自地球电影免费观看 | 精品一区二区免费观看| 一区二区日韩欧美中文字幕 | 夜夜骑夜夜射夜夜干| 久久久亚洲精品成人影院| 美女视频免费永久观看网站| 亚洲一级一片aⅴ在线观看| av视频免费观看在线观看| 啦啦啦中文免费视频观看日本| 久久久久久久久久久免费av| 久久久久精品性色| 免费观看a级毛片全部| 在线观看人妻少妇| 国产视频首页在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲精品日本国产第一区| 日本与韩国留学比较| 老熟女久久久| 黑人猛操日本美女一级片| 人妻制服诱惑在线中文字幕| 国产精品熟女久久久久浪| 在线观看免费日韩欧美大片 | 欧美国产精品一级二级三级| 国产成人精品福利久久| 国产片特级美女逼逼视频| 日本与韩国留学比较| 国产乱人偷精品视频| √禁漫天堂资源中文www| 久久鲁丝午夜福利片| 国产亚洲精品久久久com| 这个男人来自地球电影免费观看 | 日韩av免费高清视频| 在线精品无人区一区二区三| 精品少妇内射三级| 母亲3免费完整高清在线观看 | av在线播放精品| 日本-黄色视频高清免费观看| 麻豆精品久久久久久蜜桃| 亚洲国产av新网站| 考比视频在线观看| 精品少妇内射三级| 日韩av在线免费看完整版不卡| 纯流量卡能插随身wifi吗| 亚洲国产欧美在线一区| 一个人看视频在线观看www免费| av电影中文网址| 日产精品乱码卡一卡2卡三| 亚洲国产精品国产精品| 激情五月婷婷亚洲| 免费高清在线观看日韩| 秋霞在线观看毛片| 日韩制服骚丝袜av| 国产精品人妻久久久久久| 精品卡一卡二卡四卡免费| 男的添女的下面高潮视频| 中文字幕亚洲精品专区| 精品久久久久久久久av| 性高湖久久久久久久久免费观看| 国产精品99久久99久久久不卡 | 精品一区二区三卡| 国产高清有码在线观看视频| 亚洲精品aⅴ在线观看| 久久久久久久大尺度免费视频| xxx大片免费视频| 十八禁网站网址无遮挡| 中文乱码字字幕精品一区二区三区| 欧美人与性动交α欧美精品济南到 | 自拍欧美九色日韩亚洲蝌蚪91| 一个人看视频在线观看www免费| 99久国产av精品国产电影| av播播在线观看一区| 国产精品麻豆人妻色哟哟久久| 欧美日韩视频高清一区二区三区二| 日本-黄色视频高清免费观看| 蜜桃国产av成人99| 午夜激情久久久久久久| 五月开心婷婷网| 国产精品久久久久久精品电影小说| 成人黄色视频免费在线看| 成人综合一区亚洲| 91久久精品国产一区二区成人| 精品人妻偷拍中文字幕| 久久精品国产亚洲av天美| 国产视频首页在线观看| 熟女电影av网| 999精品在线视频| 大片电影免费在线观看免费| 女性生殖器流出的白浆| 99热6这里只有精品| 桃花免费在线播放| 亚洲国产精品999| 国内精品宾馆在线| 久久精品国产自在天天线| h视频一区二区三区| 毛片一级片免费看久久久久| 婷婷色麻豆天堂久久| 日本91视频免费播放| 欧美日韩视频高清一区二区三区二| 五月开心婷婷网| 99久久精品国产国产毛片| 99精国产麻豆久久婷婷| 最新中文字幕久久久久| 性色av一级| 欧美人与性动交α欧美精品济南到 | 中文字幕免费在线视频6| 亚洲国产av影院在线观看| 丝袜喷水一区| 色哟哟·www| 亚洲国产精品999| av卡一久久| 热re99久久国产66热| 成人亚洲精品一区在线观看| 人妻人人澡人人爽人人| 永久免费av网站大全| 国产视频内射| 国产国语露脸激情在线看| 国产高清三级在线| 欧美国产精品一级二级三级| 久久久久久久久久成人| 大片免费播放器 马上看| 寂寞人妻少妇视频99o| 午夜激情福利司机影院| 五月伊人婷婷丁香| 精品亚洲成国产av| 色哟哟·www| 亚洲在久久综合| 一级毛片我不卡| 成人毛片a级毛片在线播放| 在线观看www视频免费| 精品久久久久久久久亚洲| 日本91视频免费播放| 精品一区在线观看国产| 黄色配什么色好看| 亚洲av免费高清在线观看| 国产极品天堂在线| 亚洲av不卡在线观看| 久久久久视频综合| 欧美激情 高清一区二区三区| 亚洲成人手机| 永久网站在线| 国产亚洲一区二区精品| 我的女老师完整版在线观看|