Contrasting responses of multiple insect taxa to common heathland management regimes and old-growth successional stages

Maintaining heathlands in early successional stages to sustain heather ( Calluna vulgaris ) is a common, large-scale management practice in Europe. However, allowing patches of long-term natural vegetation development may increase habitat heterogeneity benefitting insects, but empirical evidence is sparse. We investigated how old-growth heathland (> 30 years abandonment) affect species richness and composition of bees (Anthophila), crane flies (Tipuloidea), ground beetles (Carabidae), hoverflies (Syrphidae) and rove beetles (Staphylinidae) in relation to their hygropreference. Adult insects, vegetation and edaphic explanatory variables were collected in old-growth, managed and wet sites and compared in four lowland heathland locations in Denmark. We found 299 species including 24 nationally red-listed. Species composition differed between managed, old-growth and wet heathland for all taxa. Indicator species and richness analyses showed a predominance of xerophilic bee species in managed heathland. Old-growth heathland showed a predominance of mesophilic indicator species, and higher richness of mesophilic crane flies and of hygrophilic ground and rove beetles compared to managed heathland. Wet heathland was generally dominated by hygrophilic species. Soil moisture, bare soil and vegetation height density were important drivers explaining the contrasting responses in richness and composition between heathland types. Our results demonstrate that heathland management focusing solely on early successional vegetation stages may homogenize insect communities. We suggest that management practices should focus on improving structural vegetation heterogeneity. This can be achieved through management regimes that reset the succession and expose bare soil, but also by allowing patches of old-growth vegetation stages to develop and by conserving existing ones.


