环境科学  2021, Vol. 42 Issue (7): 3232-3241   PDF    
青岛市农区地下水硝态氮污染来源解析
寇馨月1, 丁军军1, 李玉中1,2, 毛丽丽1, 李巧珍1, 徐春英1, 郑欠1, 庄姗1     
1. 中国农业科学院农业环境与可持续发展研究所, 农业部旱作节水重点实验室, 北京 100081;
2. 中国农业科学院环境稳定同位素实验室, 北京 100081
摘要: 为了提高作物产量,肥料大量投入在农业种植区日益普遍,导致了农区地下水硝态氮(NO3--N)污染.农业面源污染是地下水硝态氮污染的主要原因.为了保障饮用水安全,明确农区硝态氮污染的来源是十分必要的.本研究分别于2009年和2019年在青岛农区随机选取35个采样点,借助反距离加权法(IDW)对硝态氮含量进行空间分布分析,通过测定氮、氧同位素进行溯源,运用SIAR模型量化污染源的贡献率.结果表明,青岛市地下水硝态氮含量(平均值)由2009年的38.49 mg·L-1降低为2019年的22.37 mg·L-1,但仍高于世界卫生组织(WHO)规定的饮用水中硝态氮的最大允许含量.2009年和2019年硝态氮含量都呈现由南向北逐渐增加的趋势,南部污染轻,北部污染重.δ15N-NO3-δ18O-NO3-的交叉图显示青岛市地下水硝态氮主要来源是化肥、土壤氮、粪肥和污水.水同位素表明降水是青岛市地下水的主要来源.贝叶斯混合模型(SIAR模型)表明污染源贡献率为:粪肥和污水(47.42%)>土壤氮(27.80%)>化肥(14.35%)>大气氮沉降(10.43%).从2009~2019年青岛市地下水质量得到了改善,但硝态氮污染状况仍不容忽视,应根据硝态氮污染来源,有针对性地防治以确保农区饮用水安全和农业的可持续发展.
关键词: 氮同位素      氧同位素      硝态氮来源      SIAR模型      青岛农区     
Identifying the Sources of Groudwater NO3--N in Agricultural Region of Qingdao
KOU Xin-yue1 , DING Jun-jun1 , LI Yu-zhong1,2 , MAO Li-li1 , LI Qiao-zhen1 , XU Chun-ying1 , ZHENG Qian1 , ZHUANG Shan1     
1. Key Laboratory of Dryland Agriculture, Ministry of Agriculture, Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081, China;
2. Environmental Stable Isotope Laboratory, Chinese Academy of Agricultural Sciences, Beijing 100081, China
Abstract: To increase crops yields, applying large amounts of fertilizers has become increasingly common in agricultural regions, resulting in NO3--N groundwater pollution. Agricultural non-point pollution is the main source of groundwater NO3--N pollution. To ensure drinking water safety and quality, it is crucial to clarify the sources of NO3--N pollution in agricultural regions. In this study, 35 sampling sites were randomly selected in the Qingdao agricultural area in 2009 and 2019. The spatial distribution of NO3--N concentration was analyzed by the inverse distance weighting method (IDW). The nitrogen and oxygen isotopes were used as a tool to trace sources of NO3--N and the SIAR model was used to quantify contribution proportion of pollution sources. The results showed that the concentration of NO3--N (average) in groundwater in Qingdao has been reduced from 38.49 mg·L-1 in 2009 to 22.37 mg·L-1 in 2019, but it is still higher than the maximum allowable concentration of NO3--N in drinking water set by the World Health Organization (WHO). The NO3--N concentration gradually increased from south to north both in 2009 and 2019. The cross diagram of δ15N-NO3- and δ18O-NO3- show that the main sources of NO3--N in groundwater in Qingdao are chemical fertilizers, soil nitrogen, and manure and sewage. Water isotopes indicate that precipitation was the main source of groundwater in Qingdao. The SIAR model results indicated that the contribution of each source ranked as follows: manure and sewage (47.42%) > soil nitrogen (27.80%) > chemical fertilizer (14.32%) > atmospheric nitrogen depositions (10.43%). From 2009 to 2019, the quality of groundwater in Qingdao has been improved, but NO3--N pollution still cannot be ignored. According to the results, prevention and control should be made to ensure the safety of drinking water and the sustainable development of agriculture.
Key words: δ15N      δ18O      NO3--N source      SIAR model      agricultural area in Qingdao     

随着农业的发展以及“乡村振兴”政策的实行, 农业生产快速发展带来了一系列环境问题[1~3].为了追求高产过量施用化学肥料和有机肥料, 造成了农业面源污染[4], 土壤中大量硝态氮残留物渗入地下水体, 导致地下水硝态氮污染.人体饮用NO3--N含量超标的地下水, 会患血红蛋白血症[5, 6]、肝病和癌症等疾病[7]; 而在水体中的NO3--N会导致许多环境和生态问题, 如水体的富营养化[8]以及河流生态系统中的物种灭绝[9].农区地下水中的硝态氮污染已成为迫切关注的问题[10].因此, 明确地下水硝态氮污染状况, 追溯和量化硝态氮污染源对于防治农区硝态氮、保障饮用水安全和促进农业可持续发展具有重要意义.

传统的溯源方法是通过调查当地的土地利用类型、地面耕作方式并结合地下水的硝态氮的含量特征确定其污染源.得到的结果多为间接和不确定的, 较粗糙, 从而导致溯源的结果不太准确[3, 11, 12].由于不同来源的硝态氮具有不同的N、O同位素组成[4, 13], 目前对硝态氮污染进行溯源是利用NO3-中的δ15 N和δ18 O[3, 14~20], 由于每种作用引起的分馏程度不同[21], 因此应考虑同位素分馏的因素.双同位素与传统方法[22]相结合可以得到相对准确的溯源结果.为了估算硝态氮源的比例贡献, SIAR模型已广泛用于量化硝态氮来源的贡献率[23~25].

