Next Article in Journal
Effects of Two-Phase Flow of Water and Air on Shallow Slope Failures Induced by Rainfall: Insights from Slope Stability Assessment at a Regional Scale
Next Article in Special Issue
Responsible Water Reuse Needs an Interdisciplinary Approach to Balance Risks and Benefits
Previous Article in Journal
Indigenous Perspectives on Water Security in Saskatchewan, Canada
Previous Article in Special Issue
Water Security Assessment of Groundwater Quality in an Anthropized Rural Area from the Atlantic Forest Biome in Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach to Use Load Duration Curves to Evaluate Water Quality: A Study in the Doce River Basin, Brazil

1
Department of Agricultural Engineering, Federal University of Vicosa (UFV), Viçosa-MG 36570-900, Brazil
2
Minas Gerais Water Management Institute (IGAM), Belo Horizonte-MG 30000-000, Brazil
3
Department of Water Resources, Vale do Rio Verde University (Unincor), Três Corações-MG 37417-150, Brazil
*
Author to whom correspondence should be addressed.
Water 2020, 12(3), 811; https://doi.org/10.3390/w12030811
Submission received: 19 February 2020 / Revised: 9 March 2020 / Accepted: 10 March 2020 / Published: 14 March 2020
(This article belongs to the Special Issue Water Security and Governance in Catchments)

Abstract

:
Although water availability depends both on qualitative and quantitative aspects, most studies focus only on one of these. Therefore, the goal here is to relate water quality and quantity with the construction of Load Duration Curves (LDC) and to estimate E. coli load patterns in different flow conditions, seasons, and positions of two sub-basins of the Doce watershed (Brazil): Piracicaba and Piranga. A novel methodology is proposed in which the Burr XII distribution is adjusted to the LDC to compare all observed loads to their respective Total Maximum Daily Load (TMDL), allowing the estimation of the relative difference (RD) between these. Higher values of RD were observed for low flows for the Piracicaba basin, more urbanized, where point sources of pollution are the primary concern, reaching up to 99% of needed load reduction. In the Piranga basin, more agricultural, there was a broader RD variation, from 9% to 97% load reduction needed, which is an evidence of point sources of pollution combined with non-point sources. The new methodology can be used to estimate the load reduction of any pollutant and can be used by environmental agencies to identify effective practices to minimize and control pollution in different locations of the basins.

Graphical Abstract

1. Introduction

Water availability depends on the quantitative and qualitative assessments of this resource. Therefore, to ensure a sustainable future, these evaluations must be taken altogether. Nevertheless, in spite of this pre-requisite, most studies only cover the quantification and/or prevision of streamflow [1,2] or the characterization and/or modeling of water quality [3,4,5].
Watershed models consist of tools that can be used to assess water quantity and quality simultaneously. However, these approaches often require a diversity of input data and information to run, which can be viewed as a limitation for their use in the developing countries where data is fragmented, uncertain, and barely available [6,7,8,9]. The limited portability of watershed models can be overcome by simpler quantity-quality methods, which are less dependent on flow and transport equations and their adjustment parameters, such as Load Duration Curves (LDC).
An LDC is a graphical representation of water quality, namely observed and total maximum daily loads. The TMDL is an upper threshold of a predefined water quality criterion. For example, the TMDL for Escherichia coli (E. coli) can be computed considering a concentration of 1000 MPN dL−1 (Most Probable Number) for most Brazilian rivers, according to the Environment National Council (CONAMA), which is multiplied by the Flow Duration Curve (FDC) of any given gauge. The setup of LDCs for observed and total maximum daily loads requires the existence of streamflow gauges close or in the same location of water quality gauges because it depends on the generation of the FDC. Thus, when observed loads are plotted with their respective TMDL curve, it is possible to understand the water quality standard, being in conformity when the observations are plotted below the TMDL curve.
The advantage of using LDCs instead of watershed models relies on their capability to make evident the links between streamflow conditions and the expected water quality. In the sequel, the use of LDCs can help to implement more successful efforts to improve water quality in the catchment, especially in the cases where flow regimes are typically characterized by recurrent pollution [10,11,12]. Although the analysis of concentrations can be used to describe water quality in a river, the analysis of loads allows a direct comparison between concentration and streamflow. Thus, the analysis of loads can improve watershed planning through the provision of a better description of water quality concerns [13]. For instance, the LDCs can easily portray the streamflow classes in disagreement with the regulations, and therefore, shed light on best management practices to improve water quality [14,15]. In that context, pollution events related to shorter permanence (i.e., high flows) are usually associated with diffuse pollution sources, whereas events related to higher permanence (i.e., low flows) are commonly linked to point source pollution [11,13,16].
In Brazil, there have been some attempts to use LDCs, but the frequent non-synchronous streamflow and water quality sampling represented a limitation [17,18,19,20]. For instance, FDCs were developed for rivers in the state of São Paulo, Minas Gerais and Paraná (SP) but considered the monthly average of streamflow in the estimation of the load [18,20]. The monthly average makes the day-to-day flow variation smoother, thereby, hampering the comparison between the loading capacity and the observed pollutant load into the river.
A relevant application of LDCs refers to the quantification of a potential load reduction per hydrological regime. The United States Environmental Protection Agency (USEPA) set up conventional ranges for streamflow regimes, namely 0–10% for the high flows, 10–40% for moist conditions, 60%–90% for mid-range conditions, 60%–90% for dry conditions and 90%–100% for low flows. To obtain the needed load reduction, i.e., for E. coli contamination, the middle flow exceedance percentile (from the LDC) is multiplied by the 90th percentile of sampled observations in each flow regime [21,22]. This approach is inherently a setback because the load variation represented in the observations is not well-captured by this technique and will be lost. The extreme hydrological conditions, i.e., high (0–10%) and low flows (90%–100%), are also not accounted as critical conditions for load reduction. It is worth noting that the exclusion of these marginal fringes can compromise the analysis when the observed data is limited, which is a common situation in developing countries.
These limitations inspired the quantification of a potential load reduction per hydrological regime based on a new approach. This would become the purpose of this study, detailed in the forthcoming sentences. With the new approach, a trend line of observed loads is estimated using the LOESS smoothing technique [23], while the calculated TMDLs are fitted to a non-linear line using the Burr XII distribution [24,25]. Having defined these two tendencies, they can be visually compared to identify impaired pollutant loads in a straightforward manner. They can also be used to quantify load reductions as relative differences between the two curves. While using this new method, the percent of load reduction for each streamflow regime would be a summary of all samples collected in that regime, instead of standing on the 90th percentile concentration or on a single load value.
The new method was tested in two watersheds from the Brazilian state of Minas Gerais, namely the Piracicaba and Piranga River basins. In these catchments, many water bodies are highly concentrated in E. coli, which is evidence of fecal contamination and the potential presence of pathogens. This fact has motivated the selection of E. coli as the pollutant in this load assessment study. The proposed methodology allows the comparison of each observed load to the maximum allowable load for the same exceedance permanence of flow, unlike the existing methodology used by the USEPA, thus providing a more representative evaluation of water quality. Thus, besides the presentation of a new method to quantify a required load reduction, the aim here was to increase our understanding of E. coli load variation and thus optimize resource allocation and selection of appropriate best management practices to reduce the percentage of impairment in the basins of study.

2. Materials and Methods

2.1. Study Area

