Next Article in Journal
Poria Acid, Triterpenoids Extracted from Poria cocos, Inhibits the Invasion and Metastasis of Gastric Cancer Cells
Previous Article in Journal
Exploring the Synthetic Potential of γ-Lactam Derivatives Obtained from a Multicomponent Reaction—Applications as Antiproliferative Agents
Previous Article in Special Issue
The Use of a Solid Bismuth Microelectrode for Vanadium Quantification by Adsorptive Stripping Voltammetry in Environmental Water Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment in Monitoring of Water Analysis of a Brazilian River

1
Postgraduate Programme in Metrology, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 22451-900, Brazil
2
Veiga de Almeida University, Rio de Janeiro 20271-020, Brazil
3
SIG Consultoria e Assessoria Ltda, Rio de Janeiro 22745-004, Brazil
4
Telos Soluções Ambientais, São Paulo 13049-322, Brazil
5
Logistics, Operational Planning and Control, Measurement and Product Inventory Management, PETROBRAS S.A., Rio de Janeiro 20231-030, Brazil
*
Author to whom correspondence should be addressed.
Submission received: 20 May 2022 / Revised: 28 May 2022 / Accepted: 31 May 2022 / Published: 6 June 2022
(This article belongs to the Special Issue Frontiers in Analytical Methods for Water Analysis)

Abstract

:
This study aimed to introduce non-parametric tests and guard bands to assess the compliance of some river water properties with Brazilian environmental regulations. Due to the heterogeneity of the measurands pH, Biochemical Oxygen Demand (BOD), manganese molar concentration, and Escherichia coli, which could be wrongly treated as outliers, as well as the non-Gaussian data, robust methods were used to calculate the measurement uncertainty. Next, based on guard bands, the compliance assessment was evaluated using this previous uncertainty information. For these four measurands, partial overlaps between their uncertainties and the specification limit could generate doubts about compliance. The non-parametric approach for calculating the uncertainty connected to the guard bands concept classified pH and BOD as “conform”, with a risk to the consumer of up to 4.0% and 4.9%, respectively; in contrast, manganese molar concentration and Escherichia coli were “not conform”, with a risk to the consumer of up to 25% and 7.4%, respectively. The methodology proposed was satisfactory because it considered the natural heterogeneity of data with non-Gaussian behavior instead of wrongly excluding outliers. In an unprecedented way, two connected statistical approaches shed light on the measurement uncertainty in compliance assessment of water analysis.

Graphical Abstract

1. Introduction

Due to the increase of the world population and the unrestrained occupation of industry and housing on riverside lands, a severe problem has arisen in the environment: river pollution [1]. The three significant sources of pollution can be considered as derived from industry (chemical waste products that are discharged into rivers), agriculture (fertilizers and pesticides), and domestic activity (people without environmental awareness that throw rubbish into rivers or from cities without proper sanitary wastewater treatment) [2].
Monitoring and evaluating the quality of surface and groundwater are key factors for the proper management of water resources [3]. The characterization and analysis of trends in watersheds are essential for various management activities, such as planning, granting, charging, and framing of watercourses [4].
The concern with river water analysis is worldwide and can be highlighted in the most varied aspects, such as evaluation of fecal coliform (FC) concentration in river water in Kyrgyzstan, Central Asia [5]; development of a new ensemble machine-learning model for predicting monthly water quality at the Lam Tsuen River in Hong Kong [6]; application of soft computing to predict water quality in wetland [7]; exploitation of the water quality of the Surma River by applying hydrochemical and multivariate statistical methods and also with the help of Water Quality Index (WQI) analysis [8]; spatial and seasonal evaluation of the variations of surface water quality in the reservoir using multivariate statistical techniques along with the Trophic State Index and Trophic State Index deviation [9]; assessment of water quality in the polluted stretch using a cluster analysis during pre- and COVID-19 lockdown of the Tawi River Basin, Jammu, North India [10]; performance of an ecotoxicological evaluation of treated wastewater and river water in the catchment scale using a battery of biotests in Poland [11]; optimal operation of an Iranian surface water resources system in terms of quantity and quality simultaneously by using a Non-Dominated Sorting Genetic Algorithm-II (NSGA-II) algorithm [12]; evaluation and classification of the surface water quality in Dong Thap province, Vietnam, using set-pair analysis and national WQI methods [13]; use of radical terraces for erosion control and water quality improvement in Rwanda [14]; detection of phthalates in drinking and river water samples [15]; evaluation of the suitability of an Arabian surface water [16].
Specifically, concerning Brazilian rivers, some recent studies show the importance of this monitoring and evaluation, such as the use of complex network modeling to evaluate the influence of the parameters on the water quality of rivers in the interior of São Paulo State [17]; the variation of water quality in a hydrographic basin in Paraná [18]; the monitoring and modeling for management of the Piabanha River rehabilitation [19]; the development and application of a new method applied to river samples from a populated tropical urban area [20]; the evaluation of the influence of the water quality seasonality in a neotropical river in Bahia [21]; the evaluation of the influence of the land on the variability of water quality in a State of Sergipe river basin [22]; the development and optimization of an alternative analytical method using an ultrasonic bath to extract metals from Doce River sediment samples in the State of Minas Gerais [23].
The water analysis assessment in Brazilian rivers is based on the maximum allowable limits established by the Brazilian National Environmental Council (CONAMA) Resolution 357/2005 [24]. However, not always; it is a trivial task since the temporal heterogeneity of the samples must be treated instead of considering them as outliers. This issue can be potentialized when there is a partial overlap between the most probable value accompanied with its uncertainty and the specification limit, which could generate doubts about compliance with the specification. To solve this problem, this study aimed to innovatively introduce non-parametric tests connected to guard bands to assess the conformity of river water analysis concerning Brazilian environmental regulations.

