Skip to main content

Effect of COD:N ratio on biological nitrogen removal using full-scale step-feed in municipal wastewater treatment plants


This study investigated the effect of low and high chemical oxygen demand (COD):N ratios on biological nitrogen removal and microbial distributions in full-scale step-feed (SF) municipal wastewater treatment plants (WWTPs) in Thailand (SF1) and Taiwan (SF2). The SF1 WWTP had a low COD:N (4:1) ratio, a long solids retention time (SRT) (> 60 d), and low dissolved oxygen (DO) conditions (0.2 mg L− 1 in anoxic tank and 0.9 mg L− 1 in aerobic tank). The total nitrogen (TN) removal efficiency was 48%. The SF2 WWTP had a high COD:N (10:1) ratio, a short SRT (7 d), and high DO (0.6 mg L− 1 in anoxic tank and 1.8 mg L− 1 in aerobic tank). The TN removal efficiency was 61%. The nitrification and denitrification rates from these two plants were inadequate. Using a quantitative polymerase chain reaction (qPCR) technique, the populations of ammonium oxidizing bacteria (AOB) and ammonium oxidizing archaea were quantified. Measurement of ammonia monooxygenase (amoA) gene abundances identified these AOB: Nitrosomonas sp., Nitrosospira sp., Nitrosoccus sp. and Zoogloea sp. Higher amounts of the archaeal-amoA gene were found with long SRT, lower DO and COD:N ratios. Abundance of Nitrobacter sp. was slightly higher than Nitrospira sp. at the SF1, while abundance of Nitrobacter sp. was two orders of magnitude greater than Nitrospira sp. at the SF2. More denitrifying bacteria were of the nirS-type than the nirK-type, especially at higher COD:N ratio. Most bacteria belong to the phyla Acidobacteria, Actinobacteria Bacteroidetes, Chloroflexi, Proteobacteria. The results from this work showed that insufficient carbon sources at the SF1 and high DO concentration in anoxic tank of SF2 adversely affected nitrogen removal efficiencies. In further research work, advanced techniques on the next generation sequencing with different variable regions should be recommended in full-scale WWTPs.


Increases in water pollution are usually related to growing urban populations. Efficient removal of nitrogen in wastewater treatment plants (WWTPs) is essential to avoid downstream eutrophication which adversely affects not only animal but also human health globally. Nitrogen is removed from wastewaters with physical methods (air stripping), chemical methods (ion exchange), biological treatment (nitrification and denitrification processes), and/or combinations of these. Biological treatment processes are dominant over all other physical and chemical methods and are attractive because of relative low costs [1].

The most popular domestic wastewater treatment system for large communities is activated sludge process with plug flow configuration. However, with some site-specific conditions, existing processes or equipment and demand for high biological nitrogen removal efficiency, a modification of plug flow with step-feed is recommended. Dividing a reactor tank into anoxic and aerobic zones and/or using step-feed configuration are commonly recommended for improving nitrogen removal [2]. However, not all step-feed configurations require pre-anoxic process. The step-feed process has many advantages over conventional activated sludge processes, including more uniform distribution of oxygen demand, superior ability to handle peak wet–weather flows, and flexible operation. Step-feed systems can often achieve treatment objectives with smaller bioreactor volumes [3], and the process will often achieve low effluent total inorganic nitrogen concentrations [4].

The key factors affecting full-scale step-feed WWTPs are dissolved oxygen (DO), solids retention time (SRT), and hydraulic retention time (HRT). Operation with internal recycle reduced total nitrogen (TN) concentrations between 5 and 8 mg L− 1 with SRT between 3 and 15 d, total HRT 3–5 h, anoxic zone time 0.5–1.5 h, aerobic zone time 2.5–3.5 h, aerobic tank DO 2–3 mg L− 1 and anoxic tank DO ≤0.2 mg L− 1 [5]. Wang and Chen [6] reported TN removal efficiency of > 64% in a full-scale step-feed system operated at DO of < 0.25 mg L− 1 in anoxic tank and 2 mg L− 1 in aerobic tank with temperature between 11.9–23.3 °C. Moreover, in a pilot-scale step-feed system investigated by Ge et al. [7], TN removal efficiencies varied from 75 to 86% with different SRTs at 10–15 d (75%), 10–12 d (82%) and 8–10 d (86%) under DO concentrations of 1.2–2.0 mg L− 1 in aerobic zones.

