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

    鄱陽湖濕地苔草(Carex)景觀變化及其水文響應(yīng)?

    2017-07-31 20:03:44周云凱白秀玲寧立新
    湖泊科學(xué) 2017年4期
    關(guān)鍵詞:苔草蘆荻鄱陽湖

    周云凱,白秀玲??,寧立新

    (1:河南大學(xué)資源與環(huán)境研究所,開封475004)(2:河南大學(xué)環(huán)境與規(guī)劃學(xué)院,開封475004)

    鄱陽湖濕地苔草(Carex)景觀變化及其水文響應(yīng)?

    周云凱1,2,白秀玲1,2??,寧立新1,2

    (1:河南大學(xué)資源與環(huán)境研究所,開封475004)
    (2:河南大學(xué)環(huán)境與規(guī)劃學(xué)院,開封475004)

    淡水濕地水文過程控制著濕地植被景觀的形成與演變.基于Landsat TM/ETM+遙感影像數(shù)據(jù),利用決策樹分類法提取鄱陽湖濕地1992、1999、2006、2012年4期景觀信息,通過景觀格局指數(shù)、轉(zhuǎn)移矩陣和質(zhì)心遷移法對(duì)苔草景觀的空間變化及其與水文過程關(guān)系進(jìn)行分析.結(jié)果表明:研究期間,鄱陽湖濕地秋、冬季苔草景觀分布面積受到水位和退水過程的影響,低水位年的中灘位緩慢退水與低灘位快速出露更有利于苔草景觀的擴(kuò)張;苔草景觀的空間格局與水位關(guān)系密切,在低水位年份,低位洲灘提前出露,苔草景觀分布高程較低,部分低位光灘被苔草所侵占,原苔草分布的部分洲灘轉(zhuǎn)變?yōu)樘J荻景觀,景觀的破碎化程度較重,在高水位年份,低位洲灘長(zhǎng)期被水體淹沒,苔草景觀分布高程相對(duì)較高,部分蘆荻分布區(qū)被苔草所侵占,而原苔草景觀的部分區(qū)域轉(zhuǎn)變?yōu)樗w和光灘,由于該期間苔草主要集中分布在湖周和入湖河口地帶的高位洲灘上,其景觀破碎化程度較輕;水位年際間的升降變化會(huì)影響苔草景觀質(zhì)心位置,年均水位上升引起景觀質(zhì)心發(fā)生向湖岸方向推進(jìn),而年均水位下降則會(huì)導(dǎo)致苔草景觀質(zhì)心向湖心方向轉(zhuǎn)移.

    苔草景觀;景觀格局;水位變化;響應(yīng);鄱陽湖濕地

    濕地景觀格局是各種生態(tài)過程在不同尺度上綜合作用的結(jié)果[1].其中,水文過程是淡水濕地景觀格局形成與演變更為重要的控制性因素,主導(dǎo)了濕地植被基本的空間分布格局[2-3].不同的水文過程導(dǎo)致土壤水分梯度變化,引起濕地植物的分布范圍和群落結(jié)構(gòu)出現(xiàn)空間差異,形成不同的植被景觀類型,并在空間上表現(xiàn)為一定的景觀格局[4].在淡水濕地生態(tài)系統(tǒng)中,植被景觀往往表現(xiàn)出沿水分梯度呈現(xiàn)較為明顯的帶狀分布格局[5-6].而水文過程年際間的變動(dòng),又會(huì)影響濕地植物群落的穩(wěn)定性和物種多樣性,當(dāng)變動(dòng)超過植物的耐受范圍,則會(huì)引起植物群落發(fā)生正向或逆向演替,從而導(dǎo)致植被景觀格局發(fā)生改變[7-10].

    鄱陽湖是我國(guó)第一大淡水湖,也是長(zhǎng)江中游重要的通江湖泊之一.受長(zhǎng)江和“五河”(贛江、撫河、信江、饒河、修水)來水的雙重影響,湖泊水位年內(nèi)變化巨大,呈現(xiàn)豐水期與枯水期交替出現(xiàn)的獨(dú)特水文節(jié)律[11],這種節(jié)律性水位變化,使洲灘各高程帶的出露—淹沒過程不一,形成不同的濕地水文環(huán)境和與之相適應(yīng)的植物群落帶.此外,由于水位年際間的波動(dòng)性和多變性,使?jié)竦卦兴沫h(huán)境遭到破壞,引起植物群落結(jié)構(gòu)、分布范圍及分布高程發(fā)生變化[12].尤其近十幾年來,受氣候因素和三峽工程運(yùn)行的影響,江湖關(guān)系發(fā)生顯著變化[13],導(dǎo)致鄱陽湖水位較常年偏低,枯水期提前、延長(zhǎng),植被景觀的空間格局隨之發(fā)生明顯改變[14].針對(duì)鄱陽湖濕地植被,早期更多地基于野外實(shí)地調(diào)查和遙感解譯,研究鄱陽湖濕地植被類型、分布及其演化規(guī)律[15-18],近年來,許多學(xué)者開始關(guān)注鄱陽湖濕地植被演化與水文過程之間的關(guān)系,并取得了一定的研究進(jìn)展[19-22],但研究的廣度和深度還有待進(jìn)一步推進(jìn).本文基于景觀生態(tài)學(xué)角度,利用遙感影像數(shù)據(jù)研究鄱陽湖濕地典型植被——苔草(Carex)景觀的空間格局變化及其對(duì)水文變化的響應(yīng),研究結(jié)果對(duì)于進(jìn)一步深入理解濕地植被景觀格局演變具有一定的理論價(jià)值,同時(shí)對(duì)于鄱陽湖濕地生態(tài)環(huán)境保護(hù)也具有重要的現(xiàn)實(shí)意義.

    1 材料與方法

    1.1 研究區(qū)域

    鄱陽湖濕地是我國(guó)最大的淡水湖濕地,也是世界重要濕地之一,地處江西省北部,長(zhǎng)江中下游南岸,水域面積在夏季豐水期達(dá)到3700 km2以上,而在秋、冬季枯水期縮減至1000 km2以下,呈現(xiàn)“洪水一片、枯水一線”的獨(dú)特自然景觀[20].水位的周期性漲落形成了大面積水陸交替的草洲、泥灘和沙灘.在全湖尺度上,洲灘植被沿水分梯度總體呈較為明顯的條帶狀分布格局,由湖岸向湖心方向分布的主要優(yōu)勢(shì)植物群落為蘆葦群落(Phragmites)、南荻群落(Triarrhena)、苔草群落(Carex)、虉草群落(Phalaris)和沉水植物群落.蘆葦與南荻群落多分布于16~18m高程的高位洲灘,每年3月前后開始萌芽,9-10月開花,11-12月開始枯萎;苔草與虉草群落多分布于12~16m高程的中低位洲灘[15],其中,苔草群落是鄱陽湖濕地分布最廣、面積最大的植被群落類型[23-24],分為春草和秋草,春草一般2-5月為萌發(fā)生長(zhǎng)期,豐水期被湖水淹沒,大量死亡或進(jìn)入休眠狀態(tài),秋草則在秋季退水期洲灘出露后開始生長(zhǎng),一般9-12月為萌發(fā)生長(zhǎng)期,次年1-2月枯萎[23];虉草每年2月萌芽生長(zhǎng),至豐水期基本完成其生命周期,秋季洲灘出露后,只有少數(shù)虉草還可繼續(xù)生長(zhǎng);沉水植物群落分布在12m高程以下且長(zhǎng)期處于淹水狀態(tài)的低洼地帶[15].在中等尺度下,由于受微地形和水分條件的影響,各植被類型在灘地上交錯(cuò)出現(xiàn),構(gòu)成濕地群落復(fù)合體結(jié)構(gòu)和水平鑲嵌結(jié)構(gòu)[23].

    依據(jù)星子站水位數(shù)據(jù),1992-2012年間鄱陽湖濕地多年月均水位最高值在18 m左右,明顯高于苔草群落的分布高程,因此,本研究以鄱陽湖區(qū)DEM生成的18m等高線和湖區(qū)圩堤界線為基礎(chǔ),提取與鄱陽湖連成一體并受鄱陽湖洪水影響的區(qū)域,作為鄱陽湖濕地范圍(圖1).

    1.2 數(shù)據(jù)來源與處理

    研究使用的數(shù)據(jù)包括遙感數(shù)據(jù)和非遙感數(shù)據(jù)2種.由于鄱陽湖水位年內(nèi)變化劇烈,在高水位時(shí),大部分洲灘植被為湖水所淹沒,而低水位時(shí),洲灘植被出露水面,為避免水淹的干擾,在遙感影像選取時(shí),盡量選擇水位較低且苔草等植被生長(zhǎng)較為旺盛的秋、冬季影像,該時(shí)段鄱陽湖已基本進(jìn)入枯水期,苔草等植被幾乎全部出露水面,且此時(shí)苔草植被正值成熟期,與其他植被影像特征差異較為明顯,另外,在秋、冬季鄱陽湖區(qū)云雨較少,影像質(zhì)量較高,便于進(jìn)行遙感影像解譯,因此,遙感數(shù)據(jù)最終選取1992年11月28日、1999年12月10日、2006年11月3日和2012年10月26日的Landsat TM/ETM+影像,其軌道號(hào)為121/40,空間分辨率為30 m.非遙感數(shù)據(jù)包括鄱陽湖區(qū)1∶50000地形圖和DEM、圩堤數(shù)據(jù)以及野外考察采集的點(diǎn)位數(shù)據(jù),作為研究區(qū)域確定和景觀信息提取的輔助數(shù)據(jù);另外,還包括鄱陽湖水位數(shù)據(jù),鄱陽湖共設(shè)有康山、吳城、都昌、星子、湖口5個(gè)基本水文站,其中星子水文站位于鄱陽湖入江水道上段左岸,能夠較為全面地反映鄱陽湖在江湖作用層面以及大湖面的水情變化[11],故選星子水文站觀測(cè)水位(吳淞高程,下同)代表鄱陽湖水位.

    遙感影像處理首先借助ENVI 4.7軟件對(duì)遙感影像進(jìn)行波段合成;以1∶50000地形圖作為參考圖像對(duì)影像進(jìn)行幾何精校正;利用確定的鄱陽湖濕地范圍作為掩膜對(duì)合成后影像進(jìn)行裁剪;依據(jù)影像的光譜特征及計(jì)算得到的NDVI值,采用決策樹分類法對(duì)秋、冬季枯水期鄱陽湖濕地遙感影像進(jìn)行解譯.鄱陽湖濕地景觀可分為非植被景觀和植被景觀.非植被景觀包括水體、泥灘和沙灘;植被景觀主要包括沉水植物群落、中低位洲灘的苔草與虉草群落、高位洲灘的蘆葦與南荻群落.其中,沉水植被沒于湖中,在TM/ETM+遙感影像中難以與水體區(qū)分,并與水體區(qū)域疊合,故分類時(shí)未單獨(dú)劃分,而將其歸入到水體景觀中;泥灘與沙灘合并為光灘景觀;苔草、南荻、蘆葦?shù)戎脖皇芹蛾柡竦厍铩⒍咀钪饕膬?yōu)勢(shì)植被類型,呈不規(guī)則的帶狀或片狀分布,面積較大,其他濕生植被(如虉草等)在秋、冬季分布面積較小,且常與優(yōu)勢(shì)群落相混雜[18].另外,由于TM/ETM+影像的分辨率為30 m,對(duì)于面積較小的植被類型難以準(zhǔn)確解譯出來.基于上述考慮,根據(jù)各景觀的影像特征和植被的生長(zhǎng)發(fā)育節(jié)律不同,本研究將鄱陽湖濕地景觀分為水體景觀(包括水體與沉水植物群落)、光灘景觀(包括泥灘和沙灘)、苔草景觀(以苔草為主,包括虉草、藜蒿在內(nèi)的中低位洲灘植物群落)和蘆荻景觀(以蘆葦與南荻為主的高位洲灘植物群落)4種類型;利用野外采集的點(diǎn)位數(shù)據(jù)和目視判讀數(shù)據(jù)進(jìn)行解譯精度評(píng)價(jià),并通過人工目視糾正錯(cuò)分區(qū)域,確保遙感影像解譯精度達(dá)到80%以上;最終形成鄱陽湖濕地景觀類型圖(圖2).

    圖1 鄱陽湖濕地范圍Fig.1 Location of the Lake Poyang Wetland

    1.3 研究方法

    利用景觀格局指數(shù)法研究苔草景觀結(jié)構(gòu)變化.基于斑塊類型水平上選取景觀類型面積、景觀類型面積百分比、平均斑塊面積、斑塊數(shù)目、斑塊密度、平均形狀指數(shù)和平均分維數(shù)等指標(biāo),利用ArcGIS和FRAGSTATS 4.0軟件提取各景觀格局指數(shù),從分布面積、破碎化、形狀方面研究苔草景觀結(jié)構(gòu)變化.各指標(biāo)計(jì)算方法及生態(tài)意義參考文獻(xiàn)[25]和Fragstats用戶指南.

    利用轉(zhuǎn)移矩陣法研究苔草景觀與其他景觀類型間的相互轉(zhuǎn)換.通過ARCGIS空間分析模塊對(duì)1992、1999、2006和2012年4期遙感影像解譯圖進(jìn)行空間疊加分析,屬性數(shù)據(jù)整理后得到1992-1999、1999-2006和2006-2012年3個(gè)時(shí)段苔草景觀與其他景觀類型相互轉(zhuǎn)化的面積.轉(zhuǎn)移矩陣的表達(dá)式[26]為:

    式中,S為面積,n為轉(zhuǎn)移前后景觀類型數(shù),i、j(i,j=1、2、…、n)分別為轉(zhuǎn)移前、轉(zhuǎn)移后的景觀類型,Sij為前一期的i類景觀轉(zhuǎn)變成后一期的j類景觀的面積.矩陣中每一行元素代表前期i類景觀向后期各類景觀的流出信息,每一列元素代表后期的j類景觀從前期的各類景觀的流入信息.

    圖2 1992、1999、2006和2012年鄱陽湖濕地景觀類型Fig.2 Landscape types of Lake PoyangWetland in 1992,1999,2006 and 2012

    利用分布質(zhì)心遷移法研究苔草景觀空間分布的變動(dòng)情況.借助ARCGIS的空間分析工具生成1992、1999、2006和2012年4個(gè)時(shí)期苔草各斑塊的分布質(zhì)心,并生成各斑塊質(zhì)心的經(jīng)緯度,運(yùn)用公式(2)計(jì)算得到各期苔草景觀分布質(zhì)心,通過對(duì)景觀分布質(zhì)心的比較,反映不同時(shí)期鄱陽湖濕地苔草景觀在空間上的變化特征.質(zhì)心坐標(biāo)的計(jì)算公式[27]為:

    式中,X和Y分別為按面積加權(quán)得到的質(zhì)心經(jīng)度和緯度,Xi和Yi分別為苔草景觀第i個(gè)斑塊的經(jīng)度和緯度,Ci為苔草景觀第i個(gè)斑塊的面積,n為苔草景觀的斑塊總數(shù)目.

    2 結(jié)果與分析

    2.1 苔草景觀格局變化

    2.1.1 苔草景觀結(jié)構(gòu)變化 研究期間,鄱陽湖濕地苔草景觀面積經(jīng)歷了先逐漸增大后減小的變化過程(表1),最大分布面積為2006年的1019.7 km2,占研究區(qū)面積的33.73%,與1992年相比,苔草景觀面積增加了423.22 km2,擴(kuò)張近71%;2012年苔草景觀面積為737.72 km2,較2006年減少了281.98 km2.斑塊數(shù)、斑塊密度與面積的變化趨勢(shì)一致,而平均斑塊面積表現(xiàn)為先逐漸減小后增大的趨勢(shì),這表明苔草景觀的破碎化程度也經(jīng)歷了先增大后減小的變化過程.從景觀形狀來看,4期的平均分維數(shù)和平均形狀指數(shù)無明顯變化趨勢(shì),基本保持穩(wěn)定,說明苔草景觀在自然狀態(tài)下形狀的不規(guī)則性變化較小.

    表1 1992-2012年苔草景觀格局指數(shù)Tab.1 Changes of class-level indexes of Carex landscape during 1992-2012

    2.1.2 苔草景觀轉(zhuǎn)移變化 苔草景觀除了面積發(fā)生變化外,其分布區(qū)域也在研究期間發(fā)生明顯改變,這意味著苔草景觀與其他景觀類型之間進(jìn)行著較為頻繁的轉(zhuǎn)入與轉(zhuǎn)出過程.1992-1999年間,苔草景觀主要與蘆荻以及光灘景觀發(fā)生相互轉(zhuǎn)換,而與水體間的轉(zhuǎn)換極少(表2).該階段苔草景觀轉(zhuǎn)出的主要方向?yàn)樘J荻景觀,有114.94 km2的苔草景觀轉(zhuǎn)變?yōu)樘J荻景觀,占轉(zhuǎn)出總面積的70%,轉(zhuǎn)為光灘景觀部分的面積占轉(zhuǎn)出總面積的27.32%;而轉(zhuǎn)入部分也主要來自蘆荻景觀,共有142.82 km2的蘆荻景觀轉(zhuǎn)變?yōu)樘Σ菥坝^,占轉(zhuǎn)入量的52.96%,而光灘景觀有124.85 km2轉(zhuǎn)變?yōu)樘Σ菥坝^,占轉(zhuǎn)入量的46.30%.在1999-2006年間,苔草景觀面積進(jìn)一步擴(kuò)大,達(dá)到1019.70 km2,該階段苔草景觀轉(zhuǎn)出方向仍以蘆荻景觀為主,占轉(zhuǎn)出總面積的76.92%,轉(zhuǎn)為光灘景觀部分占轉(zhuǎn)出總面積的15.67%;而轉(zhuǎn)入部分則與前一階段不同,光灘成為該階段的主要轉(zhuǎn)入景觀類型,占轉(zhuǎn)入總面積的57.66%,而蘆荻景觀的轉(zhuǎn)入部分占轉(zhuǎn)入總面積的28.32%,此外,也有66.17 km2的水體景觀轉(zhuǎn)變?yōu)樘Σ菥坝^,占轉(zhuǎn)入總面積的14.02%,明顯高于上一階段.在2006-2012年間,苔草景觀面積大幅減少,有506.39 km2的苔草景觀轉(zhuǎn)出為其他景觀類型,其中有297.31 km2轉(zhuǎn)變?yōu)樗w景觀,占轉(zhuǎn)出總面積的58.71%,是該階段的主要轉(zhuǎn)出方向,其次是光灘,轉(zhuǎn)移面積為126.14 km2,占24.91%;與此同時(shí),也有部分其他景觀轉(zhuǎn)入為苔草景觀,最主要的轉(zhuǎn)入來源為蘆荻景觀,轉(zhuǎn)入面積為155.22 km2,占轉(zhuǎn)入總面積的67.17%,另外,還有29.46%來自于光灘景觀,而水體景觀轉(zhuǎn)入部分極少,僅占轉(zhuǎn)入部分的1.37%.

    表2 3個(gè)時(shí)段苔草與其他景觀類型之間的轉(zhuǎn)移面積(km2)Tab.2 Transfer among Carex and other landscape types in three different periods

    上述分析可以看出,研究期間苔草景觀與其他景觀類型之間進(jìn)行著頻繁的相互轉(zhuǎn)換,但不同時(shí)段其轉(zhuǎn)移方向及轉(zhuǎn)移量存在較大差異.1992-1999和1999-2006年期間,苔草景觀主要轉(zhuǎn)出為蘆荻景觀,2006-2012年間則主要轉(zhuǎn)出為水體景觀;而3個(gè)時(shí)段的轉(zhuǎn)入部分主要來自蘆荻和光灘景觀,占總轉(zhuǎn)入面積的86%以上.

    2.1.3 苔草景觀質(zhì)心遷移變化 研究期間,苔草景觀分布質(zhì)心均位于大湖區(qū)的西南部,但不同時(shí)期,分布質(zhì)心的遷移方向有所不同.在1992-1999年間,鄱陽湖濕地苔草景觀分布質(zhì)心向西北方向偏移了739.87 m,整體向贛江中支入湖河口形成的三角洲方向推進(jìn),質(zhì)心點(diǎn)高程由1992年的12.94m略升高到1999年的13.01 m;1999-2006年間,苔草景觀分布質(zhì)心又向東偏南方向偏移了1527.19 m,年均偏移量(218.17 m/a)明顯大于前一階段(105.70m/a),該時(shí)期質(zhì)心點(diǎn)向遠(yuǎn)離湖岸方向推進(jìn),質(zhì)心點(diǎn)高程明顯降低,降低了0.36 m;而在2006-2012年間,分布質(zhì)心向南遷移477.38 m,年均偏移量(79.56 m/a)最小,質(zhì)心點(diǎn)高程升高了0.12 m,這一階段分布質(zhì)心略表現(xiàn)出向湖岸方向推進(jìn)(圖3).

    圖3 1992-2012年苔草景觀分布質(zhì)心遷移變化Fig.3 Centriods shifting of Carex landscape during 1992-2012

    2.2 鄱陽湖水位變化

    1992-2012年間,鄱陽湖年均水位波動(dòng)頻繁(圖4),多年平均水位為13.20m,最大值出現(xiàn)在1998年,為15.61 m,最小值出現(xiàn)在2011年,為10.96 m,年均水位最大差值為4.65 m.從多年變化曲線來看,鄱陽湖年均水位整體表現(xiàn)出一定的下降趨勢(shì).2006年為相對(duì)枯水年,年均水位為11.55 m,低于多年平均值1.65 m;1999年為相對(duì)的豐水年,年均水位達(dá)到14.12 m,高于多年平均值0.92m;而1992和2012年的年均水位也高于多年平均值,分別為13.33和13.74 m.

    圖4 1992-2012年鄱陽湖星子站年均水位變化Fig.4 Annualwater level change of Xingzi Station of Lake PoyangWetland during 1992-2012

    由于受入湖河流及長(zhǎng)江徑流過程的影響,不同年份鄱陽湖水文變化過程有所差異.研究期間,水位年內(nèi)變化雖然均呈單峰模式(圖5),但2006年水位變化過程最為平緩,最高水位為16.72m,最低水位為7.84m,年內(nèi)變幅為8.88m;而1999年水位變化過程最為劇烈,最高水位達(dá)到21.71 m,最低水位為7.83 m,年內(nèi)變幅達(dá)到13.88 m,1992和2012年的水位變化過程則介于前兩者之間,年內(nèi)變幅分別為12.41和11.84m.若以年內(nèi)最高水位為界將鄱陽湖水位過程分為漲水和退水2個(gè)階段,則2006年退水時(shí)間相對(duì)提前(6月21日),退水過程最為緩慢,日均水位降低0.05 m;2012年退水時(shí)間相對(duì)滯后(8月13日),退水過程相對(duì)較快,日均水位降低0.07 m;而1999年自7月21日開始退水,退水階段水位日均降低值最大,為0.08m/d.

    圖5 1992、1999、2006和2012年鄱陽湖水位變化過程Fig.5Water level changes of Lake Poyang in 1992,1999,2006 and 2012

    表3 不同高程帶苔草景觀面積百分比構(gòu)成Tab.3 The area percentages of Carex landscape in different elevation zones

    2.3 苔草景觀對(duì)水文變化的響應(yīng)

    2.3.1 苔草景觀分布與結(jié)構(gòu)對(duì)水文變化的響應(yīng)鄱陽湖濕地苔草景觀的空間分布主要受地貌及水文條件影響.對(duì)于秋季苔草而言,地貌高程與濕地退水過程的影響更為重要.從苔草的分布高程來看(表3),研究的4個(gè)年份中,苔草均分布在16 m高程以下,其中12~14 m高程帶是苔草最為主要的分布區(qū),但不同年份各高程帶上苔草面積占比不同,反映了苔草主要分布高程的變化.在枯水年(2006年),12 m以下高程的苔草面積占比高于其他年份,而14~16 m高程的面積占比低于其他高程,說明枯水年的低水位引起苔草分布高程出現(xiàn)一定程度的下移.

    表4 退水過程不同階段經(jīng)歷時(shí)間Tab.4 Duration of different stages in the process ofwater recession

    從苔草景觀分布面積來看(表1),枯水年(2006年)苔草面積遠(yuǎn)高于其他年份,但最小面積并未出現(xiàn)在豐水年(1999年),而是出現(xiàn)在相對(duì)的平水年(1992年),這說明水位僅是影響苔草景觀分布的重要影響因子之一,苔草景觀的分布還受到其他水文因子的影響.通過對(duì)4個(gè)年份退水過程分析發(fā)現(xiàn),不同年份退水過程各階段歷時(shí)差異明顯(表4),在苔草分布的主要高程區(qū)間,水位由16m降到14m,以苔草分布面積最大的2006年歷時(shí)最長(zhǎng)(36 d),其次為2012年(29 d),而苔草分布面積最小的1992年歷時(shí)最短(14 d);當(dāng)水位由14m繼續(xù)下降到12m,2006年歷時(shí)最短(8 d),其次為2012年(18 d),而1992年歷時(shí)最長(zhǎng)(41 d).結(jié)合4個(gè)年份的苔草景觀面積可發(fā)現(xiàn),中灘位(14~16m高程帶)的退水歷時(shí)長(zhǎng)且低灘位(12~14m高程帶)退水歷時(shí)短,苔草景觀面積則大,反之,則小.據(jù)此可初步推測(cè),中灘位的緩慢退水和低灘位的快速出露聯(lián)合作用下的水文經(jīng)歷更有利于苔草的萌芽生長(zhǎng)與擴(kuò)張,形成面積較大的苔草景觀;相反,中灘位的快速退水和低灘位的緩慢出露聯(lián)合作用下的水文經(jīng)歷則會(huì)對(duì)苔草的萌發(fā)與分布產(chǎn)生抑制作用,這一定程度上反映退水過程對(duì)苔草景觀分布具有重要影響.

    與此同時(shí),隨著水位的年際間升降,苔草景觀空間結(jié)構(gòu)也發(fā)生相應(yīng)變化.在低水位年份(2006年),鄱陽湖濕地的洲灘大面積出露,尤其在“五河”河口及湖心地帶,由河流沖淤作用形成諸多大小不一的洲灘出露水面,在一定的高程帶上形成了苔草適生的水文環(huán)境,苔草開始萌芽、生長(zhǎng),苔草空間上呈較為明顯的片狀分布格局,苔草斑塊數(shù)目較多,而平均斑塊面積較小,苔草景觀的破碎化程度相對(duì)較重;而在高水位年份(1999年),由于河口與湖中部分洲灘被水體所淹沒,苔草的分布區(qū)域主要集中在湖周及河口地帶高程相對(duì)較高的洲灘上,苔草景觀空間上呈現(xiàn)出條帶狀分布格局,苔草分布的斑塊數(shù)目較少,平均斑塊面積較大,景觀破碎化程度相對(duì)較輕.而斑塊平均形狀指數(shù)與平均分維數(shù)對(duì)水位變化并不敏感,苔草斑塊始終在自然狀態(tài)下呈現(xiàn)出不規(guī)則狀.

    2.3.2 苔草景觀轉(zhuǎn)移對(duì)水文變化的響應(yīng) 由于水位年際間的變動(dòng),苔草分布區(qū)域及分布高程發(fā)生改變,導(dǎo)致苔草景觀與其他景觀類型間不斷進(jìn)行轉(zhuǎn)換.當(dāng)處于相對(duì)的高水位年份(如1999、2012年),苔草分布高程有所上移,侵占先前的蘆荻景觀分布區(qū),然由于濕地洲灘地貌的復(fù)雜性以及水文情勢(shì)的多變性,部分原苔草生長(zhǎng)區(qū)域形成更有利于蘆荻生長(zhǎng)的水文環(huán)境,使苔草景觀被蘆荻景觀所取代;此外,還有部分原苔草分布的低位洲灘由于長(zhǎng)期的水淹而不利于苔草自身的生存,轉(zhuǎn)變成為光灘和水體景觀;而處于相對(duì)的低水位年份(2006年),苔草的分布高程有所下移,部分光灘被苔草所侵占,原苔草分布的高位灘地則轉(zhuǎn)變?yōu)樘J荻景觀,此外,受局地特殊的水文與地貌條件影響,也有部分蘆荻景觀轉(zhuǎn)變?yōu)樘Σ菥坝^.

    2.3.3 苔草景觀質(zhì)心對(duì)水文變化的響應(yīng) 隨著苔草分布區(qū)域與分布高程的變化,苔草景觀的分布質(zhì)心也不斷發(fā)生遷移.從景觀質(zhì)心遷移方向與水位變化情況可看出,苔草景觀質(zhì)心遷移方向與水位升降緊密相關(guān).在水位上升時(shí)期(1992-1999、2006-2012年),苔草的主要分布區(qū)域向湖岸帶和入湖河流河口區(qū)的較高灘位轉(zhuǎn)移,其景觀質(zhì)心會(huì)發(fā)生向湖岸方向的遷移;而在水位降低時(shí)期(1999-2006年),苔草分布高程下降,湖中出露的灘地、入湖河口以及湖岸帶的中低灘位被苔草侵占,其主要分布區(qū)域向湖心低位洲灘方向擴(kuò)張,從而引起分布質(zhì)心也隨之向湖心方向遷移.

    3 結(jié)語

    利用1992、1999、2006和2012年4期Landsat TM/ETM+遙感影像數(shù)據(jù),從景觀生態(tài)學(xué)角度出發(fā),分析了鄱陽湖濕地典型植被—苔草景觀的空間變化,并對(duì)其與水文過程變化之間的關(guān)系進(jìn)行初步探討,主要得到以下結(jié)論:

    1)1992-2012年間,鄱陽湖濕地年均水位波動(dòng)頻繁,然整體略呈下降趨勢(shì).不同年份濕地水文變化過程差異明顯,其中2006年水位變化最為平緩,退水過程最慢;而1999年水位變化最為劇烈,退水過程最快.

    2)鄱陽湖濕地秋、冬季苔草景觀分布主要受地貌高程和水文過程影響.苔草景觀分布在16 m高程以下,而12~14 m高程帶是其最主要的分布區(qū)間;苔草景觀分布面積受到水位的影響,但退水過程影響更為重要,尤其是中灘位(14~16m高程帶)與低灘位(12~14m高程帶)退水過程的聯(lián)合作用對(duì)苔草景觀的分布起到促進(jìn)或抑制作用.

    3)苔草景觀空間格局與水位密切相關(guān),在低水位年份,苔草景觀的破碎化程度較重,景觀質(zhì)心向湖心方向遷移;在高水位年份,其景觀破碎化程度較輕,景觀質(zhì)心會(huì)向湖岸方向推進(jìn).

    鄱陽湖濕地植被景觀格局的變化受多種因素影響,雖然水文過程起到主要的控制作用,但其他因素(如泥沙沖淤、土壤環(huán)境、人為開墾、火燒、放牧等)的影響也不容忽視.目前,淡水濕地植被景觀格局與水文變化關(guān)系研究雖取得一定進(jìn)展,但研究的深度還十分有限,未來還需更為長(zhǎng)期的系統(tǒng)觀測(cè)與模擬試驗(yàn)才能定量地揭示兩者之間的作用機(jī)理;同時(shí),也需要加強(qiáng)其他因素在濕地植被景觀格局變化中的作用研究,從而進(jìn)一步豐富和完善濕地植被格局演變理論.

    致謝:感謝中國(guó)科學(xué)院鄱陽湖湖泊濕地觀測(cè)研究站在野外點(diǎn)位數(shù)據(jù)采集時(shí)提供的幫助;同時(shí)感謝中國(guó)科學(xué)院南京地理與湖泊研究所陳宇煒研究員,孫占東、賴錫軍副研究員提供的水位數(shù)據(jù).

    [1] Bai Junhong,Ouyang Hua,Yang Zhifeng et al.Changes in wetland landscape patterns:A review.Progress in Geography,2005,24(4):36-45.[白軍紅,歐陽華,楊志峰等.濕地景觀格局變化研究進(jìn)展.地理科學(xué)進(jìn)展,2005,24(4):36-45.]

    [2] Huang Yilong,F(xiàn)u Bojie,Chen Liding.Advances in ecohydrological process research.Acta Ecologica Sinica,2003,23(3):580-587.[黃奕龍,傅伯杰,陳利頂.生態(tài)水文過程研究進(jìn)展.生態(tài)學(xué)報(bào),2003,23(3):580-587.]

    [3] Li Xu,Xie Yonghong,Huang Jishan et al.Research progresses on the formationmechanism of vegetation distribution pattern in wetland.Wetland Science,2009,7(3):280-288.[李旭,謝永宏,黃繼山等.濕地植被格局成因研究進(jìn)展.濕地科學(xué),2009,7(3):280-288.]

    [4] Cao Qing,Huang Senkai,Huang Binbin et al.Research advances in lakeside zone landscape pattern and hydrological process.Pearl River,2014,(3):82-86.[曹慶,黃森開,黃彬彬等.水文過程對(duì)湖濱帶景觀格局影響研究進(jìn)展.人民珠江,2014,(3):82-86.]

    [5] Bornman TG,Adams JB,Bate GC.Environmental factors controlling the vegetation zonation patterns and distribution of vegetation types in the Olifants Estuary,South Africa.South African Journal ofBotany,2008,74(4):685-695.

    [6] Huang Qun,Jiang Jiahu,Lai Xijun etal.Changes of landscape structure in Dongting Lakewetlandsand the evaluation on impacts from operation of the Three Gorges project.Resources and Environment in the Yangtze Basin,2013,22(7):923-927.[黃群,姜加虎,賴錫軍等.洞庭湖濕地景觀格局變化以及三峽工程蓄水對(duì)其影響.長(zhǎng)江流域資源與環(huán)境,2013,22(7):923-927.]

    [7] Bacon PE,Stone C,Binns DL et al.Relationships between water availability and Eucalyptus camaldtdensis growth in a riparian forest.Journal ofHydrology,1993,150:541-561.

    [8] Beauchamp VB,Stromberg JC.Flow regulation of the Verde River,Arizona encourages Tamarix recruitment buthasminimal effect on Populus and Salix stand density.Wetlands,2007,27(2):381-389.

    [9] Yao Xin,Yang Guishan,Wan Rongrong et al.Impact of water level change on wetland vegetation of rivers and lakes.J Lake Sci,2014,26(6):813-821.DOI:10.18307/2014.0601.[姚鑫,楊桂山,萬榮榮等.水位變化對(duì)河流、湖泊濕地植被的影響.湖泊科學(xué),2014,26(6):813-821.]

    [10] Li Shengnan,Wang Genxu,DengWei.Research advances in wetland landscape pattern and hydrological process.Chinese Journal ofEcology,2008,27(6):1012-1020.[李勝男,王根緒,鄧偉.濕地景觀格局與水文過程研究進(jìn)展.生態(tài)學(xué)雜志,2008,27(6):1012-1020.]

    [11] DaiXue,Wan Rongrong,Yang Guishan etal.Temporal variation of hydrological rhythm in Poyang Lake and the associated water exchange with the Changjiang River.Scientia Geographica Sinica,2014,24(12):1488-1496.[戴雪,萬榮榮,楊桂山等.鄱陽湖水文節(jié)律變化及其與江湖水量交換的關(guān)系.地理科學(xué),2014,24(12):1488-1496.]

    [12] Yu Li,He Longhua,Zhang Qi etal.Landsat-TM data based study on dynamic changes of typicalwetlands of Poyang Lake. Remote Sensing Information,2010,(6):48-52.[余莉,何隆華,張奇等.基于Landsat-TM影像的鄱陽湖典型濕地動(dòng)態(tài)變化研究.遙感信息,2010,(6):48-52.]

    [13] Dai X,Wan RR,Yang GS.Non-stationarywater-level fluctuation in China's Poyang Lake and its interactionswith Yangtze River.Journal ofGeographical Sciences,2015,25(3):274-288.

    [14] Hu Zhenpeng,Ge Gang,Liu Chenglin.Cause analysis and early warning for wetland vegetation degradation in Poyang Lake.Resources and Environment in the Yangtze Basin,2015,24(3):381-386.[胡振鵬,葛剛,劉成林.鄱陽湖濕地植被退化原因分析及其預(yù)警.長(zhǎng)江流域資源與環(huán)境,2015,24(3):381-386.]

    [15] Zhang Quanjun,Yu Xiubo,Hu Binhua.Research on the characteristics of plant communities in the Poyang NanjiWetlands,China.Resources Science,2013,35(1):42-49.[張全軍,于秀波,胡斌華.鄱陽湖南磯濕地植物群落分布特征研究.資源科學(xué),2013,35(1):42-49.]

    [16] Liu Xiaoli,Ding Mingjun,Li Guicai et al.Research on variation ofwetland plant communities along elevation gradient around Poyang Lake.Yangtze River,2013,44(5):82-86.[劉肖利,丁明軍,李貴才等.鄱陽湖濕地植物群落沿高程梯度變化特征研究.人民長(zhǎng)江,2013,44(5):82-86.]

    [17] Ge Gang,Zhao Anna,Zhong Yiyong et al.Patterns of dominant populations of plants in Islets of Poyang Lake.WetlandScience,2011,9(1):19-25.[葛剛,趙安娜,鐘義勇等.鄱陽湖洲灘優(yōu)勢(shì)植物種群的分布格局.濕地科學(xué),2011,9(1):19-25.]

    [18] Wu Jiandong,Liu Guanhua,Jin Jiefeng et al.Structure analysis of beach vegetation in Poyang Lake in autumn.Jiangxi Science,2010,28(4):549-554.[吳建東,劉觀華,金杰峰等.鄱陽湖秋季洲灘植物種類結(jié)構(gòu)分析.江西科學(xué),2010,28(4):549-554.]

    [19] Dai X,Wan RR,Yang GS et al.Responses of wetland vegetation in Poyang Lake,China to water-level fluctuations. Hydrobiologia,2016,773:35-47.

    [20] Liu Xuying,Guan Yanning,Guo Shan et al.Response on wetland vegetation distribution to hydrology regularity based on harmonic-time series analysis.JLake Sci,2016,28(1):195-206.DOI:10.18307/2016.0123.[劉旭穎,關(guān)燕寧,郭杉等.基于時(shí)間序列諧波分析的鄱陽湖濕地植被分布與水位變化響應(yīng).湖泊科學(xué),2016,28(1):195-206.]

    [21] Zhang Meng,Ni Leyi,Xu Jun etal.Annual dynamics of the wetland plants community in Poyang Lake in response towater-level variation.Research of Environmental Sciences,2013,26(10):1057-1063.[張萌,倪樂意,徐軍等.鄱陽湖草灘濕地植物群落響應(yīng)水位變化的周年動(dòng)態(tài)特征分析.環(huán)境科學(xué)研究,2013,26(10):1057-1063.]

    [22] Zhang Lili,Yin Junxian,Jiang Yunzhong et al.Relationship between hydrological conditions and vegetation communities in Poyang Lake national nature reserve of China.Advances in Water Science,2012,23(6):769-776.[張麗麗,殷峻暹,蔣云鐘等.鄱陽湖自然保護(hù)區(qū)濕地植被群落與水文情勢(shì)關(guān)系.水科學(xué)進(jìn)展,2012,23(6):769-776.]

    [23] Hu Zhenpeng,Ge Gang,Liu Chenglin etal.Structure of Lake PoyangWetland plantsecosystem and influence of lakewater level for the structure.Research ofEnvironmental Sciences,2010,19(6):597-605.[胡振鵬,葛剛,劉成林等.鄱陽湖濕地植物生態(tài)系統(tǒng)結(jié)構(gòu)及湖水位對(duì)其影響研究.長(zhǎng)江流域資源與環(huán)境,2010,19(6):597-605.]

    [24] Wu Qin,Yao Bo,Zhu Lili et al.Seasonal variation in plant biomass of Carex Cinerascens and its carbon fixation assessment in a typical Poyang Lakemarshland.Research ofEnvironmental Sciences,2012,21(2):215-219.[吳琴,堯波,朱麗麗等.鄱陽湖典型苔草濕地生物量季節(jié)變化及固碳功能評(píng)價(jià).長(zhǎng)江流域資源與環(huán)境,2012,21(2):215-219.]

    [25] Li Cong,Li Xiaowen,Zheng Yu et al.Evolution of wetland landscape pattern in Hengshui National Nature Reserve.Resources Science,2008,30(10):1571-1578.[黎聰,李曉文,鄭鈺等.衡水湖國(guó)家級(jí)自然保護(hù)區(qū)濕地景觀格局演變分析.資源科學(xué),2008,30(10):1571-1578.]

    [26] Qiao Weifeng,Sheng Yehua,F(xiàn)ang Bin etal.Land use change informationmining in highly urbanized area based on transfermatrix:A case study of Suzhou,Jiangsu Province.Geographical Research,2013,32(8):1497-1507.[喬偉峰,盛業(yè)華,方斌等.基于轉(zhuǎn)移矩陣的高度城市化區(qū)域土地利用演變信息挖掘——以江蘇省蘇州市為例.地理研究,2013,32(8):1497-1507.]

    [27] Gong Zhaoning,Zhang Yiran,Gong Huili et al.Evolution ofwetland landscape pattern and its driving factors in Beijing. Acta Geographica Sinica,2011,66(1):77-88.[宮兆寧,張翼然,宮輝力等.北京濕地景觀格局演變特征與驅(qū)動(dòng)機(jī)制分析.地理學(xué)報(bào),2011,66(1):77-88.]

    Landscape pattern changes of Carex and its response to water level in Lake Poyang Wetland

    ZHOU Yunkai1,2,BAIXiuling1,2??&NING Lixin1,2
    (1:Institute ofNatural Resources and Environment,Henan University,Kaifeng 475004,P.R.China)
    (2:College of Environment and Planning,Henan University,Kaifeng 475004,P.R.China)

    The formation and evolution of vegetation landscape are controlled by the hydrological processes in freshwater wetlands. Based on the Landsat TM/ETM+remote sensing images,the landscape information on Lake PoyangWetland in 1992,1999,2006 and 2012 was extracted by using decision tree classification.The spatial variation of Carex landscape was analyzed by means of landscape pattern index,transfermatrix and centriod shiftingmethods.Meanwhile,the relationship between Carex landscape and hydrological processeswas also preliminarily discussed.Results showed that the distribution and area of Carex landscape in autumn were influenced by the water level and recession process of Lake PoyangWetland during the study period,the slow recession of intermediate-elevation(14-16m)areas and fastemersion of low-elevation(12-14m)areas play amore important role in the formation and expansion of Carex landscape than thewater level.The spatial pattern Carex landscape in Lake PoyangWetland was closely related to the annualmean water level.In the year of low annualmean water level,the distribution elevation of Carex landscape is relatively low,some exposed mudflats were occupied by Carex landscape,and the original distribution areas of Carex landscape were partly replaced by Phragmites and Triarrhena landscape,the fragmentation degree of Carex landscape was relatively heavy. While in the yearwith high annualmean water level,many low bottomlandswere submerged,the elevation of Carex distributed was relatively high.Many regions of Phragmites and Triarrhena landscape were occupied by Carex,while the original distribution areas of Carex landscape became partlymudflatorwater landscape types.Affected by the risingwater level,Carex landscapemainly shifted and centered on the high zones of estuary of input rivers and lakeshore,thus the fragmentation degree of Carex landscape was relatively light.In addition,water level inter-annual fluctuation will affect the centroid position of Carex landscape,with centriodshifting to the lakeshore during the rising water level;vice verse to the lake center during the falling water level falling.

    Carex;landscape pattern;water level change;response;Lake Poyang Wetland

    DOI 10.18307/2017.0410

    ?2017 by Journal of Lake Sciences

    ?國(guó)家自然科學(xué)基金項(xiàng)目(41101089,41371450)資助.2016-03-18收稿;2016-10-14收修改稿.周云凱(1980~),男,博士,副教授;E-mail:ykzhou2009@126.com.

    ??通信作者;E-mail:bxling@163.com.

    猜你喜歡
    苔草蘆荻鄱陽湖
    ABSTRACTS
    鄱陽湖水系之潦河
    蘆荻
    夕照高郵湖
    《鄱陽湖生態(tài)系列插畫》
    水邊的蘆荻
    文苑(2018年22期)2018-11-19 02:54:30
    水邊的蘆荻
    鄱陽湖好風(fēng)光
    老友(2017年4期)2017-02-09 00:26:04
    灰化苔草葉片鉛脅迫的光譜響應(yīng)研究
    苔草特征特性及防治技術(shù)
    欧美一级a爱片免费观看看| 国产精品人妻久久久影院| 欧美一区二区精品小视频在线| 中文亚洲av片在线观看爽| 午夜福利18| 国产欧美日韩精品亚洲av| 两个人视频免费观看高清| 国产淫片久久久久久久久| 尤物成人国产欧美一区二区三区| 亚洲,欧美,日韩| 亚洲专区中文字幕在线| 丰满的人妻完整版| av视频在线观看入口| 麻豆精品久久久久久蜜桃| 日日夜夜操网爽| h日本视频在线播放| av天堂在线播放| 久久久久久久久久黄片| 亚洲欧美清纯卡通| 九色国产91popny在线| 亚洲av五月六月丁香网| 丝袜美腿在线中文| 美女大奶头视频| 亚洲精品粉嫩美女一区| 高清毛片免费观看视频网站| 真实男女啪啪啪动态图| 久久精品91蜜桃| 国产精品久久久久久精品电影| 综合色av麻豆| 亚洲成人久久爱视频| 精品久久久噜噜| 成人特级黄色片久久久久久久| 亚洲欧美日韩无卡精品| 女人十人毛片免费观看3o分钟| 免费在线观看成人毛片| 日日摸夜夜添夜夜添小说| 久久久久免费精品人妻一区二区| 日日摸夜夜添夜夜添小说| 丰满乱子伦码专区| 黄色女人牲交| 免费不卡的大黄色大毛片视频在线观看 | www.色视频.com| 欧美性感艳星| 麻豆精品久久久久久蜜桃| 色综合婷婷激情| 国产一区二区三区视频了| 日本 av在线| 亚洲av熟女| 国产视频内射| 看免费成人av毛片| avwww免费| 岛国在线免费视频观看| 久久久久久久久久久丰满 | 波多野结衣高清作品| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 亚洲精品一卡2卡三卡4卡5卡| 国产真实伦视频高清在线观看 | 又黄又爽又免费观看的视频| 一a级毛片在线观看| 极品教师在线免费播放| 色视频www国产| 亚洲国产日韩欧美精品在线观看| 午夜影院日韩av| 亚洲av日韩精品久久久久久密| 一本一本综合久久| 韩国av一区二区三区四区| 午夜福利高清视频| 国内毛片毛片毛片毛片毛片| 春色校园在线视频观看| 国产精品精品国产色婷婷| av在线观看视频网站免费| 国产三级中文精品| 国内少妇人妻偷人精品xxx网站| 国产精品自产拍在线观看55亚洲| 亚洲最大成人av| 欧美xxxx性猛交bbbb| 天堂av国产一区二区熟女人妻| 动漫黄色视频在线观看| 精品乱码久久久久久99久播| 在现免费观看毛片| 人人妻人人看人人澡| 亚洲av中文av极速乱 | 亚洲avbb在线观看| 欧美不卡视频在线免费观看| 香蕉av资源在线| 99riav亚洲国产免费| 亚洲avbb在线观看| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 97超级碰碰碰精品色视频在线观看| 国产激情偷乱视频一区二区| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 久久久国产成人免费| 天堂影院成人在线观看| 欧美日韩黄片免| 亚洲av免费在线观看| 一本一本综合久久| 国产精品人妻久久久久久| 国产精品久久久久久亚洲av鲁大| 成人特级av手机在线观看| 很黄的视频免费| 在线国产一区二区在线| 国产精品福利在线免费观看| 男女视频在线观看网站免费| 久久精品国产鲁丝片午夜精品 | 在线免费十八禁| 免费一级毛片在线播放高清视频| 午夜福利高清视频| 岛国在线免费视频观看| 午夜精品一区二区三区免费看| 麻豆久久精品国产亚洲av| 麻豆国产97在线/欧美| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 亚洲乱码一区二区免费版| 亚洲aⅴ乱码一区二区在线播放| 又爽又黄无遮挡网站| 日韩精品中文字幕看吧| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 成人高潮视频无遮挡免费网站| 少妇人妻一区二区三区视频| www.www免费av| 欧美一级a爱片免费观看看| 嫩草影视91久久| 欧美一区二区亚洲| 精品人妻一区二区三区麻豆 | 麻豆久久精品国产亚洲av| 91久久精品国产一区二区三区| 热99re8久久精品国产| a级一级毛片免费在线观看| 日日啪夜夜撸| 精品久久久久久久末码| 国产精品一区二区免费欧美| 日韩一本色道免费dvd| 欧美性感艳星| 国产色爽女视频免费观看| 91久久精品电影网| 91在线精品国自产拍蜜月| 亚洲图色成人| 亚洲 国产 在线| 综合色av麻豆| av黄色大香蕉| 神马国产精品三级电影在线观看| 亚洲av第一区精品v没综合| 极品教师在线视频| 日本黄大片高清| 乱人视频在线观看| 亚洲不卡免费看| 精品一区二区三区视频在线观看免费| 亚洲无线在线观看| 免费黄网站久久成人精品| 亚洲中文字幕日韩| 久久久久久大精品| 直男gayav资源| 身体一侧抽搐| 国产精品女同一区二区软件 | 久久久久久久久久黄片| 三级毛片av免费| 亚洲欧美精品综合久久99| 色综合婷婷激情| 人妻夜夜爽99麻豆av| 九九在线视频观看精品| 精品久久久久久久久久免费视频| 99热精品在线国产| 老司机深夜福利视频在线观看| 久久久久久久久久成人| 三级国产精品欧美在线观看| 亚洲精品成人久久久久久| 十八禁网站免费在线| 久久久久久大精品| 男人的好看免费观看在线视频| 麻豆精品久久久久久蜜桃| 啦啦啦韩国在线观看视频| 亚洲av五月六月丁香网| 91av网一区二区| 在线观看美女被高潮喷水网站| 国产单亲对白刺激| 欧美区成人在线视频| 久久久久久久午夜电影| 免费黄网站久久成人精品| 成人毛片a级毛片在线播放| 偷拍熟女少妇极品色| 日本在线视频免费播放| 国产三级中文精品| 毛片女人毛片| 国产伦一二天堂av在线观看| 国产精品国产三级国产av玫瑰| 一区二区三区免费毛片| 中国美白少妇内射xxxbb| 热99在线观看视频| 亚洲成a人片在线一区二区| 伦理电影大哥的女人| 波多野结衣高清作品| 国产单亲对白刺激| 国产午夜福利久久久久久| 亚洲性夜色夜夜综合| 亚洲欧美精品综合久久99| 国内精品美女久久久久久| 深爱激情五月婷婷| 免费观看精品视频网站| 最新中文字幕久久久久| 免费在线观看成人毛片| 两个人的视频大全免费| 色精品久久人妻99蜜桃| 亚洲国产日韩欧美精品在线观看| 亚洲av美国av| 欧美中文日本在线观看视频| 欧美一区二区精品小视频在线| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 国产色爽女视频免费观看| 亚洲中文字幕日韩| 精品国产三级普通话版| 免费在线观看成人毛片| 在线观看免费视频日本深夜| 极品教师在线免费播放| 九色国产91popny在线| 亚洲自拍偷在线| 高清毛片免费观看视频网站| 天天一区二区日本电影三级| 欧美一区二区国产精品久久精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区中文字幕在线| 毛片一级片免费看久久久久 | 久久久久久九九精品二区国产| 少妇裸体淫交视频免费看高清| 别揉我奶头 嗯啊视频| 永久网站在线| 亚洲国产精品成人综合色| 赤兔流量卡办理| 国产av在哪里看| 波多野结衣高清作品| 精品久久久久久久久亚洲 | 一区二区三区四区激情视频 | 在线a可以看的网站| 久久久久久伊人网av| 午夜福利在线观看吧| 麻豆精品久久久久久蜜桃| 欧美高清成人免费视频www| 亚洲欧美日韩无卡精品| 成人午夜高清在线视频| 亚洲avbb在线观看| 亚洲国产精品合色在线| 18禁黄网站禁片免费观看直播| 99九九线精品视频在线观看视频| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 日本与韩国留学比较| 国产伦人伦偷精品视频| 色av中文字幕| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 日本三级黄在线观看| 免费观看在线日韩| 国产一区二区亚洲精品在线观看| 最近在线观看免费完整版| 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区三区四区免费观看 | 一本精品99久久精品77| 亚洲av二区三区四区| 久久久久久伊人网av| 亚洲最大成人手机在线| 国产精品电影一区二区三区| 国产精品久久久久久亚洲av鲁大| 中文字幕熟女人妻在线| 最新在线观看一区二区三区| 日韩在线高清观看一区二区三区 | 午夜日韩欧美国产| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| а√天堂www在线а√下载| 亚洲精品456在线播放app | 一进一出抽搐gif免费好疼| 久久精品国产清高在天天线| 欧美一区二区国产精品久久精品| 少妇人妻精品综合一区二区 | 长腿黑丝高跟| 亚洲人与动物交配视频| 在线播放无遮挡| 国内揄拍国产精品人妻在线| 给我免费播放毛片高清在线观看| 国产色婷婷99| 久久国产乱子免费精品| 一区二区三区免费毛片| 99精品在免费线老司机午夜| 中文字幕av成人在线电影| 九九热线精品视视频播放| 国产一区二区三区在线臀色熟女| 日韩精品有码人妻一区| 毛片一级片免费看久久久久 | 久久久国产成人精品二区| 国产成人一区二区在线| 色噜噜av男人的天堂激情| 久久久久国内视频| 国内精品久久久久久久电影| 小蜜桃在线观看免费完整版高清| 国产av在哪里看| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 精品一区二区三区人妻视频| 亚洲内射少妇av| 99在线视频只有这里精品首页| aaaaa片日本免费| av专区在线播放| 久久热精品热| 嫩草影院精品99| 99热这里只有精品一区| 精品不卡国产一区二区三区| 亚洲成人精品中文字幕电影| 色哟哟·www| 乱码一卡2卡4卡精品| 女同久久另类99精品国产91| 亚州av有码| 国产美女午夜福利| 中文字幕久久专区| 成人国产一区最新在线观看| 色综合站精品国产| 亚洲av不卡在线观看| 成年女人看的毛片在线观看| 久久久久久国产a免费观看| 欧美成人免费av一区二区三区| 国产精品98久久久久久宅男小说| 国产精品免费一区二区三区在线| 69人妻影院| 尤物成人国产欧美一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲一级一片aⅴ在线观看| 亚洲成人久久爱视频| 成人国产综合亚洲| 日本色播在线视频| 在线观看美女被高潮喷水网站| xxxwww97欧美| 天天一区二区日本电影三级| 久久中文看片网| x7x7x7水蜜桃| 亚洲精品亚洲一区二区| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播放欧美日韩| 免费大片18禁| av黄色大香蕉| 日韩欧美精品v在线| 国产成人aa在线观看| 免费在线观看影片大全网站| 亚洲欧美清纯卡通| 不卡一级毛片| 亚洲av免费高清在线观看| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看| 亚洲午夜理论影院| a级毛片a级免费在线| 精品国内亚洲2022精品成人| 琪琪午夜伦伦电影理论片6080| 成人美女网站在线观看视频| 亚洲国产欧洲综合997久久,| 成人欧美大片| 国产精品日韩av在线免费观看| 久久国产精品人妻蜜桃| 国产精品,欧美在线| 亚洲第一电影网av| 日韩欧美三级三区| 日本一本二区三区精品| 哪里可以看免费的av片| 欧美潮喷喷水| 露出奶头的视频| 国产欧美日韩精品一区二区| 露出奶头的视频| 国产av麻豆久久久久久久| 波多野结衣高清作品| 国内精品久久久久精免费| 亚洲精品色激情综合| 成人美女网站在线观看视频| 色哟哟·www| 亚洲av中文字字幕乱码综合| 久久久精品欧美日韩精品| 免费av观看视频| 亚洲无线在线观看| 伦精品一区二区三区| 国产精品精品国产色婷婷| 91狼人影院| 极品教师在线免费播放| 成人永久免费在线观看视频| 俄罗斯特黄特色一大片| 国产精品亚洲美女久久久| 午夜亚洲福利在线播放| 国产精品乱码一区二三区的特点| 噜噜噜噜噜久久久久久91| netflix在线观看网站| 99九九线精品视频在线观看视频| 日韩一区二区视频免费看| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 国产精品人妻久久久影院| 日韩大尺度精品在线看网址| 国产精品三级大全| 国产午夜精品论理片| 日本爱情动作片www.在线观看 | 老熟妇仑乱视频hdxx| 精品欧美国产一区二区三| 99久久久亚洲精品蜜臀av| 亚洲人与动物交配视频| 国产一区二区三区av在线 | 亚洲最大成人手机在线| 国产女主播在线喷水免费视频网站 | 99riav亚洲国产免费| 免费av不卡在线播放| 俺也久久电影网| 精品一区二区三区人妻视频| 中文字幕精品亚洲无线码一区| 国产高清视频在线播放一区| 午夜爱爱视频在线播放| 午夜免费成人在线视频| 国产亚洲欧美98| 日本与韩国留学比较| 如何舔出高潮| 国产91精品成人一区二区三区| 中文在线观看免费www的网站| 香蕉av资源在线| 国产成人a区在线观看| 欧美精品国产亚洲| 小蜜桃在线观看免费完整版高清| 国产精品精品国产色婷婷| 伊人久久精品亚洲午夜| 制服丝袜大香蕉在线| 精品国内亚洲2022精品成人| 99热这里只有精品一区| 国产亚洲精品久久久久久毛片| 亚洲一级一片aⅴ在线观看| 麻豆精品久久久久久蜜桃| 亚洲成人免费电影在线观看| 亚洲最大成人av| 联通29元200g的流量卡| 日韩欧美一区二区三区在线观看| 国产成人av教育| 波野结衣二区三区在线| 国产伦人伦偷精品视频| 国产日本99.免费观看| 欧美最新免费一区二区三区| a在线观看视频网站| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区精品| 成人av一区二区三区在线看| 一级黄片播放器| 亚洲美女黄片视频| 欧美丝袜亚洲另类 | 能在线免费观看的黄片| 欧美日本亚洲视频在线播放| 日本三级黄在线观看| 深夜精品福利| 午夜免费男女啪啪视频观看 | 51国产日韩欧美| 国产一区二区三区av在线 | 非洲黑人性xxxx精品又粗又长| 日本五十路高清| 男人狂女人下面高潮的视频| 久久婷婷人人爽人人干人人爱| 久久久久久久久久久丰满 | 91麻豆av在线| 精品国产三级普通话版| 国产精品一区二区免费欧美| 丰满的人妻完整版| 国产色婷婷99| 欧美日韩乱码在线| 婷婷精品国产亚洲av| 国产精品嫩草影院av在线观看 | 亚洲五月天丁香| 一区二区三区激情视频| 三级毛片av免费| 精品人妻1区二区| 一级av片app| 亚洲乱码一区二区免费版| 在线免费观看不下载黄p国产 | 丝袜美腿在线中文| 国产精品野战在线观看| 久久久精品欧美日韩精品| 国产精品国产高清国产av| 日本一二三区视频观看| 国产三级中文精品| 九九爱精品视频在线观看| 国产欧美日韩精品亚洲av| 国产在线男女| 99久久无色码亚洲精品果冻| 日韩欧美 国产精品| 国产单亲对白刺激| 国产免费av片在线观看野外av| 国产aⅴ精品一区二区三区波| 国产av不卡久久| 国产免费男女视频| 一边摸一边抽搐一进一小说| 国产爱豆传媒在线观看| 亚洲欧美精品综合久久99| 国产精品,欧美在线| 18+在线观看网站| 亚洲欧美日韩卡通动漫| 人妻夜夜爽99麻豆av| 国产视频内射| 我的女老师完整版在线观看| netflix在线观看网站| 亚洲五月天丁香| 日韩欧美精品v在线| 国内久久婷婷六月综合欲色啪| 一区福利在线观看| 久久精品国产99精品国产亚洲性色| 神马国产精品三级电影在线观看| 久久久国产成人免费| 精品不卡国产一区二区三区| 精品久久久久久久末码| 又紧又爽又黄一区二区| 男女视频在线观看网站免费| 亚洲经典国产精华液单| 国产成人aa在线观看| 欧美bdsm另类| 啦啦啦观看免费观看视频高清| 欧美日韩中文字幕国产精品一区二区三区| 亚洲性久久影院| 欧美xxxx黑人xx丫x性爽| 男女啪啪激烈高潮av片| 亚洲精品乱码久久久v下载方式| 国产高清视频在线观看网站| 国产免费男女视频| 一本久久中文字幕| 村上凉子中文字幕在线| 天堂影院成人在线观看| 国产精品精品国产色婷婷| 欧美区成人在线视频| 免费在线观看成人毛片| 我的女老师完整版在线观看| 亚洲性久久影院| 午夜福利视频1000在线观看| 最好的美女福利视频网| 色尼玛亚洲综合影院| 国产黄a三级三级三级人| 国产精品亚洲美女久久久| 日本一二三区视频观看| 国产一区二区三区av在线 | videossex国产| 日韩一本色道免费dvd| 色av中文字幕| 国产精品日韩av在线免费观看| 少妇熟女aⅴ在线视频| 99在线视频只有这里精品首页| 成人av一区二区三区在线看| 春色校园在线视频观看| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 欧美最新免费一区二区三区| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 亚洲人成网站在线播| 十八禁网站免费在线| 欧美日本视频| 国产成人影院久久av| 国产精品久久久久久亚洲av鲁大| 国产久久久一区二区三区| 国产一区二区激情短视频| 欧美日韩瑟瑟在线播放| 日韩亚洲欧美综合| 美女免费视频网站| 中文字幕久久专区| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区| 午夜激情福利司机影院| 国产91精品成人一区二区三区| 18+在线观看网站| 高清毛片免费观看视频网站| 精品人妻视频免费看| 99热精品在线国产| 性插视频无遮挡在线免费观看| 国内少妇人妻偷人精品xxx网站| 国产亚洲欧美98| 国产精品av视频在线免费观看| 亚洲在线观看片| 久久久久久久精品吃奶| 老熟妇仑乱视频hdxx| 一卡2卡三卡四卡精品乱码亚洲| 亚洲自偷自拍三级| 欧美激情久久久久久爽电影| 乱系列少妇在线播放| 丰满的人妻完整版| 亚洲成av人片在线播放无| 禁无遮挡网站| 欧美bdsm另类| 国产精品人妻久久久久久| 禁无遮挡网站| 男人舔女人下体高潮全视频| 国产欧美日韩精品一区二区| 成人无遮挡网站| 深夜a级毛片| av专区在线播放| 性色avwww在线观看| 美女大奶头视频| 真人做人爱边吃奶动态| 一个人观看的视频www高清免费观看| 国产爱豆传媒在线观看| 中文字幕人妻熟人妻熟丝袜美| 黄色视频,在线免费观看| 在线天堂最新版资源| 女的被弄到高潮叫床怎么办 | 欧美精品啪啪一区二区三区| 午夜精品在线福利| 日韩中字成人| 麻豆一二三区av精品| 美女cb高潮喷水在线观看| 色噜噜av男人的天堂激情| 亚洲av日韩精品久久久久久密| 国产乱人伦免费视频| 亚洲国产色片| 久久99热这里只有精品18| 搞女人的毛片| 精品人妻熟女av久视频| 美女黄网站色视频| 亚洲最大成人中文| 中文字幕精品亚洲无线码一区|