Next Article in Journal
Effects of Commonly Occurring Metal Ions on Hydroxyapatite Crystallization for Phosphorus Recovery from Wastewater
Previous Article in Journal
Estimating Hydraulic Conductivity Fields in Composite Fan Delta Using Vertical Electrical Sounding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing Water Pollution Potential in Life Cycle Impact Assessment Based on Bacterial Growth and Water Quality Models

1
State Key Laboratory of Urban Water Resource and Environment, School of Environment, Harbin Institute of Technology, Harbin 150090, China
2
College of Resources and Environment, Northeast Agricultural University, Harbin 150030, China
*
Authors to whom correspondence should be addressed.
Water 2018, 10(11), 1621; https://doi.org/10.3390/w10111621
Submission received: 6 October 2018 / Revised: 31 October 2018 / Accepted: 4 November 2018 / Published: 9 November 2018
(This article belongs to the Section Wastewater Treatment and Reuse)

Abstract

:
For the life cycle assessment (LCA) of wastewater management, eutrophication is considered the most relevant factor. However, eutrophication is not the only pathway through which wastewater influences the environment, and merely characterizing eutrophication potential is not sufficient for the LCA framework to reflect the influence of wastewater. This study defines the Bacterial Depletion of Oxygen (BDO)—a new impact category that represents the oxygen depleting potential caused by the growth of microorganisms—and characterization models and characterization factors are developed for the application of BDO. Water quality models (both one- and two-dimensional) are incorporated into the BDO characterization models so that the LCA framework includes some spatially differentiated factors, and can be used to estimate the direct impact of wastewater on receiving environment (IBDO value). Based on three case studies, this study demonstrates how the BDO category can be applied for the evaluation of wastewater management. Results show that increases in the downstream distance and self-purification coefficients reduce the IBDO value, whereas the increase in water velocity raises the IBDO value. Future integration of the BDO category with water quality models must link the dilution effect of water bodies, the environmental carrying capacity of receiving water, and the distribution of water pollutants in eutrophication and bacterial oxygen depletion.

1. Introduction

Increasingly sophisticated discharge requirements and the requirement to avoid shifting the burden from one environmental impact to another have necessitated the application of scientific tools to the decision-making process concerning wastewater management. Life cycle assessment (LCA) has been widely used for the environmental evaluation of wastewater management [1,2,3,4,5,6,7,8,9,10], and eutrophication is considered the most relevant impact category [7]. Accordingly, most of the methodological LCA progress in regards to wastewater management has focused on the development of a eutrophication impact category [11,12,13,14,15]. These methods enable the examination of the potential impacts of eutrophication in freshwater and marine ecosystems associated with the development of fate and transport models [16].
Eutrophication is not the only type of water pollution. The depletion of oxygen due to the growth of microorganisms is another important pathway through which wastewater influences the environment, especially for wastewater that contains a large amount of organics pollutants [17]. Specifically, when wastewater flows into water bodies without treatment, organic pollutants such as carbohydrates, protein, and organic acid can encourage the growth of aerobic bacteria, consuming vast amounts of dissolved oxygen (DO). If the oxygen supply from air and photosynthesis is lower than the oxygen consumption, the DO concentration decreases. Once the concentration approaches zero, the growth of anaerobic bacteria is stimulated and the water conditions deteriorate, killing aquatic species and causing serious water pollution associated with the release of unpleasant smell and depreciated value. Deterioration of water quality due to organic enrichment is a major environmental issue in industrialization and urbanization. In China, water pollution caused by microorganism growth has been found in many contaminated rivers [18,19], and reducing this kind of pollution has become the objective of many environmental policies introduced by the Chinese government [20,21].
In general, LCA results have a high level of abstraction without referring to any specific situations [22]. Generalizable LCA results are justifiable for environmental impacts with large-scale characteristics such as global warming [23]. For emissions that directly impact the regional environment, however, it is crucial for the LCA framework to incorporate spatially differentiated factors to obtain more representative results. This is especially important when conducting LCA for wastewater-related issues because the emissions of water pollutants directly affect regional water bodies [24,25]. In previous LCA studies, some regional characterization factors or models for eutrophication have been developed, including spatially explicit factors at the lake- or stream-level [11,12,13], spatially differentiated factors for different emission sources, and a fate model for phosphorus emissions on the continent and country scales [14,15]. However, most studies did not consider the transport and transformation of water pollutants, including organic matter, nitrogen, and phosphorus, in the receiving environment. When released into water bodies, water pollutants move and undergo processes that often involve simultaneous biological, chemical, and physical reactions [17]. Due to this variety of mechanisms, the pollutant concentration varies in different parts of water bodies, with different potential impacts. Such situations have rarely been included in current calculation LCA frameworks. This highlights the need to incorporate the fate of water pollutants in LCA frameworks for obtaining estimation results that represent the direct impact of water pollutants on different parts of water bodies. Water quality models (WQM) are effective tools for simulating and predicting the fate and behavior of pollutants in water environments [26]. WQM include some geological, chemical, or physical factors related to the water environment, and can be used to identify the levels in different parts of a water environment.
The goal of this study wastwo2-fold. First, a new impact category is defined for LCA system, the Bacterial Depletion of Oxygen (BDO), to characterize the oxygen depleting potential that is encouraged by the growth of microorganisms due to the release of water pollutants. Following the general framework for characterization, characterization models for the BDO category are developed. Second, WQM are incorporated into the characterization models with a BDO category. The LCA system includes the transport and transformation of water pollutants in the receiving environment, and can be used to estimate the direct impact of organic pollutants on the different parts of the receiving environment.

2. Methods and Materials

2.1. Bacterial Depletion of Oxygen

