1. Introduction
The motivation for energy neutrality, resource recovery, and decarbonization in the water sector has led to the exploration of novel wastewater treatment schemes towards carbon capture, energy recovery, and efficient nutrient removal and recovery with less energy/external carbon requirements and carbon emissions (Ren and Pagilla 2022, Wan et al. 2016). At present, several energy-efficient and even energy-positive technologies, such as shortcut nitrogen (N) removal processes (Lackner et al. 2014, Ma et al. 2016), coupled aerobic-anoxic nitrous decomposition operation (CANDO) (Scherson et al. 2014), microalgae cultivation (Chen et al. 2015), and anaerobic membrane bioreactors (Smith et al. 2012), have been proposed to replace traditional biological N removal processes (i.e., nitrification-denitrification) in municipal wastewater treatment plants (WWTPs) that are typically energy intensive and require additional carbon. Among them, the adsorption/bio-oxidation (A/B) process, where A-stage diverts most organic matter to anaerobic digestion for carbon capture and bioenergy recovery, while the subsequent B-stage aims for carbon- and energy-efficient N removal, is one of the promising engineering innovations that can meet the needs of energy and nutrient management (Guthi et al. 2022, Liu et al. 2019). Specifically, at the B-stage, the partial nitrification/anammox (PN/A) process, which relies on the cooperation of ammonium-oxidizing bacteria (AOB) and anaerobic AOB (AnAOB) and the suppression of nitrite-oxidizing bacteria (NOB) for completely autotrophic N removal via nitrite pathway, has attracted extensive attentions due to its largely reduced demands for aeration energy, external carbon, and waste sludge treatment (Cao et al. 2017, Ma et al. 2016). Side-stream PN/A has been widely studied and employed in many WWTPs for high-strength digestor centrate treatment, while mainstream PN/A application is still limited and faces challenges in stability (Cao et al. 2017, Lackner et al. 2014, Li et al. 2020b). Alternatively, partial denitrification/anammox (PDN/A) was demonstrated as a more practically viable route for mainstream anammox, although more studies are still needed for its full-scale demonstration and optimization (Ali et al. 2021, Zhang et al. 2019).
However, most of the shortcut N processes, such as PN/A or PDN/A, focused primarily on N removal via nitrite, and did not incorporate and/or optimize for biological phosphorus (P) removal, particularly enhanced biological P removal (EBPR). To date, only a few studies have attempted to incorporate mainstream PN/A with EBPR (Cao et al. 2016, Jia et al. 2023), which is usually due to the presumed incompatibility between these two processes (Roots et al. 2020, Yang et al. 2017). For instance, a favorable carbon-to-phosphorus (C/P) ratio and alternating anaerobic-aerobic/anoxic conditions are required for the enrichment of polyphosphate (polyP) accumulating organisms (PAOs) or denitrifying PAOs (DPAOs) performing EBPR, which are hardly satisfied in the B-stage PN/A process with limited influent organics from A-stage and deficiency/absence of deep anaerobic zone. The supplemental carbon sources added to the B-stage influent for EBPR, however, would potentially compromise NOB suppression, undermining the stability of shortcut N removal performance (Jia et al. 2023, Li et al. 2020b). Recently, the successful implementation of side-stream EBPR (S2EBPR) processes, which involve a side-stream anaerobic reactor for simultaneous PAO enrichment and volatile fatty acid (VFA) generation via hydrolysis and fermentation of return activated sludge (RAS) and/or mixed liquor, has promised the potential to enable stable EBPR performance with more flexible configurations without direct dependence on influent carbon (Barnard et al. 2017, Gu et al. 2019, Onnis-Hayden et al. 2020, Wang et al. 2019). Compared to conventional EBPR, S2EBPR offers unique conditions that favor PAOs over GAOs, which holds the key to EBPR performance and stability, by introducing a more complex composition of VFAs (i.e., a mixture of acetate and propionate) that leads to higher PAO diversity (Onnis-Hayden et al. 2020, Srinivasan et al. 2021, Wang et al. 2019), and consequently, different P release/uptake kinetics. In addition, extended anaerobic conditions in S2EBPR further give PAOs competitive advantages due to their distinct maintenance and delayed decay associated with the versatility of metabolic pathways and the possession of multiple intracellular polymers (i.e., PHA, glycogen, and polyP) (Li et al. 2020a). Furthermore, S2EBPR allows for more flexible carbon supplement and improved carbon utilization efficiency, and thus EBPR can be achieved for systems under previously considered unfavorable C/P ratios (Sabba et al. 2023, Wang et al. 2019). Therefore, we hypothesized that the S2EBPR process could potentially be incorporated into the A/B process for achieving shortcut N removal and efficient P removal/recovery simultaneously.
Recently, a full-scale pilot testing with the incorporation of S2EBPR into an A/B process, defined as an integrated system combining A-stage high-rate activated sludge (HRAS) with B-stage partial (de)nitrification/anammox and side-stream EBPR (i.e., HRAS-P(D)N/A-S2EBPR), was conducted at the Chesapeake-Elizabeth Sewage Treatment Plant (Virginia Beach, VA, USA) to assess the feasibility of this newly established system (
Figure 1). To achieve nitrite accumulation, a series of continuous stirred-tank reactors (CSTRs) were intermittently aerated at the same interval and duration via ammonia versus NO
x–-N (nitrite+nitrate) (AvN) control strategy (Regmi et al. 2014). To enhance P removal, the S2EBPR reactor was operated in a side-stream RAS fermentation with supplemental carbon (SSRC) configuration (Wang et al. 2019), with A-stage wasted activated sludge (WAS) fermentate as a supplement. To the best of our knowledge, this was the first demonstration and comprehensive investigation of a pilot-scale HRAS-P(D)N/A-S2EBPR system for mainstream shortcut N removal and side-stream EBPR treating real mainstream municipal wastewater. The unique configuration and operating conditions are expected to exert distinct selection forces from conventional EBPR on the microbial ecology of the system. However, several fundamental issues related to processing design and operation, as well as the underlying mechanisms, remain unresolved regarding the coupling of A/B process with S2EBPR, including 1) How would operating conditions in the S2EBPR reactor (e.g., sludge retention time (SRT), hydraulic retention time (HRT), oxidation-reduction potential (ORP), etc.), which differ from the conventional EBPR, affect AOB/NOB growth and activity? 2) What effects would intermittently aerated CSTRs at varying dissolved oxygen (DO) concentrations have on PAOs/DPAOs and their competitors, glycogen accumulating organisms (GAOs)/denitrifying GAOs (DGAOs)? 3) Is there any effect on the functionally relevant microorganisms, particularly PAOs/DPAOs, GAOs/DGAOs, AOB, NOB, and denitrifying bacteria (DNB), when the A-stage WAS fermentate was introduced into the S2EBPR reactor as the additional carbon load?
Therefore, this study aimed to comprehensively and systematically evaluate the feasibility, microbial ecology, and relevant mechanisms involved in the incorporation of S2EBPR into the novel HRAS-P(D)N/A-S2EBPR system. In addition to routine monitoring of influent and effluent water quality, mass balances of carbon, N, and P, metabolic activities batch tests, and microbial community analyses via 16S rRNA gene amplicon sequencing and single-cell Raman micro-spectroscopy (SCRS) were performed. The outcome will improve our fundamental understanding of this unique HRAS-P(D)N/A-S2EBPR process, the associated microbial community structure and phenotypes, and elucidate the key factors and mechanisms that govern process optimization.
2. Materials and methods
2.1. Pilot plant setup and operating conditions
The pilot plant is located at the Chesapeake-Elizabeth Sewage Treatment Plant (Virginia Beach, VA, USA) operated by Hampton Roads Sanitation District (HRSD). Details of this pilot plant were reported by Klaus (2019). The plant was operated in an A/B process (
Figure 1), receiving pre-screened and temperature-controlled municipal wastewater at 20℃ and a flow rate of 3.2 m
3/d (Printz 2019). The A-stage consists of three HRAS reactors in series with a total volume of 511 L, followed by an intermediate clarifier. The supernatant of the A-stage clarifier flows into the B-stage as its influent. The B-stage comprises a 53 L anaerobic reactor and four CSTRs in series with a total volume of 600 L, operated under AvN intermittent aeration control. A 174 L side-stream anaerobic reactor (S2EBPR) was introduced in B-stage. A fraction of RAS (~20–30%) from the B-stage clarifier was split into S2EBPR for sludge fermentation and PAO enrichment, with A-stage WAS fermentate as supplemental carbon. The S2EBPR is intermittently mixed for a duration of 1 min every 3 h, yielding an estimated side-stream SRT (SRT
SS) of ~7–10 h and a side-stream HRT (HRT
SS) of ~3.7 h. The effluent from the B-stage clarifier then passes through a 340 L anammox moving bed biofilm reactor (MBBR) loaded with K3 biofilm carriers (AnoxKaldnes, Sweden) for accomplishing shortcut N removal. The whole process is named HRAS-P(D)N/A-S2EBPR, while the focused processes in this study, including mainstream anaerobic reactor and CSTRs for shortcut N removal via PN and/or PDN, coupled with the side-stream anaerobic reactor for EBPR, is defined as the P(D)N-S2EBPR system.
The pilot testing period for this study was conducted from June 13th, 2019 to March 26th, 2020. The long-term operation was divided into three phases based on the major operational differences in total carbon loading level and corresponding performance: Phase 1 (start-up at low carbon load, Day 1–68), Phase 2 (nitrite accumulation at medium carbon load, Day 69–194), and Phase 3 (nitrite accumulation and EBPR at high carbon load, Day 195–288). Detailed influent and operating parameters of the system were summarized in
Table S1. The performance of the system during the pilot testing period was routinely monitored by measuring chemical oxygen demand (COD), soluble COD (sCOD), total phosphorus (TP), orthophosphate (PO
43–-P), total Kjeldahl nitrogen (TKN), nitrate-nitrogen (NO
3–-N), nitrite-nitrogen (NO
2–-N), ammonium-nitrogen (NH
4+-N), etc. in the influent, effluent, and each zone of the system. The sampling parameters, sites, frequency, and dates are summarized in
Table S2. The nitrite accumulation ratio (NAR) was calculated by the ratio of effluent NO
2–-N to the sum of effluent NO
2–-N and NO
3–-N.
2.2. Metabolic activity batch tests
To evaluate the EBPR activity, anaerobic-aerobic batch tests were conducted on-site weekly following the previously described protocol (Gu et al. 2008). Briefly, fresh activated sludge withdrawn from the 4th aerobic tank at B-stage (CSTR-4) was initially aerated for 1 h to remove residual organics. Sodium acetate (HAc) of ~100 mg COD/L was then added and the sludge was continuously purged with nitrogen gas to maintain an anaerobic condition (DO < 0.01 mg/L). After 2 h of anaerobic conditions, the sludge was bubbled with air to obtain an aerobic condition (DO > 2.0 mg/L) for 3 h. The pH and temperature in the batch test were maintained at 7.0 and 20°C, respectively. Samples were collected every 30 min, filtered by 0.45 μm filter, and then analyzed for sCOD, PO43–-P, NH4+-N, NO2–-N, and NO3–-N, etc. To investigate the effect of different carbon sources on PAO activity, activated sludge samples and A-stage WAS fermentate (fermentate) in the P(D)N-S2EBPR system were collected at the end of the pilot testing and immediately shipped to Cornell lab. Either fermentate or HAc was utilized as the carbon source and added at the beginning of the anaerobic-aerobic batch test. Simultaneously, the presence of denitrifying PAOs (DPAOs) was evaluated by carrying out an anaerobic-anoxic batch test with fermentate and NO3–-N as carbon source and electron acceptor, respectively. DPAO activity and fraction of DPAOs over total PAOs were calculated according to previous literature (Ferrentino et al. 2018, Hu et al. 2003).
To evaluate the activities of nitrifying bacteria (i.e., AOB and NOB), ex-situ nitrification batch tests were routinely carried out. Briefly, ~4 L pre-aerated activated sludge from CSTR-4 was spiked with ammonium chloride at a concentration of 25 mg N/L. The pH and DO in the batch test were maintained at ~7.0 and ~3–5 mg/L, respectively. To evaluate the activities of ordinary denitrifiers (i.e., DNB), ~2 L RAS from B-stage clarifier was mixed with A-stage effluent of ~100–150 mg COD/L and NO3–-N of 20 mg N/L. Samples were taken every 15 min for a total of 1 h, and then analyzed for NH4+-N, NO2–-N, NO3–-N, and sCOD concentrations. AOB, NOB, and DNB activities were calculated by linear regression as the NOx–-N production rate, NO3–-N production rate, and NO3–-N reduction rate, respectively.
2.3. Microbial community analysis
The B-stage microbial community was periodically analyzed via 16S rRNA gene amplicon sequencing following the method described by Wang et al. (2019). Briefly, genomic DNA was extracted using the DNeasy PowerSoil Pro Kit (Qiagen, Germantown, MD, USA). The extracted DNA was amplified targeting the V4 region of the 16S rRNA gene using the primers 515F (5'-GTGCCAGCMGCCGCGGTAA-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3') and the amplicons were sequenced on the Illumina MiSeq platform using the paired-end (2 × 250 bp) option in the Microbial Analyses, Resources, and Service (MARS) facility at the University of Connecticut (Storrs, CT, USA). Raw reads were preprocessed following the standard operating procedure by Mothur software (Schloss et al. 2009). High-quality reads were then clustered at a 97% similarity to obtain the operational taxonomic units (OTUs) and annotated based on the Silva database (Quast et al. 2012). Non-metric multidimensional scaling (NMDS) with distance bray was plotted for the microbial community analysis using the meta MDS function from the vegan package (R Core Team 2013).
2.4. Single-cell Raman micro-spectroscopy-based phenotyping
To quantify PAOs and GAOs based on their phenotypic characteristics in the microbial community, activated sludge samples were collected at the start of anaerobic phase (t = 0 min), end of anaerobic phase (t = 120 min), and end of aerobic phase (t = 300 min) of the anaerobic-aerobic batch tests and prepared for single-cell Raman micro-spectroscopy (SCRS) analysis, as described previously (Majed et al. 2012, Wang et al. 2021). The detailed SCRS analysis protocol related to sample pretreatment, data acquisition, and processing was described in the Supplementary Information.
2.5. Chemical analyses
The concentrations of COD, PO43–-P, NO3–-N, NO2–-N, NH4+-N, mixed liquor suspended solids (MLSS), and mixed liquor volatile suspended solids (MLVSS) were analyzed according to the Standard Methods (APHA. 2012) or USEPA methods or using HACH kits (HACH, Loveland, CO, USA). VFAs, including acetic, propionic, n-butyric, isobutyric, n-valeric, isovaleric, and caproic acids, were analyzed using high-performance liquid chromatography (HPLC). Glycogen was extracted (0.9M HCl and 3 h of digestion time) and determined using a liquid chromatography-mass spectrometer (LC-MS) (Thermo Scientific, Waltham, MA, USA). Poly-β-hydroxyalkanoates (PHAs), including poly-β-hydroxybutyrate (PHB), poly-β-hydroxyvalerate (PHV), and poly-β-hydroxy-2-methylvalerate (PH2MV), were extracted (3% sulfuric acid and 3 h of digestion time) and determined by a gas chromatography-mass spectrometry (GC-MS) (Agilent, Santa Clara, CA, USA) (Lanham et al. 2013).
2.6. Statistics
All the tests were performed at least in duplicate. The one-way analysis of variance (ANOVA) and the least significant difference (LSD) post hoc tests were used to test the significance among different samples using Matlab R2014a (Mathworks, Natick, MA, USA).
Figure 1.
Schematic flow diagrams of the pilot plant run in a HRAS-P(D)N/A-S2EBPR configuration and the P(D)N-S2EBPR system (marked by the red dashed box) at the Chesapeake-Elizabeth Sewage Treatment Plant (Virginia Beach, VA, USA). HRAS: high-rate activated sludge; CSTRs: continuous stirred-tank reactors; S2EBPR: side-stream biological P remover; RAS: return activated sludge; WAS: waste activated sludge; MBBR: moving bed biofilm reactor.
Figure 1.
Schematic flow diagrams of the pilot plant run in a HRAS-P(D)N/A-S2EBPR configuration and the P(D)N-S2EBPR system (marked by the red dashed box) at the Chesapeake-Elizabeth Sewage Treatment Plant (Virginia Beach, VA, USA). HRAS: high-rate activated sludge; CSTRs: continuous stirred-tank reactors; S2EBPR: side-stream biological P remover; RAS: return activated sludge; WAS: waste activated sludge; MBBR: moving bed biofilm reactor.
Figure 2.
Biological P removal performance during the pilot testing: (A) influent PO43–-P, effluent PO43–-P, and PO43–-P removal efficiency; and (B) cumulative relative frequency distribution for effluent PO43–-P during different phases. The effluent PO43–-P data shown in the summary table are in mg/L.
Figure 2.
Biological P removal performance during the pilot testing: (A) influent PO43–-P, effluent PO43–-P, and PO43–-P removal efficiency; and (B) cumulative relative frequency distribution for effluent PO43–-P during different phases. The effluent PO43–-P data shown in the summary table are in mg/L.
Figure 3.
Biological N removal performance during the pilot testing: influent TIN, effluent TIN, and TIN removal efficiency. TIN: total inorganic nitrogen (sum of NO3–-N, NO2–-N and NH4+-N).
Figure 3.
Biological N removal performance during the pilot testing: influent TIN, effluent TIN, and TIN removal efficiency. TIN: total inorganic nitrogen (sum of NO3–-N, NO2–-N and NH4+-N).
Figure 4.
Fate of PO43–-P in (A) Phase 1, (B) Phase 2, and (C) Phase 3. Value with a plus or minus sign represents the PO43–-P release and uptake, respectively. Value in brackets represents the percentage of PO43–-P in each unit compared to the total input. AN: anaerobic reactor; AE: intermittently aerated continuous stirred-tank reactors; S2EBPR: side-stream biological P remover; RAS: return activated sludge; WAS: calculated P removal via waste activated sludge. Units are in g/d.
Figure 4.
Fate of PO43–-P in (A) Phase 1, (B) Phase 2, and (C) Phase 3. Value with a plus or minus sign represents the PO43–-P release and uptake, respectively. Value in brackets represents the percentage of PO43–-P in each unit compared to the total input. AN: anaerobic reactor; AE: intermittently aerated continuous stirred-tank reactors; S2EBPR: side-stream biological P remover; RAS: return activated sludge; WAS: calculated P removal via waste activated sludge. Units are in g/d.
Figure 5.
The 16S rRNA gene amplicon sequencing-determined relative abundance of key functionally relevant microorganisms during the pilot testing: (A) known PAOs (i.e., Accumulibacter, Tetrasphaera, and Dechloromonas) and total known GAOs; (B) Nitrosomonas AOB, Nitrospira NOB, and known denitrifiers.
Figure 5.
The 16S rRNA gene amplicon sequencing-determined relative abundance of key functionally relevant microorganisms during the pilot testing: (A) known PAOs (i.e., Accumulibacter, Tetrasphaera, and Dechloromonas) and total known GAOs; (B) Nitrosomonas AOB, Nitrospira NOB, and known denitrifiers.
Figure 6.
The single-cell Raman spectroscopy-determined EBPR phenotypes within a P release and uptake batch test for each phase (error bars represent the standard error): (A) relative abundance of PAOs, GAOs and PPOs; (B) polyP, PHA and glycogen contents in PAOs; and (D) PHA content in PAOs, GAOs and PPOs. PPOs: PHA-producing organisms.
Figure 6.
The single-cell Raman spectroscopy-determined EBPR phenotypes within a P release and uptake batch test for each phase (error bars represent the standard error): (A) relative abundance of PAOs, GAOs and PPOs; (B) polyP, PHA and glycogen contents in PAOs; and (D) PHA content in PAOs, GAOs and PPOs. PPOs: PHA-producing organisms.
Table 1.
Specific kinetic rates and stoichiometric ratios observed in the ex situ P release and uptake batch tests during the pilot testing period, and comparison with other EBPR systems reported in the literature.
Table 1.
Specific kinetic rates and stoichiometric ratios observed in the ex situ P release and uptake batch tests during the pilot testing period, and comparison with other EBPR systems reported in the literature.
Parameter |
Phase 1 |
Phase 2 |
Phase 3 |
S2EBPR systemsa
|
EBPR systemsb
|
Prel [mg P/(g VSS·h)] |
7.7±2.6 |
7.2±2.4 |
18.8±7.4 |
2.9–7.0 |
3.2–31.9 |
Pup [mg P/(g VSS·h)] |
3.9±1.1 |
3.2±1.3 |
8.3±6.0 |
0.6–2.6 |
0.7–19.2 |
Pup/Prel
|
0.5±0.1 |
0.4±0.1 |
0.4±1.1 |
0.2–0.5 |
0.2–0.7 |
HAcup [mg HAc/(g VSS·h)] |
14.4±4.1 |
10.3±3.5 |
29.1±11.1 |
7.7–24.9 |
9.0–47.0 |
P/HAc (P-mol/C-mol) |
0.5±0.1 |
0.7±0.1 |
0.6±0.1 |
0.2–1.3 |
0.1–1.1 |
Table 2.
Specific kinetic rates and stoichiometric ratios observed in the ex situ P release and uptake batch tests using different carbon sources and electron acceptors with the activated sludge collected at the end of Phase 3, and comparison with other full-scale EBPR systems reported in the literature.
Table 2.
Specific kinetic rates and stoichiometric ratios observed in the ex situ P release and uptake batch tests using different carbon sources and electron acceptors with the activated sludge collected at the end of Phase 3, and comparison with other full-scale EBPR systems reported in the literature.
Carbon source |
This study |
Previous EBPR studies |
Fermentate |
HAc |
Real wastewater |
HAc |
Anaerobic results |
|
|
|
|
P release (mg P/g VSS) |
10.9 |
8.6 |
3.4–16.7a
|
5.1–24.3b
|
Prel [mg P/(g VSS·h)] |
6.6 |
7.8 |
– |
3.2–31.9c
|
VFAup [mg VFA/(g VSS·h)] |
38.2 |
37.9 |
– |
9.0–47.0d,e
|
P/VFA (P-mol/C-mol) |
0.16 |
0.25 |
0.63–1.00a
|
0.11–1.30c,f
|
PHA/VFA (C-mol/C-mol) |
0.21 |
0.23 |
1.20–1.39a
|
0.63–2.10a,c
|
PHV/PHB (C-mol/C-mol) |
3.16 |
0.11 |
0.63–0.86a
|
0.00–0.26g,h
|
Anoxic results |
|
|
|
|
P uptake (mg P/g VSS) |
2.5 |
– |
– |
– |
Pup [mg P/(g VSS·h)] |
1.2 |
– |
– |
0.0–5.9d
|
NRR [mg N/(g VSS·h)] |
1.3 |
– |
– |
1.0–10.0i,j
|
P/PHA (P-mol/C-mol) |
0.37 |
– |
– |
0.46e
|
Aerobic results |
|
|
|
|
P uptake (mg P/g VSS) |
8.6 |
7.7 |
3.6–18.2a
|
4.1–15.0b
|
Pup [mg P/(g VSS·h)] |
4.4 |
3.8 |
– |
0.7–19.2c,d
|
P/PHA (P-mol/C-mol) |
0.75 |
0.82 |
0.71–0.90a
|
0.20–3.68c,k
|