15Soil mesofauna communities are hyperdiverse and critical for ecosystem functioning. 16 However, our knowledge on spatial structure and underlying processes of community 17 assembly for soil arthropods is scarce, hampered by limited empirical data on species 18 diversity and turnover. We implement a high-throughput-sequencing approach to 19 generate comparative data for thousands of arthropods at three hierarchical levels: 20 genetic, species and supra-specific lineages. A joint analysis of the spatial arrangement 21 across these levels can reveal the predominant processes driving the variation in 22 biological assemblages at the local scale. This multi-hierarchical approach was performed 23 using haplotype-level-COI metabarcoding of entire communities of mites, springtails and 24 beetles from three Iberian mountain regions. Tens of thousands of specimens were 25 2 extracted from deep and superficial soil layers and produced comparative 26 phylogeographic data for >1000 co-distributed species and nearly 3000 haplotypes. Local 27 assemblages were highly distinctive between grasslands and forests, and within each of 28 them showed strong spatial structures and high endemicity at the scale of a few kilometres 29 or less. The local distance-decay patterns were self-similar for the haplotypes and higher 30 hierarchical entities, and this fractal structure was very similar in all three regions, 31 pointing to a significant role of dispersal limitation driving the local-scale community 32 assembly. Our results from whole-community metabarcoding provide unprecedented 33 insight into how dispersal limitations constrain mesofauna community structure within 34 local spatial settings over evolutionary timescales. If generalized across wider areas, the 35 high turnover and endemicity in the soil locally may indicate extremely high richness 36 globally, challenging our current estimations of total arthropod-diversity on Earth. 37 38 39