In uence of Selected Environmental Parameters on Rove Beetle (Coleoptera: Staphylinidae) Communities in Central European Floodplain Forests

Slavomír Stašiov Technical University of Zvolen: Technicka Univerzita vo Zvolene Juraj Litavský (  litavskyjuraj@gmail.com ) Comenius University in Bratislava Faculty of Natural Sciences: Univerzita Komenskeho v Bratislave Prirodovedecka fakulta https://orcid.org/0000-0003-1229-1098 Oto Majzlan Comenius University in Bratislava Faculty of Natural Sciences: Univerzita Komenskeho v Bratislave Prirodovedecka fakulta Marek Svitok Technical University of Zvolen: Technicka Univerzita vo Zvolene Peter Fedor Comenius University in Bratislava Faculty of Natural Sciences: Univerzita Komenskeho v Bratislave Prirodovedecka fakulta


Introduction
The alluvial landscape is characterized by a diverse composition of habitats. Its structure consists of a mosaic of wetlands, waterways, oxbow lakes, meadows and various types of oodplain forests. Their dynamics is predominantly affected by water regime (Kontriš 1981 undisputedly dynamic complexes with a high degree of biodiversity (Hughes 2007). For this reason, they are considered as priority biotopes in nature conservation.
Management and conservation of European wetlands requires a sound understanding of both extraordinary species-environment and species-species relationships in the communities, guilds etc. (Fedor 2004;Derunkov 2007Derunkov , 2009). Moreover, insects (e.g. King  Up to now 1,783 species of Staphylinoidea have been reliably known in Slovakian fauna (Zahradník 2017) while in Serbia 883 species have been recorded until now (Stančić 2013; Löbl and Löbl 2015; Majzlan and Litavský 2017). The Staphylinidae, or rove beetles, belong to the largest and most biologically diverse beetle families. The world fauna consists of more than 63,495 described species ). Representatives of the family are found among euedaphic species in deep soil layers, on the soil surface, on the vegetation up to the canopy of forests, and from marine habitats in the eulittoral zone of the sea up to high mountain zones (Thayer 2005). Moreover, the variety of feeding habits is astonishing (Irmler and Lipkow 2018). Most rove beetles are predatory as adults and larvae, feeding on other, smaller arthropods. Within the family, however, you'll nd rove beetles that specialize on a diet of fungal spores, others that eat pollen, and still others that feed on the regurgitated food from ants (Hadley 2020). Because of good knowledge of ecological demands of most Central European species of rove beetles and their distribution in all semi-natural and man-made habitats, many authors consider these organisms to be very suitable bioindicators of anthropogenic environmental changes (Boháč 1990(Boháč , 1999 Irmler and Lipkow 2018). Their research revealed the in uence of several factors on the abundance and structure of rove beetle communities in this type of environment (e.g. moisture, temperature and pH of soil, canopy cover, amount of dead biomass on soil surface, forest stand area and its circuit, habitat type, vegetation composition, fragmentation of landscape, ood characteristics, anthropogenic impact and others). However, the results of some studies contradict each other, and the signi cance of the in uence of individual factors on the formation of rove beetle communities, which was found by various authors, is also different. Therefore, we decided to test the impact of some controversial factors, as well as hitherto unexplored environmental parameters on the structure of rove beetle communities in oodplain forests.
Parallel analyses of staphylinid taxocoenoses in two similar ecosystems in different regions have inspired us to evaluate the impact of a wide range of environmental characteristics on their structure. Therefore, we compared the similar types of habitats in alluvia of three rivers, Danube (Slovakia), Tisa and Begej (Serbia), using 24 different parameters ( Table 2) and their in uence on staphylinid communities.

Study area
The research was carried out at eight sites situated in the alluvium of the three rivers: the Danube in Slovakia (S1-S5) and the Tisa (S6, S7) and the Begej in Serbia (S8). S1 -a fragment of oodplain forest Salici-Populetum. S4 -a shrub-meadow ecotone, located in the oodplain forest Biskupické luhy.
S5 -a hardwood oodplain forest located next to the oxbow lake Biskupické rameno.
S6 -a softwood oodplain forest located in the proximity of the oxbow lake of the Tisa River in the inundation area.
S8 -a hardwood oodplain forest located within the SNR "Carska bara" in the alluvium of the Begej River.
Additional details on the study area are provided in Litavský et al. (2021b).

Sampling
The research was performed from February 2015 until November 2016. The rove beetles were captured by pitfall trapping (Stašiov 2015), using 0.5 l plastic cups (9 cm in opening diameter) with and 1% formaldehyde as a xation uid (approximately one third of the volume). At each study site, ve traps were installed in a line in the distances of 5 m between neighbouring cups. These transects were placed in the middle of monitored sites. Traps were emptied and re lled approximately in half-month intervals. The contents of ve traps at each site and date were put together to represent a composite sample. The material was sorted in the laboratory and the rove beetles subsequently identi ed to the species level according to Lohse (1964Lohse ( , 1989, Assing and Schülke (2012). Specimens were xed in 75% ethyl alcohol and deposited at the Faculty of Natural Sciences of the Comenius University in Bratislava. A list of species is organized in accordance with (Newton 2019).
Concerning environmental properties, we decided to monitor the characteristics that could in uence the composition of staphylinid communities. Results about these attributes that we recorded, such as stand canopy of individual layers, species richness of individual layers, anthropogenic impact (four levels), fragment's size in which we placed pitfall traps, distance of located traps from another different fragment Due to relatively low number of sampling sites and high number of environmental variables, we did not build parsimonious models describing diversity and composition patterns, but we adopted an exploratory approach and investigated all species-environment relationships.
In order to investigate patterns of summary community characteristics (total dynamic activity, species richness, Shannon diversity and evenness), we calculated the matrix of pairwise Spearman correlation coe cients between environmental variables and community characteristics.
Principal coordinate analysis (PCoA) based on Bray-Curtis distance was used to visualize similarities in composition of rove beetle communities. The vectors of environmental variables were projected into ordinations in the directions of their maximal correlations with the con guration of sites. Statistical signi cance of tted environmental vectors was assessed using permutation tests (10,000 permutations).
This exploratory approach allowed us to investigate all species-environment relationships without the risk of missing any important information caused by removal of highly correlated variables. Indeed, the results of such a heuristic analysis should be considered carefully and regarded as a process of generation of new hypotheses rather than a generalization to wider populations.
The analyses were performed in R (R Core Team 2016) using libraries Hmisc (Harrell 2016) and vegan (Oksanen et al. 2016 recorded at a maximum of 5 sites. Apart from the Drusilla (Drusilla) canaliculata, Tasgius (Rayacheila) melanarius and Philonthus (Philonthus) spinipes, the most frequented species were recorded only at the Slovak sites. Species Omalium caesum, Oxyporus maxillosus (only in Slovakia) and Staphylinus caesareus (in Slovakia and in Serbia) were recorded at four sites. Sixty-eight (56 from Slovakian and 14 from Serbian sites) of rove beetles were recorded on one site only for the entire period of research, of which 37 were represented by one individual only. The largest number of individuals refers to the study site S1 (265 ind.) and site S5 (252 ind.), the highest species richness appeared at S5 (40 species) and S3 (39 species). The least individuals as well as species were caught on all Serbian sites.
The highest value of the Shannon's index of rove beetle communities was calculated for the site S2 (3.05) and the highest value of equitability index refers to the site S8 (0.80). The lowest values of Shannon index as well as equitability index of rove beetle communities were recorded at site S1 (2.34/0.30) ( Table  1).
The relationships between community characteristics and environmental variables were evaluated using non-parametric correlation analysis. Total dynamic activity of staphylinids showed signi cant positive relationship with number of plant species in E 2 vegetation layer ( Table 2). Species richness was signi cantly positively correlated with number of plant species in E 2 and soil pH, and negatively with relative H content of soil. We did not nd any signi cant correlation between Shannon diversity, but evenness was negatively linked with number of plant species in E 2 .

Discussion
The results of our research showed that the wetland habitats of the Danube, Begej and Tisa Rivers could provide suitable conditions for rove beetle communities with high species richness. Also, according to Kolesnikova et al. (2016), rove beetles are characterized by high species diversity in wetlands and alluvial forest. According to Derunkov (2009), wetland ecosystems in general, especially oodplain forests are key habitats for many rare species, including postglacial relicts. According to this author, the staphylinid community composition in oodplain forests is de ned by the eurytopic forest species. This opinion is also supported by the dominant occurrence of Ocypus mus and O. olens recorded at studied sites (38.8% from all caught rove beetles), because both species are typical for the moderately humid or dry forests (Šustek 1992;Massolo et al. 2006).
Similarly, high biodiversity of rove beetle communities in wetlands has been con rmed by several other authors. For instance, Derunkov (2007)  Among other factors, the rove beetles are also related to the quality of leaf litter and soil. An important property of soil is the content of N and P. Our research revealed that the largest number of individuals were recorded at the site S1 (265 ind.) and site S5 (252 ind.) (Table 1). They were characterized by a relatively high content of these nutrients in soil (N (S1) -0.29% weighted, P (S1) -36.0 mg.kg -1 ; N (S5) -0.29 % weighted, P (S5) -25.8 mg.kg -1 ) in comparison with other studied sites (within the overall range N: 0.15-0.36% weighted, P: 5.2-84.3 mg.kg -1 ). The content of these nutrients in soil and in the litter probably also had an effect on equitability of rove beetle communities, because the highest values equitability index was also recorded at site S8 (0.80) ( Table 1) with the highest P content in the soil (84.3 mg.kg -1 ) and also in the litter layer (407 mg.kg -1 ; within the overall range 85-407 mg.kg -1 ) and with the second highest N content in the soil (0.31% weighted) within the studied sites. The content of N and P in soil may have an indirect impact (through food resources of rove beetles -e.g. fungi, saprophagous) on the food available for rove beetles. Soil humus rich in nitrogen and phosphorus is an attractive source for fungi and saprophagous organisms. Therefore, such material is decomposed by activity of decomposers and saprophagous species more quickly (Wittich 1942(Wittich , 1943. The importance of nitrogen and phosphorus as the main elements determining the animal production and food resources for invertebrates, including rove beetles, has been reported by Dunger (1958). During the same research, we con rmed at the same study sites that the total dynamic activity of harvestmen showed signi cant positive relationship with proportion of N in the leaf litter, while the community equitability was negatively related to the H content of the soil (Litavský et al. 2018). As for the weevils, total dynamic activity of weevils showed a signi cant positive relationship with relative content of H in soil, while the species richness was negatively related to relative content of H and N in soil (Litavský et al. 2021b).
An important property of the soil is also its pH. We found that pH was signi cantly positively correlated with species richness of rove beetle communities (Table 2). Soil pH probably had a positive effect on species richness of rove beetle communities, because the highest species richness appeared at the sites S5 (40 species) and S3 (39 species) (Table 1), on which the highest (S3 -7.7) and the second highest pH of the soil (S5 -7.6) (within the overall range 5.9-7.7) was recorded. The pH of soil and the leaf litter probably also had a positive effect on the Shannon's index of rove beetle communities, because its highest value was calculated for the site S2 (3.05) that was characterized by the highest pH of the leaf litter (6.4; within the overall range 5.3-6.4) and the second highest pH of the soil (7.6). Moreover, according to Irmler and Lipkow (2018), the acidity of substrate also plays a role in the distribution of staphylinid species.
In terms of vegetation properties, the number of plant species in shrub vegetation layer had the most signi cant effect on rove beetle communities. This property was signi cantly positively correlated with the total dynamic activity and species richness of rove beetle communities and negatively linked with evenness of these communities ( Table 2). Rove beetle communities were probably also affected by the number of plant species in herb and tree layer, as well as stand canopy cover of shrub vegetation layer.
This assumption suggests the nding that the highest values of the Shannon's index of rove beetle communities (3.05) ( Table 1) was calculated for the site S2 with the highest number of plant species in tree vegetation layer (6 spp; within the overall range 0-6 spp). On the other hand, the second lowest number of plant species in herb vegetation layer (18 spp; within the overall range 16-32 spp) was found at the already mentioned sites S1 and S5 with the largest number of recorded individuals of staphylinids.
In addition, the site S1 was characterised also by lowest stand canopy of shrub vegetation layer (7%; within the overall range 7-65%). So, plant species richness of E 1 had the opposite effect on the rove beetle communities than plant species richness of E 2 and E 3 . It remains questionable whether the high total dynamic activity of rove beetle communities at the S1 site was signi cantly affected by the low number of plant species in E 1 or the low canopy cover of E 2 . The in uence of vegetation structure on rove beetle communities was also pointed out by other authors. For example, Klimaszewski et al. (2018) claim that broadleaved and mixed forests provide for staphylinids richer and moister ground litter, which is a more suitable habitat for these beetles. According to Rose (2001) and Klimaszewski et al. (2018), canopy cover is the most important factor in uencing the structure of rove beetle communities. During the same investigation, we found out that the total dynamic activity of weevils showed a signi cant positive relationship with number of plant species in shrub vegetation layer, while the species richness of weevils was signi cantly positively related to the number of plant species in herb and negatively to stand canopy of tree vegetation layers (Litavský et al. 2021b). We also found that, the harvestmen community equitability was positively related with species richness of plant communities in shrub layer (Litavský et al. 2018).
The obtained results also indicate the in uence of other evaluated factors on rove beetle communities. For example, thickness of litter layer probably had a positive effect on the total dynamic activity as well as the species richness of the staphylinids. This is indicated by the fact that the already mentioned pair of sites S1 and S5, with the largest number of recorded individuals (Table 1) also had the thickest (S5 -4.2 cm) and the second thickest litter layer (S1 -3.5 cm) within all studied sites (within the overall range 1.0-4.2 cm). At the site S5, with the thickest litter layer, we recorded most species of rove beetles ( Table  1). The in uence of the litter layer on rove beetle communities was also pointed out by other authors. According to Klimaszewski et al. (2018), moist litter is a very important habitat for most rove beetles, because their small bodies are prone to desiccation. Removal of the litter layer and exposure of mineral soil likely reduce soil moisture and habitat suitability for staphylinids. Irmler and Lipkow (2018) pointed that litter type is the most important factor determining the distribution of rove beetles. According to these authors, several staphylinid species responded positively to high amounts of litter fall. The mentioned highest value of the Shannon's index of rove beetle communities, recorded at the site S2 may also point to the effect of stand age on these communities. In addition to the highest number of plant species in the tree vegetation layer, this site was also characterized by the oldest vegetation in tree layer (70 years; within the overall range 10-70 years). According to Pohl et al. (2008), many staphylinid species require continuous, mature, or old growth stands. Another important factor that probably shaped the structure of studied rove beetle communities was anthropogenic impact. Namely, the sites S3 and S5 with recorded highest species richness (Table 1), in addition to the relatively high pH values of the soil, were also characterised by the least anthropogenic impact (degree 1) within all studied sites (within the overall range degrees 1-4). The in uence of anthropogenic impact on rove beetle communities was pointed out by several authors. For example, Pohl et al. (2008) pointed out that their species assemblages effectively re ect the extent of natural or human impact upon forest ecosystems. According to Klimaszewski et al. (2018), in many cases rove beetle species can be characterized as forest specialists with strong a nity for intact forest stands that have not been recently disturbed by forest management. Marcelino et al. (2016) studied the distribution and genetic variability of rove beetles in anthropogenically in uenced insular landscapes. They found that rove beetle richness was associated with anthropogenic in uence and habitat type, increasing from less to more anthropogenic impacted habitats.
Our research has revealed that wetland ecosystems in general, especially oodplain forests provide suitable conditions for staphylinid taxocoenoses with high species diversity. This diversity is conditioned by the fact that most species prefer moist habitats (Babenko et al. 2015;Irmler and Lipkow 2018;Klimaszewski et al. 2018) and also by diverse structure of microhabitats that oodplain forests provide to the rove beetles. The conservation and protection of these habitats is therefore important not only in terms of preserving rove beetle communities, but also in terms of maintaining high biodiversity and ecological stability in the landscape.

Conclusions
In summary, we examined how rove beetle communities vary across different habitats of river wetlands. During our investigation in the different stands of oodplain forests in Slovakia and Serbia, we found out that the total dynamic activity increased with species richness of shrub vegetation layer and the number of rove beetles was positively linked with the number of plant species in E 2 and soil pH. We also found that Ocypus (Ocypus) olens and Ocypus (Pseudocypus) mus preferred stands with higher soil pH and that these species could be used as suitable bioindicators to assess changes in landscape structure caused by human activity resulting in soil acidi cation. Therefore, more information about the importance of these relationships might be helpful in further elucidation of how rove beetle community metrics respond functionally to soil, vegetation, and microclimatic conditions, and how these conditions vary across various types of wetland habitats.   Tables   Table 1 Total number of rove beetle specimens and the diversity measures of rove beetle communities recorded at the study sites S1-S8   table 1