This study defines an independent impact category named the Bacterial Depletion of Oxygen (BDO), which is in parallel with other impact categories in LCA, such as eutrophication, global warming, or acidification. The BDO category is defined to represent the oxygen depleting potential caused by the growth of microorganisms due to the release of water pollutants. Based on the general characterization LCA framework [27], the characterization model for BDO category is described as:
I B D O ( S ) = p B D O p × m p
where S is the area affected by the substance p, p refers to water pollutants that can stimulate the growth of bacteria, I B D O ( S ) refers to the characterization results of the BDO impact category, m p is the emitted amount of substance p, and B D O p is the characterization factors, representing the potential contribution of substance p to bacterial growth.
Characterization factors ( B D O p ) are calculated by aggregating the potential contribution of water pollutants to the growth of microorganisms. In this study, organic matter (measured as the chemical oxygen demand, or COD), nitrogen (measured as the total nitrogen, TN) and phosphorus (measured as the total phosphorus, TP) are considered the major water pollutants (substance p) involved in the depletion of oxygen in receiving water bodies. The formula describing the characterization factors is:
B D O p = v p M p v r e f M r e f
where v p and v r e f represent the potential contribution to the bacterial growth by one mole of the substance p and the reference substance, respectively; and M p and M r e f refer to the molar mass of the substance p and the reference substance, respectively. There are three types of reference substance p: oxygen (O2), nitrate (NO3), and phosphate (PO43−). When characterizing the potential growth of bacteria in a LCA, choosing O2 can provide an intuitive measure of oxygen depletion (kg O2 eq.), whereas selecting NO3 (kg NO3 eq.) or phosphate PO43− (kg PO43− eq.) can contribute to a direct comparison between the newly defined BDO category and the currently available eutrophication category [27]. Of note, the identification of substance references would minimally influence the characterization of environmental impact in LCA. Therefore, the potential contribution of reference substance ( v r e f ) is usually assigned a value of 1 [28].
The determination of v p ( v T N , v T P , and v C O D ) is based on the identification of the characterization state. Given the scientific background of the eutrophication category, the average chemical composition of aquatic organisms (C106H263O110N16P) is assumed to be representative of average biomass [27]. This study determines the growth of bacteria with the characterization state of the BDO category, and v p is calculated on the basis of the stoichiometry and the chemical composition of the bacterial biomass (CnHaObNc) [29]:
C n H a O b N c + ( 2 n + 0.5 a 1.5 c b 2 ) O 2 n C O 2 + c N H 3 + ( a 3 c 2 ) H 2 O
where one mole of the bacterial biomass ( C n H a O b N c ) requires c moles of N and (2n + 0.5a − 1.5cb)/2 moles of chemical oxygen demand (COD) and v p is calculated as:
v T N = 1 c v C O D = 1 ( 2 n + 0.5 a 1.5 c b 2 ) v T P = 0

2.2. Incorporation of Water Quality Models into the LCA Framework

To include the fate of water pollutants in the LCA framework, this study incorporates water quality models (WQM) into the characterization model of the BDO category. Two kinds of WQM were selected in this study: the one-dimensional (1D) model and the two-dimensional (2D) model. In this study, the BDO characterization model coupled with 1D is called the one-dimensional coupling characterization model (1DCCM), and the BDO characterization model coupled with 2D is called the two-dimensional coupling characterization model (2DCCM).

2.2.1. Descriptions of 1D and 2D

The development of 1D is based primarily on the Streeter-Phelps oxygen-sag equation and the steady-state assumption that constant parameters are required, such as the flow rate of waste [30]. 1D is often used for small-sized rivers [31]. In the 1D model, the concentration of water pollutants is assumed to vary along the longitudinal orientation, while all waste constituents are assumed to be completely mixed vertically and laterally within any cross section [32]. The equations describing 1D are as follows:
c t + u c x = E 2 c x 2 k c
z t + 1 B Q x = q t
Q t + x ( Q 2 A ) + g A ( z x + Q | Q | A ) = 0
where Equation (4) is the differential equation for 1D, Equation (5) is the continuity equation of flow, Equation (6) is the momentum equation for the longitudinal orientation, c is the pollutant concentration, u is the velocity of the receiving water, E is the longitudinal dispersion coefficient, t is the time required for the pollutant to move over one cross-section, x is the longitudinal distance that pollutant travels through the cross-section, k is the self-purification coefficient of the pollutant, z is the water level, Q is the rate of flow, B is the width of the receiving water, A is the cross-sectional area, and g is the gravitational acceleration, q t is the flow rate arising from the side-river.
Based on Equations (4)–(6), the concentration of pollutants along the longitudinal orientation can be described as:
C ( x ) = m A u exp ( k x u )
where C(x) is the pollutant concentration at the cross-section, which is x meters away from the emission point, and m is the emitted amount of the pollutant.
The development of 2D is based upon the finite-difference solution of the vertically averaged equations of fluid motion [33,34]. 2D describes time-varying lateral and longitudinal distributions of water pollutants in a water body, and is often used for large water bodies including large-sized rivers, reservoirs, lakes, and estuaries [35]. The use of 2D depends on the assumption that water quality varies along the longitudinal and lateral axes of a water-body, and all waste constituents are assumed to be completely mixed vertically [34]. The equations describing 2D are as follows:
c t + u x c x + u y c y = M x 2 c x 2 + M y 2 c y 2 k c
ρ z t + ρ H u x x + ρ H u y y = 0
ρ H u x t + x ( ρ H u x u x ) + y ( ρ H u y u x ) = ρ g H z x + x ( H γ e f f u x x ) + y ( H γ e f f u x y ) τ h x
ρ H u y t + x ( ρ H u x u y ) + y ( ρ H u y u y ) = ρ g H z y + x ( H γ e f f u y x ) + y ( H γ e f f u y y ) τ h y
where Equation (8) is the differential equation of 2D; Equation (9) is the continuity equation of flow; Equations (10) and (11) are the momentum equation for the longitudinal orientation and lateral orientation, respectively; c is the pollutant concentration; k is the self-purification coefficient of the pollutant; t is the time required for the pollutant to transfer; ux and uy are water velocities along the longitudinal and lateral orientation, respectively; M x and M y are longitudinal and horizontal dispersion coefficients, respectively; x is the longitudinal distance from emission point to downstream impact point; y is the lateral distance from emission point to downstream impact point; H is the average depth of receiving water; g is the gravitational acceleration; ρ is the water density; z is the water level; τ h x is the momentum coefficient along x orientation and τ h y is the momentum coefficient along y orientation; γ e f f is the coefficient of virtual viscosity.
Based on Equations (8)–(11), the concentration of pollutants at the downstream impact point (x,y) can be described as:
C ( x , y ) = C p × Q p × exp ( k x 86400   u x ) × 1 H × π M y x u x × [ exp ( u x y 2 4 M y x ) + exp ( u x ( 2 B y ) 2 4 M y x ) ]
where C (x,y) is the pollutant concentration at the downstream impact point (x,y), Q p is the flow rate of the emission point, and C p is the pollutants concentration at the emission point.

