We propose a supersymmetric gradient flow equation in the four-dimensional Wess-Zumino model. The flow is constructed in two ways. One is based on the off-shell component fields and the other is based on the superfield formalism, in which the same result is provided. The obtained flow is supersymmetric because the flow time derivative and the supersymmetry transformation commute with each other. Solving the equation, we find that it has a damping oscillation with the flow time for nonzero mass, which is different from the Yang-Mills flow. The on-shell flow equation is also discussed.