Master Water Modeling &
Constructed Wetland Design

A complete interactive course to understand the science, design principles, and real-world application of AquaSai's Multi-Stage Recirculating (MSR) constructed wetland systems.

8
Modules
40+
Quiz Questions
3
Live Calculators
3
Case Studies

What You'll Master

This course takes you from first principles to professional practice. You will understand why constructed wetlands work at a molecular and ecosystem level, how to design them using verified scientific models, and where AquaSai has applied them in the real world.

Module 1
CW Fundamentals
Types, MSR stages, ecological basis
Module 2
Hydraulics & Flow
Darcy's Law, HLR, HRT, porosity
Module 3
Treatment Science
P-k-C*, Arrhenius, Monod kinetics
Module 4
Plant Ecology
Species, rhizosphere, polyculture
Module 5
Layout & MSR Design
HSSF, VSSF, FWS, sizing methods
Module 6
Case Studies
ETBAR, Bagmati, Koh Phangan
Module 7
Compliance & Ops
EU regulations, monitoring, O&M
Module 8
Final Assessment
Master-level quiz with certificate

The AquaSai Ecosystem

AquaSai is a nature-based water treatment venture anchored in the MSR (Multi-Stage Recirculating) Constructed Wetland — a modular, solar-powered, IoT-monitored system that achieves advanced wastewater treatment using engineered plant communities, stratified substrate, and controlled hydraulics.

Unlike conventional treatment plants that consume energy and chemicals, AquaSai systems generate ecosystem services: they sequester carbon, support biodiversity, produce biomass for biogas, and integrate into the local Water–Energy–Food nexus.

🔬

Science-First

Every design parameter is derived from peer-reviewed models: Kadlec & Wallace 2009, HYDRUS-CW2D, CWM1, Arrhenius kinetics.

☀️

Off-Grid Ready

Solar-powered pumping, IoT monitoring, minimal OPEX. Designed for remote and peri-urban communities.

🌏

Globally Deployed

Koh Phangan (Thailand), Kathmandu Bagmati (Nepal), ETBAR Tenerife (Spain). Each adapts MSR to local ecology.

📊

Verified Performance

Pre-validated with HYDRUS-CW2D simulations. Target BOD ≤ 10 mg/L, TSS ≤ 10 mg/L, E.coli ≤ 1000 CFU/100mL.

Did You Know

An AquaSai MSR constructed wetland can achieve BOD ≤ 10 mg/L and E.coli ≤ 1,000 CFU/100 mL using only plants, sunlight, gravity and gravel — with zero chemical dosing and a fraction of the energy of a conventional activated-sludge plant.

AquaSai Field Note

This course distils methods proven across three live AquaSai deployments — the Koh Phangan resort wetland (Thailand), the Bagmati river-edge system (Kathmandu, Nepal) and the ETBAR Ensenada Pelada coastal pilot (Tenerife, Spain).

— AquaSai engineering team
Try It LiveModel a full MSR train in your browser with AquaSai's design tools.
🛠️ MSR Tool🌊 HYDRUS Simulator

🌐 AquaSai MSR Ecosystem — Interactive Knowledge Map

Hover over any node to explore. Click to reveal details about each system component.

Core System
Technology
Science
Co-benefits
Projects
MSR Cutaway
Module 1
Constructed Wetland Fundamentals
Understand the ecological and engineering basis of nature-based water treatment, and how AquaSai's MSR system integrates all treatment mechanisms in synergy.
🎯 Learning Objectives
  • Distinguish between FWS, HSSF, VSSF and MSR constructed wetland typologies
  • Explain the six treatment mechanisms operating simultaneously in AquaSai MSR
  • Describe how recirculation links nitrification and denitrification across stages
  • Define primary, secondary and tertiary treatment roles in a multi-stage system

What is a Constructed Wetland?

A Constructed Wetland (CW) is an engineered shallow basin planted with aquatic macrophytes, filled with porous substrate (gravel, sand, coconut coir), through which wastewater flows and is treated by physical, chemical, and biological mechanisms working simultaneously.

Unlike a natural wetland that is preserved for biodiversity, a CW is deliberately designed and sized to achieve specific treatment targets — BOD removal, nitrification, denitrification, phosphorus sorption — validated against mathematical models.

