环境科学  2021, Vol. 42 Issue (5): 2469-2479   PDF    
长期施用化肥和有机肥对稻田土壤重金属及其有效性的影响
夏文建, 张丽芳, 刘增兵, 张文学, 蓝贤瑾, 刘秀梅, 刘佳, 刘光荣, 李祖章, 王萍     
江西省农业科学院土壤肥料与资源环境研究所, 国家红壤改良工程技术研究中心, 国家农业环境宜春观测实验站, 南昌 330200
摘要: 为研究长期施肥下红壤稻田土壤重金属累积及有效态分布特征,探讨其主要影响因子,选择持续了35年的长期定位试验的不施肥(CK)、磷钾化肥(PK)、氮磷钾化肥(NPK)和不同比例有机肥化肥配施(M1NPK、M2NPK和M3NPK)等6个处理,分析了土壤化学指标、重金属全量(THM)和有效态的分布特征,以及重金属有效态含量与化学指标和全量之间的关系.结果表明,长期施肥改变了土壤化学性质,与CK处理相比,PK处理显著提高了土壤速效磷(AP)和速效钾(AK)含量,NPK处理显著提高了土壤有机质(SOM)、阳离子交换量(CEC)、AP和AK含量,有机肥处理(M1NPK、M2NPK和M3NPK)则大幅提高了SOM、CEC、AP、AK和硝态氮(NO3--N)含量;化肥处理(PK和NPK)土壤重金属全量总体上变化不大,有机肥处理土壤Cu、Zn和Cd的全量显著提高;化肥处理显著提高了土壤有效态Cr和As含量,有机肥处理显著提高了土壤有效态Cu、Zn、Cr、Cd、As和Fe的含量;土壤有效态Cu、Zn、Cr、Cd、As、Fe与SOM、CEC、AP、NO3--N含量之间呈显著正相关,另外有效态Zn和Cd还与土壤pH呈显著正相关,而有效态Pb与土壤pH、SOM、CEC和NO3--N含量之间呈显著负相关;土壤Cu、Zn和Cd有效态含量与全量之间显著正相关,而Cr、Fe有效态含量与全量之间显著负相关;冗余分析发现,SOM和pH分别解释了重金属全量方差变异的80.7%和5.5%,土壤CEC、AP和pH分别解释了重金属有效态方差变异的81.1%、4.9%和3.3%,达到显著水平;进一步通过偏最小二乘路径模型(PLS-PM)分析,发现THM、CEC和AP对土壤重金属有效态影响的路径系数分别为0.459、0.417和0.293,达到显著水平.长期施用猪粪等有机肥显著改善了土壤化学性质,并影响了土壤重金属有效性,其中土壤CEC和AP可能起到了关键的调控作用.
关键词: 长期施肥      红壤稻田      重金属      有效性      PLS-PM分析     
Effects of Long-Term Application of Chemical Fertilizers and Organic Fertilizers on Heavy Metals and Their Availability in Reddish Paddy Soil
XIA Wen-jian , ZHANG Li-fang , LIU Zeng-bing , ZHANG Wen-xue , LAN Xian-jin , LIU Xiu-mei , LIU Jia , LIU Guang-rong , LI Zu-zhang , WANG Ping     
Institute of Soil Fertilizer and Resource Environment, Jiangxi Academy of Agricultural Sciences, National Engineering and Technology Research Center for Red Soil Improvement, Yichun Scientific Observing and Experimental Station of Agro-environment, Nanchang 330200, China
Abstract: To investigate the effects of long-term fertilization on the accumulation and availability of heavy metals in reddish paddy soil and to analyze the major influencing factors, soil samples were collected after the later rice was harvested in 2018 from a long-term fertilization field experiment that began in 1984. Six treatments were selected, namely CK (control without fertilization), PK (P and K fertilizer), and NPK (N, P, and K fertilizer), and different proportions of organic fertilizer plus chemical fertilizer (M1NPK: 30%M+70%NPK; M2NPK: 50%M+50%NPK; and M3NPK: 70%M+30%NPK), soil chemical properties, total and available heavy metal contents, and the relationships between the available forms of heavy metals and soil chemical parameters and total heavy metals (THM) were analyzed. The results showed that ① long-term fertilization changed the soil chemical properties; compared with those of CK, PK significantly increased the contents of soil available phosphorus (AP) and available potassium (AK), NPK significantly increased the soil organic matter (SOM), cation exchange capacity (CEC), AP, and AK, and the organic fertilizer treatments significantly increased the contents of SOM, CEC, AP, AK, and nitrate (NO3--N). ② There were slight variations in the THM contents under the chemical fertilizer treatments (PK and NPK), whereas the organic fertilizer treatments significantly increased the total contents of Cu, Zn, and Cd. ③ The chemical fertilizer treatments significantly increased the available Cr and As, whereas the organic fertilizer treatments significantly increased the available Cu, Zn, Cr, Cd, As, and Fe. ④ There were significant positive correlations between the available Cu, Zn, Cr, Cd, As, and Fe and the SOM, CEC, AP, and NO3--N. In addition, the available Zn and Cd were significantly positively correlated with the soil pH, whereas the available Pb was significantly negatively correlated with soil pH, SOM, CEC, and NO3--N. ⑤ There were significant positive correlations between the available and total contents of Cu, Zn, and Cd, whereas there were significant negative correlations between the available and total contents of Cr and Fe. ⑥ Redundancy analysis showed that SOM and pH accounted for 80.7% and 5.5% of the variation in THM, whereas the soil CEC, AP, and pH accounted for 81.1%, 4.9%, and 3.3% of the variation in the available heavy metals, respectively. ⑦ The partial least squares path model analysis showed that the path coefficients of the THM, CEC, and AP on the available state of heavy metals were 0.459, 0.417, and 0.293, respectively. Long-term application of organic manure, such as pig manure, significantly improved the soil chemical properties and affected the availability of heavy metals, and soil CEC and AP may play key roles in regulation.
Key words: long-term fertilization      reddish paddy soil      heavy metals      availability      PLS-PM analysis     

