Introduction

The area of forests in Poland increased between 1945 and 2012 from 21 to 30% (i.e. from 6,470,000 ha to 9,164,000 ha). Such growth occurred primarily through the afforestation of post-agricultural lands (Kaliszewski 2016), accounting for 55% of afforestation—1,477,000 ha (Krawczyk 2015). Afforested post-agricultural land differs in physical and chemical characteristics from typical forest land. The most important differences include a higher pH and a lower C:N ratio. This results from long-term agricultural use—repeated ploughing and fertilization. Ploughing is responsible for the transformation of the soil profile and fertilization for the chemical changes (Gorzelak 1996). Additional differences are generated by the fact that these lands have been deprived for many years of the influence of tree root systems, which stabilize the soil structure and enrich it biologically (Thoms et al. 2010). The greater biodiversity of bacteria and lower diversity of fungi in soils of post-agricultural lands may be triggered by higher availability of nutrients (Delgado-Baquerizo et al. 2017). Planting of trees on post-agricultural land shapes the succession process of herbaceous plants and shrubs. This process is long-lasting and depends on the fragmentation of habitats and the vicinity of inhabited stands (Ciurzycki et al. 2021). The diversity of different plant functional groups may be positively correlated with soil functions such as microbial biomass carbon, basal respiration, and carbon sequestration. Thus, with the ongoing succession of these groups on post-agricultural land, the parameters may improve (Heydari et al. 2020). Both land use history and soil type are, therefore, a set of abiotic factors that shape the living conditions of trees and soil organisms. The importance of soil biota in soil remediation processes is based on their contribution to the decomposition of dead organic matter. Fungi, bacteria, nematodes (Nematoda), springtails (Collembola) and soil mites (Acari) should be considered as the main groups of such organisms. Good bioindicators of changes in the soil environment are free-living mesostigmatid mites (Mesostigmata). They occur in large numbers in the soil and forest litter, lead a relatively sedentary lifestyle (Kamczyc et al. 2019) and are sensitive to pollution and degradation of forest areas (Manu et al. 2019). They hunt primarily on springtails, nematodes, potworms (Enchytraeidae) or insect larvae (Insecta) (Koehler 1999; Wissuwa et al. 2012). The structure of mesostigmatid mite communities depends on the particular conditions of habitat type, such as vegetation, age of the trees or human impact and abiotic factors. Mesostigmatid mites highlight the degree of anthropogenization in ecosystems, but this effect is especially pronounced in forests (Călugăr, 2021). Mesostigmatid mite communities are also subject to seasonal changes in species composition and abundance concerning temperature and precipitation (Salmane 2000; Kamczyc et al. 2022). Environmental conditions and seasonal dynamics are also expressed through changes in the relationships of different mesostigmatid developmental instars, which include adult (males, females) and three juvenile instars (deutonymphs, protonymphs and larvae) (Urbanowski et al. 2021).

The aim of the study was to recognize differences in seasonal changes in the abundance, species richness and diversity of Mesostigmata mite assemblages inhabiting pure forests (Betula pendula Roth, Tilia cordata Mill., Pinus sylvestris L.) growing on post-agricultural lands. Considering previous studies on post-agricultural lands (Scheu and Schulz 1996; Gormsen et al. 2006; Gawęda et al. 2021), we hypothesized that abundance, species richness, diversity, and relative proportion of different instars (including sex ratio) in mite communities would respond in different ways to (1) season (summer and autumn), (2) tree species, and (3) soil properties. We expected that various forest litter types would affect soil mite assemblages differently, and that effect would change with season.

Materials and methods

Site description and experimental design

The study site was located in the Opole Forest District (S–W Poland), where coniferous forests cover 60% of the forest area. The forests are dominated by mesic (57% of area), moist (40%), and marsh sites (~ 3%). Considering the soil environment, the main groups of soils are Arenosols (rusty soils)—45.3%, Podzols (podzolic soils)—23.0%, and Gleysols (gleyic soils)—16.4% (Forest Management Plan for Opole Forest Division, 2014). The vegetation season lasts 227 days. Mean annual precipitation reaches 603 mm, while mean annual temperature is 8 °C. The duration of winter is 60–70 days, and the number of days with snow cover is 58–65 days. In turn, the duration of summer is 90–100 days (Malica et al. 2024).