Information on TN removal for full-scale step-feed municipal WWTP specifically for low and high chemical oxygen demand (COD):N ratios and various DO concentrations is rare in the literature. For this reason, our work was focused on two full-scale step-feed WWTPs in Bangkok, Thailand and Taipei, Taiwan. These two full- scale WWTPs with similar configurations were selected because of low and high COD:N ratios in influent. The definitions for high and low COD:N ratios of wastewaters are > 4.3:1 and ≤ 4.3:1, respectively. The study compared the efficiencies of nitrogen removal from these step-feed WWTPs, and different observations due to design parameters and operating conditions were explained. In addition, the abundance of microbial communities in these full-scale WWTPs were investigated and discussed. The results from this work could be applied to step-feed WWTPs in either country to solve carbon limitation when treating low COD:N wastewater and/or reduce aeration energy by using low-DO processes for improving biological nitrogen removal efficiencies.

Materials and methods

Wastewater treatment systems

Two underground full-scale municipal step-feed WWTPs were selected from the downtown area of two capital cities, Bangkok, Thailand, (SF1) and Taipei, Taiwan, (SF2). All wastewaters samples were collected and analyzed over an entire year (2018–2019). Both plants were designed for removal of both organic matter and nitrogen with reaction tanks consisting of anoxic-aerobic zones. The SF1 (see Fig. 1a) had four feed points to four anoxic and four aerobic tanks. However, due to low flow conditions into the SF1 system, only two feed points were operated and rotated with another two feed points. The SF2 (see Fig. 1b) had three feed points to three anoxic tanks and three large aeration tanks (each large aeration tank was divided to four small aeration tanks) due to high flow of the system. These two WWTPs were built underground because of land limitations in these dense capital cities. The above ground areas of these two plants were used as recreation and education centers.

Fig. 1

Diagrams of the full-scale step feed systems in a Thailand and b Taiwan

Analytical methods

Influent and effluent samples from each full-scale step-feed WWTP in this work were collected monthly during 2018–2019. Characteristics of these samples were measured by using the method described in Standard Methods for the Examination of Water and Wastewater (2005). Mixed liquor suspended solids (MLSS) and mixed liquor volatile suspended solids (MLVSS) from each reactor tank were also analyzed by following the method described in Standard Methods. Temperature and pH were immediately measured in the field.

Nitrification and denitrification rates

To calculate the nitrification rates with various COD:N ratios, the concentration of ammonium-nitrogen (NH4+-N) in influent and effluent was determined. The nitrification rate was defined based on the NH4+-N removal as shown in Eq. (1) [8],

$$ {\upgamma}_{\mathrm{nitrification}}={\mathrm{Q}}_{\mathrm{in}}\times \frac{\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{in}\mathrm{fluent}}-\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{effluent}}}{{\mathrm{V}}_{\mathrm{reactor}}\times \mathrm{VS}{\mathrm{S}}_{\mathrm{nitrifying}\ \mathrm{in}\ \mathrm{reactor}}}, $$

where γnitrification is the nitrification rate (d− 1), Qin is flow rate (m3 d− 1), Vreactors is volume of reactors (m3), and VSSnitrifying in reactor is MLVSS of nitrifying organisms in reactor (mg L− 1).

The nitrifying organisms in the reactor is calculated based on MLVSS using Eqs. (2) and (3),

$$ \mathrm{MLVS}{\mathrm{S}}_{\mathrm{nitrifying}\ \mathrm{in}\ \mathrm{reactor}}={\mathrm{f}}_{\mathrm{N}}\times \mathrm{VS}{\mathrm{S}}_{\mathrm{reactors}}, $$
$$ {\mathrm{f}}_{\mathrm{N}}=\frac{0.16\times \left(\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{influent}}-\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{effluent}}\right)}{0.6\times \left(\mathrm{BO}{\mathrm{D}}_{\mathrm{influent}}-\mathrm{BO}{\mathrm{D}}_{\mathrm{effluent}}\right)+0.16\times \left(\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{influent}}-\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{effluent}}\right)}, $$

where fN is the fraction of nitrifying organisms presenting in the mixed liquor of a step-feed system. This fraction of nitrifying organisms can be estimated using Eq. (3). BODinfluent and BODeffluent are the concentrations of biochemical oxygen demand in influents and effluents (mg L− 1), respectively. The denitrification rate was defined as Eq. (4),

