Soil Physics with HYDRUS: Chapter 3

41
Soil Physics with HYDRUS: Chapter 3 Equations, Tables, and Figures

description

Soil Physics with HYDRUS: Chapter 3. Equations, Tables, and Figures. 3.1 – 3.5. 3.6 – 3.9. 3.10 – 3.13. 3.14 – 3.17. 3.18 – 3.20. 3.21 – 3.25. 3.26 – 3.30. 3.31 – 3.34. 3.35 – 3.39. 3.40. - PowerPoint PPT Presentation

Transcript of Soil Physics with HYDRUS: Chapter 3

Page 1: Soil Physics with HYDRUS: Chapter 3

Soil Physics with HYDRUS: Chapter 3

Equations, Tables, and Figures

Page 2: Soil Physics with HYDRUS: Chapter 3

3.1 – 3.54

8wR g HQL

LHK

AQ

s

w sHJ KL

dzdHKJ sw

zHKJ sw

Page 3: Soil Physics with HYDRUS: Chapter 3

3.6 – 3.92 1

w s sH H b LJ K KL L

2 1w s s

H H h zJ K KL z

s sh z b LK Kz L

bh zL

gk

K ls

Page 4: Soil Physics with HYDRUS: Chapter 3

3.10 – 3.13

1

1

N

jj

eff Nj

j j

LK

LK

w effHJ KL

1

1

N

j jj

eff N

jj

L KK

L

x sdzJ Kdx

Page 5: Soil Physics with HYDRUS: Chapter 3

3.14 – 3.172

2max 0

s

Rsz zK

( )( ) ( ) ( )

( ) 1

wH h z h zJ K h K h K hz z z zhK hz

( ) hsK h K e

( )m

ss

K K

Page 6: Soil Physics with HYDRUS: Chapter 3

3.18 – 3.20( )

1

sN

KK hha

2

2

( )01

( )0

( )

e

e

e

e

e

SdSh S

le s e

dSh S

K S K S

2

( )0

21

( )0

( )

e

e

e

e

e

SdSh S

le s e

dSh S

K S K S

Page 7: Soil Physics with HYDRUS: Chapter 3

3.21 – 3.251 2/( ) l

e s eK S K S

21/( ) 1 (1 )l m me s e eK S K S S

1 1/ , 1m n n

1 1 2 2

1 2 1 2

21/ 1/

1 2 1 1 2 2

21 1 2 2

( ) 1 (1 ) 1 (1 )( )

( )

m m m mle e e e

e s

w S w S w S w SK S K

w w

0

0

( )

( ) ( )i

h

hc

i

K h dh

K h K h

Page 8: Soil Physics with HYDRUS: Chapter 3

3.26 – 3.301

c

1 2 3 4 5 2 1( ) 4 2 4 2 ... 2 43

b

n n na

xy x dx y y y y y y y y

2 3

( )

( )

s a

as a

K h K h h

hK h K h h

h

wJ K h2

1

1 2

1( )

h

wh

dh z zJK h

Page 9: Soil Physics with HYDRUS: Chapter 3

3.31 – 3.34 1 / /

1 tanh ( )s ssi K i KKh a L z

i a

0

04( ) 1

s w z

c

i J

K hr

22 c

s sH HQ K rG G

2

2 21 cs s

H Hi KG r G r

Page 10: Soil Physics with HYDRUS: Chapter 3

3.35 – 3.392 3

1 2 3 41

2H H HG A A A Ar r r

1 2

ws

J LKH H

2 1w s

H HdbJ Kdt L

0

1 1

lnsb LLK

t b L

2

2s sG

Hi KG r

Page 11: Soil Physics with HYDRUS: Chapter 3

3.40

21sinh 1

2G

H r rr H H

G

Page 12: Soil Physics with HYDRUS: Chapter 3

TABLE 3.1Brooks and Corey (1964) soil hydraulic parameters (Ks, θr, θs, ha, and λ) for 11 textural classes (Rawls et al., 1982). Macroscopic capillary length (λc) calculated using Equation 3.21 and Equation 3.25.

Textural class Ks θr θs ha λ λc

(cm d-1)

(cm3 cm-3) (cm3 cm-3) (cm) (-) (cm)

Sand 504.0 0.020 0.417 -7.26 0.592 9.88

Loamy sand 146.6 0.035 0.401 -8.69 0.474 12.28

