Population Density and Diversity of Millipedes in Four Habitat Classes: Comparison Concerning Vegetation Type and Soil Characteristics

Our study was conducted in the Valle Nuevo National Park and included four habitat classes: tussock grass (Sabapa), pine forest (Pinoc), broadleaf forest (Boslat), and agricultural ecosystem (Ecoag). We had two main objectives: to comparatively describe millipede communities and to determine the relationships between population density/diversity and soil physicochemical variables. The research was cross-sectional and non-manipulative, with a descriptive and correlational scope; sampling followed a stratified systematic design, with eight transects and 32 quadrats of 1 m2, covering 21.7 km. We found a sandy loam soil with an extremely acidic pH. The highest population density of millipedes was recorded in Sabapa, and the lowest in Ecoag. The highest alpha diversity was shared between Boslat (Margalef = 1.72) and Pinoc (Shannon = 2.53); Sabapa and Boslat showed the highest Jaccard similarity (0.56). The null hypothesis test using the weighted Shannon index revealed a statistically significant difference in diversity between the Boslat–Sabapa and Pinoc–Sabapa pairs. Two of the species recorded highly significant indicator values (IndVal) for two habitat classes. We found significant correlations (p < 0.05) between various soil physicochemical variables and millipede density and diversity.