摘 要：為了解萬山汞礦區(qū)彈琴蛙（Nidirana adenopleura）組織重金屬富集特征及其關(guān)聯(lián)，對(duì)萬山汞礦區(qū)典型農(nóng)田生境內(nèi)彈琴蛙組織（肝臟、皮膚、股骨、肌肉和趾骨）中鉻（Cr）、錳（Mn）、鎳（Ni）、銅（Cu）、鋅（Zn）、鎘（Cd）、鉛（Pb）、砷（As）、硒（Se）和汞（Hg）這10種重金屬元素的生物富集特征差異進(jìn)行分析，結(jié)合非負(fù)矩陣分解分析法識(shí)別重金屬元素的來源及貢獻(xiàn)，并通過灰色關(guān)聯(lián)度分析法分析彈琴蛙組織重金屬元素間的關(guān)聯(lián)性。結(jié)果表明：重金屬元素在萬山汞礦區(qū)彈琴蛙體內(nèi)富集程度不同，平均質(zhì)量分?jǐn)?shù)由大到小依次為Zn、Cu、Mn、Pb、Se、Cr、Hg、Cd、Ni、As；Cu、Se、Cr、Ni、Cd和Hg主要富集在肝臟組織中，Mn、Zn和Pb主要富集在股骨和趾骨組織中，As主要富集在肌肉組織中。基于非負(fù)矩陣分解和相關(guān)分析，Hg、Cu、Cr、Pb 和Se主要源于汞礦開采，As主要來源于農(nóng)業(yè)活動(dòng)，Mn主要來源于錳礦冶煉，Ni、Zn和Cd來源于多種復(fù)合途徑?？傊?，兩棲動(dòng)物趾骨具有較好的生物指示作用，可作為一種非致命性的檢測(cè)方法，用于動(dòng)物重金屬污染的研究。研究結(jié)果可為萬山汞礦區(qū)環(huán)境污染防控和野生動(dòng)物保護(hù)提供科學(xué)依據(jù)。

關(guān)鍵詞：重金屬；生物富集；彈琴蛙；汞礦

中圖分類號(hào)：Q89；X171

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：2310 - 1490（2024）- 03 - 0596 - 11

DOI：10.12375/ysdwxb.20240316

重金屬對(duì)環(huán)境的污染被認(rèn)為是生態(tài)學(xué)研究的關(guān)鍵問題之一［1］，會(huì)對(duì)野生動(dòng)物及環(huán)境造成嚴(yán)重危害［2］，威脅生態(tài)安全。采礦活動(dòng)是環(huán)境重金屬最重要的人為污染源之一，會(huì)導(dǎo)致周邊生態(tài)系統(tǒng)有較高的污染風(fēng)險(xiǎn)［3］。萬山汞礦曾是我國最大的汞礦，已有600多年的開采歷史［4］，已于2004年關(guān)閉［5］。大規(guī)模的汞礦開采和冶煉活動(dòng)已對(duì)周邊河流［6］、土壤［4］、植被［7］和大氣［8］等多個(gè)環(huán)境領(lǐng)域造成不同程度的重金屬污染，而域內(nèi)產(chǎn)的大米是居民健康風(fēng)險(xiǎn)的主要來源［9］。

目前，關(guān)于野生動(dòng)物重金屬污染的研究集中在兩棲動(dòng)物方面，特別是無尾兩棲動(dòng)物，這是因?yàn)樗鼈儗?duì)環(huán)境變化較為敏感，是水生和陸生生態(tài)系統(tǒng)食物鏈的關(guān)鍵組成部分［10］，環(huán)境污染被認(rèn)為是全球兩棲動(dòng)物數(shù)量下降的主要原因［11］。兩棲動(dòng)物具有特殊的生活史及高滲透性皮膚，環(huán)境重金屬等污染物可通過皮膚、呼吸和食物進(jìn)入其體內(nèi)，導(dǎo)致其更易受到環(huán)境重金屬的污染［3］。這些因素表明，兩棲動(dòng)物可能是較好的指示生物，越來越多的研究利用兩棲動(dòng)物進(jìn)行重金屬等污染物的監(jiān)測(cè)［3，12］。

彈琴蛙（Nidirana adenopleura）屬蛙科（Ranidae）琴蛙屬（Nidirana），廣泛分布于我國重慶、云南、貴州、安徽和浙江等南方地區(qū)，多生活在農(nóng)作區(qū)及其附近水塘內(nèi)，種群數(shù)量多［13］。農(nóng)業(yè)活動(dòng)同樣也是環(huán)境重金屬的污染來源之一，農(nóng)田長期使用農(nóng)藥、殺蟲劑和除草劑，會(huì)導(dǎo)致銅（Cu）、鋅（Zn）、鎳（Ni）和鎘（Cd）等多種重金屬的積累［5］。已有研究表明，萬山汞礦區(qū)農(nóng)田土壤重金屬為重度污染［14］，農(nóng)產(chǎn)品受到不同程度的重金屬污染［9］，但在該區(qū)域重金屬開采造成污染的背景下，關(guān)于農(nóng)業(yè)活動(dòng)等人為因素導(dǎo)致的重金屬復(fù)合污染生境中的野生動(dòng)物重金屬復(fù)合污染的研究較少。因此，本研究以萬山汞礦區(qū)農(nóng)田生境中的彈琴蛙為研究對(duì)象，在分析彈琴蛙組織重金屬元素富集特征的基礎(chǔ)上，結(jié)合非負(fù)矩陣分解（nonnegativematrix factorization，NMF）分析法來識(shí)別重金屬元素的來源及貢獻(xiàn)，并通過灰色關(guān)聯(lián)度分析法（greyrelational analysis，GRA）分析彈琴蛙組織重金屬的關(guān)聯(lián)性，以探討兩棲動(dòng)物趾骨作為重金屬污染的生物指示作用。

