Prediction of saturated area formation is important for hydrologic modeling of watersheds with shallow, highly pervious soils. This simulation study examines the relative importance of hillslope properties and rainfall rate on the evolution of saturated source areas during wetting. The study focuses on homogeneous, nonconvergent hillslope topography with constant slope. The two-dimensional, variably saturated ground water model VS2D <Healy, 1990; Lappala et al. 1993> is used to simulate saturated area formation under the action of steady uniform rainfall. The study methodology systematically varies four hillslope properties, depth to impervious layer, slope length, slope angle, and average saturated hydraulic conductivity, and the rainfall rate. Results indicate that the fraction of the hillslope length that is saturated at equilibrium is a function of a parameter &PHgr;, which is defined as the rainfall rate multiplied by the slope length, divided by the slope angle, soil thickness, and saturated hydraulic conductivity. The temporal evolution of saturated area is analyzed in terms of the equilibrium time. Nonlinearity in the unsaturated zone and differences in hillslope properties result in a nonunique time to equilibrium relation. The time to equilibrium is maximum when factors that tend to cause surface saturation are in approximate balance with those that tend to dissipate surface saturation. The temporal evolution of surface saturation during wetting follows a wide variety of trajectories. Equivalent hillslope soil water storage ratios on slopes with different properties can result in a wide range of surface saturation conditions. Hillslopes with shallower soils and smaller slope angles are most susceptible to saturated area formation. However, sudden changes in surface saturation are possible on steep hillslopes when 1<&PHgr;<6. ¿ 2000 American Geophysical Union