施用化肥和有机肥是提高粮食产量的重要措施[1], 但由于肥料中含有少量重金属元素, 特别是有机肥中重金属经常超标[2, 3], 长期施用化肥或有机肥可能造成土壤重金属累积[2~6].我国南方地区畜禽粪便是土壤重金属的主要来源之一[7], 但施用有机肥同时也是改良土壤结构, 减少重金属活性的重要措施[8~10].有机肥可以提高土壤中有机物质含量, 通过有机物的表面羧基、羟基等官能团与重金属结合[11, 12], 或者通过提高土壤pH[8], 从而降低土壤重金属有效态.薛毅等[13]在湖南紫泥田的研究发现, 通过连续4a施用商品有机肥, 可以降低土壤Cd有效态从而降低稻米Cd含量; 刘灿等[6]在四川钙质水稻土中的研究发现, 有机肥及有机肥与化肥配施后Cu和Zn全量显著提高, 而有效态则是Cu和Cd显著提高, 有效态Zn差异不大; 何其辉等[14]通过盆栽试验研究发现, 施用猪粪使红黄泥Cu、Pb和Cd有效态含量显著降低了22.1%、13.4%和20.0%, 但对河沙泥Cu、Pb和Cd有效态影响不显著、Zn有效态显著上升, 而商品有机肥显著提高了红黄泥有效态Zn的含量, 同时显著提高了河沙泥Zn、Pb和Cd有效态含量; 王开峰等[15]分析湖南省的7个稻田长期定位试验发现, 中、高量有机肥明显提高了土壤Zn、Cu和Cd的全量及有效态含量.可见, 施用有机肥容易造成土壤重金属累积已达成共识, 但有机肥对重金属有效态的影响结论存在较大分歧, 这与有机肥的种类与用量、土壤环境、种植制度和重金属的类型等多方面的因素有关.同样, 长期化肥处理也可能导致重金属全量和活性增加.陈芳等[16]在河南潮土的长期定位试验研究发现, 土壤中As、Hg、Cd和Pb含量随耕作年限呈增长趋势, 磷肥可能是重金属增加的主要原因; 而刘灿等[6]的研究发现长期施用化肥土壤重金属全量变化不明显, 磷肥显著提高了有效Fe含量, 降低有效Mn含量.总体上, 现有研究主要围绕化肥或有机肥本身带入重金属的累积风险, 关于长期施肥造成土壤化学性质变化, 从而影响重金属有效态及相关机制研究较薄弱.

一般认为重金属全量对有效性具有较大影响[17], 但长期施肥条件下重金属有效态含量与全量的相关性因元素而异[6].常同举等[18]在紫色水稻土长期不同耕作方式的定位试验中发现土壤Zn有效态与全量呈显著正相关关系, 而Cu、Pb和Cd有效态与全量之间的相关性不显著; 李双异等[19]的研究发现Cu、Zn和Cd有效态含量与全量呈极显著正相关关系, Pb有效态与全量呈显著负相关; 刘灿等[6]在钙质水稻土长期施肥研究中发现Mn和Cu有效态与全量呈显著正相关关系, 而Fe、Zn、Pb和Cd有效态与全量之间的相关性不显著.长期施肥影响土壤理化性质, 从而影响重金属有效性, 土壤pH、有机质、CEC、有效氮和有效磷等均会影响重金属有效性[20], 但不同环境条件下影响重金属有效态的主要因子存在较大差异.长期施用化肥和有机肥对土壤重金属污染风险, 以及重金属有效性的关键影响因素尚不清楚.

为了研究长期施肥对红壤稻田重金属有效态的影响及土壤环境作用机制, 本文利用持续了35a的长期定位试验, 分析了Cu、Zn、Cr、Cd、Pb、As、Hg、Fe和Mn全量和有效态含量, 通过相关分析、主成分分析、冗余分析和偏最小二乘路径模型(PLS-PM)分析土壤重金属分布特征, 以及重金属有效态与环境因子之间的关系, 通过解析长期施肥土壤重金属有效态的关键影响因素, 以期为红壤稻田合理施肥减少环境风险提供理论参考.

1 材料与方法 1.1 研究区概况

长期定位试验于1984年开始[21], 地点位于江西省农业科学院试验基地(N28.57°, E115.94°), 该区域处于中亚热带, 海拔高度25 m, 年均温17.5℃, 无霜期280 d, 降雨量约1 600 mm, 土壤为第四纪酸性红壤发育的中潴黄泥田.试验前土壤耕层基础肥力状况为:土壤pH 6.5, 有机质25.6 g·kg-1, 全氮1.36 g·kg-1, 全磷0.49 g·kg-1, 速效氮81.6 mg·kg-1, 速效磷20.8 mg·kg-1, 速效钾35.0 mg·kg-1, 阳离子交换量(CEC)7.54 cmol·kg-1.当地土壤Cu、Zn、Cr、Pb、Cd、As、Hg、Fe和Mn含量[22]分别为47.0、115、72.0、43.0、0.171、11.0、0.30、29.9×103和475 mg·kg-1.

1.2 试验设计

本研究选取长期定位试验中的6个处理, 分别为:不施肥(CK)、施用磷钾化肥(PK)、施用氮磷钾化肥(NPK)、70%化肥+30%有机肥(M1NPK)、50%化肥+50%有机肥(M2NPK)、30%化肥+70%有机肥(M3NPK).试验常年种植模式为早稻-晚稻-冬闲, 小区面积33.3 m2, 水泥梗分隔, 独立排灌.早稻季施用纯N 150 kg·hm-2, 晚稻季施用纯N 180 kg·hm-2, 早稻和晚稻各施P2O560 kg·hm-2, K2O 150 kg·hm-2.各处理肥料用量详见表 1.肥料品种氮肥为尿素, 磷肥为过磷酸钙(含P2O5 12%), 钾肥为氯化钾(含K2O 60%).有机肥早稻为紫云英, 晚稻为猪粪.猪粪含水率85%, 养分含量N 0.45%、P2O5 0.19%和K2O 0.60%; 紫云英含水率88%, 养分含量N 0.30%、P2O50.08%和K2O 0.23%.猪粪中Cu、Zn、Cr、Pb、Cd、As、Hg、Fe和Mn含量分别为110、280、2.35、1.16、0.27、0.48、0.09、391和93.2 mg·kg-1, 紫云英中分别为1.3、3.8、0.26、0.12、0.04、0.03、0.001、16.2和21.4 mg·kg-1.尿素中Cu、Zn、Fe和Mn含量分别为0.13、2.31、10.2和0.21 mg·kg-1, 其他重金属未检出; 氯化钾中Cu、Zn、Fe和Mn含量分别为1.31、5.8、106和8.42 mg·kg-1, 其他重金属未检出; 过磷酸钙Cu、Zn、Cr、Pb、Cd、As、Hg、Fe和Mn含量分别为5.80、22.8、1.87、0.36、0.33、0.01、0.01、431和23.2mg·kg-1.其中磷肥和有机肥全作基肥, 氮肥50%作基肥, 25%作分蘖肥, 25%作幼穗分化肥.钾肥全作追肥, 50%作分蘖肥, 50%作幼穗分化肥.所有小区的播种、移栽、灌溉和打药等日常管理措施与当地习惯相同.

