基于光谱分析探讨磷酸与焦磷酸改性生物炭的磷素形态
王旭阳1, 孙涛1, 朱新萍1, 唐光木2, 贾宏涛1,*, 徐万里2
1.新疆农业大学资源与环境学院, 新疆 乌鲁木齐 830052
2.新疆农业科学院土壤肥料与农业节水研究所, 新疆 乌鲁木齐 830091
*通讯作者 e-mail: jht@xjau.edu.cn

作者简介: 王旭阳, 1996年生, 新疆农业大学资源与环境学院硕士研究生 e-mail: 1515375795@qq.com

摘要

通过磷酸(H3PO4)和焦磷酸(H4P2O7)对生物炭改性能够使其更适于农业应用。 探明H3PO4和H4P2O7改性生物炭的P赋存形态与结合方式, 将有助于揭示其表面P的生物有效性。 以麦秆生物炭(WBC)与棉秆生物炭(CBC)为原料, 分别通过H3PO4和H4P2O7制备了H3PO4改性生物炭(P-WBC和P-CBC)和H4P2O7改性生物炭(PA-WBC和PA-CBC)。 利用拉曼光谱(Raman)与扫描电镜能谱(SEM-EDS)对改性生物炭结构与P分布变化进行表征, 采用傅里叶红外光谱(FTIR)和X射线光电子能谱(XPS)探究改性生物炭表面P结合方式, 并结合Hedley磷分级方法与可见分光光度法, 定量分析改性前后生物炭中P形态及含量变化。 结果表明, H3PO4和H4P2O7改性后生物炭 IG/ ID值增大, 石墨化结构增强, 形成了含P颗粒状结构。 H3PO4和H4P2O7改性促进了生物炭表面羧基(—COOH)、 P—O—P和P—H等酸性官能团与含P基团的形成, 且H3PO4改性生物炭和H4P2O7改性生物炭表面官能团种类相似。 XPS结果显示, 与WBC和CBC相比, 改性处理中的O(1 s)峰相对含量显著增加了13.15%~32.44%, P(2 s)峰相对含量显著增加了18.54%~27.02%( p<0.05)。 反褶积分峰将P(2 s)与O(1 s)分为C—P—O, C—O—P, O=P—O C=O与(或) P=O C—O—C与(或)P—O—C和P—O—P六类。 较H3PO4改性而言, H4P2O7改性能够促进更多C—O—P, O=P—O C—O—C与(或)P—O—C和P—O—P键的形成。 改性也使得生物炭中总P含量显著增加, 且PA-WBC和PA-CBC中P含量显著高于P-WBC和P-CBC。 与WBC和CBC相比, 改性处理中活性P含量显著提高2.36~14.77 g·kg-1, 稳定态P含量显著降低0.06~0.17 g·kg-1( p<0.05)。 与P-WBC和P-CBC相比, PA-WBC和PA-CBC的活性P、 中等活性P分别显著增加了5.27~15.66和0.53~0.64 g·kg-1, 稳定态P含量减少了0.03~0.34 g·kg-1( p<0.05)。 H3PO4和H4P2O7改性改变了P在生物炭表面的结合方式, 同时增加了P的活性。 H3PO4和H4P2O7改性生物炭间, 不同形态P含量和结合方式的差异对进一步探究P的生物有效性具有重要意义。

关键词: 生物炭; 改性; 磷形态; 光谱分析
中图分类号:O657.3 文献标志码:A
Phosphorus Species of Biochar Modified by Phosphoric Acid and Pyrophosphoric Acid Based on Spectral Analysis
WANG Xu-yang1, SUN Tao1, ZHU Xin-ping1, TANG Guang-mu2, JIA Hong-tao1,*, XU Wan-li2
1. College of Resources and Environment Science, Xinjiang Agricultural University, Urumqi 830052, China
2. Institute of Soil and Fertilizer & Agricultural Sparing Water, Xinjiang Academy of Agricultural Sciences, Urumqi 830091, China
*Corresponding author
Abstract

