Abstract. Multiple discrete-energy ion bands observed by the Polar satellite in the inner magnetosphere on 9 February 1998 were investigated by means of particle simulation with a realistic model of the convection electric field. The multiple bands appeared in the energy vs. L spectrum in the 1-100 keV range when Polar traveled in the heart of the ring current along the outbound and inbound paths. We performed particle tracing, and simulated the energy vs. L spectra of proton fluxes under the dipole magnetic field, the corotation electric field, and the realistic convection electric field model with its parameters depending on the solar wind data. Simulated spectra are shown to agree well with the observed ones. A better agreement is achieved when we rotate the convection electric potential eastward by 2 h in MLT and we change the distribution function in time in the nearEarth magnetotail. It is concluded that the multiple bands are likely produced by two processes for this particular event, that is, changes in the convection electric field (for >3 keV protons) and changes in the distribution function in the nearEarth magnetotail (for <3 keV protons).