环境科学  2016, Vol. 37 Issue (4): 1379-1388   PDF    
重庆南山表层岩溶泉与地下河三氮运移及氮通量估算
张远瞩1, 贺秋芳2,3 , 蒋勇军1, 李勇1    
1. 西南大学地理科学学院, 重庆 400715;
2. 国土资源部/广西岩溶动力学重点实验室, 桂林 541004;
3. 中国地质科学院岩溶地质研究所, 桂林 541004
摘要: 岩溶区特殊的地质构造使地下水系统存在多重水流. 为研究三氮在表层岩溶带水流和地下河水流中的特征和运移方式,选取城市化进程中的重庆南山老龙洞地下河流域表层岩溶泉和地下河按月采样,分析水化学特征,结合SPSS的相关性分析,认为两种水流三氮特征差异很大. 表层岩溶泉三氮质量浓度月变化小,受降水、污水影响较小,NO3--N>>NH4+-N; 地下河水三氮质量浓度月变化大,受降水、污水影响较大,旱季NH4+-N > NO3--N,雨季NO3--N > NH4+-N. 表层岩溶泉DIN主要来源于与农业活动有关的非点源污染,在土壤、表层岩溶带中经氨化、吸附、硝化作用等过程以扩散流形式运移至出露地表. 地下河DIN随管道流、扩散流运移,并以管道流为主. DIN全年有点源的工业、生活污染物经落水洞、裂隙、溶隙进入. 雨季同时存在占更大比例的非点源污染物,或由降雨产生的表面坡流、壤中流、表层岩溶带水流经落水洞进入,或经深裂隙、溶隙下渗进入,硝化作用明显. 流域DIN输出通量为56.05kg ·(hm2 ·a)-1,其中NH4+-N、NO3--N分别占46.03%、52.51%; 根据径流分割法估算出点源污染、非点源污染的贡献率分别为25.08%、74.92%.
关键词: 表层岩溶泉     岩溶地下河     三氮运移     氮通量     重庆南山    
Characteristics and Transport Patterns of Ammonia, Nitrites, Nitrates and Inorganic Nitrogen Flux at Epikarst Springs and a Subterranean Stream in Nanshan, Chongqing
ZHANG Yuan-zhu1, HE Qiu-fang2,3 , JIANG Yong-jun1, LI Yong1    
1. School of Geographical Sciences, Southwest University, Chongqing 400715, China;
2. Karst Dynamics Laboratory, Ministry of Land and Resources & Guangxi, Guilin 541004, China;
3. Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004, China
Abstract: In a karst groundwater system, it develops complex multiple flows because of its special geological structure and unique physical patterns of aquifers. In order to investigate the characteristics and transport patterns of ammonia, nitrite and nitrate in epikarst water and subterranean stream, the water samples were collected monthly in a fast-urbanizing karst region. The results showed distinctive characteristics of three forms of inorganic nitrogen. The concentration of inorganic nitrogen was stable in the epikarst water while it was fluctuant in the subterranean stream. Epikarst water was less affected by rainfall and sewage compared with subterranean stream. In epikarst water, the nitrate concentration was much higher than the ammonia concentration. Dissolved inorganic nitrogen, mainly from non-point source pollution related to agricultural activities, passed in and out of the epikarst water based on a series of physical, chemical and biological processes in the epikarst zone, such as ammonification, adsorption and nitrification. On the contrary, subterranean stream showed a result of NH4+-N >NO3--N in dry seasons and NO3--N >NH4+-N in rainy seasons. This can be due to the fact that sanitary and industrial sewage flowed into subterranean river through sinkholes, fissures and grikes in dry season. Dissolved inorganic nitrogen in subterranean river was mainly from the non-point source pollution in wet season. Non-point source pollutants entered into subterranean water by two transport ways, one by penetration along with vadose flow through fissures and grikes, and the other by conduit flow through sinkholes from the surface runoff, soil water flow and epikarst flow. The export flux of DIN was 56.05 kg ·(hm2 ·a)-1, and NH4+-N and NO3--N accounted for 46.03% and 52.51%, respectively. The contributions of point-source pollution and non point-source pollution to the export flux of DIN were 25.08% and 74.92%, respectively, based on run-off division method.
Key words: epikarst spring     subterranean stream     transport patterns of nitrogen     inorganic nitrogen flux     Nanshan,Chongqing    


地下水是干旱缺水的西南岩溶区社会经济发展的重要支撑,但随着城市的发展和工农业活动的加强,岩溶地下水水质不断恶化,三氮成为常见污染物. 三氮即氨氮(NH4+-N)、 硝酸盐氮(NO3--N)和亚硝酸盐氮(NO2--N),合为溶解无机氮(DIN). 目前国内外地下水三氮运移的研究主要集中在非岩溶区,对三氮在水土界面、 包气带、 含水层中迁移转化机理、 过程及规律进行了定性、 定量研究[1, 2, 3, 4, 5, 6, 7, 8, 9, 10].

