环境科学学报  2017, Vol. 37 Issue (2): 566-575
黄土高塬沟壑区氮素流失的时空变异特征    [PDF全文]
林雪青, 李志 , 李志恒, 韦祎, 季翔    
西北农林科技大学资源环境学院, 杨凌 712100
摘要: 为了探讨黄土高塬沟壑区水体氮污染的时空变化情况,选取黑河流域(泾河支流)为研究区域,测定2013-2014年枯水期和汛期流域地表水和地下水中主要离子及NO3--N和NH4+-N的浓度并进行分析.结果表明,黑河流域枯水期水化学主要为Na+K-Cl-SO4型,汛期主要为Ca-Mg-HCO3型.枯水期及汛期阳离子均主要为Na+,阴离子在枯水期主要为SO42-而汛期则转变为HCO3-.汛期NO3--N浓度普遍大于枯水期,平均值分别为2.37和1.63 mg·L-1;且空间分布不均衡,地表水中的浓度为:上游(1.35 mg·L-1)>中游(1.05 mg·L-1)>下游(0.93 mg·L-1),而地下水中的浓度为:下游(3.84 mg·L-1)>中游(2.54 mg·L-1)>上游(2.35 mg·L-1).NH4+-N在时间分布上没有明显的规律,汛期及枯水期变化不大,空间分布特征与NO3--N类似,但其整体浓度较低,在0.11 mg·L-1左右波动,较为稳定均且未超过IV类水标准.水体中NO3--N不仅来自于农田氮肥的施入等人类活动,还可能来自于酸性降雨.地表水的NO3--N污染程度存在空间差异,上游污染程度大于中、下游,而超过70%的地下水水质属于良好,对当地饮用水安全暂不造成威胁.
关键词: 水化学     硝态氮     氨态氮     地表水     地下水     时空变异    
Spatio-temporal variations of nitrogen loss in the tableland-gully region of the Loess Plateau
LIN Xueqing, LI Zhi , LI Zhiheng, WEI Yi, JI Xiang    
College of Natural Resources and Environment, Northwest A & F University, Yangling 712100
Received 6 May 2016; received in revised from 5 Jul 2016; accepted 22 Jul 2016
Supported by the National Natural Science Foundation of China (No.51179161), the Science and Technology Research and Development Program of Shaanxi Province (No.2013KJXX-18) and the the Fundamental Research Funds for Central Universities (No.2452015105)
Biography: LIN Xueqing (1992-), female, E-mail:18302941034@163.com
*Corresponding author: LI Zhi, E-mail:lizhibox@126.com
Abstract: To investigate the spatiotemporal variations of the nitrogen pollution in the tableland-gully region of the Loess Plateau, the surface water and groundwater in the Heihe watershed (tributary of the Jing River catchment) have been sampled in the dry and flood seasons for the period of 2013-2014 to determine the concentration of NO3--N and NH4+-N. The major water type is Na+K-Cl-SO4 during the dry seasons and Ca-Mg-HCO3 during the flood seasons. Na+ is the major cation for both dry and flood season, while the major anions are SO42- in the dry seasons and HCO3- in the flood seasons. The loss of nitrogen is mainly in the form of NO3--N, and the concentration of NO3--N in the flood seasons is higher than that in the dry seasons, which is respectively 2.37 and 1.63 mg·L-1. Spatial variations exist in the concentration of NO3--N. For surface water, it is the highest in the upstream (1.35 mg·L-1), intermediate in the midstream (1.05 mg·L-1) and the lowest in the downstream (0.93 mg·L-1); however, for groundwater, it is the lowest in the upstream (2.35 mg·L-1), intermediate in the midstream (2.54 mg·L-1), and the highest in the downstream (3.84 mg·L-1). The pollution from NO3--N was further quantified by comparing the concentration with environmental standards. For surface water, the pollution in the upstream was greater than that in the midstream and downstream; for groundwater, more than 70% of the groundwater samples are good enough to be drinking water. The NO3--N is possibly from fertilizer application and acid precipitation. The concentration of NH4+-N is spatiotemporally invariant and it fluctuates at around 0.11 mg·L-1.
Key words: hydro-chemistry     nitrate nitrogen     ammonia nitrogen     surface water     groundwater     spatiotemporal variations    
1 引言(Introduction)

随着点源污染逐渐得到控制, 非点源污染问题日益加重, 成为水环境污染的重要污染源, 而农田氮素流失又是非点源污染的重要表现形式.自1950年起, 我国农业地区氮肥用量急剧增加, 大量氮肥施入土壤加之不合理的管理造成其有效利用率仅为30%~40%, 使得剩余的氮肥主要以NO3--N形式随降雨径流流失, 进而造成地表水污染或在土壤深剖面中累积淋失至地下水.

由于氮素随降雨径流从农田向河流、湖泊、水库等水体输移, 不仅降低了土壤肥力和化肥效率, 而且造成藻类大量繁殖、富集, 出现富营养化现象, 大量消耗水中的溶解氧并危及地表水体的环境安全(金洁等, 2005Tian et al., 2007Zhao et al., 2009).通过农田输入湖泊的氮量占输入湖泊总氮量的7%~35%(Ranjith et al., 2006), 降雨径流及农业排水携带的硝态氮是造成地表水污染的主要原因.我国已有63.3%的湖泊水体达到富营养化(王晓燕等, 2004), 两大淡水湖-巢湖和太湖都存在严重的水体富营养化现象.因此, 严格控制氮素的非点源污染在水环境保护方面具有重要的意义.

