Wastewater treatment systems selection inside watersheds by using multiobjective analysis

In order to determine required wastewater treatment plants (WWTP’s) efficiencies, several methodologies have been proposed to assist the sewage treatment systems process selection, usually aiming to achieve economic objectives. However, in water resources problems solving, search involves multiple and conflicting objectives. This work used, for Pardo’s river basin (Espírito Santo State, Brazil), water quality simulation model, optimization technique and multiobjective analysis for selecting sewage treatment systems, taking into account Brazilian environmental effluents Biochemical Oxygen Demand standards, imposed by CONAMA Resolution 430/2011. After proposing an optimization model that allowed the estimation of minimum organic matter removal efficiencies, a multiobjective analysis was applied by using the ELECTRE III method, which selected sewage treatment systems for localities responsible for domestic wastewater disposal in Pardo’s river basin. The results indicated that keeping the effluents standards can produce treatment plants overestimation and inadequate allocated financial resources distribution for sewage treatment plants implementation. Treatment systems such as Submerged Aerated Biofilters (with Nitrification), Wetlands and Facultative Ponds were the most selected, considering different contexts appreciated in this work.


INTRODUCTION
In the accelerated urban development context in various parts of the world, wastewater treatment plants (WWTP) are key technologies for meeting restrictive environmental and public health standards.Required WWTP's treatment levels depend on receiving bodies self-purification capacities, water uses downstream discharges points and discharges characteristics and conditions (JORDÃO; PESSOA, 2014).However, scarcity of water, increasing number of alternatives available for effluents treatment and growing pressure related with projects technical, environmental, economic and social balances are aspects that turn difficult WWTPs installing process decisions (CASTILLO et al., 2017), especially if analyzed within a watershed where different water bodies potential polluters are involved.
In Brazil, beyond National Water Resources Policy (BRASIL, 1997) watersheds water resources management implementation, the Basic Sanitation Law (BRASIL, 2007) adopted watersheds as reference units for actions planning.However, multiple combination of factors, related with each geographic space, involving physiographic, institutional, socio-cultural and economic aspects, requires important analytical and methodological efforts to address the sanitation issue with focus on water resources protection (ANA, 2017).
In general, several methodologies have been developed for producing tools to assist sewage treatment systems selection processes (SOUZA;FORSTER, 1996;SOUZA;CORDEIRO;SILVA, 2009), differentiating themselves, as Hamouda, Anderson and Huck (2009) suggest, in three different approaches: (1) technical design, which considers information on treatment alternatives efficiency and effectiveness; (2) technical and economic analysis, which usually uses cost as a decision variable; and (3) system analysis, or "appropriate technology analysis", through which technological choice represents different technical, economic, social, environmental and cultural criteria point of convergence and harmonization.
Associated with the technical and economic approach, Burn and Yulianti (2001), Cho, Seok Sung and Ryong Ha (2004), Aras, Togan and Berkun (2007) and Han et al. (2012) works considered watercourses self-purification capacities in effluent allocation problems, using water quality simulation models and optimization techniques.In this type of problem, objective functions usually consist of minimizing costs and water quality degradation, in order to determine the removal level required for each effluent source to maintain or reach satisfactory water quality, in economically feasible and efficient conditions.
On the other hand, Nikoo, Kerachian and Niksokhan (2012), Mahjouri and Bizhani-Manzar (2013), Cho and Lee (2014), Lei et al. (2015), Ashtiani, Niksokhan and Jamshidi (2015), Nikoo, Beiglou and Mahjouri (2016), Valory, Reis and Mendonça (2016) and Fantin, Reis and Mendonça (2017) works, in addition to minimizing treatment costs, incorporated into the wastewater management problem the "Polluter Pays Principle", which, in turn, expresses the notion that as greater the responsibility for pollution, greater should be the amount paid by the polluter.Thus, these studies considered costs distribution fairness principles, either in the search for water courses self-purification capacities adequate apportionment, or in the search for adequate costs distribution, also through water quality modeling and optimization techniques.Among the different techniques used to solve optimization problems, the water resources planning and management area has seen considerable increase in the development and application of several types of evolutionary algorithms, given their potential for solving nonlinear, non-convex, multimodal and discrete problems.Among the evolutionary algorithms, Genetic Algorithms (GA) have been frequently used (NICKLOW et al., 2010;GARCIA, 2011).
Although costs are relevant in the water resources decisionmaking process, the search for problem solving involves multiple and conflicting objectives.Thus, under the appropriate technologies approach, multiobjective analysis is often used to assist the sewage treatment systems choice process, given the wide variety of technologies available and the various factors involved in the decision-making process (MOLINOS-SENANTE et al., 2015).The multiobjective analysis constitutes a set of techniques that aid the decision-making process, taking into account the influential factors complexity and society interests (SOUZA; CORDEIRO; SILVA, 2009).The different perspectives, values and preferences of those responsible and impacted by the decisions are important factors that add complexity to the decision-making process (MARTTUNEN; LIENERT; BELTON, 2017).
Technology selection for a given location alone is a highly complex problem.Different areas appropriate technologies choice (analyzed together in a multi-layer basin with different water bodies presenting different assimilation capacities) further elevates the difficulty level.Since downstream watercourse quality is affected by upstream conditions, wastewater treatment facilities cannot be considered isolated from each other (RAHM et al., 2013).
In this perspective, the present work has studied water quality simulation model, optimization technique and multiobjective analysis use to conform sewage treatment systems selection process, within watersheds.It should be noted that, among the methodologies that have been developed to assist the decision-making process regarding wastewater planning and management it is not observed the use of water quality modeling and optimization techniques together with the multiobjective analysis methods, given that, while the former two are usually associated with the technical Bringer et al.

