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

    內(nèi)蒙古高原3大典型湖泊水化學(xué)特征及其控制因素分析

    2015-03-12 07:41:40吳用史小紅趙勝男林濤濤馬軍

    吳用,史小紅*,趙勝男,林濤濤,馬軍

    1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古 呼和浩特 010018;2. 長(zhǎng)江水利委員會(huì)長(zhǎng)江三峽水文水資源勘測(cè)局,湖北 宜昌 443000;3. 河套灌區(qū)排水事業(yè)管理局總排干溝排水站,內(nèi)蒙古 巴彥淖爾 014400

    內(nèi)蒙古高原3大典型湖泊水化學(xué)特征及其控制因素分析

    吳用1,史小紅1*,趙勝男1,林濤濤2,馬軍3

    1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古 呼和浩特 010018;2. 長(zhǎng)江水利委員會(huì)長(zhǎng)江三峽水文水資源勘測(cè)局,湖北 宜昌 443000;3. 河套灌區(qū)排水事業(yè)管理局總排干溝排水站,內(nèi)蒙古 巴彥淖爾 014400

    摘要:內(nèi)蒙古高原地處我國(guó)北方半干旱與干旱氣候區(qū),經(jīng)度跨度大,由東向西氣候呈現(xiàn)半干旱區(qū)向干旱區(qū)過渡特征。區(qū)內(nèi)湖泊多屬內(nèi)陸型湖泊,湖泊水化學(xué)特征受地理位置和氣候環(huán)境影響顯著。文章在空間上從東向西依次選擇呼倫湖、達(dá)里湖和烏梁素海為研究對(duì)象,通過對(duì)不同氣候區(qū)域和不同類型的湖泊水體中主要離子組成進(jìn)行分析對(duì)比研究,揭示區(qū)域內(nèi)湖泊水化學(xué)特征及其控制因素。結(jié)果表明,內(nèi)蒙古高原區(qū)內(nèi)湖泊水化學(xué)類型單一,陽離子以Na+主導(dǎo),陰離子以Cl-主導(dǎo),但各湖泊水體中主要離子質(zhì)量濃度差異較大。呼倫湖、達(dá)里湖和烏梁素海3個(gè)湖泊Na+分別占其陽離子摩爾總數(shù)的76.6%、97.2%和66.0%;Cl-分別占其陰離子摩爾總數(shù)的49.6%、55.5%和45.2%。湖泊水體中Na+和SO42-、Cl-的變異系數(shù)相對(duì)較小,在水體中質(zhì)量濃度相對(duì)穩(wěn)定。借助Piper三角圖,分析出呼倫湖和烏梁素海陰離子分布在SO42--Cl-線上,陽離子分布在Na++K+線上;達(dá)里湖的陰離子分布在HCO3--Cl-線上,陽離子分布在Na++K+線上,說明區(qū)內(nèi)湖泊受蒸發(fā)結(jié)晶和巖石風(fēng)化共同作用影響。湖泊水體中主要離子質(zhì)量濃度的空間變化特征受徑流,補(bǔ)排關(guān)系的影響,結(jié)合Gibbs模型圖進(jìn)一步說明了3個(gè)湖泊受蒸發(fā)濃縮作用和人類活動(dòng)影響。湖泊水化學(xué)特征的研究,對(duì)于保護(hù)區(qū)域內(nèi)湖泊環(huán)境具有重要意義。

    關(guān)鍵詞:內(nèi)蒙古高原湖泊;干旱與半干旱氣候區(qū);主要離子組成;水化學(xué)特征

    引用格式:吳用,史小紅,趙勝男,林濤濤,馬軍. 內(nèi)蒙古高原3大典型湖泊水化學(xué)特征及其控制因素分析[J]. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(7): 1202-1208.

    WU Yong, SHI Xiaohong, ZHAO Shengnan, LIN Taotao, MA Jun. Major Ion Chemistry and Influencing Factors of Three Typical Lakes in Inner Mongolia Plateau [J]. Ecology and Environmental Sciences, 2015, 24(7): 1202-1208.

    湖泊水體中的化學(xué)離子被視為天然的“示蹤劑”,能夠反映湖泊水循環(huán)過程,流域水-巖作用以及湖泊演化歷史(周嘉欣等,2014)。研究水體中主要離子的組成特征,可用于識(shí)別水體化學(xué)組成的基本過程。國(guó)內(nèi)外已有不少學(xué)者通過對(duì)各地主要湖泊的水化學(xué)組成和特征進(jìn)行了分析研究,研究結(jié)果表明湖泊水化學(xué)組成受到大氣沉降、巖石風(fēng)化、水體生物和人類活動(dòng)的綜合影響(王蘇民等,1995)。

    我國(guó)湖泊水化學(xué)特征的研究大多集中在長(zhǎng)江中下游湖泊群,西南巖溶湖、青藏高原湖區(qū)和新疆湖區(qū),研究成果較多(胡春華等,2011;白占國(guó)等,1998;侯韶華等,2009)。處于我國(guó)西北干旱和半干旱的湖區(qū),生態(tài)環(huán)境脆弱,對(duì)全球氣候變化敏感,區(qū)域內(nèi)的湖泊具有較高的生態(tài)價(jià)值和經(jīng)濟(jì)價(jià)值,對(duì)區(qū)域內(nèi)湖泊水化學(xué)特征和鹽化過程的深入研究非常重要。近年來,學(xué)者們對(duì)干旱區(qū)和半干旱區(qū)的高原湖泊進(jìn)行了相關(guān)研究,如王鵬等(2013)研究了青藏高原打加芒錯(cuò)湖水的水化學(xué)組成特征及其鹽化過程,王海雷等(2010)對(duì)高原上400多個(gè)湖泊的水化學(xué)成分與湖水鹽度的相關(guān)性進(jìn)行了綜合分析,鄭綿平等(2010)分析了青藏高原鹽湖水化學(xué)類型及其地帶分布,指出不同的水化學(xué)類型具有不同的礦物組合特征。內(nèi)蒙古高原的湖泊研究中,甄志磊(2014)利用同位素研究了達(dá)里湖與地下水的補(bǔ)給關(guān)系,汪敬忠等(2013)結(jié)合氫、氧同位素和水化學(xué)特征研究了河套平原地下水對(duì)烏梁素海礦化度的影響;曾海鰲等(2010)人分析了蒙新高原湖區(qū)各湖泊水體礦化度的影響因素。針對(duì)內(nèi)蒙古地區(qū)的湖泊,其空間上氣候條件差異顯著,湖泊類型和地理環(huán)境復(fù)雜,需要研究區(qū)內(nèi)典型湖泊水化學(xué)特征并進(jìn)行對(duì)比分析,為進(jìn)一步了解內(nèi)蒙古地區(qū)各個(gè)湖泊的水環(huán)境特征提供依據(jù)。

    內(nèi)蒙古高原地處內(nèi)陸,區(qū)內(nèi)湖泊大多屬內(nèi)陸湖泊,氣候干旱少雨,地表徑流補(bǔ)給不豐,湖泊逐漸萎縮咸化。加之,高原橫跨經(jīng)度大,從東向西氣候由半干旱向干旱過渡,蒸發(fā)量由1500 mm逐漸增大2400 mm,降雨量由256 mm逐漸減小至140 mm,蒸發(fā)量約是降雨量的5~10倍。由東向西選擇區(qū)內(nèi)3個(gè)典型湖泊:呼倫湖、達(dá)里湖和烏梁素海,湖泊類型分別為草原湖泊,巖溶堵塞湖和灌區(qū)湖泊。文中從水化學(xué)的角度出發(fā),系統(tǒng)分析不同類型湖泊(呼倫湖、達(dá)里湖和烏梁素海)主離子水化學(xué)特征及其分布規(guī)律,并試圖對(duì)各湖泊水化學(xué)成因及演化過程進(jìn)行探討,以期為內(nèi)蒙古高原干旱與半干旱湖泊的生態(tài)環(huán)境治理提供基礎(chǔ)資料和科學(xué)依據(jù)。

    1 研究區(qū)概況

    呼倫湖也稱達(dá)賚湖,地理坐標(biāo)為:116°58′~117°48′E,48°33′~49°20′N,坐落于呼倫貝爾大草原中部,是典型的草原湖泊,也是內(nèi)蒙古最大的湖泊,蓄水量為138×108m3,湖水面積2339 km2。呼倫湖主要補(bǔ)給水源為克魯倫河、烏爾遜河,達(dá)蘭鄂羅木河(新開河)。其中新開河為吞吐性河流,海拉爾河水大時(shí),順此河流入呼倫湖,呼倫湖水大時(shí)又順此河流向額爾古納河(李暢游等,2011;王志杰等,2012)。達(dá)里湖位于赤峰市克什克騰旗西北部,其地理坐標(biāo)為:116°22′~117°00′E,43°11′~43°27′N,是一個(gè)巖溶堵塞湖,蓄水量為15.53×108m3,湖水面積224.6 km2。湖水主要依賴湖面降水和地表地下徑流補(bǔ)給。入湖河流主要有貢格爾河、沙里河、亮子河及浩來河,沒有排泄出口(甄志磊等,2014)。烏梁素海位于烏拉特前旗境內(nèi),地理坐標(biāo)為:108°43′~108°57′E,40°47′~41°03′N,是黃河改道后形成的牛軛湖和典型的灌區(qū)湖泊,是河套灌區(qū)排灌水系的重要組成部分,蓄水量為2.64×108m3,湖泊現(xiàn)有水域面積285.38 km2(汪敬忠等,2013;付緒金等,2013)。內(nèi)蒙古高原的內(nèi)陸湖泊,蒸發(fā)強(qiáng)度較大以致超過湖水的補(bǔ)給量,湖泊水體因不斷被濃縮而發(fā)育成閉流類的咸水湖和鹽湖。本研究選取的呼倫湖、達(dá)里湖和烏梁素海,是內(nèi)蒙古高原主要和具有代表性的湖泊,內(nèi)蒙古高原湖區(qū)如圖1。

    2 樣品采集與分析方法

    2012年6月至8月,內(nèi)蒙古高原湖泊都處于非冰封期。對(duì)3個(gè)湖泊的水質(zhì)狀況進(jìn)行依次調(diào)查采樣。其中:呼倫湖采樣點(diǎn)9個(gè),達(dá)里湖采樣點(diǎn)9個(gè),烏梁素海采樣點(diǎn)11個(gè),采樣深度均為水面以下0.6 m處。采樣點(diǎn)編號(hào)和位置如圖2?,F(xiàn)場(chǎng)采用梅特勒(SG8-ELK,SG9-ELK)多參數(shù)水質(zhì)分析儀測(cè)定湖水礦化度,鹽度、pH值和電導(dǎo)率;HCO3-質(zhì)量濃度采用酸堿指示劑滴定法(國(guó)家環(huán)境保護(hù)總局,2002);陰、陽離子測(cè)定前,水樣經(jīng)0.45 μm的玻璃纖維濾膜過濾后,使用潔凈聚氟乙烯瓶封裝后取回水樣,在實(shí)驗(yàn)室采用離子色譜法進(jìn)行分析測(cè)定,所采用儀器為Dionex ICS-1100型號(hào)的色譜儀。

    圖1 內(nèi)蒙古高原湖泊區(qū)域概況圖Fig. 1 The locations of the lakes in the region of Inner Mongolia plateau

    圖2 研究區(qū)采樣點(diǎn)Fig. 2 The samples in study areas

    3 結(jié)果與討論

    3.1湖泊水化學(xué)指標(biāo)統(tǒng)計(jì)

    總?cè)芙庑怨腆w(TDS)是水化學(xué)成分的重要指標(biāo),一般以溶解于水中的主要離子總量表示,受降水量、蒸發(fā)量、流域巖石成分等因素影響。

    3個(gè)典型湖泊水化學(xué)指標(biāo)統(tǒng)計(jì)結(jié)果(表1)表明,呼倫湖的TDS最小,為1520 mg·L-1,達(dá)里湖最大為5990 mg·L-1,烏梁素海居中為2250 mg·L-1。分析原因,呼倫湖是內(nèi)蒙古東部草原型湖泊,屬中溫帶半干旱偏濕地區(qū),多年平均年降水量為319 mm,年蒸發(fā)量為1400~1900 mm,劇烈的蒸發(fā)導(dǎo)致呼倫湖水體的TDS值高于濕潤(rùn)地區(qū)湖泊,如太湖、鄱陽湖和洞庭湖,其TDS值分別為345.86、99.04 和165.48 mg·L-1(蔡永久等,2013)。但呼倫湖TDS值要低于干旱區(qū)和半干旱區(qū)的達(dá)里湖和烏梁素海,這是由于克魯倫河和烏爾遜河兩大河的補(bǔ)給,年補(bǔ)給量達(dá)11.94×108m3,因而湖泊水體的TDS值偏低。達(dá)里湖是內(nèi)蒙古高原干旱地區(qū)巖溶堵塞湖,湖泊水沒有排泄出口,湖泊水的排泄主要依靠劇烈的蒸發(fā)。河流帶到湖中的鹽分不能排出,鹽分在湖中不斷累積,導(dǎo)致湖水TDS值偏高。烏梁素海地處內(nèi)蒙古中西部,其屬于黃河改道形成的牛軛湖,水量的補(bǔ)給來源主要靠黃河水和河套灌區(qū)農(nóng)田退水,與河套灌區(qū)的灌溉措施密切相關(guān),湖泊換水周期較快,湖泊水的水化學(xué)類型和水質(zhì)取決黃河水和農(nóng)田退水。加之同時(shí)受干旱區(qū)氣候和灌區(qū)內(nèi)人類活動(dòng),企業(yè)、工廠的污廢水的影響,湖泊水體的TDS值相對(duì)較高。

    表1 3個(gè)湖泊主要離子質(zhì)量濃度對(duì)比分析Table 1 Analysis the content of main ions is of their discrepancy in each of lakes mg·L-1

    3個(gè)湖泊同屬內(nèi)蒙古高原干旱區(qū)和半干旱區(qū),受氣候條件的影響,區(qū)域內(nèi)蒸發(fā)量遠(yuǎn)大于降水量,湖泊湖水的排泄很大一部分靠蒸發(fā),因此,3大湖泊都有向微咸水和咸水湖泊開始轉(zhuǎn)化。

    3.2湖泊水化學(xué)特征

    3.2.1湖泊水體水化學(xué)類型

    呼倫湖TDS平均值為1520 mg·L-1,屬于微咸水湖泊。陽離子當(dāng)量濃度大小依次為Na+>Mg2+> K+>Ca2+,湖水優(yōu)勢(shì)陽離子為Na+,占陽離子總摩爾數(shù)的76.6%;陰離子當(dāng)量濃度大小依次為Cl->SO42-> HCO3-,水中優(yōu)勢(shì)陰離子為Cl-、SO42-,分別占陰離子總摩爾數(shù)的49.6%、29.3%。湖泊水的水化學(xué)類型為Cl-·SO42--Na+。

    達(dá)里湖TDS平均值為5990 mg·L-1,屬于咸水湖泊。陽離子當(dāng)量濃度大小依次為Na+>Mg2+>K+> Ca2+,湖水優(yōu)勢(shì)陽離子為Na+,占陽離子總摩爾數(shù)的97.5%;陰離子當(dāng)量濃度大小依次為Cl->HCO3-> SO42-,水中優(yōu)勢(shì)陰離子為Cl-、HCO3-,分別占陰離子總摩爾數(shù)的55.5%、35.3%。湖泊水的水化學(xué)類型為·Cl-·HCO3--Na+。

    烏梁素海TDS平均值為2250 mg·L-1,屬于微咸水湖泊。陽離子當(dāng)量濃度大小依次為Na+>Mg2+> Ca2+>K+,湖水優(yōu)勢(shì)陽離子為Na+,占陽離子總摩爾數(shù)的66.0%;陰離子當(dāng)量濃度大小依次為Cl-> HCO3->SO42-,水中優(yōu)勢(shì)陰離子為Cl-、HCO3-,分別占陰離子總摩爾數(shù)的45.2%、37.5%。湖泊水的水化學(xué)類型為Cl-·HCO3--Na+。

    3個(gè)湖泊的水化學(xué)類型表明,其水化學(xué)組成相對(duì)單一,陽離子為鈉型,陰離子為氯化物型。是處于內(nèi)蒙古高原干旱半干旱區(qū)湖泊的共同點(diǎn),因氣候環(huán)境和水汽來源相近,使得湖泊水水中離子化學(xué)組成和水化學(xué)類型存在一定的相似性。然而水體中各離子摩爾濃度差別較大,主要受水文條件和地質(zhì)環(huán)境的影響。呼倫湖湖盆地層主要是中生界地層及上第三系、第四系地層,巖性基本為沙巖、泥巖、礫巖夾煤系地層。近年來人類活動(dòng),湖區(qū)北面扎賚諾爾礦區(qū)的煤礦工業(yè)得到大力發(fā)展,工業(yè)污廢水的大量排放,夾雜著大量的硫化物進(jìn)入呼倫湖流域,從而導(dǎo)致湖水中SO42-是除Cl-外質(zhì)量濃度最大的陰離子(王蘇民等,1995)。達(dá)里湖存在烏日根塔拉-西拉沐倫東西向大斷裂上,是由玄武巖的噴發(fā)構(gòu)成的堰塞湖,沒有排泄出口。周邊地層中分布著大量的巖漿巖和變質(zhì)巖,每年達(dá)里湖的河川徑流量補(bǔ)給量為0.44×108m3,蒸發(fā)量為2.48×108m3,其虧損水量依靠地下水的補(bǔ)給,水流對(duì)巖層的淋濾作用導(dǎo)致湖泊水體中碳酸鹽的濃度要比同樣是封閉湖泊的呼倫湖要高,湖水中HCO3-離子質(zhì)量濃度僅次于Cl-離子(張郝哲等,2012)。烏梁素海是黃河改道形成的河跡湖,水源主要靠黃河灌溉水和農(nóng)田退水補(bǔ)給,依靠黃河水的大量補(bǔ)給,水體交換條件較好,水體地球化學(xué)環(huán)境以還原性為主,因此湖泊水體中HCO3-質(zhì)量濃度也相對(duì)較高(崔鳳麗等,2013)。3個(gè)湖泊水化學(xué)組成成分的成因各不相同,充分說明內(nèi)蒙古高原內(nèi)不同湖泊類型其水化學(xué)特征有差異,需對(duì)其各自水化學(xué)特征進(jìn)行詳細(xì)的分析。

    3.2.2不同湖泊水體主離子組成變化

    采用通用的水化學(xué)分類方法,將毫克當(dāng)量百分?jǐn)?shù)大于25%的離子參與水化學(xué)類型的分類。利用Gw-chart軟件(USGS開發(fā))繪出陽離子Ca2+-Mg2+-(Na++K+)組成的三角圖和陰離子SO42--Cl--(CO32-+HCO3-)組成的三角圖來表明不同水體的化學(xué)組成特征。湖泊演化一般趨勢(shì)是從淡水湖到咸水湖再到鹽湖,直至干涸結(jié)束演化,其水化學(xué)類型也遵循一定的鹽化規(guī)律:即一般為碳酸鹽型-硫酸鈉型-硫酸鎂型-氯化物型(陸瑩等,2010)。

    陰陽離子三角形圖可表示水體中溶質(zhì)載荷主要離子的相對(duì)豐度和分布特征,從而揭示不同巖石風(fēng)化對(duì)水中總?cè)苜|(zhì)成分的相對(duì)貢獻(xiàn)率。在陰離子HCO3--SO42--Cl-組成的三角組圖中,在碳酸鹽巖地區(qū)的水以HCO3-為主導(dǎo),組分點(diǎn)落在HCO3-一端;在蒸發(fā)鹽巖地區(qū)的水SO42-、Cl-含量較高,其組分點(diǎn)會(huì)落在SO42--Cl-線上,遠(yuǎn)離HCO3-一端。在陽離子Ca2+-Mg2+-(Na++K+)組成的三角圖中也有類似的反映,主要受碳酸鹽巖影響的水體,其組分點(diǎn)會(huì)落在Ca2+-Mg2+線上,靠近Ca2+端元;主要受蒸發(fā)鹽巖影響的水體,其組分點(diǎn)偏向(Na++K+)端元(王瑞久,1983)。

    圖3為呼倫湖、達(dá)里湖和烏梁素海的主要陰陽離子組成三角圖,從中可以看出呼倫湖和烏梁素海的陰離子分布在三角圖SO42--Cl-線上,陽離子分布在三角圖的Na++K+線上;達(dá)里湖的陰離子分布在三角圖HCO3--Cl-線上,陽離子分布在三角圖的Na++K+線上。通過計(jì)算各離子的變異系數(shù),表2得知:呼倫湖和烏梁素海水體中主要陰陽離子Na+和SO42-、Cl-的變異系數(shù)相對(duì)較小,達(dá)里湖水體中主要陰陽離子Na+和Cl-的變異系數(shù)相對(duì)較小,進(jìn)一步說明他們?cè)谒w中的質(zhì)量濃度相對(duì)穩(wěn)定。可以推斷呼倫湖和烏梁素海主要受蒸發(fā)鹽巖控制,達(dá)里湖受碳酸鹽巖和蒸發(fā)鹽巖的共同作用,湖泊離子組成與氣候條件和湖泊地質(zhì)環(huán)境密不可分(葉宏萌等,2010)。

    圖3 3個(gè)湖泊主離子組成Piper圖Fig. 3 The Piper triangle chart for main ions of three lakes

    表3 3個(gè)湖泊主離子現(xiàn)狀與歷史資料對(duì)比分析Table 3 The status quo icons concentration percent compare with history in the three lakes %

    表2 呼倫湖、達(dá)里湖和烏梁素海湖泊水化學(xué)特征值*Table 2 The water chemical characteristics values in Hulun Lake, Dali Lake and Ulansuhai Lake mg·L-1

    對(duì)比20世紀(jì)80~90年代資料(表3),呼倫湖各離子質(zhì)量濃度變化劇烈。近幾十年來,由于流域降雨的減少,河川徑流補(bǔ)給持續(xù)減少,湖泊水體受蒸發(fā)主控,水體中Cl-的增加。煤礦工業(yè)的大力發(fā)展,對(duì)湖泊環(huán)境造成了一定的影響,從而導(dǎo)致水體中SO42-的增多。達(dá)里湖是巖溶堵塞湖,水環(huán)境條件較為封閉,湖泊水體中各離子質(zhì)量濃度變化不大,受蒸發(fā)鹽巖控制,湖水中Na+、K+、和Cl-質(zhì)量濃度變大。取決達(dá)里湖湖泊水文循環(huán)過程,地下徑流大量補(bǔ)給湖水,達(dá)里湖碳酸鹽巖的作用要大于蒸發(fā)鹽巖。烏梁素海雖然依賴黃河水的補(bǔ)給,但同樣受蒸發(fā)鹽巖影響,流域內(nèi)蒸發(fā)量是降水量的10倍左右,水體中Na+、K+增加。

    Piper圖上可以看出,各個(gè)湖泊的離子質(zhì)量濃度在空間上也存在差異。通過繪制呼倫湖、達(dá)里湖和烏梁素??臻g上各離子質(zhì)量濃度變化圖(圖4),并結(jié)合各湖泊空間采樣點(diǎn)的位置(圖2)分析,呼倫湖從北向南湖水中各離子質(zhì)量濃度呈現(xiàn)先增大后減小的過程。在湖泊南部有兩條河流進(jìn)行補(bǔ)給,水動(dòng)力條件較好,水體交換強(qiáng)烈,各離子質(zhì)量濃度較小;湖北部是新開河,為吐納型河流,補(bǔ)給情況下各離子質(zhì)量濃度會(huì)被稀釋。達(dá)里湖空間上,離子變化幅度較小,主要原因是達(dá)里湖沒有外泄出口,湖水交換主要依靠地下水的補(bǔ)給和蒸發(fā)排泄。其全湖整體上由北向南出現(xiàn)降低趨勢(shì)。湖中部和南部有3條河流,水體交換能力較北部活躍,河流的淡水稀釋了入湖口周圍的離子質(zhì)量濃度。烏梁素海的主要的補(bǔ)給口是總排干,湖泊南部只有唯一的泄水渠烏毛計(jì)。湖泊水體從北向南各離子質(zhì)量濃度整體增大,由于南部水動(dòng)力條件較北部偏弱,沒有大量的水流交換,加上強(qiáng)烈的蒸發(fā),流進(jìn)湖泊的鹽分在湖泊南部逐漸積累,造成南部湖泊水體例子質(zhì)量濃度增大。同時(shí)進(jìn)一步說明蒸發(fā)依然是烏梁素海湖泊水體的主要排泄方式之一。

    綜上,不同湖泊水體的主要陰陽離子組成及水化學(xué)類型存在顯著差別,離子組成與水化學(xué)類型呈顯著的區(qū)域性。與歷史資料的對(duì)比分析,現(xiàn)狀下湖泊水體的離子組成的變化主要受氣候條件和人類活動(dòng)的綜合影響。同一湖泊水體中主要陰陽離子組成及水化學(xué)類型存在空間上的差異,受湖泊本身水動(dòng)力條件和補(bǔ)給排泄的綜合影響。

    圖4 3個(gè)湖泊各離子質(zhì)量濃度空間變化Fig. 4 The space changes of content of main ions for three lakes

    3.3不同湖泊水體離子來源及控制因素

    陸地水溶解鹽的可能來源包括大氣攜帶的海鹽(循環(huán)鹽)成分、可溶性巖石(蒸發(fā)巖、硅酸鹽、碳酸鹽、硫化物)的風(fēng)化物和人類活動(dòng)產(chǎn)生的污染物(Berner et al.,1987)。Gibbs通過對(duì)世界雨水、河水、湖泊等地表水化學(xué)組分分析,利用ρ(TDS) 與ρ(Na+)/[ρ(Na+)+ρ(Ca2+)]圖和ρ(TDS)與ρ(Cl-)/ [ρ(Cl-)+ρ(HCO3-)]圖,直觀地反映出水體主要組分的類型:降水控制類型、巖石風(fēng)化類型和蒸發(fā)-濃縮類型(Gibbs,1970;朱海勇等,2013)。該圖把地表水化學(xué)組分劃分為3個(gè)控制端元:巖石風(fēng)化主控端元,大氣降水主控端元和蒸發(fā)-結(jié)晶主控端元。

    將3個(gè)湖泊的水化學(xué)數(shù)據(jù)繪于Gibbs圖中,從圖5看出3個(gè)湖泊的水化學(xué)主控端元落在蒸發(fā)-結(jié)晶端元,TDS大于1000 mg·L-1,ρ(Na+)/[ρ(Na+)+ρ(Ca2+)]比值接近于1,ρ(Cl-)/[ρ(Cl-)+ρ(HCO3-)]比值在0.2~0.9之間,說明湖泊水體分布在蒸發(fā)作用很強(qiáng)的干旱區(qū)域,符合內(nèi)蒙古高原湖泊的特征,受干旱與半干旱氣候的影響。同時(shí)通過計(jì)算呼倫湖、達(dá)里湖和烏梁素海的Cl/Na值,分別為0.54,0.94和1.11,均都小于世界平均海水比值(Cl/Na=1.15)(葉宏萌等,2010),符合內(nèi)陸湖泊的特點(diǎn),受大氣環(huán)流影響較小。但3個(gè)湖泊受干旱和半干旱氣候影響的程度不同,與上文離子組成特征相符。由圖5還可知,3個(gè)湖泊部分區(qū)域還均受到不同程度的人為活動(dòng)的干擾,部分采樣點(diǎn)落在了Gibbs模型外。

    圖5 3個(gè)湖泊Gibbs模型圖示Fig. 5 The Gibbs model for three lakes

    結(jié)合表4中和圖5,分別分析3個(gè)湖泊落在的Gibbs模型外方的采樣點(diǎn),通過對(duì)比發(fā)現(xiàn),呼倫湖D11和G2兩點(diǎn),分別離達(dá)賚湖漁場(chǎng)雙山子分場(chǎng)和大沙圈漁場(chǎng)較近,兩個(gè)漁場(chǎng)人為活動(dòng)在呼倫湖流域相對(duì)較大,對(duì)湖泊水體的影響也就相對(duì)突出;達(dá)里湖的A7和D6兩個(gè)點(diǎn),分別與達(dá)里湖自然保護(hù)區(qū)的達(dá)來諾爾漁場(chǎng)、達(dá)里湖度假中心較近,人為活動(dòng)產(chǎn)生的生活污廢水對(duì)湖泊水體的離子組成造成了一定的影響;烏梁素海落在Gibbs圖外的點(diǎn)主要分布在總排干入湖口處J13、K12、L11點(diǎn)和湖的東南部烏梁素海漁場(chǎng)附近P11、Q10點(diǎn),也說明灌區(qū)內(nèi)的生活污廢水和企業(yè)工礦廢水通過總排干和排污口排入烏梁素海,一定程度上改變了湖泊水體的水化學(xué)組分。

    4 結(jié)論

    表4 3個(gè)湖泊采樣點(diǎn)對(duì)應(yīng)的ρ(Na+)/[ρ(Na+)+ρ(Ca2+)]和ρ(Cl-)/[ρ(Cl-)+ρ(HCO?)]比值Table 4 The value of ρ (Na+)/[ρ(Na+)+ρ(Ca2+)] and ρ(Cl-)/[ρ(Cl-)+ρ(HCO)] for three lakes

    表4 3個(gè)湖泊采樣點(diǎn)對(duì)應(yīng)的ρ(Na+)/[ρ(Na+)+ρ(Ca2+)]和ρ(Cl-)/[ρ(Cl-)+ρ(HCO?)]比值Table 4 The value of ρ (Na+)/[ρ(Na+)+ρ(Ca2+)] and ρ(Cl-)/[ρ(Cl-)+ρ(HCO)] for three lakes

    呼倫湖   達(dá)里湖烏梁素海點(diǎn)號(hào) ρ(Na+)/[ρ(Na+)+ρ(Ca2+)] ρ(Cl-)/[ρ(Cl-)+ρ(HCO3-)] 點(diǎn)號(hào) ρ(Na+)/[ρ(Na+)+ρ(Ca2+)]ρ(Cl-)/[ρ(Cl-)+ρ(HCO3-)] 點(diǎn)號(hào) ρ(Na+)/[ρ(Na+)+ρ(Ca2+)] ρ(Cl-)/[ρ(Cl-)+ρ(HCO3-)] A10  0.98  0.38  A7  1.00  0.48  J13 0.87  0.27 B9  0.98  0.33  B4  1.00  0.51  K12 0.90  0.20 D11  0.99  0.91  C5  1.00  0.70  L11 0.87  0.18 F9  0.99  0.87  C7  1.00  0.43  L15 0.89  0.83 G2  0.99  0.76  D4  1.00  0.51  M14 0.88  0.56 G8  0.99  0.90  D6  1.00  0.41  N13 0.89  0.51 H3  0.99  0.42  E1  1.00  0.51  O10 0.89  0.77 E2  1.00  0.50  P9 0.89  0.39 E5  1.00  0.36  P11 0.89  0.31 Q10 0.91  0.65 T5 0.95  0.42

    處于干旱-半干旱地區(qū)的內(nèi)蒙古高原湖泊,由于地理氣候的影響,蒸發(fā)強(qiáng)烈,降雨稀少,蒸發(fā)是降雨的10~20倍,導(dǎo)致該區(qū)域湖泊由淡水湖向微咸水和咸水湖泊轉(zhuǎn)化,水化學(xué)類型呼倫湖為Cl-·SO42--Na+,達(dá)里湖為Cl-·HCO3--Na+,烏梁素海為Cl-·HCO3--Na+。

    分析湖泊水中離子的組成發(fā)現(xiàn),湖泊水中各離子的質(zhì)量濃度呈顯著性的區(qū)域差異,達(dá)里湖各離子質(zhì)量濃度顯著高于呼倫湖和烏梁素海,與湖區(qū)所在的地理環(huán)境和水文條件密切相關(guān)。對(duì)同一湖區(qū),空間上湖泊水體離子質(zhì)量濃度也有差異,主要與湖泊水體的水動(dòng)力條件和補(bǔ)給排泄關(guān)系聯(lián)系緊密。

    與歷史資料的對(duì)比,并結(jié)合Piper解析了湖泊主要離子變化原因。利用Gibbs圖進(jìn)一步分析,3個(gè)湖泊水化學(xué)主控端元落在蒸發(fā)-結(jié)晶端元,同時(shí)3個(gè)湖泊采樣點(diǎn)均有部分落在Gibbs模型外,說明3個(gè)湖泊個(gè)別區(qū)域已受到了不同程度的人為活動(dòng)影響。

    參考文獻(xiàn):

    BERNER E K, BERNER R A. 1987. The global water cycle [J]. Englewood Cliffs: Prentice-Hall.

    GIBBS R J. 1970. Mechanisms controlling world water chemistr [J]. Science, 170(3962): 1088-1090.

    白占國(guó), 萬國(guó)江. 1998. 喀斯特區(qū)域的水化學(xué)不穩(wěn)定性——以黔中地區(qū)為例[J]. 水科學(xué)進(jìn)展, 12(4): 325-332.

    蔡永久, 姜加虎, 張路, 等. 2013. 長(zhǎng)江中下游湖群大型底棲動(dòng)物群落結(jié)構(gòu)及影響因素[J]. 生態(tài)學(xué)報(bào), 23(16): 4987-4999.

    曾海鰲, 吳敬祿. 2010. 蒙新高原湖泊水質(zhì)狀況及變化特征. 湖泊科學(xué)[J]. 22(6): 882-887.

    崔鳳麗, 李暢游, 史小紅, 等. 2013. 烏梁素海主要例子季節(jié)變化特征分析[J]. 干旱區(qū)資源與環(huán)境, 27(8): 137-142.

    付緒金, 賈克力, 史小紅, 等. 2013. 烏梁素海沉積物腐殖質(zhì)的組成及分布特征[J]. 湖泊科學(xué), 25(4): 489-496.

    國(guó)家環(huán)境保護(hù)總局. 2002. 水和廢水監(jiān)測(cè)分析方法[M]. 中國(guó)環(huán)境科學(xué)出版社.

    侯韶華, 徐海, 安芷生. 2009. 青海湖流域水化學(xué)主離子特征及控制因素初探[J]. 地球與環(huán)境, 37(1): 11-18.

    胡春華, 周文斌, 夏思奇. 2011. 鄱陽湖流域水化學(xué)主離子特征及其來源分析[J]. 環(huán)境化學(xué), 30(9): 1620-1626.

    李暢游, 孫標(biāo), 高占義, 等. 2011. 呼倫湖多波段遙感水深反演模型研究[J]. 水利學(xué)報(bào), 42(12): 1423-1431.

    陸瑩, 王乃昂, 李貴鵬, 等. 2010. 巴丹吉林沙漠湖泊水化學(xué)空間分布特征[J]. 湖泊科學(xué), 22(5): 774-782.

    汪敬忠, 吳敬祿, 曾海鰲, 等. 2013. 內(nèi)蒙古河套平原水體同位素及水化學(xué)特征[J]. 地球科學(xué)與環(huán)境學(xué)報(bào), 35(4): 104-111.

    王海雷, 鄭綿平. 2010. 青藏高原湖泊水化學(xué)與鹽度的相關(guān)性初步研究[J]. 地質(zhì)學(xué)報(bào), 54(10): 1517-1522.

    王鵬, 尚英男, 沈立成, 等. 2013. 青藏高原淡水湖泊水化學(xué)組成特征及其演化[J]. 環(huán)境科學(xué), 34(3): 874-881.

    王瑞久. 1983. 三線圖解及其水文地質(zhì)解釋[J]. 工程勘探, (6): 6-11.

    王蘇民. 1995. 呼倫湖古湖泊學(xué)研究[M]. 中國(guó)科學(xué)技術(shù)大學(xué)出版社.

    王志杰, 李暢游, 李衛(wèi)平, 等. 2012. 內(nèi)蒙古呼倫湖水量平衡計(jì)算與分析[J]. 湖泊科學(xué), 24(2): 273-281.

    葉宏萌, 袁旭音, 葛敏霞, 等. 2010. 太湖北部流域水化學(xué)特征及控制因素[J]. 生態(tài)環(huán)境學(xué)報(bào), 19(1): 23-27.

    張郝哲, 田明中, 郭婧, 等. 2012. 基于RS和GIS的內(nèi)蒙古達(dá)里諾爾湖1990─2010年動(dòng)態(tài)監(jiān)測(cè)[J]. 干旱區(qū)資源與環(huán)境, 26(10): 41-46.

    甄志磊, 李暢游, 李文寶, 等. 2014. 內(nèi)蒙古達(dá)里諾爾湖流域地表水和地下水環(huán)境同位素特征及補(bǔ)給關(guān)系[J]. 湖泊科學(xué), 26(6): 916-922.

    鄭綿平, 劉喜方. 2010. 青藏高原鹽湖水化學(xué)及其礦物組合特征[J]. 地質(zhì)學(xué)報(bào), 84(11): 1585-1600.

    周嘉欣, 丁永建, 曾國(guó)雄, 等. 2014. 疏勒河上游地表水水化學(xué)主離子特征及其控制因素[J]. 環(huán)境科學(xué), 35(9): 3315-3324.

    朱海勇, 陳永金, 劉加珍, 等. 2013. 塔里木河中下游地下水化學(xué)及其演變特征分析[J]. 干旱區(qū)地理, 36(1): 8-18.

    Major Ion Chemistry and Influencing Factors of Three Typical Lakes in Inner Mongolia Plateau

    WU Yong1, SHI Xiaohong1, ZHAO Shengnan1, LIN Taotao2, MA Jun3
    1. Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, Huhhot 010018, China; 2. Three Gorges Hydrology and Water Resources Survey Bureau, Bureau of Hydrology, Yangtze River Water Conservancy Commission, Yichang, 443000, China; 3. Honggebo Drainage Station of Hetao Irrigation District Management Bureau of Inner Mongolia, Bayannur 014400, China

    Abstract:Inner Mongolia plateau is located in the semi-arid and arid climate zone of northern China. The climate exhibits both semi-arid and arid characteristics as it spans across longitudes from East to West. Lakes in the region are mostly inland lakes. Their water chemistry characteristics are greatly influenced by geographical and climate. Three lakes, Hulun Lake, Dali Lake and Ulansuhai Lake, were analyzed for their main ions and the discrepancies, the chemical characteristics of the water, and their control factors. The results showed that the water chemical type, cation by Na+leading and anion by Cl-leading, was single by the geographical and climate. However, the main ions content was very different in these lakes. The percent of Na+in the total number of mole of cationic was 76.6%, 97.2% and 66.0%, while the percent of Cl-in the total number of mole of anion was 49.6%, 55.5% and 45.2%, respectively for three lakes. The coefficient of variation of Na+and SO42-, Cl-was relatively low, it showed the ion content was stable in the water. Based on the Piper triangle chart analyzed ions consist, the anion distributed in the line of SO42--Cl-and the cation in the line of Na++K+in Hulun lake and Ulansuhai Lake. But the anion was in the line of HCO3--Cl-and the cation in the line of Na++K+in Dali Lake. It shows that the water chemistry characteristics are affected by evaporation and crystallization and rock weathering. The spatial changes of the content of the main ions are also affected by the relationship between runoff and drainage, and are combined with the Gibbs model to further explain the different degree of human activity for these three typical lakes. Researching the lakes’ chemistry characteristics has a great significance for protecting lakes’ environment.

    Key words:lakes in the Inner Mongolia plateau; arid and semi-arid regions; mainly composed of ions; water chemistry characteristics

    收稿日期:2015-04-29

    *通信作者。E-mail: imaushixiaohong@163.com

    作者簡(jiǎn)介:吳用(1990年生),男,博士,主要從事流域水污染控制研究。E-mail: wuyong526@126.com

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51269016;51269017;51339002;51409288);內(nèi)蒙古農(nóng)業(yè)大學(xué)優(yōu)秀青年科學(xué)基金項(xiàng)目(2014XYQ-10)

    中圖分類號(hào):X131.2

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1674-5906(2015)07-1202-07

    DOI:10.16258/j.cnki.1674-5906.2015.07.019

    五月开心婷婷网| 老司机亚洲免费影院| 看免费av毛片| 精品久久国产蜜桃| 亚洲精品国产色婷婷电影| 狠狠婷婷综合久久久久久88av| 人妻少妇偷人精品九色| 三上悠亚av全集在线观看| 在线观看免费视频网站a站| 两个人免费观看高清视频| 侵犯人妻中文字幕一二三四区| 天堂俺去俺来也www色官网| 亚洲国产成人一精品久久久| 久久人人97超碰香蕉20202| 国产成人欧美| 日本欧美视频一区| 18禁动态无遮挡网站| 99久久人妻综合| 亚洲少妇的诱惑av| 亚洲精品一二三| 岛国毛片在线播放| 蜜桃在线观看..| 另类亚洲欧美激情| av视频免费观看在线观看| 久久久亚洲精品成人影院| 成年人免费黄色播放视频| 最近最新中文字幕大全免费视频 | 午夜免费鲁丝| 日韩免费高清中文字幕av| 卡戴珊不雅视频在线播放| av在线播放精品| a 毛片基地| 国产黄色视频一区二区在线观看| 亚洲国产欧美日韩在线播放| 日本91视频免费播放| 日日啪夜夜爽| 十八禁高潮呻吟视频| 久久婷婷青草| 国产有黄有色有爽视频| 看免费成人av毛片| 国产爽快片一区二区三区| 在线 av 中文字幕| 欧美 亚洲 国产 日韩一| 又黄又粗又硬又大视频| 欧美国产精品va在线观看不卡| 五月天丁香电影| 国产片特级美女逼逼视频| 国产精品嫩草影院av在线观看| 午夜福利影视在线免费观看| 国产精品熟女久久久久浪| 丰满少妇做爰视频| 久久国产亚洲av麻豆专区| 免费观看性生交大片5| 国产免费一区二区三区四区乱码| 9热在线视频观看99| 国产欧美日韩一区二区三区在线| 岛国毛片在线播放| 亚洲伊人久久精品综合| 久久国内精品自在自线图片| 久久国产精品男人的天堂亚洲 | 2021少妇久久久久久久久久久| 日韩制服丝袜自拍偷拍| 国产有黄有色有爽视频| xxx大片免费视频| 好男人视频免费观看在线| 欧美成人午夜免费资源| 黑丝袜美女国产一区| 精品一区二区免费观看| 一级毛片我不卡| 久久精品国产a三级三级三级| 18+在线观看网站| 一级片免费观看大全| 18禁裸乳无遮挡动漫免费视频| 大片免费播放器 马上看| 男人舔女人的私密视频| 久久国产精品大桥未久av| 国产极品粉嫩免费观看在线| 丝袜美足系列| 日本91视频免费播放| 欧美日韩一区二区视频在线观看视频在线| 在线观看三级黄色| 曰老女人黄片| 国产不卡av网站在线观看| 亚洲av国产av综合av卡| 国产成人91sexporn| 国产1区2区3区精品| 欧美成人精品欧美一级黄| 亚洲性久久影院| 欧美日韩成人在线一区二区| 啦啦啦视频在线资源免费观看| 在线亚洲精品国产二区图片欧美| 看免费成人av毛片| 十分钟在线观看高清视频www| 欧美日韩国产mv在线观看视频| 亚洲av在线观看美女高潮| 亚洲人成77777在线视频| 国产精品秋霞免费鲁丝片| 美女福利国产在线| 一个人免费看片子| av网站免费在线观看视频| 国产亚洲午夜精品一区二区久久| 三级国产精品片| 亚洲性久久影院| av免费观看日本| 一二三四中文在线观看免费高清| 丝袜喷水一区| av网站免费在线观看视频| 丰满乱子伦码专区| 欧美3d第一页| 欧美日韩综合久久久久久| 色婷婷av一区二区三区视频| a 毛片基地| 免费黄频网站在线观看国产| 亚洲国产毛片av蜜桃av| 成人国语在线视频| 一级毛片电影观看| a级毛片在线看网站| 亚洲精品乱码久久久久久按摩| 美女国产视频在线观看| √禁漫天堂资源中文www| 国产精品欧美亚洲77777| 国产成人免费无遮挡视频| 日本av手机在线免费观看| 国产在视频线精品| 日韩一区二区三区影片| 高清毛片免费看| 久久狼人影院| 人人妻人人爽人人添夜夜欢视频| 免费观看在线日韩| 国产黄频视频在线观看| 国产精品秋霞免费鲁丝片| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久精品精品| 高清欧美精品videossex| 999精品在线视频| 免费少妇av软件| 亚洲在久久综合| 国产黄色免费在线视频| 日韩伦理黄色片| 热re99久久国产66热| 久久久a久久爽久久v久久| 成年美女黄网站色视频大全免费| 欧美性感艳星| 欧美成人精品欧美一级黄| 大陆偷拍与自拍| 啦啦啦在线观看免费高清www| 熟女人妻精品中文字幕| 久久毛片免费看一区二区三区| 插逼视频在线观看| 久久精品久久久久久噜噜老黄| 免费在线观看黄色视频的| 有码 亚洲区| 少妇猛男粗大的猛烈进出视频| 丝袜脚勾引网站| 高清毛片免费看| 伦精品一区二区三区| 午夜福利网站1000一区二区三区| 婷婷色综合大香蕉| 又粗又硬又长又爽又黄的视频| 免费观看无遮挡的男女| 建设人人有责人人尽责人人享有的| 中文字幕亚洲精品专区| 欧美97在线视频| 老司机影院毛片| 99香蕉大伊视频| 狂野欧美激情性xxxx在线观看| a级毛片黄视频| 亚洲精品一二三| 亚洲国产av新网站| 亚洲欧美一区二区三区国产| 国产在线免费精品| 亚洲av福利一区| 免费黄色在线免费观看| 在线天堂中文资源库| 亚洲伊人久久精品综合| 亚洲五月色婷婷综合| 看免费成人av毛片| 欧美精品国产亚洲| 中文字幕制服av| 熟女电影av网| 精品久久国产蜜桃| 老熟女久久久| 两性夫妻黄色片 | 精品亚洲成国产av| 国产永久视频网站| av卡一久久| 久久这里有精品视频免费| 大香蕉久久成人网| 国产在视频线精品| 99re6热这里在线精品视频| 亚洲美女视频黄频| 亚洲人与动物交配视频| 人人妻人人澡人人爽人人夜夜| 国产男人的电影天堂91| 久久久久久久大尺度免费视频| 寂寞人妻少妇视频99o| 爱豆传媒免费全集在线观看| 日本av手机在线免费观看| av播播在线观看一区| 内地一区二区视频在线| 超碰97精品在线观看| h视频一区二区三区| 老司机影院毛片| 26uuu在线亚洲综合色| 搡老乐熟女国产| 高清欧美精品videossex| 九九爱精品视频在线观看| 国产精品国产av在线观看| 大陆偷拍与自拍| 少妇被粗大猛烈的视频| 久久97久久精品| 在线观看www视频免费| 国产精品人妻久久久久久| av网站免费在线观看视频| xxx大片免费视频| 亚洲综合精品二区| 亚洲精品美女久久av网站| 一级片'在线观看视频| 飞空精品影院首页| 亚洲一区二区三区欧美精品| 欧美日韩国产mv在线观看视频| 成人午夜精彩视频在线观看| 草草在线视频免费看| 中文乱码字字幕精品一区二区三区| 亚洲久久久国产精品| 少妇的丰满在线观看| 少妇人妻精品综合一区二区| 大陆偷拍与自拍| av天堂久久9| 日日啪夜夜爽| 26uuu在线亚洲综合色| 久久久久精品性色| 免费高清在线观看日韩| 国产成人av激情在线播放| 欧美精品高潮呻吟av久久| 在线观看免费日韩欧美大片| 欧美 亚洲 国产 日韩一| 99re6热这里在线精品视频| 亚洲av福利一区| 欧美日韩亚洲高清精品| 亚洲精品aⅴ在线观看| 亚洲,一卡二卡三卡| 亚洲丝袜综合中文字幕| 亚洲伊人久久精品综合| 国产精品蜜桃在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲综合精品二区| 国产精品久久久久成人av| 蜜桃在线观看..| 成人国语在线视频| 如何舔出高潮| 这个男人来自地球电影免费观看 | 99久久中文字幕三级久久日本| 侵犯人妻中文字幕一二三四区| 久久久久久久久久人人人人人人| 伦理电影免费视频| 在线观看一区二区三区激情| 欧美日韩视频精品一区| av女优亚洲男人天堂| 久久久精品免费免费高清| 最新中文字幕久久久久| 丝袜在线中文字幕| 视频在线观看一区二区三区| 亚洲精品日本国产第一区| 人妻一区二区av| 免费高清在线观看视频在线观看| 人妻人人澡人人爽人人| 免费大片黄手机在线观看| 国国产精品蜜臀av免费| 少妇精品久久久久久久| 哪个播放器可以免费观看大片| 日韩免费高清中文字幕av| 日韩,欧美,国产一区二区三区| 亚洲国产欧美在线一区| videos熟女内射| 亚洲天堂av无毛| 人妻少妇偷人精品九色| 久久女婷五月综合色啪小说| 综合色丁香网| 深夜精品福利| 99久久中文字幕三级久久日本| 最近的中文字幕免费完整| 亚洲精品456在线播放app| 在线观看免费日韩欧美大片| 成人黄色视频免费在线看| 久久久久视频综合| 丝袜脚勾引网站| 在线亚洲精品国产二区图片欧美| 90打野战视频偷拍视频| 777米奇影视久久| 黄色视频在线播放观看不卡| av在线老鸭窝| 久久久亚洲精品成人影院| 两个人看的免费小视频| 亚洲av综合色区一区| 国产成人精品福利久久| 美女国产视频在线观看| 免费日韩欧美在线观看| 国产免费又黄又爽又色| 久久久亚洲精品成人影院| 在线观看美女被高潮喷水网站| 国产爽快片一区二区三区| av.在线天堂| 成人影院久久| 婷婷成人精品国产| 日本黄色日本黄色录像| 国产国语露脸激情在线看| 午夜福利影视在线免费观看| 少妇被粗大的猛进出69影院 | 国产精品人妻久久久影院| 午夜91福利影院| 欧美 日韩 精品 国产| 人人妻人人澡人人爽人人夜夜| 国产午夜精品一二区理论片| 精品人妻在线不人妻| 欧美最新免费一区二区三区| 91国产中文字幕| 制服人妻中文乱码| 99热这里只有是精品在线观看| 亚洲中文av在线| 精品午夜福利在线看| 99热这里只有是精品在线观看| 看十八女毛片水多多多| 亚洲综合色网址| 亚洲五月色婷婷综合| 久久久久久久久久久免费av| 精品久久国产蜜桃| 久久久久久人妻| 高清黄色对白视频在线免费看| av天堂久久9| 中文字幕免费在线视频6| 久久女婷五月综合色啪小说| 大片电影免费在线观看免费| 丝袜脚勾引网站| 看免费成人av毛片| 精品久久蜜臀av无| 性色av一级| 日本黄大片高清| 国产日韩一区二区三区精品不卡| 少妇人妻精品综合一区二区| 久久精品人人爽人人爽视色| 亚洲一级一片aⅴ在线观看| 水蜜桃什么品种好| 五月开心婷婷网| 精品久久国产蜜桃| 亚洲精品一区蜜桃| 亚洲激情五月婷婷啪啪| 91在线精品国自产拍蜜月| 国产日韩欧美视频二区| 国产精品女同一区二区软件| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 国产永久视频网站| 啦啦啦中文免费视频观看日本| av黄色大香蕉| 男的添女的下面高潮视频| 免费观看性生交大片5| 在线观看一区二区三区激情| 91国产中文字幕| 国产在视频线精品| 国产精品不卡视频一区二区| 免费不卡的大黄色大毛片视频在线观看| 秋霞在线观看毛片| 岛国毛片在线播放| 午夜激情久久久久久久| 久久久久久人人人人人| 成年美女黄网站色视频大全免费| 久久热在线av| 老女人水多毛片| freevideosex欧美| 国产成人a∨麻豆精品| 新久久久久国产一级毛片| 最新的欧美精品一区二区| 国产1区2区3区精品| 人妻一区二区av| 国产高清国产精品国产三级| 在线精品无人区一区二区三| 看免费av毛片| 国产精品一区二区在线观看99| 亚洲 欧美一区二区三区| 不卡视频在线观看欧美| 精品一区在线观看国产| 在线观看免费视频网站a站| 一本色道久久久久久精品综合| 老熟女久久久| 日韩大片免费观看网站| 午夜福利视频在线观看免费| 色吧在线观看| 国产又爽黄色视频| 少妇的逼好多水| 午夜福利乱码中文字幕| 亚洲内射少妇av| 日韩电影二区| 一级毛片电影观看| 深夜精品福利| 日本爱情动作片www.在线观看| 亚洲成av片中文字幕在线观看 | 精品第一国产精品| 亚洲性久久影院| 国产精品女同一区二区软件| 日本黄色日本黄色录像| 国产成人a∨麻豆精品| videosex国产| 亚洲高清免费不卡视频| 中文字幕另类日韩欧美亚洲嫩草| 国产在视频线精品| 狠狠婷婷综合久久久久久88av| 一本久久精品| 高清黄色对白视频在线免费看| 国产老妇伦熟女老妇高清| 国产女主播在线喷水免费视频网站| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 妹子高潮喷水视频| 一区在线观看完整版| 九九爱精品视频在线观看| www.色视频.com| 国产亚洲精品第一综合不卡 | 80岁老熟妇乱子伦牲交| 久久热在线av| 色94色欧美一区二区| 99九九在线精品视频| 国产有黄有色有爽视频| 97人妻天天添夜夜摸| 亚洲精品乱码久久久久久按摩| 99热6这里只有精品| 亚洲色图 男人天堂 中文字幕 | av在线观看视频网站免费| 亚洲人与动物交配视频| 91aial.com中文字幕在线观看| 美女视频免费永久观看网站| 免费在线观看黄色视频的| 极品少妇高潮喷水抽搐| 午夜激情久久久久久久| 97超碰精品成人国产| 最近中文字幕高清免费大全6| 国产在线视频一区二区| 久久午夜福利片| 全区人妻精品视频| 男人爽女人下面视频在线观看| a级毛片在线看网站| 人妻系列 视频| 亚洲熟女精品中文字幕| 久久精品国产鲁丝片午夜精品| 日本vs欧美在线观看视频| 久久国产精品男人的天堂亚洲 | 女人久久www免费人成看片| 美女福利国产在线| 亚洲第一区二区三区不卡| www日本在线高清视频| 国产永久视频网站| 亚洲成人一二三区av| 最新中文字幕久久久久| 亚洲第一av免费看| 狂野欧美激情性xxxx在线观看| 亚洲成av片中文字幕在线观看 | 国产精品久久久久久精品电影小说| 亚洲激情五月婷婷啪啪| 国产欧美日韩一区二区三区在线| 日韩大片免费观看网站| 两个人免费观看高清视频| 黄色怎么调成土黄色| 亚洲中文av在线| 午夜福利影视在线免费观看| 亚洲精品自拍成人| 激情视频va一区二区三区| 91aial.com中文字幕在线观看| 国产深夜福利视频在线观看| 中文精品一卡2卡3卡4更新| 亚洲精品456在线播放app| 丝袜人妻中文字幕| 另类精品久久| 日韩不卡一区二区三区视频在线| 国产成人精品一,二区| 亚洲内射少妇av| 免费在线观看完整版高清| 人人妻人人添人人爽欧美一区卜| 成人亚洲精品一区在线观看| 大香蕉97超碰在线| 丝袜人妻中文字幕| 熟妇人妻不卡中文字幕| 免费久久久久久久精品成人欧美视频 | 亚洲中文av在线| 国产在视频线精品| 天天躁夜夜躁狠狠久久av| 少妇人妻精品综合一区二区| 久久精品国产自在天天线| 热re99久久精品国产66热6| 精品国产一区二区久久| 午夜影院在线不卡| 丰满迷人的少妇在线观看| 伦精品一区二区三区| 午夜福利视频精品| 熟女电影av网| 99热这里只有是精品在线观看| 亚洲精品,欧美精品| 丝袜在线中文字幕| 亚洲精品乱久久久久久| 国产精品蜜桃在线观看| 亚洲色图 男人天堂 中文字幕 | 极品少妇高潮喷水抽搐| 国产麻豆69| 飞空精品影院首页| 黄片无遮挡物在线观看| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 精品卡一卡二卡四卡免费| 久久久久久久久久成人| av一本久久久久| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 婷婷色av中文字幕| 国产黄频视频在线观看| 韩国高清视频一区二区三区| 日产精品乱码卡一卡2卡三| 九色亚洲精品在线播放| 一级,二级,三级黄色视频| 国产精品一区www在线观看| 欧美亚洲日本最大视频资源| 精品酒店卫生间| 久久热在线av| 久久毛片免费看一区二区三区| 国产男人的电影天堂91| 丁香六月天网| 国产欧美日韩综合在线一区二区| 国产精品久久久久久精品古装| 国产 一区精品| 欧美 日韩 精品 国产| 少妇被粗大猛烈的视频| 婷婷色av中文字幕| 久久久亚洲精品成人影院| 99热6这里只有精品| 亚洲精品视频女| 99热6这里只有精品| av播播在线观看一区| 久久久精品区二区三区| 美女大奶头黄色视频| 9热在线视频观看99| 国产免费一级a男人的天堂| 国产精品女同一区二区软件| 欧美老熟妇乱子伦牲交| 美女xxoo啪啪120秒动态图| 人妻 亚洲 视频| 多毛熟女@视频| 亚洲精品久久成人aⅴ小说| xxx大片免费视频| 国产69精品久久久久777片| 久久精品国产自在天天线| 男人舔女人的私密视频| 中文欧美无线码| 建设人人有责人人尽责人人享有的| 视频在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 免费观看av网站的网址| 久久久久精品性色| 国产福利在线免费观看视频| 国产片特级美女逼逼视频| 亚洲一码二码三码区别大吗| 色婷婷久久久亚洲欧美| 香蕉丝袜av| 欧美精品亚洲一区二区| 亚洲精品乱码久久久久久按摩| 麻豆乱淫一区二区| av免费观看日本| 成年av动漫网址| 99久久中文字幕三级久久日本| 亚洲精品一区蜜桃| 日本av免费视频播放| 婷婷成人精品国产| 久久这里只有精品19| 亚洲丝袜综合中文字幕| 一级毛片电影观看| kizo精华| 亚洲精品久久午夜乱码| 九草在线视频观看| videossex国产| 人人妻人人澡人人爽人人夜夜| 免费日韩欧美在线观看| 肉色欧美久久久久久久蜜桃| 高清欧美精品videossex| av网站免费在线观看视频| 激情视频va一区二区三区| 最黄视频免费看| 亚洲美女搞黄在线观看| 一级毛片 在线播放| 夜夜爽夜夜爽视频| 欧美国产精品一级二级三级| 考比视频在线观看| 久久久久久久精品精品| 天堂中文最新版在线下载| 国产精品一区www在线观看| 高清黄色对白视频在线免费看| 高清不卡的av网站| 久久久欧美国产精品| 精品国产乱码久久久久久小说| 国产无遮挡羞羞视频在线观看| 视频在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲国产精品专区欧美| 久久精品熟女亚洲av麻豆精品| 国产精品女同一区二区软件| 亚洲成色77777| 中文乱码字字幕精品一区二区三区| av在线观看视频网站免费| 高清毛片免费看| 2018国产大陆天天弄谢| 捣出白浆h1v1| 咕卡用的链子| 人成视频在线观看免费观看| 啦啦啦啦在线视频资源| 一区二区av电影网| 亚洲丝袜综合中文字幕| 久久99热这里只频精品6学生| 国产有黄有色有爽视频| 国产日韩欧美视频二区|