2.2.2. Development of 1DCCM and 2DCCM

This study defines a time-varying mass parameter for the BDO characterization model, and Equation (1) can be described as follows:
I B D O ( x ) = p [ B D O p × m p , t ( x ) ]
I B D O ( x , y ) = p [ B D O p × m p , t ( x , y ) ]
where m p , t ( x ) and m p , t ( x , y ) are the masses of the pollutant that travel through a cross-section (x) or an impact point (x,y) at a given time t, respectively; m p , t ( x ) is determined by the flow rate of receiving water and the C ( x ) of 1D; and m p , t ( x , y ) is determined by the flow rate of receiving water and the C ( x , y ) of 2D.
The 1DCCM was developed by combing Equations (7) and (13), and can be described as follows:
I B D O ( x ) = p [ B D O p × m p , t ( x ) × exp ( k p x u ) ]
where I B D O ( x ) is the characterization result of the BDO category, representing the oxygen depleting potential of the emitted water pollutants at the cross-section that is x meters away from the emission point; COD and TN are considered the major water pollutants (substance p); and t is the time required for water pollutants to reach the cross-section, determined by x and u.
The 2DCCM was developed by combing Equations (12) and (14), and can be described as follows:
I B D O ( x , y ) = { p [ B D O p × m p , t ( x , y ) × exp ( k p x u ) ] } × [ 1 ( H π M y x u x ) ] × [ exp ( u x y 2 4 M y x ) + exp ( u x ( 2 B y ) 2 4 M y x ) ]
where I B D O ( x , y ) is the characterization result of the BDO category, which represents the oxygen depleting potential of the emitted water pollutants at the impact point (x,y).

2.3. Identification of Atrazine and Its Metabolites

We first demonstrates how to calculate the characterization factors for the BDO category. The growth of bacteria (mixed cultures and pure cultures) are determined as the characterization state. The potential contribution by one mole of substance p ( v p ) to the bacterial growth w calculated by incorporating the bacterial empirical formula into Equation (3). The selected bacterial empirical formula included C5H7O2N, C7H12O4N, C9H15O5N, C9H16O5N, C4.9H9.4O2.9N, C4.7H7.7O2.1N, C4.9H9O3N, C5H8.8O3.2N, C4.1H6.8O2.2N, C5.1H8.5O2.5N and C5.3H9.2O2.5N (for mixed culture), as well as the C5H8O2N, C5H8.33O0.81N, C4H8O2N, C4.17H7.42O1.38N, C4.54H7.91O1.95N, C4.17H7.21O1.79N, C4.16H8O1.25N, and C3.85H6.69O1.78N (for pure cultures). All empirical equations were determined from the technical book. The COD and TN were identified for the substance p; therefore, he v C O D and v T N could be determined. TP was excluded because its minimal contribution identified by Equation (3). The characterization factors of the BDO category were calculated by introducing v C O D and v T N into Equation (2).
We then investigated the application of 1DCCM and 2DCCM. 1DCCM was applied to a small-sized river, the Ashihe river located in China. One river reach of the Ashihe river was selected in this study: the reach from Xiquanyan to Maan mountain. The impact of a single parameter on the characterization results of BDO, I B D O ( x ) , was explored, and the parameters included the downstream cross-section that is x meters away from the emission point, the self-purification coefficient of COD (kCOD), and the water velocity (u). We explored the combined impact of these parameters on I B D O ( x ) . 2DCCM was applied to a large-sized river, the Songhua river, located in China. The Harbin mainstay of the Songhua river was selected, which is 16.25 km long, 450 m wide, and 2.3 m deep. We explored the change in I B D O ( x ) along the longitudinal orientation at different lateral sections. Six lateral sections of the Harbin mainstay were selected: y = 0 m, 50 m, 100 m, 150 m, 200 m, and 250 m. Both the rivers were assumed to be affected by an emission point. The flow rate of wastewater was 10,000 m3/day, and the COD of raw wastewater was 259.2 mg/L and TN was 28.7 mg/L.
We finally conducted a complete LCA analysis using the CML method (the most widely used LCA method developed by Leiden University) on a full-scale wastewater treatment plant (WWTP) using the BDO category and 1DCCM. The WWTP is located near the bank of the Ashihe river, and the major treatment technology is cyclic activated sludge technique, operated at capacity of 10,000 m3/day. The components of raw wastewater entering the WWTP included a COD of 259.2 mg/L, TN of 28.7 mg/L, TP of 3.4 mg/L, and suspended solids of 187.4 mg/L. The goal of the LCA analysis was to investigate the environmental impact of different wastewater treatment levels. Three scenarios were defined for the LCA analysis based on the three wastewater treatment levels: scenario 1 (tertiary treatment), scenario 2 (intermediate treatment), and scenario 3 (basic treatment). The 10,000 m3 of wastewater emitted per day was chosen as the functional unit. The operational stage of the WWTP was evaluated, and the system boundaries included sewage treatment, electricity production, chemical manufacturing and transportation, and waste-activated sludge processing. The BDO category (kg NO3 eq.) was used in parallel with the other impact categories of the CML method, and these impact categories included acidification (A, kg SO2 eq.), freshwater aquatic ecotoxicity (FAET, kg 1,4-DCB eq.), human toxicity (HT, kg 1,4-DCB eq.), ozone depletion (OD, kg CFC-11 eq.), photochemical oxidation (PO, kg ethylene eq.), global warming (GW, kg CO2 eq.), abiotic depletion of fossil fuels (ADF, MJ), and abiotic depletion of elements (ADE, kg antimony eq.). We determined five cross-sections for the LCA analysis: section 1 (100 m away from the emission point), section 2 (20,000 m away), section 3 (40,000 m away), section 4 (60,000 m away), and section 5 (80,000 m away).

