Geodiversity is recognized as one of the most important drivers of ecosystem characteristics and biodiversity globally. However, in the northern Neotropics, the contribution of highly diverse landscapes, environmental conditions, and geological history in structuring large-scale patterns of aquatic environments and aquatic species associations remains poorly understood. We evaluated the relationships among geodiversity, limnological conditions, and freshwater ostracodes from southern Mexico to Nicaragua. A cluster analysis (CA), based on geological, geochemical, mineralogical, and water-column physical and chemical characteristics of 76 aquatic ecosystems (karst, volcanic, tectonic) revealed two main limnological regions: (1) karst plateaus of the Yucatán Peninsula and northern Guatemala, and (2) volcanic terrains of the Guatemalan highlands, mid-elevation sites in El Salvador and Honduras, and the Nicaraguan lowlands. In addition, seven subregions were recognized, demonstrating a high heterogeneity of aquatic environments. Principal component analysis (PCA) identified water chemistry (ionic composition) and mineralogy as most influential for aquatic ecosystem classification. Multi-parametric analyses, based on biological data, revealed that ostracode species associations represent disjunct faunas. Five species associations, distributed according to limnological regions, were recognized. Structural equation modeling (SEM) revealed that geodiversity explains limnological patterns of the study area. Limnology further explained species composition, but not species richness. The influence of conductivity and elevation were individually evaluated in SEM and were statistically significant for ostracode species composition, though not for species richness. We conclude that geodiversity has a central influence on the limnological conditions of aquatic systems, which in turn influence ostracode species composition in lakes of the northern Neotropical region.
