劉浪,吳恒,許先鵬,方向陽,王玥琳

摘要：為探究滇中高原哀牢山不同林分類型土壤重金屬的分布特征及對(duì)其潛在生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)，選取哀牢山4種林分（常綠闊葉林、高山櫟林、華山松林、滇油杉林）土壤為研究對(duì)象，采集林分土壤，對(duì)該土壤進(jìn)行7種重金屬（Zn、Cu、Mn、Cr、Pb、As、Ni）、有機(jī)質(zhì)含量和pH進(jìn)行檢測(cè)分析。以云南省土壤背景值和國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（Ⅰ級(jí)）作為參比值，采用Hakanson 潛在生態(tài)危害指數(shù)法（Potential ecologicalrisk index）評(píng)價(jià)研究區(qū)不同林分土壤重金屬潛在生態(tài)危害效應(yīng)。結(jié)果表明，研究區(qū)內(nèi)土壤中Zn的含量最高，其次是Mn，其余5種都在5 mg/kg以下，pH均為酸性，有機(jī)質(zhì)的含量常綠闊葉林最高，滇油杉林、華山松林次之，高山櫟林最低，變異系數(shù)為18.32% ～ 37.26%。常綠闊葉林土壤中Pb和pH為顯著負(fù)相關(guān)（P＜0.05），與有機(jī)質(zhì)顯著正相關(guān)（P＜0.05）；高山櫟林Pb和Ni極顯著負(fù)相關(guān)（P＜0.01），As和有機(jī)質(zhì)顯著正相關(guān)（P＜0.05）；華山松林Mn與Pb顯著正相關(guān)（P＜0.05）；滇油杉林Zn和Pb、Ni、有機(jī)質(zhì)為負(fù)相關(guān)外，其余4種重金屬和pH均為正相關(guān)，且相關(guān)性均不顯著（P＞0.05）；有機(jī)質(zhì)與Zn的相關(guān)性較強(qiáng)可聚類為一組。4種林分土壤重金屬的單種和多種土壤重金屬潛在生態(tài)危害都處于輕微風(fēng)險(xiǎn)狀態(tài)，單種重金屬潛在生態(tài)危害程度平均值由大到小排序?yàn)锳s、Zn、Ni、Pb、Cu、Cr、Mn，多種土壤重金屬潛在生態(tài)危害程度由大到小排序?yàn)槿A山松林、滇油杉林、常綠闊葉林、高山櫟林。

關(guān)鍵詞：滇中高原；哀牢山；森林土壤；重金屬含量；生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

中圖分類號(hào)：S714;X53;X826文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1006-8023（2024）03-0056-10

Heavy Metal Pollution Characteristics and Ecological Risk Assessment?of Soil in Different Forest Stands of Ailao Mountain in the?Central Yunnan Plateau

LIU Lang1， WU Heng1， XU Xianpeng1， FANG Xiangyang2， WANG Yuelin3*

（1.Southwest Survey and Planning Institute of National Forestry and Grassland Administration， Kunming 650031， China;?2.Kunming Forest Fire Prevention and Control and Forest Grass Information Center， Kunming 650506， China;?3.College of Tourism and Urban-Rural Planning， Chengdu University of Technology， Chengdu 610059， China）

Abstract：In order to investigate the distribution characteristics of heavy metals in soil of different forest types in Ailao Mountain on the central Yunnan Plateau and evaluate their potential ecological risks， four forest types （evergreen broad-leaved forest， alpine oak forest， Pinus armandii forest， and Yunnan oil fir forest） were selected as the research object. The forest soil was collected to detect and analyze the content of seven heavy metals （Zn， Cu， Mn， Cr， Pb， As， Ni）， organic matter， and pH of the soil. Using the soil background value of Yunnan Province and the national soil environmental quality standard （Level I） as reference values， the Hakanson potential ecological hazard index method was used to evaluate the potential ecological hazard effects of heavy metals in soil of different forest stands in the study area. The results showed that the content of Zn in the soil of the study area was the highest， followed by Mn， and the other five species were all below 5 mg/kg. The pH was acidic， and the content of organic matter was highest in evergreen broad-leaved forest， followed by Yunnan oil fir forest and Pinus armandii forest， and the lowest in alpine oak forest， with a coefficient of variation?between 18.32% and 37.26%. There was a significant negative correlation between Pb and pH in the soil of evergreen broad-leaved forest （P<0.05）， and a significant positive correlation with organic matter （P<0.05）. Pb and Ni were significantly negatively correlated （P<0.01）， while As and organic matter were significantly positively correlated （P<0.05） in the alpine oak fores. There was a significant positive correlation between Mn and Pb in Pinus armandii forest （P<0.05）. Except for the negative correlation between Zn and Pb， Ni， and organic matter in the Yunnan oil fir forest， the other four heavy metals and pH values were all positive correlated， and the correlation was not significant （P>0.05）. The correlation between organic matter and Zn was strong and can be clustered into a group. The potential ecological hazards of single and multiple heavy metals in the soil of four forest stands were all in a mild risk state. The average degree of potential ecological hazards of single heavy metals was ranked from high to low as As， Zn， Ni， Pb， Cu， Cr， and Mn， while the potential ecological hazards of multiple soil heavy metals were ranked from high to low as Pinus armandii forest， Yunnan oil fir forest， evergreen broad-leaved forest and alpine oak forest.

