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

    秦嶺南坡中段主要森林群落類型劃分及環(huán)境梯度解釋

    2016-09-19 06:38:46王宇超周亞福王得祥陜西省西安植物園陜西西安7006西北農(nóng)林科技大學(xué)林學(xué)院陜西楊凌7200
    關(guān)鍵詞:樣地坡度速效

    王宇超,周亞福,王得祥. 陜西省西安植物園,陜西 西安 7006;2. 西北農(nóng)林科技大學(xué)林學(xué)院,陜西 楊凌 7200

    秦嶺南坡中段主要森林群落類型劃分及環(huán)境梯度解釋

    王宇超1,周亞福1,王得祥2*
    1. 陜西省西安植物園,陜西 西安 710061;2. 西北農(nóng)林科技大學(xué)林學(xué)院,陜西 楊凌 712100

    該研究旨在揭示秦嶺大熊貓主要棲息地植物群落分布規(guī)律及其與環(huán)境的對(duì)應(yīng)關(guān)系。在研究區(qū)設(shè)置2條樣帶,每條樣帶各設(shè)置40~50個(gè)樣地,共計(jì)93個(gè)。通過(guò)樣地調(diào)查,利用雙向聚類法(Two-way Clustering Method)、DCA、DCCA對(duì)秦嶺大熊貓(Ailuropoda melanoleuca)棲息地森林群落進(jìn)行數(shù)量分類和排序。(1)TWINSPAN將研究區(qū)93個(gè)樣地分成太白紅杉(Larix chinensis)、巴山冷杉(Abies fargesii)、鐵杉(Tsuga chinensis)、華山松(Pinus armandii)、油松(Pinus tabulaeformi)、銳齒槲櫟(Quercus aliena var. acuteserrata)等34個(gè)群落類型,并通過(guò)去趨勢(shì)對(duì)應(yīng)分析(DCA)對(duì)93個(gè)樣地進(jìn)行排序,結(jié)果表明不同類型樣地呈現(xiàn)聚集分布,充分驗(yàn)證了TWINSPAN對(duì)群落類型劃分的科學(xué)性。(2)通過(guò)對(duì)93個(gè)樣地除趨勢(shì)典范對(duì)應(yīng)分析(DCCA)和Monte Carlo顯著性檢驗(yàn)表明,在分析12個(gè)環(huán)境因子中,制約秦嶺南坡大熊貓棲息地群落類型、植物種分布格局的主要因素是海拔、坡度、速效磷與全磷。海拔、坡度、速效磷和全磷能夠很好地解釋群落樣地及優(yōu)勢(shì)種在環(huán)境梯度上的分布,同時(shí)表現(xiàn)出樣方間在種類和環(huán)境因子組成上的相似性。

    秦嶺;南坡;雙向聚類法;去趨勢(shì)對(duì)應(yīng)分析;除趨勢(shì)典范對(duì)應(yīng)分析

    引用格式:王宇超, 周亞福, 王得祥. 秦嶺南坡中段主要森林群落類型劃分及環(huán)境梯度解釋[J]. 生態(tài)環(huán)境學(xué)報(bào), 2016, 25(6): 965-972.

    WANG Yuchao, ZHOU Yafu, WANG Dexiang. The Quantitative Classification and Environmental Interpretation of Forest Communities in the Middle Area of South Slope of Qinling Mountains [J]. Ecology and Environmental Sciences, 2016, 25(6): 965-972.

    目前,國(guó)內(nèi)針對(duì)山地植被進(jìn)行的研究,無(wú)論是從宏觀尺度還是微觀研究,都開(kāi)展得比較深入和透徹。其中,山地植物群落格局與環(huán)境關(guān)系研究是當(dāng)前植被研究的熱點(diǎn)問(wèn)題之一。許多研究認(rèn)為植被分布地形環(huán)境(徐廣平等,2005;張先平等,2006)、群落的演替階段(王世雄等,2010)、外界干擾程度(李晉鵬等,2007)、地表和植被狀況(王國(guó)宏等,2001;張玲等,2007)、土壤狀況(龔志蓮等,2015)等因素都可能對(duì)植物群落分布格局產(chǎn)生影響。任何植物群落的存在,都與其所在的環(huán)境條件有著密切的相關(guān)性。因此,清楚認(rèn)識(shí)植物群落分布格局和生態(tài)環(huán)境之間的關(guān)系是植被生態(tài)研究的基礎(chǔ)與前提。然而,數(shù)量分析為客觀、準(zhǔn)確地群落類型劃分及其與環(huán)境之間的生態(tài)關(guān)系提供了科學(xué)、合理、有效的途徑,己成為植被生態(tài)學(xué)研究的重要內(nèi)容之一(張新時(shí),1991;江洪等,1994)。其中,數(shù)量分類和排序是最主要和應(yīng)用最廣泛的多元統(tǒng)計(jì)方法,是全面認(rèn)識(shí)研究區(qū)域內(nèi)植被分布格局特征的重要手段(張峰等,2003)。目前,數(shù)量分類和排序已被較為廣泛地應(yīng)用在森林(Veazquez,1994;王晶等,2016)、濕地(李思陽(yáng)等,2016)、草甸(宋愛(ài)云等,2006)、灌叢(孫菊等,2009)等群落類型研究之中。數(shù)量分類和排序已經(jīng)在植被生態(tài)學(xué)中發(fā)展成為研究群落生態(tài)關(guān)系的重要數(shù)量方法(Russelt,1991;Pamer,1993)。

    佛坪和長(zhǎng)青國(guó)家級(jí)自然保護(hù)區(qū)位于秦嶺中段南坡,是秦嶺重要核心區(qū)域,同時(shí)也是大熊貓Ailuropoda melanoleuca等珍稀動(dòng)物的主要棲息地。其植被空間格局直接影響到大熊貓等野生動(dòng)物生存空間,所以對(duì)該研究區(qū)植物群落進(jìn)行研究顯得尤為重要。然而,近些年關(guān)于該研究區(qū)植被分布格局及其與環(huán)境關(guān)系的研究還較少,對(duì)植被分布格局及其與環(huán)境關(guān)系的認(rèn)識(shí)還不是很清楚。本研究結(jié)合我國(guó)山地植物群落數(shù)量分類和排序的相關(guān)研究,利用TWINSPAN數(shù)量分類和DCA、DCCA排序,全面深入探討秦嶺大熊貓棲息地地形、土壤等環(huán)境因子與山地植被類型及其分布的關(guān)系,并對(duì)植被類型作出環(huán)境解釋,以期為保護(hù)棲息地的森林植被,充分發(fā)揮其強(qiáng)大的生態(tài)系統(tǒng)服務(wù)功能,實(shí)現(xiàn)保護(hù)區(qū)森林生態(tài)系統(tǒng)可持續(xù)發(fā)展和管理提供理論指導(dǎo),為制定保護(hù)措施提供科學(xué)依據(jù)。

    1 研究區(qū)概況與研究方法

    1.1研究區(qū)概況

    佛坪、長(zhǎng)青國(guó)家級(jí)自然保護(hù)區(qū)位于秦嶺中段南坡,地理位置為東經(jīng)107°19′~107°55′、北緯33°17′~33°46′,總面積59146 hm2,森林覆蓋率達(dá)90%以上。研究區(qū)為亞熱帶向暖溫帶的過(guò)渡區(qū),海拔在900~3071 m之間,年均溫為11.4~15 ℃,無(wú)霜期為225~239 d,年降水量813~1129.6 mm。研究區(qū)主要以保護(hù)大熊貓、朱鹮 Nipponia nippon、羚牛Budorcas taxicolor、金絲猴 Rhinopithecus roxellana為主。該區(qū)植被類型豐富,主要有針葉林、闊葉林、竹林、草甸等植被類型。主要分布植物群落優(yōu)勢(shì)種有太白紅杉Larix chinensis、巴山冷杉Abies fargesii、華山松Pinus armandii、油松Pinus tabulaeformis、銳齒槲櫟 Quercus aliena var. acutiserrata、巴山木竹Bashania fargesii、秦嶺箭竹Fargesia qinlingensis等。

    1.2研究方法

    1.2.1樣地設(shè)置

    在查閱相關(guān)資料和踏查的基礎(chǔ)上,采用典型取樣法對(duì)研究區(qū)內(nèi)海拔900~3071 m之間的山地喬木植物群落進(jìn)行野外調(diào)查。在每個(gè)保護(hù)區(qū)內(nèi)各設(shè)置 1條樣帶,在每條樣帶上按海拔每升高 100 m設(shè)置1~4個(gè)典型樣地。2條樣帶共設(shè)置93個(gè)樣地。

    喬木樣地:面積500 m2(20 m×25 m),對(duì)樣地內(nèi)喬木層(胸徑大于4 cm)植株進(jìn)行每木檢尺,在每個(gè)樣地內(nèi)設(shè)置5個(gè)5 m×5 m的灌木樣方和1 m×1 m的草本樣方,并在每個(gè)樣地沿對(duì)角線設(shè)置25個(gè)1 m×1 m灌草頻度調(diào)查樣方。所有樣地調(diào)查記錄內(nèi)容包括植物種類、高度、蓋度、胸徑、頻度、生活力等指標(biāo),同時(shí)記錄各群落的綜合特征和生境特征,主要包括:經(jīng)緯度、海拔、群落蓋度、層間植物種類、坡位、坡度、坡向和腐殖質(zhì)層厚度等。

    1.2.2土壤采集及測(cè)定

    在每個(gè)樣地四角及中央設(shè)置5個(gè)土壤剖面,每個(gè)剖面分3層取土。分別是0~20、20~40、40~60 cm,帶回室內(nèi)分析其養(yǎng)分含量。測(cè)試指標(biāo)有:全氮、全磷、全鉀、有效氮、有效磷、有效鉀、有機(jī)質(zhì)、pH值。土樣處理以及測(cè)試方法參照鮑士旦等(1999)實(shí)驗(yàn)方法進(jìn)行。

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

    (1)物種重要值計(jì)算(王宇超等,2012)

    喬木重要值=(相對(duì)高度+相對(duì)密度+相對(duì)優(yōu)勢(shì)度)/300;

    灌木及草本重要值=(相對(duì)密度+相對(duì)頻度+相對(duì)蓋度)/300;

    (2)植物群落數(shù)量分類與排序

    本文采用雙向聚類法TWINSPAN(Two-Way Indicators Species Analysis)對(duì)群落進(jìn)行分類,去趨勢(shì)對(duì)應(yīng)分析 DCA(Detrended Correspondence Analysis)、除趨勢(shì)典范對(duì)應(yīng)分析(Detrended Canonical Correspondence Analysis,DCCA)進(jìn)行排序。

    1.2.4坡度、坡向的標(biāo)準(zhǔn)化

    野外測(cè)量時(shí),坡度以度為單位記錄(0°~90°)。坡向的數(shù)據(jù)用方位角表示(0°~360°)。為了便于分析,將坡度和坡向通過(guò)如下公式進(jìn)行標(biāo)準(zhǔn)化:式中,SP為百分比坡度,S為坡度,SP的取值范圍為(0, 200);AT為標(biāo)準(zhǔn)化坡向,A為坡向,AT的取值范圍為(0, 2),最大值2代表濕冷的東北坡,最小值0則代表干熱的西南坡。

    1.2.5數(shù)據(jù)統(tǒng)計(jì)與分析

    數(shù)據(jù)分析利用TWINSPAN 2.3、Canoco 4.5、Excel 2007等軟件。

    2 結(jié)果分析

    2.1TWINSPAN數(shù)量分類

    基于樣地調(diào)查數(shù)據(jù)生成的物種重要值矩陣,利用雙向聚類法(Two-way Clustering Method)對(duì)93個(gè)樣地進(jìn)行數(shù)量分類,該方法采用TWINSPAN軟件包進(jìn)行。本文根據(jù)物種重要值選用的假種分級(jí)為7級(jí):0.02~0.05,0.05~0.1,0.1~0.15,0.15~0.2,0.2~0.35,0.35~0.7,>0.7。列表中的最大物種數(shù)為120個(gè),用來(lái)劃分每一組樣地個(gè)數(shù)最小值為 5,最大劃分分級(jí)水平為 9,每次劃分的最多區(qū)別中數(shù)目為8。樣方分類時(shí)為了突出優(yōu)勢(shì)種的作用,對(duì)7級(jí)“假種”的權(quán)重分別?。?、1、1、1、2、3和4。

    93個(gè)樣地 658種植物進(jìn)行數(shù)量分類結(jié)果見(jiàn)圖1。其中,圖 1的下方和右方分別顯示對(duì)樣方和物種的分化水平和類型。從圖1可以看出TWINSPAN 將93個(gè)樣地劃分成39個(gè)群落類型。由于本研究只是從植物群系方面進(jìn)行劃分,所以結(jié)合野外調(diào)查數(shù)據(jù),對(duì)一些劃分過(guò)細(xì)的樣地進(jìn)行合并。樣地46、65、17、18、19、25的建群種均為銳齒槲櫟,將其合并成為一個(gè)群系;樣地12、15、24、69、45、58、51的建群種均為油松,將其合并成一個(gè)群系,樣地59、10、11、26、28、29、53建群種為板栗,將其合并成一個(gè)群系。通過(guò)合并將 93個(gè)樣地劃分成太白紅杉、巴山冷杉、油松、華山松等 34個(gè)群系,具體名錄如表1所示。

    表1 秦嶺南坡中段主要喬木群落類型Table 1 Tree community types in the middle area of south slope of Qinling Mountains

    圖1 秦嶺南坡中段植物群落TWINSPAN分類結(jié)果Fig. 1 The TWINSPAN classification for forest communities

    2.2植物群落DCA排序

    2.2.1森林植物群落DCA排序

    利用DCA對(duì)秦嶺南坡中段93個(gè)森林喬木群落樣方進(jìn)行排序,由表2可知,4個(gè)排序軸的前兩個(gè)軸特征值較大且與環(huán)境有著較高相關(guān)性,具有一定生態(tài)意義,而后兩個(gè)排序軸特征值較小且與環(huán)境之間相關(guān)性甚微,其生態(tài)意義不明確。所以本文根據(jù)前兩個(gè)排序軸作出二維排序圖(圖2),由圖2可知,TWINSPAN分類所得到各個(gè)群落類型在二維的排序圖上各自有一定的分布范圍,基本都呈現(xiàn)出相同植物群落類型聚集在同一個(gè)范圍內(nèi),DCA排序結(jié)果基本印證TWINSPAN的分類結(jié)果。這也進(jìn)一步說(shuō)明了DCA能夠較好地反映各個(gè)群落之間以及群落與環(huán)境之間的關(guān)系。

    從排序軸來(lái)看,第一軸基本上反映各植物群落所在環(huán)境的海拔梯度變化,即DCA第一軸從左到右,海拔逐漸增大;第二軸基本上反映各植物群落所在環(huán)境的濕度變化,即DCA第二軸從上到下,林下濕度逐步增加。各個(gè)群落在排序軸上呈現(xiàn)規(guī)律分布,林下環(huán)境較為干燥、分布海拔較低的栓皮櫟Quercus variabilis植物群落位于排序軸的左上方,林下環(huán)境極其濕潤(rùn),分布在溝谷中山楠群落位于排序軸下方,遠(yuǎn)離其它群落類型。水杉 Metasequoia glyptostroboides林是人工林,分布海拔較低,常常栽植在河道兩側(cè)空曠地上,其林下環(huán)境較為干燥,所以水杉林和山核桃Carya cathayensis 群落位于排序軸的上方。巴山冷杉、牛皮樺Betula albo-sinensis var. septentrionalis、太白紅杉高海拔分布植物群落位于排序軸右側(cè)。

    表2 群落DCA排序軸特征值及種-環(huán)境相關(guān)系數(shù)Table 2 Characteristic values and species-environment correlation coefficients on DCA axis forest communities

    圖2 秦嶺南坡大熊貓棲息地森林植物群落93個(gè)樣地DCA二維散布圖Fig. 2 A two-dimensional scatter plot of DCA ordination for forest communities

    從圖2可知闊雜、油松、銳齒槲櫟、白樺、板栗、短柄枹櫟Quercus glandulifera分布在林下中度濕度、中低海拔區(qū)域,位于排序軸中間靠左的位置,各個(gè)群落類型分布較為密集,相互之間有一定交錯(cuò)分布。說(shuō)明這幾個(gè)群落類型分布生境相近,在物種組成上有一定的相似性,該環(huán)境區(qū)域是整個(gè)研究區(qū)植物群落類型分布最多,物種最多的區(qū)域。

    2.3植物群落分布與環(huán)境之間的關(guān)系

    2.3.1群落樣地分布與環(huán)境因子的DCCA排序

    本研究以秦嶺南坡中段主要喬木植物群落物種組成和土壤、地形指標(biāo)所組成的環(huán)境因子矩陣進(jìn)行DCCA分析,并與DCA排序結(jié)果進(jìn)行分析對(duì)比。

    從表3可知,第1軸特征值最大(0.701),第2軸次之(0.359),而其它 2個(gè)軸的特征值較小,說(shuō)明前兩個(gè)軸包含較多信息量,故采用第1軸、第2軸數(shù)據(jù)作DCCA二維排序圖(圖3)。圖中箭頭代表各個(gè)環(huán)境因子,這是DCCA排序的特點(diǎn)。箭頭所處象限,代表著環(huán)境因子與排序軸間的正負(fù)相關(guān)性;箭頭連線的長(zhǎng)短表示植物群落的分布與該環(huán)境因子相關(guān)性的大小,箭頭連線在排序中的斜率表示環(huán)境因子與排序軸相關(guān)性的大小,箭頭所處的象限表示環(huán)境因子與排序軸之間相關(guān)性的正負(fù)。

    表3 DCCA排序軸特征值及種-環(huán)境因子相關(guān)系數(shù)Table 3 Characteristic values and species-environment correlation coefficients on DCCA axis

    從圖3和表4中可以看出,第一軸基本上反映各植物群落樣地所在環(huán)境的海拔梯度變化趨勢(shì)(海拔與DCCA第一軸典范系數(shù)0.9648,相關(guān)系數(shù)為0.9925),即沿DCCA第1軸從左到右,海拔逐步增加;第二軸基本上表現(xiàn)出植物群落所在環(huán)境的坡度、全磷、速效磷的變化趨勢(shì)(全磷與第二軸典范系數(shù)為0.6610,相關(guān)系數(shù)為0.7421;速效磷與第二軸的典范系數(shù)為0.5760,相關(guān)系數(shù)為0.6466;坡度與第二軸的典范系數(shù)為-0.4459,相關(guān)系數(shù)為-0.5006),即沿DCCA第二軸從下到上,坡度漸緩、坡向越向陽(yáng),土壤全磷、速效磷含量逐步增加。因此,海拔、全磷、速效磷和坡度是本研究 12個(gè)環(huán)境因子中影響植物群落分布的幾個(gè)主導(dǎo)因子。

    表4 環(huán)境因子與DCCA排序軸的典范系數(shù)和相關(guān)系數(shù)Table 4 Canonical coefficients and correlation coefficients of environmental factors with DCA axis

    圖3 93個(gè)樣地與12個(gè)環(huán)境因子的DCCA二維排序圖Fig. 3 Two-dimensional DCCA ordination diagram of 93 samples with 12 environmental factors

    2.3.2主要優(yōu)勢(shì)種分布與環(huán)境因子的DCCA排序

    通過(guò)分析主要優(yōu)勢(shì)種的DCCA二維排序圖(圖4)發(fā)現(xiàn),群落優(yōu)勢(shì)種分布格局所揭示的環(huán)境梯度與樣地群落的分布梯度有著很大相似性。沿著Axis1軸從左至右,喬木層中分布在低海拔的栓皮櫟、廟臺(tái)槭 Acer miaotaiense、楓楊 Pterocarya stenoptera等物種位于排序軸最左側(cè)。隨著海拔進(jìn)一步升高,進(jìn)入整個(gè)研究區(qū)中海拔區(qū)域,從圖4也可以看出這個(gè)區(qū)域優(yōu)勢(shì)種數(shù)量明顯高于兩端,該區(qū)域分布喬木有油松、銳齒槲櫟、白樺、亮葉樺Betula luminifera等,這些樹(shù)種分布在各個(gè)環(huán)境因子中心位置,對(duì)環(huán)境要求相對(duì)較低,所以它們是整個(gè)研究區(qū)分布最為廣泛的幾個(gè)樹(shù)種。高海拔太白紅杉、巴山冷杉位于排序圖的最右側(cè),與其它喬木層相隔較遠(yuǎn),說(shuō)明其生境極其嚴(yán)酷,以致其它物種無(wú)法在此生長(zhǎng)。灌木層和草本層也表現(xiàn)出與喬木層相同的變化規(guī)律。其中,巴山木竹和秦嶺箭竹作為大熊貓的主要食材,它們分別分布在中海拔區(qū)域和高海拔區(qū)域,這與現(xiàn)實(shí)調(diào)查結(jié)果相吻合。通過(guò)上面分析說(shuō)明,Axis1軸從左至右,海拔逐漸增大,環(huán)境溫度逐步降低,物種對(duì)熱量需求逐步減少。從圖4也可以看出,Axis2軸與土壤速效磷以及坡度有著密切相關(guān)關(guān)系,從上至下,坡度逐步增加,土壤速效磷逐步減少。山核桃、白樺、太白楊等位于較平緩山地上,土壤中速效磷的含量較為豐富。秦嶺冷杉 Abies chensiensis、鐵杉 Tsuga chinensis、杈葉槭 Acer robustum等位于坡度較大山地上,土壤速效磷含量較低。通過(guò)分析圖4發(fā)現(xiàn),DCCA二維排序圖所顯示的植物群落優(yōu)勢(shì)種分布情況與實(shí)際調(diào)查記錄相一致,說(shuō)明DCCA二維排序圖能較好反應(yīng)主要優(yōu)勢(shì)種與環(huán)境梯度之間的關(guān)系。因此,海拔、坡度和速效磷是影響優(yōu)勢(shì)種分布格局的主要因素,這與群落樣地分布和環(huán)境梯度關(guān)系相似。

    圖4 主要優(yōu)勢(shì)種與12個(gè)環(huán)境因子的DCCA二維排序圖Fig. 4 Two-dimensional DCCA ordination diagram of main dominant species with 12 environmental factors

    2.3.3Monte Carlo檢驗(yàn)

    從DCCA第一排序軸Monte Carlo檢驗(yàn)結(jié)果(表5)可以看出,無(wú)論是對(duì)于群落樣地,還是對(duì)于優(yōu)勢(shì)種分布,海拔、坡度、速效磷和全磷的F值都達(dá)到了顯著水平(P<0.01),因此,海拔、坡度、速效磷和全磷可以很好地解釋植物群落樣地和優(yōu)勢(shì)物種的分布。該檢驗(yàn)結(jié)果充分證明前面二維排序圖所顯示的海拔、坡度、速效磷和全磷是決定植物群落和優(yōu)勢(shì)種分布的主要影響因素。

    3 結(jié)論與討論

    分別采用雙向聚類法(TWINSPAN)、去趨勢(shì)對(duì)應(yīng)分析(DCA)和除趨勢(shì)典范對(duì)應(yīng)分析(DCCA)對(duì)秦嶺南坡中段主要喬木群落進(jìn)行分類與排序,結(jié)果表明TWINSPAN雙向聚類分析較好地將不同的群落類型進(jìn)行分類,研究區(qū) 93個(gè)群落樣方被分成34個(gè)群落類型。本文利用去趨勢(shì)對(duì)應(yīng)分析(DCA)對(duì)樣地進(jìn)行排序,結(jié)果不同類型樣地呈現(xiàn)聚集分布,其結(jié)果充分驗(yàn)證了TWINSPAN對(duì)群落類型劃分的科學(xué)性。利用DCA對(duì)TWINSPAN分類結(jié)果進(jìn)行驗(yàn)證,在以往研究中都得到普遍體現(xiàn)(王得祥等,2009)。

    表5 環(huán)境因子與DCCA第一排序軸的Monte Carlo顯著性檢驗(yàn)Table 5 Monte Carlo permutation tests of significance of selected environmental variables with the first canonical axis from DCCA ordinations

    從廣泛的意義上講,植物種、植物群落的分布格局是在不同尺度上各種氣候、地形、土壤等環(huán)境因子綜合表現(xiàn)的結(jié)果。在區(qū)域尺度上講,氣候、植物區(qū)系決定特定的植被類型;景觀尺度上,如海拔這個(gè)因子,使局部環(huán)境的水熱條件發(fā)生改變,導(dǎo)致植被類型發(fā)生變化,從而決定植被分布格局;而微生境、土壤養(yǎng)分等決定最終的植物群落類型(Burke,2001)。本文結(jié)果表明,在景觀和土壤養(yǎng)分這兩個(gè)尺度上,制約秦嶺南坡大熊貓棲息地群落類型、植物種分布格局的主要因素是海拔、坡度、速效磷與全磷。海拔與坡度作為光、熱、水條件綜合表現(xiàn),其格局變化也就是光、熱、水的空間再次分配,從而使得植物群落組成以及植物種的分布格局發(fā)生變化(張峰等,2003;代力民等,2008)。速效磷和全磷作為土壤營(yíng)養(yǎng)元素重要成分之一,對(duì)植物生長(zhǎng)有著重要作用,直接影響著植物群落的組成與生理活力,決定著生態(tài)系統(tǒng)的結(jié)構(gòu)、功能和生產(chǎn)力水平(Robertson et al.,1981;Vitousek et al.,1989)。綜上,海拔、坡度、速效磷和全磷格局變化是光、熱、水、土壤營(yíng)養(yǎng)的格局變化,從而導(dǎo)致植物群落類型及物種分布格局發(fā)生變化。

    本文通過(guò)Monte Carlo顯著性檢驗(yàn)表明,海拔、坡度、速效磷和全磷能夠很好地解釋群落樣地及優(yōu)勢(shì)種在環(huán)境梯度上的分布。同時(shí)發(fā)現(xiàn)DCA、DCCA排序圖不僅反映樣方間在種類組成上的相似性,而且也反映樣方間在環(huán)境因子組成上的相似性,而這兩種相似性往往相互聯(lián)系(代力民等,2008)。一般來(lái)講,種類組成接近的植物群落,其環(huán)境因子組成也較接近,這是由植物種、植物群落和環(huán)境因子之間相互作用的生態(tài)關(guān)系所決定的。

    BURKE A. 2001. Classification and ordination of plant communities of the Naukluft Mountains Namibia [J]. Journal of Vegetation Science, 12(1): 53-60.

    PAMER M W. 1993. Putting things in even better order: The advantages of canonical correspondence analysis [J]. Ecology, 74(8): 2215-2230.

    ROBERTSON G P, VITOUSEK P M. 1981. Nitrification potentials in primary and secondary Succession [J]. Ecology, 62(2): 376-387.

    RUSSELT SMITH J. 1991. Classification, species richness, and environmental relations of monsoon rain forest in northem Australia [J]. Journal of Vegetation Science, 2(2): 259-278.

    VEAZQUEZ A. 1994. Multivariate analysis of the vegetation of the volcanoes Taloc and Pelado, Mexico [J]. Journal of Vegetation Science,5(2): 263-270.

    VITOUSEK P M, MATSON P A, VANCKVE K. 1989. Nitrogen availability and nitrification during succession, primary, secondary and old field seres [J]. Plant and Soil, 115(2): 229-239.

    代力民, 唐立娜, 曹玉明, 等. 2008. 遼東山區(qū)生態(tài)土地分類中的植物群落數(shù)量分析[J]. 林業(yè)科學(xué), 44(3): 6-12.

    鮑士旦, 江榮風(fēng), 楊超光, 等. 1999. 土壤農(nóng)化分析[M]. 北京: 中國(guó)農(nóng)業(yè)出版社: 14-114.

    龔志蓮, 李勇. 2015. 川西南不同植被木本物種多樣性與土壤有機(jī)碳積累[J]. 生態(tài)環(huán)境學(xué)報(bào), 24(9): 1447-1452.

    江洪, 黃建輝, 陳靈芝, 等. 1994. 東靈山植物群落的排序、數(shù)量分類與環(huán)境解釋[J]. 植物學(xué)報(bào), 36(7): 539-551.

    李晉鵬, 上官鐵梁, 孟東平, 等. 2007. 山西呂梁山南段植物群落的數(shù)量分類和排序研究[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 13(5): 615-619.

    李思陽(yáng), 張娟, 姚玲愛(ài), 等. 2016. 西枝江流域浮游植物群落結(jié)構(gòu)特征與主要環(huán)境因子的關(guān)系研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 36(6): 1939-1947.

    宋愛(ài)云, 劉世榮, 史作, 等. 2006. 臥龍自然保護(hù)區(qū)亞高山草甸植物群落物種多樣性研究[J]. 林業(yè)科學(xué)研究, 19(6): 767-772.

    王得祥, 藺雨陽(yáng), 雷瑞德, 等. 2009. 秦嶺山地天然油松林群落結(jié)構(gòu)特征和數(shù)量分類研究[J]. 西北植物學(xué)報(bào), 29(5): 1016-1025.

    王國(guó)宏, 楊利民. 2001. 祁連山北坡中段森林植被梯度分析及環(huán)境解釋[J]. 植物生態(tài)學(xué)報(bào), 25(6): 733-740.

    王晶, 張欽弟, 許強(qiáng), 等. 2016. 山西龐泉溝銀露梅群落物種多樣性的海拔格局[J]. 植物學(xué)報(bào), 51(3): 335-342.

    王世雄, 王孝安, 李國(guó)慶, 等. 2010. 陜西子午嶺植物群落演替過(guò)程中物種多樣性變化與環(huán)境解釋[J]. 生態(tài)學(xué)報(bào), 30(6): 1638-1647.

    王宇超, 王得祥, 胡有寧, 等. 2012. 陜西黑河上游主要天然林類型及物種多樣性特征研究[J]. 西北農(nóng)林科大學(xué)學(xué)報(bào)(自然科學(xué)版), 40(7): 106-112.

    孫菊, 李秀珍, 胡遠(yuǎn)滿, 等. 2009. 大興安嶺溝谷凍土濕地植物群落分類、物種多樣性和物種分布梯度[J]. 應(yīng)用生態(tài)學(xué)報(bào), 20(9): 2049-2056.徐廣平, 張德罡, 蒲小鵬, 等. 2005. 東祁連山高寒草地植物群落二維極點(diǎn)排序[J]. 草原與草坪, 3: 42-45.

    張峰, 張金屯. 2003. 歷山自然保護(hù)區(qū)豬尾溝森林群落植被格局及環(huán)境解釋[J]. 生態(tài)學(xué)報(bào), 23(3): 421-427.

    張玲, 袁曉穎, 張東來(lái). 2007. 大小興安嶺過(guò)渡區(qū)木本植物群落數(shù)量分類與排序[J]. 東北林業(yè)大學(xué)學(xué)報(bào), 35(9): 49-51.

    張先平, 王孟本, 佘波, 等. 2006. 龐泉溝國(guó)家自然保護(hù)區(qū)森林群落的數(shù)量分類和排序[J]. 生態(tài)學(xué)報(bào), 26(3): 754-761.

    張新時(shí). 1991. 西藏阿里植物群落的間接梯度分析、數(shù)量分類與環(huán)境解釋[J]. 植物生態(tài)學(xué)報(bào), 15(2): 101-113.

    The Quantitative Classification and Environmental Interpretation of Forest Communities in the Middle Area of South Slope of Qinling Mountains

    WANG Yuchao1, ZHOU Yafu1, WANG Dexiang2*
    1. Xi'an Botanical Garden of Shaanxi Province, Xi'an 710061, China;2. College of forestry, Northwest A & F Unversity, Yangling 712100, China

    In order to revealed the plant community distribution pattern and the corresponding relation with environment. Tow line transects were set between the altitude of 900~3 071 m in the field, and 93 standard plots (25 m×20 m) were investigated on transects. Vegetation community, growth status and distribution environment were recorded and soil samples were collected in the standard plots. Based on analysis of plots and soil data, the main forests in giant panda habitat were classified and ordering by using the two-way indicator species analysis (TWINSPAN), Detrended Corresponding Analysis (DCA) and Detrended Canonical Correspondence Analysis (DCCA). The results showed that: (1) 93 standard plots were classified into 34 plant communities by Two-way Indicators Species Analysis (TWINSPAN). The result of ordering by DCA showed that the standard plots fit to aggregation distribution and verified the scientificity of TWINSPAN. (2) The restrictive factors of community type and species distribution were altitude, slope, rapidly-available phosphorus and total phosphorus at the landscape level by the test of significance between DCCA and Monte Carlo. Factors of altitude, slope, rapidly-available phosphorus and total phosphorus can explain distribution of communities and dominant species at the landscape level, and the species and environmental factors show great similarities between different plots.

    Qinling Mountains; south slope; TWINSPAN; DCA; DCCA

    10.16258/j.cnki.1674-5906.2016.06.008

    Q948; X171.1

    A

    1674-5906(2016)06-0965-08

    陜西省科學(xué)院科技計(jì)劃項(xiàng)目(2014k-25);陜西省科學(xué)院科技重點(diǎn)項(xiàng)目(2013K-01);國(guó)家科技基礎(chǔ)性工作專項(xiàng)重點(diǎn)項(xiàng)目(2007FY110800)

    王宇超(1978年生),男,博士,主要從事森林生態(tài)以及群落生態(tài)學(xué)研究。Email: wangyuchao2000@126.com

    王得祥(1966年生),男,教授,博士生導(dǎo)師,主要從事森林生態(tài)及森林健康可持續(xù)經(jīng)營(yíng)研究。Email: wangdx66@126.com

    2016-06-08

    猜你喜歡
    樣地坡度速效
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    昆明市主要綠化樹(shù)種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    不同殺螨劑對(duì)柑桔紅蜘蛛速效與持效性防治田間試驗(yàn)
    土壤速效鉀測(cè)定方法的影響因素
    綠色科技(2016年16期)2016-10-11 06:53:51
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    CT和MR對(duì)人上脛腓關(guān)節(jié)面坡度的比較研究
    施用不同有機(jī)物料對(duì)山地果園紅壤速效氮磷鉀含量的影響
    淫妇啪啪啪对白视频| 亚洲精品在线美女| 午夜激情欧美在线| 国产av不卡久久| 午夜两性在线视频| 国产亚洲精品av在线| 精品国产乱码久久久久久男人| 女同久久另类99精品国产91| 啦啦啦观看免费观看视频高清| 欧美国产日韩亚洲一区| 成人三级做爰电影| 嫁个100分男人电影在线观看| 90打野战视频偷拍视频| 97超视频在线观看视频| 精品国产三级普通话版| 国产三级中文精品| 人人妻,人人澡人人爽秒播| 免费在线观看影片大全网站| 精品福利观看| 少妇熟女aⅴ在线视频| 亚洲国产精品合色在线| 亚洲专区字幕在线| 国产亚洲精品综合一区在线观看| 亚洲中文日韩欧美视频| 国产亚洲精品一区二区www| 搡老熟女国产l中国老女人| 搡老熟女国产l中国老女人| 亚洲国产精品久久男人天堂| 日韩欧美在线二视频| 一卡2卡三卡四卡精品乱码亚洲| 欧美另类亚洲清纯唯美| 成人三级做爰电影| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产99精品国产亚洲性色| 一进一出抽搐动态| 一进一出抽搐动态| 日韩有码中文字幕| 久久久国产成人免费| 人妻久久中文字幕网| 亚洲五月天丁香| 亚洲七黄色美女视频| 一个人免费在线观看电影 | 成在线人永久免费视频| 搡老妇女老女人老熟妇| 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 日本一二三区视频观看| 国产精品,欧美在线| 国产激情久久老熟女| 免费看日本二区| 99久国产av精品| 久久久久久久久免费视频了| 亚洲天堂国产精品一区在线| e午夜精品久久久久久久| 亚洲精品456在线播放app | 国产一区二区三区视频了| 一级毛片高清免费大全| 精品无人区乱码1区二区| 中文字幕熟女人妻在线| 午夜福利18| 亚洲中文日韩欧美视频| 婷婷丁香在线五月| 性色av乱码一区二区三区2| 午夜影院日韩av| 中文字幕高清在线视频| 嫁个100分男人电影在线观看| 一进一出好大好爽视频| 亚洲成人免费电影在线观看| 人妻丰满熟妇av一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产成人aa在线观看| 久久久色成人| 国产亚洲av嫩草精品影院| 成人av在线播放网站| 99精品久久久久人妻精品| x7x7x7水蜜桃| 欧美色视频一区免费| 男插女下体视频免费在线播放| 好男人在线观看高清免费视频| 精品一区二区三区视频在线 | 99riav亚洲国产免费| 国产成人影院久久av| 亚洲国产色片| 99久久国产精品久久久| 99久久精品国产亚洲精品| 日韩欧美在线二视频| 两个人看的免费小视频| xxx96com| xxx96com| 一级毛片精品| 国产免费男女视频| 黄色日韩在线| 成人av一区二区三区在线看| 免费看a级黄色片| 国产蜜桃级精品一区二区三区| 亚洲av美国av| 综合色av麻豆| 91久久精品国产一区二区成人 | www日本黄色视频网| 在线永久观看黄色视频| 99在线人妻在线中文字幕| 国产精品女同一区二区软件 | 99riav亚洲国产免费| 午夜福利视频1000在线观看| 18禁黄网站禁片午夜丰满| 成熟少妇高潮喷水视频| 亚洲天堂国产精品一区在线| 亚洲欧美日韩高清在线视频| 99久久99久久久精品蜜桃| 国产乱人伦免费视频| 久久精品91蜜桃| 亚洲自拍偷在线| 在线永久观看黄色视频| 一个人免费在线观看的高清视频| 18禁美女被吸乳视频| 老鸭窝网址在线观看| 亚洲人成电影免费在线| 成年免费大片在线观看| 久久久久精品国产欧美久久久| 人人妻人人看人人澡| 日韩精品中文字幕看吧| 日韩免费av在线播放| 久久天堂一区二区三区四区| 在线观看午夜福利视频| 欧美xxxx黑人xx丫x性爽| 别揉我奶头~嗯~啊~动态视频| 国内精品久久久久久久电影| 免费看日本二区| 精品一区二区三区av网在线观看| 人人妻人人看人人澡| 超碰成人久久| 三级男女做爰猛烈吃奶摸视频| 999久久久国产精品视频| 欧美av亚洲av综合av国产av| 中文字幕人妻丝袜一区二区| 露出奶头的视频| 久久亚洲真实| 操出白浆在线播放| 国产69精品久久久久777片 | 巨乳人妻的诱惑在线观看| 亚洲激情在线av| 亚洲性夜色夜夜综合| 日本一二三区视频观看| 色视频www国产| 中文字幕最新亚洲高清| 91在线观看av| 一个人看视频在线观看www免费 | 国产亚洲精品av在线| 欧美日韩综合久久久久久 | 亚洲欧美日韩高清在线视频| 在线国产一区二区在线| 国产精品99久久99久久久不卡| 人妻夜夜爽99麻豆av| 村上凉子中文字幕在线| 国产精品,欧美在线| 99国产综合亚洲精品| 免费搜索国产男女视频| 1024手机看黄色片| 在线看三级毛片| 久久热在线av| 免费在线观看视频国产中文字幕亚洲| 亚洲中文字幕一区二区三区有码在线看 | 欧美激情久久久久久爽电影| 国产精品亚洲一级av第二区| 亚洲 国产 在线| 久久亚洲精品不卡| 天堂√8在线中文| 亚洲七黄色美女视频| 亚洲国产色片| 欧洲精品卡2卡3卡4卡5卡区| 小蜜桃在线观看免费完整版高清| 国产毛片a区久久久久| 日韩人妻高清精品专区| 精品国产三级普通话版| 午夜福利免费观看在线| 亚洲国产欧美一区二区综合| 精品一区二区三区视频在线 | 三级毛片av免费| 亚洲av电影不卡..在线观看| 波多野结衣巨乳人妻| 一个人看的www免费观看视频| 国产爱豆传媒在线观看| 久久婷婷人人爽人人干人人爱| 免费高清视频大片| а√天堂www在线а√下载| 欧美色视频一区免费| 亚洲国产精品久久男人天堂| 麻豆一二三区av精品| 精品不卡国产一区二区三区| 午夜福利18| 丰满人妻一区二区三区视频av | 好男人电影高清在线观看| 亚洲激情在线av| 少妇裸体淫交视频免费看高清| 国产精华一区二区三区| 成年女人毛片免费观看观看9| 日日夜夜操网爽| 变态另类成人亚洲欧美熟女| 久久欧美精品欧美久久欧美| 午夜免费观看网址| 日日干狠狠操夜夜爽| 在线免费观看不下载黄p国产 | 成人一区二区视频在线观看| 国产精品一区二区精品视频观看| 午夜福利18| 岛国视频午夜一区免费看| 岛国在线免费视频观看| 在线观看日韩欧美| 国产精品久久电影中文字幕| 成人精品一区二区免费| 色哟哟哟哟哟哟| 久久亚洲精品不卡| 在线观看美女被高潮喷水网站 | netflix在线观看网站| 一区二区三区激情视频| www日本黄色视频网| 国产精品av久久久久免费| a级毛片在线看网站| 99久久精品国产亚洲精品| 看片在线看免费视频| 久久香蕉精品热| 97人妻精品一区二区三区麻豆| 久久香蕉国产精品| 性色avwww在线观看| 精品久久蜜臀av无| 欧美乱妇无乱码| 性欧美人与动物交配| 免费一级毛片在线播放高清视频| 午夜两性在线视频| www日本黄色视频网| 精品久久久久久久毛片微露脸| 精品熟女少妇八av免费久了| 精品午夜福利视频在线观看一区| 真人做人爱边吃奶动态| 亚洲国产高清在线一区二区三| 搡老岳熟女国产| 一个人免费在线观看的高清视频| 久久久国产成人免费| 久久国产精品人妻蜜桃| 亚洲色图av天堂| 91av网一区二区| 久久久水蜜桃国产精品网| 国模一区二区三区四区视频 | 久久人人精品亚洲av| 亚洲国产欧美人成| 午夜精品在线福利| 午夜久久久久精精品| 亚洲午夜精品一区,二区,三区| 18禁黄网站禁片免费观看直播| 国产一区二区激情短视频| 国产美女午夜福利| 久久伊人香网站| 久久精品亚洲精品国产色婷小说| 亚洲成人久久性| 久久中文看片网| 亚洲黑人精品在线| 亚洲欧美精品综合久久99| 91九色精品人成在线观看| 成人国产一区最新在线观看| 国产亚洲欧美在线一区二区| 99在线视频只有这里精品首页| 最近最新中文字幕大全电影3| 色老头精品视频在线观看| 精品久久久久久,| 亚洲国产欧美人成| 一本一本综合久久| 两个人看的免费小视频| 一个人看的www免费观看视频| 国产不卡一卡二| 日本a在线网址| 狂野欧美激情性xxxx| 特级一级黄色大片| 久久国产乱子伦精品免费另类| or卡值多少钱| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看| 1000部很黄的大片| 国产av不卡久久| 狂野欧美白嫩少妇大欣赏| 一级毛片高清免费大全| 此物有八面人人有两片| 成人特级av手机在线观看| 欧美色视频一区免费| 国产久久久一区二区三区| 国语自产精品视频在线第100页| 最新中文字幕久久久久| 亚洲av熟女| 免费在线观看成人毛片| 国产精品人妻久久久久久| 91久久精品电影网| 激情 狠狠 欧美| 99热6这里只有精品| 亚洲欧美成人精品一区二区| 亚洲精品乱码久久久久久按摩| 亚洲av免费在线观看| 亚洲国产色片| 亚洲美女搞黄在线观看| 在线观看66精品国产| 国产精品av视频在线免费观看| 91久久精品国产一区二区成人| 欧美极品一区二区三区四区| 亚洲国产精品合色在线| 最新中文字幕久久久久| 丰满少妇做爰视频| 天堂影院成人在线观看| 三级国产精品欧美在线观看| АⅤ资源中文在线天堂| 成人性生交大片免费视频hd| 看片在线看免费视频| 人人妻人人澡欧美一区二区| 久久6这里有精品| 99久久人妻综合| 高清午夜精品一区二区三区| 国产成人精品婷婷| 三级国产精品片| 黄色配什么色好看| 搡女人真爽免费视频火全软件| 一级毛片我不卡| 日本三级黄在线观看| 日本免费在线观看一区| 欧美一区二区亚洲| 高清在线视频一区二区三区 | 精品久久国产蜜桃| 欧美xxxx性猛交bbbb| 欧美日韩国产亚洲二区| www日本黄色视频网| 欧美成人一区二区免费高清观看| 久久人妻av系列| 亚洲精品aⅴ在线观看| 国国产精品蜜臀av免费| 久久久国产成人精品二区| 色综合色国产| av免费在线看不卡| 在线免费十八禁| 日日摸夜夜添夜夜添av毛片| 亚洲精品成人久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲三级黄色毛片| 久久久午夜欧美精品| 搡老妇女老女人老熟妇| 国产在视频线在精品| 91久久精品电影网| 最近视频中文字幕2019在线8| 国产精品伦人一区二区| 夫妻性生交免费视频一级片| 国产成年人精品一区二区| 久久人人爽人人片av| 内射极品少妇av片p| 国产精品一及| 水蜜桃什么品种好| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 亚洲欧美一区二区三区国产| 尾随美女入室| 日日撸夜夜添| 国产精品美女特级片免费视频播放器| 久久鲁丝午夜福利片| 国产极品精品免费视频能看的| 日韩欧美精品v在线| 亚洲四区av| 欧美+日韩+精品| 日韩三级伦理在线观看| 国产69精品久久久久777片| 小说图片视频综合网站| 日本猛色少妇xxxxx猛交久久| 建设人人有责人人尽责人人享有的 | 久久久精品欧美日韩精品| a级毛色黄片| 欧美区成人在线视频| 亚洲18禁久久av| 欧美丝袜亚洲另类| 午夜视频国产福利| 亚洲国产精品sss在线观看| 日日摸夜夜添夜夜爱| 国产综合懂色| 婷婷色麻豆天堂久久 | 久久久久久久久中文| av免费观看日本| 午夜老司机福利剧场| 国产成人午夜福利电影在线观看| 亚洲美女视频黄频| 亚洲国产色片| 亚洲国产日韩欧美精品在线观看| 国产精品av视频在线免费观看| 2021少妇久久久久久久久久久| av专区在线播放| 村上凉子中文字幕在线| 国产在视频线精品| 中文资源天堂在线| 国产一区亚洲一区在线观看| 日韩一区二区三区影片| 久久久久国产网址| 国产亚洲91精品色在线| 久久久国产成人精品二区| 亚洲中文字幕一区二区三区有码在线看| 中国国产av一级| 99久久成人亚洲精品观看| 91午夜精品亚洲一区二区三区| 天天躁日日操中文字幕| 国产探花在线观看一区二区| 日日啪夜夜撸| 一区二区三区免费毛片| 亚洲人成网站在线观看播放| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 国产伦在线观看视频一区| 国产成人免费观看mmmm| 国产精品国产三级国产专区5o | 18禁在线播放成人免费| 亚洲国产精品国产精品| 日本免费一区二区三区高清不卡| 视频中文字幕在线观看| 亚洲在线自拍视频| 欧美又色又爽又黄视频| 午夜福利在线观看吧| 成人午夜精彩视频在线观看| 最新中文字幕久久久久| 熟女人妻精品中文字幕| 三级经典国产精品| 国产男人的电影天堂91| 日韩高清综合在线| 最近最新中文字幕大全电影3| 亚洲色图av天堂| 亚洲成人中文字幕在线播放| 色综合站精品国产| 乱系列少妇在线播放| 亚洲三级黄色毛片| 成人毛片60女人毛片免费| 七月丁香在线播放| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 精品少妇黑人巨大在线播放 | 国产亚洲最大av| 人妻系列 视频| 国产成人一区二区在线| 久久99蜜桃精品久久| 国产精品三级大全| 精品欧美国产一区二区三| 久久久a久久爽久久v久久| 久久精品91蜜桃| 日韩大片免费观看网站 | 欧美色视频一区免费| 精品久久久噜噜| 国产免费一级a男人的天堂| 欧美精品国产亚洲| a级毛色黄片| 一级黄色大片毛片| 国产精品女同一区二区软件| 国产综合懂色| 成人三级黄色视频| 日韩 亚洲 欧美在线| 日本黄色视频三级网站网址| 一个人免费在线观看电影| 亚洲三级黄色毛片| 久久这里只有精品中国| 99热6这里只有精品| 日本免费在线观看一区| 91在线精品国自产拍蜜月| 99热精品在线国产| 久久久久久久久大av| 国产午夜精品一二区理论片| 欧美又色又爽又黄视频| 免费搜索国产男女视频| 国产伦理片在线播放av一区| av免费在线看不卡| 欧美极品一区二区三区四区| 亚洲av中文av极速乱| 日韩亚洲欧美综合| 国产黄片美女视频| 国产黄片视频在线免费观看| 日日摸夜夜添夜夜添av毛片| 亚洲国产欧洲综合997久久,| 久久草成人影院| 只有这里有精品99| 国产美女午夜福利| 成年女人看的毛片在线观看| 欧美成人a在线观看| 中文资源天堂在线| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久久丰满| 日本午夜av视频| 国产黄色视频一区二区在线观看 | 91久久精品国产一区二区成人| 欧美日韩国产亚洲二区| 床上黄色一级片| 色综合站精品国产| 直男gayav资源| 欧美bdsm另类| 免费看av在线观看网站| 午夜福利高清视频| 国产精华一区二区三区| 亚洲成人av在线免费| 18禁在线播放成人免费| 国产在线男女| 精品久久久久久久久av| 天堂影院成人在线观看| 国内揄拍国产精品人妻在线| 久久久国产成人精品二区| 2022亚洲国产成人精品| 国产精品美女特级片免费视频播放器| 全区人妻精品视频| 国产精品久久视频播放| 精品少妇黑人巨大在线播放 | 久久精品熟女亚洲av麻豆精品 | 欧美精品一区二区大全| 99久国产av精品国产电影| 三级男女做爰猛烈吃奶摸视频| 超碰av人人做人人爽久久| 日本一本二区三区精品| 国产精品无大码| 一个人免费在线观看电影| 欧美成人免费av一区二区三区| 国产成人freesex在线| 亚洲伊人久久精品综合 | 午夜免费激情av| 天天一区二区日本电影三级| 丰满人妻一区二区三区视频av| 免费在线观看成人毛片| 国产午夜福利久久久久久| 美女xxoo啪啪120秒动态图| 中文字幕人妻熟人妻熟丝袜美| 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| 欧美+日韩+精品| 91久久精品国产一区二区三区| 国产激情偷乱视频一区二区| 日韩三级伦理在线观看| 亚洲国产精品久久男人天堂| 大话2 男鬼变身卡| 美女xxoo啪啪120秒动态图| 欧美高清性xxxxhd video| 一区二区三区高清视频在线| 国产精品女同一区二区软件| 国产av一区在线观看免费| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 国产综合懂色| 日韩欧美精品v在线| 男人的好看免费观看在线视频| 在线播放无遮挡| 色吧在线观看| 午夜激情欧美在线| 在线观看66精品国产| 国产男人的电影天堂91| 又爽又黄无遮挡网站| 欧美三级亚洲精品| 亚洲av一区综合| 午夜久久久久精精品| 免费看a级黄色片| 2021少妇久久久久久久久久久| 91狼人影院| 免费电影在线观看免费观看| 国产精品国产三级国产专区5o | 国产av不卡久久| 久久久a久久爽久久v久久| 黑人高潮一二区| 一级毛片aaaaaa免费看小| 色哟哟·www| 亚洲av中文字字幕乱码综合| 亚洲国产欧美人成| 少妇的逼水好多| 色综合站精品国产| 国产又色又爽无遮挡免| 中国国产av一级| 午夜视频国产福利| 综合色av麻豆| 欧美高清成人免费视频www| 亚洲av电影在线观看一区二区三区 | 久久精品人妻少妇| 午夜激情欧美在线| 2022亚洲国产成人精品| 色视频www国产| 3wmmmm亚洲av在线观看| 亚洲久久久久久中文字幕| www.色视频.com| 乱码一卡2卡4卡精品| 欧美精品国产亚洲| 成人一区二区视频在线观看| 国国产精品蜜臀av免费| 亚洲成av人片在线播放无| 搡老妇女老女人老熟妇| 久久久久国产网址| 一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 亚洲欧美日韩高清专用| 中文字幕av成人在线电影| 三级国产精品片| 日本色播在线视频| 国产在线一区二区三区精 | 69av精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 一级毛片我不卡| 国产免费一级a男人的天堂| 午夜爱爱视频在线播放| 又粗又爽又猛毛片免费看| 国产精品一二三区在线看| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| 亚洲综合色惰| 亚洲av成人精品一区久久| 亚洲国产欧美人成| 久久久久久久久久成人| 99国产精品一区二区蜜桃av| .国产精品久久| 欧美色视频一区免费| 久久久午夜欧美精品| 国产69精品久久久久777片| eeuss影院久久| 国产大屁股一区二区在线视频| 国产亚洲91精品色在线| 午夜视频国产福利| 久久99热这里只频精品6学生 | 国产精品久久久久久精品电影| 精品人妻一区二区三区麻豆| 少妇的逼好多水|