3/14
and economic approach, the latter is related to the appropriate technology analysis approach.
The selection process addressed in this study considered Biochemical Oxygen Demand conditions established for effluents release through CONAMA Resolution 430/2011 (BRASIL, 2011), aiming to evaluate these restrictions imposition consequences in the context of the search for adequate water bodies self-purification capacities distribution among polluters.Although the methods adopted here are intended to be generic and applicable to any watershed, in this study there were employed in the Pardo watershed, watercourse located in Espírito Santo State, Brazil.

STUDY AREA
The proposed methodology was applied to the Pardo river watershed (Figure 1).Pardo is a Braço Norte Esquerdo river affluent and this is an Itapemirim river tributary.Itapemirim is the most important watercourse in the Espírito Santo State southern portion.Pardo River main spring located in Ibatiba municipality (ES), at an altitude of approximately 1,244 meters, and its mouth is localized at an about 400 meters altitude.In its watershed are located parts of Ibatiba, Irupi, Iúna, Muniz Freire and Lajinha municipalities -the latter located in Minas Gerais State.Pardo river watershed comprises about 611 square kilometers drainage area.Pardo river length corresponds to approximately 57.9 kilometer.
Pardo watershed water system is composed by Pardo river and three tributaries of greater relevance (São José river, Pardinho river and Ribeirão Perdição creek).As may be observed in Figure 2   Wastewater treatment systems selection inside watersheds by using multiobjective analysis 4/14

MATERIAL AND METHODS
The methodology used in this study presents as main objective to subsidize the decision-making process associated to domestic sewage treatment systems selection, in the context of watersheds.In order to do so, it involves (1) minimum sewage treatment efficiencies required for each effluent disposal point determination stage; and (2) possible treatment technologies, from a preference structure, hierarchical stage.

Minimum sewage treatment efficiency determination
Aiming at reducing wastewater treatment costs, the self-purification concept was incorporated into the wastewater treatment systems choice process.Thus, a water quality model and an optimization technique were applied in order to determine the minimum sewage treatment efficiencies that, within the watershed, would allow respect to the established water receiving bodies environmental standards.

Water quality simulation model
The water quality simulation was performed by the computational model developed by Valory, Reis and Mendonça (2016), in the MATLAB software environment, adapted to the Pardo river watershed.This model encompasses the physical mixing interactions and biological reactions that characterize water body natural self-purification process, reproducing QUAL-UFMG model, originally presented and discussed by Von Sperling (2007), mathematical formulations and conceptual and computational structures,.
The raw effluents kinetic constants, hydrodynamic information and organic loads considered in the study (Table 1) were obtained from Calmon et al. (2016), who proposed a methodology to support surface watercourses water intended uses classification processes, with application in the same watershed considered in this work.Calmon et al. (2016) also established flow velocity and depth adjustment equations, based on Thomann and Mueller (1987).The K D coefficient was evaluated by the authors according to the water body hydraulic characteristics (depth and flow), considering the formulation proposed by EPA (1985) and Thomann and Mueller (1987).The kinetic atmospheric reaeration constant (K2) was estimated from the original expression defined by O'Connor and Dobbins in 1958, considering the relations between flow and velocity and between flow and depth.
For BOD concentrations simulation, only the deoxygenation phenomena associated with organic matter oxidation and the external fuzzy loads contributions were considered.For DO concentrations simulation, atmospheric reaeration and deoxygenation produced by organic matter oxidation were considered.
According to Von Sperling (2014a), the organic matter sedimentation disregard is in favor of results safety, since the BOD model does not consider the BOD decrease coming from it.In addition, photosynthesis and respiration processes present approximately equal rates, and can be equivalent for water quality simulations that assume permanent flow.The author also indicates that, with respect to the lack of consideration associated with benthic demand, the studies for quantification of the benthic demand coefficient are laborious and expensive, demanding large quantity of field surveys and laboratory analyzes.The use of average data from the literature is difficult to generalize, due to the complexity of all factors interaction that may influence benthic demand (sludge composition and texture, organic matter percentage, age, depth, temperature, DO, among others).

