In the pursuit of reducing carbon footprint and in view of increasing energy prices, energy efficiency is more important than ever before. Batch-wise operated aerobic granular sludge reactors consume up to 50% less energy compared to conventional activated sludge systems because pumping energy is reduced and mixing equipment is not needed. Further energy reduction efforts should therefore target aeration energy requirements. The alpha factor is an important factor influencing the oxygen transfer efficiency, however the dynamic behaviour of alpha has hardly been investigated in general and never for an aerobic granular sludge reactor. This study showed that alpha increases during the aeration phase of a cycle due to the influence of different process parameters. Through a data analysis study of 175 batch cycles of the Prototype Nereda® installation in Utrecht over the summer and winter period of 2020–2021, the exchange ratio and temperature were identified as the main influencing factors on the rate of increase of alpha in a batch cycle. A higher exchange ratio was related to a slower increase in alpha over the aeration phase, while a higher temperature was related to a faster increase in alpha. Moreover, alpha was characterized by a same minimal value at the beginning of every aeration phase, which could be explained by the adsorption of soluble biodegradable organic carbon described by a Langmuir adsorption model. Two mathematical models, a decreasing exponential and a first order model, were set up to unravel the dynamic behaviour of alpha. Both models were discussed in view of their practical implications for the design and performance optimization of aerobic granular sludge reactors and other batch-wise operated aerobic wastewater treatment systems.@en

The aerobic granular sludge (AGS) process treats wastewater with a significantly lower footprint and energy consumption compared to conventional activated sludge systems. Nevertheless, there is still potential for optimizing its performance, and mathematical models are most valuable tools to this end. Aeration energy consumption deserves particular attention, as it is the largest remaining operating cost for AGS systems. Batch-wisely operated reactors show an increasing oxygen transfer efficiency during aeration, which translates into a dynamic alpha factor. However, the dynamic nature of alpha is neglected in current models. The impact of this simplification on the operating performance was addressed for the first time in this study. Through the development of a novel 1-D biofilm reactor model, calibrated to a full-scale AGS plant, it was shown that the alpha dynamics affect both model structure and calibration, as well as the process performance. The description of the dynamic nature of alpha through the empirical relationship with the soluble biodegradable organic carbon required the addition of the state variable representing soluble slowly biodegradable organic carbon (SCB) to the biokinetic ASM2d model. Simulation results showed that alpha dynamics significantly influences simultaneous nitrification and denitrification and therefore need to be included in mathematical models to optimize AGS process performance. Different process variables such as volume exchange ratio, aeration capacity and granule size can be manipulated to improve reactor design and performance. The practical application of these new insights were discussed regarding the optimization of AGS systems, as well as other batch-wisely operated aerobic wastewater treatment systems.@en

This chapter provides a review of the models available for estimating the production and emission of methane from wastewater collection and treatment systems. The details of a number of mechanistic models as well as the simplified empirical models have been summarized. Their limitations have been identified and general methods for calibration and validation have been presented.@en

Partial nitration-anammox is a resource-efficient technology for nitrogen removal from wastewater. However, the advantages of this nitrogen removal technology are challenged by the emission of N2O, a potent greenhouse gas. In this study, a granular sludge one-stage partial nitritation-anammox reactor comprising granules and flocs was run for 337 days in the presence of influent organics to investigate its effect on N removal and N2O emissions. Besides, the effect of aeration control strategies and flocs removal was investigated as well. The interpretation of the experimental results was complemented with modelling and simulation. The presence of influent organics (1 g COD g-1 N) helped to suppress NOB and significantly reduced the overall N2O emissions while having no significant effect on anammox activity. Besides, long-term monitoring of the reactor indicated that constant airflow rate control resulted in more stable effluent quality and less N2O emissions than DO control. Still, floc removal reduced N2O emissions at DO control but increased N2O emissions at constant airflow rate. Furthermore, anammox bacteria could significantly reduce N2O production during heterotrophic denitrification, likely via competition for NO with heterotrophs. Overall, this study demonstrated that the presence of influent organics together with proper aeration control strategies and floc management could significantly reduce the N2O emissions without compromising nitrogen removal efficiency during one-stage partial nitritation-anammox processes.@en

