Wastewater Treatment Plant Assessment by Quantifying the Carbon and Water Footprint

: In the context of e ﬃ cient and sustainable management of the elements of the urban water cycle as an aim of the Water Framework Directive (WFD), the evaluation of indicators such as the water footprint (WF) and the carbon footprint (CF) in a wastewater treatment plant (WWTP) provides a quantiﬁcation of the environmental impact, both negative and positive, which implies its exploitation. In this study, in addition to WF and CF quantiﬁcation, a joint evaluation of both indicators was conducted. Consumption is indicated by the blue water footprint (WF Blue ) and emissions by CF. Both are related to the operational grey water footprint ( ∆ WF G,mef ) in two ratios, WFR and CFR. In this way, the water consumed and gases emitted are measured according to the reduction range of the pollutant load of the discharge. The results for four WWTPs show operational scenarios for better management in accordance with the WFD.


Introduction
Directive 2000/60/EC, known as the Water Framework Directive (WFD), establishes a legal basis for the protection of water bodies in all European Union member states. To this purpose, WFD challenges water resource managers to work on the evaluation and sustainable management of the water cycle, in which wastewater treatment plants (WWTPs) have a crucial role [1]. Sustainability assessment of WWTPs through indicators such as those collected in the environmental footprint family [2] is vital to identify the impact of WWTPs on the environment and to improve the process towards cleaner, and more sustainable production [3]. Footprint assessment methods clearly and easily identify, describe and quantify the complexity of water-related impact associated with wastewater management [4].
The water footprint (WF) and the carbon footprint (CF) are the most commonly used indicators to assess the sustainability of WWTPs [5]. The concept of WF responds to the need for an indicator of freshwater consumption [6]. WF measures the volume of water consumed directly and indirectly to manufacture a product or provide a service [7]. It is a multidimensional indicator with three sub-indicators: blue water footprint (WF Blue ) refers to the consumption of surface and groundwater; green water footprint (WF Green ) measures the consumption of rainwater stored in the soil, and grey water footprint (WF Grey ) refers to pollution and defines the volume of fresh water needed to assimilate the pollutant load based on the required water quality standards [7]. The Water Footprint Network (WFN) has developed the Water Footprint Assessment (WFA) tool, which is the most widely used methodology for the evaluation of WF [7]. The introduction of the ISO 14046 standard [8] resulted in other methods to evaluate WF based on life cycle assessment (LCA) to quantify potential environmental impacts related to water. In the recently published Available Water Remaining (AWARE) method [9], the water consumption impact indicator helps to assess the water scarcity footprint according to In contrast, Morera et al.'s [17] study did not consider the upstream concentration of the discharge point (c nat ). Gu et al. [18] quantified the reduction in the greywater footprint from the influent (c r ) and effluent concentrations of the WWTP. Gómez-Llanos et al. [19] considered the operational grey water footprint in addition to the effluent and influent concentration from the WWTP and the maximum concentration allowed.
The objective of this paper is to develop and apply a joint methodology for the evaluation of WF and CF in four WWTPs located in La Vera (Cáceres, Spain). Previous studies in this regard are scarce, particularly those applied to WWTPs. The WF evaluation is based on the WFA, particularly the operational grey water footprint [19], and CF is estimated as an indicator of global warming in LCA. In this way, we quantify the merged relationships between water quality and the resources necessary for their treatment. The efficiency of this process is evaluated through two new ratios in which the WF and CF are related: the water footprint ratio (WFR) and the carbon footprint ratio (CFR). These ratios, defined in this study, in addition to assessing efficiency, allow comparisons between different plants subjected to various circumstances.
Therefore, this study contributes to the development of the application of footprint assessment methodologies, in this case WF and CF, in a WWTP. The proposed methodology, based on the two ratios, enables quantifying the connection between consumption or impact as defined by the WF Blue and CF, and the operational management margin necessary to comply with the regulations, considering the working concentrations of the WWTP and the values of the receiving environment, both of which are established by the regulations as well as those existing in the receiving water body.

