DETERMINING STOICHIOMETRIC PARAMETERS USING AN LSS RESPIROMETER FOR MATHEMATICAL MODELING OF CONSTRUCTED WETLANDS

Similar documents
Modelling of Wastewater Treatment Plants

USING NUMERICAL SIMULATION SOFTWARE FOR IMPROVING WASTEWATER TREATMENT EFFICIENCY

SIMULATION AND CALIBRATION OF A FULL- SCALE SEQUENCING BATCH REACTOR FOR WASTEWATER TREATMENT

NEW BIOLOGICAL PHOSPHORUS REMOVAL CONCEPT SUCCESSFULLY APPLIED IN A T-DITCH PROCESS WASTEWATER TREATMENT PLANT

Figure Trickling Filter

Feedforward aeration control of a Biocos wastewater treatment plant

ENHANCING THE PERFORMANCE OF OXIDATION DITCHES. Larry W. Moore, Ph.D., P.E., DEE Professor of Environmental Engineering The University of Memphis

Examples of Studies conducted by

Removal of High C and N Contents in Synthetic Wastewater Using Internal Circulation of Anaerobic and Anoxic/Oxic Activated Sludge Processes

Effect of the start-up length on the biological nutrient removal process

336098: DYNAMIC MODELLING AND SIMULATION OF ANAEROBIC DIGESTER FOR HIGH ORGANIC STRENGTH WASTE

19. AEROBIC SECONDARY TREATMENT OF WASTEWATER

Biofilm Reactor Technology and Design

Simulation of a wastewater treatment plant receiving industrial effluents #

Calibration of a Bio-Kinetic Model to Simulate Microalgae Growth

Sanitary and Environmental Engineering I (4 th Year Civil)

Activated sludge modeling of Vale Faro Waste Water Treatment Plant

Water and Wastewater Engineering Dr. Ligy Philip Department of Civil Engineering Indian Institute of Technology, Madras

Sustainable Energy Generation in Microbial Fuel Cell Catalyzed with Bacillus Subtilis Species

SIMULTANEOUS ESTIMATION OF NITRIFICATION/DENITRIFICATION KINETICS AND INFLUENT NITROGEN LOAD USING ORP AND DO DYNAMICS

Journal of Biotechnology

THE INFLUENCE OF CARRIER SIZE AND SHAPE IN THE MOVING BED BIOFILM PROCESS

Kinetic models for nitrification inhibition by ammonium and nitrite in a suspended and an immobilised biomass systems

LABORATORY PROJECT: IMPACT OF VARIABLE LOADING ON A FIXED FILM REACTOR VERSUS A SUSPENDED GROWTH REACTOR

Experimental determination of the heterotroph anoxic yield in anoxic-aerobic activated sludge systems treating municipal wastewater

IWA Publishing 2012 Water Practice & Technology Vol 7 No 3 doi: /wpt

DETERMINATION OF BIOLOGICAL KINETIC CONSTANTS USING RESPIROMETRY FOR THE WATER9 AIR EMISSIONS MODEL

EXPERIMENTAL DETERMINATION OF THE HETEROTROPH ANOXIC YIELD IN ANOXIC-AEROBIC ACTIVATED SLUDGE SYSTEMS TREATING MUNICIPAL WASTEWATER

19. AEROBIC SECONDARY TREATMENT OF WASTEWATER

A General Description of the Activated Sludge Model No. 1 (ASM1)

Supplemental Information for

Modeling Reactions Between Activated Sludge Fractions and Ozone to Optimize Biosolids Reduction Processes

Mechanistic modelling of pulp and paper mill wastewater treatment plants

Module 19 : Aerobic Secondary Treatment Of Wastewater. Lecture 24 : Aerobic Secondary Treatment Of Wastewater

Study of Kinetic coefficients of a Membrane Bioreactor (MBR) for municipal wastewater treatment

Application of MBR for the treatment of textile wastewater

Biological Short-Cut Nitrogen Removal from Anaerobic Digestate in a Demonstration Sequencing Batch Reactor

Nitrogen removal during the treatment of wastewater and landfill leachate in a biological aerated filter

Troubleshooting Activated Sludge Processes. PNCWA - Southeast Idaho Operators Section Pocatello, ID February 11, 2016 Jim Goodley, P.E.

WASTEWATER TREATMENT PLANT OPTIMIZATION USING A DYNAMIC MODEL APPROACH

American Water College 2010

Augmentation of self-purification capacity of sewer pipe by immobilizing microbes on the pipe surface

