安吉¹, 李婷^1,3, 张世熔², 郑雯文¹, 狄黎平³, 郭鹏¹
1. 四川农业大学资源学院, 成都 611130;
2. 四川农业大学环境学院, 成都 611130;
3. Center for Spatial Information Science and System, George Mason University, Fairfax, VA22030, USA
收稿日期: 2018-02-08; 修回日期: 2018-03-17; 录用日期: 2018-03-17
基金项目: 国家自然科学基金项目（No.41601311）；四川省教育厅重点项目（No.17ZA0308）
作者简介: 安吉(1995-), 女, E-mail:sc_anji@163.com
通讯作者（责任作者）: 李婷, E-mail:lting1804@gmail.com

摘要: 城市化导致城市面积扩张和土地利用方式改变，从而影响了土壤重金属空间分布.本文以成都-温江这一典型"中心城-郊区-卫星城"梯度带为研究区，探讨绿地表层土壤重金属Cu、Pb和Zn的空间分布特征及其影响因素，以期为城市土壤重金属污染治理提供帮助.结果表明，研究区内土壤Cu、Pb和Zn含量均值分别为47.77、39.98和89.65 mg·kg^-1，其中土壤Cu和Pb含量均值超过土壤背景值.3种重金属含量在"中心城-郊区-卫星城"梯度带呈现先降低后增加的趋势.方差分析表明，区内不同绿地类型土壤Cu和Zn含量差异显著（p < 0.05），Pb含量在各个绿地类型下无显著性差异（p>0.05）；回归分析显示3种重金属含量随距城市中心距离增加而极显著降低（p < 0.01）；缓冲区分析表明，土壤Cu、Pb和Zn含量受局部范围内的建筑密度和道路密度影响，且随局部范围内建筑密度和道路密度的增加而显著增加（p < 0.05）.逐步回归分析显示，土壤Cu含量主要受建筑密度影响，Pb含量主要受道路密度和绿地类型影响，Zn含量则受绿地类型、道路密度和建筑密度的共同影响.研究结果显示城市化造成土壤重金属在一定程度上累积，距城市中心距离、绿地类型、建筑密度和道路密度均会影响土壤重金属的空间分布.
关键词:城市化土壤重金属空间分布建筑密度道路密度
Spatial distribution of soil heavy metal in green land along a "center city-suburb-satellite city" gradient
AN Ji¹, LI Ting^1,3 , ZHANG Shirong², ZHENG Wenwen¹, DI Liping³, GUO Peng¹
1. College of Resources, Sichuan Agricultural University, Chengdu 611130;
2. College of Environment, Sichuan Agricultural University, Chengdu 611130;
3. Center for Spatial Information Science and System, George Mason University, Fairfax, VA 22030, USA
Received 8 February 2018; received in revised from 17 March 2018; accepted 17 March 2018
Supported by the National Natural Science Foundation of China (No.41601311) and the Key Projects of Education Department of Sichuan (No.17ZA0308)
Biography: AN Ji (1995—), female, E-mail: sc_anji@163.com
^*Corresponding author: LI Ting, E-mail:lting1804@gmail.com
Abstract: Urbanization results in rapid expansion of construction land area and changes of land use types over the past decades, which greatly influences the spatial distribution of soil heavy metals around the city area. In this paper, the spatial distributions of soil heavy metals including Cu, Pb and Zn were studied in a typical "center city-suburban-satellite city" gradient zone of western Chengdu City of Sichuan Province, China, aimed to provide useful information for effectively prevention of soil heavy metal pollution. The results showed that the average concentrations of Cu, Pb and Zn in soils were 47.77, 39.98 and 89.65 mg·kg^-1, respectively, and the values of Cu and Pb concentrations in soil were higher than the background values. The concentrations of all the three soil heavy metals decreased from center city to suburban, and then increased from suburban to satellite city. The concentrations of soil Cu and Zn showed significant differences among different green space types (p < 0.05), while the concentrations of soil Pb showed no significant difference (p>0.05). The results of regression analysis indicated that the concentrations of the three heavy metals in soils significantly declined with the increasing distance from the city center (p < 0.01). The buffer zone analyses illustrated that the concentrations of Cu, Pb and Zn in the soils were positively correlated with the road density and building density in local area (p < 0.05). Stepwise regression analyses exhibited that road density was the determinants governing soil Cu concentrations, road density and different green space types were the main controlling factors for soil Pb, and the concentrations of soil Zn were together controlled by the above three factors. In conclusion, urbanization has led to the accumulation of heavy metals, and the distance from urban centers, different green space types, construction density and road density determined the spatial distribution of heavy metals in the process of urbanization.
Key words: urbanizationsoil heavy metalspatial distributionconstruction densityroad density
1 引言(Introduction)城市化促使城市人口增多、城市面积扩张, 导致城市用地需求增加, 农业用地向城市用地转变(Bren et al., 2017).该种用地方式的转变必然会导致土壤重金属分布产生差异, 造成土壤重金属不同程度的累积, 并通过大气、水体、食物链直接或间接地影响人体健康, 同时对城市生态环境产生了巨大的影响(柳云龙等, 2015；Necat A?ca et al., 2015；Bai et al., 2017).因此, 了解城市绿地土壤重金属的分布特征及其影响因素, 可为城市土壤重金属污染治理, 土地资源的可持续利用提供科学依据(方淑波等, 2016).
人类活动、土地利用方式和距城市中心距离等是导致土壤重金属空间分布差异的主要因素(张连科等, 2016；郭伟等, 2013；Wu et al., 2008；杨蕊等, 2016).如Pouyat和McDonnell(1991)、郭广慧等(2015)的研究表明, 城市中心土壤重金属含量明显高于乡村地区.庄腾飞等(2012)研究表明上海市表层土壤中Cu、Zn、Pb在国家级开发园区和一些居民区、商业区富集, 并以该区域为中心向其他地区辐射状递减；刘亚纳等(2016)研究发现, 不同功能区土壤重金属含量具有明显差异, 工业区污染最为严重, 其次是城市主干道和商业区, 城乡结合处污染程度最轻.已有研究中以市中心为直线距离作为城乡梯度的度量方法得到了广泛运用(李风平等, 2008；柳云龙等, 2012), 然而随着城市化的不断发展, 在大城市地区的影响范围之内, 且与中心城具有一定距离, 衍生了承担一部分城市功能的新兴城镇, 通常被称为“卫星城”, 从而形成“中心城-郊区-卫星城”梯度带(MacGregor-Fors, 2011；李万峰, 2014).新兴梯度带的产生必然会导致区域内土壤重金属含量及分布异质性, 然而迄今对于该种梯度带下土壤重金属分布的异质性及其变化机制仍缺乏必要的了解.
本文以“中心城-郊区-卫星城”典型梯度带“成都-温江”为研究区域, 分析该梯度带上绿地土壤重金属Cu、Pb和Zn的空间分布特征, 并采用距城市中心距离、建筑密度和道路密度等城市化量化指标, 探究城市化对土壤重金属的影响；并通过分析居住绿地、公园绿地、交通绿地、工业绿地、附属绿地土壤中重金属的含量差异了解城市土壤重金属的污染状况, 为城市土地资源的持续利用和污染治理提供科学依据.
2 材料与方法(Materials and methods)2.1 研究区概况研究区位于成都市西部, 东经102°54′~104°53′、北纬30°05′~31°26′, 区域面积约为315 km².属亚热带湿润季风气候, 年平均气温16.2 ℃, 年均降水量918 mm, 区内地势平坦.成都市是西南地区重要的科技、贸易、金融中心.近20年, 成都市快速的经济发展和城市化, 导致城市不断向临域扩张, 扩张方向主要是向成都市西、西北、西南方向(Peng et al., 2015).朝西的城市扩张, 形成了温江区卫星城, 导致成都市中心到温江区之间形成“市中心-郊区-卫星城”梯度带.研究区内有一环、二环、三环和四环4条主要环城道路, 根据城市化程度, 将三环内区域划为中心城, 三环到郫温大道区域划为郊区, 温江城区为卫星城.
2.2 土样采集与分析根据成都市西部梯度分布情况, 将研究区划分为中心城、郊区和卫星城3个梯度带；再按照《城市绿地分类》对研究区域城市绿地划分为：公园绿地、居住绿地、工业绿地、道路绿地和附属绿地.根据不同环带中的各绿地类型, 确定77个采样区, 针对采样区的大小选择不同植被类型(草地、灌木、乔木)的土壤样点, 在每个样点区域10 m²范围内随机选取5处, 用土钻取表层土壤(0~20 cm), 分别将单个地块区域内的5个子样混匀, 作为该样点剖面上的混合样品.共采集公园绿地土壤62个, 附属绿地(学校、公共区域绿地)40个, 居住绿地37个, 工业绿地17个, 交通绿地21个, 共计采样点177个(图 1).野外采集后去除样品中的草根、石块、垃圾杂质等, 样品避光自然风干, 磨碎、过塞备用.采用HNO₃-HF-HClO₄对制备的土壤样品进行消煮, 采用火焰原子吸收分光光度法检测消解液中铜(Cu)、锌(Zn)、铅(Pb)的含量(邱孟龙等, 2015).检测过程所用试剂均为优级纯.检测过程加入国家标准土壤参比物质(GSS-3)进行质量控制, 各元素的回收率均在100%±5%.
图 1(Fig. 1)