Types of Constructed Wetlands

🌊 FWS — Free Water Surface

Water flows above substrate, exposed to atmosphere. Excellent for tertiary polishing, wildlife habitat. Lower capital, moderate BOD removal.

🔄 HSSF — Horizontal Subsurface Flow

Water flows horizontally through gravel below surface. Anaerobic/anoxic zones → excellent denitrification. Most common CW type worldwide.

⬇️ VSSF — Vertical Subsurface Flow

Water loaded from top, percolates down. Intermittent dosing creates aerobic conditions → excellent nitrification. Energy-efficient via gravity.

🔀 MSR — Multi-Stage Recirculating

AquaSai's integrated design: combines HSSF + VSSF + FWS + recirculation loop in sequence. Achieves what no single type can match alone.

The AquaSai MSR — 4-Stage Treatment Train

The MSR is the heart of every AquaSai project. Flow the wastewater through all four stages by clicking each box:

⚗️
Stage 1
Primary Settlement
Anaerobic pre-treatment, TSS removal, BOD reduction ≥45%
🌿
Stage 2
HSSF Wetland
Planted horizontal bed; anoxic/aerobic mosaic; BOD, NH₄, TSS removal
⬇️
Stage 3
VSSF + Recirc
Vertical flow → nitrification; recirculation loop enhances denitrification
🏊
Stage 4
Polishing Pond
FWS pond; algae, UV-C, final pathogen kill; E.coli ≤ 1000 CFU/100mL

Why MSR Out-Performs Single-Stage Systems

Each treatment mechanism has conditions under which it thrives — anaerobic bacteria need zero oxygen, nitrifiers need high oxygen, denitrifiers need anoxic conditions with available carbon. A single bed cannot simultaneously optimize all of them.

The MSR resolves this by staging: each cell operates in its optimal redox window, and the recirculation loop allows partially nitrified effluent to re-enter the HSSF, providing the nitrate that denitrifiers need and the carbon that nitrifiers lack.

Constructed Wetland Types — Performance Comparison
Normalised treatment efficiency across 6 key parameters (0–100, higher = better). Source: Kadlec & Wallace 2009; Vymazal 2014.
FWS Pond
HSSF Wetland
VSSF Wetland
AquaSai MSR

🗺️ Constructed Wetland Type Comparison — Knowledge Map

Hover each node to see hydraulic and treatment characteristics. All four types combine in the AquaSai MSR.

CW Family
HSSF
VSSF
FWS
MSR (Combined)

📝 Module 1 Quiz

1. Which constructed wetland type creates the most aerobic conditions, making it best for nitrification?
A
FWS — Free Water Surface
B
HSSF — Horizontal Subsurface Flow
C
VSSF — Vertical Subsurface Flow
D
Pond only (FWS Stage 4)
2. The MSR recirculation loop primarily benefits which removal process?
A
TSS sedimentation
B
Denitrification (NO₃ → N₂)
C
Phosphorus precipitation
D
UV disinfection
3. What is the BOD removal efficiency typically expected from Stage 1 (primary settlement/anaerobic pre-treatment)?
A
10–15%
B
40–55%
C
80–90%
D
95–99%
Did You Know

A single reed (Phragmites australis) can transport oxygen down to its roots and leak it into the substrate — creating aerobic micro-zones millimetres from anoxic ones. This radial oxygen loss is why one AquaSai MSR bed can host nitrifiers and denitrifiers side by side.

AquaSai Field Note

Every AquaSai MSR sequences HSSF → VSSF → FWS with a recirculation loop, so each cell operates in its optimal redox window instead of compromising in a single bed. The cutaway renders above are drawn from real AquaSai project schematics.

— AquaSai MSR design library
Try It LiveExplore the 4-stage MSR train interactively and size each cell.
🛠️ Open MSR Tool
Key Takeaways — Module 1
  • MSR is a multi-stage reed bed system combining HSSF, VSSF, FWS and recirculation
  • VSSF cells create the aerobic conditions needed for nitrification
  • The recirculation loop enables denitrification by returning nitrate-rich effluent to the anoxic HSSF stage
  • Stage 1 pre-treatment (primary settlement) removes 40–55% BOD before the wetland beds