1 研究區(qū)概況

萬山汞礦位于貴州省銅仁市東南部，呈現(xiàn)出亞熱帶濕潤季風(fēng)氣候［4］。主要河流下溪河和敖寨河均源于汞礦區(qū)，下游并入瓦屋河［5］。萬山汞礦地勢(shì)較高，且處于喀斯特溶巖地區(qū)，重金屬污染物仍可通過徑流、降水和風(fēng)等方式輸送到周邊區(qū)域［3］，對(duì)下游農(nóng)田構(gòu)成潛在環(huán)境污染風(fēng)險(xiǎn)［6］。瓦屋侗族鄉(xiāng)位于萬山汞礦東北部下游區(qū)域，境內(nèi)主要河流為瓦屋河，是重要的農(nóng)業(yè)區(qū)，其灌溉水主要由河流提供［4］。本研究選擇瓦屋侗族鄉(xiāng)典型稻田生境作為彈琴蛙的取樣區(qū)。

2 材料與方法

2. 1 樣品采集與處理

在瓦屋侗族鄉(xiāng)周邊選擇3處稻田作為彈琴蛙取樣區(qū)，在每處取樣區(qū)隨機(jī)捕捉體型相近的雄性彈琴蛙成體（體長55 mm左右）［13］各5只，共捕捉彈琴蛙15只，具體取樣位置見圖1。捕捉的彈琴蛙個(gè)體平均體長為（55. 50±0. 58） mm，平均體質(zhì)量為（12. 37±0. 80） g，帶回實(shí)驗(yàn)室解剖并取其肝臟、皮膚、股骨和肌肉，依次用自來水、去離子水充分洗滌，放入烘箱烘干至恒質(zhì)量，研磨粉碎待用。趾骨選擇右后肢第5根腳趾，用蒸餾水沖洗后，剪刀去除皮膚和肌肉組織，然后使用30%的過氧化氫清洗結(jié)膜組織，再用去離子水充分洗滌趾骨樣品，最后放入烘箱烘干至恒質(zhì)量［15］。本研究通過了銅仁學(xué)院科學(xué)研究倫理委員會(huì)的審查（審批同意書編號(hào)：TU20231005）。

2. 2 樣品測(cè)定與分析

根據(jù)萬山汞礦區(qū)環(huán)境重金屬污染特征，選擇10種重金屬（類金屬）元素進(jìn)行分析。鉻（Cr）、錳（Mn）、Ni、Cu、Zn、Cd、鉛（Pb）和砷（As）采用電感耦合等離子體質(zhì)譜法（inductively coupled plasma mass spectrometry，ICP-MS）測(cè)定，硒（Se）和汞（Hg）采用原子熒光光譜法（atomic fluorescence spectrometry，AFS）測(cè)定。采用生物標(biāo)準(zhǔn)物質(zhì)（扇貝，GB10024）監(jiān)測(cè)實(shí)驗(yàn)過程，每組樣品設(shè)置空白對(duì)照樣，相對(duì)標(biāo)準(zhǔn)偏差控制在5%以內(nèi)。

2. 3 生物富集系數(shù)

生物富集系數(shù)（bioaccumulation factor，BAF）是生物體或其特定組織中的化學(xué)化合物濃度與其在周圍環(huán)境中的濃度比值［15］，可反映生物體對(duì)環(huán)境重金屬的富集能力［16］，計(jì)算公式

BAF = Cb /Cf 。

式中：Cb為彈琴蛙各組織中的重金屬含量，Cf為環(huán)境中的土壤重金屬含量。在本研究中，組織和土壤重金屬均以干質(zhì)量計(jì)，mg/kg。如果BAFgt;1. 00，則表明具有潛在的富集風(fēng)險(xiǎn)。

2. 4 土壤重金屬背景值

萬山汞礦區(qū)土壤重金屬污染已有較多研究，本研究以最新研究成果作為研究區(qū)土壤環(huán)境重金屬背景值。Cr、Mn、Ni、Cu、Zn、Cd、Pb、As和Hg的質(zhì)量分?jǐn)?shù)來源于Liu 等［4］的研究，Se 的質(zhì)量分?jǐn)?shù)來源于Chang等［17］的研究（表1）。

2. 5 重金屬來源分析模型

利用R 4. 3. 1軟件中的NMF包計(jì)算并分析彈琴蛙體內(nèi)重金屬元素的來源及貢獻(xiàn)，分別對(duì)設(shè)置為3～8個(gè)因子數(shù)量［18］的模型進(jìn)行分析，并根據(jù)殘差減少率和因子解釋程度確定最適合該模型的因子數(shù)量。最終，本研究采用3因子模型的計(jì)算結(jié)果。