3. Results

3.1. Characterization Factors of BDO

The characterization factors for the BDO category were calculated (Table 1). The values of 1, 4.4286, and 6.7857 for B D O T N (characterization factor of TN) were obtained for both mixed cultures and pure cultures when selecting O2, NO3, and PO43− as the reference substance, respectively. All the selected bacterial empirical equations produced the same number of nitrogen atoms and the same potential contribution of substance N to the bacterial growth ( v T N = 1 ). The values of B D O C O D (characterization factor of COD) depended on the identification of the bacterial empirical formula, varying with the different environmental conditions and reference substances. For the following study, the selection of the environmental conditions or the reference substances had little influence. Accordingly, the averaged characterization factors (0.3759 of B D O C O D and 4.4286 of B D O T N ) were used with the reference substance being NO3.

3.2. Application of 1DCCM and 2DCCM

Based on 1DCCM, we explored the impact of single parameters on the characterization results of the BDO category. The selected parameters included the downstream cross-section distance (x), self-purification coefficient of COD (kCOD), and water velocity (u).
As shown in Figure 1, the characterization results (IBDO) showed decreasing tendencies with increasing kCOD and x, whereas an increasing tendency with increasing u was observed. Temperature was a main contributor to the change in the self-purifying coefficient: the higher the temperature, the higher the kCOD. In general, the metabolic activity of microorganisms in the receiving water is enhanced by the increase in water temperature. This indicated that with the increase in water temperature, the self-purifying capacity of the receiving water would be stronger (leading to the increase in kCOD), and more water pollutants could be degraded before being transported to downstream cross-sections (leading to the decrease in IBDO).
With the increase in water velocity (u), a shorter period would be required or water pollutants to be self-purified by the receiving water; therefore, more water pollutants could be left to impact the downstream cross-section (leading to the increase in IBDO). Accordingly, with the increase in downstream cross-section distance (x), the receiving water would have more time to self-purify, leading to a lower environmental impact on the downstream cross-section (a decrease in IBDO).
In reality, factors such as self-purifying coefficients (k) and water velocity (u) vary with the change in environmental conditions, and these factors would impact the receiving water’s self-purifying capacity synthetically but not separately. In our case (the Harbin mainstay of the Songhua river), all three parameters (u, kCOD, and kTN) changed during various periods. For example, from March to April, the values of these parameters increased. As described above, the increased water velocity (u) had a tendency to increase IBDO, whereas the increased self-purifying coefficients (kCOD and kTN) had a tendency to decrease IBDO.
Based on the 1DCCM, it was possible to investigate the combined impact of the three increased parameters on the IBDO. As shown in Figure 2, the value of IBDO in April (1882 kgNO3) was higher than the value in March (1626 kgNO3), indicating that the decreasing tendency driven by kCOD and kTN was overwhelmed by the increasing tendency driven by u.

3.3. Results of IBDO Obtained from 2DCCM

With 2DCCM, more detailed information about the migration and transportation of water pollutants in the receiving water could be converted into the LCA characterization results. Figure 3 describes the longitudinal (x-axis) gradients of IBDO at different horizontal (y-axis) river sections. Results of IBDO at the horizontal section of y = 0 m showed a reducing tendency ranging from 1956 kgNO3 (x = 0 m, y = 0 m) to 723 kgNO3 (x = 0 m, y = 7000 m) and stabilized at 558 ± 69 kgNO3 for the rest of the longitudinal distances. However, at the horizontal sections of y = 50 m and y = 100 m, the characterization results showed increasing tendencies, with IBDO stabilizing at 453 ± 42 kgNO3and 231 ± 9 kgNO3, respectively. For the other horizontal sections of y = 150 m, y = 200 m, and y = 250 m, characterization results were almost quantitatively negligible. One possible reason for this finding might be that few water pollutants migrated to these sections in this study case.

3.4. Comparison of Total Environmental Impact and IBDO for Different Scenarios

From scenario 3 to scenario 2, reduction of IBDO and decrease in total impact were evident for all downstream sections (Figure 4). This means that upgrading the WWTP from basic treatment to intermediate treatment would be beneficial for both the receiving water and the environment overall.
With respect to the escalation from intermediate treatment to tertiary treatment, scenario 1 showed lower IBDO results than scenario 2 at all downstream locations, but the magnitude of reduction decreased gradually from DWP-1 (95.16 kgNO3) to DWP-5 (26.56 kgNO3). We compared total impact between scenario 1 and scenario 2, with scenario 1 being lower than scenario 2 at DWP-1 but higher at DWP-5. These results indicated that when examining the environmental implications of upgrading to tertiary treatment from intermediate treatment, the upgrade appeared beneficial if the evaluated downstream location was near the emission point, whereas the upgrade was probably not favorable if the evaluated downstream was far from the emission point.

