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

    山丘區(qū)中小河流邊灘植被分布驅(qū)動因子及響應(yīng)關(guān)系

    2017-06-06 12:00:34余根聽夏繼紅畢利東汪穎俊林立懷曹偉杰伊紫函
    中國水土保持科學(xué) 2017年2期
    關(guān)鍵詞:灘地植被因子

    余根聽,夏繼紅?,畢利東,汪穎俊,林立懷,曹偉杰,伊紫函

    (1.河海大學(xué)水利水電學(xué)院,210098,南京;2.龍游縣水利局,324400,浙江龍游)

    山丘區(qū)中小河流邊灘植被分布驅(qū)動因子及響應(yīng)關(guān)系

    余根聽1,夏繼紅1?,畢利東1,汪穎俊2,林立懷1,曹偉杰1,伊紫函1

    (1.河海大學(xué)水利水電學(xué)院,210098,南京;2.龍游縣水利局,324400,浙江龍游)

    河岸帶灘地植被是維持河流系統(tǒng)健康的重要屏障。以浙江省龍游縣靈山港邊灘為研究對象,運用聚類分析和RDA約束性冗余分析等方法,分析中小河流邊灘植被分布特征及其主要驅(qū)動因子,并建立各驅(qū)動因子對植被分布的響應(yīng)關(guān)系。結(jié)果表明:靈山港灘地植被種類豐富多樣,共有62種植被(被子植物61種,蕨類植物1種),其植被群落結(jié)構(gòu)包括9種類型;灘地相對水面高差因子、灘地形態(tài)系數(shù)和水文特性是植被分布空間差異性的主要影響因子,各因子貢獻率分別為37.50%、27.50%和16.82%,各因子對植被分布的驅(qū)動響應(yīng)機制不同。其中:在灘地相對水平面高差因子作用下,植被從灘地外緣到灘地內(nèi)部,耐水性由強到弱,豐富性由低到高,生物量由小到大;當(dāng)邊緣線發(fā)育系數(shù)SDI值處于2.0~4.0之間,且灘地狹長指數(shù)Pe/Pa值在0.12~0.3之間時,至少有60%的保證率使得灘地植被多樣性維持在1.03~1.96的高水平;水文特性能夠影響植被的形態(tài)特征和灘地的生境條件,水文作用下,灘地外圍砂粒質(zhì)量分?jǐn)?shù)隨著離水距離增大而逐漸減小,植被類型也會隨著離水距離的增加,從喜濕耐沖型到中生型過渡。掌握灘地植被的分布特征及其影響因子,以期為河道生態(tài)修復(fù),提供理論依據(jù)。

    山丘區(qū); 中小河流; 灘地; 植被分布特征; 驅(qū)動因子

    灘地植被是發(fā)揮灘地功能的重要因素,植被覆蓋的變化是生態(tài)環(huán)境變化的體現(xiàn)[1],一定寬度的灘地植被,可以通過過濾、滲透、吸收、滯留和沉積等機械、化學(xué)和生物功能效應(yīng),使進入地表和地下水的污染物毒性及污染程度降低,起到控制非點源污染、改善河流水文狀況和水質(zhì)的作用[2-3]。植被還可以通過改變微彎河岸擾動壓強大小,影響河岸帶灘地側(cè)向潛流交換和物質(zhì)能量傳遞,進而影響河流新陳代謝[4-5],并且植被分布下的灘地生態(tài)邊緣效應(yīng)、緩沖邊緣效應(yīng)和社會經(jīng)濟邊緣效應(yīng),對河道的建設(shè)管理和開發(fā)利用均具有深刻的影響[6]。近年來,由于人為采砂和侵占等干擾,灘地植被資源減少,功能嚴(yán)重退化,這一問題已引起社會各界的認(rèn)識和重視。例如:高敏等[7]通過二維數(shù)學(xué)模型,對植被化的漫灘水流進行流速分布研究;王月容等[8]對灘地植被緩沖帶進行劃分及功能和生態(tài)重建;閆靜等[9]、陳正兵等[10]分析了不同植物種類和不同排列方式對水流結(jié)構(gòu)的影響效應(yīng);Dwire等[11]、Burton等[12]分析不同環(huán)境梯度及城市化下的植被種類組成及結(jié)構(gòu)效應(yīng);Moffatt等[13]分析植被群落的種類組成、多樣性及周圍土地利用類型。但是,以往學(xué)者的研究多是圍繞植被對水流結(jié)構(gòu)的影響方面,或是圍繞研究所在地的植被基本情況及功能闡述方面,鮮有對植被的影響因子強度效應(yīng)進行具體分析?;诖耍P者以龍游縣靈山港灘地為研究對象,分析植被分布的特點及其主要驅(qū)動因子,深入探討不同因子對植被分布的響應(yīng)關(guān)系,以期為灘地生態(tài)修復(fù)中,植被的選擇與布置提供理論依據(jù)。

    1 研究區(qū)概況

    浙江省龍游縣位于E119°01′41″~119°19′52″,N28°44′10″~29°17′15″之間,地處浙江省西部金衢盆地,屬浙中丘陵盆地。靈山港源于海拔1 265 m的遂昌縣和尚嶺,在沐塵鄉(xiāng)馬戍口流入龍游縣境內(nèi),由南向北貫穿龍游縣中部,于龍游縣城匯入衢江,是錢塘江上游衢江右岸的一條重要支流。靈山港流域?qū)賮啛釒Ъ撅L(fēng)氣候區(qū),雨量豐沛,四季變化明顯。多年平均氣溫17.3 ℃,年平均日照時間1 966 h,多年平均降水量為1 666.4 mm。河道在龍游縣境內(nèi)流域面積367.6 km2,主流長43.79 km,平均比降2.45‰,是典型的山丘區(qū)中小河流。河道沿線灘地資源豐富,據(jù)調(diào)查統(tǒng)計,2014年河道沿線灘地總面積為1.66 km2。

    2 材料與方法

    2.1 監(jiān)測點布設(shè)與取樣

    灘地植被通過布置監(jiān)測樣方進行調(diào)查。研究將靈山港分為上、中、下游與河口4個區(qū)段,每個區(qū)段內(nèi)選取3個典型灘地(河口段除外),上游區(qū)段包括溪口四橋灘地(L1)、江潭灘地(L2)和下徐橋灘地(L3);中游區(qū)段包括寺下灘地(L4)、梅村灘地(L5)和周村灘地(L6);下游區(qū)段包括姜席堰灘地(L7)、上揚村灘地(L8)和高鐵橋灘地(L9);河口段為彩虹橋灘地(L10)。典型灘地分布見圖1(a)。每個典型灘地內(nèi)布設(shè)一定數(shù)量的喬木樣方(10 m×10 m)、灌木樣方(5 m×5 m)和草本樣方(1 m×1 m)。草本樣方布置于灘頭、灘中和灘尾斷面上,每個斷面上按照離水邊緣距離大小,每隔7.5 m布置1個草本樣方;同時,通過野外流速儀,監(jiān)測典型灘地邊緣水流流速。植被監(jiān)測樣方和流速監(jiān)測點布置見圖1(b)。

    L1:溪口四橋灘地,L2:江潭灘地,L3:下徐橋灘地,L4寺下灘地,L5:梅村灘地,L6:周村灘地,L7:姜習(xí)堰灘地,L8:上揚村灘地,L9:高鐵橋灘地,L10:彩虹橋灘地。以下類同。L1: Xikousiqiao floodplain, L2: Jiangtan floodplain, L3: Xiaxuqiao floodplain, L4: Sixia floodplain, L5: Meicun floodplain, L6: Zhoucun floodplain, L7: Jiangxiyan floodplain, L8: Shangyangcun floodplain, L9: Gaotieqiao floodplain, and L10: Caihongqiao floodplain. The same below.圖1 監(jiān)測點斷面及監(jiān)測點樣方布設(shè)Fig.1 Map of monitoring sections and plots

    對喬木與灌木,主要記錄樣方內(nèi)其數(shù)量、種類、樹高、胸徑和冠幅等參數(shù);對草本植被,主要通過摘取、分類和分裝,清點株數(shù),并帶回實驗室,經(jīng)稱重、殺青、烘干和再稱重等處理,獲得植被鮮質(zhì)量和干質(zhì)量。取樣調(diào)查分別于2015年4月、2015年10月和2016年5月進行。

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

    2.2.1 主要植被參數(shù)及處理 根據(jù)監(jiān)測數(shù)據(jù),分別計算各灘地監(jiān)測斷面草本植被生物多樣性指數(shù)[14]和生物量[15],并用均值代表灘地植被指標(biāo)參數(shù)。對于喬木與灌木,主要計算喬灌木的空間測度指數(shù)(混交度)[16]、分布均勻度(角尺度)[17]和葉面積指數(shù)(LAI)[18]等參數(shù)。為了獲取植被空間變異規(guī)律,采用半變異函數(shù)[19]進行分析,描述植被空間連續(xù)變異,反映植被不同距離測算值之間的變化。

    2.2.2 環(huán)境因子數(shù)據(jù)處理 針對灘地的特點,主要選取灘地相對水面高差(ButtElev)、土壤吸濕系數(shù)(HygrCoef)、灘地離水距離(Distance)、灘地形態(tài)系數(shù)(Shapindx)及水文特性(HydrChar)5個環(huán)境因子,通過現(xiàn)場監(jiān)測及Google earth測量,獲得相關(guān)因子的數(shù)據(jù)。其中:ButtElev選取灘地植被樣方所在位置與河道水平面之間的垂直距離表示;HygrCoef選取現(xiàn)場樣方位置表層20 cm土層測算;Distance選取監(jiān)測樣方與水邊的水平距離表示見圖1(b);Shapindx選取歐式幾何形態(tài)指標(biāo)SDI及灘地狹長性指數(shù)Pe/Pa(短軸/長軸)表示;HydrChar選取弗勞德數(shù)(Fr)與紊動能(Eε)表示。

    3 結(jié)果與分析

    3.1 植被區(qū)系構(gòu)成及分布特征

    3.1.1 總體特征 現(xiàn)場調(diào)查發(fā)現(xiàn),靈山港沿線灘地共有61種被子植物,隸屬于25科56屬。其中,絕大部分為草本植物,主要有菊科、禾本科、蓼科、傘形科、莧科、石竹科、莎草科、蕁麻科、豆科、玄參科、唇形科和茜草科等共58種,隸屬于23科54屬,菊科和禾本科分別占14.52%和17.74%,蓼科和傘形科分別占8.06%和6.45%,莧科和石竹科均占4.84%,莎草科、蕁麻科、豆科、玄參科、唇形科和茜草科均占3.23%。僅發(fā)現(xiàn)1科1屬蕨類植物,為陰石蕨(Humatarepens(L. f.) Diels)。另外,發(fā)現(xiàn)少量的核桃楸(JuglansmandshuricaMaxim.)、竹子(Bambusoideae)以及楊樹(PopulusL.)等木本植物,共3科3屬。

    3.1.2 草本群落結(jié)構(gòu)特點 靈山港灘地植被主要為草本植物,植被種類豐富,其群落結(jié)構(gòu)存在較大差異。為準(zhǔn)確劃分草本植被的群落結(jié)構(gòu),應(yīng)用SPSS,對典型灘地內(nèi)草本植物種類進行聚類分析(表1)。由表1可知,靈山港灘地草本植被可以分為9類群落結(jié)構(gòu)。從分布頻率看,沿線灘地草本植物中出現(xiàn)頻率較高的為I、II和III類群落;其中,I、II類群落主要為間接分布種,III類群落為廣布種。可見,灘地植被中廣布種和間接分布種最多,其次是內(nèi)部分布種;另外,緊鄰灘地外緣分布物種種類較少,主要為水蓼(Polygonumhydropiper)、積雪草(Centellaasiatica(L.) Urban)、沿階草(OphiopogonbodinieriLevl.)等。由此可知,草本植物分布和出現(xiàn)頻率存在較大差異,通過現(xiàn)場觀測發(fā)現(xiàn),產(chǎn)生這一現(xiàn)象的主要原因是受灘地幾何形態(tài)、地形、地貌和水流條件等因素的影響。

    表1 灘地草本植被空間分布類型

    3.1.3 林分空間結(jié)構(gòu)特點 靈山港沿線灘地中,木本植物有其獨特的林分結(jié)構(gòu)?,F(xiàn)場觀測發(fā)現(xiàn),木本植物主要出現(xiàn)于溪口四橋灘地(L1)、寺下灘地(L4)、周村灘地(L6)和彩虹橋灘地(L10)這4個灘地中。各典型灘地中,木本植物的林分指標(biāo)計算值如表2所示,木本植被個體胸徑、株高和LAI指數(shù)都不同,其周邊的草本植被類型也存在較大差異。灘地的木本LAI指數(shù)分別為30.11、28.04、34.74和44.84,基本呈現(xiàn)增長趨勢,2項擬合度較高,相關(guān)系數(shù)達到0.991 4。其中,L4灘地的木本LAI指數(shù)略有減小,這是由于灘地破碎化程度較高,使植被分布密度受到極大的限制。從混交度和角尺度上看,混交度均小于0.05,木本植被單一化程度明顯,角尺度雖有些許波動,但二者均值均小于0.475,都為均勻分布類型??傮w而言,靈山港灘地木本植被產(chǎn)量較為豐富,分布類型單一,有一定的局部小氣候效應(yīng)特征。

    表2 林分空間結(jié)構(gòu)特征

    3.2 驅(qū)動因子對植被分布的響應(yīng)關(guān)系

    3.2.1 灘地植被空間關(guān)聯(lián)性分析及關(guān)鍵驅(qū)動因子的確定 應(yīng)用半變異函數(shù)來研究植被多樣性的空間異質(zhì)性,并確定其影響范圍大小,計算結(jié)果由表3可知,L1、L5和L10灘地自相關(guān)性比較大,分別為0.266、0.222和0.253,表明隨機性因子(局地氣候、氣象條件和局地水文特性等)在空間變異中占據(jù)主導(dǎo)地位;其余灘地自相關(guān)性較小,表明結(jié)構(gòu)因子(地形地貌變化、地質(zhì)構(gòu)造和灘地形態(tài)等)在空間變異中起主導(dǎo)作用?;诖耍瑢`山港灘地植被空間分布的驅(qū)動因子,及其敏感性大小進行深入分析。

    靈山港灘地植被空間分布受多個因子影響。對土壤吸濕系數(shù)、灘地離水距離、灘地形態(tài)系數(shù)、灘地相對水面高差因子和水文特性,進行KMO和Bartlett球形度檢驗,得KMO>0.7、P<0.05,表明因子選擇具有結(jié)構(gòu)效度,可做因子分析。通過Canoco軟件RDA約束性冗余分析方法,分別計算各因子對植被分布的貢獻率大小,通過貢獻率,反映各因子對植被分布的驅(qū)動程度,各因子貢獻率結(jié)果見圖2,為了可觀性較好,圖中只展示了部分表1中的物種。由圖2 RDA雙序圖可知,相對水面高差因子、形態(tài)系數(shù)和水文特性對植被分布的貢獻率分別為37.50%、27.50%和16.82%,三者總體貢獻率達到81.82%。土壤吸濕系數(shù)和灘地離水距離的貢獻率分別為10.56%和7.62%,二者總體貢獻率為18.18%。圖中藍色箭頭代表草本植被種類,其長度在因子變量上的投影值代表物種與環(huán)境因子之間相關(guān)性的大小,若投影方向與因子正方向相同,表示物種與環(huán)境因子之間是正相關(guān),反之為負(fù)相關(guān)。從貢獻率和投影長度看,主要驅(qū)動因子對植被分布的敏感性程度順序為相對水面高差因子>形態(tài)系數(shù)>水文特性。

    表3 灘地植被分布自相關(guān)性表

    注:Cn:塊金方差反映樣點間多樣性的變化程度;a:變程值用來衡量兩樣本間的相關(guān)性,距離d

    1: 狗牙根Cynodon dactylon(Linn.)Pers. 2: 棒頭草 Polypogon fugax Nees ex Steud. 3: 大畫眉草Eragrostis cilianensis (All.) Link ex Vignolo-Lutati. 5: 看麥娘 Alopecurus aequalis Sobol. 6: 菵草Beckmannia syzigachne(Steud.)Fern. 7: 刺蓼Polygonum senticosum (Meisn.) Franch. et Sav. 10: 小飛蓬 Conyza canadensis (L.) Cronq. 11: 鼠麯 Gnaphalium affine D. Don. 13: 苣荬菜 Sonchus arvensis L. 16: 積雪草Centella asiatica (L.) Urban. 19: 水蓼Polygonum hydropiper. 21: 空心蓮子草 Alternanthera Philoxeroides(Mart.)Griseb. 22: 刺莧 Amaranthus spinosus. 24: 毒芹 Cicuta virosa L. 26: 毛豆 Glycine max. 27: 陰石蕨 Humata repens (L. f.) Diels. 28: 毛茛 Ranunculus japonicus Thunb. 29: 通泉草 Mazus japonicus (Thunb.) O. Kuntze. 35: 雪見草 Salvia plebeia R. Br. 36: 旋覆花 Inula japonica Thunb. 37: 假稻 Leersia japonica (Makino) Honda. 38: 鴨跖草 Commelina communis. 40: 野老鶴草 Geranium carolinianum L. 41: 黃鵪菜Youngia japonica. 43: 簇生卷耳Cerastium caespitosum Gilib. The same below. ButtElev:Elevation factor; Shapindx:Shape coefficient; HydrChar:Hydrological characteristics; HygrCoef:Soil moisture absorption coefficient; Distance:Distance away from the water. 圖2 物種與環(huán)境因子的RDA貢獻率雙序圖Fig.2 RDA double sequence diagram of species and environmental factors

    3.2.2 植被分布對灘地相對水面高差因子變化的響應(yīng)關(guān)系 灘地相對水面高差因子是影響植被分布最顯著的因子,以靈山港上游溪口四橋灘地(L1)、中游周村灘地(L6)、下游高鐵橋灘地(L9)以及衢江口彩虹橋灘地(L10)為例,說明不同區(qū)段灘地內(nèi),相對水面高差因子對植被分布的影響。4個灘地內(nèi),植被分布的相對水平面高差云圖見圖3~圖6。

    圖3 上游區(qū)段植被空間分布Fig.3 Spatial distribution of vegetation in the upper reaches

    圖4 中游區(qū)段植被空間分布Fig.4 Spatial distribution of vegetation in the middle reaches

    圖5 下游區(qū)段植被空間分布Fig.5 Spatial distribution of vegetation in the lower reaches

    圖6 河口植被空間分布Fig.6 Spatial distribution of vegetation in the estuary

    可見,I、II、IV、VI、VIII、IX類植被主要生長在高差較大的灘地內(nèi)部,滿足其排水和向陽條件,III、V、VII類植被主要生長在高差較小的臨水邊,滿足其水分需求較高的特點。各類植被的分布特點也存在一定的差異。I類植被(通泉草、牛毛氈、鼠麯、小飛蓬和狗牙根等)廣泛分布于高差0.25~1.2 m內(nèi),且易形成大片群落。主要是因為該區(qū)間內(nèi),土壤受水流作用,土層松動,沙性明顯,水分條件及排水較好,適宜I類植被生長。II類植被(小花糖芥、細(xì)風(fēng)輪菜和旋覆花等)分布范圍較廣,主要分布在高差>1.8 m的位置。IV類植被(艾草、野豌豆和北美車前等)大多分布于高差1.8~2.1 m內(nèi),少部分會出現(xiàn)在臨水邊0.36~0.54 m處。V和VII類植被(沿階草、積雪草、水蓼等)主要分布于高差<0.9 m的臨水邊,主要是因為積雪草葉片中的機械組織不發(fā)達,抗旱能力極差,是陰生濕生植物,而蓼科和莎草科等陽生濕生植被,根系不發(fā)達,沒有根毛,但根與莖之間有通氣的組織,以保證取得充足的氧氣,生活在陽光充足、土壤水分飽和的沼澤地區(qū)或湖邊。III類植被(棒頭草、假稻、鴨跖草和酸模葉蓼等)會與V和VII類植被集群分布。VI類植被(刺蓼、陰石蕨和野菊等)主要分布于高差0.4~1.25 m(刺蓼、山大顏和毒芹)、1.8~2.4 m(陰石蕨和野菊)的位置。VIII和IX類植被主要分布于高差1.8~2.0 m的位置,在灘地內(nèi)緣有少量出現(xiàn)。

    由圖7和圖8可知,灘地植被從外緣到內(nèi)緣(高差從小到大),生物量在逐漸增大,且植被多樣性隨著灘地相對水面高度的增加而增加,對應(yīng)的植被擴散斑塊(圓域面積)大小在逐漸擴增。總體而言,地形高度的差異,可以作為植被種類的分界線,灘地植被從灘地外緣到灘地內(nèi)部的分布規(guī)律為:耐水性由強到弱,豐富性由低到高,生物量由小到大。

    圖7 灘地外緣至內(nèi)緣植被生物量變化Fig.7 Biomass variation from outer to inner edge of floodplain

    圖8 不同相對水面高度下的植被多樣性Fig.8 Vegetation Shannon-Wiener under different relative elevation of water

    3.2.3 植被分布對形態(tài)系數(shù)的響應(yīng)關(guān)系 灘地總體形態(tài)變遷是影響植被多樣性的另一個重要因素,穩(wěn)定的灘地生境是植被分布特征優(yōu)化,以及發(fā)揮生態(tài)效益最大化的重要保障。由圖9可知,隨著灘地Pe/Pa指數(shù)的增大,植被覆蓋率也隨之增大。從圖3~6亦可看出,短寬和窄長灘地對等高線分布及灘地斷面植被容納量有明顯的影響。選擇灘地邊緣線發(fā)育系數(shù)SDI(即Richardson所提出的緊湊度C0的倒數(shù),SDI=1時,灘地形態(tài)近似圓形)和灘地橫縱徑(短長軸)比值Pe/Pa(避免SDI較大時,灘地出現(xiàn)過分窄長的情況),來反映灘地形態(tài)的彎曲和狹長。Pe/Pa和SDI與植被多樣性響應(yīng)關(guān)系分別見圖9和圖10。

    圖9 Pe/Pa與植被覆蓋度和多樣性的響應(yīng)關(guān)系Fig.9 Response of coverage and Shannon-Wiener to Pe/Pa

    圖10 SDI與植被多樣性的響應(yīng)關(guān)系Fig.10 Response of Shannon-Wiener to SDI

    由圖10可知,靈山港灘地SDI系數(shù)>2的占灘地總數(shù)的70%,說明灘地緊湊度不大,岸線幾何形狀較為復(fù)雜。在自然情況下,灘地植被多樣性隨著SDI系數(shù)的變化而變化,兩者呈現(xiàn)同樣的趨勢效應(yīng),節(jié)點變化同一性極高。當(dāng)SDI值處于2.0~4.0之間時,灘地平均植被多樣性>1出現(xiàn)的頻率占全局60%以上。圖9顯示,灘地植被多樣性隨著Pe/Pa的增大而減小;由Q-Q檢測可知,各點分布于直線附近,且K-S檢驗的漸近顯著性系數(shù)為0.580>0.05。因此,靈山港灘地的形態(tài)特征服從正態(tài)分布規(guī)律,灘地整體形態(tài)穩(wěn)定性較好,而且,當(dāng)Pe/Pa值在0.12~0.3之間時,植被多樣性維持在1.03~1.96之間的次數(shù)出現(xiàn)6次占60%。即至少有60%的保證率使得SDI系數(shù)處于2.0~4.0之間,且灘地狹長指數(shù)Pe/Pa值在0.12~0.3之間時,灘地植被多樣性較高。

    3.2.4 植被分布對水文特性的響應(yīng)關(guān)系 從植被空間變異特征(圖11)可知,距離水邊10 m以內(nèi),指數(shù)的平均差異性Δr(h)約為0.100,空間變異幅度較小,這主要是由于灘地臨水邊緣小區(qū)域受水流影響較大,植被物種差異性不明顯,多為喜濕耐沖型植被,將這一區(qū)域稱為低變幅區(qū)。距離水邊10~25 m時,指數(shù)的平均差異性Δr(h)達0.206,空間變異幅度較大,這主要是由于微地貌的變化,使得灘地淹沒所需水流深度增加,這一因素的限制,在平水位情況下,水力沖蝕削弱,主要生長中生型植被,植被種類豐富,將這一區(qū)域稱為高變幅區(qū)。距離水邊>25 m時,指數(shù)的平均差異性Δr(h)僅為0.052左右,無明顯的空間變異性(寺下、上揚村受人為干擾,變異性依舊較大),此時結(jié)構(gòu)因子正在逐漸替代水文特性成為主導(dǎo)因子,且由表3可知,變程值a均>25 m;因此,將這一區(qū)域稱為平穩(wěn)區(qū),主要為喬灌草結(jié)構(gòu)。

    圖11 植被空間相關(guān)性曲線Fig.11 Vegetation spatial correlation curve

    將低變幅區(qū)與高變幅區(qū),即離水邊緣25 m以內(nèi)的植被帶,作為植被分布對水文特性響應(yīng)關(guān)系的研究區(qū)域。各區(qū)段灘地外緣植被帶土壤性質(zhì)、植被多樣性和生物量的變化見表4。

    由表4可知,同一灘地低變幅區(qū)土壤中砂粒質(zhì)量分?jǐn)?shù)比高變幅區(qū)的大,且低變幅區(qū)的土壤質(zhì)地類別中,壤質(zhì)砂土出現(xiàn)頻率極高,達到80%。究其原因,土壤顆粒組成受水文特性的影響極大,水流沖刷和淘蝕灘地岸線,土壤中的粉粒和黏粒等細(xì)粒物質(zhì)大量流失,越靠近河流的淹沒帶,越容易受到其影響,致使粉黏粒質(zhì)量分?jǐn)?shù)下降,砂粒質(zhì)量分?jǐn)?shù)增加。水流作用改變了植物根系與土壤之間物質(zhì)交換屬性,使得植被的空間分布也存在異質(zhì)性。

    表4 植被分布對水文特性的響應(yīng)

    從植被形態(tài)可塑性角度看,為了適應(yīng)高流速、大紊動的水流條件,灘地外緣植被的葉片和桿莖均為輕柔和狹長,能夠順?biāo)鞣较騼A伏,以克服水流的拖拽力及紊動卷攜,從而保持在該條件下的生存能力。表現(xiàn)為在灘地外緣25 m帶寬范圍內(nèi),蓼科和棒頭草出現(xiàn)的頻率均占40%,菵草和狗牙根出現(xiàn)的頻率均占30%,沿階草出現(xiàn)的頻率占10%,均為喜濕耐沖型植被,與內(nèi)緣植被差異性較高。從植被生長屬性角度看,靈山港灘地基本呈現(xiàn)上、下游區(qū)段植被多樣性較高,生物量較低,而中游區(qū)段植被多樣性較低,生物量較高的現(xiàn)象。這是由于上、下游水流適度的擾動及短期水位變動,為低多樣性集群物種的存活率創(chuàng)造了條件,在主要植被類型基礎(chǔ)上,新增細(xì)風(fēng)輪菜、通泉草、小飛蓬、鼠麯、小花糖芥、鴨跖草和積雪草等植被物種,使得植被帶寬內(nèi)多樣性指數(shù)較大;但由于物種莖稈及根系較小,生物量相對較低。而中游在集中的木本植被阻流緩沖和局地小氣候效應(yīng)影響下,植被帶寬內(nèi)水文效應(yīng)薄弱,灘地生境穩(wěn)定,種類單一化明顯,多為陰石蕨、沿階草和毒芹等,使得植被多樣性較?。坏捎邶嫶蟮男K根及高含水量,使得生物量較大。

    4 結(jié)論

    1)由植被空間關(guān)聯(lián)性分析,將靈山港灘地分為低變幅區(qū)、高變幅區(qū)和平穩(wěn)區(qū)3個區(qū)域。其中,植被分布受隨機因子與結(jié)構(gòu)因子的共同影響,在灘地外緣至岸邊,植被類型表現(xiàn)為喜濕耐沖型植被、中生植被、中生植被+喬木+灌木的過渡。

    2)山丘區(qū)河道灘地植被分布的關(guān)鍵驅(qū)動因子為灘地相對水面高差因子、灘地形態(tài)系數(shù)和水文特性。從整體范圍上看,植被分布對關(guān)鍵驅(qū)動因子敏感性大小依次為灘地相對水面高差因子>灘地形態(tài)系數(shù)>水文特性。

    3)地形高度的差異可以作為植被種類的分界線,灘地植被的分布從灘地外緣到灘地內(nèi)部,耐水性由強到弱,豐富性由低到高,生物量由小到大。

    4)灘地不可能無限發(fā)育,有其自身的極限展長與極限展寬。當(dāng)SDI系數(shù)處于2.0~4.0之間,且灘地狹長指數(shù)Pe/Pa值在0.12~0.3之間時,至少有60%的保證率,使得灘地植被多樣性能維持在較高水平。

    5)河道水文特性會從不同程度影響植被的空間關(guān)聯(lián)性,如植被形態(tài)特征與生境條件。水文特性作用下,常表現(xiàn)為灘地外圍植被帶砂粒質(zhì)量分?jǐn)?shù)隨離水距離增大而減少。外圍植被帶中,蓼科、棒頭草、菵草、狗牙根和沿階草等,對水文環(huán)境適應(yīng)能力極佳,可用于生態(tài)修復(fù)建設(shè)過程中優(yōu)先選擇的對象。

    [1] 廖建華,李丹勛,王興奎,等. 長江上游植被覆蓋的時空分異季節(jié)變化及其驅(qū)動因子研究[J].環(huán)境科學(xué)學(xué)報,2009,29(5): 1103. LIAO Jianhua, LI Danxun, WANG Xingkui, et al. The spatial-temporal distribution of seasonal vegetation changes and their driving forces in the upper reaches of the Yangtze River[J]. Acta Scientiae Circumstantiae, 2009, 29(5):1103.

    [2] 郭二輝,孫然好,陳利頂. 河岸植被緩沖帶主要生態(tài)服務(wù)功能研究的現(xiàn)狀與展望[J].生態(tài)學(xué)雜志,2011,30(8):1830. GUO Erhui, SUN Ranhao, CHEN Liding. Main ecological service functions in riparian vegetation buffer zone: research progress and prospects[J]. Chinese Journal of Ecology, 2011, 30(8):1830.

    [3] 孫彭成,高建恩,王顯文,等. 柳枝稷植被過濾帶攔污增效試驗初步研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2016,35(2):314. SUN Pengcheng, GAO Jianen, WANG Xianwen, et al. Effectiveness of switchgrass vegetative filter strip in intercepting pollutants and promoting plant biomass[J]. Journal of Agro- Environment Science, 2016, 35(2):314.

    [4] 林俊強,嚴(yán)忠民,夏繼紅. 微彎河岸沿線擾動壓強分布特性試驗[J].水科學(xué)進展,2013,24(6):855. LIN Junqiang, YAN Zhongmin, Xia Jihong. Laboratory experiments on pressure distribution along sinuous riverbanks[J]. Advances in Water Science, 2013, 24(6):855.

    [5] 夏繼紅,陳永明,王為木,等. 河岸帶潛流層動態(tài)過程與生態(tài)修復(fù)[J].水科學(xué)進展,2013,24(4):589. XIA Jihong, CHEN Yongmin, WANG Weimu, et al. Dynamic processes and ecological restoration of hyporheic layer in riparian zone[J]. Advances in Water Science, 2013, 24(4):589.

    [6] 夏繼紅,林俊強,姚莉,等. 河岸帶的邊緣結(jié)構(gòu)特征與邊緣效應(yīng)[J].河海大學(xué)學(xué)報(自然科學(xué)版),2010,38(2):215. XIA Jihong, LIN Junqiang, YAO Li, et al. Edge structure and edge effect of riparian zones[J]. Journal of Hohai University(Natural Sciences),2010,38(2):215.

    [7] 高敏,槐文信,趙明登. 灘地植被化的順直復(fù)式河道漫灘水流計算[J].節(jié)水灌溉,2008(6):21. GAO Min, HUAI Wenxin, ZHAO Mingdeng. Overbank flow calculation for straight compound channel with vegetated floodplain[J].Water Saving Irrigation,2008(6):21.

    [8] 王月容,周志翔,周金星,等. 長江灘地植被緩沖帶類型及功能與生態(tài)重建[J].湖北林業(yè)科技, 2010(4):5. WANG Yuerong, ZHOU Zhixiang, ZHOU Jinxing, et al. Type,function and ecological reconstruction of vegetated buffer strips in Yangtze River Beach[J]. Hubei Forestry Science and Technology, 2010(4):5.

    [9] 閆靜,唐洪武,田志軍,等. 植物對明渠流速分布影響的試驗研究[J].水利水運工程學(xué)報,2011(4):138. YAN Jing, TANG Hongwu, TIAN Zhijun, et al. Experimental study on the influence of vegetation on the velocity distribution of open channel flows[J]. Hydro-Science and Engineering, 2011(4):138.

    [10] 陳正兵,江春波.灘地植被對河道水流影響[J].清華大學(xué)學(xué)報(自然科學(xué)版), 2012,52(6):804. CHEN Zhengbing, JIANG Chunbo. Effect of floodplain vegetation on river hydrodynamics [J].Journal of Tsinghua University (Science and Technology), 2012,52(6):804.

    [11] DWIRE K A, KAUFFMAN J B, BROOKSHIRE E N J, et al. Plant biomass and species composition along an environmental gradient in montane riparian meadows[J]. Oecologia, 2004, 139(2):309.

    [12] BURTON M L, SAMUELSON L J. Influence of urbanization on riparian forest diversity and structure in the Georgia Piedmont, US[J]. Plant Ecology, 2008, 195(1):99.

    [13] MOFFATT S F, MCLACHLAN S M, KENKEL N C. Impacts of land use on riparian forest along an urban-rural gradient in southern Manitoba[J]. Plant Ecology, 2004, 174(174):119.

    [14] 李悅,馬溪平,李法云,等.細(xì)河河岸帶植物多樣性研究[J].廣東農(nóng)業(yè)科學(xué),2011(19):131. LI Yue, MA Xiping, LI Fayun, et al. Study on plant species diversity of riparian zone along Xi river[J].Guangdong Agricultural Sciences,2011(19):131.

    [15] 孟偉,張遠,渠曉東等.河流生態(tài)調(diào)查技術(shù)方法[M].北京:科學(xué)出版社,2011:65. MENG Wei, ZHANG Yuan, QU Xiaodong, et al.Technical methods of river ecological investigation[M]. Beijing: Science Press, 2011:65.

    [16] 惠剛盈,胡艷波,趙中華. 基于相鄰木關(guān)系的樹種分隔程度空間測度方法[J].北京林業(yè)大學(xué)學(xué)報,2008,30(4):131. HUI Gangying, HU Yanbo, ZHAO Zhonghua. Evaluating tree species segregation based on neighborhood spatial relationships[J]. Journal of Beijing Forestry University, 2008, 30(4):131.

    [17] 惠剛盈,Gadow K V,胡艷波.林分空間結(jié)構(gòu)參數(shù)角尺度的標(biāo)準(zhǔn)角選擇[J].林業(yè)科學(xué)研究,2004,17(6):687. HUI Gangying, Gadow K V, HU Yanbo. Theoptimum standard angle of the uniform angle index[J]. Forest Research, 2004, 17(6):687.

    [18] 麻雪艷,周廣勝. 春玉米最大葉面積指數(shù)的確定方法及其應(yīng)用[J].生態(tài)學(xué)報,2013,33(8):2596. MA Xueyan, ZHOU Guangsheng. Method of determining the maximum leaf area index of spring maize and its application[J]. Acta Ecologica Sinica, 2013, 33(8):2596.

    [19] 李斌.青藏高原植被時空分布規(guī)律及其影響因素研究[D]北京:中國地質(zhì)大學(xué)出版社,2016:62. LI Bin. Study of vegetation’s spatial and temporal distribution and influencing factors on the Tibetan Plateau[D]. Beijing: China University of Geosciences,2016:62.

    Vegetation distribution and its driven-forces on the floodplains of small and medium rivers in hilly area

    YU Genting1, XIA Jihong1, BI Lidong1, WANG Yingjun2, LIN Lihuai1, CAO Weijie1, YI Zihan1

    (1.College of Water Conservancy and Hydropower, Hohai University, 210098, Nanjing, China; 2.Water Conservancy Bureau, 324400, Longyou, Zhejiang, China)

    [Background] Vegetation in riparian zone is an important buffer to sustain the health of a river system. It is a typical ecotone characterized as edge effect, and plays a significant role in the conservation of river system health. Therefore, it is vital to know how the vegetation in floodplain distributes and what the driven-forces of the distribution are. [Methods] Taking the Lingshan River, in Longyou County, Zhejiang Province, as a typical example of small and medium rivers in hilly area, the characteristics of the vegetation distribution and its driven-forces were studied. After investigating the vegetation species in the quadrat sampling, the data was calculated and analyzed using the Clustering Analysis (CA) in SPSS software, and the Redundancy Analysis (RDA) in Canoco5 software and through which the composition and cluster numbers of vegetation species, the cosine of the angle at the origin of the two vectors, the contribution rate of driven-forces were calculated. On the basis of GPS site survey and the maps downloaded from Google Earth, the elevation and shape coefficient were analyzed by Digital Elevation Model (DEM) in Surfer software. Then combined with the regression analysis, the relationship between the vegetation and driven-forces were obtained. [Results] In the floodplain of Lingshan River, there were 62 vegetation species (61 kinds ofAngiospermaesand one kind ofPteridophyta) and the vegetation community structure was mainly divided into 9 classes. The major three driven-forces of the vegetation distribution were elevation factor (ButtElev), shape coefficient (Shapindx) and hydrological characteristics (HydrChar), and their contribution rates were 37.50%, 27.50%, 16.82%, respectively. Particularly, under the influence of ButtElev, from outer to inner edge of floodplain, the moisture capacity varied from strong to weak, the diversity varied from low to high, and the biomass varied from small to large. There was at least 60 percent of the guarantee rate making the vegetation diversity between 1.03 and 1.96 when the SDI coefficient was between 2.0 and 4.0, and the Pe/Pa was between 0.12 and 0.3. The hydrological characteristics affected the morphological plasticity of vegetation and the habitat conditions of the floodplain. The sand content in the outer zone of the floodplain decreased with the increase of the distance from the water, and the anti-scour ability and moisture capacity of the vegetation decreased from strong to weak. [Conclusions] The elevation factor, shape coefficient and hydrological characteristics are the key driven-forces of vegetation distribution in small and medium rivers in hilly area. When we design ecological restoration projects, the three factors should be focused on, and the suitable measures should be employed in order to control and optimize the three factors.

    hilly area; small and medium rivers; floodplain; vegetation distribution; driven-forces

    2016-09-19

    2017-03-23

    項目名稱: 國家自然科學(xué)基金“蜿蜒型河岸帶潛流層水動力學(xué)機制及溶質(zhì)運移規(guī)律研究”(41471069);浙江省水利科技項目“龍游縣中小河流灘地時空演化機理及生態(tài)修復(fù)技術(shù)研究”(RC1527)

    余根聽(1992—),男,碩士研究生。主要研究方向:河岸帶生態(tài)機理。E-mail: yugenting@hhu.edu.cn

    ?通信作者簡介: 夏繼紅(1970—),男,教授,博士生導(dǎo)師。主要研究方向:河岸帶生態(tài)機理。E-mail: syjhxia@hhu.edu.cn

    X171

    A

    2096-2673(2017)02-0051-11

    10.16843/j.sswc.2017.02.007

    猜你喜歡
    灘地植被因子
    浙江省典型河流灘地時空演變分析*
    洞庭湖區(qū)灘地及稻田洼地資源分布研究
    基于植被復(fù)綠技術(shù)的孔植試驗及應(yīng)用
    因子von Neumann代數(shù)上的非線性ξ-Jordan*-三重可導(dǎo)映射
    一些關(guān)于無窮多個素因子的問題
    影響因子
    影響因子
    山丘區(qū)中小河流灘地時空演化與擾動因子分析
    綠色植被在溯溪旅游中的應(yīng)用
    基于形態(tài)指數(shù)的山丘區(qū)中小河流灘地分類方法及演變分析
    国产精品久久电影中文字幕 | 精品国内亚洲2022精品成人 | 久久久久网色| 亚洲综合色网址| 男女之事视频高清在线观看| 亚洲五月色婷婷综合| 99精国产麻豆久久婷婷| 蜜桃国产av成人99| 成人18禁高潮啪啪吃奶动态图| a级毛片在线看网站| 欧美激情 高清一区二区三区| 12—13女人毛片做爰片一| 91大片在线观看| 日韩大码丰满熟妇| 精品人妻熟女毛片av久久网站| 婷婷成人精品国产| 国产成人av激情在线播放| 亚洲熟女精品中文字幕| 国产欧美日韩精品亚洲av| 女人被躁到高潮嗷嗷叫费观| 搡老熟女国产l中国老女人| 亚洲国产中文字幕在线视频| 国产男女超爽视频在线观看| 丝袜喷水一区| 精品少妇一区二区三区视频日本电影| 黑人猛操日本美女一级片| 日本黄色视频三级网站网址 | 1024视频免费在线观看| 精品亚洲成国产av| 欧美人与性动交α欧美软件| 国产日韩一区二区三区精品不卡| 男人操女人黄网站| 国产伦理片在线播放av一区| 99国产极品粉嫩在线观看| 亚洲av成人一区二区三| 成人国语在线视频| 成人手机av| 国产av精品麻豆| 亚洲欧美色中文字幕在线| 亚洲熟女毛片儿| 啦啦啦中文免费视频观看日本| 99riav亚洲国产免费| 老司机影院毛片| 精品视频人人做人人爽| 日韩中文字幕欧美一区二区| 丁香六月欧美| 精品欧美一区二区三区在线| 国产精品久久久久久人妻精品电影 | 欧美精品亚洲一区二区| 中文字幕最新亚洲高清| 国产成人影院久久av| 国产成人精品在线电影| 天天影视国产精品| 在线观看66精品国产| 99精品在免费线老司机午夜| 亚洲第一欧美日韩一区二区三区 | 欧美精品一区二区免费开放| 国产一区二区激情短视频| 亚洲精品在线观看二区| 99九九在线精品视频| 国产亚洲午夜精品一区二区久久| 欧美日韩视频精品一区| 久久中文字幕人妻熟女| 两个人免费观看高清视频| 9热在线视频观看99| 9色porny在线观看| 午夜精品久久久久久毛片777| 五月开心婷婷网| 另类亚洲欧美激情| 搡老岳熟女国产| 另类精品久久| 中文字幕精品免费在线观看视频| 成人18禁在线播放| 亚洲国产av新网站| 国产xxxxx性猛交| 午夜福利视频精品| 精品少妇内射三级| 两个人看的免费小视频| 女人精品久久久久毛片| 国产黄色免费在线视频| 国产又色又爽无遮挡免费看| 久久精品熟女亚洲av麻豆精品| 午夜精品国产一区二区电影| 999精品在线视频| 久久热在线av| 欧美黄色淫秽网站| 超色免费av| 人人妻人人澡人人爽人人夜夜| 精品人妻在线不人妻| 丝袜在线中文字幕| 欧美国产精品va在线观看不卡| svipshipincom国产片| 美女主播在线视频| 18在线观看网站| 露出奶头的视频| kizo精华| tube8黄色片| 久久久久久免费高清国产稀缺| 亚洲一区中文字幕在线| 美女福利国产在线| 中文字幕最新亚洲高清| 成人国产av品久久久| 老鸭窝网址在线观看| 亚洲精品国产精品久久久不卡| 搡老岳熟女国产| 99国产极品粉嫩在线观看| 亚洲中文字幕日韩| 我要看黄色一级片免费的| 50天的宝宝边吃奶边哭怎么回事| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人免费观看mmmm| 成年人黄色毛片网站| 交换朋友夫妻互换小说| 岛国毛片在线播放| 亚洲欧美色中文字幕在线| 成人av一区二区三区在线看| 不卡av一区二区三区| 国产一区二区 视频在线| 国产亚洲欧美精品永久| 亚洲精品自拍成人| 日韩欧美免费精品| 国产亚洲欧美在线一区二区| 日韩欧美三级三区| 国产精品香港三级国产av潘金莲| 国产日韩欧美亚洲二区| 成人特级黄色片久久久久久久 | 久久久精品区二区三区| 久久国产精品影院| 在线 av 中文字幕| 超碰成人久久| 91九色精品人成在线观看| 久久人人爽av亚洲精品天堂| 久久久国产一区二区| 日本wwww免费看| 两性夫妻黄色片| 免费看十八禁软件| 欧美激情极品国产一区二区三区| 久久精品熟女亚洲av麻豆精品| 精品国产亚洲在线| 777久久人妻少妇嫩草av网站| 一级毛片电影观看| 国产成人一区二区三区免费视频网站| 狠狠婷婷综合久久久久久88av| 99国产精品免费福利视频| 色在线成人网| 欧美日韩亚洲高清精品| a级毛片黄视频| 亚洲自偷自拍图片 自拍| 美女福利国产在线| 青青草视频在线视频观看| 国产激情久久老熟女| 亚洲成国产人片在线观看| 午夜精品久久久久久毛片777| 中亚洲国语对白在线视频| 欧美在线一区亚洲| 视频在线观看一区二区三区| 一边摸一边抽搐一进一小说 | 90打野战视频偷拍视频| 免费高清在线观看日韩| 亚洲熟妇熟女久久| 国产91精品成人一区二区三区 | 亚洲五月婷婷丁香| a级毛片在线看网站| 亚洲av第一区精品v没综合| 国产午夜精品久久久久久| 18禁国产床啪视频网站| 久久天堂一区二区三区四区| 亚洲国产看品久久| 无遮挡黄片免费观看| 久久婷婷成人综合色麻豆| 激情在线观看视频在线高清 | 91麻豆精品激情在线观看国产 | 一个人免费在线观看的高清视频| 日韩欧美三级三区| 两个人免费观看高清视频| 日本欧美视频一区| tube8黄色片| 久久精品熟女亚洲av麻豆精品| 日本五十路高清| 亚洲成人免费av在线播放| 男女之事视频高清在线观看| 国产精品.久久久| 国产精品免费一区二区三区在线 | 国产精品一区二区免费欧美| 久久影院123| 一区二区三区国产精品乱码| 亚洲精品中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 男女床上黄色一级片免费看| 人成视频在线观看免费观看| 高清黄色对白视频在线免费看| 欧美日韩中文字幕国产精品一区二区三区 | 国产色视频综合| 欧美乱妇无乱码| 女人精品久久久久毛片| 久久精品国产99精品国产亚洲性色 | 久久婷婷成人综合色麻豆| 精品国产亚洲在线| 欧美亚洲日本最大视频资源| 亚洲全国av大片| 亚洲精品一二三| 一进一出好大好爽视频| 欧美变态另类bdsm刘玥| 久久精品熟女亚洲av麻豆精品| 日韩 欧美 亚洲 中文字幕| 日韩成人在线观看一区二区三区| 一个人免费看片子| 国产一区二区三区综合在线观看| www.999成人在线观看| 天天操日日干夜夜撸| 中亚洲国语对白在线视频| 国产在线免费精品| 国产免费现黄频在线看| 亚洲精品粉嫩美女一区| 欧美在线一区亚洲| 欧美精品一区二区免费开放| 午夜成年电影在线免费观看| 久久99热这里只频精品6学生| 九色亚洲精品在线播放| 婷婷成人精品国产| h视频一区二区三区| 高清黄色对白视频在线免费看| 午夜精品久久久久久毛片777| 国产精品亚洲av一区麻豆| 日日爽夜夜爽网站| 婷婷成人精品国产| 99国产精品免费福利视频| 麻豆av在线久日| 国产精品免费大片| 激情在线观看视频在线高清 | 精品国产国语对白av| 女同久久另类99精品国产91| 亚洲少妇的诱惑av| 欧美午夜高清在线| 天天添夜夜摸| 亚洲人成电影观看| 成年女人毛片免费观看观看9 | 超色免费av| 国产单亲对白刺激| 看免费av毛片| 九色亚洲精品在线播放| 国产野战对白在线观看| 亚洲久久久国产精品| 69精品国产乱码久久久| 国产成人av激情在线播放| 亚洲精品在线美女| 日韩视频在线欧美| 欧美激情极品国产一区二区三区| 免费看十八禁软件| 一本—道久久a久久精品蜜桃钙片| 精品一区二区三区av网在线观看 | 999久久久国产精品视频| 人妻 亚洲 视频| 一夜夜www| 这个男人来自地球电影免费观看| 国产区一区二久久| 脱女人内裤的视频| 久久久久网色| 国产野战对白在线观看| 欧美老熟妇乱子伦牲交| 亚洲av第一区精品v没综合| 免费人妻精品一区二区三区视频| 日韩中文字幕视频在线看片| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品电影小说| 日韩制服丝袜自拍偷拍| 日韩欧美一区视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 69精品国产乱码久久久| 多毛熟女@视频| 精品一区二区三区av网在线观看 | 我要看黄色一级片免费的| 777久久人妻少妇嫩草av网站| a级毛片黄视频| 夜夜骑夜夜射夜夜干| 水蜜桃什么品种好| 丰满饥渴人妻一区二区三| 老熟妇仑乱视频hdxx| 国产在线一区二区三区精| 欧美在线黄色| 精品一区二区三卡| 91麻豆精品激情在线观看国产 | 国产精品秋霞免费鲁丝片| 大片电影免费在线观看免费| 欧美国产精品一级二级三级| 久久国产精品人妻蜜桃| 久久99热这里只频精品6学生| 女人爽到高潮嗷嗷叫在线视频| 激情在线观看视频在线高清 | 香蕉久久夜色| 一二三四在线观看免费中文在| 免费在线观看影片大全网站| 自线自在国产av| 十八禁高潮呻吟视频| 国产不卡av网站在线观看| 亚洲熟妇熟女久久| 国产成人精品久久二区二区91| 制服诱惑二区| 桃花免费在线播放| 男女高潮啪啪啪动态图| 色视频在线一区二区三区| 亚洲自偷自拍图片 自拍| 蜜桃国产av成人99| 美女高潮喷水抽搐中文字幕| 亚洲精华国产精华精| 性少妇av在线| 久久 成人 亚洲| 妹子高潮喷水视频| 在线观看免费视频网站a站| 18禁美女被吸乳视频| 国产国语露脸激情在线看| 精品午夜福利视频在线观看一区 | 无限看片的www在线观看| 老司机靠b影院| 午夜福利乱码中文字幕| 国产男女超爽视频在线观看| 亚洲色图 男人天堂 中文字幕| 久久久国产一区二区| 午夜日韩欧美国产| 国产一区二区三区视频了| a级片在线免费高清观看视频| 国产免费福利视频在线观看| 国产精品久久电影中文字幕 | 国产精品 国内视频| 老司机午夜十八禁免费视频| 欧美大码av| www日本在线高清视频| 亚洲av成人不卡在线观看播放网| 国产伦理片在线播放av一区| 69av精品久久久久久 | 亚洲中文字幕日韩| 十八禁网站免费在线| 又大又爽又粗| 午夜免费鲁丝| 精品亚洲成a人片在线观看| 一二三四社区在线视频社区8| 久久亚洲真实| 成人18禁高潮啪啪吃奶动态图| 不卡一级毛片| 亚洲成人国产一区在线观看| 亚洲一码二码三码区别大吗| 久久久精品国产亚洲av高清涩受| 久久99热这里只频精品6学生| 侵犯人妻中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 精品免费久久久久久久清纯 | 成人手机av| 多毛熟女@视频| 国产精品熟女久久久久浪| 丁香六月天网| 日韩视频一区二区在线观看| 操美女的视频在线观看| 国产亚洲精品一区二区www | 美女国产高潮福利片在线看| 久久久国产欧美日韩av| 国产精品一区二区在线观看99| 欧美日韩福利视频一区二区| 淫妇啪啪啪对白视频| 在线观看66精品国产| 久久狼人影院| 在线看a的网站| 日本vs欧美在线观看视频| 日本av手机在线免费观看| 亚洲欧美精品综合一区二区三区| 搡老岳熟女国产| 亚洲精品在线美女| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美一区二区三区在线观看 | 亚洲视频免费观看视频| 国产av又大| 一本色道久久久久久精品综合| 人人妻人人爽人人添夜夜欢视频| 亚洲美女黄片视频| 丁香欧美五月| 麻豆av在线久日| 国产一区二区 视频在线| 女人被躁到高潮嗷嗷叫费观| av超薄肉色丝袜交足视频| 黄色a级毛片大全视频| 丁香六月天网| 两性夫妻黄色片| av不卡在线播放| 欧美黄色片欧美黄色片| 久久久久国产一级毛片高清牌| 亚洲五月色婷婷综合| 黄色视频不卡| 成人手机av| 一边摸一边抽搐一进一出视频| 免费在线观看完整版高清| 精品久久久久久电影网| 亚洲一区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 侵犯人妻中文字幕一二三四区| 国产精品九九99| www.精华液| 亚洲专区中文字幕在线| 亚洲七黄色美女视频| 国产在线精品亚洲第一网站| 久久久国产欧美日韩av| 国产成人影院久久av| av电影中文网址| 精品少妇黑人巨大在线播放| 欧美一级毛片孕妇| 悠悠久久av| 老鸭窝网址在线观看| 欧美乱码精品一区二区三区| a级片在线免费高清观看视频| 国产高清国产精品国产三级| 99在线人妻在线中文字幕 | 精品国产超薄肉色丝袜足j| 成年人黄色毛片网站| 亚洲精品国产区一区二| 99久久国产精品久久久| 无遮挡黄片免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 精品一品国产午夜福利视频| 国产aⅴ精品一区二区三区波| 一进一出好大好爽视频| 精品午夜福利视频在线观看一区 | 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 十八禁网站免费在线| 精品国产一区二区久久| 曰老女人黄片| 99国产精品免费福利视频| 国产单亲对白刺激| 黄色 视频免费看| 老鸭窝网址在线观看| 丰满少妇做爰视频| 80岁老熟妇乱子伦牲交| 亚洲成a人片在线一区二区| 在线观看www视频免费| 高清黄色对白视频在线免费看| 男男h啪啪无遮挡| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品影院久久| 亚洲天堂av无毛| 国产精品免费视频内射| 欧美日韩亚洲综合一区二区三区_| xxxhd国产人妻xxx| 国产精品久久久人人做人人爽| 国产区一区二久久| 中文字幕制服av| 亚洲一区二区三区欧美精品| 欧美亚洲 丝袜 人妻 在线| 母亲3免费完整高清在线观看| 亚洲精品久久午夜乱码| 人妻 亚洲 视频| 成人国语在线视频| 久久久久久久国产电影| 亚洲人成电影观看| 午夜激情av网站| 中文字幕另类日韩欧美亚洲嫩草| 看免费av毛片| 久久精品国产亚洲av香蕉五月 | 69精品国产乱码久久久| 欧美成人午夜精品| 国产成人精品久久二区二区91| 搡老熟女国产l中国老女人| 成人亚洲精品一区在线观看| 精品国内亚洲2022精品成人 | 免费黄频网站在线观看国产| 国产精品美女特级片免费视频播放器 | 国产高清videossex| 国产精品香港三级国产av潘金莲| av又黄又爽大尺度在线免费看| 波多野结衣av一区二区av| 亚洲成国产人片在线观看| 精品一区二区三卡| 9色porny在线观看| 人人妻人人澡人人爽人人夜夜| 超碰成人久久| 99riav亚洲国产免费| 国产高清视频在线播放一区| 国产成人免费观看mmmm| 精品人妻在线不人妻| 国产在线精品亚洲第一网站| 亚洲性夜色夜夜综合| 国产老妇伦熟女老妇高清| 嫁个100分男人电影在线观看| 亚洲国产成人一精品久久久| 在线av久久热| 欧美老熟妇乱子伦牲交| 精品国产国语对白av| 欧美精品人与动牲交sv欧美| 蜜桃国产av成人99| 久久精品国产99精品国产亚洲性色 | cao死你这个sao货| 999久久久精品免费观看国产| aaaaa片日本免费| 精品国产一区二区三区四区第35| 国产在线精品亚洲第一网站| 极品教师在线免费播放| 九色亚洲精品在线播放| 在线观看免费午夜福利视频| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全电影3 | 99国产精品一区二区蜜桃av | 日本黄色日本黄色录像| 成人永久免费在线观看视频 | 大型黄色视频在线免费观看| 欧美成人免费av一区二区三区 | 精品久久久久久久毛片微露脸| 欧美激情久久久久久爽电影 | 97在线人人人人妻| www.熟女人妻精品国产| 午夜福利免费观看在线| 黄色成人免费大全| 国产欧美日韩综合在线一区二区| 日韩精品免费视频一区二区三区| xxxhd国产人妻xxx| 男女免费视频国产| 欧美亚洲日本最大视频资源| 久久精品国产99精品国产亚洲性色 | 18在线观看网站| 韩国精品一区二区三区| 多毛熟女@视频| 啦啦啦 在线观看视频| av不卡在线播放| 欧美精品一区二区免费开放| 国产日韩欧美亚洲二区| 搡老乐熟女国产| 国产成人精品在线电影| 少妇 在线观看| 日本wwww免费看| 日韩精品免费视频一区二区三区| 丁香六月欧美| 人成视频在线观看免费观看| 亚洲专区中文字幕在线| 夫妻午夜视频| 视频区欧美日本亚洲| 成年版毛片免费区| 亚洲一卡2卡3卡4卡5卡精品中文| 巨乳人妻的诱惑在线观看| 精品少妇久久久久久888优播| 中文字幕人妻丝袜制服| 日本黄色视频三级网站网址 | 无限看片的www在线观看| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美亚洲二区| 国产在线精品亚洲第一网站| 男女之事视频高清在线观看| 国产无遮挡羞羞视频在线观看| 首页视频小说图片口味搜索| 一本一本久久a久久精品综合妖精| 国产成+人综合+亚洲专区| 嫁个100分男人电影在线观看| 国产亚洲精品第一综合不卡| 久久国产精品人妻蜜桃| 黄频高清免费视频| 久久国产精品影院| 久久午夜亚洲精品久久| 久久久久久久大尺度免费视频| 三级毛片av免费| 国产亚洲午夜精品一区二区久久| 欧美黑人精品巨大| 伦理电影免费视频| 国产一区二区三区综合在线观看| 亚洲中文日韩欧美视频| 狠狠婷婷综合久久久久久88av| 欧美激情久久久久久爽电影 | 亚洲国产av影院在线观看| 欧美精品高潮呻吟av久久| 国产激情久久老熟女| 91麻豆精品激情在线观看国产 | 国产熟女午夜一区二区三区| 午夜免费成人在线视频| 老司机福利观看| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 国产在线免费精品| 又黄又粗又硬又大视频| 久久久久久久久久久久大奶| 亚洲三区欧美一区| 制服诱惑二区| 成人av一区二区三区在线看| 国产成人一区二区三区免费视频网站| 日本五十路高清| 国产成人欧美在线观看 | 高清视频免费观看一区二区| 亚洲精品av麻豆狂野| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| a在线观看视频网站| 国产在线一区二区三区精| 国产精品免费大片| 男女床上黄色一级片免费看| 大片免费播放器 马上看| 亚洲七黄色美女视频| 麻豆乱淫一区二区| 久久久国产成人免费| 91成人精品电影| 亚洲avbb在线观看| 亚洲精品一二三| 真人做人爱边吃奶动态| 性色av乱码一区二区三区2| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| 日本五十路高清| 国产精品久久电影中文字幕 | 日韩视频在线欧美| 久久热在线av| 亚洲性夜色夜夜综合| 777久久人妻少妇嫩草av网站| 免费人妻精品一区二区三区视频| 欧美av亚洲av综合av国产av| 国产1区2区3区精品| 黄片大片在线免费观看| 日本欧美视频一区| 久久国产亚洲av麻豆专区| 超色免费av| 国产不卡一卡二| 黑人操中国人逼视频| 欧美日韩视频精品一区|