2. 6 灰色關(guān)聯(lián)度分析（GRA）

GRA是基于灰色系統(tǒng)的多因素分析方法［19］，能夠描述各因素關(guān)系的強(qiáng)弱、大小和次序［20?21］。本研究通過計(jì)算灰色關(guān)聯(lián)度來分析彈琴蛙趾骨與肝臟、股骨、肌肉和皮膚重金屬元素含量的關(guān)聯(lián)程度。關(guān)聯(lián)度值越接近1. 00，關(guān)聯(lián)度越強(qiáng)，當(dāng)該值≥0. 50時(shí)，通常表明關(guān)聯(lián)度顯著［21］。

2. 7 數(shù)據(jù)處理與統(tǒng)計(jì)分析

在彈琴蛙組織重金屬檢測(cè)過程中，部分樣品在檢測(cè)限下未被檢出。針對(duì)這種情況，若某種組織50%以上的樣品檢測(cè)結(jié)果高于檢測(cè)限，則被認(rèn)為是“檢測(cè)到”。濃度值低于檢測(cè)限，但認(rèn)為是“檢測(cè)到”的組織樣品，被替換為每種元素檢測(cè)限值的一半，而檢測(cè)頻率低于50%的結(jié)果不進(jìn)行多變量統(tǒng)計(jì)分析［22］。本研究中的部分?jǐn)?shù)據(jù)不符合正態(tài)分布，利用非參數(shù)Kruskal-Wallis檢驗(yàn)法比較重金屬元素含量的差異。采用Spearman相關(guān)秩次檢驗(yàn)進(jìn)行組織間、元素間的相關(guān)性分析，以plt;0. 05為有顯著差異。以上數(shù)據(jù)分析在SPSS 27. 0軟件中進(jìn)行，采用Origin 2022作圖。

3 結(jié)果與分析

3. 1 組織重金屬質(zhì)量分?jǐn)?shù)及分布

通過分析彈琴蛙組織樣品的檢測(cè)頻率，發(fā)現(xiàn)Cr、Mn、Ni、Cu、Zn、Cd、Pb和Hg這8種重金屬元素在肝臟、股骨、肌肉、皮膚和趾骨內(nèi)均被檢測(cè)到，As在肝臟、股骨、肌肉和趾骨內(nèi)被檢測(cè)到，Se在肝臟、股骨、肌肉和皮膚內(nèi)被檢測(cè)到（表2）。

重金屬元素在萬山汞礦區(qū)彈琴蛙體內(nèi)蓄積程度不同，彈琴蛙體內(nèi)重金屬平均質(zhì)量分?jǐn)?shù)由大到小依次為Zn、Cu、Mn、Pb、Se、Cr、Hg、Cd、Ni、As。Zn在彈琴蛙體內(nèi)平均質(zhì)量分?jǐn)?shù)最高，達(dá)（581. 75±8. 36）mg/kg，Cu 和Mn 次之，平均質(zhì)量分?jǐn)?shù)分別為（384. 58±18. 68）、（308. 51±19. 41）mg/kg。在非必需元素中，Pb平均質(zhì)量分?jǐn)?shù)相對(duì)較高，為（14. 29±1. 53）mg/kg，可能對(duì)彈琴蛙的健康產(chǎn)生一定影響。As的平均質(zhì)量分?jǐn)?shù)在所有元素中相對(duì)最低，僅為（0. 08±0. 01）mg/kg（ 圖2）。

通過對(duì)不同組織重金屬元素平均質(zhì)量分?jǐn)?shù)的分析發(fā)現(xiàn)，重金屬元素在彈琴蛙不同組織間的質(zhì)量分?jǐn)?shù)均存在顯著差異（plt;0. 05），肝臟是重金屬元素的主要蓄積部位。Cu、Se、Cr、Ni、Cd和Hg六種元素主要蓄積在肝臟組織中，且質(zhì)量分?jǐn)?shù)顯著高于其他組織（plt;0. 05）。Mn、Zn和Pb三種元素主要蓄積在股骨和趾骨中，且質(zhì)量分?jǐn)?shù)顯著高于其他組織（plt;0. 05）。As主要蓄積在肌肉組織中，且質(zhì)量分?jǐn)?shù)顯著高于其他組織（plt;0. 05）（圖3）。

3. 2 重金屬生物富集特征

為評(píng)估彈琴蛙對(duì)環(huán)境重金屬元素的生物富集特征及富集風(fēng)險(xiǎn)，分別計(jì)算其各組織對(duì)土壤中不同重金屬元素的生物富集系數(shù)（圖4）。彈琴蛙對(duì)土壤中重金屬元素的富集作用較低，各組織對(duì)土壤中Cr、Mn、Ni、As、Pb 和Hg 六種元素的富集系數(shù)均小于1. 00，表明具有較低的富集風(fēng)險(xiǎn)。肝臟對(duì)土壤中Cu、Se和Cd三種元素的富集風(fēng)險(xiǎn)較高，特別是Cu的富集系數(shù)達(dá)14. 45。彈琴蛙股骨、皮膚和趾骨對(duì)Zn存在一定的富集作用，但富集風(fēng)險(xiǎn)相對(duì)較低，富集系數(shù)均在1. 50以下。

