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

    應(yīng)用概率物種敏感度分布法研究太湖銅水生生物水質(zhì)基準(zhǔn)

    2015-06-27 05:54:12侯俊趙芊淵王超王沛芳苗令占呂博文
    生態(tài)毒理學(xué)報(bào) 2015年1期
    關(guān)鍵詞:水生太湖敏感度

    侯俊,趙芊淵,王超,王沛芳,苗令占,呂博文

    1. 河海大學(xué)教育部淺水湖泊綜合治理與資源開發(fā)重點(diǎn)實(shí)驗(yàn)室,南京 210098 2. 河海大學(xué)環(huán)境學(xué)院,南京 210098

    應(yīng)用概率物種敏感度分布法研究太湖銅水生生物水質(zhì)基準(zhǔn)

    侯俊1,2,*,趙芊淵1,2,王超1,2,王沛芳1,2,苗令占1,2,呂博文1,2

    1. 河海大學(xué)教育部淺水湖泊綜合治理與資源開發(fā)重點(diǎn)實(shí)驗(yàn)室,南京 210098 2. 河海大學(xué)環(huán)境學(xué)院,南京 210098

    運(yùn)用概率物種敏感度分布法獲得太湖水體中銅的急性水質(zhì)基準(zhǔn)值為14.57 μg·L-1,慢性水質(zhì)基準(zhǔn)值為3.26 μg·L-1;不同類別物種敏感性存在差異,無脊椎動物較脊椎動物更敏感,甲殼類敏感性大于魚類。概率物種敏感度分布法與傳統(tǒng)的物種敏感度分布法相比,更全面合理地考慮多種毒性效應(yīng),曲線擬合效果好,受數(shù)據(jù)量大小影響較小,結(jié)果更加穩(wěn)定。研究結(jié)論可為銅水質(zhì)標(biāo)準(zhǔn)的修訂和太湖流域水環(huán)境管理提供技術(shù)支持。

    銅;水生生物;急性毒性;慢性毒性;水質(zhì)基準(zhǔn);太湖;概率物種敏感度分布

    銅是生命體必需的微量元素,同時(shí)也是水體中重金屬污染的主要元素之一,水體中的銅主要來自大氣沉降、農(nóng)業(yè)徑流,為控制藻類水華而使用的硫酸銅以及含銅工業(yè)污廢水的直接排放。通常情況下,銅對水生生物的毒性大于它對人類和其他陸生生物的毒性[1]。為保護(hù)水體中水生生物,許多國家和地區(qū)都已開展了銅水生生物基準(zhǔn)研究,如美國、歐盟、加拿大、澳大利亞、荷蘭和新西蘭等,采用的方法主要是評價(jià)因子法、毒性百分?jǐn)?shù)排序法和物種敏感度分布法,基于水生生物基準(zhǔn)制定和頒布了一系列環(huán)境質(zhì)量標(biāo)準(zhǔn)[1-2]。我國銅的水生生物基準(zhǔn)的研究還比較少,水質(zhì)標(biāo)準(zhǔn)主要是參照國外水質(zhì)基準(zhǔn)或者標(biāo)準(zhǔn)制定的,例如我國頒布的GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中Ⅰ類和Ⅱ類Cu水質(zhì)標(biāo)準(zhǔn)限值分別為10 μg·L-1和1 000 μg·L-1[1],因此,我國水生生物的保護(hù)還缺乏可靠的科學(xué)依據(jù)。

    目前,國際上推導(dǎo)水質(zhì)基準(zhǔn)的主流方法為物種敏感度分布法(species sensitivity distribution, SSD),該方法最初是由Kooijman提出,后來很多學(xué)者對其進(jìn)行了改進(jìn)。該理論認(rèn)為:不同門類生物由于生活史、生理構(gòu)造、行為特征和地理分布等不同而產(chǎn)生了差異性,其在毒理學(xué)上反映為不同物種對污染物有不同的劑量-效應(yīng)響應(yīng)關(guān)系,即在結(jié)構(gòu)復(fù)雜的生態(tài)系統(tǒng)中,不同物種對某一脅迫因素(如有毒化學(xué)物質(zhì))的敏感性(如無觀察效應(yīng)濃度NOEC或半數(shù)致死濃度LC50等毒性數(shù)據(jù))服從一定的分布,而可獲得的毒性數(shù)據(jù)被認(rèn)為是來自于這個(gè)分布的樣本,可用來估算該分布的參數(shù)。前提假設(shè)是這些物種的選擇均具有隨機(jī)性,且能夠代表給定生態(tài)系統(tǒng)的群落結(jié)構(gòu)[3-10]。

    然而,在物種敏感度分布法(SSD)的方法理論中,對于曲線的擬合,目前研究還無法表明某一特定的分布模型適用于所有情況下毒性數(shù)據(jù)集的擬合[5],因而導(dǎo)致SSD法曲線擬合具有較大主觀性、擬合效果差、結(jié)果不穩(wěn)定等不足。為了解決上述問題,本文考慮一種新的物種敏感度曲線擬合方法,即概率物種敏感度分布法(probabilistic species sensitivity distribution, PSSD)。

    太湖是我國第三大淡水湖,在區(qū)域經(jīng)濟(jì)和社會發(fā)展中具有舉足輕重的地位,銅對太湖水生生物的毒性效應(yīng)引起了人們的廣泛關(guān)注[1,11-12]。余海靜等[13]在2012年對太湖水體采樣分析得到太湖水體總Cu含量為1.298~4.485 μg·L-1,平均值為2.22 μg·L-1,最小值位于東部湖區(qū),最大值位于北部湖區(qū),而且Cu在太湖沉積物中也有一定的富集。本研究結(jié)合太湖生物區(qū)系特點(diǎn),選取太湖水體中代表物種,應(yīng)用概率物種敏感度分布法推導(dǎo)太湖銅的水生生物水質(zhì)基準(zhǔn),以期為我國淡水水體中重金屬基準(zhǔn)的建立、重金屬污染監(jiān)測評價(jià)和防控治理提供參考。

    1 研究方法(Methodology)

    1.1 概率物種敏感度分布法(PSSD)曲線構(gòu)建

    先確定單個(gè)物種對某種污染物的敏感度,然后將這些單個(gè)物種的敏感度組合在一起得到概率物種敏感度分布(PSSD),具體構(gòu)建過程見圖1。使用的軟件是R語言(R Development Core Team 2008)。

    (1) 收集整理所有可獲得的所要研究的生態(tài)系統(tǒng)的毒性數(shù)據(jù)。將某一單個(gè)物種的所有符合要求的毒性終點(diǎn)值組合,以獲得該物種敏感度的概率密度函數(shù)(probability density function, PDF),這樣考慮了單個(gè)物種所有毒性數(shù)據(jù),最大限度地利用已獲得的所有有效數(shù)據(jù),而不是僅僅考慮一個(gè)確定的幾何均值。

    ① 對于某一物種對某一污染物只有一個(gè)符合要求的毒性終點(diǎn)值(a)的情況,我們考慮三角形概率密度分布(triangular PDF),構(gòu)建三角形概率密度分布時(shí),使用變異系數(shù)(coefficient of variation, CV)來獲得可能的最小值(a')和最大值(a")。美國環(huán)保局技術(shù)導(dǎo)則(US EPA 2000)提出30%的變異系數(shù)能夠大致代表實(shí)驗(yàn)室內(nèi)部平均差異[14],本文中,我們假定單個(gè)毒性終點(diǎn)以及復(fù)合分布邊界值的變異系數(shù)均為30%。如圖1a。

    ② 對于某一物種對某一污染物有2個(gè)符合要求的毒性終點(diǎn)值(a和b)的情況,我們考慮梯形分布(trapezoidal distribution)。這種梯形分布以均勻分布為基礎(chǔ),這2個(gè)毒性值作為均勻分布的區(qū)間端點(diǎn),每邊再補(bǔ)加上由變異系數(shù)確定的可能最值( a'和b')構(gòu)成的三角形分布。如圖1b。

    ③ 對于有多個(gè)毒性終點(diǎn)值的情況,毒性數(shù)據(jù)的最低值(a)和第二低值(b)構(gòu)成第1個(gè)區(qū)間的2個(gè)端點(diǎn),第二低值(b)和第三低值(c)構(gòu)成第2個(gè)區(qū)間的2個(gè)端點(diǎn),依此類推,有n個(gè)毒性終點(diǎn)值的情況,就有(n-1)個(gè)均勻分布區(qū)間,每邊再補(bǔ)加上由變異系數(shù)確定的可能最值(a'和c')構(gòu)成的三角形分布。不同長度的區(qū)間,概率是一樣的。也就是說在第1個(gè)區(qū)間和在第n個(gè)區(qū)間上取值的概率是相等的。如圖1c。

    (2) 將這些單個(gè)的PDF轉(zhuǎn)化成該物種的PDF。從這些單個(gè)的PDF中隨機(jī)抽取同樣數(shù)量的值(本文中隨機(jī)抽樣100 000個(gè)值,取值過少會使得結(jié)果不穩(wěn)定,取100 000,穩(wěn)定性能達(dá)到很好),即在每個(gè)小區(qū)間隨機(jī)取100 000個(gè)值(基于均勻分布),對于只有一個(gè)毒性值的情況,在該區(qū)間隨機(jī)取100 000個(gè)值(基于三角形分布),對于2個(gè)及2個(gè)以上毒性值的情況中,兩邊的小區(qū)間隨機(jī)取3 000個(gè)值(基于三角形分布),將這些值組合得到該物種的概率密度分布,即單一物種的概率密度分布PDF。如圖1d。

    圖1 概率物種敏感度分布(PSSD)的構(gòu)建過程Fig. 1 The development procedure of probabilistic species sensitivity distribution (PSSD)注:(a)物種只有1個(gè)毒性值的情況,a為該物種的毒性值,a’為可能最小值,a’=a×(1-CV)=0.7a,a’’為可能最大值,a’’=a×(1+CV)=1.3a;(b)物種有2個(gè)毒性值的情況,a、b為該物種的2個(gè)毒性值,a’、b’為可能最小值和可能最大值,a’=0.7a,b’=1.3b;(c)物種有3個(gè)毒性值的情況,a、b和c為該物種的3個(gè)毒性值,a’、c’為可能最小值和可能最大值,a’=0.7a,c’=1.3c;(d)物種有5個(gè)毒性值的情況,a、b、c、d和e為該物種的5個(gè)毒性值,a’、e’為可能最小值和可能最大值,a’=0.7a,e’=1.3e;(e)PSSD曲線。Note: (a) species with only one single toxic value (a), the possible minimum is a’, a’=a×(1-CV)=0.7a, the possible maximum is a’’, a’’=a×(1+CV)=1.3a; (b) species with two toxic values (a and b), the possible minimum and maximum are a' and b', a'=0.7a, b'=1.3b; (c) species with three toxic values (a, b and c), the possible minimum and maximum are a' and c', a'=0.7a, c'=1.3c; (d) species with five toxic values (a, b, c, d and e), the possible minimum and maximum are a' and e', a'=0.7a, e'=1.3e; (e) PSSD curves.

    (3) 將所有物種的概率密度分布組合成所要研究的環(huán)境體系的PSSD。如圖1e所示,從每一單個(gè)物種的概率密度分布PDF中隨機(jī)抽取同樣數(shù)量的值(本文中此處隨機(jī)抽取10 000個(gè)值),將這些值組合得到所要研究的環(huán)境體系的概率物種敏感度分布PSSD。圖1e中黑色粗實(shí)線PSSD即是由物種1、物種2、物種3、物種4和物種5的PDF組合成的[15-17]。

    1.2 毒性數(shù)據(jù)獲取

    毒性數(shù)據(jù)主要來自美國環(huán)保局ECOTOX數(shù)據(jù)庫(http: / /cfpub.epa.gov /ecotox /)和中國知網(wǎng)(http: / /www. cnki. net /)收錄的文獻(xiàn),以及Elsevier、Wiley Online Library、SpringerLink、IWA(International Water Association)等上發(fā)表的文獻(xiàn)。

    將所獲得的毒性數(shù)據(jù)按照以下原則篩選:盡量選擇太湖物種,江浙地區(qū)以及長江廣泛分布的物種也可以視作太湖物種;毒性試驗(yàn)方法與相關(guān)標(biāo)準(zhǔn)測試方法一致,如經(jīng)濟(jì)合作與發(fā)展組織或美國材料與試驗(yàn)協(xié)會等發(fā)布的毒性試驗(yàn)方法;實(shí)驗(yàn)用水為淡水或自來水,而實(shí)驗(yàn)用水為海水和去離子/蒸餾水的數(shù)據(jù)不可用;水生生物的急性毒性試驗(yàn)選擇24~96 h的LC50或EC50,水蚤或其他枝角類動物,選擇24 h的LC50或EC50,蚊蟲或其他昆蟲的毒性效應(yīng)終點(diǎn)以48 h的EC50或LC50表示,其他生物毒性效應(yīng)終點(diǎn)以96 h的EC50或LC50表示,如未獲得相應(yīng)24 h和48 h的毒性值,則選擇96 h的毒性值;水生生物慢性毒性終點(diǎn)使用最低可見效應(yīng)濃度LOEC和NOEC[18-20]。

    1.3 物種分組和數(shù)據(jù)處理

    為了從生態(tài)系統(tǒng)不同層次研究,將毒性數(shù)據(jù)根據(jù)生物營養(yǎng)級水平分類考慮,將收集的數(shù)據(jù)分3種情況考慮:①全部物種;②全部物種分為脊椎動物和無脊椎動物;③脊椎動物分為魚類和兩棲類(甲魚等爬行動物歸為兩棲類),無脊椎動物分為甲殼類、昆蟲類和其他無脊椎動物。獲得符合要求Cu的48個(gè)物種的急性毒性數(shù)據(jù),隸屬于8門23科34屬,見表1,符合3門8科的要求[18],太湖水生生物區(qū)系中魚類以鯉科為主,底棲動物中較多的是搖蚊科和顫蚓科,浮游動物中較多的是臂尾輪蟲科和溞科,浮游植物以藍(lán)藻門、綠藻門和硅藻門居多[,21-,22],所選物種符合太湖生物區(qū)系特征,可用來推導(dǎo)太湖Cu水質(zhì)基準(zhǔn)。一般PSSD曲線擬合數(shù)據(jù)量控制在10以上[18],甲殼類和魚類的物種數(shù)分別為10和17,可以單獨(dú)繪制PSSD曲線。表1為Cu急性毒性數(shù)據(jù)統(tǒng)計(jì)匯總。獲得符合要求Cu的6個(gè)物種的慢性毒性數(shù)據(jù),見表2。

    表1 不同物種類別的銅急性毒性數(shù)據(jù)統(tǒng)計(jì)匯總

    表2 銅的慢性毒性數(shù)據(jù)

    注:NOEC為無觀察效應(yīng)濃度,LOEC最低可見效應(yīng)濃度。

    Note: NOEC stands for no observed effect concentration; LOEC stands for lowest observed effect concentration.

    1.4 HC5和水質(zhì)基準(zhǔn)

    PSSD曲線上累積概率為5%對應(yīng)的毒性值即為HC5值,編輯R代碼,由R軟件直接輸出。部分代碼如下:

    C=sort(log(qq)) #濃度排序

    HC5,log=C[length(C)*0.05] #累積概率為5%處的對數(shù)濃度

    HC5=exp(HC5,log) #累積概率為5%處的濃度

    用急性毒性數(shù)據(jù)推導(dǎo)短期危害濃度(short term hazardous concentration, STHC5),慢性毒性數(shù)據(jù)推導(dǎo)長期危害濃度(long term hazardous concentration, LTHC5),數(shù)值上等于HC5。通常慢性毒性數(shù)據(jù)不足以構(gòu)建PSSD曲線,根據(jù)以下公式計(jì)算長期危害濃度:

    LTHC5=HC5,急性/FACR

    (1)

    式中,F(xiàn)ACR為最終急慢性比率,計(jì)算FACR需獲得至少3個(gè)科的水生生物急慢性比率(ACR),其中至少有1種是魚類,至少有1種是無脊椎動物,至少有1種是急性敏感淡水物種。由于未獲得所需物種的急慢性比率,本文采用美國環(huán)保局1995年水質(zhì)基準(zhǔn)報(bào)道大型溞(Daphnia magna)、鉤蝦(Gammarus pseudolimnaeus)和黑頭呆魚(Pimephales promelas)的ACR分別為2.418、3.297和11.2,得到Cu的FACR為4.470[58]。

    2 結(jié)果(Results)

    對符合要求的毒性數(shù)據(jù)用R軟件進(jìn)行數(shù)據(jù)處理,得到全部物種PSSD曲線,如圖2所示,圖中實(shí)線為全部物種的PSSD,圖中也繪出了原數(shù)據(jù)和幾何平均值,分別用空心和實(shí)心小圓表示。由PSSD曲線得到全部物種的HC5為14.57 μg·L-1,即短期危害濃度STHE5為14.57 μg·L-1,慢性毒性數(shù)據(jù)較少,無法擬合PSSD曲線,根據(jù)公式(1)得到長期危害濃度LTHE5為3.26 μg·L-1。

    圖2 銅的全部物種的概率物種敏感度分布曲線和物種敏感度分布曲線Fig. 2 Probabilistic species sensitivity distribution curve and species sensitivity distribution curve for total species in response to copper

    3 討論(Discussion)

    3.1 概率物種敏感度分布法與傳統(tǒng)的物種敏感度分布法比較

    傳統(tǒng)物種敏感度分布法是將篩選出的毒性數(shù)據(jù)(一個(gè)物種有多個(gè)毒性值,取幾何平均值)按從小到大的順序排列,并標(biāo)出相應(yīng)的序數(shù),求出對應(yīng)的累積概率P(P=R/(N+1),R為序數(shù),N為毒性數(shù)據(jù)總數(shù))。以累積概率為縱坐標(biāo),濃度值(或者濃度值的對數(shù))為橫坐標(biāo)作圖,選擇適當(dāng)?shù)那€擬合。由于并沒有一個(gè)特定的分布模型適用于所有毒性數(shù)據(jù)集的擬合,一般可先使用各種模型擬合進(jìn)行嘗試,然后根據(jù)決定系數(shù)、殘差平方和、F值和圖像等選擇最佳擬合曲線[4-9]。圖2中的虛線即為全部物種SSD曲線,選用quasibinomial(link=logit)擬合。

    本文的PSSD與傳統(tǒng)的SSD曲線擬合相比,用到了物種更多有效毒性數(shù)據(jù),從而考慮了水生生物對污染物多種不同響應(yīng),因?yàn)橹亟饘巽~對水生生物毒性效應(yīng)是多方面的,影響水生生物生長發(fā)育、新陳代謝、繁殖和種群數(shù)量等,不同季節(jié)、齡期等(毒性數(shù)據(jù)不充足)的水生生物對銅的敏感性也存在差異,環(huán)境中其他物質(zhì)的存在也會影響水生生物對銅的響應(yīng),制定水生生物水質(zhì)基準(zhǔn)也要從多方面對水生生物進(jìn)行保護(hù)。曲線形式由所獲得的有效毒性數(shù)據(jù)來決定,盡量考慮更多毒性效應(yīng),將這些確定的毒性數(shù)據(jù)轉(zhuǎn)化為概率信息(假設(shè)這些不同毒性效應(yīng)值以均等隨機(jī)概率出現(xiàn)),基于概率來利用數(shù)據(jù),擬合過程根據(jù)實(shí)際情況考慮三角形分布和均勻分布。圖2中全部物種的曲線擬合,易見PSSD的擬合效果好于傳統(tǒng)的SSD,幾乎所有單個(gè)值以及多個(gè)值的幾何平均值都在曲線上或是附近,偏離很小,較好地考慮了多種毒性效應(yīng),基于PSSD得到的水生生物水質(zhì)基準(zhǔn)能更全面合理地保護(hù)水生生物。傳統(tǒng)的SSD對于某單一物種有多個(gè)符合要求的毒性數(shù)據(jù)時(shí),使用單一確定的幾何平均值,如果這些多個(gè)值中有某一個(gè)值(也就是某種毒性效應(yīng)值)和其他值(其他各種毒性效應(yīng)值)相差很大,那么最終的幾何平均值(最終毒性值)則會偏向這個(gè)偏離很大的值,而PSSD通過均等的概率考慮所有值(不同的毒性效應(yīng)),使得偏離很大的離異值對最終結(jié)果不再起決定作用[1,5,8,15],得到的結(jié)果更合理。

    傳統(tǒng)的SSD法對擬合曲線的選擇也存在很大的主觀性,選擇不同的擬合曲線,所得到的結(jié)果也不同,在數(shù)據(jù)量較多的情況下,差別不大(圖2中實(shí)線和虛線分別為全部物種的PSSD和SSD);而在數(shù)據(jù)量較少的情況下,選不同曲線擬合得到的結(jié)果差異就比較明顯,大多數(shù)的情況下我們的數(shù)據(jù)量都不大,很容易出現(xiàn)這種情況。對甲殼類的毒性數(shù)據(jù)進(jìn)行傳統(tǒng)的SSD擬合,基于SSD(圖3)選擇Logistic擬合得到的HC5為6.83 μg·L-1,而基于SSD選擇Gompertz擬合得到的HC5為14.9 μg·L-1,另外2個(gè)曲線由于取值范圍原因沒能得到HC5值,基于傳統(tǒng)的SSD得到的2個(gè)結(jié)果相差一倍多,而圖3中4種模型的擬合效果綜合來看是不相上下的。而我們對甲殼類的毒性數(shù)據(jù)采用PSSD擬合(圖4b)得到HC5為1.30 μg·L-1。另外,PSSD曲線對于曲線兩端尾部數(shù)據(jù)也都擬合得很好,這樣我們所需要的處于曲線下端的HC5值就更加精確;而圖3采用4種不同模型分別擬合甲殼類SSD,從圖中能看出SSD曲線下端有較大差異,所以選不同模型獲得的HC5差異較大。

    PSSD較傳統(tǒng)的SSD更加穩(wěn)定,F(xiàn)adri等[15]的研究結(jié)果顯示應(yīng)用PSSD對三氯生(triclosan)毒性數(shù)據(jù)擬合,使用23個(gè)物種的毒性數(shù)據(jù)擬合得到的HC5為0.63 μg·L-1,增加23個(gè)新物種,對這46個(gè)物種的毒性數(shù)據(jù)進(jìn)行PSSD擬合獲得的HC5為0.66 μg·L-1,兩者及其接近;而Capdevielle等由同樣的23個(gè)物種毒性數(shù)據(jù)使用傳統(tǒng)的SSD選擇Log-logistic擬合得到的HC5為1.55 μg·L-1,Lyndall等對同樣的46個(gè)物種毒性數(shù)據(jù)使用傳統(tǒng)的SSD選用非參數(shù)的bootstrap得到HC5為0.8 μg·L-1,兩者相差一倍;46個(gè)物種基于傳統(tǒng)的SSD得到的0.8 μg·L-1與PSSD得到的0.63 μg·L-1和0.66 μg·L-1都相差較少。這一研究結(jié)果表明傳統(tǒng)的SSD更容易受到數(shù)據(jù)量大小的影響,PSSD相對而言更加穩(wěn)定,尤其是對于數(shù)據(jù)量較少時(shí)也能較好地?cái)M合數(shù)據(jù)。另外這也能說明PSSD能夠隨著數(shù)據(jù)量而不斷更新完善,而傳統(tǒng)的SSD增減數(shù)據(jù)量會有很大的變化。

    綜上可知,概率物種敏感度分布法較之傳統(tǒng)的物種敏感度分布法有較多優(yōu)勢,可以替代傳統(tǒng)的物種敏感度分布法或者作為補(bǔ)充。

    3.2 不同物種類別物種敏感度差異比較

    對銅的毒性數(shù)據(jù)進(jìn)行物種分類研究,脊椎動物有25種物種,無脊椎動物16種,魚類17種,甲殼類10種,可以單獨(dú)擬合PSSD曲線;兩棲類、昆蟲類和其他無脊椎動物物種數(shù)均少于10,無法擬合PSSD曲線。繪出不同物種類別PSSD曲線,如圖4所示,

    圖3 不同模型擬合銅的甲殼類物種敏感度分布曲線Fig. 3 Species sensitivity distribution curves for crustaceans in response to copper were simulated by different models

    顯然不同類別物種的PSSD曲線形狀存在較大差異。根據(jù)PSSD曲線得到脊椎動物、無脊椎動物、魚類和甲殼類的HC5分別為34.75 μg·L-1、0.18 μg·L-1、24.62 μg·L-1和1.30 μg·L-1,由表1得毒性均值分別為:389.982 μg·L-1、169.372 μg·L-1、384.013 μg·L-1和203.260 μg·L-1,比較HC5或毒性均值均可得水生生物不同物種類別對銅毒性敏感性大小:無脊椎動物>脊椎動物,甲殼類>魚類。脊椎動物和無脊椎動物敏感度相差顯著,脊椎動物的HC5是無脊椎動物的190多倍;魚類和甲殼類相比,相差近20倍。脊椎動物與無脊椎動物相比,所處營養(yǎng)等級更高,生理構(gòu)造更加復(fù)雜,體內(nèi)解毒機(jī)制更加完善,對銅的耐受性更強(qiáng),敏感性因而較低;脊椎動物在物種個(gè)體形體上也較無脊椎動物大,能夠積累更多的銅離子,也能使得其敏感性較低。另外,銅對不同物種的致毒機(jī)理也是不同,銅所作用的酶不同、不同物種對銅的吸收速率不同以及物種所處生命階段不同等也影響敏感性[59]。

    3.3 國內(nèi)外銅水生生物水質(zhì)基準(zhǔn)比較

    查閱文獻(xiàn)獲得部分國內(nèi)外淡水水體中銅的水生生物水質(zhì)基準(zhǔn)[1,60-63],具體見表3,將本研究所得的太湖銅水質(zhì)基準(zhǔn)與國內(nèi)外相關(guān)基準(zhǔn)進(jìn)行比較分析。孔祥臻等[60]研究得到我國重金屬銅對淡水水生物的急性HC5為4.5 μg·L-1,吳豐昌等[1]根據(jù)物種敏感度分布法得到我國銅的淡水生物水質(zhì)基準(zhǔn)急性基準(zhǔn)為30 μg·L-1,慢性基準(zhǔn)為9.44 μg·L-1,本研究得到太湖銅對全部物種的急性HC5為14.57 μg·L-1,慢性HC5為3.26 μg·L-1,高于孔祥臻等[60]的研究結(jié)果,低于吳豐昌等[1]的研究結(jié)果,由于太湖是我國淡水湖泊典型代表湖泊,太湖水體中的物種和我國淡水湖泊中的物種極為相似,因此得到的結(jié)果和我國淡水水質(zhì)基準(zhǔn)相差不顯著;由于推導(dǎo)過程使用的方法不同,我國淡水湖泊中的物種較太湖豐富,推導(dǎo)水質(zhì)基準(zhǔn)所選擇的物種存在差異,得到的結(jié)果不同。研究結(jié)果高于澳大利亞基于物種敏感度分布法得到銅的淡水水質(zhì)基準(zhǔn)1.4 μg·L-1[62],可能是水體中物種存在差異以及氣候水質(zhì)不同,另外本研究使用的概率物種敏感度分布法與澳大利亞推導(dǎo)使用的物種敏感度分布法也有所區(qū)別。吳豐昌等[1]基于毒性百分?jǐn)?shù)排序法獲得我國淡水水體中銅的基準(zhǔn)最大濃度和基準(zhǔn)連續(xù)濃度分別為9.10 μg·L-1和5.63 μg·L-1,與本研究的結(jié)果較為接近;美國毒性百分?jǐn)?shù)排序法得到的基準(zhǔn)最大濃度和基準(zhǔn)連續(xù)濃度分別為13.0 μg·L-1和9.00 μg·L-1[63],該結(jié)果與本研究結(jié)果也比較接近。本研究所得太湖銅水質(zhì)基準(zhǔn)與國內(nèi)外同類研究結(jié)果相差不大,沒有數(shù)量級的差異,產(chǎn)生差異的原因是多方面的,例如水質(zhì)基準(zhǔn)明顯的區(qū)域性,推導(dǎo)基準(zhǔn)的方法的不同等。

    圖4 銅的魚類、甲殼類、脊椎動物和無脊椎動物的概率物種敏感度分布曲線注:(a)魚類;(b)甲殼類;(c)脊椎動物;(d)無脊椎動物。Fig. 4 Probabilistic species sensitivity distribution curves for fishes, crustaceans, vertebrates and invertebrates in response to copperNote: (a) Fish; (b) Crustacea; (c) Vertebrate; (d) Invertebrate.

    所獲得物種毒性值中只有2個(gè)物種的急性毒性值小于基準(zhǔn)值(蚤狀溞1.3 μg·L-1和顫蚓0.16 μg·L-1),從而所得基準(zhǔn)值能保護(hù)太湖水體中95%以上的物種免受銅危害。余海靜等[13]采樣分析得到太湖水體總Cu為2.22 μg·L-1,總體來看,不會對水生生物造成短期和長期危害,但是北部部分湖區(qū)濃度高于4 μg·L-1,另外,北部湖區(qū)分布著太湖的主要入湖河流,接納無錫、宜興等城市的工業(yè)和生活污水,這些廢水中含有大量重金屬,再加上重金屬難以排出體外,有著生物累積效應(yīng)和放大效應(yīng),長期來看會對太湖水生生物產(chǎn)生危害,建議有關(guān)部門應(yīng)給予重視并采取相應(yīng)措施。

    綜上所述,可以得出以下結(jié)論:1)應(yīng)用概率物種敏感度分布法得到太湖銅的急性水質(zhì)基準(zhǔn)值為14.57 μg·L-1,慢性水質(zhì)基準(zhǔn)值為3.26 μg·L-1,并且與國內(nèi)外類似研究結(jié)果相差不大,所得到的基準(zhǔn)能夠保護(hù)太湖95%以上的物種。根據(jù)所得到的基準(zhǔn)值,總體來看,太湖水體中Cu不會對水生生物造成危害,但是北部部分湖區(qū)水生生物會受到Cu的長期危害,有關(guān)部門需給予重視并采取相應(yīng)措施。

    2)不同物種類別物種敏感性存在一定的差異,脊椎動物、無脊椎動物、魚類和甲殼類的HC5分別為34.75 μg·L-1、0.18 μg·L-1、24.62 μg·L-1和1.30 μg·L-1,敏感性大小順序?yàn)椋簾o脊椎動物>脊椎動物,甲殼類>魚類。

    表3 太湖銅水質(zhì)基準(zhǔn)與國內(nèi)外相關(guān)基準(zhǔn)比較

    注:STHC5表示短期危害濃度;LTHC5表示長期危害濃度;CMC表示基準(zhǔn)最大濃度;CCC表示基準(zhǔn)連續(xù)濃度。

    Note: STHC5stands for short term hazardous concentration; LTHC5stands for long term hazardous concentration; CMC stands for criteria maximum concentration; CCC stands for criteria continuous concentration.

    3)將概率物種敏感度分布法與傳統(tǒng)的物種敏感度分布法進(jìn)行比較分析,概率物種敏感度分布法有諸多優(yōu)越性:PSSD考慮了水生生物對銅多種不同響應(yīng),能對水生生物多方面保護(hù);曲線擬合效果更好;SSD擬合曲線的選擇有很大的主觀性,選不同曲線得到的HC5有一定差異,PSSD避免了這點(diǎn),且得到的HC5更精確;PSSD受數(shù)據(jù)量大小影響相對更小,因而結(jié)果也更為穩(wěn)定。

    [1] 吳豐昌, 馮承蓮, 曹宇靜, 等. 我國銅的淡水生物水質(zhì)基準(zhǔn)研究[J]. 生態(tài)毒理學(xué)報(bào), 2011, 6(6): 617-628

    Wu F C, Feng C L, Cao Y J, et al. Aquatic life ambient freshwater quality criteria for copper in China [J]. Asian Journal of Ecotoxicology, 2011, 6(6): 617-628 (in Chinese)

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

    Feng C L, Wu F C, Zhao X L, et al. Water quality criteria research and progress [J]. Science China: Earth Science, 2012, 42(5): 646-656 (in Chinese)

    [3] 王俊能, 胡習(xí)邦, 李玉炫, 等. 物種敏感性分布在鉈污染淡水生物生態(tài)風(fēng)險(xiǎn)評估中的應(yīng)用[J]. 安全與環(huán)境學(xué)報(bào), 2012, 12(6): 134-139

    Wang J N, Hu X B, Li Y X, et al. Assessing ecological risk of sudden thallium pollution accident to freshwater organisms by species sensitivity distributions [J]. Journal of Safety and Environment, 2012, 12(6): 134-139 (in Chinese)

    [4] 李玉爽, 吳豐昌, 崔曉勇, 等. 中國苯的淡水水質(zhì)基準(zhǔn)研究[J]. 生態(tài)學(xué)雜志, 2012, 31(4): 908-915

    Li Y S, Wu F C, Cui X Y, et al. Freshwater quality criteria for benzene in China [J]. Chinese Journal of Ecology, 2012, 31(4): 908-915 (in Chinese)

    [5] 張瑞卿, 吳豐昌, 李會仙, 等. 應(yīng)用物種敏感度分布法研究中國無機(jī)汞的水生生物水質(zhì)基準(zhǔn)[J]. 環(huán)境科學(xué)學(xué)報(bào), 2012, 32(2): 440-449

    Zhang R Q, Wu F C, Li H X, et al. Deriving aquatic water quality criteria for inorganic mercury in China by species sensitivity distributions [J]. Acta Scientiae Circumstantiae, 2012, 32(2): 440-449 (in Chinese)

    [6] Wu F C, Meng W, Zhao X L, et al. China embarking on development of its own national water quality criteria system [J]. Environmental Science & Technology, 2010, 44(21): 7992-7993

    [7] 吳豐昌, 孟偉, 宋永會, 等. 中國湖泊水環(huán)境基準(zhǔn)的研究進(jìn)展[J]. 環(huán)境科學(xué)學(xué)報(bào), 2008, 28(12): 2385-2393

    Wu F C, Meng W, Song Y H, et al. Research progress in lake water quality criteria in China [J]. Acta Scientiae Circumstantiae, 2008, 28(12): 2385-2393 (in Chinese)

    [8] 吳豐昌, 馮承蓮, 曹宇靜, 等. 鋅對淡水生物的毒性特征與水質(zhì)基準(zhǔn)的研究[J]. 生態(tài)毒理學(xué)報(bào), 2011, 6(4): 367-382

    Wu F C, Feng C L, Cao Y J, et al. Toxicity characteristic of zinc to freshwater biota and its water quality criteria [J]. Asian Journal of Ecotoxicology, 2011, 6(4): 367-382 (in Chinese)

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

    Wu F C, Feng C L, Zhang R Q, et al. Derivation of water quality criteria for representative water-body pollutants in China [J]. Science China: Earth Science, 2012, 42(5): 665-672 (in Chinese)

    [10] 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

    [11] 秦延文, 張雷, 鄭丙輝, 等. 太湖表層沉積物重金屬賦存形態(tài)分析及污染特征[J]. 環(huán)境科學(xué), 2012, 33(12): 4291-4299

    Qin Y W, Zhang L, Zheng B H, et al. Speciation and pollution characteristics of heavy metals in the sediment of Taihu Lake [J]. Environmental Science, 2012, 33(12): 4291-4299 (in Chinese)

    [12] 戴秀麗, 孫成. 太湖沉積物中重金屬污染狀況及分布特征探討[J]. 上海環(huán)境科學(xué), 2001, 20(2): 71-74

    Dai X L, Sun C. The characteristics of heavy metals distribution and pollution in sediment from Lake Taihu [J]. Shanghai Environmental Sciences, 2001, 20(2): 71-74 (in Chinese)

    [13] 余海靜, 張深, 鄒國防, 等. 生物配體模型預(yù)測太湖水體中Cu的形態(tài)分布和生物有效性[J]. 環(huán)境化學(xué), 2014, 33(7): 1107-1114

    Yu H J, Zhang S, Zou G F, et al. Predicting copper speciation and bio-availability in Taihu Lake using Biotic Ligand Model (BLM) [J]. Environmental Chemistry, 2014, 33(7): 1107-1114 (in Chinese)

    [14] United States Environmental Protection Agency (US EPA). Understanding and Accounting for Method Variability in Whole Effluent Toxicity Applications under the National Pollutant Discharge Elimination System Program [R]. Washington DC: US EPA, 2000

    [15] Fadri G, Bernd N. A probabilistic method for species sensitivity distributions taking into account the inherent uncertainty and variability of effects to estimate environmental risk [J]. Integrated Environmental Assessment and Management, 2012, 9(1): 79-86

    [16] Fadri G, Roland W S, Bernd N. Probabilistic material flow modeling for assessing the environmental exposure to compounds: Methodology and an application to engineered nano-TiO2particles [J]. Environmental Modelling & Software, 2010, 25: 320-332

    [17] Hayashi T I, Kashiwagi N. A Bayesian approach to probabilistic ecological risk assessment: Risk comparison of nine toxic substances in Tokyo surface waters [J]. Environmental Science and Pollution Research, 2011, 18: 365-375

    [18] United States Environmental Protection Agency (US EPA). Guidelines for Deriving Numerical National Water Quality Criteria for the Protection of Aquatic Organisms and Their Uses [R]. (PB 85-227049). Washington DC: US EPA, NTIS, 1985

    [19] Klimisch H J, Andreae M, Tillmann U. A systematic approach for evaluating the quality of experimental toxicological and ecotoxicological data [J]. Regulatory Toxicology and Pharmacology, 1997, 25(1): 1-5

    [20] Organization for Economic Co-operation and Development. Report of the OECD Workshop on Extrapolation of Laboratory Aquatic Toxicity Data to the Real Environment [R]. Paris: OECD, 1992

    [21] 蘇海磊, 吳豐昌, 李會仙, 等. 太湖生物區(qū)系研究及與北美五大湖的比較[J]. 環(huán)境科學(xué)研究, 2011, 24(12): 1346-1354

    Su H L, Wu F C, Li H C, et al. Aquatic biota of Taihu Lake and comparison with those of the North American great lakes [J]. Research of Environmental Sciences, 2011, 24(12): 1346-1354 (in Chinese)

    [22] 蘇海磊, 吳豐昌, 李會仙, 等. 我國水生生物水質(zhì)基準(zhǔn)推導(dǎo)的物種選擇[J]. 環(huán)境科學(xué)學(xué)報(bào), 2012, 25(5): 506-511

    Su H L, Wu F C, Li H C, et al. Species selection for deriving water quality criteria for protection of aquatic organisms in China [J]. Research of Environmental Sciences, 2012, 25(5): 506-511 (in Chinese)

    [23] United States Environmental Protection Agency. ECOTOX Database [DB/OL]. http://cfpub.epa.gov/ecotox/

    [24] 劉偉, 吳孝兵, 趙娟. 重金屬Cu2+對錦鯽和日本沼蝦的急性毒性研究[J]. 資源開發(fā)與市場, 2008, 24(10): 868-870

    Liu W, Wu X B, Zhao J. Study on acute toxicity of Cu2+to crucian carp and Japan marsh shrimp [J]. Resource Development and Market, 2008, 24(10): 868-870 (in Chinese)

    [25] 劉曉旭, 施蔡雷, 曹慧, 等. 銅和阿特拉津?qū)ΞT江彩鯉的聯(lián)合毒性研究[J]. 杭州師范大學(xué)學(xué)報(bào), 2010, 9(4): 263-267

    Liu X X, Shi C L, Cao H, et al. Joint toxicity of Cu2+and atrazine on C. carpio var. color [J]. Journal of Hangzhou Normal University, 2010, 9(4): 263-267 (in Chinese)

    [26] 楊麗華, 方展強(qiáng), 鄭文彪. 重金屬對鯽魚的生物毒性及安全濃度評價(jià)[J]. 華南師范大學(xué)學(xué)報(bào), 2003(2): 101-106

    Yang L H, Fang Z Q, Zheng W B. Safety assessment and acute toxicity of heavy metals to crucian Carassius auratus [J]. Journal of South China Normal University, 2003(2): 101-106 (in Chinese)

    [27] 劉大勝, 孔強(qiáng), 王大榜, 等. 金屬離子魚類急性研究及對環(huán)境標(biāo)準(zhǔn)修訂的啟示意義[J]. 環(huán)境科學(xué)與管理, 2010, 35(4): 38-42

    Liu D S, Kong Q, Wang D B, et al. Fish acute toxicity of metal ions and its inspiration for environmental standards revision [J]. Environmental Science and Management, 2010, 35(4): 38-42 (in Chinese)

    [28] 陳萬光, 屈菊平, 鄧平平, 等. Cu2+,Pb2+對高體鳑鲏幼魚的急性毒性試驗(yàn)[J]. 江蘇農(nóng)業(yè)科學(xué), 2010(4): 243-244

    Chen W G, Qu J P, Deng P P, et al. Acute toxicity test of Cu2+and Pb2+on juvenile Rhodeus ocellatus [J]. Journal of Jiangsu Agriculture, 2010(4): 243-244 (in Chinese)

    [29] 楊建華, 宋維彥. 3種重金屬離子對中華鳑鲏魚的急性毒性及安全濃度研究[J]. 安徽農(nóng)業(yè)科學(xué), 2010, 38(23): 12481-12482, 12485

    Yang J H, Song W Y. Study on acute toxicity and safe concentration of 3 heavy metal ions to Rhodens sinensis gunther [J]. Journal of Anhui Agriculture, 2010, 38(23): 12481-12482, 12485 (in Chinese)

    [30] 何志強(qiáng), 鄭永華. 重金屬(銅、鋅和汞)對中華倒刺鲃生物毒性效應(yīng)的研究[D]. 重慶: 西南大學(xué), 2008

    He Z Q, Zheng Y H. Biotoxicity effect studies on heavy metals(copper, zinc and mercury) in Spinibarbus sinensis [D]. Chongqing: Southwest University, 2008 (in Chinese)

    [31] 王偉莉, 焦聰穎, 閆振廣, 等. 水體硬度對銅和鎘生物毒性的影響[J]. 環(huán)境工程技術(shù)學(xué)報(bào), 2013, 3(3): 272-278

    Wang W L, Jiao C Y, Yan Z G, et al. Effects of water hardness on ecotoxicity of cadmium and copper to aquatic organisms [J]. Journal of Environmental Engineering Technology, 2012, 3(3): 272-278 (in Chinese)

    [32] 程霄玲, 鄭永華, 唐洪玉, 等. Cu2+、Zn2+、Cd2+對厚頜魴幼魚的聯(lián)合致毒效應(yīng)研究[J]. 淡水漁業(yè), 2009, 39(2): 54-59

    Chen X L, Zheng Y H, Tang H Y, et al. Study on the joint toxicity of Cu2+, Zn2+, Cd2+to the juvenile of Megalobrama pellegrini [J]. Freshwater Fisheries, 2009, 39(2): 54-59 (in Chinese)

    [33] 葉素蘭, 余治平. Cu2+、Pb2+、Cd2+、Cr6+對鳙胚胎和仔魚的急性致毒效應(yīng)[J]. 水產(chǎn)科學(xué), 2009, 28(5): 263-267

    Ye S L, Yu Z P. Acute toxicity of Cu2+, Pb2+, Cd2+and Cr6+to embryos and larvae of bighead carp Aristichthys nobilis [J]. Fisheries Science, 2009, 28(5): 263-267 (in Chinese)

    [34] 張穎. 重金屬對蛇鮈的急性毒性研究[J]. 西昌學(xué)院學(xué)報(bào), 2007, 21(1): 20-23

    Zhang Y. Acute toxicity of heavy metals to Longnose gudgeon [J]. Journal of Xichang College, 2007, 21(1): 20-23 (in Chinese)

    [35] 朱國霞, 季延濱, 孫學(xué)亮. Cu2+, Fe2+及銅鐵合劑對血鸚鵡幼魚的急性毒性試驗(yàn)[J]. 中國水產(chǎn), 2013(8): 56-58

    Zhu G X, Ji Y B, Sun X L. Acute toxicity test of Cu2+, Fe2+and copper iron micture on juvenile blood parrotfish [J]. Journal of Fishery of China, 2013(8): 56-58 (in Chinese)

    [36] 姚志峰, 章龍珍, 莊平, 等. 銅對中華鱘幼魚的急性毒性及對肝臟抗氧化酶活性的影響[J]. 中國水產(chǎn)科學(xué), 2010, 17(4): 731-738

    Yao Z F, Zhang L Z, Zhuang P, et al. Effects of antioxidant enzyme in liver and acute toxicity of Cu2+on juvenile Chinese sturgeon [J]. Journal of Fishery Sciences of China, 2010, 17(4): 731-738 (in Chinese)

    [37] 楊帆, 鄒容. 3種重金屬對瓦氏黃顙魚的急性毒性試驗(yàn)[J]. 水產(chǎn)科學(xué), 2008, 27(7): 368-369

    Yang F, Zou R. Acute toxicity of three heavy metals to yellow catfish Pelteobagrus vachelli [J]. Fisheries Science, 2008, 27(7): 368-369 (in Chinese)

    [38] 朱旭, 陳雷, 于俊海, 等. Cu2+對泥鰍急性毒性試驗(yàn)[J]. 遼寧化工, 2012, 41(9): 880-881

    Zhu X, Chen L, Yu J H, et al. Acute toxicity test of Cu2+for loach [J]. Liaoning Chemical Industry, 2012, 41(9): 880-881 (in Chinese)

    [39] 楊亞琴, 賈秀英. Cu2+、Zn2+和Cd2+對蟾蜍蝌蚪的聯(lián)合毒性[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2006, 12(3): 356-359

    Yang Y Q, Jia X Y. Joint toxicity of Cu2+、Zn2+and Cd2+to tadpole of Bufo bufo gargarizans [J]. Chinese Journal of Applied & Environmental Biology, 2006, 12(3): 356-359 (in Chinese)

    [40] 王愛民. 四種重金屬對綠蟾蜍蝌蚪的急性毒性研究[J]. 新疆大學(xué)學(xué)報(bào), 1990, 7(1): 60-64

    Wang A M. Acute toxicity of 4 heavy metals on the tadpoles of Bufo viridis laurenti [J]. Journal of Xinjiang University, 1990, 7(1): 60-64 (in Chinese)

    [41] 石戈, 王健鑫, 王日昕. Cu2+對中國林蛙蝌蚪生長發(fā)育的毒性效應(yīng)[J]. 吉林師范大學(xué)學(xué)報(bào), 2007(3): 71-73

    Shi G, Wang J X, Wang R X. Toxic effects of Cu2+on Rana chensinensis tadpole growth and development [J]. Journal of Jilin Normal University, 2007(3): 71-73 (in Chinese)

    [42] 姚丹, 萬琳燕, 耿寶榮, 等. Cu2+對日本林蛙蝌蚪的急性毒性研究[J]. 福建師范大學(xué)學(xué)報(bào), 2004, 20(4): 117-120

    Yao D, Wan L Y, Geng B R, et al. Acute toxicity of Cu2+to Rana japonica tadpoles [J]. Journal of Fujian Normal University, 2004, 20(4): 117-120 (in Chinese)

    [43] 黃帆, 郭正元, 徐珍, 等. 霸螨靈與Cu2+對蝌蚪的急性毒性試驗(yàn)研究[J]. 云南環(huán)境科學(xué), 2006, 25(2): 1-3

    Huang F, Guo Z Y, Xu Z, et al. Study on acute toxicity of Fenpyroximate and copper to tadpole [J]. Yunnan Environmental Science, 2006, 25(2): 1-3 (in Chinese)

    [44] 汪學(xué)英, 盧祥云, 李春梅, 等. 重金屬離子對黑斑蛙胚胎及蝌蚪的毒性影響[J]. 四川動物, 2001, 20(2): 59-61

    Wang X Y, Lu X Y, Li C M, et al. Toxicity of heavy metal irons to embryos and larvae of Rana nigromaculata [J]. Sichuan Journal of Zoology, 2001, 20(2): 59-61 (in Chinese)

    [45] 楊再福. 銅和汞對蝌蚪聯(lián)合毒性的影響[J]. 農(nóng)業(yè)環(huán)境保護(hù), 2001, 20(5): 370-371

    Yang Z F. Joint-toxicity of copper and hydrargyrum on tadpole [J]. Agro-environmental Protection, 2001, 20(5): 370-371 (in Chinese)

    [46] 王志錚, 呂敢堂, 施建軍, 等. 4種重金屬離子對中華鱉(Trionyx sinensis)稚鱉的急性毒性效應(yīng)[J]. 海洋與湖沼, 2009, 40(6): 745-752

    Wang Z Z, Lv G T, Shi J J, et al. Acute toxic effects of four heavy metals on Trionyx sinensis juveniles [J]. Oceanologia et Limnologia Sinica, 2009, 40(6): 745-752 (in Chinese)

    [47] 杞桑, 吳棉國, 張弦. 三種枝角類的實(shí)驗(yàn)生物學(xué)和急性毒性試驗(yàn)[J]. 暨南大學(xué)學(xué)報(bào), 1991, 12(3): 73-78

    Qi S, Wu J G, Zhang X. The biological characteristics of three cladocerans and acute toxic test [J]. Journal of Jinan University, 1991, 12(3): 73-78 (in Chinese)

    [48] 劉存歧, 安通偉, 張亞娟, 等. Cu2+對日本沼蝦的毒性研究[J]. 安徽農(nóng)業(yè)科學(xué), 2008, 36(28): 12285-12286

    Liu C Q, An T W, Zhang Y J, et al. Toxicity study of Cu2+on Macrobrachium nipponense [J]. Journal of Anhui Agriculture Science, 2008, 36(28): 12285-12286 (in Chinese)

    [49] 江敏, 臧維玲, 姚慶禎, 等. 四種重金屬對羅氏沼蝦仔蝦的毒性作用[J]. 上海水產(chǎn)大學(xué)學(xué)報(bào), 2002, 11(3): 203-207

    Jiang M, Zang W L, Yao Q Z, et al. The toxicity of four heavy metals on Macrobrachium rosenbergii postlarva [J]. Journal of Shanghai Fisheries University, 2002, 11(3): 203-207 (in Chinese)

    [50] 董學(xué)興, 呂林蘭, 王愛民, 等. Cu2+和Cd2+對克氏原螯蝦幼蝦的毒性效應(yīng)研究[J]. 水生態(tài)學(xué)雜志, 2010, 3(3): 90-93

    Dong X X, Lv L L, Wang A M, et al. Study on the acute toxicity of Cu2+and Cd2+acting on Procambarus clarkia juvenile [J]. Journal of Hydroecology, 2010, 3(3): 90-93 (in Chinese)

    [51] 米靜潔, 袁慧, 王蘭. 銅、鎘、鉻、鉛對河南華溪蟹的急性毒性作用[J]. 安徽農(nóng)業(yè)科學(xué), 2008, 36(17): 7273-7274, 7321

    Mi J J, Yuan H, Wang L. Study on the acute toxic effects of and on Sinopotamon henanense [J]. Journal of Anhui Agriculture Science, 2008, 36(17): 7273-7274, 7321 (in Chinese)

    [52] 趙含英, 楊家新, 陸正和, 等. Cu2+對萼花臂尾輪蟲的毒性影響[J]. 南京師大學(xué)報(bào), 2002, 25(4): 81-85, 90

    Zhao H Y, Yang J X, Lu Z H, et al. The toxic effect of Cu2+on rotifer Brachionus calyciflorus [J]. Journal of Nanjing Normal University, 2002, 25(4): 81-85, 90 (in Chinese)

    [53] 潘林. Pb2+和Cu2+對水螅的急性毒性作用以及對SOD活性的影響[J]. 高師理科學(xué)刊, 2010, 30(2): 72-75, 85

    Pan L. Effect of Pb2+and Cu2+on Hydra. sp. acute toxicity and SOD activity [J]. Journal of Science of Teachers’ College and University, 2010, 30(2): 72-75, 85 (in Chinese)

    [54] 陳德輝, 楊劍慧, 章宗涉, 等. 藻類的接觸致毒-恢復(fù)培養(yǎng)及其對化學(xué)物質(zhì)的安全評價(jià)[J]. 水生生物學(xué)報(bào), 1999, 23(5): 455-459

    Chen D H, Yang J H, Zhang Z S, et al. Algal toxicity test with exposure-recovery and its application to the safety assessmen of chemicals [J]. Acta Hydrobiologica Sinica, 1999, 23(5): 455-459 (in Chinese)

    [55] 劉益浩, 王超. 斜生柵藻與重金屬的相互作用研究[D]. 南京: 河海大學(xué), 2007

    Liu Y H, Wang C. Interactions between Scenedesmus Obliquus and heavy metals [D]. Nanjing: Hohai University, 2007 (in Chinese)

    [56] 閆海, 潘綱, 霍潤蘭. 銅、鋅和錳抑制月形藻生長的毒性效應(yīng)[J]. 環(huán)境科學(xué)學(xué)報(bào), 2001, 21(3): 328-332

    Yan H, Pan G, Huo R L. Oxic effects of copper, zinc and manganese on the inhibition of the growth of Closterium lunula [J]. Acta Scientiae Circumstantiae, 2001, 21(3): 328-332 (in Chinese)

    [57] 潘永全, 張?jiān)泼? 韓志剛, 等. 紅劍魚、孔雀魚、食蚊魚亞慢性毒性試驗(yàn)[J]. 四川動物, 2007, 26(3): 674-678

    Pan Y Q, Zhang Y M, Han Z G, et al. Subchronic toxicity test for Xiphophorus helleri, Lebistes reticulates and Gambusia affinis [J]. Sichuan Journal of Zoology, 2007, 26(3): 674-678 (in Chinese)

    [58] United States Environmental Protection Agency (US EPA). Great lakes water quality initiative criteria documents for the protection of aquatic life in ambient water [R]. (EPA-820-B-95-004). Washington DC: US EPA, Office of Water, 2009

    [59] 易秀. 銅對生物有機(jī)體的毒性及抗性機(jī)理[J]. 農(nóng)業(yè)環(huán)境保護(hù), 1997, 16(4): 187-189

    Yi X. Toxicity of copper on biological organisms and resistance mechanism of biological organisms [J]. Agro-environmental Protection, 1997, 16(4): 187-189 (in Chinese)

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

    Kong X Z, He W, Qin N, et al. Assessing acute ecological risks of heavy metals to freshwater organisms by species sensitivity distributions [J]. China Environmental Science, 2011, 31(9): 1555-1562 (in Chinese)

    [61] Canadian Council of Ministers of the Environment. Protocol for the derivation of water quality guidelines for the protection of aquatic life [R].Winnipeg: Canadian Council of Ministers of the Environment, 1999

    [62] Australia and New Zealand Environment and Conservation Council and Agriculture and Resource Management Council of Australia and New Zealand. Australia and New Zealand guidelines for fresh and marine water quality [R]. Canberra: ANZECC and ARMCANZ, 2000

    [63] United States Environmental Protection Agency (US EPA). National Recommended Water Quality Criteria [R]. Washington DC: Office of Water, Office of Science and Technology, 2009

    Deriving Aquatic Water Quality Criteria for Copper in Taihu Lake by Probabilistic Species Sensitivity Distributions

    Hou Jun1,2,*, Zhao Qianyuan1,2, Wang Chao1,2, Wang Peifang1,2, Miao Lingzhan1,2, Lv Bowen1,2

    1. Key Laboratory of Integrated Regulation and Resource Development on Shallow Lakes, Ministry of Education, Hohai University, Nanjing 210098, China 2. College of Environment, Hohai University, Nanjing 210098, China

    2 November 2014 accepted 24 November 2014

    Acute and chronic water quality criteria of copper for Taihu are derived to be 14.57 μg·L-1and 3.26 μg·L-1by the method of probabilistic species sensitivity distribution, respectively. Moreover, the results show that invertebrates are more sensitive than vertebrates, and crustaceans are more sensitive than fishes. Our method is superior to the traditional species sensitivity distribution. In our method, the toxicity effects are considered fully and reasonably, and data are fitted more efficiently. The amount of data exhibits little effect on the results. Furthermore, the method exhibits higher stability. The results can be served as reference of water quality standard of copper in Taihu. It can provide technical support for water environment management.

    copper; aquatic life; acute toxicity; chronic toxicity; water quality criteria; Taihu; probabilistic species sensitivity distribution

    國家杰出青年基金項(xiàng)目(51225901);國家自然科學(xué)基金(51479047;41430751;51209069);國家十二五水專項(xiàng)課題(2012ZX07101-008);江蘇省杰出青年基金項(xiàng)目(BK2012037);教育部創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃(IRT13061)

    侯俊(1979-),男,博士,副研究員,主要研究方向?yàn)樗Y源保護(hù)與生態(tài)修復(fù),E-mail: hhuhjyhj@126.com;

    *通訊作者(Corresponding author),E-mail: hhuhjyhj@126.com

    10.7524/AJE.1673-5897.20141102002

    2014-11-02 錄用日期:2014-11-24

    1673-5897(2015)1-191-13

    X171.5

    A

    侯俊, 趙芊淵, 王超, 等. 應(yīng)用概率物種敏感度分布法研究太湖銅水生生物水質(zhì)基準(zhǔn)[J]. 生態(tài)毒理學(xué)報(bào), 2015, 10(1): 191-203

    Hou J, Zhao Q Y, Wang C, et al. Deriving aquatic water quality criteria for copper in Taihu Lake by probabilistic species sensitivity distributions [J]. Asian Journal of Ecotoxicology, 2015, 10(1): 191-203 (in Chinese)

    猜你喜歡
    水生太湖敏感度
    全體外預(yù)應(yīng)力節(jié)段梁動力特性對于接縫的敏感度研究
    電視臺記者新聞敏感度培養(yǎng)策略
    新聞傳播(2018年10期)2018-08-16 02:10:16
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    在京韓國留學(xué)生跨文化敏感度實(shí)證研究
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    三氯生對4種水生生物的急性毒性研究
    太湖一角
    Diodes高性能汽車霍爾效應(yīng)閉鎖提供多種敏感度選擇
    依托科技創(chuàng)新 打造現(xiàn)代水生蔬菜產(chǎn)業(yè)
    長江蔬菜(2014年1期)2014-03-11 15:09:48
    欧美绝顶高潮抽搐喷水| 免费人成在线观看视频色| 成人av在线播放网站| 欧美成人性av电影在线观看| a级毛片免费高清观看在线播放| 丰满的人妻完整版| 三级国产精品欧美在线观看| 亚洲性久久影院| 日本一二三区视频观看| 99久久精品国产国产毛片| 最近最新免费中文字幕在线| 国产成人影院久久av| av天堂中文字幕网| 亚洲美女搞黄在线观看 | 日韩精品有码人妻一区| av在线蜜桃| 色5月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 久久久国产成人免费| 在线观看美女被高潮喷水网站| 精品午夜福利在线看| 国产精品亚洲美女久久久| 日韩欧美在线二视频| 欧美一区二区亚洲| 亚洲一级一片aⅴ在线观看| 麻豆成人av在线观看| 国产伦一二天堂av在线观看| 真实男女啪啪啪动态图| 亚洲欧美激情综合另类| 日本 欧美在线| 国产一区二区激情短视频| 黄片wwwwww| 国产伦精品一区二区三区视频9| 男女下面进入的视频免费午夜| 日本-黄色视频高清免费观看| 少妇熟女aⅴ在线视频| 人妻丰满熟妇av一区二区三区| 老司机午夜福利在线观看视频| 在线免费十八禁| 日韩欧美精品免费久久| 制服丝袜大香蕉在线| 国产精品人妻久久久影院| 国产午夜精品久久久久久一区二区三区 | 日本免费一区二区三区高清不卡| 国产伦在线观看视频一区| 国产精品一区二区免费欧美| 色尼玛亚洲综合影院| 99热这里只有精品一区| 亚洲久久久久久中文字幕| 人妻久久中文字幕网| x7x7x7水蜜桃| 非洲黑人性xxxx精品又粗又长| 99久久精品国产国产毛片| 久久久久国内视频| 99九九线精品视频在线观看视频| 亚洲av免费在线观看| eeuss影院久久| 亚洲精品亚洲一区二区| 不卡一级毛片| 全区人妻精品视频| 欧美日韩瑟瑟在线播放| 在线国产一区二区在线| 日日夜夜操网爽| 又爽又黄a免费视频| 日本在线视频免费播放| 国内精品久久久久久久电影| 国产美女午夜福利| 毛片女人毛片| 嫁个100分男人电影在线观看| 制服丝袜大香蕉在线| 亚洲最大成人av| 久久久久久国产a免费观看| 一进一出好大好爽视频| 免费看a级黄色片| 神马国产精品三级电影在线观看| 亚洲电影在线观看av| 又紧又爽又黄一区二区| 欧美人与善性xxx| 国产在线男女| 成人精品一区二区免费| 大又大粗又爽又黄少妇毛片口| 国产 一区 欧美 日韩| 色视频www国产| 高清毛片免费观看视频网站| 亚洲av一区综合| АⅤ资源中文在线天堂| 午夜福利欧美成人| 搡女人真爽免费视频火全软件 | 欧美高清性xxxxhd video| xxxwww97欧美| 最近在线观看免费完整版| 久久天躁狠狠躁夜夜2o2o| 极品教师在线免费播放| 亚洲国产欧洲综合997久久,| 88av欧美| 国产成年人精品一区二区| 亚洲av美国av| 色噜噜av男人的天堂激情| 色在线成人网| 亚洲av中文av极速乱 | 天美传媒精品一区二区| 欧美精品国产亚洲| 欧美黑人欧美精品刺激| 天堂动漫精品| 国产精品av视频在线免费观看| 五月伊人婷婷丁香| av在线观看视频网站免费| 日韩欧美国产一区二区入口| 亚洲中文字幕日韩| 成年版毛片免费区| 亚洲欧美日韩卡通动漫| 久久精品国产自在天天线| 中文字幕av成人在线电影| 三级国产精品欧美在线观看| 国产精品久久久久久久电影| 少妇高潮的动态图| .国产精品久久| 免费无遮挡裸体视频| 欧美日韩国产亚洲二区| 日韩欧美在线乱码| 小蜜桃在线观看免费完整版高清| 国产探花在线观看一区二区| 亚洲三级黄色毛片| 在线国产一区二区在线| 人妻制服诱惑在线中文字幕| 亚洲午夜理论影院| 亚洲国产日韩欧美精品在线观看| 熟妇人妻久久中文字幕3abv| 日韩国内少妇激情av| 亚洲国产日韩欧美精品在线观看| 91在线观看av| 成人永久免费在线观看视频| 欧美精品啪啪一区二区三区| 色综合色国产| 精品一区二区三区视频在线观看免费| 日本黄大片高清| 国产精品亚洲一级av第二区| 中文亚洲av片在线观看爽| 永久网站在线| 欧美三级亚洲精品| 精品一区二区三区人妻视频| 国产高清视频在线播放一区| 日韩一本色道免费dvd| 国产真实伦视频高清在线观看 | 精品久久久久久久末码| 久久6这里有精品| 99久久无色码亚洲精品果冻| 嫩草影院入口| 午夜老司机福利剧场| 国产精品久久久久久av不卡| 亚洲精华国产精华精| 校园人妻丝袜中文字幕| 久久天躁狠狠躁夜夜2o2o| 别揉我奶头~嗯~啊~动态视频| 日本黄色片子视频| 国产av一区在线观看免费| 国产综合懂色| 免费在线观看日本一区| 国产欧美日韩精品一区二区| 日韩欧美在线乱码| 国产精品国产高清国产av| 亚洲精品一卡2卡三卡4卡5卡| 久久香蕉精品热| 欧美三级亚洲精品| 国产一级毛片七仙女欲春2| 嫩草影视91久久| 精品久久久久久,| 热99在线观看视频| 国产成年人精品一区二区| 噜噜噜噜噜久久久久久91| av视频在线观看入口| 久久午夜福利片| 人妻少妇偷人精品九色| 一级黄片播放器| 亚洲欧美精品综合久久99| 亚洲男人的天堂狠狠| 乱码一卡2卡4卡精品| 看黄色毛片网站| 欧美日韩国产亚洲二区| 天天一区二区日本电影三级| 欧美不卡视频在线免费观看| 国内精品美女久久久久久| 亚洲精品亚洲一区二区| 一级a爱片免费观看的视频| 色吧在线观看| 成人特级av手机在线观看| 嫩草影院新地址| 久久久久九九精品影院| 在线免费十八禁| 国产精品久久久久久久久免| 九九爱精品视频在线观看| 乱码一卡2卡4卡精品| 真人做人爱边吃奶动态| 国产免费男女视频| 亚洲经典国产精华液单| 999久久久精品免费观看国产| 91久久精品电影网| 久久人人精品亚洲av| 人人妻人人澡欧美一区二区| 美女免费视频网站| 欧美bdsm另类| 美女高潮喷水抽搐中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 国产一区二区三区在线臀色熟女| 欧美一区二区亚洲| 男女视频在线观看网站免费| 一本久久中文字幕| 午夜精品久久久久久毛片777| 在线看三级毛片| 精品久久久久久成人av| 最近最新中文字幕大全电影3| 男女那种视频在线观看| 免费人成在线观看视频色| 天天一区二区日本电影三级| 亚洲三级黄色毛片| 少妇熟女aⅴ在线视频| 悠悠久久av| 久久人人精品亚洲av| 搡女人真爽免费视频火全软件 | 久久久久久久久久成人| 国产精品1区2区在线观看.| 三级国产精品欧美在线观看| 男人舔奶头视频| 久久精品夜夜夜夜夜久久蜜豆| 狂野欧美激情性xxxx在线观看| 亚洲精华国产精华精| 日本撒尿小便嘘嘘汇集6| 男女那种视频在线观看| 成年版毛片免费区| 男女啪啪激烈高潮av片| 男女下面进入的视频免费午夜| 深夜精品福利| 69人妻影院| 人妻少妇偷人精品九色| 最近最新中文字幕大全电影3| 色5月婷婷丁香| 国产午夜福利久久久久久| 欧美日本亚洲视频在线播放| 亚洲国产欧洲综合997久久,| 欧美色视频一区免费| 久久久精品大字幕| 日韩精品青青久久久久久| 男女做爰动态图高潮gif福利片| 亚洲无线在线观看| 97热精品久久久久久| 久久精品国产99精品国产亚洲性色| 十八禁国产超污无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 中文字幕av成人在线电影| 国产91精品成人一区二区三区| 欧美成人性av电影在线观看| 在线看三级毛片| 亚洲三级黄色毛片| 黄色配什么色好看| 免费看a级黄色片| 免费高清视频大片| 精品久久久久久久末码| 成年免费大片在线观看| 精品人妻一区二区三区麻豆 | 免费观看的影片在线观看| 久久久久性生活片| 国产色爽女视频免费观看| 欧美不卡视频在线免费观看| 亚洲专区中文字幕在线| 我的女老师完整版在线观看| 精品免费久久久久久久清纯| 日韩 亚洲 欧美在线| 蜜桃久久精品国产亚洲av| 欧美人与善性xxx| 久久热精品热| 一进一出好大好爽视频| 最近中文字幕高清免费大全6 | 久久精品国产鲁丝片午夜精品 | a级毛片a级免费在线| 日本 欧美在线| 蜜桃亚洲精品一区二区三区| 日韩在线高清观看一区二区三区 | avwww免费| 欧美性猛交黑人性爽| 男人的好看免费观看在线视频| 精品午夜福利视频在线观看一区| 国产午夜精品论理片| 欧美又色又爽又黄视频| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 亚洲美女搞黄在线观看 | 国产精品电影一区二区三区| 国语自产精品视频在线第100页| 亚洲真实伦在线观看| 国产黄a三级三级三级人| 最后的刺客免费高清国语| 欧美日韩亚洲国产一区二区在线观看| 91狼人影院| 亚洲av中文av极速乱 | 99热这里只有精品一区| 亚洲,欧美,日韩| 午夜福利在线在线| 变态另类丝袜制服| 在线观看午夜福利视频| 亚洲久久久久久中文字幕| 黄色一级大片看看| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 18禁黄网站禁片午夜丰满| 99热精品在线国产| 免费av不卡在线播放| 在线观看一区二区三区| 一个人免费在线观看电影| 夜夜夜夜夜久久久久| 真实男女啪啪啪动态图| 九色国产91popny在线| 日日干狠狠操夜夜爽| 精品人妻1区二区| 乱码一卡2卡4卡精品| 免费高清视频大片| 男女边吃奶边做爰视频| 黄色视频,在线免费观看| 97超级碰碰碰精品色视频在线观看| 国产激情偷乱视频一区二区| 欧美一区二区精品小视频在线| 国产精品,欧美在线| 亚洲精品久久国产高清桃花| 亚洲精品乱码久久久v下载方式| 在线观看舔阴道视频| 亚州av有码| 我要搜黄色片| 久久精品国产亚洲网站| 亚洲午夜理论影院| 午夜免费成人在线视频| 乱人视频在线观看| 精品无人区乱码1区二区| 欧美激情国产日韩精品一区| 国产精品,欧美在线| 国产日本99.免费观看| 三级男女做爰猛烈吃奶摸视频| 亚洲最大成人av| 又爽又黄无遮挡网站| 国产激情偷乱视频一区二区| 一区二区三区激情视频| 99久久精品热视频| 在线观看午夜福利视频| 国产综合懂色| 亚洲精品一区av在线观看| 在线观看美女被高潮喷水网站| 精品人妻熟女av久视频| 国产高清三级在线| 99久久久亚洲精品蜜臀av| 波多野结衣高清作品| 毛片一级片免费看久久久久 | 十八禁国产超污无遮挡网站| av.在线天堂| 久久这里只有精品中国| 一个人观看的视频www高清免费观看| 日本三级黄在线观看| 国产高清不卡午夜福利| 亚洲av免费在线观看| 在现免费观看毛片| 亚洲中文日韩欧美视频| 日韩欧美国产一区二区入口| 国产精品自产拍在线观看55亚洲| 女人被狂操c到高潮| 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 岛国在线免费视频观看| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看| 丰满人妻一区二区三区视频av| 俺也久久电影网| 国产激情偷乱视频一区二区| 人人妻,人人澡人人爽秒播| 国产精品女同一区二区软件 | 国产精品无大码| 丝袜美腿在线中文| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 欧美性猛交黑人性爽| 一级黄片播放器| 免费看光身美女| 国产精品av视频在线免费观看| 给我免费播放毛片高清在线观看| 亚洲三级黄色毛片| 欧美日韩黄片免| 欧美性猛交黑人性爽| 国产精品伦人一区二区| 婷婷精品国产亚洲av| 在线免费观看的www视频| 真人一进一出gif抽搐免费| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| 亚洲黑人精品在线| 国产精品一区二区免费欧美| 日韩人妻高清精品专区| 俄罗斯特黄特色一大片| 欧美高清性xxxxhd video| 日本黄色视频三级网站网址| 国产精品人妻久久久影院| 日韩 亚洲 欧美在线| 91狼人影院| 少妇被粗大猛烈的视频| 成年人黄色毛片网站| 国产精品人妻久久久久久| av国产免费在线观看| 免费观看人在逋| 亚洲avbb在线观看| 亚洲人成网站在线播放欧美日韩| 极品教师在线免费播放| 国产亚洲精品综合一区在线观看| 国模一区二区三区四区视频| 国产淫片久久久久久久久| 美女高潮的动态| 久久人妻av系列| 日韩精品中文字幕看吧| 搡老岳熟女国产| 少妇猛男粗大的猛烈进出视频 | 国产欧美日韩精品一区二区| 欧美日本视频| 999久久久精品免费观看国产| 3wmmmm亚洲av在线观看| 久久久久性生活片| 欧美又色又爽又黄视频| 国产中年淑女户外野战色| 18+在线观看网站| 亚洲经典国产精华液单| 精品一区二区三区人妻视频| 亚洲精品在线观看二区| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 成年女人毛片免费观看观看9| 国产精品人妻久久久久久| 91久久精品电影网| 亚洲av.av天堂| 制服丝袜大香蕉在线| 日韩欧美国产在线观看| 美女黄网站色视频| 国产色婷婷99| 亚洲真实伦在线观看| 久久中文看片网| 看黄色毛片网站| 美女大奶头视频| 免费看美女性在线毛片视频| 男插女下体视频免费在线播放| 老女人水多毛片| 韩国av一区二区三区四区| 国产黄色小视频在线观看| 一区福利在线观看| 精品国产三级普通话版| 国产精品自产拍在线观看55亚洲| 亚洲欧美日韩高清专用| 淫妇啪啪啪对白视频| 老熟妇仑乱视频hdxx| 国产一区二区在线av高清观看| 在线观看午夜福利视频| 老司机午夜福利在线观看视频| 国产精品一区二区免费欧美| 久久精品国产亚洲av涩爱 | 欧美日韩瑟瑟在线播放| 午夜久久久久精精品| 国内精品一区二区在线观看| 99久久精品热视频| 日韩欧美国产一区二区入口| 少妇丰满av| 国模一区二区三区四区视频| 啦啦啦观看免费观看视频高清| 我的老师免费观看完整版| 国产一区二区在线av高清观看| 高清毛片免费观看视频网站| 成人欧美大片| 国内精品美女久久久久久| 日韩欧美三级三区| 桃红色精品国产亚洲av| 婷婷精品国产亚洲av| 午夜免费激情av| 亚洲成人久久性| 亚洲人成伊人成综合网2020| 俺也久久电影网| 午夜久久久久精精品| 真人做人爱边吃奶动态| 亚洲三级黄色毛片| 伦理电影大哥的女人| 黄色女人牲交| 观看美女的网站| 午夜爱爱视频在线播放| 久久久久久久亚洲中文字幕| 美女黄网站色视频| 国产伦精品一区二区三区四那| 综合色av麻豆| 亚洲va日本ⅴa欧美va伊人久久| 国产精品av视频在线免费观看| 国产亚洲91精品色在线| 中亚洲国语对白在线视频| 欧美激情国产日韩精品一区| 国国产精品蜜臀av免费| 国产大屁股一区二区在线视频| 黄片wwwwww| 色视频www国产| 亚洲美女黄片视频| 天堂影院成人在线观看| 久久精品综合一区二区三区| 亚洲色图av天堂| 午夜精品在线福利| 精品久久久久久,| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 1024手机看黄色片| 在线天堂最新版资源| 中国美女看黄片| 亚洲色图av天堂| 男女视频在线观看网站免费| 不卡视频在线观看欧美| 亚洲性夜色夜夜综合| 黄色欧美视频在线观看| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 我的女老师完整版在线观看| 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 亚洲专区中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av涩爱 | 国产黄a三级三级三级人| 特级一级黄色大片| 日本一二三区视频观看| 亚洲真实伦在线观看| 国产在视频线在精品| 国产精品日韩av在线免费观看| 少妇熟女aⅴ在线视频| 亚洲四区av| 无遮挡黄片免费观看| 久久人妻av系列| 免费大片18禁| 免费观看在线日韩| 欧美中文日本在线观看视频| 亚洲精品色激情综合| 中国美白少妇内射xxxbb| 精品久久久久久久末码| 亚洲专区中文字幕在线| 国产精品永久免费网站| 国产精品爽爽va在线观看网站| 成年版毛片免费区| 乱系列少妇在线播放| avwww免费| 国产精品电影一区二区三区| 啪啪无遮挡十八禁网站| 男女啪啪激烈高潮av片| 日韩欧美在线二视频| 久久久久久久午夜电影| 国语自产精品视频在线第100页| 成人av在线播放网站| eeuss影院久久| 人妻丰满熟妇av一区二区三区| 欧美日韩乱码在线| 999久久久精品免费观看国产| 国产精品久久久久久久久免| 国产男人的电影天堂91| 国内精品一区二区在线观看| 可以在线观看的亚洲视频| 最新中文字幕久久久久| 天堂av国产一区二区熟女人妻| 国产精品久久久久久av不卡| 嫩草影视91久久| 久久久久久久久久久丰满 | 久久久精品大字幕| 校园人妻丝袜中文字幕| 两个人的视频大全免费| 亚洲美女搞黄在线观看 | 在线看三级毛片| 免费不卡的大黄色大毛片视频在线观看 | 成人二区视频| 国产av不卡久久| 免费无遮挡裸体视频| 18禁黄网站禁片午夜丰满| 亚洲美女视频黄频| 午夜免费男女啪啪视频观看 | 国产成人a区在线观看| 欧美成人性av电影在线观看| 国产高清视频在线播放一区| 色尼玛亚洲综合影院| 久久久精品欧美日韩精品| 国产熟女欧美一区二区| 国产精品伦人一区二区| 女生性感内裤真人,穿戴方法视频| 国产白丝娇喘喷水9色精品| 非洲黑人性xxxx精品又粗又长| 久久中文看片网| 日韩,欧美,国产一区二区三区 | 国产一区二区在线观看日韩| 国产在线男女| 美女cb高潮喷水在线观看| eeuss影院久久| 午夜日韩欧美国产| 亚洲va在线va天堂va国产| 国产午夜精品久久久久久一区二区三区 | 成熟少妇高潮喷水视频| 亚洲精品国产成人久久av| 动漫黄色视频在线观看| 中文字幕免费在线视频6| 亚洲最大成人av| 夜夜看夜夜爽夜夜摸| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲av香蕉五月| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩无卡精品| 国产欧美日韩精品一区二区| 日韩强制内射视频| 又爽又黄a免费视频| 久久人人爽人人爽人人片va| 在线播放无遮挡| 亚洲七黄色美女视频| 99久久九九国产精品国产免费| 亚洲黑人精品在线| 欧美区成人在线视频| 男人和女人高潮做爰伦理| 男人舔奶头视频|