At the beginning of our study, three pure stands of different tree species (Pinus sylvestris L., Tilia cordata Mill. and Betula pendula Roth) were selected (Fig. 1). Every stand was located on post-agricultural land characterized by sandy, acidic or slightly acidic soils, with no calcium carbonates and low soil organic matter content in upper soil horizons (rusty soil—Arenosol). The distance between stands was at least 50 m. Finally, eight circular study plots (0.02 ha) were established in total. The study site characteristics were done in summer 2021 to note all plant species growing on each site. On each study plot, we measured: (1) the total number of trees (N); (2) the height of the trees (m); (3) the dimension at breast height—DBH (cm). Moreover, all vascular plants were recorded and the cover of each species in each layer was estimated using the seven-level Braun-Blanquet abundance scale (Table 4 in Appendix 1).

Fig. 1
figure 1

Study sites on post-agricultural lands (July 2021) in the Opole Forest District: A Betula pendula stand, B Pinus sylvestris stand, C Tilia cordata stand

Soil analyses

The soil environment was characterized by the description of soil pit (at least to a depth of 1 m) located in the middle of the chosen stands. Subsequently, we described the reference soil groups based on the IUSS WRB (2015) and took soil and litter (mineral and organic layer separately) samples (ca. 500 g in total) from every soil horizon for further laboratory studies (16 samples in total). Soil samples were collected to determine: (1) soil texture (Casagrande’s aerometric method modified by Prószyński); (2) pH of soil and litter—measured in distilled water (potentiometric method); (3) content of organic carbon (Corg%—Tiurin’s method); (4) soil organic matter (SOM%—loss on ignition method); (5) the total nitrogen content (N%—Kjeldahl’s method); (6) C to N ratio; (7) elements such as Mg, Na, Ca, P, K (%) as well as Fe, Cd, Pb, Mn, Zn, and Cu (mg/kg) using atomic absorption spectrometry analysis using the AAS Varian 55B spectrometer (Peña et al. 2016). Moreover, we collected samples of intact soil for analysis of bulk density (BD; g/cm3)—core method, and actual soil moisture (SM; %)—gravimetric method (Lityński et al. 1976). The analyses were conducted in the laboratory at the Poznań University of Life Sciences (Table 5 in Appendix 2).

Mesostigmata mites investigation

Sampling was carried out four times: in July and October 2021 and in July and October 2022. Collecting of the soil samples (288 samples in total = eight plots × nine samples per plot × four samplings) was conducted on each study plot using metal soil corer (ø 5 cm) to the depth of 15 cm. In the laboratory, samples were placed on the Berlese-Tullgren apparatus, consisting of a funnel, a light bulb (40W), a strainer (mesh size 2 mm), and containers with 75% ethanol solution. Edaphon was extracted from the soil samples for at least seven days (until the collected samples were dry). Mesostigmatid mite individuals were selected from the alcohol solution under a stereomicroscope and placed in the Hoyer’s medium on slides. The detailed diagnosis took place under a compound microscope Zeiss Axio Scope.A1. All individuals were classified into the appropriate instar, as well as species level or higher taxonomic unit using identification keys (Micherdziński 1969; Karg 1971, 1993; Gwiazdowicz 2007; Mašán & Halliday 2014).

Data analysis