3. 3 重金屬來源分析

為分析彈琴蛙體內(nèi)重金屬元素的來源，對(duì)10種重金屬元素的含量進(jìn)行Spearman 相關(guān)性分析（圖5）。當(dāng)元素間存在顯著正相關(guān)關(guān)系時(shí)，可能具有相同或者相似的來源。Hg—Pb（0. 59）和Hg—Se（0. 52）呈顯著正相關(guān)關(guān)系（plt;0. 05），可能Pb和Se與萬山汞礦開采有關(guān)。Pb—Cd（0. 56）、Pb—Ni（0. 51）和Cu—Cr（0. 66）呈顯著正相關(guān)關(guān)系（plt;0. 05），可能這些重金屬元素存在共同或相似的來源。As—Cu（-0. 69）呈顯著負(fù)相關(guān)關(guān)系（plt;0. 05），可能這兩種元素環(huán)境來源不同。

基于NMF進(jìn)一步對(duì)彈琴蛙體內(nèi)重金屬元素的環(huán)境來源進(jìn)行解析，最適合模型的因子數(shù)量為3個(gè)，各重金屬元素的來源占比情況見圖6。

因子1的貢獻(xiàn)率為18. 67%，載荷較高的重金屬元素為Hg（70. 87%）、Cu（67. 97%）、Cr（66. 85%）、Pb（55. 89%）和Se（51. 12%）。因子2 的貢獻(xiàn)率為32. 76%，載荷較高的重金屬元素為As（50. 35%）。因子3的貢獻(xiàn)率為48. 57%，載荷較高的重金屬元素為Mn（61. 92%）。Ni、Zn和Cd在3種因子下的載荷均不高，可能這3種元素具有不同或者多種來源。

3. 4 組織重金屬的關(guān)聯(lián)性

除Se以外，其他9種重金屬元素在趾骨組織中均被檢測(cè)到，對(duì)趾骨與其他4種組織重金屬元素間的相關(guān)性進(jìn)行分析，以了解趾骨與其他組織重金屬元素的關(guān)聯(lián)性（圖7）。趾骨與其他組織重金屬元素的相關(guān)系數(shù)較大，且絕大部分為正相關(guān)關(guān)系。趾骨與肝臟組織中As和Hg的含量、與股骨組織中Mn、Pb和Hg的含量、與肌肉組織中Cr和Hg的含量以及與皮膚組織中Hg 的含量均呈顯著正相關(guān)關(guān)系（plt;0. 05）。

進(jìn)一步利用GRA 分析彈琴蛙趾骨與肝臟、股骨、肌肉和皮膚重金屬元素含量的關(guān)聯(lián)程度（圖8）。趾骨中Se未被檢測(cè)到，皮膚中As未被檢測(cè)到，未進(jìn)行分析。趾骨與其他組織重金屬元素含量關(guān)聯(lián)程度較高、關(guān)聯(lián)性較強(qiáng)，關(guān)聯(lián)度均大于0. 50。趾骨與肝臟中Cd的關(guān)聯(lián)度相對(duì)最高，為0. 88；趾骨與股骨中Cu的關(guān)聯(lián)度相對(duì)最低，為0. 60。

4 討論

4. 1 彈琴蛙重金屬富集特征

重金屬元素在彈琴蛙體內(nèi)富集特征存在一定差異，Zn平均質(zhì)量分?jǐn)?shù)相對(duì)最高，動(dòng)物體內(nèi)Zn是僅次于Fe的第二大微量元素，這與已有的兩棲動(dòng)物研究結(jié)果［23］一致，可能與萬山汞礦區(qū)土壤中Zn屬中度污染［24］有關(guān)。非必需元素Pb的平均質(zhì)量分?jǐn)?shù)相對(duì)較高，但彈琴蛙對(duì)土壤中Pb的富集風(fēng)險(xiǎn)較低，Pb可能來源于其他環(huán)境介質(zhì)。Hg作為萬山汞礦區(qū)土壤主要的重金屬污染物，雖存在極強(qiáng)的生態(tài)風(fēng)險(xiǎn)［25］，但對(duì)彈琴蛙的富集風(fēng)險(xiǎn)較低，這可能與區(qū)域內(nèi)Hg的大氣沉降來源有關(guān)［26］，本研究中彈琴蛙的取樣區(qū)距萬山汞礦存在一定距離，對(duì)其影響相對(duì)較小。As在彈琴蛙各組織內(nèi)的質(zhì)量分?jǐn)?shù)均相對(duì)較低，且具有較低的富集風(fēng)險(xiǎn)，這可能是萬山汞礦區(qū)As污染程度不重，屬清潔狀態(tài)［26］所致。

