POD-based Monte Carlo approach for the solution of regional scale groundwater flow driven by randomly distributed recharge