图 1 研究区样点分布图 Fig. 1Sample distribution in the study area

2.3 距城市中心距离计算本研究以成都市中心、卫星城(温江城区)作为城市中心, 基于ArcGIS10.2.2计算每个样点到城市中心的距离.中心城、郊区区域样点以成都市中心为城市中心点计算距离, 卫星城区域样点以温江城市中心为城市中心点计算距离.
2.4 建筑密度与道路密度计算基于ArcGIS1 0.2.2运用缓冲区分析计算建筑密度和道路密度.本文以采样点为中心建立0.1、0.3、0.5、0.7、1.0、1.3和1.5 km半径的圆形缓冲区, 计算7个不同尺度下缓冲区内的建筑密度, 同时采用专家评分法, 将高速公路、主要道路、普通道路和乡村道路的权重分别设置为0.4、0.3、0.2和0.1, 以此计算各缓冲区内道路密度(陈青松等, 2016；Liu et al., 2016).
2.5 统计分析描述性统计分析、线性回归分析和逐步回归分析通过SPSS22完成.其中不同绿地类型的土壤重金属含量差异显著性分析采用单因子方差分析(Oneway-AVOVA)和新复极差法(Duncan).逐步回归分析法是通过自变量的显著性选取最优变量的过程, 即定量判定城市化因子解释土壤重金属含量变异的能力(唐梦迎等, 2017), 本研究中影响土壤重金属空间部分的因素包括绿地类型、距城市中心距离、建筑密度和道路密度.其中绿地类型为定性变量, 将居住绿地作为参照变量, 引入D₁、D₂、D₃和D₄ 4个虚拟变量对其进行重新编码, 使其参与逐步回归分析.D₁表示附属绿地系数为1, 其他绿地系数为0；D₂, 交通绿地系数为1, 其他绿地系数为0；D₃, 工业绿地系数为1, 其他绿地系数为0；D₄, 公园绿地系数为1, 其他绿地系数为0.采用Originpro9.0和ArcGIS10.2进行重金属空间特征和影响因素分析以及图形制作.
3 结果与讨论(Results and discussion)3.1 研究区土壤重金属含量的统计特征K-S检验结果表明, 土壤Cu、Pb、Zn呈对数正态分布.土壤Cu、Pb、Zn含量平均值分别为47.77、39.98、89.65 mg·kg^-1, 均未超过国家二级标准(土壤环境质量标准, GB 15618—1995), 但有5个土样Cu含量超过国家二级标准, 存在点源污染.土壤Cu和Pb的平均值均超过姚廷伸(1987)提出的成都平原土壤背景值32.80和25.70 mg·kg^-1, 采样点中有27%的土样Zn含量超过土壤背景值98.30 mg·kg^-1, 说明Cu、Pb、Zn在表层土壤都存在不同程度的富集现象, 并且Cu和Pb的富集更为明显.变异系数反映了总体样本中各采样点的平均变异程度, 能在一定程度反映样品受人为影响的程度.Cu、Pb、Zn的变异系数分别为0.46、0.28、0.34, 根据Wilding提出的变异系数等级分类(Wilding, 1985), 土壤Cu属于高度变异, Zn和Pb属于中等变异, 这与吕建树等(2012)对日照市研究相似, 说明3种重金属均易受到外界环境的干扰, 其中Cu受外界环境影响最大, 这种变异可能由于受到城市发展中工业、交通、商业等人为活动影响(钟晓兰等, 2007).
表 1(Table 1)

