水體硬度對(duì)鋅的水質(zhì)基準(zhǔn)及生態(tài)風(fēng)險(xiǎn)評(píng)估的影響

郭 淵， 楊明儒， 何 佳,2*， 秦 寧， 吳豐昌

1.中國(guó)環(huán)境科學(xué)研究院， 環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100012

2.北京師范大學(xué)水科學(xué)研究院， 北京 100875

3.北京科技大學(xué)能源與環(huán)境工程學(xué)院， 北京 100083

鋅是一種生物體所必需的微量元素，對(duì)維持生物體正常的生長(zhǎng)發(fā)育有著十分重要的作用[1-2]. 但鋅暴露超過(guò)一定劑量會(huì)對(duì)魚類[3]、溞類和藻類[4]等水生動(dòng)、植物產(chǎn)生一定的急性、慢性毒性作用，抑制動(dòng)植物的生長(zhǎng)、繁殖，破壞動(dòng)植物體內(nèi)酶的活性，甚至導(dǎo)致死亡. 因此，許多國(guó)家頒布了鋅的淡水水生生物水質(zhì)基準(zhǔn)，如美國(guó)規(guī)定鋅的長(zhǎng)期和短期水質(zhì)基準(zhǔn)(硬度為100 mg/L)均為120 μg/L[5]；加拿大頒布鋅的長(zhǎng)期水質(zhì)基準(zhǔn)值為30 μg/L[6]；澳大利亞于2000年頒布鋅的長(zhǎng)期水質(zhì)基準(zhǔn)值為8 μg/L[7]. 在不考慮硬度影響的條件下，鋅對(duì)我國(guó)本土物種的長(zhǎng)期和短期水質(zhì)基準(zhǔn)值分別為89.7～102.33 μg/L[8]和25.03～34.5 μg/L[1]. 此外，許多國(guó)家針對(duì)飲用水安全制定了鋅飲用水水質(zhì)標(biāo)準(zhǔn)，如美國(guó)規(guī)定二級(jí)飲用水中鋅濃度的標(biāo)準(zhǔn)值為5 mg/L；加拿大則規(guī)定鋅的濃度不超過(guò)5 mg/L；我國(guó)《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》(GB 5749—2006)中規(guī)定，城鄉(xiāng)各類集中式供水的生活飲用水中鋅濃度的限值為1 mg/L.

水環(huán)境因子如pH[9]、溫度[10]和硬度[11]等，對(duì)鋅的毒性會(huì)產(chǎn)生較大影響. 有研究表明，鋁[12]、鎳[13]、銀[14]和鋅[15]等重金屬對(duì)水生生物的毒性均受硬度影響，呈現(xiàn)隨硬度升高毒性降低的趨勢(shì). 已有研究[16]表明，我國(guó)典型流域地表水硬度變化范圍為27.6～1 007.30 mg/L，幾何均值為112.10 mg/L，中位數(shù)為119.88 mg/L. 該研究對(duì)鋅的水質(zhì)基準(zhǔn)進(jìn)行推導(dǎo)時(shí)充分考慮了硬度對(duì)其毒性的影響，因此對(duì)鋅的毒性數(shù)據(jù)進(jìn)行了不同梯度的硬度校正，并推導(dǎo)了相應(yīng)硬度條件下鋅保護(hù)淡水水生生物的短期和長(zhǎng)期水質(zhì)基準(zhǔn).