肝臟在有毒物質(zhì)解毒中發(fā)揮核心作用［27］，重金屬元素易富集于肝臟［28］。在本研究中，肝臟是彈琴蛙重金屬元素的主要蓄積部位，導(dǎo)致肝臟對(duì)土壤中Cu、Se和Cd這3種元素具有較高的富集風(fēng)險(xiǎn)，但是重金屬元素在各組織中的富集程度也存在一定差異［29］。Mn、Zn和Pb這3種元素主要富集在彈琴蛙的股骨和趾骨中，這是由于Zn對(duì)骨骼生長具有重要作用，骨作為鈣化組織［30］對(duì)Pb的富集作用較強(qiáng)［31］、對(duì)Mn的富集系數(shù)較高［32］所致。As主要富集在肌肉組織中，這與澤陸蛙（Fejervarya multistriata）重金屬富集的研究結(jié)果［23］一致。

4. 2 重金屬污染源解析

萬山汞礦周邊重金屬污染問題突出［26］。在本研究中，彈琴蛙取樣區(qū)位于萬山汞礦區(qū)下游的瓦屋侗族鄉(xiāng)周邊，重金屬元素由于水力傳輸、大氣沉降等的原因輸入，導(dǎo)致該區(qū)仍受到多種重金屬的污染［4］。該區(qū)汞礦石表現(xiàn)出高Cu、Hg和Pb的特點(diǎn)［33］，Se是汞礦床礦物基質(zhì)中重要的共存元素［34］。因此，基于NMF和相關(guān)分析結(jié)果，Hg、Cu、Cr、Pb和Se這5種元素具有較高的同源性，可能主要源于萬山悠久的汞礦開采活動(dòng)。As受人為干擾影響大［35］，彈琴蛙個(gè)體采集于稻田，農(nóng)藥、殺蟲劑和除草劑的長期使用易造成As的積累［36］，As可能主要來源于農(nóng)業(yè)活動(dòng)。司前大壩周邊的錳礦冶煉導(dǎo)致瓦屋侗族鄉(xiāng)的Mn污染，這一區(qū)域農(nóng)產(chǎn)品中Mn的危害指數(shù)高于Hg［9］，導(dǎo)致Mn的貢獻(xiàn)率較高，Mn可能主要來源于錳礦冶煉活動(dòng)。Ni、Cd 與Pb 均呈顯著正相關(guān)關(guān)系（plt;0. 05），但在NMF中，Ni、Zn和Cd這3種元素在3種因子下的載荷均不高，這3種元素可能來源于多種復(fù)合途徑。

4. 3 趾骨的生物指示作用

兩棲動(dòng)物在水生生態(tài)系統(tǒng)中分布廣泛，在許多研究中被認(rèn)為是各種污染物的指示生物［1，37］。在本研究中，10種重金屬元素在彈琴蛙絕大部分組織中均被檢測(cè)到，表明作為區(qū)域優(yōu)勢(shì)種的彈琴蛙同樣可作為監(jiān)測(cè)重金屬污染的指示生物。

目前，在兩棲動(dòng)物重金屬污染的研究中多采用組織破壞性取樣技術(shù)，如取動(dòng)物的肝臟、腎臟和肌肉等［32，38?39］，但這些取樣方法均需要以殺死動(dòng)物為代價(jià)，這與動(dòng)物保護(hù)目標(biāo)相沖突［40］。本研究發(fā)現(xiàn)，彈琴蛙趾骨與肝臟、股骨、肌肉和皮膚重金屬元素含量絕大部分呈正相關(guān)關(guān)系，趾骨與其他組織的關(guān)聯(lián)程度較高、關(guān)聯(lián)性較強(qiáng)，表明趾骨具有較好的生物指示作用。因此，在今后的研究中沒必要為監(jiān)測(cè)重金屬污染而殺死這些兩棲動(dòng)物，以兩棲動(dòng)物趾骨為重金屬分析材料，用于動(dòng)物重金屬污染的研究，是簡單、安全和非致命性的檢測(cè)方法，趾骨可作為評(píng)估區(qū)域內(nèi)兩棲動(dòng)物體內(nèi)主要重金屬富集特征及環(huán)境污染的生物指示物。

5 結(jié)論

（1）重金屬元素在萬山汞礦區(qū)彈琴蛙體內(nèi)的富集程度不同，其平均質(zhì)量分?jǐn)?shù)由大到小依次為Zn、Cu、Mn、Pb、Se、Cr、Hg、Cd、Ni和As。

（2）Cu、Se、Cr、Ni、Cd和Hg主要富集在肝臟組織中，Mn、Zn和Pb主要富集在股骨和趾骨組織中，As主要富集在肌肉組織中。肝臟對(duì)土壤中Cu、Se和Cd這3種元素的富集風(fēng)險(xiǎn)較高，對(duì)其他重金屬元素的富集風(fēng)險(xiǎn)相對(duì)較低。

（3）基于NMF 和相關(guān)分析，Hg、Cu、Cr、Pb 和Se主要源于汞礦開采，As主要來源于農(nóng)業(yè)活動(dòng)，Mn主要來源于錳礦冶煉，Ni、Zn 和Cd 來源于多種復(fù)合途徑。

（4）兩棲動(dòng)物趾骨具有較好的生物指示作用，可作為一種非致命性的檢測(cè)樣本。

參考文獻(xiàn)：

［1］ CHEBOTINA M Y， GUSEVA V P， BERZIN D L， et al. Accumulation

