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

    鎘的淡水水生生物水質(zhì)基準(zhǔn)研究

    2011-09-23 02:13:46吳豐昌曹宇靜李會(huì)仙張瑞卿馮承蓮閆振廣
    環(huán)境科學(xué)研究 2011年2期
    關(guān)鍵詞:基準(zhǔn)值水生敏感度

    吳豐昌,孟 偉*,曹宇靜,3,李會(huì)仙,張瑞卿,2,馮承蓮,閆振廣

    1.中國(guó)環(huán)境科學(xué)研究院,國(guó)家環(huán)境保護(hù)湖泊污染控制重點(diǎn)實(shí)驗(yàn)室,北京 100012

    2.中國(guó)科學(xué)院廣州地球化學(xué)研究所,廣東 廣州 510640

    3.中國(guó)鐵道科學(xué)研究院,節(jié)能環(huán)保勞衛(wèi)研究所,北京 100081

    鎘的淡水水生生物水質(zhì)基準(zhǔn)研究

    吳豐昌1,孟 偉1*,曹宇靜1,3,李會(huì)仙1,張瑞卿1,2,馮承蓮1,閆振廣1

    1.中國(guó)環(huán)境科學(xué)研究院,國(guó)家環(huán)境保護(hù)湖泊污染控制重點(diǎn)實(shí)驗(yàn)室,北京 100012

    2.中國(guó)科學(xué)院廣州地球化學(xué)研究所,廣東 廣州 510640

    3.中國(guó)鐵道科學(xué)研究院,節(jié)能環(huán)保勞衛(wèi)研究所,北京 100081

    鎘是一種有毒重金屬,具有高毒性、難降解和易殘留等特點(diǎn),會(huì)對(duì)水生生物及水生態(tài)系統(tǒng)產(chǎn)生有害影響.為有效控制鎘給水生生物帶來的不利影響,亟需開展鎘的水生生物基準(zhǔn)研究,為水質(zhì)標(biāo)準(zhǔn)的制訂提供依據(jù).以我國(guó)淡水生態(tài)系統(tǒng)及其生物區(qū)系為保護(hù)對(duì)象,結(jié)合大量國(guó)內(nèi)外文獻(xiàn)報(bào)道的鎘對(duì)我國(guó)淡水生物區(qū)系中代表物種的毒理學(xué)數(shù)據(jù),運(yùn)用評(píng)價(jià)因子法、毒性百分?jǐn)?shù)排序法和物種敏感度分布法等當(dāng)前國(guó)際上廣泛使用的基準(zhǔn)推導(dǎo)方法,研究我國(guó)淡水中鎘的基準(zhǔn)值及其推導(dǎo)過程.結(jié)果表明:評(píng)價(jià)因子法得出的基準(zhǔn)值為單值,其值為0.15μg/L;毒性百分?jǐn)?shù)排序法得出的基準(zhǔn)值包括基準(zhǔn)最大濃度和基準(zhǔn)連續(xù)濃度,二者分別為7.30和0.12μg/L;物種敏感度分布法得出的基準(zhǔn)值分為短期危險(xiǎn)濃度和長(zhǎng)期危險(xiǎn)濃度,二者分別為32.50和0.46μg/L.比較了3種方法的優(yōu)缺點(diǎn),以及與國(guó)內(nèi)外已有研究基準(zhǔn)值之間的差異及形成原因,分析了影響鎘的水生生物基準(zhǔn)的關(guān)鍵因素.

    鎘;水生生物基準(zhǔn);評(píng)價(jià)因子法;毒性百分?jǐn)?shù)排序法;物種敏感度分布法

    Abstract:Cadmium is a poisonous heavy metal which is toxic,hard to degrade and easy to reside.It can cause adverse effects on aquatic organisms and ecosystems.In order to control effectively the adverse effectswhich Cd might bring to aquatic life in Chinese freshwaters,it is urgent for China to derive regional aquatic life criteria for Cd,which could provide a basis for the establishment and revision of water quality standards.In this study,all available toxicity data of Cd to Chinese representative species in freshwater were collected in order to protect the freshwater ecosystem and biota system.Three widely used criteria derivation methods concerning the assessment factormethod,toxicity percentile rank method and species sensitivity distribution method were used to derive aquatic life criteria for Cd.Meanwhile,the criteria values for freshwater in China and the derivation process were studied and compared among the three methods. The results showed that for the assessment factormethod,the criteria of freshwater Cd was expressed by one value,which was 0.15 μg/L;for toxicity percentile rank method,the criteria included criteria maximum concentration and criteria continuous concentration,which were 7.30μg/L and 0.12 μg/L, respectively;for species sensitivity distributionmethod,the criteria of short term hazardous concentration and long term hazardous concentration were 32.50μg/L and 0.46 μg/L,respectively.Finally,this study analyzed the advantages and disadvantages of these three methods and compared the difference of the Cd criteria in this research and other existing reference values in China and abroad.In addition,the possible reasons which caused this difference and the key factors which influenced the aquatic life criteria were also discussed.

    Keywords:Cd;aquatic life criteria;assessment factor method;toxicity percentile rank method;species sensitivity distribution method

    水生生物基準(zhǔn)是保護(hù)水生生物水質(zhì)基準(zhǔn)的簡(jiǎn)稱,是指水環(huán)境中的污染物對(duì)水生生物及其使用功能不產(chǎn)生長(zhǎng)期和短期不良或有害效應(yīng)的最大允許濃度,它是水質(zhì)基準(zhǔn)的核心組成部分之一[1].水生生物基準(zhǔn)以棲息于水生態(tài)系統(tǒng)中的水生生物為保護(hù)目標(biāo),旨在維護(hù)水生態(tài)系統(tǒng)結(jié)構(gòu)和功能的完整性,以及生物多樣性,是國(guó)家制訂水質(zhì)標(biāo)準(zhǔn)、評(píng)價(jià)水質(zhì)和進(jìn)行水質(zhì)管理的科學(xué)依據(jù)[2-3].水生態(tài)系統(tǒng)具有一定的區(qū)域性和層次性,區(qū)域性水環(huán)境特征如 pH,硬度,堿度,溫度,溶解氧和有機(jī)質(zhì)等多種因素都會(huì)影響污染物在水環(huán)境中的生物有效性、生物吸收和富集過程[4-6],進(jìn)而影響到該物質(zhì)的水生生物基準(zhǔn).另外,不同國(guó)家/區(qū)域的水生生物區(qū)系存在一定的差異,由于物種敏感度差異性,同一污染物的基準(zhǔn)值也可能存在差異[7].目前我國(guó)的水質(zhì)標(biāo)準(zhǔn)主要是參考和借鑒國(guó)外基準(zhǔn)和標(biāo)準(zhǔn),但是國(guó)外的基準(zhǔn)是基于國(guó)外的生態(tài)系統(tǒng)特征和國(guó)情推導(dǎo)出來的,有的適用我國(guó),有的可能不適用.所以,有必要開展我國(guó)水質(zhì)基準(zhǔn)的研究工作,為水質(zhì)標(biāo)準(zhǔn)的制訂提供依據(jù).

    鎘是一種有毒重金屬元素,它具有高毒性、難降解、易殘留等特點(diǎn),會(huì)影響水生生物胚胎的發(fā)育、幼體的存活以及成體的繁殖[8-9].鑒于鎘的毒性,世界衛(wèi)生組織于1972年將鎘列為第3位優(yōu)先研究的食品污染物,聯(lián)合國(guó)環(huán)境規(guī)劃署于1974年將其定為重點(diǎn)污染物,美國(guó)國(guó)家環(huán)境保護(hù)局(US EPA)于1986年將鎘列為制定水質(zhì)基準(zhǔn)時(shí)的優(yōu)先控制污染物.目前,許多國(guó)家頒布了鎘的淡水水生生物基準(zhǔn),如美國(guó)已頒布了7次鎘的淡水水生生物基準(zhǔn),加拿大和澳大利亞也分別于1996年和2000年頒布了鎘的淡水水生生物基準(zhǔn).而我國(guó)水質(zhì)基準(zhǔn)的研究相對(duì)滯后.閆振廣等[10]參照 US EPA推薦的水生生物基準(zhǔn)推導(dǎo)方法——毒性百分?jǐn)?shù)排序法,利用國(guó)內(nèi)外報(bào)道的我國(guó)一些水生生物物種的毒理學(xué)數(shù)據(jù),研究了我國(guó)鎘的淡水水生生物基準(zhǔn)體系;曹宇靜等[11]同樣利用毒性百分?jǐn)?shù)排序法,主要采用我國(guó)學(xué)者報(bào)道的以我國(guó)淡水水生生物作為受試物種的鎘的毒理學(xué)數(shù)據(jù),探討了基于我國(guó)區(qū)域特點(diǎn)的鎘淡水水生生物基準(zhǔn).為了更加全面、系統(tǒng)地研究鎘的淡水水生生物基準(zhǔn)的理論和方法,選取我國(guó)淡水生物區(qū)系中代表性水生生物,調(diào)研了國(guó)內(nèi)外以這些物種作為受試物種的毒理學(xué)數(shù)據(jù),然后用國(guó)際上常用的3種方法分別研究了我國(guó)淡水中鎘的水生生物基準(zhǔn)值和推導(dǎo)方法,并探討這3種方法的優(yōu)缺點(diǎn)及其與國(guó)外基準(zhǔn)的差異性,以及影響水質(zhì)基準(zhǔn)值的關(guān)鍵因素.

    1 生物區(qū)系中代表物種的選擇

    我國(guó)淡水環(huán)境特別是湖泊種類繁多、分布廣泛.淡水中的水生生物受所處的地理位置、氣候條件、水文特點(diǎn)、理化性質(zhì)和人類活動(dòng)等多種因素綜合影響,具有明顯的時(shí)空差異性.據(jù)調(diào)查[12],我國(guó)常見的水生浮游植物約有200個(gè)屬,高等植物約有84種,浮游動(dòng)物約有1 000種(其中原生動(dòng)物360余種、輪蟲類250余種、枝角類130余種、橈足類200余種),底棲動(dòng)物主要有3個(gè)門(環(huán)節(jié)動(dòng)物門、軟體動(dòng)物門和節(jié)肢動(dòng)物門),魚類約有34科(其中鯉科在極大多數(shù)湖泊中占物種總數(shù)的50%以上).由于水生生物種類繁多,在選擇受試物種時(shí)不可能涵蓋所有生物,因此按照以下原則選擇受試物種:①要充分考慮物種的多樣性,尤其是急性毒性試驗(yàn),受試物種要盡可能涵蓋魚類、底棲類和浮游類等生物;②受試物種能反映我國(guó)的生物區(qū)系特征,以我國(guó)水生態(tài)系統(tǒng)中的代表性生物為優(yōu)選對(duì)象;③應(yīng)包含有重要經(jīng)濟(jì)價(jià)值或娛樂用途的物種.同時(shí),由于不同的推導(dǎo)方法對(duì)物種的要求有所不同,如評(píng)價(jià)因子法選擇受試物種時(shí)在遵循以上原則的基礎(chǔ)上,僅要求篩選出對(duì)某種污染物最敏感的物種即可[13];而毒性百分?jǐn)?shù)排序法在選擇受試物種時(shí)除要滿足上述原則外,必須至少涵蓋3門8科的生物[1].該研究根據(jù)推導(dǎo)方法對(duì)物種的數(shù)據(jù)要求,篩選出受試物種如表1所示.

    表1 推導(dǎo)我國(guó)水生生物基準(zhǔn)擬采用的受試物種Table 1 The tested species of our country's biota system used in deriving aquatic life criteria

    2 試驗(yàn)理化參數(shù)的確定

    為保證試驗(yàn)數(shù)據(jù)的可靠性和準(zhǔn)確性,在試驗(yàn)過程中應(yīng)嚴(yán)格控制以下理化參數(shù)[14]:①試驗(yàn)溫度.視受試物種而定,一般情況下應(yīng)采用受試生物的最適生長(zhǎng)溫度.②溶解氧濃度.應(yīng)維持在其飽和濃度的60% ~105%之間.③光周期.一般情況下,光暗周期比為16 h∶8 h.④pH以中性為宜,通常應(yīng)在7~8之間.⑤水質(zhì)參數(shù).試驗(yàn)所用稀釋水的各項(xiàng)參數(shù)須保持恒定,ρ(顆粒物)≤20 mg/L,ρ(總有機(jī)碳)≤2 mg/L,ρ(CODCr)≤5 mg/L,ρ(非離子氨)≤1μg/L,ρ(殘留氯)<3μg/L,ρ(總有機(jī)磷農(nóng)藥)≤50 ng/L,ρ(多氯聯(lián)苯)≤50 ng/L,ρ(有機(jī)氯)≤25 ng/L.試驗(yàn)理化參數(shù)應(yīng)視具體的污染物和受試物種而定,且在試驗(yàn)過程中應(yīng)測(cè)定試驗(yàn)用水的水質(zhì)參數(shù),如硬度、溶解氧、顆粒物濃度和總有機(jī)碳濃度等,以便在適當(dāng)?shù)臅r(shí)候?qū)鶞?zhǔn)值進(jìn)行校正.

    3 毒理數(shù)據(jù)的篩選

    3.1 選擇依據(jù)

    不同的試驗(yàn)類型對(duì)毒理數(shù)據(jù)的要求不同,應(yīng)分別予以考慮.急性毒性試驗(yàn)數(shù)據(jù)應(yīng)建立在能夠反映待測(cè)物質(zhì)對(duì)受試物種的急性嚴(yán)重不利影響的終點(diǎn)之上,其選擇依據(jù)[1]:①當(dāng)受試物種為水蚤或其他水蚤類動(dòng)物時(shí),應(yīng)使用齡期小于24 h的生物進(jìn)行試驗(yàn);當(dāng)受試物種為蚊類時(shí),應(yīng)使用其第2代或第3代幼蟲進(jìn)行試驗(yàn).試驗(yàn)結(jié)果應(yīng)以48 hEC50或LC50表示,如果受試生物在試驗(yàn)結(jié)束時(shí)沒有出現(xiàn)反?,F(xiàn)象,也可采用大于48 h的EC50或LC50;②當(dāng)受試物種為其他生物時(shí),應(yīng)以 96 hEC50或LC50表示.如果試驗(yàn)暴露時(shí)間在48~96 h之間,則采用試驗(yàn)結(jié)束時(shí)的LC50或EC50.

    慢性毒性試驗(yàn)數(shù)據(jù)應(yīng)建立在能夠反映待測(cè)物質(zhì)在不同暴露時(shí)間內(nèi)對(duì)受試物種產(chǎn)生的慢性不利影響的終點(diǎn)之上,其選擇依據(jù):①在整個(gè)生命周期試驗(yàn)中,當(dāng)受試物種為水蚤和糠蝦時(shí),應(yīng)使用齡期小于24 h的幼體進(jìn)行試驗(yàn),且暴露時(shí)間應(yīng)大于21 d;當(dāng)受試物種為魚類或其他物種時(shí),應(yīng)使用胚胎或齡期小于48 h的生物進(jìn)行試驗(yàn),試驗(yàn)至少要持續(xù)到孵化出下一代24 d后才可結(jié)束.②在部分生命周期試驗(yàn)中,當(dāng)受試物種為魚類時(shí),應(yīng)使用處于幼齡期的生物(至少要先于性腺發(fā)育前60 d)進(jìn)行試驗(yàn),且試驗(yàn)至少要持續(xù)到孵化出下一代24 d后才可結(jié)束;③在早期生命周期試驗(yàn)中,當(dāng)受試物種為魚類時(shí),應(yīng)使用剛完成受精的個(gè)體進(jìn)行試驗(yàn),且試驗(yàn)持續(xù)時(shí)間一般為28~32 d.慢性毒性試驗(yàn)結(jié)果應(yīng)以試驗(yàn)結(jié)束時(shí)的最大無觀察效應(yīng)濃度(NOEC)和最低觀察效應(yīng)濃度(LOEC)的幾何平均值表示.

    植物毒性試驗(yàn)數(shù)據(jù)應(yīng)建立在能夠反映待測(cè)物質(zhì)對(duì)受試物種總的不利影響的終點(diǎn)之上,其選擇依據(jù)為:①當(dāng)受試物種為藻類時(shí),試驗(yàn)結(jié)果應(yīng)以96 hLC50或EC50表示;②當(dāng)受試物種為水生維管束植物時(shí),試驗(yàn)結(jié)果應(yīng)用長(zhǎng)期的LC50或EC50表示.

    3.2 鎘的毒理數(shù)據(jù)的篩選

    根據(jù)我國(guó)淡水生物區(qū)系受試物種的選取原則,通過文獻(xiàn)調(diào)研,獲得了大量鎘對(duì)代表物種的毒理效應(yīng)數(shù)據(jù).該研究中鎘的毒性數(shù)據(jù)主要來源于中國(guó)知網(wǎng)(http://www.cnki.com/)和US EPA的毒理學(xué)數(shù)據(jù)庫(kù)(http://cfpub.epa.gov/ecotox/),數(shù)據(jù)收集截止到2009年12月.這些數(shù)據(jù)可以分為動(dòng)物急性毒性數(shù)據(jù)、動(dòng)物慢性毒性數(shù)據(jù)、植物毒性數(shù)據(jù)和生物富集數(shù)據(jù).其中,共有18個(gè)物種的急性毒性數(shù)據(jù)滿足推導(dǎo)基準(zhǔn)的要求,它們分屬于5門11科16屬(見表2).有15個(gè)物種的慢性毒性數(shù)據(jù)滿足推導(dǎo)基準(zhǔn)的要求,它們分屬于5門13科14屬(見表3).由于鎘對(duì)水生植物的毒性研究較少,國(guó)內(nèi)僅有少數(shù)學(xué)者研究了鎘對(duì)水生植物的生長(zhǎng)及生理毒性[54-56],且有的沒有建立鎘的劑量 -效應(yīng)關(guān)系,有的沒有給出鎘對(duì)水生植物的96 h或長(zhǎng)期的LC50或EC50值,因此,無法獲得滿足推導(dǎo)基準(zhǔn)要求的植物毒性數(shù)據(jù).有關(guān)鎘的生物富集研究主要集中于魚類、節(jié)肢動(dòng)物及軟體動(dòng)物的肌肉和內(nèi)臟組織[19,57-58].比較發(fā)現(xiàn),孟曉紅[57]測(cè)得的鯽魚肌肉組織對(duì)鎘的生物富集系數(shù)最小,為4;褚武英[58]測(cè)得的浦東圓田螺肌肉組織對(duì)鎘的生物富集系數(shù)最大,為8 125.關(guān)于鎘的最大允許組織濃度,國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局頒布了w(鎘)在不同水產(chǎn)品中的標(biāo)準(zhǔn)限值〔《農(nóng)產(chǎn)品質(zhì)量安全:無公害水產(chǎn)品安全要求》(GB 18406—2001)〕,其中甲殼類中為0.5 mg/kg,魚類中為0.1 mg/kg,貝類中為1.0 mg/kg,軟體類中為1.0 mg/kg[59].

    表2 鎘對(duì)淡水動(dòng)物的種平均急性毒性值和屬平均急性毒性值Table 2 The calculation of species mean acute value and genusmean acute value of cadmium to freshwater animal

    表3 鎘對(duì)淡水動(dòng)物的種平均慢性毒性值和屬平均慢性毒性值Table 3 The calculation of species mean chronic value and genusmean chronic value of cadm ium to freshwater animal

    4 水生生物基準(zhǔn)推導(dǎo)方法

    制訂水生生物基準(zhǔn)的核心是構(gòu)建推導(dǎo)水生生物基準(zhǔn)的理論與方法學(xué)體系,世界上許多國(guó)家和地區(qū)開展了水生生物基準(zhǔn)推導(dǎo)方法的研究工作,如美國(guó)、加拿大、歐盟和荷蘭等.目前國(guó)際上普遍采用的方法有評(píng)價(jià)因子法、物種敏感度分布法和毒性百分?jǐn)?shù)排序法3種,不同的推導(dǎo)方法原理不同,所需的基礎(chǔ)毒性數(shù)據(jù)量和基準(zhǔn)值表達(dá)方式也不相同.

    4.1 評(píng)價(jià)因子法

    評(píng)價(jià)因子法是世界上最早用于制訂水生生物基準(zhǔn)的一種方法.該方法基于化學(xué)物質(zhì)效應(yīng)評(píng)價(jià)的長(zhǎng)期經(jīng)驗(yàn),用敏感生物的毒性數(shù)據(jù)乘以相應(yīng)的評(píng)價(jià)因子或是帶入相應(yīng)的經(jīng)驗(yàn)公式中來定值,得出的基準(zhǔn)值為單值,并且以此作為污染物在任何情況都不得超過的濃度閾值[13].該方法簡(jiǎn)單易行,所需基礎(chǔ)數(shù)據(jù)較少,在毒性數(shù)據(jù)偏少的情況下,評(píng)價(jià)因子法因其通用性而被廣泛使用.評(píng)價(jià)因子法的有效性和評(píng)價(jià)因子的適用性在某種程度上主要依賴于敏感生物的毒性值.我國(guó)常見的基于急性毒性數(shù)據(jù)的經(jīng)驗(yàn)公式如下[60]:

    水質(zhì)基準(zhǔn)(安全濃度)=(24 hLC50×0.3)/

    水質(zhì)基準(zhǔn)(安全濃度)=(48 hLC50×0.3)/

    水質(zhì)基準(zhǔn)(安全濃度)=(96 hLC50)/AF(3)其中,式(3)應(yīng)用最為普遍,式中24 h,48 h和96 hLC50分別為24 h,48 h和96 h半致死濃度;AF為評(píng)價(jià)因子,無量綱.

    許多國(guó)家應(yīng)用評(píng)價(jià)因子法推導(dǎo)水生生物基準(zhǔn),只是對(duì)于評(píng)價(jià)因子的定值有所不同.加拿大規(guī)定[13]:基于敏感生物的急性毒性值推導(dǎo)水生生物基準(zhǔn)時(shí),持久性污染物的評(píng)價(jià)因子為100,非持久性污染物的評(píng)價(jià)因子為20;基于敏感生物的慢性毒性值推導(dǎo)水生生物基準(zhǔn)時(shí),所有污染物的評(píng)價(jià)因子均為10.我國(guó)的評(píng)價(jià)因子取10~100,其中對(duì)于易分解、低殘留的污染物,評(píng)價(jià)因子取10~20;對(duì)于穩(wěn)定的、易在水生生物體內(nèi)富集的污染物,評(píng)價(jià)因子取20~100[60].

    4.2 毒性百分?jǐn)?shù)排序法

    毒性百分?jǐn)?shù)排序法是隨著對(duì)污染物物理化學(xué)特性的了解和生態(tài)毒理學(xué)等學(xué)科的發(fā)展而逐步建立起來的,是 US EPA 1985年推導(dǎo)水生生物基準(zhǔn)的標(biāo)準(zhǔn)方法[1].該方法同時(shí)考慮了污染物的急性和慢性毒性效應(yīng),其得出的基準(zhǔn)值包括基準(zhǔn)最大濃度(CMC)和基準(zhǔn)連續(xù)濃度(CCC),其中CMC考慮的是污染物對(duì)水生動(dòng)物的急性毒性效應(yīng),是最終急性毒性值(FAV,mg/L)的1/2;CCC考慮的是污染物對(duì)水生動(dòng)物的慢性毒性效應(yīng)、對(duì)水生植物的毒性效應(yīng)以及污染物的生物富集效應(yīng),它取最終慢性毒性值(FCV,mg/L)、最終植物值(FPV,mg/L)和最終殘留值(FRV,mg/L)中的最小者.

    4.2.1 最終急性毒性值(FAV)

    FAV需要至少來自3門8科物種的急性毒性數(shù)據(jù),如果可以獲得足夠的數(shù)據(jù)且這些數(shù)據(jù)符合相關(guān)要求,則按以下步驟計(jì)算:

    第1步,計(jì)算各物種的種平均急性毒性值和屬平均急性毒性值.

    第2步,將屬平均急性毒性值從高到低排列,并且給其分配等級(jí)R,最小值的等級(jí)為1,最大值的等級(jí)為N(N為屬的個(gè)數(shù)).

    第3步,計(jì)算屬平均急性毒性值的累積概率(P):P=R/(N+1).

    第4步,選擇4個(gè)累積概率接近0.05的屬平均急性毒性值,用所選擇的屬平均急性毒性值和它們的累積概率計(jì)算最終毒性急性值.

    式中,s,L,A為計(jì)算過程中采用的符號(hào),沒有特殊含義;GMAV為屬平均急性毒性值,mg/L.

    4.2.2 最終慢性毒性值(FCV)

    如果可以獲得3門8科物種的慢性毒性數(shù)據(jù),最終慢性毒性值可以參考最終急性毒性值的方法計(jì)算.如果數(shù)據(jù)量不足,則最終慢性毒性值計(jì)算為:

    式中,F(xiàn)ACR為最終急性/慢性比率,無量綱.

    4.2.3 最終植物值(FPV)

    最終植物值等于水生植物毒性試驗(yàn)(通常是用藻類所做的96 h試驗(yàn)或者是用水生維管束植物所做的慢性試驗(yàn))結(jié)果中的最小值.由于植物對(duì)鎘的敏感性通常遠(yuǎn)遠(yuǎn)低于動(dòng)物,因此在很多情況下,可以不計(jì)算該值.

    4.2.4 最終殘留值(FRV)

    計(jì)算最終殘留值需要獲得該物質(zhì)的最大允許組織濃度(它是有關(guān)部門對(duì)魚油、魚類和貝類的可食用部分的管理水平)和生物富集系數(shù)(通常采用獲得數(shù)據(jù)中的最大值),然后按式(9)計(jì)算:

    式中,MPTC為最大允許組織濃度,mg/kg;BCF為生物富集系數(shù),L/kg.

    4.3 物種敏感度分布法

    物種敏感度分布理論最初是由歐美科學(xué)家于20世紀(jì)70年代提出的[61],主要用于生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),后來一些學(xué)者對(duì)該理論進(jìn)行了深入研究,為其在環(huán)境質(zhì)量基準(zhǔn)制訂中的應(yīng)用奠定了基礎(chǔ)[62-63].物種敏感度理論認(rèn)為:不同門類的生物由于生活史、生理構(gòu)造、行為特征和地理分布等不同而產(chǎn)生了差異性,其在毒理學(xué)上反映為不同的物種對(duì)污染物有不同的劑量 -效應(yīng)響應(yīng)關(guān)系,即不同的生物對(duì)同一污染物存在敏感性差異并遵循一定的概率分布模型.在獲得所需的毒性數(shù)據(jù)后,根據(jù)毒性數(shù)據(jù)的頻數(shù)分布擬合出某種概率分布函數(shù),即物種敏感度分布曲線,基準(zhǔn)值即為物種敏感度分布曲線上累積概率為X%時(shí)對(duì)應(yīng)的毒性值.根據(jù)不同的毒性數(shù)據(jù),該研究分別采用急性和慢性毒性數(shù)據(jù)擬合物種敏感度分布曲線,將得出的基準(zhǔn)值以短期危險(xiǎn)濃度(STHCX)和長(zhǎng)期危險(xiǎn)濃度(LTHCX)來表示.使用該方法推導(dǎo)水生生物基準(zhǔn)的一般步驟為:

    第1步,計(jì)算各物種的種平均急性/慢性毒性值.對(duì)于篩選出來的每個(gè)物種的毒性效應(yīng)終點(diǎn),種平均急性/慢性毒性值等于一個(gè)物種所有符合基準(zhǔn)推導(dǎo)要求的急性/慢性毒性值的幾何平均值.

    第2步,計(jì)算各物種平均急性/慢性毒性值的累積概率.將所有物種的種平均急性/慢性毒性值按從大到小的順序進(jìn)行排列,并且給其分配等級(jí)為I,最小的等級(jí)為1,最大的等級(jí)為N,計(jì)算公式如下:

    第3步,根據(jù)物種的毒性數(shù)據(jù)分布情況選擇適當(dāng)?shù)哪P蛿M合物種的急性或慢性敏感度分布曲線.推薦模型[63-66]有波爾Ⅲ模型(Burr TypeⅢ)、邏輯斯蒂累積密度模型(Logistic CDF)、對(duì)數(shù)正態(tài)累積密度模型(Lognormal CDF)、韋布爾累積密度模型(Weibull CDF)、蒙特卡羅模型(Monte Carlo)、高斯模型(Gaussian)、龔珀資模型(Gompertz)、指數(shù)增長(zhǎng)模型(Exponential Growth)和S型模型(Sigmoid)等.大量研究[67-69]表明,沒有任何一個(gè)模型對(duì)所有數(shù)據(jù)集都有很好的擬合效果.因此,應(yīng)根據(jù)區(qū)域生物區(qū)系毒性數(shù)據(jù)的實(shí)際分布情況確定擬合模型.

    第4步,利用選定模型擬合物種敏感度分布曲線.在確定毒性數(shù)據(jù)分布參數(shù)和選定模型后,利用物種的種平均急性/慢性毒性值或它們的對(duì)數(shù)值,及其相應(yīng)的累積概率,擬合物種的敏感度分布曲線,并且評(píng)價(jià)各模型的擬合度.

    第5步,確定短期危險(xiǎn)濃度(STHCX)和長(zhǎng)期危險(xiǎn)濃度(LTHCX).這2個(gè)值分別等于對(duì)應(yīng)的急性和慢性物種敏感度分布曲線上累積概率X%對(duì)應(yīng)的毒性值.一般歐美將X取為5,則STHC5和LTHC5指的是影響不超過5%的物種,即可以保護(hù)95%以上的物種時(shí)對(duì)應(yīng)的急性濃度和慢性濃度.最后通過與物種的種平均急性/慢性毒性值相比較,判斷各模型得出的基準(zhǔn)值的準(zhǔn)確度,最終選擇擬合度和準(zhǔn)確度均較高的值作為基準(zhǔn)值.

    5 結(jié)果

    5.1 評(píng)價(jià)因子法得出的鎘的水生生物基準(zhǔn)

    該研究擬采用水生生物的急性毒性數(shù)據(jù)計(jì)算鎘的基準(zhǔn)值.由表2可知,在所有受試生物中,羅氏沼蝦對(duì)鎘最為敏感,其96 hLC50為0.015 mg/L,根據(jù)評(píng)價(jià)因子選擇原則,鎘屬于穩(wěn)定的、易在水生生物體內(nèi)富集的污染物,且殘留時(shí)間較長(zhǎng)(其半衰期長(zhǎng)達(dá) 10~35 a)[8],對(duì)水生生物危害較大.因此,選擇100作為鎘的評(píng)價(jià)因子,將這2個(gè)值代入式(3),得出鎘的基準(zhǔn)值為0.15μg/L.

    5.2 毒性百分?jǐn)?shù)排序法得出的鎘的水生生物基準(zhǔn)

    5.2.1 最終急性毒性值(FAV)

    根據(jù)表2中的數(shù)據(jù)和最終急性毒性值的計(jì)算方法,將所需數(shù)據(jù)代入式(4)~(7)中,得出淡水中鎘的最終急性毒性值為14.60μg/L.

    5.2.2 最終慢性毒性值(FCV)

    根據(jù)表3中的數(shù)據(jù)和計(jì)算最終慢性毒性值的方法,可以采用和最終急性毒性值同樣的方法計(jì)算,將所需數(shù)據(jù)代入式(4)~(7)中,僅需將式中的屬平均急性毒性值(GMAV)換成屬平均慢性毒性值(GMCV)即可,得出鎘的最終慢性毒性值為0.21μg/L

    5.2.3 最終植物值(FPV)

    由于無法獲得滿足要求的植物毒性數(shù)據(jù),而且大量的研究表明,水生植物對(duì)污染物的敏感性遠(yuǎn)遠(yuǎn)低于水生動(dòng)物的敏感性,最終植物值對(duì)鎘的基準(zhǔn)值影響不大,因此該研究沒有計(jì)算鎘的最終植物值.

    5.2.4 最終殘留值(FRV)

    在所有可獲得的鎘的生物富集系數(shù)中,浦東圓田螺肌肉組織對(duì)鎘的生物富集系數(shù)最大,為8 125,由于圓田螺屬于軟體動(dòng)物,其最大允許組織濃度為1.0 mg/kg,將這2個(gè)值代入式(9),得出鎘的最終殘留值為0.12μg/L.

    毒性百分?jǐn)?shù)排序法得出的基準(zhǔn)值以基準(zhǔn)最大濃度(CMC)和基準(zhǔn)連續(xù)濃度(CCC)表示.根據(jù)5.2.1~5.2.4節(jié)的計(jì)算結(jié)果,得出鎘的二者分別為7.30和0.12μg/L.

    5.3 物種敏感度分布法得出的鎘的水生生物基準(zhǔn)

    5.3.1 短期危險(xiǎn)濃度

    表4 不同模型擬合鎘的急性物種敏感度分布曲線結(jié)果Table 4 The results of cadmium's acute species sensitivity distribution curves fitted by differentmodels

    根據(jù)物種敏感度分布的推導(dǎo)方法,首先采用Origin 8.0對(duì)表2中18個(gè)物種的種平均急性毒性值進(jìn)行對(duì)數(shù)正態(tài)分布檢驗(yàn),其中 S-W檢驗(yàn)所得的顯著性水平為P=0.203,K-S檢驗(yàn)所得的顯著性水平為P=0.927,均大于0.05,因此數(shù)據(jù)符合對(duì)數(shù)正態(tài)分布.利用 SigmaPlot 10.0擬合鎘的急性物種敏感度分布曲線,以各物種的種平均急性毒性值的對(duì)數(shù)值為橫坐標(biāo),以其對(duì)應(yīng)的累積概率為縱坐標(biāo),觀察這些坐標(biāo)點(diǎn)的分布趨勢(shì),采用了Sigmoid, Gaussian, Lorentzian, Gompertz 和Exponential Grow th等5個(gè)模型分別擬合了鎘的急性物種敏感度分布曲線,擬合結(jié)果見表4和圖1.

    圖1 采用不同模型擬合鎘的急性物種敏感度分布曲線Fig.1 The acute species sensitivity distribution curves of differentmodels

    從表4可知,各模型的擬合度相差不大,其中Sigmoid模型的擬合度最高,其決定系數(shù)(R2)為0.995,通過將它們和表2中各物種的種平均急性毒性值進(jìn)行比較,發(fā)現(xiàn)Sigmoid模型得出的基準(zhǔn)值準(zhǔn)確度較高.因此,該研究采用Sigmoid模型得出的32.50μg/L作為鎘的短期危險(xiǎn)濃度(STHC5).

    5.3.2 長(zhǎng)期危險(xiǎn)濃度

    長(zhǎng)期危險(xiǎn)濃度的推導(dǎo)方法跟短期危險(xiǎn)濃度類似,首先采用 Origin8.0對(duì)表3中15個(gè)物種的種平均慢性毒性值進(jìn)行對(duì)數(shù)正態(tài)分布檢驗(yàn),其中S-W檢驗(yàn)所得的顯著性水平為P=0.197,K-S檢驗(yàn)所得的顯著性水平為P=0.856,均大于0.05,符合對(duì)數(shù)正態(tài)分布.利用Sigma Plot 10.0擬合鎘的慢性物種敏感度分布曲線,以各物種的種平均慢性毒性值的對(duì)數(shù)值為橫坐標(biāo),以其對(duì)應(yīng)的累積概率為縱坐標(biāo),選擇適當(dāng)?shù)哪P瓦M(jìn)行擬合.采用Sigmoid,Gaussian和Exponential Growth 3個(gè)模型擬合了鎘的慢性物種敏感度分布曲線,擬合結(jié)果見表5和圖2.從表5可知,Sigmoid模型(4個(gè)參數(shù))的擬合度最高,其R2為0.954;各模型的擬合度相差不大,但是各模型得出的長(zhǎng)期危險(xiǎn)濃度(LTHC5)相差較大,將它們和表3中各物種的種平均慢性毒性值進(jìn)行比較發(fā)現(xiàn),采用Sigmoid模型(4個(gè)參數(shù))擬合所得的0.46μg/L準(zhǔn)確度較高.因此,該研究采用0.46μg/L作為鎘的長(zhǎng)期危險(xiǎn)濃度(LTHC5).

    表5 不同模型擬合鎘的慢性物種敏感度分布曲線結(jié)果Table 5 The results of cadmium's chronic species sensitivity distribution curves fitted by differentmodels

    圖2 不同模型擬合的鎘的慢性物種敏感度分布曲線Fig.2 The chronic species sensitivity distribution curves of differentmodels

    比較幾種模型的擬合結(jié)果發(fā)現(xiàn),用Sigmoid模型擬合物種的毒性數(shù)據(jù)的決定系數(shù)最大,最適合用于敏感度分布法進(jìn)行數(shù)據(jù)擬合,可以作為基準(zhǔn)推導(dǎo)的標(biāo)準(zhǔn)參考方法.

    6 討論

    該研究在調(diào)研我國(guó)水生生物區(qū)系中代表物種的基礎(chǔ)上,結(jié)合國(guó)內(nèi)外對(duì)以這些代表物種作為受試物種的毒理學(xué)試驗(yàn)數(shù)據(jù),采用目前國(guó)際上普遍使用的推導(dǎo)水生生物基準(zhǔn)的3種方法——評(píng)價(jià)因子法、毒性百分?jǐn)?shù)排序法和物種敏感度分布法,得出了以我國(guó)淡水生態(tài)系統(tǒng)為保護(hù)目標(biāo)的鎘的水生生物基準(zhǔn).需要說明的是,該研究得出的基準(zhǔn)值是針對(duì)總鎘的,在具體應(yīng)用時(shí),可根據(jù)實(shí)際水體的硬度將總鎘轉(zhuǎn)化為可溶性鎘,轉(zhuǎn)化系數(shù)可以參照 US EPA提出的轉(zhuǎn)換系數(shù)(CF)[70].該研究和其他國(guó)家基準(zhǔn)值比較見表6.

    表6 該研究和其他國(guó)家采用不同方法得出的基準(zhǔn)值對(duì)比Table 6 The criteria of this study and other countries used different derivingmethods

    從表6可知,3種方法得出的基準(zhǔn)值之間有差距,評(píng)價(jià)因子法得出的基準(zhǔn)值為單值,其余2種方法均為雙值,即鎘分別在短期暴露和長(zhǎng)期暴露下不得超過的濃度閾值,即短期濃度限值和長(zhǎng)期濃度限值.這3種方法均有優(yōu)缺點(diǎn).評(píng)價(jià)因子法的優(yōu)點(diǎn)在于所需基礎(chǔ)數(shù)據(jù)少、計(jì)算簡(jiǎn)單易行,當(dāng)生物的毒性數(shù)據(jù)量不能滿足其他推導(dǎo)方法時(shí),該方法是好的選擇.但存在一定的缺點(diǎn):①該方法屬于經(jīng)驗(yàn)法,不同國(guó)家對(duì)同一物質(zhì)評(píng)價(jià)因子的規(guī)定上往往相差較大,很難判斷哪個(gè)比較符合實(shí)際情況;②該方法得出的基準(zhǔn)值主要依賴于敏感生物的毒性值,如果最敏感生物的毒性值測(cè)量不準(zhǔn)確,所選生物并非最敏感生物時(shí),計(jì)算出的基準(zhǔn)值可能會(huì)有偏差;③該方法沒有考慮物種之間的相互關(guān)系以及污染物的生物富集效應(yīng),也缺乏數(shù)理統(tǒng)計(jì)理論的支持.上述缺點(diǎn)使得用該方法推導(dǎo)出的基準(zhǔn)值缺乏說服力,所以一般不宜采用.毒性百分?jǐn)?shù)排序法和物種敏感度分布法共同的優(yōu)點(diǎn)在于:①2種方法分別考慮了污染物的急性和慢性毒性效應(yīng),得出的基準(zhǔn)值為雙值,更加符合污染物的實(shí)際分布規(guī)律,能夠?yàn)樗锾峁┻m當(dāng)?shù)谋Wo(hù);②2種方法將測(cè)試生物按生物分類學(xué)進(jìn)行了歸類,考慮了同一物種內(nèi)不同生物或同一屬內(nèi)不同物種對(duì)污染物的劑量-效應(yīng)響應(yīng)關(guān)系的相似性,比較符合生物學(xué)規(guī)律;③2種方法都結(jié)合了數(shù)理統(tǒng)計(jì)理論推導(dǎo)基準(zhǔn)值,其中毒性百分?jǐn)?shù)排序法認(rèn)為敏感屬的毒性數(shù)據(jù)符合對(duì)數(shù) -三角分布,物種敏感度分布法認(rèn)為敏感物種的毒性數(shù)據(jù)符合對(duì)數(shù) -正態(tài)分布,具有統(tǒng)計(jì)學(xué)意義.然而,這2種方法也都有各自的不足之處,其中毒性百分?jǐn)?shù)排序法雖然計(jì)算了各物種和屬的毒性數(shù)據(jù),但最終用于計(jì)算基準(zhǔn)值的只是累積概率接近0.05的4個(gè)屬的毒性數(shù)據(jù),這使得該方法在一定程度上存在與評(píng)價(jià)因子法類似的缺點(diǎn),即最終基準(zhǔn)值很大程度上依賴于敏感物種的數(shù)據(jù).因此,利用該方法得出的基準(zhǔn)值是否具有代表性值得考慮.物種敏感度分布法需要進(jìn)一步考慮污染物在生物體內(nèi)的富集效應(yīng),特別是那些生物富集能力較強(qiáng)的污染物(如汞和鎘等);其次,使用的模型不同,擬合出的物種敏感度分布曲線不同,因而得出的基準(zhǔn)值可能也不同,所以不同區(qū)域基準(zhǔn)可能有特定的模型推導(dǎo)方法[71-72].因此,毒性百分?jǐn)?shù)排序法和物種敏感度分布法在推導(dǎo)基準(zhǔn)時(shí)明顯優(yōu)于評(píng)價(jià)因子法.該研究中采用這2種方法分別得出的鎘的短期濃度限值和長(zhǎng)期濃度限值在同一個(gè)數(shù)量級(jí)范圍內(nèi),但彼此之間差距較大,物種敏感度分布法得出的基準(zhǔn)值是毒性百分?jǐn)?shù)排序法的4倍左右.將這2種方法得出的基準(zhǔn)值和鎘的毒性數(shù)據(jù)相比較發(fā)現(xiàn),毒性百分?jǐn)?shù)排序法得出的基準(zhǔn)值可以保護(hù)更大范圍的水生生物,因此該研究認(rèn)為該方法得出的基準(zhǔn)值作為保護(hù)我國(guó)淡水水生生物的基準(zhǔn)值可能比較恰當(dāng),其基準(zhǔn)最大濃度(CMC)和基準(zhǔn)連續(xù)濃度(CCC)分別為7.30和0.12μg/L.

    由表6可知,該研究和其他國(guó)家的基準(zhǔn)值存在較大差異,這進(jìn)一步說明我國(guó)不應(yīng)該照搬國(guó)外的水質(zhì)基準(zhǔn),而應(yīng)該以我國(guó)生物區(qū)系中代表生物的毒理學(xué)數(shù)據(jù)為基礎(chǔ)[73].綜合分析,造成這種差異的主要原因有:①各國(guó)的生物區(qū)系不同.不同的生物區(qū)系中有不同的敏感物種,不同物種對(duì)同種污染物的敏感性不同.水質(zhì)基準(zhǔn)通常是基于保護(hù)95%物種為依據(jù)確定的,但是,5%敏感物種毒性數(shù)據(jù)對(duì)基準(zhǔn)值的最終確定具有決定作用.當(dāng)然敏感物種分布和種屬也是區(qū)域生物區(qū)系和水生態(tài)系統(tǒng)特征之一.推導(dǎo)我國(guó)鎘的水質(zhì)基準(zhǔn)時(shí)選用的最敏感的4個(gè)屬分別為石斑魚屬、水蚤屬、光殼蝦屬和沼蝦屬;而美國(guó)選用的最敏感的4個(gè)屬為鮭屬、鱸魚屬、紅點(diǎn)鮭屬和斑鱒魚屬(見表7).這4個(gè)鎘的敏感性物種屬毒性數(shù)據(jù)是決定水質(zhì)基準(zhǔn)的重要因素之一.②水質(zhì)參數(shù)不同,許多水質(zhì)參數(shù)如 pH,硬度,溫度,溶解氧和溶解有機(jī)質(zhì)等均會(huì)影響污染物的毒性,如鎘的毒性受硬度的影響較大[73-74].該研究未涉及這部分內(nèi)容,需要在以后的工作中進(jìn)一步考慮這些因素的影響.

    表7 中國(guó)和美國(guó)計(jì)算鎘的基準(zhǔn)時(shí)采用的最敏感屬對(duì)比Table 7 The comparison of the most sensitive genus which were used separately by China and America in calculating cadmium's criteria

    7 結(jié)論

    a.在調(diào)研我國(guó)水生生物區(qū)系中代表物種的基礎(chǔ)上,采用國(guó)際上普遍使用的3種方法分別推導(dǎo)了鎘的淡水水生生物基準(zhǔn).采用評(píng)價(jià)因子法得出的我國(guó)鎘的淡水水生生物基準(zhǔn)是單值,為 0.15 μg/L.采用毒性百分?jǐn)?shù)排序法得出基準(zhǔn)最大濃度為7.30μg/L,基準(zhǔn)連續(xù)濃度為0.12μg/L.采用物種敏感度分布法得出的鎘的短期危險(xiǎn)濃度為32.50μg/L,長(zhǎng)期危險(xiǎn)濃度為0.46μg/L.

    b.比較3種方法得出的基準(zhǔn)值發(fā)現(xiàn),毒性百分?jǐn)?shù)排序法得出的基準(zhǔn)值較其他2種方法更為恰當(dāng);與其他國(guó)家鎘的基準(zhǔn)值比較,發(fā)現(xiàn)存在一定的差異,主要是由各國(guó)生物區(qū)系存在的差異造成的.

    [1]US EPA.Guidelines for deriving numerical national water quality criteria for the protection of aquatic organisms and their uses[R].Washington DC:Office of Research and Development,1985:1-57.

    [2]孟偉,張遠(yuǎn),鄭丙輝.水環(huán)境質(zhì)量基準(zhǔn)、標(biāo)準(zhǔn)與流域水污染物總量控制策略[J].環(huán)境科學(xué)研究,2006,19(3):1-6.

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

    [4]吳豐昌,萬國(guó)江,蔡玉蓉.沉積物-水界面的生物地球化學(xué)作用[J].地球科學(xué)進(jìn)展,1996,11(2):191-197

    [5]WU F C,XING B S.Natural organic matter and its significance in the environment[M].Beijing:Chinese Science Publisher,2009:135-237.

    [6]吳豐昌.天然有機(jī)質(zhì)及其與污染物的相互作用[M].北京:科學(xué)出版社,2010:212-245.

    [7]孟偉,吳豐昌.水質(zhì)基準(zhǔn)的理論與方法學(xué)導(dǎo)論[M].北京:科學(xué)出版社,2010:23-45.

    [8]張翠,翟毓秀,寧勁松,等.鎘在水生動(dòng)物體內(nèi)的研究概況[J].水產(chǎn)科學(xué),2007,26(8):465-470.

    [9]趙紅霞,詹勇,許梓榮.重金屬對(duì)水生動(dòng)物毒性的研究進(jìn)展:(一)[J].內(nèi)陸水產(chǎn),2003,28(1):38-40.

    [10]閆振廣,孟偉,劉征濤,等.我國(guó)淡水水生生物鎘基準(zhǔn)研究[J].環(huán)境科學(xué)學(xué)報(bào),2009,29(11):2393-2406.

    [11]曹宇靜,吳豐昌.淡水中重金屬鎘的水質(zhì)基準(zhǔn)制定[J].安徽農(nóng)業(yè)科學(xué),2010,38(3):1378-1380,1437.

    [12]金相燦.中國(guó)湖泊環(huán)境[M].北京:海洋出版社,1995:142-234.

    [13]Canadian Council of Ministers of the Environment.Protocol for the derivation of water quality guidelines for the projection of aquatic life[R].Winnipeg,Manitoba:Canadian Council of Ministers of the Environment,1999:1-10.

    [14]US EPA.40CFR part 797 Environmental effects testing guidelines[R].Washington DC:Office of Science and Technology,1987:1-35.

    [15]吳本富.Cd2+和Pb2+重金屬離子對(duì)4種水生動(dòng)物的毒性研究[D].蕪湖:安徽師范大學(xué),2007:6-14.

    [16]付榮恕,杜作濱.鉛、鎘污染對(duì)水絲蚓的急性毒性效應(yīng)[J].山東師范大學(xué)學(xué)報(bào):自然科學(xué)版,2008,23(4):93-95.

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

    [18]閆沁遠(yuǎn),祁建紅.Cd對(duì)孝感地區(qū)鯽魚的急性毒性及安全濃度試驗(yàn)[J].湛江師范學(xué)院學(xué)報(bào),2008,29(3):77-80.

    [19]楊麗華.重金屬(鎘、銅、鋅和鉻)對(duì)鯽魚的生物毒性研究[D].廣州:華南師范大學(xué),2003:27-33.

    [20]陳芳,周啟星.城市土-水界面污染流條件下加樂麝香和鎘對(duì)鯽魚的聯(lián)合毒性[J].浙江大學(xué)學(xué)報(bào):農(nóng)業(yè)與生命科學(xué)版,2009,35(2):228-236.

    [21]王銀秋,張迎梅,趙東芹.重金屬鎘、鉛、鋅對(duì)鯽魚和泥鰍的毒性[J].甘肅科學(xué)學(xué)報(bào),2003,15(1):35-38.

    [22]侯麗萍,馬廣智.鎘與鋅對(duì)草魚種的急性毒性和聯(lián)合毒性研究[J].淡水漁業(yè),2002,32(3):44-46.

    [23]柏世軍.水環(huán)境鎘對(duì)羅非魚的毒性作用和機(jī)理探討[D].杭州:浙江大學(xué),2006:53-55.

    [24]王瑞龍,馬廣智,方展強(qiáng).銅、鎘、鋅對(duì)唐魚的急性毒性和安全濃度評(píng)價(jià)[J].水產(chǎn)科學(xué),2006,25(3):117-120.

    [25]劉曉旭,施蔡雷,賈秀英.Cu2+,Cd2+對(duì)甌江彩鯉的急性毒性研究[J].杭州師范大學(xué)學(xué)報(bào):自然科學(xué)版,2009,8(4): 304-307.

    [26]周輝明,吳志強(qiáng),袁樂洋,等.三種重金屬對(duì)鯉魚幼魚的毒性和積累[J].南昌大學(xué)學(xué)報(bào):理科版,2005,29(3):292-295.

    [27]黃敏毅,張育輝.鎘對(duì)中國(guó)林蛙蝌蚪生長(zhǎng)發(fā)育的毒性效應(yīng)[J].生態(tài)學(xué)雜志,2006,25(5):535-540.

    [28]趙華斌,胡好遠(yuǎn),王瑩,等.Pb2+和 Cd2+對(duì)水螅毒性作用的初步研究[J].安徽師范大學(xué)學(xué)報(bào):自然科學(xué)版,2005,28 (4):447-450.

    [29]畢蕾,紀(jì)文繡,唐金玉,等.不同水溫條件下重金屬對(duì)三角帆蚌幼蚌的急性致毒效應(yīng)[J].安徽農(nóng)業(yè)科學(xué),2009,37 (14):6468-6471.

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

    [31]何斌,陳先均,李孟均.Cu2+,Hg2+,Cd2+,Pb2+對(duì)淡水石斑魚仔魚的急性毒性研究[J].水產(chǎn)養(yǎng)殖,2006,27(2):1-3.

    [32]陳芳,周啟星.模擬城市徑流中加樂麝香和鎘對(duì)大型水蚤的毒性效應(yīng)[J].中國(guó)環(huán)境科學(xué),2009b,29(1):58-62.

    [33]黃寧昌.鎘對(duì)溞狀蚤的毒性試驗(yàn)[J].水產(chǎn)養(yǎng)殖,1994,15 (4):19-22.

    [34]陳孝煊,吳志新,操玉濤,等.紅鰲鰲蝦蝦苗對(duì)四種重金屬離子的耐受性[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào),1999,18(5):476-478.

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

    [36]戴習(xí)林,臧維林,楊鴻山,等.Cu2+,Zn2+,Cd2+對(duì)羅氏沼蝦幼蝦的毒性作用[J].上海水產(chǎn)大學(xué)學(xué)報(bào),2001,10(4): 298-300.

    [37]KRAAK M H S,TOUSSAINT M,LAVY D,et al.Short-term effects of metals on the filtration rate of the Zebra Mussel Dreissena polymorpha[J].Environ Pollut,1994,84(2):139-143.

    [38]KRAAK M H S,LAVY D,TOUSSAINT M,et al.Toxicity of heavy metals to the Zebra Mussel(Dreissena polymorpha)[M]//NALEPA T F,SCHLOESSER D M.Zebra Mussels: biology,impacts,and control(Chapter 29).Boca Raton:Lewis Publishers,1993:491-502.

    [39]曾麗璇,陳桂珠,吳宏海.重金屬鎘和銅對(duì)河蜆呼吸和排泄的毒性研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2007,26(1):175-178.

    [40]賈秀英.鎘對(duì)泥鰍幼魚的急性和亞急性研究[J].環(huán)境污染與防治,2001,23(5):227-228.

    [41]ONER M,ATLI G,CANLIM.Effects of metal(Ag,Cd,Cr,Cu,Zn)exposures on some enzymatic and non-enzymatic indicators in the liver of Oreochromis niloticus[J].Bulletin of Environmental Contamination and Toxicology,2009,2(3): 317-321.

    [42]SILVESTRE F,DIERICK J F,DUMONT V,et al.Differential protein expression profiles in anterior gills of Eriocheir sinensis during acclimation to cadmium[J].Aquatic Toxicology,2006,76(1):46-58.

    [43]GAMA F J L,CASTELLANOS P M E,SARMA S S S,et al. Effect of pulsed exposure to heavy metals(Copper and Cadmium)on some population variables of Brachionus calyciflorus Pallas(Rotifera:Brachionidae:Monogononta)[J]. Hydrobiologia,2007,593(1):201-208.

    [44]NIEDERLEHNER B.Cadmium toxicity to a cladoceran[M]. Virginia:Virginia Polytechnic Inst.and State University,1984: 9-130.

    [45]SUEDEL B C,RODGERS J J H,DEAVER E.Experimental factors that may affect toxicity of cadmium to freshwaterorganisms[J].Arch Environ Contam Toxicol,1997,33(2): 188-193.

    [46]GUILHERMINO L,SOBRAL O,CHASTINET C,et al.A Daphnia magna first-brood chronic test:an alternative to the conventional21-day chronic bioassay?[J].Ecotoxicol Environ Safety,1999,42(1):67-74.

    [47]CAVAS T,GARANKO N N,ARKHIPCHUK V V.Induction of micronuclei and binuclei in blood,gill and liver cells of fishes subchronically exposed to cadmium chloride and copper sulphate[J].Food and Chemical Toxicology,2005,43(4): 569-574.

    [48]ARKHIPCHUK V V,GARANKO N N.Using the nucleolar biomarker and the micronucleus test on in vivo fish fin cells[J].Ecotoxicol Environ Safety,2005,62(1):42-52.

    [49]GOETTL J P J,DAVIES P H,SINLEY J R.Water pollution studies[M].Colorado:Colorado Division of W ildlife,Boulder,CO,1976:68-75.

    [50]DAVIES P H,GORMANW C,CARLSON C A,et al.Effect of hardness on bioavailability and toxicity of cadmium to Rainbow Trout[J].Chemical Speciation and Bioavailabilty,1993,5 (2):67-77.

    [51]MEBANE C A,HENNESSY D P,DILLON F S.Developing acute-to-chronic toxicity ratios for Lead,Cadmium,and Zinc using rainbow trout,a mayfly,and a midge[J].Water Air Soil Pollut,2007,188(4):41-66.

    [52]HOLDWAY D A,LOK K,SEMAAN M.The acute and chronic toxicity of cadmium and zinc to two hydra species[J]. Environmental Toxicology,2001,16(6):557-565.

    [53]TOLLETT V D,BENVENUTTI E L,DEER L A,et al. Differential toxicity to Cd,Pb,and Cu in dragonfly larvae (Insecta:Odonata)[J].Arch Environ Contam Toxicol,2009,56(1):77-84.

    [54]王模善,趙鐵銘.重金屬鎘對(duì)沉水植物毒性效應(yīng)的研究[J].西南大學(xué)學(xué)報(bào):自然科學(xué)版,2008,30(4):128-134.

    [55]邱昌恩,況琪軍,畢永紅,等.Cd2+對(duì)綠球藻生長(zhǎng)及生理特性的影響研究[J].水生生物學(xué)報(bào),2007,31(1):142-145.

    [56]周宏,項(xiàng)斯瑞.重金屬銅、鋅、鉛、鎘對(duì)小形月牙藻生長(zhǎng)及亞顯微結(jié)構(gòu)的影響[J].杭州大學(xué)學(xué)報(bào):自然科學(xué)版,1998,25 (2):85-92.

    [57]孟曉紅.金屬鎘在魚體中的生物富集作用[J].廣東微量元素科學(xué),1997,3(1):8-11.

    [58]褚武英.幾種重金屬在淡水動(dòng)物體內(nèi)的富集及其毒理研究[D].上海:上海師范大學(xué),2003:23-29.

    [59]國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 18406—2001農(nóng)產(chǎn)品質(zhì)量安全:無公害水產(chǎn)品安全要求[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2001.

    [60]周永欣,王士達(dá),夏宜琤.水生生物與環(huán)境保護(hù)[M].北京:科學(xué)出版社,1983:10-50.

    [61]KOOIJMAN S.A safety factor for LC50values allowing for differences in sensitivity among species[J].Water Res,1987,21(3):269-276.

    [62]VAN STRAALEN N M,DENNEMAN C A J.Ecotoxicological evaluation of soil quality criteria[J].Ecotoxicol Environ Safety,1989,18(3):241-251.

    [63]ALDENBERG T,SLOB W.Confidence limits for hazardous concentrations based on logistically distributed NOEC toxicity data[J].Ecotoxicol Environ Safety,1993,25(1):48-63.

    [64]Organisation for Economic Co-operation and Development. Guidance document for aquatic effects assessment[R].Paris: OECD Environment Monographs No 92. Organisation for Economic Co-operation and Development,1995:5-85.

    [65]Environment Canada.Guidance document on statisticalmethods for environmental toxicity tests[R].Ottawa:Environment Protection Series Report EPS 1/RM/46,2005:21-30.

    [66]雷炳莉,金小偉,黃圣彪,等.太湖流域3種氯酚類化合物水質(zhì)基準(zhǔn)的探討[J].生態(tài)毒理學(xué)報(bào),2009,4(1):40-49.

    [67]SMITH E P,CAIRNS J J.Extrapolation methods for setting ecological standards for water quality:statistical and ecological concerns[J].Ecotoxicology,1993,2(3):203-219.

    [68]SHAO Q.Estimation for hazardous concentrations based on NOEC data:an alternative approach[J].Envirometrics,2000,11(5):583-595.

    [69]WARNE M S.Derivation of the Australian and New Zealand water quality guidelines for toxicants[J].Australasian Journal of Ecotoxicology,2001,7(2):123-136.

    [70]US EPA.National recommended water quality criteria[R]. Washington DC:Office of Water,office of Science and Technology,2002:6-29.

    [71]WU F C,TANOUE E.Isolation and partial characterization of dissolved copper-comp lexing ligands in streamwaters[J]. Environ Sci Technol,2001,35:3646-3652.

    [72]吳豐昌,孟偉,張瑞卿,等.保護(hù)淡水水生生物硝基苯水質(zhì)基準(zhǔn)研究[J].環(huán)境科學(xué)研究,2011,24(1):1-10.

    [73]WU F C,MENG W,ZHAO X l,et al.China embarking on development of its own national water quality criteria system[J].Environ Sci Technol,2010,44(21):7792-7793.

    [74]WU F C,TANOUE E.Molecular mass distribution and fluorescence characteristics of dissolved organic ligands for copper(Ⅱ)in Lake Biwa,Japan[J].Organic Geochemistry,2001,32:11-20.

    Derivation of Aquatic Life Water Quality Criteria for Cadm ium in Freshwater in China

    WU Feng-chang1,MENG Wei1,CAO Yu-jing1,3,LI Hui-xian1,ZHANG Rui-qing1,2,F(xiàn)ENG Cheng-lian1,YAN Zhen-guang1

    1.State Environmental Protection Key Laboratory for Lake Pollution Control,Chinese Research Academy of Environmental Sciences,Beijing 100012,China
    2.Guangzhou Institute of Geochemistry,Chinese Academy of Sciences,Guangzhou 510640,China
    3.China Academy of Railway Sciences,Energy Saving& Environmental Protection& Occupational Safety and Health Research Institute,Beijing 100081,China

    X-651

    A

    1001-6929(2011)02-0172-13

    2010-09-30

    2010-10-19

    國(guó)家 重 點(diǎn) 基 礎(chǔ) 研 究 發(fā) 展 計(jì) 劃 (973)項(xiàng) 目(2008CB418200)

    吳豐昌 (1965-),男,浙江衢州人,研究員,博士,博導(dǎo),主要從事環(huán)境地球化學(xué)和水質(zhì)基準(zhǔn)研究,wufengchang@vip.sk leg.cn.

    *責(zé)任作者,孟偉(1956-),男,山東青島人,中國(guó)工程院院士,博士,博導(dǎo),主要從事流域水環(huán)境管理研究,mengwei@craes.org.cn

    猜你喜歡
    基準(zhǔn)值水生敏感度
    河北省啤酒行業(yè)清潔生產(chǎn)水平分析
    基于5G用戶體驗(yàn)的業(yè)務(wù)質(zhì)量?jī)?yōu)化模型研究及其應(yīng)用
    全體外預(yù)應(yīng)力節(jié)段梁動(dòng)力特性對(duì)于接縫的敏感度研究
    電視臺(tái)記者新聞敏感度培養(yǎng)策略
    新聞傳播(2018年10期)2018-08-16 02:10:16
    一種基于改進(jìn)差分的測(cè)井?dāng)?shù)據(jù)可逆變長(zhǎng)碼壓縮方法
    在京韓國(guó)留學(xué)生跨文化敏感度實(shí)證研究
    三氯生對(duì)4種水生生物的急性毒性研究
    Diodes高性能汽車霍爾效應(yīng)閉鎖提供多種敏感度選擇
    依托科技創(chuàng)新 打造現(xiàn)代水生蔬菜產(chǎn)業(yè)
    跟我進(jìn)城的兄弟
    椰城(2013年4期)2013-11-16 06:31:11
    亚洲av成人一区二区三| 69精品国产乱码久久久| 妹子高潮喷水视频| 国产精品亚洲av一区麻豆| 99国产精品一区二区蜜桃av| 国产精品久久久久久亚洲av鲁大| 真人做人爱边吃奶动态| 国内毛片毛片毛片毛片毛片| av中文乱码字幕在线| 久久久久久亚洲精品国产蜜桃av| 99国产综合亚洲精品| 女人高潮潮喷娇喘18禁视频| 成人精品一区二区免费| 午夜老司机福利片| 88av欧美| 老司机深夜福利视频在线观看| 欧美国产日韩亚洲一区| 制服诱惑二区| 亚洲中文av在线| 亚洲精品久久成人aⅴ小说| 欧美av亚洲av综合av国产av| 天天一区二区日本电影三级 | 亚洲天堂国产精品一区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 51午夜福利影视在线观看| 免费av毛片视频| 久久人人97超碰香蕉20202| 久久国产精品影院| 精品国产国语对白av| 啦啦啦观看免费观看视频高清 | 欧美日韩精品网址| 一级毛片高清免费大全| 亚洲av电影不卡..在线观看| 免费在线观看亚洲国产| 电影成人av| 黄频高清免费视频| 国产欧美日韩精品亚洲av| 亚洲自拍偷在线| 日韩欧美国产一区二区入口| 又紧又爽又黄一区二区| 热re99久久国产66热| 成熟少妇高潮喷水视频| 亚洲人成电影免费在线| 免费人成视频x8x8入口观看| 十八禁网站免费在线| 一区福利在线观看| 一二三四在线观看免费中文在| 窝窝影院91人妻| 亚洲av电影不卡..在线观看| 一本久久中文字幕| 俄罗斯特黄特色一大片| 国产一区二区激情短视频| 精品国产一区二区三区四区第35| 一进一出抽搐gif免费好疼| 99国产精品一区二区三区| 黄片播放在线免费| 亚洲第一欧美日韩一区二区三区| 男人的好看免费观看在线视频 | 美女 人体艺术 gogo| 日韩av在线大香蕉| 宅男免费午夜| 一级a爱片免费观看的视频| 老熟妇仑乱视频hdxx| 黄色成人免费大全| 色播在线永久视频| 久久人妻av系列| 亚洲国产中文字幕在线视频| 一级作爱视频免费观看| 1024视频免费在线观看| 国产精品av久久久久免费| 国产精品久久久人人做人人爽| 最好的美女福利视频网| 一级黄色大片毛片| 免费一级毛片在线播放高清视频 | 亚洲精品一卡2卡三卡4卡5卡| 天堂影院成人在线观看| 午夜精品久久久久久毛片777| 午夜免费激情av| 高清黄色对白视频在线免费看| 亚洲性夜色夜夜综合| 99re在线观看精品视频| 亚洲精品在线美女| 亚洲最大成人中文| 窝窝影院91人妻| 啦啦啦韩国在线观看视频| 亚洲欧洲精品一区二区精品久久久| 制服诱惑二区| 成人手机av| 久久久久国内视频| 欧美在线一区亚洲| 久久久水蜜桃国产精品网| 欧美在线黄色| 亚洲美女黄片视频| 97人妻天天添夜夜摸| 高潮久久久久久久久久久不卡| 日本三级黄在线观看| 人人妻人人澡人人看| 这个男人来自地球电影免费观看| 免费看美女性在线毛片视频| 淫秽高清视频在线观看| 亚洲欧美日韩另类电影网站| 久久精品成人免费网站| 电影成人av| 一夜夜www| 国产av精品麻豆| 长腿黑丝高跟| 琪琪午夜伦伦电影理论片6080| 男女下面进入的视频免费午夜 | 18禁黄网站禁片午夜丰满| 欧美成狂野欧美在线观看| 99riav亚洲国产免费| 精品国产乱码久久久久久男人| 久久中文字幕人妻熟女| 午夜福利免费观看在线| 变态另类丝袜制服| 男女做爰动态图高潮gif福利片 | 欧美乱码精品一区二区三区| 免费搜索国产男女视频| 亚洲精品久久成人aⅴ小说| 婷婷精品国产亚洲av在线| 国产高清videossex| 十分钟在线观看高清视频www| 亚洲欧美精品综合久久99| 69精品国产乱码久久久| 午夜精品在线福利| 国产精华一区二区三区| 亚洲九九香蕉| 亚洲av日韩精品久久久久久密| 亚洲男人天堂网一区| 丰满的人妻完整版| 久久精品国产亚洲av香蕉五月| 日本 av在线| 欧美午夜高清在线| 男人舔女人下体高潮全视频| 亚洲国产看品久久| 亚洲国产看品久久| 黑人巨大精品欧美一区二区蜜桃| 国产一级毛片七仙女欲春2 | 免费在线观看日本一区| 国产av一区在线观看免费| 女人精品久久久久毛片| 老司机午夜十八禁免费视频| 别揉我奶头~嗯~啊~动态视频| 国产成人精品久久二区二区91| 91av网站免费观看| 免费高清在线观看日韩| 黄色视频不卡| 国产欧美日韩一区二区三| 一区福利在线观看| 午夜福利一区二区在线看| 国产真人三级小视频在线观看| 人人妻人人澡欧美一区二区 | 久久久久久久久久久久大奶| 精品久久久久久久毛片微露脸| 视频区欧美日本亚洲| 十八禁网站免费在线| 午夜免费观看网址| 亚洲男人天堂网一区| 亚洲男人天堂网一区| 1024香蕉在线观看| 国产精品久久久久久精品电影 | 老司机在亚洲福利影院| 超碰成人久久| 日韩精品中文字幕看吧| 国产精品99久久99久久久不卡| 无限看片的www在线观看| 精品久久久精品久久久| 熟妇人妻久久中文字幕3abv| 国产亚洲精品一区二区www| 淫妇啪啪啪对白视频| 成人欧美大片| 精品国产国语对白av| 两性午夜刺激爽爽歪歪视频在线观看 | 国产熟女午夜一区二区三区| 国产成人精品久久二区二区91| 久久久久久大精品| 琪琪午夜伦伦电影理论片6080| 日韩 欧美 亚洲 中文字幕| 精品国产一区二区三区四区第35| 日韩欧美一区视频在线观看| 国产三级在线视频| 99re在线观看精品视频| 亚洲精品一区av在线观看| 男人舔女人的私密视频| 精品人妻1区二区| 久久国产精品人妻蜜桃| 国产精华一区二区三区| www.www免费av| 母亲3免费完整高清在线观看| 18美女黄网站色大片免费观看| 日本a在线网址| a级毛片在线看网站| 国产国语露脸激情在线看| 91麻豆av在线| 日韩有码中文字幕| 三级毛片av免费| 成熟少妇高潮喷水视频| 亚洲欧洲精品一区二区精品久久久| 免费在线观看日本一区| 一本综合久久免费| 国产精品98久久久久久宅男小说| 精品不卡国产一区二区三区| 国产精品久久久久久亚洲av鲁大| 午夜影院日韩av| 日韩精品免费视频一区二区三区| 国产高清激情床上av| 亚洲激情在线av| 国产精品九九99| 丝袜美腿诱惑在线| 最近最新中文字幕大全免费视频| 性色av乱码一区二区三区2| 国产成人精品无人区| 中国美女看黄片| 男人操女人黄网站| 亚洲免费av在线视频| 国产成人av教育| 国产精品亚洲av一区麻豆| av天堂久久9| 国产99久久九九免费精品| 国产精品电影一区二区三区| 久久国产亚洲av麻豆专区| 欧美乱色亚洲激情| www日本在线高清视频| 国语自产精品视频在线第100页| 欧美日韩黄片免| 久久久久亚洲av毛片大全| 欧美黑人精品巨大| 免费在线观看完整版高清| 视频区欧美日本亚洲| 久久人妻熟女aⅴ| 制服诱惑二区| 人人妻人人澡人人看| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看 | 老汉色av国产亚洲站长工具| 成在线人永久免费视频| 男人的好看免费观看在线视频 | 国产一级毛片七仙女欲春2 | 久久婷婷成人综合色麻豆| 韩国精品一区二区三区| 国产又色又爽无遮挡免费看| 在线天堂中文资源库| 亚洲精品中文字幕一二三四区| 电影成人av| 人人妻人人澡人人看| 亚洲免费av在线视频| 久久影院123| 国产免费男女视频| 国产99白浆流出| 好男人电影高清在线观看| 欧美在线黄色| 美女高潮喷水抽搐中文字幕| 亚洲无线在线观看| 国产亚洲精品久久久久5区| 一个人观看的视频www高清免费观看 | 怎么达到女性高潮| 免费在线观看视频国产中文字幕亚洲| 亚洲自偷自拍图片 自拍| 一区二区三区激情视频| 亚洲美女黄片视频| 婷婷六月久久综合丁香| 成人亚洲精品一区在线观看| 1024视频免费在线观看| 国产欧美日韩一区二区精品| 操美女的视频在线观看| 级片在线观看| 欧美绝顶高潮抽搐喷水| 一级a爱视频在线免费观看| 大型黄色视频在线免费观看| 国产精品久久久人人做人人爽| 制服丝袜大香蕉在线| 给我免费播放毛片高清在线观看| 亚洲免费av在线视频| 淫妇啪啪啪对白视频| 看片在线看免费视频| 少妇 在线观看| 欧美日本亚洲视频在线播放| 在线永久观看黄色视频| 成在线人永久免费视频| 久9热在线精品视频| 人妻久久中文字幕网| 精品国内亚洲2022精品成人| 日韩国内少妇激情av| 久久婷婷成人综合色麻豆| 一级a爱片免费观看的视频| 国产亚洲欧美在线一区二区| 最近最新免费中文字幕在线| 色av中文字幕| 日韩欧美在线二视频| 久久国产精品人妻蜜桃| 女性被躁到高潮视频| 国产男靠女视频免费网站| 亚洲成人国产一区在线观看| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 91老司机精品| www.自偷自拍.com| 90打野战视频偷拍视频| 久久精品成人免费网站| 国产区一区二久久| 亚洲在线自拍视频| 午夜两性在线视频| 日韩精品免费视频一区二区三区| 国内毛片毛片毛片毛片毛片| 午夜老司机福利片| 日韩大尺度精品在线看网址 | 久久中文字幕人妻熟女| 男人的好看免费观看在线视频 | 日韩成人在线观看一区二区三区| 亚洲精品粉嫩美女一区| 亚洲人成电影观看| 自拍欧美九色日韩亚洲蝌蚪91| 女人爽到高潮嗷嗷叫在线视频| 欧美午夜高清在线| 桃色一区二区三区在线观看| 超碰成人久久| 亚洲熟妇中文字幕五十中出| 这个男人来自地球电影免费观看| 国产精品久久电影中文字幕| 夜夜看夜夜爽夜夜摸| 亚洲精品国产精品久久久不卡| 亚洲成国产人片在线观看| 久久国产精品影院| 亚洲精品国产精品久久久不卡| www.www免费av| √禁漫天堂资源中文www| 一本大道久久a久久精品| 亚洲伊人色综图| 国产午夜精品久久久久久| 成人手机av| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 国产又色又爽无遮挡免费看| 成人免费观看视频高清| 在线观看日韩欧美| 黄色a级毛片大全视频| 国产乱人伦免费视频| 午夜福利一区二区在线看| av视频在线观看入口| 久久精品国产99精品国产亚洲性色 | 国内精品久久久久久久电影| 我的亚洲天堂| 纯流量卡能插随身wifi吗| 亚洲av电影在线进入| 极品人妻少妇av视频| 国产伦一二天堂av在线观看| 巨乳人妻的诱惑在线观看| 久久久久久久精品吃奶| 成人18禁在线播放| 欧美精品亚洲一区二区| 18禁国产床啪视频网站| 成年人黄色毛片网站| 久久 成人 亚洲| 日本五十路高清| 午夜亚洲福利在线播放| 多毛熟女@视频| 又紧又爽又黄一区二区| 亚洲熟妇中文字幕五十中出| 好男人在线观看高清免费视频 | 国产精品一区二区免费欧美| av免费在线观看网站| 国产蜜桃级精品一区二区三区| 一本久久中文字幕| 午夜免费鲁丝| 国产亚洲精品第一综合不卡| 日本免费a在线| 精品午夜福利视频在线观看一区| 中文字幕精品免费在线观看视频| 啦啦啦免费观看视频1| 国产精品 国内视频| 在线播放国产精品三级| 在线观看免费视频日本深夜| 宅男免费午夜| 久9热在线精品视频| 搡老熟女国产l中国老女人| 天天躁狠狠躁夜夜躁狠狠躁| 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 国产野战对白在线观看| 黄频高清免费视频| av天堂久久9| 久久婷婷成人综合色麻豆| АⅤ资源中文在线天堂| 国产97色在线日韩免费| 激情在线观看视频在线高清| 久久中文看片网| 久久久国产成人免费| 久久精品亚洲精品国产色婷小说| 午夜福利一区二区在线看| 国产欧美日韩一区二区精品| 香蕉久久夜色| avwww免费| 久9热在线精品视频| 色播在线永久视频| 丁香欧美五月| 大香蕉久久成人网| 热re99久久国产66热| 国产男靠女视频免费网站| 91av网站免费观看| 久久久久精品国产欧美久久久| 精品国产超薄肉色丝袜足j| 久久精品影院6| 日韩一卡2卡3卡4卡2021年| 在线观看舔阴道视频| 精品一品国产午夜福利视频| 一区二区三区精品91| 亚洲精品av麻豆狂野| 91成人精品电影| 啪啪无遮挡十八禁网站| 免费观看精品视频网站| 欧美亚洲日本最大视频资源| av在线天堂中文字幕| 中文字幕最新亚洲高清| 91字幕亚洲| 亚洲一码二码三码区别大吗| 怎么达到女性高潮| 波多野结衣av一区二区av| 给我免费播放毛片高清在线观看| 在线观看舔阴道视频| 欧美老熟妇乱子伦牲交| 狂野欧美激情性xxxx| 亚洲av片天天在线观看| 嫩草影视91久久| 69av精品久久久久久| 两个人看的免费小视频| 此物有八面人人有两片| 久久久久久久精品吃奶| 精品国产乱子伦一区二区三区| 91成人精品电影| 午夜免费激情av| 精品久久久精品久久久| 午夜免费观看网址| 别揉我奶头~嗯~啊~动态视频| 91在线观看av| 国产私拍福利视频在线观看| 人人妻人人澡人人看| 国产精品1区2区在线观看.| 变态另类丝袜制服| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区三区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区在线不卡| 最好的美女福利视频网| 黄片大片在线免费观看| 熟妇人妻久久中文字幕3abv| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 亚洲专区国产一区二区| 久久人妻福利社区极品人妻图片| 国内毛片毛片毛片毛片毛片| 色婷婷久久久亚洲欧美| 久久久久国产精品人妻aⅴ院| 久久香蕉精品热| 亚洲成人精品中文字幕电影| 狠狠狠狠99中文字幕| 天堂影院成人在线观看| 免费看a级黄色片| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 啦啦啦免费观看视频1| 国产熟女午夜一区二区三区| 欧美激情久久久久久爽电影 | 久99久视频精品免费| 久久精品亚洲熟妇少妇任你| 91精品三级在线观看| 人妻久久中文字幕网| 99国产极品粉嫩在线观看| 在线观看免费日韩欧美大片| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出| 中出人妻视频一区二区| 精品一区二区三区视频在线观看免费| 成人免费观看视频高清| 亚洲精品国产色婷婷电影| 亚洲自拍偷在线| 别揉我奶头~嗯~啊~动态视频| 757午夜福利合集在线观看| 美女高潮到喷水免费观看| 久久久久精品国产欧美久久久| 每晚都被弄得嗷嗷叫到高潮| 久久人妻福利社区极品人妻图片| 久久影院123| 成人三级黄色视频| 国产精品美女特级片免费视频播放器 | 亚洲成av人片免费观看| 日韩国内少妇激情av| 国产一区二区在线av高清观看| 人妻久久中文字幕网| 老汉色∧v一级毛片| 一级片免费观看大全| 9191精品国产免费久久| 亚洲,欧美精品.| 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| 日本精品一区二区三区蜜桃| 久久久久国产一级毛片高清牌| 伊人久久大香线蕉亚洲五| 视频在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品日韩av在线免费观看 | 国产精品久久久久久亚洲av鲁大| 欧美在线一区亚洲| 精品国产乱子伦一区二区三区| 两个人免费观看高清视频| 日本撒尿小便嘘嘘汇集6| 18禁观看日本| 男女下面插进去视频免费观看| 看黄色毛片网站| 黄片小视频在线播放| 国产成人欧美在线观看| 日韩精品青青久久久久久| av视频在线观看入口| 亚洲精品国产色婷婷电影| 久久国产精品男人的天堂亚洲| a级毛片在线看网站| 人人澡人人妻人| 男女做爰动态图高潮gif福利片 | 日日爽夜夜爽网站| 1024香蕉在线观看| 国产成人啪精品午夜网站| 一级a爱片免费观看的视频| 久久久久久亚洲精品国产蜜桃av| 午夜两性在线视频| 电影成人av| 在线十欧美十亚洲十日本专区| 精品久久久久久久毛片微露脸| 可以在线观看毛片的网站| 91精品三级在线观看| 成人18禁高潮啪啪吃奶动态图| 免费看a级黄色片| 老司机午夜十八禁免费视频| 纯流量卡能插随身wifi吗| 国产精品久久视频播放| 国产亚洲精品综合一区在线观看 | 日韩有码中文字幕| 色哟哟哟哟哟哟| 亚洲色图 男人天堂 中文字幕| 欧美精品啪啪一区二区三区| 久久中文字幕一级| 男男h啪啪无遮挡| 亚洲精品国产精品久久久不卡| 看免费av毛片| 国产精品二区激情视频| 黑人欧美特级aaaaaa片| 日韩精品免费视频一区二区三区| 精品日产1卡2卡| 欧美一区二区精品小视频在线| 久久午夜亚洲精品久久| 老汉色∧v一级毛片| av福利片在线| 不卡一级毛片| 国产精华一区二区三区| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 成人国产综合亚洲| 怎么达到女性高潮| 他把我摸到了高潮在线观看| 亚洲第一电影网av| 日韩欧美一区视频在线观看| 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| 成年女人毛片免费观看观看9| 精品久久久久久,| 国产亚洲欧美98| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 久久久国产欧美日韩av| 亚洲成国产人片在线观看| 精品国产一区二区三区四区第35| 成熟少妇高潮喷水视频| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 少妇被粗大的猛进出69影院| 日本在线视频免费播放| 亚洲视频免费观看视频| 九色国产91popny在线| 岛国视频午夜一区免费看| 精品久久久久久久久久免费视频| 日本精品一区二区三区蜜桃| 国产成人av激情在线播放| 国产欧美日韩精品亚洲av| 日日夜夜操网爽| 国产一区二区在线av高清观看| 脱女人内裤的视频| 久久久久久国产a免费观看| 久久久久久久久免费视频了| 国产精品一区二区三区四区久久 | 日本一区二区免费在线视频| 嫩草影视91久久| 首页视频小说图片口味搜索| 搡老岳熟女国产| 亚洲国产中文字幕在线视频| 黄色女人牲交| 久久热在线av| 国产一区二区三区综合在线观看| 久久国产亚洲av麻豆专区| 久久久精品欧美日韩精品| 亚洲专区字幕在线| 69av精品久久久久久| 女同久久另类99精品国产91| 国产色视频综合| 999精品在线视频| 国内精品久久久久久久电影| 亚洲国产毛片av蜜桃av| 午夜久久久在线观看| 国产一卡二卡三卡精品| 黄色a级毛片大全视频| 国产一区二区在线av高清观看| 男男h啪啪无遮挡| 精品国产乱码久久久久久男人| www日本在线高清视频| 亚洲精品久久国产高清桃花| 一区二区三区激情视频|