Proposed Methodology
As noted above, the assessment of the WF is considered the WFA methodology [7], and CF is evaluated as an indicator of the impact of global warming established in LCA [11,12].
According to both methodologies, it is necessary to define the physical and temporal boundaries of the studied systems and functional unit. As shown in Figure 1, the edges of the supply chain studied include different WWTP processes: the water line, made up of pre-treatment, primary, secondary or biological treatment, and possible tertiary treatment; the sludge line, where primary and secondary sludges are treated, including a possible biogas line. In each step, water is consumed or used, and GHG emissions are produced, either directly or indirectly. Figure 1 shows these consumption and emission sources; for instance, the indirect water consumption associated with energy consumption, which is included in WF Blue or direct GHG emissions produced in biological treatment, which is considered in CF. The functional unit adopted is the volume of treated water during one month of operation.
In the four WWTPs studied, the treatment of wastewater is a conventional biological process (activated sludge), and tertiary treatment is not necessary. The sludge treatment is thickening and involves physical centrifugation, without stabilisation based on aerobic or anaerobic digestion; for instance, the gas line does not exist. The final sludges are applied outside the plants in agricultural parcels. In the four WWTPs studied, the treatment of wastewater is a conventional biological process (activated sludge), and tertiary treatment is not necessary. The sludge treatment is thickening and involves physical centrifugation, without stabilisation based on aerobic or anaerobic digestion; for instance, the gas line does not exist. The final sludges are applied outside the plants in agricultural parcels.

Blue Water Footprint as a Consumption Indicator
WFBlue, according to the WFA Manual [7], is defined as the volume of freshwater consumed directly and indirectly in all processes of the supply chain of a product or service; that is, the volume of water that is incorporated into the final product or by-products, the water evaporated during the process and the lost flow [7].
In the case of a WWTP, WFBlue for the period of time considered responds to Equation (1): where, for the period of time considered (t), WFBlue (m 3 /t) is the total blue water footprint, WFElect (m 3 /t) the blue water footprint due to electricity consumption, WFChem (m 3 /t) the blue water footprint due to the consumption of products derived from the risk of purification, VolumeSludge (m 3 /t) is the volume of freshwater incorporated into the sludge produced by the operation of the treatment plant and VolumeEvap (m 3 /t) is the volume of water evaporated in the biological reactors, decanters and storm tanks of the secondary WWTP process. WFElect was evaluated with the conversion factor of ECOINVENT 3.0 used by Morera et al. [17], Gómez-Llanos et al. [19] and Martínez-Alcalá et al. [26]. A lack of data on chemicals consumed prevented assessment of WFChem. As for the WWTP sludge considered in VolumeSludge, the WFBlue calculated in this by-product is the volume of water embedded in said product; that is, the volume of water present in the extracted sludge when it is dehydrated. In our WWTPs, the dehydration phase was carried out by the centrifugation process, and the percentage of dry matter in the final removed   [7], is defined as the volume of freshwater consumed directly and indirectly in all processes of the supply chain of a product or service; that is, the volume of water that is incorporated into the final product or by-products, the water evaporated during the process and the lost flow [7].
In the case of a WWTP, WF Blue for the period of time considered responds to Equation (1): where, for the period of time considered (t), WF Blue (m 3 /t) is the total blue water footprint, WF Elect (m 3 /t) the blue water footprint due to electricity consumption, WF Chem (m 3 /t) the blue water footprint due to the consumption of products derived from the risk of purification, Volume Sludge (m 3 /t) is the volume of freshwater incorporated into the sludge produced by the operation of the treatment plant and Volume Evap (m 3 /t) is the volume of water evaporated in the biological reactors, decanters and storm tanks of the secondary WWTP process. WF Elect was evaluated with the conversion factor of ECOINVENT 3.0 used by Morera et al. [17], Gómez-Llanos et al. [19] and Martínez-Alcalá et al. [26]. A lack of data on chemicals consumed prevented assessment of WF Chem . As for the WWTP sludge considered in Volume Sludge , the WF Blue calculated in this by-product is the volume of water embedded in said product; that is, the volume of water present in the extracted sludge when it is dehydrated. In our WWTPs, the dehydration phase was carried out by the centrifugation process, and the percentage of dry matter in the final removed sludge is, on average, 23%. That is, the average percentage of water content of the sludge is 77%. Volume Evap was estimated based on the Penman-Monteith method, applied to quantify natural evaporation from open water surfaces [27].