Gas bubbles are introduced in water to absorb or strip volatile substances in a variety of unit operations, for example during (waste)water treatment. To calculate the transfer rate of substances between the liquid phase and the gas phase, different assumptions have been made in literature regarding the gas phase composition and hydraulic pressure, which both vary along the reactor height. In this study, analytical expressions were derived for the total (macroscopic) liquid-gas transfer rate, using either the complete gradients of the mole fraction and pressure (comprehensive approach) or a uniform value, for one or both of them. Simulations with the comprehensive model were performed to understand the effect of the type of volatile substance and of the reactor design and operating conditions on the total liquid-gas transfer rate. These effects were found to be highly interactive and often non-linear. Next, the simulation results of the comprehensive model were compared with those from models that assume either a uniform mole fraction or a uniform pressure in the complete reactor volume. This illustrated that for soluble substances, the mole fraction gradient strongly affects the total liquid-gas transfer rate, while the pressure gradient became only important under operating conditions that promote stripping (i.e., for a high concentration in the liquid phase and low concentration in the inlet gas). For very poorly soluble substances, the pressure became more important under conditions that promote absorption. These results on the importance of the mole fraction and pressure gradients remained equally valid when explicitly considering a typical variation of the volumetric overall transfer coefficient (KLa) along the reactor height. Finally, a simple and fast procedure was made available through a spreadsheet to select appropriate simplifying assumptions in reactor or plant-wide models. By applying the procedure to oxygen (O2), carbon dioxide (CO2), methane (CH4), nitrous oxide (N2O) and nitrogen gas (N2) in an aerobic biological wastewater treatment reactor, it was demonstrated that some common simplifications can lead to significant errors, for which corrections were proposed.@en

Gas bubbles are introduced in water to absorb or strip volatile substances in a variety of unit operations, for example during (waste)water treatment. To calculate the transfer rate of substances between the liquid phase and the gas phase, different assumptions have been made in literature regarding the gas phase composition and hydraulic pressure, which both vary along the reactor height. In this study, analytical expressions were derived for the total (macroscopic) liquid-gas transfer rate, using either the complete gradients of the mole fraction and pressure (comprehensive approach) or a uniform value, for one or both of them. Simulations with the comprehensive model were performed to understand the effect of the type of volatile substance and of the reactor design and operating conditions on the total liquid-gas transfer rate. These effects were found to be highly interactive and often non-linear. Next, the simulation results of the comprehensive model were compared with those from models that assume either a uniform mole fraction or a uniform pressure in the complete reactor volume. This illustrated that for soluble substances, the mole fraction gradient strongly affects the total liquid-gas transfer rate, while the pressure gradient became only important under operating conditions that promote stripping (i.e., for a high concentration in the liquid phase and low concentration in the inlet gas). For very poorly soluble substances, the pressure became more important under conditions that promote absorption. These results on the importance of the mole fraction and pressure gradients remained equally valid when explicitly considering a typical variation of the volumetric overall transfer coefficient (KLa) along the reactor height. Finally, a simple and fast procedure was made available through a spreadsheet to select appropriate simplifying assumptions in reactor or plant-wide models. By applying the procedure to oxygen (O2), carbon dioxide (CO2), methane (CH4), nitrous oxide (N2O) and nitrogen gas (N2) in an aerobic biological wastewater treatment reactor, it was demonstrated that some common simplifications can lead to significant errors, for which corrections were proposed.@en