They can be applied to agriculture by changing the properties of biochar by phosphoric acid (H3PO4) and pyrophosphoric acid (H4P2O7) . It is helpful to revealing the bio-availability of P on its surface for identifying P occurrence form and binding mode of changed biochar for H3PO4 and H4P2O7. This paper adopts wheat stalk biochar (WBC) and cotton stalk biochar (CBC) as raw materials. Meanwhile, this paper prepares changed H3PO4 (P-WBC and P-CBC) and changed biochar of H4P2O7 (PA-WBC and PA-CBC) by H3PO4 and H4P2O7 respectively. This paper adopts Raman spectroscopy (Raman) and scanning electron microscopy (SEM) to characterize the structure and P distribution of the changed biochar. In addition, this paper adopts Fourier transform infrared spectroscopy (FTIR) and X-ray photoelectron spectroscopy (XPS) to investigate the P-binding mode of changed biochar surface. Meanwhile, this paper quantitatively analyzes the changes of P form and content in biochar before and after changing by combining the Hedley phosphorus classification method with visible spectro-photometry. The result shows that the IG/ ID value of biochar increases and the graphitization structure enhances to form the P-containing granular structure after changing the properties of H3PO4 and H4P2O7. It has promoted the formation of carboxyl (—COOH), P—O—P and P—H acid functional groups on the surface of biochar for changes of H3PO4 and H4P2O7. They are similar to functional groups on the surface of changed biochar of H3PO4 and changed biochar of H4P2O7. XPS result shows that it increases by 13.15%~32.44% significantly by compared with WBC and CBC for the relative content of O1s peak in the changed treatment. Meanwhile, it also shows that it increases by 18.54%~27.02% significantly for the relative content of the P(2 s) peak ( p<0.05). They have divided P(2 s) and O(1 s) into C—P—O, C—O—P, O=P—O C=O and (or) P=O C—O—C and (or) P—O—C and P—O—P for the deconvolution integral peaks. It can promote the formation of C—O—P, O=P—O C—O—C and/or P—O—C and P—O—P bonds for the changed properties of H4P2O7 by comparing with the changed properties of H3PO4. It also significantly increases the total P content in biochar for the changed properties. Meanwhile, it is significantly higher than that in P—WBC and P—CBC for the P content in PA-WBC and PA-CBC. The active P content in the changed treatment significantly increases by 2.36~14.77 g·kg-1 compared with WBC and CBC. In addition, the stable P content significantly decreases by 0.06~0.17g·kg-1 ( p<0.05). The active P and moderate active P of PA-WBC and PA-CBC significantly increase by 5.27~15.66 and 0.53~0.64 g·kg-1 respectively. The stable P content decreases by 0.03~0.34 g·kg-1 ( p<0.05). In conclusion, it has changed the binding mode of P on biochar surface and increased the activity of P for the changed properties of H3PO4 and H4P2O7. The difference in the content and binding mode of P in different forms between H3PO4 and H4P2O7 modified biochar are of great significance for further exploring the bioarailability of P.

Keyword: Biochar; Modification; Phosphorus species; Spectral analysis
引言

生物炭是一种环境友好型材料, 在土壤培肥、 废水处理等多个领域得到广泛应用。 而生物炭常导致环境pH值升高并引入大量阳离子, 使得其在一些环境中的使用受到了严格限制[1]。 采取物理或化学的方法对生物炭进行改性已经成为减少生物炭在应用中受限的关键环节。 磷(P)较小的原子半径和价电子可极大的调节炭材料性质[2], 因此正磷酸(H3PO4)与焦磷酸(H4P2O7)是对生物炭化学改性的良好选择。 H3PO4和H4P2O7在降低生物炭pH值的同时增加生物炭中速效磷含量[3], 使H3PO4和H4P2O7改性生物炭具有作为新型酸性磷肥的潜力。 了解H3PO4和H4P2O7改性生物表面P形态与结合方式对于深入探究其生物可利用性具有重要意义。

P是影响作物生长的重要营养元素, 土壤中P的形态与含量是决定其植物有效性的关键因素。 长期以来, 为准确定量样品中P各形态含量, 多采用顺序提取的方法。 但顺序提取法不能确切反映P的结合方式, 在预测P形态结构演化方面存在偏差。 光谱分析在P结合方式的分析方面得到了广泛应用。 扫描电子显微镜-能谱仪(SEM-EDS)能够反映样品表面形貌及P素的表面分布情况。 X射线光电子能谱(XPS)通过X射线光源, 获得电子结合能, 为P分子结构和原子价态分析提供依据[4]。 拉曼光谱(Raman)和傅里叶变换红光谱(FTIR)等多种光谱技术的联用, 能够在一定程度上实现互补, 从中获取分子结构变化信息[5]。 但光谱分析技术多局限于半定量分析, 对不同P结合方式的定量分析涉及较少。 前人对P的研究方法多集中于顺序提取法或光谱分析方面, 关于综合运用顺序提取法和光谱分析技术, 分析P结合方式与形态间关系的研究少有报道。

