Incorporating auxiliary geophysical data into ground-water flow parameter estimation

The incorporation of auxiliary data into ground-water flow parameter estimation is a challenging task which can ultimately result in a better site characterization. In this study a maximum likelihood estimation procedure has been applied to the joint identification of the parameters of the aquifer transmissivity random field, and the parameters of the linear regression between the logarithm of transmissivity and the logarithm of the electrical transverse formation factor (TF), determined from surface geoelectrical methods (Vertical Electrical Sounding or V.E.S.). This approach is basically a co-kriging technique applied to the transmissivity and transverse formation factor random fields, but it avoids the independent estimation of the cross-covariances and the secondary variable covariance. The procedure needs some direct well data for transmissivity and a (usually larger) number of V.E.S. measurements which have to be in part at a distance from the well locations in order to provide useful information. The algorithm determines the characteristics of the local (site dependent) transmissivity-transverse formation factor relationship and utilizes this auxiliary information for a geostatistical transmissivity field estimation. The methodology is tested on a real field scenario: a fractured aquifer impacted by landfill leachate contamination. The use of the formation factor in place of the raw resistivity of the subsoil layers accounts for possible effects of clay and contaminant concentration on pore-water resistivity. The information provided by the V.E.S. can add, to some extent, to the understanding of the aquifer characteristics and vulnerability However, the specificity of each site has to be fully understood for an effective application of the present procedure. It seems unlikely that geoelectric data can differentiate between transmissivity values differing by less than two or three orders of magnitude.