• <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ī)律研究
    狠狠婷婷综合久久久久久88av| 国产成人aa在线观看| 亚洲人成77777在线视频| 成人亚洲精品一区在线观看| 久久女婷五月综合色啪小说| 成人亚洲精品一区在线观看| 国产一区二区三区综合在线观看| 免费看不卡的av| 91精品国产国语对白视频| 国产精品女同一区二区软件| 十分钟在线观看高清视频www| av网站在线播放免费| 一级毛片电影观看| 国产97色在线日韩免费| 日韩中文字幕视频在线看片| 免费不卡的大黄色大毛片视频在线观看| av片东京热男人的天堂| 国产精品亚洲av一区麻豆 | 人妻一区二区av| 你懂的网址亚洲精品在线观看| 成人手机av| 久久久精品94久久精品| 亚洲婷婷狠狠爱综合网| 亚洲国产精品国产精品| 黄色一级大片看看| 国产精品秋霞免费鲁丝片| 久久99热这里只频精品6学生| 99热全是精品| 久久av网站| 成人毛片60女人毛片免费| 色婷婷av一区二区三区视频| 久久综合国产亚洲精品| 国产黄色视频一区二区在线观看| 久久毛片免费看一区二区三区| 肉色欧美久久久久久久蜜桃| 亚洲av成人精品一二三区| 精品第一国产精品| 久久精品久久久久久久性| 999久久久国产精品视频| 一区福利在线观看| 黄色配什么色好看| 大香蕉久久成人网| 久久人人爽av亚洲精品天堂| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 亚洲成色77777| 国产一区二区三区av在线| 制服人妻中文乱码| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 国产白丝娇喘喷水9色精品| 老鸭窝网址在线观看| 可以免费在线观看a视频的电影网站 | 少妇精品久久久久久久| 国产又爽黄色视频| 亚洲人成电影观看| 9色porny在线观看| 秋霞在线观看毛片| 尾随美女入室| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜爱| 少妇 在线观看| 亚洲欧美色中文字幕在线| 精品人妻熟女毛片av久久网站| 中文天堂在线官网| 亚洲综合色惰| 久久久国产欧美日韩av| 亚洲国产精品国产精品| 久久久欧美国产精品| 人妻一区二区av| 在线观看免费视频网站a站| 成人午夜精彩视频在线观看| 永久网站在线| 亚洲精品视频女| 日韩精品有码人妻一区| 两个人看的免费小视频| 久久久久精品久久久久真实原创| av国产久精品久网站免费入址| 在现免费观看毛片| 久久久久久免费高清国产稀缺| 欧美 日韩 精品 国产| 乱人伦中国视频| 少妇被粗大的猛进出69影院| 亚洲av在线观看美女高潮| 99re6热这里在线精品视频| 亚洲激情五月婷婷啪啪| 街头女战士在线观看网站| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 免费在线观看视频国产中文字幕亚洲 | 1024视频免费在线观看| 久久亚洲国产成人精品v| 亚洲美女搞黄在线观看| 秋霞伦理黄片| 最近中文字幕2019免费版| 精品国产露脸久久av麻豆| 国产精品久久久久久精品古装| 日韩在线高清观看一区二区三区| 国产精品久久久久久久久免| 国产又色又爽无遮挡免| 国产xxxxx性猛交| 亚洲精品一二三| 久热久热在线精品观看| 一区二区三区精品91| 亚洲精品,欧美精品| 视频在线观看一区二区三区| 国产熟女欧美一区二区| 久久99蜜桃精品久久| 一本—道久久a久久精品蜜桃钙片| 久久影院123| 国产精品久久久久成人av| 国产片特级美女逼逼视频| 妹子高潮喷水视频| 久久精品国产鲁丝片午夜精品| 国产av一区二区精品久久| 高清视频免费观看一区二区| 国产成人一区二区在线| 亚洲精品美女久久久久99蜜臀 | 国精品久久久久久国模美| 一级毛片我不卡| 国语对白做爰xxxⅹ性视频网站| 国产爽快片一区二区三区| 伊人久久大香线蕉亚洲五| 国产伦理片在线播放av一区| 国产 一区精品| 青青草视频在线视频观看| 日韩中文字幕欧美一区二区 | 尾随美女入室| 色吧在线观看| 性色av一级| 亚洲精品美女久久久久99蜜臀 | 成人黄色视频免费在线看| 成人漫画全彩无遮挡| 国产亚洲一区二区精品| 看十八女毛片水多多多| 大码成人一级视频| 亚洲精品日本国产第一区| 宅男免费午夜| 国产成人91sexporn| 春色校园在线视频观看| 精品国产一区二区三区四区第35| 最新的欧美精品一区二区| 亚洲av电影在线观看一区二区三区| 激情视频va一区二区三区| 亚洲国产日韩一区二区| 精品国产乱码久久久久久小说| 亚洲,一卡二卡三卡| 日日啪夜夜爽| 只有这里有精品99| av天堂久久9| 在现免费观看毛片| 黄色一级大片看看| 一级毛片电影观看| 男女午夜视频在线观看| 黄片无遮挡物在线观看| 中文字幕av电影在线播放| 国产精品蜜桃在线观看| 国产有黄有色有爽视频| 天美传媒精品一区二区| 欧美av亚洲av综合av国产av | 国产精品三级大全| 乱人伦中国视频| 交换朋友夫妻互换小说| 丝袜美足系列| 一区福利在线观看| 搡女人真爽免费视频火全软件| 老司机影院成人| 欧美变态另类bdsm刘玥| 色婷婷av一区二区三区视频| 最近的中文字幕免费完整| 国产毛片在线视频| 爱豆传媒免费全集在线观看| 热re99久久精品国产66热6| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 女人被躁到高潮嗷嗷叫费观| 亚洲激情五月婷婷啪啪| 亚洲精品国产av成人精品| 国产成人午夜福利电影在线观看| 久久久久久伊人网av| 国产精品久久久久久精品电影小说| 黄片小视频在线播放| 精品久久久久久电影网| 亚洲av综合色区一区| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 国产欧美日韩综合在线一区二区| av女优亚洲男人天堂| 高清视频免费观看一区二区| 国产激情久久老熟女| 亚洲精品在线美女| 大话2 男鬼变身卡| 精品少妇久久久久久888优播| 日日啪夜夜爽| 欧美人与性动交α欧美精品济南到 | 国产爽快片一区二区三区| 最近最新中文字幕大全免费视频 | 2018国产大陆天天弄谢| 少妇 在线观看| 午夜福利在线免费观看网站| 亚洲精品国产色婷婷电影| 蜜桃国产av成人99| 亚洲欧美精品自产自拍| 99久久综合免费| 少妇被粗大的猛进出69影院| 日韩伦理黄色片| av不卡在线播放| 99国产精品免费福利视频| 亚洲av.av天堂| 天堂俺去俺来也www色官网| 日本黄色日本黄色录像| 2018国产大陆天天弄谢| 精品亚洲成国产av| 91aial.com中文字幕在线观看| 99九九在线精品视频| 亚洲人成网站在线观看播放| 搡女人真爽免费视频火全软件| 亚洲欧美精品综合一区二区三区 | 久久午夜综合久久蜜桃| 亚洲欧美色中文字幕在线| 亚洲 欧美一区二区三区| 丰满饥渴人妻一区二区三| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 岛国毛片在线播放| 欧美精品人与动牲交sv欧美| 成人国产av品久久久| 91国产中文字幕| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 欧美+日韩+精品| 日本av免费视频播放| 丰满迷人的少妇在线观看| 99热全是精品| 亚洲精华国产精华液的使用体验| 欧美精品人与动牲交sv欧美| 有码 亚洲区| 亚洲欧美清纯卡通| 国产亚洲精品第一综合不卡| 亚洲成色77777| 搡女人真爽免费视频火全软件| 热re99久久国产66热| 欧美成人午夜免费资源| 性高湖久久久久久久久免费观看| 婷婷色麻豆天堂久久| 男女高潮啪啪啪动态图| 好男人视频免费观看在线| av线在线观看网站| 免费观看在线日韩| 精品少妇内射三级| 亚洲三级黄色毛片| 爱豆传媒免费全集在线观看| 十八禁高潮呻吟视频| 亚洲精品久久久久久婷婷小说| 春色校园在线视频观看| 99精国产麻豆久久婷婷| 日本爱情动作片www.在线观看| 久久这里有精品视频免费| 亚洲av.av天堂| 成人二区视频| 只有这里有精品99| 精品少妇一区二区三区视频日本电影 | 午夜福利乱码中文字幕| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久久久免| 街头女战士在线观看网站| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区 | 香蕉精品网在线| 交换朋友夫妻互换小说| 女性被躁到高潮视频| 看免费成人av毛片| 男女国产视频网站| 国产精品香港三级国产av潘金莲 | 街头女战士在线观看网站| a级毛片黄视频| 最近中文字幕2019免费版| 久久久久久久精品精品| 欧美激情 高清一区二区三区| 777米奇影视久久| av线在线观看网站| 少妇人妻 视频| 亚洲精品成人av观看孕妇| 熟女电影av网| 国产免费一区二区三区四区乱码| 黄色一级大片看看| 国精品久久久久久国模美| 1024香蕉在线观看| 国产精品欧美亚洲77777| 亚洲av在线观看美女高潮| 中文欧美无线码| 亚洲国产精品成人久久小说| 亚洲第一区二区三区不卡| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 美女xxoo啪啪120秒动态图| 精品一区二区三卡| 少妇的丰满在线观看| 丝瓜视频免费看黄片| 亚洲美女黄色视频免费看| 黄色视频在线播放观看不卡| 成人毛片60女人毛片免费| 精品久久久精品久久久| 人体艺术视频欧美日本| 一二三四中文在线观看免费高清| 国产探花极品一区二区| 在线观看三级黄色| 欧美在线黄色| 日本vs欧美在线观看视频| 999精品在线视频| 2022亚洲国产成人精品| 少妇 在线观看| 午夜福利网站1000一区二区三区| 欧美 亚洲 国产 日韩一| 黑丝袜美女国产一区| 啦啦啦中文免费视频观看日本| 亚洲第一青青草原| 久久ye,这里只有精品| 在线观看免费视频网站a站| 国产精品人妻久久久影院| 十八禁高潮呻吟视频| 久久精品久久久久久噜噜老黄| 国产精品久久久久久精品古装| 亚洲精品久久久久久婷婷小说| 午夜激情av网站| 男女免费视频国产| 热99国产精品久久久久久7| 叶爱在线成人免费视频播放| 精品国产国语对白av| 亚洲中文av在线| 黄片播放在线免费| 国产伦理片在线播放av一区| 国产综合精华液| 国产av一区二区精品久久| 熟女少妇亚洲综合色aaa.| 日本欧美视频一区| 2018国产大陆天天弄谢| 精品人妻在线不人妻| 国产亚洲欧美精品永久| 国产xxxxx性猛交| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 99热国产这里只有精品6| 女人被躁到高潮嗷嗷叫费观| 最新中文字幕久久久久| 亚洲国产精品国产精品| 亚洲av男天堂| 国产亚洲av片在线观看秒播厂| 极品人妻少妇av视频| 韩国精品一区二区三区| 七月丁香在线播放| 国产亚洲av片在线观看秒播厂| 久久久久久伊人网av| 亚洲伊人色综图| 久久精品国产自在天天线| 国产精品久久久久久精品古装| 91aial.com中文字幕在线观看| 天堂俺去俺来也www色官网| 超色免费av| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 亚洲欧美一区二区三区国产| 日韩av不卡免费在线播放| a级毛片黄视频| 永久免费av网站大全| 99久国产av精品国产电影| 午夜福利在线观看免费完整高清在| 亚洲国产欧美网| 久久久久国产网址| 免费观看性生交大片5| 亚洲精华国产精华液的使用体验| 人人澡人人妻人| 18+在线观看网站| 搡老乐熟女国产| 另类亚洲欧美激情| 成人影院久久| 久久这里有精品视频免费| 美女国产视频在线观看| 老司机亚洲免费影院| 9热在线视频观看99| 在线观看免费视频网站a站| 丰满饥渴人妻一区二区三| 欧美黄色片欧美黄色片| 国产成人免费观看mmmm| 免费黄网站久久成人精品| 只有这里有精品99| 欧美日韩精品网址| 国产男女超爽视频在线观看| 久久久精品94久久精品| 久久狼人影院| 色婷婷av一区二区三区视频| 秋霞在线观看毛片| 男女下面插进去视频免费观看| av国产久精品久网站免费入址| 80岁老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 久久精品久久久久久久性| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区国产| 纵有疾风起免费观看全集完整版| 午夜福利视频精品| 国产亚洲精品第一综合不卡| 最新中文字幕久久久久| 成人二区视频| 国产黄频视频在线观看| 中文字幕精品免费在线观看视频| 热99久久久久精品小说推荐| 色哟哟·www| 国产xxxxx性猛交| 91aial.com中文字幕在线观看| 在线 av 中文字幕| 女人久久www免费人成看片| 亚洲综合色网址| 欧美日韩亚洲高清精品| 边亲边吃奶的免费视频| 秋霞在线观看毛片| 免费观看a级毛片全部| 建设人人有责人人尽责人人享有的| 最新的欧美精品一区二区| 成人二区视频| 极品人妻少妇av视频| 国产精品麻豆人妻色哟哟久久| 爱豆传媒免费全集在线观看| 少妇的逼水好多| 久久久久人妻精品一区果冻| 人人妻人人爽人人添夜夜欢视频| 高清不卡的av网站| 在线天堂最新版资源| 午夜老司机福利剧场| 精品国产乱码久久久久久男人| 热99国产精品久久久久久7| 久久精品人人爽人人爽视色| 好男人视频免费观看在线| 国产有黄有色有爽视频| 人妻少妇偷人精品九色| 一级,二级,三级黄色视频| 人妻系列 视频| 日本午夜av视频| 在线亚洲精品国产二区图片欧美| 亚洲av.av天堂| 美女午夜性视频免费| 黑人猛操日本美女一级片| 亚洲第一青青草原| 91精品国产国语对白视频| 最近中文字幕2019免费版| 免费日韩欧美在线观看| 日韩一本色道免费dvd| 丰满饥渴人妻一区二区三| 欧美日韩一级在线毛片| 丝袜喷水一区| 国产亚洲午夜精品一区二区久久| 久久久欧美国产精品| 美女xxoo啪啪120秒动态图| 九草在线视频观看| 精品国产国语对白av| 亚洲欧美精品自产自拍| 韩国av在线不卡| 69精品国产乱码久久久| 国产野战对白在线观看| 欧美国产精品一级二级三级| 久久久久精品人妻al黑| 嫩草影院入口| 老司机影院毛片| 色婷婷久久久亚洲欧美| 欧美亚洲 丝袜 人妻 在线| 晚上一个人看的免费电影| 亚洲美女视频黄频| 精品一区二区三卡| 久久久久久久久免费视频了| 90打野战视频偷拍视频| 欧美日韩国产mv在线观看视频| 国产av码专区亚洲av| 亚洲精品日本国产第一区| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av蜜桃| 搡女人真爽免费视频火全软件| 九九爱精品视频在线观看| 精品久久蜜臀av无| 久久久欧美国产精品| 一区二区三区四区激情视频| 成人亚洲欧美一区二区av| 曰老女人黄片| 亚洲天堂av无毛| 十八禁高潮呻吟视频| 国产男女内射视频| 不卡视频在线观看欧美| 精品国产国语对白av| 亚洲av欧美aⅴ国产| 婷婷色av中文字幕| 十分钟在线观看高清视频www| 久久精品aⅴ一区二区三区四区 | 天堂俺去俺来也www色官网| 亚洲国产av新网站| 考比视频在线观看| 美女福利国产在线| 精品人妻偷拍中文字幕| 天天躁夜夜躁狠狠躁躁| 国产国语露脸激情在线看| 国产精品久久久av美女十八| 黄片播放在线免费| 中文字幕人妻丝袜制服| 丰满饥渴人妻一区二区三| 国产成人91sexporn| av视频免费观看在线观看| 久久 成人 亚洲| 国产成人精品无人区| 国产精品一二三区在线看| 日韩一卡2卡3卡4卡2021年| 久久精品熟女亚洲av麻豆精品| 多毛熟女@视频| av有码第一页| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| 97在线人人人人妻| 精品一区在线观看国产| 国产xxxxx性猛交| 久久久久久久精品精品| 你懂的网址亚洲精品在线观看| 麻豆精品久久久久久蜜桃| 亚洲综合色惰| 尾随美女入室| 中文字幕精品免费在线观看视频| 青草久久国产| 午夜影院在线不卡| 看免费av毛片| 老司机影院毛片| 亚洲精品第二区| 亚洲色图 男人天堂 中文字幕| 丰满饥渴人妻一区二区三| 只有这里有精品99| 久久久久久久大尺度免费视频| 久久99蜜桃精品久久| 国产精品人妻久久久影院| 久久人人97超碰香蕉20202| 亚洲综合色惰| 丝袜喷水一区| 99久久中文字幕三级久久日本| 免费黄网站久久成人精品| 国产欧美日韩综合在线一区二区| 少妇人妻精品综合一区二区| 久久精品熟女亚洲av麻豆精品| 久久精品久久久久久久性| 亚洲熟女精品中文字幕| 香蕉精品网在线| 久久精品亚洲av国产电影网| 亚洲精品第二区| 大香蕉久久成人网| 亚洲美女搞黄在线观看| 亚洲精品国产色婷婷电影| 99精国产麻豆久久婷婷| 九色亚洲精品在线播放| 日本vs欧美在线观看视频| 久久这里有精品视频免费| av卡一久久| 高清黄色对白视频在线免费看| 国产精品一二三区在线看| 成年女人毛片免费观看观看9 | 美国免费a级毛片| 1024视频免费在线观看| videos熟女内射| 综合色丁香网| 交换朋友夫妻互换小说| 亚洲图色成人| 啦啦啦在线观看免费高清www| 亚洲第一青青草原| 国产精品熟女久久久久浪| 欧美bdsm另类| 秋霞伦理黄片| 亚洲国产最新在线播放| √禁漫天堂资源中文www| 黄色配什么色好看| 97精品久久久久久久久久精品| 亚洲av男天堂| 精品人妻熟女毛片av久久网站| 国产精品熟女久久久久浪| 国产欧美日韩综合在线一区二区| 久久毛片免费看一区二区三区| 日韩免费高清中文字幕av| 青草久久国产| 久久精品熟女亚洲av麻豆精品| 嫩草影院入口| 欧美精品一区二区免费开放| 国产成人精品一,二区| 91在线精品国自产拍蜜月| 男女下面插进去视频免费观看| www.自偷自拍.com| 久久热在线av| 日本-黄色视频高清免费观看| 日韩一区二区三区影片| 99久久中文字幕三级久久日本| 国产精品二区激情视频| 五月天丁香电影| 亚洲国产av影院在线观看| kizo精华| 五月天丁香电影| 日韩制服骚丝袜av| 91精品三级在线观看| a级毛片黄视频| 国产成人91sexporn| 老汉色av国产亚洲站长工具| 国产精品免费视频内射| 18禁裸乳无遮挡动漫免费视频| 岛国毛片在线播放| a 毛片基地| 日韩中字成人| 亚洲人成网站在线观看播放| 超碰成人久久| 搡女人真爽免费视频火全软件| 人妻一区二区av| 777久久人妻少妇嫩草av网站| av免费在线看不卡| 十八禁网站网址无遮挡| 免费黄色在线免费观看| 十分钟在线观看高清视频www|