Low Earth orbit (LEO) satellite constellations have emerged as an effective alternative for the provision of high-accuracy positioning, navigation and timing (PNT) solutions which are based on high-precision orbit and clock information. Determining an orbit with high precision is dependent on the number and distribution of ground tracking stations. Therefore, it is important to investigate methodologies that can ensure the adequate observing coverage of LEO navigation constellations. In this study, an evolutionary algorithm is applied to optimize the number and deployment of ground stations for tracking LEO constellations. According to the distribution area, two schemes of study are analyzed: (a) global deployment—the ground stations are deployed throughout the globe; (b) regional deployment—a selected region is used for deployment. For global deployment, the optimization objectives are focused on the ground station and observing rate for k-heavy observing coverage (HC), while the sole objective for the regional deployment scheme is the satellite position dilution of precision (SPDOP). It is shown that a deployment of 95 ground stations is optimal for achieving 3-HC with an observing rate of 97.37% and 4-HC with an observing rate of 92.01%. For regional distribution, 15, 20 and 25 ground stations are used for three optimal configurations of SPDOP at 2.058, 1.399 and 1.330, respectively. The results are significantly enhanced using intersatellite links for SPDOP evaluation, from 2.058, 1.399 and 1.330 to 0.439, 0.422 and 0.409, with 15, 20 and 25 ground stations, respectively.