2. Materials and Methods

This study shows the importance of monitoring the analyzed biochemical parameters in water analyses and describes how to connect robust methods and use uncertainty information in compliance assessment. This section describes an overview of the importance of quality control and how the analysis of each biochemical parameter is described.

2.1. Biochemical Parameters

This section describes an overview of the importance of quality control and how the analysis of each biochemical parameter is described.

2.1.1. pH

Measurement by pH is significant for evaluating the environmental conditions of water collections. There is a range of pH values suitable to lifeforms, but it is usually considered a good water quality indicator if pH lies between 5 and 9 [25].
The samples were collected and analyzed in situ within two hours, using proper pH meter equipment, i.e., a meter that was able to perform the measurements in accordance with the reference [26]. The calibration was based on the supplier orientation by means of an indicating (glass) electrode and a reference electrode using National Institute of Standards and Technology (NIST) buffers.

2.1.2. Biochemical Oxygen Demand (BOD)

Biochemical oxygen demand (BOD) testing determines the relative oxygen requirements of wastewaters, effluents, and polluted waters. This parameter is related to water quality and is widely used to evaluate the environmental health of natural water collections, such as rivers and lakes. It is also important in assessing the performance of wastewater treatment facilities. The standard method recommends filling an airtight bottle of the particularized size with the sample until it overflows and incubating it at the specified temperature (in the dark at 20 °C) and time (5 days). BOD was calculated based on the difference between the initial and final measurements of dissolved oxygen [27]. This Standard Method recommended using the dilution technique for higher concentrations in more polluted water or industrial wastewater.

2.1.3. Manganese molar concentration

Manganese is a common element on Earth’s crust, occurring in the form of mineral. It is the twelfth most abundant element and the fourth most abundant metal for commercial use. The average manganese content in the Earth’s crust is around 0.1 %, and in groundwaters, it is less than 0.1 mg L−1. Manganese generally is present in iron minerals due to its chemical similarity and is found in oceans, fresh waters, and soils. In aqueous matrices, the common aqueous species are Mn2+ and Mn4+. Raised manganese levels can cause stains in domestic utensils, but on the other hand, this metal is very useful for plant and animal survival. The United Nations Food and Agriculture Organization considers that in irrigation waters, its concentration cannot exceed 0.2 mg L−1.
The analytical method used is Inductively Coupled Plasma-Mass Spectrometry (ICP-MS). This method is based on the measurement of ions generated by a radiofrequency inductively coupled plasma. After being nebulized inside the equipment, chemical species are transported as an aerosol by an argon flow that reached the plasma torch. After this, the ions produced are appropriately introduced into a mass spectrometer [28].

2.1.4. Escherichia coli

Detection and evaluation of coliforms are usual indicators of a water’s quality and suitability for domestic, industrial, or other uses. Escherichia coli is used as an indicator of environmental contamination in freshwater systems since this organism is rife in human and animal excrements and not generally found in other ecosystems. Escherichia coli could be easily detected by its ability to ferment glucose, and it may be easier isolated than other known gastrointestinal pathogens.
The samples collected for this work were submitted to enzyme–substrate tests that use hydrolyzable chromogenic and fluorogenic substrates for simultaneously detecting enzymes produced by total coliforms and Escherichia coli. In this method, total coliform bacteria produce the enzyme β-d-galactosidase. Once this enzyme is released in a medium containing chromogenic substrate, this substrate is cleaved to release chromogen. Conversely, most Escherichia coli strains produce the other similar enzyme, the β-glucuronidase, which cleaves a fluorogenic substrate in the medium to release fluorogen. Therefore, by observing whether chromogen or fluorogen is released, it is possible to confirm the presence of Escherichia coli. The release of chromogen indicates that coliform bacteria are present, and the release of fluorogen indicates that Escherichia coli are present [29].
In this study, the Escherichia coli density was reported as a most probable number (MPN), based upon a serial dilution of the test sample.

2.2. Robust Statistical Methods

Robust or non-parametric statistical methods can be used to describe position and precision measures of results that deviate from normal behavior, as they do not always require the treatment of outliers. In general, robust methods should be used in preference to methods that exclude results wrongly considered as outliers [30].
The median (Med), the scaled median absolute deviation (MADe), and the normalized interquartile range (nIQR) are allowed as simple estimators for data sets with normal behavior and that depart from normality. In these cases, Med, MADe, and nIQR have a greater variance than the mean and standard deviation when applied to data with approximately normal distribution; that is, a more conservative approach [31]. However, these estimators can be used with robustness in both situations. Furthermore, the literature emphasizes that different conclusions can be reached by choosing between MADe, Equation (1) and nIQR, Equation (2) [32].
MADe ( x ) = s * MADe = 1.483   Med | x i Med ( x ) |
nIQR ( x ) = s * nIQR = 0.7413 × ( q 3 q 1 )
where q1 and q3 are the first and third quartiles, respectively.
The value to be evaluated, ( x * ) , came from a robust statistic calculated using median and robust standard deviation, s*, Equations (1) and (2). For p observations, the robust standard uncertainty of this value, evaluation of type A, u ( x * ) , could be estimated by Equation (3):
u ( x * ) = 1.25 × s * p
Based on ISO 13528 [31], when using Equation (3), the median value and the robust standard deviation, Equations (1) and (2) were determined from several samples; the uncertainty of median value could be assumed to include the effects of uncertainty concerning heterogeneity. The factor 1.25 (Equation (3)) was based on the variability of the median as a trial to estimate the mean, and it is recommended when testing results have a non-Gaussian behavior. Finally, the expanded uncertainty, U ( x * ) , was calculated by multiplying u ( x * ) by the coverage factor (k) for a confidence level of 95%.