地下水资源在北方干旱半干旱地区农牧业供水等中占有重要地位, 地下水安全直接影响到我国人民的生活质量(Kumazawa et al., 2002).20世纪60年代, 我国就有关于地下水NO3--N污染的报道, 农业地区大量氮肥的施入, 使得土壤深剖面含有大量NO3--N, 约有20%~50%随土壤水分向下淋溶至地下水, 对地下水造成污染.世界卫生组织规定地下水饮用水的NO3--N浓度标准为10 mg·L-1, 我国大部分地区都超标甚至严重超标.地下水中过量的硝酸盐会导致高铁血红蛋白症及消化系统癌症的发生, 对人类健康造成极大危害(Arauzo et al., 2015Laftouhi et al., 2003).并且土壤和地下水中的氮素污染还可能带来农产品氮浓度的增加, 从而导致粮食、蔬菜、瓜果等氮污染, 并经由饮食途经进一步危害人体健康(郭开秀等, 2011).

我国对于流域氮素流失的研究多集中在南方温暖湿润、降雨量较高的地区, 包括巢湖、太湖及九龙江, 而对于较为缺水干旱的黄土高原地区的研究甚少.黄土高原地区严重缺水, 降雨是其水资源的主要补给来源, 并且随着农村地区农业发展, 氮肥施入量逐渐增加, 氮素随降雨和土壤淋溶流失越发严重(Ran et al., 2006陆海明等, 2008), 当前的污染状况亟需进行评估.因此, 本文以黄土高塬沟壑区的典型代表黑河流域为研究对象, 对2013-2014年枯水期和汛期地表水及地下水中NO3--N和NH4+-N浓度进行分析, 探究黑河流域氮素的时空变化规律, 以期为该流域氮素污染防治提供科学依据并丰富流域氮素流失的相关研究.

2 材料与方法(Materials and methods) 2.1 研究区域概况

黑河流域位于黄土高原中南部(图 1), 陕西省长武县与甘肃省华亭县之间, 海拔1000~2500 m.流域面积1506 hm2, 呈东西走向, 形状狭长, 地势西高东低, 地貌属黄土高塬沟壑区, 包括塬面和沟壑两大单元, 是黄土高塬沟壑区的典型代表.属暖温带半湿润大陆性季风气候, 夏季炎热多雨, 冬季寒冷干燥, 年平均气温9.4 ℃, 年平均降水量522~608 mm, 多集中在7-9月.绝大部分生产可利用地都分布在塬面和梁顶上, 土壤主要是黑垆土和黄墡土.土地利用方式主要为耕地、草地和林地, 分别占流域面积的58%、36%和5%(李志等, 20082010).而农业用地的化肥施用量很高, 以长武县为例, 麦地年均化肥施用量为750~1500 kg·hm-2 (曹群虎等, 2009), 而苹果地施肥量更大为4000~10000 kg·hm-2 (朱亚莉等, 2015), 远超过国际安全施肥上限225 kg·hm-2(马俊杰等, 2003).同时, 该研究区水土流失严重, 年径流量57.7~35.4 m3·s-1, 侵蚀模数174.9~696.7 t·hm-2·a-1(胡宏祥等, 2009).施肥量大和水土流失严重的共同影响下会导致养分流失严重(Maila et al., 2004), 并且有关该地区的水体氮素污染情况尚不明确, 因此, 非常有必要研究该流域氮素时空分布情况.

图 1 黑河流域采样点分布 Fig. 1 Sampling points in the Heihe basin
2.2 样品采集及测定

为充分考虑降水径流关系对氮素淋失的影响, 分析了长期监测的降水、蒸散和径流特征(图 2), 发现6月蒸散量最大而降雨量较小, 导致径流量最小, 是枯水期的典型代表;7-9月降雨量最大, 占全年55%以上, 而蒸散量较小, 致使径流量最大, 具有典型的汛期特征.因此, 采样时间确定为6月和9月, 分别代表枯水期和汛期.4次采样的具体时间分别为2013年6月29日、2013年9月1日、2014年6月29日、2014年9月20日.

图 2 黑河流域水文气象特征 Fig. 2 Hydrometeorological characteristics in the Heihe basin

尽管样品采集选择了典型的月份, 日期的选择仍可能存在不确定性, 但同期的稳定同位素分析可为采样的代表性提供支持.2年的同位素组成分析表明, 地表水同位素受降水的影响比较明显, 而地下水同位素组成较为稳定;同位素的空间分布上, 枯水期和汛期, 地表水均在中、下游富集, 而地下水从上到下呈贫化状态(宋献芳等, 2009).这些稳定的同位素组成特征表明, 黑河流域具有稳定的水循环特征, 因此, 采样能反映该区的水文特征, 具有一定的代表性.

沿黑河流域均匀采集地表水(河水和沟道径流)和地下水(井水和泉水)共160个水样, 河水即黑河河水, 沟道径流即在河道两侧的径流, 井水为河道两侧的井水, 泉水为基岩裂隙水.每次采样在同一位置或附近采集, 即1个样点采集4次水样.用GPS测定各点的经纬度及高程, 同时用pH计和电导率仪现场测定水样的pH值、电导率(EC)及水温(T).采集的水样置于清洗过的250 mL的聚乙烯瓶中, 4 ℃保存, 上机前先过0.45 μm的微孔滤膜去除悬浮颗粒物, 用连续流动分析仪测定样品中硝态氮(NO3--N)和氨态氮(NH4+-N);离子色谱仪测定Cl-、SO42-;ICP测定K+、Ca2+、Na+、Mg2+;采用甲基橙指示剂稀硫酸滴定计算HCO3-.

