Dissolved Oxygen Sag Curve Calculator

The Streeter-Phelps equation is a first-order kinetic model that predicts dissolved oxygen (DO) depletion and recovery in streams following organic pollutant discharge. It balances deoxygenation (BOD-driven oxygen consumption) and reaeration (atmospheric oxygen transfer), yielding a DO sag curve. Despite its simplifications—such as assuming steady flow, constant temperature, and single-point BOD loading—it remains widely adopted due to its transparency, regulatory acceptance (e.g., US EPA WQX documentation, ISO 5667-14:2022 for water quality assessment), and utility in screening-level TMDL development and preliminary permit evaluations. Its analytical solution enables rapid scenario testing, making it foundational for training, compliance checks, and integration into more complex models like QUAL2K as a benchmark.

Select k_d and k_a using site-specific calibration where possible. k_d typically ranges from 0.1–0.5 day⁻¹ at 20°C and depends on wastewater composition, temperature, and microbial activity; apply the Arrhenius correction (θ = 1.047) for non-standard temperatures per ASTM D5259-18. k_a is strongly influenced by turbulence, depth, and velocity—empirical formulas (e.g., O’Connor–Dobbins, Churchill) or field tracer studies are preferred over textbook defaults. US EPA’s WASP and HEC-RAS recommend deriving k_a from hydraulic geometry and wind data. Default values (e.g., k_d = 0.2, k_a = 0.3 day⁻¹) should only serve as initial estimates and must be validated against DO profile measurements to avoid underestimating critical deficit or misplacing the critical point.

The classic Streeter-Phelps model assumes a single instantaneous BOD load and steady, uniform flow—so direct application to systems with multiple discharges or unsteady flow violates its core assumptions and risks significant error. For such cases, use segmented modeling: compute sequential sag curves with mass-balanced L₀ and DO₀ at each confluence, or adopt numerical tools like QUAL2K or CE-QUAL-W2 that support dynamic loading and variable hydraulics. Regulatory guidance (e.g., US EPA Technical Support Document for TMDLs, 2021) explicitly requires superposition or distributed loading approaches when >1 major outfall exists. Always document assumptions and perform sensitivity analysis on flow-weighted BOD inputs to quantify uncertainty in the predicted minimum DO.

For coldwater fisheries (e.g., trout), the U.S. EPA recommends a minimum 7-day average DO ≥ 6.0 mg/L, with no measurement below 5.0 mg/L (40 CFR Part 131; also reflected in state standards like Oregon DEQ OAR 340-041-0027). The Streeter-Phelps sag curve identifies the critical distance (x_c) and corresponding minimum DO (DO_min), which must exceed these thresholds to demonstrate compliance. If DO_min falls below criteria, the model helps prioritize mitigation—e.g., reducing L₀ via pretreatment or enhancing k_a via aeration weirs. Note: Some jurisdictions (e.g., California Water Code § 13050) require modeling at design low-flow (7Q10) conditions, not mean flow, to ensure protective worst-case evaluation.

Temperature critically affects both DO_s and reaction rates. DO_s decreases nonlinearly with rising temperature (e.g., from 14.6 mg/L at 0°C to 7.6 mg/L at 30°C); use the APHA Standard Methods 2540-C equation or USGS coefficients—not fixed defaults—to compute DO_s accurately. Simultaneously, k_d and k_a increase with temperature per the Arrhenius relationship (θ ≈ 1.047 for k_d, 1.024 for k_a). Using a static DO_s (e.g., 9.1 mg/L at 20°C) while modeling summer conditions will overpredict DO and underestimate deficit severity. Always pair temperature-corrected DO_s with adjusted rate constants—and validate with field DO profiles across seasons per ISO 5667-3:2012 sampling protocols.

Yes—the Streeter-Phelps model is explicitly endorsed for preliminary and screening-level analyses in key regulatory frameworks, including US EPA’s TMDL Technical Support Document (2021) and many state NPDES permit technical guidance manuals (e.g., Texas TCEQ RG-422). However, its use is often conditional: it must be applied at design low-flow (e.g., 7Q10), calibrated to local data, and supplemented with uncertainty analysis. For complex systems or contested permits, agencies increasingly require verification via calibrated numerical models (e.g., WASP, MIKE HYDRO River) or field monitoring. Documentation must cite assumptions, data sources, and validation metrics per ASTM D5259-18 to satisfy peer review requirements in formal submissions.

A mathematically predicted critical point upstream indicates model instability caused by violating the fundamental assumption that k_a > k_d. When k_a ≤ k_d, the reaeration cannot offset deoxygenation, leading to monotonic DO decline and no true minimum downstream—hence the spurious upstream solution. This commonly arises from using overly conservative k_a (e.g., 0.1 day⁻¹ in deep, slow rivers) or inflated k_d (e.g., >0.4 day⁻¹ without temperature correction). Remediate by recalibrating k_a using hydraulic data (e.g., stream power or shear velocity methods per USGS Techniques and Methods 4-A6), verifying BOD₅-to-L₀ conversion (typically L₀ = BOD₅ / (1 − e^(−k_d·5))), and ensuring units and temperature corrections are consistent. Field DO profiling is essential to confirm physical plausibility.