Die Brooks Corey Funktion stellt eine Beziehung zwischen Bodenfeuchte oder Sättigung einerseits und der Saugspannung andererseits her:

$$ \frac{\theta-\theta_r}{\theta_s-\theta_r} = {\left[\frac{\psi_0}{\psi}\right]}^\lambda mit den Parametern: \begin{tabular}{ll} \theta & Bodenfeuchte \\ \theta_s & Sättigungs-Bodenfeuchte \\ \theta_r & Rest-Bodenfeuchte \\ \psi_0 & Lufteintrittspunkt \\ \psi & Saugspannung, Kapillarspannung \\ \end{tabular} $$

Die ungesättigte Leitfähigkeit wird nach Brooks-Corey durch folgende Beziehung dargestellt:

$$ \frac{K({\theta})}{K_s} = {\left[{\frac{\theta-\theta_r}{\theta_s-\theta_r}}\right]}^{3+2*\lambda} mit den Parametern: \begin{tabular}{ll} K_\theta & ungesättigte Bodenfeuchte \\ K_s & gesättigte Bodenfeuchte \\ \lambda & Porengrößenindex \\ \end{tabular} $$

| Brooks-Corey.R

# Brooks-Corey (1966)
th <- seq(0.05,0.4,0.005) # Bodenfeuchte
tr <- 0.05            # Restfeuchte
tm <- 0.40            # maximale Feuchte
Ks <- 1E-5            # hydraul. Leitfähigkeit in m/s
ip <- 0.1           # pore size index
 
Se <- (th-tr)/(tm-tr)
 
ps0<- 0.015            # Lufteintrittspunkt: cm
psa<- seq(0.015,1E+4)# Saugspannung in cm
 
S  <- (ps0/psa)^ip
 
plot(S,psa)