Next Article in Journal
Exploring Financial Agglomeration and the Impact of Environmental Regulation on the Efficiency of the Green Economy: Fresh Evidence from 30 Regions in China
Next Article in Special Issue
New Insight into the Performance and Self-Defensive Responses of the Algal–Bacterial Granular Sludge Process under Cr(VI)-Induced Stress
Previous Article in Journal
Sustainable Tourism Development and Innovation: Recent Advances and Challenges
Previous Article in Special Issue
Application of Phosphogypsum and Organic Amendment for Bioremediation of Degraded Soil in Tunisia Oasis: Targeting Circular Economy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Key Basic Parameters Involved in Carbon Emissions in Full-Scale Wastewater Treatment Plants

1
Shanghai Investigation, Design & Research Institute Co., Ltd., Shanghai 200335, China
2
School of Environment, Harbin Institute of Technology, Harbin 150090, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(9), 7225; https://doi.org/10.3390/su15097225
Submission received: 13 March 2023 / Revised: 13 April 2023 / Accepted: 24 April 2023 / Published: 26 April 2023
(This article belongs to the Special Issue Wastewater Treatment Technology and Environmental Impact Assessment)

Abstract

:
In this study, carbon emissions in three full-scale wastewater treatment plants were determined by the emission factor method. Moreover, the correlation between basic parameters (influent water parameters and pollutant removal efficiency) and carbon emissions was examined via a structural equation model (SEM). The results showed a significant variation in the total carbon emission intensity of plants over time. The average total carbon emission intensity of plants A, B and C were 0.314, 0.404 and 0.363 kg eqCO2/m³, respectively. Meanwhile, the indirect carbon emission caused by energy and chemical agent consumption accounts for the majority of total carbon emissions (about 85%). Generally, statistical analysis results show that carbon emission intensity is positively correlated with pollutant removal efficiency. Notably, RTN showed the highest positive correlation with Eind, followed by RTN > RCODCr > RTP > TN > RNH3-N > NH3-N > TP. Moreover, capacity showed the greatest negative contribution to Eind, followed by CODCr. In contrast, the positive contribution to Edir was followed by the sequence of RTN > RCODCr > TN > RNH3-N > NH3-N. Notably, CODCr showed a significantly negative correlation with Edir, while TP and its removal showed little correlation with Edir.

1. Introduction

The ongoing emission of anthropogenic greenhouse gas (GHG) is triggering changes in many climate hazards that can impact humanity [1]. Wastewater treatment plants (WWTPs), as essential units of the urban water system, can contribute nearly 1~2% of the total global anthropogenic carbon emission [2]. More critically, carbon emissions continually increase due to the increased discharge of pollutants [3]. Generally, carbon emissions in WWTPs can be divided into direct carbon emissions and indirect carbon emissions [4]. Direct carbon emissions are mainly the GHGs (CO2, CH4 and N2O) discharged during biological reaction processes. In addition, indirect carbon emissions mainly refer to the GHG discharged during the production of consumed energy and chemical agents in wastewater treatment. Briefly, GHG emissions from WWTPs are major contributors to the overall anthropogenic GHG emission, which should be taken seriously [5].
Currently, carbon emission reduction has become an urgently prioritized objective in the wastewater industry [6,7]. Many effective measures such as appropriate process selection, technologies with low carbon emission and precision management have been applied in order to pursue environmental sustainability [8,9,10]. However, the essential prerequisite of carbon emission reduction strategies, which involves the identification of key parameters involved in carbon emissions, still requires further comprehensive research. Generally, the carbon emission of WWTPs varies greatly among influent quality, the discharge standard of pollutants and operation conditions [11]. Moreover, the first two conditions, influent conditions and discharge standards, comprise the criteria of working condition adjustments and process selection [12,13]. The study of the influence of these conditions on carbon emissions is conducive to adjusting operation strategies. Generally, low influent CODCr and TN with a relatively stable C/N ratio demonstrated high carbon emission intensity. By evaluating the carbon emissions of 50 WWTPs in Shanghai, Jiarui, Xi et al. pointed out that the lowest carbon emission was obtained when influent CODCr was 150~250 mg/L and NH3-N was 15~25 mg/L [14]. In contrast, stricter discharge limits led to a higher emission intensity [15]. Highly influent nutrients consume more oxygen and chemical agents, while poorly influent nutrients need extra carbon sources to support the growth of microorganisms, which can remove nitrogen. In particular, the influent C/N (CODCr/TN) ratio is considered one of the most significant parameters because it can markedly affect carbon emissions from both nitrification and denitrification processes [16]. By adjusting the influent C/N, Chen et al. pointed out that an influent C/N ratio of 10 would be optimal for simultaneously achieving relatively higher pollutant removal efficiency and lower GHG emissions in constructed wetlands [17]. However, there are still few comprehensive studies on the correlation between basic parameters (influent water parameters and pollutant removal efficiency) and carbon emissions in full-scale plants.
Usually, the same parameter shows different contributions to different types of carbon emission intensity. Direct emission intensity shows no significant difference between different WWTP scale groups, while indirect carbon emission intensity shows a significant scale effect. The WWTPS at small scales always obtains higher indirect emissions [14]. Meanwhile, some researchers pointed out that the improved scheduling of the influent load can reduce the energy costs of the wastewater treatment plant [18]. Xi Jiarui et al. pointed out that the indirect emission intensity of WWTPs with low CODCr and NH3-N concentrations was twice that of other WWTPs [14]. Notably, the previous study mainly focused on the influence of influent conditions and discharge standards on carbon emissions. Unfortunately, the influence of different processes is not properly excluded.
In addition, the treatment process also has a significant impact on carbon emission intensity. The highest emission intensity is usually reported in the process using membranes with higher energy and consumed chemical agents; for example, 0.79 kgCO2-eq/m3 was reported by a previous study at the Shenzhen MBR plant [19]. Via the remolding process (anaerobic fixed-film MBR reverse osmosis–chlorination process), the carbon emission of the municipal wastewater reclamation could be 0.31 kgCO2-eq/m3 [20]. Nguyen et al. pointed out that the GHG emission intensity of AAO was one-fifth that of SBR in the WWTPs of Australia [21]. In another study, AAO and oxidation ditch processes showed a lower carbon emission per ton of water in WWTPs in Shanghai [14]. Notably, the technologies mostly used in WWTPs in China are AAO and oxidation ditches, which account for over 50% of the existing WWTPs [22]. With the implementation of the Class A Discharge Standard of Pollutants (GB18918-2002) [23] for municipal wastewater treatment plants (WWTPs), treatment processes have been upgraded and rebuilt in many WWTPs in China. An AAO oxidation ditch with larger influent loads, improved TN removal performance and lower energy costs is a typical selected treatment process for implementing upgrades [24]. However, there are few studies on the carbon emission characteristics of the AAO oxidation ditch process. A more comprehensive and in-depth study should be carried out.
In this study, data from three full-scale WWTPs that adopted the AAO oxidation ditch process were collected for conducting operation performance evaluation and carbon emission characteristics analyses. Moreover, structural equation models (SEMs) were adopted to analyze the correlation between carbon emission and the main influencing factors, including capacity, influent water quality and pollutant removal efficiency. This study provides a valuable reference to identify the key parameters involved in the carbon emission of WWTPs.

