Empirical Bayesian Kriging, a Robust Method for Spatial Data Interpolation of a Large Groundwater Quality Dataset from the Western Netherlands

No single spatial interpolation method reigns supreme for modelling the precise spatial distribution of groundwater quality data. This study addresses this challenge by evaluating and comparing several commonly used geostatistical methods: Local Polynomial Interpolation (LPI), Ordinary Kriging (OK), Simple Kriging (SK), Universal Kriging (UK), and Empirical Bayesian Kriging (EBK). We applied these methods to a vast dataset of 3033 groundwater records encompassing a substantial area (11,100 km²) in the coastal lowlands of the western Netherlands. To our knowledge, no prior research has investigated these interpolation methods in this specific hydrogeological setting, exhibiting a range of groundwater qualities, from fresh to saline, often anoxic, with high natural concentrations of PO₄ and NH₄. The prediction performance of the interpolation methods was assessed through statistical indicators such as root means square error. The findings indicated that EBK outperforms the other geostatistical methods in forecasting groundwater quality for the five variables considered: Cl, SO₄, Fe, PO₄, and NH₄. In contrast, SK performed worst for the species except for SO₄. We recommend not using SK to interpolate groundwater quality species unless the data exhibit low spatial variation, high sample density, or evenly distributed sampling.