岩溶区因其特殊的地质构造使地下水系统与非岩溶区差异巨大,是一个存在裂隙、 孔隙、 管道和洞穴等多重水流的非常特殊的混合系统,其含水介质具有高度异质性. 岩溶含水层从上至下包括上覆土壤层、 表层岩溶带、 饱水带[11, 12]. 大气降水一部分经土壤层渗透进入表层岩溶带,经孔隙、 裂隙、 溶隙,以扩散流形式运移,部分以表层岩溶泉形式出露地表; 一部分降水在地表形成坡流,经落水洞或竖井以集中流形式直接进入地下河,最终在地表出露形成岩溶泉[13],二者都是岩溶区重要水资源,但其补给(扩散流和管道流)、 径流(非饱和带和饱水带)和排泄(扩散流和管道流)途径均不相同,三氮的特征及转化运移规律也应有较大差异.

地表与含水层的连通性是影响污染物迁移的重要因素:以扩散流形式运移的表层岩溶水具有渗透率低[14]、 水流速度慢、 水物理化学指标较为稳定的特点[15],其中污染物将经历长期的储存、 运移过程再出露地表[16, 17]; 相反,与地表高度连通的管道流具有渗透率高[14]、 水流速度快、 水物理化学指标极不稳定的特点[15],污染物通常随暴雨直接经落水洞进入地下水系统[18],在含水层中快速、 远距离运移[19, 20],几乎没有机会发生化学变化和生物降解作用[21, 22],自净能量有限[23]. 一些学者对NO3-在表层岩溶带[24]和地下河[25, 26]的运移进行研究,并对地下河中NO3-的运移进行模拟[27]. 但污染物中DIN还以NH4+、 NO2-等形式进入地下水,同一区域内,表层岩溶泉和地下河中三氮的季节动态变化特征不清楚,特别是对二者迁移转化规律的对比研究较少. 因此,本研究通过监测城市化进程中的重庆南山地下水系统的表层岩溶泉和地下河地球化学指标月动态变化,分析三氮在这两类地下水的不同季节动态变化特性,揭示三氮在两类地下水中不同的储存、 转化和运移规律,并在此基础上分析点源、 非点源污染物对两类地下水三氮的贡献. 这将对岩溶区不同储存形式的地下水三氮污染的治理和合理开发具有积极的指导意义.

1 研究区概况

研究区位于重庆市南岸区与巴南区境内,东经106°34′44″~106°36′48″,北纬29°28′35″~29°33′13″(图 1),流域面积13.09km2,处于川东平行岭谷区南温泉背斜,为一典型的岩溶槽谷,槽谷为近南北走向,地势北高南低. 背斜核部是下三叠统嘉陵江组(T1j)灰岩; 两翼为中三叠统雷口坡组(T2l)碳酸盐岩和上三叠统须家河组(T3xj)砂岩,石灰岩分布面积约11.3 km2. 背斜核部纵张裂隙发育形成含水层,两翼带状分布的砂岩为相对隔水层. 地下水补给主要来源于降水,由表层岩溶泉和地下河出口排泄. 核部碳酸盐岩遭溶蚀后形成洼地,洼地有较厚的覆盖层,边缘常有表层岩溶裂隙泉出露. 槽谷底部落水洞成串珠状分布,发育着一条自NNE向SSW流动的地下河,总长约6 km,最大流量达13 000 L ·s-1,最小流量达0.6 L ·s-1,常流量约50~80 L ·s-1.

图 1 研究区水文地质示意 Fig. 1 Schematic hydrogeological map of Laolongdong underground water system

该区属亚热带湿润气候,年均降水1 100 mm,主要集中在4~10月. 多年平均气温18.7℃. 强烈的岩溶作用造成该区生境破碎,快速的城市化进程更使土地利用类型日益破碎,林地和草地分别占总面积的47.19%和4.18%,建设用地和农业用地分别占27.22%和8.23%(图 2). 城市化也使原有排污管道不能满足要求,大量生活、 工业废水直接排入或被雨水冲入洼地,或通过与岩溶洞穴管道有较好水力联系的落水洞[28]直接进入地下河系统. 研究区生猪、 鸡鸭养殖基地,产生大量排泄物并露天堆放. 以蔬菜种植为主的农业长期大量施用人畜粪便及尿素、 碳铵、 复合化肥等. 生活污水、 工业废水、 人畜粪便、 化肥是研究区三氮的主要来源.

图 2 研究区土地利用类型 Fig. 2 Land uses of Laolongdong underground water system
2 研究方法 2.1 在线数据监测