本文以山东省的经济中心, “一带一路”新亚欧大陆桥经济走廊主要节点城市“青岛”为研究对象, 旨在探究青岛市农区地下水硝态氮污染的总体状况, 明确近十年前后青岛市地下水硝态氮污染的对比变化及时空分布特征, 依据水化学特征和氮氧同位素对青岛农业种植区进行了硝态氮溯源, 并利用贝叶斯混合模型量化污染来源的贡献率.本研究实现了由定性到定量, 以期为青岛市防治地下水硝态氮污染提供了坚实的理论依据、对于保障该地区居民的饮水安全具有现实意义和对促进农业的可持续发展具有长远意义.

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

青岛市(119°30′~121°00′E, 35°35′~37°09′N), 地处山东半岛东南, 东、南濒临黄海, 东北与烟台市毗邻, 西与潍坊市相连, 西南与日照市接壤; 总面积为11 282 km2.青岛为海滨丘陵城市, 地势东高西低, 南北两侧隆起, 中间低洼.温带季风性气候, 由于临海带有海洋性特征.青岛共有大小河流224条, 均为季风区雨源型, 多为独立入海的山溪性小河.农作物耕作面积68.37万hm2, 其中粮食作物耕作面积47.50万hm2, 蔬菜耕作面积10.75万hm2.青岛市蔬菜产量大约是粮食产量的2倍[26].

1.2 采样及样品测定

青岛市采样点的设置根据均匀布点、局部加密的原则, 于2009年11月和2019年4月在青岛各地水井采集水样35个, 覆盖青岛所有市(区、县), 每个采样点均用GPS定位, 记录采样井位置、深度和农田利用信息等, 采样点示意见图 1.用0.45 μm滤膜过滤后, 将滤液储存在聚乙烯瓶中, 并在-20℃冰箱中冷冻保存.用流动注射分析仪测量水样的NO3--N含量.使用多参数仪测量pH和电导率(EC).用离子色谱仪(940 professional IC Vario)测量Cl-的含量.δ15 N-NO3-δ18 O-NO3-的测定采用反硝化细菌法[27, 28], 于中国农业科学院农业环境与可持续发展研究所环境稳定同位素实验室完成.使用中国农业科学院环境稳定同位素实验室的超高精度液态水和水蒸气同位素分析仪(PicarroL1115-I, 美国)测定δD-H2 O和δ18 O-H2 O. δD-H2 O和δ18 O-H2 O的分析精度分别为0.50‰和0.10‰.

数据来源于国家科技基础条件平台:国家地球系统科学数据中心(http://www.geodata.cn) 图 1 研究区土地利用及地下水采样点分布示意 Fig. 1 Distribution of land use and groundwater sampling sites in the study area

1.3 δ值和SIAR模型

稳定的同位素比以δ表示如下[29, 30]

(1)

式中, R样品R标准分别是样本和标准品的15 N/14 N或18 O/16 O比, 标准分别是大气氮(N2)和维也纳标准平均海水(VSMOW).

目前, 贝叶斯混合同位素模型(SIAR)已广泛用于量化地下水硝态氮来源的贡献率[1, 31~33], 且能够计算多于3个污染来源的贡献率, 同时也能考虑到同位素分馏作用产生的影响, 许多学者使用SIAR成功估算了多种硝态氮源的贡献率[5, 34, 35], 该模型表示为[36]

(2)
(3)
(4)
(5)

式中,Xij表示混合物i同位素jδ值;pk表示来源k的比例; qjk是同位素j在第k种源中的浓度; Sjk表示第k个来源的第j种同位素的δ值, 服从均值为μ方差为ω的正态分布; Cjk表示第k个来源的j同位素的分馏系数, 服从均值为λ方差为τ的正态分布; εij是残余误差, 表示其他各个混合物间无法量化的方差, 其均值和标准差σ在通常情况下均为0.

2 结果与分析 2.1 水化学

对青岛市35个地下水样品进行了测量, 其理化结果见表 1.地下水的pH值为碱性(7. 67~10.17), 平均值为8.20. EC值在224.40~4 648.10 μS ·cm-1之间变化, 平均值为1 101.70 μS ·cm-1, EC值高, 表明受到了硝态氮污染[25]; Cl-值在18.85~1 111.24 mg ·L-1之间变化, 平均值为159.97 mg ·L-1.

表 1 青岛农区地下水样品的水质特征1) Table 1 Water quality characteristics of groundwater samples in Qingdao agricultural area

2.2 青岛市农区地下水硝态氮时间变化

硝态氮含量是评估水质的重要指标.世界卫生组织(WHO)规定饮用水中硝态氮的最大允许含量为11.30 mg ·L-1[37].《饮用水卫生标准》(GB 5749-2006)规定, 饮用水硝态氮的限值为10 mg ·L-1, 地下水源硝态氮的限值为20 mg ·L-1, 这也是中国饮用水(GB/T 14848-2017)的地下水水质标准.中国饮用水标准(GB/T 14848-2017)指出, 当硝态氮的含量大于3 mg ·L-1时, 表明它受到人为因素的影响[38].水中硝态氮的含量大小与水体的质量密切相关, 可以根据水体中的硝态氮含量来评估地下水的质量[39, 40].根据以上标准, 为了更好地明确水质分类, 本研究将地下水水质分为以下5个等级(表 2).

表 2 青岛市农区地下水质量分类1) Table 2 Groundwater quality classification in Qingdao agricultural area