In municipal wastewater treatment plants, some dissolved methane can enter the aerobic bioreactors. This greenhouse gas originates from sewers and return flows from anaerobic sludge treatment. In well-mixed conventional activated sludge reactors, methane emissions are largely avoided because methane oxidizing bacteria consume a large fraction, even without optimizing for this purpose. In this work, the fate of dissolved methane is studied in aerobic granular sludge reactors, as they become increasingly popular. The influence of the characteristic design and operating conditions of these reactors are studied with a mathematical model with apparent conversion kinetics and stripping: the separation of feeding and aeration in time, a higher substrate transport resistance, a high retention time of granular biomass and a taller water column. Even for a best-case scenario combining an unrealistically low intragranule substrate transport resistance, a high retention time, a tall reactor, an extremely high influent methane concentration and no oxygen limitation, the methane conversion efficiency was only 12% when feeding and aeration were separated in time, which is lower than for continuous activated sludge reactors under typical conditions. A more rigorous model was used to confirm the limited conversion, considering the multi-species and multi-substrate biofilm kinetics, anoxic methane consumers and the high substrate concentration at the bottom during upward plug flow feeding. The observed limited methane conversion is mainly due to the high concentration that accumulates during unaerated feeding phases, which favours stripping more than conversion in the subsequent aeration phase. Based on these findings, strategies were proposed to mitigate methane emissions from wastewater treatment plants with sequentially operated reactors.@en

In municipal wastewater treatment plants, some dissolved methane can enter the aerobic bioreactors. This greenhouse gas originates from sewers and return flows from anaerobic sludge treatment. In well-mixed conventional activated sludge reactors, methane emissions are largely avoided because methane oxidizing bacteria consume a large fraction, even without optimizing for this purpose. In this work, the fate of dissolved methane is studied in aerobic granular sludge reactors, as they become increasingly popular. The influence of the characteristic design and operating conditions of these reactors are studied with a mathematical model with apparent conversion kinetics and stripping: the separation of feeding and aeration in time, a higher substrate transport resistance, a high retention time of granular biomass and a taller water column. Even for a best-case scenario combining an unrealistically low intragranule substrate transport resistance, a high retention time, a tall reactor, an extremely high influent methane concentration and no oxygen limitation, the methane conversion efficiency was only 12% when feeding and aeration were separated in time, which is lower than for continuous activated sludge reactors under typical conditions. A more rigorous model was used to confirm the limited conversion, considering the multi-species and multi-substrate biofilm kinetics, anoxic methane consumers and the high substrate concentration at the bottom during upward plug flow feeding. The observed limited methane conversion is mainly due to the high concentration that accumulates during unaerated feeding phases, which favours stripping more than conversion in the subsequent aeration phase. Based on these findings, strategies were proposed to mitigate methane emissions from wastewater treatment plants with sequentially operated reactors.@en

In municipal wastewater treatment plants, some dissolved methane can enter the aerobic bioreactors. This greenhouse gas originates from sewers and return flows from anaerobic sludge treatment. In well-mixed conventional activated sludge reactors, methane emissions are largely avoided because methane oxidizing bacteria consume a large fraction, even without optimizing for this purpose. In this work, the fate of dissolved methane is studied in aerobic granular sludge reactors, as they become increasingly popular. The influence of the characteristic design and operating conditions of these reactors are studied with a mathematical model with apparent conversion kinetics and stripping: the separation of feeding and aeration in time, a higher substrate transport resistance, a high retention time of granular biomass and a taller water column. Even for a best-case scenario combining an unrealistically low intragranule substrate transport resistance, a high retention time, a tall reactor, an extremely high influent methane concentration and no oxygen limitation, the methane conversion efficiency was only 12% when feeding and aeration were separated in time, which is lower than for continuous activated sludge reactors under typical conditions. A more rigorous model was used to confirm the limited conversion, considering the multi-species and multi-substrate biofilm kinetics, anoxic methane consumers and the high substrate concentration at the bottom during upward plug flow feeding. The observed limited methane conversion is mainly due to the high concentration that accumulates during unaerated feeding phases, which favours stripping more than conversion in the subsequent aeration phase. Based on these findings, strategies were proposed to mitigate methane emissions from wastewater treatment plants with sequentially operated reactors.@en