将美国ONSET公司生产的HOBO水位温度自动记录仪(精度0.3 cm)安装在地下河出口处水文站,水文站与修建的长约10 m,宽2.1 m,高1.5 m的矩形明渠相通,15 min自动记录一次水位,运用经验公式将水位换算成流量; 将HOBO小型气象站(精度0.2 mm)安装在地下河出口处楼顶,30 min自动记录一次降雨量(图 1).

2.2 样品采集及测试

将青岛崂山电子仪器总厂有限公司生产的SYC-3型冷藏型降水、 降尘自动采样器安装在气象站边,自动采集降水,收集桶直径300 mm. 降水时自动打开(灵敏度0.2 mm ·h-1),降水停止10 min 后会自动关闭,每次降水结束后进行现场检测和样品采集. 研究期间除2014年2、 12月因降水量小未能收集到外,其余各月均采集到降水样品,共16场降水,能较好地代表全年降水特征. 选取水面离地表深度均小于1.5 m的表层岩溶泉形成的井水W1、 W2、 W3、 W4(图 1)(W1、 W2位于农业居民区,W3、 W4位于农业菜地)、 地下河出口G、 岩溶洞穴内深部裂隙水F、 污水点S1(城镇生活污水)、 S2(农业排水及生活污水),于2014年2月至2015年1月按月监测并采集水样. 现场监测(包括降水)使用美国HCAH公司的HQ340d多参数水质分析仪测定各水体的pH、 水温(T)、 电导率、 溶解氧(DO)(精度分别为0.01pH单位、 0.1℃、 1 μS ·cm-1、 0.01 mg ·L-1); 用D850水质分析仪测定NH4+、 NO2-质量浓度(精度为0.01mg ·L-1和0.001 mg ·L-1). 采集的水样(包括降水)经0.45 μm微孔滤膜过滤后,置入清洗过的500 mL的聚乙烯瓶中,用于NO3-测试. 所有样品取样后,立即密封保存于便携式冰袋中,12 h内运至室内放入4℃冰箱中并尽快测试. NO3-质量浓度用紫外分光光度法,采用日本岛津公司生产的UV2450紫外分光光度计测试(检测精度优于0.01 mg ·L-1),在西南大学地球化学与同位素实验室完成.

2.3 数据分析处理

4个表层岩溶泉(W1、 W2、 W3、 W4)数据算术平均后得到表层岩溶泉各项理化指标,地下河各指标按月流量加权[公式(1)]计算年均值. 降水各理化指标根据降水量加权平均,按公式(2)和(3)分别求得月均值和年均值. 城镇生活污水(S1)和农业区排水沟(S2)均以人畜粪便为主要氮源,具有相似的组成和变化规律,将其理化指标算术平均后作为污水指标(表 1). 采用Excel和SPSS 19.0软件进行数据计算分析.

式中,CG年为地下河T(℃)、 pH、 电导率(μS ·cm-1)、 DO(mg ·L-1)、 三氮质量浓度(mg ·L-1)的年平均值; i为1~12月; Ci为第i月的T、 pH、 EC、 DO、 三氮质量浓度; Qi为第i月流量(m3).
式中,CR月为降水T(℃)、 pH、 电导率(μS ·cm-1)、 DO(mg ·L-1)、 三氮质量浓度(mg ·L-1)的月平均值; i为各月降雨次数; Ci为第i次降雨的T、 pH、 电导率、 DO、 三氮质量浓度; Pi为第i次降雨的降雨量(mm).
式中,CR年为降水T(℃)、 pH、 电导率(μS ·cm-1)、 DO(mg ·L-1)、 三氮质量浓度(mg ·L-1)的年平均值; i为1~12月; Ci为第i月降雨T、 pH、 电导率、 DO、 三氮质量浓度的月平均值; Pi为第i月降雨量(mm).

表 1 部分水文地球化学特征指标统计 Table 1 Statistical analysis of hydrogeochemical features (partly)
3 结果与分析 3.1 表层岩溶泉、 地下河三氮质量浓度月变化规律

从年均值看,表层岩溶泉电导率略大于地下河(表 1),因为表层岩溶水经过了较长的存储时间,与围岩接触更充分,能溶解更多碳酸盐岩. 地下河DIN质量浓度较表层岩溶泉大,证明有更多的氮素进入地下河. 地下河DO小于表层岩溶泉,成为一个相对缺氧的环境,对NH4+的硝化作用将产生一定的抑制作用,这与地下河具有较大质量浓度的NH4+、 表层岩溶泉具有较大质量浓度NO3-结果一致[图 3(a)、 3(b)、 3(e)]. 表层岩溶泉NH4+-N/NO3--N几乎接近于0[图 3(f)],说明NH4+-N在表层岩溶带迁移过程中已经过硝化过程,完成了向NO3--N的转化. 地下河DIN的主要赋存形式为NH4+-N,其质量浓度是地下水Ⅴ类水下限(0.5mg ·L-1)的7.6倍,NO3--N平均质量浓度低于表层岩溶泉,NH4+-N/NO3--N远远大于表层岩溶泉[图 3(f)],且最大、 最小值出现时间与表层岩溶泉不一致,说明两种地下水三氮的转化运移方式不同.

