天津地區(qū)地表水中多環(huán)芳烴的生態(tài)風(fēng)險(xiǎn)

天津地區(qū)地表水中多環(huán)芳烴的生態(tài)風(fēng)險(xiǎn)基金項(xiàng)目：北京大學(xué)政基金，國家自然科學(xué)基金（40031010，40024101）資助作者簡介：石璇（1982—），女，本科生；*通訊聯(lián)系人石璇楊宇徐福留劉文新陶澍（北京大學(xué)環(huán)境學(xué)院,地表過程分析與模擬教育部重點(diǎn)實(shí)驗(yàn)室,北京100871）摘要:以天津地區(qū)46個(gè)水樣中8種多環(huán)芳烴的監(jiān)測濃度及其對(duì)6至38種水生生物的LC50為基礎(chǔ)資料,用重疊面積和聯(lián)合概率曲線兩種概率風(fēng)險(xiǎn)評(píng)價(jià)的方法分析這些多環(huán)芳烴的生態(tài)風(fēng)險(xiǎn)。結(jié)果表明,所研究的化合物中蒽的風(fēng)險(xiǎn)最大。進(jìn)一步,基于“等效系數(shù)”概念,計(jì)算各樣點(diǎn)多環(huán)芳烴的等效總濃度,用概率風(fēng)險(xiǎn)評(píng)價(jià)方法分析了8種多環(huán)芳烴聯(lián)合作用的總生態(tài)風(fēng)險(xiǎn),并與風(fēng)險(xiǎn)商結(jié)果對(duì)比。結(jié)果顯示,總生態(tài)風(fēng)險(xiǎn)顯著高于任一多環(huán)芳烴的單獨(dú)作用。計(jì)算采用的兩種劑量－效應(yīng)關(guān)系的不同假設(shè)對(duì)最終結(jié)果影響不大。關(guān)鍵詞:多環(huán)芳烴；生態(tài)風(fēng)險(xiǎn)；地表水；概率風(fēng)險(xiǎn)評(píng)價(jià)；水生生物EcologicalriskassessmentofpolycyclicaromatichydrocarbonsinsurfacewaterfromTianjinShiXuan,YangYu,XuFuliu,LiuWenxin,TaoShu(MOELaboratoryforEarthSurfaceProcessesAnalysisandSimulation,CollegeofEnvironmentalSciences,PekingUniversity,Beijing100871)Abstract:Basedontheobservedconcentrationsof46samplesofsurfacewaterfromTianjinandLC50valuesofanumberofaquaticorganism,probabilisticriskapproachwasappliedtoaddressecologicalriskof8PAHstolocalaquaticecosystem.Bothjointprobabilitycurveandoverlappingareaofconcentrationandtoxicitydistributionwerecalculatedforthispurpose.TheresultsindicatedthattheecologicalriskoftheindividualPAHswassmallwithanthracenebeingthemosttoxiccompound.Thejointprobabilityriskofthe8PAHswascalculatedbasedonthesumofequivalentconcentrationsderivedfromcalculatedequivalencefactorsandcomparedwithconventionalapproachofhazardquotient.Thejointriskwasmuchhigherthanthatoftheindividualcompounds.Thecalculatedriskwasslightlyinfluencedbythetypeofdose-responserelationships.Keywords:PAHs;ecologicalrisk;surfacewater;probabilityriskassessment,aquaticorganism多環(huán)芳烴對(duì)生物體的危害已得到廣泛證實(shí),其毒性多在光誘導(dǎo)條件下發(fā)生[1,2]。雖然火山活動(dòng)、火災(zāi)和生物合成等自然過程也會(huì)產(chǎn)生多環(huán)芳烴,但化石燃料燃燒無疑是環(huán)境中多環(huán)芳烴的最重要來源[1]。由于大量燃煤和其它工業(yè)污染,天津地區(qū)地表水體遭受多環(huán)芳烴嚴(yán)重污染[3],這些污染必然導(dǎo)致對(duì)生態(tài)系統(tǒng)的威脅。生態(tài)風(fēng)險(xiǎn)分析側(cè)重生態(tài)系統(tǒng)水平的污染效應(yīng),是定量研究有毒污染物生態(tài)危害的重要手段[4]。目前用于風(fēng)險(xiǎn)分析的方法很多,既有簡單的風(fēng)險(xiǎn)系數(shù)的計(jì)算[5],也有基于概率意義的風(fēng)險(xiǎn)分析[6]。概率風(fēng)險(xiǎn)分析中最常見的方法是將表征化合物暴露濃度和毒性參數(shù)的概率密度曲線置于同一座標(biāo)系下,計(jì)算其重疊部分面積,據(jù)此表征生物受不利影響的概率[7]。聯(lián)合概率曲線則以毒性響應(yīng)累積概率（橫軸）和暴露濃度超過相應(yīng)影響邊界濃度的概率（縱軸）作圖表征特定化合物的生態(tài)風(fēng)險(xiǎn)[8]。這些方法被廣泛應(yīng)用于單一污染物的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。對(duì)于多種污染物共同作用導(dǎo)致的生態(tài)危害,大多采用推廣風(fēng)險(xiǎn)系數(shù)的方法[9]。推廣的處理方法基于如下假設(shè):即環(huán)境中存在多種污染物時(shí),如果各種污染物致毒機(jī)理相同,其總效果表現(xiàn)為單一污染物獨(dú)立作用效應(yīng)的簡單疊加[10],這種疊加可以通過風(fēng)險(xiǎn)系數(shù)直接相加[911]或者對(duì)等效濃度進(jìn)行簡單加和[1213]的方法實(shí)現(xiàn)。對(duì)于作用機(jī)制不同,效應(yīng)果無法疊加的污染物,往往權(quán)宜地取單種風(fēng)險(xiǎn)系數(shù)的最大值表現(xiàn)他們的共同作用[914]。與利用簡單閾值或指數(shù)的風(fēng)險(xiǎn)分析相比,概率風(fēng)險(xiǎn)分析更接近真實(shí)情況。本研究的目的在于根據(jù)濃度－響應(yīng)定量關(guān)系的基本形式,推導(dǎo)出基于等效系數(shù)概念的對(duì)多種污染物生態(tài)危害的概率評(píng)價(jià)方法,并應(yīng)用該方法評(píng)價(jià)天津地區(qū)地表水中8種多環(huán)芳烴對(duì)水生生態(tài)系統(tǒng)的危害。1研究方法1.1數(shù)據(jù)獲取在天津地區(qū)主要河流46個(gè)樣點(diǎn)采集了水樣,分析了其中16種優(yōu)控多環(huán)芳烴含量,具體樣點(diǎn)及分析方法參見相關(guān)文獻(xiàn)[4]。根據(jù)毒性數(shù)據(jù)的可獲取性選擇了8種代表性化合物進(jìn)行分析。表1列舉了這8種多環(huán)芳烴在天津地表水體中含量的統(tǒng)計(jì)量。表1天津地表水中8種多環(huán)芳烴含量(g/l)Table1Concentrationsofthe8PAHsinsurfacewaterfromTianjin(g/l)多環(huán)芳烴蒽苯并(a)芘熒蒽芘芴菲萘二氫苊樣本量4646464646464646算術(shù)均值0.0059510.0016680.0068860.0089080.023430.035180.11220.00191標(biāo)準(zhǔn)差0.007020.0018220.0068990.0095880.059270.037740.18660.002434從美國環(huán)保署毒性數(shù)據(jù)庫（/ecotox）和Verschueren[15]中收集到蒽、苯并(a)芘、熒蒽、芘、芴、菲、萘、二氫苊等8種多環(huán)芳烴對(duì)水生生物的急性毒性數(shù)據(jù)。生物種包括綠藻（Selenastrumcapricornutum）、硅藻（Skeletonemacostatum）、五帶淡水庫蚊（Culexquinquefasciatus）、埃及伊蚊（Aedesaegypti）、水蚤（Daphniamagna）、紅鱸（Cyprinodonvariegatus）、鯰魚（Ictaluruspunctatus）、鱒魚（Oncorhynchusmykiss）、黑頭呆魚（Pimephalespromelas）、鉤蝦（Gammarusannulatus）、豹蛙（Ranapipiens）、蠑螈（Spottedsalamander）、蛤（Mulinialateralis）、淡水蝦（Palaemonetes）、蝸牛（Physaheterostropha）等生物。毒性數(shù)據(jù)主要取24~96h的LC50。由于所研究的主要多環(huán)芳烴對(duì)生物的致毒機(jī)理相似,同時(shí)作用效應(yīng)疊加[211],故有可能利用等效系數(shù)表征它們的聯(lián)合作用。表28種多環(huán)芳烴對(duì)水生生物急性毒性數(shù)據(jù)統(tǒng)計(jì)量(g/l)Table2Statisticsofacutetoxicitydataofthe8PAHstoaquaticorganism(g/l)多環(huán)芳烴蒽苯并(a)芘熒蒽芘芴菲萘二氫苊樣本量13113868183114(LC50)min0.01626.646340110240(LC50)max3601200045000200010000042402200001605001.2概率風(fēng)險(xiǎn)分析方法對(duì)于單一污染物,根據(jù)上述資料分別導(dǎo)出區(qū)域暴露濃度和生物耐受水平兩個(gè)概率密度函數(shù),用Matlab計(jì)算兩函數(shù)曲線重疊部分的面積,并求得聯(lián)合概率曲線。對(duì)于多種污染物,采用等效系數(shù)概念[1213],根據(jù)不同污染物導(dǎo)致相同生物死亡率的濃度關(guān)系,將其它化合物的暴露濃度折算成苯并(a)芘等效濃度,加和后按照單一污染物風(fēng)險(xiǎn)分析的方法進(jìn)行評(píng)價(jià)。由于文獻(xiàn)中報(bào)道毒性參數(shù)多為半致死濃度,而多環(huán)芳烴對(duì)水生生物毒性的劑量－響應(yīng)關(guān)系可能為“S”型或指數(shù)型。為對(duì)等效濃度進(jìn)行疊加,有必要根據(jù)響應(yīng)曲線類型反推濃度與LC50間定量關(guān)系。對(duì)于“S”型響應(yīng),按照概率單位的方法可以轉(zhuǎn)換為直線[131617],從而得出以下方程: (1)式中C表示特定多環(huán)芳烴的實(shí)際濃度,C等效代表與此濃度該化合物毒性相當(dāng)?shù)谋讲?a)芘濃度（即等效濃度）,LC50e和LC50i分別表示苯并(a)芘和該化合物對(duì)水生生物的半致死濃度（均值）。對(duì)于簡單指數(shù)響應(yīng)關(guān)系[18],濃度與半致死濃度關(guān)系為: (2)式中符號(hào)與式(1)相同。在這種假設(shè)下,各種多環(huán)芳烴有確定的等效系數(shù),即C等效與C的比值（表三）。表38種多環(huán)芳烴的等效系數(shù)Table3Equivalencefactorofthe8PAHs多環(huán)芳烴蒽芘熒蒽苯并(a)芘菲芴二氫苊萘等效系數(shù)18.983.03110.34等效系數(shù)可以在一定程度上反映各種多環(huán)芳烴的毒性差別,等效系數(shù)越大,毒性越強(qiáng)。所研究的8種多環(huán)芳烴中蒽對(duì)水生生物的毒性最強(qiáng),萘的毒性最弱,苯并(a)芘的毒性居中。在根據(jù)兩種不同方式進(jìn)行上述轉(zhuǎn)換后,有可能根據(jù)8種多環(huán)芳烴實(shí)測濃度求出等效濃度。等效濃度將多種多環(huán)芳烴的危害歸于統(tǒng)一尺度下,可在加和得到總有效濃度后進(jìn)行多污染物總概率風(fēng)險(xiǎn)分析[1314]。同時(shí)可與基于總有效濃度的風(fēng)險(xiǎn)商（等效濃度和苯并(a)芘對(duì)生物L(fēng)C50平均值之比）進(jìn)行對(duì)比。2結(jié)果與討論2.1天津地表水中多環(huán)芳烴含量和水生生物對(duì)多環(huán)芳烴耐受水平的分布特征采用偏度-峰度系數(shù)對(duì)本研究收集的水相多環(huán)芳烴濃度和多環(huán)芳烴對(duì)水生生物毒性的半致死濃度（LC50）原始數(shù)據(jù)及其對(duì)數(shù)變換數(shù)據(jù)進(jìn)行正態(tài)檢驗(yàn)。結(jié)果表明,天津地表水中8種多環(huán)芳烴含量及大多數(shù)半致死濃度為對(duì)數(shù)正態(tài)分布。由于部分毒性參數(shù)數(shù)據(jù)量很?。ū?）,故采用對(duì)數(shù)正態(tài)分布表征全部數(shù)據(jù)。由此得到的概率分布的一般形式為: (3)式中f(C)為概率密度函數(shù),a和b為描述分布函數(shù)的特征參數(shù),C和分別代表多環(huán)芳烴濃度和取對(duì)數(shù)后的均值。有關(guān)具體參數(shù)列于表4。表4天津地區(qū)地表水中多環(huán)芳烴暴露濃度和毒性數(shù)據(jù)的分布參數(shù)(g/l）Table4DistributionparametersforexposureconcentrationsandtoxicitydataofPAHsinsurfacewaterfromTianjin(g/l)多環(huán)芳烴菲芘熒蒽蒽萘二氫苊苯并(a)芘芴濃度均值-3.56-4.74-4.95-5.28-3.05-5.13-6.62-4.84濃度標(biāo)準(zhǔn)差0.800.850.790.951.251.541.111.31a（濃度）0.500.470.510.420.320.260.360.30b（濃度）-0.78-0.69-0.81-0.55-0.32-0.21-0.40-0.29毒性均值6.294.095.202.268.507.085.206.58毒性標(biāo)準(zhǔn)差1.421.892.502.821.190.653.322.33a（毒性）0.240.330.620.120.17b（毒性）-0.25-0.14-0.08-0.06-0.35-1.20-0.05-0.09從表中數(shù)據(jù)可以看出,天津地表水中萘的濃度最高,苯并(a)芘濃度最低；生物對(duì)蒽最敏感,對(duì)萘耐受性最強(qiáng)。就暴露濃度和半致死濃度差別而言,半致死濃度總體上遠(yuǎn)高于環(huán)境中的暴露濃度,其中蒽的差別最小,二氫苊的差別最大。盡管如此,這些污染物對(duì)水生生物仍然具有概率意義上的毒性。2.2單一多環(huán)芳烴的生態(tài)風(fēng)險(xiǎn)分別對(duì)8種多環(huán)芳烴進(jìn)行獨(dú)立概率風(fēng)險(xiǎn)分析。將暴露濃度和毒性數(shù)據(jù)分布置于同一座標(biāo)系下,計(jì)算分布的重疊面積,據(jù)此考察生物受影響比例。表5給出了8種多環(huán)芳烴濃度分布和毒性分布的重疊面積。表5天津地區(qū)地表水中多環(huán)芳烴濃度分布與對(duì)水生生物毒性分布的重疊面積Table5OverlappingareaofthedistributioncurvesbetweenthetoxicitiestoaquaticorganismandexposureconcentrationsofPAHsinsurfacewaterfromTianjin多環(huán)芳烴蒽苯并(a)芘熒蒽芘芴菲萘二氫苊重疊面積7.11026.61033.41032.41031.61031.81055.31064.1108表中化合物由左至右按照風(fēng)險(xiǎn)大小順次排列。可以看出,所研究的8種多環(huán)芳烴造成的風(fēng)險(xiǎn)都比較小。影響相對(duì)較大的為暴露濃度與半致死濃度差別最小的蒽和苯并(a)芘。以受影響生物比例為橫標(biāo),暴露濃度超過響應(yīng)邊界的概率為縱標(biāo),得到的聯(lián)合概率曲線可以更直觀地反映各化合物毒性[7]。圖1即為蒽、芴和萘的聯(lián)合概率曲線。圖1 天津地區(qū)地表水中蒽、芴和萘對(duì)水生生物毒性的聯(lián)合概率曲線Figure1Jointprobabilitycurvesoftoxicityofanthracene,fluorene,naphthalenetoaquaticorganisminTianjin聯(lián)合概率曲線的位置反映了污染物生態(tài)風(fēng)險(xiǎn)的大小。曲線越靠近坐標(biāo)軸,風(fēng)險(xiǎn)越小[78]。其反映的毒性風(fēng)險(xiǎn)順序與表5列舉結(jié)果一致。從三種代表性化合物的聯(lián)合概率曲線可以看出,所研究多環(huán)芳烴風(fēng)險(xiǎn)差別很大（坐標(biāo)橫軸尺度不同）。盡管萘的濃度很高,但由于其毒性較弱,毒性風(fēng)險(xiǎn)低于多數(shù)其它化合物。蒽的相對(duì)高風(fēng)險(xiǎn)主要與其相對(duì)較強(qiáng)的毒性有關(guān),其它毒性較強(qiáng)的化合物由于暴露濃度很小而造成的風(fēng)險(xiǎn)較低。2.3多環(huán)芳烴的總生態(tài)風(fēng)險(xiǎn)根據(jù)式（1）和（2）,在“S”響應(yīng)和指數(shù)響應(yīng)兩種假設(shè)下將8種多環(huán)芳烴的實(shí)測濃度轉(zhuǎn)換為等效濃度后求和,并據(jù)此進(jìn)行總生態(tài)風(fēng)險(xiǎn)分析。由于缺少關(guān)于濃度與毒性響應(yīng)確切關(guān)系的信息,分別在S型響應(yīng)和指數(shù)響應(yīng)假設(shè)條件下進(jìn)行分析。兩種響應(yīng)假設(shè)條件下得到的等效濃度（C等）和均服從對(duì)數(shù)正態(tài)分布,求得的概率密度函數(shù)分別為: (4)和 (5)這些分布與苯并(a)芘對(duì)水生生物L(fēng)C50分布的重疊面積分別為0.17（S響應(yīng)假設(shè)）和0.14（指數(shù)響應(yīng)假設(shè)）,均高于所有8種多環(huán)芳烴獨(dú)立作用下的重疊面積甚至8種化合物重疊面積之和（表5）,清晰地表現(xiàn)出聯(lián)合毒效應(yīng)。兩者的聯(lián)合概率曲線在圖2中給出。為便于比較,圖中還給出了單一風(fēng)險(xiǎn)最高的化合物（蒽）的獨(dú)立作用下的聯(lián)合概率曲線（圖中虛線）。顯然,無論采用哪種濃度－毒性響應(yīng)假設(shè),疊加影響顯著高于危害最大化合物的獨(dú)立作用,且S響應(yīng)假設(shè)下的風(fēng)險(xiǎn)略高于指數(shù)響應(yīng)假設(shè)的計(jì)算結(jié)果。這與重疊面積計(jì)算結(jié)果一致。為與常用非概率風(fēng)險(xiǎn)分析方法比較,計(jì)算了基于總有效濃度的各樣點(diǎn)風(fēng)險(xiǎn)商。兩種假設(shè)條件下的平均風(fēng)險(xiǎn)商分別為0.00097和0.0016。圖3給出了這些風(fēng)險(xiǎn)商的統(tǒng)計(jì)分布。圖2 天津地區(qū)地表水中8種多環(huán)芳烴等效毒性的聯(lián)合概率曲線Figure2Jointprobabilitycurveofequivalenttoxicityof8PAHsinsurfacewaterfromTianjin圖3 基于天津地區(qū)地表水中8種多環(huán)芳烴總等效濃度的風(fēng)險(xiǎn)商（經(jīng)對(duì)數(shù)變換）的統(tǒng)計(jì)分布Figure3Distributionsofhazardquotientsbasedontotalequivalentconcentrationsof8PAHsinsurfacewaterfromTianjin經(jīng)對(duì)數(shù)變換后的風(fēng)險(xiǎn)商服從正態(tài)分布。計(jì)算風(fēng)險(xiǎn)系數(shù)幾何均值分別為0.0055和0.0057。它們反映了污染物對(duì)生物危害的程度。暴露濃度對(duì)LC50的倍數(shù)越高,生物受危害越嚴(yán)重。計(jì)算表明,風(fēng)險(xiǎn)系數(shù)集中分布于數(shù)值很小的區(qū)間上,說明生物受危害較輕。與之不同,概率風(fēng)險(xiǎn)評(píng)價(jià)反映了受危害生物所占比例。重疊面積越大,聯(lián)合概率曲線離坐標(biāo)軸越遠(yuǎn),生物受害的比例越大。兩種假設(shè)下重疊面積都比較小,聯(lián)合概率曲線也接近坐標(biāo)軸。兩種方法從不同的角度說明天津地區(qū)地表水中多環(huán)芳烴對(duì)水生生物的危害較小。在方法學(xué)上,針對(duì)多種多環(huán)芳烴的概率風(fēng)險(xiǎn)評(píng)價(jià)忽略了各種多環(huán)芳烴對(duì)不同生物毒性強(qiáng)弱關(guān)系的差別；而風(fēng)險(xiǎn)商則忽略了不同生物耐受性的差別。不同折算方法對(duì)總等效濃度分布和各種風(fēng)險(xiǎn)分析結(jié)果的影響都不大。將兩種假設(shè)下總等效濃度的分布置于同一座標(biāo)系下得到的重疊面積為0.90,不重疊部分的面積顯著大于概率風(fēng)險(xiǎn)分析中重疊面積的差別（0.17-0.14=0.03）,這也說明此種概率風(fēng)險(xiǎn)分析方法對(duì)暴露濃度分布的變化不太敏感?；谥笖?shù)響應(yīng)假設(shè)下的總風(fēng)險(xiǎn)系數(shù),相當(dāng)于單種污染物風(fēng)險(xiǎn)系數(shù)的加和,這種方法的有效性與其他研究的結(jié)論一致[911]。3結(jié)論天津地表水中主要多環(huán)芳烴中蒽對(duì)水生生物毒性首屈一指。8種主要多環(huán)芳烴的聯(lián)合效果顯著高于任何單一化合物。不同劑量－響應(yīng)假設(shè)得到的結(jié)果差異不大。概率風(fēng)險(xiǎn)評(píng)價(jià)和風(fēng)險(xiǎn)商的計(jì)算從不同的角度反映了污染物的生態(tài)風(fēng)險(xiǎn)。參考文獻(xiàn)[1] 趙云英,馬永安.天然環(huán)境中多環(huán)芳烴的遷移轉(zhuǎn)化及其對(duì)生態(tài)環(huán)境的影響.海洋環(huán)境科學(xué),1998,17:68~72.[2] Erickson,R.J.;Ankley,G.T.;DeFoe,D.L.;Kosian,P.A.;Makynen,E.A.1999AdditivetoxicityofbinarymixturesofphototoxicpolycyclicaromatichydrocarbonstotheOligochaeteLumbriculusvariegatesToxicol.Appl.Pharmacol.154,1999,97-105.[3] 何新春,天津地區(qū)河水有機(jī)氯農(nóng)藥殘留與分布,碩士學(xué)位論文,北京大學(xué)環(huán)境學(xué)院,北京,2003.[4] Solomon,K.R.,Sibley,P.Newconceptsinecologicalriskassessment:wheredowegofromhere?Mar.Pollut.Bull.,2002,44:279~285.[5] Staples,C.A.,Davis,J.W..Anexaminationofthephysicalproperties,fate,ecotoxicityandpotentialenvironmentalrisksforaseriesofpropyleneglycolethers.Chemosphere,2002,49:61~73.[6]Duvall,S.E.,Barron,M.G..AscreeninglevelprobabilisticriskassessmentofmercuryinFloridaEvergladesfoodwebs,Ecotoxicol.Environ.Safe.,2000,47:298~305.[7]王喜龍,徐福留,李本綱等.天津污灌區(qū)苯并(a)芘、熒蒽和菲生態(tài)毒性的風(fēng)險(xiǎn)表征.城市環(huán)境與城市生態(tài).2002,15:10~12.[8]Wang,X.L.,Tao,S.,Dawson,R.W.,Xu,F.L..CharacterizingandcomparingrisksofpolycyclicaromatichydrocarbonsinaTianjinwastewaterirrigatedarea.Environ.Res.,2002,90:201~206.[9]Logan,D.T.;Wilson,H.T.Anecologicalriskassessmentmethodforspeciesexposedtocontaminantmixtures,Environ.Toxicol.Chem.14,1995,351-359.[10]孟紫強(qiáng),《環(huán)境毒理學(xué)》中國環(huán)境科學(xué)出版社,北京,2000,pp.45.[11] Swartz,R.C.;Schults,D.W.;Ozretich,R.J.;Lamberson,J.O.;Cole,F.A.;DeWitt,T.H;Redmond,M.S;Ferraro,S.P1995ΣPAH:Amodeltopredictthetoxicityofpolynucleararomatichydrocarbonmixturesinfield-collectedsediments,Environ.Toxicol.Chem.14,1995,1977-1987[12] Petry,T.;Schmid,P.;Schlatter,C.Theuseoftoxicequivalencyfactorsinassessingoccupationalandenvironmentalhealthriskassociatedwithexposuretoairbornemixturesofpolycyclicaromatichydrocarbons(PAHs),Chemosphere,32,1996,639-648.[13] Tsai,P.J.;Shieh,H.Y.;Lee,W.J.;Lai,S.O.Health-riskassessmentforworkersexposedtopolycyclicaromatichydrocarbons(PAHs)inacarbonblackmanufacturingindustry,Sci,TotalEnviron,278,2001,137-150,[14]Soldán,P,Toxicriskofsurfacewaterpollution—sixyearsofexperience,Environ,Int.28,2003,677-682. [15] Verschueren,K.HandbookofEnvironmentalDataonOrganicChemicals.VanNostrand-Reinhold,NewYork,2001.[16] 徐鏡波,《環(huán)境毒理學(xué)》東北師范大學(xué)出版社,長春2000,pp.28.[17] 劉筱嫻,《醫(yī)學(xué)統(tǒng)計(jì)學(xué)》科學(xué)出版社,北京2000,pp.157.[18] Mukhtasor;RehanSadiq;TahirHusain;BrainVeitch;NeilBoseAcuteecologicalriskassociatedwithsootdeposition:aPersianGulfcasestudy