天津地区土壤中萘的生态风险分析杨宇,石璇,徐福留,陶澍(北京大学环境学院,地表过程分析与模拟教育部重点实验室,北京100871)摘要:以天津地区表土188个样点的实测萘浓度以及萘对10种代表性土壤生物的半致死浓度为基础资料,尝试用三种方式评价了该地区土壤中萘的生态风险。并进行了方法学比较。研究结果发现,暴露浓度和毒性效应概率密度函数的重叠面积为6.4%,联合概率曲线非常贴近左、下两轴,准确计算的暴露浓度超过LC50的概率仅为1.67%。这些结果均说明该地区土壤萘的生态风险极低,仅严重污染的局部地区可能发生最敏感生物种受害。用于进行风险评价的三种方法各有利弊,从不同角度反映了上述事实。关键词:生态风险;概率风险评价;土壤;萘RiskanalysisofnaphthaleneinTianjinterrestrialecosystemYangYu,ShiXuan,XuFuliu,TaoShu(CollegeofEnvironmentalSciences,,LaboratoryforEarthSurfaceProcesses,PekingUniversity,Beijing100871,China)Abstract:ThreeapproacheswereappliedandevaluatedforprobabilisticriskassessmentofnaphthaleneinaterrestrialecosysteminTianjin,China,basedontheobservednaphthaleneconcentrationof188topsoilsamplesfromtheareaandLC50sofnaphthalenetotentypicalsoilspeciesfromtheliterature.ItwasfoundthattheoverlappingareaofthetwoprobabilitydensityfunctionsofconcentrationandLC50was6.4%,thejointprobabilitycurvebendtowardsandveryclosetothebottomandleftaxis,andthecalculatedprobabilitythatexposureconcentrationexceedsLC50ofvariousspecieswasaslowas1.67%,allindicatingaverymuchacceptableriskofnaphthalenetotheterrestrialecosystemandonlyverysensitivespeciesorindividualanimalsarethreatenbylocalizedextremelyhighconcentration.Thethreeapproachesrevealedsimilarresultsfromdifferentviewpoints.Keywords:ecologicalrisk;probability;riskassessment;soil;naphthalene1.前言作为多环芳烃中迁移活性最强的化合物,萘是天津地区影响范围最大的持久性有机污染物之一。近年监测结果发现该地区表土中萘含量高达275276g/kg[1]。与其他多环芳烃相似,其主要来源包括化石燃料燃烧和污水灌溉[2]。虽然目前并未发现萘具有致癌作用,但其对生态系统的影响早有定论,对多种生物具有毒性[3]。风险分析是定量研究有毒污染物生态危害的重要手段,因而受到各国环境研究者和决策者的高度重视。目前用于风险分析的方法很多既有简单的阈值比较[4],也有具有概率意义的风险分析[5-6]。尽管前者在风险管理中使用方便,但由于忽视了污染物在环境中浓度的大幅度变异以及不同生物种耐受能力的极大差别,简单的阈值很难在生态系统水平对污染的生态风险做出合理的判断,被推荐在生态风险评价的筛选阶段使用[7-8].概率风险分析的方法同时考虑生物耐受性和污染物浓度两方面的变异,可以对污染的生态风险做出整体评价。对生态风险进行有概率意义的评估可以采用不同方法。最常见的是将表征化合物基金项目:北京大学政基金,国家自然科学基金(40031010,40024101)作者简介:杨宇(1984-),男,汉族,本科生暴露浓度和毒性参数的概率密度曲线置于同一坐标系下,据此考察两个密度函数的重叠,并计算其重叠部分面积表述生物受不利影响的概率[9-10]。这样的做法实际上基于浓度分布中心趋势在数值上低于毒性参数分布的隐含假设,当事实上毒性参数分布中心低于浓度分布中心时,浓度分布正移,重叠面积反而减小,简单根据重叠面积的判断已失去原有意义。方法的缺陷显而易见。联合概率曲线则以毒性响应累积概率(横标)和暴露浓度超过相应影响边界浓度的概率(纵标)作图表征特定化合物生态风险[11]。由此得到的曲线在二维平面上的位置反映了风险水平的高低,曲线越靠近坐标轴,风险相应越小。虽然联合概率曲线比简单计算分布重叠面积要合理得多,但仍然缺乏非常确切的概率意义[8]。从理论上说,精确计算整个浓度-响应分布区间风险系数的概率密度函数是比较合理的概率风险分析。这样的计算可以通过蒙特卡罗模拟手段实现12。可以利用这一...