湖泊科学   2022, Vol. 34 Issue (5): 1413-1427.  DOI: 10.18307/2022.0501
0

综述

引用本文 [复制中英文]

李影, 秦丽欢, 雷秋良, 罗加法, 杜新忠, 闫铁柱, 刘宏斌, 小流域农业面源污染监测断面设置与污染物通量估算研究进展. 湖泊科学, 2022, 34(5): 1413-1427. DOI: 10.18307/2022.0501
[复制中文]
Li Ying, Qin Lihuan, Lei Qiuliang, Luo Jiafa, Du Xinzhong, Yan Tiezhu, Liu Hongbin. Review on agricultural non-point source pollution monitoring sections layout and pollutant loading estimation in small watershed. Journal of Lake Sciences, 2022, 34(5): 1413-1427. DOI: 10.18307/2022.0501
[复制英文]

基金项目

国家自然科学基金项目(42007095,U20A20114,42177325)和中国农业科学院科技创新工程协同创新项目(CAAS-GXAAS-XTCX2019026-03)联合资助

通信作者

雷秋良, E-mail:leiqiuliang@caas.cn

文章历史

2021-11-16 收稿
2022-01-14 收修改稿

码上扫一扫

小流域农业面源污染监测断面设置与污染物通量估算研究进展
李影1,2,3 , 秦丽欢4 , 雷秋良1 , 罗加法5 , 杜新忠1 , 闫铁柱1 , 刘宏斌1     
(1: 中国农业科学院农业资源与农业区划研究所, 农业农村部面源污染控制重点实验室, 北京 100081)
(2: 中国科学院地理科学与资源研究所, 资源与环境信息系统国家重点实验室, 北京 100101)
(3: 中国科学院大学,北京 100049)
(4: 北京市生态环境监测中心,北京 100048)
(5: 新西兰农业科学院,汉密尔顿 3240)
摘要:流域尺度面源污染的监测是系统认识农业面源污染的发生、迁移及转化过程, 并对其进行有效控制的重要基础. 当前, 田块尺度的面源污染监测方法比较成熟, 而流域尺度的监测, 尤其是监测断面的布设及采样频率设置等方面的研究较少. 本文详细梳理了小流域采样断面布设、采样频率优化和河流断面通量估算3方面的主要进展. 1)从小流域样点布设来看, 包括常规监测采样设计、针对性监测采样设计和融合前两种类型的监测采样设计3种类型, 点位布设方法上有遗传算法、模糊逻辑法、熵值法和模型法等, 样品采集方式包括随机采样、复合采样、综合采样以及连续采样4种类型, 其中复合采样应用广泛;2)从采样频率来看, 1~2周一次的采样频率即可精准获取污染负荷通量, 若需要进一步提高精度, 可在水文/水质变异大的时期提高监测频率及在特殊断面加密布点;3)在通量估算上, 当前主要的计算方法有平均法、插值法和回归/曲线法3类方法, 其中流量加权的浓度估计法、插值算法和LOADEST法是简便且精确的方法, 方法的选择上也可根据不同时期流域污染源特征进一步优化.
关键词小流域    面源污染    监测点布设    监测频率    通量估算    
Review on agricultural non-point source pollution monitoring sections layout and pollutant loading estimation in small watershed
Li Ying1,2,3 , Qin Lihuan4 , Lei Qiuliang1 , Luo Jiafa5 , Du Xinzhong1 , Yan Tiezhu1 , Liu Hongbin1     
(1: Key Laboratory of Nonpoint Source Pollution Control, Ministry of Agriculture and Rural Affairs, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, P.R.China)
(2: State Key Laboratory of Resources and Environmental Information System, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, P.R.China)
(3: University of Chinese Academy of Sciences, Beijing 100049, P.R.China)
(4: Beijing Municipal Ecological and Environmental Monitoring Center, Beijing 100048, P.R.China)
(5: AgResearch Limited, Ruakura Research Centre, Hamilton 3240, New Zealand)
Abstract: Watershed-scale non-point source pollution monitoring is an important basis for systematically understanding the occurrence, transportation and transformation processes of agricultural non-point source pollution and its effective control. Presently methods monitoring non-point source pollution at field scale are relatively mature, while methods applied at watershed scale, particularly the layout of monitoring section and sampling frequency setting, are rarely studied. This paper reviewed the main progress of sampling section layout, sampling frequency optimization and loading estimation of river sections. From the perspective of sampling layout in small watershed, there are three types of monitoring design including probabilistic design, targeted design and both. Sampling points were located with genetic algorithm, fuzzy logic method, entropy and modelling method, etc. Sample collection methods include random sampling, compound sampling, comprehensive sampling and continuous sampling. Compound sampling is widely used. Concerning the sampling frequency, once every 1-2 weeks is sufficient to obtain accurate pollution loading. If higher accuracy is necessary, the monitoring frequency and special sections can be increased when hydrology/water quality variation is large. For loading estimation, average method, interpolation method and regression or curve method were widely used. Concentration estimation algorithm by flow-weighted, interpolation and LOADEST are simple and exact. Selection of methods could be further optimized based on the characteristics of pollution sources in different periods of the basin.
Keywords: Small watershed    non-point source pollution    layout of monitoring points    monitoring frequency    loading estimation    

农业面源污染指在农业生产和农村生活区域, 氮、磷等营养盐及其他污染物受水力驱动以随机、分散、无组织方式进入受纳水体引起的水质恶化[1]. 农业面源污染对水体污染的贡献率在世界大部分国家超过了50 % [2]. 我国作为农业大国, 同样面临农业面源污染这一重大环境问题[3]. 两次全国污染源普查数据显示, 2007年我国农业源总氮、总磷排放量分别占排放总量的57.2 % 和67.4 %, 到2017年数值有所降低, 但仍占到46.5 % 和67.2 % [4-5], 这说明农业面源是我国当前乃至未来一段时期的重要污染源.

监测、规划、实施、评估等一系列系统性方法的综合使用是农业面源污染有效控制的前提和基础, 其中监测是极其重要的一环. 各国近年来开展了许多水质监测项目用于流域管理. 这些项目的制定一般遵循以下几个目的:1)摸清当前的水质现状[6]; 2)解析点源和面源污染[6]; 3)为流域模型提供数据支持[7-8]; 4)在确定边界条件和监测时间的前提下对河流的长期污染趋势进行评估[9-10]; 5)评估最佳管理措施(best management practices, BMPs)的有效性, 为政策和管理措施的制定提供指导[9-11].

相较于森林/草地主导的流域, 农业小流域受灌溉和施肥等农业活动的影响, 流域内河流的水文水质有较高的时空变异性[12]. 此外, 在农业主导的流域, 施肥期往往和高流量期(晚春季节或早夏)耦合, 导致季节性的通量估算偏差变大[13-16]. 在管道排水的流域(如城市和现代农牧区)[17-19], 由于雨污的收集设施较为完善, 可将样点布设在管道节点及出口处. 而传统农业主导的流域, 由于其污染物的广域分散性、运输途径的多样无序性和水文水质的高时空变异性, 采样断面的布设需要考虑的因素更为复杂[20]. 《流域农业面源污染监测技术规范》(NY/T 3824—2020)中规定了流域农业面源污染监测技术的监测断面设置与采样、监测指标及方法、流域农业面源污染的结果表达与质量控制等要求[1]. 尽管目前面源污染的控制取得了极大的进步, 但如何使面源管控措施兼顾经济性和有效性, 还需要进一步的研究.

河流断面污染物浓度测量及通量估算可以表征各类污染物对河流水质的影响, 是各类模型校验的主要数据源, 也是进行流域主要污染源和污染区识别的直接依据. 流域研究方案的制定, 是在了解流域基本情况的基础上, 对流域水质采样断面的位置、采样的时间和频率进行确定[21]. 本文从农业小流域水质监测的采样点的布设、采样频率的确定和通量估算方法的选择入手, 梳理了国内外研究进展, 旨在为我国的小流域农业面源污染监测设定和最佳通量估算方法的选择提供参考.

