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

    三峽水庫枯水期補(bǔ)水調(diào)度對洞庭湖越冬白鶴(Grus leucogeranus)攝食棲息地的影響*

    2022-06-30 07:06:54戴凌全湯正陽戴會超任玉峰常曼琪毛勁喬
    湖泊科學(xué) 2022年4期
    關(guān)鍵詞:三峽水庫白鶴洞庭湖

    戴凌全,王 煜,湯正陽,戴會超,,任玉峰,常曼琪,吳 倩,毛勁喬

    (1:中國長江三峽集團(tuán)有限公司智慧長江與水電科學(xué)湖北省重點(diǎn)實(shí)驗室,宜昌 443133)(2:中國長江三峽集團(tuán)有限公司長江生態(tài)環(huán)境工程研究中心,北京 100038)(3:三峽大學(xué)三峽水庫生態(tài)系統(tǒng)湖北省野外科學(xué)觀測研究站,宜昌 443002)(4:河海大學(xué)水利水電學(xué)院,南京 210098)

    洞庭湖處于東亞-澳大利西亞遷飛區(qū)上,廣袤的濕地是越冬水鳥的理想棲息地,可為水鳥提供充足的食物資源以保證其生存的基本條件[1-2].區(qū)域內(nèi)水鳥具有種類多、數(shù)量大、珍稀瀕危程度高等特點(diǎn),目前已知越冬水鳥共有7目12科53種,其中國家一級保護(hù)鳥類3種,二級保護(hù)鳥類4 種[3].鶴類是洞庭湖區(qū)重要的瀕危越冬水鳥,目前存在白鶴(Grusleucogeranus)、白頭鶴(G.monacha)、白枕鶴(G.vipio)、灰鶴(G.grus)和蓑羽鶴(G.virgo)5種鶴類.對以沉水植物為主要食物來源的鶴類而言,湖泊水深是影響鶴類攝食可及性的決定性因子之一.鶴類作為大型水禽,其攝食對水深變化非常敏感,高水深下鶴類的攝食范圍、頻率將有不同程度的降低,而過低水深條件下鶴類喜好的食物易枯死,且泥灘表面基質(zhì)干燥也妨礙鶴類正常取食.此外,水深變化節(jié)律直接影響到洲灘濕地淹沒和出露日期及持續(xù)時間,對鶴類食物資源的產(chǎn)量和分布產(chǎn)生直接影響,進(jìn)而對鶴類的數(shù)量和分布構(gòu)成較大影響[4].進(jìn)入21世紀(jì)以來,受氣候變化、水利工程運(yùn)行、楊樹(Populusdeltoides)大面積種植等多重因素的影響,洞庭湖水文節(jié)律發(fā)生了明顯變化[5-9],枯水期湖泊水位下降、水域面積縮小、調(diào)蓄能力減弱,對越冬水鳥生境造成不利影響[10-12],因此采取有效措施保障越冬水鳥水位需求已成為棲息地保護(hù)工作中的重點(diǎn)和難點(diǎn)問題[13].維持洞庭湖適宜水位可采取工程方案和非工程方案,其中工程方案包括城陵磯建閘、三口建閘等[14-15],但工程方案是否對洞庭湖生態(tài)產(chǎn)生不利影響,目前仍在論證階段[16-17].非工程方案包括保護(hù)區(qū)核心區(qū)楊樹的清理,截止2020年底核心區(qū)楊樹清理完畢,同時清理跡地的生態(tài)修復(fù)工作也正在逐步開展,主要包括封溝育洲、恢復(fù)洲灘植被、種植沉水植物等.此外,隨著三峽水庫調(diào)度規(guī)程的逐步完善,通過水庫優(yōu)化調(diào)控對洞庭湖補(bǔ)水的非工程方案正受到各方的高度重視,然而由于三峽水庫不同出庫流量條件下洞庭湖珍稀水鳥攝食棲息地面積演變規(guī)律尚不明確,因此,量化三峽水庫不同出庫流量對水鳥攝食棲息生境的影響是開展補(bǔ)水調(diào)度的關(guān)鍵環(huán)節(jié).

    洞庭湖水鳥越冬棲息地主要分為兩類:休息棲息地和攝食棲息地[18],水鳥找到合適的生存環(huán)境時,不僅要考慮環(huán)境是否適合休息嬉戲,更要考慮能否方便的獲得食物資源[19].國內(nèi)外學(xué)者對于休息棲息地開展了卓有成效的研究,多通過遙感解譯獲得不同水位下的生境面積,如關(guān)蕾等以衛(wèi)星影像為基礎(chǔ),運(yùn)用遙感技術(shù)解譯洞庭湖內(nèi)泥灘、草洲等水鳥休息棲息地分布狀況[20],周延等采用衛(wèi)星圖片提取東洞庭湖不同時期濕地景觀格局?jǐn)?shù)據(jù),分析了鶴類種群動態(tài)與濕地景觀格局之間的關(guān)聯(lián)[21],梁婕等在洞庭湖適宜生態(tài)水位的基礎(chǔ)上,基于中分辨率成像光譜儀(MODIS)影像數(shù)據(jù)分析了洞庭湖水體濕地和洲灘濕地的面積變化趨勢[22].洞庭湖濕地為越冬水鳥提供了豐富的食物和休息場所,從食性分析來看,鶴類的主要食物為沉水植物的冬芽,湖泊水深是鶴類攝食可及性的決定性因子之一[23].棲息地面積可以反映生物平面生存空間,但無法表征不同區(qū)域水深空間分布特征及獲取食物的概率,因此在評判湖泊棲息地質(zhì)量時,不僅要關(guān)注面積的大小,還應(yīng)考慮不同生物對水深的需求[24-26].相比于遙感解譯的不足,考慮目標(biāo)物種攝食對水深需求的物理棲息地模型(PHABSIM模型)可獲得不同水文條件下大量的生境水動力特征數(shù)據(jù),從而提高棲息地生境變化規(guī)律的確定性,但目前采用物理棲息地模型從水鳥攝食可及性所需的水動力條件角度確定洞庭湖候鳥棲息地面積方面的成果還較少.

    本文以越冬白鶴重要的棲息地洞庭湖為研究區(qū)域,構(gòu)建白鶴攝食水動力生境需求的物理棲息地模型,闡明三峽水庫枯水期不同補(bǔ)水流量下洞庭湖白鶴攝食棲息地面積變化規(guī)律.研究成果有助于更深入了解白鶴的攝食行為和水動力因子的關(guān)系,可為改善洞庭湖白鶴攝食棲息地的三峽水庫枯水期補(bǔ)水調(diào)度方案的制定提供理論基礎(chǔ),進(jìn)而為珍稀候鳥的保護(hù)和管理提供參考資料.

    1 研究區(qū)域概況

    洞庭湖屬于季節(jié)性湖泊,4月水位開始上漲,7-8月達(dá)到最高,最大水域面積超過2500 km2,11月-次年3月為枯水期,最小水域面積不足500 km2,城陵磯站多年最大水位變幅達(dá)18.77 m[27].當(dāng)城陵磯水位降至28.0 m時,自西南向東北逐漸形成一個傾斜的水面,當(dāng)城陵磯水位低于24.0 m時,洞庭湖形成西洞庭湖、南洞庭湖(包含橫嶺湖)、東洞庭湖3個不同的湖區(qū)[28],研究區(qū)域如圖1所示.

    圖1 洞庭湖位置及自然保護(hù)區(qū)核心區(qū)示意

    洞庭湖年內(nèi)巨大的水位差異為洲灘濕地的發(fā)育提供了良好的水文條件,為越冬水鳥提供了豐富的食物來源,每年冬季有大量來自西伯利亞和蒙古等地的水鳥來此越冬[29].為維持洞庭湖適宜的濕地面積,共建立了4個不同級別的自然保護(hù)區(qū)[30]:東洞庭湖國家級自然保護(hù)區(qū)、西洞庭湖國家級自然保護(hù)區(qū)、南洞庭湖省級自然保護(hù)區(qū)和橫嶺湖省級自然保護(hù)區(qū)(圖1),保護(hù)對象均以濕地生態(tài)系統(tǒng)和珍稀水禽為主,保護(hù)區(qū)基本情況如表1所示.

    表1 洞庭湖自然保護(hù)區(qū)基本情況

    2 材料與方法

    2.1 物理棲息地模型

    物理棲息地模型(physical habitat simulation model,PHABSIM模型)分為水動力模型(hydrodynamic model)和棲息地模型(habitat model)兩部分,水動力模型計算棲息地水位、水深及流速分布,棲息地模型將目標(biāo)物種對棲息地的偏好性和生境因子之間的關(guān)系進(jìn)行量化,根據(jù)不同流量下計算得到的水動力分布特征找出對應(yīng)的棲息地生境適宜度指數(shù)(habitat suitability index,HSI),賦予各自權(quán)重后計算出可使用的棲息地面積.與其他方法相比,PHABSIM模型考慮了生物本身對物理生境的要求,被認(rèn)為是生態(tài)流量計算中最具科學(xué)依據(jù)的方法之一.具體步驟如下:

    1)收集目標(biāo)物種棲息時間、地點(diǎn)、數(shù)量以及對應(yīng)的水動力條件;2)遴選關(guān)鍵生境因子,將棲息地適宜性指標(biāo)與目標(biāo)物種生境的水動力影響因子關(guān)聯(lián),量化目標(biāo)物種對棲息地水動力因素的偏好程度,繪制目標(biāo)物種水動力生境適宜性曲線;3)建立目標(biāo)物種棲息地的水動力學(xué)模型,以不同流量為邊界條件模擬目標(biāo)物種棲息地水動力分布的空間特征;4)將目標(biāo)物種生境適宜性曲線輸入棲息地模型,得到影響因子的棲息地適宜度指數(shù),結(jié)合水動力模型模擬的水動力分布特征計算棲息地的加權(quán)可利用面積,計算方法如下:

    (1)

    式中,WUA為目標(biāo)物種的加權(quán)可利用面積,km2;n為網(wǎng)格單元的數(shù)量;HSIi是第i個網(wǎng)格單元的棲息地生境適宜度指數(shù);Hi是第i個網(wǎng)格單元水體的水深,m;Ai是第i個網(wǎng)格單元水域的面積,km2.

    2.2 白鶴攝食生境適宜性曲線的構(gòu)建

    攝食生境適宜度曲線的構(gòu)建采用棲息地模型,適宜度曲線以生境因子的數(shù)值為橫坐標(biāo),以目標(biāo)物種對此生境因子的適宜度為縱坐標(biāo),建立目標(biāo)物種對生境因子的偏好與生境因子之間關(guān)系的連續(xù)曲線,定量描述物種在該條件下的生存質(zhì)量.目標(biāo)物種對生境因子的偏好設(shè)定在0~1之間,曲線峰值代表生物對該因子的最適宜值[31].白鶴又名西伯利亞鶴(Siberian crane),全球數(shù)量僅有3500~4000只,已被世界自然保護(hù)聯(lián)盟(IUCN)紅色名錄列為極瀕等級.白鶴僅在西伯利亞靠近北極圈附近繁殖[32],在長江中下游及以南地區(qū)越冬,其中在中國越冬的白鶴數(shù)量約占全球的90%[6,33],長江中下游地區(qū)的通江湖泊鄱陽湖是白鶴主要的越冬棲息地,而洞庭湖也是白鶴重要的棲息地之一[34-35].根據(jù)自然保護(hù)區(qū)管理局、公開權(quán)威報道及已有文獻(xiàn)數(shù)據(jù)統(tǒng)計[36],每年在洞庭湖均能發(fā)現(xiàn)不同數(shù)量的白鶴(表2),根據(jù)對白鶴出現(xiàn)時間的統(tǒng)計,白鶴越冬時間為10月-次年4月,一般在1月中旬達(dá)到峰值,因此研究時間集中在1月中旬.

    表2 洞庭湖2006-2020年白鶴數(shù)量變化及分布

    白鶴攝食棲息地選擇與水深、食物數(shù)量密切相關(guān),其中水深不但決定攝食棲息地是否可以被白鶴利用,還能影響白鶴的取食行為和取食能耗,因此選擇水深作為白鶴攝食的關(guān)鍵水動力因子.根據(jù)對白鶴食性的觀察和分析,沉水植物苦草(Vallisnerianatans)的冬芽是白鶴越冬的主要食物來源.苦草隸屬于沉水植物的水鱉科,是多年生、蓮座型沉水草本植物.苦草是洞庭湖濕地分布范圍最廣、面積最大的典型沉水植物種類之一,主要分布在淺水洼地[37].白鶴取食時通過喙部高頻率挖掘底泥尋覓冬芽,對于基質(zhì)干燥的區(qū)域,白鶴需消耗太多體力或無法取食.根據(jù)白鶴身體各部位的平均長度和取食時喙、頭、頸沒入水中的比例估算白鶴取食適宜的水深是20~30 cm,當(dāng)水深超過30 cm時取食效率與水深呈反比,當(dāng)水深超過60 cm時白鶴發(fā)現(xiàn)食物非常困難且無法落腳[24],根據(jù)白鶴攝食水動力條件[38-39],制定洞庭湖白鶴攝食的水深適宜度曲線(圖2).

    圖2 洞庭湖白鶴攝食水深-適宜度曲線

    2.3 白鶴棲息地水動力模型的構(gòu)建與驗證

    為分析白鶴棲息地水動力分布特征,基于MIKE Flood構(gòu)建了涵蓋長江干流、三口河系、洞庭湖及其四水尾閭河段的江湖一體化耦合水動力模型.長江干流、三口河系、四水尾閭河段采用一維模型,洞庭湖及江湖交匯區(qū)采用二維模型,長江干流河段5~10 km布設(shè)一個斷面,三口河系及四水尾閭河段平均10 km布設(shè)一個斷面.根據(jù)洞庭湖歷史水文資料確定湖區(qū)最高水位淹沒邊界,據(jù)此設(shè)定二維模型的計算范圍.根據(jù)邊界和地形數(shù)據(jù)劃分網(wǎng)格,洞庭湖枯水期呈現(xiàn)明顯的“河相”,水流只在河道深槽中流動,大部分洲灘濕地顯露[40],為了精確反映洞庭湖特殊的地形,提高模型的精度,二維模型采用非結(jié)構(gòu)不規(guī)則的三角形網(wǎng)格對湖區(qū)進(jìn)行網(wǎng)格劃分,網(wǎng)格大小根據(jù)地形的空間分布等確定,對湖區(qū)主河槽網(wǎng)格進(jìn)行加密.地形邊界采用長江干流、三口、四水尾閭河道地形(2018年)和洞庭湖1∶5000 實(shí)測水下地形(2011年),洞庭湖地形及網(wǎng)格劃分示意圖如圖 3所示,可見模型所用網(wǎng)格清晰地描述了河槽(藍(lán)色)、灘地(淺黃色)、島嶼(紅色)等地形的變化.根據(jù)上述原則并兼顧計算時間,網(wǎng)格按最大面積不超過0.1 km2生成,最終生成網(wǎng)格數(shù)量49164個,節(jié)點(diǎn)數(shù)量27189個.

    圖3 洞庭湖地形及網(wǎng)格劃分示意圖

    由于研究時段集中在枯水期1月份,因此采用2020年1月1日-31日實(shí)測的長江干流及洞庭湖水位數(shù)據(jù)對水動力模型的參數(shù)進(jìn)行率定和驗證.模型邊界包括水文邊界和地形邊界,其中水文邊界包括三峽水庫出庫流量、清江流量(高壩洲站)、四水入湖控制站流量(湘潭、桃江、桃源和石龜山),上游邊界采用各水文站點(diǎn)日平均流量過程(圖4),下游邊界條件采用2010-2020年螺山站水位-流量關(guān)系曲線(圖 5).因此只需給定三峽水庫下泄流量、清江流量、四水流量即可計算出整個區(qū)域內(nèi)部的水位、流量、流速等水動力參數(shù)的時空變化過程.綜合考慮模型穩(wěn)定性及計算效率等因素,設(shè)定一維模型時間步長為10 s,二維模型最大計算步長為5 s、最小計算步長為0.1 s,步長可根據(jù)流量與地形條件自動調(diào)整.

    圖4 上游邊界流量過程

    圖5 螺山站水位-流量關(guān)系曲線

    模型的主要參數(shù)為曼寧糙率系數(shù),糙率是水動力學(xué)模型最重要且較為敏感的計算參數(shù),反映了河湖底部及岸坡、湖區(qū)地表影響水流阻力的綜合系數(shù).以長江干流及洞庭湖主要站點(diǎn)水位為特征值率定模型,長江干流選取宜昌、枝城、沙市、監(jiān)利為率定點(diǎn),洞庭湖選取城陵磯、鹿角、楊柳潭、南咀為率定點(diǎn)(圖3).主要水文站點(diǎn)水位率定結(jié)果如圖6所示,從圖可以看出各站點(diǎn)水位計算結(jié)果與實(shí)測水位基本一致,一維模型計算的長江干流4個水文站的水位模擬結(jié)果與實(shí)測水位吻合較好,計算值與實(shí)測值的平均絕對誤差(mean absolute error,MAE)為0.18 m,最大絕對誤差為0.29 m,相對誤差(relative error,RE)不超過5.0%.洞庭湖出口城陵磯站和西洞庭湖南咀站模擬效果較好,而東洞庭湖鹿角站、南洞庭湖楊柳潭站模擬效果略差,分析其

    圖6 長江干流及洞庭湖水文站點(diǎn)實(shí)測水位與計算水位對比

    原因是城陵磯水位受長江干流水位影響明顯,干流流量在枯水期保持穩(wěn)定,南咀站水位主要受澧水及三口之一的松滋河影響,在枯水期流量較小且穩(wěn)定,而洞庭湖內(nèi)部站點(diǎn)受四水影響更為復(fù)雜,因此離干流較遠(yuǎn)的鹿角站和楊柳潭站模擬效果相對較差,但相對誤差也控制在8.0%以內(nèi).總體來看所建立的水動力模型的計算誤差均在可接受范圍內(nèi),滿足計算要求.通過率定和驗證,河道糙率取值在0.023~0.029之間,湖底糙率在0.029~0.041之間.

    3 結(jié)果與討論

    根據(jù)《三峽(正常運(yùn)行期)-葛洲壩水利樞紐梯級調(diào)度規(guī)程》(2019年修訂版)規(guī)定:1月三峽水庫增大下泄流量,水庫水位從正常蓄水位175 m開始下降,6 月上旬降至防洪汛限水位145 m.為分析三峽水庫不同出庫流量對洞庭湖水位及白鶴攝食棲息地的影響,分別計算水庫蓄水前后不同代表性出庫流量條件下的棲息地加權(quán)可利用面積.根據(jù)1988-2002年三峽水庫蓄水前宜昌站1月中旬流量統(tǒng)計:最小流量為3740 m3/s(1997年),最大流量為5280 m3/s(2001年),平均流量為4300 m3/s.2010年三峽水庫首次蓄水至正常蓄水位175 m,根據(jù)2011-2020年出庫流量統(tǒng)計,1月中旬出庫流量最小值為5960 m3/s(2015年),最大值為9227.5 m3/s(2019年),平均出庫流量為7000 m3/s,據(jù)此設(shè)定三峽水庫運(yùn)行后可能出庫流量范圍為5500~9500 m3/s,并以500 m3/s為流量間隔作為模擬出庫流量.根據(jù)1988-2002年清江及四水1月中旬流量統(tǒng)計,清江、湘江、資水、沅江、澧水1月中旬最小流量分別為37.52、636.14、188.74、520.53、158.92 m3/s,平均流量分別為77.18、1000.80、446.40、754.60、228.50 m3/s,最大流量分別為95.42、1630.15、786.57、1213.62、295.97 m3/s.取3種典型工況進(jìn)行對比展示:(1)長江干流宜昌站流量為3740 m3/s,代表三峽水庫運(yùn)行前1月中旬最小出庫流量,對應(yīng)的清江及四水流量取1月中旬最小流量;(2)長江干流宜昌站流量為5280 m3/s,代表三峽水庫運(yùn)行前1月中旬最大出庫流量以及三峽水庫運(yùn)行后1月中旬最小出庫流量,對應(yīng)的清江及四水流量取1月中旬同頻率流量;(3)三峽水庫出庫流量為9500 m3/s,代表三峽水庫運(yùn)行后1月中旬最大出庫流量,對應(yīng)的清江及四水流量取1月中旬最大流量.

    3.1 三峽水庫不同出庫流量下洞庭湖水位空間分布格局

    三峽水庫不同典型工況出庫流量下洞庭湖水位的空間分布如圖7所示,當(dāng)水庫出庫流量為3740 m3/s時,城陵磯、鹿角、楊柳潭、南咀的水位分別為20.31、21.63、27.43、31.58 m.當(dāng)水庫出庫流量為5280 m3/s時,城陵磯、鹿角、楊柳潭、南咀的水位分別為20.64、22.26、27.64、31.28 m.當(dāng)出庫流量為9500 m3/s時,城陵磯、鹿角、楊柳潭、南咀的水位分別為21.64、22.46、27.94、31.67 m.從水位增幅的空間分布來看,三峽水庫運(yùn)行后增加下泄流量對洞庭湖水位的提高主要集中在東洞庭湖北部,城陵磯水位最大提高1.33 m,而鹿角、楊柳潭、南咀水位最大分別提高0.83、0.51、0.39m.由此可見,枯水期1月中旬三峽水庫增加出庫流量時,洞庭湖水位比水庫運(yùn)行前有所升高,但受洞庭湖西高東低地形的影響,水庫增加出庫流量對東洞庭湖北部的水位提升起到明顯的作用,而對于南洞庭湖楊柳潭站和西洞庭湖的南咀站作用較小,因此可以認(rèn)為三峽水庫枯水期增加下泄流量對于維持東洞庭湖水位作用明顯,但對南洞庭湖和西洞庭湖而言作用有限.

    圖7 三峽水庫不同典型出庫流量下洞庭湖水位空間分布

    3.2 洞庭湖水深變化對白鶴攝食棲息地適宜度指數(shù)的影響

    圖8顯示了三峽水庫不同典型出庫流量下洞庭湖白鶴攝食棲息地適宜度指數(shù)空間分布特征,其中顏色越接近于紅色表示適宜度越大,越接近藍(lán)色表示適宜度越小.三峽水庫運(yùn)行前長江宜昌站流量為 3740 m3/s,洞庭湖白鶴攝食棲息地平均水深為1.03 m,不同區(qū)域的平均適宜度指數(shù)為0.034;當(dāng)水庫出庫流量為5280 m3/s時,攝食棲息地平均水深為1.21 m,不同區(qū)域的適宜度指數(shù)較小,平均值為0.040;當(dāng)出庫流量為9500 m3/s時,棲息地平均水深為1.50 m,水深適宜度指數(shù)平均值達(dá)到0.043.從空間分布來看,適宜度指數(shù)較高的區(qū)域主要分布在南洞庭湖南部和西洞庭湖南部,此區(qū)域地勢平坦,形成了大面積的淺水區(qū)域,而東洞庭湖在枯水期呈“河相”,水流只在河槽中流動,水深增加對棲息地適宜度指數(shù)影響較小.

    圖8 三峽水庫不同典型出庫流量下洞庭湖適宜度指數(shù)空間分布

    3.3 三峽水庫不同出庫流量對白鶴攝食棲息地面積的影響

    適宜度指數(shù)直接影響白鶴攝食棲息地加權(quán)可利用面積的大小,基于構(gòu)建的PHABSIM模型得到不同出庫流量對應(yīng)的洞庭湖白鶴攝食棲息地加權(quán)可利用面積如圖9所示.由圖可知三峽水庫運(yùn)行后出庫流量由小變大的過程中,洞庭湖白鶴攝食棲息地的加權(quán)可利用面積總體呈現(xiàn)緩慢增大的趨勢.三峽水庫出庫流量為3740 m3/s時,棲息地加權(quán)可利用面積101.40 km2;三峽水庫出庫流量為5280 m3/s時,棲息地加權(quán)可利用面積113.44 km2;當(dāng)出庫流量在6000~9500 m3/s之間時,加權(quán)可利用面積逐步增大;出庫流量為9500 m3/s,加權(quán)可利用面積為121.84 km2.由此可見三峽水庫運(yùn)行后洞庭湖白鶴攝食棲息地面積最大增加20.44 km2,對應(yīng)的增幅為20.16%,說明水庫運(yùn)行后洞庭湖枯水期白鶴潛在的攝食棲息地適宜生境面積隨著三峽水庫出庫流量的增加呈增大的趨勢,從而印證了三峽水庫補(bǔ)水調(diào)度的積極作用.

    圖9 三峽水庫出庫流量與棲息地加權(quán)可利用面積的關(guān)系

    洞庭湖不同湖區(qū)攝食棲息地加權(quán)可利用面積對三峽水庫出庫流量的響應(yīng)關(guān)系如圖10所示,不同湖區(qū)加權(quán)可利用面積總體保持相對穩(wěn)定并呈逐步增大的趨勢.三峽水庫蓄水前當(dāng)長江宜昌站流量為3740 m3/s時,東洞庭湖攝食棲息地加權(quán)可利用面積僅為7.49 km2,當(dāng)三峽水庫出庫流量為5280 m3/s,攝食棲息地加權(quán)可利用面積持續(xù)增大,達(dá)到9.28 km2,隨著出庫流量的繼續(xù)增大,加權(quán)可利用面積逐步增大,當(dāng)流量為9500 m3/s 時,加權(quán)可利用面積達(dá)到9.86 km2.對于南洞庭湖(含橫嶺湖),加權(quán)可利用面積變化趨勢與東洞庭湖一致,當(dāng)流量為3740 m3/s時,南洞庭湖白鶴攝食棲息地加權(quán)可利用面積為47.37 km2,當(dāng)出庫流量為9500 m3/s時,棲息地加權(quán)可利用面積達(dá)到最大,為60.34 km2.西洞庭主要由沅江、澧水及三口河系的松滋河匯聚而成,隨著三峽水庫出庫流量的增加,松滋口分流量增加幅度較小,根據(jù)水動力模型計算結(jié)果:出庫流量由5280 m3/s增加至9500 m3/s時,松滋河分流量增加20.8 m3/s,達(dá)到63.2 m3/s,白鶴攝食所需的加權(quán)可利用面積呈緩慢增大的趨勢.由圖 10也可看出,洞庭湖白鶴攝食棲息地加權(quán)可利用面積主要分布在南洞庭湖和西洞庭湖,兩者之和占總面積的90%以上,從白鶴數(shù)量分布來看,主要分布在南洞庭湖和西洞庭湖濕地,也正說明廣袤的淺水區(qū)域為白鶴攝食提供了良好的條件.

    圖10 洞庭湖不同湖區(qū)白鶴攝食棲息地面積變化趨勢

    4 結(jié)論

    本文以洞庭湖珍稀越冬水鳥——白鶴作為指示性候鳥,以水深作為白鶴攝食的關(guān)鍵生境因子,構(gòu)建針對洞庭湖枯水期越冬白鶴攝食的物理棲息地模型,對三峽水庫運(yùn)行前后不同出庫流量條件下洞庭湖枯水期白鶴潛在生境面積變化進(jìn)行模擬.

    1)三峽水庫運(yùn)行后1月中旬出庫流量范圍在5280~9500 m3/s,洞庭湖城陵磯、鹿角、楊柳潭、南咀水位最大提高分別為1.33、0.83、0.51、0.39 m,說明三峽水庫增加出庫流量對洞庭湖水位的影響主要集中在東洞庭湖北部,對其它湖區(qū)影響相對較小.

    2)三峽水庫出庫流量由小變大的過程中,白鶴攝食棲息地的加權(quán)可利用面積總體呈現(xiàn)逐步增大的趨勢,棲息地面積維持在101.40~121.84 km2之間,最大值與最小值差值為20.44 km2,對應(yīng)的增幅為20.16%.

    3)不同湖區(qū)加權(quán)可利用面積隨著三峽出庫流量的增加均呈逐漸增大的趨勢,洞庭湖白鶴攝食棲息地主要分布在南洞庭湖和西洞庭湖,兩者之和占總面積的90%以上.

    維持洞庭湖枯水期適宜水位的措施除了三峽水庫補(bǔ)水調(diào)度外,還應(yīng)在楊樹清理工作完成后持續(xù)開展?jié)竦厣鷳B(tài)恢復(fù).此外,白鶴在數(shù)量還受歷史棲息地選擇、忠誠度、食物資源、人類干擾、湖床底質(zhì)等多種因素共同影響,水深僅是白鶴攝食棲息地的限制性生境因子之一,后續(xù)的研究將沿著“水文、氣候等脅迫-棲息地生境(食物、攝食可及性)-水鳥數(shù)量及分布”這一思路,構(gòu)建多種關(guān)鍵生境因子共同作用下的棲息地模型,詳細(xì)全面評估棲息地質(zhì)量的演變特征,可為通過三峽水庫生態(tài)補(bǔ)水調(diào)度改善白鶴等更多洞庭湖越冬水鳥的棲息地生境提供重要的科學(xué)依據(jù).

    猜你喜歡
    三峽水庫白鶴洞庭湖
    白鶴之約
    心聲歌刊(2023年2期)2023-07-16 05:56:02
    獨(dú)孤求敗白鶴灘
    哈哈畫報(2022年1期)2022-04-19 11:27:20
    洞庭湖
    白鶴的真愛
    黃河之聲(2021年7期)2021-06-21 13:59:24
    三峽水庫迎戰(zhàn)長江2020年第2號洪水
    中國水利(2020年14期)2020-08-02 16:19:30
    輕松松聊漢語 洞庭湖
    金橋(2020年8期)2020-05-22 06:22:54
    好一個洞庭湖
    洞庭湖的麋鹿
    三峽水庫澎溪河水-氣界面CO2、CH4擴(kuò)散通量晝夜動態(tài)初探
    “白鶴”飛得高靠的抓“四好”
    中國火炬(2013年3期)2013-07-24 14:15:07
    亚洲精品在线观看二区| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 色av中文字幕| 99国产综合亚洲精品| 黄色女人牲交| 一进一出好大好爽视频| 精品人妻1区二区| 国产午夜福利久久久久久| 丝袜在线中文字幕| 婷婷精品国产亚洲av在线| 久久久国产欧美日韩av| 国产精品久久久av美女十八| av欧美777| 身体一侧抽搐| 国产亚洲精品一区二区www| xxx96com| 人妻久久中文字幕网| 国产亚洲欧美精品永久| 欧美另类亚洲清纯唯美| 最近最新中文字幕大全电影3 | 熟女少妇亚洲综合色aaa.| 精品国产一区二区久久| tocl精华| 久久性视频一级片| 啪啪无遮挡十八禁网站| 日韩精品青青久久久久久| 亚洲欧美激情在线| www日本在线高清视频| 天天添夜夜摸| 欧美老熟妇乱子伦牲交| 国产真人三级小视频在线观看| 天天躁夜夜躁狠狠躁躁| 久久久国产成人精品二区| 国产精品美女特级片免费视频播放器 | 亚洲成人精品中文字幕电影| 国产成人系列免费观看| 国产高清视频在线播放一区| 久久久久国产精品人妻aⅴ院| 啦啦啦观看免费观看视频高清 | 12—13女人毛片做爰片一| 妹子高潮喷水视频| 久久久久久人人人人人| 乱人伦中国视频| av天堂在线播放| 日韩精品免费视频一区二区三区| 精品乱码久久久久久99久播| 久久久国产精品麻豆| 欧美激情高清一区二区三区| 精品欧美国产一区二区三| 婷婷丁香在线五月| 欧美一级a爱片免费观看看 | 9热在线视频观看99| 精品国产乱子伦一区二区三区| 中文字幕久久专区| 国产成人免费无遮挡视频| 精品高清国产在线一区| 欧美日本视频| 夜夜看夜夜爽夜夜摸| 在线观看一区二区三区| 国产成人精品无人区| 欧美丝袜亚洲另类 | 国产精品自产拍在线观看55亚洲| 国产国语露脸激情在线看| 日韩欧美在线二视频| 精品一区二区三区av网在线观看| svipshipincom国产片| 波多野结衣一区麻豆| 嫁个100分男人电影在线观看| 两性夫妻黄色片| 国产欧美日韩一区二区精品| 人人妻人人澡人人看| 正在播放国产对白刺激| 成人av一区二区三区在线看| videosex国产| 午夜老司机福利片| 国产高清videossex| 国产激情欧美一区二区| 50天的宝宝边吃奶边哭怎么回事| 午夜福利成人在线免费观看| 夜夜夜夜夜久久久久| 日本vs欧美在线观看视频| 国产成人精品无人区| 视频在线观看一区二区三区| 精品欧美国产一区二区三| 色av中文字幕| 欧美日韩精品网址| 中文字幕人妻丝袜一区二区| 精品无人区乱码1区二区| 91九色精品人成在线观看| 欧美+亚洲+日韩+国产| 亚洲男人天堂网一区| 国产精品精品国产色婷婷| 欧美日韩乱码在线| АⅤ资源中文在线天堂| 在线永久观看黄色视频| 亚洲成av片中文字幕在线观看| 成人免费观看视频高清| 日韩欧美免费精品| 午夜福利视频1000在线观看 | 1024香蕉在线观看| cao死你这个sao货| 在线观看日韩欧美| 亚洲成国产人片在线观看| 精品国产一区二区久久| 久久久精品国产亚洲av高清涩受| 人人妻人人澡人人看| 日本三级黄在线观看| 精品午夜福利视频在线观看一区| 国产人伦9x9x在线观看| 天天一区二区日本电影三级 | 成人国语在线视频| 丝袜美足系列| 黄网站色视频无遮挡免费观看| 亚洲专区中文字幕在线| 麻豆久久精品国产亚洲av| 日韩免费av在线播放| 最好的美女福利视频网| 国产又爽黄色视频| 又大又爽又粗| 国产精华一区二区三区| 久久久国产成人免费| 性少妇av在线| 无限看片的www在线观看| 欧美在线黄色| 国产成人精品在线电影| 一级,二级,三级黄色视频| 久久中文字幕人妻熟女| 国产精品 国内视频| 亚洲中文av在线| 亚洲成av片中文字幕在线观看| www日本在线高清视频| 91九色精品人成在线观看| 国产99久久九九免费精品| 精品久久久久久久久久免费视频| √禁漫天堂资源中文www| 久久人妻av系列| 久久人人爽av亚洲精品天堂| 窝窝影院91人妻| 免费在线观看亚洲国产| 嫩草影院精品99| avwww免费| 亚洲第一av免费看| 日本 av在线| 久久精品影院6| 女性生殖器流出的白浆| 男女之事视频高清在线观看| 日韩av在线大香蕉| 久久亚洲真实| 制服诱惑二区| 亚洲免费av在线视频| 两个人看的免费小视频| netflix在线观看网站| 在线十欧美十亚洲十日本专区| 亚洲欧美精品综合一区二区三区| 成年版毛片免费区| 欧美日韩黄片免| 99国产极品粉嫩在线观看| 久久精品aⅴ一区二区三区四区| www.999成人在线观看| 亚洲性夜色夜夜综合| 久久亚洲真实| 又大又爽又粗| 可以免费在线观看a视频的电影网站| 怎么达到女性高潮| 又紧又爽又黄一区二区| 色哟哟哟哟哟哟| 99国产极品粉嫩在线观看| 少妇被粗大的猛进出69影院| 亚洲激情在线av| 午夜福利免费观看在线| 国产成人精品无人区| a在线观看视频网站| 精品欧美国产一区二区三| 精品日产1卡2卡| 精品久久久久久成人av| 自拍欧美九色日韩亚洲蝌蚪91| xxx96com| 女人爽到高潮嗷嗷叫在线视频| 777久久人妻少妇嫩草av网站| 精品国产超薄肉色丝袜足j| 亚洲精品在线观看二区| 久久中文看片网| 成人免费观看视频高清| 69av精品久久久久久| 精品欧美国产一区二区三| 国产精品秋霞免费鲁丝片| e午夜精品久久久久久久| 老司机午夜福利在线观看视频| 亚洲欧洲精品一区二区精品久久久| 九色国产91popny在线| 在线观看舔阴道视频| 中文亚洲av片在线观看爽| 亚洲欧美日韩无卡精品| 亚洲中文av在线| 嫁个100分男人电影在线观看| 精品国产一区二区久久| 国产精品 欧美亚洲| 美国免费a级毛片| 少妇粗大呻吟视频| 欧美黑人欧美精品刺激| 日韩视频一区二区在线观看| 中国美女看黄片| 亚洲成a人片在线一区二区| 久久热在线av| 欧美绝顶高潮抽搐喷水| 精品国产超薄肉色丝袜足j| 97人妻精品一区二区三区麻豆 | 精品一区二区三区av网在线观看| 超碰成人久久| 此物有八面人人有两片| 国产精品久久久久久精品电影 | 国产一区二区三区综合在线观看| 老鸭窝网址在线观看| 非洲黑人性xxxx精品又粗又长| 国产日韩一区二区三区精品不卡| 亚洲一区中文字幕在线| 午夜福利高清视频| 男人舔女人的私密视频| 日韩欧美一区视频在线观看| 久久久国产精品麻豆| 一本久久中文字幕| 少妇的丰满在线观看| 黄色a级毛片大全视频| 久久精品亚洲熟妇少妇任你| 每晚都被弄得嗷嗷叫到高潮| 每晚都被弄得嗷嗷叫到高潮| 村上凉子中文字幕在线| 亚洲av成人av| 亚洲国产欧美一区二区综合| 很黄的视频免费| 久久国产亚洲av麻豆专区| 黄色 视频免费看| 一二三四社区在线视频社区8| 人人妻人人澡人人看| 91麻豆精品激情在线观看国产| 淫秽高清视频在线观看| 三级毛片av免费| 亚洲中文日韩欧美视频| 一进一出抽搐gif免费好疼| 天堂动漫精品| 在线观看日韩欧美| 黄色丝袜av网址大全| 国产极品粉嫩免费观看在线| 国产人伦9x9x在线观看| 777久久人妻少妇嫩草av网站| 午夜福利18| 美女高潮到喷水免费观看| 这个男人来自地球电影免费观看| 久久草成人影院| 亚洲午夜理论影院| 国产视频一区二区在线看| 亚洲精品久久成人aⅴ小说| 午夜亚洲福利在线播放| 国产成人欧美在线观看| 日日夜夜操网爽| 精品国产超薄肉色丝袜足j| 免费av毛片视频| 日日干狠狠操夜夜爽| 十分钟在线观看高清视频www| 真人做人爱边吃奶动态| 少妇粗大呻吟视频| 天堂动漫精品| 亚洲人成电影观看| 国产精品一区二区在线不卡| 国产精品1区2区在线观看.| 亚洲精品久久成人aⅴ小说| а√天堂www在线а√下载| 很黄的视频免费| 在线国产一区二区在线| 国产欧美日韩一区二区三区在线| 女生性感内裤真人,穿戴方法视频| 大香蕉久久成人网| 欧美国产日韩亚洲一区| 一边摸一边抽搐一进一出视频| av在线播放免费不卡| 国产成人精品久久二区二区91| 成人三级做爰电影| 欧美国产精品va在线观看不卡| 国产精品久久久久久亚洲av鲁大| 亚洲无线在线观看| 可以在线观看的亚洲视频| 精品高清国产在线一区| 成人特级黄色片久久久久久久| 国产91精品成人一区二区三区| 视频在线观看一区二区三区| 日本在线视频免费播放| 侵犯人妻中文字幕一二三四区| 亚洲精品久久成人aⅴ小说| 最近最新免费中文字幕在线| 欧美另类亚洲清纯唯美| 岛国在线观看网站| 丁香欧美五月| 欧美成人性av电影在线观看| 老司机福利观看| 欧美不卡视频在线免费观看 | 在线观看免费午夜福利视频| 夜夜躁狠狠躁天天躁| 老司机在亚洲福利影院| 国产成人系列免费观看| 久久人妻熟女aⅴ| 精品久久久久久成人av| 成人国产综合亚洲| 国产一区在线观看成人免费| 欧美成人性av电影在线观看| www.精华液| 精品国产亚洲在线| 88av欧美| 麻豆国产av国片精品| 可以在线观看毛片的网站| 久久久久久人人人人人| 国产亚洲欧美在线一区二区| 亚洲人成伊人成综合网2020| 久久精品亚洲精品国产色婷小说| 丁香欧美五月| 一边摸一边抽搐一进一小说| 国产精品 欧美亚洲| 亚洲三区欧美一区| 黄片小视频在线播放| 大香蕉久久成人网| 亚洲中文日韩欧美视频| 亚洲国产精品sss在线观看| 国产成人精品无人区| 久久久久国产一级毛片高清牌| 成人永久免费在线观看视频| 久久久国产欧美日韩av| 一个人观看的视频www高清免费观看 | 999久久久精品免费观看国产| 变态另类丝袜制服| 18禁黄网站禁片午夜丰满| 9191精品国产免费久久| 久久久国产欧美日韩av| 18禁美女被吸乳视频| 国产精品久久电影中文字幕| 亚洲欧美精品综合一区二区三区| 88av欧美| 精品无人区乱码1区二区| 欧美绝顶高潮抽搐喷水| 男人舔女人下体高潮全视频| 国产色视频综合| 亚洲av成人一区二区三| 搡老岳熟女国产| 国产亚洲av嫩草精品影院| 午夜日韩欧美国产| 中文字幕人成人乱码亚洲影| 久久精品91无色码中文字幕| 淫妇啪啪啪对白视频| 乱人伦中国视频| 黄频高清免费视频| xxx96com| 久久人人97超碰香蕉20202| 90打野战视频偷拍视频| 91在线观看av| 国产又爽黄色视频| 国产成人影院久久av| 精品一区二区三区av网在线观看| 大型黄色视频在线免费观看| 久久久久久久精品吃奶| av天堂久久9| 神马国产精品三级电影在线观看 | 国产视频一区二区在线看| 久久香蕉国产精品| 大码成人一级视频| 亚洲五月色婷婷综合| 啦啦啦 在线观看视频| 18禁黄网站禁片午夜丰满| 国产亚洲精品一区二区www| 久热这里只有精品99| 国产欧美日韩综合在线一区二区| 免费在线观看视频国产中文字幕亚洲| 欧美黑人精品巨大| 一区二区三区精品91| 欧美一级a爱片免费观看看 | 国产亚洲欧美精品永久| 亚洲精品中文字幕在线视频| 91成年电影在线观看| 亚洲av电影不卡..在线观看| 日韩中文字幕欧美一区二区| 露出奶头的视频| 国产av一区在线观看免费| 午夜福利高清视频| 久久精品国产亚洲av高清一级| 国产精品一区二区免费欧美| 最好的美女福利视频网| 国产亚洲精品av在线| 精品人妻在线不人妻| 非洲黑人性xxxx精品又粗又长| 亚洲最大成人中文| 中文字幕人妻丝袜一区二区| 一区在线观看完整版| 欧美性长视频在线观看| 深夜精品福利| 国产亚洲欧美98| 精品少妇一区二区三区视频日本电影| or卡值多少钱| 久久精品国产99精品国产亚洲性色 | 无限看片的www在线观看| 天堂动漫精品| 午夜福利成人在线免费观看| 18禁观看日本| 久久久久久大精品| 正在播放国产对白刺激| 乱人伦中国视频| 又大又爽又粗| 久久 成人 亚洲| 免费av毛片视频| 老司机午夜十八禁免费视频| 麻豆成人av在线观看| 久久久国产精品麻豆| 久久热在线av| 亚洲精品一区av在线观看| 国产精品 国内视频| 99国产精品一区二区三区| 欧美日本中文国产一区发布| 熟女少妇亚洲综合色aaa.| 日韩精品青青久久久久久| 亚洲aⅴ乱码一区二区在线播放 | 久久中文字幕人妻熟女| 黑人欧美特级aaaaaa片| 女人被躁到高潮嗷嗷叫费观| 亚洲精华国产精华精| 9191精品国产免费久久| 男女做爰动态图高潮gif福利片 | 午夜福利影视在线免费观看| 国产91精品成人一区二区三区| 女人精品久久久久毛片| 一级黄色大片毛片| 国产精品九九99| 国产欧美日韩一区二区精品| 午夜福利影视在线免费观看| 男人舔女人下体高潮全视频| 欧美成人性av电影在线观看| 色播在线永久视频| 国产精品,欧美在线| 女人被狂操c到高潮| 国产一级毛片七仙女欲春2 | 欧美在线一区亚洲| 一级毛片高清免费大全| 国产精品免费一区二区三区在线| 国内精品久久久久久久电影| 成人特级黄色片久久久久久久| 国产伦人伦偷精品视频| 成人欧美大片| 国产精品香港三级国产av潘金莲| 青草久久国产| 欧美国产精品va在线观看不卡| 久久香蕉国产精品| 级片在线观看| 神马国产精品三级电影在线观看 | 亚洲一码二码三码区别大吗| 亚洲一码二码三码区别大吗| 黑人操中国人逼视频| 久久天堂一区二区三区四区| 欧美成人免费av一区二区三区| 精品电影一区二区在线| 一边摸一边抽搐一进一小说| 女性被躁到高潮视频| 成在线人永久免费视频| 国产激情欧美一区二区| 一进一出好大好爽视频| 女人被狂操c到高潮| 悠悠久久av| 一边摸一边抽搐一进一小说| 亚洲视频免费观看视频| 日韩一卡2卡3卡4卡2021年| 亚洲成人久久性| 午夜免费激情av| 老司机在亚洲福利影院| 免费女性裸体啪啪无遮挡网站| 91老司机精品| 久久天堂一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 女人高潮潮喷娇喘18禁视频| 自拍欧美九色日韩亚洲蝌蚪91| 精品电影一区二区在线| 在线观看免费日韩欧美大片| 手机成人av网站| 少妇裸体淫交视频免费看高清 | 在线观看免费视频网站a站| 热99re8久久精品国产| 久久精品91蜜桃| 69av精品久久久久久| 中文字幕高清在线视频| 99国产精品免费福利视频| 国产精品亚洲av一区麻豆| 亚洲伊人色综图| 日本一区二区免费在线视频| 久久久久久免费高清国产稀缺| 男女下面插进去视频免费观看| 免费在线观看视频国产中文字幕亚洲| 黄色视频不卡| 最近最新中文字幕大全电影3 | 日韩中文字幕欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲美女久久久| 国产成人精品在线电影| 免费观看精品视频网站| 日韩欧美三级三区| 国产精品一区二区在线不卡| 久久亚洲精品不卡| 成年女人毛片免费观看观看9| 亚洲少妇的诱惑av| 亚洲色图 男人天堂 中文字幕| 亚洲国产高清在线一区二区三 | 国产精品二区激情视频| 久久人妻av系列| 宅男免费午夜| 天堂影院成人在线观看| 亚洲精品久久成人aⅴ小说| 国产激情久久老熟女| 亚洲国产精品999在线| 18禁黄网站禁片午夜丰满| 12—13女人毛片做爰片一| 欧美成狂野欧美在线观看| 不卡一级毛片| 嫁个100分男人电影在线观看| 久热爱精品视频在线9| 亚洲精品中文字幕一二三四区| 国产伦一二天堂av在线观看| 大陆偷拍与自拍| 欧美乱码精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 亚洲精品av麻豆狂野| 琪琪午夜伦伦电影理论片6080| 午夜免费鲁丝| 国产精品一区二区三区四区久久 | 欧美成狂野欧美在线观看| 久久久久久免费高清国产稀缺| 一区二区三区精品91| 999精品在线视频| 午夜精品久久久久久毛片777| 身体一侧抽搐| tocl精华| 亚洲精品中文字幕在线视频| 国产av又大| 丝袜在线中文字幕| 69精品国产乱码久久久| 国产一区二区三区视频了| 后天国语完整版免费观看| 欧美日韩黄片免| 制服丝袜大香蕉在线| 又黄又粗又硬又大视频| 色综合站精品国产| 中文字幕高清在线视频| videosex国产| 亚洲国产欧美网| or卡值多少钱| 一级毛片高清免费大全| 90打野战视频偷拍视频| 精品卡一卡二卡四卡免费| 亚洲av日韩精品久久久久久密| 一级毛片高清免费大全| 亚洲性夜色夜夜综合| 午夜福利影视在线免费观看| 欧美成人性av电影在线观看| 欧美精品啪啪一区二区三区| 国产精品 国内视频| 久久香蕉精品热| www.www免费av| 亚洲电影在线观看av| 久久九九热精品免费| 午夜视频精品福利| 久久这里只有精品19| 久久国产乱子伦精品免费另类| 日韩一卡2卡3卡4卡2021年| 此物有八面人人有两片| 国内精品久久久久久久电影| 亚洲自拍偷在线| 久久久久国产一级毛片高清牌| 美女午夜性视频免费| 日本在线视频免费播放| 欧美av亚洲av综合av国产av| 99精品久久久久人妻精品| 12—13女人毛片做爰片一| 男人操女人黄网站| 国产精品99久久99久久久不卡| 免费高清在线观看日韩| 一区二区三区国产精品乱码| 亚洲九九香蕉| 午夜福利成人在线免费观看| 久久狼人影院| 亚洲精品国产一区二区精华液| 精品国产美女av久久久久小说| 高潮久久久久久久久久久不卡| 19禁男女啪啪无遮挡网站| 亚洲av片天天在线观看| 香蕉丝袜av| 亚洲一区中文字幕在线| 母亲3免费完整高清在线观看| 亚洲成人免费电影在线观看| 午夜精品久久久久久毛片777| 亚洲伊人色综图| 婷婷丁香在线五月| 91麻豆精品激情在线观看国产| 我的亚洲天堂| av有码第一页| 香蕉久久夜色| av视频在线观看入口| 又紧又爽又黄一区二区| 国产色视频综合| 麻豆久久精品国产亚洲av| 免费在线观看亚洲国产| 丁香欧美五月| 亚洲国产精品久久男人天堂| 十八禁网站免费在线| 9191精品国产免费久久| 色播亚洲综合网| 久久久久久免费高清国产稀缺| 欧美日韩福利视频一区二区| 久久久久亚洲av毛片大全| 久久草成人影院| 18禁国产床啪视频网站| 久久 成人 亚洲|