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

    近60年鄱陽湖東部湖灣水文連通變化及其對濕地植物與候鳥的影響*

    2023-01-13 07:21:58熊麗黎劉建新李寬意吳召仕郭玉銀王仕剛鄧燕青歐陽千林
    湖泊科學 2023年1期
    關(guān)鍵詞:湖灣苦草連通性

    熊麗黎,劉建新,李寬意,閔 翔,吳召仕,郭玉銀,王仕剛,鄧燕青,歐陽千林

    (1:中國科學院南京地理與湖泊研究所,南京 210008) (2:江西省水文監(jiān)測中心,南昌 330002) (3:江西省鄱陽湖水文生態(tài)監(jiān)測研究重點實驗室,南昌 330002) (4:中國科學院大學,北京 100864) (5:鄱陽湖水文水資源監(jiān)測中心,南昌 330002)

    水文連通指以水為介質(zhì)的物質(zhì)、能量及生物在水文循環(huán)要素內(nèi)部或要素之間的傳輸轉(zhuǎn)移[1-2],在湖泊洪水調(diào)蓄、水質(zhì)凈化和生物保護等方面發(fā)揮重要作用。水文連通性變化與水文情勢相互作用,進而影響生境面積、湖泊植被靜態(tài)格局和水生生物動態(tài)遷移、覓食和繁殖特征等[3-4]。維持連通性是河湖水系管理和健康維護的重要要求[5]。然而近半個世紀以來,在氣候變化和人類適應性水文調(diào)控的影響下,濕地局部水文連通性嚴重降低,造成濕地面積急劇萎縮和功能嚴重退化,損害水域生態(tài)完整性,制約流域經(jīng)濟、社會、環(huán)境的協(xié)同發(fā)展[4,6-7]。Grill等構(gòu)建了河流連通狀態(tài)指數(shù)(connectivity status index),結(jié)果表明大壩阻隔及水庫調(diào)節(jié)破壞了全球大約48.2%的河段水文連通性,影響了河流及河岸生態(tài)系統(tǒng)的自我修復能力[8]。因此,水文連通研究對于濕地生態(tài)環(huán)境評估與健康維護等顯得尤為重要。

    鄱陽湖是長江流域典型的吞吐型、季節(jié)性湖泊,水文特征季節(jié)性差異明顯。漲水時,主湖區(qū)水流將營養(yǎng)物質(zhì)輸送至洲灘、碟形湖及湖濱帶,水生動物也進入洪泛區(qū)覓食、繁殖甚至避難,退水時歸槽水流帶走腐殖質(zhì),魚類等水生動物回歸主湖區(qū),有效維持了洲灘、碟形湖及湖濱帶生物群落結(jié)構(gòu)的穩(wěn)定性[9]。碟形湖物種多樣性豐富,是鄱陽湖80%以上候鳥(如白枕鶴(Ciconiaalba)、東方白鸛(Ciconiaboyciana)等)的棲息地和覓食場所[10-11]。濕地植物群落沿高程的分布特征是在多種環(huán)境因子綜合作用下形成的,水文過程通常被認為是主導因子[12]。主湖區(qū)與碟形湖的水文連通變化迅速,對水位、淹沒面積等的影響十分明顯[13-14]。水文連通差異決定了碟形湖植被分布和生長特征的不同[15]。2003年以后,長江對鄱陽湖的倒灌量平均下降了近90%,主湖區(qū)枯水期提前且延長、枯水位進一步降低,導致洲灘濕地的出露時間延長[16]。根據(jù)1983年和2013年兩次鄱陽湖科考成果,雖然因分布上下限拓寬各1 m,薹草群落面積增幅約為68.9%,但鄱陽湖沉水植物面積整體縮減37.7%,群落物種由5~8種降至3~5種,苦草(Vallisnerianatans)替代竹葉眼子菜(Potamogetonwrightii)成為優(yōu)勢物種[17]。

    為了有效反映季節(jié)性洪泛湖泊水文連通狀況,學者從不同角度提出了多種水文連通度的定量分析方法,如與水文特征相關(guān)的輸移時間、換水周期等,與水動力特征相關(guān)的水沙通量、水流速率等,以及與生物相關(guān)的生物遷移能力等[18-20]。針對鄱陽湖季節(jié)性洪泛過程形成的河道-洲灘-碟形湖之間的水文連通過程,已有研究采用統(tǒng)計方法,引入水文連通函數(shù)(connectivity function)來表征湖區(qū)空間上縱向和橫向水文連通程度的變化,并有效反映了湖區(qū)水動力學條件、碟形湖水面擴張和萎縮的程度[21-22]。結(jié)合水面積、水深、流速等水文參數(shù),Li等基于湖泊洪泛水文水動力學,考慮連通性與生態(tài)環(huán)境指標之間的耦聯(lián)關(guān)系,重新定義了水文連通性的概念及其意義[23]?;谶B通的定義,圍繞水鳥、魚類、浮游藻類和大型底棲動物等,Tan等進一步探討了鄱陽湖洪泛區(qū)水文連通與生態(tài)環(huán)境要素之間的聯(lián)系,提升了對洪泛區(qū)水文連通及其生態(tài)環(huán)境服務(wù)功能的認識[24]。

    目前,相關(guān)研究主要是基于系統(tǒng)和宏觀視角,側(cè)重于水文連通的變化過程和影響機理,研究方式多基于遙感數(shù)據(jù)或者模型模擬[25-26]。然而,針對大湖洪泛的典型區(qū)域,水文連通性的長期變化特征以及有關(guān)水文連通的影響研究尚需加強,基于長時間序列實測數(shù)據(jù)以及地理地形測量資料,聚焦不同時期的水文連通改變程度及對鄱陽湖典型洪泛系統(tǒng)水生生境的影響工作仍較為少見。同時1962-2002年鄱陽湖主湖區(qū)水位呈周期性變化且略有上升趨勢,但受江湖關(guān)系變化的影響,2003年以后水位呈現(xiàn)顯著下降變化且波動較為劇烈[15]。與南磯山等其他洪泛區(qū)不同的是,東部湖灣枯水期與主湖區(qū)存在水文連通出口僅有一處,為相對獨立的水生態(tài)系統(tǒng),水量交換特征與碟形湖類似,便于分析水文連通性改變程度以及潛在影響,且因地形和距離的差異而存在生態(tài)差異。因此,本文分析對比2003年前后鄱陽湖東部湖灣和主湖區(qū)的水文連通性變化及其影響,為進一步探討碟形湖和洪泛區(qū)的水生態(tài)系統(tǒng)演變提供科學依據(jù)。主要研究目標為:(1)基于近60年實測水位數(shù)據(jù)以及近年湖泊地理地形測量資數(shù)據(jù)分析2003年前后主湖區(qū)水文變異程度及其生態(tài)效應影響;(2)從連通天數(shù)、發(fā)生時間等方面定量反映與主湖區(qū)的水文連通變化程度;(3)探討水文連通改變對濕地植物分布及候鳥棲息的潛在影響,提出相關(guān)保護對策。對東部湖灣這一典型區(qū)域開展水文連通性變化原因及其對濕地植物、候鳥的影響研究,有利于揭示鄱陽湖洪泛區(qū)水文連通性及其對濕地植被、候鳥影響的長期變化趨勢,為其水生態(tài)保護提供依據(jù)。

    1 研究區(qū)與方法

    1.1 研究區(qū)及數(shù)據(jù)

    鄱陽湖東部湖灣,面積約為580 km2,占鄱陽湖面積的17.7%,系湖區(qū)最大湖汊(圖1)。區(qū)域分布有企湖、珠池湖、白池湖、漢池湖等多個碟形湖,這些碟形湖的水量交換也受東部湖灣水文連通性影響。東部湖灣月平均水深變化范圍為4~10 m,每年3-9月與主湖區(qū)保持良好的水文連通且存在著周期性的水量交換,9月中旬以后東部湖灣流向主湖區(qū)徑流量驟減,10月中旬至次年2月與主湖區(qū)沒有水量交換,僅在虎頭渡口有徑流流出。東部湖灣生長有薹草、蘆葦?shù)葷裆参镆约爸袢~眼子菜、苦草、黑藻(Hydrillaverticillata)等沉水植物,是鄱陽湖國家鯽魚鯉魚產(chǎn)卵場保護區(qū)和鄱陽湖都昌候鳥省級自然保護區(qū)。近年來,多次現(xiàn)場調(diào)查發(fā)現(xiàn)東部湖灣植被群落發(fā)生較大改變[21]。

    圖1 鄱陽湖東部湖灣位置及遙感影像圖Fig.1 Location of theeastern bay in Lake Poyang and its remote sensing images

    本文選取棠蔭水位站和塘美水位站進行計算分析。塘美站位于東部湖灣,建于2018年7月;棠蔭站位于鄱陽湖主湖區(qū)鄰近東部湖灣,建于1961年11月。數(shù)據(jù)序列選取棠蔭站1962年1月1日-2019年10月22日和塘美站2018年8月8日-2020年10月22日的逐日水位數(shù)據(jù)。利用塘美和棠蔭兩站同時期水位資料分析兩站水位關(guān)系,基于棠蔭站水位數(shù)據(jù)分析主湖區(qū)水文變異程度。鄱陽湖湖底地形采用2010年1∶1萬實測高程數(shù)據(jù)。數(shù)據(jù)來源于江西省水文監(jiān)測中心。

    1.2 主要方法

    水文變異會導致生物棲息地發(fā)生改變,Richter等[27]提出的IHA (indicators of hydrologic alteration)/RVA (range of variability approach) 法可定量評價生態(tài)環(huán)境水文過程,在河流湖泊水文變異及其生態(tài)效益評估中得到了廣泛運用[27-30]。本文選擇IHA/RVA法分析評價主湖區(qū)水位變異程度及潛在生態(tài)影響,進一步分析水位變異對東部湖灣和主湖區(qū)的水文連通關(guān)系的影響。從適宜生境面積、連通天數(shù)以及非連通期發(fā)生時間3個方面量化表征水文連通性的改變程度。其中,適宜生境面積表征濕地植被適宜生境的大小,連通天數(shù)可反映水量與主湖區(qū)的交換頻次,非連通發(fā)生時間則在一定程度上體現(xiàn)越冬候鳥取食時間及范圍。

    1.2.1 IHA/RVA法 因鄱陽湖區(qū)不存在斷流情況,本文采用的IHA方法包括月平均水位、年極端水位、極端水位出現(xiàn)時間、高低水位頻率和持續(xù)時間、水位變化率5組32個水文指標[30],可反映水文變化對濕地植物群落分布、生物覓食和棲息地塑造的影響程度(表1)。

    表1 IHA評價指標及對濕地生態(tài)系統(tǒng)的影響Tab.1 Index of IHA method and their ecological impact on wetland

    通常,RVA閾值上下限為參照時段IHA指標發(fā)生率的75%和25%,該范圍可認為能滿足水生生物生態(tài)需水的變動要求。湖泊水位變化程度可采用分項指標改變度(Ai)和整體改變度(A0)來表示,計算公式見式(1)和(2)。分項指標改變度指每個指標的改變程度,而整體改變度綜合反映水位改變程度及其影響。當指標改變度小于33%時,屬于影響程度較小的低度改變;當改變度大于33%并小于67%時,屬于具有一定影響的中度改變;當改變度大于67%時,屬于影響程度較大的高度改變。

    (1)

    (2)

    式中,Ai為第i個指標的水文改變度;ei和e0分別為水位變化時段指標在RVA閾值的上下限之間的實際年數(shù)和預期年數(shù)。預期年數(shù)為閾值上下限的中間值(50%)與變化時段總年數(shù)(17年)之積,本文為8.5年。

    1.2.2 適宜生境面積與連通天數(shù) 利用ArcGIS 10.2從鄱陽湖東部湖灣DEM數(shù)據(jù)中提取不同高程的水面面積,建立東部湖灣水位-面積關(guān)系。東部湖灣濕生植物和沉水植物的優(yōu)勢種分別為薹草和苦草[31],是鴻雁(Ansercygnoides)、東方白鸛等越冬候鳥的主要食物來源[32]。薹草適宜生境主要取決于適宜淹沒時長和適宜水深,分別為71~225 d和0.16~1.16 m[33-34];而苦草和大多數(shù)候鳥覓食的適宜水深分別為0~1.0 m和0.2~0.4 m。根據(jù)東部湖灣水位-面積關(guān)系曲線,利用不同水位差計算不同水位情況下的水深情況,分析薹草、苦草生長以及候鳥覓食的適宜水深范圍對應的生境面積,建立東部湖灣水位與薹草、苦草生長以及候鳥覓食生境面積的關(guān)系,以反映東部湖灣適宜生境面積隨水位的變化情況。此外,對于薹草適宜生境面積,在適宜淹沒水深對應的薹草生境面積基礎(chǔ)上,進一步考慮其對適宜淹沒時長的影響。

    通常,每年3-6月鄱陽湖碟形湖與主湖區(qū)逐漸連通,7-9月保持高水位完全連通,10月后水位下降形成獨立的子湖體。對比棠蔭站和塘美站逐日水位變化,當水位高于閾值12.5 m時,東部湖灣與主湖區(qū)處于“連通”狀態(tài),兩者水位保持同步變化;反之,當水位低于該閾值時,主湖區(qū)和東部湖灣處于“非連通”狀態(tài),兩者水位變化存在著明顯差異。當東部湖灣與主湖區(qū)處于“連通”狀態(tài)時,棠蔭站和塘美站水位變化相關(guān)性顯著,可通過建立兩者水位相關(guān)關(guān)系方程,利用棠蔭站長序列的水位數(shù)據(jù)延長塘美站水位數(shù)據(jù),為東部湖灣薹草、苦草生長以及候鳥覓食適宜生境面積的分析提供數(shù)據(jù)支持。根據(jù)東部湖灣與主湖區(qū)連通時的水位閾值,可統(tǒng)計獲得研究時段內(nèi)東部湖灣處于連通狀態(tài)時各月份的平均天數(shù),即月平均連通天數(shù)。

    (3)

    式中,L為一年的天數(shù);Ωz轉(zhuǎn)化為弧度,表示非連通期發(fā)生時間,弧度0、2π分別表示1月1日和12月31日。對于某時段,年內(nèi)第l日(儒略日)處于非連通狀態(tài)共計nl天,該時段非連通狀態(tài)的平均發(fā)生時間和平均合成長度的計算公式為:

    (4)

    (5)

    其中,

    (6)

    2 結(jié)果與分析

    2.1 水位與水文連通關(guān)系

    2.1.1 東部湖灣水位-面積關(guān)系 以塘美站為依據(jù)分析鄱陽湖東部湖灣水位與水面面積及適宜生境面積的關(guān)系(圖2)。水面面積隨著塘美站水位的升高而增加,但其增加值先增加后減少且大于15.0 m以后變化平緩,水面面積最大增加值為100.3 km2,相應的塘美站水位為13.5 m。不同水位薹草、苦草生長和候鳥覓食的生境面積呈先增加后減少的變化。當塘美站水位為13.5、14.4和13.6 m時,東部湖灣適宜薹草、苦草生長和候鳥覓食的生境面積達到最大,分別為187.9、194.4和40.1 km2,占相應水位水面面積的58.7%、49.3%和13.2%。在塘美站水位低于12.8、14.0和12.4 m時,薹草、苦草生長和候鳥覓食的適宜生境面積受東部湖灣水位變化影響較大;但水位高于15.6、15.8和13.2 m時,適宜水深對應的水面面積變化平緩,此時水位的變化對東部湖灣生境的影響較小。

    圖2 東部湖灣水位-水面面積曲線(a)及水位-適宜生境面積曲線(b)Fig.2 Curves of water level and area (a), water level and suitable habitat area (b) of the eastern bay

    2.1.2 水文連通關(guān)系 根據(jù)2018年8月8日-2020年10月22日棠蔭站和塘美站逐日水位數(shù)據(jù)變化過程(圖3a)和東部湖灣水位面積關(guān)系,分析2018-2020年東部湖灣適宜生境面積變化過程(圖3b)。從棠蔭站和塘美站逐日水位變化結(jié)合實際調(diào)查可判斷,棠蔭站水位12.5 m為東部湖灣與主湖區(qū)水文連通的閾值。當棠蔭站水位高于12.5 m時,棠蔭站和塘美站水位保持同步變化,東部湖灣與鄱陽湖主湖區(qū)保持全面連通,屬“連通”狀態(tài);當棠蔭站水位低于12.5 m時,棠蔭站和塘美站水位差異明顯,東部湖灣脫離主湖區(qū)獨立成子湖,屬“非連通”狀態(tài)。通常,每年3-9月東部湖灣與主湖區(qū)處于“連通”狀態(tài),年內(nèi)其它月份多為“非連通”狀態(tài)。當塘美站水位在13~16 m范圍變化時,薹草和苦草生長適宜生境面積較大;而當塘美站水位小于13 m或大于16 m時,薹草和苦草生長的適宜生境面積急劇下降,而薹草主要受淹沒時長不足或過長的影響,在東部湖灣的適宜生境面積消失。

    圖3 2018-2020年東部湖灣水位(a)與適宜生境面積(b)Fig.3 Daily water level (a) and suitable habitat area (b) of the eastern bay from 2018 to 2020

    根據(jù)棠蔭站和塘美站逐日水位變化,分析“連通”、“非連通”兩種狀態(tài)東部湖灣與主湖區(qū)的水位相關(guān)性,如圖4所示。在水文連通狀態(tài),棠蔭站和塘美站年內(nèi)水位變化范圍較大,約為9.0 m,但兩者水位差始終保持在0.30 m左右,決定系數(shù)(R2)達0.997,可用于構(gòu)建兩站水位的函數(shù)關(guān)系。非連通狀態(tài)通常處于枯水期,棠蔭站水位主要受五河來水的影響,而塘美站水位則受到上游漳田河和潼津河來水影響,兩者水位R2為0.667。非連通狀態(tài)棠蔭站和塘美站的平均水位分別為13.0和11.5 m,水位變化范圍分別為12.5~13.5 m和10.7~12.5 m,為分析枯水期適宜生境面積變化提供依據(jù)。

    圖4 2018-2020年棠蔭站和塘美站水位差異(a)以及連通狀態(tài)(b)和非連通狀態(tài)(c)的相關(guān)關(guān)系Fig.4 Differencess of daily water level (a), their correlation relationship under hydrological connectivity (b) and non-connectivity (c) between Tangyin Station and Tangmei Station from 2018 to 2020

    2.2 主湖區(qū)水位變異程度

    鄱陽湖主湖區(qū)年際年內(nèi)水位變化較大,1962-2019年棠蔭站年最高水位和最低水位平均值分別為19.2和11.3 m,2003年以后平均降低約0.7 m,水位差范圍由4.4~10.9 m變?yōu)?.6~10.1 m(圖5)。以1962-2002年為基準期,分析2003-2019年主湖區(qū)IHA指標改變程度(表2)。第1組月平均水位改變度為27.8%,2003年以后月平均水位都有不同程度的下降,12項指標中多數(shù)變化以低改變程度為主,其中 10-12 月為中改變程度;第2組年極端水位改變度為33.4%,11項極端水位指標為中低改變程度,其中年均1、3、7以及90 d最小值4項指標均為中改變程度,2004、2007和2008年年均7 d最低水位均下降至10 m以下,不同時長的最大值指標為低改變程度,但其下降幅度多數(shù)為1 m以上,處于16.5~18.5 m之間;第3組年極端水位出現(xiàn)時間改變度為65.8%,有90%以上的年份超過RVA上限值67%,平均提前約15 d,屬于高改變度,而最高水位出現(xiàn)時間平均推后約5 d,屬于低改變程度;第4組高低水位頻率和持續(xù)時間改變度為30.6%,年高、低水位頻率為中改變程度,相應地持續(xù)時間則為低改變程度,表現(xiàn)為低水位持續(xù)時間延長,平均增加6.5 d,而高水位持續(xù)時間縮短,平均減少6 d;第5組水位變化率改變度為77.1%,3項指標變異程度較高,其中年平均水位增長率和年水位逆轉(zhuǎn)次數(shù)為高改變程度,而平均水位減少率為中改變程度。從32個IHA指標可知,低、中、高改變程度分別為19、10和3個,2003年以后鄱陽湖主湖區(qū)水位降低集中在枯水期的10-12月,同時存在著湖區(qū)低水位進一步降低且年水位波動頻率增加的變化。水位整體改變度為40.2%,屬于中改變程度,可導致洪泛區(qū)與主湖區(qū)的水文連通關(guān)系減弱。

    圖5 1962-2019年棠蔭站年水位序列(a)及IHA指標(b)Fig.5 Water level series (a) and its IHA (b) of Tangyin Station from 1962 to 2019

    表2 2003-2019年棠蔭站IHA指標改變情況Tab.2 Summary on IHA of Tangyin Station from 2003 to 2019

    2.3 連通狀況及適宜生境面積變化

    2.3.1 連通天數(shù)與適宜生境面積 分析1962-2002年和2003-2019年不同時段的平均水文連通天數(shù)變化(表3),2003年前后東部湖灣與主湖區(qū)年內(nèi)不同時段水文連通時長均有不同程度的減少。1962-2002年東部湖灣與主湖區(qū)年均水文連通天數(shù)為291.0 d,而2003-2019年水文連通天數(shù)為244.8 d,減少了約46.2 d,變化幅度為15.9%。對于不同季度,10-12月水文連通天數(shù)的變化最大,縮短了27.5 d,變化幅度達46.8%,其次1-3月水文連通天數(shù)縮短了約8.3 d,變化幅度為16.4%,這主要由于10-12月和1-3月受長江與鄱陽湖關(guān)系的改變,鄱陽湖主湖區(qū)枯水期提前且低水時段拉長;4-6月和7-9月為豐水期,受流域“五河”來水以及長江水位頂托影響,兩個時段鄱陽湖水位變化較小,連通天數(shù)縮短了5 d左右,變化幅度約為6%。2003年前后東部湖灣與主湖區(qū)年內(nèi)不同時段水文連通時長均有不同程度的減少。

    表3 1962-2019年東部湖灣與主湖區(qū)連通天數(shù)Tab.3 Number of days under hydrological connectivity between the eastern bay and the main lake of Lake Poyang from 1962 to 2019

    對于薹草和苦草生長,分析1962-2002年和2003-2019年兩個時期3-10月逐月適宜生境面積情況變化(表4)。對比2003年前后,3-10月東部湖灣平均水位下降了0.6 m,除了6月份水位不變以外,其他7個月東部湖灣平均水位均呈減少的變化,特別是9和10月,分別減少1.2和1.0 m。對于薹草適宜生境而言,2003-2019年較1962-2002年適宜水深相應的水面面積增加值最大的月份為9月,增加幅度為73.9%;而減少的月份為7月,降幅為4.1%。對于苦草適宜生境而言,2003-2019年較1962-2002年適宜水深相應的水面面積增加值最大的月份依次為9和10月,增幅分別為39.2%和21.4%;而減少的月份為3月和6月,降幅分別為10.4%和10.5%。2003年前后東部湖灣平均水位有所下降,薹草和苦草等典型濕地植物的分布將向湖中心轉(zhuǎn)移,但其適宜水深對應的生境面積有所增加。

    表4 1962-2019年東部湖灣薹草和苦草月適宜生境面積Tab.4 Monthly suitable habitat area for Carex cinerascens and Vallisneria natans in the eastern bay from 1962 to 2019

    2.3.2 非連通天數(shù)與發(fā)生時間 分析1962-2002年和2003-2019年兩個時期非連通天數(shù)和發(fā)生時間變化(表5,圖6),結(jié)果表明東部湖灣與主湖區(qū)的非連通狀態(tài)發(fā)生了明顯的改變。1962-2002年東部湖灣與主湖區(qū)平均年非連通天數(shù)為74.2 d,而2003-2019年為119.5 d,平均年非連通天數(shù)明顯增加,增加幅度達61.2%,其中2011年非連通天數(shù)達257 d,占全年的70%以上;1962-2002年東部湖灣非連通期平均角度為6.2°(1月6日),而 2003-2019年為-6.9° (12月24日),平均提前了約13 d,且2003年以后非連通期開始時間明顯提前了1個月,大約提前至10月中旬;1962-2002年東部湖灣與主湖區(qū)非連通期從11月上旬至2月上旬,而2003-2019年則從10月中旬至2月中旬,兩個時期非連通期合成長度分別為0.81和0.64,東部湖灣非連通狀態(tài)在年內(nèi)的時間跨度更長。

    表5 1962-2019年東部湖灣非連通狀態(tài)的變化Tab.5 Results on the non-connectivity between the eastern bay and main lake of Lake Poyang from 1962 to 2019

    圖6 1962-2019年東部湖灣非連通期天數(shù)(a)及方向(b)Fig.6 Days (a) and directions (b) under the non-connectivity between the eastern bay and main lake of Lake Poyang from 1962 to 2019

    3 討論

    3.1 水文連通性變化的主要原因

    本研究認為鄱陽湖東部湖灣近年來與主湖區(qū)水文連通性呈下降趨勢,尤其2003年以來連通時間縮短,這主要受主湖區(qū)水位變化的影響。2003年以后主湖區(qū)枯水期提前且時間延長,枯水期水位下降,與前人關(guān)于鄱陽湖水位變化的研究[17]基本一致,表現(xiàn)為10-12月平均水位、年極端水位、年最低水位出現(xiàn)時間、高低水位頻率以及水位變化率等水位指標發(fā)生中高度程度的改變。主湖區(qū)水位主要受到江湖關(guān)系、湖泊地形、氣候變化等影響[36-38]。根據(jù)1973-2011年多時相衛(wèi)星遙感數(shù)據(jù)分析結(jié)果,三峽工程蓄水運行引發(fā)的長江與鄱陽湖關(guān)系變化是近20年鄱陽湖汛末水位突變的重要因素之一[39],此外,2000年以后(特別是2003-2016年間)大規(guī)模的采砂造成湖盆局部形態(tài)發(fā)生改變[40],導致北部湖區(qū)湖底大幅度下切變深、局部擴寬,10和11月棠蔭站平均水位分別下降2.81和2.05 m[37,40-42]。氣候變化導致2003-2012年多年平均降雨量較1956-2002年偏少3.62%[38]。

    三峽工程的運行以及采砂、氣候變化不同程度改變了主湖區(qū)水位變化,從而引起洪泛區(qū)與主湖區(qū)的連通性,而東部湖灣作為受主湖區(qū)水位影響最為密切的區(qū)域,其水文連通性變化更為敏感且影響深遠。

    3.2 水文連通性變化對濕地植物和候鳥生境的影響

    隨著水文連通性的下降,東部湖灣水面面積呈現(xiàn)萎縮趨勢,植被及候鳥群落結(jié)構(gòu)及分布特征發(fā)生較大變化。2003年以后,水文連通性變化導致草洲提前出露,苦草、黑藻、竹葉眼子菜等沉水植物種類和數(shù)量有所減少,但薹草等濕生植物移向湖心且分布有增加的趨勢,改變了濕地植物的優(yōu)勢物種和群落結(jié)構(gòu)[43],鄱陽湖沉水植物分布上限降低1 m[17]。近年來,東部湖灣濕地植物的適宜生境面積整體上升,但對于沉水植物和濕生植物的影響各不相同,尤其在9-10月份恰好是薹草秋草生長期以及苦草的果期[42]。鄱陽湖汛末水位突變導致9-10月份鄱陽湖枯水提前、枯水水位下降[38-39],將減少淹沒時間,適宜生長面積增加,有利于薹草等濕生植物生長,使其面積呈上升趨勢。對于苦草來說,3和6月分別是其萌發(fā)期和開花期,適宜生長面積減少影響其種群生長,而9-10月適宜生長面積增加,但是汛末突變同時也使得水位下降幅度的增加,導致其原生長空間被侵占,新增的生長空間不利于新的個體萌發(fā),整體壓縮了沉水植物的生存面積,導致東部湖灣沉水植物群落密度和生物量反而下降[43]。該變化規(guī)律與鄱陽湖濕地洲灘植被整體變化趨勢一致[17,44-47]。

    濕地植物的豐富性和可取食性是白鶴、鴻雁、東方白鸛等越冬候鳥選擇棲息地的重要條件。薹草群落是鄱陽湖濕地分布最廣、面積最大的植被群落類型,它的嫩芽是越冬候鳥的重要食物來源。薹草是東部湖灣重要的濕地植物,水文連通性變化導致薹草等濕生植物移向碟形湖湖心,可能會導致鄱陽湖越冬候鳥種群和棲息地發(fā)生變化,以薹草嫩芽為主要食物的雁類(如豆雁(Anserfabalis)和鴻雁)種群增加明顯[48]。每年10月主湖區(qū)適宜鳥類覓食的草洲尚未出露,淺水區(qū)的面積也相對較小,11月星子水位下降至10 m左右,洲灘和碟形湖與主湖區(qū)水文連通性逐漸變?nèi)?,水陸過渡帶面積越來越大,12月至翌年3月份,星子水位逐步下降至6 m左右,碟形湖內(nèi)沉水植被帶達到白鶴等食塊莖鳥類的適宜覓食水深范圍,大量鳥類集中覓食[49]。對棲息地質(zhì)量變化較為敏感的種群(如鶴類)覓食棲息地發(fā)生了轉(zhuǎn)移,鄱陽湖越冬白鶴原來僅在碟形湖淺水區(qū)覓食,2010年首次發(fā)現(xiàn)在草洲取食下江萎陵菜的塊莖和老鴉瓣的球莖。此外,近些年也發(fā)現(xiàn)越冬白鶴在湖區(qū)周邊藕塘、稻田覓食的現(xiàn)象[50]。

    3.3 對鄱陽湖管理的啟示

    在自然和人類活動的雙重作用下,水文連通發(fā)生變化,使?jié)竦卦械纳鷳B(tài)平衡被打破[51]。洪泛區(qū)是鄱陽湖重要的生態(tài)系統(tǒng),其水文連通性的減弱與面積的萎縮將帶來一系列的生態(tài)問題,如沉水植物的減少以及鳥類活動區(qū)域的轉(zhuǎn)變等,必須引起高度重視。水利工程建設(shè)、采砂等人類活動以及氣候變化引起鄱陽湖主湖區(qū)水位變化是東部湖灣水文連通性變化的重要影響因素。近年來,采砂活動得到較好的控制,而水利工程調(diào)度要考慮氣候變化對水位的影響,還應充分考慮在3、6月沉水植物重要生長期以及秋冬季節(jié)的水位調(diào)控,并減少水位波動幅度,以進一步保護生態(tài)平衡。

    4 結(jié)論

    水文連通性是湖泊水文情勢變化的重要組成部分,水文連通狀況及改變程度很大程度上影響了鄱陽湖碟形湖及湖汊濕地植物群落和候鳥棲息地的分布特征。本文分析1962-2019年鄱陽湖主湖區(qū)水位的變異程度,定量探究東部湖灣與主湖區(qū)的連通天數(shù)和發(fā)生時間等水文連通特征的變化,討論鄱陽湖水文連通性改變對濕地植物、候鳥覓食的潛在影響。

    1)與1962-2002年相比,2003年以后鄱陽湖主湖區(qū)枯水期水位普遍降低,持續(xù)時間延長。從月平均水位、年極端水位、年極端水位出現(xiàn)時間、高低水位頻率和持續(xù)時間、水位變化率5組32個水文指標分析,水位中高度改變指標為13個,整體改變度為40.2%,屬于中度改變,對濕生植物和沉水植物的分布、候鳥越冬覓食棲息將產(chǎn)生一定的不利影響。

    2)2003年以后東部湖灣與主湖區(qū)的水文連通發(fā)生了較為明顯的改變,非連通期明顯提前并延長。與1962-2002年比較,2003-2019年東部湖灣與主湖區(qū)連通天數(shù)減少了46.2 d,變化幅度為15.9%,其中10-12月連通天數(shù)的變化最大,變化幅度達46.8%;與主湖區(qū)的非連通期天數(shù)有大幅度增加,發(fā)生時間提前1個月且年內(nèi)時間分布較為離散,可能導致薹草和苦草等濕地植物移向湖心且生長周期提前,越冬鳥類的覓食區(qū)域轉(zhuǎn)移至主湖區(qū)的其他淺灘。

    3)水文連通性改變將破壞鄱陽湖洪泛區(qū)現(xiàn)有的生態(tài)平衡,需要引起高度重視。采砂和水利工程建設(shè)是重要的影響因素,需要加強管理,科學調(diào)控,氣候變化的影響也不容忽視。本文受篇幅和資料限制,對于東部湖灣濕地植被、候鳥研究主要基于文獻數(shù)據(jù),更加深入的研究還需要進一步調(diào)查,另外除水位變異影響之外,濕地植被還受到土壤、水動力條件、營養(yǎng)物質(zhì)[52-53]等的制約,可將其作為進一步研究的對象。對于濕地植物和候鳥對適宜生境及維持適宜生境的水文連通需求、如何調(diào)控以達到水文連通要求等方面量化研究還需要進一步加強。

    猜你喜歡
    湖灣苦草連通性
    苦草對水體氮磷吸收與釋放作用的探析
    偏序集及其相關(guān)拓撲的連通性?
    北部太湖DO、pH時空變化及其與葉綠素a的相關(guān)性分析
    洱海湖灣水生植物群落結(jié)構(gòu)及時空分布特征研究
    這里的農(nóng)民愛種草,還能賺錢
    至愛(2020年10期)2020-11-20 02:41:06
    基于近自然結(jié)構(gòu)的湖灣生態(tài)修復工程設(shè)計
    寫真地理(2020年24期)2020-10-20 08:58:35
    擬莫比烏斯映射與擬度量空間的連通性
    不同類型底棲藻對養(yǎng)殖廢水中苦草生長的影響
    河道-灘區(qū)系統(tǒng)連通性評價研究
    高穩(wěn)定被動群集車聯(lián)網(wǎng)連通性研究
    通信學報(2016年11期)2016-08-16 03:20:04
    高清日韩中文字幕在线| 麻豆成人av在线观看| 日韩欧美国产在线观看| 亚洲最大成人手机在线| 欧美日韩中文字幕国产精品一区二区三区| 此物有八面人人有两片| 老司机深夜福利视频在线观看| 老司机午夜福利在线观看视频| www.www免费av| 噜噜噜噜噜久久久久久91| 国产精品无大码| 老女人水多毛片| 国产精品美女特级片免费视频播放器| 淫秽高清视频在线观看| 一a级毛片在线观看| 国产精品1区2区在线观看.| 真实男女啪啪啪动态图| 国产视频内射| 一边摸一边抽搐一进一小说| 精品久久国产蜜桃| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 一区二区三区免费毛片| 亚洲国产欧洲综合997久久,| 国产精品福利在线免费观看| 亚洲色图av天堂| 美女被艹到高潮喷水动态| 蜜桃亚洲精品一区二区三区| 久久国内精品自在自线图片| 搡老妇女老女人老熟妇| 亚洲av美国av| 少妇人妻一区二区三区视频| 久久久精品大字幕| 可以在线观看的亚洲视频| 两个人视频免费观看高清| 亚洲18禁久久av| 久久久久久久久大av| 少妇高潮的动态图| 精品久久久久久久久久免费视频| 性欧美人与动物交配| 免费观看的影片在线观看| 俄罗斯特黄特色一大片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲一区高清亚洲精品| 性色avwww在线观看| 免费看美女性在线毛片视频| 久久久久久久久大av| av在线观看视频网站免费| 亚洲精品一卡2卡三卡4卡5卡| 国产在线男女| 成人二区视频| 成人无遮挡网站| 啪啪无遮挡十八禁网站| 久久精品国产自在天天线| 亚洲 国产 在线| 九九在线视频观看精品| 中国美白少妇内射xxxbb| 又黄又爽又刺激的免费视频.| 最后的刺客免费高清国语| 久久中文看片网| 日韩精品中文字幕看吧| 少妇的逼水好多| 精品一区二区三区av网在线观看| 特大巨黑吊av在线直播| 乱系列少妇在线播放| 最后的刺客免费高清国语| 天美传媒精品一区二区| 精品久久久久久久久av| 久久久久久伊人网av| av天堂在线播放| 亚洲成人久久性| 黄色日韩在线| 伦精品一区二区三区| 国内精品久久久久精免费| 日韩欧美在线二视频| 综合色av麻豆| 在线观看一区二区三区| 日韩精品青青久久久久久| 一个人看的www免费观看视频| 色综合婷婷激情| 九色国产91popny在线| 成人亚洲精品av一区二区| 精品无人区乱码1区二区| 精品久久久久久久末码| 综合色av麻豆| 久久人人爽人人爽人人片va| 3wmmmm亚洲av在线观看| 成人毛片a级毛片在线播放| 在线观看舔阴道视频| 国产色婷婷99| av黄色大香蕉| 内地一区二区视频在线| 国内精品久久久久久久电影| 可以在线观看的亚洲视频| 免费观看的影片在线观看| 亚洲精华国产精华精| 午夜亚洲福利在线播放| 3wmmmm亚洲av在线观看| 麻豆久久精品国产亚洲av| 麻豆成人av在线观看| 国产白丝娇喘喷水9色精品| 在线播放无遮挡| 麻豆一二三区av精品| 成人综合一区亚洲| 男女那种视频在线观看| 免费黄网站久久成人精品| 成人亚洲精品av一区二区| 18禁在线播放成人免费| av在线蜜桃| 永久网站在线| 999久久久精品免费观看国产| 免费av不卡在线播放| 欧美日本视频| 赤兔流量卡办理| 亚洲三级黄色毛片| 亚洲在线自拍视频| 黄色丝袜av网址大全| 色av中文字幕| 蜜桃久久精品国产亚洲av| 亚洲久久久久久中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产人妻一区二区三区在| 久久久久九九精品影院| videossex国产| 国产亚洲精品久久久com| 亚洲av一区综合| 哪里可以看免费的av片| 欧美精品国产亚洲| 亚洲av成人av| 日本精品一区二区三区蜜桃| 丰满人妻一区二区三区视频av| 一进一出好大好爽视频| 久久99热这里只有精品18| 丰满的人妻完整版| 他把我摸到了高潮在线观看| 欧美bdsm另类| 最新中文字幕久久久久| 男人舔女人下体高潮全视频| 久久精品夜夜夜夜夜久久蜜豆| av国产免费在线观看| 亚洲国产精品sss在线观看| 最好的美女福利视频网| 99国产极品粉嫩在线观看| av在线观看视频网站免费| 欧美日本亚洲视频在线播放| 嫩草影视91久久| 欧美日本亚洲视频在线播放| 色综合色国产| 两性午夜刺激爽爽歪歪视频在线观看| 韩国av一区二区三区四区| 看片在线看免费视频| 午夜视频国产福利| 看黄色毛片网站| 欧美最黄视频在线播放免费| 国产精品电影一区二区三区| 九九爱精品视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产av一区在线观看免费| 又爽又黄a免费视频| 国产aⅴ精品一区二区三区波| 高清在线国产一区| 成人美女网站在线观看视频| 精品久久久噜噜| a级一级毛片免费在线观看| 免费人成在线观看视频色| 美女被艹到高潮喷水动态| 欧美日本视频| 国产成人一区二区在线| 搞女人的毛片| 欧美激情久久久久久爽电影| 国产精品一区二区三区四区免费观看 | 99久久久亚洲精品蜜臀av| 伦理电影大哥的女人| 黄色日韩在线| 国内精品美女久久久久久| 91久久精品电影网| 欧美又色又爽又黄视频| 久久久午夜欧美精品| 深夜a级毛片| 99久久精品热视频| 亚洲一级一片aⅴ在线观看| 国产欧美日韩精品亚洲av| 88av欧美| 亚洲国产欧洲综合997久久,| 久久精品综合一区二区三区| 深爱激情五月婷婷| 国产高清有码在线观看视频| 国产精品久久视频播放| 男人的好看免费观看在线视频| 毛片女人毛片| 在线免费十八禁| 美女 人体艺术 gogo| x7x7x7水蜜桃| 国产在视频线在精品| 色综合婷婷激情| 少妇熟女aⅴ在线视频| 欧美最黄视频在线播放免费| 日本一本二区三区精品| 亚洲av中文字字幕乱码综合| 欧美+日韩+精品| 免费无遮挡裸体视频| 国产午夜福利久久久久久| 亚洲精品亚洲一区二区| 色在线成人网| 丰满的人妻完整版| 国产精品三级大全| 日韩强制内射视频| 亚洲熟妇熟女久久| a级毛片a级免费在线| 欧美激情在线99| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久com| 中文字幕精品亚洲无线码一区| 91麻豆av在线| 男女那种视频在线观看| 亚洲午夜理论影院| 少妇熟女aⅴ在线视频| 国产熟女欧美一区二区| 国产高清视频在线播放一区| 人人妻人人澡欧美一区二区| 黄色女人牲交| 国产一区二区激情短视频| 欧美在线一区亚洲| 国产精品一区二区三区四区久久| 亚洲综合色惰| 99国产精品一区二区蜜桃av| 中国美女看黄片| 久久午夜亚洲精品久久| 国产成人影院久久av| 成人性生交大片免费视频hd| 免费看日本二区| 国产在线精品亚洲第一网站| 亚洲成av人片在线播放无| 老熟妇乱子伦视频在线观看| 不卡一级毛片| 精品人妻视频免费看| 久久久久久国产a免费观看| 亚洲av一区综合| 亚洲va在线va天堂va国产| 亚洲欧美激情综合另类| 亚洲精品乱码久久久v下载方式| 久久精品国产亚洲av天美| 国产在线男女| 国产精品1区2区在线观看.| 日本黄大片高清| 午夜福利18| 美女高潮喷水抽搐中文字幕| 制服丝袜大香蕉在线| 高清毛片免费观看视频网站| bbb黄色大片| 亚洲图色成人| 成人欧美大片| av天堂中文字幕网| 999久久久精品免费观看国产| 亚洲专区国产一区二区| 午夜福利高清视频| 国内揄拍国产精品人妻在线| 长腿黑丝高跟| 免费观看人在逋| 好男人在线观看高清免费视频| 女人被狂操c到高潮| 亚洲在线观看片| 久久精品国产亚洲av涩爱 | 国产一区二区三区在线臀色熟女| 精品久久久久久久末码| 亚洲人成伊人成综合网2020| 中文字幕人妻熟人妻熟丝袜美| 欧美高清成人免费视频www| 看免费成人av毛片| av在线蜜桃| 久久亚洲精品不卡| 狂野欧美激情性xxxx在线观看| 亚洲五月天丁香| 韩国av在线不卡| 级片在线观看| 免费看光身美女| 久久热精品热| 精品一区二区免费观看| 一级毛片久久久久久久久女| 1024手机看黄色片| 热99re8久久精品国产| 两人在一起打扑克的视频| 日韩欧美国产在线观看| 国产亚洲精品av在线| 国内精品美女久久久久久| 97人妻精品一区二区三区麻豆| 国产精品久久久久久精品电影| av在线天堂中文字幕| 久久国内精品自在自线图片| 婷婷丁香在线五月| 欧美潮喷喷水| 特级一级黄色大片| 全区人妻精品视频| 少妇丰满av| 亚洲最大成人av| 男人的好看免费观看在线视频| 久久精品国产99精品国产亚洲性色| av在线老鸭窝| 免费无遮挡裸体视频| 最好的美女福利视频网| 国产精品一区二区三区四区久久| 男女视频在线观看网站免费| xxxwww97欧美| 日韩一区二区视频免费看| 国产亚洲精品久久久久久毛片| 欧美日韩精品成人综合77777| 丰满乱子伦码专区| 国产精华一区二区三区| 蜜桃久久精品国产亚洲av| 中文字幕av在线有码专区| 中文亚洲av片在线观看爽| 91av网一区二区| 亚州av有码| 国产av不卡久久| 夜夜爽天天搞| 99热这里只有是精品在线观看| 欧美一级a爱片免费观看看| 波野结衣二区三区在线| 久久国产乱子免费精品| 久久精品久久久久久噜噜老黄 | 久久精品国产自在天天线| 色综合站精品国产| 国产亚洲av嫩草精品影院| 欧美日韩亚洲国产一区二区在线观看| 看黄色毛片网站| 国内精品美女久久久久久| 99久久精品热视频| av视频在线观看入口| 女人被狂操c到高潮| 日日摸夜夜添夜夜添小说| 在线a可以看的网站| 少妇丰满av| 观看美女的网站| 国产麻豆成人av免费视频| 国产成人a区在线观看| 国产精品98久久久久久宅男小说| 亚洲精品日韩av片在线观看| 日本a在线网址| 国产视频内射| av天堂中文字幕网| 亚洲精品成人久久久久久| 天堂动漫精品| 嫁个100分男人电影在线观看| 成人亚洲精品av一区二区| 久久久久免费精品人妻一区二区| 亚洲经典国产精华液单| 两个人的视频大全免费| 久久99热这里只有精品18| 搡老妇女老女人老熟妇| 两人在一起打扑克的视频| 午夜视频国产福利| 欧美日韩黄片免| 搡女人真爽免费视频火全软件 | 韩国av一区二区三区四区| 91久久精品电影网| 天堂动漫精品| 国产精品嫩草影院av在线观看 | 亚洲,欧美,日韩| 日韩人妻高清精品专区| 深夜精品福利| 日韩欧美精品免费久久| 亚洲avbb在线观看| 久久精品国产亚洲av涩爱 | 欧美成人性av电影在线观看| 午夜激情福利司机影院| 99热精品在线国产| 国产视频一区二区在线看| 草草在线视频免费看| 女人十人毛片免费观看3o分钟| 亚洲国产精品成人综合色| 午夜福利视频1000在线观看| 热99re8久久精品国产| 欧美另类亚洲清纯唯美| 国产精品国产高清国产av| 搡老岳熟女国产| 欧美黑人欧美精品刺激| 美女cb高潮喷水在线观看| 又爽又黄a免费视频| 国产老妇女一区| 久久精品国产亚洲av香蕉五月| 国产精品女同一区二区软件 | 婷婷精品国产亚洲av| 午夜福利欧美成人| 国产高清有码在线观看视频| 国产单亲对白刺激| 日韩国内少妇激情av| 亚洲自拍偷在线| 偷拍熟女少妇极品色| 日本黄色片子视频| 在线播放无遮挡| 18禁黄网站禁片免费观看直播| 亚洲av第一区精品v没综合| 性欧美人与动物交配| 欧美精品啪啪一区二区三区| 少妇丰满av| 亚洲五月天丁香| 色尼玛亚洲综合影院| av福利片在线观看| 国产精品国产三级国产av玫瑰| 亚洲av免费高清在线观看| 亚洲三级黄色毛片| 22中文网久久字幕| 女生性感内裤真人,穿戴方法视频| 成人性生交大片免费视频hd| 亚洲中文日韩欧美视频| 亚洲电影在线观看av| 亚洲18禁久久av| 国产成人av教育| 欧美日韩国产亚洲二区| 午夜福利18| 最近最新免费中文字幕在线| 中文字幕高清在线视频| www.色视频.com| 日本a在线网址| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线播放欧美日韩| 国产精品人妻久久久影院| 国产美女午夜福利| 午夜视频国产福利| 国产精品久久久久久亚洲av鲁大| 国产精品精品国产色婷婷| 亚洲avbb在线观看| 亚洲国产高清在线一区二区三| 亚洲内射少妇av| 97超级碰碰碰精品色视频在线观看| 国产免费男女视频| 草草在线视频免费看| 禁无遮挡网站| 日韩欧美在线乱码| 欧美日韩精品成人综合77777| 啪啪无遮挡十八禁网站| 尤物成人国产欧美一区二区三区| 动漫黄色视频在线观看| 色哟哟·www| 99久国产av精品| 国产视频一区二区在线看| 长腿黑丝高跟| 久久精品国产99精品国产亚洲性色| 网址你懂的国产日韩在线| 88av欧美| 桃色一区二区三区在线观看| 成年女人永久免费观看视频| 久久精品国产亚洲网站| 看免费成人av毛片| 亚洲精品色激情综合| 久久精品国产亚洲网站| 白带黄色成豆腐渣| 亚洲av美国av| 女人被狂操c到高潮| 国产视频一区二区在线看| 色综合站精品国产| 特大巨黑吊av在线直播| 久久精品91蜜桃| 亚洲无线在线观看| 精品午夜福利视频在线观看一区| 一区二区三区激情视频| 久久精品综合一区二区三区| 极品教师在线视频| 色哟哟哟哟哟哟| 国产黄色小视频在线观看| 亚洲图色成人| 97碰自拍视频| 欧美xxxx黑人xx丫x性爽| 亚洲中文日韩欧美视频| 99久久精品热视频| 免费在线观看成人毛片| 亚洲av二区三区四区| 欧美性猛交黑人性爽| 99久久成人亚洲精品观看| 成人三级黄色视频| 亚洲成人久久爱视频| 看十八女毛片水多多多| 美女黄网站色视频| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| АⅤ资源中文在线天堂| 又爽又黄无遮挡网站| 久久精品91蜜桃| 亚洲七黄色美女视频| 国产高清三级在线| 2021天堂中文幕一二区在线观| a在线观看视频网站| 国产免费男女视频| 国产美女午夜福利| 日本熟妇午夜| 亚洲国产欧美人成| 九色国产91popny在线| av国产免费在线观看| 在线免费观看不下载黄p国产 | 国产av不卡久久| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 日日摸夜夜添夜夜添av毛片 | 亚洲人成网站在线播| 欧美一区二区精品小视频在线| 夜夜爽天天搞| 午夜日韩欧美国产| 国产精品久久久久久精品电影| 久久久久性生活片| 动漫黄色视频在线观看| 99国产精品一区二区蜜桃av| 国产精品野战在线观看| 国产蜜桃级精品一区二区三区| 久久久久性生活片| 国产午夜福利久久久久久| 一区二区三区四区激情视频 | 美女黄网站色视频| 久久草成人影院| 在线观看av片永久免费下载| 中国美女看黄片| 老司机深夜福利视频在线观看| 免费高清视频大片| 日本在线视频免费播放| 午夜影院日韩av| 亚洲av免费高清在线观看| 成人国产综合亚洲| 国产一区二区三区在线臀色熟女| 成人美女网站在线观看视频| avwww免费| 成人鲁丝片一二三区免费| 欧美一区二区精品小视频在线| 国产精品久久视频播放| 特级一级黄色大片| 丰满乱子伦码专区| 国产精品一区二区免费欧美| 国产不卡一卡二| 国产亚洲欧美98| 久久热精品热| 免费观看的影片在线观看| 亚洲色图av天堂| 免费搜索国产男女视频| 久久精品国产亚洲av涩爱 | 精品一区二区三区av网在线观看| 一a级毛片在线观看| 亚洲成a人片在线一区二区| 22中文网久久字幕| 免费观看的影片在线观看| 白带黄色成豆腐渣| 嫩草影院新地址| 亚洲成av人片在线播放无| 变态另类丝袜制服| 又黄又爽又刺激的免费视频.| 久久精品夜夜夜夜夜久久蜜豆| 夜夜夜夜夜久久久久| 国产爱豆传媒在线观看| 男插女下体视频免费在线播放| 国产乱人伦免费视频| 校园人妻丝袜中文字幕| 欧美高清成人免费视频www| 亚洲熟妇熟女久久| 又黄又爽又刺激的免费视频.| 日本 av在线| 欧美成人免费av一区二区三区| 国产av在哪里看| 18禁黄网站禁片午夜丰满| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 国产白丝娇喘喷水9色精品| 69人妻影院| 国产精品99久久久久久久久| 亚洲精品国产成人久久av| av专区在线播放| 十八禁网站免费在线| 观看免费一级毛片| 亚洲综合色惰| 日本三级黄在线观看| 午夜a级毛片| 成人亚洲精品av一区二区| 亚洲精品日韩av片在线观看| 桃色一区二区三区在线观看| 此物有八面人人有两片| 简卡轻食公司| 亚洲成人久久性| 99在线视频只有这里精品首页| 国产午夜福利久久久久久| 99riav亚洲国产免费| 欧美日韩综合久久久久久 | 成人鲁丝片一二三区免费| 亚洲 国产 在线| 亚洲成人久久爱视频| 悠悠久久av| 亚洲欧美日韩高清专用| 亚洲第一区二区三区不卡| www.www免费av| 亚洲经典国产精华液单| 69av精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 国内毛片毛片毛片毛片毛片| 亚洲国产日韩欧美精品在线观看| 久久精品综合一区二区三区| 午夜日韩欧美国产| 久久久久久大精品| av在线老鸭窝| 哪里可以看免费的av片| 日韩在线高清观看一区二区三区 | 午夜福利在线观看吧| 国产白丝娇喘喷水9色精品| 色播亚洲综合网| 精品福利观看| 嫩草影院入口| 免费av不卡在线播放| 久久久久国内视频| 久久香蕉精品热| 午夜福利视频1000在线观看| 国产精品久久久久久久电影| 男人的好看免费观看在线视频| 精品不卡国产一区二区三区| 性色avwww在线观看| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 我的老师免费观看完整版| 在线免费十八禁| 亚洲 国产 在线| 韩国av在线不卡|