表 1 研究区样点土壤重金属含量统计特征表 Table 1 Statistical characteristics of soil heavy metal contents in sample plots

表 1 研究区样点土壤重金属含量统计特征表 Table 1 Statistical characteristics of soil heavy metal contents in sample plots 重金属 最大值/ (mg·kg-1) 最小值/ (mg·kg-1) 均值±标准差/ (mg·kg-1) 背景值/ (mg·kg-1) 超过背景 值占比 国家二级标准/ (mg·kg-1) 偏度 峰度 变异系数 Cu 144.84 26.14 47.77±22.10 32.8 88% 100 2.07 4.73 0.46 Pb 96.51 17.60 39.98±11.28 25.7 93% 300 1.12 6.35 0.28 Zn 202.78 56.26 89.65±30.16 98.3 27% 250 0.78 1.34 0.34

3.2 土壤重金属空间分布特征“中心城-郊区-卫星城”梯度带土壤重金属含量的空间分布存在差异(图 2).土壤Cu、Pb、Zn 3种重金属含量的分布总体均呈东西部高、中部低的趋势, 且成都市中心高于温江城区.中心城到郊区土壤重金属含量呈递减趋势, Pouyat和McDonnell(1991)的研究也表明纽约城市土壤Cu、Ni、Pb的含量从市中区到乡村浓度逐渐降低, 说明城市化影响土壤重金属的空间分布.
图 2(Fig. 2)