Keywords：Central Yunnan Plateau; Ailao Mountain; forest soil; heavy metal content; ecological risk assessment

0引言

森林生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的重要組成部分，為人類的生產(chǎn)、生活提供賴以生存的條件，也是陸地上最大的碳源碳匯，在氣候調(diào)節(jié)、水土保持、水源涵養(yǎng)和空氣凈化等方面起著非常重要的作用[1-2]。森林下方土壤質(zhì)量的好壞關(guān)乎森林生態(tài)系統(tǒng)的健康持久發(fā)展，其生物及非生物因素都是影響植物群落生存的重要環(huán)境因子，是森林生態(tài)系統(tǒng)不可缺少的重要部分[3-4]。目前隨著工業(yè)的迅速發(fā)展，大量重金屬通過大氣和水體等途徑不斷排放到環(huán)境中，通過降雨、空氣沉降及植物攔截等作用最終不斷輸入到森林土壤中，并且相關(guān)研究表明土壤雖然可容納環(huán)境系統(tǒng)中70%以上的重金屬[5]，但土壤中的重金屬具有長時(shí)間殘留、強(qiáng)隱蔽性和毒性大等特點(diǎn)[6]。森林土壤重金屬含量達(dá)到一定閾值將會(huì)對(duì)生態(tài)系統(tǒng)中的其他生物造成一定威脅，不利于其他生物生存，減少生物多樣性，從而影響整個(gè)生態(tài)系統(tǒng)的穩(wěn)定性，并通過林產(chǎn)品以及地下水等形式危害人類健康[7]。有研究發(fā)現(xiàn)，不同植物對(duì)重金屬吸附和累積的能力有所不同，林地土壤重金屬含量也有所不同，森林群落對(duì)重金屬具有明顯的消減作用[8-9]，一方面森林植被生物量較大，生命周期長，大量的重金屬可長期存儲(chǔ)在體內(nèi)，另一方面植物發(fā)達(dá)的根系能夠深入土層吸附更深層土壤中的重金屬[10]。

哀牢山國家級(jí)自然保護(hù)區(qū)境內(nèi)分布著滇中最具代表性的中山濕性常綠闊葉林以及高山櫟林、滇油杉林等多種典型林分[11]，而目前研究主要集中在森林植被的土壤營養(yǎng)元素[12]、土壤水分動(dòng)態(tài)[13]、枯落物[14]和碳儲(chǔ)量[15]等方面，對(duì)于土壤重金屬污染方面的研究還未見報(bào)道。而該研究區(qū)對(duì)于本地的水源涵養(yǎng)、生物多樣性等方面起著重要的生態(tài)功能[16]，因此很有必要對(duì)該地區(qū)不同林分類型土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)分析評(píng)價(jià)，掌握研究區(qū)森林土壤重金屬分布特征、污染程度和潛在風(fēng)險(xiǎn)，為森林生態(tài)安全提供科學(xué)保障。

1研究區(qū)概況

哀牢山國家級(jí)自然保護(hù)區(qū)位于云南省中部，地處云貴高原、橫斷山地和青藏高原三大自然地理區(qū)域的接合部。地理位置為23.35°～24.44° E，100.54°～101.30° N，海拔1 300～2 750 m，受西南季風(fēng)的影響，干濕季分明，年平均氣溫11.3 ℃，年平均降水量為1 947 mm，年均蒸發(fā)量1 192 mm，全年日照時(shí)間2 380 h。研究區(qū)內(nèi)由于海拔較高，氣候垂直分布明顯，形成的植被帶分布明顯，從山體下部到上部大致可以劃分為3個(gè)植被帶，分別為山體中下部的亞熱帶季風(fēng)常綠闊葉林、山體中部的中山濕性常綠闊葉林和山體頂部的苔蘚矮林。獨(dú)特的氣候特點(diǎn)蘊(yùn)育了我國目前面積最大以云南特有植物物種為優(yōu)勢(shì)的亞熱帶常綠闊葉林生長在哀牢山，使其保存著處于原始狀態(tài)、完整而穩(wěn)定的亞熱帶山地森林生態(tài)系統(tǒng)。

2研究方法

2.1樣品的采集及處理