1 采样位置和采样方式的确定

河流断面监测是获取流域水质基本特征, 制定管理决策的重要基础手段[22]. 具有代表性的样本采集, 对水质的准确表征和管理至关重要, 在流域决策中也具有重要意义[23-24]. 采样位置的选择可以影响监测的结果, 是能否达到监测目标的决定因素[25]. 水样的采集和分析需要巨大的经济投入, 这是我国河流水质监测点布设数量的主要限制因素[26]. 如何在资金预算有限的前提下, 合理布设水质监测断面, 优化各监测指标的采样频率, 是亟需解决的问题.

1.1 流域监测类型

流域监测可根据监测需求划分为3种类型:第1类是常规监测采样设计或概要设计(probabilistic design), 即采用随机的方式在流域进行采样点设置, 该方式主要为了能够对流域有较为公正客观的认识和评价; 第2类是针对性监测设计(targeted design), 此种方式是较为常用的设计方式, 是根据对具体问题的认知或即将发生事件的认知进行采样点设计的监测类型, 比如研究最佳管理措施的应用效果; 第3类是融合前两种类型的监测采样设计(combination of targeted and probabilistic). 针对性监测设计的主要分类及优缺点的比较见表 1[27].

表 1 监测设计的分类及特点[27] Tab. 1 Classification and characteristics of monitoring design[27]
1.2 样点布设方法

自1970s以来, 地表水质监测网的布设便受到了广泛关注, 针对水质监测断面的设置也开展了很多研究[28]. 采样断面一般分为背景断面和控制断面两类:背景断面应布设在基本未受人类活动干扰的河流源头位置; 控制断面应布设于顺直河段、河床稳定、水流集中、无浅滩处, 避开死水区、回水区、排水口处[1]. 背景断面的布设考虑因素较为简单, 但控制断面的布设需要对流域特征和监测目的进行综合考虑, 样点布设的合理性是后续进行数据分析的前提[28-29], 如果采样点的代表性不足,采样频率、数据分析和呈现的特征就变得无关紧要[23].

最初水质样本的采集只考虑便捷性, 如在桥下布设采样点[30], 随着地理信息系统技术和数学算法的发展, 水质采样点的布设方法有了很大的进步, 如遗传算法[31]、模糊逻辑法[32]、熵值法[33]和模型法[34]等. 下面介绍几种常用的水质采样点设置和优化方法.

1.2.1 河网拓扑法

该方法将流域上游最小的不分枝支流指定为一级支流, 将仅由一级支流组成的河流定义为二级支流, 从而确定河流的等级. 河流等级确定后, 通过寻找质心的方式找到采样点的位置, 并通过计算流域出口处河网的连接点数目获得流域采样点个数[35-36]. Asadi等考虑了每个采样点的潜在污染,并将此方法在伊朗进行了应用[37]. 该方法可有效地识别潜在的污染源[38], 简单且对历史水质数据的需求较低, 适用于河流分枝不多的大流域, 如果分枝的数量过多(如分枝超过100时)或没有分枝时, 采用此种方法会存在困难[23-24].

1.2.2 模糊逻辑法

该方法主要根据模糊函数计算的值对监测站进行排序, 选取耕地面积百分比、建成面积百分比、面源污染量、绿化覆盖率、滑坡面积比、山坡过度利用率、水质监测站密度等指标, 采用加权法对子流域进行评分, 在得分较高的子流域增加水质监测站, 该理论结合多准则分析提高了数据分类和排序的客观性[32]. 该方法在中国台湾、伊朗等地区得到了很好的应用[39]. 模糊理论在复杂决策过程的情况下具有以下优势:1)理论规则简单, 易解释和理解; 2)是一种高级分类方法, 通过模糊隶属函数处理信息的不确定性; 3)适用于需要通过区分等级之间的实际差异来对优先级进行排名的决策过程[39]. 模糊方法被认为比传统的分类模型(如自然断裂法)更客观、更准确, 为了节省预算, 在模糊理论法指定为一级的子流域内新建水质监测站更为可行[32].

1.2.3 遗传算法

该方法定义了适应度函数, 该函数为河流系统的代表性、水质标准、用水监督、污染源监测和水质变化检查5个对开发监测系统至关重要的标准的线性组合. 利用地理信息数据对适应度函数进行一系列计算, 得出适合度等级, 并对世代数、种群大小、交叉和变异概率等主要参数进行敏感性分析, 以确定良好的适应度水平和最优解的收敛性. Park等[31]利用遗传算法和地理信息系统相结合的方法, 设计了一个有效的大型河流水质监测网络. 虽然用遗传算法得到的监测网的优化设计优于传统的主观设计[22], 但还需要进一步的研究来定义适应度函数, 使得这些函数可以代表选择标准之间的综合关系或优先级.

1.2.4 熵值理论

该方法利用熵值理论, 在考虑系统中最大不确定性和最小冗余信息的前提下, 选择监测网中采样点和采样频率的最佳组合. 先利用河流水质模拟模型生成河流沿岸几个潜在监测点的水质数据, 再利用熵理论对含有冗余信息的台站进行剔除. 该方法在伊朗南部的Karoon河进行了应用[33], 可以方便地对现有的监测网络进行评价和修正, 但要求现有采样网络有足够数据, 以获得具有统计显著性的结果[40].

1.2.5 模型法

TWSST(targeted watershed site selection tool)模型由美国威斯康辛州自然资源部研发[34], 该模型进行监测点位布设的原理是先利用秩相关检验和典型相关找到水质和参数的关系, 进而通过K值聚类分析对河网进行分组, 使得每组的各水质参数绝对差异化, 从而得到流域监测点位的分布. 该模型可以采用先验方法进行有效的流域监测设计, 也可以通过后验方法确定监测位置所代表的空间范围.

河流物理、化学以及生物特性受到自然和人为因素的干扰, 流域的监测断面位置和监测断面密度需要能够反映这些自然和人为因素, 这就要求设计者对流域内影响因素的分布和未来的变化趋势有清晰的了解[23, 41]. 因此, 流域采样点的空间布设, 需要结合地理和人为要素, 充分利用地统计分析方法的优势, 结合多变量统计方法和优化算法, 对采样点位置进行确定[38].

1.3 采样方式

采样的方式可以分为4种类型:随机采样、复合采样、综合采样以及连续采样. 其中, 随机采样指的是在某个地方一次性采样, 并不适用于需要考虑采样对象的空间异质性的采样, 如分层的湖泊、水质随时间的变化(如降雨事件的影响); 复合采样方式是一系列随机采样的集合, 可以分为时间加权和流量加权两种采样类型, 时间加权指在固定时间间隔进行采样, 流量加权指流量到一定的峰值进行采样; 综合采样指考虑水体水质空间变化, 将不同位置的采样点混合后作为原始样品的采样方式; 连续采样则需借助自动监测设备进行, 一般用于溶解氧、pH、电导率等指标的监测. 表 2列出了不同采样方式能实现的研究目标, 从表中可以看出复合采样适合大多数研究[27, 42].

表 2 采样方式及其主要特征 Tab. 2 Sampling methods and its main characteristics

不同的研究目标对监测的频率、时间和数据分析强度的要求不同. 监测频率受研究对象的影响较大, 如措施的有效性监测需要比趋势监测采用更高的采样频次, 而趋势监测则需要更长的监测时间[43]. 表 2列出了不同监测类型所需要的监测频率、持续时间和数据分析强度[42].

1.4 采样频率的确定与调整