Optimization model
Aiming to estimate the minimum sewage treatment efficiencies and in the search for an adequate apportionment of water bodies autodepuration capacity among polluters, an optimization model that imposes the minimization of an inequality measure between treatments was proposed.This model is based on Santoro, Reis and Mendonça (2016) work conclusions.The authors evaluated different optimization models applied to the determination of sewage treatment efficiencies in the Pardo river watershed.Some of the models tested by these authors did not present adequate performances, due to the existence of multiple sewage discharges presenting very different raw sewage loads and water bodies presenting different assimilation capacities.Thus, it was opted for an optimization model proposition that would include the different environmental conditions existing in a watershed (Equation 1).
( ) in the Equation 1: COe(i): i-th effluent disposal point considered in the watershed raw organic load, in mg/s; BODr(i): BOD concentration in the water stream immediately before the i-th effluent disposal point, in mg/L; E(i): i-th effluent treatment system Bringer et al.

5/14
BOD removal efficiency; Qr(i): watercourse flow immediately before the i-th effluent disposal point, in L/s.The proposed model presents as main purpose to consider the recipient water bodies condition immediately upstream the raw domestic sewage discharges.Thus, the ( ) COe .BODr E .Qr relationships indicate that streams with higher flow rates present greater sewage assimilation capacity, while water courses presenting higher BOD concentrations have lower assimilation capacity, given their lower quality conditions upstream the releases.Based on these considerations, the lower the watercourse assimilation capacity and the greater the organic load discharged, the greater the required sewage treatment efficiency for that wastewater discharge point.
The model purpose is to minimize the sum of said relationships differences between effluent disposal points, which, therefore, establishes an inequality measure between final disposal points.
The objective function was subjected to constraint inequalities (inequalities from 2 to 13) based on both sewage treatment efficiency limits and quality standards established by Brazilian CONAMA Resolution 357/2005(BRASIL, 2005) for Class II watercourses and on sewage discharge conditions with respect to BOD (CONAMA Resolution 430/2011).Thus, the constraints that conformed the optimization models established three different final effluents disposal analysis conditions: -Effluent Discharge Condition 1: -Effluent Discharge Condition 2: -Effluent Discharge Condition 3: It is important to highlight that Inequality 13, besides allowing watercourse self-purification capacity exploration, aims to guarantee non-negative estimated sewage treatment efficiencies values.