利用GPS定位根據(jù)哀牢山森林的不同海拔，分別選取常綠闊葉林、高山櫟林、華山松林和滇油杉林4個(gè)植物群落，見表1，每個(gè)群落中按照10 m×10 m的網(wǎng)格面積布設(shè)4個(gè)典型樣方，共布設(shè)16個(gè)典型樣方。每個(gè)樣方按照上、中、下坡位各設(shè)置采樣點(diǎn)1個(gè)，共布設(shè)采樣點(diǎn)48個(gè)。每個(gè)采樣點(diǎn)分別按0～20、20～40、40～60 cm，3個(gè)層次采集各土壤剖面樣品。采樣時(shí)按照由上到下的次序分層采樣，每層采集500 g左右，隨后將同一樣地同一坡度同一層土壤均勻混合，采用四分法收集500 g左右，剔除土壤樣品中的石頭和植物的根、葉等雜物，放入標(biāo)好號(hào)的聚乙烯塑料袋密封，帶回試驗(yàn)室通風(fēng)處進(jìn)行自然風(fēng)干，風(fēng)干后磨碎過孔徑為100目的尼龍網(wǎng)篩，放入密封袋待測(cè)。

2.2分析方法

2.2.1樣品化學(xué)分析方法

土壤重金屬Zn、Cu、Mn、Cr、Pb、As、Ni采用酸溶法后用美譜達(dá)V-1800原子吸收分光光度計(jì)測(cè)定，pH用pHS-3C酸度計(jì)測(cè)定，土壤有機(jī)質(zhì)采用重鉻酸鉀氧化外加熱法，每份樣品設(shè)置3個(gè)重復(fù)，最終取3次試驗(yàn)分析的平均值為該樣點(diǎn)的最終分析值。

2.2.2土壤重金屬潛在生態(tài)危害評(píng)價(jià)方法

本研究選用Hakanson 提出的潛在生態(tài)危害指數(shù)法（Potential ecologicalrisk index）[17]，以云南省土壤背景值[18]和國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（Ⅰ級(jí)）[19]作為參比值，見表2，確定每種林分類型土壤重金屬潛在生態(tài)危害效應(yīng)及空間分布特性，評(píng)價(jià)玉溪哀牢山常綠闊葉林、高山櫟林、華山松林和滇油杉林4個(gè)林分的土壤重金屬潛在生態(tài)危害效應(yīng)。

單種重金屬潛在生態(tài)危害系數(shù)（Ei）的計(jì)算公式為

Ei=T ir×CiC0i。（1）

式中：T ir為重金屬的毒性系數(shù)；Ci為第 i 種重金屬的實(shí)測(cè)含量，mg/kg；Coi為參比值，mg/kg，即土壤背景值。

本研究以云南省土壤背景值和國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（Ⅰ級(jí)）作為參比值。毒性系數(shù)根據(jù)徐爭(zhēng)啟等[20]計(jì)算的重金屬毒性系數(shù)（T ir）：Zn=Mn=1、Cr=2、Cu=Ni=Pb=5、As=10。

多種重金屬潛在生態(tài)危害指數(shù)（RI）的計(jì)算公式為

RI=∑mi=1Ei。（2）

根據(jù)Ei和RI ，參照沉積物（土壤）中重金屬潛在生態(tài)危害系數(shù)、生態(tài)危害指數(shù)和污染程度的關(guān)系[20]， 采用丁振華等[21]給出7種重金屬的潛在生態(tài)危害指數(shù)RI 劃分的潛在生態(tài)危害分級(jí)標(biāo)準(zhǔn)，見表3。

2.2.3數(shù)據(jù)處理方法

本研究采用Excel2013、SPSS20.0統(tǒng)計(jì)軟件進(jìn)行數(shù)據(jù)分析，聚類分析在R 4.0.5中利用“pheatmap”包[22]完成分析。

3結(jié)果與分析

3.1不同林分類型森林土壤重金屬含量及其分布特征

哀牢山4種林分類型下土壤的pH、有機(jī)質(zhì)與重金屬含量見表4。7種土壤重金屬中，Zn的含量最高，其次是Mn，其余5種都在5 mg/kg以下，Pb的含量最低。Mn在常綠闊葉林中變幅最大，且變異系數(shù)最大，常綠闊葉林中最大值是最小值的6倍多，且4種林分其平均含量由大到小為：常綠闊葉林、滇油杉林、華山松林、高山櫟林，說明Mn在常綠闊葉林土壤中含量分布最不均勻。每種林分土壤中Pb含量變幅最大，其平均含量滇油杉林最高，高山櫟林次之，常綠闊葉林和華山松林相同，且在后2種林分土壤中分布較不均勻，變異系數(shù)都在60%以上。As在4種林分土壤中變幅最小，華山松林土壤As的平均含量相對(duì)最高，其余3種差別不大?？偟膩砜矗嵊蜕剂滞寥乐?種重金屬含量變異系數(shù)最大，其分布最不均勻，其次是常綠闊葉林、華山松林，高山櫟林最小。4種林分類型土壤的pH均為酸性，差異相對(duì)較小，pH的變異系數(shù)值均在8.11%以下。有機(jī)質(zhì)的含量不同林分差異較大，其平均含量常綠闊葉林最高，滇油杉林、華山松林次之，高山櫟林最低，最大值為最小值的3倍多，且各林分類型下變異系數(shù)為18.32%～37.26%。