The study was developed in two neighbor watersheds, the Piranga and Piracicaba basins, within the Doce river basin (MG), illustrated in Figure 1. This basin received worldwide attention when a tailing dam collapsed in 2015, destroying a small city and compromising the water quality from the basin headwaters towards the sea [26].
Figure 1 displays the location of the basins, the gauges, and also the boundaries of cities where the gauges are located. The position of the gauges and the different drainage areas throughout the basin allow the assessment of the water quality patterns. These areas, as well as the names of the cities and their population, are depicted in Table 1.
The basins are mostly occupied by not too populous cities. Thus, given the lack of resources from these cities, most of them lack wastewater treatment plants, which causes most of their generated sewage to be loaded into the rivers without treatment. This situation, associated with pasture and swine activities in the basins, has caused E.coli contamination to be one of their core water quality problems [27].
Regarding the loading capacities, CONAMA establishes that the loadings must respect the natural capacity of the river, which means that, after receiving the load, the water quality parameters of the river must meet its class of use. There is no maximum concentration of E. coli established in the loading legislation. Industries or cities that intend to load pollutants into rivers must apply for licenses in the responsible environmental agency and respect a list of minimum requirements, which must be fulfilled unless an unusual situation is experienced. However, given the lack of monitoring agents, situations in which this legislation is not met are frequent, compromising the water quality in many Brazilian water bodies [28].

2.2. Data Used in the Study

To understand the nexus between E.coli contamination and the hydrological patterns in the basins, the maximum allowed E. coli load was compared to the observed E.coli load in various water quality gauges. The maximum allowed load combines the water quality criteria set up by the CONAMA, which is 1000 MPN dL−1, with the flow rate at the streamflow gauge, obtained with the development of an FDC. By multiplying the FDC with the water quality criteria, it is then possible to obtain the LDC, which relates the maximum allowed load to the percentage that can be observed in time (and flow permanence), the Total Maximum Daily Load (TMDL) [29,30]. Larger allowable loads are associated with a shorter permanence, since these are associated with higher flows when the dilution capacity of the river is more elevated. On the other hand, smaller allowable loads are associated with smaller streamflows.
Because it is necessary to measure water quality and quantity at the same location in order to estimate the TMDL, streamflow, and water quality gauges approximated from each other were selected for the study. Three and four pairs of stations were selected for Piracicaba and Piranga basins, respectfully, with data from 1997 to 2018. The MG state holds four water quality campaigns throughout the year to account for three different periods: Wet, dry, and intermediary. The wet season goes from January to March, the dry season from July to September, and the intermediary seasons from April to June and October to December. During the wet and dry seasons, 51 water quality parameters are sampled, whereas 19 parameters are sampled in the intermediary seasons. E. coli concentrations are measured in all campaigns.
In this state, the assessment of water availability and water quality occur separately. Therefore, the gauges are not installed in a way to minimize the distance between then. As a consequence, there are limited locations in which it is possible to apply an LDC-like methodology, and the amount of data is frequently limited.
First, the normalized 7-day, 10-year low flow (q7,10), in m3 s−1 km−2, at the selected gauges was calculated per month. This analysis improves the understanding of the hydrological dynamics in the gauges and their potential relationship with water quality. The q7,10 was selected because it is the low flow used by the state of Minas Gerais for water resources management and planning, and it was normalized by the catchment’s upstream drainage area to estimate water production, and therefore, allow their comparison [31].
For the scope of this study, extreme outliers were evaluated and eliminated from the dataset. Inferior and superior extreme samples can be identified using box-plots in combination with Equations (1), and (2) [32], respectively:
L E i n f = Q 1 3 ( Q 3 Q 1 )
L E s u p = Q 3 + 3 ( Q 3 Q 1 )
where LEinf and LEsup are the inferior and superior limits of the box-plot for the identification of extreme outliers.
The seven time-series used in this study, which represent stream flows in the main watercourses of their respective basins, were obtained from the HidroWeb, which is operated by the Brazilian Water Agency (ANA). The FDCs were generated with the Hydrology Plus software [33].

2.3. Load Estimation

The TMDL was calculated for all gauges (LoadTMDL), to understand their water quality dynamics, and then compared with the observed E. coli load monitored in the water quality gauges (LoadOBS). When plotted together, a water quality violation is characterized when LoadOBS is higher than LoadTMDL. LoadTMDL and LoadOBS were calculated according to Equations (3) and (4), respectively:
L o a d T M D L = C M A X Q F D C CF
L o a d O B S = C O B S Q O B S CF
where CMAX is the maximum allowable concentration of E. coli, determined according to Class 2 of the CONAMA regulation, which gives 1000 MPN dL−1; QFDC is the streamflow obtained with the flow duration curve of the gauge; COBS is the observed concentration of E. coli sampled in the water quality gauges; QOBS is the streamflow observed in the day in which the water quality was sampled, and CF is a conversion factor used to estimate the load in MPN day−1 (864,000 × 103).
The Loess smoothing technique [23] was applied to estimate the tendency of LoadOBS samples in relation to the LoadTMDL estimated with the LDC. This technique highlights the tendency of observed data in relation to the maximum observed load in all gauges. Thus, it is possible to understand the expected water quality for different flow permanences. The Loess method consists of a locally weighted regression that, for each value in the x-axis, a fitted value in the y-axis is generated, the so-called xk. According to this method, the fitted values consider all values in the x-axis, but their weight varies according to the proximity of xk. The result of the smoothing technique is that it generates curves with different parameters for each value in the x-axis.
The smoothing analysis helps the visual perception of observed loads in the basin in comparison to the water quality criteria, instead of comparing the sampled points with the LoadTMDL. Thus, in areas where the smoothing curve is above the TMDL, it is likely that the water quality will be impaired for the flows related to that permanence.

2.4. Burr XII Curve

The main novelty of the proposed methodology consists in the identification of the needed load reduction. To quantify the pollutant reduction for each observation and then for all classes of flow, the LDC representing the TMDL was fitted to the Burr XII distribution [24,25]. In this case, it is possible to predict the observed load for all the x-values in which the observed E. coli was sampled, unlike the LDC method proposed by USEPA, whereas the observed load is computed based on a single concentration and the middle permanence flow of each flow regime. The Burr XII is a flexible double-powered distribution and was selected because it has a precise fit for empirical data, thus representing with accuracy the extremes of an FDC, and consequently an LDC, for the scope of this study. The Burr XII distribution can also be adjusted to different streamflow regimes, such as for perennial, intermittent and ephemeral rivers. Equation (5) describes the Burr XII equation,
L o a d p = λ [ { 1 ( p / τ ) β } / β ] α
where Loadp is the observed load arranged in descending order; p is the percentage of exceedance, τ is the cease-to-flow percentile, which is the ratio between the days in which the flow was larger than zero to the total number of observations of the series. In this study, τ = 100 because all rivers are perennial; λ is the scale parameter and β and α are related to the shape of the curve.
If the curve is not adjusted, the quantification of the difference between observed and the maximum allowed E. coli load is computed through approximation, i.e, the method proposed by the USEPA, in which the observed load is set as the multiplication of the median streamflow for each flow range by the 90th percentile concentration of a given flow regime, whereas the maximum allowable is the medium load of that regime. This allowable load is compared to the existing load, in order to identify the ranges by which the pollutant contamination is more expressive [22].
In the present study, the R software was used to fit the permanence values which the observed E. coli was sampled. The equation parameters were determined by an iteration process. The fitted equations allowed to estimate the load related to each percentage in which there were sampled observations. In turn, this allowed the estimation of total E. coli contamination and its percentage observed in each flow regime: High, mid-range and low, corresponding to a percentile of 0%–30%, 30%–70%, and 70%–100%, respectively [34]. The r-squared statistics was used to evaluate the goodness-of-fit of the fitted distribution.