对2009年和2019年收集的35个水样进行的综合评估表明(表 2), 2009年Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ类水的比例分别为0%、8.57%、14.29%、20 %和57.14%.饮用水占22.86%, 污染水占77.14%.57.14%的地下水样品严重超标, 91.43%超过了WHO规定的标准.2019年Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ类水的比例分别为5.71%、5.71%、34.29%、25.71%和28.57%.饮用水占45.71%, 污染水占54.28%.大约28.57%的地下水样品严重超标.85.71%超过了WHO规定的标准.

与2009年相比(表 2), Ⅰ类水增加了, Ⅱ类水减少了, Ⅲ类水增加了, Ⅳ类增加了, Ⅴ类水减少了.饮用水的比例提高, 但是水质不高, 大多数位于10~20 mg ·L-1之间.总体而言, 2019年青岛市硝态氮污染水平低于2009年, 表明地下水水质有所改善(图 2).但是, 污染水的比例仍然很高, 地下水的状况仍不容忽视.

图 2 2009年和2019年青岛市农区地下水硝态氮空间分布 Fig. 2 Spatial distribution of groundwater NO3--N in Qingdao agricultural area in 2009 and 2019

2.3 青岛市农区地下水硝态氮空间变化

2009年和2019年青岛市农区都受到了硝态氮的污染, 2009年污染更为严重.2009年、2019年青岛市农区地下水中的NO3--N含量范围分别为6.85~161 mg ·L-1和2.62~46.20 mg ·L-1, 平均值分别为38.49 mg ·L-1和22.37 mg ·L-1.

在2009年和2019年(图 2), 青岛市硝态氮含量都呈现由南向北逐渐增加的趋势, 南部污染轻, 北部污染重.2009年, 平度、莱西、即墨、城阳、胶州北部和胶南东南部属于严重超标, 污染严重, 大部分区域硝态氮含量超过了30 mg ·L-1.2019年青岛市大部分区域硝态氮含量低于30 mg ·L-1符合Ⅳ类水的标准, 适用于农业和工业用水, 胶南、胶州南部和即墨东部达到了最低饮用水标准.

2009~2019年硝态氮含量存在明显的空间变化, 青岛市硝态氮污染总体降低, 局部仍保持不变或增加.但是, 2019年平均值为22.37 mg ·L-1, 仍然超过了最低饮用水20 mg ·L-1的标准(GB/T 14848-2017).总体上, 青岛的地下水硝态氮由38.49 mg ·L-1降至22.37 mg ·L-1硝态氮污染减轻, 水质得到了改善.

2.4 青岛市农区地下水同位素特征 2.4.1 水体氮氧同位素

水中硝态氮的来源通常有降水、化肥、土壤氮、污水和粪肥[41~43].由于粪肥和污水的氮、氧同位素的范围重叠且同位素组成相似, 因此将粪肥和污水视为同一污染源[44~46].本研究的δ15 N-NO3-值范围为0.47‰ ~29.03‰, 平均值为9.38‰. δ18 O-NO3-值范围为-2.54‰ ~37.02‰, 平均值为10.51‰.图 3表明δ15 N-NO3-分布呈现北部高、南部低的特点, 大部分地区低于15‰, 崂山区东部最高, 胶南市南部、胶州市北部、平度市和莱西区的δ15 N-NO3-值较低.δ18 O-NO3-值呈现北部和南部低、中部高的特点, 平度市西南部、胶州西南部和胶南西北部出现部分高值, 莱西区、崂山区大部和胶南市南部的δ18 O-NO3-值较低.

图 3 青岛市农区地下水δ15 N-NO3-δ18 O-NO3-的空间分布 Fig. 3 Spatial distribution of groundwater δ15 N-NO3- and δ18 O-NO3- in Qingdao agricultural area

2.4.2 水体氢氧同位素

表 1所示, 地下水样品的δD-H2 O值范围为-66.55‰ ~-42.24‰, 平均值为-57.57‰.地下水样品的δ18 O-H2 O值在-9.94‰ ~-6.62‰之间, 平均值为-9.06‰.δD-H2 O和δ18 O-H2 O关系(图 4)可以追踪水体的补给[47].根据全球降水同位素网络(GNIP, IAEA)提供的青岛附近站点的数据计算得出当地气象水位线(LMWL)δD=6.29δ18 O-H2 O-3.63.如图 4所示, 研究区域的地下水样品的δD-H2 O和δ18 O-H2 O值位于全球气象水位线(GMWL)δD=8.1δ18 O-H2 O+10.8[48]和当地气象水位线(LMWL)附近, 表明降水是研究区域地下水的主要来源.地下水的δD-H2 O和δ18 O-H2 O值略高于GMWL和LMWL, 表明补给过程中地下水经历了蒸发[5, 32].δD-H2 O和δ18 O-H2 O的拟合回归方程为y=7.16x+7.34(R2=0.84), 斜率和截距都小于GMWL, 再次证明了蒸发的影响[4].

图 4 δD-H2 O和δ18 O-H2 O的关系 Fig. 4 Relationship between δD-H2 O and δ18 O-H2 O

3 讨论 3.1 青岛市农区地下水硝态氮硝化与反硝化过程的判别

δ15 N-NO3-δ18 O-NO3-值主要受到硝态氮来源的影响, 但生物介导的过程(矿化、固氮、硝化和反硝化)是导致氮和氧同位素分馏的重要因素[20].反硝化过程会影响δ15 N和δ18 O的值, 从而改变溯源结果.因此, 判断反硝化是否发生是确定硝态氮污染源的重要前提[49].先前的研究表明, 微生物反硝化将增加水样品中的δ15 N-NO3-δ18 O-NO3-值, 降低硝态氮含量, 并产生1.3 ∶1~2.1 ∶1的δ15 N-NO3-δ18 O-NO3-比值[50, 51].如图 7所示, 1.3 ∶1~2.1 ∶1之间只有5个采样点, 所以青岛的地下水反硝化不明显, 在研究过程中反硝化作用可以忽略.