Introduction
European heathlands are vulnerable semi-natural ecosystems of high conservation value protected under the European Union Habitats Directive (Evans 2006).Management of heathlands is necessary to sustain the cover and growth of ericoid dwarf shrubs and to prevent nutrient accumulation and succession towards forest.This is achieved through management regimes such as prescribed burning, livestock grazing or machine cutting, which deplete nutrients and favor rejuvenation and regeneration of heather (Calluna vulgaris) (Hobbs 1984;Webb 1998).Heathlands that are kept in their early successional vegetation stages are often associated with a higher cover of Calluna, but often develop into homogenous even-aged sites with uniform canopies and low structural diversity (Gimingham 1970).While most heathland management practices have focused on maintaining early successional stages, a growing body of literature shows negative consequences of a homogenous vegetation structure on heathland arthropods (Usher 1992;Bell et al. 2001;Webb et al. 2010;WallisDeVries et al. 2016;Hansen et al. 2020).
Heathlands contain a diverse and unique arthropod fauna with varied ecological requirements and life cycles (Kirby 1992;Usher 1992).Vegetation composition and structure are important drivers of variation in heathland arthropod assemblages.The cover and growth-phases of C. vulgaris (pioneer, building, mature and degenerate, total 30-40 years) are associated with different species assemblages of flies, ground beetles, butterflies and moths, millipedes, phytophagous insects, spiders and springtails (Miller 1975;Gimingham 1985;Usher 1992;Haysom and Coulson 1998).A high graminoid cover in heathlands promotes species richness of grasshoppers (Schirmel et al. 2011), while bare ground patches affect the species composition and community structure of bees, wasps and ground beetles (Gregory and Wright, 2005;Cameron and Leather 2012).An evaluation of the habitat requirements of UK heathland priority species, of which the majority are arthropods, found that shelter (topography or vegetation), bare ground, grasses, and shrubs/trees are essential habitats (Webb et al. 2010) contrasting the common heathland management aiming at a high and homogeneous Calluna cover.
The effects of heathland management on vegetation structure and composition and thus arthropod assemblages are complex and depend on several factors such as management regime and intensity, land-use history and abiotic factors (Gimingham 1975;Riis-Nielsen et al. 1991;Lake et al. 2001;Newton et al. 2009;Kepfer-Rojas et al. 2014;Schellenberg and Bergmeier 2021).Common management regimes such as prescribed burning aim to restart the life cycle of aging heather and prevent colonization of trees (Gimingham 1975).Prescribed burning removes the canopy and accumulated litter and exposes bare ground, which can be beneficial to a wide range of xerophilic epigeal and burrowing insects (Kirby 1992;Schirmel et al. 2012;Buchholz et al. 2013;Bargmann et al. 2015).Livestock grazing aims at increasing structural diversity, and maintaining open dwarf shrub vegetation with a high diversity of plant species while eliminating woody species (Lake et al. 2001).The effect of grazing on heathland arthropod assemblages is taxa specific and largely determined by grazing intensity (WallisDeVries et al. 2016).Intensive heathland management regimes may further affect soil and plant chemical properties with detrimental effects on the insect fauna (De Bruyn et al. 2001;Vogels et al. 2017).
Besides the effects of management, underlying abiotic factors also structure heathland ecosystems and may in turn affect insect communities.A key factor is soil moisture, and insect communities usually differ markedly between wet and dry heathlands (Usher 1992;De Bruyn et al. 2001;Mantilla-Contreras et al. 2011;WallisDeVries et al. 2016).
Wet heathland have often been unmanaged for decades or even centuries, and contain a large number of specialized hygrophilic arthropod species of high conservation value (Spitzer and Danks 2006).In large-scale intensively managed heathland locations, longterm unmanaged sites are not common, but allowing dwarf shrub-dominated heathlands to develop naturally through long periods of time (henceforth referred to as old-growth heathland) may benefit certain insect taxa.Over time, old-growth heathlands develop a complex vegetation structure and composition (Riis-Nielsen et al. 2005;Ransijn et al. 2015), including the different growth-phases of Calluna (Gimingham 1985).
The aim of this study was to investigate how old-growth heathland affect species richness and composition of five different insect taxa.Adult insects, vegetation and edaphic explanatory variables were collected in old-growth, managed and wet sites and compared in four lowland heathland locations in Denmark.Wet heathland was included as an outgroup, enabling us to compare the species assemblages to those of managed and oldgrowth heathland, as we expected old-growth heathland invertebrate communities to show more resemblance towards wet heathland.We focused on bees (Hymenoptera: Apoidea: Anthophila), crane flies (Diptera: Tipuloidea), ground beetles (Coleoptera: Carabidae), hoverflies (Diptera: Syrphidae) and rove beetles (Coleoptera: Staphylinidae), which are species rich families/clades (> 250 species per group in Denmark) commonly found in heathlands.Insect taxa were further divided into groups based on their hygropreferences.We expected that xerophilic species would prefer managed heathland because of the sparse vegetation and bare soil, mesophilic species would prefer old-growth heathland because of the denser vegetation and more varied conditions, and hygrophilic species would prefer wet heathland and to some extend old-growth heathland because of similarity in vegetation structure among habitats.To our knowledge, multi-taxa responses across heathland management regimes and old-growth unmanaged heathlands have not been previously reported.
We hypothesized that: (1) insect richness and species composition differ between heathland types (managed, old-growth and wet) and are linked to insect hygropreferences, with xerophilic species preferring managed heathland, mesophilic species preferring old-growth heathland, and hygrophilic species preferring wet heathland and to some extent also oldgrowth heathland; (2) differences in insect richness and composition between heathland types can be explained by differences in vegetation and edaphic conditions; (3) old-growth heathland sites in combination with managed sites provide heterogeneous habitat conditions that increase overall diversity of insects.

Study locations and experimental design
The study was conducted in four heathland locations in Jutland, Denmark (Fig. 1).The locations represent some of the largest inland heaths in Denmark, and were chosen based on historical information regarding past and present management regimes and long continuity.The heathlands are located within 70 km of each other, and are protected within the Natura 2000 network.The mean annual temperature is 8.8 °C and mean annual precipitation is 950 mm.The total annual nitrogen deposition is 11-15 kgN/ha (DCE 2020).
All heathland locations included (1) a dry old-growth heathland site primarily dominated by Calluna vulgaris or other dwarf shrubs (2) sites with different management regimes and (3) a wet heathland site (habitat code 4010) with Erica tetralix, Molinia 1 3 caerulea, Sphagnum species and different dwarf shrub species (Table 1).The managed sites have a previous history of intensive heathland management through the past decades, and have primarily been managed by prescribed burning, grazing or harvesting, which are among the most common management regimes in Denmark (Buttenschøn and Schmidt, 2015;Riis-Nielsen et al. 1991).
At each site, 8 plots were established randomly with the following criteria: Plots were placed with a minimum of 50 m to the border of the site, and with a minimum of 50 m between plots (Fig. 1) to avoid depletion of pitfall trap catches (Digweed et al. 1995).In addition, four extra plots were established in each site with pitfall traps installed.These traps were active in the same periods as in the other 8 plots, but were only used in the dataset if data from the 8 plots were incomplete, e.g.due to trampling of animals or flooding.

