传统城市化带来的城市内涝、水环境污染与地下水缺乏补给等严峻问题，人们通过反思与研究，提出了许多应对之策，如美国的低影响开发(LID)、英国的可持续城市排水系统(SUDS)、澳大利亚的水敏感城市设计(WSUD)等^[1]，均取得了良好效果。中国根据自身的水文、地理条件，提出了建设自然积存、自然渗透、自然净化的“海绵城市”的要求。在中国的海绵城市建设中，雨水花园、生物滞留设施等雨水径流集中入渗设施应用较多。这类措施主要利用入渗性能和污染物吸附性能较好的改良材料作为填料，使降雨径流携带的污染物经填料吸附、过滤、离子交换、微生物降解等物理、化学、生物作用得到自然净化，然后下渗补给地下水或收集再次利用，具有减缓城市内涝、净化水质与涵养地下水资源的综合功能。由于城市雨水径流存在严重的面源污染，且集中入渗的水量负荷与污染负荷强度大，长期集中入渗会造成填料堵塞，污染物吸附饱和等众多问题，降低设施运行效率，缩短设施运行寿命，这些都是众多专家非常关注的科学问题^[2]。

目前，国内外研究大多针对雨水花园、生物滞留设施等LID设施本身的结构、效果、影响机制与模拟研究展开^[3-5]。如王璐等^[6]将雨水花园与工程隔盐技术相结合，通过室内模拟和正交试验方法，提出适于上海滨海盐碱地区的3种雨水花园结构模式。Gurung等^[7]研究表明雨水花园对降雨径流中氮、磷、重金属具有较好的去除效果。罗鸣等^[8]基于下凹绿地和雨水花园构建了SWMM模型，模拟1年一遇、2年一遇、5年一遇重现期暴雨时的出水口径流过程及各污染物负荷总量，表明下凹绿地和雨水花园对径流和水质均有一定的削减作用，且雨水花园作用效果更显著。对低影响开发(LID)模式下雨水径流集中入渗对土壤的影响开展了一些研究，如Kim等^[9]研究了雨水花园中重金属在土壤中的迁移特性，证实了其会对土壤产生污染；Gunawardena等^[10]论述了城市雨水径流中常用的雨水入渗技术和城市雨水中主要的污染物，并对这些污染物入渗后的去向进行了调查；Tedoldi等^[11]认为SUDS的普遍应用会对土壤造成影响，并通过文献回顾分析了解决问题的方法。Tedoldi^[11]和Xie^[12]研究了流体中携带的污染物在土壤环境中的迁移行为和分布规律。但是，近年来，雨水花园集中入渗雨水径流对土壤影响过程的量化分析缺乏系统研究，影响了这类措施的合理使用。

本文对运行8~9 a的雨水花园土壤进行多次取样分析，重点研究黄土地区雨水径流集中入渗对土壤氮、磷、重金属含量的影响过程，以期为黄土地区雨水径流集中入渗措施的合理配置与推广应用提供科学依据和理论支撑，促进中国海绵城市的健康发展。

1 研究区概况

试验雨水花园位于陕西省西安市西安理工大学校园内。西安市位于黄河流域中部的陕西关中盆地，属于大陆性季风气候，平均气温13.3 ℃，冷热干湿，四季分明，冬季干冷，春季干燥，夏季湿热，秋季多雨，属半湿润气候区。根据西安市1951—2008年(58 a)的降雨统计资料，多年平均降雨580.2 mm，2017年研究区降水总量为642.9 mm，降雨量年际变化相差较大，年最大年降雨量达903.2 mm(1983年)。西安市是缺水城市，人均水资源占有量仅为全国的1/6，且地下水超采严重。西安市所在地区的黄土层极为丰厚，黄土入渗性能非常好，且为入渗、滞留自然降雨径流提供了天然的屏障。故以LID设施为单个海绵体，用于排除和蓄存雨水径流是海绵城市建设过程中重要的措施，为缓解城市内涝灾害和补充地下水资源匮缺提供良好的天然条件。本研究涉及两个雨水花园RD₁和RD₂。

雨水花园RD₁(图 1a)建成于2011年，用于处理办公楼屋面雨水径流，屋面面积，即汇水面积约为605 m²，雨水花园面积为30.24 m²，汇流比为20∶1。上部蓄水层15 cm，下部结构层填充西安市本地黄土20 cm，渗透系数大约为2.346 m/d，结构层以下为原状土层(未受扰动)，花园底部不做防渗处理，雨水径流经土壤基质直接入渗补给地下水。入流口安装45°三角堰，溢流口安装30°三角堰。雨水花园RD₁东侧30 m处有一垃圾场，主要用于收集校园内生活垃圾以及枯枝落叶等，垃圾场底部为水泥混凝土硬质地面，垃圾渗滤液不会发生垂向或水平渗漏，主要是清理不及时带来的少量垃圾碎屑沉积物随径流雨水冲刷而进入花园内。