图 5 δ18 O-H2 O和δ18 O-NO3-的关系 Fig. 5 Relationship between δ18 O-H2 O and δ18 O-NO3-

图 6 青岛市农区氯离子空间分布及Cl-和NO3-/Cl-的关系 Fig. 6 Spatial distribution of chloride ions and the relationship between Cl- and NO3-/Cl- in Qingdao agricultural area

图 7 典型硝态氮来源氮氧同位素范围以及青岛农区δ15 N-NO3-δ18 O-NO3-关系[20] Fig. 7 Range of nitrogen and oxygen isotopes of typical NO3--N sources and the relationship between δ15 N-NO3- and δ18 O-NO3- in Qingdao agricultural region

从理论上讲, 微生物硝化产生的NO3- 1/3来自大气中的O2, 2/3来自H2 O中的O.微生物硝化产生的δ18 O-NO3-值可通过以下公式计算[43, 47, 52, 53]:

(6)

其中, 大气O2中的δ18 O值为23.5‰, H2 O中的O中的δ18 O值为实测值.

硝化过程可以使用δ18 O-NO3-值确定[54].根据硝化作用计算出的δ18 O-NO3-理论值约为2‰[54, 55].如图 5所示, 大约一半的样品接近理论线[29], 且位于-10‰ ~10‰之间, 这表明青岛市地下水中的硝态氮主要来自微生物硝化作用的合成肥料[1].一些样品的δ18 O-NO3-值高于理论值, 这可能是由于δ18 O-NO3-值受到大气来源影响[20, 29], 或者是混合了其他硝态氮来源所致[32].

3.2 青岛市农区地下水NO3--N源解析

氯离子是变化较小保守示踪剂[22], 且氯离子含量和硝态氮含量的空间分布规律类似[图 2(b)图 6(a)], 可用于鉴定硝态氮来源.在两种土地利用类型中, 灌溉水中氯离子含量高(216.89 mg ·L-1), 饮用水中氯离子含量低(92.38 mg ·L-1), 这可能是由于在农田种植施肥引起的.Cl-和NO3-/Cl-的关系图也可以用来鉴别污染来源:高NO3-/Cl-比值低Cl-含量表示水样中的硝态氮主要来源于化肥, 低NO3-/Cl-比值高Cl-含量表明水样中的硝态氮主要来源于生活污水和粪肥[56, 57].由图 6可知, 青岛农区地下水中的硝态氮更多是低NO3-/Cl-比值高Cl-含量, 表明其主要来自于生活污水和粪肥.部分点具有中等的NO3-/Cl-和Cl-含量, 表明其来自于混合污染[25].

表 3总结了前人对δ15 N-NO3-δ18 O-NO3-的研究范围.根据表 3, 本研究中的δ15 N-NO3-δ18 O-NO3-的范围为:大气氮沉降(δ15 N: -15‰ ~10‰, δ18 O: 58‰ ~98‰), 铵肥(δ15 N: -5.90‰ ~6.60‰, δ18 O: -15‰ ~18‰), 化肥(δ15 N: -7.40‰ ~ 6.80‰, δ18 O: 18‰ ~22‰), 土壤氮(δ15 N-NO3-: 3‰ ~9‰, δ18 O: -15‰ ~18‰), 粪肥和污水(δ15 N: 4‰ ~25‰, δ18 O-NO3-: -15‰ ~18‰).

表 3 硝态氮来源的范围/‰ Table 3 Range of NO3--N sources/‰

青岛农区硝态氮同位素特征值以及典型硝态氮来源氮氧同位素范围如图 7所示.δ15 N-NO3-δ18 O-NO3-的交叉表明(图 7), 该研究区的大部分地下水样品均处于在污染源之内, 主要分布在化肥、土壤氮、粪肥和污水的范围中, 可能是这几种污染源混合作用导致青岛农区的硝态氮污染.从图 7知道, 落在土壤氮范围内的点较多, 说明多年的施肥使土地中积累了硝态氮.只有几个点超出了潜在污染源的范围, 落在化肥和降水之间, 这可能是化肥和降水的混合污染造成的.δ15 N-NO3-δ18 O-NO3-值没有落在降水范围内, 表明降水不是地下水中硝态氮的主要来源.按照地下水的用途划分, 大多灌溉水位于化肥和土壤氮的范围内, 饮用水位于粪肥和污水范围内, 这与灌溉水来源多位于大田, 施用化肥较多有关, 而饮用水来源多取样于农户家中, 受到种植蔬菜、养殖或排放生活污水, 所以以落在粪肥和污水范围内更多.按照纬度划分, 青岛市北部硝态氮污染更为严重, 南部较轻, 北部的点落在化肥和土壤氮的范围更多, 南部的点大多落在粪肥和污水范围内.

3.3 青岛市农区硝态氮来源的量化

采用贝叶斯混合模型估算青岛农区地下水硝态氮源的比例贡献.2个同位素(δ15 N和δ18 O)4个来源(化肥、土壤氮、粪肥及污水和大气氮沉降).此外, 可以忽略不计反硝化作用, 将分馏系数Cjk设置为0.不同硝态氮来源的均值和标准差来自于文献[33, 64, 65].