因此, 本文基于H3PO4和H4P2O7改性的麦秆与棉秆生物炭。 通过SEM-EDS, Raman, FTIR和XPS综合探究P的结构演化。 并利用P顺序分馏提取法对生物炭中的P形态及有效性进行分析。 拟揭示H3PO4和H4P2O7改性生物炭表面P结合方式与有效性间的关系, 以期为酸性生物炭的应用提供基础数据。

1 实验部分
1.1 生物炭的制备

试验设置麦秆生物炭(WBC)与棉秆生物炭(CBC)两种生物炭种类处理, 生物炭均在温度为400 ℃的厌氧环境下热解2 h制备, 实验材料由新疆农业科学院提供。 称取一定量原始生物炭(WBC、 CBC)洗净、 烘干、 粉碎后过10目筛后待用。 将浓度为70 Wt%的H3PO4于石英坩埚中, 在205 ℃马弗炉中脱水缩合1.5 h至液体浓稠, 制备H4P2O7。 将原始生物炭以20:1(mL· g-1)的比例分别于H3PO4和H4P2O7中浸渍, 用橡胶塞密封容器后, 放置在25 ℃恒温震荡箱内, 120 r· min-1下震荡24 h。 使用去离子水冲洗样品至滤液pH稳定, 将冲洗完成的样品放入65 ℃烘箱中干燥至恒重, 制备出H3PO4改性生物炭(P-WBC, P-CBC)和H4P2O7改性生物炭(PA-WBC, PA-CBC)。

1.2 测定方法

各处理C, H, N及O元素采用元素分析仪(PerkinElmer2400II, 美国)测定。 使用扫描电镜能谱(SEM-EDS, FEI F50-Burker 6-10 XFlash, 美国)观察生物炭表面形貌特征及P素分布。 采取拉曼光谱(Raman, DXR2XI, 美国)分析生物炭改性后的表面缺失程度, 测试条件为: 波长488 nm, 激光光斑1 mm, 扫描波数范围200~1 700 cm-1, 分辨率1 cm-1。 通过傅里叶红外光谱(FT/IR-4100, 日本)在波数为4 000~400 cm-1条件下进行生物炭表面官能团测定。 生物炭的化学成分采用ThermoFicher ESCALAB 250Xi型X射线光电子能谱(XPS)仪进行分析, 全谱测试条件: 通能160 eV, 15 kV, 5 mA; 步长为1 eV, 测试面积为700 μ m× 300 μ m。 分谱测试条件: 通能20 eV, 15 kV, 10 mA, 步长为0.05 eV。 通过XPSPEAK 4.1对光谱进行分析。 生物炭P形态采用Hedley磷分级体系, 采取连续浸提的方法[6], P各形态定量分析使用可见分光光度计测定(Spectrumlab 22pc, 中国)。

2 结果与讨论
2.1 生物炭基本性质

2.1.1 元素分析

各处理基本化学性质如表1所示。 H3PO4和H4P2O7改性均使WBC和CBC处理中的碳(C)元素含量显著减少(p< 0.05), 这可能归因于酸对碳结构边缘的破坏。 与WBC和CBC相比, P-WBC与PA-WBC处理氢(H)含量分别显著增加了0.27%和0.71%, P-CBC与PA-CBC处理H含量分别增加了0.99%和1.09%(p< 0.05)。 相较于P-WBC和P-CBC, PA-CBC与PA-WBC处理灰分含量显著增加, 较高的灰分含量可能是造成H4P2O7改性生物炭中C含量显著降低的重要原因[7]。 改性后麦秆与棉秆生物炭中总磷含量分别显著提高了10.90~17.24和15.99~32.51 g· kg-1, 并且H4P2O7改性处理显著高于H3PO4改性处理(p< 0.05), 这说明P被成功负载至生物炭表面, 且提升量较大。

表1 生物炭基本性质 Table 1 Basic properties of biochar

2.1.2 拉曼光谱(Raman)分析