$$ {\upgamma}_{\mathrm{denitrification}}=\frac{{\mathrm{Q}}_{\mathrm{in}}\times \left[\left({\mathrm{TKN}}_{\mathrm{in}\mathrm{fluent}}-\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{effluent}}\left)-\right(\mathrm{N}{\mathrm{O}}_3^{-}-{\mathrm{N}}_{\mathrm{effluent}}\right)\right]}{\left({\mathrm{V}}_{\mathrm{reactors}}\right)\times \left(\mathrm{VS}{\mathrm{S}}_{\mathrm{reactors}}\right)}, $$

where γdenitrification is the denitrification rate (d− 1) and [TKNinfluent – (NH4+-N)influent] can be substituted with organic nitrogen in the influent. The TN removal (%) was calculated using Eq. (5).

$$ {\mathrm{TN}}_{\mathrm{removal}}\left(\%\right)=\frac{\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{influent}}-\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{effluent}}}{\mathrm{N}{\mathrm{H}}_4^{+}-{\mathrm{N}}_{\mathrm{influent}}}\times 100, $$

Microbial communities analysis

Sludge samples from SF1 and SF2 and were taken from both anoxic and aerobic tanks for the analysis of microbial communities. The nitrifying bacterial communities were identified through analysis of ammonia monooxygenase (amoA) gene abundances of ammonium oxidizing bacteria (AOB) and ammonium oxidizing archaea (AOA). The 16S rDNA target gene of Nitrospira (NSR) and Nitrobacter (Nitro) was used to determine nitrite oxidizing bacteria (NOB) abundance. The functional targeted gene of nirK and nirS genes were used as molecular markers for denitrifying bacterial (DNB) abundances.

DNA extraction and polymerase chain reaction (PCR) amplification

1 mL of the samples were taken for DNA extraction following the manufacturer’s method using FavorPrep™ soil DNA isolation mini kit (Favogen® Biotech Corp, Taiwan). The PCR protocol and oligonucleotide primers for quantitative PCR (qPCR) and denaturing gradient gel electrophoresis are shown in Table S1 in Supplemental Materials.

qPCR of functional and 16S rDNA genes

The qPCR mixture contained 10 μL of SYBR Green by Luna® Universal qPCR Master Mix (New England Biolabs, MA, USA), 10 pmol of each primer, 1 μL of DNA template (~ 10–20 ng μL− 1) and nuclease-free water up 20 μL per reaction. Each sample and the standard series (100–108) were performed in triplicate on CFX96 Touch™ Real-Time PCR detection systems (Bio-Rad Laboratories, CA, USA) and the results were active based on the correlation coefficient of the standard curve (R2 = 0.995).

Denaturing gradient gel electrophoresis (DGGE) fingerprints

The PCR mixture contained 10X Ex Taq™ buffer, 5 units μL− 1 TaKaRa Ex Taq™, 2.5 mM dNTP Mixture, 10 pmol of each primer, 1 μL of DNA template (~ 10–20 ng μL− 1) and nuclease–free water up to 25 μL per reaction. Each sample was completed on T100™ Thermal cycler (Bio-Rad Laboratories, CA and USA). 15 μL of each PCR product was loaded into individual lanes of a DGGE gel of 8% (W/V) acrylamide gel with 35–60% (EUB) and 35–50% (AOB), and 6% (W/V) acrylamide gel with 20–50% (AOA) denaturing gradients. Electrophoresis was performed for 16 h at 58 °C with a constant voltage at 80 V in 1X TAE buffer. Each DGGE band was excised with a scalpel, DNA fragment was eluted from the band by milli-Q water overnight in a refrigerator, followed by PCR with the same primer without attached CG-clamp. Representative sequences were aligned against the National Center for Biotechnology Information database using Basic Local Alignment Search Tool. In this work, nirK and nirS genes were used because they are typically contained in denitrifying bacteria, but are structurally different from nitrite reductase (distinguish the copper-dependent nitrite reductase and cytochrome cd1-containing nitrite reductase).

Results and discussion

Key operating conditions of the full-scale step-feed WWTPs and wastewater quality

The key average operating parameters of the SF1 and SF2 WWTPs are shown in Table 1. The characteristics of influent and effluent of each SF are shown in Table 2. The BOD:TN and COD:N ratios of SF1 were 2:1 and 4:1, respectively. The BOD:TN and COD:N ratios of SF2 were 4.6:1 and 10:1, respectively. The wastewater treatment loading rate, BOD:TN and COD:N ratios of SF1 were significantly lower than those in SF2. Both SF1 and SF2 were able to remove SS, COD, and BOD well, but not nitrogen and phosphorus. At SF1, there was not enough carbon source (low BOD in the influent) for denitrifying bacteria as electron donor. For this reason, this insufficient carbon source would affect on denitrification process.