(a)表层岩溶泉三氮组成月变化; (b)地下河三氮组成月变化; (c)污水三氮组成月变化; (d)降水三氮组成月变化; (e)表层岩溶泉、 地下河、 污水、 雨水DIN月变化; (f)表层岩溶泉、 地下河、 污水、 雨水NH4+-N/NO3--N月变化图 3 DIN组成及月变化 Fig. 3 Composition and monthly variation of dissolved inorganic nitrogen

从季节变化看,表层岩溶泉和地下河水中DIN质量浓度的最大值均出现在旱季,最小值均出现在雨季[图 3(e)],体现出降水对地下水DIN的稀释作用. 但二者也存在显著差异[图 3(a)、 3(b)、 3(e)]. 表层岩溶泉相对稳定,无明显季节差异,而地下河有明显季节变化规律:旱季地下河NH4+-N>NO3--N,证明地下河中有大量未经硝化作用的新鲜污染物输入. 雨季NO3--N>NH4+-N,NH4+-N/NO3--N远远小于旱季[图 3(f)]. 说明地下河三氮在旱季和雨季具有不同的来源与转化运移规律.

地下河中NO2--N质量浓度较高,监测的12月中有8个月超过国家地下水Ⅴ类水下限(0.1 mg ·L-1),月变化较大,最大值出现时间与NO3--N一致,体现出地下河NO3--N、 NO2--N在运移转换过程中具有一定的联系.

3.2 地下河DIN通量估算

通过水位温度自动记录仪获取15 min一次的水位数据,再根据水位-流量关系计算出15 min一次的流量,累积得到每月的总流量,用公式(4)[29]求出的三氮年输出量.

式中,W为地下河三氮年输出量(kg),C为三氮月质量浓度(mg ·L-1),Q为月流量(m3).

1.表层岩溶带; 2.岩溶管道; 3.落水洞; 4.土壤层; 5.壤中流; 6.管道流; 7.表层岩溶水流; 8.表面坡流; 9.表层岩溶泉; 10.岩溶大泉; 11.裂隙; 12.溶隙; 13.污水裂隙、 溶隙流; 14.污水地表流图 4 表层岩溶泉、 地下河DIN运移路径概念模型 Fig. 4 Conceptual model of transport patterns of dissolved inorganic nitrogen at epikarst springs and subterranean stream

在流域面积为13.09 km2的研究区内,地下河DIN输出通量为56.05 kg ·(hm2 ·a)-1,远大于典型农耕区的寨底地下河的7.94 kg ·(hm2 ·a)-1[29]. 其中,NH4+-N、 NO3--N输出通量分别占总输出量的46.03%、 52.51%,其他少量以NO2--N的形式输出,而寨底地下河NO3--N输出量占DIN的99.25%. 体现出城市化进程使地下河的污染更严重,来源更复杂. DIN输出主要集中在降水量大的3~9月,占总输出量的82%以上,而并非DIN浓度高的旱季. 输出量最大的月份出现在降大暴雨9月(图 5),说明在雨季有更多氮素通过不同路径进入地下河.

图 5 地下河基流、 DIN点源输入分割 Fig. 5 Separation of base flow and point source of DIN
4 讨论 4.1 表层岩溶泉三氮的来源及转化运移

表层岩溶泉三氮质量浓度相对稳定,变异系数低(表 1),因其主要在裂隙、 溶隙中以扩散流的形式运移,运移速度慢、 时间长,对三氮有一定的调节作用,并为生物化学反应提供了机会. 相反,降水、 污水变异系数高,DIN的主要赋存形式和月变化规律与表层岩溶泉也差异很大[图 3(a)、 3(c)、 3(d)、 3(e)],说明污水不是表层岩溶泉三氮的主要来源,降水也与表层岩溶泉三氮无直接联系. 因此,与农业活动相关的人畜粪便和化肥可能成为表层岩溶泉三氮的主要来源.

