Non-point source contaminations of aquifers are complex to model, predict and manage. This study uses numerical and stochastic methods to address the importance of key sources of spatial variability. We show that heterogeneity in recharge and contaminant loading does not significantly impact management metrics and could be simplified. Also, homogenizing physical properties has more impact on predictions, but can provide useful information on concentration statistics in a regional analysis.