This paper proposes an Integrated Monte Carlo Methodology (IMCM) to solve the parameter estimation problem in water quality models. The methodology is based on Bayesian approach and Markov Chain Monte Carlo techniques and it operates by means of four modules: Markov Chain Monte Carlo (MCMC), Moving Feasible Ranges (MFR), Statistical Analysis of the Joint Posterior Distribution (SAD) and Uncertainty Propagation Analysis (UPA). The main innovation of the new proposal lies in the combination of MCMC and MFR modules which provides the joint posterior distribution of the calibrated parameters following the classical Bayesian approach. While MCMC module, based on Shuffled Complex Evolution Metropolis (SCEM-UA) algorithm, is specially designed to sample complex joint posterior shapes within certain parameter ranges, the MFR readjusts these ranges until the coverage of the feasible parameter space is guaranteed. Once the joint posterior distribution is properly defined, the SAD provides the parameter statistics and the UPA performs an analysis of the uncertainty propagation through the model. The possibilities of the new proposal have been tested on the basis of a simple model featuring different activated sludge batch experiments. IMCM has been implemented in Matlab and it is prepared to be easily connected to any software package. (C) 2010 Elsevier B.V. All rights reserved.