表 1 不同试验处理每年肥料总用量 Table 1 Annual total amount of fertilizers for different treatments

1.3 研究方法

2018年晚稻收获后, 采用土钻每小区S型采集0~20 cm耕层土壤样品.按四分法分出一半鲜样分析土壤硝态氮和铵态氮, 剩余一半样品摊匀风干, 磨细过筛分析土壤化学性质.

土壤化学性质采用常规试验方法[23], 具体为:土壤pH(水∶土2.5∶1)采用酸度计电位法; 土壤有机质(SOM)采用重铬酸钾外加热法测定; 速效磷(AP)用Olsen法; 速效钾(AK)用1 mol·L-1 NH4OAc浸提-火焰光度法; 硝态氮(NO3--N)和铵态氮(NH4+-N)采用2 mol·L-1的KCl溶液提取, 流动分析仪测定; 阳离子交换量(CEC)采用乙酸铵(NH4OAc)交换法测定; 土壤重金属全量(THM)采用HCl-HNO3-HClO4消煮, 重金属有效态采用稀酸法提取[17, 24, 25], 消煮液和提取液中Cu、Zn、Cr、Pb、Cd、Fe和Mn含量采用原子吸收分光光度计测定, As和Hg含量采用原子荧光光度计测定.土壤重金属全量和有效态含量分析过程中分别采用国家标准样GSS-22及ASA-5a进行质量控制.

1.4 数据处理

本试验数据利用Excel 2016进行整理、统计分析、计算和绘图, 运用SPSS 17.0进行单因素方差分析和差异显著性检验(LSD法, P < 0.05); 采用R(4.0.0)中的“corrplot”包进行相关性分析, “pheatmap”包进行热图聚类分析, “plspm”包进行偏最小二乘路径回归分析(PLS-PM).重金属的全量和有效态分类根据聚类树进行划分, 聚类之前数据进行归一化处理, 距离算法采用欧几里德距离, 聚类采用最大距离法, 通过层次聚类获得聚类树.采用Canoco 5.0软件进行主成分分析(PCA)和冗余分析(RDA).

2 结果与分析 2.1 长期施肥对土壤化学性质的影响

经过35 a的长期不同施肥处理, 土壤化学性质发生了显著变化(表 2).与CK处理相比, PK处理显著提高了土壤AP和AK含量, 分别增加了747%和90.5%; NPK处理显著提高了SOM、CEC、AP和AK含量, 分别增加了18.8%、13.3%、343%和17.7%; 有机肥处理(M1NPK、M2NPK和M3NPK)显著提高了SOM、CEC、AP、AK和NO3--N含量, 分别增加了33.1%~56.4%、31.6%~69.3%、860%~1 057%、16.9%~38.4%和341%~958%, 并且随着有机肥用量增加呈增加趋势.与NPK处理相比, 有机肥处理(M1NPK、M2NPK和M3NPK)显著提高了SOM、CEC、AP和NO3--N含量, 高量有机肥(M3NPK)还显著提高了土壤pH值和AK含量.土壤NH4+-N含量范围为4.80~8.26 mg·kg-1, 不同处理之间差异不显著.

表 2 长期施肥对土壤化学性质的影响1) Table 2 Soil chemical properties under long-term fertilizations

2.2 长期施肥对土壤重金属全量和有效态的影响

长期施肥影响了土壤重金属全量和有效态, 不同元素之间存在较大差异(图 1).与CK相比, 化肥处理(PK和NPK)土壤重金属全量变化不大, 有机肥处理(M1NPK、M2NPK、M3NPK)显著提高了土壤Cu、Zn和Cd的含量, 但对As、Hg、Mn和Pb无显著影响, 而Cr和Fe有所下降.根据我国农田土壤标准GB15618-2018, 有机肥处理(M1NPK、M2NPK和M3NPK)土壤铜(56.5~82.0 mg·kg-1)轻微超标, 高量有机肥M3NPK处理土壤镉(0.486 mg·kg-1)轻微超标, 其他元素均未超标.Heatmap聚类结果显示, 不同施肥处理可以分为3类:①CK和PK; ②NPK; ③M1NPK、M2NPK和M3NPK, 不同重金属全量可以划分为3类:①Cu、Zn、Cd和Hg; ②As和Mn; ③Pb、Cr和Fe [图 1(a)].

T.Cu、T.Zn、T.Cr、T.Pb、T.Cd、T.As、T.Hg、T.Fe和T.Mn分别表示Cu、Zn、Cr、Pb、Cd、As、Hg、Fe和Mn全量; A.Cu、A.Zn、A.Cr、A.Pb、A.Cd、A.As、A.Hg、A.Fe和A.Mn分别表示Cu、Zn、Cr、Pb、Cd、As、Hg、Fe和Mn有效态含量; 图中的数值代表重金属全量或有效态含量, 其中Fe全量单位为g·kg-1, Hg有效态单位为μg·kg-1, 其他为mg·kg-1; 颜色表示同一列重金属全量或有效态不同处理之间的相对值, 同一列不同字母表示处理间差异显著(P < 0.05) 图 1 长期施肥土壤重金属全量和有效态的heatmap聚类分析 Fig. 1 Cluster heatmap of the total and available contents of heavy metals in reddish paddy soil under long-term fertilization

与CK相比, 化肥处理(PK和NPK)土壤有效态Cr和As含量显著上升, 其他重金属有效态含量呈略微增加趋势; 有机肥处理(M1NPK、M2NPK和M3NPK)显著提高了土壤有效态Cu、Zn、Cr、Cd、As和Fe的含量.聚类结果显示, 不同施肥处理可以分为4类:①M3NPK; ②M1NPK和M2NPK; ③CK和NPK; ④PK, 不同重金属元素有效态可以划分为3类:①Cu、Zn、Cd、Cr、Fe和As; ②Mn; ③Pb和Hg[图 1(b)].

