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

    遼河東西支流與干流水體光學特性和組分對比研究

    2015-08-30 03:00:41邵田田宋開山李思佳中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所中國科學院濕地生態(tài)與環(huán)境重點實驗室吉林長春3002河南大學黃河文明與可持續(xù)發(fā)展研究中心暨黃河文明傳承與現(xiàn)代文明建設(shè)河南省協(xié)同創(chuàng)新中心河南開封7500河南大學環(huán)境與規(guī)劃學院河南開封7500中國科學院地理科學與資源研究所北京000
    中國環(huán)境科學 2015年9期
    關(guān)鍵詞:遼河流域遼河干流

    邵田田,宋開山,丁 智,李思佳,趙 瑩,張 柏(.中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所,中國科學院濕地生態(tài)與環(huán)境重點實驗室,吉林 長春 3002;2.河南大學黃河文明與可持續(xù)發(fā)展研究中心暨黃河文明傳承與現(xiàn)代文明建設(shè)河南省協(xié)同創(chuàng)新中心,河南 開封 7500;3.河南大學環(huán)境與規(guī)劃學院,河南 開封 7500;.中國科學院地理科學與資源研究所,北京 000)

    遼河東西支流與干流水體光學特性和組分對比研究

    邵田田1,2,3,宋開山1*,丁 智4,李思佳1,趙 瑩1,張 柏1(1.中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所,中國科學院濕地生態(tài)與環(huán)境重點實驗室,吉林 長春 130102;2.河南大學黃河文明與可持續(xù)發(fā)展研究中心暨黃河文明傳承與現(xiàn)代文明建設(shè)河南省協(xié)同創(chuàng)新中心,河南 開封 475001;3.河南大學環(huán)境與規(guī)劃學院,河南 開封 475001;4.中國科學院地理科學與資源研究所,北京 100101)

    通過對東、西遼河和遼河干流有色溶解有機物(CDOM)吸收特性和熒光特征以及溶解有機碳(DOC)濃度的分析,對比研究不同子流域CDOM的光學特性差異,并分析影響CDOM與DOC濃度的主要因素.對比研究發(fā)現(xiàn),遼河干流CDOM與DOC相關(guān)性(R = 0.89,P < 0.01)高于西遼河(R = 0.81)與東遼河(R = 0.75).總懸浮物濃度(TSM)和總堿度與不同子流域CDOM/DOC相關(guān)性較高.不同流域內(nèi)CDOM與DOC的濃度變化較大,西遼河 CDOM 和 DOC濃度明顯高于東遼河和遼河干流.流域內(nèi)土地利用、氣候條件對兩者濃度的影響顯著,林地和CDOM及DOC濃度存在顯著負相關(guān)(R = -0.41, -0.56),即林地面積越大,CDOM和DOC濃度越低;而農(nóng)田與兩者存在正相關(guān)(R = 0.40, 0.32),農(nóng)田面積越大,相應(yīng)的CDOM和DOC濃度越高;降雨量與CDOM和DOC呈現(xiàn)明顯的負相關(guān)(R = -0.53, -0.38),說明降雨對河流CDOM的稀釋作用較大.對于CDOM熒光強度Fn(355)與CDOM濃度來說,在每個子流域兩者均存在較好的相關(guān)性,特別是西遼河流域(R = 0.96, P <0.01).CDOM的三維熒光光譜發(fā)現(xiàn),東、西遼河及遼河干流均表現(xiàn)出較強的類腐殖酸熒光峰(A峰和C峰),同時遼河干流表現(xiàn)出很強的類蛋白質(zhì)峰(T峰).CDOM的吸收斜率(S)以及基于熒光光譜的FI和HIX表明東、西遼河和遼河干流CDOM的主要來源是外源高等植物輸入的大分子量DOM,但是,西遼河和遼河干流CDOM的分子量要小于東遼河CDOM組成物質(zhì).

    東遼河;西遼河;CDOM;DOC;土地利用;熒光

    有機溶解有機物(DOM),在水環(huán)境物質(zhì)循環(huán)、能量運輸以及水生生態(tài)系統(tǒng)中起到非常重要的作用,一般使用溶解有機碳(DOC)來表征DOM的濃度[1-3].有色溶解有機物(CDOM)是DOM的重要組成部分,它能夠強烈吸收紫外輻射,限制UV-B輻射在水體中的傳播,保護水中生態(tài)系統(tǒng);同時影響紫外輻射在水中的傳輸并決定水體固有光學特性,從而限制了遙感在地表水文化學研究方面的應(yīng)用[4-8].

    CDOM擁有較為獨特的光譜吸收特性及熒光特征,被廣泛的用于研究海洋、港灣、河口、湖泊、河流的 DOM 組成和來源的研究中[9-12].河流連接著陸地和海洋兩大生態(tài)系統(tǒng),是陸源物質(zhì)向海輸送的重要途徑,對海岸帶及鄰近的海洋系統(tǒng)影響深遠.而河流的DOM含量不僅反映著流域的自然狀態(tài),人類活動對流域地表的改造作用也會在河流的DOM含量中得以體現(xiàn)[13].河流DOM 主要來源于外源物質(zhì)的輸入,因此流域特征與 DOM 的濃度和特征聯(lián)系緊密[14-15].流域內(nèi)土地利用類型的變化嚴重影響著水中CDOM的濃度和光譜特征.另外,流域內(nèi)的氣候條件、水文條件也會對CDOM的特征產(chǎn)生很大影響,影響著河流有機碳含量向海的輸入[16].

    遼河流域面積為22.9×104km2,流經(jīng)河北、內(nèi)蒙古、吉林、遼寧等省份,流域特征的復(fù)雜性,土地利用類型的多樣性,氣候條件的多變性都對該區(qū)的CDOM和DOC產(chǎn)生重要影響.另外,水系流經(jīng)沈陽、撫順、鞍山、本溪、遼陽等大中型工業(yè)城市,隨著工農(nóng)業(yè)和城市的迅速發(fā)展,人為因素對遼河流域水環(huán)境的影響越來越大,所產(chǎn)生的危害也逐年增加,使遼河處于較為嚴重的污染狀態(tài)[17-18],從而影響水體中 CDOM 的光學吸收特性和DOC的濃度.因此本文利用CDOM的光譜特征以及DOC濃度來研究各子流域內(nèi)CDOM的吸收特性及其與相應(yīng) DOC的相關(guān)分析、CDOM吸收特性和熒光特征的相關(guān)關(guān)系;各子流域內(nèi),對CDOM和DOC產(chǎn)生影響的環(huán)境因素進行綜合分析以及東、西遼河及遼河干流 CDOM的組成特征和來源分析.

    1 材料與方法

    1.1研究區(qū)概況

    遼河流域位于中國東北地區(qū)南部(116°30′~125°47′E, 38°43′~45°N),北與松花江流域接壤,南與渤海灣相接,中部為遼河平原及松遼平原的一部分.流域由 2個獨立水系組成:一為東、西遼河水系,二是渾河、太子河水系.本文的研究區(qū)域主要集中于第一水系,即東、西遼河,于福德店匯流為遼河干流,經(jīng)雙臺子河由盤山入海,全長1.39×103km,其中干流長 516km.西遼河是遼河的最大支流,容納 4條支流:西拉木倫河、老哈河、新開河和教來河;東遼河發(fā)源于長白山系吉林哈達嶺支脈薩哈嶺,是遼河東部第一大支流,由遼河掌、渭津河、拉津河等70多條河流匯聚而成.本研究區(qū)氣候?qū)儆诎霛駶?半干旱的溫帶氣候區(qū),年平均溫度由南向北降低(4~7℃),降水量則自東向西遞減(350~1000mm).遼河流域上游山丘區(qū),多為黃白土和風沙土,水土流失嚴重,植被覆蓋度低于30%,是中國東北地區(qū)風沙干旱較為嚴重的地區(qū). 1.2 野外樣品的采集

    表1 采樣點編號及位置點Table 1 Numbers and locations of sampling sites

    圖1 采樣點位置分布以及研究區(qū)土地利用Fig.4 Location of sampling stations and the land use/cover of Liaohe River watershed

    于2012年10月對西遼河、東遼河以及遼河干流進行采樣,采集22組;并于2013年10月對西遼河進行補充采樣(2組).其中西遼河 5組數(shù)據(jù),東遼河9組數(shù)據(jù),遼河干流10組數(shù)據(jù),共計24組數(shù)據(jù),采樣點具體位置如表1和圖1所示.在采樣過程中,每個樣點采集 2500mL水樣并立刻放在車載冰箱里冷藏,送至實驗室測量其光學參數(shù)以及水質(zhì)參數(shù),整個測試過程需要2~3d.

    1.3CDOM吸收光譜指數(shù)

    CDOM吸收系數(shù)的測定是通過0.22μm的微孔聚碳酸酯膜過濾水樣后,利用 UV-2600紫外分光光度計測得200~800nm的吸光度,計算得到各波長的吸收系數(shù),通過校正得到用于分析的吸收系數(shù).本文采用aCDOM(355)作為CDOM的濃度[19].

    S值是 CDOM 吸收曲線的斜率,可以反映CDOM 組成分子的大小,表征 CDOM的來源[20-21].本研究選擇275~295nm,采用線性擬合方法進行S值的求算[22].

    CDOM的三維熒光數(shù)據(jù)(EMMs)通過日立F-7000熒 光 光 度 計 (Hitachi Fluorescence spectrophotometer F-7000)測定獲得,掃描光譜進行儀器自動校正.EMMs主要受到瑞利散射、拉曼散射以及內(nèi)部濾波器的影響.其中,拉曼散射可以通過扣除空白平均值的方法進行消除;瑞利散射主要是將發(fā)射波長≤激發(fā)波長+5nm以及發(fā)射波長≥激發(fā)波長+300nm的2個區(qū)域內(nèi)的數(shù)值全部用0值來取代;對儀器的激發(fā)與發(fā)射光譜的校正主要是去除儀器內(nèi)部濾波器的影響,從而得到CDOM真實熒光光譜[23-24].熒光強度采用同步測定的硫酸奎寧作參照,在激發(fā)波長 355nm、發(fā)射波長450nm處的熒光值作為該點的熒光強度,用QSU單位表示[25].另外,溶液中的有機污染物、離子強度等會對CDOM的熒光強度產(chǎn)生影響[26-28],也有學者認為,自然河流中的離子強度對河流CDOM的三維熒光光譜并無明顯影響[29].本研究為淡水河流,因此未考慮離子強度等對CDOM熒光的影響.

    腐殖化指數(shù)[13](HIX)以及熒光指數(shù)[24](FI)基于EMMs來計算.其中HIX為在254nm激發(fā),發(fā)射波長435 ~ 480nm與300 ~ 345nm的波段內(nèi)的熒光強度積分值的比值,HIX < 2表示不存在腐殖化的物質(zhì),而HIX > 10則表示存在典型的富里酸物質(zhì);FI為在370nm處激發(fā),發(fā)射波長450nm和500nm的比值,即450/500[24].當FI < 1.3表示水體中主要以陸源的、高等植物DOM為主,而FI值在1.8左右表示主要以微生物所產(chǎn)生DOM為主.

    1.4 其他參數(shù)的測定

    葉綠素 a(Chla)濃度通過 90%的丙酮溶液萃取,并用UV-2600紫外分光光度計分別測得630, 647,664,750nm處的吸光度,進而計算 Chla的濃度[30-31].總懸浮顆粒物濃度(TSM)的測定采用稱重法獲得;總氮(TN)通過堿性過硫酸鉀紫外分光光度法測得;總磷(TP)采用鉬藍分光光度法測得,具體步驟詳見文獻[31].電導由電導率儀測定(上海雷磁,DDS-307A);濁度通過紫外分光光度計測得;總堿度采用國標甲基橙滴定法獲得.溶解有機碳(DOC)濃度是通過總?cè)芙馓迹═DC)減去溶解無機碳(DIC)得到,而TDC和DIC則通過總有機碳分析儀(日本島津公司,TOC-VCPN)測得.

    2 結(jié)果

    2.1CDOM吸收特征以及與DOC的關(guān)系

    通過對東、西遼河與遼河干流CDOM的吸收光譜進行分析發(fā)現(xiàn)所有樣品的吸收存在共同點,即在700nm之后基本為0,而在較短波段(280~500nm)呈指數(shù)增長(如圖2).雖然3個流域CDOM的吸收光譜趨勢基本一致,但數(shù)值差別較大,基本趨勢是西遼河>遼河干流>東遼河.位于西遼河流域的2個樣點(S23、S24)的CDOM吸收特征與東遼河的類似(圖 2a),在較短波段隨波長增加減小的較為緩慢.

    對各子流域CDOM與DOC分別進行回歸分析發(fā)現(xiàn),DOC濃度隨著aCDOM(355)的增加而增大,西遼河DOC濃度最高,遼河干流次之,東遼河的最低(圖3).可以看出,每個子流域的aCDOM(355)與 DOC都存在線性相關(guān)性,其中遼河干流的相關(guān)性大于西遼河和東遼河.對所有采樣點數(shù)據(jù)進行相關(guān)分析發(fā)現(xiàn),其決定系數(shù)略小于遼河干流,但明顯高于東、西遼河.總的來看,遼河干流數(shù)據(jù)較為規(guī)律,RMSE為0.63,小于其他2個流域;西遼河數(shù)據(jù)較為離散,主要是因為位于西遼河上游的 2點(表1,圖2,3)較為不一致;而東遼河流域數(shù)據(jù)整體分布較為離散.

    通過對基于CDOM吸收所計算的兩個參量S275-295統(tǒng)計分析(表 2)可以看出,大部分數(shù)據(jù)的S275-295在0.0170~0.0194之間,僅僅有3個點(S06、S07和S24)的值較小,而這3個采樣點均位于水庫的上下游;從 S275-295的均值來看,西遼河CDOM略大于遼河干流,而東遼河則小于西遼河與遼河干流CDOM的斜率.

    圖2 東、西遼河及遼河干流CDOM吸收特征Fig.4 Characteristic of CDOM absorption in different sub-watersheds, the West Liaohe River watershed, East Liaohe River watershed and the Main Stem Liaohe River watershed

    圖3 αCDOM(355)與DOC濃度的相關(guān)性分析Fig.4 Correlation between αCDOM(355) and DOC concentration in different sub-watersheds

    表2 CDOM的吸收和熒光光譜指數(shù)的統(tǒng)計Table 1 Spectroscopic indices based on CDOM absorption and fluorescence

    2.2CDOM的熒光特征及其與吸收的關(guān)系

    2.2.1CDOM的熒光強度與吸收的關(guān)系 CDOM 的化學特性復(fù)雜,組成成分很多,但其與熒光強度有很強的相關(guān)性.本研究將各子流域的數(shù)據(jù)分別對Fn(355)和aCDOM(355)進行回歸分析,如圖4.西遼河aCDOM(355)和熒光強度的相關(guān)性最好,東遼河次之,遼河干流相對來說R較低.同時將3個流域的數(shù)據(jù)進行分析,其相關(guān)性則高于各子流域兩者的相關(guān)性.

    2.2.2CDOM的熒光特征 本文對24組數(shù)據(jù)的三維熒光光譜進行計算分析,得到樣品的2個指數(shù)(表2)和熒光峰值(表3).根據(jù)表2中基于EMMs所計算的兩個指數(shù)HIX和FI可以看出,西遼河HIX的均值為6.91,明顯大于遼河干流(6.37)和東遼河(4.44);而FI則基本相等.熒光峰主要包括2大類,類腐殖質(zhì)熒光峰(紫外區(qū)類腐殖酸熒光峰A,Ex/Em=230~270/440~455;可見光區(qū)類腐殖酸熒光峰 C,Ex/Em=305~340/410~450)和類蛋白質(zhì)熒光峰(類酪氨酸熒光峰 B,Ex/Em=220~230(270~280)/300~320;類色氨酸熒光峰 T,Ex/Em=220~230(270~280)/320~360).表 3表明各子流域水樣表現(xiàn)出較強的類腐殖酸熒光峰,而遼河干流同時也表現(xiàn)出較強的類蛋白質(zhì)熒光峰.

    圖4 aCDOM(355)與Fn (355)的相關(guān)性分析Fig.4 Correlation between aCDOM(355) and Fn (355)

    3 討論

    3.1CDOM和DOC濃度及相關(guān)性的影響因素

    3.1.1水質(zhì)參數(shù)對CDOM和DOC的影響 影響CDOM與DOC相關(guān)性的因素較多,但如果CDOM在DOC中的比例較恒定或者CDOM與無色溶解有機物(UDOM)呈比例時,兩者就會呈現(xiàn)較好的關(guān)系[31].因此,本研究引入 CDOM 在DOC中所占的比例,即CDOM/DOC,用于分析兩者相關(guān)性的影響因素,各流域的統(tǒng)計值如表4所示.可以看出,3流域中遼河干流CDOM/DOC變化較小,西遼河兩者的比例波動較大,東遼河CDOM/DOC的比率變化雖低于西遼河,但也波動明顯,與圖3所得到結(jié)果較為一致.而所有數(shù)據(jù)DOC與CDOM相關(guān)性比較好的原因是在較大的流域內(nèi),DOC濃度存在一定的梯度,遼河干流較高的相關(guān)性提升了總體數(shù)據(jù)的相關(guān)性.

    表4 CDOM/DOC比值統(tǒng)計Table 1 CDOM/DOC in the three watersheds

    東遼河CDOM/DOC與所有水質(zhì)參數(shù)均不存在顯著相關(guān)性,這說明在東遼河流域,影響CDOM與DOC濃度相關(guān)性的因素不確定,不同的水質(zhì)參數(shù)對兩者關(guān)系產(chǎn)生不同的影響,從而導致圖 3a中數(shù)據(jù)較為離散,相關(guān)性較弱.遼河干流CDOM/DOC與總堿度存在顯著正相關(guān),與葉綠素a、總懸浮物負相關(guān),說明遼河干流影響CDOM 與DOC相關(guān)性的決定因素是總堿度、葉綠素a和總懸浮物濃度.西遼河CDOM/DOC與TN、TP、總堿度存在顯著正相關(guān),與總懸浮物正相關(guān),與其他2個子流域不同的是,在西遼河流域TN和TP 對CDOM與DOC相關(guān)性影響較大.對比表5中不同子流域各參數(shù)的對數(shù)值與 CDOM/DOC相關(guān)性以及圖 3可以發(fā)現(xiàn),遼河干流和所有數(shù)據(jù)DOC與 CDOM相關(guān)性較好的主要原因是與TSM和總堿度存在明顯的相關(guān)性,西遼河流域并沒有發(fā)現(xiàn)很好的相關(guān)性.主要原因在于西遼河的S23(老哈河)CDOM和DOC變化與西遼河流域其他點的非同步性,其CDOM/DOC為0.14,小于均值 0.56.因此,除水質(zhì)參數(shù)外,其它因素也會對DOC與CDOM的相關(guān)性造成影響.

    3.1.2不同土地類型對CDOM及DOC的影響河水流經(jīng)不同的土地類型所攜帶的 CDOM和DOC濃度是不相同的,流經(jīng)濕地的河流DOM濃度一般較高,而流經(jīng)森林和農(nóng)田的河流DOM在不同的研究區(qū)域表現(xiàn)出不同的濃度,大小不一致.如歐洲中部的河流,流經(jīng)農(nóng)田區(qū)的DOM要高于流經(jīng)森林的濃度[33];而對于加拿大安大略州中南部的河流,流經(jīng)森林的 DOM濃度要高于農(nóng)田的濃度[16].因此,不同地區(qū)的土地利用類型對河流DOM濃度的影響是不一致的.本文基于DEM數(shù)據(jù)采用SWAT模型對遼河流域進行子流域的劃分,相應(yīng)的落在本研究區(qū)的與采樣點對應(yīng)的子流域有12個(圖1).在進行子流域的大小與DOC及CDOM濃度分析時發(fā)現(xiàn)子流域面積與兩者并不存在顯著相關(guān)性(CDOM: R = -0.37; DOC: R = -0.27),這與 Gareber等[33]研究結(jié)果一致.通過子流域各土地利用類型面積的統(tǒng)計發(fā)現(xiàn),除子流域1和 2之外,其余子流域的主要利用類型均為農(nóng)田,而且面積超過子流域面積的60%.由于子流域主要土地利用類型較為單一,而且面積差距較大,因此并不能突顯出子流域土地利用對 DOC及CDOM的作用.

    表5 水質(zhì)參數(shù)對CDOM/DOC的影響Table 1 Impacts of water quality parameters on CDOM/DOC

    為了進一步分析土地利用類型對CDOM及DOC的影響,本文對采樣點做了以5km為半徑的緩沖區(qū),在排除子流域內(nèi)水流路徑等水文過程影響的前提下分析緩沖區(qū)內(nèi)土地利用對該采樣點CDOM和DOC濃度的影響[33].本文將各采樣點所對應(yīng)的緩沖區(qū)的土地利用類型與 CDOM和DOC濃度進行相關(guān)分析(表6).可以看出,林地、農(nóng)田兩種類型對CDOM和DOC濃度有較大影響,其中林地與DOC存在顯著負相關(guān),與CDOM負相關(guān),而農(nóng)田與兩者均存在正相關(guān).由此可見,農(nóng)田面積越大,CDOM和DOC濃度越高;相反的,林地面積越大,DOC和CDOM濃度則越小,這與前人研究結(jié)果一致[33].出現(xiàn)這一現(xiàn)象的主要原因是,對于本研究區(qū),農(nóng)業(yè)土壤中的有機碳是水中DOC和CDOM的主要來源.農(nóng)業(yè)有機化肥的使用、農(nóng)田耕作方式等的影響使土壤有機碳轉(zhuǎn)變?yōu)槿芙庥袡C碳[33-34],存儲于水體中,使得農(nóng)田與河流的DOC和CDOM存在正相關(guān)關(guān)系;而森林因為有植被覆蓋,水土保持較農(nóng)田區(qū)好,土壤中的有機碳對河流溶解有機物的輸入較少,因而其CDOM 和 DOC濃度較低,相應(yīng)的也可以在草地與CDOM和DOC的相關(guān)性上發(fā)現(xiàn)此規(guī)律.另外,草地與CDOM/DOC負相關(guān),農(nóng)田與之則是正相關(guān),而水體面積與之則是顯著負相關(guān).從而說明,在小面積范圍內(nèi)草地、農(nóng)田和水體對CDOM和DOC相關(guān)性的影響較為突出.

    表6 土地利用類型與DOC、CDOM和Fn (355)的相關(guān)性Table 1 Correlations between land use/cover types and DOC, CDOM and Fn (355)

    另外,土地利用類型的影響也解釋了老哈河CDOM和DOC變化與西遼河流域其他點的非同步性.S23位于內(nèi)蒙古赤峰市紅山水庫下游, 其 CDOM的濃度受水庫影響較大,水庫水體更新較慢,所受日照時間較流動水體長,因而光化學反應(yīng)較為強烈,從而導致CDOM吸收衰減較慢,濃度較低[35];其周圍環(huán)境主要是蘆葦濕地,而濕地 DOC濃度較其他土地利用類型影響下的DOC濃度高[16,33].因此,若將點S23排除,西遼河流域DOC與CDOM的相關(guān)性則明顯提高,R = 0.99.同時也解釋了老哈河與西拉沐倫河的CDOM吸收特性與東遼河的類似而與西遼河其他樣點不同(圖 2a)的主要原因在于兩個采樣點所處的地理位置,與東遼河類似,比較靠近水源地,周圍土地利用類型主要是草原和林地(圖 1),耕地較少,受到人為擾動(城市污染、農(nóng)業(yè)污染等)較小,水體保持較為原始的狀態(tài),因而CDOM的濃度均較低.

    3.1.3降雨量對CDOM及DOC的影響 根據(jù)整個遼河流域氣象站點降雨量數(shù)據(jù)對遼河流域進行插值,從而獲得每個采樣點的降雨量數(shù)據(jù),分析其與CDOM和DOC濃度的相關(guān)性(圖5).本文水樣采集時間為10月上旬,利用9月份的降雨量數(shù)據(jù)與CDOM和DOC進行分析.需要注意的是,S23和S24為2013年的補充樣點,因此采用2013年 9月份降雨量數(shù)據(jù)進行分析.可以看出,CDOM和DOC的濃度隨降雨量的增加而降低,降雨量與 CDOM濃度存在顯著負相關(guān)(R = -0.53),與 DOC負相關(guān)(R = -0.38),這與 Huang等[34]的研究結(jié)果一致.遼河流域降雨量由西向東逐漸增加,而由南向北則遞減.CDOM和DOC濃度則是西遼河大于遼河干流,東遼河最低.濃度較高的西遼河與濃度較低的東遼河交匯后形成濃度低于西遼河而高于東遼河的干流濃度.西遼河濃度較高的原因主要有兩點:第一,西遼河流域植被覆蓋度較低,風沙較大,土壤侵蝕較為嚴重,土壤中的有機物更容易進入河流形成溶解有機物;第二,西遼河流域降雨較少,減輕了對水中溶解有機物的稀釋作用,同時流域內(nèi)較高的潛在蒸散量增加了對水中溶解有機物濃縮作用[35].另外,降雨量與CDOM濃度的負相關(guān)性也說明了降雨的稀釋作用要大于地表徑流所攜帶的 DOM入河流的作用.

    圖5 降雨量與CDOM和DOC濃度的關(guān)系Fig.4 Relationships between precipitation and CDOM and DOC

    綜上分析,影響CDOM與DOC濃度的因子較多,包括氣象數(shù)據(jù)、土地利用數(shù)據(jù)和水質(zhì)參數(shù)數(shù)據(jù).不同流域內(nèi),其影響因子也存在很大差別.

    3.1.4CDOM的組分與來源分析 由表2可知,所有樣點的FI均小于1.3,這表明遼河水體主要以陸源的、高等植物 DOM為主,而以微生物降解所產(chǎn)生的小分子DOM較少,并且遼河FI值小于 McKnight等[23]所研究的 Missouri河(1.5)、Ogeechee河(1.4)、Ohio河(1.5)以及Yakima河(1.4 - 1.5).CDOM的HIX均大于2,說明所有樣點均存在腐殖化的物質(zhì),西遼河S09的HIX為11.22說明該樣點存在典型的富里酸物質(zhì).另外,西遼河與遼河干流的HIX均值明顯高于東遼河,說明了東遼河CDOM的腐殖化程度相對較低,從而也說明東遼河CDOM的分子量相對較大.另外東遼河S275-295均值明顯低于西遼河與遼河干流,進一步說明東遼河 CDOM主要是由大分子物質(zhì)組成.需要注意的是,西遼河S23和S24的HIX值與東遼河較為相近,主要原因是它們均位于較為靠近源頭的位置,腐殖化程度較低.另外,HIX與CDOM濃度存在顯著正相關(guān),說明CDOM濃度越高,其包含的腐殖化物質(zhì)越多.

    由表 3可以看出,類蛋白質(zhì)峰是遼河干流CDOM的重要組成部分,這是與東、西遼河所不同的.一般認為B、T峰通常反映的是生物降解來源的色氨酸和酪氨酸所形成的熒光峰值,而C、A峰則反映的是外源輸入的腐殖酸和富里酸形成的熒光峰值[37-38].因此,對于遼河水體,外源輸入是DOM的主要來源,這與根據(jù)FI得到結(jié)果一致.另外,也有部分研究認為浮游植物生長過程中也會產(chǎn)生 C、A峰等類腐殖質(zhì)熒光[39],但是在本研究中葉綠素含量與A峰、C峰不存在任何相關(guān)性,因此可以基本排除類腐殖質(zhì)熒光峰主要是由浮游植物生長產(chǎn)生的可能.還有學者研究認為,河流CDOM一般僅僅表現(xiàn)出類腐殖酸峰[40-41],但是如果受人類活動污染影響較大的話則會在3D-EEMs上表現(xiàn)出很強的類蛋白質(zhì)峰[37,41-42].遼河流經(jīng)沈陽、遼陽以及鐵嶺等城市,受到較強的人為干擾,干流接納了城市的生活污水、生產(chǎn)廢水等以及各支流所攜帶的污染物,因此在三維熒光光譜上表現(xiàn)出很強的類蛋白質(zhì)峰;而東、西遼河受到人為干擾相對較小,因此蛋白質(zhì)峰不明顯.另外,本文雖對CDOM和DOC產(chǎn)生影響的各項因素進行分析,但也主要是針對水質(zhì)參數(shù)、土地類型和降雨量,而流域的其他特征還需要進一步的細化與分析,這也是以后要研究的重點.

    4 結(jié)論

    4.1西遼河CDOM和DOC濃度明顯高于東遼河和遼河干流,遼河干流CDOM與DOC的相關(guān)性明顯高于東、西遼河.

    4.2TSM和總堿度與 CDOM/DOC顯著相關(guān);林地和農(nóng)田面積對CDOM和DOC濃度影響較大,林地面積越大,兩者濃度越低,農(nóng)田面積越大,兩者濃度越高;降雨量與CDOM和DOC呈現(xiàn)明顯的負相關(guān).Fn(355)與 CDOM在每個子流域均存在較好的相關(guān)性,特別是西遼河流域.

    4.3通過CDOM的吸收斜率S和基于熒光光譜的FI、HIX發(fā)現(xiàn)東、西遼河和遼河干流CDOM的主要來源是外源高等植物 DOM的輸入,組成也傾向于較大的分子量物質(zhì),但西遼河和遼河干流相對與東遼河CDOM組成物質(zhì)的分子量要小.

    4.4通過三維熒光光譜分析發(fā)現(xiàn),東、西遼河及遼河干流均表現(xiàn)出較為強烈的類腐殖酸熒光峰,而遼河干流同時也表現(xiàn)出較強的類蛋白質(zhì)熒光峰,這主要與人類活動影響較大、河流攜帶較多的污染物有關(guān).

    [1] Benner R. Chemical composition and reactivity of dissolved organic matter [C]. D.A. Hansell and C.A. Carlson (eds.). Biogeochemistry of marine dissolved organic matter. Academic Press, 2002,59-90.

    [2] Spencer RGM, Stubbins A, Hernes PJ, et al. Photochemical degradation of dissolved organic matter and dissolved lignin phenols from the Congo River [J]. Journal of Geophysical Research: Biogeosciences (2005-2012), 2009,114(G3).

    [3] Zhang YL, Zhang EL, Yin Y, et al. Characteristics and sources of chromophoric dissolved organic matter in lakes of the Yungui Plateau, China, differing in trophic state and altitude [J]. Limnology and Oceanography, 2010,55(6):2645-2659.

    [4] Cory N, Buffam I, Laudon H, et al. Landscape control of stream water aluminum in a boreal catchment during spring flood [J]. Environmental Science and Technology, 2006,40(11):3494-3500.

    [5] Zhu WN, Tian YQ, Yu Q, et al. Using Hyperion imagery to monitor the spatial and temporal distribution of colored dissolved organic matter in estuarine and coastal regions [J]. Remote Sensing of Environment, 2013,134:342-354.

    [6] Stedmon C A, Markager S, Tranvik L, et al. Photochemical production of ammonium and transformation of dissolved organic matter in the Baltic Sea [J]. Marine Chemistry, 2007,104(3/4):227-240.

    [7] Tranvik L, Downing J A, Cotner J B, et al. Lakes and reservoirs as regulators of carbon cycling and climate [J]. Limnology and Oceanography, 2009,54(6):2298-2314.

    [8] Del Vecchio R, Blogh N V. Spatial and seasonal distribution of chromophoric dissolved organic matter and dissolved organic carbon in the Middle Atlantic Bight [J]. Marine Chemistry, 2004, 89(1):169-187.

    [9] Chen Z, Hu C, Conmy RN, et al. Colored dissolved organic matter in Tampa Bay, Florida [J]. Marine Chemistry, 2007,104(1):98-109.

    [10] Carstea E M, Baker A, Pavelescu G, et al. Continuous fluorescence assessment of organic matter variability on the Bournbrook River, Birmingham, UK [J]. Hydrological Processes, 2009,23(13):1937-1946.

    [11] Siegel H, Stottmeister I, Reiβmann J, et al. Siak River System—East-Sumatra: Characterisation of sources, estuarine processes, and discharge into the Malacca Strait [J]. Journal of Marine Systems, 2009,77(1):148-159.

    [12] Singh S, D'Sa EJ, Swenson EM. Chromophoric dissolved organic matter (CDOM) variability in Barataria Basin using excitationemission matrix (EEM) fluorescence and parallel factor analysis (PARAFAC) [J]. Science of the Total Environment, 2010,408(16):3211-3222.

    [13] Zsolnay A, Baigar E, Jimenez M, et al. Differentiating with fluorescence spectroscopy the sources of dissolved organic matter in soils subjected to drying [J]. Chemosphere, 1999,38(1):45-50.

    [14] Aitkenhead-Peterson J A, McDowell W H, Neff J C, et al. Sources, production, and regulation of allochthonous dissolved organic matter inputs to surface waters (pp. 71-91) [M]. San Diego: Academic Press., 2003.

    [15] Williams C J, Yamashita Y, Wilson H F, et al. Unraveling the role of land use and microbial activity in shaping dissolved organic matter characteristics in stream ecosystems [J]. Limnology and Oceanography, 2010,55(3):1159.

    [16] Wilson H F, Xenopoulos M A. Ecosystem and seasonal control of stream dissolved organic carbon along a gradient of land use [J]. Ecosystems, 2008,11(4):555-568.

    [17] 惠秀娟,楊 濤,李法云.遼寧省遼河水生態(tài)系統(tǒng)健康評價 [J].應(yīng)用生態(tài)學報, 2011,22(1):181-188.

    [18] 王雪蕾,蔡明勇,鐘部卿,等.遼河流域非點源污染空間特征遙感解析 [J]. 環(huán)境科學, 2013,34(10):3788-3796.

    [19] 殷 燕,吳志旭,張運林,等.新安江水庫夏季CDOM吸收光譜特征及來源分析 [J]. 中國環(huán)境科學, 2014(12):3207-3214.

    [20] Coble P G, Del Castillo C E, Avril B. Distribution and optical properties of CDOM in the Arabian Sea during the 1995 Southwest Monsoon [J]. Deep-Sea Research II, 1998,45(10/11):2195-2223.

    [21] Del Castillo C E, Coble P J, Morel J M, et al. Analysis of the optical properties of the Orinoco River plume by absorption and fluorescence spectroscopy [J]. Marine Chemistry, 1999,66(1):35-51.

    [22] Helms J R, Stubbins A, Ritchie J D, et al. Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and photobleaching of chromophoric dissolved organic matter [J]. Limnology and Oceanography, 2008,53(3):955-969.

    [23] McKnight D M, Boyer E W, Westerhoff P K, et al. Spectrofluorometric characterization of dissolved organic matter for indication of precursor organic material and aromaticity [J]. Limnology and Oceanography, 2001,46(1):38-48.

    [24] Zhang Y L, Yin Y, Feng L Q, et al. Characterizing chromophoric dissolved organic matter in Lake Tianmuhu and its catchment basin using excitation-emission matrix fluorescence and parallel factor analysis [J]. Water Research, 2011,45(16):5110-5122.

    [25] Wada S, Aoki MN, Tsuchiya Y, et al. Quantitative and qualitative analyses of dissolved organic matter released from Ecklonia cava Kjellman, in Oura Bay, Shimoda, Izu Peninsula, Japan [J]. Journal of Experimental Marine Biology and Ecology, 2007,349(2):344-358.

    [26] Ryan D K, Weber J H. Fluorescence quenching titration for determination of complexing capacities and stability constants of fulvic acid [J]. Analytical Chemistry, 1982,54(6):986-990.

    [27] Wu F C, Cai Y C, Evans D, et al. Complexation between Hg (II)and dissolved organic matter in stream waters: an application of fluorescence spectroscopy [J]. Biogeochemistry, 2004,71(3):339-351.

    [28] Bai Y C, Wu F C, Wan G J, et al. Ultraviolet absorbance titration for the determination of conditional stability constants of Hg (II)and dissolved organic matter [J]. Chinese Journal of Geochemistry, 2008,27(1):46-52.

    [29] 韓宇超,郭衛(wèi)東.河口區(qū)有色溶解有機物(CDOM)三維熒光光譜的影響因素 [J]. 環(huán)境科學學報, 2008,28(8):1646-1653.

    [30] Song K S, Wang Z M, Blackwell J, et al. Water quality monitoring using Landsat Themate Mapper data with empirical algorithms in Chagan Lake, China [J]. Journal of Applied Remote Sensing, 2011,5(1):053506-053506-16.

    [31] Song K S, Zang S Y, Zhao Y et al. Spatiotemporal characterization of dissolved carbon for inland waters in semihumid/semiarid region, China [J]. Hydrology and Earth System Sciences Earth, 2013,10(5):6559-6597.

    [32] 姜廣甲,馬榮華,段洪濤.利用 CDOM吸收系數(shù)估算太湖水體表層DOC濃度 [J]. 環(huán)境科學, 2012,33(7):2335-2343.

    [33] Graeber D, Gelbrecht J, Pusch M T, et al. Agriculture has changed the amount and composition of dissolved organic matter in Central European headwater streams [J]. Science of the Total Environment, 2012,438:435-446.

    [34] Huang W, Chen R F. Sources and transformations of chromophoric dissolved organic matter in the Neponset River Watershed [J]. Journal of Geophysical Research: Biogeosciences (2005-2012), 2009,114(G4).

    [35] Larson J H, Frost P C, Zheng Z, et al. Effects of upstream lakes on dissolved organic matter in streams [J]. Limnology and Oceanography, 2007,52(1):60-69.

    [36] 張偉科,楊艷昭,封志明,等.西遼河流域潛在蒸散量時空格局[J]. 地理科學進展, 2010,29(9):1115-1120.

    [37] 傅平青,劉叢強,吳豐昌.溶解有機質(zhì)的三維熒光光譜特征研究[J]. 光譜學與光譜分析, 2005,25(12):2024-2028.

    [38] 郭衛(wèi)東,程遠月,吳 芳.海洋熒光溶解有機物研究進展 [J]. 海洋通報, 2007,26(1):98-106.

    [39] 任保衛(wèi),趙衛(wèi)紅,王江濤.膠州灣圍隔實驗中溶解有機物三維熒光特征 [J]. 環(huán)境科學, 2007,28(4):712-718.

    [40] Baker A. Fluorescence excitation-emission matrix characterization of some sewage-impacted rivers [J]. Environmental Science and Technology, 2001,35(5):948-953.

    [41] Baker A. Fluorescence excitation-emission matrix characterization of river waters impacted by a tissue mill effluent [J]. Environmental Science and Technology, 2002,36(7):1377-1382.

    [42] 邵田田,趙 瑩,宋開山,等.遼河下游CDOM吸收與熒光特性的季節(jié)變化研究 [J]. 環(huán)境科學, 2014,35(10):3755-3763.

    致謝:感謝師弟王銘在數(shù)據(jù)實驗室分析中的幫助以及師弟馬建行在數(shù)據(jù)處理中給予的支持和幫助.

    Comparative characterization of optically active components in waters from the east/west tributaries and the stem of Liaohe River.

    SHAO Tian-tian1,2,3, SONG Kai-shan1*, DING Zhi4, LI Si-jia1, ZHAO Ying1, Zhang Bai1(1.Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Key Laboratory of Wetland Ecology and Environment, Changchun 130102, China;2.Key Research Institute of Yellow River Civilization and Sustainable Development and Collaborative Innovation Center on Yellow River Civilization of Henan Province, Henan University, Kaifeng 475001, China;3.College of Environment and Planning, Henan University, Kaifeng 475001, China;4. Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China).

    China Environmental Science, 2015,35(9):2794~2804

    Comparative analyses about the differences of optical characteristics and composition of riverine waters in the East/West tributaries and the stem of Liaohe River were conducted and the influenced factors impacting on chromophoric dissolved organic matter (CDOM) and dissolve organic carbon (DOC) in different sub-watersheds were examined based on the absorption and fluorescent optical properties. Relationship between CDOM absorption and DOC concentration in the stem of Liaohe River exhibited the best (R = 0.89, P < 0.01) compared with that in the East Liaohe River (R = 0.75) and the West Liaohe River (R = 0.81). Close relationships between two water qualities (total suspended matter [TSM] and total alkalinity) and ratios of CDOM and DOC (CDOM/DOC) were obtained in this investigation. CDOM and DOC concentrations showed great variability in different sub-watersheds. CDOM and DOC concentrations in West Liaohe River were higher than that in East and the stem of Liaohe River. Land use/cover and climate conditions exerted great influences on CDOM and DOC. Significant negative correlations were obtained between forest area and CDOM and DOC,with the R = -0.41 and -0.56respectively. Positive correlations were observed between cropland area and CDOM and DOC (R = 0.40 for CDOM and R =0.32 for DOC). Furthermore, negative correlations were obtained between precipitations and CDOM and DOC, with the R = -0.53 and -0.38 respectively, which demonstrated that dilute effect of rainfall exhibited great influence on CDOM and DOC concentrations. Moreover, relationships between Fn (355) and CDOM absorption in different sub-watersheds showed strong correlations, especially in the West Liaohe River (R = 0.96, P < 0.01). All the samples exhibited fairly strong humic-like fluorophore (fluorophore A and fluorophore C) in the fluorescence excitation-emission matrices (EEMs). Meantime, samples in the main stem of Liaohe River showed strong protein-like fluorophore (fluorophore T). Main sources of CDOM in the studied rivers were terrestrial (allochthonous) origins, which showed the high weight molecular substances based on the S, FI and HIX values. However, molecular weights in West Liaohe River and the main stem of Liaohe River were inclined to be lower than that in the East Liaohe River.

    East Liaohe River;West Liaohe River;chromophoric dissolved organic matter (CDOM);DOC;Land use/cover;Fluorescence

    X522

    A

    1000-6293(2015)09-2794-11

    2015-01-23

    國家重點基礎(chǔ)研究發(fā)展規(guī)劃項目(2013CB430401);國家自然科學基金項目(41471290)

    *責任作者, 研究員, songks@iga.ac.cn

    邵田田(1986-),山東濟寧人,博士研究生,主要研究方向為水體光學特性.發(fā)表論文6篇.

    猜你喜歡
    遼河流域遼河干流
    遼河口
    松花江干流哈爾濱江段封凍
    中國水利(2020年23期)2020-12-24 06:24:40
    傾聽
    遼河(2020年8期)2020-09-02 14:25:25
    遼河文訊
    遼河(2018年5期)2018-11-15 03:39:26
    基于ELM算法的遼河沈陽段水污染評價
    遼寧省遼河流域石油煉制排放限值的制定
    遼河流域水生態(tài)管理指標體系構(gòu)建
    遼河流域水文特性分析
    新疆額敏河干流防洪工程施工分析
    甘肅黑河干流中游濕地現(xiàn)狀及保護對策
    国产美女午夜福利| 美女内射精品一级片tv| 一区二区三区乱码不卡18| 少妇熟女欧美另类| 久久久久久久亚洲中文字幕| 国产av码专区亚洲av| 国产真实伦视频高清在线观看| 久久久久久久久久久丰满| 赤兔流量卡办理| 精品熟女少妇av免费看| 99久久精品热视频| 王馨瑶露胸无遮挡在线观看| 午夜精品国产一区二区电影 | 男人舔奶头视频| 丰满人妻一区二区三区视频av| 成人漫画全彩无遮挡| 黄色视频在线播放观看不卡| 国产探花在线观看一区二区| 国产免费一级a男人的天堂| 精品国产三级普通话版| 国产精品久久久久久av不卡| 蜜桃久久精品国产亚洲av| 国产毛片在线视频| 免费高清在线观看视频在线观看| 国产亚洲最大av| 色婷婷久久久亚洲欧美| 精品久久久噜噜| 成人亚洲欧美一区二区av| 伦理电影大哥的女人| 成年免费大片在线观看| 人人妻人人爽人人添夜夜欢视频 | 一级毛片电影观看| 久久精品国产自在天天线| 简卡轻食公司| 久久精品久久精品一区二区三区| av在线观看视频网站免费| 成人毛片a级毛片在线播放| 我的女老师完整版在线观看| 高清欧美精品videossex| 国产欧美日韩一区二区三区在线 | 在线观看免费高清a一片| 国产精品国产三级国产专区5o| tube8黄色片| .国产精品久久| 精品久久久久久久久亚洲| 在线观看国产h片| 男人爽女人下面视频在线观看| 亚洲国产精品国产精品| 全区人妻精品视频| 久久人人爽人人爽人人片va| 午夜爱爱视频在线播放| 国产成人午夜福利电影在线观看| av在线观看视频网站免费| 成人特级av手机在线观看| 国产片特级美女逼逼视频| 免费黄频网站在线观看国产| 亚洲av日韩在线播放| 免费电影在线观看免费观看| 国产精品嫩草影院av在线观看| 乱码一卡2卡4卡精品| 黄色日韩在线| 国产精品国产三级国产av玫瑰| 精品久久久久久久末码| 伊人久久国产一区二区| 午夜福利在线观看免费完整高清在| 国产成人91sexporn| av卡一久久| 精品99又大又爽又粗少妇毛片| 国产探花极品一区二区| 亚洲欧洲国产日韩| 一本久久精品| 男人添女人高潮全过程视频| 日韩一本色道免费dvd| 国产成人福利小说| 国产精品国产三级国产专区5o| 国产在线男女| 91aial.com中文字幕在线观看| 欧美性猛交╳xxx乱大交人| 麻豆乱淫一区二区| 亚洲欧美日韩东京热| 欧美日韩在线观看h| 91久久精品国产一区二区三区| 国产成人aa在线观看| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 男人狂女人下面高潮的视频| 人人妻人人澡人人爽人人夜夜| 国内揄拍国产精品人妻在线| 国产在视频线精品| 人妻 亚洲 视频| 精品人妻偷拍中文字幕| 噜噜噜噜噜久久久久久91| 久久精品夜色国产| 一个人看视频在线观看www免费| 久久99蜜桃精品久久| 深爱激情五月婷婷| 极品教师在线视频| 美女脱内裤让男人舔精品视频| 美女国产视频在线观看| 亚洲国产精品成人久久小说| 国产伦理片在线播放av一区| 少妇人妻一区二区三区视频| av在线播放精品| 成人亚洲精品av一区二区| 大香蕉久久网| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看| 黄色欧美视频在线观看| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| 欧美高清性xxxxhd video| .国产精品久久| av免费在线看不卡| 最近的中文字幕免费完整| 国产精品国产av在线观看| 久久久午夜欧美精品| 久久久精品欧美日韩精品| 成年免费大片在线观看| 波野结衣二区三区在线| 国产黄a三级三级三级人| 精品一区二区三区视频在线| 色5月婷婷丁香| 插逼视频在线观看| 色吧在线观看| 美女cb高潮喷水在线观看| 一级黄片播放器| 久久久久网色| 国产精品人妻久久久影院| 国产毛片a区久久久久| 日日啪夜夜撸| 久久久精品94久久精品| 国产亚洲精品久久久com| 在线亚洲精品国产二区图片欧美 | 交换朋友夫妻互换小说| 免费av毛片视频| 欧美zozozo另类| 色视频www国产| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 亚洲怡红院男人天堂| 精品少妇久久久久久888优播| 久久久久精品性色| 18禁在线无遮挡免费观看视频| 婷婷色综合www| 青春草国产在线视频| 男女无遮挡免费网站观看| .国产精品久久| 黄片无遮挡物在线观看| 最新中文字幕久久久久| 偷拍熟女少妇极品色| 午夜福利视频1000在线观看| 国产91av在线免费观看| 免费av观看视频| 精品熟女少妇av免费看| 精品久久久久久电影网| 视频中文字幕在线观看| av在线观看视频网站免费| 国产成人精品婷婷| 久久精品久久精品一区二区三区| 亚洲欧美日韩另类电影网站 | 久久久久久久大尺度免费视频| 日韩国内少妇激情av| 国产高潮美女av| 嫩草影院精品99| 亚洲美女搞黄在线观看| 男人舔奶头视频| 久久6这里有精品| 国产综合懂色| 亚洲激情五月婷婷啪啪| 欧美成人一区二区免费高清观看| 色5月婷婷丁香| 久久久成人免费电影| 国产淫片久久久久久久久| 91狼人影院| 少妇人妻久久综合中文| 热re99久久精品国产66热6| 国产高清不卡午夜福利| 日韩强制内射视频| 国产日韩欧美亚洲二区| 街头女战士在线观看网站| 好男人视频免费观看在线| 日韩欧美一区视频在线观看 | 免费高清在线观看视频在线观看| 免费观看的影片在线观看| kizo精华| 中国三级夫妇交换| 国产v大片淫在线免费观看| 日本免费在线观看一区| av国产精品久久久久影院| 99热这里只有是精品50| 又大又黄又爽视频免费| 久久99热6这里只有精品| 亚洲成人av在线免费| 久久久久久国产a免费观看| 午夜精品国产一区二区电影 | 插阴视频在线观看视频| 十八禁网站网址无遮挡 | 国产精品秋霞免费鲁丝片| 乱系列少妇在线播放| 亚洲欧美日韩东京热| 精品午夜福利在线看| 免费看av在线观看网站| 成人毛片60女人毛片免费| 黄色怎么调成土黄色| 亚洲人成网站高清观看| 成人国产麻豆网| 一级二级三级毛片免费看| 91久久精品国产一区二区成人| 99热国产这里只有精品6| 免费电影在线观看免费观看| videossex国产| 在线观看人妻少妇| 肉色欧美久久久久久久蜜桃 | 狂野欧美激情性bbbbbb| 伦精品一区二区三区| 久久热精品热| 亚洲精品色激情综合| 王馨瑶露胸无遮挡在线观看| 99热这里只有精品一区| 岛国毛片在线播放| 国产真实伦视频高清在线观看| 亚洲av电影在线观看一区二区三区 | 狂野欧美白嫩少妇大欣赏| 欧美潮喷喷水| 午夜福利网站1000一区二区三区| 丰满人妻一区二区三区视频av| 插阴视频在线观看视频| 精品人妻一区二区三区麻豆| av福利片在线观看| 久久精品久久久久久噜噜老黄| 久久6这里有精品| 成人二区视频| 尾随美女入室| 婷婷色综合大香蕉| 少妇人妻精品综合一区二区| 你懂的网址亚洲精品在线观看| 亚洲av一区综合| 黄色欧美视频在线观看| 久久久久精品久久久久真实原创| 美女高潮的动态| 国产精品一区二区性色av| 国精品久久久久久国模美| 成年女人在线观看亚洲视频 | 成人一区二区视频在线观看| 亚洲成人中文字幕在线播放| 卡戴珊不雅视频在线播放| 高清在线视频一区二区三区| av在线蜜桃| 综合色av麻豆| 极品教师在线视频| 夫妻性生交免费视频一级片| 大片电影免费在线观看免费| 人妻少妇偷人精品九色| 亚洲最大成人中文| 久久久成人免费电影| 久久久久精品性色| 男女下面进入的视频免费午夜| 一级片'在线观看视频| a级毛色黄片| 欧美最新免费一区二区三区| 亚洲自偷自拍三级| 亚洲人与动物交配视频| 中国美白少妇内射xxxbb| 五月伊人婷婷丁香| 日韩免费高清中文字幕av| 日韩一区二区视频免费看| 欧美日韩在线观看h| 亚洲精品影视一区二区三区av| 午夜精品一区二区三区免费看| 国产有黄有色有爽视频| 国产乱人视频| 一级毛片久久久久久久久女| 国产一级毛片在线| 国产成人一区二区在线| 亚洲欧美精品专区久久| 久久精品综合一区二区三区| 少妇熟女欧美另类| 内射极品少妇av片p| 少妇丰满av| 日日摸夜夜添夜夜爱| 国产久久久一区二区三区| 免费av不卡在线播放| 久久久久久久久大av| 亚洲av免费高清在线观看| 免费人成在线观看视频色| 亚洲激情五月婷婷啪啪| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 男人舔奶头视频| 国产成人免费无遮挡视频| 中文字幕人妻熟人妻熟丝袜美| 免费观看av网站的网址| 我的老师免费观看完整版| 成人国产av品久久久| 午夜视频国产福利| 男女啪啪激烈高潮av片| 视频中文字幕在线观看| 亚洲美女搞黄在线观看| 国产午夜福利久久久久久| 夫妻午夜视频| 天堂俺去俺来也www色官网| 天天一区二区日本电影三级| 国产成人91sexporn| 香蕉精品网在线| av在线观看视频网站免费| 亚洲成人久久爱视频| 夜夜看夜夜爽夜夜摸| 男男h啪啪无遮挡| 在线精品无人区一区二区三 | 久久久久性生活片| av网站免费在线观看视频| 禁无遮挡网站| 亚洲电影在线观看av| 好男人视频免费观看在线| 成人高潮视频无遮挡免费网站| 一级二级三级毛片免费看| 国产在线男女| 91久久精品国产一区二区三区| 联通29元200g的流量卡| 国产av不卡久久| 观看免费一级毛片| 中国三级夫妇交换| 免费黄频网站在线观看国产| 99久久中文字幕三级久久日本| 亚洲精品乱码久久久v下载方式| 亚洲欧美日韩无卡精品| 欧美丝袜亚洲另类| 免费看a级黄色片| 少妇人妻精品综合一区二区| 久久久午夜欧美精品| 国产美女午夜福利| 亚洲av国产av综合av卡| 久久精品国产亚洲网站| 蜜臀久久99精品久久宅男| 亚洲最大成人手机在线| 人妻夜夜爽99麻豆av| 亚洲最大成人手机在线| 两个人的视频大全免费| 国产精品99久久99久久久不卡 | 18禁动态无遮挡网站| 亚洲精品,欧美精品| 韩国av在线不卡| 午夜日本视频在线| 国产精品一二三区在线看| 男人舔奶头视频| 欧美xxxx性猛交bbbb| av网站免费在线观看视频| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 精品熟女少妇av免费看| 男人爽女人下面视频在线观看| 日本-黄色视频高清免费观看| 一本久久精品| 身体一侧抽搐| 精品人妻熟女av久视频| 22中文网久久字幕| 久久人人爽人人片av| 男人舔奶头视频| 免费黄色在线免费观看| av在线观看视频网站免费| 三级国产精品片| 免费看a级黄色片| 韩国高清视频一区二区三区| 欧美日韩视频精品一区| 久久久久久久大尺度免费视频| 日本午夜av视频| 国产精品蜜桃在线观看| 亚洲国产av新网站| 中文在线观看免费www的网站| 亚洲av男天堂| 国产日韩欧美在线精品| 精品久久久久久久久av| 秋霞伦理黄片| 亚洲av电影在线观看一区二区三区 | 中国三级夫妇交换| 青春草视频在线免费观看| 又粗又硬又长又爽又黄的视频| 大又大粗又爽又黄少妇毛片口| 国产国拍精品亚洲av在线观看| 精品视频人人做人人爽| 丰满乱子伦码专区| av在线亚洲专区| 国产高清国产精品国产三级 | 亚洲天堂av无毛| 嘟嘟电影网在线观看| 亚洲欧美日韩无卡精品| 熟妇人妻不卡中文字幕| 国产精品一及| 久久精品熟女亚洲av麻豆精品| 一级a做视频免费观看| 日韩欧美一区视频在线观看 | 又粗又硬又长又爽又黄的视频| 日韩电影二区| 舔av片在线| 亚洲,欧美,日韩| 嫩草影院新地址| 人妻一区二区av| 久久精品久久精品一区二区三区| 免费av不卡在线播放| 国语对白做爰xxxⅹ性视频网站| 成年女人看的毛片在线观看| 国产美女午夜福利| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 亚洲国产av新网站| 亚洲在线观看片| 国产亚洲午夜精品一区二区久久 | 毛片一级片免费看久久久久| 国产av码专区亚洲av| 国产精品熟女久久久久浪| 午夜福利高清视频| 欧美3d第一页| 国产伦精品一区二区三区视频9| 97在线人人人人妻| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| 观看美女的网站| 国产成人精品福利久久| 美女cb高潮喷水在线观看| 成人毛片60女人毛片免费| 51国产日韩欧美| 久久99蜜桃精品久久| 亚洲国产成人一精品久久久| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 嫩草影院精品99| 国产一区二区三区综合在线观看 | tube8黄色片| 亚洲内射少妇av| 18+在线观看网站| av在线天堂中文字幕| 亚洲欧美日韩东京热| 嫩草影院入口| 成人毛片60女人毛片免费| 秋霞在线观看毛片| 男女下面进入的视频免费午夜| 中文字幕人妻熟人妻熟丝袜美| 中文精品一卡2卡3卡4更新| 一级爰片在线观看| 亚洲国产av新网站| 91精品国产九色| h日本视频在线播放| 久久久国产一区二区| 九九爱精品视频在线观看| 久久99精品国语久久久| 在线 av 中文字幕| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 成人美女网站在线观看视频| 男女下面进入的视频免费午夜| 色视频www国产| 日本黄大片高清| 18禁在线播放成人免费| 国产色爽女视频免费观看| av国产免费在线观看| 少妇人妻精品综合一区二区| 成人亚洲精品一区在线观看 | 国产av不卡久久| av在线天堂中文字幕| 丝袜脚勾引网站| 天天一区二区日本电影三级| 舔av片在线| 国产成人午夜福利电影在线观看| 日本av手机在线免费观看| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 日韩av不卡免费在线播放| 美女高潮的动态| 国产精品av视频在线免费观看| 男女那种视频在线观看| 男女国产视频网站| 成人高潮视频无遮挡免费网站| 国语对白做爰xxxⅹ性视频网站| 91久久精品电影网| 亚洲精品456在线播放app| 亚洲综合色惰| 赤兔流量卡办理| 久久6这里有精品| 免费观看av网站的网址| 国产精品国产三级国产专区5o| 男女下面进入的视频免费午夜| 又黄又爽又刺激的免费视频.| 97在线视频观看| 不卡视频在线观看欧美| 亚洲欧美清纯卡通| 午夜福利在线在线| 国产在线一区二区三区精| 亚洲内射少妇av| 十八禁网站网址无遮挡 | 亚洲精品乱久久久久久| 日韩成人伦理影院| 欧美xxxx黑人xx丫x性爽| 中文天堂在线官网| 特大巨黑吊av在线直播| av又黄又爽大尺度在线免费看| 免费av观看视频| 卡戴珊不雅视频在线播放| 成人毛片60女人毛片免费| 又爽又黄无遮挡网站| 亚洲av.av天堂| 另类亚洲欧美激情| 男男h啪啪无遮挡| 91久久精品电影网| 午夜精品国产一区二区电影 | 亚洲国产精品国产精品| 97在线视频观看| 99久久九九国产精品国产免费| 日韩亚洲欧美综合| 国产真实伦视频高清在线观看| 中文在线观看免费www的网站| 综合色av麻豆| 国产高清不卡午夜福利| 99热6这里只有精品| 草草在线视频免费看| av专区在线播放| 女人十人毛片免费观看3o分钟| 国产精品一区二区性色av| 国产一级毛片在线| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久精品免费免费高清| 日日摸夜夜添夜夜爱| 永久免费av网站大全| 特大巨黑吊av在线直播| 亚洲内射少妇av| 九草在线视频观看| 在线免费十八禁| 黄片wwwwww| 男女无遮挡免费网站观看| 我要看日韩黄色一级片| 欧美人与善性xxx| 国产av国产精品国产| 波野结衣二区三区在线| 最后的刺客免费高清国语| 国产成人精品久久久久久| 国产精品一区www在线观看| 蜜臀久久99精品久久宅男| 国产精品伦人一区二区| 亚洲美女视频黄频| av在线app专区| 亚州av有码| 秋霞伦理黄片| 偷拍熟女少妇极品色| 特级一级黄色大片| 亚洲色图综合在线观看| 九草在线视频观看| 六月丁香七月| 欧美最新免费一区二区三区| 69人妻影院| 免费av不卡在线播放| 最近手机中文字幕大全| 寂寞人妻少妇视频99o| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 久久人人爽人人片av| av播播在线观看一区| 女的被弄到高潮叫床怎么办| 最新中文字幕久久久久| 成年av动漫网址| 国产精品久久久久久精品电影| 成人鲁丝片一二三区免费| 国产永久视频网站| 乱系列少妇在线播放| 久久久久久九九精品二区国产| 国模一区二区三区四区视频| 大香蕉久久网| 性色avwww在线观看| 免费观看av网站的网址| 人妻少妇偷人精品九色| 久久久午夜欧美精品| 欧美成人精品欧美一级黄| 免费观看a级毛片全部| 国产在线一区二区三区精| 国产亚洲91精品色在线| 人妻一区二区av| 日韩亚洲欧美综合| 少妇人妻 视频| 97超碰精品成人国产| 男插女下体视频免费在线播放| 欧美日韩视频精品一区| 成人欧美大片| 欧美日韩国产mv在线观看视频 | 一级毛片久久久久久久久女| av在线蜜桃| 别揉我奶头 嗯啊视频| 欧美性感艳星| 国产女主播在线喷水免费视频网站| 亚洲人成网站在线播| 久久ye,这里只有精品| 精品熟女少妇av免费看| 亚洲精品国产av蜜桃| 搡老乐熟女国产| 成年女人看的毛片在线观看| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 成年女人看的毛片在线观看| 久久久久久久国产电影| 久久久精品欧美日韩精品| 制服丝袜香蕉在线| 内地一区二区视频在线| 欧美性感艳星| 成人毛片a级毛片在线播放| 欧美人与善性xxx| 国产高清国产精品国产三级 | 在线观看美女被高潮喷水网站| 欧美日本视频| 最近中文字幕高清免费大全6| 美女脱内裤让男人舔精品视频| 好男人在线观看高清免费视频| 久久99精品国语久久久| 成年av动漫网址| 国产男人的电影天堂91| 中国美白少妇内射xxxbb| 王馨瑶露胸无遮挡在线观看|