Water table elevations are usually sampled in space using piezometric measurements that are unfortunately expensive to obtain and are thus scarce over space. Most of the time,
piezometric data are sparsely distributed over large areas, thus providing limited direct information about the level of the corresponding water table. As a consequence, there is a real need for approaches that are able at the same time to (1) provide spatial predictions at unsampled locations and (2) enable the user to account for all potentially available secondary information sources that are in some way related to water table elevations. In
this paper, a recently developed Bayesian data fusion (BDF) framework is applied to the problem of water table spatial mapping. After a brief presentation of the underlying
theory, specific assumptions are made and discussed to account for a digital elevation model and for the geometry of a corresponding river network. On the basis of a data set for the Dijle basin in the north part of Belgium, the suggested model is then implemented and results are compared to those of standard techniques such as ordinary kriging and cokriging. Respective accuracies and precisions of these estimators are finally evaluated using a ‘‘leave-one-out’’ cross-validation procedure. Although the BDF methodology was illustrated here for the integration of only two secondary information sources (namely, a digital elevation model and the geometry of a river network), the method can be applied for incorporating an arbitrary number of secondary information sources, thus opening new avenues for the important topic of data integration in a spatial mapping context.