在监测项目设计过程中, 设计者对监测位置的布设给予了高度关注, 但对监测频率的关注却并不充分, Nguyen等评估了59项研究, 只有1项研究考虑了监测频率和监测成本[44]. 当前, 流量已可以用现场仪器连续或近连续测量(水位-流量传感器可以通过水位和流量的特定关系将测量数据转换为流量数据)[45-46], 测量精度也较高. 例如, 美国地质勘探局的7个测站的流量的测量效果较好, 测量误差仅为5 %~10 % [47]. 而水质指标的浓度测量由于人力和财力的限制大部分还是离散采样(如1次/周或1次/月)[48-50], 这导致在通量估算过程中, 样品收集带来的不确定性较高, 达到了累计不确定性的±4 %~±50 % [51]. 虽然有研究表明采样频率越高, 对污染物动态变化过程的描述越精确[52-53], 但是高频率水质监测需要耗费较多的人力物力, 往往不太现实. Kovács等研究发现将采样频率从1次/d降低到1次/3 d, 可以节省大约50 % 的成本[54]. 因此, 合理地设置采样频率可以减少样品采集和测试中产生的成本消耗. 国内外对采样频率对通量估计影响的研究总结见表 3.

表 3 采样频率和算法对流域通量估算影响研究总结 Tab. 3 Summary of the influence of sampling frequency and algorithm on fluxes estimation
1.4.1 采样频率的确定

虽有部分学者提出采用时间序列自相关分析结合数据分布特征来设定较为现实合理的监测频率[49, 55], 但多数仍是采用偏差分析的方法比较不同采样频率相对于最高采样频率下平均浓度的误差值, 从而得到满足一定误差范围的合理采样频率[52, 56]. 通常采样频率降低时, 检测具有统计学意义趋势的能力会下降[57], 通量估计的准确度也会随着采样频率的降低而降低, 并且采样频率的降低对不同的水质指标造成的误差结果不同[58]. 河流水文监测的不确定性易受水文特征和流域特征的影响. 水文特征包括:水文响应关系、浓度-流量相关关系、基流指数(base flow index,BFI)和变率指数(flashiness index,FI)等[49, 59-61]; 流域特征包括人口密度、流域尺度等[62-63]. 因此, 采样频率的设定取决于流域特征[62], 尤其是水文条件[49]和浓度行为[64].

对Loire河的研究结果表明, 当采样频率定为1次/月时, 总氮(TN)年通量估算的均方根误差(RMSE)范围为6.6 %~20.4 %, 硝态氮(NO3--N)为4 %~26 %、溶解态活性磷(SRP)为5 %~33 %、颗粒态磷(PP)为6 %~18 %, 总磷(TP)为13 %~42 %; 此外, TN的年通量估算精度(standard deviation, SD)为8.8 %, NO3--N的SD为13 %, TP的SD为20 %, SRP的SD为26 %, PP的SD为34 % [65]. 因此, 对于不同指标, 最佳采样频率的设定略有不同. 在大型流域, 采样频率对误差的影响不大, 但对精度有近线性影响, 若要达到10 % 精度, NO3--N采样间隔为15天, SRP和TP采样间隔为10天, SS采样间隔为5天[65]; 而在管道排水的流域, SRP的频率设置为13~26 h、NO3--N设置为2.7~17.5天, 可以达到±10 % 的准确性[66]; 相较于森林流域, 农业流域的水文水质时空变异性较高, 因此准确地进行通量估算更为困难[12], 在面积为217 km2的山地农业小流域, 1次/2周的采样频率可以使TN年通量的RMSE控制在5 % 左右, 但这个频率却不足以分析通量的周期或季节特征(4个季节的平均RMSE高达17 %)[67].

1.4.2 采样频率的优化

浓度与流量的复杂关系将会对通量估算产生影响[12, 68], 在高流量事件中, 流量和浓度之间通常存在强相关关系[16, 64], 水位上升期间, 物质浓度可能达到最大, 在第一次冲刷之后则有可能产生稀释效应, 浓度逐渐降低[69-70], 大量悬浮沉积物的运输可能只发生在高流量时期[71], 高流量事件贡献了营养物的主要年通量[46, 72]. 因此, 基于水文变异的采样方案对通量估算来说是最准确的[63], 并且河流监测采样应尽可能多地捕捉高流量事件, 至少累积涵盖年流量的80 %~85 % [46, 73-74], 这就要求在暴雨径流等事件发生时进行加密观测[56].

相较于自然植被覆盖的流域, 农业流域由于灌溉和施肥等农业活动的影响, 河流水文水质的时空变异性较高[12], 且存在特定的季节特征, 导致季节性的通量估算变异变大[13-16]. 此外, 在较小的流域或者流域上游水文响应关系较强, 水文特征的时空变异较大, 通量估算的偏差也较大[49, 62, 75-77]. 在这类站点, 低频率的采样不能捕获短时内发生的水文的快速响应[61, 75, 78].

由于变异因素的存在, 如研究区域独特的特征和特殊的河流水文状况、不同的算法、不同的污染物和研究持续时间等, 现有研究结果很难被普适化[48, 79]. 对于时间自相关性强的污染物序列, 其浓度变化稳定, 可在保证精度的前提下, 减少频率来降低成本; 而对于污染物时间变异性较大的序列, 若要反映真实的水质变化过程, 则需要提高采样频率[49, 55]. 因此, 需要及时进行数据分析, 对于时间序列变化相对不稳定的指标进行加密监测和频率调整. 《流域农业面源污染监测技术规范》(NY/T 3824—2020)中规定了加密监测时段, 分别为农业面源污染剧增期(耕作、施肥等农事活动密集, 且发生径流排水并引起控制断面水位明显变化的时段)和人口剧增时段(如重要节假日或旅游季节). 此外, 规范中还规定汛期监测频率应高于非汛期[1].

2 断面通量估算方法

估算算法对通量估计的影响不如采样频率大, Vries和Klavers[102]在荷兰一个研究区的研究表明, 河流污染物的通量估算首先应考虑采样策略, 其次是算法. 针对当前连续/高频的流量数据和离散/低频的水质数据, 主要形成了3类计算通量的方法:平均法、插值法和回归/曲线法[58, 62, 82, 103-104]. 平均法使用实测数据的平均值作为时段数据的代表值; 插值法假定数据(浓度和流量)由瞬时采样数据表示; 回归(外推)方法取决于经验关系和水文水质的变化特征[65].

2.1 常用估算方法

河流监测断面的污染物通量是断面流量和浓度的函数, 污染物通量的变化与降雨-产流过程有很强的相关关系[50, 105]. 监测断面污染物浓度与流量有正相关、负相关和无关3种关系[106]. 当前各学者采用实测断面平均浓度、样品时间平均浓度、断面瞬时流量、采样期间平均流量、采样代表时段平均流量等概念构造了多种监测断面通量估算方法. 在此, 对应用较多的8种通量估算方法进行了分析比较, 各种方法的公式及物理意义如表 4所示.

表 4 通量估算方法* Tab. 4 Fluxes estimation methods

从环境水力学角度对表中的估算公式进行分析, 可将监测断面的流量及平均浓度表达为时间平均的形式[93-94]

$ Q(t)=\frac{1}{T} \int_T Q(t) \mathrm{d} t+Q^{\prime \prime}(t)=Q_{\mathrm{a}}+Q^{\prime \prime}(t) $ (1)
$ C(t)=\frac{1}{T} \int_T C(t) \mathrm{d} t+C^{\prime \prime}(t)=C_{\mathrm{a}}+C^{\prime \prime}(t) $ (2)

式中, Qa为时段平均流量; Ca为时段平均浓度; T为估算时段; Q″和C″分别为流量和浓度时均距平值. 时段通量可以通过式(3)表达:

$ W=\int_T F(t) \mathrm{d} t=\int_T Q(t) \cdot C(t) \mathrm{d} t=Q_a \cdot C_a \cdot T+\int_T Q^{\prime \prime} \cdot C^{\prime \prime} \mathrm{d} t \approx \sum\limits_{i=1}^n Q_i \cdot C_i \Delta t_i= \\ Q_{\mathrm{a}} \cdot C_{\mathrm{a}} \cdot T+\sum\limits_{i=1}^n Q^{\prime \prime}{ }_i \cdot C_i^{\prime \prime} \Delta t_i $ (3)

在8种算法中(表 4), 算法A和D使用时间加权平均浓度值, 在高流量时具有稀释效果; 而算法B和C使用部分流量加权浓度数据; 算法E、F和G则是应用通量的基本定义[110]. 这些算法通过不同的方程来关联浓度和流量数据, 估计给定时间段内的通量. 将表 4中各方法与式(3)进行对比, 可以看到方法A、D实际上只包含了(3)中的第1项(对流项), 而忽略了第2项(时均离散项), 在使用过程中会反映出时间加权平均浓度对流量加权平均浓度的高估或低估; 方法A和D的差别在于方法A采用的是离散的实测流量平均, 方法D采用的是时段的平均流量; 方法B、C、E、F和G则两项都包括, 方法B和C的差别在于采用离散的流量平均还是连续的流量平均; 方法E则采用了流量加权平均浓度, 与时段平均流量的乘积作为时段通量; 方法F采用插值方法, 将离散浓度数据插值为高频浓度数据, 与监测获得的高频流量数据的乘积作为时段通量; 方法G利用浓度与流量的相关关系, 用流量预测浓度, 结合监测获得的高频流量数据得到时段通量; 方法H利用连续的流量监测和离散的浓度监测, 用流量和时间预测瞬时通量, 瞬时通量加和得到时段通量. 显然, 方法A~G均是先得到对应的流量和浓度数据, 再进行通量的估算; 方法H是直接以通量与流量或事件的回归关系得到瞬时通量, 要求流量数据具有连续性; 方法A、D仅适用于推断流域断面(断面流速均匀)的平均污染通量的估算; 方法B适合点源占优的情况; 方法C、E适合面源贡献较大的通量估算; 方法F适用污染物排放短时变异不强的情况; 方法G和H要求流量与浓度或通量具有较好的回归关系(表 5).

表 5 时段通量估算方法的应用取向分析[94] Tab. 5 Application orientation analysis of the estimation methods for period-fluxes[94]
2.2 估算方法的选择

由于不同方法的计算结果存在较大差异, 目前对估算方法的适用性研究是当前流域污染物通量研究的热点之一[110-112]. 研究表明, 由于污染物类型的差别, 在统计时段通量时(如年通量)对于不同的污染物应采用不同的估算方法[95, 97]. 在A~F这7种方法中, 基于流量的浓度估计算法(E)和插值算法(F)是大多数水质参数通量估算的首选算法[47, 65-66, 91]. 虽然对于NO3--N通量的估算采用方法E较为合适[58, 66], 但方法F的变异系数最小[47]. 同时, 方法F对于TN和TP的估算效果好于方法E[82]. 方法D对于TP(尤其是PP)的估算效果最好[65], 但在高流量时可能会由于稀释效应产生高估[82, 91]. 总体来看, 线性插值法(F)的表现效果是7种方法中最好的[65-66]. 由于方法H所采用的LOADEST模型在计算时会运用Akaike信息准则(AIC)和Schwarz后验概率准则(SPPC)对模型进行优化、检验, 因此, 在预测通量和观测趋势方面表现良好, 但普遍倾向于高估[113]. 基于方法H的通量估算的研究多采用低频的水质数据, 并且采样频率对结果误差的影响较小[100]. 此外, 对于G和H这类回归模型的估算效果不仅随水质参数的不同而不同, 也强烈依赖降雨事件样品所占的比例[114].

通常情况下, 使用污染物的年通量来评估当前水质的年内状况是无法实现的[6], 季节通量的估算往往不如年度通量的估算准确[60]. 污染源类型和丰/枯水期的不同也会影响估算方法的表现效果[92, 111], 夏季发生的暴雨事件会使污染物浓度产生更大的变异性, 导致通量估算的不确定性大于其他3个季节[66-67]. 此外, 在面源占优的情况下, 丰水期的监测采用较高的监测频率对提高长时间通量的估算精度是必要的[94, 115]. 因此, 有研究将丰水期和枯水期分别进行分析, 发现采用时段通量平均浓度与时段平均流量之积(方法E), 分别计算丰、平、枯季的污染通量, 然后相加求和来估算年通量, 是最简便精确的方法[93]; 也有学者通过比较不同水文阶段污染物通量估算结果, 提出枯水期采用以点源占优的方法B, 丰水期采用强调时段总径流量以及以面源占优的方法E, 平水期则结合相关的水文特征进行方法的选择[92]; 如果采样频率低且污染物输出与流量的相关性较好, 方法H是可以保证估算精准度的首选方法[100]. 此外, 也可以根据不同时期流域污染的点源和面源特征进行方法的选择[94, 96-97](表 5).

流域土地利用特征也会对通量估算效果产生影响. 由于农业活动具有明显的时空变异特征, 农业流域污染物输出的时空变异性也更高, 因此, 准确地进行通量估算要比森林流域更具挑战性[12, 16]; 此外, 通量估算方法的表现效果也受采样频率的影响, 一般采样频率越低, 表现效果会越差[67]. 由此可以看出导致通量测算结果不确定性的因素, 不仅包括通量估算方法, 还包括水文特征[80, 105]、污染源特征(点源或面源类型)[94, 115]、采样策略[116-117]、采样时间[98]和流域大小[82]等因素.

3 结语

由于经济和地理因素, 当前我国农业面源污染监测断面设置数量较少, 尤其是偏远地区, 监测断面数量和监测频率有待提高, 监测断面布设和采样策略的制定也亟需科学化和标准化. 经济欠发达地区往往对先进技术和设备的投入不足, 这也在一定程度上影响到监测工作的质量和水平. 随着3S技术、信息技术和化学分析技术的发展, 我国农业面源污染监测在朝着标准化、自动化和智能化的方向发展. 未来还需建立科学的监测体系和标准, 加强人员培训, 改造和引进监测设备和技术, 促进监测站的标准化、规范化和自动化运营.

本文通过文献回顾的方式, 就农业小流域采样点布设、采样频率优化和断面通量估算方法3个主要环节进行了详细的分析和介绍. 虽然各流域的监测方案应在充分考虑监测目的、流域特征、水文特征和人类活动的前提下制定, 但通过文献回顾仍能发现一些共性的规律, 可为今后的监测项目提供参考.

4 参考文献

