Modification of Anaerobic Digestion Model No. 1 for modeling anaerobic digestion of cattle manure with changing solids retention time

In this study, Anaerobic Digestion Model No.1 (ADM1) was modified to incorporate changes in biochemical parameters due to solids retention time (SRT) variations. Cattle manure (CM) and thermally hydrolyzed CM were selected for testing. Continuous anaerobic digestion reactors were operated under different SRT conditions ranging from 6.6 to 36.0 days for both samples. The biochemical parameters (kch, kli, kpr, um,ac, um,bu, um,pro, um,va, Kac, Kbu, Kpro, and Kva) for each SRT condition were determined. To modify ADM1, the equations obtained through linear regression were substituted into biochemical parameters as a function of SRT. The modified ADM1 demonstrated superior accuracy compared with conventional ADM1. This study implies the feasibility of optimizing biochemical parameters for modeling in response to changes in environmental variables.

Increasing aggregate size reduces single-cell organic carbon incorporation by hydrogel-embedded wetland microbes

Microbial degradation of organic carbon in sediments is impacted by the availability of oxygen and substrates for growth. To better understand how particle size and redox zonation impact microbial organic carbon incorporation, techniques that maintain spatial information are necessary to quantify elemental cycling at the microscale. In this study, we produced hydrogel microspheres of various diameters (100, 250, and 500 μm) and inoculated them with an aerobic heterotrophic bacterium isolated from a freshwater wetland (Flavobacterium sp.), and in a second experiment with a microbial community from an urban lacustrine wetland. The hydrogel-embedded microbial populations were incubated with 13C-labeled substrates to quantify organic carbon incorporation into biomass via nanoSIMS. Additionally, luminescent nanosensors enabled spatially explicit measurements of oxygen concentrations inside the microspheres. The experimental data were then incorporated into a reactive-transport model to project long-term steady-state conditions. Smaller (100 μm) particles exhibited the highest microbial cell-specific growth per volume, but also showed higher absolute activity near the surface compared to the larger particles (250 and 500 μm). The experimental results and computational models demonstrate that organic carbon availability was not high enough to allow steep oxygen gradients and as a result, all particle sizes remained well-oxygenated. Our study provides a foundational framework for future studies investigating spatially dependent microbial activity in aggregates using isotopically labeled substrates to quantify growth.

The ASM2d model with two-step nitrification can better simulate biological nutrient removal systems enriched with complete ammonia oxidizing bacteria (comammox Nitrospira)

The discovery of comammox Nitrospira, a complete ammonia-oxidizing microorganism belonging to the genus Nitrospira, has brought new insights into the nitrification process in wastewater treatment plants (WWTPs). The applicability of Activated Sludge Model No. 2 d with one-step nitrification (ASM2d-OSN) or two-step nitrification (ASM2d-TSN) for the simulation of the biological nutrient removal (BNR) processes of a full-scale WWTP in the presence of comammox Nitrospira was studied. Microbial analysis and kinetic parameter measurements showed comammox Nitrospira was enriched in the BNR system operated under low dissolved oxygen (DO) and long sludge retention time (SRT). The relative abundance of Nitrospira under the conditions of stage I (DO = 0.5 mg/L, SRT = 60 d) was about twice of that under stage II conditions (DO = 4.0 mg/L, SRT = 26 d), and the copy number of the comammox amoA gene for stage I was 33 times higher than that for stage II. Compared to the ASM2d-OSN model, the ASM2d-TSN model simulated the performance of the WWTP under stage I conditions better, and the Theil inequality coefficient values of all the tested water quality parameters were lower than using ASM2d-OSN. These results indicate that an ASM2d model with two-step nitrification is a better choice for the simulation of WWTPs with the presence of comammox.

A critical review on microbial ecology in the novel biological nitrogen removal process: Dynamic balance of complex functional microbes for nitrogen removal

