SUPPORTING INFORMATION Effects of management on aquatic tree-hole communities in temperate forests are mediated by detritus amount and water chemistry Martin M. Gossner, Peggy Lade, Anja Rohland, Nora Sichardt, Tiemo Kahl, Jürgen Bauhus, Wolfgang W. Weisser, Jana S. Petermann Table of Contents Appendix S1: Details on study plots and occurring tree-holes……………………………. 2 Appendix S2: Details on tree-hole mapping, sampling, and species identification………. 7 Appendix S3: Details on data analysis……………………………………………………. 13 Appendix S4: Details on water properties of sampled tree-holes………………………. 17 Appendix S5: Details on sampled dendrolimnotoxenes…………………………………... 18 Appendix S6: Community composition in tree-holes…………………………………… 20 Appendix S7: Detailed results of Structural Equation Models……………………………. 21 Appendix S8: Tree-hole richness………………………………………………………….. 26 Appendix S9: Diversity partitioning………………………………………………………. 27 Appendix S10: Detailed Acknowledgements……………………………………………... 28 References…………………………………………………………………………………. 29 1 SUPPORTING INFORMATION Appendix S1: Details on study plots and occurring tree-holes In the Schwäbische Alb (SA) 29 and in the Hainich-Dün (HD) 24 experimental plots (100 m x 100 m) were studied and represent a range of forest management intensities, from unmanaged beech forests to differently managed beech- and conifer forests (Fig. S1-1, Table S1-1). Unmanaged forests were set-aside 20-70 years ago, but all were influenced by humans at some point, since no virgin forests exist in Central Europe. In the Schwäbische Alb, the “unmanaged” forests (five plots) were formerly used as pasture woodland (approx. until the beginning of the 20th century) resulting in very old beech trees (200 years) with large crowns and low branching. The structural importance of these trees for biodiversity led to the protection of the stands, but small interventions such as removal of spruce trees to protect the old beech trees are still practiced in these forests. In the HD, unmanaged forests (six stands) are located within the National Park Hainich which covers an area of 7,600 ha and was founded in 1997. The core area of the national park was a protection zone of a military training range from 1945 until 1990 and is characterized by uneven-aged broad-leaved forests with tree ages of up to 250 years. The managed forests in our study represent different silvicultural systems: 1) in selection cutting forests only individual trees or small groups of trees are harvested according to their target diameter and stem quality, resulting in forests with an uneven age-structure (only in HD; 6 plots); 2) even-aged forests, where trees are from one cohort, are typically regenerated either using small clear cuts or strip cuts (coniferous forests) or through extended regeneration phases using shelterwood systems. Within age-class forests different developmental stages were considered (see Table S1-1): immature beech stands (thickets and stands of pole-sized trees; SA: 9 plots, HD: 5 plots), early-mature and mature beech stands (SA: 9, HD: 7), early-mature and mature spruce stands (Schwäbische Alb: 6). Previous 2 SUPPORTING INFORMATION studies have shown that plant diversity increases with forest management intensity in these forests due to increased disturbance (Boch et al. 2013). Fig. S1-1: Map showing the location of studied plots in the two regions Schwäbische Alb and Hainich-Dün. Different forest management types are indicated by different symbols. Grey shaded areas in the two regions represent forests, white areas predominantly represent arable land, but also urban areas and water bodies. The mean distance between the two regions is 321km (MIN: 297km, MAX: 344km) and among plots within region it is 7.3km (MIN: 0.3km; MAX: 19km) in the SA and 16.4km (MIN: 0.3km; MAX: 36km) in Hainich Dün. 3 SUPPORTING INFORMATION Table S1-1: The number of tree-holes occurring in the study plots (1ha) of the Schwäbische Alb and Hainch-Dün region, separated by stratum (canopy >2m height, ground <2m height), type (pan/rot hole) and formation type. AC=Age-class forest Canopy Pan Study plot Schwäbische Alb AEW01 AEW02 AEW03 AEW04 AEW05 AEW06 AEW07 AEW08 AEW09 AEW13 AEW14 AEW15 AEW16 AEW17 AEW22 AEW23 AEW26 AEW27 Ground Pan Rot hole Total Rot hole Main tree species Management type Forked Forked Branch Stem Forked Forked Root stem branch break hole stem branch hole Branch Stem Stump break hole Spruce Spruce Spruce Beech Beech Beech Beech Beech Beech Spruce Spruce Beech Beech Beech Beech Beech Beech Beech Age class, immature timber Age class, immature timber Age class, immature timber Age class, pole wood Age class, mature timber Age class, immature timber Unmanaged, mature timber Unmanaged, mature timber Unmanaged, mature timber Age class, mature timber Age class, mature timber Age class, thicket Age class, thicket Age class, pole wood Age class, immature timber Age class, mature timber Age class, pole wood Age class, pole wood 0 0 0 0 11 20 5 6 9 0 0 0 0 0 10 24 0 0 0 0 0 0 0 0 1 2 7 0 0 0 0 0 0 0 0 0 5 0 4 0 3 0 2 5 1 0 0 0 0 0 2 7 0 0 0 0 0 0 0 0 15 11 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 16 4 14 16 3 0 0 2 2 0 4 12 0 0 0 0 0 0 0 0 4 6 2 0 0 0 0 0 0 0 0 0 0 13 3 0 8 0 0 0 2 12 8 0 2 0 0 0 0 0 5 13 7 0 38 24 43 50 31 12 8 2 4 0 16 43 0 0 4 SUPPORTING INFORMATION AEW30 AEW32 AEW37 AEW38 AEW42 AEW43 AEW45 AEW47 AEW48 AEW49 AEW50 Beech Spruce Beech Beech Beech Beech Beech Beech Beech Beech Beech Age class, mixed immature timber Age class, immature timber Age class, thicket Age class, pole wood Age class, immature timber Age class, mature timber Age class, thicket Age class, mixed immature timber Age class, mixed immature timber Unmanaged, mature timber Unmanaged, mature timber 12 0 0 0 11 8 0 10 9 9 13 0 0 0 0 8 1 0 0 0 4 6 0 0 0 0 0 0 0 0 0 7 11 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 11 2 0 0 0 0 0 7 8 0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 16 0 13 8 19 9 0 10 9 28 38 HainichDün HEW04 HEW05 HEW06 HEW08 HEW10 HEW15 HEW16 HEW17 HEW18 HEW19 HEW21 HEW22 HEW26 HEW28 HEW32 HEW34 HEW37 Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Beech Age class, thicket Age class, immature timber Age class, mature timber Selection cutting Unmanaged, mature timber Age class, thicket Age class, pole wood Age class, pole wood Age class, pole wood Age class, immature timber Age class, mature timber Age class, mature timber Selection cutting Selection cutting Selection cutting Unmanaged, mature timber Unmanaged, mature timber 3 25 13 1 9 4 2 3 12 12 11 8 1 0 3 10 15 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 2 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 9 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 2 0 0 0 0 0 4 1 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 0 0 0 0 0 0 0 0 1 6 6 2 5 6 0 1 0 2 8 29 5 5 1 14 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 1 2 1 0 0 1 0 0 0 3 1 1 0 1 0 0 1 0 0 0 1 0 1 0 0 0 0 0 0 2 1 0 0 0 5 34 20 4 16 16 3 11 13 17 21 39 8 7 4 27 40 5 SUPPORTING INFORMATION HEW38 HEW39 HEW46 HEW47 HEW48 HEW49 HEW50 Total Beech Beech Beech Beech Beech Beech Beech Unmanaged, mature timber Unmanaged, mature timber Age class, immature timber Age class, mature timber Selection cutting Selection cutting Unmanaged, mature timber 14 14 10 7 6 2 6 4 1 5 2 0 0 1 8 14 4 0 0 0 4 1 9 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 18 12 5 5 7 6 2 0 0 0 0 0 0 0 0 6 1 0 0 0 4 0 0 0 0 0 0 0 45 56 25 14 13 8 18 348 65 91 16 10 10 262 12 37 59 910 6 SUPPORTING INFORMATION Appendix S2: Details on tree-hole mapping, sampling, and species identification Mapping of water-filled tree-holes For assessing the availability of canopy tree-holes we first used binoculars to assess the occurrence of potential water bodies in forked stems and branches, branch breaks and stem holes (see Fig. S2-1). In total, we mapped 910 tree-holes (Table S1-1). We then used a 15 meter long telescopic pole with a camera attached at the top (tree-top camera) for detailed assessment, in particular to check if the potential tree-hole was indeed water-filled and hence a potential habitat for a local community (Fig. S2-2). This allowed the assessment of treeholes at heights of up to 17m when adding the height of the observer. Our approach thus did not allow us to examine tree-holes higher than 17m in detail. Based on assessment with binoculars we estimated that only 3% of the tree-holes were at heights above 17m. Since the probability of desiccation increases with height because tree-holes located higher in the canopy are usually smaller and more exposed to sunlight and wind (Yanoviak 1999) the large majority of organisms can be reasonably expected to dwell in lower tree-holes. We thus consider our sampling to be representative of the tree-hole communities in these plots. Using the tree-top camera, each tree-hole was evaluated on a display attached to the telescopic pole and a stick attached to the top was used to measure the depth of tree-holes. The shape and surface area of each tree-hole was estimated and the maximum water volume was assessed using equations for volumes of cones, hemispheres or prisms depending on the shape of the tree-hole. For those tree-holes that were mapped but whose communities were not sampled, the experience from sampled tree-holes (see below) was used to arrive at a reliable estimation of water volume. Only tree-holes with a maximum water volume of >50ml were used to calculate the density of tree-holes (per 1-ha plot). Additionally, the tree-hole type (pan, maintaining an unbroken bark lining vs. rot hole, penetrating through to the wood of the tree; Kitching 1971), tree-hole origin (forked stem, forked branch, root hole, branch 7 SUPPORTING INFORMATION break, stem hole, stump), and the height above ground was noted. Water-filled holes did in rare cases occur in deadwood logs and these were assigned to the category “stem hole”. Treeholes were classified as “canopy tree-holes” when found higher than 2m above ground and “ground tree-holes” when found below 2m above the ground, following classifications by Yanoviak and Fincke (2005) and Kitching (1971). Fig. S2-1: Selection of tree-holes sampled in the present study. Stem forks are a common type found in managed beech forests (Fagus sylvatica) (a-c). Sometimes birds (i.e. here the common blackbird Turdus merula) build their nests above this type of tree-hole (b). Other common types of tree-holes in the studied forests were branch forks (pans) in beech (d), root holes (pans) in beech (e, k) and hornbean (Carpinus betulus) (j), and stumps (rot holes) in spruce (Picea abies; h). Rot holes, such as branch breaks (e.g. in beech; f), stem holes (e.g. in lime trees, Tilia spp.; g) and tree-hole in lying dead wood (i.e. in lying hornbean trunk; i), were much rarer, particularly in managed forests. 8 SUPPORTING INFORMATION Fig. S2-2: Mapping and sampling of tree-holes. a-c) A tree-top camera was used for mapping tree-holes in the canopy. d) A photo of each tree-hole seen with the tree-top camera was taken. e/i) Sampling of tree-holes in the canopy of (e) a managed age-class forest and (i) an unmanaged former woodland pasture in the SA (i) was done using single rope climbing technique. f) Sampling of a tree-hole at the base of a beech tree trunk (root hole). g/h) Treeholes in a branch fork with numerous larvae of the syrphid Myatropa florea with their long posterior breathing tubes erected to the water surface. 9 SUPPORTING INFORMATION Sampling of tree-hole communities The sampling of selected tree-hole communities was conducted from June 23 - July 20 in the SA and from June 10 -20, 2011 in the HD. Tree-hole communities at ground level were sampled in a total of 14 plots in the SA and 24 plots in the HD across the whole range of forest management intensities (Table S1-1). Due to time restrictions, canopy tree-holes were sampled only in 9 and 18 plots, respectively. In each plot we randomly selected tree-holes with a volume of 50ml (if not available smaller tree-holes were sampled), pans at ground level (seven in SA and two at HD) and pans in the canopy (five in SA and one in HD). Additionally, rot holes were sampled where available. In one plot of the HD and twelve plots of the SA, fewer than the targeted number of tree-holes occurred and thus only these were sampled. In two plots of the SA and one plot of the HD, tree-holes located higher than the mapping limit of 17m were included because of an insufficient number of tree-holes at lower heights. All sampled canopy tree-holes were reached by using the single-rope climbing technique (Perry 1978). The ropes were installed using a crossbow. The sampling of these tree-holes followed the protocol suggested by Yanoviak and Fincke (2005). The tree-hole was photographed and the surrounding bark was checked for insect larvae and pupae. Subsequently, pH, temperature, and oxygen saturation were measured by a mobile sensor (see paragraph “Assessment of water properties”) as these measures need to be taken prior to disturbing the system. Additionally, a 15-ml water sample was taken from the free water body with the help of a 10-ml syringe. The water was filtered using a 0.45-μm filter extension (PVDF diaphragm) and the filtrate was collected in a 15-ml falcon tube, immediately transferred to a cooler and frozen at -20 °C in the evening of each day. Subsequently, all water and detritus of the tree-hole were quantitatively removed and transferred to sampling jars. The tree-holes were repeatedly flushed with a laboratory wash bottle filled with rain water and subsequently checked for insect using a flashlight to ensure 10 SUPPORTING INFORMATION complete sampling. Then, tree parameters (diameter at breast height, DBH) and parameters of the tree-hole, such as exposition and height on the tree, type (pan/rot hole), origin (forked stem, forked branch, root hole, branch break, stem hole, stump), shape, depth, length and width were noted. The latter measurements were used to calculate the opening area the maximum volume of non-sampled tree-holes. The maximum volume was measured by filling the tree-hole with rain water. Finally, some leaves from around the tree were added to facilitate the recovery of the system. In the laboratory, the detritus of each water-filled tree-hole was separated into coarse detritus (leaves, flower parts etc.) and fine detritus using a sieve of 0.5mm mesh size and dried at room temperature. Assessment of water properties Oxygen saturation, pH, and water temperature were measured in the field with mobile electrodes (Oxi 330, pH 330, WTW GmbH). Nitrate- (NO3-), ammonium- (NH4+), and phosphate (PO34-) contents (mg/l) were measured by photometric determination of concentrations. For the determination of nitrate content, a subsample of 1 ml of each sample was diluted, if necessary, and mixed with 50 μl of a 5-% Resorcin solution and 1.3 ml of concentrated acid sulfur. After mixing (vortex) and heat dissipation, an absorbance measurement was taken in a 1.5 ml cuvette at 360 nm. The colorimetrical determination of ammonium 1ml of each water sample was charged with 0.5 ml natrium-phenolate-solution, 0.25 ml natrium-nitroprussid- and natrium-hypochlorite-NaOH-solution. The test tubes were subsequently stored in a water quench in the dark at 30°C for 30 minutes. Then, the absorbance measurement was taken in a 1.5-ml cuvette at 630 nm. For the measurement of phosphate concentration 2.5 ml of each water sample were charged with 100 μl Molybdat- 11 SUPPORTING INFORMATION reagent solution and 25 μl ascorbic acid –solution and mixed subsequently. After a settlement period of 20 minutes an absorbance measurement was taken in a 1.5-ml cuvette at 880 nm. Due to technical constraints, temperature, oxygen content, and phosphate could only be assessed in the HD. Keys used for insect identification Insect larvae and pupae were identified to family level following common insect larval identification keys (Stehr 1987; Nilsson 1996; Stresemann et al. 2005). For genus and species identification in Diptera we used for Ceratopogonidae Cranston (1982), Lindegaard (1997), and Pankratova (1970) for Chironomidae, Nilsson (1997), for Culicidae Cranston et al. (1987), Mohrig (1969), and Utrio (1976), for Muscidae Rozkošný and Gregor (2004), for Syrphidae Dixon (1960), Hartley (1961), and Rotheray (1993), for Psychodidae Nilsson (1997), and for Tipulidae Nilsson (1997). For species identification in the Coleoptera family Scirtidae the key of Klausnitzer (1996) was used. 12 SUPPORTING INFORMATION Appendix S3: Details on data analysis Calculation of tree-hole diversity Tree-hole density and the number of tree-hole types per plot were assessed using all 910 mapped tree-holes. We calculated the number of tree-hole types as follows. This method has previously been used to calculate the richness of dead wood structures (Siitonen et al. 2000). We included all tree-hole characteristics that were assessed for the mapped tree-holes (see above); tree species, exposition, direction of opening, height, type, type of origin, shape, opening area, depth, maximum water volume. The number of tree-hole parameters was reduced based on ecological meaningfulness and then by removing highly correlated parameter (Pearson’s correlation coefficient >0.4). Our final set of parameters for the variable tree-hole diversity comprised tree-hole origin, tree species, direction of opening (top, top-side, side), height (≤ 0.1 m, < 10 m, ≥ 10), and area of opening (≤ 20 cm², < 50 cm², ≥ 50 cm²). Analyses of causal relationships with Structural Equation Models To analyze causal relationships between forest management intensity, regional and local factors and tree-hole communities we used structural equation modelling (SEM). Variables and causal chains in our initial model were based on evidence from previous tree-hole studies (e.g. Kitching 1971; Schmidl, Sulzer & Kitching 2008). Here, we used a model-generating approach, where one tentative initial model is simplified until the model fit cannot be improved (Grace 2006). Analyses were conducted with the package lavaan in R (Rosseel 2014) which fits models using Maximum Likelihood based on covariance matrices. The hypothesized full models are shown in Figure S3-1. In the structural equation models, we used the region (SA, HD) as a grouping factor and merged chemical and physical water properties into a single composite variable, determined by pH, ammonium and nitrate concentration (mg/l). Because other chemical properties were only available from the HD 13 SUPPORTING INFORMATION region we excluded these from the models. We additionally calculated a separate model for the HD region including all available water properties. The results were qualitatively similar and our overall conclusions remained unchanged. Thus, we present only the model including both regions. All variables were tested for normality and homoscedasticity prior to analyses and transformed where necessary. For ammonium, detritus amount, and insect abundance and species richness an ln-transformation and for nitrate a square-root transformation was applied and the density of tree-holes per plot, the height of the tree-hole above ground, and maximum water volume were ln-transformed. We used two axes of the NMDS (see above) as response variable in our SEM (Tylianakis et al. 2008; Blakely & Didham 2010). For a measure of model fit in SEMs, the overall model p-value was used following Grace (2006, p.128f). For models without an appropriate fit (model p-value <0.05), we first used the modification indices to evaluate if important paths were missing from the model (MI values larger than 3.84), included those and subsequently deleted non-significant variables in a stepwise procedure, starting with those showing the highest p-values in both regions. We followed this procedure until an adequate model fit was achieved. To calculate the total effect (including all indirect effects) we multiplied path coefficients of each individual path and summed up all paths between forest management intensity and species abundance and richness, respectively. 14 SUPPORTING INFORMATION Figure S3-1: Hypothesized full model for the indirect effects of forest management on abundance and species richness (top) and community composition, measured as scores along the first and second axis of a Nonmetric Multi Dimensional Scaling (NMDS) based on BrayCurtis distances (bottom). The composite variable “water properties” was not measured directly. 15 SUPPORTING INFORMATION Appendix S4: Details on water properties of sampled tree-holes Table S4-1: Overview over the measured properties of the sampled tree-holes. Certain properties (O2 concentration, temperature, phosphate) were only measured in the HainichDün. Schwäbische Alb Hainich-Dün Property N Mean SE Min Max N Mean SE Min Max O2 [%] 0 NA NA NA NA 56 10.9 0.915 2.00 36.0 Temperature [°C] 0 NA NA NA NA 64 15.7 0.502 7.00 25.0 pH* 68 7.08 6.69/7.49 4.93 8.85 66 6.70 6.30/7.08 4.10 8.50 Ammonium [mg/l] 96 1.27 0.190 0.00 15.0 50 4.68 1.24 0.00 57.0 Nitrate [mg/l] 96 0.29 0.073 0.00 5.05 50 22.6 1.95 0.68 72.0 Phosphate [mg/l] 0 NA NA NA NA 47 1.79 0.22 0.00 6.34 Max. water volume [ml] 112 982 188 6.00 16620 87 1094 184 5.00 10192 Coarse detritus [ml] 112 155 41.4 0.00 4030 87 427 73.7 0.00 3530 Fine detritus [ml] 112 125 28.1 0.00 1570 87 142 35.4 0.00 1892 Total detritus [ml] 112 280 63.2 0.00 5600 87 566 98.0 0.00 4982 Proportion detritus** 112 0.343 0.025 0.00 1.67 87 0.580 0.078 0.00 5.00 Height [m]*** 112 0.554 0.03 34.0 87 3.89 0.593 0.10 18.3 4.37 * because pH is a logarithmic value median and quartiles are given ** proportion of detritus relative to maximum water volume *** height above ground 16 SUPPORTING INFORMATION Appendix S5: Details on sampled dendrolimnotoxenes A total of 41,146 dendrolimnetoxenes were found in the studied tree-holes. These species might not establish persistent populations there. They consisted mainly of oribatid mites (please note that abundance of mites was not measured in the Schwäbische Alb). Dendrolimnetoxenes comprised species from at least eight different orders (Tab. S5-1). Araneae were observed occasionally in water-filled tree-holes near ground level, but might have escaped during sampling. Similarly, Dermaptera were observed occasionally, but only up to a height of five meters. Acari, Collembola, Diplopoda, Gastropoda, Isopoda, and Oligochaeta (Lumbricidae, Enchytraeidae) frequently occurred in tree-holes and with the exception of Enchytraeidae (only at ground level) up to 15 meters above ground. Dendrolimnotoxenes were excluded from all analyses. Table S5-1: Taxonomic groups of dendrolimnetoxenes found in the tree-holes of the present study. Dendrolimnetoxenes are species that are present in tree-holes, but might not establish persistent populations there according to the definition given by Kitching (1971; 2000). All adult Lumbricidae were identified to species level. In the SA 112 and in the Hainich-Dün 87 tree-hole communities were sampled. The frequency gives the proportion of tree-holes in which a particular species occurred. NA=not assessed. Order Family Species Individuals (total) Schw. Alb Hainich Frequency [%] Schw. Alb Hainich Araneae 0 2 0 2.30 Collembola 5 15 4.46 9.20 Dermaptera 2 0 0.89 0 Diplopoda 16 4 11.61 1.15 Gastropoda 7 137 2.68 40.23 31 34 13.39 10.34 47 24 22.32 8.05 Isopoda Oligochaeta Lumbricidae total 17 SUPPORTING INFORMATION . Lumbricidae Allolobophoridella eiseni 11 3 9.82 2.30 Lumbricidae Not identified 36 21 16.96 8.05 8 9 2.68 5.75 NA 40,734 NA 75.86 Enchytraeidae ‚Acari‘ mainly Oribatidae 18 SUPPORTING INFORMATION Appendix S6: Community composition in tree-holes The change in community composition as a function of forest management is illustrated in the NMDS plots (Fig. S6-1). They show a higher relative abundance of filter feeders (green) in tree-holes of plots under higher forest management intensity. Fig. S6-1: Nonmetric multidimensional scaling (NMDS) of tree-hole communities (circles) and individual species (numbered crosses referring to species IDs Fig. 1 of the main manuscript) in the SA (stress value 0.185) and Hainich Dün (0.185). Forest management intensity is fitted post hoc as a vector (black arrow) to illustrate the location of plots with high or low forest management intensity in the ordination diagram. Green=filter feeders (FF), orange=fine detritus feeders (FD), blue=macrosaprophages (MS), red=predators (P). All sampled tree-holes with a full set of environmental variables were included (Schwäbische Alb: 112, Hainich-Dün: 87). 19 SUPPORTING INFORMATION Appendix S7: Detailed results from Structural Equation Models Table S7-1: Results from structural equation models on the indirect effects of forest management on insect communities in water-filled tree-holes. Covariances, R-squares, and total effects are given. For standardized path coefficients of regressions and overall model fit see Figure 5 and 6 of the main manuscript and Figure S7-1. Region was included as grouping factor in the Abundance-Richness and Community composition model. In the Hainich model additionally measured water properties and distance to the next tree-hole were included. Abundance-Richness Model Schw. Alb Covariances: Plant species richness ~~Insect richness ~~NMDS1 ~~NMDS2 Tree-hole density ~~Plant species richness ~~NMDS1 ~~NMDS2 NMDS1~~NMDS2 pH ~~Ammonium ~~Nitrate ~~Phosphate ~~O2 ~~ForMI 0.044 Community composition model Hainich-Dün n.s. Schw. Alb Hainich-Dün Hainich-Dün 0.349 -0.048 -0.051 -0.138 0.251 * 0.688 0.341 ** 2.414 -0.011 Hainich model -1.140 * -0.130 0.136 -0.280 0.173 0.105 0.226 * -0.236 -0.108 -0.167 0.130 0.257 0.343 * ** 0.136 0.337 -0.012 -0.182 * 0.111 0.143 0.447 -0.030 -0.268 * 20 SUPPORTING INFORMATION ~~Height Ammonium ~~Nitrate ~~Phosphate ~~O2 ~~ForMI ~~Height Nitrate ~~Phosphate ~~O2 ~~ForMI ~~Height ForMI~~Height 0.273 * 1.873 0.277 * 0.306 0.287 * 0.431 0.325 * 0.108 0.004 0.200 -0.284 0.148 0.144 -1.132 2.363 * 0.537 1.312 -0.321 * 0.007 0.192 -0.299 0.175 0.147 -0.159 0.355 * 0.066 0.224 -0.042 * 0.221 0.167 -0.071 -0.296 * -0.066 -0.107 0.041 Phosphate ~~O2 -0.003 ~~ForMI -0.362 O2~~ForMI -0.114 R-square: Water properties Tree-hole density Distance to next tree-hole Plant species richness Volume Detritus Insect abundance Insect richness NMDS1 NMDS2 1 0.558 1 0.330 1 0.575 1 0.331 0.099 0.042 0.529 0.578 0.532 0.010 0.005 0.589 0.393 0.616 0.102 0.048 0.517 0.009 0.009 0.572 0.374 0.254 0.384 0.262 * 0.325 0.214 0.174 0.072 0.559 0.532 0.571 21 SUPPORTING INFORMATION Total effects: Overall forest management effect on Insect abundance Insect species richness NMDS1 NMDS2 -3.021** -2.984** -0.151* -0.164* 0.042 n.s. 0.292* 22 SUPPORTING INFORMATION Figure S7-1: Indirect effects of forest management intensity on the abundance and species richness of insect communities in water-filled tree-holes of the Hainich-Dün region. In addition to the variables used in the model across regions distance to the next tree-hole (Distance) and additional water properties (Phosphate, O2 content) were included in the model 23 SUPPORTING INFORMATION which were not available for Schwäbische Alb. Top: Full hypothesized model; Bottom: Results from a structural equation model carried out with the lavaan package in R. Black solid arrows indicate paths of the final model. Values show standardized path coefficients for the two regions (Chi²= 66.179). Model structure did differ significantly from observed data (p=0.027) indicating low model fit. 24 SUPPORTING INFORMATION Appendix S8: Number of tree-hole types Fig. S8-1: Relationship of tree-hole density and number of tree-hole types per 1-ha plot (F1,36=28.64, p<0.001; adjusted R²=0.428). Insets show the number of tree-holes (bottom right; F1,51=36.88, p<0.001; adjusted R²=0.408) and number of tree-hole types (top left; F1,36=32.28, p<0.001; adjusted R²=0.458) as function of forest management intensity. Please note that the number of tree-hole types was only assessed in a subset of plots. 25 SUPPORTING INFORMATION Appendix S9: Diversity partitioning Table S9-1: Results from stand-level analyses on the differences of mean α-diversity, diversity, γ-diversity, and the proportion of -diversity (based on a multiplicative diversity partitioning of species richness) between low (ForMI <0.85; Schwäbische Alb 2 plots, Hainich-Dün 11 plots), medium (0.85x1.20; 5, 8) and high (>1.20; 7, 5) forest management intensity. Results are based on linear models. Df F p Direction of effect α-diversity Region ForMI Residuals 1 2 33 5.081 5.239 0.031 0.028 -diversity Region ForMI Residuals 1 2 34 0.263 0.342 0.612 0.563 γ-diversity Region ForMI Residuals 1 2 34 2.110 17.980 0.155 <0.001 ↓ Proportion of -diversity Region ForMI Residuals 1 1 49 0.337 5.082 0.565 0.031 Hainich-Dün>Schwäbische Alb ↓ ↑ 26 SUPPORTING INFORMATION Appendix S10: Detailed Acknowledgements We are grateful to Steffen Boch, Jörg Müller, Stephanie A. Socher, and Markus Fischer for providing the data on plant species richness, Brenda Guidetti for assistance in combining the data of the two regions, Norbert Höser for the identification of Lumbricidae, Reiner Bark and colleagues for their engagement in building and optimizing the tree-top camera, Kirsten Küsel for providing lab facilities for the analyses of water properties, Ulrich Kern for the graphical illustrations, and two anonymous reviewers for their valuable comments which greatly improved the manuscript. We thank the managers of the three Exploratories, Kirsten ReichelJung, Swen Renner, Katrin Hartwich, Sonja Gockel, Kerstin Wiesner, and Martin Gorke for their work in maintaining the plot and project infrastructure; Christiane Fischer and Simone Pfeiffer for giving support through the central office, Michael Owonibi for managing the central data base, and Markus Fischer, Eduard Linsenmair, Dominik Hessenmöller, Jens Nieschulze, Daniel Prati, Ingo Schöning, François Buscot, Ernst-Detlef Schulze and the late Elisabeth Kalko for their role in setting up the Biodiversity Exploratories project. The work was funded by the DFG Priority Program 1374 "Infrastructure-Biodiversity-Exploratories" (DFG-WE 3081/21-1). Field work permits were issued by the responsible state environmental offices of Baden-Württemberg and Thüringen. 27 SUPPORTING INFORMATION References Blakely, T.J. & Didham, R.K. (2010) Disentangling the mechanistic drivers of ecosystem-size effects on species diversity. Journal of Animal Ecology, 79, 1204-1214. Boch, S., Prati, D., Müller, J., Socher, S., Baumbach, H., Buscot, F. et al. (2013) High plant species richness indicates management-related disturbances rather than the conservation status of forests. Basic and Applied Ecology, 14, 496-505. Cranston, P.S. (1982) A key to the larvae of the British Orthocladiinae (Chironomidae). Scientific Publications of the Freshwater Biological Association, 45, 1-152. Cranston, P.S., Ramsdale, C., Snow, K. & White, G. (1987) Keys to the adults, male hypopygia, fourth-instar larvae and pupae of the British mosquitoes (Culicidae) with notes on their ecology and medical importance. Scientific Publications of the Freshwater Biological Association, 48, 1-152. Dixon, T. (1960) Key to and descriptions of the third instar larvae of some species of Syrphidae (Diptera) occurring in Britain. Transactions of the Royal Entomological Society of London, 112, 345-379. Grace, J. (2006) Structural equation modeling and natural systems. Cambridge University Press, Cambridge. Hartley, J. (1961) A taxonomic account of the larvae of some British Syrphidae. Proceedings of the Zoological Society of London, pp. 505-573. Wiley Online Library. Kitching, R.L. (1971) An ecological study of water-filled tree-holes and their position in the woodland ecosystem. Journal of Animal Ecology, 40, 281-302. Kitching, R.L. (2000) Food webs and container habitats: the natural history and ecology of phytotelmata. Cambridge University Press, Cambridge. Klausnitzer, B. (1996) Die Larven der Käfer Mitteleuropas: L6: Polyphaga: Teil 5. Gustav Fischer, Jena. 28 SUPPORTING INFORMATION Lindegaard, C. (1997) Diptera Chironomidae, non-biting midges. Aquatic Insects of North Europe, 2, 265-294. Mohrig, W. (1969) Die Culiciden Deutschlands: Untersuchungen zur Taxonomie, Biologie und Ökologie der einheimischen Stechmücken. Parasitologische Schriflenreihe, 18, 1260. Nilsson, A.N. (1996) Aquatic insects of North Europe: a taxonomic handbook. Volume 1: Ephemeroptera, Plecoptera, Heteroptera, Neuroptera, Megaloptera, Coleoptera, Trichoptera, Lepidoptera. Apollo Books. Nilsson, A.N. (1997) Aquatic insects of North Europe: a taxonomic handbook. Volume 2: Odonata, Diptera. Apollo Books. Pankratova, V.Y. (1970) Key to the larvae of the genera of the sub-family Orthocladiinae from Larvae and Pupae of midges of the sub-family Orthocladiinae of the fauna of the USSR. FBA Translations (New Series), 54, 51-55. Perry, D.R. (1978) A method of access into the crowns of emergent and canopy trees. Biotropica, 10, 155-157. Rosseel, Y. (2014) Package ‘lavaan’: Latent Variable Analysis Rotheray, G.E. (1993) Colour guide to hoverfly larvae (Diptera, Syrphidae) in Britain and Europe. Dipterists Digest, 9, 1-155. Rozkošný, R. & Gregor, F. (2004) Insecta, Diptera, Muscidae. Elsevier, München. Schmidl, J., Sulzer, P. & Kitching, R.L. (2008) The insect assemblage in water filled treeholes in a European temperate deciduous forest: community composition reflects structural, trophic and physicochemical factors. Hydrobiologia, 598, 285-303. Siitonen, J., Martikainen, P., Punttila, P. & Rauh, J. (2000) Coarse woody debris and stand characteristics in mature managed and old-growth boreal mesic forests in southern Finland. Forest Ecology and Management, 128, 211-225. Stehr, F.W. (1987) Immature insects [Vol. 1]. Kendall/Hunt Publishing Co. 29 SUPPORTING INFORMATION Stresemann, E., Klausnitzer, B., Hannemann, H.-J. & Senglaub, K. (2005) StresemannExkursionsfauna von Deutschland, Band 2: Wirbellose: Insekten. Spektrum Akademischer Verlag. Tylianakis, J.M., Rand, T.A., Kahmen, A., Klein, A.M., Buchmann, N., Perner, J. & Tscharntke, T. (2008) Resource heterogeneity moderates the biodiversity-function relationship in real world ecosystems. PLoS Biology, 6, 947-956. Utrio, P. (1976) Identification key to Finnish mosquito larvae (Diptera, Culicidae). Annales Agriculturae Fenniae (Finland), 15, 128-136. Yanoviak, S.P. (1999) Community structure in water-filled tree holes of Panama: effects of hole height and size. Selbyana, 20, 106-115. Yanoviak, S.P. & Fincke, O.M. (2005) Sampling methods for water-filled tree holes and their artificial analogues. Insect Sampling in Forest Ecosystems (ed. S.R. Leather), pp. 168185. Blackwell. 30
© Copyright 2024 ExpyDoc