Breakthrough and terminal head loss are the main parameters that determine the performance of rapid sand filters. Carman-Kozeny and Ergun equations can be applied to estimate head loss, but can only be applied to clean filter beds. Elaborated models are needed to predict head loss in dirty filters. In this study, a neuro-fuzzy modeling approach was proposed to estimate head loss in dirty filters. Hydraulic loading rate, influent iron concentration, bed porosity, and operating time were selected as input variables. Various types of membership functions were tried. Two rule-base generation methods-subtractive clustering and grid partition were used for a first-order, Sugeno-type inference system. Using 11 rules and the grid-partition method, an optimum rule base set was developed and the lowest root mean squared error (RMSE) was obtained.