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

    環(huán)境變化對瀾滄江出境水文過程變異的影響

    2024-05-13 20:00:17王樂揚張建云寧忠瑞賈雨凡孫高霞王國慶
    水科學進展 2024年2期
    關鍵詞:人類活動徑流量氣候變化

    王樂揚 張建云 寧忠瑞 賈雨凡 孫高霞 王國慶

    摘要:變化環(huán)境改變河流水文節(jié)律,影響流域水資源的開發(fā)利用;瀾滄江-湄公河是位于東南亞源于中國西南的跨境河流,科學厘清瀾滄江徑流變化原因是實現(xiàn)瀾滄江-湄公河流域跨境水利益共享的重要基礎工作?;?961—2020年允景洪水文站的實測徑流量資料,采用數(shù)理統(tǒng)計與水文模擬相結合的研究途徑,診斷徑流序列的變異特征,定量分析氣候變化和人類活動對徑流及其年內分配的影響。結果表明:① 1961—2020年允景洪站實測年徑流量呈顯著性減少趨勢,其中,1—5月徑流量呈增加趨勢,6—12月徑流量呈減少趨勢;實測年徑流量在2008年發(fā)生突變。② RCCC-WBM模型能夠較好地模擬基準期(1961—2008年)的天然徑流量過程,模擬相對誤差小于1%;模擬的人類活動影響期間(2009—2020年)天然徑流量約為371.8 mm,較基準期偏少6.8%,這是氣候變化引起的。③ 水利工程建設運行是瀾滄江最主要的人類活動,是水庫蓄水運行期(2009—2013年)徑流減少的主要原因;氣候變化對水庫穩(wěn)定正常運行期(2014—2020年)徑流減少貢獻更大,約為61.7%。④ 水庫運行調度使得2014—2020年汛期(6—11月)徑流量減少44.3%,非汛期(12月至翌年5月)徑流量增加134.3%;瀾滄江干流水庫對徑流的削峰補枯調節(jié)作用在一定程度上保障了湄公河中下游的防洪安全與供水安全。

    關鍵詞:氣候變化;人類活動;徑流量;RCCC-WBM模型;歸因識別;瀾滄江流域

    中圖分類號:TV121.2

    文獻標志碼:A

    文章編號:1001-5791(2024)02-0313-12

    收稿日期:2023-10-08;網絡出版日期:2024-01-31

    網絡出版地址:https://link.cnki.net/urlid/32.1309.P.20240131.1028.002

    基金項目:國家自然科學基金資助項目(92047301;52121006)

    作者簡介:王樂揚(2000—),男,山東成武人,博士研究生,主要從事流域生態(tài)水文模型、水文水資源等方面的研究。

    E-mail:yueyangwang413@163.com

    通信作者:王國慶,E-mail:gqwang@nhri.cn

    氣候變化是流域水循環(huán)變化的直接驅動力,全球變暖影響下,蒸散發(fā)、土壤水和徑流量時空分布格局發(fā)生明顯變化[1-3;隨著經濟社會的快速發(fā)展,人類活動對區(qū)域水資源的影響愈加顯著[4-6??茖W識別不同驅動要素對流域水文的影響是國際水文科學協(xié)會(IAHS)2013—2022年計劃“Panta Rhei—everything flows”的重要主題[7-8,變化環(huán)境下的水資源問題成為多學科交叉研究的核心內容[9-10。

    在全球地緣政治經濟格局頻繁變動與淡水危機雙重作用下,國際河流的水資源變化直接關系到流域內各國家水資源開發(fā)利用的競爭關系,也影響著國家地緣安全[11-12。瀾滄江-湄公河是發(fā)源于中國青海省,流經緬甸、老撾、泰國、柬埔寨和越南的一條國際河流。變化環(huán)境下瀾滄江-湄公河流域不同區(qū)間徑流發(fā)生了較為明顯的改變,未來氣候變化下流域水文干旱情勢也更為嚴峻[13-15。位于中國境內的瀾滄江是湄公河的上游,其出境水量一直是下游流域國家關注的焦點,科學厘清環(huán)境變化對瀾滄江出境水文過程的影響,是解決瀾滄江-湄公河流域水事爭端、支撐中國跨境河流外交的重要基礎工作[16-17。

    徑流變化由多種驅動要素共同影響造成,歸因分析是研究各要素貢獻水平的一種方法。宋曉猛等[18將歸因分析定義為在一定可信度條件下,通過數(shù)學方法量化不同驅動因素對某系統(tǒng)變量演變的相對貢獻率,并可通過某些方式檢測出變量的變化過程或趨勢。近些年,水文學者在徑流變化歸因量化方法等方面做了許多有益探索,總體可劃分為4類:① 分項調查法,對水文過程的影響因子進行單獨調查,在此基礎上進行徑流還原計算,該方法概念清晰,但人力資源投入大[19-20;② 統(tǒng)計分析法,該方法對水文氣象觀測資料要求較高,其獲得可靠結論的前提是樣本充分,即連續(xù)觀測年份足夠多,一般需要30 a以上的資料,這在一定程度上限制了該方法的應用[21-22;③ 試驗流域法,該方法通過平行試驗流域對比分析,直觀明了地揭示流域治理對水文過程的影響,然而該方法試驗周期長,且很難找到相似的試驗流域,且大尺度流域不適用[23;④ 流域水文模型法,運用水文模型定量分析不同驅動因子的貢獻率,簡便且具有一定物理基礎,應用相對廣泛[24。Rsnen等[25基于1960—2014年的實測資料利用統(tǒng)計分析方法研究了水電站運行調度對瀾滄江河川徑流的影響,發(fā)現(xiàn)水電站運行對徑流量及其年內分配均會產生一定的影響;然而,該研究忽略了氣候要素變化對徑流的影響。Fan等[26采用水文模擬方法評估了降水、氣溫對湄公河上游和怒江流域徑流的影響,研究發(fā)現(xiàn)降水對這2個流域徑流的影響更為明顯,特別是汛期降水的影響更為顯著。流域水文情勢也在一定程度上受到水利工程以及其他人類活動的影響,基于徑流和氣候要素之間的相互作用研究對人類活動的影響考慮相對欠缺;因此,科學量化解析變化環(huán)境下不同驅動要素對年徑流及其月尺度徑流變化歸因亟待加強[27-28。

    隨著全球經濟一體化與可持續(xù)發(fā)展深度融合,人們越發(fā)關注國際跨境河流水量分配和管理利用問題[29-30。瀾滄江-湄公河是東南亞最大的跨境河流,位于中國境內的瀾滄江出境徑流量變化會在一定程度上影響到中下游的航運、灌溉及其他水資源利用,因此,瀾滄江出境水量已成為下游流域國家關注的重點和國際爭端的起因[31。本文在分析瀾滄江流域環(huán)境變化的基礎上,定量解析氣候變化和人類活動對瀾滄江流域徑流過程的影響,以期為實現(xiàn)瀾滄江-湄公河流域跨境水資源利益共享及解決水事爭端提供支撐。

    1 資料與方法

    1.1 流域概況及資料來源

    瀾滄江位于湄公河上游,發(fā)源于青海省唐古拉山,界于94°E—102°E、21°N—34°N,由北向南呈條帶狀,流經青海省、西藏自治區(qū)和云南省,在西雙版納州勐臘縣關累鎮(zhèn)西南流出中國國境,中國境內主干流總長度為2 130 km,允景洪水文站以上流域面積為141 779 km2。圖1給出了瀾滄江流域水系及水文氣象站點位置示意。

    瀾滄江流域由北向南跨越13個緯度,氣溫由北向南遞增,縱跨不同的氣候帶。其中,位于青海南部的源頭地區(qū)屬于高寒氣候區(qū),位于西藏的中游上段屬于高原溫帶氣候區(qū),中游滇西北區(qū)屬亞熱帶,下游滇西南地區(qū)丘陵和盆地交錯分布,屬于亞熱帶或熱帶氣候。瀾滄江流域支流甚多,流域面積大于1 000 km2的支流有41條,較大支流主要分布在上游和下游(圖1)。瀾滄江河流天然落差大,約為4 583 m,允景洪站多年平均徑流量約540億m3,流域水資源及水能資源均非常豐富。

    本文收集整理了允景洪站逐月徑流量資料,瀾滄江流域內及附近周邊17個氣象站點氣溫、降水資料以及12個雨量站點的降水資料,其中氣象、水文資料的序列為1961—2020年。由于瀾滄江流域面積較大,以0.5°分辨率進行網格劃分,采用分布式計算模擬流域水文過程;基于氣象站點、雨量站點的降水資料,采用線性插值方法將站點數(shù)據(jù)插值到網格[32。從中國科學院資源環(huán)境科學與數(shù)據(jù)中心收集了瀾滄江流域土地利用資料,從湄公河委員會官方網站、中國華能瀾滄江水電集團網站和開放發(fā)展湄公河計劃網站收集了瀾滄江截止到2020年的已建水電工程資料。

    1.2 研究方法

    徑流變化是氣候變化和人類活動兩大驅動因子綜合作用的結果?;谒哪M解析不同驅動因子對徑流影響包括4個主要步驟[19,24:① 采用Mann-Kenall方法、有序聚類方法等數(shù)理統(tǒng)計方法診斷徑流系列的趨勢性及突變特征,將徑流劃分為弱人類活動影響的“天然階段”和“人類活動影響階段”,將“天然階段”作為基準期;② 利用天然時期的水文氣象資料率定模型參數(shù),驗證流域水文模型在研究流域的適應性,只有流域水文模擬對天然徑流過程具有較好的模擬效果,才可用于歸因量化識別;③ 采用人類活動影響時期的氣象資料驅動率定的流域水文模型,模擬“人類活動影響階段”的天然徑流量過程;④ 以人類活動影響期間還原的天然徑流量與基準時期天然徑流量之差作為氣候變化對徑流的影響量,以人類活動影響期間還原的天然徑流量與相應時期的實測徑流量之差作為人類活動對于徑流的影響量。

    Mann-Kendall(M-K)檢測方法是一種非參數(shù)序列趨勢診斷方法[33-34,通過計算M-K標準化統(tǒng)計量(UFK值)確定序列變化的趨勢及其顯著性。當統(tǒng)計量為正或負時說明序列呈現(xiàn)增加或減少趨勢,當統(tǒng)計量介于給定顯著水平α=0.05的臨界值±1.96區(qū)間內,說明系列為非顯著性變化趨勢;反之,則變化趨勢顯著。有序聚類方法是常用的序列突變診斷方法,其原理是以某一時間點()將系列劃分為不同系列,最優(yōu)分割點*使得2個系列的離差平方和最小[34。

    RCCC-WBM模型是由水利部應對氣候變化研究中心團隊研發(fā)并逐步完善的大尺度流域水文模型[35-36。該模型考慮了降雨徑流、融雪徑流和地下徑流3種徑流,能夠較好地模擬高寒、濕暖等不同氣候區(qū)域的水文過程。模型輸入包括逐時段降水量、蒸發(fā)能力和氣溫,其中,蒸發(fā)能力一般采用E601實測值,也可采用Penman-Monteith等公式進行估算[37。針對大尺度流域,模型在對流域進行網格劃分和網格產流計算的基礎上,采用VIC模型的匯流計算方案進行大尺度流域水文過程模擬[37。選用Nash-Sutcliffe模型效率系數(shù)(ENS)和模擬總量相對誤差(ER)為目標函數(shù)進行參數(shù)率定,若ENS越接近于0,同時ER越接近于1,說明對徑流過程的模擬效果越好[38

    2 結果與討論

    2.1 瀾滄江流域氣候及環(huán)境變化

    2.1.1 土地利用類型變化

    基于中國科學院資源環(huán)境科學與數(shù)據(jù)中心1980年、2000年和2015年3期土地利用資料,分析了3個時期瀾滄江流域耕地、林地、灌木、草地、水體、未利用地和建設用地的面積及不同土地利用類型的轉移變化(圖2)。

    由圖2可以看出:① 流域主要土地利用類型為草地、林地和灌木3類;草地面積最大,廣泛分布于流域上游,占流域總面積的43%(1980年);林地次之,主要分布于流域中下游地區(qū),占流域總面積的28%(1980年);灌木面積次于林地面積,占流域總面積的14%(1980年)。② 草地、林地和未利用地面積變化較大,但其相對變化較小,均不到0.7%;相比而言,水體和建設用地面積相對變化較大,1980—2015年,水體面積減少15%,建設用地面積增加1.1倍,但二者面積變化較小,不到流域總面積的1%。③ 從不同土地利用類型轉變面積來看,耕地、林地、草地和灌木用地轉移變化最大,且主要是四者之間相互轉化。譬如,1980—2000年,林地多轉變?yōu)楦亍⒉莸睾凸嗄镜?,同時,耕地、草地和灌木地也是轉變?yōu)榱值氐闹匾M成??傮w來看,在過去幾十年,瀾滄江流域盡管存在不同土地利用類型的轉化,但不同土地利用類型面積變化相對較小。

    2.1.2 水能資源開發(fā)程度變化

    瀾滄江流域具有十分豐富的水能資源,蘊藏著約3 656萬kW的水能,其中干流上約2 545 kW,為中國重點開發(fā)的13個水電基地之一。瀾滄江流域內已修建有大朝山(2003年投產運行)、景洪(2008年投產運行)、小灣(2009年投產運行)、功果橋(2011年投產運行)、糯扎渡(2012年投產運行)等大型水電站,其中,小灣和糯扎渡調節(jié)庫容較大[39。圖3給出了瀾滄江截至到2020年已建水電工程的累積庫容??梢钥闯觯?008年之前瀾滄江處于低程度開發(fā),累積庫容約18億m3,不到年徑流量的4%;之后,隨著景洪、小灣等水電站工程的建設運行,2008—2013年累積庫容增加明顯,后期累積庫容相對穩(wěn)定。2020年已建水電站工程庫容接近多年平均年徑流量,水電工程的徑流調節(jié)能力顯著增強。

    2.1.3 氣溫和降水演變

    采用M-K檢測方法診斷了流域平均氣溫和降水演變趨勢的顯著性。圖4給出了瀾滄江流域1961—2020年氣溫、降水演變過程及其趨勢顯著性診斷結果。

    瀾滄江流域多年平均氣溫為13.1 ℃,多年平均年降水量為1 094.2 mm。由圖4可以看出,1961—2020年期間,氣溫呈現(xiàn)顯著的上升趨勢,線性趨勢率為0.27 ℃/(10 a)。降水以自然波動為主,其中,1975—1982年為枯水段,8 a中只有2個年份降水量略多于多年均值;1993—2002年降水相對較豐,10 a中只有1994年、1997年降水量略低于多年均值;2009—2018年又是一個持續(xù)少水階段,10 a中只有2010年和2016年降水略多于多年均值。趨勢顯著性診斷結果表明:各月氣溫均呈現(xiàn)顯著性上升趨勢;不同月份降水量均為非顯著性變化趨勢,其中,1月、3—5月和9月降水量呈現(xiàn)增加趨勢,6—8月降水量為減少趨勢,其余月份降水量基本不存在趨勢性變化;年和汛期降水量均呈現(xiàn)減少趨勢,相比而言,汛期降水量減少相對明顯。

    瀾滄江流域上下游氣候差異較大,降水自南向北具有遞減趨勢,下游降水量較大,部分站點年降水量超過1 800 mm,上游降水相對偏少,個別站點降水量不足500 mm;氣溫也呈現(xiàn)自南向北遞減態(tài)勢。在瀾滄江上、中、下游不同區(qū)域選取昌都、德欽、思茅3個典型氣象站點,診斷分析了1961—2020年氣溫和降水的變化趨勢特征(表1)。

    統(tǒng)計結果表明,昌都、德欽和思茅3個站多年平均降水量分別為483.2、633.4和1 488.1mm,上、中游的昌都站和德欽站年均氣溫均在10 ℃以下,而下游思茅站年均氣溫約為18.5 ℃,遠高于中上游的2個站點。由表1可以看出,1961—2020年3站均呈現(xiàn)顯著升溫態(tài)勢,中、下游的德欽站和思茅站線性升溫幅度更大,均超過0.4 ℃/(10 a);就降水的年代際趨勢來看,上游昌都站呈現(xiàn)非顯著性增加趨勢,中、下游的2個代表站降水量呈現(xiàn)非顯著性減少趨勢,其中,下游的思茅站線性遞減率更大,約為-25.10 mm/(10 a)。

    2.2 變化環(huán)境下徑流過程變異特征

    瀾滄江允景洪站多年平均徑流深為381.6 mm,變化環(huán)境下實測徑流量過程發(fā)生一定的變化。圖5(a)給出了瀾滄江允景洪站1961—2020年實測徑流量距平及其累積過程,采用M-K檢測方法診斷了月徑流量演變趨勢的顯著性(圖5(b))。

    由圖5(a)可以看出:允景洪站實測徑流在2000年之前以自然波動為主,其中,20世紀60年代水量偏豐,多數(shù)年份高于多年均值;2000年之后,徑流量相對偏少,只有2001年、2004年、2007年和2017年4個年份徑流量高于多年均值,其中2008—2013年實測徑流量持續(xù)偏少,且偏少幅度較大,接近-100 mm。圖5(b)趨勢診斷結果表明:1—5月徑流量呈現(xiàn)增加趨勢,6—12月徑流量呈現(xiàn)減少趨勢,其中,只有1月、2月、6月、12月徑流量為顯著性變化趨勢;統(tǒng)計結果表明,1961—2020年期間汛期(6—11月)和年徑流量均為顯著性減少趨勢。

    基于聚類分析方法對年徑流量突變診斷結果表明,實測徑流量在2008年發(fā)生了較為明顯的突變,結合瀾滄江流域內人類活動狀況(圖3),將1961—2008年作為人類活動影響相對較小的天然基準期,將2009—2020年作為人類活動影響時期。統(tǒng)計結果表明,2009—2020年實測徑流量較基準期偏少21.9%,其中,2009—2013年和2014—2020年期間分別偏少34.1%和13.2%,2009—2013年實測徑流偏少較多,這與圖5(a)徑流量距平給出的結果一致。

    根據(jù)瀾滄江流域降水和徑流的年內分布特征,將6—11月和12月至翌年5月分別定義為汛期和非汛期,圖6給出了允景洪站2個時期實測徑流量較基準期徑流量變化的年內分配過程。可以看出,2009—2013年非汛期徑流量較基準期變化不大,但2014—2020年期間非汛期徑流量較基準期明顯偏大;同時,也可發(fā)現(xiàn)2個時期汛期徑流量均較基準期明顯偏小,特別是主汛期7—10月徑流偏小幅度較大。

    2.3 模型參數(shù)率定及天然徑流過程模擬

    利用1961—2008年的逐月氣象水文資料率定RCCC-WBM模型參數(shù)并驗證模型的可靠性。為避免初始狀態(tài)變量對模擬效果的影響,將1961年作為模型預熱期,將1962—1990年作為率定期率定模型參數(shù),將1991—2008年作為驗證期驗證模型的可靠性。圖7給出了1961—2008年允景洪站實測與模擬徑流量年內分配過程。

    由圖7可以看出,1961—2008年瀾滄江允景洪站實測與模擬月徑流量過程擬合良好,7—8月模擬徑流量略微偏大,10—11月模擬徑流量略微偏小。統(tǒng)計結果表明,年徑流量模擬誤差為0.7%,率定期和驗證期的ENS分別為0.87和0.85,ER分別為-3.1%和1.4%,說明RCCC-WBM模型對允景洪站月徑流量過程具有較好的模擬能力。

    保持模型參數(shù)不變,利用1961—2020年的實測氣象資料驅動RCCC-WBM模型,模擬天然狀態(tài)下流域的徑流量過程。圖8給出了1961—2020年瀾滄江允景洪站天然徑流量模擬偏差及其累積過程??梢钥闯?,2008年之前天然徑流量模擬偏差正負相間且變幅較小,多在±10.0 mm之間,2008年之后模擬值偏大,個別年份甚至偏大100.0 mm以上。模擬偏差的累積過程在2008年之前圍繞0軸線上下波動,2008年之后由于模擬值偏大,模擬偏差累積過程持續(xù)升高,在2014年之后有所放緩。統(tǒng)計結果表明,2009—2020年模擬的天然徑流量為371.8 mm,盡管較同期實測徑流量偏大,但仍較1961—2008年偏少6.8%。

    2.4 環(huán)境變化對徑流變異的影響

    基于實測徑流距平演變特征及突變診斷結果和圖3分析結果,以1961—2008年為基準期,將人類活動影響期分段定義為水庫蓄水運行期(2009—2013年)和穩(wěn)定調度運行期(2014—2020年);根據(jù)實測徑流和模擬的天然徑流過程分析氣候變化和人類活動對徑流量的影響(表2)。

    由表2可以看出:① 與基準期相比,2009—2013年和2014—2020年實測徑流量分別減少136.2、52.8 mm,氣候變化和人類活動均使得2個影響時期徑流量減少;② 2009—2013年期間,氣候變化和人類活動的影響導致徑流量分別減少19.9、116.3 mm,占徑流減少總量的14.6%和85.4%,人類活動是徑流減少的主要原因;③ 2014—2020年期間,氣候變化和人類活動的影響導致徑流量分別減少32.6、20.2 mm,占該時期徑流減少總量的61.7%和38.3%,氣候變化是該時期徑流減少的主要驅動因素;④ 對比2個影響時期,氣候變化對徑流的影響量由-19.9 mm變化到-32.6 mm,影響幅度有增大趨勢;而人類活動對徑流的影響量由-116.3 mm變化到-20.2 mm,影響幅度大幅度減小。

    環(huán)境變化不僅影響徑流量,而且影響徑流的年內分配。圖9給出了氣候變化和人類活動對允景洪站2009—2013年和2014—2020年2個時期徑流年內分配的影響。

    由圖9可以看出:① 2009—2013年期間,氣候變化導致5月、7月和11月徑流略微增加,其他月份徑流減少,9月減少幅度最大,約為-10.1 mm。在2014—2020年期間,氣候變化因素導致4—6月和11月徑流量增加,但增加幅度較小,均不到1.5 mm;氣候變化使得其他月份徑流減小,相比而言,8月減少幅度最大,約為-15.8 mm,其次為9月和7月,減少幅度分別為-8.8、-4.2 mm。② 水庫建設運行調度是瀾滄江流域最主要的人類活動;在水庫蓄水運行期(2014—2020年),人類活動使得汛期徑流量減少42.2%,非汛期徑流量增加6.4%;在水庫穩(wěn)定調度運行期(2014—2020年),人類活動的影響使得徑流量減少44.3%,非汛期徑流量增加134.3%。③ 總體來看,氣候變化對主汛期的8—9月徑流影響幅度較大,而對非汛期月份徑流影響幅度較小;水利工程運行調度在蓄水階段對汛期影響較大,盡管在穩(wěn)定調度運行階段對年徑流量的影響較小,但削峰增枯的作用明顯。

    2.5 討論

    氣溫升高導致蒸發(fā)能力和流域蒸發(fā)損失增大,從而導致徑流減少;降水是河川徑流的主要來源,降水的增減直接影響徑流的豐枯。瀾滄江流域2009—2013年和2014—2020年氣溫較1961—2008年分別升高0.97、1.03 ℃,降水量分別減少11.8、23.5 mm,這是徑流減少的原因之一,其中,后一個時期氣溫升高及降水減少幅度均大于前一個時期,這是氣候變化影響量幅度增大的重要原因。

    下墊面變化、工農業(yè)發(fā)展及水電工程開發(fā)是瀾滄江流域主要的人類活動結果。前述土地利用分析結果表明,瀾滄江流域過去40 a盡管土地利用類型發(fā)生了變化,但變化較小;基于本研究收集到的用水資料統(tǒng)計得到,2010年和2020年瀾滄江流域用水總量分別為27.6、26.7億m3,約占允景洪站年徑流量的5%,流域內用水量總體變化不大,并且占徑流總量的比率較低,因此,土地利用變化和用水量變化不會對流域水文過程產生較大的影響。然而,瀾滄江干流修建了一系列水電工程,特別是2008—2014年期間,已建水庫的累積庫容大幅度增加,而實測徑流也在2008年前后發(fā)生突變,因此,水利工程建設與運行是流域的主要人類活動,且與徑流變化密切相關。

    圖10給出了瀾滄江流域徑流突變前后3個時段的降水—徑流關系,可以看出:基準期1961—2008年的降水—徑流關系點群最高,2009—2013年的點群最低,2014—2020年的點群大多位于基準期點群之中,但位置相對略低。由此說明,同樣降水條件下基準期的徑流量最大,其次為2014—2020年,2009—2013年的徑流量最小。結合圖3不難推斷,2009—2013年期間是瀾滄江水庫蓄水運行期,該時期水庫蓄水調度對徑流量的影響相對較大,其影響量約占該時期天然徑流量的30%;2014—2020年是水庫穩(wěn)定調度期,水庫穩(wěn)定調度對徑流量的影響相對較小,影響量約占該時期天然徑流量的5.3%。

    湄公河流域雨熱同期,降水主要集中在汛期6—10月,非汛期1—5月降水相對偏少,而非汛期正是作物生長季節(jié),需水量大[13-14。基于本研究收集的氣象資料統(tǒng)計得到,1961—2008年,70%以上的徑流量集中在汛期6—10月,非汛期徑流量不到年徑流量的30%;而在2009—2013年和2014—2020年,由于水庫的調節(jié)作用,非汛期徑流量占比增加至年徑流量的40.3%和61.6%。瀾滄江干流水庫的修建和運行盡管對徑流量有一定的影響,但其對徑流的調節(jié)作用對于保障湄公河中下游的農業(yè)灌溉起到至關重要的作用。

    3 結論

    基于瀾滄江流域水文氣象、水利工程、工農業(yè)用水和下墊面等資料,分析了流域氣候與徑流的演變特征;采用流域水文模擬與數(shù)理統(tǒng)計分析相結合的途徑,量化了氣候要素和人類活動對允景洪站徑流變化的影響。主要結論如下:

    (1) 1961—2020年,瀾滄江流域氣溫呈現(xiàn)顯著上升趨勢,降水以豐枯交替波動為主;變化環(huán)境影響下,允景洪站實測年和汛期徑流量出現(xiàn)顯著減少趨勢,1—5月徑流量出現(xiàn)增加趨勢,6—12月徑流量出現(xiàn)減少趨勢;實測年徑流量在2008年發(fā)生突變,2009—2020年實測徑流量較1961—2008年徑流量減少21.9%。

    (2) RCCC-WBM模型對允景洪站天然徑流量具有較好的模擬效果,率定期(1961—1990年)和驗證期(1991—2008年)的Nash-Sutcliffe模型效率系數(shù)分別為0.87和0.85,模擬相對誤差分別為-3.1%和1.4%,可見該模型可以用于人類活動影響時期天然徑流量的模擬還原。

    (3) 人類活動是2009—2013年徑流減少的主要原因,其貢獻占比約為85.4%;氣候變化是2014—2020年徑流減少的主要驅動因素,其貢獻占比約為61.7%;氣候變化對徑流的影響幅度有增大趨勢,人類活動對2個階段徑流的影響幅度顯著減小,由-116.3 mm(2009—2013年)變化到-20.2 mm(2014—2020年)。

    (4) 以水電工程為主的人類活動對允景洪站徑流量的影響主要體現(xiàn)為增枯削峰作用,使得非汛期1—5月和12月徑流增加,6—11月徑流減少。這一方面增大了湄公河中下游旱季灌溉用水量,另一方面通過攔蓄汛期徑流,減輕了湄公河中下游防汛抗洪壓力。

    參考文獻:

    [1]IPCC.Climate change 2021:the physical science basis[M].Cambridge:Cambridge University Press,2021.

    [2]KIESEL J,GERICKE A,RATHJENS H,et al.Climate change impacts on ecologically relevant hydrological indicators in three catchments in three European ecoregions[J].Ecological Engineering,2019,127:404-416.

    [3]張建云,王國慶,金君良,等.1956—2018年中國江河徑流演變及其變化特征[J].水科學進展,2020,31(2):153-161.(ZHANG J Y,WANG G Q,JIN J L,et al.Evolution and variation characteristics of the recorded runoff for the major rivers in China during 1956—2018[J].Advances in Water Science,2020,31(2):153-161.(in Chinese))

    [4]FENG X M,F(xiàn)U B J,PIAO S L,et al.Revegetation in China′s Loess Plateau is approaching sustainable water resource limits[J].Nature Climate Change,2016,6(11):1019-1022.

    [5]李艷忠,劉昌明,劉小莽,等.植被恢復工程對黃河中游土地利用/覆被變化的影響[J].自然資源學報,2016,31(12):2005-2020.(LI Y Z,LIU C M,LIU X M,et al.Impact of the grain for green project on the land use/cover change in the Middle Yellow River[J].Journal of Natural Resources,2016,31(12):2005-2020.(in Chinese))

    [6]王樂揚,李清洲,王金星,等.變化環(huán)境下近60年來中國北方江河實測徑流量及其年內分配變化特征[J].華北水利水電大學學報(自然科學版),2020,41(2):36-42.(WANG Y Y,LI Q Z,WANG J X,et al.The variation characteristics of recorded runoff and its annual distribution in North China during the recent 60 years in the context of environment change[J].Journal of North China University of Water Resources and Electric Power (Natural Science Edition),2020,41(2):36-42.(in Chinese))

    [7]MONTANARI A,YOUNG G,SAVENIJE H H G,et al.“Panta Rhei:everything flows”:change in hydrology and society:the IAHS scientific decade 2013—2022[J].Hydrological Sciences Journal,2013,58(6):1256-1275.

    [8]BLSCHL G,BIERKENS M,CHAMBEL A,et al.Twenty-three unsolved problems in hydrology (UPH):a community perspective[J].Hydrological Sciences Journal,2019,64:1141-1158.

    [9]EVARISTO J,MCDONNELL J J.Global analysis of streamflow response to forest management[J].Nature,2019,570(7762):455-461.

    [10]POKHREL Y,F(xiàn)ELFELANI F,SATOH Y,et al.Global terrestrial water storage and drought severity under climate change[J].Nature Climate Change,2021,11:226-233.

    [11]鐘蘇娟,毛熙彥,黃賢金.地緣安全視角下的中國國際河流水資源開發(fā)利用[J].世界地理研究,2022,31(3):466-477.(ZHONG S J,MAO X Y,HUANG X J.Water resources exploitation and utilization of international rivers in China from the perspective of geo-security[J].World Regional Studies,2022,31(3):466-477.(in Chinese))

    [12]許紫月,馬凱,袁旭,等.跨境流域水-能源-糧食關聯(lián)關系研究進展與展望[J].地理科學,2023,43(8):1442-1450.(XU Z Y,MA K,YUAN X,et al.Progress and prospect of water-energy-food nexus in the transboundary river basins[J].Scientia Geographica Sinica,2023,43(8):1442-1450.(in Chinese))

    [13]李楊,王婕,唐雄朋,等.變化環(huán)境下瀾滄江-湄公河流域徑流演變及年內分配特征[J].水土保持研究,2021,28(3):141-148.(LI Y,WANG J,TANG X P,et al.Characteristics of runoff variation and seasonal pattern of Lancang-Mekong River in the context of environment change[J].Research of Soil and Water Conservation,2021,28(3):141-148.(in Chinese))

    [14]李昱,席佳,張弛,等.氣候變化對瀾湄流域氣象水文干旱時空特性的影響[J].水科學進展,2021,32(4):508-519.(LI Y,XI J,ZHANG C,et al.Impact of climate change on the spatio-temporal characteristics of meteorological and hydrological drought over the Lancang-Mekong River basin[J].Advances in Water Science,2021,32(4):508-519.(in Chinese))

    [15]談曉珊,王婕,唐雄朋,等.1960—2012年瀾滄江-湄公河流域氣候變化趨勢及不同區(qū)間的徑流響應[J].水資源與水工程學報,2020,31(4):1-8.(TAN X S,WANG J,TANG X P,et al.Variation trends of climate change and hydrological responses of individual regions in Lancang-Mekong River basin from 1960—2012[J].Journal of Water Resources and Water Engineering,2020,31(4):1-8.(in Chinese))

    [16]孫周亮,劉艷麗,劉冀,等.基于博弈論組合賦權法的瀾滄江-湄公河水量分配[J].水資源與水工程學報,2020,31(1):1-5.(SUN Z L,LIU Y L,LIU J,et al.Water allocation of Lancang-Mekong River based on game theory combination weighting method[J].Journal of Water Resources and Water Engineering,2020,31(1):1-5.(in Chinese))

    [17]施國良,薛航.基于XGBoost算法的跨境河流水資源爭端影響因素分析[J].水資源保護,2023,39(6):137-144.(SHI G L,XUE H.XGBoost-based analysis of factors affecting water resources disputes in transboundary rivers[J].Water Resources Protection,2023,39(6):137-144.(in Chinese))

    [18]宋曉猛,張建云,占車生,等.氣候變化和人類活動對水文循環(huán)影響研究進展[J].水利學報,2013,44(7):779-790.(SONG X M,ZHANG J Y,ZHAN C S,et al.Review for impacts of climate change and human activities on water cycle[J].Journal of Hydraulic Engineering,2013,44(7):779-790.(in Chinese)

    [19]張建云,王國慶.河川徑流變化及歸因定量識別[M].北京:科學出版社,2014.(ZHANG J Y,WANG G Q.Quantitative identification of river runoff change and attribution[M].Beijing:Science Press,2014.(in Chinese))

    [20]劉曉燕,王富貴,楊勝天,等.黃土丘陵溝壑區(qū)水平梯田減沙作用研究[J].水利學報,2014,45(7):793-800.(LIU X Y,WANG F G,YANG S T,et al.Sediment reduction effect of level terrace in the hilly-gully region in the Loess Plateau[J].Journal of Hydraulic Engineering,2014,45(7):793-800.(in Chinese))

    [21]謝平,吳林倩,吳子怡,等.西南諸河地表水資源演變的基因圖譜[J].水科學進展,2022,33(3):416-428.(XIE P,WU L Q,WU Z Y,et al.Gene map of water resources evolution of Southwest rivers in China[J].Advances in Water Science,2022,33(3):416-428.(in Chinese))

    [22]DEY P,MISHRA A.Separating the impacts of climate change and human activities on streamflow:a review of methodologies and critical assumptions[J].Journal of Hydrology,2017,548:278-290.

    [23]BOSCH J M,HEWLETT J D.A review of catchment experiments to determine the effect of vegetation changes on water yield and evapotranspiration[J].Journal of Hydrology,1982,55(1/2/3/4):3-23.

    [24]王國慶,張建云,管曉祥,等.中國主要江河徑流變化成因定量分析[J].水科學進展,2020,31(3):313-323.(WANG G Q,ZHANG J Y,GUAN X X,et al.Quantifying attribution of runoff change for major rivers in China[J].Advances in Water Science,2020,31(3):313-323.(in Chinese))

    [25]R?S?NEN T A,SOMETH P,LAURI H,et al.Observed river discharge changes due to hydropower operations in the Upper Mekong basin[J].Journal of Hydrology,2017,545:28-41.

    [26]FAN H,HE D M.Temperature and precipitation variability and its effects on streamflow in the upstream regions of the Lancang-Mekong and Nu-Salween rivers[J].Journal of Hydrometeorology,2015,16(5):2248-2263.

    [27]WU J W,MIAO C Y,WANG Y M,et al.Contribution analysis of the long-term changes in seasonal runoff on the Loess Plateau,China,using eight Budyko-based methods[J].Journal of Hydrology,2017,545:263-275.

    [28]HAN Z Y,LONG D,F(xiàn)ANG Y,et al.Impacts of climate change and human activities on the flow regime of the dammed Lancang River in Southwest China[J].Journal of Hydrology,2019,570:96-105.

    [29]劉宗瑞,周海煒,趙又霖.國際河流跨境洪災防治的合作特征及演進趨勢:基于防洪合作協(xié)議的內容分析[J].中國人口·資源與環(huán)境,2015,25(12):78-84.(LIU Z R,ZHOU H W,ZHAO Y L.Cooperation characters and evolution trends of international rivers′ cross-border flood prevention based on the content analysis of flood prevention cooperation agreements[J].China Population,Resources and Environment,2015,25(12):78-84.(in Chinese))

    [30]何大明,劉恒,馮彥,等.全球變化下跨境水資源理論與方法研究展望[J].水科學進展,2016,27(6):928-934.(HE D M,LIU H,F(xiàn)ENG Y,et al.Perspective on theories and methods study of transboundary water resources under the global change[J].Advances in Water Science,2016,27(6):928-934.(in Chinese))

    [31]王濤,劉承良,杜德斌.1948—2018年國際河流跨境水沖突的時空演化規(guī)律[J].地理學報,2021,76(7):1792-1809.(WANG T,LIU C L,DU D B.Spatio-temporal dynamics of international freshwater conflict events and relations from 1948 to 2018[J].Acta Geographica Sinica,2021,76(7):1792-1809.(in Chinese))

    [32]張瑋瑋,張眉,吳楊,等.復雜地形下浙江夏季氣候要素空間插值方法評價[J].干旱氣象,2020,38(4):674-682.(ZHANG W W,ZHANG M,WU Y,et al.Evaluation of spatial interpolation method about climatic elements in summer in Zhejiang Province under complex topography[J].Journal of Arid Meteorology,2020,38(4):674-682.(in Chinese))

    [33]管曉祥,張建云,鞠琴,等.多種方法在水文關鍵要素一致性檢驗中的比較[J].華北水利水電大學學報(自然科學版),2018,39(2):51-56.(GUAN X X,ZHANG J Y,JU Q,et al.Comparison of consistency testing for key hydrological elements by using multiple statistical methods[J].Journal of North China University of Water Resources and Electric Power (Natural Science Edition),2018,39(2):51-56.(in Chinese))

    [34]吳雷,許有鵬,王強,等.長三角地區(qū)夏季降水結構演變及其非平穩(wěn)性[J].水科學進展,2022,33(5):730-742.(WU L,XU Y P,WANG Q,et al.Evolution and non-stationary characteristics of summer precipitation structure over the Yangtze River Delta[J].Advances in Water Science,2022,33(5):730-742.(in Chinese))

    [35]趙建華,王國慶,張建云,等.RCCC-WBM水量平衡模型在北方典型流域的適用性研究[J].水文,2018,38(2):25-29,14.(ZHAO J H,WANG G Q,ZHANG J Y,et al.Suitability of RCCC-WBM model to typical river catchments in North China[J].Journal of China Hydrology,2018,38(2):25-29,14.(in Chinese))

    [36]WANG G Q,YANG Q L,BAO Z X.RCCC-WBM model and its application for impact assessment of climate change on water resource variability[M].Macao:Unique Network Printing Fty Ltd,2022.

    [37]鞠琴,高慧濱,王國慶,等.基于能量平衡原理的潛在蒸散發(fā)模型構建[J].水科學進展,2022,33(5):794-804.(JU Q,GAO H B,WANG G Q,et al.Modeling potential evapotranspiration based on energy balance[J].Advances in Water Science,2022,33(5):794-804.(in Chinese))

    [38]NASH J E,SUTCLIFFE J V.River flow forecasting through conceptual models part:Ⅰ:a discussion of principles[J].Journal of Hydrology,1970,10(3):282-290.

    [39]陸穎,何大明,何開為,等.瀾滄江-湄公河水電站單位發(fā)電耗水量分析計算[J].水科學進展,2018,29(3):415-423.(LU Y,HE D M,HE K W,et al.Analysis and calculation of water consumption of hydropower plants in Lancang-Mekong River basin[J].Advances in Water Science,2018,29(3):415-423.(in Chinese))

    Impact of environmental changes on variability of outlet hydrological process of

    the Lancang River basin

    The study is financially supported by the National Natural Science Foundation of China (No.92047301;No.52121006).

    WANG Yueyang1,2,ZHANG Jianyun2,3,NING Zhongrui2,4,JIA Yufan2,4,SUN Gaoxia3,5,WANG Guoqing3,5

    (1. College of Hydrology and Water Resources,Hohai University,Nanjing 210098,China;

    2. Yangtze Institute for Conservation and

    Development,Nanjing 210098,China;

    3. National Key Laboratory of Water Disaster Prevention,Nanjing Hydraulic Research

    Institute,Nanjing 210029,China;

    4. Cooperative Innovation Center for Water Safety and Hydro Science,Hohai University,

    Nanjing 210098,China;

    5. Research Center for Climate Change,Ministry of Water Resources,Nanjing 210029,China)

    Abstract:Environment change alters hydrological rhythms and affects utilization of water resources for a river basin.The Lancang-Mekong River is a transboundary river in Southeast Asia originated in the Southwestern China.It is of significance to scientifically quantifying runoff change attribution of the Langcang River for supporting water benefits share between countries which the Lancang-Mekong River runs across.Based on the recorded runoff data series from 1961 to 2020 gauged at the Yunjinghong hydrometric station on the Langcang River,the variability characteristics of runoff series and the impact of climate change and human activities on runoff change as well as its seasonal distribution were investigated by using statistical methods and hydrological simulation methods.The results show that:① The measured annual runoff at Yunjinghong station showed a significant decrease trend from 1961 to 2020,with an increase in runoff during the period from January to May and a decrease in runoff in period from June to December.The abrupt change of the measured annual runoff series was occurred in 2008.② The RCCC-WBM model performs well for monthly natural discharge simulation in the baseline period of 1961—2008 with relative error of less than 1.0%.The simulated natural annual runoff in the human-disturbance period of 2009—2020 is about 371.8 mm,which is 6.8% less than that in baseline period.The reduction is induced by climate change.③ Reservoirs construction and operation is the utmost human activity in the Langcang River.Human activities are the main reason of runoff decrease in the reservoir impoundment period of 2009—2013,while climate change is principal driver of runoff reduction in in the reservoir operation period of 2014—2020,which contribution is about 61.7% of total runoff reduction.④ The reservoirs operation makes runoff decrease by 44.3% in flood season from June to November and increase by 134.3% in dry season from December to May for period of 2014—2020.The regulation function of reservoirs on the Langcang River for decreasing flood discharge and increasing low flow has been guaranteeing security of water supply and flood control for the middle and lower reaches of the Mekong River.

    Key words:climate change;human activities;runoff;RCCC-WBM model;attribution analysis on runoff change;Lancang River basin

    猜你喜歡
    人類活動徑流量氣候變化
    《應對氣候變化報告(2022)》發(fā)布
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    水文比擬法在計算河川徑流量時的修正
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    人為因素主導下海岸帶生態(tài)系統(tǒng)變遷探討
    綠色科技(2015年6期)2015-08-05 18:01:22
    山地表生地球化學研究現(xiàn)狀與展望
    應對氣候變化需要打通“網關”
    太陽能(2015年7期)2015-04-12 06:49:50
    SCS模型在紅壤土坡地降雨徑流量估算中的應用
    亚洲精品456在线播放app | 国产精品久久久久久久久免 | av在线天堂中文字幕| 亚洲美女黄片视频| 国产精品99久久99久久久不卡| 757午夜福利合集在线观看| 久久久久免费精品人妻一区二区| 看黄色毛片网站| 国产伦在线观看视频一区| 伊人久久精品亚洲午夜| 婷婷亚洲欧美| 国产精品久久电影中文字幕| 亚洲熟妇熟女久久| 久久久久久久久大av| 久久香蕉国产精品| 亚洲av成人av| 少妇的逼好多水| 一二三四社区在线视频社区8| 91麻豆精品激情在线观看国产| www.色视频.com| 3wmmmm亚洲av在线观看| 一本综合久久免费| 亚洲av成人av| 波多野结衣高清无吗| 国产精品99久久99久久久不卡| 免费观看精品视频网站| 亚洲天堂国产精品一区在线| 最新在线观看一区二区三区| e午夜精品久久久久久久| 亚洲人成电影免费在线| 久久天躁狠狠躁夜夜2o2o| 日韩av在线大香蕉| 亚洲黑人精品在线| 午夜福利在线观看吧| 高清毛片免费观看视频网站| 久久性视频一级片| 精品欧美国产一区二区三| 久久性视频一级片| 天美传媒精品一区二区| 亚洲乱码一区二区免费版| 亚洲成人久久性| 婷婷精品国产亚洲av| 内射极品少妇av片p| 久久久久免费精品人妻一区二区| 成人精品一区二区免费| 午夜免费男女啪啪视频观看 | 国产aⅴ精品一区二区三区波| e午夜精品久久久久久久| 中国美女看黄片| 极品教师在线免费播放| 丰满乱子伦码专区| 成年女人永久免费观看视频| 国产精品免费一区二区三区在线| 久久这里只有精品中国| 最新中文字幕久久久久| 精品国产美女av久久久久小说| 床上黄色一级片| www.www免费av| 狂野欧美激情性xxxx| 亚洲av电影不卡..在线观看| 欧美乱妇无乱码| 亚洲精品粉嫩美女一区| 他把我摸到了高潮在线观看| 搡女人真爽免费视频火全软件 | 精品人妻一区二区三区麻豆 | 麻豆成人午夜福利视频| 午夜影院日韩av| 亚洲成人精品中文字幕电影| 国内揄拍国产精品人妻在线| 亚洲精品日韩av片在线观看 | 桃红色精品国产亚洲av| 亚洲国产欧洲综合997久久,| 国产精品香港三级国产av潘金莲| 国产老妇女一区| 熟女少妇亚洲综合色aaa.| 国产69精品久久久久777片| 少妇人妻精品综合一区二区 | 亚洲avbb在线观看| 成人特级av手机在线观看| 亚洲男人的天堂狠狠| 在线免费观看的www视频| 婷婷六月久久综合丁香| 日韩成人在线观看一区二区三区| 此物有八面人人有两片| 国产高潮美女av| 亚洲精品在线观看二区| 婷婷精品国产亚洲av在线| 国产三级在线视频| 国产午夜精品久久久久久一区二区三区 | 黄色女人牲交| 男人舔女人下体高潮全视频| 伊人久久大香线蕉亚洲五| 麻豆成人av在线观看| 天天一区二区日本电影三级| 在线播放无遮挡| 色综合婷婷激情| 激情在线观看视频在线高清| 搡女人真爽免费视频火全软件 | 国产精华一区二区三区| 亚洲黑人精品在线| 国产亚洲欧美在线一区二区| www.色视频.com| 十八禁网站免费在线| 亚洲七黄色美女视频| 国产黄色小视频在线观看| or卡值多少钱| 国产精品 国内视频| 精品99又大又爽又粗少妇毛片 | 亚洲人成网站高清观看| 亚洲精品在线美女| 日韩 欧美 亚洲 中文字幕| 五月玫瑰六月丁香| 亚洲精品美女久久久久99蜜臀| 亚洲成人精品中文字幕电影| 日韩精品青青久久久久久| 久久国产乱子伦精品免费另类| 一进一出好大好爽视频| 哪里可以看免费的av片| 亚洲,欧美精品.| 九色成人免费人妻av| 国产精品爽爽va在线观看网站| 久久精品91无色码中文字幕| 夜夜夜夜夜久久久久| 国产精品 欧美亚洲| 99久久成人亚洲精品观看| 国产69精品久久久久777片| 一级黄片播放器| 国产视频内射| 中文字幕精品亚洲无线码一区| 操出白浆在线播放| 精品日产1卡2卡| 国产精品亚洲美女久久久| 午夜福利18| 久久久精品欧美日韩精品| 久久性视频一级片| 亚洲男人的天堂狠狠| 精品久久久久久久久久免费视频| 国内精品久久久久精免费| 99热这里只有精品一区| 中出人妻视频一区二区| 日本一本二区三区精品| 久久精品国产亚洲av涩爱 | 亚洲欧美日韩卡通动漫| 亚洲男人的天堂狠狠| 三级男女做爰猛烈吃奶摸视频| 一a级毛片在线观看| 久久久久久大精品| 69人妻影院| 亚洲av熟女| 久久6这里有精品| 一个人免费在线观看电影| 免费搜索国产男女视频| 成年人黄色毛片网站| 一本一本综合久久| 99久久精品国产亚洲精品| 亚洲无线观看免费| av国产免费在线观看| 久久亚洲真实| 可以在线观看毛片的网站| 日本五十路高清| 91在线精品国自产拍蜜月 | eeuss影院久久| 露出奶头的视频| 日本一二三区视频观看| 午夜老司机福利剧场| 亚洲av熟女| 国产精品久久视频播放| 在线观看午夜福利视频| 麻豆成人午夜福利视频| 亚洲精品456在线播放app | 99在线视频只有这里精品首页| 国产欧美日韩精品一区二区| 亚洲乱码一区二区免费版| 久久久国产精品麻豆| 最近最新免费中文字幕在线| 国产亚洲欧美98| 午夜精品久久久久久毛片777| 男插女下体视频免费在线播放| 亚洲五月婷婷丁香| 黄色片一级片一级黄色片| 国产v大片淫在线免费观看| 国内久久婷婷六月综合欲色啪| 久久精品国产综合久久久| 亚洲成人免费电影在线观看| 琪琪午夜伦伦电影理论片6080| 最新在线观看一区二区三区| 成年女人看的毛片在线观看| 俄罗斯特黄特色一大片| 欧美大码av| 激情在线观看视频在线高清| 精品久久久久久,| 九九久久精品国产亚洲av麻豆| 美女黄网站色视频| 久久久久久久午夜电影| 国产精品精品国产色婷婷| 性色avwww在线观看| 亚洲第一欧美日韩一区二区三区| 人妻丰满熟妇av一区二区三区| 熟女电影av网| 黄色视频,在线免费观看| 内地一区二区视频在线| av在线蜜桃| 亚洲激情在线av| 一夜夜www| 欧美成人一区二区免费高清观看| a在线观看视频网站| 国产美女午夜福利| 黄色丝袜av网址大全| av专区在线播放| 免费看光身美女| 又粗又爽又猛毛片免费看| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色| 91在线观看av| 成人特级黄色片久久久久久久| 亚洲电影在线观看av| 身体一侧抽搐| 亚洲成人精品中文字幕电影| 久久香蕉国产精品| 成人无遮挡网站| 亚洲成a人片在线一区二区| 成人精品一区二区免费| www.色视频.com| 黄色片一级片一级黄色片| 国产精品一区二区三区四区免费观看 | 激情在线观看视频在线高清| 精品人妻偷拍中文字幕| 国产高清激情床上av| 久久精品综合一区二区三区| 国产一区二区激情短视频| 亚洲精品粉嫩美女一区| 欧美成人性av电影在线观看| 美女cb高潮喷水在线观看| 亚洲成a人片在线一区二区| 久久人妻av系列| 韩国av一区二区三区四区| 午夜福利高清视频| av视频在线观看入口| 日韩大尺度精品在线看网址| 久久精品国产自在天天线| 欧美乱妇无乱码| 免费观看精品视频网站| 精品国产超薄肉色丝袜足j| 好男人电影高清在线观看| 国产成人系列免费观看| АⅤ资源中文在线天堂| 免费看美女性在线毛片视频| av女优亚洲男人天堂| 天堂影院成人在线观看| 伊人久久精品亚洲午夜| 丁香欧美五月| 欧美日韩福利视频一区二区| 亚洲人成网站在线播放欧美日韩| 天美传媒精品一区二区| 小蜜桃在线观看免费完整版高清| 亚洲国产精品sss在线观看| 成人午夜高清在线视频| x7x7x7水蜜桃| 亚洲不卡免费看| a级一级毛片免费在线观看| 特级一级黄色大片| 好男人电影高清在线观看| 久久香蕉国产精品| 母亲3免费完整高清在线观看| 久久久国产成人精品二区| 少妇高潮的动态图| 国产精品自产拍在线观看55亚洲| 欧美最黄视频在线播放免费| 久久中文看片网| 日本黄色片子视频| 国产探花极品一区二区| 午夜福利18| 成年女人永久免费观看视频| 日韩高清综合在线| 亚洲avbb在线观看| 国产熟女xx| 国产午夜精品久久久久久一区二区三区 | 最近最新中文字幕大全电影3| av中文乱码字幕在线| 青草久久国产| 欧美日韩综合久久久久久 | 9191精品国产免费久久| 国内少妇人妻偷人精品xxx网站| 中文字幕人成人乱码亚洲影| 欧美黑人欧美精品刺激| 搡老妇女老女人老熟妇| 97超级碰碰碰精品色视频在线观看| 日韩欧美在线乱码| 精品电影一区二区在线| 国产高清视频在线观看网站| 国产成人aa在线观看| 三级男女做爰猛烈吃奶摸视频| 叶爱在线成人免费视频播放| 18+在线观看网站| 国产97色在线日韩免费| 伊人久久大香线蕉亚洲五| 一本久久中文字幕| 黄色成人免费大全| 69人妻影院| 天堂网av新在线| 91字幕亚洲| 国产伦人伦偷精品视频| 国产美女午夜福利| 国产色婷婷99| 国产精品国产高清国产av| 好男人电影高清在线观看| 国产黄a三级三级三级人| 一级作爱视频免费观看| 国产高清videossex| av天堂中文字幕网| 国产精品,欧美在线| 欧美黄色淫秽网站| 精品人妻1区二区| 午夜老司机福利剧场| 香蕉久久夜色| 婷婷六月久久综合丁香| 中国美女看黄片| 亚洲男人的天堂狠狠| 香蕉久久夜色| 熟女电影av网| 国内揄拍国产精品人妻在线| 国产伦人伦偷精品视频| 又爽又黄无遮挡网站| 成人高潮视频无遮挡免费网站| 日本 av在线| 中亚洲国语对白在线视频| 国产成+人综合+亚洲专区| 国产一区二区激情短视频| 黑人欧美特级aaaaaa片| 国产欧美日韩一区二区三| 久久草成人影院| 九九在线视频观看精品| 波多野结衣高清作品| 嫁个100分男人电影在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品999在线| 女生性感内裤真人,穿戴方法视频| 国产成人aa在线观看| 午夜日韩欧美国产| av欧美777| 女人被狂操c到高潮| 精品一区二区三区人妻视频| 午夜免费激情av| 国产精品98久久久久久宅男小说| 一个人免费在线观看电影| 黑人欧美特级aaaaaa片| av中文乱码字幕在线| 国产一区二区三区在线臀色熟女| 精品熟女少妇八av免费久了| 人妻丰满熟妇av一区二区三区| 12—13女人毛片做爰片一| АⅤ资源中文在线天堂| 99在线视频只有这里精品首页| 久久久久久久久中文| 每晚都被弄得嗷嗷叫到高潮| 嫩草影院精品99| 少妇的丰满在线观看| 天堂av国产一区二区熟女人妻| 91av网一区二区| aaaaa片日本免费| 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 少妇的逼水好多| 亚洲自拍偷在线| 18禁在线播放成人免费| 51国产日韩欧美| 在线播放无遮挡| 三级国产精品欧美在线观看| 国产亚洲欧美在线一区二区| 精品午夜福利视频在线观看一区| 欧美区成人在线视频| 国产91精品成人一区二区三区| av女优亚洲男人天堂| 欧美激情久久久久久爽电影| 舔av片在线| 中文字幕av在线有码专区| 亚洲内射少妇av| av黄色大香蕉| 国产av在哪里看| 母亲3免费完整高清在线观看| 精品人妻一区二区三区麻豆 | 欧美日韩乱码在线| 国产成人aa在线观看| 久久久久久久久中文| 午夜免费成人在线视频| 国产精品 国内视频| 午夜免费男女啪啪视频观看 | 国产色爽女视频免费观看| 波多野结衣高清无吗| 免费人成视频x8x8入口观看| 久久国产精品人妻蜜桃| 国产极品精品免费视频能看的| 国产免费av片在线观看野外av| 久久久久国内视频| 国产在视频线在精品| 又粗又爽又猛毛片免费看| 91久久精品国产一区二区成人 | 亚洲中文日韩欧美视频| 欧美av亚洲av综合av国产av| 亚洲五月天丁香| 免费看光身美女| 69人妻影院| 又黄又爽又免费观看的视频| 色精品久久人妻99蜜桃| av天堂中文字幕网| 久久性视频一级片| 99久久久亚洲精品蜜臀av| 亚洲最大成人手机在线| 三级国产精品欧美在线观看| 18禁黄网站禁片午夜丰满| 国产99白浆流出| 欧美激情在线99| 波多野结衣巨乳人妻| 日本与韩国留学比较| 天美传媒精品一区二区| 亚洲不卡免费看| 精品国产亚洲在线| 黄色视频,在线免费观看| 亚洲不卡免费看| 国产爱豆传媒在线观看| 国产伦人伦偷精品视频| 国产久久久一区二区三区| 欧美日本亚洲视频在线播放| 狂野欧美激情性xxxx| 成人三级黄色视频| 亚洲不卡免费看| 一进一出好大好爽视频| 人妻丰满熟妇av一区二区三区| 天堂√8在线中文| 制服人妻中文乱码| 3wmmmm亚洲av在线观看| 国产美女午夜福利| 特级一级黄色大片| 99久久精品热视频| tocl精华| 亚洲中文日韩欧美视频| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 一本一本综合久久| 少妇高潮的动态图| 国产一区二区三区在线臀色熟女| 最近最新中文字幕大全电影3| 色在线成人网| 午夜福利视频1000在线观看| 亚洲国产精品合色在线| 亚洲精品美女久久久久99蜜臀| 高清在线国产一区| 十八禁网站免费在线| 精品久久久久久久人妻蜜臀av| 最近视频中文字幕2019在线8| 国产精品99久久久久久久久| 少妇裸体淫交视频免费看高清| 国产高清videossex| 精品久久久久久,| 欧美成人a在线观看| 国产精品野战在线观看| www日本黄色视频网| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 久久人人精品亚洲av| 国产亚洲精品久久久com| 久久久久国产精品人妻aⅴ院| 亚洲,欧美精品.| 熟女人妻精品中文字幕| 欧美日韩瑟瑟在线播放| 天天添夜夜摸| 国产精品一区二区三区四区免费观看 | 91久久精品国产一区二区成人 | 国产视频内射| 99热精品在线国产| 亚洲天堂国产精品一区在线| 最新中文字幕久久久久| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利视频1000在线观看| 久9热在线精品视频| 欧美日韩福利视频一区二区| 亚洲第一电影网av| 欧美精品啪啪一区二区三区| av天堂在线播放| 国产精品精品国产色婷婷| 最近在线观看免费完整版| 成人永久免费在线观看视频| 国产三级在线视频| 久久精品91无色码中文字幕| 欧美一级a爱片免费观看看| 国产高潮美女av| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 国产伦在线观看视频一区| 亚洲精品日韩av片在线观看 | 国产精品三级大全| 少妇人妻精品综合一区二区 | 欧美日韩综合久久久久久 | 51午夜福利影视在线观看| 老司机福利观看| 国产精品98久久久久久宅男小说| 免费高清视频大片| 青草久久国产| 成人18禁在线播放| 色尼玛亚洲综合影院| 男插女下体视频免费在线播放| 18禁在线播放成人免费| 国产高清视频在线观看网站| 美女黄网站色视频| av国产免费在线观看| 婷婷精品国产亚洲av| 亚洲精品日韩av片在线观看 | 99久久精品国产亚洲精品| 九九在线视频观看精品| 国产一区二区在线av高清观看| 亚洲精品456在线播放app | 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美激情综合另类| 香蕉久久夜色| 久久精品国产亚洲av香蕉五月| 欧美日韩一级在线毛片| 又黄又爽又免费观看的视频| 亚洲精品色激情综合| 日本黄色片子视频| 老司机在亚洲福利影院| 中文字幕人妻丝袜一区二区| 国产日本99.免费观看| 精品不卡国产一区二区三区| 九色成人免费人妻av| 黄色片一级片一级黄色片| 看黄色毛片网站| 日韩欧美一区二区三区在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久亚洲精品不卡| 桃色一区二区三区在线观看| 亚洲 欧美 日韩 在线 免费| 国产亚洲欧美在线一区二区| 欧美一级毛片孕妇| 亚洲电影在线观看av| av在线天堂中文字幕| 国产又黄又爽又无遮挡在线| 亚洲av成人不卡在线观看播放网| 日韩高清综合在线| 国产亚洲精品一区二区www| 色吧在线观看| 精品日产1卡2卡| 精品一区二区三区视频在线观看免费| 国产黄色小视频在线观看| 啦啦啦观看免费观看视频高清| 91av网一区二区| 成人av在线播放网站| 老司机在亚洲福利影院| 亚洲性夜色夜夜综合| 亚洲avbb在线观看| 婷婷精品国产亚洲av在线| 日本一二三区视频观看| 91在线观看av| 波多野结衣高清作品| 欧美极品一区二区三区四区| 不卡一级毛片| 免费大片18禁| 亚洲成av人片在线播放无| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式 | 久久精品91无色码中文字幕| 少妇熟女aⅴ在线视频| 亚洲18禁久久av| 99久久成人亚洲精品观看| 色在线成人网| 国产中年淑女户外野战色| 日本成人三级电影网站| 一夜夜www| 真人做人爱边吃奶动态| 91麻豆精品激情在线观看国产| 免费看日本二区| 一区二区三区高清视频在线| а√天堂www在线а√下载| 一级毛片高清免费大全| 在线免费观看的www视频| 国产精品久久久久久久久免 | 校园春色视频在线观看| av天堂在线播放| 三级毛片av免费| 国产99白浆流出| 十八禁人妻一区二区| 久久久久国产精品人妻aⅴ院| 欧美性猛交╳xxx乱大交人| 午夜精品一区二区三区免费看| 精品国产美女av久久久久小说| 国产单亲对白刺激| 久久久久亚洲av毛片大全| 一级作爱视频免费观看| 国产一级毛片七仙女欲春2| a在线观看视频网站| 一个人看视频在线观看www免费 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久久末码| 好看av亚洲va欧美ⅴa在| 亚洲熟妇熟女久久| 国产在线精品亚洲第一网站| 最新中文字幕久久久久| 亚洲五月婷婷丁香| 国产伦在线观看视频一区| 欧美日韩中文字幕国产精品一区二区三区| 久久草成人影院| 一级黄色大片毛片| 国产伦人伦偷精品视频| 国产激情偷乱视频一区二区| 在线免费观看的www视频| 国产综合懂色| 国产精品久久久人人做人人爽| 国产成人欧美在线观看| 中文字幕av成人在线电影| 免费看十八禁软件| 亚洲色图av天堂| 成人午夜高清在线视频| 国产熟女xx| 噜噜噜噜噜久久久久久91| 18+在线观看网站| 国产亚洲欧美在线一区二区|