近年來(lái)，在黃河、長(zhǎng)江、松花江等流域中均有高含量的鋅檢出，其中，在黃河水體中鋅的暴露濃度平均值高達(dá)360.07 μg/L[17-18]，對(duì)水生態(tài)環(huán)境造成了較大的威脅[19-20]. 為了切實(shí)保護(hù)淡水水生生物的生態(tài)安全，對(duì)實(shí)際水體中鋅的生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)估十分必要[21-22]. 基于已有研究[23-24]篩選出國(guó)內(nèi)八大流域中鋅的暴露濃度，該研究在充分考慮硬度對(duì)鋅毒性影響的作用下，對(duì)鋅在國(guó)內(nèi)八大流域中的生態(tài)風(fēng)險(xiǎn)進(jìn)行了評(píng)估，以期更科學(xué)、準(zhǔn)確地評(píng)估鋅在流域不同水化學(xué)條件下的生態(tài)風(fēng)險(xiǎn)，為流域水生態(tài)安全區(qū)域化精細(xì)管理提供科學(xué)依據(jù).

1 材料與方法

1.1 毒性數(shù)據(jù)的篩選

該研究所選數(shù)據(jù)均來(lái)自Ecotox毒理數(shù)據(jù)庫(kù)(https://cfpub.epa.gov/ecotox)及文獻(xiàn)檢索數(shù)據(jù)庫(kù)(中國(guó)知網(wǎng)和Web of Science). 篩選出的受試生物主要包括我國(guó)本土物種、國(guó)際通用物種(在我國(guó)境內(nèi)水體中廣泛分布的物種)和引進(jìn)物種(具有重要經(jīng)濟(jì)價(jià)值或娛樂(lè)用途)[25]，均能反映我國(guó)淡水生物區(qū)系特征，能夠充分代表水體中不同生態(tài)營(yíng)養(yǎng)級(jí)別. 毒性數(shù)據(jù)篩選過(guò)程遵循以下原則：①剔除不含水硬度或Ca2+、Mg2+濃度的數(shù)據(jù)；②納入試驗(yàn)用水為標(biāo)準(zhǔn)稀釋水的毒性數(shù)據(jù)；③剔除未設(shè)置對(duì)照組試驗(yàn)及對(duì)照組不符的毒性數(shù)據(jù)；④同一物種同一毒性效應(yīng)測(cè)試終點(diǎn)試驗(yàn)數(shù)據(jù)相差10倍以上時(shí)，剔除離群值；⑤急性毒性試驗(yàn)暴露時(shí)間為1～4 d，終點(diǎn)為L(zhǎng)C50、EC50或IC50；⑥慢性毒性試驗(yàn)的暴露時(shí)間在21 d以上，終點(diǎn)為NOEC、LOEC等慢性終點(diǎn).

1.2 物種敏感度分布法

物種敏感度分布法(species sensitivity distribution，SSD)利用統(tǒng)計(jì)軟件(MATLAB R2016b)對(duì)毒理學(xué)數(shù)據(jù)進(jìn)行模型擬合，構(gòu)建物種的敏感度分布曲線，通過(guò)有限物種的可接受水平代表整個(gè)生態(tài)系統(tǒng)，是目前國(guó)際上較為成熟的水質(zhì)基準(zhǔn)推導(dǎo)方法之一[26]. 該基準(zhǔn)推導(dǎo)是通過(guò)將物種的急性、慢性毒性值或其對(duì)數(shù)值從小到大進(jìn)行排序，確定其毒性秩次(R)，依據(jù)式(1)計(jì)算物種的累積頻率(FR).

FR=f/(f+1)×100%

(1)

式中：FR為累積頻率，%；f為頻數(shù)，表示毒性秩次(R)對(duì)應(yīng)的物種數(shù)，個(gè).

對(duì)所獲得的累積頻率與毒性值進(jìn)行SSD模型擬合(包括正態(tài)分布、對(duì)數(shù)正態(tài)分布、邏輯斯諦分布、對(duì)數(shù)邏輯斯諦分布)，依據(jù)模型擬合的決定系數(shù)(r2)、均方根(RMSE)、殘差平方和(SSE)以及Kolmogorov-Smirnov(K-S)檢驗(yàn)系數(shù)選取最優(yōu)模型，r2越大，RMSE及SSE越小，數(shù)據(jù)擬合程度越高. 確定最優(yōu)擬合模型后計(jì)算其HC5值，除以評(píng)估因子[27-28](物種數(shù)大于15時(shí)，一般取值為2)即可確定污染物質(zhì)最終的淡水水生生物水質(zhì)基準(zhǔn).

1.3 生態(tài)風(fēng)險(xiǎn)評(píng)估

該研究所選的風(fēng)險(xiǎn)評(píng)估方法為概率密度重疊面積法和聯(lián)合概率分步法. 概率密度重疊法是將毒理數(shù)據(jù)和暴露濃度的概率密度曲線建立在同一坐標(biāo)系下，對(duì)兩條曲線重疊部分面積進(jìn)行計(jì)算，以重疊部分面積表征生物受到不利影響的概率. 重疊的面積越大，風(fēng)險(xiǎn)水平越高，反之則風(fēng)險(xiǎn)水平越小[29]. 暴露濃度概率密度曲線與毒性數(shù)據(jù)曲線交點(diǎn)表示生物耐受性水平的累積概率與暴露濃度的反累積概率之和取極小值. 交點(diǎn)的左側(cè)暴露濃度反累積概率降低較快，生物耐受性水平累積概率增加較慢，暴露濃度反累積概率與生物耐受性水平累積概率之和逐漸減小，交點(diǎn)右側(cè)則與之相反.

聯(lián)合概率分布法是將表征污染物的暴露濃度和效應(yīng)濃度的累積概率曲線建立在同一坐標(biāo)體系下，并通過(guò)概率單位轉(zhuǎn)換對(duì)概率曲線進(jìn)行直線轉(zhuǎn)換. 轉(zhuǎn)化獲得的累積概率直線上的點(diǎn)表示在特定濃度下受污染物危害的水生生物比例，暴露濃度累積概率直線上的點(diǎn)則表示暴露濃度超過(guò)此特定濃度的概率. 因此，聯(lián)合概率分布法反映了不同損害水平下暴露濃度超過(guò)安全閾值的概率[29].

通過(guò)以上兩種風(fēng)險(xiǎn)評(píng)估方法分別對(duì)各流域的鋅污染風(fēng)險(xiǎn)狀況及不同硬度條件下鋅污染風(fēng)險(xiǎn)的差別進(jìn)行分析，并采用Origin 9.1軟件繪圖.

2 結(jié)果與討論

2.1 鋅的短期水質(zhì)基準(zhǔn)

2.1.1鋅的未校正硬度短期水質(zhì)基準(zhǔn)

經(jīng)過(guò)數(shù)據(jù)篩選獲得鋅的急性毒性數(shù)據(jù)共214條，包括8門40科77個(gè)物種(見(jiàn)表1)，其中本土物種61種、國(guó)際通用物種5種以及引進(jìn)物種10種，均為中國(guó)本土及廣泛分布物種. 覆蓋魚類(一種硬骨魚綱鯉科、一種硬骨魚綱非鯉科)、甲殼類、環(huán)節(jié)類、水生植物以及上述未涉及的其他門類動(dòng)物，滿足HJ 837—2017《淡水水生生物水質(zhì)基準(zhǔn)制定技術(shù)指南》對(duì)水質(zhì)基準(zhǔn)推導(dǎo)的物種要求[41]. 鋅的毒性濃度范圍為52～70 975 μg/L.

表1 鋅對(duì)淡水水生生物的急性毒性

續(xù)表1

基于鋅對(duì)水生生物的急性毒性數(shù)據(jù)，計(jì)算種平均急性毒性值(見(jiàn)表1)，并使用4種模型進(jìn)行SSD模型擬合，通過(guò)r2、RMSE、SSE和P值(K-S檢驗(yàn))進(jìn)行擬合效果評(píng)價(jià)，結(jié)果如表2所示. 擬合最優(yōu)模型為邏輯斯蒂分布模型，鋅的急性毒性HC5計(jì)算結(jié)果為143.34 μg/L，除以評(píng)估因子后，未校正硬度條件下短期水質(zhì)基準(zhǔn)值為71.67 μg/L.

表2 未校正硬度條件下鋅的急性毒性物種敏感度分布模型參數(shù)

2.1.2鋅的硬度校正后短期水質(zhì)基準(zhǔn)

對(duì)鋅的急性毒性值和水體硬度值進(jìn)行線性擬合，獲得線性相關(guān)關(guān)系：lgCA=0.62lgH+2.06，其中l(wèi)gCA為對(duì)數(shù)急性毒性幾何均值，lgH為對(duì)數(shù)硬度值，P<0.05，具有顯著線性相關(guān)關(guān)系.

對(duì)急性毒性數(shù)據(jù)進(jìn)行水體硬度校正，分別獲得水體硬度(H)為50、100、150、200、250和350 mg/L時(shí)的毒性值，基于4種模型分別進(jìn)行SSD擬合，通過(guò)r2、RMSE、SSE、P值(K-S檢驗(yàn))的統(tǒng)計(jì)學(xué)驗(yàn)證分析(見(jiàn)表3)，發(fā)現(xiàn)邏輯斯蒂分布模型SSD曲線擬合最優(yōu). 在水體硬度為50、100、150、200、250和350 mg/L的條件下，鋅的急性毒性HC5分別為106.78、164.57、211.95、253.62、291.51和359.63 μg/L，相應(yīng)的短期水質(zhì)基準(zhǔn)值分別為53.39、82.29、105.98、126.81、145.76和179.82 μg/L.

表3 校正硬度條件下鋅的急性毒性物種敏感度分布模型參數(shù)

2.2 鋅的長(zhǎng)期水質(zhì)基準(zhǔn)

2.2.1鋅的未校正硬度長(zhǎng)期水質(zhì)基準(zhǔn)

經(jīng)過(guò)數(shù)據(jù)篩選獲得鋅對(duì)水生生物慢性數(shù)據(jù)共38條，包括6門15科20個(gè)物種，其中本土物種15種、國(guó)際通用物種4種和引進(jìn)物種1種，滿足HJ 831—2017《淡水水生生物水質(zhì)基準(zhǔn)制定技術(shù)指南》對(duì)水質(zhì)基準(zhǔn)推導(dǎo)的物種要求[41]. 鋅的毒性濃度范圍為30～2 660 μg/L.

基于鋅對(duì)水生生物的慢性毒性數(shù)據(jù)，計(jì)算種平均慢性毒性值(見(jiàn)表4)，并使用4種模型進(jìn)行SSD模型擬合(見(jiàn)表5)，結(jié)果表明正態(tài)分布模型擬合結(jié)果最優(yōu)，基于正態(tài)分布模型計(jì)算獲得鋅的慢性毒性HC5為23.80 μg/L，除以評(píng)估因子后，未校正硬度條件下長(zhǎng)期基準(zhǔn)值為11.90 μg/L.

表4 鋅對(duì)淡水生物的慢性毒性

表5 未校正水體硬度條件下鋅的慢性物種敏感度分布模型參數(shù)

2.2.2鋅的硬度校正后慢性水質(zhì)基準(zhǔn)

對(duì)鋅的慢性毒性和水體硬度進(jìn)行線性擬合，得到線性相關(guān)關(guān)系：lgCC=0.29lgH+1.37，其中l(wèi)gCC為對(duì)數(shù)慢性毒性幾何均值，P<0.05，具有顯著線性相關(guān)關(guān)系.

對(duì)慢性毒性數(shù)據(jù)進(jìn)行水體硬度校正，分別獲得水體硬度(H)為50、100、150、200、250和350 mg/L時(shí)的毒性值，基于4種模型分別進(jìn)行SSD擬合，通過(guò)r2、RMSE、SSE、P值(K-S檢驗(yàn))統(tǒng)計(jì)學(xué)驗(yàn)證分析(見(jiàn)表6)，發(fā)現(xiàn)正態(tài)分布模型為最優(yōu)擬合曲線. 在水體硬度為50、100、150、200、250和350 mg/L的條件下，鋅的慢性毒性HC5分別為20.49、25.03、28.13、30.57、32.60和34.37 μg/L，相應(yīng)的長(zhǎng)期水質(zhì)基準(zhǔn)值為10.25、12.52、14.07、15.29、16.30和17.19 μg/L.

表6 校正硬度條件下鋅的慢性毒性物種敏感度分布模型參數(shù)

比較分析該研究與其他國(guó)家或地區(qū)推導(dǎo)的鋅水質(zhì)基準(zhǔn)值發(fā)現(xiàn)，在100 mg/L硬度條件下，美國(guó)基于毒性百分?jǐn)?shù)法制定的短期和長(zhǎng)期水質(zhì)基準(zhǔn)值均為120 μg/L(長(zhǎng)期基準(zhǔn)值基于急慢性比推導(dǎo)獲得)，該研究短期基準(zhǔn)值(82.29 μg/L)相對(duì)較低，美國(guó)用于基準(zhǔn)推導(dǎo)中最敏感的4個(gè)物種(賴爪網(wǎng)紋溞、模糊網(wǎng)紋溞、條紋狼鱸和長(zhǎng)鰭鱥)鋅毒性分別為50.7、174.1、119.4和227.8 μg/L[5]，而我國(guó)鋅毒性最敏感的前4個(gè)物種(劍水溞、溞狀鉤蝦、近頭狀偽蹄型藻和羅氏沼蝦)毒性分別為52、77、80、82 μg/L，相比美國(guó)更為敏感. 在50 mg/L硬度下，加拿大基于評(píng)價(jià)因子法制定的長(zhǎng)期基準(zhǔn)值為30 μg/L，該研究長(zhǎng)期基準(zhǔn)值(10.25 μg/L)低于加拿大相同硬度水平下的基準(zhǔn)值. 在150 mg/L硬度條件下，澳大利亞和新西蘭基于物種敏感度分布法制定的長(zhǎng)期基準(zhǔn)值為8 μg/L，該研究長(zhǎng)期基準(zhǔn)值為14.07 μg/L，高于澳大利亞和新西蘭相同硬度水平下的基準(zhǔn)值. 造成這一差異主要原因是：①不同國(guó)家在進(jìn)行水質(zhì)基準(zhǔn)推導(dǎo)時(shí)各國(guó)的生物區(qū)系不盡相同; ②不同的生物區(qū)系中存在不同的敏感物種; ③使用方法、受試生物和水化學(xué)條件不同.