图 2 土壤重金属含量随经度变化 Fig. 2Variation of soil heavy metal content with longitude

土壤Zn和Pb的含量在温江城区明显高于绕城路到三环郊区, 但是低于成都市中心.这可能由于随着温江城市化发展, 人类活动增加导致重金属累积, 然而温江属于成都市卫星城, 城市兴起时间晚于成都市中心, 人类活动对土壤重金属的干扰时间相对较短, 因此表现出温江城区土壤重金属含量较成都市中心低.Cu的含量在温江城区没有表现出明显高于郊区, 这可能是因为土壤中Cu含量主要受自然和人为因素共同影响(Zhou et al., 2016), 温江城区土壤与郊区土壤母质相同, 而温江城市化程度低, 人类干扰强度小、时间短, 因此Cu含量在温江城区和郊区没有表现出显著差异.但是在温江城区和郊区之间有3个样点Cu含量偏高, 郊区内有一个样点Zn含量偏高, 这是因为这些样点属于工业绿地土壤, 采集于以铜像制造为主的金属制品厂和有色金属制厂, 受到工业排放影响, 土壤Cu和Zn含量明显高于其他样点, 说明Cu和Zn存在点源污染.
3.3 土壤重金属空间分布的影响因素3.3.1 不同绿地类型对土壤重金属的影响研究区不同绿地类型下土壤重金属的含量具有差异(表 2).工业绿地土壤Cu含量与交通绿地、公园绿地具有显著性差异(p < 0.05), 其含量高低为工业绿地>附属绿地>居住绿地>公园绿地>交通绿地, 说明工业对绿地土壤Cu的含量有较大的影响.不同绿地类型中土壤Pb含量无显著性差异(p>0.05), 交通绿地土壤Pb含量最高, 公园绿地最低.这是因为Pb主要来源于交通污染源和生活废弃物的堆放(郭伟等, 2013), 而汽车尾气排放颗粒物中Pb含量较高的缘故(方晰等, 2010).居住绿地、附属绿地、工业绿地土壤中Zn含量差异不显著(p>0.05), 但显著高于交通绿地、公园绿地(p < 0.05).工业绿地土壤中Zn含量最高, 说明Zn的来源主要是工业排放等人类活动(Liang et al., 2017), 公园绿地Zn含量最低, 究其原因是因为公园绿地土壤环境质量较好, 受工业污染小, 因此公园绿地Zn含量最低.
表 2(Table 2)

表 2 不同绿地类型土壤重金属含量分布特征 Table 2 Distribution characteristics of soil heavy metals in different types of green land

表 2 不同绿地类型土壤重金属含量分布特征 Table 2 Distribution characteristics of soil heavy metals in different types of green land 重金属 含量/(mg·kg-1) 居住绿地 附属绿地 交通绿地 工业绿地 公园绿地 Cu 50.98±19.92abc 56.05±27.89ab 38.67±6.69c 63.95±38.13a 41.28±14.24bc Pb 38.98±10.45ab 37.12±17.39ab 48.03±7.96a 38.88±5.73ab 39.07±9.46ab Zn 107.98±34.99a 110.96±25.81a 71.54±7.09b 113.21±32.33a 69.90±8.20b 注：同一行中相同字母表示差异不显著, 不同字母表示差异显著(p < 0.05).

3.3.2 距城市中心的距离对土壤重金属的影响研究区土壤重金属与距城市中心距离回归分析结果表明(图 3), 土壤Cu、Pb、Zn含量随距城市中心距离增加而显著降低(p < 0.01), 且随距城市中心距离增加, 土壤Cu、Pb、Zn含量降低趋势明显减小.随着距城市中心距离的增加, 城市建设年龄越小, 土壤重金属Cu、Pb、Zn含量逐渐减少, 吴绍华等(2011)研究也表明, 宜兴市土壤Cd、Cu、Pb和Zn含量沿城市扩张方向均表现出明显的递减趋势, 土壤重金属含量随着城市建设时间的增长逐渐增高.Wang等(2017)的研究也表明, 北京城市土壤中的Pb、Cu和Zn浓度与城市年龄线性相关.而道路密度、建筑密度、人口密度等存在的差异是导致重金属元素空间分布的梯度差异的主要原因(吴绍华等, 2011).在距城市中心距离10000 m左右, 存在一个样点Cu和Zn含量异常高于同等距离下的其他样点值, 这是因为该样点采集于郊区的金属制品厂绿地, 受到工业污染影响, 说明工业分布也会影响土壤重金属的空间分布.
图 3(Fig. 3)