Optimization tecnique
For the proposed optimization model application it was utilized the Genetic Algorithm (GA) technique.Originally established by John Holland, in 1975, at the University of Michigan, and popularized by one of his students, David Goldberg, the GAs belong to the class of probabilistic search and optimization methods that attempt to direct the search to space regions where it is likely that solutions are (TANOMARU, 1995).They were inspired by the species evolution mechanisms, including population genetics processes, survival and individuals adaptation.
Lacerda and Carvalho (1999) point out that the first step of GA is the generation of an "initial chromosomes population", which is formed by a random set of chromosomes that represents problem possible solutions.During the "evolutionary process", this population is evaluated and each chromosome receives a note, called "aptitude", which reflects the quality of the solution it represents.In general, the fittest chromosomes are selected and the least discarded (Darwinism).The selected "individuals" can undergo changes in their fundamental characteristics through "crossover" and "mutation" operators, generating offspring for the next generation.This process is repeated until a satisfactory solution is found.
Thus, the basic idea of GAs is to treat possible solutions to the problem as "individuals" of a "population", which will "evolve" at each iteration or "generation" (POZO, 2005), been strongly dependent on the genetic operators performance (GREFENSTETTE, 1986;GEN;CHENG, 1996apud LEIRAS, 2010).
The AG parameters operators and values used in this work (Table 2) reproduced those used by Valory, Reis and Mendonça (2016).When estimating sewage treatment efficiencies for hypothetical effluent releases in the upper portion of the Santa Maria da Vitória river watershed (Espírito Santo, Brazil), these authors performed tests with the main GA parameters and their respective more commonly used values in water quality management problems.In this work stage the functions available in the Toolbox "Optimization" of MATLAB software were used.ELECTRE III makes peer-to-peer alternatives comparisons in order to accept, reject or, more often, evaluate the credibility of the statement "alternative 'a' is at least as good as alternative 'b'", taking into account the indifference and preference thresholds defined for each criterion, the degree or importance coefficients (weights) assigned to each criterion and the possible relative comparability difficulties of two alternatives when the 'a' alternative is significantly better than the 'b' alternative in certain criteria, but much worse than the 'b' alternative in at least one of the criteria.
Thus, for each criterion, two indices are calculated.One expresses the extent to which the alternatives performances ' a ' and ' b ' are in agreement with the statement 'a' is higher than 'b' (concordance index), while the other expresses exactly the opposite -to what extent the performances of the alternatives oppose this statement (discordance index).The concordance indices associated with the criteria analysis weights give rise to a global concordance matrix, which expresses the extent to which all criteria performances are in agreement with the statement "'a' is greater than 'b'".Subsequently, a credibility matrix is established, which consists of the global concordance index weakened by the discordance indices, comprehensively expressing to what extent "'a' is greater than 'b'".A classification algorithm is then used to obtain the final alternatives ordering.
Subsequent topics detail the steps in the method execution, considering the sewage treatment system selection applied to the Pardo river watershed problem.

Alternatives performances
Considering that ELECTRE III makes peer-to-peer sewage treatment alternatives comparisons by means of systems performances analysis in relation to each criterion, quantitative and qualitative comparisons between the sewage treatment systems data provided by Von Sperling (2014b) were made.
The (g (a)) alternatives performance quantitatively compared by the aforementioned author were determined by using the linear values normalization technique, as performed by Hunt (2013) when establishing a multicriteria decision support model applied to small municipalities wastewater treatment systems selection.For Cr01 to Cr03 criteria the normalization equation took Equation 14 form.For Cr04 to Cr09 criteria values normalization was performed according to Equation 15.

Weights, preference, indifference and veto thresholds definition
ELECTRE III is a technique that uses anticipated preferences articulation.Hence, its application depends on the definition by the decision maker (s) of the degrees or importance coefficients attributed to the different criteria (weights) and α and β coefficients values, necessary for the preference, indifference and veto thresholds calculation.The three thresholds can be defined and calculated, as following: -The indifference threshold ( ( ) . q q q g a α β = + ) corresponds to the greatest difference between the performances of two alternatives compatible with a situation of indifference; -The preference threshold ( ( ) ) corresponds to the smallest performance difference of two alternatives from which the best performing alternative is strictly preferred; -The veto threshold ( ( ) ) is the smallest performance difference of two alternatives from which it cannot be accepted that the worst of the two alternatives (with respect to a given criterion) can be considered as good as the best alternative, even if its performance in all other criteria are better.
As it was not possible local preference structures identification (fundamental to the decision-making process), the values attributed by Cordeiro Netto, Souza and Lopes Júnior (2001) were adapted when analyzing alternatives for anaerobic reactors effluents post-treatment by using ELECTRE III as multiobjective analysis technique.It should be noted that many of the criteria adopted by them were more recently used in sewage treatment systems selection processes employing multiobjective analysis techniques (KARIMI et al., 2011;KALBAR et al., 2012;HUNT, 2013;CORNELLI (2014); MOLINOS-SENANTE et al., 2015).Thus, for Cr01 to Cr09 criteria, the alpha and beta values required to calculate the indifference and preference thresholds were, respectively, 0.05 and 0.00; and 0.15 and 0.00.For Cr10 to Cr17 criteria, the referred values were, respectively, 0.00 and 1.00; and 0.00 and 3.00.Considering that in developing countries costs are important aspects in sewage treatment systems selection processes, the veto threshold was considered only for criteria Cr05 and Cr06, so that, in a pairwise alternatives comparison, an expensive implantation and operation treatment system could not be considered as good as a much more economical.In this case, alpha and beta assumed values, respectively, 0.99 and 0.00.