2.3. Use of Uncertainty Information in Compliance Assessment

Conformity assessment can be understood as any process that evaluates if standards and specified requirements are met. In conformity assessment, a measurement result accompanied with its uncertainty can be useful to decide whether an item of interest meets a requirement. This requirement can be evaluated based on one or two tolerance limits [33].
Once the acceptance range of permissible measured values of a measurand is defined, the risks of false acceptance/rejection decisions related to measurement uncertainty can reach a balance to minimize the costs associated with these improper decisions.
Typical scenarios arise when measurement results accompanied by their measurement uncertainties are used to assess compliance with an upper specification limit. When the measured value plus its uncertainty is above or below the control limit is clear. However, in other situations, false decisions can be made due to the partial overlap of the uncertainty bands about the specification limits [34].
A good and current alternative for conformity assessment is the concept of decision rules based on guard bands that are based on measurement uncertainty. Based on this premise, a rejection zone can be outlined as specification limit L plus a value g (called the guard band). This amount was selected in such a way that for a measurement result equal to or greater than L + g, the false rejection probability is less than or equal to α; that is, a low probability that the allowable limit has not been exceeded. Figure 1 illustrates the consumer’s risk, which can be defined as the probability of a false acceptance (i.e., accepting an item or lot outside the specification limits).
An excellent option for applying this approach is the Monte Carlo method (MCM). In a simple way, the basis of Monte Carlo simulations is to randomly select a number with a known probability distribution function and repeat this routine n times [35], using the median and the robust standard uncertainty. Subsequently, histograms can represent graphically this probability density function, and the guard bands help us to define if the item of interest meets a specified requirement. Finally, Bayes’ theorem can be used to estimate the risk of the consumer [36].
The use of uncertainty information in compliance assessment has been published in original research and studies in areas covering the latest in basic and applied research in analytical sciences related to the risk of false conformity assessment applied to automotive fuel analysis employing a multiparameter approach [37]; connected to a data reconciliation approach solving the discrepancy between producers’ and consumers’ measurements [38]; to study the probability of false conforming in customs control of denatured alcohols [39]; to shed light on the importance of analytical chemistry concerning industrial practices [40]; applied to drug and medicine analyses [41]; use of measurement uncertainty in the conformity/non-conformity assessment of pharmaceutical products [42]. Finally, it is worth highlighting an innovative study that discussed the environmental behavior concerning measurement uncertainty of concentrations of pollutants using the guard bands approach [43].

3. Experimental

To characterize the compliance assessment of the water analysis, samples were collected at the same locations, varying over time, at coordinates S 22°18′23.4″, W 41°49′21.5″, in 2020 from January to December, on the right margin of Macaé River, Brazil, Figure 2. The Macaé River basin is about 1710 km2, covering six municipalities in Rio de Janeiro state [44]. In the 1940s, the analyzed section of the Macaé River was rectified to facilitate the drainage of swamps in the region. This process changed the behavior of the river, its depth, and the speed of its waters.
The sampling and analytical strategy is based on the environmental monitoring of the Macaé River for use in a thermal power plant electrical system. The number of analyses was 52 for pH (every seven days), 94 for BDO (every four days), 72 (every five days) for Manganese molar concentration, and 60 for Escherichia coli (every six days). All the samples were collected at the same point by a system of stainless-steel bucket, separated into proper sample bottles, and sent to the laboratory in a thermal shipping box.
Every sample was labeled and sent to a laboratory except pH measurements, which were analyzed in situ using a Sanxin SX836 portable pH meter, in accordance with the preservation procedures for each parameter. The samples from the Macaé River Margin were carried to the laboratory in ice-cooled boxes to avoid degradation. After standard preservation procedures, the samples were sent to an accredited laboratory for the BOD, Mn, and Escherichia coli measurements following the selected measurement procedure for each measurand by ISO/IEC 17025 [45], located in Rio de Janeiro City, 190 km from the sample point. This laboratory has procedures to monitor the validity of results, such as the use of certified reference materials, use of control charts, participation in intralaboratory and interlaboratory comparisons, etc.

4. Results

The measurement uncertainty is considered as the square root of the quadratic sum of the analytical uncertainty and sampling uncertainty. The approach used in this study considered both the analytical variability and the temporal variability, without stratifying each uncertainty source. The raw data (before any statistical treatment) are available in Table 1.

4.1. Preliminary Statistical Evaluation

Data from Table 1 were treated, and they do not follow a normal distribution based on the Kolmogorov–Smirnov test; in addition, several outliers are met [46].

4.2. Statistical Treatment

Here, Equations (1)–(3) were applied to all data from Table 1 without excluding any outliers since this study has considered that the variability of the results is due to the temporal heterogeneity of the sampling target. The results of the statistical treatment are presented in Table 2. In a traditional approach, i.e., always considering that data are normally distributed, and there are no outliers, U ( x ) SD is calculated as   k × SD / p , where SD is the standard deviation. There are partial overlaps between the central measures (medians) plus expanded uncertainties and their respective specification limits. It can generate doubts concerning the compliance assessment.

