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

    南北過(guò)渡帶典型流域徑流變化歸因?qū)Ρ妊芯?/h1>
    2020-12-30 12:55:30師潤(rùn)田鵬趙廣舉
    人民黃河 2020年12期
    關(guān)鍵詞:過(guò)渡帶下墊面源區(qū)

    師潤(rùn) 田鵬 趙廣舉

    摘 要:為探明中國(guó)南北過(guò)渡帶不同河川徑流變化的差異,采用線性趨勢(shì)回歸法與累積距平法對(duì)比分析了1975—2015年南北過(guò)渡帶秦嶺南麓漢江源武侯站與北麓渭河源武山站的徑流變化特征,采用基于Budyko水熱平衡方程的彈性系數(shù)法開展徑流變化的歸因識(shí)別。研究發(fā)現(xiàn):漢江源區(qū)徑流呈微弱減少趨勢(shì),渭河源區(qū)徑流呈極顯著減少趨勢(shì),分別于1990年和1993年發(fā)生突變;兩流域年徑流量均對(duì)氣候變化較為敏感,特別是漢江源區(qū),其徑流對(duì)降水變化表現(xiàn)出較高的敏感性;降水、潛在蒸散發(fā)和下墊面變化對(duì)漢江源徑流變化的貢獻(xiàn)率分別為65.23%、6.02%和28.75%,對(duì)渭河源徑流變化的貢獻(xiàn)率分別為24.51%、26.62%和48.87%。研究結(jié)果表明渭河源區(qū)徑流對(duì)人類活動(dòng)影響更加敏感。

    關(guān)鍵詞:南北過(guò)渡帶;徑流變化;Budyko水熱平衡方程;氣候變化;人類活動(dòng);漢江源區(qū);渭河源區(qū)

    中圖分類號(hào):P333;TV121 文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.12.006

    Abstract:This study attempted to identify the differences in runoff changes of two rivers in the south and north of Chinas north-south transitional zone. The linear regression test and the cumulative anomaly method were applied to detect temporal variations in runoff at Wuhou and Wushan stations in the sources of Hanjiang River and Weihe River from 1975 to 2015. The elasticity coefficient method based on the Budyko hydrothermal balance equation was applied to quantify the attributions of climate change and land surface changes. The results show that the annual runoff at Wuhou in Hanjiang River shows a slightly decreasing trend, whereas runoff at Wushan Station shows a significant decreasing trend. Abrupt changes are detected in 1990 and 1993 respectively. Runoff in both watersheds shows high sensitivity to climate change, especially in the source of the Hanjiang River, its runoff shows high sensitivity to the changes of precipitation. The changes in precipitation, potential evapotranspiration and land surface contribute 65.23%, 6.02% and 28.75% respectively to the runoff changes in the source of the Hanjiang River and 24.51%, 26.62% and 48.87% respectively in the source of the Weihe River. Furthermore, the runoff in the source of the Weihe River is highly influenced by human activities comparing to that in the source of Hanjiang River.

    Key words: northsouth transition zone; runoff changes; Budyko equation; climate change; human activity; the source of the Hanjiang River; the source of the Weihe River

    以秦嶺—大巴山為主體的中國(guó)南北過(guò)渡帶,地跨溫帶和亞熱帶,是中國(guó)南北氣候和地理的天然界線[1]。南北過(guò)渡帶是許多河流的源頭,包括長(zhǎng)江流域的漢江、嘉陵江、岷江、沱江及黃河流域的渭河、伊洛河等[2]。該地區(qū)水資源豐富,其中漢江是陜西省引漢濟(jì)渭和國(guó)家南水北調(diào)中線重要引水工程的優(yōu)質(zhì)水源地[3-4],跨地區(qū)調(diào)水保障了沿線和受水地區(qū)社會(huì)生產(chǎn)和人民生活基本用水。受氣候變化和人類活動(dòng)的影響,南北過(guò)渡帶地區(qū)的河川徑流量日益減少,導(dǎo)致水資源時(shí)空分布不均,研究南北過(guò)渡帶河川徑流變化及其驅(qū)動(dòng)因素對(duì)緩解區(qū)域水資源供需矛盾、促進(jìn)區(qū)域社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展具有重要意義。

    氣候波動(dòng)和人類活動(dòng)是導(dǎo)致河川徑流變化的主要因素[5-6],有關(guān)學(xué)者在南北過(guò)渡帶地區(qū)開展了相關(guān)的研究,如柯新月等[7]采用徑流還原法分析灞河和洵河流域1956—2011年的實(shí)測(cè)徑流量,發(fā)現(xiàn)人類活動(dòng)是導(dǎo)致兩流域徑流減少的主要原因,貢獻(xiàn)率超過(guò)了75%;李凌程等[8]定量評(píng)估了南水北調(diào)中線工程漢江上游流域徑流變化的主要驅(qū)動(dòng)因素,發(fā)現(xiàn)該流域1991—1999年徑流減少主要是受氣候因素驅(qū)動(dòng),而2000—2008年人類活動(dòng)影響逐步增強(qiáng)且占主導(dǎo)地位。大量研究基于經(jīng)驗(yàn)統(tǒng)計(jì)法[9-12]和水文模型法[13-15],分析了不同流域徑流變化特征與驅(qū)動(dòng)因素,但針對(duì)南北過(guò)渡帶不同流域河川徑流變化的對(duì)比研究相對(duì)較少。筆者選擇南北過(guò)渡帶的漢江源頭武侯水文站以上流域和渭河源頭武山水文站以上流域?yàn)檠芯繀^(qū),采用多種統(tǒng)計(jì)方法辨析實(shí)測(cè)徑流量變化趨勢(shì)及突變特征,基于Budyko水熱平衡方程識(shí)別氣候因子對(duì)徑流變化的貢獻(xiàn)率,辨識(shí)南北過(guò)渡帶典型流域徑流變化主要原因,以期為區(qū)域水資源可持續(xù)利用提供科學(xué)依據(jù)。

    1 研究區(qū)概況與數(shù)據(jù)獲取

    1.1 研究區(qū)概況

    漢江源流域位于漢江上游干流武侯水文站以上(東經(jīng)106°00′—106°50′,北緯32°30′—33°40′),流域呈帶狀,面積為3 092 km2,多年(1975—2015年)平均徑流量為10.45億m3。漢江發(fā)源于陜西省寧強(qiáng)縣秦嶺南麓,漢江源流域水系(見圖1(b))主要由干流玉帶河和支流沮水河組成。玉帶河全長(zhǎng)101.1 km,河道平均比降為8.2%,流域面積為831 km2,占流域總面積的26.88%;沮水河長(zhǎng)130 km,第二支流較多,流域面積為1 747 km2,占流域總面積的56.50%。整個(gè)流域?qū)俦眮啛釒瘻貪駶?rùn)季風(fēng)氣候區(qū),降水均勻且豐沛,多年平均降水量為853.16 mm,多年平均氣溫為12.7 ℃,多年平均日照時(shí)數(shù)為1 585.5 h。

    渭河源流域位于黃河中游渭河干流武山水文站以上(東經(jīng)104°00′—105°00′ ,北緯34°20′—35°30′),流域呈葉片狀,總面積為8 080 km2,多年(1975—2015年)平均徑流量為4.80億m3。渭河發(fā)源于甘肅省渭源縣鳥鼠山,渭河源流域水系(見圖1(c))主要由左岸秦祁河、咸河和右岸榜沙河3大支流組成。秦祁河和咸河全長(zhǎng)分別為67.5、69.0 km,平均比降分別為6.0%、5.3%,流域面積分別為858、1 159 km2,左岸水少質(zhì)差,洪水含沙量大。榜沙河流經(jīng)秦嶺北麓山地,河道總長(zhǎng)102.6 km,支流頗多,流域面積為3 597 km2,年徑流量為4.39億m3,右岸徑流來(lái)源以降水為主,植被較好,水量大,含沙量小。整個(gè)流域地處秦嶺北麓西段黃土高原邊緣復(fù)合地帶,受大陸性季風(fēng)氣候影響,光熱資源豐富,降水極不均勻,多年平均降水量為468.52 mm,多年平均氣溫為5.8 ℃,多年平均日照時(shí)數(shù)為2 348.4 h。

    1.2 數(shù)據(jù)來(lái)源與處理

    根據(jù)漢江源區(qū)和渭河源區(qū)水文站的分布特點(diǎn),并考慮數(shù)據(jù)的準(zhǔn)確性、完整性和連續(xù)性,選取武侯水文站和武山水文站及其周邊8個(gè)具有代表性的氣象站(略陽(yáng)、勉縣、留壩、寧強(qiáng)和渭源、漳縣、隴西、武山),如圖1所示。武侯站設(shè)立在陜西省漢中市勉縣武侯鎮(zhèn),武山站于1974年7月1日由原車家川站下遷至甘肅省天水市武山縣。收集整理兩個(gè)水文站1975—2015年徑流數(shù)據(jù),研究南北過(guò)渡帶局部地區(qū)1975—2015年徑流演變特征,并選取各氣象站同時(shí)間序列的逐日氣侯數(shù)據(jù)(包括降水、最高和最低以及平均氣溫、日照時(shí)數(shù)、平均風(fēng)速、相對(duì)濕度),采用算術(shù)平均法計(jì)算流域平均降水量,采用被廣泛應(yīng)用的FAO修正PenmanMonteith公式[16]計(jì)算流域年潛在蒸散發(fā)量。水文數(shù)據(jù)來(lái)源于《長(zhǎng)江水文年鑒》《黃河水文年鑒》和《黃河泥沙公報(bào)》,氣象數(shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)共享服務(wù)網(wǎng)(http://data.cma.cn)。

    2 方 法

    2.1 趨勢(shì)及突變分析

    采用線性趨勢(shì)回歸法[17]和累積距平法分析年徑流、降水和潛在蒸散發(fā)序列的變化趨勢(shì)。線性回歸法較為簡(jiǎn)單,在此不作贅述。累積距平法可由曲線直觀反映序列的變化趨勢(shì),結(jié)合實(shí)測(cè)值,檢驗(yàn)序列的演變趨勢(shì)及突變點(diǎn)[18-19]。

    2.2 徑流變化影響因子敏感性分析

    基于Budyko水熱平衡方程探究徑流變化的主要原因已成為研究熱點(diǎn),其中,彈性系數(shù)法最具代表性,與數(shù)理統(tǒng)計(jì)方法相比,其物理意義更明顯,所用數(shù)據(jù)更易獲取[20-21]。

    2.3 降水量、潛在蒸散發(fā)量、下墊面變化對(duì)徑流變化的貢獻(xiàn)

    根據(jù)突變點(diǎn)對(duì)應(yīng)年份將研究序列劃分為兩個(gè)時(shí)段,徑流變化(dR)為后一時(shí)段的多年平均徑流深R-2減去前一時(shí)段的多年平均徑流深R-1。

    3 結(jié)果與分析

    3.1 南北過(guò)渡帶典型流域水文情勢(shì)變化分析

    采用線性趨勢(shì)回歸法分析武侯站和武山站1975—2015年年徑流量變化趨勢(shì),結(jié)果如圖2所示。武侯站年徑流量呈微弱減少趨勢(shì),年均減少0.070億m3;而武山站年徑流量減少趨勢(shì)極顯著(α<0.01)(見表1),年均減少0.128億m3。武侯站的年徑流量出現(xiàn)多個(gè)峰值,振蕩幅度較大,較為離散,極值差異顯著,變異系數(shù)CV=0.57,多年平均徑流量為10.45億m3。武山站的年徑流量減少平緩,數(shù)據(jù)離散程度不大,無(wú)明顯峰值,極值差異不大,變異系數(shù)CV=0.49,多年平均徑流量為4.80億m3,約為武侯站的1/2。總體上,兩站1975—1992年年徑流量的振蕩幅度均高于1993—2015年的。

    進(jìn)一步判斷漢江源區(qū)和渭河源區(qū)氣象時(shí)間序列的趨勢(shì)顯著性,結(jié)果見表2。漢江源區(qū)的年降水量、年潛在蒸散發(fā)量序列趨勢(shì)檢驗(yàn)統(tǒng)計(jì)值為-0.54和-0.29,渭河源區(qū)的檢驗(yàn)統(tǒng)計(jì)值為-0.75和4.13。這表明漢江源區(qū)年降水量、年潛在蒸散發(fā)量序列呈不顯著減少趨勢(shì),渭河源區(qū)年降水量序列呈不顯著減少趨勢(shì),而年潛在蒸散發(fā)量序列呈極顯著增加趨勢(shì)。

    綜合分析過(guò)渡帶南北氣象水文要素趨勢(shì)變化規(guī)律,其中降水量和潛在蒸散發(fā)量分別有增加和減少趨勢(shì),而徑流量均呈減少趨勢(shì),這意味著氣候變化不是兩流域徑流量減少的唯一影響因素,人類活動(dòng)造成的下墊面改變也可能是秦嶺南北麓徑流量減少的重要原因,需要進(jìn)一步深入辨析各因子對(duì)徑流量變化的貢獻(xiàn)。

    3.2 南北過(guò)渡帶典型流域徑流突變及其與降水關(guān)系分析

    采用累積距平曲線分析兩個(gè)流域的年徑流量階段變化特征,結(jié)果如圖3所示,秦嶺南麓漢江源區(qū)的徑流量、降水量累積距平曲線均在1990年發(fā)生了明顯轉(zhuǎn)折,而北麓渭河源區(qū)的兩曲線明顯轉(zhuǎn)折分別發(fā)生在1993年和1992年。將1990年、1993年對(duì)應(yīng)的轉(zhuǎn)折點(diǎn)作為兩流域徑流量變化的突變點(diǎn),據(jù)此劃分徑流量變化的基準(zhǔn)期(人類活動(dòng)影響較弱時(shí)段)和變化期(人類活動(dòng)影響劇烈時(shí)段)。漢江源區(qū)基準(zhǔn)期(1975—1990年)的年均徑流量為13.10億m3,變化期(1991—2015年)為8.75億m3,變化期比基準(zhǔn)期減少了33.21%;渭河源區(qū)變化期(1994—2015年)年均徑流量為3.16億m3,較基準(zhǔn)期(1975—1993年)的6.70億m3減少了52.83%。過(guò)渡帶南北徑流突變時(shí)間基本一致,與馬新萍等[31]、劉登偉等[32]研究結(jié)果相同。漢江源區(qū)徑流減少量較渭河源區(qū)大,但渭河源區(qū)徑流相對(duì)減少量大于漢江源區(qū)。

    徑流突變前,漢江源區(qū)的徑流量和降水量累積距平曲線較為貼合,突變后兩者更為貼合,而渭河源區(qū)則相反,說(shuō)明突變以后,漢江源區(qū)徑流—降水關(guān)系未發(fā)生顯著變化,且受降水影響較突變前更甚,渭河源區(qū)的徑流—降水相關(guān)性明顯下降。1975—2015年漢江源區(qū)的徑流量和降水量累積距平變化具有較好的同步性和一致性,總體呈W形上下波動(dòng);而渭河源區(qū)的徑流量與降水量累積距平曲線變化趨勢(shì)總體一致,呈鐘形分布。這表明漢江源區(qū)的徑流—降水相關(guān)性較高,其徑流很大程度上受降水影響,但渭河源區(qū)的徑流除受降水變化影響外,其他影響因素同樣不可忽視。

    3.3 南北過(guò)渡帶典型流域徑流變化歸因分析

    3.3.1 對(duì)氣候因子和下墊面變化的敏感性分析

    漢江源區(qū)和渭河源區(qū)各時(shí)期的氣候因子P、PET和下墊面特征參數(shù)ω及其彈性系數(shù)見表3。與基準(zhǔn)期相比,兩個(gè)源區(qū)變化期的P和R均減小,PET和ω均增大,且渭河源區(qū)更加明顯;兩個(gè)源區(qū)的干旱指數(shù)PET/P均增大,徑流系數(shù)R/P均減小,R與P呈正相關(guān),但與PET、ω呈負(fù)相關(guān);漢江源區(qū)和渭河源區(qū)的干旱指數(shù)均值分別為0.83和1.72,渭河源區(qū)干旱指數(shù)>1>漢江源區(qū)干旱指數(shù),表明漢江源區(qū)的氣候相對(duì)濕潤(rùn),渭河源區(qū)氣候相對(duì)干旱。

    漢江源區(qū)徑流對(duì)P、PET和ω的彈性系數(shù)大小關(guān)系為|εP|>|εPET|>|εω|,其范圍分別為1.87~2.16、-1.16~-0.87和-0.74~-0.61,即P、PET或ω每增加1%將分別導(dǎo)致1.87%~2.16%的徑流增加,0.87%~1.16%或0.61%~0.74%的徑流減少。渭河源區(qū)三者大小關(guān)系為|εP|>|εω|>|εPET|,其范圍分別為2.60~3.18、-2.35~-1.63和-2.18~-1.60,即P、ω或PET每增加1%將分別導(dǎo)致2.60%~3.18%的徑流增加,1.63%~2.35%或1.60%~2.18%的徑流減少。另外,漢江源區(qū)|εP|與|εPET|的和小于渭河源區(qū),這表明兩流域均對(duì)P的變化最為敏感,且渭河源區(qū)徑流對(duì)氣候變化的敏感性高于漢江源區(qū)。

    表3同時(shí)反映了各時(shí)期徑流對(duì)氣候因子和下墊面敏感性的變化。變化期漢江源區(qū)的|εP|、|εPET|較基準(zhǔn)期的增量均大于|εω|的增量,而渭河源區(qū)|εω|的增量大于|εP|、|εPET|的,表明南北過(guò)渡帶兩典型流域徑流對(duì)P、PET和ω的變化的敏感性正在發(fā)生變化,其中漢江源區(qū)徑流對(duì)氣候變化的敏感性在不斷增強(qiáng),對(duì)下墊面變化的敏感性逐漸減弱,而渭河源區(qū)徑流對(duì)各因子的敏感性均逐漸增強(qiáng),尤其是對(duì)下墊面變化的敏感性。

    3.3.2 徑流變化的歸因識(shí)別

    兩流域P、PET和ω對(duì)徑流減少的貢獻(xiàn)見表4。與基準(zhǔn)期相比,在變化期,南北過(guò)渡帶兩個(gè)典型流域的P和R均減小,ω和PET均增大。圖4(d)表明漢江源區(qū)降水量對(duì)徑流變化的貢獻(xiàn)率達(dá)到了65.23%,人類活動(dòng)引起的下墊面特征變化的貢獻(xiàn)率次之,為28.75%,潛在蒸散發(fā)量的貢獻(xiàn)率僅占6.02%;在渭河源區(qū),人類活動(dòng)的貢獻(xiàn)率最大,為48.87%,由P和PET引起的氣候變化貢獻(xiàn)率約占1/2,這與左德鵬等[33]研究結(jié)果相同,且PET的貢獻(xiàn)略大于P。

    以上研究表明,漢江源區(qū)降水是造成徑流變化的主導(dǎo)因子,且貢獻(xiàn)率為65.23%,而渭河源區(qū)徑流變化與人類活動(dòng)密切相關(guān),下墊面變化對(duì)徑流減少的貢獻(xiàn)率接近50%。由此可知,南北過(guò)渡帶典型流域間徑流變化差異主導(dǎo)因子不同,對(duì)徑流減少貢獻(xiàn)也不同。

    4 討 論

    對(duì)比南北過(guò)渡帶秦嶺南北麓河川徑流變化發(fā)現(xiàn),渭河源區(qū)徑流量減少顯著,而漢江源區(qū)徑流量未發(fā)生明顯減少。兩個(gè)流域的徑流量在20世紀(jì)90年代初發(fā)生突變,該時(shí)期處于枯水期,其主要原因是北半球大尺度的氣候波動(dòng)[31]。

    對(duì)比而言,漢江源區(qū)的徑流量遠(yuǎn)大于渭河源區(qū),但渭河源區(qū)徑流量減少趨勢(shì)顯著性高于漢江源區(qū)。一方面降水是徑流的主要來(lái)源,過(guò)渡帶秦嶺南麓的降水量明顯高于北麓,潛在蒸散發(fā)量南麓低于北麓;另一方面漢江源區(qū)多山地陡坡,土層淺薄,徑流可快速匯集至出水口,茂密的植被有效減少了蒸散發(fā)消耗,而位于黃土高原邊緣的渭河源區(qū)地形平緩開闊,尤其是渭河左岸地區(qū),其土壤質(zhì)地疏松多孔不易蓄水,導(dǎo)致徑流匯集緩慢,且被深厚的黃土覆蓋,徑流匯集過(guò)程中下滲和蒸散發(fā)消耗較多。兩個(gè)源區(qū)氣候變化亦不同,漢江源區(qū)氣候因子變化趨勢(shì)不顯著,渭河源區(qū)潛在蒸散發(fā)量呈極顯著增大,降水呈不顯著減少且年潛在蒸散發(fā)量遠(yuǎn)大于該區(qū)年降水量,其徑流變化受蒸發(fā)能力直接影響[33]。

    河川徑流量的減少不僅受氣候波動(dòng)影響,還與人類活動(dòng)引起的下墊面變化關(guān)系密切。據(jù)統(tǒng)計(jì),1985—2008年漢江源區(qū)森林覆蓋率提高了22.5%。退耕還林全面實(shí)施后,漢江源區(qū)2000—2012年的土地利用變化較整個(gè)漢江流域相對(duì)緩慢,其間林地面積增加約3%,草地、水域及其他用地變化很小。此外,漢江源區(qū)無(wú)大型水利工程,少量中、小型水利工程對(duì)該流域水文的影響并不顯著。渭河流域?qū)嵤┝舜竺娣e的坡改梯、退耕還林還草等水土保持措施,使得下墊面環(huán)境發(fā)生顯著變化。這些措施一方面攔蓄地表徑流,增加了降水的入滲和蒸發(fā),另一方面隨著林地面積增加,林冠層和地被物對(duì)雨水的截持量明顯增加,使得徑流減少[34]。截至2009年底,武山縣總計(jì)退耕還林35 933 hm2。除受水土保持措施和退耕還林還草工程的影響,農(nóng)業(yè)灌溉引蓄水等水利水保工程(淤地壩、水庫(kù)、水窖)對(duì)徑流的減少起到很大的作用[11]。1996年起陜甘地區(qū)實(shí)施“集雨窖灌”工程,每年集蓄雨水約600萬(wàn)m3,對(duì)比漢江源區(qū),渭河源區(qū)土地利用覆被變化、農(nóng)業(yè)灌溉、水土保持與水利工程實(shí)施及快速城市化對(duì)河川徑流變化的影響更為劇烈,是該區(qū)徑流減少的主要因素。

    5 結(jié) 論

    以南北過(guò)渡帶秦嶺南麓漢江源區(qū)和北麓渭河源區(qū)為研究區(qū),采用線性趨勢(shì)回歸和累積距平等方法分析其徑流變化特征,運(yùn)用彈性系數(shù)法量化氣候和人類活動(dòng)對(duì)南北過(guò)渡帶秦嶺南、北麓徑流變化的貢獻(xiàn),主要結(jié)論如下。

    (1)南北過(guò)渡帶的漢江源區(qū)和渭河源區(qū)1975—2015年徑流量均呈減少趨勢(shì),且渭河源區(qū)減少趨勢(shì)更加顯著。漢江源區(qū)年徑流量減少速率為0.070億m3/a,渭河源區(qū)為0.128億m3/a;受氣候干濕交替影響,兩流域年徑流量分別在1990年和1993年發(fā)生突變,突變后年均徑流量分別減少了4.35億m3和3.54億m3。

    (2)過(guò)渡帶南北氣候差異明顯,渭河源區(qū)干旱指數(shù)>1>漢江源區(qū)干旱指數(shù)。采用彈性系數(shù)法計(jì)算發(fā)現(xiàn)兩流域?qū)邓兓顬槊舾校己釉磪^(qū)徑流對(duì)氣候變化的敏感性高于漢江源區(qū)。漢江源區(qū)徑流對(duì)氣候和下墊面變化的敏感性分別呈不顯著增強(qiáng)和減弱趨勢(shì),而渭河源區(qū)對(duì)兩者的敏感性均呈顯著增強(qiáng)趨勢(shì),尤其是下墊面變化。

    (3)漢江源區(qū)降水、潛在蒸散發(fā)、下墊面特征對(duì)徑流變化貢獻(xiàn)率分別為65.23%、6.02%、28.75%,而渭河源區(qū)三者的貢獻(xiàn)率分別為24.51%、26.62%、48.87%。由此可見,漢江源區(qū)徑流受氣候變化影響更為顯著,而渭河源區(qū)人類活動(dòng)的貢獻(xiàn)率接近50%,原因一方面是大規(guī)模的水土保持措施顯著減少了地表徑流,另一方面是農(nóng)業(yè)灌溉、雨水集蓄等貢獻(xiàn)較大。

    參考文獻(xiàn):

    [1] 趙芳,張百平,朱連奇,等.秦巴山地垂直帶譜結(jié)構(gòu)的空間分異與暖溫帶—亞熱帶界線問(wèn)題[J].地理學(xué)報(bào),2019,74(5):889-901.

    [2] 張百平.中國(guó)南北過(guò)渡帶研究的十大科學(xué)問(wèn)題[J].地理科學(xué)進(jìn)展,2019,38(3):305-311.

    [3] 王浩,劉家宏.引漢濟(jì)渭工程在國(guó)家水資源戰(zhàn)略布局中的作用[J].中國(guó)水利,2015(14):47-50.

    [4] 閆寶偉,郭生練,肖義.南水北調(diào)中線水源區(qū)與受水區(qū)降水豐枯遭遇研究[J].水利學(xué)報(bào),2007,38(10):1178-1185.

    [5] MILLIMAN J D, FARNSWORTH K L, JONES P D, et al. Climatic and Anthropogenic Factors Affecting River Discharge to the Global Ocean, 1951-2000[J]. Global and Planetary Change, 2008,62(3-4):187-194.

    [6] 王國(guó)慶,張建云,劉九夫,等.氣候變化對(duì)水文水資源影響研究綜述[J].中國(guó)水利,2008(2):55-58.

    [7] 柯新月,汪妮.秦嶺南北典型流域徑流變化規(guī)律的對(duì)比研究[J].西安理工大學(xué)學(xué)報(bào),2019,35(4):452-458.

    [8] 李凌程,張利平,夏軍,等.氣候波動(dòng)和人類活動(dòng)對(duì)南水北調(diào)中線工程典型流域徑流影響的定量評(píng)估[J].氣候變化研究進(jìn)展,2014,10(2):118-126.

    [9] CHEN Z S, CHEN Y N, LI B F. Quantifying the Effects of Climate Variability and Human Activities on Runoff for Kaidu River Basin in Arid Region of Northwest China[J]. Theoretical & Applied Climatology, 2013,111(3-4):537-545.

    [10] LI E H, MU X M, ZHAO G J, et al. Variation of Runoff and Precipitation in the HekouLongmen Region of the Yellow River Based on Elasticity Analysis[J]. The Scientific World Journal, 2014(3):929858.

    [11] ZHAO G J, MU X M, TIAN P, et al. Climate Changes and Their Impacts on Water Resources in Semiarid Regions: a Case Study of the Wei River Basin, China[J]. Hydrological Processes,2013, 27(26):3852-3863.

    [12] 冉大川,吳永紅,李雪梅,等.河龍區(qū)間近期人類活動(dòng)減水減沙貢獻(xiàn)率分析[J].人民黃河,2012,34(2):84-86.

    [13] 宋曉猛,張建云,占車生,等.氣候變化和人類活動(dòng)對(duì)水文循環(huán)影響研究進(jìn)展[J].水利學(xué)報(bào),2013,44(7):779-790.

    [14] 葉許春,張奇,劉健,等.氣候變化和人類活動(dòng)對(duì)鄱陽(yáng)湖流域徑流變化的影響研究[J].冰川凍土,2009,31(5):835-842.

    [15] 董磊華,熊立華,于坤霞,等.氣候變化與人類活動(dòng)對(duì)水文影響的研究進(jìn)展[J].水科學(xué)進(jìn)展,2012,23(2):278-285.

    [16] ALLEN R G, PEREIRA L S, RAES D, et al. Crop EvapotranspirationGuidelines for Computing Crop Water Requirements[G]. FAO Irrigation and Drainage. Rome: FAO, 1998:17-30.

    [17] 胡慧杰,崔凱,曹茜,等.黃河近百年徑流演變特征分析[J].人民黃河,2019,41(9):14-19.

    [18] 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.

    [19] 田小靖,趙廣舉,穆興民,等.水文序列突變點(diǎn)識(shí)別方法比較研究[J].泥沙研究,2019,44(2):33-40.

    [20] 張麗梅,趙廣舉,穆興民,等.基于Budyko假設(shè)的渭河徑流變化歸因識(shí)別[J].生態(tài)學(xué)報(bào),2017,38(21):7607-7617.

    [21] NING T T, LI Z, LIU W Z. Separating the Impacts of Climate Change and Land Surface Alteration on Runoff Reduction in the Jing River Catchment of China[J]. Catena,2016(147):80-86.

    [22] 張丹,梁康,聶茸,等.基于Budyko假設(shè)的流域蒸散發(fā)估算及其對(duì)氣候與下墊面的敏感性分析[J].資源科學(xué),2016,38(6):1140-1148.

    [23] 孫悅,李棟梁,朱擁軍.渭河徑流變化及其對(duì)氣候變化與人類活動(dòng)的響應(yīng)研究進(jìn)展[J].干旱氣象,2013,31(2):396-405.

    [24] CHOUDHURY B. Evaluation of an Empirical Equation for Annual Evaporation Using Field Observations and Results from a Biophysical Model[J]. Journal of Hydrology,1999,216(1-2):99-110.

    [25] YANG H B, YANG D W, LEI Z D, et al. New Analytical Derivation of the Mean Annual WaterEnergy Balance Equation[J]. Water Resources Research,2008,44(3):131-139.

    [26] BUDYKO M I. Climate and Life[M]. New York: Academic Press, 1974:55-80.

    [27] XU X Y, YANG D W, YANG H B, et al. Attribution Analysis Based on the Budyko Hypothesis for Detecting the Dominant Cause of Runoff Decline in Haihe Basin[J]. Journal of Hydrology,2014,510:530-540.

    [28] 楊大文,張樹磊,徐翔宇.基于水熱耦合平衡方程的黃河流域徑流變化歸因分析[J].中國(guó)科學(xué):技術(shù)科學(xué),2015,45(10):22-32.

    [29] 孫福寶,楊大文,劉志雨,等.基于Budyko假設(shè)的黃河流域水熱耦合平衡規(guī)律研究[J].水利學(xué)報(bào),2007,38(4):409-416.

    [30] SCHAAKE J C. From Climate to Flow[M]. New York: John Wiley and Sons Inc, 1990:65-94.

    [31] 馬新萍,白紅英,馮海鵬,等.52 a來(lái)秦嶺南北徑流變化對(duì)比及影響因素[J].干旱區(qū)地理,2013,36(6):1032-1040.

    [32] 劉登偉,延軍平.秦嶺南北徑流變化特征對(duì)比分析[J].干旱區(qū)資源與環(huán)境,2008,22(11):62-67.

    [33] 左德鵬,徐宗學(xué),隋彩虹,等.氣候變化和人類活動(dòng)對(duì)渭河流域徑流的影響[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,49(2):115-123.

    [34] 劉曉燕,劉昌明,楊勝天,等.基于遙感的黃土高原林草植被變化對(duì)河川徑流的影響分析[J].地理學(xué)報(bào),2014,69(11):1595-1603.

    【責(zé)任編輯 張 帥】

    猜你喜歡
    過(guò)渡帶下墊面源區(qū)
    不同下墊面對(duì)氣溫的影響
    冬小麥蒸散源區(qū)代表性分析
    核磁共振方法研究油水過(guò)渡帶驅(qū)替特征
    北京與成都城市下墊面閃電時(shí)空分布特征對(duì)比研究
    小興安嶺森林濕地不同過(guò)渡帶土壤微生物群落結(jié)構(gòu)研究
    森林工程(2018年3期)2018-06-26 03:40:44
    流域下墊面變化對(duì)潮白河密云水庫(kù)上游徑流影響分析
    渭河源區(qū)徑流量變化特征及趨勢(shì)分析
    下墊面變化對(duì)徑流及洪水影響分析
    《小斷塊油藏油水過(guò)渡帶提高采收率技術(shù)研究》通過(guò)中檢
    川西坳陷中段盆山過(guò)渡帶構(gòu)造變形特征及其演化
    斷塊油氣田(2012年6期)2012-03-25 09:53:54

    在线观看免费高清a一片| 欧美成人一区二区免费高清观看| 久久精品久久精品一区二区三区| 亚洲av一区综合| 日韩不卡一区二区三区视频在线| 久久久久久久大尺度免费视频| 大香蕉97超碰在线| 欧美激情在线99| 国产欧美另类精品又又久久亚洲欧美| 少妇丰满av| 精华霜和精华液先用哪个| 亚洲av国产av综合av卡| 大话2 男鬼变身卡| 肉色欧美久久久久久久蜜桃 | 久久久国产一区二区| 亚洲内射少妇av| 有码 亚洲区| 亚洲最大成人av| 好男人在线观看高清免费视频| 国产91av在线免费观看| 国产欧美亚洲国产| 亚洲国产av新网站| 丝袜美腿在线中文| 日本一二三区视频观看| 啦啦啦在线观看免费高清www| 日本三级黄在线观看| 亚洲精品日本国产第一区| 日本爱情动作片www.在线观看| 舔av片在线| 免费av毛片视频| 亚洲伊人久久精品综合| 一区二区av电影网| 少妇人妻 视频| 亚洲国产av新网站| 亚洲精品色激情综合| 少妇的逼好多水| 国产精品爽爽va在线观看网站| 日韩强制内射视频| 国产高清国产精品国产三级 | 亚洲欧洲日产国产| 亚洲精品第二区| 午夜免费观看性视频| 亚洲av日韩在线播放| 丰满人妻一区二区三区视频av| 18禁在线无遮挡免费观看视频| 欧美极品一区二区三区四区| av天堂中文字幕网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久欧美国产精品| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久午夜电影| 性色av一级| 欧美老熟妇乱子伦牲交| 日本色播在线视频| 一个人观看的视频www高清免费观看| 国产亚洲一区二区精品| 日本熟妇午夜| www.色视频.com| 大片免费播放器 马上看| 中文字幕免费在线视频6| 国产一区二区三区综合在线观看 | 我的女老师完整版在线观看| av网站免费在线观看视频| 亚洲精品乱码久久久v下载方式| 直男gayav资源| 免费人成在线观看视频色| 欧美+日韩+精品| 大香蕉久久网| 国产老妇女一区| 亚洲精品自拍成人| 日韩一本色道免费dvd| 亚洲精品自拍成人| 少妇人妻一区二区三区视频| 国内揄拍国产精品人妻在线| 九九久久精品国产亚洲av麻豆| 亚洲精品色激情综合| 男女国产视频网站| 国产成人a∨麻豆精品| 色视频在线一区二区三区| 国产国拍精品亚洲av在线观看| 一级片'在线观看视频| 国产成人福利小说| 韩国av在线不卡| av在线亚洲专区| 2021天堂中文幕一二区在线观| 久久ye,这里只有精品| 亚洲精品第二区| 97超碰精品成人国产| 亚洲真实伦在线观看| 亚洲av电影在线观看一区二区三区 | 亚洲精品成人av观看孕妇| 久久久精品免费免费高清| 久久久久久九九精品二区国产| 国产亚洲午夜精品一区二区久久 | 97人妻精品一区二区三区麻豆| 嫩草影院新地址| 特大巨黑吊av在线直播| 精品国产露脸久久av麻豆| 成人国产av品久久久| 亚洲成人av在线免费| www.色视频.com| 国产亚洲精品久久久com| 国产一区有黄有色的免费视频| 成年女人看的毛片在线观看| 夫妻午夜视频| 免费看日本二区| 九草在线视频观看| av免费在线看不卡| 亚洲人与动物交配视频| 欧美激情久久久久久爽电影| 精品国产三级普通话版| 毛片女人毛片| 午夜视频国产福利| 你懂的网址亚洲精品在线观看| 黄片wwwwww| 一区二区三区四区激情视频| 久久久a久久爽久久v久久| 国产伦理片在线播放av一区| 在线观看av片永久免费下载| 我要看日韩黄色一级片| 毛片女人毛片| 人妻 亚洲 视频| 婷婷色麻豆天堂久久| 国产综合精华液| 香蕉精品网在线| 女人久久www免费人成看片| 精品久久久久久电影网| 国产人妻一区二区三区在| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 国产在视频线精品| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 嫩草影院入口| 日韩电影二区| 日韩中字成人| 日韩亚洲欧美综合| 国产伦理片在线播放av一区| 亚洲av.av天堂| 国产片特级美女逼逼视频| 亚洲av中文字字幕乱码综合| 黄片wwwwww| 18禁裸乳无遮挡免费网站照片| 联通29元200g的流量卡| 欧美日韩精品成人综合77777| 丰满乱子伦码专区| 国产精品精品国产色婷婷| 免费观看av网站的网址| 2021天堂中文幕一二区在线观| 国产成年人精品一区二区| 亚洲精品影视一区二区三区av| 狂野欧美白嫩少妇大欣赏| 在线免费观看不下载黄p国产| 欧美丝袜亚洲另类| 涩涩av久久男人的天堂| 久久久久久久久久久免费av| 五月伊人婷婷丁香| 超碰av人人做人人爽久久| 国产精品久久久久久精品古装| 亚洲无线观看免费| 三级国产精品欧美在线观看| 一本色道久久久久久精品综合| 在线观看国产h片| 美女cb高潮喷水在线观看| 国产欧美日韩一区二区三区在线 | 国产色婷婷99| 精品少妇黑人巨大在线播放| 日韩人妻高清精品专区| 国产午夜精品久久久久久一区二区三区| 日韩制服骚丝袜av| 亚洲欧洲日产国产| 在线观看免费高清a一片| 亚洲欧美日韩另类电影网站 | 美女高潮的动态| 精品视频人人做人人爽| 亚洲第一区二区三区不卡| 麻豆国产97在线/欧美| 亚洲精品国产av成人精品| 成年免费大片在线观看| 亚洲丝袜综合中文字幕| 黄色一级大片看看| 神马国产精品三级电影在线观看| 久久久久久伊人网av| 国产高清有码在线观看视频| 国产亚洲精品久久久com| 成年av动漫网址| 亚洲成色77777| 久久人人爽人人爽人人片va| 女的被弄到高潮叫床怎么办| 国产爱豆传媒在线观看| 一级毛片我不卡| 亚洲av中文字字幕乱码综合| 色吧在线观看| 久久99热6这里只有精品| 人体艺术视频欧美日本| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 国产极品天堂在线| 最近手机中文字幕大全| 伊人久久精品亚洲午夜| 男女国产视频网站| av黄色大香蕉| 成人黄色视频免费在线看| 亚洲精品日韩av片在线观看| 亚洲va在线va天堂va国产| 亚洲精品一二三| 亚洲怡红院男人天堂| 欧美xxxx性猛交bbbb| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 亚洲经典国产精华液单| 青青草视频在线视频观看| 日日啪夜夜爽| 一二三四中文在线观看免费高清| a级毛色黄片| 在线观看人妻少妇| 成人高潮视频无遮挡免费网站| 日韩一本色道免费dvd| 国产av不卡久久| 久久精品夜色国产| 18+在线观看网站| 亚洲电影在线观看av| 国产亚洲5aaaaa淫片| 久久久成人免费电影| 91在线精品国自产拍蜜月| 婷婷色av中文字幕| 日日摸夜夜添夜夜添av毛片| 日日啪夜夜撸| 美女xxoo啪啪120秒动态图| 在线亚洲精品国产二区图片欧美 | 好男人在线观看高清免费视频| 日本熟妇午夜| 搡老乐熟女国产| 亚洲国产精品国产精品| 精品久久久久久久久亚洲| 中文乱码字字幕精品一区二区三区| 国产一区二区三区综合在线观看 | 亚洲欧美一区二区三区国产| 人人妻人人澡人人爽人人夜夜| 男女国产视频网站| 国产色婷婷99| 搡女人真爽免费视频火全软件| 亚洲欧美日韩卡通动漫| 久热这里只有精品99| 久久精品久久久久久噜噜老黄| 久久精品夜色国产| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 日韩欧美一区视频在线观看 | 日本午夜av视频| 久久久久精品性色| 制服丝袜香蕉在线| 精品久久久久久久末码| 国产欧美亚洲国产| 噜噜噜噜噜久久久久久91| 精品人妻一区二区三区麻豆| 精品久久久久久久久亚洲| 狂野欧美白嫩少妇大欣赏| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 在线播放无遮挡| 亚洲av福利一区| 国产91av在线免费观看| av国产免费在线观看| 国产黄色视频一区二区在线观看| 欧美日韩综合久久久久久| 观看美女的网站| 免费人成在线观看视频色| 亚洲综合色惰| 99久久人妻综合| 国产真实伦视频高清在线观看| 精品一区二区三区视频在线| 亚洲电影在线观看av| 91久久精品国产一区二区成人| 波多野结衣巨乳人妻| 18禁裸乳无遮挡动漫免费视频 | xxx大片免费视频| 午夜精品国产一区二区电影 | 老司机影院成人| 色综合色国产| 91在线精品国自产拍蜜月| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 香蕉精品网在线| 午夜福利在线在线| 97精品久久久久久久久久精品| 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 久久久久国产网址| 午夜亚洲福利在线播放| 99re6热这里在线精品视频| 亚洲国产欧美人成| 国产伦精品一区二区三区四那| 日韩大片免费观看网站| 国产精品久久久久久久久免| 国产老妇女一区| 国产男人的电影天堂91| 一区二区三区乱码不卡18| 欧美成人精品欧美一级黄| 久久精品国产亚洲av涩爱| 五月伊人婷婷丁香| 欧美xxxx黑人xx丫x性爽| 亚洲天堂国产精品一区在线| 欧美最新免费一区二区三区| 亚洲国产av新网站| 五月天丁香电影| 国产成人一区二区在线| 亚洲美女视频黄频| 亚洲精品日本国产第一区| 亚洲精品一二三| www.色视频.com| 国产精品嫩草影院av在线观看| 乱码一卡2卡4卡精品| 国产爱豆传媒在线观看| 校园人妻丝袜中文字幕| 内射极品少妇av片p| 国产黄a三级三级三级人| 国产一区有黄有色的免费视频| 草草在线视频免费看| 特大巨黑吊av在线直播| 看非洲黑人一级黄片| 欧美3d第一页| 亚洲欧美日韩无卡精品| 超碰97精品在线观看| 亚洲真实伦在线观看| 最近2019中文字幕mv第一页| 天天躁夜夜躁狠狠久久av| 好男人在线观看高清免费视频| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 国产高清不卡午夜福利| 久久韩国三级中文字幕| .国产精品久久| 免费观看无遮挡的男女| 看十八女毛片水多多多| 亚洲av一区综合| 国产日韩欧美亚洲二区| 白带黄色成豆腐渣| 熟女人妻精品中文字幕| 人人妻人人爽人人添夜夜欢视频 | 一级毛片久久久久久久久女| 国产 一区 欧美 日韩| 日产精品乱码卡一卡2卡三| 国产成人免费无遮挡视频| 国产综合懂色| 国产精品国产三级专区第一集| 女人十人毛片免费观看3o分钟| 一级爰片在线观看| 自拍偷自拍亚洲精品老妇| 亚洲婷婷狠狠爱综合网| 亚洲欧美成人精品一区二区| 国产精品久久久久久精品电影| 99热全是精品| 只有这里有精品99| 建设人人有责人人尽责人人享有的 | 内射极品少妇av片p| 美女xxoo啪啪120秒动态图| 日韩在线高清观看一区二区三区| 美女xxoo啪啪120秒动态图| 在线播放无遮挡| 激情 狠狠 欧美| 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品乱久久久久久| 欧美97在线视频| 在线观看三级黄色| 国产精品蜜桃在线观看| 少妇人妻 视频| 国产久久久一区二区三区| 视频区图区小说| 精品99又大又爽又粗少妇毛片| 国内揄拍国产精品人妻在线| av在线app专区| 日本wwww免费看| 亚洲图色成人| 亚洲av成人精品一二三区| 亚洲精品乱码久久久v下载方式| 国产精品秋霞免费鲁丝片| 亚洲av在线观看美女高潮| 日本wwww免费看| 国产在线一区二区三区精| 99热网站在线观看| 国产一区有黄有色的免费视频| 国产欧美另类精品又又久久亚洲欧美| 欧美潮喷喷水| 国产亚洲最大av| 日韩在线高清观看一区二区三区| 午夜亚洲福利在线播放| 国产精品.久久久| 一个人观看的视频www高清免费观看| 网址你懂的国产日韩在线| 性插视频无遮挡在线免费观看| av天堂中文字幕网| 日韩亚洲欧美综合| 男女边摸边吃奶| 久久人人爽人人片av| 国产亚洲5aaaaa淫片| 久久人人爽人人片av| 在线天堂最新版资源| 欧美日韩在线观看h| 亚洲真实伦在线观看| 色吧在线观看| 丝袜美腿在线中文| 热re99久久精品国产66热6| 秋霞伦理黄片| 精品久久久久久久久亚洲| 国产精品久久久久久精品电影小说 | 国产精品99久久99久久久不卡 | a级毛片免费高清观看在线播放| 小蜜桃在线观看免费完整版高清| 男女国产视频网站| 久久精品国产亚洲av涩爱| 日日摸夜夜添夜夜爱| 三级男女做爰猛烈吃奶摸视频| 男插女下体视频免费在线播放| 亚洲第一区二区三区不卡| 麻豆成人午夜福利视频| 免费看a级黄色片| 国产精品伦人一区二区| 精品国产三级普通话版| 国产成人a区在线观看| 国产成人福利小说| 伦精品一区二区三区| 日日摸夜夜添夜夜爱| 最近中文字幕2019免费版| 国产亚洲精品久久久com| 男的添女的下面高潮视频| 99热6这里只有精品| 少妇人妻一区二区三区视频| 国产在线一区二区三区精| 黑人高潮一二区| 日韩一区二区视频免费看| 听说在线观看完整版免费高清| 久久久久久久久久成人| 一区二区三区免费毛片| 欧美精品一区二区大全| 国产视频内射| 久久久久精品久久久久真实原创| 日本一本二区三区精品| 噜噜噜噜噜久久久久久91| 伊人久久精品亚洲午夜| 18禁裸乳无遮挡免费网站照片| 国产免费一级a男人的天堂| 国产黄色视频一区二区在线观看| 精品人妻视频免费看| 国产91av在线免费观看| av又黄又爽大尺度在线免费看| 一级毛片aaaaaa免费看小| 丰满乱子伦码专区| 欧美高清成人免费视频www| 久久99热这里只频精品6学生| 久久女婷五月综合色啪小说 | 一个人看的www免费观看视频| 亚洲va在线va天堂va国产| 成人毛片a级毛片在线播放| 青春草国产在线视频| 在线观看一区二区三区| 黄色欧美视频在线观看| 超碰av人人做人人爽久久| 国产黄片美女视频| 国产片特级美女逼逼视频| 亚洲av国产av综合av卡| 亚洲av不卡在线观看| 亚洲国产av新网站| 欧美潮喷喷水| 免费观看性生交大片5| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影小说 | 3wmmmm亚洲av在线观看| 国产亚洲av片在线观看秒播厂| 日韩av免费高清视频| av在线亚洲专区| 欧美zozozo另类| 亚洲第一区二区三区不卡| 免费在线观看成人毛片| 日韩精品有码人妻一区| 天堂网av新在线| 一级二级三级毛片免费看| 国产黄片美女视频| 国产 精品1| 婷婷色综合大香蕉| 精品久久久噜噜| 国产精品精品国产色婷婷| 白带黄色成豆腐渣| 简卡轻食公司| 国产一区二区在线观看日韩| 国产免费又黄又爽又色| 我要看日韩黄色一级片| h日本视频在线播放| 麻豆成人午夜福利视频| 亚洲精品成人久久久久久| 真实男女啪啪啪动态图| 99九九线精品视频在线观看视频| 国产高清不卡午夜福利| 97精品久久久久久久久久精品| 看非洲黑人一级黄片| 天堂网av新在线| 超碰av人人做人人爽久久| 人体艺术视频欧美日本| 人妻系列 视频| 亚洲欧洲国产日韩| 免费人成在线观看视频色| 校园人妻丝袜中文字幕| 99视频精品全部免费 在线| 日本欧美国产在线视频| 人妻夜夜爽99麻豆av| 在线 av 中文字幕| kizo精华| 国产 一区精品| 国国产精品蜜臀av免费| 午夜免费男女啪啪视频观看| 久久人人爽av亚洲精品天堂 | 欧美人与善性xxx| 校园人妻丝袜中文字幕| 女人久久www免费人成看片| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 日本欧美国产在线视频| 亚洲欧美一区二区三区黑人 | 欧美日韩一区二区视频在线观看视频在线 | 国产精品99久久99久久久不卡 | 菩萨蛮人人尽说江南好唐韦庄| 伊人久久国产一区二区| 免费电影在线观看免费观看| 青春草视频在线免费观看| 久久久久久国产a免费观看| 国产免费又黄又爽又色| 亚洲欧洲日产国产| 日韩一区二区视频免费看| 国产精品三级大全| 日韩欧美一区视频在线观看 | 中文在线观看免费www的网站| 丰满人妻一区二区三区视频av| 亚洲av成人精品一区久久| 一级片'在线观看视频| 久热这里只有精品99| 久久精品综合一区二区三区| 国产探花极品一区二区| 国产亚洲精品久久久com| freevideosex欧美| av线在线观看网站| 亚洲精品国产成人久久av| 国产男女超爽视频在线观看| 久久人人爽人人片av| 日日撸夜夜添| 欧美人与善性xxx| 亚洲av免费高清在线观看| 免费观看a级毛片全部| 夜夜爽夜夜爽视频| 丝袜美腿在线中文| 亚洲精品亚洲一区二区| 中文字幕人妻熟人妻熟丝袜美| 一级二级三级毛片免费看| 亚洲电影在线观看av| 嘟嘟电影网在线观看| 2018国产大陆天天弄谢| 国产毛片a区久久久久| 国产人妻一区二区三区在| 搡女人真爽免费视频火全软件| 国产精品国产av在线观看| 欧美日韩国产mv在线观看视频 | 久久ye,这里只有精品| 国产熟女欧美一区二区| 一级av片app| 一边亲一边摸免费视频| 亚洲经典国产精华液单| 麻豆精品久久久久久蜜桃| 亚州av有码| 亚洲国产日韩一区二区| 下体分泌物呈黄色| 精品一区二区免费观看| 欧美日本视频| 大香蕉97超碰在线| 26uuu在线亚洲综合色| 成人一区二区视频在线观看| 欧美zozozo另类| 亚洲精品自拍成人| 老女人水多毛片| 国产黄片视频在线免费观看| 在线观看av片永久免费下载| av黄色大香蕉| 校园人妻丝袜中文字幕| 精品人妻偷拍中文字幕| 亚洲国产精品国产精品| 色视频在线一区二区三区| av网站免费在线观看视频| 久久精品国产a三级三级三级| 深爱激情五月婷婷| 麻豆成人av视频| 卡戴珊不雅视频在线播放| 亚洲精品乱码久久久v下载方式| 交换朋友夫妻互换小说| 在线免费观看不下载黄p国产| 少妇丰满av| 久久久a久久爽久久v久久| 国国产精品蜜臀av免费| 国产成人一区二区在线| 毛片一级片免费看久久久久| 国产精品久久久久久久电影| 大香蕉97超碰在线| 边亲边吃奶的免费视频| 美女高潮的动态| 亚洲美女视频黄频| 午夜精品一区二区三区免费看| 国产精品爽爽va在线观看网站| 九草在线视频观看| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 人妻系列 视频| 国产一区二区在线观看日韩| 老师上课跳d突然被开到最大视频| 少妇的逼好多水| 久久精品国产亚洲av涩爱| 日韩免费高清中文字幕av|