The last 85,000 years were characterized by high climate and environmental variability on the Yucatán Peninsula. Heinrich stadials are examples of abrupt climate transitions that involved shifts in regional temperatures and moisture availability. Thus, they serve as natural experiments to evaluate the contrasting responses of aquatic and terrestrial ecosystems. We used ostracodes and pollen preserved in a 75.9-m-long sediment core (PI-6, ~85 ka) recovered from Lake Petén Itzá, Guatemala, to assess the magnitude and velocity of community responses. Ostracodes are sensitive to changes in water temperature and conductivity. Vegetation responds to shifts in temperature and the ratio of evaporation to precipitation. Ostracodes display larger and more rapid community changes than does vegetation. Heinrich Stadial 5-1 (HS5-1) was cold and dry and is associated with lower ostracode and vegetation species richness and diversity. In contrast, the slightly warmer and dry conditions during HS6 and HS5a are reflected in higher ostracode species richness and diversity. Our paleoecological study revealed the greatest ecological turnover for ostracodes occurred from 62.5 to 51.0 ka; for pollen, it was at the Pleistocene/Holocene transition. Future studies should use various climate and environmental indicators from lake and marine sediment records to further explore late glacial paleoclimate causes and effects in the northern neotropics.

Three new species of the genus Cypretta Vávra, 1985 are described from Southern Mexico, representing the first record of the genus in the country. Cypretta campechensis n. sp. is closely related to Florida and North Carolina species such as C. nigra Furtos, 1936, C. brevisaepta Furtos, 1934 and C. bilicis Furtos, 1936. Cypretta spinosa n. sp. is related to the North and South American species C. intonsa Furtos, 1936 and C. vivacis Würdig & Pinto, 1993, but also to some Aus-tralian and South-East Asian Islands ones, such as C. obfuscata Victor & Fernando, 1981 and C. raciborskii (Grochmalicki, 1915). The last species described herein, Cypretta maya n. sp., is closely related to the South-East Asia islands, Australian and north Indian species, in particular to C. hirsuta Henry, 1923, C. longidactyla Victor & Fernando, 1981, and C. patialensis Battish, 1982.

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.

We evaluated how ranges of four endemic and non-endemic aquatic ostracode species changed in response to long-term (glacial–interglacial cycles) and abrupt climate fluctuations during the last 155 kyr in the northern Neotropical region. We employed two complementary approaches, fossil records and species distribution models (SDMs). Fossil assemblages were obtained from sediment cores PI-1, PI-2, PI-6 and Petén-Itzá 22-VIII-99 from the Petén Itzá Scientific Drilling Project, Lake Petén Itzá, Guatemala. To obtain a spatially resolved pattern of (past) species distribution, a downscaling cascade is employed. SDMs were reconstructed for the last interglacial (∼120 ka), the last glacial maximum (∼22 ka) and the middle Holocene (∼6 ka). During glacial and interglacial cycles and marine isotope stages (MISs), modelled paleo-distributions and paleo-records show the nearly continuous presence of endemic and non-endemic species in the region, suggesting negligible effects of long-term climate variations on aquatic niche stability. During periods of abrupt ecological disruption such as Heinrich Stadial 1 (HS1), endemic species were resilient, remaining within their current areas of distribution. Non-endemic species, however, proved to be more sensitive. Modelled paleo-distributions suggest that the geographic range of non-endemic species changed, moving southward into Central America. Due to the uncertainties involved in the downscaling from the global numerical to the highly resolved regional geospatial statistical modelling, results can be seen as a benchmark for future studies using similar approaches. Given relatively moderate temperature decreases in Lake Petén Itzá waters (∼5 ^∘C) and the persistence of some aquatic ecosystems even during periods of severe drying in HS1, our data suggest (1) the existence of micro-refugia and/or (2) continuous interaction between central metapopulations and surrounding populations, enabling aquatic taxa to survive climate fluctuations in the northern Neotropical region.

The Yucatán Peninsula is a karstic region, rich in subterranean environments with a diverse crustacean stygobiont fauna. In order to gain insights into the biological evolution of the subterranean environments of this region, we evaluated the ostracode species composition of caves and cenotes in five independent sampling campaigns (2008, 2013, 2017–2019). Using morphometric analyses, we evaluated inter-population morphological variability; using molecular analysis based on mitochondrial COI and nuclear 18S rDNA, we evaluated genetic differentiation in selected species. The observed fauna is composed of 20 (epigean) species, presenting a lack of strict stygobionts. Morphometric analyses discriminated up to three morphotypes in each of the three most abundant species: Cytheridella ilosvayi, Alicenula sp. and Cypridopsis vidua. High intraspecific morphological variability was found either in shape or size. Phylogenetic analysis based on COI demonstrated the existence of three lineages on C. ilosvayi, with high support (>0.9). The 18S rDNA sequences were identical among individuals of different populations. A lack of congruence between the genetic markers precluded us from postulating speciation in subterranean environments. It is likely that Late Pleistocene—Early Holocene climate variation related to sea level and precipitation was forcing agent for epigean ostracode dominance in subterranean environments of the Peninsula.

In the Yucatán Peninsula, anthropogenic activities such as urbanization and the final disposal of solids and wastewater critically impact aquatic systems. Here, we evaluated the anthropogenic-related environmental alteration of Lake La Sabana, located in the northern limits of one of the main cities of the Mexican Caribbean. We evaluated lake water quality, using physical, chemical, and microbiological indicators, and heavy metals in surficial sediment and fish tissue to evaluate the potential environmental risk. Multivariate analyses revealed that environmental conditions in La Sabana are spatially and temporally heterogeneous. Medium to bad water quality was determined within different basins by the National Sanitation Foundation water quality index, related to the degree of anthropogenic influence at each zone. The central-south zones displayed critical microbiological values largely exceeding national standards. Heavy metals in sediment and fish tissue such as Zn and Hg were relatively low, but Hg concentrations threaten the ecological environment. Incipient wastewater treatment and its final disposal in La Sabana are mainly responsible for the changes in the trophic status and availability of nutrients, which in turn may have promoted changes in the biological structure and aquatic plant invasions. Lake La Sabana can be considered a model of the potential and sequential effects of anthropogenic alterations in the oligotrophic karst tropical aquatic systems in the Yucatan Peninsula.