Flow plan
Module 2
Hydraulics & Flow
Master the equations governing how water moves through porous substrate — the foundation of every wetland sizing calculation.
⏱ ~7 min read
🎯 Learning Objectives
  • Apply Darcy's Law (Q = K·i·A) to calculate flow velocity through porous substrate
  • Compute HLR and HRT for any given system geometry and flow
  • Predict how porosity, bed depth and surface area affect hydraulic retention time
  • Identify causes and remediation strategies for short-circuiting and clogging

Darcy's Law — The Foundation of Porous Flow

All subsurface flow in constructed wetlands is governed by Darcy's Law, which relates the volumetric flow rate through a porous medium to the hydraulic gradient and medium properties:

Darcy's Law — Porous Media Flow
$$Q = K \cdot A_c \cdot \frac{\Delta h}{L}$$
$Q$ = volumetric flow rate [m³/d]  |  $K$ = hydraulic conductivity [m/d]  |  $A_c$ = cross-sectional area [m²]  |  $\Delta h / L$ = hydraulic gradient [–]
AquaSai gravel: K ≈ 10,000–100,000 m/d (new bed); drops to 500–2,000 m/d after biofilm development.

For HSSF wetlands, the hydraulic conductivity K of the gravel substrate typically ranges from 10,000–100,000 m/d. However, biological clogging (biofilm growth) progressively reduces K over time — a key long-term design consideration that MSR addresses through VSSF pulse dosing that aerates and prevents clogging.

Hydraulic Loading Rate (HLR)

HLR — Surface Loading Rate
HLR = Q / As
HLR [m³/m²·d = m/d], Q = daily flow [m³/d], As = surface area [m²]

HLR is the most important hydraulic design parameter. It determines contact time, velocity through substrate, and is the core input for the P-k-C* treatment model. Typical AquaSai MSR HLR: 0.05–0.15 m/d (conservative, ensuring full treatment).

Hydraulic Retention Time (HRT)

HRT — Contact Time in the System
HRT = (As · d · n) / Q
HRT [days], d = depth [m], n = porosity [–], Q = flow [m³/d]

Porosity (n) of the substrate is the fraction of void space available for water storage. Typical values: gravel n=0.35–0.40, coconut coir n=0.50–0.65. A higher porosity means more water is in contact with biofilm-covered surfaces for longer — better treatment at the same area.

🔢 Interactive Hydraulics Calculator

80 m³/d
800 m²
0.6 m
0.38
HLR
0.10
m/d
HRT
2.28
days
Pore Volume
182.4
Loading Rate
OK
HRT vs. Surface Area — Trade-off Visualisation
For Q = 80 m³/d, d = 0.6 m, n = 0.38. Move the sliders above to update this chart in real time.

Key Design Ranges (AquaSai Standard)

ParameterHSSF CellVSSF CellPolishing Pond
HLR (m/d)0.05–0.150.08–0.250.10–0.30
HRT (days)2–51–33–7
Depth (m)0.4–0.80.5–0.80.5–1.5
Porosity0.35–0.420.38–0.45N/A (open water)
L:W ratio3:1 – 10:11:1 – 2:11:1 – 3:1

📝 Module 2 Quiz

1. A wetland has Q = 100 m³/d and a surface area of 1,000 m². What is the HLR?
A
0.01 m/d
B
0.10 m/d
C
10.0 m/d
D
100 m/d
2. If Q = 50 m³/d, As = 600 m², depth = 0.6 m, porosity = 0.40 — what is the HRT?
A
0.6 days
B
1.4 days
C
2.88 days
D
7.2 days
3. Which substrate property causes HLR performance to degrade over time in HSSF systems?
A
Increased porosity from plant root growth
B
Increased temperature in the substrate
C
Biological clogging reducing hydraulic conductivity K
D
Evapotranspiration losses increasing concentration
Did You Know

Darcy's Law (1856) — written for sand filters in Dijon — still governs how water creeps through the gravel of a modern constructed wetland. AquaSai sizes every HSSF bed so the hydraulic gradient never exceeds the substrate's conductivity, preventing surface breakout and short-circuiting.