Sampling of insects
Insects were sampled recurrently from early June to early August of 2017 using sweep nets (crane flies and hoverflies), pan traps (bees and hoverflies), and pitfall traps (ground beetles and rove beetles).Sweep net sampling was conducted in a circle (r: 15 m) around the plot center for a total of 10 min shared by 2-3 collectors at three occasions during the summer (Online Appendix A) using a 38 cm diameter sweep net.Sweep net sampling was initiated  from the periphery of the circle to ensure least possible disturbance of the fauna.Sweep net catches were emptied in bags containing 70% ethanol.Sampling was conducted from 9 AM to 6 PM in dry weather and wind speed < 8 m/s.Three coloured pan traps (radius: 7.5 cm, volume: 500 ml) in blue, white and yellow, were placed on top of the vegetation or at ground level in each plot and filled with 50% propylene glycol and 50% soap water.Pan traps were active between 8 and 16 days in three periods (Online Appendix A).A pitfall trap was installed in each plot, consisting of a 50 cl plastic cup (diameter: 10 cm) with propylene glycol and a roof to prevent flooding.Pitfall traps were active between 12 and 30 days in three periods (Online Appendix A).

Species identification
Bees were identified to species using identification keys for the different bee genera (Amiet 1996;Amiet et al. 1999Amiet et al. , 2001Amiet et al. , 2004Amiet et al. , 2007;;Bogusch and Straka 2012;Schmid-Egger and Scheuchl 1997;Falk and Lewington 2015).Due to the issue of workers having identical appearance within the species complex Bombus terrestris, B. lucorum, B. cryptarum and B. magnus, specimens were pooled together and treated as one species (Bombus terrestris complex Group).Honey bees (Apis mellifera) only occur commercially in Denmark and were excluded from the dataset.Species identification was performed using the key by Stubbs and Kramer (2016) for crane flies; Luff (2007) for ground beetles; Bartsch et al. (2009) for hoverflies, with supplementary information on morphological characters from Van Veen (2010).Females of the genus Sphaerophoria were treated as Sphaerophoria female Group, since morphological identification is uncertain.Similarly, no distinction between species of Melanostoma scalare and Melanostoma mellinum was made and they were treated as Melanostoma Group.Rove beetles were identified using the keys by Hansen (1951;1952;1954) as well as Assing and Schülke (2012).Species in the morphological similar subfamily Aleocharinae were separated in two groups, one with the genus Dinarda, and one with the rest of the Aleocharinae, of which most belonged to a single species, Drusilla canaliculata.Following identification, species were assigned to xerophiles, mesophiles and hygrophiles based on their hygropreference.Bees were assigned according to Cederberg et al. (2022); crane flies according to information contained in the Catalogue of the Craneflies of the World (Oosterbroek 2022); ground beetles according to Lindroth (1985;1986); hoverflies according to information contained in Syrph the Net Database (Speight et al. 2000); and rove beetles according to Hansen et al. (2018).Red-listed species were assigned according to the Danish Red List (Moeslund et al. 2019), and expert assessments were used to evaluate the status of crane flies and rove beetles, which are not included in the Danish Red List.