研究区表层岩溶泉大多出露于核部洼地边缘,发育有较厚的覆盖层(主要为粉质黏土),其土地利用方式主要是农业菜地和分散的农业居住用地(图 2). 菜地中化肥、 粪肥大量使用,养殖场牲畜粪便未经处理,直接露天堆放. 化肥和动物粪便中的大量氮化物在降水驱动下进入土壤,被微生物氨化为NH4+-N,NH4+-N一部分被植物吸收,一部分经微生物硝化作用转化为NO3--N. 未被硝化的NH4+-N因研究区碱性的粉质黏土对带正电的NH4+具有很强的吸附作用[30],使其不能随入渗水进入表层岩溶带,造成表层岩溶泉中很小质量浓度的NH4+-N. 硝化作用形成的NO3-因带负电而不易被碱性的粉质黏土吸附,植物吸收后过量部分随入渗水进入表层岩溶带,经过表层岩溶带中裂隙、 溶隙的调节后以表层岩溶泉出露地表. 因此,经过一系列物理、 化学、 生物作用后,表层岩溶泉中无机氮几乎以NO3--N形式存在,并相对稳定.

可见,表层岩溶泉中三氮主要来源于与农业活动有关的非点源污染,4月因春季大量施用肥料使DIN质量浓度明显升高[图 3(a)、 3(e)]也证明了农业活动对表层岩溶泉的影响. 化肥、 人畜粪便等在土壤中经氨化、 硝化细菌作用转化为NO3-,下渗进入表层岩溶带后,在岩溶裂隙和溶隙中随扩散流运移(图 4).

4.2 地下河三氮的来源及转化运移 4.2.1 降水、 污水对地下河三氮的影响

大气降水入渗补给是地下水获得补充水量的主要方式,也是可溶态氮向下迁移的主要载体和驱动力[31, 32]. 研究区全年降水量1 273.36 mm,主要集中在3~9月,地下河对降水的响应与表层岩溶泉有很大不同. 降水时地下河流量迅速增加,降水结束后退峰过程相对较长,表现出典型的南方岩溶地下河管道流的水文动态变化特征(图 5),说明该区属管道与裂隙组合的地下水系统,地表与地下河通过管道直接相通,同时有深部裂隙和溶隙的调节[33],因此,地下河中扩散流和管道流并存,并以管道流为主.

地下河三氮质量浓度对降水三氮质量浓度有积极响应,除降水量突然增加的3、 9月外,二者具有相同的变化趋势[图 3(b)、 3(d)、 3(e)]. 同时,地下河三氮质量浓度对降水量也积极响应,具有雨季降低旱季升高的规律. 降水量与地下河DIN输出量的关系更加直接:将地下河每月DIN输出量与降水量用SPSS软件进行线性回归,设降水量为X,地下河DIN输出量为Y,那么,Y=0.056X+0.215,相关系数R=0.90,其相关性在α<0.01水平上达到显著正相关.

地下河、 污水三氮质量浓度也具有较一致变化趋势. 首先,二者都以NH4+-N为主要赋存形式; 其次,都对降水敏感响应[图 3(b)、 3(c)、 3(e)],从2014年8月至2015年1月二者的NH4+-N几乎同步变化. 污水和地下河中NO3--N和NO2--N都具有旱季质量浓度低、 雨季质量浓度高的规律.

因此,以管道流为主的地下河水与降水、 污水联系紧密,降水、 污水均可能通过管道与地下河直接相通. 与表层岩溶泉不同,污水可能成为地下河三氮的来源.

4.2.2 枯水期地下河三氮来源与转化运移

岩溶区降水进入地下河的路径主要有:表面坡流、 壤中流、 表层岩溶带水流、 落水洞水流、 深溶隙水流、 深裂隙水流[13, 34],其中,表面坡流、 壤中流通过落水洞直接进入,表层岩溶带水流部分经落水洞进入,部分经深溶隙和深裂隙进入. 因此,降水进入地下河的路径可以合并为:落水洞(管道流)、 深溶隙和深裂隙(扩散流),具体路径因降水的特征不同而有所差异,被降水驱动的污染物的运移路径也因此而不同. 可见,不同降水条件下,进入地下河的氮源及其路径和转化方式也不同.

2014年2月、 12月和2015年1月,月降水量均小于20 mm,未产生10 mm以上降水,日均径流量仅为55L ·s-1,且3个月均无明显变化(图 5). 这期间降水少,缺乏运移非点源污染物的驱动力和载体,大气、 地表、 土壤中含氮物质不能随降水进入地下河,生活、 工业污水等点污染源成为该段时间地下河三氮的主要来源[35, 36, 37]. 3个月中因没有降水的“稀释作用”,地下河DIN平均质量浓度高达10.35 mg ·L-1. 同时,低降水导致的低DO浓度使微生物的硝化作用被抑制,使NH4+-N占DIN的85.25%~94.70%,NH4+-N/NO3--N全年最高[图 3(f)]. 枯水期地下河水和污水三氮组成相似,并且二者的NO3--N质量浓度几乎一致(均小于1 mg ·L-1)[图 3(b)、 3(c)、 3(e)],进一步证明枯水期地下河三氮主要来源于点源污染物: 生活污水、 工业废水. 其运移路径主要有两种:一是通过管道引入落水洞直接排入地下河,这是主要路径. 二是排污管、 渠泄漏后经其下裂隙、 溶隙渗入地下河(图 4),三氮经这两种路径输入地下河的位置相对固定、 污染物量较稳定.