2.5. Estimation of the Needed E. coli Load Reduction

There are different methods to estimate pollutant reduction according to the flow regimes. Here, the reduction, or Relative Difference (RD), is calculated for all observations and compares the observed load with the maximum load, which is obtained with the adjusted curve of the FDC. The RD between the actual sampled load and the fitted load was calculated according to Equation (6),
RD = ( L o a d O B S L o a d M A X L o a d O B S ) 100
where LoadOBS corresponds to the observed load of the water quality gauges, while the LoadMAX represents the maximum allowed load estimated from the fitted TMDL curve.
The Relative Difference (RD) that represents the difference between the observed load and the maximum accepted load considering Class 2 of the CONAMA 357/2005 legislation [35] was calculated for all observed samples. The positive values indicate that there is a need for E. coli load reduction, whereas negative values indicate that the load in the section is not causing the river section to be impaired. In the format of Equation (6), positive values range from 0% to 100%, indicating the percentage of reduction required to meet the water quality criteria. Negative values, on the other hand, can be smaller than negative 100%, depending on the magnitude of LoadOBS in relation to the LoadMAX. Nevertheless, these are always classified as “no-load reduction needed”.
To facilitate the interpretation, the RD was also divided into three groups, representing three classes of flow regime: High (HF), mid-range (MRF), and low flows (LF). For each class, the number of impaired observations was quantified, and then it was possible to estimate the magnitude of the needed reduction for each flow regime with the use of box-plots where the load variation in each class of flow was represented. The RD observations were also divided into seasons, which represent the frequency in which the water quality samples are collected, in order to identify when the pollution is more expressive during the year, and if it matches the expected flow (high, medium or low).
The main limitation of this approach, however, is that the goodness-of-fit of the Burr XII curve to the FDC will determine the validity of the RD estimates, and hence the magnitude of the needed reduction. Additionally, the margin of safety (MOS) and allocation of loads that characterize the final aspects of the LDC methodology proposed by the USEPA were not covered here.

3. Results

3.1. Water Quality and Streamflow Patterns

Prior to the estimation of the E. coli load, the pollutant concentration and streamflow patterns throughout the studied basins were investigated to increase the knowledge about their relationship with the existing load. Figure 2 illustrates the box-plots of the stations involved in the study for the Piracicaba and Piranga basins (Figure 2a,b), as well as the normalized 7-day, 10-year low flow (q7,10) of the selected gauges per month for the Piracicaba and Piranga basins (Figure 2c,d).
In Figure 2, the dashed line represents the maximum pollutant concentration allowed in Class 2 rivers, 1000 MPN dL−1. After the outlier is removed, according to Equations (1) and (2), the number of E. coli samples in the Piracicaba basin was 80, 81 and 79 for gauges RD025 (56610000), RD029 (56659998), and RD031 (56696000), respectively. For the Piranga basin, the number of samples was 72, 73, 114, and 108 for gauges RD001 (56028000), RD007 (56075000), RD013 (56110005), and RD023 (56539000), respectively.
In most stations, the water quality was impaired for more than half of the samples, in particular stations RD025 and RD029 in the Piracicaba basin and station RD013 in the Piranga basin, in which cases more than 75 percent of the data were above the maximum allowed concentration established by the legislation.
With regards to the hydrological patterns, which summarize the historical flowrates from 1989 to 2018, it is clear that these are similar: Streamflow increases from November to March and decreases from March to October. In a study developed in the Piranga basin, streamflow regionalization techniques were applied to estimate the water availability variation throughout the network, and it was observed that, in some points of the hydrography, the difference between dry and wet seasons reaches almost 120% [31].
Given the nature of E. coli data, characterized by a wide range of values, smaller observations commonly get overwhelmed with larger ones. Therefore, log10 transformation was used in the analysis of the data patterns and generation of graphs, whereas the Burr XII was adjusted to untransformed data, and so was performed the RD computation.

3.2. Piracicaba Basin

Figure 3 illustrates the sampled E. coli load, in comparison with the maximum allowable E. coli load (TMDL) for the stations RD025, RD029, and RD031, generated with the development of LDC (Equation (3)). The Loess curve was also represented in the map to indicate the tendency of the sampled E. coli load in comparison to the points representing the TMDL, whereas the grey area surrounding the curve represents the 95% confidence interval for the smoothing curve. The Loess covers all frequencies of permanence for which E. coli data were sampled. Thus, considering that, for some gauges, sample collection could not cover the entire frequency range (0–100%), the Loess may also not cover all the range.
In the Piracicaba basin, there is an impairment for E. coli for all gauges, at some flow range, and E. coli contamination is higher as the drainage area decreases once there is a larger distance between the Loess and the TMDL curves in upstream regions. For gauges RD025 and RD029, a similar pattern is recognized: the tendency line is distant from the TMDL and the percentage of impairment in these gauges is also higher, 90% for both RD025 and RD029 in comparison to 68.8% in gauge RD031.
Figure 4 illustrates the adjusted TMDL curve according to the Burr XII distribution considering the flow duration points that represent the maximum loading capacity, as presented in Figure 3 for the gauges in the Piracicaba basin. It also illustrates a summary of the samples (a box-plot of the observed E. coli) in comparison to the TMDL curve for the high (HF), mid-range (MRF) and low flows (LF).
The box-plots allow recognizing how extensive the impairment of each flow regime is in relation to the TMDL curve. In contrast to Figure 3, Figure 4 allows the pattern of each flow regime to be interpreted and understood. For example, as observed from the tendency curves, gauges RD025 and RD029 have similar E. coli contamination patterns throughout the three regimes of flow (Figure 4a,c). Gauge RD031, on the other hand, have higher E.coli impairment for high flows, decreasing to almost none for mid-range flows (considering the median), and increasing for low flows but in a smaller magnitude, compared to the observed for high flows.
Overall, fewer observations were collected in high flow regimes, which is expected given the smaller permanence in time [36]. In gauges RD025 and RD029, the higher percentage of impairment is observed in mid-range conditions, while for gauge RD031 a smaller impairment was observed for this regime. Equations (7)–(9) illustrate the best-fit TMDL equations for gauges RD025, RD029, and RD031, respectively, as were illustrated in Figure 4.
L o a d R D 025 = 1.11 × 10 13 [ { 1 ( p 100 ) 1.03 × 10 1 } / 1.03 × 10 1 ] 1.89 ( R 2 = 0.995 )
L o a d R D 029 = 2.16 × 10 13 [ { 1 ( p 100 ) 2.34 × 10 2 } / 2.34 × 10 2 ] 1.75 ( R 2 = 0.995 )
L o a d R D 031 = 1.19 × 10 13 [ { 1 ( p 100 ) 2.06 × 10 1 } / 2.06 × 10 1 ] 3.14 ( R 2 = 0.994 )
Figure 5 illustrates the RD between observed E. coli samples and the adjusted TMDL curve for each flow regime. When this information is combined with Figure 4, it is possible to assess the magnitude of the impairment in each flow class, and therefore an estimate for the needed E. coli load reduction.
In all gauges in the Piracicaba basin, the RD, estimated by the median, was between 50% and 100%, indicating the need to reduce this pollutant in all flow regimes. Additionally, it was observed in gauges RD025 and RD029 that there was a higher range between the first and the third quantile for HFs, decreasing for MRFs and again increasing for LFs. It can have been caused by the constant sewage load into the rivers, which have a higher impact on the LF. The variation of RD in HF is probably due to the variability of diffuse sources of pollution.
Figure 6 represents the RD divided into the seasons: January to March (JFM), April to June (AMJ), July to September (JAS) and October to December (OND). Although, the samples were divided into groups of three months, most samples were mainly collected in the first month within the season: January, April, July, and October, whereas fewer were collected in the second (February, May, August and November), and none during the last (March, June, September and December).
As observed in Figure 5, the higher E. coli load reduction in the Piracicaba basin is necessary for the LF. For gauge RD031, although the Loess indicated that the E. coli observations were similar to the TMDL, the needed reduction is still high for all flow regimes. It indicates that, although there is a higher variation within the samples, also observed in Figure 5c, there are more observations that do not match the water quality criteria.
In Figure 6, a smaller median is observed in the first semester of the year that represents mainly wet conditions. In all gauges, a higher RD was observed for the months with lower flows, July and October, where the RD almost reached 100% of reduction.
In RD031, although higher values of reduction are needed, the Loess curve was closer to the TMDL curve than for gauges RD025 and RD029. This happens because the reduction is measured by the median, considering the non-parametric distribution of the E. coli samples. Although there are many samples below the TMDL curve, the samples above it are many more, leading the needed reduction to higher levels. Table 2 illustrates the RD in the Piracicaba basin, according to the different flow regimes and seasons.

