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

    近50年氣候和人類活動對錫林河流域徑流的影響

    2015-02-02 02:31:44焦瑋朱仲元宋小園隋佳碩席小康董志兵
    中國水土保持科學(xué) 2015年6期
    關(guān)鍵詞:人類活動徑流量貢獻(xiàn)率

    焦瑋,朱仲元,宋小園,隋佳碩,席小康,董志兵

    (內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,010018,呼和浩特)

    ?

    近50年氣候和人類活動對錫林河流域徑流的影響

    焦瑋,朱仲元?,宋小園,隋佳碩,席小康,董志兵

    (內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,010018,呼和浩特)

    摘要:氣候變化和人類活動是影響河川徑流變化的2大主要因素?;?963—2011年錫林浩特水文站的實(shí)測水文數(shù)據(jù)資料,采用累積距平法(CA)對錫林河近50年徑流年際變化序列的趨勢進(jìn)行分析,確定出2個(gè)拐點(diǎn)年份為1985和1999年,其中1999年為突變年,運(yùn)用修正后的的累積量斜率變化率比較分析法(SCRCQ)定量計(jì)算了氣候及人類活動對徑流變化的貢獻(xiàn)率。以初始研究階段T1(1963—1985年)為基準(zhǔn)期,不考慮潛在蒸散量變化的影響,T2(1986—1999年)和T3(2000—2011年)時(shí)期的降水量減少對于徑流量減少的貢獻(xiàn)率分別為26.30%和57.36%,人類活動對于徑流量減少的貢獻(xiàn)率分別為73.70%和42.64%;如果考慮潛在蒸散量變化,T2、T3時(shí)段的潛在蒸散量減少對于徑流量減少產(chǎn)生的負(fù)貢獻(xiàn)率分別為-23.47%和-52.25%,而人類活動對錫林河流域徑流變化的貢獻(xiàn)率分別增至97.17%和94.89%。顯然,人類活動是該區(qū)間錫林河流域徑流量減少的最重要的驅(qū)動因素。

    關(guān)鍵詞:徑流量; 氣候變化; 人類活動; 貢獻(xiàn)率; 累積距平法; 累積量變化率比較分析法; 錫林河流域

    項(xiàng)目名稱: 國家自然科學(xué)基金“內(nèi)蒙古典型草原水文過程及其擾動與觸發(fā)草地退化的水文臨界條件實(shí)驗(yàn)與模擬研究”(51269014);內(nèi)蒙古自治區(qū)自然科學(xué)基金“天然榆樹疏林草地蒸散發(fā)及其與環(huán)境因子的關(guān)系研究”(2013MS0607);內(nèi)蒙古自治區(qū)科技計(jì)劃項(xiàng)目

    氣候變化和人類活動是影響河川徑流變化的2大主要因素[1]。氣候變化在一定程度上可以導(dǎo)致水資源時(shí)空上的重新分布,甚至?xí)黾雍闈场⒏珊禐?zāi)害發(fā)生的可能性,進(jìn)而影響到區(qū)域水資源的可持續(xù)利用[2-4]。氣候變化對于河川徑流量的影響主要體現(xiàn)在降水量和氣溫的變化。大量研究認(rèn)為降水量的變化是影響產(chǎn)流的主導(dǎo)因素,往往不考慮氣溫因子,然而氣溫變化會導(dǎo)致潛在蒸散量變化,從而影響產(chǎn)流量的大小。在全球變暖的大背景下,我國的區(qū)域氣候逐漸升溫,水資源問題頻現(xiàn),河流徑流量減少[5-7]。近幾十年來,隨著人類社會的進(jìn)步,人口增加、社會經(jīng)濟(jì)快速發(fā)展、城市化進(jìn)程加快,加劇了水資源供需矛盾,水短缺問題日益凸顯;因此,學(xué)術(shù)界廣泛認(rèn)為在研究流域河川徑流量變化時(shí)應(yīng)當(dāng)是氣候和人類活動二者疊加后的效應(yīng)。

    國內(nèi)學(xué)者對于流域徑流變化方面的研究已逐步從定性分析過渡到定量解析的階段[8-12],對于徑流貢獻(xiàn)率分解也進(jìn)行了許多探討[13-17]。我國北方干旱半干旱地區(qū)水資源時(shí)空分布變異性大,人均占有量不均勻,生態(tài)環(huán)境問題頻現(xiàn),導(dǎo)致生態(tài)系統(tǒng)退化嚴(yán)重,在草原型內(nèi)陸河流域地區(qū)尤為顯著[18]。研究此區(qū)域的徑流年際變化情況,定量分解有關(guān)影響因子的貢獻(xiàn)率,為將來流域水資源可持續(xù)利用和有效配置提供參考依據(jù),對于認(rèn)識變化環(huán)境下內(nèi)陸河流域的水文循環(huán)過程演變規(guī)律具有重要意義。

    1研究區(qū)概況

    內(nèi)蒙古自治區(qū)境內(nèi)內(nèi)陸河流域面積約54.5萬km2,其中有69條河流的流域面積超過1 000 km2。錫林河干流全長175 km,流域面積1.1萬多km2,是內(nèi)蒙古典型的草原型內(nèi)陸河。河流發(fā)源于赤峰市克什克騰旗寶爾圖,自東向西流經(jīng)赤峰市和錫林郭勒盟,最后流入查干諾爾沼澤,屬查干淖爾水系。近年來,隨著城市新區(qū)的建立,該流域承載著城市工農(nóng)業(yè)用水及景觀用水的功能,錫林河生態(tài)濕地周邊大力發(fā)展建設(shè),開采利用地下水及地表水,錫林河水庫建成后致使下游斷流,直接影響錫林河流域的徑流變化。錫林河流域位于E 115°32′~117°16′,N 43°26′~44°39′之間,屬中溫帶半干旱大陸性季風(fēng)氣候。春季風(fēng)大干燥,蒸發(fā)強(qiáng)度大;夏季炎熱,雨水多集中在6—8月;秋季涼爽短暫;冬季漫長寒冷,風(fēng)雪交加[19-21]。根據(jù)錫林浩特氣象站1955—2013年資料統(tǒng)計(jì),年均氣溫1.7 ℃,多年平均降雨量282.5 mm,多年平均蒸發(fā)能力為1 862.89 mm。

    2數(shù)據(jù)來源

    圖1 錫林河流域及相關(guān)氣象站水文站位置Fig.1 Xilin River Basin and location of hydro-meteorological stations

    筆者選取錫林浩特水文站1963—2011年的日徑流序列實(shí)測數(shù)據(jù)。錫林浩特水文站設(shè)立于1957年7月,系國家基本站點(diǎn),位于錫林河流域中游,E 116°10′,N 43°49′,控制流域面積3 852 km2。由于該站水文觀測資料時(shí)間序列較長且完整,控制流域出口,因此選擇該站實(shí)測徑流量數(shù)據(jù)作為分析數(shù)據(jù)。研究區(qū)降水?dāng)?shù)據(jù)根據(jù)流域內(nèi)錫林浩特市氣象站,周邊2個(gè)氣象站(阿巴嘎旗和林西縣)1963—2011年的年降水量實(shí)測數(shù)據(jù),利用ArcGIS平臺Kriging方法差值計(jì)算得到。各臺站的年降水?dāng)?shù)據(jù)收集自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http:∥cdc.cma.gov.cn/home.do),各站位置見圖1。

    3研究方法

    累積距平法(CA)[22]是通過觀察差積曲線判斷數(shù)據(jù)點(diǎn)離散程度和變化趨勢的一種非線性統(tǒng)計(jì)方法:

    (1)

    該方法可以判斷離散數(shù)據(jù)與其均值之間離散幅度的大小,本文使用SPSS軟件將序列進(jìn)行量綱歸一化處理后對其進(jìn)行累積距平,若得到的累積距平值>0,表明離散數(shù)據(jù)>均值,反之則<均值,如果差積曲線包括以上2部分,則可用來確定變化趨勢的拐點(diǎn)。CA法可以直接反映不同時(shí)期徑流量的變化進(jìn)程,判斷降水量和徑流量的突變年份,同時(shí)避免了近年來廣泛使用的降水量—徑流量雙累積曲線法(DMC)判斷突變拐點(diǎn)時(shí)的人為因素的局陷性[10],采用累積距平法診斷出的突變點(diǎn)較為清晰,雙累積量在拐點(diǎn)前后的線性擬合程度較好。以確定出的突變年份為界分為基準(zhǔn)期和措施期,與基準(zhǔn)期相比,在措施期人類活動對于徑流的影響擾動較大。

    分析累積徑流量、累積降水量及累積潛在蒸散量趨勢變化時(shí)采用統(tǒng)計(jì)分析方法,定量分解降水量和其他影響因素對河川徑流量的貢獻(xiàn)率應(yīng)用累積量變化率比較分析法(SCRCQ)[10]。該方法可以更為客觀地分析各因子的貢獻(xiàn)率,引入累積量旨在降低實(shí)測年際資料的影響,在影響因子以及徑流量變化不明顯的情形,利用SCRCQ法仍能達(dá)到定量分析影響程度的目的。

    SCRCQ法的原理是,將降水量和其他影響因素的總和定義為100%;各影響因素對因變量的影響程度等于其在時(shí)間序列上的累積斜率變化率與因變量累積斜率變化率的比值。因其具有良好適應(yīng)性的特點(diǎn),本文具體采用修正后的SCRCQ法來研究氣候變化和人類活動對于徑流量的影響程度,計(jì)算過程如下:

    (2)

    式中:CPi為降水量因子對徑流量變化的貢獻(xiàn)率,%;SRi2和SRi1為拐點(diǎn)年份前后累積徑流量-年份線性關(guān)系式的斜率,108m3/a;SPi2和SPi1為累積降水量-年份線性關(guān)系式的斜率,mm/a。

    同理,可求得

    (3)

    式中:CEi為蒸散量(氣溫)因子對徑流量變化的貢獻(xiàn)率,%;SEi2和SEi1為拐點(diǎn)年份前后累積蒸散量-年份線性關(guān)系式的斜率,mm/a。

    在求得上述氣候因子的貢獻(xiàn)率后,根據(jù)水量平衡原理,可得到

    CHi=100-CPi-CEi。

    (4)

    式中CHi為人類活動對徑流量減少的貢獻(xiàn)率,%。

    圖2 錫林河流域1963—2011年徑流量和降水量變化特征Fig.2 Variations in annual runoff and precipitation at the Xilin River Basin from 1963 to 2011

    由于研究區(qū)缺乏長系列持續(xù)觀測的蒸散量數(shù)據(jù),本文采用L. Zhang等[23]提出的張氏公式計(jì)算潛在蒸散量來代替實(shí)測蒸散量。

    (5)

    式中:E為潛在蒸散量,mm;p為年均降水量,mm;q為森林覆蓋率,%。

    采用該方法的前提是在要區(qū)分人類活動和自然因素2方面的影響;因此,排除人工林地只需搜集天然林所占研究區(qū)的面積比例,即天然林覆蓋度[11]。由于研究區(qū)缺乏詳盡的天然林覆蓋度的相關(guān)資料,本文采用中科院遙感所制作的中國區(qū)的GLC2000土地覆蓋數(shù)據(jù),在ArcGIS平臺將GLC2000矢量數(shù)據(jù)轉(zhuǎn)化成grid柵格數(shù)據(jù),根據(jù)研究區(qū)的流域外邊界圖裁剪得到錫林河流域2000年的土地覆蓋數(shù)據(jù),提取常綠針葉林、落葉闊葉林和灌叢的面積,進(jìn)而計(jì)算得到研究區(qū)的森林覆蓋率為3.89%。已有研究[24-26]表明錫林河流域范圍內(nèi)雖然灌木林增幅較大;但因?yàn)槠渥陨碓济娣e很小,所以森林覆蓋率變化不明顯。因此,可近似認(rèn)為在本文研究期內(nèi)流域的森林覆蓋率大致保持該水平不變。

    4結(jié)果與分析

    4.1徑流量和降水量年際變化特征

    錫林河流域在1963—2011年的徑流量年際變化特征如圖2所示。年徑流量在波動中整體以11萬m3/a的速度減少,其中1963—1979年變幅較小且走勢平穩(wěn),1979—1985年有較大幅度下降,1985—2000年間波動幅度較大,1993年達(dá)到最高值,后期總體呈現(xiàn)下降的趨勢,雖有小幅回升但程度有限。

    流域的年均降水量整體以0.963 mm/a的速度減少,19世紀(jì)60年代初至70年代末呈增長態(tài)勢(圖2),到80年代及以后的降水量變化幅度加大,總體具有減少的趨勢。由徑流量和降水量的5點(diǎn)滑動平均值變化可知在80年代中期前后二者均發(fā)生明顯變化而且整體趨勢呈現(xiàn)良好的同步性。

    4.2降水量和徑流量拐點(diǎn)年份

    如圖3所示,降水量-徑流量雙累積曲線(DMC)中雙累積量之間具有非常好的線性擬合關(guān)系(R2>0.98,P<0.000 1),判斷拐點(diǎn)年份為1979年和1999年。可以看到20世紀(jì)80年代前后雙累積曲線的拐點(diǎn)并不明顯,這對拐點(diǎn)年份的選擇產(chǎn)生了不確定因素。

    圖3  錫林河流域年降水量-徑流量雙累積曲線 變化趨勢及拐點(diǎn)年份Fig.3  Changing trend of the double cumulative curve  of precipitation and runoff showing the years  of inflexion points at the Xilin River Basin

    圖4 錫林河流域量綱歸一化序列累積距平曲線Fig.4 Cumulative anomalous curve of standardized annual runoff(a) and precipitation variations(b) at the Xilin River Basin

    如圖4(a)中的徑流量量綱歸一化序列累積距平曲線所示,1963—1979年徑流量的年際變化有小幅波動,1979—1985年間年徑流的累積距平分布小幅度下降,1985年是典型枯水年,1993—2006豐水期,1999年徑流量累積距平達(dá)到峰值,前后總體發(fā)生了增加到減少的突變。圖4(b)顯示在整個(gè)時(shí)間序列錫林河流域的年降水量的拐點(diǎn)年份與同期年徑流量累積距平表現(xiàn)出相似性,1979年由增加轉(zhuǎn)變?yōu)闇p少,以1985年附近為界由減少轉(zhuǎn)為增加,1998年發(fā)生了由增加趨勢到減少趨勢的突變。

    綜上所述,確定徑流量和降水量的拐點(diǎn)年份分別為1985和1999年,其中1999年為突變年。在拐點(diǎn)1985年之前人類活動較為輕微;因此氣候因素的影響是導(dǎo)致徑流量變化的主要因素,將1963—1985年確定為基準(zhǔn)期。1985年之后是氣候因子和人類活動共同作用的措施期。徑流量在1999年發(fā)生突變說明受氣候和人類活動的影響程度加大。

    4.3氣候和人類活動對錫林河流域徑流量減少的貢獻(xiàn)率分析

    以上述確定的拐點(diǎn)年份為界,將各變量累積量的變化過程劃分為3個(gè)時(shí)期1963—1985年(T1),1986—1999年(T2)和2000—2011年(T3),分別對年份與累積徑流量、累積降水量進(jìn)行線性擬合,其相關(guān)線性回歸方程見圖5,相關(guān)程度非常高。根據(jù)SCRCQ法中相關(guān)系數(shù)的獲取原則,將各線性回歸方程的斜率提取至表1。T1為基準(zhǔn)期,之后各時(shí)期的變量變化率都是與基準(zhǔn)期相比,若變化率>0,趨勢減少;<0,則增加。

    3個(gè)時(shí)期的年份-累積徑流量的線性擬合方程的斜率分別為0.178/(108m3·a-1),0.234/(108m3·a-1),0.114/(108m3·a-1)(表1),與基準(zhǔn)期T1相比,T2時(shí)段徑流量減少了31.65%,T3時(shí)段徑流量增長了35.75%;T3與T2相比徑流量增加了51.20%,增幅較為明顯。根據(jù)公式(2)~(4)分別估算各影響因素對于徑流量減少的貢獻(xiàn)率,T2和T3時(shí)段的貢獻(xiàn)率均是與基準(zhǔn)期作比較得到,結(jié)果見表2。在不考慮潛在蒸散量的影響時(shí),措施期T2(1985—1999)和T3(2000—2011)時(shí)段降水量減少對于徑流量減少的貢獻(xiàn)率分別為26.30%和57.36%,人類活動對于徑流減少的貢獻(xiàn)率分別為73.70%和42.64%,可見在T3時(shí)期,人類活動對徑流量影響程度最大。如果考慮潛在蒸散量變化,在措施期潛在蒸散量的減少對徑流變化產(chǎn)生的負(fù)貢獻(xiàn)率分別為-23.47%和-52.25%,人類活動的貢獻(xiàn)率分別增至97.17%和94.89%。顯而易見,人類活動是該流域地表徑流量縮減的最重要驅(qū)動因素。

    SD為標(biāo)準(zhǔn)差。SD is standard deviation圖5 錫林河流域年份與累積徑流量和累積降水量之間的關(guān)系Fig.5 Changes of variations in cumulative runoff and precipitation over years at the Xilin River Basin

    表1 不同時(shí)期年份與各累積量斜率變化率

    Note:SRishows that the slopes of the linear relationship between year and cumulative runoff.SPishows that the slopes of the linear relationship between year and cumulative precipitation.SEishows that the slopes of the linear relationship between year and cumulative potential evapotranspiration.

    表2不同時(shí)期氣候因子和人類活動對徑流量變化的貢獻(xiàn)率

    Tab.2Quantitative impact of climatic factors and human activities on runoff change at the Xilin River Basin during different periods

    時(shí)期PeriodsCPi/%CEi/%(CEi+CHi)/%CHi/%T1————T226.30-23.4773.7097.17T357.36-52.2542.6494.89

    Note:CPiis the contribution rate of precipitation factor of runoff change;CEiis the contribution rate of evaporation factor of runoff change;CHiis the contribution rate of human activities on the reduction of runoff.

    4.4人類活動因素分析

    結(jié)果表明:錫林河流域徑流變化主要受人類活動影響,下面對不同人類活動的影響作進(jìn)一步分析。以流域所在的錫林郭勒盟為例,分別從人口及載畜量變化、生態(tài)保護(hù)措施以及水庫建設(shè)等方面對此進(jìn)行討論。

    錫林河流域自20世紀(jì)80年代以來,隨著人口急劇增長,城市化進(jìn)程加快,城鎮(zhèn)用地面積和有效灌溉面積均有不同程度的增加[27]。據(jù)統(tǒng)計(jì),錫林郭勒盟人口從1960年的50.34萬人增至2011年的103.31萬人,在進(jìn)入20世紀(jì)80年代及21世紀(jì)之后增速較快(圖6(a))。錫林郭勒盟有效灌溉面積從1988年的0.17萬hm2增至2011年的3.94萬hm2,增長約22倍,2000年以后大幅增長(圖6(b))。研究區(qū)城鎮(zhèn)用地面積20年間增加36.35 km2。

    已有研究[28]表明,錫林河流域自1985年以來其實(shí)際載畜量遠(yuǎn)遠(yuǎn)大于正常水平,超出天然草地的承載能力。統(tǒng)計(jì)顯示,流域范圍內(nèi)牲畜總量在近20年間增加約2倍。過度放牧和超載養(yǎng)畜等不合理的人類活動會引起流域植被的退化,降低地表對降雨徑流的緩沖攔蓄作用,加速水土流失[29]。2000年開始,錫林郭勒盟開展了一系列的生態(tài)保護(hù)工作,包括“退耕還林還草”“圍封禁牧”“天然林保護(hù)工程”“京津風(fēng)沙源治理”等水土保持措施。統(tǒng)計(jì)[30]顯示,錫盟水土流失治理面積從2000年的1.49萬hm2增長至2011年的2.69萬hm2。

    根據(jù)朝博等[31]對內(nèi)蒙古草原水資源開發(fā)利用狀況的研究,為改善草原牧區(qū)生態(tài)環(huán)境,國家實(shí)施了大量生態(tài)恢復(fù)工程,其中還林、造林工程均為水資源“大用戶”,消耗大量水資源。據(jù)水行政主管部門統(tǒng)計(jì),2005年末錫盟機(jī)電井供水量2.69億m3,其中51%用于農(nóng)業(yè)灌溉,19%用于工業(yè)生產(chǎn)。退牧還草工程使部分牧民轉(zhuǎn)變?yōu)檗r(nóng)民,開墾草原種植作物,表面減輕了草場的負(fù)擔(dān),實(shí)則加重了對水資源的需求。王軍等[20]96對錫林河流域水資源狀況的研究表明,該地區(qū)煤炭資源豐富,大量工礦企業(yè)以地下水作為供水水源,造成地下水位下降,市區(qū)地下水開采已接近警戒線[32],流域南部已有降落漏斗形成。2010 年流域水資源開發(fā)利用總量占可利用總量的93.19%,基于目前需水現(xiàn)狀,到2020年錫林河流域?qū)⒂? 513.85萬m3的缺水口,缺水程度達(dá)19.47%。

    水庫對調(diào)節(jié)河流徑流量有著獨(dú)特的功能,2012年全盟已建成大、中、小型水庫35座,總庫容5.2億m3,其中錫林河水庫是一座以防洪為主、兼顧工業(yè)供水及灌溉、養(yǎng)殖和旅游等綜合利用的中型水庫。部分?jǐn)r蓄滯留的地表水用以生態(tài)工程和農(nóng)業(yè)灌溉以及景觀觀賞,這些水量部分通過水循環(huán)作用蒸散到大氣中[29]185,從而間接作用其徑流量的變化。

    錫林河流域從20世紀(jì)80年代中期開始頻繁人類活動,人口和載畜量快速增長,與1985年徑流量出現(xiàn)拐點(diǎn)相一致;進(jìn)入21世紀(jì)后,人口擴(kuò)張和工業(yè)發(fā)展,農(nóng)業(yè)灌溉以及水利工程的修建等人類活動共同影響流域徑流,水土保持工作的展開使人類活動對徑流的影響有所減輕。

    數(shù)據(jù)來源于錫林郭勒盟統(tǒng)計(jì)局。data in the graph from Xilingol league area Bureau of Statistics. 圖6 錫林郭勒盟部分人類活動變化過程Fig.6 Variations of some human activities in Xilin Gol League

    5結(jié)論與討論

    1)1963—2011年間,錫林河流域徑流量和降水量分別以11萬m3/a和0.963 mm/a的速度減少,在20世紀(jì)80年代中期前后二者均發(fā)生明顯變化而且整體趨勢呈現(xiàn)良好的同步性。

    2)本研究采用累積距平法(CA)診斷拐點(diǎn)年份,確定徑流量和降水量的突變年份為1999年,拐點(diǎn)年份為1985年和1999年。

    3)以T1(1963—1985)為基準(zhǔn)期(人類活動輕微),在不考慮潛在蒸散量變化的影響時(shí),措施期T2(1985—1999)和T3(2000—2011)時(shí)段降水量減少對于徑流變化的貢獻(xiàn)率分別為26.30%和57.36%,人類活動對于徑流減少的貢獻(xiàn)率分別為73.70%和42.64%;如果考慮潛在蒸散量變化,在T2和T3時(shí)期潛在蒸散量變化對徑流減少產(chǎn)生的負(fù)貢獻(xiàn)率分別為-23.47%和-52.25%,而人類活動的貢獻(xiàn)率分別增至97.17%和94.89%。

    4)人類活動是導(dǎo)致錫林河流域徑流量減少的最重要的驅(qū)動因素,其中人口及載畜量變化,工業(yè)發(fā)展與農(nóng)業(yè)灌溉以及水利工程的修建等因子共同影響流域徑流,水土保持工作的展開使人類活動對徑流的影響有所減輕。

    前人對于錫林河流域徑流量變化的影響因素已做過一些研究,但楊力哲等[26]307沒有細(xì)致地分析徑流在年際上的突變特征,僅判斷出徑流量累積值發(fā)生偏移的年代,未查明拐點(diǎn)年份,只定性的分析了氣候因素、人類活動及下墊面變化等潛在因素對錫林河徑流產(chǎn)生的影響;而本研究通過累積距平法確認(rèn)出的拐點(diǎn)年份與楊力哲的研究結(jié)果一致,并且通過劃分不同時(shí)段對徑流變化的影響因素進(jìn)行了定量解析,結(jié)果更具說服力。本文分解出的人類活動影響的貢獻(xiàn)率包含了所有類型的人類活動對徑流變化的影響,由于各因子間的相互作用十分復(fù)雜;因此,如何進(jìn)一步定量區(qū)分各因子對徑流變化的影響有待深入研究。

    參考文獻(xiàn)6

    [1]Zhan Chesheng, Zeng Sidong, Jiang Shanshan, et al. An integrated approach for partitioning the effect of climate change and human activities on surface runoff [J]. Water Resources Management, 2014, 28(11):3843-3858

    [2]吳紹洪, 黃季焜, 劉燕華, 等. 氣候變化對中國的影響利弊[J].中國人口·資源與環(huán)境,2014,24(1):7-13

    [3]王剛, 肖偉華, 路獻(xiàn)品, 等. 氣候變化對旱澇事件影響研究進(jìn)展[J].災(zāi)害學(xué),2014,29(2):142-148

    [4]夏軍, 劉春蓁, 任國玉. 氣候變化對我國水資源影響研究面臨的機(jī)遇與挑戰(zhàn)[J].地球科學(xué)進(jìn)展,2011,26(1):1-12

    [5]邢子強(qiáng), 嚴(yán)登華, 魯帆, 等. 人類活動對流域旱澇事件影響研究進(jìn)展[J].自然資源學(xué)報(bào),2013,28(6):1070-1082

    [6]夏軍, 邱冰, 潘興瑤, 等. 氣候變化影響下水資源脆弱性評估方法及其應(yīng)用[J].地球科學(xué)進(jìn)展,2012,27(4):443-451

    [7]陳亞寧, 楊青, 羅毅, 等. 西北干旱區(qū)水資源問題研究思考[J]. 干旱區(qū)地理,2012,35(1):1-9

    [8]趙曉坤, 王隨繼, 范小黎. 1954—1993年間窟野河徑流量變化趨勢及其影響因素分析[J]. 水資源與水工程學(xué)報(bào),2010,21(05):32-36

    [9]李志, 劉文兆, 鄭粉莉, 等. 黃土塬區(qū)氣候變化和人類活動對徑流的影響[J]. 生態(tài)學(xué)報(bào),2010,30(9):2379-2386

    [10] 王隨繼, 閆云霞, 顏明, 等. 皇甫川流域降水和人類活動對徑流量變化的貢獻(xiàn)率分析:累積量斜率變化率比較方法的提出及應(yīng)用[J]. 地理學(xué)報(bào),2012,67(3):388-397

    [11] 王隨繼,李玲,顏明. 氣候和人類活動對黃河中游區(qū)間產(chǎn)流量變化的貢獻(xiàn)率[J]. 地理研究,2013,32(3):395-402

    [12] 張調(diào)風(fēng), 朱西德, 王永劍, 等. 氣候變化和人類活動對湟水河流域徑流量影響的定量評估[J]. 資源科學(xué),2014,36(11):2256-2262

    [13] 王寶鑒, 宋連春, 張強(qiáng), 等. 石羊河流域水資源對氣候變暖的響應(yīng)及對生態(tài)環(huán)境的影響[J]. 地球科學(xué)進(jìn)展,2007,22(7):730-737

    [14] 李加強(qiáng), 陳亞寧, 李衛(wèi)紅, 等. 天山北坡中小河川降水與徑流變化特征:以精河為例[J]. 干旱區(qū)地理,2010,33(4):615-622

    [15] 陳忠升, 陳亞寧, 李衛(wèi)紅, 等. 塔里木河干流徑流損耗及其人類活動影響強(qiáng)度變化[J]. 地理學(xué)報(bào),2011,66(1):89-98

    [16] 鄧振鏞, 張強(qiáng), 王潤元, 等. 河西內(nèi)陸河徑流對氣候變化的響應(yīng)及其流域適應(yīng)性水資源管理研究[J]. 冰川凍土,2013,35(5):1267-1275

    [17] 張曉曉, 張鈺, 徐浩杰, 等. 河西走廊三大內(nèi)陸河流域出山徑流變化特征及其影響因素分析[J]. 干旱區(qū)資源與環(huán)境,2014,28(4):66-72

    [18] 李春梅,高素華. 我國北方干旱半干旱地區(qū)水資源演變規(guī)律及其供需狀況評價(jià)[J]. 水土保持學(xué)報(bào),2002,16(2):68-71

    [19] Tong C,Wu J,Yong S,et al. A landscape-scale assessment of steppe degradation in the Xilin River Basin, Inner Mongolia, China [J].Journal of Arid Environments,2004,59(1):133-149

    [20] 王軍, 李和平, 趙淑銀, 等.錫林河流域水資源評價(jià)與開發(fā)利用潛力分析研究[J]. 水資源與水工程學(xué)報(bào),2011,22(4):95-102

    [21] 楊淇越, 吳錦奎,丁永建, 等. 錫林河流域地表水和淺層地下水的穩(wěn)定同位素研究[J]. 冰川凍土,2009,31(5):850-856

    [22] Ran L S, Wang S J, Fan X L. Channel change at Toudaoguai station and its responses to the operation of upstream reservoirs in the upper Yellow River[J]. Journal of Geographical Sciences, 2010, 20(2):231-247

    [23] Zhang L, Dawes W R, Walker G R. Response of mean annual evapotranspiration to vegetation changes at catchment scale[J]. Water Resources Research, 2001, 37(3):701-708

    [24] 李月臣, 劉春霞. 錫林浩特市土地利用/覆蓋變化研究[J]. 水土保持通報(bào),2006,4(2):60-64

    [25] 徐廣才, 康慕誼, 李亞飛. 錫林郭勒盟土地利用變化及驅(qū)動力分析[J]. 資源科學(xué),2011,33(4):690-697

    [26] 楊立哲, 錢虹, 郝璐. 錫林河近50年徑流變化特征及其影響因素分析[J]. 草業(yè)科學(xué),2015,32(3):303-310

    [27] 內(nèi)蒙古自治區(qū)統(tǒng)計(jì)局.內(nèi)蒙古統(tǒng)計(jì)年鑒[M].北京:中國統(tǒng)計(jì)出版社,2012:493-517

    [28] 姜曄, 畢曉麗, 黃建輝, 等. 內(nèi)蒙古錫林河流域植被退化的格局及驅(qū)動力分析[J]. 植物生態(tài)學(xué)報(bào),2010,34(10):1132-1141

    [29] 李玉蘭, 賈國華. 錫林河水庫水土流失現(xiàn)狀及成因分析[J]. 安徽農(nóng)學(xué)通報(bào)(上半月刊),2012,18(21):185-188

    [30] 錫林郭勒盟黨史地方志編纂委員會辦公室. 錫林郭勒年鑒[M].呼和浩特:內(nèi)蒙古出版集團(tuán)內(nèi)蒙古文化出版社,2012:168

    [31] 朝博, 烏云,烏恩. 氣候變化背景下內(nèi)蒙古草原水資源保護(hù)與可持續(xù)利用[J].中國草地學(xué)報(bào),2012.34(5):99-106

    [32] 趙為耕. 淺談錫盟地下水的可持續(xù)利用[J]. 內(nèi)蒙古水利,2014(1):99-100

    (責(zé)任編輯:程云郭雪芳)

    Impacts of climate change and human activities on runoff yield of the

    Xilin River Basin over nearly 50 years

    Jiao Wei, Zhu Zhongyuan, Song Xiaoyuan, Sui Jiashuo, Xi Xiaokang, Dong Zhibing

    (Hydraulics and Civil Engineering College,Inner Mongolia Agricultural University, 010018, Hohhot, China)

    Abstract:Climate change and human activities are two main factors influencing the runoff change. To some extent, climate change can lead to the spatial-temporal redistribution of water resources. The effect of climate change on river runoff yield is mainly reflected in the change of precipitation and temperature. With the increase of population, rapid development of economy and society, and the accelerated process of urbanization, a serious contradiction arise between supply and demand of water resources, and the question of water resource shortage is increasingly concerned. A prevailing thought is that research on river runoff changes should take into account the combination of climate change and human activities. In this study, Xilin River, a typical inland river basin in arid and semi-arid region of Inner Mongolia, was selected as our research object. Based on the observed data of Xilinhot Hydrological Station during 1963—2011, we used the method of cumulative anomaly (CA) to analyze the runoff sequence interannual trend of Xilin River over nearly 50 years, in order to determine the years of inflection point, i.e., 1985 and 1999. Using the revised slope change ratio of cumulative quantity (SCRCQ), we calculated quantitatively the contribution rate of the precipitation, potential evapotranspiration and human activities on runoff. With years of inflection point as a boundary, change process of cumulants of the variables was divided into three periods, 1963—1985 (T1), 1986—1999 (T2) and 2000—2011 (T3). The linear regression equation was adopted to analyze the relationship between accumulated runoff and time. The slope of the linear regression equation was obtained, and then the rate of change was calculated. Compared with the benchmark stage T1, runoff during the T2period decreased by 31.65%, and that during T3period increased by 35.75%. If not considering the effect of evapotranspiration in the area, the reduced precipitation contributed to the runoff reduction by respectively 26.30% and 57.36% respectively during the period of T2and T3, and the contribution rate of human activity to the runoff reduction was 73.70% and 42.64% respectively, with the initial research phase T1(1963—1985) as a benchmark. On the contrary, regarding to the influence of potential evapotranspiration, the contributions of reduced potential evapotranspiration to the runoff yield change were -23.47% and -52.25% during the period of T2and T3respectively, and the contributions of human activities to the runoff yield change were increased to 97.17% and 94.89% respectively. Obviously, human activities were the dominant factor influencing the runoff yield changes of the Xilin River Basin. The human activities include population explosion, livestock grazing, industrial development, agricultural irrigation and water conservancy projects, and the influence of human activities on runoff is mitigated because of soil and water conservation projects. Since the interaction between various factors is very complex, how to further quantitatively measure influence of each factor on runoff changes need intensive study.

    Keywords:runoff yield; climate change; human activities; contribution rate; cumulative anomaly (CA); slope change ratio of cumulative quantity; Xilin River Basin

    通信作者?簡介: 朱仲元(1956—),男,教授,博士。主要研究方向:水文水資源與草原生態(tài)。E-mail: nmgzzy@tom.com

    作者簡介:第一 焦瑋(1992—),女,碩士研究生。主要研究方向:水資源利用與保護(hù)。E-mail: jw731011410@126.com

    收稿日期:2015-05-15修回日期: 2015-11-03

    中圖分類號:K903

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1672-3007(2015)06-0012-08

    猜你喜歡
    人類活動徑流量貢獻(xiàn)率
    一種通用的裝備體系貢獻(xiàn)率評估框架
    關(guān)于裝備體系貢獻(xiàn)率研究的幾點(diǎn)思考
    水文比擬法在計(jì)算河川徑流量時(shí)的修正
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    В первой половине 2016 года вклад потребления в рост китайской экономики достиг 73,4 процента
    中亞信息(2016年10期)2016-02-13 02:32:45
    人為因素主導(dǎo)下海岸帶生態(tài)系統(tǒng)變遷探討
    綠色科技(2015年6期)2015-08-05 18:01:22
    山地表生地球化學(xué)研究現(xiàn)狀與展望
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    資江流域徑流量演變規(guī)律研究
    日韩大码丰满熟妇| 狂野欧美激情性xxxx| 亚洲中文av在线| 亚洲精品国产区一区二| 大片免费播放器 马上看| 狂野欧美激情性xxxx| 高清黄色对白视频在线免费看| 另类精品久久| 好男人视频免费观看在线| 一本—道久久a久久精品蜜桃钙片| 免费女性裸体啪啪无遮挡网站| 亚洲成人av在线免费| h视频一区二区三区| 女的被弄到高潮叫床怎么办| 女性被躁到高潮视频| 精品少妇内射三级| 国产精品 欧美亚洲| 桃花免费在线播放| 久久久精品国产亚洲av高清涩受| 日韩av在线免费看完整版不卡| 人体艺术视频欧美日本| 免费少妇av软件| 亚洲精品一区蜜桃| 黄色视频在线播放观看不卡| 欧美人与性动交α欧美软件| 亚洲精品国产一区二区精华液| 中文欧美无线码| 国产有黄有色有爽视频| 国产男女超爽视频在线观看| 黄片小视频在线播放| 久久久久久久国产电影| av不卡在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲美女视频黄频| 亚洲人成网站在线观看播放| 精品一区二区三区av网在线观看 | 亚洲欧美一区二区三区国产| av在线观看视频网站免费| 亚洲图色成人| 亚洲欧洲精品一区二区精品久久久 | av在线app专区| 大码成人一级视频| 久久久久精品人妻al黑| 777久久人妻少妇嫩草av网站| av福利片在线| 巨乳人妻的诱惑在线观看| 欧美日韩av久久| 男男h啪啪无遮挡| 9色porny在线观看| 久久久精品国产亚洲av高清涩受| 国产一区二区三区综合在线观看| 老熟女久久久| 国产老妇伦熟女老妇高清| 9色porny在线观看| 爱豆传媒免费全集在线观看| 中文天堂在线官网| 人妻一区二区av| 99久久人妻综合| 建设人人有责人人尽责人人享有的| 国产免费视频播放在线视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品aⅴ在线观看| 性高湖久久久久久久久免费观看| 欧美激情 高清一区二区三区| 男女床上黄色一级片免费看| 亚洲美女搞黄在线观看| 国产日韩欧美视频二区| 高清欧美精品videossex| 中国国产av一级| 最新在线观看一区二区三区 | 国产乱人偷精品视频| 亚洲av综合色区一区| 免费日韩欧美在线观看| 一区二区三区四区激情视频| 夫妻午夜视频| 搡老乐熟女国产| 视频区图区小说| 丝袜美腿诱惑在线| 亚洲精品久久成人aⅴ小说| 一级黄片播放器| 亚洲国产欧美网| 亚洲国产中文字幕在线视频| 亚洲熟女精品中文字幕| 国产男女超爽视频在线观看| 两个人免费观看高清视频| 亚洲视频免费观看视频| 在线观看人妻少妇| 999久久久国产精品视频| 中文字幕制服av| 日日啪夜夜爽| 精品人妻在线不人妻| 亚洲欧美色中文字幕在线| 国产深夜福利视频在线观看| 国产精品一二三区在线看| 在线免费观看不下载黄p国产| av卡一久久| 少妇猛男粗大的猛烈进出视频| 啦啦啦中文免费视频观看日本| 久久久久精品国产欧美久久久 | 永久免费av网站大全| 亚洲国产精品999| 亚洲国产精品国产精品| 伊人亚洲综合成人网| 在线亚洲精品国产二区图片欧美| 免费日韩欧美在线观看| netflix在线观看网站| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲| 国产黄色视频一区二区在线观看| 久久韩国三级中文字幕| 别揉我奶头~嗯~啊~动态视频 | 日日摸夜夜添夜夜爱| 亚洲精品自拍成人| 久久免费观看电影| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区黑人| 国产精品.久久久| 不卡视频在线观看欧美| 尾随美女入室| 五月开心婷婷网| 欧美日韩综合久久久久久| 天堂中文最新版在线下载| 伊人亚洲综合成人网| 亚洲欧美精品自产自拍| 精品国产露脸久久av麻豆| a级毛片在线看网站| 在线观看免费日韩欧美大片| 免费观看人在逋| 午夜免费男女啪啪视频观看| 久久精品国产a三级三级三级| 日韩 亚洲 欧美在线| 久久韩国三级中文字幕| av在线播放精品| 搡老乐熟女国产| 成人三级做爰电影| 国产成人91sexporn| 哪个播放器可以免费观看大片| 久久久精品国产亚洲av高清涩受| 亚洲一卡2卡3卡4卡5卡精品中文| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产最新在线播放| 国产在视频线精品| 青青草视频在线视频观看| 亚洲av电影在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲婷婷狠狠爱综合网| 免费日韩欧美在线观看| 国产乱人偷精品视频| av有码第一页| 成年美女黄网站色视频大全免费| 在线天堂中文资源库| 捣出白浆h1v1| 男女之事视频高清在线观看 | av线在线观看网站| www.自偷自拍.com| 涩涩av久久男人的天堂| 久久精品久久精品一区二区三区| 欧美精品av麻豆av| 嫩草影院入口| 97在线人人人人妻| 亚洲精品一二三| 亚洲精品美女久久av网站| 久久天堂一区二区三区四区| 人妻一区二区av| videosex国产| 国产成人一区二区在线| 成人午夜精彩视频在线观看| 午夜影院在线不卡| 国产精品麻豆人妻色哟哟久久| 亚洲图色成人| 久久人人97超碰香蕉20202| 啦啦啦啦在线视频资源| 亚洲欧美激情在线| 侵犯人妻中文字幕一二三四区| 日韩一区二区三区影片| 亚洲欧美一区二区三区国产| 亚洲精品中文字幕在线视频| 国产亚洲精品第一综合不卡| 国产伦理片在线播放av一区| 九草在线视频观看| 99精国产麻豆久久婷婷| 欧美在线一区亚洲| 日本色播在线视频| av福利片在线| 国产一区亚洲一区在线观看| 一边亲一边摸免费视频| 熟女少妇亚洲综合色aaa.| 亚洲国产精品一区三区| 亚洲自偷自拍图片 自拍| 亚洲精品成人av观看孕妇| 成年av动漫网址| 亚洲色图综合在线观看| av国产久精品久网站免费入址| 欧美人与性动交α欧美软件| av女优亚洲男人天堂| 亚洲精品第二区| 日韩 亚洲 欧美在线| 一本大道久久a久久精品| 色播在线永久视频| 亚洲国产欧美网| netflix在线观看网站| 丰满少妇做爰视频| 嫩草影院入口| 欧美日韩国产mv在线观看视频| 精品卡一卡二卡四卡免费| 色播在线永久视频| 在线观看www视频免费| 天天躁夜夜躁狠狠躁躁| 久久国产精品男人的天堂亚洲| 不卡视频在线观看欧美| 老司机影院毛片| 久久人妻熟女aⅴ| 中文欧美无线码| 18禁裸乳无遮挡动漫免费视频| 久久久久久久久免费视频了| 天美传媒精品一区二区| 国精品久久久久久国模美| 777久久人妻少妇嫩草av网站| 丰满饥渴人妻一区二区三| 少妇精品久久久久久久| 精品一区二区三卡| 国产又色又爽无遮挡免| 亚洲成人国产一区在线观看 | 极品少妇高潮喷水抽搐| 如何舔出高潮| 18禁观看日本| 午夜精品国产一区二区电影| 中文字幕制服av| 亚洲av成人不卡在线观看播放网 | 欧美精品人与动牲交sv欧美| 一级,二级,三级黄色视频| 最新的欧美精品一区二区| 香蕉国产在线看| 少妇被粗大猛烈的视频| 亚洲美女视频黄频| 国产成人免费观看mmmm| 可以免费在线观看a视频的电影网站 | 操出白浆在线播放| 夫妻午夜视频| 男男h啪啪无遮挡| av女优亚洲男人天堂| 美女视频免费永久观看网站| 一本一本久久a久久精品综合妖精| 只有这里有精品99| 在线观看国产h片| 丁香六月欧美| 久久久久久免费高清国产稀缺| 亚洲在久久综合| 叶爱在线成人免费视频播放| 一区在线观看完整版| www.熟女人妻精品国产| 一区二区三区乱码不卡18| 婷婷成人精品国产| 亚洲第一av免费看| 国产在线免费精品| 麻豆av在线久日| 99热网站在线观看| 日本爱情动作片www.在线观看| 亚洲一区中文字幕在线| 成年av动漫网址| 国产精品国产av在线观看| 国产高清国产精品国产三级| 国产日韩一区二区三区精品不卡| 王馨瑶露胸无遮挡在线观看| 9热在线视频观看99| 欧美日韩亚洲高清精品| 久久久国产欧美日韩av| 国产在线视频一区二区| 久久久亚洲精品成人影院| 人成视频在线观看免费观看| 只有这里有精品99| 最近最新中文字幕免费大全7| 国产成人av激情在线播放| 国产成人欧美在线观看 | av在线播放精品| 在现免费观看毛片| 悠悠久久av| 久久久精品免费免费高清| 波野结衣二区三区在线| 日韩av不卡免费在线播放| 少妇人妻 视频| 一级片'在线观看视频| 精品卡一卡二卡四卡免费| 女的被弄到高潮叫床怎么办| 熟女av电影| 色精品久久人妻99蜜桃| 成人国语在线视频| 久久青草综合色| 国产亚洲av高清不卡| 日本欧美视频一区| 丝袜美腿诱惑在线| 精品国产一区二区久久| 人人澡人人妻人| 最近2019中文字幕mv第一页| 一级爰片在线观看| 亚洲成人av在线免费| 久久这里只有精品19| 亚洲情色 制服丝袜| a级毛片黄视频| 欧美精品高潮呻吟av久久| 成人免费观看视频高清| 大香蕉久久成人网| 亚洲国产中文字幕在线视频| 天天躁日日躁夜夜躁夜夜| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看| 国产精品久久久久久人妻精品电影 | 亚洲成人国产一区在线观看 | 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到| 国产片内射在线| 在现免费观看毛片| 午夜影院在线不卡| 欧美日韩综合久久久久久| 成年女人毛片免费观看观看9 | 国产乱来视频区| 久久毛片免费看一区二区三区| 久久99一区二区三区| 午夜激情久久久久久久| 99九九在线精品视频| 亚洲精品国产av蜜桃| 一本一本久久a久久精品综合妖精| 亚洲国产欧美一区二区综合| 中国国产av一级| 国产亚洲一区二区精品| 久久婷婷青草| 国产成人精品久久二区二区91 | 高清不卡的av网站| 欧美亚洲日本最大视频资源| 捣出白浆h1v1| 日本午夜av视频| 久久久久久久精品精品| 黄网站色视频无遮挡免费观看| 人妻 亚洲 视频| 免费人妻精品一区二区三区视频| 大码成人一级视频| 欧美乱码精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 国产高清不卡午夜福利| 国产精品成人在线| 熟女少妇亚洲综合色aaa.| 一级毛片我不卡| 国产乱人偷精品视频| 日韩,欧美,国产一区二区三区| a级毛片在线看网站| 日韩大片免费观看网站| 一区二区三区乱码不卡18| 精品人妻一区二区三区麻豆| 午夜福利免费观看在线| bbb黄色大片| 国产成人免费无遮挡视频| 伊人久久大香线蕉亚洲五| 人体艺术视频欧美日本| 777米奇影视久久| 久久97久久精品| 在线观看www视频免费| 久久久欧美国产精品| 亚洲国产欧美一区二区综合| 成人影院久久| 亚洲三区欧美一区| 精品少妇黑人巨大在线播放| 男女免费视频国产| 狂野欧美激情性xxxx| 亚洲四区av| av视频免费观看在线观看| 亚洲精品日韩在线中文字幕| 一级毛片电影观看| 亚洲国产精品成人久久小说| 中文字幕制服av| 丝袜喷水一区| 1024香蕉在线观看| 国产一区有黄有色的免费视频| 嫩草影视91久久| 捣出白浆h1v1| 中文欧美无线码| 视频区图区小说| 亚洲国产中文字幕在线视频| 一区二区三区四区激情视频| 在线观看www视频免费| 两个人免费观看高清视频| 亚洲成色77777| 一区在线观看完整版| 欧美国产精品一级二级三级| 久久99热这里只频精品6学生| 一级a爱视频在线免费观看| 亚洲国产欧美日韩在线播放| 日本欧美视频一区| 亚洲久久久国产精品| 亚洲美女视频黄频| 国产精品久久久久成人av| 亚洲国产精品999| 国产成人91sexporn| 欧美精品亚洲一区二区| 成年动漫av网址| 99香蕉大伊视频| 2021少妇久久久久久久久久久| a级片在线免费高清观看视频| 七月丁香在线播放| 99热全是精品| 午夜免费男女啪啪视频观看| 18在线观看网站| 亚洲色图综合在线观看| 国产一区有黄有色的免费视频| 亚洲av国产av综合av卡| 日本av手机在线免费观看| 飞空精品影院首页| 国产一区二区激情短视频 | 国产99久久九九免费精品| 在线亚洲精品国产二区图片欧美| 国产一区二区激情短视频 | 青春草国产在线视频| 日本黄色日本黄色录像| 成人手机av| 美女福利国产在线| 久久亚洲国产成人精品v| 一区二区三区精品91| 日日摸夜夜添夜夜爱| 国产xxxxx性猛交| 亚洲精品久久久久久婷婷小说| 久久97久久精品| 亚洲成人免费av在线播放| 只有这里有精品99| 国产精品嫩草影院av在线观看| 久久久精品免费免费高清| 亚洲av国产av综合av卡| 久久 成人 亚洲| tube8黄色片| 91aial.com中文字幕在线观看| 亚洲成人手机| 男女边吃奶边做爰视频| 女人精品久久久久毛片| 97精品久久久久久久久久精品| 男女床上黄色一级片免费看| 国产成人91sexporn| av女优亚洲男人天堂| 久久国产精品大桥未久av| 成人漫画全彩无遮挡| 日韩欧美一区视频在线观看| 亚洲国产av新网站| 久热爱精品视频在线9| 狠狠婷婷综合久久久久久88av| 在线观看国产h片| 一个人免费看片子| 色精品久久人妻99蜜桃| 欧美乱码精品一区二区三区| 免费在线观看完整版高清| 国产成人欧美在线观看 | 制服丝袜香蕉在线| 国产精品免费大片| 激情视频va一区二区三区| 久久久久久久久久久免费av| 一本色道久久久久久精品综合| 国产 一区精品| 欧美日韩综合久久久久久| 男女边吃奶边做爰视频| 飞空精品影院首页| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 不卡av一区二区三区| 成人亚洲精品一区在线观看| 成人18禁高潮啪啪吃奶动态图| 久久人人97超碰香蕉20202| 国产成人精品无人区| 人人澡人人妻人| 欧美精品高潮呻吟av久久| 久久性视频一级片| 午夜91福利影院| 1024香蕉在线观看| 极品人妻少妇av视频| 最近2019中文字幕mv第一页| 国产成人精品福利久久| 大码成人一级视频| 麻豆乱淫一区二区| 成人亚洲欧美一区二区av| 免费av中文字幕在线| 黄色 视频免费看| 大话2 男鬼变身卡| 亚洲精品久久午夜乱码| 一个人免费看片子| 视频区图区小说| 在线观看www视频免费| 人成视频在线观看免费观看| 国产成人午夜福利电影在线观看| 视频在线观看一区二区三区| 亚洲欧美一区二区三区久久| 日韩精品免费视频一区二区三区| 亚洲精品视频女| 欧美国产精品va在线观看不卡| 久久女婷五月综合色啪小说| 免费观看性生交大片5| 欧美黄色片欧美黄色片| 午夜久久久在线观看| 亚洲国产毛片av蜜桃av| 日韩一卡2卡3卡4卡2021年| 免费看av在线观看网站| 1024视频免费在线观看| 久久精品国产综合久久久| 熟女av电影| 美女扒开内裤让男人捅视频| 97在线人人人人妻| 中文字幕人妻熟女乱码| 悠悠久久av| 亚洲精品视频女| 国产精品久久久久成人av| 黄色视频不卡| 丰满迷人的少妇在线观看| 免费黄频网站在线观看国产| 成人黄色视频免费在线看| 亚洲欧美精品综合一区二区三区| 波野结衣二区三区在线| av一本久久久久| 国产成人系列免费观看| 青草久久国产| 老司机亚洲免费影院| 精品国产露脸久久av麻豆| 综合色丁香网| 欧美黑人精品巨大| 色婷婷av一区二区三区视频| 久久久久久久大尺度免费视频| 最新的欧美精品一区二区| 国产精品熟女久久久久浪| 久久精品人人爽人人爽视色| 久久青草综合色| 国产成人一区二区在线| √禁漫天堂资源中文www| 免费少妇av软件| 久久久精品国产亚洲av高清涩受| 热re99久久精品国产66热6| av网站免费在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片| 国产亚洲精品第一综合不卡| 午夜福利网站1000一区二区三区| 久久精品国产综合久久久| 国产一卡二卡三卡精品 | 哪个播放器可以免费观看大片| 一二三四在线观看免费中文在| 黄频高清免费视频| 日韩欧美精品免费久久| 九草在线视频观看| 亚洲av电影在线观看一区二区三区| 国产日韩欧美在线精品| 波多野结衣一区麻豆| 国产精品99久久99久久久不卡 | 校园人妻丝袜中文字幕| 在线观看一区二区三区激情| 精品国产乱码久久久久久小说| 国产乱来视频区| 亚洲av日韩在线播放| 亚洲成人手机| av有码第一页| 高清不卡的av网站| www.av在线官网国产| 丁香六月欧美| 19禁男女啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 1024香蕉在线观看| 免费久久久久久久精品成人欧美视频| 极品少妇高潮喷水抽搐| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 青草久久国产| 欧美亚洲日本最大视频资源| av免费观看日本| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 热99国产精品久久久久久7| 国产av码专区亚洲av| 少妇猛男粗大的猛烈进出视频| 人人妻人人爽人人添夜夜欢视频| 精品一区在线观看国产| 欧美在线黄色| 各种免费的搞黄视频| 国产在线一区二区三区精| 中文字幕人妻熟女乱码| 飞空精品影院首页| 两性夫妻黄色片| 妹子高潮喷水视频| 街头女战士在线观看网站| 欧美成人午夜精品| 亚洲美女黄色视频免费看| 超碰97精品在线观看| 两性夫妻黄色片| 飞空精品影院首页| 少妇被粗大猛烈的视频| 精品亚洲乱码少妇综合久久| 日韩不卡一区二区三区视频在线| 中文字幕亚洲精品专区| 国产国语露脸激情在线看| 亚洲av在线观看美女高潮| 少妇人妻久久综合中文| 欧美av亚洲av综合av国产av | 黄网站色视频无遮挡免费观看| 国产黄频视频在线观看| 下体分泌物呈黄色| 电影成人av| 久久久久国产精品人妻一区二区| 十八禁人妻一区二区| 韩国av在线不卡| 亚洲美女黄色视频免费看| 欧美成人精品欧美一级黄| 久久这里只有精品19| 日本91视频免费播放| 爱豆传媒免费全集在线观看| 国产色婷婷99| 一边亲一边摸免费视频| 两个人免费观看高清视频| av一本久久久久| 国产一卡二卡三卡精品 | 国产成人a∨麻豆精品| 欧美激情极品国产一区二区三区| 日本爱情动作片www.在线观看|