4.2.3 平水期、 丰水期地下河三氮来源与转化运移

研究期间,降水主要集中在3~9月,但10、 11月仍然有日降水量大于10 mm的降水,使地下河流量发生明显变化. 因此,不再将平水期和丰水期进行区别分析. 这期间地下河DIN输出量远远大于枯水期(图 5),说明有更多的三氮输入,其三氮组成也与枯水期很不相同. 随着降水的增加,NH4+-N质量浓度和占DIN比例大大降低,NO3--N却显著升高[图 3(b)],可能因为:第一,旱季存储有高NO3--N的表层岩溶带水在降水驱动下,经深溶隙、 深裂隙进入地下河[26],现场观察到降水后洞穴中裂隙水明显增加也证明了这个假设. 同时,降水产生的高NO3--N低NH4+-N的壤中流、 表层岩溶带水流经落水洞进入地下河. 第二,雨水带入的溶解氧使硝化细菌在有利条件下将地下水一部分NH4+-N转化为NO3--N. 其中间产物NO2--N的增加就是这一过程的证明[图 3(b)、 3(e)].

9月120 mm的大暴雨后地下河DIN质量浓度迅速上升,出现了大流量下的高DIN值[图 3(b)、 3(e)],这显然与降水的“稀释效应”不符,说明有大量三氮进入地下河,并以NO3--N为主要形式存在. 这与Tucker在美国Kentucky岩溶区研究结果一致[38]. 当月深裂隙水流DIN质量浓度为全年最低的7.93 mg ·L-1,体现出降水的“稀释效应”,明显低于地下河的9.35 mg ·L-1,说明高DIN质量浓度不是由表层岩溶带经深裂隙、 溶隙下渗产生. 现场观察表明,大暴雨后流域内积水现象突出,排水沟逸出、 大量露天堆放牲畜粪便被冲走,一些从未被淹没积累大量氮素的地方被淹,形成的表面坡流具有较高质量浓度的DIN,通过落水洞直接补给地下河,造成了地下河的高DIN质量浓度. 大量表面坡流引入的新鲜氧气更充分,加剧了硝化作用,出现NO3--N质量浓度最高值和NH4+-N质量浓度最低值[图 3(b)、 3(e)]. 9月大暴雨将地表储存的大量氮素清洗后造成了其后的10月地下河DIN质量浓度和输出量的最低值(图 5).

因此,平水期、 丰水期三氮来源和输入转化方式有:第一,与旱季时相同的生活、 工业废水点源污染通过相对固定的落水洞、 岩溶裂隙、 溶隙进入地下河; 第二,降水使表层岩溶带暂时饱和,高硝氮、 低氨氮的表层岩溶带水经深裂隙、 深溶隙(扩散流)继续向下运移进入地下河[25, 26]; 第三,降水形成的表面坡流、 壤中流、 表层岩溶带水流将地表氮素溶解后经落水洞(管道流)进入地下河(图 4). 可见,该时期地下水三氮来源于点源和非点源污染,并发生强烈硝化作用.

4.3 地下河中点源非点源氮的输入量估算

研究区地下河中DIN主要来源于四方面. 一是生活、 工业废水的点源污染; 二是大气氮沉降; 三是农业使用的化肥、 农肥和养殖业产生的牲畜粪便; 四是土壤有机氮的分解. 后三者均以非点源形式进入地下河. 降水是非点源污染物的载体和重要驱动力,但点源污染物不因降水的变化而不同,排放量相对稳定,成为枯水季节水质污染的主要来源[35],可根据水文特征由地下河基流推求[36]. 根据研究区地下河水文特征,用水平线分割法将流量稳定的枯水期1、 2、 12月流量分割为来自深层地下水补给的基流[39],用公式(5)[37]求得地下河稳定的月点源污染输入量. 根据公式(6)计算出研究年56.05 kg ·(hm2 ·a)-1的DIN通量中,点源污染、 非点源污染分别为14.06 kg ·(hm2 ·a)-1、 41.99 kg ·(hm2 ·a)-1,贡献率分别为25.08%、 74.92%,丰水期和平水期以非点源污染为主(图 5).

式中,Wg月为月点污染源输入DIN量(kg),C枯月为1、 2、 12月DIN质量浓度(mg ·L-1),C枯月为1、 2、 12月地下河流量(m3).
式中,WT为地下河DIN输出量(kg); Wu为非点源污染物输入DIN量(kg).