2.3 水體硬度對(duì)鋅水質(zhì)基準(zhǔn)的影響差異

對(duì)比硬度校正前后鋅的水質(zhì)基準(zhǔn)發(fā)現(xiàn)，除硬度校正為50 mg/L時(shí)基準(zhǔn)值低于校正前外，硬度校正為100～350 mg/L時(shí)，急性、慢性基準(zhǔn)值均大于校正前，且基準(zhǔn)值隨著硬度的升高而升高. 當(dāng)硬度校正為100 mg/L時(shí)，急性基準(zhǔn)值相比于校正前有所增加，而慢性基準(zhǔn)值校正前后差距不大；當(dāng)硬度校正為350 mg/L時(shí)，急性、慢性基準(zhǔn)值均是未校正硬度時(shí)的兩倍左右. 這說(shuō)明硬度對(duì)鋅的毒性產(chǎn)生了較大影響，水體中鋅毒性隨著硬度的增加而降低.

硬度對(duì)鋅的水生生物毒性產(chǎn)生影響的主要原因是：Ca2+、Mg2+的存在影響了金屬的化學(xué)形態(tài)進(jìn)而對(duì)金屬的毒性產(chǎn)生了影響[42]；同時(shí)，Ca2+、Mg2+也會(huì)與水生生物細(xì)胞上的金屬吸收點(diǎn)發(fā)生點(diǎn)位競(jìng)爭(zhēng). 以魚為例，Ca2+是與魚鰓上皮細(xì)胞緊密聯(lián)合的鰓膜穩(wěn)定劑，但Ca2+與Zn2+之間對(duì)于魚鰓上的金屬吸收點(diǎn)存在一定的點(diǎn)位競(jìng)爭(zhēng)，當(dāng)水體中Ca2+濃度增加時(shí)，魚鰓通過(guò)吸收外部Ca2+而降低表面滲透性，從而導(dǎo)致對(duì)Zn2+的吸收減少，因此較高的水體硬度會(huì)對(duì)鋅的毒性有一定的緩解作用[39].