图 3 距城市中心距离与土壤重金属Cu、Pb、Zn含量回归分析 Fig. 3Correlation between distance from the urban center and soil Cu, Zn, Pb concentrations

3.3.3 建筑密度对土壤重金属分布的影响不同缓冲区内建筑密度与绿地土壤Cu、Pb、Zn含量线性回归分析结果表明(表 3), 土壤Cu含量随各个缓冲区内建筑密度增加而显著增加(p < 0.01), 其中, 1 km范围内的建筑密度对土壤Cu含量变异的独立解释能力最高(R²=0.170).土壤Pb含量随0.5~1.5 km缓冲区内建筑密度增加而极显著增加(p < 0.01), 且1 km范围内的建筑密度对土壤Pb含量变异的独立解释能力最高(R²=0.104).土壤Zn含量随各个缓冲区内建筑密度增加而显著增加(p < 0.05), 且1.3 km范围内的建筑密度对土壤Zn含量变异的独立解释能力最高(R²=0.116).建筑密度的增加必然会导致城市绿地减少, 人工封闭地面的增加, 一方面城市土壤封闭材料会蓄积重金属(魏宗强等, 2014)；另一方面建筑密度的增加会促使人口密度相应增加, 人类活动增强, 燃煤排放和生活废弃物堆放随之增加, 造成土壤重金属累积(戴彬等, 2015；银燕春等, 2015), 因此土壤重金属含量随建筑密度增加而显著增加.
表 3(Table 3)

表 3 不同距离缓冲区建筑密度与重金属回归分析 Table 3 Building density and heavy metal regression analysis under different distance buffer

表 3 不同距离缓冲区建筑密度与重金属回归分析 Table 3 Building density and heavy metal regression analysis under different distance buffer 重金属 缓冲半 径/km 拟合函数 R2 P 0.1 Y=0.305X+37.113 0.085 0.009 0.3 Y=0.316X+32.078 0.092 < 0.001 0.5 Y=0.374X+27.985 0.133 < 0.001 Cu 0.7 Y=0.378X+27.317 0.136 < 0.001 1 Y=0.421X+24.701 0.170 < 0.001 1.3 Y=0.391X+22.660 0.146 < 0.001 1.5 Y=0.246X+22.276 0.131 < 0.001 0.1 Y=0.089X+37.637 0.008 0.172 0.3 Y=0.180X+34.871 0.024 0.056 0.5 Y=0.270X+31.737 0.065 0.004 Pb 0.7 Y=0.293X+30.862 0.078 0.002 1 Y=0.335X+29.369 0.104 < 0.001 1.3 Y=0.275X+29.622 0.067 0.003 1.5 Y=0.273X+29.005 0.066 0.003 0.1 Y=0.198X+79.663 0.031 0.032 0.3 Y=0.234X+72.721 0.047 0.010 0.5 Y=0.240X+70.750 0.05 0.008 Zn 0.7 Y=0.208X+72.750 0.035 0.022 1 Y=0.212X+72.007 0.037 0.020 1.3 Y=0.315X+56.358 0.116 < 0.001 1.5 Y=0.339X+55.372 0.108 < 0.001

3.3.4 道路密度对土壤重金属分布的影响道路的增加, 一方面导致人工封闭土壤增加；另一方面交通量增加, 汽车尾气排放和轮胎磨损粉尘增加, 造成重金属污染(魏宗强等, 2014；De et al., 2016；Liu et al., 2016).已有研究表明, 道路密度显著影响土壤重金属空间分布, 并且呈现出尺度效应(方淑波等, 2015).
研究区内土壤重金属Cu、Pb、Zn含量与不同缓冲区(0.1~1.5 km)内道路密度线性回归分析表明(表 4), 土壤Cu和Zn含量随各个缓冲区内道路密度增加而显著增加(p < 0.05).其中, 1 km缓冲区内道路密度对土壤Cu和Zn含量变异的独立解释能力最高(Cu, R²=0.119; Zn, R²=0.110).土壤Pb含量随0.3~1.5 km缓冲区内道路密度增加而显著增加(p < 0.05), 且1.5 km缓冲区内道路密度对土壤Pb含量变异的独立解释能力最高(R²=0.111).研究结果表明, 城市区域土壤Cu、Pb、Zn含量受局部范围内的道路密度的影响, 不同重金属含量受到不同范围内道路密度影响, 进一步了解道路密度对重金属含量的影响范围及原因, 对城市土壤管理及重金属污染治理具有重要意义.
表 4(Table 4)

表 4 不同距离缓冲区道路密度与重金属回归分析 Table 4 Road density and heavy metal regression analysis under different distance buffer