Table 1 Key average parameters of two full-scale step feed (SF) WWTPs
Table 2 Characteristics of wastewaters from each full-scale step feed WWTP

The process performances on COD, NH4+-N, NO2- N, and NO3- N concentrations profile of SF1 and SF2 WWTPs are shown in Fig. 2. Although BOD:N and COD:N ratios of SF1 were significantly lower than SF2, the NH4+-N removal efficiency of SF1 (> 88%) was higher than that of SF2 (59%). Average temperature of wastewater at SF1 was higher than average temperature of wastewater at SF2 (Table 2). Higher average temperature and longer SRT for SF1 could be significant factors promoting AOB activities. Although the DO concentration in the aerobic tank of SF1 was quite low, the DO was sufficient for adequate nitrification. It is also noted that the concentration of nitrifier communities at SF1 was significantly higher than that for SF2 (see subsequent Section microbial communities AOB and NOB populations and communities).

Fig. 2

Process performances on COD concentrations profile from a SF1 and from b SF2, and nitrogen (NH4+-N, NO2- N, and NO3-N) concentrations profile from c SF1 and from d SF2

Nitrification and denitrification rates

Overall biological nitrogen removal in SF1 and SF2 was determined by calculating nitrification and denitrification rates in the aerobic and anoxic tanks. These two rates should not be the same value [9]. In this work the nitrification and denitrification rates of SF1 and SF2 were significantly different (Table 3). There are several possible explanations for this inequality. First, Thai sewage piping combines wastewater and rainwater that occurs all seasons, diluting the Thai influent BOD and SS to very low levels (< 50 mg L− 1 and < 90 mg L− 1, respectively) [10]. Second, the MLVSS:MLSS ratio of SF1 was only 0.45–0.55 compared with that of SF2 (0.8–0.82) due to longer SRT in both anoxic and aerobic tanks in SF1. For this reason, when MLVSS is used to calculate biomass, inaccurate higher estimations of microorganisms would result. The significant difference of MLVSS:MLSS ratio between SF1 and SF2 might also be due to the absence of primary clarifier in SF1. The main purpose of a primary clarifier is to remove solids and particulates. Third, other factors affecting the growth of nitrifiers and denitrifiers would include DO concentration, SRT duration and temperature. Especially important was maintenance of appropriate DO concentration (< 0.2 mg L− 1) in the anoxic phase [6].

Table 3 Nitrification and denitrification rates and TN removal efficiencies of both SF1 and SF2

In this study, it was shown that the longer SRT (> 60 d) of SF1 promoted TN removal efficiency (48%) although the COD:N at this plant was quite low. Davies et al. [11] reported that longer SRTs improved nitrification and denitrification, resulting in high TN removal efficiency. At SF2, the operation was normal with sufficient carbon presence (high COD:N ratio) but TN removal efficiency (only 61%) was not much better than that of SF1. Moreover, the nitrification and denitrification rates were only 1.23 and 0.12 g NH4+-N g− 1 VSS d− 1, respectively. The main reason for the low nitrogen removal efficiency and differing microbial processes was excessive DO in anoxic tank (0.6 mg L− 1). Other investigators [9, 12, 13] have stated that high DO concentration enhances nitrification rates in aerobic tank while not increasing denitrification rates in anoxic tank. Maintaining lower DO concentrations (< 1.0 mg L− 1) throughout an entire year in aerobic tank adversely affected the nitrification process (TN removal efficiency only 55%). Meng et al. [14] reported that TN removal was increased to 78% by increasing DO concentrations to > 1.0 mg L− 1. Wang and Chen [6] demonstrated that a simultaneous nitrification-denitrification (anoxic-aerobic) process could result in TN removal efficiency of 57%. They also suggested that DO concentration in anoxic tank should be < 0.25 mg L− 1 and aeration should be reduced when the DO concentration exceeds 2 mg L− 1 in the aerobic tank. This current study showed that at both SF1 and SF2, DO levels between 0.2 and 0.6 mg L− 1 in anoxic tank could be postulated to impact on the denitrification rate.

Microbial communities by using qPCR and DGGE technique

AOB and archaea (AOA) populations and communities

