We study the eigenvalue problem of the Laplacian on a thin domain under a mixed boundary condition by using a pure variational approach. We obtain the detailed asymptotic behavior of such eigenvalues when the domain becomes thinner and shrinks to a surface, and we establish the characteristic geometric dependency. Furthermore, we also present an application to a certain bifurcation problem arising in population dynamics.