5 结论

(1)在人类活动影响下,表层岩溶泉和地下河均含有较高质量浓度的氮化合物,但其组成形式有很大差异. 表层岩溶泉以NO3--N为主要赋存形式,占DIN比均在98.15%以上,并相对稳定; 地下河旱季以NH4+-N为主要赋存形式,占DIN比78.85%以上,而雨季以NO3--N为主.

(2)表层岩溶泉DIN主要来源于农业、 养殖业的肥料、 人畜粪便等非点源污染. 污染物通过氨化作用、 植物吸收、 硝化作用、 吸附作用等一系列物理、 化学、 生物作用后经裂隙、 溶隙,随扩散流运移并出露地表. 地下河水DIN随扩散流和管道流运移,并以管道流为主. 其来源和运移方式因季节不同而有所差异. 枯水期主要是来源于点源的工业、 生活污染物经落水洞、 裂隙、 溶隙进入; 平水期、 丰水期DIN来源于点源污染和非点源污染,并以非点源污染物为主,包括表层岩溶带饱和后经深裂隙和后继续下渗进入地下河(扩散流)和降雨产生的地表径流冲刷污染物、 壤中流、 表层岩溶带流从落水洞进入地下河(管道流)两种形式.

(3)全年流域DIN经地下河流出73 365.23 kg,输出通量为56.05 kg ·(hm2 ·a)-1,其中NH4+-N、 NO3--N分别占总输出量的46.03%、 52.51%; 点源污染物、 非点源污染物的贡献率分别为25.08%、 74.92%.

