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

    氣候因素對(duì)阿克蘇河徑流變化影響的定量評(píng)估

    2021-01-27 01:44:38李虹彬劉亞婷王衛(wèi)光龐麥田
    灌溉排水學(xué)報(bào) 2021年1期
    關(guān)鍵詞:徑流量水文徑流

    李虹彬,劉亞婷,王衛(wèi)光,2*,龐麥田

    氣候因素對(duì)阿克蘇河徑流變化影響的定量評(píng)估

    李虹彬1,劉亞婷1,王衛(wèi)光1,2*,龐麥田1

    (1.河海大學(xué) 水文水資源學(xué)院,南京 210098;2.河海大學(xué) 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,南京 210098)

    【目的】探究阿克蘇河徑流變化的歸因。【方法】采用Mann-Kendall-Sneyers檢驗(yàn)估計(jì)流域水文、氣象數(shù)據(jù)的變化趨勢和突變情況,并試圖用3種改進(jìn)的彈性系數(shù)法評(píng)估氣候要素(尤其是溫度)對(duì)阿克蘇河徑流變化的影響?!窘Y(jié)果】徑流和溫度分別以11.73 m3/(s·10 a)和0.26 ℃/10 a的速率顯著遞增。降雨也呈顯著上升趨勢,增幅為14.31 mm/10 a的速率遞增,而潛在蒸散發(fā)呈顯著下降趨勢(21.89 mm/10 a),徑流的突變點(diǎn)為1993年。因此,采用1960—1993年的多年平均徑流和氣候因素作為衡量變化的基準(zhǔn)值。1994—2010年,徑流量相對(duì)1960—1993年增加了23.1%。敏感性分析表明,溫度升高對(duì)徑流變化的相對(duì)貢獻(xiàn)率為46.7%,是徑流量增加的最主要因素,而降水量和潛在蒸散發(fā)的作用分別占徑流量變化的20.5%和24.9%;其他非氣候要素(主要是人類活動(dòng))對(duì)徑流變化的相對(duì)貢獻(xiàn)率約為7.9%?!窘Y(jié)論】氣候變化是阿克蘇河徑流量變化的主要原因,其中溫度的貢獻(xiàn)率最大。

    徑流變化;Mann-Kendall-Sneyers檢驗(yàn);彈性系數(shù);氣候變化

    0 引 言

    【研究意義】近年來,全球氣候的變暖加速了水文循環(huán),加劇了極端水文事件的發(fā)生[1]。隨著全球平均氣溫的升高而導(dǎo)致的水資源短缺問題加劇,氣候變化對(duì)水文循環(huán)的影響已經(jīng)成為近年來水文學(xué)家們關(guān)注的熱點(diǎn)問題[2-3]。由于干旱區(qū)脆弱的水文生態(tài)環(huán)境,氣候變化帶來的水資源短缺等一系列問題尤為突出。此外,我國西北干旱區(qū)的平均氣溫以0.33~0.39 ℃/10 a的速率上升,高于同一時(shí)期的我國平均0.249 ℃/10 a的上升速率,遠(yuǎn)大于全球平均0.13 ℃/10 a的升溫速率[4-5]。這一現(xiàn)象進(jìn)一步加劇了我國西北干旱區(qū)以冰川融雪補(bǔ)給為主的水資源系統(tǒng)的脆弱性和不確定性,同時(shí)也使得該地區(qū)綠洲經(jīng)濟(jì)與沙漠生態(tài)之間的水資源供需矛盾變得更加尖銳[6]。又由于我國干旱區(qū)的徑流補(bǔ)給來源主要是降水與冰雪融水,降雨和溫度的耦合作用將顯著影響徑流的穩(wěn)定性,而作為反映大氣蒸發(fā)能力指標(biāo)的潛在蒸散發(fā)也是影響氣候變化的主要?dú)夂蜃兞?,因此,在干旱區(qū)定量評(píng)估降雨、溫度,潛在蒸散發(fā)等氣候變量以及其他影響因素(主要是人類活動(dòng))對(duì)徑流變化的影響對(duì)旱區(qū)水資源可持續(xù)發(fā)展的管理和布局具有重要意義。

    【研究進(jìn)展】目前,定量探討氣候變化對(duì)徑流影響的研究中,多數(shù)都是基于以下3種方法進(jìn)行:經(jīng)驗(yàn)統(tǒng)計(jì)法,基于敏感性的彈性系數(shù)法和水文模型[7]。江善虎等[8]利用雙累積曲線法確定了氣候變化和人類活動(dòng)對(duì)老哈河流域徑流的影響。黨素珍等[9]通過Budyko框架下的彈性系數(shù)法定量研究了氣候變化和人類活動(dòng)對(duì)涇河上游徑流變化的影響。張艷霞等[10]通過雙累積曲線法、累積量斜率變化率比較法以及彈性分析法定量研究了氣候變化和人類活動(dòng)對(duì)錫林河徑流量影響的貢獻(xiàn)率。王國慶等[11]采用水文模型模擬的途徑定量探討了山西省嵐河流域河川徑流變化受氣候變化和人類活動(dòng)的相對(duì)影響。可以發(fā)現(xiàn),現(xiàn)有的對(duì)徑流變化影響要素的研究多集中在濕潤半濕潤地區(qū),而干旱地區(qū)的此類研究還較少。此外,經(jīng)驗(yàn)統(tǒng)計(jì)法雖然易于使用,但需要非常詳盡的水文氣象資料,耗費(fèi)人力物力,并且對(duì)不同要素之間的book=114,ebook=2相互作用難以有效的考慮[11]。水文模型物理基礎(chǔ)完備,但是模型結(jié)構(gòu)和計(jì)算參數(shù)的不確定性都會(huì)對(duì)結(jié)果造成明顯的誤差影響[12]。而基于模型的氣候彈性法雖然在物理上是合理的,但是對(duì)模型參數(shù)率定的要求較高,并且結(jié)果的準(zhǔn)確性也不很理想[13]。而干旱區(qū)水文、氣象站點(diǎn)稀少且分布不均勻,資料不詳盡,代表性較差,此外,由于干旱區(qū)高溫少雨的氣候條件和脆弱的水文生態(tài)環(huán)境,模型的參數(shù)難以得到有效的率定[14]。因此,現(xiàn)有方法在干旱區(qū)的使用仍具有一定的局限性。而基于非參數(shù)的彈性系數(shù)法,可以直接使用長時(shí)期的直接觀測資料進(jìn)行氣候變化和人類活動(dòng)對(duì)徑流變化影響的定量評(píng)估,其穩(wěn)健性和結(jié)果的合理性都具有一定的保障[15]?!厩腥朦c(diǎn)】然而上述方法都難以直接考慮溫度對(duì)徑流變化的影響,以往考慮溫度變化對(duì)徑流的影響都是通過潛在蒸散發(fā)對(duì)徑流的影響而間接體現(xiàn)。而隨著全球變暖的加劇,溫度對(duì)河川徑流變化的影響將變得更加劇烈,尤其是在寒旱區(qū),其地表徑流主要來自上游冰川融雪來水補(bǔ)給,溫度對(duì)該地區(qū)徑流變化的作用尤為強(qiáng)烈[16],定量評(píng)估以溫度為主的氣候要素對(duì)徑流變化的影響,將有助于加深對(duì)寒旱區(qū)水文過程的理解以及對(duì)未來水資源量的預(yù)測。Li等[17]提出了一種改進(jìn)的彈性系數(shù)法,為考慮溫度對(duì)徑流變化的影響提供了一種行之有效的思路。

    位于我國西北干旱區(qū)的阿克蘇河是塔里木河的主要支流,也是目前唯一與塔里木河干流常年保持自然水力聯(lián)系的源流,其補(bǔ)給水量約占塔里木河干流來水量的73%[18]。同時(shí),阿克蘇河起源于天山山脈,是典型的冰雪融水補(bǔ)給型河流。阿克蘇河徑流量的變化將顯著影響塔里木河干流的來水量,進(jìn)而顯著影響干流的綠洲農(nóng)業(yè)經(jīng)濟(jì)發(fā)展和生態(tài)保護(hù)[19]。近年來,圍繞阿克蘇河流域徑流變化特征以及其影響因素已進(jìn)行了一系列的深入研究,阿克蘇河徑流量逐年增加的現(xiàn)象被廣泛地發(fā)現(xiàn),特別是在20世紀(jì)90年代出現(xiàn)了由偏枯向偏豐的顯著轉(zhuǎn)變[20-22],與此同時(shí),該流域氣候條件發(fā)生了由“暖干”向“暖濕”的明顯轉(zhuǎn)變,徑流趨勢與氣候變化形成了明顯的響應(yīng)關(guān)系,證明氣候要素(主要是降雨、溫度等)是當(dāng)前階段造成徑流變化的最重要原因[2,14,19]。然而,目前多數(shù)研究仍采用滑動(dòng)平均或多項(xiàng)式、模型模擬、差積曲線分析和小波分析等方法定性擬合水文、氣候變化趨勢、分析氣候特征對(duì)徑流變化的影響[20],而鮮有通過定量分析來探究阿克蘇河徑流對(duì)氣候變化的響應(yīng)情況的研究。【擬解決問題】因此,在全球變暖日益加劇的情況下,利用非參數(shù)Mann-Kendall-Sneyers檢驗(yàn)對(duì)徑流序列進(jìn)行突變點(diǎn)檢驗(yàn),并利用Mann-Kendall趨勢檢驗(yàn)分析氣候變量的變化趨勢,繼而利用基于LI等[17]、ZHENG等[12]所提出的改進(jìn)的基于非參數(shù)的氣候彈性方法定量評(píng)估以溫度為主的氣候要素對(duì)阿克蘇河徑流變化的影響。

    1 材料與方法

    1.1 研究區(qū)概況

    阿克蘇河流域位于天山中段西部南麓地區(qū),塔里木河盆地西北部(75°35′—80°59′E,40°17′—42°27′N),地勢西北高,東南低,垂直地帶分布顯著,地形由南向北、由西向東逐漸下降。流域面積約5×108km2,塔里木河的主要補(bǔ)給來源[20]。該地區(qū)屬于大陸性干旱氣候,氣候干旱,降水稀少,蒸發(fā)量大,日照時(shí)間長,光熱資源豐富[21]。徑流主要來源為大氣降水、冰川和季節(jié)性積雪融水以及少量的地下水補(bǔ)給[22]。年極端最高溫度和最低溫度分別為40.2 ℃和27.6 ℃,年平均日照時(shí)間為2 850 h。庫馬拉克河和托什干河是其2條主要支流,在喀拉都維匯合形成阿克蘇河,山區(qū)是阿克蘇河的主要產(chǎn)流區(qū),復(fù)雜的氣候條件導(dǎo)致了徑流的不規(guī)則波動(dòng)[22]。本研究中選取了流域內(nèi)位于山區(qū)的3個(gè)具代表性的水文站(表1,圖1):沙里桂蘭克、協(xié)和拉、西大橋和各站1956—2011年的實(shí)測徑流資料以及流域內(nèi)自1956—2011年的逐日降雨、日平均氣溫、日最高、最低氣溫、平均相對(duì)濕度、平均風(fēng)速、日照時(shí)間等氣象資料,綜合后分別求算各水文、氣象指標(biāo)的年平均值并基于Penman-Monteith綜合法計(jì)算流域內(nèi)潛在蒸散發(fā)值。

    表1 阿克蘇河流域水文站概況 Table 1 General situation of hydrological stations in Aksu River basin

    圖1 研究區(qū)位置及水文站分布 Fig.1 Location of study area and distribution of hydrological stations

    book=117,ebook=3

    1.2 研究方法

    1.2.1 彈性理論方法

    彈性理論可以用來檢測自變量對(duì)因變量的影響[17]。為了定量評(píng)估氣候因素對(duì)徑流變化的影響,首先假設(shè)徑流變化歸因于氣候因素和非氣候因素相關(guān)的變化(如人類活動(dòng)和其他因素)可以近似如下:

    式中:Q、C和N分別代表徑流,氣候要素(本研究中包括(溫度(T)、降雨(P)和潛在蒸散發(fā)(PET)),和非氣候要素的變化量,f'代表徑流的變化率。

    式中:Ci分別代表各氣候要素(溫度(T)、降雨(P)、潛在蒸散發(fā)(PET))。本研究中假設(shè)各要素對(duì)徑流變化的影響作用相互獨(dú)立,那么式(1)、式(2)可以改寫為:

    式中:QCi表示由要素Ci變化所引起的徑流變化量,該式可用于區(qū)分徑流變化中氣候和非氣候要素的影響。

    根據(jù)Schaake[23]的定義,徑流的氣候彈性可寫為:

    式中:εi為徑流對(duì)氣候要素變化的彈性系數(shù)。基于式(6),可將式(4)、式(5)改寫為:

    遵循以前的研究結(jié)果[17]和上述公式,可依據(jù)下式直接根據(jù)觀測值確定徑流和氣候要素間的關(guān)系:

    式中:Qj和Cij表示年氣候要素值和徑流量(j=1956,1957,…,2011);Qb和Cib是徑流和氣候要素的基準(zhǔn)值(后文將有具體說明)。因此,?i可被看作(Qj-Qb)/Qb和(Cij-Cib)/Cib的線性回歸系數(shù)。為了定量識(shí)別1956—2011年由氣候要素所引起的徑流變化,徑流變化率的計(jì)算式為:

    式中:ΔCi和ΔQ是1956—2011年氣候要素和徑流量的改變量;Φi是氣候要素對(duì)徑流變化的相對(duì)貢獻(xiàn)率。為了驗(yàn)證本研究結(jié)果的合理性,同時(shí)選用了2種具有較高有效性和穩(wěn)定性的基于非參數(shù)的彈性方法來探究年徑流對(duì)氣候要素變化的敏感性[12,24],計(jì)算式:

    式中:為多年平均氣候要素值;Xi為年平均氣候要素值;為多年平均徑流;Qi為年平均徑流。

    1.2.2 突變點(diǎn)監(jiān)測方法

    非參數(shù)Mann-Kendall-Sneyers檢驗(yàn)被廣泛應(yīng)用于確定氣候因子和徑流的階躍變化點(diǎn)的出現(xiàn)[25-26]。x1,?,xn代表數(shù)據(jù)的時(shí)間序列。

    式中:tk為mi的累積序列;mi為序列中滿足后項(xiàng)大于前項(xiàng)條件的累積次數(shù)。在時(shí)間序列服從隨機(jī)獨(dú)立條件假設(shè)的條件下,統(tǒng)計(jì)量的均值和方差計(jì)算式為:

    最后定義統(tǒng)計(jì)量:

    為了獲取統(tǒng)計(jì)量UBK,將徑流數(shù)據(jù)時(shí)間序列按逆序排列,重復(fù)式(12)、式(13)、式(14),規(guī)定UB1=0,在給定顯著水平下,若UFK和UBK曲線在置信區(qū)間中存在交點(diǎn),則該交點(diǎn)即為理論突變點(diǎn)。

    1.2.3 趨勢性檢驗(yàn)

    Mann-Kendall非參數(shù)趨勢檢驗(yàn),由于計(jì)算簡單且不受少數(shù)異常值的干擾,適用于水文和氣象等非正態(tài)分布的數(shù)據(jù)檢測[27]。在本研究中應(yīng)用Mann-Kendall法檢測研究區(qū)內(nèi)徑流、溫度、降雨、潛在蒸散發(fā)的顯著性趨勢。

    2 結(jié)果與分析

    2.1 徑流和氣候要素趨勢性分析

    圖2為1956—2011年阿克蘇河的年徑流量、溫度、降水和潛在蒸散量的變化。Mann-Kendall檢驗(yàn)結(jié)果表明,徑流和溫度分別以11.73 m3/(s·10 a)和0.26 ℃/10 a的速率呈顯著的遞增趨勢(P<0.01)。降雨以14.31 mm/10 a的速率遞增(P<0.05),而潛在蒸散發(fā)卻以21.89 mm/10 a的速率顯著遞減(P<0.01)。

    運(yùn)用Mann-Kendall-Sneyers檢驗(yàn)方法,對(duì)1956—2011年阿克蘇河的徑流量進(jìn)行突變點(diǎn)檢驗(yàn)。年徑流量序列突變點(diǎn)的計(jì)算結(jié)果見圖3。UFk和UBk曲線的交點(diǎn)表明徑流突變點(diǎn)在95%的置信水平上發(fā)生在1993年。該檢驗(yàn)結(jié)果與以前阿克蘇河地區(qū)的研究結(jié)book=0,ebook=4果基本一致[19,21]。隨后為了進(jìn)一步研究氣候要素對(duì)徑流變化的影響,對(duì)降雨、溫度、潛在蒸散發(fā)都進(jìn)行了相同的檢驗(yàn)。檢驗(yàn)結(jié)果表明,溫度(1993年)、降雨(1991年)和潛在蒸散發(fā)(1994年)的突變點(diǎn)都發(fā)生在20世紀(jì)90年代初,也就是說,在1994年之后,由于氣候要素的影響,徑流量發(fā)生了顯著的變化。因此,本研究將1956—1993年的多年平均徑流量和氣候要素作為衡量徑流與氣候要素變化的基準(zhǔn)值。計(jì)算結(jié)果表明,阿克蘇河1994—2011年的徑流量、溫度和降雨分別相對(duì)于1956—1993年增長了23.1%、32.4%和22.5%,而潛在蒸散發(fā)減少了7.6%(表2)。

    表2 阿克蘇河徑流與氣候要素變化情況(1956—2011年) Table 2 Changes of runoff and climate factors in Aksu River (1956—2011)

    圖2 年徑流量與氣候要素變化趨勢 Fig.2 Trends in streamflow and climate factors

    圖3 阿克蘇河徑流M-K-S檢驗(yàn) Fig.3 The M-K-S test of annual streamflow for Aksu River

    2.2 徑流敏感性系數(shù)

    圖4為溫度、降雨、潛在蒸散發(fā)和徑流量之間的線性回歸關(guān)系。阿克蘇河1956—2011年徑流量與溫度、降雨變化率之間顯著正相關(guān)(圖4(a)、圖4(c)),潛在蒸散發(fā)變化率與徑流量之間顯著負(fù)相關(guān)(圖4(b))。由線性回歸方程的系數(shù)即為所求的氣候彈性,由圖4可得到徑流相對(duì)于溫度(εT)、降雨(εP)和潛在蒸散發(fā)(εPET)的彈性分別為1.138、0.270 8和-1.003,這意味著溫度、降雨和潛在蒸散發(fā)每增加1%將導(dǎo)致徑流增加1.138%、0.270 8%和 -1.003%。阿克蘇河徑流變化對(duì)溫度變化的響應(yīng)最為敏感,潛在蒸散發(fā)次之,對(duì)降雨量的變化響應(yīng)最不敏感。此外,由回歸關(guān)系所得的氣候彈性結(jié)果與式(10)、式(11)計(jì)算所得的結(jié)果基本一致(表3),側(cè)面驗(yàn)證了結(jié)果的可靠性。

    2.3 氣候要素對(duì)徑流變化的影響

    基于氣候變化速率和徑流對(duì)降雨量、溫度和潛book=0,ebook=5在蒸散發(fā)變化的彈性系數(shù),利用式(9)可以定量評(píng)估氣候因素對(duì)徑流變化的相對(duì)影響(表3)。由表3可知,溫度的變化對(duì)阿克蘇河1956—2011年來徑流變化的影響最大,推測這應(yīng)該與阿克蘇河上游冰川融雪水的徑流補(bǔ)給比例較大有關(guān);此外,該地區(qū)溫度上升速度為0.26 ℃/10 a,高于同期全國平均增溫速率(0.249 ℃/10 a)[4-5]。特別是相對(duì)于溫度呈下降趨勢的西南地區(qū)[28],溫度對(duì)阿克蘇河徑流變化的影響程度更加劇烈。

    與此同時(shí),由表3、表4可以看出,3種氣候彈性方法(式(8)、式(10)、式(22))的最終計(jì)算結(jié)果基本一致,這意味著本研究所采用的彈性系數(shù)的計(jì)算方法在阿克蘇河地區(qū)的應(yīng)用結(jié)果是基本可信的,本研究的結(jié)論具有一定的合理性和可靠性。以3種計(jì)算方法的均值作為最終的結(jié)果,溫度、降水和潛在蒸散發(fā)對(duì)阿克蘇河徑流變化的相對(duì)貢獻(xiàn)率分別為46.7%、20.5%和24.9%,這表明其他因素(如人類活動(dòng)等)對(duì)徑流變化的相對(duì)貢獻(xiàn)大約為7.9%,說明氣候要素是阿克蘇河徑流變化的主導(dǎo)因素。

    表3 基于不同彈性系數(shù)法的氣候因素對(duì)徑流變化的影響結(jié)果 Table 3 Influence of climate factors on streamflow changes based on different elastic coefficient methods

    表4 基于不同氣候彈性系數(shù)法的氣候要素相對(duì)貢獻(xiàn)率 Table 4 Relative contribution rate of climate elements based on different climate elasticity coefficient methods

    圖4 阿克蘇河1956—2011年徑流變化率與 氣候要素變化率的回歸關(guān)系 Fig.4 The regression relationships between the proportional change of streamflow and climate factors for Aksu River from 1956 to 2011

    book=0,ebook=5

    從目前的研究來看,阿克蘇河流域溫度的變化是徑流變化中最主要的影響要素(46.7%),而非氣候因素(如人類活動(dòng)等)對(duì)徑流變化的影響不大,1956—2011年來徑流的變化只有約7.9%的變化是由非氣候要素變化所引起的。這與該區(qū)域以前的研究基本一致[17,20],但這一結(jié)論存在著一定的矛盾,非氣候要素(主要指人類活動(dòng))與氣候變化之間并不是完全割裂的2類過程,而是存在一定的相互作用關(guān)系[1,7,12]。以人類活動(dòng)為例,人類活動(dòng)的不斷發(fā)展顯著影響了地球上的水文與能量循環(huán)以及水熱平衡,一定程度上促使了氣候變化的發(fā)生與發(fā)展[24,29],反之,氣候的變化也會(huì)影響人類活動(dòng)的舉措與劇烈程book=120,ebook=6度[10-11]。因此,本文中溫度、潛在蒸散發(fā)等氣候要素對(duì)徑流變化的影響中理應(yīng)包含一部分人類活動(dòng)所產(chǎn)生的影響,但本研究主要旨在探究徑流變化對(duì)氣候要素變化的定量響應(yīng)情況,氣候變化與人類活動(dòng)之間的耦合作用的深入探討暫不是本文的研究范圍,尚有待在后續(xù)研究中做進(jìn)一步的分析與討論。

    此外,盡管阿克蘇河上游來水量不斷增加,但塔里木河干流徑流量在1956—2011年來卻呈遞減的傾向,并且人類活動(dòng)是干流徑流量減少的主導(dǎo)因素[1,2,16]。這與本研究的結(jié)果存在一定的沖突,究其原因,造成這種現(xiàn)象的原因主要有以下幾點(diǎn):

    1)阿克蘇河冰川融雪水補(bǔ)給徑流比例較大(約42%),阿克蘇河又處于我國氣溫變化的敏感地帶,溫度的上升加速了上游山區(qū)的冰雪融化,繼而增加了河道的冰雪融水補(bǔ)給。此外,研究區(qū)氣溫呈明顯的上升趨勢,上升速度為0.26 ℃/10 a,高于同期全球平均水平(0.13 ℃/10 a),蔣艷等[22]分析指出,隨著阿克蘇河流域內(nèi)氣溫的逐漸增加,以冰川融雪為主要補(bǔ)給來源的阿克蘇河流域山區(qū)徑流量顯著的上升。

    2)本研究選取的水文站位于山區(qū)和出山口,對(duì)山區(qū)徑流變化上具有較好的代表性,而山區(qū)的人類活動(dòng)強(qiáng)度較弱,因此以人類活動(dòng)為主的非氣候要素對(duì)徑流變化的影響不大。此外,山區(qū)受冰雪融水補(bǔ)給量大,水文環(huán)境較原始,因此對(duì)氣候變化的響應(yīng)情況更為敏感。

    3)人類活動(dòng)對(duì)徑流的影響主要包括直接的水資源開發(fā)利用和對(duì)下墊面條件的間接影響兩方面[30]。阿克蘇河作為新疆重要的棉花和糧食生產(chǎn)基地,人類活動(dòng)對(duì)徑流的影響主要體現(xiàn)在灌溉面積的激增所帶來的水資源量的劇烈消耗[6,16]。但阿克蘇河灌區(qū)主要分布在下游平原地區(qū),對(duì)山區(qū)的開采利用還比較少,隨著人類活動(dòng)不斷地加劇,尤其是阿克蘇河流域灌溉農(nóng)業(yè)的飛速發(fā)展,灌溉面積將不斷向上游擴(kuò)張,流域上下游的耗水量都呈現(xiàn)出顯著的遞增趨勢,人類活動(dòng)將會(huì)對(duì)徑流變化產(chǎn)生越來越重要的影響。此外,由于阿克蘇河下游綠洲農(nóng)業(yè)的擴(kuò)張,以及高耗水作物棉花的廣泛種植,大量消耗了上游增加的來水量[6],從而導(dǎo)致補(bǔ)給干流的水量減少,結(jié)合干流愈發(fā)劇烈的人類活動(dòng)的影響,使得塔里木河干流的水量逐年減少。

    4)由于非參數(shù)方法評(píng)估氣候彈性是假設(shè)每一種氣候要素之間都是相互獨(dú)立互不干擾的,然而,事實(shí)上每種氣候要素或是其他的相關(guān)因素之間都是相互聯(lián)系的,存在著一定的相互作用,因此該類方法仍具有一定的不確定性。同時(shí),觀測值可能受不確定因素的影響進(jìn)而影響計(jì)算結(jié)果,也會(huì)對(duì)結(jié)果造成一定的不確定性影響。如何降低評(píng)估過程的不確定性仍需要在今后的研究中進(jìn)一步探索。

    4 結(jié) 論

    1)徑流和溫度分別以11.73 m3/(s·10 a)和0.26 ℃/10 a的速率顯著遞增(P<0.01)。降雨量以14.31 mm/10 a的速率遞增(P<0.05),而潛在蒸散發(fā)以21.89 mm/10 a的速率呈顯著的遞減趨勢(P<0.01)。突變檢驗(yàn)結(jié)果表明阿克蘇河水文、氣象要素突變點(diǎn)都發(fā)生在20世紀(jì)90年代初。

    2)阿克蘇河徑流變化對(duì)溫度變化的響應(yīng)最為敏感,潛在蒸散發(fā)次之,對(duì)降雨量的變化的敏感性最弱。氣候變化是阿克蘇河徑流變化的主導(dǎo)原因,但隨著流域人類活動(dòng)的加劇,人類活動(dòng)將對(duì)徑流變化產(chǎn)生越來越重要的影響。

    [1] CHEN Y N, LI Z, FAN Y T, et al. Progress and prospects of climate change impacts on hydrology in the arid region of northwest China[J]. Environmental Research, 2015, 139: 11-19.

    [2] 陳亞寧, 徐長春, 郝興明, 等. 新疆塔里木河流域近50a氣候變化及其對(duì)徑流的影響[J]. 冰川凍土, 2008, 30(6): 921-929. CHEN Yaning, XU Changchun, HAO Xingming, et al. Fifty-year climate change and its effect on annual runoff in the Tarim river basin, China[J]. Journal of Glaciology and Geocryology, 2008, 30(6): 921-929.

    [3] LUO M, LIU T, MENG F H, et al. Identifying climate change impacts on water resources in Xinjiang, China[J]. Science of the Total Environment, 2019, 676: 613-626.

    [4] LI B F, CHEN Y N, CHEN Z S, et al. Variations of temperature and precipitation of snowmelt period and its effect on runoff in the mountainous areas of Northwest China[J]. Journal of Geographical Sciences, 2013, 23(1): 17-30.

    [5] 張雪芹, 孫楊, 毛煒嶧, 等. 中國干旱區(qū)氣溫變化對(duì)全球變暖的區(qū)域響應(yīng)[J]. 干旱區(qū)研究, 2010, 27(4): 592-599. ZHANG Xueqin, SUN Yang, MAO Weiyi, et al. Regional response of temperature change in the arid regions of China to global warming[J]. Arid Zone Research, 2010, 27(4): 592-599.

    [6] THEVS N, PENG H Y, ROZI A, et al. Water allocation and water consumption of irrigated agriculture and natural vegetation in the Aksu-Tarim river basin, Xinjiang, China[J]. Journal of Arid Environments, 2015, 112: 87-97.

    [7] WANG X X. Advances in separating effects of climate variability and human activity on stream discharge: An overview[J]. Advances in Water Resources, 2014, 71: 209-218.

    [8] 江善虎, 任立良, 雍斌, 等. 氣候變化和人類活動(dòng)對(duì)老哈河流域徑流的book=121,ebook=7影響[J]. 水資源保護(hù), 2010, 26(6): 1-4, 15. JIANG Shanhu, REN Liliang, YONG Bin, et al. The effects of climate variability and human activities on runoff from the Laohahe River Basin in Northern China[J]. Water Resources Protection, 2010, 26(6): 1-4, 15.

    [9] 黨素珍, 董國濤, 蔣曉輝, 等. 氣候變化和人類活動(dòng)對(duì)涇河上游徑流的影響[J]. 南水北調(diào)與水利科技, 2014, 12(4): 30-34. DANG Suzhen, DONG Guotao, JIANG Xiaohui, et al. Impacts of climate change and human activities on runoff in the upper reaches of Jinghe river[J]. South-to-North Water Transfers and Water Science & Technology, 2014, 12(4): 30-34.

    [10] 張艷霞. 錫林河流域徑流量對(duì)氣候變化與人類活動(dòng)雙重作用的定量響應(yīng)研究[D]. 呼和浩特: 內(nèi)蒙古大學(xué), 2018. ZHANG Yanxia. Quantitative response of runoff in Xilin river basin to the dual effects of climate change and human activities[D]. Hohhot: Inner Mongolia University, 2018.

    [11] 王國慶, 張建云, 鮑振鑫, 等. 人類活動(dòng)和氣候變化對(duì)嵐河流域河川徑流的影響[J]. 灌溉排水學(xué)報(bào), 2019, 38(6): 113-118. WANG Guoqing, ZHANG Jianyun, BAO Zhenxin, et al. Change in runoff within the Lan river basin due to climate change and human activities[J]. Journal of Irrigation and Drainage, 2019, 38(6): 113-118.

    [12] ZHENG H X, ZHANG L, ZHU R R, et al. Responses of streamflow to climate and land surface change in the headwaters of the Yellow River Basin[J]. Water Resources Research, 2009, 45(7): W00A19.

    [13] 陳瑩, 陳興偉. 基于彈性分析的閩江流域徑流演變與氣候變化關(guān)系[J]. 福建師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 27(6): 101-105. CHEN Ying, CHEN Xingwei. Relationship between the changes of runoff and climate in the Min river basin based on elasticity analysis[J]. Journal of Fujian Normal University (Natural Science Edition), 2011, 27(6): 101-105.

    [14] 王志成, 方功煥, 張輝, 等. 基于高空與地面觀測的阿克蘇河流域氣候水文要素變化分析[J]. 氣候變化研究進(jìn)展, 2018, 14(1): 1-10. WANG Zhicheng, FANG Gonghuan, ZHANG Hui, et al. Analysis of hydrometeorological variations based on the sounding and near-surface observations in Aksu River Basin[J]. Climate Change Research, 2018, 14(1): 1-10.

    [15] SANKARASUBRAMANIAN A, VOGEL R M, LIMBRUNNER J F. Climate elasticity of streamflow in the United States[J]. Water Resources Research, 2001, 37(6): 1 771-1 781.

    [16] YU Y, DE YU R, CHEN X, et al. Agricultural water allocation strategies along the oasis of Tarim River in Northwest China[J]. Agricultural Water Management, 2017, 187: 24-36.

    [17] LI B F, CHEN Y N, XIONG H G. Quantitatively evaluating the effects of climate factors on runoff change for Aksu River in northwestern China[J]. Theoretical and Applied Climatology, 2016, 123(1/2): 97-105.

    [18] 茍樹生, 張雄文, 王彥國, 等. 近50年來塔里木河干流水量、水質(zhì)及耗水分析[J]. 干旱區(qū)研究, 2010, 27(6): 861-870. GOU Shusheng, ZHANG Xiongwen, WANG Yanguo, et al. Analysis on runoff volumes, water quality and water consumption of the Tarim river in recent 50 years[J]. Arid Zone Research, 2010, 27(6): 861-870.

    [19] 黃國強(qiáng). 近50年來阿克蘇河徑流演變特征分析[J]. 地下水, 2019, 41(5): 132-135. HUANG Guoqiang. Analysis on the evolution characteristics of runoff in the Aksu river in the past 50 years[J]. Ground Water, 2019, 41(5): 132-135.

    [20] 柏玲, 陳忠升, 王充, 等. 西北干旱區(qū)阿克蘇河徑流對(duì)氣候波動(dòng)的多尺度響應(yīng)[J]. 地理科學(xué), 2017, 37(5): 799-806. BAI Ling, CHEN Zhongsheng, WANG Chong, et al. Multi-scale response of runoff to climate change in the Aksu river basin in arid area of the northwestern China[J]. Scientia Geographica Sinica, 2017, 37(5): 799-806.

    [21] FAN Y T, CHEN Y N, LI W H. Increasing precipitation and baseflow in Aksu River since the 1950s[J]. Quaternary International, 2014, 336: 26-34.

    [22] 蔣艷, 周成虎, 程維明. 阿克蘇河流域徑流補(bǔ)給及徑流變化特征分析[J]. 自然資源學(xué)報(bào), 2005, 20(1): 27-34. JIANG Yan, ZHOU Chenghu, CHENG Weiming. Analysis on runoff supply and variation characteristics of Aksu drainage basin[J]. Journal of Natural Resources, 2005, 20(1): 27-34.

    [23] SCHAAKE J C. From climate to flow in Climate Change and U.S. Water Resources[M]. New York: John Wiley, 1990:177-206.

    [24] AHN K H, MERWADE V. Quantifying the relative impact of climate and human activities on streamflow[J]. Journal of Hydrology, 2014, 515: 257-266.

    [25] 張應(yīng)華, 宋獻(xiàn)方. 水文氣象序列趨勢分析與變異診斷的方法及其對(duì)比[J]. 干旱區(qū)地理, 2015, 38(4): 652-665. ZHANG Yinghua, SONG Xianfang. Techniques of abrupt change detection and trends analysis in hydroclimatic time-series: Advances and evaluation[J]. Arid Land Geography, 2015, 38(4): 652-665.

    [26] 陳俊旭, 趙紅玲, 呂燕, 等. 紅河流域徑流演變分異及其特征[J]. 水電能源科學(xué), 2018, 36(10): 22-26. CHEN Junxu, ZHAO Hongling, LYU Yan, et al. Study on the runoff variation and its characteristics of red river basin[J]. Water Resources and Power, 2018, 36(10): 22-26.

    [27] 彭兆亮, 何斌, 王國利, 等. 碧流河流域氣候變化對(duì)水資源量的影響研究[J]. 南水北調(diào)與水利科技, 2010, 8(5): 76-79. PENG Zhaoliang, HE Bin, WANG Guoli, et al. A study on characteristics of climate change and its influence on water resources in resent 30 years in Biliuhe river basin[J]. South-to-North Water Transfers and Water Science & Technology, 2010, 8(5): 76-79.

    [28] 任國玉, 徐銘志, 初子瑩, 等. 近54年中國地面氣溫變化[J]. 氣候與環(huán)境研究, 2005, 10(4): 717-727. REN Guoyu, XU Mingzhi, CHU Ziying, et al. Changes of surface air temperature in China during 1951—2004[J]. Climatic and Environmental Research, 2005, 10(4): 717-727.

    book=122,ebook=8[29] 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:362-375.

    [30] 劉啟寧, 辛卓航, 韓建旭, 等. 變化環(huán)境下東北半干旱地區(qū)徑流演變規(guī)律分析:以洮兒河流域?yàn)槔齕J]. 水力發(fā)電學(xué)報(bào), 2020, 39(5): 51-63. LIU Qining, XIN Zhuohang, HAN Jianxu, et al. Analysis on evolving trends of streamflow in the semi-arid region in Northeast China under changing environment. A case study of Taoer River basin[J]. Journal of Hydroelectric Engineering, 2020, 39(5): 51-63.

    Assessing the Impact of Meteorological Factors on Streamflow in Aksu River

    LI Hongbin1, LIU Yating1, WANG Weiguang1,2*, PANG Maitian1

    (1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China; 2.State Key Laboratory of Hydrology-water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China)

    【Background】The effects of potential climate change on water flow in rivers have been well documented but with the focus putting on precipitation and potential evapotranspiration (PET). There is a lack of study on the impact of individual meteorological factors.【Objective】Taking Aksu river in Xinjiang as an example, this paper investigated the relationship between water flow rate in the river and some meteorological factors.【Method】The trends and abrupt changes in hydrological and meteorological data measured from the basin were analyzed using the Mann-Kendall-Sneyers test and Mann-Kendall nonparametric trend test; three improved elastic coefficient methods were used to evaluate the influence of the meteorological factors, particularly temperature, on water flow rate in the river.【Result】Both water flow rate in river and temperature in the basin had increased at rates of 11.73 (m3/s)/10a and 0.26 ℃/10a, respectively, at significant level (P<0.01). The precipitation trended up at a rate of 14.31 mm/10a at significant level (P<0.05), while the potential evapotranspiration had been in decline at a rate of 21.89 mm/10a at significant level (P<0.01). The test results showed that the temperature changed abruptly in 1993, the precipitation in 1991 and the potential evapotranspiration in the early 1990s. Water flow rate in the river underwent an abrupt change in 1993, indicating that water in the river has changed significantly since 1994 due to climatic impact. A multi-year average water flow rate and climate factors from 1960 to 1993 were thus used as the baseline to measure the change in water flow after 1993. From 1994 to 2010, streamflow, precipitation and temperature increased by 23.1%, 22.5% and 32.4% respectively, while the potential evapotranspiration reduced by 7.6%, compared to the averages over 1960 to 1993. The linear relationship between temperature, precipitation, potential evapotranspiration and streamflow passed the statistical test at significant level of 1%, indicating that the linear regression coefficients reflect the relationship between surface runoff and climate factors and that it can be used to deduce the climate elasticity. Sensitivity analysis showed that the elasticity of the streamflow to the temperature (εT) and the precipitation (εP) was 1.138 and 0.270 8 respectively, meaning that every 1% increase in the temperature and precipitation will lead to 1.138% and 0.270 8% increase in streamflow respectively. The elasticity coefficient of the streamflow to the potential evapotranspiration was -1.003(εPET), indicating that a 1% decrease in potential evapotranspiration would lead to a 1.003% rise in streamflow. Water flow in the river is most sensitive to temperature and least to precipitation, with the potential evapotranspiration between. It was also found that the relative contribution of temperature to streamflow was 46.7%, while the contribution of precipitation and potential evapotranspiration was 20.5% and 24.9% respectively. Non-climatic factors, such as anthropogenic activities, contributed only 7.9% to the streamflow change.【Conclusion】Our analysis showed that climate change is currently the main factor behind the streamflow change in Aksu river, despite the increasing influence of the continuous intensified anthropogenic activities at upper and lower reaches of the basin.

    streamflow variation; Mann-Kendall-Sneyers test; elastic coefficient; climate change

    P333.1

    A

    10.13522/j.cnki.ggps.2020287

    1672 - 3317(2021)01 - 0115 - 08

    2020-05-30

    國家自然科學(xué)基金項(xiàng)目(51979071);國家級(jí)大學(xué)生創(chuàng)新訓(xùn)練項(xiàng)目(202010294009)

    李虹彬(1998-),男。學(xué)士,主要從事水文物理規(guī)律模擬研究。 E-mail: lihongbin_hhu@163.com

    王衛(wèi)光(1979-),男。教授,博士生導(dǎo)師,主要從事水文物理規(guī)律模擬研究。E-mail: wangweiguang2016@126.com

    李虹彬, 劉亞婷, 王衛(wèi)光, 等. 氣候因素對(duì)阿克蘇河徑流變化影響的定量評(píng)估[J]. 灌溉排水學(xué)報(bào), 2021, 40(1): 115-122.

    LI Hongbin, LIU Yating, WANG Weiguang, et al. Assessing the Impact of Meteorological Factors on Streamflow in Aksu River[J]. Journal of Irrigation and Drainage, 2021, 40(1): 115-122.

    責(zé)任編輯:趙宇龍

    猜你喜歡
    徑流量水文徑流
    2022年《中國水文年報(bào)》發(fā)布
    水文
    水文水資源管理
    水文比擬法在計(jì)算河川徑流量時(shí)的修正
    水文
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    探秘“大徑流”
    攻克“大徑流”
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    資江流域徑流量演變規(guī)律研究
    18禁黄网站禁片午夜丰满| 天堂网av新在线| www日本黄色视频网| 麻豆av噜噜一区二区三区| 最后的刺客免费高清国语| 国产成人福利小说| www.色视频.com| 成人毛片a级毛片在线播放| 少妇熟女aⅴ在线视频| 国产单亲对白刺激| 少妇裸体淫交视频免费看高清| 日本 av在线| 中文亚洲av片在线观看爽| 欧美极品一区二区三区四区| 97超视频在线观看视频| 日本五十路高清| 九九在线视频观看精品| 成人高潮视频无遮挡免费网站| 国产高清三级在线| 性欧美人与动物交配| 成人无遮挡网站| 麻豆国产97在线/欧美| 嫩草影视91久久| 国产精品乱码一区二三区的特点| 啦啦啦韩国在线观看视频| 午夜激情欧美在线| 色综合亚洲欧美另类图片| 亚洲真实伦在线观看| 国产高潮美女av| 内射极品少妇av片p| 久久伊人香网站| 亚洲精品久久国产高清桃花| 淫秽高清视频在线观看| 丰满的人妻完整版| 深夜a级毛片| 国产老妇女一区| 高清在线国产一区| 99久久精品国产亚洲精品| 亚洲五月天丁香| 亚洲无线在线观看| www.色视频.com| 午夜免费成人在线视频| 日韩中文字幕欧美一区二区| 国内精品一区二区在线观看| 波野结衣二区三区在线| 小蜜桃在线观看免费完整版高清| 美女 人体艺术 gogo| 女人十人毛片免费观看3o分钟| 亚洲av一区综合| 亚洲av成人av| 精品人妻视频免费看| 亚洲专区中文字幕在线| 亚洲人成网站在线播放欧美日韩| 国产亚洲精品综合一区在线观看| 亚洲人成网站高清观看| 久久久久精品国产欧美久久久| 成人鲁丝片一二三区免费| 日本黄色片子视频| xxxwww97欧美| 久9热在线精品视频| 一本一本综合久久| 美女xxoo啪啪120秒动态图 | 亚洲内射少妇av| 日韩欧美一区二区三区在线观看| 亚洲真实伦在线观看| 成年女人看的毛片在线观看| 免费无遮挡裸体视频| a级一级毛片免费在线观看| 麻豆国产av国片精品| 中国美女看黄片| а√天堂www在线а√下载| a级毛片免费高清观看在线播放| 国产国拍精品亚洲av在线观看| 亚洲成人免费电影在线观看| 亚洲中文字幕一区二区三区有码在线看| 在线a可以看的网站| 成人高潮视频无遮挡免费网站| .国产精品久久| 99久久99久久久精品蜜桃| 首页视频小说图片口味搜索| 一本久久中文字幕| 久久天躁狠狠躁夜夜2o2o| av视频在线观看入口| 69人妻影院| 日本一二三区视频观看| 欧美成人一区二区免费高清观看| 丝袜美腿在线中文| 亚洲黑人精品在线| 性色av乱码一区二区三区2| 久久亚洲真实| 久久久久国内视频| 色哟哟·www| 精品一区二区免费观看| 最新中文字幕久久久久| 国产三级黄色录像| 国产精品久久视频播放| 男插女下体视频免费在线播放| 嫩草影院新地址| 免费高清视频大片| 少妇丰满av| 热99在线观看视频| 不卡一级毛片| 欧美日韩黄片免| 一区二区三区激情视频| 精品一区二区三区视频在线观看免费| 亚洲电影在线观看av| 亚洲精品在线美女| 欧美最新免费一区二区三区 | 国产精品98久久久久久宅男小说| 亚洲 欧美 日韩 在线 免费| 亚洲人成网站在线播放欧美日韩| 99国产极品粉嫩在线观看| 亚洲专区中文字幕在线| 大型黄色视频在线免费观看| 少妇熟女aⅴ在线视频| 国产aⅴ精品一区二区三区波| 欧美极品一区二区三区四区| 如何舔出高潮| 亚洲第一欧美日韩一区二区三区| 亚洲avbb在线观看| 国产aⅴ精品一区二区三区波| 啦啦啦韩国在线观看视频| 在线观看免费视频日本深夜| 91久久精品电影网| 精品国产亚洲在线| 亚洲电影在线观看av| 亚洲一区二区三区色噜噜| 91午夜精品亚洲一区二区三区 | 国内精品一区二区在线观看| 伊人久久精品亚洲午夜| 99精品在免费线老司机午夜| 淫妇啪啪啪对白视频| 午夜久久久久精精品| 国产久久久一区二区三区| 免费人成在线观看视频色| 国产av麻豆久久久久久久| 久久久精品大字幕| 看免费av毛片| 成人欧美大片| 欧美乱色亚洲激情| 免费一级毛片在线播放高清视频| 精品久久久久久久久久久久久| 国产成人欧美在线观看| 精品欧美国产一区二区三| 在线免费观看不下载黄p国产 | 九九热线精品视视频播放| 成年女人永久免费观看视频| 国产精品精品国产色婷婷| 听说在线观看完整版免费高清| 最近中文字幕高清免费大全6 | 国产亚洲精品久久久久久毛片| 91狼人影院| 午夜久久久久精精品| 又爽又黄无遮挡网站| 蜜桃久久精品国产亚洲av| 亚洲在线观看片| 窝窝影院91人妻| 成人高潮视频无遮挡免费网站| 日本黄色视频三级网站网址| 日本黄色视频三级网站网址| 久久热精品热| 在线观看免费视频日本深夜| 国产精品国产高清国产av| 欧美绝顶高潮抽搐喷水| 网址你懂的国产日韩在线| 亚洲va日本ⅴa欧美va伊人久久| 午夜视频国产福利| 18禁黄网站禁片免费观看直播| 久久久久亚洲av毛片大全| 中文字幕av成人在线电影| 一区福利在线观看| 长腿黑丝高跟| 免费大片18禁| 亚洲av第一区精品v没综合| 两个人视频免费观看高清| 亚洲aⅴ乱码一区二区在线播放| 免费搜索国产男女视频| 99热只有精品国产| 制服丝袜大香蕉在线| 青草久久国产| 国产精品嫩草影院av在线观看 | 他把我摸到了高潮在线观看| 有码 亚洲区| 久久国产乱子伦精品免费另类| 老熟妇仑乱视频hdxx| 亚洲经典国产精华液单 | 欧美性感艳星| 亚洲性夜色夜夜综合| 亚洲aⅴ乱码一区二区在线播放| 成年人黄色毛片网站| 国产精品嫩草影院av在线观看 | 亚洲片人在线观看| 日本黄色视频三级网站网址| 欧美精品啪啪一区二区三区| 亚洲最大成人中文| 国产精品98久久久久久宅男小说| 成年女人看的毛片在线观看| 日韩欧美国产一区二区入口| 日韩人妻高清精品专区| 亚洲美女黄片视频| 自拍偷自拍亚洲精品老妇| 久久久久久久久大av| 99热只有精品国产| 一级毛片久久久久久久久女| 成人无遮挡网站| 最后的刺客免费高清国语| 99久久无色码亚洲精品果冻| a级毛片免费高清观看在线播放| 久久精品夜夜夜夜夜久久蜜豆| 搡老岳熟女国产| 黄片小视频在线播放| 久久久精品欧美日韩精品| 欧美日韩中文字幕国产精品一区二区三区| 99在线人妻在线中文字幕| 俄罗斯特黄特色一大片| 欧美xxxx黑人xx丫x性爽| 欧美日韩中文字幕国产精品一区二区三区| 国产一区二区三区在线臀色熟女| 最近中文字幕高清免费大全6 | 日本五十路高清| 久久国产乱子伦精品免费另类| av中文乱码字幕在线| 精品人妻视频免费看| www.色视频.com| 国产欧美日韩一区二区三| 女生性感内裤真人,穿戴方法视频| 午夜福利免费观看在线| 噜噜噜噜噜久久久久久91| 国产高清三级在线| 亚洲av免费高清在线观看| www.www免费av| 日本一二三区视频观看| 一级黄片播放器| 国产精品美女特级片免费视频播放器| 精品久久久久久久久av| 免费电影在线观看免费观看| 日韩欧美三级三区| 国产精品精品国产色婷婷| 嫩草影院入口| 久久99热6这里只有精品| 久久久久精品国产欧美久久久| 国产在线男女| 波野结衣二区三区在线| 91午夜精品亚洲一区二区三区 | 色综合婷婷激情| 99久久99久久久精品蜜桃| 欧美+日韩+精品| 熟妇人妻久久中文字幕3abv| .国产精品久久| 一进一出抽搐gif免费好疼| 国产免费一级a男人的天堂| 国产精品影院久久| av国产免费在线观看| 日韩精品青青久久久久久| 成人国产一区最新在线观看| 国内少妇人妻偷人精品xxx网站| 欧美日本视频| 成人精品一区二区免费| 老司机福利观看| 精品无人区乱码1区二区| 性欧美人与动物交配| 国产精品伦人一区二区| 中文字幕久久专区| 人妻制服诱惑在线中文字幕| 亚洲成人免费电影在线观看| 午夜福利欧美成人| 极品教师在线免费播放| 真人一进一出gif抽搐免费| 国产色爽女视频免费观看| 人人妻人人看人人澡| 好看av亚洲va欧美ⅴa在| 少妇熟女aⅴ在线视频| 日韩免费av在线播放| 欧美黄色片欧美黄色片| 成人午夜高清在线视频| 亚洲精华国产精华精| 久久久久性生活片| 亚洲人成网站在线播| 在线国产一区二区在线| av黄色大香蕉| 别揉我奶头 嗯啊视频| 日本在线视频免费播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区二区三区不卡视频| 三级国产精品欧美在线观看| 日韩人妻高清精品专区| 日韩欧美三级三区| 午夜老司机福利剧场| 精品日产1卡2卡| 国产在线男女| 婷婷精品国产亚洲av在线| 国产主播在线观看一区二区| 免费av毛片视频| 欧美性猛交黑人性爽| 午夜福利在线观看吧| 三级毛片av免费| 国产久久久一区二区三区| 国产欧美日韩精品一区二区| 婷婷亚洲欧美| 日韩欧美国产一区二区入口| 在线观看舔阴道视频| 亚洲三级黄色毛片| av福利片在线观看| 亚洲精品亚洲一区二区| 国产精品一及| 一级毛片久久久久久久久女| 色播亚洲综合网| 国产亚洲av嫩草精品影院| 欧美色欧美亚洲另类二区| 国产淫片久久久久久久久 | 久久久久性生活片| 99久久精品国产亚洲精品| 免费看日本二区| 亚洲美女搞黄在线观看 | 国产高清视频在线观看网站| 两个人视频免费观看高清| 女人被狂操c到高潮| 国产午夜精品论理片| 12—13女人毛片做爰片一| 国产主播在线观看一区二区| 色哟哟·www| 男插女下体视频免费在线播放| 最新在线观看一区二区三区| 国产成人a区在线观看| 欧美激情国产日韩精品一区| 一区二区三区高清视频在线| 欧美+日韩+精品| 亚洲精品在线观看二区| 午夜福利免费观看在线| 一进一出抽搐动态| 国产色婷婷99| 波野结衣二区三区在线| 2021天堂中文幕一二区在线观| 欧美性猛交黑人性爽| 嫩草影院入口| 欧美成狂野欧美在线观看| 成人鲁丝片一二三区免费| 日日夜夜操网爽| 悠悠久久av| 99久久精品热视频| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 日日夜夜操网爽| 日韩人妻高清精品专区| 亚洲精品粉嫩美女一区| 熟女人妻精品中文字幕| 天堂av国产一区二区熟女人妻| 一区二区三区高清视频在线| 欧美精品国产亚洲| 亚洲成av人片免费观看| 最近最新中文字幕大全电影3| 欧美高清成人免费视频www| 日本免费一区二区三区高清不卡| 国产av不卡久久| 国产麻豆成人av免费视频| 国产欧美日韩精品一区二区| 精品人妻偷拍中文字幕| 亚洲精品粉嫩美女一区| 少妇人妻精品综合一区二区 | 久久人人爽人人爽人人片va | 欧美国产日韩亚洲一区| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| av视频在线观看入口| 精品无人区乱码1区二区| 中文字幕高清在线视频| av在线蜜桃| 午夜福利在线在线| 亚洲无线在线观看| 日韩大尺度精品在线看网址| 亚洲欧美日韩高清在线视频| 欧美激情国产日韩精品一区| 精品人妻偷拍中文字幕| 精品午夜福利在线看| 欧美色视频一区免费| 看免费av毛片| 每晚都被弄得嗷嗷叫到高潮| 天堂√8在线中文| 亚洲人与动物交配视频| 欧美激情国产日韩精品一区| 久久精品国产自在天天线| 一区二区三区免费毛片| 国产在视频线在精品| av福利片在线观看| 亚洲乱码一区二区免费版| 俄罗斯特黄特色一大片| 精品欧美国产一区二区三| 成人精品一区二区免费| 久久精品人妻少妇| 噜噜噜噜噜久久久久久91| 久久99热6这里只有精品| 成人高潮视频无遮挡免费网站| 免费电影在线观看免费观看| 热99re8久久精品国产| 人人妻人人看人人澡| 久久国产乱子伦精品免费另类| 极品教师在线免费播放| 成人亚洲精品av一区二区| 欧美最新免费一区二区三区 | 俄罗斯特黄特色一大片| 精品99又大又爽又粗少妇毛片 | 国产成+人综合+亚洲专区| 神马国产精品三级电影在线观看| 69人妻影院| 又爽又黄无遮挡网站| 色5月婷婷丁香| 日韩欧美国产一区二区入口| 噜噜噜噜噜久久久久久91| 成人av在线播放网站| 亚洲av中文字字幕乱码综合| 深夜a级毛片| 最后的刺客免费高清国语| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 麻豆国产97在线/欧美| 首页视频小说图片口味搜索| 99久久九九国产精品国产免费| 国产精品一及| bbb黄色大片| 天天一区二区日本电影三级| 国产私拍福利视频在线观看| 99在线视频只有这里精品首页| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头 嗯啊视频| 亚洲第一区二区三区不卡| 亚洲av美国av| 亚洲av日韩精品久久久久久密| a在线观看视频网站| 看免费av毛片| 欧美又色又爽又黄视频| avwww免费| 最近在线观看免费完整版| 亚洲无线观看免费| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 国产精品野战在线观看| 每晚都被弄得嗷嗷叫到高潮| 嫩草影院精品99| 国产v大片淫在线免费观看| 国产精品久久久久久久电影| 舔av片在线| 女人被狂操c到高潮| 亚洲第一区二区三区不卡| 成人毛片a级毛片在线播放| 久久久色成人| 成人鲁丝片一二三区免费| 又紧又爽又黄一区二区| a级毛片a级免费在线| 黄色丝袜av网址大全| 国产野战对白在线观看| 给我免费播放毛片高清在线观看| 九色国产91popny在线| 嫩草影视91久久| 午夜福利在线观看吧| av专区在线播放| 中文字幕人成人乱码亚洲影| 日本一二三区视频观看| 特级一级黄色大片| 精品乱码久久久久久99久播| 色视频www国产| 国产一级毛片七仙女欲春2| 欧美性猛交╳xxx乱大交人| 精品人妻一区二区三区麻豆 | 久久精品国产亚洲av香蕉五月| 麻豆国产av国片精品| 国产精品乱码一区二三区的特点| 欧美黑人欧美精品刺激| 国产69精品久久久久777片| 免费在线观看影片大全网站| 91久久精品国产一区二区成人| 中文资源天堂在线| 国产成人福利小说| 一卡2卡三卡四卡精品乱码亚洲| 精品人妻偷拍中文字幕| 午夜福利免费观看在线| 成年女人毛片免费观看观看9| 久久久久久久久中文| 2021天堂中文幕一二区在线观| 一区二区三区免费毛片| 中文字幕高清在线视频| 99久久精品热视频| 亚洲成a人片在线一区二区| 日韩中文字幕欧美一区二区| 啪啪无遮挡十八禁网站| 久久久久久久久中文| 精品午夜福利在线看| 亚洲人成网站高清观看| 亚洲欧美日韩无卡精品| 男人和女人高潮做爰伦理| 亚洲av二区三区四区| 久久精品综合一区二区三区| 亚洲久久久久久中文字幕| 欧美又色又爽又黄视频| 午夜福利视频1000在线观看| 欧美乱色亚洲激情| 噜噜噜噜噜久久久久久91| 日本一本二区三区精品| 男女下面进入的视频免费午夜| 一进一出抽搐动态| 亚洲中文字幕一区二区三区有码在线看| 69人妻影院| 免费av观看视频| 亚洲熟妇中文字幕五十中出| 亚洲不卡免费看| 大型黄色视频在线免费观看| 日韩欧美精品免费久久 | 亚洲av美国av| 真实男女啪啪啪动态图| 欧美日韩综合久久久久久 | 久久久久久九九精品二区国产| 一区二区三区免费毛片| 亚洲最大成人中文| 偷拍熟女少妇极品色| 18+在线观看网站| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 亚洲av成人不卡在线观看播放网| 久久九九热精品免费| 熟女电影av网| 日本 欧美在线| 看十八女毛片水多多多| av视频在线观看入口| 国产精品亚洲美女久久久| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久99热这里只有精品18| 婷婷丁香在线五月| 国产午夜精品久久久久久一区二区三区 | 久久久久久久久久黄片| 尤物成人国产欧美一区二区三区| 丁香欧美五月| 51国产日韩欧美| 日韩免费av在线播放| 1024手机看黄色片| 在线观看66精品国产| 欧美成人a在线观看| 亚洲av五月六月丁香网| 成人午夜高清在线视频| 亚洲最大成人中文| 亚洲av熟女| 又黄又爽又免费观看的视频| 精品免费久久久久久久清纯| 欧美绝顶高潮抽搐喷水| 亚洲内射少妇av| 看片在线看免费视频| 国产精品,欧美在线| 在线播放国产精品三级| 欧美bdsm另类| 欧美色欧美亚洲另类二区| av在线观看视频网站免费| 亚洲人成伊人成综合网2020| 成人毛片a级毛片在线播放| 最近在线观看免费完整版| 国产精品一区二区性色av| 色播亚洲综合网| 亚洲成av人片在线播放无| 丰满人妻一区二区三区视频av| 久久精品国产99精品国产亚洲性色| 中文字幕av成人在线电影| 最近在线观看免费完整版| 免费在线观看影片大全网站| 欧美成人性av电影在线观看| 久久久久久久亚洲中文字幕 | 国模一区二区三区四区视频| 免费在线观看日本一区| 国产黄a三级三级三级人| 日日夜夜操网爽| 99在线人妻在线中文字幕| 中文字幕久久专区| 国产一区二区在线av高清观看| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 真人做人爱边吃奶动态| 欧美日韩中文字幕国产精品一区二区三区| 麻豆一二三区av精品| 舔av片在线| 亚洲欧美日韩东京热| а√天堂www在线а√下载| 人妻制服诱惑在线中文字幕| 欧美zozozo另类| 国产乱人伦免费视频| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 国内精品一区二区在线观看| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片 | 精品久久久久久久久av| 日日摸夜夜添夜夜添小说| 露出奶头的视频| 好看av亚洲va欧美ⅴa在| 少妇裸体淫交视频免费看高清| 少妇熟女aⅴ在线视频| 亚洲精华国产精华精| 亚洲中文字幕一区二区三区有码在线看| 国产中年淑女户外野战色| 久久国产乱子免费精品| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影| 精品久久久久久久久久免费视频| 九色成人免费人妻av| 91狼人影院| 欧美激情在线99| 动漫黄色视频在线观看| 简卡轻食公司| 亚洲美女黄片视频| 99热这里只有精品一区| 丰满人妻一区二区三区视频av| 脱女人内裤的视频| 国内毛片毛片毛片毛片毛片| 又黄又爽又免费观看的视频| 久久午夜福利片|