Abstract. A linear estimator, cokriging, was applied to estimate hydraulic conductivity, using pressure head, solute concentration, and solute arrival time measurements in a hypothetical, heterogeneous vadose zone under steady state infiltrations at different degrees of saturation. Covariances and cross-covariances required by the estimator were determined by a first-order approximation in which sensitivity matrices were calculated using an adjoint state method, The effectiveness of the pressure, concentration, and arrival time measurements for the estimator were then evaluated using two statistical criteria, L • and L 2 norms, i.e., the average absolute error and the mean square error of the estimated conductivity field. Results of our analysis showed that pressure head measurements at steady state flow provided the best estimation of hydraulic conductivity among the three types of measurements. In addition, head measurements of flow near saturation were found more useful for estimating conductivity than those at low saturations. The arrival time measurements do not have any significant advantage over concentration. Factors such as variability, linearity, and ergodicity were discussed to explain advantage and limitation of each type of data set. Finally, to take advantage of different types of data set (e.g., head and concentration), a computationally efficient estimation approach was developed to combine them sequentially to estimate the hydraulic conductivity field. The conductivity field estimated by using this sequential approach proves to be better than all the previous estimates, using one type of data set alone.
IntroductionDuring the past decades, the cokriging technique has been applied extensively to many studies of subsurface hydrology. timates of the heterogeneous conductivity field. On the other hand, Yeh et al. [1995, 1996] proposed iterative cokriging techniques for nonlinear systems in which the requirement of unbiasedness and minimum variance were imposed during each iteration. By incorporating the nonlinear relationship between head and conductivity in groundwater flow systems, the estimated conductivity field revealed more detailed heterogeneity than using the linear model and more closely resembled the true field. This iterative approach was further extended to unsaturated flow by Zhang and Yeh [1997] to estimate parameters for unsaturated hydraulic conductivity in the vadose zone. Similarly, a quasi-linear geostatistical approach was presented by Kitanidis [1995] in an attempt to incorporate the nonlinear relationship between the parameter and secondary information of the subsurface flow system.In spite of the advance in parameter estimation techniques, a practical question remains regarding what kind of measurements (e.g., head, concentration, or other variables) are most useful for estimating the hydrological parameters. Several attempts to address this issue were carried out in the past. For example, Yeh and Zhang [1996] found that pressure head measurements can improve estimates of s...