SIAR模型结果显示(图 8), 粪肥和污水(47.42%)>土壤氮(27.80%)>化肥(14.35%)>大气氮沉降(10.43%), 地下水硝态氮主要来源于粪肥和污水, 生活污水和粪肥因其接近的氮氧同位素值, 目前定量模型难区分开, 所以粪肥和污水的贡献率高.化肥的贡献率比较低, 可能是由于化肥氮素更多地储存在土壤氮中, 这是值得今后研究的问题.在饮用水源中, 粪肥和污水(42.38%)>土壤氮(29.57%)>化肥(16.45%)>大气氮沉降(11.60%).饮用水中不同来源的贡献率和青岛农区总的结果类似.在灌溉水源中, 粪肥和污水(38.58%)>土壤氮(32.79%)>化肥(19.38%)>大气氮沉降(9.25%).灌溉水中粪肥和污水的贡献率有所降低, 土壤氮和化肥的贡献率上升, 这可能是多年的农业种植导致氮素的持续污染所致.根据SIAR模型的结果, 青岛市应减少粪肥的使用和污水的排放, 以使青岛农区的硝态氮污染水平能够迅速降低.本研究的结果也提示在追求经济发展注重效益的同时也该关注环境问题, 这样才可以使农业可持续发展不再成为奢望.

图 8 不同污染源的贡献率 Fig. 8 Contribution proportions of various sources

4 结论

(1) 2009年达到可饮用水(Ⅲ类)标准的占22.86%, 57.14%的地下水样品严重超标(Ⅴ类), 91.43%超过了WHO规定的标准. 2019年达到可饮用水(Ⅲ类)标准的占45.71%, 28.57%的地下水样品严重超标(Ⅴ类), 85.71%超过了WHO规定的标准.总体而言, 2019年青岛市农区硝态氮含量污染水平低于2009年, 地下水水质有所改善.

(2) 青岛农区硝态氮含量呈现由南向北逐渐增加的趋势, 南部污染轻, 北部污染重.

(3) 青岛农区硝态氮主要来自于化肥、土壤氮、粪肥和污水, 可能是这几种污染源混合污染.研究区域主要发生了硝化作用, 反硝化过程不是主要因素.

(4) SIAR模型表明, 硝态氮来源的贡献率为:粪肥和污水(47.42%)>土壤氮(27.80%)>化肥(14.35%)>大气氮沉降(10.43%).应在实际农业种植过程中减少粪肥和生活污水的使用, 注意硝态氮污染的防治.氮循环伴随着复杂的分馏, 这影响了可追溯性的准确性.未来的研究应仔细分析分馏以准确得出结果.