The novel biological nitrogen removal process has been extensively studied for its high nitrogen removal efficiency, energy efficiency, and greenness. A successful novel biological nitrogen removal process has a stable microecological equilibrium and benign interactions between the various functional bacteria. However, changes in the external environment can easily disrupt the dynamic balance of the microecology and affect the activity of functional bacteria in the novel biological nitrogen removal process. Therefore, this review focuses on the microecology in existing the novel biological nitrogen removal process, including the growth characteristics of functional microorganisms and their interactions, together with the effects of different influencing factors on the evolution of microbial communities. This provides ideas for achieving a stable dynamic balance of the microecology in a novel biological nitrogen removal process. Furthermore, to investigate deeply the mechanisms of microbial interactions in novel biological nitrogen removal process, this review also focuses on the influence of quorum sensing (QS) systems on nitrogen removal microbes, regulated by which bacteria secrete acyl homoserine lactones (AHLs) as signaling molecules to regulate microbial ecology in the novel biological nitrogen removal process. However, the mechanisms of action of AHLs on the regulation of functional bacteria have not been fully determined and the composition of QS system circuits requires further investigation. Meanwhile, it is necessary to further apply molecular analysis techniques and the theory of systems ecology in the future to enhance the exploration of microbial species and ecological niches, providing a deeper scientific basis for the development of a novel biological nitrogen removal process.

Mathematical modelling of an intermittent anoxic/aerobic MBBR: Estimation of nitrification rates and energy savings

This study aimed at modelling the performance of a novel MBBR configuration, named A/O-MBBR, comprised of a pre-anoxic reactor, with an HRT of 4.5 h, coupled with an intermittent anoxic/aerobic MBBR (HRT = 6.8 h). The lab-scale system was fed with municipal wastewater with an average influent Total Ammonia Nitrogen (TAN) and total COD (TCOD) concentrations of 46 mg of TAN-N L^-1 and 310 mg TCOD L^-1. During the whole experimental period, TAN removal efficiency was always higher than 96%; denitrification was also very effective, achieving nitrate and nitrite concentrations in the effluent both lower than 5 mg NO_x-N L^-1 on average. Moreover, TCOD average removal efficiency was equal to 85%. Modelling was performed to investigate the nitrification efficacy enhancement; to this aim, a biofilm model was developed, adopting the equations for mixed-culture biofilms and the Activated Model Sludge n°1 (ASM1) for the biological processes rates. The model allowed to determine the maximum uptake rate for autotrophic growth (μ_A was 2.5 d^-1) and the semisaturation constant (K_OA was 0.2 mg O₂ L^-1), suggesting that the nitrification process was 3-fold faster than average and very effective at low oxygen concentrations. The model estimated that about 85% of TAN was removed by the biofilm and only the remaining part by suspended biomass in the bulk liquid. Finally, it was assessed that the A/O-MBBR configuration allowed for a 45-60% savings of the energy requirement compared to a Benchmark WWTP layout.

The r/K selection theory and its application in biological wastewater treatment processes

Understanding the characteristics of functional organisms is the key to managing and updating biological processes for wastewater treatment. This review, for the first time, systematically characterized two typical types of strategists in wastewater treatment ecosystems via the r/K selection theory and provided novel strategies for selectively enriching microbial community. Functional organisms involved in nitrification (e.g., Nitrosomonas and Nitrosococcus), anammox (Candidatus Brocadia), and methanogenesis (Methanosarcinaceae) are identified as r-strategists with fast growth capacities and low substrate affinities. These r-strategists can achieve high pollutant removal loading rates. On the other hand, other organisms such as Nitrosospira spp., Candidatus Kuenenia, and Methanosaetaceae, are characterized as K-strategists with slow growth rates but high substrate affinities, which can decrease the pollutant concentration to low levels. More importantly, K-strategists may play crucial roles in the biodegradation of recalcitrant organic pollutants. The food-to-microorganism ratio, mass transfer, cell size, and biomass morphology are the key factors determining the selection of r-/K-strategists. These factors can be related with operating parameters (e.g., solids and hydraulic retention time), biomass morphology (biofilm or granules), and operating modes (continuous-flow or sequencing batch), etc., to achieve the efficient acclimation of targeted r-/K-strategists. For practical applications, the concept of substrate flux was put forward to further benefit the selective enrichment of r-/K-strategists, fulfilling effective management and improvement of engineered pollution control bioprocesses. Finally, the future perspectives regarding the development of the r/K selection theory in wastewater treatment processes were discussed.