3.3. Piranga Basin

Figure 7 illustrates the sampled E. coli load in relation to the points representing the maximum allowable E. coli load (TMDL) in the stations RD001, RD007, RD013, and RD023, generated with the development of LDCs. The Loess curve indicates the tendency of sampled E. coli loads in comparison to the points representing the TMDL in the Piranga basin, within the 95% confidence interval.
In the Piranga basin, a similar pattern is observed for gauges RD001 and RD007 (Figure 7a,b). The impairment is more frequent in streamflow classes of higher flows, and it decreases for mid-flows. The main difference between these gauges is that the tendency line of the observed data is constant for low flows for the RD001 gauge, whereas it continues to decrease for the RD007 gauge. Nevertheless, within the 95% confidence interval (the gray area surrounding the smoothing curve) the TMDL and the Loess curve of the observed data are only different for high flows, from 0–25%.
The RD023 gauge, which is located near the mouth of the basin, is associated with the larger drainage area when compared to the other gauges and is inside a protected area. Unlike the other gauges, where the most sampled points were above the TMDL curve, in RD023 the impairment with the water quality standards is more expected from the high flow conditions to the mid-flow conditions. This situation might have happened because the constant load, generated from cities, is smaller than the loading capacity of the river in low flow conditions. In this condition, it is also expected that the E. coli load from upland locations in the basin has decayed before reaching this gauge. On the other hand, during high flows, the streamflow is higher, resulting in a higher velocity of the contaminant from bigger cities, as Ponte Nova, which caused the E. coli to reach gauge RD023.
Figure 8 illustrated the adjusted TMDL curve according to the Burr XII distribution considering the flow duration points that represent the maximum loading capacity for the Piranga basin, as well as box-plots of the observed E. coli, sampled points and the TMDL curve for the assessed flow regimes.
Considering the median of all observations collected in each flow regime, the pattern observed in Figure 7 persists in Figure 8. The loads for high and mid-flow conditions are similar for gauges RD001 and RD007, but for low flows they tend to maintain the same pattern of mid-flow conditions just for gauge RD001 and decrease for gauge RD007. For gauge RD013, a reduction in E. coli contamination is observed as the flow decreases. Finally, for gauge RD023 the main contamination is observed when flow conditions are high, whereas 75% of all observations are below the TMDL curve for mid-range and low flows.
Equations (10)–(13) illustrate the best-fit TMDL equation for gauges RD001, RD007, RD013, and RD023, respectively.
L o a d R D 001 = 1.36 × 10 13 [ { 1 ( p 100 ) 7.74 × 10 2 } / 7.74 × 10 2 ] 1.09   ( R 2 = 0.993 )
L o a d R D 007 = 5.26 × 10 13 [ { 1 ( p 100 ) 2.23 × 10 1 } / 2.23 × 10 1 ] 7.37   ( R 2 = 0.996 )
L o a d R D 013 = 7.07 × 10 13 [ { 1 ( p 100 ) 4.75 × 10 1 } / 4.75 × 10 1 ] 5.82   ( R 2 = 0.996 )
L o a d R D 023 = 1.55 × 10 14 [ { 1 ( p 100 ) 1.38 × 10 1 } / 1.38 × 10 1 ] 8.81 × 10 1   ( R 2 = 0.989 )
Figure 9 depicts the relative difference (RD) between observed E. coli samples and the adjusted TMDL curve for each flow regime in the Piranga basin. As with Figure 5 and Figure 6, the combined use of Figure 9 and Figure 10 allows the magnitude of impairment to be assessed that is in each flow class and in each gauge, giving an idea of the need for E. coli load reduction.
In gauge RD001, there is a higher need for reduction in the LFs, whereas in gauge RD007 the reduction needed in the LFs and HFs is similar. The RD is consistently high in gauge RD013. This is a consequence of sewage load from Ponte Nova. In gauge RD023, a higher reduction is needed in the HFs. The variance in the gauges is a consequence of a smaller influence of point sources of pollution, meaning the sewage load in the rivers, as was observed for the Piracicaba basin. The negative values are obtained when the observed load is smaller than the TMDL, and higher negative magnitudes indicate a smaller LoadOBS (Equation (6)), which explains why, for some gauges, the box-plot is not entirely on the plot (i.e., Figure 9d). Thus, it is possible to affirm that lower values of E. coli load (in relation to the TMDL curve) were observed in gauge RD023, followed by gauge RD007 and then RD001. Gauge RD013 presented only positive values.
Figure 10 illustrates the RD variation in the Piranga basin throughout the seasons. This figure shows higher load variances, with larger differences between the first and third quantiles, unlike the situation observed in the Piracicaba basin. As perceived from previous results, a higher pollutant reduction is needed in gauge RD013, which showed higher RD for all classes of flow throughout the year. It also shows that a higher E. coli reduction is needed in dry months (July and October) in comparison to wet months (January and April). In gauge RD023, the RD is positive only in wet months, indicating a need for load reduction in diffuse sources of pollution. Table 3 illustrates the RD in the Piranga basin, according to the different flow regimes and different seasons.
In general, the variation is more significant in this basin than in the Piracicaba basin. Besides, no-load reduction is needed in gauge RD023 during the dry season, but this observation is not valid for the other gauges.

4. Discussion

4.1. E. coli Concentration and Streamflow Patterns