参考文献
[1] 明晓星, 杨平恒, 谢世友, 等. 金佛山世界遗产地岩溶地下河系统硝酸盐来源与转化[J]. 湖泊科学, 2019, 31(5): 1299-1309.
Ming X X, Yang P H, Xie S Y, et al. Sources and transformations of nitrate of the subterranean river system in Jinfoshan Karst World Heritage[J]. Journal of Lake Sciences, 2019, 31(5): 1299-1309.
[2] Ding J T, Xi B D, Gao R T, et al. Identifying diffused nitrate sources in a stream in an agricultural field using a dual isotopic approach[J]. Science of the Total Environment, 2014, 484: 10-18. DOI:10.1016/j.scitotenv.2014.03.018
[3] 赵然, 韩志伟, 申春华, 等. 典型岩溶地下河流域水体中硝酸盐源解析[J]. 环境科学, 2020, 41(6): 2664-2670.
Zhao R, Han Z W, Shen C H, et al. Identifying nitrate sources in a typical karst underground river basin[J]. Environmental Science, 2020, 41(6): 2664-2670. DOI:10.3969/j.issn.1000-6923.2020.06.037
[4] 金赞芳, 胡晶, 吴爱静, 等. 基于多同位素的不同土地利用区域水体硝酸盐源解析[J]. 环境科学, 2021, 42(4): 1696-1705.
Jin Z F, Hu J, Wu A J, et al. Identify the nitrate sources in different land areas based on multipleisotopes[J]. Environmental Science, 2021, 42(4): 1696-1705.
[5] Zhang Q Q, Wang H W. Assessment of sources and transformation of nitrate in the alluvial-pluvial fan region of north China using a multi-isotope approach[J]. Journal of Environmental Sciences, 2020, 89: 9-22. DOI:10.1016/j.jes.2019.09.021
[6] Pasten-Zapata E, Ledesma-Ruiz R, Harter T, et al. Assessment of sources and fate of nitrate in shallow groundwater of an agricultural area by using a multi-tracer approach[J]. Science of the Total Environment, 2014, 470-471: 855-864. DOI:10.1016/j.scitotenv.2013.10.043
[7] Zhai Y Z, Zhao X B, Teng Y G, et al. Groundwater nitrate pollution and human health risk assessment by using HHRA model in an agricultural area, NE China[J]. Ecotoxicology and Environmental Safety, 2017, 137: 130-142. DOI:10.1016/j.ecoenv.2016.11.010
[8] Wang Y J, Peng J F, Cao X F, et al. Isotopic and chemical evidence for nitrate sources and transformation processes in a plateau lake basin in Southwest China[J]. Science of the Total Environment, 2020, 711. DOI:10.1016/j.scitotenv.2019.134856
[9] Gautam S, Iqbal M Z. Using stable isotopes of nitrogen to study its source and transformation in a heavily farmed watershed[J]. Environmental Earth Sciences, 2010, 60(1): 11-20. DOI:10.1007/s12665-009-0165-7
[10] Jin Z X, Wang J F, Chen J G, et al. Identifying the sources of nitrate in a small watershed using δ15N-δ18O isotopes of nitrate in the Kelan Reservoir, Guangxi, China[J]. Agriculture, Ecosystems & Environment, 2020, 297. DOI:10.1016/j.agee.2020.106936
[11] 李玉中, 贾小妨, 徐春英, 等. 山东省地下水硝酸盐溯源研究[J]. 生态环境学报, 2013, 22(8): 1401-1407.
Li Y Z, Jia X F, Xu C Y, et al. A study on the source tracing of groundwater nitrate in Shandong province[J]. Ecology and Environmental Sciences, 2013, 22(8): 1401-1407. DOI:10.3969/j.issn.1674-5906.2013.08.021
[12] 盛婷, 杨平恒, 谢国文, 等. 基于δ15N和δ18O的农业区地下河硝酸盐污染来源[J]. 环境科学, 2018, 39(10): 4547-4555.
Sheng T, Yang P H, Xie G W, et al. Nitrate-nitrogen pollution sources of an underground river in karst agricultural area using 15N and 18O isotope technique[J]. Environmental Science, 2018, 39(10): 4547-4555.
[13] Kaown D, Koh D C, Mayer B, et al. Identification of nitrate and sulfate sources in groundwater using dual stable isotope approaches for an agricultural area with different land use (Chuncheon, mid-eastern Korea)[J]. Agriculture, Ecosystems & Environment, 2009, 132(3-4): 223-231.
[14] Blarasin M, Cabrera A, Matiatos I, et al. Comparative evaluation of urban versus agricultural nitrate sources and sinks in an unconfined aquifer by isotopic and multivariate analyses[J]. Science of the Total Environment, 2020, 741. DOI:10.1016/j.scitotenv.2020.140374
[15] 徐志伟, 张心昱, 于贵瑞, 等. 中国水体硝酸盐氮氧双稳定同位素溯源研究进展[J]. 环境科学, 2014, 35(8): 3230-3238.
Xu Z W, Zhang X Y, Yu G R, et al. Review of dual stable isotope technique for nitrate source identification in surface-and groundwater in China[J]. Environmental Science, 2014, 35(8): 3230-3238.
[16] Yang P H, Li Y, Groves C, et al. Coupled hydrogeochemical evaluation of a vulnerable karst aquifer impacted by septic effluent in a protected natural area[J]. Science of the Total Environment, 2020, 658: 1475-1484.
[17] Liu T, Wang F, Michalski G, et al. Using 15N, 17O, and 18O to determine nitrate sources in the Yellow River, China[J]. Environmental Science & Technology, 2013, 47(23): 13412-13421.
[18] Yu Y L, Jin Z, Chu G C, et al. Effects of valley reshaping and damming on surface and groundwater nitrate on the Chinese Loess Plateau[J]. Journal of Hydrology, 2020, 584. DOI:10.1016/j.jhydrol.2020.124702
[19] Li C, Li S L, Yue F J, et al. Identification of sources and transformations of nitrate in the Xijiang River using nitrate isotopes and Bayesian model[J]. Science of the Total Environment, 2019, 646: 801-810. DOI:10.1016/j.scitotenv.2018.07.345
[20] Kendall C, Elliott E M, Wankel S D. Tracing anthropogenic inputs of nitrogen to ecosystems[A]. In: Michener R, Lajtha K (Eds. ). Stable Isotopes in Ecology and Environmental Science (2nd ed. )[M]. Malden: Blackwell Publishing Ltd, 2007. 375-449.
[21] Hosono T, Tokunaga T, Kagabu M, et al. The use of δ15N and δ18O tracers with an understanding of groundwater flow dynamics for evaluating the origins and attenuation mechanisms of nitrate pollution[J]. Water Research, 2013, 47(8): 2661-2675. DOI:10.1016/j.watres.2013.02.020
[22] Meghdadi A, Javar N. Quantification of spatial and seasonal variations in the proportional contribution of nitrate sources using a multi-isotope approach and Bayesian isotope mixing model[J]. Environmental Pollution, 2018, 235: 207-222. DOI:10.1016/j.envpol.2017.12.078
[23] Liu J, Shen Z Y, Yan T Z, et al. Source identification and impact of landscape pattern on riverine nitrogen pollution in a typical urbanized watershed, Beijing, China[J]. Science of the Total Environment, 2018, 628-629: 1296-1307. DOI:10.1016/j.scitotenv.2018.02.161
[24] Yue F J, Li S L, Waldron S, et al. Rainfall and conduit drainage combine to accelerate nitrate loss from a karst agroecosystem: Insights from stable isotope tracing and high-frequency nitrate sensing[J]. Water Research, 2020, 186. DOI:10.1016/j.watres.2020.116388
[25] 金赞芳, 张文辽, 郑奇, 等. 氮氧同位素联合稳定同位素模型解析水源地氮源[J]. 环境科学, 2018, 39(5): 2039-2047.
Jin Z F, Zhang W L, Zheng Q, et al. Contribution of nitrogen sources in water sources by combining nitrogen and oxygen isotopes and SIAR[J]. Environmental Science, 2018, 39(5): 2039-2047.
[26] 青岛市统计局, 国家统计局青岛调查队. 青岛统计年鉴[M]. 北京: 中国统计出版社, 2020.
[27] Casciotti K L, Sigman D M, Hastings M G, et al. Measurement of the oxygen isotopic composition of nitrate in seawater and freshwater using the denitrifier method[J]. Analytical Chemistry, 2002, 74(19): 4905-4912. DOI:10.1021/ac020113w
[28] 徐春英, 李玉中, 郝卫平, 等. 反硝化细菌法结合痕量气体分析仪/同位素比质谱仪分析水体硝酸盐氮同位素组成[J]. 分析化学, 2012, 40(9): 1360-1365.
Xu C Y, Li Y Z, Hao W P, et al. Analysis of nitrogen isotopic composition of nitrate in water by denitrifier method and trace-gas/isotope ratio mass spectrometry[J]. Chinese Journal of Analytical Chemistry, 2012, 40(9): 1360-1365.
[29] Kendall C, McDonnell J J. Isotope tracers in catchment hydrology[M]. Amsterdam: Elsevier, 1998: 519-576.
[30] Kohl D H, Shearer G B, Commoner B. Fertilizer nitrogen: contribution to nitrate in surface water in a corn belt watershed[J]. Science, 1971. DOI:10.1126/science.174.4016.1331
[31] 傅雪梅, 孙源媛, 苏婧, 等. 基于水化学和氮氧双同位素的地下水硝酸盐源解析[J]. 中国环境科学, 2019, 39(9): 3951-3958.
Fu X M, Sun Y Y, Su J, et al. Source of nitrate in groundwater based on hydrochemical and dual stable isotopes[J]. China Environmental Science, 2019, 39(9): 3951-3958. DOI:10.3969/j.issn.1000-6923.2019.09.042
[32] Fadhullah W, Yaccob N S, Syakir M I, et al. Nitrate sources and processes in the surface water of a tropical reservoir by stable isotopes and mixing model[J]. Science of the Total Environment, 2020, 700. DOI:10.1016/j.scitotenv.2019.134517
[33] 王静, 叶寅, 王允青, 等. 利用氮氧同位素示踪技术解析巢湖支流店埠河硝酸盐污染源[J]. 水利学报, 2017, 48(10): 1195-1205.
Wang J, Ye Y, Wang Y Q, et al. Using δ15N and δ18O values to identify sources of nitrate in the Dianbu River in the Chaohu Lake Basin[J]. Journal of Hydraulic Engineering, 2017, 48(10): 1195-1205.
[34] Xue D M, De Baets B, Van Cleemput O, et al. Classification of nitrate polluting activities through clustering of isotope mixing model outputs[J]. Journal of Environmental Quality, 2013, 42(5): 1486-1497. DOI:10.2134/jeq2012.0456
[35] 张妍, 张秋英, 李发东, 等. 基于稳定同位素和贝叶斯模型的引黄灌区地下水硝酸盐污染源解析[J]. 中国生态农业学报, 2019, 27(3): 484-493.
Zhang Y, Zhang Q Y, Li F D, et al. Source identification of nitrate contamination of groundwater in Yellow River Irrigation Districts using stable isotopes and Bayesian model[J]. Chinese Journal of Eco-Agriculture, 2019, 27(3): 484-493.
[36] Parnell A C, Inger R, Bearhop S, et al. Source partitioning using stable isotopes: coping with too much variation[J]. PLoS One, 2010, 5(3). DOI:10.1371/journal.pone.0009672
[37] World Health Organization (WHO). Guidelines for drinking-water quality (4th ed.)[M]. Geneva: WHO, 2011.
[38] Babiker I S, Mohamed M A A, Terao H, et al. Assessment of groundwater contamination by nitrate leaching from intensive vegetable cultivation using geographical information system[J]. Environment International, 2004, 29(8): 1009-1017. DOI:10.1016/S0160-4120(03)00095-3
[39] 刘宏斌, 雷宝坤, 张云贵, 等. 北京市顺义区地下水硝态氮污染的现状与评价[J]. 植物营养与肥料学报, 2001, 7(4): 385-390.
Liu H B, Lei B K, Zhang Y G, et al. Investigation and evaluation on nitrate pollution in groundwater of Shunyi district[J]. Plant Nutrition and Fertilizer Science, 2001, 7(4): 385-390. DOI:10.3321/j.issn:1008-505X.2001.04.005
[40] 贾小妨. 山东省地下水硝酸盐分布规律及溯源研究[D]. 北京: 中国农业科学院, 2010.
Jia X F. Research on distributed pattern and sources of nitrate in groundwater of Shandong province[D]. Beijing: Chinese Academy of Agricultural Sciences, 2010.
[41] Zhang H, Xu Y, Cheng S Q, et al. Application of the dual-isotope approach and Bayesian isotope mixing model to identify nitrate in groundwater of a multiple land-use area in Chengdu Plain, China[J]. Science of the Total Environment, 2020, 717. DOI:10.1016/j.scitotenv.2020.137134
[42] Xue D M, Botte J, De Baets B, et al. Present limitations and future prospects of stable isotope methods for nitrate source identification in surface-and groundwater[J]. Water Research, 2009, 43(5): 1159-1170. DOI:10.1016/j.watres.2008.12.048
[43] Zhang Q Q, Wang H W, Wang L. Tracing nitrate pollution sources and transformations in the over-exploited groundwater region of north China using stable isotopes[J]. Journal of Contaminant Hydrology, 2018, 218: 1-9. DOI:10.1016/j.jconhyd.2018.06.001
[44] Yu L, Zheng T Y, Zheng X L, et al. Nitrate source apportionment in groundwater using Bayesian isotope mixing model based on nitrogen isotope fractionation[J]. Science of the Total Environment, 2020, 718. DOI:10.1016/j.scitotenv.2020.137242
[45] Biddau R, Cidu R, Da Pelo S, et al. Source and fate of nitrate in contaminated groundwater systems: assessing spatial and temporal variations by hydrogeochemistry and multiple stable isotope tools[J]. Science of the Total Environment, 2019, 647: 1121-1136. DOI:10.1016/j.scitotenv.2018.08.007
[46] 王开然, 郭芳, 姜光辉, 等. 15N和18O在桂林岩溶水氮污染源示踪中的应用[J]. 中国环境科学, 2014, 34(9): 2223-2230.
Wang K R, Guo F, Jiang G H, et al. Application of 15N and 18O to nitrogen pollution source in karst water in Eastern Guilin[J]. China Environmental Science, 2014, 34(9): 2223-2230.
[47] 殷超, 杨海全, 陈敬安, 等. 基于水化学和氮氧同位素的贵州草海丰水期水体硝酸盐来源辨析[J]. 湖泊科学, 2020, 32(4): 989-998.
Yin Y, Yang H Q, Chen J A, et al. Tracing nitrate sources with dual isotopes and hydrochemical characteristics during wet season in Lake Caohai, Guizhou Province[J]. Journal of Lake Sciences, 2020, 32(4): 989-998.
[48] Craig H. Isotopic variations in meteoric waters[J]. Science, 1961, 133(3465): 1702-1703. DOI:10.1126/science.133.3465.1702
[49] Nisi B, Vaselli O, Delgado Huertas A, et al. Dissolved nitrates in the groundwater of the Cecina Plain (Tuscany, Central-Western Italy): clues from the isotopic signature of NO3-[J]. Applied Geochemistry, 2013, 34: 38-52. DOI:10.1016/j.apgeochem.2012.12.004
[50] Panno S V, Kelly W R, Martinsek A T, et al. Estimating background and threshold nitrate concentrations using probability graphs[J]. Ground Water, 2006, 44(5): 697-709. DOI:10.1111/j.1745-6584.2006.00240.x
[51] Nikolenko O, Jurado A, Borges A V, et al. Isotopic composition of nitrogen species in groundwater under agricultural areas: a review[J]. Science of the Total Environment, 2018, 621: 1415-1432. DOI:10.1016/j.scitotenv.2017.10.086
[52] Durka W, Schulze E D, Gebauer G, et al. Effects of forest decline on uptake and leaching of deposited nitrate determined from 15N and 18O measurements[J]. Nature, 1994, 372(6508): 765-767. DOI:10.1038/372765a0
[53] Li S L, Liu C Q, Li J, et al. Assessment of the sources of nitrate in the changjiang river, china using a nitrogen and oxygen isotopic approach[J]. Environmental Science & Technology, 2010, 44(5): 1573-1578.
[54] Ogrinc N, Tamse S, Zavadlav S, et al. Evaluation of geochemical processes and nitrate pollution sources at the Ljubljansko polje aquifer (Slovenia): a stable isotope perspective[J]. Science of the Total Environment, 2019, 646: 1588-1600. DOI:10.1016/j.scitotenv.2018.07.245
[55] 李艳利, 杨梓睿, 尹希杰, 等. 太子河下游河流硝酸盐来源及其迁移转化过程[J]. 环境科学, 2019, 39(3): 1076-1084.
Li Y L, Yang Z R, Yin X J, et al. Identification of nitrate sources and the fate of nitrate in downstream areas: a case study in the Taizi River basin[J]. Environmental Science, 2019, 39(3): 1076-1084.
[56] Liu C Q, Li S L, Lang Y C, et al. Using δ15N-and δ18O-values to identify nitrate sources in karst ground water, Guiyang, Southwest China[J]. Environmental Science & Technology, 2006, 40(22): 6928-6933.
[57] Yue F J, Liu C Q, Li S L, et al. Analysis of δ15N and δ18O to identify nitrate sources and transformations in Songhua River, Northeast China[J]. Journal of Hydrology, 2014, 519: 329-339. DOI:10.1016/j.jhydrol.2014.07.026
[58] Nejatijahromi Z, Nassery H R, Hosono T, et al. Groundwater nitrate contamination in an area using urban wastewaters for agricultural irrigation under arid climate condition, southeast of Tehran, Iran[J]. Agricultural Water Management, 2019, 221: 397-414. DOI:10.1016/j.agwat.2019.04.015
[59] Jin Z F, Qin X, Chen L X, et al. Using dual isotopes to evaluate sources and transformations of nitrate in the West Lake watershed, eastern China[J]. Journal of Contaminant Hydrology, 2015, 177-178: 64-75. DOI:10.1016/j.jconhyd.2015.02.008
[60] Heaton T H E, Stuart M E, Sapiano M, et al. An isotope study of the sources of nitrate in Malta's groundwater[J]. Journal of Hydrology, 2012, 414-415: 244-254. DOI:10.1016/j.jhydrol.2011.10.037
[61] Zhang Y, Shi P, Li F D, et al. Quantification of nitrate sources and fates in rivers in an irrigated agricultural area using environmental isotopes and a Bayesian isotope mixing model[J]. Chemosphere, 2018, 208: 493-501. DOI:10.1016/j.chemosphere.2018.05.164
[62] Oelmann Y, Kreutziger Y, Bol R, et al. Nitrate leaching in soil: tracing the NO3- sources with the help of stable N and O isotopes[J]. Soil Biology and Biochemistry, 2007, 39(12): 3024-3033. DOI:10.1016/j.soilbio.2007.05.036
[63] Xing M, Liu W G. Using dual isotopes to identify sources and transformations of nitrogen in water catchments with different land uses, Loess Plateau of China[J]. Environmental Science and Pollution Research, 2016, 23(1): 388-401. DOI:10.1007/s11356-015-5268-y
[64] 李玉中, 王利民, 徐春英, 等. 地下水硝酸盐污染治理研究[M]. 北京: 中国农业科学技术出版社, 2014.
[65] Zhang Q Q, Wang X K, Sun F X, et al. Assessment of temporal and spatial differences of source apportionment of nitrate in an urban river in China, using δ15N and δ18O values and an isotope mixing model[J]. Environmental Science and Pollution Research, 2015, 22(24): 20226-20233. DOI:10.1007/s11356-015-5674-1