Inorganic carbon limitation during nitrogen conversions in sponge-bed trickling filters for mainstream treatment of anaerobic effluent

Anaerobic sewage treatment is a proven technology in warm climate regions, and sponge-bed trickling filters (SBTFs) are an important post-treatment technology to remove residual organic carbon and nitrogen. Even though SBTFs can achieve a reasonably good effluent quality, further process optimization is hampered by a lack of mechanistic understanding of the factors influencing nitrogen removal, notably when it comes to mainstream anaerobically treated sewage. In this study, the factors that control the performance of SBTFs following anaerobic (i.e., UASB) reactors for sewage treatment were investigated. A demo-scale SBTF fed with anaerobically pre-treated sewage was monitored for 300 days, showing a median nitrification efficiency of 79% and a median total nitrogen removal efficiency of 26%. Heterotrophic denitrification was limited by the low organic carbon content of the anaerobic effluent. It was demonstrated that nitrification was impaired by a lack of inorganic carbon rather than by alkalinity limitation. To properly describe inorganic carbon limitation in models, bicarbonate was added as a state variable and sigmoidal kinetics were applied. The resulting model was able to capture the overall long-term experimental behaviour. There was no nitrite accumulation, which indicated that nitrite oxidizing bacteria were little or less affected by the inorganic carbon limitation. Overall, this study indicated the vital role of influent characteristics and operating conditions concerning nitrogen conversions in SBTFs treating anaerobic effluent, thus facilitating further process optimization.

Bifurcation Analysis of an Impulsive System Describing Partial Nitritation and Anammox in a Hybrid Reactor

Low-energy nitrogen removal under mainstream conditions is a technology that has received significant attention in recent years as the water industry drives toward long-term sustainability goals. Simultaneous partial nitritation-Anammox (PN/A) is one process that can provide substantial energy reduction and lower sludge yields. Mathematical modeling of the PN/A process offers engineers insights into the operating conditions necessary to maximize its potential. Laureni et al. (Laureni et al. Water Res. 2019, 14) have recently published a simplified mechanistic model of the process operated as a sequencing batch reactor that investigated the effect of three key operating parameters on performance (Anammox biofilm activity, dissolved oxygen concentration and fraction of solids wasted). The analysis of the model was limited, however, to simulation with relatively few discrete parameter sets. Here, we demonstrate through the use of bifurcation theory applied to an impulsive dynamical system that the parameter space can be partitioned into regions in which the system converges to different fixed points that represent different outcomes: either the washout of nitrite-oxidizing bacteria or their survival. Mapping process performance data onto these spaces allows engineers to target suitable operating regimes for specific objectives. Here, for example, we note that the nitrogen removal efficiency is maximized close to the curve that separates the regions in parameter space where nitrite-oxidizing bacteria washout from the region in which they survive. Further, control of solids washout and Anammox biofilm activity can also reduce oxygen requirements while maintaining an appropriate hydraulic retention time. The approach taken is significant given the possibility for using such a methodology for models of increasing complexity. This will enable engineers to probe the entire parameter space of systems of higher dimension and realism in a consistent manner.

A novel technology with precise oxygen-input control: Application of the partial nitritation-anammox process