[1]
Technical specification for agricultural non-point source pollution monitoring at the watershed scale. NY/T 3824-2020. [流域农业面源污染监测技术规范. NY/T 3824—2020. ]
[2]
Agency EP ed. National water quality inventory: Report to congress. Washington, DC (United States): Environmental Protection Agency, 2017.
[3]
Xue BL, Zhang HW, Wang GQ et al. Evaluating the risks of spatial and temporal changes in nonpoint source pollution in a Chinese River Basin. Science of the Total Environment, 2022, 807: 151726. DOI:10.1016/j.scitotenv.2021.151726
[4]
Ministry of Ecology and Environment of the People's Republic of China, National Bureau of Statistics, Ministry of Agriculture and Rural Affairs of the People's Republic of China eds. Bulletin of the second national survey of pollution sources. https://www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202006/t20200610_783547.html, 2020-06-09. [中华人民共和国生态环境部, 国家统计局, 中华人民共和国农业农村部. 第二次全国污染普查公报. ]
[5]
Ministry of Ecology and Environment of the People's Republic of China, National Bureau of Statistics, Ministry of Agriculture and Rural Affairs of the People's Republic of China eds. Bulletin of the second national survey of pollution sources. https://www.mee.gov.cn/gkml/hbb/bgg/201002/t20100210_185698.htm, 2010-02-06. [环境保护部, 国家统计局, 农业部. 第一次全国污染普查公报. ]
[6]
Roygard J, McArthur K, Clark M. Diffuse contributions dominate over point sources of soluble nutrients in two sub-catchments of the Manawatu River, New Zealand. New Zealand Journal of Marine and Freshwater Research, 2012, 46(2): 219-241. DOI:10.1080/00288330.2011.632425
[7]
Gassman PW, Reyes MR, Green CH et al. The soil and water assessment tool: Historical development, applications, and future research directions. Transactions of the ASABE, 2007, 50(4): 1211-1250. DOI:10.13031/2013.23637
[8]
Ullrich A, Volk M. Influence of different nitrate-N monitoring strategies on load estimation as a base for model calibration and evaluation. Environmental Monitoring and Assessment, 2010, 171(1/2/3/4): 513-527. DOI:10.1007/s10661-009-1296-8
[9]
Brauer N, O'Geen AT, Dahlgren RA. Temporal variability in water quality of agricultural tailwaters: Implications for water quality monitoring. Agricultural Water Management, 2009, 96(6): 1001-1009. DOI:10.1016/j.agwat.2009.01.011
[10]
King KW, Harmel RD. Considerations in selecting a water quality sampling strategy. Transactions of the ASAE, 2003, 46(1): 63-73. DOI:10.13031/2013.12549
[11]
Lennartz B, Tiemeyer B, de Rooij G et al. Artificially drained catchments—from monitoring studies towards management approaches. Vadose Zone Journal, 2010, 9(1): 1-3. DOI:10.2136/vzj2009.0149
[12]
Kelly PT, Vanni MJ, Renwick WH. Assessing uncertainty in annual nitrogen, phosphorus, and suspended sediment load estimates in three agricultural streams using a 21-year dataset. Environmental Monitoring and Assessment, 2018, 190(2): 91. DOI:10.1007/s10661-018-6470-4
[13]
Richards RP, Baker DB. Trends in water quality in LEASEQ rivers and streams (northwestern Ohio), 1975-1995. Journal of Environmental Quality, 2002, 31(1): 90-96. DOI:10.2134/jeq2002.9000
[14]
Royer TV, David MB, Gentry LE. Timing of riverine export of nitrate and phosphorus from agricultural watersheds in Illinois: Implications for reducing nutrient loading to the Mississippi River. Environmental Science & Technology, 2006, 40(13): 4126-4131. DOI:10.1021/es052573n
[15]
Renwick WH, Vanni MJ, Zhang QY et al. Water quality trends and changing agricultural practices in a midwest US watershed, 1994-2006. Journal of Environmental Quality, 2008, 37(5): 1862-1874. DOI:10.2134/jeq2007.0401
[16]
Carpenter SR, Booth EG, Kucharik CJ et al. Extreme daily loads: Role in annual phosphorus input to a north temperate lake. Aquatic Sciences, 2015, 77(1): 71-79. DOI:10.1007/s00027-014-0364-5
[17]
Steidl J, Lischeid G, Engelke C et al. The curse of the past—What can tile drain effluent tell us about arable field management?. Agriculture, Ecosystems & Environment, 2022, 326: 107787. DOI:10.1016/j.agee.2021.107787
[18]
Barkle G, Stenger R, Moorhead B et al. The importance of the hydrological pathways in exporting nitrogen from grazed artificially drained land. Journal of Hydrology, 2021, 597: 126218. DOI:10.1016/j.jhydrol.2021.126218
[19]
Liu WL, Youssef MA, Birgand FP et al. Processes and mechanisms controlling nitrate dynamics in an artificially drained field: Insights from high-frequency water quality measurements. Agricultural Water Management, 2020, 232: 106032. DOI:10.1016/j.agwat.2020.106032
[20]
Wang YG, Wang HY, Zheng YL et al. Advances in research methods and control technologies of agricultural non-point source pollution: A review. Chinese Journal of Agricultural Resources and Regional Planning, 2021, 42(1): 25-33. [王一格, 王海燕, 郑永林等. 农业面源污染研究方法与控制技术研究进展. 中国农业资源与区划, 2021, 42(1): 25-33. DOI:10.7621/cjarrp.1005-9121.20210104]
[21]
Agency EP ed. Monitoring guidance for determining the effectiveness of nonpoint source controls. Washington, DC: Environmental Protection Agency, 2016.
[22]
Behmel S, Damour M, Ludwig R et al. Water quality monitoring strategies—A review and future perspectives. Science of the Total Environment, 2016, 571: 1312-1329. DOI:10.1016/j.scitotenv.2016.06.235
[23]
Alilou H, Moghaddam Nia A, Keshtkar H et al. A cost-effective and efficient framework to determine water quality monitoring network locations. Science of the Total Environment, 2018, 624: 283-293. DOI:10.1016/j.scitotenv.2017.12.121
[24]
Do HT, Lo SL, Chiueh PT et al. Design of sampling locations for mountainous river monitoring. Environmental Modelling & Software, 2012, 27. DOI:10.1016/j.envsoft.2011.09.007
[25]
da Luz N, Kumpel E. Evaluating the impact of sampling design on drinking water quality monitoring program outcomes. Water Research, 2020, 185: 116217. DOI:10.1016/j.watres.2020.116217
[26]
Yang X, Liu Q, Luo X et al. Spatial regression and prediction of water quality in a watershed with complex pollution sources. Scientific Reports, 2017, 7: 8318. DOI:10.1038/s41598-017-08254-w
[27]
Dressing SA, Meals DW, Harcum JB et al eds. Monitoring and evaluating nonpoint source watershed projects. Washington, DC: U.S. Environmental Protection Agency, 2016.
[28]
Karamouz M, Kerachian R, Akhbari M et al. Design of river water quality monitoring networks: A case study. Environmental Modeling & Assessment, 2008, 14(6): 705-714. DOI:10.1007/s10666-008-9172-4
[29]
Sanders TG, Ward RC, Loftis JC et al. Design of networks for monitoring water quality. Colorado: Water Resouces Publications, 1983.
[30]
Strobl RO, Robillard PD. Network design for water quality monitoring of surface freshwaters: A review. Journal of Environmental Management, 2008, 87(4): 639-648. DOI:10.1016/j.jenvman.2007.03.001
[31]
Park SY, Choi JH, Wang S et al. Design of a water quality monitoring network in a large river system using the genetic algorithm. Ecological Modelling, 2006, 199(3): 289-297. DOI:10.1016/j.ecolmodel.2006.06.002
[32]
Chang CL, Lin YT. A water quality monitoring network design using fuzzy theory and multiple criteria analysis. Environmental Monitoring and Assessment, 2014, 186(10): 6459-6469. DOI:10.1007/s10661-014-3867-6
[33]
Karamouz M, Nokhandan AK, Kerachian R et al. Design of on-line river water quality monitoring systems using the entropy theory: A case study. Environmental Monitoring and Assessment, 2009, 155(1/2/3/4): 63-81. DOI:10.1007/s10661-008-0418-z
[34]
Wisconsin Department of Natural Resources, Optimizing a monitoring design in targeted watersheds: The targeted watershed site selection tool. Madison: Wisconsin Department of Natural Resources, 2015.
[35]
Sharp WE. A topologically optimum water-sampling plan for rivers and streams. Water Resources Research, 1971, 7(6): 1641-1646. DOI:10.1029/WR007i006p01641
[36]
Horton RE. Erosional development of streams and their drainage basins: hydrophysical approach to quantitative morphology. Geological Society of America Bulletin, 1945, 56(3): 275-370. DOI:10.1177/030913339501900406
[37]
Asadi A, Moghaddam NA, Bakhtiari EB et al. An integrated approach for prioritization of river water quality sampling points using modified Sanders, analytic network process, and hydrodynamic modeling. Environmental Monitoring and Assessment, 2021, 193(8): 1-15. DOI:10.1007/s10661-021-09272-y
[38]
Jiang JP, Tang SJ, Han DW et al. A comprehensive review on the design and optimization of surface water quality monitoring networks. Environmental Modelling & Software, 2020, 132: 104792. DOI:10.1016/j.envsoft.2020.104792
[39]
Alilou H, Rahmati O, Singh VP et al. Evaluation of watershed health using Fuzzy-ANP approach considering geo-environmental and topo-hydrological criteria. Journal of Environmental Management, 2019, 232: 22-36. DOI:10.1016/j.jenvman.2018.11.019
[40]
Li J, Tian LQ, Wang YH et al. Optimal sampling strategy of water quality monitoring at high dynamic lakes: A remote sensing and spatial simulated annealing integrated approach. Science of the Total Environment, 2021, 777: 146113. DOI:10.1016/j.scitotenv.2021.146113
[41]
Camara M, Jamil NR, Abdullah AFB et al. Economic and efficiency based optimisation of water quality monitoring network for land use impact assessment. Science of the Total Environment, 2020, 737: 139800. DOI:10.1016/j.scitotenv.2020.139800
[42]
MacDonald LH, Smart AW, Wissmar RC. Monitoring guidelines to evaluate effects of forestry activities on streams in the Pacific Northwes and Alaska, 1991.
[43]
Gibson GR, Barbour MT, Stribling JB et al eds. Biological criteria: Technical guidance for streams and small rivers. Washington, DC (United States): Environmental Protection Agency, 1996.
[44]
Nguyen TH, Helm B, Hettiarachchi H et al. The selection of design methods for river water quality monitoring networks: A review. Environmental Earth Sciences, 2019, 78(3): 1-17. DOI:10.1007/s12665-019-8110-x
[45]
Horowitz AJ. An evaluation of sediment rating curves for estimating suspended sediment concentrations for subsequent flux calculations. Hydrological Processes, 2003, 17(17): 3387-3409. DOI:10.1002/hyp.1299
[46]
Horowitz AJ. Determining annual suspended sediment and sediment-associated trace element and nutrient fluxes. Science of the Total Environment, 2008, 400(1/2/3): 315-343. DOI:10.1016/j.scitotenv.2008.04.022
[47]
Jiang Y, Frankenberger JR, Bowling LC et al. Quantification of uncertainty in estimated nitrate-N loads in agricultural watersheds. Journal of Hydrology, 2014, 519: 106-116. DOI:10.1016/j.jhydrol.2014.06.027
[48]
Cassidy R, Jordan P. Limitations of instantaneous water quality sampling in surface-water catchments: Comparison with near-continuous phosphorus time-series data. Journal of Hydrology, 2011, 405(1/2): 182-193. DOI:10.1016/j.jhydrol.2011.05.020
[49]
Jones AS, Horsburgh JS, Mesner NO et al. Influence of sampling frequency on estimation of annual total phosphorus and total suspended solids Loads1. JAWRA Journal of the American Water Resources Association, 2012, 48(6): 1258-1275. DOI:10.1111/j.1752-1688.2012.00684.x
[50]
Pohle I, Baggaley N, Palarea-Albaladejo J et al. A framework for assessing concentration-discharge catchment behavior from low-frequency water quality data. Water Resources Research, 2021, 57(9): e2021WR029692. DOI:10.1029/2021WR029692
[51]
Harmel RD, Cooper RJ, Slade RM et al. Cumulative uncertainty in measured streamflow and water quality data for small watersheds. Transactions of the ASABE, 2006, 49(3): 689-701. DOI:10.13031/2013.20488
[52]
Skarbøvik E, Stålnacke P, Bogen J et al. Impact of sampling frequency on mean concentrations and estimated loads of suspended sediment in a Norwegian River: Implications for water management. Science of the Total Environment, 2012, 433: 462-471. DOI:10.1016/j.scitotenv.2012.06.072
[53]
Lu XX, Li SY, He M et al. Seasonal changes of nutrient fluxes in the Upper Changjiang Basin: An example of the Longchuanjiang River, China. Journal of Hydrology, 2011, 405(3/4): 344-351. DOI:10.1016/j.jhydrol.2011.05.032
[54]
Kovács J, Korponai J, Székely Kovács I et al. Introducing sampling frequency estimation using variograms in water research with the example of nutrient loads in the Kis-Balaton Water Protection System (W Hungary). Ecological Engineering, 2012, 42: 237-243. DOI:10.1016/j.ecoleng.2012.02.004
[55]
Alewell C, Lischeid G, Hell U et al. High temporal resolution of ion fluxes in semi-natural ecosystems-gain of information or waste of resources?. Biogeochemistry, 2004, 69(1): 19-35. DOI:10.1023/b:biog.0000031029.46798.7f
[56]
Li B, Yang GS, Wan RR et al. Temporal variability of water quality in Poyang Lake outlet and the associated water level fluctuations: A water quality sampling revelation. Resources and Environment in the Yangtze Basin, 2017, 26(2): 289-296. [李冰, 杨桂山, 万荣荣等. 鄱阳湖出流水质2004—2014年变化及其对水位变化的响应: 对水质监测频率的启示. 长江流域资源与环境, 2017, 26(2): 289-296. DOI:10.11870/cjlyzyyhj201702015]
[57]
da Luz N, Tobiason JE, Kumpel E. Water quality monitoring with purpose: Using a novel framework and leveraging long-term data. Science of the Total Environment, 2022, 818: 151729. DOI:10.1016/j.scitotenv.2021.151729
[58]
Birgand F, Faucheux C, Gruau G et al. Uncertainties in assessing annual nitrate loads and concentration indicators: Part 1. impact of sampling frequency and load estimation algorithms. Transactions of the ASABE, 2010, 53(2): 437-446. DOI:10.13031/2013.29584
[59]
Moatar F, Meybeck M, Raymond S et al. River flux uncertainties predicted by hydrological variability and riverine material behaviour. Hydrological Processes, 2013, 27(25): 3535-3546. DOI:10.1002/hyp.9464
[60]
Kerr JG, Eimers MC, Yao HX. Estimating stream solute loads from fixed frequency sampling regimes: The importance of considering multiple solutes and seasonal fluxes in the design of long-term stream monitoring networks. Hydrological Processes, 2016, 30(10): 1521-1535. DOI:10.1002/hyp.10733
[61]
Stelzer RS, Likens GE. Effects of sampling frequency on estimates of dissolved silica export by streams: The role of hydrological variability and concentration-discharge relationships. Water Resources Research, 2006, 42(7): W07415. DOI:10.1029/2005WR004615
[62]
Johnes PJ. Uncertainties in annual riverine phosphorus load estimation: Impact of load estimation methodology, sampling frequency, baseflow index and catchment population density. Journal of Hydrology, 2007, 332(1/2): 241-258. DOI:10.1016/j.jhydrol.2006.07.006
[63]
Horowitz AJ, Clarke RT, Merten GH. The effects of sample scheduling and sample numbers on estimates of the annual fluxes of suspended sediment in fluvial systems. Hydrological Processes, 2015, 29(4): 531-543. DOI:10.1002/hyp.10172
[64]
Littlewood IG. Hydrological regimes, sampling strategies, and assessment of errors in mass load estimates for United Kingdom rivers. Environment International, 1995, 21(2): 211-220. DOI:10.1016/0160-4120(95)00011-9
[65]
Moatar F, Meybeck M. Compared performances of different algorithms for estimating annual nutrient loads discharged by the eutrophic River Loire. Hydrological Processes, 2005, 19(2): 429-444. DOI:10.1002/hyp.5541
[66]
Williams MR, King KW, Macrae ML et al. Uncertainty in nutrient loads from tile-drained landscapes: Effect of sampling frequency, calculation algorithm, and compositing strategy. Journal of Hydrology, 2015, 530: 306-316. DOI:10.1016/j.jhydrol.2015.09.060
[67]
Li Y, Yen H, Harmel RD et al. Effects of sampling strategies and estimation algorithms on total nitrogen load determination in a small agricultural headwater watershed. Journal of Hydrology, 2019, 579: 124114. DOI:10.1016/j.jhydrol.2019.124114
[68]
Aulenbach BT, Burns DA, Shanley JB et al. Approaches to stream solute load estimation for solutes with varying dynamics from five diverse small watersheds. Ecosphere, 2016, 7(6): e01298. DOI:10.1002/ecs2.1298
[69]
Vanni MJ, Renwick WH, Headworth JL et al. Dissolved and particulate nutrient flux from three adjacent agricultural watersheds: A five-year study. Biogeochemistry, 2001, 54(1): 85-114. DOI:10.1023/A:1010681229460
[70]
Jiang R, Woli KP, Kuramochi K et al. Hydrological process controls on nitrogen export during storm events in an agricultural watershed. Soil Science and Plant Nutrition, 2010, 56(1): 72-85. DOI:10.1111/j.1747-0765.2010.00456.x
[71]
Gao P, Josefson M. Temporal variations of suspended sediment transport in Oneida Creek watershed, central New York. Journal of Hydrology, 2012, 426/427: 17-27. DOI:10.1016/j.jhydrol.2012.01.012
[72]
Horowitz AJ. A review of selected inorganic surface water quality-monitoring practices: Are we really measuring what we think, and if so, are we doing it right?. Environmental Science & Technology, 2013, 47(6): 2471-2486. DOI:10.1021/es304058q
[73]
Hooper RP, Aulenbach BT, Kelly VJ. The national stream quality accounting network: A flux-based approach to monitoring the water quality of large rivers. Hydrological Processes, 2001, 15(7): 1089-1106. DOI:10.1002/hyp.205
[74]
Gonzalez-Nicolas A, Schwientek M, Sinsbeck M et al. Characterization of export regimes in concentration-discharge plots via an advanced time-series model and event-based sampling strategies. Water, 2021, 13(13): 1723. DOI:10.3390/w13131723
[75]
Coynel A, Schäfer J, Hurtrez JE et al. Sampling frequency and accuracy of SPM flux estimates in two contrasted drainage basins. Science of the Total Environment, 2004, 330(1/2/3): 233-247. DOI:10.1016/j.scitotenv.2004.04.003
[76]
Duvert C, Gratiot N, Némery J et al. Sub-daily variability of suspended sediment fluxes in small mountainous catchments & ndash implications for community-based river monitoring. Hydrology and Earth System Sciences, 2011, 15(3): 703-713. DOI:10.5194/hess-15-703-2011
[77]
Yang TT, Zhang L, Yue Y et al. Optimal estimates for dissolved and suspended particulate material fluxes in the Yangtze River, China. Environmental Science and Pollution Research, 2021, 28(30): 41337-41350. DOI:10.1007/s11356-021-13581-7
[78]
Moatar F, Person G, Meybeck M et al. The influence of contrasting suspended particulate matter transport regimes on the bias and precision of flux estimates. Science of the Total Environment, 2006, 370(2/3): 515-531. DOI:10.1016/j.scitotenv.2006.07.029
[79]
Defew LH, May L, Heal KV. Uncertainties in estimated phosphorus loads as a function of different sampling frequencies and common calculation methods. Marine and Freshwater Research, 2013, 64(5): 373. DOI:10.1071/mf12097
[80]
Richards RP, Holloway J. Monte Carlo studies of sampling strategies for estimating tributary loads. Water Resources Research, 1987, 23(10): 1939-1948. DOI:10.1029/WR023i010p01939
[81]
Rekolainen S, Posch M, Kämäri J et al. Evaluation of the accuracy and precision of annual phosphorus load estimates from two agricultural basins in Finland. Journal of Hydrology, 1991, 128(1/2/3/4): 237-255. DOI:10.1016/0022-1694(91)90140-D
[82]
Kronvang B, Bruhn AJ. Choice of sampling strategy and estimation method for calculating nitrogen and phosphorus transport in small lowland streams. Hydrological Processes, 1996, 10(11): 1483-1501. DOI:10.1002/(SICI)1099-1085(199611)10:111483:AID-HYP386>3.0.CO;2-Y
[83]
Phillips JM, Webb BW, Walling DE et al. Estimating the suspended sediment loads of rivers in the LOIS study area using infrequent samples. Hydrological Processes, 1999, 13(7): 1035-1050. DOI:10.1002/(SICI)1099-1085(199905)13:71035:AID-HYP788>3.0.CO;2-K
[84]
Guo YP, Markus M, Demissie M. Uncertainty of nitrate-N load computations for agricultural watersheds. Water Resources Research, 2002, 38(10): 3-1. DOI:10.1029/2001WR001149
[85]
Moatar F, Meybeck M. Riverine fluxes of pollutants: Towards predictions of uncertainties by flux duration indicators. Comptes Rendus Geoscience, 2007, 339(6): 367-382. DOI:10.1016/j.crte.2007.05.001
[86]
Bowes MJ, Smith JT, Neal C. The value of high-resolution nutrient monitoring: A case study of the River Frome, Dorset, UK. Journal of Hydrology, 2009, 378(1/2): 82-96. DOI:10.1016/j.jhydrol.2009.09.015
[87]
Tiemeyer B, Kahle P, Lennartz B. Designing monitoring programs for artificially drained catchments. Vadose Zone Journal, 2010, 9(1): 14-24. DOI:10.2136/vzj2008.0181
[88]
Birgand F, Appelboom TW, Chescheir GM et al. Estimating nitrogen, phosphorus, and carbon fluxes in forested and mixed-use watersheds of the lower coastal plain of north Carolina: Uncertainties associated with infrequent sampling. Transactions of the ASABE, 2011, 54(6): 2099-2110. DOI:10.13031/2013.40668
[89]
Sergeant CJ, Nagorski S. The implications of monitoring frequency for describing riverine water quality regimes. River Research and Applications, 2015, 31(5): 602-610. DOI:10.1002/rra.2767
[90]
Reynolds KN, Loecke TD, Burgin AJ et al. Optimizing sampling strategies for riverine nitrate using high-frequency data in agricultural watersheds. Environmental Science & Technology, 2016, 50(12): 6406-6414. DOI:10.1021/acs.est.5b05423
[91]
Elwan A, Singh R, Patterson M et al. Influence of sampling frequency and load calculation methods on quantification of annual river nutrient and suspended solids loads. Environmental Monitoring and Assessment, 2018, 190(2): 78. DOI:10.1007/s10661-017-6444-y
[92]
Tan H, Zhao WJ, Wang L. Study on the estimation of pollutant fluxes for Dongliao River based on hydrological periods. Journal of Anhui Agricultural Sciences, 2012, 40(1): 337-339. [谭恒, 赵文晋, 王伦. 东辽河分水期污染物通量估算研究. 安徽农业科学, 2012, 40(1): 337-339. DOI:10.3969/j.issn.0517-6611.2012.01.132]
[93]
Li HE, Li C. Comparision of calculation methods of pollution load for control section of river. Journal of Water Resources and Water Engineering, 2013, 24(2): 1-4. [李怀恩, 李层. 河流控制断面污染负荷计算方法比较. 水资源与水工程学报, 2013, 24(2): 1-4.]
[94]
Fu G. Analysis of the estimation methods for the river pollutant fluxes(Ⅰ): Comparison and analysis of the estimation methods of period fluxes. Research of Environmental Sciences, 2003, 16(1): 1-4. [富国. 河流污染物通量估算方法分析(Ⅰ)——时段通量估算方法比较分析. 环境科学研究, 2003, 16(1): 1-4. DOI:10.13198/j.res.2003.01.3.fug.001]
[95]
Hao CL, Deng YX, Wang YH et al. Study on the selection and error analysis of riverine pollutant flux estimation methods. Acta Scientiae Circumstantiae, 2012, 32(7): 1670-1676. [郝晨林, 邓义祥, 汪永辉等. 河流污染物通量估算方法筛选及误差分析. 环境科学学报, 2012, 32(7): 1670-1676. DOI:10.13671/j.hjkxxb.2012.07.032]
[96]
Wang H. Annual pollutants flux calculation of water quality monitoring sections in mainstream of the Huaihe River. Water Resources Protection, 2004, 20(6): 37-39. [王晖. 淮河干流水质断面污染物年通量估算. 水资源保护, 2004, 20(6): 37-39. DOI:10.3969/j.issn.1004-6933.2004.06.010]
[97]
Wang WP, Hong HS, Zhang YZ et al. Preliminary estimate for the contaminations fluxes from Jiulong River to the sea. Marine Environmental Science, 2006, 25(2): 45-47, 57. [王卫平, 洪华生, 张玉珍等. 九龙江污染物入海通量初步估算. 海洋环境科学, 2006, 25(2): 45-47, 57. DOI:10.3969/j.issn.1007-6336.2006.02.012]
[98]
Wang JL, Jiang GQ, Ren XW et al. Investigation of monitoring strategies estimation methods for watershed pollutant fluxes. Environmental Protection of Xinjiang, 2011, 33(2): 1-7, 18. [王君丽, 姜国强, 任秀文等. 流域污染物通量测算方法研究. 新疆环境保护, 2011, 33(2): 1-7, 18. DOI:10.3969/j.issn.1008-2301.2011.02.001]
[99]
Cai YQ, Tang XY, Zhang JQ et al. Study on optimization of the monitoring frequency of nitrogen and phosphorus loss flux in the river of a small hilly catchment. Environmental Pollution & Control, 2021, 43(5): 614-619, 658. [蔡燕秋, 唐翔宇, 张建强等. 丘陵小流域氮磷流失通量监测频次的优化研究. 环境污染与防治, 2021, 43(5): 614-619, 658.]
[100]
Cheng GW, Du ZP, Yan CA et al. Impacts of water quality monitoring frequency and extreme climate on the estimation of nitrogen and phosphorus fluxes in river of Plateau Lake basin. Acta Scientiae Circumstantiae, 2020, 40(11): 3982-3989. [程国微, 杜展鹏, 严长安等. 水质监测频率与极端气候对高原湖泊入湖河流氮磷通量估算的影响. 环境科学学报, 2020, 40(11): 3982-3989. DOI:10.13671/j.hjkxxb.2020.0155]
[101]
Su JJ, Li XY, Wu Z. Comparison of different methods estimating annual sediment loads in river cross sections based on irregularly measured data. China Environmental Science, 2017, 37(1): 218-228. [苏静君, 李叙勇, 吴震. 利用实测值估算断面泥沙年负荷的方法比较. 中国环境科学, 2017, 37(1): 218-228.]
[102]
Vries AD, Klavers HC. Riverine flues of pollutants: Monitoring strategy first, calculation methods second. European Water Management, 1994, 4(2): 12-17.
[103]
Quilbé R, Rousseau AN, Duchemin M et al. Selecting a calculation method to estimate sediment and nutrient loads in streams: Application to the Beaurivage River (Québec, Canada). Journal of Hydrology, 2006, 326(1/2/3/4): 295-310. DOI:10.1016/j.jhydrol.2005.11.008
[104]
Preston SD, Bierman Jr VJ, Silliman SE. An evaluation of methods for the estimation of tributary mass loads. Water Resources Research, 1989, 25(6): 1379-1389. DOI:10.1029/WR025i006p01379
[105]
Gorski G, Zimmer MA. Hydrologic regimes drive nitrate export behavior in human-impacted watersheds. Hydrology and Earth System Sciences, 2021, 25(3): 1333-1345. DOI:10.5194/hess-25-1333-2021
[106]
Ren Y. Study on surface water quality assessment and pollution load of Ebinur Lake Watershed based on SPARROW model[Dissertation]. Urumqi: Xinjiang University, 2017. [任岩. 基于SPARROW模型的艾比湖流域地表水水质评价及污染负荷研究[学位论文]. 乌鲁木齐: 新疆大学, 2017. ] ]
[107]
Shih G, Abtew W, Obeysekera J. Accuracy of nutrient runoff load calculations using time-composite sampling. Transactions of the Asae American Society of Agricultural Engineers, 1994, 37(2): 419-429. DOI:10.13031/2013.28093
[108]
Littlewood IG. Estimating contaminant loads in rivers: A review. Wallingford: Institute of Hydrology, 1992.
[109]
Littlewood IG, Watts CD, Custance JM. Systematic application of United Kingdom River flow and quality databases for estimating annual river mass loads (1975-1994). Science of the Total Environment, 1998, 210/211: 21-40. DOI:10.1016/S0048-9697(98)00042-4
[110]
Webb BW, Phillips JM, Walling DE et al. Load estimation methodologies for British rivers and their relevance to the LOIS RACS(R) programme. Science of the Total Environment, 1997, 194/19: 379-389. DOI:10.1016/S0048-9697(96)05377-6
[111]
Zhang D, Duan H, Yang HX et al. Study on the estimation of river pollutant flux in rural area. China Biogas, 2015, 33(3): 95-98. [张丹, 段慧, 杨洪霞等. 农村河流断面污染物通量估算研究. 中国沼气, 2015, 33(3): 95-98.]
[112]
Marinović Ruždjak A, Ruždjak D. Evaluation of river water quality variations using multivariate statistical techniques. Environmental Monitoring and Assessment, 2015, 187(4): 1-14. DOI:10.1007/s10661-015-4393-x
[113]
Jha B, Jha MK. Rating curve estimation of surface water quality data using LOADEST. Journal of Environmental Protection, 2013, 4(8): 849-856. DOI:10.4236/jep.2013.48099
[114]
Park YS, Engel BA. Analysis for regression model behavior by sampling strategy for annual pollutant load estimation. Journal of Environmental Quality, 2015, 44(6): 1843-1851. DOI:10.2134/jeq2015.03.0137
[115]
Fu G, Lei K. Analysis of the estimation methods for the river pollutant fluxes(Ⅱ): Error judgment of time-averaged or section-averaged dispersion fluxes. Research of Environmental Sciences, 2003, 16(1): 5-9, 42. [富国, 雷坤. 河流污染物通量估算方法分析(Ⅱ)——时空平均离散通量误差判断. 环境科学研究, 2003, 16(1): 5-9, 42. DOI:10.13198/j.res.2003.01.7.fug.002]
[116]
Aulenbach BT, Hooper RP. The composite method: An improved method for stream-water solute load estimation. Hydrological Processes, 2006, 20(14): 3029-3047. DOI:10.1002/hyp.6147
[117]
Snelder TH, McDowell RW, Fraser CE. Estimation of catchment nutrient loads in New Zealand using monthly water quality monitoring data. JAWRA Journal of the American Water Resources Association, 2017, 53(1): 158-178. DOI:10.1111/1752-1688.12492