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

    基于多水文改變指標評價東江流域河流流態(tài)變化及其對生物多樣性的影響

    2016-11-28 08:56:09顧西輝孔冬冬劉劍宇
    生態(tài)學報 2016年19期
    關鍵詞:東江流態(tài)赤字

    顧西輝,張 強,*,孔冬冬,王 月,劉劍宇

    1 中山大學水資源與環(huán)境系,廣州 510275 2 中山大學華南地區(qū)水循環(huán)與水安全廣東省普通高校重點實驗室,廣州 510275 3 中山大學廣東省城市化與地理環(huán)境空間模擬重點實驗室,廣州 510275

    ?

    基于多水文改變指標評價東江流域河流流態(tài)變化及其對生物多樣性的影響

    顧西輝1,2,3,張 強1,2,3,*,孔冬冬1,2,3,王 月1,2,3,劉劍宇1,2,3

    1 中山大學水資源與環(huán)境系,廣州 510275 2 中山大學華南地區(qū)水循環(huán)與水安全廣東省普通高校重點實驗室,廣州 510275 3 中山大學廣東省城市化與地理環(huán)境空間模擬重點實驗室,廣州 510275

    新豐江、楓樹壩和白盆珠3座大型水庫的建立對東江流域河道流量和河流流態(tài)過程有了較大改變,威脅河道下游生態(tài)系統(tǒng)的健康?;趶V義指標生態(tài)剩余和生態(tài)赤字評價了東江流域受水庫影響后流域生態(tài)需水需求目標總的盈余和缺失變化過程,基于IHA32指標計算的Do和DHRAM評價了水庫對下游河段河流水文過程總的改變程度以及威脅河道生態(tài)系統(tǒng)健康的風險性大小,并進一步分析了對河道生物多樣性的影響。研究結果如下:(1)水庫對流量歷時曲線(FDC)有顯著影響,曲線上部下降,尾部上升,尤其體現(xiàn)在秋季和冬季。降水對年與夏季生態(tài)剩余影響較大,水庫對各季節(jié)生態(tài)剩余和生態(tài)赤字均有較大影響:秋季和冬季生態(tài)赤字幾乎為0,生態(tài)剩余顯著增加。生態(tài)剩余和生態(tài)赤字與大部分IHA中32個指標具有很強的相關性,可作為衡量東江流域年和季節(jié)徑流變化的生態(tài)指標。(2)龍川、河源、嶺下和博羅4站點總體改變程度分別為58.48%、54.04%、54.32%和52.47%。河流流態(tài)變化導致總季節(jié)生態(tài)剩余增加并維持在較高水平,進一步引起河流生物多樣性下降,并維持在較低水平。龍川和河源兩站河流流態(tài)的變化對河流生態(tài)系統(tǒng)造成了高風險性影響,嶺下和博羅兩站則為中等風險。

    生態(tài)剩余;生態(tài)赤字;IHA;水文改變;東江流域

    河流流態(tài)和河道流量因跟水資源緊密相關,所以是水循環(huán)過程中最受關注的元素,并且是保持河道生態(tài)系統(tǒng)完整性的核心要素,因此研究河流流態(tài)變化特征對于流域水資源管理和生態(tài)系統(tǒng)保護具有重要意義。氣候變化和人類活動導致水資源時空分布特征[1]及河流流態(tài)[2]發(fā)生變異,通過水庫或大壩導致下游河道徑流失去了原有的天然變化特征[3],并進一步影響河道生態(tài)系統(tǒng)。因此在進行水資源開發(fā)與利用時,有必要考慮徑流變化的生態(tài)效應。為了評價河流流態(tài)變化過程和對生態(tài)系統(tǒng)的影響,常常需要指標來量化水庫或大壩對河流水文改變程度。Olden和Poff總結了170多個水文指標,發(fā)現(xiàn)這些指標過于冗余且具有明顯相關性[4]。Richter等總結并歸納了33個水文改變指標(IHA,Indicators of Hydrologic Alteration)來充分反映河流流態(tài)年內及年際變化特征[5],并被廣泛使用[6- 7]。盡管相比170多個水文指標,IHA中的33個指標已有較大簡化,但指標之間相關性問題仍然沒有得到較好解決[8],且不利于水庫等基于生態(tài)河流流態(tài)的運行調度[9]。為解決上述問題,本文采用基于流量歷時曲線(Flow Duration Curve, FDC)的生態(tài)赤字、生態(tài)剩余等[10]無維度指標評價河道流量需求總的缺失和盈余,并進一步利用基于IHA32指標計算的Do(Degree of hydrologic alteration)[11]和DHRAM(Dundee Hydrological Regime Assessment Method)[12]計算水庫對河流流態(tài)總體改變程度以及造成的河道生態(tài)系統(tǒng)的風險性大小。生態(tài)赤字、生態(tài)剩余、Do和DHRAM等可以從總體角度分析河流流態(tài)變化過程,利于流域水資源管理。

    東江流域內建有新豐江、楓樹壩及白盆珠三大水庫,且長期為廣州、深圳等珠三角重要城市及香港供水。近幾十年來,由于氣候變化和人類活動影響,尤其是流域內三座控制型水庫的影響,河道流量過程已有較大改變。以往的研究主要集中在河道流量年或者季節(jié)性趨勢變化,較少有研究河流流態(tài)改變程度以及對河道生態(tài)系統(tǒng)的影響,尤其是3座大型水庫對下游水文過程改變引起相關生態(tài)的風險性評價。因此本文研究目標如下:(1)利用生態(tài)赤字和生態(tài)剩余兩個無維度指標分析水庫影響下河道生態(tài)流量需求總的缺失和盈余變化過程;(2)基于IHA32指標計算Do和DHRAM評價水庫對下游水文過程總的改變程度以及造成的河道生態(tài)系統(tǒng)風險性大??;(3)進一步分析水庫引起流域水文改變后,對河道生物多樣性的影響。

    1 研究區(qū)域與數(shù)據(jù)

    東江流域是珠江流域的重要支流,發(fā)源于江西省,主要河段位于廣東省,全長562 km,集水面積為35340 km2。東江流域為廣東省主要城市及香港等主要水源區(qū),水資源得到高強度開發(fā)。另外,流域內建有3座控制性大型水庫[13],顯著改變了東江流域河流流態(tài)、增加了水文過程復雜度[13]。

    本文選取東江流域主要控制性水文站點(龍川、河源、嶺下和博羅4站)[13]1954—2009年日流量數(shù)據(jù)進行分析,數(shù)據(jù)來源于廣東省水文局,無缺測,并且經(jīng)過系統(tǒng)整編,質量可靠。水文序列按照水庫建成時間,劃分為受水庫影響前序列(天然序列)和受水庫影響后序列(改變序列)。龍川站和河源站分別位于楓樹壩和新豐江水庫下游[13],因此分別以兩個水庫建成時間1974、1962年作為分割點;嶺下站和博羅站均位于楓樹壩和新豐江水庫下游[13],張強等發(fā)現(xiàn)嶺下和博羅兩站徑流在1973年有一個改變[14],因此以楓樹壩水庫建成時間1974年作為分割點。另外選取分布在整個東江流域29個降水站點1959—2009年日降水數(shù)據(jù),數(shù)據(jù)有少量缺測。

    2 研究方法

    2.1 生態(tài)剩余和生態(tài)赤字

    Vogel等人于2007年提出生態(tài)剩余和生態(tài)赤字兩個廣義指標來評價河道徑流生態(tài)機制[10]。生態(tài)剩余和生態(tài)赤字均以流量歷時曲線(FDC)為基礎。FDC由選擇時間段的日流量數(shù)據(jù)構造,衡量流量超過一個給定值的時間歷時百分比。一段時間內日流量數(shù)據(jù)Qi由大到小排列,其超過概率為[10]:

    pi=i/(n+1)

    (1)

    式中,i為秩次,n為日流量觀測值Qi的樣本大小。FDC可以表述為Qi是pi的函數(shù)。日流量序列既可以構造年尺度FDC也可以構造季節(jié)尺度FDC。東江流域1954—2009年日流量數(shù)據(jù),以水庫建成時間為分割點,分割點前的序列代表天然機制徑流,分別構造分割點前序列每一年年FDC和季節(jié)FDC,然后求得25%和75%分位數(shù)的年FDC和季節(jié)FDC,作為河流生態(tài)系統(tǒng)保護范圍。一個給定年的年FDC或者季節(jié)FDC高于75%分位數(shù)FDC,兩個曲線圍成的面積定義為生態(tài)剩余;一個給定年的年FDC或者季節(jié)FDC低于25%分位數(shù)FDC,兩個曲線圍成的面積定義為生態(tài)赤字[10]。生態(tài)剩余和生態(tài)赤字統(tǒng)一定義為生態(tài)徑流(eco-flow)指標。

    2.2 水文改變程度評價

    IHA32個指標用于量化河流流態(tài)變化特征[15]。分割前序列每個參數(shù)的范圍被分為3類,將默認分類分位數(shù)由33%、67%調整為25%和75%,從而和FDC分類標準一致。用主成分分析從IHA32個指標中識別與生態(tài)最相關的水文指標,定量評價改變后河流流態(tài)與天然河流流態(tài)的偏差,水文改變程度計算如下[11]:

    (2)

    式中,Di為第i個指標的水文改變程度,No,i為改變后徑流序列IHA值在25%—75%分類范圍內的年數(shù),Ne相應期望年數(shù)(Ne=P×NT,P為50%,NT改變后徑流序列的總年數(shù))基于每個指標的改變程度Di,改變后序列總的水文改變程度計算如下[11]:

    (3)

    Black等人于2005年提出了另一個廣義的水文改變指標DHRAM來衡量人類活動對水文機制改變的程度和范圍[12]?;贗HA32個指標改變前后變化百分比,DHRAM計算一個分數(shù)值(變化范圍為0—30),分數(shù)值越大,河流流態(tài)改變程度越大,生態(tài)系統(tǒng)遭受破壞的風險越大。Do和DHRAM共同作為衡量總的水文改變程度的指標,評價水文機制改變對河流生態(tài)系統(tǒng)的風險。

    2.3 生物多樣性影響評價

    SI(Shannon Index,Shannon指數(shù))是運用最廣泛的評價生物多樣性指標[16]:

    (4)

    式中,pi表示群落屬于第i個物種的比例。SI越大表示生物多樣性越豐富。Yang等人用GP(Genetic Programming)算法,基于IHA32個指標建立了SI與水文指標的最佳擬合關系[17]:

    (5)

    式中,Dmin表示最小一天流量的“Julian”日期,Min3、Min7分別表示為最小3d流量和最小7天流量,Max3表示為最大3d流量,Q3、Q5分別表示為3、5月流量,Rrate表示為連續(xù)日流量之間正差異均值。由于缺乏東江流域河道生物群落和種類數(shù)量數(shù)據(jù),無法直接計算SI指標,因此公式(5)SI與水文指標構建的關系式,則可以初步粗略預測水庫建立后對河道生物多樣性的影響程度,但是還需要在以后的生物群落和種類數(shù)據(jù)實際測量中進行進一步驗證。

    3 研究結果

    3.1 生態(tài)徑流指標變化

    圖1為建庫前后年和季節(jié)歷年FDC散點分布特征。由圖1看出,年和夏季FDC建庫前后FDC分布范圍較一致,建庫后高流量和低流量能較好覆蓋建庫前高流量和低流量出現(xiàn)區(qū)域。對于春、秋以及冬季,建庫前后高流量和低流量出現(xiàn)范圍具有較大差別,尤其以秋、冬季最為明顯。秋、冬季建庫后高流量量級和次數(shù)明顯降低,低流量量級和次數(shù)明顯增加,這必然導致秋、冬季由低流量引起的生態(tài)剩余增加,生態(tài)赤字減小。由FDC的變化可以初步判斷出水庫對年和季節(jié)生態(tài)指標影響,然而生態(tài)指標的變化特征還受降水的影響(圖2,圖3)。

    圖2為通過年FDC和季節(jié)FDC得到年和季節(jié)生態(tài)指標(生態(tài)剩余和生態(tài)赤字)及1959—2009年間年和季節(jié)降水距平的時間變化特征。從年尺度來看,年生態(tài)剩余的時間變化與降水距平較為一致,受水庫影響較少。降水增加導致年高流量增加,年FDC超過75%分位數(shù)線的部分增加,生態(tài)剩余上升,東江流域河道生態(tài)需水處于過于滿足狀態(tài)。整個東江流域20世紀70年代中期到80年代中期,降水明顯增加,生態(tài)剩余較為明顯,20世紀80年代中期之后,降水下降,生態(tài)剩余也呈下降趨勢。相比年生態(tài)剩余,年生態(tài)赤字與降水距平并不完全一致,水庫對年生態(tài)剩余的影響較大。降水減少,導致年低流量減少,年FDC低于25%分位數(shù)線的部分增加,生態(tài)赤字上升,東江流域河道需水處于過于缺乏狀態(tài)。年生態(tài)赤字與降水距平較一致的時間段集中在1959—1965、1990—2009年兩個時間段。這兩個時間段,降水明顯減少,生態(tài)赤字同時明顯增加。但是對于龍川、嶺下和博羅站,1974—1985年生態(tài)赤字的變化并未完全體現(xiàn)出降水的減少程度,這段時間由于楓樹壩和新豐江兩座大型水庫的調蓄作用,年低流量的減少程度減緩了。

    圖2 生態(tài)徑流指標與降水距平百分率的時間變化Fig.2 Changes in the annual and season eco-flow metrics and the annual precipitation anomaly in the period from 1954 to 2009紅色和綠色條形為降水距平百分率,藍色曲線為生態(tài)剩余,黑色曲線為生態(tài)赤字

    相比年生態(tài)指標變化,季節(jié)生態(tài)指標時間變化差異性較為明顯,并且變化幅度較大,穩(wěn)定性較小。季節(jié)生態(tài)剩余和生態(tài)赤字與降水距平的時間變化一致性較小,水庫的影響更加明顯(圖2)。在4個季節(jié)中,夏季生態(tài)指標與降水距平具有更好的一致性。但是對于河源站,1962—1970年夏季降水明顯增加,生態(tài)剩余卻呈顯著下降,新豐江水庫的建成對河源站高流量有顯著的影響。相比春季和夏季,枯水期秋季和冬季,東江流域生態(tài)指標受到水庫的顯著影響,降水的影響受到極大的削弱。河源、嶺下和博羅站在1962年之后,不管降水如何變化,生態(tài)赤字一直穩(wěn)定在0附近,新豐江水庫對保持秋、冬季低流量在較高的水平具有重要的作用;由于降水的減少,加上楓樹壩和新豐江兩座水庫的調節(jié)作用,秋季生態(tài)剩余在1974年之后呈穩(wěn)定下降趨勢,龍川和河源兩站冬季生態(tài)剩余則在高位震蕩。

    圖3 年和季節(jié)生態(tài)指標10a年際間變化箱型圖Fig.3 Boxplot of decadal changes in the seasonal ecosurplus and ecodeficit in the period from 1954 to 2009

    圖3展示了年和季節(jié)生態(tài)指標10a統(tǒng)計變化特征。4個站點年生態(tài)剩余均表現(xiàn)為1974年之前處于較低水平,1974年之后突然增加;年生態(tài)赤字變化幅度相比年生態(tài)剩余較小,龍川和河源兩站年生態(tài)赤字呈緩慢增加趨勢,嶺下和河源兩站1964—1993年維持在0值附近,1994—2009年顯著增加。4個站點春季生態(tài)剩余10a統(tǒng)計中值呈現(xiàn)規(guī)律的拋物線形:1954—1983年呈增加趨勢,1983—2009年呈下降趨勢;河源站春季生態(tài)赤字1954—1983呈現(xiàn)“低-高-低”的循環(huán),然后在最近的20a呈增加趨勢,表明河源站豐水期徑流有一個10a的周期,在20世紀80年代之后呈增加趨勢;龍川、嶺下和博羅在20世紀90年代之前春季生態(tài)赤字一直維持在穩(wěn)定的低水平,然后在近20a呈突然增加趨勢。

    3.2 水文整體改變程度的評價

    Do和DHRAM被用來衡量水文機制整體改變程度(表1)。DHRAM基于IHA5組33個參數(shù),將每類參數(shù)的平均均值和離差系數(shù)的變化程度分為3類(1代表改變程度最低,3代表改變程度最高),然后求出每類總的改變程度,進行總體評價,確定水文改變對河流生態(tài)系統(tǒng)造成的風險等級。從表1中可以看出,龍川站IHA1、3、5組(分別對應月徑流、極端徑流出現(xiàn)時間和徑流改變頻率)整體改變程度較大,達到中等及以上水平;河源站和嶺下站均集中在IHA1、3組整體改變程度較大,博羅站僅IHA第3組整體改變程度較大,達到中等水平。DHRAM將總體改變程度分為5級,1級代表水文機制沒有受到影響(生態(tài)系統(tǒng)無風險),5級代表受到嚴重影響(生態(tài)系統(tǒng)有嚴重風險)。從DHRAM總體打分及評價結果來看,龍川和河源兩站總得分均為13分,改變等級為4級——高風險影響;嶺下和博羅兩站總得分均為9分,改變等級為3級——中等風險影響。從Do來看,龍川站水文機制改變程度最大,達到58.48%,博羅站最小達到52.47%。與DHRAM不同的是,河源站和嶺下站水文機制改變程度幾乎相等,均在54%左右。盡管4站點之間Do值差別較小,但是均在50%以上,說明東江流域水文機制均較建庫前有較大改變。采用Do值的另一個作用是可以用于水庫調度中下游水文機制改變程度最小的調度目標,這在接下來的研究中繼續(xù)進行。

    表1 水庫建成前后水文改變的總體評價:DHRAM和Do結果

    3.3 生態(tài)徑流的變化對生物多樣性的影響

    已經(jīng)有170多個水文指標[4]公開發(fā)表,盡管IHA總結了32個水文指標[5],但是與生態(tài)最相關指標(Ecologically Relevant Hydrologic Indicators, ERHIs)的選擇依然是關鍵問題。本文用PCA從IHA32個指標中選擇ERHIs做進一步的分析。選擇特征值大于1的主成分,保證結果中每個變量至少有更多的方差,共選擇了7個主成分:PC1—PC7。因為在一個主成分下選擇多個指標,代表河流流態(tài)相同的方程成分,所以在選擇的ERHIs中,很多指標有較高的相關性例如:最小7d、最小30d和最小90d等。因此在明顯具有相關性的指標中,只選擇一個最高荷載的指標:最小7d流量。在4個站點中出現(xiàn)次數(shù)最多的指標被選為東江流域ERHIs:最小7d流量(Min7)、最大7天流量(Max7)、最小流量出現(xiàn)日期(Dmin)、最大流量出現(xiàn)日期(Dmax),高流量平均歷時(High Pulse Count, HPC)。上述ERHIs與Yang等人用GP篩選出的6個ERHIs基本吻合[17],兩者相差水文轉變數(shù)量(Number of reversals,RV)、上升率(Rise rates,Rrate)兩個指標。為了全面考慮,本文同樣將RV和Rrate兩個指標納入ERHIs中分析(圖4)。

    從圖4中可以看出,最小7d流量在1974年之后突然增加,然后維持在高位震蕩。最大7d流量由1974年之前的高位震蕩,突然減小,并且穩(wěn)定在較低水平。最小日流量出現(xiàn)時間建庫后序列震蕩幅度明顯增加,以嶺下站和河源站最顯著。最大日流量出現(xiàn)時間建庫后序列較建庫前更加穩(wěn)定,震蕩幅度顯著減小。高流量天數(shù)(HPC)建庫后均成緩慢震蕩下降趨勢。高流量量級和天數(shù)減少,直接減弱洪水對河道的沖於作用,影響底層無脊椎生物生存與繁衍;河岸邊濕地生態(tài)用水補充可能困難,濕地地下水也較難補充,并且影響濕地土壤營養(yǎng)補充,進而威脅濕地系統(tǒng)的生物多樣性。Rrate建庫后的變化很顯著,均表現(xiàn)為突然下降,然后穩(wěn)定在較低水平。RV則相反,一直呈顯著增加趨勢。ERHIs的變化特征必將引起總生態(tài)赤字和總生態(tài)剩余的變化(圖5)并且導致河道水生生態(tài)系統(tǒng)生物多樣性的變化(SI,圖6)。

    圖4 ERHIs隨時間變化特征Fig.4 The time-varying characteristics of ERHIs

    圖5 總季節(jié)生態(tài)剩余和總季節(jié)生態(tài)赤字的時間變化特征Fig.5 The time-varying characteristics of total season ecosurplus and ecodeficit圖中陰影部分表示loess擬合曲線的95%置信區(qū)間

    總季節(jié)生態(tài)剩余和總季節(jié)生態(tài)赤字就是分別將各季節(jié)生態(tài)剩余和生態(tài)赤字求和(圖5)??偧竟?jié)生態(tài)剩余和生態(tài)赤字可以反映出季節(jié)尺度生態(tài)變化在年尺度的影響,從而和年尺度的生態(tài)多樣性指標SI值對應??偧竟?jié)生態(tài)剩余和生態(tài)赤字用局部加權多項式(Loess函數(shù))回歸進行擬合,并給出了95%置信區(qū)間,以判斷其變化趨勢(圖5)。從圖5中可以看出,4個水文站點總季節(jié)生態(tài)剩余明顯高于總季節(jié)生態(tài)赤字,并且總季節(jié)生態(tài)剩余和生態(tài)赤字均具有相似的變化特征。1980年之前總季節(jié)生態(tài)剩余一直呈增加趨勢,1980年之后處于穩(wěn)定狀態(tài)。龍川站總季節(jié)生態(tài)赤字一直處于平穩(wěn)狀態(tài),變化較小,河源、嶺下和博羅三站總季節(jié)生態(tài)赤字1974年之前處于微弱的下降狀態(tài),1974年之后反而處于緩慢的上升狀態(tài),其中夏季生態(tài)赤字的增加起到了關鍵的作用(圖2)。

    圖6給出了生物多樣性指標SI的變化特征。從圖6中可以看出,龍川、嶺下和博羅3站SI值在1974年有一個明顯的轉折點:1974年之前持續(xù)下降,1974年之后緩慢上升,1990年代之后在緩慢下降。相比建庫前,龍川、嶺下和博羅三站SI值均在下降,生物多樣性減小,水庫對河流生態(tài)系統(tǒng)造成了顯著影響。河源站SI值在水庫建成后1990年代呈顯著下降趨勢, 生物多樣性受到嚴重影響,這與河源站1990年代降水顯著下降有顯著關系(圖2)。比較圖5和圖6,建庫后總季節(jié)生態(tài)剩余和生態(tài)赤字的處于穩(wěn)定性變化,SI值也重新處于穩(wěn)定狀態(tài),說明:建庫后河流生態(tài)河流流態(tài)的重新平衡,對建庫后的生物多樣性的重新處于平穩(wěn)狀態(tài)具有重要作用。

    圖6 生物多樣性指標SI時間變化特征Fig.6 The time-varying characteristics of biodiversity indicators SI黑色曲線表示loess函數(shù)擬合曲線,陰影部分表示loess擬合曲線的95%置信區(qū)間

    3.4 生態(tài)徑流指標與IHA32指標比較

    為了理解生態(tài)徑流指標和IHA32指標之間的關系,IHA32指標與各生態(tài)徑流指標的相關關系如圖7所示。大部分生態(tài)徑流指標與IHA32指標具有很強的正或者負的相關關系(圖7)。冬季月徑流(12—2月)與冬季、年以、季節(jié)總生態(tài)剩余、總生態(tài)改變以及SI具有很強的正相關關系。最小1、3、7、30d和90d流量幾乎與各生態(tài)指標(除去夏季和秋季生態(tài)剩余和生態(tài)赤字)均有較強的正相關關系。相比不同歷時的最小流量,最大1、3、7、30d和90d流量與各生態(tài)徑流指標相關關系較弱。FR(Fall rate)與LPD(Low pulse duration)與大部分生態(tài)徑流指標具有較強的負相關關系。在所有的生態(tài)徑流指標中,冬季、年以及季節(jié)總生態(tài)剩余指標、秋季和冬季生態(tài)赤字、總生態(tài)改變以及SI與IHA32指標具有較好的相關關系,能夠較好的表達IHA32指標的信息:例如冬季生態(tài)剩余增加,生態(tài)赤字減少,可以充分反映出冬季月徑流增加,不用歷時最小流量的增加等。在所有的生態(tài)指標中,夏季生態(tài)剩余和生態(tài)赤字與IHA32指標的相關關系最弱,表達IHA32指標的信息最少。通過以上的分析,生態(tài)徑流指標可以很好的反映出IHA32指標的信息,生態(tài)剩余和生態(tài)赤字可以提供水文改變的一個好的評價標準。生態(tài)剩余和生態(tài)赤字與IHA32指標的計算是相互獨立的,同時生態(tài)剩余和生態(tài)赤字的應用能夠有效的解決大量水文指標之間的冗余與相關關系,可以作為衡量水文改變的很好的指標[8,18]。

    4 結論

    用東江流域1954—2009年日流量數(shù)據(jù)計算生態(tài)剩余和生態(tài)赤字,并結合Do和DHRAM評價東江流域水文整體改變程度,分析了ERHIs、總季節(jié)生態(tài)剩余與生態(tài)赤字以及SI指標的變化,研究水文改變對生態(tài)多樣性的影響。研究結果如下:

    (1)由于水庫的影響,建庫后高流量量級和次數(shù)明顯降低,低流量量級和次數(shù)明顯增加,導致FDC上端向下移動,尾部向上移動。FDC的變化引起高流量部分低于25%FDC,產(chǎn)生生態(tài)赤字,低流量部分高于75%FDC,產(chǎn)生生態(tài)剩余。

    (2)年生態(tài)剩余的時間變化與降水距平較為一致,受水庫影響較少,年生態(tài)赤字則受水庫影響較大。夏季生態(tài)剩余和生態(tài)赤字與降水的變化較為一致,其他季節(jié)則受水庫的影響較大,尤其是秋季和冬季,建庫后生態(tài)赤字幾乎為0,生態(tài)剩余顯著增加。

    (3)龍川和河源兩站河流流態(tài)的變化對河流生態(tài)系統(tǒng)造成了高風險的影響,嶺下和博羅兩站則是中等風險(DHRAM評價),龍川、河源、嶺下和博羅4站點總體改變程度分別為58.48%、54.04%、54.32%和52.47%(Do)。河流流態(tài)變化導致總季節(jié)生態(tài)剩余增加并維持在較高水平,進一步引起河流生物多樣性下降,并維持在較低水平。

    (4)生態(tài)徑流指標(生態(tài)剩余和生態(tài)赤字)與IHA32指標具有較好的相關關系,能夠體現(xiàn)大部分IHA參數(shù)的信息。生態(tài)剩余和生態(tài)赤字能夠很好的應用與東江流域反映東江流域年和季節(jié)徑流變化。

    [1] 張強, 孫鵬, 陳喜, 陳曉宏. 1956—2000年中國地表水資源狀況: 變化特征、成因及影響. 地理科學, 2011, 31(12): 1430- 1436.

    [2] Barnett T P, Pierce D W, Hidalgo H G, Bonfils C, Santer B D, Das T, Bala G, Wood A W, Nozawa T, Mirin A A, Cayan D R, Dettinger M D. Human-Induced changes in the hydrology of the western united states. Science, 2008, 319(5866): 1080- 1083.

    [3] Poff N L, Olden J D, Merritt D M, Pepin D M. Homogenization of regional river dynamics by dams and global biodiversity implications. The National Academy of Sciences of the USA, 2007, 104(14): 5732- 5737.

    [4] Olden J D, Poff N L. Redundancy and the choice of hydrologic indices for characterizing streamflow regimes. River Research and Applications, 2003, 19(2): 101- 121.

    [5] Richter B D, Baumgartner J V, Powell J, Braun D P. A method for assessing hydrologic alteration within ecosystems. Conservation Biology, 1996, 10(4): 1163- 1174.

    [6] Zhang Q, Xiao M H, Liu C L, Singh V P. Reservoir-induced hydrological alterations and environmental flow variation in the East River, the Pearl River basin, China. Stochastic Environmental Research and Risk Assessment, 2014, 28(8): 2119- 2131.

    [7] Chen Y Q D, Yang T, Xu C Y, Zhang Q, Chen X, Hao Z C. Hydrologic alteration along the Middle and Upper East River (Dongjiang) basin, South China: a visually enhanced mining on the results of RVA method. Stochastic Environmental Research and Risk Assessment, 2010, 24(1): 9- 18.

    [8] Gao Y X, Vogel R M, Kroll C N, Poff N L, Olden J D. Development of representative indicators of hydrologic alteration. Journal of Hydrology, 2009, 374(1/2): 136- 147.

    [9] Shiau J T, Wu F C. Compromise programming methodology for determining instream flow under multiobjective water allocation criteria. Journal of the American Water Resources Association, 2006, 42(5): 1179- 1191.

    [10] Vogel R M, Sieber J, Archfield S A, Smith M P, Apse C D, Lee A H. Relations among storage, yield, and instream flow. Water Resources Research, 2007, 43(5): W05403.

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

    [12] Black A R, Rowan J S, Duck R W, Bragg O M, Clelland B E. DHRAM: a method for classifying river flow regime alterations for the EC Water Framework Directive. Aquatic Conservation: Marine and Freshwater Ecosystems, 2005, 15(5): 427- 446.

    [13] Zhou Y, Zhang Q, Li K, Chen X H. Hydrological effects of water reservoirs on hydrological processes in the East River (China) basin: complexity evaluations based on the multi-scale entropy analysis. Hydrological Processes, 2012, 26(21): 3253- 3262.

    [14] 張強, 崔瑛, 陳永勤. 水文變異條件下的東江流域生態(tài)徑流研究. 自然資源學報, 2012, 27(5): 790- 800.

    [15] The Nature Conservancy. Indicators of hydrologic alteration version 7.1 user′s manual: 2009. (2009-04). http://www.conservationgateway. org/ConservationPractices/Freshwater/EnvironmentalFlows/MethodsandTools/IndicatorsofHydrologicAlteration/Pages/IHA-Software-Download.aspx.

    [16] Kuo S R, Lin H J, Shao K T. Seasonal changes in abundance and composition of the fish assemblage in Chiku Lagoon, southwestern Taiwan. Bulletin of Marine Science, 2001, 68(1): 85- 99.

    [17] Yang Y C E, Cai X M, Herricks E E. Identification of hydrologic indicators related to fish diversity and abundance: A data mining approach for fish community analysis. Water Resources Research, 2008, 44(4): W04412.

    [18] Gao B, Yang D W, Zhao T G, Yang H B. Changes in the eco-flow metrics of the Upper Yangtze River from 1961 to 2008. Journal of Hydrology, 2012, 448- 449: 30- 38.

    Based on multiple hydrological alteration indicators evaluating the characteristics of flow regime with the impact on the diversity of hydrophily biology

    GU Xihui1,2,3, ZHANG Qiang1,2,3,*, KONG Dongdong1,2,3, WANG Yue1,2,3, LIU Jianyu1,2,3

    1DepartmentofWaterResourcesandEnvironment,SunYat-senUniversity,Guangzhou510275,China2KeyLaboratoryofWaterCycleandWaterSecurityinSouthernChinaofGuangdongHighEducationInstitute,SunYat-senUniversity,Guangzhou510275,China3GuangdongProvincialKeyLaboratoryofUrbanizationandGeo-simulation,SchoolofGeographyandPlanning,SunYat-senUniversity,Guangzhou510275,China

    Dam-induced hydrological alterations and related ecological problems have created considerable concern for hydrologists, ecologists, and policy-makers. The East River basin in China is the primary provider of water resources for mega-cities within the Pearl River Delta and meets 80% of the annual water demand of Hong Kong. Daily streamflow data from 4 hydrological stations covering the period of 1954—2009 are analyzed with respect to eco-flow regimes of the East River basin using eco-flow metrics (ecosurplus and ecodeficit) and the Indicators of Hydrologic Alteration (IHA) technique. In this study, ecodeficit and ecosurplus were analyzed to determine the ecological impact of water impoundments. In addition, Doand DHRAM (Dundee Hydrological Regime Assessment Method) were employed to evaluate the degree of alteration of hydrological regimes, and ERHIs (Ecologically Relevant Hydrologic Indicators) were used to analyze the influence of hydrological alterations on ecological diversity. Results indicate that the flow duration curve (FDC) is significantly affected by water reservoirs, which is largely mirrored by the running down of the upper curve tail and running up of the lower curve tail, especially in autumn and winter. In this way, the magnitude and frequency of high flows decrease and those of low flows increase because of reservoir regulation. Precipitation exerts greater influence on annual and summer ecosurplus, whereas seasonal ecosurplus and ecodeficit are more influenced by the hydrological regulations of the reservoirs. The ecodeficit in the autumn and winter is nearly zero, whereas the ecosurplus in these same seasons increases significantly. In addition, ecodeficit and ecosurplus in other seasons, particularly autumn and winter, are more influenced by reservoir regulation. Impacts of reservoirs on hydrological regimes and eco-flow regimes differ from one station to another owing to varying degrees of influence from the reservoirs on hydrological processes. The longer the distance between a reservoir and hydrological station, the weaker the influence the water reservoir has on hydrological processes. Because of changes in the streamflow mechanism, the fluvial ecosystems are under a high risk of degradation at the Longchuan and Heyuan stations and a moderate risk of degradation at the Lingxia and Boluo stations. The general alteration degrees of hydrology are 58.48%, 54.04%, 54.32%, and 52.47% at the Longchuan, Heyuan, Lingxia, and Boluo stations, respectively. The changes in the streamflow mechanism tend to increase seasonal ecosurplus and trigger risk of a high level of ecosurplus that causes a further decline in biodiversity. Comparing the eco-flow and IHA metrics demonstrated that combining these two groups has the potential to provide a sufficient measure of the change in the flow regime. Ecodeficit and ecosurplus can be accepted in the evaluation of alterations of hydrological processes at annual and seasonal time scales. This study provides a theoretical background for water resource management with consideration of eco-flow variations in response to reservoir regulation in other highly regulated river basins of the globe.

    ecosurplus; ecodeficit; IHA; hydrology alteration; East River

    國家杰出青年科學基金項目(51425903);國家自然科學基金重大國際合作研究項目(51320105010);香港特別行政區(qū)研究資助局項目(CUHK441313);中山大學濱海小流域自然地理綜合過程觀測與實驗平臺建設(2012年度);中央高校基本科研業(yè)務費專項資金;廣東省教育部產(chǎn)學研結合項目(2012B091100471)

    2014- 12- 12;

    日期:2016- 01- 15

    10.5846/stxb201412122480

    *通訊作者Corresponding author.E-mail: zhangq68@mail.sysu.edu.cn

    顧西輝,張強,孔冬冬,王月,劉劍宇.基于多水文改變指標評價東江流域河流流態(tài)變化及其對生物多樣性的影響.生態(tài)學報,2016,36(19):6079- 6090.

    Gu X H, Zhang Q, Kong D D, Wang Y, Liu J Y.Based on multiple hydrological alteration indicators evaluating the characteristics of flow regime with the impact on the diversity of hydrophily biology.Acta Ecologica Sinica,2016,36(19):6079- 6090.

    猜你喜歡
    東江流態(tài)赤字
    側邊機組故障對泵站前池流態(tài)的影響
    大電機技術(2022年1期)2022-03-16 06:40:24
    改進邊界條件的非恒定流模型在城市河流橡膠壩流態(tài)模擬中的應用
    萬綠水 東江請原創(chuàng)歌曲|全國征集
    嶺南音樂(2016年4期)2016-11-21 09:25:24
    東江本地早快速投產(chǎn)配套技術研究
    浙江柑橘(2016年4期)2016-03-11 20:12:58
    動態(tài)流態(tài)冰蓄冷系統(tǒng)在千級凈化廠房的應用
    機電信息(2015年3期)2015-02-27 15:54:46
    湖南東江水庫移民農(nóng)業(yè)事業(yè)單位體制改革探討
    基于TM遙感影像的河口流態(tài)信息半定量化研究
    互信赤字
    東江縱隊
    軍事歷史(1984年2期)1984-08-21 06:27:02
    一个人看视频在线观看www免费| 国产老妇伦熟女老妇高清| 国产精品久久视频播放| 色视频www国产| 中文字幕人妻熟人妻熟丝袜美| 欧美+日韩+精品| 国产精品一区二区性色av| 最新中文字幕久久久久| 天美传媒精品一区二区| 欧美色欧美亚洲另类二区| 99精品在免费线老司机午夜| 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 18禁在线播放成人免费| 老司机影院成人| 一本精品99久久精品77| 亚洲经典国产精华液单| 26uuu在线亚洲综合色| 狂野欧美激情性xxxx在线观看| 日韩亚洲欧美综合| 国产 一区 欧美 日韩| 国产白丝娇喘喷水9色精品| 欧美日韩在线观看h| 噜噜噜噜噜久久久久久91| 能在线免费观看的黄片| 国产色婷婷99| 国产极品精品免费视频能看的| www日本黄色视频网| 免费观看人在逋| 91在线精品国自产拍蜜月| 插阴视频在线观看视频| 人妻少妇偷人精品九色| 欧美3d第一页| 欧美在线一区亚洲| 日日撸夜夜添| 中文精品一卡2卡3卡4更新| 精品久久国产蜜桃| 99久久中文字幕三级久久日本| 有码 亚洲区| 观看美女的网站| 亚洲av一区综合| 免费人成视频x8x8入口观看| 一本久久中文字幕| 久久精品国产亚洲网站| 99热这里只有精品一区| 91精品一卡2卡3卡4卡| 午夜福利在线观看免费完整高清在 | 青春草视频在线免费观看| 麻豆成人av视频| 成人毛片60女人毛片免费| 亚洲真实伦在线观看| 精品一区二区免费观看| 最近视频中文字幕2019在线8| 好男人视频免费观看在线| 韩国av在线不卡| 亚洲图色成人| 在线播放国产精品三级| av在线老鸭窝| 午夜视频国产福利| 别揉我奶头 嗯啊视频| 亚洲丝袜综合中文字幕| 国产精品一区二区在线观看99 | 乱人视频在线观看| 性插视频无遮挡在线免费观看| 人妻系列 视频| 晚上一个人看的免费电影| 国产成人精品一,二区 | 久久久久久久午夜电影| 最近中文字幕高清免费大全6| 国产老妇伦熟女老妇高清| 高清毛片免费看| 日日摸夜夜添夜夜添av毛片| 国产精品人妻久久久久久| 免费av不卡在线播放| 国产高潮美女av| 久久精品夜夜夜夜夜久久蜜豆| 国产精品一区二区三区四区久久| 成人无遮挡网站| 99热全是精品| 精品国内亚洲2022精品成人| 人体艺术视频欧美日本| 国产不卡一卡二| 我要看日韩黄色一级片| 亚洲国产色片| 日本一本二区三区精品| 九草在线视频观看| 欧美一区二区精品小视频在线| 亚洲av.av天堂| 两个人视频免费观看高清| 免费在线观看成人毛片| 亚洲不卡免费看| 99久久九九国产精品国产免费| 免费观看a级毛片全部| 国产精品一区二区三区四区久久| 欧美zozozo另类| 三级毛片av免费| 国产 一区 欧美 日韩| 欧美最新免费一区二区三区| 美女大奶头视频| 欧美一区二区国产精品久久精品| 国产精华一区二区三区| 欧美成人免费av一区二区三区| 国产精品不卡视频一区二区| 亚洲欧美精品自产自拍| 精品人妻视频免费看| a级毛片a级免费在线| 亚洲国产精品成人久久小说 | 亚洲欧美日韩无卡精品| 青青草视频在线视频观看| 91久久精品国产一区二区成人| 久99久视频精品免费| 晚上一个人看的免费电影| 亚洲久久久久久中文字幕| 美女被艹到高潮喷水动态| 色综合色国产| 国产爱豆传媒在线观看| 久久婷婷人人爽人人干人人爱| or卡值多少钱| 日本一二三区视频观看| 午夜免费激情av| 欧美xxxx性猛交bbbb| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人免费av一区二区三区| 国产色婷婷99| 日本熟妇午夜| 悠悠久久av| 我的女老师完整版在线观看| 欧美在线一区亚洲| 国产69精品久久久久777片| 麻豆成人av视频| 久久亚洲精品不卡| 男插女下体视频免费在线播放| 国产色婷婷99| 国产av在哪里看| 亚洲av熟女| 午夜久久久久精精品| 麻豆久久精品国产亚洲av| 91久久精品国产一区二区三区| 免费看av在线观看网站| 国产精品一区二区在线观看99 | 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 色视频www国产| 成年版毛片免费区| 国产精品一区二区三区四区久久| 麻豆乱淫一区二区| 看免费成人av毛片| 亚洲av.av天堂| 国产成人aa在线观看| 国产精品99久久久久久久久| 九九久久精品国产亚洲av麻豆| 国产精品久久视频播放| 99热精品在线国产| 亚洲国产精品sss在线观看| av卡一久久| 99热这里只有是精品50| 又粗又爽又猛毛片免费看| 狠狠狠狠99中文字幕| 最近最新中文字幕大全电影3| 边亲边吃奶的免费视频| 天美传媒精品一区二区| 久久人妻av系列| 日韩一本色道免费dvd| 午夜亚洲福利在线播放| 天堂影院成人在线观看| 中文字幕av成人在线电影| 69人妻影院| 中文在线观看免费www的网站| 欧美一区二区精品小视频在线| 中国国产av一级| 最近2019中文字幕mv第一页| 日韩一区二区视频免费看| 全区人妻精品视频| 成人毛片60女人毛片免费| 久久热精品热| 九色成人免费人妻av| 在现免费观看毛片| 国内久久婷婷六月综合欲色啪| 亚洲,欧美,日韩| 91久久精品电影网| 长腿黑丝高跟| 性色avwww在线观看| 欧美丝袜亚洲另类| 国产亚洲av片在线观看秒播厂 | 看免费成人av毛片| 成人亚洲欧美一区二区av| 中文欧美无线码| 亚洲av成人av| 夜夜爽天天搞| 在线免费十八禁| 校园春色视频在线观看| 午夜视频国产福利| 精品熟女少妇av免费看| 91av网一区二区| 男女边吃奶边做爰视频| 日本一二三区视频观看| 免费黄网站久久成人精品| av免费观看日本| 午夜亚洲福利在线播放| 麻豆av噜噜一区二区三区| 亚洲丝袜综合中文字幕| 青春草视频在线免费观看| 18+在线观看网站| 身体一侧抽搐| 我的女老师完整版在线观看| 天天躁日日操中文字幕| 亚洲一区高清亚洲精品| 少妇的逼水好多| 12—13女人毛片做爰片一| 国产精品久久久久久精品电影| 日本爱情动作片www.在线观看| 狂野欧美白嫩少妇大欣赏| 大又大粗又爽又黄少妇毛片口| 人体艺术视频欧美日本| 久久久成人免费电影| 岛国在线免费视频观看| 中文字幕制服av| 九草在线视频观看| 国产成人福利小说| 午夜福利成人在线免费观看| 久久热精品热| 在线观看午夜福利视频| 人人妻人人看人人澡| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清专用| 精品国内亚洲2022精品成人| 日产精品乱码卡一卡2卡三| 中国美女看黄片| 国产在线精品亚洲第一网站| 国内少妇人妻偷人精品xxx网站| 91久久精品国产一区二区三区| 免费在线观看成人毛片| 欧美成人精品欧美一级黄| 男女那种视频在线观看| 我的老师免费观看完整版| 天堂√8在线中文| 乱系列少妇在线播放| 日韩一区二区三区影片| 看免费成人av毛片| 午夜a级毛片| 国产高清有码在线观看视频| 国产乱人偷精品视频| 大香蕉久久网| 日韩欧美三级三区| 男人的好看免费观看在线视频| 成人欧美大片| 在线观看一区二区三区| 国产精品美女特级片免费视频播放器| 看黄色毛片网站| 国产成人精品久久久久久| 国产精品日韩av在线免费观看| 国产精品一二三区在线看| 波野结衣二区三区在线| 国产一区亚洲一区在线观看| 一本一本综合久久| 伦精品一区二区三区| 亚洲欧美成人综合另类久久久 | 一级二级三级毛片免费看| 精品99又大又爽又粗少妇毛片| 色94色欧美一区二区| 国国产精品蜜臀av免费| 亚洲精品中文字幕在线视频| 人妻一区二区av| 免费高清在线观看视频在线观看| 91在线精品国自产拍蜜月| 久久久a久久爽久久v久久| 在线天堂最新版资源| 各种免费的搞黄视频| 久久精品夜色国产| 九草在线视频观看| 少妇的逼好多水| 超碰97精品在线观看| 久久人妻熟女aⅴ| 少妇的逼好多水| 日韩亚洲欧美综合| 免费大片黄手机在线观看| 国产av码专区亚洲av| xxx大片免费视频| 久久人人爽人人片av| 亚洲人成77777在线视频| 搡老乐熟女国产| 国产在线一区二区三区精| 亚洲av国产av综合av卡| 男男h啪啪无遮挡| 考比视频在线观看| 国产国语露脸激情在线看| 熟女人妻精品中文字幕| 五月开心婷婷网| 97在线视频观看| 99久久人妻综合| 亚洲国产精品一区二区三区在线| 午夜日本视频在线| 91精品伊人久久大香线蕉| 精品亚洲成a人片在线观看| 国产免费一级a男人的天堂| 久久久久久久精品精品| 国产亚洲av片在线观看秒播厂| 寂寞人妻少妇视频99o| 成人18禁高潮啪啪吃奶动态图 | 久久久久国产网址| 亚州av有码| 狂野欧美激情性bbbbbb| 国产探花极品一区二区| 极品人妻少妇av视频| 99久国产av精品国产电影| 久久这里有精品视频免费| 亚洲精品视频女| 久久综合国产亚洲精品| 老熟女久久久| 两个人的视频大全免费| 伦理电影免费视频| 中文字幕免费在线视频6| 永久网站在线| 久久青草综合色| 91精品一卡2卡3卡4卡| 色哟哟·www| 欧美精品高潮呻吟av久久| 极品人妻少妇av视频| 亚洲精品国产av蜜桃| 欧美日韩视频精品一区| 久久99热6这里只有精品| av网站免费在线观看视频| 国产亚洲精品久久久com| 亚洲中文av在线| 91在线精品国自产拍蜜月| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黑丝袜美女国产一区| 卡戴珊不雅视频在线播放| 亚洲,欧美,日韩| 国产精品一区二区在线不卡| 国产精品国产av在线观看| 免费久久久久久久精品成人欧美视频 | 国产成人精品无人区| 黄色毛片三级朝国网站| 久久鲁丝午夜福利片| 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 久久久国产一区二区| 欧美一级a爱片免费观看看| 一级毛片 在线播放| 亚洲av在线观看美女高潮| 少妇丰满av| 中文字幕制服av| 91国产中文字幕| 51国产日韩欧美| 中文天堂在线官网| 国产精品一区www在线观看| 我的老师免费观看完整版| 大香蕉97超碰在线| 国产免费福利视频在线观看| 亚洲精品国产av成人精品| 精品国产一区二区三区久久久樱花| 大片免费播放器 马上看| 国产成人精品福利久久| 91在线精品国自产拍蜜月| 一个人看视频在线观看www免费| 青春草亚洲视频在线观看| 国产av一区二区精品久久| 男女边吃奶边做爰视频| av视频免费观看在线观看| √禁漫天堂资源中文www| 一级毛片aaaaaa免费看小| 一级片'在线观看视频| 男人添女人高潮全过程视频| 自拍欧美九色日韩亚洲蝌蚪91| tube8黄色片| 亚洲精品乱久久久久久| 久久国产精品男人的天堂亚洲 | 国产精品嫩草影院av在线观看| 天堂中文最新版在线下载| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 五月天丁香电影| 人妻夜夜爽99麻豆av| 亚洲精品日本国产第一区| 两个人的视频大全免费| 国产精品99久久99久久久不卡 | 久久久久久人妻| 视频区图区小说| 久久精品国产亚洲av天美| 九九久久精品国产亚洲av麻豆| 亚洲av在线观看美女高潮| www.色视频.com| 久久精品国产鲁丝片午夜精品| 色婷婷av一区二区三区视频| 又黄又爽又刺激的免费视频.| 亚洲av男天堂| 午夜福利视频在线观看免费| 各种免费的搞黄视频| 日本欧美视频一区| 少妇的逼好多水| 久久精品国产亚洲av涩爱| 国产高清国产精品国产三级| 在线播放无遮挡| 久久午夜综合久久蜜桃| 99国产精品免费福利视频| 久久久久久久精品精品| 性高湖久久久久久久久免费观看| 99九九在线精品视频| 日本免费在线观看一区| 亚洲欧美中文字幕日韩二区| 亚洲一区二区三区欧美精品| 啦啦啦啦在线视频资源| 日韩伦理黄色片| 国产在线一区二区三区精| 一级毛片 在线播放| 国产成人精品一,二区| 久久人人爽人人爽人人片va| 少妇的逼水好多| 黑人欧美特级aaaaaa片| 黄色一级大片看看| 精品久久久精品久久久| 永久网站在线| 久久久久国产精品人妻一区二区| 国产黄频视频在线观看| 久久久久久久国产电影| 99久久中文字幕三级久久日本| 国产日韩欧美在线精品| 男女国产视频网站| 一个人看视频在线观看www免费| 国精品久久久久久国模美| 亚洲人成77777在线视频| 春色校园在线视频观看| 母亲3免费完整高清在线观看 | 国产精品女同一区二区软件| av有码第一页| 高清欧美精品videossex| 欧美精品一区二区大全| 麻豆精品久久久久久蜜桃| 亚洲天堂av无毛| 久久精品国产亚洲av天美| 成人国语在线视频| 亚洲精品久久成人aⅴ小说 | 涩涩av久久男人的天堂| 日韩av免费高清视频| 日韩在线高清观看一区二区三区| 久久久久久久精品精品| 黑丝袜美女国产一区| 老司机影院成人| 亚洲av国产av综合av卡| 中文字幕亚洲精品专区| 最近最新中文字幕免费大全7| 黄色欧美视频在线观看| 中国国产av一级| 看免费成人av毛片| 日韩免费高清中文字幕av| 成人毛片a级毛片在线播放| 亚洲欧美日韩另类电影网站| av福利片在线| 国产日韩欧美亚洲二区| 有码 亚洲区| 高清视频免费观看一区二区| 丝袜美足系列| 人成视频在线观看免费观看| 亚洲成人一二三区av| 亚洲怡红院男人天堂| av在线老鸭窝| av一本久久久久| 亚洲在久久综合| 毛片一级片免费看久久久久| 人妻夜夜爽99麻豆av| 国产精品国产三级国产专区5o| 少妇被粗大猛烈的视频| tube8黄色片| 丰满乱子伦码专区| 国产一区二区三区综合在线观看 | 日本午夜av视频| 国产黄色视频一区二区在线观看| 亚洲天堂av无毛| 新久久久久国产一级毛片| 国产男女内射视频| 亚洲,欧美,日韩| av国产久精品久网站免费入址| 国产女主播在线喷水免费视频网站| 人妻系列 视频| 国产欧美日韩综合在线一区二区| 免费少妇av软件| 丰满迷人的少妇在线观看| videosex国产| 女性被躁到高潮视频| 在线观看免费日韩欧美大片 | 观看美女的网站| 久久久久久久久久人人人人人人| av女优亚洲男人天堂| av播播在线观看一区| 亚洲不卡免费看| 2018国产大陆天天弄谢| 伊人久久精品亚洲午夜| 中文字幕精品免费在线观看视频 | 青春草亚洲视频在线观看| 日韩熟女老妇一区二区性免费视频| 久久国内精品自在自线图片| 国产极品天堂在线| 内地一区二区视频在线| 有码 亚洲区| 日本黄色日本黄色录像| 国产精品熟女久久久久浪| 在线观看免费日韩欧美大片 | 丰满乱子伦码专区| 日韩免费高清中文字幕av| 母亲3免费完整高清在线观看 | 久久国产亚洲av麻豆专区| 精品国产露脸久久av麻豆| 精品一区二区三区视频在线| 日本av免费视频播放| 亚洲少妇的诱惑av| 久久ye,这里只有精品| 成人亚洲精品一区在线观看| a 毛片基地| 最黄视频免费看| 一级爰片在线观看| 观看美女的网站| 肉色欧美久久久久久久蜜桃| 麻豆乱淫一区二区| 亚洲av.av天堂| 蜜桃在线观看..| 久久午夜福利片| 母亲3免费完整高清在线观看 | 人人妻人人添人人爽欧美一区卜| 夜夜看夜夜爽夜夜摸| 在线免费观看不下载黄p国产| 人人妻人人澡人人看| av专区在线播放| 一二三四中文在线观看免费高清| 人成视频在线观看免费观看| 中文字幕免费在线视频6| 日韩av免费高清视频| 成年人免费黄色播放视频| 五月天丁香电影| 日本午夜av视频| 成人综合一区亚洲| 一区二区三区四区激情视频| 我的老师免费观看完整版| 飞空精品影院首页| 国产国拍精品亚洲av在线观看| 久久鲁丝午夜福利片| 国产老妇伦熟女老妇高清| 欧美人与性动交α欧美精品济南到 | 亚洲精品久久久久久婷婷小说| 久久99热6这里只有精品| 高清黄色对白视频在线免费看| 人妻少妇偷人精品九色| a级毛片在线看网站| 日韩欧美精品免费久久| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 制服丝袜香蕉在线| 亚洲国产色片| 久久免费观看电影| 99热全是精品| 成人手机av| 成年av动漫网址| 色婷婷av一区二区三区视频| 一个人免费看片子| 精品久久久久久久久av| 成人午夜精彩视频在线观看| 亚洲国产精品一区三区| 看十八女毛片水多多多| 99久国产av精品国产电影| 22中文网久久字幕| 国产精品国产三级国产专区5o| 日本色播在线视频| 亚洲精品一区蜜桃| 亚洲国产欧美在线一区| h视频一区二区三区| 制服人妻中文乱码| 国产欧美亚洲国产| 丰满少妇做爰视频| 在线看a的网站| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片 | 日韩伦理黄色片| 高清在线视频一区二区三区| 大码成人一级视频| kizo精华| 中文字幕av电影在线播放| 女性生殖器流出的白浆| 国模一区二区三区四区视频| 亚洲av欧美aⅴ国产| 自线自在国产av| 久久久久人妻精品一区果冻| 一级毛片黄色毛片免费观看视频| 嘟嘟电影网在线观看| 午夜久久久在线观看| 久久精品人人爽人人爽视色| 亚洲欧美日韩另类电影网站| 亚洲一级一片aⅴ在线观看| 日本黄色日本黄色录像| 日韩 亚洲 欧美在线| 美女国产视频在线观看| 99re6热这里在线精品视频| 中文乱码字字幕精品一区二区三区| 美女主播在线视频| 另类亚洲欧美激情| 高清欧美精品videossex| av福利片在线| 国产精品人妻久久久久久| 大香蕉97超碰在线| 免费播放大片免费观看视频在线观看| 国产又色又爽无遮挡免| 99久久精品国产国产毛片| 久久精品国产自在天天线| 18+在线观看网站| 99久国产av精品国产电影| 久久久久精品久久久久真实原创| 亚洲精品一二三| 老司机影院成人| 日本wwww免费看| 99国产精品免费福利视频| 九色亚洲精品在线播放| 一本—道久久a久久精品蜜桃钙片| 看免费成人av毛片| 国产日韩欧美视频二区| 日韩三级伦理在线观看| videosex国产|