廣州地區(qū)典型多溴聯(lián)苯醚遷移和歸趨行為模擬

田 慧,郭 強,毛瀟萱,黃 韜,馬建民,吳軍年,張 干,李 軍,高 宏*

(1.蘭州大學資源與環(huán)境學院,甘肅 蘭州 730000；2.中國科學院廣州地球化學研究所,有機地球化學國家重點實驗室,廣東 廣州 510640)

廣州地區(qū)典型多溴聯(lián)苯醚遷移和歸趨行為模擬

田 慧1,郭 強1,毛瀟萱1,黃 韜1,馬建民1,吳軍年1,張 干2,李 軍2,高 宏1*

(1.蘭州大學資源與環(huán)境學院,甘肅 蘭州 730000；2.中國科學院廣州地球化學研究所,有機地球化學國家重點實驗室,廣東 廣州 510640)

運用多介質(zhì)逸度模型對典型PBDEs(BDE47、BDE99、BDE209)在廣州地區(qū)大氣、水體、土壤、沉積物中的濃度分布與多介質(zhì)間的遷移、歸趨進行了模擬研究并分析了 3種化合物在研究區(qū)域環(huán)境多介質(zhì)間的遷移通量,確定其在環(huán)境中的主要遷移過程;結(jié)合實際監(jiān)測數(shù)據(jù),對模型的可靠性進行驗證;以BDE47和BDE209為例,對模型的輸入?yún)?shù)進行靈敏度分析;并以BDE209為例,對模型進行不確定性分析.通過模擬濃度與實測濃度的對比,表明模型在該地區(qū)具有很好的適用性.結(jié)果表明,環(huán)境系統(tǒng)達到平衡時,BDE47、BDE99和BDE209在土壤和沉積物中的含量分別占其在環(huán)境系統(tǒng)總含量的17.73%和82.26%,14.65%和85.35%,4.81%和95.19%;PBDEs從環(huán)境系統(tǒng)中的消失途徑主要為大氣平流輸出和土壤降解;logKow和大氣平流輸入是影響化合物在環(huán)境相中濃度分布的最主要因素;不確定性分析指出BDE209在土壤相中濃度的變異系數(shù)最大.

廣州；多溴聯(lián)苯醚(PBDEs)；多介質(zhì)環(huán)境模型；靈敏度分析；不確定性分析

多溴聯(lián)苯醚(PBDEs)是一種添加型溴代阻燃劑,廣泛應用于塑料制品、電子電器、汽車和紡織纖維等產(chǎn)品的生產(chǎn)過程中,由于沒有化學鍵的束縛,因此在產(chǎn)品生產(chǎn)、使用和報廢處置過程中都有可能釋放到環(huán)境中.隨著電子產(chǎn)業(yè)的飛速發(fā)展,全球 PBDEs的需求量和消耗量不斷增加[1-2].同時,由于其持久性、半揮發(fā)性和長距離遷移性而廣泛的存在于各環(huán)境介質(zhì)中[3-9],并且在偏遠地區(qū)也發(fā)現(xiàn)了PBDEs的存在[10-11].

珠江三角洲是我國最大的計算機及其零部件、通訊器材生產(chǎn)地以及電子垃圾拆卸地.環(huán)境介質(zhì)中的PBDEs的濃度在過去以低溴BDEs為主,現(xiàn)在 BDE209的濃度已經(jīng)遠遠超出了其余PBDEs同系物之和[12].目前,PBDEs已廣泛的存在于珠江三角洲各環(huán)境介質(zhì)中[13-18].廣州市是珠江三角洲地區(qū)的工業(yè)及經(jīng)濟中心,夏季大氣中以四、五和十溴BDEs為主,其含量占∑PBDEs的80%以上,其中 BDE209、BDE47、BDE99是主要成分,尤以BDE209高于香港[19]、Michigan[20]、Indiana[20]等沒有電子垃圾回收廠的地區(qū).研究表明,廣州地區(qū)PBDEs的污染主要來自五溴和十溴工業(yè)品的使用[21].