Vegetation and edaphic variables
Vegetation analysis was conducted in the plots in a 5 m radius circle.The presence of all vascular plant species and a visual assessment of their cover including bare ground were recorded by estimates of two separate recorders.Vegetation height density was measured using a ruler in the plots at 1 m intervals across the plot in the cardinal directions, as the height where the ruler was no longer visible from a top-down perspective, thus representing a combination of height and density similar to Robel et al. (1970).Plant coverage data were assigned to functional groups (Table 2).
Soil was sampled 5 m from the plot center in the directions N, S, E and W using a soil corer with an inner diameter of 5 cm.The 4 samples were divided into two fractions, one comprising the mor-layer and another the top 10 cm of mineral soil (e-horizon) and pooled respectively.After drying at 55 °C, soil pH was measured on the mineral soil after 2 h extraction in a 0.01 M CaCl 2 solution on ratios 1:5 (soil:extractant) using a Radiometer combination-electrode GK2401 (Radiometer, Copenhagen, Denmark).Mineral soil for C:N analysis was further grinded on a Retsch Planetary Ball Mill PM400 (Retsch, Germany) at 370 rpm for 12 min.Carbon and nitrogen contents were measured using dry combustion method on a Thermo Scientific FLASH 2000 soil CN analyzer.Soil moisture content in the top layer was measured using a Theta ML3 Kit (Buch & Holm, Denmark) five meter from the plot center in the directions N, S, E and W at four different occasions during the summer and averaged.Edaphic variables are shown in Table 2.

Data analyses
Before analyses, insect species data from the three collection events were aggregated on plot level, and plots with incomplete data were omitted (Table 1).A total of 82 unspecified specimens were removed from the dataset, crane flies 16 specimens and hoverflies 66 specimens.All analyses were done with the statistical software R (R Core Team 2020), and graphs were created using the packages 'ggplot2' (Wickham 2009) and 'effects' (Fox et al. 2016).Sample completeness was evaluated on plot level, site level and heathland type level using the package 'iNEXT' (Hsieh et al. 2016).Because of low replication among the different management regimes (burned, grazed and harvested), they were treated as one group (managed) in the analyses (Table 1).This grouping was justified as we expected vegetation structure and composition to be relatively uniform among management regimes and different from old-growth and wet heathland.

Species richness
To assess the influence of heathland type (managed, old-growth and wet) on species richness of the different taxa according to their hygropreferences (xerophilic, mesophilic and hygrophilic), we performed generalized linear mixed effects models (GLMMs) using the 'lme4' package (Bates et al. 2007) or 'glmmTMB' package (Magnusson et al. 2017).
Heathland type (managed, old-growth and wet) was included as a fixed effect and heathland location as a random effect.Significance was assessed using the 'lmerTest' package (Kuznetsova et al. 2017), and in case of significance followed by a Tukey's post hoc test using the 'Multcomp' package (Hothorn et al. 2016).Assumptions of independence, homoscedasticity and normality of the residuals were assessed using the 'DHARMa' package (Hartig 2022).In most cases, a Poisson distribution for the count richness data was appropriate, but some models showed signs of underdispersion and a proper distribution was chosen.Model validation and full results are included in Online Appendix B.

Vegetation and edaphic variables
We compared if vegetation and edaphic variables differed between managed, old-growth and wet heathland.For vegetation cover, we performed GLMMs (Bates et al. 2007) with a negative binomial distribution.For edaphic variables we performed Linear mixed models (LMMs) (Bates et al. 2007) with a Gaussian distribution on log-transformed data.In case of significance, we used Tukey's post hoc test for pairwise comparisons (Table 2).
To assess the effect of vegetation and edaphic variables on species richness of the different insect taxa according to their hygropreference, we performed GLMMs using the 'lme4' package (Bates et al. 2007) or 'glmmTMB' package (Magnusson et al. 2017).Before analyses, collinearity of all measured predictor variables were assessed, and a generalized variance inflation factor (GVIF) with a threshold of 3 was used to remove collinear variables (Zuur et al. 2010).After exclusion of collinear variables, the remaining predictor variables were scaled to unit variance and included in GLMMs as fixed effects.Heathland location was used as a random effect.Model validation was done using the 'DHARMa' package (Hartig 2022).In most cases, a Poisson distribution for the count richness data was appropriate, but some models showed signs of underdispersion and a proper distribution was chosen.In addition, some models showed signs of non-linearity between dependent and independent variables.This was accounted for by visually inspecting relationships and adding a quadratic term to the predictor variable.Model validation and full results are included in Online Appendix C.