loam 62.16 0.041 0.412 -14.7 0.322 22.18

Loam 31.68 0.027 0.434 -11.1 0.220 17.79

Silt loam 16.32 0.015 0.486 -20.7 0.211 33.38

Sandy clay loam 10.32 0.068 0.330 -28.1 0.250 44.16

Clay loam 5.52 0.075 0.390 -25.9 0.194 42.27

Silty clay loam 3.60 0.040 0.432 -32.6 0.151 55.03

clay 2.88 0.109 0.321 -29.2 0.168 48.61

Silty clay 2.16 0.056 0.423 -34.2 0.127 58.96

Clay 1.44 0.090 0.385 -37.3 0.131 64.07

Page 13: Soil Physics with HYDRUS: Chapter 3

TABLE 3.2Soil properties (texture, sand, silt, clay, structure, and Ks) for the Cecil loamy sand (plot 4 from Bruce et al., 1983).

Horizon Texture Clay Silt Sand Structure Ks

% % % cm h-1

Ap Loamy sand 7 15 78 Weak, medium granular 19.19BA Clay loam 37 20 43 Weak, medium,

subangular blocky7.69

Bt1 Clay 50 20 30 Strong, medium, subangular blocky

10.73

Bt2 Clay 41 25 34 Strong, medium, subangular blocky

0.206

BC Clay loam 36 27 37 Weak, medium, subangular blocky

0.035

C Sandy clay loam

24 24 52 Massive 0.467

Page 14: Soil Physics with HYDRUS: Chapter 3

TABLE 3.3Soil texture/structure categories for estimation of macroscopic capillary length (λc), adapted from Elrick and Reynolds (1992).

Soil texture/structure category λc

(cm)Coarse and gravelly sands; may also include some highly structured soils with large cracks and/or macropores

2.8

Most structured soils from clays through loams; also includes unstructured medium and fine sands

8.3

Soils which are both fine textured (clayey) and unstructured 25Compacted, structureless, clayey materials such as landfill caps and liners, lacustrine or marine sediments, etc

100

Page 15: Soil Physics with HYDRUS: Chapter 3

TABLE 3.4Dimensionless coefficients for the polynomial (Equation 3.35) describing the geometric factor G, valid for H/r < 10 (Bosch and West, 1998)

Soil texture/structure A1 A2 A3 A4

Sand 0.079 0.516 -0.048 0.002Structured loams and clays 0.083 0.514 -0.053 0.002Unstructured clays 0.094 0.489 -0.053 0.002

Page 16: Soil Physics with HYDRUS: Chapter 3

Figure 3.1 A vertical column of soil 40.0 cm in length with 5.0 cm of water ponded on the surface and water allowed to drip from the bottom. The soil saturated hydraulic conductivity is 10.0 cm h-1

Page 17: Soil Physics with HYDRUS: Chapter 3

Figure 3.2 A vertical soil column of length L and saturated hydraulic conductivity Ks. Water is ponded to a depth b at the top of the column and water discharges freely from the bottom at atmospheric pressure.

Page 18: Soil Physics with HYDRUS: Chapter 3

Figure 3.3 A vertical soil column with an arbitrary point at height z inside the column, instead of the top of the column.

Page 19: Soil Physics with HYDRUS: Chapter 3

0.01

0.1

1

10

100

1000

10000

100000

0.001 0.01 0.1 1

Clay Fraction (g/g)

Satu

rate

d Hy

drau

lic C

ondu

ctiv

ity (c

m/d

)

Figure 3.4 Saturated hydraulic conductivity as a function of clay fraction for 324 soils from the UNSODA database (Nemes et al., 2001).

Page 20: Soil Physics with HYDRUS: Chapter 3

Figure 3.5 Saturated hydraulic conductivity as a function of depth in a Cecil series soil from North Carolina Piedmont (Schoneberger and Amoozegar, 1990).

Page 21: Soil Physics with HYDRUS: Chapter 3

Figure 3.6 A soil column with N layers of thickness L1...LN and saturated hydraulic conductivites K1...KN on the left side. An equivalent uniform column on the right with a saturated hydraulic conductivity Keff.

Page 22: Soil Physics with HYDRUS: Chapter 3