Sampled material allowed to determine mite abundance (ind.), species richness and Shannon–Wiener diversity index (H′ = − Σpiln(pi), where pi is the proportion of particular species in the mite community) for each of the plots. During analysis raw, data for abundance, species richness and diversity of mesostigmatid mite communities were presented as mean (± SE) values per sample within collection month and habitat type. In addition, for each species, indicators of dominance and constancy (the ratio of samples in which a given species occurs to all samples collected in a given variant) were determined for each habitat variant. Statistical analyses were done using R software (R version 4.1.2; R Core Team R: A Language and Environment for Statistical Computing; available online: https://www.r-project.org/). We used generalized linear models (GLM) assuming a Tweedie distribution for mesostigmatid abundance, species richness, and diversity calculated per sample and per plot. We applied predictors in models with low variance inflation factors (VIF < 5). Furthermore, we conducted Tukey post-hoc tests for each model to determine the differences in studied factors between month/habitat combinations. The statistical significance of variables used in GLM’s was calculated using t-values. Results were treated as significant when p < 0.05. In order to describe the relationship between mesostigmatid mite assemblages and habitat characteristics, including soil properties, we conducted canonical correspondence analysis (CCA) implemented in the vegan package in R. The variables included in the CCA final model were based on Akaike information criterion (AIC). Furthermore, environmental variables were selected on VIF < 2. To test the factor's significance, we conducted a permutational analysis of variance (PERMANOVA). We plotted a cumulative species number using the vegan::specaccum() function (Oksanen et al. 2018). Furthermore, we used the bipartite package to reveal the relationship between mesostigmatid mite taxa and habitats (Dormann et al. 2008). In the conducted bipartite interaction network we assumed habitats as lower-level and mite taxa as higher-level groups. The taxa level response to specific habitat was described by specialisation index d′, derived from the Kulback-Leibler distance, expressing whether a given taxon is a non-specialist or a perfect specialist (range from 0 to 1). Additionally, d′ index describes how strongly the taxa differ from a random sampling of available interacting partners (Dormann 2011).

Results

Overall results

In total, 399 individuals (266 females, 50 males and 83 juveniles) were classified as 38 taxa (33 species, five genera) (Table 1). The GLM models assuming Tweedie distribution revealed that collection month and Fe content in soil affected abundance, species richness and diversity. Moreover, mesostigmatid mite abundance was influenced by litter N content (Table 2). The canonical correspondence analysis (CCA) revealed significant differences in mite taxa distribution. Axis 1 and axis 2 explained 26.1% and 25.0% of variance in community composition, respectively (Fig. 2). Axis 1 reflects collection time, while Axis 2 reflects the tree species studied. The analysis of variance conducted for the model showed that collection month (p = 0.002), year (p = 0.002), soil Mn content (p = 0.001), soil Na content (p = 0.028) and litter thickness (p = 0.036) significantly influenced the mite communities structure (Table 3). Species of the genus Paragamasus were influenced by soil Mn content and litter thickness. In contrast, species belonging to the genus Veigaia were affected by soil Na content and were also the least sensitive to soil Mn among all other genera.

Table 1 Species composition, dominance, and constancy of occurrence of individual mite species in each type of tree stand
Table 2 Generalized linear models assuming a Tweedie distribution explaining mite abundance, species richness, and Shannon’s diversity
Fig. 2
figure 2

Canonical correspondence analysis (CCA) for Mesostigmata communities across four sampling periods (July, October; 2020, 2021) in pure stands on post-agricultural lands (Opole Forest District). Mite taxa labelled with the first two letters of the genus name and the first two letters of the species name for each mite taxon (for abbreviations of mite taxa see Table 6 in Appendix 3). BPB. pendula, PSP. sylvestris, TCT. cordata, Jul—July, Oct—October, FLt—forest litter thickness, Sna—soil Na content, Smn—soil Mn content

Table 3 The influence of environmental variables on Mesostigmata communities (PERMANOVA)

Seasonal changes in mesostigmatid mite communities

The highest abundance calculated per sample was recorded in P. sylvestris stand in October 2022 (3.19 ± 0.62 ind.), while the lowest was recorded in P. sylvestris stand in July 2022 (0.37 ± 0.17). The highest abundance in July was reported in T. cordata stand in 2021 (1.59 ± 0.50 ind.) (Fig. 3A). The Tukey post-hoc tests revealed that mesostigmatid mite abundance was significantly lower in P. sylvestris stand in July (0.57 ± 0.17 ind.) than in P. sylvestris and T. cordata stands in October (2.17 ± 0.41 and 2.15 ± 0.34, respectively). Moreover, abundance in P. sylvestris stand in October was significantly higher than in B. pendula stand in July (0.78 ± 0.18 ind.).

Fig. 3
figure 3

Mean abundance (A), species richness (B) and Shannon’s diversity (C) of mesostigmatid mites in the total communities across four sampling periods (July, October; 2020, 2021) in pure stands on post-agricultural lands (Opole Forest District). Abbreviations of variants: Oct—October, Jul—July, BPB. pendula, PSP. sylvestris, TCT. cordata, Oct_PS—October/P. sylvestris, Jul_PS—July/P. sylvestris, Oct_TC—October/T. cordata, Jul_BP—July/B. pendula, Jul_TC—July/T. cordata, Oct_BP—October/B. pendula

The highest species richness per sample was recorded in P. sylvestris stand in October 2022 (1.85 ± 0.35 species) and the lowest in P. sylvestris stand in July 2022 (0.30 ± 0.12) (Fig. 2B). Species richness of mesostigmatid mite communities was significantly lower in P. sylvestris stand in July (0.52 ± 0.15 species) than in P. sylvestris and T. cordata stands in October (1.37 ± 0.25 and 1.37 ± 0.19, respectively) (Fig. 3B).

The highest mesostigmatid mite Shannon’s diversity was found in P. sylvestris stand in October 2022 (0.52 ± 0.11), while the lowest was in T. cordata stand in July 2022 (0.05 ± 0.03). A significantly higher diversity was found in P. sylvestris stand in October (0.40 ± 0.07) compared to P. sylvestris and B. pendula stands in July (0.12 ± 0.05 and 0.14 ± 0.05, respectively) (Fig. 3C).

Mite communities structure

Most mite individuals belonged to the Parasitidae (111 ind.), Laelapidae (95) and Veigaiidae (70) families (Table 1). The most abundant species were Hypoaspis (Gaeolaelaps) aculeifer (21.6% of all recorded mites), Veigaia nemorensis (7.8%) and Trachytes aegrota (7.0%). Each of these species was most abundant in the T. cordata stand in October. Only two species occurred in each studied variant – Veigaia nemorensis and Hypoaspis aculeifer, while 15 mite species were unique for one variant (Fig. 4, Table 6 in Appendix 3). The cumulative species richness for T. cordata stand in October and July exceeded those for other examined month/habitat combinations (Fig. 5). The lowest juvenile abundance was recorded in T. cordata stand in July 2022 (0.07 ± 0.05 ind.; 2.4% of all juveniles). However, the number of juveniles was not statistically different among studied habitats and collection months (Fig. 3C). Interestingly, most females (2.19 ± 0.49 ind.; 22.2% of all females) were reported in P. sylvestris stand in October 2022, while the lowest female abundance was found in P. sylvestris stand in July of the same year (0.15 ± 0.09) (Fig. 6A). Males were the most numerous in the T. cordata stand in October 2022 (0.52 ± 0.20 ind.; 28.0% of all males). Moreover, no males were recorded in July 2022 in T. cordata and B. pendula stand (Fig. 4B). Furthermore, the most juveniles were found in P. sylvestris (0.56 ± 0.19 ind.; 18.1% of all juveniles) and T. cordata stand (0.56 ± 0.15; 18.1%) in October 2022, as well as in B. pendula stand in October 2021 (0.56 ± 0.30; 12.1%) (Fig. 6C).

Fig. 4
figure 4

Co-occurrence network for bipartite relationships between mesostigmatid mite taxa (upper boxes) and collection Month/Habitat (lower boxes). Boxes are proportional to total mite abundance, whereas ribbon width is proportional to the co-occurrence. Oct_PS—October/P. sylvestris, Jul_PS—July/P. sylvestris, Oct_TC—October/T. cordata, Jul_BP—July/B. pendula, Jul_TC—July/T. cordata, Oct_BP—October/B. pendula. For abbreviations of mite taxa see Table 6 in Appendix 3

Fig. 5
figure 5

Cumulative species richness (data presented as mean values) for each month/habitat examined. Oct_PS—October/P. sylvestris, Jul_PS—July/P. sylvestris, Oct_TC—October/T. cordata, Jul_BP—July/B. pendula, Jul_TC—July/T. cordata, Oct_BP—October/B. pendula

Fig. 6
figure 6

Mean female (A), male (B) and juvenile (C) abundance of mesostigmatid mite communities across four sampling periods (July, October; 2020, 2021) in pure stands on post-agricultural lands (Opole Forest District).. Oct—October, Jul—July, BPB. pendula, PSP. sylvestris, TCT. cordata

The largest number of H. aculeifer females was found in samples taken in October from the P. sylvestris and T. cordata stands (in both cases there were 22 individuals from two years, which accounted for 81.48% of all individuals of this species in the variant). In contrast, the largest number of H. aculeifer juveniles was found in samples taken in July from the T. cordata stand (6 ind.; 42.86%). Similarly, in the case of V.nemoresnis, the most females were found in samples collected in October from the T.cordata stand (10 ind.; 83.33%), while the most juvenile stages were reported in July in the T. cordata stand (4 ind.; 66,67%). No T. aegrota individuals were found in B. pendula stands, while the highest number of individuals of this species was reported in the T. cordata stands in October (35.71% of all T. aegrota individuals in total).

Discussion

Seasonal changes in mite communities

Generally, we found a significantly higher abundance of Mesostigmata mites in October than in July, as we expected from the results of other studies (Fujii and Takeda 2017). It alludes to the study by Kaczmarek et al. (2011), highlighting that soil moisture is the main factor that allows soil mite assemblage to develop. Indeed, in the climatic conditions of Poland, lower abundance in July may be explained by higher temperature and rapidly decreasing soil moisture in summer, which causes a decrease in mite density. With the end of summer, the temperature remains high, but soil moisture increases. Consequently, this causes an increase in soil mite density (Salmane 2000). The same relationship also applied to species richness and diversity of soil acarofauna. Similar importance of precipitation and temperature for the soil mite community structure was also revealed by Kamczyc et al. (2022), who reported that the abundance of mesostigmatid mites in decomposing litter of broadleaved and coniferous temperate tree species was positively correlated with the temperature of the sampling month and negatively correlated with the temperature of the month before sampling. The difference with our study was that the abundance shown in summer (July) was higher than in autumn (October). Species richness and diversity showed only slight changes between consecutive months of the growing season (April–October).

An interesting result of our study is the increase in the proportion of juvenile instars in summer. They are considered to be more sensitive than adult mite instars to temperature increases and humidity decreases, which is associated inter alia with a lower degree of body sclerotization (Huhta and Hänninen 2001). Adult mites should be more resistant than juveniles to water loss and less susceptible to changes in temperature and humidity (Urbanowski et al. 2021). An explanation for our result may be the overlap of the period of high mortality with the hatching of mesostigmatid juveniles. As demonstrated by Kaczmarek et al. (2011), under similar climatic conditions, there are two density peaks for mesostigmatid mites—in January and at the turn of spring and summer.

Most identified individuals belonged to the Parasitidae, Laelapidae and Veigaiidae families. Mites of the Parasitidae family are found in all types of forests, meadows, bogs, and farmland. They mainly inhabit the litter but also the humus layer. Their main prey are springtails and other mites, such as Oribatida and Astigmata (Błoszyk 2008). The Laelapidae family includes, among others, large predatory species of the Hypoaspis genus. Hypoaspis (Geolaelaps) aculeifer is a common pioneer species (Wissuwa et al. 2012), also found in large numbers on former farmland. Hypoaspis aculeifer is a species that inter alia preys on springtails (Cortet et al. 2003). Veigaiidae family contains species specialised in feeding on springtails, but may also feed on other soil fauna (Koehler 1999). They primarily inhabit forest litter in the Palearctic but are also found in meadows and agricultural land (Błoszyk 2008). Together with the mites of the Parasitidae family, they form the dominant group of predators in the Mesostigmata order (Koehler 1997; Błoszyk 2008; Skorupski 2008; Kazemi et al. 2013). Veigaia nemorensis is a common species in the litter and upper soil layers on former farmland and even on degraded or post-industrial sites (Skorupski et al. 2013; Manu et al. 2017). Veigaia nemorensis, in contrast to the general characteristics of the Veiigaidae family, mainly preys on bacterial and root-feeding nematode species (Manu et al. 2017).

Tree species impact

The differences in litter impact of tree species such as P. sylvestris, B. pendula and T. cordata on mesostigmatid mite assemblages were compared by Kamczyc et al. (2019) and Urbanowski et al. (2021). Kamczyc et al. (2019) showed that tree litter only affected abundance, while there was no influence on species richness and diversity of mesostigmatid mite assemblages. In addition, they indicated a positive effect of P. sylvestris litter comparing to T. cordata litter on soil mite density. In contrast, Urbanowski et al. (2021) showed no significant differences in the effect of P. sylvestris and B. pendula litter on the abundance, species richness and diversity of mesostigmatid mite communities. Our study did not show differences in the influence of individual tree species on the soil Mesostigmata. However, when comparing mite assemblages in summer and autumn, we observed high differences in mite abundance, species richness and diversity between the studied stands. Our results indicate that the most extreme abundance values ​​were found in P. sylvestris stand, which provided the most favourable conditions for the soil fauna in autumn and the least favourable in summer. This effect may be explained by the highest water absorption capacity of P. sylvestris litter, a coniferous species, with the greatest decreases in moisture during dry months (Zhou et al. 2018). Interestingly, Ilek et al. (2024) revealed that forest litter with a low C:N ratio and a predominant share of oak leaves achieved the greatest water storage capacity, while litter with a high C:N ratio and a predominant share of P. sylvestris needles had the lowest water storage capacity. This is consistent with our results because the C:N ratio was the lowest in the P. sylvestris stand and the highest in the T. cordata stand. Also in the B. pendula stand, mite abundance and diversity were significantly lower in summer than in P. sylvestris stand in autumn. We assume that this results from the fact that B. pendula stands are less able to protect the soil in summer from heating and drying out (Jonczak et al. 2020).

Interaction between soil parameters and mite communities

Our study revealed that Fe content in soil affected abundance, species richness and diversity of mite communities—the biodiversity decreased as soil Fe concentration increased. The negative effects of this nutrient on soil fauna were also noted by Huot et al. (2018) and Rusek and Marshall (2000). An increase in soil Fe concentration may limit the development of springtails, reduce manganese uptake, and stress soil organisms through the energy-demanding detoxification process (Rusek and Marshall 2000). Christophe et al. (2014) proved that the concentration of elements such as Fe in the rhizosphere solution and the soil environment changes with the seasons and is highest in autumn. Widowati and Sukristyonubowo (2019) showed that Fe becomes less toxic to plants when the soil pH is lower. In turn, Thoms and Gleixner (2013) observed that soil pH is lower in autumn than in summer.

Also, the collection month significantly affected the abundance, species richness and diversity of mesostigmatid mite communities. Changes in weather conditions with the ongoing seasons appear to be crucial to the density and diversity of soil fauna communities. This is not surprising because mites are ectothermic organisms, depending on climatic conditions such as temperature and precipitation (Thakur et al. 2018). This result was also confirmed by other studies conducted in central Europe (Urbanowski et al. 2021; Kamczyc et al. 2022). Additionally, species richness and diversity of mesostigmatid mites increased with litter nitrogen content, which may be associated with a higher biomass and species richness of microorganisms and nematodes (Kaneko and Salamanca 1999; Sánchez-Moreno et al. 2009; Renčo et al. 2020). Interestingly, in our study, the abundance of mites was not influenced by litter nitrogen content, which is correlated with the type of litter (Cornwell et al. 2008; Horodecki and Jagodziński 2019). It corresponds to the conclusions of studies conducted by Seniczak et al. (2018) on moss mite (Oribatida) and mesostigmatid mite communities. They reported that nutrient-rich habitats generate high species diversity but low densities of soil mites. This may mean that poorly fertile habitats have the opposite trend in mite community structure.

Conclusions

In conclusion, we revealed that different tree species growing on post-agricultural lands shape various soil conditions, which change between the summer and autumn seasons. Large fluctuations in soil mite abundance, species richness and diversity between P. sylvestris and B. pendula stands in summer and autumn may depend on tree species and season, which may constitute proof for using mixed stands during afforestation of post-agricultural lands. Seasonal change in the proportion between instars was unclear, but contrary to expected results, the proportion of young individuals increased in summer.