与CK相比, 施用化肥土壤重金属有效态占全量的比例变化较小, PK处理显著提高了有效态Cu、Cr、As和Fe的比例, NPK处理显著提高了有效态Cr、As和Hg的比例; 而有机肥处理(M1NPK、M2NPK和M3NPK)显著提高了Cu、Zn、Cr、As和Fe有效态的比例, M3NPK处理显著降低了有效态Pb的比例, M1NPK处理显著提高了有效态Hg的比例(表 3).

表 3 长期施肥下土壤重金属有效态占全量的比例1)/% Table 3 Proportion of soil available heavy metals to total contents under long-term fertilization/%

2.3 土壤重金属有效性与土壤化学性质和重金属全量的相关性

土壤有效态Cu、Zn、Cr、Cd、As、Fe与SOM、CEC、AP、NO3--N含量之间呈显著正相关, 另外有效态Zn、Cd还与土壤pH呈显著正相关, 而有效态Pb与土壤pH、SOM、CEC和NO3--N含量之间呈显著负相关.有效态Hg和Mn与土壤化学指标之间相关性不显著, 土壤AK和NH4+-N与各重金属有效态之间相关性不显著[图 2(a)].

橙色椭圆表示重金属有效态与化学指标或全量之间呈显著正相关(P < 0.05), 蓝色椭圆表示呈显著负相关(P < 0.05), 椭圆面积越小表示相关系数绝对值越大 图 2 重金属有效量与土壤化学性质和重金属全量之间的相关关系 Fig. 2 Correlations between available heavy metals and soil chemical properties and total heavy metals

土壤Cu、Zn和Cd有效态含量与全量之间显著正相关, 而Cr和Fe有效态含量与全量之间显著负相关, 其余重金属有效态含量与全量之间相关性不显著[图 2(b)].同时有效态Cu、Zn、Cd、As、Fe与全量Cu、Zn、Cd、Hg显著正相关, 有效态Pb与全量Cu、Zn、Cd和As显著负相关, 全量Cr与有效态As、Hg和Fe显著负相关, 全量Fe与有效态Cu、Cr、As和Hg显著负相关.可见, 受长期施肥影响土壤中不同重金属有效态与全量之间的关系存在差异, 并且不同类型重金属有效态与其他元素之间存在较为复杂的相关关系, 其原因还有待深入分析.

2.4 土壤重金属的主成分分析与冗余分析

通过对9种土壤重金属的全量和有效态分别进行主成分分析, 重金属全量和有效态的前两个主成分共解释了总体方差变异的93.6%和94.4%(图 3).长期不同施肥使土壤中重金属全量和有效态发生了明显变化, 从主成分分布图中可以看出, 不同施肥处理重金属的全量和有效态分布规律一致, CK处理分布在第二象限, 化肥处理(PK和NPK)分布在第二和第三象限, 而有机肥处理(M1NPK、M2NPK和M3NPK)分布在第一和第四象限, 并且随着有机肥用量增加, 样点分布由第四象限的左下方向第一象限右上方移动.

图 3 长期施肥对土壤重金属全量和有效态影响的主成分分析 Fig. 3 Principal component analysis (PCA) of total heavy metals and available heavy metals in soils under long-term fertilization

通过冗余分析发现, 土壤重金属全量和有效态与土壤化学性质间存在较强的约束关系.土壤重金属全量变异的88.1%可以由土壤pH、SOM、CEC、AP、NO3--N和NH4+-N解释, 其中SOM和pH分别解释了方差变异的80.7%和5.5%, 达到显著水平[图 4(a)]; 土壤重金属有效态变异的90.5%可以由土壤pH、SOM、CEC、AP、NO3--N和NH4+-N解释, 其中土壤CEC、AP和pH分别解释了方差变异的81.1%、4.9%和3.3%, 达到显著水平[图 4(b)].冗余分析结果, 反映了土壤重金属全量及有效态与土壤pH、SOM、CEC和AP之间有较强的关联.

图 4 土壤重金属全量和有效态与土壤化学指标的冗余分析 Fig. 4 Redundancy analysis (RDA) of the correlations between soil variables with total heavy metals and available heavy metals

2.5 土壤重金属有效性的PLS-PM分析

在前述分析的基础上, 进一步通过偏最小二乘路径模型(PLS-PM)分析了有机肥、化肥以及土壤化学性质(pH、SOM、CEC和AP)与重金属全量(THM)及重金属有效态(AHM)之间的关系, 模型的拟合优度为0.650(图 5).有机肥对SOM、CEC和AP和THM的路径系数分别为0.953、0.969、0.764和0.978, 达到显著水平, 但对土壤pH的影响不显著.化肥对土壤pH、SOM和CEC的路径系数分别为-0.436、0.277和0.128, 达到显著水平, 而对AP和THM之间路径系数不显著.CEC、AP和THM对有效态重金属影响的路径系数分别为0.417、0.293和0.459, 达到显著水平.土壤pH和SOM对重金属有效态的路径系数为负, 但不显著.有机肥和化肥主要通过影响土壤CEC和AP和THM间接影响土壤重金属的有效性, 间接效应系数分别为0.958和0.138(表 4), 可见有机肥的影响远大于化肥.

OM表示有机肥(organic manure), CF表示化肥(chemical fertilizers), THM表示重金属全量; 蓝色箭头表示路径系数为正, 红色箭头表示路径系数为负, 虚线表示路径系数不显著, 线条粗细表示路径系数大小 图 5 长期施肥对土壤重金属有效性影响的偏最小二乘路径模型(PLS-PM)分析 Fig. 5 Directed graph of the partial least squares path model (PLS-PM) of the effects of soil available heavy metals under long-term fertilization

表 4 长期施肥通过不同参数影响土壤重金属有效态的间接效应系数 Table 4 Indirect effect coefficient of long-term fertilization on soil heavy metal availability through different parameters

3 讨论 3.1 长期施肥对红壤稻田土壤重金属累积的影响

