Pollution has become a very serious threat to our environment. Monitoring pollution is the first step toward planning to save the environment. The use of differential equations of monitoring pollution has become possible. In this paper the pollution problem of three lakes with interconnecting channels has been studied. The variational iteration method has been applied to compute an approximate solution of the system of differential equations, governing on the problem. Three different types of input models: sinusoidal, impulse, and step will be considered for monitoring the pollution in the lakes. The results are compared with those obtained by Adomian decomposition method. This comparison reveals that the variational iteration method is easier to be implemented.