Grey Water Footprint as an Efficiency Indicator
In this paper, WF Grey is evaluated based on the WFA and a previous study [19]. Considering WF Grey as dilution factor to express the efficiency of the treatment, the operational grey water footprint, ∆WF G,mef , defined by Gómez-Llanos et al. [19], enables the evaluation of the WTTP operation, according to Equation (2), depending on the value of this operational index: where Q r is the volume of treated water in a period of time (volume/time), c r and c ef are the concentrations of the respective pollutant measured in the WWTP influent and effluent, respectively (mass/volume), c max is the concentration of the respective pollutant allowed by the standard (mass/volume) and c nat is the concentration of the respective pollutant measured in the receiving body of water (mass/volume). As can be seen from Figure 2, for the current scenario, where the concentration of effluent (c ef ) is lower than the concentration allowed by the standard (c max ), ∆WF G,mef results in a positive value, which indicates a higher limit in WWTP management. In a scenario with equal flow and standard concentration values, an increase in ∆WF G,mef benefits the receiving environment, since the level of required self-purification of biodegradable pollutant discharge is reduced. In the scenario with discharge over the concentration by the standard, ∆WF G,mef is a negative value, which would indicate non-compliance with the legal standard. WWTP managers are interested in rigorous compliance with the discharge requirement, especially if it reduces WF Blue . In this scenario, ∆WF G,mef is null or has positive value close to zero.
Water 2020, 12, x FOR PEER REVIEW 5 of 16 sludge is, on average, 23%. That is, the average percentage of water content of the sludge is 77%. VolumeEvap was estimated based on the Penman-Monteith method, applied to quantify natural evaporation from open water surfaces [27].

Grey Water Footprint as an Efficiency Indicator
In this paper, WFGrey is evaluated based on the WFA and a previous study [19]. Considering WFGrey as dilution factor to express the efficiency of the treatment, the operational grey water footprint, ∆WFG,mef, defined by Gómez-Llanos et al. [19], enables the evaluation of the WTTP operation, according to Equation (2), depending on the value of this operational index: where Q is the volume of treated water in a period of time (volume/time), cr and cef are the concentrations of the respective pollutant measured in the WWTP influent and effluent, respectively (mass/volume), cmax is the concentration of the respective pollutant allowed by the standard (mass/volume) and cnat is the concentration of the respective pollutant measured in the receiving body of water (mass/volume). As can be seen from Figure 2, for the current scenario, where the concentration of effluent (cef) is lower than the concentration allowed by the standard (cmax), ∆WFG,mef results in a positive value, which indicates a higher limit in WWTP management. In a scenario with equal flow and standard concentration values, an increase in ∆WFG,mef benefits the receiving environment, since the level of required self-purification of biodegradable pollutant discharge is reduced. In the scenario with discharge over the concentration by the standard, ∆WFG,mef is a negative value, which would indicate non-compliance with the legal standard. WWTP managers are interested in rigorous compliance with the discharge requirement, especially if it reduces WFBlue. In this scenario, ∆WFG,mef is null or has positive value close to zero.

Carbon Footprint Assessment in WWTP
CF is the indicator of the global warming impact category in LCA. CF is measured as CO2 equivalent mass per unit of time. To quantify CF, it is necessary to establish the sources of the GHG emissions and the corresponding conversion factors, which were obtained from the characterisation models developed by other authors. The total CF is the sum of the CF of each emission source, as indicated in Equation (3): where, for the period of time considered (t), CF (kg CO2/t) is the total carbon footprint, Emi (unit/t) is the value of emission source i and Fi (kg CO2/unit) is the conversion factor for emission source i. In a WWTP, direct emissions are produced especially in secondary treatment. In this biological treatment, the GHGs emitted are methane (CH4) and nitrous oxide (N2O). The estimate of GHGs is quantified based on the proposed methodology by the Intergovernmental Panel on Climate Change [28]. Indirect sources of GHG emissions are the electrical consumption of the plant and the treatment of the sludge obtained. The required data to estimate CF and the references for the considered emission factors are in Table 1.