Species composition
Non-Metric Multidimensional Scaling (NMDS) was used to visualize the species composition of the insect taxa between heathland types using the function metaMDS in the 'vegan' package (Oksanen et al. 2016).Plots with fewer than three specimens were omitted to achieve a clear interpretation of data.We used a Bray-Curtis dissimilarity matrix calculated for Hellinger transformed abundance data.Ordinations were set to 999 iterations and stress values below 0.20 were accepted.Stress values above 0.20 were evaluated against the distribution of stress values from independent permutations of the dataset to assess its difference from a random ordination structure (Dexter et al. 2018).To assess the influence of vegetation and edaphic variables on species composition we used the envfit function 1 3 in the 'vegan' package (Oksanen et al. 2016).This function was fitted onto NMDS axes 1 and 2 for insect community data, and statistical significance was assessed using 9999 permutations.
Permutational multivariate analysis of variance (PERMANOVA; Anderson 2001) was used to test for differences in species composition between heathland types and locations using the Adonis function in the 'vegan' package (Oksanen et al. 2016).We independently conducted a PERMANOVA for heathland type (with permutations restricted within heathland locations using the argument strata) and heathland location.The pseudo p-values were generated using 9999 permutations.In the case of a significant PERMANOVA result, pairwise comparisons were performed for heathland types with Bonferroni-corrected p-values using the 'pairwiseAdonis' package (Arbizu 2017).As PERMANOVA can be sensitive to heterogeneity in dispersions (Anderson 2013), we tested for differences in dispersions among heathland types using the betadisper function in the 'Vegan' package (Oksanen et al. 2016).
To evaluate differences in species occurrence among heathland types we further performed indicator species analyses using the multipatt function in the 'indicspecies' package (De Cáceres et al. 2010).Indicator association values (IndVal) were calculated for each species related to either managed, old-growth and wet heathland or a combination, as the square root of the product of the specificity and the fidelity.We restricted permutations within heathland locations, and assessed statistical significance using 9999 permutations.

Results
In total, 14 702 identified specimens belonging to 299 species including 5 complexes of insects were found.Species diversity comprised between 6 and 26% of the total Danish species within each taxon.Xerophiles constituted 76 species, mesophiles 139 species and hygrophiles 84 species.18 species were DK red-listed species, and 6 species were considered red-listed/rare by expert opinion (Table 3).Sample completeness was generally high for bees, ground beetles, hoverflies and rove beetles (plot level x: 0.82; site level x: 0.95; type level x: 0.99) (Online Appendix D), while crane flies had too few specimens in managed and old-growth sites to make any robust estimation.Roughly 1/3 of all species were singletons and doubletons.A full species list is found in Online Appendix E.

Richness of insect taxa
Xerophilic species richness varied among insect taxa and heathland types with plot median values ranging from 2 to 6 for bees, 0 for crane flies, 0-1 for ground beetles, 2-3 for hoverflies, and 0-2 for rove beetles (Fig. 2).Bees were the only taxon showing a significantly higher richness of xerophilic species in managed plots compared to old-growth plots, while xerophilic bees, ground beetles and rove beetles showed a significantly higher species richness in old-growth plots compared to wet plots (Fig. 2).
Mesophilic species richness was generally higher than that of xerophilic species richness with median values ranging from 3 to 4 for bees, 0-1 for crane flies, 2-7 for ground beetles, 5 for hoverflies, and 1-5 for rove beetles (Fig. 2).Crane flies were the only taxon showing a significantly higher richness of mesophilic species in old-growth plots compared to managed plots, while both mesophilic ground beetles and rove beetles showed a significantly lower species richness in wet plots compared to managed and old-growth plots (Fig. 2).
Hygrophilic species richness showed median values ranging from 0 for bees, 0-3 for crane flies, 0-2 for ground beetles, 2-3 for hoverflies, and 1-3 for rove beetles (Fig. 2).Hygrophilic ground beetle and rove beetle richness were significantly higher in old-growth plots compared to managed plots.Hygrophilic crane fly and ground beetle richness were significantly higher in wet plots compared to managed and old-growth plots, while no difference in hygrophilic rove beetle richness was found between old-growth and wet plots (Fig. 2).
In addition, the test was conducted on the number of red-listed species across all taxa, and no difference was found between heathland types (GLMM: Chi = 0.01, P = 0.99).