AquaSai Field Note

On the Bagmati river-edge system, seasonal monsoon flows swing by an order of magnitude. AquaSai uses HRT and hydraulic-loading design margins (2–5 d in HSSF, 1–3 d in VSSF) so treatment holds in both wet and dry seasons.

— AquaSai Kathmandu project data
Try It LiveRun Darcy and HRT scenarios in the HYDRUS-CW2D simulator.
🌊 HYDRUS Simulator
Key Takeaways — Module 2
  • Darcy's Law (Q = K·i·A) governs subsurface flow; K is the single most critical substrate parameter
  • HLR (m/yr) = Q / As — the primary sizing metric linking flow to treatment area
  • HRT (d) = As · d · n / Q — residence time determines treatment contact time
  • Biological clogging reduces effective K over time; conservative design uses 20–50% of clean-gravel K
Rhizosphere
Module 3
Treatment Science
Learn the mathematical models that predict pollutant removal — from the P-k-C* model through Arrhenius temperature corrections to Monod microbial kinetics.
⏱ ~10 min read
🎯 Learning Objectives
  • Apply the P-k-C* model to predict effluent quality from any inlet concentration and area
  • Use the Arrhenius equation to temperature-correct rate constants for your site climate
  • Explain the role of the irreducible background concentration C* in treatment limits
  • Interpret Monod kinetics and half-saturation constants for nitrification biofilms

The P-k-C* Model (Kadlec & Wallace 2009)

The most widely validated model for constructed wetland performance is the P-k-C* (P-order rate constant background concentration) model. For first-order kinetics (P=1), the effluent concentration is:

P-k-C* Model — First Order (Kadlec & Wallace 2009)
$$C_e = (C_0 - C^*) \cdot e^{-k_T / \text{HLR}} + C^*$$
$C_e$ = effluent [mg/L]  |  $C_0$ = influent [mg/L]  |  $C^*$ = background concentration [mg/L]
$k_T$ = Arrhenius-corrected rate constant [m/yr]  |  HLR = hydraulic loading rate [m/yr]
Worked example: C0=250 mg/L BOD, k20=37.5 m/yr, T=20°C (θ=1.000, kT=37.5), HLR=37.5 m/yr, C*=3.5. Ce = (250−3.5)·exp(−37.5/37.5)+3.5 = 246.5·e⁻¹+3.5 ≈ 246.5×0.368+3.5 ≈ 94.2 mg/L. Removal = (250−94.2)/250 = 62.3%.

The C* background concentration represents the irreducible minimum — wetlands cannot achieve zero output because internal biological processes continuously regenerate small amounts of each pollutant. Key C* values: BOD 3.5 mg/L, NH₄-N 0 mg/L, TN 1.5 mg/L, TP 0.02 mg/L.

Arrhenius Temperature Correction

Arrhenius Temperature Factor
$$k_T = k_{20} \cdot \theta^{(T-20)}$$
$k_{20}$ = rate constant at 20 °C  |  $\theta$ = temperature coefficient (1.000–1.090)  |  $T$ = water temperature [°C]
Nitrification (NH₄-N) most temperature-sensitive: θ=1.060 — performance drops ~40% at 10°C vs 20°C.

Rate Constants for Key Pollutants

Pollutantk₂₀ (m/yr)θ (–)C* (mg/L)Dominant mechanism
BOD37.51.0003.5Aerobic/anaerobic degradation
NH₄-N20.01.0600.0Nitrification (aerobic)
NO₃-N30.01.0900.0Denitrification (anoxic)
PO₄-P11.01.0100.02Sorption, precipitation
TSS50.01.0002.0Sedimentation, filtration
E.colikd=1.5/d1.0500.0UV, dessication, predation
Worked example: NH₄-N nitrification. k20=20.0, θ=1.060. At T=10°C: kT = 20.0 × 1.060^(10−20) = 20.0 × 1.060^(−10) = 20.0 / 1.791 ≈ 11.2 m/yr — 44% of the 20°C rate. At T=30°C: kT = 20.0 × 1.060^10 = 20.0 × 1.791 ≈ 35.8 m/yr — 79% higher than 20°C.