Reliable and accurate oxygen-input control, which is critical to maintaining efficient nitrogen removal performance for partial nitritation-anammox (PN-A) process, remains one of the main operational difficulties. In this study, a novel, yet simple system (a simple process for autotrophic nitrogen-removal, SPAN) with precise oxygen-input control was developed to treat ammonium-rich wastewater via PN-A process. SPAN brings oxygen to biomass by circulating water and creating water spray (shower) at the water-air interface, and effectively balances the activities of core functional microorganisms through precise oxygen-input control. The oxygen-input rate is decided by the water circulation rate and shower rate and is measurable and predictable. Therefore, the required amount of oxygen for ammonium oxidation can be precisely delivered to the biomass by adjusting the circulation rate and shower rate. The results of two parallel SPAN reactors demonstrated that during long-term operation, the required oxygen input was precisely and reliably controlled. More than 99% of NH₄⁺-N and 81% - 85% of total nitrogen were stably removed, with anammox bacteria contributing to more than 96% of total nitrogen removal. Anammox bacteria were efficiently enriched to the highest level among the key nitrogen-converting microbial groups, both in terms of abundance (8.17%) and nitrogen-conversion capacity, while ammonium oxidizing bacteria were well controlled to provide sufficient ammonium-oxidizing capacity. Nitrite oxidizing bacteria were maintained stable (relative abundance of 1.08%-1.88%) and their activity was effectively suppressed. This study provided a novel technology, SPAN, to precisely control oxygen input in PN-A system, and proved that SPAN was effective and reliable in achieving long-term high-efficiency nitrogen removal.

Nitrifying biofilms deprived of organic carbon show higher functional resilience to increases in carbon supply

In nitrifying biofilms, the organic carbon to ammonia nitrogen (C/N) supply ratio can influence resource competition between heterotrophic and nitrifying bacteria for oxygen and space. We investigated the impact of acute and chronic changes in carbon supply on inter-guild competition in two moving bed biofilm reactors (MBBR), operated with (R1) and without (R0) external organic carbon supply. The microbial and nitrifying community composition of the reactors differed significantly. Interestingly, acute increases in the dissolved organic carbon inhibited nitrification in R1 ten times more than in R0. A sustained increase in the carbon supply decreased nitrification efficiency and increased denitrification activity to a greater extent in R1, and also increased the proportion of potential denitrifiers in both bioreactors. The findings suggest that autotrophic biofilms subjected to increases in carbon supply show higher nitrification and lower denitrification activity than carbon-fed biofilms. This has significant implications for the design of nitrifying bioreactors. Specifically, efficient removal of organic matter before the nitrification unit can improve the robustness of the bioreactor to varying influent quality. Thus, maintaining a low C/N ratio is important in nitrifying biofilters when acute carbon stress is expected or when anoxic activity (e.g. denitrification or H₂S production) is undesirable, such as in recirculating aquaculture systems (RAS).

Evolutionary causes and consequences of metabolic division of labour: why anaerobes do and aerobes don't

Metabolic division of the labour of organic matter decomposition into several steps carried out by different types of microbes is typical for many anoxic - but not oxic environments. An explanation of this well-known pattern is proposed based on the combination of three key insights: (i) well-studied anoxic environments are high flux environments: they are only anoxic because their high organic matter influx leads to oxygen depletion; (ii) shorter, incomplete catabolic pathways provide the capacity for higher flux, but this capacity is only advantageous in high flux environments; (iii) longer, complete catabolic pathways have energetic happy ends but only with high redox potential electron acceptors. Thus, aerobic environments favour longer pathways. Bioreactors, in contrast, are high flux environments and therefore favour division of catabolic labour even if aeration keeps them aerobic; therefore, host strains and feeding strategies must be carefully engineered to resist this pull.

Modelling anaerobic, aerobic and partial nitritation-anammox granular sludge reactors - A review