Insect richness in relation to vegetation and edaphic variables
Overall, the measured explanatory variables differed between heathland types (Table 2).In total, three variables were removed by the GVIF and excluded from the models (Table 2).A Pearson correlation matrix is included in Online Appendix F.  NS not significant, NA not available 1 3 Xerophilic bee richness was negatively associated to veg.height density (GLMM: est = − 0.28 ± 0.11, Z = − 2.6, p = 0.01) and showed a nonlinear negative association to soil moisture (GLMM: est = − 0.28 ± 0.12, Z = − 2.3, p = 0.02) (Fig. 3), while mesophilic bee richness showed no significant association to any variable.
Richness of xerophilic, mesophilic and hygrophilic hoverflies showed no significant association to any variable.

Species composition
Species compositions of all insect taxa were significantly different between heathland types, explaining 12-28% of the variation (Table 4), whereas heathland location was significant for all taxa and explained 20%, 18%, 23% and 26% for bees, ground beetles, hoverflies and rove beetles, respectively.Crane flies were omitted from species composition analyses because of a limited number of specimens in managed and old-growth plots.Bees showed a stress value of 0.25, but the ordination was considered different from a random ordination structure (Online Appendix G).
Pairwise tests showed a significant difference in the species composition of all insect taxa in all possible comparisons, except for hoverflies which were not different between managed and old-growth heathland (Table 4; Fig. 4).However, bees, ground beetles and hoverflies violated the assumptions of similar multivariate dispersions (betadisper test) indicating that differences in β-diversity between heathland types might be causing the differences in composition between heathland types.The ordinations show that this is the case only for bees, where there is an overlap between old-growth and wet heathlands.
Several vegetation and edaphic variables were significantly associated to the NMDS axes 1 and 2 (Fig. 4).Insect assemblages in managed plots were generally associated with exposed bare soil and high pH, whereas old-growth plots were associated with other dwarf shrub and forb cover.Wet plots were associated with high soil moisture.In addition, both old-growth and managed plots were associated with the cover of Calluna, while wet and old-growth plots were associated with vegetation height density (Fig. 4).
Indicator species analysis identified 23, 18 and 19 species significantly associated to managed, old-growth and wet heathland, respectively (Table 5).The species associated to managed heathland primarily comprised xerophilic species, while old-growth heathland primarily comprised mesophilic species, and wet heathland was dominated by hygrophilic species (Table 5).Full indicator species analysis is included in Online Appendix H.