Activated Sludge Process Optimisation AERATION PERFORMANCE. Application Note XAUK WW

Determination of design criteria for UASB reactors as a wastewater pretreatment system in tropical small communities

Optimal sizing of Activated Sludge Process with ASM3

IMPROVED BIO-TOWER TECHNOLOGY MUNICIPAL SEWAGE TREATMENT

MARPAK modular biomedia WASTEWATER TREATMENT

AEROBIC GRANULATION - AN ECONOMICALLY VIABLE OPTION FOR THE TREATMENT OF WASTEWATER

TWO YEARS OF BIOLOGICAL PHOSPHORUS REMOVAL WITH AN ADVANCED MSBR SYSTEM AT THE SHENZHEN YANTIAN WASTEWATER TREATMENT PLANT

Dynamic Performance of UASB Reactor Treating Municipal Wastewater

JET-LOOP SYSTEM. The Jet-Loop System applies for aerobic wastewater treatment using ejector(s) as the mechanical system for aeration.

Design of horizontal and vertical subsurface flow constructed wetlands treating industrial wastewater

NUTRIENT OPTIMIZATION FOR PULP & PAPER WASTEWATER TREATMENT PLANTS AN OPPORTUNITY FOR MAJOR COST SAVINGS

ADVANCED CONTROL SYSTEM OF ACTIVATED SLUDGE PROCESSES USING IN- SITU AMMONIUM AND NITRATE PROBES

Domestic Waste Water (Sewage): Collection, Treatment & Disposal

Energy and chemical efficient nitrogen removal at a full-scale MBR water reuse facility

AquaNereda Aerobic Granular Sludge Technology

Comparing biofilm models for a single species biofilm system

Benchmarking Biological Nutrient Removal in Wastewater Treatment Plants: Influence of Mathematical Model Assumptions

Environmental Biotechnology Cooperative Research Centre Date submitted: March 2008 Date published: March 2011

Aeration University Advanced Concepts in Energy Efficiency

8.8 Calculation of oxygen requirements Peder Maribo

Assessment of Clogging in Constructed Wetlands

Unit Treatment Processes in Water and Wastewater Engineering

Study on the Heavy Metals Removal Efficiencies of Constructed Wetlands with Different Substrates

Two-dimensional cellular automaton model for mixed-culture biofilm

MODELLING OF THE MARIANRIDGE WASTEWATER TREATMENT PLANT

Inlet Process air and wash cycle scour air. Air grid Floor Nozzle

KINETIC ANALYSIS AND SIMULATION OF UASB ANAEROBIC TREATMENT OF A SYNTHETIC FRUIT WASTEWATER

Lagoons Operation and Management in New Brunswick

Anaerobic Reactor Technologies

IMPROVING THE DESIGN OF A FULL SCALE WASTEWATER TREATMENT PLANT WITH THE USE OF THE COMPLEX ACTIVATED SLUDGE MODEL

Activated Sludge Processes

Chapter 11. Secondary Clarifiers

Comparison of Water Quality Parameters

KEY WORDS Anaerobic lagoons ; reaction constant ; operation ; per capita sludge ; loading rate.

A. Damayanti*, Mohd Razman Salim, Zaini Ujang

Effect of the temperature on the performance of a sludge activated petrochemical wastewater treatment plant

Oxygen. Oxygen is one of the fundamental resources required by life forms on Earth. Aquatic ecosystems have a wide assortment of life forms.

Chapter 2: Conventional Wastewater Treatment (continue)

SEQUENCING BATCH REACTOR FOR ORGANICS AND NITROGEN REMOVAL: MODELING AND DESIGN OPTIMIZATION VELMURUGAN SUBRAMANIAN

Activated Sludge Treatment Extended Aeration (EA) and Sequencing Batch Reactor (SBR) Systems

Degradation of wine distillery wastewaters by the combination of aerobic biological treatment with chemical oxidation by Fenton s reagent

MATHEMATICAL MODELING OF HOUSEHOLD WASTEWATER TREATMENT BY DUCKWEED BATCH REACTOR

Fixed-Film Processes

Lecture 3 CE 433. Excerpts from Lecture notes of Professor M. Ashraf Ali, BUET.

THE USE OF ULTRASOUND TO ACCELERATE THE ANAEROBIC DIGESTION OF WASTE ACTIVATED SLUDGE. Ewa Neczaj, Joanna Lach

