天津地区土壤中萘的生态风险分析

SeriesofSelectedPapersfromChun-TsungScholars,PekingUniversity(2003)146天津地区土壤中萘的生态风险分析杨宇,石璇,徐福留,陶澍（北京大学环境学院，地表过程分析与模拟教育部重点实验室，北京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.前言作为多环芳烃中迁移活性最强的化合物，萘是天津地区影响范围最大的持久性有机污染物之一。近年监测结果发现该地区表土中萘含量高达275276g/kg[1]。与其他多环芳烃相似，其主要来源包括化石燃料燃烧和污水灌溉[2]。虽然目前并未发现萘具有致癌作用，但其对生态系统的影响早有定论,对多种生物具有毒性[3]。风险分析是定量研究有毒污染物生态危害的重要手段，因而受到各国环境研究者和决策者的高度重视。目前用于风险分析的方法很多，既有简单的阈值比较[4]，也有具有概率意义的风险分析[5-6]。尽管前者在风险管理中使用方便，但由于忽视了污染物在环境中浓度的大幅度变异以及不同生物种耐受能力的极大差别，简单的阈值很难在生态系统水平对污染的生态风险做出合理的判断，被推荐在生态风险评价的筛选阶段使用[7-8].概率风险分析的方法同时考虑生物耐受性和污染物浓度两方面的变异，可以对污染的生态风险做出整体评价。对生态风险进行有概率意义的评估可以采用不同方法。最常见的是将表征化合物基金项目：北京大学政基金，国家自然科学基金（40031010，40024101）作者简介：杨宇（1984-），男，汉族，本科生SeriesofSelectedPapersfromChun-TsungScholars,PekingUniversity(2003)147暴露浓度和毒性参数的概率密度曲线置于同一坐标系下，据此考察两个密度函数的重叠，并计算其重叠部分面积表述生物受不利影响的概率[9-10]。这样的做法实际上基于浓度分布中心趋势在数值上低于毒性参数分布的隐含假设，当事实上毒性参数分布中心低于浓度分布中心时，浓度分布正移，重叠面积反而减小，简单根据重叠面积的判断已失去原有意义。方法的缺陷显而易见。联合概率曲线则以毒性响应累积概率（横标）和暴露浓度超过相应影响边界浓度的概率（纵标）作图表征特定化合物生态风险[11]。由此得到的曲线在二维平面上的位置反映了风险水平的高低，曲线越靠近坐标轴，风险相应越小。虽然联合概率曲线比简单计算分布重叠面积要合理得多，但仍然缺乏非常确切的概率意义[8]。从理论上说，精确计算整个浓度－响应分布区间风险系数的概率密度函数是比较合理的概率风险分析。这样的计算可以通过蒙特卡罗模拟手段实现12。可以利用这一概念计算暴露浓度超过特定耐受水平的概率，并据此进行评价[12]。本研究旨在根据天津地区萘的土壤含量实测数据和文献中报道的萘对土壤生物的毒性资料，分别用上述三种方法表征萘对该地区土壤生态系统的风险，着重对三种方法进行方法学比较，为更深入地研究本地区多环芳烃和其他持久性有机污染物毒性建立必要手段。2.研究方法土壤萘含量资料来自对取自天津地区的188个表土样品的分析。采样点大致均匀分布于天津地区，平均采样密度约为1/60km2。所有样品均作平行测定。样品中多环芳烃含量均用加速溶剂仪（DionexASE-300）提取，以1:1二氯甲烷－丙酮为提取液。提取温度140C，压强1500psi，预热５min，静态提取５min。提取液用硅胶层析柱净化。样品多环芳烃（PAH16）用GC-MS（安捷伦6890GC/5973MSD）以外标定量法测定。详细分析方法及仪器条件设置参见相关文献[13]。从美国环保署毒性数据库中收集萘对土壤生物的毒性资料，主要为土壤动物7-8d暴露的半致死浓度。涉及生物、实验条件和LC50值列于表1。尽管能收集到的毒性参数不可能包括全部土壤生物，而且采用生物并非来自研究地区，但这些资料涉及到多数代表性物种，且一般而言，不同区域同一物种的耐毒差异远低于不同物种间差异，这些资料可以用于本项研究。表1土壤动物毒性参数[拉丁名称汉译[14]]土壤生物拉丁实验周期(d)LC50(ug/g)螨属亚纲acari816节腹亚目arthropleona814节足动物arthropoda811弹尾目collembola815爱胜蚓Eiseniafeida24670跳虫spider812莴苣lactuca7100中气门亚目mesostigmata813甲螨亚目Oribatida815前气门亚目prostigmata817根据上述资料分别求得区域暴露浓度和生物耐受水平的两个概率密度函数，用Matlab计算两函数曲线重叠部分的面积和暴露浓度超过生物耐受水平的概率，并得到联合概率曲线。3．结果与讨论3.1天津地区土壤萘含量及土壤动物对萘毒性耐受水平分布特征图1分别为天津表层土壤中萘含量（左）和土壤动物对萘毒性耐受水平的统计分布（右）。两者分别表现了天津地区土壤中萘的含量和对土壤生物毒性的分布特征。天津地区表土中萘的平均含量为0.275g/g，标准差为0.176g/g。经分布检验（Kolmogorov-Smirnov）后，土壤萘含量服从典SeriesofSelectedPapersfromChun-TsungScholars,PekingUniversity(2003)148型的对数正态分布。土壤生物对萘的耐受水平的数据量较小，即便作正态检验，结论的可靠性也不高，故直接作对数正态假定[11]。根据实际观测结果得到的两个变量的概率密度函数为：311352.071999.0lg1210ln011113.1)(CnapeCnapCnapf329996.1483671.1lg5050225010ln489215.0)(LCeLCLCf式中Cnap为土壤萘浓度(ug/g)，LC50则代表表述萘对土壤动物毒性的半致死浓度。3.2概率密度函数重叠面积评价将浓度和毒性耐受水平两个概率密度函数放在同一坐标系中，并考察两者重叠程度，结果见图2。图中阴影部分即为两函数重叠部分。对天津土壤中萘而言，实际暴露浓度水平总体上大大低于导致10种代表性土壤生物对萘的耐受性水平，即在同一坐标上，前者分布显著较后者偏左。因此可以根据两者重叠面积评价危害风险。从图2可以看出，两个概率密度函数重叠部分极小。计算概率仅为6.4％（图2中阴影部分）。说明天津地区萘对土壤生物的危害极低。3.3联合概率曲线评价可以综合反映致毒效应的LC50的累积函数和土壤萘暴露浓度的反累积函数（与累积函数和为1）进一步计算图3所示联合概率曲线，籍以表述土壤萘对土壤生物毒性的全面关系。图中横标为累积危害概率，纵标则代表暴露浓度高于影响边界浓度（相应的LC50）的概率。王喜龙等人（2002）曾经用类似方式计算天津地区苯并芘、荧蒽和菲三者对水生生物的毒1.51.00.503.0表土萘含量(g/g)2.00f1(Cnap)的概率1.010500.045毒性耐受水平(g/g)0.0300f2(LC50)的概率0.015图1天津地区表土萘的含量分布（左）和土壤动物对萘毒性耐受水平（右）的分布特征1.51.00.502.0浓度(g/g)0f1(Cnap)或f2(LC50)的概率1.0毒性耐受水平表土萘含量图2暴露浓度与生物耐受水平的重叠，重叠部分用阴影表示0.060.040.0201.0累积危害概率0.50.5超过影响边界的概率的概率图3天津地区土壤萘对生物毒性的联合概率曲线SeriesofSelectedPapersfromChun-TsungScholars,PekingUniversity(2003)149性，结果发现三者对天津地区水生生物的毒性远远高于本研究结果。三种PAHs浓度分布和毒性分布曲线重叠概率均在50％以上，分别达0.582、0.623和0.752。相应的联合概率曲线在位置上也很接近斜率为－1的直线。反映出高的多的毒性风险。相比之下，天津地区土壤中萘对土壤生物的毒性相对低得多，概率密度函数重叠仅5％左右，联合概率曲线也非常贴近左、下两轴（横标累积毒性概率范围远远低于1.0）。两者都反映了很低的暴露风险。3.4暴露浓度超过LC50概率的直接计算如果将土壤环境中萘的浓度和萘对不同土壤生物物种毒性的LC50视为两个独立变量，据此可以得到关于两个变量的二元概率密度函数。该函数服从二元对数正态分布。通过二重积分计算暴露浓度超过生物LC50的概率P计算公式为：00502150)(CnapCnapLCddLCfCnapfP据此求得天津地区土壤中萘浓度超过土壤生物LC50值得概率为1.67%。同样反映了相对较低的风险水平。3.5方法学比较与结果尽管三种方法可以得出非常相似的风险分析结果，但三者既有内在联系也有差别。概率密度函数是累积函数的导数，反映了累积概率增加或反累积概率下降的速率。随浓度升高，生物耐受水平的累积概率增加，而暴露浓度的反累积概率则减小。图2所示重叠的概率密度函数交于一点，该点左侧暴露浓度概率密度较大，即暴露浓度反