除水體硬度外，不同pH、DOC濃度、溫度等條件也會(huì)對(duì)鋅的毒性產(chǎn)生影響. 如在急性試驗(yàn)條件下，當(dāng)硬度較低時(shí)，低pH(酸性)導(dǎo)致鋅的毒性上升，加速了試驗(yàn)物種死亡. 而在慢性暴露的條件下，盡管生物體未死亡，但也加重了機(jī)體受損的情況[43]. 這可能是由于鋅在不同條件下形態(tài)不同，同時(shí)也受到水體中不同的pH、DOC濃度等水質(zhì)參數(shù)的影響，導(dǎo)致生物體與其結(jié)合點(diǎn)位、方式產(chǎn)生了變化，從而導(dǎo)致毒性差異[44]. 需針對(duì)不同水化學(xué)條件下鋅的水質(zhì)基準(zhǔn)開展進(jìn)一步研究，以提高鋅的水質(zhì)標(biāo)準(zhǔn)制修訂的合理性和準(zhǔn)確性.

2.4 鋅的生態(tài)風(fēng)險(xiǎn)評(píng)估

2.4.1概率密度重疊面積法生態(tài)風(fēng)險(xiǎn)評(píng)估

首先應(yīng)用概率密度重疊面積法對(duì)鋅的生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)估，八大流域水體中鋅的暴露濃度均來(lái)源于已有研究[21]. 硬度校正前后鋅在八大流域中的急性、慢性生態(tài)風(fēng)險(xiǎn)概率密度曲線如圖1所示，計(jì)算得出的急性、慢性生態(tài)風(fēng)險(xiǎn)概率值如表7所示.

