Mathematical models based on nonlinear programming are developed for operating recharge basin systems. The objective of these optimization models is to determine the operation policy that maximizes the infiltration subject to constraints for continuity, infiltration, groundwater flow, and the physical constraints describing the recharge basins. Two basic physical processes that are involved in the recharge basin systems are the infiltration process and the soil moisture redistribution process. The infiltration process is modeled using both the Philip's (1957) and the Green-Ampt (1911) infiltration models and the soil moisture redistribution process is modeled using the Darcy's equation as well as a kinematic wave model (Charbeneau, 1984). An application is presented.