雨水花园RD₂(图 1b)建成于2012年，大致为椭圆形，长抽为6 m，短轴为2 m，雨水花园面积约为9.42 m²，主要汇集路面和屋面的混合雨水，汇水面积约为140 m²，汇流比为15∶1。花园中间用隔板分割为两个面积相同的雨水花园，其中一侧做防渗处理，底部埋设穿孔排水管，排水管周围填筑15 cm砾石反滤层；另一侧不做防渗处理，雨水径流可直接入渗补给地下水。在花园入流口安装两个45°三角堰，出水口和溢流口安装30°三角堰。上部蓄水层50 cm，下部填充西安市本地黄土60 cm，底层为原状土层(未受扰动)。雨水花园结构见图 2。

2 试验方法

每次降雨开始时在入流、出流、溢流三角堰处记录各场次降雨的入流、出流瞬时流量，降雨初期一般每隔5 min采集一次进、出水水样，在降雨持续20 min后可按照每隔10 min采集1次进、出水样，降雨持续1 h后按照每隔20 min或30 min采集1次进、出水样，一般降雨水样采集数量大于8个，现场采集的水样及时放入-4 ℃冰箱内，水样分析指标包括TSS，COD，NH₄-N，NO₃-N，TN和TP，水质指标一般在3~5 d内分析完毕。

2017年4月至2019年2月在两个雨水花园RD₁和RD₂(入渗一侧)内共采集了7次土样(距离花园RD₁西侧5 m处的绿化地中采集一对照土样CK)，土壤含水率保持在11.3%~30.1%，pH值保持在7.55~7.76，呈微碱性，主要测定土壤N，P，TOC，重金属含量。为了保持取样的一致性，一般定于降雨后2~3 d内进行土样采集，采用3点(花园RD₁为同心圆上3点，花园RD₂为直线上3点)取样混合法采集土样。花园RD₁土层深度较深，分5层进行采集，分别为0—10，20—30，40—50，70—80，90—100 cm；雨水花园RD₂土层较浅(60 cm)，分3层进行土样采集，分别为0—10，20—30，40—50 cm；对照组CK分4层采集，分别为0—10，20—30，40—50，70—80 cm。样品采集后分成两部分，一部分立即进行含水率测定，另一部分晾在阴凉处，在自然状态下进行通风阴干，剔除其中的草根等杂物，将其充分混匀，采用四分法取少量阴干土壤，过2 mm筛，测定NH₄-N，NO₃-N，另一部分过1 mm筛测定TN，TP，重金属Cu，Zn，Cd，超出最小检测线值记为0。经处理后的土样放入-20 ℃冰箱内待测，两周内分析完毕。7次土样采集前降雨情况见表 1，水样和土样分析方法见表 2。重金属Cu，Zn，Cd的形态采用经典的Tessierr等连续5级提取法^[13]，包括可交换态(S₁)、碳酸盐结合态(S₂)、铁—锰氧化物结合态(S₃)、有机结合态(S₄)和残渣态(S₅)。

3 结果与分析 3.1 雨水花园入流污染物浓度和负荷分析

2016—2018年对雨水花园RD₁和RD₂监测了16~18场降雨事件，主要采集雨水花园入流水样，测定其径流中污染物含量，其径流中各污染物EMC浓度和雨水花园单位面积受纳污染负荷量见表 3。

由表 3可以看出，雨水花园RD₁入流中COD，TSS，NH₄-N，TN，TP等污染指标的平均EMC浓度较RD₂大，但NO₃-N和重金属Cu，Zn，Cd含量均小于RD₂。其中RD₁中COD的EMC浓度是RD₂的1.80倍，单位面积受纳污染负荷量是RD₂的2.44倍。说明RD₁径流受有机物污染较为严重，降雨径流进入雨水花园后将破坏土体平衡，COD的降解将会消耗土体中氧气，成为厌氧状态，影响生物生存^[14]。雨水花园RD₁东侧30 m处有一校园垃圾场，垃圾场3面设有高2 m的水泥围墙(敞开一面背离雨水花园)，受校园垃圾日常清扫、垃圾装卸等活动的影响，雨水花园RD₁汇水渠中堆积的粉尘较多，夏季降雨突然来不及清扫时，粉尘随径流会进入RD₁中，导致雨水花园RD₁径流中TSS，N，P等污染物含量较高。