Effect of Temperature Shock on Activities of Phosphorus-accumulating Organisms

CTB3365x Introduction to Water Treatment

WEFTEC 2002 RECIPROCATING SUBSURFACE TREATMENT SYSTEM KEEPS AIRPORT OUT OF THE DEEP FREEZE

Bypass filtration in cooling water systems

BIODEGRADATION OF DE-ICING COM- POUNDS IN COLUMNS SIMULATING A RANGE OF NATURAL CONDITIONS

Sludge Formation In The Activated Sludge Process: Mathematical Analysis

Hach Method Luminescence Measurement of Dissolved Oxygen (LDO ) in Water and Wastewater

BIOLOGICAL PHOSPHOROUS REMOVAL AN OPERATOR S GUIDE

Modelling of an Upflow Anaerobic Sludge Blanket reactor

Bioreactor System ERT 314. Sidang /2011

ENVE 302 Environmental Engineering Unit Processes DENITRIFICATION

Assessment of Biodegradability of Synthetic Tanning Agents Used in Leather Tanning Process

COMBINED CONTROL OF SEWER AND WASTEWATER TREATMENT PLANT UNSING IDENTIFIABLE MODELS THE SMArt CONTROL OF WASTEWATER SYSTEMS PROJECT

Transcription:

DETERMINING STOICHIOMETRIC PARAMETERS USING AN LSS RESPIROMETER FOR MATHEMATICAL MODELING OF CONSTRUCTED WETLANDS Inês Isabel da Silva Baranda July 2017 Civil Engineering Department, Instituto Superior Técnico, Universidade Técnica de Lisboa, Portugal ABSTRACT Despite being a natural treatment, constructed wetlands (from now on referred as CW) constitute an extremely complex system due to numerous physicochemical transformations that occur simultaneously. As a result, these solutions have been studied in greater detail using this model. Most mathematical models currently available are based on proposed models for activated sludge (ASM) systems, which often comprise many parameters, both kinetic and stoichiometric. An example of this is CWM1- Constructed Wetland Model 1, one of the most popular models in CW. Due to a very small number of studies dedicated to determine model parameters for constructed wetlands, modeling studies for these systems often rely on the parameters obtained for activated sludge. This study proposed to measure maximum heterotrophic growth yield (Y H ) of the biofilm from two constructed wetlands with subsurface flow using a LSS respirometer (Liquid-phase principle: Static gas, liquid Static). Tests were performed to biomass from two CW supplied with different organic load, however the study of CW feed with a lower load is in an early stage. Biofilm samples, used in respirometric tests, were obtained using 100 to 300 ml bulk volume of gravel and washing vigorously with 1.5 L of tap water. Tests had duration between 8 to 15 hours and sodium acetate was used as a soluble substrate. The oxygen uptake curve, from the CW fed with higher organic loads, only revealed consumption of added substrate. The range of the stoichiometric Y H was from 0.64 to 0.74 gcqogcqo -1 comprising the Y H range of values for activated sludge and constructed wetlands with vertical flow. Respirometric tests, carried out from CW biofilm fed with a lower load, revealed two distinct stages as observed in previous studies. First stage was associated with readily biodegradable substrate consumption and the second, associated with storage products consumption. Keywords constructed wetlands, respirometry, Y H, oxygen uptake rate curve.

1 INTRODUCTION Constructed wetlands are a sustainable solution for the biological treatment of effluents. This type of treatment takes advantage of natural degradation mechanisms in wetlands, for the pollutants removal on wastewater. The advantages of this type of solution include its optimum landscape integration and high sustainability, regarding reduced energy consumption and chemical reagents. Despite being a natural treatment, built wetlands constitute an extremely complex system due to many physical and chemical transformations that occur simultaneously. Besides that, the increased knowledge about biological degradation mechanisms of wastewater treatment processes, resulted in the elaboration of several mathematical models, currently used as instruments for design, operation and optimization of treatment infrastructures. For that reason, constructed wetlands have been studied in greater detail with the use of modeling. Most of the currently available mathematical models are based on activated sludge (ASM) or wetland constructed models with parameters imported from activated sludge. It s the case of CWM1- Constructed Wetland model 1 (Langergraber et al. 2009). This model comprises a high number of parameters, both kinetic and stoichiometric, deduced for activated sludge. The obvious lack of specific information for this type of biological treatment indicates how important are developing studies, in way to obtain kinetic and stoichiometric parameters related to the constructed wetlands. Currently, respirometric method represents an important support tool for modeling biological treatment systems, since it allows to obtain some necessary parameters for models calibration that could not be directly measured or determined by analytical methods. Respirometric consists in measuring the Dissolved Oxygen (DO) dynamics during substrate consumption by the biomass, resulting in an Oxygen Uptake Rate (OUR) curve. These curves allow the estimation of values for parameters related to the aerobic bioconversion of substrates and the identification of different consumption phases. The present work aims to understand the development of biofilm from constructed wetlands using the respirometric method to determine the maximum heterotrophic growth yield (Y H ). 2 MATERIALS AND METHODS 2.1 Respirometer set-up In this study, the respirometer used is an LSS type which works with a liquid-phase measurement principle and a static gas and static liquid operational regimen. It was developed at Instituto Superior Técnico. Used respirometer includes a main reactor opened at the top but closed at the bottom, built from an acrylic cylinder. The main reactor has inlet-outlet ports for connection of the tubes. Respirometer also includes a reading cell, without aeration, made of glass. The system also includes an air compressor