表 4 不同距离缓冲区道路密度与重金属回归分析 Table 4 Road density and heavy metal regression analysis under different distance buffer 重金属 缓冲半 径/km 拟合函数 R2 P 0.1 Y=0.206X+38.351 0.034 0.023 0.3 Y=0.204X+36.494 0.034 0.025 0.5 Y=0.260X+32.777 0.060 0.004 Cu 0.7 Y=0.352X+28.172 0.117 < 0.001 1 Y=0.355X+25.856 0.119 < 0.001 1.3 Y=0.336X+25.964 0.106 < 0.001 1.5 Y=0.343X+25.559 0.110 < 0.001 0.1 Y=0.099X+37.780 0.001 0.296 0.3 Y=0.206X+33.993 0.034 0.029 0.5 Y=0.215X+33.280 0.038 0.022 Pb 0.7 Y=0.274X+31.327 0.067 0.003 1 Y=0.291X+29.674 0.076 0.002 1.3 Y=0.324X+27.918 0.097 < 0.001 1.5 Y=0.345X+27.150 0.111 < 0.001 0.1 Y=0.256X+74.137 0.058 0.005 0.3 Y=0.230X+72.299 0.045 0.011 0.5 Y=0.292X+65.856 0.078 0.001 Zn 0.7 Y=0.341X+62.128 0.109 < 0.001 1 Y=0.342X+58.766 0.110 < 0.001 1.3 Y=0.284X+62.608 0.073 0.002 1.5 Y=0.277X+63.207 0.069 0.002

3.3.5 城市化对土壤重金属的综合影响分析结合缓冲区分析, 选取最优缓冲半径内的建筑密度、道路密度, 以及绿地类型、距城市中心距离与重金属含量进行逐步多元回归分析(表6).剔除不显著变量后, Cu的回归模型中只包含了建筑密度, 并解释了Cu在“中心城-郊区-卫星城”梯度带17%的空间变异, 说明城市土壤Cu主要受建筑密度的影响.道路密度和绿地类型一同解释了城市土壤Pb含量在梯度带16%的空间变异, 说明城市土壤Pb含量主要受道路密度的影响, 同时也受到城市格局即不同绿地类型分布的影响, 这与方淑波等(2015)的研究结果一致, Pb的浓度分布不仅响应于城乡的梯度格局, 也与城市化特殊的格局分布相关.绿地类型、道路密度和建筑密度一同解释了城市土壤Zn含量在梯度带43%的空间变异, 其中受绿地类型的影响最大, 说明Zn受人为因素影响较大(Liang et al., 2017).城市化因子解释了土壤Cu、Pb、Zn含量在梯度带16%~43%的空间变异, 说明除绿地类型、距城市中心距离、建筑密度和道路密度外, 土壤重金属的空间分布还受到其它城市化因子的影响.吴绍华等(2011)和方淑波等(2015)研究表明, 城市建设时间、人口密度、交通密度、工业活动以及反映城市化扩张过程的景观形状指数均影响土壤重金属的空间分布.
表 5(Table 5)

表 5 城市化因子与土壤重金属含量关系的逐步多元回归分析 Table 5 Stepwise multiple regression analysis of the relationships between urbanization factors and soil heavy metal

表 5 城市化因子与土壤重金属含量关系的逐步多元回归分析 Table 5 Stepwise multiple regression analysis of the relationships between urbanization factors and soil heavy metal 重金属 模型 R2 P Cu Y=45.77+0.42X1 0.170 ＜0.01 Pb Y=38.65+0.33X2+0.24D1 0.160 ＜0.01 Zn Y=105.33-0.53D4 -0.38D2+0.21X2+0.20X1 0.431 ＜0.01 注：X1, 建筑密度；X2, 道路密度；X3, 距城市中心的距离.D1表示附属绿地系数为1, 其他绿地系数为0；D2, 交通绿地系数为1, 其他绿地系数为0；D3, 工业绿地系数为1, 其他绿地系数为0；D4, 公园绿地系数为1, 其他绿地系数为0.