3 结果与分析(Results and analysis) 3.1 水化学及电导率特征

为说明黑河流域地表水及地下水的基本水化学特征, 绘制Piper三线图(图 3)进行分析.以2014年枯水期(6月)及汛期(9月)黑河流域水样为例进行分析, 黑河流域水化学特征具有时间上的变化趋势, 且不同水体变化不一.枯水期地表水中各阴、阳离子的浓度均值略大于地下水, 是地下水的1.2倍;而汛期的浓度均值则是地表水小于地下水, 是地下水的0.7倍.该研究区域水体中离子组成变化较大, 枯水期水体中阳离子浓度均值呈现:Na+(73.9 mg·L-1)>Mg2+(29.4 mg·L-1)>Ca2+(24.7 mg·L-1)>K+(2.1 mg·L-1), 阴离子浓度平均值呈现:SO42-(75.6 mg·L-1)>HCO3-(31.2 mg·L-1)>Cl-(23.3 mg·L-1).汛期HCO3-浓度明显增大, 总的来说, 阳离子浓度均值呈现:Na+ (74.8 mg·L-1)>Ca2+(57.1 mg·L-1)>Mg2+(33.3 mg·L-1)>K+(2.8 mg·L-1), 阴离子浓度均值呈现:HCO3-(316.5 mg·L-1) >SO42-(92.9 mg·L-1)>Cl-(25.1 mg·L-1).根据Piper三线图(图 3)可知, 枯水期水化学主要为Na+K-Cl-SO4类型, 汛期水化学主要为Ca-Mg-HCO3类型.

图 3 黑河流域水样Piper三线图(a.枯水期;b.汛期) Fig. 3 The Piper diagram of Heihe basin

电导率(EC)通常表征水体中阴、阳离子的浓度(刘相超等, 2008).由电导率仪测定的结果可知, 所采集水样的EC在2013年枯水期及汛期分别为710和798 μS·cm-1, 2014年分别为725和791 μS·cm-1, 与2014年阴、阳离子浓度变化情况一致.该研究区域地下水中的EC值均高于300 μS·cm-1, 最高达2610 μS·cm-1, 说明受污染程度较大.地下水中EC均值(2013年枯水期及汛期分别为876、910 μS·cm-1, 2014年分别为803、1007 μS·cm-1)高于地表水(2013年枯水期及汛期分别为584、599 μS·cm-1, 2014年分别为726、755 μS·cm-1).这表明地下水更新速率小于地表水, 滞留时间大于地表水, 一旦受到污染其治理难度远大于地表水.

3.2 流域氮素流失时间变化 3.1.1 不同水文期变化特征

黑河流域水体中的氮素以NO3--N为主, 其浓度无论汛期还是枯水期均高于NH4+-N.对于NO3--N而言, 汛期水样中的平均浓度大于枯水期(表 1), 如2013年6、9月分别为1.89、2.13 mg·L-1, 而2014年6、9月分别为1.37、2.60 mg·L-1.但NH4+-N浓度在两个时期变化不大, 6月为0.02~0.44 mg·L-1, 而9月为0.02~0.21 mg·L-1, 整体波动较小.主要是由于NH4+-N带正电, 易被土壤胶体吸附, 随降雨径流流失或通过壤中流途径淋失风险较小, 对地表水及地下水造成的影响小.汛期和枯水期的差异主要受降水量年内分布影响, 黑河流域汛期降水量占全年降水量的55%以上, 而汛期离子浓度明显高于枯水期, 说明降水径流过程对NO3--N流失有明显的季节性影响(Chen et al., 2005Guo et al., 2004张亚丽等, 2014).

表 1 2013、2014年NO3--N和NH4+-N变化 Table 1 Change of nitrate、ammonia nitrogen in 2013、2014

不同时期不同类型水样中NO3--N浓度及变化情况不同.NO3--N浓度在井水中最大, 沟道径流中最小;NO3--N浓度在河水和沟道径流中呈现枯水期 < 汛期, 而在井水中则为枯水期 > 汛期, 泉水中的NO3--N变化不一.NH4+-N浓度在不同类型水样中变化不明显, 没有明显的增加或减少, 整体较稳定且浓度低, 均在0.11 mg·L-1左右波动.根据差异性分析可知(表 2), 地表水及地下水中NO3--N浓度普遍呈显著性差异.而对于NH4+-N而言, 在4个不同时期差异均不显著.

表 2 2013、2014年地表水及地下水NO3--N和NH4+-N差异性分析 Table 2 Variability analysis of nitrate、ammonia nitrogen in 2013、2014
3.1.2 年际变化

表 1可知, 2013年河水、沟道径流、泉水、井水中NO3--N浓度均值分别为1.43、0.69、1.89、3.25 mg·L-1, 2014年分别为1.45、1.17、1.36、5.31 mg·L-1, 随时间推移, 除泉水外其他3类水样的NO3--N浓度均增大.NH4+-N在该流域年际变化不明显, 浓度普遍较低, 在2013年分别为0.11、0.11、0.09、0.13 mg·L-1, 2014年分别为0.13、0.12、0.08、0.13 mg·L-1, 4种类型的水样彼此间差异及年际变化不明显.NO3--N浓度的年际变化可能与气候有显著的关系, 如2013年和2014年降水量分别为574.0和413.2 mm, 其中, 汛期降水量的比例分别为70%和47%, 2013年较高的全年和汛期降水量可能是NO3--N浓度偏低的主要原因, 在降雨初期土壤表层NO3--N随径流流失, 造成水体中浓度迅速升高, 而随着降雨的继续, 雨量较大的情况下会对水体中NO3--N造成稀释作用使其浓度降低.