3.5. Contribution of Impact Categories on Total Environmental Impact

Using 1DCCM provided further insight into how the contribution of different impact categories for one scenario would change with different environmental conditions. Taking scenario 1 as an example (Table 2), BDO accounted for 57.37% of the total impact at DWP-1 (100 m) but the percentage decreased to 41.49% at DWP-5 (80,000 m). For the other impact categories, ADF of 6.87%, GW of 8.38%, FAET of 13.74%, and A of 7.85% were obtained at DWP-1, but ADF of 9.43%, GW of 11.51%, FAET of 18.86%, and A of 10.78% were recorded at DWP-5.

4. Discussion

Characterizing the potential bacterial depletion of oxygen is a midpoint impact assessment that explains the quantitative relationship between water pollutants, such as chemical oxygen demand (COD) and its midpoint impact. Coupled with water quality models, the developed 1DCCM and 2DCCM could better represent a more specific midpoint effect by aggregating water pollutants’ migration and transformation in the receiving water. Although the choice of aerobic or anaerobic state was insignificant in this study, in LCA practice, using anaerobic bacterial growth as the characterization state would increase the characterization result of the BDO impact category and add weight to the BDO impact category among the trade-offs across various impact categories.
The impact assessment obtained from general characterization models were generalizable, so the characterization results usually account for potential site-generic impact to the total environment. considerable work has been completed to consider regional environmental characteristics in the LCA framework [36,37,38,39,40]. By integrating water quality models into the BDO characterization model in this study, the characterization results were more site-specific in terms of the WWTP cases. More importantly, the 1DCCM and 2DCCM allow for more flexible LCA analysis, exploring one-scenario environmental implications in light of different receiving water features and comparing different scenarios’ environmental implications with a wider evaluating scope. Thus, more informative results could be obtained. Obtaining more informative results means that more information is needed. Except for the emission amount of water pollutants needed in the general eutrophication characterization models, longitudinal parameters relevant to the downstream water body were necessary for the calculation of 1DCCM, and both longitudinal and horizontal parameters were indispensable for 2DCCM. To capture the migration and transformation of water pollutants more accurately, a three-dimensional water quality model might be coupled with LCA characterization models in the future. If so, vertical parameters would be required. Collecting these parameters related to the receiving water would, to some extent, increase the complexity of performing the LCA analysis.
With the coupling models, the potential contribution of COD and TN to bacterial depletion of oxygen could be modelled and quantified. Similarly, contribution to the eutrophication potential could be addressed by the coupling models. However, the coupling models would not be able to simultaneously address the eutrophication and bacterial depletion of oxygen. To avoid repeated calculations, the amount of water pollutants that were accounted for by the eutrophication potential cannot be used again in the characterization of bacterial depletion of oxygen. Maybe a relevant scientific principle is needed to distribute the water pollutants reasonably in order for coupling models to convincingly quantify both eutrophication and BDO.
In the present work, the potential bacterial growth has been put forward as the characterization state of the evaluated downstream river sections. During the process of migration and transformation, the self-purifying mechanism could also be accounted for by the potential bacterial growth, partly due to the metabolic activity of bacteria whose growth is activated by water pollutants. Thus, bacterial growth could be treated as the force driving the removal of water pollutants before water pollutants reach the evaluated downstream water body, and could be regarded as the potential environmental impact once water pollutants have entered the evaluated downstream waterbody. There is an limitation in the application of 1DCCM or 2DCCM in full LCA analysis: there is no mechanism to handle the situation where the environmental carrying capacity is broken down so that the self-purifying mechanism fails to operate anymore before water pollutants reach the evaluated downstream waterbody. This issue highlights a future research direction for coupling models, where the potential environmental impact characterized by bacterial growth could be linked with the environmental carrying capacity of receiving water.
The self-purifying mechanism was the major basis for the development of 1DCCM and 2DCCM. Another important mechanism that was not considered here was that a greater water volume of the receiving water could help dilute the content of water pollutants, contributing to the reduction of the negative effect of water pollutants. A greater water volume was generally represented by higher water velocity, meaning that with the consideration of the effect of dilution in coupling models, a higher velocity would drive a decrease in IBDO, opposite of the tendency shown in Figure 1. As the coupling models continue to develop, it is important to synthesize the dilution effect and self-purifying effect improve the simulation of water pollutant activity in receiving water in the LCA characterization results.

5. Conclusions

In order to map the oxygen depleting potential caused by bacterial growth due to the release of water pollutants (COD and TN) and the fate of water pollutants in receiving water (transportation and transformation), we defined and developed a new impact category, the Bacterial Depletion of Oxygen (BDO) category, and its associated characterization models and factors. Characterization factors were calculated based on the chemical composition of the bacterial biomass. Coupling models (1DCCM and 2DCCM) were formed by integrating the BDO category into WQMs. By incorporating site-specific information into the evaluation process, the coupling models enabled the LCA of wastewater-related issues for actual regional situations. The impacts of parameters on modelling output (IBDO value) were investigated, and the results showed that the IBDO value was reduced by the increase in the downstream distance and self-purification coefficients but was raised by the increase in water velocity. Future development of the BDO category and coupling models need to create links with the dilution effect of the receiving environment, and include the more completed processes of oxygen depletion and water pollutant behavior.

Author Contributions

S.B. and X.Z. performed the experiments and wrote the paper. X.W. and N.R. contributed to the discussion and analyzed the experiments.

Funding