The E. coli concentration in the gauges is directly related to the proximity of urban centers, where the population is mostly concentrated [37]. In the Piracicaba basin, the stations RD025 and RD029 were installed in city surroundings, whereas the gauge RD031 is within the extent of Cel. Fabriciano city, but was installed upstream of its urban center. In the latter case, the station did not receive the raw (untreated) sewage generated in one of the most populated cities in the basin. This fact explains the larger E. coli concentration upstream of the Piracicaba basin in relation to downstream.
Additionally, although the impairment at gauge RD031 is greater than 50% (Figure 2a), if the confidence interval of the Loess curve is considered in relation to the TMDL curve, these overlap, resulting in an impairment percentage of about 20%, mainly observed in high flows. In this gauge, the E. coli contamination is probably attributed to runoff, cattle grazing, or/and the sewage from upstream cities that arrive faster in the gauge during events of high flows, not allowing the process of bacterial decay.
For the Piranga basin, the stations RD001 and RD007 are located close to the urban center of the cities, while station RD013 is located downstream from the city center. For the first stations, most sewage generated by the city is likely loaded after their location, which would result in sampled E. coli values that do not reflect the water quality in the cities. For station RD013, most sewage load generated by the city likely reaches the gauge, although some of the loaded E. coli can decay from the moment the sewage reaches the river. This affirmation is supported by the fact that in this gauge the smoothing curve of the observed samples of E. coli is distant from the TMDL curve for all permanences of streamflow. This situation is probably a consequence of the untreated sewage load from Ponte Nova, one of the most populated cities of this basin, with over 53,000 inhabitants.
The station RD023, on the other hand, is located in a preserved state park, the Doce State Park, away from the urban centers of both cities where it is located. Although a higher E. coli concentration is observed near the cities, it cannot be ruled out that this region is a pole of swine farming activities, which also contribute to a high concentration of pathogens, due to high percentages of waste deposition that are not respecting the environmental regulations [27].
Apart from station RD013, the tendency of the median is to reduce as the drainage area increase. The same pattern was observed in other basins, in which higher E. coli loading was observed to be higher in watersheds with smaller drainage areas and with a higher population density [38]. It is suggested that this situation probably arises because, in smaller basins, the pollutant source is likely to be closer to the streams, resulting in a faster transport rate, thus reducing the impact of E. coli decay. Similarly, in the surroundings of Chinese counties and evidenced that the arrangements of the cities play an important role in water quality. Thus, government initiatives, such as structural reforms and environmental regulations, impact water quality [37]. An example of the importance of government actions in order to improve water quality happened in the USA when a detergent ban and a total P discharge was imposed in order to reduce algal blooms [39].

4.2. E. coli Load Patterns in the Piracicaba Basin

As expected, a higher RD was observed for the seasons with lower flows, which can receive closer attention in the establishment of measurements to reduce the pollution load. In a study developed in urban catchments in Florida (USA), higher E. coli contamination was observed during the summer, with wet and hot conditions [38]. This difference in relation to the results presented for this basin is probably due to the difference in scale, which is higher here, resulting in a delay in the transport velocity, lessening the effect of the diffuse pollution.
In a different study, the dynamics of fecal indicators were analyzed according to different conditions of land use, season and water chemistry. It was observed that sometimes the indicators had different patterns throughout the basin [40]. E. coli, in agricultural basins, maintained high concentrations throughout the year, whereas, in urban cities, it tended to increase during the summer season and is smaller during the winter.
The combination of these analyses can improve the current water quality assessment method where only the pollutant concentration is monitored. For instance, it is possible to estimate the pollution pattern, i.e., the regime where it is more recurrent and when there are more samples away from water quality criteria. In gauge RD025, for instance, the percentage of impairment is higher for greater flows (almost 94% of impairment during HF). However, as it is observed in Figure 5a, the magnitude of pollution is higher for low flows, as the median RD of low flows in this class is more distant to x-axis (RD equivalent to zero). As a whole, E. coli contamination in the basin is a concern and is the reason why all classes of flow and sources of pollution (point and non-point) should be addressed.

4.3. E. coli Load Patterns in the Piranga Basin

In the Piranga basin, the total impairment is ampler than in the Piracicaba basin, varying from over 95% in gauge RD013 to below 10% in gauge RD023. The impairment is larger in RD001 than in RD007, while this was not possible to detect in the smoothing curve. In a study developed in different scales in Texas: Field, small watershed and river basin, E. coli concentrations were evaluated with different agricultural management and land use [41]. In this study, the authors also observed that the primary E. coli sources in rural basins are wildlife, streamflow resuspension, failing wastewater treatment plants and animal feeding operations. Additionally, in the river courses evaluated, E. coli contamination tends to be higher upstream and decreases downstream. Throughout the Piranga basin, the main sources of pollution are likely caused by cultivated and grazed lands and swine production [27], except by gauge RD013, in which contamination is caused by failing and non-existing wastewater treatment plants.
In the Piranga basin, there is a need for E. coli load reduction in the HFs in all gauges. However, in gauges RD001 and RD013, the magnitude of reduction is higher for the LFs (Table 3). This can be used as evidence to address these flows first when summed with the fact that low flows can represent critical conditions in the context of water supply.
In gauge RD001, E. coli load reduction is also needed in all seasons, but unlike in gauge RD013, a higher variation throughout the seasons is observed.
With the information of the flow regimes with the must recurrence of impairment and the months in which the problems are more aggravated, it is easier for the statewide environmental agencies to identify potential sources of pollution and thus impose regulations to minimize this load.

4.4. E. coli Load Variation in Piracicaba and Piranga Basins

In order to understand the dynamics of the observed E. coli load throughout the Doce river basin, the smoothing curves of the gauges within the Piracicaba and Piranga sub-basins were plotted together in Figure 11.
Given the smaller size of the Piracicaba basin in comparison to the Piranga, there is not a significant load variation throughout the former basin (Figure 11a). For instance, with 95% confidence, the load in gauge RD029 is equivalent to the E. coli loads in gauges RD025 and RD031 at some flow permanence. In this basin, it is also possible to observe a steeper decrease in the E. coli load in the gauge RD031 as the permanence increases in relation to the other gauges in the basin.
In the Piranga basin (Figure 11b), a broader load variation is observed. The load in gauge RD013 is higher in all permanences of flow, as observed in other studies, where the agricultural basin presented higher and constant E. coli concentrations, in relation to the urban [40]. In gauge RD001, the load in low flow conditions is higher than the loads in the gauges RD007 and RD023 for higher permanences, even considering its smaller drainage area. The difference in the load of gauge RD013, with respect to the other gauges throughout the permanences of streamflow, can be attributed to the fact that it is located in the city with a higher amount of inhabitants, over 53,000, in relation to the others, in which the population is under 7,000 inhabitants.
The Loess curve does not present a continuum decrease in its load magnitude, as expected for the TMDL curve, because it was generated for the sampled observations, which reflected the water quality of these gauges. Therefore, the decrease in permanence does not imply in the reduction of the load, once it will follow the water quality observed in each permanence of flow and, for this reason, the Loess curve can have positive and negative slopes throughout the streamflows.
On that basis, the source control actions to be taken in order to improve water quality in the studied basins should prioritize regions where it could have more impact in the basin. In the Piracicaba basin, the gauges upstream (RD025 and RD029) are likely good candidates given the smaller drainage area and a probable high percentage of point sources. In the Piranga basin, on the other hand, the gauge where the adoption of pollution control practices would have a higher impact is the RD013. In the listed gauges, the main advantage is the existence of point sources of pollution, which are more feasible to identify and further control.

5. Conclusions

