
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70311
ECE370311
ECE-2024-07-01393.R1
Applied Ecology
Ecosystem Ecology
Global Change Ecology
Soil Ecology
Research Article
Research Article
The influence of forest types including native and non‐native tree species on soil macrofauna depends on site conditions
Wenglein et al.
Wenglein Ronja https://orcid.org/0009-0006-7926-0877
1 ronja.wenglein@uni-goettingen.de

Lu Jing‐Zhong https://orcid.org/0000-0002-4051-8993
1
Scheu Stefan https://orcid.org/0000-0003-4350-9520
1 2
1 J.‐F. Blumenbach Institute of Zoology and Anthropology University of Göttingen Göttingen Germany
2 Centre of Biodiversity and Sustainable Land Use University of Göttingen Göttingen Germany
* Correspondence
Ronja Wenglein, J.‐F. Blumenbach Institute of Zoology and Anthropology, University of Göttingen, Untere Karspüle 2, Tierökologie, 37073, Göttingen, Lower Saxony, Germany.
Email: ronja.wenglein@uni-goettingen.de

18 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7031126 8 2024
09 7 2024
01 9 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

The ongoing climate change calls for managing forest ecosystems in temperate regions toward more drought‐resistant and climate‐resilient stands. Yet ecological consequences of management options such as planting non‐native tree species and mixing coniferous and deciduous tree species have been little studied, especially on soil animal communities, key in litter decomposition and pest control. Here, we investigated the taxonomic and trophic structure of soil macrofauna communities in five forest types including native European beech (Fagus sylvatica), range‐expanding Norway spruce (Picea abies) and non‐native Douglas fir (Pseudotsuga menziesii) as well as conifer‐beech mixtures across loamy and sandy sites in northern Germany. Abundance of primary decomposers (feeding predominantly on litter) was high in Douglas fir and beech forests, benefiting from less acidic soil and more favorable litter resources compared to spruce forests, while secondary decomposers (feeding predominantly on microorganisms and microbial residues) reached highest densities in spruce forests. Differences in abundance and species richness among forest types generally varied between regions and were most pronounced in Douglas fir of the sandy region. However, trophic guilds differed more between regions than between forest types, indicating that environmental factors outweigh the importance of forest type on soil macrofauna communities. The analysis of stable isotopes (δ15N and δ13C values) supported the general robustness in trophic position of macrofauna trophic guilds against variations in forest types and regions, but indicated reduced detrital shifts and food‐chain lengths in coniferous compared to European beech forests with mixtures mitigating these effects. Overall, for evaluating consequences of future forest management practices on the structure and functioning of soil animal communities, regional factors need to be considered, but in particular at loamy sites the taxonomic and trophic structure of soil macrofauna communities are resistant against changes in forest types.

Douglas fir
forest mixtures
macrofauna
non‐native species
soil fauna
stable isotopes
Deutsche Forschungsgemeinschaft 10.13039/501100001659 316045089 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:18.09.2024
Wenglein, R. , Lu, J.‐Z. , & Scheu, S. (2024). The influence of forest types including native and non‐native tree species on soil macrofauna depends on site conditions. Ecology and Evolution, 14 , e70311. 10.1002/ece3.70311
==== Body
pmc1 INTRODUCTION

Temperate forests are facing increasing decline due to climate extremes and associated bark beetle outbreaks calling for adaptation of silvicultural practices. European forests consist in large of monocultures. To enhance their resilience to increasing temperature and extreme weather events, forests are increasingly managed toward mixed forests including the plantation of non‐native tree species (Brang et al., 2014). These management practices likely are associated with changes in the diversity and functioning of above‐ and belowground animal communities. However, their consequences for the structure of soil animal communities are still little understood (Ganault et al., 2021). Especially consequences of the conversion of monocultures into mixed forests including non‐native tree species on the taxonomic and functional composition of soil animals needs further investigation. Soil animals are important drivers of litter decomposition, nutrient mineralization and soil formation (Abd El‐Wakeil, 2015; Barnes et al., 2018; Gilbert et al., 2014; Nielsen, 2019), and also serve as major agents linking belowground and aboveground components of ecosystems (Bardgett & van der Putten, 2014; Sackett et al., 2010; Scheu, 2001).

The dominant tree species of forests in central Europe, Norway spruce (Picea abies), is heavily suffering from drought and associated bark beetle outbreaks (Krejza et al., 2021). Confronted with the dieback of Norway spruce, the non‐native, more drought resistant Douglas fir (Pseudotsuga menziesii) is increasingly considered as alternative (Vitali et al., 2017). Furthermore, as Douglas fir is less acidifying the soil compared to Norway spruce and its leaf litter is more palatable for soil organism than that of Norway spruce, the humus form in Douglas fir forests is of higher quality than in Norway spruce forests which is also reflected by enhanced nitrification (Kupka et al., 2013; Podrázský et al., 2020). Nevertheless, the impact of Douglas fir on soils and soil organisms is still little known and depends on site conditions (Cremer & Prietzel, 2017; Lu & Scheu, 2021; Mrak et al., 2024). For example, microbial biomass has been shown to be significantly reduced in Doulgas fir at sandy but not all loamy sites compared to native European beech forests (Lu & Scheu, 2021). As is typical for introduced tree species, Doulas fir is colonized by only few, mostly generalist herbivore species (Schmid et al., 2014). Concerning soil animals, the abundance and activity of soil mesofauna has been found to be higher in Norway spruce than in Douglas fir forests (Kohlert & Roth, 2000; Rożen et al., 2010), whereas the abundance of macrofauna decomposers has been found to increase in Douglas fir compared to Norway spruce forests (Engel, 2001). However, little is known on soil macrofauna communities and their trophic organization in mixed stands including Douglas fir.

Since monocultures of conifers are increasingly being replaced by mixed stands of conifers and broadleaf species, in particular European beech (Fagus silvatica), changes due to mixing tree species are becoming increasingly relevant. It has been shown that the decomposition of litter in mixed spruce‐beech forests is faster than in pure spruce stands, independent of leaf/needle litter species (Berger & Berger, 2014). Mixed stands of conifers and beech are intermediate between the respective monocultures in storing soil organic carbon (Cremer et al., 2016). Furthermore, mixing conifers with beech reduces acidification (Foltran et al., 2023) and therefore may improve stand conditions. However, shifts in soil chemistry in mixed stands of European beech and Douglas fir are similar to those in mixed stands of European beech and Norway spruce (Foltran et al., 2020) but site conditions are playing an important role (Thurm & Pretzsch, 2016). Generally, increasing tree species richness has been found to increase the diversity of soil organisms (Ganault et al., 2021). However, abundance and functional diversity of ground beetles in mixed stands of beech and conifers were intermediate between the respective pure stands (Kriegel et al., 2021). Similarly, the biomass of different macrofauna groups in mixed conifer‐beech stands was intermediate between the respective pure stands (Scheu et al., 2003). More research on a wider spectrum of taxa and considering different site conditions is needed to explore if increased tree diversity and mixing of tree species increases soil fauna diversity.

Soil macrofauna species are an important component of forest ecosystems significantly contributing to leaf litter decomposition (Gongalsky, 2021; Slade & Riutta, 2012). Since macrofauna species are sensitive to variations in habitat conditions, it is important to understand their response to planting non‐native tree species such as Douglas fir and to the plantation of mixed forests (Eggleton et al., 2005; Wu & Wang, 2019). Although soil animal diversity is resilient against forest management practices, macrofauna detritivores have been shown to sensitively respond to changes in environmental conditions such as pH (Pollierer et al., 2021). Also, the abundance and diversity of litter dwelling macrofauna predators such as spiders have been shown to sensitively respond to the plantation of Douglas fir forests (Kriegel et al., 2021; Matevski & Schuldt, 2021; Matevski & Schuldt, 2023).

Changes in forest types lead to changes in habitat conditions and supply of resources likely affecting the trophic structure of soil animal communities. The analysis of natural variations in stable isotope ratios of carbon (13C/12C) and nitrogen (15N/14N) is increasingly used for characterizing trophic niches of species. The method allows insight into the trophic position of species (15N/14N ratios) as well as the use of basal resources (13C/12C ratios) (Potapov et al., 2019). Tree species affect soil organic matter turnover and are associated with different litter and root stable isotope signatures (Lorenz et al. 2020). Furthermore, litter quality may influence stable isotope signatures of soil fauna guilds as indicated by different 13C/12C ratios in coniferous compared to beech forests (Klarner et al., 2014). Also, the trophic position of generalist predators has been shown to be lower in Douglas fir than in beech forests (Wildermuth et al., 2023), but the influence of Douglas fir on the trophic structure of soil macrofauna communities has not been examined.