國外關(guān)于 PBDEs在環(huán)境中行為的研究大都集中在水生食物網(wǎng)模型[22-25];也有一些關(guān)于長距離遷移[26-27]及室內(nèi)污染物來源分布[28-29]的研究探討;此外,Palm[30]使用多介質(zhì)逸度模型對瑞典環(huán)境相中的PBDEs濃度進行了模擬.國內(nèi)研究則比較單一,鮮見關(guān)于PBDEs在環(huán)境中行為趨勢的報道.因此,本文針對廣泛存在于廣州各環(huán)境介質(zhì)中的BDE47、BDE99和BDE209進行多介質(zhì)間的遷移和歸趨模擬,更全面的了解典型PBDEs在廣州地區(qū)的行為,為評價其對廣州地區(qū)人群的潛在生態(tài)風險以及污染物的控制提供基礎理論支持.

1 研究方法

1.1 模型概況

本文基于逸度[31](Fugacity)概念,選取大氣、水體、土壤和沉積物四個環(huán)境主相(其中大氣相包括氣子相和氣溶膠子相,水體相包括水子相和懸浮顆粒物子相,土壤相包括氣子相、水子相和固體子相,沉積物相包括固體沉積物子相和水子相),在穩(wěn)態(tài)、平衡的假設條件下,結(jié)合廣州地區(qū)環(huán)境特點,建立 PBDEs的多介質(zhì)遷移歸趨模型.所涉及的過程包括化合物在相鄰環(huán)境相界面間的遷移、大氣水體平流、各環(huán)境相中的降解、地表水徑流、表層沉積物的掩埋等.模型的輸入?yún)?shù)包括化合物的物理化學性質(zhì)和研究區(qū)的環(huán)境參數(shù),輸出結(jié)果則包括各環(huán)境相及其子相的逸度f(Pa)、逸度容量 Z(mol/(m3?Pa))、化合物在各環(huán)境介質(zhì)間的遷移通量 D(mol/(h?Pa))、遷移速率N(mol/h)、以及當環(huán)境系統(tǒng)達到平衡時化合物在各環(huán)境相中的濃度C(mol/m3)等.

1.2 數(shù)據(jù)來源

模型需要輸入的參數(shù)主要包括 BDE47、BDE99、BDE209的物理化學性質(zhì)參數(shù)(蒸汽壓、辛醇-水分配系數(shù)、溶解度、在各環(huán)境介質(zhì)中的半衰期等)和研究區(qū)的環(huán)境參數(shù)(溫度、降雨量、面積、有機碳含量等),具體參數(shù)見表1和表2.

表1 PBDEs主要物理化學參數(shù)Table 1 Main physical and chemical parameters of PBDEs

Palm等[13]估算出PBDEs的人均年排放量(表3),該排放量適用于沒有生產(chǎn)PBDEs工業(yè)品工廠的地區(qū),廣州市人口2010年為1270.08萬,則據(jù)此可以計算出各污染物向環(huán)境中的直接排放量(表3).

由于廣州冬季上風向地區(qū)空氣中PBDEs濃度數(shù)據(jù)的缺失,采用廣州夏季上風向地區(qū)空氣中PBDEs的濃度作為模型中大氣平流輸入污染物的濃度.因廣州夏季盛行東或東南季風,因此選取了從東莞方向平流輸入的 BDE47、BDE99、BDE209,濃 度 依 次 為 2181.47, 1058.21, 4191.93pg/m3[21].水體平流輸入濃度參考珠江入?？谒w中 PBDEs濃度,依次為 0.04,0.001,6.37ng/L[32].

表2 研究區(qū)主要環(huán)境參數(shù)Table 2 Main environmental parameters of the study area

表3 廣州地區(qū)PBDEs排放量Table 3 Emission of PBDEs in Guangzhou

1.3 靈敏度分析和不確定性分析

為確定輸入?yún)?shù)對模型輸出結(jié)果的影響程度, 以BDE47和BDE209為例,選取除常數(shù)外的代表性輸入?yún)?shù),對參數(shù)分別進行±10%的變動,依據(jù) CS=(Y1.1-Y0.9)/(0.2×Y1.0)進行靈敏度分析,其中Y1.1、Y0.9、Y1.0分別為參數(shù)取均值1.1、0.9和1.0倍時模型的響應值,Cs為靈敏度系數(shù).本文以|CS|>0.5為篩選標準,對模型的關(guān)鍵參數(shù)進行識別.