4.3. Compliance Assessment

The uncertainty is used to evaluate the conformity/non-conformity with the specification. In this study, the value of g is 1.64 u ( x * ) for a significance level of 5%, where u ( x * ) is calculated from Equation (3).
Next, histograms with the central value (median of actual data distribution without excluding outliers), its respective uncertainty, i.e., guard band (1.64 u ( x * ) ) from Equation (3), and specification limit(s) are presented. In this study, Monte Carlo simulations were carried out with 100000 pseudorandom values for each water analysis parameter. Here, the consumer’s risk is evaluated, i.e., the risk of a false acceptance of the analyzed item.

4.3.1. pH

Figure 3a,b show the pH conformity assessment based on MADe and nIQR, respectively.
The lower and the upper acceptance limits for pH were 6.10 and 8.90 (MADe approach) and 6.19 and 8.81 (nIQR approach), respectively.

4.3.2. Biochemical Oxygen Demand (BOD)

Figure 4a,b show the BOD conformity assessment based on MADe and nIQR, respectively.
The upper acceptance limit for BOD was 4.42 mg L−1 (MADe approach) and 4.50 mg L−1 (nIQR approach), respectively.

4.3.3. Manganese Molar Concentration

Figure 5a,b show the conformity assessment of the manganese molar concentration based on MADe and nIQR, respectively.
The upper acceptance limit for Mn was 0.094 mol L−1 (MADe approach) and 0.090 mol L−1 (nIQR approach), respectively.

4.3.4. Escherichia coli

Figure 6a,b show the Escherichia coli conformity assessment based on MADe and nIQR, respectively.
The upper acceptance limit for Escherichia coli was 970 MPN per 100 mL (MADe approach) and 971 MPN per 100 mL (nIQR approach), respectively.

4.3.5. Summarizing Different Approaches

To prepare the discussion in the next section, Table 3 summarizes the comparison of the results for compliance assessment using the different methods (including the traditional approach, mean plus uncertainty arising standard deviation) and the percent consumer risk for each case.

5. Discussion

Some results showed partial overlaps between the central measures (medians) plus expanded uncertainties and their respective specification limits. It could generate doubts concerning the compliance assessment. Therefore, in this section, we discuss how our study could contribute to the solution of this problem. The robust methods were based on Equations (1)–(3), and the guard bands were calculated applying 1.64 u ( x * ) .
Concerning the first property studied, pH, the conformity test has demonstrated, with a probability greater than 95% (significance level of 5%), that the pH, 6.20 (median), conforms to the requirement based on both MADe and nIQR approaches, even though the expanded uncertainty arising from nIQR approach is almost twice the MADe approach, in what was reflected in the risk for the consumer, although all of them were less than 5%. When the standard deviation approach was considered, no partial overlap was noticed, and the compliance assessment was clearly defined. For this biochemical parameter, all approaches presented agreed among themselves.
Regarding the second parameter, BOD, the conformity test has demonstrated, with a probability greater than 95% (significance level of 5%), that the BOD, 4.50 mg L−1 (median), conforms to the requirement based on both MADe and nIQR approaches, even though the expanded uncertainty arising from the nIQR approach is greater than the MADe approach, in what was reflected in the risk for the consumer, although both were less than 5%. When the standard deviation approach was considered, no partial overlap was noticed, and the non-compliance assessment was clearly defined. For this biochemical parameter, the nIQR approach seemed to be more appropriated because it was more conservative than the MADe approach; furthermore, the statistical methods were powerful enough to handle the absence of normality and the presence of outliers. It probably occurred because approximately 50% of the results were the same, so MADe tended to zero, not being a good estimator. Although 33% of BOD data are above the specification limit, they are not “lost”. These sample data estimate the population over time and apply robust methods to calculate the position measure (median), accompanied by its dispersion measures (MADe and nIQR) that are part of the entire data set. Here, we understand that the decision-making was based on all results during 2020 that generated a median value over time considering the heterogeneity of the data, and not just on a percentage of BOD data (33%).
For the next parameter, Mn, the conformity test has demonstrated, with a probability greater than 95% (significance level of 5%), that the manganese molar concentration, 0.096 mol L−1 (median), did not comply with the requirement based on both MADe and nIQR approaches. The expanded uncertainty arising from the nIQR approach is similar to the MADe approach; however, the risk for the consumer for nIQR approach is greater than the MADe approach. When the standard deviation approach was considered, no partial overlap was noticed, and the non-compliance assessment was clearly defined.
For the last property, Escherichia coli, the conformity test has demonstrated, with a probability greater than 95% (significance level of 5%), that Escherichia coli, 1026 MPN per 100 mL (median), does not conform to the requirement based on both MADe and nIQR approaches. Both the expanded uncertainty and the consumer’s risk arising from nIQR approach are similar to the MADe approach. On the other hand, no partial overlap was noticed when the standard deviation approach was considered, and the non-compliance assessment was clearly defined.
In the cases studied, the robust standard deviations calculated by MADe and nIQR approaches, which were used in the guard bands calculations, reached the same conclusions concerning compliance/non-compliance with the Brazilian specification, probably because the robustness of these estimators against outliers was not exceeded. Based on the robustness of statistical tests and knowledge about water analysis, consistent conclusions can be reached and supported by the expertise of the professional involved in this task. This study did not aim to make a decision based on a specific analysis or a fraction of samples to assess the compliance of river water analysis concerning Brazilian environmental regulations, but rather all the analyses carried out during 2020 considered the natural non-Gaussian behavior of the data and the median below/above of the specification limit, but with this limit within the uncertainty range.