This study presents an alternative method to estimate E. coli load reduction in different locations of Piranga and Piracicaba basins, considering different flow regimes and seasons, using the concept of Load Duration Curves (LDC) and Total Maxim Daily Loads (TMDL). The results made evident that the loading of raw sewage in the rivers is one of the leading causes of E. coli contamination in all flow regimes for both basins. Thus, the first practice in the basins to control E. coli pollution should be the installation of wastewater treatment plants.
In the Piracicaba basin, the tendency is that the observed load is higher than the TMDL curve for all flow regimes, with a higher need for load reduction in low flows, in upstream gauges. On gauge RD031, on the other hand, located downstream of the basin, although there is an impairment of almost 70%, the TMDL load is equal to the observed load for about 80% of the streamflow permanence, considering the smoothing (Loess) curve.
For the Piranga basin, there is a higher variation in the E. coli contamination. Nevertheless, more efforts should focus on controlling the contamination in the higher flows (HF), because the RD was positive for all gauges under this condition, and the minimum percentage of impairment was over 65%.
The study proposes a different way to monitor water quality in Brazilian waters. Nowadays, the primary method used to measure water quality is through the measurement of pollutant concentration, not considering the streamflow related to that concentration. However, in order to improve water quality, it is essential to monitor both aspects, namely water quality and quantity, to propose techniques that are appropriate for the watershed and will control pollution.
The study brought different analyses that can be explored to improve water quality in the studied and other basins. For instance, with this methodology, it is possible to assess water quality and propose periodical goals in order to achieve a certain water standard. The method can also help define the needed efficiency in wastewater treatment plants, and thus, the definition of proper technology. Also, the identification of potential pollution sources can be facilitated with the seasonality approach. This type of approach can aid in the understanding of fate, transport, and survival of E. coli, which is a pollutant of many basins. Nevertheless, the methodology proposed can be applied to different contaminants.

Author Contributions

L.d.O.S., A.C.B. and F.F.P. conceived and designed the study. L.d.O.S. collected the data, performed the computations, analyzed the results and wrote the original draft of the manuscript. A.C.B. supervised the research work and contributed to the interpretation of the results and manuscript revision. M.C.d.M. and F.F.P. also contributed with the result interpretation and elaboration of the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Council for Scientific and Technological Development (CNPq) grant number 165565/2017-9, Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), finance code: 001, and Minas Gerais Research Foundation (FAPEMIG) grant number PPM-00911-15.

Acknowledgments