雨水花园RD₂中重金属含量较高，其中重金属Cu，Zn的EMC含量分别是RD₁的2.41，2.24倍，单位面积受纳污染负荷量分别是RD₁的1.81，1.69倍，说明重金属含量主要来自路面径流。主要是由于城市道路行驶车辆较多，汽车尾气排放、油脂泄露、车胎磨损将导致路面沉积物中重金属Zn含量升高。韦毓韬等^[15]研究表明，道路雨水径流中Cu，Zn，Cd，Pb等重金属主要受交通因素的影响。

3.2 土壤N，P，TOC含量垂向分布

雨水花园0—100 cm不同土层深度N，P含量如图 3所示。

由图 2可以看出，雨水花园和CK土壤中NH₄-N和TOC含量随土层深度增加呈逐渐递减的趋势，上层土壤中NH₄-N和TOC含量明显高于下层土壤。其中RD₁土壤中NH₄-N含量由表层下10 cm的67.69 mg/kg降低至表层下100 cm的21.75 mg/kg，降低了211.18%。这主要是由于降雨径流中NH₄-N主要以颗粒态形式存在，进入雨水花园后被拦截在表层土壤。此外，NH₄-N主要带正电荷，而土壤大多带负电荷，雨水径流中NH₄-N易吸附在土壤胶体表面，在降雨冲刷作用下不易随水分向下迁移。RD₁土壤中TOC含量由表层下10 cm的13.83 g/kg降低至表层下100 cm的5.42 g/kg，降低了155.38%。随取样深度增加，雨水花园和CK土壤中NO₃-N含量呈增加趋势，上层土壤中NO₃-N含量小于下层土壤。RD₁土壤中NO₃-N含量由表层下10 cm的65.84 mg/kg增加至表层下100 cm的70.23 mg/kg，下层土壤NO₃-N含量比上层增加了6.68%。这就说明降雨径流中NO₃-N进入土壤后发生了淋溶现象，NO₃-N随水分入渗易向下迁移。王禄等^[16]通过人工滤柱试验发现，NH₄-N主要集中表层下50 cm土壤，并随滤池深度增加呈逐渐降低的趋势，而NO₃-N随水流被快速排除试验系统。有关研究表明^[17]，硝酸盐作为阴离子不被土壤介质所吸附，在土壤/水系统中通常具有很强的移动性。与NH₄-N含量的变化规律一致，雨水花园土壤中TN含量随土层深度增加呈逐渐减小的趋势。雨水花园土壤中TP含量随土层深度增加逐渐增大，而CK土壤中TP含量随深度增加呈减小趋势。RD₁不同层深度TP含量分别为250.56，239.93，273.09，366.49，485.83 mg/kg，表层下100 cm土壤TP含量较表层下10 cm增加了93.90%。说明雨水花园土壤中可溶性磷随水分入渗发生了明显的迁移作用，具有显著的深层富集现象。因此，雨水花园、生物滞留设施不宜用于地下水位较高的地区，若确实需要应用，需做好前处理，避免硝酸盐和磷酸盐的淋溶而发生次生灾害。

3.3 土壤重金属垂向分布

将雨水花园2017年4月至2019年2月采集的7次土样中3种重金属含量分别取均值，得到雨水花园土壤重金属含量在垂向上的分布规律(图 3)。

由图 3可以看出，雨水花园和CK土壤重金属Zn和Cd大多富集在土壤表层下30 cm范围内，并且土壤Cd含量随土层深度增加呈逐渐减小的趋势。土壤中重金属Cu的分布较均匀，RD₁的90—100 cm和RD₂的40—50 cm土层中Cu含量略大于其他各层，而CK表层土壤重金属Cu含量略大于其他各层。研究表明，城市道路径流雨水中重金属Zn的主要来源于车辆轮胎磨损、发动机润滑油遗漏；重金属Cu的主要来源于金属电镀、车辆轴承及制动部件磨损；重金属Cd主要来源于车辆轮胎磨损^[18]。董微砾^[19]研究新疆乌鲁木齐市城市道路雨水中重金属Cu平均值为209.7~237.4 μg/L，Zn平均值为530.4~745.5 μg/L，Cd平均值为0.7~3.3 μg/L。Xie等^[12]通过大量研究认为Zn在市区的主要来源是车辆的轮胎磨损、制动磨损、大气干湿沉积和石油、天然气的燃烧作用。