This work shows how more variables can be monitored with a single off-gas sampler on sequentially operated than on continuously fed and aerated reactors and applies the methods to data from a full-scale aerobic granular sludge reactor as a demonstration and to obtain insight in this technology. First, liquid-gas transfer rates were calculated. Oxygen (O2) absorption and carbon dioxide (CO2) emission rates showed comparable cyclic trends due to the coupling of O2 consumption and CO2 production. Methane (CH4) emissions showed a stripping profile and nitrous oxide (N2O) emissions showed two peaks each cycle, which were attributed to different production pathways. Secondly, aeration characteristics were calculated, of which the gradual improvement within cycles was explained by surfactants degradation. Thirdly, liquid phase concentrations were estimated from off-gas measurements via a novel calculation procedure. As such, an average influent CH4 concentration of 0.7 g·m−3 was found. Fourthly, reaction rates could be estimated from off-gas data because no feeding or discharge occurred during reaction phases. The O2 consumption rate increased with increasing dissolved oxygen and decreased once nitrification was complete. Fifthly, greenhouse gas emissions could be derived, indicating a 0.06% N2O emission factor. Sixthly, off-gas gave an indication of influent characteristics. The CO2 emitted per kg COD catabolized corresponded with the TOC/COD ratio of typical wastewater organics in cycles with balanced nitrification and denitrification. High nitrogen removal efficiencies were associated with high catabolized COD/N ratios as estimated from the O2 absorption. Finally, mass balances could be closed using off-gas O2 data. As such, an observed yield of 0.27 g COD/g COD was found. All these variables could be estimated with a single sampler because aeration without feeding creates a more homogeneous off-gas composition and simplifies liquid-phase mass balances. Therefore, off-gas analyzers may have a broader application potential for sequentially operated reactors than currently acknowledged.@en

This work shows how more variables can be monitored with a single off-gas sampler on sequentially operated than on continuously fed and aerated reactors and applies the methods to data from a full-scale aerobic granular sludge reactor as a demonstration and to obtain insight in this technology. First, liquid-gas transfer rates were calculated. Oxygen (O2) absorption and carbon dioxide (CO2) emission rates showed comparable cyclic trends due to the coupling of O2 consumption and CO2 production. Methane (CH4) emissions showed a stripping profile and nitrous oxide (N2O) emissions showed two peaks each cycle, which were attributed to different production pathways. Secondly, aeration characteristics were calculated, of which the gradual improvement within cycles was explained by surfactants degradation. Thirdly, liquid phase concentrations were estimated from off-gas measurements via a novel calculation procedure. As such, an average influent CH4 concentration of 0.7 g·m−3 was found. Fourthly, reaction rates could be estimated from off-gas data because no feeding or discharge occurred during reaction phases. The O2 consumption rate increased with increasing dissolved oxygen and decreased once nitrification was complete. Fifthly, greenhouse gas emissions could be derived, indicating a 0.06% N2O emission factor. Sixthly, off-gas gave an indication of influent characteristics. The CO2 emitted per kg COD catabolized corresponded with the TOC/COD ratio of typical wastewater organics in cycles with balanced nitrification and denitrification. High nitrogen removal efficiencies were associated with high catabolized COD/N ratios as estimated from the O2 absorption. Finally, mass balances could be closed using off-gas O2 data. As such, an observed yield of 0.27 g COD/g COD was found. All these variables could be estimated with a single sampler because aeration without feeding creates a more homogeneous off-gas composition and simplifies liquid-phase mass balances. Therefore, off-gas analyzers may have a broader application potential for sequentially operated reactors than currently acknowledged.@en

This study deals with the effect of aeration control strategies on the nitrogen removal efficiency and nitrous oxide (N2O) emissions in a partial nitritation–anammox reactor with granular sludge. More specifically, dissolved oxygen (DO) control, constant airflow and effluent ammonium (NH4+) control strategies were compared through a simulation study. Particular attention was paid to the effect of flocs, which are deliberately or unavoidable present besides granules in this type of reactor. When applying DO control, DO setpoints had to be adjusted to the amount of flocs present in the reactor to maintain high nitrogen removal and reduce N2O emissions, which is difficult to realize in practice because of variable floc fractions. Constant airflow rate control could maintain a good nitrogen removal efficiency independent of the floc fraction in the reactor, but failed in N2O mitigation. Controlling aeration based on the effluent ammonium concentration results in both high nitrogen removal and relatively low N2O emissions, also in the presence of flocs. Fluctuations in floc fractions caused significant upsets in nitrogen removal and N2O emissions under DO control but had less effect at constant airflow and effluent ammonium control. Still, rapid and sharp drops in flocs led to a peak in N2O emissions at constant airflow and effluent ammonium control. Overall, effluent ammonium control reached the highest average nitrogen removal and lowest N2O emissions and consumed the lowest aeration energy under fluctuating floc concentrations.@en