We are thankful to Fernando Pacheco for his constructive and critical comments on this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, M.; Shao, Q.; Zhang, L.; Chiew, F.H. A New Regionalization Approach and Its Application to Predict Flow Duration Curve in Ungauged Basins. J. Hydrol. 2010, 389, 137–145. [Google Scholar] [CrossRef]
  2. Pruski, F.F.; de Araújo Nunes, A.; Pruski, P.L.; del Giudice Rodriguez, R. Improved Regionalization of Streamflow by Use of the Streamflow Equivalent of Precipitation as an Explanatory Variable. J. Hydrol. 2013, 476, 52–71. [Google Scholar] [CrossRef]
  3. Nejadhashemi, A.P.; Woznicki, S.A. Comparison of Four Models (STEPL, PLOAD, L-THIA, AND SWAT) in Simulating Sediment, Nitrogen, and Phosphorus Loads and Pollutant Source Areas. Trans. ASABE 2011, 54, 875–890. [Google Scholar] [CrossRef]
  4. Wu, Y.; Chen, J. Investigating the Effects of Point Source and Nonpoint Source Pollution on the Water Quality of the East River (Dongjiang) in South China. Ecol. Indic. 2013, 32, 294–304. [Google Scholar] [CrossRef]
  5. Shrestha, N.K.; Leta, O.T.; De Fraine, B.; Garcia-Armisen, T.; Ouattara, N.K.; Servais, P.; van Griensven, A.; Bauwens, W. Modelling Escherichia Coli Dynamics in the River Zenne (Belgium) Using an OpenMI Based Integrated Model. J. Hydroinf. 2014, 16, 354–374. [Google Scholar] [CrossRef] [Green Version]
  6. Beskow, S.; Mello, C.R.D.; Norton, L.D. Development, Sensitivity and Uncertainty Analysis of LASH Model. Sci. Agric. 2011, 68, 265–274. [Google Scholar] [CrossRef] [Green Version]
  7. Paiva, R.C.D.; Collischonn, W.; Bonnet, M.P.; De Gonçalves, L.G.G. On the Sources of Hydrological Prediction Uncertainty in the Amazon. Hydrol. Earth Syst. Sci. 2012, 16, 3127–3137. [Google Scholar] [CrossRef] [Green Version]
  8. De Bressiani, D.A.; Gassman, P.W.; Fernandes, J.G.; Garbossa, L.H.P.; Srinivasan, R.; Bonumá, N.B.; Mendiondo, E.M. A Review of Soil and Water Assessment Tool (SWAT) Applications in Brazil: Challenges and Prospects. Int. J. Agric. Biol. Eng. 2015, 8, 1–27. [Google Scholar] [CrossRef]
  9. De Mello, C.R.; Norton, L.D.; Pinto, L.C.; Beskow, S.; Curi, N. Agricultural Watershed Modeling: A Review for Hydrology and Soil Erosion Processes. Ciência e Agrotecnologia 2016, 40, 7–25. [Google Scholar] [CrossRef] [Green Version]
  10. Liao, H.; Krometis, L.A.H.; Cully Hession, W.; Benitez, R.; Sawyer, R.; Schaberg, E.; von Wagoner, E.; Badgley, B.D. Storm Loads of Culturable and Molecular Fecal Indicators in an Inland Urban Stream. Sci. Total Environ. 2015, 530–531, 347–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Bonta, J.V.; Cleland, B. Incorporating Natural Variability, Uncertainty, and Risk into Water Quality Evaluations Using Duration Curves. J. Am. Water Resour. Assoc. 2003, 39, 1481–1496. [Google Scholar] [CrossRef]
  12. Johnson, S.L.; Whiteaker, T.; Maidment, D.R. A Tool for Automated Load Duration Curve Creation. J. Am. Water Resour. Assoc. 2009, 45, 654–663. [Google Scholar] [CrossRef]
  13. Cleland, B.R. TMDL Development from the “Bottom up”—Part III: Duration Curves and Wet-Weather Assessments. Proceedings of the Water Environment Federation. January 2003, pp. 1740–1766. Available online: https://www.researchgate.net/publication/228822472_TMDL_Development_from_the_Bottom_Up-_PART_III_Durations_Curves_and_Wet-Weather_Assessments (accessed on 14 March 2020).
  14. Teague, A.; Bedient, P.B.; Guven, B. Targeted Application of Seasonal Load Duration Curves Using Multivariate Analysis in Two Watersheds Flowing into Lake Houston. J. Am. Water Resour. Assoc. 2011, 47, 620–634. [Google Scholar] [CrossRef]
  15. Stallard, M.A.; Otter, R.R.; Winesett, S.; Barbero, M.; Bruce, M.; Layton, A.; Bailey, F.C. A Watershed Analysis of Seasonal Concentration- and Loading-Based Results for Escherichia Coli in Inland Waters. Bull. Environ. Contam. Toxicol. 2016, 97, 838–842. [Google Scholar] [CrossRef] [PubMed]
  16. Cleland, B. TMDL Development from the “Bottom up”—Part II: Using Duration Curves to Connect the Pieces. In Proceedings of the National TMDL Science and Policy 2002 Specialty Conference, Phoenix, AZ, USA, 13–16 November 2002. [Google Scholar]
  17. Cunha, D.G.F.; Calijuri, M.C. Análise Probabilística de Ocorrência de Incompatibilidade Da Qualidade Da Água Com o Enquadramento Legal de Sistemas Aquáticos—Estudo de Caso Do Rio Pariquera-Açu (SP). Eng. Sanitária e Ambient. 2010, 15, 337–346. [Google Scholar] [CrossRef] [Green Version]
  18. Cunha, F.D.G.; Calijuri, C.; Mendiondo, E.M. Integração Entre Curvas de Permanência de Quantidade e Qualidade Da Água Como Uma Ferramenta Para a Gestão Eficiente Dos Recursos Hídricos. Eng. Sanitária e Ambient. 2012, 17, 369–376. [Google Scholar] [CrossRef] [Green Version]
  19. Teodoro, A.; Ide, C.N.; Ribeiro, M.L.; Broch, S.A.O.; Silva, J.B.D. Implementação Do Conceito Capacidade de Diluição de Efluentes No Modelo de Qualidade Da Água QUAL-UFMG: Estudo de Caso No Rio Taquarizinho (MS). Eng. Sanit. e Ambient. 2013, 18, 275–288. [Google Scholar] [CrossRef] [Green Version]
  20. Marques, L.O.D.A.; Taffarello, D.; Calijuri, M.D.C.; Mendiondo, E.M.; Ferreira, M.D.S.; Cunha, D.G.F. Phosphorus and Thermotolerant Coliforms’ Loads in Brazilian Watersheds with Limited Data: Considerations on the Integrated Analysis of Water Quality and Quantity. Rev. Bras. Recur. Hídricos 2019, 24, 1–13. [Google Scholar] [CrossRef] [Green Version]
  21. United States Environmental Protection Agency. An Approach for Using Load Duration Curves in the Development of TMDLs; USEPA: Washington, DC, USA, 2007.
  22. United States Environmental Protection Agency. Options for the Expression of Daily Loads in TMDLs; USEPA: Washington, DC, USA, 2007.
  23. Cleveland, W.S. Robust Locally Weighted Regression and Smoothing Scatterplots. J. Am. Stat. Assoc. 1979, 74, 829–836. [Google Scholar] [CrossRef]
  24. Shao, Q.; Wong, H.; Xia, J.; Ip, W.C. Models for Extremes Using the Extended Three-Parameter Burr XII System with Application to Flood Frequency Analysis. Hydrol. Sci. J. 2004, 49, 685–701. [Google Scholar] [CrossRef]
  25. Shao, Q.; Zhang, L.; Chen, Y.D.; Singh, V.P. A New Method for Modelling Flow Duration Curves and Predicting Streamflow Regimes under Altered Land-Use Conditions. Hydrol. Sci. J. 2009, 54, 606–622. [Google Scholar] [CrossRef]
  26. Fernandes, G.W.; Goulart, F.F.; Ranieri, B.D.; Coelho, M.S.; Dales, K.; Boesche, N.; Bustamante, M.; Carvalho, F.A.; Carvalho, D.C.; Dirzo, R.; et al. Deep into the Mud: Ecological and Socio-Economic Impacts of the Dam Breach in Mariana, Brazil. Nat. e Conserv. 2016, 14, 35–45. [Google Scholar] [CrossRef]
  27. De Melo, M.C.; de Mendes, R.S.; da Fonseca, M.; Aquino, F.A.; Rodrigues, A.F. Environmental Monitoring Planning in River Basins: The Case of Pig Farming in River Piranga Basin. Caminhos Geogr. 2017, 18, 454–471. [Google Scholar] [CrossRef] [Green Version]
  28. Nunes, L.G.D.P.; Oliveira, M.D.V.; de Souza, A.A.; Lopes, L.D.F.; Dias, P.C.E.S.; Nogueira, G.B.; Souza, M.A.A.D. Water Quality Comparison between a Supply Network and Household Reservoirs in One of the Oldest Cities in Brazil. Int. J. Environ. Health Res. 2019, 29, 173–180. [Google Scholar] [CrossRef]
  29. Kim, J.; Engel, B.A.; Park, Y.S.; Theller, L.; Chaubey, I.; Kong, D.S.; Lim, K.J. Development of Web-Based Load Duration Curve System for Analysis of Total Maximum Daily Load and Water Quality Characteristics in a Waterbody. J. Environ. Manag. 2012, 97, 46–55. [Google Scholar] [CrossRef]
  30. Benham, B.L.; Baffaut, C.; Zeckoski, R.W.; Mankin, K.R.; Pachepsky, Y.A.; Sadeghi, A.M.; Brannan, K.M.; Soupir, M.L.; Habersack, M.J. Modeling Bacteria Fate and Transport in Watersheds to Support TMDLs. Trans. ASABE 2006, 49, 987–1002. [Google Scholar] [CrossRef] [Green Version]
  31. De Serrano, L.O.; Ribeiro, R.B.; Borges, A.C.; Pruski, F.F. Low-Flow Seasonality and Effects on Water Availability throughout the River Network. Water Resour. Manag. 2020. [Google Scholar] [CrossRef]
  32. Oliveira, F.A. Procedimentos Para Aprimorar a Regionalização de Vazões: Estudo de Caso Na Bacia Do Rio Grande. Master’s Thesis, Federal University of Vicosa, Viçosa, Brazil, April 2008; p. 106. [Google Scholar]
  33. De Sousa, J.R.C. Hydrology Plus: Sistema Computacional Para Estudos Meteorológicos e Hidrológicos; Universidade Federal de Viçosa: Viçosa, Brazil, 2017. [Google Scholar]
  34. Shen, J.; Zhao, Y. Combined Bayesian Statistics and Load Duration Curve Method for Bacteria Nonpoint Source Loading Estimation. Water Res. 2010, 44, 77–84. [Google Scholar] [CrossRef]
  35. CONAMA, Environment National Council (Brazil). Resolution 357: Provides the Classification of Water Bodies and Environmental Guidelines for Its Setting as Well as Establishes the Conditions and Standards of Effluent Discharge, and Give Olher Measures; Ministry of the Environment: Brasilia, Brazil, 2005; p. 27.
  36. Babbar-Sebens, M.; Karthikeyan, R. Consideration of Sample Size for Estimating Contaminant Load Reductions Using Load Duration Curves. J. Hydrol. 2009, 372, 118–123. [Google Scholar] [CrossRef]
  37. Zhang, Q.; Yang, L.; Song, D. Environmental Effect of Decentralization on Water Quality near the Border of Cities: Evidence from China’s Province-Managing-County Reform. Sci. Total Environ. 2020, 708, 1–11. [Google Scholar] [CrossRef]
  38. Chin, D.A. Linking Pathogen Sources to Water Quality in Small Urban Streams. J. Environ. Eng. 2010, 136, 249–253. [Google Scholar] [CrossRef]
  39. Lebo, M.E.; Paerl, H.W.; Peierls, B.L. Evaluation of Progress in Achieving TMDL Mandated Nitrogen Reductions in the Neuse River Basin, North Carolina. Environ. Manag. 2012, 49, 253–266. [Google Scholar] [CrossRef] [PubMed]
  40. Badgley, B.D.; Steele, M.K.; Cappellin, C.; Burger, J.; Jian, J.; Neher, T.P.; Orentas, M.; Wagner, R. Fecal Indicator Dynamics at the Watershed Scale: Variable Relationships with Land Use, Season, and Water Chemistry. Sci. Total Environ. 2019, 697, 134113. [Google Scholar] [CrossRef]
  41. Harmel, R.D.; Karthikeyan, R.; Gentry, T.; Srinivasan, R. Effects of Agricultural Management, Land Use on E. Coli Concentrations in Runoff and Streamflow. Trans. ASABE 2010, 53, 1833–1841. [Google Scholar] [CrossRef]