施肥是影响土壤重金属含量变化的重要因素, 一般长期施用化肥对土壤重金属全量的影响不大[2], 但仍然有累积的风险.任顺荣等[26]的研究发现25 a长期不同施肥土壤重金属Cu、Zn、Cd、Cr、Hg和As均呈增加趋势, 氮肥对重金属影响不大, 磷肥使重金属含量有所增加, 而有机肥使Cu、Zn和Hg增幅较大.由于氮肥和钾肥中重金属元素含量较低[27], 而磷肥中重金属Cd含量相对较高[28], 长期施用磷肥的情况下, 可能会造成Cd的累积[29].我国磷矿具有镉含量低的优势, 但酸性稻田仍然需要注意镉污染风险[25~30].本研究长期不施肥的对照处理水稻产量低, 对土壤重金属的影响较小, 与当地土壤的背景值相比变化不大.试验中化肥重金属含量低, 长期施用化肥处理(PK和NPK)土壤重金属全量与对照相比变化不大, NPK处理土壤Cd和As有下降趋势, 这与作物收获带走有关[5, 6].

长期施用有机肥或有机肥与化肥配施对土壤重金属全量和有效态有较大影响, 但不同重金属元素之间差异较大[2, 31], 这与有机肥的种类、来源、用量, 以及种植制度、土壤性质和元素本身的迁移性有关[31~35].王改玲等[36]的研究发现化肥配施商品有机肥增加了土壤Cu、Zn、As和Pb的含量, 降低土壤Hg含量, 但对土壤Cd、Cr和Ni的影响不显著; 猪粪连续施用显著增加耕层土壤Cu和Zn, 而施用鸡粪显著增加土壤Cu、Zn、Cd、Cr、As和Pb[32].本研究中有机肥为猪粪, 有机肥处理(M1NPK、M2NPK和M3NPK)显著提高了土壤Cu、Zn和Cd的含量, 其他重金属差异不显著.而有研究结果显示施用有机肥显著提高了耕层土壤Cu、Zn和As含量[37], 这可能与早期猪粪中As含量较高有关.本研究中有机肥处理As的全量无明显增加, 但有效态含量显著提高, 其原因还有待进一步分析.刘荣乐等[38]的研究发现我国主要商品有机肥中重金属含量变异很大, 商品有机肥中重金属含量与所使用的有机物料中重金属含量呈显著正相关, 其中猪粪中重金属以Cu、Zn和Cd为主, 而江西省猪粪中Cu和Zn含量超过全国平均水平[2], 可能是本研究中有机肥处理土壤Cu和Zn含量显著提高的主要原因.经过35 a的长期施肥, 与我国农田土壤标准GB 15618-2018相比, 有机肥处理(M1NPK、M2NPK、M3NPK)土壤铜(56.5~82.0 mg·kg-1)轻微超标, 高量有机肥M3NPK处理土壤镉(0.486 mg·kg-1)轻微超标.同时施用有机肥处理籽粒中Cu、Zn和Cd含量有明显的增加, 其中M3NPK处理水稻籽粒Cd含量达到0.28 mg·kg-1, 超过国家食品安全限值, 与李本银等[39]的研究结果一致.因此有机肥长期施用需要考虑生态效益和环境承载力, 避免过量施用造成环境和粮食安全风险[39, 40].

本文通过对土壤重金属进行Heatmap聚类分析发现, 重金属全量不同施肥处理可以分为: ① CK和PK; ② NPK; ③ M1NPK、M2NPK和M3NPK这3类[图 1(a)], 重金属有效态不同施肥处理可以分为: ① M3NPK; ② M1NPK和M2NPK; ③ CK和NPK; ④ PK这4类[图 1(b)], 长期施用有机肥处理土壤重金属全量和有效态发生了显著变化, 该结果与主成分分析结果(图 3)一致.根据聚类树重金属全量可以划分为: ① Cu、Zn、Cd和Hg; ② As和Mn; ③Pb、Cr和Fe这3类, 施用有机肥土壤Cu、Zn和Cd含量显著增加, 而土壤Hg含量变化规律与Cu、Zn、Cd一致[图 1(a)], 可能是这4个元素聚到一起的原因.有机肥处理(M1NPK、M2NPK、M3NPK)显著提高了土壤有效态Cu、Zn、Cr、Cd、As和Fe的含量, 不同处理有效态Mn的含量差异均不显著, 而有效态Hg和Pb含量随着有机肥用量增加呈下降趋势, 因此重金属有效态可以划分为: ① Cu、Zn、Cd、Cr、Fe和As; ② Mn; ③ Pb和Hg这3类.

3.2 长期施肥对红壤稻田土壤重金属有效性的影响

施用有机肥不仅带入了大量重金属, 使土壤重金属全量分布情况有明显变化[图 1(a)], 同时重金属有效态的分布也与对照和化肥处理明显不同[图 1(b)].然而不同重金属元素有效态与全量之间的相关性并不一致[图 2(b)], 这可能与不同元素本身的性质和来源有关[6, 41], 在稻田[18]、旱地[42]和菜地[43]的研究中也发现部分重金属有效态与对应重金属全量间的相关性不高.本研究中Cu、Zn和Cd有效态与全量具有较高的相关性, 表明Cu、Zn和Cd的有效态和全量的来源都主要为有机肥带入.施用有机肥处理(M1NPK、M2NPK和M3NPK)明显提高了Cu、Zn、Cr、As和Fe有效态与全量之间的比值, 但Cd有效态与全量比值没有明显变化(表 3).长期施肥影响土壤重金属有效性, 一方面与施肥带入有关, 另一方面由于施肥能显著改变土壤化学性质, 如土壤pH、土壤溶液中离子组成、阳离子交换量等, 或肥料成分直接与土壤中金属离子发生反应, 从而影响重金属的有效性[39~43].

相关分析结果表明土壤重金属有效态主要受SOM、CEC、AP、NO3--N和pH影响, 其中Zn、Cd有效态与上述指标显著正相关, 而Pb有效态与上述指标显著负相关[图 2(a)].进一步通过RDA分析发现重金属全量主要受SOM和pH影响[图 4(a)], 而重金属有效态主要受CEC、AP和pH影响[图 4(b)].在此基础上通过PLS-PM分析发现, 长期施用化肥对重金属的间接效应系数为0.138, 主要通过AP(0.074)、CEC(0.054)和pH(0.031)影响重金属有效性; 而有机肥的间接效应系数为0.959, 主要通过重金属全量THM(0.449)、CEC(0.404)和AP(0.224)影响重金属有效性(表 4).张民等[44]的研究发现土壤重金属的全量和有效态与有机质和速效磷含量呈显著正相关, 猪粪中全磷和有机质含量较高[45], 长期施用猪粪使土壤Cu、Zn和Cd发生累积, 同时也通过影响土壤SOM、pH、CEC和AP等影响重金属有效性.稻田土壤中磷对重金属的固定和活化受到土壤pH、氧化还原反应、根际酸化效应和根表铁膜等影响[46], 稻田渗滤液中磷和重金属含量显著相关[47], 但磷与重金属之间存在复杂的源汇关系和固定活化作用, 两者之间的作用和相关机制还有待深入研究[46].通过改善土壤理化性质有助于缓解土壤重金属污染[48], 酸性土壤中碱性调节效果优于有机调节[34], 而本研究发现土壤CEC和AP可能是酸性稻田中重金属有效态调控的关键因子, 但具体的调控措施还有待进一步研究.

