A large database including temporal trends of physical, ecological and socio-economic data was developed within the EUROCAT project. The aim was to estimate the nutrient fluxes for different socio-economic scenarios at catchment and coastal zone level of the Po catchment (Northern Italy) with reference to the Water Quality Objectives reported in the Water Framework Directive (WFD 2000/60/CE) and also in Italian legislation. Emission data derived from different sources at national, regional and local levels are referred to point and non-point sources. While non-point (diffuse) sources are simply integrated into the nutrient flux model, point sources are irregularly distributed. Intensive farming activity in the Po valley is one of the main Pressure factors Driving groundwater pollution in the catchment, therefore understanding the spatial variability of groundwater nitrate concentrations is a critical issue to be considered in developing a Water Quality Management Plan. In order to use the scattered point source data as input in our biogeochemical and transport models, it was necessary to predict their values and associated uncertainty at unsampled locations. This study reports the spatial distribution and uncertainty of groundwater nitrate concentration at a test site of the Po watershed using a probabilistic approach. Our approach was based on geostatistical sequential Gaussian simulation used to yield a series of stochastic images characterized by equally probable spatial distributions of the nitrate concentration across the area. Post-processing of many simulations allowed the mapping of contaminated and uncontaminated areas and provided a model for the uncertainty in the spatial distribution of nitrate concentrations.