3.2土壤重金屬含量與pH、有機(jī)質(zhì)的相關(guān)性分析

對(duì)4種林分類型下土壤中重金屬、有機(jī)質(zhì)含量及pH進(jìn)行相關(guān)性分析，結(jié)果見表5。常綠闊葉林土壤中Zn和Cu、Ni、pH呈正相關(guān)性外，與其余4種重金屬及有機(jī)質(zhì)呈負(fù)相關(guān)，且相關(guān)性不顯著（P＞0.05）；Pb和pH表現(xiàn)出顯著負(fù)相關(guān)性（P＜0.05），與有機(jī)質(zhì)表現(xiàn)出正相關(guān)性顯著（P＜0.05）；其余重金屬之間、重金屬與有機(jī)質(zhì)和pH之間相關(guān)性均不顯著（P＞0.05）。高山櫟林土壤中Zn和Cu、Mn、Cr、Pb呈正相關(guān)性外，與其余2種重金屬及pH、有機(jī)質(zhì)呈負(fù)相關(guān)，且相關(guān)性不顯著（P＞0.05）；Pb和Ni表現(xiàn)出極顯著負(fù)相關(guān)性（P＜0.01），相關(guān)性系數(shù)為1；As和有機(jī)質(zhì)表現(xiàn)出顯著正相關(guān)性（P＜0.05），其余重金屬之間、重金屬與有機(jī)質(zhì)和pH之間相關(guān)性均不顯著（P＞0.05）。華山松林土壤中Zn和As、有機(jī)質(zhì)為正相關(guān)外，其余5種重金屬和pH均為負(fù)相關(guān)，且相關(guān)性均不顯著（P＞0.05）；Mn與Pb表現(xiàn)出顯著正相關(guān)性（P＜0.05），有機(jī)質(zhì)和pH表現(xiàn)出顯著負(fù)相關(guān)性（P＜0.05）；其余重金屬之間、重金屬與有機(jī)質(zhì)和pH之間相關(guān)性均不顯著（P＞0.05）。滇油杉林土壤中Zn和Pb、Ni、有機(jī)質(zhì)為正相關(guān)外，其余4種重金屬和pH均為正相關(guān)，且相關(guān)性均不顯著（P＞0.05）?？傮w來看，pH與重金屬之間的相關(guān)性不顯著，也無明顯規(guī)律可循；而有機(jī)質(zhì)與不同林分類型土壤中的重金屬含量均有一定的正/負(fù)相關(guān)性，且相關(guān)系數(shù)值較大，說明有機(jī)質(zhì)對(duì)重金屬可能存在一定的吸附效果，這與Brookes[23]的研究結(jié)果相似。

3.3不同林分土壤重金屬、有機(jī)質(zhì)含量及pH聚類分析

對(duì)不同林分類型土壤重金屬、有機(jī)質(zhì)含量及pH進(jìn)行層次聚類分析，得到相關(guān)性聚類熱圖（圖1）。由圖1（圖上部）可看出，最底層滇油杉林和華山松林的重金屬含量更相似可以劃分為第1組，其次是常綠闊葉林與第1組劃分成第2層次為第2組，最后是高山櫟林與第2組劃分成第3層次為第3組，可能與林分的地理空間位置相近，土壤本底值相似有關(guān)系。從重金屬、有機(jī)質(zhì)、pH看（圖1左側(cè)），除了Zn之外的6種重金屬與pH的相關(guān)性較強(qiáng)可聚類為一組，有機(jī)質(zhì)（OM）與Zn的相關(guān)性較強(qiáng)可聚類為一組，可能二者之間的相互作用較其他重金屬更為明顯。

3.4不同林分土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

由表6看出，無論用云南省土壤背景值還是國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（Ⅰ級(jí)）作為參比值，4種林分類型下每種土壤重金屬的單種重金屬潛在生態(tài)危害系數(shù)（Ei）均遠(yuǎn)遠(yuǎn)小于40，都處于輕微風(fēng)險(xiǎn)狀態(tài)，每種林分下多種土壤重金屬潛在生態(tài)危害指數(shù)（RI）也均遠(yuǎn)小于150，處于輕微風(fēng)險(xiǎn)狀態(tài)。4種林分土壤單種重金屬潛在生態(tài)危害系數(shù)（Ei）危害程度平均值由大到小排序?yàn)椋篈s、Zn、Ni、Pb、Cu、Cr、Mn，可能與As的毒性系數(shù)較大有關(guān)系，使As的潛在生態(tài)風(fēng)險(xiǎn)系數(shù)計(jì)算結(jié)果偏大。從林分類型來看，多種土壤重金屬潛在生態(tài)危害指數(shù)（RI）危害程度由大到小排序?yàn)椋喝A山松林、滇油杉林、常綠闊葉林、高山櫟林。