同時,以 BDE209在各環(huán)境介質(zhì)中濃度為例,運用蒙特卡羅方法對模型的輸出結(jié)果進行不確定性分析,對已經(jīng)篩選出的靈敏度高的參數(shù)進行隨機取值,程序重復運行10000次,給出模型的不確定性信息.本文所有程序運算均在Matlab7.0中進行.

2 結(jié)果與討論

2.1 模型驗證和濃度分布

為驗證模型的適用性,從文獻中搜集了 3種化合物在廣州地區(qū)大氣[21,19,39]、土壤[17]、沉積物[12,18]中的實測濃度,并且參考珠江入?？谒w中PBDEs濃度,將模擬濃度與實測濃度的對數(shù)值進行比較.由圖1可見∶3種化合物在各環(huán)境介質(zhì)中的絕大多數(shù)模擬濃度與實測濃度的吻合度較高, BDE99在土壤中的濃度略高于實測值,BDE209在土壤和沉積物中的濃度均接近于實測濃度的最小值,原因可能是過小的估計了BDE209的排放量,廣州是電子產(chǎn)業(yè)集中地區(qū),雖然《斯德哥爾摩公約》第4次締約方大會已經(jīng)將五溴和八溴聯(lián)苯醚列入 POPs名單,但全球十溴聯(lián)苯醚的需求量仍然很大,由此可見導致偏差的原因極有可能是Palm等[13]對BDE209排放的估算低于高速發(fā)展下的廣州的實際排放量.總體上來說,模型能較準確地反映研究區(qū)的污染狀況,因此該模型在該研究區(qū)具有一定的適用性.

由表 4可知,3種化合物在土壤和沉積物中的濃度遠遠超過了其在大氣和水體中的濃度,當研究區(qū)環(huán)境達到平衡時,BDE47在土壤和沉積物相中所占質(zhì)量比率分別為 17.73%和 82.26%, BDE99在土壤和沉積物相中所占比率分別為14.65%和85.35%,BDE209在土壤和沉積物相中所占比率分別為4.81%和95.19%,由此可見,3種同系物在土壤和沉積物中的質(zhì)量占環(huán)境中總質(zhì)量的 98%以上,從而表明土壤和沉積物是目標化合物在研究區(qū)的主要儲庫.同時,隨著溴代程度的升高,PBDEs在土壤中的富集程度會相對降低,在表層沉積物中的富集程度會增加,在實際環(huán)境中可隨著表層沉積物的沉降進入到深層的沉積物.BDE47、BDE99、BDE209在大氣氣溶膠子相和水體懸浮顆粒物子相中質(zhì)量占對應總相質(zhì)量的比例分別為39.22%、26.59%,70.13%、73.76%, 96.31%、100%,表明隨著溴代數(shù)的增加,化合物在各子相之中更傾向于向固體附著.

2.2 多介質(zhì)相間遷移與歸趨行為

由圖2可知,3種化合物的消失途徑基本相同,主要為大氣平流輸出和土壤降解,并且前者速率高出后者約1個數(shù)量級,BDE47、BDE99和BDE209的特征遷移距離分別為 811,600,209km[38],表明PBDEs易通過大氣向下風向方向進行傳輸,說明廣州可成為其下風向地區(qū)的污染源;界面間則以大氣-土壤、水體-沉積物遷移為主,表明PBDEs在進入環(huán)境后通過一系列的遷移過程,最終賦存在土壤和沉積物當中,究其原因,可能是其低水溶性和高辛醇-水分配系數(shù)(logKow)所決定的[17].通過對PBDEs進入到環(huán)境中各種行為速率的分析,基本可以推測出其主要行為∶大氣中的化合物大部分遷移出了研究的環(huán)境相,對下風向的地區(qū)造成一定程度的污染;土壤降解是污染物消失在研究區(qū)環(huán)境中的第二個主要途徑;大氣中未遷出的污染物則主要隨著干濕沉降、擴散等過程進入到土壤、水體等環(huán)境中,最終賦存在土壤和沉積物當中.

圖1 PBDEs模擬濃度與實測濃度的對比Fig.1 Comparsion between the calculated and the observed concentration of PBDEs

表4 PBDEs在環(huán)境介質(zhì)中模擬的濃度(mol/m3)Table 4 Calculated concentrations of PBDEs in environmental media(mol/m3)

