科技论文
摘 要： 稳定同位素技术是研究环境和生态系统中元素循环途径的重要方法. 稳定同位素的丰度变化反映了自然界和生物体内混合、分馏双重作用的结果，因此可作为指标计算混合物的来源贡献，或研究造成分馏的化学反应和生物代谢路径. 从20世纪中期确立稳定同位素的基础原理，70年间该技术在地球化学、环境科学、生态学、微生物学、食品科学等领域获得了大量有价值的成果. 其中，水体作为自然环境的重要组成和人类社会的重要资源，已有诸多研究涉及稳定同位素在水环境中污染物溯源、水生态系统元素迁移转化、水生生物营养来源和营养关系等方面的应用. 通过梳理常见的碳、氮、氧稳定同位素在水环境和水生态领域的研究进展，发现污染物和食物来源分析已不局限于定性识别，基于数学模型的混合物组分定量评估方法正得到越来越多的应用；同时，为开展水体脱氮强度和通量估算、水生生物营养级计算及食物网分析，精确测量18O、15N和13C的富集程度、通过试验和调研获取运算所需的基础参数都是关键步骤. 虽然在实际应用中存在待完善之处，稳定同位素技术的前沿研究仍昭示了其整体化、精细化的发展方向. 未来与计算科学的方法学进步相结合，将为水科学研究提供更有力的技术支撑.
关键词 : 稳定同位素;分馏;多元混合模型;硝酸盐溯源;食性分析;营养级;
Abstract： Stable isotope technique is an important method to study element cycle in the environment and ecosystems. The change in the abundance of stable isotopes reflects the results of mixing and fractionation in nature and organisms, therefore, it can be used as an indicator to calculate the source contribution of a mixture, or study the chemical reactions and biological pathways that cause fractionation. Since the establishment of the basic principle of stable isotopes in the middle of the last century, the technology has obtained many valuable achievements in geochemistry, environmental science, ecology, microbiology, food science and other fields during the past 70 years. Among them, there have been many studies on the application of stable isotopes to evaluate the pollutant sources, element migration and transformation, food sources and nutritional relations in aquatic environment, which is an essential component of nature andasignificant resource of human society. By summarizing the research progress of common carbon, nitrogen and oxygen stable isotopes in the field of water environment and ecology, it is found that the analysis of pollutants and food sources is not limited to qualitative identification. In fact, the quantitative assessment method of mixture components based on mathematical models has been increasingly used. Meanwhile, to evaluate the intensity and flux of denitrification, or conduct the trophic level calculation and food-web analysis, it is important to measure the enrichment degree of 18O,15N and 13C accurately, and get the appropriate basic parameters for model by experiment and survey. Although some aspects still need to be optimized for practical applications, the research front of stable isotope technology shows its development direction of integrity and fineness. Combining with methodological progress of computational science, stable isotopes will provide more powerful technical support for water science research.
Keyword： stable isotope; fractionation; multivariate mixed model; nitrate source identification; diet analysis; trophic level;
稳定同位素指核结构稳定、不发生或不易发生放射性衰变、半衰期极长的一系列核素. 1912年，英国物理学家John Thomson发现天然氖由20和22两种不同质量数的原子组成，首次证明了常态元素稳定同位素的存在. 该技术的发展始于20世纪中期的稳定同位素地球化学领域，自然界常见的氢、碳、氮、氧、硫等元素都有超过一个的稳定同位素，研究者测定各类天然物质中元素的同位素组成，分析变化规律并推测原因所在.
近年来，稳定同位素技术在地球化学、环境科学、生态学、微生物学、食品科学等领域得到了广泛利用[1,2,3]. 应用重点并非探究同位素丰度的分布特征，而是将其作为观测手段研究过程、开展溯源. 物理、化学和生物过程都可造成稳定同位素的丰度差异，反之稳定同位素化学性质相同，其变化不影响过程，因此可作为示踪剂为回溯环境、生物等多种因素对元素的共同作用提供科学有效的依据.
水是生物不可或缺的基础物质，水生态环境不仅是环境科学和生态学的研究对象，其污染状况、物种分布等信息也是水资源管理的重要支撑. 随着色谱质谱联用和傅立叶变换核磁共振技术的诞生，同位素的检测手段愈加丰富，研究成本大幅降低，已有大量研究涉及稳定同位素在水体污染物溯源、水生生物食物网研究等不同方向的应用. 聚焦研究中最常见的碳、氮、氧稳定同位素，介绍了分馏和混合的基础理论及模型，水质分析方面，以氮素和硝酸盐为例探讨了N、O稳定同位素及以IsoSource和SIAR为代表的混合模型在污染物迁移转化和溯源中的应用；水生态方面，基于C、N稳定同位素特征、浓缩系数和基线生物选择综述了该技术在营养级、生态位和食性分析中的应用.
1、稳定同位素分馏和混合的基础理论
1.1、 改变同位素丰度变化的驱动力：分馏作用
1.1.1、分馏作用的定义
同位素中子数的差异使其质量和原子间结合力不同，也导致了地球化学及生物化学过程中扩散和反应速率的差异. 最终表现为同位素分馏(isotopic fractionation)效应，即一种元素的同位素以不同比例在不同物相中分配. 通常而言，轻同位素优先参与化学反应，造成动力学非平衡分馏[4]. 而当反应已达到化学平衡，不同化合物对重同位素的竞争导致了热力学平衡分馏.
分馏效应形成了具有不同同位素丰度的物质库. 2H、13C、15N和18O等重同位素在环境中的绝对丰度低，一般接近或低于1%. 为了放大细微变化，便于比较，一般用样品和国际公认标准物质同位素比值的千分差描述稳定同位素含量，记作δ. δ值越大(越正)表示重同位素越多，越小(越负)表示轻同位素越多[5].
1.1.2 、自然分馏
自然界的物理过程和化学反应均会造成同位素分馏. 以氮素为例，挥发、扩散、水解、固氮、硝化、反硝化等反应体系具有不同的同位素分馏特征. 氨化作用造成的分馏仅为±1‰左右. 但硝化作用会在反应物里富集15N，浓缩系数在-12‰到-29‰之间；氨挥发也会导致NH4+里15N含量上升，氮的轻同位素更易于以NH3的形式挥发[6,7]. 李荣富等[8]的综述总结了15N在氮循环过程各环节中的分馏系数，并指出并不是所有环节的同位素分馏都得到了清晰阐述，如硝酸盐异化还原、厌氧氨氧化的分馏效应就暂未明确.
1.1.3 、营养分馏
生物的新陈代谢会导致同位素分馏. 呼吸和排泄优先排出营养元素的轻同位素，造成食物中的重同位素在捕食者体内富集. 研究[9,10]表明. 13C一般富集0~1‰，而15N的浓缩系数通常可达到3‰~4‰，并随着食物链进一步放大. 因此，与食物组成相近的δ13C常被用于追溯水体中捕食者的食物组成和贡献，富集程度更高的δ15N则用于开展水生态系统的营养水平研究[11,12,13,14].
1.2 、计算源贡献的理论及方法：多源混合模型
分馏和混合是稳定同位素研究中的两个重要概念. 分馏作用导致物质库之间的同位素丰度差异，而多个物质库的相互混合进一步造成环境中稳定同位素的复杂分布. 混合模型用数学方法推断混合物中各物质库的贡献率，被广泛应用于污染物定量溯源和捕食者食性分析.
作为最基础的混合模型，质量平衡模型(mass-balance model)可用于n种同位素和n+1个来源的溯源计算. 最常用的是双同位素和3来源的定量评估，可以避免单一同位素的区分度不够(如NO3-里的15N). 但也有研究[15]指出质量平衡模型未考虑不确定性的影响，在混合物具有多种来源时不能使用. 为解决多来源的问题，2003年，Phillips等[16]通过对质量平衡方程组的反复迭代提出了IsoSource模型. IsoSource适用于n种同位素和>n+1个源，获得的不是点估计，而是分布. 2008年，Parnell等人开发的SIAR(stable isotope analysis in R)模型引入了不确定性和同位素分馏的考虑，其基于狄利克雷分布和贝叶斯算法，结合似然函数和逻辑先验分布计算贡献率. 吴文欢等[17]对以上三者的适用情境和使用技巧进行了详细论述.
除质量平衡模型、IsoSource和SIAR外，IsoError、IsoEconc、IsotopeR、MixSIR等模型也在污染物溯源中有所应用，冯建祥、Hopkins等[18,19]的综述比较了它们的原理、特征和优缺点. 但在实际使用中为避免误用，还需仔细评估背景资料是否充分，科学问题、试验设计和模型选择是否合理[20].
2 、N/O双同位素在水体氮素迁移转化及溯源中的应用
2.1、反硝化脱氮
氮循环是生物地球化学循环中的重要环节，生态系统中氮的输入、转化、利用和消除都是被长期探讨的科学课题. N/O双同位素用于脱氮过程研究，提供了用稳定同位素探讨元素迁移转化的范例.
厌氧条件下反硝化菌有偏好地利用轻同位素转化成N2和N2O，造成15N和18O富集在剩余的NO3-中. 二者的同位素分馏系数比值约为2，即δ15N升高1‰，δ18O相应上升0.5‰[6]. Xue等[7]认为，实际情况下15N和18O以1.3∶1到2.1∶1之间的比例富集都可以推论反硝化的发生. 表1总结了一些分馏系数的比值. 综合脱氮前的N/O同位素比例和分馏系数之比，Panno[21]、Houlton等[22]分别试算了密西西比河和森林中因反硝化损失的氮量，为研究氮素迁移转化提供了新方法.
表1水体脱氮过程中氮氧同位素分馏系数对比
2.2、硝酸盐溯源
我国许多地方的地表水和地下水正在面临硝酸盐污染[29,30]. 六成分图、派珀图等是硝酸盐污染分析的经典方法[31]，但随着同位素方法的建立和推广，N/O双同位素被广泛用于氮源定性识别. 1998年Kendall等[6]根据大量试验结果提出了利用N/O双同位素示踪不同来源硝酸盐的经验方法，利用δ15N的差异区分来自肥料和降水中NH4+硝化、土壤氮素、有机肥与污水的硝酸盐，再通过δ18O的差异区分来自硝酸盐肥料和大气降水的硝酸盐. 此方法得到了后续研究的证实[7,32,33].
N/O同位素溯源研究通常需要结合水化学和其他同位素数据. 例如，Min等[34]用δ15N联合NO3-浓度解释了生活污水和化肥对韩国洛东江流域地下水硝酸盐的贡献；Nyilitya等[35]利用NO3-、Cl-浓度、硼(boron, B)浓度和δ11 B探讨了肯尼亚基苏木城和Kano平原地下水硝酸盐来源及转化；Wang等[36]比较了Cl-与NO3-浓度、NO3-/Cl-比值的相关性，以讨论云南程海氮素的来源和反硝化的发生；以及在赣江流域开展的研究[37]中，分析了δ15N和NO3-水平、NO3-/Cl-比值的关系. 此外，研究地点的背景信息(如土地使用情况)[31,37]、排放负荷和入水负荷[38]等也常被用作综合分析.
在氮源定性识别的基础上，能量化污染源贡献率的定量评估方法也逐渐受到青睐. 研究者将一系列数学模型用于地表水和地下水的硝酸盐来源探讨，得到氮源的相对贡献. 近年来，IsoSource和SIAR在我国各地得到了越来越多的应用(见表2).
表2 稳定同位素在硝酸盐定量溯源中的应用
3、 C/N双同位素在水生生物营养谱构建和食性研究中的应用
3.1 、营养级估算
稳定同位素在营养生态学的应用始于20世纪90年代中期. 以营养分馏为基础，聚焦水生生物体内的15N含量可反映生态系统中不同物种的营养水平. Vender等[57]发现传统食物分析和稳定同位素法得到的342种鱼类的营养级无显着性差异，证明了新方法的可靠性. 稳定同位素计算结果用连续数值标示生物在食物网中所处位置，与传统方法获得的正整数相比更能反映其在生态系统中的真实状况[58].
捕食者和食物之间较确定的δ15N差值是稳定同位素用于营养级关系计算的基础，根据实际需要选用有针对性的浓缩系数和基线生物是开展研究的前提. 通过试验或文献调研，Minagawa等[10]提出15N浓缩系数的平均值(±误差)为3.4±1.1‰，其后Post[59]和Mccutchan[60]的结论分别为3.4±0.98‰和2.2±0.18‰. Mccutchan同时发现针对不同的食物结构和捕食者部位，15N呈现出不同的富集模式. Zanden等[61]按照生物分类(鱼类和无脊椎动物)、栖息地(海水和淡水)、调查方式(实验室和野外)及食物类型(肉食和植食)比较了15N的浓缩系数，发现后两者有显着性差异(P<0.05)，但整体在3‰上下. 国内万祎和蔡德陵等[62,63,64]分别通过渤海湾水生生物网调查和鳀鱼养殖试验得出了3.8‰和2.5‰的15N浓缩系数并得到应用，但很多研究还是采用3.4‰[65,66,67].
通常情况下浮游植物的数量变动剧烈，不适合反映水域时间和空间的平均信息，因此基线生物多选用研究水域中常年存在、食性单一的浮游动物或底栖动物[68]. 如李红燕、魏虎进等[63,69]以水体浮游动物优势种中华哲水蚤、太平洋纺锤水蚤、针刺拟哲水蚤作为基线生物，以其δ15N值为基线值. 也有研究[59]认为浮游动物易受外界干扰、季节波动明显，应选择个体较大、生活周期长的螺类和双壳类，如已有实际应用的福寿螺[66]、铜锈环棱螺[65]、珠蚌[70]、栉孔扇贝[67]、翡翠贻贝[68]等. 在国内，稳定同位素技术曾被用于三峡库区、太湖等淡水水域和象山港、胶州湾、流沙湾等海域的连续营养谱构建，所调查的生物、使用的基线生物和浓缩系数等详见表3.
表3 稳定同位素在国内水域营养谱构建中的应用
3.2、生态位研究
利用δ13 C-δ15 N散点图量化营养结构和生态位的方法为食物网研究提供了新思路. 在用δ13 C-δ15 N二维坐标关联同位素含量和营养功能群的基础上[77]，提出δ15N 差值(δ15N Range, NR)，δ13C 差值(δ13C range, CR)、总面积(Total area, TA)、平均离心距离(Mean distance to centroid, CD)、平均最邻近距离(Mean nearest neighbor distance, NND)和平均最邻近距离标准差 (Standard deviation of nearest neighbor distance, SDNND)六个指标，分别表示营养层次、食源多样性、占据生态位或食物网中营养多样性的总量、营养多样性平均水平、群落的整体密度和营养生态位分布的均匀程度[78]. 此方法量化了食物网的营养结构多样性程度和冗余度，有利于评价生态系统中每个物种的功能及生态位变化.
张文博、谢斌等[79,80]使用δ13 C-δ15 N散点图研究了海洋渔场中小型消费者的C/N稳定同位素比值，分别分析了华南海陵湾、陵水湾两个水域和连云港海州湾不同季节的水生生物NR、CR、TA、CD、NND、SDNND时空变化. 盖珊珊、俞雅文等[81,82]将此方法用于不同生物的生态位研究，探讨了两种鱼类和两种蟹类生态位的宽度和重叠程度. 其中，盖珊珊等除NR、CR和TA外，还计算了标准椭圆面积(Standard Ellipse Area, SEA)和营养生态位重叠面积(Overlap Area, OA). SEA被认为比TA更好，且更不易受样本数量干扰[83]. OA量化了生态位的重叠程度，能更好地评价不同物种利用食物资源的竞争强弱.
3.3、 食物源分析
胃含物分析是开展食性研究的经典方法，通过直接解剖动物胃肠道并分析其中残留食物组成了解其食物来源. 实际应用中发现其具有一系列局限性，包括样本量小时偶然性强、不适于难以解剖的小型动物、只能反映较短时间内的摄食状况、食物残渣偏向难消化的食物类型等[84]. 与之相比，生物体内某些组织的稳定同位素具有较长的周转时间，可用来研究长生命周期内消费者对食物的代谢和吸收[85,86].
与定性判断氮源类似，可通过对比捕食者和不同来源食物的同位素值分析食物组成. 例如，当几种潜在食源的δ13C值差异显着，消费者的δ13C落在颗粒有机物和固着藻类之间，可认为这二者是主要食物[65]. 张波等[87]根据生物的δ13C值由栖息水层加深逐渐增大，提出了崂山湾鳅虎鱼不同生命阶段摄食不同深度水生生物的变化规律. 崔莹等[88]结合中华绒螯蟹各发育阶段特点和C/N同位素比值分析其洄游期食物组成，强调了结合稳定同位素与生活习性开展分析的必要性. 定量分析方面，混合模型同样被用于食物贡献率计算，表4总结了部分应用. 此时要注意通过胃含物、文献查阅等传统方法确定食源范围，避免取回的样品缺失重要的食物来源；另外尤其针对15N等营养分馏较明显的同位素，需要选取合适的分馏系数[89].
表4 稳定同位素在水生生物食性研究中的应用
4、 其他特定化合物稳定同位素技术
利用N/O同位素进行硝酸盐溯源是针对化合物(或离子)中某元素开展的稳定同位素分析(compound-specific isotope analysis, CSIA)应用之一. 在环境科学领域，CSIA也被用于卤代烃、多环芳烃等有机污染物的溯源[96]. 而在微生物系统结构和生物化学过程解析方面，Ohkochi等[3]分离纯化了日本Kaiike湖中的光合色素并测定其C/N稳定同位素组成，解释了与色素相关的自养微生物的生活深度、同化途径和生物功能. Isaji等[97]测定了意大利某盐场中营养盐、叶绿素的δ15N变化，以此探讨底栖微生物中氮素的转化途径及铵态氮在初级生产中的循环利用.
前文所述C/N同位素分析食源和营养级的研究都是将个体或组织(如肌肉)作为对象，检测同位素的整体丰度. 目前已知15N浓缩的主要原因是氨基酸代谢里的脱氨基反应有较明显的同位素分馏，导致排泄出更轻的氮同位素[10]. 因此有研究者分离作为生物标记物的氨基酸并测定其δ15N，从中获得与食物链相关的信息. 其理论基础是：生物体内以谷氨酸为代表的一类氨基酸可显着富集15N(约8‰)，而苯丙氨酸等的15N含量随食物链富集程度较弱(约0.4‰)，二者的δ15N值携带了浓缩系数和基线的双重信息[98]. 另有研究[99]认为，利用谷氨酸、脯氨酸、天冬氨酸、赖氨酸等多种氨基酸的δ15N计算营养级可能比仅采用谷氨酸和苯丙氨酸更准确. Ishikawa等[100]总结了氨基酸氮元素的CSIA方法在食物网生态学中的应用. 于灏等[101]综述了以脂肪酸为目标化合物的CSIA在海洋食物网研究中的进展.
随着进一步研究，还出现了针对化合物特定位置元素开展的稳定同位素分析(position -specific isotope analysis, PSIA)，以获得更详细的分子内同位素分布信息[102]. 例如，一氧化二氮中氮素的PSIA方法可被用于研究硝化和反硝化过程，且在一定程度上推算二者发生的比例[103]. PSIA可以解答CSIA提出的问题，加深我们对元素在环境中归趋与去向的理解.
5、结论与展望
稳定同位素技术为研究自然和生物过程中的元素迁移转化提供了新手段，其在水环境和水生态领域得到快速发展，在模型方法逐渐完善的同时积累了大量有价值的数据. 但在实际应用中，还存在以下需要深入解决的问题：
a)磷是水生态系统中的重要元素，与水域富营养化关系密切. 然而除31P外，P的其他同位素均具有放射性，无法类比C/N/O的方法研究水体中P的来源和转化. 考虑到自然界中P通常以磷酸盐形式出现，目前已出现了用PO43-中O的稳定同位素开展P循环分析与来源解析的方法. 但由于该领域起步较晚，经验不足，目前仍存在诸多问题.
b)自然界的同位素比值极其复杂，其原因部分归结于同位素分馏系数随环境条件发生变化，且往往与混合同时发生. 因此需要认识到二者的不可分割，并在研究之前深入理解该水域中影响同位素分馏的可能因素. 以SIAR为代表、基于贝叶斯算法的模型纳入了分馏效应和不确定性的考量，建议优先选用较新的模型进行计算.
c) 正确的参数传入是模型获得可靠结果的基础. 以营养级研究为例，获取适当的浓缩系数对后续估算十分关键. 然而目前在用的许多浓缩系数或为实验室的分馏试验结果，通常是特定生物饲喂单一来源食物所得；或为野外营养级富集结果，受野外环境复杂食源影响较大；或为二者平均，可能会掩盖多种因素对营养分馏的影响. 因此，模型的基础参数需结合实际情况谨慎选用，必要时开展进一步试验和讨论.
整体同位素—CSIA—PSIA的进展揭示了稳定同位素技术深入化、精细化的发展趋势. 随着计算科学的发展，新的算法让基于大量数据的复杂计算变得方便快捷. 未来的研究可能依托多种计算机模型开展，如更加优化的质量守恒模型、分子/原子间相互作用模型、食物网模型等，水体中污染物的转化归趋、及水生态系统各物种相互作用都是尺度巨大、变化众多的复杂过程，方法学的进步将极大有助于稳定同位素为水科学研究提供更加合理和明确的证据.
参考文献
[1] AMBROSE S H,DENIRO M J.The isotopic ecology of East African mammals[J] Oecologia,1986 ,69(3):395-406.
[2]钟敏用碳氮稳定同位素对大米产地溯源的研究[D].大连:大连海事大学.2013:11-13.
[3] OHKOUCHI N,NAKAJIMA Y,OKADA H,et al Biogeochemical processes in the saline meromictic Lake Kaiike ,Japan.implications from molecular isotopic evidences of photosynthetic pigments[J] Environmental Microbiology,2005,7(7):1009-1016.
[4] MECKENSTOCK R U,MORASCH B, GRIEBLER C.et al. Stable isotope fractionation analysis as a tool to monitor biodegradation in contaminated acquifers[J].Journal of Contaminant Hydrology,2004,75(3):215-255.
[5]徐军应用碳、氮稳定性同位素探讨淡水湖泊的食物网结构和营养级关系[D].武汉:中国科学院研究生院(水生生物研究所),2005:3-5.
[6] KENDALL C,MCDONNELL J J.Isotope tracers in catchment hydrology[M] Amsterdam: Elsevier Science,1998:.519-576.
[7] XUE D,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.
[8]李荣富，罗跃辉，曾洪玉,等稳定同位素技术在环境水体氮的生物地球化学循环研究中的应用[J].南京大学学报(自然科学),2016,52(1):16-26.
[9] DENIRO M J,EPSTEIN S Infuence of diet on the distribution of carbon isotopes in animals[J] Geochimica et Cosmochimica Acta,1978,42(5):495-506.
[10] MINAGAWA M,WADA E Stepwise enrichment of 15N along food chains:Further evidence and the relation between δ15N and animal age[J] Geochimica et Cosmochimica Acta,1984,.48(5)-1135-1140
[11] GUZZO M M.HAFFNE G D,SORGE S,et al Spatial and temporal variabiltis of δ13C and δ15N within lower trophic levels of a large lake:implications for estimating trophic relationships of consumers[J].Hydrobiologia,2011,675:41-53.
[12] PERSIC A.ROCHE H,RAMADE F.Stable carbon and nitrogen isotope quantitative structural assessment of dominant species from the Vaccares Lagoon trophic web (Camargue Biosphere Reserve France)[J] Estuarine, Coastal and Shelf Science, 2004,60(2):261-272.
[13] PETERSON B J,FRY B. Stable isotopes in ecosystem studies([J]. Annual Review of Ecology and Systematics, 1987,18.293-320.
[14] STOWASSER G,ATKINSON A.MCGILL R A R,et al.Food web dynamics in the Scotia Sea in summer:a stable isotope study[J]. Deep Sea Research Part l:Topical Studies in 0ceanography,.2012, 59/60:208-221.
[15] MOORE J W,EMMENS B X Incorporating uncertainty and prior information into stable isotope mixing models[J] Ecology Letters,2008,11(5).470-480.
[16] PHILLIPS D L,GREGG J W. Source partitining using stable isotopes: coping with too many sources[J]. Oecologia,2003, 136(2):261-269.
[17]吴文欢，何小娟,苏跃龙等于氮氧稳定同位素识别水体氮源的研究进展[J]环境科学与技术2016,39(8)-77-84.
[18] HOPKINS J B IERGUSON J M. Estimating the diets of animals using stable isotopes and a comprehensive Bayesian mixing model[J].PloS one ,2012.7(1):e28478.
[19]冯建祥，黄敏参黄茜，等稳定同位素在滨海湿地生态系统研究中的应用现状与前景[J].生态学杂志2013,32(4): 1065-1074.
[20] PHILLIPS D L,INGER R,BEARHOP S,et al.Best practices for use of stable isotope mixing models in food-web studies[J].Canadian Journal of Zoology,2014.92:823-835.
[21] HOULTON B Z,SIGMAN D M.HEDIN L 0. lsotopic evidence for large gaseous nitrogen losses from tropical rainforests[J] Proceedings of the National Academy of Sciences ofthe United States of America,2006, 103(23)-8745-8750.
[22] PANNO S V,HACKLEY K C,KELLY W R,et al.sotopic Evidence of nitrate sources and denitification in the Mississippi Rivellinois[J] Journal of Environmental Quality,2006,35(2):495-504.
[23] MENGIS M,SCHIF S L.HARRIS M.et al. Multiple geochemical and isotopic approaches for assessing ground water NO3- elimination in a riparian zone[J]. Groundwater, 999,37(3):448-457.
[24] FUKADA T,HISCOCK K M,DENNIS P F,et al.A dual isotope approach to identify denitrification in groundwater at a river-bank infiltration site[J].Water Research,2003,37(13):3070-3078.
[25] FUKADA T,HISCOCK K M,DENNIS P F.A dual-isotope approach to the nitrogen hydrochemistry of an urban aquifer[J] Applied Geochemisty,2004,19(5)-709-719.
[26] LEHMANN M F,REICHERT P,BERNASCONI S M,et al.Modelling nitrogen and oxygen isotope fractionation during denitification in a lacustrine redox-ransition zone[J] Geochimica et Cosmochimica Acta,2003,67(14):2529-2542.
[27] GRANGER J,SIGMAN D M,.LEHMANN M F,et al.Nitrogen and oxygen isotope fractionation during dissimilatory nitrate reduction by dentrifying bacteria[J] Limnology and Oceanography,2008,53(6)-.2533-2545.
[28] WANG Ang,FANG Yunting,CHEN Dexiang,et al.High nitrogen isotope fractionation of nitrate during denitrification in four forest soils and its implications for denitrification rateestimates[J]. Science of The Total Environment,2018,633:1078- 1088.
[29] HU Kelin,HUANG Yuangfang,LI Hong,et al .Spatial variability of shallow groundwater level, electrical conductivity and nitrate concentration,and risk assessment of nitrate contamination in North China Plain[J] Environment International,2005,31(6):896-903.
[30]张鑫,张妍.毕直磊，等.中国地表水硝酸盐分布及其来源分析[].环境科学,2020.41(4): 1594-1606.
[31]赵解春,李玉中,ICHIJI Y,等地下水硝酸盐污染来源的推断与溯源方法概述[J]中国农学通报,2010,26(18):374-378.
[32]徐志伟,张心昱，于贵瑞等中国水体硝酸盐氮氧双稳定同位素溯源研究进展[J]环境科学,2014,35(8):3230-3238.
[33]张千千，缪丽萍孙继朝,等稳定同位素识别水体硝酸盐污染来源的研究进展[J].长江流域资源与环境,2015,24(5):742-749.
[34] MIN J H,YUN S T,KIM K,et al.Nitrate contamination of lluvial groundwaters in the Nakdong River basin,Korea[J]. Geosciences Journal,2002,6(1):35-46.
[35] NYILITYA B,MUREITHI S,BOECKX P. Tracking sources and fate of groundwater nitrate in Kisumu City and Kano Plains ,Kenya[J]. Water,2020,12(2):401.
[361 WANG Yajun,PENG Jianfeng, CAO Xiaofeng.et al.lsotopic and chemical evidence for nitrate sources and transformation processes in a plateau lake basin in Southwest China
[37] WANG Peng,LIU Junzheng,Ql Shuhua,et al.Tracing sources of nitrate using water chemistry,land use and nitrogen isotopes in the Ganjiang River,China[J].Isotopes in Environmental and Health Studies,2017,53(5):539-551.
[38]袁丽君,刘广,张泽洲，等漏湖氮污染双同位素溯源与清单统计法对比研究[J].土壤2018,50(4):738-745.
[39] DEUTSCH B,MEWES M,LISKOW I,et al.Quantification of difuse nitrate inputs into a small river system using stable isotopes of oxygen and nitrogen in nitrate[J].Organic Geochemistry,2006,37(10):1333-1342.
[40]杨锐婧，冯民权汪银龙.汾河下游丰水和枯水期的河流硝酸盐污染来源特征[J.水土保持通报,2019,39(6)-211-217.
[41]王希欢，杨芳，驼娟等乌梁素海硝酸盐来源的季节性变化[J]环境科学研究,2021,34(5):1091-1098.
[42]康萍萍，许士国，禹守泉同位素溯源解析地下水库对地下水氮分布影响[J].水资源保护,2016,32(05):79-84+90.
[43]卢丽,李文莉,裴建国等基于IsoSource的桂林寨底地下河硝酸盐来源定量研究[J].地球学报,2014.35(02)-248-254.
[44]徐璐，蒋勇军,段世辉等于双同位素δ 15N-NO3-5180-NO3-和IsoSource模型的岩溶槽谷区地下水硝酸盐来源的定量示踪:以重庆龙凤和龙车岩溶地下河系统为例[J]环境科学.2020.41(8):3637-3645.
[45]盛婷基于氮氧同位素和IsoSource模型的农业区地下河硝酸盐来源研究[D].重庆西南大学,2019:41-50.
[46]刘坤华北集约农区地下水氮素来源及影响因素分析[D].北京:中国农业大学2016:30-40.
[47]傅雪梅孙源媛苏婧等基于水化学和氮氧双同位素的地下水硝酸盐源解析[J]中国环境科学2019,39(09):3951-3958.
[48]邢萌，刘卫国.浐河、灞河硝酸盐端元贡献比例一基于硝酸盐氮、 氧同位素研究[J] 地球环境学报,2016,7(01):27-36.
[49] WANG Meng,LU Baohong,WANG Jianqun,et al. Using dual isotopes and a Bayesian isotope mixing model to evaluate nitrate sources of surface water in a drinking water source watershed,East China[J].Water,2016,8(8)-355.
[50] YANG Liping.,HAN Jiangpei,XUE Jianlong,et al. Nitrate source apportionment in a subtropical watershed using Bayesian model [J] Science of The Total Environment,2013.463-464:340-347.
[51]张妍,张秋英,李发东等基于稳定同位素和贝叶斯模型的引黄灌区地下水硝酸盐污染源解析[J]中国生态农业学报(中英文).2019,27(03):484-493.
[52]彭聪,潘晓东焦友军,等多种同位素手段的硝酸盐污染源解析:以会仙湿地为例J]环境科学,2018,39(12):5410-5417.
[53] XUE Dongmei,DE BAETS B,VAN CL EEMPUT O,et al.Calsification of nitrate polluting activities through clustering of isotope mixing model outputs[J]Jourmal of Environmental Quality,2013,42(5): 1486-1497.
[54] DIVERS M,ELLIOTT E,BAIN D.Quantification of nitrate sources to an urban stream using dual nitrate isotopes[J] Environmental science and technology,2014,48(18):10580-
[55]MATIATOS 1.Nitrate source identification in groundwater of multiple land-use areas by combining isotopes and multivariate statistical analysis:.A case study of Asopos basin (CentralGreece)[J]. Science of The Total Environment,2016 541:802-814.
[56]KIM K H,YUN S T,MAYER B,et al.Quantification of nitrate sources in groundwater using hydrochemical and dual isotopic data combined with a Bayesian mixing model[J]. Agriculture, Ecosystems and Environment,2015, 199:369-381.
[57] RYBCZYNSKI S,WALTERS D,FRITZ K.et al.Comparing trophic position of stream fishes using stable isotope and gut contents analyses[J] Ecology of Freshwater Fish,2008,17:199-206.
[58]BOECKLEN W J,YARNES C T,COOK B A,et al.On the use of stable isotopes in trophic ecology[J] Annual Review of Ecology ,Evolution, and Systematics,2011 ,42(1):411-440.
[59] POST D M.Using stable isotopes to estimate trophic position: Models, methods ,and assumptions[J]. Ecology,2002,83(3):703-718.
[60] MCCUTCHAN JR J H,LEWIS JR W M.KENDALL C,et al.Variation in trophic shift for stable isotope ratios of carbon,nitrogen,and sulfur[J] Oikos,2003,102(2):378-390.
[61] ZANDEN M J V,RASMUSSEN J B. Variation in δ15N and δ13C trophic fractionation:lmplications for aquatic food web studies([J]. Limnology and Oceanography,2001.46(8).2061-2066.
[62]万神,胡建英,安立会，等利用稳定氮和碳同位素分析渤海湾食物网主要生物种的营养层次[J].科学通报2005,50(7):708-712.
[63]魏虎进,朱小明,纪雅宁,等于稳定同位素技术的象山港海洋牧场区食物网基础,与营养级的研究[].应用海洋学学报2013,32(2):250-257.
[64]蔡德陵,李红燕唐启升,等黄东海生态系统食物网连续营养谱的建立:来自碳氮稳定同位素方法的结果[J].中国科学C辑:生命科学, 2005,35(2)-123-130.
[65]李斌徐丹丹，王志坚，等三峡库区小江库湾鱼类食物网的稳定C、N同位素分析[J] 生态学报，2013,33(20):6704-6711.
[66]周正,米武娟许元钊,等.克氏原螫虾两种养殖模式的食物网结构及其食性比较[J]水生生物学报,2020,44(01):133-142.
[67]王荦，杜双成杨婷越,等.应用稳定同位素技术评价大连近岸海域食物网营养结构[J]生态学杂志,2017 36(05):1452-1457.
[68]杨国欢,侯秀琼孙省利,等.流沙湾食物网结构的初探:基于稳定同位素方法的分析结果[J]水生生物学报.2013,37(01):150-156.
[69]李红燕稳定碳、氮同位素在生态系统中的应用研究:以无定河、黄东海生态系统为例[D] 青岛中国海洋大学2004:60-66.
[70] CABANA G,RASMUSSEN J B Comparison of aquatic food chains using nitrogen isotopes[J]. Proceedings of the National Academy of Sciences of the United States of America,1996, ,93(20):10844-10847.
[71]李云凯，贡艺基于碳、氮稳定同位素技术的东太湖水生食物网结构[J]生态学杂志2014,33(06): 1534-1538.
[72]宋博，陈琳琳,闫朗,等山东东营和烟台潮间带海草床食物网结构特征[J].生物多样性,2019,27(9):984-992.
[73]麻秋云，韩东燕,刘贺，等应用稳定同位素技术构建胶州湾食物网的连续营养谱[J]生态学报2015,35(21):7207-7218.
[74]蒋日进,章守字，王凯,等枸杞岛近岸海域食物网的稳定同位素分析[J].生态学杂志，2014.33(4):930-938.
[75] CHEN Qian,LIU Yang.HO W T,et al.Use of stable isotopes to understand food webs in Macao wetlands[J] Wetlands Ecology and Management,2017,25(1):59-66.
[76]余婕.刘敏侯立军，等崇明东滩大型底栖动物食源的稳定同位素示踪[J]自然资源学报2008,23(2):319-326.
[77] ABRANTES K,SHEAVES M.Food web structure in a near-pristine mangrove area of the Australian Wet Tropics([J]. Estuarine,Coastal and Shelf Science ,2009,82(4):597-607.
[78] LAYMAN C A.ARRINGTON D A.MONTARA C G,et al.Can stable isotope ratios provide for community-wide measures of trophic structure?[J] Ecology,2007 88():42-48.
[79]张文博,黄洪辉李纯厚,等华南典型海湾主要渔业生物碳氮稳定同位素研究[J]南方水产科学,2019,15(5):9-14.
[80]谢斌,李云凯张虎等基于稳定同位素技术的海州湾海洋牧场食物网基础及营养结构的季节性变化[J] .应用生态学报,2017,28(7):2292-2298.
[81]盖珊珊赵文溪.宋静静,等小黑山鸟人工鱼礁区许氏平触和大泷六线鱼的营养生态位研究[J].生态学报,2019,39(18);6923-6931.
[82]俞雅文徐东坡.运用碳氮稳定同位素技术探究中华绒螯蟹与无齿螳臂相手蟹的营养生态位特征[J].长江流域资源与环境2019,28(9);.2092-2099.
[83] JACKSON A,INGER R,PARNELL A.et al.Comparing isotopic niche width among and within communities: SIBER - Stable Isotope Bayesian Elipses in R[J].The Jourmal of animal ecology,2011. 80:595-602.
[84]李云凯,贡艺,陈新军稳定同位素技术在头足类摄食生态学研究中的应那[J].应用生态学报,2014,25(5):1541-1546.
[85] KIDD K A,SCHINDLER D W,HESSLEIN R H.et al.Correlation between stable nitrogen isotope ratios and concentrations of organochlorines in biota from a freshwater food web[J].Science of The Total Environment,1995, 160/161:381-390.
[86] HOBSON K A,VAN WIL GENBURG S L,WASSENAAR L I,et al.A multi-isotope (613C,615N,δ2H) feather isoscape to assign Afrotropical migrant birds to origins[J]. Ecosphere,2012,3(5):article44.
[87]张波，袁伟戴芳群应用稳定同位素技术研究崂山湾夏季鱼类群落的摄食生态[J]水产学报,2016,40(4)-585-594.
[88]崔莹，吴旭干,赵峰，等基于碳、氮稳定同位素分析中华绒蛰蟹(Eriocheir sinensis)幼蟹溯河洄游期食物来源变化[J]海洋与湖沼，2019.50(4):822-829.
[89]陈展彦，武海涛王云彪,等于稳定同位素的湿地食物源判定和食物网构建研究进展[J]应用生态学报,2017 ,28(7).2389-2398.
[90] HARRIGAN P,ZIEMAN J C,MACKO S. The base of nutritional support for the gray snapper (L utjanus Griseus):An evaluation based on a combined stomach content and stableisotope analysis[J] Bulletin of Marine Science, 198944:65-77.
[91] WHITLEDGE G,RABENI C F.Energy sources and ecological role of crayfishes in an Ozark stream:Insights from stable isotopes and gut analysis[J] Canadian Journal of Fisheries and Aquatic Sciences, 1997 ,54:2555-2563.
[92]全为民长江口盐沼湿地食物网的初步研究:稳定同位素分析[D].上海县大学2007:71-89.
[93]鹿志创,田甲申,王召会等应用碳氮稳定同位素技术研究江豚(Neophocanaasiaerienalsssp sunameri食性J].生态学报2016,36(1):69-76 LU Zhichuang.TIAN Jiashen,WANG Zhaohui,et al.Using stable isotope technique to study feeding habits of the finless porpoise(Neophocaenaasiaerientalsssp sunameri)[J]. Acta Ecologica Sinica,2016,36(1):69-76.