4結(jié)論與討論

哀牢山不同林分類型下7種土壤重金屬中，Zn的含量最高，其次是Mn，其余5種都在5 mg/kg以下，Pb的含量最低，且每種重金屬在不同林分含量變異系數(shù)各不相同，在不同采樣點(diǎn)分布較不均勻，表明研究區(qū)內(nèi)土壤重金屬空間分布差異特征比較明顯，這與其他學(xué)者的研究結(jié)果一致[24-25]。并且森林土壤重金屬的含量、分布特征，受重金屬元素的來源、種類以及林地土壤的理化性質(zhì)和原有的土壤母質(zhì)影響很大[26-27]。簡單來說，土壤中重金屬既可能本身就存在于原有土壤母質(zhì)中，也有可能是由人類的生產(chǎn)生活等活動(dòng)帶入[28]?？梢钥隙ǖ氖牵煌N類的重金屬但其來源相同，說明這幾類重金屬可能存在一定的相關(guān)性[25]。分析結(jié)果顯示，不同林分類型土壤重金屬之間表現(xiàn)出的相關(guān)性各異，比如高山櫟林中Pb和Ni表現(xiàn)出極顯著負(fù)相關(guān)性（P＜0.01），相關(guān)性系數(shù)為1；As和有機(jī)質(zhì)表現(xiàn)出顯著正相關(guān)性（P＜0.05），華山松林中Mn與Pb表現(xiàn)出顯著正相關(guān)性（P＜0.05），其他林分類型下各種土壤重金屬之間相關(guān)性不顯著（P＞0.05），說明這幾種土壤重金屬的來源存在相同的伴生性[29]。

從土壤的理化性質(zhì)（有機(jī)質(zhì)含量及pH）與重金屬的關(guān)系來分析，其理化性質(zhì)對(duì)重金屬起到累積和吸附的效果[30]，土壤重金屬含量與有機(jī)質(zhì)含量呈明顯的相關(guān)性[31]，并且pH越高土壤對(duì)重金屬離子的吸附效果越明顯[32]。本研究結(jié)果顯示，常綠闊葉林中Pb與有機(jī)質(zhì)表現(xiàn)出顯著正相關(guān)性（P＜0.05），高山櫟林中As和有機(jī)質(zhì)表現(xiàn)出顯著正相關(guān)性（P＜0.05）且相關(guān)系數(shù)值較大，說明土壤中的有機(jī)質(zhì)可能對(duì)重金屬產(chǎn)生一定的吸附作用，這一研究與Brookes[23]的結(jié)果相似。常綠闊葉林中Pb和pH表現(xiàn)出顯著負(fù)相關(guān)性（P＜0.05），而其他重金屬含量與林分中pH間的相關(guān)性關(guān)系均不顯著（P＞0.05）。同時(shí)從土壤重金屬和理化性質(zhì)的層次聚類分析結(jié)果也證實(shí)，有機(jī)質(zhì)與Zn的相關(guān)性較強(qiáng)可聚類為一組，這些結(jié)果都進(jìn)一步說明了土壤理化性質(zhì)與重金屬之間關(guān)系的復(fù)雜性[24]。

以云南省土壤背景值和國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（Ⅰ級(jí)）作為參比值的單種和多種土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果均為輕微污染程度，可能與研究區(qū)屬于天然林并遠(yuǎn)離城區(qū)處于比較原始的狀態(tài)，土壤重金屬主要屬于自然來源，一是來源于巖石圈內(nèi)地球化學(xué)過程中形成的原生的土壤母質(zhì)，二是植被攔截到的重金屬塵降，三是天然土壤污染[33]，幾乎不受人類工農(nóng)業(yè)活動(dòng)的影響。4種林分土壤單種重金屬潛在生態(tài)危害系數(shù)（Ei）危害程度平均值由大到小排序?yàn)椋篈s、Zn、Ni、Pb、Cu、Cr、Mn，這樣的分析結(jié)果顯示利用重金屬潛在生態(tài)風(fēng)險(xiǎn)法，在計(jì)算的過程中加入了毒性系數(shù)（T i r），這樣的計(jì)算方式充分考慮了不同重金屬的生物毒性，從而使評(píng)價(jià)結(jié)果更綜合、全面地評(píng)價(jià)[34]。評(píng)價(jià)結(jié)果顯示As的毒性系數(shù)比其他種類的重金屬都大，與As的毒性系數(shù)最大值（10）有關(guān)系[20]，使得As的潛在生態(tài)風(fēng)險(xiǎn)系數(shù)計(jì)算結(jié)果最大。并且研究結(jié)果顯示，從林分類型來看，多種土壤重金屬潛在生態(tài)危害指數(shù)（RI）危害程度由大到小排序?yàn)椋喝A山松林、滇油杉林、常綠闊葉林、高山櫟林?？赡芘c不同林分的植被凈化作用不同有關(guān)，并且華山松林生長的坡度相較于其他林分要平穩(wěn)，對(duì)于大氣沉降帶來的重金屬更有利于存儲(chǔ)，而其他林分所處坡度較陡，降雨對(duì)土壤層起到更大強(qiáng)度的淋溶和沖刷作用[24]，土壤中重金屬含量相對(duì)較低，因此在研究區(qū)內(nèi)華山松林對(duì)土壤重金屬的吸附和存儲(chǔ)效果較其他林分明顯，直接影響到重金屬在土壤中的釋放、遷移和沉積等自然循環(huán)過程，以及重金屬在自然界的毒性強(qiáng)弱[6，35]，最終導(dǎo)致華山松林多種土壤重金屬潛在生態(tài)危害指數(shù)最大。但總體來看，由于4種天然林中重金屬總量較少或由于植被的凈化作用，研究區(qū)森林土壤重金屬污染屬于輕微污染狀態(tài)。

