Supersymmetry, a symmetry between fermions and bosons, provides a promising extension of the standard model but is still lack of experimental evidence. Recently, the interest in supersymmetry arises in the condensed matter community owing to its potential emergence at the continuous quantum phase transition. In this work, we demonstrate that 2+1D supersymmetry, relating massive Majorana and Ising fields, might emerge at the first-order quantum phase transition of the Ising magnetization by tuning a single parameter. Although the emergence of the SUSY is only allowed in a finite range of scales due to the existence of relevant masses, the scale range can be large when the masses before scaling are small. We show that the emergence of supersymmetry is accompanied by a topological phase transition for the Majorana field, where its non-zero mass changes the sign but keeps the magnitude. An experimental realization of this scenario is proposed using the surface state of a 3+1D time-reversal invariant topological superconductor with surface magnetic doping.