When modelling soil acidification at the European scale, it is inevitable that both model and data have varying degrees of associated uncertainty. The present study attempted to quantify the uncertainty in long-term forecasts of soil solution concentrations of Al{sup 3+} and No{sub 3}{sup {minus}} concentrations resulting from the uncertainty in small resolution European-scale maps and input data, using the Netherlands as a case study. Large-scale forecasts were made with a relatively simple dynamic process-oriented model, SMART2. Model outputs were considered as block median concentrations and the block areal fractions in which concentrations exceeded a critical level. Sources of uncertainty considered included (1) uncertainty in soil and vegetation maps (categorical data), and (2) uncertainty in soil and vegetation-related parameters (continuous data). The uncertainty in model outputs was quantified by an efficient two-step Monte Carlo simulation approach, which takes spatial correlation into account. The uncertainty in the input data at the European scale led to major uncertainties in the predicted Al{sup 3+} concentration. Uncertainties in the areas where the Al{sup 3+} concentration exceeded the maximum allowable concentration were much smaller. The uncertainties in soil-related parameters contributed most to the uncertainty in the Al{sup 3+} concentration, whereas the uncertainty contributed by the soilmore » and vegetation maps was negligible. For the NO{sub 3}{sup {minus}} concentration, however, the uncertainty originated mainly from the soil and vegetation maps. Evaluation of the different error sources is of great practical significance, as it identified which sources need further improvement. The present study shows that the uncertainty contribution of the different error sources depends greatly on the model output considered.« less