Carbon Footprint Assessment in WWTP
CF is the indicator of the global warming impact category in LCA. CF is measured as CO 2 equivalent mass per unit of time. To quantify CF, it is necessary to establish the sources of the GHG emissions and the corresponding conversion factors, which were obtained from the characterisation models developed by other authors. The total CF is the sum of the CF of each emission source, as indicated in Equation (3): where, for the period of time considered (t), CF (kg CO 2 /t) is the total carbon footprint, Em i (unit/t) is the value of emission source i and F i (kg CO 2 /unit) is the conversion factor for emission source i. In a WWTP, direct emissions are produced especially in secondary treatment. In this biological treatment, the GHGs emitted are methane (CH 4 ) and nitrous oxide (N 2 O). The estimate of GHGs is quantified based on the proposed methodology by the Intergovernmental Panel on Climate Change [28]. Indirect sources of GHG emissions are the electrical consumption of the plant and the treatment of the sludge obtained. The required data to estimate CF and the references for the considered emission factors are in Table 1.

Joint Assessment of WF and CF in WWTP
For the joint assessment of the WF and CF, two new relations were defined: the WF ratio (WFR) and the CF ratio (CFR).
The relationship of the blue water footprint with the operational grey water footprint [19] is defined as the water footprint ratio, expressed in Equation (4): where ∆WF G,mef is operational grey water footprint (m 3 /time) and WF Blue is blue water footprint (m 3 /time). WFR is dimensionless and establishes the relationship of the volume of water consumed directly or indirectly with the purification function in a WWTP. With this ratio, the efficiency with which the plant reduces the wastewater pollutant load is quantified up to the WWTP output parameters based on the water consumption made for it. That is the ratio of the volume of water required to carry out the purification function of the plant.
If the operational grey water footprint indicator is related to the CF, the carbon footprint ratio (CFR) is defined by Equation (5). CFR indicates the purification function in a WWTP concerning the environmental impact in terms of GHG emissions.
where ∆WF G,mef is the operational grey water footprint (m 3 /time) and CF is the carbon footprint (kg CO 2 -eq/time), and K is equal to the unit kg CO 2 -eq per cubic meter and is used so that CFR is dimensionless. Figure 3 shows a scheme with the variation in the base indicators (∆WF G,mef , WF Blue and CF) of Equations (4) and (5) for the determined period of operation in a WWTP (OR t ) and an estimate of the life period. The range of variation of ∆WF G,mef in the life period of the WWTP with proper operation could cover from ∆WFG,mef with a value of zero, with effluent concertation equal to the maximum legal concentration (c ef ≈ c max , ∆WF G,mef ≈ 0), to ∆WF G,mef , with a value of effluent concentration equal to natural concentration (c ef ≈ c nat , ∆WF G,mef > 0) for a specific contaminant in the aqueous receiving medium. The variation of the determined life period studied is actually restricted at the interval marked with OR t in Figure 3. Considering WF Blue or CF variation during the lifetime of the WWTP, the range could be defined by unknown upper and lower limits, WF Blue,min and WF Blue,Max , Water 2020, 12, 3204 7 of 16 and similarly CF Min and CF Max . This interval would be measured during a specific period of time, OR t in Figure 3. Where cef is the concentration of effluent of a determined pollutant, cmax is the concentration allowed by the standard and cnat is the concentration of respective pollutant measured in receiving body of water. Possible operational limits of WWTP in time t and for its life period, where GE, good environmental performance limit; GM, good manager performance limit; OS, operational scenario to reduce consumption; PE, poor environmental performance limit.
Throughout the operating life of the WWTP, ∆WFG,mef value a could decrease to zero, with cef equal to cmax, and value b would be the maximum of ∆WFG,mef when cef is similar to cnat. For the CF and WFBlue values, c would reach an unknown minimum value marked by the fixed consumption of the WWTP, and finally, d would be the unknown maximum value of these indicators.
In this way, WFR should be defined for the different scenarios as: In general, Equation (6) is verified: CFR would be expressed similarly to Equations (6) and (7). Therefore, for the particular time under study, as the values obtained are approximately those defined in Equation (6), the operation of the WWTP would be estimated from the environmental and managerial point of view for both the best and worst scenario. Figure 4 shows a summary of the ∆WFG,mef, WFBlue and CF values corresponding to the descriptions of different scenarios identified. Where c ef is the concentration of effluent of a determined pollutant, c max is the concentration allowed by the standard and cnat is the concentration of respective pollutant measured in receiving body of water. Possible operational limits of WWTP in time t and for its life period, where GE, good environmental performance limit; GM, good manager performance limit; OS, operational scenario to reduce consumption; PE, poor environmental performance limit.
Throughout the operating life of the WWTP, ∆WF G,mef value a could decrease to zero, with c ef equal to c max , and value b would be the maximum of ∆WF G,mef when c ef is similar to c nat . For the CF and WF Blue values, c would reach an unknown minimum value marked by the fixed consumption of the WWTP, and finally, d would be the unknown maximum value of these indicators.
In this way, WFR should be defined for the different scenarios as: In general, Equation (6) is verified: CFR would be expressed similarly to Equations (6) and (7). Therefore, for the particular time under study, as the values obtained are approximately those defined in Equation (6), the operation of the WWTP would be estimated from the environmental and managerial point of view for both the best and worst scenario. Figure 4 shows a summary of the ∆WF G,mef , WF Blue and CF values corresponding to the descriptions of different scenarios identified.
Water 2020, 12, x FOR PEER REVIEW 8 of 16 Figure 4. ∆WFG,mef, WFBlue and CF values corresponding to descriptions of different scenarios identified. GE, good environmental performance limit; GM, good manager performance limit; OS, operational scenario to reduce consumption; PE, poor environmental performance limit.
From the management point of view, the operating scenario is defined as the maximum values of both indicators, WFBlue or CF and ∆WFG,mef, where the necessary measures must be studied and established to reduce the negative impact on the environment indicated by WFBlue and CF with the least damage to the operating range. Operating with low ∆WFG,mef values, even close to zero, would be enough, since the discharge would continue to comply with the regulations. Operating with smaller WFBlue and CF values would be more appropriate. This scenario is defined as the good operating limit for managers.
From an environmental point of view, operating with high ∆WFG,mef values near the limit established by cnat is more beneficial, since WWTP discharge would require less self-purification volume. Operating with smaller WFBlue and CF values would be more appropriate. This limit is the good environmental performance limit (GE).
The poor environmental performance limit is defined with maximum WFBlue and CF values and minimum ∆WFG,mef values. This scenario indicates a negative impact on the environment characterised by WFBlue and CF, and a higher volume of self-purification required by being an adjusted value of ∆WFG,mef.

