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

    應(yīng)用物種敏感性分布評(píng)估三氯卡班對(duì)我國(guó)淡水環(huán)境的生態(tài)風(fēng)險(xiǎn)

    2017-10-24 10:23:29楊先海范德玲劉濟(jì)寧石利利
    關(guān)鍵詞:三氯淡水毒性

    汪 貞, 楊先海, 范德玲, 郭 敏, 劉濟(jì)寧, 石利利

    (環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所, 江蘇 南京 210042)

    應(yīng)用物種敏感性分布評(píng)估三氯卡班對(duì)我國(guó)淡水環(huán)境的生態(tài)風(fēng)險(xiǎn)

    汪 貞, 楊先海, 范德玲, 郭 敏, 劉濟(jì)寧①, 石利利

    (環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所, 江蘇 南京 210042)

    應(yīng)用物種敏感性分布(species sensitivity distribution,SSD)方法構(gòu)建了三氯卡班(triclocarban,TCC)對(duì)淡水生物的SSD曲線,計(jì)算了TCC對(duì)淡水生物的5%急性危害濃度(HC5)。采用對(duì)數(shù)正態(tài)(log-normal)分布模型得到的急性ρ(HC5)=3.85 μg·L-1,在評(píng)估因子(AF)取值為3、急慢性毒性比(ACR)為39.3的基礎(chǔ)上,計(jì)算了TCC的慢性預(yù)測(cè)無(wú)效應(yīng)濃度(PNEC)為32.7 ng·L-1。采用商值法對(duì)21個(gè)有數(shù)據(jù)報(bào)道的典型地表水調(diào)查點(diǎn)水體中TCC的生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)估,結(jié)果表明,TCC高風(fēng)險(xiǎn)比例為28.6%, 中風(fēng)險(xiǎn)比例為47.6%,低風(fēng)險(xiǎn)比例為23.8%。TCC對(duì)我國(guó)地表水的影響應(yīng)該引起關(guān)注。

    三氯卡班; 風(fēng)險(xiǎn)商; 物種敏感性分布; 地表水

    三氯卡班(triclocarban,TCC)是一種廣譜抗菌劑,能殺滅革蘭陽(yáng)性菌、革蘭陰性菌、真菌、酵母菌和病毒等微生物,被廣泛地應(yīng)用于洗滌劑、牙膏、肥皂、空氣清新劑、洗手液等個(gè)人護(hù)理用品以及化妝品、織物消毒劑和醫(yī)用消毒液等生產(chǎn)領(lǐng)域。相比于毒性已經(jīng)被廣泛研究認(rèn)可的三氯生(triclosan,TCS),TCC一直被認(rèn)為是無(wú)毒安全的殺菌劑,是TCS的良好替代品,自1957年問世以來(lái),被世界各國(guó)大量使用,其環(huán)境危害研究也一直沒有得到足夠重視。自2004年環(huán)境中TCC的ng·L-1級(jí)檢測(cè)方法逐漸建立與完善以來(lái)[1],越來(lái)越多的調(diào)查發(fā)現(xiàn)TCC已存在于各種環(huán)境介質(zhì)中,在許多湖泊、河流以及生物體內(nèi)都有檢出[2-4]。盡管水環(huán)境中TCC濃度相對(duì)較低,一般檢出濃度在ng·L-1~μg·L-1水平[5],但是其具有化學(xué)性質(zhì)穩(wěn)定、難降解,親脂性高,容易在生物體內(nèi)富集等特點(diǎn),可以隨著水生生物食物鏈生物放大,最終通過食物進(jìn)入人體中[6]。已有大量證據(jù)表明環(huán)境中TCC能對(duì)生態(tài)系統(tǒng)和人類健康產(chǎn)生危害[7]。TCC還有可能引發(fā)包括癌癥、生殖功能障礙和發(fā)育異常等病癥的諸多問題,是一種新型的內(nèi)分泌干擾物(EDCs)[8]。美國(guó)食品藥品監(jiān)督管理局(Food and Drug Administration,FDA)于2016年9月2日發(fā)布一項(xiàng)最終裁決令,將禁止含有TCC等19種活性成分的非處方(OTC)消費(fèi)品抗菌皂繼續(xù)上市銷售[9]。此外,TCC曾2次被提名列入美國(guó)環(huán)境保護(hù)局(United States Environmental Protection Agency,USEPA)飲用水污染物候選清單(CCL-3,CCL-4)。因此,TCC對(duì)生態(tài)系統(tǒng)和人體健康造成的潛在影響在逐漸受到關(guān)注,開展人體健康風(fēng)險(xiǎn)和水生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)等研究非常必要。

    物種敏感性分布法(species sensitivity distribution,SSD)是基于不同物種對(duì)某一環(huán)境脅迫的敏感度服從一定累積概率分布假設(shè),以統(tǒng)計(jì)分布模型來(lái)描述不同物種樣本對(duì)脅迫因素的敏感性差異,實(shí)現(xiàn)將單一物種的測(cè)試結(jié)果外推至生態(tài)系統(tǒng)的風(fēng)險(xiǎn)評(píng)估方法[10]。該方法已列入歐盟風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)導(dǎo)則[11],廣泛應(yīng)用于生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)和環(huán)境質(zhì)量基準(zhǔn)(environmental quality criteria,EQC)的制定。筆者通過構(gòu)建TCC的急性物種敏感性分布曲線,參考我國(guó)淡水環(huán)境中TCC濃度現(xiàn)狀,評(píng)估我國(guó)淡水環(huán)境中TCC生態(tài)風(fēng)險(xiǎn),以期為我國(guó)水體TCC的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)與管理提供參考。

    1 材料與方法

    1.1 方法原理與基本步驟

    物種敏感性分布法理論認(rèn)為:由于生活史、生理構(gòu)造、行為特征和地理分布等不同,不同門類的生物在毒理學(xué)上表現(xiàn)為不同物種對(duì)污染物有不同的劑量-效應(yīng)響應(yīng)關(guān)系,即在結(jié)構(gòu)復(fù)雜的生態(tài)系統(tǒng)中,不同物種對(duì)某一脅迫因素(如有毒化學(xué)品)的敏感程度服從一定的(累積)概率分布。SSD的用法包括正向(forward use)和反向(inverse use)2種。正向方法主要用于生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),即已知污染物濃度水平,通過SSD曲線計(jì)算潛在影響比例(potential affected fraction,PAF),用以表征生態(tài)系統(tǒng)或者不同類別生物的生態(tài)風(fēng)險(xiǎn);而反向方法通過確定保護(hù)一個(gè)生態(tài)系統(tǒng)中大部分物種的污染物濃度水平來(lái)制定環(huán)境質(zhì)量基準(zhǔn),一般使用5%危害濃度(hazardous concentration 5%,HC5),HC5指影響不超過5%的物種即可以保護(hù)95%以上的物種時(shí)對(duì)應(yīng)的急性濃度/慢性濃度[12]。筆者將通過反向方法推算HC5,在其基礎(chǔ)上推算預(yù)測(cè)無(wú)效應(yīng)濃度(predicted no effect concentration,PNEC),最后結(jié)合環(huán)境濃度以商值法表征我國(guó)淡水環(huán)境中TCC的生態(tài)風(fēng)險(xiǎn)。

    物種敏感度分布曲線的建立步驟:(1)毒性數(shù)據(jù)獲取;(2)數(shù)據(jù)處理;(3)曲線擬合;(4)HC5的計(jì)算。風(fēng)險(xiǎn)商的計(jì)算步驟:(1)環(huán)境濃度數(shù)據(jù)整理;(2)PNEC預(yù)測(cè);(3)風(fēng)險(xiǎn)商(risk quotient,RQ)的計(jì)算。

    1.2 毒理數(shù)據(jù)的獲取與整理

    SSD的構(gòu)建可以使用半數(shù)致死濃度(LC50)、半數(shù)效應(yīng)濃度(EC50)、無(wú)可見效應(yīng)濃度(NOEC)及最低可見效應(yīng)濃度(LOEC)等急性或慢性數(shù)據(jù),由于TCC相關(guān)的慢性毒性數(shù)據(jù)較少,故使用急性毒性數(shù)據(jù)構(gòu)建SSD。利用USEPA的ECOTOX數(shù)據(jù)庫(kù)(http:∥www.epa.gov/ecotox/)、高產(chǎn)量物質(zhì)信息系統(tǒng)(HPVIS,https:∥ofmpub.epa.gov/oppthpv/hpv-ez.html-menu)以及相關(guān)文獻(xiàn),參考表1的毒性數(shù)據(jù)篩選標(biāo)準(zhǔn),搜集TCC對(duì)水生生物的急性、慢性毒理數(shù)據(jù)。對(duì)于同一物種有多個(gè)數(shù)據(jù)的情況,采用其所有濃度數(shù)據(jù)的幾何平均值。共收集到4門14科15個(gè)淡水物種的急性毒性數(shù)據(jù)和2門4科4個(gè)淡水物種的慢性毒性數(shù)據(jù)[5,13-14](表2~3)。

    表1毒性數(shù)據(jù)篩選標(biāo)準(zhǔn)

    Table1Criteriaforselectingtoxicitydata

    數(shù)據(jù)類型暴露終點(diǎn)暴露時(shí)間/d試驗(yàn)介質(zhì)急性毒性LC50/EC50≤10淡水慢性毒性NOEC/LOEC≥20淡水

    LC50為半數(shù)致死濃度; EC50為半數(shù)效應(yīng)濃度; NOEC為無(wú)可見效應(yīng)濃度; LOEC為最低可見效應(yīng)濃度。

    1.3 數(shù)據(jù)處理

    將表2的毒性數(shù)據(jù)進(jìn)行整理,按照數(shù)值大小對(duì)物種毒性數(shù)據(jù)進(jìn)行排序,每個(gè)物種的累積概率(P)計(jì)算公式為P=i/(n+1)。其中,i為物種從小到大排序的秩;n為物種數(shù)。

    1.4 SSD曲線擬合

    目前,還沒有理論研究證明SSD屬于某一特定曲線形式,因此可選擇不同的擬合方法。常用的擬合方法包括log-normal、log-logistic、Gaussian和Burr Ⅲ 分布模型等,美國(guó)和歐盟推薦使用log-normal分布模型擬合SSD曲線,而澳大利亞和新西蘭則推薦使用Burr Ⅲ分布[15]??梢杂糜赟SD曲線擬合的軟件工具有SigmaPlot、ETX 2.0和BurrliOZ等。筆者采用SigmaPlot 13.0進(jìn)行sigmoid、Gaussian、log-normal和log-logistic分布擬合。

    1.5 HC5與PNEC的計(jì)算

    在SSD擬合曲線上對(duì)應(yīng)5%累積概率的污染物質(zhì)量濃度為HC5。根據(jù)擬合方程,利用Matlab軟件的solve函數(shù)可求出5%累積概率對(duì)應(yīng)的HC5。

    在慢性毒性數(shù)據(jù)缺乏的情況下,PNEC(CPNE)的計(jì)算方法如下:采用急性慢性毒性比(ACR,RAC),即3個(gè)或3個(gè)以上物種的急性慢性毒性比的幾何平均值(要求至少包含魚類、無(wú)脊椎動(dòng)物類和敏感淡水類生物中的1種)[16],將急性HC5轉(zhuǎn)換為慢性HC5,再參照歐盟《化學(xué)品的注冊(cè)、評(píng)估、授權(quán)和限制》(REACH)法規(guī)對(duì)危害評(píng)估因子的選擇方法[11],將評(píng)估因子(AF,FA)取值為3,即CPNE=ρ(HC5)/RAC/FA。根據(jù)表2~3收集的急性、慢性毒性數(shù)據(jù),推算得到的RAC為39.3。

    表2三氯卡班的淡水急性毒性數(shù)據(jù)

    Table2Acutetoxicitydataoftriclocarbaninfreshwater

    門科 物種終點(diǎn)暴露時(shí)間/h終點(diǎn)值/(μg·L-1)數(shù)據(jù)來(lái)源綠藻門小球藻科羊角月牙藻(Pseudokirchneriellasubcapitata)EC507217.0A柵藻科近具刺?hào)旁?Scenedesmussubspicatus)EC507220.0B節(jié)肢動(dòng)物門溞科模糊網(wǎng)紋溞(Ceriodaphniadubia)EC50483.10A溞科大型溞(Daphniamagna)EC504810.2C鉤蝦科鉤蝦(Gammarusfasciatus)EC507213.0A糠蝦科巴西尼糠蝦(Mysidopsisbahia)LC50/EC5024~9618.0C搖蚊科擬長(zhǎng)跗搖蚊(Paratanytarsusparthenogeneticus)EC504880.0A原生動(dòng)物門四膜科嗜熱四膜蟲(Tetrahymenathermophila)EC5024295文獻(xiàn)[5]脊索動(dòng)物門蛙科黑斑蛙(Pelophylaxnigromaculatus)LC5096217文獻(xiàn)[13]負(fù)子蟾科非洲爪蟾(Xenopuslaevis)LC5096252文獻(xiàn)[13]太陽(yáng)魚科藍(lán)鰓太陽(yáng)魚(Lepomismacrochirus)LC5048~9670.3C青鳉科青鳉(Oryziaslatipes)LC509685.0文獻(xiàn)[14]鯉科黑頭呆魚(Pimephalespromelas)LC509692.0A鮭科虹鱒(Oncorhynchusmykiss)LC5096133C科斑點(diǎn)叉尾(Ictaluruspunctatus)LC5048~96140A

    LC50為半數(shù)致死濃度; EC50為半數(shù)效應(yīng)濃度; A為ECOTOX數(shù)據(jù)庫(kù); B為高產(chǎn)量物質(zhì)信息系統(tǒng)(HPVIS); C為ECOTOX數(shù)據(jù)庫(kù)& HPVIS。

    表3三氯卡班的慢性毒性數(shù)據(jù)

    Table3Chronictoxicitydataoftriclocarbaninfreshwater

    門物種終點(diǎn)暴露時(shí)間/d終點(diǎn)值/(μg·L-1)范圍幾何均值數(shù)據(jù)來(lái)源節(jié)肢動(dòng)物門大型溞(Daphniamagna)NOEC/LOEC21~420.25~150.74C脊索動(dòng)物門黑頭呆魚(Pimephalespromelas)NOEC21~351.6~9.83.00C節(jié)肢動(dòng)物門巴西尼糠蝦(Mysidopsisbahia)NOEC/LOEC280.06~0.60.20B節(jié)肢動(dòng)物門擬長(zhǎng)跗搖蚊(Paratanytarsusparthenogeneticus)NOEC/LOEC21~281.28~13.624.23A

    NOEC為無(wú)可見效應(yīng)濃度; LOEC為最低可見效應(yīng)濃度; A為ECOTOX數(shù)據(jù)庫(kù); B為高產(chǎn)量物質(zhì)信息系統(tǒng)(HPVIS); C為ECOTOX數(shù)據(jù)庫(kù)& HPVIS。

    1.6 生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)是通過一個(gè)濃度閾值或者風(fēng)險(xiǎn)值來(lái)量化污染物對(duì)環(huán)境的生態(tài)危害,是通過風(fēng)險(xiǎn)表征來(lái)量化表達(dá)的。商值法、概率法(probabilistic ecological risk assessment,PERA)和安全閾值法(the margin of safety,MOS10)是常用的量化表征方法[17]。筆者采用商值法評(píng)價(jià)TCC的生態(tài)風(fēng)險(xiǎn)。商值法中RQ(QR)的計(jì)算公式為QR=CPE/CPNE,其中,CPE為預(yù)測(cè)環(huán)境濃度(PEC)。筆者研究中,用TCC的環(huán)境濃度(measured environmental concentration,MEC)代替 PEC 進(jìn)行生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。TCC的環(huán)境濃度數(shù)據(jù)主要來(lái)源于文獻(xiàn)數(shù)據(jù)和筆者所在研究團(tuán)隊(duì)的實(shí)測(cè)數(shù)據(jù),包括池塘、湖泊、河流和水庫(kù)(飲用水源地)等。

    當(dāng)QR<0.1時(shí),表示化學(xué)品對(duì)環(huán)境存在的風(fēng)險(xiǎn)較低(低風(fēng)險(xiǎn)),當(dāng)QR為0.1~1.0時(shí),表明化學(xué)品對(duì)環(huán)境存在一定風(fēng)險(xiǎn)(中風(fēng)險(xiǎn)),需要對(duì)相關(guān)風(fēng)險(xiǎn)源展開跟蹤觀察;當(dāng)QR>1.0時(shí),表明化學(xué)品對(duì)環(huán)境存在比較嚴(yán)重的風(fēng)險(xiǎn)(高風(fēng)險(xiǎn)),需要采取相應(yīng)的風(fēng)險(xiǎn)削減措施[18-19]。

    2 結(jié)果與分析

    2.1 我國(guó)不同水體中TCC的環(huán)境濃度

    根據(jù)文獻(xiàn)[2-3,20-24]以及筆者所在研究團(tuán)隊(duì)未發(fā)表的數(shù)據(jù)整理了我國(guó)不同地表水水體中TCC的環(huán)境濃度(表4)。由表4可知,我國(guó)不同類型的地表水中TCC的環(huán)境濃度差異明顯。池塘中ρ(TCC)范圍為69.9~1 110 ng·L-1,大型湖泊中ρ(TCC)為10.3~17.2 ng·L-1,河流中ρ(TCC)為河流>大型湖泊>飲用水源地。

    表4我國(guó)不同水體中三氯卡班(TCC)的環(huán)境濃度

    Table4ConcentrationsoftriclocarbaninsurfacewaterenvironmentofChina

    地表水類型調(diào)查點(diǎn)ρ(TCC)/(ng·L-1)范圍均值樣點(diǎn)數(shù)檢出率/%數(shù)據(jù)來(lái)源池塘大明湖—69.9 1100[20]白鷺湖—11101100[21]湖泊巢湖—13.33—[22]太湖—10.33—[22]鄱陽(yáng)湖—13.75—[22]洞庭湖—17.25—[22]河流小清河4231100[20]流溪河

    “—”表示無(wú)數(shù)據(jù); LOQ為檢測(cè)下限; ND表示未檢測(cè)。*表示筆者所在研究團(tuán)隊(duì)于2016年對(duì)無(wú)錫、寧波、重慶和深圳4地的飲用水源地進(jìn)行采樣調(diào)查獲得的數(shù)據(jù)。

    2.2 不同曲線方程擬合結(jié)果

    采用SigmaPlot 13.0建立了4個(gè)曲線方程,各方程的參數(shù)、校正決定系數(shù)(Radj2)和HC5見表5,相應(yīng)的擬合曲線見圖1。sigmoid、Gaussian、log-normal和logistic分布的Radj2均>0.95,盡管Gaussian分布的Radj2最大,殘差平方和(RSS)最小,但與其他分布模型差異不大。log-normal分布模型是美國(guó)和歐盟推薦使用的SSD曲線擬合模型,后續(xù)的評(píng)估將采用log-normal分布得到的HC5,即3.85 μg·L-1。

    表5擬合SSD曲線方程和HC5值

    Table5EquationforfittingSSDcurvesandHC5values

    分布函數(shù)函數(shù)參數(shù)abx0校正決定系數(shù)Radj2殘差平方和RSSρ(HC5)/(μg·L-1)sigmoid 2.02730.6827 2.58390.97130.03141.18Gaussian2.12411.56804.48510.97190.03071.54log-normal364492.51781.757913562.32520.96980.03313.85logistic1182.1669-1.8957107.99160.97150.03123.41

    2.3 PNEC預(yù)測(cè)

    采用log-normal分布獲得的ρ(HC5)=3.85 μg·L-1,采用1.5節(jié)中的公式計(jì)算的PNEC為32.7 ng·L-1。

    2.4 基于商值法的長(zhǎng)期生態(tài)風(fēng)險(xiǎn)評(píng)估

    根據(jù)國(guó)內(nèi)文獻(xiàn)以及筆者所在研究團(tuán)隊(duì)的實(shí)測(cè)數(shù)據(jù),我國(guó)水環(huán)境中TCC濃度范圍為

    圖1 不同模型擬合三氯卡班的急性物種敏感度分布曲線Fig.1 SSD curves of triclocarban relative to fitting model

    表6我國(guó)淡水環(huán)境中三氯卡班的生態(tài)風(fēng)險(xiǎn)

    Table6EcologicalriskoftriclocarbaninfreshwaterenvironmentofChina

    地表水類型調(diào)查點(diǎn)風(fēng)險(xiǎn)商風(fēng)險(xiǎn)等級(jí)池塘大明湖 2.14高白鷺湖33.90高湖泊巢湖0.41中太湖0.31中鄱陽(yáng)湖0.42中洞庭湖0.53中河流小清河12.90高流溪河0.23中石井河4.83高九龍河0.43中黃河0.18中遼河0.43中海河0.95中珠江2.11高長(zhǎng)江0.17中東江1.09高飲用水源地駱馬湖0.09低無(wú)錫0低寧波0.03低重慶0.03低深圳0低

    3 討論

    足夠的生態(tài)毒性數(shù)據(jù)和環(huán)境濃度數(shù)據(jù)是采用物種敏感性分布法開展生態(tài)風(fēng)險(xiǎn)評(píng)估的基礎(chǔ),經(jīng)濟(jì)合作與發(fā)展組織(Organization for Economic Co-operation and Development,OECD)以及澳大利亞水質(zhì)標(biāo)準(zhǔn)中推薦的最小數(shù)據(jù)量為5個(gè)[25],其他研究或政府指導(dǎo)文件中推薦的最小數(shù)據(jù)量為8~10個(gè)不等。筆者整理了4門14科15個(gè)物種的急性毒性數(shù)據(jù),符合USEPA構(gòu)建SSD曲線時(shí)3門8科生態(tài)毒性數(shù)據(jù)的要求。4類21個(gè)地表水調(diào)查點(diǎn)基本涵蓋了我國(guó)淡水生態(tài)系統(tǒng)的不同類型,環(huán)境濃度數(shù)據(jù)可以代表我國(guó)淡水生態(tài)系統(tǒng)中TCC的污染現(xiàn)狀。需要指出的是,用于構(gòu)建SSD曲線的15個(gè)物種中,青鳉(Oryziaslatipes)、黑斑蛙(Pelophylaxnigromaculatus)為本土物種,羊角月牙藻(Pseudokirchneriellasubcapitata)、近具刺?hào)旁?Scenedesmussubspicatus)、大型溞(Daphniamagna)、嗜熱四膜蟲(Tetrahymenathermophila)等物種在我國(guó)的水生態(tài)系統(tǒng)中也有分布。因此,筆者的評(píng)價(jià)結(jié)論可以為我國(guó)的淡水生態(tài)

    系統(tǒng)TCC的生態(tài)風(fēng)險(xiǎn)評(píng)估提供參考。

    在慢性毒性數(shù)據(jù)缺失的情況下,ACR法可以較準(zhǔn)確地實(shí)現(xiàn)急性HC5到慢性HC5的外推[18]。在筆者研究中,ACR是關(guān)鍵參數(shù)。不同機(jī)構(gòu)或?qū)W者在ACR的取值上有不同的規(guī)定和看法。USEPA在計(jì)算 ACR時(shí)要求采用至少3個(gè)物種的急慢性毒性比的幾何平均值,包括1種魚類、1種無(wú)脊椎動(dòng)物和1種敏感的淡水生物[16];OECD和澳大利亞推薦使用的ACR默認(rèn)值為10[26];也有一些研究采用的ACR值為100。在使用了該方法的已有研究中,不同的評(píng)估物質(zhì)具有不同的ACR值,不同的研究者使用了不同的ACR值。如芐氯菊酯和高效氰戊菊酯的ACR為11.4[27],毒死蜱、硫丹、氟樂靈和百菌清4種農(nóng)藥的ACR為7.42~34.6[28],五氯酚、2,4-二氯酚、2,4,6-三氯酚、4-硝基酚、2,4-二硝基酚和2-甲基-4,6-二硝基酚的ACR為1.65~7.58[18]。而在另外一項(xiàng)研究中,五氯酚、2-氯苯酚、2,4-二氯苯酚和2,4,6-三氯苯酚的ACR為2.77~20.49[29]。筆者推導(dǎo)的ACR為39.3,盡管數(shù)值上大于其他物質(zhì)的ACR,但數(shù)量級(jí)水平基本一致。在此基礎(chǔ)上,筆者研究得到的PNEC為32.7 ng·L-1,文獻(xiàn)[3]基于大型溞的慢性毒性(21 d NOEC為2.9 μg·L-1)以50為評(píng)估因子得到的TCC的PNEC為58 ng·L-1,兩者基本在相同的數(shù)量級(jí)水平,但筆者研究中的PNEC更加保守。從風(fēng)險(xiǎn)評(píng)估結(jié)果來(lái)看,不同類型地表水有風(fēng)險(xiǎn)級(jí)別的差異,如池塘均為高風(fēng)險(xiǎn),4大淡水湖泊(巢湖、太湖、鄱陽(yáng)湖、洞庭湖)均為中等風(fēng)險(xiǎn),有數(shù)據(jù)報(bào)道的10條河流中,4條為高風(fēng)險(xiǎn),6條為中等風(fēng)險(xiǎn),5個(gè)飲用水源地調(diào)查點(diǎn)均為低風(fēng)險(xiǎn)。由于環(huán)境中TCC主要來(lái)源于生活污水,而居民區(qū)附近的池塘、河流等水體往往成為生活污水直接排放的受納水體,小清河、石井河、珠江和東江等河流的TCC濃度數(shù)據(jù)來(lái)自濟(jì)南、廣州等大城市的市內(nèi)河段,這可能是池塘以及小清河、石井河、珠江、東江水環(huán)境中TCC均為高風(fēng)險(xiǎn)的原因。一項(xiàng)針對(duì)藻、溞、魚開展的TCC風(fēng)險(xiǎn)評(píng)估研究表明城市河流中還沒有污水處理設(shè)施覆蓋的集水區(qū),TCC的風(fēng)險(xiǎn)值大于1[14]。因此,筆者研究結(jié)果基本上與已有研究結(jié)果一致。

    值得注意的是,筆者研究中搜索的數(shù)據(jù)量(包括環(huán)境暴露水平和生物毒性數(shù)據(jù))有限,由于慢性毒性數(shù)據(jù)較少,慢性HC5是基于4個(gè)物種推導(dǎo)的ACR進(jìn)行轉(zhuǎn)換得到的,ACR取值不同可能會(huì)使推導(dǎo)的慢性HC5存在一定差異,以及評(píng)估因子AF的取值,以上過程都存在一定的不確定性。

    綜上所述,TCC對(duì)我國(guó)淡水環(huán)境的影響應(yīng)該引起環(huán)境管理部門的重視。鑒于筆者研究所得到的結(jié)論可能受本土物種毒性數(shù)據(jù)少、以急性毒性數(shù)據(jù)推導(dǎo)長(zhǎng)期生態(tài)風(fēng)險(xiǎn)、評(píng)估因子等不確定性因素的影響,建議深入研究TCC對(duì)我國(guó)本土物種的慢性毒性影響,系統(tǒng)調(diào)查TCC在我國(guó)水環(huán)境中的賦存現(xiàn)狀,以進(jìn)一步研究針對(duì)保護(hù)我國(guó)水生態(tài)系統(tǒng)與人體健康的TCC標(biāo)準(zhǔn)限值。

    [1] HALDEN R U,PAULL D H.Analysis of Triclocarban in Aquatic Samples by Liquid Chromatography Electrospray Ionization Mass Spectrometry[J].Environmental Science and Technology,2004,38(18):4849-4855.

    [2] ZHAO J L,ZHANG Q Q,CHEN F,etal.Evaluation of Triclosan and Triclocarban at River Basin Scale Using Monitoring and Modeling Tools:Implications for Controlling of Urban Domestic Sewage Discharge[J].Water Research,2013,47(1):395-405.

    [3] ZHAO J L,YING G G,LIU Y S,etal.Occurrence and Risks of Triclosan and Triclocarban in the Pearl River System,South China:From Source to the Receiving Environment[J].Journal of Hazardous Materials,2010,179(1/2/3):215-222.

    [4] GAUTAM P,CARSELLA J S,KINNEY C A.Presence and Transport of the Antimicrobials Triclocarban and Triclosan in a Wastewater-Dominated Stream and Freshwater Environment[J].Water Research,2014,48:247-256.

    [5] 白琦鋒,高禮,袁濤.低劑量三氯生和三氯卡班對(duì)嗜熱四膜蟲的毒性研究[J].環(huán)境化學(xué),2012(5):720-725.[BAI Qi-feng,GAO Li,YUAN Tao.Toxicity of Low Concentration Exposures of Triclosan and Triclocarban onTetrahymenathermophila[J].Environmental Chemistry,2012,31(5):720-725.]

    [6] WU C,SPONGBERG A L,WITTER J D.Adsorption and Degradation of Triclosan and Triclocarban in Soils and Biosolids-Amended Soils[J].Journal of Agricultural and Food Chemistry,2009,57(11):4900-4905.

    [7] 李根.環(huán)境激素三氯卡班與雌二醇對(duì)大型溞的復(fù)合干擾效應(yīng)[D].廣州:暨南大學(xué),2014.[LI Gen.The Joint Effects of Triclocarban and Estradiol onDaphniamagna[D].Guangzhou:Jinan University,2014.]

    [8] CHEN J G,AHN K C,GEE N A,etal.Triclocarban Enhances Testosterone Action:A New Type of Endocrine Disruptor?[J].Endocrinology,2008,149(3):1173-1179.

    [9] FDA.Safety and Effectiveness of Consumer Antiseptics; Topical Antimicrobial Drug Products for Over-the-Counter Human Use[J].Federal Register,2016,81(172):61106-61135.

    [10] WHEELER J R,GRIST E P M,LEUNG K M Y,etal.Species Sensitivity Distributions:Data and Model Choice[J].Marine Pollution Bulletin,2002,45:192-202.

    [11] European Chemicals Bureau.Technical Guidance Document on Risk Assessment Part Ⅱ[R].Luxembourg:Office for Official Publications of the European Communities,2003:101-118.

    [12] 胡習(xí)邦,王俊能,許振成,等.應(yīng)用物種敏感性分布評(píng)估DEHP對(duì)區(qū)域水生生態(tài)風(fēng)險(xiǎn)[J].生態(tài)環(huán)境學(xué)報(bào),2012,21(6):1082-1087.[HU Xi-bang,WANG Jun-neng,XU Zhen-cheng,etal.Assessing Aquatic Ecological Risk of DEHP by Species Sensitivity Distributions[J].Ecology and Environmental Sciences,2012,21(6):1082-1087.]

    [13] 高昆,陳娟,李圓圓,等.三氯生和三氯卡班對(duì)兩棲動(dòng)物蝌蚪的急性毒性[J].生態(tài)毒理學(xué)報(bào),2016,11(4):226-231.[GAO Kun,CHEN Juan,LI Yuan-yuan,etal.Acute Toxicity of Triclosan and Triclocarban toAmphibiantadpoles[J].Asian Journal of Ecotoxicology,2016,11(4):226-231.]

    [14] TAMURA I,KAGOTA K,YASUDA Y,etal.Ecotoxicity and Screening Level Ecotoxicological Risk Assessment of Five Antimicrobial Agents:Triclosan,Triclocarban,Resorcinol,Phenoxyethanol and P-Thymol[J].Journal of Applied Toxicology,2012,33(11):1222-1229.

    [15] 劉亞莉,謝玉為,張效偉,等.應(yīng)用物種敏感性分布評(píng)價(jià)敵敵畏對(duì)淡水生物的生態(tài)風(fēng)險(xiǎn)[J].生態(tài)毒理學(xué)報(bào),2016,11(2):531-538.[LIU Ya-li,XIE Yu-wei,ZHANG Xiao-wei,etal.Assessing Ecological Risks of Dichlorvos to Freshwater Organisms by Species Sensitivity Distribution[J].Asian Journal of Ecotoxicology,2016,11(2):531-538.]

    [16] 張怡婷,王蕾,劉濟(jì)寧,等.應(yīng)用不同毒理學(xué)終點(diǎn)評(píng)估酚類物質(zhì)的生態(tài)危害和風(fēng)險(xiǎn)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2016,32(2):326-331.[ZHANG Yi-ting,WANG Lei,LIU Ji-ning,etal.Using Different Toxicological End Points to Evaluate Phenolic Compounds for Ecological Hazard and Risk[J].Journal of Ecology and Rural Environment,2016,32(2):326-331.]

    [17] 雷炳莉,黃圣彪,王子健.生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)理論和方法[J].化學(xué)進(jìn)展,2009,21(增刊1):350-358.[LEI Bing-li,HUANG Sheng-biao,WANG Zi-jian.Theories and Methods of Ecological Risk Assessment[J].Progress in Chemistry,2009,21(Suppl.1):350-358.]

    [18] 高培.壬基酚的水質(zhì)基準(zhǔn)探討和生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[D].青島:中國(guó)海洋大學(xué),2014.[GAO Pei.Derivation of Water Quality Criteria for Nonylphenol and Its Application in Ecological Risk Assessment[D].Qingdao:Ocean University of China,2014.]

    [19] GAO P,LI Z Y,GIBSON M,etal.Ecological Risk Assessment of Nonylphenol in Coastal Waters of China Based on Species Sensitivity Distribution Model[J].Chemosphere,2014,104:113-119.

    [20] 孫靜,趙汝松,王霞,等.液液萃取-HPLC-ESI-MS法同步測(cè)定環(huán)境水樣中的三氯卡班和三氯生[J].山東輕工業(yè)學(xué)院學(xué)報(bào),2011,25(1):35-37.[SUN Jing,ZHAO Ru-song,WANG Xia,etal.Simultaneous Determination of Triclosan and Triclocarban in Environmental Water Samples by Liquid Extration Combined With HPLC-ESI-MS[J].Journal of Shandong Polytechnic University,2011,25(1):35-37.]

    [21] 朱開杭,曾慶玲,沈春花,等.SPE-HPLC法檢測(cè)環(huán)境水樣中的痕量三氯卡班[J].環(huán)境工程學(xué)報(bào),2016,10(1):503-508.[ZHU Kai-hang,ZENG Qing-ling,SHEN Chun-hua,etal.Determination of Trace Triclocarban in Environmental Water Samples by SPE-HPLC[J].Chinese Journal of Environmental Engineering,2016,10(1):503-508.]

    [22] WU C X,HUANG X L,WITTER J D,etal.Occurrence of Pharmaceuticals and Personal Care Products and Associated Environmental Risks in the Central and Lower Yangtze River,China[J].Ecotoxicology and Environmental Safety,2014,106:19-26.

    [23] Lü M,SUN Q,XU H L,etal.Occurrence and Fate of Triclosan and Triclocarban in a Subtropical River and Its Estuary[J].Marine Pollution Bulletin,2014,88(1/2):383-388.

    [24] 張芹,張圣虎,汪貞,等.駱馬湖表層水體中32種PPCPs類物質(zhì)的污染水平、分布特征及風(fēng)險(xiǎn)評(píng)估[J].環(huán)境科學(xué),2017,38(1):162-169.[ZHANG Qin,ZHANG Sheng-hu,WANG Zhen,etal.Pollution Level,Distribution Characteristics and Risk Assessment of 32 PPCPs in Surface Water of Luomahu Lake[J].Environmental Science,2017,38(1):162-169.]

    [25] HOSE G C,PJ V D B.Confirming the Species-Sensitivity Distribution Concept for Endosulfan Using Laboratory,Mesocosm,and Field Data[J].Archives of Environmental Contamination and Toxicology,2004,47(4):511-520.

    [26] ST M J,WARNE M.Derivation of the Australian and New Zealand Water Quality Guidelines for Toxicants[J].Australasian Journal of Ecotoxicology,2001,7:123-136.

    [27] 趙晨,彭書傳,陳天虎,等.巢湖流域和太湖流域沉積物中芐氯菊酯和高效氰戊菊酯的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].環(huán)境科學(xué)學(xué)報(bào),2016,36(3):1080-1091.[ZHAO Chen,PENG Shu-chuan,CHEN Tian-hu,etal.Ecological Risk Assessment of Sediment-Associated Permethrin and Esfenvalerate in Chaohu Lake and Taihu Lake Watersheds[J].Acta Scientiae Circumstantiae,2016,36(3):1080-1091.]

    [28] 梁霞,周軍英,李建宏,等.物種敏感度分布法(SSD)在農(nóng)藥水質(zhì)基準(zhǔn)推導(dǎo)中的應(yīng)用[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2015,31(3):398-405.[LIANG Xia,ZHOU Jun-ying,LI Jian-hong,etal.Application of Species Sensitivity Distribution (SSD) to Derivation of Water Quality Criteria for Pesticides[J].Journal of Ecology and Rural Environment,2015,31(3):398-405.]

    [29] 張瀚文.四溴雙酚A和4-叔丁基苯酚的多物種水生態(tài)危害評(píng)估[D].長(zhǎng)春:東北師范大學(xué),2014.[ZHANG Han-wen.Multi-Species Water Ecological Risk Assessment of TBBPA and 4-TBP[D].Changchun:Northeast Normal University,2014.]

    EcologicalRiskAssessmentofTriclocarbaninFreshWaterofChinabySpeciesSensitivityDistribution.

    WANG Zhen, YANG Xian-hai, FAN De-ling, GUO Min, LIU Ji-ning, SHI Li-li

    (Nanjing Institute of Environmental Sciences, Ministry of Environmental Protection, Nanjing 210042, China)

    The species sensitivity distributions (SSD) method was used to plot SSD curves of triclocarban (TCC) to freshwater organisms, and acute hazardous concentration for 5% of the species (HC5) of TCC was calculated. The log-normal distribution model was used to work out acuteρ(HC5)=3.85 μg·L-1, and based on AF=3 and acute to chronic ratios=39.3, chronic PNEC was worked out to be 32.7 ng·L-1. Risk quotient method was used in assessing ecological risk of triclocarban in 21 typical reported survey points of surface water, and the results of the assessment showed that the percentage for high risk, medium risk and low risk of TCC were 28.6%, 47.6% and 23.8%. The ecological risk of TCC for surface water in China deserve more attention.

    triclocarban; risk quotient; species sensitivity distributions (SSD); surface water

    2016-11-28

    科技基礎(chǔ)性工作專項(xiàng)(2015FY110900)

    ① 通信作者E-mail: ljn@nies.org

    X826

    A

    1673-4831(2017)10-0921-07

    10.11934/j.issn.1673-4831.2017.10.008

    汪貞(1982—),男,湖北云夢(mèng)人,助理研究員,博士,主要研究方向?yàn)榛瘜W(xué)品生態(tài)毒理與風(fēng)險(xiǎn)評(píng)估。E-mail: wangzhen@nies.org

    責(zé)任編輯: 李祥敏)

    猜你喜歡
    三氯淡水毒性
    不簡(jiǎn)單!一口普通的淡水蝦塘,他們竟能做到畝產(chǎn)2000多斤,獲利3萬(wàn)多/畝
    漲瘋了!碘漲50%,三氯漲超30%,溶劑漲超250%……消毒劑要漲價(jià)了
    鯨豚趣多多之它們愛淡水
    跳到海里喝淡水
    動(dòng)物之最——毒性誰(shuí)最強(qiáng)
    歐盟食品安全局:三氯蔗糖無(wú)致癌風(fēng)險(xiǎn)
    三氯生對(duì)4種水生生物的急性毒性研究
    RGD肽段連接的近紅外量子點(diǎn)對(duì)小鼠的毒性作用
    當(dāng)冷鹽水遇見溫淡水
    PM2.5中煤煙聚集物最具毒性
    人妻夜夜爽99麻豆av| 国产亚洲5aaaaa淫片| 成人亚洲精品av一区二区| 床上黄色一级片| 欧美日韩一区二区视频在线观看视频在线 | 免费观看的影片在线观看| 亚洲国产精品sss在线观看| 成人综合一区亚洲| ponron亚洲| 中文字幕久久专区| 精品99又大又爽又粗少妇毛片| 久久精品夜色国产| 偷拍熟女少妇极品色| 边亲边吃奶的免费视频| 久久精品熟女亚洲av麻豆精品 | 国产精品三级大全| 中文欧美无线码| 色网站视频免费| 精品人妻一区二区三区麻豆| 欧美zozozo另类| 亚洲熟妇中文字幕五十中出| 日韩精品青青久久久久久| 亚洲欧美成人精品一区二区| 国产精品嫩草影院av在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人a∨麻豆精品| 久久久久久久午夜电影| 美女国产视频在线观看| 亚洲无线观看免费| 男人和女人高潮做爰伦理| 久久久久久久久久黄片| 99视频精品全部免费 在线| 国产亚洲av片在线观看秒播厂 | 午夜精品在线福利| 欧美精品国产亚洲| 国产永久视频网站| 久久久久久久久中文| 国产精品国产三级国产专区5o| 18禁在线播放成人免费| 欧美日韩精品成人综合77777| 十八禁网站网址无遮挡 | 深夜a级毛片| 99视频精品全部免费 在线| 国产伦一二天堂av在线观看| 午夜视频国产福利| 国产色爽女视频免费观看| 亚洲人成网站在线播| 秋霞在线观看毛片| 亚洲国产欧美在线一区| 久久精品人妻少妇| 熟女电影av网| 97人妻精品一区二区三区麻豆| 国产91av在线免费观看| 久99久视频精品免费| 五月天丁香电影| 九九久久精品国产亚洲av麻豆| 色网站视频免费| 精品亚洲乱码少妇综合久久| 大香蕉97超碰在线| 少妇的逼好多水| 2021天堂中文幕一二区在线观| 91精品一卡2卡3卡4卡| 六月丁香七月| 精品午夜福利在线看| 国产极品天堂在线| 日韩精品青青久久久久久| 一区二区三区高清视频在线| 国产精品久久久久久久电影| 久久久久国产网址| 欧美激情国产日韩精品一区| 水蜜桃什么品种好| 亚洲精品久久午夜乱码| 97在线视频观看| 国产成年人精品一区二区| 免费av毛片视频| 老司机影院毛片| 亚洲性久久影院| 国产又色又爽无遮挡免| 亚洲不卡免费看| 热99国产精品久久久久久7| 丝袜美足系列| 久久精品久久精品一区二区三区| 男女免费视频国产| 精品人妻在线不人妻| 丰满饥渴人妻一区二区三| 美女福利国产在线| 久久久久视频综合| 亚洲一级一片aⅴ在线观看| 国产在线免费精品| 国产精品麻豆人妻色哟哟久久| 国产亚洲一区二区精品| 国产在视频线精品| 国产1区2区3区精品| 啦啦啦中文免费视频观看日本| 亚洲色图综合在线观看| 黄色怎么调成土黄色| 精品少妇久久久久久888优播| 黑人欧美特级aaaaaa片| 精品第一国产精品| 成人国产av品久久久| 亚洲av欧美aⅴ国产| 成人免费观看视频高清| 成人午夜精彩视频在线观看| 热re99久久精品国产66热6| 9191精品国产免费久久| 性色avwww在线观看| 久久精品国产综合久久久| 看免费av毛片| 最黄视频免费看| 永久网站在线| 久久精品国产亚洲av涩爱| 久久久国产欧美日韩av| 亚洲五月色婷婷综合| videossex国产| 日日爽夜夜爽网站| 五月伊人婷婷丁香| 天天操日日干夜夜撸| 男女下面插进去视频免费观看| av网站免费在线观看视频| 久久久久久久久久久久大奶| 久久99精品国语久久久| 一区在线观看完整版| 精品一区二区三区四区五区乱码 | 精品国产乱码久久久久久男人| 久久女婷五月综合色啪小说| 亚洲av电影在线进入| 久久久精品免费免费高清| 成年av动漫网址| 极品少妇高潮喷水抽搐| 日韩欧美精品免费久久| 在线观看美女被高潮喷水网站| 国产精品欧美亚洲77777| 丰满乱子伦码专区| 免费不卡的大黄色大毛片视频在线观看| 国产福利在线免费观看视频| 九草在线视频观看| 日本91视频免费播放| 蜜桃国产av成人99| 不卡av一区二区三区| 黄色毛片三级朝国网站| 久久久精品94久久精品| 日本vs欧美在线观看视频| 免费在线观看黄色视频的| 人体艺术视频欧美日本| 免费在线观看完整版高清| 十八禁网站网址无遮挡| 90打野战视频偷拍视频| 色网站视频免费| 国产福利在线免费观看视频| 国产黄色免费在线视频| 日本av免费视频播放| 美女脱内裤让男人舔精品视频| www.av在线官网国产| 成人免费观看视频高清| 人体艺术视频欧美日本| 精品福利永久在线观看| 纵有疾风起免费观看全集完整版| 亚洲中文av在线| 精品酒店卫生间| 在线精品无人区一区二区三| 国产又色又爽无遮挡免| 捣出白浆h1v1| 色网站视频免费| 永久网站在线| 成人免费观看视频高清| 亚洲,欧美精品.| 男女无遮挡免费网站观看| 涩涩av久久男人的天堂| 久热久热在线精品观看| 久久精品夜色国产| 99久久精品国产国产毛片| av在线播放精品| 国产一区二区 视频在线| 国产精品免费大片| 亚洲精品乱久久久久久| 大片电影免费在线观看免费| 久久久久精品人妻al黑| 国产深夜福利视频在线观看| 超碰97精品在线观看| 一边摸一边做爽爽视频免费| 我要看黄色一级片免费的| 啦啦啦视频在线资源免费观看| 国产综合精华液| 制服人妻中文乱码| 丰满迷人的少妇在线观看| 亚洲内射少妇av| 中文精品一卡2卡3卡4更新| 国产亚洲欧美精品永久| 综合色丁香网| 大片免费播放器 马上看| a 毛片基地| 成人二区视频| tube8黄色片| 午夜免费男女啪啪视频观看| 国产一区二区激情短视频 | 一区二区三区激情视频| 自线自在国产av| 成年动漫av网址| 在线观看人妻少妇| 国产精品久久久av美女十八| 亚洲,欧美精品.| av有码第一页| 熟女电影av网| 超色免费av| 99久久中文字幕三级久久日本| 亚洲内射少妇av| 亚洲精品av麻豆狂野| 国产亚洲最大av| 夫妻午夜视频| 色哟哟·www| 国产激情久久老熟女| 高清av免费在线| 国产精品免费大片| 日本猛色少妇xxxxx猛交久久| 高清av免费在线| 欧美中文综合在线视频| 国产亚洲欧美精品永久| 又粗又硬又长又爽又黄的视频| 777米奇影视久久| 色吧在线观看| 国产精品三级大全| 这个男人来自地球电影免费观看 | 国产探花极品一区二区| 人人妻人人澡人人爽人人夜夜| 午夜日本视频在线| 搡女人真爽免费视频火全软件| 777米奇影视久久| 日日摸夜夜添夜夜爱| 青春草亚洲视频在线观看| 亚洲精品国产一区二区精华液| 欧美另类一区| 天天躁日日躁夜夜躁夜夜| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久精品古装| 999精品在线视频| 三上悠亚av全集在线观看| 国产淫语在线视频| 精品一区二区三卡| 欧美97在线视频| 国产在线免费精品| 久久久亚洲精品成人影院| 国产一区二区三区av在线| 成人手机av| 26uuu在线亚洲综合色| 久久久精品国产亚洲av高清涩受| 久久久久久伊人网av| 18禁裸乳无遮挡动漫免费视频| 啦啦啦在线观看免费高清www| 丰满乱子伦码专区| 青青草视频在线视频观看| 日本av手机在线免费观看| 日本av手机在线免费观看| 国产精品久久久久久av不卡| 亚洲综合精品二区| 国产精品嫩草影院av在线观看| 一级a爱视频在线免费观看| 亚洲欧美日韩另类电影网站| 成年人午夜在线观看视频| 九色亚洲精品在线播放| 国产精品久久久久久久久免| 久久这里有精品视频免费| 亚洲人成77777在线视频| 亚洲av国产av综合av卡| 一本色道久久久久久精品综合| 黄片播放在线免费| 亚洲欧美精品综合一区二区三区 | 999精品在线视频| 99久久人妻综合| 欧美黄色片欧美黄色片| 母亲3免费完整高清在线观看 | 久久久久久久久免费视频了| av卡一久久| 一级黄片播放器| av在线app专区| 天天躁狠狠躁夜夜躁狠狠躁| 精品亚洲成a人片在线观看| 赤兔流量卡办理| 中文字幕人妻丝袜一区二区 | 国产xxxxx性猛交| 岛国毛片在线播放| 纵有疾风起免费观看全集完整版| 亚洲人成电影观看| 国产麻豆69| 啦啦啦在线观看免费高清www| 国产日韩一区二区三区精品不卡| 免费观看无遮挡的男女| av电影中文网址| 亚洲欧美日韩另类电影网站| 亚洲精品国产av蜜桃| 天天影视国产精品| 国产精品不卡视频一区二区| 久久久国产精品麻豆| 男女下面插进去视频免费观看| 久久99蜜桃精品久久| 看非洲黑人一级黄片| 777米奇影视久久| 人人妻人人添人人爽欧美一区卜| 丰满少妇做爰视频| 最近的中文字幕免费完整| 一级爰片在线观看| 伦精品一区二区三区| 欧美精品人与动牲交sv欧美| 亚洲一级一片aⅴ在线观看| 国产免费视频播放在线视频| 性色av一级| 久久毛片免费看一区二区三区| 天堂中文最新版在线下载| 亚洲伊人久久精品综合| 男男h啪啪无遮挡| 黄色 视频免费看| 亚洲三级黄色毛片| 免费大片黄手机在线观看| 午夜日本视频在线| 自拍欧美九色日韩亚洲蝌蚪91| 狂野欧美激情性bbbbbb| 宅男免费午夜| 久久久久久久久久久久大奶| 欧美国产精品一级二级三级| 亚洲精品第二区| 久久久欧美国产精品| 99香蕉大伊视频| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 亚洲人成网站在线观看播放| 精品国产露脸久久av麻豆| 精品国产一区二区久久| 观看美女的网站| 日韩一区二区三区影片| 黄片小视频在线播放| 母亲3免费完整高清在线观看 | 日韩av免费高清视频| 毛片一级片免费看久久久久| 日韩制服丝袜自拍偷拍| 看十八女毛片水多多多| 看非洲黑人一级黄片| 9色porny在线观看| a级毛片黄视频| 伊人久久国产一区二区| 国产精品不卡视频一区二区| 成年美女黄网站色视频大全免费| 看非洲黑人一级黄片| 欧美精品av麻豆av| 欧美人与性动交α欧美软件| 晚上一个人看的免费电影| av视频免费观看在线观看| 美女脱内裤让男人舔精品视频| 如日韩欧美国产精品一区二区三区| 国产精品人妻久久久影院| 亚洲,一卡二卡三卡| 狂野欧美激情性bbbbbb| 少妇精品久久久久久久| 看十八女毛片水多多多| 成人手机av| 国产精品欧美亚洲77777| 亚洲精品国产av成人精品| 香蕉精品网在线| 另类亚洲欧美激情| 日本免费在线观看一区| 免费在线观看黄色视频的| 成人漫画全彩无遮挡| 欧美日韩av久久| 欧美av亚洲av综合av国产av | 精品一品国产午夜福利视频| 欧美xxⅹ黑人| 性少妇av在线| 99久久综合免费| 丁香六月天网| 纵有疾风起免费观看全集完整版| 一边亲一边摸免费视频| 国产 精品1| 亚洲国产精品999| 69精品国产乱码久久久| 午夜福利视频精品| 久久久国产欧美日韩av| 国产福利在线免费观看视频| 高清av免费在线| 国产精品久久久久成人av| 久久久久久久亚洲中文字幕| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看| 日韩免费高清中文字幕av| 亚洲欧美色中文字幕在线| 最近最新中文字幕大全免费视频 | 国产 精品1| 久久国产精品男人的天堂亚洲| 欧美日韩综合久久久久久| 大码成人一级视频| 曰老女人黄片| kizo精华| 如日韩欧美国产精品一区二区三区| 久久精品久久精品一区二区三区| 性高湖久久久久久久久免费观看| 交换朋友夫妻互换小说| 久久影院123| 亚洲精品国产av蜜桃| 国产在线视频一区二区| 亚洲色图 男人天堂 中文字幕| 国产日韩一区二区三区精品不卡| 久久久久国产精品人妻一区二区| 成年女人在线观看亚洲视频| 亚洲一区二区三区欧美精品| 久久久久人妻精品一区果冻| 激情视频va一区二区三区| 国产男女超爽视频在线观看| 多毛熟女@视频| 欧美成人午夜免费资源| 满18在线观看网站| 久久久久人妻精品一区果冻| 蜜桃国产av成人99| 亚洲成人av在线免费| 日本猛色少妇xxxxx猛交久久| 国产乱来视频区| 久久鲁丝午夜福利片| 777米奇影视久久| 一二三四在线观看免费中文在| 午夜福利乱码中文字幕| 亚洲精品久久久久久婷婷小说| 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 免费女性裸体啪啪无遮挡网站| 国产视频首页在线观看| 国产欧美日韩一区二区三区在线| 久久鲁丝午夜福利片| 亚洲伊人久久精品综合| 美女国产高潮福利片在线看| 午夜福利,免费看| 欧美日韩精品成人综合77777| 波野结衣二区三区在线| 少妇人妻久久综合中文| 老司机亚洲免费影院| 成年女人在线观看亚洲视频| 美女大奶头黄色视频| 国产高清国产精品国产三级| 一区二区三区激情视频| 赤兔流量卡办理| 91国产中文字幕| 亚洲,一卡二卡三卡| 亚洲伊人色综图| 黄频高清免费视频| 国产精品免费大片| 老汉色∧v一级毛片| 精品国产一区二区三区久久久樱花| 国产97色在线日韩免费| 国产精品免费视频内射| 美女高潮到喷水免费观看| a 毛片基地| 香蕉精品网在线| 久久久久久久久久人人人人人人| 七月丁香在线播放| 少妇 在线观看| 伦理电影大哥的女人| 国产亚洲欧美精品永久| 蜜桃国产av成人99| 狠狠婷婷综合久久久久久88av| 成人18禁高潮啪啪吃奶动态图| 晚上一个人看的免费电影| 久久久久久久久久久免费av| 十分钟在线观看高清视频www| 免费日韩欧美在线观看| 考比视频在线观看| 午夜日本视频在线| 最近中文字幕高清免费大全6| 久久精品国产亚洲av天美| av.在线天堂| 国产精品 欧美亚洲| 啦啦啦啦在线视频资源| www.av在线官网国产| 七月丁香在线播放| 水蜜桃什么品种好| 久久精品国产亚洲av涩爱| 久久这里只有精品19| 成年人午夜在线观看视频| 日韩在线高清观看一区二区三区| 日本午夜av视频| 亚洲欧美一区二区三区久久| 成年动漫av网址| 如日韩欧美国产精品一区二区三区| 嫩草影院入口| 日韩精品免费视频一区二区三区| 91国产中文字幕| videossex国产| 午夜免费观看性视频| 亚洲熟女精品中文字幕| 美国免费a级毛片| 另类亚洲欧美激情| 日韩伦理黄色片| 在线 av 中文字幕| 亚洲成av片中文字幕在线观看 | 国产成人免费观看mmmm| 一边摸一边做爽爽视频免费| 美女脱内裤让男人舔精品视频| 久久久久久久久免费视频了| 免费少妇av软件| 日韩伦理黄色片| 少妇人妻精品综合一区二区| 波野结衣二区三区在线| 中国国产av一级| 欧美在线黄色| 亚洲美女视频黄频| 日韩中文字幕欧美一区二区 | 91久久精品国产一区二区三区| 精品亚洲乱码少妇综合久久| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 欧美老熟妇乱子伦牲交| 成年女人毛片免费观看观看9 | 最新中文字幕久久久久| 大香蕉久久网| 亚洲精品国产色婷婷电影| 天天躁狠狠躁夜夜躁狠狠躁| 最新中文字幕久久久久| 国产精品免费视频内射| 精品国产乱码久久久久久男人| 黄色怎么调成土黄色| h视频一区二区三区| 一级毛片 在线播放| 女人被躁到高潮嗷嗷叫费观| 一级片免费观看大全| 精品卡一卡二卡四卡免费| 日韩欧美一区视频在线观看| 免费播放大片免费观看视频在线观看| 亚洲国产精品成人久久小说| 捣出白浆h1v1| 色哟哟·www| av线在线观看网站| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 人妻一区二区av| 在线观看三级黄色| 午夜精品国产一区二区电影| 男的添女的下面高潮视频| 91精品伊人久久大香线蕉| 成人午夜精彩视频在线观看| 日本vs欧美在线观看视频| 激情视频va一区二区三区| 亚洲av日韩在线播放| 一二三四中文在线观看免费高清| 久久这里有精品视频免费| 国产一区二区在线观看av| 国产亚洲一区二区精品| 寂寞人妻少妇视频99o| 久久 成人 亚洲| 日韩在线高清观看一区二区三区| 成人国产av品久久久| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 水蜜桃什么品种好| 熟女电影av网| 一区二区三区四区激情视频| 欧美变态另类bdsm刘玥| 成人国语在线视频| 18禁观看日本| 在线观看免费日韩欧美大片| 亚洲av免费高清在线观看| 国产精品.久久久| 在线观看美女被高潮喷水网站| 一区二区日韩欧美中文字幕| 亚洲av男天堂| 国产欧美日韩一区二区三区在线| 日本色播在线视频| 一区二区三区激情视频| 国产又爽黄色视频| 看非洲黑人一级黄片| 美女视频免费永久观看网站| 欧美日韩成人在线一区二区| 亚洲精品成人av观看孕妇| 麻豆av在线久日| 少妇人妻精品综合一区二区| 男女免费视频国产| 9191精品国产免费久久| 精品少妇黑人巨大在线播放| 久久午夜综合久久蜜桃| 国产xxxxx性猛交| 亚洲美女黄色视频免费看| 下体分泌物呈黄色| 亚洲国产av影院在线观看| 久久久久久久大尺度免费视频| 91aial.com中文字幕在线观看| 少妇人妻精品综合一区二区| 中文天堂在线官网| av不卡在线播放| 久久女婷五月综合色啪小说| 一区二区三区四区激情视频| 国产探花极品一区二区| 观看美女的网站| 人人妻人人爽人人添夜夜欢视频| 在线观看国产h片| 国产成人午夜福利电影在线观看| 多毛熟女@视频| 国产精品香港三级国产av潘金莲 | a 毛片基地| 最近最新中文字幕大全免费视频 | 国产成人一区二区在线| 丰满乱子伦码专区| 精品一品国产午夜福利视频| 国产一区亚洲一区在线观看| videos熟女内射| 爱豆传媒免费全集在线观看| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| 欧美成人午夜免费资源| 人人妻人人澡人人看| 自线自在国产av| 精品久久久久久电影网| 精品人妻一区二区三区麻豆| av片东京热男人的天堂| 高清在线视频一区二区三区| 桃花免费在线播放| tube8黄色片| 熟女少妇亚洲综合色aaa.| 久久这里有精品视频免费| 久久久精品94久久精品| 亚洲视频免费观看视频| 国产成人av激情在线播放| 九草在线视频观看|