Industrial Effluent Mixing Zone Evaluation for Pulp Mill in Northern Ontario

Case Study 2: Industrial Effluent Mixing Zone Evaluation for Pulp Mill in Northern Ontario

Scenario Project Type: Environmental impact assessment (EIA) for mill expansion and updated effluent discharge permit application. Location Context: Remote boreal forest tributary of the Mattagami River, Ontario, Canada — low-gradient, cold-water system (summer avg. 14°C), high natural organic load, and Indigenous fishing rights considerations. DO naturally fluctuates between 9.0–11.2 mg/L due to cool temperatures and forest canopy shading. Constraints: Limited historical hydrological data; Indigenous knowledge indicated fish avoidance near discharge during late summer; provincial regulation requires DO ≥ 6.5 mg/L at all points within 5 km mixing zone; winter ice cover precludes field validation until spring.

Given Data

• Initial dissolved oxygen (DO₀): 10.3 mg/L (pre-discharge ambient, measured in May)
• Initial BOD (L₀): 42.0 mg/L (combined mill effluent + stormwater runoff BOD₅, conservative estimate)
• Deoxygenation rate constant (k_d): 0.12 day⁻¹ (adjusted for 14°C using Arrhenius correction from lab 20°C k_d = 0.20)
• Reaeration rate constant (k_a): 0.22 day⁻¹ (low turbulence, shallow riffle-pool morphology)
• Saturation DO (DOₛ): 10.0 mg/L (14°C, 95% saturation)
• Distance downstream (x): 5 km

Calculation Using the same Streeter–Phelps framework and calibrated travel time: Average velocity estimated from bathymetry and flow gauging (Q = 4.8 m³/s, mean depth = 1.1 m, width = 12 m) → v ≈ 0.36 m/s ≈ 31.1 km/day → t = 5 km / 31.1 km/day ≈ 0.1608 days

Compute:

• e^(−k_a·t) = e^(−0.22 × 0.1608) ≈ e^(−0.0354) ≈ 0.965
• e^(−k_d·t) = e^(−0.12 × 0.1608) ≈ e^(−0.0193) ≈ 0.981
• (DOₛ − DO₀) = 10.0 − 10.3 = −0.3
• (k_d·L₀)/(k_a − k_d) = (0.12 × 42.0)/(0.22 − 0.12) = 5.04 / 0.10 = 50.4
• Term1 = (−0.3) × 0.965 ≈ −0.289
• Term2 = 50.4 × (0.981 − 0.965) = 50.4 × 0.016 ≈ 0.806
• DO(5 km) = 10.0 − (−0.289) + 0.806 = 10.0 + 0.289 + 0.806 ≈ 11.095 → 11.10 mg/L

However, the minimum DO (sag point) occurs at t_max = ln(k_a/k_d)/(k_a − k_d) = ln(0.22/0.12)/(0.10) ≈ ln(1.833)/0.10 ≈ 0.607/0.10 ≈ 6.07 days → ~189 km downstream — far beyond the 5 km mixing zone. Thus, DO remains near saturation throughout the zone.

Result and Decision Predicted DO at 5 km = 11.10 mg/L, exceeding the 6.5 mg/L requirement by >4 mg/L. Yet community interviews and benthic macroinvertebrate surveys revealed elevated sulfide odors and low chironomid diversity within 1 km — indicating anaerobic microzones not captured by bulk BOD/DO modeling. Engineers concluded that while bulk DO compliance was satisfied, localized oxygen depletion near sediment interface required mitigation. They mandated installation of diffused aeration at the discharge pipe and revised effluent limits to include maximum sulfide and total suspended solids (TSS) — addressing root causes of benthic hypoxia.

Lesson Bulk DO modeling validates regulatory compliance but may overlook microhabitat-scale deoxygenation — integrate biological field indicators (e.g., benthic assemblages, redox potential) to detect hidden stressors beneath acceptable DO averages.