of trace elements in the marsh frog Pelophylax ridibundus in

cooling ponds of the middle Urals［J］. Inland Water Biology，

2022， 15（2）： 189-194.

［2］ 王祥生， 婁虎， 張明海. 野生動(dòng)物重金屬污染研究進(jìn)展［J］. 生

命科學(xué)研究， 2019， 23（1）： 78-86.

WANG X S， LOU H， ZHANG M H. Advances in heavy metal pollution

in wildlife［J］. Life Science Research， 2019， 23（1）： 78-

86.

［3］ MANI M， ALTUNI?IK A， GEDIK K. Bioaccumulation of trace elements

and health risk predictions in edible tissues of the marsh

frog［J］. Biological Trace Element Research， 2022， 200（10）：

4493-4504.

［4］ LIU B L， TIAN K， HE Y， et al. Dominant roles of torrential

floods and atmospheric deposition revealed by quantitative source

apportionment of potentially toxic elements in agricultural soils

around a historical mercury mine， southwest China［J］. Ecotoxi‐

cology and Environmental Safety， 2022， 242： 113854.

［5］ ADLANE B， XU Z D， XU X H， et al. Evaluation of the potential

risks of heavy metal contamination in rice paddy soils around an

abandoned Hg mine area in southwest China［J］. Acta Geochimica，

2020， 39（1）： 85-95.

［6］ 蔡敬怡， 譚科艷， 路國慧， 等. 貴州萬山廢棄礦區(qū)小流域系統(tǒng)

沉積物及懸浮物重金屬的空間分布特征［J］. 巖礦測(cè)試，

2019， 38（3）： 305-315.

CAI J Y， TAN K Y， LU G H， et al. The spatial distribution characteristics

of heavy metals in river sediments and suspended matter

in small tributaries of the abandoned Wanshan mercury mines，

Guizhou Province［J］. Rock and Mineral Analysis， 2019， 38（3）：

305-315.

［7］ TANG X C， WU X L， XIA P H， et al. Health risk assessment of

heavy metals in soils and screening of accumulating plants around

the Wanshan mercury mine in northeast Guizhou Province， China

［J］. Environmental Science and Pollution Research， 2021， 28

（35）： 48837-48850.

［8］ LIAN M， SHANG L H， DUAN Z， et al. Lidar mapping of atmospheric

atomic mercury in the Wanshan area， China［J］. Environmental

Pollution， 2018， 240： 353-358.

［9］ 湛天麗， 黃陽， 何騰兵， 等. 貴州銅仁汞礦區(qū)主要農(nóng)產(chǎn)品重金

屬污染及其健康風(fēng)險(xiǎn)評(píng)估［J］. 安全與環(huán)境學(xué)報(bào)， 2017， 17

（4）： 1524-1529.

ZHAN T L， HUANG Y， HE T B， et al. Heavy metal pollution

and health risk in major agricultural products in Tongren mercury

mine area， Guizhou［J］. Journal of Safety and Environment，

2017， 17（4）： 1524-1529.

［10］ CRUZ-SANTIAGO O， PéREZ-MALDONADO I N， GONZáLEZMILLE

D J， et al. Nondestructive biomarkers in giant toad （Rhi?

nella marina） to assess the effect of complex mixture of pollutants

in Coatzacoalcos River， Mexico［J］. Environmental Toxicology

and Pharmacology， 2021， 82： 103558.

［11］ SANTANA E， SCHIESARI L， GOMES F， et al. Morphophysiological

traits of an amphibian exposed to historical industrial pollution

in a Brazilian biodiversity hotspot［J］. Amphibia-Reptilia，

2021， 42（3）： 283-295.

［12］ ZHELEV Z， ARNAUDOVA D， TSONEV S. Genotoxicity and

erythrocyte nuclear abnormalities in Pelophylax ridibundus （Pallas，

1771） （Anura： Ranidae） in an industrial area in southern

Bulgaria： evaluation as biomarkers for ecological stress assessmen

［t J］. Acta Zoologica Bulgarica， 2022， 74（1）： 59-67.

［13］ 費(fèi)梁， 胡淑琴， 葉昌媛， 等. 中國動(dòng)物志： 兩棲綱： 下卷： 無

尾目： 蛙科［M］. 北京： 科學(xué)出版社， 2009： 1171-1175.

FEI L， HU S Q， YE C Y， et al. Fauna Sinica： Amphibia： Vol.

3： Anura： Ranidae［M］. Beijing： Science Press， 2009： 1171-

1175.

［14］ 湛天麗， 黃陽， 滕應(yīng)， 等. 貴州萬山汞礦區(qū)某農(nóng)田土壤重金

屬污染特征及來源解析［J］. 土壤通報(bào)， 2017， 48（2）： 474-

480.

ZHAN T L， HUANG Y， TENG Y， et al. Pollution characteristics

and sources of heavy metals in farmland soil in Wanshan mining

areas， Guizhou Province［J］. Chinese Journal of Soil Science，

2017， 48（2）： 474-480.

［15］ SIMON E， KUNDRáT J T， BRAUN M， et al. Impact of anthropogenic

activities on the concentration of trace elements in toe