【參考文獻(xiàn)】

[1]張禎祺，蔡惠文，張平平，等.基于GEE遙感云平臺(tái)的三江源植被碳源/匯時(shí)空變化研究[J].自然資源遙感，2023，35（1）：231-242.

ZHANG Z Q， CAI H W， ZHANG P P， et al. A GEE-based study on the temporal and spatial variations in the carbon source/sink function of vegetation in the Three-River Headwaters region[J]. Remote Sensing for Natural Resources， 2023， 35（1）： 231-242.

[2]李亦秋，鄧歐，李若霜，等.黔中喀斯特城市新區(qū)森林生態(tài)系統(tǒng)水土保持及其生態(tài)效益[J].水土保持通報(bào)，2022，42（4）：244-250.

LI Y Q， DENG O， LI R S， et al. Soil and water conservation and ecological benefits of forest ecosystem in new urban Karst areas of central Guizhou Province[J]. Bulletin of Soil and Water Conservation， 2022， 42（4）： 244-250.

[3]ILANY T， ASHTON M S， MONTAGNINI F， et al. Using agroforestry to improve soil fertility： effects of intercropping on Ilex paraguariensis （yerba mate） plantations with Araucaria angustifolia[J]. Agroforestry Systems， 2010， 80（3）： 399-409.

[4]單奇華，俞元春，張建鋒，等.城市森林土壤肥力質(zhì)量指標(biāo)篩選：以南京市為例[J].土壤，2009，41（5）：777-783.

SHAN Q H， YU Y C， ZHANG J F， et al. Determination of fertility index of urban forest soil in Nanjing city[J]. Soils， 2009， 41（5）： 777-783.

[5]NRIAGU J O. A history of global metal pollution[J]. Science， 1996， 272（5259）： 223.

[6]韓春梅，王林山，鞏宗強(qiáng)，等.土壤中重金屬形態(tài)分析及其環(huán)境學(xué)意義[J].生態(tài)學(xué)雜志，2005，24（12）：1499-1502.

HAN C M， WANG L S， GONG Z Q， et al. Chemical forms of soil heavy metals and their environmental significance[J]. Chinese Journal of Ecology， 2005， 24（12）： 1499-1502.

[7]郭笑笑，劉叢強(qiáng)，朱兆洲，等.土壤重金屬污染評(píng)價(jià)方法[J].生態(tài)學(xué)雜志，2011，30（5）：889-896.

GUO X X， LIU C Q， ZHU Z Z， et al. Evaluation methods for soil heavy metals contamination： a review[J]. Chinese Journal of Ecology， 2011， 30（5）： 889-896.

[8]盧德亮，喬璐，陳立新，等.哈爾濱市區(qū)綠地土壤重金屬污染特征及植物富集[J].林業(yè)科學(xué)，2012，48（8）：16-24.

LU D L， QIAO L， CHEN L X， et al. Soil pollution characteristics by heavy metals and the plant enrichment in green space of urban areas of Harbin[J]. Scientia Silvae Sinicae， 2012， 48（8）： 16-24.

[9]劉曉輝，呂憲國，劉惠清.溝谷地不同植被下土壤重金屬縱向分異研究[J].環(huán)境科學(xué)，2007，28（12）：2766-2770.

LIU X H， LYU X G， LIU H Q. Longitudinal variation of heavy metals contents in topsoil covered with different vegetation in gully area[J]. Environmental Science， 2007， 28（12）： 2766-2770.

[10]WANG G， SU M Y， CHEN Y H， et al. Transfer characteristics of cadmium and lead from soil to the edible parts of six vegetable species in southeastern China[J]. Environmental Pollution， 2006， 144（1）： 127-135.

[11]牛躍，劉洋，張一平，等.哀牢山植物多樣性垂直分布格局[J].西部林業(yè)科學(xué)，2008，37（4）：27-34.

