In this paper, we deal with modeling crisp input-fuzzy output data by constructing a MARS-fuzzy regression model with crisp parameters estimation and fuzzy error terms for the fuzzy data set. The proposed method is a two-phase procedure which applies the MARS technique at phase one and an optimization problem at phase two to estimate the center and fuzziness of the response variable. A realistic application of the proposed method is also presented in a hydrology engineering problem. Empirical results demonstrate that the proposed approach is more efficient and more realistic than some traditional least-squares fuzzy regression models.