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

    甘肅祁連山國(guó)家自然保護(hù)區(qū)植物群落分布格局及其與環(huán)境因子的關(guān)系

    2022-04-08 01:08:39蘇軍德趙曉冏李國(guó)霞
    廣西植物 2022年3期
    關(guān)鍵詞:植物群落海拔

    蘇軍德 趙曉冏 李國(guó)霞

    摘 要:植物群落分布格局是環(huán)境因子和人類活動(dòng)共同作用的結(jié)果,尤其是海拔梯度被認(rèn)為是植物群落分布格局的決定性因子,為探求甘肅祁連山國(guó)家自然保護(hù)區(qū)植物群落分布格局與環(huán)境因子的關(guān)系及其驅(qū)動(dòng)機(jī)制,該文在野外調(diào)查的基礎(chǔ)上,運(yùn)用數(shù)量分類和除趨勢(shì)典范對(duì)應(yīng)分析(DCCA)排序等方法,探討了研究區(qū)內(nèi)植物群落特征及其與環(huán)境因子的關(guān)系。結(jié)果表明:(1)88個(gè)樣方共記錄物種85種,隸屬30科56屬,利用雙向指示種分析法(TWINSPAN)將其分為9個(gè)植物群落。(2)9個(gè)植物群落在DCCA排序圖上聚集分布,呈現(xiàn)出較好的環(huán)境梯度,其中海拔對(duì)植物群落分布格局影響最大,其次為降水、溫度、坡度、坡向和土壤腐殖質(zhì)。(3)影響植物群落空間分布格局的變量中,環(huán)境因子占25.24%,空間因子占13.21%,空間因子和環(huán)境因子交叉作用占9.03%,群落分布格局未被空間因子和環(huán)境因子解釋部分占52.52%,這部分主要反映了人類活動(dòng)對(duì)研究區(qū)植物群落分布格局的影響。該研究成果對(duì)區(qū)域內(nèi)植被的生態(tài)恢復(fù)和生物多樣性的穩(wěn)定維持具有重要意義。

    關(guān)鍵詞: 植物群落, 祁連山國(guó)家自然保護(hù)區(qū), 分布格局, DCCA, 海拔

    中圖分類號(hào):Q948.15

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-3142(2022)03-0450-10

    Distribution pattern of plant community and its relationship

    with environmental factors in Qilian Mountains

    National Nature Reserve of Gansu Province

    SU Junde1*, ZHAO Xiaojiong2,3, LI Guoxia1

    ( 1. Gansu Vocational & Technical College of Nonferrous Metallurgy, Jinchang 737100, Gansu, China; 2. Gansu

    Academy of Ecological Environment Design and Research, Lanzhou 730000, China; 3. Nourthwest Institute of Eco-Environment

    and Resources, Chinese Academy of Sciences, Lanzhou 733000, China )

    Abstract:The distribution pattern of plant community is the result of the interaction of environmental factors and human activities, especially the elevation gradient is considered as the decisive factor of the distribution pattern of plant community. In order to explore the relationship between plant community distribution pattern and environmental factors and its driving mechanism in Qilian Mountains National Nature Reserve of Gansu Province,? the characteristics of plant community and its relationship with environmental factors in the study area were disucssed, using quantitative classification and DCCA sequencing methods on the basis of field investigation. The results were as follows:(1) There were 85 species of seed plants belonging to 30 families and 56 genera in 88 sampling plots of Qilian Mountains National Nature Reserve of Gansu Province, including twelve species of Gramineae, nine species of Compositae, eight species of Fabaceae, six species of Rosaceae and five species of Chenopodiaceae, accounting for 14.11%, 10.59%, 9.41%, 7.06% and 5.88% of total species, respectively. These plants could be classified into nine plant communities throught TWINSPAN. (2) Nine plant communities were clustered and distributed on the DCCA sequencing diagram, showing a good environmental gradient. Altitude had a great influence on the distribution pattern of plant communities, followed by precipitation, temperature, slope, slope direction and soil humus. (3) Among the variables influencing plant community spatial distribution pattern, environmental factors accounted for 25.24%, spatial factors accounted for 13.21%, space factor and environment factor interactions accounted for 9.03%, 52.52% of the community distribution pattern was not accounted by spatial factor and environmental factor section, and this part mainly reflected the human activity influence on plant community distribution pattern in the study area. The results of this paper are of great significance to the ecological restoration of vegetation and the stable maintenance of biodiversity in the region.

    Key words: plant community, Qilian Mountains National Nature Reserve of Gansu Province, distribution pattern, DCCA, altitude

    植被空間分布作為衡量森林生態(tài)系統(tǒng)結(jié)構(gòu)和功能的重要指標(biāo)(唐志紅等,2020),受海拔、溫度、降水等眾多環(huán)境因子的影響(王金蘭等,2019),能夠很好地反映物種組成和群落的結(jié)構(gòu)特征(張建貴等,2019)。但近年來(lái),隨著全球氣候變暖和人類活動(dòng)的加劇,植被群落面臨了諸如生長(zhǎng)不良、功能衰退、穩(wěn)定性變差等問(wèn)題,嚴(yán)重影響了森林生態(tài)系統(tǒng)的生態(tài)效益及區(qū)域經(jīng)濟(jì)的協(xié)調(diào)發(fā)展,同時(shí)也受到了全社會(huì)和科學(xué)界的高度關(guān)注。李金等(2020)運(yùn)用植物多樣性指數(shù)及冗余分析,深入分析了阿勒泰小克蘭河上游植物的垂直分布格局;單元琪等(2020)通過(guò)實(shí)地植物樣方調(diào)查,發(fā)現(xiàn)海拔是影響三江平原七星河流域濕地植物多樣性的重要因素;劉夢(mèng)婷等(2019)研究表明,森林采伐對(duì)新疆天山云杉林群落分布格局具有十分重要的影響。由此可見(jiàn),研究植物群落的分布格局及其與環(huán)境因子的關(guān)系,對(duì)深入了解物種的演化趨勢(shì)、生態(tài)系統(tǒng)的功能發(fā)揮以及物種分布的驅(qū)動(dòng)機(jī)制都有十分重要的意義。

    甘肅祁連山國(guó)家自然保護(hù)區(qū)地理位置特殊,具有獨(dú)特而典型的自然生態(tài)系統(tǒng)和生物區(qū)系,因其物種豐富(李艷龍等,2020),成為我國(guó)生物多樣性保護(hù)的優(yōu)先區(qū)域,也是西北地區(qū)重要的生物物種資源庫(kù)和野生動(dòng)物遷徙的重要廊道(丁文廣,2019),在維護(hù)我國(guó)西北生態(tài)安全和促進(jìn)區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展方面具有十分重要的作用。近年來(lái),眾多學(xué)者對(duì)區(qū)域內(nèi)生物多樣性(溫婭麗和張虎,2002)、植被NDVI(蘇軍德和李國(guó)霞,2019)、森林生態(tài)效益(汪有奎等,2013)等方面進(jìn)行了研究,但都基于遙感數(shù)據(jù),關(guān)于甘肅祁連山自然保護(hù)區(qū)植物群落分布及其與環(huán)境因子關(guān)系的研究還鮮有報(bào)道。本研究以甘肅祁連山國(guó)家自然保護(hù)區(qū)北坡上的植物群落為研究對(duì)象,通過(guò)實(shí)地樣方調(diào)查并結(jié)合雙向指示種分析法,深入分析區(qū)域內(nèi)植被群落空間分布格局與其環(huán)境因子的關(guān)系,并探索區(qū)域內(nèi)植被群落空間分布的驅(qū)動(dòng)機(jī)制,以期為區(qū)域內(nèi)植被的生態(tài)恢復(fù)和生物多樣性的穩(wěn)定維持奠定一定的理論基礎(chǔ)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于甘肅省祁連山國(guó)家自然保護(hù)區(qū)境內(nèi),地處99°30′ E、38°15′ N(圖1),由于遠(yuǎn)離海洋,長(zhǎng)期受西風(fēng)影響,具有大陸性高寒半濕潤(rùn)山地氣候特征(陳童堯等,2020)。年均氣溫為6 ℃,年均降水量約為550 mm,集中分布在5—9月,年蒸發(fā)量為1 234.2 mm,無(wú)霜期為110 d。區(qū)域植被主要分布在1 500~3 700 m之間,具有明顯的垂直梯帶性和陰陽(yáng)坡差異,主要有青海云杉林、灌木林和少量的祁連圓柏、樺木和山楊林。草原有草甸草原、荒漠草原和高寒草原,區(qū)域內(nèi)植被生長(zhǎng)良好,森林覆蓋率達(dá)28.8%(宋偉宏等,2019)。土壤以山地森林灰褐土、草原栗鈣土、亞高山灌叢草甸土為主。

    1.2 樣方的布設(shè)與調(diào)查

    以垂直高度100 m為間隔,于2017年8月在研究區(qū)內(nèi)不同生境和不同植被群落處設(shè)置88個(gè)30 m×30 m的樣地,每塊樣地中隨機(jī)布設(shè)1個(gè)15 m×15 m的喬木大樣方,在喬木樣方內(nèi)再沿對(duì)角線布設(shè)3個(gè)5 m×5 m的灌木樣方和隨機(jī)布設(shè)1 m×1 m的草本樣方,統(tǒng)計(jì)樣方內(nèi)喬木和灌木的株高、胸徑、株數(shù)以及草本植物的多度和蓋度,同時(shí)詳細(xì)記錄每個(gè)樣方的海拔、坡度、坡向、土壤腐殖質(zhì)厚度等環(huán)境因子,氣溫和降水?dāng)?shù)據(jù)采用寺大隆氣象站(海拔2 600 m)觀測(cè)值。共調(diào)查喬木樣方88個(gè),灌木樣方264個(gè),草本樣方337個(gè)。

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

    以重要值(IV=相對(duì)頻度+相對(duì)多度+相對(duì)顯著度或相對(duì)蓋度)(張世雄等,2020)為樣方物種信息,結(jié)合海拔、氣溫、降水、坡度等環(huán)境因子,運(yùn)用軟件WINTWIN 2.3完成雙向指示種分析(two-way indicator species analysis, TWINSPAN)(王國(guó)宏和楊利民,2001),將88個(gè)樣方劃分成不同植被群落。并在此基礎(chǔ)上,利用軟件CANOCO 5.0完成樣方物種信息的除趨勢(shì)典范對(duì)應(yīng)分析(detrended canonical correspondence analysis,DCCA)(潘婷婷等,2020)和環(huán)境-空間變量分離分析。

    2 結(jié)果與分析

    2.1 典型樣地植物群落物種組成

    本文共調(diào)查樣地88塊,記錄種子植物85種,隸屬30科56屬,其中禾本科12種、菊科9種、豆科8種、薔薇科6種、藜科5種,分別占物種總數(shù)的14.11%、10.59%、9.41%、7.06%和5.88%。植被群落由喬木層、灌木層和草本層組成,其中喬木層5種,灌木層13種,草本層67種。群落不同層片物種重要值(表1)表明,在甘肅祁連山國(guó)家自然保護(hù)區(qū)喬木層優(yōu)勢(shì)種有祁連圓柏、青海云杉、紅樺、山楊和白樺,灌木層優(yōu)勢(shì)種有鬼箭錦雞兒、狹葉錦雞兒、金露梅、白刺、紅砂、鹽爪爪、合頭草和蒙古蕕,草本層優(yōu)勢(shì)種有小甘菊、冰草、芨芨草、牛筋草、甘肅萱草、甘青蒿、垂穗披堿草、小花鳳毛菊、早熟禾和針茅。

    2.2 植被群落分類

    TWINSPAN是將樣方與物種進(jìn)行雙向聚類,以二歧式的方式對(duì)群落類型進(jìn)行分割劃分。 依據(jù)《中國(guó)植被》的分類原則,以重要值為依據(jù),對(duì)甘肅祁連山自然保護(hù)區(qū)88個(gè)樣方進(jìn)行TWINSPAN分類,并根據(jù)群落生境的指示種和優(yōu)勢(shì)種進(jìn)行命名,結(jié)果將88個(gè)樣方劃分為9個(gè)群落類型(圖2)。

    群叢Ⅰ:祁連圓柏+鬼箭錦雞兒+甘肅萱草,由喬木層、灌木層和草本層構(gòu)成,喬木層以祁連圓柏為主,郁閉度為0.27~0.52,灌木層主要有鬼箭錦雞兒和金露梅,草本層有甘肅萱草、甘青蒿等,分布于海拔3 050~3 300 m之間,包括樣方34、43、62、66、67、68、70、71、73、74、75、80、83、84、85、86、87,伴有青海云杉、金露梅、小葉金露梅、早熟禾、裂葉獨(dú)活、針茅、蒲公英、蘚生馬先蒿、小米草、肉果草、龍膽等,草本植物蓋度為82%。

    群叢Ⅱ:鬼箭錦雞兒+山生柳+青藏蓼,由灌木層和草本層構(gòu)成,灌木層主要以鬼箭錦雞兒和山生柳為主,草本層主要有青藏蓼、珠芽蓼等,分布于海拔3 190~3 700 m之間,包括樣方61、63、64、65、69、72、76,伴有青海杜鵑、甘肅萱草、高原毛茛、垂穗披堿草、早熟禾、火絨草等,草本植物蓋度為78%。

    群叢Ⅲ:銀露梅+狹葉錦雞兒+甘肅萱草,由灌木層和草本層構(gòu)成,灌木層主要以銀露梅和狹葉錦雞兒為主,草本層主要有甘肅萱草、垂穗披堿草等,分布于海拔2 350~2 600 m之間,包括樣方25、29、32、37、39、44、55、56、59、77、78、79、88,伴有高頭烏、鳳毛菊、針茅、早熟禾、唐松草、甘肅萱草、裂葉獨(dú)活、甘肅馬先蒿、雙花堇菜,草本植物蓋度為85%。

    群叢Ⅳ:青海云杉+金露梅+垂穗披堿草,由喬木層、灌木層和草本層構(gòu)成,喬木層主要以青海云杉為主,郁閉度為0.43~0.67,灌木層有金露梅等,草本層主要有垂穗披堿草、冰草等,分布于海拔2 450~2 950 m之間,包括樣方30、31、35、40、45、47、48、49、51、53、57、58、60、81、82,伴有蒲公英、車前、針茅、早熟禾、濕地繁縷,草本層植物蓋度為72%。

    群叢Ⅴ:蒙古蕕+狹葉錦雞兒+甘肅馬先蒿,由灌木層和草本層構(gòu)成,灌木層主要有蒙古蕕和狹葉錦雞兒,草本層主要有甘肅馬先蒿,分布于海拔2 360~2 750 m之間,包括樣方26、27、28、33、36、38、41、42、46、50、52、54,伴有金露梅、冰草、紅花薔薇、早熟禾、堇菜、唐松草、鳳毛菊、裂葉獨(dú)活、高頭烏,草本層植物蓋度為67%。

    群叢Ⅵ:狹葉錦雞兒+甘蒙錦雞兒+冰草,由灌木層和草本層構(gòu)成,灌木層主要以狹葉錦雞兒和甘蒙錦雞兒為主,草本層主要有冰草,分布于海拔2 300~2 360 m之間,包括樣方5、6、14、21、22、23、24,伴有甘肅馬先蒿、雙花堇菜、唐松草、裂葉獨(dú)活、拉拉藤,草本植物蓋度為58%。

    群叢Ⅶ:紅砂+合頭草+芨芨草,由灌木層和草本層構(gòu)成,灌木層主要有紅砂和合頭草,草本層以芨芨草為主,分布于海拔1 550~1 730 m之間,包括樣方2、3、9、10、17,伴有旱柳、沙棘、披堿草、冰草、早熟禾、狼毒花,草本植物蓋度為64%。

    群叢Ⅷ:鹽爪爪+紅砂+牛筋草,由灌木層和草本層構(gòu)成,灌木層主要有鹽爪爪和紅砂,草本層以牛筋草為主,分布于海拔1 500~1 800 m之間,包括樣方11、12、13、15、16、20,伴有沙棘、芨芨草、白蓮蒿、冰草、狼毒花、早熟禾,草本植物蓋度為70%。

    群叢Ⅸ:白刺+檉柳+冰草,由灌木層和草本層構(gòu)成,灌木層主要有白刺和檉柳,草本層以冰草為主,分布于海拔1 500~1 700 m之間,包括樣方1、4、7、8、18、19,伴有旱柳、狼毒花、牛筋草、早熟禾、白蓮蒿,草本植物蓋度為68%。

    2.3 植被群落的DCCA排序分析

    2.3.1 環(huán)境因子與DCCA排序軸的相關(guān)性 甘肅祁連山自然保護(hù)區(qū)植被群落DCCA排序結(jié)果(表2)表明前四個(gè)排序軸對(duì)物種和環(huán)境關(guān)系的解釋量累計(jì)為89.75%,其中第一軸為64.21%,這說(shuō)明第一軸在對(duì)群落與環(huán)境關(guān)系的解釋中占主導(dǎo)地位。蒙特卡洛檢驗(yàn)(Du et al., 2020)也同樣表明,前四個(gè)排序軸所反映的環(huán)境因子與物種呈顯著相關(guān)(P<0.01),其中第一軸蒙特卡洛檢驗(yàn)結(jié)果為F=11.875,P=0.001。

    各環(huán)境因子與第一軸呈顯著相關(guān)關(guān)系,相關(guān)性大小表現(xiàn)為海拔>年均降水>年均溫度>坡度>土壤腐殖質(zhì)>坡位>坡向,表明第一軸主要反映了由海拔變化引起的水熱生態(tài)環(huán)境、地形及土壤因子的變化。海拔、年均降水、坡度和年均溫度與第二軸呈顯著相關(guān)關(guān)系,相關(guān)性大小表現(xiàn)為年均溫度>年均降水>坡度>海拔,表明第二軸主要反映了溫度的梯度變化(表2)。

    2.3.2 植被群落樣方與環(huán)境因子的DCCA排序 甘肅祁連山自然保護(hù)區(qū)植被群落DCCA分析結(jié)果表明(表2),前四個(gè)排序軸的特征值分別為0.962、0.643、0.471和0.371,由于第一排序軸和第二排序軸特征值相對(duì)較高,反映了重要的生態(tài)學(xué)意義,因此采用前兩個(gè)排序軸數(shù)據(jù)作群落樣方與環(huán)境因子DCCA二維排序圖(圖3),結(jié)果顯示,88個(gè)樣方被分為9個(gè)群落類型,且在DCCA二維排序圖中聚集分布,反映了海拔-水-熱的梯度變化。群叢Ⅸ到群叢Ⅰ沿第一排序軸從左到右交叉分布,海拔從1 500 m到3 700 m,依次出現(xiàn)了紅砂、鹽爪爪、白刺、狹葉錦雞兒、金露梅、青海云杉、祁連圓柏、鬼箭錦雞兒等生態(tài)習(xí)性各異的植被系列,表明了海拔的梯度變化。群叢Ⅵ到群叢Ⅰ沿第二排序軸從下往上交叉分布,表明了溫度的梯度變化。

    樣方與環(huán)境因子的DCCA二維排序圖直觀地展現(xiàn)了樣方群落與環(huán)境因子的相互關(guān)系(圖3),DCCA第一軸反映海拔的梯度變化,從左往右海拔依次升高,第二軸反映溫度的梯度變化,從下往上

    溫度依次升高。群落Ⅶ(紅砂+合頭草+芨芨草)、群落Ⅷ(鹽爪爪+紅砂+牛筋草)、群落Ⅸ(白刺+檉柳+冰草)屬旱生荒漠群叢,分布在海拔相對(duì)較低、溫度較高的區(qū)域,但在第二軸上并未產(chǎn)生明顯的分化。群落Ⅰ(祁連圓柏+鬼箭錦雞兒+甘肅萱草)受海拔、坡度和降水的影響較大,分布在海拔較高、坡度較陡、降水較多的區(qū)域。群落Ⅲ(銀露梅+狹葉錦雞兒+甘肅萱草)受制于坡向和海拔的影響,主要分布在海拔相對(duì)較低、土壤腐殖質(zhì)較厚的區(qū)域。群落Ⅳ(青海云杉+金露梅+垂穗披堿草)的限制因子是海拔和溫度, 主要分布在海拔相對(duì)較低、生境干旱、土壤腐殖質(zhì)較厚的區(qū)域。群落Ⅴ(蒙古蕕+狹葉錦雞兒+甘肅馬先蒿)主要分布在海拔較高、坡度陡峭、土壤腐殖質(zhì)較薄的區(qū)域。群落Ⅱ(鬼箭錦雞兒+山生柳+青藏蓼)分布在高海拔的冷濕區(qū)域。

    2.3.3 影響群落分布格局環(huán)境因子的定量分析 DCCA分析結(jié)果表明,海拔、降水、溫度、坡度、坡向、土壤腐殖質(zhì)共同影響了研究區(qū)群落物種的分布格局,依據(jù)Borcard等(1992)的方法,對(duì)影響甘肅祁連山自然保護(hù)區(qū)群落分布格局的環(huán)境空間變量定量分離的結(jié)果顯示,環(huán)境因子對(duì)群落分布格局的解釋率占25.24%,空間因子對(duì)群落分布格局解釋率占13.21%,空間因子和環(huán)境因子交叉作用對(duì)群落分布格局的解釋率占9.03%,群落分布格局未被空間因子和環(huán)境因子解釋部分占52.52%。

    3 討論

    物種組成是反映植物群落結(jié)構(gòu)變化的重要因子(Aardt et al.,2020),甘肅祁連山國(guó)家自然保護(hù)區(qū)因其豐富的植物多樣性,在維護(hù)我國(guó)西部生態(tài)安全方面具有重要作用。本研究調(diào)查的88個(gè)樣方中,發(fā)現(xiàn)種子植物85種,隸屬30科56屬。禾本科、菊科和豆科的物種在研究區(qū)內(nèi)占主要優(yōu)勢(shì),這與唐志紅對(duì)祁連山東段植物群落特征研究結(jié)果較為一致。研究區(qū)共發(fā)現(xiàn)喬木5種,主要以青海云杉為主,且有絕對(duì)優(yōu)勢(shì),主要分布在中海拔的半陰坡和半陽(yáng)坡,在高海拔的陽(yáng)坡有少量祁連圓柏分布,在峽谷地帶分布有少量紅樺和白樺,山楊在河谷地段呈斑狀分布。由于喬木層物種在甘肅祁連山國(guó)家自然保護(hù)區(qū)北坡分布較少(張超,2020),加之其個(gè)體生長(zhǎng)緩慢,所以植被組成相對(duì)穩(wěn)定(蔣友嚴(yán)等,2017),而灌木層和草本層植被物種多樣性和豐富度卻隨著喬木層植被的演替而不斷變化(張永虎,2010)。灌木作為甘肅祁連山自然保護(hù)區(qū)北坡的主要植物類群,從低海拔到高海拔都有分布,共有13種,占植物種類的15.29%。

    植物群落的空間分布是海拔及眾多環(huán)境因子共同作用的結(jié)果,但近年來(lái),海拔被認(rèn)為是植物群落空間分布的決定性因素。劉夢(mèng)婷等(2019)認(rèn)為天山云杉群落的空間分布的主導(dǎo)因子是海拔;高輝等(2020)研究表明由海拔梯度引起的濕度和土壤含氮量是決定色季拉山物種多樣性及其空間分布的重要因素;潘婷婷等(2020)也發(fā)現(xiàn)海拔是南京北部郊野森林群落結(jié)構(gòu)和植物生長(zhǎng)的限制因子。本研究TWINSPAN分析劃分的9個(gè)群落類型,隨海拔梯度的升高,依次出現(xiàn)了超旱生荒漠類型(群落Ⅶ、Ⅷ、Ⅸ)、山地旱生類型(群落Ⅴ、Ⅵ)、山地中生類型(群落Ⅲ、Ⅳ)、山地濕生類型(群落Ⅰ)、山地高寒類型(群落Ⅱ),體現(xiàn)了甘肅祁連山國(guó)家自然保護(hù)區(qū)植被群落的垂直梯帶性分布特征。易玉媛和王紅義(2013)的研究也同樣表明,祁連山北坡植被群落隨著海拔的升高,依次呈現(xiàn)出荒漠草原植被、干性灌叢草原植被、山地森林草原植被、亞高山灌叢草甸植被、高山寒漠草甸植被的分布格局。由此可見(jiàn),海拔是影響祁連山自然保護(hù)區(qū)北坡植被群落空間分布最重要的環(huán)境因子,決定著不同地理氣候帶群落尺度上植被的空間分布。

    在區(qū)域及較大尺度上,降水和溫度是影響物種多樣性和植被分布的重要因素(劉冠成等,2018),而在景觀及小尺度上,環(huán)境因子則主要決定植被的空間分布(劉春艷等,2017)。蒙特卡洛檢驗(yàn)結(jié)果顯示,所有排序軸都通過(guò)統(tǒng)計(jì)檢驗(yàn)(F=2.531,P<0.001),這說(shuō)明所選環(huán)境因子均對(duì)研究區(qū)植被群落空間分布具有重要影響。甘肅祁連山國(guó)家自然保護(hù)區(qū)物種多樣性和植被群落分布在區(qū)域尺度上主要受降水和溫度的影響,隨著海拔的升高,降水增多,溫度降低(陳銀萍等,2013),以420 mm降水量為界,研究區(qū)植被形成了不同的演替格局。而在景觀尺度上,地形及土壤腐殖質(zhì)是植被群落空間分布的驅(qū)動(dòng)因子,隨著海拔的升高,土壤鹽分減小,有機(jī)質(zhì)及全氮含量增加(張蕊等,2014),且在不同坡度和坡向土壤肥力又有差異(賈鵬麗等,2020),南坡及坡度較緩的區(qū)域土壤肥力較高,草本及灌木主要分布在西南坡和西北坡,喬木則主要分布在北坡。由此可見(jiàn),由海拔主導(dǎo)的水熱及地形土壤因素是影響研究區(qū)植被群落分布格局的重要環(huán)境因子(Xin et al.,2020)。值得一提的是,研究區(qū)群落分布格局的環(huán)境空間變量定量分離結(jié)果顯示,群落分布格局未被空間因子和環(huán)境因子解釋部分占52.52%,這說(shuō)明除環(huán)境因子外,人類活動(dòng)(如砍伐、放牧)對(duì)區(qū)域內(nèi)植被群落的空間分布也有一定的影響(修麗娜,2014)。因此,加大對(duì)研究區(qū)現(xiàn)有森林資源的保護(hù)和宣傳教育,對(duì)甘肅祁連山國(guó)家自然保護(hù)區(qū)生態(tài)恢復(fù)和重建具有十分重要的意義。

    4結(jié)論

    (1)甘肅祁連山國(guó)家自然保護(hù)區(qū)植被群落包括種子植物85種,隸屬30科56屬,其中禾本科12種、菊科9種、豆科8種、薔薇科6種、藜科5種,分別占物種總數(shù)的14.11%、10.59%、9.41%、7.06%和5.88%。

    (2)甘肅祁連山國(guó)家自然保護(hù)區(qū)植被群落包括9個(gè)群叢類型:祁連圓柏+鬼箭錦雞兒+甘肅萱草、鬼箭錦雞兒+山生柳+青藏蓼、銀露梅+狹葉錦雞兒+甘肅萱草、青海云杉+金露梅+垂穗披堿草、蒙古蕕+狹葉錦雞兒+甘肅馬先蒿、狹葉錦雞兒+甘蒙錦雞兒+冰草、紅砂+合頭草+芨芨草、鹽爪爪+紅砂+牛筋草、白刺+檉柳+冰草。

    (3)海拔是影響甘肅祁連山自然保護(hù)區(qū)植被群落分布格局的主要因素,其次為降水、溫度、坡向、坡度和土壤腐殖質(zhì)。

    (4)影響植物群落空間分布格局的變量中,環(huán)境因子解釋了群落分布格局的25.24%,空間因子解釋了群落分布格局的13.21%,空間因子和環(huán)境因子交叉作用對(duì)群落分布格局的解釋率占9.03%,群落分布格局未被空間因子和環(huán)境因子解釋部分占52.52%。

    參考文獻(xiàn):

    AARDT ACV, CODRON D, THERON EJ, et al.,2020. Plant community structure and possible vegetation changes after drought on a granite catena in the Kruger National Park[J]. S Africa, 62(2): e1-e11.

    BORCARD D, LEGENDRE P, DRAPEAU P, 1992. Partialling out the spatial component of ecological variation[J]. Ecology, 15(73): 1045-1055.

    CHEN TY, JIA YF, WANG JN, et al., 2020. Current situation and function of soil conservation in National Nature Reserves in the Qilian mountains based on InVEST model[J]. Arid Zone Res, 37(1):150-159. [陳童堯, 賈燕鋒, 王佳楠, 等, 2020. 基于InVEST模型的祁連山國(guó)家級(jí)自然保護(hù)區(qū)土壤保持現(xiàn)狀與功能[J]. 干旱區(qū)研究, 37(1):150-159.]

    CHEN YP, YU F, YANG ZJ, et al., 2013. Size structure and dynamic of the Picea crassifolia populations on eastern Qilian mountains [J]. Guhaia, 33(2): 236-241. [陳銀萍, 于飛, 楊宗娟, 等, 2013. 東祁連山青海云杉種群大小結(jié)構(gòu)及其動(dòng)態(tài)研究[J]. 廣西植物, 33(2): 236-241.]

    DU BJ, LIU CY, YANG JX, et al., 2020. Abutment pressure distribution pattern and size optimization of coal pillar under repeated mining: A case study[J]. Arab J Geosci, 13(23):1261.

    DING WG, 2019. Innovative research on ecological environment protection system and mechanism in Qilian Mountain Nature Reserve of Gansu Province[D]. Lanzhou: Lanzhou University: 23-25. [丁文廣, 2019. 甘肅祁連山自然保護(hù)區(qū)生態(tài)環(huán)境保護(hù)體制機(jī)制創(chuàng)新研究[D]. 蘭州: 蘭州大學(xué): 23-25.]

    GAO H, LIU LJ, FANG JP, 2020.Distribution pattern and environmental interpretation of forest communities along elevation gradients in Sejila Mountain \[J\]. J Nat Sci Hunan Norm Univ, 43(5): 39-47. \[高輝,劉麗娟,方江平,2020. 色季拉山沿海拔梯度的森林群落分布格局及環(huán)境解釋\[J\]. 湖南師范大學(xué)自然科學(xué)學(xué)報(bào),43(5):39-47.\]

    JIA PL, FENG HY, LI M, et al., 2020. Soil microbial diversity of black soil under different land use patterns in northeast China[J]. Trans Chin Soc Agric Eng, 36(20): 171-178. [賈鵬麗, 馮海艷, 李淼, 等, 2020. 東北黑土區(qū)不同土地利用方式下農(nóng)田土壤微生物多樣性[J]. 農(nóng)業(yè)工程學(xué)報(bào), 36(20): 171-178.]

    JIANG YY, DU WT, HUANG J,et al., 2017. Analysis of vegetation changes in the Qilian mountains during 2000-2015 [J]. J Glac Geo, 39(5): 1130-1136. [蔣友嚴(yán), 杜文濤, 黃進(jìn), 等, 2017. 2000-2015年祁連山植被變化分析[J]. 冰川凍土, 39(5): 1130-1136.]

    LIU CY, ZHU KW, LIU JP, 2017. Evolution and predication of land cover and biodiversity function in Chongqing section of Three Gorges Reservoir area[J]. Trans Chin Soc Agric Eng, 33(19): 258-267. [劉春艷, 朱康文, 劉吉平, 2017. 三峽庫(kù)區(qū)重慶段土地覆蓋和生物多樣性功能演化及預(yù)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 33(19): 258-267.]

    LIU GC, HUANG YX, WANG QG,et al., 2018. Effects of environment factors on plant species diversity research progress[J]. Chin Agric Sci Bull, 34(13): 83-89. [劉冠成, 黃雅曦, 王慶貴, 等, 2018. 環(huán)境因子對(duì)植物物種多樣性的影響研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào), 34(13): 83-89.]

    LIU MT, WANG ZX, WANG YP,et al., 2019. Plant communities pattern of Picea tianschanica forest and their interrelations with environmental factors in Tianshan area[J]. For Res, 32(6): 90-98. [劉夢(mèng)婷, 王振錫, 王雅佩, 等, 2019. 新疆天山云杉林群落分布格局及環(huán)境解釋[J]. 林業(yè)科學(xué)研究, 32(6): 90-98.]

    LI J, XU HL, ZHAO XF, et al., 2020. Vertical distribution and environment analysis of plant communities in upper reaches of Crane River in Aletai[J]. Adm Technol Environ Monit, 32(4):10-14. [ 李金, 徐海量, 趙新風(fēng), 等, 2020. 阿勒泰小克蘭河上游植物群落的垂直分布格局及環(huán)境分析[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù), 32(4):10-14.]

    LI YL, WANG ZP, LIU HT, et al., 2020. Impact assessment of road engineering on biodiversity in Qilian Mountain National Nature Reserve of Gansu Province[J]. Prot For Sci Technol, 2(8): 44-47. [李艷龍, 王作鵬, 劉海濤, 等, 2020. 道路工程對(duì)甘肅祁連山國(guó)家級(jí)自然保護(hù)區(qū)生物多樣性影響評(píng)價(jià)[J]. 防護(hù)林科技, 2(8): 44-47.]

    PAN TT,CHEN L, YANG GD,et al., 2020. Communities species diversity and environment interpretation of Suberban Forest in Northern Nanjing[J]. J Nanjing For Univ (Nat Sci Ed), (10): 48-54. [潘婷婷, 陳林, 楊國(guó)棟, 等, 2020. 南京北部郊野森林群落物種多樣性及其環(huán)境解釋[J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), (10): 48-54.]

    SHAN YQ, YAO YL, ZHANG XX,et al.,2020. Study on wetland plant diversity and its influencing factors in the Qixiang River Basin of Sanjiang Plain[J]. Acta Ecol Sin, 40(5): 1629-1636. [單元琪, 姚允龍, 張欣欣, 等, 2020. 三江平原七星河流域濕地植物多樣性及影響因素[J]. 生態(tài)學(xué)報(bào), 40(5):1629-1636.]

    SU JD, LI GX, 2019. Spatial and temporal variation characteristics of vegetation NDVI in Qilian mountains based on GIS. Mine Surv[J], 47(4): 33-38. [蘇軍德, 李國(guó)霞, 2019. 基于GIS的祁連山區(qū)植被NDVI時(shí)空變化特征研究[J]. 礦山測(cè)量, 47(4): 33-38.]

    SONG WH, WANG LN, ZHANG JL, 2019. Temporal and spatial changes in response to climate in grasslands of the Qilian Mountain National Nature Reserve[J]. Pratac Sci, 36(9): 2233-2249. [宋偉宏, 王莉娜, 張金龍, 2019. 甘肅祁連山自然保護(hù)區(qū)草地時(shí)空變化及其對(duì)氣候的響應(yīng)[J]. 草業(yè)科學(xué), 36(9): 2233-2249.]

    TANG ZH, YU QS, LIU HJ,et al., 2020. Characteristics of alpine vegetation communityand its relationship to topographic climate factors in the eastern Qilian mountain[J]. Acta Ecol Sin, 40(1): 223-232.[唐志紅, 尉秋實(shí), 劉虎俊, 等, 2020. 祁連山東段高寒植被群落特征及其與地形氣候因子關(guān)系研究[J]. 生態(tài)學(xué)報(bào), 40(1): 223-232.]

    WEN YL, ZHANG H, 2002. Plant diversity and distribution in the north slope of Qilian mountains[J]. J Gansu For Sci Technol, (1): 37-39. [溫婭麗, 張虎, 2002. 祁連山(北坡)植物多樣性及其分布特點(diǎn)[J]. 甘肅林業(yè)科技, (1): 37-39.]

    WANG JL, CAO WX, ZHANG DG, et al., 2019. Structure and species diversity of alpine Rhododendron shrub-herb community and its response to altitude gradients in eastern Qilian mountains[J]. Grassl Turf, 39(5): 1-9. [王金蘭, 曹文俠, 張德罡, 等, 2019. 東祁連山高寒杜鵑灌叢群落結(jié)構(gòu)和物種多樣性對(duì)海拔梯度的響應(yīng)[J]. 草原與草坪, 39(5): 1-9.]

    WANG GH, YANG LM, 2001. Gradient analysis environment interpretation of woody plant communities the middle section of the northern slopes of Qilian mountain, Gansu, China[J]. Acta Phytoecol Sin, 1(6): 733-740. [王國(guó)宏, 楊利民, 2001. 祁連山北坡中段森林植被梯度分析及環(huán)境解釋[J]. 植物生態(tài)學(xué)報(bào),? 1(6): 733-740.]

    WANG YK, GUO SX, WANG J,et al., 2013. Estimation of forest ecosystem service value in the Qilian mountains National Nature Reserve in Gansu of China[J]. J Desert Res, 33(6): 1905-1911. [汪有奎, 郭生祥, 汪杰, 等, 2013. 甘肅祁連山國(guó)家級(jí)自然保護(hù)區(qū)森林生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估[J]. 中國(guó)沙漠, 33(6): 1905-1911.]

    XIN G, PENG Y, MA XH, et al., 2020. Design of distillation reactor with novel catalysts distribution pattern for n-amyl acetate synthesis in industrial scale[J]. Fuel, 280: 1-8.

    XIU LN, 2014. Spatio-temporal distribution characteristics

    of potential natural vegetation in China based on CSCS Model[D]. Lanzhou: Lanzhou University: 3-18. [修麗娜, 2014. 基于CSCS模型的中國(guó)潛在自然植被時(shí)空分布特征研究[D]. 蘭州: 蘭州大學(xué): 3-18.]

    YI YY, WANG HY. 2013. Vertical distribution of vegetation communities on the north slope of Qilian mountains[J]. Prot For Sci Technol, 10: 6-9. [易玉媛, 王紅義, 2013. 祁連山北坡植被群落垂直分布特征[J]. 防護(hù)林科技, 10: 6-9.]

    ZHANG C, 2020. Investigation andanalysis of biodiversity and natural ecosystem in Qilian Mountain National Nature Reserve[J]. For Sci Technol Inf, 52(3): 12-14. [ 張超, 2020. 祁連山國(guó)家公園生物多樣性和自然生態(tài)系統(tǒng)的特征調(diào)查分析[J]. 林業(yè)科技情報(bào), 52(3): 12-14.]

    ZHANG JG, WANG LD, YAO T, et al., 2019. Plant community structure and species diversity differences in alpine grassland in the Qilian mountains with different levels of degradation[J]. Acta Agr Sin, 28(5): 15-25. [張建貴, 王理德, 姚拓, 等, 2019. 祁連山高寒草地不同退化程度植物群落結(jié)構(gòu)與物種多樣性研究[J]. 草業(yè)學(xué)報(bào), 28(5): 15-25.]

    ZHANG R, CAO JJ, GUO RY,et al., 2014. Effect of converting cultivated land to forest-grassland on soil carbon, nitrogen and phosphoras in sub-alpine grassland region of north slope of Qilian mountains[J]. Ecol Environ Sci,23(6): 938-944. [張蕊, 曹靜娟, 郭瑞英, 等, 2014. 祁連山北坡亞高山草地退耕還林草混合植被對(duì)土壤碳氮磷的影響[J]. 生態(tài)環(huán)境學(xué)報(bào), 23(6): 938-944.]

    ZHANG SX, QIN RM, YANG XY,et al., 2020. Patterns of species diversity in herbaceous communities an different altitude gradients and their relationships with environment factors in the Luliang mountanis in Shanxi Province,China[J]. Guihaia, 40(12): 1860-1868. [張世雄, 秦瑞敏, 楊曉艷, 等, 2020. 山西呂梁山草本群落物種多樣性的海拔梯度格局及與環(huán)境因子的關(guān)系[J]. 廣西植物, 40(12): 1860-1868.]

    ZHANG YH, 2010. Vegetation distribution and succession in the eastern part of Qilian mountains[J]. Yellow River, 32(8): 94-95. [張永虎, 2010. 祁連山東段植被分布特征及演替規(guī)律[J]. 人民黃河, 32(8): 94-95.]

    (責(zé)任編輯 李 莉)

    收稿日期:2021-05-19

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(41661016);甘肅省社會(huì)科學(xué)規(guī)劃項(xiàng)目(19YB155);甘肅省人文社會(huì)科學(xué)項(xiàng)目(20ZZ11);甘肅省青年科技基金計(jì)劃(18JR3RC420);金昌市青年人才基金(2020CZA52W);甘肅省高等學(xué)校創(chuàng)新基金項(xiàng)目(2020A-297);隴原青年創(chuàng)新創(chuàng)業(yè)人才項(xiàng)目(2019-12) [Supported by? National Natural Science Foundation of China (41661016); Social Science Planning Project of Gansu Province(19YB155); General Project of Humanities and Social Science of Gansu Province(20ZZ11); Youth Science and Technology Fund of Gansu Province(18JR3RC420); Jinchang Young Talent Fund(2020CZA52W); Innovation Fund of Colleges and Universities in Gansu Province(2020A-297); Project of Longyuan Youth Innovative and Entrepreneurial Talents (2019-12)]。

    第一作者: 蘇軍德(1986-),碩士,副教授,主要從事生態(tài)環(huán)境學(xué)研究,(E-mail)sujd2009@163.com。

    猜你喜歡
    植物群落海拔
    高海拔的甘南
    南美三大高海拔球場(chǎng)
    海拔4000米的父愛(ài)
    普達(dá)措,高海拔桃源秘境
    上海外環(huán)綠帶植物群落鄰體干擾特征分析
    綠色科技(2016年21期)2016-12-27 17:06:10
    談城市濱水植物群落的構(gòu)建
    縱向嶺谷區(qū)植物群落研究
    華北土石山區(qū)植物群落土壤呼吸特性及其影響因子
    珍稀瀕危植物專類園植物配置模式研究
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    美女国产高潮福利片在线看| 大型av网站在线播放| 免费搜索国产男女视频| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 国产精品一区二区精品视频观看| 亚洲精品一区av在线观看| 久久久水蜜桃国产精品网| 亚洲欧美一区二区三区黑人| 亚洲成a人片在线一区二区| 成人18禁高潮啪啪吃奶动态图| 久久精品亚洲精品国产色婷小说| 黑人巨大精品欧美一区二区蜜桃| 麻豆久久精品国产亚洲av| 两个人免费观看高清视频| 久久久国产精品麻豆| 一夜夜www| 亚洲av美国av| 十八禁人妻一区二区| 欧美成人免费av一区二区三区| 看片在线看免费视频| 在线观看www视频免费| 美女大奶头视频| 琪琪午夜伦伦电影理论片6080| 禁无遮挡网站| 亚洲在线自拍视频| 色哟哟哟哟哟哟| 黄色女人牲交| 一边摸一边抽搐一进一小说| 亚洲天堂国产精品一区在线| 亚洲中文字幕日韩| 精品乱码久久久久久99久播| 大香蕉久久成人网| 淫妇啪啪啪对白视频| 一二三四在线观看免费中文在| 成人三级黄色视频| 狠狠狠狠99中文字幕| 亚洲欧美激情综合另类| 精品久久久精品久久久| 亚洲美女黄片视频| 男女之事视频高清在线观看| 精品午夜福利视频在线观看一区| 国产免费男女视频| 国产真人三级小视频在线观看| 亚洲av日韩精品久久久久久密| 久久午夜亚洲精品久久| 窝窝影院91人妻| 精品一区二区三区视频在线观看免费| 免费高清视频大片| 亚洲人成电影免费在线| 成人精品一区二区免费| 久久人人97超碰香蕉20202| 免费高清视频大片| 非洲黑人性xxxx精品又粗又长| 免费在线观看日本一区| 激情视频va一区二区三区| 欧美日韩亚洲综合一区二区三区_| 在线观看午夜福利视频| 多毛熟女@视频| 免费搜索国产男女视频| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产色婷婷电影| 美女国产高潮福利片在线看| 免费av毛片视频| 色播亚洲综合网| 99久久综合精品五月天人人| 9热在线视频观看99| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 欧美日韩一级在线毛片| 黑人操中国人逼视频| 午夜精品国产一区二区电影| 亚洲狠狠婷婷综合久久图片| 老司机午夜福利在线观看视频| 无限看片的www在线观看| 男女午夜视频在线观看| 日韩av在线大香蕉| 日韩欧美国产在线观看| 色哟哟哟哟哟哟| 天堂√8在线中文| 波多野结衣一区麻豆| 淫妇啪啪啪对白视频| 欧美日韩黄片免| 免费在线观看完整版高清| 久久天躁狠狠躁夜夜2o2o| 成人18禁在线播放| 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| 无人区码免费观看不卡| 99国产精品一区二区三区| 欧美av亚洲av综合av国产av| 电影成人av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区国产一区二区| 1024香蕉在线观看| 亚洲成av人片免费观看| 精品第一国产精品| 欧美亚洲日本最大视频资源| 99在线人妻在线中文字幕| 午夜免费激情av| 色尼玛亚洲综合影院| 国产精品99久久99久久久不卡| 美女大奶头视频| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 国产乱人伦免费视频| 日韩 欧美 亚洲 中文字幕| 满18在线观看网站| 高清在线国产一区| 中文亚洲av片在线观看爽| 少妇被粗大的猛进出69影院| 91国产中文字幕| 制服人妻中文乱码| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区三| 麻豆av在线久日| 午夜福利一区二区在线看| 国产亚洲精品久久久久5区| 精品欧美一区二区三区在线| 侵犯人妻中文字幕一二三四区| 91麻豆精品激情在线观看国产| 视频在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 黄片小视频在线播放| 禁无遮挡网站| 男男h啪啪无遮挡| 此物有八面人人有两片| 久久人妻福利社区极品人妻图片| 成熟少妇高潮喷水视频| 久久香蕉精品热| 亚洲精品久久成人aⅴ小说| 麻豆av在线久日| 伦理电影免费视频| 悠悠久久av| 亚洲中文av在线| 成年人黄色毛片网站| 看免费av毛片| 国产精品影院久久| 精品免费久久久久久久清纯| 在线播放国产精品三级| 日本黄色视频三级网站网址| 免费高清在线观看日韩| 亚洲第一电影网av| 18禁观看日本| 国产一区二区激情短视频| 可以在线观看的亚洲视频| 看黄色毛片网站| 99国产极品粉嫩在线观看| 久久人妻福利社区极品人妻图片| 精品久久蜜臀av无| 精品少妇一区二区三区视频日本电影| 久久久久久大精品| 美女扒开内裤让男人捅视频| 欧美一区二区精品小视频在线| 国产精品爽爽va在线观看网站 | 亚洲av第一区精品v没综合| 大码成人一级视频| 久久久久久免费高清国产稀缺| 女生性感内裤真人,穿戴方法视频| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 日韩欧美三级三区| 久久久久久亚洲精品国产蜜桃av| 一级a爱视频在线免费观看| 日韩欧美三级三区| 在线av久久热| 激情视频va一区二区三区| 国产亚洲欧美在线一区二区| 熟女少妇亚洲综合色aaa.| 亚洲中文字幕日韩| 免费女性裸体啪啪无遮挡网站| 成熟少妇高潮喷水视频| 午夜精品国产一区二区电影| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av | 亚洲一卡2卡3卡4卡5卡精品中文| 咕卡用的链子| 国产私拍福利视频在线观看| 18禁黄网站禁片午夜丰满| 一本综合久久免费| 国产成年人精品一区二区| 可以在线观看的亚洲视频| 老司机深夜福利视频在线观看| 日本 av在线| 久久午夜亚洲精品久久| 日韩精品中文字幕看吧| 欧美激情极品国产一区二区三区| 亚洲国产中文字幕在线视频| 好男人电影高清在线观看| 午夜亚洲福利在线播放| 人人澡人人妻人| 自线自在国产av| 大型av网站在线播放| 国产在线精品亚洲第一网站| 女人高潮潮喷娇喘18禁视频| 嫁个100分男人电影在线观看| 久久国产亚洲av麻豆专区| 亚洲国产精品999在线| 一级,二级,三级黄色视频| 国产精品电影一区二区三区| 欧美在线黄色| 操出白浆在线播放| 国产精品一区二区精品视频观看| 免费高清在线观看日韩| 人人澡人人妻人| 人人妻人人澡欧美一区二区 | 少妇粗大呻吟视频| 老鸭窝网址在线观看| 亚洲欧美日韩高清在线视频| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 精品久久久精品久久久| 亚洲精品美女久久久久99蜜臀| 十分钟在线观看高清视频www| 国产精品一区二区精品视频观看| 成人亚洲精品一区在线观看| 制服人妻中文乱码| 免费少妇av软件| 禁无遮挡网站| 9热在线视频观看99| 99在线人妻在线中文字幕| 亚洲第一欧美日韩一区二区三区| 999久久久国产精品视频| 啪啪无遮挡十八禁网站| 欧美乱码精品一区二区三区| 巨乳人妻的诱惑在线观看| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 久久婷婷成人综合色麻豆| 成人亚洲精品一区在线观看| 两人在一起打扑克的视频| 国产亚洲精品av在线| 日本五十路高清| 女生性感内裤真人,穿戴方法视频| 久久久国产成人精品二区| 成人特级黄色片久久久久久久| 天天添夜夜摸| 夜夜躁狠狠躁天天躁| 涩涩av久久男人的天堂| 精品欧美国产一区二区三| 亚洲专区字幕在线| 黄片小视频在线播放| 女人精品久久久久毛片| 国产高清视频在线播放一区| www.www免费av| 99久久99久久久精品蜜桃| 一进一出好大好爽视频| 国产成人欧美在线观看| 午夜福利高清视频| 嫁个100分男人电影在线观看| 午夜福利影视在线免费观看| 老汉色av国产亚洲站长工具| 国产精品久久电影中文字幕| av中文乱码字幕在线| 又大又爽又粗| www日本在线高清视频| 啦啦啦 在线观看视频| 国产激情久久老熟女| 伊人久久大香线蕉亚洲五| 一本久久中文字幕| 757午夜福利合集在线观看| 欧美激情极品国产一区二区三区| 久久人人爽av亚洲精品天堂| 在线观看免费视频网站a站| 国产成人av教育| 黑人巨大精品欧美一区二区mp4| 久久精品亚洲精品国产色婷小说| 女生性感内裤真人,穿戴方法视频| 正在播放国产对白刺激| 国产xxxxx性猛交| 成年人黄色毛片网站| 欧美不卡视频在线免费观看 | 亚洲成人精品中文字幕电影| av网站免费在线观看视频| 亚洲国产精品成人综合色| 在线观看午夜福利视频| 黄色丝袜av网址大全| 久久国产精品影院| 丰满的人妻完整版| 给我免费播放毛片高清在线观看| 欧美一级a爱片免费观看看 | 好男人在线观看高清免费视频 | 欧美激情久久久久久爽电影 | 级片在线观看| 亚洲三区欧美一区| 长腿黑丝高跟| 国内精品久久久久久久电影| 色婷婷久久久亚洲欧美| av视频免费观看在线观看| 一二三四在线观看免费中文在| 国产不卡一卡二| 成人手机av| 精品久久久久久久人妻蜜臀av | 满18在线观看网站| 91老司机精品| 久久久久久人人人人人| or卡值多少钱| 国产精品久久电影中文字幕| 熟妇人妻久久中文字幕3abv| 91字幕亚洲| 在线永久观看黄色视频| 久久精品人人爽人人爽视色| 电影成人av| 夜夜夜夜夜久久久久| 一边摸一边抽搐一进一出视频| 男女做爰动态图高潮gif福利片 | 操美女的视频在线观看| 国产精品秋霞免费鲁丝片| 操出白浆在线播放| 国产欧美日韩综合在线一区二区| 一区在线观看完整版| tocl精华| 后天国语完整版免费观看| 成年人黄色毛片网站| 成人手机av| 国产av又大| 69精品国产乱码久久久| 99香蕉大伊视频| 青草久久国产| 大香蕉久久成人网| 亚洲成av人片免费观看| 欧美乱妇无乱码| 91国产中文字幕| 好男人在线观看高清免费视频 | 精品久久久久久久毛片微露脸| 琪琪午夜伦伦电影理论片6080| 19禁男女啪啪无遮挡网站| 国产黄a三级三级三级人| 老汉色av国产亚洲站长工具| 老汉色av国产亚洲站长工具| АⅤ资源中文在线天堂| 日韩三级视频一区二区三区| 日韩三级视频一区二区三区| 国产一区二区在线av高清观看| 亚洲精品国产精品久久久不卡| 国产午夜精品久久久久久| 琪琪午夜伦伦电影理论片6080| 视频区欧美日本亚洲| 嫩草影院精品99| 日韩欧美免费精品| av视频在线观看入口| 国产精品久久久久久人妻精品电影| 亚洲一码二码三码区别大吗| 亚洲国产欧美一区二区综合| 丝袜人妻中文字幕| 午夜福利视频1000在线观看 | 免费在线观看视频国产中文字幕亚洲| 超碰成人久久| av免费在线观看网站| 亚洲欧美日韩另类电影网站| 日韩中文字幕欧美一区二区| 中文亚洲av片在线观看爽| 色哟哟哟哟哟哟| 韩国精品一区二区三区| 午夜福利18| 国产不卡一卡二| 91麻豆精品激情在线观看国产| 99国产精品免费福利视频| 久久久久久久午夜电影| 极品教师在线免费播放| 满18在线观看网站| 十八禁人妻一区二区| 欧美日韩福利视频一区二区| 久久这里只有精品19| 日韩一卡2卡3卡4卡2021年| 精品久久蜜臀av无| 黄色丝袜av网址大全| 日韩大尺度精品在线看网址 | 一区二区日韩欧美中文字幕| 宅男免费午夜| 久久精品国产99精品国产亚洲性色 | 天天添夜夜摸| 久久久久精品国产欧美久久久| 国产成人av激情在线播放| 亚洲中文日韩欧美视频| 欧美性长视频在线观看| 一边摸一边做爽爽视频免费| 99久久久亚洲精品蜜臀av| 欧美日本亚洲视频在线播放| 一区福利在线观看| 亚洲成av片中文字幕在线观看| 免费高清在线观看日韩| 人成视频在线观看免费观看| 亚洲五月婷婷丁香| av片东京热男人的天堂| 亚洲国产中文字幕在线视频| 中文字幕久久专区| 天天一区二区日本电影三级 | or卡值多少钱| 夜夜爽天天搞| 久99久视频精品免费| 人妻久久中文字幕网| 50天的宝宝边吃奶边哭怎么回事| av电影中文网址| 亚洲五月色婷婷综合| 日韩精品青青久久久久久| 国产成人系列免费观看| 精品国产美女av久久久久小说| av有码第一页| 一级片免费观看大全| 国产色视频综合| 国产精品亚洲av一区麻豆| 女警被强在线播放| 免费不卡黄色视频| 久久人妻av系列| 亚洲va日本ⅴa欧美va伊人久久| 久久精品91蜜桃| 国产亚洲精品第一综合不卡| 国产91精品成人一区二区三区| 亚洲狠狠婷婷综合久久图片| 十八禁人妻一区二区| 久久久国产成人精品二区| 国产成人av激情在线播放| 少妇熟女aⅴ在线视频| 亚洲国产欧美日韩在线播放| 美女午夜性视频免费| 国产亚洲欧美在线一区二区| 色精品久久人妻99蜜桃| 两人在一起打扑克的视频| 亚洲在线自拍视频| 成人国产综合亚洲| 神马国产精品三级电影在线观看 | av中文乱码字幕在线| 少妇裸体淫交视频免费看高清 | 久久精品亚洲精品国产色婷小说| 美女午夜性视频免费| 日本一区二区免费在线视频| 亚洲中文日韩欧美视频| 亚洲国产欧美日韩在线播放| 黄色女人牲交| 两个人视频免费观看高清| 三级毛片av免费| 啦啦啦观看免费观看视频高清 | 久久久久久大精品| 真人一进一出gif抽搐免费| 伦理电影免费视频| 又黄又爽又免费观看的视频| 国产成年人精品一区二区| 午夜福利成人在线免费观看| 亚洲黑人精品在线| 麻豆国产av国片精品| 精品卡一卡二卡四卡免费| 欧美在线黄色| 日韩高清综合在线| 国产1区2区3区精品| 一二三四在线观看免费中文在| 一级黄色大片毛片| 欧美在线黄色| 日本免费a在线| 成人手机av| 一进一出抽搐动态| 欧美国产日韩亚洲一区| 丰满人妻熟妇乱又伦精品不卡| 女性被躁到高潮视频| 最好的美女福利视频网| 久久精品人人爽人人爽视色| 久久久国产成人免费| 人人澡人人妻人| 亚洲国产高清在线一区二区三 | 国产人伦9x9x在线观看| 亚洲国产精品合色在线| 精品久久久久久,| 99国产综合亚洲精品| x7x7x7水蜜桃| 亚洲av片天天在线观看| 一进一出抽搐gif免费好疼| 日本精品一区二区三区蜜桃| 非洲黑人性xxxx精品又粗又长| 亚洲国产欧美日韩在线播放| 日韩av在线大香蕉| 亚洲情色 制服丝袜| 少妇 在线观看| 色综合亚洲欧美另类图片| 人成视频在线观看免费观看| 成人亚洲精品av一区二区| 欧美国产精品va在线观看不卡| 国产男靠女视频免费网站| 女生性感内裤真人,穿戴方法视频| 久久久国产成人免费| 国产99久久九九免费精品| 国产一区二区三区视频了| 欧美老熟妇乱子伦牲交| 久久久久九九精品影院| 天堂影院成人在线观看| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 一级,二级,三级黄色视频| 正在播放国产对白刺激| 中文亚洲av片在线观看爽| 免费在线观看完整版高清| 亚洲成av人片免费观看| 中文字幕久久专区| 18禁观看日本| 露出奶头的视频| 两人在一起打扑克的视频| 中文字幕精品免费在线观看视频| 波多野结衣av一区二区av| 精品免费久久久久久久清纯| 成熟少妇高潮喷水视频| 国产精品一区二区精品视频观看| 丰满的人妻完整版| 午夜精品在线福利| 一级毛片精品| 亚洲一区二区三区不卡视频| 一级黄色大片毛片| 美女高潮喷水抽搐中文字幕| 国产熟女午夜一区二区三区| 久久国产亚洲av麻豆专区| 人人妻,人人澡人人爽秒播| 国产精品精品国产色婷婷| 亚洲第一欧美日韩一区二区三区| 手机成人av网站| 亚洲激情在线av| 99在线视频只有这里精品首页| 一区二区三区国产精品乱码| 国产三级在线视频| 免费搜索国产男女视频| 最新美女视频免费是黄的| 老熟妇仑乱视频hdxx| 欧美日韩亚洲综合一区二区三区_| 精品日产1卡2卡| 精品不卡国产一区二区三区| 亚洲久久久国产精品| bbb黄色大片| 国产精品乱码一区二三区的特点 | 久久九九热精品免费| 成人精品一区二区免费| 97人妻精品一区二区三区麻豆 | 黑人巨大精品欧美一区二区mp4| 国产精品一区二区三区四区久久 | 涩涩av久久男人的天堂| 亚洲aⅴ乱码一区二区在线播放 | 国产一级毛片七仙女欲春2 | 亚洲国产精品999在线| 午夜福利在线观看吧| 中文字幕久久专区| 色在线成人网| a在线观看视频网站| 一级毛片高清免费大全| 精品国产一区二区三区四区第35| 十八禁人妻一区二区| 中文字幕高清在线视频| 啦啦啦免费观看视频1| 久久性视频一级片| 欧美乱码精品一区二区三区| 久久婷婷人人爽人人干人人爱 | 成人国产综合亚洲| 国产欧美日韩精品亚洲av| 欧美黄色淫秽网站| 少妇 在线观看| 成人亚洲精品一区在线观看| 精品久久久久久成人av| 午夜两性在线视频| 在线永久观看黄色视频| 久9热在线精品视频| 757午夜福利合集在线观看| 老汉色av国产亚洲站长工具| 91麻豆av在线| 久久久久久人人人人人| 啪啪无遮挡十八禁网站| 久久国产亚洲av麻豆专区| 在线观看免费视频网站a站| 亚洲人成电影免费在线| 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区 | 18禁美女被吸乳视频| 精品国产乱子伦一区二区三区| 亚洲美女黄片视频| 午夜精品在线福利| 国产精品国产高清国产av| 一二三四社区在线视频社区8| 91麻豆av在线| 天天躁夜夜躁狠狠躁躁| 麻豆av在线久日| 国产精品亚洲美女久久久| 十八禁网站免费在线| 91麻豆av在线| 嫩草影视91久久| 午夜免费激情av| 在线十欧美十亚洲十日本专区| 亚洲精品粉嫩美女一区| 午夜成年电影在线免费观看| 久久久久亚洲av毛片大全| 十八禁人妻一区二区| 在线十欧美十亚洲十日本专区| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品久久男人天堂| 神马国产精品三级电影在线观看 | 久久精品亚洲熟妇少妇任你| 久久久久久久精品吃奶| 久久国产精品影院| 视频在线观看一区二区三区| 变态另类成人亚洲欧美熟女 | 欧美绝顶高潮抽搐喷水| 成在线人永久免费视频| 亚洲国产日韩欧美精品在线观看 | 国产亚洲欧美精品永久| 人人妻人人澡欧美一区二区 | 国产成人av教育| 一区二区三区激情视频| 亚洲黑人精品在线| 国产成+人综合+亚洲专区| videosex国产| 国产99久久九九免费精品| 宅男免费午夜| 多毛熟女@视频| 日韩成人在线观看一区二区三区| bbb黄色大片| 老司机靠b影院| 国产亚洲av嫩草精品影院| 一级黄色大片毛片|