bones of the common toad （Bufo bufo）［J］. North-Western Journal

of Zoology， 2016， 13（1）： 63-69.

［16］ 劉利， 張樂， 劉云鵬， 等. 包頭南海子濕地赤膀鴨（Mareca

strepera）組織與環(huán)境中重金屬含量的相關(guān)性［J］. 生態(tài)與農(nóng)村

環(huán)境學(xué)報(bào)， 2019， 35（4）： 476-483.

LIU L， ZHANG L， LIU Y P， et al. Correlation between heavy

metal content in Mareca strepera tissues and the environment in

the Nanhaizi Wetland， Baotou City［J］. Journal of Ecology and

Rural Environment， 2019， 35（4）： 476-483.

［17］ CHANG C Y， CHEN C Y， YIN R S， et al. Bioaccumulation of

Hg in rice leaf facilitates selenium bioaccumulation in rice

（Oryza sativa L. ） leaf in the Wanshan mercury mine［J］. Environmental

Science amp; Technology， 2020， 54（6）： 3228-3236.

［18］ 王君. 基于NMF非負(fù)矩陣模型的土壤中重金屬和多環(huán)芳烴

污染源解析［J］. 廣州化工， 2021， 49（12）： 135-138； 145.

WANG J. Non-negative matrix factorization （NMF） for source

analysis of heavy metals and polycyclic aromatic hydrocarbon

（PAHs） in soil of industrial site［J］. Guangzhou Chemical Industry，

2021， 49（12）： 135-138； 145.

［19］ 蔡新成， 花立民， 葉國輝， 等. 基于標(biāo)志重捕法的高原鼢鼠

種群動(dòng)態(tài)及其環(huán)境影響因子 ［J］. 生態(tài)學(xué)報(bào)， 2024， 44（2）：

591-601.

CAI X C， HUA L M， YE G H， et al. Population characteristics

of plateau zokor（ Eospalax baileyi） and its relationship with environmental

factors［J］. Acta Ecologica Sinica， 2024， 44（2）：

591-601.

［20］ 包俊江， 武曉東， 付和平， 等. 嚙齒動(dòng)物種群數(shù)量與年降水

量的灰色關(guān)聯(lián)度分析［J］. 草業(yè)科學(xué)， 2011， 28（4）： 642-

647.

BAO J J， WU X D， FU H P， et al. The grey correlative degree

analysis of rodent population number and annual precipitation

［J］. Pratacultural Science， 2011， 28（4）： 642-647.

［21］ 侯蒙京， 高金龍， 葛靜， 等. 青藏高原東部高寒沼澤濕地動(dòng)

態(tài)變化及其驅(qū)動(dòng)因素研究［J］. 草業(yè)學(xué)報(bào)， 2020， 29（1）： 13-

27.

HOU M J， GAO J L， GE J， et al. An analysis of dynamic

changes and their driving factors in marsh wetlands in the eastern

Qinghai-Tibet Plateau［J］. Acta Prataculturae Sinica， 2020， 29

（1）： 13-27.

［22］ DOS SANTOS LIMA G， MENEGARIO A A， SUAREZ C A， et

al. Pelagic and estuarine birds as sentinels of metal（loid）s in the

South Atlantic Ocean： ecological niches as main factors acting on

bioaccumulation［J］. Environmental Pollution， 2023， 326：

121452.

［23］ 張曉東， 盧志宏， 王婷， 等. 萬山汞礦區(qū)澤陸蛙重金屬富集特

征及健康評(píng)價(jià)［J］. 環(huán)境科學(xué)與技術(shù)， 2023， 46（10）： 170-180.

ZHANG X D， LU Z H， WANG T， et al. Heavy metal enrichment

characteristics and health assessment of Fejervarya multi?

striata in Wanshan mercury mining area， Guizhou Province［J］.

Environmental Science amp; Technology， 2023， 46（10）： 170-180.

［24］ 茍?bào)w忠， 阮運(yùn)飛. 萬山汞礦區(qū)土壤重金屬污染特征及來源解

析［J］. 化工環(huán)保， 2020， 40（3）： 336-341.

GOU T Z， RUAN Y F. Characteristics and source of heavy metals

in contaminated soil in Wanshan mercury mine area［J］. Environmental

Protection of Chemical Industry， 2020， 40（3）：

336-341.

［25］ 李秀華， 趙玲， 滕應(yīng)， 等. 貴州汞礦區(qū)周邊農(nóng)田土壤汞鎘復(fù)

合污染特征空間分布及風(fēng)險(xiǎn)評(píng)估［J］. 生態(tài)環(huán)境學(xué)報(bào)， 2022，

31（8）： 1629-1636.

LI X H， ZHAO L， TENG Y， et al. Characteristics， spatial distribution

and risk assessment of combined mercury and cadmium

pollution in farmland soils surrounding mercury mining areas in

Guizhou［J］. Ecology and Environmental Sciences， 2022， 31

（8）： 1629-1636.

［26］ 胡國成， 張麗娟， 齊劍英， 等. 貴州萬山汞礦周邊土壤重金

屬污染特征及風(fēng)險(xiǎn)評(píng)價(jià)［J］. 生態(tài)環(huán)境學(xué)報(bào)， 2015， 24（5）：