4 结论(Conclusions)研究区绿地土壤3种重金属统计结果表明, Cu属于高度变异, Pb和Zn属于中等变异, 土壤中Cu和Pb平均含量高于成都土壤背景值, Zn有27%的样点超过背景值, 3种重金属平均含量均未超过国家二级标准, 但有5个土样Cu含量超过国家二级标准, 说明城市化已经造成土壤重金属在一定程度上的累积, 且存在点源污染.
“中心城-郊区-卫星城”梯度带土壤重金属Cu、Pb、Zn含量的空间分布存在差异, 在梯度带上呈现先降低后增加的趋势.卫星城兴起时间晚于中心城, 受人类活动干扰时间较短, 因此表现出卫星城土壤重金属含量低于中心城.郊区人类活动干扰强度较小, 土壤重金属含量低于城市, 但存在点源污染.随距城市中心距离增加土壤重金属含量显著降低(p < 0.01), 城市化影响土壤重金属的空间分布.
城市绿地土壤重金属含量受不同城市化因子影响.研究区不同绿地类型土壤重金属含量具有明显差异, 土壤Zn、Cu在工业绿地含量最高, Pb在交通绿地平均含量高于其他绿地.土壤Cu、Pb、Zn含量受局部范围内的建筑密度、道路密度影响, 且随道路密度、建筑密度增加而显著增加(p < 0.05).逐步回归分析显示, Cu主要受建筑密度影响, Pb主要受道路密度和绿地类型影响, Zn受绿地类型、道路密度、建筑密度影响.
在城市化进程下, 不同绿地类型、距离城市中心距离、建筑密度和道路密度在一定程度上影响“中心城-郊区-卫星城”梯度带土壤重金属的空间分布, 进一步了解城市化进程下“中心城-郊区-卫星城”梯度带土壤重金属来源和去向, 可为城市土壤重金属污染治理、土地资源的可持续利用和生态城市建设提供科学依据.

参考文献

A?ca N. 2015. Spatial distribution of heavy metal content in soils around an industrial area in Southern Turkey[J]. Arabian Journal of Geosciences, 8(2): 1111–1123.DOI:10.1007/s12517-013-1240-7

Bai Y P, Deng X Z, Jiang S J, et al. 2017. Exploring the relationship between urbanization and urban eco-efficiency:evidence from prefecture-level cities in China[J]. Journal of Cleaner Production.

Bren C D, Reitsma F, Baiocchi G, et al. 2017. Future urban land expansion and implications for global croplands[J]. Proceedings of the National Academy of Sciences of the United States of America, 114(34): 8939.DOI:10.1073/pnas.1606036114

陈青松, 李婷, 张世熔, 等. 2016. 城乡交错带土壤氮素空间分布及其影响因素[J]. 生态学报, 2016, 36(8): 2133–2141.

戴彬, 吕建树, 战金成, 等. 2015. 山东省典型工业城市土壤重金属来源、空间分布及潜在生态风险评价[J]. 环境科学, 2015, 36(2): 507–515.

De S S, Ball A S, Huynh T, et al. 2016. Metal accumulation in roadside soil in Melbourne, Australia:Effect of road age, traffic density and vehicular speed[J]. Environmental Pollution, 208: 102–109.DOI:10.1016/j.envpol.2015.09.032

方淑波, 崔曲, 庞慧焕, 等. 2016. 基于源汇理论的浦东新区土壤重金属沿城乡梯度的分布特征[J]. 生态学杂志, 2016, 35(3): 772–780.

方淑波, 崔曲, 庞慧焕., 等. 2015. 上海市城乡梯度上土壤重金属空间分布及其影响因素定量分析[J]. 环境科学学报, 2015, 35(12): 3976–3984.

方晰, 金文芬, 李开志, 等. 2010. 长沙市韶山路沿线不同绿地土壤重金属含量及其潜在生态风险[J]. 水土保持学报, 2010, 24(3): 64–70.

郭广慧, 雷梅, 乔鹏炜. 2015. 北京市城市发展中土壤重金属的空间分布[J]. 环境工程技术学报, 2015, 5(5): 424–428.

郭伟, 孙文惠, 赵仁鑫, 等. 2013. 呼和浩特市不同功能区土壤重金属污染特征及评价[J]. 环境科学, 2013, 34(4): 1561–1567.

李风平, 娄翼来, 李琪, 等. 2008. 沈阳南部沿城乡梯度土壤铅时空分布研究[J]. 土壤通报, 2008, 39(4): 935–937.

李万峰. 2014. 卫星城理论的产生、演变及对我国新型城镇化的启示[J]. 经济研究参考, 2014, 41: 4–8.

Liang J, Feng C T, Zeng G M, et al. 2017. Spatial distribution and source identification of heavy metals in surface soils in a typical coal mine city, Lianyuan, China[J]. Environmental Pollution, 225: 681–690.DOI:10.1016/j.envpol.2017.03.057

Liu A, Gunawardana C, Gunawardena J, et al. 2016. Taxonomy of factors which influence heavy metal build-up on urban road surfaces[J]. Journal of Hazardous Materials, 310: 20–29.DOI:10.1016/j.jhazmat.2016.02.026

Liu X L, Li T, Zhang S R, et al. 2016. The role of land use, construction and road on terrestrial carbon stocks in a newly urbanized area of western Chengdu, China[J]. Landscape and Urban Planning, 147: 88–95.DOI:10.1016/j.landurbplan.2015.12.001