麦秆与棉秆生物炭改性前后的Raman图谱分别见图1(a, b), 在1 600和1 350 cm-1处有两个明显的散射峰, 分别对应G峰和D峰。 D峰代表材料炭结构中的缺陷程度和无定型碳的无序程度, G峰归因于石墨中平面内的振动模式, 即反映了有序碳原子的有序度, G峰和D峰的强度比(IG/ID)可用于表示类碳材料的碳化程度, 其与IG/ID成正比[1]。 经计算WBC, P-WBC和PA-WBC的IG/ID分别为0.36, 0.40和0.41; CBC, P-CBC和PA-CBC的IG/ID分别为0.31, 0.37和0.40。 在同一改性处理下, 麦秆生物炭IG/ID均高于棉秆生物炭, 表明H3PO4和H4P2O7改性使生物炭的有序度增加, 石墨化程度增强; 在同一改性处理下, 棉秆生物炭具有更多缺陷和无定型碳。

图1 麦秆(a)与棉秆(b)生物炭改性前后的拉曼光谱
P-与PA-分别代表经过H3PO4改性和H4P2O7改性, 下同
Fig.1 Raman spectra of the wheat straw (a) and cotton stalk (b) biochars before and after modification
P- and PA- represent modified by H3PO4 and H4P2O7, respectively, Same below

2.2 生物炭P素表征

2.2.1 扫描电镜能谱(SEM-EDS)分析

通过SEM-EDS表征生物炭改性前后表面形貌及相应P素分布变化, 结果如图2(a, b)所示。 可以看出, WBC处理表面粗糙, 呈现出密集的结构特征, 这可能来自热解过程中生物质自身纤维素与盐类物质的残留[8]。 与WBC和CBC处理相比, 改性后生物炭表面P素显著增加。 改性后生物炭表面形成了P含量较高的颗粒状结构, 且PA-WBC和PA-CBC处理表面的颗粒状结构多于P-WBC和P-CBC处理, 这可能归因于磷酸盐与焦磷酸盐在生物炭表面的形成与积累, 为改性生物炭表面P素增加提供了有力证据。

图2 麦秆(a)与棉秆(b)生物炭改性前后SEM图像与相应的P素映射图像Fig.2 SEM image and the corresponding P elemental mapping images of the wheat straw (a) and cotton stalk (b) biochars before and after modification

2.2.2 傅里叶红外光谱(FTIR)分析

图3(a, b)分别为麦杆与棉杆各处理生物炭表面的官能团变化。 结果显示, 在490~2 420 cm-1之间可以观察到多个信号, 与WBC和CBC处理相比, 改性后的生物炭在490, 799, 1 085, 1 420和2 416 cm-1处的特征峰发生改变。 其中490 cm-1对应P— (OH)3的变角振动[9]、 799 cm-1对应C— P的伸缩振动、 1 085 cm-1归因于P+— O的电离和聚磷酸链中的P— O— P的对称振动[10]、 1 420 cm-1对应羧基(— COOH)基团的形成、 2 416 cm-1归因于PH3或P— H基团的形成[10]。 相较于WBC和CBC处理, H3PO4和H4P2O7改性有效促进了P— (OH)3键的形成, 破坏了WBC和CBC处理中C— P键, 并且形成了新的含P基团。 麦秆与棉秆生物炭通过H3PO4改性后的红外光谱与H4P2O7改性后的红外光谱变化相近, 这表明通过H3PO4和H4P2O7改性后的生物炭表面官能团相似。 红外光谱分析显示, 改性使生物炭表面含P基团增加。

图3 麦秆(a)与棉秆(b)生物炭改性前后的傅里叶红外光谱Fig.3 FTIR spectra of the wheat straw (a) and cotton stalk (b) biochars before and after modification

2.2.3 X射线光电子能谱(XPS)分析

图4为麦杆与棉杆各处理生物炭的XPS图谱。 可以看出, XPS总谱图大约在132, 286和534 eV处显示了三个主要特征峰, 分别对应P(2p), C(1s)和O(1s)峰。 与WBC和CBC相比, 改性后生物炭的C1s峰相对含量下降了42.22%~59.46%, 这一结果可能归因于生物炭边缘结构不稳定, 即H3PO4和H4P2O7作为一种氧化剂, 导致生物碳表面C原子被O或P取代, 或以C与酸反应生成的CO或CO2的形式丢失[10]。 改性后O(1s)峰相对强度增加了13.15%~32.14%, 这可能是由生物炭表面COOH和磷酸盐基团的形成所导致[11]。 并且P(2p)特征峰相对强度增加了2.02%~8.48%, 进一步说明H3PO4和H4P2O7改性使P素在生物炭表面成功结合, 并促进了P, O和C间化学键的形成。