對不同溴代程度的同系物進行分析∶BDE47在大氣-水體和水體-大氣之間的交換速率相差4個數(shù)量級,而 BDE99則相差不到 5個數(shù)量級,BDE209相差5個多數(shù)量級,并且3種同系物從大氣-水體凈遷移速率隨著溴代程度的增加而增加,溴代程度高的 BDE209更容易從大氣沉降.Palm[30]指出由于BDE209低的蒸汽壓,因此更易沉降到地表,但空氣平流輸出依然是其遷出研究區(qū)大氣環(huán)境相的主要途徑;吳有方等[38]的研究表明PBDEs同系物隨著溴代程度的增加,更容易形成點源污染,低溴聯(lián)苯醚進行遠距離遷移的潛力最大.而在水環(huán)境中,PBDEs則更傾向于向顆粒態(tài)分配,隨著懸浮顆粒物的沉降和再懸浮在水體和沉積物之間遷移,BDE47、BDE99、BDE209在水體和沉積物之間的遷移速率的對數(shù)值分別為 0.353,0.354,0.354log(mol/h),各同系物之間并沒有太大的差別.BDE47和 BDE99在環(huán)境中主要消失途徑相同,除了隨大氣平流遷出研究區(qū)和土壤中的降解外,還有大氣中的降解,而對于BDE209,除了大氣平流和土壤降解外,最主要的過程則為沉積物中的降解.這與Palm[30]的研究稍許有些差異,他指出PBDEs在環(huán)境中消失的途徑主要為大氣平流輸出、土壤向水體的地表徑流、土壤中的降解,但兩篇文章均認為其最主要的消失途徑為大氣的平流輸出,并且環(huán)境中 80%以上的PBDEs均通過此途徑在研究區(qū)消失,但會對周邊地區(qū)產(chǎn)生污染.

圖2 PBDEs在多環(huán)境介質(zhì)中的遷移速率Fig.2 Transfer fluxes of PBDEs in the multimedia environment

2.3 靈敏度分析

表5 靈敏度分析中的參數(shù)代碼Table 5 Codes of parameters for sensitivity analysis

以BDE47和BDE209在環(huán)境介質(zhì)中的濃度為例,對模型的輸入?yún)?shù)進行靈敏度分析,其中分析代碼見表5所示,部分分析結(jié)果見圖3,靈敏度系數(shù)大于 0表明參數(shù)對模型輸出結(jié)果具有正影響作用,反之,小于0則為負影響.

兩種化合物的大氣平流輸入對所有環(huán)境介質(zhì)中的濃度分布均有很大的影響;并且PBDEs在大氣相中的濃度分布在|Cs|>0.5的篩選標準下僅僅受平流輸入的影響,因此圖 3并未列出BDE47和BDE209在大氣相中濃度分布的靈敏度分析結(jié)果;此外,logKow對 BDE47在水體中濃度影響最大,靈敏度系數(shù)可達到-7.4,對BDE209在水體中的濃度影響最大的為懸浮顆粒物的沉降速率;對兩種化合物在土壤相中濃度影響較大的參數(shù)則主要為土壤的性質(zhì)(土壤密度、土壤深度、土壤中固體體積分數(shù)、土壤中的半衰期)和降雨;沉積物中BDE47濃度分布的靈敏參數(shù)主要為 logKow,其靈敏度系數(shù)可達 3.45,對 BDE209靈敏度影響較大的參數(shù)則主要為沉積物的性質(zhì)(水體面積、表層沉積物的深度、沉積物中固體體積分數(shù)、懸浮顆粒物的再懸浮速率、沉積物中的半衰期)以及水體平流輸入濃度等.

對所有靈敏度系數(shù)進行比較可知∶logKow、大氣平流輸入濃度是影響化合物在環(huán)境相中濃度分布的最主要因素.

圖3 BDE47和BDE209在環(huán)境介質(zhì)中濃度的靈敏度分析Fig.3 Sensitivity parameters of BDE47 and BDE209 concentrations in environmental compartments

2.4 不確定性分析

不確定性分析結(jié)果如圖 4所示.BDE209在大氣

