• <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á)的差異
    男的添女的下面高潮视频| 99视频精品全部免费 在线| 丝袜喷水一区| 在线观看人妻少妇| 国产免费又黄又爽又色| 免费黄频网站在线观看国产| 国产 一区 欧美 日韩| 在线观看人妻少妇| 亚洲图色成人| 欧美高清性xxxxhd video| 成人鲁丝片一二三区免费| 一本一本综合久久| 美女cb高潮喷水在线观看| 最近手机中文字幕大全| 午夜爱爱视频在线播放| 免费看日本二区| 久久精品国产a三级三级三级| 男人添女人高潮全过程视频| 尾随美女入室| 在线天堂最新版资源| 久久久久久九九精品二区国产| 国产一区二区三区av在线| 成人国产av品久久久| 一个人观看的视频www高清免费观看| 欧美成人a在线观看| 国产精品嫩草影院av在线观看| 欧美3d第一页| 成人二区视频| 99热这里只有是精品在线观看| 国产精品久久久久久久久免| 国产乱人偷精品视频| 天堂网av新在线| 亚洲美女搞黄在线观看| 国产 一区 欧美 日韩| 午夜免费观看性视频| 欧美日本视频| 看免费成人av毛片| 插逼视频在线观看| 两个人的视频大全免费| 国产乱人视频| 日韩伦理黄色片| 亚洲成人中文字幕在线播放| 久久久成人免费电影| 婷婷色综合大香蕉| 久久久久久九九精品二区国产| 国产精品一二三区在线看| 97精品久久久久久久久久精品| 亚洲精品第二区| 国产淫片久久久久久久久| 午夜爱爱视频在线播放| 天堂俺去俺来也www色官网| 亚洲av不卡在线观看| 尤物成人国产欧美一区二区三区| 老女人水多毛片| 国产成人精品福利久久| 少妇丰满av| 国产亚洲最大av| 亚洲国产精品成人久久小说| 一级爰片在线观看| 1000部很黄的大片| 一本色道久久久久久精品综合| 天天躁日日操中文字幕| 18禁动态无遮挡网站| 亚洲,欧美,日韩| 久久精品人妻少妇| 免费大片18禁| a级毛色黄片| 80岁老熟妇乱子伦牲交| 免费看a级黄色片| 亚洲精品国产色婷婷电影| 亚洲欧美中文字幕日韩二区| 国产在线男女| 亚洲欧美精品专区久久| 免费观看性生交大片5| 国产亚洲一区二区精品| 国产av国产精品国产| 免费黄色在线免费观看| 高清视频免费观看一区二区| 亚洲成人中文字幕在线播放| 欧美xxxx黑人xx丫x性爽| 内地一区二区视频在线| 国产有黄有色有爽视频| 少妇的逼好多水| 日韩亚洲欧美综合| 小蜜桃在线观看免费完整版高清| 国产真实伦视频高清在线观看| 久久久精品欧美日韩精品| 免费黄色在线免费观看| 伊人久久精品亚洲午夜| 久久精品国产自在天天线| 嫩草影院新地址| 久久精品夜色国产| av黄色大香蕉| 免费av不卡在线播放| 亚洲国产av新网站| 免费观看的影片在线观看| 熟妇人妻不卡中文字幕| 伦理电影大哥的女人| 男女下面进入的视频免费午夜| 久久6这里有精品| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 久久久成人免费电影| 午夜老司机福利剧场| 中文在线观看免费www的网站| 一级毛片黄色毛片免费观看视频| 亚州av有码| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 久久久午夜欧美精品| 国产视频首页在线观看| 18禁动态无遮挡网站| 少妇高潮的动态图| 国产成年人精品一区二区| 一本一本综合久久| 国产伦理片在线播放av一区| 男女下面进入的视频免费午夜| 久久6这里有精品| 九九爱精品视频在线观看| 欧美最新免费一区二区三区| 一边亲一边摸免费视频| 国产精品一区www在线观看| 亚洲丝袜综合中文字幕| 少妇人妻一区二区三区视频| 一本久久精品| 国产在线男女| 久久精品久久久久久噜噜老黄| 国产成人a∨麻豆精品| 精品亚洲乱码少妇综合久久| 国产高潮美女av| 亚洲国产欧美人成| 91精品国产九色| 女人久久www免费人成看片| www.色视频.com| 亚洲国产最新在线播放| 精品视频人人做人人爽| 亚洲av中文av极速乱| 国产精品女同一区二区软件| 亚洲图色成人| 久久精品人妻少妇| 大片免费播放器 马上看| 中文字幕人妻熟人妻熟丝袜美| 成年女人看的毛片在线观看| 久久影院123| 晚上一个人看的免费电影| 女人久久www免费人成看片| 免费av毛片视频| 日本欧美国产在线视频| 亚洲性久久影院| av国产免费在线观看| 嫩草影院入口| 久久久精品免费免费高清| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 久久久久久久国产电影| 国产精品一区二区性色av| 亚洲精品国产成人久久av| 热re99久久精品国产66热6| 九草在线视频观看| 各种免费的搞黄视频| 简卡轻食公司| 国产成人精品福利久久| 麻豆国产97在线/欧美| 麻豆乱淫一区二区| 成人亚洲欧美一区二区av| 亚洲自偷自拍三级| 极品教师在线视频| 国产成人精品久久久久久| 亚洲欧洲日产国产| 中文字幕免费在线视频6| 国产有黄有色有爽视频| 欧美一区二区亚洲| 国产成人精品一,二区| 男女边摸边吃奶| 日本av手机在线免费观看| 久久久久久久久久久丰满| 大话2 男鬼变身卡| 美女xxoo啪啪120秒动态图| 中文字幕免费在线视频6| 在线天堂最新版资源| 国国产精品蜜臀av免费| 蜜桃久久精品国产亚洲av| 国产成人精品久久久久久| 激情 狠狠 欧美| 成人亚洲精品av一区二区| 肉色欧美久久久久久久蜜桃 | 麻豆精品久久久久久蜜桃| 日韩国内少妇激情av| 精品久久久久久电影网| 亚洲人与动物交配视频| 国产男女超爽视频在线观看| 色哟哟·www| 成人二区视频| 91久久精品电影网| 精品久久久久久久末码| 中国美白少妇内射xxxbb| 少妇人妻 视频| 久久久久网色| 亚洲色图av天堂| 欧美精品国产亚洲| 久久久久久久大尺度免费视频| 一级毛片黄色毛片免费观看视频| 偷拍熟女少妇极品色| 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 一区二区三区免费毛片| 国产伦理片在线播放av一区| 黑人高潮一二区| 天堂俺去俺来也www色官网| av福利片在线观看| 一级二级三级毛片免费看| 99久久人妻综合| 久久久久久九九精品二区国产| 黄色怎么调成土黄色| videos熟女内射| 亚洲精品乱码久久久久久按摩| 一区二区三区四区激情视频| 成人综合一区亚洲| 欧美日韩视频精品一区| 亚洲国产色片| 中文字幕制服av| 你懂的网址亚洲精品在线观看| 久久精品综合一区二区三区| 老师上课跳d突然被开到最大视频| 天堂中文最新版在线下载 | 亚洲天堂国产精品一区在线| 国产精品无大码| 国产亚洲精品久久久com| 欧美日本视频| 国产av不卡久久| 91精品国产九色| av在线蜜桃| 日本爱情动作片www.在线观看| 久久久久久久亚洲中文字幕| 日韩伦理黄色片| 国产老妇伦熟女老妇高清| 免费看光身美女| 精品午夜福利在线看| 热re99久久精品国产66热6| 你懂的网址亚洲精品在线观看| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 汤姆久久久久久久影院中文字幕| 国产亚洲精品久久久com| 久久久久久久久久成人| 国产亚洲午夜精品一区二区久久 | 97超视频在线观看视频| 欧美日韩亚洲高清精品| 日韩av在线免费看完整版不卡| 天天躁日日操中文字幕| 亚洲在久久综合| 各种免费的搞黄视频| 亚洲精品中文字幕在线视频 | 精品人妻一区二区三区麻豆| 97人妻精品一区二区三区麻豆| 久久久色成人| 精品国产一区二区三区久久久樱花 | 在线精品无人区一区二区三 | 大香蕉97超碰在线| 亚洲经典国产精华液单| 黄色怎么调成土黄色| 精品久久久噜噜| eeuss影院久久| 精品人妻熟女av久视频| 高清午夜精品一区二区三区| 69av精品久久久久久| 波多野结衣巨乳人妻| 日日撸夜夜添| 蜜桃久久精品国产亚洲av| 男人舔奶头视频| 日韩免费高清中文字幕av| 欧美人与善性xxx| 一级片'在线观看视频| 中文字幕制服av| 日韩人妻高清精品专区| 嫩草影院入口| 中文在线观看免费www的网站| 亚洲欧美一区二区三区黑人 | 亚洲成人精品中文字幕电影| 亚洲综合色惰| 色婷婷久久久亚洲欧美| 国产 一区 欧美 日韩| 大陆偷拍与自拍| 视频区图区小说| 卡戴珊不雅视频在线播放| 国产爽快片一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产伦精品一区二区三区视频9| 一区二区三区乱码不卡18| 国产精品久久久久久久久免| www.av在线官网国产| 亚洲欧美成人综合另类久久久| 亚洲高清免费不卡视频| 乱码一卡2卡4卡精品| 永久网站在线| www.色视频.com| 一本一本综合久久| 精品久久久久久久人妻蜜臀av| 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 看免费成人av毛片| 校园人妻丝袜中文字幕| 另类亚洲欧美激情| xxx大片免费视频| 日本三级黄在线观看| 2021天堂中文幕一二区在线观| 亚洲国产色片| 在线看a的网站| 男男h啪啪无遮挡| 伦精品一区二区三区| 有码 亚洲区| 在线观看人妻少妇| 美女国产视频在线观看| 国产91av在线免费观看| 精品人妻视频免费看| 最近的中文字幕免费完整| 在线天堂最新版资源| 亚洲aⅴ乱码一区二区在线播放| 日韩成人伦理影院| 久久久久久国产a免费观看| 国产综合懂色| www.色视频.com| 如何舔出高潮| av免费观看日本| 久久久久性生活片| 亚洲精品久久久久久婷婷小说| 热99国产精品久久久久久7| 成年版毛片免费区| 国产成人免费无遮挡视频| 尾随美女入室| 亚洲最大成人av| 美女被艹到高潮喷水动态| 免费电影在线观看免费观看| 又爽又黄a免费视频| 黄片wwwwww| 建设人人有责人人尽责人人享有的 | 久久久欧美国产精品| 国产亚洲91精品色在线| 国产色婷婷99| 777米奇影视久久| 伊人久久国产一区二区| 亚洲国产av新网站| av线在线观看网站| 一区二区av电影网| 久久精品国产a三级三级三级| 免费观看a级毛片全部| 人妻系列 视频| 精品人妻偷拍中文字幕| 亚洲国产色片| 午夜福利高清视频| 亚洲欧美一区二区三区国产| 久久精品人妻少妇| 99热6这里只有精品| 午夜福利高清视频| 欧美xxxx性猛交bbbb| av在线蜜桃| 毛片女人毛片| 天堂中文最新版在线下载 | 亚洲在久久综合| 久久影院123| 午夜免费男女啪啪视频观看| 亚洲精品国产色婷婷电影| 可以在线观看毛片的网站| 精品午夜福利在线看| 真实男女啪啪啪动态图| 五月玫瑰六月丁香| 男女无遮挡免费网站观看| 丰满人妻一区二区三区视频av| 日本三级黄在线观看| 免费播放大片免费观看视频在线观看| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 欧美高清性xxxxhd video| 26uuu在线亚洲综合色| 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 肉色欧美久久久久久久蜜桃 | 中文天堂在线官网| 狂野欧美激情性bbbbbb| 精品人妻偷拍中文字幕| 久久人人爽人人片av| 亚洲av中文av极速乱| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 久久精品夜色国产| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 久久99精品国语久久久| 一级av片app| 亚洲天堂av无毛| 国产精品久久久久久精品古装| 国产精品无大码| 女人被狂操c到高潮| 国产真实伦视频高清在线观看| 日本黄色片子视频| 国产精品福利在线免费观看| 国产v大片淫在线免费观看| 亚洲熟女精品中文字幕| 中文字幕制服av| 欧美bdsm另类| 国产久久久一区二区三区| 国产精品秋霞免费鲁丝片| av线在线观看网站| 国产精品国产三级国产av玫瑰| 极品少妇高潮喷水抽搐| 黑人高潮一二区| 夜夜爽夜夜爽视频| 精品久久久久久久久av| 国产又色又爽无遮挡免| 视频中文字幕在线观看| 大香蕉久久网| 日韩欧美精品免费久久| 在线看a的网站| 麻豆成人午夜福利视频| 婷婷色麻豆天堂久久| 韩国av在线不卡| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 国产毛片在线视频| 精品少妇久久久久久888优播| 精品少妇黑人巨大在线播放| 最近中文字幕高清免费大全6| 在线播放无遮挡| 免费电影在线观看免费观看| 亚洲精品一区蜜桃| 色吧在线观看| 狠狠精品人妻久久久久久综合| 午夜福利高清视频| 秋霞在线观看毛片| 日韩欧美精品v在线| 欧美日韩亚洲高清精品| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 亚州av有码| 亚洲三级黄色毛片| 激情五月婷婷亚洲| 亚洲无线观看免费| 亚洲色图综合在线观看| 大陆偷拍与自拍| 99久久精品国产国产毛片| 日本黄大片高清| 男人爽女人下面视频在线观看| 久久久久久九九精品二区国产| 在线免费观看不下载黄p国产| 亚洲av中文av极速乱| 日本黄色片子视频| 免费播放大片免费观看视频在线观看| 国产日韩欧美亚洲二区| 高清av免费在线| 亚洲av中文字字幕乱码综合| 国产综合精华液| 看免费成人av毛片| 一本久久精品| 日韩一本色道免费dvd| a级毛色黄片| 精华霜和精华液先用哪个| 亚洲不卡免费看| 色吧在线观看| 人妻 亚洲 视频| 一级二级三级毛片免费看| 久久鲁丝午夜福利片| 久久午夜福利片| 免费观看在线日韩| 亚洲人与动物交配视频| 亚洲天堂av无毛| 亚洲不卡免费看| 亚州av有码| 亚洲av中文字字幕乱码综合| 日本av手机在线免费观看| 亚洲经典国产精华液单| 亚洲精品日本国产第一区| 亚洲精品视频女| 又爽又黄无遮挡网站| 亚洲av成人精品一二三区| av在线老鸭窝| 狂野欧美激情性xxxx在线观看| 99re6热这里在线精品视频| 天天躁日日操中文字幕| 人妻夜夜爽99麻豆av| 少妇熟女欧美另类| 精品国产乱码久久久久久小说| av线在线观看网站| 亚洲精品久久午夜乱码| kizo精华| 91久久精品国产一区二区成人| 国产高清不卡午夜福利| 久久久久精品性色| 国产精品无大码| 亚洲精品一区蜜桃| 精品久久久久久久久av| 十八禁网站网址无遮挡 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利视频精品| 免费看日本二区| 国产极品天堂在线| 国产精品99久久99久久久不卡 | 美女视频免费永久观看网站| 午夜福利网站1000一区二区三区| 久久国内精品自在自线图片| 亚洲av二区三区四区| 色视频在线一区二区三区| 97超碰精品成人国产| a级毛片免费高清观看在线播放| 天堂网av新在线| 亚洲精品色激情综合| 国产大屁股一区二区在线视频| 久久久久精品久久久久真实原创| 成人鲁丝片一二三区免费| 亚洲av欧美aⅴ国产| 卡戴珊不雅视频在线播放| 麻豆成人av视频| 男女边吃奶边做爰视频| 中文字幕久久专区| 五月玫瑰六月丁香| 国产色爽女视频免费观看| 久久久成人免费电影| 亚洲久久久久久中文字幕| 26uuu在线亚洲综合色| 晚上一个人看的免费电影| 日韩欧美精品v在线| 久久精品久久久久久久性| 神马国产精品三级电影在线观看| 亚洲人与动物交配视频| 国产爱豆传媒在线观看| 久久久久网色| 亚洲av日韩在线播放| 国产精品嫩草影院av在线观看| 日本-黄色视频高清免费观看| 免费观看的影片在线观看| 午夜福利在线观看免费完整高清在| 日本欧美国产在线视频| av网站免费在线观看视频| 亚洲人与动物交配视频| 亚洲av欧美aⅴ国产| 99热6这里只有精品| 亚洲av成人精品一二三区| 国产大屁股一区二区在线视频| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 秋霞在线观看毛片| 制服丝袜香蕉在线| 在现免费观看毛片| 成人亚洲欧美一区二区av| 国模一区二区三区四区视频| 夫妻性生交免费视频一级片| 日本一二三区视频观看| 国产免费福利视频在线观看| 男人和女人高潮做爰伦理| 白带黄色成豆腐渣| 少妇猛男粗大的猛烈进出视频 | 亚洲av成人精品一二三区| 欧美日韩亚洲高清精品| 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 3wmmmm亚洲av在线观看| 久久久久久久大尺度免费视频| 又大又黄又爽视频免费| 久久久久久久久久久丰满| 美女国产视频在线观看| 国产色爽女视频免费观看| 五月玫瑰六月丁香| 久久精品熟女亚洲av麻豆精品| 夫妻午夜视频| 欧美高清成人免费视频www| 美女xxoo啪啪120秒动态图| 十八禁网站网址无遮挡 | 国产精品一区二区性色av| 熟女人妻精品中文字幕| 亚洲av国产av综合av卡| 国产精品伦人一区二区| 黄片wwwwww| 亚洲国产精品国产精品| 我要看日韩黄色一级片| 99精国产麻豆久久婷婷| a级毛色黄片| 各种免费的搞黄视频| 欧美高清成人免费视频www| 精品久久久精品久久久| 国产精品成人在线| 嘟嘟电影网在线观看| 国产国拍精品亚洲av在线观看| 99re6热这里在线精品视频| 日韩亚洲欧美综合| 99热国产这里只有精品6| 亚洲三级黄色毛片| 久久久精品欧美日韩精品| 欧美+日韩+精品| 国产精品av视频在线免费观看| 成人漫画全彩无遮挡| 久久这里有精品视频免费| 美女xxoo啪啪120秒动态图| 国产伦在线观看视频一区| 亚洲欧美日韩卡通动漫| 久久97久久精品| 国产男人的电影天堂91| 伊人久久国产一区二区| 日日撸夜夜添| 国产精品国产av在线观看| 国产免费一区二区三区四区乱码| 亚洲真实伦在线观看| 99热这里只有是精品50| 啦啦啦中文免费视频观看日本| 国产高潮美女av| 久久精品国产a三级三级三级| 久久影院123| 亚洲欧美一区二区三区黑人 | 日韩精品有码人妻一区| 18+在线观看网站| 亚洲欧美清纯卡通| 亚洲熟女精品中文字幕| av专区在线播放| 一区二区三区免费毛片| 国产午夜精品久久久久久一区二区三区| 日韩视频在线欧美| 国产精品久久久久久av不卡|