WWTPs Data Inventory
The WF and CF were evaluated in four WWTPs located in the La Vera region, northeast of Extremadura (SW Spain), which serve several municipalities ( Table 2). The WWTPs are (1) Aldeanueva de la Vera and Cuacos de Yuste, (2) Arroyomolinos, Pasarón de la Vera and Tejeda de Tiétar, (3) Villanueva de la Vera and Valverde de la Vera and (4) Torremenga.
The four WWTPs were designed for different treatment capacities, specifically different equivalent population (EP) values, as shown in Table 2. The inventory data were extracted from the quarterly reports by their management company. In these reports, daily and monthly data of each WWTP are collected. The temporal boundary was the operational period from January 2014 to December 2016. The assessment was conducted considering the monthly data of the four WWTPs. By way of a summary, Table 2 shows the monthly average of the studied period of treated flow and electricity consumption, and Table 3 shows the average monthly concentrations at the entry and exit of each WWTP. From the management point of view, the operating scenario is defined as the maximum values of both indicators, WF Blue or CF and ∆WF G,mef , where the necessary measures must be studied and established to reduce the negative impact on the environment indicated by WF Blue and CF with the least damage to the operating range. Operating with low ∆WF G,mef values, even close to zero, would be enough, since the discharge would continue to comply with the regulations. Operating with smaller WF Blue and CF values would be more appropriate. This scenario is defined as the good operating limit for managers.
From an environmental point of view, operating with high ∆WF G,mef values near the limit established by c nat is more beneficial, since WWTP discharge would require less self-purification volume. Operating with smaller WF Blue and CF values would be more appropriate. This limit is the good environmental performance limit (GE).
The poor environmental performance limit is defined with maximum WF Blue and CF values and minimum ∆WF G,mef values. This scenario indicates a negative impact on the environment characterised by WF Blue and CF, and a higher volume of self-purification required by being an adjusted value of ∆WF G,mef .