Discussion
Despite the strong dependence of lowland European heathlands on management, few studies have contrasted the impact of managed heathland to old-growth succession heathland on insect communities.Our analyses across multiple insect taxa showed that species richness and composition generally differed between heathland types (managed, old-growth and wet heathland), and that differences were linked to shifts in insect hygropreferences.
The measured vegetation and edaphic variables to some extent explained the differences in species richness and composition between heathland types.However, the results were taxaspecific and contrasted among taxa.We hypothesized that xerophilic species would show a preference for managed heathland because of the sparse vegetation and bare soil.This was, however, only the case for xerophilic bees, which showed a much higher richness and number of xerophilic indicator species compared to old-growth and wet heathland.Xerophilic bee assemblages consisted predominantly of ground-nesting species, which are highly affected by soil conditions, where soil moisture (Julier and Roulston 2009;Pane and Harmon-Threatt 2017) and bare soil (Gregory and Wright 2005;Potts et al. 2005) are important factors.This corresponded with the analysis, where xerophilic bee species richness in managed heathland was explained by low soil moisture and low vegetation height density.
Surprisingly, xerophilic ground and rove beetles showed no difference in richness between managed and old-growth heathland, and only a few xerophilic ground and rove beetle indicator species were found.Xerophilic ground beetles have previously been associated with high species prevalence and many indicator species in young successional vegetation stages of dry habitats (Schirmel and Buchholz 2011;Buchholz et al. 2013).This indicates that xerophilic epigeic heathland insect communities may be able to persist after 30-120 years of abandonment, or this could simply be a spillover effect from nearby managed sites.Some rarely encountered xerophilic species were even unique to old-growth areas such as the rove beetle Platydracus latebricola found in Nørholm.It should, however, be noted that the investigated old-growth sites were primarily dominated by dwarf-shrubs with low to intermediate levels of grass and tree encroachment.Heathland dominated by grasses and trees would likely have a negative effect on xerophilic beetle communities (Desender and Turin 1989;Schirmel and Buchholz 2011;Buchholz et al. 2013).In accordance with our hypothesis, old-growth heathland showed a predominance of mesophilic indicator species, and higher richness of mesophilic crane flies and hygrophilic ground and rove beetles compared to managed heathland.Mesophilic crane fly richness showed a weak negative relationship to Calluna cover and a positive relationship to other dwarf shrub cover, which provides little biological explanation.It is thus more likely, that mesophilic crane flies simply benefit from the accumulation of organic matter and mosses in old-growth areas, as many species are considered detritivores and bryophagous (Pritchard 1983).In addition, rove beetles have a high affinity to microhabitats, especially congregated piles of various types of organic matter, e.g.litter, dung, fungi and animal nests (Betz et al. 2018;Hansen et al. 2018), and early successional heathland stages with bare soil may lack such substrates.This is in accordance with mesophilic and hygrophilic rove beetle richness showing a negative relationship to bare soil and a positive relationship to vegetation height density, respectively.
While species richness was generally low in wet heath compared to managed and oldgrowth heathland, most taxa except bees had distinct assemblages of hygrophilic species, with a comparable number of indicator species as managed and old-growth heathland.Crane flies were highly abundant and species rich in wet areas compared to dry areas, which can be attributed to the hydrophilic nature of the larvae (Pritchard 1983;Byriel et al. 2020).Some rare crane fly species findings included the peat bog-specialists Phylidorea heterogyna and Erioptera nielseni, which were only known from a few old records in Denmark and were considered threatened or even extinct (Salmela 2012), but occurred in high abundances in three of the four wet heath sites.Bee communities in wet heathland resembled communities in old-growth heathland as indicated by the NMDS plot.Contribution to similarity could include an overlap in floral ericoid resources, e.g.Vaccinium spp., which are attractive to Bombus spp.(Moquet et al. 2017) and the oligolectic Andrena lapponica, which was identified as an indicator species of old-growth and wet heathland.
As opposed to the other taxa investigated, hoverflies showed remarkably similar richness median values and no response to vegetation and edaphic variables.Species composition indicated low variation among managed and old-growth heathland, while wet heathland showed distinct assemblages.In contrast to bees, hoverflies do not have parental care and are not restricted to a certain home range after oviposition (Rotheray and Gilbert 2011), which together with their high mobility makes them more suitable for environmental evaluation at larger scales (Sommaggio 1999).Common heathland hoverfly assemblages include long-foraging aphidophagous species associated with e.g.coniferous and deciduous forests (Moquet et al. 2018), thus explaining the low affinity to local vegetation characteristics.Wet heathland may show distinct hoverfly communities as they provide breeding habitats for saprophagous and aquatic larvae (Moquet et al. 2018), which corresponded with the three hygrophilic indicator species identified from wet heathland.
Several studies have shown that heathland management regimes such as grazing, harvesting and burning influence the structure of arthropod communities (Bell et al. 2001;Bargmann et al. 2015;WallisDeVries et al. 2016;Hansen et al. 2020).Although our study did not allow for a robust comparison between different management regimes, we observed taxa-specific responses to management as inferred from the NMDS plots.Especially insect assemblages in harvested/burned plots varied from other management regimes.These plots were characterized by a high cover of bare soil and forbs, which are known to influence pollinating and soil-dwelling insects through changes in food availability, temperature and nesting conditions (Gregory and Wright 2005;Meyer et al. 2009;Cameron and Leather 2012;Schirmel et al. 2012;Moquet et al. 2018).Grazed heathland was similarly characterized by distinguished hoverfly and rove beetle communities, and none of the measured variables seemed to explain this relationship.Large grazing herbivores are known to affect insect communities directly by creating microhabitats such as dung or trampled soil (Dennis et al. 1997;van Klink et al. 2020), which are important substrates for many species of hoverflies and rove beetles e.g. the dung-associated hoverfly Eristalinus sepulchralis found as an indicator of managed heathland.