圖1 鋅在八大流域中的暴露濃度及硬度校正前后的急性、慢性毒性概率密度曲線

由表7可以看出：硬度未校正前，鋅在八大流域水體中的急性風(fēng)險(xiǎn)表現(xiàn)為黃河>長(zhǎng)江>太湖>松花江>珠江>海河>淮河>遼河，慢性風(fēng)險(xiǎn)表現(xiàn)為黃河>長(zhǎng)江>海河>太湖>珠江>松花江>淮河>遼河；硬度校正后，急性風(fēng)險(xiǎn)表現(xiàn)為黃河>長(zhǎng)江>太湖>松花江>海河>珠江>淮河>遼河，慢性風(fēng)險(xiǎn)表現(xiàn)為黃河>長(zhǎng)江>海河>珠江>太湖>松花江>淮河>遼河.

表7 八大流域水生生物鋅污染的急性、慢性風(fēng)險(xiǎn)概率

未進(jìn)行硬度校正情況下，除黃河受鋅污染的急性生態(tài)風(fēng)險(xiǎn)概率超過(guò)40%外，其余流域急性風(fēng)險(xiǎn)均低于20%，其中遼河最低，僅為2.7%. 黃河受鋅污染的慢性生態(tài)風(fēng)險(xiǎn)影響程度高達(dá)74.4%，長(zhǎng)江、海河發(fā)生慢性鋅污染風(fēng)險(xiǎn)概率也超過(guò)了40%. 當(dāng)水體硬度校正為50mg/L時(shí)，各流域鋅污染急性、慢性風(fēng)險(xiǎn)概率相比于校正前均有所升高，其中海河慢性風(fēng)險(xiǎn)從校正前的40.64%升至校正后的45%，為最大漲幅，其余流域漲幅均低于4.5%. 當(dāng)水體硬度條件校正為100～350 mg/L時(shí)，各流域急性、慢性生態(tài)風(fēng)險(xiǎn)概率相比于校正前均有所降低，且隨硬度升高風(fēng)險(xiǎn)概率呈依次遞減趨勢(shì). 急性生態(tài)風(fēng)險(xiǎn)概率隨硬度增大明顯下降，例如，在水體硬度校正為350 mg/L時(shí)，淮河、遼河和海河水體中鋅的生態(tài)風(fēng)險(xiǎn)概率相比于校正前最大降幅均在200%以上，太湖、長(zhǎng)江、珠江和松花江的風(fēng)險(xiǎn)概率降幅均接近200%. 相比于急性風(fēng)險(xiǎn)，硬度校正后慢性風(fēng)險(xiǎn)的降幅較小. 例如，在水體硬度校正為100 mg/L時(shí)，各流域慢性風(fēng)險(xiǎn)相比于校正前未發(fā)現(xiàn)顯著性差異，校正前松花江慢性風(fēng)險(xiǎn)為35.61%，校正后僅為35.10%. 但在硬度校正為350 mg/L的條件下，海河慢性風(fēng)險(xiǎn)相比校正前降幅最大(8.6%)，長(zhǎng)江、珠江、松花江、淮河和遼河降幅均在7.0%左右，黃河和太湖降幅在6.0%左右. 研究表明，鋅在各流域水體中的急性、慢性生態(tài)風(fēng)險(xiǎn)均受到水體硬度的顯著影響.

