All Tools

Watershed Explorer

Explore the water cycle from precipitation to groundwater. Calculate water budgets using SCS curve numbers, visualize aquifer cross-sections with Darcy flow, model contaminant plume migration, and analyze well drawdown with the Theis equation.

Cross-Section View

Water TableI = 32.2R = 0.3ET = 17.5P = 50 mm/dayCN = 55 | Soil: sand

Controls

Precipitation50 mm/day
Slope5 °

Results

P=I+R+ETP = I + R + ET
50.0 = 32.2 + 0.3 + 17.5 mm/day
Infiltration
32.2 mm/day
64%
Runoff
0.3 mm/day
1%
Evapotranspiration
17.5 mm/day
35%
SCS Curve Number
55

Water Budget Components

Reference Guide

Darcy's Law

Darcy's Law describes the flow of fluid through a porous medium. The specific discharge (Darcy velocity) is proportional to the hydraulic gradient.

v=Kdhdlv = -K \frac{dh}{dl}

Where KK is hydraulic conductivity (m/day), and dh/dldh/dl is the hydraulic gradient. The actual pore velocity is vs=v/nv_s = v/n where nn is porosity.

Water Budget

The water budget partitions precipitation into infiltration, surface runoff, and evapotranspiration.

P=I+R+ETP = I + R + ET

Runoff is estimated using the SCS Curve Number method, where S=25400CN254S = \frac{25400}{CN} - 254 is the maximum soil retention and Q=(P0.2S)2P+0.8SQ = \frac{(P - 0.2S)^2}{P + 0.8S} for P>0.2SP > 0.2S.

Theis Equation

The Theis equation gives the drawdown at distance r from a pumping well in a confined aquifer.

s=Q4πTW(u),u=r2S4Tts = \frac{Q}{4\pi T} W(u), \quad u = \frac{r^2 S}{4Tt}

For small uu (large time or small distance), the Cooper-Jacob approximation gives W(u)0.5772ln(u)W(u) \approx -0.5772 - \ln(u).

Contaminant Dispersion

A contaminant plume spreads by advection (bulk flow) and dispersion (mixing). The 2D Gaussian solution describes concentration in space and time.

C(x,y,t)=C04πtDLDTe(xvt)24DLty24DTtC(x,y,t) = \frac{C_0}{4\pi t \sqrt{D_L D_T}}\, e^{-\frac{(x-vt)^2}{4D_Lt} - \frac{y^2}{4D_Tt}}

The dispersion coefficients are DL=αLvD_L = \alpha_L v (longitudinal) and DT0.1DLD_T \approx 0.1 D_L (transverse), where αL\alpha_L is dispersivity.