This research was funded by the National Key R&D Program (No. 2016YFC0401105), State Key Laboratory of Urban Water Resource and Environment (Harbin Institute of Technology) (No. XQ2018002), and the National Natural Science Foundation of China (Grant No. 31870110).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alfonsin, C.; Lebrero, R.; Estrada, J.M.; Muñoz, R.; Kraakman, N.B.; Feijoo, G.; Moreira, M.T. Selection of odour removal technologies in wastewater treatment plants: A guideline based on Life Cycle Assessment. J. Environ. Manag. 2015, 149, 77–84. [Google Scholar] [CrossRef] [PubMed]
  2. Harder, R.; Heimersson, S.; Svanstrom, M.; Peters, G.M. Including pathogen risk in life cycle assessment of wastewater management. 1. Estimating the burden of disease associated with pathogens. Environ. Sci. Technol. 2014, 48, 9438–9445. [Google Scholar] [CrossRef] [PubMed]
  3. Heimersson, S.; Harder, R.; Peters, G.M.; Svanstrom, M. Including pathogen risk in life cycle assessment of wastewater management. 2. Quantitative comparison of pathogen risk to other impacts on human health. Environ. Sci. Technol. 2014, 48, 9446–9453. [Google Scholar] [CrossRef] [PubMed]
  4. Pasqualino, J.C.; Meneses, M.; Abella, M.; Castells, F. LCA as a decision support tool for the environmental improvement of the operation of a municipal wastewater treatment plant. Environ. Sci. Technol. 2009, 43, 3300–3307. [Google Scholar] [CrossRef] [PubMed]
  5. Risch, E.; Gutierrez, O.; Roux, P.; Boutin, C.; Corominas, L. Life cycle assessment of urban wastewater systems: Quantifying the relative contribution of sewer systems. Water Res. 2015, 77, 35–48. [Google Scholar] [CrossRef] [PubMed]
  6. Yoshida, H.; Clavreul, J.; Scheutz, C.; Christensen, T.H. Influence of data collection schemes on the Life Cycle Assessment of a municipal wastewater treatment plant. Water Res. 2014, 56, 292–303. [Google Scholar] [CrossRef] [PubMed]
  7. Corominas, L.; Foley, J.; Guest, J.S.; Hospido, A.; Larsen, H.F.; Morera, S.; Shaw, A. Life cycle assessment applied to wastewater treatment: State of the art. Water Res. 2013, 47, 5480–5492. [Google Scholar] [CrossRef] [PubMed]
  8. Li, Y.; Hou, X.; Zhang, W.; Xiong, W.; Wang, L.; Zhang, S.; Wang, P.; Wang, C. Integration of life cycle assessment and statistical analysis to understand the influence of rainfall on WWTPs with combined sewer systems. J. Clean. Prod. 2018, 172, 2521–2530. [Google Scholar] [CrossRef]
  9. Bai, S.; Zhu, X.; Wang, X.; Ren, N. Identify stakeholders’ understandings of life cycle assessment results on wastewater related issues. Sci. Total Environ. 2018, 622–623, 869–874. [Google Scholar] [CrossRef] [PubMed]
  10. Bai, S.; Zhao, X.; Wang, D.; Zhang, X.; Ren, N. Engaging multiple weighting approaches and Conjoint Analysis to extend results acceptance of life cycle assessment in biological wastewater treatment technologies. Bioresour. Technol. 2018, 265, 349–356. [Google Scholar] [CrossRef] [PubMed]
  11. Basset-Mens, C.; Anibar, L.; Durand, P.; van der Werf, H.M.G. Spatialised fate factors for nitrate in catchments: Modelling approach and implication for LCA results. Sci. Total Environ. 2006, 367, 367–382. [Google Scholar] [CrossRef] [PubMed]
  12. Gallego, A.; Rodriguez, L.; Hospido, A.; Moreira, M.T.; Feijoo, G. Development of regional characterization factors for aquatic eutrophication. Int. J. Life Cycle Assess. 2010, 15, 32–43. [Google Scholar] [CrossRef]
  13. Azevedo, L.B.; Henderson, A.D.; van Zelm, R.; Jolliet, O.; Huijbregts, M.A.J. Assessing the importance of spatial variability versus model choices in life cycle impact assessment: The case of freshwater eutrophication in europe. Environ. Sci. Technol. 2013, 47, 13565–13570. [Google Scholar] [CrossRef] [PubMed]
  14. Helmes, R.J.K.; Huijbregts, M.A.J.; Henderson, A.D.; Jolliet, O. Spatially explicit fate factors of phosphorous emissions to freshwater at the global scale. Int. J. Life Cycle Assess. 2012, 17, 646–654. [Google Scholar] [CrossRef]
  15. Struijs, J.; Beusen, A.; de Zwart, D.; Huijbregts, M. Characterization factors for inland water eutrophication at the damage level in life cycle impact assessment. Int. J. Life Cycle Assess. 2011, 16, 59–64. [Google Scholar] [CrossRef]
  16. Morelli, B.; Hawkins, T.R.; Niblick, B.; Henderson, A.D.; Golden, H.E.; Compton, J.E.; Cooter, E.J.; Bare, J.C. Critical review of eutrophication models for life cycle assessment. Environ. Sci. Technol. 2018, 52, 9562–9578. [Google Scholar] [CrossRef] [PubMed]
  17. Gonzalez, S.O.; Almeida, C.A.; Calderon, M.; Mallea, M.A.; Gonzalez, P. Assessment of the water self-purification capacity on a river affected by organic pollution: Application of chemometrics in spatial and temporal variations. Environ. Sci. Pollut. Res. Int. 2014, 21, 10583–10593. [Google Scholar] [CrossRef] [PubMed]
  18. Meng, W.; Su, Y.; Zheng, B. Analysis of current situation of water pollution and its control strategy for Chinese river basins. J. China Inst. Water. 2004, 4. [Google Scholar] [CrossRef]
  19. MEP. Report on the State of the Environment in China; Ministry of Environmental Protection: Beijing, China, 2014.
  20. MEP. Discharge Standard of Pollutants for Municipal Wastewater Treatment Plant (GB18918-2002); China Environment Press: Beijing, China, 2002.
  21. MEP. Surface Water Quantity Standards (GB8938-2002); China Environment Press: Beijing, China, 2002.
  22. Pizzol, M.; Weidema, B.; Brandão, M.; Osset, P. Monetary valuation in life cycle assessment: A review. J. Clean. Prod. 2015, 86, 170–179. [Google Scholar] [CrossRef]
  23. Renou, S.; Thomas, J.S.; Aoustin, E.; Pons, M.N. Influence of impact assessment methods in wastewater treatment LCA. J. Clean. Prod. 2008, 16, 1098–1105. [Google Scholar] [CrossRef]
  24. Zang, Y.W.; Li, Y.; Wang, C.; Zhang, W.L.; Xiong, W. Towards more accurate life cycle assessment of biological wastewater treatment plants: A review. J. Clean. Prod. 2015, 107, 676–692. [Google Scholar] [CrossRef]
  25. Zhao, X.; Yang, J.; Zhang, X.; Wang, L.; Ma, F. Evaluation of bioaugmentation using multiple life cycle assessment approaches: A case study of constructed wetland. Bioresour. Technol. 2017, 244, 407–415. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Q.G.; Li, S.B.; Jia, P.; Qi, C.J.; Ding, F. A review of surface water quality models. Sci. World J. 2013, 2013. [Google Scholar] [CrossRef] [PubMed]
  27. Guinée, J.B. Handbook on life cycle assessment operational guide to the ISO standards. Int. J. Life Cycle Assess. 2001, 7, 311–313. [Google Scholar] [CrossRef]
  28. Scientific Applications International Corporation (SAIC); Curran, M.A. Life-Cycle Assessment: Principles and Practice; National Risk Management Research Laboratory, Office of Research and Development, US Environmental Protection Agency: Cincinnati, OH, USA, 2006.
  29. Rittmann, B.E.; McCarty, P.L. Environmental Biotechnology: Principles and Applications; Tata McGraw-Hill Education: New York, NY, USA, 2012. [Google Scholar]
  30. Wu, J.S.; Ahlert, R.C. Applications of a Steady-State, One-Dimensional Water-Quality Model. Water Resour. Bull. 1979, 15, 660–670. [Google Scholar] [CrossRef]
  31. Lap, B.Q.; Mori, K.; Inoue, E. A one-dimensional model for water quality simulation in medium- and small-sized rivers. Paddy Water Environ. 2007, 5, 5–13. [Google Scholar] [CrossRef]
  32. Sudjono, P.; Koga, K.; Araki, H.; Booij, N. Development and application of one-dimensional water quality model for channel networks and tidal rivers. Hydraul. Eng. Softw. VI 1996, 18, 309–319. [Google Scholar]
  33. Buchak, E.M.; Edinger, J.E. User Guide for LARM2: A Longitudinal-Vertical, Time-Varying Hydrodynamic Reservoir Model; Edinger (Je) Associates Inc.: Wayne, PA, USA, 1982. [Google Scholar]
  34. Martin, J.L. Application of two-dimensional water quality model. J. Environ. Eng. 1988, 114, 317–336. [Google Scholar] [CrossRef]
  35. Covelli, P.; Marsili-Libelli, S.; Pacini, G. SWAMP: A two-dimensional hydrodynamic and quality modeling platform for shallow waters. Numer. Methords Part Differ Equ. 2002, 18, 663–687. [Google Scholar] [CrossRef]
  36. Bellekom, S.; Potting, J.; Benders, R. Feasibility of applying site-dependent impact assessment of acidification in LCA. Int. J. Life Cycle Assess. 2002, 11, 417–424. [Google Scholar] [CrossRef]
  37. Zimmermann, T. Parameterized tool for site specific LCAs of wind energy converters. Int. J. Life Cycle Assess. 2013, 18, 49–60. [Google Scholar] [CrossRef]
  38. Sleeswijk, A.W. Regional LCA in a global perspective. A basis for spatially differentiated environmental life cycle assessment. Int. J. Life Cycle Assess. 2011, 16, 106–112. [Google Scholar] [CrossRef]
  39. Boulay, A.M.; Bulle, C.; Bayart, J.B.; Deschenes, L.; Margni, M. Regional characterization of freshwater use in LCA: Modeling direct impacts on human health. Environ. Sci. Technol. 2011, 45, 8948–8957. [Google Scholar] [CrossRef] [PubMed]
  40. Kim, S.; Dale, B.E. Regional variations in greenhouse gas emissions of biobased products in the United States-corn-based ethanol and soybean oil. Int. J. Life Cycle Assess. 2009, 14, 540–546. [Google Scholar] [CrossRef]