Figure 3.7 A column consisting of a sand layer 30.0-cm thick overlying a loam layer 20.0-cm thick. The saturated hydraulic conductivities of the sand and loam layers are 21.0 and 1.32 cm h-1, respectively. Water is ponded at the top of the column to a depth of 5.0 cm and allowed to drip from the bottom of the column.

Page 23: Soil Physics with HYDRUS: Chapter 3

Figure 3.8 Geometry considered for defining the drainage equation. Water height in the ditch or drain is z0 (Warrick, 2003).

Page 24: Soil Physics with HYDRUS: Chapter 3

Figure 3.9 Examples of unsaturated hydraulic conductivity curves as a function of pressure head (top) and volumetric water content (bottom) for the sand, loam, and clay soil textural classes based on parameters in Table 3.1.

Page 25: Soil Physics with HYDRUS: Chapter 3

Figure 3.10 Boreholes and wetting fronts in a soil with a small λc and a soil with a large λc.

Page 26: Soil Physics with HYDRUS: Chapter 3

Figure 3.11 Excel spreadsheet with Brooks and Corey parameters for a loamy sand from Table 3.1.

Page 27: Soil Physics with HYDRUS: Chapter 3

Figure 3.12 The sum of the Simpson terms in column C is calculated in cell C18. The integral of K(h) is calculated in cell C19 by multiplying the sum by 10/3 (h/3) (see Equation 3.27). In cell C20, λc is calculated by dividing the integral by the difference between Ks and K(-100) (see Equation 3.25).

Page 28: Soil Physics with HYDRUS: Chapter 3

Figure 3.13 Distribution of h as a function of z for a steady low infiltration rate (0.5 cm d -1) and a high infiltration rate (1.5 cm d-1) based on Equation 3.31 (adapted from Jury and Horton, 2004).

-100

-90

-80

-70

-60

-50

-40

-30

-20

-10

0

-40 -30 -20 -10 0

Pressure Head (cm)

Dep

th (c

m)

i = 0.5 cm/di = 1.5 cm/d

Page 29: Soil Physics with HYDRUS: Chapter 3

Figure 3.14 Steady 3-D infiltration of water from a ring of radius r at the soil surface.

Page 30: Soil Physics with HYDRUS: Chapter 3

Figure 3.15 Steady 3-D infiltration of water from a borehole.

Page 31: Soil Physics with HYDRUS: Chapter 3

Figure 3.16 New Project dialog window in HYDRUS-1D.

Page 32: Soil Physics with HYDRUS: Chapter 3

Figure 3.17 Project panel in HYDRUS-1D.

Page 33: Soil Physics with HYDRUS: Chapter 3

Figure 3.18 Geometry Information dialog window in HYDRUS-1D.

Page 34: Soil Physics with HYDRUS: Chapter 3

Figure 3.19 Time Information dialog window in HYDRUS-1D.

Page 35: Soil Physics with HYDRUS: Chapter 3

Figure 3.20 Soil Hydraulic Property Model dialog window in HYDRUS-1D.

Page 36: Soil Physics with HYDRUS: Chapter 3

Figure 3.21 Water Flow Parameters dialog window in HYDRUS-1D.

Page 37: Soil Physics with HYDRUS: Chapter 3

Figure 3.22 Rosetta Lite dialog window in HYDRUS-1D.

Page 38: Soil Physics with HYDRUS: Chapter 3

Figure 3.23 Predicted parameters for a loamy sand textural class in the output panel.

Page 39: Soil Physics with HYDRUS: Chapter 3

Figure 3.24 Section of capillary tube of radius R and length L filled with water flowing in response to a pressure difference P2 - P1 (Jury and Horton, 2004).

Page 40: Soil Physics with HYDRUS: Chapter 3

0

0.2

0.4

0.6

0.8

1

0 0.2 0.4 0.6 0.8 1 1.2

V(r)/Vmax

r/R

Figure 3.25 Parabolic distribution of velocities within a cylinder. The velocity is zero at the wall where r = R and at a maximum = Vmax at the center of the cylinder where the r = 0 (Jury and Horton, 2004).

Page 41: Soil Physics with HYDRUS: Chapter 3

Figure 3.26 Excel spreadsheet for calculating pressure head distribution under steady flow in a profile of the Cecil soil with properties from Table 2.4 and Table 3.2 . Water is ponded to a depth of 0.1 cm at the surface and there are tile drains at a depth of 250 cm (pressure head equals zero).