<br />where Lp is the L value at t; ty' Because L in Eq, 47
<br />is not expressed as an explicit function of t, the value
<br />of L for given t must be determined numerically, for
<br />instance, by the Newton-Raphson method, After L is
<br />computed, the value of L will be substituted either
<br />into Eq, 38 or 44 for the computation of the f value,
<br />Note that use of Eq. 44 has a slight advantage over
<br />that of Eq, 38 because any enor involved in the
<br />evaluation of dL/dt in Eq, 38 can be avoided,
<br />However, use of either Eq, 38 or 44 requires the prior
<br />computation of the values of tp and Lp'
<br />
<br />Determination of the tp and Lp values, Before
<br />ponding (t<O;tp), the cumulative inmtration, F(t), is rt
<br />for constant rainfall intensity, r. Thus, from Eq, 38,
<br />after integration and rearrangement, one obtains at t
<br />;tp
<br />
<br />(r - K )
<br />L 0 t
<br />P - as - 00 p
<br />
<br />, , , , ' , , .(48)
<br />
<br />and from Eq. 44, at t; tp
<br />
<br />~KS(h - ~o)
<br />L -
<br />P (r - Kg)
<br />
<br />, , ' , , . .(49)
<br />
<br />Equating Eqs, 48 and 49 yields
<br />
<br />t =
<br />p
<br />
<br />(8 - 8 )(h - ~ )
<br />K s 0 0
<br />~ s (r - K )(r - K )
<br />o s
<br />
<br />, , , , ,(50)
<br />
<br />The values of Lp and tp can be determined from Eqs,
<br />49 and 50 respectively, Like Eq, 47, the fl value must
<br />be evaluated or given first in order for Eqs, 49 and 50
<br />to be useful,
<br />
<br />EvalWltion of the fl value, As defined in Eq, 42,
<br />the fl value depends on the K(z} and ",(z}
<br />distributions, which are, of course, unknown. How-
<br />ever, if one knows the K-"'relationship for the soil,
<br />the fl value can be readily evaluated from Eq, 42 and
<br />the K. '" relation, It is understood that the K.w
<br />relation is hysteretic, but will be assumed unique
<br />herein as long as the wetting process is only con.
<br />sidered, Several forms of the K.",relation have been
<br />assumed by different investigators in an attempt to
<br />best fit the field or experimental data, Because the
<br />soil after ponding (for h>O) has both the saturated
<br />and unsaturated zones, all the forms of the K.",
<br />relation must be modified to include K; K, for ",;;.0.
<br />
<br />The simplest of the available K-'" relations is a
<br />linear relationship between K and "', proposed by
<br />Richards (1931),
<br />
<br />K .. a;' + Ks
<br />
<br />for ;. :OS 0
<br />
<br />(5la)
<br />
<br />K = K
<br />s
<br />
<br />for ;. ~ 0
<br />
<br />(Sib)
<br />
<br />where "an is simply a constant.
<br />
<br />Raats (1971) and Braester (1973) among others
<br />used an exponential form in K and "'.
<br />
<br />K" K ealjJ
<br />s
<br />
<br />for 1/1 :S 0
<br />
<br />(52a)
<br />(52b)
<br />
<br />K ,., K
<br />s
<br />
<br />for I/J 2:: 0
<br />
<br />It should be noted that Eqs, 51 and 52 have only one
<br />parameter, a, needed to be determined from the K-",
<br />data,
<br />
<br />Several researchers proposed multi-parameter
<br />models which should fit the K-'" data better than the
<br />one-parameter model. For example, Gardner's (1958)
<br />original three.parameter model for "'<0 may be
<br />extended to include "'; 0 so that
<br />
<br />K -
<br />
<br />a
<br />
<br />(53a)
<br />
<br />for I/J :S 0
<br />
<br />(_~)n + b
<br />
<br />K - K
<br />s
<br />
<br />for tjI ~ 0
<br />
<br />(53b)
<br />
<br />where a, b, and n are all constants. However, the
<br />necessity of having K ; K, at "'; 0 results in the loss
<br />of freedom for Eq, 53a to choose the a value other
<br />than bI(,. Conversely, if K ; K, at '" ; 0 is merely
<br />considered as a data point in the K-", relationship, the
<br />values of a, b, and n in Eq, 53a can be determined by
<br />using a least squares optimization procedure. Mter
<br />that, one should take the value of K, equal to a/b, In
<br />view of the difficulty in integrating Eq, 53a with
<br />respe<:t to '" upon substitution of Eq, 53a into Eq, 42
<br />for the evaluation of the fl value, the following
<br />compatible form similar to Eq, 53a is proposed
<br />herein:
<br />
<br />K-
<br />
<br />a
<br />(_". + b)n
<br />
<br />(54a)
<br />
<br />for 1/J:5 0 . . . . .
<br />
<br />K - K
<br />s
<br />
<br />for ~~ 0 . , , , , (54b)
<br />
<br />By the same token, the values of a, b, and n in Eq,
<br />54a can be estimated by using a least squares
<br />optimization technique, and hence the value of K,
<br />may be taken to be equal to a/b".
<br />
<br />The values of a, b, and n determined from Eqs.
<br />53 and 54 are not all dimensionless, In order to make
<br />these dimensionally different parameters consistent
<br />(i.e" dimensionless), some investigators proposed
<br />dimensionless forms of Eq, 53 with the same number
<br />of parameters used in Eq, 53, Of them, Wei (1971)
<br />has used the most general one that is nevertheless not
<br />adopted in this study due to the same reason as given
<br />for Eq, 53,
<br />
<br />Substituting Eqs, 51, 52, and 54 into Eq, 42
<br />yields
<br />
<br />34
<br />
|