A new empirical model of plasmapause location as functions of magnetic local time and geomagnetic indices has been developed based on the observations from THEMIS mission. We use the two-year data of electron density inferred from spacecraft potential to identify the plasmapause crossings and create a database of plasmapause locations. The database is further used to build up an empirical model of plasmapause related to magnetic local time based on the equation from O’Brien and Moldwin(2003). The new model is compared with previous plasmapause location models. It is found that our newly developed model is the best in predicting plasmapause locations among the existing models. The models based on Kp and Dst indices are better than the model based on AE index, suggesting that the plasmapause location is controlled by large scale convection of the magnetosphere. A new empirical model of plasmapause location as functions of magnetic local time and geomagnetic indices has been developed based on the observations from THEMIS mission. The use of two-year data of electron density inferred from spacecraft potential to identify the plasmapause crossings and create a database of plasmapause locations. The database is further used to build up an empirical model of plasmapause related to magnetic local time based on the equation from O’Brien and Moldwin (2003). The new model is compared with previous plasmapause location models. that our newly developed model is the best in predicting plasma pause locations among the existing models. The Models based on Kp and Dst indices are better than the model based on AE index, suggesting that the plasma pause location is controlled by large scale convection of the magnetosphere.