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

    1960-2013年廣東省梅江流域徑流特征及其對降水的滯后效應

    2016-07-26 11:31:01董才文張正棟楊傳訓萬露文
    水土保持通報 2016年3期
    關鍵詞:滯后效應集中度徑流

    董才文, 張正棟, 楊傳訓, 萬露文, 張 倩

    (華南師范大學 地理科學學院, 廣東 廣州 510631)

    ?

    1960-2013年廣東省梅江流域徑流特征及其對降水的滯后效應

    董才文, 張正棟, 楊傳訓, 萬露文, 張 倩

    (華南師范大學 地理科學學院, 廣東 廣州 510631)

    摘要:[目的] 全面了解梅江流域的水文特征及其變化趨勢,為河流水量調節(jié)、防洪防汛及流域生態(tài)環(huán)境保護等方面提供決策依據(jù)。[方法] 以該地區(qū)的典型河流梅江為例,利用1960—2013年流域的徑流、降水數(shù)據(jù),結合累積距平、變異系數(shù)、集中度和集中期等數(shù)理方法,對梅江的徑流量、降水量在年代、年際、季節(jié)等時間尺度上進行分析,并用集中期來反映徑流對降水變化的滯后效應。[結果] 梅江徑流量年際變化整體呈下降趨勢,遞減率為5.79 mm/10 a;流域徑流主要集中在4—9月,季節(jié)上主要集中在夏季和春季,徑流集中期主要集中在5月下旬至7月中旬;流域內多年降水量呈波動減少趨勢,但減少幅度很小,降水集中期主要集中在5—6月;流域內年徑流量與降水量存在極大的相關性,并呈線性回歸關系(r=0.824>0.428=R0.001,54,回歸顯著水平超過0.001),流域徑流對降水具有滯后效應,多年平均滯后時間為21 d,且滯后天數(shù)隨著時間推移呈擴大趨勢。[結論] 可以推測是人類活動改變了地表覆被狀況或者是氣候的變化。

    關鍵詞:徑流; 降水; 滯后效應; 集中度; 集中期; 梅江流域

    文獻參數(shù): 董才文, 張正棟, 楊傳訓, 等.1960—2013年廣東省梅江流域徑流特征及其對降水的滯后效應[J].水土保持通報,2016,36(3):38-45.DOI:10.13961/j.cnki.stbctb.2016.03.008

    水資源是人類生產與生活活動的基礎自然資源,是生態(tài)環(huán)境建設的控制性因素,同時也是具有戰(zhàn)略性的經濟資源,是一個國家綜合國力的有機組成部分[1]。作為水文循環(huán)的重要環(huán)節(jié),河川徑流是水資源綜合開發(fā)利用、科學管理及優(yōu)化調度的重要依據(jù)。近年來,隨著全球變暖和人類活動影響的加劇,河川徑流發(fā)生了顯著的時空變化,出現(xiàn)了一系列生態(tài)環(huán)境問題:流域旱澇災害頻發(fā)、水土流失嚴重、水質惡化等,這直接影響了社會經濟的可持續(xù)發(fā)展,全球氣候變化是導致徑流發(fā)生變化的最直接原因,因此研究氣候變化與徑流的關系具有重要意義。

    國內外許多學者對這兩者的關系進行了研究,如劉昌明等[2]利用黃河干流控制流域匯水的主要水文站近40 a的逐月徑流資料,分成了4個系列進行時序分析,得出了黃河下游徑流主要來源于上中游入水,以及中下游基流量在39 a中有下降趨勢這些非常有意義的結論。馬振梅等[3]通過對西北干早地區(qū)石羊河流域氣候變化和人類活動對于徑流的影響研究表明,流域大部分集水區(qū)徑流具有明顯的減少趨勢,氣候變化對徑流變化影響的比重占64%~87%。王兆禮等[4]探討了東江流域徑流年際變化特征及其對氣候變化和植被覆蓋變化的響應,結果表明氣候變化、厄爾尼諾現(xiàn)象、太陽黑子數(shù)和NDVI都對徑流變化有一定的影響,但導致東江流域徑流變化的根本原因是氣候因子的綜合作用。凌紅波等[5]利用葉爾羌河源流區(qū)近50 a的逐月降水量與徑流量資料,借助非參數(shù)檢驗、小波分析等方法研究了葉爾羌河年徑流量與降水量的集中性及多時間尺度相關關系。國外對于這方面的研究多傾向于采用大中尺度水文模型來研究徑流量和氣候變化的關系,如Nash等[6]采用一個修正的水平衡模型研究了科羅拉多河水文要素系統(tǒng),發(fā)現(xiàn)降水是影響徑流變化最主要的因素。Milly等[7]用多個GCM集合模擬1990—1998年全球大尺度徑流分布,并指出氣候脅迫已經對20世紀70年代以后全球徑流分布產生影響,且氣候強迫信號在歐亞大陸北部和北美洲西北部的高緯地區(qū)比較顯著;Githui等[8]通過設置SWAT模型不同情景模擬,預測了肯尼亞恩佐亞流域氣候變化的徑流響應,發(fā)現(xiàn)徑流對溫度變化并不敏感,但如果未來年均降水量增加2.4%~23.2%,年徑流將增加6%~11.5%,表現(xiàn)出良好的降水—徑流線性關系。

    梅江發(fā)源于汕尾陸豐與河源紫金交界的烏突山,是廣東省沿海河流韓江上游干流段(以梅江為韓江主源)。作為典型的華南濕熱地區(qū)的河流,梅江的徑流全由降水補給,河流水位、徑流量及泥沙含量等水文特征及流域周邊環(huán)境均受氣候變化的影響。河流支流多,流域范圍廣,汛期長,年徑流量大。受季風環(huán)流、鋒面雨和臺風雨的共同影響,流域降水和徑流的年內變化和年際變化較大,流域旱澇災害頻發(fā),影響了社會經濟的發(fā)展。目前,對梅江流域的研究較少,現(xiàn)有的工作主要集中于水質研究[9-10]和水文特性的簡單分析[11],因而有必要對流域的徑流量和降水量的特征和規(guī)律進行綜合分析,并找出二者的關系。

    本文以梅江流域為研究區(qū)域,基于流域中游1960—2013年54 a的徑流量實測數(shù)據(jù)和降水量數(shù)據(jù),采用變異系數(shù)、累積距平、集中度和集中期等方法,對流域內的徑流量、降水量在多時間尺度上進行分析,深入理解徑流與降水的集中度和集中期,并用集中期來反映梅江徑流對降水的滯后效應,旨在全面了解梅江流域的水文特征及其變化趨勢,為河流水量調節(jié)、防洪防汛及流域生態(tài)環(huán)境保護等方面提供決策依據(jù),也可以豐富華南濕熱區(qū)河流徑流和降水的研究案例。這也是首次用集中度和集中期的方法分析華南濕熱地區(qū)河流的徑流和降水的關系。

    1研究區(qū)概況

    梅江流域位于廣東省東部,梅州市的中西部,是韓江主流,也是梅州市最主要的河流,上游稱琴江,流經五華縣水寨與五華河匯合后稱梅江,由西南向東北流經五華、興寧、梅縣、梅江區(qū),至大埔縣的三河壩與汀江、梅潭河匯合后稱韓江。梅江的支流主要有五華河、北琴江、寧江、程江、石窟河、松源河等,干流全長307 km,流域集雨面積為13 929 km2,河床比降為0.59‰,地理坐標為東經115°08′—116°33′,北緯23°21′—25°10′。流域屬亞熱帶季風氣候區(qū),多年平均氣溫為21.1 ℃,無霜期一般都在320 d以上,多年平均降水量是1 489.4 mm,但降水量年內分配不均,主要集中在3—9月,約占年降水總量的82.5%,年均相對濕度80%左右。流域徑流完全由降水補給,徑流量、降水量兩者的年際變化規(guī)律比較接近。連續(xù)最大4個月徑流量基本上出現(xiàn)于5—8月,約占全年徑流量的53.1%,汛期(4—9月)徑流量約占全年徑流量的73.7%。流域在地貌上表現(xiàn)為低山、丘陵、盆地相見分布,流域屬南亞熱帶赤紅壤分布區(qū),分布著赤紅壤、紅壤、山黃壤和南方山地灌叢草甸土等?,F(xiàn)有植被主要是以馬尾松為主的次生林、桉樹林和稀疏草坡以及農作物,平均植被覆蓋度約為28%。

    2數(shù)據(jù)來源與研究方法

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

    本研究選取了興寧市水口水文站1960—2013年的觀測月徑流數(shù)據(jù)和梅縣氣象站1960—2013年的月降水數(shù)據(jù),徑流數(shù)據(jù)來自廣東省水利廳,降水數(shù)據(jù)來自廣東省氣象局。水口水文站設立于1954年,是梅江的重要控制站,該站以上有琴江、五華河和興寧江等匯合,控制流域面積6 480 km2,水情變化復雜,水雨情信息對梅州市防災減災起著無可替代的關鍵作用。該水文站位于水口鎮(zhèn)水口圩,東經115°53′,北緯23°59′,海拔109 m 。由于該站處于梅江干流中游,水文觀測資料時間序列較長且系統(tǒng)全面,位于梅江與寧江的匯入口附近,因此選擇該水文站實測徑流量數(shù)據(jù)作為分析數(shù)據(jù),具有一定的代表性,可以用來分析梅江干流水文特征。梅縣氣象站位于梅州市梅縣區(qū),是國家基本氣象站,設立于1952年,地理位置為東經116°04′,北緯24°15′,海拔89.3 m。研究區(qū)四季時間定義如下:春季(3—5月),夏季(6—8月),秋季(9—11月),冬季(12月至翌年2月)。

    2.2研究方法

    本文所采用的研究方法主要包括:累積距平法、變差系數(shù)法、集中度和集中期以及相關數(shù)理統(tǒng)計方法等。

    (1) 根據(jù)梅江流域年徑流量數(shù)據(jù),采用累積距平法得到年徑流量累積距平序列,從而繪制累積距平曲線。當累積距平持續(xù)增大時,表明該時段內徑流量距平持續(xù)為正;當累積距平持續(xù)不變時,表明該時段距平持續(xù)為零即保持平均;當累積距平持續(xù)減小時,表明該時段內徑流量距平持續(xù)為負。據(jù)此,可以較為直觀而準確地確定徑流量年際變化階段[12]。

    (2) 因為氣候的季節(jié)性變化,氣象要素中的降水和氣溫都有明顯的季節(jié)性波動,在一定程度上決定了徑流年內的季節(jié)性變化。本文用徑流年內變差系數(shù) 來衡量徑流年內季節(jié)性變化程度[13]。徑流年內變差系數(shù)Cv的計算公式如下:

    (1)

    (2)

    (3) 集中度是利用逐月徑流量(降水量)反映年內徑流量(降水量)集中程度的重要指標,集中期是利用逐月徑流量(降水量)反映年內最大徑流量(降水量)出現(xiàn)時段的重要指標[14]。集中度和集中期的計算是將一年內各月的徑流量作為向量看待,月徑流量的大小為向量的長度,所處的月份為向量的方向。1—12月每月的方位角h依次為0,30°,60°,…,330°,并把每個月的徑流量分解為x和y兩個方向上的分量,則x和y方向上的向量合成:

    (3)

    (4)

    (5)

    (6)

    (7)

    式中:Ri——第i月的徑流量(m3/s);h——第i月徑流的方位角(°);Rx,Ry——X,Y方向合成向量;R——Rx和Ry的合成總向量;Cd——集中度(m3/s);D——集中期。集中度和集中期的詳細計算過程見參考文獻[15]。

    3結果與分析

    3.1流域徑流特征分析

    根據(jù)水口水文站1960—2013年這54 a的徑流數(shù)據(jù),得出梅江多年平均徑流量為153 m3/s(4.83×109m3),1960—2013年,年徑流量呈波動變化(圖1),總體上徑流變化呈下降趨勢,但趨勢不明顯(r=0.161<0.222=R0.1,54,回歸顯著水平不到0.01)。由表1可見,流域最大徑流量出現(xiàn)在1983年,為316.1 m3/s(9.97×109m3),最小值出現(xiàn)在2009年,為39.5 m3/s(1.25×10 9m3),絕對變化幅度276.6 m3/s(8.72×109m3),相對變化幅度為8。根據(jù)徑流年內變差系數(shù) 的計算公式,得出梅江年徑流量的變差系數(shù)為0.366,說明其年際變化總體上較為平穩(wěn)。

    另外,從圖1中的累積距平曲線來看,54 a來梅江流域徑流呈現(xiàn)出明顯的階段性變化特征,1961—1972,2002—2012年徑流量以減少為主,累積距平曲線呈波動式下降狀態(tài),也表明最近10 a來河流進入枯水期;1973—2001年累積距平曲線表現(xiàn)為明顯的上升趨勢,累積距平為正,徑流量以增加為主,河流為豐水期。

    圖1 廣東省梅州市梅江徑流年際變化和累積距平曲線

    集水面積/km2多年平均流量/(m3·s-1)最大水年徑流量/(m3·s-1)年份最小水年徑流量/(m3·s-1)年份年際極值比變差系數(shù)Cv14061153316.1198339.5200980.366

    從梅江流域徑流的年內分配特征來看(圖2),徑流年內分配不均勻,主要集中在4—9月,季節(jié)上主要是在夏季和春季,占年徑流量的68%以上,各個年代的徑流也主要集中在這一階段;1—3和10—12月的徑流比較小。針對各個年代的最大徑流月份而言,1960,1970和2000s最大徑流出現(xiàn)在6月,1980和2010s出現(xiàn)在5月,1990s出現(xiàn)在8月,多年平均最大徑流出現(xiàn)在6月,這主要是各個年代徑流的集中時期不同造成的,這種現(xiàn)象的出現(xiàn)與降水量在年內時間上的分配有關,廣東緯度低,處于夏季風和臺風活動的前沿,河流屬暖季泛濫性河流,徑流全部靠降水補給。由于各年間季風和臺風強弱不同,來去的遲早和停留的時間長短不等,導致河流的徑流量存在年際上的變化,具有不穩(wěn)定性。

    圖2 廣東省梅州市梅江流域徑流年內分配特征

    梅江各季節(jié)徑流量均呈波動性變化(圖3),從各個季節(jié)徑流量的線性趨勢來看,春季徑流量呈現(xiàn)下降的趨勢(線性方程斜率為-5.31 mm/10 a),但趨勢并不明顯(r=0.084<0.092=R0.5,54,回歸顯著水平不到0.05);夏季徑流量也呈緩慢下降趨勢(線性方程斜率為-4.15 mm/10 a),但變化趨勢很不明顯(r=0.055<0.092=R0.5,54,回歸顯著水平不到0.05);秋季徑流量減少趨勢明顯(線性方程斜率為-12.97 mm/10 a,r=0.283>0.263=R0.05,54,回歸顯著水平在0.05以上);冬季徑流量同樣呈現(xiàn)緩慢下降趨勢(線性方程斜率為-0.69 mm/10 a),但變化趨勢極不明顯(r=0.032<0.092=R0.5,54,回歸顯著水平不到0.05)。春夏秋冬各季節(jié)平均徑流量為別為169.6,249.2,128.8,64.3 m3/s,徑流量所占年徑流量的百分比分別集中在23.6%~35.9%,34.7%~48.3%,17.2%~25.9%,8.7%~12.2%,均說明梅江的徑流主要集中在夏季和春季,秋冬季節(jié)徑流量較小。夏季和春季所占的徑流比例達68%以上,但其變化趨勢不明顯,而秋季徑流量下降趨勢明顯,整體上梅江流域年徑流量變化不明顯。春夏秋冬各季節(jié)的變差系數(shù)為:0.482,0.343,0.457,0.449,說明春季的年際徑流波動最大,夏季的年際徑流量較穩(wěn)定。

    由于徑流變化在月份和季節(jié)上存在有明顯差異,徑流分配特征隨年代變化也存在差異。為此,進一步討論分析梅江徑流量的集中程度和集中期,有利于加深對河流水文特征的了解,這對于華南濕熱區(qū)河流的水量調節(jié)、防洪防汛及流域生態(tài)保護等方面具有現(xiàn)實意義。

    梅江流域徑流集中度和集中期的年際變化表明(圖4),集中度圍繞平均值Cd=35%上下波動,最高達60.4%,最低為19.7%,說明梅江徑流主要集中在1 a內的某幾個月內,這與圖2體現(xiàn)的徑流主要集中在4—9月相符。梅江年內徑流集中期也存在上下波動趨勢,并且從趨勢線可以看出隨時間的推移發(fā)現(xiàn)集中期逐漸向后推移,從分布的時間范圍上來看,主要集中在5月下旬至7月中旬,54 a中有43 a落在這一時間段內,占總年份的80%。其中,集中在6月1日至6月30日的有27 a,集中在7月1—16日的有12 a,分別占總年份的的50%和22.2%。

    各年代的徑流統(tǒng)計特征也表現(xiàn)出一定的差異性,如表2所示,徑流集中度處于32.9%~42.4%,較為集中;集中期的合成向量處于141°~175.88°之間,轉化為徑流集中期的時間為5月23日至6月27日,多年平均徑流集中期的時間為6月17日。

    表2 流域徑流統(tǒng)計特征

    圖3 廣東省梅州市梅江流域季節(jié)徑流量及其占年徑流量比例變化

    圖4 廣東省梅州市梅江流域徑流集中度和集中期年際變化

    3.2流域降水特征分析

    如圖5所示,1960—2013年,梅江流域內降水量呈波動減少趨勢,但減少幅度很小,平均每10 a減少2.45 mm,多年平均降水量為1 489.4 mm。其中,最大降水量為2 355.4 mm,出現(xiàn)在1983年,最小降水量為1 011.3 mm,出現(xiàn)在1991年,絕對變化幅度為478.1 mm,相對變化幅度為2.33,變差系數(shù)為0.194。與流域內的徑流量相比,降水的相對變化幅度和變差系數(shù)均較小,說明流域內降水量的年際變化不大。

    從圖5的累積距平曲線可以看出,多年來梅江流域的降水量也呈現(xiàn)出相對明顯的階段性變化特征,1961—1971,2001—2004,2006—2012年降水量以減少為主,累積距平曲線呈明顯下降狀態(tài);1972—1994年累積距平曲線呈波動式上升狀態(tài),降水量以增加為主;1995—2000年降水變化不明顯。

    圖5 廣東省梅州市梅江流域降水年際變化和累積距平曲線

    梅江流域降水量的年內分配特征顯示(圖6),降水主要集中在4—6和8月,占到年降水量的55%以上,從季節(jié)上看主要是在春季和夏季,各年代的降水變化也主要體現(xiàn)在這一階段,1—3,7和9—12月所占比例較小。各年代最大降水出現(xiàn)的月份也有差異,具體而言,1970,1980和2010 s最大降水出現(xiàn)在5月,263.3,249.5,289.9 mm,1960和2000 s出現(xiàn)在6月份, 277.8和269.6 mm,而1990 s最大降水是出現(xiàn)在8月,為212.4 mm,這說明各年代或年內降水的分配特征和集中程度存在差異。同時,比較圖2與圖6可以看出,徑流峰值落后于降水,這說明徑流對降水有滯后效應。

    圖6 廣東省梅州市梅江流域降水年內分配特征

    各月的降水量數(shù)據(jù)得到梅江流域降水集中度和集中期時間序列(圖7),集中度圍繞平均值Cd=39.5%上下波動,最高達45.3%,最低為35.6%,但集中度的趨勢基本保持不變。較高的集中度表明,流域年內降水較為集中,這在年代尺度上的各月降水量變化也能體現(xiàn)(圖6),而降水集中期有隨時間向前推移的趨勢。同時,與梅江流域的徑流集中度相比,降水集中度大于徑流集中度,表明降水波動小于徑流,降水集中期更穩(wěn)定。從降水集中期的分布時間范圍來看,主要分布在5和6月,其中,集中在5月的有30 a,集中在6月的有14 a,分別占總年份的55.6%和25.9%。

    表3中的各年代降水年內分配特征顯示,集中程度較高(Cd≥35.6%),集中期的合成向量主要處于133.57°~154.84°之間,轉化為降水集中期的時間為5月15日至6月6日,多年平均降水量集中期的時間為5月28日。

    圖7 廣東省梅州市梅江流域降水集中度和集中期年際變化

    時間段降水量/mm集中度/%集中期合成向量方向/(°)最大降水出現(xiàn)時間1960s1442.643.1152.0206031970s1488.937.0154.8406061980s1578.139.5133.5705151990s1472.535.6140.9905232000s1461.045.3146.4905292010s1499.936.4150.880602多年平均1489.439.5145.860528

    3.3徑流對降水的響應

    梅江為典型的華南地區(qū)的外流河,河流徑流完全依賴降水補給,徑流量受降水量的直接影響,因此,需要考慮降水量和徑流量的關系。采用梅江流域年徑流量和年降水量進行關聯(lián)分析表明(圖8),兩者存在極大的相關性,相關系數(shù)為0.824?;貧w趨勢分析表明,年徑流量與年降水量的關系為線性回歸關系(r=0.824>0.428=R0.001,54),兩者的回歸方程有很好的擬合度,即徑流量會隨降水量的增加而增加,反之亦然。

    圖8 廣東省梅州市梅江流域年平均徑流量與降水量回歸分析

    隨著人類水土資源開發(fā)活動的日漸深入,傳統(tǒng)的自然水循環(huán)模式已逐漸被“自然—人工”二元水循環(huán)模式所取代。人類對水土資源的開發(fā)利用直接改變了影響產匯流過程的下墊面條件,改變了水的循環(huán)和轉化過程,使得流域水文情勢發(fā)生改變[16]。徑流形成過程是從降水落到流域表面至水流匯聚到流域出口斷面的整個復雜的物理過程,可以概括為產流過程和匯流過程。徑流是流域中氣候、下墊面條件和人類活動等眾多因素綜合作用的產物,由降水所形成的徑流需要經過停蓄階段、漫流階段、河槽集流階段3個階段,各個階段因降雨量大小、下墊面因素不同而需要耗費時間。因此,流域內每次降水匯流至水文站需要一定的時間,即降水所對應在水文站觀測的水文特征具有時間差,本研究嘗試應用降水和徑流的集中期來進一步探討兩者的內在聯(lián)系。

    從圖9可以看出,徑流集中期出現(xiàn)在降水集中期以后,且波動趨勢相近。徑流集中期相對降水集中期的滯后天數(shù)大多在5~40 d內,占到整個時間段的83.3%,平均滯后21 d,說明從流域降水開始,經過蓄流、匯流至觀測水文站需21 d左右的時間。此外,從圖9滯后天數(shù)的趨勢線來看,滯后天數(shù)有增大趨勢,1960—2000s間,各年代的平均滯后天數(shù)為20.8,23,16.3,28.6和24.9 d,隨時間推移呈上升趨勢。滯后天數(shù)的增加主要是因為徑流集中期向后推遲和降水集中期的向前推移。

    圖9 廣東省梅州市梅江流域降水和徑流的集中期對比及滯后時間年際變化

    水文現(xiàn)象具有必然性和隨機性,其演變規(guī)律不僅具有趨勢性、周期性等確定性演變特點,而且也具有不確定性的特征。隨著人類活動和全球氣候變化的影響,水文現(xiàn)象的不確定性更加突顯。流域的徑流對于降水具有滯后效應,基于集中度、集中期理論和梅江中游54 a的徑流量、降水量數(shù)據(jù),計算出梅江流域多年平均滯后天數(shù)長達21 d,且滯后天數(shù)隨時間推移有增加趨勢。這其中很重要的一個原因是人類活動改變地表覆被而使產流的各個階段時間延長,從而使徑流對降水的滯后時間延長[17-18],如將坡地改為梯田、退耕還林、植樹造林和蓄水工程建設等都能起到一定的蓄水作用,增加土壤的持水量和水的滯留時間[19]。梅江流域近30,40 a來,人口大量增加,工業(yè)化和城市化進程加快,這使得流域內天然林破壞嚴重,面積急劇減少,但同時因流域水土保持和經濟發(fā)展的需要,部分耕地轉變?yōu)榱值睾筒莸?,人工?主要包括用材林、經濟林和防護林等)的面積迅速增加,總體上流域內的林地和草地面積是呈不斷增長趨勢的,這對徑流產生了影響。另一方面,梅州市政府和社會對水利的投入增加,流域水利事業(yè)發(fā)展的步伐加快,水庫山塘、溝渠等蓄水和引水工程設施建設全面提速,這也很大程度上延長了降水停蓄時間。

    4討論與結論

    4.1討 論

    (1) 流域徑流對降水的多年滯后天數(shù)為21 d,這僅是應用集中度和集中期理論獲得的初步結論,而水文現(xiàn)象具有復雜性和不確定性,這就需要進一步通過獲取野外試驗數(shù)據(jù),明確降水強度對地表徑流的影響(降水入滲速率、徑流模數(shù)等參數(shù)都會影響徑流的形成時間),進而與此結論進行比較。

    (2) 本文主要分析了降水變化對徑流變化的影響,雖然降水是影響河川徑流最基本的因素,但實際上土地利用與土地覆被變化對流域徑流的減少也有顯著影響[20],徑流對降水的滯后時間隨年際變化呈現(xiàn)增大趨勢,這可能是因為地表覆被變化(如流域內植被指數(shù)NDVI的升高)或人類對降水進行攔截引起的,還有原因可能是氣候因素的變化。因此梅江流域內土地利用/覆被的變化對河川徑流到底存在什么樣的影響?影響程度如何?流域內徑流對降水滯后的時間增大是否與該流域土地利用/覆被的變化或者氣候的變化有關?這些都需要接下來進一步的探究。

    4.2結 論

    (1) 徑流量呈波動變化,主要集中在4—9月,占年徑流量的68%左右,集中期主要分布在5月下旬至7月中旬;流域內降水量趨勢幾十年來基本保持不變,主要集中在4—6和8月,占到年降水量的55%左右,集中期主要分布在5和6月。

    (2) 梅江各季節(jié)徑流量均呈波動性變化,從各個季節(jié)徑流量的線性趨勢來看,4個季節(jié)的徑流量均呈現(xiàn)下降的趨勢,其中秋季下降趨勢明顯,其他季節(jié)變化趨勢不明顯。

    (3) 從徑流量和降水量的累積距平曲線來看,梅江流域徑流量和降水量呈現(xiàn)出明顯的階段性變化特征,1961—1972和2001—2012年,兩者均呈現(xiàn)較為明顯的減少趨勢,1972—1994年,徑流量和降水量又表現(xiàn)出增加的趨勢。這也表明近10 a多來,梅江流域來水量整體呈偏枯態(tài)勢,氣候相對干燥。

    (4) 采用相關性分析法得出,流域內年徑流量與降水量存在極大的相關性,并呈線性回歸關系(r=0.824,大于R0.001,54=0.428),這表明了流域降水量增加會使徑流量增加,反之亦然。

    (5) 根據(jù)集中度和集中期理論,梅江流域的徑流對降水具有滯后效應,54 a平均時間為21 d,且滯后天數(shù)隨時間推移呈擴大趨勢。其中主要原因可能是人類活動,如將坡地改為梯田、退耕還林和蓄水工程建設等延長了降水停蓄時間。

    [參考文獻]

    [1]陳志愷.中國水資源的可持續(xù)利用問題[J].水文,2003,23(1):1-5.

    [2]劉昌明,成立.黃河干流下游斷流的徑流序列分析[J].地理學報,2000,55(3):257-266.

    [3]Ma Zhenmei, Kang Shaozhong, Zhang Lu, et al. Analysis of impacts of climate variability and human activity on streamflow for a river basin in arid region of northwest China[J]. Journal of Hydrology, 2008,352(3):239-249.

    [4]王兆禮,陳曉宏,楊濤.近50 a東江流域徑流變化及影響因素分析[J].自然資源學報,2010,25(8):1365-1374.

    [5]凌紅波,徐海量,張青青,等.葉爾羌河年徑流量與降水量的集中性及相關性分析[J].中國沙漠,2012,32(6):1757-1764.

    [6]Nash L L, Gleick P H. Sensitivity of streamflow in the Colorado basin to climatic changes[J]. Journal of Hydrology, 1991,125(3):221-241.

    [7]Milly P C D, Dunne K A, Vecchia A V. Global pattern of trends in streamflow and water availability in a changing climate[J]. Nature, 2005,438(7066):347-350.

    [8]Githui F, Gitau W, Mutua F, et al. Climate change impact on SWAT simulated streamflow in western Kenya[J]. International Journal of Climatology, 2009,29(12):1823-1834.

    [9]李善吉.對梅江河水中微量元素的區(qū)域性研究[J].廣東微量元素科學,2004,11(1):54-56.

    [10]黃燕戀,黃立民.梅江河近十年水質變化趨勢分析[J].廣東水利水電,2011,5(5):40-42.

    [11]李建東.韓江流域中上游:梅州地區(qū)水文特性分析[J].資源與環(huán)境,2012,17(2):34-36.

    [12]馬新萍,白紅英,侯欽磊,等.1959—2010年秦嶺灞河流域徑流量變化及其影響因素分析[J].資源科學,2012,34(7):1299-1303.

    [13]劉士余,章俊霞,羅志軍,等.近50 a贛西北大坑小流域徑流對降雨的響應[J].水土保持研究,2012,19(1):19-22.

    [14]劉賢趙,李嘉竹,宿慶,等.基于集中度與集中期的徑流年內分配研究[J].地理科學,2007,27(6):791-795.

    [15]楊遠東.河川徑流年內分配的計算方法[J].地理學報,1984,39(2):218-227.

    [16]王浩,嚴登華,楊大文,等.水文學方法研究[M].北京:科學出版社,2012:183-184.

    [17]劉昌明.黃河流域水循環(huán)演變若干問題的研究[J].水科學進展,2004,15(5):608-614.

    [18]王根緒,張鈺,劉桂民,等.馬營河流域1967—2000年土地利用變化對河流徑流的影響[J].中國科學:D輯,2005,35(7):671-681.

    [19]萬榮榮,楊桂山.流域土地利用/覆被變化的水文效應及洪水響應[J].湖泊科學,2004,16(3):259-264.

    [20]朱麗,秦富倉,姚云峰,等.北京市紅門川流域森林植被/土地覆被變化的水文響應[J].生態(tài)學報,2010,30(16):4287-4294.

    收稿日期:2015-05-11修回日期:2015-08-10

    通訊作者:張正棟(1968—)男(漢族),甘肅省榆中市人,博士,教授,主要從事土地利用、景觀生態(tài)和3S應用等方面的研究。E-mail:zhangzdedu@163.com。

    文獻標識碼:A

    文章編號:1000-288X(2016)03-0038-08

    中圖分類號:P333

    Runoff Characteristics and Hysteresis to Precipitation in Meijiang River Basin in Guangdong Province During 1960-2013

    DONG Caiwen, ZHANG Zhengdong, YANG Chuanxun, WAN Luwen, ZHANG Qian

    (SchoolofGeography,SouthChinaNormalUniversity,Guangdong,Guangzhou510631,China)

    Abstract:[Objective] Taking Meijiang river basin as a case study, this paper investigates the variations of precipitation and runoff in order to provide supports for river regulation, flood control and environmental protection. [Methods] Runoff and precipitation data during 1960—2013 was collected. Mathematical methods including accumulative anomaly, variation coefficient, concentration frequency and concentration period were used to investigate the variations of runoff and precipitation at decadal, annual and seasonal time scales. The hysteresis effect of runoff to precipitation was discussed. [Results] The inter-annual changes of runoff was declined at the rate of 5.79 mm/10 a. The runoff was mainly concentrated in the summer and spring time from April to September. The fluctuation of precipitation showed a slight decreasing trend in the basin. The precipitation mainly concentrated in May to June. Correlation between annual runoff and precipitation was significant(the significant level is more than 0.001, r=0.824>0.428=R0.001,54). There was a hysteresis effect of runoff to precipitation, the average hysteresis period was about 21 days, and this hysteresis period was expected to increase over time. [Conclusion] Human activity may be the main reason that alters the status of land cover or the change of climate.

    Keywords:runoff; precipitation; hysteresis effect; concentration frequency; concentration period; Meijiang River basin

    資助項目:國家自然科學基金項目“變化環(huán)境下南方濕熱區(qū)韓江流域景觀格局演變與生態(tài)水文耦合研究”(41471147)

    第一作者:董才文(1989—),男(漢族),湖北省黃石市人,碩士研究生,研究方向為景觀生態(tài)與土地利用。E-mail:540329883@qq.com。

    猜你喜歡
    滯后效應集中度徑流
    多個大股東治理對企業(yè)投資的影響研究
    會計之友(2025年3期)2025-01-22 00:00:00
    京津冀縣域人口集中度分析
    客聯(lián)(2022年10期)2022-07-06 09:06:16
    新廣告商:廣告業(yè)周期性在弱化,而集中度在提升 精讀
    英語文摘(2020年11期)2020-02-06 08:53:24
    基于環(huán)境保護的企業(yè)社會責任與企業(yè)財務績效的關系研究
    商情(2016年32期)2017-03-04 00:54:25
    保險公司資本結構、業(yè)務集中度與再保險需求研究
    煤炭行業(yè)未來在提高集中度
    能源(2016年3期)2016-12-01 05:10:51
    城鎮(zhèn)化中人口結構變化與經濟增長的關系
    Topmodel在布哈河流域徑流模擬中的應用
    探秘“大徑流”
    攻克“大徑流”
    亚洲av.av天堂| 搡老乐熟女国产| av.在线天堂| 狂野欧美激情性bbbbbb| 中文字幕人妻丝袜制服| 免费不卡的大黄色大毛片视频在线观看| 久久国产精品男人的天堂亚洲 | av有码第一页| 午夜激情av网站| 成人国语在线视频| 久久久久久久久久久免费av| 新久久久久国产一级毛片| 国产成人a∨麻豆精品| 亚洲天堂av无毛| av黄色大香蕉| 99九九在线精品视频| 亚洲内射少妇av| 99国产精品免费福利视频| 亚洲精品国产av成人精品| 日本午夜av视频| 成人综合一区亚洲| 一级毛片我不卡| 久久久久国产精品人妻一区二区| 欧美日韩成人在线一区二区| 国产深夜福利视频在线观看| 美女xxoo啪啪120秒动态图| 免费在线观看完整版高清| 黄色配什么色好看| 国产黄色免费在线视频| 婷婷色av中文字幕| av免费在线看不卡| 欧美激情国产日韩精品一区| 久久久久网色| freevideosex欧美| 考比视频在线观看| 亚洲五月色婷婷综合| 久久久久久久精品精品| 制服人妻中文乱码| 九色成人免费人妻av| 精品国产乱码久久久久久小说| 汤姆久久久久久久影院中文字幕| 亚洲一码二码三码区别大吗| 国产男女超爽视频在线观看| 免费看光身美女| 久久久久久久久久成人| 97超碰精品成人国产| 在线观看免费视频网站a站| 久久久a久久爽久久v久久| 亚洲国产欧美日韩在线播放| 麻豆乱淫一区二区| 午夜激情av网站| 最后的刺客免费高清国语| 啦啦啦视频在线资源免费观看| 久久人人爽av亚洲精品天堂| 制服人妻中文乱码| 97在线视频观看| 精品久久蜜臀av无| 欧美国产精品一级二级三级| 伦理电影免费视频| 国产欧美另类精品又又久久亚洲欧美| 日本色播在线视频| 国产精品偷伦视频观看了| 天天躁夜夜躁狠狠久久av| 男女免费视频国产| 免费观看av网站的网址| 国产精品久久久久久久电影| 国产成人精品在线电影| 人人妻人人澡人人爽人人夜夜| 不卡视频在线观看欧美| 久久久久精品人妻al黑| 日本黄色日本黄色录像| 亚洲,欧美,日韩| 亚洲,欧美,日韩| 熟女av电影| 赤兔流量卡办理| 久久久久精品性色| 亚洲综合精品二区| 黑人高潮一二区| 最近最新中文字幕免费大全7| 我的女老师完整版在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩三级伦理在线观看| 黑人高潮一二区| 免费av中文字幕在线| 日本爱情动作片www.在线观看| 精品国产国语对白av| 热99国产精品久久久久久7| 最近的中文字幕免费完整| 建设人人有责人人尽责人人享有的| 老熟女久久久| 秋霞伦理黄片| 欧美激情国产日韩精品一区| 少妇高潮的动态图| 大话2 男鬼变身卡| 国产极品天堂在线| 国产国语露脸激情在线看| 大片电影免费在线观看免费| 啦啦啦中文免费视频观看日本| 啦啦啦视频在线资源免费观看| 欧美 亚洲 国产 日韩一| 在线免费观看不下载黄p国产| 亚洲av综合色区一区| 高清在线视频一区二区三区| 精品一品国产午夜福利视频| 久久久久久久久久成人| 90打野战视频偷拍视频| 中文天堂在线官网| 国产精品一国产av| 菩萨蛮人人尽说江南好唐韦庄| xxx大片免费视频| 久久影院123| 99久久综合免费| 女性被躁到高潮视频| 51国产日韩欧美| 国产精品久久久久久久久免| 久久精品国产亚洲av天美| 亚洲国产看品久久| 久久久久精品人妻al黑| 一级毛片黄色毛片免费观看视频| 久久久亚洲精品成人影院| 国产国语露脸激情在线看| 有码 亚洲区| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区| 美女视频免费永久观看网站| 亚洲av综合色区一区| 美女视频免费永久观看网站| 成人午夜精彩视频在线观看| 自线自在国产av| 精品视频人人做人人爽| 人人妻人人爽人人添夜夜欢视频| 亚洲人成77777在线视频| 午夜日本视频在线| 日韩制服骚丝袜av| 亚洲国产精品一区三区| 一级毛片我不卡| 免费看av在线观看网站| 日韩精品免费视频一区二区三区 | 久久久久精品人妻al黑| 高清黄色对白视频在线免费看| 少妇人妻 视频| 中文字幕最新亚洲高清| 亚洲精品国产av成人精品| 老熟女久久久| 久久久久国产网址| 国产精品 国内视频| 久久精品国产自在天天线| 九色成人免费人妻av| 国产免费一区二区三区四区乱码| 啦啦啦啦在线视频资源| 热99久久久久精品小说推荐| 制服丝袜香蕉在线| 亚洲av男天堂| 亚洲成人手机| 精品99又大又爽又粗少妇毛片| 成年人免费黄色播放视频| freevideosex欧美| 各种免费的搞黄视频| 色视频在线一区二区三区| 51国产日韩欧美| 午夜精品国产一区二区电影| 日本欧美视频一区| 十八禁网站网址无遮挡| 国产黄色免费在线视频| 亚洲国产精品专区欧美| 黄片播放在线免费| 麻豆精品久久久久久蜜桃| 伊人久久国产一区二区| 成年人免费黄色播放视频| 久久人妻熟女aⅴ| 波多野结衣一区麻豆| 一级片'在线观看视频| 中国国产av一级| 国产精品熟女久久久久浪| 在线亚洲精品国产二区图片欧美| 国产精品国产av在线观看| 免费观看av网站的网址| 搡女人真爽免费视频火全软件| 精品一区在线观看国产| 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久| 香蕉丝袜av| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 美女中出高潮动态图| 亚洲,欧美,日韩| 亚洲成人av在线免费| 免费观看av网站的网址| 日本与韩国留学比较| 精品久久国产蜜桃| 美女内射精品一级片tv| 婷婷色av中文字幕| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| √禁漫天堂资源中文www| 大香蕉97超碰在线| av黄色大香蕉| 黄色一级大片看看| 午夜av观看不卡| 精品一区在线观看国产| 国产片特级美女逼逼视频| 韩国精品一区二区三区 | 精品一区在线观看国产| 国产精品久久久久久精品电影小说| 国产精品麻豆人妻色哟哟久久| 大陆偷拍与自拍| 亚洲av成人精品一二三区| 天美传媒精品一区二区| 97精品久久久久久久久久精品| 一级爰片在线观看| 大陆偷拍与自拍| av片东京热男人的天堂| 久久青草综合色| 国产成人欧美| 午夜福利,免费看| 18在线观看网站| 亚洲伊人色综图| 国产成人欧美| 免费在线观看完整版高清| 男女下面插进去视频免费观看 | 国产无遮挡羞羞视频在线观看| 国产av一区二区精品久久| 永久免费av网站大全| 一本色道久久久久久精品综合| 精品久久久精品久久久| 亚洲美女搞黄在线观看| 亚洲精品美女久久av网站| 只有这里有精品99| 亚洲欧洲精品一区二区精品久久久 | 精品午夜福利在线看| av国产久精品久网站免费入址| 女人被躁到高潮嗷嗷叫费观| 蜜桃在线观看..| 欧美亚洲 丝袜 人妻 在线| 国产精品免费大片| 热99国产精品久久久久久7| 免费黄频网站在线观看国产| 国产一区二区激情短视频 | 新久久久久国产一级毛片| 五月开心婷婷网| 国产 一区精品| 日本av手机在线免费观看| 亚洲欧美日韩另类电影网站| 在线观看人妻少妇| av.在线天堂| 美女国产视频在线观看| 18禁在线无遮挡免费观看视频| 久久人人97超碰香蕉20202| 国产在线视频一区二区| 日韩不卡一区二区三区视频在线| 国产av国产精品国产| 免费黄网站久久成人精品| 国产高清不卡午夜福利| 亚洲av男天堂| 曰老女人黄片| videos熟女内射| 90打野战视频偷拍视频| av福利片在线| 天堂8中文在线网| 少妇熟女欧美另类| 精品视频人人做人人爽| 伊人亚洲综合成人网| 日韩av免费高清视频| 2021少妇久久久久久久久久久| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 九色成人免费人妻av| a 毛片基地| 国产精品国产av在线观看| 国产精品 国内视频| √禁漫天堂资源中文www| 青春草国产在线视频| 欧美激情极品国产一区二区三区 | 国产成人aa在线观看| 18+在线观看网站| 亚洲一区二区三区欧美精品| 久久人人爽人人爽人人片va| 亚洲国产色片| 成年女人在线观看亚洲视频| 亚洲精品视频女| 欧美国产精品一级二级三级| 国产精品久久久久久精品电影小说| 一级爰片在线观看| 国产熟女欧美一区二区| 日本wwww免费看| 卡戴珊不雅视频在线播放| 精品亚洲乱码少妇综合久久| 亚洲欧美精品自产自拍| 又粗又硬又长又爽又黄的视频| 三级国产精品片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av在线观看美女高潮| 日韩熟女老妇一区二区性免费视频| 亚洲av男天堂| 中文字幕精品免费在线观看视频 | 亚洲国产av影院在线观看| 亚洲国产最新在线播放| 涩涩av久久男人的天堂| 欧美日韩一区二区视频在线观看视频在线| 久久久国产一区二区| 国产精品一区www在线观看| 亚洲国产看品久久| 丝袜人妻中文字幕| 日韩成人av中文字幕在线观看| 亚洲精品456在线播放app| 免费大片黄手机在线观看| 日日撸夜夜添| 在现免费观看毛片| 少妇的逼水好多| av网站免费在线观看视频| 欧美精品av麻豆av| 美女视频免费永久观看网站| 亚洲精品美女久久久久99蜜臀 | 一级a做视频免费观看| 国产精品一二三区在线看| 亚洲国产日韩一区二区| 精品国产乱码久久久久久小说| 又黄又粗又硬又大视频| 国产色婷婷99| 日韩av不卡免费在线播放| 亚洲精品久久午夜乱码| av黄色大香蕉| 母亲3免费完整高清在线观看 | 制服诱惑二区| 精品少妇内射三级| 国产xxxxx性猛交| 香蕉国产在线看| 大片免费播放器 马上看| 丝袜喷水一区| 国产精品国产三级专区第一集| 卡戴珊不雅视频在线播放| 美女中出高潮动态图| 亚洲天堂av无毛| 咕卡用的链子| 亚洲国产精品999| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 两个人免费观看高清视频| 日本午夜av视频| 一区二区三区乱码不卡18| 国产精品国产av在线观看| 久久精品久久久久久噜噜老黄| 精品久久久精品久久久| 久久精品aⅴ一区二区三区四区 | 国产在线视频一区二区| 国产欧美另类精品又又久久亚洲欧美| 亚洲五月色婷婷综合| 亚洲成av片中文字幕在线观看 | 永久网站在线| 国产激情久久老熟女| 建设人人有责人人尽责人人享有的| 久久久久国产网址| 日本vs欧美在线观看视频| 水蜜桃什么品种好| 精品少妇黑人巨大在线播放| 亚洲综合色惰| 亚洲人与动物交配视频| 亚洲精品美女久久av网站| 又粗又硬又长又爽又黄的视频| 久久精品久久久久久久性| 九九在线视频观看精品| 丰满乱子伦码专区| 久久人人97超碰香蕉20202| 精品人妻偷拍中文字幕| 久久97久久精品| 男女无遮挡免费网站观看| 制服诱惑二区| videossex国产| 国产成人午夜福利电影在线观看| 草草在线视频免费看| 水蜜桃什么品种好| av有码第一页| 日韩中文字幕视频在线看片| 国产精品久久久久久精品古装| 香蕉丝袜av| 国产精品一区二区在线不卡| 老熟女久久久| 男女边摸边吃奶| 日日摸夜夜添夜夜爱| 少妇精品久久久久久久| 汤姆久久久久久久影院中文字幕| 精品卡一卡二卡四卡免费| 午夜福利,免费看| 爱豆传媒免费全集在线观看| 热re99久久精品国产66热6| 九九在线视频观看精品| xxxhd国产人妻xxx| 日韩电影二区| 两性夫妻黄色片 | av在线老鸭窝| 欧美精品一区二区免费开放| 高清欧美精品videossex| 亚洲人成网站在线观看播放| 人人妻人人澡人人看| 曰老女人黄片| 最近中文字幕2019免费版| 人人妻人人添人人爽欧美一区卜| 午夜福利视频在线观看免费| 国产 一区精品| 中文乱码字字幕精品一区二区三区| av有码第一页| 男女高潮啪啪啪动态图| 久久精品国产a三级三级三级| 男女边摸边吃奶| 婷婷色综合www| 亚洲国产日韩一区二区| √禁漫天堂资源中文www| 男人舔女人的私密视频| 成人毛片a级毛片在线播放| 高清不卡的av网站| 免费看不卡的av| 国产精品国产三级国产av玫瑰| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 欧美日韩国产mv在线观看视频| 熟女av电影| 中文精品一卡2卡3卡4更新| 午夜精品国产一区二区电影| 亚洲精品自拍成人| 亚洲av免费高清在线观看| 女人久久www免费人成看片| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品美女久久久久99蜜臀 | www日本在线高清视频| 亚洲精品美女久久久久99蜜臀 | 人妻少妇偷人精品九色| 国产成人精品福利久久| 亚洲欧洲精品一区二区精品久久久 | 免费av不卡在线播放| 国产片内射在线| 男女午夜视频在线观看 | 亚洲成人一二三区av| 日韩一本色道免费dvd| 中文字幕精品免费在线观看视频 | 亚洲色图综合在线观看| 国产淫语在线视频| 一本—道久久a久久精品蜜桃钙片| 国产欧美日韩一区二区三区在线| 18在线观看网站| 日本欧美国产在线视频| 五月开心婷婷网| 人人澡人人妻人| 中国美白少妇内射xxxbb| 国产精品人妻久久久影院| 自线自在国产av| 香蕉国产在线看| 亚洲精品第二区| 制服丝袜香蕉在线| 国产成人a∨麻豆精品| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 美女中出高潮动态图| av免费在线看不卡| 中文字幕制服av| 一区二区三区精品91| 中文乱码字字幕精品一区二区三区| 欧美最新免费一区二区三区| 国精品久久久久久国模美| 极品人妻少妇av视频| 全区人妻精品视频| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 国产一区二区在线观看av| 蜜桃国产av成人99| 丝袜人妻中文字幕| 欧美日韩av久久| 狠狠婷婷综合久久久久久88av| 啦啦啦视频在线资源免费观看| 在线精品无人区一区二区三| 一级毛片电影观看| 免费少妇av软件| 久久久久久伊人网av| 久久久久久久国产电影| 国产黄频视频在线观看| 狂野欧美激情性bbbbbb| 寂寞人妻少妇视频99o| 亚洲国产最新在线播放| 大香蕉久久网| 中国国产av一级| 日韩电影二区| 一本—道久久a久久精品蜜桃钙片| 熟妇人妻不卡中文字幕| 69精品国产乱码久久久| 啦啦啦啦在线视频资源| 汤姆久久久久久久影院中文字幕| 国产精品嫩草影院av在线观看| 久久久久久久久久久免费av| 人人妻人人添人人爽欧美一区卜| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 午夜福利乱码中文字幕| 国产亚洲欧美精品永久| videos熟女内射| 亚洲综合色惰| 久久人人97超碰香蕉20202| 五月伊人婷婷丁香| 晚上一个人看的免费电影| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 一级,二级,三级黄色视频| 97人妻天天添夜夜摸| 日韩av在线免费看完整版不卡| 九色亚洲精品在线播放| 永久免费av网站大全| 国内精品宾馆在线| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 高清视频免费观看一区二区| 国产亚洲最大av| 男人添女人高潮全过程视频| 五月玫瑰六月丁香| 国产亚洲av片在线观看秒播厂| 菩萨蛮人人尽说江南好唐韦庄| 亚洲图色成人| 国产一区二区三区av在线| 欧美3d第一页| 亚洲美女黄色视频免费看| 欧美3d第一页| 国产精品一国产av| 国产无遮挡羞羞视频在线观看| 国产精品一国产av| 欧美亚洲日本最大视频资源| 国产精品无大码| 欧美成人精品欧美一级黄| 国产精品无大码| 丰满饥渴人妻一区二区三| 国产精品国产三级专区第一集| 日韩一区二区三区影片| 22中文网久久字幕| 男人爽女人下面视频在线观看| 免费女性裸体啪啪无遮挡网站| 最近手机中文字幕大全| 国产熟女欧美一区二区| 春色校园在线视频观看| av.在线天堂| 最近手机中文字幕大全| 久久国产精品大桥未久av| 国产免费一区二区三区四区乱码| 久久久国产精品麻豆| av卡一久久| 国产一区二区激情短视频 | 亚洲美女黄色视频免费看| 成人二区视频| 婷婷色综合大香蕉| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 丁香六月天网| 99久久中文字幕三级久久日本| 亚洲伊人色综图| 中国三级夫妇交换| 91成人精品电影| 亚洲av福利一区| 在线观看一区二区三区激情| 又粗又硬又长又爽又黄的视频| 久久精品熟女亚洲av麻豆精品| 有码 亚洲区| 午夜视频国产福利| 国产精品一区二区在线不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品人妻偷拍中文字幕| 麻豆乱淫一区二区| 这个男人来自地球电影免费观看 | 久久久久视频综合| 国产日韩欧美在线精品| 少妇的逼好多水| 青春草视频在线免费观看| 男女国产视频网站| 午夜福利,免费看| xxxhd国产人妻xxx| 夫妻性生交免费视频一级片| 一本久久精品| 寂寞人妻少妇视频99o| 美女国产高潮福利片在线看| 国产永久视频网站| 黄网站色视频无遮挡免费观看| 男人操女人黄网站| 飞空精品影院首页| 菩萨蛮人人尽说江南好唐韦庄| 国产一区有黄有色的免费视频| 国产亚洲欧美精品永久| 国产av一区二区精品久久| 黑人欧美特级aaaaaa片| 少妇高潮的动态图| 高清欧美精品videossex| av网站免费在线观看视频| a级毛色黄片| 欧美日韩精品成人综合77777| 乱人伦中国视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产毛片av蜜桃av| 男女高潮啪啪啪动态图| 国产毛片在线视频| 91精品国产国语对白视频| 肉色欧美久久久久久久蜜桃| 亚洲国产精品成人久久小说| 这个男人来自地球电影免费观看 | 各种免费的搞黄视频| 精品人妻在线不人妻| 免费大片18禁| 黄色一级大片看看| 中文字幕精品免费在线观看视频 | 侵犯人妻中文字幕一二三四区| 男人操女人黄网站| 高清在线视频一区二区三区| 国产精品免费大片| 亚洲在久久综合| 免费观看无遮挡的男女| 国产熟女欧美一区二区| 最近2019中文字幕mv第一页| 久久99精品国语久久久| 久久人人97超碰香蕉20202| 黑人高潮一二区| 尾随美女入室| 王馨瑶露胸无遮挡在线观看|