Assuming that there exist at least two fermionic parameters, the classical N = 1 supersymmetric Korteweg-de Vries (SKdV) system can be transformed to some coupled bosonic systems. The boson fields in the bosonized SKdV (BSKdV) systems are defined on even Grassmann algebra. Due to the intrusion of other Grassmann parameters, the BSKdV systems are different from the usual non-suppersymmetric integrable systems, and many more abundant solution structures can be unearthed. With the help of the singularity analysis, the Painlevé property of the BSKdV system is proved and a Bäcklund transformation (BT) is found. The BT related nonlocal symmetry, we call it as residual symmetry, is used to find symmetry reduction solutions of the BSKdV system. Hinted from the symmetry reduction solutions, a more generalized but much simpler method is established to find exact solutions of the BSKdV and then the SKdV systems, which actually can be applied to any fermionic systems.

Bosonization approach to the classical supersymmetric systems is presented. By introducing the multifermionic parameters in the expansions of the superfields, the N = 1 supersymmetric KdV (sKdV) equations are transformed to a system of coupled bosonic equations. The method can be applied to any fermionic systems. By solving the coupled bosonic equations, some novel types of exact solutions can be explicitly obtained. Especially, the richness of the localized excitations of the supersymmetric integrable system are discovered. The rich multi-soliton solutions obtained here have not yet been obtained by using other methods. Unfortunately, the traditional known multi-soliton solutions can also not be obtained by the bosonization approach of this paper. Some open problems on the bosonization of the supersymmetric integrable models are proposed in the both classical and quantum levels.