4 结论

(1) 长期施肥改变了土壤化学性质, 特别是有机肥处理大幅提高了土壤SOM、CEC、AP、AK和NO3--N.长期施用化肥对土壤重金属全量影响不大, 长期施用有机肥使土壤Cu、Zn和Cd全量和有效态显著提高, 并显著提高了Cr、As和Fe有效态含量.

(2) Heatmap聚类结果显示, 重金属全量可以划分为①CK和PK; ②NPK; ③M1NPK、M2NPK和M3NPK这3类, 重金属有效态可以划分为①Cu、Zn、Cd、Cr、Fe和As; ②Mn; ③Pb和Hg这3类.

(3) 土壤CEC、AP和pH分别解释了重金属有效态方差变异的81.1%、4.9 %和3.3%; PLS-PM分析发现, CEC和AP对重金属有效态影响的路径系数分别为0.417和0.293, 达显著水平.长期施肥条件下土壤CEC和AP可能对重金属有效态起关键调控作用.

参考文献
[1] 朱兆良, 金继运. 保障我国粮食安全的肥料问题[J]. 植物营养与肥料学报, 2013, 19(2): 259-273.
Zhu Z L, Jin J Y. Fertilizer use and food security in China[J]. Plant Nutrition and Fertilizer Science, 2013, 19(2): 259-273.
[2] 穆虹宇, 庄重, 李彦明, 等. 我国畜禽粪便重金属含量特征及土壤累积风险分析[J]. 环境科学, 2020, 41(2): 986-996.
Mu H Y, Zhuang Z, Li Y M, et al. Heavy metal contents in animal manure in China and the related soil accumulation risks[J]. Environmental Science, 2020, 41(2): 986-996.
[3] 王美, 李书田. 肥料重金属含量状况及施肥对土壤和作物重金属富集的影响[J]. 植物营养与肥料学报, 2014, 20(2): 466-480.
Wang M, Li S T. Heavy metals in fertilizers and effect of the fertilization on heavy metal accumulation in soils and crops[J]. Journal of Plant Nutrition and Fertilizer, 2014, 20(2): 466-480.
[4] Uprety D, Hejcman M, Száková J, et al. Concentration of trace elements in arable soil after long-term application of organic and inorganic fertilizers[J]. Nutrient Cycling in Agroecosystems, 2009, 85(3): 241-252. DOI:10.1007/s10705-009-9263-x
[5] Ning C C, Gao P D, Wang B Q, et al. Impacts of chemical fertilizer reduction and organic amendments supplementation on soil nutrient, enzyme activity and heavy metal content[J]. Journal of Integrative Agriculture, 2017, 16(8): 1819-1831. DOI:10.1016/S2095-3119(16)61476-4
[6] 刘灿, 秦鱼生, 赵秀兰. 长期不同施肥对钙质紫色水稻土重金属累积及有效性的影响[J]. 农业环境科学学报, 2020, 39(7): 1494-1502.
Liu C, Qin Y S, Zhao X L. Long-term effect of fertilization on accumulation and availability of heavy metal in a calcareous paddy soil[J]. Journal of Agro-Environment Science, 2020, 39(7): 1494-1502.
[7] Peng H, Chen Y L, Weng L P, et al. Comparisons of heavy metal input inventory in agricultural soils in North and South China: a review[J]. Science of the Total Environment, 2019, 660: 776-786. DOI:10.1016/j.scitotenv.2019.01.066
[8] Huang Q Q, Yu Y, Wan Y N, et al. Effects of continuous fertilization on bioavailability and fractionation of cadmium in soil and its uptake by rice (Oryza sativa L.)[J]. Journal of Environmental Management, 2018, 215: 13-21.
[9] Zhen H Y, Jia L, Huang C D, et al. Long-term effects of intensive application of manure on heavy metal pollution risk in protected-field vegetable production[J]. Environmental Pollution, 2020, 263. DOI:10.1016/j.envpol.2020.114552
[10] 宁川川, 王建武, 蔡昆争. 有机肥对土壤肥力和土壤环境质量的影响研究进展[J]. 生态环境学报, 2016, 25(1): 175-181.
Ning C C, Wang J W, Cai K Z. The effects of organic fertilizers on soil fertility and soil environmental quality: a review[J]. Ecology and Environmental Sciences, 2016, 25(1): 175-181.
[11] Li F, Li Z A, Mao P, et al. Heavy metal availability, bioaccessibility, and leachability in contaminated soil: effects of pig manure and earthworms[J]. Environmental Science and Pollution Research, 2019, 26(20): 20030-20039. DOI:10.1007/s11356-018-2080-5
[12] Wan Y N, Huang Q Q, Wang Q, et al. Accumulation and bioavailability of heavy metals in an acid soil and their uptake by paddy rice under continuous application of chicken and swine manure[J]. Journal of Hazardous Materials, 2020, 384. DOI:10.1016/j.jhazmat.2019.121293
[13] 薛毅, 尹泽润, 盛浩, 等. 连续4a施有机肥降低紫泥田镉活性与稻米镉含量[J]. 环境科学, 2020, 41(4): 1880-1887.
Xue Y, Yin Z R, Sheng H, et al. Reduction of soil cadmium activity and rice cadmium content by 4-year-consecutive application of organic fertilizer[J]. Environmental Science, 2020, 41(4): 1880-1887.
[14] 何其辉, 谭长银, 曹雪莹, 等. 肥料对土壤重金属有效态及水稻幼苗重金属积累的影响[J]. 环境科学研究, 2018, 31(5): 942-951.
He Q H, Tan C Y, Cao X Y, et al. Effects of fertilizer on the availability of heavy metals in soil and its accumulation in rice seedling[J]. Research of Environmental Sciences, 2018, 31(5): 942-951.
[15] 王开峰, 彭娜, 王凯荣, 等. 长期施用有机肥对稻田土壤重金属含量及其有效性的影响[J]. 水土保持学报, 2008, 22(1): 105-108.
Wang K F, Peng N, Wang K R, et al. Effects of long-term manure fertilization on heavy metal content and its availability in paddy soil[J]. Journal of Soil and Water Conservation, 2008, 22(1): 105-108. DOI:10.3321/j.issn:1009-2242.2008.01.023
[16] 陈芳, 董元华, 安琼, 等. 长期肥料定位试验条件下土壤中重金属的含量变化[J]. 土壤, 2005, 37(3): 308-311.
Chen F, Dong YH, An Q, et al. Variation of soil heavy metal contents in a long-term fertilization experiment[J]. Soils, 2005, 37(3): 308-311. DOI:10.3321/j.issn:0253-9829.2005.03.015
[17] 钟晓兰, 周生路, 李江涛, 等. 土壤有效态Cd、Cu、Pb的分布特征及影响因素研究[J]. 地理科学, 2010, 30(2): 254-260.
Zhong X L, Zhou S L, Li J T, et al. Soil available Cd, Cu, Pb distribution characteristic and its influencing factors[J]. Scientia Geographica Sinica, 2010, 30(2): 254-260.
[18] 常同举, 崔孝强, 阮震, 等. 长期不同耕作方式对紫色水稻土重金属含量及有效性的影响[J]. 环境科学, 2014, 35(6): 2381-2391.
Chang T J, Cui X Q, Ruan Z, et al. Long-term effects of tillage methods on heavy metal accumulation and availability in purple paddy soil[J]. Environmental Science, 2014, 35(6): 2381-2391.
[19] 李双异, 刘赫, 汪景宽. 长期定位施肥对棕壤重金属全量及其有效性影响[J]. 农业环境科学学报, 2010, 29(6): 1125-1129.
Li S Y, Liu H, Wang J K. Effects of long-term located fertilization on heavy metals and their availability in brown earth[J]. Journal of Agro-Environment Science, 2010, 29(6): 1125-1129.
[20] 王美, 李书田, 马义兵, 等. 长期不同施肥措施对土壤铜、锌、镉形态及生物有效性的影响[J]. 农业环境科学学报, 2014, 33(8): 1500-1510.
Wang M, Li S T, Ma Y B, et al. Influence of different long-term fertilization practices on fractionations and bioavailability of Cu, Zn, and Cd in soils[J]. Journal of Agro-Environment Science, 2014, 33(8): 1500-1510.
[21] 冀建华, 刘秀梅, 李祖章, 等. 长期施肥对黄泥田碳和氮及氮素利用的影响[J]. 中国农业科学, 2011, 44(12): 2484-2494.
Ji J H, Liu X M, Li Z Z, et al. Effects of long-term fertilization on carbon and nitrogen in yellow clayey soil and its nitrogen utilization[J]. Scientia Agricultura Sinica, 2011, 44(12): 2484-2494. DOI:10.3864/j.issn.0578-1752.2011.12.010
[22] 毛大发, 鄢新华, 刘小兵, 等. 试论南昌—莲塘一带土壤环境地球化学特征及其环境质量[J]. 地质与勘探, 2003, 39(3): 72-77.
Mao D F, Yan X H, Liu X B, et al. Soil geochemical characteristic and environment quality in Nanchang—Liantang region[J]. Geology and Prospecting, 2003, 39(3): 72-77. DOI:10.3969/j.issn.0495-5331.2003.03.016
[23] 鲍士旦. 土壤农化分析[M]. (第三版). 北京: 中国农业出版社, 2000.
[24] 李发生, 韩梅, 熊代群, 等. 不同浸提剂对几种典型土壤中重金属有效态的浸提效率研究[J]. 农业环境科学学报, 2003, 22(6): 704-706.
Li F S, Han M, Xiong D Q, et al. Efficiency of some extractants for available heavy metals from several typical soils[J]. Journal of Agro-Environment Science, 2003, 22(6): 704-706. DOI:10.3321/j.issn:1672-2043.2003.06.016
[25] 周卫红, 张静静, 邹萌萌, 等. 土壤重金属有效态含量检测与监测现状、问题及展望[J]. 中国生态农业学报, 2017, 25(4): 605-615.
Zhou W H, Zhang J J, Zou M M, et al. The detection and monitoring of available heavy metal content in soil: a review[J]. Chinese Journal of Eco-Agriculture, 2017, 25(4): 605-615.
[26] 任顺荣, 邵玉翠, 高宝岩, 等. 长期定位施肥对土壤重金属含量的影响[J]. 水土保持学报, 2005, 19(4): 96-99.
Ren S R, Shao Y C, Gao B Y, et al. Effects of long-term located fertilization on heavy-metal content of soil[J]. Journal of Soil and Water Conservation, 2005, 19(4): 96-99. DOI:10.3321/j.issn:1009-2242.2005.04.023
[27] Raven K P, Loeppert R H. Trace element composition of fertilizers and soil amendments[J]. Journal of Environmental Quality, 1997, 26(2): 551-557.
[28] 王起超, 麻壮伟. 某些市售化肥的重金属含量水平及环境风险[J]. 农村生态环境, 2004, 20(2): 62-64.
Wang Q C, Ma Z W. Heavy metals in chemical fertilizer and environmental risks[J]. Rural Eco-Environment, 2004, 20(2): 62-64.
[29] 杨旭, 余垚, 李花粉, 等. 我国与欧美化肥重金属限量标准的比较和启示[J]. 植物营养与肥料学报, 2019, 25(1): 149-156.
Yang X, Yu Y, Li H F, et al. Comparison of heavy metal limits for chemical fertilizers in China, EU and US and enlightenments[J]. Journal of Plant Nutrition and Fertilizers, 2019, 25(1): 149-156.
[30] 马文奇, 张福锁, 张卫锋. 关乎我国资源、环境、粮食安全和可持续发展的化肥产业[J]. 资源科学, 2005, 27(3): 33-40.
Ma W Q, Zhang F S, Zhang W F. Fertilizer production and consumption and the resources, environment, food security and sustainable development in China[J]. Resources Science, 2005, 27(3): 33-40. DOI:10.3321/j.issn:1007-7588.2005.03.006
[31] Sukitprapanon T, Suddhiprakarn A, Kheoruenromne I, et al. Partitioning and potential mobilization of aluminum, arsenic, iron, and heavy metals in tropical active and post-active acid sulfate soils: Influence of long-term paddy rice cultivation[J]. Chemosphere, 2018, 197: 691-702. DOI:10.1016/j.chemosphere.2018.01.099
[32] 何梦媛, 董同喜, 茹淑华, 等. 畜禽粪便有机肥中重金属在土壤剖面中积累迁移特征及生物有效性差异[J]. 环境科学, 2017, 38(4): 1576-1586.
He M Y, Dong T X, Ru S H, et al. Accumulation and migration characteristics in soil profiles and bioavailability of heavy metals from livestock manure[J]. Environmental Science, 2017, 38(4): 1576-1586.
[33] 陈卫平, 杨阳, 谢天, 等. 中国农田土壤重金属污染防治挑战与对策[J]. 土壤学报, 2018, 55(2): 261-272.
Chen W P, Yang Y, Xie T, et al. Challenges and countermeasures for heavy metal pollution control in farmlands of China[J]. Acta Pedologica Sinica, 2018, 55(2): 261-272.
[34] 王腾飞, 谭长银, 曹雪莹, 等. 长期施肥对土壤重金属积累和有效性的影响[J]. 农业环境科学学报, 2017, 36(2): 257-263.
Wang T F, Tan C Y, Cao X Y, et al. Effects of long-term fertilization on the accumulation and availability of heavy metals in soil[J]. Journal of Agro-Environment Science, 2017, 36(2): 257-263.
[35] Lu H L, Wu Y X, Liang P X, et al. Alkaline amendments improve the health of soils degraded by metal contamination and acidification: crop performance and soil bacterial community responses[J]. Chemosphere, 2020, 257. DOI:10.1016/j.chemosphere.2020.127309
[36] 王改玲, 李立科, 郝明德, 等. 长期定位施肥对土壤重金属含量的影响及环境评价[J]. 水土保持学报, 2010, 24(3): 60-63, 70.
Wang G L, Li L K, Hao M D, et al. Effects of long-term fertilization on heavy-metal contents of soil and environmental quality evaluation[J]. Journal of Soil and Water Conservation, 2010, 24(3): 60-63, 70.
[37] 李祖章, 谢金防, 蔡华东, 等. 农田土壤承载畜禽粪便能力研究[J]. 江西农业学报, 2010, 22(8): 140-145, 149.
Li Z Z, Xie J F, Cai H D, et al. Environmental loading capacity of farmland soil for dung of livestock and poultry[J]. Acta Agriculturae Jiangxi, 2010, 22(8): 140-145, 149. DOI:10.3969/j.issn.1001-8581.2010.08.043
[38] 刘荣乐, 李书田, 王秀斌, 等. 我国商品有机肥料和有机废弃物中重金属的含量状况与分析[J]. 农业环境科学学报, 2005, 24(2): 392-397.
Liu R L, Li S T, Wang X B, et al. Contents of heavy metal in commercial organic fertilizers and organic wastes[J]. Journal of Agro-Environment Science, 2005, 24(2): 392-397. DOI:10.3321/j.issn:1672-2043.2005.02.041
[39] 李本银, 汪鹏, 吴晓晨, 等. 长期肥料试验对土壤和水稻微量元素及重金属含量的影响[J]. 土壤学报, 2009, 46(2): 281-288.
Li B Y, Wang P, Wu X C, et al. Effect of long-term fertilization experiment on concentration of micronutrients and heavy metals in soil and brown rice[J]. Acta Pedologica Sinica, 2009, 46(2): 281-288. DOI:10.3321/j.issn:0564-3929.2009.02.013
[40] 柳开楼, 李大明, 黄庆海, 等. 红壤稻田长期施用猪粪的生态效益及承载力评估[J]. 中国农业科学, 2014, 47(2): 303-313.
Liu K L, Li D M, Huang Q H, et al. Ecological benefits and environmental carrying capacities of red paddy field subjected to long-term pig manure amendments[J]. Scientia Agricultura Sinica, 2014, 47(2): 303-313.
[41] 张磊, 宋凤斌, 崔良. 化肥施用对土壤中重金属生物有效性的影响研究[J]. 中国生态农业学报, 2006, 14(4): 122-125.
Zhang L, Song F B, Cui L. Effects of application of chemical fertilizers on bioavailability of heavy metals in soil[J]. Chinese Journal of Eco-Agriculture, 2006, 14(4): 122-125.
[42] Zhou S W, Liu J, Xu M G, et al. Accumulation, availability, and uptake of heavy metals in a red soil after 22-year fertilization and cropping[J]. Environmental Science and Pollution Research, 2015, 22(19): 15154-15163. DOI:10.1007/s11356-015-4745-7
[43] Hao X Z, Zhou D M, Huang D Q, et al. Heavy metal transfer from soil to vegetable in southern Jiangsu Province, China[J]. Pedosphere, 2009, 19(3): 305-311. DOI:10.1016/S1002-0160(09)60121-1
[44] 张民, 龚子同. 我国菜园土壤中某些重金属元素的含量与分布[J]. 土壤学报, 1996, 33(1): 85-93.
Zhang M, Gong Z T. Contents and distribution of some heavy metal elements in the vegetable cultivated soils in China[J]. Acta Pedologica Sinica, 1996, 33(1): 85-93. DOI:10.3321/j.issn:0564-3929.1996.01.010
[45] 黄绍文, 唐继伟, 李春花. 我国商品有机肥和有机废弃物中重金属、养分和盐分状况[J]. 植物营养与肥料学报, 2017, 23(1): 162-173.
Huang S W, Tang J W, Li C H. Status of heavy metals, nutrients, and total salts in commercial organic fertilizers and organic wastes in China[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(1): 162-173.
[46] Seshadri B, Bolan N S, Wijesekara H, et al. Phosphorus-cadmium interactions in paddy soils[J]. Geoderma, 2016, 270: 43-59. DOI:10.1016/j.geoderma.2015.11.029
[47] Liu X P, Bi Q F, Qiu L L, et al. Increased risk of phosphorus and metal leaching from paddy soils after excessive manure application: insights from a mesocosm study[J]. Science of the Total Environment, 2019, 666: 778-785. DOI:10.1016/j.scitotenv.2019.02.072
[48] Li C, Sanchez G M, Wu Z F, et al. Spatiotemporal patterns and drivers of soil contamination with heavy metals during an intensive urbanization period (1989-2018) in southern China[J]. Environmental Pollution, 2020, 260. DOI:10.1016/j.envpol.2020.114075