🌡️ Thermal Front Analysis
ΔT = T_ocean - T_river
Thermal Front Strength = ΔT × (ρ_ocean - ρ_river) / ρ_avg
Mixed Layer Depth = 2.5 × u*² / (N² × f)
🦠 Primary Productivity
PP = P_max × (I/I_k) × exp(1 - I/I_k) × f(N)
f(N) = [NO₃] / (K_s + [NO₃])
Light Attenuation: I_z = I_0 × exp(-k × z)
🏔️ Sedimentation Dynamics
Settling Velocity: w_s = (2gr²(ρ_s - ρ_w)) / (9μ)
Deposition Rate = C × w_s × (1 - τ_b/τ_cr)
Resuspension = M × (τ_b/τ_cr - 1)
💨 Oxygen Dynamics
dO₂/dt = k_a(O₂_sat - O₂) + PP - R - SOD
O₂_sat = f(T, S, P)
Reaeration: k_a = 0.31 × U₁₀^1.24 × Sc^(-0.5)
🧪 Nutrient Cycling
dN/dt = -μ_max × N/(K_s + N) × B + R_N + D_N
Mixing Coefficient: K_z = κu*z(1 - z/h)
Advection: u∂N/∂x + v∂N/∂y + w∂N/∂z
🌊 Physical Mixing
Richardson Number: Ri = N²/(∂u/∂z)²
Stratification: N² = -(g/ρ₀)(∂ρ/∂z)
Wind Mixed Layer: h_ml = 2.5u*/(f × N)