2. Materials and Methods

2.1. Wastewater Treatment System Description

Carbon emission accounting was carried out at three full-scale tertiary municipal WWTPs in Lu’an, Anhui province. Additionally, the satellite images of studied WWTPs and general information are listed in Table S1.
The influent water in the studied plants comprised mainly domestic sewage, and the effluent can meet the first class, A (from Discharge Limits of Pollutants for Municipal Wastewater Treatment Plant, China (GB18918-2002): CODCr < 50 mg/L, BOD5 < 10 mg/L, TN < 15 mg/L and TP < 0.5 mg/L). As Figure 1 shows, the plants adopt the anaerobic/anoxic/oxic (A/A/O) oxidation ditch process for extensive nitrogen and phosphorus removal. The wastewater and externally returned active sludge first flow into an anaerobic pool, where organic matter is hydrolyzed and acidified to improve biodegradability, and phosphorus accumulation organisms (PAOs) can gain a competitive advantage in the following stage via the synthesis of poly-β-hydroxybutyrate (PHB). Secondly, wastewater flows into the oxidation ditch, and the transformation from an anoxic zone to an aerobic zone is realized by adjusting the guide wall and aeration device. Then, the solid–liquid separation of the effluent is performed in a secondary sedimentation tank. The residual nutrient will be further removed in an activated sand filter. Finally, the effluent is disinfected and discharged into the natural water body.

2.2. Carbon Emission Accounting

Carbon emission accounting is carried out based on the emission factor method, and the parameters were selected by referring to previous studies [25,26]. Carbon emission was calculated with respect to two parts: direct carbon emission and indirect carbon emission. Direct carbon emission, E d i r , is the equivalent CO2 emission that is mainly from CH4 ( E d i r C H 4 ) and N2O ( E d i r N 2 O ) emissions and from biochemical processing units in WWTP. In addition, indirect carbon emission, E i n d , refers to the equivalent CO2 emission during the production of consumed energy and chemical reagents during wastewater treatment. They can be calculated as the following equations, which are provided by the Intergovernmental Panel on Climate Change (IPCC) [25]:
E d i r = G W P C H 4 × E d i r C H 4 + G W P N 2 O × E d i r N 2 O Q
E d i r C H 4 = Q × ( C O D i ) × E F C H 4
E d i r N 2 O = Q × ( T N i ) × E F N 2 O
E i n d = E C i n d × E F C O 2
where GWP refers to the global warming potential within a certain future time period of the GHG. In this study, G W P C H 4 and G W P N 2 O are 25 and 298, separately [27]. Additionally, the adopted emission factor ( E F N 2 O , E F C H 4 and E F C O 2 ) in this study can be found in Table 1 and Table 2.

2.3. Statistical Analysis

R 4.2.1 was used for data processing and figure drawing in this study. The Pearson correlation coefficient between WWTP performance and influent parameters was calculated by using the following equation:
r ( X , Y ) = i = 1 n ( X i X ¯ ) ( Y i Y ¯ ) i = 1 n ( X i X ¯ ) 2 i = 1 n ( Y i Y ¯ ) 2
where X and Y refer to the corresponding observed variables (capacity, CODCr, TN, etc.), r(X,Y) refers to the Pearson correlation coefficient between X and Y, and subscript i represents the number of observations.
Additionally, the structural equation model (SEM) was used to analyze the relationship between basic parameters (influent water parameters and pollutant removal efficiency) and carbon emission intensity. The SEM was constructed using the lavaan 0.6–12 package in R [35], and the model parameter was estimated by the maximum likelihood method. Carbon emissions during the process are related to the capacity, influent water quality (CODCr, TN, NH3-N, TP, etc.) and their removal efficiencies (RCODCr, RTN, RNH3-N, RTP, etc.). All plausible paths between carbon emissions and influent quality were tested. The path coefficient was used to measure the degree of influence or effect between variables. The regression associations implied by the model can be represented by the following:
E i n d = γ 11 c a p a c i t y + γ 12 C O D C r + γ 13 T N + γ 14 N H 3 N + γ 15 T P + γ 16 R C O D C r + γ 17 R T N + γ 18 R N H 3 N + γ 19 R T P + α 1 + ϵ 1
E d i r = γ 21 c a p a c i t y + γ 22 C O D C r + γ 23 T N + γ 24 N H 3 N + γ 25 T P + γ 26 R C O D C r + γ 27 R T N + γ 28 R N H 3 N + γ 29 R T P + α 2 + ϵ 2
E t o t a l = γ 31 E i n d + γ 32 E d i r + α 3 + ϵ 3
where γ refers to the path coefficients, α refers to the constant terms, and ϵ refers to the error terms.
The model evaluation of SEM first involved examining whether the results of parameters estimated in the model had statistical significance, testing the significance of the path coefficient and selecting the path in the model. The chi square/degrees of freedom, non-normed fit index (NNFI), adjusted goodness-of-fit index (AGFI), normed fit index (NFI) and standardized root mean square residual (SRMR) were used to evaluate the successful fit of the model.

