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

    新疆達坂城鹽湖不同植物群落土壤呼吸研究①

    2017-07-21 05:05:48李典鵬姚美思韓東亮王麗萍胡保安王寧寧蔣大勇朱新萍賈宏濤
    土壤 2017年3期
    關鍵詞:芨芨草鳶尾土壤濕度

    李典鵬,姚美思,韓東亮,王麗萍,胡保安,王寧寧,蔣大勇,張 凱,朱新萍,賈宏濤

    ?

    新疆達坂城鹽湖不同植物群落土壤呼吸研究①

    李典鵬,姚美思,韓東亮,王麗萍,胡保安,王寧寧,蔣大勇,張 凱,朱新萍,賈宏濤*

    (新疆農(nóng)業(yè)大學草業(yè)與環(huán)境科學學院,烏魯木齊 830052)

    新疆鹽湖分布地域廣泛而不均勻。鹽湖生態(tài)系統(tǒng)碳循環(huán)是干旱區(qū)生態(tài)系統(tǒng)碳循環(huán)的重要組分,其動態(tài)變化對干旱區(qū)碳平衡有著重要的影響,然而目前對干旱區(qū)鹽湖土壤呼吸的研究相對薄弱。為探討鹽湖沿岸不同植物群落土壤呼吸及其影響因子,以新疆達坂城鹽湖為典型研究區(qū),選取沿岸小獐毛()、鳶尾()、芨芨草()、黑果枸杞()植物群落土壤為主要研究對象,以周邊農(nóng)田撂荒地為對照,利用LI-COR 8100對達坂城鹽湖沿岸4種植物群落和撂荒地土壤呼吸進行了監(jiān)測,結果表明:鹽湖沿岸土壤呼吸具有明顯的日變化。監(jiān)測時段內(nèi)土壤呼吸作用均為單峰曲線,排放通量在13:00左右最高,07:00左右出現(xiàn)最低值。不同群落土壤呼吸速率變化范圍在0.89 ~ 4.34 μmol/(m2·s)之間,撂荒地為1.77 ~ 3.48 μmol/(m2·s),表現(xiàn)為黑果枸杞>小獐毛>芨芨草>撂荒地>鳶尾。不同植物群落地下生物量影響土壤CO2排放,研究區(qū)土壤呼吸日平均累計CO2排放通量為209.28 mg/(m2·d)。土壤呼吸速率與土壤5 cm溫度呈現(xiàn)顯著正相關(<0.01),可以解釋土壤呼吸日變化的69% ~ 78%,是決定達坂城鹽湖土壤呼吸變化的主要因子。利用土壤5 cm處溫度得到小獐毛、鳶尾、芨芨草、黑果枸杞、撂荒地的10值分別為2.30、1.24、1.20、1.57、1.28,小獐毛樣地土壤呼吸對溫度更敏感。除撂荒地外各群落土壤呼吸與土壤濕度呈顯著負相關(2= 0.38 ~ 0.51)。該研究對進一步明確干旱區(qū)鹽湖生態(tài)系統(tǒng)土壤呼吸變化規(guī)律和影響因素,以及對估算區(qū)域碳平衡及“碳匯”功能具有重要意義。

    干旱區(qū);鹽湖;碳循環(huán);土壤呼吸;影響因子

    呼吸作用是生物能學的一個重要組成部分,由于土壤生物過程十分復雜,故現(xiàn)階段的土壤呼吸主要集中在研究土壤CO2釋放量[1]。土壤呼吸作為全球生態(tài)系統(tǒng)碳循環(huán)的重要環(huán)節(jié),也是全球氣候變化的關鍵生態(tài)過程。不同時間尺度上探討影響土壤呼吸特征及其敏感環(huán)境因素,對準確評估生態(tài)系統(tǒng)碳收支有著重要的意義。作為溫室氣體的主要組成部分,CO2在限制地表紅外輻射向外傳送起到重要的作用,更多的熱量被存儲在了大氣中,從而導致氣候變暖。特定條件下土壤作為“碳匯”,對緩解區(qū)域氣候變化起到了一定的作用;據(jù)估計,在全球尺度上,土壤每年向大氣中排放77 Pg C[2],土壤呼吸作為“碳源”釋放碳的速率比人類活動釋放碳的速率要大一個數(shù)量級。因此,在氣候變暖與全球碳循環(huán)之間的正反饋中,土壤呼吸是一個關鍵的過程。為了正確評價氣候變化中的不確定性,需要了解不同類型土壤及特定條件下的呼吸特征及其影響因素。

    鹽湖是濕地生態(tài)系統(tǒng)的重要組成部分,具有鹽分和水分雙重脅迫的特點,干旱區(qū)鹽湖濕地土壤呼吸的研究相對薄弱。在極端干旱的鹽生荒漠系統(tǒng),鹽堿土對大氣中的CO2有很好的吸收能力,吸收CO2量可達62 ~ 622 g/(m2·a)[3-4],其吸收量相當于植物的2 ~ 15倍,且濕潤年份日累積最大CO2吸收量是干旱年份的4倍[5]。在一定程度上說明水分能夠影響到鹽堿土的“碳匯”功能,鹽湖作為干旱區(qū)特殊的濕地,無論是從潮濕土壤有機固碳還是鹽漬化土壤無機固碳的角度來看,均應該是典型“碳匯”。但其中的核心問題取決于鹽湖周邊土壤的呼吸,即土壤CO2排放特征。

    通過對鹽湖周邊土壤呼吸的研究,明確驅動鹽湖土壤呼吸的關鍵因素,可為濕地生態(tài)系統(tǒng)在全球氣候變化中作用的準確評估提供參考數(shù)據(jù)。本研究擬以新疆達坂城鹽湖為典型研究區(qū)域,以其周邊生長不同主要植物群落立地土壤為研究對象,對土壤呼吸速率及其影響因子進行原位監(jiān)測,明確達坂城鹽湖沿岸不同植物群落土壤呼吸速率日變化特征,為進一步明確內(nèi)陸鹽湖濕地在干旱區(qū)陸地碳循環(huán)中的地位提供科學依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    達坂城鹽湖又名鹽湖、破城子湖、烏魯木齊湖等[6],位于新疆維吾爾自治區(qū)烏魯木齊市區(qū)東南45 km達坂城區(qū)西部(88°03′53″ ~ 88°12′15″E,43°21′00″ ~43°25′25″N)(圖1)。湖盆為封閉內(nèi)流盆地,周圍無常年性地表河流,靠大氣降水和地下水補給,面積35 km2,湖面海拔1 066 m。湖盆在柴窩堡山間盆地東南部,南北高山夾持,氣候干旱,年平均氣溫 6.1℃(月平均氣溫–22.1℃,7月平均氣溫 35.8℃,年較差 31.1℃),日較差12.1℃;年降水量71.8 mm,年蒸發(fā)量2 648 mm,相對濕度38%;冰凍期5個月,年平均降(積)雪150 d,最大積雪深度30 ~ 40 cm;年日照數(shù)2 500 ~ 2 600 h,全年多刮西北風,最大風力10 級,最大風速30 m/s,屬于干旱大陸性氣候。鹽湖沿岸主要植被類型為小獐毛()、鳶尾()、芨芨草()、黑果枸杞()、苔草()等。

    1.2 研究方案

    試驗于2015年7—9月進行,以鹽湖為核心,根據(jù)試驗區(qū)植物群落特點,由湖心向外依次選取小獐毛()、鳶尾()、芨芨草()、黑果枸杞()群落和撂荒地(Abandoned land)土壤作為監(jiān)測對象,研究樣地基本性質見表1。在每個群落內(nèi)隨機選取3塊樣地,每塊樣地10 m × 10 m,沿樣方對角線設置2個長期固定的底座,每月中下旬選擇2 ~ 3個相對晴朗的天氣進行土壤呼吸的監(jiān)測,從上午07:00至晚上19:00,每2小時監(jiān)測1次,監(jiān)測土壤呼吸的同時監(jiān)測環(huán)境因素指標,包括近地面氣溫,土壤5、10、15、20和25 cm處土溫,表土層含水量(0 ~ 10 cm)、氣壓、風速、光照強度等。

    1.3 研究方法

    1.3.1 土壤呼吸測定 利用土壤碳通量測量系統(tǒng)LI-COR 8100(Li-COR,Lincoln,NE,USA)測定土壤呼吸速率。每塊樣地內(nèi)設兩個內(nèi)徑為20 cm、高10 cm 的PVC管基座,將基座嵌入土中,使上端高出地表約5 cm,剪除基座環(huán)內(nèi)的植物。為減少由于嵌入基座對土壤擾動而造成短期內(nèi)土壤呼吸速率的波動,經(jīng)過24 h 的平衡后,土壤呼吸速率恢復到正常水平才進行相關指標的測定。監(jiān)測時每個基座測定1 次,即每個處理6個監(jiān)測數(shù)據(jù)。

    表1 達坂城鹽湖研究樣地概況

    1.3.2 土壤溫度、濕度測定 與土壤呼吸監(jiān)測同步,監(jiān)測土壤濕度和土壤溫度;土壤濕度采用野外采樣-室內(nèi)烘干法進行,以基座為圓心在半徑為1 m的圓形范圍內(nèi)均為土壤濕度的采樣范圍;5 cm土壤溫度采用LI-COR 8100 自帶的溫度探頭,監(jiān)測前一天在PVC管基座附近安置地溫計,同步對10、15、20、25 cm處土壤地溫進行監(jiān)測。

    1.3.3 植被調(diào)查方法 植被調(diào)查采用常規(guī)方法,由于不同植物群落的差異,所以調(diào)查方法也略有不同。小獐毛群落和撂荒地植被調(diào)查方法如下:在所選樣地附近,按一定方向設置50 m 樣線,每隔10 m 布設一個1 m×1 m 的樣方。鳶尾、芨芨草、黑果枸杞群落生物量測定,按一定方向設置50 m 樣線,則為每隔10 m 布設一個5 m×5 m 的樣方,對樣方內(nèi)地表活體植物和凋落物進行采集和分析,估計樣方內(nèi)植被密度、蓋度,用稱重法測定生物量。根系樣品采集采用根鉆法,選擇剪完草的5個樣方,采取0 ~ 5、5 ~ 10、10 ~ 20、20 ~ 30、30 ~ 50、50 ~ 70、70 ~ 100 cm 分層根系樣品。

    1.3.4 數(shù)據(jù)處理 采用單因素方差分析(ANOVA)檢驗不同群落土壤呼吸速率的顯著性。對所得數(shù)據(jù)采用 Excel 2016、SPSS 20.0 進行數(shù)據(jù)處理分析和制圖,繪圖工具為MapInfo 8.5和Adobe Photoshop CS5。采用 Van’t Hoff 模型分析土壤溫度與土壤呼吸速率的關系:

    = aebT(1)

    式中:表示土壤呼吸(μmol/(m2·s),下同),表示土壤溫度(℃,下同),a、b均為參數(shù)(下同)。

    采用線性回歸分析的方法分析土壤濕度與土壤呼吸速率的關系:

    = a+ c (2)

    式中:表示土壤濕度(%),c為參數(shù)(下同)。

    采用多性回歸分析法分析土壤水分和溫度與土壤呼吸速率的關系:

    = a+ b+ c (3)

    采用公式(4)計算土壤呼吸對溫度變化敏感程度:

    10= e10b(4)

    式中:10為土壤呼吸的溫度敏感性,即土壤溫度升高10℃時土壤呼吸速率增加或減少的倍數(shù);b為Van’t Hoff 指數(shù)模型 (= aebT) 中的溫度反應常數(shù)。

    2 結果與分析

    2.1 不同植被群落土壤呼吸日動態(tài)

    圖2為達坂城鹽湖不同植被群落土壤呼吸速率的日動態(tài),可以看出,各土壤呼吸速率的日變化均為先增后降的單峰曲線。不同植被類型土壤呼吸的日變化特征存在差異,大小為黑果枸杞>小獐毛>芨芨草>鳶尾,其中鳶尾群落土壤呼吸速率明顯低于其他群落。曲線的最高峰值均出現(xiàn)在13:00左右,小獐毛、鳶尾、芨芨草、黑果枸杞和撂荒地曲線峰值分別為4.02、2.26、3.37、4.94、3.48 μmol/(m2·s)。除撂荒地外,各個群落土壤呼吸速率最低值出現(xiàn)在19:00;撂荒地土壤呼吸速率在17:00時降低到2.05 μmol/(m2·s),隨后其呼吸速率又呈現(xiàn)增加趨勢。

    2.2 植被生物量對土壤CO2排放通量的影響

    從圖3可以看出,小獐毛、鳶尾、芨芨草、黑果枸杞群落和撂荒地地上活體生物量分別為144.56、32.70、55.21、51.47、21.9 g/m2,小獐毛群落土壤地上活體生物量與其他4個處理差異達到極顯著水平(<0.01)。小獐毛、鳶尾、芨芨草、黑果枸杞和撂荒地地表凋落物生物量分別為84.81、14.84、33.36、32.13、9.43 g/m2,小獐毛群落土壤地表凋落物與其他群落的凋落物差異達顯著水平(<0.05);鳶尾與芨芨草、芨芨草與黑果枸杞地表凋落物差異性未達到顯著水平(>0.05)。小獐毛、鳶尾、芨芨草、黑果枸杞和撂荒地土壤CO2排放總量分別為5.35、3.55、5.16、6.23、4.84 g/(m2·d);各群落和撂荒地之間土壤CO2日平均通量的差異性達到極顯著水平(<0.01)??傮w來看,活體生物量大、凋落物高的小獐毛群落土壤,其CO2排放通量大于生物量小、凋落物低的鳶尾、芨芨草群落和撂荒地,低于黑果枸杞群落。

    本文利用土壤碳通量測量系統(tǒng)LI-COR 8100測定土壤呼吸,而土壤呼吸包括了微生物呼吸和根呼吸以及無機過程,微生物呼吸與土壤中微生物數(shù)量密切相關;根呼吸與植物地下根系生物量有關;在鹽生荒漠系統(tǒng)中土壤溫度的變化改變了CO2在土壤溶液中的溶解度,催生出一個日出夜進的“無機呼吸”過程,此過程與環(huán)境溫度的變化關系密切。通過分析不同植物群落土壤0 ~ 100 cm 根系與土壤CO2日平均通量的關系,結果表明不同植物群落土壤根系生物量差異性達到極顯著水平,且各群落土壤CO2排放通量也存在顯著性差異,通過對根系生物量與CO2通量擬合的結果表明,土壤呼吸與根系生物量呈線性相關:= 0.111 7+ 3.385 8(2 = 0.64,<0.05),說明土壤呼吸也受地下根系生物量的影響。

    2.3 土壤呼吸與土壤溫度、濕度的關系

    在監(jiān)測時段內(nèi),各群落和撂荒地土壤5 cm溫度呈現(xiàn)出先增加后下降的趨勢(圖4),但不同群落土壤表層溫度峰值出現(xiàn)的時間存在差異,小獐毛、鳶尾、芨芨草群落土壤峰值出現(xiàn)在17:00,分別為32.3℃、39.2℃、42.3℃,撂荒地的峰值出現(xiàn)在15:00,為41.4℃,黑果枸杞群落土壤最高溫度出現(xiàn)在13:00,為40.8℃。除小獐毛群落10 cm土壤溫度日變化為單增曲線之外,其他監(jiān)測點10 cm土壤溫度變化均為單峰曲線;小獐毛群落10 cm土壤溫度最高值較表層峰值滯后1 h,為27.6℃。鳶尾、芨芨草群落和撂荒地10 cm土壤溫度峰值與土壤表層峰值出現(xiàn)時刻一致,黑果枸杞群落10 cm土壤溫度最高值較土壤表層溫度峰值滯后4 h出現(xiàn)。且各群落土壤10 cm以下土層溫度變幅不大。

    通過 Van’t Hoff 指數(shù)方程分別擬合土壤呼吸速率與土壤5、10、15、20、25 cm 溫度的回歸關系,并用決定系數(shù)(2) 來解釋土壤呼吸速率的變異量,結果表明土壤呼吸速率與5 cm地溫的擬合關系最為理想,即Van’t Hoff指數(shù)方程能很好地描述土壤呼吸對土壤5 cm 溫度的響應(表2)。不同群落及撂荒地土壤5 cm溫度與土壤呼吸速率具有指數(shù)關系(圖4),在> 99% 的置信區(qū)間內(nèi),土壤5 cm溫度的變化能夠解釋土壤呼吸變化的69% ~ 78%,而10 cm 土層土壤溫度的變化可以解釋土壤呼吸變化的24% ~ 44%。然而,在對15、20、25 cm土層溫度與土壤呼吸速率進行擬合時,三者決定系數(shù)與5、10 cm決定系數(shù)相比有明顯減弱,因此15 cm 以及15 cm 以下土層溫度對土壤呼吸的貢獻率可以不用考慮。5 cm 土層土壤溫度對土壤呼吸速率的決定系數(shù)高于10 cm土層,說明在鹽湖沿岸生態(tài)系統(tǒng)中,5 cm 土層土壤溫度可以作為影響土壤呼吸急劇變化的變量。

    表2 土壤呼吸的溫度敏感系數(shù)及土壤呼吸日變化與土壤溫度指數(shù)回歸關系

    土壤呼吸溫度敏感系數(shù)(10)大小順序依次為小獐毛群落>黑枸杞群落>撂荒地>鳶尾群落>芨芨草群落(表2),小獐毛群落土壤的10值為2.30,明顯高于其他4個處理,說明小獐毛群落土壤呼吸對溫度變化的響應程度大于其他群落土壤。但是土壤呼吸和土壤溫度、土壤濕度等生物環(huán)境因素有著多重關系,所以如果僅用10解釋鹽湖不同群落和撂荒地土壤呼吸的變化規(guī)律具有一定的不確定性。

    各群落土壤濕度與土壤呼吸的關系如圖6所示。不同群落土壤呼吸與土壤濕度進行擬合分析可知(表3),各群落土壤濕度與土壤呼吸速率呈顯著負相關,土壤呼吸曲線決定系數(shù)在0.38 ~ 0.51之間;撂荒地土壤濕度與土壤呼吸速率呈顯著正相關(2= 0.38)。各群落和撂荒地土壤濕度和土壤呼吸之間線性關系達到顯著水平,說明土壤濕度也是影響土壤呼吸的因子。4種植物群落土壤濕度與土壤呼吸相關性決定系數(shù)大小順序為鳶尾>小獐毛>芨芨草>黑果枸杞,決定系數(shù)的大小在一定程度上反映不同群落土壤呼吸速率對水分的依賴程度,出現(xiàn)此現(xiàn)象的原因可能與植物群落與湖邊的距離有關,距離不同,水分的脅迫程度不同。

    表3 土壤呼吸與土壤濕度的線性回歸關系

    2.4 土壤呼吸與土壤溫度及土壤濕度的線性回歸分析

    不同群落土壤呼吸與土壤表層溫度、濕度的線性回歸分析見表4。各群落和撂荒地土壤呼吸與溫度的相關系數(shù)分別為0.85、0.84、0.73、0.77、0.76,相關性顯著。在各群落和撂荒地土壤回歸方程中,土壤濕度和表層溫度對土壤呼吸的影響均存在較小的交互作用,其中芨芨草群落土壤濕度與表層溫度對土壤呼吸的交互作用大于小獐毛和鳶尾群落,在回歸方程中,土壤濕度和表層溫度均進入回歸方程(95% 置信區(qū)間),以溫度-水分-土壤呼吸雙因子模型可以解釋土壤呼吸的變化。從各群落土壤呼吸與濕度、溫度的相關系數(shù)可以看出,沿湖心向外土壤呼吸對溫度的響應呈先降低再增加的變化趨勢;不同群落土壤呼吸與土壤濕度呈顯著的負相關,由芨芨草群落到小獐毛群落土壤呼吸和含水量的負相關依次增強,即鳶尾>小獐毛>芨芨草>黑果枸杞。撂荒地土壤呼吸與土壤溫度、濕度呈正相關。

    表4 土壤呼吸與土壤溫度和土壤濕度的函數(shù)模擬

    注:R表示土壤呼吸速率,表示溫度,表示土壤濕度;** 表示達到<0.01顯著水平,* 表示達到<0.05顯著水平。

    3 討論

    3.1 不同植物群落土壤呼吸日變化特征

    鹽湖生態(tài)系統(tǒng)不同群落條件下土壤呼吸日變化呈現(xiàn)明顯單峰曲線,此結果與國內(nèi)外大部分學者研究相同[7-9],峰值均出現(xiàn)在13:00,這可能與土壤溫度和濕度有關,凌晨土壤溫度較低,且土壤濕度較高,二者同時限制土壤呼吸,導致土壤呼吸較弱。隨著氣溫的上升,土壤溫度與土壤濕度也隨之變化,同時也促使土壤呼吸速率的增加,當土壤溫度和土壤濕度在13:00達到適宜條件時,土壤酶活性較高,促進了植物根系和微生物的代謝速率[10-11],故土壤呼吸在此時出現(xiàn)監(jiān)測時段內(nèi)的峰值。13:00后土壤溫度持續(xù)上升,由于溫度升高和蒸發(fā)速率加強導致土壤濕度下降,濕度與溫度交互作用從而使土壤呼吸減弱,故13:00后土壤呼吸速率呈下降趨勢。在干旱區(qū)鹽生荒漠,馬杰等[12]對原始鹽生荒漠多枝檉柳群落土壤CO2通量觀測表明,鹽生荒漠土壤CO2排放通量呈現(xiàn)出單峰的日變化過程。相關研究結果表明,土壤呼吸日變化相對于土壤溫度的變化存在不定的提前效應[10]或滯后[13]現(xiàn)象。在本研究中,土壤呼吸速率最大值出現(xiàn)在13:00,相對于土壤表層溫度最高值出現(xiàn)提前2 h,這與Jassal等[14]的研究結果一致,其原因可能是土壤溫度和含水量的變化所導致,除此之外底物供應也可能是調(diào)節(jié)土壤呼吸日變化的一個重要因子[15]。

    3.2 植物類型對土壤呼吸的影響

    植物是生態(tài)系統(tǒng)的重要組成部分,通過植物的相關指標可以間接地反映生態(tài)系統(tǒng)的穩(wěn)定性和環(huán)境狀況。植物類型和生物量的大小是影響土壤和大氣之間進行CO2交換的因素,不同群落有不同植物類型和覆蓋度,凋落物、生物量的數(shù)量方面也有較大差異,這些都有可能對土壤的呼吸速率造成一定的影響。通過對鹽湖生態(tài)系統(tǒng)的研究發(fā)現(xiàn),建群種分別為小獐毛、鳶尾、芨芨草、黑果枸杞,單位面積的生物量沿湖心向外有逐步遞減的趨勢,小獐毛樣地地上生物量與鳶尾、芨芨草地上生物量存在顯著性差異,且CO2通量也存在顯著性差異。黑果枸杞地上生物量和凋落物均小于小獐毛,但其呼吸速率大于其他3個群落和撂荒地。地表凋落物分解速率也是影響土壤與大氣之間CO2交換速率的因素之一,本研究表明,除撂荒地外,其他群落都呈現(xiàn)土壤呼吸速率隨凋落物增加而增加的趨勢。董自紅等[16]、王憶慧等[17]通過對高山草甸草原植物群落的研究結果表明土壤呼吸與凋落物存在線性正相關,本研究結果與此相同。與同處于干旱區(qū)的艾比湖相比,達坂城鹽湖沿岸不同群落土壤呼吸速率要高于艾比湖流域[18],究其原因,是因為艾比湖流域植被類型主要以胡楊、梭梭等木本植物為主,植物覆蓋度低、植物殘體較少,導致土壤微生物活動較弱、有機質的礦化速率緩慢,從而降低了土壤CO2的排放。不同土地利用方式對土壤呼吸也有不同影響,達坂城鹽湖沿岸土壤呼吸速率與橡膠林土壤相近[19],低于種植玉米農(nóng)田土壤[20]。在不考慮環(huán)境條件的情況下,此類結果均有可能是由于不同植被類型所導致。

    3.3 環(huán)境因素對土壤呼吸的影響

    從本質上來說,土壤呼吸是一個生物化學過程,土壤溫度、濕度幾乎影響土壤呼吸的各個方面。在生物化學水平,溫度通過影響酶活性從而影響根和根際微生物的呼吸速率;在物理水平,溫度和水分相互作用影響土壤孔隙結構與外界的氣體交換。就大部分情況而言,土壤表層及亞表層溫度在時間尺度上變化較為明顯,通過分析可知土壤溫度函數(shù)都可以用于解釋土壤呼吸的動態(tài)變化。微生物呼吸作為土壤CO2排放的來源之一,其釋放速率受溫度的影響較大,并且與土壤溫度有著密切的聯(lián)系;溫度通過影響土壤中各類酶活性、微生物代謝速率、土壤動物呼吸速率和植物根系生長從而影響土壤CO2的排放[21-22]。土壤呼吸與土壤溫度的擬合方式有冪函數(shù)、線性、拋物線和指數(shù)函數(shù)等,不同研究結果擬合模型也有所不同,但指數(shù)函數(shù)擬合的方法較多[23],原因可能與溫度敏感系數(shù)的計算有關。本研究中土壤呼吸與土壤溫度呈顯著的指數(shù)相關(2分別為0.74、0.71、0.69、0.78、0.69;<0.01),且土壤呼吸與土壤5 cm溫度的指數(shù)擬合相關性大于10 cm,這與陳寶玉等[24]研究結果不一致,這可能是由于鹽湖沿岸土壤5 cm溫度在時間尺度上的變化更為明顯以及土壤有機質水平較高造成的。相關研究表明全球各生態(tài)系統(tǒng)的10的變化范圍在1.3 ~ 3.3之間波動,中值為2.40[2]。本研究中小獐毛、鳶尾、芨芨草、黑果枸杞和撂荒地樣地的10分別為2.30、1.24、1.20、1.57、1.28,中值為1.28。與前人研究成果[25]相比,鹽湖沿岸群落土壤呼吸對溫度的敏感程度更低,這可能是由環(huán)境因素導致的,鹽湖夏季的日較差較低,溫度小范圍的變化降低了土壤呼吸的溫度敏感程度,從而降低了10值。

    土壤濕度是影響土壤呼吸的另一個重要因子。土壤濕度對土壤呼吸的直接影響是通過影響根和微生物的生理過程,對土壤呼吸的間接影響是通過影響底物和氧氣的擴散。相關研究結果表明,土壤呼吸速率和土壤濕度之間存在正相關、負相關或是沒有相關性[26]。本研究中,達坂城鹽湖沿岸不同群落土壤呼吸與土壤濕度呈顯著負相關,這與高艷紅等[27]在騰格里沙漠得到的結論相同,這可能是因為土壤含水量超過了田間持水量的50%,限制了微生物活性,從而降低了微生物的異養(yǎng)呼吸。撂荒地土壤呼吸速率與土壤濕度呈正相關,能夠解釋土壤呼吸速率變化的38%,這與方晰等[28]對杉木人工林地土壤呼吸研究的結果和黃湘等[29]在塔里木河的研究結果一致。當土壤濕度在一個很小的范圍內(nèi)變化,土壤呼吸與土壤濕度之間將不會存在相關性,這是因為痕量般的水分變化可能會被其他因子所掩蓋[30]。本研究中,撂荒地地表植被稀疏、地上地下生物量較低、礫石含量高、土壤水分狀況較差,這都是限制土壤微生物活性和根呼吸的重要因子。當土壤含水量未超土壤水分閾值時,土壤微生物活性、根系呼吸也會出現(xiàn)正反饋。

    土壤CO2通量與濕度的關系比較復雜,包括了生物、物理、化學機制,并且隨著區(qū)域和時間尺度的不同而改變。土壤濕度的改變通常伴隨著土壤微生物活性、根以及溫度的變化,在干旱區(qū)還會影響到可溶性鹽分的分布。在極端條件下的鹽生荒漠系統(tǒng)中土壤經(jīng)常會出現(xiàn)“碳匯”和無機CO2交換的現(xiàn)象[31-32],目前對土壤無機CO2通量形成的機制仍沒有統(tǒng)一的定論,主要觀點有碳酸鹽系統(tǒng)的溶解/沉淀和CO2的溶解化學等。本研究區(qū)域位于達坂城鹽湖沿岸,土壤類型屬于鹽漬土,在監(jiān)測時段內(nèi),土壤未出現(xiàn)吸收CO2現(xiàn)象。

    4 結論

    在監(jiān)測時段內(nèi),不同植物群落和撂荒地土壤呼吸日變化均呈單峰曲線,土壤呼吸峰值出現(xiàn)在13:00左右。各群落間及撂荒地土壤CO2日排放總量存在顯著性差異,黑果枸杞群落土壤CO2日排放總量大于其他群落及撂荒地土壤;地上、地下生物量是影響土壤CO2排放總量的重要因子。土壤呼吸速率與土壤表層溫度存在顯著指數(shù)相關(2= 0.69 ~ 0.78),表明土壤溫度是影響土壤呼吸的環(huán)境因子;小獐毛群落土壤對溫度敏感強度大于其他土壤。不同樣地土壤呼吸與土壤濕度相關性顯著(2= 0.38 ~ 0.51),土壤濕度是決定達坂城鹽湖土壤呼吸速率的重要因子之一。

    鹽湖沿岸土壤與鹽生荒漠的土壤均為鹽漬化土壤,而二者的生態(tài)功能卻截然不同。我們的研究表明,在監(jiān)測時段內(nèi)土壤CO2未出現(xiàn)負通量,所以鹽湖周邊不同群落立地土壤表現(xiàn)為CO2的源,這可能與土壤鹽分類型以及土壤水分條件有關。要進一步準確評估達坂城鹽湖的生態(tài)系統(tǒng)碳收支狀況,仍需對各植物群落的固碳潛力進行研究,從而為鹽湖生態(tài)系統(tǒng)碳循環(huán)過程的系統(tǒng)研究提供科學數(shù)據(jù)。

    [1] 楊濤, 姜文波, 孫希春, 等. 不同植被類型高山草甸土呼吸作用的研究[J]. 土壤學報, 1989, 26(2): 186–192

    [2] Raich J W, Schlesinger W H. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate[J]. Ellus, 1992, 44(B): 81–99

    [3] 李彥, 王勤學, 馬健, 等. 鹽生荒漠地表水、熱CO2輸送的實驗研究[J]. 地理學報, 2004, 59(1): 33–39

    [4] 劉冉, 王勤學, 唐立松, 等. 鹽生荒漠地表水熱與二氧化碳通量的季節(jié)變化及驅動因素[J]. 生態(tài)學報, 2009, 29(1): 67–75

    [5] 劉冉, 李彥, 劉燕. 鹽生荒漠CO2通量對兩個極端降水年份響應[J]. 生態(tài)學雜志, 2013, 32(10): 2545–2551

    [6] 鄭喜玉, 張明剛, 徐昶, 等. 中國鹽湖志[M]. 北京: 科學出版社, 2002

    [7] 楊磊, 賈宏濤, 胡毅, 等. 水分條件對巴音布魯克高寒濕地CO2排放的影響[J].草業(yè)科學, 2015, 32(2): 168–173

    [8] 劉躍輝, 艾力·買買提依明, 楊帆, 等. 塔克拉瑪干沙漠腹地冬季土壤呼吸及其驅動因子[J]. 生態(tài)學報, 2015, 35(20): 6711–6719

    [9] 張增信, 施政, 何容, 等. 北亞熱帶次生櫟林和人工松林土壤呼吸日變化[J]. 南京林業(yè)大學學報(自然科學版), 2010, 34(1): 19–23

    [10] Bontti E E, Decant J P, Munson S M, et al. Litter decomposition in grasslands of Central North America[J]. Global Change Biology, 2009, 15(5): 1356–1363

    [11] Raich J W, Tufekcioglu A. Vegetation and soil respiration:Correlations and controls[J]. Biogechemistry, 2000, 48(1): 71–90

    [12] 馬杰, 于丹丹, 鄭新軍. 鹽生荒漠土壤CO2通量及其環(huán)境影響因素[J]. 生態(tài)學雜志, 2013, 32(10): 2532–2538

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

    [14] Jassal R, Black A, Novak M, et al. Relationship between soil CO2concentrations and forest-floor CO2effluxes[J]. Agricultural and Forest Meteorology, 2005, 130(3): 176–192

    [15] Tang J, Baldocchi D D. Spatial-temporal variation in soil respiration in an oak-grass savanna ecosystem in California and its partitioning into autotrophic and heterotrophic components[J]. Biogeochemistry, 2005, 73(1):183–207

    [16] 董自紅, 蔣平安, 賈宏濤, 等. 新疆高山草甸草原植物群落夏季CO2排放日變化分析[J].新疆農(nóng)業(yè)科學, 2007, 44(1): 1–5

    [17] 王憶慧, 龔吉蕊, 劉敏, 等. 草地利用方式對土壤呼吸和凋落物分解的影響[J]. 植物生態(tài)學報, 2015, 39(3): 239–248

    [18] 楊建軍, 呂光輝, 張燕, 等. 艾比湖流域土壤呼吸日變化及水熱因子影響[J]. 新疆農(nóng)業(yè)科學, 2009, 46(2): 223–231

    [19] 王春燕, 陳秋波, 袁坤, 等. 橡膠林土壤呼吸速率及其與土壤溫濕度的關系[J]. 土壤學報, 2013, 50(5):119–128

    [20] 張賽, 王龍昌, 張曉雨, 等. 紫色土丘陵區(qū)保護性耕作下旱地土壤呼吸及影響因素[J]. 土壤學報, 2014, 51(3):520–530

    [21] 王忠媛, 謝江波, 王玉剛, 等. 溫度對鹽土和堿土土壤無機CO2通量的影響[J]. 生態(tài)學雜志, 2013, 32(10): 2525–2531

    [22] 全權, 張震, 何念鵬, 等. 短期氮添加對東靈山三種森林土壤呼吸的影響[J]. 生態(tài)學雜志, 2015, 34(3): 797–804

    [23] Luo Y Q, Wan S Q, Hui D F, et al. Acclimatization of soil respiration to warming in a tall grass prairie[J]. Nature, 2001, 413: 622–625

    [24] 陳寶玉, 劉世榮, 葛建平, 等. 川西高山針葉林土壤呼吸速率與不同土層溫度的關系[J]. 應用生態(tài)學報, 2007,18(6):1219–1224

    [25] 張盼弟, 王旭, 陳寶瑞, 等. 呼倫貝爾地區(qū)貝加爾針茅草甸草原土壤CO2排放特征[J]. 應用生態(tài)學報, 2014, 25(2): 387–393

    [26] 孔雨光, 張金池, 王因花, 等. 蘇北淤泥質海岸典型防護林地土壤呼吸及其溫度敏感性[J]. 生態(tài)學報, 2009, 29(8): 4084–4093

    [27] 高艷紅, 張志山, 劉立超, 等. 騰格里沙漠植被重建對土壤呼吸的影響. 土壤學報, 2009, 46(4):627–633

    [28] 方晰, 田大倫, 項文化. 杉木人工林林地土壤CO2釋放量及其影響因子的研究. 林業(yè)科學, 2005, 41(2): 1–7

    [29] 黃湘, 李衛(wèi)紅, 陳亞寧, 等. 塔里木河下游荒漠河岸林群落土壤呼吸及其影響因子[J]. 生態(tài)學報, 2007, 27(5): 1952–1959

    [30] 李虎, 邱建軍, 王立剛. 農(nóng)田土壤呼吸特征及根呼吸貢獻的模擬分析[J]. 農(nóng)業(yè)工程學報, 2008, 24(4): 14–20

    [31] Liu J B, Feng W, Zhang Y Q. Abiotic CO2exchange between soil and atmosphere and its response to temperature[J]. Environmental Earth Sciences, 2015, 73: 2463–2471

    [32] Emma L, Angela M, Detweiler L T. Assessing the role of alkaline soils on the carbon cycle at a playa site[J]. Environmental Earth Sciences, 2013, 70(3): 1047–1056

    Soil Respirations Under Different Plant Communities in Dabancheng Salt Lake of Xinjiang

    LI Dianpeng, YAO Meisi, HAN Dongliang, WANG Liping, HU Bao’an, WANG Ningning, JIANG Dayong, ZHANG Kai, ZHU Xinping, JIA Hongtao*

    (College of Grassland and Environment Sciences, Xinjiang Agricultural University, Urumqi 830052, China)

    The distribution of salt lakes in Xinjiang is broad but uneven. The carbon cycle of the salt lake ecosystem is an important component and has an important impact on carbon equilibrium in the ecosystem of the arid region. However, the research on soil respiration in salt lakes in the arid region is relative weak. A study was conducted in order to explore soil respiration and its influential factors of different plant communities along the Salt Lake in Dabancheng of Xinjiang.,,,were taken as the main plant objects with the abandoned farmland around as the control, soil respiration were measured by LI-COR 8100 (Li-COR, Lincoln, NE, USA). The results showed soil respiration of Salt Lake coastal had obvious diurnal variation and appeared in the pattern of single peak, the highest point occurred at 13:00 while the lowest occurred at around 07:00. The rates of soil respiration ranged from 1.77 to 3.48 μmol/(m2·s) in CK and from 0.89 to 4.34 μmol/(m2·s) under different communities, with an order of. Underground biomass influenced soil CO2emissions, the daily average CO2emission flux of soil respiration in the study region was 209.28 mg/(m2·d). Soil temperature at 5 cm depth had significant correlation with soil respiration rate, which is the dominant factor in determining the diurnal variation of soil respiration and could explain 69%-78% of it. The10values calculated based on soil temperature at 5 cm depth were 2.30, 1.24, 1.20, 1.57, and 1.28, soil respiration was most sensitive to soil temperature under. Except for abandoned land, there were significantly negative correlation (2= 0.38-0.51) between soil respiration and soil moisture. The above results disclosed further soil respiration variation and influential factors of Slat Lake ecosystem in the arid region, which is helpful in estimating carbon equilibrium and sequestration.

    Arid region; Salt take; Carbon cycle; Soil respiration; Influential factors

    10.13758/j.cnki.tr.2017.03.028

    S154.1

    A

    國家自然科學基金項目(31560171)和國家大學生創(chuàng)新訓練計劃項目(201510758004)資助。

    (hongtaojia@126.com)

    李典鵬(1992—),男,湖南新華人,在讀本科生,主要研究領域為土壤生態(tài)系統(tǒng)物質循環(huán)。E-mail: LLdp05120@126.com

    猜你喜歡
    芨芨草鳶尾土壤濕度
    “人間彩虹”鳶尾
    竹子和芨芨草
    鳶尾,只綻放一天的彩虹女神
    中老年保健(2021年7期)2021-08-22 07:44:28
    土壤濕度傳感器在園林綠化灌溉上的應用初探
    鳶尾素與惡性腫瘤相關研究進展
    甘肅科技(2020年20期)2020-04-13 00:30:48
    基于51單片機控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    鳶尾苷元在兔體內(nèi)的藥動學
    中成藥(2017年10期)2017-11-16 00:49:54
    甘肅省甘谷縣第三中學芨芨草文學社作品小輯
    芨芨草文學社
    大作文(2016年11期)2016-05-30 00:25:54
    四川盆地土壤濕度時空分布及影響因子分析
    狂野欧美激情性bbbbbb| 男女之事视频高清在线观看| 久久99热这里只频精品6学生| 欧美精品一区二区免费开放| 99国产精品一区二区三区| 99精品久久久久人妻精品| 国产在线观看jvid| 真人做人爱边吃奶动态| 午夜福利在线免费观看网站| 午夜福利视频在线观看免费| 亚洲一区中文字幕在线| 成人国产av品久久久| 午夜久久久在线观看| 亚洲一码二码三码区别大吗| 日本精品一区二区三区蜜桃| 飞空精品影院首页| 国精品久久久久久国模美| 国产欧美日韩一区二区精品| 成人亚洲精品一区在线观看| 日韩欧美国产一区二区入口| 高清欧美精品videossex| 成人三级做爰电影| 欧美在线一区亚洲| 国产一区二区激情短视频 | 久久久久国内视频| 日韩欧美一区二区三区在线观看 | 99久久99久久久精品蜜桃| 午夜精品久久久久久毛片777| 色94色欧美一区二区| 91字幕亚洲| 久久国产精品人妻蜜桃| 香蕉国产在线看| 亚洲国产日韩一区二区| 999久久久精品免费观看国产| 黄片播放在线免费| tocl精华| 久久精品人人爽人人爽视色| 成年av动漫网址| 午夜福利视频精品| 老司机福利观看| 久久久久国内视频| 99国产精品免费福利视频| 国产主播在线观看一区二区| 精品人妻一区二区三区麻豆| 免费在线观看日本一区| 人人妻人人澡人人看| 国产免费现黄频在线看| 国产亚洲av片在线观看秒播厂| 欧美日韩黄片免| 日韩免费高清中文字幕av| 搡老乐熟女国产| 美女扒开内裤让男人捅视频| 99精品久久久久人妻精品| 成年美女黄网站色视频大全免费| 日韩欧美一区二区三区在线观看 | 国产不卡av网站在线观看| 国产深夜福利视频在线观看| 黄色视频在线播放观看不卡| 一区二区三区四区激情视频| 亚洲国产中文字幕在线视频| 操出白浆在线播放| 最近最新中文字幕大全免费视频| www.999成人在线观看| 十八禁高潮呻吟视频| 亚洲专区国产一区二区| 午夜福利影视在线免费观看| 久久久久久久精品精品| 各种免费的搞黄视频| 两性夫妻黄色片| 国产精品一区二区在线观看99| 亚洲国产精品999| 国产一区有黄有色的免费视频| 一边摸一边做爽爽视频免费| 免费一级毛片在线播放高清视频 | 秋霞在线观看毛片| 久久综合国产亚洲精品| 青草久久国产| 18禁观看日本| e午夜精品久久久久久久| 国产成人精品无人区| 久久性视频一级片| 极品人妻少妇av视频| 成人18禁高潮啪啪吃奶动态图| 国产亚洲欧美精品永久| 一进一出抽搐动态| 又大又爽又粗| 美女脱内裤让男人舔精品视频| 国产日韩欧美在线精品| 999久久久精品免费观看国产| 免费日韩欧美在线观看| av在线播放精品| 爱豆传媒免费全集在线观看| 日本av免费视频播放| 亚洲av日韩在线播放| 午夜福利在线免费观看网站| h视频一区二区三区| 伊人久久大香线蕉亚洲五| 欧美人与性动交α欧美软件| 一级片'在线观看视频| 欧美日韩福利视频一区二区| 亚洲人成77777在线视频| 91成人精品电影| 亚洲男人天堂网一区| 亚洲精品av麻豆狂野| 日本猛色少妇xxxxx猛交久久| videos熟女内射| 亚洲欧洲精品一区二区精品久久久| 国产熟女午夜一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美黄色片欧美黄色片| 少妇精品久久久久久久| 999精品在线视频| 麻豆乱淫一区二区| 国产精品一区二区精品视频观看| 99九九在线精品视频| 真人做人爱边吃奶动态| 国产99久久九九免费精品| 国产精品一区二区在线不卡| 操美女的视频在线观看| 欧美日韩一级在线毛片| 如日韩欧美国产精品一区二区三区| 国产黄色免费在线视频| 精品福利永久在线观看| cao死你这个sao货| 视频在线观看一区二区三区| 亚洲伊人久久精品综合| 久久这里只有精品19| 一进一出抽搐动态| 大香蕉久久成人网| 18在线观看网站| 久久久久国内视频| cao死你这个sao货| tocl精华| 69精品国产乱码久久久| 波多野结衣一区麻豆| 乱人伦中国视频| 日韩中文字幕欧美一区二区| 老司机亚洲免费影院| 人妻人人澡人人爽人人| 热re99久久国产66热| 精品少妇内射三级| 亚洲欧美一区二区三区久久| 午夜91福利影院| 日韩免费高清中文字幕av| 午夜日韩欧美国产| 9191精品国产免费久久| 91精品伊人久久大香线蕉| 黄色怎么调成土黄色| 亚洲av电影在线观看一区二区三区| 法律面前人人平等表现在哪些方面 | 久久精品熟女亚洲av麻豆精品| 天天躁夜夜躁狠狠躁躁| 999精品在线视频| 丝袜在线中文字幕| 成人av一区二区三区在线看 | 欧美xxⅹ黑人| videosex国产| 亚洲人成77777在线视频| 欧美黑人欧美精品刺激| 日韩电影二区| 天天躁狠狠躁夜夜躁狠狠躁| 成人国产一区最新在线观看| 另类精品久久| 搡老乐熟女国产| 国产国语露脸激情在线看| 性高湖久久久久久久久免费观看| 岛国毛片在线播放| 国产又爽黄色视频| 国产精品香港三级国产av潘金莲| 在线观看免费视频网站a站| 国产主播在线观看一区二区| 一区福利在线观看| 一区二区日韩欧美中文字幕| 欧美在线黄色| 男女无遮挡免费网站观看| 欧美变态另类bdsm刘玥| 精品一区二区三卡| 欧美97在线视频| 美女中出高潮动态图| 国产欧美日韩一区二区三 | av有码第一页| 成在线人永久免费视频| 香蕉丝袜av| 亚洲国产精品一区二区三区在线| 欧美精品人与动牲交sv欧美| 最近最新免费中文字幕在线| a在线观看视频网站| 久久性视频一级片| 日本猛色少妇xxxxx猛交久久| 一级片'在线观看视频| 中文精品一卡2卡3卡4更新| 高清在线国产一区| 国产99久久九九免费精品| 看免费av毛片| 日日爽夜夜爽网站| 午夜91福利影院| 老司机亚洲免费影院| 秋霞在线观看毛片| 秋霞在线观看毛片| 老熟妇仑乱视频hdxx| 男女午夜视频在线观看| 久久精品aⅴ一区二区三区四区| 久久精品熟女亚洲av麻豆精品| 在线观看免费视频网站a站| 国产淫语在线视频| 美女高潮喷水抽搐中文字幕| 国产精品二区激情视频| 国产亚洲欧美在线一区二区| 日本a在线网址| 久久性视频一级片| 三级毛片av免费| 夫妻午夜视频| 午夜免费鲁丝| 亚洲欧美清纯卡通| 又紧又爽又黄一区二区| 纯流量卡能插随身wifi吗| 天天影视国产精品| 男女床上黄色一级片免费看| 99精国产麻豆久久婷婷| 精品一区二区三区四区五区乱码| 国产片内射在线| 天堂8中文在线网| 精品卡一卡二卡四卡免费| 性高湖久久久久久久久免费观看| 日韩欧美免费精品| 国产色视频综合| 日韩制服骚丝袜av| svipshipincom国产片| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 丝袜喷水一区| 国产精品二区激情视频| 交换朋友夫妻互换小说| 精品国产乱码久久久久久小说| 婷婷丁香在线五月| 啦啦啦中文免费视频观看日本| 男人爽女人下面视频在线观看| 丝袜在线中文字幕| 啦啦啦免费观看视频1| 高清视频免费观看一区二区| 国产有黄有色有爽视频| 黄片大片在线免费观看| 国产亚洲欧美精品永久| 我的亚洲天堂| 色老头精品视频在线观看| 飞空精品影院首页| 成年动漫av网址| 一本一本久久a久久精品综合妖精| 国产男女超爽视频在线观看| 欧美亚洲 丝袜 人妻 在线| 一进一出抽搐动态| 午夜激情久久久久久久| 中文字幕av电影在线播放| 欧美精品亚洲一区二区| 丝瓜视频免费看黄片| 国产伦理片在线播放av一区| 国产亚洲精品久久久久5区| 久久久国产成人免费| 中文字幕最新亚洲高清| 在线十欧美十亚洲十日本专区| 成年av动漫网址| 久久亚洲精品不卡| 国产av精品麻豆| 国产一区二区三区av在线| 99热国产这里只有精品6| 久久精品亚洲av国产电影网| 亚洲第一欧美日韩一区二区三区 | 欧美日韩视频精品一区| 亚洲国产av新网站| 国产一区二区 视频在线| 少妇猛男粗大的猛烈进出视频| 久久久精品区二区三区| 午夜免费成人在线视频| av欧美777| 他把我摸到了高潮在线观看 | 十八禁高潮呻吟视频| 午夜福利一区二区在线看| 国产一区二区三区综合在线观看| 午夜福利在线免费观看网站| 黄片大片在线免费观看| 亚洲免费av在线视频| 久久热在线av| 午夜福利影视在线免费观看| 一二三四在线观看免费中文在| 老熟妇仑乱视频hdxx| 欧美午夜高清在线| 欧美激情高清一区二区三区| 免费人妻精品一区二区三区视频| 老司机午夜十八禁免费视频| 91大片在线观看| 咕卡用的链子| 久久久久久免费高清国产稀缺| 无遮挡黄片免费观看| 在线精品无人区一区二区三| 日韩人妻精品一区2区三区| 色婷婷av一区二区三区视频| 久久综合国产亚洲精品| a级毛片黄视频| 91成人精品电影| 十八禁网站免费在线| 亚洲黑人精品在线| 亚洲精品乱久久久久久| av国产精品久久久久影院| 国产视频一区二区在线看| 亚洲精品久久成人aⅴ小说| 在线av久久热| 在线观看免费高清a一片| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三 | 国产精品秋霞免费鲁丝片| 淫妇啪啪啪对白视频 | 精品人妻一区二区三区麻豆| 高清视频免费观看一区二区| 欧美xxⅹ黑人| av片东京热男人的天堂| xxxhd国产人妻xxx| 一级a爱视频在线免费观看| 美女中出高潮动态图| 在线看a的网站| 狂野欧美激情性xxxx| 一本综合久久免费| 久久久久久久国产电影| 精品乱码久久久久久99久播| 男女免费视频国产| 亚洲国产欧美在线一区| 9热在线视频观看99| 久久人人爽人人片av| 日韩电影二区| 天天操日日干夜夜撸| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av男天堂| 亚洲中文日韩欧美视频| 水蜜桃什么品种好| 久久久久国产精品人妻一区二区| 午夜福利视频精品| 亚洲伊人久久精品综合| 汤姆久久久久久久影院中文字幕| 国产又色又爽无遮挡免| 亚洲精品一卡2卡三卡4卡5卡 | 久久久国产欧美日韩av| 99国产精品免费福利视频| 国产色视频综合| 亚洲情色 制服丝袜| 黄网站色视频无遮挡免费观看| 国产av一区二区精品久久| 性色av乱码一区二区三区2| 在线av久久热| 久久精品亚洲熟妇少妇任你| 精品亚洲成a人片在线观看| 黑丝袜美女国产一区| 夜夜夜夜夜久久久久| 一本大道久久a久久精品| 亚洲欧美色中文字幕在线| 香蕉丝袜av| 首页视频小说图片口味搜索| 老司机深夜福利视频在线观看 | 国产av又大| 久久99热这里只频精品6学生| 不卡一级毛片| 亚洲av美国av| 免费高清在线观看视频在线观看| 国产人伦9x9x在线观看| 一级a爱视频在线免费观看| 两个人免费观看高清视频| 久久精品成人免费网站| √禁漫天堂资源中文www| 一边摸一边抽搐一进一出视频| 欧美黑人精品巨大| a级毛片在线看网站| 91大片在线观看| 成人av一区二区三区在线看 | 91老司机精品| 超碰97精品在线观看| 老汉色av国产亚洲站长工具| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 亚洲精品国产色婷婷电影| 一本综合久久免费| 国产av一区二区精品久久| 亚洲av男天堂| 国产野战对白在线观看| 大码成人一级视频| 久久精品aⅴ一区二区三区四区| 欧美另类一区| 国产91精品成人一区二区三区 | 一区二区三区激情视频| 免费观看a级毛片全部| 飞空精品影院首页| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻一区二区三区麻豆| 91麻豆av在线| 在线永久观看黄色视频| 色视频在线一区二区三区| 亚洲国产日韩一区二区| 久久精品国产综合久久久| 精品国产乱码久久久久久男人| 欧美 亚洲 国产 日韩一| 免费观看av网站的网址| 欧美精品亚洲一区二区| 久久精品亚洲av国产电影网| 午夜精品国产一区二区电影| 国产区一区二久久| 成年动漫av网址| 国产片内射在线| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 久久精品人人爽人人爽视色| 我要看黄色一级片免费的| 最近最新中文字幕大全免费视频| 精品少妇一区二区三区视频日本电影| 亚洲熟女精品中文字幕| 久久人人97超碰香蕉20202| 免费不卡黄色视频| 国产淫语在线视频| 9色porny在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产片内射在线| 午夜福利,免费看| 窝窝影院91人妻| 91成年电影在线观看| 一本综合久久免费| 亚洲国产欧美一区二区综合| 一区二区三区激情视频| 人人妻人人澡人人看| 最近中文字幕2019免费版| 美国免费a级毛片| 天堂中文最新版在线下载| 久久久久网色| 免费高清在线观看视频在线观看| 狠狠精品人妻久久久久久综合| av片东京热男人的天堂| 午夜福利乱码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线免费观看网站| 色老头精品视频在线观看| 久久精品成人免费网站| 亚洲国产成人一精品久久久| 女人精品久久久久毛片| 人妻 亚洲 视频| 亚洲免费av在线视频| 日韩电影二区| 50天的宝宝边吃奶边哭怎么回事| 久久久国产精品麻豆| 夜夜骑夜夜射夜夜干| 伊人亚洲综合成人网| 性色av一级| 99久久国产精品久久久| 精品国产一区二区三区四区第35| 久久热在线av| 久久国产精品大桥未久av| 女性被躁到高潮视频| 亚洲精品自拍成人| 国产欧美日韩一区二区三区在线| 成人国产av品久久久| 国产精品九九99| 在线观看www视频免费| 丝袜人妻中文字幕| 黑人巨大精品欧美一区二区mp4| 制服诱惑二区| 亚洲成人免费av在线播放| 多毛熟女@视频| 精品少妇黑人巨大在线播放| 亚洲av片天天在线观看| 国产一区二区三区在线臀色熟女 | 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 三级毛片av免费| 国产xxxxx性猛交| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区综合在线观看| 又黄又粗又硬又大视频| 超色免费av| 男人舔女人的私密视频| 777米奇影视久久| 久久亚洲国产成人精品v| 午夜福利视频精品| 精品福利永久在线观看| 少妇粗大呻吟视频| 亚洲国产欧美在线一区| 欧美黑人欧美精品刺激| 一区二区三区乱码不卡18| 精品福利永久在线观看| 欧美 日韩 精品 国产| 一本色道久久久久久精品综合| 国产精品国产av在线观看| xxxhd国产人妻xxx| 精品免费久久久久久久清纯 | 欧美日韩av久久| 成人手机av| av在线app专区| 丰满人妻熟妇乱又伦精品不卡| 日韩三级视频一区二区三区| 不卡av一区二区三区| 一级,二级,三级黄色视频| 午夜视频精品福利| 一级片'在线观看视频| 中文字幕人妻丝袜一区二区| 又黄又粗又硬又大视频| 91国产中文字幕| 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 视频区欧美日本亚洲| 下体分泌物呈黄色| 久久香蕉激情| 久久国产精品男人的天堂亚洲| 黄片播放在线免费| 狠狠婷婷综合久久久久久88av| 日韩一区二区三区影片| av片东京热男人的天堂| 宅男免费午夜| 午夜91福利影院| 国产黄色免费在线视频| 伊人久久大香线蕉亚洲五| 免费久久久久久久精品成人欧美视频| www.av在线官网国产| 久久精品亚洲熟妇少妇任你| 天堂8中文在线网| 精品国产一区二区三区久久久樱花| 一级a爱视频在线免费观看| 女警被强在线播放| 人妻 亚洲 视频| 在线观看免费高清a一片| 人成视频在线观看免费观看| 纯流量卡能插随身wifi吗| 大片电影免费在线观看免费| 久久久水蜜桃国产精品网| 亚洲自偷自拍图片 自拍| 国产色视频综合| 欧美av亚洲av综合av国产av| 50天的宝宝边吃奶边哭怎么回事| 天堂中文最新版在线下载| 婷婷色av中文字幕| 国产男人的电影天堂91| 19禁男女啪啪无遮挡网站| 曰老女人黄片| 热99re8久久精品国产| 国产av一区二区精品久久| 热re99久久国产66热| 1024视频免费在线观看| 精品人妻一区二区三区麻豆| 丰满迷人的少妇在线观看| 91av网站免费观看| 欧美日韩黄片免| 婷婷丁香在线五月| 亚洲精品国产精品久久久不卡| 下体分泌物呈黄色| 一区二区av电影网| 18禁国产床啪视频网站| 无限看片的www在线观看| 黑丝袜美女国产一区| 秋霞在线观看毛片| 欧美黄色淫秽网站| 国产91精品成人一区二区三区 | 国产日韩欧美在线精品| 亚洲 欧美一区二区三区| 大型av网站在线播放| 午夜视频精品福利| 多毛熟女@视频| 老司机影院毛片| 亚洲中文字幕日韩| 亚洲第一欧美日韩一区二区三区 | 最近最新中文字幕大全免费视频| 中文字幕制服av| 一级毛片电影观看| 丝袜美腿诱惑在线| 国产一级毛片在线| 亚洲第一欧美日韩一区二区三区 | 亚洲av日韩精品久久久久久密| 亚洲专区国产一区二区| 亚洲av男天堂| 一二三四社区在线视频社区8| 国产精品久久久久久精品电影小说| 人人妻人人澡人人看| kizo精华| 免费在线观看视频国产中文字幕亚洲 | 亚洲av成人一区二区三| 国产黄色免费在线视频| av在线老鸭窝| 窝窝影院91人妻| 亚洲激情五月婷婷啪啪| 国产亚洲av片在线观看秒播厂| 久久久久久久精品精品| 一级片免费观看大全| 国产日韩欧美亚洲二区| 亚洲精品久久午夜乱码| 亚洲国产欧美一区二区综合| 精品人妻一区二区三区麻豆| 高潮久久久久久久久久久不卡| 69av精品久久久久久 | 91av网站免费观看| 女人被躁到高潮嗷嗷叫费观| 久久久精品94久久精品| 香蕉丝袜av| 在线亚洲精品国产二区图片欧美| 国产高清videossex| 少妇裸体淫交视频免费看高清 | 后天国语完整版免费观看| 美女午夜性视频免费| 亚洲av片天天在线观看| 亚洲美女黄色视频免费看| h视频一区二区三区| 在线观看舔阴道视频| 国产熟女午夜一区二区三区| 日韩欧美一区二区三区在线观看 | e午夜精品久久久久久久| 精品人妻熟女毛片av久久网站| 亚洲一区中文字幕在线| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕人妻丝袜一区二区| 欧美精品亚洲一区二区| 一本—道久久a久久精品蜜桃钙片| 老汉色av国产亚洲站长工具| 久久国产精品影院| 91麻豆av在线| 成在线人永久免费视频|