刘亚纳, 朱书法, 魏学锋, 等. 2016. 河南洛阳市不同功能区土壤重金属污染特征及评价[J]. 环境科学, 2016, 37(6): 2322–2328.

柳云龙, 章立佳, 庄腾飞, 等. 2015. "城郊乡"梯度下土壤Cu、Zn、Pb含量的空间变异特征[J]. 长江流域资源与环境, 2015, 24(7): 1207–1213.

柳云龙, 章立佳, 韩晓非, 等. 2012. 上海城市样带土壤重金属空间变异特征及污染评价[J]. 环境科学, 2012, 33(2): 599–605.

吕建树, 张祖陆, 刘洋, 等. 2012. 日照市土壤重金属来源解析及环境风险评价[J]. 地理学报, 2012, 67(7): 109–122.

Macgregor-Fors I. 2011. Misconceptions or misunderstandings? On the standardization of basic terms and definitions in urban ecology[J]. Landscape & Urban Planning, 100(4): 347–349.

Peng W F, Wang G J, Zhou J M, et al. 2015. Studies on the temporal and spatial variations of urban expansion in Chengdu, western China, from 1978 to 2010[J]. Sustainable Cities & Society, 17: 141–150.

Pouyat R V, Mcdonnell M J. 1991. Heavy metal accumulations in forest soils along an urban- rural gradient in Southeastern New York, USA[J]. Water Air & Soil Pollution, 57(1): 797–807.

邱孟龙, 李芳柏, 王琦, 等. 2015. 工业发达城市区域耕地土壤重金属时空变异与来源变化[J]. 农业工程学报, 2015, 31(2): 298–305.

唐梦迎, 丁建丽, 夏楠, 等. 2017. 干旱区典型绿洲土壤有机质含量分布特征及其影响因素[J]. 土壤学报, 2017, 54(3): 759–766.

Wang M E, Liu R, Chen W P, et al. 2018. Effects of urbanization on heavy metal accumulation in surface soils, Beijing[J]. Journal of Environmental Sciences, 64: 328–334.DOI:10.1016/j.jes.2016.11.026

魏宗强, 颜晓, 吴绍华, 等. 2014. 人工封闭对城市土壤功能的影响研究进展[J]. 生态环境学报, 2014, 23(4): 710–715.

Wilding L P. 1985. Spatial variability: its documentation, accommodation and implication to soil surveys[A]. In: Nielsen D R, Bouma J (Eds. ). Soil spatial variability[M]. Wageningen: PUDOC publishers. 166-194

Wu S H, Zhou S L, Yang D Z, et al. 2008. Spatial distribution and sources of soil heavy metals in the outskirts of Yixing City, Jiangsu Province, China[J]. Chinese Science Bulletin, 53(1): 188–198.

吴绍华, 周生路, 潘贤章, 等. 2011. 城市扩张过程对土壤重金属积累影响的定量分离[J]. 土壤学报, 2011, 48(3): 496–505.DOI:10.11766/trxb201006220258

杨蕊, 李小平, 王继文, 等. 2016. 西宁市城市土壤重金属分布特征及其环境风险[J]. 生态学杂志, 2016, 35(6): 1531–1538.

姚廷伸. 1987. 川西平原农业经济自然区冲积性水稻土及粮食作物(小麦、大米)中十一种化学元素背景值的研究[J]. 四川师范大学学报(自然科学版), 1987(增刊): 1–67.

银燕春, 王莉淋, 肖鸿, 等. 2015. 成都市区、城郊和农村生活垃圾重金属污染特性及来源[J]. 环境工程学报, 2015, 9(1): 393–400.

张连科, 张花娟, 黄学敏, 等. 2016. 包头市不同功能区土壤重金属污染评价[J]. 水土保持研究, 2016, 23(2): 352–356.

钟晓兰, 周生路, 李江涛, 等. 2007. 长江三角洲地区土壤重金属污染的空间变异特征——以江苏省太仓市为例[J]. 土壤学报, 2007, 44(1): 33–40.DOI:10.11766/200508260106

Zhou J, Feng K, Li Y J, et al. 2016. Factorial Kriging analysis and sources of heavy metals in soils of different land-use types in the Yangtze River Delta of Eastern China[J]. Environmental Science and Pollution Research, 23(15): 14957–14967.DOI:10.1007/s11356-016-6619-z

庄腾飞, 柳云龙, 章立佳, 等. 2012. 上海城郊土壤重金属含量的空间变异与分布特征——以上海市闵行区为例[J]. 长江流域资源与环境, 2012, 21(z1): 99–104.

中华人民共和国环境保护部. 1995. GB15618-1995土壤环境质量标准[S]. 北京: 中华人民共和国环境保护部