6. Conclusions

This study aimed to use robust methods connected to the guard bands to assess the conformity when there is a partial overlap of the uncertainty bands about the specification limits of some biochemical parameters of the river water analysis. The methodology was considered suitable and promising since it respected the particularity of data sets.
The temporal heterogeneity of the samples was treated as a factor concerning the behavior of the samples and not as outliers. Moreover, robust methods connected to the concept of guard bands were beneficial for evaluating the risk of false conformity when there is a partial overlap between the most probable value accompanied by its uncertainty and the specification limit.
For future works, we propose applying this methodology to other water quality parameters and areas whose data are expected to have similar behavior, e.g., the classification of waste from steel companies and air quality. In addition, we recommend considering the spatial heterogeneity of the sampling target, as recently proposed by [47], and segregating the contribution from the sampling step and the analytical one.

Author Contributions

Conceptualization, E.C.d.O.; methodology, E.C.d.O.; software, E.C.d.O.; validation, L.P.B., V.F.S., M.B. and E.C.d.O.; formal analysis, L.P.B., V.F.S., M.B. and E.C.d.O.; investigation, L.P.B., V.F.S., M.B. and E.C.d.O.; resources, L.P.B., V.F.S., M.B. and E.C.d.O.; data curation, L.P.B., V.F.S., M.B. and E.C.d.O.; writing—original draft preparation, L.P.B., V.F.S., M.B. and E.C.d.O.; writing—review and editing, L.P.B., V.F.S., M.B. and E.C.d.O.; visualization, L.P.B., V.F.S., M.B. and E.C.d.O.; supervision, E.C.d.O.; project administration, E.C.d.O.; funding acquisition, L.P.B. and E.C.d.O. All authors have read and agreed to the published version of the manuscript.

Funding

