• <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在线| 欧美日韩视频精品一区| 亚洲五月天丁香| 久久人妻av系列| 黄色怎么调成土黄色| 久久午夜亚洲精品久久| 精品国产一区二区久久| 99国产精品免费福利视频| 成人特级黄色片久久久久久久| 免费观看精品视频网站| 黄色a级毛片大全视频| 亚洲欧美日韩高清在线视频| 人人妻人人添人人爽欧美一区卜| 亚洲熟女毛片儿| 99re在线观看精品视频| 亚洲精品中文字幕一二三四区| 很黄的视频免费| 窝窝影院91人妻| 黑人巨大精品欧美一区二区mp4| 国产不卡一卡二| 香蕉丝袜av| 亚洲欧美激情综合另类| 一本大道久久a久久精品| 人人妻,人人澡人人爽秒播| 亚洲av成人av| 99热国产这里只有精品6| 久久热在线av| 丝袜在线中文字幕| 90打野战视频偷拍视频| 久久久久久亚洲精品国产蜜桃av| 亚洲av片天天在线观看| 精品人妻1区二区| 美女高潮到喷水免费观看| 午夜精品在线福利| 一级毛片精品| av中文乱码字幕在线| 国产精品一区二区在线观看99| 精品国产国语对白av| 精品视频人人做人人爽| 美女福利国产在线| 精品欧美一区二区三区在线| 在线观看免费午夜福利视频| 欧美激情高清一区二区三区| 国产91精品成人一区二区三区| 欧美黄色片欧美黄色片| 久久天堂一区二区三区四区| 俄罗斯特黄特色一大片| 久久久久久亚洲精品国产蜜桃av| 女性被躁到高潮视频| 国产成人系列免费观看| 国产精品 国内视频| 一级片'在线观看视频| 91大片在线观看| 丰满迷人的少妇在线观看| av天堂久久9| 最近最新免费中文字幕在线| 欧美成狂野欧美在线观看| 国产亚洲欧美98| 正在播放国产对白刺激| 正在播放国产对白刺激| 怎么达到女性高潮| 精品一区二区三区av网在线观看| 亚洲精品久久午夜乱码| 成人永久免费在线观看视频| 99国产精品一区二区三区| 757午夜福利合集在线观看| 视频区图区小说| 亚洲熟女精品中文字幕| 久久久精品国产亚洲av高清涩受| 国产精品九九99| 国产男女超爽视频在线观看| 久久这里只有精品19| 丝袜美足系列| 久久狼人影院| 18在线观看网站| 中文字幕最新亚洲高清| 老熟妇仑乱视频hdxx| 亚洲精品国产精品久久久不卡| 免费在线观看黄色视频的| 人妻丰满熟妇av一区二区三区 | 亚洲国产毛片av蜜桃av| 久久亚洲真实| 久久影院123| 久久香蕉精品热| 久久精品国产清高在天天线| 男人操女人黄网站| 91精品国产国语对白视频| 国产精品影院久久| 国精品久久久久久国模美| e午夜精品久久久久久久| 久久久久精品国产欧美久久久| 国产在线精品亚洲第一网站| 夫妻午夜视频| 黄色a级毛片大全视频| 激情在线观看视频在线高清 | 日本a在线网址| 亚洲av成人av| 啦啦啦免费观看视频1| 99热只有精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 高清在线国产一区| www.自偷自拍.com| 亚洲aⅴ乱码一区二区在线播放 | 天堂俺去俺来也www色官网| 美女高潮喷水抽搐中文字幕| 99精品久久久久人妻精品| 一级毛片女人18水好多| 久久久国产欧美日韩av| 免费看十八禁软件| 美女福利国产在线| 国产主播在线观看一区二区| 一级作爱视频免费观看| 欧美 亚洲 国产 日韩一| 美女午夜性视频免费| 久久久久久久久久久久大奶| 国产激情欧美一区二区| 精品国产一区二区三区四区第35| 亚洲av第一区精品v没综合| 亚洲精品国产区一区二| 久久热在线av| 女人爽到高潮嗷嗷叫在线视频| 久久ye,这里只有精品| 国产日韩一区二区三区精品不卡| 一级毛片女人18水好多| 久久人人爽av亚洲精品天堂| 亚洲精品国产色婷婷电影| 9色porny在线观看| 高清欧美精品videossex| 国产成人欧美| 久久久水蜜桃国产精品网| 视频区图区小说| av不卡在线播放| 欧美黑人精品巨大| 国产欧美日韩综合在线一区二区| 咕卡用的链子| 丰满的人妻完整版| 国产麻豆69| 午夜91福利影院| 久久国产精品男人的天堂亚洲| 丝袜人妻中文字幕| 午夜日韩欧美国产| 日韩欧美一区二区三区在线观看 | 亚洲一区中文字幕在线| 亚洲久久久国产精品| 黄色片一级片一级黄色片| 99国产精品一区二区三区| 一区二区日韩欧美中文字幕| 别揉我奶头~嗯~啊~动态视频| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 视频区图区小说| 欧美精品啪啪一区二区三区| www.999成人在线观看| 黄色片一级片一级黄色片| 亚洲熟妇中文字幕五十中出 | 免费av中文字幕在线| 男女高潮啪啪啪动态图| 国产高清videossex| 久久久久视频综合| 亚洲精品久久午夜乱码| 女性被躁到高潮视频| 欧美亚洲日本最大视频资源| 国产有黄有色有爽视频| 夜夜躁狠狠躁天天躁| 国产精品影院久久| 国产黄色免费在线视频| 午夜福利影视在线免费观看| 99re6热这里在线精品视频| 99热网站在线观看| 精品一区二区三区视频在线观看免费 | 老司机靠b影院| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 自线自在国产av| 在线看a的网站| 50天的宝宝边吃奶边哭怎么回事| 天天躁狠狠躁夜夜躁狠狠躁| 如日韩欧美国产精品一区二区三区| 午夜日韩欧美国产| 人人妻人人爽人人添夜夜欢视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品香港三级国产av潘金莲| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区mp4| 日韩熟女老妇一区二区性免费视频| 欧美丝袜亚洲另类 | 国产麻豆69| 国产精品久久电影中文字幕 | 中文字幕高清在线视频| 欧美日韩精品网址| av中文乱码字幕在线| 久久久水蜜桃国产精品网| 在线观看午夜福利视频| 国产男女内射视频| 国产人伦9x9x在线观看| 国产午夜精品久久久久久| 国产成人精品久久二区二区91| 麻豆国产av国片精品| 高清在线国产一区| 国产精品久久久久久精品古装| 亚洲人成77777在线视频| 成人三级做爰电影| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 怎么达到女性高潮| 亚洲成人免费av在线播放| 99re6热这里在线精品视频| av天堂在线播放| 国产乱人伦免费视频| 亚洲全国av大片| 日本欧美视频一区| 亚洲美女黄片视频| 免费在线观看影片大全网站| 国产真人三级小视频在线观看| 老司机影院毛片| 亚洲人成电影免费在线| 国产高清国产精品国产三级| 亚洲色图综合在线观看| 精品国产亚洲在线| 久久久精品免费免费高清| 午夜福利在线免费观看网站| 一级作爱视频免费观看| 国产精品1区2区在线观看. | 日韩一卡2卡3卡4卡2021年| 变态另类成人亚洲欧美熟女 | 在线天堂中文资源库| 欧美成人免费av一区二区三区 | 亚洲av欧美aⅴ国产| 国精品久久久久久国模美| 久久亚洲精品不卡| 成人国语在线视频| 大型黄色视频在线免费观看| 久久性视频一级片| 水蜜桃什么品种好| 18禁美女被吸乳视频| xxxhd国产人妻xxx| 欧美黄色片欧美黄色片| 人妻久久中文字幕网| 飞空精品影院首页| 最新的欧美精品一区二区| 成人特级黄色片久久久久久久| 麻豆av在线久日| 色综合欧美亚洲国产小说| 欧美乱色亚洲激情| 国产免费av片在线观看野外av| 后天国语完整版免费观看| 丁香六月欧美| 国产精品99久久99久久久不卡| 国产日韩欧美亚洲二区| 黑丝袜美女国产一区| 夜夜夜夜夜久久久久| 少妇猛男粗大的猛烈进出视频| 免费在线观看影片大全网站| 99精品欧美一区二区三区四区| 可以免费在线观看a视频的电影网站| 欧美日韩黄片免| 在线看a的网站| 亚洲成人手机| 五月开心婷婷网| 男女高潮啪啪啪动态图| av欧美777| 国产成人免费观看mmmm| 老熟女久久久| 视频区图区小说| 妹子高潮喷水视频| 嫁个100分男人电影在线观看| 美国免费a级毛片| 久久久久精品人妻al黑| 国产aⅴ精品一区二区三区波| 国产xxxxx性猛交| 亚洲人成伊人成综合网2020| 成人av一区二区三区在线看| 久久性视频一级片| 日韩中文字幕欧美一区二区| 亚洲精品自拍成人| 9热在线视频观看99| 俄罗斯特黄特色一大片| 91在线观看av| 国产人伦9x9x在线观看| 久久九九热精品免费| 久久久国产成人精品二区 | 午夜精品久久久久久毛片777| 欧美最黄视频在线播放免费 | 老司机亚洲免费影院| 国产欧美日韩精品亚洲av| www.熟女人妻精品国产| 久热这里只有精品99| 久久久久视频综合| 亚洲国产毛片av蜜桃av| 国产97色在线日韩免费| 国产男靠女视频免费网站| 午夜精品在线福利| 欧美日韩黄片免| 亚洲精品自拍成人| 在线观看66精品国产| 久久精品国产亚洲av高清一级| 在线国产一区二区在线| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情 高清一区二区三区| 人妻久久中文字幕网| 涩涩av久久男人的天堂| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 国产av又大| 变态另类成人亚洲欧美熟女 | 人人妻人人澡人人爽人人夜夜| 9色porny在线观看| 深夜精品福利| 久久久久久久久免费视频了| 99精品欧美一区二区三区四区| svipshipincom国产片| 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| 亚洲伊人色综图| 欧美av亚洲av综合av国产av| av在线播放免费不卡| 夜夜爽天天搞| 午夜久久久在线观看| 精品福利观看| 国产精品香港三级国产av潘金莲| 久久人妻熟女aⅴ| 91国产中文字幕| 欧美大码av| 亚洲全国av大片| 免费观看人在逋| 国产成人欧美在线观看 | 成年人黄色毛片网站| 老司机在亚洲福利影院| 国产精品 国内视频| 美女国产高潮福利片在线看| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 中文字幕人妻丝袜制服| 757午夜福利合集在线观看| 中文字幕人妻丝袜制服| 老司机亚洲免费影院| 黄频高清免费视频| 精品久久久精品久久久| 久久人妻熟女aⅴ| 美国免费a级毛片| 黑丝袜美女国产一区| 在线国产一区二区在线| 久久中文字幕人妻熟女| 午夜福利影视在线免费观看| 欧美老熟妇乱子伦牲交| 久久午夜综合久久蜜桃| 欧美激情 高清一区二区三区| 正在播放国产对白刺激| 亚洲av美国av| av国产精品久久久久影院| 久久久久精品人妻al黑| 日日摸夜夜添夜夜添小说| www.999成人在线观看| 黄色a级毛片大全视频| 人人妻人人爽人人添夜夜欢视频| 免费在线观看日本一区| 久久久久久久午夜电影 | 男人舔女人的私密视频| 国产精品秋霞免费鲁丝片| 怎么达到女性高潮| 成人18禁高潮啪啪吃奶动态图| 亚洲午夜精品一区,二区,三区| 日本五十路高清| 丝袜美腿诱惑在线| 中出人妻视频一区二区| 午夜视频精品福利| 最近最新中文字幕大全免费视频| 亚洲精品久久午夜乱码| 91麻豆精品激情在线观看国产 | 夜夜爽天天搞| 免费在线观看日本一区| 精品亚洲成国产av| 午夜精品在线福利| 高清欧美精品videossex| 午夜免费观看网址| 亚洲在线自拍视频| 黄色视频不卡| 9191精品国产免费久久| 精品免费久久久久久久清纯 | 高清毛片免费观看视频网站 | 新久久久久国产一级毛片| 久久这里只有精品19| 日本五十路高清| 十分钟在线观看高清视频www| 天天躁狠狠躁夜夜躁狠狠躁| 免费久久久久久久精品成人欧美视频| 精品国产美女av久久久久小说| 正在播放国产对白刺激| 91精品三级在线观看| 我的亚洲天堂| 日韩熟女老妇一区二区性免费视频| 国产成人免费无遮挡视频| 久久精品人人爽人人爽视色| 一本大道久久a久久精品| а√天堂www在线а√下载 | 一级作爱视频免费观看| 午夜免费成人在线视频| 精品久久蜜臀av无| 丰满的人妻完整版| 免费在线观看日本一区| 日韩视频一区二区在线观看| 国产免费男女视频| 亚洲欧美激情综合另类| 久久精品亚洲熟妇少妇任你| 日本欧美视频一区| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 一区福利在线观看| tocl精华| 久久久久国产精品人妻aⅴ院 | 日韩欧美在线二视频 | 亚洲欧美激情综合另类| 超碰97精品在线观看| 国产精品 国内视频| 亚洲专区字幕在线| 亚洲熟妇中文字幕五十中出 | 热re99久久国产66热| 国产高清激情床上av| 最新在线观看一区二区三区| 一级黄色大片毛片| 老司机影院毛片| 亚洲 欧美一区二区三区| 80岁老熟妇乱子伦牲交| 在线av久久热| 男人舔女人的私密视频| 一区在线观看完整版| 男人舔女人的私密视频| 中文字幕色久视频| 欧美午夜高清在线| 丁香六月欧美| 国产伦人伦偷精品视频| 亚洲欧美日韩另类电影网站| 一级作爱视频免费观看| 午夜福利视频在线观看免费| 午夜久久久在线观看| 久久中文看片网| 午夜免费鲁丝| 叶爱在线成人免费视频播放| 曰老女人黄片| 亚洲五月婷婷丁香| 亚洲全国av大片| 18禁裸乳无遮挡动漫免费视频| 午夜福利乱码中文字幕| 夜夜躁狠狠躁天天躁| 午夜福利免费观看在线| 岛国在线观看网站| 搡老岳熟女国产| 人妻 亚洲 视频| 亚洲免费av在线视频| 亚洲精品成人av观看孕妇| a级片在线免费高清观看视频| 国产在线一区二区三区精| 操出白浆在线播放| 国产成人系列免费观看| 国产高清videossex| 欧美黑人欧美精品刺激| 色综合婷婷激情| 亚洲国产精品一区二区三区在线| 亚洲三区欧美一区| 视频区图区小说| 亚洲国产毛片av蜜桃av| 欧美精品av麻豆av| 国产一区二区激情短视频| 久久精品亚洲精品国产色婷小说| 1024视频免费在线观看| 中文字幕制服av| 国产成人免费观看mmmm| 国产精品久久久人人做人人爽| xxxhd国产人妻xxx| 91av网站免费观看| 国产精品亚洲一级av第二区| 热99re8久久精品国产| 日本欧美视频一区| 久久久久久久久免费视频了| 国产精品自产拍在线观看55亚洲 | 一级片免费观看大全| 999久久久精品免费观看国产| 亚洲人成伊人成综合网2020| 在线观看免费视频网站a站| 国产成人av激情在线播放| 露出奶头的视频| 欧美不卡视频在线免费观看 | 啦啦啦在线免费观看视频4| 香蕉丝袜av| 亚洲欧美日韩另类电影网站| 18禁黄网站禁片午夜丰满| 日本欧美视频一区| 我的亚洲天堂| 国产伦人伦偷精品视频| 丰满饥渴人妻一区二区三| 日本wwww免费看| 亚洲片人在线观看| 日日爽夜夜爽网站| 亚洲一码二码三码区别大吗| 国产精品永久免费网站| 成人精品一区二区免费| 女性生殖器流出的白浆| 啦啦啦视频在线资源免费观看| 午夜福利欧美成人| 亚洲免费av在线视频| 午夜免费成人在线视频| 十分钟在线观看高清视频www| 91字幕亚洲| bbb黄色大片| 精品国产亚洲在线| 制服诱惑二区| 亚洲三区欧美一区| 又大又爽又粗| 黑人巨大精品欧美一区二区mp4| 又大又爽又粗| 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 国产亚洲一区二区精品| 国产精品成人在线| 欧美亚洲日本最大视频资源| 精品一区二区三区四区五区乱码| 日韩欧美一区二区三区在线观看 | 国产高清激情床上av| 啦啦啦视频在线资源免费观看| 久久久久久久精品吃奶| 国产欧美日韩一区二区三区在线| 国产精品免费大片| 亚洲av日韩精品久久久久久密| 97人妻天天添夜夜摸| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 美女福利国产在线| 亚洲一区二区三区欧美精品| 午夜免费成人在线视频| 91大片在线观看| 丝袜美腿诱惑在线| 天堂俺去俺来也www色官网| 首页视频小说图片口味搜索| 中文字幕人妻丝袜一区二区| 老司机亚洲免费影院| 久久香蕉激情| av在线播放免费不卡| 黄色毛片三级朝国网站| cao死你这个sao货| 午夜视频精品福利| 多毛熟女@视频| av有码第一页| 一级a爱视频在线免费观看| 亚洲人成伊人成综合网2020| 宅男免费午夜| 欧美 亚洲 国产 日韩一| 亚洲熟女精品中文字幕| 国产淫语在线视频| 亚洲中文日韩欧美视频| a在线观看视频网站| 久久久国产欧美日韩av| 精品国产乱子伦一区二区三区| 高清黄色对白视频在线免费看| 午夜精品国产一区二区电影| 免费av中文字幕在线| 老司机影院毛片| 建设人人有责人人尽责人人享有的| 在线av久久热| 老汉色av国产亚洲站长工具| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色片一级片一级黄色片| 久久 成人 亚洲| videos熟女内射| 国产精品久久电影中文字幕 | 午夜福利免费观看在线| 免费一级毛片在线播放高清视频 | 啦啦啦视频在线资源免费观看| 国产成人精品在线电影| 亚洲专区字幕在线| 欧美精品人与动牲交sv欧美| 黄网站色视频无遮挡免费观看| 一级作爱视频免费观看| 日本五十路高清| 成人三级做爰电影| 80岁老熟妇乱子伦牲交| 搡老熟女国产l中国老女人| 欧美乱妇无乱码| 日本欧美视频一区| 99国产精品99久久久久| 一边摸一边抽搐一进一小说 | 黑人猛操日本美女一级片| 色94色欧美一区二区| 建设人人有责人人尽责人人享有的| 中出人妻视频一区二区| 久久久久久久午夜电影 | 亚洲第一青青草原| 9191精品国产免费久久| 久久久久久久久久久久大奶| 妹子高潮喷水视频| 亚洲,欧美精品.| 视频区图区小说| 精品久久久久久,| 免费久久久久久久精品成人欧美视频| 国产91精品成人一区二区三区| 淫妇啪啪啪对白视频| 视频区欧美日本亚洲| 国产免费男女视频| 久久天躁狠狠躁夜夜2o2o| 在线观看www视频免费| 国产成人一区二区三区免费视频网站| 色在线成人网| 亚洲国产中文字幕在线视频| 99热只有精品国产| netflix在线观看网站| av网站免费在线观看视频| 老汉色∧v一级毛片| 最新在线观看一区二区三区| 久久精品国产亚洲av香蕉五月 | 久久久精品免费免费高清| 国产在视频线精品| 一进一出抽搐动态|