⚗️ Live P-k-C* Treatment Calculator

250 mg/L
36 m/yr
22 °C
kT (adjusted)
38.8
m/yr
Effluent Ce
29.4
mg/L
Removal
88.2
%
EU Target
Pollutant Removal Curve — Effluent vs. HLR
P-k-C* first-order model at T = 22 °C. Higher HLR (faster loading) → worse effluent. Select pollutant above to update.
BOD
NH₄-N
NO₃-N
PO₄-P
TSS

Monod Kinetics — Microbial Growth Rates

Within the biofilm on substrate grains, microbial reaction rates follow Monod kinetics — a saturation function where growth rate approaches a maximum as substrate (food) becomes abundant:

Monod Equation
$$\mu = \mu_{\max} \cdot \frac{S}{K_s + S}$$
$\mu$ = specific growth rate [d⁻¹]  |  $\mu_{\max}$ = maximum growth rate  |  $S$ = substrate concentration [mg/L]  |  $K_s$ = half-saturation constant [mg/L]
When $S \gg K_s$: zero-order kinetics. When $S \ll K_s$: first-order — this is why P-k-C* works at low effluent concentrations.

When substrate is abundant (S >> Ks), bacteria grow at their maximum rate — zero-order kinetics. When substrate is scarce (S << Ks), growth is proportional to substrate — first-order kinetics. This is why the P-k-C* model (first-order) works so well at low concentrations typical of polishing stages.

🧬 Biochemical Process Web — Pollutant Transformation Pathways

Interactive map of the seven major transformation processes in a constructed wetland. Hover each node for the governing equation and key organisms.

Process Hub
Aerobic
Anoxic/Anaerobic
Chemical
Physical

📝 Module 3 Quiz

1. In the P-k-C* model, what does C* represent?
A
The influent concentration before treatment
B
The concentration at which treatment becomes 50% efficient
C
The background (irreducible minimum) effluent concentration
D
The critical saturation concentration for biofilm growth
2. NH₄-N has θ = 1.060. If k₂₀ = 20.0 m/yr, what is kT at 30°C?
A
20.0 m/yr (θ = 1.000 for NH₄)
B
30.0 m/yr
C
≈ 35.8 m/yr
D
≈ 43.0 m/yr
3. Which pollutant has the highest k₂₀ rate constant (fastest removal) in constructed wetlands?
A
BOD (k₂₀ = 37.5)
B
NO₃-N (k₂₀ = 30.0)
C
NH₄-N (k₂₀ = 20.0)
D
TSS (k₂₀ = 50.0)
Did You Know

The rhizosphere — the thin film around plant roots — can hold 10⁸–10⁹ bacteria per gram of substrate. This biofilm, not the plant itself, does most of the pollutant breakdown in an AquaSai wetland; the plant's job is to feed and oxygenate it.

AquaSai Field Note

AquaSai calibrates the P-k-C* rate constant with Arrhenius temperature correction for each climate — tropical Koh Phangan runs warm and fast; a Tenerife winter runs cooler, so beds are sized with extra area to hold the same effluent target year-round.

— AquaSai treatment modelling notes
Try It LiveTest P-k-C* and Monod kinetics against your own influent numbers.
🛠️ MSR Tool📄 Modelling Article
Key Takeaways — Module 3
  • The P-k-C* model (Kadlec & Wallace) is the industry standard for wetland treatment prediction
  • C* (background concentration) sets an irreducible treatment floor — BOD 3.5, TN 1.5, TP 0.02 mg/L
  • Arrhenius correction (θ^(T−20)) is essential in variable climates — nitrification drops ~40% at 10°C
  • Monod kinetics explain substrate-limited growth; at C >> Ks, removal approaches maximum rate μmax
Botanical plate
Module 4
Plant Ecology & Selection
Learn which plants to choose, why polyculture outperforms monoculture, and how the rhizosphere creates the treatment engine within each wetland cell.
⏱ ~6 min read
🎯 Learning Objectives
  • Explain how macrophyte roots enhance treatment via radial oxygen loss (ROL)
  • Identify the five primary AquaSai polyculture species and their specific functional roles
  • Select appropriate species for HSSF, VSSF and FWS wetland stages
  • Explain why polyculture outperforms monoculture in resilience and year-round performance