2.4.2聯(lián)合概率分布法生態(tài)風(fēng)險(xiǎn)評(píng)估

概率密度函數(shù)重疊面積表征方法得到的結(jié)果直接反映了污染物在各流域水體中的生態(tài)風(fēng)險(xiǎn)狀況，聯(lián)合概率分布法則可以根據(jù)受影響的物種比例來(lái)確定保護(hù)水平. 由于聯(lián)合概率分步法充分考慮了環(huán)境暴露濃度和毒性數(shù)據(jù)的不確定性和可變性. 因此選擇應(yīng)用聯(lián)合概率分布法，計(jì)算了鋅在八大流域水體中的暴露濃度以及不同硬度條件下超過(guò)保護(hù)95%生物不受影響的急性和慢性HC5值的超標(biāo)率. 如圖2所示，暴露濃度對(duì)應(yīng)的累積概率值越小，超標(biāo)情況越嚴(yán)重，鋅對(duì)水生生物造成的生態(tài)風(fēng)險(xiǎn)越嚴(yán)重.

圖2 鋅在八大流域中的暴露濃度及硬度校正前后的急性、慢性毒性聯(lián)合概率分布

如表8所示：硬度校正前，鋅在各流域中對(duì)急性HC5超標(biāo)率表現(xiàn)為黃河>長(zhǎng)江>太湖>海河>松花江>珠江>淮河>遼河，慢性HC5超標(biāo)率表現(xiàn)為遼河>黃河>長(zhǎng)江>海河>珠江>松花江>太湖>淮河；硬度校正后，鋅在各流域中對(duì)急性HC5超標(biāo)率表現(xiàn)為黃河>長(zhǎng)江>太湖>海河>珠江>松花江>淮河>遼河，慢性HC5超標(biāo)率表現(xiàn)為遼河>黃河>海河>長(zhǎng)江>淮河>珠江>松花江>太湖.