(HAILEA V-20, China), a hydraulic pump (Flojet, China), a mechanical stirrer (BS VELP, Italy) and a dissolved oxygen sensor (ProODO, YSI, USA), which also gives temperature readings. The Figure 1, shows schematically all the elements that compose the respirometer system. Figure 1 LSS respirometer used during the experimental phase Air compressor is connected through a rubber tube to the main reactor bottom, in order to provide necessary aeration to the suspended biomass. In this respirometric system, liquid circulates between reactor and reading cell, where the probe is inserted. Connection between these two elements is made by a rubber and a cap, which prevents oxygen to be transferred between the cell and the outside. DO sensor is connected to a computer, where readings are registered. The selected DO sensor uses a luminescence-based measurement principle, allowing less fouling problems. Used pump, that makes liquid recirculate between reactor and reading cell, is also connected to a computer, allowing the control of pumping cycle times. For this experiment, OUR tests were conducted with a cycle, including a 3-minute idle period with no recirculation (static liquid) and a 2-minute pumping period. Since, during the idle period, there is no renewal of dissolved oxygen, the decrease of dissolved oxygen results only from consumption by the biomass. Mechanical stirrer is used to promote a continuous biomass suspension in the main reactor, while air compressor aeration keeps dissolved oxygen levels close to saturation, so that DO is not a limiting factor in the OUR measurement. 2.2 Biofilm source and sampling This study was carried out with biofilm samples, taken from two horizontal subsurface flow, constructed wetland installed in Hydraulics Laboratory of Instituto Superior Técnico, in Lisbon,

Portugal. Both constructed wetlands were planted with Scirpus holoschoenus, on a gravel media with a porosity of 30%. The main difference between constructed wetland 2 (CW 2) and constructed wetland 9 (CW 9) is the way that they are fed. The CW 2 is fed with synthetic sewage with a Chemical Oxygen Demand (COD) content of 3838 mgo 2 L -1. Feeding was done at a rate of 10 L per day, five days a week, with a rest period during weekends. The CW 9 was fed 5 times a week with quantities much lower than that CW 2, resulting from diluted leftovers from synthetic sewage. Although, it was not quantified in precise terms the COD, is estimated that the mass load is approximately 10 go 2 / m² day. To obtain each biofilm sample, a volume between 100 and 300 ml of gravel was carefully removed from constructed wetland center. In order to detach the biofilm, 1,5 L of tap water at room temperature was added to the gravel sample followed by a vigorous manual shaking period, after which most of the biofilm was observed to be in suspension. Finally this volume is placed in the main reactor by adding tap water to a final volume of approximately 2.5L. 2.3 Respirometric test Respirometric test allows to determine the OUR curve corresponding to the oxidation of a specific amount of rapidly biodegradable substrate. The OUR curve obtained corresponds to the line that best fits the values of DO concentration, and which minimizes distances of measured values along time, in idle periods. After placing the sample in main reactor, the biofilm suspension was aerated until endogenous respiration level was attained. This procedure ensures that oxygen consumption, during the test, was due only to the oxidation of the added substrate. Once endogenous respiration level was attained, were added 200mL of a sodium acetate solution corresponding to a carbonaceous substrate, which will allow measuring the oxygen consumption capacity of biofilm. Acetate concentration, in the added solution, was adjusted in order to give COD levels in reactor, between 83 and 519 mgl -1. This values range is based on previous studies, where several concentrations of acetate were tested inside the reactor. COD content acetate solution was measured using cuvette tests (25-1500 mg/l,merck). Respirometric ends when endogenous respiration is again attained. Depending on the amount of substrate added, the tests had durations between 8 and 15 hours. In order to avoid interference from nitrification, a nitrification inhibitor (NTH 166 600, WTW) was added to biomass suspension at the beginning of each respirometric test. 2.4 Calculations for Y S e Y STO In real waste water treatment, biofilm is exposed to variable concentrations of different substrates. It means that microorganisms experienced two contrasting situations: feast conditions, when external