Wastewater treatment processes with granular sludge are compact and are becoming increasingly popular. Interest has been accompanied by the development of mathematical models. This contribution simultaneously reviews available models in the scientific literature for anaerobic, aerobic and partial nitritation-anammox granular sludge reactors because they comprise common phenomena (e.g. liquid, gas and granule transport) and thus pose similar challenges. Many of the publications were found to have no clearly defined goal. The importance of a goal is stressed because it determines the appropriate model complexity and helps other potential users to find a suitable model in the vast amount of literature. Secondly, a wide variety was found in the model features. This review explains the chosen modelling assumptions based on the different reactor types and goals wherever possible, but some assumptions appeared to be habitual within fields of research, without clear reason. We therefore suggest further research to more clearly define the range of operational conditions and goals for which certain simplifying assumptions can be made, e.g. when intragranule solute transport can be lumped in apparent kinetics and when biofilm models are needed, which explicitly calculate substrate concentration gradients inside granules. Furthermore, research is needed to better mechanistically understand detachment, removal of influent particulate matter and changes in the mixing behaviour inside anaerobic systems, before these phenomena can be adequately incorporated in models. Finally, it is suggested to perform full-scale model validation studies for aerobic and anammox reactors. A spreadsheet in the supplementary information provides an overview of the features in the 167 reviewed models.

The future of WRRF modelling - outlook and challenges

The wastewater industry is currently facing dramatic changes, shifting away from energy-intensive wastewater treatment towards low-energy, sustainable technologies capable of achieving energy positive operation and resource recovery. The latter will shift the focus of the wastewater industry to how one could manage and extract resources from the wastewater, as opposed to the conventional paradigm of treatment. Debatable questions arise: can the more complex models be calibrated, or will additional unknowns be introduced? After almost 30 years using well-known International Water Association (IWA) models, should the community move to other components, processes, or model structures like 'black box' models, computational fluid dynamics techniques, etc.? Can new data sources - e.g. on-line sensor data, chemical and molecular analyses, new analytical techniques, off-gas analysis - keep up with the increasing process complexity? Are different methods for data management, data reconciliation, and fault detection mature enough for coping with such a large amount of information? Are the available calibration techniques able to cope with such complex models? This paper describes the thoughts and opinions collected during the closing session of the 6th IWA/WEF Water Resource Recovery Modelling Seminar 2018. It presents a concerted and collective effort by individuals from many different sectors of the wastewater industry to offer past and present insights, as well as an outlook into the future of wastewater modelling.

Modelling aerobic granular sludge reactors through apparent half-saturation coefficients

During biological wastewater treatment, substrates undergo simultaneous diffusion and reactions inside microbial aggregates, creating microscale spatial substrate gradients and limiting the macroscale reaction rates. For flocculent and anaerobic granular sludge, this rate-limiting effect of diffusion is often lumped in model parameters, like the half-saturation coefficients of Monod kinetics in activated sludge models (ASM). Yet, an explicit description of the reaction-diffusion process with biofilm models is more common for aerobic granular sludge. This work investigates whether apparent half-saturation coefficients could have applications for aerobic granular sludge as well and examines the implications of this simplification. To this end, the macroscopic reaction rates predicted with a one-dimensional biofilm (1D) model were fitted with Monod kinetics. The results showed that the macroscale rates could indeed be described using apparent kinetics, at the very least over a time scale where the microbial population distribution stays fixed. However, the coefficients were sensitive to changes in the microbial population distribution, which can be affected by long-term changes in operating conditions. Also the activity of organisms that compete for the same substrates affect the parameter value. Be that as it may, apparent kinetics also depend on the operating conditions for flocculent and anaerobic granular sludge, but they have still been used successfully for design and optimization. Therefore, the last section of this work illustrates that they may also have applications for aerobic granular sludge. A simple model for ammonium removal using apparent half-saturation coefficients for oxygen and ammonium is applied to a full-scale reactor, taking advantage of the batch-wise operation and on-line monitoring data for regular recalibration.

High-Rate Partial Nitritation of Municipal Wastewater after Psychrophilic Anaerobic Pretreatment