中濃度范圍為2603～4048pg/m3,水體中濃度范圍為 0.61～1.46ng/L,土壤中濃度范圍為 1.62～4.56ng/g,沉積物中濃度范圍為48.34～116.78ng/g. BDE209在大氣、水體、土壤和沉積物中模擬濃度分布的變異系數(shù)分別為5.25、10.44、12.07、11.05,其中,土壤的變異系數(shù)最大,這可能與進行不確定性分析的參數(shù)的靈敏度和不同參數(shù)間靈敏度的相對平衡有關(guān)系.

圖4 BDE209模擬計算濃度頻數(shù)分布Fig.4 Frequency distribution of BDE209 calculated concentration

3 結(jié)論

3.1 模型模擬的3種PBDEs同系物的濃度與研究區(qū)實測濃度的吻合度較高,表明該模型適用于研究區(qū)PBDEs的模擬.通過質(zhì)量分布比例確定廣州地區(qū)所研究的3種BDEs同系物的最終儲庫為土壤和沉積物,占各自總環(huán)境質(zhì)量的98%以上.

3.2 PBDEs在環(huán)境相中的消失途徑主要為大氣平流輸出和在土壤中的降解,可能成為下風向城市的污染源;界面間則以大氣-土壤、水體-沉積物的遷移為主.

3.3 logKow和大氣的平流輸入對所有環(huán)境相中PBDEs濃度均有較大的影響,并且針對各個具體的環(huán)境介質(zhì),參數(shù)的靈敏度系數(shù)大小均不相同,其中對土壤相和沉積物相靈敏度影響因素較多.

3.4 BDE209在大氣、水體、土壤和沉積物中不確定性分析得出各介質(zhì)中的模擬濃度范圍分別為 2603～4048pg/m3、0.61～1.46ng/L、 1.62～4.56ng/g和48.34～116.78ng/g.土壤中濃度的相對離散程度大于其余3介質(zhì)中濃度的離散程度.

[1] 葉細標,傅 華.多溴聯(lián)苯醚的環(huán)境暴露及健康危害 [J]. 環(huán)境與職業(yè)醫(yī)學, 2007,24(1)：95-101.

[2] de Wit CA. An overview of brominated flame retardants in the environment [J]. Chemosphere, 2002,46(5)：583-624.

[3] 胡永彪,李英明,耿大瑋,等.北京冬季大氣中多溴聯(lián)苯醚的污染水平和分布特征 [J]. 中國環(huán)境科學, 2013,33(1)：9-13.

[4] 楊 雪,劉大錳,孫俊玲.北京市大氣 PM2.5中多溴聯(lián)苯醚污染水平與來源分析 [J]. 現(xiàn)代地質(zhì), 2010,24(2)：355-361.

[5] 謝 薇,孫天河,韓 鵬,等.青藏高原中部地區(qū)表層土壤中多溴聯(lián)苯醚的分布特征及影響因素 [J]. 現(xiàn)代地質(zhì), 2012,26(5)：917-925.

[6] 楊永亮,潘 靜,李 悅,等.青島近岸沉積物中持久性有機污染物多氯萘和多溴聯(lián)苯醚 [J]. 科學通報, 2003,48(21)：2244-2251.

[7] 周變紅,張承中,蔣君麗,等.西安市大氣中多溴聯(lián)苯醚的季節(jié)變化特征及人體暴露評估 [J]. 中國環(huán)境科學, 2012,32(9)：1591-1596.

[8] 陳社軍,田 密,羅 勇,等.多溴聯(lián)苯醚在電子廢棄物拆卸地沉積物中的垂直分布 [J]. 環(huán)境科學, 2010,31(12)：3088-3092.

[9] 劉芃巖,高 麗,趙雅嫻,等.分散液相微萃取氣相色譜/氣相色譜質(zhì)譜法測定白洋淀水中多溴聯(lián)苯醚 [J]. 分析化學, 2010,38(4)：498-502.

[10] Hites RA. Polybrominated diphenyl ethers in the environment and in people： a meta-analysis of concentrations [J]. Environmental Science and Technology, 2004,38(4)：945-956.

[11] Wang X M, Ding X, Mai B X, et al. Polybrominated diphenyl ethers in airborne particulates collected during a research expedition from the Bohai Sea to the Arctic [J]. Environmentalscience and technology, 2005,39(20)：7803-7809.

