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

    高寒草地土壤理化性質(zhì)特征及質(zhì)量評(píng)價(jià)

    2022-05-11 10:44:40楊沖王文穎劉攀周華坤毛旭鋒
    草原與草坪 2022年1期
    關(guān)鍵詞:差異質(zhì)量

    楊沖,王文穎,劉攀,周華坤,毛旭鋒

    (1.青海師范大學(xué)地理科學(xué)學(xué)院,青海 西寧 810008;2.青海師范大學(xué)生命科學(xué)學(xué)院,青海 西寧 810008;3.中國(guó)科學(xué)院西北高原生物研究所,青海 西寧 810008)

    土壤作為一種重要的自然資源,是植物生長(zhǎng)、發(fā)育的物質(zhì)基礎(chǔ),能夠控制和調(diào)節(jié)植物生長(zhǎng)的生態(tài)過程,是陸地生態(tài)系統(tǒng)的重要組成部分,也是人類賴以生存的物質(zhì)基礎(chǔ),其質(zhì)量狀況直接關(guān)系到人類的健康和社會(huì)經(jīng)濟(jì)的可持續(xù)性發(fā)展[1-4]。土壤質(zhì)量是土壤的固有屬性,指一種特定類型的土壤在自然或有管理的生態(tài)系統(tǒng)邊界內(nèi)發(fā)揮作用的能力,是土壤維持生態(tài)系統(tǒng)生物生產(chǎn)力、保護(hù)環(huán)境質(zhì)量和促進(jìn)動(dòng)植物健康能力的綜合衡量指標(biāo),能反映自然因素和人類活動(dòng)對(duì)土壤的影響[5-8]。土壤的形成、分布、人為干擾以及其所處的綜合自然環(huán)境等因素均能導(dǎo)致土壤質(zhì)量的時(shí)空分異格局[9-10]。目前,土壤質(zhì)量評(píng)價(jià)的相關(guān)研究已經(jīng)取得了豐碩成果,但仍沒有一個(gè)普適的、統(tǒng)一的評(píng)價(jià)標(biāo)準(zhǔn)[11],不同的研究區(qū)域、不同的評(píng)價(jià)目的和對(duì)土壤不同功能的側(cè)重決定了評(píng)價(jià)指標(biāo)的差異,不同的評(píng)價(jià)方法也會(huì)對(duì)土壤質(zhì)量評(píng)價(jià)結(jié)果產(chǎn)生顯著影響[12]。常用的土壤質(zhì)量評(píng)價(jià)方法包括土壤質(zhì)量卡片及監(jiān)測(cè)系統(tǒng)[13]、土壤質(zhì)量指數(shù)法[14]、主成分分析法[15]、模糊關(guān)聯(lián)法[16]、動(dòng)態(tài)土壤質(zhì)量模型[17]、管理評(píng)估法[18]等。土壤質(zhì)量受多個(gè)因素影響,且各個(gè)因素之間存在一定的相關(guān)性,致使反映土壤質(zhì)量狀況的若干指標(biāo)之間存在信息重疊[19]。主成分分析就是把多個(gè)指標(biāo)化為少數(shù)幾個(gè)綜合指標(biāo)的一種統(tǒng)計(jì)分析方法,將多個(gè)影響土壤質(zhì)量的因素進(jìn)行降維分析,提取主成分,弱化變量之間的自相關(guān)引起的誤差,被廣泛應(yīng)用在土壤質(zhì)量的評(píng)價(jià)研究中[20]。

    黃河源區(qū)是青藏高原的重要組成部分,是黃河上游主要的產(chǎn)流區(qū)、水源涵養(yǎng)區(qū)和水源補(bǔ)給區(qū),提供了約三分之一的黃河干流水量[21-22],是我國(guó)重要的生態(tài)安全屏障。黃河源區(qū)植被主要以草地為主,高寒草原、高寒草甸和沼澤草甸是黃河源區(qū)3種主要的天然草地類型,占該區(qū)可利用草地面積的80%以上[23],人工草地作為黃河源區(qū)草地恢復(fù)的一項(xiàng)重要措施,有助于提高植被生產(chǎn)力、改善土壤理化性質(zhì),提高土壤質(zhì)量,同時(shí)對(duì)畜牧業(yè)的發(fā)展和草地生態(tài)的恢復(fù)有重要作用,已成為研究高寒草地生態(tài)問題的重要組成部分[24-25]。為了保護(hù)黃河源區(qū)的生態(tài)環(huán)境以及草地資源,需要了解黃河源區(qū)主要草地生態(tài)系統(tǒng)的土壤理化性質(zhì)特征并評(píng)價(jià)其土壤質(zhì)量的狀況,進(jìn)而才能制定科學(xué)合理的高寒草地利用、保護(hù)、恢復(fù)措施。本試驗(yàn)以黃河源區(qū)不同類型高寒草地為研究對(duì)象,在總結(jié)和參考前人研究的基礎(chǔ)上,選擇12個(gè)土壤指標(biāo)(土壤含水量、pH、容重、有機(jī)碳、全碳、氨氮、硝氮、全氮、速效磷、全磷、速效鉀、全鉀)對(duì)土壤理化性質(zhì)特征進(jìn)行研究,運(yùn)用主成分分析方法對(duì)土壤質(zhì)量進(jìn)行評(píng)價(jià),以期能夠客觀、全面地反映該區(qū)土壤狀況的真實(shí)面貌,為黃河源區(qū)高寒草地的合理利用、人工草地的建植與改良以及畜牧業(yè)的可持續(xù)發(fā)展提供理論依據(jù)。

    1 材料和方法

    1.1 研究區(qū)域概況

    研究地點(diǎn)位于青海省果洛藏族自治州(圖1),地處黃河源區(qū),地理坐標(biāo)E 97°54′~101°50′,N 32°31′~35°40′,境內(nèi)平均海拔4 200 m,是“中華水塔”的重要組成部分。一年只有冷暖兩季,沒有四季之分,冷季持續(xù)時(shí)間長(zhǎng)達(dá)7~8個(gè)月,全年無絕對(duì)無霜期,年平均氣溫-4 ℃~2 ℃,年降水量235.0~974.6 mm,屬高寒半濕潤(rùn)和半干旱氣候區(qū)。境內(nèi)植被類型主要以高寒草原、高寒草甸、高寒灌叢和藏嵩草沼澤草甸為主,局部地區(qū)有少量林地、耕地和人工草地分布,畜牧業(yè)為當(dāng)?shù)刂饕a(chǎn)方式,放牧家畜主要為藏綿羊和牦牛。

    圖1 研究區(qū)示意圖Fig.1 Schematic diagram of the study area

    1.2 樣地設(shè)置

    本研究中,共設(shè)置高寒草原(G1)、高寒草甸(G2)、沼澤草甸(G3)、人工草地(G4)等4種不同類型的高寒草地(表1)。其中人工草地是高寒草甸極度退化后,人工翻耕種植垂穗披堿草的草地。播種時(shí)用磷酸二胺和羊板糞作基肥,磷酸二胺施用量150~300 kg/hm2,后期無施肥措施。人工草地建植后第1年到第2年的返青期絕對(duì)禁牧,之后作為冬季牧場(chǎng)。每個(gè)草地類型包括3個(gè)10 m×10 m的重復(fù)樣地,重復(fù)樣地彼此間盡可能使它們有相近的地形、植被和土壤類型(樣地之間空間距離1~3 km,避免假重復(fù)),共計(jì)12個(gè)樣地。G1草地在禁牧區(qū)(禁牧年限為7年),無放牧。G2、G3、G4草地為冬季牧場(chǎng)。G1草地在瑪多縣,G2,G3,G4草地在瑪沁縣。

    表1 樣地信息

    1.3 樣品采集與處理

    2017年8月中旬,在每個(gè)類型樣地采用隨機(jī)步程法設(shè)置3個(gè)樣方,對(duì)每個(gè)樣方用全球定位系統(tǒng)(GPS)進(jìn)行定位,在樣方上用土鉆采集土壤樣品,測(cè)定深度為20 cm,每10 cm為1層,重復(fù)3次,將3次重復(fù)分層混和為一個(gè)混合樣,然后將土樣分層分別用自封袋封裝帶回實(shí)驗(yàn)室,用孔徑2 mm土篩剔除植物根系和石礫等雜物,將篩分出來的土樣風(fēng)干后過1 mm和0.15 mm土篩進(jìn)一步分樣,然后進(jìn)行土壤理化性質(zhì)測(cè)定。另在該樣方做土壤剖面,用鋁土盒和容積100 cm3環(huán)刀在每層中間取樣,每個(gè)樣方3個(gè)重復(fù),做好標(biāo)記用自封袋封裝帶回實(shí)驗(yàn)室,分別用于測(cè)定土壤含水量和容重。

    1.4 分析方法

    土壤含水量和容重用烘干法測(cè)定,土壤pH用酸度計(jì)法測(cè)定,土壤全氮、全碳、有機(jī)碳(酸熏處理)通過Elementar Vario ELⅢ元素分析儀測(cè)定,土壤氨氮、硝氮、速效磷、全磷通過CleverChem Anna全自動(dòng)間斷化學(xué)分析儀測(cè)定,土壤速效鉀、全鉀通過島津A6300原子吸收儀測(cè)定。

    1.5 數(shù)據(jù)處理及評(píng)價(jià)方法

    數(shù)據(jù)記錄用Excel 2007完成,作圖用ArcGIS完成,數(shù)據(jù)統(tǒng)計(jì)分析及評(píng)價(jià)用SPSS 21完成,用平均值和標(biāo)準(zhǔn)差表示測(cè)定結(jié)果,用方差分析(ANOVA,LSD)、獨(dú)立樣本T檢驗(yàn)來比較不同處理間差異顯著性。主成分分析方法[26]:(1) KMO和Bartlett 球形度檢驗(yàn),判定選取指標(biāo)是否可以進(jìn)行主成分分析。(2)用標(biāo)準(zhǔn)化法消除不同指標(biāo)間由于量綱不同而造成的數(shù)量級(jí)上的差異。(3)提取特征值大于1,累計(jì)方差貢獻(xiàn)率大于90%的主成分。(4)主成分特征向量與相應(yīng)指標(biāo)的乘積求和,求得主成分方程。(5)將標(biāo)準(zhǔn)化后的值代入各主成分方程得到各主成分得分,主成分得分與相應(yīng)的主成分貢獻(xiàn)率的乘積求和為綜合得分。

    2 結(jié)果與分析

    2.1 不同類型高寒草地土壤理化性質(zhì)特征的變化

    2.1.1 土壤容重、pH和含水量變化 在0~10 cm土層,不同草地類型土壤容重為0.78~1.35 g/cm3,具體表現(xiàn)為G1>G4>G2>G3,G2草地與其他草地之間差異顯著(P<0.05),G3草地與其他草地之間差異顯著(P<0.05);土壤pH值在各草地為5.6~7.8,具體表現(xiàn)為G1>G4>G2>G3,各草地之間均差異顯著(P<0.05);土壤含水量在各草地的變化為7.6%~87.6%,具體表現(xiàn)為G3>G2>G4>G1,各草地之間均差異顯著(P<0.05)。

    在10~20 cm土層,不同草地類型土壤容重0.93~1.57 g/cm3,具體表現(xiàn)為G1>G2>G4>G3,G1草地與其他草地之間差異顯著(P<0.05),G3草地與其他草地之間差異顯著(P<0.05);土壤pH值在各草地為5.8~7.9,具體表現(xiàn)為G1>G4>G2>G3,G1草地與其他草地之間差異顯著(P<0.05),G4草地與其他草地之間差異顯著(P<0.05);土壤含水量在各草地為9.6%~77.7%,具體表現(xiàn)為G3>G2>G4>G1,G2草地與其他草地之間差異顯著(P<0.05),G3草地與其他草地之間差異顯著(P<0.05)。

    在0~20 cm土層,土壤容重在各草地為0.86~1.46 g/cm3,具體表現(xiàn)為G1>G4>G2>G3,G1草地與其他草地之間差異顯著(P<0.05),G3草地與其他草地之間差異顯著(P<0.05);土壤pH值在各草地為5.7~7.85,具體表現(xiàn)為G1>G4>G2>G3,各草地之間均差異顯著(P<0.05);土壤含水量在各草地為8.6%~82.65%,具體表現(xiàn)為G3>G2>G4>G1,各草地之間均差異顯著(P<0.05)。

    在同一草地的不同土壤深度,土壤容重隨土壤深度的增加呈增大趨勢(shì),其中在G1、G2草地,10~20 cm土層的土壤容重顯著高于0~10 cm土層(P<0.05);土壤pH值隨土壤深度的增加呈增大趨勢(shì),但在各草地不同土壤深度均差異不顯著;土壤含水量在G1草地隨土壤深度的增加呈增大趨勢(shì),在其他草地隨土壤深度的增加呈減小趨勢(shì),其中在G2草地0~10 cm土層的土壤含水量顯著高于10~20 cm土層(P<0.05),其他草地差異不顯著(圖2)。

    圖2 土壤容重、pH和含水量Fig.2 Soil bulk density,pH and water contentin soil

    2.1.2 土壤有機(jī)碳、全碳含量變化 在0~10 cm土層,土壤有機(jī)碳含量在各草地為12.06~161.53 g/kg,在10~20 cm土層,土壤有機(jī)碳含量在各草地為8~131.76 g/kg,在0~20 cm土層,土壤有機(jī)碳含量在各草地為10.03~146.65 g/kg,具體表現(xiàn)均為G3>G2>G4>G1,其中G1與G4草地之間差異不顯著,G2與其他草地之間差異顯著(P<0.05),G3與其他草地之間差異顯著(P<0.05)。在同一草地的不同土壤深度,土壤有機(jī)碳和全碳含量隨著土壤深度的增加有減少的趨勢(shì),各草地0~10 cm土層和10~20 cm土層的土壤有機(jī)碳、全碳含量均差異不顯著(圖3)。

    圖3 土壤有機(jī)碳、全碳含量Fig.3 Contents oforganic carbon and total carbonin soil

    2.1.3 土壤氨氮、硝氮和全氮含量變化 在0~10 cm土層,土壤氨氮含量在各草地為5.32~50.10 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G1草地與其他草地之間差異均顯著(P<0.05),G4草地與其他草地之間差異均顯著(P<0.05);土壤硝氮含量在各草地為5.31~15.71 mg/kg,具體表現(xiàn)為G4>G3>G1>G2,G2與G3草地之間差異顯著(P<0.05),G4與其他草地之間差異均顯著(P<0.05);土壤全氮含量在各草地為1.38~11.51 g/kg,具體表現(xiàn)為G3>G2>G4>G1,G1與G2草地之間差異顯著(P<0.05);G3與其他草地之間均差異顯著(P<0.05)(圖4)。

    在10~20 cm土層,土壤氨氮含量在各草地為3.33~41.55 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G1草地與其他草地之間差異均顯著(P<0.05),G4草地與其他草地之間差異均顯著(P<0.05);土壤硝氮含量在各草地為3.42~13.22 mg/kg,具體表現(xiàn)為G4>G3>G1>G2,G2與G3間差異顯著(P<0.05),G4與其他草地之間差異均顯著(P<0.05);土壤全氮含量在各草地為1.33~8.58 g/kg,具體表現(xiàn)為G3>G2>G4>G1,G1與G2草地之間差異顯著(P<0.05);G3與其他草地之間均差異顯著(P<0.05)(圖4)。

    在0~20 cm土層,土壤氨氮含量在各草地為4.33~45.74 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G1草地與其他草地之間差異均顯著(P<0.05),G4草地與其他草地之間差異均顯著(P<0.05);土壤硝氮含量在各草地為4.37~14.47 mg/kg,具體表現(xiàn)為G4>G3>G1>G2,G2與G3草地之間差異顯著(P<0.05),G4與其他草地之間均差異顯著(P<0.05);土壤全氮含量在各草地為1.36~10.05 g/kg,具體表現(xiàn)為G3>G2>G4>G1,各草地之間均差異顯著(P<0.05)。

    圖4 土壤氨氮、硝氮、全氮含量Fig.4 Contents of ammonia nitrogen,nitrate nitrogen and total nitrogen in soil

    在同一草地的不同土壤深度,土壤氨氮、硝氮和全氮含量隨著土壤深度的增加均有減少的趨勢(shì),其中G1草地土壤氨氮含量在0~10 cm土層和10~20 cm土層之間差異顯著(P<0.05);G2草地土壤全氮含量在0~10 cm土層和10~20 cm土層之間差異顯著(P<0.05);G3草地土壤氨氮和硝氮含量在0~10 cm土層和10~20 cm土層之間差異顯著(P<0.05)。

    2.1.4 土壤速效磷、全磷含量變化 在0~10 cm土層,土壤速效磷含量在各草地為4.73~11.25 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G1、G2草地與G3、G4草地之間差異顯著(P<0.05);土壤全磷含量在各草地為0.43~1.72 g/kg,具體表現(xiàn)為G3>G2>G4>G1,各草地之間差異顯著(P<0.05)(圖5)。

    在10~20 cm土層,土壤速效磷含量在各草地為2.41~6.79 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G3與G4草地之間差異顯著(P<0.05),G1與其他草地之間差異均顯著(P<0.05);土壤全磷含量在各草地為0.45~1.49 g/kg,具體表現(xiàn)為G3>G2>G4>G1,G1、G2草地分別與G3、G4草地之間差異顯著(P<0.05)(圖5)。

    在0~20 cm土層,土壤速效磷含量在各草地為3.57~9.02 mg/kg,具體表現(xiàn)為G3>G2>G4>G1,G1與其他草地之間差異均顯著(P<0.05),G4與其他草地之間差異均顯著(P<0.05);土壤全磷含量在各草地的變化為0.44~1.61 g/kg,具體表現(xiàn)為G3>G2>G4>G1,各草地之間差異顯著(P<0.05)。

    在同一草地的不同土壤深度,土壤速效磷和全磷含量隨著土壤深度的增加均有減少的趨勢(shì),其中G1、G2、G3草地土壤速效磷含量在0~10 cm土層和10~20 cm土層之間差異顯著(P<0.05);各草地土壤全磷含量在0~10 cm土層和10~20 cm土層之間差異不顯著(圖5)。

    圖5 土壤速效磷、全磷含量Fig.5 Contents ofavailable phosphorus and total phosphorusin soil

    2.1.5 土壤速效鉀、全鉀含量 在0~10 cm土層,土壤速效鉀含量在各草地的變化為101.15~216.32 mg/kg,具體表現(xiàn)為G4>G2>G1>G3,其中G3與其他草地之間差異均顯著(P<0.05);土壤全鉀含量在各草地的變化為8.23~17.25 g/kg,具體表現(xiàn)為G4>G1>G2>G3,G2與G4間差異顯著(P<0.05),G3與其他草地間均差異顯著(P<0.05)(圖6)。

    在10~20 cm土層,土壤速效鉀含量在各草地的變化為82.45~154.96 mg/kg,具體表現(xiàn)為G4>G1>G2>G3,G3與G4草地之間差異顯著(P<0.05);土壤全鉀含量在各草地的變化為7.24~18.76 g/kg,具體表現(xiàn)為G4>G2>G1>G3,G3與其他草地之間差異均顯著(P<0.05),G4與其他草地之間均差異顯著(P<0.05)。

    在0~20 cm土層,土壤速效鉀含量在各草地為91.8~185.64 mg/kg,土壤全鉀含量各草地為7.74~18.01 g/kg,兩者變化趨勢(shì)相同,具體表現(xiàn)為G4>G1>G2>G3,G3與其他草地之間差異均顯著(P<0.05),G4與其他草地之間差異均顯著(P<0.05)。

    在同一草地的不同土壤深度,土壤速效鉀含量隨著土壤深度的增加均有減少的趨勢(shì),土壤全鉀含量隨著土壤深度的增加變化趨勢(shì)不明顯;其中G1和G4草地土壤速效鉀含量在0~10 cm土層和10~20 cm土層之間差異顯著(P<0.05);各草地土壤全鉀含量在0~10 cm土層和10~20 cm土層之間差異均不顯著(圖6)。

    圖6 速效鉀、全鉀含量Fig.6 Contents ofavailable potassium and total potassium in soil

    2.2 不同高寒草地0~20 cm土壤理化成分之間的相關(guān)性分析

    土壤pH值與土壤容重呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤含水量、有機(jī)碳、全碳、氨氮、全氮、速效磷、全磷含量之間呈極顯著負(fù)相關(guān)關(guān)系(P<0.01),與土壤全鉀含量之間呈顯著正相關(guān)關(guān)系(P<0.05);土壤容重與土壤含水量、有機(jī)碳、全碳、氨氮、全氮、速效磷、全磷含量之間呈極顯著負(fù)相關(guān)關(guān)系(P<0.01),與土壤速效鉀含量呈顯著正相關(guān)關(guān)系(P<0.05),與土壤全鉀含量呈極顯著正相關(guān)關(guān)系(P<0.01);土壤含水量與土壤有機(jī)碳、全碳、氨氮、全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤速效鉀含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05),與土壤全鉀含量呈極顯著負(fù)相關(guān)關(guān)系(P<0.01);土壤有機(jī)碳含量與土壤全碳、氨氮、全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤速效鉀含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05),與土壤全鉀含量呈極顯著負(fù)相關(guān)關(guān)系(P<0.01);土壤全碳含量與土壤氨氮、全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤速效鉀含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05),與土壤全鉀含量呈極顯著負(fù)相關(guān)關(guān)系(P<0.01);土壤氨氮含量與土壤全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01);土壤硝氮含量與其他土壤理化成分之間相關(guān)性均不顯著;土壤全氮含量與土壤速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤速效鉀含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05),與土壤全鉀含量呈極顯著負(fù)相關(guān)關(guān)系(P<0.01);土壤速效磷含量與土壤全磷含量呈極顯著正相關(guān)關(guān)系(P<0.01),土壤速效鉀含量與土壤全鉀含量呈極顯著正相關(guān)關(guān)系(P<0.01);土壤全磷含量與土壤全鉀含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05)(表2)。

    表2 土壤理化指標(biāo)之間的相關(guān)性

    2.3 基于主成分分析的不同高寒草地0~20 cm土壤質(zhì)量綜合評(píng)價(jià)

    2.3.1 KMO和Bartlett球形度檢驗(yàn) KMO檢驗(yàn)用于檢查變量間的相關(guān)性和偏相關(guān)性,Bartlett球形度檢驗(yàn)用于檢驗(yàn)相關(guān)陣中各變量間的相關(guān)性。KMO統(tǒng)計(jì)量的取值在0和1之間,當(dāng)所有變量間的簡(jiǎn)單相關(guān)系數(shù)平方和遠(yuǎn)遠(yuǎn)大于偏相關(guān)系數(shù)平方和時(shí),KMO值接近1,KMO值越接近于1,意味著變量間的相關(guān)性越強(qiáng),原有變量越適合作因子分析;當(dāng)所有變量間的簡(jiǎn)單相關(guān)系數(shù)平方和接近0時(shí),KMO值接近0,KMO值越接近于0,意味著變量間的相關(guān)性越弱,原有變量越不適合作因子分析。KMO度量標(biāo)準(zhǔn):0.9以上表示非常適合;0.8表示適合;0.7表示一般;0.6表示不太適合;0.5以下表示極不適合。將土壤含水量、pH值、容重、土壤有機(jī)碳、全碳、氨氮、硝氮、全氮、速效磷、全磷、速效鉀、全鉀含量等12個(gè)土壤指標(biāo)經(jīng)KMO和Bartlett 檢驗(yàn),KMO值為0.810,Bartlett的球形度檢驗(yàn)的相伴概率P<0.01(極顯著水平),說明本研究選取的土壤質(zhì)量評(píng)價(jià)指標(biāo)間存在較強(qiáng)的相關(guān)性,采用主成分分析法來評(píng)價(jià)各樣地的土壤質(zhì)量狀況是可行的(表3)。

    表3 KMO和Bartlett 球形度檢驗(yàn)

    2.3.2 計(jì)算貢獻(xiàn)率和主成分提取 對(duì)主成分進(jìn)行提取時(shí),依據(jù)主成分特征值大于1的原則。從表4結(jié)果共提取2個(gè)主成分,第1主成分對(duì)總方差的貢獻(xiàn)率為78.037%,第2主成分對(duì)總方差的貢獻(xiàn)率為13.783%,累積貢獻(xiàn)率為91.820%,說明這2個(gè)主成分代表了該試驗(yàn)中原始數(shù)據(jù)91.820%的信息。

    表4 方差分解主成分提取分析

    2.3.3 計(jì)算綜合得分并排序 應(yīng)用主成分分析法對(duì)不同高寒草地土壤質(zhì)量進(jìn)行綜合評(píng)價(jià),就是利用主成分綜合得分的大小來評(píng)價(jià)土壤質(zhì)量的高低,主成分綜合得分越大土壤質(zhì)量越高,反之則越低。在0~10、10~20以及0~20 cm土層,土壤質(zhì)量高低排序均為G3>G2>G4>G1;各草地0~10 cm土層土壤質(zhì)量均高于10~20和0~20 cm土層,0~20 cm土層土壤質(zhì)量高于10~20 cm土層(表5)。

    表5 不同高寒草地土壤質(zhì)量綜合得分及排序

    3 討論

    本研究中,4種不同類型高寒草地的土壤理化性狀存在差異。其中,沼澤草甸的土壤容重最小,可能的原因是藏嵩草沼澤化草甸土壤環(huán)境有明顯的高濕和低溫特征,且根系發(fā)達(dá),土壤中存在大量未分解或半分解的有機(jī)殘?bào)w,使土壤容重變小[27-28]。沼澤草甸和高寒草甸的土壤pH值比較低,可能是因?yàn)橥寥赖膒H值與土壤含水量之間存在較強(qiáng)的負(fù)相關(guān)關(guān)系[29],而沼澤草甸和高寒草甸有著較高的土壤含水量。沼澤草甸和高寒草甸兩處草地的土壤含水量較高,可能是因?yàn)閮商帢拥匚挥诤庸鹊屯萏?,降水較多,并且均有河流經(jīng)過,且兩處草地植物根系發(fā)達(dá),使得土壤孔隙度增大,土壤毛管持水量上升,保水力增加,另外,兩處草地植被茂盛,一定程度上阻止了土壤表層水分的蒸發(fā)。沼澤草甸土壤有機(jī)碳和全碳含量最高,可能是因?yàn)檎訚赏了诌^多,濕生植物生長(zhǎng)旺盛,植物根系十分發(fā)達(dá),根系脫落物和溢泌物較多,且土壤處于缺氧狀態(tài),有機(jī)物質(zhì)分解緩慢,土壤有機(jī)碳含量較高[30];高寒草原的土壤有機(jī)碳含量較低,一方面可能是因?yàn)楦吆菰寥赖纳郴潭容^高,土壤通透性及含氧量相對(duì)較高,土壤呼吸作用強(qiáng)烈,導(dǎo)致土壤有機(jī)碳有降低的趨勢(shì)[31],另一方面土壤有機(jī)質(zhì)主要來源于植物地上部分的凋落物及地下的根系,植物每年都有大量有機(jī)物質(zhì)進(jìn)入土壤,并通過微生物的活動(dòng),逐漸變?yōu)橥寥烙袡C(jī)碳,而高寒草原植物有著相對(duì)較低的地上和地下生物量。沼澤草甸土壤全氮含量顯著高于其他草地,可能是因?yàn)橥寥廊?5%來源于土壤有機(jī)質(zhì)[32],而沼澤草甸的土壤有機(jī)質(zhì)含量最高。土壤速效磷、全磷含量在各草地均呈現(xiàn)出普遍匱乏狀態(tài),這可能是因?yàn)辄S河源區(qū)地勢(shì)高、溫度較低,較低的氣溫抑制了有機(jī)物分解和礦化,造成土壤磷含量較低;沼澤草甸由于水分大,有助于速效磷的形成,其速效磷含量最高,這與慕軍鵬[33]在若爾蓋的研究結(jié)果一致。沼澤草甸土壤速效鉀和全鉀含量較低,可能是因?yàn)殁浕衔镆话愣疾粨]發(fā),但在水中有較高的溶解度,它的代換量比磷大,容易從土壤膠體上代換出來[34],而沼澤草甸的土壤水分過多,容易造成鉀元素淋失,本研究結(jié)果與劉玉萍[35]在青海湖區(qū)的研究結(jié)果一致;此外,NH4+和K+具有相似的化學(xué)性質(zhì)、電負(fù)性和相近的離子半徑,在同時(shí)存在這兩種離子的環(huán)境中,K+和NH4+的吸收存在拮抗作用,NH4+的存在會(huì)抑制K+的釋放,常觀察到NH4+和K+在晶格中會(huì)共同占據(jù)某一原子位置并形成相互取代的現(xiàn)象[36],而沼澤草甸的氨氮含量較高,可能也是沼澤草甸速效鉀和全鉀含量較低的原因。

    本研究中,高寒草地的不同土層之間土壤理化性狀存在差異。其中,土壤容重隨土壤深度的增加呈增大趨勢(shì),這與前人[37-38]的研究結(jié)果一致,可能的原因是由于植被對(duì)土壤的保護(hù)作用,使得土壤表層受水力和風(fēng)力及凍融作用的侵蝕作用較小,且植被的增加相應(yīng)增加了凋落物的輸入,使得土壤表層的有機(jī)質(zhì)和其他植物殘?bào)w等腐殖質(zhì)含量比較多,土壤容重較小。同一樣地不同土層間pH值沒有顯著差異,可能是由于土壤的緩沖特性使得土壤pH在一定的土壤深度上總體上保持穩(wěn)定[39]。高寒草原土壤含水量隨土壤深度的增加呈增加趨勢(shì),在其他高寒草地隨土壤深度的增加呈降低趨勢(shì),可能的原因是土壤含水量不僅與降水關(guān)系密切,而且還受蒸發(fā)、太陽輻射等因素的制約;在一定的時(shí)期內(nèi),高寒草原的降水較少,蒸發(fā)旺盛,導(dǎo)致土壤上層含水量低于土壤下層;沼澤草甸和高寒草甸兩處草地的土壤含水量較高,可能是因?yàn)閮商帢拥匚挥诤庸鹊屯萏?,降水較多,并且均有河流經(jīng)過,且兩處草地植物根系發(fā)達(dá),使得土壤孔隙度增大,土壤毛管持水量上升,保水力增加,另外,兩處草地植被茂盛,一定程度上阻止了土壤表層水分的蒸發(fā)。土壤氨氮、硝氮、全氮含量在不同土層間的變化趨勢(shì)一致,均隨著土層深度的增加呈減少趨勢(shì),原因可能是土壤淋溶過程使氮素在土壤不同深度含量不同,根系和地表凋落物中碳氮比和微生物區(qū)系也會(huì)使土壤氮分解不同[40-41]。土壤速效磷、全磷、速效鉀、全鉀含量均隨著土壤深度的增加呈降低趨勢(shì),但速效磷和速效鉀含量在各土層之間變化較大,全磷和全鉀含量在不同土層之間變化不大,相對(duì)比較穩(wěn)定,可能是因?yàn)槿缀腿浐恐饕c成土母質(zhì)關(guān)系較大。

    相關(guān)性分析的結(jié)果表明,土壤有機(jī)碳含量與土壤全碳、氨氮、全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),土壤全碳含量與土壤氨氮、全氮、速效磷、全磷含量之間呈極顯著正相關(guān)關(guān)系(P<0.01),此結(jié)果與已有的研究結(jié)果[42-43]相一致,土壤碳氮通常高度耦合,另外,土壤中的粘土礦物、碳酸鹽和有機(jī)質(zhì)能緩慢固定來自基巖的磷,反過來,較高的土壤磷含量可進(jìn)一步加強(qiáng)固氮作用和有機(jī)質(zhì)的積累過程[44],最終導(dǎo)致三者之間的強(qiáng)相關(guān)關(guān)系。

    土壤質(zhì)量狀況數(shù)值化的綜合評(píng)價(jià),能較好地反映土壤質(zhì)量的實(shí)際情況。土壤質(zhì)量的評(píng)價(jià)體系和方法很難建立統(tǒng)一的標(biāo)準(zhǔn),研究區(qū)域和尺度的不同,側(cè)重的土壤功能不同,決定了選取評(píng)價(jià)指標(biāo)和方法的差異。本研究通過主成分分析對(duì)不同類型高寒草地土壤質(zhì)量進(jìn)行了評(píng)價(jià),選取的12個(gè)土壤指標(biāo)經(jīng)KMO和Bartlett 檢驗(yàn),KMO值為0.810,Bartlett的球形度檢驗(yàn)的相伴概率P<0.01,說明本研究選取的土壤質(zhì)量評(píng)價(jià)指標(biāo)間存在較強(qiáng)的相關(guān)性,共提取2個(gè)主成分,這2個(gè)主成分代表了該試驗(yàn)中原始數(shù)據(jù)91.820%的信息,說明采用主成分分析法來評(píng)價(jià)各草地的土壤質(zhì)量狀況是準(zhǔn)確可行的。評(píng)價(jià)結(jié)果顯示,土壤質(zhì)量高低依次為沼澤草甸>高寒草甸>人工草地>高寒草原,各草地 0~10 cm土層土壤質(zhì)量均高于10~20和0~20 cm土層,此結(jié)果與前人的研究結(jié)果[45-47]相一致。高寒草原的土壤質(zhì)量最差,可能是因?yàn)楦吆菰寥赖耐翆虞^薄,一般在10~25 cm左右,缺少淀積層,發(fā)育程度弱,土壤受水分限制嚴(yán)重并且有沙化趨勢(shì),植被稀疏且根系不發(fā)達(dá),土壤有機(jī)質(zhì)和其他成分含量較少,導(dǎo)致土壤養(yǎng)分較低,土壤質(zhì)量差。高寒草甸土壤質(zhì)量較好,可能是因?yàn)橄啾扔诟吆菰?,高寒草甸土壤A層會(huì)形成 10~15 cm的腐殖層,有利于營(yíng)養(yǎng)物質(zhì)的形成和轉(zhuǎn)換,提高了土壤質(zhì)量。沼澤草甸土壤質(zhì)量最好,可能是因?yàn)檎訚刹莸橹脖簧a(chǎn)力較高,根系發(fā)達(dá),有機(jī)物質(zhì)輸入多,會(huì)形成一層厚厚的毛氈狀腐質(zhì)化土層,土壤礦化能力較強(qiáng),土壤養(yǎng)分含量較高。退化高寒草甸建植人工草地土壤質(zhì)量低于高寒草甸,可能是因?yàn)楦吆莸厣鷳B(tài)系統(tǒng)是一個(gè)脆弱的生態(tài)系統(tǒng),發(fā)生嚴(yán)重退化后,建植人工草地盡管短期內(nèi)可有效恢復(fù)地上植被,但土壤質(zhì)量的恢復(fù)則是一個(gè)長(zhǎng)期的過程。各草地0~10 cm土層土壤質(zhì)量均高于10~20 cm土層,即表層土壤質(zhì)量最好,這與表層土壤碳氮磷鉀總含量和速效養(yǎng)分含量高密切相關(guān),也與植被養(yǎng)分輸入多密切相關(guān),另外也說明一旦表層剝蝕,土壤質(zhì)量將嚴(yán)重下降,因此要預(yù)防水土流失以及過載過牧對(duì)表層土壤的破壞。

    4 結(jié)論

    在0~20 cm土層土壤特征方面,土壤含水量、土壤有機(jī)碳、全碳、氨氮、全氮、速效磷以及全磷含量高度耦合,具體表現(xiàn)均為沼澤草甸>高寒草甸>人工草地>高寒草原;土壤硝氮含量表現(xiàn)為人工草地>沼澤草甸>高寒草原>高寒草甸;土壤pH值與土壤容重呈極顯著正相關(guān)關(guān)系,具體表現(xiàn)均為高寒草原>人工草地>高寒草甸>沼澤草甸;土壤速效鉀與全鉀含量表現(xiàn)均為人工草地>高寒草原>高寒草甸>沼澤草甸。

    基于主成分分析的評(píng)價(jià)結(jié)果顯示,在0~10、10~20及0~20 cm土層,土壤質(zhì)量高低排序均為沼澤草甸>高寒草甸>人工草地>高寒草原;各草地均為0~10 cm土層土壤質(zhì)量最高。

    猜你喜歡
    差異質(zhì)量
    相似與差異
    音樂探索(2022年2期)2022-05-30 21:01:37
    “質(zhì)量”知識(shí)鞏固
    質(zhì)量守恒定律考什么
    找句子差異
    做夢(mèng)導(dǎo)致睡眠質(zhì)量差嗎
    DL/T 868—2014與NB/T 47014—2011主要差異比較與分析
    生物為什么會(huì)有差異?
    關(guān)于質(zhì)量的快速Q(mào)&A
    質(zhì)量投訴超六成
    汽車觀察(2016年3期)2016-02-28 13:16:26
    M1型、M2型巨噬細(xì)胞及腫瘤相關(guān)巨噬細(xì)胞中miR-146a表達(dá)的差異
    亚洲激情五月婷婷啪啪| 国产熟女欧美一区二区| 97精品久久久久久久久久精品| 少妇的逼好多水| 99久久中文字幕三级久久日本| 久久久亚洲精品成人影院| 欧美激情国产日韩精品一区| 亚洲天堂国产精品一区在线| 麻豆成人午夜福利视频| 亚洲国产欧美在线一区| 只有这里有精品99| 成人毛片a级毛片在线播放| videos熟女内射| 性色avwww在线观看| 亚洲内射少妇av| 精华霜和精华液先用哪个| 国产午夜福利久久久久久| 久久久精品免费免费高清| 亚洲天堂国产精品一区在线| 国产中年淑女户外野战色| 啦啦啦中文免费视频观看日本| 亚洲最大成人中文| 国产精品国产三级国产专区5o| 高清欧美精品videossex| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美 国产精品| 国产精品99久久久久久久久| 欧美成人精品欧美一级黄| 99久国产av精品国产电影| 激情五月婷婷亚洲| 黄片无遮挡物在线观看| 99视频精品全部免费 在线| 日本熟妇午夜| 国产毛片在线视频| 大又大粗又爽又黄少妇毛片口| 新久久久久国产一级毛片| 欧美丝袜亚洲另类| 女人十人毛片免费观看3o分钟| 嫩草影院新地址| 欧美日韩视频高清一区二区三区二| 国产一区亚洲一区在线观看| 亚洲精品亚洲一区二区| 黄色一级大片看看| 国产亚洲av片在线观看秒播厂| 少妇熟女欧美另类| 在线免费观看不下载黄p国产| 亚洲自偷自拍三级| 久久久久精品久久久久真实原创| 国产免费一级a男人的天堂| 少妇高潮的动态图| 秋霞伦理黄片| av线在线观看网站| 亚洲人与动物交配视频| 国模一区二区三区四区视频| 国产高清国产精品国产三级 | 国产高清三级在线| 久久精品国产亚洲av天美| 特大巨黑吊av在线直播| 建设人人有责人人尽责人人享有的 | 97精品久久久久久久久久精品| 国产精品久久久久久精品电影小说 | 黄片无遮挡物在线观看| 2022亚洲国产成人精品| 日本黄大片高清| 精品国产一区二区三区久久久樱花 | 黄片无遮挡物在线观看| 22中文网久久字幕| 亚洲精品,欧美精品| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩一区二区视频在线观看视频在线 | 国产精品一二三区在线看| 久久国内精品自在自线图片| 青青草视频在线视频观看| 国产男人的电影天堂91| 中文欧美无线码| 嘟嘟电影网在线观看| 成人漫画全彩无遮挡| 一级毛片久久久久久久久女| eeuss影院久久| 精品人妻偷拍中文字幕| av福利片在线观看| 中文资源天堂在线| 在线播放无遮挡| 赤兔流量卡办理| 国产成人午夜福利电影在线观看| 成年女人在线观看亚洲视频 | 亚洲图色成人| 毛片一级片免费看久久久久| 高清av免费在线| 国产精品偷伦视频观看了| 亚洲不卡免费看| 欧美xxxx性猛交bbbb| 精品一区二区三区视频在线| 久久影院123| 久久人人爽人人爽人人片va| 欧美日韩在线观看h| 国产白丝娇喘喷水9色精品| 熟女av电影| 亚洲三级黄色毛片| 国产精品一区www在线观看| 日韩欧美精品免费久久| 国产黄片视频在线免费观看| 亚洲精品久久午夜乱码| 国产高清国产精品国产三级 | 男的添女的下面高潮视频| 十八禁网站网址无遮挡 | 国产亚洲91精品色在线| 国产一区二区三区av在线| 国产精品一区二区性色av| 天堂中文最新版在线下载 | 97人妻精品一区二区三区麻豆| 日日啪夜夜撸| 我的老师免费观看完整版| 欧美精品一区二区大全| 亚洲av二区三区四区| .国产精品久久| 久久久久久久久久久免费av| 国产高清三级在线| 天堂网av新在线| 亚洲欧美日韩另类电影网站 | 国产伦在线观看视频一区| 不卡视频在线观看欧美| 欧美日本视频| 久久亚洲国产成人精品v| 另类亚洲欧美激情| 精品久久久久久久人妻蜜臀av| a级一级毛片免费在线观看| 亚洲人成网站高清观看| 新久久久久国产一级毛片| 熟妇人妻不卡中文字幕| 亚洲精品国产av蜜桃| 嘟嘟电影网在线观看| 亚洲美女搞黄在线观看| 极品少妇高潮喷水抽搐| 亚洲av成人精品一区久久| 免费人成在线观看视频色| 80岁老熟妇乱子伦牲交| 免费大片黄手机在线观看| 亚洲色图av天堂| 精品国产一区二区三区久久久樱花 | 久久99热这里只有精品18| 自拍欧美九色日韩亚洲蝌蚪91 | 99热全是精品| 精品久久久久久久末码| 亚洲综合精品二区| 99久国产av精品国产电影| 男人和女人高潮做爰伦理| 国产人妻一区二区三区在| 老司机影院成人| 国模一区二区三区四区视频| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美精品免费久久| 亚洲成人久久爱视频| 亚洲成人久久爱视频| 国产69精品久久久久777片| 久热久热在线精品观看| 欧美激情在线99| 美女cb高潮喷水在线观看| 免费在线观看成人毛片| 亚洲人与动物交配视频| 欧美日韩亚洲高清精品| 一级黄片播放器| 国产美女午夜福利| 久久久久国产精品人妻一区二区| 在线免费观看不下载黄p国产| 性色av一级| 亚洲精品中文字幕在线视频 | 中文在线观看免费www的网站| 高清毛片免费看| 赤兔流量卡办理| 少妇高潮的动态图| 纵有疾风起免费观看全集完整版| 国产精品偷伦视频观看了| 在线观看一区二区三区| 国产国拍精品亚洲av在线观看| 欧美少妇被猛烈插入视频| 欧美少妇被猛烈插入视频| 久久精品久久精品一区二区三区| 日本欧美国产在线视频| 天堂俺去俺来也www色官网| 亚洲成人一二三区av| 国产一区有黄有色的免费视频| 男人添女人高潮全过程视频| 免费少妇av软件| 成年女人在线观看亚洲视频 | av免费在线看不卡| 亚洲自偷自拍三级| 日本爱情动作片www.在线观看| 国精品久久久久久国模美| 欧美成人一区二区免费高清观看| 免费观看的影片在线观看| 亚洲欧美一区二区三区国产| 超碰97精品在线观看| 国产91av在线免费观看| 亚洲av成人精品一区久久| 亚洲精华国产精华液的使用体验| 不卡视频在线观看欧美| 国产成人精品久久久久久| 亚洲精品日韩av片在线观看| 爱豆传媒免费全集在线观看| 99九九线精品视频在线观看视频| 欧美成人a在线观看| 国产淫片久久久久久久久| 亚洲欧美日韩另类电影网站 | 在线观看免费高清a一片| 少妇人妻精品综合一区二区| 免费看日本二区| 欧美一级a爱片免费观看看| 中文字幕久久专区| 一个人看视频在线观看www免费| 久久亚洲国产成人精品v| 国产成人精品久久久久久| 成人毛片60女人毛片免费| 欧美精品人与动牲交sv欧美| 亚洲精品自拍成人| 3wmmmm亚洲av在线观看| 国产 一区精品| 身体一侧抽搐| 夫妻性生交免费视频一级片| 久久女婷五月综合色啪小说 | 亚洲成人精品中文字幕电影| 大码成人一级视频| av卡一久久| 69av精品久久久久久| 亚洲最大成人手机在线| 久久精品国产亚洲网站| 久久久a久久爽久久v久久| 国产成人午夜福利电影在线观看| 欧美激情久久久久久爽电影| 三级国产精品欧美在线观看| 黄片无遮挡物在线观看| 久久久久久伊人网av| 日韩强制内射视频| 色播亚洲综合网| av女优亚洲男人天堂| 涩涩av久久男人的天堂| 在线播放无遮挡| 国产有黄有色有爽视频| 蜜臀久久99精品久久宅男| 亚洲欧美成人综合另类久久久| 高清av免费在线| 最近手机中文字幕大全| 国产精品久久久久久精品古装| 成人一区二区视频在线观看| 久久久亚洲精品成人影院| 综合色av麻豆| 亚洲精品第二区| 观看免费一级毛片| 国产精品一区二区三区四区免费观看| 日韩人妻高清精品专区| 亚洲精品日韩av片在线观看| 一级a做视频免费观看| 天天躁日日操中文字幕| 国产亚洲精品久久久com| 女的被弄到高潮叫床怎么办| 久久久精品免费免费高清| 国产精品一区二区在线观看99| 自拍偷自拍亚洲精品老妇| 中文精品一卡2卡3卡4更新| 啦啦啦中文免费视频观看日本| 午夜激情福利司机影院| 日日啪夜夜爽| 十八禁网站网址无遮挡 | 最后的刺客免费高清国语| 春色校园在线视频观看| 女人被狂操c到高潮| 在线观看一区二区三区激情| 日韩国内少妇激情av| 天天躁夜夜躁狠狠久久av| 成人国产av品久久久| 亚洲怡红院男人天堂| 国国产精品蜜臀av免费| 美女内射精品一级片tv| 日本黄色片子视频| 午夜视频国产福利| 国内精品美女久久久久久| 少妇人妻 视频| 少妇裸体淫交视频免费看高清| 国产精品一区www在线观看| 最近最新中文字幕大全电影3| 白带黄色成豆腐渣| 免费观看a级毛片全部| 久久久久久久久久成人| 不卡视频在线观看欧美| 日日啪夜夜撸| 九色成人免费人妻av| av国产久精品久网站免费入址| 九九爱精品视频在线观看| 欧美+日韩+精品| 亚洲精品,欧美精品| 国产一区二区亚洲精品在线观看| 高清欧美精品videossex| 欧美高清性xxxxhd video| 国产成人freesex在线| 欧美日韩在线观看h| 噜噜噜噜噜久久久久久91| 2021少妇久久久久久久久久久| 国产在线一区二区三区精| 一个人看的www免费观看视频| 久久精品久久久久久噜噜老黄| 一区二区三区四区激情视频| 国产亚洲一区二区精品| 成人午夜精彩视频在线观看| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 国产中年淑女户外野战色| 亚洲丝袜综合中文字幕| 三级国产精品欧美在线观看| 日本wwww免费看| 亚洲高清免费不卡视频| 亚洲精品一区蜜桃| 国产成人精品福利久久| 成人黄色视频免费在线看| 男女国产视频网站| 在现免费观看毛片| 久热这里只有精品99| 成人毛片60女人毛片免费| 欧美另类一区| 日韩亚洲欧美综合| 免费黄频网站在线观看国产| 色网站视频免费| 国产亚洲91精品色在线| 日韩欧美精品免费久久| 成年免费大片在线观看| 熟妇人妻不卡中文字幕| 青春草视频在线免费观看| 小蜜桃在线观看免费完整版高清| 边亲边吃奶的免费视频| 亚洲自偷自拍三级| 亚洲欧美清纯卡通| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产日韩一区二区| av在线观看视频网站免费| 国国产精品蜜臀av免费| 国产老妇伦熟女老妇高清| av在线老鸭窝| 99视频精品全部免费 在线| 亚洲欧美成人综合另类久久久| 青春草亚洲视频在线观看| 国产亚洲一区二区精品| 白带黄色成豆腐渣| 亚洲美女视频黄频| 插阴视频在线观看视频| 国产亚洲av片在线观看秒播厂| 精品熟女少妇av免费看| 亚洲内射少妇av| 51国产日韩欧美| 99re6热这里在线精品视频| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 日日啪夜夜爽| 麻豆久久精品国产亚洲av| 自拍欧美九色日韩亚洲蝌蚪91 | 在线观看人妻少妇| 亚洲精品亚洲一区二区| 日韩成人av中文字幕在线观看| 熟妇人妻不卡中文字幕| 欧美区成人在线视频| 精品一区在线观看国产| 国产成人aa在线观看| 国产综合懂色| 制服丝袜香蕉在线| 麻豆国产97在线/欧美| 狂野欧美白嫩少妇大欣赏| 国产毛片在线视频| 国产一区二区在线观看日韩| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 久久久久国产网址| 只有这里有精品99| 亚洲av一区综合| 亚洲国产成人一精品久久久| www.av在线官网国产| 在线观看一区二区三区| 一区二区三区四区激情视频| 晚上一个人看的免费电影| 久久久久国产网址| 日韩精品有码人妻一区| 久久久久久久久久久免费av| 街头女战士在线观看网站| 黄色怎么调成土黄色| 国产精品久久久久久久久免| 黄色配什么色好看| 精品国产三级普通话版| 亚洲婷婷狠狠爱综合网| 又黄又爽又刺激的免费视频.| 美女脱内裤让男人舔精品视频| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 联通29元200g的流量卡| 日韩一本色道免费dvd| 亚洲最大成人av| 欧美 日韩 精品 国产| h日本视频在线播放| 九色成人免费人妻av| 亚洲精品成人av观看孕妇| h日本视频在线播放| 国产黄a三级三级三级人| 蜜桃亚洲精品一区二区三区| 国产白丝娇喘喷水9色精品| 国产一区亚洲一区在线观看| 99热全是精品| 久久久久久久久久久免费av| 五月玫瑰六月丁香| 国产中年淑女户外野战色| 亚洲国产日韩一区二区| 国产精品av视频在线免费观看| 五月天丁香电影| 国产精品不卡视频一区二区| 免费看a级黄色片| 97热精品久久久久久| 18+在线观看网站| 精品国产三级普通话版| 久久久色成人| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 女人十人毛片免费观看3o分钟| 老女人水多毛片| 日韩av不卡免费在线播放| 啦啦啦中文免费视频观看日本| 国产成人精品一,二区| 国产极品天堂在线| 性插视频无遮挡在线免费观看| 午夜免费男女啪啪视频观看| 免费看av在线观看网站| 男男h啪啪无遮挡| 精品视频人人做人人爽| 婷婷色综合大香蕉| 一本久久精品| av专区在线播放| 日韩中字成人| 下体分泌物呈黄色| 日本一本二区三区精品| 免费观看在线日韩| 国产精品女同一区二区软件| 最近最新中文字幕大全电影3| 天天躁夜夜躁狠狠久久av| 精品少妇久久久久久888优播| 直男gayav资源| 在线观看av片永久免费下载| 久久97久久精品| 亚洲精品一区蜜桃| 街头女战士在线观看网站| 亚洲精品国产av成人精品| 亚洲欧美日韩东京热| 国产男人的电影天堂91| 国产精品久久久久久久久免| 天天一区二区日本电影三级| 国产午夜精品久久久久久一区二区三区| 国语对白做爰xxxⅹ性视频网站| 欧美丝袜亚洲另类| 国产精品国产av在线观看| 日韩中字成人| 简卡轻食公司| 黄色日韩在线| 搡老乐熟女国产| 亚洲人成网站在线播| 99视频精品全部免费 在线| 久久久精品94久久精品| 蜜桃亚洲精品一区二区三区| 人体艺术视频欧美日本| 亚洲最大成人av| 亚洲av在线观看美女高潮| 天天躁夜夜躁狠狠久久av| 欧美激情在线99| 日韩电影二区| 少妇被粗大猛烈的视频| 国内精品宾馆在线| 搞女人的毛片| 亚洲国产欧美在线一区| 免费看av在线观看网站| 欧美日韩一区二区视频在线观看视频在线 | 高清视频免费观看一区二区| 亚洲国产欧美在线一区| 在线亚洲精品国产二区图片欧美 | 成人毛片60女人毛片免费| 亚洲欧美日韩卡通动漫| 波多野结衣巨乳人妻| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 人体艺术视频欧美日本| 日本av手机在线免费观看| 亚洲精品一二三| 精品久久国产蜜桃| 久久久a久久爽久久v久久| 亚洲国产av新网站| 性色av一级| 成人毛片a级毛片在线播放| 国产成人免费无遮挡视频| 久久这里有精品视频免费| av又黄又爽大尺度在线免费看| 国产精品久久久久久久久免| 婷婷色麻豆天堂久久| 一个人观看的视频www高清免费观看| 80岁老熟妇乱子伦牲交| www.色视频.com| 天堂俺去俺来也www色官网| 国产爱豆传媒在线观看| 日本一二三区视频观看| 国产老妇女一区| 国产免费视频播放在线视频| 九草在线视频观看| 观看美女的网站| 久热这里只有精品99| 国产欧美另类精品又又久久亚洲欧美| 国产精品蜜桃在线观看| 国产黄色免费在线视频| 国产片特级美女逼逼视频| 最近中文字幕2019免费版| 视频中文字幕在线观看| 联通29元200g的流量卡| 欧美bdsm另类| a级毛色黄片| 国产又色又爽无遮挡免| 国内揄拍国产精品人妻在线| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 性插视频无遮挡在线免费观看| 午夜免费男女啪啪视频观看| 97超视频在线观看视频| 午夜精品国产一区二区电影 | 亚洲精品亚洲一区二区| 热99国产精品久久久久久7| 看十八女毛片水多多多| 国产精品精品国产色婷婷| 国产日韩欧美在线精品| 亚洲美女视频黄频| 亚洲精品中文字幕在线视频 | 欧美精品国产亚洲| 久久久久久久久久久丰满| 国产一区二区三区av在线| 最近的中文字幕免费完整| 国产午夜福利久久久久久| 精品久久久久久久末码| 日韩,欧美,国产一区二区三区| 国产一区二区亚洲精品在线观看| 观看免费一级毛片| 美女主播在线视频| 青春草国产在线视频| 亚洲精品视频女| 99热全是精品| 日韩人妻高清精品专区| 老司机影院成人| 国内精品美女久久久久久| 一本久久精品| 人体艺术视频欧美日本| videossex国产| 网址你懂的国产日韩在线| 欧美激情在线99| 日韩,欧美,国产一区二区三区| 国产精品精品国产色婷婷| a级一级毛片免费在线观看| 一级a做视频免费观看| 国产免费又黄又爽又色| 中国三级夫妇交换| 在线观看人妻少妇| 99热网站在线观看| 久久精品国产鲁丝片午夜精品| 26uuu在线亚洲综合色| 91精品国产九色| av播播在线观看一区| 青春草国产在线视频| 神马国产精品三级电影在线观看| 国产成人aa在线观看| 久久精品人妻少妇| 久久久色成人| 人妻夜夜爽99麻豆av| 在现免费观看毛片| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 亚洲第一区二区三区不卡| 亚洲最大成人手机在线| 亚洲电影在线观看av| 男插女下体视频免费在线播放| 精品国产一区二区三区久久久樱花 | 边亲边吃奶的免费视频| 国产乱来视频区| 又大又黄又爽视频免费| 亚洲美女搞黄在线观看| 精品国产三级普通话版| 狠狠精品人妻久久久久久综合| 欧美日韩视频精品一区| 欧美极品一区二区三区四区| 一级毛片 在线播放| av专区在线播放| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 在线观看美女被高潮喷水网站| 人人妻人人爽人人添夜夜欢视频 | 99热全是精品| 久久99蜜桃精品久久| 婷婷色av中文字幕| 黄片wwwwww| 狠狠精品人妻久久久久久综合| 亚洲美女视频黄频| 日日啪夜夜爽| 一本色道久久久久久精品综合| 国产一区有黄有色的免费视频| 777米奇影视久久| 国产欧美日韩精品一区二区| 国国产精品蜜臀av免费| 亚洲欧美成人精品一区二区| 欧美一区二区亚洲| 天天一区二区日本电影三级| 一二三四中文在线观看免费高清| 五月开心婷婷网| 性插视频无遮挡在线免费观看| 美女脱内裤让男人舔精品视频| 禁无遮挡网站| 交换朋友夫妻互换小说| 国产乱来视频区| 中文欧美无线码| 国产精品久久久久久久久免| 成人毛片a级毛片在线播放| 大片免费播放器 马上看|