3. Results

3.1. Characteristics of Influent Quality

The average daily concentrations of four crucial water quality parameters that correspond to the wastewater influent quality at the studied WWTPs in 2021 are shown in Figure S1. As Figure S1 shows, the monthly median influent CODCr varied from 73.5 to 175.0 mg/L in Plant A, from 88.0 to 170.0 mg/L in Plant B and from 87.0 to 145.0 mg/L in Plant C, which is close to the influent water quality of WWTPs in China reported in previous studies [36,37]. Notably, the sources of wastewater influents are mainly municipal sewage, as well as a possible mixture of stormwater and surface water. In July and August, the CODCr concentration reached an obvious “valley” and the capacity showed a dramatic increase (Figure S1a), which coincided with the food season in Lu’an.
The monthly median influent TN varied from 22.5 to 37.0 mg/L in Plant A, from 17.0 to 38.5 mg/L in Plant B and from 27.0 to 41.5 mg/L in Plant C (Figure S1b). Meanwhile, the monthly median influent NH3-N varied from 18.5 to 34.3 mg/L in Plant A, from 11.8 to 18.9 mg/L in Plant B and from 22.7 to 34.8 mg/L in Plant C (Figure S1c). In this study, NH3-N in influent wastewater accounted for about 48~90%, which is consistent with a previous study [37]. The high proportion of NH3-N in influent wastewater is closely related to the expansion of human activities and industrial processes.
Additionally, the monthly median influent TP was around 2.7~4.6 mg/L in Plant A, 2.3~5.2 mg/L in Plant B and 2.4~4.2 mg/L in Plant C (Figure S1d). Notably, the monthly median capacity varied from 3210 to 4000 m3/h in Plant A, 748 to 1255 m3/h in Plant B and 1880 to 2160 m3/h in Plant C, which showed a different trend compared with the influent water quality. In accordance with CODCr, the TN, NH3-N and TP exhibited their minimum value in August. The obvious “valley” is mainly for the rainfall effects on inflow and infiltration in wastewater treatment systems [38], and mixed rainwater may greatly influence the subsequent wastewater treatment process. The impact of influent wastewater quality on carbon emissions will be discussed in detail in Section 3.4.

3.2. The Correlation between WWTP Performance and Influent Parameters

The major pollutant removal rate is shown in Figure S2. The monthly median CODCr removal rate differed from 82.0% to 92.5% in Plant A, from 87.5% to 92.5% in Plant B and from 88.0% to 94.0% in Plant C (Figure S2a). The monthly median TN removal rate varied from 64.5% to 80% in Plant A, from 68.0% to 85.0% in Plant B and from 63.0% to 74.8% in Plant C (Figure S2b). Additionally, the monthly median NH3-N removal rate differed from 90.5% to 99.3% in Plant A, from 92.5% to 98.5% in Plant B and from 96.8% to 99.6% in Plant C (Figure S2c). In addition, the monthly median TP removal rate differed from 89.0% to 94.2% in Plant A, from 91.6% to 97.5% in Plant B and from 89.2% to 93.5% in Plant C (Figure S2d). Notably, the major nutrient removal rate was consistent with the influent CODCr. In the influent “valley” (July and August), nutrient removal became unstable and exhibited poor performance.
The correlation between the removal rate and influent parameters is shown in Figure 2. Apparently, there is a negative correlation between capacity and influent water quality, which is further evidence of the rainfall effects on inflow and infiltration in wastewater treatment systems. Additionally, the influent TN showed a significant negative correlation with the removal of TN (RTN). This is mainly because the higher influent TN needs more carbon sources for denitrification, which is the main denitrification path in this study. Thus, RTN will decrease if there is no sufficient carbon source in the influent wastewater or the external carbon source is not obtained in a timely and precise manner. Interestingly, the TP removal rate (RTP) increased with the increase in influent TP, which was mainly due to the fluctuation of influent C/P from 23 to 33, and the variation in influent phosphorus concentrations had little effect on biological phosphorus removal efficiency. With the subsequent chemical phosphorus removal and the increase in phosphorus mass concentration, the mass concentration of iron–phosphorus precipitates in the effluent increases, and the chance of dissolved phosphorus being complexed by precipitate adsorption increases greatly, which will improve the phosphorus removal rate.
In addition, there was a negative correlation between RTP and RTN, which is mainly for the competition between phosphate-accumulating organisms (PAOs) and denitrifying bacteria for carbon sources [39]. If the carbon source, which is the electron donor during denitrification, is insufficient due to the interference of biological phosphorus removal, various denitrifying enzymes will compete for electrons, leading to a decrease in denitrification efficiency and the accumulation of N2O [40].

3.3. Carbon Emission Characteristics of WWTPs

