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

    多因子影響下銅的水質(zhì)基準(zhǔn)及生態(tài)風(fēng)險(xiǎn)

    2022-07-18 14:13:48李亞兵劉紅玲
    中國(guó)環(huán)境科學(xué) 2022年7期
    關(guān)鍵詞:多因子基準(zhǔn)值水生

    劉 娜,李亞兵,劉紅玲,2*

    多因子影響下銅的水質(zhì)基準(zhǔn)及生態(tài)風(fēng)險(xiǎn)

    劉 娜1,李亞兵1,劉紅玲1,2*

    (1.南京大學(xué)環(huán)境學(xué)院,污染控制與資源化研究國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023;2.江蘇省生態(tài)環(huán)境保護(hù)化學(xué)品安全與健康風(fēng)險(xiǎn)研究重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023)

    為探究特定水環(huán)境中多因子影響下的銅水質(zhì)基準(zhǔn)及生態(tài)風(fēng)險(xiǎn),采用物種權(quán)重敏感度分布法、水效應(yīng)比法和生物配體模型推導(dǎo)保護(hù)太湖水生生物銅的水質(zhì)基準(zhǔn).根據(jù)推導(dǎo)結(jié)果,推薦采用最大濃度基準(zhǔn)值(CMC) 1.43 μg/L和持續(xù)濃度基準(zhǔn)值(CCC) 1.33 μg/L.結(jié)合水效應(yīng)比法和生物配體模型,采用聯(lián)合概率法評(píng)估太湖銅的生態(tài)風(fēng)險(xiǎn).結(jié)果表明,兩種方法下豐水期5%水生生物受到銅慢性毒性影響的概率分別為23.43%和39.43%,而未考慮多因子影響的風(fēng)險(xiǎn)概率為85.01%,高估了太湖銅生態(tài)風(fēng)險(xiǎn).可見(jiàn),水環(huán)境多因子對(duì)水質(zhì)基準(zhǔn)和生態(tài)風(fēng)險(xiǎn)的影響不容忽視,我國(guó)目前使用的銅標(biāo)準(zhǔn)可能無(wú)法保護(hù)特定區(qū)域的水生生物.考慮多因子影響可提高基準(zhǔn)值推導(dǎo)和風(fēng)險(xiǎn)評(píng)估的科學(xué)性,避免“過(guò)保護(hù)”和“欠保護(hù)”現(xiàn)象.

    銅;多因子影響;水質(zhì)基準(zhǔn);生態(tài)風(fēng)險(xiǎn);物種權(quán)重敏感度分布法;水效應(yīng)比法;生物配體模型

    銅是生物必需的微量元素,但攝入過(guò)量會(huì)危害生物體,尤其對(duì)水生生物會(huì)產(chǎn)生不同水平的毒性效應(yīng)[1-3].已有研究評(píng)估了多種金屬對(duì)水生生物的毒性,銅的潛在危害最大[4].最初研究關(guān)注銅本身給水生生物帶來(lái)的危害,近年來(lái)大量研究表明生物因素(不同物種、不同生命階段等)和區(qū)域水環(huán)境特征(溫度、硬度、溶解氧、pH值和有機(jī)質(zhì)等)能對(duì)銅的生物有效性及毒性產(chǎn)生影響[5-7],這也推動(dòng)了方法學(xué)的發(fā)展.其中,生物配體模型(BLM)考慮顆粒吸附、絡(luò)合作用、無(wú)機(jī)配體結(jié)合、競(jìng)爭(zhēng)活性位點(diǎn)等多個(gè)影響因素,能較好地預(yù)測(cè)重金屬對(duì)水生生物的毒性[8];物種權(quán)重敏感度分布法(WSSD)和水效應(yīng)比法(WER)則通過(guò)鏈接真實(shí)水環(huán)境,對(duì)毒性值進(jìn)行校正[9-10].三種方法也逐漸應(yīng)用于水質(zhì)基準(zhǔn)的研究中,但都集中于單一方法的應(yīng)用,缺乏對(duì)不同方法的適用性討論.如利用BLM對(duì)金屬水質(zhì)基準(zhǔn)進(jìn)行修訂[11];通過(guò)WSSD法得到校驗(yàn)后太湖鉛的水質(zhì)基準(zhǔn)[12];利用WER法所得銅水質(zhì)基準(zhǔn)值較U.S. EPA國(guó)家標(biāo)準(zhǔn)值更為嚴(yán)格[13].另一方面,我國(guó)現(xiàn)行水質(zhì)標(biāo)準(zhǔn)采用的是統(tǒng)一標(biāo)準(zhǔn)值,缺乏以保護(hù)水生生物為目標(biāo)的銅水質(zhì)標(biāo)準(zhǔn)[14].因此, 探討不同方法的適用性,研究多因子影響下銅的水質(zhì)基準(zhǔn)和生態(tài)風(fēng)險(xiǎn)十分必要.

    太湖作為我國(guó)最大的淡水湖之一,也是受銅等重金屬污染最嚴(yán)重的水域之一[15].本文以太湖為例,在收集銅的急慢性毒性數(shù)據(jù)基礎(chǔ)上,補(bǔ)充了6種太湖本土生物急慢性毒性數(shù)據(jù),考慮生物類群和水化學(xué)等多因子對(duì)水質(zhì)基準(zhǔn)的影響,采用WSSD法、WER以及BLM模型法,得到基于多因子影響的保護(hù)太湖水生生物銅的水質(zhì)基準(zhǔn).進(jìn)一步基于WER法和BLM模型法,采用聯(lián)合概率法對(duì)太湖豐枯水期銅的生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià),并對(duì)不同方法進(jìn)行討論.以期為太湖銅水質(zhì)基準(zhǔn)的制定提供參考,同時(shí)為多因子影響下我國(guó)保護(hù)水生生物重金屬類水質(zhì)標(biāo)準(zhǔn)的制定及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)提供方法借鑒,支撐不同水域污染下“因地制宜”式的控制與管理.

    1 材料與方法

    1.1 水樣采集及參數(shù)測(cè)定

    為推導(dǎo)太湖銅水質(zhì)基準(zhǔn),采集太湖重污染區(qū)陳東港、大浦港、分水大橋、分水橋和平臺(tái)山5 個(gè)代表點(diǎn)水樣,如圖1(a).溫度等參數(shù)及DOC分別采用YSI多參數(shù)水質(zhì)儀(6600V2-4, Ecosense, Ohio, USA)和島津TOC-5000A總有機(jī)碳測(cè)定儀(Shimadzu,京都,日本)進(jìn)行測(cè)定.為測(cè)定太湖Cu2+的環(huán)境濃度以評(píng)價(jià)太湖銅的水生生態(tài)風(fēng)險(xiǎn),考慮水生態(tài)功能區(qū),選擇太湖52個(gè)點(diǎn)位,分別于豐水期(2019年9月)和枯水期(2020年1月)依據(jù)HJ/T91相關(guān)規(guī)定[16]進(jìn)行采樣.采樣點(diǎn)位如圖1(b),其中湖體37個(gè)(綠色點(diǎn)為9個(gè)國(guó)控?cái)嗝?,湖口及周邊河流15個(gè)(紅色點(diǎn)).水樣采集后經(jīng)0.45μm水性微孔濾膜過(guò)濾,濾液中加入適量硝酸使其硝酸含量達(dá)1%作為測(cè)定樣品.銅的濃度利用電感耦合等離子體質(zhì)譜(ICP-MS (NexION 300X))測(cè)定.測(cè)定的標(biāo)準(zhǔn)曲線具有良好的線性(2=0.998),檢出限為0.18μg/L,重復(fù)樣品的標(biāo)準(zhǔn)偏差RSD<5.50%,回收率在90.5%~110.1%,本方法準(zhǔn)確性和穩(wěn)定性較好.

    1.2 太湖銅的水質(zhì)基準(zhǔn)

    1.2.1 毒性數(shù)據(jù)的收集 通過(guò)U.S.EPA的ECOTOX數(shù)據(jù)庫(kù)(http://cfpub.epa.gov/ecotox/)搜集太湖本土水生生物銅(Cu2+)的毒性數(shù)據(jù),數(shù)據(jù)篩選原則如下:

    (1) 剔除非太湖本土生物毒性數(shù)據(jù).

    (2) 生物毒性測(cè)試符合國(guó)標(biāo)[17]或OECD[18]及U.S. EPA的標(biāo)準(zhǔn)(包括受試生物年齡及健康狀態(tài)、空白對(duì)照的要求、測(cè)試時(shí)長(zhǎng)、生物培養(yǎng)、溫度、光照、喂食與否等).

    (3) 測(cè)試暴露介質(zhì)為淡水,測(cè)試地點(diǎn)為實(shí)驗(yàn)室,毒性數(shù)據(jù)包含明確毒性終點(diǎn)、毒性單位、暴露時(shí)間及暴露方式.

    (4) 對(duì)于急性毒性數(shù)據(jù),毒性終點(diǎn)選擇半數(shù)致死濃度(LC50)或半數(shù)有效應(yīng)濃度(EC50).優(yōu)先選擇暴露96h的魚類和兩棲類的LC50,水生無(wú)脊椎動(dòng)物96h(水溞類48h)的LC50,水生植物37d及水生藻類£96h的EC50.對(duì)于慢性毒性數(shù)據(jù),選擇NOEC(無(wú)觀察效應(yīng)濃度)或LOEC(最低效應(yīng)濃度)為毒性終點(diǎn).毒性指標(biāo)選擇短生命周期生物>4d及長(zhǎng)生命周期生物>7d的亞致死效應(yīng),優(yōu)先選擇標(biāo)準(zhǔn)方法獲得的大型溞21d毒性數(shù)據(jù)[18].

    1.2.2 急性毒性試驗(yàn) 對(duì)太湖6種本土水生生物(大型溞、青鳉魚、青魚、黃顙魚仔魚、林蛙蝌蚪、斜生柵藻)進(jìn)行急性毒性試驗(yàn)以補(bǔ)充毒性數(shù)據(jù),并利用二步外推法補(bǔ)充慢性毒性數(shù)據(jù)[19].對(duì)大型溞和青鳉進(jìn)行原水毒性試驗(yàn)以利用WER推導(dǎo)水質(zhì)基準(zhǔn).

    (1)急性試驗(yàn)試劑及生物.金屬銅(CuCl2) 標(biāo)準(zhǔn)溶液購(gòu)自百靈威,曝氣時(shí)間>3d的曝氣水作為稀釋水.其水質(zhì)參數(shù)為pH值(8.00±0.14)、堿度(81.25±2.17) mg(CaCO3)/L、硬度(125.50±4.95) mg(CaCO3)/L、總有機(jī)碳TOC 2.55mg(C)/L.試驗(yàn)所用魚和蝌蚪購(gòu)自江蘇溧水淡水漁業(yè)研究所,斜生柵藻購(gòu)自中科院武漢水生所,大型溞為實(shí)驗(yàn)室馴養(yǎng)生物.受試生物于實(shí)驗(yàn)室馴養(yǎng)兩周以上,期間每日喂食并定時(shí)換水,斜生柵藻則定期接種更新培養(yǎng)基并保證生物死亡率£10%.

    (2)急性毒性試驗(yàn).急性毒性試驗(yàn)依據(jù)OECD和U.S.EPA技術(shù)指南進(jìn)行.選擇大小相近、健康的個(gè)體進(jìn)行隨機(jī)分組,將暴露液配置為5個(gè)濃度梯度并設(shè)置稀釋水空白對(duì)照,暴露液每天更新一次,每濃度4 組平行,每平行10個(gè)生物.除大型溞和斜生柵藻暴露48h外,其余受試生物暴露時(shí)長(zhǎng)均為96h.實(shí)驗(yàn)條件按照國(guó)標(biāo)或OECD標(biāo)準(zhǔn).每24h觀察記錄生物死亡及異常情況,并將死亡個(gè)體及時(shí)取出,實(shí)驗(yàn)過(guò)程中不喂食.對(duì)于藻類則取對(duì)數(shù)生長(zhǎng)期的斜生柵藻,每24h用分光光度法測(cè)定藻液濃度.

    1.2.3 太湖銅的水質(zhì)基準(zhǔn)推導(dǎo) 利用WSSD、WER和BLM法推導(dǎo)太湖銅水質(zhì)基準(zhǔn).

    傳統(tǒng)物種敏感度分布法(SSD),通常利用Weibull法(式(1))計(jì)算累積概率[20].本文在SSD基礎(chǔ)上,利用WSSD法推導(dǎo)太湖銅的水質(zhì)基準(zhǔn).考慮生物類群間的能量傳遞,將水生生物分為植物(P)、脊椎動(dòng)物(V)及無(wú)脊椎動(dòng)物(I).根據(jù)對(duì)太湖21世紀(jì)的水生生物的調(diào)查統(tǒng)計(jì),太湖水生生物共524 種,植物、脊椎和無(wú)脊椎動(dòng)物的物種數(shù)分別為353、60和111 種,占比分別為67.4%、11.4%和21.2%[21-25].根據(jù)式(2)~(5)賦予權(quán)重以計(jì)算累積概率P(式(1)~(5)中,P指累積概率,為排序,為樣本量,為太湖物種數(shù),p、v和i分別為太湖中植物、脊椎和無(wú)脊椎動(dòng)物數(shù)量,、和分別為可獲得的太湖植物、脊椎和無(wú)脊椎動(dòng)物相應(yīng)的毒性數(shù)據(jù)量,式(5)中的P或v或I取決于排序的物種類別,擬合WSSD曲線,得出急慢性HC5值(影響5%物種的危害濃度),最后除以評(píng)價(jià)因子(AF=2)[26],得到太湖銅的急慢性基準(zhǔn)值分別以最大濃度基準(zhǔn)值(CMC)和持續(xù)濃度基準(zhǔn)值(CCC)表示.

    WER法基于經(jīng)驗(yàn)[27],推導(dǎo)因?qū)嶒?yàn)室與真實(shí)環(huán)境間的差異而導(dǎo)致的毒性比值WER,從而對(duì)實(shí)驗(yàn)數(shù)據(jù)推導(dǎo)的水質(zhì)基準(zhǔn)值進(jìn)行修正;BLM法則是從機(jī)制出發(fā)[28-29],反映多因子對(duì)生物有效性和毒性的影響.本文通過(guò)2 種本土生物(大型溞、青鳉)原位水和實(shí)驗(yàn)水的急性毒性試驗(yàn),利用WER方法獲得WER實(shí)際校正后的銅基準(zhǔn)值.測(cè)定陳東港等5個(gè)太湖代表性點(diǎn)位BLM模型所需參數(shù):pH值、溫度、溶解性有機(jī)碳(DOC)、主要陽(yáng)離子(Ca2+,Mg2+,Na+和K+)、陰離子(SO42-和Cl-)和堿度.以BLM模型中WQC計(jì)算模塊得到BLM法太湖銅的CMC及CCC值.

    1.3 太湖銅的水生生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    聯(lián)合概率法(JPC)能夠?qū)⒍拘詳?shù)據(jù)和化合物的環(huán)境濃度相結(jié)合[30-31],反映環(huán)境濃度超過(guò)影響一定百分比物種的概率.本文利用WER和BLM模型法得出兩種水效應(yīng)比WER實(shí)際和WERBLM,以校正毒性數(shù)據(jù),結(jié)合太湖不同水期52個(gè)點(diǎn)位Cu2+的暴露濃度,通過(guò)JPC反映太湖銅水生生態(tài)風(fēng)險(xiǎn).在水生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)中,通常以保護(hù)95%的水生生物作為保護(hù)目標(biāo)[30],故本文將太湖銅的水生生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表征為“太湖銅濃度超過(guò)影響5%生物的概率”.

    1.4 統(tǒng)計(jì)學(xué)分析

    補(bǔ)充的急性毒性數(shù)據(jù)LC50及95%的置信區(qū)間以GraphPad Prism 8中的Sigmoidal dose-response模型獲得.慢性數(shù)據(jù)運(yùn)用Microsoft Office Excel 2019外推獲得,最終急慢性毒性數(shù)據(jù)采用Kolmogorov- Smirnov進(jìn)行正態(tài)性檢驗(yàn),SSD及WSSD曲線均采用GraphPad Prism 8中l(wèi)og-Normal方法擬合,以R2評(píng)價(jià)擬合優(yōu)度.聯(lián)合概率曲線由GraphPad Prism 8軟件繪制,BLM由Windward Environmental提供(http: //www.windwardenv.com/biotic-ligand-model).

    2 結(jié)果與討論

    2.1 急慢性數(shù)據(jù)補(bǔ)充結(jié)果

    補(bǔ)充的急慢性毒性數(shù)據(jù)如表1.對(duì)銅最敏感的水生生物為斜生柵藻,最耐受物種為青鳉等脊椎動(dòng)物.結(jié)合搜集的毒性數(shù)據(jù),最終銅的急慢性毒性數(shù)據(jù)見(jiàn)表2.共獲得銅的22種太湖水生生物急性毒性數(shù)據(jù),包括5門14 科;16種太湖水生生物慢性毒性數(shù)據(jù),分屬于5門11科,滿足保護(hù)水生生物水質(zhì)基準(zhǔn)推導(dǎo)要求[32].

    2.2 WSSD推導(dǎo)太湖銅水質(zhì)基準(zhǔn)結(jié)果

    由銅對(duì)水生生物的急慢性毒性數(shù)據(jù),得到未考慮因子影響的太湖銅急慢性傳統(tǒng)SSD曲線如圖2所示.整體而言,對(duì)于銅的毒性,藻類是最敏感生物,其次是無(wú)脊椎動(dòng)物,最不敏感類群為魚類.根據(jù)傳統(tǒng)SSD,太湖銅的急慢性HC5值分別為5.70μg/L和3.44μg/L.考慮太湖水生生物組成,得到太湖銅急慢性物種權(quán)重敏感度分布曲線WSSD如圖3.WSSD與傳統(tǒng)SSD法識(shí)別出的太湖銅敏感物種與耐受物種一致,但閾值上有一定差異.

    通過(guò)SSD和WSSD獲得的HC5值,推導(dǎo)出未考慮(SSD)和考慮生物影響(WSSD)的太湖銅水質(zhì)基準(zhǔn)(WQC)如表3.其中,WSSD推導(dǎo)的CMC和CCC分別為1.43μg/L和1.33μg/L,急性基準(zhǔn)值為傳統(tǒng)SSD推導(dǎo)值CMC 2.85μg/L的一半,慢性基準(zhǔn)值也較SSD的CCC 1.72μg/L低,可見(jiàn)考慮物種組成的WSSD比傳統(tǒng)SSD法對(duì)特定區(qū)域水生生物的保護(hù)更為嚴(yán)格.本文的基準(zhǔn)值與Shi等[33]利用WSSD推導(dǎo)的CMC 5.3μg/L和CCC 3.7μg/L略有差異,可能是由于補(bǔ)充了敏感物種藻類等生物的毒性數(shù)據(jù)以及在賦予權(quán)重時(shí)對(duì)太湖物種組成調(diào)研時(shí)間不同造成的.

    表1 太湖6種本土水生生物的銅急性毒性數(shù)據(jù)及預(yù)測(cè)慢性毒性數(shù)據(jù)

    表2 太湖水生生物銅的急慢性毒性數(shù)據(jù)

    續(xù)表2

    注:*為本研究補(bǔ)充的數(shù)據(jù).

    我國(guó)目前的水質(zhì)標(biāo)準(zhǔn)主要參考美國(guó)等發(fā)達(dá)國(guó)家的水質(zhì)標(biāo)準(zhǔn)[34],使用最多的水質(zhì)基準(zhǔn)推導(dǎo)方法是SSD法,而我國(guó)水域眾多,與其他國(guó)家水生物種差異大,國(guó)內(nèi)不同水域中物種及生態(tài)狀況也有差異.相關(guān)研究表明不同物種對(duì)銅的敏感性存在顯著性差異,且相較而言無(wú)脊椎動(dòng)物對(duì)銅的敏感度更高[35],與本文的結(jié)論一致.此外,國(guó)內(nèi)不同水域的水質(zhì)基準(zhǔn)差異顯著,甚至可達(dá)6倍[36].由此可見(jiàn),在推導(dǎo)特定水域的水質(zhì)基準(zhǔn)時(shí),考慮物種的影響,結(jié)合研究區(qū)域特定的本土生物是非常必要的,WSSD法比傳統(tǒng)SSD更能真實(shí)地反映研究區(qū)域化合物對(duì)水生生物潛在的影響,體現(xiàn)同一化合物對(duì)不同區(qū)系水生生物的差異,從而針對(duì)性地支撐特定區(qū)域水生生物的保護(hù).

    表3 未考慮和考慮太湖水生生物組成的銅基準(zhǔn)閾值

    2.3 WER和BLM法推導(dǎo)太湖銅水質(zhì)基準(zhǔn)結(jié)果

    2.3.1 WER法推導(dǎo)結(jié)果 根據(jù)太湖2種本土生物(大型溞、青鳉)的3個(gè)點(diǎn)位的原位水及實(shí)驗(yàn)水急性毒性試驗(yàn)結(jié)果,得出實(shí)際水效應(yīng)比WER實(shí)際為 2.55 (表4),最終推導(dǎo)出WER法的太湖銅水質(zhì)基準(zhǔn)CMC和CCC分別為7.27μg/L和4.39μg/L.

    表4 WER實(shí)際的推導(dǎo)結(jié)果

    2.3.2 BLM推導(dǎo)結(jié)果 BLM模型輸入?yún)?shù)如表5,包括點(diǎn)位(陳東港CD、大浦港DP、分水大橋FSD、分水橋FS和平臺(tái)山PT和實(shí)驗(yàn)室LAB)的水樣溫度、pH值、Cu2+濃度、有機(jī)物、主要陰陽(yáng)離子含量及堿度.其中,HA和缺省值2-分別采用推薦值10%和10-10mg/L[37-38],主要陰離子則采用2008~ 2014年報(bào)道的濃度均值[29,39],得到BLM模型法推導(dǎo)太湖銅的瞬時(shí)水質(zhì)基準(zhǔn)值及模型水效應(yīng)比見(jiàn)表6,CMC和CCC分別為42.57μg/L和26.44μg/L,結(jié)果與Zhang等[40]利用BLM模型得出的太湖水質(zhì)基準(zhǔn)CMC 32.2μg/L近似.BLM模型較為全面的考慮了影響生物有效性的因素(競(jìng)爭(zhēng)、吸附解析等)[41],能夠一定程度上反映不同區(qū)域的差異.Chen等利用BLM模型得到瀾滄江上中下游銅的急性水質(zhì)基準(zhǔn)CMC分別為22.39, 31.20和26.79μg/L[42].可見(jiàn),對(duì)于不同水環(huán)境的基準(zhǔn)值有所差異,BLM法應(yīng)用于水質(zhì)基準(zhǔn)中,能夠?qū)⑺h(huán)境因子對(duì)基準(zhǔn)的影響納入考慮,在結(jié)果中體現(xiàn)出時(shí)空差異性.

    表5 BLM模型輸入?yún)?shù)

    注:點(diǎn)位見(jiàn)2.3.2節(jié).

    WER法和BLM模型法所得水效應(yīng)比WER實(shí)際2.55和WERBLM2.0值相近,急慢性基準(zhǔn)值相差6倍以內(nèi),此差異是由于兩種方法特點(diǎn)、內(nèi)在機(jī)制不同所造成的.兩種方法從不同層面鏈接了野外環(huán)境,在區(qū)域性水質(zhì)基準(zhǔn)的推導(dǎo)中,具有較強(qiáng)的針對(duì)性.

    由于水域污染日益嚴(yán)重,水體富營(yíng)養(yǎng)化加劇,目前太湖生態(tài)狀況正在惡化.根據(jù)調(diào)查結(jié)果,太湖水生物種數(shù)大量減少,其中魚類和大型浮游動(dòng)物數(shù)量急劇下降,耐污物種也在增加[21,43].因此,基于WSSD、WER和BLM三種方法結(jié)果,為最大程度保護(hù)太湖水生生物的安全,本文推薦采用WSSD法推導(dǎo)的CMC 1.43μg/L和CCC 1.33μg/L作為保護(hù)太湖水生生物銅的水質(zhì)基準(zhǔn).WER和BLM法則可根據(jù)太湖污染及生態(tài)狀況的變化,進(jìn)行互補(bǔ)性應(yīng)用.

    表6 BLM得出的水質(zhì)基準(zhǔn)及模型水效應(yīng)比結(jié)果

    2.4 基于多因子的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

    利用WERBLM和WER實(shí)際校正后的銅慢性毒性數(shù)據(jù)及枯豐水期太湖銅的濃度分別見(jiàn)表7和表8.

    表7 WER校正后的太湖銅慢性毒性

    結(jié)合校正后的毒性數(shù)據(jù),得到JPC曲線如圖4.對(duì)于太湖豐水期(圖4,(a)),未考慮多因子影響的評(píng)估結(jié)果為“太湖銅對(duì)5%水生生物造成慢性危害的概率為85.01%”,而通過(guò)兩種WER校正后的評(píng)估風(fēng)險(xiǎn)概率分別為39.43%和23.43%,遠(yuǎn)低于傳統(tǒng)方法的評(píng)估結(jié)果,枯水期(圖4,(b))的結(jié)果類似.此外,相比傳統(tǒng)方法,考慮多因子的風(fēng)險(xiǎn)評(píng)估識(shí)別出銅在豐水期的風(fēng)險(xiǎn)明顯高于枯水期,體現(xiàn)了太湖銅生態(tài)風(fēng)險(xiǎn)的時(shí)間差異性.由此可見(jiàn),在評(píng)估太湖銅的水生生態(tài)風(fēng)險(xiǎn)時(shí),忽略水化學(xué)等多因子的影響,很可能會(huì)使得評(píng)估風(fēng)險(xiǎn)偏高,評(píng)估不確定性大,且“掩蓋”風(fēng)險(xiǎn)在特定區(qū)域的時(shí)間差異性,而鏈接研究區(qū)域的真實(shí)環(huán)境,能夠更為客觀地評(píng)估化合物風(fēng)險(xiǎn),指示不同時(shí)間特定區(qū)域的生態(tài)風(fēng)險(xiǎn)狀況,更好地支撐化合物的管控以及水生生態(tài)的保護(hù).

    表8 太湖52 個(gè)點(diǎn)位枯豐水期銅離子濃度

    3 結(jié)論

    3.1 WSSD,WER和BLM模型3 種方法得出環(huán)境多因子影響下,保護(hù)太湖水生生物銅的CMC分別為1.43、7.27和42.57μg/L,CCC分別為1.33、4.39和26.44μg/L.根據(jù)太湖現(xiàn)狀,為最大程度保護(hù)太湖水生生物安全,本文推薦采用WSSD法推導(dǎo)的CMC 1.43μg/L和CCC 1.33μg/L作為保護(hù)太湖水生生物銅的水質(zhì)基準(zhǔn).

    3.2 結(jié)合WER和BLM法,得出太湖豐水期銅的慢性生態(tài)風(fēng)險(xiǎn)分別為23.43%和39.43%,枯水期分別為0.64%和5.40%,傳統(tǒng)方法則高估了風(fēng)險(xiǎn)(85.01%和82.09%),太湖銅的水生生物風(fēng)險(xiǎn)不容忽視,尤其需要關(guān)注豐水季節(jié).

    3.3 在水質(zhì)基準(zhǔn)和生態(tài)風(fēng)險(xiǎn)評(píng)估中,傳統(tǒng)方法會(huì)導(dǎo)致結(jié)果的偏高,考慮區(qū)域性生物組成和水環(huán)境因子的影響十分必要.

    3.4 WSSD,WER和BLM法各有其應(yīng)用要求與特點(diǎn),對(duì)于我國(guó)重金屬類保護(hù)水生生物水質(zhì)基準(zhǔn)的制定及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),需根據(jù)實(shí)際情況進(jìn)行方法的擇優(yōu)或互補(bǔ)性應(yīng)用,達(dá)到“因地制宜”式保護(hù)水生生物的目標(biāo).

    [1] Chain F J J, Finlayson S, Crease T, et al. Variation in transcriptional responses to copper exposure across Daphnia pulex lineages [J]. Aquatic Toxicology, 2019,210:85-97.

    [2] Prato E, Parlapiano I, Biandolino F, et al. Sublethal effects of copper on some biological traits of the amphipod Gammarus aequicauda reared under laboratory conditions [J]. Chemosphere, 2013,93(6): 1015-1022.

    [3] Rodgher S, Lombardi A T, Mel?o M G, et al. Change in life cycle parameters and feeding rate of(Crustacea, Cladocera) exposure to dietary copper [J]. Ecotoxicology, 2008,17(8):826-833.

    [4] Donnachie R L, Johnson A C, Moeckel C, et al. Using risk-ranking of metals to identify which poses the greatest threat to freshwater organisms in the UK [J]. Environmental Pollution, 2014,194:17-23.

    [5] 廖 偉,劉大慶,馮承蓮,等.不同生長(zhǎng)階段斑馬魚對(duì)Cu2+的毒性響應(yīng)差異 [J]. 環(huán)境科學(xué)研究, 2020,33(3):626-633.

    Liao W, Liu D Q, Feng C L, et al. Difference in toxicity response of Zebrafish to Cu2+at different life stages [J].Research of Environmental Sciences, 2020,33(3):626-633.

    [6] 王 振,金小偉,王子健,等.銅對(duì)水生生物的毒性:類群特異性敏感度分析 [J]. 生態(tài)毒理學(xué)報(bào), 2014,9(4):640-646.

    Wang Z, Jin X W, Wang Z J, et al. Taxon-specific sensitivity differences of copper to aquatic organisms [J]. Asian Journal of Ecotoxicology, 2014,9(4):640-646.

    [7] Rogevich E C, Hoang T C, Rand G M, et al. The effects of water quality and age on the acute toxicity of copper to the Florida applesnail (Pomacea paludosa) [J]. Archives of Environmental Contamination and Toxicology, 2008,54(4):690-696.

    [8] Lock K, Van E H, Karel A C, et al. Development of a biotic ligand model (BLM) predicting nickel toxicity to barley () [J]. Chemosphere, 2007,66(7):1346-1352.

    [9] Carlson A R. Guidelines for deriving numerical aquatic site- specific water quality criteria bymodifying national criteria [M]. Environmental Research Laboratory, Office of Research and Development, US Environmental Protection Agency, 1984.

    [10] Chen Y, Yu S, Tang S, et al. Site-specific water quality criteria for aquatic ecosystems: A case study of pentachlorophenol for Tai Lake, China [J]. Science of The Total Environment, 2016,541:65-73.

    [11] Feng C L, Wu F C, Zheng B H, et al. Biotic ligand models for metals-a practical application in the revision of water quality standards in China [J]. Environmental Science & Technology, 2012, 46(20):10877-10878.

    [12] 孫雪華,孫 成,劉紅玲.考慮物種權(quán)重校驗(yàn)保護(hù)太湖水生生物的鉛基準(zhǔn)[J]. 環(huán)境化學(xué), 2020,39(6):1578-1589.

    Sun X H, Sun C, Liu H L. Weighted species sensitivity distribution method to derive site-specific quality criteria of lead for protection of aquatic life in Tai Lake [J]. Environmental Chemistry, 2020,39(6): 1578-1589.

    [13] Arnold W R, Warren-Hicks W J. Assessment of aquatic ecological risk and site-specific criteria of copper in san francisco bay, California, USA [J]. Integrated Environmental Assessment and Management, 2007,3(1):32-48.

    [14] 蘇海磊,李信茹,陶艷茹,等.美國(guó)水質(zhì)標(biāo)準(zhǔn)制定研究及其對(duì)中國(guó)的借鑒意義 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2021,30(11):2267-2274.

    Su H L, Li X R, Tao Y R, et al. Study on water quality standard formulation of USA and its revelation to China [J]. Ecology and Environmental Sciences, 2021,30(11):2267-2274.

    [15] 周巧巧,任 勃,李有志,等.中國(guó)河湖水體重金屬污染趨勢(shì)及來(lái)源解析 [J]. 環(huán)境化學(xué), 2020,39(8):2044-2054.

    Zhou Q Q, Ren B, Li Y Z, et al. Trends and sources of dissolved heavy metal pollution in water of rivers and lakes in China [J]. Environmental Chemistry, 2020,39(8):2044-2054.

    [16] 齊文啟,連 軍,孫宗光,等.《地表水和污水監(jiān)測(cè)技術(shù)規(guī)范》(HJ/T91-2002)的相關(guān)技術(shù)說(shuō)明 [J]. 中國(guó)環(huán)境監(jiān)測(cè), 2006,(1):54-57.

    Qi W Q, Lian J, Sun C G, et al. Technical explanations for “Technical specifications requirements for monitoring of surface water and waste water” (HJ/T 91- 2002) [J]. Environmental Monitoring in China, 2006, (1):54-57.

    [17] 國(guó)家環(huán)境保護(hù)局科技標(biāo)準(zhǔn)司.水質(zhì)-物質(zhì)對(duì)淡水魚(斑馬魚)急性毒性測(cè)定方法(GB/T 13267-1991) [S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社, 1992.

    MEPPRC. Water quality-determination of the acute toxicity of substance to freshwater fish (Brachydanio rerio Hamilton-Buchanan) [S]. Beijing: China Standards Press, 1992.

    [18] OECD. OECD guideline for testing of chemicals-sp., acute immobilisation test. Guideline [R]. Paris, France: OECD Environment Directorate, 2004.

    [19] Giesy J P, Graney R L. Recent developments in and intercomparisons of acute and chronic bioassays and bioindicators [J]. Hydrobiologia, 1989,188(1):21-60.

    [20] 王業(yè)耀,張鈴松,孟凡生,等.水生生物水質(zhì)基準(zhǔn)研究進(jìn)展及建立我國(guó)氨氮水質(zhì)基準(zhǔn)的探討[J]. 南水北調(diào)與水利科技, 2012,10(5):108-113.

    Wang Y Y, Zhang L S, Meng F S, et al. Research progress on water quality criteria for aquatic life and establishment of water quality criteria of ammonia nitrogen in China [J]. South-to-North Water Transfer and Water Science & Technology, 2012,10(5):108-113.

    [21] 陳立僑,劉 影,楊再福,等.太湖生態(tài)系統(tǒng)的演變與可持續(xù)發(fā)展[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003,(4):99-106.

    Chen L J, Liu Y, Yang Z F, et al. Ecological succession and sustainables development in Tai Lake [J]. Journal of East China Normal University (Natural Science), 2003,(4):99-106.

    [22] Cai Y J, Gong Z J, Qin B Q, et al. Community structure and diversity of macrozoobenthos in Tai Lake, a large shallow eutrophic lake in China [J] . Biodiversity Science, 2010,18(1):50-59.

    [23] 鄭丙輝,田自強(qiáng),張 雷,等.太湖西岸湖濱帶水生生物分布特征及水質(zhì)營(yíng)養(yǎng)狀況[J]. 生態(tài)學(xué)報(bào), 2007,(10):4214-4223.

    Zhen B H, Tian Z Q, Zhang L, et al. The characteristics of the Hydrobios' distribution and the analysis of water quality along the west shore of Tai Lake [J]. Acta Ecologica Sinica, 2007,27(10):4214-4223.

    [24] Zhu Y Q. Ichtyoloical survey of Tai Lake during 2002~2003 [J]. Journal of Lake Sciences, 2004,(2):120-124.

    [25] Xu H, Cai Y J, Tang X M, et al. Community structure of macrozoobenthos and the evaluation of water environment in Tai Lake [J]. Journal of Lake Sciences, 2015,27(5):840-852.

    [26] Kooijman S A L M. A safety factor for LC50 values allowing for differences in sensitivity among species [J]. Water Research, 1987, 21(3):269-276.

    [27] Welsh P G, Lipton J, Chapman G A, et al. Evaluation of water-effect ratio methodology for establishing site-specific water quality criteria [J]. Environmental Toxicology and Chemistry, 2000,19(6):1616-1623.

    [28] Paquin P R, Gorsuch J W, Apte S, et al. The biotic ligand model: A historical overview [J]. Comparative Biochemistry and Physiology C-Toxicology & Pharmacology, 2002,133(1/2):3-35.

    [29] Xie R L, Zhou Q X. Methodology of water quality criteria in foreign countries and its future researching prospects [J]. World Sci-Tech R&D, 2012,34(6):939-944.

    [30] Solomon K, Giesy J, Jones P. Probabilistic risk assessment of agrochemicals in the environment [J]. Crop Protection, 2000,19(8-10): 649-655.

    [31] Wang X L, Xu F L, Li B G, et al. Risk assessment of ecotoxicity of benzo(a)pyrene, fluorethene and phenanthrene in Tianjin wastewater irrigated area [J]. Urban Environment and Urban Ecology, 2002, 15(4):10-12.

    [32] 馮承蓮,吳豐昌,趙曉麗,等.水質(zhì)基準(zhǔn)研究與進(jìn)展 [J]. 中國(guó)科學(xué):地球科學(xué), 2012,42(5):646-656.

    Feng C L, Wu F C, Zhao X L, et al. Research and development of quality criterion [J]. Scientia Sinica (Terrae), 2012,42(5):646-656.

    [33] Shi R, Yang C, Su R, et al. Weighted species sensitivity distribution method to derive site-specific quality criteria for copper in Tai Lake, China [J]. Environmental Science and Pollution Research, 2014,21 (22):12968-12978.

    [34] 吳豐昌,馮承蓮,張瑞卿,等.我國(guó)典型污染物水質(zhì)基準(zhǔn)研究 [J]. 中國(guó)科學(xué):地球科學(xué), 2012,42(5):665-672.

    Wu F C, Wu F C, Zhang R Q, et al. Research on water quality benchmark of typical pollutants in my country [J]. Scientia Sinica (Terrae), 2012,42(5):665-672.

    [35] 孔祥臻,何 偉,秦 寧,等.重金屬對(duì)淡水生物生態(tài)風(fēng)險(xiǎn)的物種敏感性分布評(píng)估 [J]. 中國(guó)環(huán)境科學(xué), 2011,31(9):1555-1562.

    Kong X Z, He W, Qin N T, et al. Species sensitivity analysis of heavy metals to freshwater organisms [J]. China Environmental Science, 2011,31(9):1555-1562.

    [36] Wang Y Z, Yan Z G, Zhan X, et al. Preliminary aquatic life criteria development and ecological risk assessment of ammonia in seven major basins in China [J]. Research of Environmental Sciences, 2016, 29(1):77-83.

    [37] The Biotic Ligand Model Windows Interface, Version 2.2.3: User’s Guide and Reference Manual [M]. HydroQual, Inc, Mahwah, NJ, April 2005.

    [38] Gondek J C, Gensemer R W, Claytor C A, et al. Framework for derivation of water quality criteria using the biotic ligand model: Copper as a case study [J]. Integrated Environmental Assessment and Management, 2018,14(6):736-749.

    [39] Ye H M, Yuan X Y, Ge M X, et al. Water chemistry characteristics and controlling factors in the northern rivers in the Taihu Basin [J]. Ecology and Environmental Sciences, 2010,19(1):23-27.

    [40] Zhang Y H, Zang W C, Qin L M, et al. Water quality criteria for copper based on the BLM approach in the freshwater in China [J]. Plos One, 2017,12(2):e0170105.

    [41] Wang C Y, Chen H, An L H, et al. An updated review on biotic ligand model in predicting metal bioavailability in surface waters [J]. Enuivonmental Science and Technology, 2011,34(8):75-80.

    [42] 陳 莎.瀾滄江銅的水質(zhì)基準(zhǔn)與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究 [D]. 昆明:昆明理工大學(xué), 2014.

    Chen S. Research on water quality benchmark and ecological risk assessment of copper in Lancang River [D]. Kunming: Kunming University of Science and Technology, 2014.

    [43] Fan L X. Historical evolution of water ecological setting in Tai Lake [J]. Journal of Lake Sciences, 1996,8(4):297-304.

    致謝:本研究的采樣測(cè)定工作由王遵堯老師協(xié)助完成,在此表示感謝.

    The water quality criteria and ecological risks of copper under the influence of multiple factors.

    LIU Na1, LI Ya-bing1, LIU Hong-ling1,2*

    (1.State Key Laborotory of Pollution Control & Resource Reuse, School of the Environment, Nanjing University, Nanjing 210023, China;2.Jiangsu Province Ecology and Environment Protection Key Laboratory of Chemical Safety and Health Risk, Nanjing 210023, China)., 2022,42(7):3353~3361

    In order to explore the impact of multiple factors on water quality criteria and ecological risks, the weighted species sensitivity distribution (WSSD), water effect ratio (WER) and biological ligand model (BLM) were used to derive the water quality criteria of copper for protecting aquatic organisms in Tai Lake. The criterion maximum concentration (CMC) of 1.43μg/L and the criterion continuous concentration (CCC) of 1.33μg/L were recommended as the water quality criteria of copper in Tai Lake. Joint probability curve (JPC) combining WER and BLM were used to further assess the ecological risks. There were 23.43% (WER) and 39.43% (BLM) probability of affecting 5% aquatic organisms during wet season in Tai Lake, respectively, while the traditional method overestimated the risks of copper in Tai Lake with a probability of 85.01%. These findings clearly suggest that the influence of multiple factors in the aquatic environment on water quality criteria and ecological risks cannot be ignored, currently environmental quality standard of copper adopted in China may not be able to protect aquatic organisms in specific areas such as Tai Lake. Considering multiple factors will obtain more reasonable water quality criteria and ecological risks assessment results, and avoid the phenomenon of “over-protection” or “under-protection”.

    copper;influence of multiple factors;water quality criteria;ecological risks;weighted species sensitivity distribution (WSSD);water effect ratio(WER);biological ligand model(BLM)

    X524

    A

    1000-6923(2022)07-3353-09

    劉 娜(1996-),女,江西吉安人,南京大學(xué)碩士研究生,研究方向?yàn)樗|(zhì)基準(zhǔn)及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).

    2021-12-10

    國(guó)家自然科學(xué)基金資助項(xiàng)目(22176095,21677073);國(guó)家科技重大專項(xiàng)(2018ZX07208001);國(guó)家重點(diǎn)研發(fā)項(xiàng)目(2018YFC1801505)

    * 責(zé)任作者, 副教授,hlliu@nju.edu.cn

    猜你喜歡
    多因子基準(zhǔn)值水生
    河北省啤酒行業(yè)清潔生產(chǎn)水平分析
    基于5G用戶體驗(yàn)的業(yè)務(wù)質(zhì)量?jī)?yōu)化模型研究及其應(yīng)用
    基于BP神經(jīng)網(wǎng)絡(luò)的多因子洪水分類研究
    基于打分法的多因子量化選股策略研究
    一種基于改進(jìn)差分的測(cè)井?dāng)?shù)據(jù)可逆變長(zhǎng)碼壓縮方法
    基于多因子的ZigBee安全認(rèn)證機(jī)制
    三氯生對(duì)4種水生生物的急性毒性研究
    依托科技創(chuàng)新 打造現(xiàn)代水生蔬菜產(chǎn)業(yè)
    跟我進(jìn)城的兄弟
    椰城(2013年4期)2013-11-16 06:31:11
    “多因子 全方位”加強(qiáng)對(duì)高校家庭經(jīng)濟(jì)困難學(xué)生引導(dǎo)
    一边亲一边摸免费视频| 国精品久久久久久国模美| 亚洲欧美中文字幕日韩二区| 亚洲国产精品国产精品| 黄色配什么色好看| 这个男人来自地球电影免费观看 | 久久韩国三级中文字幕| 日本欧美国产在线视频| 男女无遮挡免费网站观看| 精品国产露脸久久av麻豆| 久久久久精品久久久久真实原创| 黄片播放在线免费| 97精品久久久久久久久久精品| 久久精品国产a三级三级三级| 在线天堂中文资源库| 国产黄频视频在线观看| 黄色怎么调成土黄色| 亚洲一区中文字幕在线| 国产日韩欧美在线精品| 国产成人精品无人区| 一本色道久久久久久精品综合| 亚洲成人手机| videos熟女内射| 亚洲成人av在线免费| 日本色播在线视频| 国产成人a∨麻豆精品| 国产精品蜜桃在线观看| 捣出白浆h1v1| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 韩国高清视频一区二区三区| 狠狠精品人妻久久久久久综合| 国产毛片在线视频| 国产精品一国产av| 亚洲精品av麻豆狂野| 国产成人欧美| 国产精品一区二区在线观看99| 欧美bdsm另类| 亚洲国产最新在线播放| 国产成人精品一,二区| 精品国产露脸久久av麻豆| 色婷婷av一区二区三区视频| 中文字幕精品免费在线观看视频| 精品视频人人做人人爽| 日韩人妻精品一区2区三区| 国产激情久久老熟女| 国产片特级美女逼逼视频| 亚洲欧美色中文字幕在线| 免费女性裸体啪啪无遮挡网站| av片东京热男人的天堂| www日本在线高清视频| 在线观看人妻少妇| 边亲边吃奶的免费视频| 青春草亚洲视频在线观看| 亚洲五月色婷婷综合| 两个人免费观看高清视频| 永久免费av网站大全| 啦啦啦在线观看免费高清www| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 日本免费在线观看一区| 七月丁香在线播放| 亚洲国产色片| 久久久欧美国产精品| 这个男人来自地球电影免费观看 | 欧美成人午夜精品| 亚洲久久久国产精品| 如日韩欧美国产精品一区二区三区| 国产一区二区 视频在线| 日韩中文字幕欧美一区二区 | 久久婷婷青草| 久久99精品国语久久久| 久久 成人 亚洲| 男人舔女人的私密视频| 日本-黄色视频高清免费观看| 18禁裸乳无遮挡动漫免费视频| 夜夜骑夜夜射夜夜干| 午夜福利视频在线观看免费| 边亲边吃奶的免费视频| 青青草视频在线视频观看| 成人二区视频| 亚洲人成77777在线视频| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 日本91视频免费播放| 一级毛片黄色毛片免费观看视频| 国产精品三级大全| 另类精品久久| 老司机影院毛片| 亚洲av.av天堂| 日韩免费高清中文字幕av| 少妇精品久久久久久久| 亚洲精品视频女| 在线观看国产h片| 一区二区三区激情视频| 久久久久久久久久久免费av| 少妇的丰满在线观看| 日本欧美国产在线视频| 日韩在线高清观看一区二区三区| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 一边摸一边做爽爽视频免费| 欧美日本中文国产一区发布| 极品人妻少妇av视频| 亚洲精品国产av成人精品| 久久久久久久亚洲中文字幕| 18在线观看网站| 观看美女的网站| 晚上一个人看的免费电影| 久久久久久久久久久免费av| 日韩av免费高清视频| 国产黄色视频一区二区在线观看| 亚洲国产日韩一区二区| 国产一区亚洲一区在线观看| 麻豆乱淫一区二区| 老司机影院毛片| 大码成人一级视频| 少妇人妻精品综合一区二区| 国产成人精品福利久久| 久久久久国产一级毛片高清牌| 亚洲一区中文字幕在线| 不卡视频在线观看欧美| 中国三级夫妇交换| 亚洲精品第二区| 日本欧美国产在线视频| 久久狼人影院| 国产 一区精品| 亚洲成人手机| 天天影视国产精品| 日韩一本色道免费dvd| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| 免费观看av网站的网址| 丰满迷人的少妇在线观看| 女人被躁到高潮嗷嗷叫费观| 人人妻人人澡人人爽人人夜夜| 一级毛片 在线播放| 国产老妇伦熟女老妇高清| 国产精品国产av在线观看| 亚洲熟女精品中文字幕| 超碰成人久久| 亚洲男人天堂网一区| 最新的欧美精品一区二区| 亚洲av.av天堂| 国产片特级美女逼逼视频| 婷婷成人精品国产| 亚洲国产日韩一区二区| 高清av免费在线| 国产精品一二三区在线看| 五月开心婷婷网| 18禁裸乳无遮挡动漫免费视频| 国产成人aa在线观看| kizo精华| 精品一区二区免费观看| 97精品久久久久久久久久精品| 老汉色av国产亚洲站长工具| 自拍欧美九色日韩亚洲蝌蚪91| 99热全是精品| 国产有黄有色有爽视频| 久久人人爽av亚洲精品天堂| 国产日韩欧美亚洲二区| 黄片无遮挡物在线观看| 午夜老司机福利剧场| 男女无遮挡免费网站观看| 中文欧美无线码| 亚洲欧美一区二区三区久久| 国产精品国产av在线观看| 看非洲黑人一级黄片| 中文字幕精品免费在线观看视频| 人妻一区二区av| 亚洲美女视频黄频| 麻豆精品久久久久久蜜桃| 国产白丝娇喘喷水9色精品| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产专区5o| 看十八女毛片水多多多| 国产亚洲最大av| 纯流量卡能插随身wifi吗| 亚洲第一区二区三区不卡| 久久ye,这里只有精品| 免费黄网站久久成人精品| 天堂中文最新版在线下载| 国产人伦9x9x在线观看 | av网站在线播放免费| 妹子高潮喷水视频| 国产av精品麻豆| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 国产成人精品久久二区二区91 | 99久久精品国产国产毛片| 熟女av电影| 日本wwww免费看| 人妻少妇偷人精品九色| 美女脱内裤让男人舔精品视频| 丰满乱子伦码专区| 自线自在国产av| 亚洲一码二码三码区别大吗| 高清不卡的av网站| av视频免费观看在线观看| 少妇被粗大猛烈的视频| 亚洲五月色婷婷综合| 夫妻性生交免费视频一级片| 国产综合精华液| 久久久国产欧美日韩av| 欧美日韩精品成人综合77777| 欧美人与性动交α欧美软件| 国产在线免费精品| 免费av中文字幕在线| 在线亚洲精品国产二区图片欧美| 纵有疾风起免费观看全集完整版| 久久精品久久精品一区二区三区| 在线观看免费视频网站a站| 日本午夜av视频| 老司机亚洲免费影院| 亚洲久久久国产精品| 大码成人一级视频| 精品亚洲乱码少妇综合久久| 日韩av不卡免费在线播放| 国产成人a∨麻豆精品| 久热这里只有精品99| 26uuu在线亚洲综合色| 久久久久精品性色| 精品国产露脸久久av麻豆| 美女xxoo啪啪120秒动态图| 日日啪夜夜爽| freevideosex欧美| 超碰成人久久| 亚洲人成77777在线视频| av天堂久久9| 久久久国产欧美日韩av| 又大又黄又爽视频免费| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 国产免费又黄又爽又色| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 亚洲精品自拍成人| 亚洲av免费高清在线观看| 国产在线视频一区二区| 中国国产av一级| 国产野战对白在线观看| 久久精品久久久久久噜噜老黄| 人妻少妇偷人精品九色| 午夜福利在线观看免费完整高清在| 人人澡人人妻人| 丁香六月天网| 纵有疾风起免费观看全集完整版| 中文字幕最新亚洲高清| 热re99久久精品国产66热6| 深夜精品福利| 国产极品粉嫩免费观看在线| 国精品久久久久久国模美| 亚洲少妇的诱惑av| 成人毛片60女人毛片免费| 丰满乱子伦码专区| 免费看av在线观看网站| 十八禁网站网址无遮挡| 亚洲欧洲国产日韩| 人妻人人澡人人爽人人| 免费少妇av软件| 久久久亚洲精品成人影院| 丰满迷人的少妇在线观看| 国产精品蜜桃在线观看| 国产男女超爽视频在线观看| 女人精品久久久久毛片| 美女国产高潮福利片在线看| 日韩在线高清观看一区二区三区| 日本色播在线视频| 天天躁日日躁夜夜躁夜夜| 欧美97在线视频| 国产黄色免费在线视频| 亚洲精品久久午夜乱码| 欧美97在线视频| 精品一区二区免费观看| 飞空精品影院首页| 18在线观看网站| 亚洲精品日本国产第一区| 国产成人精品在线电影| 亚洲色图 男人天堂 中文字幕| 老司机亚洲免费影院| 欧美最新免费一区二区三区| 午夜福利视频精品| 欧美激情极品国产一区二区三区| 久久综合国产亚洲精品| 欧美激情 高清一区二区三区| 在线观看三级黄色| 欧美精品高潮呻吟av久久| 亚洲欧美成人综合另类久久久| 欧美成人午夜精品| freevideosex欧美| 黄色视频在线播放观看不卡| 亚洲伊人色综图| 午夜福利在线免费观看网站| 成年女人毛片免费观看观看9 | 亚洲精品第二区| 免费在线观看完整版高清| av福利片在线| 午夜老司机福利剧场| 在线观看免费视频网站a站| 街头女战士在线观看网站| 卡戴珊不雅视频在线播放| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品古装| 看十八女毛片水多多多| 爱豆传媒免费全集在线观看| 日韩一区二区三区影片| 久久久久久久久久久免费av| 亚洲综合色网址| 亚洲精品第二区| 色94色欧美一区二区| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网| 丝袜美腿诱惑在线| av天堂久久9| 久久精品国产亚洲av天美| 亚洲一码二码三码区别大吗| 精品少妇久久久久久888优播| 亚洲伊人色综图| 男的添女的下面高潮视频| 少妇猛男粗大的猛烈进出视频| 免费在线观看完整版高清| 久久久精品94久久精品| 国产精品国产av在线观看| 熟女少妇亚洲综合色aaa.| 涩涩av久久男人的天堂| 日日撸夜夜添| 伊人久久国产一区二区| 一区二区av电影网| 久久久久久久大尺度免费视频| 在线观看免费日韩欧美大片| 午夜激情av网站| 男女免费视频国产| 欧美国产精品va在线观看不卡| 国产成人精品在线电影| 久久精品国产综合久久久| 亚洲经典国产精华液单| 久久国产亚洲av麻豆专区| 国产精品二区激情视频| 我的亚洲天堂| 久久鲁丝午夜福利片| 精品第一国产精品| 午夜福利,免费看| 男女高潮啪啪啪动态图| 精品视频人人做人人爽| 多毛熟女@视频| 国产免费又黄又爽又色| 亚洲av福利一区| 一级爰片在线观看| 亚洲综合色网址| 欧美在线黄色| 国产精品三级大全| 日韩一区二区视频免费看| 高清av免费在线| 男人添女人高潮全过程视频| 午夜久久久在线观看| 巨乳人妻的诱惑在线观看| 久久久久视频综合| 亚洲三级黄色毛片| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| 黑人猛操日本美女一级片| 久久精品熟女亚洲av麻豆精品| 欧美精品国产亚洲| 精品视频人人做人人爽| 激情五月婷婷亚洲| 亚洲国产精品999| 丰满迷人的少妇在线观看| 欧美黄色片欧美黄色片| 精品少妇久久久久久888优播| 深夜精品福利| 人人澡人人妻人| 亚洲久久久国产精品| 我的亚洲天堂| av有码第一页| 看免费av毛片| 色94色欧美一区二区| 大片免费播放器 马上看| 国产白丝娇喘喷水9色精品| 男男h啪啪无遮挡| 精品酒店卫生间| 国产免费视频播放在线视频| 免费看不卡的av| 色吧在线观看| 久热久热在线精品观看| 国产精品成人在线| 热re99久久精品国产66热6| 午夜av观看不卡| 飞空精品影院首页| 亚洲国产av影院在线观看| 亚洲国产精品国产精品| 成年女人毛片免费观看观看9 | 精品久久久精品久久久| 大话2 男鬼变身卡| 午夜激情av网站| 日韩av在线免费看完整版不卡| 水蜜桃什么品种好| 丝袜人妻中文字幕| 国产精品久久久久久久久免| 久久这里只有精品19| 久久久久久久国产电影| 亚洲第一区二区三区不卡| 亚洲国产精品一区二区三区在线| 制服人妻中文乱码| 国产xxxxx性猛交| 日韩中字成人| 少妇熟女欧美另类| 高清黄色对白视频在线免费看| 久久精品国产综合久久久| 18禁国产床啪视频网站| 中文字幕色久视频| 天天躁夜夜躁狠狠躁躁| 母亲3免费完整高清在线观看 | 国产黄色视频一区二区在线观看| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 亚洲欧美精品自产自拍| 一级毛片我不卡| 国产精品一区二区在线观看99| 国产免费视频播放在线视频| 国产欧美亚洲国产| 人人妻人人添人人爽欧美一区卜| 亚洲男人天堂网一区| 两个人看的免费小视频| 国产精品亚洲av一区麻豆 | 亚洲伊人久久精品综合| 飞空精品影院首页| 国产精品一国产av| 在线观看www视频免费| 久久99精品国语久久久| 亚洲四区av| 国产精品av久久久久免费| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀 | 黄片小视频在线播放| 精品久久久久久电影网| 亚洲精品成人av观看孕妇| 熟女电影av网| 国产视频首页在线观看| 国产一区有黄有色的免费视频| 亚洲av成人精品一二三区| 男人爽女人下面视频在线观看| 国产白丝娇喘喷水9色精品| 人人妻人人添人人爽欧美一区卜| 熟女av电影| 青春草亚洲视频在线观看| 亚洲视频免费观看视频| 国产一区亚洲一区在线观看| 久久人人97超碰香蕉20202| 制服诱惑二区| 激情视频va一区二区三区| 天天操日日干夜夜撸| 一区二区三区四区激情视频| 久久婷婷青草| 成人黄色视频免费在线看| 91精品国产国语对白视频| 国产一区有黄有色的免费视频| 免费女性裸体啪啪无遮挡网站| 精品久久久精品久久久| 国产白丝娇喘喷水9色精品| 久久青草综合色| 欧美97在线视频| 看非洲黑人一级黄片| 免费黄色在线免费观看| 满18在线观看网站| 青青草视频在线视频观看| 1024香蕉在线观看| 男女边吃奶边做爰视频| 精品视频人人做人人爽| 99久久中文字幕三级久久日本| 国产精品蜜桃在线观看| 国产精品成人在线| 国产精品免费视频内射| 欧美精品一区二区大全| 欧美 亚洲 国产 日韩一| 美女脱内裤让男人舔精品视频| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 国产日韩一区二区三区精品不卡| 午夜免费观看性视频| 一边亲一边摸免费视频| 日本欧美视频一区| 国产精品亚洲av一区麻豆 | 精品亚洲成a人片在线观看| 99香蕉大伊视频| 如何舔出高潮| 欧美日韩一级在线毛片| 国产成人欧美| 香蕉国产在线看| 国产成人欧美| 欧美人与性动交α欧美精品济南到 | 亚洲成人一二三区av| 免费观看在线日韩| 卡戴珊不雅视频在线播放| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 亚洲精品一二三| 国产激情久久老熟女| 久久精品国产亚洲av天美| 免费高清在线观看视频在线观看| 亚洲欧美清纯卡通| 精品人妻熟女毛片av久久网站| 亚洲欧洲日产国产| av线在线观看网站| 丝袜人妻中文字幕| 侵犯人妻中文字幕一二三四区| 视频在线观看一区二区三区| 丝袜脚勾引网站| 中文天堂在线官网| 一级黄片播放器| 一区二区三区精品91| 另类精品久久| 日韩精品有码人妻一区| 极品少妇高潮喷水抽搐| 亚洲综合色惰| 老女人水多毛片| 大片电影免费在线观看免费| 久热这里只有精品99| 免费av中文字幕在线| 亚洲国产毛片av蜜桃av| 2018国产大陆天天弄谢| 国产乱人偷精品视频| 在线观看免费日韩欧美大片| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 亚洲av免费高清在线观看| 伊人亚洲综合成人网| 少妇 在线观看| 亚洲一码二码三码区别大吗| 亚洲四区av| 香蕉丝袜av| 9191精品国产免费久久| 国产精品av久久久久免费| 在现免费观看毛片| 日韩一本色道免费dvd| 麻豆乱淫一区二区| freevideosex欧美| 国产精品久久久久久精品古装| 国产国语露脸激情在线看| 老女人水多毛片| 国产亚洲午夜精品一区二区久久| 国产午夜精品一二区理论片| 亚洲欧洲精品一区二区精品久久久 | 精品久久蜜臀av无| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 国产在线免费精品| 欧美日韩国产mv在线观看视频| av.在线天堂| 亚洲欧美日韩另类电影网站| av免费观看日本| 国产乱人偷精品视频| 欧美精品国产亚洲| 久久毛片免费看一区二区三区| 午夜激情久久久久久久| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频| 美女中出高潮动态图| 国产激情久久老熟女| 日本91视频免费播放| 制服诱惑二区| 久久久久久免费高清国产稀缺| 一区二区三区精品91| 在线亚洲精品国产二区图片欧美| 日韩制服骚丝袜av| 自线自在国产av| 亚洲一级一片aⅴ在线观看| 美女午夜性视频免费| 国产日韩一区二区三区精品不卡| 久久精品国产亚洲av涩爱| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 超碰成人久久| 免费观看a级毛片全部| www.自偷自拍.com| 国产一区二区在线观看av| 国产野战对白在线观看| 久久久久国产一级毛片高清牌| 国产av一区二区精品久久| 制服人妻中文乱码| 伦精品一区二区三区| 日本av手机在线免费观看| 欧美激情极品国产一区二区三区| 久久 成人 亚洲| 国产高清不卡午夜福利| 亚洲综合精品二区| 天堂中文最新版在线下载| 国产精品99久久99久久久不卡 | 久久精品国产亚洲av高清一级| 欧美国产精品va在线观看不卡| 精品福利永久在线观看| 久久人人爽人人片av| 日韩中字成人| 女人高潮潮喷娇喘18禁视频| 国产一区二区在线观看av| 国产成人欧美| 国产 一区精品| 亚洲av在线观看美女高潮| 看非洲黑人一级黄片| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 如何舔出高潮| 亚洲国产av影院在线观看| 亚洲三级黄色毛片| 国产亚洲最大av| 99国产精品免费福利视频| 在线观看免费高清a一片| 国产精品av久久久久免费| 国产精品久久久久久精品电影小说| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 亚洲天堂av无毛| 亚洲少妇的诱惑av|