ABSTRACT We investigate two dark energy cosmological models (i.e., the ΛCDM and ϕCDM models) with massive neutrinos assuming two different neutrino mass hierarchies in both the spatially flat and non-flat scenarios, where in the ϕCDM model the scalar field possesses an inverse power-law potential, V(ϕ) ∝ ϕ −α (α > 0). Cosmic microwave background data from Planck 2015, baryon acoustic oscillation data from 6dFGS, SDSS-MGS, BOSS-LOWZ and BOSS CMASS-DR11, the joint light-curve analysis compilation of SNe Ia apparent magnitude observations, and the Hubble Space Telescope H 0 prior, are jointly employed to constrain the model parameters. We first determine constraints assuming three species of degenerate massive neutrinos. In the spatially flat (non-flat) ΛCDM model, the sum of neutrino masses is bounded as Σm ν < 0.165(0.299) eV at 95% confidence level (CL). Correspondingly, in the flat (non-flat) ϕCDM model, we find Σm ν < 0.164(0.301) eV at 95% CL. The inclusion of spatial curvature as a free parameter results in a significant broadening of confidence regions for Σm ν and other parameters. In the scenario where the total neutrino mass is dominated by the heaviest neutrino mass eigenstate, we obtain similar conclusions to those obtained in the degenerate neutrino mass scenario. In addition, the results show that the bounds on Σm ν based on two different neutrino mass hierarchies have insignificant differences in the spatially flat case for both the ΛCDM and ϕCDM models; however, the corresponding differences are larger in the non-flat case.