Many aquatic species have restricted dispersal capabilities, making them the most vulnerable organisms to climate change and land use change patterns. These factors deplete Nymphaea species’ suitable habitats, threatening their populations and survival. In addition, the species are poorly documented, which may indicate how scarce they are or will become. Members of Nymphaea are ecologically important as well as having cultural and economic value, making them of conservation interest. Therefore, using the maximum entropy (MaxEnt) approach, climatic variables, land use, and presence points were modeled for seven Nymphaea species in South America, using three general circulation models (CCSM4, HADGEM2-AO, and MIROC5) and in two representative concentration pathways (RCPs 4.5 and 8.5) and two scenarios (2050 and 2070). Our results indicated that mean diurnal range (bio2), precipitation of the wettest month (bio13), temperature seasonality (bio15), and land use (dom_lu) were the main influencing factors. For all species, suitable areas were concentrated east of Brazil, and they were variable in northern parts of the continent. Besides, inconsistent expansion and contraction of suitable habitats were noticed among the species. For example, N. amazonum, N. rudgeana, and N. lasiophylla future habitat expansions declined and habitat contraction increased, while for N. ampla and N. jamesoniana, both future habitat expansion and contraction increased, and for N. pulchella and N. rudgeana it varied in the RCPs. Moreover, the largest projected suitable habitats were projected outside protected areas, characterized by high human impacts, despite our analysis indicating no significant change between protected and unprotected areas in suitable habitat change. Finally, understanding how climate change and land use affect species distribution is critical to developing conservation measures for aquatic species.