3.2 氮素空间变化特征

按照经纬度及高程将研究区域整体划分为上、中、下游, 2013-2014年黑河流域由上游到下游NO3--N浓度变化趋势如图 4所示.2013年河水(图 4a)和沟道径流(图 4b)中NO3--N浓度沿上游至下游均有减小趋势, 相反泉水(图 4c)和井水(图 4d)中NO3--N浓度逐渐增大, 2014年的变化情况与2013年类似.NH4+-N浓度的变化趋势(图未给出)与NO3--N一致, 在整个研究区域浓度小, 波动不明显.2013-2014年NO3--N浓度最大值均出现在上游井水, 分别为32.85和17.10 mg·L-1, 均值分别为9.27和5.40 mg·L-1;最小值在上游沟道径流, 分别为0.14和0.02 mg·L-1, 均值分别为1.09和0.47 mg·L-1.NH4+-N浓度没有明显规律.

图 4 2013-2014年黑河流域NO3--N空间分布图 Fig. 4 Spatial variation in nitrate concentrations in 2013、2014 along Heihe basin

利用变异系数对黑河流域上、中、下游进行变异程度定量, 可知2013年各河段NO3--N浓度为中等变异(0.1 < CV < 1.0), 井水NO3--N浓度为在上、中游变异程度最大;而NH4+-N变异程度较低, 多属弱变异.2014年NO3--N多属于中等变异, 但上游河水和径流达到强变异(CV≥1.0);NH4+-N浓度变异程度最小值和最大值分别在下游和中游泉水, 分别为0.11和1.08.由此可知, 不同时期和不同空间位置下, 水样的NO3--N和NH4+-N变异程度不一, 2014年的变异程度较2013年强.变异程度越强, 表明水体中NO3--N浓度受空间位置影响更大;而对于低浓度低变异程度的NH4+-N来说, 空间位置的影响可以忽略.造成空间差异的主要原因是黄土高塬沟壑区塬高沟深, 黑河流域上游塬面少沟道多, 地势起伏大, 水对土壤中的氮素的冲刷能力更强, 造成黑河流域氮素流失空间不均衡.对2013-2014年上、中、下游地表水及地下水中NO3--N浓度进行单因素方差分析, 结果(表 3)表明, 地表水中NO3--N浓度无空间差异, 但地下水上、中游与下游间差异显著.而不同河段的均值表明, 地表水中NO3--N浓度呈现:上游(2.84 mg·L-1)>中游(2.17 mg·L-1)>下游(1.67 mg·L-1), 而地下水呈现:下游(7.68 mg·L-1)>中游(4.07 mg·L-1)>上游(3.10 mg·L-1).

表 3 2013、2014年地表水及地下水NO3--N浓度均值空间变化 Table 3 Spatial variation of nitrate nitrogen in 2013、2014
3.3 主要离子来源分析

利用阴、阳离子之间的相关性分析, 可以大致判断各离子来源及其一致性.人类活动产物的特征主要富含K、Ca、S、Cl和N, 其中, K、Ca、S、Cl同时又是岩石/土壤风化作用的产物(Stallard et al., 1981).因此, N在水体的主要存在形式NO3--N及其浓度变化则能指示人为活动对水体化学特征的影响.根据表 4可知, 无论枯水期还是汛期NH4+-N与其他离子之间相关性均不显著, 受外界影响较小.汛期流域水体中NO3--N浓度明显增大, 其受人为因素影响较大, 农业化肥和牲畜粪便会随着降雨及地下径流污染地表水和地下水, 成为该流域水体氮素主要来源(Nolan et al., 2000).同样, NO3--N与HCO3-也具有显著相关性, 反映二者有共同来源, 表明水体中NO3--N不仅来自于农田氮肥的施入等人类活动, 还可能来自于酸性降雨.K+、Na+、Mg2+都显示出与SO42-具有显著的相关性, 这表明阳离子可能主要来源于硫酸盐.Cl-与SO42-也表现出显著的相关性, 这表明Cl-对水体的酸化和水质变化有一定影响(张芳等, 2013).SO42-与Ca2+也表现出显著的相关性, 表明SO42-应该主要来源于石膏(CaSO4·H2O)等蒸发岩的溶解(耿金等, 2013).同时, 该地区上游洗煤严重, 因此, 水体中Cl-和SO42-也可能来自于工业污染源.汛期各离子之间的相关系数明显显著于枯水期, 表明这些离子相互影响或有相同的来源.

表 4 黑河流域离子相关性分析 Table 4 Correlation analysis of ions of Heihe basin
3.4 氮素污染状况评估

按照《中华人民共和国地表水环境质量标准》(GB-3838-2002)中总氮的标准, 将采集的2013-2014年河水和径流水样分类(图 5), 可知黑河流域上游地表水水样氮污染程度较为均匀, 每类水质所占比率较为接近.中、下游主要达到Ⅳ类水标准, 分别约占37%和34%.上游水样达到Ⅴ类及劣Ⅴ类水标准所占比率较中、下游大, 约为22%.说明黑河流域由上游至下游NO3--N流失逐渐稀释减小, 使得到达下游时减缓了NO3--N污染.根据国家地下水硝、氨态氮质量标准(GB/T 14848-93)可知, 地下水中NO3--N浓度绝大部分未超过国家饮用水标准10 mg·L-1, 超过70%的水质属于良好(≤5 mg·L-1), 对当地饮用水安全暂不会造成威胁.NH4+-N浓度达到Ⅲ类水标准所占比例最高(表 5), 在2013、2014年分别为96%和77%, 并且所有水样的NH4+-N均没有超过0.5 mg·L-1,即未超过Ⅳ类水标准.