WWTPs Data Inventory
The WF and CF were evaluated in four WWTPs located in the La Vera region, northeast of Extremadura (SW Spain), which serve several municipalities ( Table 2). The WWTPs are (1) Aldeanueva de la Vera and Cuacos de Yuste, (2) Arroyomolinos, Pasarón de la Vera and Tejeda de Tiétar, (3) Villanueva de la Vera and Valverde de la Vera and (4) Torremenga. The four WWTPs were designed for different treatment capacities, specifically different equivalent population (EP) values, as shown in Table 2. The inventory data were extracted from the quarterly reports by their management company. In these reports, daily and monthly data of each WWTP are collected. The temporal boundary was the operational period from January 2014 to December 2016. The assessment was conducted considering the monthly data of the four WWTPs. By way of a summary, Table 2 shows the monthly average of the studied period of treated flow and electricity consumption, and Table 3 shows the average monthly concentrations at the entry and exit of each WWTP. The maximum concentration (c max ) values for the main pollutants analysed (biochemical oxygen demand at five days (BOD 5 ), total nitrogen (TN) and total phosphorus (TP)) were obtained from Directive 91/271/EEC and its transposition to Spanish legislation in Royal Decree-Law 11/1995. The Royal Decree-Law establishes the standards applicable to urban wastewater treatment and the limit concentrations for BOD 5 , total nitrogen and total phosphorus, shown in Table 3. Table 4 shows the average values of WF Blue per cubic metre of treated water for the four WWTPs in the studied period (2014 to 2016). These values represent the freshwater consumed in the operation phase of each WWTP to purify the discharge. As can be seen in Table 4, the four plants had their highest consumption in 2015, as indicated by WF Blue . Torremenga WWTP had its highest WF Blue value in 2015 (1.37 m 3 ·month −1 /m 3 on average, as shown in Table 4). On average, Torremenga is the plant with the most water consumption despite being the smallest (933 PE).