Here we analyze the effect of Douglas fir compared to Norway spruce in monoculture and in mixture with European beech on macrofauna communities and their tropic niches. Three guilds spanning a wide range of trophic levels including primary decomposers, secondary decomposers and predators are investigated. We aimed at investigation the response of abundance, species richness and guild structure of macrofauna to the plantation of non‐native Douglas fir monocultures and mixtures with European beech, in comparison to monocultures of Norway spruce and European beech as well as mixtures of European beech and Norway spruce. To investigate the generality of the response, we studied forest stands differing in site conditions, i.e. loamy and sandy sites. Shifts in trophic niches of macrofauna communities were studied using natural variations in stable isotope ratios of carbon and nitrogen to investigate differences between forest types and sites.

We hypothesize that (1) the abundance and species richness of macrofauna guilds decreases from beech to spruce to Douglas fir forests, with the differences being more pronounced in sandy compared to loamy sites, (2) in mixed forests of spruce and beech as well as Douglas fir and beech the abundance of macrofauna is intermediate between the respective monocultures but species richness is higher, (3) macrofauna communities of Douglas fir and Douglas fir/beech mixtures predominantly comprise generalist species and are similar to those in beech/spruce mixtures, and (4) the trophic position of macrofauna species (15N/14N ratios) is little affected by forest type, while the use of basal resources (13C/12C ratios) differs among stand types as well as between sandy and loamy sites.

2 MATERIAL AND METHODS

2.1 Study sites

We investigated a total of eight study sites each composed of five forest types (quintets) located in northern Germany grouped into four sandy and four loamy sites. The four loamy sites are characterized by more fertile partly podsolic Cambisols and Luvisols. The four sandy sites are characterized by Podzols on out‐washed sand. The mean annual precipitation at the four loamy sites is 821–1029 mm while at the four sandy sites it is 627–746 mm. More detailed information on site characteristics is given in Ammer et al. (2020) and Foltran et al. (2020). The five different forest types at each site included monocultures of European beech (F. silvatica), Douglas fir (P. menziesii) and Norway spruce (P. abies), as well as Douglas fir/European beech and Norway spruce/European beech mixtures. The average age of trees was 80 years (Glatthorn et al., 2023). The distance between forest type stands within sites ranged from 76 m to 4600 m, the sites were between 5 and 190 km apart. Plots of 50 m × 50 m were established in each stand. In pure stands, more than 90% of total basal area comprised the target tree species; in mixed stands, the basal area of beech ranged 33%–53% and of conifers 53–60%.

2.2 Animal extraction and identification

Samples were taken between November 2017 and January 2018. Two soil cores of 20 cm diameter were taken at a distance of approximately 10 m per plot using a metal cylinder and separated into litter layer and two soil layers, 0–5 and 5–10 cm soil depth. Temporal variations in soil fauna biomass are only moderate in temperate forests (Schaefer, 1990), likely due to litter and soil buffering climate extremes. Macrofauna was extracted from these layers by heat using a high‐gradient extractor for seven to 10 days (Kempson et al., 1963). Animals were collected in 50% diethylene glycol, filtered through 200 μm mesh and transferred into 70% ethanol. Macrofauna were sorted into broad taxonomic groups (Araneida, Chilopoda, Coleoptera, Coleoptera larvae, Diplopoda, Diptera larvae, Isopoda, Lumbricidae, Pseudoscorpionida and rest). Chilopoda, Coleoptera, Diplopoda and Isopoda were identified to species level if possible and Coleoptera larvae to family level using keys of Eason (1964), Klausnitzer (1978), Hopkin (1991), Lohse and Lucht (1999), Lompe (2002), Bährmann (2008) and Decker et al. (2024). Araneidae and Lumbricidae were not determined to species due to partial degrading during storage. The two samples per plot were pooled at each layer and species richness and abundance were calculated from data pooled across layers. Juveniles were only included into community analysis.

The animal taxa were assigned to trophic guilds based on litter corrected δ15N values (see Table S1) with corrections from literature for the taxa of Carabidae larvae, Geophilidae and Lithobidae (Bonato et al., 2021). Species without available δ15N values were ascribed to trophic guilds according to literature (Lawrence and Newton 1980; Arnett, 2002; Irmler et al., 2018; Schaefer, 2018).

2.3 Stable isotope analysis

Natural variations in bulk stable isotope values of 15N and 13C were used to quantify trophic niches of macrofauna taxa. If possible one individual per plot was measured for Aleocharinae larvae, Cantheridae larvae, Chilopoda, Diplopoda, Elateridae larvae, Isopoda and Staphylinidae larvae, with one individual per species kept as voucher. From the rest of Coleoptera larvae and Coleoptera adults, a maximum of three individuals per species was analyzed. Individuals were freeze dried overnight and then stored in a desiccator with silica gel. Samples >1 mg were ground and homogenized using mortar and pestle, and a subsample was taken for stable isotope analysis. The dry weight of all individuals used for stable isotope analysis was measured and used for biomass calculations. Samples were weighed into tin capsules, and stable isotope values were measured using a coupled system of an elemental analyzer (Flash 2000, Thermo Fisher Scientific, Massachusetts, USA) and a mass spectrometer (Delta V Advantage, Thermo Electron, Bremen, Germany). If individual samples were below 100 μg dry weight, a modified setup for small sample size was used (Langel & Dyckmans, 2014). Atmospheric nitrogen and Vienna PeeDee belemnite were used as primary standard, acetanilide (C8H9NO, Merck, Darmstadt) was used as internal standard.

Natural variations in stable isotope ratios of carbon and nitrogen (δX) were expressed as δX (‰) = (Rsample‐Rstandard)/Rstandard х 1000, with R being the ratio between the heavy and light isotopes (13C/12C or 15N/14N). Animal stable isotope values were calibrated to that of leaf litter at plot level (Klarner et al., 2014); stable isotope values of litter were taken from Lu et al. (2022). Litter corrected δ13C and δ15N values, denoted as Δ13C and Δ15N values, were used for statistical analysis.

2.4 Statistical analysis

All analyses were performed in R 4.4.0 (R Core Team, 2024). We used linear mixed models (LMMs) for macrofauna community and generalized linear mixed models (GLMMs) to analyze the effect of forest type (Douglas fir, Douglas fir/European beech, European beech, Norway spruce/European beech, Norway spruce), region (sandy, loamy) and their interaction on macrofauna guilds. Generalized linear mixed models with a ‘poisson’ distribution were selected as the best fit for the data. The eight sites were included as random effect. The response variables included abundance and richness of the respective guilds. Specimens from the three layers per soil core were summed as one community to focus on the factors of interest. To check for normality and homoscedasticity, model residuals were plotted.

To analyze the species‐based community structure of the different guilds, nonmetric multidimensional scaling (NMDS) and permutational multivariate analyses of variance (PERMANOVA) were used with Bray‐Curtis dissimilarities. Only taxa occurring in more than one plot were included into the analysis. Environmental variables (soil pH, soil carbon content, C/N ratio, water content, PLFA proportion of Gram+ and Gram− bacteria and fungi, litter mass; Lu & Scheu, 2021) that may influence macrofauna community structure were selected through stepwise forward selection and fitted into NMDS using the envfit function.

Isotopic metrics (Layman et al. 2007, Cucherousset & Villéger 2015) were used to analyze trophic structure of the macrofauna and guilds. Values were scaled between 0 and 1 to place species in a two‐dimensional space of Δ15N and Δ13C values. For guilds, to reach the minimum number of replicates (>2) for the analysis, loamy and sandy sites were combined, resulting in a total of 10 plots for secondary decomposers and predators and 9 plots for primary decomposers. For primary decomposers, beech at loamy sites was excluded from the analysis as not enough samples were available to calculate multidimensional metrics. Single dimensional metrics included minimum, maximum and range calculated as biomass weighted means of Δ15N and Δ13C values. Five biomass weighted multidimensional isotopic metrics (Isotopic Divergence (IDiv), Isotopic Dispersion (IDis), Isotpic Evenness (IEve), Isotopic Uniqueness (IUni) and Isotopic Richness (IRic)) were calculated. IDiv measures the distance between species and the center of the convex hull area. Values close to 0 indicate that extreme values are rare in the community, values close to 1 represent a community dominated by extreme values. IDis combines IDiv and the convex hull area. IDis equals 0 if all species have the same isotopic values and tends to 1 when contrasting stable isotope values are abundant and far from the center of gravity. IEve quantifies the species distribution in the stable isotope space. IEve is close to 1 when community values are evenly distributed in space, while IEve tends to 0 when species are packed in clusters. IUni evaluates the closeness of isotopic values across the community. IUni values tend to 0 when many species share the same isotopic position and IUni tends to 1 when species are unique in their isotopic values. IRic represents the level of trophic diversity of all species in the community based on isotopic niches scaled from 0 to 1. More detailed information on the isotopic metrics is given in Cucherousset and Villéger (2015). Linear models were fitted for metrics and guilds with forest type and region as factors. The post hoc HSD.test function was used for inspecting differences among forest types and regions. For LMM, the lmer and glmer function of the ‘lme4’ package was used (Bates et al., 2015), and for NMDS and PERMANOVA the metaNMDS and adonis2 function of the ‘vegan’ package was used (Oksanen et al. 2022). For pairwise comparisons, the glht function of the ‘multcomp’ package was used (Hothorn et al., 2008). Graphics were generated with ‘ggplot2’ (Wickham, 2016).

