同一物种对同一污染物的敏感度可以相差100倍——拟长腹剑水蚤的镉半数致死浓度仅19μg/L,而天津厚蟹高达50,000μg/L。这种跨物种的毒性差异,正是物种敏感度分布法存在的前提。生态环境部2026版镉基准报告基于78个海洋物种的198条急性毒性数据,用正态分布模型拟合出了一条从最敏感到最耐受物种的连续分布曲线,从曲线尾部外推的HC₅=70μg/L成为短期基准定值的核心依据。
SSD法的核心逻辑:从单物种到全群落
SSD法假设不同物种对同一污染物的敏感性服从某种统计分布(常用正态分布或逻辑斯谛分布)。通过拟合分布曲线,可以外推保护X%物种的浓度阈值(HCₓ)。该方法的核心优势在于将保护对象从"单一敏感物种"扩展到"群落水平的物种多样性",符合生态保护的实质目标。HJ1260—2022规定中国海洋生物水质基准推导必须使用SSD法,且受试物种需涵盖3个营养级至少5科8种——中国2026版镉基准短期78个物种、51个科,远超最低要求。
SSD法对数据质量的要求极为严格:每条毒性数据必须是物种水平的效应浓度(EC₅₀/LC₅₀等),效应终点必须具有生态相关性(生长、繁殖、存活),组织、细胞或分子水平的亚个体效应数据一律剔除。这也是2,426条原始数据最终仅254条进入SSD拟合的根本原因——数据质量比数据数量更关键。
同效应急/慢性值计算的数学规范
短期数据(ATV)按效应类别分为生长类(EC₅₀)和存活类(LC₅₀)。分物种计算各效应类别的几何均值得到AVE,取两类中较小值作为最敏感AVE纳入模型。公式为:AVEᵢₖ = ⁿ√(ATVᵢₖ₁ × ATVᵢₖ₂ × … × ATVᵢₖₙ)。78个物种中仅文蛤和真赤鲷同时获得两类效应AVE——文蛤生长类84μg/L、存活类127μg/L,取84;真赤鲷生长类800μg/L、存活类8,042μg/L,取800。其余物种单类数据直接纳入。
长期数据(CTV)按生长、繁殖、存活三类别处理。对于NOEC和LOEC配对数据,先计算MATC=√(NOEC×LOEC),再与其他慢性指标(EC₁₀、EC₂₀、EC₅₀)按优先序(EC₁₀>EC₂₀>MATC>NOEC>LOEC>EC₅₀>LC₅₀)选取并计算几何均值得到CVE。26个物种中仅长牡蛎和蒙古裸腹溞获得多类CVE,取最小CVE纳入模型。长期最敏感数据为紫贻贝CVE=6.25μg/L(繁殖类),远低于同物种的存活类数据。
正态分布vs逻辑斯谛分布:模型选择的统计验证
将最敏感AVE和CVE取常用对数后从小到大排序,计算累积频率F_R=Σfᵢ/(N+1)×100%,以lg毒性值为自变量x、累积频率F_R为因变量y,分别拟合正态分布和逻辑斯谛分布模型。拟合软件为"国家生态环境基准计算软件物种敏感度分布法(1.0版)"。
评价参数包括RMSE(均方根误差,越小越好)和P值(Anderson-Darling检验,需>0.05)。短期:正态分布RMSE=0.0414、P>0.05,逻辑斯谛分布RMSE=0.0420、P>0.05;长期:正态分布RMSE=0.0425、P>0.05,逻辑斯谛分布RMSE=0.0459、P>0.05。两种模型均通过统计检验,但正态分布RMSE在急性和慢性尺度上均更优,确定为正态分布为最优拟合模型。拟合曲线显示,短期最敏感10%物种的AVE集中在19~354.7μg/L,长期最敏感10%物种的CVE集中在6.25~215.7μg/L,曲线尾部拟合良好、无异常外推风险。
HC₅定值与评估因子的数学逻辑
由正态分布模型外推各累积频率对应的物种危害浓度:短期HC₅=70.06μg/L、HC₁₀=139.54μg/L、HC₂₅=441.08μg/L、HC₅₀=1,583μg/L、HC₇₅=5,687μg/L、HC₉₀=17,979μg/L、HC₉₅=23,578μg/L;长期HC₅=8.30μg/L、HC₁₀=18.24μg/L、HC₂₅=67.90μg/L、HC₅₀=292.7μg/L、HC₇₅=1,261μg/L、HC₉₀=4,610μg/L、HC₉₅=10,318μg/L。
评估因子取值逻辑:有效毒性数据物种数>15时取2,≤15时取3。短期78个物种和长期26个物种均>15,故SAF=LAF=2。SWQC=70.06/2≈35μg/L,LWQC=8.30/2≈4.2μg/L。评估因子取2而非1的核心原因是SSD拟合基于实验室单物种数据,存在"实验室到野外"的外推不确定性,即使大样本模型也无法完全消除这种不确定性,因此保留至少2倍安全系数。
海水水质参数对SSD的潜在影响
盐度通过络合作用降低Cd²⁺相对含量——随盐度增加,Cd²⁺转化为CdCl₂和CdCl₃⁺,后者不易通过生物膜从而降低生物有效性。但文献数据无法建立统一校正模型:甲壳类(7个物种)呈"盐度降低、毒性升高"的一致趋势,鱼类(2个物种)呈"先降后升"的非单调响应。温度影响代谢速率进而改变毒性效应,但仅筛选到5个物种的24条急性数据和1个物种的3条慢性数据,不足以支撑建模。
目前全球无国家或国际组织发布海水水质参数的毒性校正技术方法指南,本次SSD拟合未进行盐度和温度校正,这是基准推导中最大的不确定性来源。未来随着更多不同盐度/温度梯度下的毒性数据积累,可能建立校正模型以提升SSD在不同海域的适用性。