WF Assessment
The year with the lowest water consumption of all WWTPs was 2016, except Cuacos, where the lowest water consumption was in 2014. The average value for the three studied years for the WWTPs is 0.78 m 3 per metre of treated water.
The purification function of plants is defined by the indicator ∆WF G,mef . In general, for the four WWTPs studied and three pollutant parameters analysed, the study scenario was that in which those parameters are above the values established by existing legislation. Figure 5 shows average values of ∆WF G,mef per cubic metre of treated water to three pollutant parameters analysed. In general, TN and TP concentrations show lower values than BOD 5 in the three years for the four WWTPs. For example, BOD 5 ∆WFG,mef at Cuacos was 9.41 m 3 per cubic metre of water treated in 2014, versus to TN and TP with ∆WF G,mef values of 2.42 and 6.62 m 3 per cubic metre cubic of water treated, respectively. Therefore, TN and TP concentrations have a more limited reduction margin with regard to the values established by the standard, similar to what was observed by Chen et al. [5]. The purification function of plants is defined by the indicator ∆WFG,mef. In general, for the four WWTPs studied and three pollutant parameters analysed, the study scenario was that in which those parameters are above the values established by existing legislation. Figure 5 shows average values of ∆WFG,mef per cubic metre of treated water to three pollutant parameters analysed. In general, TN and TP concentrations show lower values than BOD5 in the three years for the four WWTPs. For example, BOD5 ∆WFG,mef at Cuacos was 9.41 m 3 per cubic metre of water treated in 2014, versus to TN and TP with ∆WFG,mef values of 2.42 and 6.62 m 3 per cubic metre cubic of water treated, respectively. Therefore, TN and TP concentrations have a more limited reduction margin with regard to the values established by the standard, similar to what was observed by Chen et al. [5].
The highest ∆WFG,mef values for BOD5 were in 2015. Figure 5 also shows that, in general,  The highest ∆WF G,mef values for BOD 5 were in 2015. Figure 5 also shows that, in general, ∆WF G,mef increased significantly in 2015 in the Cuacos, Tejeda and Villanueva WWTPs. On the other hand, ∆WF G,mef decreased in 2014 in the Tejeda, Villanueva and Torremenga WWTPs, which indicates operations with a lower operational margin. Table 5 shows the results obtained for CF per cubic metre of treated water in the four WWTPs. As can be appreciated, the year with the highest GHG emissions was 2015 for all WWTPs; for example, 0.52 kg CO 2 /m 3 for Cuacos. Tejeda had the lowest CF values in the three years under study. The average value for the three studied years for the WWTPs is 0.28 kg CO 2 /m 3 , similar to that obtained by Lorenzo-Toja et al. [24], who analysed 22 treatment plants in Spain by LCA, getting a value of 0.278 kg CO 2 /m 3 for plants with biological treatment of active sludge.  Table 5 shows the results obtained for CF per cubic metre of treated water in the four WWTPs. As can be appreciated, the year with the highest GHG emissions was 2015 for all WWTPs; for example, 0.52 kg CO2/m 3 for Cuacos. Tejeda had the lowest CF values in the three years under study. The average value for the three studied years for the WWTPs is 0.28 kg CO2/m 3 , similar to that obtained by Lorenzo-Toja et al. [24], who analysed 22 treatment plants in Spain by LCA, getting a value of 0.278 kg CO2 /m 3 for plants with biological treatment of active sludge.  As can be seen in Figure 6, the performance of WFR and CFR is similar for TP and TP, respectively. This is a consequence of the similar evolution of annual values of CF and WF Blue . As can be seen in Tables 4 and 5, both indicators showed similar growth, with a steeper slope for WF Blue , and the CF value is approximately constant. CF and WF Blue have a strong relationship with electricity consumption. Similar to the results reported by Gómez-Llanos et al. [19] and Chen et al. [5], the main factor in WF Blue assessment is electricity consumption. However, the primary emission source in CF assessment of the four WWTPs studied is CH 4 , with electricity consumption as the second factor, similar to Hospido et al. [22]. Figure 6a,c show the relationships among values of TN ∆WF G,mef and WF Blue and CF, respectively. The scenario of 2014 was from an environmental point of view, not only by lower water consumption as indicated by WF Blue , but also by the highest value of TN ∆WF G,mef . In contrast, 2016 could be considered a beneficial scenario from the business point of view, with a WF Blue value close to the minimum and rigorous compliance with the legal standards. The operational scenario (OS) indicates the reduction in water consumption, as indicated for 2015. Figure 6b,d indicate the relationships among values of TP ∆WF G,mef and WF Blue and CF, respectively. For this contaminant, 2016 could be considered a beneficial scenario from the business point of view. However, a good environmental scenario is not identified, since the highest TP ∆WF G,mef value was produced with the highest WF Blue and CF value in 2015.