The different types of carbon emissions of the studied WWTPs in Lu’an were accounted, and the results are shown in Figure 3. The total carbon emission intensity ( E t o t a l ) reached a relative peak in the summer, which is due to the sudden increase in influent load and relatively low influent nutrients during the rainy season. The average monthly total carbon emission intensity is about 0.273~0.334 kg eqCO2/m³ in Plant A, 0.301~0.514 kg eqCO2/m³ in Plant B and 0.322~0.413 kg eqCO2/m³ in Plant C, which is close to previous studies [14,41]. Plant B, with a lower capacity (20,000 m3/d), exhibits a higher total carbon emission intensity, which should be due to the scale effect [42].
As Figure 3 shows, major carbon emission is caused by electricity consumption ( E i n d e l e ), which accounts for 44.3~77.7% in Plant A, 43.7~83.4% in Plant B and 36.4~70.7% in Plant C. Then, carbon emissions introduced by chemical agent consumption account for 4.9~35.6% in Plant A, 2.9~40.6% in Plant B and 12.9~46.9% in Plant C. It is important to note that there is an excessive variation in chemical agent consumption relative to influent water quality and organic removal efficiency, which indicates a further improvement in the management of chemical agent dosing. That is, the refined management of chemical agent dosing will be one of the priorities of carbon emission reduction in these WWTPs. Collectively, the indirect carbon emission accounts for 79.4~88.4% in Plant A, 81.7~92.1% in Plant B and 81.7~87.7% in Plant C.
While direct carbon emission accounts for 11.6~20.6% in Plant A, 7.8~18.3% in Plant B and 12.2~18.3% in Plant C, it should be noted that direct carbon emissions due to N2O emissions, E d i r N 2 O , comprise a major part of direct carbon emissions. This is mainly related to the adopted treatment process, which is mainly based on the aerobic environment.
Notably, the proportion of different types of carbon emissions also shows clear differences over time. E i n d e l e comprises a relatively large portion during the winter. At lower ambient temperatures, the activity of microorganisms decreased significantly, and increasing dissolved oxygen (DO) to maintain stable performance is the most common and preferred strategy [43], which leads to more power consumption. Meanwhile, carbon emissions caused by chemicals consumption ( E i n d c h e m i c a l s ) account for 2.9~46.9%. The high chemical cost mainly occurs in the summer and autumn seasons, which can be attributed to the low C/N and TP in influent water [44]. Significantly, Plant B exhibited far more carbon emissions (0.404 kg eqCO2/m³) in September than in other months. This is mainly attributed to the high hydraulic load (~130% design capacity). In addition, the high hydraulic loading rate is a common problem in WWTPs in China [36]. Under a higher hydraulic load, WWTPs would have little room for regulation and cannot effectively cope with changes in wastewater influent quality and capacity, eventually leading to more carbon emissions.

3.4. SEM Analysis on Basic Parameters Involved in Carbon Emissions

Regarding the wide range of parameters related to carbon emissions, SEM was used to identify the degree of correlation. SEM (Figure 4) shows the influence of the influent water quality and pollutant removal rate on carbon emissions. The non-normed fit index (NNFI) = 0.907, adjusted goodness-of-fit index (AGFI) = 0.997 and chi square/degrees of freedom = 1~2 indicate a good fit of the model to the pure data [45,46].
In line with previous studies, it can be observed that a significantly strong relationship (0.85 **) exists between indirect carbon emission (Eind) and total carbon emission (Etotal), which indicates that indirect carbon emissions account for large proportions of the total carbon emissions of WWTPs [47]. Notably, Edir and Eind have a significant positive inter-relationship. That is, when Eind is higher, Edir is usually higher as well. An improper operational control strategy results in useless energy and chemical agent consumption besides more greenhouse gas emissions. For instance, excessive dissolved oxygen (DO) requires more electricity, which results in more Eind. Moreover, it can inhibit the production and activity of denitrification enzymes, and nitrous oxide reductase (NOR) is more sensitive to the fluctuations of DO than other denitrification enzymes, which leads to increased nitrous oxide accumulation [48].
Collectively, influent parameters have different impacts on direct carbon emissions (Edir), indirect carbon emissions (Eind) and total carbon emissions (Etotal). Higher concentrations of TN result in higher carbon emission intensity, which involves both Eind and Edir. In contrast, higher concentrations of CODCr during the study period are more beneficial for pollutant removal [49], which decreases Eind and Edir. Capacity has a significant negative effect on Eind and Etotal, which should be attributed to the scale effect [42,50]. In this study, Plant B shows a higher Etotal with a significant fluctuation under the lowest capacity (Figure 3). Notably, TP and its removal RTP are mainly related to the chemical agent’s dosage during tertiary treatments, which is an important part of Eind. Meanwhile, they showed little correlation with Edir in the SEM analysis. Notably, phosphorus removal is linked to nitrogen removal (Figure 2). Murnleitner et al. pointed out that there is competition between phosphate-accumulating organisms (PAOs) and denitrifying bacteria for carbon sources [39], which will further impact direct carbon emissions. All concerned influent parameters and pollutant removal efficiency have a direct positive contribution to carbon emissions during wastewater treatment processes, with the exception of capacity. Influent parameters show different influences on direct and indirect carbon emissions. Among these, capacity is significantly negatively associated with Eind (−0.54 **) and has no significant correlation with Edir. This indicates that capacity is a key factor affecting the carbon emissions of WWTPs due to the scale effect. Conclusively, RTN and RCODCr are singled out as the significant positive influential factors for Edir and Eind, while CODCr shows a negative influence on Edir and Eind. The positive contribution to Edir followed the sequence of RTN > RCODCr > TN > RNH3-N > NH3-N, while the positive contribution to Eind followed the sequence of RTN > RCODCr > RTP > TN > RNH3-N > NH3-N > TP.

4. Discussion

Generally, the average influent concentrations of the CODCr, NH3-N, TN and TP of the studied WWTPs in 2021 were 124.13, 22.07, 29.69 and 3.40 mg/L, respectively. In this study, the average removal rates of CODCr, NH3-N, TN and TP were 88.7%, 97.2%, 71.9% and 93.0%, respectively, which showed good performance for pollutant removal. Additionally, the average total carbon emission intensity of the AAO oxidation ditch was about 0.314~0.404 kg eqCO2/m³, which is slightly above the total carbon emission intensity of the oxidation ditch (0.318 kg eqCO2/m³) in the previous study [14]. This is mainly due to the seriously insufficient carbon source in wastewater influents with a C/N ratio of around 4.2, which is lower than the minimum value in the study by Sun et al. [37]. The insufficient carbon source is crucial to BNR processes as CODCr acts as a limiting factor for phosphorus release and denitrification. This leads to additional carbon sources for nitrogen and phosphorus removal, which introduces higher carbon emissions caused by chemical agent costs (0.04 kg eqCO2/m³).
Although there have been many studies focused on carbon neutrality and energy self-sufficiency in WWTPs [51], most WWTPs still require considerable energy and chemical agents to remove pollutants. Among them, RTN shows the highest positive correlation with Edir and Eind. Nitrogen removal has been known as an important challenge for WWTPs’ performance [52]. Moreover, there is difficulty in coping with nitrogen removal during the WWTPs’ carbon emission reduction. The removal of nitrogen in an efficient and low-carbon manner is an urgent problem that needs to be overcome. Additionally, phosphorus removal has also been another key challenge for WWTPs. In our study, RTP is mainly positively related to Eind but shows little correlation with Edir.
The SEM results indicated that there was a significant correlation between pollutant removal and carbon emission. Guo et al. pointed out that there should be a trade-off between carbon emission and the pollutant removal of activated sludge processes, which needs further research [53]. In particular, the removal of TN has the greatest impact on carbon emissions (Figure 4), and the adoption of low carbon emission nitrogen removal processes is also one of the important ways to realize carbon emission reductions in the wastewater industry [54]. Notably, the capacity also shows a significant scale effect on the carbon emission per ton of water, which was also proved in a previous study on 217 wastewater treatment plants in the Valencia region [42]. In general, the AAO oxidation ditch has a high application potential, but its application conditions still need to be further studied according to influent conditions and pollutant removal requirements.