Figure 1. Identification and location of the studied basins, with the plot of streamflow and water quality gauges selected for this study.
Figure 1. Identification and location of the studied basins, with the plot of streamflow and water quality gauges selected for this study.
Water 12 00811 g001
Figure 2. E. coli concentrations in the (a) Piracicaba, and (b) Piranga basins (the dashed line represents the maximum allowed concentration in Class 2 rivers), and monthly streamflow variation in (c) Piracicaba and (d) Piranga basins.
Figure 2. E. coli concentrations in the (a) Piracicaba, and (b) Piranga basins (the dashed line represents the maximum allowed concentration in Class 2 rivers), and monthly streamflow variation in (c) Piracicaba and (d) Piranga basins.
Water 12 00811 g002
Figure 3. Comparison between load duration curves (TMDL) and Loess curve (representing the tendency of E. coli samples) in the Piracicaba river basin for gauges: (a) RD025/56610000, (b) RD029/56659998 and (c) RD031/56696000.
Figure 3. Comparison between load duration curves (TMDL) and Loess curve (representing the tendency of E. coli samples) in the Piracicaba river basin for gauges: (a) RD025/56610000, (b) RD029/56659998 and (c) RD031/56696000.
Water 12 00811 g003
Figure 4. E. coli samples and box-plot per class of flow in relation to the TMDL curve for gauges: (a,b) RD025/56610000, (c,d) RD029/56659998 and (e,f) RD031/56696000. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Figure 4. E. coli samples and box-plot per class of flow in relation to the TMDL curve for gauges: (a,b) RD025/56610000, (c,d) RD029/56659998 and (e,f) RD031/56696000. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Water 12 00811 g004
Figure 5. Percentage of reduction per class of flow for gauges (a) RD025/56610000, (b) RD029/56659998, and (c) RD031/56696000. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Figure 5. Percentage of reduction per class of flow for gauges (a) RD025/56610000, (b) RD029/56659998, and (c) RD031/56696000. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Water 12 00811 g005
Figure 6. Percentage of reduction per season of flow for gauges (a) RD025/56610000, (b) RD029/56659998, and (c) RD031/56696000. JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.
Figure 6. Percentage of reduction per season of flow for gauges (a) RD025/56610000, (b) RD029/56659998, and (c) RD031/56696000. JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.
Water 12 00811 g006
Figure 7. Comparison between load duration curves (TMDL) and Loess curve (representing the tendency of E. coli samples) in the Piranga river basin: Gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and (d) 56539000/RD023.
Figure 7. Comparison between load duration curves (TMDL) and Loess curve (representing the tendency of E. coli samples) in the Piranga river basin: Gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and (d) 56539000/RD023.
Water 12 00811 g007
Figure 8. E. coli samples and box-plots per class of flow in relation to the TMDL curve for gauge (a,b) 56028000/RD001, (c,d) 56075000/RD007, (e,f) 56110005/RD013 and (g,h) 56539000/RD023. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Figure 8. E. coli samples and box-plots per class of flow in relation to the TMDL curve for gauge (a,b) 56028000/RD001, (c,d) 56075000/RD007, (e,f) 56110005/RD013 and (g,h) 56539000/RD023. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Water 12 00811 g008
Figure 9. Percentage of reduction per class of flow for gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and, (d) 56539000/RD023. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Figure 9. Percentage of reduction per class of flow for gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and, (d) 56539000/RD023. HF, high flows, MRF, Mid-range flows, LF, low-flows.
Water 12 00811 g009
Figure 10. Percentage of reduction per season for gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and, (d) 56539000/RD023. JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.
Figure 10. Percentage of reduction per season for gauge (a) 56028000/RD001, (b) 56075000/RD007, (c) 56110005/RD013 and, (d) 56539000/RD023. JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.
Water 12 00811 g010
Figure 11. Tendency curves (Loess) of observed points in the gauges within the (a) Piracicaba and (b) Piranga basins.
Figure 11. Tendency curves (Loess) of observed points in the gauges within the (a) Piracicaba and (b) Piranga basins.
Water 12 00811 g011
Table 1. List of water quality and streamflow gauges used in the study, with the upstream drainage area. Information is also provided for the cities where the gauges are located, namely city name and their population.
Table 1. List of water quality and streamflow gauges used in the study, with the upstream drainage area. Information is also provided for the cities where the gauges are located, namely city name and their population.
BasinQuality GaugeStreamflow GaugeDrainage Area (km2)CityPopulation
PiracicabaRD025566100001160Rio Piracicaba11,614
RD029566599983060Nova Era15,837
RD031566960005270Timóteo/Cel. Fabriciano85,888 and 106,945
PirangaRD001560280001400Piranga6,156
RD007560750004260Porto Firme5,081
RD013561100056230Ponte Nova53,169
RD0235653900015,900Córrego Novo/Marliéria2,020 and 2,924
Table 2. Needed E. coli reduction (RD) per regime of flow and per season in the Piracicaba basin.
Table 2. Needed E. coli reduction (RD) per regime of flow and per season in the Piracicaba basin.
GaugeHF (%)MRF (%)LF (%)
56610000/RD02585.394.499.7
56659998/RD02980.794.499.4
56696000/RD03190.095.399.9
JFM (%)AMJ (%)JAS (%)OND (%)
56610000/RD02590.593.8598.399.6
56659998/RD02986.694.4 98.1 97.9
56696000/RD03191.196.1 98.9 99.8
Note: HF, high flows, MRF, Mid-range flows, LF, low-flows, JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.
Table 3. Needed E. coli reduction (RD) per regime of flow and per season in the Piranga basin.
Table 3. Needed E. coli reduction (RD) per regime of flow and per season in the Piranga basin.
GaugeHF (%)MRF (%)LF (%)
56028000/RD00167.558.593.5
56075000/RD00765.9-62.4
56110005/RD01393.594.497.4
56539000/RD02339.1--
JFM (%)AMJ (%)JAS (%)OND (%)
56028000/RD00174.277.784.284.0
56075000/RD00759.439.862.459.4
56110005/RD01395.795.697.396.4
56539000/RD02323.99.1--
Note: HF, high flows, MRF, Mid-range flows, LF, low-flows, JFM, January to March; AMJ, April to June; JAS, July to September; OND, October to December.

Share and Cite

MDPI and ACS Style

Serrano, L.d.O.; Borges, A.C.; Pruski, F.F.; Melo, M.C.d. A New Approach to Use Load Duration Curves to Evaluate Water Quality: A Study in the Doce River Basin, Brazil. Water 2020, 12, 811. https://doi.org/10.3390/w12030811

AMA Style

Serrano LdO, Borges AC, Pruski FF, Melo MCd. A New Approach to Use Load Duration Curves to Evaluate Water Quality: A Study in the Doce River Basin, Brazil. Water. 2020; 12(3):811. https://doi.org/10.3390/w12030811

Chicago/Turabian Style

Serrano, Ligia de Oliveira, Alisson Carraro Borges, Fernando Falco Pruski, and Marília Carvalho de Melo. 2020. "A New Approach to Use Load Duration Curves to Evaluate Water Quality: A Study in the Doce River Basin, Brazil" Water 12, no. 3: 811. https://doi.org/10.3390/w12030811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop