In this study we develop a first-order, nonstationary stochastic model for steady state, unsaturated flow in randomly heterogeneous media. The model is applicable to the entire domain of a bounded vadose zone, unlike most of the existing stochastic models. Because of its nonstationarity, we solve it by the numerical technique of finite differences, which renders the flexibility in handling diff...