substrate is available in excess, and famine conditions, in the absence of external substrate. In order to assure biomass response on this variability, Gujer et al. (1999) proposed the ASM3 model that introduces a phenomenon of carbon storage sources, observed by several researchers. This model assumes that substrate is first stored before being used for biofilm growth. However, experimental data have demonstrated the existence of simultaneous storage and growth processes during feast conditions. In the famine phase, microorganisms consume previously created storage products, which function as a sort of reserve. Based on this, the SSAG (Simultaneous Storage And Growth) model, proposed as a complement to the ASM3 model, (Hoque et al. 2010). In this study, identifying the period of readily biodegradable organic matter (S s ) consumption and the period of storage products consumption (that ends when OUR settles back at endogenous level), was performed visually, following methodology proposed by Andreottola et al. (2007) and Ortigara et al. (2011), as represented in Figure 2, which is a graphical representation of DO consumption profiles associated with these two phases: the first one corresponds to a period in which the rapidly biodegradable substrate is consumed, and the second to a period in which the storage products are consumed. Figure 2 Graphical visualization of used method to distinguish between initial (YS) and storage (YSTO) yield values Besides that, Figure 2 also indicates the calculation procedure for two yield values, according to a simplified method proposed by Ortigara et al. (2011). The initial growth yield value (Y S ) was obtained using equation 1 below. t f t0 Y S = 1 ΔO 2 acetate = 1 [OUR(t) OUR endogenous(t) OUR STO(t)]dt S S S S (eq.1) Where S S is the initial concentration of added substrate (acetate - COD, mgo 2 L -1 ), ΔO 2 acetate will be the amount of oxygen consumed for the oxidation of added substrate, OUR (t) is the oxygen consumption rate at the end of the test, OUR endogenous (t) is the endogenous respiration level, OUR STO is the oxygen consumption rate at the beginning of famine phase, t 0 is the instant substrate is added and t f is the instant that all the added substrate has been consumed.

The storage yield was then calculated according to equation 2. Y STO = 1 ΔO 2 STO S S (eq.2) Where S S is the initial concentration of added substrate (acetate - CDO, mgo 2 L -1 ) and ΔO 2 STO is the oxygen consumption associated to storage The growth yield coefficient, Y H refers to the total growth of the heterotrophic biomass by aerobic conditions. If the OUR curve do not consumption of stored products, Y H is equal to Y S. 2.5 Temperature corretion Temperature is a parameter with a strong influence on microorganism s metabolism, affecting organic matter s oxidation rate. In general, higher temperatures accelerate the biological decomposition reactions of organic matters, thus influencing its removal efficiency by the organisms. Because respirometry interprets and evaluates the metabolism of microorganisms, OUR curves, obtained at room temperature, are affected by daily temperature variations. As such, it isn t easy to accurately identify endogenous respiration. In order to eliminate the influence of temperature variations, OUR curve was corrected by adopting a reference temperature of 20 C. So, in respirometry tests performed, temperature in the reading cell was continuously registered and its correction was performed using equation 3 indicated in Ortigara et al. (2011). OUR 20 C = α OUR T (T 20 C) (α = 1.08) (3) Where OUR 20 C is the oxygen consumption rate corrected to the temperature of 20 C (mgo 2 /Lh), OUR T is the oxygen consumption rate at a temperature T (mgo 2 /Lh), T is the temperature measured at each value of OUR ( C) and α is the Arrhenius constant. 3 Results and Discussion Respirometry tests results, performed in CW2, revealed an OUR profile with only one phase related to the external substrate consumption. After acetate addition, a high OUR peak can be identified. After this initial high value, OUR decreases down to the endogenous level. During this decreasing period, OUR curve shows only one slope. Note that the endogenous OUR level increases from its initial value. That result probably occurs from biomass growth.

