HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Figures Only Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Basile, A. Articles by Randazzo, L. Search for Related Content PubMed Articles by Basile, A. Articles by Randazzo, L. GeoRef GeoRef Citation Agricola Articles by Basile, A. Articles by Randazzo, L. Related Collections Hydraulic Conductivity/Relative Permeability Hysteresis Water Retention/Capillary Pressure

ORIGINAL RESEARCH

Scaling Approach to Deduce Field Unsaturated Hydraulic Properties and Behavior from Laboratory Measurements on Small Cores

^a Institute for Mediterranean Agricultural and Forestry systems (ISAFOM), National Research Council (CNR), Ercolano (NA), Italy
^b Department for Agro-Forestry Systems Management (DITEC), University of Basilicata, Potenza, Italy
^c Department of Agricultural Engineering and Agronomy (DIAAT), University of Naples "Federico II", Portici (NA), Italy
^* Corresponding author (acoppola{at}unibas.it)

Received 31 October 2005.

	ABSTRACT

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Hydraulic properties should be determined at the scale of the process modeled. The methods to hydraulically characterize a soil in situ remain extremely difficult to implement, requiring measurements of water content and pressure head with adequate time–depth resolution. We recently proposed a method that reduced the number of field measurements required for complete field hydraulic characterization. In this paper, we extend the previous method to cases where the only available information consists of laboratory hydraulic properties, along with a measurement of the maximum water content in the field, whose value depends on the way the wetting of the porous medium is performed. Specifically, the field hydraulic parameters were estimated with a scaling procedure accounting for the ratio between the total porosity estimated by the maximum water content in the laboratory and the partial porosity effectively involved in the field, as estimated by the field-measured maximum water content. The scaling method was evaluated with data from four soils (three volcanic sandy loam soils and one silty clay loam soil). Soil hydraulic properties were measured both in situ by a field internal drainage test and in the laboratory in an evaporation experiment (Wind's method). Scaling-based hydraulic properties were also compared with those estimated by applying a simplified method based on the unit-gradient water flow assumption where only water content measurements performed during the internal drainage test were of concern. The hydraulic properties estimated with the scaling method were compared with the measured ones and with those from the unit-gradient method in terms of the Relative Mean Error (RME) and Relative Root Square Mean Error (RRMSE). The scaling method proved to be especially effective when applied to the three sandy loam soils, where scaled retention and hydraulic conductivity curves to a large extent reproduced those measured in the field. For the silty clay loam soil, appropriate results were observed only for the water retention curve, while poorer scaled hydraulic conductivity was obtained.

Key Words: REV • representative elementary volume

	INTRODUCTION

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
EFFICIENT MANAGEMENT of soil and water necessitates characterization of soil hydraulic behavior across large areas. This requires that soil water–pressure head, (h), and hydraulic conductivity–water content, k(), relationships be specified at the scale of concern.

Laboratory measurements are more straightforward than those in the field, but their validity is essentially related to the sample sizes, which have to appropriately represent the heterogeneity of the medium being studied. A large body of literature discusses the validity of laboratory hydraulic characterization in inferring field hydrological behavior (van Genuchten et al., 1999).

The most rigorous field method to determine k() and (h) is the instantaneous profile method (Watson, 1966). Indeed, at the field scale these functions are generally obtained by monitoring a redistribution process along a profile of limited dimensions in space (pedon scale) through direct measurements of both soil water content and soil water pressure head profiles.

This technique is cumbersome and time-consuming. Therefore, two main simplifications have been largely applied: (i) considering approximations about the flow regimes during redistribution or (ii) drawing field hydraulic variables partly from lab measurements.

The first condition has produced several approximate methods, mainly based on the assumption of unit hydraulic head gradients, that are useful to characterize hydraulic relationships by measuring only water contents during redistribution (Libardi et al., 1980; Sisson et al., 1980; Chong et al., 1981). The method proposed by Libardi et al. (1980) is additionally based on the assumptions of a linear water content–storage relationship and an exponential k() relationship. Following the analysis provided by Flühler et al. (1976), Sisson and van Genuchten (1991) assumed the adequacy of unit-gradient water flow models given that the instantaneous profile method in its classical form would suffer from the uncertainties arising from noises, especially those observed in pressure head measurements. Actually, the malfunction of even one tensiometer introduces significant errors into the calculation of hydraulic properties. Their study showed that the data of the instantaneous profile method applications could be used in a Cauchy problem to find hydraulic parameters. However, though attractive, such simplified methods should be applied with caution because the founding assumptions were only rarely validated on critical field measurements. As stated by Kutilek and Nielsen (1994), "the accuracy of field-measured estimates of hydraulic properties will increase as the number of simplifying assumptions decreases."

In the second approach, the field method requires the direct measurement of either soil water content or soil water pressure head profiles and the indirect determination of the other variables through the separately lab-determined water retention curves. For example, Shouse et al. (1992a, 1992b) estimated the soil water contents in fields by coupling pressure head measured in the field and laboratory (h). Ahuja and Williams (1991) and Rockhold et al. (1996) estimated (h) in the field from measurements on soil cores. In the latter methods complications arise because in situ determined hydraulic properties often disagree with those determined on undisturbed samples collected at the same site, in a way that field hydraulic properties cannot be straightforwardly drawn from laboratory ones. As shown in Fig. 10 of the paper by Pachepsky et al. (2004), if one considers the UNSODA database (Leij et al., 1996), at the same pressure heads, field water contents are almost always markedly lower than laboratory values. Actually, what is generally observed is that the field retention curves are significantly lower than laboratory ones (Basile et al., 2003). Basile et al. (2003) also argued that the laboratory desorption curves represent the maximum values that may occur in the field.

Laboratory hydraulic conductivities at h = 0, k₀, are generally observed to be different from those measured in the field. It has been amply shown (Eching et al., 1994; Rockhold et al., 1996; Basile, 2004, unpublished data, among others) that the laboratory determined k₀ can be one or more orders of magnitude greater than those measured in the field. Consequently, erroneous estimates of field absolute hydraulic conductivities can arise when using the lab k₀ to match field relative hydraulic conductivities.

At present, there are still great difficulties in providing quantitative explanations and formalization of differences in laboratory vs. field hydraulic properties. Recent efforts (Basile et al., 2003) aimed to propose criteria to systematically manage the observed dissimilarities between field and laboratory behavior in a way that laboratory information could be converted to the field.

The main practical purpose for which the formulation was developed was to reduce the number of field measurements required for a complete field hydraulic characterization.

The objective of this study was to formulate a scaling approach that allows the method proposed by Basile et al. (2003) to be extended to cases when only laboratory hydraulic properties are available. This objective is crucial since field-scale water flow predictions are frequently made with only laboratory or small-scale data available.

The appropriateness of the proposed approach was evaluated by comparing the scaling-based predictions of hydraulic properties with those from the classical analysis of the instantaneous profile method (Watson, 1966) and those from the unit-gradient simplification proposed by Sisson and van Genuchten (1991). A numerical analysis was also provided to establish the applicability conditions of the method depending on the ratio between field/lab ₀ values, as well as on the ratio between field or lab values and the laboratory n parameter.

	THEORY

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
In their study, Basile et al. (2003) determined the hydraulic properties of the soils at the pedon and laboratory scales. Soil water retention was described by the unimodal (h) relationship proposed by van Genuchten (1980) and expressed here in terms of the effective saturation, S_e, as follows:

[1]

with S_e = ( – _r)/(₀ – _r), _r and ₀ being the residual water content and the water content at h = 0, respectively, and in which (cm^–1), n, and m are curve-fitting parameters.

Mualem's expression was used to calculate relative hydraulic conductivity, k_r (Mualem, 1976). Assuming m = 1 – 1/n, van Genuchten (1980) obtained a closed-form analytical solution to predict k_r at a specified volumetric water content:

[2]

in which k₀ is the hydraulic conductivity measured at ₀, and is a parameter which accounts for the dependence of the tortuosity on the water content. In this study was fixed at a value of 0.5.

For all the horizons a discrepancy was found between the water content values at h = 0 obtained with the two methods, with a tendency for the water retention values to be higher for the laboratory curves. The two-parameter vectors differed only in parameter ₀ and . The latter was higher for the laboratory, as expected given that the larger pores, nearly all filled with water in the laboratory tests, empty at relatively lower pressure head values, h. Once the larger pores have been emptied, at a value of h approximately equal to 1/, the fraction of pores that empties in response to a given variation in the pressure head h tends to become the same both for the field and laboratory samples, as indicated by comparable values of parameter n that were observed for the two curves.

The authors ascribed such behavior to the nearly complete saturation (i.e., a higher fraction of air removed during wetting) that could be achieved in the laboratory and interpreted such differences as a result of different hysteretic paths being observed, arising in turn from different wetting procedures. This hypothesis was reinforced by observations of nonhysteretic k() curves when field and laboratory measurements were compared, which, in the van Genuchten–Mualem framework, may only be expected when similar n values are also obtained (see Eq. [2]). In order that hydraulic conductivity curves, k(), almost overlapping, may be derived from k_r() curves arising from different lab–field retention curves, the matching factors for relative conductivity in the field and in the laboratory have to be different. They proposed an expression (Eq. [7] in Basile et al., 2003) to derive the field k₀ from the lab one, which for applications can be put in the following form, accounting for the differences in the whole range of the field hydraulic conductivity curve:

[3]

where ^Fk_r is derived from laboratory n and by parameters ₀ and determined in the field in the early phases of the field test. In practice, Eq. [3] permits scaling of the ^Lk₀ values by the ratio between the integrals of the k_r functions. In Eq. [3] and in the following, labels L and F are used to signify laboratory and field, respectively.

In this study the field hydraulic properties are deduced from the laboratory ones and from parameter ₀ determined in a field test. If the theoretical framework holds, a simple scaling procedure of the laboratory retention curve can be applied to derive the field retention curve. With the (h) parameters, the k() field curve is estimated through Mualem's model by using a k₀ matching factor, lower than the lab one, obtained as described in Eq. [3]. A schematic view of the proposed scaling procedure is given in Fig. 1 , the formalization of which is given in Eq. [4]:

[4]

where h_c is a critical pressure head taken on the lab water retention curve and corresponding to the ₀ on the field curve, while S_e is the effective saturation.

View larger version (9K): [in this window] [in a new window]   Fig.1. Schematic of the scaling procedure: (A) laboratory (dashed line) and field (solid line) measured water retention curves, (B) effective saturation (Se) of the measured laboratory curve; (C) laboratory-measured (dashed line) and field-scaled (bold solid line) probability density functions (PDFs), (D) laboratory- (dashed line) and field-measured (solid line) water retention curves, and field-scaled retention curve (bold solid line). The critical pressure head, hc, is the value on the laboratory curves corresponding to F0.

1. The porous media being compared are the same, a condition in principle applying each time the lab sample is actually a representative elementary volume (REV) for the pedon scale.
   
2. Due to the wetting procedure, in the field the porous system is only partially filled with water at h = 0.
   
3. The part of the porous system not involved in the field wetting is in the range of larger pores.
   

Figure 1A shows a typical case of lab–field differences in the retention curve. The S_e(h) curve in Fig. 1B corresponds to the (h) curve in the laboratory (Fig. 1A, dashed line). We recall that the actual saturation, S_e, is considered a cumulative distribution function of pore size with a density function f(h) that may be expressed by the following equation (Durner, 1994):

[5]

The f(h) curve corresponding to Fig. 1B is depicted in Fig. 1C (dashed line). Assumptions 1, 2, and 3 lead one to consider that the fraction of the total porous system (as detected in the laboratory) actually involved in the field corresponds to the shaded side of the lab pore size distribution in Fig. 1C. Through Eq. [4], the shaded part of the lab distribution transforms to the estimated field distribution (solid line in Fig. 1C) to give a ^FS_e(h) = 1 at a water content ^F₀ and at h = 0. In practice, this corresponds to shifting the part of the lab curve from ^LS_e(^F₀) to ^LS_e = 0 to the left along the h axis of h – h_c and upward along the S_e axis of ^LS_e(h)/^LS_e(h_c). Finally, the scaled S_e(h) converts to the field (lab-scaled) (h) curve in Fig. 1D (dashed line).

This scaling formulation is something like the approach of the most commonly used hysteresis scaling procedures (Kool and Parker, 1987; Luckner et al., 1989).

Once scaled ^FS_e(h) are estimated according to Eq. [4], the relevant value of can be determined by an optimization procedure.

	MATERIALS AND METHODS

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
The study concerned four soil profiles whose main soil physical and textural characteristics are reported in Table 1. The soils are located in alluvial plains (Piano Campano and Piana del Volturno, Campania region, Southern Italy). Three of these (site 2, site 3 and site 15) are sandy loam soils including volcanic sediments and ash fall. These volcanic soils are very fertile and intensively cultivated. The VIT soil, located in the Volturno plain, is a silty clay loam soil with vertic characteristics.

View larger version (69K): [in this window] [in a new window]   Fig. 2. Schematic of instantaneous profile method for the (a) infiltration and (b) redistribution stages.

The water retention curve for each horizon was obtained by coupling the measurements of and h performed at the same depth and time. The unsaturated hydraulic conductivity curve was obtained by the algorithm suggested by Watson (1966), hereafter summarized following Kutilek and Nielsen (1994).

Water content (z, t) measurements allow calculation of the water stored W_t at time t in the profile between the soil surface and depth z:

[6]

The average flux density during a time interval t = t₂ – t₁ is (t) = –W/t, W being the water loss from the soil layer (0, z). Substituting the average flux density in the finite difference form of the Darcy equation, = –k()H/z, yields

[7]

where the gradients of the total potential head, H (=h + z), are calculated from h(z,t) measurements, and is related to the mean of h in H.

Sisson et al. (1980) estimated hydraulic conductivity without the need for calculating hydraulic gradients and water flux densities. Assuming unit-gradient water flow conditions and a relatively homogeneous soil profile, the Richards equation reduces to

[8]

For a initial condition of (z,0) = _in(z), the solution of the Eq. [8], known as a Cauchy problem (Lax, 1972), is

[9]

that is, the ratio of the depth to time at which water contents are measured. When the profile is saturated initially, _in equals the water content at saturation.

Sisson and van Genuchten (1991) estimated the parameters in Eq. [2] by equating the experimental values of the derivative of the k() (Eq. [9]) to the derivative of the Eq. [2]:

[10]

with A = 1 – S_e^1/m, in a parameter optimization procedure. According to the approach proposed by Sisson and van Genuchten (1991), the parameter in Eq. [1] (and thus the retention curve) is also estimated by using an experimental (h) point that, in our study, was taken to be h = –60 cm.

Laboratory Method
In the laboratory, the (h) and k() relationships were determined by means of an evaporation experiment (Wind, 1968) according to the calculus procedure proposed by Tamari et al. (1993). Soil samples ( = 86 mm, h = 150 mm) were collected in duplicate at the end of the field experiment in the middle of each horizon (Fig. 2b). In the laboratory, the samples, after preliminary air drying, were saturated from the bottom and the saturated hydraulic conductivity measured by a permeameter (Reynolds and Elrick, 2002). Then, the sample drying process was conducted by starting from hydrostatic equilibrium conditions, setting a zero flux at the bottom surface, and at the same time, at appropriately preset time intervals, monitoring the weight of the whole sample and the pressure heads at two different depths by means of tensiometers. The equipment and data acquisition were completely automated. Measurements were interrupted when air entered the upper tensiometer. In relation to the texture of the different samples, each test lasted between 3 and 5 d.

An iterative procedure was applied for estimating the water retention curve. We started with a "guess" water retention curve, with initial estimates of parameters in Eq. [1], to convert pressure head data at any one time into estimated water contents at the tensiometer depths., The average water content of the soil sample was then calculated and compared with the average measured water content obtained from the sample weight. If significant differences between calculated and measured average water content were observed, a new parameter set was estimated. The water retention curve was optimized using an iterative procedure, by minimizing the squared differences between calculated and measured average water content. After convergence of the algorithm, water contents at the tensiometer depths were calculated from the measured pressure heads and the estimated water retention curve. The distributions of the measured pressure heads, h(z, t), and of the estimated water contents, (z, t), were used to calculate the hydraulic conductivity according to the instantaneous profile method previously described.

Moreover, three points—at 3, 6, and 12 bars for Site 2 and Site 15—and four points—at 2, 4, 12, and 15 bars for Site VIT—of the water retention curve were determined by means of the pressure plate apparatus.

	RESULTS AND DISCUSSION

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
In this section, details will be given on (i) comparisons between field and laboratory measurements of hydraulic properties, to evaluate the scaling framework; (ii) comparison between field measured hydraulic properties and those deduced by applying the scaling procedure. The goodness of the scaling procedure will be evaluated by comparing relevant estimates with those derived by the unit-gradient method; and (iii) numerical analysis of the applicability conditions of the scaling procedure in terms of van Genuchten model parameters, to quantify the extent of the ratios between field and laboratory measured hydraulic parameters for the scaling procedure to be expected to be valid.

Hereafter, labels WR and HC will be used to signify water retention and hydraulic conductivity, respectively, with the superscripts F for field curves and L for laboratory curves. Moreover, subscripts m, s, and ug will be used to identify measured, scaled, and unit-gradient curves and parameters, respectively.

Comparing Measured Field and Laboratory Hydraulic Properties
The water retention values from field (circles) and laboratory (triangles) experiments for the Ap horizons of the investigated sites are plotted in Fig. 3 . They were interpolated by applying the Eq. [1] with m = 1 – 1/n. Relevant parameters are shown in Table 2, along with parameters pertaining to the other horizons not shown herein. The graph shows that for all horizons there is a discrepancy between the water content values at h = 0 obtained with the two methods, with water retention values being higher for the laboratory curves. For lower pressure head values the curves tend to converge progressively. The gradual convergence of the laboratory and field retention curves past the inflection point is always reflected in similar values of parameter n. Analogous considerations can be made for the other horizons on the basis of the value of the parameters shown in Table 2.

View larger version (10K): [in this window] [in a new window]   Fig. 3. Soil water retention curves of the Ap horizons. Triangles and circles refer to laboratory and field experimental data, respectively; lines are for relevant curve-fittings by Eq. [1].

View this table: [in this window] [in a new window]   Table 2. Parameters of the water retention curve (Eq. [1]). The field parameters were obtained by fixing n at the laboratory value, 0 at the measured value in the field and setting m = 1 – 1/n. Scaled and unit gradient parameters were obtained according to the procedures described in the text.

View larger version (10K): [in this window] [in a new window]   Fig. 4. Hydraulic conductivity curves of the Ap horizons. Triangles and circles refer to laboratory and field experimental data, respectively; dashed lines refer to the application of Mualem's model (Eq. [2]) to laboratory water retention data.

Comparison of Measured, Scaled and Unit-Gradient Hydraulic Properties
Figure 5 shows the results of applying the scaling procedure proposed (Eq. [4]) to derive the field retention curves from the lab ones. In the graph the scaled field retention curves (^FWR_s, solid lines) are compared with the field-measured curves (^FWR_m, circles). The relevant parameters are given in Table 2. Note the satisfactory, sometime very precise, overlapping of the reconstructed curves on the retention curves actually measured in the field, even for site VIT.

View larger version (14K): [in this window] [in a new window]   Fig. 5. Field soil water retention curves of the Ap horizons. Solid lines and circles refer to scaled and measured field water retention curve, respectively. Dashed lines refer to the application of the unit-gradient method.

View larger version (11K): [in this window] [in a new window]   Fig. 6. Field hydraulic conductivity curves of the Ap horizons. Circles refer to field measured data. Solid (FHCs(1)) and dashed (FHCs(2)) lines refer to the scaling method with kr() determined from measured and scaled water retention curves, respectively. Dotted lines refer to the application of the unit-gradient method.

These results seem to confirm that the hysteresis scheme hypothesized holds for explaining differences between lab and field measurements, as only in this case does the schematization leading to Eq. [3] and [4] hold. This also implies that the laboratory samples are actually REVs for the corresponding horizons at the pedon scale. In fact, the main condition for which the scaling procedure we proposed held is that we are dealing with the same porous medium.

In the same graphs (Fig. 5 and 6), the unit-gradient water retention and hydraulic conductivity curves are also shown (^FWR_ug in Fig. 5 and ^FHC_ug in Fig. 6). The relevant parameters are given in Table 2.

Although the Sisson et al. (1980) unit-gradient approach holds strictly only for uniform profiles (Sites 3 and VIT), it was applied (less rigorously) also to situations where the profile exhibited some weak layering (Sites 2 and 15).

What should be noted is that the scaling procedure always yielded results comparable to those from the unit-gradient approach, although in practice calling only for laboratory measurements. For Site 15, the ^FWR_ug behavior at lower h has to be attributed to the lack of information in the dry range during the field drainage experiment.

Difference-Based Statistics
By reasoning from analogy, similar results can be drawn by observing the difference-based statistics provided in Table 3, both for the water retention and hydraulic conductivity, for all the horizons investigated.

[11]

[12]

where N is the number of data points, E_i and M_i are the estimated and measured water contents or conductivities, respectively, and is the average value of the measured data. Depending on the goal of the comparison and on the curves to be compared, the values and meanings of the estimated (E_i) and the measured (M_i) water contents will change according to the scheme reported in Table 3. The E_i values are calculated from the relevant parameters at the same values of pressure head and water content of the M_i, for the water retention and hydraulic conductivity curve, respectively.

Relative mean error percentages reflect the average deviation between estimated and measured hydraulic data and were calculated to give information on the relative size of systematic errors as their positive (negative) values indicate the percentage of overestimation (underestimation) on the average.

The RRMSE is used to describe accuracy encompassing both random and systematic errors. It can be especially useful to provide information about the relative size of the differences between estimated and measured water retentions and hydraulic conductivities by normalizing the differences by the mean of measurements.

Every RME index of the ^LWR_m/^FWR_m comparison is positive, confirming for all the horizons investigated that the laboratory soil water retention curves are positioned above the field ones according to the framework reported. The other comparisons by means of the RME index for the water retentions do not show any trend. The scaled method for the hydraulic conductivity, ^FHC_s seems to be site-dependent, overestimating at Sites 2 and 3 and underestimating at Site 15. Instead, the unit-gradient method, ^FHC_ug, underestimates for all soils except for the Ap horizon of Site 2.

Figure 7 illustrates graphs of the RRMSE for the ^LWR_m–^FWR_m, ^FWR_s–^FWR_m, ^FWR_ug–^FWR_m for water retention series, and for the ^LHC_m–^FHC_m, ^FHC_s(2)–^FHC_m, ^FHC_ug–^FHC_m for hydraulic conductivity series.

View larger version (20K): [in this window] [in a new window]   Fig. 7. Relative Root Square Mean Error (RRMSE) indexes for water retention and hydraulic conductivity.

As for the F_s–F_m series, what should first be noted is that the scaling procedure provides estimated water retention curves differing by <10% from the measured ones (slightly higher for horizons Ap and 2Ab of Site 15) and frequently improved when compared to the unit-gradient method (see F_ug–F_m indexes) (Site 2 A, Site 2 Bw, Site 3 Bw, Site 15 A, Site 15 2Ab). Concerning the corresponding hydraulic conductivity estimations, differences are frequently very low, with values higher than 10% only for Site 3. Also considering the hydraulic conductivity predictions, the F_s series proved to be in a number of cases better than the F_ug series (Site 2 Ap, Site 2 A, Site 15 A, Site 15 2Ab). However, one should note that enhanced water retention predictions do not directly mean better predictions of hydraulic conductivities.

Overall, the scaling procedure provides results comparable to those from the unit-gradient method, with the difference that the former just requires measurement in the field of the water content at h = 0, while the latter calls for field water content profile measurements over time and at least one (h) experimental point for estimating the relevant (see Table 4 for the comparison of the main requirements of each applied method).

Both in the tables and in the graphs, difference-based statistics are also shown for the Ap horizon of Site VIT. It is worth noting that Site VIT is a structured silty-clay loam soil, with a larger REV for hydraulic properties (at least 10³ cm³) than that required for well-sorted soils, especially for hydraulic conductivity (Kutilek and Nielsen, 1994). For such a soil, the samples to be analyzed in the laboratory, given their sizes, were not expected to represent the REV of the soil horizon at the pedon scale. This is the classic case where the schematization we proposed does not apply because the responses observed in the lab and in the field refer to two different porous media, each with its own hydraulic properties. From the RRMSE index graphs, larger discrepancies between estimated and measured field hydraulic conductivities emerged (see also Fig. 6),although with acceptable water retention curves (see Fig. 5).

It could also be argued that acceptable water retentions and poor hydraulic conductivity estimations are largely to be attributed to the sample size we used, which would suffice as a representative volume for water retention but would be inadequate as a REV for hydraulic conductivity.

However, a better estimation of hydraulic conductivity was obtained by the unit-gradient method, since it is in any case based on water content measured directly in the field.

Numerical Analysis on the Applicability Conditions of the Scaling Procedure
Figures 8a , 8b, and 8c convey information on the applicability conditions of the proposed scaling procedure in terms of van Genuchten parameters. The graph shows trends of the ^F/^L ratio as related to the ^F₀/^L₀ ratio for different values of shape parameters ^L and n and for a starting value of ^L₀ = 0.5, although it should be pointed out that the plots are independent of the starting value of ^L₀. Each line of any plot provides the value of ^F/^L for a given ^F₀/^L₀ ratio and n value so that the scaling procedure holds. In a sense, it might be said that when the ^F/^L ratio adjusts according to the plots as those in Fig. 5, the selected lab sample size may well represent the REV for horizon at the pedon scale.

View larger version (35K): [in this window] [in a new window]   Fig. 8. Trends of F/L ratio as related to the F0/L0 ratio for different values of shape parameters L and n and for a starting value of L0 = 0.5 (Plots a, b, c). Plot d is the same as Plot b with Fn = Ln x 1.1.

Additionally, to give an example of the errors one can make by applying the proposed scaling procedure when the hysteresis framework does not perfectly hold and, consequently, the ^Ln parameter does not represent the ^Fn, Fig. 8d describes analogous relationships for a case of ^Fn = ^Ln x 1.1. In other words, Fig. 8d provides the ^F/^L ratios necessary for the laboratory and the scaled field retention curves to overlap, for the case of ^Fn values 10% higher than ^Ln. Differences between ^Fn–^Ln might occur either due to the inability of the lab sample size to encompass the REV for the horizon or due to measurement errors, or might even be the result of the ^Ln being estimated by a pedotransfer function. Of course, this circumstance is only of speculative concern, as ^Fn should be known in order for such a graph to be of practical interest, a condition that the proposed scaling procedure aims above all to overcome. As an example, while applying the scaling procedure we proposed, for ^Ln = 1.3 and for a ^F₀/^L₀ ratio equal to 0.8, one would select an ^F/^L ratio of 0.37 from Fig. 8b. By contrast, to predict the correct field retention curve for a real ^Fn 10% higher than the ^Ln, one should use a ^F/^L ratio of 0.20 from Fig. 8d, thus with an error in estimating ^F of about 100%. From the same graph in Fig. 7 one might easily note that the higher the ^Ln value, the smaller the error in estimating ^F, thus suggesting that the robustness of the proposed scaling procedure increases with the ^Ln values. Analogous considerations can be made for k₀ predictions (results not shown).

	CONCLUSIONS

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
The procedure we proposed, derived in the framework of a hysteretic interpretation of the differences between lab and field hydraulic properties, allows reproduction of both the field-measured retention and hydraulic conductivity with adequate approximation for many purposes. Despite the expected differences between the measured and scaling-based field hydraulic properties, the scaling procedure provided hydraulic property predictions whose accuracy was comparable to that obtained by the unit-gradient approximation, although there existed some cases where either the (h) or the k() curve was better described by the unit-gradient method. However, it should be noted that the unit-gradient method requires measurements over time of the soil water content along the profile while our scaling method needs just one soil water content measurement in the field at h = 0 and laboratory characterization. Moreover, at least in principle, there are no restrictions in the applicability of the proposed scaling method (i.e., a priori assumptions on the gradient or the type of process). The only constraint is that the soil sample has to be large enough to be the REV for the horizon at the pedon scale.

Within the scaling framework proposed, it is important to match the lab and field experimental hydraulic conductivity by different matching points, depending on the initial and boundary conditions applied. The field k₀, which is lower than its lab counterpart, can be estimated from the ^Fk₀ formula (Eq. [3]) derived in the theoretical framework. When this value is used to match Mualem's relative hydraulic conductivity estimated from the retention curve scaled according to the proposed procedure, field-estimated hydraulic conductivities roughly reproduce the measured ones.

The concept that different initial and boundary conditions result in different hysteresis curves being measured may also be extended given that, in practice, different measurement techniques may result in different initial and boundary conditions being applied, such that the hydraulic behavior we actually measure depends on the technique used, rather than being an intrinsic attribute of soils.

The proposed procedure provides a systematic technique for scaling large laboratory data sets to the field, either deriving from direct measurements or estimated by pedotransfer functions, which too frequently are incorrectly used for field applications.

	REFERENCES

TOP EXECUTIVE SUMMARY ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 

• Ahuja, L.R., and R.D. Williams. 1991. Scaling water characteristics and hydraulic conductivity based on Gregson-Hector-McGrown approach. Soil Sci. Soc. Am. J. 55:308–319.[Abstract/Free Full Text]
• Basile, A., G. Ciollaro, and A. Coppola. 2003. Hysteresis in soil-water characteristics as a key to interpreting comparisons of laboratory and field measured hydraulic properties. Water Resour. Res. 39(12):1355. doi:10.1029/2003WR002432.[CrossRef]
• Chong, S.K., R.E. Green, and L.R. Ahuja. 1981. Simple in situ determination of hydraulic conductivity by power function descriptions of drainage. Water Resour. Res. 17:1109–1114.
• Durner, W. 1994. Hydraulic conductivity estimation for soils with heterogeneous pore structure. Water Resour. Res. 30(2):211–223.[CrossRef]
• Eching, S.O., J.W. Hopmans, and W.W. Wallender. 1994. Estimation of in situ unsaturated soil hydraulic functions from scaled cumulative drainage data. Water Resour. Res. 30:2387–2394.[CrossRef]
• Flühler, H., M.S. Ardakani, and L.H. Stolzy. 1976. Error propagation in determining hydraulic conductivities from successive water content and pressure head profiles. Soil Sci. Soc. Am. J. 40:830–836.[Abstract/Free Full Text]
• Jørgensen, S.E., L. Kamp-Nielsen, T. Christensen, J. Windolf-Nielsen, and B. Westergaard. 1986. Validation of a prognosis based upon a eutrophication model. Ecol. Modell. 35:165–182.
• Kool, J.B., and J.C. Parker. 1987. Development and evaluation of closed-form expressions for hysteretic soil hydraulic properties. Water Resour. Res. 23:105–114.
• Kutilek, M., and D.R. Nielsen. 1994. Soil hydrology. Catena Verlag, Geoscience Publ., Reiskirchen, Germany.
• Lax, P.D. 1972. The formation and decay of shock waves. Am. Math. Mon. 79:227–241.[Medline]
• Leij, F.J., W.J. Alves, and M.Th. van Genuchten. 1996. The UNSODA Unsaturated Soil Hydraulic Database: User's manual version 1.0. Tech. EPA/600/R-96/095. Natl. Risk Management Lab., Cincinnati, OH.
• Libardi, P.L., K. Reichardt, D.R. Nielsen, and J.W. Biggar. 1980. Simple field methods for estimating soil hydraulic conductivity. Soil Sci. Soc. Am. J. 44:3–7.
• Luckner, L., M.Th. van Genuchten, and D.R. Nielsen. 1989. A consistent set of parametric models for the two-phase flow of immiscible fluids in the subsurface. Water Resour. Res. 25:2187–2193.
• Mualem, Y. 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12:513–522.[CrossRef]
• Pachepsky, Y.A., K.R.J. Smettem, J. Vanderborght, M. Herbst, H. Vereecken, and J.H.M. Wosten. 2004. Reality and fiction of models and data in soil hydrology. In R.A. Feddes et al. (ed.) Unsaturated-zone modeling. Kluwer Academic Publishers, Dordrecht, The Netherlands.
• Reynolds, W.D., and D.E. Elrick. 2002. The soil solution phase. Falling head soil core (tank) method. p. 809–812. In J.H. Dane and G.C. Topp (ed.) Methods of soil analysis. SSSA Book Ser. 5. SSSA, Madison, WI.
• Rockhold, M.L., R.E. Rossi, and R.G. Hills. 1996. Application of similar media scaling and conditional simulation for modelling water flow and tritium transport at the Las Cruces Trench Site. Water Resour. Res. 32:595–609.[CrossRef]
• Shouse, P.J., J.B. Sisson, G. de Rooij, J.A. Jobes, and M.Th. van Genuchten. 1992a. Application of fixed gradient methods for estimating soil hydraulic conductivity. p. 675–684. In M.Th. van Genuchten et al. (ed.) Proceedings of the International Workshop on Indirect Methods for Estimating the Hydraulic Properties of Unsaturated Soils. University of California Press, Riverside.
• Shouse, P.J., J.B. Sisson, T.R. Ellsworth, and J.A. Jobes. 1992b. Estimating in situ unsaturated hydraulic properties of vertically heterogeneous soils. Soil Sci. Soc. Am. J. 56:1673–1679.[Abstract/Free Full Text]
• Sisson, J.B., A.H. Ferguson, and M.Th. van Genuchten. 1980. Simple method for predicting drainage from field plots. Soil Sci. Soc. Am. J. 44:1147–1152.[Abstract/Free Full Text]
• Sisson, J., and M.Th. van Genuchten. 1991. An improved analysis of gravity drainage experiments for estimating the unsaturated soil hydraulic functions. Water Resour. Res. 27:569–575.[CrossRef]
• Tamari, S., L. Bruckler, J. Halbertsma, and J. Chadoeuf. 1993. A simple method for determining soil hydraulic properties in the laboratory. Soil Sci. Soc. Am. J. 57:642–651.[Abstract/Free Full Text]
• van Genuchten, M.Th. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44:892–898.[Abstract/Free Full Text]
• van Genuchten, M.Th., F.J. Leij, and L. Wu (ed.) 1999. Characterization and measurement of the hydraulic properties of unsaturated porous media. Parts 1 & 2. Proceedings of the International Workshop on Characterization and Measurement of the Hydraulic Properties of Unsaturated Porous Media, Riverside, CA.
• Warrick, A.W. 1998. Spatial variability. p. 655–676. In Environmental soil physics. D. Hillel (ed.) Academic Press, New York.
• Watson, K.K. 1966. An instantaneous profile method for determining the hydraulic conductivity of unsaturated porous materials. Water Resour. Res. 2:709–715.
• Williams, R.D., L.R. Ahuja, and J.W. Naney. 1992. Comparison of methods to estimate soil water characteristics from soil texture and limited data. Soil Sci. 153:172–184.
• Wind, G.P. 1968. Capillary conductivity data estimated by a simple method. p. 181–191. In P.E. Rijtema and H. Wassink (ed.) Water in the unsaturated zone. Vol. 1. Int. Assoc. Sci. Hydrol., UNESCO, Paris.

This Article Abstract Figures Only Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Basile, A. Articles by Randazzo, L. Search for Related Content PubMed Articles by Basile, A. Articles by Randazzo, L. GeoRef GeoRef Citation Agricola Articles by Basile, A. Articles by Randazzo, L. Related Collections