The authors thank for the financial support provided by the Brazilian funding agencies CNPq, CAPES, FINEP and FAPERJ. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data was obtained from the Telos Soluções Ambientais and are available from the authors with the permission of them.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mateo-Sagasta, J.; Zadeh, S.M.; Turral, H.; Burke, J. Water Pollution from Agriculture: A Global Review. Executive Summary; FAO: Rome, Italy; International Water Management Institute (IWMI): Colombo, Sri Lanka, 2017; p. 35. [Google Scholar]
  2. Akhtar, N.; Syakir Ishak, M.I.; Bhawani, S.A.; Umar, K. Various Natural and Anthropogenic Factors Responsible for Water Quality Degradation: A Review. Water 2021, 13, 2660. [Google Scholar] [CrossRef]
  3. Klamt, R.A.; Lobo, E.A.; da Costa, A.B.; Delevati, D.M. Evaluation of water resource preservation areas in the Hydrographical Basin of Andreas Stream, RS, Brazil, using environmental monitoring programs. Rev. Ambiente Água 2019, 14, e2307. [Google Scholar] [CrossRef]
  4. Thornton, J.A.; Harding, W.R.; Slawski, T.M.; Lin, H. Monitoring and Evaluation: The Foundation for Lake and Reservoir Management. Earth 2022, 3, 115–124. [Google Scholar] [CrossRef]
  5. Li, Y.; Ma, L.; Li, Y.; Abdyzhapar Uulu, S.; Abuduwaili, J. Exploration of the driving factors and distribution of fecal coliform in rivers under a traditional agro-pastoral economy in Kyrgyzstan, Central Asia. Chemosphere 2022, 286, 131700. [Google Scholar] [CrossRef]
  6. Asadollah, S.B.H.S.; Sharafati, A.; Mott, D.; Yaseen, Z.M. River water quality index prediction and uncertainty analysis: A comparative study of machine learning models. J. Environ. Chem. Eng. 2021, 9, 104599. [Google Scholar] [CrossRef]
  7. Pham, Q.B.; Mohammadpour, R.; Linh, N.T.T.; Mohajane, M.; Pourjasem, A.; Sammen, S.S.; Anh, D.T.; Nam, V.T. Application of soft computing to predict water quality in wetland. Environ. Sci. Pollut. Res. 2021, 28, 185–200. [Google Scholar] [CrossRef]
  8. Howladar, M.F.; Chakma, E.; Koley, N.J.; Islam, S.; Numanbakth, M.A.A.; Ahmed, Z.; Chowdhury, T.R.; Akter, S. The water quality and pollution sources assessment of Surma river, Bangladesh using, hydrochemical, multivariate statistical and water quality index methods. Groundw. Sustain. Dev. 2021, 12, 100523. [Google Scholar] [CrossRef]
  9. Mamun, M.; Kim, J.Y.; An, K.-G. Multivariate Statistical Analysis of Water Quality and Trophic State in an Artificial Dam Reservoir. Water 2021, 13, 186. [Google Scholar] [CrossRef]
  10. Kour, G.; Kothari, R.; Dhar, S.; Pathania, D.; Tyagi, V.V. Impact assessment on water quality in the polluted stretch using a cluster analysis during pre- and COVID-19 lockdown of Tawi river basin, Jammu, North India: An environment resiliency. Energ. Ecol. Environ. 2021; in press. [Google Scholar] [CrossRef]
  11. Szklarek, S.; Kiedrzynska, E.; Kiedrzynski, M.; Mankiewicz-Boczek, J.; Mitsch, W.J.; Zalewski, M. Comparing ecotoxicological and physicochemical indicators of municipal wastewater effluent and river water quality in a Baltic Sea catchment in Poland. Ecol. Indic. 2021, 126, 107611. [Google Scholar] [CrossRef]
  12. Babamiri, O.; Marofi, S. A multi-objective simulation–optimization approach for water resource planning of reservoir–river systems based on a coupled quantity–quality model. Environ. Earth Sci. 2021, 80, 389. [Google Scholar] [CrossRef]
  13. Giao, N.T.; Nhien, H.T.H.; Anh, P.K.; Ni, D.V. Use of Complex Network Modelling to Assess the Influence of the Parameters on Water Quality of Rivers. Environ. Monit. Assess 2021, 193, 319. [Google Scholar] [CrossRef] [PubMed]
  14. Uwacu, R.A.; Habanabakize, E.; Adamowski, J.; Schwinghamer, T.D. Using radical terraces for erosion control and water quality improvement in Rwanda: A case study in Sebeya catchment. Environ. Dev. 2021, 39, 100649. [Google Scholar] [CrossRef]
  15. Yang, D.; Yang, Y.; Li, Y.; Yan, S.; Chen, Y.; Wang, J.; Xiao, J.; Sun, C. Dispersive Liquid—Liquid Microextraction Based on Solidification of Floating Organic Drop Combined with High Performance Liquid Chromatography for Analysis of 15 Phthalates in Water. J. AOAC Int. 2019, 102, 942–951. [Google Scholar] [CrossRef] [PubMed]
  16. Al-Othman, A.A. Evaluation of the suitability of surface water from Riyadh Mainstream Saudi Arabia for a variety of uses. Arab. J. Chem. 2019, 12, 2104–2110. [Google Scholar] [CrossRef] [Green Version]
  17. Gayer, F.A.M.; de Angelis, D.F.; de Angelis, A.F.; Poletti, E.C.C. Use of Complex Network Modelling to Assess the Influence of the Parameters on Water Quality of Rivers. Water Air Soil Pollut. 2021, 232, 324. [Google Scholar] [CrossRef]
  18. Godoy, R.F.B.; Crisiogiovanni, E.L.; Trevisan, E.; Radomski, F.A.D. Spatial and temporal variation of water quality in a watershed in center-west Paraná, Brazil. Water Supply 2021, 21, 1734–1788. [Google Scholar] [CrossRef]
  19. Costa, D.A.; Silva Junior, L.C.S.; Azevedo, J.P.S.; Santos, M.A.; Assumpção, R.S.F.V. From Monitoring and Modeling to Management: How to Improve Water Quality in Brazilian Rivers? A Case Study: Piabanha River Watershed. Water 2021, 13, 176. [Google Scholar] [CrossRef]
  20. Stoll, G.C.; Carreira, R.S.; Massone, C.G. Polychlorinated biphenyls (PCBs) in water: Method development and application to river samples from a populated tropical urban area. Anal. Bioanal. Chem. 2020, 412, 2477–2486. [Google Scholar] [CrossRef]
  21. Aboim, I.L.; Gomes, D.F.; Junior, P.O.M. Phytoplankton response to water quality seasonality in a Brazilian neotropical river. Environ. Monit. Assess 2020, 192, 70. [Google Scholar] [CrossRef]
  22. Cruz, M.A.S.; Gonçalves, A.A.; Aragão, R.; Amorim, J.R.A.; Mota, P.V.M.; Srinivasan, V.S.; Garcia, C.A.B.; Figueiredo, E.E. Spatial and seasonal variability of the water quality characteristics of a river in Northeast Brazil. Environ. Earth Sci. 2019, 78, 68. [Google Scholar] [CrossRef]
  23. Mimura, A.M.S.; Oliveira, M.A.L.; Ciminelli, V.S.T.; Silva, J.C.J. Optimization of Ultrasound-Assisted Extraction of Cr, Cu, Zn, Cd, and Pb from Sediment, Followed by FAAS and GFAAS Analysis. J. AOAC Int. 2016, 99, 252–259. [Google Scholar] [CrossRef] [PubMed]
  24. Brazilian National Environment Council (CONAMA). Resolution no. 357. Provides the Classification of Water Bodies and Environmental Guidelines for Their Framework, as Well as Stablishes the Conditions and Standards for Effluents Discharge; Brazilian National Environment Council (CONAMA): Brasília, Brazil, 2005; pp. 1–27. [Google Scholar]
  25. APHA; AWWA; WPCF. 4500-h+ ph. In Standard Methods For the Examination of Water and Wastewater, 23rd ed.; American Public Health Association (APHA): Washington, DC, USA, 2018. [Google Scholar] [CrossRef]
  26. Rajan, A.; Al Raisi, A.; Thankamony, R.; Perumal, P.; Al Hosani, S. Eutrophication sources, impacts and management: A case study from Abu Dhabi. Aquat. Ecosyst. Health Manag. 2020, 23, 175–186. [Google Scholar] [CrossRef]
  27. APHA; AWWA; WPCF. Standard Methods For the Examination of Water and Wastewater, 23rd ed.; American Public Health Association (APHA): Washington, DC, USA, 2019. [Google Scholar] [CrossRef]
  28. APHA; AWWA; WPCF. 3500-mn manganese. In Standard Methods For the Examination of Water and Wastewater, 23rd ed.; American Public Health Association (APHA): Washington, DC, USA, 2018. [Google Scholar] [CrossRef]
  29. APHA; AWWA; WPCF. Standard Methods For the Examination of Water and Wastewater; American Public Health Association (APHA): Washington, DC, USA, 2017. [Google Scholar] [CrossRef]
  30. Bueno, A.F.; Ozaki, D.A.F.; Barbosa, E.; dos Santos, E.E.; Gomes, M.M.; Iida, P.H.; dos Santos, S.C.A.; de Sá, F.G.; de Oliveira, E.C. Validation of the performance of process stream analyzer systems with nonparametric behavior. Accredit. Qual. Assur. 2014, 19, 185–193. [Google Scholar] [CrossRef]
  31. ISO 13528:2015; Statistical Methods for Use in Proficiency Testing by Interlaboratory Comparisons. International Organization for Standardization: Geneva, Switzerland, 2015.
  32. Holleben, C.; Oliveira, E.C. Variations in the classification of laboratory performance due to different protocols for use in proficiency testing programs. Measurement 2020, 152, 107354. [Google Scholar] [CrossRef]
  33. JCGM 106: 2012; Evaluation of Measurement Data—The Role of Measurement Uncertainty in Conformity Assessment. Joint Committee for Guides in Metrology: Sèvres, France, 2012. Available online: www.bipm.org/en/publications/guides/gum.html (accessed on 26 May 2022).
  34. Ellison, S.L.R.; Williams, A. (Eds.) EURACHEM/CITAC Guide: Use of Uncertainty Information in Compliance Assessment, 2nd ed.; Construction Innovation and Technology Application Centre: Hong Kong, China, 2021; Available online: https://www.eurachem.org/index.php/publications/guides (accessed on 14 September 2021).
  35. De Oliveira, E.C.; Biazon, C.L.; Moreira, R.M.; Filho, P.L.S. Uncertainty evaluation in the determination of oil and grease content in produced water by colorimetric method using Monte Carlo Simulation. Pet. Sci. Technol. 2019, 37, 436–442. [Google Scholar] [CrossRef]
  36. Toczek, W.; Smulko, J. Risk Analysis by a Probabilistic Model of the Measurement Process. Sensors 2021, 21, 2053. [Google Scholar] [CrossRef]
  37. De Oliveira, E.C.; Lourenço, F.R. Risk of false conformity assessment applied to automotive fuel analysis: A multiparameter approach. Chemosphere 2021, 263, 128265. [Google Scholar] [CrossRef]
  38. De Oliveira, E.C.; Lourenço, F.R. Data reconciliation applied to the conformity assessment of fuel products. Fuel 2021, 300, 120936. [Google Scholar] [CrossRef]
  39. Kuselman, I.; Pennecchi, F.; da Silva, R.J.N.B.; Hibbert, D.B. Conformity assessment of multicomponent materials or objects: Risk of false decisions due to measurement uncertainty—A case study of denatured alcohols. Talanta 2017, 164, 189–195. [Google Scholar] [CrossRef]
  40. De Oliveira, E.C. Use of Measurement Uncertainty in Compliance Assessment with Regulatory Limits. Braz. J. Anal. Chem. 2020, 7, 1–2. [Google Scholar] [CrossRef]
  41. Separovic, L.; Simabukuro, R.S.; Couto, A.R.; Bertanha, M.L.G.; Dias, F.R.S.; Sano, A.Y.; Caffaro, A.M.; Lourenço, F.R. Measurement Uncertainty and Conformity Assessment Applied to Drug and Medicine Analyses—A Review. Crit. Rev. Anal. Chem. 2021; in press. [Google Scholar] [CrossRef] [PubMed]
  42. Simabukuro, R.; Jeong, N.A.; Lourenço, F.R. Application of Measurement Uncertainty on Conformity Assessment in Pharmaceutical Drug Products. J. AOAC Int. 2021, 104, 585–591. [Google Scholar] [CrossRef] [PubMed]
  43. Pennecchi, F.R.; Kuselman, I.; da Silva, R.J.N.B.; Hibbert, D.B. Risk of a false decision on conformity of an environmental compartment due to measurement uncertainty of concentrations of two or more pollutants. Chemosphere 2018, 202, 165–176. [Google Scholar] [CrossRef] [PubMed]
  44. Junior, J.L.; Tavares, L.P.S.; Kalas, F.A.; Rodrigues, P.P.G.W.; Wasserman, J.C.A. Reservoir implantation for flood dampening in the Macaé River basin using the Mohid Land model. Ciência Nat. 2019, 41, e29. [Google Scholar] [CrossRef]
  45. ISO/IEC 17025:2017; General Requirements for the Competence of Testing and Calibration Laboratories. International Organization for Standardization: Geneva, Switzerland, 2017.
  46. Oliveira, E.C.; Orlando, A.F.; Ferreira, A.L.S.; Chaves, C.E.O. Comparison of different approaches for detection and treatment of outliers in meter proving factors determination. Flow Meas Instrum. 2016, 48, 29–35. [Google Scholar] [CrossRef]
  47. Borges, C.; Palma, C.; Dadamos, T.; Bettencourt da Silva, R.J.N. Evaluation of seawater composition in a vast area from the Monte Carlo simulation of georeferenced information in a Bayesian framework. Chemosphere 2021, 263, 128036. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Acceptance and rejection zones for simultaneous upper and lower limits based on guard bands.
