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

    祁連山區(qū)主要植被類(lèi)型下土壤團(tuán)聚體變化特征

    2020-11-30 14:33:18于文竹趙恒策
    關(guān)鍵詞:土壤結(jié)構(gòu)祁連山草甸

    魏 霞,賀 燕,魏 寧,于文竹,崔 霞,趙恒策

    祁連山區(qū)主要植被類(lèi)型下土壤團(tuán)聚體變化特征

    魏 霞1,賀 燕1,魏 寧2,于文竹1,崔 霞1,趙恒策1

    (1. 蘭州大學(xué)資源環(huán)境學(xué)院,蘭州 730000;2. 西北農(nóng)林科技大學(xué)理學(xué)院,楊陵 712100)

    土壤團(tuán)聚體是反映土壤結(jié)構(gòu)穩(wěn)定性、肥力和質(zhì)量狀況的重要指標(biāo),與侵蝕過(guò)程、水土流失、環(huán)境質(zhì)量密切相關(guān),研究土壤的理化性質(zhì),對(duì)保護(hù)土壤資源、提高生產(chǎn)率、維護(hù)土壤生態(tài)系統(tǒng)平衡具有重要意義。以青藏高原北緣的祁連山區(qū)為研究對(duì)象,采集該區(qū)4種主要植被類(lèi)型——荒漠、草原、草甸及灌叢的土壤,分析了不同植被類(lèi)型的土壤團(tuán)聚體指標(biāo),水穩(wěn)定性團(tuán)聚體質(zhì)量百分?jǐn)?shù)(percentage of water-stable aggregates,WSA)、平均質(zhì)量直徑(mean weightdiameter,MWD)、幾何平均直徑(geometric mean diameter,GMD)、結(jié)構(gòu)體破壞率(aggregate destruction rate,PAD)、平均重量比表面積(mean weight soil specific area,MWSSA)和分形維數(shù)隨海拔高度和土壤深度的變化趨勢(shì)。結(jié)果表明:0~30 cm土層WSA、MWD、GMD、MWSSA 依次:荒漠<草原<草甸<灌叢,>30~40 cm土層植被類(lèi)型對(duì)團(tuán)聚體無(wú)顯著影響(>0.05);隨土壤深度的增加,土壤團(tuán)聚體穩(wěn)定性和團(tuán)聚度逐漸降低,土壤結(jié)構(gòu)趨于惡化,草甸帶表層土壤團(tuán)聚體穩(wěn)定性顯著高于深層土壤(0.05);隨海拔高度增大,在海拔1 692~2 800 m土壤團(tuán)聚體逐漸穩(wěn)定,土壤結(jié)構(gòu)改善,在海拔2 800~3 639 m土壤團(tuán)聚體穩(wěn)定性逐漸降低,土壤結(jié)構(gòu)趨于惡化。WSA、MWD和GMD受1~4 mm粒級(jí)主導(dǎo)作用,分形維數(shù)主要受0.038~0.25 mm粒級(jí)的影響,MWSSA不能準(zhǔn)確分析該地區(qū)的團(tuán)聚體水穩(wěn)性。

    土壤;團(tuán)聚體;植被;海拔;祁連山

    0 引 言

    土壤團(tuán)聚體是細(xì)小土粒經(jīng)過(guò)有機(jī)和礦物質(zhì)的膠結(jié)、離子凝聚及團(tuán)聚作用而形成的,是直徑<10 mm的土壤結(jié)構(gòu)單位[1]。它不僅對(duì)土壤水分運(yùn)移和儲(chǔ)存、土壤孔隙狀況、固碳、生物活動(dòng)、土壤抗蝕能力、根系生長(zhǎng)等方面至關(guān)重要[2-4],而且能夠通過(guò)抵抗雨滴打擊和地表徑流侵蝕作用,防止土壤養(yǎng)分流失[5],有效改善土壤功能和環(huán)境[6],是評(píng)價(jià)土體安全、土壤質(zhì)量和抗蝕性能的關(guān)鍵指標(biāo)[7]。

    許多學(xué)者以水穩(wěn)定性團(tuán)聚體質(zhì)量百分?jǐn)?shù)(percentage of water-stable aggregates,WSA)、平均質(zhì)量直徑(mean weight diameter,MWD)、幾何平均直徑(geometric mean diameter,GMD)作為土壤團(tuán)聚體評(píng)價(jià)指標(biāo)[8-10]。WSA直接表征土壤團(tuán)聚體穩(wěn)定性,但不足以反映其他粒級(jí)對(duì)團(tuán)聚體的影響;MWD、GMD表征土壤團(tuán)聚體的團(tuán)聚度和大小分布狀況,其值越大,土壤平均粒徑的團(tuán)聚度越高,穩(wěn)定性越好[10]。平均重量比表面積(mean weight soil specific area,MWSSA)和結(jié)構(gòu)體破壞率(percentage of aggregate destruction,PAD)也是研究團(tuán)聚體特征的有效工具,分別反映土壤團(tuán)聚體的分散性和其受到侵蝕后崩解、破壞的程度[11],但MWSSA在反映土壤團(tuán)聚體質(zhì)量的準(zhǔn)確性和靈敏性仍不清楚。蘇靜等[12]通過(guò)比較黃土高原土壤分形維數(shù)(fractal dimension,)、WSA和MWD 等團(tuán)聚體評(píng)價(jià)方法,發(fā)現(xiàn)評(píng)價(jià)土壤結(jié)構(gòu)較好,還能反映土壤理化性質(zhì)和肥力狀況。目前,國(guó)內(nèi)外學(xué)者主要在黃土高原、紫色丘陵及南方紅壤等易侵蝕區(qū),從土地利用模式、施肥方式、植被恢復(fù)、植被群落類(lèi)型等對(duì)土壤團(tuán)聚體特征的影響等方面進(jìn)行了大量的研究,在土壤團(tuán)聚體分布規(guī)律、土壤穩(wěn)定性、土壤結(jié)構(gòu)改善方面取得了豐富的研究成果[5-9]。劉雷等[5, 9-10]研究了黃土高原土壤團(tuán)聚體穩(wěn)定性的變化特征及其影響因素,指出MWD空間變化受有機(jī)質(zhì)、土壤質(zhì)地、根系質(zhì)量密度、氣候和植被類(lèi)型等因素的綜合影響。于海艷等[8]探討了不同林分類(lèi)型團(tuán)聚體穩(wěn)定性隨土壤深度的變化特征,研究表明團(tuán)聚體穩(wěn)定性隨土壤深度增加而逐漸降低。謝錦升等[6-7, 11]研究了植被恢復(fù)、植物籬等對(duì)土壤團(tuán)聚體穩(wěn)定性影響,提出植物籬有利于增強(qiáng)團(tuán)聚體穩(wěn)定性,團(tuán)聚體是反映土壤質(zhì)量、抗蝕性能的重要指標(biāo)。

    作為生態(tài)安全屏障,祁連山區(qū)受氣候變暖和人為破壞的雙重影響,土壤侵蝕加劇,水土流失嚴(yán)重,土壤理化性質(zhì)改變[13],是中國(guó)極強(qiáng)度生態(tài)脆弱區(qū)之一[14]。祁連山區(qū)地貌特征復(fù)雜,水熱條件差異顯著,表現(xiàn)為從東到西土壤質(zhì)地粗化,植被帶從溫性草原帶過(guò)渡到高寒荒漠帶,年降雨量從東到西遞減并隨海拔升高而增大[15]。土壤理化性質(zhì)、植被分布和降雨量等空間變化可能引起土壤結(jié)構(gòu)和團(tuán)聚體穩(wěn)定性發(fā)生變化。因此研究祁連山區(qū)不同植被類(lèi)型下土壤團(tuán)聚體空間變異,具有生產(chǎn)價(jià)值和理論意義。但由于祁連山區(qū)獨(dú)特的自然地理?xiàng)l件和氣候條件,有關(guān)祁連山區(qū)的團(tuán)聚體研究鮮有報(bào)道。鑒于此,本文以“絲綢之路經(jīng)濟(jì)帶”重點(diǎn)保護(hù)區(qū)——祁連山區(qū)為對(duì)象,通過(guò)分析土壤團(tuán)聚體WSA、MWD、GMD、MWSSA、PAD和等指標(biāo)的變化趨勢(shì),研究祁連山區(qū)不同植被類(lèi)型、不同海拔、不同土壤深度下土壤團(tuán)聚體變化特征,以期為該區(qū)土壤質(zhì)量改良提供理論依據(jù),為水土保持、生態(tài)恢復(fù)治理工作提供參考。

    1 材料與方法

    1.1 研究區(qū)概況

    祁連山地處于青藏高原東北緣(35°48′~40°05′ N,93°18′~103°54′ E),地勢(shì)西高東低,海拔為1 623~5 767 m,大部分在3 500~5 000 m[15],總面積1.84×105km2,草原和草甸占總土地面積的60%以上[16]。祁連山區(qū)是西北地區(qū)主要水源地,也是“絲綢之路經(jīng)濟(jì)帶”重點(diǎn)保護(hù)區(qū)[14]。該區(qū)屬大陸性高寒半干旱山地氣候區(qū),年平均氣溫為0.6 ℃,年降雨量為400~700 mm。祁連山區(qū)植被分布具有明顯的水平地帶性,自東向西依次為溫性草原、溫帶針葉林、寒溫性針葉林、高寒灌叢、高寒草甸、高寒草原、高寒荒漠[15]。與植被分布相對(duì)應(yīng),該區(qū)主要土壤類(lèi)型有黑氈土、黑鈣土、高寒草甸土、棕漠土、冷鈣土和栗鈣土等。

    1.2 樣品采集與分析

    在2018年6-8月對(duì)祁連山區(qū)野外考察的基礎(chǔ)上選取典型樣地,樣地主要選擇地勢(shì)平坦且能代表荒漠、草原、草甸和灌叢4種植被分布的區(qū)域,共設(shè)置了56個(gè)采樣點(diǎn)(圖1),每個(gè)采樣點(diǎn)采用土鉆(直徑為10 cm)取0~10、>10~20、>20~30、>30~40 cm深度的土壤樣品。每個(gè)樣地按“十字型”布設(shè)3個(gè)樣方,設(shè)置3個(gè)重復(fù),將土壤剖面相同土層的重復(fù)樣品混合均勻裝袋。

    將已編號(hào)的原狀土樣帶回實(shí)驗(yàn)室進(jìn)行前期處理,放置通風(fēng)陰涼處自然風(fēng)干,并除去動(dòng)植物殘?bào)w、石礫等。土壤有機(jī)碳(soil organic carbon,SOC)采用重鉻酸鉀外加熱法測(cè)定。水穩(wěn)性團(tuán)聚體含量的測(cè)定綜合參考Kemper等[17-18]所采用的濕篩法和分離方法,稱取過(guò)6 mm篩孔的風(fēng)干土樣25 g,放入TTF-100型團(tuán)聚體分析儀(浙江省上虞市舜戈實(shí)驗(yàn)儀器廠)的套篩(自上而下為4、2、1、0.25和0.038 mm孔徑)中,將套篩置于清水中快速浸潤(rùn)10 min,再以30次/min的頻率振動(dòng)10 min,套篩里的土樣用布氏漏斗抽濾,于105 ℃下烘干稱質(zhì)量。將各烘干土樣掃入裝有5 g/L六偏磷酸鈉的鋁盒中,再放進(jìn)恒溫?fù)u床內(nèi)以150 rpm震蕩12~16 h,然后用對(duì)應(yīng)孔徑的篩子和布氏漏斗抽濾,于105 ℃下烘干稱質(zhì)量,則測(cè)得各級(jí)孔徑團(tuán)聚體質(zhì)量。土壤團(tuán)聚體穩(wěn)定性指標(biāo)WSA、MWD、GMD計(jì)算公式如下:

    平均重量比表面積(MWSSA,cm2/g)表征土壤團(tuán)聚體的分散性,計(jì)算公式如(4):

    式中為2.65 g/cm3。

    根據(jù)楊培嶺等[19]的推導(dǎo)公式計(jì)算分形維數(shù)為

    土壤結(jié)構(gòu)結(jié)構(gòu)體破壞率(PAD,%)計(jì)算公式如下:

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

    采用Microsoft Excel 2010和IBM SPSS Statistics 20統(tǒng)計(jì)軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。采用單因素方差分析(one-way ANOVA)和最小顯著差異法(least-significant difference,LSD)比較相同土壤深度不同植被類(lèi)型間的差異;采用雙因素方差分析比較不同土壤深度間的差異(=0.05);采用Pearson相關(guān)分析對(duì)土壤團(tuán)聚指標(biāo)、有機(jī)碳和各粒級(jí)團(tuán)聚體百分含量進(jìn)行相關(guān)性評(píng)價(jià)。

    2 結(jié)果與分析

    2.1 不同植被類(lèi)型土壤團(tuán)聚體變化差異

    土壤WSA、MWD、GMD、MWSSA、PAD和的最大值與最小值之比分別為3.57、5.89、8.24、44.12、28.23和2.40倍(表1),研究區(qū)土壤團(tuán)聚體空間變化大,分布不均。變異系數(shù)(coefficient of variation,CV)可以表征空間變異程度的大小。根據(jù)Nielsen等[20]的分級(jí)標(biāo)準(zhǔn),CV≤10%表示弱變異性,10%<CV<100%表示中等變異性,CV≥100%表示強(qiáng)變異性。土壤WSA、MWD、GMD、MWSSA、PAD和的變異系數(shù)分別為17.86%、36.66%、38.23%、62.50%、45.38%和10.91%,表明祁連山區(qū)土壤團(tuán)聚體穩(wěn)定性的空間變異性為中等變異性。

    由圖2可知,同土層不同植被類(lèi)型下團(tuán)聚體指標(biāo)存在差異。0~10 cm土層各植被類(lèi)型間WSA平均值依次:荒漠<草原(或草甸)<灌叢,>10~20 cm土層草甸的WSA顯著高于荒漠(<0.05),>20~40 cm土層各植被類(lèi)型間WSA無(wú)顯著差異。0~30 cm土層各植被類(lèi)型間MWD平均值和GMD平均值都體現(xiàn):荒漠<草原<草甸<灌叢,>30~40 cm土層各植被類(lèi)型間MWD、GMD無(wú)顯著差異。0~40 cm各土層荒漠土壤MWSSA顯著低于草甸、灌叢(<0.05),而草甸與灌叢同土層MWSSA無(wú)顯著差異。0~10 cm土層土壤PAD平均值荒漠顯著高于草原、灌叢、草甸,>10~30 cm土層草甸PAD顯著低于荒漠、草原和灌叢(<0.05),>30~40 cm土層各植被類(lèi)型間PAD無(wú)顯著差異。圖2f可知,0~10 cm土層平均值荒漠和草原顯著高于灌叢和草甸,>10~20 cm土層草甸顯著低于荒漠(<0.05),>20~40 cm土層各植被類(lèi)型間無(wú)顯著差異。

    表1 祁連山區(qū)土壤團(tuán)聚體指標(biāo)統(tǒng)計(jì)

    注:小寫(xiě)字母表示同土層不同植被類(lèi)型間顯著差異(P<0.05),大寫(xiě)字母表示同植被類(lèi)型下不同土層間顯著差異(P<0.05)。

    2.2 不同植被類(lèi)型團(tuán)聚體隨土壤深度變化特征

    雙因素方差分析表明(表2),土壤深度對(duì)WSA有顯著影響(<0.05),對(duì)土壤MWD、GMD、PAD和存在極顯著影響(<0.01),但對(duì)土壤MWSSA無(wú)顯著影響(>0.05)。植被類(lèi)型與土壤深度的交互效應(yīng)對(duì)土壤MWD和GMD有極顯著影響(<0.01),而對(duì)其他團(tuán)聚體指標(biāo)無(wú)顯著影響(>0.05)。

    表2 不同植被類(lèi)型及土壤深度下土壤團(tuán)聚體指標(biāo)雙因素方差分析

    注:*,<0.05;**,<0.01;***,<0.001。

    Note: *,<0.05; **,<0.01; ***,<0.001.

    由圖2可知,同一植被類(lèi)型不同土壤深度下團(tuán)聚體變化規(guī)律不統(tǒng)一。隨土壤深度的增加,WSA逐漸減小,僅草甸帶>30~40與0~10、>10~20、>20~30 cm土層間WSA存在顯著差異(<0.05);隨土壤深度的增加,MWD、GMD均逐漸減小,草甸、灌叢帶>30~40和0~10、>10~20 cm土層間MWD、GMD差異顯著(<0.05);各土層間MWSSA無(wú)顯著差異;隨土壤深度的增加,PAD逐漸增大,在草原和草甸帶0~10 cm土層PAD顯著低于>20~30、>30~40 cm土層(<0.05);隨土壤深度的增加,逐漸增大,草甸帶0~10、>10~20和>20~30、>30~40 cm土層間差異顯著(<0.05)。

    2.3 不同植被類(lèi)型團(tuán)聚體隨海拔變化特征

    因表層采集土壤樣品較下層多,本文選用不同植被類(lèi)型表層(0~10 cm)土壤研究團(tuán)聚體隨海拔變化特征。采集土樣的海拔為1 692~3 639 m(圖3),1 692~3 200 m低海拔上發(fā)育高寒荒漠和草原,在3 000~3 639 m高海拔上發(fā)育著草甸和灌叢。

    圖3 祁連山區(qū)不同植被類(lèi)型不同指標(biāo)沿海拔變化趨勢(shì)

    由圖3a可知,在海拔1 692~3 639 m,WSA隨海拔升高明顯的先逐漸增大后減小,在海拔2 800 m處最大,荒漠土壤WSA隨海拔()變化的擬合方程為WSA=?4×10-5x+0.225?218.73(2=0.76,<0.01)。MWD、GMD隨海拔梯度變化趨勢(shì)基本一致(圖3b和圖 3c),在海拔1 692~2 800 m,MWD、GMD均隨海拔的升高逐漸增大,在海拔2 800~3 639 m,MWD、GMD隨海拔的升高而逐漸減小,在草甸和灌叢帶,MWD和GMD隨海拔增大而顯著降低(2=0.44,0.01;2=0.51,<0.01)。MWSSA變化有所不同(圖3d),與海拔梯度不相關(guān)。PAD變化趨勢(shì)與MWD、GMD相反(圖3e),在海拔1 692~2 800 m,PAD隨海拔的升高而減小,在海拔2 800~3 639 m,PAD隨海拔的升高逐漸增大,草原帶土壤PAD隨海拔升高而顯著降低(2=0.71,<0.01),回歸方程為PAD=?2×10-5x?0.131+233.46。由圖3f可知,在海拔1 692~2 800 m,隨海拔的升高而減小,在海拔2 800~3 639 m,隨海拔的升高逐漸增大,荒漠土壤與海拔顯著相關(guān)(R=0.57,<0.01)。

    2.4 土壤團(tuán)聚體指標(biāo)與各粒級(jí)的相關(guān)分析

    土壤是由細(xì)顆粒組成的,具有自相似性結(jié)構(gòu)和分形特征,土壤水穩(wěn)定性團(tuán)聚體指標(biāo)受各粒級(jí)直接影響。由表3可知,各粒級(jí)團(tuán)聚體含量(>4、>2~4、>1~2、>0.25~1 mm)、有機(jī)碳與WSA、MWD、GMD和MWSSA極顯著正相關(guān)(<0.01),與PAD、極顯著負(fù)相關(guān)(<0.01),用相關(guān)系數(shù)來(lái)評(píng)價(jià)土壤團(tuán)聚體粒級(jí)對(duì)土壤團(tuán)聚體的貢獻(xiàn),>2~4、>1~2 mm團(tuán)聚體含量對(duì)土壤團(tuán)聚體指標(biāo)貢獻(xiàn)較大,0.038~0.25 mm與極顯著正相關(guān)(<0.01),相關(guān)系數(shù)為0.501。說(shuō)明,土壤大團(tuán)粒結(jié)構(gòu)中1~4 mm粒級(jí)團(tuán)聚體是土壤團(tuán)聚體穩(wěn)定性、團(tuán)聚度的主導(dǎo)因素,小團(tuán)粒結(jié)構(gòu)中0.038~0.25 mm粒級(jí)團(tuán)聚體是土壤結(jié)構(gòu)的主要影響因素。MWD與GMD極顯著正相關(guān)(2=0.95),MWSSA與WSA、MWD、GMD、相關(guān)關(guān)系不顯著,說(shuō)明MWSSA不適用于研究水穩(wěn)性團(tuán)聚體特征。

    表3 土壤團(tuán)聚體各指標(biāo)與團(tuán)聚體粒級(jí)、有機(jī)質(zhì)的Pearson相關(guān)分析

    注(Note):=161。

    3 討 論

    研究區(qū)土壤團(tuán)聚體指標(biāo)均屬于中等變異,變異系數(shù)最小,為10.91%。0~10 cm土層各植被類(lèi)型WSA、MWD、GMD平均值依次:荒漠<草原(或草甸)<灌叢,PAD、平均值荒漠顯著高于草原、灌叢、草甸,說(shuō)明土壤團(tuán)聚體穩(wěn)定性及團(tuán)聚度依次:荒漠<草原(或草甸)<灌叢,土壤結(jié)構(gòu)穩(wěn)定性荒漠小于草原、灌叢、草甸,自東向西團(tuán)聚體穩(wěn)定性降低,草甸相比于荒漠、草原有助于增加土壤團(tuán)聚體,提高土壤穩(wěn)定性和團(tuán)聚度,改善土壤結(jié)構(gòu)。這是因?yàn)閃SA、MWD、GMD與植被覆蓋度、有機(jī)質(zhì)正相關(guān)[8-9],與有機(jī)碳負(fù)相關(guān)(表3),0~20 cm土層荒漠土壤團(tuán)聚體穩(wěn)定性顯著低于草甸、灌叢的(<0.05),而0~20 cm荒漠土壤PAD和顯著高于草甸、灌叢的(<0.05),原因是草甸、灌叢土壤在植被保護(hù)和有機(jī)質(zhì)膠結(jié)作用下,利于土壤團(tuán)聚體形成,降低了雨滴打擊作用引起的土壤團(tuán)聚體崩解,改善土壤質(zhì)量,而西南荒漠帶土壤受風(fēng)蝕作用,1~4 mm 粒級(jí)大團(tuán)聚體破碎,含碳量少的0.038~0.25 mm 粒級(jí)小團(tuán)聚體增多,土壤結(jié)構(gòu)穩(wěn)定性差。

    土壤深度變化對(duì)土壤團(tuán)聚體指標(biāo)存在顯著影響(MWSSA除外)。隨土壤深度增加,WSA、MWD和GMD均逐漸減小,而PAD和逐漸增大,且0~30 cm土層內(nèi)荒漠和草甸、灌叢團(tuán)聚體指標(biāo)差異顯著(>20~30 cm土層WSA和除外),>30~40 cm土層各植被類(lèi)型間團(tuán)聚體指標(biāo)無(wú)顯著差異。說(shuō)明隨土壤深度的增加,土壤團(tuán)聚體穩(wěn)定性和團(tuán)聚度逐漸減小,土壤結(jié)構(gòu)趨于惡化,植被類(lèi)型對(duì)土壤團(tuán)聚體的影響比土層更大,但在>30~40 cm土層內(nèi)對(duì)土壤團(tuán)聚體無(wú)顯著作用。這與于海艷等[8, 21]研究結(jié)果一致,WSA、MWD隨土壤剖面逐漸降低,PAD、變化趨勢(shì)與此相反。主要因?yàn)橥寥缊F(tuán)聚體穩(wěn)定性(MWD)和有機(jī)碳(表3)、根系質(zhì)量密度顯著正相關(guān)[9],祁連山草甸、灌叢地下植物量及根系主要分布在0~20 cm處[22-23],土壤有機(jī)碳富集在土壤表面,隨土壤深度的增加而降低,30 cm土層以下趨于穩(wěn)定[24-25]。而MWSSA變化趨勢(shì)不顯著,是因?yàn)樗焕趫F(tuán)聚體水穩(wěn)定性表達(dá),與王潤(rùn)澤等[11]研究結(jié)果一致。

    隨海拔的升高,在海拔1 692~2 800 m,WSA、MWD、GMD逐漸增大,PAD和逐漸減小,海拔2 800~3 639 m,WSA、MWD、GMD逐漸減小,PAD和逐漸增大。說(shuō)明低于海拔2 800 m土壤團(tuán)聚體隨海拔的升高逐漸穩(wěn)定,土壤結(jié)構(gòu)和質(zhì)地得到改善,高于2 800 m海拔土壤團(tuán)聚體隨海拔的升高逐漸被破壞,土壤結(jié)構(gòu)趨于惡化。這一變化趨勢(shì)是受微生物、有機(jī)質(zhì)和植被覆蓋度等共同影響,祁連山生物量、碳密度均隨海拔的增加呈先增加后降低趨勢(shì),在海拔2 800~2 950 m達(dá)到最大值,2 800~3 300 m海拔梯度內(nèi)生物量波動(dòng)降低,有機(jī)碳密度在海拔3 300 m急劇降低[26-27]。此外,該區(qū)在3 000~3 639 m高海拔處,土壤有機(jī)質(zhì)和植被覆蓋度較高,土壤條件較好,但有機(jī)質(zhì)只有在微生物的吸附粘結(jié)作用下具有團(tuán)聚能力[2],并不是高海拔寒區(qū)團(tuán)聚體形成的制約因素,該海拔梯度內(nèi)生物量逐漸降低是土壤團(tuán)聚體穩(wěn)定性降低和結(jié)構(gòu)惡化的原因。

    4 結(jié) 論

    1)祁連山區(qū)土壤團(tuán)聚體各指標(biāo)變化范圍大,且土壤團(tuán)聚體空間變異性屬于中等變異,這主要受植被類(lèi)型、土壤深度、海拔與有機(jī)碳的影響。植被類(lèi)型對(duì)土壤團(tuán)聚體有顯著影響(<0.01),0~10 cm土層水穩(wěn)定性團(tuán)聚體質(zhì)量百分?jǐn)?shù)(percentage of water-stable aggregates,WSA)、平均質(zhì)量直徑(mean weight diameter,MWD)、幾何平均直徑(geometric mean diameter,GMD)、平均重量比表面積(mean weight soil specific area,MWSSA)的平均值依次:荒漠<草原或草甸<灌叢,結(jié)構(gòu)體破壞率(percentage of aggregate destruction,PAD)和分形維數(shù)(fractal dimension,)平均值荒漠顯著高于草原、灌叢、草甸,草甸顯著提高(<0.05)土壤團(tuán)聚體的穩(wěn)定性,改善土壤結(jié)構(gòu),>30~40 cm土層植被類(lèi)型對(duì)土壤團(tuán)聚體無(wú)顯著影響。

    2)土壤深度變化對(duì)團(tuán)聚體有顯著影響,隨著土壤深度增大,WSA、MWD和GMD均逐漸減小,而PAD和值逐漸增大,草甸、灌叢帶0~10 cm土層團(tuán)聚體穩(wěn)定性顯著高于>30~40 cm,說(shuō)明土壤團(tuán)聚度、團(tuán)聚體及結(jié)構(gòu)的穩(wěn)定性隨土壤深度的增加逐漸減小,表層土壤團(tuán)聚體穩(wěn)定性顯著高于深層。

    3)隨海拔升高,團(tuán)聚體穩(wěn)定性和團(tuán)聚度呈先增大后降低趨勢(shì),土壤結(jié)構(gòu)穩(wěn)定性先降低后增大,海拔2 800 m為土壤團(tuán)聚體穩(wěn)定性的分界點(diǎn)。針對(duì)不同海拔土壤性質(zhì),調(diào)整植被模式,改良土壤。

    4)土壤團(tuán)聚體穩(wěn)定性受1~4 mm粒級(jí)大團(tuán)聚體的主導(dǎo)作用,與有機(jī)碳顯著正相關(guān);土壤結(jié)構(gòu)主要受0.038~0.25 mm粒級(jí)小團(tuán)聚體的影響。MWSSA不適用于研究水穩(wěn)性團(tuán)聚體特征,研究結(jié)果對(duì)于祁連山區(qū)生態(tài)恢復(fù)治理和環(huán)境評(píng)價(jià)具有重要意義。

    [1]Bronick C J, Lal R. Soil structure and management: A review[J]. Geoderma, 2005, 124(1/2): 0-22.

    [2]王清奎,汪思龍. 土壤團(tuán)聚體形成與穩(wěn)定機(jī)制及影響因素[J]. 土壤通報(bào),2005,36(3):415-421.

    Wang Qingkui, Wang Silong. Forming and stable mechanism of soil aggregate and influencing factors[J]. Chinese Journal of Soil Science, 2005, 36(3): 415-421. (in Chinese with English abstract)

    [3]He Y, Xu C, Gu F, et al. Soil aggregate stability improves greatly in response to soil water dynamics under natural rains in long-term organic fertilization[J]. Soil & Tillage Research, 2018, 184: 281-290.

    [4]A Cerdà. Aggregate stability against water forces under different climates on agriculture land and scrubland in southern Bolivia[J]. Soil & Tillage Research, 2000, 57(3): 159-166.

    [5]劉雷,安韶山,黃華偉. 應(yīng)用Le Bissonnais 法研究黃土丘陵區(qū)植被類(lèi)型對(duì)土壤團(tuán)聚體穩(wěn)定性的影響[J]. 生態(tài)學(xué)報(bào),2013,33(20):6670-6680.

    Liu Lei, An Shaoshan, Huang Huawei. Application of Le Bissonnais method to study soil aggregate stability under different vegetation on the loess plateau[J]. Acta Ecological Sinica, 2013, 33(20): 6670-6680. (in Chinese with English abstract)

    [6]謝錦升,楊玉盛,陳光水,等. 植被恢復(fù)對(duì)退化紅壤團(tuán)聚體穩(wěn)定性及碳分布的影響[J]. 生態(tài)學(xué)報(bào),2008,28(2):702-709.

    Xie Jinsheng, Yang Yusheng, Chen Guangshui, et al. Effects of vegetation restoration on water stability and organic carbon distribution in aggregates of degraded red soil in subtropics of China[J]. Acta Ecological Sinica, 2008, 28(2): 702-709. (in Chinese with English abstract)

    [7]蒲玉琳,林超文,謝德體,等. 植物籬—農(nóng)作坡地土壤團(tuán)聚體組成和穩(wěn)定性特征[J]. 應(yīng)用生態(tài)學(xué)報(bào),2013,24(1):122-128.

    Pu Yulin, Lin Chaowen, Xie Deti, et al. Composition and stability of soil aggregates in hedgerow-crop slope land[J]. Chinese Journal of Applied Ecology, 2013, 24(1): 122-128. (in Chinese with English abstract)

    [8]于海艷,宮汝寧,周婭,等. 北京八達(dá)嶺地區(qū)4種人工林土壤團(tuán)聚體穩(wěn)定性及有機(jī)碳特征[J]. 水土保持學(xué)報(bào),2015,29(5):162-166.

    Yu Haiyan, Gong Runing, Zhou Ya, et al. Characteristics of soil aggregate stability and soil organic carbon under four typical artificial plantations in Beijing Badaling Mountain area[J]. Journal of Soil and Water Conservation, 2015, 29(5): 162-166. (in Chinese with English abstract)

    [9]耿韌,張光輝,洪大林,等. 黃土高原農(nóng)地草地林地土壤團(tuán)聚體穩(wěn)定性沿降水梯度的變化特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(3):141-148.

    Geng Ren, Zhang Guanghui, Hong Dalin, et al. Variation characteristics of aggregate stability of cropland, grassland and woodland along precipitation gradient in Loess Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2019, 35(3): 141-148. (in Chinese with English abstract)

    [10]Yea L, Tana W, Fang L, et al. Spatial analysis of soil aggregate stability in a small catchment of the Loess Plateau, China: I. Spatial variability[J]. Soil & Tillage Research, 2018, 179: 71-81.

    [11]王潤(rùn)澤,諶蕓,李鐵,等. 紫色土區(qū)植物籬籬前淤積帶土壤團(tuán)聚體穩(wěn)定性特征研究[J]. 水土保持學(xué)報(bào),2018,32(2):210-216.

    Wang Runze, Chen Yun, Li Tie, et al. Study on stability of soil aggregates in the sedimentation zone in front of hedgerows in purple soil area[J]. Journal of Soil and Water Conservation, 2018, 32(2): 210-216. (in Chinese with English abstract)

    [12]蘇靜,趙世偉. 土壤團(tuán)聚體穩(wěn)定性評(píng)價(jià)方法比較[J]. 水土保持通報(bào),2009,29(5):114-117.

    Su Jing, Zhao Shiwei. Comparison of the analysis methods for soil aggregate stability[J]. Bulletin of Soil and Water Conservation, 2009, 29(5): 114-117. (in Chinese with English abstract)

    [13]Yuan H, Hou F. Grazing intensity and soil depth effects on soil properties in alpine meadow pastures of Qilian Mountain in northwest China[J]. Acta Agriculturae Scandinavica, Section B: Soil & Plant Science, 2015, 65(3): 222-232.

    [14]王濤,高峰,王寶,等. 祁連山生態(tài)保護(hù)與修復(fù)的現(xiàn)狀問(wèn)題與建議[J]. 冰川凍土,2017,39(2):229-234.

    Wang Tao, Gao Feng, Wang Bao, et al. Status and suggestions on ecological protection and restoration of Qilian Mountains[J]. Journal of Glaciology and Geocryology, 2017, 39(2): 229-234. (in Chinese with English abstract)

    [15]張富廣. 氣候變化背景下祁連山高寒草甸上界現(xiàn)代分布變化研究[D]. 甘肅:蘭州大學(xué),2018.

    Zhang Fuguang. Modern Distribution Changes in the Upper Boundary of Alpine Meadow Belt in the Qilian Mountains With Climate Warming[D]. Gansu: Lanzhou University, 2018. (in Chinese with English abstract)

    [16]He Z, Zhao W, Liu H, et al. The response of soil moisture to rainfall event size in subalpine grassland and meadows in a semi-arid mountain range: A case study in northwestern China’s Qilian Mountains[J]. Journal of Hydrology, 2012(420/421): 183-190.

    [17]Kemper W D, Rosenau R C, Klute A. Aggregate stability and size distribution[J]. Methods of Soil Analysis, 1986(1): 425-442.

    [18]Elliott E T. Aggregate structure and carbon, nitrogen, and phosphorus in native and cultivated soils[J]. Soil Science Society of America Journal, 1986, 50(3): 627-633.

    [19]楊培嶺,羅遠(yuǎn)培,石元春. 用粒徑的重量分布表征的土壤分形特征[J]. 科學(xué)通報(bào),1993,38(20):1896-1899.

    Yang Peiling, Luo Yuanpei, Shi Yuanchun. Soil fractal character token by particle-mass distribution[J]. Chinese Science Bulletion, 1993, 38(20): 1896-1899. (in Chinese with English abstract)

    [20]Nielsen D R, Bouma J. Soil spatial variability[C]// Proceedings of a Workshop of the ISSS and the SSSA, Las Vegas, USA. 30th November to 1st December, 1984. Pudoc, Wageningen: The Netherlands,1985: 243.

    [21]馬帥,趙世偉,李婷,等. 子午嶺林區(qū)植被自然恢復(fù)下土壤剖面團(tuán)聚體特征研究[J]. 水土保持學(xué)報(bào),2011,25(2):157-161.

    Ma Shuai, Zhao Shiwei, Li Ting, et al. Study of the character of soil aggregates under vegetation restoration in Ziwu Mountain[J]. Journal of Soil and Water Conservation, 2011, 25(2): 157-161. (in Chinese with English abstract)

    [22]于應(yīng)文,胡自治,張德罡,等. 金露梅株叢的植物量及其空間分布的研究[J]. 甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào),1999,34(3):237-242.

    Yu Yingwen, Hu Zizhi, Zhang Degang, et al. Study on phytomass and its distribution of Potentilla fruticosa[J]. Journal of Gansu Agricultural University, 1999, 34(3): 237-242. (in Chinese with English abstract)

    [23]史麗麗,趙成章,樊潔平,等. 祁連山地甘肅臭草斑塊土壤水分與植被蓋度空間格局[J]. 生態(tài)學(xué)雜志,2013,32(2):285-291.

    Shi Lili, Zhao Chengzhang, Fan Jieping, et al. Analysis of the soil moisture and Melica przewalskyi patches in the degraded alpine grassland of Qilian Mountains, China[J]. Chinese Journal of Ecology, 2013, 32(2): 285-291. (in Chinese with English abstract)

    [24]齊鵬,劉賢德,趙維俊,等. 祁連山中段青海云杉林土壤養(yǎng)分特征[J]. 山地學(xué)報(bào),2015,33(5):538-545.

    Qi Peng, Liu Xiande, Zhao Weijun, et al. Soil nutrient characteristics of Picea crassifolia forest in the middle segment of Qilian Mountains[J]. Mountain Research, 2015, 33(5): 538-545. (in Chinese with English abstract)

    [25]李占斌,周波,馬田田,等. 黃土丘陵區(qū)生態(tài)治理對(duì)土壤碳氮磷及其化學(xué)計(jì)量特征的影響[J]. 水土保持學(xué)報(bào),2017,31(6):312-318.

    Lin Zhanbin, Zhou Bo, Ma Tiantian, et al. Effects of ecological management on characteristics of soil carbon, nitrogen, phosphorus and their stoichiometry in Loess Hilly Region, China[J]. Journal of Soil and Water Conservation, 2017, 31(6): 312-318. (in Chinese with English abstract)

    [26]賈文雄,劉亞榮,張禹舜,等. 祁連山草甸草原物種多樣性和生物量與氣候要素的關(guān)系[J]. 干旱區(qū)研究,2015,32(6):1167-1172.

    Jia Wenxiong, Liu Yarong, Zhang Yushun, et al. Species diversity and biomass of meadow steppe in Qilian Mountains and their relationships with climate factors[J]. Arid Zone Research, 2015, 32(6): 1167-1172. (in Chinese with English abstract)

    [27]曾立雄,雷蕾,王曉榮,等. 海拔梯度對(duì)祁連山青海云杉林喬木層和土壤層碳密度的影響[J]. 生態(tài)學(xué)報(bào),2018,38(20):7168-7177.

    Zeng Lixiong, Lei Lei, Wang Xiaorong, et al. Effect of altitudinal variation on carbon density in arbor layer and soil layer of Picea crassifolia forest in Qilian Mountains[J]. Acta Ecologica Sinica, 2018, 38(20): 7168-7177. (in Chinese with English abstract)

    Variation characteristics of soil aggregates under main vegetation types in Qilian Mountainous areas

    Wei Xia1, He Yan1, Wei Ning2, Yu Wenzhu1, Cui Xia1, Zhao Hengce1

    (1.,,730000,; 2.,712100,)

    Soil aggregates reflect the stability of soil structure, soil fertility and quality, and they are closely related to the environmental quality, soil and water loss, and soil erosion. Although vegetation types, soil depths, and altitude likely have important effects on soil properties and soil aggregates, few studies have concentrated on the relationships between their interactive effects and soil aggregate stability on the regional scale. In order to study the variation characteristics of soil aggregates in the Qilian Mountains. In this study, 4 vegetation type zones (desert, steppe, meadow, and shrub) were selected in the Qilian Mountains (northwest China) as the subjects. The stability and size distribution of soil aggregates was measured by the method of wet-sieving. The characteristics of soil aggregate index were analyzed at different soil depths and different altitudes, including the percentage of water-stable aggregate (WSA), mean weight diameter (MWD), geometric mean diameter (GMD), mean weight soil specific area (MWSSA), aggregate processing damage rate ( PAD) and fractal dimension (). The results showed that the coefficients of variation of soil aggregate index under the 4 vegetation types were high (ranging from 10.91% to 62.50%), which indicated a moderate spatial variability. Among the 4 grassland types, WSA, MWD, GMD and MWSSA all showed the same increased in order: desert < steppe or meadow 30-40 cm depth, no significant differences in aggregate stability index were found among different vegetation types (>0.05). The reason why the aggregate stability of desert was lowest could be explained by the fact that the lower soil organic carbon (SOC) and biomass of desert, and the small aggregates couldn’t be formed. Along the vertical direction of the soil profiles from aboveground to underground, with the increase of soil depth (0-40 cm), WSA, MWD, GMD, and MWSSA gradually all decreased. Whereas PAD andbehaved opposite to them, and with the increase in soil depth, they gradually increased, which indicated that the stability and aggregation degree of soil aggregates decreased with soil depth. This could be explained as that the SOC appeared enrichment phenomenon in the 0-20 cm layer and decreased with soil depth. The soil aggregate index all tended to be significantly higher in surface layer than in lower layers only at the meadow (<0.05), but no significant difference in MWSSA (>0.05) were found among all soil depths. With increasing altitude, the aggregate stability (WSA, MWD, and GMD) increased gradually from 1 692 m, reached a peak at 2 800 m, and then decreased quickly. But the trends were reverse for the PAD and, and there was no significant correlation between MWSSA and altitude under different vegetation types (>0.05), which suggested that their distribution clearly showed unimodal patterns across all the altitude. This indicated that the stability and aggregation degree of soil aggregates first increased and then decreased with increasing altitude, and the maximum point appeared at elevation of 2 800 m. In addition, correlation analysis showed that WSA, MWD, and GMD were mainly affected by 1-4 mm,was mainly affected by 0.038-0.25 mm, and MWSSA couldn't accurately express the characteristics of water-stable aggregates. This research will guide the practice of reducing soil erosion for the different conditions and different vegetation types, and results have great significance for controlling grassland degradation, promoting soil structure stability and the sustainable development of grassland animal husbandry.

    soils; aggregates; vegetations; altitude; Qilian Mountains

    魏 霞,賀 燕,魏 寧,于文竹,崔 霞,趙恒策. 祁連山區(qū)主要植被類(lèi)型下土壤團(tuán)聚體變化特征[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(2):148-155.doi:10.11975/j.issn.1002-6819.2020.02.018 http://www.tcsae.org

    Wei Xia, He Yan, Wei Ning, Yu Wenzhu, Cui Xia, Zhao Hengce. Variation characteristics of soil aggregates under main vegetation types in Qilian Mountainous areas[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(2): 148-155. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.02.018 http://www.tcsae.org

    2019-05-07

    2019-12-10

    國(guó)家自然科學(xué)基金項(xiàng)目(51679115);蘭州大學(xué)西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金及蘭州大學(xué)中央高校基本科研業(yè)務(wù)費(fèi)專(zhuān)項(xiàng)基金(lzujbky-2018-kb01)

    魏 霞,博士后,副教授,主要從事土壤侵蝕與水土保持、水文學(xué)及水資源學(xué)等方面的研究。Email:weix@lzu.edu.cn

    10.11975/j.issn.1002-6819.2020.02.018

    S157

    A

    1002-6819(2020)-02-0148-08

    猜你喜歡
    土壤結(jié)構(gòu)祁連山草甸
    高山草甸
    圖志
    發(fā)展(2024年1期)2024-04-25 03:45:30
    The world's narrowest river
    山地草甸
    祁連山下
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    祁連山草原:如夢(mèng)如幻近高天
    祁連山
    黃河之聲(2016年20期)2016-02-21 11:55:33
    論太原市櫻花栽培情況
    土壤結(jié)構(gòu)三維可視化研究
    地球(2015年6期)2015-03-31 07:03:10
    欧美精品啪啪一区二区三区| 成人特级黄色片久久久久久久| 少妇人妻精品综合一区二区 | 少妇熟女aⅴ在线视频| 国产单亲对白刺激| 亚洲黑人精品在线| 国产91精品成人一区二区三区| 国产乱人视频| 国产探花极品一区二区| 中文字幕人妻熟人妻熟丝袜美 | 欧美日韩国产亚洲二区| 免费看日本二区| 桃红色精品国产亚洲av| 国产国拍精品亚洲av在线观看 | 国产av不卡久久| 网址你懂的国产日韩在线| 两人在一起打扑克的视频| 久久久久国内视频| 日韩欧美三级三区| 午夜影院日韩av| 国产精品久久久久久久电影 | 久久婷婷人人爽人人干人人爱| 精品电影一区二区在线| 99久久九九国产精品国产免费| 国产av一区在线观看免费| 嫁个100分男人电影在线观看| 成年女人永久免费观看视频| 色噜噜av男人的天堂激情| 日韩中文字幕欧美一区二区| 免费人成在线观看视频色| 成人鲁丝片一二三区免费| 国模一区二区三区四区视频| 国产黄a三级三级三级人| 99国产极品粉嫩在线观看| 老司机在亚洲福利影院| 国产精品电影一区二区三区| 精品熟女少妇八av免费久了| 中文字幕熟女人妻在线| 99久久精品热视频| 在线天堂最新版资源| 免费一级毛片在线播放高清视频| 在线观看日韩欧美| 欧美激情在线99| 精品久久久久久久毛片微露脸| 久久久成人免费电影| 精品久久久久久久久久免费视频| www.色视频.com| 国产伦在线观看视频一区| 久久久国产精品麻豆| 亚洲av五月六月丁香网| 国产又黄又爽又无遮挡在线| 91久久精品电影网| 熟女人妻精品中文字幕| 丰满乱子伦码专区| 少妇丰满av| 老熟妇仑乱视频hdxx| 国产毛片a区久久久久| av专区在线播放| 亚洲精品乱码久久久v下载方式 | 亚洲国产日韩欧美精品在线观看 | 久久久久亚洲av毛片大全| 白带黄色成豆腐渣| 午夜免费男女啪啪视频观看 | 国产日本99.免费观看| 禁无遮挡网站| 18禁美女被吸乳视频| 国产精品一及| 村上凉子中文字幕在线| 亚洲欧美日韩东京热| 亚洲精品乱码久久久v下载方式 | 香蕉久久夜色| 制服人妻中文乱码| 色综合欧美亚洲国产小说| 亚洲最大成人手机在线| 叶爱在线成人免费视频播放| 日韩欧美一区二区三区在线观看| 免费一级毛片在线播放高清视频| 欧美另类亚洲清纯唯美| 久久久久性生活片| 一本久久中文字幕| 亚洲无线在线观看| 成人三级黄色视频| 成人特级黄色片久久久久久久| 国产伦在线观看视频一区| 啦啦啦韩国在线观看视频| 制服人妻中文乱码| 欧美成人一区二区免费高清观看| 淫秽高清视频在线观看| 国产蜜桃级精品一区二区三区| 国产野战对白在线观看| 成年女人看的毛片在线观看| www日本在线高清视频| 老汉色∧v一级毛片| 97人妻精品一区二区三区麻豆| 婷婷精品国产亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 男女那种视频在线观看| 免费在线观看成人毛片| 99精品久久久久人妻精品| 18禁裸乳无遮挡免费网站照片| 欧美一区二区亚洲| 国产黄a三级三级三级人| 亚洲,欧美精品.| 最新中文字幕久久久久| 女人高潮潮喷娇喘18禁视频| 国产毛片a区久久久久| 内射极品少妇av片p| 成人欧美大片| 精品久久久久久成人av| 国产日本99.免费观看| 国产亚洲av嫩草精品影院| 午夜亚洲福利在线播放| 久久久精品欧美日韩精品| 亚洲18禁久久av| 热99在线观看视频| 无遮挡黄片免费观看| 国产99白浆流出| 国产毛片a区久久久久| 在线国产一区二区在线| 精品久久久久久久久久久久久| 国产精品久久电影中文字幕| 欧美中文日本在线观看视频| 成人性生交大片免费视频hd| 国产精品久久视频播放| 深爱激情五月婷婷| 琪琪午夜伦伦电影理论片6080| 亚洲av成人av| 久久精品国产自在天天线| 男女午夜视频在线观看| 精品不卡国产一区二区三区| 久久久久国产精品人妻aⅴ院| 在线国产一区二区在线| 国产欧美日韩精品一区二区| 日韩大尺度精品在线看网址| 亚洲国产精品久久男人天堂| 精品午夜福利视频在线观看一区| 可以在线观看的亚洲视频| 欧美黑人欧美精品刺激| 久久久久国产精品人妻aⅴ院| 国产激情偷乱视频一区二区| 午夜久久久久精精品| xxxwww97欧美| 欧美色欧美亚洲另类二区| 97超级碰碰碰精品色视频在线观看| 99riav亚洲国产免费| 国产精华一区二区三区| 日日干狠狠操夜夜爽| 成年女人永久免费观看视频| 少妇的逼好多水| 一区福利在线观看| 国产精品影院久久| 悠悠久久av| 三级国产精品欧美在线观看| 日本黄大片高清| 亚洲中文字幕一区二区三区有码在线看| 亚洲 欧美 日韩 在线 免费| 国产视频一区二区在线看| av福利片在线观看| 国产老妇女一区| 国产精品久久久久久久电影 | 露出奶头的视频| 麻豆国产97在线/欧美| 久久久久久国产a免费观看| 成人特级av手机在线观看| АⅤ资源中文在线天堂| 一区福利在线观看| 久久久成人免费电影| 欧美黑人欧美精品刺激| 长腿黑丝高跟| 国产精品自产拍在线观看55亚洲| 岛国视频午夜一区免费看| 午夜免费观看网址| 久久精品夜夜夜夜夜久久蜜豆| 男女下面进入的视频免费午夜| 伊人久久大香线蕉亚洲五| 久久精品人妻少妇| 美女免费视频网站| eeuss影院久久| 国产蜜桃级精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 男女做爰动态图高潮gif福利片| 久久99热这里只有精品18| 十八禁网站免费在线| 久久久久久大精品| 天美传媒精品一区二区| 香蕉丝袜av| 超碰av人人做人人爽久久 | 欧美日韩福利视频一区二区| 国产精品 国内视频| 一夜夜www| 色播亚洲综合网| 日韩欧美免费精品| 两个人视频免费观看高清| a在线观看视频网站| 欧美成狂野欧美在线观看| 好男人在线观看高清免费视频| 非洲黑人性xxxx精品又粗又长| 黄色片一级片一级黄色片| 中出人妻视频一区二区| 国产精品永久免费网站| 国产探花在线观看一区二区| 男女那种视频在线观看| 国产激情欧美一区二区| 亚洲国产欧美网| 国产精品1区2区在线观看.| 午夜福利在线观看免费完整高清在 | xxx96com| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 欧美日韩综合久久久久久 | 日本黄色视频三级网站网址| 日韩中文字幕欧美一区二区| 日韩人妻高清精品专区| 婷婷精品国产亚洲av在线| 啦啦啦免费观看视频1| 久久99热这里只有精品18| 国产成人福利小说| 香蕉丝袜av| 一个人观看的视频www高清免费观看| 18禁裸乳无遮挡免费网站照片| 国产真实乱freesex| 国产麻豆成人av免费视频| www国产在线视频色| 欧美国产日韩亚洲一区| 国内精品久久久久精免费| 91九色精品人成在线观看| 欧美日韩精品网址| 人妻夜夜爽99麻豆av| 国产野战对白在线观看| 在线观看一区二区三区| 国产成年人精品一区二区| 国产精品,欧美在线| 欧美性猛交黑人性爽| 少妇丰满av| 舔av片在线| 老司机午夜福利在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 色噜噜av男人的天堂激情| 在线观看日韩欧美| 国产亚洲精品av在线| 日本与韩国留学比较| 高清毛片免费观看视频网站| 成人永久免费在线观看视频| 国产一区二区激情短视频| 亚洲国产中文字幕在线视频| 人人妻人人澡欧美一区二区| 18禁黄网站禁片午夜丰满| 国产精品1区2区在线观看.| 露出奶头的视频| 制服人妻中文乱码| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 国内少妇人妻偷人精品xxx网站| 男女视频在线观看网站免费| 欧美绝顶高潮抽搐喷水| 嫩草影院精品99| 国产高清三级在线| 他把我摸到了高潮在线观看| 特大巨黑吊av在线直播| 日本五十路高清| 亚洲最大成人手机在线| 国产精品亚洲一级av第二区| 岛国在线免费视频观看| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 精品日产1卡2卡| 男人舔奶头视频| 男插女下体视频免费在线播放| 一本综合久久免费| АⅤ资源中文在线天堂| 精品久久久久久久久久免费视频| 国产精品久久久久久精品电影| 日本黄色片子视频| 人人妻,人人澡人人爽秒播| 男插女下体视频免费在线播放| 成人高潮视频无遮挡免费网站| 99在线人妻在线中文字幕| 色尼玛亚洲综合影院| 精品人妻偷拍中文字幕| 十八禁人妻一区二区| 伊人久久精品亚洲午夜| 国产一区二区三区在线臀色熟女| 午夜福利在线在线| 一进一出抽搐动态| 成人特级黄色片久久久久久久| 欧美成人性av电影在线观看| 国产午夜精品论理片| 变态另类丝袜制服| 中文字幕久久专区| 国产精品,欧美在线| 人妻久久中文字幕网| 神马国产精品三级电影在线观看| 日韩有码中文字幕| 国产高清视频在线播放一区| tocl精华| 欧美日韩黄片免| 色播亚洲综合网| 欧美乱色亚洲激情| 精品久久久久久久毛片微露脸| 国内精品久久久久精免费| 草草在线视频免费看| 夜夜夜夜夜久久久久| 悠悠久久av| 我要搜黄色片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 狂野欧美激情性xxxx| 国产在线精品亚洲第一网站| 一级黄色大片毛片| 男人舔女人下体高潮全视频| 俄罗斯特黄特色一大片| 国产精品电影一区二区三区| 长腿黑丝高跟| 看免费av毛片| 国产一区二区在线观看日韩 | 手机成人av网站| 嫩草影院精品99| 丁香六月欧美| 国产一区二区亚洲精品在线观看| 国产精品电影一区二区三区| 天堂√8在线中文| 精品久久久久久久久久久久久| www日本黄色视频网| 桃色一区二区三区在线观看| 此物有八面人人有两片| 狂野欧美白嫩少妇大欣赏| 久久精品91蜜桃| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 国产精品三级大全| 国产精品电影一区二区三区| 国产成人福利小说| 日韩大尺度精品在线看网址| 变态另类丝袜制服| 国产激情欧美一区二区| 亚洲国产欧美人成| 免费观看精品视频网站| АⅤ资源中文在线天堂| 国产成人a区在线观看| 男女做爰动态图高潮gif福利片| 中出人妻视频一区二区| 日韩 欧美 亚洲 中文字幕| 成人鲁丝片一二三区免费| 又粗又爽又猛毛片免费看| av女优亚洲男人天堂| 精品一区二区三区视频在线观看免费| 桃红色精品国产亚洲av| 成年女人看的毛片在线观看| 欧美成狂野欧美在线观看| 在线天堂最新版资源| 欧美成人一区二区免费高清观看| 在线观看66精品国产| 国产精品久久久久久人妻精品电影| 亚洲一区二区三区不卡视频| 偷拍熟女少妇极品色| 亚洲国产日韩欧美精品在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品精品国产色婷婷| 亚洲av中文字字幕乱码综合| 露出奶头的视频| 亚洲av中文字字幕乱码综合| 精品久久久久久久毛片微露脸| 国产伦精品一区二区三区四那| 十八禁网站免费在线| 国产精品一区二区免费欧美| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 免费高清视频大片| 亚洲av中文字字幕乱码综合| 精品久久久久久久毛片微露脸| 久久久久久久久中文| 成人特级av手机在线观看| 麻豆一二三区av精品| 麻豆国产97在线/欧美| av专区在线播放| 床上黄色一级片| 啦啦啦韩国在线观看视频| 国产av不卡久久| 精品国内亚洲2022精品成人| 国产一区二区激情短视频| 窝窝影院91人妻| 国产精品三级大全| 亚洲成av人片在线播放无| 精品久久久久久久毛片微露脸| 在线国产一区二区在线| 制服丝袜大香蕉在线| 91在线观看av| 午夜亚洲福利在线播放| 久久久久国产精品人妻aⅴ院| 欧美zozozo另类| 亚洲最大成人中文| 国产精品 国内视频| 久久国产精品影院| 久久九九热精品免费| 不卡一级毛片| 51国产日韩欧美| 国产精品久久久久久亚洲av鲁大| 女人被狂操c到高潮| 国产野战对白在线观看| 亚洲av电影不卡..在线观看| www.色视频.com| 最近最新中文字幕大全电影3| 99热只有精品国产| 国产激情欧美一区二区| 嫩草影院精品99| 国产成+人综合+亚洲专区| 天美传媒精品一区二区| 国产精品98久久久久久宅男小说| 中国美女看黄片| 成人永久免费在线观看视频| 欧美性猛交╳xxx乱大交人| 欧美最黄视频在线播放免费| 天天添夜夜摸| 日韩精品中文字幕看吧| 搡老妇女老女人老熟妇| 亚洲专区国产一区二区| 一本一本综合久久| 午夜福利在线在线| 老司机福利观看| 国产蜜桃级精品一区二区三区| 亚洲天堂国产精品一区在线| 国产成人福利小说| 欧美日韩国产亚洲二区| 一个人观看的视频www高清免费观看| 黄色成人免费大全| 桃红色精品国产亚洲av| 91字幕亚洲| 色视频www国产| eeuss影院久久| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 日韩大尺度精品在线看网址| 嫩草影院精品99| 内射极品少妇av片p| 一个人观看的视频www高清免费观看| 亚洲专区国产一区二区| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 在线免费观看不下载黄p国产 | 国产精品综合久久久久久久免费| xxxwww97欧美| 国产高清激情床上av| 国产av一区在线观看免费| 亚洲熟妇中文字幕五十中出| aaaaa片日本免费| 亚洲专区国产一区二区| 麻豆国产av国片精品| 欧美又色又爽又黄视频| 中出人妻视频一区二区| 免费电影在线观看免费观看| 看片在线看免费视频| 亚洲av一区综合| 亚洲av五月六月丁香网| 日韩亚洲欧美综合| 2021天堂中文幕一二区在线观| 国产精品久久久久久久久免 | 国产国拍精品亚洲av在线观看 | 十八禁网站免费在线| 午夜免费男女啪啪视频观看 | 亚洲 欧美 日韩 在线 免费| 99精品久久久久人妻精品| av中文乱码字幕在线| 人妻久久中文字幕网| 免费电影在线观看免费观看| 精品久久久久久久人妻蜜臀av| 欧美日韩乱码在线| 女同久久另类99精品国产91| 51国产日韩欧美| 丁香欧美五月| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 又黄又爽又免费观看的视频| 久久九九热精品免费| 亚洲av成人不卡在线观看播放网| 亚洲av不卡在线观看| 久久精品国产清高在天天线| 国产激情欧美一区二区| 搡老妇女老女人老熟妇| 亚洲国产高清在线一区二区三| 国产成人av教育| av国产免费在线观看| 天堂影院成人在线观看| 成人特级av手机在线观看| 国产熟女xx| or卡值多少钱| 9191精品国产免费久久| 欧美日韩中文字幕国产精品一区二区三区| 可以在线观看毛片的网站| 亚洲国产欧美人成| 中文字幕高清在线视频| 国产精品自产拍在线观看55亚洲| 久久精品国产综合久久久| 极品教师在线免费播放| 国产伦一二天堂av在线观看| av在线天堂中文字幕| 一个人看视频在线观看www免费 | 成人国产综合亚洲| 叶爱在线成人免费视频播放| 国产探花极品一区二区| 成人午夜高清在线视频| 757午夜福利合集在线观看| 两个人视频免费观看高清| 岛国在线观看网站| 欧美激情在线99| 欧美激情久久久久久爽电影| 亚洲欧美日韩卡通动漫| 一区二区三区国产精品乱码| 国产亚洲欧美在线一区二区| 久久久久久久亚洲中文字幕 | 动漫黄色视频在线观看| 男插女下体视频免费在线播放| 亚洲人成网站在线播放欧美日韩| 91久久精品国产一区二区成人 | 婷婷六月久久综合丁香| 精品久久久久久成人av| 男女下面进入的视频免费午夜| 欧美极品一区二区三区四区| 少妇丰满av| 国产精品 欧美亚洲| 热99在线观看视频| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 亚洲人成网站在线播| avwww免费| 亚洲不卡免费看| 草草在线视频免费看| 在线播放国产精品三级| 免费av毛片视频| xxx96com| 波多野结衣高清作品| 国产色婷婷99| 亚洲av成人av| 亚洲中文日韩欧美视频| 露出奶头的视频| 一本综合久久免费| 午夜免费激情av| 国产男靠女视频免费网站| 一本精品99久久精品77| 国内毛片毛片毛片毛片毛片| 精品不卡国产一区二区三区| 亚洲成a人片在线一区二区| 国产精品久久久久久精品电影| 18+在线观看网站| 老鸭窝网址在线观看| 伊人久久大香线蕉亚洲五| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品国产精品久久久不卡| 欧美精品啪啪一区二区三区| 久久久国产成人免费| 91av网一区二区| 午夜精品一区二区三区免费看| 日本a在线网址| 国产极品精品免费视频能看的| 亚洲精品456在线播放app | 免费av毛片视频| 18+在线观看网站| 国模一区二区三区四区视频| 男女午夜视频在线观看| 精品一区二区三区av网在线观看| 男女之事视频高清在线观看| 一区二区三区高清视频在线| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产精品麻豆| 小说图片视频综合网站| 色视频www国产| 波多野结衣高清作品| 国产精品国产高清国产av| 熟女人妻精品中文字幕| 搡女人真爽免费视频火全软件 | 欧美黑人欧美精品刺激| 亚洲欧美日韩东京热| 欧美日韩一级在线毛片| 男女下面进入的视频免费午夜| 波多野结衣巨乳人妻| 一边摸一边抽搐一进一小说| 久久国产乱子伦精品免费另类| 国内精品久久久久精免费| 欧美zozozo另类| 日韩欧美国产一区二区入口| 欧美大码av| 亚洲av第一区精品v没综合| 欧美一区二区精品小视频在线| 欧美极品一区二区三区四区| 久久久久久大精品| 全区人妻精品视频| 精品日产1卡2卡| 亚洲无线在线观看| 中文亚洲av片在线观看爽| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩精品亚洲av| 亚洲色图av天堂| 美女 人体艺术 gogo| 久久中文看片网| 亚洲国产精品合色在线| 国产中年淑女户外野战色| 亚洲,欧美精品.| 岛国在线观看网站| 美女 人体艺术 gogo| 深夜精品福利| 欧美日韩精品网址| 在线a可以看的网站| 首页视频小说图片口味搜索| 19禁男女啪啪无遮挡网站| 精品不卡国产一区二区三区| 国产单亲对白刺激| 国产亚洲欧美在线一区二区| 精品不卡国产一区二区三区| 蜜桃亚洲精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲一级av第二区| 亚洲欧美精品综合久久99| 19禁男女啪啪无遮挡网站| 国产视频一区二区在线看| 99久国产av精品| 亚洲av电影在线进入| 动漫黄色视频在线观看|