NIU Y， LIU Y， ZHANG Y P， et al. Vertical distribution patterns of plant species diversity in Ailao Mountains[J]. Journal of West China Forestry Science， 2008， 37（4）： 27-34.

[12]羅鈺穎，朱玉婷，唐青，等.哀牢山不同海拔表層土壤理化性質(zhì)變化規(guī)律研究[J].西部林業(yè)科學(xué)，2020，49（1）：140-148.

LUO Y Y， ZHU Y T， TANG Q， et al. Variation pattern of physicochemical properties in different altitudes of Ailao Mountain[J]. Journal of West China Forestry Science， 2020， 49（1）： 140-148.

[13]鞏合德，張一平，劉玉洪，等.哀牢山常綠闊葉林土壤水分動(dòng)態(tài)變化[J].東北林業(yè)大學(xué)學(xué)報(bào)，2008，36（1）：53-54.

GONG H D， ZHANG Y P， LIU Y H， et al. Dynamics of soil moisture in evergreen broad-leaved forest in Ailao Mountain[J]. Journal of Northeast Forestry University， 2008， 36（1）： 53-54.

[14]馮文婷.哀牢山常綠闊葉林土壤簡單糖類和氨基酸研究：地上凋落物、地下碳輸入和植物養(yǎng)分吸收對(duì)其的影響[D].西雙版納：中國科學(xué)院研究生院（西雙版納熱帶植物園），2007.

FENG W T. Soil simple sugars and amino acids in an evergreen forest of Ailao Mountains[D]. Xishuangbanna： Xishuangbanna Tropical Botanical Garden， Chinese Academy of Sciences， 2007.

[15]張鵬超，張一平，楊國平，等.哀牢山亞熱帶常綠闊葉林喬木碳儲(chǔ)量及固碳增量[J].生態(tài)學(xué)雜志，2010，29（6）：1047-1053.

ZHANG P C， ZHANG Y P， YANG G P， et al. Carbon storage and sequestration of tree layer in subtropical evergreen broadleaf forests in Ailao Mountain of Yunnan[J]. Chinese Journal of Ecology， 2010， 29（6）： 1047-1053.

[16]李晉，陶晶，晁增華，等.哀牢山自然保護(hù)區(qū)資源保護(hù)價(jià)值評(píng)估[J].內(nèi)蒙古林業(yè)調(diào)查設(shè)計(jì)，2021，44（5）：44-47， 57.

LI J， TAO J， CHAO Z H， et al. Evaluation of resource conservation value in Ailao Mountain Nature Reserve[J]. Inner Mongolia Forestry Investigation and Design， 2021， 44（5）： 44-47， 57.

[17]XU Z Q， NI S J， TUO X G， et al. Calculation of heavymetal's toxicity coefficient in the evaluation of potential ecological risk index[J]. Environmental Science and Technology， 2008， 31（2）：112-115.

[18]國家環(huán)境保護(hù)局主持，中國環(huán)境監(jiān)測(cè)總站主編.中國土壤元素背景值[M].北京：中國環(huán)境科學(xué)出版社，1990.

Chaired by the State Environmental Protection Bureau of the People's Republic of China， Chief Editor of China Environmental Monitoring Station. Background values of soil elements in China[M]. Beijing： China Environmental Science Press， 1990.

[19]中華人民共和國生態(tài)環(huán)境部，國家市場(chǎng)監(jiān)督管理總局. 土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)（試行）： GB 15618—2018[S].北京：中國環(huán)境出版社集團(tuán)， 2018.

Ministry of Ecology and Environment of the People's Republic of China， State Administration for Market Regulation. Soil environmental quality- Risk control standard for soil contamination of agricultural land： GB 15618—2018 [S]. Beijing： China Environmental Press Group， 2018.

[20]徐爭(zhēng)啟，倪師軍，庹先國，等.潛在生態(tài)危害指數(shù)法評(píng)價(jià)中重金屬毒性系數(shù)計(jì)算[J].環(huán)境科學(xué)與技術(shù)，2008，31（2）：112-115.

XU Z Q， NI S J， TUO X G， et al. Calculation of heavy metals' toxicity coefficient in the evaluation of potential ecological risk index[J]. Environmental Science & Technology， 2008， 31（2）： 112-115.

[21]丁振華，賈洪武，劉彩娥，等.黃浦江沉積物重金屬的污染及評(píng)價(jià)[J].環(huán)境科學(xué)與技術(shù)，2006，29（2）：64-66，119.

DING Z H， JIA H W， LIU C E， et al. Pollution and assessment of heavy metals in Huangpu River sediments[J]. Environmental Science & Technology， 2006， 29（2）： 64-66， 119.

[22]KOLDE R. pheatmap： Pretty heatmaps. R package version 1.0.12[EB/OL]. https：//CRAN. R-project. org/package=pheatmap. 2019.