表 5 2013、2014年地下水NO3--N、NH4+-N频率分布数据统计表 Table 5 Nitrate Ammonia nitrogen concentrations and frequencies in 2013、2014

图 5 上、中、下游地表水水质变化 Fig. 5 Variations of surface water along heihe watershed
4 讨论(Discussion)

黑河流域汛期前后, 不同水体中阳离子均以Na+为主, 阴离子以HCO3-的浓度有明显增大, 其他离子变化程度较小.水体中化学成分的主要来源有大气沉降、岩石风化、蒸发-结晶及人类活动影响.Piper三线图不但反映了河水的化学组成, 而且还可以区分不同风化源区的物质组成, 从而辨别其控制端元(叶宏萌等, 2010).纯碳酸岩的风化物质以HCO3-为主, 以蒸发岩风化产物为主的数据点落在Cl--SO42-一端, 阳离子受碳酸岩风化产物影响落在Ca2+-Mg2+, 靠近Ca2+端, 白云岩风化产物应落在Ca2+-Mg2+线中间, 蒸发岩矿物风化产物落在Na++K+(韩贵林等, 2000侯韶华等, 2009).根据Piper三线图可初步判断黑河流域水体主要受碳酸岩和蒸发岩矿物风化产物的影响, 这与我国河流主要受碳酸岩和蒸发岩情况一致(Ming et al., 1982).枯水期水体阳离子中Na+占优势, 这主要是由于河水的蒸发浓缩作用造成.汛期Ca2+和Na+为优势离子, 地表水中Na+浓度较在枯水期降低, 主要是由于降雨造成稀释, 而地下水中Ca2+和Na+浓度明显增大.阴离子中造成HCO3-浓度增大的现象可能是因为富含HCO3-的水体对其进行补给, 其次有CO2进入水体.对于作为地下水的基岩裂隙泉水, 人为因素造成的影响可能性很小, 主要是由于自然因素造成的, 酸性降雨对其影响较大.由于水岩相互作用和蒸发浓缩作用使得不同时期水体中离子组成不一, 水化学类型也随之变化.枯水期水化学类型主要为Na+K-Cl-SO4型, 汛期主要为Ca-Mg-HCO3型, 雨水的水化学类型为Ca-HCO3, 流域水化学特征受降雨影响较大.同时, 人类活动如施肥、生活污水等造成化学离子的输入及降雨对水体的离子组成及水化学类型也有很大影响, 这与该地区近十几年来耕地大幅转变为苹果园的土地利用方式变化背景一致.

黑河流域4种水样的NO3--N平均浓度存在地下水大于地表水的规律, 如井水(4.43 mg·L-1)>泉水(1.63 mg·L-1)>河水(1.44 mg·L-1)>沟道径流(0.93 mg·L-1), 在枯水期含有高浓度NO3--N的浅层地下水对地表水的补给将会对地表水环境造成威胁.NH4+-N浓度普遍较小,在0.11 mg·L-1左右, 可知可溶态氮素主要以NO3--N的形式流失进入水体造成污染, 与以往的研究结果一致(陆海明等, 2008Sunndar et al., 2009).EC值变化规律与NO3--N一致, 同样为地下水 > 地表水, 金赞芳等(2004)的研究表明, 未受污染的地下水EC<300 μS·cm-1, 而EC值大于1000 μS·cm-1地区的NO3--N均值也大于10 mg·L-1.黑河流域中水体的EC值多介于300~1000 μS·cm-1, 虽然NO3--N浓度并未超过10 mg·L-1, 但通过EC值可知, 该区域的水环境安全需要高度重视.

本研究区降雨量低且较为集中, 氮素流失时期较南方湿润区短(Lu et al., 2007Cao et al., 2006韩梅等, 2014), 具有明显的季节性差异.汛期地表水(河水、沟道径流)中的NO3--N浓度普遍大于枯水期, 河水和沟道径流中NO3--N浓度分别增加了3.2倍和2.5倍.张继宗等(2009)在太湖流域的研究同样反映出汛期氮素流失较非汛期大.降雨加剧了氮素的流失程度, 使地表中的氮肥或其他形式的氮素随降雨流失到河水或进入径流(Haro et al., 2010Peter et al., 2003).有研究表明, 降雨总量增加1%, 土壤侵蚀就会增加1.7%(曹崇法等, 2001), 降雨多发季节雨水较易形成地表径流, 造成严重的水土和水肥流失, 从而导致水体中营养物质浓度急剧升高(Albert et al., 1978).

汛期井水中NO3--N浓度普遍较枯水期减小, 2013年由3.56 mg·L-1降至3.03 mg·L-1, 2014年由5.40 mg·L-1降至5.24 mg·L-1, 可能是由于受到了浓度较低的水源补给, 如降水通过快速通道直接补给(程立平等, 2012Sugimoto et al., 2003), 或者是含有低浓度NO3--N的河水通过土壤对地下水的补给导致其浓度降低.泉水没有明显的规律, 主要是由于泉水主要来自承压水层, 受外界影响较小, 且采样点普遍远离农田, 不受施肥等人为因素影响.

