\sin \theta_{ij}$, and $c_{ij} = \cos \theta_{ij}$. Here we use the units of $\hbar = c = 1$.

As shown in [@af03], there are two MSW resonances and three RSF resonances. The resonance densities are written as $$\begin{aligned}
\rho_{res}({\rm MSW-H}) &=&
\frac{m_u |\Delta m^2_{31}| \cos2\theta_{13}}{2\sqrt{2} G_F E_\nu}
\frac{1}{Y_e} , \\
\rho_{res}({\rm MSW-L}) &=&
\frac{m_u \Delta m^2_{21} \cos2\theta_{12}}{2\sqrt{2} G_F E_\nu}
\frac{1}{Y_e} , \\
\rho_{res}({\rm RSF-H}) &=&
\frac{m_u |\Delta m^2_{31}| \cos2\theta_{13}}{2\sqrt{2} G_F E_\nu}
\frac{1}{|1-2Y_e|} , \\
\rho_{res}({\rm RSF-L}) &=&
\frac{m_u \Delta m^2_{21} \cos2\theta_{12}}{2\sqrt{2} G_F E_\nu}
\frac{1}{|1-2Y_e|} , \\
\rho_{res}({\rm RSF-X}) &=&
\frac{m_u |\Delta m^2_{31}| \cos^2\theta_{13}}{2\sqrt{2} G_F E_\nu}
\frac{1}{1-Y_e} ,\end{aligned}$$ where suffixes \`\`H” and \`\`L” correspond to the high and low resonance densities deduced using $|\Delta m^2_{31}|$ and $\Delta m^2_{21}$, respectively. When $\sin^22\theta_{13}$ is larger than a critical value, the RSF-X resonance changes to the RSF-E resonance [@af03]. The RSF-E resonance appears in place of RSF-X resonance owing to a nonlinear dependence of the effective squared masses on the density. We do not distinguish the RSF-E resonance from the RSF-X and use a common term, the RSF-X resonance, in this study.

Most of neutrino oscillation parameters have been determined precisely by recent neutrino experiments [@sno04; @sk06; @kl08]. We use the squared mass differences 