表8 聯(lián)合概率分布法評(píng)估鋅在水體中暴露濃度對(duì)HC5的超標(biāo)率

未校正硬度條件下，黃河、遼河、黃河、長(zhǎng)江和海河水體中鋅的慢性HC5超標(biāo)率均高于50%，其中，遼河不存在急性HC5超標(biāo)情況，但慢性HC5超標(biāo)率在70%以上；而黃河水體中鋅的急性、慢性超標(biāo)率均相對(duì)較高，急性和慢性超標(biāo)率分別為42%和71%. 當(dāng)水體硬度校正為50 mg/L時(shí)，各流域水體中鋅的急性、慢性超標(biāo)率均有所上升，其中黃河慢性HC5超標(biāo)率高達(dá)73%. 當(dāng)水體硬度校正為100～350 mg/L時(shí)，各流域的急性和慢性HC5超標(biāo)率則隨著硬度的增加而呈降低趨勢(shì). 其中，當(dāng)水體硬度校正為350 mg/L時(shí)，太湖水體中鋅的急性HC5超標(biāo)率降低了50%以上，長(zhǎng)江的急性超標(biāo)率降低了75%以上，而珠江、松花江、海河的急性超標(biāo)率均從6%降至1%，降幅達(dá)600%. 鋅在各流域水體中的慢性超標(biāo)率情況類似，其中遼河的慢性超標(biāo)率降低了60%左右.