WFR and CFR Analysis
Next, the WFR and CFR values representing the different operating limit scenarios are assessed based on the analysis carried out in Figure 6 for each WWTP. Table 6 shows the rates for each year and WWTP. The highlighted cells with a coloured background show some of the situations of better or worse operation of WWTPs from the environmental or managerial point of view. Table 6. Annual average values of WFR and CFR for four WWTPs, 2014-2016. Good environmental performance limit (GE) is indicated by green cells, good manager performance limit (GM) by yellow cells, operational scenario to reduce consumption (OS) by orange cells, and poor environmental performance limit (PE) by grey cells. The values represented in Figure 6 have been written in red. As can be seen in Table 6, from an environmental point of view, the values corresponding to proper operation give the highest WFR and CFR values for TN and TP in 2014-2016. There are situations where this scenario has not been identified: WFR and CFR of TP at Cuacos and WFR and CFR of TP and TN at Torremenga. In both cases, the highest value cannot be clearly categorised as an excellent environmental limit, since WF Blue and CF were not near the minimum in 2014. From the management point of view, the cases where the better operation of the WWTP is observed the tendency is to generate intermediate values of WFR and CFR for TN and TP.

Cuacos
Regarding BOD 5 , the highest value of WFR and CFR are linked to the operational scenario to reduce consumption (OS), and both indicators, WF Blue and CF and ∆WF G,mef , are the maximum values for the period studied. This situation was investigated to evaluate the possibility of decreasing the environmental impact in correlation with the allowed diminution in the operating margin. In the four WWTPs studied, the OS was in 2015.
Poor environmental limit was obtained with the lowest rates of WFR and CFR for TP and TP in the period analysed (Table 6). In cases where the better operation of the WWTP is observed according to the manager, the tendency is the generation of intermediate WFR and CFR values for TN and TP and a minimum BOD 5 value (Table 6).

Conclusions
A joint evaluation of water and carbon footprint in a municipal wastewater treatment service of highlights the following aspects: • The operation and exploitation of WWTPs cause environmental impact, which is quantified by WF Blue , at 0.78 m 3 per cubic metre of water treated on average, and CF, at 0.28 kg CO 2 /m 3 for the four WWTPs. However, WWTPs have a positive impact on the environment; the reduction in the pollution load in the discharge as indicated by ∆WF G,mef is better than the strict reduction established by the regulations in the four WWTPs in the period studied.

•
In the four WWTPs, the values of ∆WF G,mef revealed a high operational margin in the elimination of BOD 5 . Instead, there is a lower operational range of the total nitrogen and total phosphorous pollutant parameters. For example, BOD 5 ∆WF G,mef at Cuacos was 9.41 m 3 per cubic metre of water treated in 2014, compared to ∆WF G,mef values of 2.42 and 6.62 m 3 per cubic metre cubic of water treated for TN and TP, respectively.

•
The proposed WFR and CFR ratios indicate optimal operational scenarios from an environmental or managerial point of view. The highest values of WFR and CFR for TN and TP belong to the best scenario for the environment. In cases where better operation of the WWTP is observed, according to the manager, the tendency is to generate intermediate values of WFR and CFR for TN and TP.

•
Poor environmental scenarios, taking into account the criteria of managers, show the lowest rates of WFR and CFR for TN and TP in the period of time analysed. • Regarding BOD 5 , the highest WFR and CFR values are found in a scenario where the possibility of decreasing the environmental impact in correlation with the allowed reduction in the operation range was investigated. In the four WWTPs studied, this scenario happened in 2015.
The application limitations of the proposed methodology will be derived mainly from the compilation of the necessary data, as well as from the estimation of the factors associated with the blue water footprint evaluation linked to energy consumption or the calculation of the carbon footprint.
This study, together with previous works, establishes a new framework for multidisciplinary decision-making that combines water use, consumption and environmental impacts. Regarding the evaluation of footprints in a WWTP, the future development of the water footprint, mainly the gray water footprint, plays a relevant role in the study of the implication of a WWTP in the receiving environment and its capacity for self-purification.

Conflicts of Interest:
The authors declare no conflict of interest. BOD 5 Biochemical