[23]BROOKES P C. The use of microbial parameters in monitoring soil pollution by heavy metals[J]. Biology and Fertility of Soils， 1995， 19（4）： 269-279.

[24]莊玉婷，馮嘉儀，儲(chǔ)雙雙，等.粵西地區(qū)不同林分類型土壤重金屬含量及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2018，39（5）：25-31.

ZHUANG Y T， FENG J Y， CHU S S， et al. Contents and ecology risk assessments of heavy metals in different forest soils in West Guangdong[J]. Journal of South China Agricultural University， 2018， 39（5）： 25-31.

[25]簡毅，張健，楊萬勤，等.岷江下游小型集水區(qū)不同林分土壤重金屬污染特征及其生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].生態(tài)環(huán)境學(xué)報(bào)，2015，24（9）：1526-1533.

JIAN Y， ZHANG J， YANG W Q， et al. Heavy metal pollution characteristics and ecological risk assessment of soils in different forest stands of a small catchment at the downstream of Minjiang River[J]. Ecology and Environmental Sciences， 2015， 24（9）： 1526-1533.

[26]MANTA D S， ANGELONE M， BELLANCA A， et al. Heavy metals in urban soils： a case study from the city of Palermo （Sicily）， Italy[J]. Science of the Total Environment， 2002， 300（1/2/3）： 229-243.

[27]楚春暉，佘宇晨，佘濟(jì)云，等.亞熱帶不同森林類型的土壤重金屬空間分布特征及其潛在生態(tài)風(fēng)險(xiǎn)[J].水土保持學(xué)報(bào)，2014，28（5）：258-263.

CHU C H， SHE Y C， SHE J Y， et al. Spatial distribution and ecological risk assessment of heavy metals in soil from different forests in subtropics[J]. Journal of Soil and Water Conservation， 2014， 28（5）： 258-263.

[28]潘勇軍，陳步峰，肖以華，等.廣州市城市森林土壤重金屬污染狀況及其評(píng)價(jià)[J].生態(tài)環(huán)境，2008，17（1）：210-215.

PAN Y J， CHEN B F， XIAO Y H， et al. Heavy metal pollution status and evaluation of urban forest soils in Guangzhou[J]. Ecology and Environment， 2008， 17（1）： 210-215.

[29]雷康平，黃鋮，朱榮莉，等.涼山州土壤重金屬與粳稻中重金屬相關(guān)性評(píng)價(jià)[J].農(nóng)業(yè)研究與應(yīng)用，2022，35（2）：62-69.

LEI K P， HUANG C， ZHU R L， et al. Correlation evaluation of heavy metals in soil and Japonica rice in Liangshan prefecture[J]. Agricultural Research and Application， 2022， 35（2）： 62-69.

[30]VANDECASTEELE B， QUATAERT P， TACK F M G. Uptake of Cd， Zn and Mn by willow increases during terrestrialisation of initially ponded polluted sediments[J]. Science of the Total Environment， 2007， 380（1/2/3）： 133-143.

[31]許超，夏北成，何石媚，等.大寶山礦山下游地區(qū)稻田土壤重金屬含量特征[J].中山大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，47（3）：122-127.

XU C， XIA B C， HE S M， et al. Characteristics of Cd， Zn， Pb， Cu content of paddy soils in the lower stream of Dabaoshan area， Guangdong[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni， 2008， 47（3）： 122-127.

[32]喬冬梅，齊學(xué)斌，龐鴻賓，等.不同pH值對(duì)重金屬Pb2+形態(tài)的影響研究[J].水土保持學(xué)報(bào)，2010，24（6）：173-176.

QIAO D M， QI X B， PANG H B， et al. The effect of pH level on existing form of Pb2+[J]. Journal of Soil and Water Conservation， 2010， 24（6）： 173-176.

[33]徐鴻志，常江.安徽省主要土壤重金屬污染評(píng)價(jià)及其評(píng)價(jià)方法研究[J].土壤通報(bào)，2008，39（2）：411-415.

XU H Z， CHANG J. Evaluation and methods of heavy metal pollution in main soils of Anhui Province[J]. Chinese Journal of Soil Science， 2008， 39（2）： 411-415.

[34]徐燕，李淑芹，郭書海，等.土壤重金屬污染評(píng)價(jià)方法的比較[J].安徽農(nóng)業(yè)科學(xué)，2008，36（11）：4615-4617.

XU Y， LI S Q， GUO S H， et al. Comparison of assessment methods of heavy metal pollution in soil[J]. Journal of Anhui Agricultural Sciences， 2008， 36（11）： 4615-4617.

[35]錢進(jìn)，王子健，單孝全.土壤中微量金屬元素的植物可給性研究進(jìn)展[J].環(huán)境科學(xué)，1995，16（6）：73-75.

QIAN J， WANG Z J， SHAN X Q. Progress in the investigation on plant availability of soil trace metals[J]. Chinese Journal of Enviromental Science， 1995， 16（6）： 73-75.