Due to the great success of the CYclone Global Navigation Satellite System (CYGNSS) mission, the follow-on GNSS Reflectometry (GNSS-R) missions are being planned. In the perceivable future, signal sources for GNSS-R missions can originate from multiple global navigation satellite systems (GNSSs) including Global Positioning System (GPS), Galileo, GLONASS, and BeiDou. On the other hand, to facilitate the operational capability for sensing ocean, land, and ice features globally, multi-satellite low Earth orbit (LEO) constellations with global coverage and high spatio-temporal resolutions should be considered in the design of the follow-on GNSS-R constellation. In the present study, the particle swarm optimization (PSO) algorithm was applied to seek the optimal configuration parameters of 2D-lattice flower constellations (2D-LFCs) composed of 8, 24, 60, and 120 satellites, respectively, for global GNSS-R observations, and the fitness function was defined as the length of the time for the percentage coverage of the reflection observations reaches 90% of the globe. The configuration parameters for the optimal constellations are presented, and the performances of the optimal constellations for GNSS-R observations including the visited and the revisited coverages, and the spatial and temporal distributions of the reflections were further compared. Although the results showed that all four optimized constellations could observe GNSS reflections with proper temporal and spatial distributions, we recommend the optimal 24-and 60-satellite 2D-LFCs for future GNSS-R missions, taking into account both the performance and efficiency for the deployment of the GNSS-R missions.