图4 麦秆(a)与棉秆(b)生物炭改性前后的XPS光谱及特征峰相对含量Fig.4 XPS spectra and the relative content of characteristic peaks of wheat straw (a) and cotton stalk (b) biochars before and after modification

为进一步了解各生物炭中P的结合形式, 对P(2p)与O(1s)的高分辨率特征峰进行了反褶积分析, 得到各组分的相对含量。 反褶积分结果分别如图5(a— d)所示。 P(2p)峰约在133.1 eV(C— P— O), 134.0 eV(C— O— P)和135.0 eV(O=P— O)出现3处主峰[12]。 O(1s)特征峰被反褶积分为3个峰, 分别是位于约531.00 eV的C=O与(或) P=O峰, 约532.20 eV的C— O— C与(或)P— O— C峰, 以及位于约533.5 eV的P— O— P峰[10]。 分峰结果显示, 在WBC和CBC处理中, 未能观察到P(2p)特征峰。 改性后, PA-WBC与PA-CBC处理中C— O— P相对含量最高, 分别达到55.81%和50.19%, 高于P-WBC与P-CBC处理17.58%和10.73%。 与H3PO4改性处理相比, PA-WBC与PA-CBC处理中C— P— O的含量分别降低了22.87%和25.06%。 P元素位于第三周期, 其共价半径(107± 3 pm)显著高于C(73 pm), 使得P更倾向于形成sp3构型[10], 加之C原子易被P原子替代[2], 使得H3PO4和H4P2O7促进了更多由O桥接的P与生物炭中的C结合。 C— O— P键中O— P的存在也证明了在生物炭表面有磷酸盐与焦磷酸盐结构的形成[13]

图5 麦秆与棉秆生物炭改性前后的高分别率P(2p)(a, b)和O(1s)(c, d)图谱Fig.5 High-resolution P(2p) (a, b) and O(1s) (c, d) spectra of the wheat straw and cotton stalk biochars before and after modification

O(1s)特征峰的反褶积结果显示, 与WBC和CBC相比, P-WBC, PA-WBC, P-CBC和PA-CBC处理中P— O— P相对含量提高了15.98%~25.29%, PA-WBC和PA-CBC处理中P— O— P相对含量分别高于P-WBC和P-CBC处理3.64%和3.66%。 C=O与(或) P=O峰相对含量减少了9.86%~32.86%。 H3PO4和H4P2O7改性使C=O与(或) P=O和C— O— C与(或)P— O— C发生了正化学位移, 结合能增加, 这意味着电荷从C原子向P原子转移, 形成了有关P的化学键[10]。 通过H3PO4和H4P2O7改性后, 生物炭形成大量含P基团, 促进了电荷从C原子向P原子转移。 由此可推测, 改性后生物炭中的O原子更多趋向于与P原子的结合。 有研究发现, 与C— P键相比, O— P体系的稳定性较差, 键能表现为最弱, O— P键易发生断裂[3], 从而释放含P基团。 PA-WBC与PA-CBC处理中O— P的结合方式相对含量较高, 这可能造成其能够释放更多的含P基团, 进而增加P的活性。 XPS分析结果证明, H3PO4和H4P2O7改性丰富了生物炭表面P的结合方式, H4P2O7改性处理中P的活性可能更强。

2.2.4 生物炭P形态变化