Figure 1. Characterization results for BDO impact category in response to the change of single factor including kCOD, x, and u (1DCCM).
Figure 1. Characterization results for BDO impact category in response to the change of single factor including kCOD, x, and u (1DCCM).
Water 10 01621 g001
Figure 2. Characterization results for BDO impact category in response to the combined change in water velocity (u) and self-purification coefficients (kCOD, kTN) using 1DCCM for the Ashihe river for different months.
Figure 2. Characterization results for BDO impact category in response to the combined change in water velocity (u) and self-purification coefficients (kCOD, kTN) using 1DCCM for the Ashihe river for different months.
Water 10 01621 g002
Figure 3. Longitudinal gradients of characterization results for the BDO impact category at different horizontal sections using 2DCCM for the Songhua river, China. Six horizontal sections were defined: y = 0 m, 50 m, 100 m, 150 m, 200 m, and y = 250 m.
Figure 3. Longitudinal gradients of characterization results for the BDO impact category at different horizontal sections using 2DCCM for the Songhua river, China. Six horizontal sections were defined: y = 0 m, 50 m, 100 m, 150 m, 200 m, and y = 250 m.
Water 10 01621 g003
Figure 4. Comparison of total environmental impact and BDO characterization results among different scenarios at different longitudinal parts of receiving water using 1DCCM in a full LCA analysis relevant to WWTP. Five downstream water body areas are defined here: DWP-1 (distance of 100 m from downstream water body part to emission point), DWP-2 (20,000 m), DWP-3 (40,000 m), DWP-4 (60,000 m), and DWP-5 (80,000 m).
Figure 4. Comparison of total environmental impact and BDO characterization results among different scenarios at different longitudinal parts of receiving water using 1DCCM in a full LCA analysis relevant to WWTP. Five downstream water body areas are defined here: DWP-1 (distance of 100 m from downstream water body part to emission point), DWP-2 (20,000 m), DWP-3 (40,000 m), DWP-4 (60,000 m), and DWP-5 (80,000 m).
Water 10 01621 g004
Table 1. Characterization factors for COD and TN with respect to the impact category of BDO based on different environmental conditions and reference substances (O2, NO3, and PO43−).
Table 1. Characterization factors for COD and TN with respect to the impact category of BDO based on different environmental conditions and reference substances (O2, NO3, and PO43−).
Empirical FormulaEnvironmental ConditionMolecular Weight v C O D v T N B D O C O D   ( O 2 ) B D O T N   ( O 2 ) B D O C O D   ( NO 3 ) B D O T N   ( NO 3 ) B D O C O D   ( PO 4 3 ) B D O T N   ( PO 4 3 )
C5H7O2NMixed culture1130.210.210.38754.42860.59376.7857
C7H12O4NMixed culture1740.1410.1410.26724.42860.40946.7857
C9H15O5NMixed culture2170.1010.1010.20394.42860.31256.7857
C9H16O5NMixed culture2180.1010.1010.19874.42860.30446.7857
C4.9H9.4O2.9NMixed culture1290.1810.1810.34904.42860.53496.7857
C4.7H7.7O2.1NMixed culture1120.2110.2110.40154.42860.61526.7857
C4.9H9O3NMixed culture1300.2010.2010.39544.42860.60586.7857
C5H8.8O3.2NMixed culture1340.2110.2110.39944.42860.61216.7857
C4.1H6.8O2.2NMixed culture1050.2510.2510.49054.42860.75156.7857
C5.1H8.5O2.5NMixed culture1240.1910.1910.37084.42860.56816.7857
C5.3H9.2O2.5NMixed culture1270.1810.1810.34754.42860.53256.7857
C5H8O2NPure culture1140.1910.1910.36904.42860.56546.7857
C5H8.33O0.81NPure culture950.1710.1710.3269 4.42860.50086.7857
C4H8O2NPure culture1020.2410.2410.4559 4.42860.69856.7857
C4.17H7.42O1.38NPure culture940.2210.2210.4226 4.42860.64746.7857
C4.54H7.91O1.95NPure culture1080.2110.2110.4043 4.42860.61946.7857
C4.17H7.21O1.79NPure culture1000.2310.2310.4477 4.42860.68606.7857
C4.16H8O1.25NPure culture920.2110.2110.4049 4.42860.62046.7857
C3.85H6.69O1.78NPure culture950.2610.2610.4990 4.42860.76466.7857
Maximum 2180.2510.2510.49904.42860.76466.7857
Minimum 920.1010.1010.19874.42860.30456.7857
Average 1250.1910.1910.37594.42860.57596.7857
Table 2. Contribution of different impact categories to total impact of each scenario at 3 different downstream: 100 m, 40,000, and 80,000 m.
Table 2. Contribution of different impact categories to total impact of each scenario at 3 different downstream: 100 m, 40,000, and 80,000 m.
Scenario-1Scenario-2Scenario-3
100 m40,000 m80,000 m100 m40,000 m80,000 m100 m40,000 m80,000 m
ADE1.61%1.74%2.21%0.86%1.29%1.75%0.38%0.61%0.90%
ADF6.87%7.41%9.43%3.74%5.59%7.58%2.29%3.67%5.38%
GW8.38%9.05%11.51%4.80%7.16%9.72%2.84%4.57%6.68%
OD0.09%0.09%0.12%0.05%0.07%0.10%0.02%0.03%0.05%
HT1.76%1.90%2.41%0.95%1.41%1.92%0.54%0.87%1.27%
FAET13.74%14.83%18.86%7.29%10.88%14.76%4.26%6.84%10.01%
PO2.34%2.30%3.21%1.20%1.91%0.59%0.78%1.26%1.84%
A7.85%8.47%10.78%4.17%6.23%8.45%2.57%4.13%6.04%
BDO57.37%54.01%41.49%76.85%65.48%53.17%86.31%78.02%67.86%

Share and Cite

MDPI and ACS Style

Bai, S.; Wang, X.; Zhao, X.; Ren, N. Characterizing Water Pollution Potential in Life Cycle Impact Assessment Based on Bacterial Growth and Water Quality Models. Water 2018, 10, 1621. https://doi.org/10.3390/w10111621

AMA Style

Bai S, Wang X, Zhao X, Ren N. Characterizing Water Pollution Potential in Life Cycle Impact Assessment Based on Bacterial Growth and Water Quality Models. Water. 2018; 10(11):1621. https://doi.org/10.3390/w10111621

Chicago/Turabian Style

Bai, Shunwen, Xiuheng Wang, Xinyue Zhao, and Nanqi Ren. 2018. "Characterizing Water Pollution Potential in Life Cycle Impact Assessment Based on Bacterial Growth and Water Quality Models" Water 10, no. 11: 1621. https://doi.org/10.3390/w10111621

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