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

    生物效應(yīng)比(BER)技術(shù)預(yù)測我國水生生物基準(zhǔn)探討

    2016-06-07 06:21:49王曉南閆振廣劉征濤王婉華中國環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室國家環(huán)境保護化學(xué)品生態(tài)效應(yīng)與風(fēng)險評估重點實驗室北京000中海石油環(huán)保服務(wù)有限公司天津30045
    中國環(huán)境科學(xué) 2016年1期
    關(guān)鍵詞:美國中國

    王曉南,閆振廣,劉征濤*,王婉華,張 聰(.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室,國家環(huán)境保護化學(xué)品生態(tài)效應(yīng)與風(fēng)險評估重點實驗室,北京 000;.中海石油環(huán)保服務(wù)有限公司,天津 30045)

    ?

    生物效應(yīng)比(BER)技術(shù)預(yù)測我國水生生物基準(zhǔn)探討

    王曉南1,閆振廣1,劉征濤1*,王婉華1,張 聰2(1.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室,國家環(huán)境保護化學(xué)品生態(tài)效應(yīng)與風(fēng)險評估重點實驗室,北京 100012;2.中海石油環(huán)保服務(wù)有限公司,天津 300452)

    摘要:基于中國與美國間生物敏感性的差異,對生物效應(yīng)比(BER)技術(shù)預(yù)測我國水生生物基準(zhǔn)值進行了研究.首先,對中美共有的且急性生物毒性數(shù)據(jù)符合3門8科的污染物進行篩選,然后,依據(jù)兩國物種的敏感性與代表性建立基于不同生物組合的BER技術(shù),最后,對不同生物組合BER技術(shù)預(yù)測水生生物急性基準(zhǔn)值與實測基準(zhǔn)值進行比較與分析,篩選出預(yù)測效果較好的BER技術(shù).結(jié)果表明:共篩選出9個中美共有且毒性數(shù)據(jù)豐富的污染物(As(III)、Cr(VI)、Hg、Cu、Zn、Pb、對硫磷、毒死蜱和三丁基錫),依據(jù)本土生物毒性數(shù)據(jù)推導(dǎo)出該9種污染物的我國水生生物急性基準(zhǔn)值分別為201.72、2.64、0.74、1.32、55.83、92.25、0.12、0.36和0.38μg/L.此外,在對7種生物組合方式BER技術(shù)的預(yù)測效果進行分析比較的基礎(chǔ)上,初步提出基于同屬或科的生物組合方式的BER技術(shù)可較好地對9種污染物的本土水生生物急性基準(zhǔn)值進行預(yù)測.研究結(jié)果可為在本土生物毒性數(shù)據(jù)缺乏時充分利用現(xiàn)有毒性數(shù)據(jù)或進行少量毒性試驗的基礎(chǔ)上對我國水生生物急性基準(zhǔn)值的預(yù)測提供幫助.

    關(guān)鍵詞:生物效應(yīng)比;水生生物基準(zhǔn);預(yù)測基準(zhǔn)值;中國;美國

    * 責(zé)任作者, 研究員, liuzt@craes.org.cn

    目前,我國正在系統(tǒng)地開展環(huán)境基準(zhǔn)技術(shù)方法研究[1-4],包括:水生生物基準(zhǔn)、生態(tài)學(xué)基準(zhǔn)、沉積物基準(zhǔn)、土壤環(huán)境基準(zhǔn)等.其中,在水生生物基準(zhǔn)技術(shù)方法研究方面,國內(nèi)學(xué)者對水生生物基準(zhǔn)推導(dǎo)過程中最少生物毒性數(shù)據(jù)需求[5]、水質(zhì)參數(shù)歸一化[6-7]、水生生物基準(zhǔn)受試物種的篩選[8-11]等技術(shù)方法進行了研究.在水生生物基準(zhǔn)值推導(dǎo)方面,國內(nèi)學(xué)者一般進行污染物對本土生物的生態(tài)毒理學(xué)試驗或進行生物毒性數(shù)據(jù)搜集,然后依據(jù)所得到的急慢性毒性數(shù)據(jù),采用相應(yīng)的方法推導(dǎo)這些污染物的水生生物基準(zhǔn)值,如:依據(jù)搜集的毒性數(shù)據(jù)先后推導(dǎo)了硝基苯[12]、鎘[13]、氨氮[14]、鋅[15]等污染物的水生生物基準(zhǔn)值;其他學(xué)者則進行污染物的生態(tài)毒理學(xué)試驗并結(jié)合搜集的毒性數(shù)據(jù)先后推導(dǎo)了氯酚類化合物[16]、三氯生[17]、硝基苯[18]、四溴雙酚-A[19]等污染物的水生生物基準(zhǔn)值.通過不斷地研究和發(fā)展,我國水生生物基準(zhǔn)研究工作逐步從對國外的相關(guān)技術(shù)方法進行引進和學(xué)習(xí),發(fā)展到對我國本土水生生物基準(zhǔn)方法學(xué)和技術(shù)體系進行探索和研究.

    在推導(dǎo)本土水生生物基準(zhǔn)時,美國[20]、澳大利亞[21]等國家推薦使用本土物種的毒性數(shù)據(jù),并且有學(xué)者對使用非本土物種毒性數(shù)據(jù)推導(dǎo)出的基準(zhǔn)值能否保護本土水生生物提出了質(zhì)疑[22].因此,各國在推導(dǎo)本國水生生物基準(zhǔn)時,基本均面臨本土生物毒性數(shù)據(jù)缺乏的問題.即使在水生生物基準(zhǔn)發(fā)展較早的美國,至今也才推薦了58種污染物(化學(xué)的、物理的、生物的)的水生生物基準(zhǔn)值.在我國,由于水生生物基準(zhǔn)研究工作起步較晚,各項相關(guān)研究尚不完善,并且我國也存在本土物種的生態(tài)毒理學(xué)數(shù)據(jù)普遍較缺乏的問題.因此,在推導(dǎo)水生生物基準(zhǔn)時,如何利用好現(xiàn)有的生物毒性數(shù)據(jù),如何解決本土毒性數(shù)據(jù)缺乏的問題,成為水生生物基準(zhǔn)研究中急需解決的重要問題.

    在推導(dǎo)水生生物基準(zhǔn)時,直接采用他國生物毒性數(shù)據(jù)或基準(zhǔn)值存在爭議[22],然而,如果將各類污染物的水生生物基準(zhǔn)都從頭開始研究、測試、推導(dǎo),短期內(nèi)難以發(fā)揮有效的科技支撐作用.本文以中國與美國為例,利用中美物種間的毒性效應(yīng)差異比,對生物毒性效應(yīng)比(BER)技術(shù)預(yù)測水生生物基準(zhǔn)值進行探討.設(shè)定BER等于“我國某一類群生物的毒性值”除以“美國同一類群生物的毒性值”,我國預(yù)測的水生生物基準(zhǔn)值等于美國基準(zhǔn)值乘以BER.該方法的探索,以期在本土生物毒性數(shù)據(jù)缺乏時,為我國水生生物基準(zhǔn)的研究與推導(dǎo)提供有用的技術(shù)信息.

    1 材料與方法

    1.1 生物毒性數(shù)據(jù)的搜集與篩選原則

    水生生物毒性數(shù)據(jù)的主要來源包括:中國知網(wǎng)CNKI、ECOTOX(http://cfpub.epa.gov/ ecotox/)等數(shù)據(jù)庫以及公開發(fā)表的文獻.

    數(shù)據(jù)搜集完成后,參照美國水生生物基準(zhǔn)指南[20],對搜集的生物毒性數(shù)據(jù)進行初步篩選.毒性數(shù)據(jù)的初步篩選原則包括:測試信息不完全、有疑點的數(shù)據(jù)(未設(shè)立對照組、對照組大量死亡等)、性質(zhì)不穩(wěn)定(易揮發(fā)、易水解)而未采用流水式試驗等的毒性數(shù)據(jù)不可用.完成初步篩選后,需進一步對毒性數(shù)據(jù)進行篩選,二次篩選原則包括:浮游甲殼類生物的急性毒性試驗終點為48-h LC50或EC50,魚、貝、蝦、蟹等生物的試驗終點是96-h LC50或EC50,優(yōu)先選擇流水式試驗結(jié)果及對試驗溶液濃度有監(jiān)控的毒性數(shù)據(jù);所有毒性數(shù)據(jù)都要求有明確的測試終點、測試時間及對測試階段或指標(biāo)的詳細(xì)描述.當(dāng)同一個物種或同一個終點有多個毒性值可用時,使用該多個值的幾何平均值SMAV(種平均急性值)[20].

    1.2 中美共有污染物的篩選

    US EPA至今共推薦了58種污染物(生物、化學(xué)、物理)的水生生物基準(zhǔn)值,由于對美國水生生物分類的認(rèn)識有限,本文在進行BER技術(shù)研究時,只采用了US EPA公開發(fā)布的58種污染物的毒性數(shù)據(jù)作為美國物種毒性數(shù)據(jù).為了進一步進行中美生物毒性效應(yīng)比BER的分析研究,依據(jù)1.1節(jié)的搜集和篩選原則,結(jié)合《中國動物志》[23]和中國生物物種名錄(http: data.sp2000.cn/2011_ cnnode_c/info_checklist.php),搜集和篩選該58種污染物的我國本土生物毒性數(shù)據(jù),然后,依據(jù)“三門八科”[20]最少生物毒性數(shù)據(jù)需求原則,篩選出中美共有的且毒性數(shù)據(jù)比較充足的污染物.

    1.3 BER技術(shù)預(yù)測水生生物急性基準(zhǔn)初步分析

    通過先前對水質(zhì)基準(zhǔn)受試生物篩選[8-11,24]的研究,結(jié)合US EPA物種種間關(guān)系估算模型[25]的研究,發(fā)現(xiàn)同類物種(如同目、同科、同屬、同種)對同個污染物的敏感性比較相似,以至US EPA在其物種種間關(guān)系估算模型中指出:要優(yōu)先選用生物學(xué)距離≤3(US EPA規(guī)定同目生物學(xué)距離為3、同科為2、同屬為1)[25]的替代物種進行估算來獲取更為可靠的毒性值.因此,從US EPA現(xiàn)有的水生生物基準(zhǔn)文件中可以得到不同生物分類對污染物的敏感性情況,借鑒US EPA物種種間關(guān)系估算方法、水效應(yīng)比技術(shù)[26],利用不同生物組合方式的中美物種間的敏感性差異,來對我國水生生物急性基準(zhǔn)進行初步估算和預(yù)測.BER技術(shù)的意義在于可充分利用現(xiàn)有的本土生物毒性數(shù)據(jù)或者進行少量的生物毒性試驗,對我國水生生物基準(zhǔn)進行估算和預(yù)測.

    針對1.2節(jié)篩選的中美共有典型污染物,選擇兩國同類的并且敏感的生物進行組合,依據(jù)不同生物組合進行BER(式1)計算,然后采用BER技術(shù)預(yù)測我國的HC5(5%物種危害濃度)或FAV(最終急性值)(式2),并對預(yù)測HC5(FAV)值和實測HC5(FAV)值進行數(shù)量級差異分析.前人研究提出,如果預(yù)測值和實測值間差異在一個數(shù)量級以內(nèi)則認(rèn)為差異較小、預(yù)測較好[27-32],超過一個數(shù)量級則認(rèn)為差異較大、預(yù)測較差.最后,依據(jù)不同生物組合BER技術(shù)的預(yù)測效果,初步提出我國BER預(yù)測技術(shù).

    式中:我國生物毒性數(shù)據(jù)為我國不同生物組合方式物種急性毒性值的幾何平均值,美國生物毒性數(shù)據(jù)為美國不同生物組合方式物種急性毒性值的幾何平均值.

    式中:美國水生生物基準(zhǔn)值為US EPA公布推薦的相關(guān)污染物的基準(zhǔn)值.

    1.4 實測水生生物基準(zhǔn)值的推導(dǎo)

    依據(jù)US EPA水生生物基準(zhǔn)指南[20],對1.2節(jié)篩選出的各污染物的GMAV(屬急性平均值)進行計算和排序(排序第幾為第R,最低的為1,最高的為N),采用下列公式(3~6)計算CMC(急性基準(zhǔn)值),式中P為累積概率P = R/(N + 1),L為選擇P值最接近0.05的4個GMAV和相應(yīng)的P值進行計算,FAV為最終急性值,可用于計算CMC值.

    2 結(jié)果與分析

    2.1 中美共有污染物的生物毒性數(shù)據(jù)

    依據(jù)1.2節(jié)篩選原則,最終篩選出9個中美共有且本土生物毒性數(shù)據(jù)豐富的污染物,包括:As(III)、Cr(VI)、Hg、Cu、Zn、Pb、對硫磷、毒死蜱和三丁基錫(TBT).受篇幅限制,只列出毒死蜱、三丁基錫的本土生物毒性數(shù)據(jù)和文獻來源(表1、表2),其中Hg、Cr(VI)、As(III)和對硫磷的數(shù)據(jù)見前人研究[24],Cu、Zn、Pb的數(shù)據(jù)參考前研究[30,33-34].美國共有污染物的生物毒性數(shù)據(jù)可在US EPA基準(zhǔn)文件[35]中查找.

    表1 毒死蜱的本土生物毒性數(shù)據(jù)Table 1 Toxicity data of native species for chlorpyrifos

    表2 三丁基錫的本土生物毒性數(shù)據(jù)Table 2 Toxicity data of native species for TBT

    2.2 不同生物組合方式的BER技術(shù)

    依據(jù)中美典型污染物現(xiàn)有的生物毒性數(shù)據(jù)和生物敏感性分布,結(jié)合我國生物分類情況,選出7種不同生物組合方式的BER技術(shù)并計算BER 值,然后對我國水生生物基準(zhǔn)值進行預(yù)測并對預(yù)測效果進行分析.具體組合方式如下:

    (1) 組合1:我國大型溞+鯉魚+鯽魚+羽搖蚊/美國基準(zhǔn)文件前4種[35].

    該組合關(guān)注的是中國重要物種和美國敏感性物種間差異:

    大型溞:美國MTDR推薦了一種浮游甲殼類生物[20],大型溞是其中的一種,并且該物種在我國分布廣泛且毒性數(shù)據(jù)充足;

    鯉魚、鯽魚:美國MTDR推薦了兩種魚類[20],我國魚類以鯉科為主,鯉科魚中以鯉魚、鯽魚分布廣泛且毒性數(shù)據(jù)較充足;

    羽搖蚊:美國MTDR推薦的在水生生物基準(zhǔn)推導(dǎo)過程中至少需要一種昆蟲類生物的毒性數(shù)據(jù)[20],其中羽搖蚊在我國分布廣泛且常被用于生態(tài)毒理學(xué)試驗.

    美國基準(zhǔn)文件前4種物種的名稱和生物分類可在文獻[35]中查找.

    (2) 組合2:與美國基準(zhǔn)文件前4屬同類的我國4屬/美國基準(zhǔn)文件前4屬,各污染物的屬組合結(jié)果見表3,我國缺少毒死蜱相應(yīng)4屬的毒性數(shù)據(jù),因此未對其進行研究.

    該組合關(guān)注的是中美敏感屬間的差異,在計算水生生物基準(zhǔn)值時,US EPA關(guān)注前4屬敏感性生物[20],它們對基準(zhǔn)值的大小將產(chǎn)生最直接的影響.

    表3 我國的4個敏感生物屬Table 3 The 4 sensitive genera in China

    (3) 組合3:我國4種與美國基準(zhǔn)文件前4個物種同種的生物(同時考慮我國生物分類情況、本土基準(zhǔn)受試生物篩選結(jié)果[8-11]和現(xiàn)有毒性數(shù)據(jù)情況)/美國基準(zhǔn)文件前4種,各污染物的物種組合結(jié)果見表4,無同種生物時,用同屬的其他物種替代.該組合關(guān)注的是中美敏感物種間的差異.

    表4 我國的4種敏感生物物種Table 4 The 4 sensitive species in China

    (4) 組合4:我國溞科大型溞+蚤狀溞+鋸頂?shù)皖~溞+模糊網(wǎng)紋溞/美國基準(zhǔn)文件中同類的4個溞科物種.

    該組合關(guān)注的是中美敏感性溞科物種間的差異,由2.1節(jié)本土生物毒性數(shù)據(jù)和美國生物毒性數(shù)據(jù)排序可知,淡水枝角類生物如溞科的溞屬、低額溞屬、網(wǎng)紋溞屬的物種對各類污染物的敏感性均排在前列,因此,本文嘗試對我國4種枝角類溞科敏感性生物(大型溞、蚤狀溞、模糊網(wǎng)紋溞、鋸頂?shù)最~溞,沒有毒性數(shù)據(jù)時可用同屬或同科的物種替代)和美國基準(zhǔn)文件中相似的4個溞科物種(沒有同種的數(shù)據(jù)時可用同屬同科的數(shù)據(jù)代替,同屬同科數(shù)據(jù)不足4個時可用與我國同種的數(shù)據(jù)重復(fù)一次計算;沒有溞科數(shù)據(jù)時可用文件中敏感的甲殼綱其他物種代替)(表5)進行BER分析.

    表5 美國基準(zhǔn)文件中的4種生物物種Table 5 The 4species in US EPA criteria documents

    (5) 組合5:指定4屬中各1種生物,同時考慮我國本土基準(zhǔn)受試生物篩選結(jié)果/美國基準(zhǔn)文件中同類的4屬生物.

    該組合關(guān)注的是中國重要物種和美國同類物種間差異,我國指定的4屬為甲殼綱溞科溞屬的大型溞、硬骨魚綱鯉科鯉屬的鯉魚、硬骨魚綱鯉科鯽屬的鯽魚和昆蟲綱搖蚊科搖蚊屬的羽搖蚊幼蟲.美國同類的4屬生物優(yōu)先選擇與我國同種生物,沒有同種生物時,依次以同屬、同科、同綱的其他生物代替,如果數(shù)據(jù)還是不足,以所缺屬的其他生物重復(fù)一次計算,各污染物的物種組合結(jié)果見表6.

    (6) 組合6:我國的2種鯉科魚/美國的2種鯉科魚.

    該組合關(guān)注的是中國重要鯉科物種和美國同類物種間差異,我國指定的2種鯉科魚為分布廣泛、毒性數(shù)據(jù)較充足的硬骨魚綱鯉科鯉屬的鯉魚、硬骨魚綱鯉科鯽屬的鯽魚.美國的2種鯉科魚優(yōu)先選擇與我國同種生物(鯉魚、鯽魚),沒有同種生物時,依次以同屬、同科的其他生物代替,如果數(shù)據(jù)還是不足,以所缺屬的其他生物重復(fù)一次計算,各污染物的物種組合結(jié)果見表7.

    表6 美國基準(zhǔn)文件中的4屬生物物種Table 6 Species of the 4genera in US EPA criteria documents

    表7 美國基準(zhǔn)文件中的2種鯉科魚Table 7 The 2cyprinid fish in US EPA criteria documents

    (7) 組合7:我國的2種鯉科魚/美國2種鮭科魚.

    該組合關(guān)注的是中國重要鯉科物種和美國重要鮭科物種間差異,本文中我國指定的2種鯉科魚為分布廣泛、毒性數(shù)據(jù)較充足的硬骨魚綱鯉科鯉屬的鯉魚、硬骨魚綱鯉科鯽屬的鯽魚.美國魚類以鮭科魚為主,鮭科魚毒性數(shù)據(jù)豐富,2種鮭科魚的選擇較多,為使分析比較更為可靠,本文優(yōu)先選用數(shù)據(jù)充足的硬骨魚綱鮭科的虹鱒和溪紅點鮭.本文對中美2種鯉科鮭科魚進行BER計算,探索中美兩國各自主要魚類物種用于BER技術(shù)的可行性.

    2.3 實測本土水生生物基準(zhǔn)值

    依據(jù)1.4節(jié)水生生物基準(zhǔn)值推導(dǎo)方法,結(jié)合2.1節(jié)搜集篩選到的本土水生生物毒性數(shù)據(jù),計算9種共有污染物As(III)、Cr(VI)、Hg、Cu、Zn、Pb、對硫磷、毒死蜱和三丁基錫的本土水生生物急性基準(zhǔn)值(表8).

    表8 9種污染物的本土水生生物基準(zhǔn)實測值Table 8 Measured aquatic life criteria for 9pollutants

    2.4 預(yù)測本土水生生物基準(zhǔn)值

    采用不同生物組合方式BER技術(shù)預(yù)測與估算本土水生生物基準(zhǔn)值的方法見1.3節(jié),以生物組合1的預(yù)測方法為例:組合1的BER值=大型溞、鯉魚、鯽魚與羽搖蚊毒性數(shù)據(jù)的幾何平均值/美國基準(zhǔn)文件前4種毒性數(shù)據(jù)的幾何平均值,預(yù)測本土水生生物基準(zhǔn)值=BER×美國水生生物基準(zhǔn)值.

    2.5 BER技術(shù)預(yù)測本土水生生物基準(zhǔn)初步分析

    本文對上述7種生物組合方式的BER技術(shù)進行了探索和分析,對基于實測毒性數(shù)據(jù)的基準(zhǔn)值和基于BER技術(shù)所預(yù)測的基準(zhǔn)值進行了推導(dǎo)和差異分析,依據(jù)不同生物組合的預(yù)測效果,初步提出我國BER預(yù)測技術(shù).7種生物組合方式的BER技術(shù)對9種污染物水生生物的預(yù)測效果如下:

    (1) 組合1:我國大型溞+鯉魚+鯽魚+羽搖蚊/美國基準(zhǔn)文件前4種的BER技術(shù)預(yù)測效果見圖1,從圖1可知,生物組合方式1所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在1個數(shù)量級以內(nèi)的有1 個,超過1個數(shù)量級的有8個,一般認(rèn)為差異在1個數(shù)量級以內(nèi)的預(yù)測效果較好,超過1個數(shù)量級則預(yù)測效果較差[24,27-32],可見生物組合方式1的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測存在較大的不確定性.

    圖1 生物組合1預(yù)測基準(zhǔn)值與實際值的比較Fig.1 Comparison of measured criteria and predicted criteria based on BER method of biological group 1圖中橫縱坐標(biāo)值為對實際基準(zhǔn)值和預(yù)測基準(zhǔn)值進行對數(shù)處理lg (μg/L)后的值,下圖同

    圖2 生物組合2預(yù)測基準(zhǔn)值與實際值的比較Fig.2 Comparison of measured criteria and predicted criteria based on BER method of biological group 2

    (2) 組合2:與US EPA基準(zhǔn)文件前4屬同類的我國4屬/美國基準(zhǔn)文件前4屬的BER技術(shù)的預(yù)測效果見圖2.從圖2可知,生物組合方式2所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在0.5個數(shù)量級以內(nèi)的有5個,1個數(shù)量級以內(nèi)的有3個,差異在1個數(shù)量級以內(nèi)的認(rèn)為預(yù)測效果較好[24,27-32],可見與生物組合方式1相比,生物組合方式2的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測較好.

    (3) 組合3:我國4種與美國基準(zhǔn)文件前4個物種同種的生物(同時考慮我國生物分類情況、本土基準(zhǔn)受試生物篩選結(jié)果和現(xiàn)有毒性數(shù)據(jù)情況)/美國基準(zhǔn)文件前4種的BER技術(shù)預(yù)測效果見圖3.從圖3可知,生物組合方式3所預(yù)測的基準(zhǔn)值與實測基準(zhǔn)值差異在1個數(shù)量級以內(nèi)的有8 個,超過1個數(shù)量級的有1個,差異在1個數(shù)量級以內(nèi)的認(rèn)為預(yù)測效果較好[24,27-32],超過1個數(shù)量級則預(yù)測效果較差,可見生物組合方式3和2的預(yù)測效果相似,均比生物組合方式1的BER技術(shù)的預(yù)測效果好,但存在一定的不確定性.

    圖3 生物組合3預(yù)測基準(zhǔn)值與實際值的比較Fig.3 Comparison of measured criteria and predicted criteria based on BER method of biological group 3

    (4) 組合4:我國溞科大型溞+蚤狀溞+鋸頂?shù)皖~溞+模糊網(wǎng)紋溞/美國基準(zhǔn)文件中同類的4個溞科物種的BER技術(shù)預(yù)測效果見圖4.從圖4可知,生物組合方式4所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在0.5個數(shù)量級以內(nèi)的有5個,1個數(shù)量級以內(nèi)的有4個,差異在1個數(shù)量級以內(nèi)則認(rèn)為預(yù)測效果較好[24,27-32],可見生物組合方式4的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測效果較好.

    (5) 組合5:我國指定4屬中各1種生物,同時考慮我國本土基準(zhǔn)受試生物篩選結(jié)果/美國同類的4屬生物的BER技術(shù)預(yù)測效果見圖5.從圖5可知,生物組合方式5所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在0.5個數(shù)量級以內(nèi)的有7個,1個數(shù)量級以內(nèi)的有2個,差異在1個數(shù)量級以內(nèi)則認(rèn)為預(yù)測效果較好[24,27-32],可見生物組合方式5的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測效果較好.

    圖4 生物組合4預(yù)測基準(zhǔn)值與實際值的比較Fig.4 Comparison of measured criteria and predicted criteria based on BER method of biological group 4

    圖5 生物組合5預(yù)測基準(zhǔn)值與實際值的比較Fig.5 Comparison of measured criteria and predicted criteria based on BER method of biological group 5

    (6) 組合6:我國的2種鯉科魚/美國的2種鯉科魚的BER技術(shù)預(yù)測效果見圖6.從圖6可知,生物組合方式6所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在1個數(shù)量級以內(nèi)的有7個,超過1個數(shù)量級的有2個,差異在1個數(shù)量級以內(nèi)則認(rèn)為預(yù)測效果較好[24,27-32],超過1個數(shù)量級則預(yù)測效果較差,可見生物組合方式6的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測效果存在一定的不確定性.

    圖6 生物組合6預(yù)測基準(zhǔn)值與實際值的比較Fig.6 Comparison of measured criteria and predicted criteria based on BER method of biological group 6

    (7) 組合7:我國的2種鯉科魚/美國的2種鮭科魚的BER技術(shù)預(yù)測效果見圖7.從圖7可知,生物組合方式7所預(yù)測的基準(zhǔn)值與實測推導(dǎo)的基準(zhǔn)值差異在1個數(shù)量級以內(nèi)的有8個,超過1個數(shù)量級的有1個,差異在1個數(shù)量級以內(nèi)則認(rèn)為預(yù)測效果較好[24,27-32],超過1個數(shù)量級則預(yù)測效果較差,可見生物組合方式7的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測效果比生物組合方式6相對要好,仍存在一定的不確定性.

    圖7 生物組合7預(yù)測基準(zhǔn)值與實際值的比較Fig.7 Comparison of measured criteria and predicted criteria based on BER method of biological group 7

    從上述的7種生物組合方式的BER技術(shù)對9種污染物水生生物基準(zhǔn)值預(yù)測效果來看,基于科、屬的生物組合方式2、4、5的BER技術(shù)可較好地對水生生物基準(zhǔn)值進行預(yù)測,多數(shù)預(yù)測值的差異在1個數(shù)量級內(nèi),甚至在0.5個數(shù)量級內(nèi),表明預(yù)測效果較好[24,27-32],而基于種的生物組合方式1、3、6、7在預(yù)測基準(zhǔn)值時存在一定的不確定性.

    基于科、屬的BER預(yù)測效果較好可能是因為:基于科、屬的生物組合方式包含了4種不同分類群的生物毒性數(shù)據(jù),生物分類比較全面,具有一定的代表性和穩(wěn)定性.荷蘭在推導(dǎo)其國家水生生物基準(zhǔn)時,規(guī)定當(dāng)存在4種不同生物分類群的毒性數(shù)據(jù)時,即可推導(dǎo)水生生物基準(zhǔn)值[67].US EPA規(guī)定在推導(dǎo)各州水生生物基準(zhǔn)值時,需要至少4種不同生物分類群的物種毒性數(shù)據(jù)[68].此外,美國ICE技術(shù)也指出在預(yù)測水生生物基準(zhǔn)值時,同科、同屬替代—預(yù)測物種對的預(yù)測效果較好[25].通過比較分析可知,本文基于科、屬的生物組合方式2、4、5的BER技術(shù)對水生生物基準(zhǔn)值的預(yù)測效果較好.而基于特定種的生物組合方式6、7的BER技術(shù)由于物種數(shù)量較少(2種),以致預(yù)測效果不確定性較大.基于不同生物分類間比較的組合方式1的BER技術(shù),由于中美間進行比較的生物分類差別較大,以致預(yù)測效果很差.基于與美國前4種物種同種的BER技術(shù)(組合方式3),由于兩國間特有物種的差異,以致其預(yù)測效果也存在一定的不確定性.

    2.6 不確定性分析

    由于美國至今僅推薦了58種污染物的水生生物基準(zhǔn)值,我國本土生物毒性數(shù)據(jù)也存在缺乏的問題,因此,通過篩選后本文進行了9種污染物的BER技術(shù)的探索與分析,污染物數(shù)量較少可能會帶來水生生物基準(zhǔn)值估算的一些不確定性.采用7種生物組合方式來對中美水生生物基準(zhǔn)預(yù)測與估算進行探索,多樣的生物組合方式可一定程度上減小BER估算預(yù)測的不確定性.

    3 結(jié)論

    3.1 在US EPA推薦的58種污染物的基礎(chǔ)上,結(jié)合“三門八科”生物毒性數(shù)據(jù)需求原則,最終篩選出9個中美共有且本土生物毒性數(shù)據(jù)較豐富的污染物,包括:As(III)、Cr(VI)、Hg、Cu、Zn、Pb、對硫磷、毒死蜱和三丁基錫(TBT);依據(jù)本土生物毒性數(shù)據(jù),推導(dǎo)該9種污染物的本土水生生物急性基準(zhǔn)值分別為:201.72、2.64、0.74、1.32、55.83、92.25、0.12、0.36和0.38 μg/L.

    3.2 提出了7種本土生物組合方式的BER技術(shù),并對基于實測毒性數(shù)據(jù)的基準(zhǔn)值和基于BER技術(shù)所預(yù)測的基準(zhǔn)值進行了差異分析,依據(jù)不同生物組合的預(yù)測效果,初步提出基于科或?qū)俚纳锝M合方式2(與US EPA前4屬同類的我國4屬/美國基準(zhǔn)文件前4屬)、生物組合方式4(我國溞科大型溞+蚤狀溞+鋸頂?shù)皖~溞+模糊網(wǎng)紋溞/美國基準(zhǔn)文件中類似的4個溞科物種)和生物組合方式5(我國指定4屬中各1種生物/美國同類的4屬生物)的BER技術(shù)可較好地對9種污染物的本土水生生物基準(zhǔn)值進行預(yù)測.研究結(jié)果表明在充分地利用現(xiàn)有生物毒性數(shù)據(jù)或進行較少的生物毒性試驗(4種)的基礎(chǔ)上可對我國水生生物基準(zhǔn)值進行預(yù)測.

    參考文獻:

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

    [2] 劉征濤.水環(huán)境質(zhì)量基準(zhǔn)方法與應(yīng)用 [M].北京:科學(xué)出版社, 2012.

    [3] 閆振廣,余若禎,焦聰穎,等.水質(zhì)基準(zhǔn)方法學(xué)中若干關(guān)鍵技術(shù)探討 [J].環(huán)境科學(xué)研究, 2012,25(4):397-403.

    [4] 王曉南,劉征濤,王婉華,等.重金屬鉻(VI)的生態(tài)毒性及其土壤環(huán)境基準(zhǔn) [J].環(huán)境科學(xué), 2014,35(8):3155-3161.

    [5] 劉征濤,王曉南,閆振廣,等.“三門六科”水質(zhì)基準(zhǔn)最少毒性數(shù)據(jù)需求原則 [J].環(huán)境科學(xué)研究, 2012,25(12):1364-1369.

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

    [7] 閆振廣,孟 偉,劉征濤,等.我國淡水生物氨氮基準(zhǔn)研究 [J].環(huán)境科學(xué), 2011,32(6):1564-1570.

    [8] 王曉南,鄭 欣,閆振廣,等.水質(zhì)基準(zhǔn)魚類受試生物篩選 [J].環(huán)境科學(xué)研究, 2014,27(4):341-348.

    [9] 鄭 欣,閆振廣,王曉南,等.水質(zhì)基準(zhǔn)甲殼類受試生物篩選 [J].環(huán)境科學(xué)研究, 2014,27(4):356-364.

    [10] 蔡 靳,閆振廣,何 麗,等.水質(zhì)基準(zhǔn)兩棲類受試生物篩選 [J].環(huán)境科學(xué)研究, 2014,27(4):349-355.

    [11] 王偉莉,閆振廣,劉征濤,等.水質(zhì)基準(zhǔn)本土環(huán)節(jié)動物與水生昆蟲受試生物篩選 [J].環(huán)境科學(xué)研究, 2014,27(4):365-372.

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

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

    [14] 閆振廣,孟 偉,劉征濤,等.我國淡水生物氨氮基準(zhǔn)研究 [J].環(huán)境科學(xué), 2011,32(6):1564-1570.

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

    [16] Jin X, Zha J, Xu Y, et al.Derivation of aquatic predicted no-effect concentration (PNEC) for 2, 4-dichlorophenol: Comparing native species data with non-native species data [J].Chemosphere, 2011,84(10):1506-1511.

    [17] Wang X, Liu Z, Yan Z, et al.Development of aquatic life criteria for triclosan and comparison of the sensitivity between native and non-native species [J].Journal of hazardous materials, 2013,260: 1017-1022.

    [18] Yan Z, Zhang Z, Wang H, et al.Development of aquatic life criteria for nitrobenzene in China [J].Environmental Pollution, 2012,162:86-90.

    [19] Yang S, Yan Z, Xu F, et al.Development of freshwater aquatic life criteria for tetrabromobisphenol A in China [J].Environmental Pollution, 2012,169:59-63.

    [20] US EPA.Guidelines for deriving numerical national water quality criteria for the protection of aquatic organisms and their uses (PB 85-227049) [R].Duluth, Minnesota, Narragansett, Rhode Island; Corvallis, Oregon, USA: Office of Research and Development, Environmental Research Laboratories, U.S.Environmental Protection Agency, 1985.

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

    [22] Davies P, Cook L, Goenarso D.Sublethal responses to pesticides of several species of Australian freshwater fish and crustaceans and rainbow trout [J].Environmental Toxicology and Chemistry, 1994,13(8):1341-1354.

    [23] 中國科學(xué)院.中國動物志 [M].北京:科學(xué)出版社, 1979.

    [24] Wang X, Yan Z, Liu Z, et al.Comparison of species sensitivity distributions for species from China and the USA [J].Environmental Science and Pollution Research, 2014,21(1): 168-176.

    [25] USEPA.Web-based interspecies correlation estimation (Web-ICE) for acute toxicity: User manual [M].Gulf Breeze:USEPA, Office of Research and Development, 2013.

    [26] 閆振廣,孟 偉,劉征濤,等.我國典型流域鎘水質(zhì)基準(zhǔn)研究 [J].環(huán)境科學(xué)研究, 2010,23(10):1221-1228.

    [27] Dyer S, Versteeg D, Belanger S, et al.Comparison of species sensitivity distributions derived from interspecies correlation models to distributions used to derive water quality criteria [J].Environmental Science and Technology, 2008,42(8):3076-3083.

    [28] Dyer S, Versteeg D, Belanger S, et al.Interspecies correlationestimates predict protective environmental concentrations [J].Environmental Science and Technology, 2006,40(9):3102-3111.

    [29] Blum D, Speece R.Determining chemical toxicity to aquatic species [J].Environmental Science and Technology, 1990,24(3): 284-293.

    [30] Feng C, Wu F, Dyer S, et al.Derivation of freshwater quality criteria for zinc using interspecies correlation estimation models to protect aquatic life in China [J].Chemosphere, 2013,90(3): 1177-1183.

    [31] Wu F, Mu Y, Chang H, et al.Predicting water quality criteria for protecting aquatic life from physicochemical properties of metals or metalloids [J].Environmental Science and Technology, 2012,47(1): 446-453.

    [32] Awkerman J, Raimondo S, Barron M.Development of species sensitivity distributions for wildlife using interspecies toxicity correlation models [J].Environmental Science and Technology, 2008,42(9):3447-3452.

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

    [34] 閆振廣,何 麗,高 富,等.鉛水生生物水質(zhì)基準(zhǔn)研究與初步應(yīng)用 [J].環(huán)境科學(xué)與技術(shù), 2014,37(4):31-37.

    [35] US EPA.1995Updates: water quality criteria documents for the protection of aquatic life in ambient water [R].Office of Water, United States Environmental Protection Agency, 1996.

    [36] 馬繼華,蔣耀培,劉 泉,等.“毒死蜱”對南美白對蝦的急性毒性影響 [J].水產(chǎn)科技情報, 2009,(4):192-194.

    [37] 葛士林,曹傳旺,王志英.3種農(nóng)藥對紅裸須搖蚊體內(nèi)蛋白質(zhì)質(zhì)量分?jǐn)?shù)和AChE活性的影響 [J].東北林業(yè)大學(xué)學(xué)報, 2011,39(1): 108-109.

    [38] 譚亞軍.四種殺蟲劑對大型蚤(D.magna)的慢性毒性 [D].浙江大學(xué), 2004.

    [39] 徐 燕,邊文杰,李少南,等.5種殺蟲劑在水-人工底泥系統(tǒng)中對大型溞的急性毒性及其比較 [J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011,30(5): 855-859.

    [40] 吳長興.毒死蜱和氟蟲腈的環(huán)境毒理與風(fēng)險 [D].杭州:浙江大學(xué)環(huán)境與資源學(xué)院, 2010.

    [41] 來有鵬,劉賢金,余向陽,等.農(nóng)藥對河蝦的急性毒性及其谷光甘肽轉(zhuǎn)移酶(GSTs)的影響 [J].農(nóng)藥, 2008,47(11):820-822.

    [42] 丁正峰,薛 暉,王曉豐,等.毒死蜱(CPF)對克氏原螯蝦的急性毒性及組織病理觀察 [J].生態(tài)與農(nóng)村環(huán)境學(xué)報, 2012,28(4): 462-467.

    [43] 李 康.阿特拉津和毒死蜱對中華絨螯蟹(Eriocheir sinensis)的毒性效應(yīng)研究 [D].華東師范大學(xué), 2005.

    [44] 邊文杰,徐燕,李少南,等.毒死蜱與魚肝微粒體P450的相互作用[J].浙江農(nóng)業(yè)學(xué)報, 2011,23(4):776-781.

    [45] 張 莉,顧曉軍,劉 新.魚類肝臟微粒體對于硫代磷酸酯類殺蟲劑活化代謝率的種間差異 [J].生態(tài)毒理學(xué)報, 2010,5(6): 809-816.

    [46] 李常健,駱 鷹,楊錦兀,等.阿維菌素·毒死蜱對草魚的毒性效應(yīng)研究 [J].安徽農(nóng)業(yè)科學(xué), 2011,39(6):3365-3367.

    [47] 丁正峰,薛 暉,王曉豐,等.毒死蜱(CPF)對河川沙塘鱧幼魚的急性毒性 [J].中國水產(chǎn)科學(xué), 2012,19(3):528-535.

    [48] 劉 麗.毒死蜱對鯉魚的免疫毒性研究 [J].河南師范大學(xué), 2013.

    [49] 馬軍國,李效宇.毒死蜱對尖膀胱螺的急性毒性試驗 [J].湖北農(nóng)業(yè)科學(xué), 2012,51(19):4330-4332.

    [50] 馬軍國.白鰱細(xì)胞色素P450(CYP3A137)的克隆及其表達研究[J].河南師范大學(xué), 2013.

    [51] 夏曉華,張林霞,王莉娜,等.毒死蜱對大鱗副泥鰍的急性毒性及遺傳毒性研究 [J].河南農(nóng)業(yè)科學(xué), 2013,1:129-131.

    [52] 王 娜,劉莉莉,孫凱峰,等.有機磷農(nóng)藥的構(gòu)效關(guān)系及其對浮游生物的毒性效應(yīng) [J].生態(tài)環(huán)境學(xué)報, 2012,21(1):118-123.

    [53] 鄭 璐.稻田排水中毒死蜱對鯽魚毒性效應(yīng)的研究 [D].吉林農(nóng)業(yè)大學(xué), 2013.

    [54] 夏錦瑜,王冬蘭,張志勇,等.毒死蜱和丙溴磷對鯽魚的毒性效應(yīng)及其在魚體中的富集 [J].江蘇農(nóng)業(yè)學(xué)報, 2010,26(5):1088-1092.

    [55] 侯方浩,余向陽,趙于丁,等.毒死蜱對錦鯽性腺的影響及其在魚組織中的富集 [J].江蘇農(nóng)業(yè)學(xué)報, 2009,25(1):188-191.

    [56] 尹曉輝,朱國念,鄧建宇,等.毒死蜱對中華蟾蜍胚胎發(fā)育毒性及生態(tài)毒理效應(yīng) [C].第六屆全國環(huán)境化學(xué)大會暨環(huán)境科學(xué)儀器與分析儀器展覽會摘要集, 2011.

    [57] 沈 堅,趙 穎,李少南,等.三種常用農(nóng)藥對環(huán)棱螺,圓田螺和河蜆的急性毒性研究 [J].農(nóng)藥學(xué)學(xué)報, 2013,15(5):559-566.

    [58] 宋志慧,陳天乙,馬維琦.有機錫化合物對蘿卜螺的毒性作用 [J].環(huán)境科學(xué), 1997,4:83-84.

    [59] 王珊珊,馮 流.有機錫化合物的毒性效應(yīng)及其影響因素 [J].安全與環(huán)境學(xué)報, 2005,5(3):12-15.

    [60] 張寶東,張毓琪.氯化三丁基錫對鋰魚的急性毒性 [J].交通環(huán)保:水運版, 1995,16(1):18-19.

    [61] 譚 燕,李遠(yuǎn)友,杜永兵,等.三丁基錫對孔雀魚的毒性效應(yīng)研究[J].環(huán)境化學(xué), 2006,25(3):315-319.

    [62] 宋志慧,陳天乙,劉如冰.三丁基錫對搖蚊幼蟲的毒性作用 [J].環(huán)境科學(xué), 1998,19(2):87-88.

    [63] 王 蕾.三丁基錫對非洲爪蟾和紅鯽的毒性效應(yīng) [J].華東師范大學(xué), 2008.

    [64] 朱春華,薛海波,李廣麗.三丁基錫(TBT)對羅氏沼蝦的毒性效應(yīng)[J].生態(tài)學(xué)雜志, 2011,30(3):527-532.

    [65] 陳天乙,鄭素平.有機錫化合物對水生無脊椎動物的毒性研究[J].環(huán)境化學(xué), 1994,13(3):266-268.

    [66] 陳天乙,譚元生,田 明.有機錫化合物對底棲動物的急性毒性研究 [J].環(huán)境科學(xué), 1994,15(5):63-64.

    [67] RIVM.Guidance document on deriving environmental risk limits in the Netherlands.Report No.601501012 [R].Bilthoven,the Netherlands: National Institute of Public Health and the Environment, 2001.

    [68] USEPA.Water quality standards handbook [M].Washington DC: Office of Water, 1994.

    Study of using biological effect ratio (BER) method to predict the aquatic life criteria in China.


    WANG Xiao-nan1, YAN Zhen-guang1, LIU Zheng-tao1*, WANG Wan-hua1, ZHANG Cong2(1.State Key Laboratory of Environmental Criteria and Risk Assessment, State Environmental Protection Key Laboratory of Ecological Effects and Risk Assessment of Chemicals, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;2.China Offshore Environmental Services Co., Ltd., Tianjin 300452, China).China Environmental Science, 2016,36(1):276~285

    Abstract:In this study, predicted aquatic life criteria using biological effect ratio (BER) method based on the difference of species sensitivity between China and the USA was investigated.First, pollutants with acute toxicity data from 3Phyla and 8Families that both in China and the USA were selected.Second, the different biological groups were developed for BER method based on the sensitivity and representation of species of the two countries.Third, effective BER method was selected based on the comparison of predicted criteria maximum concentration (CMC) and the measured CMC.Results showed that: 9pollutants of As(III), Cr(VI), Hg, Cu, Zn, Pb, parathion, chlorpyrifos and TBT both in China and the USA were selected, and the measured CMCs for protecting the Chinese native aquatic species were derived to be 201.72, 2.64, 0.74, 1.32, 55.83, 92.25, 0.12, 0.36 and 0.38μg/L, respectively.Moreover, comparison of predicted criteria and the measured criteria of 7different biological groups for BER method showed that BER method of biological groups based on the same genera or family could predict the CMCs of 9pollutants in China well.The result of this study could provide useful information for predicting CMC making full use of the existed toxicity data or just carrying out less toxicity test when toxicity data of native species is lacking.

    Key words:biological effect ratio (BER);aquatic life criteria;predicted criteria;China;USA

    中圖分類號:X171.5

    文獻標(biāo)識碼:A

    文章編號:1000-6923(2016)01-0276-10

    收稿日期:2015-05-26

    基金項目:國家水體污染與治理科技重大專項(2012ZX07501-003-06);環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室自由探索基金(SKLECRA2014OFP06)

    作者簡介:王曉南(1986-),男,山西運城人,助理研究員,博士,主要從事生態(tài)毒理學(xué)、環(huán)境基準(zhǔn)與風(fēng)險評估研究.發(fā)表論文20余篇.

    猜你喜歡
    美國中國
    當(dāng)代美國恐怖影片中的哥特元素
    我國警衛(wèi)反恐怖斗爭對策研究
    科技視界(2016年21期)2016-10-17 18:54:07
    美國商業(yè)銀行衍生品業(yè)務(wù)的特點分析
    商(2016年27期)2016-10-17 05:55:39
    淺析當(dāng)代中國社會發(fā)展的代價問題
    科技視界(2016年20期)2016-09-29 14:11:15
    肯尼·格雷特,爵士的“中國”調(diào)子
    美國黑人影片中的文化認(rèn)同焦慮
    英媒:“中國”成美國網(wǎng)絡(luò)威脅敏感詞
    久久九九热精品免费| 欧美乱色亚洲激情| 一进一出抽搐gif免费好疼 | 韩国av一区二区三区四区| 国产一区二区三区综合在线观看| 超色免费av| 久久久国产成人精品二区 | 国产成人精品在线电影| 美女 人体艺术 gogo| 久久人妻av系列| 老司机在亚洲福利影院| 美女大奶头视频| 淫秽高清视频在线观看| 免费不卡黄色视频| 高清欧美精品videossex| 12—13女人毛片做爰片一| 少妇被粗大的猛进出69影院| 亚洲一码二码三码区别大吗| www.精华液| 9191精品国产免费久久| 亚洲全国av大片| 亚洲熟妇中文字幕五十中出 | 日韩免费av在线播放| 欧美激情 高清一区二区三区| 免费搜索国产男女视频| 男女高潮啪啪啪动态图| 精品少妇一区二区三区视频日本电影| 亚洲激情在线av| 久久亚洲真实| 成人影院久久| 精品国产亚洲在线| 精品乱码久久久久久99久播| 亚洲色图 男人天堂 中文字幕| 久久久久精品国产欧美久久久| 手机成人av网站| 美女大奶头视频| 91av网站免费观看| 久久久国产一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 无遮挡黄片免费观看| 久久久水蜜桃国产精品网| 久久精品亚洲av国产电影网| 精品国产乱子伦一区二区三区| 久久精品91蜜桃| 国产主播在线观看一区二区| 十八禁网站免费在线| 一进一出好大好爽视频| 午夜日韩欧美国产| 久久热在线av| netflix在线观看网站| 精品午夜福利视频在线观看一区| 亚洲avbb在线观看| 91老司机精品| 免费在线观看视频国产中文字幕亚洲| 亚洲精品成人av观看孕妇| 在线观看www视频免费| 精品人妻1区二区| 他把我摸到了高潮在线观看| 午夜福利在线免费观看网站| 一个人观看的视频www高清免费观看 | 久久伊人香网站| 精品国内亚洲2022精品成人| 99国产精品一区二区三区| 国产亚洲精品久久久久久毛片| 9热在线视频观看99| 亚洲国产精品合色在线| 久久香蕉精品热| 亚洲av成人不卡在线观看播放网| 国产精品香港三级国产av潘金莲| av在线天堂中文字幕 | 日本一区二区免费在线视频| 一级片'在线观看视频| 夜夜躁狠狠躁天天躁| 欧美 亚洲 国产 日韩一| 亚洲成人久久性| 国产成年人精品一区二区 | 中文亚洲av片在线观看爽| 欧美午夜高清在线| 18禁黄网站禁片午夜丰满| 久久 成人 亚洲| 午夜精品国产一区二区电影| av在线天堂中文字幕 | 精品一区二区三区视频在线观看免费 | 叶爱在线成人免费视频播放| 国产97色在线日韩免费| 又紧又爽又黄一区二区| 99在线人妻在线中文字幕| 欧美色视频一区免费| 夜夜躁狠狠躁天天躁| 999精品在线视频| 久久久水蜜桃国产精品网| 18禁黄网站禁片午夜丰满| 亚洲av成人一区二区三| 在线观看一区二区三区激情| 中亚洲国语对白在线视频| 国产精品一区二区在线不卡| 亚洲中文字幕日韩| 免费高清在线观看日韩| 夜夜夜夜夜久久久久| 91九色精品人成在线观看| 亚洲精品一区av在线观看| 淫秽高清视频在线观看| 国产激情久久老熟女| 欧美日韩av久久| 日本黄色视频三级网站网址| 久久狼人影院| 国产成人精品在线电影| a级毛片在线看网站| 成人黄色视频免费在线看| 亚洲成a人片在线一区二区| 91麻豆av在线| 国产99久久九九免费精品| 日韩大尺度精品在线看网址 | 亚洲成国产人片在线观看| 老司机午夜十八禁免费视频| 人人澡人人妻人| 在线国产一区二区在线| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看| 满18在线观看网站| 99香蕉大伊视频| 久久精品国产99精品国产亚洲性色 | 国产精品二区激情视频| 亚洲人成电影观看| 在线国产一区二区在线| av视频免费观看在线观看| www.www免费av| 国产片内射在线| 日本黄色视频三级网站网址| 黄色视频,在线免费观看| 亚洲第一欧美日韩一区二区三区| 激情视频va一区二区三区| 国产精品一区二区三区四区久久 | 日本黄色日本黄色录像| 少妇的丰满在线观看| 老司机亚洲免费影院| 日韩精品免费视频一区二区三区| 中文欧美无线码| 脱女人内裤的视频| 亚洲第一av免费看| 亚洲一区二区三区欧美精品| 可以在线观看毛片的网站| 丝袜在线中文字幕| 日本a在线网址| 亚洲欧美精品综合久久99| 又黄又爽又免费观看的视频| 亚洲熟妇中文字幕五十中出 | 在线观看一区二区三区| 69精品国产乱码久久久| 美国免费a级毛片| 日韩精品免费视频一区二区三区| 久久久国产成人免费| av国产精品久久久久影院| av天堂在线播放| 亚洲精品中文字幕一二三四区| 欧美成人免费av一区二区三区| 午夜激情av网站| 亚洲第一青青草原| 欧美日韩福利视频一区二区| 色婷婷久久久亚洲欧美| 日本精品一区二区三区蜜桃| av片东京热男人的天堂| 国产视频一区二区在线看| 亚洲精品美女久久久久99蜜臀| 免费在线观看黄色视频的| 热re99久久国产66热| 国产男靠女视频免费网站| 免费女性裸体啪啪无遮挡网站| 人人妻人人澡人人看| netflix在线观看网站| 亚洲国产精品999在线| 久久久久久久久久久久大奶| 日日夜夜操网爽| 麻豆久久精品国产亚洲av | 国产熟女午夜一区二区三区| 久久精品国产99精品国产亚洲性色 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲视频免费观看视频| 日本黄色日本黄色录像| 亚洲黑人精品在线| e午夜精品久久久久久久| 精品一区二区三区av网在线观看| 老司机深夜福利视频在线观看| 一进一出抽搐gif免费好疼 | 精品无人区乱码1区二区| 成人影院久久| www日本在线高清视频| 嫩草影视91久久| 日本精品一区二区三区蜜桃| x7x7x7水蜜桃| 一级a爱视频在线免费观看| 一二三四社区在线视频社区8| 国产精品美女特级片免费视频播放器 | 日韩免费高清中文字幕av| av欧美777| 国产av又大| 亚洲av成人一区二区三| 欧美激情 高清一区二区三区| 国产欧美日韩综合在线一区二区| 女同久久另类99精品国产91| 国产精品av久久久久免费| 校园春色视频在线观看| 国产成人欧美| 人成视频在线观看免费观看| 一级毛片高清免费大全| 宅男免费午夜| 久久草成人影院| 成人特级黄色片久久久久久久| 成人免费观看视频高清| 色综合欧美亚洲国产小说| 亚洲精品成人av观看孕妇| 80岁老熟妇乱子伦牲交| 日本vs欧美在线观看视频| 变态另类成人亚洲欧美熟女 | 男人操女人黄网站| 好看av亚洲va欧美ⅴa在| 亚洲第一av免费看| 久99久视频精品免费| 日日摸夜夜添夜夜添小说| 久久人人精品亚洲av| 最近最新中文字幕大全免费视频| 久久人人爽av亚洲精品天堂| 亚洲国产精品999在线| 欧美精品一区二区免费开放| 757午夜福利合集在线观看| 久久中文看片网| 99热国产这里只有精品6| 不卡av一区二区三区| 久热这里只有精品99| 欧美日韩一级在线毛片| 老司机亚洲免费影院| 美女午夜性视频免费| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| av欧美777| 欧美精品一区二区免费开放| 女性被躁到高潮视频| 久久中文看片网| 国产成人一区二区三区免费视频网站| 亚洲色图av天堂| 亚洲精品美女久久av网站| 大码成人一级视频| 一进一出抽搐gif免费好疼 | 91字幕亚洲| 12—13女人毛片做爰片一| 99国产精品一区二区蜜桃av| 精品国产乱子伦一区二区三区| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美日韩瑟瑟在线播放| 国内毛片毛片毛片毛片毛片| 夜夜躁狠狠躁天天躁| 黑人巨大精品欧美一区二区蜜桃| 视频区欧美日本亚洲| 日韩大码丰满熟妇| 级片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲片人在线观看| 男人的好看免费观看在线视频 | 久久热在线av| 日韩大尺度精品在线看网址 | 亚洲专区字幕在线| 色在线成人网| 国产av精品麻豆| av片东京热男人的天堂| 一本大道久久a久久精品| 99久久人妻综合| 欧美日韩中文字幕国产精品一区二区三区 | 精品一区二区三区视频在线观看免费 | 天堂√8在线中文| 香蕉国产在线看| 亚洲在线自拍视频| 久久亚洲精品不卡| 美女 人体艺术 gogo| 黄色成人免费大全| 国产伦人伦偷精品视频| 免费观看精品视频网站| 在线永久观看黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文日韩欧美视频| 国产无遮挡羞羞视频在线观看| 女生性感内裤真人,穿戴方法视频| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产区一区二| 这个男人来自地球电影免费观看| 色老头精品视频在线观看| 久久狼人影院| 18禁黄网站禁片午夜丰满| 少妇 在线观看| 久久人妻福利社区极品人妻图片| 久久国产亚洲av麻豆专区| 国产三级在线视频| 在线观看午夜福利视频| 又黄又爽又免费观看的视频| 女性生殖器流出的白浆| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 天堂动漫精品| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 亚洲免费av在线视频| 999久久久国产精品视频| 18禁裸乳无遮挡免费网站照片 | 精品福利永久在线观看| 女性被躁到高潮视频| 大陆偷拍与自拍| 日日夜夜操网爽| 免费在线观看完整版高清| 美女扒开内裤让男人捅视频| 国产xxxxx性猛交| 国产成人一区二区三区免费视频网站| 精品电影一区二区在线| 波多野结衣一区麻豆| 麻豆成人av在线观看| 99国产综合亚洲精品| 日日爽夜夜爽网站| 亚洲熟妇中文字幕五十中出 | 女性生殖器流出的白浆| 国产精品久久电影中文字幕| 丝袜在线中文字幕| 在线免费观看的www视频| 天天躁夜夜躁狠狠躁躁| 欧美激情 高清一区二区三区| 欧美色视频一区免费| 看免费av毛片| 久久婷婷成人综合色麻豆| 久久精品影院6| 桃红色精品国产亚洲av| 久久中文字幕人妻熟女| 国产欧美日韩一区二区精品| 亚洲情色 制服丝袜| 欧美日韩精品网址| 18禁国产床啪视频网站| 91国产中文字幕| 中文字幕高清在线视频| 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 热re99久久国产66热| 日韩精品免费视频一区二区三区| 成年人黄色毛片网站| 久久精品亚洲精品国产色婷小说| 男女做爰动态图高潮gif福利片 | 国产免费av片在线观看野外av| 香蕉国产在线看| 欧美日韩国产mv在线观看视频| 一a级毛片在线观看| 欧美日韩视频精品一区| 9191精品国产免费久久| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 久久久久国产精品人妻aⅴ院| 国产有黄有色有爽视频| 亚洲一区二区三区不卡视频| 日韩欧美免费精品| 搡老熟女国产l中国老女人| 淫秽高清视频在线观看| netflix在线观看网站| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| 国产一区二区三区在线臀色熟女 | 亚洲精品中文字幕一二三四区| 久久精品国产综合久久久| 亚洲精品国产精品久久久不卡| 国产在线精品亚洲第一网站| 国产精品二区激情视频| 日日夜夜操网爽| 亚洲熟妇熟女久久| 看片在线看免费视频| 久久精品影院6| 亚洲欧美日韩高清在线视频| www.999成人在线观看| 亚洲情色 制服丝袜| 日韩欧美三级三区| 制服人妻中文乱码| 在线免费观看的www视频| 欧美大码av| 亚洲激情在线av| 色尼玛亚洲综合影院| 在线观看一区二区三区| 日韩欧美在线二视频| 一a级毛片在线观看| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区视频在线观看| 天堂影院成人在线观看| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 国产精品九九99| 国产99久久九九免费精品| 亚洲成av片中文字幕在线观看| 久久人人爽av亚洲精品天堂| 亚洲成av片中文字幕在线观看| 最新美女视频免费是黄的| 精品一品国产午夜福利视频| 国产高清videossex| 一区福利在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人精品巨大| 黄色丝袜av网址大全| 国产蜜桃级精品一区二区三区| 欧美另类亚洲清纯唯美| 免费av毛片视频| 日韩欧美在线二视频| 精品一品国产午夜福利视频| 色尼玛亚洲综合影院| 国产欧美日韩一区二区三| 欧美大码av| a级片在线免费高清观看视频| 搡老乐熟女国产| 91麻豆精品激情在线观看国产 | 国产一区二区在线av高清观看| 国产成人免费无遮挡视频| 天堂√8在线中文| 国产精品电影一区二区三区| av有码第一页| 亚洲中文日韩欧美视频| 91av网站免费观看| a级毛片黄视频| 日韩人妻精品一区2区三区| 亚洲国产精品合色在线| 国产精品 国内视频| 曰老女人黄片| 色精品久久人妻99蜜桃| 久久久久久免费高清国产稀缺| 国产精品成人在线| 欧美日韩亚洲高清精品| 日本一区二区免费在线视频| 神马国产精品三级电影在线观看 | 亚洲人成网站在线播放欧美日韩| 亚洲中文av在线| 精品高清国产在线一区| 国产精品av久久久久免费| 国产伦一二天堂av在线观看| 中文字幕人妻丝袜制服| 久久精品亚洲熟妇少妇任你| 久久99一区二区三区| 亚洲av五月六月丁香网| 超色免费av| 亚洲精品中文字幕在线视频| 久久香蕉激情| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 美女大奶头视频| 法律面前人人平等表现在哪些方面| 91麻豆精品激情在线观看国产 | av免费在线观看网站| 免费人成视频x8x8入口观看| 一区二区三区激情视频| av福利片在线| 村上凉子中文字幕在线| 老司机午夜福利在线观看视频| 国产一区在线观看成人免费| 精品久久蜜臀av无| 精品国产美女av久久久久小说| 国产aⅴ精品一区二区三区波| 欧美黑人欧美精品刺激| 真人做人爱边吃奶动态| 香蕉丝袜av| 欧美国产精品va在线观看不卡| 亚洲国产精品一区二区三区在线| 国产乱人伦免费视频| 丝袜人妻中文字幕| ponron亚洲| 女警被强在线播放| 国产野战对白在线观看| 欧美另类亚洲清纯唯美| 88av欧美| 真人做人爱边吃奶动态| 亚洲av第一区精品v没综合| 欧美精品亚洲一区二区| 大香蕉久久成人网| 好男人电影高清在线观看| 久久久国产成人精品二区 | 国产高清国产精品国产三级| 亚洲av美国av| 精品卡一卡二卡四卡免费| 日日干狠狠操夜夜爽| 久热爱精品视频在线9| 中文欧美无线码| 99久久久亚洲精品蜜臀av| 精品国产一区二区久久| 亚洲欧美日韩无卡精品| 熟女少妇亚洲综合色aaa.| 久久精品国产清高在天天线| 国产99久久九九免费精品| 99国产综合亚洲精品| 九色亚洲精品在线播放| 香蕉国产在线看| 99国产极品粉嫩在线观看| 在线av久久热| 在线观看免费视频网站a站| 久久久久久久久免费视频了| 亚洲av电影在线进入| 免费日韩欧美在线观看| 麻豆av在线久日| 国产一区二区三区视频了| 亚洲色图 男人天堂 中文字幕| 在线观看免费视频日本深夜| 久久中文字幕人妻熟女| 亚洲国产欧美一区二区综合| 精品卡一卡二卡四卡免费| 亚洲中文av在线| 久久影院123| 国产精品影院久久| 精品国产一区二区久久| 老熟妇仑乱视频hdxx| 制服诱惑二区| 久久精品91蜜桃| 91国产中文字幕| 在线观看www视频免费| 满18在线观看网站| 精品国内亚洲2022精品成人| 久久天堂一区二区三区四区| 99久久综合精品五月天人人| 18禁观看日本| 真人一进一出gif抽搐免费| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品一区二区三区在线| 国产熟女xx| 成人免费观看视频高清| 亚洲色图综合在线观看| 亚洲成av片中文字幕在线观看| 最近最新中文字幕大全电影3 | 国产精品久久久久久人妻精品电影| 午夜91福利影院| 国产精品九九99| 中文字幕最新亚洲高清| 欧美日韩亚洲高清精品| 母亲3免费完整高清在线观看| 国产视频一区二区在线看| 久久这里只有精品19| 天天影视国产精品| 欧美久久黑人一区二区| 久久精品人人爽人人爽视色| 手机成人av网站| 国产色视频综合| 亚洲伊人色综图| √禁漫天堂资源中文www| 99国产精品免费福利视频| 精品久久久精品久久久| 免费看十八禁软件| www.自偷自拍.com| 亚洲男人的天堂狠狠| 亚洲成av片中文字幕在线观看| 又黄又粗又硬又大视频| 国产成人影院久久av| 日韩免费av在线播放| 久久精品国产99精品国产亚洲性色 | 看黄色毛片网站| 人成视频在线观看免费观看| 黄色 视频免费看| 女人被狂操c到高潮| 88av欧美| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| av中文乱码字幕在线| 男女高潮啪啪啪动态图| 19禁男女啪啪无遮挡网站| 久久精品亚洲熟妇少妇任你| 欧美精品一区二区免费开放| 日韩成人在线观看一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲av电影在线进入| 97人妻天天添夜夜摸| 操出白浆在线播放| av在线播放免费不卡| 欧美最黄视频在线播放免费 | 啦啦啦在线免费观看视频4| 久久精品国产亚洲av高清一级| 热re99久久精品国产66热6| 日韩一卡2卡3卡4卡2021年| 国产精品久久久人人做人人爽| 精品国产国语对白av| 又大又爽又粗| 嫩草影视91久久| 日本三级黄在线观看| 啦啦啦 在线观看视频| 在线十欧美十亚洲十日本专区| 如日韩欧美国产精品一区二区三区| 亚洲自拍偷在线| 久久国产精品人妻蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美激情在线| 最新在线观看一区二区三区| 他把我摸到了高潮在线观看| 亚洲午夜理论影院| 日韩大码丰满熟妇| 乱人伦中国视频| 国产高清国产精品国产三级| 丝袜美腿诱惑在线| 国产亚洲精品久久久久久毛片| 99久久精品国产亚洲精品| 1024香蕉在线观看| 老司机深夜福利视频在线观看| 欧美色视频一区免费| 99re在线观看精品视频| 国产激情欧美一区二区| 精品福利永久在线观看| 男人舔女人的私密视频| √禁漫天堂资源中文www| 黑丝袜美女国产一区| 一边摸一边做爽爽视频免费| 国产亚洲欧美98| 亚洲成人国产一区在线观看| 国产伦一二天堂av在线观看| 国产色视频综合| 美女大奶头视频| 久久 成人 亚洲| 久久精品成人免费网站| 亚洲伊人色综图| 国内久久婷婷六月综合欲色啪| 亚洲狠狠婷婷综合久久图片| 亚洲熟女毛片儿| 女人被狂操c到高潮| 亚洲精品粉嫩美女一区| 午夜两性在线视频|