[1] In this paper we derive analytical solutions to statistical moments for transient saturated flow in two-dimensional, bounded, randomly heterogeneous porous media. By perturbation expansions, we first derive partial differential equations governing the zerothorder head h and the first-order head term h, where orders are in terms of the standard deviation of the log transmissivity. We then so...