3 RESULTS

A total of 3003 individuals including 80 juveniles were analyzed. Primary decomposers comprised 183 individuals and 21 taxa, while secondary decomposers comprised 1909 individuals and 48 taxa and predators comprised 911 individuals and 56 taxa. Elateridae larvae was the most abundant taxon accounting for 25% of total individuals followed by Cantharidae larvae accounting for 19% of total individuals. Fourty‐three taxa were only found once and 21 taxa were only found twice.

Abundance of total macrofauna differed significantly between forest types but this varied between regions; at sandy sites the abundance was highest in Norway spruce forests and declined to European beech to Doulgas fir forests, with the respective mixed forests being intermediate, whereas at loamy sites the abundance was similar in each of the forest types (Figure 1a, Table 1). Overall, macrofauna species richness did not differ significantly between regions but among forest types; it was highest in beech and lower in Douglas fir and spruce forests, with mixtures being in between the pure stands (Figure 1b, Table 1).

FIGURE 1 Abundance (ind/m2) (a) and species richness (number of species per sample) (b) of total macrofauna in five forest types [Douglas fir (Dou), Douglas fir‐beech mixture (DouBe), beech (Be), spruce‐beech mixture (SpBe), spruce (Sp)] and two regions (sandy, loamy), means ± SE.

TABLE 1 Generalized linear mixed effects models with Poission distribution of the effect of forest type, region and their interaction on the abundance and species richness of total macrofauna and three macrofauna guilds (primary decomposers, secondary decomposers and predators).

Abundance		Total macrofauna	Primary decomposers	Secondary decomposers	Predator	
	df	X 2	p‐value	X 2	p‐value	X 2	p‐value	X 2	p‐value	
Intercept	1	1533.09	<.001	7.30	.007	410.74	<.001	68.99	<.001	
Forest type (F)	4	29.79	<.001	20.78	<.001	273.23	<.001	73.96	<.001	
Region (R)	1	5.80	.016	0.09	.768	5.52	.019	0.22	.635	
F x R	4	112.44	<.001	15.40	.004	77.22	<.001	47.04	<.001	
Species richness		Total macrofauna	Primary decomposer	Secondary decomposers	Predator	
	df	X 2	p‐value	X 2	p‐value	X 2	p‐value	X 2	p‐value	
Intercept	1	361.69	<.001	0.25	.618			12.76	<.001	
Forest type (F)	4	10.29	.036	6.05	.195	1.15	0.887	9.40	.052	
Region (R)	1	0.48	.488	0.09	.763	0.16	0.691	2.06	.152	
F x R	4	8.44	.077	3.45	.486	1.86	0.7614	16.30	.003	
Note: X 2‐ and p‐values of type llI Anovas; significant effects are given in bold.

Abbreviation: df, degrees of freedom.

3.1 Abundance and species richness of trophic guilds

Across forest types and regions, the abundance of the three trophic guilds differed significantly (ANOVA: F2,14 = 38.6, p < .001); it increased from primary decomposers (74 ± 14 ind/m2; mean ± SE) to predators (367 ± 47 ind/m2) to secondary decomposers (770 ± 77 ind/m2). However, in each of the trophic guilds the abundance varied with forest type, and this depended on region (significant forest type × region interaction; Table 1). In primary decomposers at the sandy sites, the density was highest in beech and Douglas fir‐beech mixtures and lowest in Douglas fir and spruce monocultures, whereas at the loamy sites the density was lowest in pure spruce and beech‐spruce mixtures and similar in the other forest types (Figure 2a). Contrasting the pattern in primary decomposers, the abundance of secondary decomposers was generally highest in pure spruce and spruce‐beech mixtures than in the other forest types, with this being more pronounced at the sandy than at the loamy sites (Figure 2b). Similar to primary and secondary decomposers, the abundance of predators was generally higher at the sandy than at the loamy sites (Figure 2c), and similar to secondary decomposers, the abundance was low in Douglas fir and beech‐Douglas fir mixtures. However, at the sandy sites the abundance of predators was similarly high in beech, spruce and beech spruce mixtures, whereas at the loamy sites it was at similar low level in spruce and beech spruce mixtures and higher in beech and beech‐Douglas fir mixtures.

FIGURE 2 Abundance (ind/m2) and species richness (number of species per sample) for three trophic guilds (primary decomposers, secondary decomposers, predators) in five forest types [Douglas fir (Dou), Douglas fir‐beech mixture (DouBe), beech (Be), spruce‐beech mixture (SpBe), spruce (Sp)] in two regions (sandy, loamy), means ± SE.

Species richness of the three trophic guilds generally differed much less between forest types and regions than abundance (Table 1, Figure 2d–f). The numbers were generally low with on average 1.5 ± 0.5 primary decomposer, 7.1 ± 2.0 secondary decomposer and 5.3 ± 1.3 predator species per sample. Only the richness of predators varied among forest types but this depended on region (significant forest type × region interaction; Table 1, Figure 2f). In beech, Douglas fir and Douglas fir‐beech mixture it was higher at the loamy than at the sandy sites, whereas in spruce and spruce‐beech mixtures it was lower at the loamy than at the sandy sites.

3.2 Community composition including trophic guilds

Macrofauna community composition differed significantly among forest types (PERMANOVA; R 2 = .173, p = .001) and between regions (PERMANOVA; R 2 = .105, p = .001), but the interaction was not significant. The first axis mainly sperated the two regions and the second axis mainly the forest types (Figure 3). Separation of the regions was closely associated with Gram+ bacteria being more abundant at sandy sites. Among forerst types spruce and beech differed most with Douglas fir being intermediate. Seperation of forest types was closely associated with soil pH and (in opposite direction) with soil carbon.

FIGURE 3 Nonmetric multidimensional scaling (NMDS) reflecting the relationship between macrofauna species in five forest types [Douglas fir (Dou), Douglas fir‐beech mixture (DouBe), beech (Be), spruce‐beech mixture (SpBe), spruce (Sp)] at two regions (sandy and loamy) and environmental factors including soil pH (pH), soil carbon (C), litter mass (Litter), Gram+ bacteria (Gram+) and Gram− bacteria (Gram‐). Ellipes represent 95% confidence ranges for Be (green), DouBe (orange), Dou (brown), SpBe (dark blue) and Sp (light blue).

Primary decomposer community composition differed significantly between regions (PERMANOVA; R 2 = .128, p = .001) but not between forest types (Figure 4a). The first axis separated the two regions and the second axis the forest types (although not significant). The communities at the sandy sites were associated with high proportion of Gram+ bacteria, whereas those at the loamy sites with high pH.

FIGURE 4 Nonmetric multidimensional scaling (NMDS) reflecting the relationship between (a) primary decomposers, (b) secondary decomposers and (c) predators in five forest types (Dou = Douglas fir, DouBe = Douglas fir‐ beech mixture, Be = beech, SpBe = spruce‐beech mixture and Sp = spruce). Environmental factors include soil pH (pH), soil carbon (C), litter mass (Litter), Gram+ bacteria (Gram+) and Gram− bacteria (Gram‐). Ellipes represent 95% confidence ranges for Be (green), DouBe (orange), Dou (brown), SpBe (dark blue) and Sp (light blue).

Secondary decomposer community composition significantly differed between both forest types and regions (PERMANOVA; R 2 = .213, p = .002 and R 2 = .08, p = .001, respectively). The communities at the sandy sites again were associated with high proportion of Gram+ bacteria and those at the loamy sites with high pH (Figure 4b). The second axis mainly separated Douglas fir‐beech from spruce‐beech forests with the former being associated with higher amounts of litter and the latter with higher pH and higher proportion of Gram− bacteria.