Partial nitritation/anammox can provide energy-efficient nitrogen removal from the main stream of municipal wastewater. The main bottleneck is the growth of nitrite oxidizing bacteria (NOB) at low temperatures (<15 °C). To produce effluent suitable for anammox, real municipal wastewater after anaerobic pretreatment was treated by enriched ammonium oxidizing bacteria (AOB) in suspended sludge SBR at 12 °C. NOB were continually washed out using aerobic duration control strategy (ADCS). Solids retention time was set to 9-16 days. Using this approach, average ammonia conversion higher than 57% at high oxidation rate of 0.4 ± 0.1 kg-N kg-VSS^-1 d^-1 was achieved for more than 100 days. Nitrite accumulation (N-NO₂^-/N-NO_X) of 92% was maintained. Thus, consistently small amounts of present NOB were efficiently suppressed. Our mathematical model explained how ADCS enhanced the inhibition of NOB growth via NH₃ and HNO₂. This approach will produce effluent suitable for anammox even under winter conditions in mild climates.

Intermittent Aeration Suppresses Nitrite-Oxidizing Bacteria in Membrane-Aerated Biofilms: A Model-Based Explanation

Autotrophic ammonium oxidation in membrane-aerated biofilm reactors (MABRs) can make treatment of ammonium-rich wastewaters more energy-efficient, especially within the context of short-cut ammonium removal. The challenge is to exclusively enrich ammonium-oxidizing bacteria (AOB). To achieve nitritation, strategies to suppress nitrite-oxidizing bacteria (NOB) are needed, which are ideally grounded on an understanding of underlying mechanisms. In this study, a nitrifying MABR was operated under intermittent aeration. During eight months of operation, AOB dominated, while NOB were suppressed. On the basis of dissolved oxygen (DO), ammonium, nitrite, and nitrate profiles within the biofilm and in the bulk, a 1-dimensional nitrifying biofilm model was developed and calibrated. The model was utilized to explore the potential mechanisms of NOB suppression associated with intermittent aeration, considering DO limitation, direct pH effects on enzymatic activities, and indirect pH effects on activity via substrate speciation. The model predicted strong periodic shifts in the spatial gradients of DO, pH, free ammonia, and free nitrous acid, associated with aerated and nonaerated phases. NOB suppression during intermittent aeration was mostly explained by periodic inhibition caused by free ammonia due to periodic transient pH upshifts. Dissolved oxygen limitation did not govern NOB suppression. Different intermittent aeration strategies were then evaluated for nitritation success in intermittently aerated MABRs: both aeration intermittency and duration were effective control parameters.

Biofilm Thickness Influences Biodiversity in Nitrifying MBBRs-Implications on Micropollutant Removal

In biofilm systems for wastewater treatment (e.g., moving bed biofilms reactors-MBBRs) biofilm thickness is typically not under direct control. Nevertheless, biofilm thickness is likely to have a profound effect on the microbial diversity and activity, as a result of diffusion limitation and thus substrate penetration in the biofilm. In this study, we investigated the impact of biofilm thickness on nitrification and on the removal of more than 20 organic micropollutants in laboratory-scale nitrifying MBBRs. We used novel carriers (Z-carriers, AnoxKaldnes) that allowed controlling biofilm thickness at 50, 200, 300, 400, and 500 μm. The impact of biofilm thickness on microbial community was assessed via 16S rRNA gene amplicon sequencing and ammonia monooxygenase (amoA) abundance quantification through quantitative PCR (qPCR). Results from batch experiments and microbial analysis showed that (i) the thickest biofilm (500 μm) presented the highest specific biotransformation rate constants (kbio, L g(-1) d(-1)) for 14 out of 22 micropollutants; (ii) biofilm thickness positively associated with biodiversity, which was suggested as the main factor for the observed enhancement of kbio; (iii) the thinnest biofilm (50 μm) exhibited the highest nitrification rate (gN d(-1) g(-1)), amoA gene abundance and kbio values for some of the most recalcitrant micropollutants (i.e., diclofenac and targeted sulfonamides). Although thin biofilms favored nitrification activity and the removal of some micropollutants, treatment systems based on thicker biofilms should be considered to enhance the elimination of a broad spectrum of micropollutants.