Why Plants Matter

Wetland plants are not merely decorative — they are the biological engine of treatment. Their roots pump oxygen into the anoxic substrate (radial oxygen loss, ROL), creating alternating aerobic microsites where nitrifiers thrive and anaerobic zones where denitrifiers work. Their rhizosphere (root zone) surface provides up to 10× more biofilm attachment area per gram of substrate than bare gravel.

Plant uptake of nitrogen and phosphorus is real but secondary — it contributes roughly 10–25% of nutrient removal. The primary mechanism is microbial, catalysed by the conditions the plant roots create.

AquaSai Core Species — Click each plant to learn more

Cattail
Common Cattail
Typha latifolia
HSSFVSSFN removalBiomass
Canna
Canna Lily
Canna indica
HSSFP uptakeAesthetic
Water Iris
Water Iris
Iris pseudacorus
FWSMetalsAesthetic
Vetiver
Vetiver Grass
Chrysopogon zizanioides
VSSFTSSErosion
Water Hyacinth
Water Hyacinth
Eichhornia crassipes
FWSN + PManaged
Duckweed
Duckweed
Lemna minor
FWSN + PProtein feed

Polyculture: The AquaSai Principle

AquaSai systems always use polyculture — deliberately mixing 3–8 species per cell. Research and field observation confirm that polycultures consistently outperform monocultures:

🔬 Deeper Redox Zonation

Different root architectures (deep tap roots, shallow lateral mats) create more complex redox gradients than any monoculture.

🛡️ Resilience & Redundancy

Disease, pest pressure, or temperature extremes rarely impact all species simultaneously. Functional continuity is maintained.

🌸 Ecosystem Services

Diverse canopies provide wildlife habitat, aesthetic value, educational potential — integral to AquaSai's co-benefit model.

🌡️ Seasonal Complementarity

Different species peak in summer vs. winter, maintaining year-round treatment performance even in temperate climates.

Did You Know

Vetiver (Chrysopogon zizanioides) roots can reach 3–4 m deep in a single season — anchoring banks, mining nutrients and creating vast biofilm surface. AquaSai pairs it with canna, iris and reeds so root depths and bloom seasons complement rather than compete.

AquaSai Field Note

The polyculture plates shown here mirror the planting architecture AquaSai used in Nepal — a layered assemblage that keeps the bed productive across the year and doubles as a botanical amenity rather than an eyesore.

— AquaSai plant-ecology library
Try It LiveMatch species to your climate zone and treatment target.
🛠️ MSR Tool

📝 Module 4 Quiz

1. What is "Radial Oxygen Loss" (ROL) and why does it matter in a HSSF wetland?
A
Oxygen lost from above-ground stems to the atmosphere — reduces plant growth
B
Oxygen consumed by root respiration — creates anoxic zones
C
Oxygen pumped from aerenchyma into the root zone — creates aerobic microsites for nitrifiers
D
Oxygen released by algae in the polishing pond
2. What fraction of nutrient removal is typically attributed to direct plant uptake vs. microbial processes?
A
75–90% uptake, 10–25% microbial
B
10–25% uptake, 75–90% microbial
C
Exactly 50% each
D
0% uptake — plants are structural only
3. Which AquaSai plant species is particularly noted for TSS and sediment filtration in VSSF cells?
A
Eichhornia crassipes (Water Hyacinth)
B
Lemna minor (Duckweed)
C
Iris pseudacorus (Water Iris)
D
Chrysopogon zizanioides (Vetiver Grass)
Wetland Plant Species — Functional Performance Matrix
Relative effectiveness for key wetland functions (0–100 scale, expert consensus). Source: Vymazal 2013, Stefanakis 2019.
Key Takeaways — Module 4
  • Plant selection should balance treatment function, climate resilience and local biodiversity
  • Phragmites australis dominates globally but native alternatives (Typha, Juncus) offer ecological co-benefits
  • Root oxygen release (ROL) from VSSF plants creates aerobic micro-zones critical for nitrification
  • Macrophyte biomass and litter decomposition are significant internal nutrient cycles — harvest may be needed