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

    未來氣候變化下嫩江流域生態(tài)流量與生物多樣性演變預估

    2025-03-22 00:00:00馬秋生陳末孫佳豪章清松吳燕鋒
    水生態(tài)學雜志 2025年2期

    摘要:分析嫩江流域未來氣候變化下的生態(tài)流量變化和生物多樣性演變,為流域水資源管理和水生態(tài)系統(tǒng)保護提供科學依據(jù)。以嫩江流域為研究對象,以大賚水文站近55年的日徑流量數(shù)據(jù)和基于CMIP6氣候情景數(shù)據(jù),利用率定后的 HYDROTEL模型,在模擬未來日徑流量的基礎上,采用水文變異指標法評估未來生態(tài)流量、環(huán)境流量和生物多樣性的變化。研究表明,與歷史時期比較,在未來SSP126、SSP370和SSP585這3種氣候情景下,嫩江流域水文指標將發(fā)生較大改變,來水不確定性增大,生態(tài)流量穩(wěn)定性降低, SSP126和SSP370情境下總改變程度分別為69%和75%。環(huán)境流枯水流量事件和高脈沖流量事件頻率在不同情景下均有增加,SSP126情景下增加3.23%和0.54%,SSP370情景下增加7.56%和0.38%,SSP585情景下增加8.01%和1.02%。嫩江流域環(huán)境流組成中大部分流量過程歸入枯水流量和高流量脈沖事件,流量變化范圍縮小,組成趨于單一。與歷史時期比較(生物多樣性指數(shù)均值為0.29),未來氣候變化下,嫩江流域生物多樣性呈增加趨勢,在SPP126、SSP137和SSP585情景下生物多樣性指數(shù)均值分別為0.39、0.36和0.37。

    關鍵詞:生態(tài)流量;環(huán)境流量;氣候變化;生物多樣性;嫩江流域

    中圖分類號:TV211.1;X143" " " " 文獻標志碼:A" " " " 文章編號:1674-3075(2026)02-0082-09

    生態(tài)流量研究始于1940年,最初主要集中在探索魚類生長繁殖與河流流量的關系,以緩解水利工程建設和水資源開發(fā)利用對漁業(yè)生產(chǎn)的影響(Acreman et al,2004)。隨后,河流生態(tài)需水理論開始萌芽,學者們陸續(xù)提出了河道內(nèi)流量、最小流量以及環(huán)境流量等概念,生態(tài)流量研究受到了管理部門和學術領域的普遍關注。近年來,流域尺度的生態(tài)流量逐漸引起學者們的關注,河流生態(tài)流量保障逐漸成為流域水資源管理的重要內(nèi)容之一(Poff et al,2013;邴建平等,2023;張輝等,2022)。流域生態(tài)流量是指確保流域生態(tài)系統(tǒng)穩(wěn)定健康且滿足人類發(fā)展需要的流量,其更關注河流流量與生態(tài)系統(tǒng)保護、維持和完整性恢復的關系(田震等,2023)。經(jīng)過不斷探索和研究,河流生態(tài)流量的評價方法也不斷發(fā)展,如變異范圍法、水文時間變異檢驗法和Tennant法等。有效評估流域水文情勢及河流和環(huán)境流量的變化,對維護流域生態(tài)系統(tǒng)穩(wěn)定健康有著重要意義。

    河流生態(tài)流量涵蓋水文情勢和環(huán)境流量2個方面。在氣候變化和高強度的人類活動影響下,全球水文系統(tǒng)發(fā)生了明顯改變,導致河流的水文情勢、環(huán)境流量組分和生物多樣性的改變,流域生態(tài)系統(tǒng)完整性受到的影響日益加?。ㄍ醮箅p等,2023)。由于生物多樣性與河流生態(tài)系統(tǒng)的多功能性密切相關(王凱等,2022;Sun et al,2017),開展氣候變化下流域生態(tài)流量和生物多樣性影響評估,成為保護流域生態(tài)系統(tǒng)的新需求,也越來越受到國內(nèi)外學者的關注。劉綠柳等(2021)研究了氣候變化下黃河流域生態(tài)流量的變化,發(fā)現(xiàn)在未來氣候下,人類活動將會對黃河流域產(chǎn)生徑流量減少等負面影響。李澤君等(2021)研究了變化環(huán)境下韓江生態(tài)流量演變特征,發(fā)現(xiàn)未來氣候下生態(tài)流量對氣候變化的響應將更加敏感。綜合國內(nèi)外研究,目前缺乏對未來多種氣候變化情景下流域水文演變和生態(tài)效應的探討。已有研究多集中在單一氣候條件下的模擬分析,對未來氣候變化下生態(tài)流量、環(huán)境流量和生物多樣性變化的研究較為有限。此外,更加復雜多變的氣候條件會給流域生態(tài)流量、環(huán)境流量變化及生物多樣帶來更多不確定性,因而迫切需要深入探究多氣候變化情景下生態(tài)流量和生物多樣性的演變規(guī)律。

    為此,本文選擇嫩江流域為研究區(qū),利用嫩江流域出口控制斷面水文站——大賚水文站近55年的日徑流量數(shù)據(jù)和CMIP6氣候情景數(shù)據(jù),采用HYDROTEL模型模擬未來日徑流量數(shù)據(jù),在計算水文改變程度指標、環(huán)境流指標和生物多樣性指標的基礎上,綜合分析嫩江流域未來氣候變化下的生態(tài)流量變化和生物多樣性演變,以期為未來氣候下嫩江流域的水資源管理和水生態(tài)系統(tǒng)保護提供科學依據(jù)。

    1" "材料與方法

    1.1" "研究區(qū)概況

    嫩江是松花江的北源,發(fā)源于內(nèi)蒙古自治區(qū)境內(nèi)大興安嶺-伊勒呼里山的中段南側(cè)。嫩江干流流經(jīng)黑龍江省與內(nèi)蒙古自治區(qū)、吉林省的交界地帶,最后在肇源縣茂興鎮(zhèn)當權村三岔河與松花江南源匯合為松花江。嫩江全長1 370 km,流域面積29.7萬km2,多年平均流量823.4 m3/s,年均徑流總量200億m3。嫩江流域徑流量的年內(nèi)和年際變化較大,每年的1、2、3、11、12月是枯水和枯冰期;4、5、6、10月是平水期;7—9月份是豐水期。此外,嫩江流域也是多種水生生物和野生動物的棲息地,魚類資源十分豐富。嫩江流域區(qū)域位置如圖1。

    1.2" "研究數(shù)據(jù)

    選取嫩江流域大賚水文站1963—2018年的日徑流量作為歷史時期數(shù)據(jù),以第6次耦合模式比較計劃(CMIP6)氣候情景數(shù)據(jù)為驅(qū)動(趙宗慈等,2016),采用HYDROTEL模型模擬未來氣候變化下2025—2100年的日徑流量。

    選取CMIP6中的不同共享社會經(jīng)濟路徑(shared socioeconomic pathways, SSP)(張麗霞等,2019)中的SSP126、SSP370和SSP585情景數(shù)據(jù),以此驅(qū)動HYDROTEL模型(周天軍等,2019),研究3種氣候變化情景下嫩江流域的生態(tài)流量和生物多樣性的變化。3種情景模式詳細信息見表1。

    1.3" "研究方法

    基于歷時和模擬的日徑流量數(shù)據(jù),計算水文情勢改變指標、環(huán)境流指標和生物多樣性指數(shù),對比分析歷史時期和不同SSPs情景下模擬的3個指標的差異性。從水文情勢改變指標和環(huán)境流指標分析生態(tài)流量的變化,從生物多樣性指標視角分析生物多樣性的演變(周毅等,2017a)。

    1.3.1" "HYDROTEL" "水文模型 利用率定和驗證后的HYDROTEL模型在嫩江流域開展未來氣候變化下日徑流模擬。HYDROTEL是一個分布式水文模型,主要由氣象數(shù)據(jù)插值模塊、積雪模塊、凍土模塊、蒸散發(fā)模塊、垂向水量收支平衡模塊、陸面產(chǎn)匯流模塊和河道徑流模塊組成。吳燕鋒等(2019)在嫩江流域構建了耦合濕地模塊的HYDROTEL水文模型,在率定和驗證期間的納什系數(shù)分別達到0.74和0.72,克林效率系數(shù)分別達到0.79和0.78,顯示出該模型在嫩江流域有較好的適用性。本研究基于吳燕鋒等(2019)已構建的模型,以CMIP6的3種排放情景數(shù)據(jù)為驅(qū)動,開展未來氣候變化下日徑流模擬。

    1.3.2" "水文情勢改變指標" "水文情勢指標可以表征流域流量基本特征與生態(tài)環(huán)境的關系。Richer等(1998)總結(jié)了33個水文改變指標(indicators of hydrologic alteration,IHA),用于反映河流流態(tài)年內(nèi)及年際特征變化,在河流生態(tài)流量管理和天然水文情勢保障方面應用廣泛。IHA 涵蓋了 5 組共 33 個水文指標(Zhang et al,2014;Chen et al,2010),這些指標體現(xiàn)了徑流量的量級、出現(xiàn)時間、持續(xù)時間、頻率等水文特征。由于嫩江流域并未出現(xiàn)日徑流量斷流的現(xiàn)象,因此無需考慮“斷流天數(shù)”這一水文指標。其余32個水文指標如表2。

    采用32個IHA指標反映河流的水文情勢狀況,并從流量、時間、頻率、延時以及變化率5個方面來反映河道水文情勢。為了量化水文指標的改變程度,在IHA的基礎上引入變動范圍法(RVA)(米瑋潔等,2024;Jiang et al,2014),將未來時期各指標發(fā)生頻率的75%和25%作為滿足流域生態(tài)需求的變動范圍,其改變度計算公式(Shiau et al,2007)為:

    [Dm=Nm?NεNε×100%]" " " " " " " " " " " ①

    式中:Dm為第m個IHA指標的改變程度;Nm為第m個IHA指標在變化后實際落于RVA閾值內(nèi)的觀測年數(shù);Nε為第m個IHA指標在變化后預期落于RVA閾值內(nèi)的年數(shù)。

    當Dm值為0~33%時為無改變或者低改變;33%~67%時為中度改變;67%~100%時則為高度改變。

    1.3.3" "環(huán)境流指標" "環(huán)境流指標可以表征流域流量變化過程中對生態(tài)環(huán)境的影響。環(huán)境流組成包括枯水流量、特枯流量、高流量脈沖、小洪水和大洪水5種流量事件,一共34個環(huán)境流指標(Richter,2007)。選擇歷史時期和未來時期前后環(huán)境流指標離散系數(shù)和偏差系數(shù)對嫩江流域環(huán)境流組成進行分析。

    離散系數(shù)(CD)可以反映與均值的偏離程度,其計算公式(周毅等,2017a)如下:

    [CD=H?L/M]" " " " " " " " " " " " "②

    式中:CD為離散系數(shù);H、L、M分別代表水文序列的第75、25和50百分位數(shù)。

    偏差系數(shù)(DF)可以反映兩個時期前后各指標相較于影響前的偏差(周毅等,2017a),計算公式為:

    [DF=D0?D1/D0]" " " " " " " " " " "③

    式中:[DF]為偏差系數(shù);D0和D1分別代表干擾前后各環(huán)境流指標值。

    1.3.4" "生物多樣性評價指標" "生物多樣性指標(SI)可以反映流域生態(tài)環(huán)境變化下生物多樣性和豐富度。為了使生物多樣性與IHA指標建立聯(lián)系,Yang等(2008)運用遺傳規(guī)劃法,選擇IHA指標作為水文變量,構建了生物多樣性指標與IHA中7個指標的公式(程俊翔等,2018;周毅等,2017a)。

    SI = [Dmin/Min7+DminQ3+Q5+Min3+2Max3]+Rrate" ④

    式中:SI為生物多樣性指標,其值越大表示生物多樣性越豐富;Dmin為最小流量的儒略日;Q3、Q5分別為3月和5月的月平均流量, m3/s;Min3、Min7分別為年最小3 d和7 d平均流量, m3/s;Max3為年最大3 d平均流量,m3/s;Rrate為漲水速率。

    2" "結(jié)果與分析

    2.1" "未來氣候下水文情勢變化

    利用IHA-RVA法計算了在SSP126、SSP370和SSP585情景下大賚水文站32種水文指標改變程度(圖2)。在SSP126情景模式下,嫩江流域共有10個水文指標發(fā)生低度改變,6個中度改變,16個高度改變,總改變度為69%;在SSP370情景模式下,嫩江流域共有8個水文指標發(fā)生低度改變,8個中度改變,16個高度改變,總改變度為75%;在SSP585情景模式下,嫩江流域共有11個水文指標發(fā)生低度改變,8個中度改變,13個高度改變,總改變度為66%。

    在SSP126情景模式下,嫩江流域有11個月平均流量呈增加趨勢,只有9月是下降趨勢;在SSP370和SSP585情景模式下,嫩江流域皆有10個月平均流量呈增加趨勢,只有7和10月是下降趨勢。

    2.2" "未來氣候下環(huán)境流量變化

    SSP126情景模式下,嫩江流域大洪水事件減少,枯水流量和高流量脈沖事件增多,無特枯流量事件發(fā)生;SSP370情景模式下,流域大洪水事件發(fā)生頻率變化很微小,枯水流量事件和高流量脈沖事件發(fā)生頻率增加,小洪水事件增加,沒有特枯流量事件發(fā)生;SSP585情景模式下,大洪水事件增加,且會出現(xiàn)1次極大的峰值,小洪水事件發(fā)生頻率降低,枯水流量和高流量脈沖事件發(fā)生頻率增加,沒有特枯流量事件發(fā)生(圖3)。在SSP126情景模式下,嫩江特枯流量事件為0,枯水流量、高流量脈沖、小洪水事件分別增加3.23%、0.54%、5.73%,大洪水事件減少0.89%;SSP370 模式下,特枯流量事件為 0,枯水流量增加7.56%,高流量脈沖、小洪水事件分別增加0.38%、4.39%,大洪水事件減少1.33%。SSP585 模式下,特枯流量事件為0,枯水流量增加 8.01%,高流量脈沖、小洪水、大洪水事件分別增加1.02%、0.25%和0.42%。

    由IHA及偏離系數(shù)、偏差系數(shù)計算可知(圖3),未來氣候下,環(huán)境流趨于單一,流量變化范圍縮小,枯水流量和高流量脈沖事件集中。SSP126模式下,1、4、5、10月的枯水和高流量脈沖事件的頻率及持續(xù)時間變化較大;SSP370模式下,1、5、10、11、12月的相關指標變化較大;SSP585模式下,1、5、10、11月的相關指標變化較大。

    2.3" "未來氣候下生物多樣性變化

    SPP126、SSP370和SSP585情景下,嫩江流域SI平均值分別為0.39、0.36和0.37,高于歷史時期平均值(0.29)。在3種情景模式下,受氣候徑流影響,嫩江流域生物多樣性指數(shù)波動較大,均值都比歷史時期大,但總體增長幅度較小(圖4)。從SI平均指數(shù)看,3種情景模式對嫩江流域的影響程度為SSP126gt;SSP585gt;SSP370(圖5)。

    3" "討論

    氣候變化會影響河道生態(tài)流量,進而影響生物多樣性的復雜性和豐富性。在未來氣候的3種情景模式下,嫩江流域水文指標都發(fā)生了50%以上程度的改變,月平均流量減少,來水不確定性增大,嫩江流域生態(tài)流量的穩(wěn)定性降低。周毅等(2017b)認為,黃河流域水文情勢指標發(fā)生較大改變,會對流域生態(tài)系統(tǒng)產(chǎn)生不利影響。在這種情境下,水文情勢變化對水生生物的影響也將進一步增強。水質(zhì)、水溫等條件會隨水文情勢改變,進而影響流域水環(huán)境狀況,影響水生生物的生存和繁殖(李建等,2013;Gooseff et al,2017)。因此,未來氣候變化影響下,嫩江流域的水生生態(tài)環(huán)境將發(fā)生顯著變化,進一步改變嫩江流域河道生物多樣性的組成。

    環(huán)境流組成的變化會影響生物棲息地環(huán)境,尤其是對處于食物鏈頂端的魚類影響最為明顯。在未來環(huán)境3種不同情景模式下,嫩江流域的環(huán)境流組成發(fā)生了明顯的變化。其中,枯水流量和高流量脈沖事件的頻率和持續(xù)時間變化程度較大。這是由于降水的持續(xù)時間、頻次和強度在很大的程度上決定了區(qū)域或流域河道流量的大小及其演變。Wu等(2022;2023)發(fā)現(xiàn),未來氣候變化下嫩江流域極端降水事件的頻次和強度將大幅增加,其中水文干旱和洪水的頻次、持續(xù)時間也會進一步增加。因此,不斷增多的極端水文事件可能是導致枯水流量和高流量脈沖事件的頻率與持續(xù)時間發(fā)生明顯變化的主要原因。環(huán)境流組成趨于單一化會對下游生態(tài)環(huán)境產(chǎn)生不利影響,阻礙魚類產(chǎn)卵 (魏卿等,2019)。此外,流域水生生物的健康也與環(huán)境流事件的發(fā)生時間、頻率和持續(xù)時間緊密相關(王鴻翔等,2019;方春明等,2014)。因此,嫩江流域環(huán)境流中枯水流量事件和高流量脈沖事件增多,會降低流域水域溶解營養(yǎng)物質(zhì)的能力,進而改變水生生物的生態(tài)環(huán)境,影響水體營養(yǎng)物質(zhì)的遷移和轉(zhuǎn)化,最終影響魚類的繁殖和棲息。

    在未來環(huán)境3種不同情景模式下,嫩江流域的生物多樣性指數(shù)增加,但增幅較小。氣候變化對生物多樣性影響顯著,涉及物種的分布、生理、行為、物候、遺傳、入侵和滅絕等方面(Bellard et al,2012)。此外,嫩江流域的物種豐富度和均勻度將提高,可能導致流域的物種組成發(fā)生變化。然而,由于生態(tài)流量的改變,對氣候變化敏感或不適應的物種會大幅度減少或消失,而適應或抗逆性強的物種則會增加或擴散,最終改變嫩江流域生物多樣性的格局。

    科學合理的水文調(diào)控可以在很大程度上緩解氣候變化導致的水文指標改變所帶來的負面影響。一方面,借助現(xiàn)有的基礎設施(如干流和支流的水庫大壩),充分發(fā)揮其水文調(diào)蓄功能,優(yōu)化水資源調(diào)度,最大程度保障維系河道生態(tài)系統(tǒng)健康的生態(tài)流量;另一方面,充分運用基于自然的水資源解決方案(如利用大量的沼澤濕地)(Wu et al,2021),依托國家河湖復蘇工程,實現(xiàn)基于綠色-灰色基礎設施的嫩江流域水資源綜合管理,滿足流域生態(tài)系統(tǒng)健康發(fā)展所需的水文情勢。

    探究未來氣候變化下的流域生態(tài)流量變化,并對生物多樣性進行預估,可為流域的水資源管理和水生態(tài)系統(tǒng)保護提供科學依據(jù),具有一定的現(xiàn)實意義。本文對未來氣候的變化情景僅基于HYDROTEL模型模擬的未來日均流量,較少考慮極端或突發(fā)氣候情景,這些問題對生態(tài)流量和生物多樣性的影響值得進一步研究和分析。此外,SI指數(shù)主要用于表征魚類的生物多樣性,無法涵蓋微生物、藻類、水生高等植物、各種無脊椎動物和脊椎動物等其他水生生物。因此,后續(xù)需要采用實際采樣與數(shù)值模型相結(jié)合的方法,系統(tǒng)分析生態(tài)流量變化對水生生物多樣性的影響機制。

    參考文獻

    邴建平, 鄧鵬鑫, 吳智, 等, 2023. 贛江流域生態(tài)流量與地表水資源可利用量研究[J]. 人民長江, 54(2): 127-131, 170.

    BING J P, DENG P X, WU Z, et al, 2023. Ecological flow and available surface water resources in Ganjiang River Basin[J]. Yangtze River, 54(2): 127-131, 170.

    程俊翔, 徐力剛, 姜加虎, 2018. 生態(tài)最相關水文指標的優(yōu)選及其在洞庭湖環(huán)境流量估算中的應用[J]. 湖泊科學, 30(5): 1235-1245.

    CHENG J X, XU L G, JIANG J H, 2018. Optimal selection of the most ecologically relevant hydrologic indicators and its application for environmental flow calculation in Lake Dongting[J]. Journal of Lake Sciences, 30(5): 1235-1245.

    方春明, 胡春宏, 陳緒堅, 2014. 三峽水庫運用對荊江三口分流及洞庭湖的影響[J]. 水利學報, 45(1): 36-41.

    FANG C M, HU C H, CHEN X J, 2014. Impacts of Three Georges Reservoir’s operation on outflow of the three outlets of Jingjiang River and Dongting Lake[J]. Journal of Hydraulic Engineering, 45(1): 36-41.

    李建, 夏自強, 戴會超, 等, 2013. 三峽初期蓄水對典型魚類棲息地適宜性的影響[J]. 水利學報, 44(8): 892-900.

    LI J, XIA Z Q, DAI H C, et al, 2013. Effect of the Three Gorges Reservoir initial filling on downstream habitat suitability of the typical fishes[J]. Journal of Hydraulic Engineering, 44(8): 892-900.

    李澤君, 黃本勝, 邱靜, 等, 2021. 變化環(huán)境下韓江生態(tài)流量演變特征分析[J]. 水資源保護, 37(5): 22-29.

    LI Z J, HUANG B S, QIU J, et al, 2021. Analysis on evolution characteristics of ecological flow of Hanjiang River under changing environment[J]. Water Resources Protection, 37(5): 22-29.

    劉綠柳, 魏麟驍, 徐影, 等, 2021. 氣候變化對黃河流域生態(tài)徑流影響預估[J]. 水科學進展, 32(6): 824-833.

    LIU L L, WEI L X, XU Y, et al, 2021. Projection of climate change impacts on ecological flow in the Yellow River Basin[J]. Advances in Water Science, 32(6): 824-833.

    米瑋潔, 張舜, 胡俊, 等, 2024. 黃河內(nèi)蒙古河段水文情勢變化及生態(tài)效應初探[J]. 水生態(tài)學雜志, 45(1): 42-49.

    MI W J, ZHANG S, HU J, et al, 2024. Ecological effects of flow regime alteration in Inner Mongolia reach of the Yellow River[J]. Journal of Hydroecology, 45(1): 42-49.

    田震, 敖偲成, 李先福, 等, 2023. 高黎貢山南段河流生態(tài)系統(tǒng)健康評價[J]. 水生態(tài)學雜志, 44(1): 25-33.

    TIAN Z, AO S C, LI X F, et al, 2023. Health assessment of river ecosystem of the southern gaoligong mountains[J]. Journal of Hydroecology, 44(1): 25-33.

    王大雙, 鄒進, 劉磊, 等, 2023. 基于改進RVA法的云龍水庫流域生態(tài)調(diào)度研究[J]. 水文, 43(1): 41-46.

    WANG D S, ZOU J, LIU L, et al, 2023. Study on ecological regulation of Yunlong Reservoir basin based on improved RVA method[J]. Journal of China Hydrology, 43(1): 41-46.

    王鴻翔, 李萌萌, 馬志軍, 等, 2019. 長江荊江三口環(huán)境流變化研究[J]. 水力發(fā)電, 45(5): 1-6.

    王凱, 王聰, 馮曉明, 等, 2022. 生物多樣性與生態(tài)系統(tǒng)多功能性的關系研究進展[J]. 生態(tài)學報, 42(1): 11-23.

    WANG K, WANG C, FENG X M, et al, 2022. Research progress on the relationship between biodiversity and ecosystem multifunctionality[J]. Acta Ecologica Sinica, 42(1): 11-23.

    魏卿, 薛聯(lián)青, 張敏, 等, 2019. 淮河流域環(huán)境流變化及其對洪澤湖魚類棲息地的生態(tài)影響[J]. 水資源保護, 35(4): 89-94.

    WEI Q, XUE L Q, ZHANG M, et al, 2019. Changes of environmental flow in Huaihe River Basin and its ecological impact on fish habitat in Hongze Lake[J]. Water Resources Protection, 35(4): 89-94.

    吳燕鋒, 章光新, 齊鵬, 等, 2019. 耦合濕地模塊的流域水文模型模擬效率評價[J]. 水科學進展, 30(3): 326-336.

    WU Y F, ZHANG G X, QI P, et al, 2019. Integration of wetland modules into the watershed hydrological model: assessment of simulation accuracy[J]. Advances in Water Science, 30(3): 326-336.

    張輝, 曾晨軍, 李婷, 等, 2022. 基于四大家魚產(chǎn)卵需求的漢江中下游生態(tài)流量研究[J]. 水生態(tài)學雜志, 43(3): 1-8.

    ZHANG H, ZENG C J, LI T, et al, 2022. Ecological flow in the mid-lower Hanjiang River based on spawning demands of the four major Chinese carps[J]. Journal of Hydroecology, 43(3): 1-8.

    張麗霞, 陳曉龍, 辛曉歌, 2019. CMIP6情景模式比較計劃(ScenarioMIP)概況與評述[J]. 氣候變化研究進展, 15(5): 519-525.

    ZHANG L X, CHEN X L, XIN X G, 2019. Short commentary on CMIP6 scenario model intercomparison project(ScenarioMIP)[J]. Climate Change Research, 15(5): 519-525.

    趙宗慈, 羅勇, 黃建斌, 2016. CMIP6的設計[J]. 氣候變化研究進展, 12(3): 258-260.

    ZHAO Z C, LUO Y, HUANG J B, 2016. Design of cmip6[J]. Climate Change Research, 12(3): 258-260.

    周天軍, 鄒立維, 陳曉龍, 2019. 第六次國際耦合模式比較計劃(CMIP6)評述[J]. 氣候變化研究進展, 15(5): 445-456.

    ZHOU T J, ZOU L W, CHEN X L, 2019. Commentary on the coupled model intercomparison project phase 6(CMIP6)[J]. Climate Change Research, 15(5): 445-456.

    周毅, 崔同, 高滿, 等, 2017. 考慮不同水文年及IHA指標相關性的水文特征評估方法[J]. 水文, 37(2): 20-25.

    ZHOU Y, CUI T, GAO M, et al, 2017. A new method to assess hydrological regime alterations considering hydrological year types and correlation among IHA[J]. Journal of China Hydrology, 37(2): 20-25.

    周毅, 崔同, 鄭鑫, 等, 2017. 基于IHA的黃河源區(qū)水文情勢及環(huán)境流變化[J]. 人民黃河, 39(7): 61-64, 70.

    ZHOU Y, CUI T, ZHENG X, et al, 2017. Hydrological regime and environmental flow change in the source region of Yellow River based on IHA[J]. Yellow River, 39(7): 61-64, 70.

    ACREMAN M C, DUNBAR M J, 2004. Defining environmental river flow requirements-a review[J]. Hydrology and Earth System Sciences, 8(5): 861-876.

    BELLARD C, BERTELSMEIER C, LEADLEY P, et al, 2012. Impacts of climate change on the future of biodiversity[J]. Ecology Letters, 15(4): 365-377.

    CHEN Y D, YANG T, XU C Y, et al, 2010. Hydrologic alteration along the middle and upper east river (Dongjiang) basin, South China: a visually enhanced mining on the results of RVA method[J]. Stochastic Environmental Research and Risk Assessment, 24(1): 9-18.

    GOOSEFF M N, WLOSTOWSKI A, MCKNIGHT D M, et al, 2017. Hydrologic connectivity and implications for ecosystem processes - Lessons from naked watersheds[J]. Geomorphology, 277: 63-71.

    HUNTINGTON T G, 2006. Evidence for intensification of the global water cycle: review and synthesis[J]. Journal of Hydrology, 319(1/2/3/4): 83-95.

    JIANG L Z, BAN X, WANG X L, et al, 2014. Assessment of hydrologic alterations caused by the Three Gorges dam in the middle and lower reaches of Yangtze River, China[J]. Water, 6(5): 1419-1434.

    POFF N L, MATTHEWS J H, 2013. Environmental flows in the anthropocence: past progress and future prospects[J]. Current Opinion in Environmental Sustainability, 5(6): 667-675.

    RICHTER B D, BAUMGARTNER J V, BRAUN D P, et al, 1998. A spatial assessment of hydrologic alteration within a river network[J]. Regulated Rivers: Research amp; Management, 14(4): 329-340.

    RICHTER B D, THOMAS G A, 2007. Restoring environmental flows by modifying dam operations[J]. Ecology and Society, 12: art12.

    SHIAU J T, WU F C, 2007. Pareto-optimal solutions for environmental flow schemes incorporating the intra-annual and interannual variability of the natural flow regime[J]. Water Resources Research, 43(6): W06433.

    SUN G, HALLEMA D, ASBJORNSEN H, 2017. Ecohydrological processes and ecosystem services in the anthropocene: a review[J]. Ecological Processes, 6(1): 35.

    WU Y F, SUN J X, JUN XU Y, et al, 2022. Projection of future hydrometeorological extremes and wetland flood mitigation services with different global warming levels: a case study in the Nenjiang river basin[J]. Ecological Indicators, 140: 108987.

    WU Y F, ZHANG G X, XU Y J, et al, 2021. River damming reduces wetland function in regulating flow[J]. Journal of Water Resources Planning and Management, 147(10): 05021014.

    WU Y F, SUN J X, HU B T, et al, 2023. Can the combining of wetlands with reservoir operation reduce the risk of future floods and droughts?[J]. Hydrology and Earth System Sciences, 27(14): 2725-2745.

    YANG Y E, CAI X M, HERRICKS E E, 2008. Identification of hydrologic indicators related to fish diversity and abundance: a data mining approach for fish community analysis[J]. Water Resources Research, 44(4): W04412.

    ZHANG Q, XIAO M Z, LIU C L, et al, 2014. Reservoir-induced hydrological alterations and environmental flow variation in the East River, the Pearl River basin, China[J]. Stochastic Environmental Research and Risk Assessment, 28(8): 2119-2131.

    (責任編輯" "鄭金秀)

    Changes in the Ecological Flow and Biodiversity of Nenjiang River Basin Under Future Climate Scenarios

    MA Qiusheng1,2, CHEN Mo1,2, SUN Jiahao1,2, ZHANG Qingsong1,2, WU Yanfeng3

    (1. School of Hydraulic and Electric Power, Heilongjiang University, Harbin" "150080, P.R. China;

    2. Institute of Groundwater in Cold Region, Heilongjiang, Harbin 150080, P.R. China;

    3. Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences,Changchun" "130102, P.R. China)

    Abstract:Climate change alters runoff, which subsequently effects ecological flow and biodiversity. As a result, the assessment of climate change impacts on ecological flow and biodiversity in river basins is now required for protecting river basin ecosystems. However, the extent of these changes remains largely unknown and has increasingly attracted the attention of researchers at home and abroad. In this study, the Nenjiang River basin was selected for study, and we comprehensively analyzed changes in ecological flow, environmental flow and biodiversity over the 75 year period 2025-2100 using the hydrological variability index method. Based on daily runoff data from Dalai Hydrological Station from 1963 to 2018, our calibrated HYDROTEL model was used to simulate the daily runoff under the CMIP6 climate scenarios (SSP126, SSP370 and SSP585). Changes in the indicators of hydrological alteration, environmental flow, and biodiversity were then calculated and used to comprehensively assess the changes in ecological flow and the evolution of biodiversity in the Nenjiang River basin under the three future climate scenarios. The hydrological indicators changed markedly under the future climate scenarios compared with the historical period: the uncertainty of water inflow increased, and the stability of ecological flows decreased. The changes in hydrological indicators were 69% under scenario SSP126 and 75% under scenario SSP370. The frequency of low-flow events and high pulse flow events of ambient streamflow increased by 3.23% and 0.54% under scenario SSP126, by 7.56% and 0.38% under scenario SSP370, and by 8.01% and 1.02% under scenario SSP585 scenario. Most of the flow processes in Nenjiang River basin were classified as dry flow and high flow pulse events, and the range of flow change decreased, presenting a gradual homogenization of environmental flow. Compared to the historical period, with a mean biodiversity index of 0.29, biodiversity in the Nenjiang River basin increased under future climate scenarios, with mean biodiversity indices of 0.39, 0.36 and 0.37, respectively, under climate scenarios SPP126, SSP137 and SSP585. This study provides a scientific basis for water resource management and aquatic ecosystem protection in the Nenjiang River basin as climate changes.

    Key words: ecological flow; environmental flow; climate change; biodiversity; Nenjiang River basin

    基金項目:國家重點研發(fā)計劃(2021YFC3200203);國家自然科學基金(421010514,42207088)。

    作者簡介:馬秋生,1997年生,男,碩士研究生,主要從事水文與水資源方面研究。E-mail:846842247@qq.com

    通信作者:陳末,1989年生,博士,副教授,主要從事水文與水資源方面研究。E-mail:chenmocc12@163.com

    亚洲人与动物交配视频| 精华霜和精华液先用哪个| 日本黄色片子视频| 亚洲av成人av| 一区二区三区激情视频| 别揉我奶头 嗯啊视频| 国产av不卡久久| 午夜福利高清视频| 国产精品一区二区性色av| 国产成人a区在线观看| 最好的美女福利视频网| 又爽又黄a免费视频| 午夜免费激情av| 能在线免费观看的黄片| 桃色一区二区三区在线观看| 级片在线观看| 免费av观看视频| 午夜激情福利司机影院| 最近视频中文字幕2019在线8| 国产女主播在线喷水免费视频网站 | 国产女主播在线喷水免费视频网站 | 我的女老师完整版在线观看| a级毛片免费高清观看在线播放| 九九热线精品视视频播放| 欧美zozozo另类| 精品人妻熟女av久视频| 亚洲精品456在线播放app | 成年版毛片免费区| 国内精品一区二区在线观看| 午夜免费男女啪啪视频观看 | 欧美另类亚洲清纯唯美| xxxwww97欧美| 热99re8久久精品国产| 99riav亚洲国产免费| 69av精品久久久久久| 日韩欧美免费精品| 舔av片在线| 亚洲熟妇熟女久久| 国产精品,欧美在线| 男女做爰动态图高潮gif福利片| 亚洲精品一区av在线观看| 免费观看在线日韩| 在线观看66精品国产| 国产精品一区二区三区四区免费观看 | 亚洲国产精品成人综合色| 国产伦人伦偷精品视频| 男女之事视频高清在线观看| 制服丝袜大香蕉在线| 欧美性感艳星| 欧美一区二区亚洲| 男女之事视频高清在线观看| 日本免费a在线| 天美传媒精品一区二区| 精品免费久久久久久久清纯| 搡老妇女老女人老熟妇| 久久九九热精品免费| 国产日本99.免费观看| 天堂影院成人在线观看| 成人欧美大片| 精品福利观看| 又粗又爽又猛毛片免费看| 99久久精品热视频| 丰满乱子伦码专区| 少妇猛男粗大的猛烈进出视频 | 黄色女人牲交| 2021天堂中文幕一二区在线观| 日本黄大片高清| 亚洲经典国产精华液单| 中国美白少妇内射xxxbb| 精品国内亚洲2022精品成人| 午夜久久久久精精品| 欧美成人免费av一区二区三区| x7x7x7水蜜桃| 精品欧美国产一区二区三| 男女视频在线观看网站免费| 噜噜噜噜噜久久久久久91| 长腿黑丝高跟| 欧美三级亚洲精品| 不卡视频在线观看欧美| 51国产日韩欧美| 欧美色视频一区免费| www.www免费av| 免费人成视频x8x8入口观看| 欧美性感艳星| 嫩草影院精品99| 女生性感内裤真人,穿戴方法视频| 免费不卡的大黄色大毛片视频在线观看 | 嫩草影视91久久| 97碰自拍视频| 神马国产精品三级电影在线观看| 日韩精品中文字幕看吧| 中文字幕久久专区| 免费人成视频x8x8入口观看| 国产精品,欧美在线| 神马国产精品三级电影在线观看| 亚洲第一区二区三区不卡| 欧美一区二区精品小视频在线| 亚洲国产高清在线一区二区三| 亚洲中文字幕日韩| 午夜福利在线观看免费完整高清在 | 日本-黄色视频高清免费观看| 亚洲人成网站在线播放欧美日韩| 久久精品国产亚洲av涩爱 | 国产精品,欧美在线| 日韩中字成人| avwww免费| 99久久精品热视频| 88av欧美| 免费高清视频大片| 亚洲内射少妇av| 琪琪午夜伦伦电影理论片6080| 国产一区二区在线av高清观看| 丰满乱子伦码专区| 亚洲欧美日韩卡通动漫| 一本久久中文字幕| 久久久国产成人免费| 国产成人av教育| 女的被弄到高潮叫床怎么办 | 久久欧美精品欧美久久欧美| 深夜a级毛片| 成人国产麻豆网| 国产精品,欧美在线| 欧美一区二区国产精品久久精品| 国产黄色小视频在线观看| 色精品久久人妻99蜜桃| 日本一二三区视频观看| 能在线免费观看的黄片| 女人被狂操c到高潮| 老女人水多毛片| 欧美高清成人免费视频www| 麻豆成人av在线观看| 12—13女人毛片做爰片一| 波多野结衣高清无吗| 国产视频一区二区在线看| 国产午夜精品论理片| 在线观看舔阴道视频| 成人精品一区二区免费| 国产高清三级在线| 色av中文字幕| 日本撒尿小便嘘嘘汇集6| 好男人在线观看高清免费视频| 麻豆国产av国片精品| 国产精品久久视频播放| www.www免费av| aaaaa片日本免费| 在线播放国产精品三级| 国产aⅴ精品一区二区三区波| 国产综合懂色| 又黄又爽又刺激的免费视频.| 男人狂女人下面高潮的视频| 女人被狂操c到高潮| 亚洲av不卡在线观看| 蜜桃久久精品国产亚洲av| 国产av一区在线观看免费| 老师上课跳d突然被开到最大视频| 国产 一区 欧美 日韩| 欧美一级a爱片免费观看看| 亚洲av中文字字幕乱码综合| 中文在线观看免费www的网站| 精品久久久久久成人av| 久久精品国产99精品国产亚洲性色| 乱系列少妇在线播放| 国产精华一区二区三区| 亚洲精华国产精华精| 国产免费视频播放在线视频| 97热精品久久久久久| 国产在视频线精品| 久久久久性生活片| 内地一区二区视频在线| tube8黄色片| 久热这里只有精品99| 91狼人影院| 亚洲色图av天堂| 欧美日韩一区二区视频在线观看视频在线| 国产女主播在线喷水免费视频网站| 亚洲精品久久久久久婷婷小说| 3wmmmm亚洲av在线观看| 日韩三级伦理在线观看| 我要看日韩黄色一级片| 女人十人毛片免费观看3o分钟| 久久国内精品自在自线图片| 水蜜桃什么品种好| 91精品国产九色| 国产日韩欧美在线精品| 91精品一卡2卡3卡4卡| 精品熟女少妇av免费看| 国产 一区精品| 国产真实伦视频高清在线观看| 人妻夜夜爽99麻豆av| 欧美xxxx黑人xx丫x性爽| 视频中文字幕在线观看| 精品久久久久久久末码| 国产人妻一区二区三区在| 免费看日本二区| tube8黄色片| 91狼人影院| 大陆偷拍与自拍| 大话2 男鬼变身卡| 亚洲美女黄色视频免费看| 黑人高潮一二区| 亚洲高清免费不卡视频| 特大巨黑吊av在线直播| 人人妻人人看人人澡| 身体一侧抽搐| xxx大片免费视频| 女人十人毛片免费观看3o分钟| 毛片女人毛片| 亚洲人与动物交配视频| 特大巨黑吊av在线直播| 中文字幕人妻熟人妻熟丝袜美| 80岁老熟妇乱子伦牲交| 欧美+日韩+精品| 网址你懂的国产日韩在线| 日本-黄色视频高清免费观看| 亚洲精品视频女| 最新中文字幕久久久久| 亚洲av.av天堂| 国国产精品蜜臀av免费| 色视频在线一区二区三区| 久久国产亚洲av麻豆专区| 女性被躁到高潮视频| 亚洲成人一二三区av| 中文在线观看免费www的网站| 国产精品国产av在线观看| 三级经典国产精品| 欧美精品人与动牲交sv欧美| 亚洲天堂av无毛| 一级毛片 在线播放| 一级二级三级毛片免费看| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩另类电影网站 | 国产深夜福利视频在线观看| 午夜免费观看性视频| 在线精品无人区一区二区三 | 卡戴珊不雅视频在线播放| 久热久热在线精品观看| 三级经典国产精品| 国产成人免费观看mmmm| 大码成人一级视频| 亚洲精品国产色婷婷电影| 卡戴珊不雅视频在线播放| 国产 精品1| 丝瓜视频免费看黄片| 小蜜桃在线观看免费完整版高清| 成人高潮视频无遮挡免费网站| 熟女电影av网| 精品人妻视频免费看| 国产亚洲欧美精品永久| 亚洲色图综合在线观看| 久久久色成人| 欧美 日韩 精品 国产| 国产一区二区三区综合在线观看 | 99久久人妻综合| 高清在线视频一区二区三区| 老师上课跳d突然被开到最大视频| 国产一区有黄有色的免费视频| 国产在线男女| 亚洲欧美日韩东京热| 99视频精品全部免费 在线| 纵有疾风起免费观看全集完整版| freevideosex欧美| 精品国产乱码久久久久久小说| 少妇人妻久久综合中文| 国产亚洲欧美精品永久| 国语对白做爰xxxⅹ性视频网站| 亚洲成人一二三区av| 国产成人一区二区在线| av在线蜜桃| 亚洲欧洲国产日韩| 欧美日韩视频高清一区二区三区二| 国产精品久久久久成人av| 女的被弄到高潮叫床怎么办| 成人一区二区视频在线观看| 黄色日韩在线| 亚洲真实伦在线观看| 国产亚洲午夜精品一区二区久久| 97热精品久久久久久| 91精品国产九色| 18+在线观看网站| 久久久久久久久久人人人人人人| 91aial.com中文字幕在线观看| 一级a做视频免费观看| 日本一二三区视频观看| 亚洲av综合色区一区| h视频一区二区三区| 在线天堂最新版资源| 欧美精品一区二区免费开放| 免费av中文字幕在线| 久久久色成人| 久久久a久久爽久久v久久| 国产黄片美女视频| 日本午夜av视频| 久久影院123| 欧美精品亚洲一区二区| 成人美女网站在线观看视频| 纵有疾风起免费观看全集完整版| 久久久国产一区二区| 精品久久久精品久久久| 成年av动漫网址| 国产成人freesex在线| freevideosex欧美| 人人妻人人看人人澡| 午夜福利网站1000一区二区三区| 一区二区av电影网| 丰满少妇做爰视频| 久久久精品免费免费高清| 亚洲成人中文字幕在线播放| 在线观看免费高清a一片| 成人美女网站在线观看视频| 国产黄频视频在线观看| 欧美xxxx黑人xx丫x性爽| 99久久人妻综合| 午夜福利在线观看免费完整高清在| 亚洲精品国产成人久久av| 在线亚洲精品国产二区图片欧美 | 毛片女人毛片| 日韩欧美一区视频在线观看 | av在线观看视频网站免费| 欧美zozozo另类| 国产 一区 欧美 日韩| 男女啪啪激烈高潮av片| 亚洲,一卡二卡三卡| 国产视频内射| 伦理电影大哥的女人| 亚洲精品乱久久久久久| 99热6这里只有精品| 日本黄色日本黄色录像| 精品久久久久久电影网| 在线免费十八禁| 亚洲,欧美,日韩| 免费播放大片免费观看视频在线观看| 日韩欧美精品免费久久| 成年免费大片在线观看| 欧美变态另类bdsm刘玥| 欧美精品人与动牲交sv欧美| 亚洲自偷自拍三级| 中文字幕精品免费在线观看视频 | 久久精品久久久久久久性| 国产成人午夜福利电影在线观看| 久久久久久九九精品二区国产| av又黄又爽大尺度在线免费看| 一本久久精品| av在线app专区| 亚洲成人av在线免费| 精品久久久久久电影网| 欧美成人a在线观看| 亚洲精品中文字幕在线视频 | 18禁在线无遮挡免费观看视频| 成人二区视频| 内射极品少妇av片p| 男人舔奶头视频| 自拍欧美九色日韩亚洲蝌蚪91 | 免费久久久久久久精品成人欧美视频 | 成人美女网站在线观看视频| 亚洲成色77777| 亚洲国产精品专区欧美| 毛片一级片免费看久久久久| 久久久久久人妻| 欧美最新免费一区二区三区| www.色视频.com| 乱系列少妇在线播放| 亚洲av中文字字幕乱码综合| 久久久亚洲精品成人影院| 色吧在线观看| 王馨瑶露胸无遮挡在线观看| 久久精品国产鲁丝片午夜精品| 色网站视频免费| 日本vs欧美在线观看视频 | 大片电影免费在线观看免费| 女人十人毛片免费观看3o分钟| 久热久热在线精品观看| 国产成人freesex在线| 性高湖久久久久久久久免费观看| 一区二区av电影网| 美女cb高潮喷水在线观看| 黄色欧美视频在线观看| 男的添女的下面高潮视频| 噜噜噜噜噜久久久久久91| 日韩一区二区三区影片| 国产精品99久久久久久久久| 精品国产乱码久久久久久小说| 天天躁日日操中文字幕| 日日啪夜夜爽| 日本-黄色视频高清免费观看| 久久久久视频综合| 亚洲欧美中文字幕日韩二区| 午夜免费观看性视频| 久久99精品国语久久久| 视频中文字幕在线观看| 99热这里只有精品一区| 一级毛片黄色毛片免费观看视频| 免费观看性生交大片5| 水蜜桃什么品种好| 成人高潮视频无遮挡免费网站| 高清日韩中文字幕在线| 秋霞在线观看毛片| 嘟嘟电影网在线观看| 亚洲久久久国产精品| 日本wwww免费看| 国内精品宾馆在线| 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩无卡精品| av卡一久久| 精品久久久噜噜| 亚洲精品亚洲一区二区| 久久久久久久精品精品| 欧美成人午夜免费资源| 女人久久www免费人成看片| 丰满人妻一区二区三区视频av| 99热全是精品| 三级经典国产精品| 成人特级av手机在线观看| 干丝袜人妻中文字幕| av天堂中文字幕网| 国产国拍精品亚洲av在线观看| 亚洲精品成人av观看孕妇| 日韩,欧美,国产一区二区三区| 一本色道久久久久久精品综合| 晚上一个人看的免费电影| 久久热精品热| 亚洲精品日本国产第一区| 欧美+日韩+精品| 久久99热这里只有精品18| 亚洲va在线va天堂va国产| 国产精品福利在线免费观看| 一级毛片电影观看| 性色av一级| 夫妻性生交免费视频一级片| 一级片'在线观看视频| 超碰97精品在线观看| .国产精品久久| 噜噜噜噜噜久久久久久91| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 高清av免费在线| 一级av片app| 1000部很黄的大片| 国产亚洲午夜精品一区二区久久| 99热这里只有是精品在线观看| 久久女婷五月综合色啪小说| 国产男女内射视频| 久热久热在线精品观看| 午夜福利在线观看免费完整高清在| 久久综合国产亚洲精品| 精品亚洲成a人片在线观看 | 日韩欧美精品免费久久| 成人午夜精彩视频在线观看| 99热国产这里只有精品6| 久久久久网色| 亚洲欧美日韩东京热| 啦啦啦在线观看免费高清www| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 一本色道久久久久久精品综合| 五月天丁香电影| 国产一区二区三区综合在线观看 | 婷婷色综合大香蕉| 少妇的逼水好多| 中国美白少妇内射xxxbb| 街头女战士在线观看网站| 日韩欧美精品免费久久| 91精品一卡2卡3卡4卡| 久久久久久久久大av| 国产精品国产三级国产专区5o| 一级毛片黄色毛片免费观看视频| 国产综合精华液| 一级毛片aaaaaa免费看小| 99热国产这里只有精品6| 最后的刺客免费高清国语| 成年av动漫网址| 亚洲精品aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 另类亚洲欧美激情| 国产精品福利在线免费观看| 丰满乱子伦码专区| 亚洲图色成人| 一级毛片黄色毛片免费观看视频| 18禁裸乳无遮挡动漫免费视频| 亚洲成人一二三区av| 寂寞人妻少妇视频99o| 高清不卡的av网站| 人妻 亚洲 视频| 亚洲一级一片aⅴ在线观看| 成人二区视频| 午夜福利网站1000一区二区三区| 天天躁日日操中文字幕| 蜜桃在线观看..| 爱豆传媒免费全集在线观看| 久久久久久久国产电影| 亚洲精品乱码久久久久久按摩| videossex国产| 国产毛片在线视频| 国产精品三级大全| 直男gayav资源| 自拍偷自拍亚洲精品老妇| 一本一本综合久久| 国产精品一区二区三区四区免费观看| 黑人猛操日本美女一级片| 日日啪夜夜爽| 国产 一区精品| 欧美精品亚洲一区二区| 精品国产露脸久久av麻豆| 久久精品久久精品一区二区三区| 午夜免费观看性视频| 日本午夜av视频| 97精品久久久久久久久久精品| 99视频精品全部免费 在线| 狂野欧美激情性xxxx在线观看| 少妇的逼水好多| 国产欧美另类精品又又久久亚洲欧美| 一级黄片播放器| 人体艺术视频欧美日本| 亚洲国产精品一区三区| av在线老鸭窝| 亚洲成人中文字幕在线播放| 老司机影院成人| 黄色一级大片看看| 免费av中文字幕在线| 亚洲美女搞黄在线观看| 日本黄大片高清| 国精品久久久久久国模美| 在线观看免费视频网站a站| 搡女人真爽免费视频火全软件| 亚洲人与动物交配视频| 成人亚洲精品一区在线观看 | 国国产精品蜜臀av免费| 亚洲美女搞黄在线观看| 久久久午夜欧美精品| 一本久久精品| 国产成人freesex在线| 亚洲国产欧美人成| 亚洲精品一区蜜桃| 亚洲欧美精品专区久久| 日韩亚洲欧美综合| 看十八女毛片水多多多| 久久影院123| 青春草视频在线免费观看| 少妇人妻 视频| 国产精品一二三区在线看| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 久久久精品免费免费高清| 国产成人91sexporn| 免费看日本二区| 这个男人来自地球电影免费观看 | 亚洲aⅴ乱码一区二区在线播放| 国产在线一区二区三区精| 免费人妻精品一区二区三区视频| 色婷婷久久久亚洲欧美| 亚洲av二区三区四区| 美女xxoo啪啪120秒动态图| 99热这里只有是精品50| 久久国内精品自在自线图片| 国产成人精品福利久久| 99九九线精品视频在线观看视频| 美女中出高潮动态图| 美女国产视频在线观看| 特大巨黑吊av在线直播| 这个男人来自地球电影免费观看 | 国产欧美亚洲国产| 国产伦精品一区二区三区视频9| 99热全是精品| av.在线天堂| 欧美国产精品一级二级三级 | 黄色视频在线播放观看不卡| 精品国产露脸久久av麻豆| 免费观看无遮挡的男女| 久久av网站| 91久久精品国产一区二区三区| 久久久久久久亚洲中文字幕| 久久99热这里只有精品18| 如何舔出高潮| 国产精品福利在线免费观看| 亚洲人与动物交配视频| 亚洲av日韩在线播放| 王馨瑶露胸无遮挡在线观看| 国产熟女欧美一区二区| 国产精品嫩草影院av在线观看| 我的老师免费观看完整版| av不卡在线播放| 亚洲成人av在线免费| av在线播放精品| 最近中文字幕2019免费版| 亚洲国产最新在线播放| 天堂中文最新版在线下载| 国产 精品1| 国产真实伦视频高清在线观看| 久久久久久久国产电影| 亚洲欧美成人综合另类久久久| 久久久久久久久久久丰满| 亚洲av成人精品一二三区| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 国产精品一二三区在线看| 街头女战士在线观看网站| 麻豆乱淫一区二区| 免费高清在线观看视频在线观看| 麻豆精品久久久久久蜜桃| 亚洲精品,欧美精品| 国产精品国产三级国产专区5o| 亚洲精品自拍成人| 最近2019中文字幕mv第一页| 在线观看av片永久免费下载| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩东京热| 午夜免费鲁丝| 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 亚洲无线观看免费| 亚洲国产精品一区三区| 国产成人一区二区在线| 午夜免费鲁丝| 国产精品99久久久久久久久| 日韩成人av中文字幕在线观看| 不卡视频在线观看欧美|