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

    不同經(jīng)營方式及生境對大興安嶺高緯度林區(qū)生長盛季森林土壤CO2、CH4、N2O通量的影響

    2015-12-06 06:45:30任樂馬秀枝范雪松
    生態(tài)環(huán)境學報 2015年3期
    關鍵詞:興安土壤濕度落葉松

    任樂,馬秀枝,范雪松

    1. 內蒙古農業(yè)大學生態(tài)環(huán)境學院,內蒙古 呼和浩特 100093;2. 內蒙古農業(yè)大學林學院,內蒙古 呼和浩特 100093

    不同經(jīng)營方式及生境對大興安嶺高緯度林區(qū)生長盛季森林土壤CO2、CH4、N2O通量的影響

    任樂1,馬秀枝2*,范雪松1

    1. 內蒙古農業(yè)大學生態(tài)環(huán)境學院,內蒙古 呼和浩特 100093;2. 內蒙古農業(yè)大學林學院,內蒙古 呼和浩特 100093

    為了確定不同森林經(jīng)營利用方式(原始林、皆伐林、漸伐林)及不同生境(有林地、無林地、濕地)對于森林土壤-大氣主要溫室氣體通量的影響,于2013年生長盛季(7─8月),在內蒙古大興安嶺興安落葉松(Larix gmelinii)林區(qū)選取了16塊樣地,采用靜態(tài)箱/氣相色譜法進行了土壤-大氣溫室氣體通量的原位觀測。氣體監(jiān)測的同時,測定了土壤濕度和各層土壤溫度以及各樣地土壤的物理、化學性質。實驗結果表明:大興安嶺興安落葉松林區(qū)生長盛季 CH4通量平均值為C(-133.6±62.3) μg·m-2·h-1。沼澤濕地(PD2和PD7)CH4表現(xiàn)為釋放,CH4通量平均值為C(130.0±41.0) μg·m-2·h-1,其它14個樣地土壤CH4為吸收,吸收通量變化范圍為C(-242.0~-42.6) μg·m-2·h-1,變異系數(shù)為47%。所有樣地CO2的地-氣通量變化范圍是 C(94.9~1164.0) μg·m-2·h-1,平均值為 C(671.3±324.3) μg·m-2·h-1,變異系數(shù)是 48%。所有樣地 N2O通量的范圍是N(1.2~21.6) μg·m-2·h-1,平均值為N(7.3±6.1) μg·m-2·h-1,變異系數(shù)為84%。通過數(shù)據(jù)分析顯示,大興安嶺林區(qū)生長盛季不同經(jīng)營方式及不同生境對森林土壤-大氣CH4、CO2、N2O通量均無顯著相關性。CH4通量與土壤各層溫度及各層濕度并未發(fā)現(xiàn)有顯著相關性,僅顯示CH4通量與土壤有機質(SOC)呈正相關關系、土壤全氮(TN)呈負相關關系;CO2通量與土壤5 cm處溫度顯著相關(P=0.042 3),與土壤0 cm及10 cm處溫度無顯著相關性,CO2通量與土壤各層溫度呈正相關關系,與腐殖質層、0~10 cm 以及10~20 cm土壤濕度呈負相關關系,與SOC、TN呈負相關關系;N2O通量與土壤各層溫度、濕度及SOC均呈正相關關系,與TN呈負相關關系。

    興安落葉松林;CO2、CH4、N2O;不同經(jīng)營方式;不同生境

    陸地生態(tài)系統(tǒng)是溫室氣體CH4、CO2和N2O主要的源和匯(Chapius-Lardy等,2007)。森林作為陸地生態(tài)系統(tǒng)的重要組成部分,占陸地總面積的1/3(4.1×109hm2),在我國森林面積占土地面積的23.7%(齊玉春等,2002)。目前國內外有關森林生態(tài)系統(tǒng)土壤溫室氣體通量的研究大多集中在熱帶、亞熱帶及溫帶地區(qū)(孫向陽和徐化成,2001;沙麗清等,2004;莫江明等,2006)。大興安嶺北部林區(qū)位于寒溫帶,是我國面積最大的原始林區(qū)。興安落葉松(Larix gmelinii)是該區(qū)主要的地帶性植被,林分覆蓋面積為1.56×107hm2,占全國森林面積的13.2%(蔣延玲和周廣勝,1999),對溫室氣體的排放和吸收有重要的影響作用。

    眾所周知土壤中的氣體的產生和消耗都是由于微生物活動的結果,但是土壤與大氣之間的通量卻與土壤理化性質參數(shù)有關,比如N含量、碳氮比、溫度以及含水量(Smith等,2003;Lilly等,2009)。研究發(fā)現(xiàn)成熟云杉林(Picea abies(L.)Karst.)中,在影響土壤CO2排放的所有非生物因素中,土壤濕度變化是影響土壤 CO2排放變化的極其重要因素(Subke等,2003)。還有研究指出溫度是影響土壤CO2排放的最主要因素之一,兩者之間有較明顯的關系,許多的實驗結果表明溫度與 CO2排放有顯著的指數(shù)函數(shù)相關關系(Raich和 Schlesinger,1992),其中5和10 cm土層溫度對土壤CO2排放影響最大。

    Li等(2008)的研究表明,在不同緯度、不同海拔地區(qū),由于土壤溫度和濕度的差異,以及不同植被條件的影響,土壤地表CO2通量的變化范圍很大。肖冬梅等(2004)和傅民杰等(2009)的研究也表明,不同緯度地區(qū)N2O的釋放也有顯著差異,其平均通量和累積釋放量均隨緯度升高而下降。受土壤溫度和水分的影響,森林土壤CH4、CO2和N2O通量往往呈現(xiàn)明顯的季節(jié)性變化。

    因此,為了更好的了解大興安嶺林區(qū)不同森林經(jīng)營方式及不同生境下土壤-大氣CH4、CO2和N2O通量的空間變異情況,本研究于2013年7月19日至8月15日期間,對內蒙古大興安嶺高緯度興安落葉松(Larix gmelinii)林區(qū)土壤開展了CH4、CO2和N2O通量的原位觀測研究,確定大興安嶺高緯度興安落葉松林區(qū)土壤-大氣溫室氣體通量的空間變異,為該地區(qū)溫室氣體通量的總體核算提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗地處內蒙古海拉爾市根河市,依托內蒙古大興安嶺森林生態(tài)系統(tǒng)國家野外科學觀測研究站試驗區(qū)(50°49′~50°51′N,121°30′~121°31′E),試驗區(qū)面積1.1萬 hm2,其中原始林3200 hm2。該地區(qū)屬寒溫帶半濕潤氣候區(qū)。冬季(平均氣溫<10 ℃)長達9個月,夏季(平均氣溫≥22 ℃)不超過1個月。≥10 ℃年積溫1403 ℃,全年最高氣溫40 ℃,極端最低氣溫-50 ℃,全年平均溫度為-5.4 ℃。年降水量450~550 mm,其中60%集中在7月、8月。9月末至次年5月初為降雪期,降雪厚度20~40 cm,降雪量占全年降水總量的 12%。全年地表蒸發(fā)量800~1200 mm,年均日照2594 h,全年無霜期80 d。土壤凍結狀態(tài)從9月下旬到次年的5月初,凍結時間長達7個多月。該區(qū)土壤以棕色針葉林土為主,土壤表層有較厚的枯枝落葉層,厚度達5~8 cm。土壤pH為4.5~6.5,腐殖質層一般在10 cm左右。由于當?shù)睾錆駶?,苔蘚枯枝落葉層發(fā)育較厚,滯水性強,使土壤在生長季處于濕潤狀態(tài),灰化現(xiàn)象不明顯,土層30~40 cm石礫較多。興安落葉松林郁閉度為0.8。

    1.2 實驗設置

    于2013年生長季7月19日─8月15日,在興安落葉林3種森林經(jīng)營利用方式下,又分別按照不同坡度和坡位共選取了 16塊樣地,分別涵蓋該林區(qū)的原始林、漸伐林及皆伐林3種不同經(jīng)營方式,人工林、無林地及沼澤3種不同生境。其中包括原始林4塊樣地,即PD1、LPL-Y、LPL-M1和LPL-N;漸伐林兩塊樣地,PD3和PD4;皆伐林6塊樣地,PD5、PD6、PD8、HL-M、HL-N和HL-O;人工林1塊樣地,PD9;無林地1塊樣地,PD10;沼澤濕地2塊樣地,PD2和PD7,樣地具體信息見表1。

    表1 興安落葉松林樣地基本信息Table 1 The basic information of sample plots in Larix gmelinii Forest

    在面積為10 m×10 m的樣地中,選擇3個間隔相等的點,并且要求3個點能代表該樣地的基本情況。在每個樣方內,于試驗開始前一周將采樣箱基座插入采樣點,盡量不破壞原有植被,以便使人為擾動對微量氣體交換的影響降到最低。各樣地基本情況介紹如下:

    興安落葉松研究區(qū)內選擇3種不同森林經(jīng)營利用方式,即原始林、漸伐林、皆伐林以及3種不同的生境,即人工林、無林地和沼澤濕地內分別設置樣地。進行土壤-大氣溫室氣體通量(CO2、CH4、N2O)的原位觀測。

    原始林:喬木層結構簡單,主要為興安落葉松林,有時混生少量白樺,郁閉度0.6~0.8。灌木層以杜香(palustre)占絕對優(yōu)勢,間有聚集分布的杜斯越橘和越橘(Vaccinium vitisidaed),高40~50 cm,層蓋度76%以上,因杜香等灌木根系在土壤中形成較厚的盤結層,草本植物難以生長,只有少數(shù)的小葉章(Deyeuxia -anguttifolia)、紅花鹿蹄草(incamatd)、苔草(Carex spp)等。蘚被層較發(fā)達,蓋度40%以上,厚度5~15 cm,海拔838 m。

    漸伐林:喬木層結構簡單,主要為興安落葉松林,有時混生少量白禪,郁閉度 0.6。灌木層以杜香為主,間有聚集分布的杜斯越橘和越橘,高10~20 cm,層蓋度95%以上,因杜香等灌木根系在土壤中形成較厚的盤結層,草本植物難以生長,只有少數(shù)的小葉章、紅花鹿蹄草、笞草等。蘚被層較發(fā)達,蓋度40%以上,厚度5~15 cm。海拔822 m。

    皆伐林:喬木層結構簡單,主要為興安落葉松林,有時混生少量白樺,郁閉度0.7。灌木層以杜香占絕對優(yōu)勢,間有聚集分布的杜斯越橘和越橘,高10~20 cm,層蓋度85%,因杜香等灌木根系在土壤中形成較厚的盤結層,草本植物難以生長,只有少數(shù)的小葉章、紅花鹿蹄草、苔草等。蘚被層較發(fā)達,蓋度40%以上,厚度5~15 cm。海拔825 m。

    人工林:喬木層結構簡單,主要為興安落葉松林,有時泥生少量白禪,郁閉度 0.6。灌木層以杜香占絕對優(yōu)勢,間有聚集分布的杜斯越橘和越橘,高40~50 cm,層蓋度60%以上,因杜香等灌木根系在土壤中形成較厚的盤結層,草本植物難以生長,只有少數(shù)的小葉章、紅花鹿蹄草、苔草等。蘚被層較發(fā)達,蓋度40%以上,厚度5~15 cm。海拔781 m。

    無林地:過度采伐跡地。沒有喬灌木。以草本為主,蓋度30%。海拔734 m

    沼澤地:主要河谷低洼地區(qū),長期積水,喬木樹種少,只有凸處少有落葉松林,無灌木,草本層,塔頭苔草為主,蓋度80%,海拔734 m。

    1.3 樣品采集和氣體通量測定

    采用靜態(tài)箱氣相色譜法,靜態(tài)箱由厚 2.0 mm的不銹鋼板制成,箱體規(guī)格為40 cm×40 cm ×40 cm,包括地箱和頂箱兩部分。地箱邊緣設有水槽,于每次觀測前,首先在水槽中加入適量的水,后將頂箱置于地箱上,上下箱體用水槽中的水密封。采用30 min罩箱時間,即每個采樣箱分別在罩箱后的 0、10、20和30 min抽取氣體樣品。采樣容器為100 mL帶三通閥的醫(yī)用注射器,將注射器與箱體一側的三通閥相連,抽取30~60 mL氣體樣品放入氣袋,同時記錄采樣時間。氣體樣品帶回實驗室后,一周之內用安捷特 7890A氣相色譜儀分析 CO2、CH4和N2O濃度。于7─8月期間選取晴好天氣的9:00─11:00時段采集氣體,采集氣體的同時測定土壤濕度、土壤表層、5、10、15、20和25 cm處溫度,利用鋁盒法測定表層土壤(0~20 cm)含水量。同時記錄氣壓和箱內溫度,其他氣象資料來自大興安嶺生態(tài)站的自動和人工氣象站。其中PD1-PD10為1號凍土丼到 10號凍土丼,LPL-Y、LPL-M1和LPL-N分別為杜香幼齡林、杜香中齡林和杜香近熟林,HL-M、HL-N和HL-O分別為草類中齡林、草類近熟林和草類過熟林。

    通量是指單位時間通過某單位面積界面輸送的物理量。氣體交換通量(F)計算公式:

    式中:ρ為箱內氣體密度;Δm和ΔC分別為Δt時間內箱內氣體質量和混合比濃度的變化;A、V、H分別為采樣箱的底面積、體積和氣室高度;ΔC/Δt為箱內氣體濃度變化。當F為負值時表示吸收,F(xiàn)為正值時表示排放。氣體通量計算過程中,通過公式中引入箱內溫度和氣壓值,對氣體濃度進行校正。

    1.4 土壤物理、化學性質測定

    SOC(土壤有機質)的測定:取5~15 cm土層土壤,采用重鉻酸鉀容量法,在外加熱的條件下(油浴溫度為180 ℃,沸騰5 min),用一定濃度的重鉻酸鉀-硫酸溶液氧化土壤有機質,剩余的重鉻酸鉀用硫酸亞鐵來滴定,從所消耗的重鉻酸鉀量,計算有機碳的質量分數(shù)。

    NT(土壤總氮)的測定:樣品在加速劑的參與下,用濃硫酸消煮時,各種含氮有機化合物,經(jīng)過復雜的高溫分解反應,轉化為氨與硫酸結合成硫酸銨。堿化后蒸餾出來的氨用硼酸吸收,以標準酸溶液滴定,求得土壤全氮的質量分數(shù)。

    堿解氮的測定:在擴散皿中,用 1.0 mol·L-1NaOH水解土壤,使易水解態(tài)氮堿解轉化為NH3,NH3擴散后為H3BO3所吸收。H3BO3吸收液中的再用標準算滴定,由此計算土壤中堿解氮的質量分數(shù)。

    土壤濕度:采用鋁盒法測定土壤含水量。用鋁盒取樣,并稱重,用烘箱在105 ℃下烘至恒重,計算樣品中損失的質量與烘干土質量的比例,即得土壤含水量。

    1.5 數(shù)據(jù)分析

    利用SAS 9.0統(tǒng)計軟件做相關性分析,利用多元回歸分析和方差分析測定CH4、CO2和N2O通量與土壤參數(shù)的相關性分析,以及CH4、CO2和N2O通量在不同經(jīng)營方式下的差異性分析。

    2 結果分析

    2.1 不同生境下土壤CH4通量

    在所選擇的所有樣地中(除沼澤濕地),總的CH4通量平均估計值為 C(-133.6±62.3) μg·m-2·h-1。14個樣地(除了PD2和PD7)CH4通量變化范圍為C(-242.0~-42.6) μg·m-2·h-1,變異系數(shù)為47%。在沼澤濕地 CH4通量平均值為 C(130.0±41.0) μg·m-2·h-1(表2),并且甲烷表現(xiàn)為釋放。變異系數(shù)范圍為6%~83%。由統(tǒng)計分析得出甲烷通量與土壤濕度、溫度以及其他土壤參數(shù)沒有顯著相關性。

    表2 興安落葉松林各樣地土壤-大氣溫室氣體通量Table 2 Flux of greenhouse gas between soil and atmosphere in Larix gmelinii Forest

    通過對數(shù)據(jù)的整合,在原始林、漸伐林和皆伐林3種不同森林經(jīng)營方式中,原始林CH4通量平均值為C(-99.3±60.0) μg·m-2·h-1,漸伐林CH4通量平均值為C(-159.1±25.1) μg·m-2·h-1,皆伐林CH4通量平均值為C(-170±55.0) μg·m-2·h-1(表4),由統(tǒng)計分析得出原始林、漸伐林和皆伐林之間沒有顯著性差異(P=0.1734)。

    表3 興安落葉松林樣地土壤性質Table 3 Soil properties of sample plots in Larix gmelinii Forest

    表4 興安落葉松林各樣地土壤-大氣溫室氣體通量Table 4 Effect of forest management on soil GHG flux in Larix gmelinii Forest

    在人工林、無林地和沼澤3種不同生境中,無林地CH4通量平均值為C(-44.4±31.9) μg·m-2·h-1,人工林CH4通量平均值為C(-91.1±5.9) μg·m-2·h-1,沼澤CH4通量平均值為C(130±41.0) μg·m-2·h-1(表4),沼澤中甲烷表現(xiàn)為釋放,而無林地和人工林表現(xiàn)為吸收,不同生境中的甲烷通量并未表現(xiàn)出顯著性差異(P=0.2004)。

    2.2 不同生境下土壤CO2通量

    在所有樣地里 CO2的通量范圍是C(94.9~1164.0) μg·m-2·h-1,平均值為C(671.3±324.3) μg·m-2·h-1,變異系數(shù)是 48%。變異系數(shù)范圍在5%~66%。在皆伐林6號凍土丼測到的CO2通量最高,達到 C(1164.0±416.0) μg·m-2·h-1,而在沼澤濕地中 2號凍土丼測到的 CO2通量最少,為C(94.4±44.8) μg·m-2·h-1(表2)。對于所有調查樣地中,CO2通量與土壤濕度沒有顯著相關性,僅表現(xiàn)為隨濕度的增加CO2通量而減少的趨勢(圖1)。隨溫度的升高,CO2通量呈增加趨勢(圖2)。CO2通量與土壤5 cm處溫度顯著相關(P=0.0423)。其他土壤參數(shù)與土壤CO2的釋放也并沒有發(fā)現(xiàn)顯著相關性。

    圖1 CO2通量與土壤濕度的關系Fig. 1 Correlation between soil moisture and CO2emissions

    圖2 CO2通量與土壤溫度的關系Fig. 2 Correlation between soil temperature and CO2emissions

    不同森林經(jīng)營方式下,原始林CO2通量平均值為C(535.0±245) μg·m-2·h-1,漸伐林CO2通量平均值為C(685.9±237.2) μg·m-2·h-1,皆伐林CO2通量平均值為C(828.3±219.2) μg·m-2·h-1,皆伐林CO2通量達到最高,而原始林CO2通量最低(表4),由統(tǒng)計分析得出原始林、漸伐林和皆伐林3種不同森林經(jīng)營方式下 CO2通量沒有顯著性差異(P=0.3675)。

    不同生境下的無林地、人工林和沼澤的CO2通量平均值分別為:C(474.4±89.3) μg·m-2·h-1、C(446.8±125.2) μg·m-2·h-1,C(114.3±27.5) μg·m-2·h-1(表 4),經(jīng)統(tǒng)計分析表明無林地、人工林和沼澤 3種不同生境下 CO2通量沒有顯著性差異(P=0.0791)。

    2.3 不同生境下土壤N2O通量

    各樣地土壤 N2O 通量的變化范圍是N(1.2~21.6) μg·m-2·h-1,平均值為 N(7.3±6.1) μg·m-2·h-1(表2),變異系數(shù)為84%,變異系數(shù)范圍為8%~133%。其中樣地PD5的變異性達到最高。經(jīng)統(tǒng)計分析表明,N2O通量與土壤溫度、土壤濕度沒有顯著相關性,N2O通量與土壤0~10、10~20 cm溫度及濕度呈正相關變化趨勢。N2O通量與SOC、TN也無顯著相關性,N2O通量與SOC呈正相關變化趨勢(圖5),與TN呈負相關變化趨勢(圖6)。

    在原始林、漸伐林和皆伐林3種不同經(jīng)營方式中,原始林N2O通量為N(5.0±3.8) μg·m-2·h-1,漸伐林N2O通量為N(4.2±2.1) μg·m-2·h-1,皆伐林N2O通量為N(8.7±6.5) μg·m-2·h-1。皆伐林N2O通量值最高。原始林和漸伐林N2O通量值相近(表4)。由統(tǒng)計分析得出,原始林、漸伐林和皆伐林3種不同經(jīng)營方式對N2O通量沒有顯著性差異(P=0.4569)。

    無林地 N2O 通量平均值為 N(18.7±20.7) μg·m-2·h-1,人工林 N2O 通量平均值為 N(3.3±3.0) μg·m-2·h-1,沼澤 N2O 通量平均值為 N(9.0±4.0) μg·m-2·h-1,顯然在無林地、人工林和沼澤3種不同生境下無林地的 N2O通量高于人工林和沼澤(表4)。統(tǒng)計分析結果表明,無林地、人工林和沼澤 3種不同生境對 N2O通量沒有顯著性差異(P=0.3418)。

    2.4 氣體通量與土壤性質、環(huán)境因子之間的相關分析

    本實驗的土壤為典型森林土壤,為酸性土壤,pH的范圍為4.94~6.01,SOC在33.04~99.15 g·kg-1,SOC質量分數(shù)最低的是PD4,而質量分數(shù)最高的樣地為LPL-N。TN的范圍則差異較大,在2.03~22.64 g·kg-1,其中HL-M質量分數(shù)最低,PD4的質量分數(shù)卻高達22.64 g·kg-1。由于SOC和TN范圍的差異,造成各樣地碳氮比范圍也較大,為 1.5~27.7。各樣地堿解氮變化范圍為:0.16~0.39,各個樣地間堿解氮的含量差異較小,其中LPLY樣地中堿解氮含量最高,而PD1樣地中堿解氮最低(表2)。在實驗期間由于降雨較頻繁,所以造成土壤濕度較大。氣體通量與土壤的SOC、TN、pH以及堿解氮之間通過SAS軟件分析并無顯著相關性,但是各氣體通量隨各土壤性質呈現(xiàn)一定的變化趨勢。CH4通量隨SOC的增加呈增加趨勢,隨TN的增加呈減少趨勢。CO2通量隨SOC增加呈現(xiàn)減少趨勢(圖3),隨TN的增加也呈減少趨勢(圖4)。N2O通量與SOC含量呈線性正相關(圖5),N2O通量與TN含量呈負相關(圖6)。各通量與pH和堿解氮并無任何相關關系。

    圖3 CO2通量與SOC的關系Fig. 3 Correlation between SOC and CO2emissions

    圖4 CO2通量與TN的關系Fig. 4 Correlation between TN and CO2emissions

    原始林、漸伐林和皆伐林3種不同森林經(jīng)營方式中,原始林的SOC為73.7 g·kg-1,高于漸伐林及皆伐林,而漸伐林的 SOC達到最低值,為 50.2g·kg-1。原始林、漸伐林和皆伐林TN值分別為:7.3、13.6、6.4 g·kg-1。所以原始林的碳氮比達到最大值,為 10.1,漸伐林的碳氮比最小,僅為 3.7,皆伐林的碳氮比與原始林相近,為 9.2。原始林、漸伐林和皆伐林3種不同經(jīng)營方式的pH值相近,分別為5.44、5.37、5.62。由于本次野外調查期間降雨頻繁,導致該區(qū)的濕度在各個樣地WFPS均很大,并且彼此間濕度差異很?。ū?)。由統(tǒng)計分析結果表明,在本實驗下的原始林、漸伐林和皆伐林3種不同經(jīng)營方式下土壤的SOC、TN、pH以及堿解氮沒有顯著性差異,不同經(jīng)營方式與CH4通量、CO2通量及N2O通量無顯著性差異。

    圖5 N2O通量與SOC的關系Fig. 5 Correlation between SOC and N2O emissions

    圖6 N2O通量與TN的關系Fig. 6 Correlation between TN and N2O emissions

    3 討論

    3.1 土壤性質對土壤溫室氣體的影響

    很多研究表明土壤濕度和溫度對CO2通量有重要作用。在我們的野外調查中,CO2與溫度呈正相關,隨溫度的增加 CO2通量也呈增加趨勢。但是CO2通量僅與 5 cm深度土壤溫度存在顯著相關性(P=0.0423),與王飛(2013)研究的有關對杜香-興安落葉松林土壤呼吸速率與土壤溫度之間呈極顯著的相關關系,而且5 cm土壤溫度與土壤呼吸之間的相關性要好于10 cm土深的土壤溫度,也就是說土壤深度5 cm的地溫能夠更好反映溫度對土壤呼吸的影響這一結果一致。賈淑霞(2009)對興安落葉松林的研究也得出一致的結果。在其余土壤層次的溫度沒有顯著相關性。還有許多的實驗結果表明溫度與CO2排放有顯著的指數(shù)函數(shù)相關關系,其中10 cm和5 cm土層溫度對土壤CO2排放影響最大。馬秀枝(2012)等的研究也得出CO2通量與土壤溫度呈顯著正相關關系。

    地表 CO2通量與土壤濕度之間呈現(xiàn)線性負相關,但是相關性不顯著,解釋率較低,該結論與張磊(2013)的研究結果一致。陳全勝等(2004)對高寒草地表層土壤(0~20 cm)水分研究結果也表明水分狀況對土壤呼吸有一定程度的影響,從而影響 CO2通量。同時,在 Inclan(2007)與 Dilustro(2005)的研究中分別顯示,土壤含水率對土壤地表 CO2通量的影響存在一個臨界值(12% cm3H20cm-3soil),當土壤濕度低于這個臨界值時,地表CO2通量與土壤濕度成正的線性關系,高于臨界值時,則影響不明顯。由于在該研究時間段內,該地區(qū)降雨頻繁導致土壤濕度偏大,可能導致土壤濕度高于了某一臨界值而導致與CO2通量的關系成負相關且沒有顯著相關性。CO2通量顯示與SOC、TN呈負相關關系。謝軍飛和李玉娥(2002)研究表明pH值存在最適范圍,超出這個范圍時,會導致溫室氣體的排放大幅度減少,而本研究各樣地pH值相近,所以并沒有發(fā)現(xiàn)pH值對溫室氣體通量造成影響。

    在整個研究區(qū)內各個樣地中,PD2 和 PD7為沼澤濕地,僅有這兩個樣地是CH4的源,其他樣地則均為CH4的匯。與陳全勝(2004)研究得到的土壤含水量的增大,CH4由吸收轉為排放研究結果一致。在研究中,CH4通量沒有顯示出與土壤溫度和濕度有任何相關性。而徐星凱和周禮愷(1999)研究表明,土壤水分含量明顯影響CH4通量,并且隨土壤水分增加,CH4通量呈增加趨勢。本研究并沒有得出水分與CH4通量的相關關系可能是由于土壤濕度超出臨界值,對CH4通量的影響不顯著。CH4通量隨SOC的增加呈增加趨勢,隨TN的增加呈減少趨勢。

    通過對14個樣地(除PD2和PD7外),N2O通量隨0~10、10~20 cm的土壤濕度的增加呈增加趨勢,但是卻沒有顯示顯著的相關性。這一研究結果與Christensen(2008)等研究的N2O通量隨土壤水分的增加而增加,而高水分含量條件下N2O的產生并非與土壤含水量成正比這一結果一致。N2O通量在土壤5 cm和10 cm處溫度顯示隨溫度升高通量增加的趨勢,同樣地,沒有顯著相關性。該研究結果劉碩等(2013)研究得出的針葉林和針闊混交林土壤 N2O的排放通量隨溫度和土壤含水量升高而增加結果一致。N2O通量與TN質量分數(shù)呈負相關,與Zhisheng Yao等(2010)研究的結果相反。而N2O通量與SOC質量分數(shù)呈線性正相關。而各個樣地的CO2通量與N2O通量變化趨勢一致。這一研究結果與Zhisheng Yao等(2010)研究結果一致。

    4 結論

    (1)CH4通量在7─8月間與土壤溫度、土壤濕度以及土壤性質均未發(fā)現(xiàn)有顯著相關關系,僅表現(xiàn)為隨CH4通量隨SOC的增加呈增加趨勢,隨TN的增加呈減少趨勢。

    (2)CO2通量在7─8月間與土壤濕度沒有顯著相關性,僅表現(xiàn)為隨濕度的增加CO2通量而減少的趨勢。隨溫度的升高,CO2通量呈增加趨勢。CO2通量與土壤5 cm處溫度顯著相關(P=0.0423)。其他土壤參數(shù)與土壤CO2的釋放也并沒有發(fā)現(xiàn)顯著相關性,僅表現(xiàn)為隨SOC的增加而減少,隨TN的增加而減少。

    (3)N2O通量在7─8月間隨0~10、10~20 cm的土壤濕度的增加呈增加趨勢,但是卻沒有顯示顯著的相關性。N2O通量在土壤5 cm和10 cm處溫度顯示隨溫度升高通量增加的趨勢,同樣地,沒有顯著相關性。N2O通量與TN呈負相關,而N2O通量與SOC呈線性正相關。而各個樣地的CO2通量與N2O通量變化趨勢一致。

    (4)不同森林經(jīng)營方式及不同生境在 7─8月間對CH4通量、CO2通量以及N2O通量沒有顯著性差異,不同森林經(jīng)營方式及不同生境,土壤性質即SOC、TN、堿解氮以及pH與CH4通量、CO2通量以及N2O通量沒有顯著相關性。

    CHAPIUS-LARDY L, WRAGE N, METAY A, et al. 2007. Soils, a sink for N2O[J]. Global Change Biol, 13(1): 1-17.

    CHRISTENSEN T R, FRIBORG T, JOHANSSON M. 2008. Trace gas budgets of high arctic permafrost regions [C] Proceedings, NICOP: 251-256.

    DILUSTRO J J, COLLINS B, DUNCAN L, et al. 2005. Moisture and soil texture effects on Soil CO2efflux components in southeastern mixed pine forests[J]. Forest Ecology and Management, 204(1): 87-97.

    INCLAN R, DELA TORRE D, BENITO M, et al. 2007. Soil CO2efflux in a mixed pine-oak forest in Valsain (Central Spain)[J]. Scientific World Journal, 7(s1): 166-174.

    LILLY A, BALL B C, MCTAGGART I P, et al. 2009. Spatial modeling of nitrous oxide emissions at the national scale using soil, climate and land use information[J]. Global Change Biol, 15(9): 2321-2332.

    LI Y L, OTIENO D, OWEN K, et al. 2008. Temporal variability in soil CO2emission in an orchard forest ecosystem[J]. Pedosphere, 18(3): 273-283.

    SMITH K A, BALL T, CONEN F, et al. 2003. Exchange of greenhouse gases between soil and atmosphere: interactions of soil physical factors and biological processes[J]. European Journal of Soil Science, 54(4): 779-791.

    SUBKE J A, REICHSTEIN M, TENHUNEN J D. 2003. Explaining temporal variation in soil CO2efflux in a mature spruce forest in Southern Germany [J]. Soil Biology & Biochemistry, 35(11): 1467-1483.

    RAICH J W, SCHLESINGER W H. 1992. The Global carbon dioxide flux in soil respiration and its relationship to vegetation [J]. Tellus, 44B(2): 81-99.

    ZHISHENG YAO, BENJAMIN WOLF, WEIWEI CHEN, et al. 2010. Spatial variability of N2O, CH4and CO2fluxes within the Xilin River catchment of Inner Mongolia, China: a soil core study[J]. Plant and Soil, 331(1-2): 341-359.

    鮑士旦. 2000. 土壤農化分析[M]. 中國農業(yè)出版社.

    陳全勝, 李凌浩, 韓興國, 等. 2004. 典型溫帶草原群落土壤呼吸溫度敏感性與土壤水分的關系[J]. 生態(tài)學報, 24(4): 831-836.

    傅民杰, 王傳寬, 王穎, 等. 2009. 氣候暖化對解凍期不同緯度興安落葉松林土壤氧化亞氮釋放的影響[J]. 應用生態(tài)學報, 20(7): 1635-1642.

    蔣延玲, 周廣勝. 1999. 中國主要森林生態(tài)系統(tǒng)公益的評估[J]. 植物生態(tài)學報, 23(5): 426-432.

    賈淑霞. 2009. 落葉松和水曲柳人工林土壤、根系和土壤微生物呼吸研究[D]. 哈爾濱: 東北林業(yè)大學: 48-49.

    劉碩, 李玉娥, 萬運帆, 等. 2013. 溫度和土壤含水量對溫帶森林土壤溫室氣體排放的影響[J]. 生態(tài)環(huán)境學報, 22(7): 1093-1098.

    莫江明, 方運霆, 李德軍, 等. 2006. 鼎湖山主要森林土壤 CO2排放和CH4吸收特征[J]. 廣西植物, 26(2): 142-147.

    馬秀枝, 張秋良, 李長生, 等. 2012. 寒溫帶興安落葉松林土壤溫室氣體通量的時間變異[J]. 應用生態(tài)學報, 23(8): 2149-2156.

    齊玉春, 羅輯, 董云社, 等. 2002. 貢嘎山山地暗針葉林帶森林土壤溫室氣體N2O和CH4排放研究[J]. 中國科技(D輯), 32(11): 934-941.

    孫向陽, 徐化成. 2001. 北京低山區(qū)兩種人工林土壤中N2O排放通量的研究[J]. 林業(yè)科學, 37(5): 57-63.

    沙麗清, 鄭征, 唐建維, 等. 2004. 西雙版納熱帶季節(jié)雨林的土壤呼吸研究[J]. 中國科學(D輯), 34(2): 167-174.

    王飛. 2013. 興安落葉松天然林碳密度與碳平衡研究[D]. 內蒙古呼和浩特: 內蒙古農業(yè)大學: 103-104.

    謝軍飛, 李玉娥. 2002. 農田土壤溫室氣體排放機理與影響因素研究進展[J]. 中國農業(yè)氣象, 23(4): 47-52.

    肖冬梅, 王淼, 姬蘭柱, 等. 2009. 長白山闊葉紅松林土壤N2O排放通量的變化特征[J]. 生態(tài)學雜志, 23(5): 46-52.

    徐星凱, 周禮愷. 1999. 土壤源CH4氧化的主要影響因子與減排措施[J].生態(tài)農業(yè)研究, 7: 18-22.

    張磊. 2013. 亞熱帶杉木林土壤溫室氣體排放對氮磷添加的響應[D]. 西南大學: 22-23.

    Effect of N2O, CH4and CO2Fluxes in the Thriving Season of Larix gmelinii Forest in Cold-temperate Zone under Different Forest Management and Topographic Condition

    REN Le1, MA XiuZhi2*, FAN Xuesong1
    1. College of Ecology and Environmental Science, Inner Mongolia Agricultural University, Huhhot 010019, China ; 2. College of Forest, Inner Mongolia Agricultural University, Huhhot 010019, China

    In order to identify the effects of different forest management (virgin forest, clear-cutting forest and successive cutting forest) and topographic conditions (forestland, non-forest land, marshland) on soil-atmosphere exchange of greenhouse gases, nitrous oxide (N2O), methane (CH4) and carbon dioxide (CO2) fluxes of soil cores from 16 representative sites were determined in situ of Larix gmelinii Forest of Inner Mongolia, and the relative soil properties such as soil moisture, soil temperature in each soil layer was measured simultaneously. The soil average CH4fluxes in the thriving growing season was (-133.6±62.3) μg·m-2·h-1. Strong CH4emission was found in two wetlands (PD2 and PD7), whereas CH4showed soil sink in other 14 sample sites, with a range from -242.0 to -42.6 μg·m-2·h-1and CV of 47%.CO2emissions from all sites ranged from 94.9 to 1164.0 μg C·m-2·h-1, with a mean value of (671.3±324.3) μg C·m-2·h-1and a CV of 48%. Soil N2O emissions across all of the investigated sites ranged from 1.2 to 21.6 μg N·m-2·h-1, with a mean of (7.3±6.1) μg N·m-2·h-1and a CV of 84%. Different forest management and different habitats have no obvious effects on soil N2O, CH4and CO2fluxes. There is no correlation was found between soil CH4flux and soil temperature and soil moisture at each depths, but CH4flux is positively correlated with soc, but negatively with soil TN. CO2flux is negatively related with soil moisture in each layer, and positively correlated with soil temperature at 5 cm depth(P=0.042 3). Positive correlation was found between N2O flux and soil moisture, temperature, and SOC, but negatively with soil TN.

    Larix gmelinii forest; CO2; CH4; N2O; different forest management; different habitat

    10.16258/j.cnki.1674-5906.2015.03.003

    Q948

    A

    1674-5906(2015)03-0378-09

    任樂,馬秀枝,范雪松. 不同經(jīng)營方式及生境對大興安嶺高緯度林區(qū)生長盛季森林土壤CO2、CH4、N2O通量的影響[J]. 生態(tài)環(huán)境學報, 2015, 24(3): 378-386.

    REN Le, MA XiuZhi, FAN Xuesong. Effect of N2O, CH4and CO2Fluxes in the Thriving Season of Larix gmelinii Forest in Cold-temperate Zone under Different Forest Management and Topographic Condition [J]. Ecology and Environmental Sciences, 2015, 24(3): 378-386.

    國家自然科學基金項目(31160117;31260119)

    任樂(1990年生),女,碩士研究生,專業(yè)方向為森林生態(tài)學。E-mail: 623349271@qq.com *通信作者:馬秀枝(1974年生),女,副教授,主要從事生態(tài)系統(tǒng)碳循環(huán)以及溫室氣體相關研究。E-mail: luckmxy@imau.edu.cn

    2014-11-26

    猜你喜歡
    興安土壤濕度落葉松
    山西落葉松雜交良種逾10萬畝
    土壤濕度傳感器在園林綠化灌溉上的應用初探
    興安行吟(組詩)
    北極光(2019年12期)2020-01-18 06:22:16
    興安四月樹
    文苑(2019年22期)2019-12-07 05:29:32
    基于51單片機控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    長白落葉松離體再生體系的建立
    追夢興安
    草原歌聲(2018年2期)2018-12-03 08:14:12
    神奇的落葉松提取物
    四川盆地土壤濕度時空分布及影響因子分析
    中國不同氣候區(qū)土壤濕度特征及其氣候響應
    高清在线视频一区二区三区| 97人妻天天添夜夜摸| 久久久久久久国产电影| 亚洲久久久国产精品| 男的添女的下面高潮视频| 免费播放大片免费观看视频在线观看| 日韩成人av中文字幕在线观看| 久久久久精品性色| 精品一品国产午夜福利视频| 美女大奶头黄色视频| 成人国语在线视频| 王馨瑶露胸无遮挡在线观看| 久久久精品区二区三区| 亚洲精华国产精华液的使用体验| 汤姆久久久久久久影院中文字幕| av天堂久久9| 人妻少妇偷人精品九色| 在线观看免费视频网站a站| 狠狠婷婷综合久久久久久88av| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 亚洲欧美一区二区三区国产| 亚洲人成网站在线观看播放| 18在线观看网站| 免费久久久久久久精品成人欧美视频 | videosex国产| 亚洲av成人精品一二三区| 亚洲精品久久午夜乱码| 国产精品蜜桃在线观看| 另类亚洲欧美激情| 免费高清在线观看视频在线观看| 天美传媒精品一区二区| 欧美日韩亚洲高清精品| 十八禁网站网址无遮挡| 国产国语露脸激情在线看| 日本av手机在线免费观看| 国产精品一二三区在线看| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 免费观看性生交大片5| 久久久精品区二区三区| 久久av网站| 久久免费观看电影| 两个人看的免费小视频| 性色av一级| 亚洲欧美成人精品一区二区| 国产一区二区三区av在线| 国产日韩一区二区三区精品不卡| 女的被弄到高潮叫床怎么办| 亚洲久久久国产精品| 国产精品一国产av| 天堂8中文在线网| 99热这里只有是精品在线观看| 国产成人午夜福利电影在线观看| 街头女战士在线观看网站| 久久午夜综合久久蜜桃| 妹子高潮喷水视频| 日韩成人伦理影院| 午夜av观看不卡| 青青草视频在线视频观看| 人妻系列 视频| 热99久久久久精品小说推荐| 欧美日韩视频精品一区| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品古装| 亚洲,欧美精品.| 国产极品天堂在线| 久久青草综合色| 久热久热在线精品观看| 国产免费现黄频在线看| 国产男人的电影天堂91| 午夜激情av网站| 日韩精品免费视频一区二区三区 | 91精品伊人久久大香线蕉| 黑人巨大精品欧美一区二区蜜桃 | 曰老女人黄片| 国产又色又爽无遮挡免| 久久99蜜桃精品久久| 欧美性感艳星| 免费黄网站久久成人精品| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 王馨瑶露胸无遮挡在线观看| h视频一区二区三区| 精品人妻偷拍中文字幕| av天堂久久9| 热re99久久精品国产66热6| 久久久久久久久久成人| 日本wwww免费看| videossex国产| tube8黄色片| 久久韩国三级中文字幕| 国产一区有黄有色的免费视频| 永久网站在线| 久久ye,这里只有精品| 国产乱来视频区| 伦精品一区二区三区| 成人黄色视频免费在线看| 卡戴珊不雅视频在线播放| 日韩制服丝袜自拍偷拍| 最近最新中文字幕大全免费视频 | 搡女人真爽免费视频火全软件| av一本久久久久| 夜夜骑夜夜射夜夜干| 一级爰片在线观看| 日韩av在线免费看完整版不卡| 亚洲精品国产色婷婷电影| 国产日韩欧美视频二区| 亚洲成色77777| 99热网站在线观看| 国产成人aa在线观看| kizo精华| 国产成人91sexporn| 夜夜骑夜夜射夜夜干| 在线观看免费视频网站a站| 蜜臀久久99精品久久宅男| 视频中文字幕在线观看| 日韩制服丝袜自拍偷拍| 日韩av不卡免费在线播放| 久久影院123| 中文字幕制服av| 亚洲精品美女久久av网站| 免费播放大片免费观看视频在线观看| 在线 av 中文字幕| 亚洲精品美女久久久久99蜜臀 | 夜夜爽夜夜爽视频| 国产精品成人在线| 性色avwww在线观看| 在线观看人妻少妇| 日本av免费视频播放| 免费久久久久久久精品成人欧美视频 | 成年人午夜在线观看视频| 国产无遮挡羞羞视频在线观看| 免费大片黄手机在线观看| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 夜夜爽夜夜爽视频| 日韩中字成人| 中文乱码字字幕精品一区二区三区| 亚洲国产看品久久| 亚洲色图 男人天堂 中文字幕 | 亚洲久久久国产精品| 久久这里只有精品19| 亚洲av免费高清在线观看| 两个人免费观看高清视频| 成人黄色视频免费在线看| 亚洲国产最新在线播放| a级片在线免费高清观看视频| 哪个播放器可以免费观看大片| 国产免费一区二区三区四区乱码| 久热这里只有精品99| 日本与韩国留学比较| 久久韩国三级中文字幕| 久久这里有精品视频免费| 两性夫妻黄色片 | 欧美日本中文国产一区发布| 日本wwww免费看| 五月天丁香电影| 国产精品人妻久久久久久| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 不卡视频在线观看欧美| 777米奇影视久久| 日韩一区二区三区影片| 国产精品久久久久久久电影| 黑人欧美特级aaaaaa片| 国产免费又黄又爽又色| 在线观看免费日韩欧美大片| 国产在视频线精品| 18禁国产床啪视频网站| av又黄又爽大尺度在线免费看| 韩国av在线不卡| 国产在线免费精品| 久久ye,这里只有精品| 熟妇人妻不卡中文字幕| 日韩人妻精品一区2区三区| 日韩成人av中文字幕在线观看| videos熟女内射| 日韩一本色道免费dvd| xxxhd国产人妻xxx| h视频一区二区三区| 美女主播在线视频| 在线免费观看不下载黄p国产| 日本黄色日本黄色录像| 国产男女内射视频| 亚洲图色成人| 免费人妻精品一区二区三区视频| 久久午夜福利片| 亚洲精品久久成人aⅴ小说| 国产精品无大码| 99九九在线精品视频| 秋霞伦理黄片| 中文字幕人妻熟女乱码| 狂野欧美激情性bbbbbb| 国产av精品麻豆| 免费av中文字幕在线| 久久久久久久大尺度免费视频| 日本免费在线观看一区| 看免费av毛片| 宅男免费午夜| 在线观看三级黄色| 亚洲精品美女久久av网站| 亚洲国产av新网站| 9热在线视频观看99| 少妇高潮的动态图| 一级a做视频免费观看| 欧美成人精品欧美一级黄| 久久人妻熟女aⅴ| 亚洲情色 制服丝袜| 91成人精品电影| 91aial.com中文字幕在线观看| 丝袜美足系列| 热99久久久久精品小说推荐| 精品久久久久久电影网| 国产69精品久久久久777片| 国产av精品麻豆| 精品第一国产精品| 久久国产精品男人的天堂亚洲 | 国产日韩欧美亚洲二区| 成人黄色视频免费在线看| 在线 av 中文字幕| 精品国产一区二区三区久久久樱花| 男女高潮啪啪啪动态图| 精品亚洲成a人片在线观看| 蜜臀久久99精品久久宅男| 亚洲精品国产av成人精品| 欧美老熟妇乱子伦牲交| 日韩av免费高清视频| 欧美 日韩 精品 国产| 久久久久久久久久久久大奶| 美国免费a级毛片| 宅男免费午夜| 国产精品欧美亚洲77777| 成年av动漫网址| 美女大奶头黄色视频| 亚洲成人手机| 高清视频免费观看一区二区| 一级a做视频免费观看| 乱人伦中国视频| 这个男人来自地球电影免费观看 | 日本av手机在线免费观看| 美女中出高潮动态图| 婷婷色综合大香蕉| 欧美3d第一页| 欧美精品人与动牲交sv欧美| 午夜老司机福利剧场| 深夜精品福利| 男女高潮啪啪啪动态图| 免费人成在线观看视频色| 99热国产这里只有精品6| 精品国产一区二区三区四区第35| 日韩伦理黄色片| 欧美日韩av久久| 成人亚洲欧美一区二区av| 青青草视频在线视频观看| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 妹子高潮喷水视频| 午夜日本视频在线| 97人妻天天添夜夜摸| 欧美日韩综合久久久久久| 国产综合精华液| 欧美97在线视频| 国产一区二区在线观看日韩| 一级黄片播放器| 国产精品久久久久久久电影| 亚洲,一卡二卡三卡| 赤兔流量卡办理| 一级毛片 在线播放| 晚上一个人看的免费电影| 精品福利永久在线观看| 午夜老司机福利剧场| 捣出白浆h1v1| 色视频在线一区二区三区| 日韩一区二区视频免费看| av卡一久久| kizo精华| 欧美精品av麻豆av| 亚洲丝袜综合中文字幕| 国产亚洲欧美精品永久| 亚洲美女搞黄在线观看| 亚洲熟女精品中文字幕| 一边摸一边做爽爽视频免费| 日韩av不卡免费在线播放| 久久精品国产综合久久久 | 精品午夜福利在线看| 欧美亚洲日本最大视频资源| 亚洲人成77777在线视频| 久久人人爽人人爽人人片va| 97在线人人人人妻| 精品国产乱码久久久久久小说| 国产乱来视频区| 欧美精品人与动牲交sv欧美| 丝袜脚勾引网站| 九九在线视频观看精品| 啦啦啦视频在线资源免费观看| 日韩成人伦理影院| 在线观看国产h片| 天天影视国产精品| 一本色道久久久久久精品综合| 如何舔出高潮| 黄片无遮挡物在线观看| 国产精品99久久99久久久不卡 | av有码第一页| 黄色毛片三级朝国网站| 国产一区二区三区综合在线观看 | 亚洲精品美女久久av网站| 在线观看国产h片| 午夜免费观看性视频| 大码成人一级视频| 丁香六月天网| 日韩大片免费观看网站| 亚洲国产精品国产精品| 一区二区av电影网| 亚洲精品中文字幕在线视频| 成人手机av| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 婷婷色综合www| 欧美日韩综合久久久久久| 国产免费视频播放在线视频| 性高湖久久久久久久久免费观看| 男男h啪啪无遮挡| 久久久久精品人妻al黑| 校园人妻丝袜中文字幕| 99国产精品免费福利视频| 精品一区在线观看国产| av视频免费观看在线观看| 精品久久蜜臀av无| 乱码一卡2卡4卡精品| 日韩精品有码人妻一区| 飞空精品影院首页| av免费观看日本| 热re99久久精品国产66热6| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 欧美精品一区二区大全| 国产精品99久久99久久久不卡 | 日本免费在线观看一区| 人妻人人澡人人爽人人| 精品熟女少妇av免费看| 在线 av 中文字幕| 国产永久视频网站| 亚洲欧美日韩另类电影网站| 亚洲av中文av极速乱| 下体分泌物呈黄色| 亚洲一区二区三区欧美精品| a级毛色黄片| 91午夜精品亚洲一区二区三区| 五月天丁香电影| 亚洲欧美一区二区三区黑人 | 午夜影院在线不卡| 2018国产大陆天天弄谢| 男人添女人高潮全过程视频| 国产国拍精品亚洲av在线观看| 国产精品久久久av美女十八| 亚洲av福利一区| 一个人免费看片子| 大码成人一级视频| 日韩欧美一区视频在线观看| 黑人猛操日本美女一级片| 久久久精品免费免费高清| 国产亚洲午夜精品一区二区久久| 80岁老熟妇乱子伦牲交| 91久久精品国产一区二区三区| 日韩精品免费视频一区二区三区 | 自线自在国产av| 日韩中字成人| 2021少妇久久久久久久久久久| 在线天堂最新版资源| 成人无遮挡网站| 99国产综合亚洲精品| 国产免费视频播放在线视频| 黑丝袜美女国产一区| 午夜老司机福利剧场| 精品国产一区二区三区久久久樱花| av国产久精品久网站免费入址| 黑丝袜美女国产一区| 一区二区三区乱码不卡18| 国产熟女欧美一区二区| www日本在线高清视频| 成年美女黄网站色视频大全免费| 人妻少妇偷人精品九色| 超色免费av| 性色av一级| 国产 一区精品| 女的被弄到高潮叫床怎么办| 日韩制服骚丝袜av| 精品亚洲成国产av| 18禁裸乳无遮挡动漫免费视频| 少妇精品久久久久久久| 五月玫瑰六月丁香| 亚洲丝袜综合中文字幕| 国产日韩一区二区三区精品不卡| 久久鲁丝午夜福利片| 久久青草综合色| av有码第一页| 春色校园在线视频观看| 久久久久久久久久久久大奶| a 毛片基地| 日韩人妻精品一区2区三区| 少妇 在线观看| 国产极品天堂在线| 欧美日韩亚洲高清精品| 久久这里只有精品19| 精品亚洲成a人片在线观看| 日韩中文字幕视频在线看片| h视频一区二区三区| 亚洲精品国产色婷婷电影| 亚洲国产精品999| 人人澡人人妻人| 九草在线视频观看| 日韩免费高清中文字幕av| 久久99一区二区三区| 精品一区二区三区视频在线| 欧美日韩成人在线一区二区| 在线亚洲精品国产二区图片欧美| 日日撸夜夜添| 日韩一区二区视频免费看| 搡老乐熟女国产| 亚洲国产精品999| 亚洲色图综合在线观看| 老司机影院毛片| 国内精品宾馆在线| freevideosex欧美| 在线亚洲精品国产二区图片欧美| 一级毛片 在线播放| 国产又爽黄色视频| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 久久青草综合色| 免费高清在线观看日韩| a级毛片黄视频| 男女边摸边吃奶| 久久精品夜色国产| 免费观看a级毛片全部| 精品99又大又爽又粗少妇毛片| 精品一区二区三区四区五区乱码 | 国产精品秋霞免费鲁丝片| 一个人免费看片子| 男女国产视频网站| 看免费成人av毛片| 一级毛片我不卡| 交换朋友夫妻互换小说| 妹子高潮喷水视频| 色婷婷av一区二区三区视频| 成人综合一区亚洲| 久久综合国产亚洲精品| 成人亚洲欧美一区二区av| 在线 av 中文字幕| 欧美精品一区二区免费开放| 伊人久久国产一区二区| 国产激情久久老熟女| 亚洲精品一二三| 精品少妇久久久久久888优播| 亚洲人成77777在线视频| 老司机亚洲免费影院| 亚洲精品视频女| 欧美成人精品欧美一级黄| 七月丁香在线播放| 精品国产一区二区三区四区第35| 亚洲情色 制服丝袜| 国产免费现黄频在线看| 亚洲国产精品成人久久小说| 国产精品一国产av| 亚洲伊人色综图| 久久久久视频综合| 热99久久久久精品小说推荐| 久久国产亚洲av麻豆专区| 成人无遮挡网站| 久久女婷五月综合色啪小说| 春色校园在线视频观看| 亚洲精品久久久久久婷婷小说| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 日韩视频在线欧美| 色视频在线一区二区三区| 国产高清不卡午夜福利| 51国产日韩欧美| 欧美日韩视频精品一区| 国产色爽女视频免费观看| 男女国产视频网站| 国产白丝娇喘喷水9色精品| 色婷婷久久久亚洲欧美| 制服诱惑二区| 十八禁高潮呻吟视频| 高清av免费在线| 久久久久久久久久成人| 亚洲欧美清纯卡通| 只有这里有精品99| av在线app专区| 一级毛片电影观看| 大香蕉久久网| 下体分泌物呈黄色| 欧美bdsm另类| 欧美精品高潮呻吟av久久| 久久午夜综合久久蜜桃| 国产高清不卡午夜福利| 侵犯人妻中文字幕一二三四区| 欧美日韩视频精品一区| 亚洲国产毛片av蜜桃av| 国产综合精华液| 精品国产露脸久久av麻豆| 欧美精品亚洲一区二区| 五月伊人婷婷丁香| 亚洲国产看品久久| 亚洲色图综合在线观看| 午夜福利在线观看免费完整高清在| 99热这里只有是精品在线观看| 成人国产麻豆网| h视频一区二区三区| 亚洲,一卡二卡三卡| 男女免费视频国产| 超碰97精品在线观看| 18禁裸乳无遮挡动漫免费视频| kizo精华| 精品久久久精品久久久| 久久99精品国语久久久| 午夜福利视频在线观看免费| 国产精品女同一区二区软件| 免费女性裸体啪啪无遮挡网站| 秋霞在线观看毛片| 国产成人免费观看mmmm| 蜜桃在线观看..| 毛片一级片免费看久久久久| 精品人妻偷拍中文字幕| 韩国精品一区二区三区 | 成人免费观看视频高清| 亚洲欧美一区二区三区国产| 国产欧美日韩综合在线一区二区| 国产色婷婷99| 黄色怎么调成土黄色| 欧美xxⅹ黑人| 晚上一个人看的免费电影| 大片电影免费在线观看免费| 精品人妻熟女毛片av久久网站| 高清欧美精品videossex| 欧美精品一区二区免费开放| 午夜影院在线不卡| 亚洲色图综合在线观看| 美女大奶头黄色视频| 免费黄网站久久成人精品| 日本猛色少妇xxxxx猛交久久| 五月开心婷婷网| 下体分泌物呈黄色| 久久这里只有精品19| 免费看不卡的av| 精品国产国语对白av| 91午夜精品亚洲一区二区三区| 丝瓜视频免费看黄片| 欧美精品av麻豆av| 国产精品久久久久久久电影| 国产精品久久久久久精品古装| 一本—道久久a久久精品蜜桃钙片| 97超碰精品成人国产| videossex国产| 亚洲精品,欧美精品| 久久久久久久亚洲中文字幕| 一区二区三区四区激情视频| 边亲边吃奶的免费视频| 日韩制服丝袜自拍偷拍| 久久午夜综合久久蜜桃| 亚洲第一区二区三区不卡| 飞空精品影院首页| 夫妻午夜视频| 美女脱内裤让男人舔精品视频| 在线观看www视频免费| 777米奇影视久久| 国产一区二区激情短视频 | 中国美白少妇内射xxxbb| 午夜老司机福利剧场| 亚洲美女黄色视频免费看| 久久久久久久亚洲中文字幕| 国精品久久久久久国模美| 男女边摸边吃奶| 在线观看人妻少妇| 啦啦啦中文免费视频观看日本| 欧美精品人与动牲交sv欧美| 十八禁高潮呻吟视频| 国产成人av激情在线播放| 午夜福利在线观看免费完整高清在| 女的被弄到高潮叫床怎么办| 激情五月婷婷亚洲| 你懂的网址亚洲精品在线观看| 色视频在线一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲美女搞黄在线观看| 久久热在线av| 日本黄大片高清| 高清毛片免费看| 9色porny在线观看| 人成视频在线观看免费观看| 欧美日韩亚洲高清精品| 伦理电影免费视频| 永久网站在线| 午夜精品国产一区二区电影| 亚洲欧洲精品一区二区精品久久久 | 好男人视频免费观看在线| 国产精品成人在线| 日韩精品免费视频一区二区三区 | 久久鲁丝午夜福利片| 51国产日韩欧美| 日韩三级伦理在线观看| 国产精品.久久久| 一边摸一边做爽爽视频免费| 777米奇影视久久| 国产精品一区二区在线观看99| 免费看av在线观看网站| 久久久国产一区二区| 久久国产精品大桥未久av| 国产成人欧美| 97在线视频观看| 午夜福利在线观看免费完整高清在| 精品99又大又爽又粗少妇毛片| 丰满迷人的少妇在线观看|