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

    半干旱區(qū)濕地干草原交錯帶判定及植物群落與環(huán)境關系的典范對應分析

    2016-02-09 12:31:11李海濤劉小丹張克斌西拉杰馬默
    中國水土保持科學 2016年6期
    關鍵詞:排序環(huán)境

    李海濤,劉小丹,張克斌?,西拉杰·馬默

    (1.北京林業(yè)大學水土保持學院,100083,北京;2.河北環(huán)境工程學院,066102,河北秦皇島)

    半干旱區(qū)濕地干草原交錯帶判定及植物群落與環(huán)境關系的典范對應分析

    李海濤1,劉小丹2,張克斌1?,西拉杰·馬默1

    (1.北京林業(yè)大學水土保持學院,100083,北京;2.河北環(huán)境工程學院,066102,河北秦皇島)

    為探究半干旱區(qū)濕地干草原生態(tài)交錯帶的邊界及影響其植被空間分布格局的主要環(huán)境因子,以寧夏哈巴湖荒漠濕地自然保護區(qū)內(nèi)的四兒灘濕地為例,采用樣線調(diào)查法獲得植被與環(huán)境數(shù)據(jù):首先利用游動分割窗技術對研究區(qū)進行交錯帶定量判定,將其分為旱生帶、交錯帶、濕生帶并計算沿不同生境梯度β多樣性變化趨勢;然后,采用典范對應分析法對研究區(qū)植被數(shù)據(jù)矩陣和環(huán)境數(shù)據(jù)矩陣進行排序,同時進行排序軸與土壤環(huán)境因子間的秩相關分析。結果表明:采用游動分割窗技術判定四兒灘濕地干草原交錯帶最適窗口寬度為6~8個樣方;不同生境帶β多樣性指數(shù)存在明顯梯度變化;典范對應分析結果顯示,植被在CCA排序圖上的分布與游動分割窗技術的劃分結果相吻合,含鹽量(質(zhì)量分數(shù))、含水量是影響四兒灘植被空間分布的主要環(huán)境因子,對含鹽量和含水量有較高適應值的物種集中分布于土壤鹽漬化嚴重的濕生帶。研究結果將為四兒灘濕地的恢復和保護提供參考。

    四兒灘濕地;生態(tài)交錯帶;游動分割窗技術;典范對應分析

    生態(tài)交錯帶是一個氣候、環(huán)境敏感地帶,對其周圍環(huán)境變化具有提前指示作用。群落交錯帶的定量判定是研究其生態(tài)學過程的基礎,特別是在斑塊尺度格局和動態(tài)、群落尺度和種群統(tǒng)計等研究中具有重要意義[1]。植被是研究交錯帶的主要載體,因此研究生態(tài)交錯帶,需探討影響植被分布的主要環(huán)境因子。國內(nèi)外學者分別對生態(tài)交錯帶和植被分布與環(huán)境關系的對應分析進行了大量研究,但將兩者結合分析的研究極少。于大炮等[2]利用游動分割窗技術對長白山高山林線交錯區(qū)進行了定量判定,并能準確地反映交錯帶植被的變異。郝占慶等[34]通過樣地調(diào)查結合典范對應分析探討了長白山北坡草本和木本植物分布與環(huán)境因子的對應關系。賈寶全等[5]通過植被蓋度隨著與綠洲距離的變化發(fā)現(xiàn)植被蓋度的變化為綠洲附近與交錯帶外圍靠近流沙的地方較低,而位于流沙和綠洲之間的部分相應較高。朱軍濤等[6]研究額濟納荒漠綠洲的植物群落分布狀況與地下水環(huán)境的對應關系,得出地下水環(huán)境的礦化度、pH、水埋深等因子對植被群落產(chǎn)生顯著影響。C.J.Burk[7]發(fā)現(xiàn)新英格蘭內(nèi)陸與淡水沼澤交錯區(qū)的維管植物多樣性呈現(xiàn)濕地與水域交錯區(qū)向沼澤中部逐漸增大的趨勢。

    我國西部內(nèi)陸濕地主要分布于新疆、甘肅和寧夏等地,占我國濕地總面積的10%左右,而且西北地區(qū)濕地以鹽湖和沼澤濕地為主,主要分布于江河源頭地區(qū)、綠洲和內(nèi)陸湖濱等生態(tài)環(huán)境脆弱地帶,一經(jīng)破壞則很難恢復[8]。四兒灘濕地作為寧夏哈巴湖荒漠濕地保護區(qū)內(nèi)最大的一塊濕地,對其連續(xù)11年的采樣調(diào)查發(fā)現(xiàn),濕地面積不斷萎縮,濕地功能逐漸退化。本文擬通過游動分割窗技術首先對其進行生境帶的劃分,采用典范對應分析探討影響研究區(qū)鹽堿濕地植被分布的主要環(huán)境因子,分析其萎縮、退化原因,以期為哈巴湖荒漠濕地的保護和恢復提供理論依據(jù),同時為保護干旱、半干旱區(qū)濕地和維持其生態(tài)系統(tǒng)穩(wěn)定性提供參考。

    1 研究區(qū)概況

    鹽池縣位于寧夏回族自治區(qū),南靠黃土高原,北與毛烏素沙地相連,地形上自南向北是從黃土丘陵區(qū)向鄂爾多斯臺地過渡地帶。在氣候上,屬典型的大陸性季風氣候,降雨少,年降水量346.9 mm;日照時長達2804.9 h,年均無霜期190 d,年均氣溫9.3℃;植被屬于干草原向荒漠的過渡地帶,植被類型有灌叢、草原、草甸、沼澤、荒漠等植被類型;土壤類型主要為灰鈣土。

    四兒灘濕地位于鹽池縣城以南8 km處,E107°24′, N37°43′,面積大約為1 420 hm2,海拔為1 320~1 370 m,隸屬于哈巴湖荒漠濕地國家級自然保護區(qū)。濕地水域面積主要受區(qū)域降水量影響,季節(jié)性積水特征明顯,土壤鹽漬化嚴重,屬于我國干旱、半干旱區(qū)的典型內(nèi)陸濕地。

    2 研究方法

    2.1 調(diào)查與采樣

    鹽池四兒灘濕地是國家荒漠化定位監(jiān)測項目的濕地定位監(jiān)測點。筆者于2015年7月(植物生長季)對四兒灘濕地進行野外植被調(diào)查。樣線布設從濕地核心區(qū)固定的水泥樁開始向東、西、南、北4個方向輻射4條樣線,考慮到調(diào)查區(qū)域都是草本植被群落,沿樣線每隔20 m隨機布設1 m×1 m的小樣方,直至旱生帶植被群落結構不發(fā)生變化時樣方停止。每個樣方調(diào)查內(nèi)容包括植物名稱、株數(shù)、高度、蓋度、生物量等指標。由于樣方間隔距離較小,每隔1個樣方取土樣,取土深度為0~30 cm,土樣測定內(nèi)容包括速效磷、速效鉀、全氮、有機質(zhì)、含水量和電導率。有研究指出土壤含鹽量(質(zhì)量分數(shù))與電導率呈顯著正相關,所以直接采用電導率表示土壤鹽漬化程度是可行的[9]。

    2.2 生態(tài)交錯帶判定

    生態(tài)交錯帶的判定基于研究區(qū)的植被重要值,利用多個(偶數(shù))樣方組成的窗體沿樣線逐個樣方滑動,采用平方歐氏距離(Dn)計算2個1/2窗體A和B間的相異系數(shù)[10]。

    式中:Dn為窗體為n時的平方歐氏距離,XiA和XiB分別對應參數(shù)為i時A、B半窗體的值,m為窗口變量數(shù)。

    2.3 β多樣性指標

    β多樣性也被稱為生境間的多樣性,強調(diào)的是沿環(huán)境梯度不同生境群落之間物種組成的相異性,β多樣性越大,表明不同群落之間的共有種越少[11]。計算四兒灘濕地沿濕地中心輻射出的4條樣線β多樣性指標變化可以進一步檢驗游動分割窗技術的準確性。

    二元屬性數(shù)據(jù)測度包括:

    數(shù)量數(shù)據(jù)測度法為Bray-curtis指數(shù)

    式中:g(H)、l(H)分別為沿生境梯度H增加和失去的物種數(shù)。a、b為2群落的物種數(shù),c為2群落共有的物種數(shù),j為2群落共有種中重要值較小者之和。

    2.4 數(shù)據(jù)處理

    采用游動分割窗技術對樣地進行生態(tài)交錯帶的定量判定,根據(jù)所得結果對研究區(qū)進行濕生帶、交錯帶、旱生帶3帶劃分;然后計算出4條樣線從濕地中心沿3個生境梯度β多樣性指標的變化(表1);最后在每條樣線上的不同生境帶選擇3個典型樣方,以調(diào)查的32種草本物種重要值組成36×32的植被數(shù)據(jù)矩陣,以測得的土壤理化性質(zhì)組成36×6的環(huán)境數(shù)據(jù)矩陣,利用Canoco for windows 4.5進行典范對應分析。

    3 結果與分析

    3.1 適宜窗口寬度的選擇

    圖1是以四兒灘濕地南樣線上的植物重要值為指標,分別采用4、6、8、10個樣方窗口寬度計算的平方歐氏距離在樣帶系列上的分布圖。由圖1可知,當窗口寬度為4個樣方時,平方歐氏距離出現(xiàn)多波峰噪聲干擾,并且峰值陡上陡下,說明以4個樣方作為窗口寬度太小不能準確判斷交錯帶。當窗口寬度為6~10個樣方時,平方歐式距離變化趨于一致,但窗口寬度為6個樣方時,還是存在2半窗口之間過度不平穩(wěn)的現(xiàn)象。窗口寬度為10個樣方時,平方歐式距離變化不明顯,并且削弱了最大峰值。綜合比較可知,南樣線選擇8個樣方作為窗口寬度最為合理。3.2 交錯帶判定結果分析

    圖1 南樣線不同窗口寬度平方歐式距離(Dn)與樣帶系列Fig.1 Squared euclidean distance(Dn)under different window width of south transect line

    根據(jù)上文對于南樣線適宜窗口寬度的選擇方法判定其余3條樣線,可知東樣線的適宜窗口寬度為6個樣方,西樣線和北樣線依舊為8個樣方。圖2為四兒灘東、西、北3條樣線分別在適宜窗口寬度下,平方歐氏距離在樣帶系列上的分布圖。由游動分割窗技術的峰值理論可以明顯判斷出各樣線交錯帶位置和寬度[12]。由圖2可知各樣線峰值較高說明交錯帶過渡明顯,不同樣線的生態(tài)交錯帶寬度不同。結合野外調(diào)查可知,以東樣線為例:濕生帶(樣方1~6)只出現(xiàn)鹽爪爪(Kalidium foliatum)、堿蓬(Suaeda glauca)等濕生植物,交錯帶(樣方7~11)不只出現(xiàn)芨芨草(Achnatherum splendens)等濕生植物、還有草瑞香(Diarthron inifalium)等交錯帶特有種及蒙古蟲實(Corispermum mongolicum)等草原植物。旱生帶(樣方12~21)主要出現(xiàn)了刺沙蓬(Salsola ruthenica)、豬毛蒿(Artemisia scoparia)等耐干旱的植物。這也說明游動分割窗技術能夠準確判定交錯帶的位置。

    3.3 不同生境梯度植物群落β多樣性變化

    由表1可以看出4條樣線Jaccard指數(shù)、Sorensen指數(shù)和Bray-curtis 指數(shù)都是交錯帶旱生帶>濕生帶交錯帶的生境梯度變化,而Cody指數(shù)則是濕生帶交錯帶>交錯帶旱生帶,說明四兒灘濕地濕生帶與交錯帶共有種數(shù)少,交錯帶與旱生帶共有種數(shù)相對較多更替速率大。不同樣線間的相同生境梯度的β多樣性指數(shù)差異不大,從濕地中心向外輻射的4條樣線都表現(xiàn)出一致的規(guī)律,表明四兒灘濕地確實存在明顯的生境帶梯度。

    圖2 東樣線6個樣方、西和北樣線8個樣方為窗口寬度的平方歐式距離與樣帶系列Fig.2 Squared euclidean distance(Dn)under 6 or 8 quadrats as window width at east or west and north transect line

    3.4 四兒灘濕地植被典范對應分析

    為探究前文得出的生境帶分化原因,對四兒灘濕地植被與環(huán)境因子進行典范對應分析。CCA二維排序軸是一種通過特征根值反映相應的排序軸所包涵的物種重要值矩陣中的信息量大小的生態(tài)梯度抽象軸[13],相關系數(shù)表示抽象軸與真實環(huán)境梯度的相關性的大小[14]。由表2可得,前兩軸特征根值之和占總特征根值的68.8%,而排序軸1,物種與環(huán)境的相關系數(shù)高達0.975,累計貢獻率為50.9,反映了排序的絕大部分信息,排序結果可靠。

    表1 2015年四兒灘濕地不同生境梯度植物群落β多樣性Tab.1 Plant community diversity of different habitat gradient in Siertan Wetland

    表2 典范對應排序軸特征值及解釋比例Tab.2 Eigenvalues and explanation percentage of CCA axes

    由圖3可知,電導率和含水量箭頭連線最長且都與物種第1排序軸和環(huán)境第1排序軸呈極顯著正相關,相關系數(shù)分別為0.868 4、0.917 0和0.841 6、0.891 4(表3),說明它們在第1排序軸上的位置很大程度上反映了其生態(tài)特點,是影響四兒灘濕地植被分布的主要環(huán)境因子。相比較而言速效磷、速效鉀、有機質(zhì)箭頭連線較短與第1排序軸正相關性較小,對植被分布影響較小,而N元素雖箭頭連線較長,但與第1排序軸呈不顯著負相關。從圖3還可以看出電導率和含水量與第1軸方向基本一致,即表明第1軸從左到右,土壤含水量和含鹽量逐漸增大。沿第1排序軸從右到左,將圖3各物種垂直投影于環(huán)境變量箭頭線(或延長線)上,根據(jù)物種距離變量箭頭的相對位置,可得出種1(鹽爪爪)、種2 (堿蓬)和種5(芨芨草)在含鹽量和含水量較大時具有最適值,分布于鹽漬化程度較高的濕生帶;而種23(苦豆子)、24(駱駝蓬)、25(老瓜頭)等物種投影于兩環(huán)境變量延長線上,最適值較低分布于旱生帶,種7(白刺),6(草瑞香)最適值適中,主要分布于交錯帶。對與前文通過游動分割窗技術的3帶劃分及野外調(diào)查的植被分布結果相吻合。

    表3 物種、環(huán)境因子前兩個CCA排序軸及環(huán)境因子間的相關系數(shù)Tab.3 Correlation coefficients for species axis 1 and axis 2,environmental factors axis 1 and axis 2 and environmental factors

    圖3 四兒灘濕地植被與環(huán)境因子的CCA二維排序圖Fig.3 CCA two-dimensional ordination diagram of vegetation and environmental factors in Siertan Wetland

    4 結論與討論

    游動分割窗技術主要通過窗體滑動求得相鄰取樣區(qū)間的滑動平均相異系數(shù),消除因個別樣方取樣偏差或環(huán)境突變造成的單個兩兩樣方間相異系數(shù)的陡然變化引起的交錯帶誤判[15]。采用游動分割窗技術過程中,必須首先探討最適窗口寬度:窗口選擇過小,滑動平均就沒有意義并產(chǎn)生多波峰干擾判斷;若窗口寬度選擇太大,容易削弱波峰。對四兒灘濕地草本重要值計算分析發(fā)現(xiàn)窗口寬度選擇為6~8個窗口寬度時交錯帶邊界清晰可見。然而田曉玲等[16]在研究晉西林草復合系統(tǒng)交錯區(qū)中認為8~10個窗口寬度為宜;石培禮等[15]在四川岷江冷杉林線交錯帶的研究中得出最適窗口寬度應為8~12。綜合分析可知由于樣地尺度和研究對象的不同可能導致最適窗口寬度的差異,如研究海拔梯度變化引起的高山林線交錯帶和復合農(nóng)林系統(tǒng)時,由于布設樣線跨度大,樣方數(shù)量多,增大窗口寬度可以消減多波峰干擾。四兒灘濕地4條樣線沿不同生境梯度上β多樣性表現(xiàn)出較大差異,出現(xiàn)明顯物種更替,體現(xiàn)在:濕生帶明顯集中分布耐濕、耐鹽堿的植被;交錯帶植被種類多,生物多樣性豐富,主要以以濕生、鹽中生、中生的植被為主;旱生帶以中生、旱中生的草原植被為主。

    四兒灘濕地從中心向外輻射范圍不足1 km,屬于小規(guī)模斑塊狀濕地,各環(huán)境因子空間變異值不大。不同于郝占慶等[3]對長白山北坡草本植物分布與環(huán)境關系的研究,區(qū)域跨度大,坡度、海拔等環(huán)境因子大幅變化導致土壤因子的變化顯著、植被分布空間差異明顯;但是,典范對應分析法在對應分析的迭代過程中結合了多元回歸分析,將每次得到的排序值均與環(huán)境因子進行多元線性回歸,使環(huán)境因子對排序結果的影響更為直觀準確[17],所以同樣適用于小尺度的生態(tài)交錯區(qū)。四兒灘濕地植被空間分布格局主要源于地表水存在有邊緣向中心遞增的水文梯度,受地勢影響土壤鹽分也主要集中在濕地中心。CCA結果也表明相比較其他營養(yǎng)元素,含水量和含鹽量是影響四兒灘植被分布的主要因素,即對水分和鹽度的耐受力決定了植被的空間分布格局。

    由于干旱、半干旱地區(qū)濕地受土壤沙化和土地鹽漬化的影響,使得濕地植被分布呈斑塊狀和片狀,濕地規(guī)模也受到限制,對環(huán)境響應更加敏感[18]。不同于其他大規(guī)模典型濕地,西北內(nèi)陸濕地受降雨量影響較大,季節(jié)積水現(xiàn)象明顯,鹽堿化程度較高,容易產(chǎn)生逆向演替逐漸退化。鑒于西北地區(qū)濕地多分布于生態(tài)環(huán)境敏感地帶的特點及生態(tài)交錯帶的提前指示作用,對四兒灘濕地生態(tài)交錯帶的判定和影響其植被分布的研究結論對于荒漠濕地的保護意義重大。今后保護區(qū)的工作重點應是保護濕地不受人為破壞、補充濕地入水量、改善鹽漬化土地,維護濕地植物物種多樣性,使?jié)竦刂脖豢臻g分布格局朝正向演替。

    [1] 王曉,張克斌,楊曉暉,等.半干旱區(qū)濕地干草原交錯帶邊界判定及其變化[J].生態(tài)學報,2012,32 (16):5122.Wang Xiao,Zhang Kebin,Yang Xiaohui,et al.Research on boundary definition and changes of wetland-dry grassland[J].Acta Ecologica Sinica,2012,32(16):5122. (in Chinese)

    [2] 于大炮,唐立娜,王紹先,等.長白山北坡植被垂直帶群落交錯區(qū)的定量判定[J].應用生態(tài)學報,2004,15 (10):1763. Yu Dapao,Tang Lina,Wang Shaoxian,et al.Quantitative methodologies for ecotone determination on north slope of Changbai Mountains[J].Chinese Journal of Applied Ecology,2004,15(10):1763.(in Chinese)

    [3] 郝占慶,郭水良.長白山北坡草本植物分布與環(huán)境關系的典范對應分析[J].生態(tài)學報,2003,23(10): 2002. Hao Zhanqing,Guo Shuiliang.Canonical correspondence analysis on relationship of herbs with their environments on northern slope of Changbai mountain[J].Acta Ecologica Sinica,2003,23(10):2002.(in Chinese)

    [4] 郝占慶,郭水良,葉吉.長白山北坡木本植物分布與環(huán)境關系的典范對應分析[J].植物生態(tài)學報,2003, 27(6):737. Hao Zhanqing,Guo Shuiliang,Ye Ji.Canonical correspondence analysis on relationship of woody plants with their environments on the northern slope of Changbai mountain[J].Chinese Journal of Applied Ecology, 2003,27(6):737.(in Chinese)

    [5] 賈寶全,慈龍駿,蔡體久,等.綠洲荒漠交錯帶環(huán)境特征初步研究[J].應用生態(tài)學報,2002,13(9):1104. Jia Baoquan,Ci Longjun,Cai Tijiu,et al.Preliminary research on environmental characteristics of oasis-desert ecotone[J].Chinese Journal of Applied Ecology,2002, 13(9):1104.(in Chinese)

    [6] 朱軍濤,于靜潔,王平,等.額濟納荒漠綠洲植物群落的數(shù)量分類及其與地下水環(huán)境的關系分析[J].植物生態(tài)學報,2011,35(5):480. Zhu Juntao,Yu Jingjie,Wang Ping,et al.Quantitative classification and analysis of relationships between plant communities and their groundwater environment in the Ejin Desert Oasis of China[J].Chinese Journal of Applied Ecology,2011,35(5):480.(in Chinese)

    [7] Burk C J.A four-year analysis of vegetation following an oil spill in a freshwater marsh[J].Journal of Applied E-cology,1977,323(1):117.

    [8] 朱剛,高會軍,曾光.西北內(nèi)陸干旱區(qū)濕地景觀格局變化:以新疆準噶爾內(nèi)流區(qū)為例[J].干旱區(qū)資源與環(huán)境,2014,28(8):77. Zhu Gang,Gao Huijun,Zeng Guang.Changes of wetland landscape pattern of the arid inland in northwest of Chinaa case of the inland water of Junggar in Xinjiang[J]. Journal of Arid Land Resources and Environment,2014, 28(8):77.(in Chinese)

    [9] 張瑜斌,鄧愛英,莊鐵誠,等.潮間帶土壤鹽度與電導率的關系[J].生態(tài)環(huán)境,2003,12(2):165. Zhang Yubin,Deng Aiying,Zhuang Tiecheng,et al.Relation between soil salinity in intertidal zone and electric conductivity[J].Ecology and Environment,2003,12 (2):165.(in Chinese)

    [10]石培禮,李文華.生態(tài)交錯帶的定量判定[J].生態(tài)學報,2002,22(4):589. Shi Peili,Li Wenhua.Quantitative methodologies for ecotone determination[J].Acta Ecologica Sinica,2002,22 (4):589.(in Chinese)

    [11]彭思羿,胡廣,于明堅.千島湖島嶼維管植物β多樣性及其影響因素[J].生態(tài)學報,2014,34(14):3867. Peng Siyi,Hu Guang,Yu Mingjian.Beta diversity of vascular plant and its influencing factors on islands in the Thousand Island Lake[J].Acta Ecologica Sinica,2014, 34(14):3867.(in Chinese)

    [12]時項鋒,許冬梅,邱開陽,等.游動分割窗技術在景觀界面影響域判定中的應用:以牛枝子黑沙蒿群落界面為例[J].草業(yè)科學,2010,27(4):31. Shi Xiangfeng,Xu Dongmei,Qiu Kaiyang,et al.Application of moving split-window technique in quantitative methodologies for edge influence of landscape boundary:a case study of Lespedeza potaninii-Artemisia ordosica community[J].Pratacultural Science,2010,27(4):31. (in Chinese)

    [13]余敏,周志勇,康峰峰,等.山西靈空山小蛇溝林下草本層植物群落梯度分析及環(huán)境解釋[J].植物生態(tài)學報,2013,37(5):378. Yu Min,Zhou Zhiyong,Kang Fengfeng,et al.Gradient analysis and environmental interpretation of understory herb-layer communities in Xiaoshegou of Lingkong Mountain,Shanxi,China[J].Chinese Journal of Plant Ecology,2013,37(5):378.(in Chinese)

    [14]張玲,王承義.大小興安嶺過渡區(qū)闊葉紅松林次生演替階段群落多樣性指數(shù)與環(huán)境因子關系[J].森林工程,2014,30(5):3. Zhang Ling,Wang Chengyi.Relationship between species diversity and environmental factors of plant community at different succession stages in ecotone of Great Xing' an and Xiaoxing'an mountains[J].Forest Engineering, 2014,30(5):3.(in Chinese)

    [15]石培禮,劉興良.游動分割窗技術在生態(tài)交錯帶定量判定中的應用:以四川巴郎山岷江冷杉林線為例[J].植物生態(tài)學報,2002,26(2):190.Shi Peili,Liu Xingliang.The application of moving splitwindow technique in determining ecotone:a case study of abies faxoniana timberline in Balang mountain in Sichuan province[J].Acta Phytoecologica Sinica,2002,26(2): 190.(in Chinese)

    [16]田曉玲,畢華興,云雷,等.晉西黃土區(qū)林草復合系統(tǒng)草本植物多樣性特征[J].北京林業(yè)大學學報,2011, 33(1):66. Tian Xiaoling,Bi Huaxing,Yun Lei,et al.Herbaceous plant diversity of silvopastoral system in the loess region of western Shanxi province,northern China[J].Journal of Beijing Forestry University,2011,33(1):66.(in Chinese)

    [17]楊思存,逄煥成,王成寶,等.基于典范對應分析的甘肅引黃灌區(qū)土壤鹽漬化特征研究[J].中國農(nóng)業(yè)科學, 2014,47(1):102.

    Yang Sicun,Pang Huancheng,Wang Chengbao,et al. Characterization of soil salinization based on canonical correspondence analysis method in Gansu Yellow River irrigation district of northwest China[J].Scientia Agricultura Sinica,2014,47(1):102.(in Chinese)

    [18]李瑞,張克斌,劉云芳,等.西北半干旱區(qū)濕地生態(tài)系統(tǒng)植物群落空間分布特征研究[J].北京林業(yè)大學學報,2008,30(1):11. Li Rui,Zhang Kebin,Liu Yunfang,et al.Plant community spatial distribution pattern of wetland ecosystem in a semi-arid region of northwestern China[J].Journal of Beijing Forestry University,2008,30(1):11.(in Chinese)

    Ecotone definition between wetland and dry grassland in semi-arid region and the canonical correspondence analysis of plant community in relation to environment factors

    Li Haitao1,Liu Xiaodan2,Zhang Kebin1,Siraj Mammo1

    (1.College of Soil and Water Conservation,Beijing Forestry University,100083,Beijing,China; 2.Hebei University of Environmental Engineering,066102,Qinhuangdao,Hebei,China)

    [Background]Siertan Wetland is the largest wetland of Habahu Desert Wetland Nature Reserve in Yanchi County,Ningxia Autonomous Region.Through the annual vegetation investigation which has been carried out for 11 years,we found that the area of Siertan Wetland is shrinking year by year.Many researchers concluded that the ecotone is very sensitive to the change of environment,thus it can be used to predict the change of the ecological environment in advance.This study is to quantitatively determine the ecological ecotone boundary and the main environmental factors that affect the spatial distribution pattern of the species.[Methods]Firstly we used moving split-window technique to quantitatively determine the ecological ecotone.Use of the moving split-window technique required that:we chose an appropriate even number of representative samples as a whole window,then divided the whole window in half and calculated their dissimilarity coefficient.We continuously moved the window right to the next sampling point and repeated the operation above again,until the right window reached the final sampling points.We then draw a scatter diagram of the dissimilarity coefficient changes by quadrats.Further,according to the result of determination we divided the research area into wetland regions,ecotone and arid region,and analyzed characteristics of beta diversity in the three different habitats of the four different sample lines.Finally,Canonical Correspondence Analysis(CCA)was used to determine which environmental factors contributed most to the distribution of the vegetation.[Results] The results showed that ecotone was determined clearly by moving split-window technique under 6-8 quadrats,the optimal window width of east transect line was 6 quadrates,the optimal window width of west,south and north transect line was 8 quadrates.Beta diversity index showed significant differences in different habitats.The result of CCA indicated that the total salt content and moisture content were the main environmental factors that affected the vegetation distribution in Siertan Wetland.The correlation of the two environmental factors and the first species axis were respectively 0.868 4 and 0.917 0,the correlation of the two environmental factors and the first environment axis were respectively 0.841 6 and 0.891 4.[Conclusions]According to the dimension of investigatory scale and the visibility of ecotone transitivity,appropriate window width can be chosen especially in the selection of appropriate smallest window width.Species with higher fitness for total salt content and moisture content variability are mainly distributed in wet region with serious soil salinization.The results of the study will provide valuable reference for Siertan Wetland's restoration and protection.The future work of Habahu Nature Reserve administration should be focused on protecting Siertan Wetland from human being destruction,ensuring adequate water supply,improving the salinization land,and maintaining plant species diversity in wetland.

    Siertan Wetland;ecotone;moving split-window technique;canonical correspondence analysis

    Q948

    A

    1672-3007(2016)06-0070-08

    10.16843/j.sswc.2016.06.009

    2016 06 27

    2016 09 19

    項目名稱:“十三五”國家重點研發(fā)計劃“生態(tài)畜牧業(yè)與沙化草地治理技術集成及產(chǎn)業(yè)示范”(2016YFC0500908);國家自然科學基金“半荒漠濕地植被群落時空分布特征研究”(30771764);國家林業(yè)局荒漠化定位監(jiān)測項目“鹽池荒漠化定位監(jiān)測”(660550)

    李海濤(1991—)男,碩士研究生。主要研究方向:水土保持與荒漠化防治。Email:lhtloyal@bjfu.edu.cn

    ?通信作者簡介:張克斌(1957—),男,教授,博士生導師。主要研究方向:荒漠化防治與監(jiān)測及干旱區(qū)環(huán)境管理。Email:ctccd@126.com

    猜你喜歡
    排序環(huán)境
    排排序
    排序不等式
    長期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學習的虛擬仿真環(huán)境
    恐怖排序
    孕期遠離容易致畸的環(huán)境
    不能改變環(huán)境,那就改變心境
    節(jié)日排序
    環(huán)境
    孕期遠離容易致畸的環(huán)境
    晚上一个人看的免费电影| 一区在线观看完整版| 最近最新中文字幕免费大全7| 9色porny在线观看| 久久影院123| 久久精品国产综合久久久| 欧美成人午夜精品| 欧美黄色片欧美黄色片| 超色免费av| 中文精品一卡2卡3卡4更新| 久久午夜福利片| 国产一级毛片在线| 不卡视频在线观看欧美| 久久99精品国语久久久| 久久久久久久精品精品| 狂野欧美激情性bbbbbb| 女人被躁到高潮嗷嗷叫费观| 国产成人精品一,二区| 18禁国产床啪视频网站| 国产成人精品在线电影| 日韩制服骚丝袜av| 夜夜骑夜夜射夜夜干| 高清在线视频一区二区三区| 精品国产露脸久久av麻豆| 亚洲国产av新网站| 最近手机中文字幕大全| 精品一区二区三卡| 纵有疾风起免费观看全集完整版| 国产成人精品婷婷| 亚洲欧美精品综合一区二区三区 | 波多野结衣一区麻豆| 欧美日韩亚洲国产一区二区在线观看 | 在线观看www视频免费| 美女中出高潮动态图| 免费女性裸体啪啪无遮挡网站| 久久这里只有精品19| 亚洲色图综合在线观看| 美女视频免费永久观看网站| 亚洲经典国产精华液单| 最近最新中文字幕免费大全7| 午夜老司机福利剧场| av线在线观看网站| 永久网站在线| 精品午夜福利在线看| 春色校园在线视频观看| 久久婷婷青草| 日韩av在线免费看完整版不卡| 欧美中文综合在线视频| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| av免费观看日本| av国产精品久久久久影院| 婷婷成人精品国产| 九草在线视频观看| 美女国产视频在线观看| 春色校园在线视频观看| 国产午夜精品一二区理论片| 欧美精品一区二区免费开放| 好男人视频免费观看在线| 久久精品国产亚洲av天美| 有码 亚洲区| 国产精品偷伦视频观看了| 免费观看性生交大片5| 伦理电影大哥的女人| 午夜福利影视在线免费观看| 精品国产露脸久久av麻豆| 在现免费观看毛片| 日本色播在线视频| 国产亚洲一区二区精品| 只有这里有精品99| 少妇人妻精品综合一区二区| 亚洲av男天堂| 欧美+日韩+精品| 久久韩国三级中文字幕| 国产一区二区三区av在线| 国产日韩欧美在线精品| 亚洲婷婷狠狠爱综合网| 成年女人毛片免费观看观看9 | 卡戴珊不雅视频在线播放| 菩萨蛮人人尽说江南好唐韦庄| 老鸭窝网址在线观看| 日韩一区二区视频免费看| 人妻系列 视频| 99久久综合免费| 可以免费在线观看a视频的电影网站 | 亚洲国产精品成人久久小说| 国产亚洲午夜精品一区二区久久| 日产精品乱码卡一卡2卡三| 最近中文字幕2019免费版| 国产精品国产av在线观看| 女人被躁到高潮嗷嗷叫费观| 午夜福利视频精品| 高清av免费在线| 久久99蜜桃精品久久| 91精品三级在线观看| 在线亚洲精品国产二区图片欧美| 美国免费a级毛片| 天堂8中文在线网| 亚洲一级一片aⅴ在线观看| 香蕉丝袜av| 日本91视频免费播放| 亚洲综合精品二区| 国产亚洲最大av| 制服诱惑二区| 国产有黄有色有爽视频| 九九爱精品视频在线观看| 日本爱情动作片www.在线观看| 午夜日本视频在线| 午夜日韩欧美国产| av福利片在线| 女人久久www免费人成看片| 欧美成人精品欧美一级黄| 日日爽夜夜爽网站| 欧美精品一区二区免费开放| 国产成人精品在线电影| tube8黄色片| a级毛片黄视频| av卡一久久| 一区二区av电影网| 精品国产乱码久久久久久男人| 美女午夜性视频免费| 国产乱来视频区| 中文字幕精品免费在线观看视频| 免费久久久久久久精品成人欧美视频| 在线观看免费日韩欧美大片| 免费女性裸体啪啪无遮挡网站| 18+在线观看网站| 欧美日韩亚洲高清精品| 免费在线观看视频国产中文字幕亚洲 | 欧美精品亚洲一区二区| 我的亚洲天堂| 精品国产一区二区久久| 高清黄色对白视频在线免费看| 免费看不卡的av| 亚洲精品中文字幕在线视频| 热re99久久国产66热| 国产男人的电影天堂91| 如日韩欧美国产精品一区二区三区| 精品酒店卫生间| 一级爰片在线观看| 久久精品久久久久久久性| 久久国内精品自在自线图片| 丰满迷人的少妇在线观看| 中文字幕色久视频| 亚洲精品,欧美精品| 久久99一区二区三区| 99热全是精品| 国产成人91sexporn| 人体艺术视频欧美日本| 久久精品国产亚洲av涩爱| 黑人巨大精品欧美一区二区蜜桃| 一个人免费看片子| 欧美激情极品国产一区二区三区| 男的添女的下面高潮视频| 免费在线观看黄色视频的| 少妇的丰满在线观看| 久久综合国产亚洲精品| 一级毛片我不卡| 日韩中文字幕欧美一区二区 | 亚洲av电影在线观看一区二区三区| 午夜老司机福利剧场| 一区二区三区激情视频| 老女人水多毛片| 国产探花极品一区二区| 男的添女的下面高潮视频| 爱豆传媒免费全集在线观看| 我的亚洲天堂| 久久久久国产精品人妻一区二区| 在线亚洲精品国产二区图片欧美| 热99久久久久精品小说推荐| 99国产综合亚洲精品| 亚洲欧美成人综合另类久久久| 黑人欧美特级aaaaaa片| 日韩 亚洲 欧美在线| 国产野战对白在线观看| 久久狼人影院| 最近中文字幕2019免费版| 波多野结衣一区麻豆| 97在线视频观看| 9191精品国产免费久久| 国产无遮挡羞羞视频在线观看| www.自偷自拍.com| 亚洲精品久久成人aⅴ小说| 女人久久www免费人成看片| 成人黄色视频免费在线看| 午夜激情久久久久久久| 99九九在线精品视频| 丝袜在线中文字幕| 久久久久久久亚洲中文字幕| 黄色配什么色好看| 这个男人来自地球电影免费观看 | 国产97色在线日韩免费| 色播在线永久视频| 韩国高清视频一区二区三区| 美女午夜性视频免费| 日本免费在线观看一区| 少妇人妻 视频| 午夜免费鲁丝| 人妻少妇偷人精品九色| 精品亚洲乱码少妇综合久久| 捣出白浆h1v1| 边亲边吃奶的免费视频| 亚洲av.av天堂| 久久久久国产一级毛片高清牌| 国产综合精华液| 久久久精品区二区三区| 巨乳人妻的诱惑在线观看| 欧美激情极品国产一区二区三区| 91午夜精品亚洲一区二区三区| 爱豆传媒免费全集在线观看| 免费观看无遮挡的男女| 精品国产超薄肉色丝袜足j| 久久国产精品男人的天堂亚洲| 女性生殖器流出的白浆| 国产极品天堂在线| 国产成人精品婷婷| 亚洲综合精品二区| 男人舔女人的私密视频| 1024香蕉在线观看| 日韩人妻精品一区2区三区| 久久久久久伊人网av| 卡戴珊不雅视频在线播放| 亚洲三区欧美一区| 久久精品熟女亚洲av麻豆精品| 国产成人精品在线电影| 男人操女人黄网站| 国产精品女同一区二区软件| a 毛片基地| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品一区在线观看| 国产熟女午夜一区二区三区| 国产精品一区二区在线不卡| 欧美日韩精品网址| 国产精品麻豆人妻色哟哟久久| 亚洲经典国产精华液单| av免费观看日本| 美女xxoo啪啪120秒动态图| 亚洲欧美一区二区三区国产| 国产亚洲av片在线观看秒播厂| 精品酒店卫生间| 亚洲欧洲国产日韩| 久久久久人妻精品一区果冻| 中文乱码字字幕精品一区二区三区| 日本91视频免费播放| 欧美日韩综合久久久久久| 国产一区二区三区av在线| 欧美少妇被猛烈插入视频| a 毛片基地| 赤兔流量卡办理| 精品福利永久在线观看| 日韩三级伦理在线观看| 精品国产国语对白av| 飞空精品影院首页| 久久久久精品人妻al黑| 国产黄色免费在线视频| 国产精品 欧美亚洲| 最近2019中文字幕mv第一页| 国产xxxxx性猛交| 亚洲第一区二区三区不卡| 中国三级夫妇交换| 美国免费a级毛片| 十八禁高潮呻吟视频| 国产男女超爽视频在线观看| 国产成人91sexporn| 国产精品久久久久成人av| 午夜福利视频在线观看免费| 你懂的网址亚洲精品在线观看| 国产成人免费观看mmmm| 大片免费播放器 马上看| 国产爽快片一区二区三区| 亚洲成色77777| 观看av在线不卡| 精品少妇久久久久久888优播| 丝袜在线中文字幕| 人人妻人人爽人人添夜夜欢视频| 成人手机av| 亚洲色图综合在线观看| 香蕉精品网在线| 久久精品国产亚洲av涩爱| 美女主播在线视频| 免费观看无遮挡的男女| 熟妇人妻不卡中文字幕| 性高湖久久久久久久久免费观看| 免费播放大片免费观看视频在线观看| 男男h啪啪无遮挡| 一本—道久久a久久精品蜜桃钙片| av天堂久久9| 国产精品国产三级国产专区5o| 亚洲国产精品999| 亚洲人成电影观看| 美女xxoo啪啪120秒动态图| 欧美日韩亚洲高清精品| 国产一区二区激情短视频 | 欧美人与性动交α欧美软件| 欧美精品高潮呻吟av久久| 永久网站在线| 国产一区亚洲一区在线观看| 国产精品熟女久久久久浪| 一区福利在线观看| 男女啪啪激烈高潮av片| 伦理电影免费视频| 麻豆精品久久久久久蜜桃| 国产精品二区激情视频| 日本午夜av视频| 热re99久久国产66热| 久久久久久免费高清国产稀缺| 国产探花极品一区二区| 两性夫妻黄色片| 国产综合精华液| av在线app专区| 免费观看在线日韩| 国产激情久久老熟女| 亚洲三区欧美一区| 国产欧美亚洲国产| 成人国产av品久久久| 亚洲人成77777在线视频| 搡女人真爽免费视频火全软件| 亚洲精品国产av成人精品| 国产又爽黄色视频| 日韩av在线免费看完整版不卡| 波野结衣二区三区在线| 免费在线观看黄色视频的| 好男人视频免费观看在线| 91成人精品电影| 边亲边吃奶的免费视频| 国产男女内射视频| 性色av一级| 女性生殖器流出的白浆| 桃花免费在线播放| a 毛片基地| 高清视频免费观看一区二区| 少妇人妻 视频| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 日韩人妻精品一区2区三区| 久久ye,这里只有精品| 黄色 视频免费看| 免费大片黄手机在线观看| 69精品国产乱码久久久| 精品国产国语对白av| 国产免费又黄又爽又色| 亚洲一区二区三区欧美精品| 最近中文字幕高清免费大全6| 女人被躁到高潮嗷嗷叫费观| 成人亚洲精品一区在线观看| 一本—道久久a久久精品蜜桃钙片| 在线免费观看不下载黄p国产| 午夜福利,免费看| 又大又黄又爽视频免费| 精品久久久精品久久久| 日韩视频在线欧美| 婷婷色麻豆天堂久久| 中文字幕av电影在线播放| 伊人久久国产一区二区| 精品国产乱码久久久久久男人| 有码 亚洲区| 久久久久久免费高清国产稀缺| 日韩熟女老妇一区二区性免费视频| 成人国产av品久久久| 欧美日韩一级在线毛片| 王馨瑶露胸无遮挡在线观看| 亚洲欧美色中文字幕在线| 香蕉国产在线看| 美国免费a级毛片| 国产伦理片在线播放av一区| 菩萨蛮人人尽说江南好唐韦庄| 曰老女人黄片| 黄片小视频在线播放| 亚洲国产av新网站| 老司机亚洲免费影院| 99国产精品免费福利视频| 最近最新中文字幕免费大全7| 在线免费观看不下载黄p国产| 精品国产一区二区三区久久久樱花| 亚洲欧洲日产国产| 美女视频免费永久观看网站| 久久99一区二区三区| 久久国产亚洲av麻豆专区| 成年女人在线观看亚洲视频| 国产熟女午夜一区二区三区| 99久久精品国产国产毛片| 宅男免费午夜| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说| 国产片内射在线| 久久久久久久亚洲中文字幕| 亚洲美女黄色视频免费看| av电影中文网址| 考比视频在线观看| 少妇被粗大猛烈的视频| 97人妻天天添夜夜摸| 看十八女毛片水多多多| 国产片内射在线| 免费黄网站久久成人精品| 男的添女的下面高潮视频| 免费观看无遮挡的男女| h视频一区二区三区| 亚洲av成人精品一二三区| 婷婷色综合www| 亚洲图色成人| 日本-黄色视频高清免费观看| 搡老乐熟女国产| 男女午夜视频在线观看| 女性生殖器流出的白浆| 观看av在线不卡| av电影中文网址| 精品亚洲成a人片在线观看| 边亲边吃奶的免费视频| 久久99热这里只频精品6学生| 国产精品久久久av美女十八| 国产黄频视频在线观看| 人妻人人澡人人爽人人| 亚洲国产精品999| 精品午夜福利在线看| 男人舔女人的私密视频| 日韩av不卡免费在线播放| 丝袜美足系列| av视频免费观看在线观看| 搡女人真爽免费视频火全软件| 亚洲欧美精品综合一区二区三区 | 亚洲美女视频黄频| 一级a爱视频在线免费观看| 丝袜美腿诱惑在线| 成人漫画全彩无遮挡| tube8黄色片| 免费观看a级毛片全部| 丝袜脚勾引网站| 90打野战视频偷拍视频| 亚洲欧洲精品一区二区精品久久久 | 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| 亚洲经典国产精华液单| 国产高清国产精品国产三级| 国产野战对白在线观看| 免费观看a级毛片全部| 亚洲欧美一区二区三区国产| 90打野战视频偷拍视频| 国产亚洲最大av| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 好男人视频免费观看在线| 国产伦理片在线播放av一区| 久久久久视频综合| 国产成人欧美| 18在线观看网站| 日韩免费高清中文字幕av| 人妻一区二区av| 高清av免费在线| av在线app专区| 国产精品一国产av| 亚洲一区中文字幕在线| 久久国产精品大桥未久av| 观看美女的网站| 秋霞伦理黄片| 最近最新中文字幕免费大全7| 免费看av在线观看网站| 男女下面插进去视频免费观看| 丝袜美足系列| 欧美激情 高清一区二区三区| 七月丁香在线播放| 一边摸一边做爽爽视频免费| 夫妻性生交免费视频一级片| 国产精品秋霞免费鲁丝片| av又黄又爽大尺度在线免费看| 国产精品蜜桃在线观看| 男女边摸边吃奶| 免费在线观看视频国产中文字幕亚洲 | 精品亚洲成国产av| 欧美日韩一级在线毛片| 日韩 亚洲 欧美在线| 永久网站在线| 亚洲国产最新在线播放| 亚洲精品美女久久av网站| 免费av中文字幕在线| 精品亚洲成国产av| 亚洲欧美成人综合另类久久久| 亚洲av综合色区一区| 一级片免费观看大全| 我要看黄色一级片免费的| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 国产av一区二区精品久久| 亚洲精品乱久久久久久| 日韩制服丝袜自拍偷拍| 大香蕉久久成人网| 在线 av 中文字幕| 丝袜美腿诱惑在线| 亚洲,一卡二卡三卡| 久久久久久久亚洲中文字幕| 精品久久久久久电影网| 18禁动态无遮挡网站| 街头女战士在线观看网站| 在线 av 中文字幕| 国产有黄有色有爽视频| 免费观看无遮挡的男女| 亚洲成av片中文字幕在线观看 | 五月天丁香电影| 黑丝袜美女国产一区| 好男人视频免费观看在线| 亚洲国产日韩一区二区| 久久人人爽av亚洲精品天堂| 国产成人午夜福利电影在线观看| 国产精品人妻久久久影院| 国产亚洲精品第一综合不卡| 精品一区二区免费观看| 亚洲三区欧美一区| 男女边摸边吃奶| 久久久久精品性色| 欧美日韩综合久久久久久| 精品国产国语对白av| 亚洲成人av在线免费| 我的亚洲天堂| 高清视频免费观看一区二区| 久久热在线av| 18+在线观看网站| 国产亚洲午夜精品一区二区久久| 久久久久久久亚洲中文字幕| 女人精品久久久久毛片| 99久久人妻综合| 日本wwww免费看| 青春草国产在线视频| 欧美日韩综合久久久久久| 国产精品 欧美亚洲| 亚洲欧美清纯卡通| 久久久久久久亚洲中文字幕| 精品国产乱码久久久久久男人| 日本vs欧美在线观看视频| 啦啦啦视频在线资源免费观看| 国产精品一二三区在线看| 人妻 亚洲 视频| 男女边摸边吃奶| 大片免费播放器 马上看| 成年人免费黄色播放视频| 国产精品女同一区二区软件| 考比视频在线观看| 精品一区在线观看国产| 国产精品国产av在线观看| 一级片'在线观看视频| 亚洲国产毛片av蜜桃av| 久久久久久伊人网av| 欧美成人午夜免费资源| 亚洲美女视频黄频| 久久精品国产亚洲av天美| 亚洲色图 男人天堂 中文字幕| 超碰97精品在线观看| 91国产中文字幕| 性色av一级| 如日韩欧美国产精品一区二区三区| 久久精品久久久久久久性| 国产精品偷伦视频观看了| 建设人人有责人人尽责人人享有的| 亚洲精品自拍成人| 水蜜桃什么品种好| av网站免费在线观看视频| 97人妻天天添夜夜摸| 亚洲激情五月婷婷啪啪| 一边摸一边做爽爽视频免费| 两性夫妻黄色片| 老汉色∧v一级毛片| 欧美日韩精品成人综合77777| 三级国产精品片| 亚洲熟女精品中文字幕| 婷婷色av中文字幕| 69精品国产乱码久久久| 伊人久久大香线蕉亚洲五| 亚洲欧美成人综合另类久久久| 亚洲国产精品国产精品| 色吧在线观看| 久久国产精品男人的天堂亚洲| 日韩av免费高清视频| 热99久久久久精品小说推荐| 亚洲情色 制服丝袜| 亚洲久久久国产精品| 中文字幕av电影在线播放| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 色网站视频免费| 9191精品国产免费久久| 免费观看av网站的网址| 免费日韩欧美在线观看| 亚洲精品aⅴ在线观看| 韩国av在线不卡| 日本av免费视频播放| 丝袜在线中文字幕| 日本wwww免费看| 在线观看一区二区三区激情| 日本猛色少妇xxxxx猛交久久| 亚洲精华国产精华液的使用体验| 久久99热这里只频精品6学生| 久久女婷五月综合色啪小说| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 青春草国产在线视频| 亚洲熟女精品中文字幕| 日韩制服丝袜自拍偷拍| 成人国语在线视频| 日韩一区二区视频免费看| 美女xxoo啪啪120秒动态图| 成年人免费黄色播放视频| 91精品国产国语对白视频| 国产av码专区亚洲av| 欧美日韩av久久| 黄片无遮挡物在线观看| 丝瓜视频免费看黄片| 高清av免费在线| 免费看不卡的av| 久久久久久久精品精品| 久久精品人人爽人人爽视色| 午夜福利影视在线免费观看| 国产福利在线免费观看视频| 在线免费观看不下载黄p国产| 欧美日韩精品成人综合77777| tube8黄色片| 久久ye,这里只有精品|