Figure 3 OUR profile obtained for CW2 Biofilm from CW2 was already studied in the past, so, this study made a comparison with the previous one. In 2014, CW were fed with a COD content of 757 mgo 2 L -1. As Figure 3 reveals, OUR profile no longer presents two phases as it was in the past. Comparing with carried out tests in the same CW, the most significant change is the absence of evidence of storage. This difference in behavior is mainly due to the increase in the amount of organic load of the influent solution with which the CW was fed. In 2014, when the CW was fed with a lower amount of organic matter, the available substrate amount limited the growth of biofilm, in some periods. As such, microorganisms capable of producing and storing polymers had a competitive advantage, compared to microorganisms without such ability, and tests performed revealed the existence of storage phenomenon. This mechanism allowed them to survive after food exhaustion. However, when CW was fed with an amount of organic matter about five times higher (January 2016), storage was no longer a fundamental mechanism for microorganism s survival. Although biofilm was also subject to dynamic conditions of substrate availability, the duration of absence periods was lower, considering higher feed concentration in the medium culture. The concentration change of the affluent solution probably favored the microorganisms with direct substrate use to the detriment of storage solutions, influencing the biodegradation process of organic matter.

Figure 4 Comparison with respirometric test performed in the past Respirometry results tests, performed in CW9, revealed an OUR profile with two different phases (Figure 6). The first one is related to the external substrate consumption (feast conditions), while the second one corresponds to the consumption of the previously stored products (famine conditions). After acetate addition, a high OUR peak can be identified. After reaching the maximum level, OUR drops to a level higher than the endogenous respiration, followed by a gradual decrease down to endogenous level. During this decreasing period, two slopes can be identified: an initial higher slope, related to external substrate consumption, and a second lower slope, related to the previously stored products consumption. Figure 5 OUR profile obtained for CW9

4 Conclusions The respirometric tests performed with biofilm from constructed wetlands feed with higher organic load (3838 mgo 2 L -1 ), didn t showed evidence consumption of the previously stored products but only direct consumption of the external substrate. However, constructed wetlands fed with lower organic load, presented two distinct phases. The first is related to the external substrate consumption, while the second one corresponds to the consumption of the products previously stored. The ranges of values obtained for Y S and Y STO comprise the range of values found in literature for proposed models of activated sludge systems as well as the value obtained by Ortigara et al. (2011) for a vertical flow constructed wetland. It was also found that biofilm sample, required for a respirometric test, had a lower volume in the case of constructed wetlands, fed with higher organic load. This fact is probably due to a lower development of the biofilm in constructed wetlands fed with less organic load, resulting from the less availability of substrate. In this way, the importance of properly determining the substrate / biomass ratio in each CW was identified because this will allow a better understanding of the development of the biofilm through the stoichiometric coefficients. This study will continue to determine the transition conditions between the two types of observed behavior, namely the direct consumption of substrate and the consumption of substrate followed by consumption of storage products. The results obtained in this study are relevant for application in several studies, namely modeling, since they allow a better understanding of the systems under study, they also make possible the optimization of their treatment efficiency. References Andreottola, G., Oliveira, E., Foladori, P., Peterlini, R., & Ziglio, G. (2007). Método Respirométrico Para O Monitoramento de Processos Biológicos. Water Science and Technology. Guisasola, a, Sin, G., Baeza, J. a, Carrera, J., and Vanrolleghem, P. a. (2005). Limitations of ASM1 and ASM3: a comparison based on batch oxygen uptake rate profiles from different CD fullscale wastewater treatment plants. Water Science and Technology, 52(10-11), 69 77. Gujer, W., Henze, M., Mino, T., & van Loosdrecht, M. (1999). Activated sludge model no. 3. Water Science and Technology, 39 (1), 183 193. Henze, M., Gujer, W., Mino, T., and Loosdrecht, V. (2000). Activated Sludge Models ASM1, ASM2, ASM2d AND ASM3. IWA Publishing.

Langergraber, G., Rousseau, D., Garcıa, J., & Mena, J. (2009). CWM1: a general model to describe biokinetic processes in subsurface flow constructed wetlands. Water Science & Technology, 59(9):1687-97. Ortigara, A., Foladori, P., & Andreottola, G. (2011). Kinetics of heterotrophic biomass and storage mechanism in wetland cores measured by respirometry. Water Science and Technology. Pisoeiro, J., Galvão, A., Pinheiro, H., Ferreira, F., & Matos, J. (2017). Determining stoichiometric parameters of detached biomass from a HSSF-CW using respirometry. Elsevier Ecological Engineering Volume 98, 388 393.