綜上，通過(guò)兩種生態(tài)風(fēng)險(xiǎn)評(píng)估方法所獲得的結(jié)論基本一致，即當(dāng)水體硬度為50 mg/L時(shí)，相比于硬度校正前，各流域水體中發(fā)生急性、慢性鋅污染的生態(tài)風(fēng)險(xiǎn)概率和HC5超標(biāo)率均有所升高；當(dāng)水體硬度為100～350 mg/L時(shí)，相較于硬度校正前，各流域水體中發(fā)生急性、慢性鋅污染的生態(tài)風(fēng)險(xiǎn)概率和HC5超標(biāo)率均有所降低，并且隨著水體硬度的增加而降低. 以上研究表明，開展不同流域水體生態(tài)風(fēng)險(xiǎn)評(píng)估時(shí)，應(yīng)充分考慮不同流域水體硬度的梯度變化，科學(xué)、合理、有效地保護(hù)流域生態(tài)系統(tǒng)安全.

3 結(jié)論

a) 通過(guò)物種敏感度分布法推導(dǎo)了我國(guó)淡水水生生物在水體硬度校正前的長(zhǎng)期、短期水質(zhì)基準(zhǔn)值分別為11.90和71.67 μg/L. 在水體硬度校正為50、100、150、200、250和350 mg/L的條件下，鋅的短期基準(zhǔn)值分別為53.39、82.29、105.98、126.81、145.76和179.82 μg/L，鋅的長(zhǎng)期基準(zhǔn)值分別為10.25、12.52、14.07、15.29、16.30和17.19 μg/L. 水體硬度對(duì)水質(zhì)基準(zhǔn)具有顯著性影響，建議在鋅的水質(zhì)基準(zhǔn)和標(biāo)準(zhǔn)制修訂過(guò)程中，應(yīng)充分考慮水體硬度對(duì)淡水水生生物毒性的影響.

b) 八大流域水體受鋅污染的急性、慢性生態(tài)風(fēng)險(xiǎn)和HC5超標(biāo)率均受到水體硬度的顯著影響，生態(tài)風(fēng)險(xiǎn)和HC5超標(biāo)率隨水體硬度的升高均呈逐步下降趨勢(shì). 考慮到我國(guó)水體流域范圍大，水體硬度南北差別較大，建議基于不同水體硬度條件開展區(qū)域性生態(tài)風(fēng)險(xiǎn)評(píng)估，為流域水生態(tài)系統(tǒng)的精細(xì)化管理提供技術(shù)支撐.