Concordance and discordance indices calculation
Thresholds determination makes it possible to calculate concordance and discordance indices.Let (a,b) be a pair of alternative sewage treatment systems.The concordance index C j (a,b), which aims to measure whether "alternative 'a' is at least as good as the alternative 'b' " in criterion g j , is calculated by Equation 16.
The discordance index is calculated according to Equation 17.

Global concordance index calculation
In order to express to what extent the performances of all the criteria are in agreement with the assertion "a is greater than b", the global concordance index C(a,b) (Equation 18) are  = , whatever the importance relative to the criterion, j w .Thus, the credibility index can be defined according to Equation 19.

Ranking algorithm
Based on the alternatives credibility indexes, the final ELECTRE III stage consists in the application of a classification algorithm that, through two complete pre-orders results intersection, provides the final sewage treatment systems classification considered in this work.
These two pre-orders are constructed in different ways.The first one is obtained in a decreasing way, starting with the best alternative and ending with the worst attribution -called descending distillation.The second is obtained in ascending order, starting with the worst classified alternative and ending with the assignment of the best -ascending distillation.Dias, Figueira and Roy (2006) present in detail the procedure performed by ELECTRE III-IV software to establish these pre-orders, also presenting a method application example, to facilitate classification algorithm understanding.

RESULTS AND DISCUSSION
The methodology application results for the Pardo river watershed will be presented in three stages.Firstly, the study area receiving raw sewage releases watercourses assimilation capacities will be discussed.Next, the sewage treatment efficiencies obtained through combined use of the water quality model with the optimization technique estimated for each locality will be presented.Finally, the sewage treatment systems selected for the Pardo river watershed will be presented and discussed, named from this point onwards as WWTP 01 (Ibatiba), WWTP 02 (Irupi), WWTP 03 (Iúna), WWTP 04 (Santíssima Trindade) and WWTP 05 (Nossa Senhora das Graças).

Pardo river watershed water quality modeling considering raw effluents final disposal
Using the water quality model developed in the MATLAB computational environment, Pardo river watershed three main watercourses DO and BOD concentration profiles were obtained considering final raw sewage disposal (Figures 3, 4 and 5).Table 4 shows, for each locality, the organic load discharged in   Bringer et al. 9/14 the watercourse (CO e ), the water flow immediately upstream the sewage discharge (Q r ), and the BOD concentrations in the water stream immediately upstream the discharge point (BOD r ) and at the discharge point (BOD r_e ).
From generated concentration profiles analysis and based on CONAMA Resolution 357/2005, even for the estimated Pardo river watershed raw sewage inflow, DO concentrations remained above the 5 mg/L limit, established for Brazilian Class 2 watercourses.However, with regard to BOD, only Santíssima Trindade and Nossa Senhora das Graças communities sewage discharges did confer on the receiving water body (Ribeirão Perdição creek) quality characteristics, DO and BOD, satisfying to Class 2 limit, function of the smaller organic load and the conditions assumed for watercourses sections located upstream the discharge points .
While Irupi wastewater discharges caused a BOD concentration increase slightly higher than 11.0 mg/L (11.2 mg/L), at Pardinho River kilometer 5, the sewage from Ibatiba produced an approximately 16.0 mg/L increase, at Rio Pardo kilometer 16.3.These localities are the main Pardo river watershed pollution sources.
Due to the lack of domestic sewage inflow, the São José stream affluence, at Rio Pardo Km 23.3, produced considerable dilution effect of the sewage produced by Ibatiba population.Similarly, the Pardinho river affluence, at Km 30.8, and the Ribeirão Perdição creek affluence, at kilometer 37, caused small BOD concentrations reductions and DO concentrations increases in the Pardo River.The Pardo River BOD peak, located between kilometers 30 and 40 (Figure 3) resulted from Iúna sewage discharge, responsible for an approximately 3.0 mg/L Biochemical Oxygen Demand increase.