图6反映了不同处理中各形态P素的相对含量, Hedley根据植物对各形态P素吸收利用的难易程度, 将其分为活性P(H2O-Pi, NaHCO3-Pi和NaHCO3-Po)、 中等活性P(NaOH-Pi和NaOH-Po)、 稳定态P(HCl-P)和残留态P(Rs-P)[14]。 与WBC相比, P-WBC与PA-WBC中活性P显著提高了2.36~5.91 g· kg-1, 中等活性P显著提高了0.47~0.85 g· kg-1, 稳定态P降低了0.14~0.17 g· kg-1(p< 0.05)。 与CBC相比, P— CBC与PA-CBC中活性P显著提高了4.27~14.78 g· kg-1, 中等活性P显著提高了0.49~0.90 g· kg-1, 稳定态P显著降低了0.89~1.03 g· kg-1(p< 0.05)。 其中, 稳定态P的减少可归因于稳定态P向活性态P的转化[15]。 总体来看, WBC和CBC中P素有效性相对较低, 活性P占总P含量的41.21%与41.13%。 改性后, P-WBC和PA-WBC处理中, 活性P相对含量分别增加了42.05%与42.33%; P-CBC和PA-CBC处理中, 活性P相对含量分别提高了44.05%与47.81%。 由此可见, 改性后生物炭中P的活性提高, 且H4P2O7改性处理中的活性P含量最高, 这与XPS分析结果相一致。

图6 麦秆与棉秆生物炭改性前后磷的相对含量变化Fig.6 Changes of phosphorus relative content of the wheat straw (a) and cotton stalk (b) biochars before and after modification

3 结论

(1)通过H3PO4和H4P2O7改性后, 麦秆与棉秆生物炭C含量显著减少, H含量显著增加, P素大量负载, 生物炭结构有序度增加(p< 0.05)。

(2)通过H3PO4和H4P2O7改性后, 麦秆与棉秆生物炭表面形成了含P颗粒状结构, 酸性官能团与含P基团增加, 形成了P(2p)峰。 与H3PO4改性处理相比, H4P2O7改性生物炭中C— P— O和C=O与(或) P=O的相对含量减少, C— O— P, O=P— O C— O— C与(或)P— O— C和P— O— P的相对含量增加。 改性生物炭中形成了活性较高的P结合方式。

(3)通过H3PO4和H4P2O7改性后, 麦秆与棉秆生物炭中活性P含量与中等活性P含量均显著增加。 H4P2O7改性生物炭中活性P与中等活性P含量显著高于H3PO4改性生物炭(p< 0.05)。 生物炭改性后P的活性增强。

参考文献
[1] Qian T T, Zhang X S, Hu J Y, et al. Chemosphere, 2013, 93(9): 3. [本文引用:2]
[2] Bi Z H, Huo L, Kong Q Q, et al. ACS Applied Materials & Interfaces, 2019, 3(24): 23. [本文引用:2]
[3] Valero-Romero M J, García-Mateos F J, Rodríguez-Mirasol J, et al. Fuel Processing Technology, 2017, 11(14): 119. [本文引用:2]
[4] Bulusheva L G, Arkhipov V E, Popov K M, et al. Materials, 2020, 13(5): 1173. [本文引用:1]
[5] Wang R, Ma Y Y, Zhang D, et al. E3S Web of Conferences, 2020, 145(2): 77. [本文引用:1]
[6] Hedley M J, Stewart J W B, ChauhanB S. Soil Science Society of America Journal, 1982, 46(5): 973. [本文引用:1]
[7] Wei J, Tu C, Yuan G D, et al. Environmental Pollution, 2019, 4(251): 59. [本文引用:1]
[8] Sun T, Xu Y M, Sun Y B, et al. Environmental Pollution, 2020, 8(271): 6. [本文引用:1]
[9] WENG Shi-fu, XU Yi-zhuang(翁诗甫, 徐怡庄). Fourier Transform Infrared Spectroscopy Analysis(傅里叶变换红外光谱分析). Third edition(3版). Beijing: Chemical Industry Press(北京: 化学工业出版社), 2016. 412. [本文引用:1]
[10] Liu Y T, Li K X, Pu L T, et al. Journal of Materials Chemistry A, 2015, 3(42): 21152. [本文引用:6]
[11] Wei J, Tu C, Yuan G D, et al. Environmental Pollution, 2019, 4(251): 58. [本文引用:1]
[12] Qian T T, Yang Q, Jun F C D, et al. Chemical Engineering Journal, 2018, 11(15): 1577. [本文引用:1]
[13] Liu Q, Li D, Chen H, et al. Bioresource Technology, 2021, 329: 8. [本文引用:1]
[14] Taskin M B, Kadioglu Y K, Sahin O, et al. Communications in Soil Science and Plant Analysis, 2019, 50(13): 5. [本文引用:1]
[15] Teng Z, Zhu J, Shao W, et al. Journal of Environmental Management, 2020, 256(15): 7. [本文引用:1]