Study limitations
Because of low replication among the different management regimes (burned, grazed and harvested), we were unable to elucidate the effect of individual management regimes on insect taxa diversity.Therefore, the results presented here should be viewed as an initial exploration of the mechanisms separating insect communities in managed vs. old-growth succession heathland.A balanced dataset including the different management regimes would likely provide a more nuanced depiction of insect diversity and associated drivers across managed and old-growth heathland.Despite this limitation, our study demonstrated a clear effect of the different heathland types on the investigated insect taxa, highlighting the importance of old-growth successional stages in structuring insect communities.Future studies should disentangle the effects of individual management regimes on the investigated insect taxa, as it would be of great conservation value.
The sampling conducted in this study was done at three different occasions from early June to early August.This interval was chosen as we expected it would represent the main peak in heathland insect diversity across all taxa.However, the interval may have been too short to cover the full phenology of all the insect taxa investigated.Sampling in the spring and autumn would unquestionably have yielded additional species, but we were limited by the number of samples that we would be able to sort and identify.

Implications for conservation management
The high number of insect species recorded as well as the large proportion of nationally red-listed species (~ 8% of all identified species) highlight the significance of heathlands in insect conservation and the value of preserving this semi-natural habitat.With our multitaxa approach, we found supporting evidence for the hypothesis that old-growth heathland sites in combination with managed heathland sites provide heterogeneous habitat conditions that increase overall diversity of insects.Heathland management focusing solely on maintaining early successional stages with a high cover of even-aged Calluna may be beneficial to certain xerophilic insect species, but may limit a large proportion of the mesophilic and hygrophilic heathland insect fauna, which require more buffered conditions provided by the dense vegetation and humid microclimate.
Our study challenges the general conception that heathland succession deteriorate xerophilic insect communities and endangered species as previously reported for ground beetles (Buchholz et al. 2013) and butterflies (Schirmel and Buchholz 2011).Hence, we found no difference in red-listed species between managed, old-growth and wet heathland, and xerophilic species seemed to persist to some extent in old-growth heathland sites.In addition, no association between red-list status and hygropreference was found in this study with 10 xerophilic, 9 mesophilic and 5 hygrophilic species being red-listed.In contrast to Schirmel and Buchholz (2011) and Buchholz et al. (2013) this study did not investigate succession from dry heathland to grass-invaded and tree-dominated sites, which likely explain the discrepancies.Additional studies are needed to understand to what extent xerophilic and endangered heathland species may persist following heathland succession.
Overall, our results demonstrate that heathland management focusing on early successional vegetation stages may homogenize insect communities, emphasizing the need for rethinking current management practices.We suggest that management practices should focus on improving structural vegetation heterogeneity.This can be achieved through management regimes that reset the succession and expose bare soil, but also by allowing patches of old-growth vegetation stages to develop and by conserving existing ones.

Fig. 2
Fig. 2 Species richness of xerophilic (Xero), mesophilic (Meso) and hygrophilic (Hygro) insect taxa in managed, old-growth and wet plots.Black dots indicate outliers.The letters indicate statistically significant difference at p < 0.05 between managed, old-growth and wet heathland for the different hygropreferences.NS not significant, NA not available

Table 1
Heathland locations, sites and typesYear in brackets refer to year of management abandonment.The no. of complete plots for the different taxa in the different heathland sites used in analyses are shown

Table 2
Vegetation and edaphic variables measured in the heathland plots (n = 126)Mean and SE values are shown for the different heathland types.Variables were included in richness models and NMDS graphs following a generalized variance inflation factor (GVIF) with a threshold of 3. Different letters indicate significant differences between heathland types (Tukey's test, p < 0.05)

Table 3
Diversity and abundance of insect taxa in the four heathlands Red-listed species in brackets were based on expert assessments

Table 4
Results from the PERMANOVA and pairwise PERMANOVA showing the effect of heathland type (man: Managed.old: Old-growth.wet: Wet) on species composition of the different insect taxa Distance: Bray-Curtis.Based on 9999 permutations a Violation of similar multivariate dispersions

Table 5
Summary of the IndVal indicator species analysesResults show the total no. of species significantly associated to managed, old-growth, and wet heathland or the combination