This study deals with the effect of aeration control strategies on the nitrogen removal efficiency and nitrous oxide (N2O) emissions in a partial nitritation–anammox reactor with granular sludge. More specifically, dissolved oxygen (DO) control, constant airflow and effluent ammonium (NH4+) control strategies were compared through a simulation study. Particular attention was paid to the effect of flocs, which are deliberately or unavoidable present besides granules in this type of reactor. When applying DO control, DO setpoints had to be adjusted to the amount of flocs present in the reactor to maintain high nitrogen removal and reduce N2O emissions, which is difficult to realize in practice because of variable floc fractions. Constant airflow rate control could maintain a good nitrogen removal efficiency independent of the floc fraction in the reactor, but failed in N2O mitigation. Controlling aeration based on the effluent ammonium concentration results in both high nitrogen removal and relatively low N2O emissions, also in the presence of flocs. Fluctuations in floc fractions caused significant upsets in nitrogen removal and N2O emissions under DO control but had less effect at constant airflow and effluent ammonium control. Still, rapid and sharp drops in flocs led to a peak in N2O emissions at constant airflow and effluent ammonium control. Overall, effluent ammonium control reached the highest average nitrogen removal and lowest N2O emissions and consumed the lowest aeration energy under fluctuating floc concentrations.@en

Multivariate statistical analysis was applied to investigate the dependencies and underlying patterns between N2O emissions and online operational variables (dissolved oxygen and nitrogen component concentrations, temperature and influent flow-rate) during biological nitrogen removal from wastewater. The system under study was a full-scale reactor, for which hourly sensor data were available. The 15-month long monitoring campaign was divided into 10 sub-periods based on the profile of N2O emissions, using Binary Segmentation. The dependencies between operating variables and N2O emissions fluctuated according to Spearman's rank correlation. The correlation between N2O emissions and nitrite concentrations ranged between 0.51 and 0.78. Correlation >0.7 between N2O emissions and nitrate concentrations was observed at sub-periods with average temperature lower than 12 °C. Hierarchical k-means clustering and principal component analysis linked N2O emission peaks with precipitation events and ammonium concentrations higher than 2 mg/L, especially in sub-periods characterized by low N2O fluxes. Additionally, the highest ranges of measured N2O fluxes belonged to clusters corresponding with NO3-N concentration less than 1 mg/L in the upstream plug-flow reactor (middle of oxic zone), indicating slow nitrification rates. The results showed that the range of N2O emissions partially depends on the prior behavior of the system. The principal component analysis validated the findings from the clustering analysis and showed that ammonium, nitrate, nitrite and temperature explained a considerable percentage of the variance in the system for the majority of the sub-periods. The applied statistical methods, linked the different ranges of emissions with the system variables, provided insights on the effect of operating conditions on N2O emissions in each sub-period and can be integrated into N2O emissions data processing at wastewater treatment plants.@en

Multivariate statistical analysis was applied to investigate the dependencies and underlying patterns between N2O emissions and online operational variables (dissolved oxygen and nitrogen component concentrations, temperature and influent flow-rate) during biological nitrogen removal from wastewater. The system under study was a full-scale reactor, for which hourly sensor data were available. The 15-month long monitoring campaign was divided into 10 sub-periods based on the profile of N2O emissions, using Binary Segmentation. The dependencies between operating variables and N2O emissions fluctuated according to Spearman's rank correlation. The correlation between N2O emissions and nitrite concentrations ranged between 0.51 and 0.78. Correlation >0.7 between N2O emissions and nitrate concentrations was observed at sub-periods with average temperature lower than 12 °C. Hierarchical k-means clustering and principal component analysis linked N2O emission peaks with precipitation events and ammonium concentrations higher than 2 mg/L, especially in sub-periods characterized by low N2O fluxes. Additionally, the highest ranges of measured N2O fluxes belonged to clusters corresponding with NO3-N concentration less than 1 mg/L in the upstream plug-flow reactor (middle of oxic zone), indicating slow nitrification rates. The results showed that the range of N2O emissions partially depends on the prior behavior of the system. The principal component analysis validated the findings from the clustering analysis and showed that ammonium, nitrate, nitrite and temperature explained a considerable percentage of the variance in the system for the majority of the sub-periods. The applied statistical methods, linked the different ranges of emissions with the system variables, provided insights on the effect of operating conditions on N2O emissions in each sub-period and can be integrated into N2O emissions data processing at wastewater treatment plants.@en

