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

    高寒草甸退化草地植被與土壤因子關(guān)系冗余分析

    2018-05-10 07:23:51劉育紅魏衛(wèi)東楊元武
    西北農(nóng)業(yè)學(xué)報 2018年4期
    關(guān)鍵詞:草甸土壤溫度全氮

    劉育紅 ,魏衛(wèi)東,楊元武,張 英

    (青海大學(xué) 農(nóng)牧學(xué)院,西寧 810016)

    高寒草甸在惡劣的環(huán)境因素和不合理的人類利用共同驅(qū)動下出現(xiàn)退化問題[1],使得其生態(tài)環(huán)境惡化,生態(tài)、生產(chǎn)功能受到較大影響,區(qū)域生態(tài)保護(hù)和社會發(fā)展功能被削弱[2]。由此,研究人員針對高寒草甸退化、退化草地土壤特征等問題開展相應(yīng)研究,并取得一定成果。一段時間以來,研究主要集中在退化草地植物群落及特征與土壤性質(zhì)或土壤環(huán)境關(guān)系方面[3-5];不同退化階段土壤化學(xué)性質(zhì)變化方面[6-8],退化草甸土壤有機(jī)碳、有機(jī)碳組分、土壤碳儲量、土壤微生物量碳方面[9-13];退化草地土壤粒徑等物理性質(zhì)方面[14-17];退化草地土壤酶活性、土壤呼吸方面[18-19]。但是,這些研究側(cè)重探討的是退化高寒草甸植物群落以及土壤中碳、氮、磷、土壤顆粒、土壤含水量等因子單獨或少量幾個因子結(jié)合后的變化特征,而將退化草地土壤多個理化因子相結(jié)合,深入探討退化草地植被與土壤理化因子間關(guān)系的研究較少。從高寒草甸不同退化程度草地植被與草地土壤理化因子關(guān)系角度看,存在的主要科學(xué)問題有:退化草地植被與土壤理化因子間呈現(xiàn)怎樣的關(guān)系,退化草地植被與土壤理化因子間關(guān)系的密切程度,哪些土壤因子能夠更敏感地反映退化草地植被的變化。基于這些科學(xué)問題,選擇具有典型高寒草甸特征的研究樣地,進(jìn)行草地退化現(xiàn)狀調(diào)查、植物群落特征觀測、土壤樣品采集及測定,同時運用數(shù)量生態(tài)學(xué)中的冗余分析(redundancy analysis,RDA)方法,研究高寒草甸退化草地植被與草地土壤因子的內(nèi)在聯(lián)系,以期為高寒草甸退化草地的相關(guān)研究和生態(tài)治理提供理論依據(jù),也為高寒地區(qū)草地群落、物種與環(huán)境間關(guān)系的研究提供實例。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)在青海省果洛藏族自治州(以下簡稱果洛州)瑪沁縣(MQ)、甘德縣(GD)、達(dá)日縣(DR)。果洛州(N32°21′~35°45′,E96°56′~101°45′)總面積7.6萬 km2,可利用草地面積585.6萬 hm2,其中高寒草甸占 56.3%;全州平均海拔4 200 m,大氣含氧量僅為海平面的60%;果洛州屬高原大陸型氣候,氣溫低、溫差大、輻射強(qiáng)、年均溫-4 ℃,年均降水量513.2 mm,年均蒸發(fā)量1 462.4 mm,年均日照時數(shù)2 260.3 h,全年無絕對無霜期;果洛州境內(nèi)河流縱橫、湖泊眾多、黃河自西至東橫穿全境;境內(nèi)高寒草甸草地植物群落中莎草科高山嵩草(Kobresiapygmaea)、矮嵩草(Kobresiahumilis)等為優(yōu)勢種,伴生種以禾本科(Gramineae)、菊科(Compositae)、龍膽科(Gentianaceae)、薔薇科(Rosaceae)、豆科(Leguminosae) 等植物為主;高寒草甸草地土壤類型以高山草甸土分布最多。

    1.2 樣地設(shè)置

    在具有典型高寒草甸特征的果洛州瑪沁縣、甘德縣、達(dá)日縣天然草地設(shè)置研究樣地(表1),利用空間分布代替時間演替的方法研究高寒草甸退化草地植被與土壤因子的關(guān)系。依據(jù)任繼周[20]、趙新全[21]關(guān)于草地退化程度劃分的方法,結(jié)合草地類型、土壤侵蝕現(xiàn)狀、嚙齒動物危害等情況綜合評價樣地的退化程度。研究樣地處于陽坡,坡度11°~21°;為高山嵩草草地型,主要植物為高山嵩草、矮嵩草、線葉嵩草(Kobresiacapillifolia)、垂穗披堿草(Elymusnutans)等;土壤為高山草甸土;分為4種退化程度,即:未退化(un-degradation,UD)、輕度退化(light degradation,LD)、中度退化(moderate degradation,MD)和重度退化(heavy degradation,HD)。每個研究樣地大小為50 m×50 m,在每個樣地內(nèi)進(jìn)行植物群落數(shù)量特征觀測時隨機(jī)設(shè)置觀測樣方,重復(fù)6次。

    1.3 測定項目與方法

    2016年8月,在各研究樣地隨機(jī)設(shè)置6個1 m2的觀測樣方,對每個樣方進(jìn)行植物群落特征調(diào)查。包括植被總蓋度(多人目測平均法)、植物物種數(shù)及不同物種蓋度、高度(自地表至植株頂端自然高度,每樣方測定20株,不足20株的物種按實際株數(shù)測定)、地上生物量(分物種齊地面刈割稱鮮質(zhì)量)等;重要值作為退化草地植物優(yōu)勢度指標(biāo)。

    重要值=(相對蓋度+相對高度+相對生物量)/3。

    完成植物群落特征調(diào)查后,在每個樣方內(nèi)土鉆采集0~30 cm土層的土樣,同一樣地不同樣方土樣混合為1個土壤樣品,帶回實驗室自然風(fēng)干后用于測定土壤性質(zhì)。在采集土樣的同時,環(huán)刀法測定土壤容質(zhì)量,測溫計測定土壤溫度,水分測定儀(TDR 300)測定土壤含水量。土壤有機(jī)碳采用重鉻酸鉀容量法[22]測定;土壤全氮采用半微量凱氏定氮法、土壤全磷采用鉬銻抗比色法[23]測定;土壤有效氮采用比色法、土壤有效磷采用碳酸氫鈉浸提鉬銻抗比色法、土壤全鉀和土壤有效鉀采用火焰光度法[22]測定;采用比重計法[24]測定粘粒比例(粒徑<0.002 mm)和砂粒比例(0.05 mm<粒徑<2 mm);土壤pH采用水土體積比2.5∶1 的電位法[22]測定。

    1.4 數(shù)據(jù)分析方法

    應(yīng)用Canoco 4.5軟件基于線性模型對退化草地植被-土壤因子進(jìn)行冗余分析[25]。其中物種數(shù)據(jù)為26種植物構(gòu)建的26×20維重要值矩陣;退化草地土壤主要理化因子構(gòu)建13×20維土壤環(huán)境因子矩陣,包括有機(jī)碳(soil organic carbon,SOC)、全氮(total nitrogen,TN)、有效氮(available nitrogen,AN)、全磷(total phosphorus,TP)、有效磷(available phosphorus,AP)、全鉀(total potassium,TK)、有效鉀(available potassium,AK)、容質(zhì)量(bulk density,BD)、土壤含水量(soil water content,SWC)、土壤溫度(soil temperature,ST)、粘粒比例(clay,以下簡稱粘粒)、砂粒比例(sand,以下簡稱沙粒)、pH共13個指標(biāo)。進(jìn)行冗余分析時,對退化草地植被數(shù)據(jù)進(jìn)行中心化和標(biāo)準(zhǔn)化,排序軸特征值采用Monte Carlo permutation test檢驗顯著性;通過前向選擇法(forward selection)Monte Carlo permutation test隨機(jī)置換999次檢驗13個土壤因子邊際作用(marginal effects)[26],以此反映土壤因子對退化草地植被影響的顯著性,并按照其特征值(eigenvalues)進(jìn)行重要性排序;Cano Draw繪制植被-土壤因子RDA二維排序圖;SPSS 20.0進(jìn)行Spearman相關(guān)分析。

    2 結(jié)果與分析

    2.1 高寒草甸退化草地植物群落組成

    植被是高寒草甸生態(tài)系統(tǒng)的重要組成部分。隨著高寒草甸退化程度的變化,植物種類、物種重要值均發(fā)生相應(yīng)變化。在退化程度不同的研究樣地,物種總數(shù)呈LD>MD>UD>HD的變化趨勢;UD樣地,莎草科高山嵩草、矮嵩草、青藏苔草和禾本科垂穗披堿草具有較高優(yōu)勢度,雜類草及豆科植物優(yōu)勢度較低且種類相對較少;LD、MD樣地雜類草和豆科植物開始大量出現(xiàn),植物種類最多;HD樣地高山嵩草、青藏苔草等退出,黃帚橐吾、全緣葉綠絨蒿等雜類草成為優(yōu)勢種。莎草科植物重要值總體隨草地退化程度的加劇而減小,禾本科垂穗披堿草重要值也具同樣規(guī)律,說明隨著草地退化,優(yōu)良牧草在植物群落中的優(yōu)勢地位逐漸消失;雜類草中美麗風(fēng)毛菊、乳白香青、釘柱委陵菜等重要值在中度退化時較大,黃帚橐吾、全緣葉綠絨蒿、矮生忍冬、細(xì)葉亞菊等則在重度退化時較大,反映出隨著草地退化,雜類草在植物群落中逐漸占優(yōu)勢,草地生產(chǎn)功能隨之降低乃至喪失(表2)。

    2.2 高寒草甸植被與草地退化的關(guān)系

    利用研究樣地植被-土壤因子數(shù)據(jù)矩陣進(jìn)行冗余分析發(fā)現(xiàn),RDA排序圖中不同退化程度樣地分布在4個區(qū)域。樣地號為1、5、6、12、13的退化程度為UD,2、7、8、14、15為LD, 3、4、9、16、17為MD;樣地號為10、11、18、19、20為HD。說明RDA排序圖第1排序軸主要反映退化程度,自左向右草地退化程度加劇。另外,沿第1排序軸水平方向UD、LD樣地位置臨近,并與MD、HD樣地位置相對較遠(yuǎn),說明不同退化程度間存在差異(圖1)。

    從植物在退化草地的分布格局看,26種植物呈明顯的區(qū)域分布。物種S1、S2、S3、S4、S5、S12、S14集中分布在UD樣地,與UD樣地關(guān)聯(lián)度較高;物種S6、S7、S8、S9、S11、S13、S22、S23、S25、S26集中分布在MD樣地,與MD樣地關(guān)聯(lián)度較高;物種S10、S15、S16、S17、S18、S19、S20、S21、S24則集中分布在HD樣地,與HD樣地關(guān)聯(lián)度較高(圖1)。沿第1排序軸,物種由分布在UD樣地的高山嵩草、矮嵩草等向分布在HD樣地的黃帚橐吾等轉(zhuǎn)變,說明植物物種分布不同,植物物種分布與草地退化程度相關(guān)。

    表2 高寒草甸不同退化程度草地植物物種組成及重要值Table 2 Species composition and important value of alpine meadow plants on different degraded grassland

    注:“/” 表示該物種未出現(xiàn)。

    Note:“/” Species did not appear.

    綜上所述,當(dāng)高寒草甸草地退化加劇時,植被中莎草科嵩草屬植物在群落中的優(yōu)勢地位逐漸被菊科、豆科的毒雜草植物取代,黃帚橐吾、細(xì)葉亞菊、黃花棘豆等毒雜草是表征草地重度退化的重要物種。

    2.3 高寒草甸退化草地植被與土壤因子的關(guān)系

    對研究樣地植被-土壤因子數(shù)據(jù)進(jìn)行冗余分析,結(jié)果表明:第1、第2排序軸特征值分別為0.504、0.240,2個排序軸共解釋74.5%的草地植被變化和83.2%的植被-土壤因子關(guān)系。說明第1、第2排序軸較好地反映退化草地植被與土壤因子間的關(guān)系。退化草地植被與土壤因子相關(guān)系數(shù)在第1、第2排序軸分別為0.995、0.982,說明退化草地植被與土壤因子關(guān)系密切;Monte Carlo permutation test也表明第1排序軸及所有排序軸反映的退化草地植被與土壤因子間呈極顯著相關(guān)關(guān)系(P<0.01)(表3)。

    冗余分析結(jié)果顯示,土壤溫度、容質(zhì)量與第1排序軸正相關(guān),相關(guān)系數(shù)分別為0.947、0.915;土壤有效氮、土壤含水量、全氮、全磷與第1排序軸負(fù)相關(guān),相關(guān)系數(shù)分別為-0.899、-0.881、-0.880、-0.818,說明第1排序軸主要反映土壤溫度、容質(zhì)量、土壤含水量以及全氮、全磷等因子的綜合變化,沿第1排序軸自左到右,土壤溫度升高、容質(zhì)量增大、有效氮、全氮、全磷含量降低、土壤含水量減少。土壤粘粒、有效鉀與第2排序軸正相關(guān),相關(guān)系數(shù)分別為0.851、0.804;土壤沙粒與第2排序軸負(fù)相關(guān),相關(guān)系數(shù)為-0.896,說明第2排序軸主要反映土壤機(jī)械組成的綜合變化,沿第2排序軸自下向上,粘粒比例、有效鉀含量增加,沙粒比例減少(表4、圖2)。

    對13個土壤因子邊際作用檢驗后發(fā)現(xiàn),不同土壤因子對草地植被分布的影響不同。依據(jù)土壤因子特征值排序為土壤溫度>容質(zhì)量>有效氮>全氮>土壤含水量>全磷>有機(jī)碳>有效鉀>粘粒>沙粒>有效磷>pH>全鉀,說明土壤溫度、有效氮等理化指標(biāo)對不同退化程度草地植被的分布更重要(表4)。

    圖中數(shù)字代表樣地號,S1-S26分別代表不同植物物種 Numbers represent the plots, S1-S26 represents the different plant species

    圖1 研究樣地及植被RDA二維排序圖Fig.1 RDA two-dimensional ordination diagram of plots and vegetation

    圖2反映退化草地植被-土壤因子冗余分析二維排序結(jié)果。結(jié)合圖2、圖1可以看出,UD樣地中有效氮、全氮、土壤含水量、有機(jī)碳、全磷較其他土壤因子對物種S1、S2、S3、S4、S5、S12、S14的影響更大;LD、MD樣地中有效鉀、粘粒、有效磷、pH較其他土壤因子對物種S6、S7、S8、S9、S11、S13、S22、S23、S25、S26的影響更大;HD樣地中,土壤溫度、容質(zhì)量、沙粒、全鉀較其他土壤因子對物種S10、S15、S16、S17、S18、S19、S20、S21、S24的影響更大。說明隨著草地退化程度的加劇,在土壤有效氮、全氮等因子水平減小,土壤溫度、容質(zhì)量等因子水平增大的背景下,草地植被由莎草科、禾本科逐步過渡到雜類草,最終演替為毒雜草群落。同時,當(dāng)土壤有效氮、全氮、土壤含水量等水平較高,土壤溫度、容質(zhì)量等水平較低時,有利于高山嵩草、垂穗披堿草等高寒草甸優(yōu)勢種植物在群落中的分布;反之,則有利于黃帚橐吾、矮生忍冬等毒雜草植物在群落中的分布。說明草地植被演替與土壤因子變化相關(guān)。

    不同植物與不同土壤因子間的相關(guān)性分析表明:莎草科高山嵩草、青藏苔草,禾本科垂穗披堿草與有機(jī)碳、全氮、有效氮、全磷、土壤含水量極顯著正相關(guān)(P<0.01),與容質(zhì)量、土壤溫度極顯著負(fù)相關(guān)(P<0.01);雜類草中,細(xì)葉亞菊、唐古特虎耳草等與有機(jī)碳、全氮、有效氮、全磷、土壤含水量極顯著負(fù)相關(guān)(P<0.01),與容質(zhì)量、土壤溫度極顯著正相關(guān)(P<0.01);豆科植物中黃花棘豆同樣與有機(jī)碳、全氮、有效氮、全磷、土壤含水量極顯著負(fù)相關(guān)(P<0.01),而與容質(zhì)量、土壤溫度極顯著正相關(guān)(P<0.01)(表5)。進(jìn)一步說明不同土壤因子對不同植物分布的影響不同。其次,與研究樣地植被存在極顯著、顯著相關(guān)較多的土壤因子是全氮、有效氮、全磷、有機(jī)碳、容質(zhì)量、土壤溫度、土壤含水量,這一結(jié)果與冗余分析中對土壤因子邊際作用檢驗并按照特征值排序的結(jié)果一致(表5)。

    圖2 退化草地植被與土壤因子RDA二維排序圖Fig.2 Degraded grassland vegetation and RDA two-dimensional ordination diagram of soil factors

    3 結(jié)論與討論

    青藏高原高寒草甸由于受獨特的氣候、地理等因素綜合作用,生態(tài)系統(tǒng)結(jié)構(gòu)相對簡單,易受環(huán)境和人為因素的擾動而發(fā)生草地退化,導(dǎo)致草地植被與土壤因子間相互響應(yīng)與反饋。同時,由于草地植被物種間存在的種間、種內(nèi)關(guān)系,以及土壤因子間的互作,共同使退化草地植被與土壤因子的關(guān)系以復(fù)雜的動態(tài)變化形式存在。

    本研究基于冗余分析方法對高寒草甸退化草地植被與土壤因子的關(guān)系進(jìn)行分析,將退化草地研究樣地、樣地植被、樣地土壤因子有機(jī)結(jié)合,通過RDA二維排序圖直觀揭示樣地、物種、土壤因子間的關(guān)系,主要表現(xiàn)為:高寒草甸群落中,植物優(yōu)勢度在草地退化演替格局中發(fā)生變化,莎草科、禾本科植物優(yōu)勢度總體隨草地退化程度的加劇而降低,雜類草、豆科植物優(yōu)勢度總體隨退化程度的加劇而逐漸增加,植被在不同退化程度草地進(jìn)行的這種演替,使得高寒草甸的生產(chǎn)功能降低甚至喪失。利用退化草地植被—土壤因子數(shù)據(jù)矩陣進(jìn)行冗余分析能夠較好反映草地植被與土壤因子的關(guān)系。在RDA二維排序圖中,第1排序軸反映高寒草甸退化程度,沿第1排序軸草地退化程度加劇,土壤溫度升高、容質(zhì)量增大、有效氮、全氮、全磷含量降低、土壤含水量減少。第2排序軸則主要反映土壤機(jī)械組成的綜合變化。第1、第2排序軸能夠解釋74.5%的退化草地植被變化和83.2%的植被—土壤因子關(guān)系且反映出退化草地植被與土壤因子極顯著相關(guān)(P<0.01)。不同土壤因子在草地植被分布時的重要性不同,土壤有效氮、全氮等含量較高,土壤溫度、容質(zhì)量水平較低時,有利于高山嵩草、矮嵩草、垂穗披堿草等高寒草甸優(yōu)勢種植物的分布;雜類草及豆科植物在土壤溫度、容質(zhì)量等因子處于較高水平且土壤有效氮、土壤含水量等因子處于較低水平時更易在群落中分布。13個土壤因子的邊際作用檢驗也表明,土壤溫度、容質(zhì)量、有效氮、全氮、全磷等與退化草地植被間的關(guān)系更密切。

    冗余分析應(yīng)用于草地植被與土壤因子間、環(huán)境與土壤因子關(guān)系的研究已有報道。王興等[27]在荒漠草原棄耕地對草地土壤與植被間關(guān)系進(jìn)行RDA分析,認(rèn)為土壤碳酸鈣表征土壤沙化程度和退化階段,棄耕地植物分布格局是在土壤退化狀態(tài)約束下物種擴(kuò)散競爭的結(jié)果;白曉航等[28]將RDA分析用于亞高山草甸群落特征與環(huán)境因子間關(guān)系研究中,認(rèn)為群落分布格局與坡向、土壤溫度、海拔等因子極顯著相關(guān)(P<0.01),坡向、土壤溫度等因子對植被類型的形成具有關(guān)鍵作用;酈威等[29]對遼寧南部太子河流域河岸帶土壤理化特征與環(huán)境因子關(guān)系進(jìn)行冗余分析,認(rèn)為海拔、年降水量等是影響土壤空間差異的主導(dǎo)因子。由此看出,雖然不同研究人員在不同地區(qū)針對不同研究對象的研究結(jié)論不同,但是植被與土壤因子間的相關(guān)性均集中體現(xiàn)在土壤化學(xué)性質(zhì)、地形地貌、氣象因素等方面。而本研究在選取影響植被分布的因子時,側(cè)重選擇土壤物理、化學(xué)性質(zhì),而坡向、海拔、降水量等地形、氣象因素未納入分析中,應(yīng)在今后的研究中給予補(bǔ)充。

    表5 高寒草甸退化草地植被與土壤因子相關(guān)性Table 5 Correlation between vegetation of degraded alpine grassland and soil factors

    注:++表示極顯著正相關(guān)(P<0.01),+表示顯著正相關(guān)(P<0.05),--表示極顯著負(fù)相關(guān)(P<0.01),-表示顯著負(fù)相關(guān)(P<0.05)。

    Note: ++ Very significantly positive correlation(P<0.01),+ significantly positive correlation(P<0.05), --very significantly negative correlation(P<0.01), -significantly negative correlation(P<0.05).

    冗余分析應(yīng)用于植被與環(huán)境間關(guān)系的研究也存在不足。高寒草甸植被退化與土壤退化互為因果,植被退化與土壤退化耦合導(dǎo)致草地退化。在高寒草甸退化進(jìn)程中植被退化與土壤退化存在時間差,由于土壤具有相對較強(qiáng)的抗退化能力,草地植被異質(zhì)化出現(xiàn)一段時間后,才引起土壤中某些性質(zhì)的異質(zhì)化。但是在退化高寒草甸應(yīng)用冗余分析時,不能反映不同退化程度下植被與土壤因子的時間順序問題。

    目前,關(guān)于植被與土壤因子間關(guān)系研究更多的是探討植被群落多樣性、群落特性等與土壤環(huán)境因子間的相關(guān)性[30-33]。趙景學(xué)等[34]研究藏北高寒植被群落物種多樣性與土壤環(huán)境因子之間的關(guān)系,發(fā)現(xiàn)高寒草地土壤有機(jī)質(zhì)、全氮、有效氮、全磷與多樣性指數(shù)及豐富度指數(shù)呈極顯著正相關(guān)(P<0.01);保婭等[35]研究高山嵩草群落生物多樣性與土壤養(yǎng)分的相關(guān)性,發(fā)現(xiàn)土壤全氮、全鉀、pH對Simpson 指數(shù)影響較大,全氮、pH對Shannon-Wiener指數(shù)影響較大,說明土壤全鉀、pH影響植物物種在群落中的分布,與本研究得到的植被與pH、全鉀間未表現(xiàn)顯著相關(guān)的結(jié)論不一致,可能的原因在于高山嵩草草甸空間分布范圍廣、高山草甸土空間異質(zhì)性大、樣地退化程度不一致所導(dǎo)致。羅亞勇等[36]也以高寒草甸退化草地為對象,研究草甸退化過程中植物群落物種多樣性、生產(chǎn)力與土壤特性的關(guān)系,結(jié)果表明,群落物種多樣性和生產(chǎn)力與土壤有機(jī)碳、全氮、土壤含水量等呈極顯著正相關(guān)(P<0.01),這一結(jié)論中認(rèn)為群落生產(chǎn)力與諸多土壤因子正相關(guān),與本研究得到的莎草科、禾本科植物的分布與有機(jī)碳、全氮、土壤含水量等因子間極顯著正相關(guān)(P<0.01)的結(jié)論實質(zhì)是一致的。

    目前,在高寒草甸退化草地植物物種與土壤環(huán)境因子研究領(lǐng)域,尚未見相關(guān)文獻(xiàn)報道,但有研究人員在亞熱帶稻田棄耕濕地方面開展類似研究,反映出丁香蓼(Ludwigiaprostrata)—水竹葉(Murdanniatriquetra)物種組主要分布在土壤有效鉀含量較高區(qū)域,小白酒草(Conyzacanadensis)—雙穗雀稗(Paspalumpaspaloides)物種組的分布則與土壤有效磷含量、pH正相關(guān)[37]。

    本研究發(fā)現(xiàn),莎草科、禾本科植物與土壤全氮、土壤含水量等因子極顯著正相關(guān)(P<0.01),與容質(zhì)量、土壤溫度間極顯著負(fù)相關(guān)(P<0.01),雜類草中的很多物種與全氮、土壤含水量間極顯著負(fù)相關(guān)(P<0.01),而與容質(zhì)量、土壤溫度間極顯著正相關(guān)(P<0.01)。這一研究結(jié)果僅僅表明在土壤全氮、土壤含水量等因子相對較高,土壤容質(zhì)量等相對較低時更有利于莎草科、禾本科植物的分布與生長;而當(dāng)土壤全氮、土壤含水量等因子含量相對較低,土壤容質(zhì)量相對較高時對雜類草分布和生長的影響要小于莎草科等植物,并不代表雜類草就適合在較低氮磷含量條件下分布與生長。從高寒草甸退化實際看,當(dāng)草地土壤養(yǎng)分匱乏、土壤退化嚴(yán)重時,草地植物群落中所有植物都受到極為不利的影響,植被總蓋度往往低于20%[38],群落地上部生物量大幅度減少,而這種狀態(tài)在添加氮、磷等養(yǎng)分后可以得到改善[39]。

    高寒草甸生態(tài)系統(tǒng)在青藏高原空間異質(zhì)性高,草地退化情況以及土壤因子變化極其復(fù)雜,在今后開展關(guān)于高寒草地植被與土壤因子相關(guān)性研究時,應(yīng)以更多的研究樣地為樣本,納入更多的土壤因子指標(biāo),如土壤微生物、土壤活性有機(jī)碳、土壤基礎(chǔ)呼吸等,以期在更廣的維度反映自然規(guī)律,深刻揭示高寒草甸植物與土壤因子間的關(guān)系,揭示草地退化時植被退化、土壤退化的發(fā)生機(jī)理,為研究和治理高寒草甸生態(tài)系統(tǒng)問題提供科學(xué)依據(jù)。

    參考文獻(xiàn)Reference:

    [1] 周 麗,張德罡,贠旭江,等.退化高寒草甸植被與土壤特征[J].草業(yè)科學(xué),2016,33(11):2196-2201.

    ZHOU L,ZHANG D G,YUN X J,etal.The vegetation and soil characteristics of degraded alpine meadow[J].PataculturalScience,2016,33(11):2196-2201.

    [2] 張法偉,王軍邦,林 麗,等.青藏高原高寒嵩草草甸植被群落特征對退化演替的響應(yīng)[J].中國農(nóng)業(yè)氣象,2014,35(5):504-510.

    ZHANG F W,WANG J B,LIN L,etal.Response of plant community of alpineKobresiameadow to degradation succession in Qinghai-Tibetan plateau[J].ChineseJournalofAgrometeorology,2014,35(5):504-510.

    [3] 林 麗,李以康,張法偉,等.高寒矮嵩草群落退化演替系列氮磷生態(tài)化學(xué)計量學(xué)特征[J].生態(tài)學(xué)報,2013,33(17):5245-5251.

    LIN L,LI Y K,ZHANG F W,etal.Soil nitrogen and phosphorus stoichiometry in a degradation series ofKobresiahumulismeadows in the Tibetan Plateau[J].ActaEcologicaSinica,2013,33(17):5245-5251.

    [4] 宗 寧,石培禮,蔣 婧,等.淺耕對西藏高原退化草甸土壤和植物群落特征的影響[J].草業(yè)科學(xué),2014,31(1):8-14.

    ZONG N,SHI P L,JIANG J,etal.Effects of shallow plowing on soil and plant community of degraded alpine meadow in Tibetan Plateau[J].PrataculturalScience,2014,31(1):8-14.

    [5] 石紅霄,侯向陽,師尚禮,等.高山嵩草草甸初級生產(chǎn)力、多樣性與土壤因子的關(guān)系[J].草業(yè)學(xué)報,2015,24(10):40-47.

    SHI H X,HOU X Y,SHI SH L,etal.Relationships between plant diversity,soil property and productivity in an alpine meadow[J].ActaPrataculturaeSinica,2015,24(10):40-47.

    [6] 張生楹,張德罡,柳小妮,等.東祁連山不同退化程度高寒草甸土壤養(yǎng)分特征研究[J].草業(yè)科學(xué),2012,29(7):1028-1032.

    ZHANG SH Y,ZHANG D G,LIU X N,etal.Soil nutrient characteristics of alpine meadow at different degradation degrees in eastern Qilian Mountains[J].PrataculturalScience,2012,29(7):1028-1032.

    [7] 陳淑燕,張德罡.不同退化階段高寒草甸草地土壤鉀素的變化分析[J].草原與草坪,2013,33(3):74-77.

    CHEN SH Y,ZHANG D G.Potassium variation in alpine meadow soil at different degeneration stages[J].GrasslandandTurf,2013,33(3):74-77.

    [8] 羅亞勇,張 宇,張靜輝,等.不同退化階段高寒草甸土壤化學(xué)計量特征[J].生態(tài)學(xué)雜志,2012,31(2):254-260.

    LUO Y Y,ZHANG Y,ZHANG J H,etal.Soil stoichiometry characteristics of alpine meadow at its different degradation stages[J].ChineseJournalofEcology,2012,31(2):254-260.

    [9] 劉淑麗,林 麗,張法偉,等.放牧季節(jié)及退化程度對高寒草甸土壤有機(jī)碳的影響[J].草業(yè)科學(xué),2016,33(1):11-18.

    LIU SH L,LIN L,ZHANG F W,etal.Effects of grazing season and degradation degree on the soil organic carbon in alpine meadow[J].PrataculturalScience,2016,33(1):11-18.

    [10] 劉育紅,魏衛(wèi)東,溫小成,等.退化高寒草甸土壤有機(jī)碳組分特征[J].西北農(nóng)業(yè)學(xué)報,2015,24(2):168-174.

    LIU Y H,WEI W D,WEN X CH,etal.Characteristics of soil organic carbon fractions in alpine meadow with different degradation[J].ActaAgriculturaeBoreali-occidentalisSinica,2015,24(2):168-174.

    [11] 劉淑麗,林 麗,杜巖功,等.青海省高寒草甸不同退化階段土壤無機(jī)碳分異特征[J].生態(tài)學(xué)雜志,2014,33(5):1290-1296.

    LIU SH L,LIN L,DU Y G,etal.The variation of soil inorganic carbon of alpine meadow at different degeneration stages in Qinghai province[J].ChineseJournalofEcology,2014,33(5):1290-1296.

    [12] 林 麗,張法偉,李以康,等.高寒矮嵩草草甸退化過程土壤碳氮儲量及C/N化學(xué)計量學(xué)特征[J].中國草地學(xué)報,2012,34(3):42-47.

    LIN L,ZHANG F W,LI Y K,etal.The soil carbon and nitrogen storage and C/N metrological characteristics of chemistry inKobresiahumilismeadow in degradation succession stages[J].ChineseJournalofGrassland,2012,34(3):42-47.

    [13] 魏衛(wèi)東,劉育紅.不同退化程度高寒草地土壤微生物量碳特征分析[J].西北農(nóng)業(yè)學(xué)報,2014,23(2):205-210.

    WEI W D,LIU Y H.Characteristics analysis of soil microbial biomass carbon on different degraded degree on alpine grassland[J].ActaAgriculturaeBoreali-occidentalisSinica,2014,23(2):205-210.

    [14] 孫 哲,王一博,劉國華,等.基于多重分形理論的多年凍土區(qū)高寒草甸退化過程中土壤粒徑分析[J].冰川凍土,2015,37(4):980-990.

    SUN ZH,WANG Y B,LIU G H,etal.Heterogeneity analysis of soil particle size distribution in the process of degradation of alpine meadow in the permafrost regions based on multifractal theory[J].JournalofGlaciologyandGeocryology,2015,37(4):980-990.

    [15] 魏茂宏,林慧龍.江河源區(qū)高寒草甸退化序列土壤粒徑分布及其分形維數(shù)[J].應(yīng)用生態(tài)學(xué),2014,25(3):679-686.

    WEI M H,LIN H L.Soil particle size distribution and its fractal dimension among degradation sequences of the alpine meadow in the source region of the Yangtze and Yellow River,Qinghai-Tibetan Plateau,China[J].ChineseJournalofAppliedEcology,2014,25(3):679-686.

    [16] 尤全剛,薛 嫻,彭 飛,等.高寒草甸草地退化對土壤水熱性質(zhì)的影響及其環(huán)境效應(yīng)[J].中國沙漠,2015,35(5):1183-1192.

    YOU Q G,XUE L,PENG F,etal.Alpine meadow degradation effect on soil thermal and hydraulic properties and its environmental impacts[J].JournalofDesertResearch,2015,35(5):1183-1192.

    [17] 張 泉,劉詠梅,楊勤科,等.祁連山退化高寒草甸土壤水分空間變異特征分析[J].冰川凍土,2014,36(1):88-94.

    ZHANG Q,LIU Y M,YANG Q K,etal.Analysis of the spatial variability of soil moisture in degrading alpine meadow in the Qilian Mountains[J].JournalofGlaciologyandGeocryology,2014,36(1):88-94.

    [18] 胡 雷,王長庭,王根緒,等.三江源區(qū)不同退化演替階段高寒草甸土壤酶活性和微生物群落結(jié)構(gòu)的變化[J].草業(yè)學(xué)報,2014,23(3):8-19.

    HU L,WANG CH T,WANG G X,etal.Changes in the activities of soil enzymes and microbial community structure at different degradation successional stages of alpine meadows in the headwater region of Three Rivers,China[J].ActaPrataculturaeSinica,2014,23(3):8-19.

    [19] 魏衛(wèi)東,劉育紅.三江源區(qū)高寒草甸退化對土壤呼吸的影響[J].黑龍江畜牧獸醫(yī),2014(19):102-105.

    WEI W D,LIU Y H.Effects of alpine meadow degradation on soil respiration in the source area of Three Major Rivers of China [J].HeilongjiangAnimalScienceandVeterinaryMedicine,2014(19):102-105.

    [20] 任繼周.草業(yè)科學(xué)研究方法[M].北京:中國農(nóng)業(yè)出版社,1998:95.

    REN J ZH.Research Methods of Pratacultural Science[M].Beijing:China Agriculture Press,1998:3-15.

    [21] 趙新全.三江源區(qū)退化草地生態(tài)系統(tǒng)恢復(fù)與可持續(xù)管理[M].北京:科學(xué)出版社,2011:127-128.

    ZHAO X Q.Ecosystem Restoration and Sustainable Management of Degraded Grassland in the Source Area of Three Major Rivers of China [M].Beijing:Science Press,2011:127-128.

    [22] 張甘霖,龔子同.土壤調(diào)查實驗室分析方法[M].北京:科學(xué)出版社,2012:38-80.

    ZHANG G L,GONG Z T.Soil Survey Laboratory Methods[M].Beijing:Science Press,2012:38-80.

    [23] 鮑士旦.土壤農(nóng)化分析[M].北京:中國農(nóng)業(yè)出版社,1999:30-38.

    BAO SH D.Soil and Agricultural Chemistry Analysis[M].Beijing:China Agriculture Press,1999:30-38.

    [24] 陳麗瓊.比重計法測定土壤顆粒組成的研究[J].環(huán)境科學(xué)導(dǎo)刊,2010,29(4):97-99.

    CHEN L Q.Research on structure of soil particle by hydrometer method[J].EnvironmentalScienceSurvey,2010,29(4):97-99.

    [25] FRANKLIN J,WISER S K,DRAKE D,etal.Environment,disturbance history and rain forest composition across the islands of Tonga,Western Polynesia [J].JournalofVegetationScience,2006,17:233-244.

    [26] LEGENDRE P,LEGENDRE L.Numerical Ecology [M].Amsterdam:Elsevier,2012:126.

    [27] 王 興,宋乃平,楊新國,等.荒漠草原棄耕地恢復(fù)草地土壤與植被的RDA分析[J].草業(yè)學(xué)報,2014,23(2):90-97.

    WANG X,SONG N P,YNAG X G,etal.Redundancy analysis of soil vegetation of recovered grassland on abandoned land in the desert steppe[J].ActaPrataculturaeSinica,2014,23(2):90-97.

    [28] 白曉航,張金屯,曹 科,等.小五臺山亞高山草甸的群落特征及物種多樣性[J].草業(yè)科學(xué),2016,33(12):2538-2543.

    BAI X H,ZHANG J T,CAO K,etal.Community characteristics and species diversity of subalpine meadows in Xiaowutai Mountain[J].PrataculturalScience,2016,33(12):2538-2543.

    [29] 酈 威,盧振蘭,孔維靜,等.太子河河岸帶土壤理化性質(zhì)特征及其與環(huán)境因子的關(guān)系[J].生態(tài)科學(xué),2013,32(1):90-97.

    LI W,LU ZH L,KONG W J,etal.Physical and chemical properties of riparian soils in the Taizi river basin and their relationship with environmental factors[J].EcologicalScience,2013,32(1):90-97.

    [30] 李立新,陳英智.退化演替對高山草地植被和土壤理化特性影響[J].水土保持通報,2015,35(4):40-43.

    LI L X,CHEN Y ZH.Effect of degradation succession on vegetation and soil physical and chemical characteristics in alpine grassland [J].BulletinofSoilandWaterConservation,2015,35(4):40-43.

    [31] 許加星,徐力剛,姜加虎.鄱陽湖典型洲灘植物群落結(jié)構(gòu)變化及其與土壤養(yǎng)分的關(guān)系[J].濕地科學(xué),2013,11(2):186-191.

    XU J X,XU L G,JIANG J H.Change of vegetation community structure and the relationship between it and soil nutrients in typical beaches in Poyang lake area[J].WetlandScience,2013,11(2):186-191.

    [32] 張發(fā)偉,郭竹筠,李以康,等.青海湖芨芨草干草原植被退化對土壤溫濕特征的影響[J].干旱區(qū)研究,2013,30(2):219-225.

    ZHANG F W,GUO ZH J,LI Y K,etal.Effect of soil temperature and volumetric moisture content on vegetation degeneration over theAchnatherumsplendenssteppe around the Qinghai Lake[J].AridZoneResearch,2013,30(2):219-225.

    [33] 劉 蓉,張衛(wèi)國,江小雷,等.垂穗披堿草群落退化演替的植被特性及其與土壤性狀的相關(guān)性研究[J].草業(yè)科學(xué),2010,27(10):96-103.

    LIU R,ZHANG W G,JIANG X L,etal.Study on the characteristics of degradation succession ofElymusnutanscommunity and its correlation to soil properties[J].PrataculturalScience,2010,27(10):96-103.

    [34] 趙景學(xué),曲廣鵬,多吉頓珠,等.藏北高寒植被群落物種多樣性與土壤環(huán)境因子的關(guān)系[J].干旱區(qū)資源與環(huán)境,2011,25(6):105-108.

    ZHAO J X,QU G P,DUO JIDUNZHU,etal.Relationship between species diversity and soil factors of alpine grasslands in north Tibet[J].JournalofAridLandResourcesandEnvironment,2011,25(6):105-108.

    [35] 保 婭,達(dá)哇卓瑪.高山嵩草群落生物多樣性與土壤養(yǎng)分的相關(guān)性研究[J].青海草業(yè),2015,24(2):12-15.

    BAO Y,DAWAZHUMA.Study on relationship betweenKobresiapygmaeacommunity biodiversity and soil nutrition[J].QinghaiPratacultural,2015,24(2):12-15.

    [36] 羅亞勇,孟慶濤,張靜輝,等.青藏高原東緣高寒草甸退化過程中植物群落物種多樣性、生產(chǎn)力與土壤特性的關(guān)系[J].冰川凍土,2014,36(5):1298-1305.

    LUO Y Y,MENG Q T,ZHANG J H,etal.Species diversity and biomass in relation to soil properties of alpine meadows in the eastern Tibetan Plateau in different degradation stages[J].JournalofGlaciologyandGeocryology,2014,36(5):1298-1305.

    [37] 彭 億,李裕元,李忠武,等.亞熱帶稻田棄耕濕地土壤因子對植物群落結(jié)構(gòu)的影響[J].應(yīng)用生態(tài)學(xué)報,2009,20(7):1543-1550.

    PENG Y,LI Y Y,LI ZH W,etal.Effects of soil factors on vegetation community structure in an abandoned subtropical paddy wetland[J].ChineseJournalofAppliedEcology,2009,20(7):1543-1550.

    [38] 辛有俊,杜鐵瑛.青海天然草地退化及恢復(fù)研究[M].西寧:青海人民出版社,2013,53.

    XIN Y J,DU T Y.Natural Grassland Degradation and Its Rehabilitation in Qinghai[M].Xining:Qinghai People Press,2013:53.

    [39] 楊元武,周華坤,李希來,等.高寒草甸物種多樣性和生產(chǎn)力對樣方添加的初期響應(yīng)[J].西北農(nóng)業(yè)學(xué)報,2017,26(2):159-166.

    YANG Y W,ZHOU H K,LI X L,etal.Initial response of species diversity and productivity to nutrients addition on alpine mesdow[J].ActaAgriculturaeBoreali-occidentalisSinica,2017,26(2):159-166.

    猜你喜歡
    草甸土壤溫度全氮
    高山草甸
    The world's narrowest river
    山地草甸
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    遼東山區(qū)3種人工林土壤呼吸對土壤溫度和土壤水分的響應(yīng)
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    豐鎮(zhèn)市農(nóng)田土壤有機(jī)質(zhì)與全氮含量關(guān)系分析
    不同土地利用方式對黒壚土有機(jī)質(zhì)和全氮分布規(guī)律的影響
    土壤與作物(2015年3期)2015-12-08 00:47:01
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    紅堿淖流域濕地土壤全氮含量及分布特征研究
    亚洲av.av天堂| 亚洲欧洲日产国产| 交换朋友夫妻互换小说| 久久精品国产a三级三级三级| 黄色怎么调成土黄色| 日韩av在线免费看完整版不卡| 免费少妇av软件| 熟女人妻精品中文字幕| 大码成人一级视频| 91午夜精品亚洲一区二区三区| 91精品一卡2卡3卡4卡| 国产精品久久久久久精品古装| 18禁在线播放成人免费| 国产 一区精品| 免费高清在线观看日韩| 亚洲av在线观看美女高潮| 一本大道久久a久久精品| 欧美精品国产亚洲| 久久热精品热| 国产又色又爽无遮挡免| 国产精品久久久久久久久免| 毛片一级片免费看久久久久| 亚州av有码| 欧美日韩一区二区视频在线观看视频在线| 中文字幕制服av| 国产av国产精品国产| 亚洲av男天堂| 亚洲av欧美aⅴ国产| 91精品三级在线观看| 国产精品麻豆人妻色哟哟久久| 中文字幕精品免费在线观看视频 | 亚洲第一区二区三区不卡| 精品人妻在线不人妻| 国产一级毛片在线| 一区二区三区免费毛片| 中文精品一卡2卡3卡4更新| 男女边摸边吃奶| 岛国毛片在线播放| 精品视频人人做人人爽| 麻豆乱淫一区二区| 国产成人免费无遮挡视频| 成人国语在线视频| 只有这里有精品99| 成人毛片60女人毛片免费| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三卡| 涩涩av久久男人的天堂| 欧美xxxx性猛交bbbb| 亚洲综合精品二区| 黄色怎么调成土黄色| 国产精品欧美亚洲77777| 精品国产一区二区久久| 免费人妻精品一区二区三区视频| 一级毛片黄色毛片免费观看视频| 另类亚洲欧美激情| 久久精品国产鲁丝片午夜精品| 三级国产精品片| 简卡轻食公司| 亚洲国产最新在线播放| 99九九在线精品视频| 最后的刺客免费高清国语| 国产男人的电影天堂91| 观看av在线不卡| 免费观看av网站的网址| 性高湖久久久久久久久免费观看| 在线播放无遮挡| 亚洲欧美色中文字幕在线| 如何舔出高潮| 国产黄色视频一区二区在线观看| 丝袜美足系列| 亚洲欧洲精品一区二区精品久久久 | 一本色道久久久久久精品综合| 亚洲色图综合在线观看| 亚洲天堂av无毛| 在线看a的网站| 国产精品人妻久久久久久| 男女无遮挡免费网站观看| 亚洲国产欧美日韩在线播放| a级毛片黄视频| 黄色毛片三级朝国网站| 国产精品三级大全| 亚洲天堂av无毛| 天天躁夜夜躁狠狠久久av| 国产欧美亚洲国产| 在线观看一区二区三区激情| 熟女人妻精品中文字幕| 一级二级三级毛片免费看| 日韩伦理黄色片| 亚洲精品aⅴ在线观看| 中国美白少妇内射xxxbb| 纯流量卡能插随身wifi吗| 美女cb高潮喷水在线观看| 激情五月婷婷亚洲| 26uuu在线亚洲综合色| 国产女主播在线喷水免费视频网站| 伊人久久精品亚洲午夜| 欧美精品高潮呻吟av久久| 妹子高潮喷水视频| 91精品国产国语对白视频| 各种免费的搞黄视频| 哪个播放器可以免费观看大片| 免费观看在线日韩| 亚洲欧美成人精品一区二区| 热re99久久精品国产66热6| 国语对白做爰xxxⅹ性视频网站| av又黄又爽大尺度在线免费看| 91午夜精品亚洲一区二区三区| 国产成人freesex在线| 国产亚洲精品第一综合不卡 | 国产亚洲最大av| 18在线观看网站| 国产男女超爽视频在线观看| 免费人成在线观看视频色| 老司机影院成人| 国产高清有码在线观看视频| 日日撸夜夜添| 亚洲一区二区三区欧美精品| 成人亚洲欧美一区二区av| 一区二区日韩欧美中文字幕 | 国产色婷婷99| 免费黄色在线免费观看| 伦理电影免费视频| 亚洲精品一区蜜桃| 国产高清不卡午夜福利| a级毛片在线看网站| 18禁裸乳无遮挡动漫免费视频| 黑丝袜美女国产一区| 免费看光身美女| 国产欧美亚洲国产| 大香蕉久久网| 成年人免费黄色播放视频| 91精品三级在线观看| 熟女av电影| 在线播放无遮挡| 国产伦理片在线播放av一区| 免费看不卡的av| 伦精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 国产精品不卡视频一区二区| 色5月婷婷丁香| 日日摸夜夜添夜夜添av毛片| 国产午夜精品一二区理论片| 国产精品.久久久| 国产在线免费精品| 亚洲性久久影院| 国产精品99久久久久久久久| 日日摸夜夜添夜夜添av毛片| 亚洲精华国产精华液的使用体验| 在线观看免费视频网站a站| 国产日韩欧美在线精品| 久久狼人影院| 久久久久久久亚洲中文字幕| 亚洲欧美色中文字幕在线| 乱码一卡2卡4卡精品| 高清在线视频一区二区三区| 亚洲av在线观看美女高潮| 久久久久久久国产电影| 蜜桃国产av成人99| 日本黄色片子视频| 亚洲,一卡二卡三卡| 9色porny在线观看| 成人18禁高潮啪啪吃奶动态图 | 天堂8中文在线网| 伦理电影免费视频| 夫妻午夜视频| 一级片'在线观看视频| 久久久久久久久大av| 亚洲伊人久久精品综合| 少妇人妻久久综合中文| 日韩伦理黄色片| 精品国产一区二区久久| 中国三级夫妇交换| 亚洲内射少妇av| 久久ye,这里只有精品| 国产综合精华液| 色视频在线一区二区三区| 久久青草综合色| videossex国产| av不卡在线播放| 精品人妻熟女av久视频| 少妇 在线观看| 99热全是精品| 国产精品一二三区在线看| 一边摸一边做爽爽视频免费| 少妇 在线观看| 热re99久久精品国产66热6| 99热这里只有是精品在线观看| 在线免费观看不下载黄p国产| 日韩一区二区三区影片| 国产精品一二三区在线看| 性色av一级| 国产精品久久久久久久久免| 女的被弄到高潮叫床怎么办| 亚洲国产欧美日韩在线播放| 亚洲不卡免费看| 黑人巨大精品欧美一区二区蜜桃 | 精品人妻熟女av久视频| 亚洲欧美中文字幕日韩二区| 国产精品嫩草影院av在线观看| 亚洲精品日本国产第一区| 在线 av 中文字幕| 日韩免费高清中文字幕av| 在线 av 中文字幕| 亚洲欧美中文字幕日韩二区| 在线观看www视频免费| 欧美日韩亚洲高清精品| 涩涩av久久男人的天堂| 日韩精品免费视频一区二区三区 | 亚洲欧洲精品一区二区精品久久久 | 国产日韩欧美视频二区| 丝袜在线中文字幕| 久久精品国产鲁丝片午夜精品| 天天影视国产精品| 一级毛片aaaaaa免费看小| 亚洲欧美成人精品一区二区| 亚洲中文av在线| 日本黄大片高清| 一级,二级,三级黄色视频| 哪个播放器可以免费观看大片| 亚洲伊人久久精品综合| 美女国产视频在线观看| 又粗又硬又长又爽又黄的视频| 91久久精品国产一区二区三区| 久久 成人 亚洲| 人妻人人澡人人爽人人| 亚洲精品日本国产第一区| 国产精品国产三级国产专区5o| 妹子高潮喷水视频| 久久久久久人妻| 大又大粗又爽又黄少妇毛片口| 久久国产精品男人的天堂亚洲 | 午夜精品国产一区二区电影| 狠狠精品人妻久久久久久综合| 丁香六月天网| 80岁老熟妇乱子伦牲交| 欧美激情 高清一区二区三区| 精品人妻熟女毛片av久久网站| 午夜福利网站1000一区二区三区| 午夜日本视频在线| freevideosex欧美| 亚洲熟女精品中文字幕| 精品久久久久久久久亚洲| 少妇人妻久久综合中文| 欧美另类一区| 卡戴珊不雅视频在线播放| 国产精品久久久久成人av| 日本黄色片子视频| 内地一区二区视频在线| 久久这里有精品视频免费| 少妇高潮的动态图| 欧美精品人与动牲交sv欧美| 日韩中字成人| 亚洲欧洲日产国产| 国产免费一区二区三区四区乱码| 亚洲国产精品成人久久小说| 女的被弄到高潮叫床怎么办| 看免费成人av毛片| 国产精品99久久久久久久久| 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 日韩制服骚丝袜av| 久久ye,这里只有精品| 一区二区三区乱码不卡18| 日韩一区二区三区影片| 精品亚洲成a人片在线观看| 美女中出高潮动态图| www.av在线官网国产| 男女免费视频国产| 免费播放大片免费观看视频在线观看| 国国产精品蜜臀av免费| 91精品伊人久久大香线蕉| 日韩 亚洲 欧美在线| 乱人伦中国视频| 91精品三级在线观看| 22中文网久久字幕| 97超视频在线观看视频| 一区二区三区四区激情视频| 亚洲欧美中文字幕日韩二区| 欧美日韩亚洲高清精品| 久久精品熟女亚洲av麻豆精品| 国产白丝娇喘喷水9色精品| 一区二区三区精品91| 国产精品人妻久久久影院| av国产精品久久久久影院| 狂野欧美白嫩少妇大欣赏| 在线观看人妻少妇| 久久久久国产网址| 日韩av在线免费看完整版不卡| 日韩不卡一区二区三区视频在线| 亚洲精品国产色婷婷电影| 精品久久久精品久久久| 妹子高潮喷水视频| 欧美97在线视频| 黄色一级大片看看| 91国产中文字幕| 午夜久久久在线观看| 我要看黄色一级片免费的| 免费黄网站久久成人精品| 99re6热这里在线精品视频| 亚洲精品日韩在线中文字幕| 国产成人精品婷婷| 春色校园在线视频观看| 国产欧美日韩一区二区三区在线 | 日韩三级伦理在线观看| 亚洲成人手机| 热re99久久精品国产66热6| 日韩一区二区视频免费看| 国产男女超爽视频在线观看| 少妇丰满av| 久久精品久久精品一区二区三区| 99热全是精品| 另类亚洲欧美激情| 久久免费观看电影| 国产亚洲午夜精品一区二区久久| 2021少妇久久久久久久久久久| 精品人妻一区二区三区麻豆| 欧美最新免费一区二区三区| 又大又黄又爽视频免费| 国产黄片视频在线免费观看| 另类精品久久| 色5月婷婷丁香| 男女边摸边吃奶| 黑人猛操日本美女一级片| 人妻少妇偷人精品九色| 精品视频人人做人人爽| 午夜福利在线观看免费完整高清在| 美女国产视频在线观看| 美女福利国产在线| 亚洲精品456在线播放app| 欧美老熟妇乱子伦牲交| 国产精品 国内视频| 国产精品人妻久久久影院| 一区二区三区免费毛片| av线在线观看网站| 成人免费观看视频高清| 美女主播在线视频| 校园人妻丝袜中文字幕| 秋霞在线观看毛片| 80岁老熟妇乱子伦牲交| 国产成人精品无人区| 各种免费的搞黄视频| 国语对白做爰xxxⅹ性视频网站| 三上悠亚av全集在线观看| 亚洲性久久影院| 亚洲精品久久久久久婷婷小说| 国产亚洲最大av| 国产精品久久久久久精品电影小说| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花| 一区二区av电影网| 亚洲成人手机| av黄色大香蕉| 99精国产麻豆久久婷婷| 色吧在线观看| 欧美 亚洲 国产 日韩一| 国产精品麻豆人妻色哟哟久久| 免费观看在线日韩| 亚洲国产欧美在线一区| 99久国产av精品国产电影| kizo精华| 欧美人与性动交α欧美精品济南到 | 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 看十八女毛片水多多多| 黄片播放在线免费| 国产精品国产三级国产av玫瑰| 国产亚洲最大av| 热99国产精品久久久久久7| 久久精品国产鲁丝片午夜精品| 亚洲av二区三区四区| 国产极品粉嫩免费观看在线 | 美女国产高潮福利片在线看| 在线观看一区二区三区激情| 国产白丝娇喘喷水9色精品| 成年美女黄网站色视频大全免费 | 美女国产高潮福利片在线看| 亚洲av成人精品一区久久| 人体艺术视频欧美日本| 日本wwww免费看| av卡一久久| 亚洲av在线观看美女高潮| videosex国产| 看非洲黑人一级黄片| 午夜免费鲁丝| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 欧美日韩在线观看h| 久久久精品94久久精品| 亚洲精品第二区| 国产精品人妻久久久久久| 国产一区二区三区av在线| 日本爱情动作片www.在线观看| 亚洲人成网站在线播| www.av在线官网国产| 日本色播在线视频| 亚洲,一卡二卡三卡| 久久久久久人妻| 性高湖久久久久久久久免费观看| 3wmmmm亚洲av在线观看| 久久国产亚洲av麻豆专区| 自线自在国产av| 欧美最新免费一区二区三区| 亚洲怡红院男人天堂| 大片免费播放器 马上看| 纯流量卡能插随身wifi吗| 看十八女毛片水多多多| 少妇人妻久久综合中文| 人妻人人澡人人爽人人| 在线观看一区二区三区激情| 成人毛片60女人毛片免费| 狂野欧美激情性xxxx在线观看| 国产视频内射| 国产在线免费精品| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费高清a一片| 亚洲少妇的诱惑av| 国产成人精品久久久久久| 男女高潮啪啪啪动态图| 麻豆乱淫一区二区| 精品少妇黑人巨大在线播放| 大香蕉久久成人网| 青青草视频在线视频观看| 欧美日韩亚洲高清精品| 久久国产精品大桥未久av| 亚洲av福利一区| 99热6这里只有精品| 伦精品一区二区三区| 欧美日韩视频精品一区| 高清不卡的av网站| 精品人妻熟女av久视频| 日本91视频免费播放| videosex国产| 国产黄频视频在线观看| 人人妻人人添人人爽欧美一区卜| 欧美成人午夜免费资源| 国产免费福利视频在线观看| 国产高清国产精品国产三级| 成人二区视频| 色94色欧美一区二区| 日韩av免费高清视频| 欧美丝袜亚洲另类| 黑丝袜美女国产一区| 亚洲综合色网址| 久久99热6这里只有精品| 亚洲av二区三区四区| 少妇猛男粗大的猛烈进出视频| 婷婷成人精品国产| 欧美最新免费一区二区三区| 亚洲内射少妇av| 建设人人有责人人尽责人人享有的| 女性生殖器流出的白浆| xxx大片免费视频| 建设人人有责人人尽责人人享有的| 亚洲久久久国产精品| 亚洲怡红院男人天堂| 在现免费观看毛片| 97超碰精品成人国产| 少妇人妻 视频| 高清欧美精品videossex| 亚洲,欧美,日韩| 少妇人妻 视频| 亚洲欧美中文字幕日韩二区| 国产一区二区三区综合在线观看 | 又大又黄又爽视频免费| 国内精品宾馆在线| 秋霞伦理黄片| 一区二区av电影网| 简卡轻食公司| 少妇高潮的动态图| 国产男女内射视频| 欧美成人精品欧美一级黄| 免费高清在线观看日韩| 有码 亚洲区| 久久久久久人妻| 一级,二级,三级黄色视频| 校园人妻丝袜中文字幕| 在线观看一区二区三区激情| 日韩伦理黄色片| 国产视频内射| 欧美最新免费一区二区三区| a级毛片免费高清观看在线播放| 亚洲精品国产av蜜桃| 三级国产精品欧美在线观看| 亚洲精品久久午夜乱码| 日本欧美视频一区| 69精品国产乱码久久久| 五月天丁香电影| 少妇高潮的动态图| 一级毛片电影观看| 老熟女久久久| 亚洲激情五月婷婷啪啪| 精品国产国语对白av| 国产精品不卡视频一区二区| 色婷婷av一区二区三区视频| 51国产日韩欧美| videosex国产| av天堂久久9| 91成人精品电影| av播播在线观看一区| 夜夜骑夜夜射夜夜干| 91精品三级在线观看| 亚洲av福利一区| 国产精品久久久久久久久免| 纵有疾风起免费观看全集完整版| 纯流量卡能插随身wifi吗| 三级国产精品欧美在线观看| 91久久精品国产一区二区成人| 日韩在线高清观看一区二区三区| 色吧在线观看| 亚洲精品国产av蜜桃| 男女啪啪激烈高潮av片| videossex国产| 午夜福利,免费看| 国产一区亚洲一区在线观看| 妹子高潮喷水视频| 美女大奶头黄色视频| 欧美日韩国产mv在线观看视频| 亚洲人成网站在线观看播放| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 精品人妻偷拍中文字幕| 午夜福利视频精品| 亚洲人成网站在线观看播放| 久久人人爽人人爽人人片va| 亚洲少妇的诱惑av| 久久精品国产自在天天线| 一级黄片播放器| 国产成人精品久久久久久| 毛片一级片免费看久久久久| 亚洲久久久国产精品| 欧美国产精品一级二级三级| 乱码一卡2卡4卡精品| 免费播放大片免费观看视频在线观看| 亚洲精品成人av观看孕妇| 天天操日日干夜夜撸| 国产在视频线精品| 欧美最新免费一区二区三区| 欧美日韩视频高清一区二区三区二| 久久婷婷青草| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 国产 一区精品| 99热这里只有是精品在线观看| 人成视频在线观看免费观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| 精品少妇久久久久久888优播| av视频免费观看在线观看| 色吧在线观看| 国产在线一区二区三区精| 精品国产乱码久久久久久小说| 欧美3d第一页| 国产精品国产三级专区第一集| 国产精品久久久久成人av| 久久国产精品大桥未久av| 久久青草综合色| 午夜福利,免费看| 久久人人爽av亚洲精品天堂| 97超碰精品成人国产| av一本久久久久| 热99久久久久精品小说推荐| 97在线视频观看| 日韩视频在线欧美| 国产成人精品婷婷| 精品一区二区三卡| xxx大片免费视频| 一级爰片在线观看| 大码成人一级视频| 亚洲av成人精品一区久久| 日本vs欧美在线观看视频| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 亚洲精品456在线播放app| 成人手机av| 欧美3d第一页| 亚洲天堂av无毛| 飞空精品影院首页| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久小说| 纯流量卡能插随身wifi吗| 精品国产乱码久久久久久小说| 嘟嘟电影网在线观看| 亚洲熟女精品中文字幕| 欧美 日韩 精品 国产| 看非洲黑人一级黄片| 成年av动漫网址| 国产国语露脸激情在线看| 黑人巨大精品欧美一区二区蜜桃 | 午夜av观看不卡| 亚洲av不卡在线观看| 久久久久精品久久久久真实原创| 九九爱精品视频在线观看| 精品久久久久久电影网| 国产深夜福利视频在线观看| www.色视频.com| 国产一区有黄有色的免费视频| 啦啦啦在线观看免费高清www| 大又大粗又爽又黄少妇毛片口| 少妇的逼好多水| 国产av一区二区精品久久| 亚洲精品国产色婷婷电影| 久久亚洲国产成人精品v| 人人妻人人添人人爽欧美一区卜| 国产探花极品一区二区| 日韩成人av中文字幕在线观看| 国产在线视频一区二区| 国产av国产精品国产| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图 | 一区二区三区四区激情视频| 日韩一区二区视频免费看| 日产精品乱码卡一卡2卡三| 一本一本综合久久| 精品酒店卫生间|