Minimum required BOD removal efficiencies
Water quality modeling combined with optimization technique provided the necessary treatment efficiencies, for each locality, for the different sewage disposal conditions evaluated (Table 5).
As the main Pardo river watershed pollution source, Ibatiba sewage must be submitted to a treatment process that may guarantee a minimum 90% BOD removal efficiency, regardless of the final disposal condition evaluated.When considering the BOD discharge restrictions established by CONAMA Resolution 430/2011, Irupi sewage should be submitted to the same treatment level required for Ibatiba.As Irupi discharges the third largest amount of organic load in the watershed in a Pardinho River section that presents a lower dilution capacity, for the third sewage discharge condition (condition that considers the use of self-purification capacity for sewage assimilation) it is necessary an 81.1% organic matter removal efficiency for the raw sewage generated by its inhabitants.
As for Irupi, discharge conditions 1 and 2, for Iúna raw sewage there were required the same treatment efficiencies required for Ibatiba's.However, when considering the prospect of Pardo River's self-purification capacity use, the simulations indicated the need to remove only 27.9% of the organic matter contained in the raw sewage.It is important to note that secondary sewage treatment systems present as minimal BOD removal efficiency approximately 30%, according to treatment systems operating characteristics presented by Von Sperling (2014b).
Considering that Santíssima Trindade and Nossa Senhora das Graças communities raw effluents discharges do not result in Ribeirão Perdição creek violation of Class 2 DO and BOD limits, the simulations results always pointed out the smallest possible BOD removal values.When the minimum sewage treatment efficiency was limited to 60% (discharge condition 1), the optimization model indicated the need to install 60% efficiency sewage treatment plants in both locations.When assumed 120 mg / L maximum BOD (discharge condition 2), 70% treatment efficiencies were required (minimum efficiency required for 120 mg / L BOD, considering the adopted raw sewage characteristics).For the third discharge condition, for which the minimum sewage treatment efficiency could be 1%, the model indicated the need to remove only 3.9% and 2.3% of the BOD from the sewage generated by Santíssima Trindade and Nossa Senhora das Graças inhabitants.Thus, the results associated with discharge conditions 1 and 2 demonstrate that the restrictions imposed by CONAMA Resolution 430/2011, which ignore water bodies self-purification capacities use possibility (60% minimum BOD removal efficiencies or maximum 120 mg / L BOD), would lead to a significant increase in effluent treatment efficiencies, leading to the selection of more robust sewage treatment systems and, possibly, inadequate distribution of the generally limited financial resources allocated to wastewater treatment plants implementation and operation.
It can be stated that, in general, the water bodies effluents assimilation capacities were adequately distributed among the polluters.The water flow and quality conditions immediately prior to receiving sewage had a direct influence on the estimated treatment efficiencies, as sought by the objective function established Wastewater treatment systems selection inside watersheds by using multiobjective analysis 10/14 in Equation 1.This perspective can be observed from the analysis of the results presented for Ibatiba and Iúna, when considered the sewage discharge third condition.These locations, although presenting very close organic loads (Table 4), required very distinct treatment efficiencies (90.0 and 27.9%, respectively), since Pardo river assimilation capacity at the Iúna sewage discharge point was considerably higher than that of the same watercourse when receiving Ibatiba discharge.It should be noted that in Santoro, Reis and Mendonça (2016) work -whose optimization models did not consider the condition of the recipient bodies immediately upstream raw domestic sewage inputs -the estimated treatment efficiency for Iúna reached 94%.
Although the results obtained for discharge condition 3 indicate the need of less than 30% BOD removal efficiencies for Iúna, Santíssima Trindade and Nossa Senhora das Graças sewage, sewage treatment systems were also proposed for these locations, as presented in the subsequent topic.