参考文献
[1] 阮晓红, 王超, 朱亮, 等. 氮在饱和土壤层中迁移转化特征研究[J]. 河海大学学报, 1996, 24 (2): 51-55.
[2] Jansons V, Busmanis P, Dzalbe I, et al. Catchment and drainage field nitrogen balances and nitrogen loss in three agriculturally influenced Latvian watersheds[J]. European Journal of Agronomy, 2003, 20 (1-2): 173-179.
[3] Flipo N, Even S, Poulin M, et al. Modeling nitrate fluxes at the catchment scale using the integrated tool CAWAQS[J]. Science of the Total Environment, 2007, 375 (1-3): 69-79.
[4] Salazar O, Wesström I, Youssef M A, et al. Evaluation of the DRAINMOD-N II model for predicting nitrogen losses in a loamy sand under cultivation in south-east Sweden[J]. Agricultural Water Management, 2009, 96 (2): 267-281.
[5] 吴迪. 红壤丘陵区地下水"三氮"迁移转化数值模拟研究[D]. 阜新: 辽宁工程技术大学, 2011.
[6] 李想. 千烟洲试验区地下水中"三氮"的迁移转化规律及数值模拟研究[D]. 阜新: 辽宁工程技术大学, 2011.
[7] 周美英. 地下水中"三氮"运移的随机模拟研究——以沈阳市为例[D]. 北京: 中国地质大学, 2012.
[8] Hinkle S R, Tesoriero A J. Nitrogen speciation and trends, and prediction of denitrification extent, in shallow US groundwater[J]. Journal of Hydrology, 2013, 509 : 343-353.
[9] Zhang Q Q, Sun J C, Liu J T, et al. Driving mechanism and sources of groundwater nitrate contamination in the rapidly urbanized region of South China[J]. Journal of Contaminant Hydrology, 2015, 182 : 221-230.
[10] 刘丽雅, 何江涛, 王俊杰. 浑河傍河区地下水氮污染来源贡献研究[J]. 岩土工程技术, 2015, 29 (2): 59-64.
[11] Aquilina L, Ladouche B, Dörfliger N. Water storage and transfer in the epikarst of karstic systems during high flow periods[J]. Journal of Hydrology, 2006, 327 (3-4): 472-485.
[12] Gaouzi F Z J E, Sebilo M, Ribstein P, et al. Using δ15N and δ18 O values to identify sources of nitrate in karstic springs in the Paris basin (France)[J]. Applied Geochemistry, 2013, 35 : 230-243.
[13] Ford D, Williams P. Karst hydrogeology and geomorphology[M]. Chichester: John Wiley & Sons Ltd, 2007. 104-208.
[14] Atkinson T C. Diffuse flow and conduit flow in limestone terrain in the Mendip Hills, Somerset (Great Britain)[J]. Journal of Hydrology, 1977, 35 (1-2): 93-110.
[15] Jacobson R L, Langmuir D. Controls on the quality variations of some carbonate spring waters[J]. Journal of Hydrology, 1974, 23 (3-4): 247-265.
[16] Fretwell B A, Burgess W G, Barker J A, et al. Redistribution of contaminants by a fluctuating water table in a micro-porous, double-porosity aquifer: field observations and model simulations[J]. Journal of Contaminant Hydrology, 2005, 78 (1-2): 27-52.
[17] Mahler B J, Garner B D, Musgrove M, et al. Recent (2003-05) water quality of Barton Springs, Austin, Texas, with emphasis on factors affecting variability[R]. Virginia: United States Geological Survey, 2006.
[18] Morasch B. Occurrence and dynamics of micropollutants in a karst aquifer[J]. Environmental Pollution, 2013, 173 : 133-137.
[19] Mahler B, Massei N. Anthropogenic contaminants as tracers in an urbanizing karst aquifer[J]. Journal of Contaminant Hydrology, 2007, 91 (1-2): 81-106.
[20] Vesper D J, Loop C M, White W B. Contaminant transport in karst aquifers[J]. Speleogenesis and Evolution of Karst Aquifers, 2003, 1 (2): 1-11.
[21] Katz B G. Nitrate contamination in Karst groundwater[A]. In: Culver D, White W (eds.). Encyclopedia of Caves (2nd ed.)[M]. Amsterdam, The Netherlands: Elsevier Science, 2004. 415-418.
[22] Einsiedl F, Radke M, Maloszewski P. Occurrence and transport of pharmaceuticals in a karst groundwater system affected by domestic wastewater treatment plants[J]. Journal of Contaminant Hydrology, 2010, 117 (1-4): 26-36.
[23] Magal E, Arbel Y, Caspi S, et al. Determination of pollution and recovery time of karst springs, an example from a carbonate aquifer in Israel[J]. Journal of Contaminant Hydrology, 2013, 145 : 26-36.
[24] 林明珠. 不同土地覆被下表层岩溶带水文地球化学特征及其碳汇效应研究[D]. 重庆: 西南大学, 2012.
[25] Mahler B J, Valdes D, Musgrove M, et al. Nutrient dynamics as indicators of karst processes: comparison of the Chalk aquifer (Normandy, France) and the Edwards aquifer (Texas, U.S.A.)[J]. Journal of Contaminant Hydrology, 2008, 98 (1-2): 36-49.
[26] 杨平恒, 袁道先, 任幼蓉, 等. 川东平行岭谷区典型岩溶含水系统中NO3-的存储和运移[J]. 环境科学, 2012, 33 (9): 3124-3131.
[27] Xu Z X, Hu B X, Davis H, et al. Simulating long term nitrate-N contamination processes in the Woodville Karst Plain using CFPv2 with UMT3D[J]. Journal of Hydrology, 2015, 524 : 72-88.
[28] 蓝家程, 杨平恒, 任坤, 等. 重庆老龙洞地下河流域氮、磷及微生物污染调查研究[J]. 环境科学, 2014, 35 (4): 1297-1303.
[29] 谢运球, 玉宏, 裴建国, 等. 桂林寨底地下河农业系统氮流失估算[J]. 农业环境科学学报, 2015, 34 (10): 1979-1984.
[30] 许可. 地下水中三氮污染物研究综述与展望[J]. 工程勘察, 2012, 40 (8): 32-35.
[31] 杨维, 郭毓, 王晓华, 等. 氮素在包气带与饱水层迁移转化的实验研究[J]. 环境科学研究, 2008, 21 (3): 69-75.
[32] 王开然, 郭芳, 姜光辉, 等. 桂林峰林平原区岩溶含水层氮污染空间分布特征[J]. 环境科学研究, 2013, 26 (3): 281-286.
[33] 任坤, 杨平恒, 江泽利, 等. 降雨期间岩溶城镇区地下河水重金属变化特征及来源解析[J]. 环境科学, 2015, 36 (4): 1270-1276.
[34] 袁道光, 蔡桂鸿. 岩溶环境学[M]. 重庆: 重庆出版社, 1988.
[35] 李怀恩. 估算非点源污染负荷的平均浓度法及其应用[J]. 环境科学学报, 2000, 20 (4): 397-400.
[36] 陈友媛, 惠二青, 金春姬, 等. 非点源污染负荷的水文估算方法[J]. 环境科学研究, 2003, 16 (1): 10-13.
[37] 蔡明, 李怀恩, 刘晓军. 非点源污染负荷估算方法研究[J]. 人民黄河, 2007, 29 (7): 36-37, 39.
[38] Tucker N L. Nonpoint agricultural pollution in a karst aquifer: lost river groundwater drainage basin, warren county, kentucky[D]. Bowling Green: Western Kentucky University, 1982.
[39] 徐磊磊, 刘敬林, 金昌杰, 等. 水文过程的基流分割方法研究进展[J]. 应用生态学报, 2011, 22 (11): 3073-3080.