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

    白洋淀沉積物重金屬潛在生態(tài)風險及生物可利用性分析*

    2022-11-09 03:35:04古小治陳開寧
    湖泊科學 2022年6期
    關鍵詞:淀區(qū)白洋淀點位

    劉 成,黃 蔚,古小治,張 雷, 陳開寧,2**

    (1:中國科學院南京地理與湖泊研究所,湖泊與環(huán)境國家重點實驗室,南京 210008) (2:蘇州科技大學,江蘇水處理技術與材料協(xié)同創(chuàng)新中心,蘇州 215009)

    近二十多年來,我國在河流、湖泊等水環(huán)境治理中取得了諸多重要進展,尤其是在水體富營養(yǎng)化控制方面取得了顯著成效,全國地表水環(huán)境質量不斷提升,水質不達標斷面比例逐漸降低,多個富營養(yǎng)化湖泊水質狀況總體好轉[1-2]. 然而,根據(jù)Huang等[3]對我國142個湖庫的研究表明,雖然富營養(yǎng)化狀況總體逐漸改善,但水體重金屬污染依然較為突出,部分區(qū)域甚至呈現(xiàn)加重趨勢. 湖泊水體中重金屬大多以顆粒物形式沉降至沉積物中[4],由此導致沉積物中重金屬不斷富集,成為水體重金屬重要的“匯”. 而在沉積物-水界面環(huán)境發(fā)生改變時,重金屬又可能以離子擴散、再懸浮等形式向上覆水體中釋放[5-6],導致沉積物成為水體重金屬污染重要的“源”. 因此,在湖泊水體重金屬遷移轉化過程中,沉積物同時扮演著重金屬的匯和源的角色,并不斷發(fā)生著轉變.

    沉積物中重金屬的源-匯轉換現(xiàn)象在淺水湖泊中尤為突出,一方面,淺水湖泊沉積物-水界面環(huán)境極易在人為和自然因素下隨季節(jié)發(fā)生轉變,如藻類和水生植物大量降解期間水體及沉積物氧消耗加劇[7-8],由此可能導致一些Fe/Mn氧化物結合態(tài)金屬溶解釋放[9];另一方面,淺水湖泊長期存在的風浪、底棲等擾動過程導致大量沉積物再懸浮,進入水體的沉積物顆粒由還原環(huán)境向氧化環(huán)境轉變,從而可能導致一些有機物及硫化物結合態(tài)金屬釋放[10]. 此外,淺水湖泊相比深水湖泊更易發(fā)生富營養(yǎng)化,我國淺水湖泊大多處于中營養(yǎng)至富營養(yǎng)狀態(tài)[11],在富營養(yǎng)化發(fā)展過程中,藻類和水生植物大量繁殖、降解導致水體及沉積物中有機質含量上升,近三十年來我國多個湖庫呈現(xiàn)有機質污染不斷加重趨勢[12],不斷升高的有機質進一步加劇了沉積物顆粒與重金屬的耦合和釋放過程[13]. 淺水湖泊中不斷發(fā)生的上述沉降、再懸浮、釋放等過程以及不斷發(fā)展的富營養(yǎng)化進程使得其沉積物重金屬更加容易被生物利用,從而威脅水生態(tài)安全. 因此,在對重金屬污染狀況進行評估的同時,還需要對其生物可利用性進行更加深入地研究.

    使用潛在生態(tài)風險指數(shù)法等基于背景值的評估方法可以較為直觀地反應沉積物重金屬污染狀況[14],而沉積物重金屬對水環(huán)境的影響及對水生生物的威脅則需要進一步研究. 如通過對沉積物重金屬形態(tài)的解析,可以揭示其中可交換態(tài)和易受氧化-還原環(huán)境影響的結合態(tài)賦存特征及釋放潛力[9]. 此外,湖泊水體及沉積物中多種常見有害重金屬(主要包括Cd、Cu、Ni、Pb、Zn等)可與還原態(tài)硫結合,從而降低其生物可利用性,減小對水生態(tài)的威脅. 因而,研究人員進一步提出基于酸可揮發(fā)性硫(acid volatile sulfide, AVS)和同步可提取金屬(simultaneously extracted metals, SEM)的重金屬生物可利用性評價方法,該方法可用于評價與沉積物重金屬生物可利用性直接相關的底棲動物可能受到的潛在毒性影響,在國內外沉積物重金屬生物可利用性研究中獲得了大量應用[15-17].

    白洋淀是華北平原最大的淺水湖泊濕地,素有“華北之腎”之稱,對區(qū)域氣候調節(jié)、水利調控、生物多樣性保護等具有重要的作用. 自1960s以來,由于經(jīng)濟發(fā)展和人類活動的影響,白洋淀水生態(tài)環(huán)境質量逐漸退化[18]. 研究人員針對白洋淀水生態(tài)環(huán)境退化及恢復做了大量工作,包括水環(huán)境質量[19]、生物資源[20]、內源污染[21]、有機污染[22]等方面. 在重金屬污染方面,多位研究人員對白洋淀水體及沉積物重金屬分布特征及潛在生態(tài)風險等進行了研究,如高秋生等[23]和白紅軍等[24]的研究均表明,白洋淀沉積物中As和Cd污染相對較為突出,且北部燒車淀等區(qū)域污染較重;李必才等[25]認為沉積物具有一定的Pb和Zn污染;趙鈺等[26]則發(fā)現(xiàn)沉積物中Cd的可交換及碳酸鹽結合態(tài)和Fe/Mn氧化物結合態(tài)比例較高;薛培英等[27]的研究進一步揭示As在間隙水中濃度較高. 以上研究為白洋淀沉積物重金屬研究及整治提供了豐富的數(shù)據(jù)和堅實的支撐. 然而,目前對白洋淀沉積物重金屬的研究遠少于氮、磷營養(yǎng)鹽等方面的研究[18],已有研究多集中于部分重金屬的污染分析,仍需要在此基礎上更加深入地分析多種重金屬污染賦存總量、形態(tài)及其對水環(huán)境的影響,剖析典型污染重金屬的生物可利用性及其對水生態(tài)的潛在威脅. 因此,本研究在對白洋淀多個區(qū)域沉積物重金屬表層及垂向分布和潛在生態(tài)風險研究的基礎上,結合重金屬形態(tài)和AVS-SEM分析,進一步解析其生物可利用性,以期為白洋淀沉積物重金屬研究及整治提供進一步支撐.

    圖1 白洋淀沉積物柱樣采樣點位分布Fig.1 Distribution of the sediment column sampling sites in Lake Baiyangdian

    1 材料與方法

    1.1 柱狀樣品采集及預處理

    為研究白洋淀沉積物重金屬潛在生態(tài)風險及生物可利用性,于2017年5月在白洋淀主要水域設置了11個采樣點位,采樣點位所在區(qū)域包括燒車淀區(qū)域(B1~B4)、郭里口區(qū)域(B5、B6)、王家寨區(qū)域(B7、B8)和南部淀區(qū)(B9、B10、B11),點位具體分布狀況見圖1. 使用柱狀重力采樣器(Φ90 mm×500 mm)原位采集了每個點位柱狀沉積物樣品,柱樣深度均達到20 cm以上. 柱樣采集后立即密封,并保持原狀轉運至實驗室進行分析.

    沉積物柱樣轉運至實驗室后,按照0~2、2~4、4~6、6~8、8~10、10~15、15~20 cm間隔對柱狀樣品進行切割分層,分層后的樣品立即置于聚乙烯樣品袋混勻. 隨后,立即取部分0~2 cm新鮮樣品置于50 mL離心管,離心獲取表層沉積物間隙水樣品,使用0.45 μm濾膜過濾后,置于10 mL離心管,加一滴濃硝酸(優(yōu)級純)保存?zhèn)錅y. 此外,分取0~2、2~4、4~6、6~8、8~10 cm部分新鮮樣品立即在厭氧條件下分析AVS和SEM. 所有剩余分層樣品置于真空冷凍干燥機進行干燥,干燥后的沉積物樣品使用瑪瑙研缽研磨,并過0.15 mm尼龍分樣篩備測.

    1.2 分析測試方法

    1.2.1 AVS和SEM分析 對表層10 cm沉積物樣品按照2 cm分層后,分析了其中AVS和SEM含量. 首先使用冷擴散法分析了AVS含量[28]:取新鮮表層泥樣(5±0.5) g加入具塞敞口瓶中,將裝有5 mL堿性乙酸鋅的比色管(10 mL)置于瓶中,使用氮氣吹脫形成無氧環(huán)境,加入6 mol/L鹽酸15 mL和0.1 mol/L抗壞血酸溶液2 mL,充分提取16 h后,使用對氨基二甲基苯胺光度法分析堿性乙酸鋅提取液中的S2-[29],進而獲得AVS含量. 將沉積物浸提后的溶液經(jīng)0.45 μm濾膜過濾后,分析其重金屬(Cd、Cu、Ni、Pb、Zn)含量得到SEM含量.

    1.2.2 沉積物重金屬形態(tài)分級 參照BCR金屬形態(tài)分級方法[30]對表層2 cm沉積物樣品重金屬形態(tài)進行了分析,將As、Cd、Cr、Cu、Ni、Pb、Zn 7種常見有害重金屬形態(tài)分級為:F1,可交換態(tài)及碳酸鹽結合態(tài)(0.11 mol/L 醋酸提取);F2,鐵/錳氧化物結合態(tài)(0.1 mol/L 鹽酸羥胺提取);F3,有機物及硫化物結合態(tài)(1 mol/L醋酸銨提取);F4,殘渣態(tài). 其中,F(xiàn)1、F2和F3形態(tài)總體上為可轉化態(tài),分別在酸性、還原性和氧化性增強時有可能溶解釋放至間隙水及上覆水中[9],從而增強金屬的生物可利用性,F(xiàn)4形態(tài)則較為穩(wěn)定.

    1.2.3 沉積物和水樣重金屬及總有機碳(TOC)分析 經(jīng)凍干過篩后的沉積物樣品置于PTFE消解器皿中,使用鹽酸、氫氟酸、硝酸、高氯酸消解后,用電感耦合等離子體光譜儀(ICP-AES,PerkinElmer DV4300,USA)分析Fe、Mn、Ti、V和Zn,檢測限分別為5、0.5、1、2和2 mg/kg;使用電感耦合等離子體質譜儀(ICP-MS; Agilent 7700x, USA)分析Cr、Co、Ni、Cu、As、Mo、Cd、Sb、Tl和Pb,檢測限分別為0.1、0.01、0.05、0.02、0.1、0.05、0.01、0.05、0.02和0.01 mg/kg. 沉積物中Hg的測定使用Hydra-c型全自動測汞儀(Teledyne Leeman Labs,USA)完成(檢測限0.005 mg/kg). 共計分析了沉積物中16種重金屬(其中As為類金屬),并著重分析了As、Cd、Cr、Cu、Hg、Ni、Pb、Zn這8種常見有害重金屬在間隙水中的濃度分布(ICP-MS,水樣中各金屬檢測限分別為0.05、0.005、0.05、0.01、0.02、0.03、0.01、0.1 μg/L). 研究區(qū)域沉積物及間隙水中各重金屬含量均高于所用方法檢測限. 沉積物重金屬分析過程中,每個消解批次隨即抽取10%的樣品進行平行樣質量控制,平行樣相對偏差確保在5%以內,大于5%時對整個批次樣品重新進行消解及分析. 每個消解批次使用水系沉積物成分分析標準物質(GSD-9)保證精確度,并同步處理2份空白樣品. 以上分析方法參考美國環(huán)境保護署(USEPA)的標準方法(Method 200.7、6020B)[31-32],并得到中國合格評定國家認可委員會認可(CNAS L 1628),樣品分析由中國科學院南京地理與湖泊研究所湖泊與環(huán)境國家重點實驗室公共技術服務中心完成. 此外,使用重鉻酸鉀容量法對沉積物TOC進行了分析[29].

    1.3 沉積物重金屬風險及生物可利用性分析

    1.3.1 潛在生態(tài)風險分析 潛在生態(tài)風險指數(shù)(RI)分析參考H?kanson[14]于1980年提出的方法,基于該方法,結合我國《土壤環(huán)境質量:農(nóng)用地土壤污染風險管控標準(試行)》(GB 15618—2018)和《農(nóng)用污泥污染物控制標準》(GB 4284-2018)等土壤及沉積物重金屬污染整治過程中常用參考標準,著重分析了我國沉積物及土壤中8種常見有害重金屬(As、Cd、Cr、Cu、Hg、Ni、Pb和Zn)的潛在生態(tài)風險,分析過程中使用河北省土壤重金屬背景值[33]作為參考. 沉積物重金屬RI值范圍所對應的潛在生態(tài)風險為:RI<150,低潛在生態(tài)風險; 150≤RI<300,中潛在生態(tài)風險; 300≤RI<600,高潛在生態(tài)風險; 600≤RI<1200,很高潛在生態(tài)風險;RI≥1200,極高潛在生態(tài)風險.

    1.3.2 生物可利用性分析 一方面,針對前述8種有害重金屬在沉積物中的賦存形態(tài)及在間隙水中的分布進行了解析,研究重金屬潛在釋放形態(tài)特征;另一方面,針對這些有害重金屬中易與還原態(tài)硫結合的Cd、Cu、Ni、Pb、Zn等[15],使用ΣSEM-AVS、ΣSEM/AVS及ΣSEMx,oc(ΣSEM-AVS/fOC)值評估其生物可利用性[15,34],其中,ΣSEM-AVS<0 μmol/g(dw)或ΣSEM/AVS<1表示沉積物重金屬無毒性;ΣSEM-AVS>0 μmol/g(dw)或ΣSEM/AVS>1表示沉積物重金屬可能對底棲生物或其他生物具有一定的潛在毒性;ΣSEM-AVS>5 μmol/g(dw)或ΣSEM/AVS>2表示沉積物重金屬具有顯著毒性;此外,鑒于有機質與重金屬的耦合特征,當ΣSEMx,oc(ΣSEM-AVS/fOC)超過147.5 μmol/g(OC)時表示沉積物中重金屬具有較大生物可利用性.

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

    點位圖由Surfer軟件(Golden software, Golden, CO, USA)進行繪制;數(shù)據(jù)圖的繪制和Spearman相關性圖制作使用Origin 2021(Origin Lab, Northampton, MA, USA)軟件完成;使用Kruskal-Wallis單因素方差分析方法對沉積物重金屬剖面差異性進行分析,使用聚類分析法(平方歐氏距離)對重金屬關系及來源進行統(tǒng)計分析,以上統(tǒng)計分析使用SPSS 26(IBM, New York, NY, USA)軟件完成.

    2 結果與分析

    2.1 沉積物重金屬含量空間分布及垂向變化特征

    對各點位表層10 cm沉積物中重金屬含量均值進行分析發(fā)現(xiàn)(附表Ⅰ),各金屬在沉積物中含量均值高低次序依次為Fe (29630.50 mg/kg)>Ti (3213.07 mg/kg)>Mn (539.44 mg/kg)>Zn (104.01 mg/kg)>V (76.63 mg/kg)>Cr (52.60 mg/kg)>Cu (43.49 mg/kg)>Ni (35.83 mg/kg)>Pb (26.75 mg/kg)>Co (10.32 mg/kg)>As (8.96 mg/kg)>Mo (2.06 mg/kg)>Sb (1.57 mg/kg) >Tl (0.43 mg/kg)>Cd (0.31 mg/kg)>Hg (0.16 mg/kg). 其中,Cr、Ni、Cu、Zn、Cd、Hg、Pb、Fe、V、Co等在沉積物中的含量總體呈現(xiàn)出明顯的自北向南逐漸降低趨勢,Mo、Sb呈現(xiàn)出自北向南升高的趨勢,而As、Mn、Ti、Tl在各點位的分布則無明顯變化趨勢. 通過與區(qū)域重金屬背景值[33]對比發(fā)現(xiàn),研究淀區(qū)Cr、As、Mn、Ti、Co、Tl含量均值低于各自背景值,而這幾種金屬變異系數(shù)(CV)均低于0.20. 以上分布趨勢表明,Cr、As、Mn、Ti、Co、Tl等低于或接近區(qū)域背景值的金屬在各淀區(qū)沉積物中分布總體較為均勻,其余均值超過相應背景值的各金屬CV多大于0.20,呈現(xiàn)較大空間異質性. 各金屬空間CV與其和背景值的偏差比例間呈現(xiàn)顯著相關性(P<0.01,R2=0.760),可見,重金屬在沉積物中較大的空間異質性主要由不同淀區(qū)污染輸入的差異所致. 其中,北部燒車淀區(qū)域(B1~B4)沉積物重金屬含量明顯高于其他淀區(qū),尤其是位于白溝引河入湖區(qū)域的B1~B3點位,可能由白溝引河入淀污染所致,河流輸入的金屬污染物隨著顆粒物不斷吸附、沉降至河口區(qū)域[35],導致該區(qū)域各金屬含量顯著高于其他區(qū)域.

    對20 cm內各點位、各重金屬的剖面分布特征進一步進行分析,結果發(fā)現(xiàn),各重金屬在不同點位間的垂向剖面分布存在顯著差異(P<0.01),同樣存在較大異質性. 其中,較為典型的8種有害重金屬(As、Cd、Cr、Cu、Hg、Ni、Pb、Zn)剖面分布狀況如附圖I所示. 由圖可見,各區(qū)域重金屬含量自表層2 cm往下至20 cm 變化趨勢相差較大,部分呈現(xiàn)出常見的自表層往下逐漸降低趨勢,同時,亦有部分點位沉積物重金屬含量自表層往下總體相近甚至上升. 表層10 cm內重金屬含量總體相近,大多變化發(fā)生在10 cm以下的沉積物中,部分點位表層0~10 cm重金屬含量呈現(xiàn)出逐漸上升的趨勢,如B2、B3點位Cr、Cu、Zn、As、Cd、Hg、Pb等金屬含量.

    2.2 沉積物重金屬潛在生態(tài)風險分析

    使用潛在生態(tài)風險指數(shù)法對8種常見有害重金屬(As、Cd、Cr、Cu、Hg、Ni、Pb和Zn)在沉積物中的生態(tài)風險進行分析,結果表明(附圖Ⅱ),表層 10 cm沉積物重金屬RI均值高低次序依次為:B3 (1431.1)>B7 (262.2)>B6 (261.2)>B2 (239.9)>B11 (206.5)>B1 (153.3)>B8 (152.3)>B9 (148.5)>B5 (146.6)>B4 (103.3)>B10 (100.2). 從各金屬對RI的平均貢獻看,單個金屬風險值占RI的比例高低次序為:Hg (43.5%)>Cd (41.2%)>Cu (4.3%)>As (3.7%)>Ni (3.1%)>Pb (2.7%)>Cr (0.9%)>Zn (0.6%). 總體上,北部的燒車淀、郭里口、王家寨等區(qū)域RI值高于南部淀區(qū),與各金屬含量趨勢相近,但由于各金屬毒性響應系數(shù)的差異,導致RI趨勢有一定的波動. 其中,除位于燒車淀核心區(qū)域的B3點位外,其余各點位表層10 cm沉積物重金屬RI均值均小于300,生態(tài)風險較低,僅在B6點位的4 cm和6 cm處和B7點位的10 cm 處具有中等風險,RI分別為440.7、314.8和384.3. B3點位沉積物重金屬RI均值達到極高水平,Hg和Cd兩種金屬的RI分別達到1110.9和247.8,占比分別達到70.7%和21.5%,是潛在生態(tài)風險較高的2種金屬. 結合前述重金屬含量分析結果,該點位各金屬含量均較高,尤其是對潛在生態(tài)風險貢獻較大的Hg和Cd 2種金屬,含量顯著高于其他點位,且2種金屬的毒性響應系數(shù)較高[14],由此導致該點位潛在生態(tài)風險尤為突出.

    2.3 沉積物重金屬形態(tài)特征

    表層沉積物通常對氧化還原狀況變化較為敏感,當表層沉積物被氧化時,可能導致金屬硫化物的氧化;而當表層沉積物由于沉降顆粒有機質快速降解等導致耗氧及還原性加強時,則可能導致鐵/錳氧化物結合態(tài)金屬被還原[9]. 這些過程均可能導致表層沉積物中重金屬溶解至間隙水,進而釋放至上覆水[6,36-37],而重金屬的賦存形態(tài)變化則在此過程中具有重要影響. 因此,為了進一步研究重金屬在沉積物中的賦存形態(tài),對表層2 cm典型有害重金屬進行了形態(tài)分析. 其中,Cr主要以殘渣態(tài)為主(圖2),可轉化態(tài)(F1、F2和F3)比例總和最高在30%左右;Ni、Cu、Zn、As、Cd和Pb的可轉化態(tài)比例較高,最高達到90%,而這3種形態(tài)又分別在酸性、還原性和氧化性增強時有可能溶解釋放至間隙水及上覆水中,具有較高生物可利用潛力[9]. 此外,從空間分布上看,Cr、Ni、Cu、Zn可轉化態(tài)比例呈現(xiàn)出明顯的自北往南逐漸降低的趨勢,Cd、Pb的可交換態(tài)及碳酸鹽結合態(tài)(F1)比例也呈現(xiàn)出這一趨勢. 由此可見,北部燒車淀區(qū)域污染輸入的重金屬不僅含量較高,且可轉化態(tài)比例較高. 這些較高比例的可轉化態(tài)重金屬可能更加容易受到沉積物-水系統(tǒng)環(huán)境變化的影響,從而溶解釋放,由此導致表層沉積物中最新沉積的重金屬具有更高的生物可利用潛力.

    圖2 沉積物典型重金屬形態(tài)特征Fig.2 Fraction characteristics of typical heavy metals in the sediment

    圖3 白洋淀沉積物間隙水重金屬濃度Fig.3 Concentrations of heavy metals in the pore water of Lake Baiyangdian

    2.4 間隙水重金屬濃度分布特征

    鑒于表層沉積物在沉積物-水界面重金屬遷移轉化過程中的關鍵作用,對表層沉積物間隙水中As、Cd、Cr、Cu、Hg、Ni、Pb、Zn 8種常見有害重金屬濃度進行分析,結果表明(圖3),間隙水各金屬濃度高低次序依次為:As((17.07±0.23) μg/L)>Hg((2.39±0.94) μg/L)>Ni((2.04±1.48) μg/L)>Cu((1.97±0.92) μg/L)>Zn((1.45±0.75) μg/L)>Cr((0.42±0.23) μg/L)>Pb((0.20±0.12) μg/L)>Cd((0.02±0.01) μg/L). 間隙水中Hg濃度約為地表水V類標準限值(1 μg/L)(GB 3838-2002)的3倍左右,其余金屬均值均在地表水I類限值濃度以下. 其中,濃度較高的As和Hg在沉積物中含量相對低于其他多個金屬,Hg在沉積物中含量遠低于其他金屬,As含量也顯著低于區(qū)域背景值(附表Ⅰ). 可見,白洋淀沉積物中As和Hg與其他金屬相比更加容易溶解至間隙水中,從而可能增加其生物可利用性. 然而,除了Hg外,其他金屬在間隙水中的濃度相對仍處于較低水平,間隙水中溶解態(tài)金屬的風險總體較低. 如Cu的濃度遠低于10 μg/L,根據(jù)Strom等[13]的研究結果,對底棲生物基本無暴露風險.

    2.5 沉積物AVS及SEM分布特征

    圖4 白洋淀沉積物酸可揮發(fā)性硫(AVS)分布特征Fig.4 Distribution characteristics of acid volatile sulfide (AVS) in the sediments of Lake Baiyangdian

    對SEM的分析結果表明(圖5),與沉積物重金屬含量分布相似,燒車淀區(qū)域B1~B3點位沉積物中ΣSEM含量顯著高于其他點位,表層10 cm均值分別為3.30、3.19和6.25 μmol/g(dw),B3點位8 cm處最高達到了8.81 μmol/g(dw). 其余點位ΣSEM含量均在2 μmol/g(dw)以下,各淀區(qū)均值為(2.23±1.53) μmol/g(dw),略高于太湖沉積物ΣSEM均值(1.65 μmol/g(dw))[39]. 其中,各金屬對ΣSEM的貢獻率大小依次為:SEMZn(56.0%)>SEMNi(20.1%)> SEMCu(19.3%)> SEMPb(4.4%)> SEMCd(0.1%). 與太湖、巢湖等長江中下游湖泊相似的是[39],對沉積物重金屬潛在生態(tài)風險貢獻較大的Cd對ΣSEM的貢獻卻較小,在一定程度上降低了Cd的生物有效性. 此外,沉積物中ΣSEM含量大多低于AVS含量,有利于還原態(tài)硫對金屬的固定[15].

    圖5 白洋淀沉積物同步可提取金屬(SEM)分布特征Fig.5 Distribution characteristics of simultaneously extracted metals (SEM) in the sediments of Lake Baiyangdian

    3 討論

    3.1 淺水湖泊高有機質含量沉積物重金屬賦存形態(tài)對生物可利用性的影響

    根據(jù)對沉積物重金屬含量及潛在生態(tài)風險的分析,雖然多個重金屬呈現(xiàn)明顯的污染特征,但除了B3點位外,其余點位潛在生態(tài)風險總體較低. 重金屬潛在生態(tài)風險雖然在一定程度上反映了沉積物中重金屬的污染及生態(tài)危害,但其基于區(qū)域背景值的計算方法往往忽視了重金屬的生物有效性及毒性. 因此,在重金屬含量及潛在生態(tài)風險研究的基礎上,進一步研究其生物有效性,更能全面地反映重金屬污染及其對水生態(tài)環(huán)境健康的威脅. 前述對重金屬形態(tài)特征的分析結果表明,多個金屬可轉化態(tài)含量較高,具有較大的生物可利用潛力. 淺水湖泊通常水柱交換性較好,總體上水柱至水底邊界層區(qū)域氧化條件較好,因此,可氧化態(tài)(有機物及硫化物結合態(tài))金屬的生物可利用性受到了更多關注[45]. 在此過程中,與上覆水柱直接接觸的表層沉積物尤為重要,其中重金屬也更容易被生物利用[10]. 對白洋淀表層10 cm沉積物AVS和SEM的分析表明,沉積物中ΣSEM值大多低于AVS(2.5節(jié)),因此,ΣSEM-AVS值大多低于0,由此導致ΣSEMx,oc大多為負值,遠低于147.5 μmol/g(OC). 僅在B3、B6、B10、B11點位部分層位沉積物中發(fā)現(xiàn)ΣSEM-AVS>0,但均小于5 μmol/g(dw)(圖6),說明大部分區(qū)域沉積物重金屬尚無明顯毒性. 在ΣSEM-AVS>0的點位中,位于郭里口區(qū)域的B6點位ΣSEM/AVS值較高,范圍在0.79~9.05之間,表層10 cm均值為5.05;其ΣSEMx,oc均值為17.48 μmol/g(OC),遠低于147.5 μmol/g(OC). 以上分析表明,基于AVS-SEM方法可以發(fā)現(xiàn),各研究淀區(qū)中僅郭里口區(qū)域B6點位沉積物重金屬具有一定的生物毒性,且并不突出. 造成這一現(xiàn)象的主要原因是白洋淀沉積物中AVS含量較高. 根據(jù)對白洋淀沉積物TOC含量的分析結果(附表I),表層10 cm沉積物中TOC含量顯著高于巢湖、太湖等典型富營養(yǎng)化湖泊中的含量均值[45-46]. 從沉積物TOC平面分布看,并未發(fā)現(xiàn)大部分重金屬所呈現(xiàn)的北部高于南部的現(xiàn)象. 由于白洋淀水生植物生長茂盛,植物衰亡后大量碎屑沉降至沉積物中[47],導致沉積物TOC含量總體較高,可能主要來源于內生性植物碎屑的沉降,在一些河口區(qū)域可能會部分來源于陸源有機質輸入[48]. 因此,分析的16種重金屬中,僅As、Cd、Hg、Pb、Mo和Sb與TOC含量呈顯著正相關(P<0.01,附圖Ⅲ),大多金屬并未呈現(xiàn)出其他湖泊中常見的與TOC顯著相關性[45]. 沉積物中如此高的內生性TOC含量使得表層沉積物長期存在著大量有機質的降解過程,從而加劇沉積物-水界面處氧的消耗[49]. 這一過程進而導致表層沉積物中長期存在著高含量的還原態(tài)硫,并與金屬結合,降低其生物有效性[15]. 此外,白洋淀內部河淀縱橫交錯,有143個淀泊且大多水域面積較小[18],難以出現(xiàn)太湖、巢湖等大型淺水湖泊中常見的頻繁水動力擾動及沉積物再懸浮現(xiàn)象[50]. 重金屬與厭氧沉積物中還原態(tài)硫結合后,經(jīng)再懸浮氧化釋放的現(xiàn)象應沒有其他大型淺水湖泊突出. 因此,對于白洋淀高有機質含量的沉積物,SEM及有機物和硫化物結合態(tài)金屬的生物可利用性總體不突出.

    圖6 沉積物ΣSEM-AVS和ΣSEM/AVS值分布特征Fig.6 Distribution characteristics of ΣSEM-AVS and ΣSEM/AVS in the sediment

    前述2.3節(jié)的分析結果表明,沉積物中Cd、Pb等金屬的可交換態(tài)及碳酸鹽結合態(tài)(F1)和鐵/錳氧化物結合態(tài)(F2)比例較高. 其中,F(xiàn)2在還原態(tài)沉積物中很可能會發(fā)生溶解釋放,但由于沉積物中AVS含量較高,Cd、Cu、Ni、Pb、Zn等F2形態(tài)較高的金屬即使發(fā)生了溶解也會被還原態(tài)硫快速固定,因而這幾種金屬在間隙水中的濃度并不突出. 將沉積物中易與還原態(tài)硫結合的幾種金屬各形態(tài)及總量與AVS、SEM含量及金屬在間隙水中的濃度進行了相關性分析(附圖IV),并未發(fā)現(xiàn)AVS含量與各金屬有機物及硫化物結合態(tài)(F3)之間顯著的相關性,其原因可能是由于沉積物中AVS含量較高,遠超出SEM含量. 然而,各金屬SEM含量與F3形態(tài)之間均呈現(xiàn)出顯著的相關性(P<0.01),說明還原環(huán)境下硫化物對金屬的固定作用,從而降低其生物可利用性. 此外,As不僅F2形態(tài)比例較高,且在還原態(tài)條件下容易溶解并釋放[51],這一現(xiàn)象在其他厭氧沉積物中也時常發(fā)生[4]. 由此導致沉積物間隙水中As濃度較高,與薛培英等[27]的研究結果一致. 在高有機質導致的還原態(tài)界面環(huán)境下,缺少了鐵/錳氧化物的固定作用,Hg也更加容易溶解釋放至間隙水和上覆水中[52]. 因此,雖然As和Hg在沉積物中的總量并不突出,但在間隙水中的濃度較高,從而可能具有更高的生物可利用性. 白洋淀沉積物有機質多來源于水生植物腐爛降解[47],在植物快速降解期間,水體會發(fā)生顯著的pH和溶解氧降低過程[53],進而加劇鐵/錳氧化物結合態(tài)以及可交換態(tài)和碳酸鹽結合態(tài)金屬的溶解,在這些環(huán)境過程變化下,As和Hg等易溶解釋放的金屬更應得到重視.

    3.2 沉積物重金屬來源及其對生物可利用性的影響

    對沉積物中各金屬進行了聚類分析和相關性分析,結果顯示(附圖Ⅲ),除Fe外,其余金屬均呈現(xiàn)出明顯的同源性,各重金屬污染來源總體相近. 相關性分析結果也顯示,大部分金屬之間呈現(xiàn)顯著正相關關系(P<0.05),其中,Mo、Sb和Ti與其他金屬相關性不顯著,甚至出現(xiàn)負相關. Mo和Sb這2種金屬含量呈現(xiàn)自北向南逐漸上升的趨勢,與其他大部分金屬趨勢相反;Ti在各點位間無明顯趨勢,且含量多低于背景值,污染相對較輕. 根據(jù)前述重金屬含量分析結果,研究淀區(qū)沉積物大部分重金屬含量呈現(xiàn)出自北往南逐漸降低的趨勢,且各金屬在空間分布的差異性主要由污染所致而非背景值差異(2.1節(jié)). 已有研究也表明[19,21,23],北部燒車淀、郭里口和王家寨等淀區(qū)營養(yǎng)鹽、重金屬等各類污染狀況總體要重于南部淀區(qū). 因而,各金屬主要污染來源應集中于北部各淀區(qū)附近. 根據(jù)白洋淀周邊水系分布及入淀流量特征[18],大部分入淀河流分布在北部,且自1960s以來,由于多個水利工程的實施以及區(qū)域內水資源短缺影響,導致多條河流斷流或干涸,實際入淀河流僅有6條,多集中于北部,如白溝引河和府河. 北部白溝引河、府河等河流是主要的入淀污染來源,由此導致北部淀區(qū)承接著大部分入淀污染[18]. 因此,大部分入淀污染隨北部河道輸入,導致北部燒車淀等淀區(qū)重金屬含量總體顯著高于南部淀區(qū). 以上分析即表明,北部入淀污染是白洋淀重金屬污染的主要來源. 北部燒車淀區(qū)域附近是重金屬含量最高的區(qū)域,尤其是位于區(qū)域核心的B3點位,重金屬含量及潛在生態(tài)風險均較高,這一方面是由于入淀河流導致的顆粒污染物輸入以及水動力因素的改變,導致該區(qū)域重金屬污染沉降較大,與已有多個河口區(qū)域研究結果相近[54-55];另一方面,根據(jù)已有研究結果[23-24],該區(qū)域內還存在一定的生活污水、漁業(yè)和養(yǎng)殖業(yè)歷史污染,均可能加劇沉積物重金屬的污染.

    北部入淀污染不僅導致燒車淀等區(qū)域沉積物中重金屬含量較高,且其中可轉化態(tài)重金屬的比例亦顯著高于南部淀區(qū). 前述2.3節(jié)的結果表明,Cr、Ni、Cu、Zn等多個金屬的可轉化態(tài)呈現(xiàn)自北往南逐漸降低的趨勢. 其中,潛在生態(tài)風險貢獻率較高的Cd可交換態(tài)及碳酸鹽結合態(tài)比例尤其高,且亦呈現(xiàn)出北部高于南部的趨勢. 重金屬形態(tài)分析所針對的表層沉積物主要為近期沉積的水體顆粒污染物,已有研究表明[56],這些沉降顆粒物通常粒徑較小,多小于63 μm,極易吸附重金屬,且其中吸附的重金屬比歷史沉積中的重金屬具有更高的生物可利用性. 本研究中對表層沉積物金屬形態(tài)的分析也反映出最新沉積的重金屬可轉化態(tài)比例大多比較高,且可轉化態(tài)比例由北向南逐漸降低的趨勢與重金屬含量變化趨勢一致. 由此表明,北部入淀污染不僅是淀區(qū)重金屬污染的主要來源,更顯著增加了表層沉積物重金屬生物可利用潛力. 毛欣等[57]的研究表明,白洋淀的沉積速率約為1 cm/a,在1990年之前的沉積物中污染累積一直處于接近背景值的水平,大部分污染累積發(fā)生在2000年以后,即本研究所關注的20 cm深度范圍內的沉積物. 根據(jù)2.1節(jié)重金屬剖面分布特征的研究結果,表層10 cm內各重金屬變化總體較小,甚至出現(xiàn)表層往下升高的趨勢. 這一趨勢表明,近10年來白洋淀沉積物重金屬富集污染總體趨于穩(wěn)定并逐漸降低的趨勢,在一定程度上反映了整個流域近年來各類污染控制措施不斷完善,從而削減了重金屬污染的匯入及在沉積物中的富集. 而北部淀區(qū)近10年來周邊污染的輸入不僅增加了該區(qū)域沉積物重金屬污染累積量,也是提升該區(qū)域重金屬生物可利用性的主要原因.

    4 結論

    對白洋淀多個淀區(qū)沉積物重金屬含量、分布、潛在生態(tài)風險及生物可利用性進行了分析,結果發(fā)現(xiàn):

    1)沉積物重金屬含量呈現(xiàn)自北向南逐漸降低的趨勢,北部燒車淀區(qū)域是主要的重金屬污染富集區(qū)域. 研究的16種重金屬中,大多呈現(xiàn)出由污染輸入差異所導致的空間分布異質性,大部分污染累積發(fā)生在2000年以后,近10年來沉積物重金屬污染總體趨于穩(wěn)定并逐漸降低的趨勢.

    2)沉積物重金屬潛在生態(tài)風險總體較低,呈現(xiàn)出自北往南生態(tài)風險總體降低的趨勢. 北部燒車淀區(qū)域部分點位潛在生態(tài)風險較高,RI值最高達到1431.1,Hg和Cd是對潛在生態(tài)貢獻最高的2種金屬,在各淀區(qū)平均貢獻率分別達到43.5%和41.2%.

    3)各重金屬污染來源相近,北部入淀污染不僅是重金屬富集總量較高的主要原因,也提升了該區(qū)域重金屬生物可利用潛力. 沉積物可轉化態(tài)重金屬比例大多在30%~90%之間,且各金屬可轉化態(tài)比例總體呈現(xiàn)出自北往南逐漸降低的趨勢.

    4)高有機質含量導致沉積物中AVS含量遠高于SEM,AVS和ΣSEM均值分別為(10.59±6.37)和(2.23±1.53) μmol/g(dw). Cd、Cu、Ni、Pb、Zn等金屬由于還原態(tài)硫的固定而生物可利用性較低. 然而,As和Hg在這樣的環(huán)境下更容易溶解和釋放,在間隙水中的濃度分別達到(17.07±0.23)和(2.39±0.94) μg/L,是潛在生物可利用性相對較高的金屬,在白洋淀這樣高有機質含量的沉積物中應給予更多關注.

    5 附錄

    附圖Ⅰ~Ⅳ和附表Ⅰ見電子版(DOI: 10.18307/2022.0614).

    猜你喜歡
    淀區(qū)白洋淀點位
    1年期國債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    秋韻白洋淀
    河北水利(2024年1期)2024-04-12 07:29:57
    白洋淀生態(tài)環(huán)境危機與治理的歷史考察
    機器人快速示教方法及示教點位姿變換的研究
    美麗河北之白洋淀
    雄安新區(qū)白洋淀冰封期水體污染特征及水質評價*
    湖泊科學(2020年4期)2020-07-19 13:02:34
    雄安新區(qū)首個永久性水利工程通水
    河北畫報(2019年2期)2019-04-03 07:06:08
    美麗的白洋淀
    從白洋淀逃荒到北京
    垂直面內建立基線的特殊點位高程測量法
    国产免费男女视频| 国产单亲对白刺激| 国产精品影院久久| 岛国在线观看网站| 国产在线精品亚洲第一网站| 性少妇av在线| 亚洲一码二码三码区别大吗| 亚洲久久久国产精品| 美女大奶头视频| 亚洲熟女毛片儿| 中文字幕精品免费在线观看视频| 国产激情欧美一区二区| e午夜精品久久久久久久| 18美女黄网站色大片免费观看| 天堂√8在线中文| 国产又色又爽无遮挡免费看| 欧美+亚洲+日韩+国产| 亚洲精品国产区一区二| 又黄又粗又硬又大视频| 亚洲在线自拍视频| 欧美日韩福利视频一区二区| 欧美日韩亚洲综合一区二区三区_| 久久精品aⅴ一区二区三区四区| 新久久久久国产一级毛片| av免费在线观看网站| 精品久久久久久久久久免费视频 | 国产亚洲欧美精品永久| 别揉我奶头~嗯~啊~动态视频| 香蕉丝袜av| 免费高清视频大片| 亚洲欧美日韩无卡精品| 老司机靠b影院| 曰老女人黄片| 日韩免费高清中文字幕av| 在线天堂中文资源库| 91成人精品电影| 妹子高潮喷水视频| 又黄又爽又免费观看的视频| 9191精品国产免费久久| 黄色怎么调成土黄色| 国产日韩一区二区三区精品不卡| 亚洲人成伊人成综合网2020| 别揉我奶头~嗯~啊~动态视频| 亚洲人成电影观看| 夜夜爽天天搞| 亚洲一码二码三码区别大吗| 成人免费观看视频高清| 男女床上黄色一级片免费看| 99精品欧美一区二区三区四区| 80岁老熟妇乱子伦牲交| 国产av在哪里看| www.自偷自拍.com| 亚洲成人免费电影在线观看| 亚洲国产精品合色在线| 久久久久国产精品人妻aⅴ院| 9191精品国产免费久久| 看片在线看免费视频| 又黄又粗又硬又大视频| 亚洲av片天天在线观看| 夜夜躁狠狠躁天天躁| 国产激情久久老熟女| 欧美日韩亚洲国产一区二区在线观看| 精品国产一区二区久久| 99在线视频只有这里精品首页| 成人黄色视频免费在线看| 国产av又大| 久久性视频一级片| 久久精品成人免费网站| 国产精品久久视频播放| 男女下面插进去视频免费观看| 国产午夜精品久久久久久| 久久香蕉国产精品| 十八禁网站免费在线| 日韩欧美三级三区| 午夜福利一区二区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 波多野结衣av一区二区av| 村上凉子中文字幕在线| 国产亚洲精品第一综合不卡| 亚洲少妇的诱惑av| 精品电影一区二区在线| 欧美日韩黄片免| 午夜精品国产一区二区电影| 97碰自拍视频| 欧美+亚洲+日韩+国产| 91成人精品电影| 国产在线观看jvid| 一a级毛片在线观看| 国产主播在线观看一区二区| 婷婷六月久久综合丁香| 亚洲精品一二三| 国产真人三级小视频在线观看| 最好的美女福利视频网| 中文亚洲av片在线观看爽| 欧美精品亚洲一区二区| 老司机亚洲免费影院| 欧美人与性动交α欧美软件| 老司机深夜福利视频在线观看| 久久精品亚洲精品国产色婷小说| 一区二区三区国产精品乱码| 亚洲精品国产色婷婷电影| 精品午夜福利视频在线观看一区| 国产精品亚洲av一区麻豆| 精品久久久久久电影网| 精品久久久久久电影网| 欧美色视频一区免费| 这个男人来自地球电影免费观看| 国产成人精品无人区| 69av精品久久久久久| 最近最新中文字幕大全免费视频| 国产精品九九99| 亚洲av电影在线进入| 国产一区二区三区在线臀色熟女 | 免费在线观看黄色视频的| 欧美色视频一区免费| 久久九九热精品免费| 中文欧美无线码| 久久久国产成人免费| 免费在线观看亚洲国产| 国产麻豆69| 精品无人区乱码1区二区| 男人的好看免费观看在线视频 | 久久中文字幕一级| 亚洲熟妇中文字幕五十中出 | 国产精品亚洲av一区麻豆| 国产区一区二久久| 首页视频小说图片口味搜索| 欧美中文日本在线观看视频| 女人被狂操c到高潮| 免费av中文字幕在线| www.精华液| 国产aⅴ精品一区二区三区波| 日日夜夜操网爽| 国产熟女xx| 一边摸一边抽搐一进一出视频| 又黄又爽又免费观看的视频| 国产精品1区2区在线观看.| 黄网站色视频无遮挡免费观看| 身体一侧抽搐| 亚洲精品一区av在线观看| 精品一区二区三卡| 丰满的人妻完整版| 国产一区二区激情短视频| 久久中文字幕人妻熟女| 最近最新免费中文字幕在线| 黄片大片在线免费观看| 黄色片一级片一级黄色片| 国产亚洲欧美98| 人妻久久中文字幕网| 亚洲欧洲精品一区二区精品久久久| 两个人看的免费小视频| 99国产精品99久久久久| 岛国在线观看网站| 国产亚洲av高清不卡| 国产99白浆流出| 热re99久久精品国产66热6| 日本免费一区二区三区高清不卡 | 妹子高潮喷水视频| 亚洲人成网站在线播放欧美日韩| 国产亚洲av高清不卡| 一边摸一边做爽爽视频免费| 精品一区二区三区视频在线观看免费 | 色综合欧美亚洲国产小说| 日本vs欧美在线观看视频| 狂野欧美激情性xxxx| 久久久国产一区二区| svipshipincom国产片| 欧洲精品卡2卡3卡4卡5卡区| 91大片在线观看| 真人做人爱边吃奶动态| av电影中文网址| 亚洲一区中文字幕在线| 少妇 在线观看| 国产无遮挡羞羞视频在线观看| 一进一出抽搐gif免费好疼 | 69精品国产乱码久久久| 免费不卡黄色视频| 国产精品亚洲av一区麻豆| 成人亚洲精品av一区二区 | 一进一出抽搐gif免费好疼 | 中国美女看黄片| 男人的好看免费观看在线视频 | 国产一区二区三区在线臀色熟女 | 亚洲在线自拍视频| 韩国av一区二区三区四区| 黄色成人免费大全| 亚洲五月色婷婷综合| 女人被躁到高潮嗷嗷叫费观| 国产男靠女视频免费网站| 一进一出好大好爽视频| 99在线人妻在线中文字幕| 男女之事视频高清在线观看| 成年人黄色毛片网站| 亚洲成国产人片在线观看| 久久久久精品国产欧美久久久| 中文字幕av电影在线播放| 久久精品亚洲熟妇少妇任你| 91成人精品电影| 国产精品久久久av美女十八| 久久狼人影院| 18禁国产床啪视频网站| 国产欧美日韩精品亚洲av| 国产亚洲精品久久久久5区| cao死你这个sao货| 在线观看午夜福利视频| 亚洲人成电影免费在线| 国产激情久久老熟女| 一级毛片女人18水好多| 亚洲一区二区三区不卡视频| 亚洲中文日韩欧美视频| xxxhd国产人妻xxx| 国产aⅴ精品一区二区三区波| 高潮久久久久久久久久久不卡| 久久国产亚洲av麻豆专区| 亚洲va日本ⅴa欧美va伊人久久| 精品国产超薄肉色丝袜足j| 亚洲成人免费av在线播放| 精品人妻在线不人妻| 亚洲,欧美精品.| 51午夜福利影视在线观看| 精品熟女少妇八av免费久了| 午夜久久久在线观看| 免费观看精品视频网站| 亚洲国产精品合色在线| 成人国产一区最新在线观看| 国产主播在线观看一区二区| 久久久国产精品麻豆| 老鸭窝网址在线观看| 99热国产这里只有精品6| 熟女少妇亚洲综合色aaa.| 脱女人内裤的视频| 国产精品98久久久久久宅男小说| 成人av一区二区三区在线看| av视频免费观看在线观看| 国产三级黄色录像| 女性生殖器流出的白浆| 久久精品91无色码中文字幕| 69精品国产乱码久久久| 人妻久久中文字幕网| 国产精品野战在线观看 | 国产精品二区激情视频| 日韩视频一区二区在线观看| 国产成人影院久久av| 69精品国产乱码久久久| 国产激情久久老熟女| 欧美日韩国产mv在线观看视频| 亚洲三区欧美一区| 女人高潮潮喷娇喘18禁视频| 中文字幕人妻丝袜一区二区| 少妇裸体淫交视频免费看高清 | 亚洲国产看品久久| 久久久精品欧美日韩精品| 91国产中文字幕| 91成人精品电影| 欧美老熟妇乱子伦牲交| 性色av乱码一区二区三区2| 在线观看免费日韩欧美大片| www.精华液| 亚洲自偷自拍图片 自拍| 久久久久久免费高清国产稀缺| 日韩精品免费视频一区二区三区| 久久久久久久午夜电影 | 成人永久免费在线观看视频| 亚洲人成网站在线播放欧美日韩| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲中文av在线| 18禁裸乳无遮挡免费网站照片 | 国产免费男女视频| 丝袜在线中文字幕| 久热这里只有精品99| 99国产精品免费福利视频| 日韩精品免费视频一区二区三区| 一级毛片高清免费大全| 人人妻人人爽人人添夜夜欢视频| 亚洲精品美女久久久久99蜜臀| 一级,二级,三级黄色视频| av片东京热男人的天堂| 亚洲色图综合在线观看| 久久精品91蜜桃| 黄色女人牲交| 国产一区二区三区综合在线观看| 日本一区二区免费在线视频| 亚洲第一欧美日韩一区二区三区| 精品日产1卡2卡| 国产精品乱码一区二三区的特点 | 99精品久久久久人妻精品| 91老司机精品| 搡老岳熟女国产| 青草久久国产| 一级黄色大片毛片| 91av网站免费观看| 人人澡人人妻人| 精品一区二区三卡| 性欧美人与动物交配| 亚洲av五月六月丁香网| 黑丝袜美女国产一区| 欧美久久黑人一区二区| 精品国产一区二区三区四区第35| 国产精品亚洲av一区麻豆| 丝袜在线中文字幕| 一级,二级,三级黄色视频| 久久精品亚洲精品国产色婷小说| 国产一区二区三区在线臀色熟女 | 88av欧美| 欧美一级毛片孕妇| 韩国精品一区二区三区| 大香蕉久久成人网| 亚洲人成网站在线播放欧美日韩| 啪啪无遮挡十八禁网站| 欧美国产精品va在线观看不卡| 国产一区二区在线av高清观看| 免费av中文字幕在线| 精品久久久久久,| 精品熟女少妇八av免费久了| 老汉色∧v一级毛片| 一级,二级,三级黄色视频| 国产精品乱码一区二三区的特点 | 夜夜夜夜夜久久久久| 12—13女人毛片做爰片一| 视频在线观看一区二区三区| 9热在线视频观看99| 久久中文字幕一级| 精品国产国语对白av| 亚洲久久久国产精品| 成人影院久久| 国产精品国产av在线观看| 大码成人一级视频| 美女大奶头视频| 国产极品粉嫩免费观看在线| 国产精品综合久久久久久久免费 | 性少妇av在线| 亚洲av成人不卡在线观看播放网| 在线观看免费午夜福利视频| 亚洲国产欧美日韩在线播放| 青草久久国产| 国产一区二区三区视频了| 久久亚洲真实| 美女国产高潮福利片在线看| 一级毛片精品| 18禁国产床啪视频网站| 国产片内射在线| 国产精品香港三级国产av潘金莲| 老司机福利观看| 最好的美女福利视频网| 夜夜夜夜夜久久久久| 乱人伦中国视频| 亚洲人成电影免费在线| 成人av一区二区三区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 色哟哟哟哟哟哟| 一级,二级,三级黄色视频| 亚洲午夜精品一区,二区,三区| 精品免费久久久久久久清纯| 热99国产精品久久久久久7| 亚洲人成77777在线视频| 可以在线观看毛片的网站| 亚洲久久久国产精品| 欧美久久黑人一区二区| 99国产精品免费福利视频| 国产精品一区二区在线不卡| 18禁国产床啪视频网站| 国产99久久九九免费精品| 国产91精品成人一区二区三区| www.精华液| 在线观看一区二区三区| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 两人在一起打扑克的视频| 久久欧美精品欧美久久欧美| 性少妇av在线| 日韩成人在线观看一区二区三区| 国产精品亚洲av一区麻豆| 亚洲成人免费av在线播放| 国产麻豆69| 亚洲专区中文字幕在线| 精品国产一区二区久久| 神马国产精品三级电影在线观看 | 美女午夜性视频免费| 国产成人欧美在线观看| 国产在线观看jvid| 人妻久久中文字幕网| 丁香欧美五月| 亚洲精品国产一区二区精华液| 国产不卡一卡二| 大码成人一级视频| 久久精品aⅴ一区二区三区四区| 一级片'在线观看视频| 在线观看一区二区三区激情| 成人亚洲精品av一区二区 | 人人澡人人妻人| 中文字幕人妻丝袜一区二区| 亚洲色图综合在线观看| 亚洲av日韩精品久久久久久密| 中文字幕av电影在线播放| 亚洲精品粉嫩美女一区| 日日摸夜夜添夜夜添小说| 国产精品九九99| 在线观看免费高清a一片| 国产亚洲欧美在线一区二区| 美女福利国产在线| 大陆偷拍与自拍| 韩国av一区二区三区四区| 亚洲九九香蕉| 免费日韩欧美在线观看| 亚洲精品中文字幕一二三四区| 国产精品乱码一区二三区的特点 | 多毛熟女@视频| 一进一出好大好爽视频| 国产成人av激情在线播放| 亚洲精品美女久久久久99蜜臀| 黄色怎么调成土黄色| 国产欧美日韩精品亚洲av| 久久 成人 亚洲| 久久久久久大精品| 国产精品野战在线观看 | 欧美黄色淫秽网站| 亚洲va日本ⅴa欧美va伊人久久| 高清毛片免费观看视频网站 | 岛国在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 欧美精品啪啪一区二区三区| 男人舔女人下体高潮全视频| 免费搜索国产男女视频| 一级毛片女人18水好多| 欧美日本亚洲视频在线播放| 桃红色精品国产亚洲av| 美女大奶头视频| 婷婷丁香在线五月| 男女做爰动态图高潮gif福利片 | 高潮久久久久久久久久久不卡| 国产成年人精品一区二区 | 日本a在线网址| a级毛片黄视频| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 国产欧美日韩一区二区精品| 国产av在哪里看| 精品久久蜜臀av无| 精品国产美女av久久久久小说| xxx96com| 国产成人精品在线电影| 亚洲va日本ⅴa欧美va伊人久久| 老鸭窝网址在线观看| 亚洲精品成人av观看孕妇| 黄片大片在线免费观看| 18禁观看日本| 国产区一区二久久| 午夜a级毛片| 制服人妻中文乱码| 99在线人妻在线中文字幕| 在线观看免费高清a一片| 91av网站免费观看| 一级a爱视频在线免费观看| 男人的好看免费观看在线视频 | 国产单亲对白刺激| 一进一出好大好爽视频| 中文字幕av电影在线播放| 久久中文字幕一级| 久久国产精品影院| 99精品久久久久人妻精品| 国产野战对白在线观看| 国产精品 欧美亚洲| 一二三四社区在线视频社区8| 免费一级毛片在线播放高清视频 | 精品久久蜜臀av无| 黄片小视频在线播放| 免费一级毛片在线播放高清视频 | 国产亚洲欧美精品永久| 亚洲激情在线av| 中文字幕av电影在线播放| 中文字幕人妻丝袜制服| 精品少妇一区二区三区视频日本电影| 69精品国产乱码久久久| 91成人精品电影| 免费久久久久久久精品成人欧美视频| 人人妻人人爽人人添夜夜欢视频| 欧美日韩国产mv在线观看视频| 国产成人精品无人区| 日本黄色视频三级网站网址| 日韩一卡2卡3卡4卡2021年| 美女午夜性视频免费| 9色porny在线观看| 俄罗斯特黄特色一大片| 水蜜桃什么品种好| 一级作爱视频免费观看| 波多野结衣av一区二区av| 日本三级黄在线观看| 亚洲久久久国产精品| 精品久久久精品久久久| 最近最新免费中文字幕在线| 国产精品综合久久久久久久免费 | 黄片播放在线免费| 国产有黄有色有爽视频| 久久香蕉精品热| 国产麻豆69| 在线观看免费视频网站a站| 淫妇啪啪啪对白视频| 亚洲性夜色夜夜综合| av免费在线观看网站| 久久香蕉激情| 精品久久久久久,| 制服人妻中文乱码| 极品教师在线免费播放| 国产精品98久久久久久宅男小说| 欧美日韩一级在线毛片| av在线播放免费不卡| 久久人人爽av亚洲精品天堂| 久久中文字幕人妻熟女| 亚洲中文字幕日韩| 亚洲成a人片在线一区二区| 精品人妻在线不人妻| 一级毛片女人18水好多| 最近最新中文字幕大全免费视频| 超碰97精品在线观看| 性欧美人与动物交配| 在线免费观看的www视频| 色综合婷婷激情| 中亚洲国语对白在线视频| 天天躁夜夜躁狠狠躁躁| 亚洲avbb在线观看| 在线天堂中文资源库| 热99国产精品久久久久久7| 欧美精品一区二区免费开放| 国产亚洲精品第一综合不卡| 制服人妻中文乱码| 香蕉丝袜av| 丰满饥渴人妻一区二区三| 97超级碰碰碰精品色视频在线观看| 在线永久观看黄色视频| 男男h啪啪无遮挡| 日韩av在线大香蕉| 女人被躁到高潮嗷嗷叫费观| 自拍欧美九色日韩亚洲蝌蚪91| 深夜精品福利| 超碰成人久久| 久9热在线精品视频| 国产又爽黄色视频| 窝窝影院91人妻| 欧美精品亚洲一区二区| 国产麻豆69| 欧美午夜高清在线| 在线观看舔阴道视频| 丝袜人妻中文字幕| 午夜免费成人在线视频| 欧美大码av| 久久中文字幕人妻熟女| 精品一区二区三区视频在线观看免费 | 看黄色毛片网站| 一个人观看的视频www高清免费观看 | 女同久久另类99精品国产91| 两性夫妻黄色片| 午夜影院日韩av| 欧美久久黑人一区二区| 18禁美女被吸乳视频| 999精品在线视频| 欧美日韩亚洲国产一区二区在线观看| 99在线人妻在线中文字幕| 欧美日韩av久久| 国产麻豆69| 亚洲,欧美精品.| 又黄又粗又硬又大视频| 国产av又大| 久久国产精品男人的天堂亚洲| 欧美最黄视频在线播放免费 | 高清欧美精品videossex| 两性夫妻黄色片| 91成人精品电影| 黄色怎么调成土黄色| 69av精品久久久久久| 精品电影一区二区在线| 国产精品 欧美亚洲| 极品教师在线免费播放| 身体一侧抽搐| 欧美另类亚洲清纯唯美| 91成年电影在线观看| 亚洲人成网站在线播放欧美日韩| 性欧美人与动物交配| 日韩三级视频一区二区三区| 777久久人妻少妇嫩草av网站| 女警被强在线播放| 国产91精品成人一区二区三区| 国产免费男女视频| 别揉我奶头~嗯~啊~动态视频| 男人舔女人下体高潮全视频| 亚洲午夜理论影院| 中文欧美无线码| av欧美777| 亚洲午夜理论影院| 如日韩欧美国产精品一区二区三区| 久久狼人影院| 男男h啪啪无遮挡| 国产精华一区二区三区| 久久狼人影院| 欧美日韩国产mv在线观看视频| 制服诱惑二区| 婷婷六月久久综合丁香| 精品卡一卡二卡四卡免费| 国产色视频综合| 啦啦啦免费观看视频1| 中文字幕av电影在线播放| 18禁裸乳无遮挡免费网站照片 | 日韩精品中文字幕看吧| 三级毛片av免费| 黄片小视频在线播放| 欧美激情久久久久久爽电影 | 国产亚洲精品一区二区www| 一区在线观看完整版| 两个人免费观看高清视频| 在线观看舔阴道视频| 69精品国产乱码久久久| 一进一出好大好爽视频| 免费日韩欧美在线观看| 久久久精品国产亚洲av高清涩受| 人人妻人人爽人人添夜夜欢视频| 精品福利观看|