This work demonstrates the usefulness of non-linear data reconciliation to evaluate available measurements and estimate unmeasured variables for a full-scale partial nitritation (SHARON) reactor for the treatment of wastewater with high ammonium concentrations. Despite a lack of measured data, the bilinear approach of formulating system constraints allowed to satisfy the requirements for data reconciliation and gross error detection, leading to a balanced data set and the estimation of unmeasured variables.@en

This work demonstrates the usefulness of non-linear data reconciliation to evaluate available measurements and estimate unmeasured variables for a full-scale partial nitritation (SHARON) reactor for the treatment of wastewater with high ammonium concentrations. Despite a lack of measured data, the bilinear approach of formulating system constraints allowed to satisfy the requirements for data reconciliation and gross error detection, leading to a balanced data set and the estimation of unmeasured variables.@en

This work demonstrates the usefulness of non-linear data reconciliation to evaluate available measurements and estimate unmeasured variables for a full-scale partial nitritation (SHARON) reactor for the treatment of wastewater with high ammonium concentrations. Despite a lack of measured data, the bilinear approach of formulating system constraints allowed to satisfy the requirements for data reconciliation and gross error detection, leading to a balanced data set and the estimation of unmeasured variables.@en

This study deals with the potential and the limitations of dynamic models for describing and predicting nitrous oxide (N2O) emissions associated with biological nitrogen removal from wastewater. The results of a three-week monitoring campaign on a full-scale partial nitritation reactor were reproduced through a state-of-the-art model including different biological N2O formation pathways. The partial nitritation reactor under study was a SHARON reactor treating the effluent from a municipal wastewater treatment plant sludge digester. A qualitative and quantitative comparison between experimental data and simulation results was performed to identify N2O formation pathways as well as for model identification. Heterotrophic denitrifying bacteria and ammonium oxidizing bacteria (AOB) were responsible for N2O formation under anoxic conditions, whereas under aerated conditions the AOB were the most important N2O producers. Relative to previously proposed models, hydroxylamine (NH2OH) had to be included as a state variable in the AOB conversions in order to describe potential N2O formation by AOB under anoxic conditions. An oxygen inhibition term in the corresponding reaction kinetics was required to fairly represent the relative contribution of the different AOB pathways for N2O production. Nevertheless, quantitative prediction of N2O emissions with models remains a challenge, which is discussed.@en

This study deals with the potential and the limitations of dynamic models for describing and predicting nitrous oxide (N2O) emissions associated with biological nitrogen removal from wastewater. The results of a three-week monitoring campaign on a full-scale partial nitritation reactor were reproduced through a state-of-the-art model including different biological N2O formation pathways. The partial nitritation reactor under study was a SHARON reactor treating the effluent from a municipal wastewater treatment plant sludge digester. A qualitative and quantitative comparison between experimental data and simulation results was performed to identify N2O formation pathways as well as for model identification. Heterotrophic denitrifying bacteria and ammonium oxidizing bacteria (AOB) were responsible for N2O formation under anoxic conditions, whereas under aerated conditions the AOB were the most important N2O producers. Relative to previously proposed models, hydroxylamine (NH2OH) had to be included as a state variable in the AOB conversions in order to describe potential N2O formation by AOB under anoxic conditions. An oxygen inhibition term in the corresponding reaction kinetics was required to fairly represent the relative contribution of the different AOB pathways for N2O production. Nevertheless, quantitative prediction of N2O emissions with models remains a challenge, which is discussed.@en