[12] 陳社軍.珠江三角洲河流、河口及鄰近南海海域和長江三角洲主要水體沉積物中的多溴聯(lián)苯醚 [D]. 廣州：中國科學院廣州地球化學研究所,2006.

[13] Palm A, Cousins I T, Mackay D, et al. Assessing the environmental fate of chemicals of emerging concern： a case study of the polybrominated diphenyl ethers [J]. Environmental Pollution, 2002,117(2)：195-213.

[14] 黃玉妹,陳來國,文麗君,等.廣州市室內(nèi)塵土中多溴聯(lián)苯醚的分布特點及來源 [J]. 中國環(huán)境科學, 2009,29(11)：1147-1152.

[15] 陳多宏,畢新慧,鐘流舉,等.典型電子垃圾拆解區(qū)大氣中多溴聯(lián)苯醚(PBDEs)的季節(jié)變化特征 [J]. 生態(tài)環(huán)境學報, 2010,19(5)：1068-1072.

[16] 羅孝俊,余 梅,麥碧嫻.多溴聯(lián)苯醚(PBDEs)在珠江水體中的分布于分配 [J]. 科學通報, 2008,53(2)：141-146.

[17] 鄒夢瑤,龔 劍,冉 勇.珠江三角洲流域土壤多溴聯(lián)苯醚(PBDEs)的分布及環(huán)境行為 [J]. 生態(tài)環(huán)境學報, 2009,18(1)：122-127.

[18] 陳社軍,麥碧嫻,曾永平.珠江三角洲及南海北部海域表層沉積物中多溴聯(lián)苯醚的分布特征 [J]. 環(huán)境科學學報, 2005,25(9)：1265-1271.

[19] Li J, Liu X, Yu L-L, et al. Comparing polybrominated diphenyl ethers (PBDEs) in airborne particles in Guangzhou and Hong Kong： sources, seasonal variations and inland outflow [J]. Journal of Environmental Monitoring, 2009,11(6)：1185-1191.

[20] Hoh E, Hites R A. Brominated flame retardants in the atmosphere of the East-Central United States [J]. Environmental Science and Technology, 2005,39(20)：7794—7802.

[21] 陳來國.廣州市夏季大氣中多氯聯(lián)苯和多溴聯(lián)苯醚的含量及組成對比 [J]. 環(huán)境科學學報, 2008,28(1)：150-159.

[22] Bhavsar S P, Gandhi N, Gewurtz S B, et al. Fate of PBDEs in juvenile lake trout estimated using a dynamic multichemical fish model [J]. Environmental Science and Technology, 2008,42(10)：3724-3731.

[23] Gandhi N, Bhavsar S P, Gewurtz S B, et al. Development of a multichemical food web model： Application to PBDEs in Lake Ellasj?en, Bear Island, Norway [J]. Environmental Science and Technology, 2006,40(15)：4714-4721.

[24] Mattila T J, Verta M. Modeling the importance of biota and black carbon as vectors of polybrominated diphenyl ethers (PBDEs) in the Baltic Sea ecosystem [J]. Environmental Science and Technology, 2008,42(13)：4831-4836.

[25] Wu J-P, Luo X-J, Zhang Y, et al. Biomagnification of polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls in a highly contaminated freshwater food web from South China [J]. Environmental Pollution, 2009,157(3)：904-909.

[26] Gouin T, Harner T. Modelling the environmental fate of the polybrominated diphenyl ethers [J]. Environment International, 2003,29(6)：717-724.

[27] Wania F, Dugani C B. Assessing the long-range transport potential of polybrominated diphenyl ethers： A comparison of four multimedia models [J]. Environmental Toxicology and Chemistry, 2003,22(6)：1252-1261.

[28] Jones-Otazo H A, Clarke J P, Diamond M L, et al. Is house dust the missing exposure pathway for PBDEs? An analysis of the urban fate and human exposure to PBDEs [J]. Environmental Science and Technology, 2005,39(14)：5121-5130.

[29] Zhang X, Diamond M L, Ibarra C, et al. Multimedia modeling of polybrominated diphenyl ether emissions and fate indoors [J]. Environmental Science and Technology, 2009,43(8)：2845-2850.

[30] Palm A.The Environmental fate of polybrominated diphenyl ethers in the centre of Stockholm： Assessment using a multimedia fugacity model [R]. Swedish, 2001