5. Conclusions

The carbon emissions of WWTPs are rather complicated due to the variety of GHGs, their generation and the wide range of influencing parameters. The main objective of this study is to identify key basic parameters involved in full-scale WWTPs, such as capacity, influent water quality and pollutant removal efficiency. This study is inspiring in a certain sense with respect to forming clear judgments and thus creating appropriate designs, operations and the management of carbon emission reductions in wastewater treatment processes. The primary conclusions have been summarized as follows:
  • In 2021, the average removal rates of CODCr, NH3-N, TN and TP were 86.1%, 97.4%, 71.0% and 91.5% in Plant A; 89.6%, 96.1%, 75.8% and 95.1% in Plant B; and 90.3%, 98.1%, 69.0% and 92.5% in Plant C. These showed good performances for the AAO oxidation ditch with respect to pollutant removal.
  • Carbon emissions during wastewater treatment processes mainly consist of indirect carbon emissions (~90%). The SEM results show that influent CODCr and TN and their removal should be key indicators related to carbon emissions. For domestic sewage, a higher influent organic matter concentration helps reduce energy and agent consumption and greenhouse gas emissions.
  • The SEM results indicated that the positive contribution to Eind followed the sequence of RTN > RCODCr > RTP > RNH3-N > TN > NH3-N > TP. Notably, capacity showed a significant negative contribution to Eind. Additionally, capacity showed the highest negative correlation with Eind, followed by CODCr, while the contribution to Edir followed the sequence of RTN > RCODCr > TN > RNH3-N > NH3-N. Notably, CODCr showed a significantly negative correlation with Edir.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/su15097225/s1, Table S1. Images of the WWTPs (from Tiandi Maps) and general information. Figure S1. The water quality parameters variation of influent water: (a) CODCr, (b) TN, (c) NH3-N, (d) TP and (e) capacity. Figure S2. Removal rate of major pollutants: (a) CODCr, (b) TN, (c) NH3-N and (d) TP.

Author Contributions

