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

    水文變異條件下鄱陽湖流域的生態(tài)流量

    2015-01-18 07:39:45劉劍宇顧西輝
    生態(tài)學(xué)報 2015年16期
    關(guān)鍵詞:生態(tài)

    劉劍宇, 張 強(qiáng),3,*, 顧西輝

    1 中山大學(xué)水資源與環(huán)境系, 廣州 510275

    2 中山大學(xué)華南地區(qū)水循環(huán)與水安全廣東省普通高校重點(diǎn)實驗室, 廣州 510275

    3 宿州學(xué)院地球科學(xué)與工程學(xué)院, 宿州 234000

    水文變異條件下鄱陽湖流域的生態(tài)流量

    劉劍宇1,2, 張 強(qiáng)1,2,3,*, 顧西輝1,2

    1 中山大學(xué)水資源與環(huán)境系, 廣州 510275

    2 中山大學(xué)華南地區(qū)水循環(huán)與水安全廣東省普通高校重點(diǎn)實驗室, 廣州 510275

    3 宿州學(xué)院地球科學(xué)與工程學(xué)院, 宿州 234000

    受氣候變化和人類活動綜合影響,鄱陽湖流域水文狀況發(fā)生變異。河流生態(tài)系統(tǒng)適應(yīng)了變異前的水文狀況,變異后勢必會影響當(dāng)?shù)厣鷳B(tài)系統(tǒng)?;诖耍捎?種變異檢測方法對水文變異進(jìn)行綜合診斷,闡明水文變異原因。在此基礎(chǔ)上,采用15種概率分布函數(shù)分別擬合5站各月變異前日流量序列,最終確定5站點(diǎn)各月最優(yōu)分布函數(shù)及所對應(yīng)的概率密度最大處的流量,即得河道內(nèi)生態(tài)流量。研究表明:(1)撫河于1962年發(fā)生弱變異,贛江、修河于1968年發(fā)生中變異,信江、饒河于1991年發(fā)生弱變異;(2)變異后,贛江、信江、饒河、修河生態(tài)需水滿足率平均上升11%,撫河生態(tài)需水滿足率下降32%;(3)水文變異增加提高生態(tài)需水滿足率,水利工程建設(shè)降低年均生態(tài)需水滿足率、提高干季生態(tài)需水滿足率。高森林覆蓋率提高干季生態(tài)需水滿足率,對年均生態(tài)需水滿足率影響不明顯。研究結(jié)果為鄱陽湖流域水資源管理及區(qū)域水資源規(guī)劃與配置提供重要科學(xué)依據(jù)。

    鄱陽湖流域; 生態(tài)流量; 水文變異

    在全球氣候變化和人類活動共同作用下,河流水文過程發(fā)生顯著變異[1],改變了流域生態(tài)系統(tǒng)已適應(yīng)了的變異前的水文狀態(tài)[2]。為科學(xué)、合理計算河流生態(tài)需水,保障河流生態(tài)系統(tǒng)良性發(fā)展,有必要在水文變異的基礎(chǔ)上對河流生態(tài)需水進(jìn)行研究。崔瑛等[3]回顧了國內(nèi)外生態(tài)需水的研究進(jìn)展,認(rèn)為生態(tài)需水尚無統(tǒng)一認(rèn)可的概念。李捷等[4]提出逐月頻率計算法,通過與Tennant法對比得出該方法的優(yōu)缺點(diǎn)及適用范圍。張華等[5]運(yùn)用水文平衡原理構(gòu)建湖泊生態(tài)需水模型,估算了東居延海不同湖面面積下的生態(tài)需水量。李劍鋒等[6]提出考慮水文變異的河道內(nèi)生態(tài)需水計算方法,對黃河干流各站的月平均流量序列進(jìn)行生態(tài)需水計算。但相關(guān)研究尚未在鄱陽湖流域開展。鄱陽湖濕地是我國最大的淡水湖生態(tài)濕地,對維系區(qū)域和國家生態(tài)安全具有重要作用。受人類活動與氣候變化影響,該流域氣象水文過程發(fā)生顯著變異,變異后勢必會對當(dāng)?shù)匾约佰蛾柡^(qū)生態(tài)系統(tǒng)造成不同程度的影響。國務(wù)院于2009年正式批復(fù)《鄱陽湖生態(tài)經(jīng)濟(jì)區(qū)規(guī)劃》,這標(biāo)志著鄱陽湖生態(tài)經(jīng)濟(jì)區(qū)正式上升為國家戰(zhàn)略。因此,研究鄱陽湖流域生態(tài)需水更具有實踐意義。

    本文對鄱陽湖流域水文變異情況下,探討:1)采用M-K檢驗法、累積距平法、有序聚類法等八種變異診斷方法對變異點(diǎn)進(jìn)行系統(tǒng)檢測,增加變異檢驗的可信度;2)選用Log-Logistic分布、Weibull分布、Pearson-III型分布等15種概率分布函數(shù)分別對各站點(diǎn)各月份變異前的日流量序列進(jìn)行擬合分析;3)在擬合優(yōu)度檢驗方面,采用3種擬合優(yōu)度檢驗方法對概率分布函數(shù)進(jìn)行擬合優(yōu)度分析,更好的反映站點(diǎn)水文特征。并通過對變異前后生態(tài)需水滿足率的變化,分析討論水文變異對河道內(nèi)生態(tài)需水的影響。該項研究對科學(xué)分析計算河道內(nèi)生態(tài)需水具有重要理論意義,對鄱陽湖流域水資源管理具有重要參考價值。

    1 研究區(qū)概況與數(shù)據(jù)來源

    圖1 鄱陽湖流域、主要水文站點(diǎn)、水庫、灌區(qū)位置示意圖 Fig.1 Location of Poyang Lake,hydrological stations, water reservoirs and irrigation areas in the Poyang Lake Basin

    鄱陽湖是我國第一大淡水湖,流域面積16.22萬km2,占長江流域面積的9%。鄱陽湖入湖水量主要受贛江、撫河、信江、饒河、修河五大支流(以下簡稱“五河”)影響,構(gòu)成山江湖一體的核心-邊緣結(jié)構(gòu)體系[7](圖1)。鄱陽湖地處亞熱帶季風(fēng)氣候區(qū),降水量季節(jié)分配不均,流域大部分地區(qū)4—6月降水總量占全年降水總量的45%—50%。鄱陽湖流域共建有大中型水庫 208 座,總庫容 204.2 億 m3,10個大中型引水灌溉渠系,其中大于6 670 hm2以上的灌區(qū)有7個(圖1)。

    本文數(shù)據(jù)為鄱陽湖流域水系“五河”主要控制站的日流量(表1),數(shù)據(jù)來源于江西省水文局。部分缺失數(shù)據(jù)通過與相鄰的水文站水文序列建立回歸關(guān)系進(jìn)行插補(bǔ)(R2>0.8)。

    2 研究方法

    受氣候變化和人類活動的雙重影響,水文序列往往發(fā)生變異,使水文序列分布發(fā)生改變。河流生態(tài)系統(tǒng)適應(yīng)了變異前的水文狀態(tài),變異后勢必影響當(dāng)?shù)厣鷳B(tài)系統(tǒng)?;诖?,本文首先對研究區(qū)域河流水文序列進(jìn)行變異檢驗分析,再結(jié)合生物系統(tǒng)基礎(chǔ)理論對鄱陽湖流域生態(tài)需水進(jìn)行研究。

    表1 鄱陽湖流域水文站流量數(shù)據(jù)

    2.1 水文變異檢驗方法

    變異點(diǎn)檢驗方法眾多,李劍鋒等[6]使用滑動秩和檢驗法對黃河干流水文變異進(jìn)行了分析,馬嵐等[8]利用重新標(biāo)度極差分析法對石羊河下游徑流變異進(jìn)行了檢測。考慮到單一檢測方法有其適用流域的局限性,本文采用Hurst系數(shù)法對變異程度進(jìn)行初步分析,并選用M-K檢驗法、累積距平法、有序聚類法、Lee-heg法、滑動T檢驗法、滑動F檢驗法、滑動游程檢驗法和滑動秩和檢驗法八種方法對水文變異進(jìn)行綜合診斷。

    2.2 生態(tài)需水計算方法

    若水文變異診斷存在變異點(diǎn),則認(rèn)為水文序列的總體分布不一致,不符合水文資料的一致性要求??梢哉J(rèn)為流域生態(tài)環(huán)境適應(yīng)了變異前的水文狀態(tài),因此計算河道內(nèi)生態(tài)流量時,只考慮變異前的水文序列。若水文序列不存在變異點(diǎn),則采用整個水文序列計算生態(tài)需水。

    流量作為河流生態(tài)系統(tǒng)最重要的環(huán)境因子之一,可以作為河流生態(tài)環(huán)境的重要考察指標(biāo)。生態(tài)適宜性理論[9]認(rèn)為,生物在適宜的環(huán)境中,生物數(shù)量最多,生長最好;隨著生態(tài)因子偏離適宜值,生物生長繁殖速度減慢;當(dāng)超過生物耐受區(qū)間,生物數(shù)量即會減少,直至消亡。長期的自然選擇促使生物適應(yīng)出現(xiàn)頻率較高的環(huán)境因子。因此,本文取概率密度最大處流量作為適宜生物生長繁殖的生態(tài)流量。Green等[10]研究指出,任一種概率分布函數(shù)只可能對某種分布或者某個子樣本容量的檢驗效果較好,并不存在一種占絕對優(yōu)勢的概率分布函數(shù)。采用Log-Logistic分布、Weibull分布、Pearson-III型分布等[11]15種分布函數(shù)對5個站點(diǎn)變異前各月流量序列進(jìn)行擬合,并采用Kolmogorov-Smirnov(K-S)、Anderson Darling(A-D)和Chi-Squared(C-S)[12]3種檢驗方法進(jìn)行擬合優(yōu)度檢驗分析。

    3 變異點(diǎn)檢驗及成因分析

    3.1 變異點(diǎn)檢測

    采用Hurst系數(shù)法對變異程度作初步檢驗,虎山、李家渡、梅港、外洲、萬家埠三站Hurst系數(shù)分別為0.63、0.51、0.62、0.71、0.76,根據(jù)謝平等[13]基于Hurst系數(shù)對變異程度的劃分,虎山、李家渡、梅港三站為弱變異,外洲、萬家埠兩站為中變異。采用M-K檢驗法、累積距平法、有序聚類法等方法[14]對具體變異點(diǎn)進(jìn)行綜合檢測,將可能變異點(diǎn)中權(quán)重最高的年份作為最終變異點(diǎn)(表2)。

    表2 變異點(diǎn)綜合診斷結(jié)果

    虎山、梅港和李家渡三站的年徑流量均發(fā)生了弱變異,變異點(diǎn)分別是1991、1991、1962年;萬家埠和外洲兩站年徑流量發(fā)生了中變異,變異年份同為1968年。各站點(diǎn)水文序列變異點(diǎn)被多種方法同時檢測出,增加了變異點(diǎn)存在的可信度。

    3.2 水文變異成因分析

    整個鄱陽湖流域近50年來氣候變化明顯。鄱陽湖流域降雨量呈增加趨勢,蒸發(fā)量呈下降趨勢,降雨量在20世紀(jì)60年代中后期以及90年代初發(fā)生突變增加,蒸發(fā)量明顯減小[15]。運(yùn)用一元線性法對徑流序列進(jìn)行趨勢分析,發(fā)現(xiàn)信江、饒河、贛江、修河四河徑流序列都有不同程度的增加趨勢,這與鄱陽湖流域降水量突變增加、蒸發(fā)量突變減少基本相符,水利工程等會引起徑流量減小,因此,氣候變化是上述四河徑流突變增加的主要原因。相反,撫河徑流序列的線性系數(shù)為-1.318,趨勢減少,如果徑流突變由氣候變化主導(dǎo),那徑流變異應(yīng)該是突變增加,因此,撫河水文變異主要受人類活動的影響。

    各站點(diǎn)變異具體成因分析

    (1)信江、饒河 鄱陽湖流域年降水量在1990年突變增加,夏季降水量和暴雨頻率在1992年突變增加[16],與梅港、虎山兩站水文序列1991年變異時間點(diǎn)基本一致。饒河從水庫數(shù)量、規(guī)模和級別上都落后于其他流域,水利工程對降雨徑流的調(diào)蓄較少,徑流量受水利工程影響較小。因此,氣候變化是引起信江、饒河兩河水文變異(增加)的主要原因。

    (2)贛江、修河 東亞夏季風(fēng)在20世紀(jì)60年代中后期發(fā)生了一次突變,由正常偏強(qiáng)夏季風(fēng)轉(zhuǎn)為弱夏季風(fēng)[17]。夏季風(fēng)減弱,南方降水增多,平均氣溫降低,蒸發(fā)量減少,導(dǎo)致年降水量在60年代中后期發(fā)生突變。年降水量的增加和蒸發(fā)量的減少是引起兩河徑流變異增加。這與本文檢測出的贛江、修河徑流變異時間大致相同。

    (3)撫河 李家渡上游有贛撫平原灌區(qū)(1958—1960)、洪門水庫(1958—1961)。贛撫平原水利工程1958年動工,1960年秋開始收益,東西灌渠設(shè)計引水流量分別為60m3/s、107 m3/s,占撫河多年平均流量的42%;洪門水庫為江西四大水庫之一,始建于1958年,1961年開始蓄水,流域面積2736km2,最大庫容12億m3。本文檢驗出的變異時間與洪門水庫和贛撫平原灌渠建成時間相近,表明人類活動是導(dǎo)致?lián)岷铀淖儺惖闹饕颉?/p>

    4 結(jié)果分析

    4.1 生態(tài)需水計算結(jié)果

    根據(jù)變異點(diǎn)發(fā)生年份將水文序列分段,選用變異前序列計算河道內(nèi)生態(tài)需水。將概率密度最大處流量作為河道內(nèi)生態(tài)流量,求概率密度最大處流量。以贛江外洲站4月份為例,選用上述3種檢驗方法對15種分布函數(shù)進(jìn)行擬合優(yōu)度分析,綜合檢驗結(jié)果表明Johnson SB分布為外洲站4月份水文序列最優(yōu)概率分布函數(shù)(表3)。同時,以該月為例比較了幾個概率分布函數(shù)的理論概率分布曲線與經(jīng)驗概率分布(圖2),認(rèn)為Johnson SB法與經(jīng)驗概率分布吻合更好。

    表3 外洲站4月逐日流量的15種概率分布函數(shù)擬合優(yōu)度檢驗表

    圖2 外洲站4月逐日流量的理論與經(jīng)驗概率分布曲線及相應(yīng)的概率差異圖Fig.2 Probability distribution functions and Probability Difference for the respectively daily flow at the April of Waizhou station

    采用相同方法對其它各站點(diǎn)站逐月日流量序列進(jìn)行擬合優(yōu)度檢驗分析,得到相應(yīng)水文序列最優(yōu)分布函數(shù)與生態(tài)流量(表4)。

    選用15種分布函數(shù)對5個站點(diǎn)變異前水文序列進(jìn)行系統(tǒng)分析,并采用3種檢驗方法對擬合優(yōu)度綜合檢驗。檢驗分析發(fā)現(xiàn),不同子流域、不同月份最優(yōu)分布函數(shù)存在差異。因此,按不同流域、不同月份水文序列分別進(jìn)行概率分析是必要的。

    在擬合優(yōu)度綜合檢驗過程中發(fā)現(xiàn),某種分布函數(shù)在一種擬合優(yōu)度檢驗方法下是最優(yōu)的,但另外一種方法檢驗效果卻非常差。例如對外洲站5月份水文序列分布函數(shù)的擬合檢驗, K-S檢驗廣義Pareto分布為最優(yōu),A-D檢驗為最劣,C-S檢驗拒絕。故選用A-D檢驗、C-S檢驗最優(yōu)的,K-S檢驗第二的Johnson SB分布作為外洲站5月份最優(yōu)分布函數(shù)。因此,采用多種方法對概率分布函數(shù)擬合優(yōu)度進(jìn)行綜合檢驗更為合理。

    4.2 與其它方法對比

    將本文所采用的方法與最小月平均流量法、逐月最小流量法等常用方法比較。由于篇幅限制,只列出各方法所求得的年均生態(tài)流量(表5)。最小月平均流量法與最枯月頻率法所給出的生態(tài)流量各月相同,不能反映生態(tài)需水的月份差異;除外洲站外,各站生態(tài)流量的Tennant法[18]等級基本處于差水平。逐月最小流量法與本文方法都能反映生態(tài)需水的月份差異,但最小月平均流量法考慮的是最小月平均流量這一極端條件,因此,該方法下各月生態(tài)流量均小于本文生態(tài)流量。逐月頻率計算法下生態(tài)流量與本文生態(tài)流量相當(dāng),逐月頻率計算法在一定的保證率下對生態(tài)需水進(jìn)行計算,其保證率的設(shè)定至今沒有統(tǒng)一定論,計算結(jié)果主觀性較強(qiáng)。本文方法計算頻率最大處流量,計算結(jié)果在Tennant法等級中屬于好或中水平,可以滿足水生生物的正常需求,其計算結(jié)果具有確定性與合理性。

    表4 5站點(diǎn)各月份水文序列最優(yōu)分布函數(shù)及其對應(yīng)的生態(tài)流量

    表5 與其它生態(tài)需水計算方法比較

    5 討論

    如果實測流量大于生態(tài)流量,則認(rèn)為該流量滿足生態(tài)需水的要求。各月生態(tài)需水得以滿足的日數(shù)與序列相應(yīng)月份總?cè)諗?shù)之比為生態(tài)需水滿足率[19]。通過對各流域生態(tài)需水滿足率的對比分析來探討氣候變化與人類活動對生態(tài)需水的影響。

    5.1 變異前后年生態(tài)需水滿足率變化

    除撫河李家渡外,其它四站水文變異后生態(tài)需水滿足率均上升,平均上升幅度為11%,河流生態(tài)系統(tǒng)的正常需水得到更好保障(表6)。河流徑流量增加使得變異后生態(tài)需水滿足率提高。鄱陽湖流域1961—2003年饒河、信江、贛江中下游、信江中下游地區(qū)降水量呈增加趨勢,暴雨、特大暴雨頻次增加明顯對年降水量增加貢獻(xiàn)最大,同時流域內(nèi)蒸發(fā)量的減少一定程度上也增加了河川徑流[20]。變異后蒸發(fā)量的減少,降水量的增加是徑流量增加、生態(tài)需水滿足率提高的主要原因。

    撫河變異前,各月生態(tài)需水滿足率絕大部分在70%以上,說明撫河變異前70%以上的時間河流流量滿足生態(tài)需水的要求,能保障河流生態(tài)系統(tǒng)的正常需水。撫河流域水文變異后,各月生態(tài)需水滿足率均下降,近半數(shù)月份生態(tài)需水滿足率不到50%,平均生態(tài)需水滿足率下降32%。撫河徑流量的大量減少導(dǎo)致變異后生態(tài)需水滿足率降低。沿線灌區(qū)從撫河干支流上大量引水用于工農(nóng)業(yè)生產(chǎn),造成河道內(nèi)徑流量大幅度減少,使撫河生態(tài)需水滿足率變化方向與其它四河相反。李家渡生態(tài)需水滿足率下降與沿途灌渠大量引水有直接關(guān)系。

    表6 水文變異前后生態(tài)需水滿足率

    5.2 變異前后干濕季生態(tài)需水滿足率變化差異

    鄱陽湖流域4—6月是全年雨量最大的3個月份,所以把徑流量最大的4、5、6月定義為濕季,把徑流量最小的11、12月及次年1月定義為干季,對比分析變異后干濕季生態(tài)需水滿足率的差異。饒河虎山、贛江外洲、修河萬家埠3站干季生態(tài)需水滿足率高于濕季,撫河李家渡、信江梅港兩站生態(tài)需水滿足率干季小于濕季。流域森林植被對徑流變化有重要影響,森林可以起到蓄水、保水、保土作用,調(diào)節(jié)年際間枯水的流量,并且顯著提高干季徑流量。江西從20世紀(jì)80年代開始實行山江湖工程、退耕還林工程、中德造林工程、長(珠)防林工程以及各種人造林工程,大量植樹造林使江西省的森林覆蓋率迅速提高到60.1%,居全國第二位。饒河流域森林覆蓋率達(dá)67.8%,修河流域為64.4%,贛江流域63.6%,信江流域為54.3%,撫河流域為53.0%。

    饒河虎山站干季生態(tài)需水滿足率與年生態(tài)需水滿足率均為“五河”最高。饒河流域第一產(chǎn)業(yè)只占12.1%,農(nóng)業(yè)用地在“五河”中最少,土地利用類型以林地為主,饒河流域森林覆蓋率在“五河”中最高,占流域面積的67.8%,森林覆蓋率明顯增加干季徑流量,增大干季生態(tài)需水滿足率。饒河上游婺源、德興一帶是江西省的三大暴雨中心之一,該流域水資源量相當(dāng)?shù)刎S富。

    贛江外洲站年生態(tài)需水滿足率僅高于撫河,干季生態(tài)需水滿足率高出濕季39.1%,為“五河”干濕季生態(tài)需水滿足率差別最大的河流。水庫對徑流量變化具有重要的調(diào)節(jié)作用,使下游河道枯水流量相對穩(wěn)定,徑流量年內(nèi)分配較天然均勻,在一定程度上緩解干旱。贛江流域內(nèi)有水庫107座,庫容達(dá)75.13億m3,流域面積、水庫數(shù)量和容量均居五河之首,尤其是1990年修建的贛江流域最大水庫萬安水庫,其總庫容達(dá)22.16億m3,對枯水的調(diào)節(jié)發(fā)揮及其重要的作用。

    修河流域森林覆蓋率居鄱陽湖流域第二位,但干季生態(tài)需水滿足率僅高出濕季2.4%,這主要是由于萬家埠上游的潦河流域只有7座中型水庫,水庫以灌溉型為主,流域內(nèi)水庫數(shù)量、規(guī)模和級別均落后于其它流域,水庫調(diào)蓄功能小于其它流域。

    而撫河李家渡站、信江梅港站干季生態(tài)需水滿足率均低于濕季,分別低8.5%、3.5%。撫河流域耕地資源豐富,相應(yīng)的林業(yè)用地較少。撫河與信江兩流域森林覆蓋率遠(yuǎn)小于其它流域,森林調(diào)節(jié)能力有限是造成干季生態(tài)需水滿足率小于濕季的主要原因。撫河干、濕兩季生態(tài)需水滿足率的差距大于信江,一方面是撫河流域森林覆蓋率小于信江流域,另一方面撫河流域水庫18座遠(yuǎn)小于信江流域的35座,撫河流域森林與水庫的調(diào)節(jié)能力均小于信江。

    6 結(jié)論

    (1)受氣候變化和人類活動雙重影響,河流水文狀況發(fā)生不同程度變異?;⑸健⒚犯鬯男蛄芯l(fā)生弱變異,變異年份為1991年,變異主要原因為90年代初降水量的突變增加。外洲、萬家埠兩站水文序列發(fā)生中變異,變異年份為1968年,變異主要是受20世紀(jì)60年代中期夏季風(fēng)變化影響。撫河流域徑流于1962年發(fā)生弱變異,變異原因是受水利工程建設(shè)影響。

    (2)鄱陽湖流域水文變異提高了生態(tài)需水滿足率。另外,氣候變化與人類活動對生態(tài)需水影響的程度均不同,水利工程降低年均生態(tài)需水滿足率,而提高干季生態(tài)需水滿足率;高森林覆蓋率提高干季生態(tài)需水滿足率,而對年均生態(tài)需水滿足率影響不明顯;修河、饒河、信江、贛江徑流變異增加使年生態(tài)需水滿足率上升,撫河徑流變異減小使年生態(tài)需水滿足率大幅度減小。高森林覆蓋率與水利工程的大量建設(shè)使饒河虎山、贛江外洲、修河萬家埠3站干季生態(tài)需水滿足率高于濕季;低森林覆蓋率使信江流域、撫河流域干季生態(tài)需水滿足率小于濕季。

    (3)水文序列變異點(diǎn)檢測符合水文資料一致性要求,選用最優(yōu)分布函數(shù)所對應(yīng)流量作為河流的生態(tài)流量,符合水生生物生長需求。與Tennant法、最小月流量法、逐月頻率計算法等方法比較,本文計算結(jié)果更具確定性與合理性。本文研究對鄱陽湖流域各支流水資源管理與規(guī)劃提供重要科學(xué)依據(jù)。

    [1] Zhang Q, Xu C Y, Zhang Z X, Chen Y Q D. Changes of atmospheric water vapor budget in the Pearl River basin and possible implications for hydrological cycle. Theoretical and Applied Climatology, 2010, 102(1/2): 185- 195.

    [2] Ziegler A D, Sheffield J, Maurer E P, Nijssen B, Wood E F, Lettenmaier D P. Detection of intensification in global and continental-scale hydrological cycles: temporal scale of evaluation. Journal of Climate, 2003, 16(3): 535- 547.

    [3] 崔瑛, 張強(qiáng), 陳曉宏, 江濤. 生態(tài)需水理論與方法研究進(jìn)展. 湖泊科學(xué), 2010, 22(4): 465- 480.

    [4] 李捷, 夏自強(qiáng), 馬廣慧, 郭利丹. 河流生態(tài)徑流計算的逐月頻率計算法. 生態(tài)學(xué)報, 2007, 27(7): 2916- 2921.

    [5] 張華, 張?zhí)m, 趙傳燕. 極端干旱區(qū)尾閭湖生態(tài)需水估算——以東居延海為例. 生態(tài)學(xué)報, 2014, 34(8): 2102- 2108.

    [6] 李劍鋒, 張強(qiáng), 陳曉宏, 江濤. 考慮水文變異的黃河干流河道內(nèi)生態(tài)需水研究. 地理學(xué)報, 2011, 66(1): 99- 110.

    [7] 崔麗娟, 趙欣勝. 鄱陽湖濕地生態(tài)能值分析研究. 生態(tài)學(xué)報, 2004, 24(7): 1480- 1485.

    [8] 馬嵐, 魏曉妹. 石羊河下游年徑流序列的變異點(diǎn)分析. 干旱地區(qū)農(nóng)業(yè)研究, 2006, 24(2): 174- 177.

    [9] Mackenzie A, Ball A S, Virdee S R. BIOS Instant Notes in Ecology (Instant Notes Series). 2nd ed. UK: BIOS Scientific Publishers Ltd., 2001: 20- 63.

    [10] Green J R, Hegazy Y A S. Powerful modified-EDF goodness-of-fit tests. Journal of the American Statistical Association, 1976, 71(353): 204- 209.

    [11] 屠其璞, 丁裕國. 氣象應(yīng)用概率統(tǒng)計學(xué). 北京: 氣象出版社, 1984.

    [12] Sürücü B. A power comparison and simulation study of goodness-of-fit tests. Computers & Mathematics with Applications, 2008, 56(6): 1617- 1625.

    [13] 謝平, 雷紅富, 陳廣才, 李晶. 基于Hurst系數(shù)的流域降雨時空變異分析方法. 水文, 2008, 28(5): 6- 10.

    [14] 謝平, 陳廣才, 雷紅富, 武方圓. 水文變異診斷系統(tǒng). 水力發(fā)電學(xué)報, 2010, 29(1): 85-91.

    [15] 郭華, 蘇布達(dá), 王艷君, 姜彤. 鄱陽湖流域1955- 2002年徑流系數(shù)變化趨勢及其與氣候因子的關(guān)系. 湖泊科學(xué), 2007, 19(2): 163- 169.

    [16] Zhang Q, Peng J T, Xu C Y, Singh V P. Spatiotemporal variations of precipitation regimes across Yangtze River basin, China. Theoretical and Applied Climatology, 2014, 115(3/4): 703- 712.

    [17] 郭其蘊(yùn), 蔡靜寧, 邵雪梅, 沙萬英. 1873—2000年東亞夏季風(fēng)變化的研究. 大氣科學(xué), 2004, 28(2): 206- 215.

    [18] Tennant D L. Instream flow regimens for fish, wildlife, recreation, and related environmental resources. American Fisheries Society, 1976, 1(4): 6- 10.

    [19] 張強(qiáng), 李劍鋒, 陳曉宏, 江濤. 水文變異下的黃河流域生態(tài)流量. 生態(tài)學(xué)報, 2011, 31(17): 4826- 4834.

    [20] 孫鵬, 張強(qiáng), 陳曉宏. 鄱陽湖流域枯水徑流演變特征、成因與影響. 地理研究, 2011, 30(9): 1702- 1712.

    Evaluation of ecological flow with considerations of hydrological alterations in the Poyang Lake basin

    LIU Jianyu1,2, ZHANG Qiang1,2,3,*, GU Xihui1,2

    1DepartmentofWaterResourcesandEnvironment,SunYat-senUniversity,Guangzhou510275,China2KeyLaboratoryofWaterCycleandWaterSecurityinSouthernChinaofGuangdongHighEducationInstitute,SunYat-SenUniversity,Guangzhou510275,China3SchoolofEarthSciencesandEngineering,SuzhouUniversity,Suzhou234000,China

    Fluvial hydrological processes within the Poyang Lake basin are significantly altered under the influence of human activities and climate changes. Poyang Lake is the largest freshwater lake and plays an important role in conservation of biological diversity and also in flood mitigation in the Yangtze River basin. In this case, ecological instream streamflow is also altered and the ecological environment is potentially affected under the influences of altered hydrological processes. With consideration of hydrological alterations, re-evaluation of ecological instream streamflow will be of great scientific and practical merit in terms of water resource management and conservation of ecological environment. We use the Hurst coefficient method to preliminary analyze the degree of alterations in five major tributaries of the Poyang lake basin. Then we use eight mutation testing methods to comprehensively investigate change points. We use 15 kinds of probability distribution functions to fit respectively daily flow for each month before occurrence of the hydrological alterations. Finally, the monthly optimal distribution functions and corresponding streamflow with the largest probability are determined and computed, and accepted as the ecological instream streamflow. Results indicate that: (1) main causes behind the hydrological alterations of the Gan, Xin, Rao and Xiu Rivers are the climate changes such as the precipitation changes observed in this study. However, the Hu River is mainly affected by human activities; (2) due to large amounts of water diversion irrigation system, hydrological processes in the Hu River were significantly altered after 1962. The hydrological alterations of the Gan and Xiu Rivers occurred in 1968. The East Asian Summer Monsoon led to anincrease of precipitation during the mid-1960s; The Xin, Rao Rivers were dominated by significant hydrological alterations in 1991. Increase of precipitation in the early 1990s is the main cause behind this hydrological alterations that occurred in the 1960s;(3) The increase of forest coverage and the construction of water conservancy are expected to improve the satisfaction rate of the ecological water requirement in the dry season. The satisfactory rate of ecological water requirement in the dry season is higher than that in the wet season in the Rao, Gan, and Xiu Rivers. However, the Fu and Xin Rivers do not follow this trend. The forestation and increased vegetation coverage can greatly alter the spatiotemporal distribution of water or runoff and can significantly increase the runoff during the dry seasons. Woodland has an important influence on the changes of runoff. Due to the large forest coverage, the ecological instream flow is greatly satisfied during dry seasons in the Rao River when compared to the other four rivers. A difference is found in the Gan River in terms of water requirements of ecological instream flow between wet and dry seasons. The number of reservoirs in the Gan River basin is high when compared to the other four tributaries of the Poyang Lake basin that were considered in this study. Forest coverage in the Xiu River basin accounts for 64.4% of the total area, which is ranked the second largest in the Poyang Lake. However, the difference between wet and dry seasons is small. The quantity and scale of the reservoir in the Xiu River basin falls far behind others. Therefore, the regulation activities of the reservoirs in Xiu River basin is not evident than other river basins. On the contrary, the satisfactory rate of ecological water requirements for the dry season is below that for the wet season in the Fu and Gan River. The rate of forest coverage in those two basins is far less than in other basins. The ability to regulate the forest is limited. These results will provide an important scientific basis for the planning and management of water resources within the Poyang Lake basin under a changing environment.

    Poyang Lake basin; ecological streamflow; hydrological alterations

    鄱陽湖濕地與流域研究教育部重點(diǎn)實驗室(江西師范大學(xué))主任開放基金資助項目(ZK2013006); 安徽省高校引進(jìn)“領(lǐng)軍人才”專項項目

    2014- 04- 08;

    日期:2014- 09- 30

    10.5846/stxb201404080664

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

    劉劍宇, 張強(qiáng), 顧西輝.水文變異條件下鄱陽湖流域的生態(tài)流量.生態(tài)學(xué)報,2015,35(16):5477- 5485.

    Liu J Y, Zhang Q, Gu X H.Evaluation of ecological flow with considerations of hydrological alterations in the Poyang Lake basin.Acta Ecologica Sinica,2015,35(16):5477- 5485.

    猜你喜歡
    生態(tài)
    “生態(tài)養(yǎng)生”娛晚年
    住進(jìn)呆萌生態(tài)房
    生態(tài)之旅
    生態(tài)之旅
    生態(tài)之旅
    大營鎮(zhèn)生態(tài)雞
    貴茶(2019年3期)2019-12-02 01:46:32
    生態(tài)之旅
    潤豐達(dá) 微平衡生態(tài)肥
    茶葉通訊(2017年2期)2017-07-18 11:38:40
    生態(tài)保護(hù) 有你有我
    “知”與“信”:《逃逸行為》的生態(tài)自我
    成年人黄色毛片网站| 在线观看免费视频日本深夜| 一本久久中文字幕| 精品人妻一区二区三区麻豆 | 亚洲av电影不卡..在线观看| 欧美bdsm另类| 最新在线观看一区二区三区| 美女黄网站色视频| 午夜福利在线观看免费完整高清在 | 一夜夜www| 国产一区二区亚洲精品在线观看| 精品一区二区免费观看| 欧美日本亚洲视频在线播放| 少妇熟女aⅴ在线视频| 午夜精品久久久久久毛片777| а√天堂www在线а√下载| 亚洲性夜色夜夜综合| 免费电影在线观看免费观看| 少妇被粗大猛烈的视频| 永久网站在线| 国产视频一区二区在线看| 亚洲国产精品合色在线| 人人妻人人看人人澡| 在线免费十八禁| 国产亚洲91精品色在线| 成人毛片a级毛片在线播放| 九九爱精品视频在线观看| 日本 av在线| 国产黄片美女视频| 亚洲在线观看片| 老司机福利观看| 99国产精品一区二区蜜桃av| 亚洲av电影不卡..在线观看| 99久久九九国产精品国产免费| av天堂在线播放| 99在线人妻在线中文字幕| 国模一区二区三区四区视频| 国产精品永久免费网站| 又紧又爽又黄一区二区| 日日撸夜夜添| 毛片女人毛片| 人妻丰满熟妇av一区二区三区| .国产精品久久| 午夜福利欧美成人| eeuss影院久久| 亚洲欧美日韩高清在线视频| 人妻夜夜爽99麻豆av| 三级毛片av免费| 深爱激情五月婷婷| 精品一区二区免费观看| 午夜福利在线在线| 久久久久久国产a免费观看| 久久午夜亚洲精品久久| 99精品在免费线老司机午夜| 日韩一本色道免费dvd| 久久精品影院6| 18禁裸乳无遮挡免费网站照片| 免费黄网站久久成人精品| 熟妇人妻久久中文字幕3abv| 五月伊人婷婷丁香| 亚洲精品影视一区二区三区av| 国产麻豆成人av免费视频| 国产色婷婷99| 97碰自拍视频| 俄罗斯特黄特色一大片| 99久国产av精品| 国产精品亚洲美女久久久| 91午夜精品亚洲一区二区三区 | 欧美色欧美亚洲另类二区| 成人欧美大片| 日韩精品青青久久久久久| 女生性感内裤真人,穿戴方法视频| 久久久色成人| 少妇人妻精品综合一区二区 | 成年人黄色毛片网站| 成人国产麻豆网| 国内少妇人妻偷人精品xxx网站| 成人永久免费在线观看视频| 精华霜和精华液先用哪个| 一级毛片久久久久久久久女| 国产精品久久久久久久久免| 一个人看的www免费观看视频| 精华霜和精华液先用哪个| 99九九线精品视频在线观看视频| 成人鲁丝片一二三区免费| 韩国av一区二区三区四区| 欧美+亚洲+日韩+国产| 国产精品人妻久久久影院| 免费观看的影片在线观看| 精品人妻一区二区三区麻豆 | 亚洲电影在线观看av| 色吧在线观看| 成人性生交大片免费视频hd| 亚洲午夜理论影院| 三级男女做爰猛烈吃奶摸视频| 黄色丝袜av网址大全| 91午夜精品亚洲一区二区三区 | av在线观看视频网站免费| 22中文网久久字幕| 国产伦一二天堂av在线观看| 欧美高清成人免费视频www| 网址你懂的国产日韩在线| 色哟哟哟哟哟哟| 午夜影院日韩av| av在线观看视频网站免费| 久久热精品热| av国产免费在线观看| 国产精品人妻久久久影院| 久久久成人免费电影| 亚洲精品粉嫩美女一区| 超碰av人人做人人爽久久| 男插女下体视频免费在线播放| 99热这里只有精品一区| 国产精品不卡视频一区二区| 悠悠久久av| 一进一出好大好爽视频| 观看美女的网站| 日韩 亚洲 欧美在线| 久久久久久久久久久丰满 | 日韩中文字幕欧美一区二区| 欧美+亚洲+日韩+国产| 精华霜和精华液先用哪个| 成人国产麻豆网| 国产亚洲91精品色在线| 黄片wwwwww| 国产精品乱码一区二三区的特点| 亚洲一区二区三区色噜噜| 99久久精品热视频| 狂野欧美激情性xxxx在线观看| 亚洲欧美清纯卡通| 国产69精品久久久久777片| 久久精品国产鲁丝片午夜精品 | eeuss影院久久| 国产精品人妻久久久影院| 欧美一区二区亚洲| 一边摸一边抽搐一进一小说| 日本黄色片子视频| 国产aⅴ精品一区二区三区波| 国内精品一区二区在线观看| 婷婷色综合大香蕉| 99riav亚洲国产免费| 国产探花在线观看一区二区| eeuss影院久久| www.www免费av| 啦啦啦观看免费观看视频高清| 亚洲精品亚洲一区二区| 丰满人妻一区二区三区视频av| 精华霜和精华液先用哪个| 免费看光身美女| 久久香蕉精品热| 久久人人精品亚洲av| 美女黄网站色视频| 亚洲最大成人av| 亚洲四区av| 成人性生交大片免费视频hd| 亚洲不卡免费看| 国产高清不卡午夜福利| 亚洲无线在线观看| 成人亚洲精品av一区二区| 99视频精品全部免费 在线| 国产男人的电影天堂91| 国产69精品久久久久777片| 亚洲色图av天堂| 国产精品一区二区三区四区久久| 国产v大片淫在线免费观看| 又紧又爽又黄一区二区| 18禁在线播放成人免费| 又黄又爽又刺激的免费视频.| 又粗又爽又猛毛片免费看| 一夜夜www| 乱码一卡2卡4卡精品| 久久热精品热| 毛片一级片免费看久久久久 | 亚洲中文字幕日韩| 天美传媒精品一区二区| 中文字幕久久专区| 午夜福利欧美成人| 一区二区三区免费毛片| 桃色一区二区三区在线观看| 亚洲美女搞黄在线观看 | 精品久久久久久,| 婷婷六月久久综合丁香| 欧美黑人欧美精品刺激| 女人十人毛片免费观看3o分钟| 中亚洲国语对白在线视频| 国产单亲对白刺激| 免费无遮挡裸体视频| 国产精品一区二区三区四区久久| 91麻豆av在线| 波野结衣二区三区在线| 国国产精品蜜臀av免费| 极品教师在线视频| 又粗又爽又猛毛片免费看| 色播亚洲综合网| 国产精品一区二区免费欧美| 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| 久久久久国产精品人妻aⅴ院| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 久久香蕉精品热| 嫩草影院精品99| 黄色丝袜av网址大全| 亚洲 国产 在线| 国产精品久久久久久av不卡| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 少妇人妻精品综合一区二区 | 成人鲁丝片一二三区免费| 韩国av一区二区三区四区| 美女大奶头视频| 亚洲aⅴ乱码一区二区在线播放| 久久人妻av系列| 在线天堂最新版资源| 美女黄网站色视频| 嫁个100分男人电影在线观看| 欧美极品一区二区三区四区| 看十八女毛片水多多多| 最好的美女福利视频网| 亚洲18禁久久av| 久久久久九九精品影院| 亚洲成a人片在线一区二区| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 男人和女人高潮做爰伦理| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av| 97超级碰碰碰精品色视频在线观看| .国产精品久久| 九九热线精品视视频播放| 国产久久久一区二区三区| 日韩人妻高清精品专区| 亚洲四区av| 亚洲欧美日韩东京热| 国产精品人妻久久久影院| 国产精品久久久久久久久免| 99国产极品粉嫩在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲国产日韩欧美精品在线观看| 久久人人爽人人爽人人片va| 色av中文字幕| 久久人人精品亚洲av| 亚洲av免费高清在线观看| 日韩,欧美,国产一区二区三区 | 伊人久久精品亚洲午夜| 成人特级av手机在线观看| 在线播放国产精品三级| av中文乱码字幕在线| 春色校园在线视频观看| 成年版毛片免费区| 狂野欧美激情性xxxx在线观看| 精品久久久久久久久久久久久| 中文字幕av在线有码专区| 色噜噜av男人的天堂激情| 欧美xxxx性猛交bbbb| av天堂中文字幕网| 亚洲aⅴ乱码一区二区在线播放| 真人做人爱边吃奶动态| 欧美精品国产亚洲| 日韩欧美在线二视频| 亚洲精品一区av在线观看| 国内精品宾馆在线| 亚洲人成网站在线播放欧美日韩| 成人午夜高清在线视频| 午夜日韩欧美国产| 国产一区二区激情短视频| 国产精品久久久久久久久免| 亚洲精品一区av在线观看| 一本一本综合久久| 国产熟女欧美一区二区| 中亚洲国语对白在线视频| 尾随美女入室| 成人特级av手机在线观看| 天堂影院成人在线观看| 色噜噜av男人的天堂激情| 热99re8久久精品国产| 99视频精品全部免费 在线| 成人特级av手机在线观看| 国产色爽女视频免费观看| 亚洲精品日韩av片在线观看| xxxwww97欧美| 成人国产一区最新在线观看| 国内精品宾馆在线| 熟女电影av网| 精品久久久久久久末码| 色综合色国产| 91在线精品国自产拍蜜月| av专区在线播放| 亚洲av.av天堂| 无人区码免费观看不卡| 国产精品99久久久久久久久| 精品人妻一区二区三区麻豆 | 亚洲av美国av| 啪啪无遮挡十八禁网站| 99久国产av精品| 伦理电影大哥的女人| 男女做爰动态图高潮gif福利片| 国产精品久久久久久av不卡| 国产精品一区二区免费欧美| 日本色播在线视频| 成年免费大片在线观看| 久久久午夜欧美精品| 精品欧美国产一区二区三| 国产黄片美女视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲午夜理论影院| 成人特级av手机在线观看| 亚洲av成人精品一区久久| 国产主播在线观看一区二区| 国产精品久久电影中文字幕| avwww免费| 又爽又黄无遮挡网站| 男人舔奶头视频| 在线免费十八禁| 亚洲国产高清在线一区二区三| 久久亚洲真实| 天天一区二区日本电影三级| 69人妻影院| 欧美日本视频| 俄罗斯特黄特色一大片| 91狼人影院| 亚洲av日韩精品久久久久久密| 国内精品久久久久精免费| 国产91精品成人一区二区三区| 亚洲久久久久久中文字幕| 美女大奶头视频| 国产视频内射| 国语自产精品视频在线第100页| 男人舔女人下体高潮全视频| 悠悠久久av| 亚洲av免费在线观看| 99精品久久久久人妻精品| 久久久久久久亚洲中文字幕| 久久久久久九九精品二区国产| 乱人视频在线观看| 两人在一起打扑克的视频| 亚洲中文字幕一区二区三区有码在线看| 中文字幕熟女人妻在线| 99久久精品一区二区三区| 亚洲av成人av| 日日夜夜操网爽| 五月玫瑰六月丁香| 大又大粗又爽又黄少妇毛片口| 97超视频在线观看视频| 99热这里只有是精品50| 亚洲国产欧洲综合997久久,| 大又大粗又爽又黄少妇毛片口| 美女高潮喷水抽搐中文字幕| 午夜精品一区二区三区免费看| 日韩欧美三级三区| 亚洲欧美日韩高清在线视频| 少妇的逼好多水| 黄色欧美视频在线观看| 精品久久久久久久久av| 成人精品一区二区免费| 亚洲狠狠婷婷综合久久图片| 色综合婷婷激情| 国产一区二区在线av高清观看| 成人综合一区亚洲| 精品久久久久久久久av| 亚洲中文日韩欧美视频| 亚洲av日韩精品久久久久久密| 99热这里只有是精品50| 精品午夜福利在线看| 国产aⅴ精品一区二区三区波| 香蕉av资源在线| 亚洲美女搞黄在线观看 | 国产精品三级大全| 亚洲av中文字字幕乱码综合| 狂野欧美激情性xxxx在线观看| 中文字幕高清在线视频| 极品教师在线免费播放| 内射极品少妇av片p| 国产伦人伦偷精品视频| 成人国产综合亚洲| 国产精品久久视频播放| 一区二区三区激情视频| 1024手机看黄色片| 嫩草影视91久久| 精品无人区乱码1区二区| 日日啪夜夜撸| 成人二区视频| 精品国产三级普通话版| 中文在线观看免费www的网站| 99热网站在线观看| 亚洲在线观看片| 国产亚洲91精品色在线| 午夜视频国产福利| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 欧美成人性av电影在线观看| 婷婷亚洲欧美| 99riav亚洲国产免费| 国产老妇女一区| 最近中文字幕高清免费大全6 | 九九热线精品视视频播放| 赤兔流量卡办理| 久久久精品欧美日韩精品| 尾随美女入室| 狠狠狠狠99中文字幕| 少妇高潮的动态图| 性色avwww在线观看| 久久精品91蜜桃| 国产精华一区二区三区| 久久久久久九九精品二区国产| 欧美xxxx性猛交bbbb| 国产av在哪里看| 干丝袜人妻中文字幕| 美女免费视频网站| 日韩欧美在线二视频| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 美女黄网站色视频| 亚洲欧美日韩高清在线视频| 淫秽高清视频在线观看| 亚洲最大成人中文| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 最近最新免费中文字幕在线| 麻豆成人av在线观看| 性色avwww在线观看| 欧美日韩黄片免| 国产精品美女特级片免费视频播放器| 亚洲自拍偷在线| 成人鲁丝片一二三区免费| 日本色播在线视频| 一本精品99久久精品77| 精品久久国产蜜桃| 韩国av一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 亚洲精华国产精华液的使用体验 | 一区二区三区激情视频| 国产淫片久久久久久久久| 日韩欧美国产一区二区入口| 精品午夜福利视频在线观看一区| 国产精品一及| 亚洲av中文字字幕乱码综合| 如何舔出高潮| 亚洲精品成人久久久久久| 日本三级黄在线观看| 99热只有精品国产| 免费一级毛片在线播放高清视频| 黄色视频,在线免费观看| 亚洲午夜理论影院| 成人午夜高清在线视频| 久久香蕉精品热| 简卡轻食公司| 亚洲成人精品中文字幕电影| 亚洲国产精品久久男人天堂| 亚洲18禁久久av| 午夜免费男女啪啪视频观看 | 老司机福利观看| 高清日韩中文字幕在线| 美女被艹到高潮喷水动态| 黄色视频,在线免费观看| 日韩精品中文字幕看吧| 男插女下体视频免费在线播放| 午夜福利高清视频| 又爽又黄无遮挡网站| 最新中文字幕久久久久| 十八禁国产超污无遮挡网站| 丰满乱子伦码专区| 亚洲无线观看免费| 亚洲国产精品成人综合色| 欧美不卡视频在线免费观看| 欧美在线一区亚洲| 99热6这里只有精品| 国产成人av教育| 成熟少妇高潮喷水视频| 日韩欧美国产一区二区入口| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 成人三级黄色视频| 18禁黄网站禁片午夜丰满| 白带黄色成豆腐渣| 男女边吃奶边做爰视频| 99热这里只有精品一区| 国产午夜精品久久久久久一区二区三区 | 日本熟妇午夜| 国产免费av片在线观看野外av| 91麻豆av在线| 12—13女人毛片做爰片一| 欧美成人性av电影在线观看| 亚洲性久久影院| 亚洲av免费在线观看| 99久国产av精品| 欧美日本亚洲视频在线播放| 两人在一起打扑克的视频| 久久这里只有精品中国| 久久精品国产鲁丝片午夜精品 | 99热这里只有是精品50| 高清日韩中文字幕在线| 男女那种视频在线观看| 婷婷色综合大香蕉| 美女高潮的动态| 国产精品福利在线免费观看| 国产av一区在线观看免费| 99热网站在线观看| 人人妻人人看人人澡| 中文字幕熟女人妻在线| 欧美高清成人免费视频www| 夜夜爽天天搞| 欧美激情国产日韩精品一区| 淫秽高清视频在线观看| 给我免费播放毛片高清在线观看| 免费看av在线观看网站| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 很黄的视频免费| 亚洲色图av天堂| 黄色视频,在线免费观看| 国产又黄又爽又无遮挡在线| 精品人妻视频免费看| 亚洲精品久久国产高清桃花| 黄色一级大片看看| 国产精品爽爽va在线观看网站| 精品久久久久久久久久免费视频| 自拍偷自拍亚洲精品老妇| 日本爱情动作片www.在线观看 | 久久久久九九精品影院| 久久精品国产亚洲av涩爱 | 真人做人爱边吃奶动态| 国产精品,欧美在线| 国产亚洲精品av在线| 天美传媒精品一区二区| 午夜福利在线观看免费完整高清在 | 国产精品99久久久久久久久| 一区二区三区免费毛片| 国内精品美女久久久久久| 亚洲一区高清亚洲精品| 少妇丰满av| 欧美最黄视频在线播放免费| 狂野欧美激情性xxxx在线观看| 变态另类成人亚洲欧美熟女| 免费观看的影片在线观看| 国产精品无大码| 国产精品久久电影中文字幕| 99在线人妻在线中文字幕| 国产探花在线观看一区二区| 女人十人毛片免费观看3o分钟| 看十八女毛片水多多多| 国产美女午夜福利| 国产男靠女视频免费网站| 乱码一卡2卡4卡精品| 免费高清视频大片| 日本熟妇午夜| 少妇裸体淫交视频免费看高清| 国产熟女欧美一区二区| 能在线免费观看的黄片| 国产主播在线观看一区二区| 日韩欧美在线二视频| 色哟哟·www| 日韩,欧美,国产一区二区三区 | 国内久久婷婷六月综合欲色啪| 国产黄a三级三级三级人| 亚洲午夜理论影院| 天美传媒精品一区二区| 亚洲成人久久性| 两人在一起打扑克的视频| 亚州av有码| 午夜免费男女啪啪视频观看 | 亚洲欧美激情综合另类| 少妇被粗大猛烈的视频| 成人精品一区二区免费| 狠狠狠狠99中文字幕| 身体一侧抽搐| 国产中年淑女户外野战色| 高清日韩中文字幕在线| 啦啦啦韩国在线观看视频| 久久亚洲真实| 不卡一级毛片| 国产亚洲欧美98| 国产亚洲91精品色在线| 亚洲国产日韩欧美精品在线观看| 给我免费播放毛片高清在线观看| 久久久午夜欧美精品| 十八禁国产超污无遮挡网站| 长腿黑丝高跟| h日本视频在线播放| 日本a在线网址| 国产高潮美女av| 国产69精品久久久久777片| 国产男人的电影天堂91| 日韩在线高清观看一区二区三区 | 色综合婷婷激情| 久久久久国内视频| 国产高潮美女av| 99久久成人亚洲精品观看| 国产精品三级大全| 人人妻人人看人人澡| 丝袜美腿在线中文| 精品人妻1区二区| 国产久久久一区二区三区| 极品教师在线免费播放| 麻豆久久精品国产亚洲av| 蜜桃久久精品国产亚洲av| 日本精品一区二区三区蜜桃| 老熟妇乱子伦视频在线观看| 久久久久性生活片| 国产精品一区二区三区四区久久| 亚洲,欧美,日韩| 久99久视频精品免费| 国产精品无大码| 天堂av国产一区二区熟女人妻| 18禁裸乳无遮挡免费网站照片| 又粗又爽又猛毛片免费看| 男女那种视频在线观看| 69av精品久久久久久| 国产精品一区二区性色av| 日本免费一区二区三区高清不卡| 乱码一卡2卡4卡精品| 欧美激情在线99| 人人妻人人澡欧美一区二区| 99久久成人亚洲精品观看| 亚洲七黄色美女视频|