-
镉是一种有毒的金属元素。据《全国土壤污染状况调查公报》,镉的点位超标率达7.0%,是点位超标率最高的污染物之一。土壤中的镉可通过食物链进入人体,进而对人的身体健康造成严重的损害[1],因此亟待对镉污染土壤进行修复。
植物修复是指利用植物吸收土壤中的重金属,最终清除土壤中重金属的一种环境修复技术[2],具有经济、绿色、效果好等特点,适合大面积污染土地的修复。目前,多采用草本植物进行修复,例如鬼针草Bidens pilosa在镉质量分数为2.66 mg·kg−1的土壤生长60 d后,富集系数为4.16,对土壤中镉的去除率为4.3%~6.2%[3]。伴矿景天Sedum plumbizincicola在镉质量分数为0.55和1.85 mg·kg−1的土壤上生长至3 m时,修复效率分别为37.15%和21.82%[4]。但是草本植物存在生物量小、收割成本高、生物质后处理困难等不足,至今没有大规模的推广应用。研究表明:很多乔木可以积累重金属。张永超等[5]发现:白榆Ulmus pumila在镉质量分数为10.00 mg·kg−1的土壤中生长至7 m时,根、茎和叶的镉质量分数分别达11.2、2.42和4.74 mg·kg−1[5]。尽管乔木中重金属质量分数可能低于超积累植物,但是由于生物量大,其积累总量可能远高于超积累植物,而且重金属一旦进入乔木中,便可永久性的积累,不影响其工业用途。与超积累草本植物相比,木本植物尤其是乔木具有生物量大,不进入食物链,可以持续修复等优点,在重金属污染土壤修复中具有较大潜力[6]。然而,关于乔木吸收镉的研究大多针对器官水平,对其结构部位的研究鲜有报道。
本研究选择了在污染土壤栽植的泡桐Paulownia fortunei、楸树Catalpa bungei、悬铃木Platanus acerifolia、黑杨Populus nigra 和垂柳Salix babylonica等5种树种,分析了镉在其不同部位的分布特征,比较了不同部位的镉富集系数和镉的积累总量,旨在为合理植树造林奠定理论基础,同时也为中国镉污染土壤进行植物修复的树种选择提供依据。
-
研究区位于河南省某污染的工业用地。该地属暖温带季风气候,四季分明,气候温和,光、热、水资源丰富,年平均气温为15.5 ℃,年平均降水量为567.9 mm。拥有丰富的动植物资源及矿产资源。该地区分布有泡桐、楸树、悬铃木、黑杨和垂柳人工林纯林(于2013年前后栽植),树龄为10 a,株行距为4 m×5 m。本次采样时间为2023年7月21日。采样地泡桐、楸树、悬铃木、黑杨和垂柳土壤的镉质量分数分别为1.91、2.90、4.04、10.27和36.37 mg·kg−1。
-
随机选择污染土壤上栽植的泡桐、楸树、悬铃木、黑杨和垂柳各5株,用围尺测得胸径。采集的植物样品分为叶片、叶柄、枝皮、枝材、茎皮、茎材、根皮和根材8个部分,分别装入塑封袋。茎皮和茎材的采样高度与胸径高度一致,茎皮取样大小为4 cm×4 cm。枝皮、枝材和叶的采样高度为距地面3~5 m处的树冠。从树木东、南、西、北各个方向采集0~20 cm土壤样品,混合后装入样品袋。
-
用自来水冲洗掉附着于植物样品表面的泥土及其他污染物,然后用去离子水冲洗3遍,晾干。将样品放入烘箱105 ℃杀青30 min后,75 ℃ 下烘干至恒量。植物样品剪碎后,先用样品粉碎机初步粉碎,再用球磨机进一步磨碎,过0.149 mm筛,用来测定植物样品中镉质量分数(由于垂柳叶柄较小,难以获取,因此不作测定);土壤样品自然风干,拣出砖头、碎石、杂草等,用木棒碾碎过0.149 mm 筛,用于土壤镉质量分数测定。
土壤样品采用盐酸-硝酸-氢氟酸-高氯酸全消解法进行消解,使土壤样品中的镉全部进入试液,用石墨炉(WFX-200)测定土壤中的镉质量分数。植物样品用硝酸-过氧化氢(m硝酸∶m过氧化氢 = 10∶1)消解,经消解后植物样品用电感耦合等离子体发射光谱仪(Avio 200) 测定镉质量分数。
-
植物对镉的富集系数(F)计算公式如下:F = Ctarget/Csoil。其中:Ctarget为目标部位镉质量分数;Csoil为土壤镉质量分数。
-
植物对镉的分配计算公式如下:R = mpart/mtoatal。其中:R为相对含量;mpart为植物某部位的镉积累量; mtotal为植物中镉的积累总量。
-
采用异速生长方程估算各树种不同部位及整体生物量[7−10]。叶部、枝部、茎部和根部由于分成2个部分检测其镉质量分数,因此对两者取平均算得叶部、枝部、茎部和根部的平均镉质量分数,再估算出各树种不同部位及整体镉积累量。
-
采用SPSS 23进行数据处理及单因素方差分析,采用Origin 2022作图。
-
比较不同树种叶部的镉质量分数发现(图1A):垂柳叶片的镉质量分数最高,达10.52 mg·kg−1;其次是黑杨,达7.59 mg·kg−1;泡桐叶片的镉质量分数最低,为0.66 mg·kg−1。黑杨叶柄的镉最高,达7.45 mg·kg−1;其次是楸树,达2.04 mg·kg−1;泡桐叶柄的镉质量分数最低,为0.23 mg·kg−1。对比叶片和叶柄的镉质量分数发现:泡桐叶片镉质量分数显著大于叶柄(P<0.05),其他树木叶片和叶柄之间无显著差异。比较不同树木叶部的富集系数发现(图1B):黑杨叶片富集系数最高,达1.00;其次是楸树,达0.71;垂柳叶片的富集系数最小,为0.29。黑杨叶柄的富集系数最高,达0.98;其次是楸树,达0.50;泡桐叶柄的富集系数最小,为0.13。对比叶片和叶柄的富集系数发现:泡桐、楸树和黑杨叶片的富集系数均大于叶柄,悬铃木叶柄的富集系数大于叶片。泡桐叶片富集系数显著高于叶柄(P<0.05),悬铃木、黑杨的叶片与叶柄富集系数差异不显著。
-
比较5种树种枝部的镉质量分数发现(图2A):垂柳枝皮的镉质量分数最高,达16.08 mg·kg−1;其次是黑杨,达10.51 mg·kg−1;泡桐枝皮的镉质量分数最低,为1.24 mg·kg−1。垂柳枝材的镉质量分数最高,达2.55 mg·kg−1;泡桐枝材的镉质量分数最低,为0.19 mg·kg−1。对比枝皮和枝材的镉质量分数发现:所有树种枝皮的镉质量分数显著高于枝材的镉质量分数(P<0.05)。比较不同树种枝部的富集系数发现(图2B):黑杨枝皮的富集系数最高,达1.38;其次是楸树,达1.24;垂柳枝皮的富集系数最低,为0.44。黑杨枝材的富集系数最高,达0.32,其次是楸树、泡桐和悬铃木,分别为0.11、0.11和0.09;垂柳枝材的富集系数最低,为0.07。对比枝皮和枝材的富集系数发现:所有树种枝皮的富集系数显著高于枝材(P<0.05)。
-
比较5种树种茎部的镉质量分数发现(图3A):垂柳茎皮的镉质量分数最高,达35.30 mg·kg−1;其次是黑杨,达5.73 mg·kg−1;泡桐茎皮的镉质量分数最低,为0.12 mg·kg−1。垂柳茎材的镉质量分数最高,达2.08 mg·kg−1;其次是黑杨,达1.43 mg·kg−1;泡桐茎材的镉质量分数最低,为0.12 mg·kg−1。对比茎皮和茎材的镉质量分数发现:楸树、黑杨和垂柳茎皮的镉质量分数显著高于茎材的镉质量分数(P<0.05),泡桐和悬铃木茎皮与茎材的镉质量分数差异不显著。比较不同树种茎部的富集系数发现(图3B):垂柳茎皮的富集系数最高,达0.97;其次是楸树,达0.82;悬铃木茎皮的富集系数最低,为0.12。黑杨茎材的富集系数最高,达0.19,其次是楸树、泡桐和悬铃木,分别达0.10、0.07和0.06;垂柳茎材的富集系数最低,为0.05。对比茎皮与茎材的富集系数发现:除悬铃木外,所有树种茎皮的富集系数显著高于茎材(P<0.05)。
-
比较5种树种根部的镉质量分数发现(图4A):垂柳根皮的镉质量分数最高,达23.21 mg·kg−1;其次是黑杨,达12.01 mg·kg−1;泡桐根皮的镉质量分数最低,为0.60 mg·kg−1。垂柳根材的镉质量分数最高,达9.76 mg·kg−1;其次是黑杨,达4.50 mg·kg−1,悬铃木根材的镉质量分数最低,为0.39 mg·kg−1。对比根皮和根材的镉质量分数发现:楸树、悬铃木、黑杨和垂柳根皮的镉质量分数显著高于根材的镉质量分数(P<0.05),泡桐根皮与根材的镉质量分数差异不显著。比较不同树种根部的富集系数发现(图4B):黑杨根皮的富集系数最高,达1.58;其次是垂柳,达0.64;泡桐根皮的富集系数最低,为0.34。黑杨根材的富集系数最高,达0.59,其次是垂柳和泡桐,分别达0.27、0.25;悬铃木根材的富集系数最低,为0.10。对比根皮与根材的富集系数发现:除泡桐外,所有树种根皮的富集系数均显著高于根材的富集系数(P<0.05)。
-
利用异速生长方程计算出每株树的叶部、枝部、茎部、根部以及整株的生物量(表1),并进一步计算得到不同部位的积累量及积累总量(表2),黑杨和垂柳不同部位镉积累量从大到小依次表现为茎部、根部、枝部、叶部;楸树不同部位镉积累量从大到小依次表现为茎部、枝部、根部、叶部;泡桐和悬铃木不同部位镉积累量从大到小依次表现为枝部、茎部、根部、叶部。
表 1 不同树种不同部位及整株的生物量
Table 1. Biomass of different parts of the tree and the whole plant
树种 叶部生物量/kg 枝部生物量/kg 茎部生物量/kg 根部生物量/kg 整株生物量/kg 泡桐 6.22±0.77 a 28.59±4.73 a 31.54±6.16 bc 17.11±2.88 b 71.72±13.08 b 楸树 2.59±0.63 b 6.86±1.89 c 22.38±6.10 c 10.42±2.74 b 42.25±11.13 b 悬铃木 6.93±2.10 a 15.30±5.51 b 79.63±29.16 a 27.39±7.74 a 181.28±68.93 a 黑杨 2.95±0.66 b 11.32±3.91 bc 48.76±16.96 b 14.83±2.92 b 79.36±24.69 b 垂柳 1.79±0.38 b 7.03±2.14 c 31.64±10.55 bc 12.78±3.64 b 65.15±20.08 b 说明:不同字母表示不同树种间相同部位及整株的生物量差异显著 (P<0.05)。 表 2 不同树种不同部位及整株镉积累量
Table 2. Cd accumulation in different parts and whole trees
树种 叶部 枝部 茎部 根部 积累总量/mg 平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg泡桐 0.45 2.80±0.34 c 0.72 20.59±3.40 b 0.34 10.72±2.09 c 0.52 8.89±1.50 c 43.01±7.34 c 楸树 2.05 5.31±1.29 c 2.52 17.29±4.76 b 1.33 29.77±8.12 c 1.08 11.25±2.96 c 63.62±16.85 c 悬铃木 1.71 11.86±3.60 b 2.25 34.42±12.40 b 0.37 29.46±10.79 c 1.00 27.39±7.74 c 103.13±34.52 c 黑杨 7.52 22.19±4.99 a 6.48 73.36±25.37 a 3.58 174.57±60.71 b 8.25 122.38±24.09 b 392.51±115.15 b 垂柳 10.52 18.82±4.05 a 9.32 65.48±19.94 a 18.69 591.29±197.21 a 16.48 210.68±59.93 a 886.28±281.11 a 说明:不同字母表示不同树种间相同部位镉积累量及积累总量差异显著(P<0.05)。 比较镉在5种树种不同部位中镉的相对含量发现(图5):泡桐不同部位镉的相对含量从大到小依次为枝部、茎部、根部、叶部;楸树不同部位镉的相对含量从大到小依次为茎部、枝部、根部、叶部;悬铃木不同部位镉的相对含量从大到小依次为茎部、根部、枝部、叶部;黑杨不同部位镉的相对含量从大到小依次为枝部、根部、茎部、叶部;垂柳不同部位镉的相对含量从大到小依次为茎部、根部、枝部、叶部。总体上看,楸树、悬铃木、垂柳均表现为茎部镉的相对含量最高,叶部镉的相对含量最低。
-
本研究表明:不同树种对重金属的富集系数不同,黑杨叶片和叶柄的富集系数均为最大;黑杨枝皮和枝材的富集系数均最大;垂柳茎皮的富集系数最大,黑杨茎材的富集系数最大;黑杨根皮和根材的富集系数均最大。由此可见, 5种供试树种对镉的吸收特征不同。在这些树种中,黑杨和垂柳均有较高的镉积累量,这与前人研究结果类似[11−12],且黑杨和垂柳具有易成活、生长快、生态价值高等优势,因此在镉污染土壤修复中具有较大的推广应用价值。
积累量在重金属积累的植物中相当可观。镉超积累植物龙葵 Solanum nigrum在镉质量分数为2.00 mg·kg−1的土壤中生长86 d后积累量达0.025 mg·株−1[13];镉超积累植物商陆 Phytolacca acinosa在镉质量分数为5.00 mg·kg−1的土壤中生长60 d后地上部(镉的主要储存部分)的积累量达7.58 mg·株−1[14];镉超积累植物东南景天 Sedum alfredii在镉质量分数为2.47 mg·kg−1的土壤中生长210 d后地上部的积累量达9.60 μg·株−1[15]。树木相比超积累植物拥有更长的生长周期、更大的生物量和更高的单株镉积累量,因此具有较大的镉污染修复潜力。在树木吸收镉的研究方面,前人也做过类似的研究。徐爱春[16]研究发现:旱柳 Salix matsudana在镉质量分数为20 mg·L−1的营养液中培养一段时间后,积累量达25.81 mg·株−1。周洁等[17]研究杂交柳(耳柳Salix aurita×银柳Salix argyracea)在镉质量分数为5.28 mg·kg−1的田间土壤中生长1 a后,镉的最大积累量达73.66 mg·株−1。乔木对镉污染土壤修复后再回收利用具有一定的经济效益,比如重金属的回收、植物能源及建材方面的利用[18]。此外,乔木具有多年生、持续修复、不进入食物链的优势,因此在镉污染土壤中植树是一种较为理想的修复技术。
本研究发现:所有树种的根、茎和枝的镉质量分数均表现为树皮部大于木材部,即从大到小依次为根皮、根材、茎皮、茎材、枝皮、枝材。可见,乔木具有优先把镉积累到皮部的特性,这种现象在其他研究中也有发现。例如,唐丽清等[19]研究发现:北京市台基厂大街行道树国槐 Sophora japonica的树皮中镉质量分数高于树干。RODRÍGUEZ等[20]研究表明:西班牙发电厂附近的地中海松Pinus halepensis树皮中镉质量分数高于树干。田胜尼等[21]研究指出:镉质量分数在腺柳 Salix chaenomeloides根皮部最高,达1 438.919 μg·kg−1,乔木中最低,为228.065 μg·kg−1。这可能是乔木的解毒机制导致的。镉主要通过土壤溶液进入根系,并从木质部转移到树皮。由于树皮代谢相对旺盛[22],且树皮中富含有助于与二价离子结合的萜类、软木脂、脂肪酸酶、氨基酸等物质,这些化学物质与镉离子能够结合,生成稳定的络合物,达到解毒目的。而树干中主要为纤维素、半纤维素和木质素,与镉离子能够结合的物质较少,所以只有少部分镉积累在树干中[23]。
-
对河南某处镉污染土壤上栽植的5种树种调查发现:所有树木根皮、茎皮和枝皮均显著大于其对应的根材、茎材和枝材的镉质量分数,这可能是乔木本身的解毒机制造成的。在这些树种中,黑杨和垂柳均有较高的镉积累量,其中黑杨多数部位(除茎皮外)的富集系数均大于其他树种,因此在镉污染土壤修复中具有较大的优势。此外,采用乔木进行镉污染土壤修复时,木材中镉质量分数可达到较高水平,因此在后期木材的加工使用过程中应注意其对人体的健康风险。
Cd accumulation characteristics of different greening tree species
-
摘要:
目的 乔木对重金属具有富集能力,但是目前对其积累特征尚不明确,一定程度上限制了乔木修复重金属污染土壤技术的推广应用。探讨5种乔木绿化树种对镉的富集特征。 方法 分析了镉污染土壤上栽植的泡桐Paulownia fortunei、楸树Catalpa bungei、悬铃木Platanus acerifolia、黑杨Populus nigra 和垂柳Salix babylonica等树木不同部位对镉的富集特征,并分析了不同树种对镉的积累总量。 结果 供试树种对镉均具有一定的积累能力,单株对镉的积累量为43.01~886.28 mg·株−1,其中垂柳对镉的积累总量最大,达886.28 mg·株−1;其次是黑杨,达392.51 mg·株−1。不同树种不同部位对镉的富集能力也不同,泡桐、楸树、悬铃木、黑杨和垂柳不同部位的富集系数分别为0.07~0.56、0.10~1.24、0.06~1.04、0.32~1.58、0.06~0.97。其中黑杨多数部位(除茎皮外)的富集系数均大于其他树种。所有树种根皮、茎皮和枝皮的镉质量分数均显著(P<0.05)大于其对应的根材、茎材和枝材。所有树种不同部位镉质量分数与土壤镉质量分数均呈极显著正相关(P<0.01)。 结论 乔木主要将镉积累到枝部和茎部,其中泡桐、楸树和悬铃木主要积累在枝部,黑杨和垂柳主要积累在茎部。此外,黑杨和垂柳具有较高的镉积累量,是今后镉污染土壤修复中的首选树种。图5表2参23 Abstract:Objective While trees can accumulate heavy metals, their accumulation characteristics are currently unclear; this, to some degree, restricts the application of the technology of using trees to remediate heavy metal-contaminated soil. This study was designed to explore the cadmium (Cd) accumulation capacities of five greening tree species. Method The Cd accumulation characteristics of different parts of various greening trees (Paulownia fortunei, Catalpa bungei, Platanus acerifolia, Populus nigra, and Salix babylonica) planted in Cd-contaminated soil were analyzed, and the total amounts of Cd in the trees were determined. Result All the tested trees had a certain Cd accumulation capacity, and the cumulative amount of Cd per tree ranged from 43.01 to 886.28 mg·plant−1. Among the trees, S. babylonica had the highest accumulation amount of Cd (886.28 mg·plant−1), followed by Populus nigra (392.51 mg·plant−1). The Cd accumulation ability varied among different tree species and parts, with concentration coefficients ranging from 0.07 to 0.56, 0.10 to 1.24, 0.06 to 1.04, 0.32 to 1.58, and 0.06 to 0.97 for different parts of Paulownia fortunei, C. bungei, Platanus acerifolia, Populus nigra, and S. babylonica, respectively. The concentration coefficients of most parts (except for the stem bark) of Populus nigra were higher than those of the other trees. The Cd concentrations in the root bark, stem bark, and branch bark of all trees were significantly higher than those in their corresponding root wood, stem wood, and branch wood (P<0.05). A highly significant positive correlation was found between the Cd contents in different parts of trees and those in the soil (P<0.01). Conclusion The tested tree species accumulated Cd mainly in their branches and stems. Paulownia fortunei, C. bungei, and Platanus acerifolia accumulated Cd mainly in their branches, while Populus nigra and S. babylonica accumulated Cd mainly in their stems. Populus nigra and S. babylonica displayed higher Cd concentrations, which are the preferred tree species for the remediation of Cd-contaminated soil. [Ch, 5 fig. 2 tab. 23 ref.] -
Key words:
- trees /
- cadmium /
- total amount of accumulation /
- concentration coefficient /
- soil remediation
-
除黑龙江、吉林、内蒙古和新疆4省(区)外,中国其余省(区、市)都有竹类生长和分布。据第8次(2009−2013年)森林资源清查,中国竹林面积达601 万hm2,占中国森林面积的3%,占世界竹林面积约20%[1]。在廊道旁、公园、古寺庙、风景区等地方种植竹子以增加景观优质性,是园林配置的一部分。竹林分布广,面积大,因此需要考虑竹林保护与防火等问题。有些竹种具有一定的防火能力,被作为防火植物使用,如毛竹Phyllostachys edulis、雷竹Ph. praecox等[2-3]。对于竹子的燃烧性能方面国内外研究较少,但是在森林可燃物的燃烧性和抗火性方面,国内外都进行了大量研究[4-9]。李树华等[10]认为:在火灾危险带种植刚竹Ph. sulphurea等植物可以减缓火势蔓延。钟安建等[11]对南昌城区15种园林树种的抗火性进行研究,认为珊瑚树Viburnum odoratissimum抗火性能最强,桂花Osmanthus fragrans抗火性能最差。金钱荣等[12]将木荷Schima superba选为防火功能较强的行道绿化树种。李世友等[13]对20种园林绿化植物的鲜枝叶进行燃烧试验及燃烧性排序。何忠华等[14]对12种园林树种的抗火性进行了综合评价,认为乐昌含笑Michelia chapensis抗火性最强。森林植物叶燃烧性研究方法可以为竹叶研究提供借鉴。张雨瑶等[15]对11种园林木本植物的新叶片和2种对比植物老活叶片进行了垂直燃烧实验,认为鹅掌楸Liriodendron chinense等燃烧性较强。氧指数试验法主要用于测定聚合材料的阻燃性能,如对于各种纺织品[16]、玻璃纤维增强塑料[17]、聚氯乙烯(PVC)管[18]、橡胶[19]等阻燃性能的测定,在森林可燃物研究方面的应用较少。本研究对17种园林竹鲜叶进行燃烧性比较,旨在分析园林竹鲜叶的易燃性差异,为竹林保护与防火提供依据。
1. 材料与方法
1.1 样品采集
以17种园林竹为研究对象(括号中数字为样品代号),车筒竹Bambusa sinospinosa (1)、慈竹Neosinocalamus affinis (2)、灰金竹Phyllostachys nigra var. henonis (3)、灰香竹Chimonocalamus pallens (4)、料慈竹B. distegia (5)、龙竹Dendrocalamus giganteus (6)、绵竹B. intermedia (7)、青皮竹B. textilis (8)、沙罗单竹Schizostachyum funghomii (9)、秀叶箭竹Fargesia yuanjiangensis (10)、小佛肚竹B. ventricosa (11)、孝顺竹B. multiplex (12)、野龙竹D. semiscandens (13)、椅子竹D. bambusoides (14)、油竹B. surrecta (15)、云南甜龙竹D. hamiltonii (16)、紫竹Ph. nigra (17)。以5种常见易燃园林绿化用木本植物的老叶作对比,即阴香Cinnamomum burmanni (18)、桂花Osmanthus fragrans (19)、滇润楠Machilus yunnanensis (20)、蓝桉Eucalyptus globulus (21)、云南樟C. glanduliferum (22)。所有植物均栽植于西南林业大学校园内。由于新叶含水率呈动态变化,而老叶含水率相对稳定且易燃,故选老叶为实验样品。取叶时,选多株、不同枝条上外形和大小相似、质量相近的多片竹叶,于防火期采集健康的完整分枝,立刻带回实验室。
1.2 实验方法
采集同枝条上的老叶,分为2组,分别进行燃烧实验和含水率测定。燃烧实验前测定鲜叶质量、叶脉长度并在白纸上勾绘出鲜叶外形,实验在高浓度医用氧条件下进行,点火气体为丙烷气。将竹叶叶尖朝上、叶柄朝下放入试件夹中,点火器火焰长度为10~15 mm,从上朝下点火,用秒表记录竹叶燃烧时间。每种鲜叶重复6次实验。含水率(H)测定采用105 ℃烘干恒量法,取相对含水率。实验采用JF-3型氧指数测定仪进行。
1.3 数据获取及计算方法
叶片单位面积质量(W)、绝对线速率(V1)、绝对面积损失速率(V2)、绝对质量损失速率(V3)、相对线速率(V4)、相对面积损失速率(V5)和相对质量损失速率(V6)参照李世友等[6]、张雨瑶等[15]、郑永波等[20]和苏文静等[21]方法进行。
1.4 数据处理
运用SPSS 18.0软件,以平均V1、V2、V3、V4、V5、V6等6个指标进行因子分析,得到22种植物鲜叶的燃烧性能得分并排序。根据燃烧性能得分,应用聚类分析法划分等级。采用因子分析法对数据进行标准化处理,通过KMO值和Bartlett球体检验提取公因子,利用旋转法使因子变量更具有可解释性,计算因子变量得分。
2. 结果与分析
2.1 含水率、单位面积质量及燃烧速率
由表1可知:22种植物鲜叶的含水率和单位面积质量均差别较大,5种木本植物鲜叶单位面积质量均大于竹叶。单位面积质量最小、含水率较小的秀叶箭竹,燃烧速率最大。单位面积质量最大、含水率较大的云南樟,燃烧速率最小。含水率最大、单位面积质量较小的椅子竹,燃烧速率较小。含水率最小、单位面积质量中等的车筒竹,燃烧速率接近最大值。由此可见:鲜叶燃烧速率与单位面积质量、平均含水率有关。
表 1 22种植物鲜叶的含水率、单位面积质量及燃烧速率Table 1 Moisture content, mass per unit area and burning rate of fresh leaves of 22 plants speices代号 H/% W/(g·m−2) 绝对燃烧速率 相对燃烧速率 V1/(cm·s−1) V2/(cm2·s−1) V3/(g·s−1) V4/(%·s−1) V5/(%·s−1) V6/(%·s−1) 1 40.51 146 1.079 0.863 0.012 8.667 8.667 8.667 2 56.99 56 1.349 2.087 0.012 8.566 8.566 8.566 3 51.86 104 0.552 0.541 0.005 6.882 6.882 6.882 4 55.12 104 0.697 0.402 0.004 6.954 6.954 6.954 5 43.84 116 0.642 1.283 0.014 3.140 3.140 3.140 6 58.91 96 0.421 1.320 0.012 2.140 2.140 2.140 7 53.07 92 0.425 0.644 0.006 3.450 3.450 3.450 8 56.73 58 0.981 1.316 0.008 6.820 6.820 6.820 9 42.93 66 1.245 1.486 0.010 7.600 7.600 7.600 10 44.08 53 1.194 1.058 0.006 9.450 9.450 9.450 11 44.27 87 1.171 1.885 0.016 7.583 7.583 7.583 12 46.07 70 0.849 1.040 0.007 7.040 7.040 7.040 13 56.83 70 0.858 2.177 0.016 4.200 4.200 4.200 14 58.97 72 0.520 0.737 0.005 4.700 4.700 4.700 15 55.34 102 0.216 0.345 0.004 1.500 1.500 1.500 16 58.79 84 0.546 1.579 0.013 2.660 2.660 2.660 17 43.15 94 0.521 0.604 0.006 5.140 5.140 5.140 18 52.36 190 0.330 0.977 0.018 2.967 2.967 2.967 19 47.55 322 0.316 0.747 0.037 4.200 3.020 4.400 20 49.12 230 0.173 0.486 0.011 1.767 1.767 1.767 21 46.93 483 0.146 0.228 0.011 0.883 0.883 0.883 22 52.21 185 0.118 0.544 0.010 0.983 1.000 0.983 2.2 数据的标准化处理
由于所获得数据数值不同,单位不同,无法进行比较和计算,因此需要进行无量纲化处理。使用SPSS软件对数据进行标准化处理,结果如表2所示。
表 2 22种植物鲜叶燃烧性评价指标无量纲化后得分Table 2 Fresh leaf combustibility of 22 plants species evaluation index dimensionless points代号 V1 V2 V3 V4 V5 V6 1 1.108 53 −0.270 53 0.134 17 1.396 29 1.403 15 1.393 99 2 1.809 85 1.895 67 0.134 17 1.359 08 1.366 27 1.356 76 3 −0.260 34 −0.840 40 −0.849 73 0.738 69 0.751 40 0.735 89 4 0.116 30 −1.086 40 −0.990 29 0.765 21 0.777 69 0.762 44 5 −0.026 57 0.472 77 0.415 28 −0.639 89 −0.614 88 −0.643 72 6 −0.600 61 0.538 25 0.134 17 −1.008 29 −0.980 00 −1.012 41 7 −0.590 22 −0.658 11 −0.709 17 −0.525 68 −0.501 69 −0.529 43 8 0.853 98 0.531 17 −0.428 06 0.715 85 0.728 77 0.713 03 9 1.539 71 0.832 03 −0.146 95 1.003 20 1.013 56 1.000 61 10 1.407 24 0.074 57 −0.709 17 1.684 76 1.689 04 1.682 67 11 1.347 50 1.538 17 0.696 40 0.996 94 1.007 35 0.994 34 12 0.511 11 0.042 72 −0.568 62 0.796 90 0.809 09 0.794 14 13 0.534 49 2.054 95 0.696 40 −0.249 38 −0.227 85 −0.252 92 14 −0.343 46 −0.493 52 −0.849 73 −0.065 17 −0.045 29 −0.068 58 15 −1.133 09 −1.187 28 −0.990 29 −1.244 07 −1.213 68 −1.248 36 16 −0.275 92 0.996 62 0.274 72 −0.816 72 −0.790 14 −0.820 69 17 −0.340 86 −0.728 90 −0.709 17 0.096 92 0.115 36 0.093 65 18 −0.836 98 −0.068 78 0.977 51 −0.703 62 −0.678 05 −0.707 50 19 −0.873 34 −0.475 83 3.648 09 −0.249 38 −0.658 70 −0.179 18 20 −1.244 78 −0.937 74 −0.006 39 −1.145 71 −1.116 19 −1.149 92 21 −1.314 91 −1.394 34 −0.006 39 −1.471 38 −1.438 96 −1.475 84 22 −1.387 64 −0.835 09 −0.146 95 −1.434 54 −1.396 24 −1.438 97 2.3 KMO值和Bartlett球体检验
因子分析法并不能适用于任何情况,只有当样品数量大于评价指标数量时,才能得出KMO值和Bartlett球体检验结果,判断原始数据是否能够进行因子分析。对标准化后的数据进行KMO值和Bartlett球体检验,结果KMO值为0.625>0.500,Bartlett检验接近0,说明指标具有相关性,适合做因子分析。
2.4 公因子提取
由表3可知:特征值大于1的公因子有2个,累积方差贡献率达到了89.623%,因此可用来描述22种园林植物鲜叶的燃烧性。
表 3 解释的总方差表Table 3 Interpretation of the total variance table成分 初始特征值 提取载荷平方和 旋转载荷平方和 总计 方差百分比/% 累积/% 总计 方差百分比/% 累积/% 总计 方差百分比/% 累积/% 1 4.114 68.567 68.567 4.114 68.567 68.567 4.084 68.070 68.070 2 1.263 21.056 89.623 1.263 21.056 89.623 1.293 21.553 89.623 3 0.591 9.843 99.466 4 0.031 0.515 99.981 5 0.001 0.019 100.000 6 1.887×10−8 3.145×10−7 100.000 2.5 因子旋转
采用最大方差法(varimax)进行因子旋转,目的是使公因子的相对负荷的方差之和最大,且保持原公共因子的正交性和公共方差总和不变。使每个因子的最大载荷变量数量最小,以简化对因子的解释。利用SPSS软件进行旋转,得到表4因子载荷矩阵。主成分1在绝对线速率(V1)、相对线速率(V4)、相对面积损失速率(V5)、相对质量损失速率(V6)上的载荷系数较大,体现了燃烧性能(f1)。主成分2在绝对面积损失速率(V2)、绝对质量损失速率(V3)的载荷系数较大,体现了燃烧性能(f2)。
表 4 旋转后因子载荷矩阵Table 4 Rotated factor load matrix评价指标 主成分 评价指标 主成分 1 2 1 2 V1 0.949 0.227 V4 0.981 −0.014 V2 0.480 0.718 V5 0.990 −0.600 V3 −0.224 0.850 V6 0.979 −0.007 2.6 计算因子得分
运用SPSS软件得出因子得分系数矩阵(表5),因子得分模型可表示为:
表 5 因子得分系数矩阵Table 5 Component score coefficient matrix评价指标 主成分 评价指标 主成分 1 2 1 2 V1 0.224 0.125 V4 0.245 −0.066 V2 0.079 0.538 V5 0.250 −0.102 V3 −0.103 0.680 V6 0.244 −0.060 f1=0.224x1+0.079x2−0.103x3+0.245x4+0.250x5+0.244x6;
f2=0.125x1+0.538x2+0.680x3−0.066x4−0.102x5−0.060x6。
其中:xi为V1~V6的数值标准化后的数据,将表2的相关变量相应的代入上式中即得到22种植物鲜叶燃烧性公因子得分。再以各公因子的方差百分比作为权数计算22种植物鲜叶燃烧性综合评价得分。计算公式为:
F= λ1f1+ λ2f2=0.685 67f1+0.210 56f2。
其中:F为22种植物鲜叶的燃烧性能得分,λi为第i个公因子的方差百分比。得分大于0,说明该植物鲜叶的燃烧性能大于22种植物鲜叶燃烧性能的平均水平,反之则比较差;得分越高代表燃烧性能越好。
各植物鲜叶燃烧性能的最后得分及排名如表6所示。由表6可知:5种木本植物得分均小于0,且有2种燃烧性能得分排名最后,说明5种木本植物鲜叶的燃烧性能均低于平均水平。22种植物鲜叶的燃烧性能从大到小的顺序依次为慈竹、小佛肚竹、秀叶箭竹、沙罗单竹、车筒竹、青皮竹、孝顺竹、野龙竹、灰香竹、灰金竹、桂花、料慈竹、紫竹、椅子竹、云南甜龙竹、阴香、龙竹、绵竹、滇润楠、油竹、云南樟、蓝桉。其中,慈竹得分最高,说明最易燃,油竹得分最低,说明最难燃,但较云南樟、蓝桉易燃。具体来看,油竹的含水率较大、单位面积质量较大,在17种竹类中得分最低。秀叶箭竹含水率较小、单位面积质量最小,得分排在前列。蓝桉含水率较大、单位面积质量最大,得分排在最后。进一步说明了鲜叶的燃烧速率与单位面积质量、平均含水率有关。
表 6 22种植物鲜叶的燃烧性能得分及排序Table 6 Combustibility property score and rank of fresh leaves of 22 plants species代号 f1 f2 F 排序 代号 f1 f2 F 排序 1 1.245 11 −0.234 50 0.804 5 12 0.767 30 −0.482 73 0.424 7 2 1.546 36 1.026 31 1.276 1 13 0.031 05 1.700 27 0.379 8 3 0.510 85 −1.232 08 0.091 10 14 −0.072 37 −0.873 29 −0.233 14 4 0.609 57 −1.418 80 0.119 9 15 −1.157 88 −1.172 78 −1.041 20 5 −0.478 60 0.676 99 −0.186 12 16 −0.608 56 0.871 97 −0.234 15 6 −0.844 15 0.532 83 −0.467 17 17 0.014 46 −0.940 74 −0.188 13 7 −0.494 31 −0.792 39 −0.506 18 18 −0.807 66 0.681 30 −0.410 16 8 0.808 52 −0.063 17 0.541 6 19 −0.878 34 2.210 90 −0.137 11 9 1.168 40 0.310 35 0.866 4 20 −1.192 01 −0.405 70 −0.903 19 10 1.638 70 −0.650 80 0.987 3 21 −1.483 82 −0.585 92 −1.141 22 11 1.089 91 1.240 98 1.009 2 22 −1.412 53 −0.399 00 −1.053 21 2.7 聚类分析
应用SPSS软件对17种竹叶的燃烧性能得分进行聚类分析,由图1所示:17种园林竹鲜叶的燃烧性划为易燃和较易燃2个等级。其中,慈竹、小佛肚竹、秀叶箭竹、沙罗单竹、车筒竹、青皮竹、孝顺竹、野龙竹、灰香竹、灰金竹等易燃;料慈竹、紫竹、椅子竹、云南甜龙竹、龙竹、绵竹、油竹等较易燃。
3. 结论与讨论
对17种园林竹和5种易燃木本植物鲜叶燃烧性6个指标的因子分析可知:各植物得分差距较大,最高分与最低分之间相差2.417,说明22种植物鲜叶的燃烧性差距较大。与5种园林木本植物相比,竹叶均为易燃叶。料慈竹、椅子竹、云南甜龙竹、龙竹、紫竹、绵竹和油竹的鲜叶燃烧性能相对较低,尤其是油竹,比桂花、阴香和滇润楠还难燃。绝对线速率(V1)、相对线速率(V4)、相对面积损失速率(V5)和相对质量损失速率(V6)对其燃烧性影响较大。基于17种竹的燃烧性能得分,SPSS聚类分析将其划为易燃和较易燃2个等级,其中易燃竹种10种,较易燃竹种7种。
鲜叶的燃烧性受自身理化性质和生态学、生物学特性等多因素的综合影响。昆明地区旱季降雨稀少,园林竹浇水较为频繁,浇水周期、浇水量和浇水次数对竹叶的含水率造成一定影响。施肥也会影响竹子生理性能。研究表明施氮肥会提高大豆Glycine max的脂肪含量[22-23];不同磷含量培养液处理下植株幼苗的株高、茎叶生物量和总生物量差异极其显著[24];不同磷源处理下云南松Pinus yunnanensis幼苗体内磷含量明显不同[25]。施肥对植物化学成份的影响一定程度上也影响其燃烧性。本研究中的竹叶样品采自竹下较低部位;竹子受自身生长因素及光照等外部因素影响,不同空间部位的竹叶生长发育不均衡,也会导致竹叶不同的理化性质和生态学特性。以后的研究中,要尽量减少人工经营措施对实验取样的干扰,并且考虑不同空间部位对竹叶的作用,使样品更具有代表性。本研究根据竹叶的燃烧速率来分析燃烧性,而没有分析理化性质、生态学特性等对燃烧性的影响。因此,以上鲜叶的燃烧性排序及分类是在特定条件下得出的,能否适用于其他条件还需要进一步验证。
-
表 1 不同树种不同部位及整株的生物量
Table 1. Biomass of different parts of the tree and the whole plant
树种 叶部生物量/kg 枝部生物量/kg 茎部生物量/kg 根部生物量/kg 整株生物量/kg 泡桐 6.22±0.77 a 28.59±4.73 a 31.54±6.16 bc 17.11±2.88 b 71.72±13.08 b 楸树 2.59±0.63 b 6.86±1.89 c 22.38±6.10 c 10.42±2.74 b 42.25±11.13 b 悬铃木 6.93±2.10 a 15.30±5.51 b 79.63±29.16 a 27.39±7.74 a 181.28±68.93 a 黑杨 2.95±0.66 b 11.32±3.91 bc 48.76±16.96 b 14.83±2.92 b 79.36±24.69 b 垂柳 1.79±0.38 b 7.03±2.14 c 31.64±10.55 bc 12.78±3.64 b 65.15±20.08 b 说明:不同字母表示不同树种间相同部位及整株的生物量差异显著 (P<0.05)。 表 2 不同树种不同部位及整株镉积累量
Table 2. Cd accumulation in different parts and whole trees
树种 叶部 枝部 茎部 根部 积累总量/mg 平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg平均镉质量分
数/(mg·kg−1)积累
量/mg泡桐 0.45 2.80±0.34 c 0.72 20.59±3.40 b 0.34 10.72±2.09 c 0.52 8.89±1.50 c 43.01±7.34 c 楸树 2.05 5.31±1.29 c 2.52 17.29±4.76 b 1.33 29.77±8.12 c 1.08 11.25±2.96 c 63.62±16.85 c 悬铃木 1.71 11.86±3.60 b 2.25 34.42±12.40 b 0.37 29.46±10.79 c 1.00 27.39±7.74 c 103.13±34.52 c 黑杨 7.52 22.19±4.99 a 6.48 73.36±25.37 a 3.58 174.57±60.71 b 8.25 122.38±24.09 b 392.51±115.15 b 垂柳 10.52 18.82±4.05 a 9.32 65.48±19.94 a 18.69 591.29±197.21 a 16.48 210.68±59.93 a 886.28±281.11 a 说明:不同字母表示不同树种间相同部位镉积累量及积累总量差异显著(P<0.05)。 -
[1] 李杨, 丁振军, 王星蒙. 葫芦岛地区土壤镉元素空间含量分布特征及污染评价研究[J]. 环境科学与管理, 2022, 47(10): 166 − 170. LI Yang, DING Zhenjun, WANG Xingmeng. Content distribution characteristics and pollution evaluation of cadmium in soil, Huludao area [J]. Environmental Science and Management, 2022, 47(10): 166 − 170. [2] 顾继光, 周启星, 王新. 土壤重金属污染的治理途径及其研究进展[J]. 应用基础与工程科学学报, 2003, 11(2): 143 − 151. GU Jiguang, ZHOU Qixing, WANG Xin. Reused path of heavy metal pollution in soils and its research advance [J]. Journal of Basic Science and Engineering, 2003, 11(2): 143 − 151. [3] 张云霞, 周浪, 肖乃川, 等. 鬼针草(Bidens pilosa L. )对镉污染农田的修复潜力[J]. 生态学报, 2020, 40(16): 5805 − 5813. ZHANG Yunxia, ZHOU Lang, XIAO Naichuan, et al. Remediation potential of B. pilosa L. in cadmium-contaminated farmland [J]. Acta Ecologica Sinica, 2020, 40(16): 5805 − 5813. [4] 柏佳, 谭长吟, 曹雪莹, 等. 3种有机酸对伴矿景天修复效率及土壤微生物数量的影响[J]. 水土保持学报, 2020, 34(2): 318 − 324. BAI Jia, TAN Changyin, CAO Xueying, et al. Effect of three organic acids on the remediation efficiency of Sedum plumbizincicola and soil microbial quantity [J]. Journal of Soil and Water Conservation, 2020, 34(2): 318 − 324. [5] 张永超, 陈存根, 钟发明, 等. 重金属铅、镉在白榆中分布规律和累积特性研究[J]. 西北林学院学报, 2011, 26(5): 6 − 11. ZHANG Yongchao, CHEN Cungen, ZHONG Faming, et al. Distribution and accumulation of lead and cadmium in Ulmus pumila [J]. Journal of Northwest Forestry University, 2011, 26(5): 6 − 11. [6] 曲豪杰, 张涵丹, 谭志超, 等. 纳塔栎和柳叶栎对铅锌矿区污染土壤的修复潜力分析: 田间试验[J]. 生态学报, 2020, 40(19): 6969 − 6981. QU Haojie, ZHANG Handan, TAN Zhichao, et al. The remediation potential of Quercus nuttallii and Quercus phellos to heavy metalpolluted soil in lead-zinc mining area: a field test [J]. Acta Ecologica Sinica, 2020, 40(19): 6969 − 6981. [7] 李培坤, 牛帅, 王雪颖, 等. 河南农业大学文化路校区校园植物调查与分析[J]. 河南农业大学学报, 2017, 51(2): 243 − 249. LI Peikun, NIU Shuai, WANG Xueying, et al. lnvestigation and analysis of campus plants in Wenhua Road Campus of Henan Agricultural University [J]. Journal of Henan Agricultural University, 2017, 51(2): 243 − 249. [8] 张欣. 基于树木年轮信息的北京园林树木碳固定研究[D]. 保定: 河北农业大学, 2021. ZHANG Xin. Research on Landscape Trees Carbon Fixation in Beijing Based on Tree Rings Information [D]. Baoding: Hebei Agricultural University, 2021. [9] 罗云建, 王效科, 逯非. 中国主要林木生物量模型手册[M]. 北京: 中国林业出版社, 2015. LUO Yunjian, WANG Xiaoke, LU Fei. Comprehensive Database of Biomass Regressions for China’s Tree Species [M]. Beijing: China Forestry Publishing House, 2015. [10] 周国逸, 尹光彩, 唐旭利, 等. 中国森林生态系统碳储量——生物量方程[M]. 北京: 科学出版社, 2018. ZHOU Guoyi, YIN Guangcai, TANG Xuli, et al. Carbon Storage-Biomass Equation of Forest Ecosystems in China [M]. Beijing: Science Press, 2018. [11] 张晓丽, 翟飞飞, 李伟, 等. 27个柳树无性系对镉的吸收分配特性[J]. 林业科学, 2017, 53(4): 9 − 17. ZHANG Xiaoli, ZHAI Feifei, LI Wei, et al. Characteristics of cadmium absorption and distribution in 27 willow clones [J]. Scientia Silvae Sinicae, 2017, 53(4): 9 − 17. [12] 郑青雅, 马云璐, 戎恭, 等. 2年生杨树吸收和富集土壤镉特性研究[J]. 安徽农业大学学报, 2023, 50(2): 297 − 302. ZHENG Qingya, MA Yunlu, RONG Gong, et al. Study on the uptake and enrichment of heavy metal Cd in two-year old poplar trees [J]. Journal of Anhui Agricultural University, 2023, 50(2): 297 − 302. [13] 李志贤, 冯涛, 陈章, 等. 镉胁迫对龙葵镉的吸收积累及生理响应的影响[J]. 水土保持学报, 2017, 31(5): 328 − 333. LI Zhixian, FENG Tao, CHEN Zhang, et al. Effects of different levels Cd stress on Cd uptake and physiological response of Solanum nigrum L. [J]. Journal of Soil and Water Conservation, 2017, 31(5): 328 − 333. [14] 聂发辉, 吴彩斌, 吴双桃. 商陆对镉的富集特征[J]. 浙江林学院学报, 2006, 23(4): 400 − 405. NIE Fahui, WU Caibin, WU Shuangtao. Characteristics of Phytolacca acinosa accumulating Cd [J]. Journal of Zhejiang Forestry College, 2006, 23(4): 400 − 405. [15] 刘波, 陈绩, 马嘉伟, 等. 交流电场和水分对褪色柳-东南景天混栽修复镉污染土壤的影响[J]. 浙江农林大学学报, 2021, 38(6): 1238 − 1244. LIU Bo, CHEN Ji, MA Jiawei, et al. Effects of AC electric field and soil moisture on phytoremediation of cadmium contaminated soil by mixed planting of Salix discolor and Sedum alfredii [J]. Journal of Zhejiang A&F University, 2021, 38(6): 1238 − 1244. [16] 徐爱春. 柳树无性系镉积累和生理变化规律研究[D]. 北京: 中国林业科学研究院, 2007. XU Aichun. Study on Cd Accumulation and Physiological Change Rule in Clones of Salix [D]. Beijing: China Academy of Forestry, 2007. [17] 周洁, 陈庆生, 王保松, 等. 镉超富集柳树对镉吸收、积累及分配特性的田间试验研究[J]. 江苏林业科技, 2023, 50(1): 1 − 6. ZHOU Jie, CHEN Qingsheng, WANG Baosong, et al. Characteristics of cadmium absorption, accumulation and distribution in cadmium-tolerant willows of field test [J]. Journal of Jiangsu Forestry Science &Technology, 2023, 50(1): 1 − 6. [18] 钟斌, 陈俊任, 彭丹莉, 等. 速生林木对重金属污染土壤植物修复技术研究进展[J]. 浙江农林大学学报, 2016, 33(5): 899 − 909. ZHONG Bin, CHEN Junren, PENG Danli, et al. Research progress of heavy metal phytoremediation technology of fast-growing forest trees in soil [J]. Journal of Zhejiang A&F University, 2016, 33(5): 899 − 909. [19] 唐丽清, 邱尔发, 韩玉丽, 等. 不同径级国槐行道树重金属富集效能比较[J]. 生态学报, 2015, 35(16): 5353 − 5363. TANG Liqing, QIU Erfa, HAN Yuli, et al. Variation of heavy metal enrichment efficiency in roadside trees of Sophora japonica L. with different diameters at breast height [J]. Acta Ecologica Sinica, 2015, 35(16): 5353 − 5363. [20] RODRÍGUEZ M, JOSE A, GUTIERREZ, et al. Wood and bark of Pinus halepensis as archives of heavy metal pollution in the Mediterranean Region [J]. Environmental Pollution, 2018, 239: 438 − 447. [21] 田胜尼, 张静, 孙庆业, 等. 铜尾矿自然定居腺柳对重金属吸收及分布的研究[J]. 农业环境科学学报, 2013, 32(9): 1771 − 1777. TIAN Shengni, ZHANG Jing, SUN Qingye, et al. Heavy metals absorption and distribution by Salix chaenomeloides settled naturally on the copper tailings [J]. Journal of Agro-Environment Science, 2013, 32(9): 1771 − 1777. [22] 王荣芬, 邱尔发, 唐丽清, 等. 行道树毛白杨树干中重金属元素分布[J]. 生态学报, 2014, 34(15): 4212 − 4222. WANG Rongfa, QIU Erfa, TANG Liqing, et al. Heavy metals distribution in roadside tree trunks under a heavy urban transportation condition [J]. Acta Ecologica Sinica, 2014, 34(15): 4212 − 4222. [23] GONZALEZ-OREJA J A, ROZAS M A, ALKORTA I, et al. Dendroremediation of heavy metal polluted soils [J]. Reviews on Environmental Health, 2008, 23(3): 223 − 234. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20230630