A model for the transient one-dimensional moisture movement in the saturated-unsaturated zone using a finite difference method is developed. Hysteresis in the soil water retention is incorporated. The model considers layered geologic formations. Monte Carlo simulation, together with the nearest neighbor model is used. Outputs of the model include pressure head, water content, and the water table elevation. Two Monte Carlo simulations of 100 realizations each are made for a 12-day simulation period with different input values. The simulation results show that the S.D. of the outputs increases with an increase in the input, the S.D. of the log K$$. The model is applied to predict a long term water table fluctuation, and the predicted water table agress well with the observed one.