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

    雅魯藏布江流域水化學(xué)時(shí)空變化及其控制因素

    2018-11-28 02:56:34劉佳駒趙雨順郭懷成
    中國環(huán)境科學(xué) 2018年11期
    關(guān)鍵詞:平水雅魯藏布江豐水期

    劉佳駒,趙雨順,黃 香,郭懷成*

    ?

    雅魯藏布江流域水化學(xué)時(shí)空變化及其控制因素

    劉佳駒1,趙雨順2,黃 香3,郭懷成1*

    (1.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,北京 100871;2.天津大學(xué)環(huán)境科學(xué)與工程學(xué)院,天津 300072;3.西藏大學(xué)化學(xué)與環(huán)境學(xué)院,西藏 拉薩 850000)

    為研究雅魯藏布江全流域水化學(xué)特征及其控制因素,2016年先后采集雅魯藏布江干流及其支流豐、平、枯三個(gè)水期水樣212個(gè),綜合運(yùn)用數(shù)理統(tǒng)計(jì)、Piper三線圖、Gibbs 模型和離子比等方法,分析了全流域河水的水文地球化學(xué)特征,并探討了流域的水化學(xué)演化規(guī)律.結(jié)果表明,河水中陽離子均以 Ca2+、Na+和Mg2+為主;陰離子以HCO3-和 SO42-為主;TDS平均值為204.51mg/L,礦化度較低;水化學(xué)類型以 HCO3·SO4(SO4·HCO3)-Ca·Mg(Mg·Ca)型水為主;雅魯藏布江主干河流內(nèi)的主要離子含量變化存在波動(dòng)且季節(jié)性變化顯著.空間上,河水中絕大多數(shù)離子遵循著先增大后減少的趨勢;水化學(xué)樣品均分布在 Gibbs 模型左中部,說明該流域水化學(xué)離子組成受巖石風(fēng)化作用控制;主成分分析及相關(guān)分析表明雅魯藏布江流域水化學(xué)組分受人為影響作用控制,方解石、白云巖的風(fēng)化及硫酸溶解也起到非常重要的作用.流域三個(gè)水期的重金屬元素的含量絕大多數(shù)能滿足地表水I類水體的要求.

    水化學(xué);雅魯藏布江;時(shí)空變化;水質(zhì)

    青藏高原被稱為亞洲水塔.其為世界近1/3人口的生活、生產(chǎn)提供淡水資源[1].由于地形、地質(zhì)和氣候等自然條件以及人類生產(chǎn)、生活等的影響,水質(zhì)受多因素、多層次和多尺度相互綜合影響,水質(zhì)的變化會(huì)對(duì)水生生態(tài)系統(tǒng)產(chǎn)生一系列的生態(tài)效應(yīng),能在很大程度上反映流域的變化特征[2].作為世界的“第三極”,青藏高原在通過陸地和海洋能量交換發(fā)生的元素地球化學(xué)循環(huán)中占有非常重要的比例[1].研究河流和湖泊水循環(huán)和水化學(xué)特征以及水資源和環(huán)境的保護(hù)具有重要的科學(xué)意義.此外,在水化學(xué)研究中揭示最主要離子的來源以及其控制因素非常重要[3].

    近年來,眾多國內(nèi)外學(xué)者開始對(duì)雅魯藏布江流域展開了較為深入的研究并取得了豐碩的成果,主要分自然影響與人為影響兩方面:隨著全球氣候變暖, 青藏高原的冰川面積減少,導(dǎo)致的水文水生態(tài)環(huán)境也在發(fā)生著變化,進(jìn)而影響了高原河流的水化學(xué)性質(zhì)[4-5].Adachi等[6]對(duì)高原東部雅魯藏布江流域主要化學(xué)組成的研究表明,巖石風(fēng)化風(fēng)化和侵蝕過程對(duì)流域水化學(xué)性質(zhì)具有強(qiáng)烈的影響;沈大軍等[7]對(duì)雅魯藏布江流域水資源進(jìn)行了評(píng)價(jià),認(rèn)為雅魯藏布江外流水系河水泥沙含量較小,河水礦化度較小,水質(zhì)良好;李紅敬等[8]通過對(duì)雅魯藏布江的水質(zhì)分析發(fā)現(xiàn),在6~10月,溶氧值低、水溫較高,12~2月,溶氧值高、水溫較低、8月懸浮物含量高,下游較中游的pH值和懸浮物含量低;Huang等[9]對(duì)雅魯藏布江全流域的水質(zhì)進(jìn)行了評(píng)估,并基于統(tǒng)計(jì)學(xué)對(duì)水化學(xué)的控制因素進(jìn)行了分析.一些學(xué)者的研究表明,雅魯藏布江的水化學(xué)特征變化一定程度上受到人類活動(dòng)的影響[10]:由于青藏高原是西風(fēng)、印度洋季風(fēng)和東亞季風(fēng)的交匯點(diǎn),大氣環(huán)流也導(dǎo)致了印度和尼泊爾等南亞國家的工農(nóng)業(yè)污染物跨界傳輸.高原上的人為活動(dòng)和氣候變化都將影響雅魯藏布江的水化學(xué)過程,并對(duì)青藏高原南部和印度和孟加拉國下游地區(qū)的人類生產(chǎn)生活構(gòu)成風(fēng)險(xiǎn)[8].Huang等[11]對(duì)雅魯藏布江流域水質(zhì)研究表明, Cd,Hg,Ni和Cu等重金屬含量相對(duì)較高,污染被認(rèn)為是由當(dāng)?shù)夭傻V業(yè)造成的.在鹽湖中,還發(fā)現(xiàn)了高濃度的重金屬元素,如鉛和砷.這可能是由于地?zé)岬茸匀贿^程引發(fā)[10].總的來說,雅魯藏布江水化學(xué)特性還沒有經(jīng)過系統(tǒng)深入研究,對(duì)全流域、季節(jié)性變化的研究較為缺乏.相對(duì)于我國其它流域,有關(guān)雅魯藏布江流域水化學(xué)特征、化學(xué)風(fēng)化過程,尤其是季節(jié)性變化及影響因素的研究尚不夠全面、細(xì)致、深入.因此,有必要繼續(xù)開展和深入研究雅魯藏布江全流域化學(xué)風(fēng)化和元素地球化學(xué)循環(huán)過程,為將來探索“第三極”以及全球氣候環(huán)境變化提供基礎(chǔ).

    本文基于2016年對(duì)我國西藏境內(nèi)雅魯藏布江全流域干支流三個(gè)水文期的河水采樣,結(jié)合已有的流域環(huán)境背景資料, 分析報(bào)道了水體樣品中主要陰陽離子的含量及其時(shí)空分布特征,探討了水化學(xué)的主控因素,進(jìn)一步研究雅江水質(zhì)和“亞洲水塔”的水質(zhì)變遷提供基礎(chǔ)數(shù)據(jù),為了解流域水環(huán)境和水資源提供科學(xué)依據(jù).

    1 研究區(qū)域概況及樣品采集與分析

    1.1 研究區(qū)概況

    圖1 青藏高原雅魯藏布江流域和采樣斷面分布

    研究區(qū)地處西藏的南部,岡底斯山脈、念青唐古拉山脈以南,喜馬拉雅山脈以北的地區(qū).雅魯藏布江發(fā)源于西藏自治區(qū)西南部喜馬拉雅山北麓的杰馬央宗冰川,自西向東流經(jīng)仲巴、日喀則、澤當(dāng)?shù)鹊?在米林縣派區(qū)附近,形成一大拐彎,折向南流,經(jīng)中國巴昔卡出境,流入印度洋的孟加拉灣.流域面積為24萬km2,沿雅魯藏布江干流谷地,海拔高度為2800~4500m[12].

    雅魯藏布江流域下游地區(qū)高溫多雨.降水的年際變化較小,年內(nèi)分配很不均勻.年均降水量從河流上游200mm左右至最下游的5000mm左右.全年平均降水主要集中在7~9月,達(dá)到全年降水量的50%~80%.年均最高溫度出現(xiàn)在6月,下游地區(qū)多出現(xiàn)在7月,年均溫度最低出現(xiàn)在1月.降水最多的月份,其冰雪融水補(bǔ)給河流的水量也大.此外,該流域還具有枯水期水量較大、懸移質(zhì)泥沙含量少、河水溫度低、河水礦化度小、總硬度低等特點(diǎn)[13].

    雅魯藏布江有5條流域面積超過10000km2的支流,如,帕隆藏布、年楚河、尼洋河、拉薩河、多雄藏布.尼洋河是雅魯藏布江的一級(jí)支流,尼洋河所處地區(qū)的氣候溫和濕潤.年均氣溫8℃左右,年降水量 600~900mm,年內(nèi)徑流量相對(duì)集中,主要集中在每年的6~9月.拉薩河位于西藏中南部全長568km,流域面積31760km2.拉薩河流域面積僅占西藏自治區(qū)總面積的2.7%,而流域內(nèi)的人口、耕地卻約占全自治區(qū)的15%,是西藏工、農(nóng)、牧業(yè)集中的地區(qū).

    1.2 水樣采集與分析

    依據(jù)《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ 494-2009)相關(guān)內(nèi)容以及雅魯藏布江流域地貌特征,并兼顧各支流匯入位置,在干流及其支流布設(shè)67個(gè)采樣點(diǎn),保證每個(gè)樣品點(diǎn)具有一定的代表性[11].采樣分別在2016年9月15日~10月8日(豐水期),2016年11月8~11月26日(平水期)和2017年4月22日~5月7日(枯水期)分3次采集樣品.如圖1所示.每個(gè)時(shí)間段樣本在半個(gè)月內(nèi)從下游到上游連續(xù)采集.為保持采樣條件的一致性,選擇24h內(nèi)沒有降雨的天數(shù)作為采集樣本日[14].

    樣品收集之前,每個(gè)樣品瓶用硝酸洗滌,之后用蒸餾水沖洗.將水樣品收集在2L聚丙烯瓶中,然后立即通過0.45mm醋酸纖維素膜過濾器.將樣品置于干凈的自密封聚乙烯(PE)塑料袋中,并保存在4℃的冰箱中直至測試.為了確保數(shù)據(jù)的準(zhǔn)確性,所有樣品瓶在收集水樣之前用收集的濾液洗滌2~3次.為避免污染和破壞,將樣品置于冰箱中冷凍運(yùn)輸.將由水,過濾器和沉淀物組成的樣品運(yùn)送到實(shí)驗(yàn)室.將水樣在4℃的冰箱中儲(chǔ)存,并在自然風(fēng)化條件下將沉淀樣品保存在無塵實(shí)驗(yàn)室中[11].

    雅魯藏布江河水基礎(chǔ)理化性質(zhì),陰陽離子及重金屬的測量儀器及方法如表1所示.在實(shí)驗(yàn)過程中,儀器由國家黃土標(biāo)準(zhǔn)物質(zhì)(GBW07408)標(biāo)定,并由標(biāo)準(zhǔn)溶液控制.陽離子和陰離子的分析精度優(yōu)于2%.總陽離子和陰離子的總和表明與使用探針測量的TDS值(Tz++Tz-=1.02TDS,2=0.98)匹配度良好,這表明數(shù)據(jù)質(zhì)量很高[11].

    所有統(tǒng)計(jì)計(jì)算均使用Excel 2007(Microsoft Office)和spss18.0進(jìn)行.實(shí)驗(yàn)結(jié)果圖使用Origin8.0制作[11].首先,通過統(tǒng)計(jì)檢驗(yàn)確定離子元素的含量是否遵循正態(tài)分布.然后,使用Pearson相關(guān)分析評(píng)估離子元素濃度的相關(guān)關(guān)系.

    表1 基礎(chǔ)理化性質(zhì)、陰陽離子及重金屬測量儀器及方法

    2 結(jié)果和討論

    2.1 河水水化學(xué)總體特征

    2.1.1 基礎(chǔ)理化性質(zhì) 通過對(duì)常規(guī)水化學(xué)參數(shù)進(jìn)行分析,雅魯藏布江流域水樣品pH值、T、EC、DO、TDS變化范圍及均值如表2所示.雅魯藏布江流域三期水體的水溫為0.0~20.5℃,平水期平均溫度最低(4.2℃);雅魯藏布江及其支流的pH值變化幅度較小,介于7.97~9.03之間,平均值為8.46,整體呈弱堿性,碳酸鹽系統(tǒng)是天然水中優(yōu)良的緩沖系統(tǒng),它對(duì)避免天然水 pH 值急劇變化起緩沖作用.豐水季節(jié)河水 pH 值(8.6)略低于枯水期(8.7)及平水期(8.8),這可能的原因是豐水期降雨導(dǎo)致水體的pH值減小;3個(gè)水期TDS變化很大,13.30~451.00mg/L,平均值為183.95mg/L,高于世界河流的平均值(120mg/L)[15](表2);電導(dǎo)率(EC)可以反映水體的離子強(qiáng)度,雅魯藏布江流域三期水體的電導(dǎo)率變化范圍均值為18.5~632.0μS/cm,均值為252.3μS/cm;在溫度和鹽度的影響下,雅魯藏布江及其支流的DO在枯水期最高(均值為10.85mg/L),在豐水期最低(均值為8.1mg/L)(表2),這可能是水體中DO的含量與空氣中的氧的分壓,水溫有相關(guān),枯水期的溫度水體較低,因此導(dǎo)致DO含量高.雅魯藏布江流域主要河道的陽離子濃度依次為Ca2+>Na+>Mg2+> K+.Ca2+為是優(yōu)勢陽離子,占陽離子總量的31.3%.陰離子濃度順序?yàn)?HCO3->SO42->Cl-.其中HCO3-和SO42-為主要陰離子,分別占陰離子總量的72.35%和14.32%.Ca2+和SO42-的平均濃度分別是全球河流平均濃度的3~4以上[15](表3).

    表2 基礎(chǔ)理化性質(zhì)

    表3 雅魯藏布江及世界大河主要離子含量統(tǒng)計(jì)(mg/L)

    2.1.2 水化學(xué)類型 河流水化學(xué)特征主要受天然補(bǔ)給來源、化學(xué)風(fēng)化以及人類活動(dòng)的影響.決定河流水化學(xué)特征的主要過程是降水對(duì)地表組成物質(zhì)的淋溶過程、主支流交匯處的混合過程、懸浮物的吸附和解析過程以及蒸發(fā)濃縮而導(dǎo)致的變質(zhì)等過程.Piper三線圖可以用來分析河水水化學(xué)成分的演化規(guī)律,結(jié)果如圖2所示.雅魯藏布江及其支流的樣品離子含量分布接近CO32-+HCO3-軸分布(圖2),陰離子主要由HCO3-和SO42-組成.在陽離子三角圖中,采樣點(diǎn)位于Ca2++Mg2+軸附近,表明該區(qū)域的陽離子主要為Ca2+和Mg2+,這可能反映了流域水化學(xué)性質(zhì)主要受碳酸鹽風(fēng)化.流域以SO4·HCO3(HCO3·SO4)- Ca·Mg(Mg·Ca)型為主.在雅魯藏布江及其支流中,80%以上的Ca2+, Mg2+和HCO3-產(chǎn)生于碳酸鹽風(fēng)化過程中.在河流的中上部,尤其是年楚河附近,Na+的比例高達(dá)29%,Cl-高達(dá)31%.說明蒸發(fā)巖的溶解是其主要控制因素.

    圖2 雅魯藏布江干流及其支流河水陰陽離子三角圖

    圖3 雅魯藏布江流域豐、平、枯3季水體重金屬濃度箱線圖

    2.1.3 重金屬分析 青藏高原近30a的金屬采礦業(yè)影響了該地區(qū)的地表水質(zhì)量.重金屬含量的增加會(huì)對(duì)漁業(yè)和飲用水安全產(chǎn)生影響.如圖3所示雅魯藏布江流域在豐、平、枯3個(gè)水期的重金屬濃度的箱線圖.研究表明,Al(均值為62.65mg/L)和Fe(均值為60.8mg/L)濃度較高.Cu,Zn,Ag,Cd和Cr的濃度普遍較低,符合中華人民共和國地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)的I類水標(biāo)準(zhǔn). 雅魯藏布江流域三個(gè)水期的As元素的含量大部分能滿足地表水I類水體的要求.但是,羊巴井地?zé)崽锵掠螛悠吩谟昙竞涂菟诘腁s濃度分別為170.5ug/L和144.0ug/L,超過了地表水環(huán)境類水V標(biāo)準(zhǔn).青藏高原沉積物中As的濃度是中國地殼豐度的11.6倍.地?zé)峄顒?dòng)釋放的氣液導(dǎo)致As含量高.其含量的超標(biāo)可能會(huì)影響該地區(qū)居民的飲水安全[18].90.2%水樣中汞的濃度達(dá)到地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)的I類水標(biāo)準(zhǔn).上游瑪玖藏布采樣點(diǎn)As濃度極高(枯水期為181.2ug/L,豐水期為217.58ug/L,平水期為182.1ug/L),這可能與采礦業(yè)以及As金屬的熔煉有關(guān)[19].但雅魯藏布江干流及其支流重金屬濃度季節(jié)變化不明顯.

    2.2 河水主要離子時(shí)空變化特征

    圖4 雅魯藏布江干流豐,平,枯三季主要離子濃度的變化

    徑流補(bǔ)給來源及氣候的顯著差異是影響季風(fēng)區(qū)河流水文、水化學(xué)季節(jié)性特征的重要因素.通過分析雅魯藏布江干流30個(gè)采樣點(diǎn)的離子濃度,結(jié)果表明,該河流的水化學(xué)特征具有顯著的季節(jié)差異,如圖4所示.雅魯藏布江干流陽、陰離子濃度變化趨勢總體表現(xiàn)為豐水期最低.豐水季節(jié)流域河水中大部分陰陽離子平均濃度(mg/L)均小于枯水期及平水期,離子濃度值下降明顯,較枯水期離子濃度減少約 20.97mg/L,較平水期減少約為50.37mg/L.(表3)陽離子變化以Ca2+為主,枯水期較平水期減小了約9.35mg/L,其余陽離子變化不明顯.流域陰離子濃度總量變化幅度遠(yuǎn)遠(yuǎn)超過了陽離子,HCO3-變化最為明顯,豐水期較平水期減少了約37.42mg/L.豐水期降水量明顯增加,而離子來源沒有明顯變化的情況下,形成的稀釋效應(yīng).相反,平水期和枯水期的流速相對(duì)較慢,水體與巖石接觸時(shí)間較長,增加了巖石的溶解.在枯水期和豐水期,雅魯藏布江干流的陽離子濃度依次為,Ca2+>Na+>Mg2 +>K+,但在平水期,表現(xiàn)為Ca2+>Mg2+>Na+>K+,這種變化在中游和下游最為明顯(圖4).在3個(gè)水期,水樣中陰離子濃度變化趨勢均為HCO3->SO42->Cl->NO3-.

    從空間上看,雅魯藏布江干流離子濃度從上游至下游含量變化存在波動(dòng),但總體呈先增后減的趨勢(圖4).最上游采樣點(diǎn)(Y17、Y16)測得的TDS濃度與整個(gè)干流相比明顯較低.這一結(jié)果可能是由于冰川和融雪稀釋了離子濃度.雅江上游地區(qū)表現(xiàn)為較為干燥的草原氣候,其蒸發(fā)量明顯超過降雨量,使得水體在一定程度上得到蒸發(fā)濃縮,大量裸地、灘涂的存在使得地表大面積暴露,導(dǎo)致當(dāng)?shù)仫L(fēng)化侵蝕作用強(qiáng)烈,大量鹽離子得以釋放進(jìn)入水體.因此,Y14和Y15采樣點(diǎn)離子濃度相比前段劇烈增加,增速最快的離子為Ca2+、Na+、HCO3-、SO42-.從源頭向下游地區(qū)過渡的過程中隨著降水增加,徑流量得以補(bǔ)充,同時(shí),灌木、針葉林、闊葉林的組合植被類型減輕了水土流失,最終導(dǎo)致了水體中的離子濃度呈現(xiàn)逐漸下降的趨勢.對(duì)比枯水季節(jié)與豐水季節(jié)流域離子空間特征分布發(fā)現(xiàn)(圖4),豐水季節(jié)陽離子總濃度值均低于枯水季節(jié).其中,上游降幅最大,約為 29.65mg/L,中游和下游降幅分別約為16.95mg/L和14.69mg/L,表明上游陽離子含量季節(jié)性差異最大,中游和下游依次降低.雅江流域從上游至下游離子濃度季節(jié)性差異呈逐步降低的趨勢.

    2.3 離子來源分析

    2.3.1 相關(guān)性分析 相關(guān)分析通常用于揭示主要離子的來源.表4顯示了雅魯藏布江流域各化學(xué)成分之間的相關(guān)矩陣.TDS與Ca2+,Mg2+,K+,Na+,HCO3-, Cl-和SO42-呈顯著正相關(guān)(<0.01),表明這些物質(zhì)都是TDS的來源.其中,TDS與Ca2+和HCO3-的相關(guān)性最顯著,相關(guān)系數(shù)分別為0.874和0.924.變化趨勢相同,統(tǒng)計(jì)學(xué)上表明Ca2+和HCO3-是TDS的主要來源.HCO3-與Ca2+和Mg2+有顯著相關(guān)性(<0.01).由于碳酸鹽和硫酸鹽礦物的風(fēng)化溶解作用,SO42-與 Ca2+、Mg2+的相關(guān)性也很強(qiáng)說明它們之間有可能來自同一來源,可能來源于碳酸鹽和硫酸鹽礦物的風(fēng)化溶解.Na+、K+都與Cl-有強(qiáng)相關(guān)性,說明它們有可能來自蒸發(fā)巖的溶解.

    表4 離子相關(guān)性分析

    2.3.2 主成分分析 為了進(jìn)一步探討雅魯藏布江流域水體主要離子與地質(zhì)背景之間的關(guān)系,了解主要離子的控制因素,應(yīng)用主成分分析法(PCA)對(duì)雅魯藏布江流域3個(gè)水期主要離子間的關(guān)系進(jìn)行分析.基于3個(gè)水期的212份水樣,通過因子分析得到3個(gè)主要公共因素.如圖5所示,3個(gè)水期的第一因子與K+,Na+,Cl-,Li,SiO2之間存在很強(qiáng)的相關(guān)性,這與高原各類蒸發(fā)巖的溶解有關(guān)[19].第二因子與Mg2+, SO42-,Ca2+和HCO3-存在較高相關(guān)性.說明方解石、白云巖的風(fēng)化及硫酸溶解的相關(guān)貢獻(xiàn)較大.其中碳酸鹽風(fēng)化作用是Ca2+的主要來源,80%~90%的溶解Ca2+和Mg2+是由雅魯藏布江流域碳酸鹽風(fēng)化產(chǎn)生的.第三個(gè)因素與Mn和Zn有較強(qiáng)相關(guān)性,表明與采礦業(yè)等人類活動(dòng)相關(guān),但人類活動(dòng)的影響并不明顯.其中流域豐水期及平水期的三個(gè)主要因素累計(jì)方差貢獻(xiàn)率為75.14%.第1,2,3因子貢獻(xiàn)率分別達(dá)到34.5%,26.2%和14.4%.枯水期的累計(jì)方差貢獻(xiàn)率為66.67%.對(duì)比平水期及豐水期可以發(fā)現(xiàn),枯水期累積方差百分比(66.67%)比平水期及豐水期(75.14%)略有降低,但離子負(fù)荷值增高,說明豐水期及平水期的主要離子來源的變化更加復(fù)雜(表5).

    圖5 雅魯藏布江流域主要離子主成分分析

    表5 雅魯藏布江流域3個(gè)水期離子和元素的主成分分析(PCA)

    圖6 雅魯藏布江流域主要離子吉布斯圖

    Y-雅魯藏布江干流,L-拉薩河,N-尼洋河,NY-年楚河,P-帕隆藏布江,D-多雄藏布江

    2.3.3 吉布斯分析 Gibbs 設(shè)計(jì)了半對(duì)數(shù)坐標(biāo)圖,縱坐標(biāo)為水樣TDS的對(duì)數(shù)坐標(biāo),橫坐標(biāo)分別為陽離子質(zhì)量濃度的比值(Na+)/(Na++Ca2+)或者陰離子質(zhì)量濃度的比值(Cl-)/(Cl-+HCO3-)[21]. Gibbs圖顯示Na+/(Na++Ca2+)(<0.38)和Cl-/(Cl-+HCO3-)(<0.21)比值較低, TDS值中等,說明雅魯藏布江干流及其支流的離子化學(xué)性質(zhì)主要是由巖石風(fēng)化控制.所采集水樣遠(yuǎn)離“降水主導(dǎo)”和“蒸發(fā)結(jié)晶”作用控制區(qū)域,且河水的Cl-/Na+平均值為0.65,遠(yuǎn)低于世界平均海水比例(Cl-/Na+=1.15),這些均表明研究區(qū)內(nèi)大氣降水和蒸發(fā)結(jié)晶對(duì)水體中離子質(zhì)量濃度的影響不是主要控制因素.

    3 結(jié)論

    雅魯藏布江流域的地表水質(zhì)呈弱堿性.河流的平均TDS為204.51mg/L,高于全球平均水平.主要陽離子是Ca2+,Na+和Mg2+,而主要陰離子是HCO3-和SO42-.流域大部分地區(qū)的水化學(xué)特征為Ca2+-HCO3-型.大部分重金屬元素均能達(dá)到中華人民共和國地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)的I類水標(biāo)準(zhǔn).然而,受到上游采礦及金屬冶煉的影響,有少數(shù)采樣點(diǎn)的As含量超標(biāo),中游羊八井地?zé)崽锔浇鼧悠稟s濃度較高,則與地?zé)峄顒?dòng)有關(guān).季節(jié)變化上,雅魯藏布江流域的水化學(xué)特征具有顯著的季節(jié)性差異,其中,豐水期的離子濃度最低,枯水期和平水期的離子濃度相對(duì)較高.在空間上,雅魯藏布江干流的離子濃度從上游到下游呈逐步降低的趨勢,支流中的離子濃度波動(dòng)顯著.主成分分析和吉布斯模型表明,雅魯藏布江流域的離子組成主要受巖石風(fēng)化作用控制.總體來說,雅魯藏布江絕大部分采樣斷面的水質(zhì)良好,極個(gè)別采樣斷面重金屬超標(biāo).在全球氣候變化的背景及人類活動(dòng)增多的雙重壓力下,雅魯藏布江的水化學(xué)特征會(huì)受到影響從而間接會(huì)對(duì)青藏高原的人類健康構(gòu)成風(fēng)險(xiǎn).本文的研究結(jié)論對(duì)進(jìn)一步研究雅江水質(zhì)和“亞洲水塔”的水質(zhì)變遷提供了基礎(chǔ),為了解流域水環(huán)境和水資源的科學(xué)管理提供基礎(chǔ)性科學(xué)依據(jù).

    [1] Li F, Zhang Y, Xu Z, et al. The impact of climate change on runoff in the southeastern Tibetan Plateau [J]. Journal of Hydrology, 2013, 505(505):188-201.

    [2] Cullaj A, Hasko A, Miho A, et al. The quality of Albanian natural waters and the human impact [J]. Environment International, 2005, 31(1):133-152.

    [3] 樂嘉祥,王德春.中國河流水化學(xué)特征[J]. 地理學(xué)報(bào), 1963,(1):3-15.

    [4] Wang S L, Jin H J, Li S X, et al. Permafrost degradation on the Qinghai-Tibet Plateau and its environmental impacts. [J]. Permafrost & Periglacial Processes, 2000,11(1):43–53.

    [5] Wu W, Xu S, Yang J, et al. Silicate weathering and CO2, consumption deduced from the seven Chinese rivers originating in the Qinghai- Tibet Plateau [J]. Chemical Geology, 2008,249(3):307-320.

    [6] Adachi K, Tainosho Y. Characterization of heavy metal particles embedded in tire dust. [J]. Environment International, 2004,30(8):1009.

    [7] 沈大軍,陳傳友.青藏高原水資源及其開發(fā)利用[J]. 自然資源學(xué)報(bào), 1996,(1):8-14.

    [8] 李紅敬,張 娜,林小濤.西藏雅魯藏布江水質(zhì)時(shí)空特征分析[J]. 河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010,38(2):126-130.

    [9] X. Huang ?, Sillanp?? M, Gjessing E T, et al. Water quality in the southern Tibetan Plateau: chemical evaluation of the Yarlung Tsangpo (Brahmaputra) [J]. River Research & Applications, 2011,27(1):113–121.

    [10] Yao T, Thompson L, Yang W, et al. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings [J]. Nature Climate Change, 2012,2(9):663-667.

    [11] Huang X, Sillanp?? M, Gjessing E T, et al. Water quality in the Tibetan Plateau: major ions and trace elements in the headwaters of four major Asian rivers [J]. Science of the Total Environment, 2009,407(24):6242

    [12] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

    [13] Huang Xiang, Sillanp??, M, Gjessing, E. T, et al. Environmental impact of mining activities on the surface water quality in Tibet: Gyama valley [J]. Science of the Total Environment, 2010,408(19):4177-4184.

    [14] Pant R R, Zhang F, Rehman F U, et al. Spatiotemporal variations of hydrogeochemistry and its controlling factors in the Gandaki River Basin, Central Himalaya Nepal [J]. Science of the Total Environment, 2017,622-623:770.

    [15] Gaillardet J, Dupré B, Louvat P, et al. Global silicate weathering and CO2, consumption rates deduced from the chemistry of large rivers [J]. Chemical Geology, 1999,159(1–4):3-30.

    [16] Sarin M M, Krishnaswami S, Dilli K, et al. Major ion chemistry of the Ganga- Brahmaputra river system: Weathering processes and fluxes to the Bay of Bengal [J]. Geochimica Et Cosmochimica Acta, 1989,53(5):997-1009.

    [17] Ahmad T, Khanna P P, Chakrapani G J, et al. Geochemical characteristics of water and sediment of the Indus river, Trans- Himalaya, India: constraints on weathering and erosion [J]. Journal of Asian Earth Sciences, 1998,16(2/3):333-346.

    [18] Thomas J, Joseph S, Thrivikramji K P, et al. Seasonal variation in major ion chemistry of a tropical mountain river, the southern Western Ghats, Kerala, India [J]. Environmental Earth Sciences, 2014,71(5): 2333-2351.

    [19] Dekov V M, Komy Z, Araújo F, et al. Chemical composition of sediments, suspended matter, river water and ground water of the Nile (Aswan-Sohag traverse). [J]. Science of the Total Environment, 1997,201(3):195.

    [20] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

    [21] Gibbs R J. Mechanisms Controlling World Water Chemistry [J]. Science, 1971,172(3985):870-2.

    [22] 劉 劍,姚治君,陳傳友.雅魯藏布江徑流變化趨勢及原因分析[J]. 自然資源學(xué)報(bào), 2007,(3):471-477.

    [23] 楊志剛,卓 瑪,路紅亞,等.1961~2010年西藏雅魯藏布江流域降水量變化特征及其對(duì)徑流的影響分析[J]. 冰川凍土, 2014,36(1):166-172.

    [24] 陳西平.雅魯藏布江水環(huán)境的探討[J]. 環(huán)境科學(xué), 1989,(5):85-89.

    [25] Diamantini E, Lutz S R, Mallucci S, et al. Driver detection of water quality trends in three large European river basins [J]. Science of the Total Environment, 2018,612:49-62.

    [26] Zayed J, Guessous A, Lambert J, et al. Estimation of annual Mn emissions from MMT source in the Canadian environment and the Mn pollution index in each province [J]. Science of the Total Environment, 2003,312(1-3):147.

    [27] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

    致謝:本研究由西藏大學(xué)及西藏自治區(qū)各級(jí)政府支持下完成,在此表示感謝.

    Spatiotemporal variations of hydrochemistry and its controlling factors in the Yarlung Tsangpo River.

    LIU Jia-ju1, ZHAO Yu-shun2, HUANG Xiang3, GUO Huai-cheng1*

    (1.College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China;2.College of Environmental Sciences and Engineering, Tianjin University, Tianjin 300072, China;3.Department of Chemistry and Environmental Science, Tibet University, Lhasa 850000, China)., 2018,38(11):4289~4297

    In order to study the characteristics of the water chemistry of the Yarlung Tsangpo River and its control factors, 212 water samples from the mainstream of the Yarlung Tsangpo River and its tributaries were collected in the three water periods in 2016. The mathematical statistics, Piper three-line map, Gibbs model, and ion ratio were analyzed. The hydrochemical characteristics of rivers in the whole basin, and the evolution of water chemistry in the basin. The results showed that the cations in the river water were dominated by Ca2+, Na+, and Mg2+; the anions were dominated by HCO3-and SO42-; the average value of TDS was 204.51mg/L, the degree of salinity is relatively low; the type of water chemistry is dominated by HCO3·SO4(SO4·HCO3)-Ca·Mg (Mg·Ca) type water; The main ion concentration changes in the main stream of Yarlung Tsangpo fluctuated, and the seasonal changes were significant. Spatially, most ions in the river follow the trend of increasing first and then decreasing. The chemical samples of water are distributed in the middle left part of the Gibbs model, indicating that the chemical composition of the watershed in the basin is controlled by rock weathering. Principal component analysis and related analysis show that the Yarlung Tsangpo River Basin The chemical composition of water is controlled by human influence, and weathering of calcite, dolomite and dissolution of sulfuric acid also play a very important role. The vast majority of heavy metals in the three watersheds of the basin meet the requirements of Class I water bodies in surface waters.

    hydrochemical;Yarlung Tsangpo River;spatio-temporal variations;river water quality

    X826

    A

    1000-6923(2018)11-4289-09

    劉佳駒(1991-),男,河北保定人,北京大學(xué)博士研究生,研究方向?yàn)榄h(huán)境規(guī)劃.發(fā)表論文3篇.

    2018-04-02

    西藏環(huán)境水化學(xué)調(diào)查及飲用水安全評(píng)價(jià)技術(shù)研究項(xiàng)目(22015FY111000)

    * 責(zé)任作者, 教授, ghc1111@pku.edu.cn

    猜你喜歡
    平水雅魯藏布江豐水期
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    牛欄江流域浮游動(dòng)物群落結(jié)構(gòu)特征及驅(qū)動(dòng)因子
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調(diào)查分析
    紹興市柯橋區(qū)平水鎮(zhèn)中心幼兒園幼兒作品發(fā)表
    五絕·平水 四支韻游王桃園賞桃花(外一首)
    陶山(2019年2期)2019-12-15 20:52:00
    金代平水縣相關(guān)問題考辨
    中國與孟加拉國在雅魯藏布江河流治理中的合作與問題探究
    商(2016年32期)2016-11-24 15:28:37
    西藏雅魯藏布江區(qū)橋梁水文計(jì)算
    西藏科技(2015年2期)2015-09-26 12:18:32
    雅魯藏布江特大橋的臨時(shí)鋼棧橋結(jié)構(gòu)計(jì)算
    西藏科技(2015年3期)2015-09-26 12:11:11
    滇池流域入湖河流水文形貌特征對(duì)豐水期大型底棲動(dòng)物群落結(jié)構(gòu)的影響
    嘟嘟电影网在线观看| 久久久久久久久久人人人人人人| 亚洲精品,欧美精品| 男女边吃奶边做爰视频| 国产老妇伦熟女老妇高清| 成年av动漫网址| 精品少妇黑人巨大在线播放| 日本午夜av视频| 久久久久人妻精品一区果冻| 在现免费观看毛片| 久久久久久久久久久丰满| 亚洲,欧美,日韩| 91午夜精品亚洲一区二区三区| 777米奇影视久久| videossex国产| 在线亚洲精品国产二区图片欧美 | 人妻一区二区av| 欧美xxⅹ黑人| 日韩制服骚丝袜av| 人妻人人澡人人爽人人| 欧美精品高潮呻吟av久久| 欧美激情极品国产一区二区三区 | 男的添女的下面高潮视频| 国产精品国产三级国产av玫瑰| 热re99久久国产66热| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| av黄色大香蕉| 你懂的网址亚洲精品在线观看| 亚洲伊人久久精品综合| 免费观看的影片在线观看| 亚洲图色成人| 亚洲电影在线观看av| 国产成人91sexporn| av卡一久久| 精品99又大又爽又粗少妇毛片| 亚洲欧美日韩东京热| 国产一区二区三区av在线| 99九九线精品视频在线观看视频| 亚洲精品中文字幕在线视频 | 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 乱人伦中国视频| 国产白丝娇喘喷水9色精品| av福利片在线| 最后的刺客免费高清国语| 丝瓜视频免费看黄片| 国产在线免费精品| 99久久中文字幕三级久久日本| 久久午夜福利片| 99久久人妻综合| 亚洲成色77777| 久久久久精品久久久久真实原创| 在线观看免费高清a一片| 国产精品蜜桃在线观看| 久久女婷五月综合色啪小说| 欧美精品国产亚洲| 狂野欧美激情性bbbbbb| 午夜福利,免费看| 国产精品国产三级专区第一集| 多毛熟女@视频| 精品人妻熟女av久视频| 久久久欧美国产精品| 亚洲精品视频女| 丰满乱子伦码专区| 亚洲精品456在线播放app| 久久久久久久久久成人| 日韩精品有码人妻一区| 婷婷色综合www| 在线观看人妻少妇| 久久久久国产网址| 男的添女的下面高潮视频| 久久精品久久久久久噜噜老黄| 国产精品欧美亚洲77777| av国产久精品久网站免费入址| 亚洲精品国产av蜜桃| 国产精品.久久久| 少妇人妻 视频| 在线观看免费高清a一片| 久久精品国产a三级三级三级| 老司机影院毛片| 美女国产视频在线观看| 国产精品一区www在线观看| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 极品教师在线视频| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 一本久久精品| 少妇的逼水好多| 国产精品福利在线免费观看| 午夜91福利影院| 精品人妻一区二区三区麻豆| 美女大奶头黄色视频| 国产亚洲午夜精品一区二区久久| 91精品国产九色| 黄色一级大片看看| 少妇的逼好多水| 久久精品夜色国产| 熟女av电影| 亚洲天堂av无毛| 高清不卡的av网站| 国产一区有黄有色的免费视频| 高清黄色对白视频在线免费看 | 十八禁网站网址无遮挡 | 黄色欧美视频在线观看| 亚洲精品乱码久久久久久按摩| 亚洲av电影在线观看一区二区三区| 久久精品久久久久久久性| 久久久久国产网址| 久久人人爽人人爽人人片va| 国产精品久久久久久精品古装| 看免费成人av毛片| 精品午夜福利在线看| a 毛片基地| 日韩在线高清观看一区二区三区| 国产精品秋霞免费鲁丝片| 日韩免费高清中文字幕av| 婷婷色av中文字幕| av天堂中文字幕网| 亚洲精品国产色婷婷电影| a级毛色黄片| 岛国毛片在线播放| 精品久久久久久久久亚洲| 亚洲高清免费不卡视频| 丁香六月天网| 男女免费视频国产| 啦啦啦中文免费视频观看日本| 乱人伦中国视频| 在线看a的网站| 夫妻午夜视频| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品一二三| 哪个播放器可以免费观看大片| 国产成人免费无遮挡视频| 新久久久久国产一级毛片| 三级国产精品欧美在线观看| 在线观看美女被高潮喷水网站| 纯流量卡能插随身wifi吗| 精品国产一区二区三区久久久樱花| 啦啦啦在线观看免费高清www| 午夜激情福利司机影院| 老司机影院毛片| 一边亲一边摸免费视频| 狠狠精品人妻久久久久久综合| 你懂的网址亚洲精品在线观看| 夫妻午夜视频| 亚洲国产精品国产精品| 国产极品天堂在线| 一本久久精品| 十八禁网站网址无遮挡 | 80岁老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 久久久欧美国产精品| 最黄视频免费看| 久久国内精品自在自线图片| 亚洲久久久国产精品| 少妇人妻久久综合中文| 韩国av在线不卡| 岛国毛片在线播放| 亚洲精品国产色婷婷电影| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 亚洲精品久久午夜乱码| 成年女人在线观看亚洲视频| 少妇人妻 视频| av女优亚洲男人天堂| 99re6热这里在线精品视频| 精品国产一区二区久久| a级毛色黄片| 免费av不卡在线播放| 午夜激情久久久久久久| 国产69精品久久久久777片| 黄色日韩在线| 丰满乱子伦码专区| 日本91视频免费播放| 午夜日本视频在线| 能在线免费看毛片的网站| av专区在线播放| 看非洲黑人一级黄片| 国产精品久久久久成人av| 蜜桃在线观看..| 人人妻人人澡人人看| 男人和女人高潮做爰伦理| 亚洲国产精品999| 亚洲伊人久久精品综合| 国产精品国产三级国产专区5o| 黄色日韩在线| 国产综合精华液| 乱码一卡2卡4卡精品| 婷婷色麻豆天堂久久| 美女福利国产在线| 人妻 亚洲 视频| 国产欧美日韩综合在线一区二区 | 成人亚洲精品一区在线观看| 成人18禁高潮啪啪吃奶动态图 | 国产av一区二区精品久久| 免费观看a级毛片全部| 久久国内精品自在自线图片| 亚洲中文av在线| 国产精品免费大片| 乱码一卡2卡4卡精品| 久久99精品国语久久久| 中文字幕精品免费在线观看视频 | 久久婷婷青草| 桃花免费在线播放| 女的被弄到高潮叫床怎么办| 亚洲精品自拍成人| 国产精品熟女久久久久浪| 国产探花极品一区二区| 精品久久久久久久久av| 狠狠精品人妻久久久久久综合| 欧美精品一区二区大全| 午夜福利影视在线免费观看| 免费看光身美女| 啦啦啦在线观看免费高清www| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 国产免费视频播放在线视频| 国产毛片在线视频| 午夜激情久久久久久久| 久久久久国产精品人妻一区二区| 国产精品一区二区三区四区免费观看| √禁漫天堂资源中文www| 麻豆成人午夜福利视频| 精品酒店卫生间| 国国产精品蜜臀av免费| 一级爰片在线观看| 九九久久精品国产亚洲av麻豆| 亚洲精品成人av观看孕妇| 夫妻性生交免费视频一级片| 国产日韩一区二区三区精品不卡 | 亚洲精华国产精华液的使用体验| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩中文字幕视频在线看片| 乱人伦中国视频| 只有这里有精品99| 久久久久网色| 国产高清不卡午夜福利| 9色porny在线观看| 赤兔流量卡办理| 国产成人免费观看mmmm| 自线自在国产av| 国产极品天堂在线| 国产精品99久久99久久久不卡 | 一区二区av电影网| 日本黄色片子视频| 日韩av在线免费看完整版不卡| 美女内射精品一级片tv| 国产成人freesex在线| 国产男女内射视频| 国产日韩欧美在线精品| av免费在线看不卡| 亚洲国产精品国产精品| 日本色播在线视频| 在线观看www视频免费| 三级国产精品欧美在线观看| 亚洲va在线va天堂va国产| 精品人妻熟女毛片av久久网站| 少妇的逼好多水| 日日啪夜夜爽| 国产伦精品一区二区三区视频9| 国产精品国产三级国产专区5o| 亚洲欧美精品自产自拍| 人妻人人澡人人爽人人| av福利片在线观看| 日日爽夜夜爽网站| 麻豆成人av视频| 国产成人精品婷婷| 国产日韩欧美在线精品| av国产精品久久久久影院| 人人妻人人看人人澡| 爱豆传媒免费全集在线观看| 精品亚洲成a人片在线观看| 日本色播在线视频| av女优亚洲男人天堂| 最近中文字幕2019免费版| 黄色怎么调成土黄色| 久久久精品免费免费高清| 色吧在线观看| 97在线人人人人妻| 国产 一区精品| 久久精品夜色国产| 久久久亚洲精品成人影院| 我要看黄色一级片免费的| 看十八女毛片水多多多| 国产综合精华液| 国产av一区二区精品久久| 桃花免费在线播放| 亚洲不卡免费看| 深夜a级毛片| 亚洲欧美日韩另类电影网站| 久久久久久久久久久久大奶| 久久久久网色| 人妻夜夜爽99麻豆av| 美女福利国产在线| 精品少妇黑人巨大在线播放| 色视频在线一区二区三区| 国产精品国产av在线观看| 爱豆传媒免费全集在线观看| 一个人免费看片子| 亚洲伊人久久精品综合| 亚洲精品一区蜜桃| 亚洲国产av新网站| 亚洲精品一区蜜桃| 国产成人91sexporn| 亚洲精品国产av蜜桃| 国产午夜精品久久久久久一区二区三区| 精品人妻熟女av久视频| 国产永久视频网站| 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 国产免费视频播放在线视频| 成人免费观看视频高清| 久久韩国三级中文字幕| 女性生殖器流出的白浆| av卡一久久| 国产欧美亚洲国产| 久久6这里有精品| 色5月婷婷丁香| 91久久精品国产一区二区三区| 国产日韩欧美在线精品| 七月丁香在线播放| 免费黄频网站在线观看国产| 丰满人妻一区二区三区视频av| 久久久久久久久久久丰满| 少妇丰满av| 男人狂女人下面高潮的视频| 亚洲美女黄色视频免费看| 一级黄片播放器| 日韩一本色道免费dvd| 男人爽女人下面视频在线观看| 全区人妻精品视频| 99久久综合免费| 国产精品.久久久| 亚洲色图综合在线观看| 看免费成人av毛片| 91精品伊人久久大香线蕉| 人妻制服诱惑在线中文字幕| 国精品久久久久久国模美| 国产在线视频一区二区| 一二三四中文在线观看免费高清| 99re6热这里在线精品视频| 丰满迷人的少妇在线观看| 国产精品国产av在线观看| √禁漫天堂资源中文www| 免费人成在线观看视频色| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产精品一区三区| 伦精品一区二区三区| 欧美性感艳星| 亚洲在久久综合| 少妇裸体淫交视频免费看高清| 亚洲国产精品999| 青春草亚洲视频在线观看| 97超碰精品成人国产| 伦理电影免费视频| 国产精品人妻久久久久久| 久久久精品94久久精品| 中国美白少妇内射xxxbb| 成人影院久久| 免费看不卡的av| 久久ye,这里只有精品| 国产成人午夜福利电影在线观看| 99九九线精品视频在线观看视频| 十八禁网站网址无遮挡 | 免费观看av网站的网址| 亚洲av免费高清在线观看| 精品熟女少妇av免费看| 曰老女人黄片| 男人舔奶头视频| 99热全是精品| 如日韩欧美国产精品一区二区三区 | 丰满人妻一区二区三区视频av| 丁香六月天网| 五月玫瑰六月丁香| 日韩欧美一区视频在线观看 | 久久久久国产网址| 18禁在线无遮挡免费观看视频| 色94色欧美一区二区| 美女内射精品一级片tv| 亚洲电影在线观看av| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 国产 精品1| 男女啪啪激烈高潮av片| 妹子高潮喷水视频| 日韩精品有码人妻一区| 中文在线观看免费www的网站| 99热网站在线观看| 精品久久久噜噜| 大片电影免费在线观看免费| 亚洲人成网站在线观看播放| 欧美日韩视频高清一区二区三区二| 精品人妻熟女毛片av久久网站| 国产高清有码在线观看视频| 少妇高潮的动态图| 一本久久精品| 欧美精品人与动牲交sv欧美| 久久99热6这里只有精品| 国产一区二区三区av在线| 超碰97精品在线观看| 欧美精品一区二区免费开放| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆| videos熟女内射| 成人国产av品久久久| 免费大片黄手机在线观看| 中文字幕av电影在线播放| 女人久久www免费人成看片| av卡一久久| 高清午夜精品一区二区三区| 久久久久精品性色| 国产黄频视频在线观看| 中文字幕亚洲精品专区| 我的女老师完整版在线观看| 日韩精品免费视频一区二区三区 | 久久久午夜欧美精品| 啦啦啦啦在线视频资源| 日韩精品有码人妻一区| 又黄又爽又刺激的免费视频.| 日韩大片免费观看网站| 秋霞伦理黄片| 成人漫画全彩无遮挡| 日本黄色日本黄色录像| av国产久精品久网站免费入址| 久久人妻熟女aⅴ| 国产 精品1| 亚洲精品,欧美精品| 女的被弄到高潮叫床怎么办| 女人精品久久久久毛片| 久热这里只有精品99| 一级二级三级毛片免费看| 日本黄色日本黄色录像| 最后的刺客免费高清国语| 国产成人精品无人区| 在线免费观看不下载黄p国产| av播播在线观看一区| 三上悠亚av全集在线观看 | 少妇的逼水好多| 欧美日韩视频高清一区二区三区二| 精品酒店卫生间| 午夜福利影视在线免费观看| 91午夜精品亚洲一区二区三区| 性色av一级| 精品国产乱码久久久久久小说| 高清在线视频一区二区三区| 在线观看三级黄色| 亚洲无线观看免费| 国产精品久久久久久av不卡| 久久99一区二区三区| 久久精品国产亚洲网站| 99热国产这里只有精品6| 日本91视频免费播放| 国产精品三级大全| 国产探花极品一区二区| 亚洲综合精品二区| a级毛片在线看网站| 男人狂女人下面高潮的视频| 爱豆传媒免费全集在线观看| 秋霞伦理黄片| 精品少妇久久久久久888优播| 波野结衣二区三区在线| 女性被躁到高潮视频| 日产精品乱码卡一卡2卡三| 男人狂女人下面高潮的视频| 我的女老师完整版在线观看| 欧美变态另类bdsm刘玥| 国精品久久久久久国模美| 国产有黄有色有爽视频| 久久人人爽人人片av| 久久精品国产自在天天线| av在线老鸭窝| 亚洲av电影在线观看一区二区三区| 黄色一级大片看看| 五月天丁香电影| 免费人妻精品一区二区三区视频| 99久久人妻综合| 2018国产大陆天天弄谢| 不卡视频在线观看欧美| 春色校园在线视频观看| 免费高清在线观看视频在线观看| 久久精品国产亚洲av天美| 久久人人爽人人爽人人片va| 内射极品少妇av片p| 久久午夜福利片| 亚洲av欧美aⅴ国产| 国产精品伦人一区二区| 国产成人免费观看mmmm| 免费看不卡的av| 卡戴珊不雅视频在线播放| 亚洲欧美一区二区三区黑人 | 如日韩欧美国产精品一区二区三区 | 国产日韩欧美视频二区| 一区二区av电影网| 国产伦精品一区二区三区四那| 天堂8中文在线网| 久久国产精品大桥未久av | 哪个播放器可以免费观看大片| 国产日韩欧美亚洲二区| 噜噜噜噜噜久久久久久91| 国产黄片美女视频| 国产精品99久久99久久久不卡 | 中文字幕av电影在线播放| 欧美xxⅹ黑人| 亚洲成人一二三区av| 七月丁香在线播放| 欧美老熟妇乱子伦牲交| 日韩欧美精品免费久久| 日日啪夜夜撸| 亚洲精品一区蜜桃| 亚洲精品国产成人久久av| 久久久久久久久久久免费av| 乱系列少妇在线播放| 美女主播在线视频| 青春草视频在线免费观看| 天堂8中文在线网| 夜夜看夜夜爽夜夜摸| 久久99热6这里只有精品| 亚洲av男天堂| 日韩在线高清观看一区二区三区| 成人无遮挡网站| 男人舔奶头视频| av线在线观看网站| 国产视频首页在线观看| 久久久久久久久久人人人人人人| 一级黄片播放器| 亚洲国产精品专区欧美| 亚洲精品一二三| 丰满少妇做爰视频| 自拍偷自拍亚洲精品老妇| 在线观看人妻少妇| 国产综合精华液| 日本-黄色视频高清免费观看| 国产亚洲一区二区精品| 各种免费的搞黄视频| 久久久久久久久久成人| 国产精品熟女久久久久浪| 国产 一区精品| 国产色婷婷99| 亚洲欧美日韩卡通动漫| 美女大奶头黄色视频| 内地一区二区视频在线| 欧美另类一区| 成人亚洲欧美一区二区av| 亚洲精品成人av观看孕妇| 亚洲综合色惰| 色5月婷婷丁香| 久久国内精品自在自线图片| 日本色播在线视频| 精品视频人人做人人爽| 男男h啪啪无遮挡| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 一级毛片久久久久久久久女| 成年美女黄网站色视频大全免费 | 亚洲国产日韩一区二区| 丝袜脚勾引网站| 日本wwww免费看| 在线播放无遮挡| 我要看黄色一级片免费的| 亚洲天堂av无毛| 久久精品熟女亚洲av麻豆精品| 久久青草综合色| 纵有疾风起免费观看全集完整版| 中文字幕制服av| 欧美国产精品一级二级三级 | 最近2019中文字幕mv第一页| 少妇人妻久久综合中文| 国产伦在线观看视频一区| 国产精品嫩草影院av在线观看| 看非洲黑人一级黄片| 欧美变态另类bdsm刘玥| 国产在视频线精品| 69精品国产乱码久久久| 国产在线男女| 午夜老司机福利剧场| 国产精品无大码| 久久久久久久亚洲中文字幕| 久久97久久精品| av网站免费在线观看视频| 久久精品国产鲁丝片午夜精品| 黑人猛操日本美女一级片| 精品国产一区二区三区久久久樱花| 国产69精品久久久久777片| 又粗又硬又长又爽又黄的视频| 成人免费观看视频高清| 国产乱来视频区| 亚洲国产欧美在线一区| www.av在线官网国产| 免费黄色在线免费观看| 一本色道久久久久久精品综合| 人妻制服诱惑在线中文字幕| 最近的中文字幕免费完整| 丰满人妻一区二区三区视频av| 日韩,欧美,国产一区二区三区| 国产黄色视频一区二区在线观看| 女人精品久久久久毛片| 亚洲欧美日韩东京热| 少妇 在线观看| 日韩中文字幕视频在线看片| 五月开心婷婷网| 国产精品国产三级国产专区5o| 在线观看三级黄色| 亚洲av二区三区四区| 狂野欧美白嫩少妇大欣赏| 97超视频在线观看视频| 少妇高潮的动态图| 青春草国产在线视频| 一级片'在线观看视频| 亚洲成人av在线免费|