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

    典型重污染水體枯水期溶解性有機(jī)質(zhì)光譜特征及其與水質(zhì)的相關(guān)性

    2023-12-13 15:22:35劉潤東顧鑫才陳重軍姜晶
    江蘇農(nóng)業(yè)學(xué)報 2023年7期

    劉潤東 顧鑫才 陳重軍 姜晶

    摘要: 利用紫外可見光譜分析、三維熒光光譜結(jié)合平行因子法研究了佛山市典型重污染水體(汾江河水系和佛山水道水系)溶解性有機(jī)質(zhì)(DOM)的來源與性質(zhì),在此基礎(chǔ)上分析了DOM的光譜特征及其與水質(zhì)的相關(guān)性。結(jié)果表明,兩條水系水質(zhì)均較差,氮、磷含量超標(biāo)嚴(yán)重,部分采樣點(diǎn)位為重污染水體。兩條水系DOM主要包括3個熒光組分,其中C1和C2組分為類腐殖質(zhì),C3組分為類蛋白質(zhì)。2021年1月的汾江河水系有大分子芳香性物質(zhì)輸入,此外其他采樣時間兩條水系DOM均表現(xiàn)為低相對分子質(zhì)量、低芳香性和高親水性特征。DOM熒光組分C1、C2、C3與總磷(TP)含量、溶解性有機(jī)碳(DOC)含量、化學(xué)需氧量(CODCr)均極顯著正相關(guān)(P<0.01),C1、C2與DO含量極顯著負(fù)相關(guān)(P<0.01),C3與DO含量顯著負(fù)相關(guān)(P<0.05),表明在重污染水體中其基本水質(zhì)指標(biāo)受到DOM輸入、微生物活動和光化學(xué)轉(zhuǎn)化等因素的綜合影響。主成分分析結(jié)果表明,兩條水系DOM與總氮(TN)、銨態(tài)氮(NH+4-N)、TP和CODCr具有較強(qiáng)的同源性,均受陸源輸入和內(nèi)源生物作用的影響,外源輸入影響較大。在重污染水體中可用DOM的光譜參數(shù)快速預(yù)測水體富營養(yǎng)化狀態(tài)。本研究結(jié)果可為重污染水環(huán)境治理修復(fù)提供理論依據(jù)。

    關(guān)鍵詞: 重污染水體;溶解性有機(jī)質(zhì);三維熒光光譜;平行因子分析;紫外可見光譜分析

    中圖分類號: X52;X832?? 文獻(xiàn)標(biāo)識碼: A?? 文章編號: 1000-4440(2023)07-1510-09

    Spectral characteristics of dissolved organic matter in typical heavy polluted water during dry season and its relationship with water quality

    LIU Run-dong1,2, GU Xin-cai1,3, CHEN Chong-jun1, JIANG Jing1

    (1.School of Environmental Science and Engineering, Suzhou University of Science and Technology, Suzhou 215009, China;2.National and Local Joint Engineering Laboratory of Municipal Sewage Resource Utilization Technology, Suzhou 215009, China;3.Institute of Agricultural Resources and Environmental Sciences, Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China)

    Abstract: Sources and properties of dissolved organic matter(DOM) in typical heavily polluted river systems (Fenjiang River system and Foshan waterway system) were studied by using ultraviolet-visible spectroscopy, three-dimensional fluorescence spectroscopy and parallel factor analysis. On this basis, the spectral characteristics of DOM and their relationships with water quality were analyzed. The results showed that the water quality of the two water systems was not optimistic, the nitrogen and phosphorus contents exceeded the standards of China, and some sampling points were heavily polluted water bodies. The DOM in the two water systems mainly included three fluorescent components, with C1 and C2 components being humic-like substances and C3 component being protein-like substances. In January 2021, macromolecular aromatic substances were input into Fenjiang River system. In addition, the DOM of the two water systems showed relatively low molecular weight,low aromaticity and high hydrophilicity during other sampling times. The fluorescence components C1, C2 and C3 of DOM were significantly positively correlated with total phosphorus (TP) content, dissolved organic carbon (DOC) content and chemical oxygen demand (CODCr) (P<0.01), C1 and C2 were significantly negatively correlated with DO content (P<0.01), and C3 was significantly negatively correlated with DO content (P<0.05), indicating that the basic water quality indices in heavily polluted water were affected by DOM input, microbial activity and photochemical transformation. The results of principal component analysis showed that the DOM of the two rivers had strong homology with total nitrogen (TN), ammonium nitrogen (NH4+-N), TP and CODCr, which were affected by terrestrial inputs and endogenous biological processes. Exogenous inputs had a greater impact. The spectral parameters of DOM can be used to quickly predict the eutrophication status of water in heavily polluted water. The results of this study can provide a theoretical basis for subsequent efficient water environment treatments.

    Key words: heavily polluted water;dissolved organic matter;three-dimensional fluorescence spectroscopy;parallel factor analysis;ultraviolet-visible spectroscopy

    溶解性有機(jī)質(zhì)(Dissolved organic matter,DOM),是由多糖、腐殖質(zhì)、富里酸和類蛋白質(zhì)等芳香族或脂肪族物質(zhì)組成的復(fù)雜混合物[1],含有豐富的碳、氮、磷等營養(yǎng)元素,能夠參與水環(huán)境中多種物理、化學(xué)和生物過程,在水體生態(tài)系統(tǒng)物質(zhì)循環(huán)、污染物遷移轉(zhuǎn)化和生物降解等方面發(fā)揮著重要作用[2]。目前,國內(nèi)外對海洋、湖泊、河流及濕地等水環(huán)境中DOM已經(jīng)有大量研究,對DOM的來源、性質(zhì)和環(huán)境行為都有系統(tǒng)報道[3-4]。對水體DOM的深入研究結(jié)果表明,對DOM的來源及其組分的異質(zhì)性和復(fù)雜性進(jìn)行全面分析,可以為后續(xù)的水環(huán)境管理及生態(tài)修復(fù)提供重要的理論支撐。城市、農(nóng)村的水體污染程度與民眾的生活息息相關(guān),也越來越受到關(guān)注。尤其是在城市化進(jìn)程加快的背景下,城市、農(nóng)村河流接納超負(fù)荷污染的現(xiàn)象普遍。同時,相比于海洋湖泊,城市和農(nóng)村河流受生活及工農(nóng)業(yè)用水的污染更嚴(yán)重,其水體中DOM組分和性質(zhì)更為復(fù)雜[5-6]。而目前針對城市和農(nóng)村重污染水體中DOM性質(zhì)的研究較少,特別是DOM性質(zhì)與水質(zhì)關(guān)系的研究尚不明確。

    然而DOM性質(zhì)與水質(zhì)的關(guān)系還受水體環(huán)境、DOM來源與組分的影響。如在富營養(yǎng)化的滇池,受到生活污水、藍(lán)藻水華等因素影響,水體有機(jī)碳氮含量高,DOM性質(zhì)與化學(xué)需氧量及總氮的關(guān)聯(lián)性很強(qiáng)。徐兵兵[7]發(fā)現(xiàn),在腐殖化程度較低的河流,類腐殖質(zhì)組分與磷的遷移轉(zhuǎn)化密切相關(guān),有色溶解性有機(jī)質(zhì)在355 nm處的吸收系數(shù)和溶解性總磷濃度能夠較好地表征類腐殖質(zhì)組分。Feng等[8]發(fā)現(xiàn)南四湖相對分子質(zhì)量高的類腐殖質(zhì)、含有色氨酸和酪氨酸的芳香蛋白質(zhì)均與水質(zhì)顯著相關(guān),而相對分子質(zhì)量低的類腐殖質(zhì)與水質(zhì)無顯著相關(guān)。三維熒光光譜結(jié)合平行因子法(EEMs-PARAFAC)以及紫外可見光譜分析(UV-Vis))是研究DOM的快速檢測手段,具有方便快捷、信息量豐富、需要水樣少,靈敏度高等優(yōu)點(diǎn)[9-10]。利用光譜技術(shù)研究城市、農(nóng)村重污染水體DOM性質(zhì)并明確DOM性質(zhì)與水質(zhì)的關(guān)聯(lián),將對未來重污染水體水環(huán)境監(jiān)測與治理起到重要作用。

    改革開放以來,佛山市汾江河和佛山水道兩岸快速城市化,但是部分城中村村民環(huán)保意識薄弱以及工業(yè)納污坑塘和污水處理設(shè)施未同步跟進(jìn),大量污水排入水體,致使河道污染嚴(yán)重,部分河段甚至達(dá)到黑臭水體標(biāo)準(zhǔn)并引起政府和群眾的強(qiáng)烈關(guān)注。因此,本研究以佛山市重污染水體汾江河和佛山水道水系為研究對象,運(yùn)用UV-Vis、EEMs-PARAFAC研究匯入兩條水體的DOM含量和光譜特征,分析DOM來源、性質(zhì)及其與水質(zhì)的關(guān)聯(lián)性,對后續(xù)的污染治理或河道修復(fù)具有指導(dǎo)意義。

    1 材料與方法

    1.1 研究區(qū)域

    研究區(qū)域如圖1所示,汾江河水系除了主要水道汾江河以外,還有6條支涌,分別為:良安截洪溝、羅村涌、良安涌主涌、良安副涌、沙徑支涌、王芝截洪溝羅村段。佛山水道水系主要水道為佛山水道,還有6條支涌,分別為:謝邊涌、水頭涌、鹽步涌、河?xùn)|大涌、東秀大涌、牛肚涌。所選采樣點(diǎn)均為主支流交匯處。參照《水和廢水監(jiān)測分析方法》(第四版)布設(shè)采樣點(diǎn),采樣點(diǎn)位置及具體采樣點(diǎn)信息分別如圖1和表1所示。

    1.2 樣品采集

    本研究的水系全長25.5 km,河寬50~200 m,水深1.5~6.0 m,河床高程變化不大。全年受季風(fēng)影響,雨量充沛,季節(jié)差異性較大,特別是冬季時,雨量明顯減少。10月至次年3月為枯水期,為更大程度上反映人類擾動對水體DOM性質(zhì)的影響,減少雨水對水體DOM性質(zhì)的影響,本研究在枯水期連續(xù)采樣,采樣點(diǎn)位置信息見表1。2020年秋季(11月)和2021年冬季(1月),于河道的中央采集表層水以下0.5 m處的水樣,并用0.45 μm孔徑的針式過濾器(水系)過濾,將所得濾液保存于經(jīng)酸洗的棕色玻璃小瓶中冷藏,帶回實(shí)驗(yàn)室進(jìn)行分析。此外,用經(jīng)過鹽酸酸洗的聚乙烯瓶取500 ml水樣用于總氮 (TN)質(zhì)量濃度、總磷 (TP) 質(zhì)量濃度、化學(xué)需氧量(CODCr)、銨態(tài)氮(NH+4-N)質(zhì)量濃度的測定?,F(xiàn)場測定水體物理指標(biāo),包括透明度(SD)和水溫。

    1.3 樣品分析

    1.3.1 常規(guī)理化指標(biāo)和營養(yǎng)鹽指標(biāo)分析 水樣常規(guī)指標(biāo)包括溶解氧含量 (DO)、pH和氧化還原電位(ORP),在現(xiàn)場用水質(zhì)多參數(shù)測量儀(HI98196)分析,透明度(SD)在現(xiàn)場用賽氏透明盤測定??偟═N)質(zhì)量濃度用堿性-過硫酸鉀消解紫外可見光度法測定,總磷(TP)質(zhì)量濃度用過硫酸鉀消解-鉬酸鹽分光光度計測定,化學(xué)需氧量(CODCr)用重鉻酸鉀法測定,銨態(tài)氮(NH+4-N)質(zhì)量濃度用納氏試劑法測定[11]。

    1.3.2 溶解性有機(jī)質(zhì)分析 水樣經(jīng)過0.45 μm水系針式過濾器過濾后用總有機(jī)碳測定儀(日本島津公司產(chǎn)品,型號TOC5000A)測定溶解性有機(jī)碳(DOC)質(zhì)量濃度。紫外可見吸收光譜通過紫外-可見分光光度計(日本島津公司產(chǎn)品,型號UV-2700)測量,波長和間隔分別設(shè)置為200~800 nm和1 nm,扣除700~800 nm的吸光度均值進(jìn)行光譜結(jié)果校正。使用配備有450 W氙燈和CCD檢測器的熒光分光光譜儀(美國Horiba公司產(chǎn)品,型號Fluorolog-3)測定DOM的三維熒光光譜(EEM)。激發(fā)光譜和發(fā)射光譜的波長分別為250~450 nm(5 nm間隔)和250~600 nm(1 nm間隔),積分時間和狹縫寬度分別設(shè)置為0.2 s和5 nm,以純水進(jìn)行空白校準(zhǔn)。

    1.4 數(shù)據(jù)分析

    紫外可見光譜的吸收系數(shù)a(λ)通過如下公式計算:

    a(λ)=2.303×A(λ)/l

    式中,A(λ)為在λ處經(jīng)過700~800 nm 波長校正后的吸光度,l為比色皿的寬度,單位m。

    以SUVA254表征DOM芳香性,SUVA254越大,芳香性越強(qiáng);以E2∶E3表征相對分子質(zhì)量,E2∶E3隨著DOM相對分子質(zhì)量增大而減小[12]。

    SUVA254=a(254)DOC

    E2∶E3=a(250)a(365)

    通過三維熒光光譜,以腐殖質(zhì)指數(shù)(HIX)、自生源指數(shù)(BIX)[13-14]和熒光指數(shù)(FI)[10]分別表征DOM的腐殖化程度、內(nèi)源物質(zhì)對DOM的相對貢獻(xiàn)和DOM的來源。HIX是激發(fā)波長為254 nm時發(fā)射波長∑(F435→F480)與∑(F300→F345)的比值,BIX是激發(fā)波長為310 nm時,發(fā)射波長380 nm與430 nm的比值,F(xiàn)I是激發(fā)波長為370 nm時,發(fā)射波長470 nm與520 nm的比值。

    采用MATLAB(2012a)中drEEM和N-way工具箱的EEM矩陣方法進(jìn)行PARAFAC分析[15]。為確保模型的準(zhǔn)確性,首先排除掉異常樣品,運(yùn)行并測試模型,采用折半分析和殘差分析來驗(yàn)證PARAFAC模型的穩(wěn)定性和確定DOM的最佳組分?jǐn)?shù)。利用SPSS20.0軟件對數(shù)據(jù)進(jìn)行統(tǒng)計學(xué)分析,采用獨(dú)立樣本t檢驗(yàn)來分析不同水系和同一水系不同月份光譜參數(shù)之間是否存在顯著性差異。以P<0.05和P<0.01表示兩個獨(dú)立樣本組在不同統(tǒng)計學(xué)水平上存在顯著性差異或者變量間顯著性相關(guān)。采用Origin2021軟件繪圖并對熒光組分和熒光指數(shù)進(jìn)行主成分分析(Principal component analysis,PCA)。

    2 結(jié)果與分析

    2.1 水質(zhì)現(xiàn)狀

    所調(diào)查水系的水質(zhì)較差(表2)。兩條水系的SD較低,均接近30 cm。不同采樣時間兩條水系DO質(zhì)量濃度較低,其中11月佛山水系DO質(zhì)量濃度僅為1.87 mg/L,按中國城市黑臭水體污染程度分級標(biāo)準(zhǔn)為輕度黑臭水體。兩條水系pH為6~8,總體上符合地表水Ⅰ~Ⅴ標(biāo)準(zhǔn)[16]。2021年1月兩條水系的氧化還原電位均明顯低于2020年11月的對應(yīng)水系,說明2021年1月兩條水系氧化性較弱。兩條水系的TN、TP質(zhì)量濃度均超地表水Ⅴ類水標(biāo)準(zhǔn),部分點(diǎn)位NH+4-N質(zhì)量濃度甚至超過15 mg/L,呈現(xiàn)重度富營化狀態(tài),NH+4-N質(zhì)量濃度遠(yuǎn)超水體自凈能力,是典型的重度黑臭水體。值得注意的是,兩條水系的CODCr總體處于地表水Ⅲ~Ⅳ水標(biāo)準(zhǔn),這說明兩條水系有較高的氮污染,反硝化作用對生物可利用碳的需求量較大,造成水體碳源不足[17]。多數(shù)采樣點(diǎn)DOC質(zhì)量濃度低于5 mg/L,與天然湖泊、河流的DOC質(zhì)量濃度接近。因此,充分理解重污染水體DOM的來源及性質(zhì),對評估DOM作為反硝化碳源的生物可利用性也有實(shí)際意義。

    2.2 紫外可見光譜參數(shù)特征

    2020年11月佛山水道水系和汾江河水系SUVA254分別為(4.35 ±0.36) L/(mg·m)和(3.99±0.42) L/(mg·m),次年1月佛山水道水系和汾江河水系SUVA254分別為(4.55±1.57) L/(mg·m)和(6.99±3.90) L/(mg·m)(圖2a)??傮w而言,2020年11月兩水系SUVA254波動不大,而次年1月兩水系SUVA254波動較大,這說明1月份輸入到這兩個水系的DOM中芳香性物質(zhì)存在差異,某些芳香性和疏水性強(qiáng)的物質(zhì)進(jìn)入水體。一般認(rèn)為,當(dāng)SUVA254>4 L/(mg·m) 時,DOM顯現(xiàn)出較高芳香性和疏水性;當(dāng)SUVA254<3 L/(mg·m)時,DOM顯現(xiàn)出較低芳香性和高親水性[18]。2021年1月汾江河水系SUVA254均值達(dá)6.99 L/(mg·m),2020年11月汾江河水系SUVA254接近4 L/(mg·m),這表明2021年1月份的汾江河水系存在高芳香性物質(zhì)的輸入。E2∶E3可表征DOM相對分子質(zhì)量變化,該比值隨DOM相對分子質(zhì)量增大而減小[19]。佛山水道水系2020年11月與次年1月E2∶E3無顯著差異,而汾江河水系2020年11月E2∶E3為6.13±0.65,顯著高于次年1月的E2∶E3(P<0.01),同時汾江河水系次年1月的E2∶E3也顯著低于佛山水道水系2020年11月和次年1月的E2∶E3(P<0.05)(圖2b),這說明2021年1月汾江河水系有相對分子質(zhì)量較大的物質(zhì)輸入,高芳香性的DOM相對分子質(zhì)量大,這與SUVA254的結(jié)果一致。

    2.3 三維熒光光譜及光譜參數(shù)特征

    將PARAFAC 分離的熒光組分與OpenFluor數(shù)據(jù)庫已公開的數(shù)據(jù)[20]進(jìn)行對比,設(shè)置激發(fā)與發(fā)射光譜相似得分參數(shù)為0.95(圖3)。C1組分(Ex/Em=325 nm/410 nm),為高芳香性的大分子類腐殖質(zhì)物質(zhì),與人為活動有關(guān)[5,20-21];C2組分(Ex/Em=380 nm/470 nm)為外源類腐殖酸物質(zhì)[5,22];C3組分(Ex/Em=280 nm/335 nm)為類蛋白熒光物質(zhì),在自然界中該組分能被微生物輕易的利用降解[22-23]。C1和C2組分的相關(guān)性較高,R2達(dá)0.98,而C1和C3,C2和C3的相關(guān)性則較低,R2均只有0.61,這說明C1和C2具有較好的同源性,而C1和C3,C2和C3的來源存在差異(圖4)。C3組分主要是內(nèi)源代謝貢獻(xiàn),一般為水體微生物代謝產(chǎn)物,本研究生活污水中分解的蛋白質(zhì)也會對C3有較大貢獻(xiàn)。C1和C3的Fmax值相當(dāng),C2的Fmax略低(圖5)??傮w而言,不同月份不同河道間均值差異性不顯著,但是部分采樣點(diǎn)位差異明顯,這跟污染負(fù)荷輸入密切相關(guān)。

    熒光指數(shù)FI用于表征DOM中腐殖質(zhì)組分的來源。FI越高表明DOM的生物來源越強(qiáng)。當(dāng)FI>1.9時,微生物代謝產(chǎn)物為DOM的主要來源,F(xiàn)I<1.4表明DOM主要受陸源輸入的影響。當(dāng)1.4<FI<1.9時,DOM成分同時受生物源和陸源的影響[10]。兩條水系FI為1.4~1.9,表明水系中DOM同時受陸源和生物源的共同影響(圖6)。HIX用于表征DOM的腐殖化程度,并且其腐殖化程度隨HIX的升高而升高[20,24]。兩條水系HIX均小于4,說明腐殖化程度較低。此外,2020年11月佛山水道水系以及2020年11月和次年1月汾江河水系BIX (>1)較高,表明水體存在較強(qiáng)的自生源特征及較高的新近DOM釋放的可能[13,20]。2020年11月和次年1月汾江河水系BIX分別為(1.06±0.05)、(1.02±0.05),生物活性較高。而次年1月佛山水道水系BIX顯著低于同時期的汾江河水系,也顯著低于佛山水道水系2020年11月的BIX(P<0.05)。結(jié)合2021年1月佛山水道水系FI顯著低于同時期的汾江河水系(P<0.05),而HIX值顯著高于2020年11月佛山水道水系(P<0.05),說明2021年1月佛山水道水系陸源輸入的DOM主要與地表徑流等高腐殖化廢水的排入有關(guān),而2021年1月汾江河水系的DOM則很有可能與生物可利用性高的生活污水有關(guān)。

    2.4 DOM熒光組分及光譜參數(shù)與水質(zhì)參數(shù)的關(guān)系

    研究DOM熒光組分與水質(zhì)參數(shù)的相關(guān)關(guān)系,有助于深刻理解城市重污染水體DOM的性質(zhì),有助于利用DOM性質(zhì)快速預(yù)測水體水質(zhì)狀況。經(jīng)過SPSS相關(guān)性分析和指標(biāo)遴選,DOM組分之間及其與水質(zhì)參數(shù)的關(guān)系如表3所示,DOC與C1、C2和C3呈極顯著正相關(guān) (P<0.01),表明DOM是水體中重要的溶解有機(jī)碳庫,其與碳在環(huán)境中的遷移轉(zhuǎn)化過程密不可分[25]。C1和C3與SUVA254呈顯著負(fù)相關(guān)(P<0.05), C2與SUVA254呈極顯著負(fù)相關(guān)(P<0.01),這說明水體中DOM呈現(xiàn)出低芳香性特征。HIX與E2∶E3呈顯著正相關(guān)(P<0.05),這說明在DOM腐殖化程度較高的采樣點(diǎn)位,其相對分子質(zhì)量較小。這與先前的研究中較高腐殖化程度DOM其芳香性和相對分子質(zhì)量往往較大的結(jié)論相反,造成這種差異的原因可能是由于樣品采集時期屬于枯水期(11月和次年1月),地表徑流和土壤淋溶對水體中DOM分子的影響較弱,而光化學(xué)轉(zhuǎn)化則顯著影響了DOM的分子結(jié)構(gòu),使其芳香性降低并由大分子物質(zhì)轉(zhuǎn)化為小分子物質(zhì)[26]。

    從熒光組分與水質(zhì)參數(shù)之間的關(guān)系看,C1、C2和C3與DO的質(zhì)量濃度呈顯著負(fù)相關(guān),其原因可能是隨著外源碳輸入的增強(qiáng),水體中碳源豐富,好氧微生物活性增強(qiáng),對碳進(jìn)行分解,同時產(chǎn)生大量類蛋白物質(zhì)從而導(dǎo)致DO質(zhì)量濃度降低[27]。C1、C2和C3與TP質(zhì)量濃度、CODCr均呈極顯著正相關(guān),這說明在城市重污染水體中,可以用DOM的組分快速預(yù)測水體營養(yǎng)狀態(tài)[28]。熒光組分(C1、C2和C3)與TP質(zhì)量濃度、CODCr極顯著相關(guān),說明DOM對TP質(zhì)量濃度和CODCr的貢獻(xiàn)較大,但值得注意的是僅C3組分與TN質(zhì)量濃度呈顯著正相關(guān),這可能與類蛋白物質(zhì)(C3)中蛋白質(zhì)、氨基酸和酰胺類等有機(jī)氮組分含量有關(guān),而NH+4-N在TN中占較大比重(表2),也可能是TN質(zhì)量濃度與C1、C2無顯著相關(guān)的原因。

    對汾江河水系水體的光譜參數(shù)與水質(zhì)參數(shù)進(jìn)行主成分分析,結(jié)果表明第一主成分和第二主成分共解釋了所選指標(biāo)59.5%方差變化,其中第一主成分與第二主成分對所選指標(biāo)的解釋率分別為42.6%和16.6%(圖7)。第一主成分軸上總氮(TN)質(zhì)量濃度、銨態(tài)氮(NH+4-N)質(zhì)量濃度、總磷(TP)質(zhì)量濃度和化學(xué)需氧量(CODCr)具有較強(qiáng)的正相關(guān)性,表明第一主成分正方向與水體富營養(yǎng)化有關(guān),同時C1、C2和C3與第一主成分荷載呈正相關(guān),表明C1、C2和C3與總氮(TN)質(zhì)量濃度、銨態(tài)氮(NH+4-N)質(zhì)量濃度、總磷(TP)質(zhì)量濃度和化學(xué)需氧量(CODCr)具有較強(qiáng)的同源性[29]。重污染水體中總氮、銨態(tài)氮、總磷和化學(xué)需氧量(CODCr)一般為外源輸入,與周邊生活污水以及工業(yè)廢水的輸入有關(guān)[30],因此水體中DOM主要考慮為外源輸入。C3、總氮(TN)和自身源指數(shù)(BIX)與第二主成分組分呈正相關(guān)性,表明第二主成分正方向與自生源特征有關(guān),較高含量的類蛋白質(zhì)組分C3和總氮(TN)為微生物的生長代謝提供了碳源,同時也再次驗(yàn)證了C3中有機(jī)氮對總氮(TN)的貢獻(xiàn)。

    3 結(jié)論

    本研究結(jié)論如下:(1)汾江河水系和佛山水道水系DOM主要包括3個熒光組分,分別為類腐殖質(zhì)C1(325 nm/410 nm)、類腐殖質(zhì)C2(380 nm/470 nm)和類蛋白質(zhì)組分C3(280 nm/335 nm)。除次年1月汾江河水系有相對分子質(zhì)量較大的芳香性物質(zhì)輸入,其他采樣時間兩條水系DOM均表現(xiàn)為低相對分子質(zhì)量、低芳香性和高親水性特征。兩條水系DOM腐殖化程度低。(2)水質(zhì)參數(shù)和DOM組成受到面源輸入、微生物活動和光化學(xué)轉(zhuǎn)化等因素的綜合影響。(3)DOM熒光組分C1、C2、C3與水質(zhì)參數(shù)TP、DOC質(zhì)量濃度、CODCr均極顯著相關(guān)(P<0.01),C1、C2與DO質(zhì)量濃度極顯著相關(guān)(P<0.01),C3與TN質(zhì)量濃度顯著相關(guān)(P<0.05)。主成分分析結(jié)果表明,DOM熒光組分與TN、NH+4-N、TP和CODCr具有較強(qiáng)的同源性,受外源輸入影響較大,即與周邊生活污水以及工業(yè)廢水的輸入有關(guān)。因此,在重污染水體中可用DOM的光譜參數(shù)快速預(yù)測水體富營養(yǎng)化狀態(tài)。

    參考文獻(xiàn):

    [1] HUANG C,ZHANG L,LI Y,et al. Carbon and nitrogen burial in a plateau lake during eutrophication and phytoplankton blooms [J]. Science of the Total Environment,2018,616/617:296-304.

    [2] REN H,F(xiàn)AN T,YAO X,et al. Investigation of the variations in dissolved organic matter properties and complexations with two typical heavy metals under the influence of biodegradation:A survey of an entire lake [J]. Science of the Total Environment,2022,806:150485.

    [3] ASMALA E,MASSICOTTE P,CARSTENSEN J. Identification of dissolved organic matter size components in freshwater and marine environments [J]. Limnology and Oceanography,2021,66(4):1381-1393.

    [4] 李帥東,張明禮,楊 浩,等. 昆明松華壩庫區(qū)表層土壤溶解性有機(jī)質(zhì)(DOM)的光譜特性 [J]. 光譜學(xué)與光譜分析,2017,37(4):1183-1188.

    [5] REN W,WU X,GE X,et al. Characteristics of dissolved organic matter in lakes with different eutrophic levels in southeastern Hubei Province,China [J]. Journal of Oceanology and Limnology,2021,39(4):1256-1276.

    [6] 何 杰,李學(xué)艷,林 欣,等. 光譜特征法辨識不同污染景觀河道中溶解性有機(jī)物的組分與來源 [J]. 環(huán)境科學(xué)學(xué)報,2021,41(3):1000-1010.

    [7] 徐兵兵. 南苕溪溶解性有機(jī)質(zhì)光譜特性及其與水質(zhì)參數(shù)的相關(guān)性研究 [J]. 環(huán)境污染與防治,2021,43(9):1144-1150,1158.

    [8] FENG L,ZHANG J,F(xiàn)AN J,et al. Tracing dissolved organic matter in inflowing rivers of Nansi Lake as a storage reservoir:Implications for water-quality control [J]. Chemosphere,2022,286:131624.

    [9] CHEN W,HABIBUL N,LIU X Y,et al. FTIR and synchronous fluorescence heterospectral two-dimensional correlation analyses on the binding characteristics of copper onto dissolved organic matter [J]. Environmental Science & Technology,2015,49(4):2052-2058.

    [10]LI Y,ZHANG Y,LI Z,et al. Characterization of colored dissolved organic matter in the northeastern South China Sea using EEMs-PARAFAC and absorption spectroscopy[J]. Journal of Sea Research,2022,180:102159.

    [11]陳丙法. 湖岸帶藍(lán)藻聚積污染物特性變化及控制研究 [D]. 蘇州:蘇州科技學(xué)院,2015.

    [12]NI M,LI S. Biodegradability of riverine dissolved organic carbon in a dry-hot valley region:Initial trophic controls and variations in chemical composition [J]. Journal of Hydrology,2019,574:430-435.

    [13]HUGUET A,VACHER L,RELEXANS S,et al. Properties of fluorescent dissolved organic matter in the Gironde Estuary [J]. Organic Geochemistry,2009,40(6):706-719.

    [14]HUR J,PARK M H,SCHLAUTMAN M A. Microbial transformation of dissolved leaf litter organic matter and its effects on selected organic matter operational descriptors [J]. Environmental Science & Technology,2009,43(7):2315-2321.

    [15]ISHII S K,BOYER T H. Behavior of reoccurring PARAFAC components in fluorescent dissolved organic matter in natural and engineered systems:a critical review [J]. Environmental Science & Technology,2012,46(4):2006-2017.

    [16]國家環(huán)境保護(hù)總局. 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn):GB 3838-2002[S]. 北京:中國標(biāo)準(zhǔn)出版社,2002.

    [17]向 衡,韓 蕓,劉 琳,等. 用于河道水反硝化脫氮補(bǔ)充碳源選擇研究 [J]. 水處理技術(shù),2013,39(5):64-68.

    [18]MATILAINEN A,GJESSING E T,LAHTINEN T,et al. An overview of the methods used in the characterisation of natural organic matter (NOM) in relation to drinking water treatment [J]. Chemosphere,2011,83:1431-1442.

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

    [20]陳丙法. 水生植物生物炭對水體重金屬作用機(jī)制研究 [D]. 北京:中國科學(xué)院大學(xué),2020.

    [21]LIU C,DU Y,YIN H,et al. Exchanges of nitrogen and phosphorus across the sediment-water interface influenced by the external suspended particulate matter and the residual matter after dredging [J]. Environmental Pollution,2019,246:207-216.

    [22]LEE S A,KIM T H,KIM G. Tracing terrestrial versus marine sources of dissolved organic carbon in a coastal bay using stable carbon isotopes [J]. Biogeosciences,2020,17(1):135-144.

    [23]GAO Z,GUEGUEN C. Size distribution of absorbing and fluorescing DOM in Beaufort Sea,Canada Basin [J]. Deep Sea Research Part I:Oceanographic Research Papers,2017,121:30-37.

    [24]SINGH S,D'SA E J,SWENSON E M. Chromophoric dissolved organic matter(CDOM) variability in Barataria Basin using excitation-emission matrix (EEM) fluorescence and parallel factor analysis(PARAFAC) [J]. Science of the Total Environment,2010,408(16):3211-3222.

    [25]許金鑫. 城市河流溶解性有機(jī)質(zhì)的光譜特征研究 [D]. 上海:上海師范大學(xué),2020.

    [26]XU H,LI Y,LIU J,et al. Photogeneration and steady-state concentration of hydroxyl radical in river and lake waters along middle-lower Yangtze region,China[J]. Water Research,2020,176:115774.

    [27]蔡文良,許曉毅,羅固源,等. 長江重慶段溶解性有機(jī)物的熒光特性分析[J]. 環(huán)境化學(xué),2012,31(7):1003-1008.

    [28]ZHANG Y,ZHOU Y,SHI K,et al. Optical properties and composition changes in chromophoric dissolved organic matter along trophic gradients:Implications for monitoring and assessing lake eutrophication [J]. Water Research,2018,131:255-263.

    [29]李曉潔,高紅杰,郭冀峰,等. 三維熒光與平行因子研究黑臭河流DOM [J]. 中國環(huán)境科學(xué),2018,38(1):311-319.

    [30]程慶霖,鄭丙輝,王圣瑞,等. 滇池水體有色溶解性有機(jī)質(zhì)(CDOM)三維熒光光譜特征[J]. 光譜學(xué)與光譜分析,2014,34(3):698-703.

    (責(zé)任編輯:成紓寒)

    久久久久久久久久成人| 亚洲av福利一区| 超碰av人人做人人爽久久| 麻豆成人午夜福利视频| 午夜精品一区二区三区免费看| 国产黄片美女视频| 一区二区三区精品91| 一级毛片aaaaaa免费看小| 国产成人精品福利久久| 精华霜和精华液先用哪个| 国产精品国产三级国产av玫瑰| 一二三四中文在线观看免费高清| 国产爽快片一区二区三区| 麻豆成人av视频| freevideosex欧美| 一级片'在线观看视频| 欧美xxxx黑人xx丫x性爽| 久久久久久久国产电影| 亚洲人成网站在线观看播放| 2022亚洲国产成人精品| 综合色av麻豆| 大香蕉97超碰在线| 日日啪夜夜爽| 波野结衣二区三区在线| 亚洲怡红院男人天堂| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 免费黄频网站在线观看国产| 日本熟妇午夜| 26uuu在线亚洲综合色| 嫩草影院新地址| 欧美丝袜亚洲另类| 麻豆成人av视频| 国产极品天堂在线| 视频中文字幕在线观看| 免费看av在线观看网站| 成年版毛片免费区| 久久久成人免费电影| 91久久精品国产一区二区成人| 香蕉精品网在线| 在线精品无人区一区二区三 | 欧美亚洲 丝袜 人妻 在线| 深夜a级毛片| 2021天堂中文幕一二区在线观| 久久国产乱子免费精品| 国产成人免费观看mmmm| 男女边吃奶边做爰视频| 国产黄片视频在线免费观看| 亚洲国产最新在线播放| 日本一二三区视频观看| 欧美xxxx性猛交bbbb| 久久久久久久大尺度免费视频| 久久久久久久久久久丰满| 夜夜爽夜夜爽视频| 三级经典国产精品| 国产一区亚洲一区在线观看| 色婷婷久久久亚洲欧美| 午夜免费男女啪啪视频观看| 在线观看av片永久免费下载| 偷拍熟女少妇极品色| 99热这里只有是精品在线观看| 97超碰精品成人国产| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 国产综合精华液| 国产精品嫩草影院av在线观看| 精品一区二区三区视频在线| 久久精品综合一区二区三区| 久久国内精品自在自线图片| 97在线人人人人妻| 网址你懂的国产日韩在线| 亚洲一区二区三区欧美精品 | 国产视频内射| 乱码一卡2卡4卡精品| 在线精品无人区一区二区三 | av国产免费在线观看| 一级毛片久久久久久久久女| 国语对白做爰xxxⅹ性视频网站| 干丝袜人妻中文字幕| 欧美潮喷喷水| 网址你懂的国产日韩在线| av在线观看视频网站免费| 我的女老师完整版在线观看| 一区二区三区四区激情视频| 水蜜桃什么品种好| 99热这里只有是精品50| 色综合色国产| 亚洲熟女精品中文字幕| 久久影院123| 婷婷色综合www| 国产精品久久久久久精品古装| 日本一本二区三区精品| 亚洲国产最新在线播放| 精品久久久精品久久久| 国产精品一二三区在线看| 国产大屁股一区二区在线视频| 亚洲无线观看免费| 中文天堂在线官网| 在线观看三级黄色| 欧美日韩国产mv在线观看视频 | 亚洲精品456在线播放app| 可以在线观看毛片的网站| 成年版毛片免费区| 国产女主播在线喷水免费视频网站| 国产 一区精品| 亚洲国产欧美在线一区| 亚洲aⅴ乱码一区二区在线播放| 亚洲四区av| 18禁动态无遮挡网站| 久久精品国产亚洲av涩爱| 日韩av不卡免费在线播放| 一级毛片 在线播放| 国产av国产精品国产| 水蜜桃什么品种好| 亚洲不卡免费看| 精品久久久久久久人妻蜜臀av| 欧美最新免费一区二区三区| 男女啪啪激烈高潮av片| 久久久成人免费电影| 嫩草影院精品99| 性色avwww在线观看| 欧美丝袜亚洲另类| videos熟女内射| 午夜视频国产福利| 永久免费av网站大全| 国产av不卡久久| 亚洲精品第二区| 麻豆成人av视频| 日本wwww免费看| 美女cb高潮喷水在线观看| 亚洲三级黄色毛片| eeuss影院久久| 亚洲国产最新在线播放| kizo精华| 国产成人免费观看mmmm| 青春草国产在线视频| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 国产精品一区www在线观看| 97超碰精品成人国产| 26uuu在线亚洲综合色| 国产在视频线精品| 小蜜桃在线观看免费完整版高清| 国产伦精品一区二区三区视频9| 亚洲精品国产av成人精品| 99久久精品一区二区三区| 自拍偷自拍亚洲精品老妇| 久久97久久精品| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 青青草视频在线视频观看| 成年av动漫网址| 欧美三级亚洲精品| 日韩成人av中文字幕在线观看| 亚洲精品乱码久久久v下载方式| 色网站视频免费| 国产亚洲91精品色在线| 国产精品伦人一区二区| 午夜福利在线观看免费完整高清在| 激情五月婷婷亚洲| 97在线人人人人妻| 日韩av在线免费看完整版不卡| 久久久久久久精品精品| 伊人久久国产一区二区| 亚洲av中文av极速乱| 亚洲色图综合在线观看| 人妻 亚洲 视频| 国产 一区精品| 一本久久精品| 欧美人与善性xxx| 国产成人免费观看mmmm| 热re99久久精品国产66热6| 黄色日韩在线| 综合色av麻豆| 老师上课跳d突然被开到最大视频| 午夜福利高清视频| 性插视频无遮挡在线免费观看| 亚洲av成人精品一区久久| 一区二区三区精品91| 欧美成人一区二区免费高清观看| av在线亚洲专区| 免费av毛片视频| 国产 精品1| 日韩,欧美,国产一区二区三区| 爱豆传媒免费全集在线观看| 日本黄色片子视频| 麻豆久久精品国产亚洲av| 久久人人爽av亚洲精品天堂 | 免费看光身美女| 97人妻精品一区二区三区麻豆| 亚洲精品久久午夜乱码| 精品久久久噜噜| 亚洲自拍偷在线| 插阴视频在线观看视频| 一个人看视频在线观看www免费| 国产精品一及| 18禁动态无遮挡网站| 三级国产精品欧美在线观看| 亚洲aⅴ乱码一区二区在线播放| 青春草视频在线免费观看| 亚洲精品日韩av片在线观看| 亚洲欧美一区二区三区国产| 国产黄片美女视频| 超碰av人人做人人爽久久| 丝袜喷水一区| 99热国产这里只有精品6| 久久久久精品性色| 美女脱内裤让男人舔精品视频| 久久影院123| videos熟女内射| av在线观看视频网站免费| 18禁裸乳无遮挡免费网站照片| kizo精华| 日韩av免费高清视频| 国产精品蜜桃在线观看| 亚洲精品日韩在线中文字幕| 好男人在线观看高清免费视频| 亚洲在久久综合| 日韩三级伦理在线观看| av线在线观看网站| 欧美日韩视频精品一区| 狂野欧美激情性bbbbbb| 夜夜看夜夜爽夜夜摸| 亚洲av一区综合| 国产精品麻豆人妻色哟哟久久| 国产黄频视频在线观看| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区黑人 | a级毛片免费高清观看在线播放| 亚洲图色成人| av免费观看日本| av女优亚洲男人天堂| 精品酒店卫生间| 97在线视频观看| 亚洲精品乱码久久久v下载方式| 男人舔奶头视频| 欧美日韩国产mv在线观看视频 | 国产精品无大码| 日韩亚洲欧美综合| 成人亚洲欧美一区二区av| 一级二级三级毛片免费看| 真实男女啪啪啪动态图| 亚洲怡红院男人天堂| 麻豆精品久久久久久蜜桃| 精品国产露脸久久av麻豆| 成人漫画全彩无遮挡| 一级爰片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲一区二区精品| 91午夜精品亚洲一区二区三区| 成人无遮挡网站| 高清视频免费观看一区二区| 成人特级av手机在线观看| 国产日韩欧美亚洲二区| 久久久久精品久久久久真实原创| 成年免费大片在线观看| 男女啪啪激烈高潮av片| 麻豆成人午夜福利视频| 天堂中文最新版在线下载 | 亚洲无线观看免费| 国产视频首页在线观看| 人妻 亚洲 视频| 中国三级夫妇交换| 97在线人人人人妻| 九九久久精品国产亚洲av麻豆| 国产精品一二三区在线看| 国产免费一级a男人的天堂| 嫩草影院入口| www.av在线官网国产| 91久久精品国产一区二区三区| 国产色爽女视频免费观看| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站 | 亚洲成人精品中文字幕电影| 亚洲精品第二区| 寂寞人妻少妇视频99o| 欧美成人a在线观看| 欧美日本视频| 精品一区在线观看国产| 免费人成在线观看视频色| 国产伦精品一区二区三区视频9| 啦啦啦啦在线视频资源| 欧美亚洲 丝袜 人妻 在线| 久久久久网色| 国产免费一级a男人的天堂| 亚洲高清免费不卡视频| 大香蕉久久网| 一区二区三区四区激情视频| 日韩在线高清观看一区二区三区| 亚洲精品国产成人久久av| 国产成人a区在线观看| 亚洲国产av新网站| 内地一区二区视频在线| 亚洲自拍偷在线| 少妇人妻精品综合一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲国产最新在线播放| 大陆偷拍与自拍| 黄色怎么调成土黄色| 18+在线观看网站| 精品午夜福利在线看| 97在线人人人人妻| 久久久久久久久大av| 国产91av在线免费观看| 免费播放大片免费观看视频在线观看| 18禁动态无遮挡网站| 2021天堂中文幕一二区在线观| 国产精品99久久99久久久不卡 | 嘟嘟电影网在线观看| 亚洲精品久久久久久婷婷小说| 久久久久精品性色| 国产爱豆传媒在线观看| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 美女内射精品一级片tv| av福利片在线观看| 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 在线观看美女被高潮喷水网站| 少妇的逼水好多| xxx大片免费视频| 又爽又黄a免费视频| 亚洲欧洲国产日韩| 亚洲不卡免费看| 99热国产这里只有精品6| 女人被狂操c到高潮| 久久久午夜欧美精品| 精品人妻熟女av久视频| 亚洲精品乱码久久久久久按摩| 亚洲国产色片| 国产乱人视频| 国产老妇伦熟女老妇高清| 插逼视频在线观看| 日韩免费高清中文字幕av| .国产精品久久| 91精品伊人久久大香线蕉| 国产成人freesex在线| 99九九线精品视频在线观看视频| 国产成人一区二区在线| 精品少妇久久久久久888优播| 人妻一区二区av| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 亚洲av福利一区| 亚洲av成人精品一区久久| 色吧在线观看| 国产一区有黄有色的免费视频| 精品一区二区免费观看| 日日啪夜夜爽| 97超视频在线观看视频| 99热这里只有是精品50| av免费在线看不卡| 国产欧美日韩一区二区三区在线 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成年女人在线观看亚洲视频 | 色视频www国产| 免费观看av网站的网址| 青青草视频在线视频观看| 天堂网av新在线| 亚洲精品成人久久久久久| 亚洲欧美中文字幕日韩二区| 黄色一级大片看看| 男女边吃奶边做爰视频| 欧美成人午夜免费资源| 男人舔奶头视频| 最后的刺客免费高清国语| 欧美+日韩+精品| 大码成人一级视频| 免费黄网站久久成人精品| 日日啪夜夜撸| 春色校园在线视频观看| 欧美xxⅹ黑人| 国产美女午夜福利| 自拍偷自拍亚洲精品老妇| 欧美丝袜亚洲另类| av免费观看日本| 中文字幕制服av| 欧美bdsm另类| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 一级爰片在线观看| 国产一区二区亚洲精品在线观看| 免费av观看视频| 亚洲av电影在线观看一区二区三区 | 欧美+日韩+精品| 国产亚洲午夜精品一区二区久久 | 亚洲欧美日韩卡通动漫| 国产黄片美女视频| 99久久精品国产国产毛片| 国内精品美女久久久久久| 国内精品宾馆在线| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 汤姆久久久久久久影院中文字幕| 免费看日本二区| 99热国产这里只有精品6| 国产久久久一区二区三区| 99热网站在线观看| 亚洲三级黄色毛片| 我要看日韩黄色一级片| 国产成人freesex在线| 欧美变态另类bdsm刘玥| 国产视频首页在线观看| 免费看光身美女| 在线观看免费高清a一片| 久久久久久久久久人人人人人人| 婷婷色综合www| 成人毛片60女人毛片免费| 简卡轻食公司| 街头女战士在线观看网站| 欧美 日韩 精品 国产| 男女下面进入的视频免费午夜| 直男gayav资源| av免费在线看不卡| 午夜免费观看性视频| 91久久精品国产一区二区三区| 亚洲无线观看免费| 久久ye,这里只有精品| 国产淫语在线视频| www.色视频.com| 亚洲成人久久爱视频| 熟女电影av网| 日韩成人av中文字幕在线观看| 国产伦精品一区二区三区四那| 亚洲国产精品专区欧美| 国产精品久久久久久久久免| 2022亚洲国产成人精品| 麻豆成人午夜福利视频| 国产精品久久久久久久电影| 男男h啪啪无遮挡| 内射极品少妇av片p| 国产精品爽爽va在线观看网站| 成人美女网站在线观看视频| 日本wwww免费看| 蜜桃亚洲精品一区二区三区| 午夜日本视频在线| 日韩亚洲欧美综合| 亚洲成人av在线免费| 色播亚洲综合网| 国产成人freesex在线| 国产精品精品国产色婷婷| 美女被艹到高潮喷水动态| 免费不卡的大黄色大毛片视频在线观看| 丝瓜视频免费看黄片| 最新中文字幕久久久久| 久热久热在线精品观看| 汤姆久久久久久久影院中文字幕| 亚洲国产精品成人久久小说| 一个人看的www免费观看视频| 天美传媒精品一区二区| 秋霞伦理黄片| 日韩三级伦理在线观看| 亚洲图色成人| 精华霜和精华液先用哪个| 久久精品夜色国产| 两个人的视频大全免费| 美女cb高潮喷水在线观看| 黄色日韩在线| 久久97久久精品| 搞女人的毛片| 中国国产av一级| 小蜜桃在线观看免费完整版高清| 久久人人爽人人爽人人片va| kizo精华| a级毛片免费高清观看在线播放| 久久精品国产亚洲av涩爱| 王馨瑶露胸无遮挡在线观看| 久久久久久九九精品二区国产| 热re99久久精品国产66热6| 草草在线视频免费看| 黄色欧美视频在线观看| 免费大片18禁| av福利片在线观看| 国产精品99久久99久久久不卡 | 高清在线视频一区二区三区| 久久久久久国产a免费观看| 伦精品一区二区三区| 18禁在线无遮挡免费观看视频| 身体一侧抽搐| 日本黄色片子视频| 高清日韩中文字幕在线| 人人妻人人看人人澡| 亚洲aⅴ乱码一区二区在线播放| 成年女人看的毛片在线观看| 青春草国产在线视频| 亚洲美女视频黄频| 国产伦精品一区二区三区四那| 精品少妇久久久久久888优播| 色网站视频免费| 国产美女午夜福利| 久久久久精品久久久久真实原创| 青春草视频在线免费观看| 亚洲av一区综合| 人体艺术视频欧美日本| 99热全是精品| 免费av毛片视频| 成人午夜精彩视频在线观看| 色播亚洲综合网| 成人毛片60女人毛片免费| 亚洲成人一二三区av| 日韩大片免费观看网站| 国产伦理片在线播放av一区| 97热精品久久久久久| 久久久久久九九精品二区国产| 蜜桃亚洲精品一区二区三区| 亚洲自偷自拍三级| 亚洲在久久综合| 禁无遮挡网站| 亚洲欧美精品专区久久| 免费黄网站久久成人精品| 亚洲内射少妇av| 人人妻人人爽人人添夜夜欢视频 | av又黄又爽大尺度在线免费看| 尾随美女入室| 国产在视频线精品| 国产精品偷伦视频观看了| 直男gayav资源| a级毛片免费高清观看在线播放| 少妇高潮的动态图| 下体分泌物呈黄色| 男女下面进入的视频免费午夜| 男女无遮挡免费网站观看| 国产黄片视频在线免费观看| 在线看a的网站| 久热久热在线精品观看| 久久人人爽人人爽人人片va| 看免费成人av毛片| 成年人午夜在线观看视频| 亚洲天堂av无毛| 国产高清三级在线| 美女xxoo啪啪120秒动态图| 国产av码专区亚洲av| a级毛色黄片| 亚洲怡红院男人天堂| 亚洲人成网站在线观看播放| 69人妻影院| 午夜福利网站1000一区二区三区| 欧美 日韩 精品 国产| 久久人人爽人人爽人人片va| 黄色一级大片看看| 狂野欧美激情性bbbbbb| 菩萨蛮人人尽说江南好唐韦庄| 人人妻人人看人人澡| 能在线免费看毛片的网站| 午夜福利在线在线| 老女人水多毛片| 在线观看美女被高潮喷水网站| 永久网站在线| 亚洲av不卡在线观看| 国产成人精品婷婷| 成人高潮视频无遮挡免费网站| 亚洲国产精品成人久久小说| 男女下面进入的视频免费午夜| 成人午夜精彩视频在线观看| videossex国产| 蜜臀久久99精品久久宅男| 老女人水多毛片| 日韩成人av中文字幕在线观看| 在线看a的网站| 中文字幕免费在线视频6| 日本猛色少妇xxxxx猛交久久| 国产毛片a区久久久久| 国产 精品1| 99久久精品热视频| 熟女av电影| 亚洲最大成人手机在线| 黄片wwwwww| 亚洲av中文字字幕乱码综合| 97热精品久久久久久| 国产精品福利在线免费观看| 黄色欧美视频在线观看| 在线观看免费高清a一片| 神马国产精品三级电影在线观看| 涩涩av久久男人的天堂| 校园人妻丝袜中文字幕| 边亲边吃奶的免费视频| 中文字幕制服av| 精品一区在线观看国产| 国产精品久久久久久久电影| 国产v大片淫在线免费观看| 色哟哟·www| 老司机影院毛片| 干丝袜人妻中文字幕| 午夜激情久久久久久久| 国产黄片美女视频| a级毛色黄片| 国产成人精品婷婷| 人人妻人人澡人人爽人人夜夜| 麻豆精品久久久久久蜜桃| 777米奇影视久久| 成年版毛片免费区| a级毛色黄片| 大片电影免费在线观看免费| 亚洲成色77777| 热re99久久精品国产66热6| 少妇的逼好多水| 99精国产麻豆久久婷婷| 国产一级毛片在线| 成人特级av手机在线观看| 日本免费在线观看一区| 欧美成人精品欧美一级黄| 欧美国产精品一级二级三级 | 女人久久www免费人成看片| 麻豆乱淫一区二区| 欧美精品国产亚洲| 成年版毛片免费区| 日本黄色片子视频| 如何舔出高潮| 成人特级av手机在线观看| 免费观看av网站的网址| 国产91av在线免费观看| 一级二级三级毛片免费看| 日韩欧美 国产精品| 国产日韩欧美亚洲二区| 久久97久久精品|