Figure 1. Acceptance and rejection zones for simultaneous upper and lower limits based on guard bands.
Molecules 27 03628 g001
Figure 2. Macaé River.
Figure 2. Macaé River.
Molecules 27 03628 g002
Figure 3. pH conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper ac-ceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Figure 3. pH conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper ac-ceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Molecules 27 03628 g003aMolecules 27 03628 g003b
Figure 4. BOD conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Figure 4. BOD conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Molecules 27 03628 g004aMolecules 27 03628 g004b
Figure 5. Mn molar concentration conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Figure 5. Mn molar concentration conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Molecules 27 03628 g005aMolecules 27 03628 g005b
Figure 6. Escherichia coli conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Figure 6. Escherichia coli conformity assessment based on (a) MADe approach and (b) nIQR approach. p(AL)—probability density at the lower acceptance limit; p(AU)—probability density at the upper acceptance limit; AL—lower acceptance limit; AU—upper acceptance limit; TL—lower tolerance limit; TU—upper tolerance limit.
Molecules 27 03628 g006aMolecules 27 03628 g006b
Table 1. Raw data of some biochemical properties of the Macaé River collected in 2020.
Table 1. Raw data of some biochemical properties of the Macaé River collected in 2020.
pH ResultsBOD Results (mg L−1)
7.205.105.236.306.235.904.4018.004.504.504.503.97
5.996.016.006.105.806.104.404.504.50101.004.503.97
5.865.865.915.935.966.0091.914.504.504.506.234.18
6.126.196.206.216.246.294.2092.3077.854.504.504.18
6.096.196.026.736.126.864.604.504.5039.184.505.72
6.977.057.087.237.307.534.804.504.504.504.505.72
6.307.506.306.306.576.344.404.504.504.505.915.91
5.908.106.066.116.106.204.504.504.504.504.504.50
6.936.947.667.65 7.804.5083.2522.372.336.23
8.90101.004.504.502.334.50
10.004.5098.974.503.486.92
3.004.504.504.503.486.92
6.9883.258.7691.914.504.50
6.984.509.174.5022.3798.87
8.764.109.1792.304.504.50
39.184.5077.854.50
Manganese results (mg L−1)Escherichia coli results (MPN * per 100 mL)
0.3760.0790.0810.0790.0960.115901113891450010581043
0.0910.2700.0920.0300.0960.115110011311048500911952
0.2800.3200.0920.0700.0960.13711261000104297019001019
0.0810.2900.0990.2100.0960.13710989641023110210341096
0.1390.0830.0730.0770.0960.0819021220111110421079926
0.0950.0830.0730.0900.0980.139500122011281056983950
0.3100.0830.3000.2300.1000.140901800998111411101118
0.0570.0850.0760.0910.1000.14010001000952965960933
0.3300.0950.1260.0910.1000.16810281121103699110941030
0.0120.0910.0780.0920.1000.168100098592211001139994
0.0890.0910.0780.2000.1130.376
0.0790.0810.0990.0960.1130.080
* Most Probable Number.
Table 2. Statistical results.
Table 2. Statistical results.
pHBOD (mg L−1)Mn (mol L−1)Escherichia coli (MPN per 100 mL)
Specification limit ‡6–95 0.1 1000
Median6.204.50.09561026
Mean6.4017.20.1281016
Standard deviation (SD)0.6229.10.081184
IQR0.8703.0950.0553142
s * nIQR 0.6452.29430.0410105
s * MADe 0.3340.14830.0248109
p52947260
U ( x ) SD0.186.10.01948
U ( x * ) nIQR0.230.600.01235
U ( x * ) MADe0.120.040.00736
g = 1.64 × 1.25 × s * nIQR p 0.180.480.01028
g = 1.64 × 1.25 × s * MADe p 0.0950.030.00629
Upper specification limit, ‡ based on Brazilian regulations.
Table 3. Summarizing and comparing different approaches.
Table 3. Summarizing and comparing different approaches.
pHBOD (mg L−1)Mn (mol L−1)Escherichia coli
(MPN per 100 mL)
ResultConsumer’s RiskResultConsumer’s RiskResultConsumer’s RiskResultConsumer’s Risk
Median and
U ( x * ) MADe
6.20   ±   0.12 0.04%4.50 ± 0.040.00%0.096 ± 0.00712.7%1026 ± 367.4%
Median and
U ( x * ) nIQR
6.20   ±   0.23 4.0%4.50 ± 0.604.9%0.096 ± 0.01225.3%1026 ± 356.9%
Mean and
U ( x ) SD
6.40   ±   0.18 0.00%17.2 ± 6.1Without
overlap
0.128 ± 0.019Without
overlap
1016 ± 4825.2%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brandão, L.P.; Silva, V.F.; Bassi, M.; de Oliveira, E.C. Risk Assessment in Monitoring of Water Analysis of a Brazilian River. Molecules 2022, 27, 3628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27113628

AMA Style

Brandão LP, Silva VF, Bassi M, de Oliveira EC. Risk Assessment in Monitoring of Water Analysis of a Brazilian River. Molecules. 2022; 27(11):3628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27113628

Chicago/Turabian Style

Brandão, Luciene Pires, Vanilson Fragoso Silva, Marcelo Bassi, and Elcio Cruz de Oliveira. 2022. "Risk Assessment in Monitoring of Water Analysis of a Brazilian River" Molecules 27, no. 11: 3628. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27113628

Article Metrics

Back to TopTop