一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法
技术领域：
：1.本发明属于农业面源污染监控
技术领域：
：，具体涉及一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法。
背景技术：
：：2.农田过量施用氮磷肥料与土壤侵蚀流失，已经成为全球河流、湖泊等地表水体面源污染的最主要贡献(ongleyetal.,2010；sunetal.,2012)。目前，大多数研究主要集中于农田尺度上可溶性氮磷流失和源头化肥减施面源污染防控方面，但对集约化农业生产驱动的土壤侵蚀对河湖水体面源污染的贡献知之甚少(fissoreetal.,2017；sarapatkaetal.,2018)。自上世纪90年代末，国内外学者开始研究水土流失引起的农业面源污染。研究的空间尺度已从径流小区观测发展到农田、乃至流域尺度，研究的理念已从单一的农田氮磷流失研究开始向坡面流失源区与河湖水库汇集区响应关系的研究发展(collinsetal.,2009；wangetal.,2018)。研究手段从氮磷浓度的离散型观测、地球化学溯源，开始向连续型的智能物联网在线监测、多环境同位素专一性溯源方面发展。3.目前，泥沙溯源辨识的方法主要有地球化学元素、frn、多环境同位素专一性溯源等方法。地球化学元素指纹技术是通过不同母质的重金属组成不同来判断不同土壤泥沙入河来源，但其不能辨识不同土地利用类型泥沙来源；散落的环境放射性核素(falloutradionuclides，frn)技术利用137cs、210pbex和7be具有不同半衰期(30年、22年和53天)，确定不同土地利用类型的土壤侵蚀速率，用于辨析河库泥沙污染物沉积的时间和泥沙的侵蚀类型来源，但不能辨识不同土地利用类型泥沙入河污染来源贡献；特定单体化合物稳定性同位素技术(compound-specificstableisotope，cssi)可以确定不同土地利用类型入河泥沙污染物的贡献百分比,其原理利用不同植物产生相同的脂肪酸作为生物标志物标记土壤，通过分析不同植物的脂肪酸稳定性同位素(δ13c)值辨识流域出口输出不同来源的泥沙贡献比例(gibbs,2008)，但不能确定流域不同土地利用类型的入河泥沙污染负荷。4.流域实时监测技术是一种通过监测降雨条件下流域出口径流和泥沙输出量来反映流域内部泥沙输出实时情况的经典方法(lietal.,2020)，但只能量化泥沙及相关养分入河湖负荷，但无法对入河泥沙来源贡献进行辨识。5.上述这些技术，都是分开应用很少联合应用，难于定量流域不同土地利用类型的入河泥沙污染负荷不能辨识不同土地利用类型泥沙来源。目前已有学者将frn-cssi联合溯源示踪解析流域泥沙污染物的土地利用来源及其输入水系的负荷通量(于寒青等，2019)，但因其技术原理存在问题(frn包括137cs、210pb和7be，分别用于测定过去50年、100年和数月的平均土壤侵蚀速率。cssi测定是次降雨事件不同土地利用类型侵蚀泥沙入河污染的贡献的百分比，所以137cs和210pb测定的长时间平均土壤侵蚀速率与cssi确定的短期降雨入河泥沙时间不匹配。尽管7be可以测定次降雨事件侵蚀速率，但只局限于面蚀，不包括沟蚀过程。因此，frn与cssi联用技术原理存在问题)，故联合frn-cssi示踪技术无法对暴雨条件下农田氮磷流失进入地表水系的通量过程进行实时定量解析。技术实现要素：6.为了解决上述技术问题，本发明提出了一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法，能简单、全面、准确地定量解析不同土地利用侵蚀土壤氮磷的入河污染负荷，为水生态治理提供科学依据。7.本发明的技术方案为：8.一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法，包括以下步骤：9.s1.实时监测流域出口的径流泥沙中氮磷含量，确定流域中泥沙氮磷总入河负荷量；10.s2.应用cssi朔源技术对流域中的入河泥沙来源进行辨识，计算不同土地利用类型泥沙来源的贡献百分比；11.s3.联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷。12.进一步地，所述步骤s1包括：13.采集入河口悬浮泥沙样品，将泥沙过滤后，一部分105℃后烘干用于计算径流中泥沙浓度，另一部分40℃烘干用于测定泥沙氮磷浓度，泥沙氮、泥沙磷分别采用凯式定氮法和hclo4-h2so4-钼蓝比色法进行测定；根据式ⅰ和式ⅱ计算泥沙氮磷总入河负荷量；[0014][0015]式ⅰ中，sci是泥沙浓度(gl-1)，si是1l水沙混合物中烘干泥沙的质量(g)；[0016][0017]式ⅱ中，pnl/ppl是泥沙氮磷总入河负荷量(kgha-1)，n是样本总数，i是次降雨事件中的单个样本数，sci是泥沙浓度(gl-1)，vi是不同样本期的径流量(m3h-1),ti是不同样本数的采样间隔时间(h)，pni/ppi是每次降雨事件中泥沙氮磷的浓度(gkg-1)；ax是子流域面积(ha)。[0018]进一步地，所述步骤s2包括：[0019]s21.流域cssi泥沙样品的采集[0020]a.不同土地利用类型样品采集：降雨径流发生后，在不同土地利用类型集水区采集0-2cm的移动泥沙样品；[0021]b.入河泥沙样品采集：降雨径流发生期间，采集河口悬浮泥沙样品；[0022]s22.室内分析[0023]采集的cssi样品一部分经酸化去除无机碳后，直接用元素分析仪测定土壤δ13c丰度；另一部分使用二氯甲烷萃取、皂化后调至酸性，使用无水硫酸钠去除水分后，再采用三氟化硼甲醇溶液将脂肪酸甲酯化后溶解在正己烷溶液中，然后利用气相色谱-红外质谱系统(gc-irms)用于分析脂肪酸甲酯化合物燃烧转化成co2中的13c值，得到脂肪酸中的δ13c丰度；[0024]s23.泥沙来源贡献比例计算[0025]选择合适的脂肪酸来确定泥沙中不同土地利用来源的碳同位素比例；基于不同土地利用类型筛选的脂肪酸δ13c和bulkδ13c的结果，应用iso-source混合模型，确定泥沙混合物中不同土地利用来源的碳同位素比例；将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用泥沙来源贡献，计算公式如式ⅰ所示：[0026][0027]式ⅲ中，sn％:泥沙第n个源土壤贡献，in:使用同位素混合模型估算的混合物中第n个源土壤的同位素比例(％)，cn％:第n个源土壤有机碳含量。[0028]进一步地，步骤s3中，所述不同土地利用侵蚀土壤氮磷的入河污染负荷根据公式ⅳ计算：[0029]源n的泥沙氮磷的入河负荷(t·ha-1)＝泥沙氮磷总入河负荷量(pnl/ppl，kgha-1)×源n对流域入河泥沙贡献的百分数sn％(％)ⅳ。[0030]本发明的有益效果为：[0031]本发明通过实时监测流域出口的径流泥沙，确定流域中泥沙氮磷总的入河负荷；利用cssi朔源技术对流域中的入河泥沙污染来源进行辨识，确定不同泥沙污染来源的贡献比例；联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷。本发明的优点是能简单、全面、准确地定量解析不同土地利用侵蚀土壤氮磷的入河污染负荷，为水生态治理提供科学依据，可有效防控输入河湖库的泥沙负荷。附图说明[0032]图1为本发明的一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法的流程图。具体实施方式[0033]下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。[0034]实施例1[0035]以某流域为例，在一次暴雨事件下，采集入河口悬浮泥沙样品，将泥沙过滤后，一部分105℃后烘干用于计算径流中泥沙浓度，另一部分40℃烘干用于测定泥沙氮磷含量，泥沙氮、泥沙磷分别采用凯式定氮法和hclo4-h2so4-钼蓝比色法进行测定，结合流域实时监测计算泥沙氮磷总入河负荷量，结果如表1所示。[0036]表1入河总的泥沙氮磷负荷[0037][0038]在雨后采集流域中来自农地、林地、河岸以及道路集水区0-2cm的移动泥沙样品，同时在暴雨过程中，在流域入河口采集悬浮泥沙样品。将这些样品带回实验室处理分析，cssi样品一部分经酸化去除无机碳后，直接用元素分析仪测定土壤δ13c丰度；另一部分使用二氯甲烷萃取、皂化后调至酸性，使用无水硫酸钠去除水分后，再采用三氟化硼甲醇溶液将脂肪酸甲酯化后溶解在正己烷溶液中，然后利用气相色谱-红外质谱系统(gc-irms)用于分析脂肪酸甲酯化合物燃烧转化成co2中的13c值，得到脂肪酸中的δ13c丰度(gibbsetal,2008)。将cssi源-汇样品进行分析，通过选择合适的脂肪酸来确定泥沙中不同土地利用来源的碳同位素比例；基于不同土地利用类型筛选的脂肪酸δ13c和bulkδ13c的结果，应用iso-source混合模型(gibbs,2008)确定泥沙混合物中不同土地利用来源的碳同位素比例；将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用类型泥沙来源贡献，结果如表2所示。[0039]表2不同土地利用类型来源泥沙的贡献百分比[0040][0041]通过实时监测技术确定流域出口入河总的泥沙氮磷负荷，利用cssi朔源技术得出不同土地利用类型来源泥沙的贡献百分比，联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷，结果如表3所示。[0042]表3某流域不同土地利用类型侵蚀土壤氮磷的入河污染负荷[0043][0044]上述说明是针对发明较佳可行实施例的详细说明，但实施例并非用以限定本发明的专利申请范围，凡本发明所提示的技术精神下所完成的同等变化或修饰变更，均应属于本发明所涵盖专利范围。当前第1页12当前第1页12
技术领域：
：1.本发明属于农业面源污染监控
技术领域：
：，具体涉及一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法。
背景技术：
：：2.农田过量施用氮磷肥料与土壤侵蚀流失，已经成为全球河流、湖泊等地表水体面源污染的最主要贡献(ongleyetal.,2010；sunetal.,2012)。目前，大多数研究主要集中于农田尺度上可溶性氮磷流失和源头化肥减施面源污染防控方面，但对集约化农业生产驱动的土壤侵蚀对河湖水体面源污染的贡献知之甚少(fissoreetal.,2017；sarapatkaetal.,2018)。自上世纪90年代末，国内外学者开始研究水土流失引起的农业面源污染。研究的空间尺度已从径流小区观测发展到农田、乃至流域尺度，研究的理念已从单一的农田氮磷流失研究开始向坡面流失源区与河湖水库汇集区响应关系的研究发展(collinsetal.,2009；wangetal.,2018)。研究手段从氮磷浓度的离散型观测、地球化学溯源，开始向连续型的智能物联网在线监测、多环境同位素专一性溯源方面发展。3.目前，泥沙溯源辨识的方法主要有地球化学元素、frn、多环境同位素专一性溯源等方法。地球化学元素指纹技术是通过不同母质的重金属组成不同来判断不同土壤泥沙入河来源，但其不能辨识不同土地利用类型泥沙来源；散落的环境放射性核素(falloutradionuclides，frn)技术利用137cs、210pbex和7be具有不同半衰期(30年、22年和53天)，确定不同土地利用类型的土壤侵蚀速率，用于辨析河库泥沙污染物沉积的时间和泥沙的侵蚀类型来源，但不能辨识不同土地利用类型泥沙入河污染来源贡献；特定单体化合物稳定性同位素技术(compound-specificstableisotope，cssi)可以确定不同土地利用类型入河泥沙污染物的贡献百分比,其原理利用不同植物产生相同的脂肪酸作为生物标志物标记土壤，通过分析不同植物的脂肪酸稳定性同位素(δ13c)值辨识流域出口输出不同来源的泥沙贡献比例(gibbs,2008)，但不能确定流域不同土地利用类型的入河泥沙污染负荷。4.流域实时监测技术是一种通过监测降雨条件下流域出口径流和泥沙输出量来反映流域内部泥沙输出实时情况的经典方法(lietal.,2020)，但只能量化泥沙及相关养分入河湖负荷，但无法对入河泥沙来源贡献进行辨识。5.上述这些技术，都是分开应用很少联合应用，难于定量流域不同土地利用类型的入河泥沙污染负荷不能辨识不同土地利用类型泥沙来源。目前已有学者将frn-cssi联合溯源示踪解析流域泥沙污染物的土地利用来源及其输入水系的负荷通量(于寒青等，2019)，但因其技术原理存在问题(frn包括137cs、210pb和7be，分别用于测定过去50年、100年和数月的平均土壤侵蚀速率。cssi测定是次降雨事件不同土地利用类型侵蚀泥沙入河污染的贡献的百分比，所以137cs和210pb测定的长时间平均土壤侵蚀速率与cssi确定的短期降雨入河泥沙时间不匹配。尽管7be可以测定次降雨事件侵蚀速率，但只局限于面蚀，不包括沟蚀过程。因此，frn与cssi联用技术原理存在问题)，故联合frn-cssi示踪技术无法对暴雨条件下农田氮磷流失进入地表水系的通量过程进行实时定量解析。技术实现要素：6.为了解决上述技术问题，本发明提出了一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法，能简单、全面、准确地定量解析不同土地利用侵蚀土壤氮磷的入河污染负荷，为水生态治理提供科学依据。7.本发明的技术方案为：8.一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法，包括以下步骤：9.s1.实时监测流域出口的径流泥沙中氮磷含量，确定流域中泥沙氮磷总入河负荷量；10.s2.应用cssi朔源技术对流域中的入河泥沙来源进行辨识，计算不同土地利用类型泥沙来源的贡献百分比；11.s3.联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷。12.进一步地，所述步骤s1包括：13.采集入河口悬浮泥沙样品，将泥沙过滤后，一部分105℃后烘干用于计算径流中泥沙浓度，另一部分40℃烘干用于测定泥沙氮磷浓度，泥沙氮、泥沙磷分别采用凯式定氮法和hclo4-h2so4-钼蓝比色法进行测定；根据式ⅰ和式ⅱ计算泥沙氮磷总入河负荷量；[0014][0015]式ⅰ中，sci是泥沙浓度(gl-1)，si是1l水沙混合物中烘干泥沙的质量(g)；[0016][0017]式ⅱ中，pnl/ppl是泥沙氮磷总入河负荷量(kgha-1)，n是样本总数，i是次降雨事件中的单个样本数，sci是泥沙浓度(gl-1)，vi是不同样本期的径流量(m3h-1),ti是不同样本数的采样间隔时间(h)，pni/ppi是每次降雨事件中泥沙氮磷的浓度(gkg-1)；ax是子流域面积(ha)。[0018]进一步地，所述步骤s2包括：[0019]s21.流域cssi泥沙样品的采集[0020]a.不同土地利用类型样品采集：降雨径流发生后，在不同土地利用类型集水区采集0-2cm的移动泥沙样品；[0021]b.入河泥沙样品采集：降雨径流发生期间，采集河口悬浮泥沙样品；[0022]s22.室内分析[0023]采集的cssi样品一部分经酸化去除无机碳后，直接用元素分析仪测定土壤δ13c丰度；另一部分使用二氯甲烷萃取、皂化后调至酸性，使用无水硫酸钠去除水分后，再采用三氟化硼甲醇溶液将脂肪酸甲酯化后溶解在正己烷溶液中，然后利用气相色谱-红外质谱系统(gc-irms)用于分析脂肪酸甲酯化合物燃烧转化成co2中的13c值，得到脂肪酸中的δ13c丰度；[0024]s23.泥沙来源贡献比例计算[0025]选择合适的脂肪酸来确定泥沙中不同土地利用来源的碳同位素比例；基于不同土地利用类型筛选的脂肪酸δ13c和bulkδ13c的结果，应用iso-source混合模型，确定泥沙混合物中不同土地利用来源的碳同位素比例；将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用泥沙来源贡献，计算公式如式ⅰ所示：[0026][0027]式ⅲ中，sn％:泥沙第n个源土壤贡献，in:使用同位素混合模型估算的混合物中第n个源土壤的同位素比例(％)，cn％:第n个源土壤有机碳含量。[0028]进一步地，步骤s3中，所述不同土地利用侵蚀土壤氮磷的入河污染负荷根据公式ⅳ计算：[0029]源n的泥沙氮磷的入河负荷(t·ha-1)＝泥沙氮磷总入河负荷量(pnl/ppl，kgha-1)×源n对流域入河泥沙贡献的百分数sn％(％)ⅳ。[0030]本发明的有益效果为：[0031]本发明通过实时监测流域出口的径流泥沙，确定流域中泥沙氮磷总的入河负荷；利用cssi朔源技术对流域中的入河泥沙污染来源进行辨识，确定不同泥沙污染来源的贡献比例；联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷。本发明的优点是能简单、全面、准确地定量解析不同土地利用侵蚀土壤氮磷的入河污染负荷，为水生态治理提供科学依据，可有效防控输入河湖库的泥沙负荷。附图说明[0032]图1为本发明的一种实时监测与cssi技术解析泥沙来源氮磷入河负荷的方法的流程图。具体实施方式[0033]下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。[0034]实施例1[0035]以某流域为例，在一次暴雨事件下，采集入河口悬浮泥沙样品，将泥沙过滤后，一部分105℃后烘干用于计算径流中泥沙浓度，另一部分40℃烘干用于测定泥沙氮磷含量，泥沙氮、泥沙磷分别采用凯式定氮法和hclo4-h2so4-钼蓝比色法进行测定，结合流域实时监测计算泥沙氮磷总入河负荷量，结果如表1所示。[0036]表1入河总的泥沙氮磷负荷[0037][0038]在雨后采集流域中来自农地、林地、河岸以及道路集水区0-2cm的移动泥沙样品，同时在暴雨过程中，在流域入河口采集悬浮泥沙样品。将这些样品带回实验室处理分析，cssi样品一部分经酸化去除无机碳后，直接用元素分析仪测定土壤δ13c丰度；另一部分使用二氯甲烷萃取、皂化后调至酸性，使用无水硫酸钠去除水分后，再采用三氟化硼甲醇溶液将脂肪酸甲酯化后溶解在正己烷溶液中，然后利用气相色谱-红外质谱系统(gc-irms)用于分析脂肪酸甲酯化合物燃烧转化成co2中的13c值，得到脂肪酸中的δ13c丰度(gibbsetal,2008)。将cssi源-汇样品进行分析，通过选择合适的脂肪酸来确定泥沙中不同土地利用来源的碳同位素比例；基于不同土地利用类型筛选的脂肪酸δ13c和bulkδ13c的结果，应用iso-source混合模型(gibbs,2008)确定泥沙混合物中不同土地利用来源的碳同位素比例；将泥沙中的有机碳含量和模型输出碳同位素比例转换成不同土地利用类型泥沙来源贡献，结果如表2所示。[0039]表2不同土地利用类型来源泥沙的贡献百分比[0040][0041]通过实时监测技术确定流域出口入河总的泥沙氮磷负荷，利用cssi朔源技术得出不同土地利用类型来源泥沙的贡献百分比，联合应用实时监测与cssi泥沙朔源，定量确定不同土地利用侵蚀土壤氮磷的入河污染负荷，结果如表3所示。[0042]表3某流域不同土地利用类型侵蚀土壤氮磷的入河污染负荷[0043][0044]上述说明是针对发明较佳可行实施例的详细说明，但实施例并非用以限定本发明的专利申请范围，凡本发明所提示的技术精神下所完成的同等变化或修饰变更，均应属于本发明所涵盖专利范围。当前第1页12当前第1页12
再多了解一些
本文用于企业家、创业者技术爱好者查询，结果仅供参考。
