Spaceborne global navigation satellite system reflectometry (GNSS-R) is an innovative bistatic radar remote sensing technique utilizing low Earth orbit (LEO) based GNSS-R instruments to acquire GNSS L-band opportunistic signals for measuring geophysical parameters. A GNSS-R LEO constellation with an optimization design for its specialized missions is very significant and necessary. However, the constellation design involves multi-parameter and multi-objective optimization, and the classical analytic solution is not capable of such a complicated issue. This study proposes a multi-objective LEO constellation design method with a genetic algorithm (GA) and presents a framework for designing two GNSS-R LEO constellations, termed "lower-latitude constellation" for typhoons and hurricanes observation in the tropics and "global constellation" for global geophysical parameter measurements. Then, the observation capability of both designed constellations is evaluated in terms of the number of reflection points, spatial coverage density, and revisit time to verify the GA efficiency in LEO constellation design. Results show that the two designed LEO constellations with high fitness function values possess optimal orbit parameter set configuration and outperform the existing CyGNSS constellations in observation performance. Compared with CyGNSS, the number of reflection points observed by the lower-latitude constellation and the global constellation increases by 38% and 45%, as well as the spatial coverage density increases by 28% and 36%. The revisit time for the lower-latitude constellation is reduced by 0.29 h, whereas the revisit time for the global constellation increases by one hour.