When sewage or industrial wastewater enters a river, microbes oxidize the organic matter and drain dissolved oxygen (DO) from the water. The Streeter-Phelps equation balances deoxygenation against atmospheric reaeration to locate the downstream critical point, telling you whether a discharge load will still meet water-quality standards.
Parameters
Ultimate BOD L₀
mg/L
Ultimate (final) BOD just downstream of the outfall — the discharge load
Saturated DO D_s
mg/L
Equilibrium DO from temperature and altitude (about 9.1 mg/L at 20 °C)
Initial DO deficit D₀
mg/L
DO deficit right after mixing (D_s − measured DO) at the outfall
Deoxygenation rate k_d
1/day
First-order rate of microbial BOD decay — typically 0.1–0.5 at 20 °C
Reaeration rate k_r
1/day
First-order atmospheric reaeration rate — set by velocity, depth and turbulence
Flow velocity v
km/day
Mean velocity converting elapsed time t into downstream distance x = v·t
Results
—
Critical time t_c (day)
—
Critical distance x_c (km)
—
Critical DO deficit D_c (mg/L)
—
Critical DO (mg/L)
—
DO at 1 day (mg/L)
—
Water quality verdict
—
River cross-section view — downstream DO sag
Upstream is on the left, downstream on the right. Colour shading shows DO concentration (deep blue = high / pale = low). A fish marker sits at the critical point.
DO sag curve — downstream distance vs DO / deficit
Downstream DO deficit D(t) [mg/L]. L₀: ultimate BOD, D₀: initial deficit, k_d: deoxygenation rate, k_r: reaeration rate. Actual DO follows DO = D_s − D.
Critical time t_c and the matching downstream critical distance x_c. v: mean velocity [km/day]. The balance between k_d and k_r sets where the sag bottom occurs.
Critical deficit and critical DO, derived from dD/dt = 0. Keeping DO_c above the regulatory standard (≈ 5 mg/L for warm-water fisheries) is the core requirement of any discharge permit.
River DO/BOD model (Streeter-Phelps)
🙋
Why do fish die when sewage enters a river? Are they killed by the dirt itself?
🎓
In the vast majority of cases it's actually oxygen starvation, not direct toxicity. Microbes in the water consume the organic matter (proteins, sugars, fats) in the wastewater and breathe up huge amounts of dissolved oxygen in the process — that's "deoxygenation". Meanwhile oxygen slowly seeps into the water from the atmosphere, called "reaeration". Wherever deoxygenation outruns reaeration, the river's DO keeps falling until it hits a minimum — the "critical point". Below about 2 mg/L most fish can no longer breathe and die.
🙋
So Streeter-Phelps is what predicts that drop? It's from 1925 — is a hundred-year-old equation still good enough?
🎓
Yes — Streeter and Phelps fit it to Ohio River data back then, and it's the textbook classic of water-quality engineering. But it still shows up first in every wastewater permit application and river standard study. Why? Because it has a closed-form analytic solution, you can see at a glance which parameter matters, and it runs instantly so it's perfect for sensitivity work. For serious projects we'd use a numerical model like QUAL2K or WASP, but engineers still use Streeter-Phelps as a sanity check on those outputs. It's also unbeatable for teaching, because the deoxygenation-versus-reaeration physics is right on the surface.
🙋
If I push the ultimate BOD L₀ slider from 30 to 100, the critical DO goes negative. How am I supposed to read that?
🎓
Good catch. Streeter-Phelps happily computes D = D_s − DO even when the load is too high, but physical DO can never go below zero. So a negative critical-DO prediction is really a "this river will go anaerobic" warning. Under anaerobic conditions a different microbial community kicks in: sulfate reducers release hydrogen sulfide (that classic rotten-egg sewer smell), methanogens make methane, and the system shifts into a regime the simple model can't describe. For engineers it's a hard red flag — BOD load must come down or treatment must be upgraded. This tool flags any critical DO below 2 mg/L as "severely polluted / anoxic risk" for that reason.
🙋
When I change velocity v, the critical time stays the same but the critical distance moves. Does that mean steep mountain streams and slow big rivers behave very differently?
🎓
Exactly. Flow velocity v only rescales time t into downstream distance x = v·t, so the chemistry itself is unchanged. But in reality, fast turbulent streams also have much larger reaeration k_r — the O'Connor-Dobbins formula gives k_r ≈ 3.93·v^0.5·H^(-1.5), with depth H. So a fast river has its critical point far downstream but only a shallow dip, while a slow lowland river has the dip close to the outfall but goes much deeper. The classic Japanese pattern is that headwater streams self-purify easily and downstream urban rivers like the lower Tama River suffer the worst sags.
🙋
So when we design a wastewater plant, we use this tool to figure out how low L₀ must be?
🎓
In practice it's one more step: first you do a mass-balance mixing calculation using wastewater flow Q_w · BOD C_w and river flow Q_r · BOD C_r to get the post-mix L₀ and D₀. You feed those into a Streeter-Phelps run and check that the predicted critical DO stays above the environmental standard (in Japan, 7.5 mg/L for class A streams, 5 mg/L for class B). If it doesn't, you compare costs of upgrading treatment (activated sludge → A2O → tertiary), of relocating the outfall to a more turbulent reach, or of converting combined sewers to separate systems. A surprisingly cheap fix is sometimes "discharge into a fast, shallow riffle" — the extra reaeration can let you skip a whole treatment stage.
Frequently Asked Questions
Published by Streeter and Phelps in 1925, this equation gives the analytical solution for how dissolved oxygen (DO) changes downstream of a wastewater outfall. It combines two first-order reactions: oxygen consumption by microbial decomposition of organic matter (deoxygenation) and atmospheric reaeration at the water surface, yielding the closed form D(t) = (k_d·L₀)/(k_r−k_d)·(exp(−k_d·t)−exp(−k_r·t)) + D₀·exp(−k_r·t). Despite being simple, it is still the first tool reached for in wastewater discharge permitting, allowable-load decisions and river quality assessment, and forms the oxygen core of advanced models such as QUAL2K.
The DO deficit D(t) first grows because deoxygenation dominates, then turns around when reaeration catches up. The maximum is the critical point at t_c = ln[(k_r/k_d)·{1 − D₀(k_r−k_d)/(k_d·L₀)}]/(k_r−k_d). DO at this point (critical DO) is the lowest oxygen value anywhere along the river, so it decides whether fish can live or the channel goes anaerobic and produces hydrogen sulfide and methane. Wastewater discharge permits are typically issued only when the predicted critical DO stays above the relevant environmental standard (e.g. 5 mg/L).
At 20 °C, deoxygenation rate k_d is roughly 0.1–0.5 day⁻¹ (typically 0.2–0.4 for urban rivers dominated by domestic sewage). Reaeration k_r depends strongly on hydraulics: 0.2–0.5 day⁻¹ for slow large rivers, 0.5–1.5 for moderately sloped streams, and 2–10 for steep reaches with riffles or waterfalls. As long as k_r > k_d the DO will recover downstream; if k_r < k_d the river cannot self-purify and stays oxygen-depleted. Temperature shifts both rates upward (θ=1.047 for k_d, 1.024 for k_r) while reducing saturated DO, which is why summer is the worst-case design condition.
Real rivers also experience sediment oxygen demand (SOD), algal photosynthesis and night-time respiration, additional oxygen draw from nitrification, flow-dependent reaeration (O'Connor-Dobbins k_r = 3.93·v^0.5·H^(-1.5)) and dilution by tributaries. The bare two-term Streeter-Phelps model is therefore not accurate enough for detailed design. It remains, however, the standard first move for (1) near-field screening at outfalls, (2) sanity-checking complex numerical models, and (3) teaching. Production design uses QUAL2K, WASP or MIKE11, but Streeter-Phelps is still how engineers confirm the results are physically plausible.
Real-world applications
Wastewater treatment plant permits and tertiary upgrades: Most national water-quality laws require that the post-mixing critical DO meet the environmental standard for the receiving stream. In Japan that means 7.5 mg/L for class A waters (upstream drinking-water sources), 5 mg/L for class B (fishable) and 2 mg/L for class C; the EU Water Framework Directive and US EPA criteria follow similar logic. When a Streeter-Phelps run predicts a critical DO below the target, the plant must add anaerobic-anoxic-oxic (A2O) stages, membrane bioreactors (MBR) or ozonation to cut BOD, nitrogen and phosphorus further before discharge.
Combined sewer overflow (CSO) impact assessment: Older cities combine stormwater and sewage in a single pipe, so heavy rain forces untreated sewage straight into the river. CSO events push L₀ above 100 mg/L for hours, and Streeter-Phelps runs reliably predict a critical DO near zero — the mechanism behind classic fish kills in Tokyo Bay tributaries and the historic Thames and Seine. Cities now build deep CSO storage tunnels (the Tokyo Metropolitan Outer Underground Discharge Channel, London's Thames Tideway) and gradually convert combined systems to separated networks, with DO modelling justifying the spend.
Industrial discharge allocation: Food, pulp & paper, brewery and livestock industries produce high-BOD effluent and must work out a permissible discharge load given the river's flow and existing background load. Streeter-Phelps lets the engineer simulate "current critical DO + my new discharge" and back-calculate the largest L₀ that still meets the standard. The same workflow underpins new plant siting, ISO 14001 environmental impact assessments and expansion permits.
Urban river restoration and visualising self-purification: River revival projects such as Cheonggyecheon in Seoul, the Sumida in Tokyo and the Thames in London deliberately raise the reaeration coefficient k_r through dredging, increased flow and engineered riffles. Showing the public a Streeter-Phelps curve where only k_r changes — and DO 10 km downstream rises by 2 mg/L thanks to a small weir under a bridge — is one of the most persuasive ways to argue for funding.
Common misconceptions and pitfalls
The first misconception is to equate "BOD below the limit" with "water quality is fine". BOD is a prediction of how much oxygen the organic load will eventually consume; it is not the current DO. Even at modest BOD, reaches with low velocity and weak reaeration, or thick sediment beds with high SOD, can run below the DO standard. BOD is a load metric; DO is what actually decides whether fish survive. Design judgement must be made on DO, which is the entire point of this tool.
The second pitfall is to ignore temperature and season. The tool lets you enter k_d, k_r and saturated DO separately, but in nature they are linked by temperature. At 25 °C in mid-summer k_d is about 1.3× its 20 °C value, k_r rises slightly, and saturated DO drops from 9.1 to 8.4 mg/L. The critical DO drops dramatically, which is the basic reason fish kills cluster in July and August every year on stressed rivers — and the same pattern shows up on the Ebro in Spain or southern US rivers in drought summers. Always evaluate the critical DO at peak summer temperatures, not at the annual average.
The third trap is to read negative DO predictions as "the equation is broken". Real DO cannot go below zero, but Streeter-Phelps is a linear model that simply asks "how much oxygen would be missing if supply were unlimited". Negative DO output is therefore not a bug — it is a flag that the river will actually go anaerobic and a different chemistry will take over. Treat it as an automatic "design fails, anoxic risk, cut the load" verdict. This tool already warns you when critical DO drops below 2 mg/L; for rigorous work the model should be capped at DO = 0 and an anaerobic decomposition model substituted beyond that point.
How to Use
Enter ultimate BOD (mg/L) — typically 150–300 mg/L for municipal sewage, 500–2000 mg/L for food processing waste
Set saturated DO (mg/L) — depends on water temperature; use 10–12 mg/L at 5°C, 8–9 mg/L at 20°C, 6–7 mg/L at 30°C
Input initial DO deficit (mg/L) — the oxygen shortfall immediately downstream of discharge
Specify deoxygenation rate constant (k₁ day⁻¹) — typically 0.1–0.5 for treated effluent, 0.3–1.0 for raw sewage
Click Calculate to find critical time, critical distance, minimum DO concentration, and water quality classification
Worked Example
Municipal WWTP discharge into a river: ultimate BOD = 180 mg/L, saturated DO = 9 mg/L, initial DO deficit = 7 mg/L, deoxygenation rate k₁ = 0.35 day⁻¹. The simulator predicts critical time t_c = 2.8 days at critical distance x_c ≈ 8.4 km (assuming 3 km/day velocity). Minimum DO drops to 2.1 mg/L, triggering a "Poor" water quality verdict. After 1 day, DO recovers to 4.3 mg/L as organic matter depletes and reaeration occurs.
Practical Notes
Critical DO often falls below 5 mg/L threshold for fish survival; check sewage treatment plant pre-discharge BOD removal efficiency
Cold-water rivers (5°C) recover slower than warm rivers because saturation DO is higher but reaeration rate k₂ is lower
Industrial discharges (tannery, brewery, slaughterhouse) typically require BOD pre-treatment to prevent anoxic zones and fish kill events
Use 24-hour BOD (BOD₅) measured value multiplied by 1.15–1.25 to estimate ultimate BOD for this model