2014年较2013年河水、沟道径流、井水中的NO3--N均有不同程度的增大, 而泉水则减小, 采样地点较为一致, 自然因素如地貌、地形、植被、土壤性质等方面变化不大, 主要受降雨及人为因素影响.2013年和2014年降雨量分别为574和413.2 mm, 而汛期降水量分别占全年的70%和47%.短期内人类活动导致的氮素输入没有大的差异, 降水径流量直接决定了水体中的NO3--N浓度, 由于2013年全年和汛期降雨量都高于2014年, 地表径流对地表水中氮素的稀释作用更大, 导致2013年汛期地表水硝态氮浓度偏低.泉水为基岩裂隙泉水, 受人为因素影响小, 主要是自然因素影响NO3--N浓度, 这一部分需要深入研究.

4个不同时期地下水中的NO3--N浓度从上游至下游均有逐渐增加的趋势, 这与王庆锁等(2011)在密云水库流域的研究结果一致.影响地下水中硝态氮空间分布特征的因素主要有2类, 氮素来源及其所处的周围环境.造成黑河流域地下水这种变化主要是由于水体在下游的富集作用及下游土壤中高浓度的NO3--N淋溶造成的, 根据离子相关性分析的结果可知, 农业化肥、牲畜粪便及酸性降雨是该地区NO3--N主要来源.该研究区域主要地形为塬面和坡地, 在中、下游地区由于高程较低, 地势较为平缓, 塬区主要集中在中、下游, 耕地面积较上游大, 受人类活动影响大, 施入土壤中的农业化肥较多, 产流系数较低, 再加之NO3--N本身不易被土壤胶体吸附的特性, 造成较多NO3--N在降雨条件下随水分向下运移并在深层土壤剖面累积进而淋失进入地下水(Di et al., 2002Hu et al., 2000).而在上游地区, 农业耕地范围较小, 主要地形为自然坡地, 氮肥施入量远远低于下游, 加之含水层渗透性一般较下游好, 地下水径流较快速, 大部分氮素被携带向下游运移, 因此,在该区域含量较低(王佳因等, 2013).地表水则相反, 自上游至下游有逐渐减小的趋势, 影响地表水中氮素浓度的主要因素为降雨-地表径流, 土地类型、植被覆盖、施肥量等都会影响地表径流中氮素浓度.该研究区域坡地所占比例远高于塬面, 坡地植被类型普遍为刺槐及荒草地且上游坡度较下游地区坡度高, 造成在降雨条件下对土壤表层的氮素冲刷作用较强, 养分流失较为严重, 也使得其随径流进入地表水中并造成污染.并且在上游区域有煤矿厂, 洗煤对河水造成污染, 导致在该区域离子浓度较大, NO3--N浓度较高.由于径流量不断增大导致下游NO3--N浓度被稀释, 或接受了其他水源的补给使得NO3--N浓度沿河流方向逐渐减小.在4种类型的水样中, 井水中NO3--N浓度最大, 井水又是该地区主要饮用水, 对当地居民的生活有一定影响, 但整体来说未超过国家规定的饮用水标准, 暂不会造成威胁, 但仍需引起高度关注.通过对黑河流域地表水及地下水中氮素流失情况进行分析, 本研究初步探讨了氮素流失的时空特征, 为流域水资源管理提供了重要的决策依据.但氮素流失和水污染状况受到多种因素影响, 各离子的相关性研究结果表明, 汛期NO3--N与Cl-显示出显著性相关性, 表明农业化肥和牲畜粪便等污染成为该流域水体氮素的主要来源.NO3--N与HCO3-同样具有显著相关性, 反映二者有共同来源, 表明水体中NO3--N不仅来自于农田氮肥的施入等人类活动, 还可能来自于酸性降雨.

5 结论(Conclusions)

1)黑河流域枯水期地表水及地下水的水化学主要为Na+K-Cl-SO4型, 汛期主要为Ca-Mg-HCO3型.汛期离子浓度较枯水期有明显增大, 阳离子主要为Na+, 阴离子在枯水期主要为SO42-, 汛期则转变为HCO3-.

2)流域水体中氮素浓度的时间变化受降水径流关系的影响.2013-2014年枯水期NO3--N浓度普遍小于汛期;同一类型水样在年际变化上明显.同一时期不同类型水样之间NO3--N浓度存在较大差异, 井水浓度最大而径流浓度最小, 但NH4+-N浓度差异不明显.

3)水体中NO3--N的空间分布不均衡, 地表水表现为:上游>中游>下游, 而地下水则为:下游>中游>上游.地表水的污染程度在不同河段有差异, 上游污染严重, 中、下游普遍达Ⅳ类水质标准.地下水NO3--N浓度在不同河段均未超过国家标准, 但地下水的EC值普遍高于地表水且大部分在300~1000 μS·cm-1之间, 说明水质状况较差.因此, 流域地表水健康及地下水安全问题应引起重视.

3)氮素主要以NO3--N的形式流失, NH4+-N浓度在整个研究区域普遍较低, 在不同时期均较稳定、波动小, 均值为0.11 mg·L-1, 大部分达Ⅲ类水质标准.