879-885.

HU G C， ZHANG L J， QI J Y， et al. Contaminant characteristics

and risk assessment of heavy metals in soils from Wanshan

mercury mine area， Guizhou Province［J］. Ecology and Environmental

Sciences， 2015， 24（5）： 879-885.

［27］ SHAFI J， MIRZA Z S， KOSOUR N， et al. Assessment of heavy

metals content in organs of edible fish species of River Ravi in

Pakistan［J］. Journal of Animal amp; Plant Sciences， 2023， 33

（1）： 201-210.

［28］ 韓瑞杰， 任逸晨， 黃濤， 等. 包頭市三類濕地中重金屬污染

程度及生物富集研究［J］. 環(huán)境工程， 2019， 37（1）： 29-34.

HAN R J， REN Y C， HUANG T， et al. Study on pollution degree

and bio-concentration of heavy metals in three types of wetlands

in Baotou， China［J］. Environmental Engineering， 2019，

37（1）： 29-34.

［29］ DING J， WANG S N， YANG W Z， et al. Tissue distribution

and association of heavy metal accumulation in a free-living resident

passerine bird tree sparrow Passer montanus［J］. Environmental

Pollution， 2023， 316： 120547.

［30］ 韓瑞杰， 求瑞娟， 任逸晨， 等. 尾礦庫濕地重金屬污染的生

物富集特征及風(fēng)險(xiǎn)評(píng)價(jià)［J］. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)， 2019， 35

（1）： 22-27.

HAN R J， QIU R J， REN Y C， et al. Bioaccumulation characteristics

and risk assessment of heavy metal pollution in tailing

wetland［J］. Journal of Ecology and Rural Environment， 2019，

35（1）： 22-27.

［31］ POGáNYOVá A， HAAS M， SOLáR J. Lead content in soil，

plants， rodents， and amphibians in the vicinity of a heating

plant’s ash waste［J］. Environmental Monitoring and Assessment，

2022， 194： 21.

［32］ ROMANOVA E B， RYABININA E S， BORYAKOV A V.

Heavy metal accumulation in the tissues and organs of Pelophy?

lax ridibundus （Pallas， 1771） and Pelophylax lessonae （Camerano，

1882） （Amphibia： Ranidae） inhabiting the waterbodies

of Nizhny Novgorod［J］. Biology Bulletin， 2021， 48（10）： 1862-

1869.

［33］ 趙賓， 李成凱， 馬驊， 等. 貴州萬山汞礦區(qū)土壤地球化學(xué)特

征［J］. 世界有色金屬， 2017（11）： 264-265.

ZHAO B， LI C K， MA H， et al. The soil geochemical characteristics

of mercury mining area in Wanshan of Guizhou Province

［J］. World Nonferrous Metals， 2017（11）： 264-265.

［34］ ZHANG H， FENG X B， LARSSEN T. Selenium speciation， distribution，

and transport in a river catchment affected by mercury

mining and smelting in Wanshan， China［J］. Applied Geochemistry，

2014， 40： 1-10.

［35］ 廖澤源， 李杰芹， 沈智杰， 等. 重慶某鐵礦區(qū)周邊耕地土壤

重金屬污染評(píng)價(jià)及來源解析［J］. 環(huán)境科學(xué)， 2024， 45（4）：

2450-2460.

LIAO Z Y， LI J Q， SHEN Z J， et al. Assessment and source

analysis of heavy metal pollution in arable land around an iron

ore mining area in Chongqing［J］. Environmental Science，

2024， 45（4）： 2450-2460.

［36］ ZHANG X W， WEI S， SUN Q Q， et al. Source identification

and spatial distribution of arsenic and heavy metals in agricultural

soil around Hunan industrial estate by positive matrix factorization

model， principle components analysis and geo statistical

analysis［J］. Ecotoxicology and Environmental Safety， 2018，

159： 354-362.

［37］ VASCONCELOS-TEIXEIRA R， TITON S C M， TITON B， et

al. Stress response， immunity， and organ mass in toads （Rhi?

nella diptycha） living in metal-contaminated areas［J］. Biological

Trace Element Research， 2022， 200（2）： 800-811.

［38］ ZHANG W Y， SUN H， SU R， et al. Fat rather than health： ecotoxic

responses of Bufo raddei to environmental heavy metal

stress during the non-breeding season［J］. Ecotoxicology and Environmental

Safety， 2022， 244： 114040.

［39］ 韓曦英， 朱健勛， 于澎. 中國林蛙體內(nèi)重金屬與環(huán)境的相關(guān)

性分析［J］. 野生動(dòng)物學(xué)報(bào)， 2018， 39（3）： 633-638.

HAN X Y， ZHU J X， YU P. Correlation analysis of heavy metals

in Rana chensinensis and their environment［J］. Chinese Journal

of Wildlife， 2018， 39（3）： 633-638.

［40］ PFLEEGER A Z， EAGLES-SMITH C A， KOWALSKI B M， et

al. From tails to toes： developing nonlethal tissue indicators of

mercury exposure in five amphibian species［J］. Ecotoxicology，

2016， 25（3）： 574-583.