The different COD:N ratios at SF1 and SF2 affected the sizes of archaeal-amoA (AOA) and bacterial-amoA (AOB) populations. Figure 3 shows that a large amount of AOA (1.0 × 105 copies mL− 1 sludge) was found with the low COD:N ratio of SF1. However, very low amounts of AOA (1.0 × 100 copies mL− 1 sludge) were found with the higher COD:N ratio of SF2. The very low amount of AOA found in SF2, but not in SF1, is due to the high DO concentration (> 1.8 mg L− 1) maintained in the aerobic tank. The quantitative results for AOA and AOB in this study were similar to that of Kayee et al. [15] who found an abundance of AOA in municipal full-scale anoxic and aerobic tanks at the Bangkok WWTP, which had low COD:N ratio (4.3:1), still higher than the COD:N ratio of SF1 in this work. In Kayee’s work, it was shown that the low DO concentration was maintained in aerobic tank of these WWTPs. Low DO concentration and longer SRT in full-scale WWTPs promoted large populations of AOA [16,17,18]. Gao et al. [16] found that a high ratio of COD:N (10.7:1) in full-scale activated sludge plants in Beijing, China, along with low DO concentration (0.5 mg L− 1) in the aerobic tank, led to large populations of AOA. It is noted that the relatively high temperature of wastewater in this work (only for SF1) could be a contributing factor for the abundance of AOA and AOB populations. Several studies reported the effects of warm climate on AOA and nitrifying community. For example, Limpiyakorn et al. [17] found high AOA and other nitrifying communities in domestic and industrial WWTPs in Thailand. Sinthusith et al. [18] reported that long SRT with high temperature (30 °C) and pH > 7 at the WWTP in Thailand was associated with the dominance of AOA amoA genes over AOB amoA genes.

Fig. 3

Microbial abundance of AOA, AOB, NOB, DNB and total bacteria (EUB) a at the SF1 with low COD:N and b the SF2 with high COD:N. Remark: BDL = below method detection limit

As indicated in Fig. 4, AOB species in the SF1 were the same as in SF2. AOB species included Nitrosomonas europaea, Nitrosomonas halophile, Nitrosospira multiformis, Nitrosospira tenuis and Zoogloea caeni via 16S rRNA of CTO primer pairs and Nitrosoccus halophilus via 16S rRNA of EUB primer pairs (see Table S2). Moreover, AOB communities present in this study were similar to those found by Shen et al. [19], who investigated the microbial community in a full-scale domestic WWTP (anoxic/oxic process). The main AOA communities at SF1 were Crenarchaeotal sp. and uncultured Thaumarchaeote. Thaumarchaeota are autotrophic and capable of performing the oxidation of NH4+ to NO2 [20, 21]. Generally, Crenarchaeotes have been found in extreme environments, such as low oxygen concentrations in aquatic systems, hot springs, and full-scale anaerobic digester systems [22,23,24]. For full scale WWTP applications, it would be advantageous to maintain conditions which support AOA and AOB communities to improve biological nitrogen removal.

Fig. 4

Phylogenetic tree of microbial communities at a the SF1 with low COD:N and b the SF2 with high COD:N

NOB populations and communities