参考文献
[${referVo.labelOrder}] Arauzo M, Martinez-Bastida J J. 2015. Environmental factors affecting diffuse nitrate pollution in the major aquifers of central Spain:groundwater vulnerability vs[J]. groundwater pollution[J]. Environmental Earth Sciences, 73(12) : 8271–8286. DOI:10.1007/s12665-014-3989-8
[${referVo.labelOrder}] Andersen H E, Kronvang B, Larsen S E. 1999. Agricultural practices and diffuse nitrogen pollution in Denmark:Empirical leaching and catchment models[J]. Water Science and Technology, 39(12) : 257–264. DOI:10.1016/S0273-1223(99)00342-X
[${referVo.labelOrder}] Alberts E E, Schuman G E, Burwell R E. 1978. Seasonal runoff losses of nitrogen and phosphorus from Missouri Valley loess watersheds[J]. Journal of Environmental Quality, 7(2) : 203–208.
[${referVo.labelOrder}] 曹群虎, 段长林, 君广斌, 等. 2009. 长武县富士果园施肥现状调查及建议[J]. 农家之友, 2009, 5(16) : 34–35.
[${referVo.labelOrder}] 蔡崇法, 丁树文, 史志华, 等. 2013. GIS支持下三峡库区典型小流域土壤养分流失量预测[J]. 水土保持学报, 2013, 15(1) : 9–12.
[${referVo.labelOrder}] Cao W Z, Hong H S, Zhang Y Z, et al. 2006. Anthropogenic nitrogen sources and exports in a village-scale catchment in Southeast China[J]. Environmental Geochemistry and Health, 28(1/2) : 45–51.
[${referVo.labelOrder}] Chen J Y, Tang C Y, Sakura Y, et al. 2005. Nitrate pollution from agriculture in different hydrogeological zones of the regional groundwater flow system in the North China Plain[J]. Hydrogeology Journal, 13(3) : 481–492. DOI:10.1007/s10040-004-0321-9
[${referVo.labelOrder}] 程立平, 刘文兆. 2012. 黄土塬区几种典型土地利用类型的土壤水稳定同位素特征[J]. 应用生态学报, 2012, 23(3) : 651–658.
[${referVo.labelOrder}] Di H J, Cameron K C. 2002. Nitrate leaching in temperate agroecosystems:sources, factors and mitigating strategies[J]. Nutrient Cycling in Agroecosystems, 64(3) : 237–256. DOI:10.1023/A:1021471531188
[${referVo.labelOrder}] 郭开秀, 姚春霞, 陈亦, 等. 2011. 上海市秋季蔬菜硝酸盐浓度及风险摄入评估[J]. 环境科学, 2011, 32(4) : 1177–1181.
[${referVo.labelOrder}] Guo H Y, Zhu J G, Wang X R, et al. 2004. Case study on nitrogen and phosphorus emissions from paddy field in Taihu Region[J]. Environmental Geochemistry and Health, 26(2) : 209–219. DOI:10.1023/B:EGAH.0000039583.71423.b4
[${referVo.labelOrder}] 耿金, 陈建生, 张时音. 2013. 赤水河上游流域水化学变化与离子成因分析[J]. 水文, 2013, 34(1) : 44–50.
[${referVo.labelOrder}] 胡宏祥, 洪天求, 刘路. 2009. 水土流失量和养分流失量的预测[J]. 环境科学研究, 2009, 9(3) : 356–361.
[${referVo.labelOrder}] 韩贵琳, 刘丛强. 2000. 贵州乌江水系的水文地球化学研究[J]. 中国岩溶, 2000, 19(1) : 35–43.
[${referVo.labelOrder}] 侯昭华, 徐海, 安芷生. 2009. 青海湖流域水化学主离子特征及控制因素初探[J]. 地球与环境, 2009, 37(1) : 11–19.
[${referVo.labelOrder}] 韩梅, 周小平, 程媛华, 等. 2014. 环太湖主要河流氮素组成特征及来源[J]. 环境科学研究, 2014, 27(12) : 1450–1457.
[${referVo.labelOrder}] Haro P S, Albert L C, Velazquez P D. 2010. Fertilizer standards for controlling groundwater nitrate pollution from agriculture:El Salobral-Los Llanos case study, Spain[J]. Journal of Hydrology, 392(3) : 174–187.
[${referVo.labelOrder}] Hu Q, Paul W, Wim V. 2000. Removal of nitrate from groundwater by cyanobacteria:quantitative assessment of factors influencing nitrate uptake[J]. Applied and Environmental Microbiology, 66(1) : 133–139. DOI:10.1128/AEM.66.1.133-139.2000
[${referVo.labelOrder}] 金洁, 杨京平. 2005. 从水环境角度探析农田氮素流失及控制对策[J]. 应用生态学报, 2005, 16(3) : 579–582.
[${referVo.labelOrder}] 金赞芳, 王飞儿, 陈英旭. 2004. 城市地下水硝酸盐污染及其成因分析[J]. 土壤学报, 2004, 41(2) : 252–258.
[${referVo.labelOrder}] Kumazawa K. 2002. Nitrogen fertilization and nitrate pollution in groundwater in Japan:Present status and measures for sustainable agriculture[J]. Nutrient Cycling in Agroecosystems, 63(2/3) : 129–137. DOI:10.1023/A:1021198721003
[${referVo.labelOrder}] Laftouhi N E, Marnik V, Mohammed J, et al. 2003. Groundwater nitrate pollution in the Essaouira Basin (Morocco)[J]. Comptes Rendus Geoscience, 335(3) : 307–317. DOI:10.1016/S1631-0713(03)00025-7
[${referVo.labelOrder}] 陆海明, 尹澄清, 王夏晖, 等. 2008. 华北半干旱区2个农业流域地表氮素流失特征的对比研究[J]. 环境科学, 2008, 29(10) : 2689–2695.
[${referVo.labelOrder}] 李志, 刘文兆, 王秋贤. 2008. 黄土塬区不同地形部位和土地利用方式对土壤物理性质的影响[J]. 应用生态学报, 2008, 19(6) : 1303–1308.
[${referVo.labelOrder}] 李志, 刘文兆, 郑粉莉. 2010. 基于CA-Markov模型的黄土塬区黑河流域土地利用变化[J]. 农业工程学报, 2010, 26(1) : 346–352.
[${referVo.labelOrder}] 刘相超, 宋献方, 夏军, 等. 2008. 怀沙河流域地表地下水化学特征与地下水资源保护[J]. 重庆交通大学学报:自然科学版, 2008, 27(2) : 305–309.
[${referVo.labelOrder}] 陆海明, 尹澄清, 王夏晖, 等. 2008. 于桥水库周边农业小流域氮素流失浓度特征[J]. 环境科学学报, 2008, 28(2) : 349–355.
[${referVo.labelOrder}] Lu H M, Yin C Q, Wang W D, et al. 2007. comparative study of nutrient transfer via surface runoff from two small agricultural catchments in north China[J]. Environmental Geology, 52(8) : 1549–1558. DOI:10.1007/s00254-006-0599-0
[${referVo.labelOrder}] 马俊杰, 牛建军. 2003. 黑河水库汇流区土壤侵蚀评判方法研究[J]. 水土保持通报, 2003, 23(4) : 24–27.
[${referVo.labelOrder}] Maila Y A, El-nahal I, Aimer M R. 2004. Seasonal variations and mechanisms of groundwater nitrate pollution in the Gaza Strip[J]. Environmental Geology, 47(1) : 84–90. DOI:10.1007/s00254-004-1136-7
[${referVo.labelOrder}] Ming H H, Stallard R F, Edmond J M. 1982. Major ion chemistry of some large Chinese rivers[J]. Natural, 298(5) : 550–553.
[${referVo.labelOrder}] Nolan B T, Stoner J D. 2000. Nutrients in groundwaters of the conterminous United States, 1992-1995[J]. Environmental Science & Technology, 34(7) : 1156–1165.
[${referVo.labelOrder}] Peter J T, Jody S B, Keith L W, et al. 2003. Nitrate in groundwaters of intensive agricultural areas in coastal Northeastern Australia[J]. Agriculture, Ecosystems & Environment, 94(1) : 49–58.
[${referVo.labelOrder}] Ranjith P U, Peter P M, Harold E G, et al. 2006. Nitrogen losses in runoff from three adjacent agricultural watersheds with claypan soils[J]. Agriculture, Ecosystems & Environment, 117(1) : 39–48.
[${referVo.labelOrder}] RAN N S. 2006. Nitrate pollution and its distribution in the groundwater of Srikakulam district, Andhra Pradesh, India[J]. Environmental Geology, 51(4) : 631–645. DOI:10.1007/s00254-006-0358-2
[${referVo.labelOrder}] Stallard R F, Edmond J M. 1981. Geochemistry of the Amazon:1[J]. Precipitation chemistry and the marine contribution to the dissolved load at the time of peak discharge[J].Journal of Geophysical Research:Oceans (1978-2012), 86(10) : 9844–9858.
[${referVo.labelOrder}] Sugimoto A, Naito D, Yanagisawa N, et al. 2003. Characteristics of soil moisture in permafrost observed in East Siberian taiga with stable isotopes of water[J]. Hydrological Processes, 17(6) : 1073–1092. DOI:10.1002/(ISSN)1099-1085
[${referVo.labelOrder}] Surindra S, Preeti B, Sushma S, et al. 2009. Nitrate contamination in groundwater of some rural areas of Rajasthan, India[J]. Journal of Hazardous Materials, 171(1) : 189–199.
[${referVo.labelOrder}] Tian Y H, Yin B, Yang L Z, et al. 2007. Nitrogen runoff and leaching losses during rice-wheat rotations in Taihu Lake region, China[J]. Pedosphere, 17(4) : 445–456. DOI:10.1016/S1002-0160(07)60054-X
[${referVo.labelOrder}] 王晓燕, 王晓峰, 汪清平, 等. 2004. 北京密云水库小流域非点源污染负荷估算[J]. 地理科学, 2004, 24(2) : 227–231.
[${referVo.labelOrder}] 王庆锁, 孙东宝, 郝卫平, 等. 2011. 密云水库流域地下水硝态氮的分布及其影响因素[J]. 土壤学报, 2011, 48(1) : 141–150.
[${referVo.labelOrder}] 叶宏萌, 袁旭音, 葛敏霞, 等. 2010. 太湖北部流域水化学特征及其控制因素[J]. 生态环境学报, 2010, 19(1) : 23–27.
[${referVo.labelOrder}] Zhao T Q, Xu H S, He Y X, et al. 2009. Agricultural non-point nitrogen pollution control function of different vegetation types in riparian wetlands:A case study in the Yellow River wetland in China[J]. Journal of Environmental Sciences, 21(7) : 933–939. DOI:10.1016/S1001-0742(08)62364-5
[${referVo.labelOrder}] 朱亚莉, 慕芳. 2015. 长武县小麦施肥现状及对策[J]. 现代农业科技, 2015, 23(4) : 220–220.
[${referVo.labelOrder}] 张亚丽, 张依章, 张远, 等. 2014. 浑河流域地表水和地下水氮污染特征研究[J]. 中国环境科学, 2014, 34(1) : 170–177.
[${referVo.labelOrder}] 张芳, 李发东, 李静, 等. 2013. 德州引黄灌区主要河系水化学空间特征分析[J]. 中国生态农业学报, 2013, 21(4) : 487–493.
[${referVo.labelOrder}] 张继宗, 张维理, 雷秋良, 等. 2009. 太湖平原农田区域地表水特征及对氮磷流失的影响[J]. 生态环境学报, 2009, 18(4) : 1497–1503.