Similar to primary decomposers, predator community composition significantly differed between regions (PERMANOVA; R 2 = .088, p = .001) but not among forest types (Figure 4c). Communities at the sandy sites again correlated with high proportion of Gram+ bacteria, whereas those at the loamy sites with a number of environmental factors including the amount of litter and soil carbon concentration.

3.3 Stable isotope metrics

The Δ13C values of the total macrofauna communities were higher in beech forests than in conifer‐only forests, and this was due to similar shifts in both minimum and maximum values (F 4,34 = 5.80, p = .001; Figure 5). In sandy sites, the Δ13C values were higher for average position (F 1,34 = 5.77, p = .021). Minimum Δ15 N values of total macrofauna communities also differed between forest types with lowest values in beech and similarly high values in spruce and Douglas fir forests (F 4,34 = 3.88, p = .011; Figure 5). By contrast, the range between minimum and maximum Δ15N values was largest in beech and lowest in spruce and Douglas fir forests (F 4,34 = 4.0, p = .009). Multidimensional metrics did not differ significantly among forest types, but Isotopic Dispersion differed between regions (Table S2).

FIGURE 5 Average position, maximum, minimum, and range between minimum and maximum Δ13C (upper panel) and Δ15N (lower panel) values of total macrofauna in Douglas fir (Dou), Douglas fir‐beech (DouBe), beech (Be), spruce‐beech (SprBe) and spruce forests (Spr); means (circles), confidence intervals and individual measurements (small dots); means marked with different letters differ significantly (p < .05, Tukey's HSD test).

In primary decomposers, none of the Δ13C metrics (average position, maximum, minimum and range) differed significantly among forest types, but between regions with average position (F 1,4 = 7.06, p = .033), maximum (F 1,4 = 8.31, p = .024) and minimum (F 1,4 = 8.63, p = .022) being 2.0, 1.6 and 2.8‰ higher at the sandy sites, respectively (Figure 6, Table S2). Secondary decomposers showed a similar pattern in Δ13C metrics, but neither differed significantly among forest types nor between regions (Table S2, Figure S1 and Data S2). Across forest types and regions, the average position, maximum, minimum and range were 2.7 ± 0.9‰, 4.9 ± 0.5‰, 0.8 ± 1.3‰, and 4.2 ± 1.1‰. In predators, minimum (F 5,5 = 6.85, p = .045) and range (F 5,5 = 31.63, p = .0028) of Δ13C values differed significantly between forest types but not between regions (Table S2); minimum values were highest in beech‐spruce mixture and lowest in spruce, and the range was largest in beech and spruce forests and smaller in the other forest types (Figure 7).

FIGURE 6 Average position, maximum, minimum and range between minimum and maximum Δ13C values of primary decomposers in loamy and sandy sites; means (circles), confidence intervals and individual measurements (small dots); means marked with different letters differ significantly (p < .05, Tukey's HSD test).

FIGURE 7 Average position, maximum, minimum, and range between minimum and maximum Δ13C (upper panel) and Δ15N (lower panel) values of predators in Douglas fir (Dou), Douglas fir‐beech (DouBe), beech (Be), spruce‐beech (SprBe) and spruce forests (Spr); means (circles) confidence intervals and individual measurements (small dots); means marked with different letters differ significantly (p < .05, Tukey's HSD test).

Δ15N values of primary and secondary decomposers generally did not differ significantly among forest types and between regions (Table S3); in primary decomposers across forest types and regions the average position, maximum, minimum and range were 0.6 ± 0.7‰, 3.3 ± 1.0‰, −2.6 ± 1.5‰ and 5.9 ± 1.9‰, respective values for secondary decomposers were 2.9 ± 1.5‰, 6.8 ± 1.3‰, −2.3 ± 3.7‰ and 9.2 ± 3.6‰. In predators, minimum Δ15N values were highest in spruce and lowest in Douglas fir and beech, with conifer‐beech mixtures being in between (ANOVA; F 5,5 = 227.41, p = .0001; Figure 7). Across forest types and regions, the average position, maximum, minimum and range were 4.5 ± 0.8‰, 8.5 ± 2.9‰, 0.2 ± 1.1‰ and 8.4 ± 2.8‰ (Figure S2).

Multidimensional metrics of each of the three guilds neither differed significantly among forest types nor between regions (Table S3, Figure S3), with the exception of isotopic evenness in primary decomposers which was significant higher at sandy than loamy sites (F 1,4 = 11.09, p = .013).

4 DISCUSSION

We investigated the response of three macrofauna trophic guilds, primary decomposers, secondary decomposers and predators, to different forest types including pure and mixed stands as well as native and introduced conifer species. Overall, differences in macrofauna abundance and species richness in the different forest types were guild specific, with only primary decomposers benefitting from Douglas fir compared to spruce forests. Secondary decomposers and predators reached maximum abundance in spruce forest and spruce‐beech mixtures especially at the sandy sites. Species richness of primary and secondary decomposers generally differed little between forest types, but at sandy sites species richness of predators was lowest in spruce forests and beech‐spruce mixtures. Communities of all three guilds differed between regions, but only the total macrofauna and secondary decomposer communities were influenced by forest type. Stable isotope analysis indicated limited effects of forest type on the trophic niches of macrofauna species of each of the trophic guilds, but strongly influenced Δ13C values and food‐chain length of the total macrofauna community.

4.1 Soil fauna guilds in different forest types and regions

Region, i.e. loamy and sandy sites, differed significantly in the abundance of trophic guilds but little in their richness, although both showed similar patterns. Typically, differences in abundance among forest types were more pronounced in the sandy than loamy sites, in line with our first hypothesis. At sandy sites, the abundance of primary decomposers in beech and Douglas fir‐beech mixed forests considerably exceeded that in spruce, spruce‐beech mixed and Douglas fir forests. As beech litter is of low food quality this is unlikely due to factors related to litter quality (Jacob et al., 2010; Thomas & Prescott, 2000) and may point to the importance of root‐derived resources. In fact, at sandy sites fine root biomass in beech forests has been shown to be larger than in coniferous forests (Lwila et al., 2021), and fine root biomass and production typically increase with nutrient deficiency (Lwila et al., 2023). In contrast to primary decomposers, secondary decomposers peaked in spruce, and at sandy sites also in beech‐spruce mixed forests, which may be related to the thick litter layer in spruce and associated supply of microbial resources. Previously, it has been shown that saprotrophic beetles prefer spruce to Douglas fir in particular at sandy sites (Gossner et al., 2016). The very low abundance of secondary decomposers in Douglas fir and Douglas fir‐beech forests may also be related to reduced soil microbial biomass associated with reduced root‐derived resources (Lu & Scheu, 2021; Lwila et al., 2021).

Typically, the abundance and species richness of total macrofauna and each of the three macrofauna trophic guilds in mixed stands was intermediate between pure coniferous and pure beech stands, supporting our second hypothesis for abundance but not species richness. This is consistent with earlier studies suggesting that mixed forests of beech and conifers increase habitat complexity and resource availability compared to conifer monocultures for a wide range of taxa including plants, fungi, ground beetles and springtails (Budde et al., 2011; Korboulewsky et al., 2021; Kriegel et al., 2021; Likulunga et al., 2021). Generally, however, differences in abundance between forest stands and also regions were considerably more pronounced than those between richness, as also reported for earthworms and mesofauna (Korboulewsky et al., 2016).

Primary decomposer abundance and richness was higher in Douglas fir than spruce forests, supporting our first and second hypotheses. Presumably, primary decomposers profit from higher food quality and the higher pH of Douglas fir compared to spruce needles (Pontégnie et al., 2005). At the sandy sites the abundance and richness of primary decomposers were highest in beech and Douglas fir‐beech mixtures indicating that potential negative effects of Douglas fir are offset by the presence of beech. Similarly, Engel (2001) found Isopoda and Diplopoda to benefit from Douglas fir compared to spruce, and David et al. (2023) found species richness of macrofauna detritivores to be increased in mixed conifer forests, supporting the positive influence of Douglas fir and its mixture with beech for primary decomposers.

In contrast to primary decomposers, the abundance of secondary decomposers was generally highest in spruce forests and similarly low in beech and Douglas fir forests. As most of the species comprised secondary decomposers, total macrofauna abundance was similar to secondary decomposers. Secondary decomposers comprised predominantly beetle larvae, in particular larvae of Elateridae and Aleocharinae, predominantly living as detritivores, but also consuming plant roots or living as predators (Samoylova & Tiunov, 2017; Wolters, 1989). Based on stable isotope analysis, we classified them as secondary decomposers and in our study they accounted for 40% of the individuals of this trophic guild. They typically reach high density in soils of low pH (Kula, 2014) and are found in high numbers in spruce forests (Schaefer, 1990; Scheu et al., 2003). This resembles the high density of secondary decomposer mesofauna, such as Collembola and Oribatida, in acidic forests with thick organic layers (Korboulewsky et al., 2016; Maraun & Scheu, 2000). Aleocharinae rove beetles favor microhabitats with high amounts of deadwood and thick humus layers (Irmler & Gürlich, 2007). Generally, larvae and adults of rove beetles live as predators and reach higher abundance in spruce than beech forests (Sührig & Schaefer, 2001) preying e.g., on microarthropods such as springtails. However, in part they also feed on fungi and this applies in particular to Aleocharinae (Irmler & Lipkow, 2018; Scheu & Falca, 2000). Although their abundance differed between forest types, species richness remained unaffected by forest type, this, however, needs to be treated with caution as larvae were only identified to subfamily level.

Predator abundance and richness were generally highest in beech forests and at the sandy sites also in spruce and beech‐spruce mixed forests. Predatory macroarthropods predominantly comprised centipedes, which typically reach high density in temperate forests and dominate the biomass of macrofauna predators (Scheu et al., 2003). A number of centipede species reach high abundance in beech dominated forests feeding on a wide spectrum of prey (Bonato et al., 2021; Ferlian et al., 2012; Potapov et al., 2022). Adult and larvae of rove beetles also contributed to the high density of macrofauna predators. Both presumably feed heavily on secondary decomposer mesofauna such as Collembola (Günther et al., 2014; Hartmann, 1979), which reach high densities in spruce and beech forests (Kohlert & Roth, 2000; Salamon et al., 2008).

4.2 Community composition

Contrasting the strong difference in total macrofauna community composition among forest types, the community composition of the three macrofauna trophic guilds mainly differed between regions but little among forest types, partly rejecting our third hypothesis. Regional environmental factors have also been shown to be more important than forest type in shaping soil mesofauna communities (Erdmann et al., 2012). NMDS ordination identified soil pH and soil organic carbon as important factors structuring the community composition of total as well as each of the three trophic guilds of macrofauna, with the separation of forest types being most pronounced in total macrofauna communities, which is similar to earlier studies (Pollierer et al., 2021). Strong separation of beech and spruce forests along pH and soil organic carbon gradients reflects soil acidification as well as litter accumulation in spruce stands, with the intermediate position of Douglas fir likely reflecting lower adidification than in spruce forests and generally more palatable litter compared to spruce and beech (Kupka et al., 2013). In addition to abiotic factors, microbial community composition, in particular Gram+ bacteria, were identified as important structuring force of total as well as each of the three trophic guilds of macrofauna. Compared to Gram− bacteria, Gram+ bacteria more heavily depend on complex litter compounds (Fanin et al., 2019; Kramer & Gleixner, 2008) and their role in separating macrofauna trophic guilds between sandy and loamy sites may reflect that litter resources play a more important role in nourishing macrofauna communities at the loamy than at the sandy sites, similar to what has been observed in Collembola along a mountain gradient (Lux et al., 2024a, 2024b). The fact that the community structure of secondary decomposers also differed among forest types, mostly between Douglas fir‐beech and spruce‐beech mixtures, may also be related to differences in microbial community structure, in particular Gram+ bacteria being associated with Douglas fir‐beech mixed forests, likely again reflecting increased litter quality in Douglas fir‐beech mixed forests compared to beech forests, as also indicated by higher abundance of primary decomposers. However, also other factors such as canopy openness not considered in our study may contribute to differences in macrofauna community structure between forest types and this needs closer attention in future studies.

4.3 Trophic structure

Generally, the trophic structure of the three trophic guilds of macrofauna varied mostly between coniferous and beech/mixed forest with little difference among forest types and between regions in the three trophic guilds, pointing to the resistance of the trophic structure of soil food webs to variations in environmental factors. The consistency of the average trophic position of each of the three trophic guilds across forest types, as indicated by Δ15N values, supports our fourth hypothesis. Total macrofauna community had lower minimum Δ15N values in beech and Douglas fir‐beech mixed forests, likely resulting from comprising more primary decomposers with low Δ15N values. Presumably, this was also responsible for the wider range of Δ15N values in beech compared to conifer forests. Among macrofauna guilds, only minimum Δ15N values of predators were significantly higher in spruce forests than in the other forest types. Potentially, predators in spruce forest feed less on low trophic level prey such as primary decomposers than in other forest types (Ferlian & Scheu, 2014; Scheu & Falca, 2000), which likely is related to low food quality of spruce needles resulting in low abundance of primary decomposers, whereas secondary decomposers such as Collembola are thriving in spruce forests (Korboulewsky et al., 2021).

The detrital shift in δ13C values, which has been identified as a typical feature of soil food webs, also applied to the studied macrofauna trophic guilds and averaged 3.2‰, which is well in the range reported previously (Potapov et al., 2019). The lower shift in δ13C values in total macrofauna in coniferous compared to beech forests may reflect the low quality of beech litter (Klarner et al., 2014). For macrofauna guilds, however, only in predators minimum and range of Δ13C values differed among forest types, with in particular spruce and beech‐spruce mixed forests differing from the other forest types. Potentially, this again is related to the thick leaf litter layer in spruce and beech‐spruce mixed forests restricting the access to prey deeper in soil relying on root‐derived resources (Ferlian & Scheu, 2014; Günther et al., 2014). The lower average, maximum and minimum Δ13C values of total macrofauna and primary decomposers at the loamy than the sand sites may point to increased consumption of more nutrient rich and less decomposed litter resources at the former sites.

5 CONCLUSIONS

Overall, soil macrofauna trophic guilds and communities varied more between regions than among forest types, with differences between forest types being most prominent in total macrofauna communities. However, in particular primary decomposers benefited from non‐native Douglas fir, whereas the opposite was true for secondary decomposers. Presumably, this reflects higher food quality of Douglas fir compared to spruce leaf litter favoring primary decomposers, and poor litter quality of spruce litter resulting in a shift toward secondary decomposers feeding on microorganisms in spruce forests. The influence of mixed forests on soil macrofauna abundance and species richness typically was intermediate between the respective pure stands suggesting that mixed forests may buffer potential detrimental effects of pure stands of both native and non‐native conifers. However, in particular at sandy sites Douglas fir detrimentally affected macrofauna trophic guilds pointing to the importance of regional factors. The strong influence of regional factors was also reflected in the distinct difference in the structure of macrofauna trophic guilds between sandy and loamy sites markedly exceeding the influence of forest types. Differences in Δ13C values of total macrofauna communities in different forest types point to the adaptation of soil macrofauna food webs to differences in resource quality and environmental factors between deciduous and conifer forest. Lower minimum Δ15N values in beech than conifer forests support this conclusion. Significant changes in the trophic position of macrofauna predators with forest type point to predators as being most sensitive to environmental changes. Overall, the results indicate that mixed forests of beech and conifers may buffer potential detrimental effects of monocultures of conifers with this also applying to non‐native Douglas fir.

AUTHOR CONTRIBUTIONS

Ronja Wenglein: Conceptualization (equal); data curation (equal); formal analysis (lead); investigation (lead); methodology (equal); project administration (supporting); software (equal); validation (equal); visualization (lead); writing – original draft (lead); writing – review and editing (equal). Jing‐Zhong Lu: Conceptualization (equal); data curation (equal); formal analysis (supporting); investigation (equal); writing – review and editing (equal). Stefan Scheu: Conceptualization (supporting); data curation (supporting); formal analysis (equal); funding acquisition (lead); investigation (supporting); methodology (lead); project administration (lead); supervision (lead); visualization (supporting); writing – original draft (supporting); writing – review and editing (equal).

CONFLICT OF INTEREST STATEMENT

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Supporting information

Data S1.

Data S2.

ACKNOWLEDGEMENTS

We acknowledge Christian Bluhm and Theodora Volovei for field and lab assistance. We thank Susanne Boening‐Klein for sorting animals and her advice on stable isotope measurements, and also credits to the Center of Stable Isotope Research and Analysis of the University Göttingen (KOSI) for stable isotope measurements. Special thanks to S. Müller for coordinating the project. The work formed part of the Research Training Group 2300 (EnriCo) funded by the German Research Foundation (DFG; grant no 316045089). Open Access funding enabled and organized by Projekt DEAL.

DATA AVAILABILITY STATEMENT

Data and codes used in the study are available in the Supplementary Material.
==== Refs
REFERENCES

Abd El‐Wakeil, K. F. (2015). Effects of terrestrial isopods (Crustacea: Oniscidea) on leaf litter decomposition processes. The Journal of Basic & Applied Zoology, 69 , 10–16. 10.1016/j.jobaz.2015.05.002
Ammer, C. , Annighöfer, P. , Balkenhol, N. , Hertel, D. , Leuschner, C. , Polle, A. , Lamersdorf, N. , Scheu, S. , & Glatthorn, J. (2020). RTG 2300‐study design, location, topography and climatic conditions of research plots in 2020. PANGAEA. 10.1594/PANGAEA.923125
Arnett, R. H. (2002). American beetles. Vol. 2, Polyphaga: Scarabaeoidea through Curculionoidea. CRC Press.
Bardgett, R. D. , & van der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature, 515 , 505–511. 10.1038/nature13855 25428498
Barnes, A. D. , Jochum, M. , Lefcheck, J. S. , Eisenhauer, N. , Scherber, C. , O'Connor, M. I. , De Ruiter, P. , & Brose, U. (2018). Energy flux: The link between multitrophic biodiversity and ecosystem functioning. Trends in Ecology & Evolution, 33 , 186–197. 10.1016/j.tree.2017.12.007 29325921
Bates, D. , Kliegl, R. , Vasishth, S. , & Baayen, H. (2015). Parsimonious mixed models. Douglas Bates and Reinhold Kliegl and Shravan Vasishth and Harald Baayen.arXiv preprint arXiv: 1506.04967 10.48550/ARXIV.1506.04967
Berger, T. W. , & Berger, P. (2014). Does mixing of beech (Fagus sylvatica) and spruce (Picea abies) litter hasten decomposition? Plant and Soil, 377 , 217–234. 10.1007/s11104-013-2001-9 24744450
Bonato, L. , Peretti, E. , Sandionigi, A. , & Bortolin, F. (2021). The diet of major predators of forest soils: A first analysis on syntopic species of Chilopoda through DNA metabarcoding. Soil Biology and Biochemistry, 158 , 108264. 10.1016/j.soilbio.2021.108264
Brang, P. , Spathelf, P. , Larsen, J. B. , Bauhus, J. , Bonc Ina, A. , Chauvin, C. , Drossler, L. , Garcia‐Guemes, C. , Heiri, C. , Kerr, G. , Lexer, M. J. , Mason, B. , Mohren, F. , Muhlethaler, U. , Nocentini, S. , & Svoboda, M. (2014). Suitability of close‐to‐nature silviculture for adapting temperate European forests to climate change. Forestry, 87 , 492–503. 10.1093/forestry/cpu018
Budde, S. , Schmidt, W. , & Weckesser, M. (2011). Auswirkungen der Einbringung von Buche (Fagus sylvatica L.) auf die Artendiversität und Naturnähe von Nadelholzbeständen in Niedersachsen.
Bährmann, R. (2008). Bestimmung wirbelloser Tiere: Bildtafeln für zoologische Bestimmungsübungen und Exkursionen (5th ed.). Spektrum Akademischer Verlag.
Cremer, M. , Kern, N. V. , & Prietzel, J. (2016). Soil organic carbon and nitrogen stocks under pure and mixed stands of European beech, Douglas fir and Norway spruce. Forest Ecology and Management, 367 , 30–40. 10.1016/j.foreco.2016.02.020
Cremer, M. , & Prietzel, J. (2017). Soil acidity and exchangeable base cation stocks under pure and mixed stands of European beech, Douglas fir and Norway spruce. Plant and Soil, 415 , 393–405. 10.1007/s11104-017-3177-1
Cucherousset, J. , & Villéger, S. (2015). Quantifying the multiple facets of isotopic diversity: New metrics for stable isotope ecology. Ecological Indicators, 56 , 152–160. 10.1016/j.ecolind.2015.03.032
David, J.‐F. , Coq, S. , Decaens, T. , Ganault, P. , Hättenschwiler, S. , & Nahmani, J. (2023). Influence of tree species richness and identity on saprophagous macroarthropod abundance and species richness in some European forest types. Forest Ecology and Management, 534 , 120862. 10.1016/j.foreco.2023.120862
Decker, P. , Allspach, A. , Wesenberg, J. , & Xylander, W. E. (2024). BODENTIER hoch 4 – Onlineportal mit App zum Erleben, Erkennen, Erfassen und Erforschen. https://bodentierhochvier.de/
Eason, T. E. H. (1964). Centipedes of the British Isles. Warne.
Eggleton, P. , Vanbergen, A. J. , Jones, D. T. , Lambert, M. C. , Rockett, C. , Hammond, P. M. , Beccaloni, J. , Marriott, D. , Ross, E. , & Giusti, A. (2005). Assemblages of soil macrofauna across a Scottish land‐use intensification gradient: Influences of habitat quality, heterogeneity and area. Journal of Applied Ecology, 42 , 1153–1164. 10.1111/j.1365-2664.2005.01090.x
Engel, K. (2001). Zersetzer und ihre leistung: buche‐fichte‐douglasie. LWF, 33 , 18–21.
Erdmann, G. , Scheu, S. , & Maraun, M. (2012). Regional factors rather than forest type drive the community structure of soil living oribatid mites (Acari, Oribatida). Experimental & Applied Acarology, 57 , 157–169. 10.1007/s10493-012-9546-9 22460402
Fanin, N. , Kardol, P. , Farrell, M. , Nilsson, M.‐C. , Gundale, M. J. , & Wardle, D. A. (2019). The ratio of gram‐positive to gram‐negative bacterial PLFA markers as an indicator of carbon availability in organic soils. Soil Biology and Biochemistry, 128 , 111–114. 10.1016/j.soilbio.2018.10.010
Ferlian, O. , & Scheu, S. (2014). Shifts in trophic interactions with forest type in soil generalist predators as indicated by complementary analyses of fatty acids and stable isotopes. Oikos, 123 , 1182–1191. 10.1111/j.1600-0706.2013.00848.x
Ferlian, O. , Scheu, S. , & Pollierer, M. M. (2012). Trophic interactions in centipedes (Chilopoda, Myriapoda) as indicated by fatty acid patterns: Variations with life stage, forest age and season. Soil Biology and Biochemistry, 52 , 33–42. 10.1016/j.soilbio.2012.04.018
Foltran, E. C. , Ammer, C. , & Lamersdorf, N. (2020). Douglas fir and Norway spruce admixtures to beech forests along in northern Germany–Are soil nutrient conditions affected? (preprint). biorxiv. 10.1101/2020.09.25.313213
Foltran, E. C. , Ammer, C. , & Lamersdorf, N. (2023). Do admixed conifers change soil nutrient conditions of European beech stands? Soil Research, 61 , 647–662. 10.1071/SR22218
Ganault, P. , Nahmani, J. , Hättenschwiler, S. , Gillespie, L. M. , David, J.‐F. , Henneron, L. , Iorio, E. , Mazzia, C. , Muys, B. , Pasquet, A. , Prada‐Salcedo, L. D. , Wambsganss, J. , & Decaëns, T. (2021). Relative importance of tree species richness, tree functional type, and microenvironment for soil macrofauna communities in European forests. Oecologia, 196 , 455–468. 10.1007/s00442-021-04931-w 33959812
Gilbert, K. J. , Fahey, T. J. , Maerz, J. C. , Sherman, R. E. , Bohlen, P. , Dombroskie, J. J. , Groffman, P. M. , & Yavitt, J. B. (2014). Exploring carbon flow through the root channel in a temperate forest soil food web. Soil Biology and Biochemistry, 76 , 45–52. 10.1016/j.soilbio.2014.05.005
Glatthorn, J. , Appleby, S. , Balkenhol, N. , Kriegel, P. , Likulunga, L. E. , Lu, J. , Matevski, D. , Polle, A. , Riebl, H. , Rivera Pérez, C. A. , Scheu, S. , Seinsche, A. , Schall, P. , Schuldt, A. , Wingender, S. , & Ammer, C. (2023). Species diversity of forest floor biota in non‐native Douglas‐fir stands is similar to that of native stands. Ecosphere, 14 , e4609. 10.1002/ecs2.4609
Gongalsky, K. B. (2021). Soil macrofauna: Study problems and perspectives. Soil Biology and Biochemistry, 159 , 108281. 10.1016/j.soilbio.2021.108281
Gossner, M. M. , Wende, B. , Levick, S. , Schall, P. , Floren, A. , Linsenmair, K. E. , Steffan‐Dewenter, I. , Schulze, E.‐D. , & Weisser, W. W. (2016). Deadwood enrichment in European forests–Which tree species should be used to promote saproxylic beetle diversity? Biological Conservation, 201 , 92–102. 10.1016/j.biocon.2016.06.032
Günther, B. , Rall, B. C. , Ferlian, O. , Scheu, S. , & Eitzinger, B. (2014). Variations in prey consumption of centipede predators in forest soils as indicated by molecular gut content analysis. Oikos, 123 , 1192–1198. 10.1111/j.1600-0706.2013.00868.x
Hartmann, P. (1979). Biologisch‐Ökologische Untersuchungen an Staphylinidenpopulationen verschiedener Ökoysteme des Solling (Dissertation). Georg‐August‐Universität.
Hothorn, T. , Bretz, F. , & Westfall, P. (2008). Simultaneous inference in general parametric models. Biometrical Journal, 50 (3 ), 346–363.18481363
Hopkin, S. (1991). A key to the woodlice of Britain and Ireland, First edition, Environmental understanding for all. Field Studies Council.
Irmler, U. , & Gürlich, S. (2007). What do rove beetles (Coleoptera: Staphylinidae) indicate for site conditions. Faunistisch‐Ökologische Mitteilungen, 8 , 439–455.
Irmler, U. , Klimaszewski, J. , & Betz, O. (2018). Introduction to the biology of rove beetles. In O. Betz , U. Irmler , & J. Klimaszewski (Eds.), Biology of rove beetles (Staphylinidae) (pp. 1–4). Springer International Publishing. 10.1007/978-3-319-70257-5_1
Irmler, U. , & Lipkow, E. (2018). Effect of environmental conditions on distribution patterns of rove beetles. In O. Betz , U. Irmler , & J. Klimaszewski (Eds.), Biology of rove beetles (Staphylinidae) (pp. 117–144). Springer International Publishing. 10.1007/978-3-319-70257-5_7
Jacob, M. , Viedenz, K. , Polle, A. , & Thomas, F. M. (2010). Leaf litter decomposition in temperate deciduous forest stands with a decreasing fraction of beech (Fagus sylvatica). Oecologia, 164 , 1083–1094. 10.1007/s00442-010-1699-9 20596729
Kempson, D. , Monte, L. , & Ghelardi, R. (1963). A new extractor for woodland litter. Pedobiologia, 3 , 1–21.
Klausnitzer, B. (1978). Ordnung Coleoptera (LARVEN). Springer Netherlands. 10.1007/978-94-009-9975-6
Klarner, B. , Ehnes, R. B. , Erdmann, G. , Eitzinger, B. , Pollierer, M. M. , Maraun, M. , & Scheu, S. (2014). Trophic shift of soil animal species with forest type as indicated by stable isotope analysis. Oikos, 123 , 1173–1181. 10.1111/j.1600-0706.2013.00939.x
Kohlert, A. , & Roth, M. (2000). The effects of exotic tree species (douglas fir: Pseudotsuga menziesii) on litter decomposing arthropods on soil surface dwelling regulators. Mitteilungen der Deutschen Gesellschaft Fpr Allgemeine Und Angewandte Entomologie, 12 , 71–74.
Korboulewsky, N. , Heiniger, C. , De Danieli, S. , & Brun, J. J. (2021). Effect of tree mixture on Collembola diversity and community structure in temperate broadleaf and coniferous forests. Forest Ecology and Management, 482 , 118876. 10.1016/j.foreco.2020.118876
Korboulewsky, N. , Perez, G. , & Chauvat, M. (2016). How tree diversity affects soil fauna diversity: A review. Soil Biology and Biochemistry, 94 , 94–106. 10.1016/j.soilbio.2015.11.024
Kramer, C. , & Gleixner, G. (2008). Soil organic matter in soil depth profiles: Distinct carbon preferences of microbial groups during carbon transformation. Soil Biology and Biochemistry, 40 , 425–433. 10.1016/j.soilbio.2007.09.016
Krejza, J. , Cienciala, E. , Světlík, J. , Bellan, M. , Noyer, E. , Horáček, P. , Štěpánek, P. , & Marek, M. V. (2021). Evidence of climate‐induced stress of Norway spruce along elevation gradient preceding the current dieback in Central Europe. Trees, 35 , 103–119. 10.1007/s00468-020-02022-6
Kriegel, P. , Matevski, D. , & Schuldt, A. (2021). Monoculture and mixture‐planting of non‐native Douglas fir alters species composition, but promotes the diversity of ground beetles in a temperate forest system. Biodiversity and Conservation, 30 , 1479–1499. 10.1007/s10531-021-02155-1
Kula, E. (2014). Revitalization liming and the response of soil fauna. Acta Universitatis Agriculturae et Silviculturae Mendelianae Brunensis, 58 , 149–158. 10.11118/actaun201058040149
Kupka, I. , Podrázský, V. , & Kubeček, J. (2013). Soil‐forming effect of Douglas fir at lower altitudes ‐ a case study. Journal of Forest Science, 59 , 345–351. 10.17221/27/2013-JFS
Langel, R. , & Dyckmans, J. (2014). Combined 13C and 15N isotope analysis on small samples using a near‐conventional elemental analyzer/isotope ratio mass spectrometer setup. Rapid Communications in Mass Spectrometry, 28 , 1019–1022. 10.1002/rcm.6878 24677523
Lux, J. , Xie, Z. , Sun, X. , Wu, D. , Pollierer, M. M. , & Scheu, S. (2024a). Trophic positions of soil micro arthropods in forests increase with elevation, but energy channels remain unchanged. Ecosphere, 15 , e4747. 10.1002/ecs2.4747
Lux, J. , Xie, Z. , Sun, X. , Wu, D. , & Scheu, S. (2024b). Trophic niches of Collembola communities change with elevation, but also with body size and life form. Oecologia, 204 , 491–504. 10.1007/s00442-023-05506-7 38265599
Lorenz, M. , Derrien, D. , Zeller, B. , Udelhoven, T. , Werner, W. , & Thiele‐Bruhn, S. (2020). The linkage of 13C and 15N soil depth gradients with C:N and O:C stoichiometry reveals tree species effects on organic matter turnover in soil. Biogeochemistry, 151 , 203–220. 10.1007/s10533-020-00721-3
Lawrence, J. F. , & Newton, A. F., Jr. (1980). Coleoptera associated with the fruiting bodies of slime molds (Myxomycetes). The Coleopterists Bulletin, 34 , 129–143.
Layman, C. A. , Arrington, D. A. , Montaña, C. G. , & Post, D. M. (2007). Can stable isotope ratios provide for community‐wide measurements of trophic structure? Ecology, 88 , 42–48. 10.1890/0012-9658(2007)88[42:CSIRPF]2.0.CO;2 17489452
Likulunga, L. E. , Rivera Pérez, C. A. , Schneider, D. , Daniel, R. , & Polle, A. (2021). Tree species composition and soil properties in pure and mixed beech‐conifer stands drive soil fungal communities. Forest Ecology and Management, 502 , 119709. 10.1016/j.foreco.2021.119709
Lohse, G. A. , & Lucht, W. H. (1999). Die Käfer Mitteleuropas, Band 14, Die Käfer Mitteleuropas. Springer Berlin Heidelberg/Springer Spektrum.
Lompe, A. (2002). Die Käfer Europas, Ein Bestimmungswerk im Internet. Käfer Europas. https://www.coleonet.de/coleo
Lu, J. , Cordes, P. H. , Maraun, M. , & Scheu, S. (2022). High consistency of trophic niches in generalist arthropod species (Oribatida, Acari) across soil depth and forest type. Ecology and Evolution, 12 , e9572. 10.1002/ece3.9572 36523514
Lu, J.‐Z. , & Scheu, S. (2021). Response of soil microbial communities to mixed beech‐conifer forests varies with site conditions. Soil Biology and Biochemistry, 155 , 108155. 10.1016/j.soilbio.2021.108155
Lwila, A. S. , Mund, M. , Ammer, C. , & Glatthorn, J. (2021). Site conditions more than species identity drive fine root biomass, morphology and spatial distribution in temperate pure and mixed forests. Forest Ecology and Management, 499 , 119581. 10.1016/j.foreco.2021.119581
Lwila, A. S. , Post‐Leon, A. , Ammer, C. , & Mund, M. (2023). Site properties, species identity, and species mixture affect fine root production, mortality, and turnover rate in pure and mixed forests of European beech, Norway spruce, and Douglas‐fir. Ecological Indicators, 147 , 109975. 10.1016/j.ecolind.2023.109975
Maraun, M. , & Scheu, S. (2000). The structure of oribatid mite communities (Acari, Oribatida): Patterns, mechanisms and implications for future research. Ecography, 23 , 374–382. 10.1111/j.1600-0587.2000.tb00294.x
Matevski, D. , & Schuldt, A. (2021). Tree species richness, tree identity and non‐native tree proportion affect arboreal spider diversity, abundance and biomass. Forest Ecology and Management, 483 , 118775. 10.1016/j.foreco.2020.118775
Matevski, D. , & Schuldt, A. (2023). Non‐native Douglas fir promotes epigeal spider density, but has a mixed effect on functional diversity. Biodiversity and Conservation, 32 , 1233–1250. 10.1007/s10531-023-02547-5
Mrak, K. , Covre‐Foltran, E. , & Lamersdorf, N. (2024). Elevated nitrate concentrations in soil solution under pure Douglas fir stands can be lowered by mixing with European beech and by site selection. Forest Ecology and Management, 564 , 122004. 10.1016/j.foreco.2024.122004
Nielsen, U. N. (2019). Soil Fauna assemblages: Global to local scales (1st ed.). Cambridge University Press. 10.1017/9781108123518
Oksanen, J. , Simpson, G. , Blanchet, F. , Kindt, R. , Legendre, P. , Minchin, P. , O'Hara, R. , Solymos, P. , Stevens, M. , Szoecs, E. , Wagner, H. , Barbour, M. , Bedward, M. , Bolker, B. , Borcard, D. , Carvalho, G. , Chirico, M. , De Caceres, M. , Durand, S. , … Weedon, J. (2022). _vegan: Community ecology package_. R package version 2.6‐4.
Podrázský, V. , Kupka, I. , & Prknová, H. (2020). Substitution of Norway spruce for Douglas‐fir: Changes of soil microbial activities as climate change induced shift in species composition – A case study. Central European Forestry Journal, 66 , 71–77. 10.2478/forj-2020-0007
Pollierer, M. M. , Klarner, B. , Ott, D. , Digel, C. , Ehnes, R. B. , Eitzinger, B. , Erdmann, G. , Brose, U. , Maraun, M. , & Scheu, S. (2021). Diversity and functional structure of soil animal communities suggest soil animal food webs to be buffered against changes in forest land use. Oecologia, 196 , 195–209. 10.1007/s00442-021-04910-1 33852071
Pontégnie, M. , de Warnaffe, G. B. , & Lebrun, P. (2005). Impacts of silvicultural practices on the structure of hemi‐edaphic macrofauna community. Pedobiologia, 49 , 199–210. 10.1016/j.pedobi.2004.09.005
Potapov, A. M. , Beaulieu, F. , Birkhofer, K. , Bluhm, S. L. , Degtyarev, M. I. , 0044evetter, M. , Goncharov, A. A. , Gongalsky, K. B. , Klarner, B. , Korobushkin, D. I. , Liebke, D. F. , Maraun, M. , Mc Donnell, R. J. , Pollierer, M. M. , Schaefer, I. , Shrubovych, J. , Semenyuk, I. I. , Sendra, A. , Tuma, J. , … Scheu, S. (2022). Feeding habits and multifunctional classification of soil‐associated consumers from protists to vertebrates. Biological Reviews, 97 , 1057–1117. 10.1111/brv.12832 35060265
Posit team . (2024). RStudio: Integrated development environment for R. Posit Software, PBC. http://www.posit.co/
Potapov, A. M. , Tiunov, A. V. , & Scheu, S. (2019). Uncovering trophic positions and food resources of soil animals using bulk natural stable isotope composition. Biological Reviews, 94 , 37–59. 10.1111/brv.12434 29920907
Rożen, A. , Sobczyk, Ł. , Liszka, K. , & Weiner, J. (2010). Soil faunal activity as measured by the bait‐lamina test in monocultures of 14 tree species in the Siemianice common‐garden experiment, Poland. Applied Soil Ecology, 45 , 160–167. 10.1016/j.apsoil.2010.03.008
Sackett, T. E. , Classen, A. T. , & Sanders, N. J. (2010). Linking soil food web structure to above‐ and belowground ecosystem processes: A meta‐analysis. Oikos, 119 , 1984–1992. 10.1111/j.1600-0706.2010.18728.x
Salamon, J.‐A. , Scheu, S. , & Schaefer, M. (2008). The Collembola community of pure and mixed stands of beech (Fagus sylvatica) and spruce (Picea abies) of different age. Pedobiologia, 51 , 385–396. 10.1016/j.pedobi.2007.10.002
Samoylova, E. S. , & Tiunov, A. V. (2017). Flexible trophic position of polyphagous wireworms (Coleoptera, Elateridae): A stable isotope study in the steppe belt of Russia. Applied Soil Ecology, 121 , 74–81. 10.1016/j.apsoil.2017.09.026
Schaefer, M. (1990). The soil fauna of a beech forest on limestone: Trophic structure and energy budget. Oecologia, 82 , 128–136. 10.1007/BF00318544 28313148
Schaefer, M. (2018). Brohmer‐Fauna von Deutschland: Ein Bestimmungsbuch unserer heimischen Tierwelt, 25., durchgesehene Auflage. ed. Quelle & Meyer Verlag.
Scheu, S. (2001). Plants and generalist predators as links between the below‐ground and above‐ground system. Basic and Applied Ecology, 2 , 3–13. 10.1078/1439-1791-00031
Scheu, S. , Albers, D. , Alphei, J. , Buryn, R. , Klages, U. , Migge, S. , Platner, C. , & Salamon, J. (2003). The soil fauna community in pure and mixed stands of beech and spruce of different age: Trophic structure and structuring forces. Oikos, 101 , 225–238. 10.1034/j.1600-0706.2003.12131.x
Scheu, S. , & Falca, M. (2000). The soil food web of two beech forests (Fagus sylvatica) of contrasting humus type: Stable isotope analysis of a macro‐ and a mesofauna‐dominated community. Oecologia, 123 , 285–296. 10.1007/s004420051015 28308733
Schmid, M. , Pautasso, M. , & Holdenrieder, O. (2014). Ecological consequences of Douglas fir (Pseudotsuga menziesii) cultivation in Europe. European Journal of Forest Research, 133 , 13–29. 10.1007/s10342-013-0745-7
Slade, E. M. , & Riutta, T. (2012). Interacting effects of leaf litter species and macrofauna on decomposition in different litter environments. Basic and Applied Ecology, 13 , 423–431. 10.1016/j.baae.2012.06.008
Sührig, A. , & Schaefer, M. (2001). Diversity of epigeic arthropod communities in mixed and pure stands of beech and spruce II: Carabidae and Staphylinidae. Verhandlungen der Gesellschaft Fur Okologie, Verhandlungen der Gesellschaft Fur Okologie, 31 , 159.
Thomas, K. D. , & Prescott, C. E. (2000). Nitrogen availability in forest floors of three tree species on the same site: The role of litter quality. Canadian Journal of Forest Research, 30 , 1698–1706. 10.1139/x00-101
Thurm, E. A. , & Pretzsch, H. (2016). Improved productivity and modified tree morphology of mixed versus pure stands of European beech (Fagus sylvatica) and Douglas‐fir (Pseudotsuga menziesii) with increasing precipitation and age. Annals of Forest Science, 73 , 1047–1061. 10.1007/s13595-016-0588-8
Vitali, V. , Büntgen, U. , & Bauhus, J. (2017). Silver fir and Douglas fir are more tolerant to extreme droughts than Norway spruce in south‐western Germany. Global Change Biology, 23 , 5108–5119. 10.1111/gcb.13774 28556403
Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis (2nd ed.. 2016. ed). Use R! Springer International Publishing: Imprint: Springer.
Wildermuth, B. , Fardiansah, R. , Matevski, D. , Lu, J.‐Z. , Kriegel, P. , Scheu, S. , & Schuldt, A. (2023). Conifers and non‐native tree species shift trophic niches of generalist arthropod predators in central European beech forests. BMC Ecology and Evolution, 23 , 3. 10.1186/s12862-023-02105-1 36737705
Wolters, V. (1989). The influence of omnivorous elaterid larvae on the microbial carbon cycle in different forest soils. Oecologia, 80 , 405–413. 10.1007/BF00379044 28312070
Wu, P. , & Wang, C. (2019). Differences in spatiotemporal dynamics between soil macrofauna and mesofauna communities in forest ecosystems: The significance for soil fauna diversity monitoring. Geoderma, 337 , 266–272. 10.1016/j.geoderma.2018.09.031