降雨径流导致了雨水花园土壤中重金属含量的增加，0—50 cm土层，雨水花园RD₁和RD₂土壤Zn含量分别为81.58，70.66，59.73，56.93，58.48 mg/kg和77.36，84.74，92.13，99.18，106.23 mg/kg，可以看出，除表层土壤(0—10 cm)外，其余各层土壤，雨水RD₂中Zn含量分别是RD₁的1.20，1.54，1.74，1.68倍，表明雨水花园土壤Zn含量受降雨汇水面类型的影响较大，汇水面为道路的雨水花园，其土壤中重金属Zn含量明显高于汇水面为屋面的雨水花园由于RD₂主要接纳路面径流，故RD₂入流水质中Zn浓度远高于RD₁。

不同汇水面对雨水花园土壤中重金属Cd的影响也不同，0—50 cm各层土壤，雨水花园RD₁土壤Cd含量分别为0.581，0.531，0.363 mg/kg，RD₂土壤重金属Cd分别为0.237，0.147，0.256 mg/kg，而CK土壤分别为0.277，0.165，0.233 mg/kg。可以看出，雨水花园RD₂和CK土壤中Cd含量差异较小，但RD₁土壤Cd含量远大于RD₂和CK，分别是RD₂和CK的2.45，3.62，1.42，2.09，3.22，1.56倍。分析其原因，主要是受RD₁东侧的垃圾场影响。雨水花园RD₁主要汇集屋面径流，降雨径流经坡屋面自由落水进入地面汇水沟后集中进入雨水花园，而垃圾场西侧正好为雨水花园部分地面汇水沟，每天校园内清理垃圾时，有大量粉尘落入东侧地面汇水沟和雨水花园RD₁，粉尘中可能含有重金属Cd，在降雨冲刷作用下进入了雨水花园，导致雨水花园RD₁中Cd含量远大于RD₂。

3.4 土壤重金属与TOC的关系

由图 4可以看出，土壤中重金属Cu和Zn与土壤有机碳(TOC)呈直线关系，R²分别达0.840，0.861，重金属Cd与TOC呈二次曲线关系，R²达0.802。TOC是通过微生物作用所形成的腐殖质、动植物残体和微生物体的合称，土壤TOC是雨水花园径流雨水中重金属的重要载体，TOC含量较高的土壤可有效削减径流雨水中的重金属，从而减少重金属对土壤中动、植物的危害，并且可防止重金属向下迁移，影响地下水水质。王腾云等^[20]通过福建沿海地区土壤—稻谷重金属含量关系与影响因素研究，发现富含有机碳的土壤条件有利于阻断稻谷对土壤Cd的吸收，降低土壤Cd污染的生态风险。谢娜等^[21]通过研究不同土地利用方式土壤有机碳变化特征及与重金属的相关性，发现典型土壤重金属Hg，Cd，Cr，As，Pb与SOC储量之间存在显著相关性，园地和建设用地最高r=0.99，林地和耕地次之，分别为r=0.87和r=0.86。可见土壤中TOC与土壤重金属具有密切相关性。

3.5 土壤重金属形态分布

可交换态和碳酸盐结合态两种形态的重金属不大稳定，随水分入渗易向下迁移，铁—锰氧化物结合态、有机结合态和残渣态在土壤中的存在形态较稳定。重金属Cu，Zn，Cd的各形态含量如图 5所示。

由图 5可以看出，雨水花园和CK土壤中重金属Cu的5种形态关系为：S₅ > S₃ > S₄ > S₂ > S₁，而Zn的5种形态关系为：S₃ > S₅ > S₂ > S₄ > S₁。故雨水花园土壤中重金属Cu和Zn主要以铁—锰氧化物结合态和残渣态形式存在，存在形态相对稳定，随水分入渗不易发生向下迁移。雨水花园和CK土壤中重金属Cd的5种存在形态与Cu和Zn明显不同，其关系为：S₁ > S₂ > S₃ > S₄ > S₅，由此可见雨水花园土壤中Cd主要以可交换态和碳酸盐结合态形式存在，这就说明Cd在土壤中很不稳定，随水分入渗有向下迁移的风险，易进入地下水。

4 结论

(1) 雨水花园汇水面类型影响降雨径流中污染物浓度，雨水花园RD₁入流中COD，TSS，N，P等污染物含量较高，而雨水花园RD₂入流中重金属含量远高于RD₁。

(2) 受降雨径流集中入渗污染影响，雨水花园0—50 cm土壤NH₄-N和TN含量较高。土壤NH₄-N，TN和TOC随土层深度增加逐渐减小，而NO₃-N和TP随土层深度增加逐渐增大，发生了明显的淋溶。

(3) 雨水花园土壤重金属大多集中在0—30 cm范围。Cu和Zn主要以铁—锰氧化物结合态和残渣态形式存在，Cd主要以可交换态和碳酸盐结合态形式存在，随水分入渗易向下迁移。雨水花园土壤重金属与TOC具有较好的拟合关系(R² > 0.8)。