In the second step of nitrification, Nitrobacter sp. and Nitrospira sp. are classically acknowledged as the most relevant NOB group in WWTPs. In Fig. 3, the copies number of Nitrospira via NSR gene were found 1.0 × 104 copies mL− 1 sludge at SF1 (low COD:N ratio), but they were present at less than 1.0 × 102 copies mL− 1 sludge at SF2 (high COD:N ratio). For Nitrobacter via Nitro gene 1.0 × 105 copies mL− 1 sludge) no significant difference was found between SF1 and SF2. Yu et al. [25] reported on fluorescence in situ hybridization (FISH) technique results in their submerged membrane bioreactors under two SRTs (30 and 90 d). The fast-growing Nitrobacter sp. was the dominant species at SRT 30 d, while the slow-growing Nitrospira was dominant at SRT 90 d. The results from this work coincided with a study by Yu et al. The amount of Nitrospira that was found in SF1 with long SRT (> 60 d) was significantly higher than the amount of Nitrospira in SF2 with short SRT (7 d). Moreover, Huang et al. [26] studied distribution of NOB communities in the full-scale WWTP by controlling DO concentration. They found that Nitrospira was dominant when low DO (< 0.9 mg L− 1) concentration was controlled, while Nitrobacter increased when DO concentrations were increased at higher than 0.9 mg L− 1. This DO fact could explain why Nitrospira was dominant at SF1 (operated with DO concentrations of 0.2–0.9 mg L− 1) and Nitrobacter was dominant at SF2 (operated with DO concentrations of 0.6–1.8 mg L− 1. Consequently, low DO conditions and long SRT would be the major operating conditions that contributed to high Nitrospira population in WWTP. It should be noted that in this work other NOB communities were not analyzed because the qPCR technique for analysis of 16S RNA would only reveal Nitrospira and Nitrobacter.

Effects of COD:N ratios on populations of DNB

The abundance of the denitrifying bacteria in this work is shown in Fig. 3. The denitrifiers are found in both anoxic and aerobic tanks of SF1 and SF2. Two gene types (nirK and nirS) were used to characterize denitrifiers. The nirK-type bacteria in both SF1 and in SF2 were found to be 104 copies mL− 1 sludge in anoxic tanks and 105 copies mL− 1 sludge in aerobic tanks. For the nirS-type denitrifiers, in SF1 anoxic and aerobic tanks, the bacteria were present in the same order of magnitude (105). However, in SF2 the nirS-type denitrifiers existed at two orders of magnitude higher (106 in anoxic tank and 107 in aerobic tank) than the nirK-type denitrifiers (104 in anoxic tank and 105 in aerobic tank). The higher population of nirS-type denitrifiers is attributed to the high COD:N ratio SF2. This observation is similar to the results from Wang et al. [27]. They found the number of nirS-type denitrifiers (104 to 105) was higher than that of nirK-type denitrifiers (103 to 104) in two full-scale WWTPs (upflow anaerobic sludge blanket and anaerobic/aerobic). Geets et al. [28] reported that nirK-type denitrifiers (105) were lower than nirS-type denitrifiers (106) in the sludge of industrial influent of anaerobic digester. They also found that nirK-type denitrifiers (106) were lower than nirS-type denitrifiers (107) in the sludge of domestic wastewater influent from hospital wastewater. From the present work and the cited studies, it is recommended that the nirS-type denitrifier growth be encouraged in full-scale WWTPs.

The bacterial species of denitrifiers at SF1 are shown to be similar to those at SF2 by PCR-DGGE technique (Fig. 3). The identified denitrifiers belonged to phyla Acidobacteria, Actinobacteria Bacteroidetes, Chloroflexi, Proteobacteria (Table S2). These five main phyla from this study are consistent with the work of Shen et al. [19]. They found these same phyla in full-scale domestic WWTPs (step-feed and anoxic/aerobic process). Wang et al. [27] also investigated the denitrifier communities in full-scale anoxic/oxic reactors using different analytical technique (next generation sequencing (NGS) by pyrosequencing and Illumina high-throughput sequencing). For this reason, Wang et al. [27] were more specific than the results from this work, identifying Genus of Thauera, Paracoccus, Hyphomicrobium, Comamonas and Azoarcus. In further research work, new and/or advanced techniques on the NGS with different variable regions should be recommended in full-scale WWTPs in order to identify specific denitrifier groups. With additional information we should understand which denitrifier communities are most effective in overall biological nitrogen removal in WWTPs.


In earlier studies it has been shown that TN removal is improved with higher COD:N ratios, longer SRTs, and low DO concentrations in anoxic zones. This work demonstrated that each variable is important in order to achieve adequate treatment in full-scale step-feed WWTPs. SF1 had longer SRT (> 60 d) and low anoxic DO (average 0.2 mg L− 1), but its low COD:N ratio (4:1) substantially impeded the denitrification portion of treatment. The COD:N ratio (10:1) at SF2 would provide enough carbon for denitrification, but the average DO (0.6 mg L− 1) in the anoxic region was too high for complete TN removal.

The distributions of archaeal and bacterial communities are dependent on operating parameters of the WWTP. At SF1 with low COD:N ratio, low DO, long SRT and high temperature the microbial abundance of AOA was greater than that of AOB. However, at SF2 with opposite parameter values, the AOB was more abundant. The predominant AOB communities were Nitrosomonas sp., Nitrosospira sp., Zoogloea sp. and Nitrosoccus sp. Higher amounts of Nitrospira were present at lower COD:N ratios. Higher amounts of Nitrobacter were found with high DO concentrations and higher COD:N. This work shows that at high COD:N ratio, nirS-type denitrifiers are more prevalent than nirS-type denitrifiers. General microbial communities belonged to phyla Acidobacteria, Actinobacteria Bacteroidetes, Chloroflexi, Proteobacteria were identified.

The results from this work, longer SRT > 30 d could be suggested as possible practical solution on issue of solving carbon limitation when treating low COD:N wastewater. Another practical solution on issue of reducing aeration energy, DO concentration in aeration tank should be maintained from 0.9 to 1.4 mg L− 1.

Availability of data and materials

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


  1. 1.

    Peng YZ, Ge SJ. Enhanced nutrient removal in three types of step feeding process from municipal wastewater. Bioresour Technol. 2011;102:6405–13.

    Article  Google Scholar 

  2. 2.

    Riffat R. Fundamentals of wastewater treatment and engineering. London: CRC Press; 2012.

    Book  Google Scholar 

  3. 3.

    Johnson BR, Daigger GT, Crawford G, Wable MV, Goodwin S. Full-scale step- feed nutrient removal systems: a comparison between theory and reality. The 76th Annual Water Environment Federation Technical Exhibition and Conference. Los Angeles; 2003.

    Google Scholar 

  4. 4.

    Amad O. Step feed BNR process achieving TN < 4 mg L−1: a case study for WSSC – Piscataway WWTP, Prince George County, Maryland. The 76th Annual Water Environment Federation Technical Exhibition and Conference. Los Angeles; 2003.

    Google Scholar 

  5. 5.

    Tchobanoglous G, Burton FL, Stensel HD. Wastewater engineering: treatment and reuse. 4th Singapore: McGraw-Hill Higher Education; 2004.

    Google Scholar 

  6. 6.

    Wang QB, Chen QW. Simultaneous denitrification and denitrifying phosphorus removal in a full-scale anoxic-oxic process without internal recycle treating low strength wastewater. J Environ Sci-China. 2016;39:175–83.

    Article  Google Scholar 

  7. 7.

    Ge SJ, Peng YZ, Lu CC, Wang SY. Practical consideration for design and optimization of the step feed process. Front Env Sci Eng. 2013;7:135–42.

    Article  Google Scholar 

  8. 8.

    Carrera J, Vicent T, Lafuente J. Effect of influent COD/ N ratio on biological nitrogen removal (BNR) from high-strength ammonium industrial wastewater. Process Biochem. 2004;39:2035–41.

    Article  Google Scholar 

  9. 9.

    Chiu YC, Lee LL, Chang CN, Chao AC. Control of carbon and ammonium ratio for simultaneous nitrification and denitrification in a sequencing batch bioreactor. Int Biodeter Biodegr. 2007;59:1–7.

    Article  Google Scholar 

  10. 10.

    Noophan P, Paopuree P, Kanlayaras K, Sirivithayapakorn S, Techkarnjanaruk S. Nitrogen removal efficiency at centralized domestic wastewater treatment plants in Bangkok, Thailand. EnvironmentAsia. 2009;2:30–5.

    Google Scholar 

  11. 11.

    Davies WJ, Le MS, Heath CR. Intensified activated sludge process with submerged membrane microfiltration. Water Sci Technol. 1998;38:421–8.

    Article  Google Scholar 

  12. 12.

    Pochana K, Keller J. Study of factors affecting simultaneous nitrification and denitrification (SND). Water Sci Technol. 1999;39:61–8.

    Article  Google Scholar 

  13. 13.

    Zhao HW, Mavinic DS, Oldham WK, Koch FA. Controlling factors for simultaneous nitrification and denitrification in a two-stage intermittent aeration process treating domestic sewage. Water Res. 1999;33:961–70.

    Article  Google Scholar 

  14. 14.

    Meng QJ, Yang FL, Liu LF, Meng FG. Effects of COD/N ratio and DO concentration on simultaneous nitrification and denitrification in an airlift internal circulation membrane bioreactor. J Environ Sci-China. 2008;20:933–9.

    Article  Google Scholar 

  15. 15.

    Kayee P, Sonthiphand P, Rongsayamanont C, Limpiyakorn T. Archaeal amoA genes outnumber bacterial amoA genes in municipal wastewater treatment plants in Bangkok (Retracted article. See vol. 72, pg. 262, 2016). Microb Ecol. 2011;62:776–88.

    Article  Google Scholar 

  16. 16.

    Gao JF, Luo X, Wu GX, Li T, Peng YZ. Quantitative analyses of the composition and abundance of ammonia-oxidizing archaea and ammonia-oxidizing bacteria in eight full-scale biological wastewater treatment plants. Bioresour Technol. 2013;138:285–96.

    Article  Google Scholar 

  17. 17.

    Limpiyakorn T, Sonthiphand P, Rongsayamanont C, Polprasert C. Abundance of amoA genes of ammonia-oxidizing archaea and bacteria in activated sludge of full-scale wastewater treatment plants. Bioresour Technol. 2011;102:3694–701.

    Article  Google Scholar 

  18. 18.

    Sinthusith N, Terada A, Hahn M, Noophan P, Munakata-Marr J, Figueroa LA. Identification and quantification of bacteria and archaea responsible for ammonia oxidation in different activated sludge of full-scale wastewater treatment plants. J Environ Sci Heal A. 2015;50:169–75.

    Article  Google Scholar 

  19. 19.

    Shen YJ, Yang DH, Wu Y, Zhang H, Zhang XX. Operation mode of a step-feed anoxic/oxic process with distribution of carbon source from anaerobic zone on nutrient removal and microbial properties. Sci Rep-UK. 2019;9:1153.

    Article  Google Scholar 

  20. 20.

    Leininger S, Urich T, Schloter M, Schwark L, Qi J, Nicol GW, et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature. 2006;442:806–9.

    Article  Google Scholar 

  21. 21.

    Santoro AE, Francis CA, de Sieyes NR, Boehm AB. Shifts in the relative abundance of ammonia-oxidizing bacteria and archaea across physicochemical gradients in a subterranean estuary. Environ Microbiol. 2008;10:1068–79.

    Article  Google Scholar 

  22. 22.

    Treusch AH, Leininger S, Kletzin A, Schuster SC, Klenk HP, Schleper C. Novel genes for nitrite reductase and Amo-related proteins indicate a role of uncultivated mesophilic Crenarchaeota in nitrogen cycling. Environ Microbiol. 2005;7:1985–95.

    Article  Google Scholar 

  23. 23.

    Martens-Habbena W, Berube PM, Urakawa H, de la Torre JR, Stahl DA. Ammonia oxidation kinetics determine niche separation of nitrifying archaea and bacteria. Nature. 2009;461:976–9.

    Article  Google Scholar 

  24. 24.

    Chouari R, Guermazi S, Sghir A. Co-occurence of Crenarchaeota, Thermoplasmata and methanogens in anaerobic sludge digesters. World J Microb Biot. 2015;31:805–12.

    Article  Google Scholar 

  25. 25.

    Yu T, Qi R, Li D, Zhang Y, Yang M. Nitrifier characteristics in submerged membrane bioreactors under different sludge retention times. Water Res. 2010;44:2823–30.

    Article  Google Scholar 

  26. 26.

    Huang ZH, Gedalanga PB, Asvapathanagul P, Olson BH. Influence of physicochemical and operational parameters on Nitrobacter and Nitrospira communities in an aerobic activated sludge bioreactor. Water Res. 2010;44:4351–8.

    Article  Google Scholar 

  27. 27.

    Wang Z, Zhang XX, Lu X, Liu B, Li Y, Long C, et al. Abundance and diversity of bacterial nitrifiers and denitrifiers and their functional genes in tannery wastewater treatment plants revealed by high-throughput sequencing. PLoS One 2014;9;e113603.

    Article  Google Scholar 

  28. 28.

    Geets J, de Cooman M, Wittebolle L, Heylen K, Vanparys B, De Vos P, et al. Real-time PCR assay for the simultaneous quantification of nitrifying and denitrifying bacteria in activated sludge. Appl Microbiol Biot. 2007;75:211–21.

    Article  Google Scholar 

Download references


This work was supported by the Thailand Science Research and Innovation (TSRI), Kasetsart University Research and Development Institute (KURDI) and Faculty of Engineering, Kasetsart University, Bangkok, Thailand (under Postdoctoral fellowship no. Post Doc.62/08/ENV) for grants that supported this research. The authors also would like to thank John Elliott (Golden, Colorado) for great help and Nimaradee Boonapatcharoen at Pilot plant development and training institute to support instrument in molecular analysis.


This work was supported by the Thailand Science Research and Innovation (TSRI), Kasetsart University Research and Development Institute (KURDI) and Faculty of Engineering, Kasetsart University, Bangkok, Thailand (under Postdoctoral fellowship no. Post Doc.62/08/ENV).

Author information




All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pongsak Noophan.

Ethics declarations

Competing interests

The authors declare they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

: Figure S1. DGGE fingerprints of (a) total bacteria, (b) amoA-AOB and (c) Arch amoA-AOA genes. Table S1. Oligonucleotide primers for PCR amplification via quantitative polymerase chain reaction (qPCR) and denaturing gradient gel electrophoresis (DGGE) techniques. Table S2. AOA, AOB and total bacteria communities.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Phanwilai, S., Noophan, P., Li, CW. et al. Effect of COD:N ratio on biological nitrogen removal using full-scale step-feed in municipal wastewater treatment plants. Sustain Environ Res 30, 24 (2020).

Download citation


  • Biological nitrogen removal
  • COD:N
  • Full-scale
  • Step-feed