Pardo river watershed sewage treatment systems selection
After obtaining the minimum BOD removal efficiencies required for each locality, those treatment alternatives presenting lower efficiencies as well as those alternatives that involved soil effluent disposal were excluded from the initial set of 30 (thirty).Thus, Table 6 presents the alternatives that have passed to the multiobjective analysis process.
It is important to note that, although the alternatives listed in Table 6 were considered, possible local technical conditions (such as terrain slope, hydrogeological characteristics and available area) could have restricted the consideration of some alternatives analyzed in this study.Nevertheless, Table 7 presents the synthesis of the sewage treatment systems selected for the five Pardo river watershed effluent disposal points, after ELECTRE III application.
From Table 7 analysis, the treatment systems selected for Ibatiba (WWTP 01) were identical (A29 -Submersed Aerated Biofilter with Nitrification) for the three final effluents disposal conditions evaluated, since the efficiency removal required for the Ibatiba remained constant (E1 = 90%).The results provided by ELECTRE III, for Ibatiba, are available in Table 8 and Figures 6 and 7.
The final sewage treatment systems ordering (Figure 7) was made by the ELECTRE III software from the results of the up and down distillations (Figure 6), calculated on the basis of credibility indices values (Table 8).As alternative A29 was ranked better than alternative A30 in the descending distillation and since both treatment systems were also classified in the ascending distillation, A29 took the first place in the ordering made by ELECTRE III, while A30 took the second place.Similarly, as in one of the distillations the alternative A30 assumed better classification than the alternatives A25 and A26, and the same classification in the other pre-order, the sewage treatment system A30 was considered superior to the treatment systems A25 and A26.On the other hand, treatment systems A25 and A26 were considered indifferent because they belonged to the same equivalence class in the two pre-orders and, therefore, assumed the third place in the final ordering.It is important to highlight that among the four alternatives evaluated for Ibatiba, the Submerged Aerated Biofilter with Nitrification resulted as the most advantageous treatment system from the economic point of view.Thus, the veto threshold incorporation (in this work considered only in the criteria related to costs) had a considerable influence on the results obtained, since it made  that the alternatives most benefited by its effect assumed at least good placement in the final ordering.
When required 90% minimum efficiencies, the sewage treatment systems selected for Irupi (WWTP 02) and Iúna (WWTP 03) were the same ones selected for Ibatiba, since the same weights and criteria were considered in the multiobjective analysis.
When required treatment efficiencies for Irupi and Iúna, 81.1% and 27.9% (WWTP's 02 and 03, respectively), the Wetland type system was recommended for both locations, presenting as advantages practically zero energy requirements, high BOD removal efficiency, simple construction, operation and maintenance, reduced implantation and operation costs, good resistance to load variations and no sludge generation.It is important to highlight that the selected alternative presented one of the lowest implantation and operation costs (criteria influenced by the veto effect) when compared to the other treatment systems that entered the multiobjective analysis, besides presenting one of the best coliforms removal performances (criterion presenting greater weight in the selection process).Nevertheless, in general, the aforementioned treatment system has good performance in most of the analysis criteria considered in the selection process, which favors its good placement in the ELECTRE III hierarchy.
It is important to note that the choice of the same system (Wetland type) for Irupi and Iúna, with considerably different treatment efficiencies (27.9 and 81.1%), indicates the need to incorporate, in the multiobjective analysis, a criterion regarding the BOD removal efficiency.Such incorporation would probably make it possible for the multiobjective analysis results to be closer to those provided by water quality simulation model and optimization technique application.
Regardless of the required BOD removal efficiencies, two treatment systems took the first place in the final order indicated by ELECTRE III for Santíssima Trindade (WWTP 04) and Nossa Senhora das Graças (WWTP 05) communities: Facultative Pond (A03) and Wetland (A13).As Facultative Ponds advantages, it is worth noting the satisfactory BOD removal efficiency, simple construction, operation and maintenance, practically zero energy requirements, sludge removal required only after periods of more than 20 years and reduced implementation and operating costs.It is important to note that for the WWTP 03 (Iúna) choice the greatest relevance multiobjective analysis criterion consisted of the coliform removal efficiency, for Trindade (WWTP 04) and Nossa Senhora das Graças (WWTP 05) the greater criteria weight in the multiobjective analysis were related with implementation and operation costs, as well as operational simplicity.Thus, even starting from the same treatment alternatives set (for WWTP 03, WWTP 04 and WWTP 05, treatment efficiencies lower than 30% were required), the results for the third sewage discharge condition were influenced by differences in the criteria weights.

CONCLUSIONS
From water quality model combined with an optimization technique and multiobjective analysis application results, the main conclusions of this work can be summarized, as follows: -   Wastewater treatment systems selection inside watersheds by using multiobjective analysis 12/14 considered as constraints on the optimization problem (discharge conditions 1 and 2), efficiencies ranged from 60 to 90%.In scenarios where these restrictions were not incorporated into the optimization model efficiencies ranged from approximately 3% to 90%.Thus, ignoring water bodies self-purification capacities may can lead to significant increase in required effluent treatment efficiencies, overestimation of treatment plants and, possibly, inadequate distribution of financial resources, usually limited, for WWTP's implementation and operation; -The application of ELECTRE III allowed to indicate the treatment systems considered most appropriate, within the analysis context, for each Pardo River watershed raw sewage discharge point.Multiobjective analysis results demonstrated that the method is sensitive to each criterion weights variations and to the adoption of veto thresholds.Whenever required removal efficiencies above 90%, ELECTRE III indicated the Submerged Aerated Biofilter system (with Nitrification) as an alternative treatment.The Wetland type treatment system was indicated for different locations and different discharge conditions.In some cases, both the Facultative Ponds and the Wetlands were selected as the first treatment option; -The sewage treatment systems indicated for the Pardo river watershed five localities are associated with the analysis context considered, resulting from the preference structure adopted.In this way, the establishment of a diverse preference structure or the incorporation of local technical restrictions, fundamental to the decision-making processes, could restrict some treatment technologies implementation.
-The use of water quality simulation model and optimization technique, prior to the application of multiobjective analysis, has proved to be relevant in the sewage treatment systems selection context, since it includes, in the appropriate technology search, watercourses self-purification capacity analysis.This approach allowed to indicate sewage treatment systems, at the river watershed level, considering the "Polluter Pays Principle" and aiming at minimizing required treatment efficiencies (by taking advantage of watercourses self-purification capacities).However, the results indicated the need to adopt, in the multiobjective analysis, the criterion associated to BOD removal efficiency (in this work, the BOD removal efficiency constituted a criterion applied to the sewage treatment systems pre-selection), to avoid selecting the same types of sewage treatment systems for locations where the required minimum efficiencies are much different.
, the Pardo River receives domestic sewage from Ibatiba and Iúna municipalities.The Pardinho river, in turn, has its quality altered by Irupi municipality sanitary sewage.Santíssima Trindade and Nossa Senhora das Graças villages discharge wastewater in the Ribeirão Perdição creek.In none of the effluent disposal points there exists Sewage Treatment Plants (SWT).

.
Wastewater treatment systems selection inside watersheds by using multiobjective analysis 8/14 calculated.The global concordance index consists of the partial concordance indexes C j (b,a) sum in each of the criteria, weighted by each criterion weights, w j .Taking into account both the global concordance index and the discordance indices for each criterion, the credibility indexes ( , ) a b σ are established, which comprehensively express to what extent the "a" alternative is considered better than "b" alternative.Credibility is merely the global concordance index weakened by the discordance indices.In discordance indices absence, This credibility value is reduced in one or more discordant criteria presence, when ( ) j D a,b > C(a,b).In accordance with the veto effect (

Figure 3 .
Figure 3. Pardo river DO and BOD concentration profiles.
Final disposal of raw sewage, regardless of the scenario considered, would not allow compliance with the water quality standards associated with BOD established for Brazilian CONAMA class 2 watercourses in long Pardo and Pardinho rivers stretches; -The use of the water quality model combined with the optimization technique allowed the appropriation of minimum BOD removal efficiencies for the different Pardo river watershed effluent disposal points, seeking the most equitable distribution of watercourses assimilation capacities among polluters.In simulations where the minimum BOD removal levels or maximum BOD concentrations indicated by Brazilian environmental legislation wereTable 8. Credibility matrix for treatment alternatives considered in the multiobjective analysis for WWTP 01 (Ibatiba) choice.alternative B credibility index lies in the line A with column B intersection cell.

Figure 6 .
Figure 6.Result of distillations carried out by ELECTRE III.

Figure 7 .
Figure 7. Final order carried out by ELECTRE III.

Table 1 .
Pardo river watershed raw effluents kinetic constants, hydrodynamic variables and organic loads.

Table 2 .
(ACHILLAs et al., 2013)s used in GA application.Considering the necessary effluent disposal sources treatment efficiencies in the watershed and aiming to incorporate into the sewage treatment systems selection process the Appropriate Technology corollaries, the final indication of the sewage treatment systems came from ELECTRE III -(ACHILLAS et al., 2013)application -using the "ELECTRE III-IV" software developed by the University of Paris -Dauphine Laboratory of Analysis and Modeling of Decision Support Systems (LAMSADE).ELECTRE is one of the European decision analysis school main and most used methods families (SOUZA; CORDEIRO; SILVA, 2009).Achillas et al. (2013)indicate that, in addition to being widely used in the field of waste management in general, ELECTRE III has been pointed out as one of the most appropriate methodologies for sewage treatment systems selection.It presents the advantage of being able to incorporate into the decision-making process a large number of quantitative and qualitative evaluation criteria(ACHILLAs et al., 2013), being therefore suitable for sewage treatment systems selection processes.

Table 3
summarizes the weights assigned to each criterion considered in this work.

Table 3 .
Weights assigned to the criteria considered in the multiobjective analysis.

Table 4 .
Waste load and characteristics of receiving watercourses.

Table 6 .
Sewage treatment alternatives considered in the multiobjective analysis.

Table 7 .
Sewage treatment alternatives selected by ELECTRE III.