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

    氣候變化影響下的流域水循環(huán):以長(zhǎng)江和黃河典型流域?yàn)槔?/h1>
    2016-11-16 03:31:04徐宗學(xué)劉瀏
    關(guān)鍵詞:太湖流域水循環(huán)徑流

    徐宗學(xué)劉瀏

    (1 北京師范大學(xué)水科學(xué)研究院,全球變化研究協(xié)同創(chuàng)新中心,北京 100875;2 中國(guó)農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,中國(guó)農(nóng)業(yè)水問題研究中心,北京 100083)

    氣候變化影響下的流域水循環(huán):以長(zhǎng)江和黃河典型流域?yàn)槔?/p>

    徐宗學(xué)1劉瀏2

    (1 北京師范大學(xué)水科學(xué)研究院,全球變化研究協(xié)同創(chuàng)新中心,北京 100875;2 中國(guó)農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,中國(guó)農(nóng)業(yè)水問題研究中心,北京 100083)

    在對(duì)IPCC提供的多種大氣環(huán)流模式(GCMs)適用性評(píng)估的基礎(chǔ)上,采用SDSM和ASD統(tǒng)計(jì)降尺度模型生成未來氣候變化情景,驅(qū)動(dòng)分布式水文模型VIC和SWAT,分別對(duì)長(zhǎng)江和黃河典型流域進(jìn)行分布式水文模擬,定量探討氣候變化對(duì)流域水循環(huán)的影響。結(jié)果表明,適應(yīng)性評(píng)估可以有效降低GCMs選擇的不確定性,統(tǒng)計(jì)降尺度方法能夠明顯改善全球氣候模式降水和氣溫輸出結(jié)果。與基準(zhǔn)期(1961—1990年)相比,未來時(shí)期(2046—2065年和2081—2100年)長(zhǎng)江下游太湖流域徑流量呈現(xiàn)微弱減少趨勢(shì),但汛期東南部徑流量增加顯著,而黃河上游流域徑流量則呈下降趨勢(shì)。研究結(jié)果可為開展我國(guó)各大流域適應(yīng)氣候變化研究提供一定的參考依據(jù)。

    氣候變化,水循環(huán),降尺度,水文模型,大氣環(huán)流模式

    0 引言

    政府間氣候變化委員會(huì)(IPCC)第五次評(píng)估報(bào)告[1]指出:全球氣候變化已是不爭(zhēng)的事實(shí),將對(duì)全球和區(qū)域水資源安全構(gòu)成嚴(yán)重威脅。我國(guó)氣候變化趨勢(shì)與全球基本一致[2-3]。氣候變化對(duì)人類、生態(tài)、環(huán)境產(chǎn)生的多尺度、全方位、多層次的影響[4-11],已使其成為威脅全球可持續(xù)性發(fā)展的危險(xiǎn)因素之一。水循環(huán)作為人類—生態(tài)—環(huán)境的紐帶,預(yù)估其對(duì)未來氣候變化的響應(yīng),已成為各國(guó)科學(xué)家、公眾和決策者共同關(guān)心的課題。

    在全球變暖背景下,中國(guó)的高溫、低溫、強(qiáng)降水、干旱等極端天氣氣候事件的頻率和強(qiáng)度存在不同的變化趨勢(shì),并有區(qū)域差異[12-13]。在過去100年中,我國(guó)主要流域水文氣候要素發(fā)生了明顯的變化,主要表現(xiàn)為大部分地區(qū)氣溫顯著升高,近50 年?yáng)|部季風(fēng)區(qū)域降水存在“南澇北旱”的趨勢(shì)變化[14-17],北方多數(shù)外流河流域徑流減少,加劇了北方水資源的供需矛盾和南方防洪抗災(zāi)的壓力[18-20];全國(guó)大部分地區(qū)潛在蒸發(fā)能力(水面蒸發(fā))明顯下降[21-24];未來的氣候變化可能對(duì)中國(guó)水資源產(chǎn)生較大影響[25-28],北方部分江河徑流量可能減少、南方江河徑流量可能增加[29],各流域年平均蒸發(fā)量可能增大[30],旱澇等災(zāi)害的出現(xiàn)頻率可能增加,進(jìn)一步加劇水資源的脆弱性與供需矛盾[29]。

    針對(duì)我國(guó)大江大河流域水循環(huán)及水資源格局變化以及南方洪澇災(zāi)害與北方水資源短缺等重大水問題,開展氣候變化背景下陸地水循環(huán)響應(yīng)機(jī)理研究,對(duì)于科學(xué)認(rèn)識(shí)全球氣候變化背景下我國(guó)陸地水循環(huán)時(shí)空演變規(guī)律,評(píng)估氣候變化對(duì)流域水循環(huán)的影響,具有重要科學(xué)意義與應(yīng)用價(jià)值。

    本文選取長(zhǎng)江和黃河典型流域分別作為我國(guó)南方、北方流域的代表,采用秩打分法對(duì)IPCC提供的多種大氣環(huán)流模式(GCMs)進(jìn)行適應(yīng)性評(píng)估,并對(duì)篩選出的GCMs輸出結(jié)果進(jìn)行降尺度處理,分別構(gòu)建適合長(zhǎng)江流域和黃河流域的未來氣候變化情景,進(jìn)而驅(qū)動(dòng)分布式水文模型,模擬研究區(qū)水循環(huán)要素的時(shí)空演變規(guī)律,評(píng)估流域水循環(huán)對(duì)未來氣候變化的響應(yīng)。

    1 GCM適應(yīng)性評(píng)估

    1.1 評(píng)估方法

    將GCMs輸出統(tǒng)計(jì)特征值與實(shí)測(cè)數(shù)據(jù)統(tǒng)計(jì)特征值的擬合程度作為目標(biāo)函數(shù),根據(jù)各個(gè)目標(biāo)函數(shù)的表現(xiàn)進(jìn)行評(píng)分,進(jìn)而綜合評(píng)價(jià)氣候模式的整體表現(xiàn)。

    以均值、變異系數(shù)(=標(biāo)準(zhǔn)差/均值)和歸一化的均方根誤差(NRMSE)表征GCMs輸出與實(shí)測(cè)變量均值和方差的吻合程度。NRMSE定義為均方根誤差與相應(yīng)標(biāo)準(zhǔn)差的比值。

    以時(shí)間序列(多年平均月序列)的Pearson相關(guān)系數(shù)和空間序列(各個(gè)站點(diǎn)氣候要素均值)的Pearson相關(guān)系數(shù)表征GCMs模擬值與實(shí)測(cè)值之間的吻合程度,可以分別用于評(píng)估模式模擬變量的年內(nèi)變化能力和空間擬合程度。

    利用Mann-Kendall方法檢測(cè)各個(gè)變量的長(zhǎng)期變化趨勢(shì)(MK Zc)及其變化幅度(MK Slope),如果GCM輸出變量的長(zhǎng)期變化趨勢(shì)(MK Zc)與實(shí)測(cè)值一致的區(qū)域沒有達(dá)到95%,則趨勢(shì)變化幅度(MK Slope)不參與評(píng)分。

    利用經(jīng)驗(yàn)正交函數(shù)(EOF)反映變量時(shí)空綜合變異的特點(diǎn),前期分析得知近20個(gè)變量的EOF第一和第二特征向量可以解釋其絕大部分的方差,因此,選用前兩個(gè)EOF的特征向量進(jìn)行秩評(píng)分。

    利用概率密度函數(shù)的兩個(gè)統(tǒng)計(jì)量BS(Brier score)和Sscore(Skill score)來評(píng)價(jià)GCMs對(duì)概率密度函數(shù)的模擬效果,BS是概率預(yù)測(cè)的均方差,Sscore用以描述模擬概率分布與實(shí)測(cè)值的重疊程度。

    以上包括了均值、變異系數(shù)、NRMSE、時(shí)間序列相關(guān)系數(shù)、空間序列相關(guān)系數(shù)、趨勢(shì)分析秩統(tǒng)計(jì)量、趨勢(shì)分析變化幅度、EOF第一特征向量、EOF第二特征向量、BS和Sscore共11個(gè)統(tǒng)計(jì)量,每個(gè)統(tǒng)計(jì)量的表現(xiàn)均作為目標(biāo)函數(shù)之一。對(duì)每個(gè)目標(biāo)函數(shù)分別進(jìn)行秩評(píng)分,根據(jù)GCMs表現(xiàn)的差異賦以0~10不同的值,計(jì)算公式如下:

    其中,xi為模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)應(yīng)的待評(píng)估統(tǒng)計(jì)量之間的相對(duì)誤差,誤差小則給的分?jǐn)?shù)也小。對(duì)于特定GCM輸出變量,最后總評(píng)分由所有統(tǒng)計(jì)量對(duì)應(yīng)的分?jǐn)?shù)求平均值得到,模擬效果越好,評(píng)分越小。該評(píng)分結(jié)果代表了GCMs輸出結(jié)果與“實(shí)測(cè)”值之間的統(tǒng)計(jì)特征擬合程度,適用于不同GCMs之間的對(duì)比分析,并不代表某一GCM的實(shí)際模擬精度。

    1.2數(shù)據(jù)

    1)GCMs輸出數(shù)據(jù)

    為盡可能多地納入全球不同國(guó)家的GCM,以及保證各GCM輸出數(shù)據(jù)的完整性和可靠度,本文選用了IPCC第四次評(píng)估報(bào)告推薦的GCMs(表1),要素包括500,700,850hPa等壓面的平均氣溫、相對(duì)濕度、緯向風(fēng)速、經(jīng)向風(fēng)速和位勢(shì)高度共15個(gè)氣象變量,以及氣溫和降水兩個(gè)地面氣候要素,時(shí)間尺度為月,所有數(shù)據(jù)均來自于IPCC數(shù)據(jù)中心,各GCM驅(qū)動(dòng)場(chǎng)等詳細(xì)信息可參見網(wǎng)站http://ipcc-ddc.cru.uea.ac.uk。所有GCMs輸出數(shù)據(jù)統(tǒng)一插值到2.5°×2.5°的網(wǎng)格分辨率,數(shù)據(jù)序列為1960—1999/2000年。

    2)ERA-40/NCEP再分析資料

    選取ERA-40和NCEP再分析資料作為各等壓面氣象變量的“實(shí)際”值,用來評(píng)估GCMs輸出結(jié)果的適用性。

    3)地面觀測(cè)資料

    利用中國(guó)氣象局提供的地表月平均氣溫和月降水格點(diǎn)數(shù)據(jù)(http://ncc.cma.gov.cn),重采樣成2.5°×2.5°網(wǎng)格分辨率,用來評(píng)估GCMs地面輸出成果。

    1.3長(zhǎng)江流域評(píng)估結(jié)果

    本研究選取的長(zhǎng)江流域GCM網(wǎng)格如圖1所示。綜合評(píng)價(jià)結(jié)果表明, FGOALS、ECHAM4、ECHAM5、HadCM3、HadGEM1和MRI在長(zhǎng)江流域的表現(xiàn)較其他氣候模式要好一些(圖2)。

    表1 氣候模式簡(jiǎn)介Table 1 Introduction of GCMs

    張?jiān)鲂诺萚31]在研究長(zhǎng)江流域水汽收支的時(shí)空變化與環(huán)流特征時(shí)指出,源自孟加拉灣的經(jīng)向水汽輸送對(duì)流域相對(duì)濕度起到重要影響,大氣緯向風(fēng)速和經(jīng)向風(fēng)速的模擬效果很大程度上影響流域內(nèi)相對(duì)濕度的時(shí)空分布。本研究中,相對(duì)濕度模擬較好的幾個(gè)GCM(ECHAM4.、FGOALS、HadCM3和HadGEM1)對(duì)緯向風(fēng)速和經(jīng)向風(fēng)速的模擬同樣都較好,這與上述得出的結(jié)論一致。而對(duì)氣溫模擬明顯要好于相對(duì)濕度等其他變量,多數(shù)研究也都表明GCM對(duì)氣溫的模擬較好。徐影等[32]指出,ECHAM4和HadCM2模式對(duì)中國(guó)地區(qū)氣候模擬最好,這兩種模式及其在其基礎(chǔ)上發(fā)展的模式在本研究中也都表現(xiàn)出良好的模擬能力。但Maxino等[33]評(píng)估澳大利亞Murray-Darling流域GCM適用性時(shí),發(fā)現(xiàn)絕大多數(shù)氣候模式的氣溫模擬值較實(shí)測(cè)值偏大。曹穎等[34]評(píng)價(jià)黃河流域GCM時(shí),氣溫模擬值均偏小。

    圖1 長(zhǎng)江流域及選取網(wǎng)格圖Fig. 1 Location of the Yangtze River Basin and selected GCM grids

    本文選取太湖流域作為典型區(qū)域開展氣候變化對(duì)長(zhǎng)江流域水循環(huán)的影響研究,并開展GCMs的適應(yīng)性評(píng)估。圖3為RS評(píng)分排序,評(píng)分分?jǐn)?shù)越低,表明模式模擬氣候的能力越強(qiáng)??梢钥闯?,BCCR模式表現(xiàn)最好。因此,選取BCCR模式進(jìn)行氣候變化影響評(píng)估,其模擬的日尺度氣候數(shù)據(jù)是在A1B情景下輸出的,分為兩個(gè)部分:當(dāng)前時(shí)期1961—2000 年(20C3M情景)和未來時(shí)期(2046—2065年和2081—2100年)。該情景是一種在經(jīng)濟(jì)發(fā)展與溫室氣體排放之間保持平衡的發(fā)展情景,較其他情景可能更為符合太湖流域未來水—經(jīng)濟(jì)—環(huán)境協(xié)調(diào)發(fā)展規(guī)劃。

    圖2 長(zhǎng)江流域GCMs綜合評(píng)分圖Fig. 2 Rank Score (RS) values of GCMs for the Yangtze River Basin

    圖3 太湖流域GCMs綜合評(píng)分Fig. 3 RS values of GCMs for the Taihu Basin

    1.4黃河流域評(píng)估結(jié)果

    本文選用了覆蓋黃河流域的25個(gè)網(wǎng)格GCMs數(shù)據(jù)(圖4),并將所有氣候模式數(shù)據(jù)以線性插值的方法統(tǒng)一插值到2.5°×2.5°分辨率的網(wǎng)格上。綜合評(píng)估結(jié)果(圖5)表明,最好的10個(gè)氣候模式為MRI、HadGEM1、INM、CSIRO30、MIROC-M、 HadCM3、BCCR、GFDL20、CGCM47和GFDL21。

    圖4 黃河流域及選取網(wǎng)格圖Fig. 4 Location of the Yellow River Basin and selected GCM grids

    圖5  黃河流域GCMs綜合評(píng)分圖Fig. 5 RS values of GCMs for the Yellow River Basin

    周曉霞等[35]在研究夏季亞洲季風(fēng)區(qū)的水汽輸送及其對(duì)中國(guó)降水的影響時(shí)指出,印度季風(fēng)帶來的緯向水汽輸送和經(jīng)向水汽輸送對(duì)流域相對(duì)濕度起到重要影響,大氣緯向風(fēng)速和經(jīng)向風(fēng)速的模擬效果很大程度上影響流域內(nèi)相對(duì)濕度的時(shí)空分布。本研究中,相對(duì)濕度模擬較好的幾個(gè)氣候模式(MRI、HadGEM1和GFDL20)等對(duì)緯向風(fēng)速和經(jīng)向風(fēng)速的模擬同樣都較好,這與上述結(jié)論是一致的。

    為探討氣候變化影響研究中由多模式及不同情景造成的不確定性,綜合考慮各模式情景數(shù)據(jù)的一致性和完整性,本研究選取MRI、CSIRO30和INM三種在黃河流域表現(xiàn)最優(yōu)的模式輸出的A2和B1情景下的數(shù)據(jù)進(jìn)行對(duì)比研究。

    2 氣候變化對(duì)長(zhǎng)江典型流域水循環(huán)的影響

    2.1未來氣候變化情景構(gòu)建

    本文作者已針對(duì)太湖流域開展了統(tǒng)計(jì)降尺度模型對(duì)比分析研究[36],即采用SDSM和ASD分別構(gòu)建太湖流域未來氣候變化情景,對(duì)比分析表明ASD對(duì)太湖流域未來氣候的模擬優(yōu)于SDSM。因此,本文采取ASD模型構(gòu)建的氣候變化情景開展研究。

    ASD降尺度結(jié)果顯示,與未來氣溫的變化相比(圖6),未來降水的變化更為復(fù)雜。

    2.2長(zhǎng)江典型流域水循環(huán)對(duì)氣候變化的響應(yīng)

    1)分布式水文模型構(gòu)建

    太湖流域地勢(shì)平坦、河網(wǎng)復(fù)雜、圩區(qū)眾多,沒有明確流域出口,并且基于掌握的90m DEM難以提取太湖流域河網(wǎng),難以進(jìn)行全流域產(chǎn)匯流模擬,因此本研究選取基于網(wǎng)格的分布式 VIC-3L模型[37-42],直接獲取太湖流域網(wǎng)格化徑流深。作者已成功構(gòu)建太湖流域VIC模型[43],可供本研究開展氣候變化影響研究。

    2)氣候變化影響定量評(píng)估

    采用ASD模型生成的A1B情景下基準(zhǔn)期以及未來時(shí)期日降水、日最高氣溫和最低氣溫的氣象數(shù)據(jù)系列,通過泰森多邊形法插值到太湖流域5km×5km網(wǎng)格,從而建立氣候強(qiáng)迫數(shù)據(jù);根據(jù)已建立的太湖流域VIC模型土壤和植被參數(shù)文件[43],使VIC模型在太湖流域的1452個(gè)網(wǎng)格上連續(xù)運(yùn)行,獨(dú)立輸出每個(gè)網(wǎng)格不同時(shí)期的日徑流深數(shù)據(jù)系列。

    太湖流域2046—2065年各網(wǎng)格相對(duì)基準(zhǔn)期的多年平均月徑流深變化的情況如圖7所示。由圖可見,太湖流域大部分地區(qū)未來時(shí)期年內(nèi)平均月徑流深呈現(xiàn)減少趨勢(shì),而各月的減少幅度不一。較基準(zhǔn)期流域徑流深減少幅度最大的地區(qū)是浙西區(qū)與杭嘉湖區(qū);而其他地區(qū)相對(duì)較小,尤其是上海地區(qū)在多數(shù)月(1—4和9—12月)都呈現(xiàn)徑流深增加的趨勢(shì)。不過,浙西區(qū)在汛期呈現(xiàn)略微增加的趨勢(shì)。

    圖6 ASD降尺度模型模擬太湖流域未來兩個(gè)時(shí)段的日最高氣溫(a—b)、日最低氣溫(c—d)和降水(e—f)及其變化情況(a、c、e為3個(gè)要素與實(shí)測(cè)值的對(duì)比,b、d、f為年內(nèi)各月要素相對(duì)基準(zhǔn)期的變化)Fig. 6 The ASD model simulated daily maximum air temperature (a-b), minimum air temperature (c-d), precipitation (e-f),and monthly changing patterns for the future two periods (a, c, e: comparisons between simulations and observations, b, d,f: monthly changes for two future periods comparing with that in the baseline period)

    2081—2100年,太湖流域大部分地區(qū)未來時(shí)期年內(nèi)平均月徑流深呈現(xiàn)減少趨勢(shì),而各月的減少幅度不一(圖略)。較基準(zhǔn)期而言,浙西區(qū)的徑流深雖有減少的趨勢(shì),不過相比2046—2065年,減少幅度有所減少。相比之下, 3—5月,流域內(nèi)徑流深減少幅度最大的是陽(yáng)澄淀區(qū);而在冬季杭嘉湖區(qū)的徑流減少幅度較大;與前一時(shí)期相似的是上海地區(qū)在多數(shù)月(1—4和10—12月)都呈現(xiàn)徑流深增加的趨勢(shì)。

    3 氣候變化對(duì)黃河典型流域水循環(huán)的影響

    3.1未來氣候變化情景構(gòu)建

    基于GCMs適應(yīng)性評(píng)估結(jié)果,本研究共選取A2和B1兩種情景下的三種GCM(MRI、CSIRO30、INM)輸出結(jié)果,數(shù)據(jù)按照2.5°×2.5°的網(wǎng)格進(jìn)行插值處理,模擬時(shí)段為基準(zhǔn)期1961—1990年、未來時(shí)期2046—2065年和2081—2100年。采用 SDSM模型進(jìn)行降尺度處理,選取12個(gè)日預(yù)報(bào)因子,覆蓋了黃河流域的25個(gè)網(wǎng)格。觀測(cè)數(shù)據(jù)包含1961—2000年間黃河流域79個(gè)測(cè)站的日降水量、平均氣溫、最高及最低氣溫的日觀測(cè)值。

    圖7 研究區(qū)A1B情景下2046—2065年多年平均月徑流深變化的空間分布Fig. 7 Spatial distribution of monthly runoff changes under A1B scenario for 2046-2065

    圖8顯示了三種GCMs在不同情景(A2、B1)和未來時(shí)期共六種情形下月最高氣溫、月最低氣溫的模擬結(jié)果。可以看出,月最高、最低氣溫在六種情形下的大部分月份里,都顯示出增長(zhǎng)趨勢(shì)。

    不同情景下,降水量的變化差異較大(圖9)??偟膩碚f,CSIRO30模擬的降水量減少幅度是三種模式中最小的,月變化降低不到25%。MRI模擬的降水量變化在三個(gè)模式中增幅最大。

    3.2黃河典型流域水循環(huán)對(duì)氣候變化的響應(yīng)

    1)分布式水文模型構(gòu)建

    本研究選取黃河上游流域(圖10)作為典型區(qū)開展氣候變化對(duì)黃流域水循環(huán)影響研究。

    黃河上游流域地勢(shì)起伏明顯,河網(wǎng)較易提取,流域出口水文站數(shù)據(jù)豐富,因此選用SWAT模型[44-45]進(jìn)行分布式模擬,生成流域出口流量過程。采用模型自帶敏感性分析模塊對(duì)SWAT模型進(jìn)行了參數(shù)敏感性分析,然后基于此對(duì)參數(shù)進(jìn)行率定和驗(yàn)證。以1971—1985年為率定期對(duì)模型參數(shù)進(jìn)行率定,1986—2000年為模型驗(yàn)證期對(duì)模型進(jìn)行驗(yàn)證。

    蘭州水文站在率定期和驗(yàn)證期的徑流模擬值與實(shí)測(cè)值的對(duì)比如圖11,總體來說,模型的評(píng)價(jià)指標(biāo)較為令人滿意。率定期納西效率系數(shù)為0.667,確定性系數(shù)為0.773;驗(yàn)證期納西效率系數(shù)為0.626,確定性系數(shù)為0.709。因此,模型較好模擬了流域內(nèi)的水量平衡。

    圖8 不同情景(A2、B1)月最高氣溫(a,b)、月最低氣溫(c,d)的模擬結(jié)果Fig. 8 Simulated monthly maximum (a, b) /minimum (c, d) air temperature under A2 and B1 scenarios

    圖9 三種GCMs(a—b:CSIRO30 ,c—d:INM,e—f:MRI)在不同情景(a,c,e:A2情景,b,d,f:B1情景)的降水量模擬值Fig. 9 Simulated monthly precipitation from three different GCMs (a-b: CSIRO30, c-d: INM, e-f: MRI) under two different scenarios (a, c, e: A2 scenario, b, d, f: B1 scenario)

    圖10 黃河上游位置圖Fig. 10 Location of the Upper Yellow River Basin

    圖11 蘭州水文站月徑流實(shí)測(cè)值與模擬值擬合曲線(a為率定期;b為驗(yàn)證期)Fig. 11 Comparison of monthly discharge between observed and simulated values at the Lanzhou hydrological station (a: calibration period, b: validation period)

    圖12 CSIRO模擬的流域各子流域多年平均徑流深在2046—2065年(a,b)和2081—2100年(c,d)較基準(zhǔn)期變化的空間分布Fig. 12 The CSIRO outputs: changes of averaged annual runoff for 2046-2065 (a, b) and 2081-2100 (c, d) compared with those in the baseline period

    2)黃河上游流域徑流量對(duì)未來氣候變化的響應(yīng)

    圖12—14分別給出了三個(gè)氣候模式降尺度情景下各子流域多年平均徑流深在2046—2065年和2081—2100年兩個(gè)時(shí)期內(nèi)較基準(zhǔn)期變化的空間分布。

    CSIRO模式相關(guān)組合情景下流域徑流深的變化幅度在-35.8~43mm,其中流域大多數(shù)地區(qū)的徑流深的變化量在-5~-20mm。整體來看,流域徑流量在不同組合情景下的變化趨勢(shì)較為一致,流域內(nèi)大部分子流域徑流量呈減少趨勢(shì)。還可以看出,兩個(gè)時(shí)期徑流深的變化較為相似;對(duì)不同情景而言,A2情景下流域徑流深呈現(xiàn)減少趨勢(shì)的面積要大于B1情景??偟膩碇v,CSIRO氣候模式降尺度情景下,流域徑流深在大部分子流域呈減少趨勢(shì);B1排放情景下,流域徑流深呈減少趨勢(shì)的子流域面積明顯少于A2情景。

    INM模式相關(guān)組合情景下流域徑流深的變化幅度在-33.9~61mm,大部分子流域徑流深變化幅度為負(fù)值,說明流域大部分子流域徑流深呈減少趨勢(shì)。流域徑流深在不同組合情景下的變化趨勢(shì)較為一致,即都表現(xiàn)為流域南部子流域內(nèi)徑流深呈減少趨勢(shì),減少幅度大于20mm,而北部子流域徑流深則較基準(zhǔn)期有增有減。

    MRI模式相關(guān)組合情景下流域徑流深的變化幅度在-32.9~69mm,其中流域內(nèi)有將超過一半流域面積的平均徑流深都呈減少趨勢(shì)。在2046—2065年,A2情景呈增加趨勢(shì)的子流域面積大于B1情景,且A2情景下減少幅度在16~20mm的流域面積大于B1情景;在2081—2100年,B1情景呈增加趨勢(shì)的子流域面積大于A2情景,且B1情景下減少幅度在16~20mm的流域面積大于A2情景??偟膩碚f,MRI模式降尺度情景下流域超過一半面積的徑流深都呈減少趨勢(shì)。

    圖13 INM模擬的流域各子流域多年平均徑流深在2046—2065年(a,b)和2081—2100年(c,d)較基準(zhǔn)期變化的空間分布Fig. 13 The INM outputs: changes of averaged annual runoff for 2046-2065 (a, b) and 2081-2100 (c, d) compared with those in the baseline period

    圖14 MRI模擬的流域各子流域多年平均徑流深在2046—2065年(a,b)和2081—2100年(c,d)較基準(zhǔn)期變化的空間分布Fig. 14 The MRI outputs: changes of averaged annual runoff for 2046-2065 (a, b) and 2081-2100 (c, d) compared with those in the baseline period

    4 結(jié)論

    水循環(huán)作為氣候系統(tǒng)的重要組成部分,其對(duì)氣候變化的反應(yīng)尤為敏感。探討氣候變化對(duì)流域水循環(huán)和水資源的影響,已經(jīng)成為研究工作的熱點(diǎn)[46]。本研究以長(zhǎng)江和黃河典型流域?yàn)檠芯繉?duì)象,通過對(duì)全球氣候模式的適應(yīng)性評(píng)估以及降尺度處理,全球氣候模式與流域水文模型的耦合,定量探討了氣候變化對(duì)流域水循環(huán)的影響。結(jié)論如下:

    1)黃河流域表現(xiàn)最好的3個(gè)氣候模式為MRI、HadGEM1和INM;長(zhǎng)江流域FGOALS、ECHAM4、ECHAM5的表現(xiàn)最好,進(jìn)一步表明我國(guó)南北區(qū)域氣候變化影響差異顯著。

    2)利用VIC模型,對(duì)長(zhǎng)江下游太湖流域水文過程進(jìn)行模擬,發(fā)現(xiàn)2046—2065年較基準(zhǔn)期徑流深減少幅度最大的地區(qū)是浙西區(qū)與杭嘉湖區(qū);浙西區(qū)在汛期呈現(xiàn)略微增加的趨勢(shì)。2081—2100年較基準(zhǔn)期而言,浙西區(qū)的徑流深雖有減少的趨勢(shì),不過相比2046—2065年減少幅度有所緩解。

    3)SWAT模型在黃河上游流域應(yīng)用得到了相對(duì)滿意結(jié)果。A2和B1情景下,黃河上游流域未來徑流量整體呈減少趨勢(shì),且2081—2100年的減少幅度較2046—2065年大。不同氣候變化組合情景下,黃河上游流域徑流量的空間變化表現(xiàn)出了一定的一致性,流域大部分區(qū)域徑流深呈減少趨勢(shì),且在紅原、久治、若爾蓋地區(qū)減少幅度最大。

    鑒于氣候變化以及流域自身水文循環(huán)的復(fù)雜性和不確定性,未來工作應(yīng)結(jié)合AR5最新研究成果,圍繞以下幾方面開展氣候變化影響研究:1)開展陸地水循環(huán)多要素時(shí)空變化的檢測(cè)與歸因研究,揭示氣候變化下水循環(huán)的響應(yīng)機(jī)理;2)研發(fā)未來水文情景變化的概率預(yù)估理論與方法,定量分析氣候變化對(duì)水資源影響的不確定性以及水資源風(fēng)險(xiǎn)管理與決策問題;3)發(fā)展新一代大尺度分布式陸地水循環(huán)模擬系統(tǒng),實(shí)現(xiàn)陸—?dú)怆p向耦合,揭示氣候變化、人類活動(dòng)與水文循環(huán)相互作用的機(jī)理;4)提出氣候變化影響下水資源脆弱性評(píng)估理論與方法,為應(yīng)對(duì)氣候變化的適應(yīng)性對(duì)策制定奠定基礎(chǔ)。

    [1]IPCC. Climate change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, NY: Cambridge University Press, 2013.

    [2]秦大河. 氣候變化的事實(shí)與影響及對(duì)策. 中國(guó)科學(xué)基金, 2003,17(1): 1-3.

    [3]秦大河, 陳振林, 羅勇, 等. 氣候變化科學(xué)的最新認(rèn)識(shí). 氣候變化研究進(jìn)展, 2007, 3(2): 63-73.

    [4]Hadson D A, Jones R G. Regional climate model simulations of present-day and future climates of southern Africa. Hadley Centre Technical Note 39. London Road, Bracknell, UK, Hadley Centre for Climate Prediction and Research, Met Office, 2002.

    [5]Varis O, Kajander T, Lemmel?, R. Climate and Water: From Climate Models to Water Resources Management and Vice Versa. Climatic Change, 2004, 66(3): 321-344.

    [6]Burn D H, Sharif M, Zhang K. Detection of trends in hydrological extremes for Canadian watersheds. Hydrological Processes, 2010,24(13): 1781-1790.

    [7]Elfert S, Bormann H. Simulated impact of past and possible future land use changes on the hydrological response of the Northern German lowland "Hunte" catchment. Journal of Hydrology, 2010,383(3-4): 245-255.

    [8]Karamouz M, Noori N, Moridi A, et al. Evaluation of floodplain variability considering impacts of climate change. Hydrological Processes, 2011, 25(1): 90-103.

    [9]Jung I W, Chang H. Assessment of future runoff trends under multiple climate change scenarios in the Willamette River Basin,Oregon, USA. Hydrological Processes, 2011, 25(2): 258-277.

    [10]Liu Q, McVicar T R. Assessing climate change induced modification of Penman potential evaporation and runoff sensitivity in a large water-limited basin. Journal of Hydrology,2012, 464: 352-362.

    [11]Marvel K, Bonfils C. Identifying external influences on global precipitation. Proceedings of the National Academy of Sciences of the United States of America (PNAS). 2013, 110(48):19301-19306.

    [12]Piao S, Ciais P, Huang Y, et al. The impacts of climate change on water resources and agriculture in China. Nature, 2010, 467: 43-51.

    [13]第二次氣候變化國(guó)家評(píng)估報(bào)告編寫委員會(huì). 第二次氣候變化國(guó)家評(píng)估報(bào)告. 北京: 科學(xué)出版社, 2011.

    [14]葉篤正, 黃榮輝. 黃河長(zhǎng)江流域旱澇規(guī)律和成因研究. 濟(jì)南: 山東科學(xué)技術(shù)出版社, 1996.

    [15]符淙斌, 安芷生, 郭維棟. 中國(guó)生存環(huán)境演變和北方干旱化趨勢(shì)預(yù)測(cè)研究(Ⅰ): 主要研究成果. 地球科學(xué)進(jìn)展, 2005, 20 (11):1168-1175.

    [16]丁一匯, 任國(guó)玉. 中國(guó)氣候變化科學(xué)概論. 北京: 氣象出版社,2008.

    [17]丁一匯, 任國(guó)玉, 石廣玉, 等. 氣候變化國(guó)家評(píng)估報(bào)告(I): 中國(guó)氣候變化的歷史和未來趨勢(shì). 氣候變化研究進(jìn)展, 2006, 2(1): 3-8.

    [18]陳志愷. 持續(xù)干旱與華北水危機(jī). 中國(guó)水利, 2002, 4: 8-11.

    [19]劉昌明. 黃河流域水循環(huán)演變?nèi)舾蓡栴}的研究. 水科學(xué)進(jìn)展,2004, 15(5): 608-614.

    [20]錢正英, 沈國(guó)舫, 石玉林. 東北地區(qū)有關(guān)水土資源配置、生態(tài)與環(huán)境保護(hù)和可持續(xù)發(fā)展的若干戰(zhàn)略問題研究(綜合卷). 北京: 科學(xué)出版社, 2007.

    [21]邱新法, 劉昌明, 曾燕. 黃河流域近40 年蒸發(fā)皿蒸發(fā)量的氣候變化特征. 自然資源學(xué)報(bào), 2003, 18(4): 437-442.

    [22]楊建平, 丁永建, 陳仁升, 等. 近40a中國(guó)北方降水量與蒸發(fā)量變化. 干旱區(qū)資源與環(huán)境, 2003,17(2): 6-11.

    [23]Liu B H, Xu M, Mark H, et al. A spatial analysis of pan evaporation trends in China, 1955-2000. Journal of Geophysical Research. 2004, 109: 1-9.

    [24]任國(guó)玉, 郭軍. 中國(guó)水面蒸發(fā)量的變化. 自然資源學(xué)報(bào), 2006,21(1): 31-44.

    [25]王國(guó)慶, 王云璋, 康玲玲. 黃河上中游徑流對(duì)氣候變化的敏感性分析. 應(yīng)用氣象學(xué)報(bào), 2002, 1: 117-121.

    [26]劉春蓁. 氣候變化對(duì)江河流量變化趨勢(shì)影響研究進(jìn)展. 地球科學(xué)進(jìn)展, 2007, 22(8): 777-783.

    [27]張建云, 劉九夫. 氣候異常對(duì)水資源影響評(píng)估分析模型. 水科學(xué)進(jìn)展, 2000, 11: 1-9.

    [28]夏軍, Thomas Tanner, 任國(guó)玉, 等. 氣候變化對(duì)中國(guó)水資源影響的適應(yīng)性評(píng)估與管理框架. 氣候變化研究進(jìn)展, 2008, 4(4):215-219.

    [29]林而達(dá), 許吟隆, 蔣金荷, 等. 氣候變化國(guó)家評(píng)估報(bào)告(Ⅱ): 氣候變化的影響與適應(yīng). 氣候變化研究進(jìn)展, 2006, 2(2): 51-56.

    [30]王守榮, 鄭水紅, 程磊. 氣候變化對(duì)西北水循環(huán)和水資源影響的研究. 氣候與環(huán)境研究, 2003, 8(1): 43-51.

    [31]張?jiān)鲂牛?姜彤, 張金池, 等. 長(zhǎng)江流域水汽收支的時(shí)空變化與環(huán)流特征. 湖泊科學(xué), 2008(6): 733-740.

    [32]徐影, 丁一匯, 趙宗慈. 近30年人類活動(dòng)對(duì)東亞地區(qū)氣候變化影響的檢測(cè)與評(píng)估. 應(yīng)用氣象學(xué)報(bào), 2002, 13(5): 513-525.

    [33]Maxino C C, Mcavaney B J, Pitman A J, et al. Ranking the AR4 climate models over the Murray-Darling Basin using simulated maximum temperature, minimum temperature and precipitation. International Journal of Climatology, 2008, 28(8): 1097-1112.

    [34]曹穎, 張光輝. 大氣環(huán)流模式在黃河流域的適用性評(píng)價(jià). 水文,2009, 29(05): 1-5.

    [35]周曉霞, 丁一匯, 王盤興. 夏季亞洲季風(fēng)區(qū)的水汽輸送及其對(duì)中國(guó)降水的影響. 氣象學(xué)報(bào), 2008, 66(1): 59-70.

    [36]初祁, 徐宗學(xué), 蔣昕昊. 兩種統(tǒng)計(jì)降尺度模型在太湖流域的應(yīng)用對(duì)比. 資源科學(xué), 2012, 34(12): 2323-2336.

    [37]袁飛, 謝正輝, 任立良, 等. 氣候變化對(duì)海河流域水文特性的影響. 水利學(xué)報(bào), 2005, 36(3): 274-279.

    [38]劉謙. VIC大尺度陸面水文模型在中國(guó)區(qū)域的應(yīng)用. 長(zhǎng)沙: 湖南大學(xué), 2004.

    [39]Su F G, Xie Z H. A model for assessing effects of climate change on runoff in China. Progress in Natural Science, 2003, 13(9):701-707.

    [40]Lohmann D, Raschke E, Nijssen B, et al. Regional scale hydrology:Ⅱ. Application of the VIC-2L model to the Weser River,Germany. Hydrological Sciences Journal, 1998, 43(1): 143-158.

    [41]Nijssen B, Lettenmaier D P, Liang X, et al. Streamflow simulation for continental-scale river basins. Water Resources Research. 1997,33: 711-724.

    [42]Abdulla F A, Lettenmaier D P, Wood E F, et al. Application of a macroscale hydrologic model to estimate the water balance of the Arkansas-Red River Basin. Journal of Geophysical Research. 1996,101(D3): 7449-7459.

    [43]劉瀏, 徐宗學(xué), 黃俊雄. 氣候變化對(duì)太湖流域徑流的影響. 北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010(3): 371-377.

    [44]Arnold J G, Srinivasan R, Muttiah R S, et al. Large area hydrologic modeling and assessment Part I: Model development. Journal of the American Water Resources Association, 1998, 34(1): 73-89.

    [45]Neitsch S L, Arnold J G, Kiniry J R, et al. Soil and water assessment tool theoretical documentation version 2005. Temple,Texas: Grassland, Soil and Water Research Laboratory, Agricultural Research Service, Blackland Research Center, Texas Agricultural Experiment Station, 2005.

    [46]徐宗學(xué), 劉瀏, 劉兆飛. 氣候變化影響下的流域水循環(huán). 北京: 科學(xué)出版社, 2015.

    Impact of Climate Change on Hydrological Cycle: Case Studies in the Yangtze River and Yellow River Basins

    Xu Zongxue1, Liu Liu2
    (1 College of Water Sciences, Beijing Normal University, Joint Center for Global Change Studies (JCGCS), Beijing 100875 2 College of Water Resources and Civil Engineering, China Agricultural University, Center for Agricultural Water Research in China, Beijing 100083)

    Based upon adaptive assessment of different GCMs recommended by IPCC, the future climate change scenarios were generated by using SDSM and ASD, respectively, and were used to drive the distributed hydrological model VIC and SWAT. The VIC was applied for simulating hydrological processes in the Taihu basin, which is selected as the typical watershed of the Yangtze River basin. The SWAT model was run for simulatimg hydrology in the upper reaches of the Yellow River basin. Then, the impact of climate change on hydrological cycle was quantitatively investigated. Results show that the methods adopted in this study for GCMs adaptive assessment and downscaling could reduce uncertainties effectively. It was detected that a decreasing trend in the upper reaches of the Yellow River basin; a slightly decreasing trend in the lower reaches, of the Yangtze River Basin,but with a significant increasing trend in the southeast of the Taihu basin during flood seasons for the future periods (2046-2065 and 2081-2100), comparing with the runoff in the baseline period (1961-1990). These results are of greatly significance for adapting climate change in different river basins for the future.

    climate change, hydrological cycle, downscaling, hydrological model, General Circulation Model (GCM)

    10.3969/j.issn.2095-1973.2016.03.017

    2015年8月10日;

    2016 年4月1日

    徐宗學(xué)(1962—),Email: zxxu@bnu.edu.cn

    資助信息: 國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2010CB428402);國(guó)家自然科學(xué)基金項(xiàng)目(51509247);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20130008120005)

    猜你喜歡
    太湖流域水循環(huán)徑流
    用袋裝水模擬水循環(huán)
    硫酸法鈦白粉水洗及水循環(huán)膜處理工藝
    化工管理(2021年7期)2021-05-13 00:46:24
    太湖流域低污染水處理技術(shù)研究
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    水循環(huán)高效礦井乏風(fēng)熱泵系統(tǒng)分析與應(yīng)用
    同煤科技(2015年4期)2015-08-21 12:51:02
    探秘“大徑流”
    攻克“大徑流”
    太湖流域:統(tǒng)籌兼顧 多管齊下 積極推進(jìn)節(jié)水型社會(huì)建設(shè)
    太陽(yáng)能水循環(huán)空調(diào)
    江埡水庫(kù)降雨徑流相關(guān)圖的建立

    国产精华一区二区三区| 中文资源天堂在线| 天堂av国产一区二区熟女人妻 | 亚洲中文字幕日韩| 18禁黄网站禁片免费观看直播| 国产成年人精品一区二区| 国产精品精品国产色婷婷| 中国美女看黄片| 久久香蕉激情| 淫妇啪啪啪对白视频| 国产精品国产高清国产av| 999久久久国产精品视频| 中文字幕人成人乱码亚洲影| 国产亚洲精品久久久久5区| 1024手机看黄色片| 日韩欧美在线二视频| 天堂影院成人在线观看| 麻豆国产av国片精品| 长腿黑丝高跟| 黄片小视频在线播放| 90打野战视频偷拍视频| 亚洲av电影在线进入| 少妇裸体淫交视频免费看高清 | 久久久久九九精品影院| 一夜夜www| 国产成人精品久久二区二区91| 亚洲精华国产精华精| 18禁裸乳无遮挡免费网站照片| 国产熟女午夜一区二区三区| xxx96com| 老熟妇乱子伦视频在线观看| 午夜福利18| 久久香蕉国产精品| 成人av在线播放网站| 男女那种视频在线观看| 国产视频一区二区在线看| 五月伊人婷婷丁香| 90打野战视频偷拍视频| 日韩有码中文字幕| 啦啦啦观看免费观看视频高清| 精品不卡国产一区二区三区| av在线播放免费不卡| 俄罗斯特黄特色一大片| e午夜精品久久久久久久| 看片在线看免费视频| 久久香蕉精品热| 国产不卡一卡二| 亚洲人成77777在线视频| 99在线视频只有这里精品首页| 久久伊人香网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av成人不卡在线观看播放网| 亚洲,欧美精品.| 18美女黄网站色大片免费观看| 老司机福利观看| 亚洲国产日韩欧美精品在线观看 | 在线观看舔阴道视频| 国产爱豆传媒在线观看 | netflix在线观看网站| 久久婷婷成人综合色麻豆| 亚洲激情在线av| 国产成人精品无人区| 久久久久久亚洲精品国产蜜桃av| 欧美乱色亚洲激情| 精华霜和精华液先用哪个| 听说在线观看完整版免费高清| 成年人黄色毛片网站| 久久久久精品国产欧美久久久| 少妇熟女aⅴ在线视频| 亚洲精品一区av在线观看| 琪琪午夜伦伦电影理论片6080| 久久九九热精品免费| 久久午夜综合久久蜜桃| 12—13女人毛片做爰片一| 亚洲男人天堂网一区| 嫁个100分男人电影在线观看| 成人亚洲精品av一区二区| 亚洲欧美日韩东京热| 亚洲成人精品中文字幕电影| 国产精品av久久久久免费| 午夜福利视频1000在线观看| 日本三级黄在线观看| 老司机深夜福利视频在线观看| 亚洲av成人不卡在线观看播放网| 亚洲狠狠婷婷综合久久图片| 欧美丝袜亚洲另类 | 国产精华一区二区三区| 精品一区二区三区av网在线观看| 日韩中文字幕欧美一区二区| 久久精品国产亚洲av香蕉五月| 国产高清有码在线观看视频 | 精品一区二区三区视频在线观看免费| 老汉色av国产亚洲站长工具| 后天国语完整版免费观看| 日韩欧美国产在线观看| 亚洲中文字幕日韩| 久久久久久大精品| 久久久久国产一级毛片高清牌| 99国产精品一区二区三区| 亚洲国产精品sss在线观看| 久久精品国产亚洲av高清一级| 欧美色欧美亚洲另类二区| 9191精品国产免费久久| 五月伊人婷婷丁香| 嫁个100分男人电影在线观看| 日本一本二区三区精品| 亚洲全国av大片| 精品久久久久久久毛片微露脸| 成人av在线播放网站| 一级毛片高清免费大全| 午夜影院日韩av| 亚洲欧美一区二区三区黑人| 日韩欧美国产在线观看| 亚洲 国产 在线| 九九热线精品视视频播放| 国内精品久久久久久久电影| 老熟妇乱子伦视频在线观看| 亚洲精品美女久久久久99蜜臀| 久久九九热精品免费| 欧美最黄视频在线播放免费| 村上凉子中文字幕在线| 色精品久久人妻99蜜桃| 一级毛片女人18水好多| 最好的美女福利视频网| 1024视频免费在线观看| 国产精品电影一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人一区二区免费高清观看 | 在线a可以看的网站| 午夜免费观看网址| 热99re8久久精品国产| 国产熟女xx| 巨乳人妻的诱惑在线观看| 99国产精品一区二区蜜桃av| 国产人伦9x9x在线观看| 成人av一区二区三区在线看| or卡值多少钱| 在线播放国产精品三级| 精品电影一区二区在线| 中文字幕熟女人妻在线| 国产精品 欧美亚洲| 久久精品91蜜桃| 亚洲熟妇熟女久久| 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利视频1000在线观看| 黑人操中国人逼视频| 亚洲精品在线观看二区| 久久久久国内视频| 国产亚洲精品综合一区在线观看 | svipshipincom国产片| 久久精品91蜜桃| 久久久久久久午夜电影| 成年女人毛片免费观看观看9| 看黄色毛片网站| 久久人人精品亚洲av| 岛国在线免费视频观看| 中文在线观看免费www的网站 | 国产视频内射| 小说图片视频综合网站| 琪琪午夜伦伦电影理论片6080| 国产区一区二久久| 在线观看66精品国产| 十八禁网站免费在线| 亚洲欧美日韩高清专用| 久久久久久久久中文| 久久国产乱子伦精品免费另类| 69av精品久久久久久| 最近最新中文字幕大全电影3| 老司机靠b影院| 听说在线观看完整版免费高清| 色av中文字幕| 亚洲全国av大片| 99在线人妻在线中文字幕| 这个男人来自地球电影免费观看| 久久久久久九九精品二区国产 | 国产一区二区三区在线臀色熟女| 欧美黑人巨大hd| 国产真实乱freesex| 欧美极品一区二区三区四区| 99国产精品99久久久久| 中文亚洲av片在线观看爽| 搡老岳熟女国产| 亚洲人成77777在线视频| 欧美色视频一区免费| 欧美最黄视频在线播放免费| 久久这里只有精品中国| 国产亚洲欧美在线一区二区| 欧美黑人巨大hd| xxxwww97欧美| 免费一级毛片在线播放高清视频| 亚洲专区国产一区二区| 久久人人精品亚洲av| 亚洲中文日韩欧美视频| 麻豆久久精品国产亚洲av| 国产主播在线观看一区二区| 叶爱在线成人免费视频播放| 1024手机看黄色片| 一级毛片女人18水好多| 亚洲一区中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 最近视频中文字幕2019在线8| 熟女少妇亚洲综合色aaa.| 成人特级黄色片久久久久久久| 国产一区二区激情短视频| 首页视频小说图片口味搜索| 又紧又爽又黄一区二区| 欧美性猛交╳xxx乱大交人| 法律面前人人平等表现在哪些方面| 一级黄色大片毛片| 色综合站精品国产| 正在播放国产对白刺激| 搡老岳熟女国产| tocl精华| 国产精品一区二区三区四区免费观看 | 美女扒开内裤让男人捅视频| 男人舔女人的私密视频| 欧美成人一区二区免费高清观看 | www.自偷自拍.com| 午夜免费激情av| 好看av亚洲va欧美ⅴa在| 美女午夜性视频免费| www.www免费av| 中文字幕av在线有码专区| 亚洲国产看品久久| 免费在线观看亚洲国产| 麻豆久久精品国产亚洲av| 日本a在线网址| 亚洲人成网站高清观看| 最好的美女福利视频网| 一级作爱视频免费观看| 国产精品乱码一区二三区的特点| 欧美黄色片欧美黄色片| 熟女少妇亚洲综合色aaa.| ponron亚洲| 啪啪无遮挡十八禁网站| 欧美乱码精品一区二区三区| 国产精品久久久人人做人人爽| 久99久视频精品免费| 又黄又爽又免费观看的视频| 88av欧美| svipshipincom国产片| 伊人久久大香线蕉亚洲五| 日韩av在线大香蕉| 亚洲成av人片在线播放无| 一进一出好大好爽视频| 成人三级做爰电影| 一进一出抽搐gif免费好疼| 久久精品国产亚洲av高清一级| 国产久久久一区二区三区| 午夜福利高清视频| 香蕉丝袜av| 最新在线观看一区二区三区| 国产亚洲精品一区二区www| 色噜噜av男人的天堂激情| 俺也久久电影网| 日日夜夜操网爽| 人人妻人人澡欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜理论影院| 中文字幕人妻丝袜一区二区| 欧洲精品卡2卡3卡4卡5卡区| 99国产极品粉嫩在线观看| 男插女下体视频免费在线播放| 一本大道久久a久久精品| 天堂动漫精品| 国产av一区在线观看免费| 亚洲真实伦在线观看| 91国产中文字幕| 老司机深夜福利视频在线观看| 丰满人妻一区二区三区视频av | 日本黄色视频三级网站网址| 国产黄a三级三级三级人| 欧洲精品卡2卡3卡4卡5卡区| 五月伊人婷婷丁香| 成年免费大片在线观看| 日韩免费av在线播放| 老鸭窝网址在线观看| 麻豆一二三区av精品| 色综合亚洲欧美另类图片| 女人被狂操c到高潮| 每晚都被弄得嗷嗷叫到高潮| 青草久久国产| 一二三四在线观看免费中文在| 国产精品av久久久久免费| 亚洲av成人精品一区久久| 美女免费视频网站| 欧美色欧美亚洲另类二区| 国产午夜福利久久久久久| 日韩av在线大香蕉| 亚洲午夜理论影院| 哪里可以看免费的av片| svipshipincom国产片| 久久国产乱子伦精品免费另类| 午夜精品久久久久久毛片777| 激情在线观看视频在线高清| 国产精品自产拍在线观看55亚洲| 国产精品电影一区二区三区| 免费在线观看完整版高清| 欧美极品一区二区三区四区| 97人妻精品一区二区三区麻豆| 久久午夜亚洲精品久久| 黑人巨大精品欧美一区二区mp4| 国产黄色小视频在线观看| 日日夜夜操网爽| 国产私拍福利视频在线观看| 可以免费在线观看a视频的电影网站| 一进一出抽搐gif免费好疼| 欧美在线一区亚洲| 特级一级黄色大片| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久久久久久久| 村上凉子中文字幕在线| 中文字幕精品亚洲无线码一区| 在线观看免费视频日本深夜| 日日干狠狠操夜夜爽| 看片在线看免费视频| 精品乱码久久久久久99久播| www.熟女人妻精品国产| 一a级毛片在线观看| 热99re8久久精品国产| 久久精品91蜜桃| 色播亚洲综合网| 女人爽到高潮嗷嗷叫在线视频| 免费高清视频大片| 欧美大码av| 在线视频色国产色| 日本一区二区免费在线视频| 三级国产精品欧美在线观看 | 伦理电影免费视频| 亚洲精品久久成人aⅴ小说| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 一本久久中文字幕| 久热爱精品视频在线9| 欧美日韩乱码在线| 国产亚洲精品一区二区www| www.999成人在线观看| www.精华液| 久久久久亚洲av毛片大全| 午夜精品在线福利| av福利片在线观看| 人成视频在线观看免费观看| 国产欧美日韩一区二区精品| 特级一级黄色大片| 叶爱在线成人免费视频播放| 亚洲午夜理论影院| 变态另类丝袜制服| 黄色a级毛片大全视频| 少妇被粗大的猛进出69影院| 午夜成年电影在线免费观看| 一区福利在线观看| 美女扒开内裤让男人捅视频| 99re在线观看精品视频| 成人国语在线视频| 久久人人精品亚洲av| 精品无人区乱码1区二区| 999久久久国产精品视频| 国产亚洲av嫩草精品影院| 欧美性长视频在线观看| 岛国在线免费视频观看| 欧美性长视频在线观看| 宅男免费午夜| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av| 日韩三级视频一区二区三区| 桃色一区二区三区在线观看| 日本成人三级电影网站| 日韩三级视频一区二区三区| 日日夜夜操网爽| 国产成年人精品一区二区| 黄色丝袜av网址大全| 在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲成人精品中文字幕电影| 亚洲一区高清亚洲精品| 精品久久久久久,| 中文字幕精品亚洲无线码一区| 成人av一区二区三区在线看| bbb黄色大片| 老司机靠b影院| 麻豆国产av国片精品| 欧美一区二区精品小视频在线| 麻豆国产av国片精品| 999精品在线视频| 久久午夜亚洲精品久久| 日韩免费av在线播放| 三级毛片av免费| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 国产在线观看jvid| √禁漫天堂资源中文www| 亚洲片人在线观看| 午夜成年电影在线免费观看| 一进一出好大好爽视频| 人妻丰满熟妇av一区二区三区| 变态另类丝袜制服| 韩国av一区二区三区四区| 在线观看美女被高潮喷水网站 | 久久久久久国产a免费观看| 成人精品一区二区免费| 精品国产乱码久久久久久男人| 国产伦在线观看视频一区| 精品免费久久久久久久清纯| 黄色成人免费大全| 成年人黄色毛片网站| 99久久久亚洲精品蜜臀av| 操出白浆在线播放| 久久久国产欧美日韩av| 成人三级做爰电影| 国产私拍福利视频在线观看| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 特级一级黄色大片| 久久精品91无色码中文字幕| 亚洲国产高清在线一区二区三| tocl精华| 午夜免费成人在线视频| 午夜福利成人在线免费观看| 香蕉av资源在线| 欧美久久黑人一区二区| 色综合站精品国产| 免费在线观看视频国产中文字幕亚洲| 亚洲黑人精品在线| 国产午夜福利久久久久久| 成人av在线播放网站| 精品久久蜜臀av无| 日韩精品青青久久久久久| 亚洲人成伊人成综合网2020| 男人舔奶头视频| 啪啪无遮挡十八禁网站| 一区二区三区激情视频| 国产黄色小视频在线观看| 国产精品一区二区三区四区免费观看 | 欧美激情久久久久久爽电影| 日韩欧美 国产精品| 久久中文字幕一级| 欧美色欧美亚洲另类二区| 波多野结衣巨乳人妻| 一区福利在线观看| 一a级毛片在线观看| 国产精品久久视频播放| 黑人欧美特级aaaaaa片| 国产精品国产高清国产av| 亚洲精品色激情综合| 九九热线精品视视频播放| 一级作爱视频免费观看| 一级毛片女人18水好多| 亚洲国产精品sss在线观看| 国产精品久久久久久人妻精品电影| 成人三级做爰电影| 91字幕亚洲| or卡值多少钱| xxxwww97欧美| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 很黄的视频免费| 欧美性长视频在线观看| 国产免费av片在线观看野外av| 亚洲欧美激情综合另类| 日日爽夜夜爽网站| 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 日韩欧美在线二视频| 亚洲中文字幕一区二区三区有码在线看 | xxx96com| 国产av一区二区精品久久| 亚洲五月婷婷丁香| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 欧美激情久久久久久爽电影| 国产91精品成人一区二区三区| 1024视频免费在线观看| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区| 欧美成狂野欧美在线观看| 亚洲专区国产一区二区| 免费高清视频大片| xxx96com| 亚洲乱码一区二区免费版| 97人妻精品一区二区三区麻豆| 欧美成人午夜精品| 国产一区二区激情短视频| 亚洲精品国产一区二区精华液| 精品国产亚洲在线| 欧美成人性av电影在线观看| 色尼玛亚洲综合影院| 免费无遮挡裸体视频| 90打野战视频偷拍视频| www国产在线视频色| 在线观看www视频免费| 长腿黑丝高跟| 成人av在线播放网站| 曰老女人黄片| 国产精品日韩av在线免费观看| 两个人看的免费小视频| 少妇粗大呻吟视频| 免费在线观看完整版高清| 亚洲欧美激情综合另类| 成人国语在线视频| 亚洲第一电影网av| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美98| 国产av在哪里看| 欧美日韩福利视频一区二区| 午夜激情av网站| 国产熟女午夜一区二区三区| 亚洲中文日韩欧美视频| АⅤ资源中文在线天堂| 亚洲专区字幕在线| 亚洲九九香蕉| 国产精品一区二区免费欧美| 最近视频中文字幕2019在线8| 日本撒尿小便嘘嘘汇集6| 五月伊人婷婷丁香| 成年人黄色毛片网站| 亚洲中文字幕一区二区三区有码在线看 | 国产1区2区3区精品| 丝袜美腿诱惑在线| 国产免费男女视频| 男插女下体视频免费在线播放| 特级一级黄色大片| 久久性视频一级片| 超碰成人久久| 少妇裸体淫交视频免费看高清 | 久久精品国产清高在天天线| 亚洲avbb在线观看| 哪里可以看免费的av片| 一级作爱视频免费观看| 在线观看一区二区三区| 热99re8久久精品国产| 一区福利在线观看| 色尼玛亚洲综合影院| 午夜福利免费观看在线| 久久久久久久久免费视频了| or卡值多少钱| 脱女人内裤的视频| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| av在线天堂中文字幕| 久久久久国产一级毛片高清牌| 99热只有精品国产| √禁漫天堂资源中文www| 亚洲精品粉嫩美女一区| 国产又黄又爽又无遮挡在线| a级毛片a级免费在线| 午夜精品一区二区三区免费看| а√天堂www在线а√下载| 精品国产亚洲在线| 久99久视频精品免费| 亚洲国产看品久久| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美日韩无卡精品| 成人18禁高潮啪啪吃奶动态图| 午夜视频精品福利| 国产69精品久久久久777片 | 久久精品91蜜桃| 男女视频在线观看网站免费 | av国产免费在线观看| svipshipincom国产片| 亚洲国产精品999在线| 正在播放国产对白刺激| av国产免费在线观看| 久久久久久久久免费视频了| 久久中文字幕人妻熟女| 一本大道久久a久久精品| 国模一区二区三区四区视频 | 波多野结衣巨乳人妻| 国产av不卡久久| 亚洲专区字幕在线| 男女之事视频高清在线观看| 亚洲国产精品合色在线| 黄色视频,在线免费观看| 又紧又爽又黄一区二区| 天堂√8在线中文| 国产亚洲精品综合一区在线观看 | 怎么达到女性高潮| 91麻豆精品激情在线观看国产| 午夜日韩欧美国产| 久久精品综合一区二区三区| 在线十欧美十亚洲十日本专区| 丝袜美腿诱惑在线| 亚洲中文字幕一区二区三区有码在线看 | 天天一区二区日本电影三级| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 18禁裸乳无遮挡免费网站照片| av福利片在线观看| 亚洲成人久久爱视频| 午夜福利高清视频| 成人三级做爰电影| 十八禁人妻一区二区| 欧美黄色淫秽网站| 少妇熟女aⅴ在线视频| 香蕉av资源在线| 欧美色视频一区免费| 无限看片的www在线观看| 国产精品香港三级国产av潘金莲| 日本一区二区免费在线视频| 99久久99久久久精品蜜桃| 国产一区在线观看成人免费| 久久久久国内视频| 国内精品久久久久精免费| 亚洲全国av大片| 亚洲中文av在线| 精品午夜福利视频在线观看一区| 色在线成人网| 给我免费播放毛片高清在线观看| 欧美日韩国产亚洲二区| 琪琪午夜伦伦电影理论片6080| 蜜桃久久精品国产亚洲av| 国产精品亚洲一级av第二区| 一个人免费在线观看电影 | 亚洲色图 男人天堂 中文字幕|