Despite significant advancements in our understanding of the relationship between soil animal life and environmental gradients, as well as their impact on ecosystem functioning, there is a growing need for quantitative research that considers the spatial and temporal heterogeneity of soil biodiversity. Such research can help to address current knowledge gaps and support the development of legislation and conservation measures aimed at preserving soil biodiversity in the future. Earthworms, one of the most important group of soil animals, show patchy horizontal distribution in response to spatial heterogeneity in soil physico-chemical parameters, food availability and biological interactions at various scales.The aim of this study was to explore whether the spatial and temporal variability and co-occurrence probability of earthworm species and ecological categories (epigeic, endogeic and anecic) can be characterized by exploiting habitat micro-topography and insolation.Field research was carried out in spring and summer at seven locations using a regular grid of 49 sampling points with 5 m inter-sample distance over an altitudinal gradient. The results of Geary's C and Moran's I indices, as well as the semivariograms calculated for each species, season and location confirmed a patchy spatial distribution. The range of semivariograms varied from 3.2 m to 32.8 m in spring, and from 0.3 m to 32.8 m in summer. Overall (all locations and both seasons), the mean range was similar for adult (8.4 m) and juvenile (9.2 m) earthworms. In spring, epigeic species had a larger range than endogeic species. Abundances of juvenile earthworms were negatively correlated with microtopographic elevation in spring (at most locations) and in summer (at all locations). The relationship between adult earthworms and microtopographic elevation could not be interpreted generally within season. Additionally, correlations between earthworm abundance and average monthly insolation were weak. Pairwise analysis of species co-occurrence revealed a random pattern for most species. In the analysis of co-occurrence of ecological categories, endogeics and epigeics had more positive probabilities than other pairs. The obtained results could be useable in soil biodiversity assessments, or creation of monitoring protocols.