Conceptualization, H.Y. and H.L.; methodology, K.G.; validation, K.G.; formal analysis, K.G.; investigation, K.G., H.Y. and Q.Z.; resources, H.Y.; data curation, K.G.; writing—original draft preparation, K.G.; writing—review and editing, H.L., Q.Z. and H.Y.; supervision, H.Y.; project administration, H.Y.; funding acquisition, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Research Project of China Three Gorges Corporation (No. 202103547), the Research Project of Shanghai Investigation, Design & Research Institute Co., Ltd. (No. 2022QT(831)-001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors greatly appreciate the financial support from “Research on carbon reduction technology based on AAMBBR combined magnetic coagulation” (2022QT(831)-001) and the research project of China Three Gorges Group Co. (Contract No. 202103547).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mora, C.; Spirandelli, D.; Franklin, E.C.; Lynham, J.; Kantar, M.B.; Miles, W.; Smith, C.Z.; Freel, K.; Moy, J.; Louis, L.V.; et al. Broad threat to humanity from cumulative climate hazards intensified by greenhouse gas emissions. Nat. Clim. Chang. 2018, 8, 1062–1071. [Google Scholar] [CrossRef]
  2. Ritchie, H.; Roser, M.; Rosado, P. CO2 and Greenhouse Gas Emissions. Available online: https://ourworldindata.org/co2-and-greenhouse-gas-emissions (accessed on 20 February 2023).
  3. Su, H.; Yi, H.; Gu, W.; Wang, Q.; Liu, B.; Zhang, B. Cost of raising discharge standards: A plant-by-plant assessment from wastewater sector in China. J. Environ. Manag. 2022, 308, 114642. [Google Scholar] [CrossRef] [PubMed]
  4. Li, L.; Wang, X.; Miao, J.; Abulimiti, A.; Jing, X.; Ren, N. Carbon neutrality of wastewater treatment—A systematic concept beyond the plant boundary. Environ. Sci. Ecotechnol. 2022, 11, 100180. [Google Scholar] [CrossRef]
  5. Du, W.-J.; Lu, J.-Y.; Hu, Y.-R.; Xiao, J.; Yang, C.; Wu, J.; Huang, B.; Cui, S.; Wang, Y.; Li, W.-W. Spatiotemporal pattern of greenhouse gas emissions in China’s wastewater sector and pathways towards carbon neutrality. Nat. Water 2023, 1, 166–175. [Google Scholar] [CrossRef]
  6. Zhang, Y.; Ni, X.; Wang, H. Visual analysis of greenhouse gas emissions from sewage treatment plants based on CiteSpace: From the perspective of bibliometrics. Environ. Sci. Pollut. Res. 2023, 30, 45555–45569. [Google Scholar] [CrossRef]
  7. You, X.; Yang, L.; Zhou, X.; Zhang, Y. Sustainability and carbon neutrality trends for microalgae-based wastewater treatment: A review. Environ. Res. 2022, 209, 112860. [Google Scholar] [CrossRef]
  8. Liang, W.; Yu, C.; Ren, H.; Geng, J.; Ding, L.; Xu, K. Minimization of nitrous oxide emission from CASS process treating low carbon source domestic wastewater: Effect of feeding strategy and aeration rate. Bioresour. Technol. 2015, 198, 172–180. [Google Scholar] [CrossRef]
  9. Wu, Z.; Duan, H.; Li, K.; Ye, L. A comprehensive carbon footprint analysis of different wastewater treatment plant configurations. Environ. Res. 2022, 214, 113818. [Google Scholar] [CrossRef]
  10. Kiyan, A.; Gheibi, M.; Akrami, M.; Moezzi, R.; Behzadian, K.; Taghavian, H. The Operation of Urban Water Treatment Plants: A Review of Smart Dashboard Frameworks. Environ. Ind. Lett. 2023, 1, 28–45. [Google Scholar] [CrossRef]
  11. Hua, H.; Jiang, S.; Yuan, Z.; Liu, X.; Zhang, Y.; Cai, Z. Advancing greenhouse gas emission factors for municipal wastewater treatment plants in China. Environ. Pollut. 2022, 295, 118648. [Google Scholar] [CrossRef]
  12. Bozkurt, H.; van Loosdrecht, M.C.M.; Gernaey, K.V.; Sin, G. Optimal WWTP process selection for treatment of domestic wastewater—A realistic full-scale retrofitting study. Chem. Eng. J. 2016, 286, 447–458. [Google Scholar] [CrossRef]
  13. Keller, J.; Hartley, K. Greenhouse gas production in wastewater treatment: Process selection is the major factor. Water Sci. Technol. 2003, 47, 43–48. [Google Scholar] [CrossRef] [PubMed]
  14. Xi, J.; Gong, H.; Zhang, Y.; Dai, X.; Chen, L. The evaluation of GHG emissions from Shanghai municipal wastewater treatment plants based on IPCC and operational data integrated methods (ODIM). Sci. Total Environ. 2021, 797, 148967. [Google Scholar] [CrossRef] [PubMed]
  15. Neethling, J.; Falk, M.W.; Reardon, D.J.; Clark, D.L.; Pramanik, A. WERF Nutrient Challenge—Nutrient Regulations, Treatment Performance, and Sustainability Collide. In Proceedings of the Water Environment Federation, Los Angeles, CA, USA, 15–19 October 2011; pp. 1–16. [Google Scholar]
  16. Zheng, M.; Zhou, N.; He, S.; Chang, F.; Zhong, J.; Xu, S.; Wang, Z.; Liu, T. Nitrous oxide (N2O) emissions from a pilot-scale oxidation ditch under different COD/N ratios, aeration rates and two shock-load conditions. J. Environ. Manag. 2021, 280, 111657. [Google Scholar] [CrossRef]
  17. Chen, X.; Zhu, H.; Yan, B.; Shutes, B.; Tian, L.; Wen, H. Optimal influent COD/N ratio for obtaining low GHG emissions and high pollutant removal efficiency in constructed wetlands. J. Clean. Prod. 2020, 267, 122003. [Google Scholar] [CrossRef]
  18. Simon-Várhelyi, M.; Cristea, V.M.; Luca, A.V. Reducing energy costs of the wastewater treatment plant by improved scheduling of the periodic influent load. J. Environ. Manag. 2020, 262, 110294. [Google Scholar] [CrossRef]
  19. Liao, X.; Tian, Y.; Gan, Y.; Ji, J. Quantifying urban wastewater treatment sector’s greenhouse gas emissions using a hybrid life cycle analysis method—An application on Shenzhen city in China. Sci. Total Environ. 2020, 745, 141176. [Google Scholar] [CrossRef]
  20. Wang, S.; Liu, H.; Gu, J.; Zhang, M.; Liu, Y. Towards carbon neutrality and water sustainability: An integrated anaerobic fixed-film MBR-reverse osmosis-chlorination process for municipal wastewater reclamation. Chemosphere 2022, 287, 132060. [Google Scholar] [CrossRef]
  21. Nguyen, T.K.L.; Ngo, H.H.; Guo, W.; Chang, S.W.; Nguyen, D.D.; Nghiem, L.D.; Liu, Y.; Ni, B.; Hai, F.I. Insight into greenhouse gases emissions from the two popular treatment technologies in municipal wastewater treatment processes. Sci. Total Environ. 2019, 671, 1302–1313. [Google Scholar] [CrossRef]
  22. Zhang, Q.H.; Yang, W.N.; Ngo, H.H.; Guo, W.S.; Jin, P.K.; Dzakpasu, M.; Yang, S.J.; Wang, Q.; Wang, X.C.; Ao, D. Current status of urban wastewater treatment plants in China. Environ. Int. 2016, 92–93, 11–22. [Google Scholar] [CrossRef]
  23. GB 18918-2002; Discharge Standard of Pollutants for Municipal Wastewater Treatment Plant. Minstry of Ecology and Environment of the Pepole’s Republic of China: Beijing, China, 2002. Available online: https://www.mee.gov.cn/ywgz/fgbz/bz/bzwb/shjbh/swrwpfbz/200307/t20030701_66529.shtml (accessed on 20 February 2023).
  24. Lou, T.; Peng, Z.; Jiang, K.; Niu, N.; Wang, J.; Liu, A. Nitrogen removal characteristics of biofilms in each area of a full-scale AAO oxidation ditch process. Chemosphere 2022, 302, 134871. [Google Scholar] [CrossRef] [PubMed]
  25. IPCC. 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories; IPCC: Geneva, Switzerland, 2019. [Google Scholar]
  26. China urban water association. Guidelines for Carbon Accounting and Emission Reduction in the Urban Water Sector, 1st ed.; China Architecture Publishing & Media Co., Ltd.: Beijing, China, 2022. [Google Scholar]
  27. IPCC; TEAP. Safeguarding the Ozone Layer and the Global Climate System; Environmental Policy Collection; IPCC: Geneva, Switzerland, 2005. [Google Scholar]
  28. Yan, X.; Li, L.; Liu, J. Characteristics of greenhouse gas emission in three full-scale wastewater treatment processes. J. Environ. Sci. 2014, 26, 256–263. [Google Scholar] [CrossRef] [PubMed]
  29. Ren, Y.G.; Wang, J.H.; Li, H.F.; Zhang, J.; Qi, P.Y.; Hu, Z. Nitrous oxide and methane emissions from different treatment processes in full-scale municipal wastewater treatment plants. Environ. Technol. 2013, 34, 2917–2927. [Google Scholar] [CrossRef] [PubMed]
  30. Masuda, S.; Sano, I.; Hojo, T.; Li, Y.-Y.; Nishimura, O. The comparison of greenhouse gas emissions in sewage treatment plants with different treatment processes. Chemosphere 2018, 193, 581–590. [Google Scholar] [CrossRef]
  31. Foley, J.; de Haas, D.; Yuan, Z.; Lant, P. Nitrous oxide generation in full-scale biological nutrient removal wastewater treatment plants. Water Res. 2010, 44, 831–844. [Google Scholar] [CrossRef]
  32. Valkova, T.; Parravicini, V.; Saracevic, E.; Tauber, J.; Svardal, K.; Krampe, J. A method to estimate the direct nitrous oxide emissions of municipal wastewater treatment plants based on the degree of nitrogen removal. J. Environ. Manag. 2021, 279, 111563. [Google Scholar] [CrossRef]
  33. Li, H.; Peng, D.; Liu, W.; Wei, J.; Wang, Z.; Wang, B. N2O generation and emission from two biological nitrogen removal plants in China. Desalination Water Treat. 2016, 57, 11800–11806. [Google Scholar] [CrossRef]
  34. Ministry of Ecology and Environment. The Baseline Emission Factor of China’s Regional Power Grid in 2019 Emission Reduction Project; Ministry of Ecology and Environment: Beijing, China, 2019. [Google Scholar]
  35. Rosseel, Y. Lavaan: An R Package for Structural Equation Modeling. J. Stat. Softw. 2012, 48, 1–36. [Google Scholar] [CrossRef]
  36. Zhang, J.; Shao, Y.; Wang, H.; Liu, G.; Qi, L.; Xu, X.; Liu, S. Current operation state of wastewater treatment plants in urban China. Environ. Res. 2021, 195, 110843. [Google Scholar] [CrossRef]
  37. Sun, Y.; Chen, Z.; Wu, G.; Wu, Q.; Zhang, F.; Niu, Z.; Hu, H.-Y. Characteristics of water quality of municipal wastewater treatment plants in China: Implications for resources utilization and management. J. Clean. Prod. 2016, 131, 1–9. [Google Scholar] [CrossRef]
  38. Cahoon, L.B.; Hanke, M.H. Rainfall effects on inflow and infiltration in wastewater treatment systems in a coastal plain region. Water Sci. Technol. 2017, 75, 1909–1921. [Google Scholar] [CrossRef] [PubMed]
  39. Murnleitner, E.; Kuba, T.; van Loosdrecht, M.C.M.; Heijnen, J.J. An integrated metabolic model for the aerobic and denitrifying biological phosphorus removal. Biotechnol. Bioeng. 1997, 54, 434–450. [Google Scholar] [CrossRef]
  40. Zumft, W.G. Cell biology and molecular basis of denitrification. Microbiol. Mol. Biol. Rev. 1997, 61, 533–616. [Google Scholar] [CrossRef] [PubMed]
  41. Zeng, S.; Chen, X.; Dong, X.; Liu, Y. Efficiency assessment of urban wastewater treatment plants in China: Considering greenhouse gas emissions. Resour. Conserv. Recycl. 2017, 120, 157–165. [Google Scholar] [CrossRef]
  42. Hernández-Chover, V.; Bellver-Domingo, Á.; Hernández-Sancho, F. Efficiency of wastewater treatment facilities: The influence of scale economies. J. Environ. Manag. 2018, 228, 77–84. [Google Scholar] [CrossRef] [PubMed]
  43. Zhou, H.; Li, X.; Xu, G.; Yu, H. Overview of strategies for enhanced treatment of municipal/domestic wastewater at low temperature. Sci. Total Environ. 2018, 643, 225–237. [Google Scholar] [CrossRef]
  44. Machat, H.; Boudokhane, C.; Roche, N.; Dhaouadi, H. Effects of C/N Ratio and DO concentration on Carbon and Nitrogen removals in a Hybrid Biological Reactor. Biochem. Eng. J. 2019, 151, 107313. [Google Scholar] [CrossRef]
  45. Hu, L.t.; Bentler, P.M. Cutoff criteria for fit indexes in covariance structure analysis: Conventional criteria versus new alternatives. Struct. Equ. Model. Multidiscip. J. 1999, 6, 1–55. [Google Scholar] [CrossRef]
  46. Hair, J.; Black, W.C.; Babin, B.J.; Anderson, R.E.; Tatham, R.L. SEM: An Introduction. In Multivariate Data Analysis: A Global Perspective, 7th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2010. [Google Scholar]
  47. Xu, R.; Li, Y.; Luo, Y.; Fang, F.; Feng, Q.; Cao, J.; Luo, J. Prediction and Evaluation of Indirect Carbon Emission from Electrical Consumption in Multiple Full-Scale Wastewater Treatment Plants via Automated Machine Learning-Based Analysis. ACS EST Eng. 2022, 3, 360–372. [Google Scholar] [CrossRef]
  48. Otte, S.; Grobben, N.G.; Robertson, L.A.; Jetten, M.S.; Kuenen, J.G. Nitrous oxide production by Alcaligenes faecalis under transient and dynamic aerobic and anaerobic conditions. Appl. Environ. Microbiol. 1996, 62, 2421–2426. [Google Scholar] [CrossRef]
  49. Benckiser, G.; Eilts, R.; Linn, A.; Lorch, H.J.; Sümer, E.; Weiske, A.; Wenzhöfer, F. N2O emissions from different cropping systems and from aerated, nitrifying and denitrifying tanks of a municipal waste water treatment plant. Biol. Fertil. Soils 1996, 23, 257–265. [Google Scholar] [CrossRef]
  50. Molinos-Senante, M.; Sala-Garrido, R.; Iftimi, A. Energy intensity modeling for wastewater treatment technologies. Sci. Total Environ. 2018, 630, 1565–1572. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, L.; Ling, J.; Lin, M. Carbon neutrality: A comprehensive bibliometric analysis. Environ. Sci. Pollut. Res. 2023, 30, 45498–45514. [Google Scholar] [CrossRef] [PubMed]
  52. McCarty, P.L. What is the Best Biological Process for Nitrogen Removal: When and Why? Environ. Sci. Technol. 2018, 52, 3835–3841. [Google Scholar] [CrossRef] [PubMed]
  53. Guo, J.; Fu, X.; Andrés Baquero, G.; Sobhani, R.; Nolasco, D.A.; Rosso, D. Trade-off between carbon emission and effluent quality of activated sludge processes under seasonal variations of wastewater temperature and mean cell retention time. Sci. Total Environ. 2016, 547, 331–344. [Google Scholar] [CrossRef] [PubMed]
  54. Cogert, K.I.; Ziels, R.M.; Winkler, M.K.H. Reducing Cost and Environmental Impact of Wastewater Treatment with Denitrifying Methanotrophs, Anammox, and Mainstream Anaerobic Treatment. Environ. Sci. Technol. 2019, 53, 12935–12944. [Google Scholar] [CrossRef]
Figure 1. Simplified wastewater treatment process diagram of studied WWTPs.
Figure 1. Simplified wastewater treatment process diagram of studied WWTPs.
Sustainability 15 07225 g001
Figure 2. The correlation heatmap between influent parameters, operation parameters and (data with p-value > 0.05 are already hidden) influent performance parameters: capacity, CODCr, TN, NH3-N and TP; removal efficiency: RCODCr, RTN, RTP and RNH3-N. “**”and ”***” suggest that P value is less than 0.01 and 0.001, respectively. And the “×” suggests that the correlation is not significant.
Figure 2. The correlation heatmap between influent parameters, operation parameters and (data with p-value > 0.05 are already hidden) influent performance parameters: capacity, CODCr, TN, NH3-N and TP; removal efficiency: RCODCr, RTN, RTP and RNH3-N. “**”and ”***” suggest that P value is less than 0.01 and 0.001, respectively. And the “×” suggests that the correlation is not significant.
Sustainability 15 07225 g002
Figure 3. The variation in total carbon emission intensity and the proportion of different types of carbon emissions. Plants A–C are three full-scale wastewater treatment plants in Anhui province and the detailed information can be found in Table S1.
Figure 3. The variation in total carbon emission intensity and the proportion of different types of carbon emissions. Plants A–C are three full-scale wastewater treatment plants in Anhui province and the detailed information can be found in Table S1.
Sustainability 15 07225 g003
Figure 4. The path diagram of SEM analysis. Single arrows denote the causal relationship from one variable to another, and double-headed arrows denote the association between two variables. Numbers adjacent to arrows are path coefficients, which represent the strength of the relationship. Red lines indicate that the correlation is positive, and dark blue lines indicate that the correlation is negative. * represents p < 0.05; ** represents p < 0.01. The SEM model fits well (chi square = 16.615; degrees of freedom (DF) = 9; standardized root mean square residual (SRMR) = 0.013).
Figure 4. The path diagram of SEM analysis. Single arrows denote the causal relationship from one variable to another, and double-headed arrows denote the association between two variables. Numbers adjacent to arrows are path coefficients, which represent the strength of the relationship. Red lines indicate that the correlation is positive, and dark blue lines indicate that the correlation is negative. * represents p < 0.05; ** represents p < 0.01. The SEM model fits well (chi square = 16.615; degrees of freedom (DF) = 9; standardized root mean square residual (SRMR) = 0.013).
Sustainability 15 07225 g004
Table 1. Description and emission factor value of GHG in previous studies.
Table 1. Description and emission factor value of GHG in previous studies.
This StudyPrevious Studies
Emission FactorValueValueDescriptionReference
E F C H 4 0.0030 kg CH4/kg CODinfluent0.00306Orbal oxidation ditch in Beijing WWTPs[28]
0.0079pre-anaerobic carrousel oxidation ditch in Jinan[29]
0.00133oxidation ditch in the Akiu sewage treatment plant in Sendai city, Japan[30]
E F N 2 O 0.0017 kg N2O/kg TNinfluent0.00571aeration oxidation ditch in Brisbane, Queensland[31]
0.00014oxidation ditch in Akiu sewage treatment plant in Sendai city, Japan[30]
0.00173Orbal oxidation ditch in Beijing WWTPs[28]
0.000295plug-flow AS tank[32]
0.00037~0.0015Orbal oxidation ditch in Xi’an, No.3 WWTP[33]
Table 2. Equivalent CO2 emission factor value for consumed energy and chemical reagent production.
Table 2. Equivalent CO2 emission factor value for consumed energy and chemical reagent production.
Consumed Energy and Chemical ReagentValueReference
electricity0.7921 kg CO2-eq/kWh[34]
polyaluminium chloride0.53 kg CO2-eq/kg[26]
sodium acetate0.623 kg CO2-eq/kg[26]
ferric chloride0.26 kg CO2-eq/kg[26]
sodium hypochlorite0.99 kg CO2-eq/kg[26]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, K.; Yang, H.; Zhao, Q.; Liu, H. Identification of Key Basic Parameters Involved in Carbon Emissions in Full-Scale Wastewater Treatment Plants. Sustainability 2023, 15, 7225. https://doi.org/10.3390/su15097225

AMA Style

Gao K, Yang H, Zhao Q, Liu H. Identification of Key Basic Parameters Involved in Carbon Emissions in Full-Scale Wastewater Treatment Plants. Sustainability. 2023; 15(9):7225. https://doi.org/10.3390/su15097225

Chicago/Turabian Style

Gao, Kuo, Hong Yang, Qingliang Zhao, and Haichen Liu. 2023. "Identification of Key Basic Parameters Involved in Carbon Emissions in Full-Scale Wastewater Treatment Plants" Sustainability 15, no. 9: 7225. https://doi.org/10.3390/su15097225

APA Style

Gao, K., Yang, H., Zhao, Q., & Liu, H. (2023). Identification of Key Basic Parameters Involved in Carbon Emissions in Full-Scale Wastewater Treatment Plants. Sustainability, 15(9), 7225. https://doi.org/10.3390/su15097225

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