• <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久久| 亚洲av在线观看美女高潮| 亚洲第一青青草原| av在线观看视频网站免费| 啦啦啦在线观看免费高清www| 99九九在线精品视频| 国产成人欧美| 九草在线视频观看| 日本一区二区免费在线视频| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| kizo精华| 亚洲精品日韩在线中文字幕| 久久ye,这里只有精品| a级毛片黄视频| 亚洲成人免费av在线播放| 色94色欧美一区二区| 日本色播在线视频| 丁香六月欧美| 乱人伦中国视频| 伦理电影大哥的女人| 一区在线观看完整版| 青春草亚洲视频在线观看| 777久久人妻少妇嫩草av网站| tube8黄色片| 成人亚洲欧美一区二区av| 又大又黄又爽视频免费| 少妇猛男粗大的猛烈进出视频| 亚洲精品中文字幕在线视频| 女人久久www免费人成看片| av免费观看日本| 天堂8中文在线网| 亚洲三区欧美一区| 老汉色av国产亚洲站长工具| 少妇人妻 视频| 国精品久久久久久国模美| 国产成人欧美| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩综合在线一区二区| 欧美日韩国产mv在线观看视频| 成年av动漫网址| 国产一区二区激情短视频 | 一级片免费观看大全| 亚洲成国产人片在线观看| 中文字幕最新亚洲高清| 亚洲欧洲日产国产| 99久久精品国产亚洲精品| 香蕉丝袜av| 叶爱在线成人免费视频播放| 在线精品无人区一区二区三| 成年动漫av网址| 欧美av亚洲av综合av国产av | 成人午夜精彩视频在线观看| 亚洲成人免费av在线播放| 国产精品香港三级国产av潘金莲 | 在线天堂中文资源库| 久久精品人人爽人人爽视色| 亚洲成人av在线免费| 男女下面插进去视频免费观看| 纵有疾风起免费观看全集完整版| 免费高清在线观看视频在线观看| 两个人看的免费小视频| 少妇猛男粗大的猛烈进出视频| 国产片特级美女逼逼视频| 另类精品久久| 蜜桃在线观看..| 国产成人精品无人区| 午夜日本视频在线| 久久精品久久久久久噜噜老黄| 黄片播放在线免费| 黄色 视频免费看| av在线老鸭窝| 亚洲国产精品国产精品| 女人爽到高潮嗷嗷叫在线视频| 日本欧美国产在线视频| 日韩免费高清中文字幕av| 高清不卡的av网站| 丁香六月欧美| 欧美乱码精品一区二区三区| 久久国产亚洲av麻豆专区| 日韩不卡一区二区三区视频在线| 男的添女的下面高潮视频| 宅男免费午夜| 日本黄色日本黄色录像| av.在线天堂| 只有这里有精品99| 亚洲av男天堂| 在线观看国产h片| 亚洲欧美一区二区三区久久| 国产精品国产三级国产专区5o| 狠狠精品人妻久久久久久综合| 欧美97在线视频| 久久久国产精品麻豆| 欧美亚洲 丝袜 人妻 在线| 亚洲精品乱久久久久久| netflix在线观看网站| 国产亚洲精品第一综合不卡| 电影成人av| 青春草亚洲视频在线观看| 国产片特级美女逼逼视频| 久久久久久免费高清国产稀缺| 国产成人欧美| 亚洲国产欧美日韩在线播放| 男人舔女人的私密视频| 丝瓜视频免费看黄片| 乱人伦中国视频| 午夜福利视频在线观看免费| 如何舔出高潮| 中文字幕亚洲精品专区| 亚洲图色成人| 免费看不卡的av| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩成人在线一区二区| 亚洲国产欧美在线一区| 91精品国产国语对白视频| 久热爱精品视频在线9| 国产亚洲av片在线观看秒播厂| 男女床上黄色一级片免费看| 国产成人精品久久二区二区91 | 国产亚洲欧美精品永久| 欧美乱码精品一区二区三区| 欧美日韩亚洲高清精品| 伦理电影大哥的女人| 国产97色在线日韩免费| 99热全是精品| 国产高清不卡午夜福利| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 精品国产乱码久久久久久男人| 亚洲精品国产一区二区精华液| 亚洲欧美激情在线| 久久女婷五月综合色啪小说| 如日韩欧美国产精品一区二区三区| 精品久久久久久电影网| 18禁裸乳无遮挡动漫免费视频| 亚洲国产av影院在线观看| av又黄又爽大尺度在线免费看| 亚洲国产中文字幕在线视频| 国产激情久久老熟女| 国产精品一区二区在线不卡| 久久99热这里只频精品6学生| 日日爽夜夜爽网站| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 高清视频免费观看一区二区| 永久免费av网站大全| 免费不卡黄色视频| 欧美日韩视频高清一区二区三区二| 欧美 亚洲 国产 日韩一| 亚洲国产精品国产精品| 女性生殖器流出的白浆| 久久精品亚洲av国产电影网| 国产激情久久老熟女| 欧美国产精品va在线观看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 黑丝袜美女国产一区| 国产一区二区在线观看av| 久久青草综合色| 一边摸一边做爽爽视频免费| 纯流量卡能插随身wifi吗| 飞空精品影院首页| 在线观看免费视频网站a站| 啦啦啦 在线观看视频| 亚洲av电影在线进入| 久久久亚洲精品成人影院| 制服诱惑二区| 各种免费的搞黄视频| 亚洲欧洲精品一区二区精品久久久 | 777米奇影视久久| 国产成人免费无遮挡视频| 麻豆乱淫一区二区| 久久狼人影院| 男女之事视频高清在线观看 | 国产一卡二卡三卡精品 | 中国三级夫妇交换| 美国免费a级毛片| 老熟女久久久| 国产高清不卡午夜福利| 精品国产一区二区久久| 丝袜美足系列| 午夜影院在线不卡| 免费高清在线观看日韩| 满18在线观看网站| 制服丝袜香蕉在线| 亚洲欧洲国产日韩| 亚洲欧洲精品一区二区精品久久久 | 婷婷色麻豆天堂久久| 看非洲黑人一级黄片| 欧美另类一区| 五月天丁香电影| 黑人巨大精品欧美一区二区蜜桃| videos熟女内射| 性高湖久久久久久久久免费观看| 午夜日本视频在线| www.熟女人妻精品国产| 下体分泌物呈黄色| 午夜福利乱码中文字幕| 母亲3免费完整高清在线观看| 777米奇影视久久| a级毛片黄视频| 日韩,欧美,国产一区二区三区| 最新在线观看一区二区三区 | 亚洲精品国产区一区二| 大话2 男鬼变身卡| 国精品久久久久久国模美| 美女中出高潮动态图| 欧美日韩视频精品一区| 国产在视频线精品| 国产精品.久久久| 夫妻午夜视频| tube8黄色片| 制服诱惑二区| 亚洲精品第二区| 最近中文字幕2019免费版| 丰满乱子伦码专区| 亚洲成人免费av在线播放| 大香蕉久久成人网| 一边摸一边抽搐一进一出视频| 久久久久精品人妻al黑| 国产高清国产精品国产三级| 欧美精品高潮呻吟av久久| 91aial.com中文字幕在线观看| 午夜福利影视在线免费观看| 91精品三级在线观看| 少妇的丰满在线观看| 国产无遮挡羞羞视频在线观看| 亚洲欧美一区二区三区黑人| av卡一久久| 国产一卡二卡三卡精品 | 天天躁夜夜躁狠狠躁躁| 日韩一本色道免费dvd| 精品卡一卡二卡四卡免费| 亚洲成人国产一区在线观看 | 不卡av一区二区三区| 欧美日本中文国产一区发布| 国产亚洲最大av| 丁香六月天网| 国产精品成人在线| 在线免费观看不下载黄p国产| 午夜福利在线免费观看网站| 日本爱情动作片www.在线观看| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 丝瓜视频免费看黄片| 欧美日韩成人在线一区二区| 国产精品一区二区精品视频观看| 久久国产精品大桥未久av| 制服人妻中文乱码| 国产av一区二区精品久久| 精品一区二区免费观看| 中文字幕亚洲精品专区| 国产精品一国产av| 侵犯人妻中文字幕一二三四区| 精品国产乱码久久久久久男人| 看免费av毛片| 人人妻人人爽人人添夜夜欢视频| 欧美精品av麻豆av| 国产精品香港三级国产av潘金莲 | 精品久久久精品久久久| 色94色欧美一区二区| 黑人欧美特级aaaaaa片| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 精品亚洲乱码少妇综合久久| 在线观看免费午夜福利视频| 啦啦啦啦在线视频资源| 又大又爽又粗| a级毛片黄视频| 黑丝袜美女国产一区| 精品少妇久久久久久888优播| 精品人妻在线不人妻| 亚洲伊人久久精品综合| 999久久久国产精品视频| 亚洲一区二区三区欧美精品| 777久久人妻少妇嫩草av网站| 国产av码专区亚洲av| xxxhd国产人妻xxx| 中文字幕人妻丝袜一区二区 | 国产欧美亚洲国产| 欧美久久黑人一区二区| 亚洲人成网站在线观看播放| 国产成人啪精品午夜网站| 亚洲国产精品国产精品| 日韩不卡一区二区三区视频在线| 亚洲欧美激情在线| 精品福利永久在线观看| www日本在线高清视频| videosex国产| 久久韩国三级中文字幕| 女性生殖器流出的白浆| 久久国产亚洲av麻豆专区| 深夜精品福利| 久久久精品94久久精品| 午夜福利一区二区在线看| 亚洲一级一片aⅴ在线观看| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 亚洲精品国产区一区二| 午夜av观看不卡| 最近2019中文字幕mv第一页| 久久久亚洲精品成人影院| 热99国产精品久久久久久7| 成人亚洲精品一区在线观看| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜一区二区 | 国产免费现黄频在线看| 美女福利国产在线| 国产一区二区激情短视频 | 亚洲欧洲日产国产| 日韩一区二区视频免费看| 各种免费的搞黄视频| 乱人伦中国视频| 精品国产乱码久久久久久小说| 在线观看www视频免费| 男女边摸边吃奶| 老司机影院成人| 你懂的网址亚洲精品在线观看| 大片免费播放器 马上看| 啦啦啦 在线观看视频| 9色porny在线观看| 久久av网站| 成人影院久久| 欧美黄色片欧美黄色片| 中国国产av一级| 久久毛片免费看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产乱人偷精品视频| 久久久久精品性色| 黄频高清免费视频| 99久国产av精品国产电影| 免费高清在线观看日韩| 精品久久久精品久久久| 狠狠婷婷综合久久久久久88av| 日日爽夜夜爽网站| 亚洲av日韩精品久久久久久密 | 在现免费观看毛片| 人成视频在线观看免费观看| 男人爽女人下面视频在线观看| 你懂的网址亚洲精品在线观看| netflix在线观看网站| 90打野战视频偷拍视频| 美女福利国产在线| 色网站视频免费| 亚洲国产欧美网| 黄色一级大片看看| 国产成人免费无遮挡视频| 看免费成人av毛片| 午夜福利视频精品| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 伦理电影大哥的女人| 日日摸夜夜添夜夜爱| 黄色视频不卡| 亚洲成人手机| 国产熟女欧美一区二区| 日本欧美国产在线视频| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 天天躁夜夜躁狠狠久久av| 亚洲欧美中文字幕日韩二区| 久久久亚洲精品成人影院| 最近2019中文字幕mv第一页| 欧美日韩成人在线一区二区| 999精品在线视频| 一区二区三区精品91| 久久人人爽人人片av| 日韩av在线免费看完整版不卡| 丝袜脚勾引网站| 久久久久久久久久久久大奶| 国产亚洲av高清不卡| 亚洲激情五月婷婷啪啪| 欧美成人午夜精品| 国产 一区精品| 欧美 日韩 精品 国产| 别揉我奶头~嗯~啊~动态视频 | 一边摸一边做爽爽视频免费| 街头女战士在线观看网站| 日韩中文字幕欧美一区二区 | 777米奇影视久久| 飞空精品影院首页| 日韩不卡一区二区三区视频在线| 天天躁夜夜躁狠狠躁躁| 久久精品亚洲av国产电影网| 老司机在亚洲福利影院| 亚洲av中文av极速乱| 久久人人爽人人片av| 在线观看免费日韩欧美大片| 老汉色∧v一级毛片| 99香蕉大伊视频| 久久久久久久国产电影| 精品久久久久久电影网| 波多野结衣av一区二区av| 国产精品亚洲av一区麻豆 | 精品第一国产精品| 99re6热这里在线精品视频| 午夜免费鲁丝| 成人亚洲精品一区在线观看| 高清欧美精品videossex| 亚洲第一区二区三区不卡| 成年人午夜在线观看视频| 欧美在线一区亚洲| 侵犯人妻中文字幕一二三四区| 国产精品国产三级专区第一集| 啦啦啦啦在线视频资源| 日韩精品有码人妻一区| 国产亚洲精品第一综合不卡| 亚洲精品国产一区二区精华液| 考比视频在线观看| 色视频在线一区二区三区| 岛国毛片在线播放| 国产在视频线精品| 久久久久久久久久久久大奶| 在线观看www视频免费| 99精品久久久久人妻精品| 亚洲美女搞黄在线观看| 日韩大片免费观看网站| 秋霞伦理黄片| 91老司机精品| 欧美最新免费一区二区三区| 老司机影院成人| 一级,二级,三级黄色视频| 日韩中文字幕欧美一区二区 | 啦啦啦视频在线资源免费观看| 成人免费观看视频高清| 欧美 亚洲 国产 日韩一| 亚洲av欧美aⅴ国产| 最近的中文字幕免费完整| 国产黄频视频在线观看| 午夜精品国产一区二区电影| 一本—道久久a久久精品蜜桃钙片| 天天影视国产精品| 丰满乱子伦码专区| 久久鲁丝午夜福利片| 91精品伊人久久大香线蕉| 精品国产乱码久久久久久男人| 国产老妇伦熟女老妇高清| 男女边摸边吃奶| 免费看av在线观看网站| 亚洲伊人久久精品综合| 男女无遮挡免费网站观看| 夫妻午夜视频| 捣出白浆h1v1| 久久精品国产综合久久久| 51午夜福利影视在线观看| 老司机在亚洲福利影院| 亚洲色图 男人天堂 中文字幕| 999久久久国产精品视频| 中国三级夫妇交换| 久久毛片免费看一区二区三区| 国产成人av激情在线播放| 色94色欧美一区二区| 欧美精品高潮呻吟av久久| 看非洲黑人一级黄片| videos熟女内射| 国产成人啪精品午夜网站| 久久av网站| 热re99久久国产66热| 你懂的网址亚洲精品在线观看| 97人妻天天添夜夜摸| 国产一区二区三区av在线| 国产深夜福利视频在线观看| 日韩av免费高清视频| 亚洲欧美成人精品一区二区| 中文字幕高清在线视频| 51午夜福利影视在线观看| 亚洲欧美日韩另类电影网站| 久久久久精品久久久久真实原创| 国产又爽黄色视频| 丰满迷人的少妇在线观看| 黄色一级大片看看| 欧美国产精品va在线观看不卡| 最新的欧美精品一区二区| 欧美亚洲日本最大视频资源| 亚洲av男天堂| 国产又色又爽无遮挡免| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av高清一级| 国产乱来视频区| 少妇人妻 视频| 精品一区二区三卡| 日韩成人av中文字幕在线观看| 免费高清在线观看日韩| 国产黄色视频一区二区在线观看| 少妇被粗大猛烈的视频| 国产日韩欧美在线精品| av天堂久久9| 亚洲中文av在线| 日本欧美国产在线视频| 91aial.com中文字幕在线观看| 精品亚洲乱码少妇综合久久| 免费av中文字幕在线| www.av在线官网国产| 99久国产av精品国产电影| 亚洲精品,欧美精品| 人人妻,人人澡人人爽秒播 | 久久97久久精品| 黄频高清免费视频| 久久精品久久久久久久性| 99热网站在线观看| 国产精品国产av在线观看| 久久韩国三级中文字幕| 美女大奶头黄色视频| 欧美亚洲日本最大视频资源| 国产一区亚洲一区在线观看| 亚洲,欧美精品.| 亚洲欧美一区二区三区久久| 又大又爽又粗| 伊人久久国产一区二区| 久久ye,这里只有精品| 久久女婷五月综合色啪小说| 一二三四在线观看免费中文在| 免费在线观看黄色视频的| 欧美97在线视频| 亚洲精品久久午夜乱码| 男女免费视频国产| 欧美成人精品欧美一级黄| 亚洲成国产人片在线观看| 欧美日韩精品网址| 久久久久精品久久久久真实原创| 午夜福利乱码中文字幕| 美女福利国产在线| 久久免费观看电影| 亚洲av成人不卡在线观看播放网 | 中文字幕色久视频| 深夜精品福利| 99久国产av精品国产电影| 日本av免费视频播放| 欧美精品高潮呻吟av久久| 黄片播放在线免费| 欧美激情高清一区二区三区 | 日韩一卡2卡3卡4卡2021年| 久久精品aⅴ一区二区三区四区| 久久99热这里只频精品6学生| 中文字幕人妻丝袜制服| 丰满饥渴人妻一区二区三| 一本色道久久久久久精品综合| 久久av网站| 又黄又粗又硬又大视频| 欧美日韩视频精品一区| 男人添女人高潮全过程视频| 高清不卡的av网站| 天天躁夜夜躁狠狠久久av| 黑人巨大精品欧美一区二区蜜桃| 男女国产视频网站| 国产av精品麻豆| 丰满饥渴人妻一区二区三| 国产成人欧美在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 美女午夜性视频免费| 国产99久久九九免费精品| 亚洲精品aⅴ在线观看| 女的被弄到高潮叫床怎么办| 一边亲一边摸免费视频| 国产一区二区激情短视频 | 无限看片的www在线观看| 青草久久国产| 国产一区有黄有色的免费视频| 黄色怎么调成土黄色| 高清黄色对白视频在线免费看| 在线亚洲精品国产二区图片欧美| 亚洲成国产人片在线观看| 精品少妇黑人巨大在线播放| av女优亚洲男人天堂| 亚洲欧洲国产日韩| 国产一卡二卡三卡精品 | 色94色欧美一区二区| 久久人人爽人人片av| 成人国语在线视频| 亚洲精品美女久久久久99蜜臀 | 精品人妻在线不人妻| 99久久99久久久精品蜜桃| 又大又爽又粗| 天天躁狠狠躁夜夜躁狠狠躁| 啦啦啦中文免费视频观看日本| 久久久久久久久久久久大奶| 伊人久久国产一区二区| 黄片播放在线免费| 亚洲一码二码三码区别大吗| 三上悠亚av全集在线观看| 欧美另类一区| 香蕉丝袜av| 欧美 日韩 精品 国产| 免费av中文字幕在线| 国产有黄有色有爽视频| 亚洲国产欧美日韩在线播放| svipshipincom国产片| 国产av码专区亚洲av| 丝袜在线中文字幕| 一级毛片电影观看| 99精国产麻豆久久婷婷| 老司机靠b影院| 免费av中文字幕在线| 精品国产乱码久久久久久男人| 飞空精品影院首页| 欧美另类一区| av.在线天堂| 狠狠精品人妻久久久久久综合| 中文字幕最新亚洲高清| 欧美激情高清一区二区三区 | 又粗又硬又长又爽又黄的视频| videosex国产| 亚洲精品在线美女| 午夜久久久在线观看| 国产av码专区亚洲av| 久久久久久久大尺度免费视频| 精品一区在线观看国产| 国产伦人伦偷精品视频| 久久久久久久国产电影| 在线观看国产h片| 制服丝袜香蕉在线|