[31] 唐納德?麥凱.環(huán)境多介質(zhì)模型-逸度方法 [M]. 北京：化學工業(yè)出版社, 2007.

[32] 管玉峰,涂秀云,吳宏海.珠江入?？谒w中多溴聯(lián)苯醚及其來源分析 [J]. 生態(tài)環(huán)境學報, 2011,20(3)：474-479.

[33] ?rn U, Eriksson L, Jakobsson E, et al. Synthesis and characterization of polybrominated diphenyl ethers-unlabelled and radiolabelled tetra-, penta-and hexa-bromodiphenyl ethers [J]. Acta Chemica Scandinavica, 1996,50(9)：802-807.

[34] Hansen B, Munn S, de Bruijn J, et al. EUR 20402EN- European Union Risk Assessment Report bis (pentabromophenyl) Ether [R]. European Commission. 2002：1-282.

[35] 中國氣象局.中國氣象科學數(shù)據(jù)共享服務網(wǎng)：中國(1971-2000年)氣候標準值[EB/OL].北京：國家氣象信息中, 2009[2009-04-05].http：//cdc.cma.gov.cn/shishi/climate.jsp?stprovid= 廣 東&station=59287.

[36] 廣州市統(tǒng)計局,廣州統(tǒng)計年鑒2010 [M]. 北京：中國統(tǒng)計出版社, 2010.

[37] 廣州市志(卷二) [M]. 廣州市地方編纂委員會.廣州：廣州出版社, 1998.

[38] 吳有方,方利江,丁中原,等.廣州地區(qū) PBDEs長距離遷移潛力和總持久性模擬 [J]. 中國環(huán)境科學, 2012,32(3)：426-432.

[39] Chen L-G, Mai B-X, Bi X-H, et al. Concentration levels, compositional profiles, and gas-particle partitioning of polybrominated diphenyl ethers in the atmosphere of an urban city in South China [J]. Environmental Science and Technology, 2006,40(4)：1190-1196.

Simulating the transfer and fate of typical PBDEs in Guangzhou.

TIAN Hui1, GUO Qiang1, MAO Xiao-xuan1, HUANG Tao1, MA Jian-min1, WU Jun-nian1, ZHANG Gan2, LI Jun2, GAO Hong1*
(1.College of Resource and Environment, Lanzhou University, Lanzhou 730000, China；2.State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China). China Environmental Science, 2014,34(3)：758～765

A multimedia fugacity model was applied to simulate the concentrations distribution, transfer and fate of typical polybrominated diphenyl ethers (BDE47、BDE99、BDE209) in air, water, soil, sediment of Guangzhou Area. Meanwhile, the transfer fluxes between different compartments were analyzed in order to infer the main transfer process. Moreover, the simulated concentrations in air, soil and sediment were compared to monitored data for validation purpose. In addition taking BDE47and BDE209as samples, input parameters of the model were tested and the key parameters were identified using sensitivity analysis method, while the uncertainty of these key parameters was estimated taking BDE209as sample. The reliability of the model was verified by the agreement between calculated and measured concentrations. When the system reaches equilibrium, soil and sediment compartments were the main reservoirs of PBDEs in Guangzhou area because of the mass fraction of BDE47、BDE99 and BDE209 in soil and sediment of the total content, respectively, were 17.73% and 82.26%, 14.65% and 85.35%, 4.81% and 95.19%;The air advection outflow and soil degradation were the major routes for PBDEs to disappear in the study area. The results of sensitivity analysis in this study also indicated that the logKow、air advection inflow had significant influence on concentrations of PBDEs in various media. Uncertainty analysis for environmental multimedia fugacity model predicted that the coefficient of variation of BDE209in the soil concentrations had the most significant compared with others.

Guangzhou；polybrominated diphenyl ethers；multimedia environmental model；sensitivity analysis；uncertainty analysis

X131

：A

：1000-6923(2014)03-0758-08

田 慧(1989-),女,山西臨汾人,蘭州大學碩士研究生,研究方向為干旱區(qū)持久性有機物區(qū)域環(huán)境行為.

2013-06-20

國家自然科學基金項目(40971267/D010903);教育部“春暉計劃”國際合作科研項目(Z2008-1-62025)

* 責任作者, 教授, honggao@lzu.edu.cn