1. Introduction
Tectonic processes are fundamental drivers of global, long-term biodiversity patterns in all ecosystems (e.g., [
1,
2]). Understanding these processes, e.g., for the functioning of forests as one of the most important and widespread, terrestrial ecosystems, is of fundamental importance for ecology and biology.
Forests are long-lived ecosystems with multiple functions: a) highly important economic factor and raw material supplier for forest-based industries, b) climate regulator, c) protector against and moderator of natural hazards, d) conservator of biodiversity, natural cycles and provider of habitats for flora and fauna, and e) provider of natural experiences and recreational areas for humans [e.g., 3,4]
In Germany, approx. 11.4 million ha are forested [
5]. With approx. 2.56 Mio ha or one third of its area, Bavaria is of the most densely forested states in Germany. In Bavaria the responsibility for the forest is divided among different owners: 57 % is privately owned (private forests; PF), 30 % is owned by the Free State Bavaria (Bavarian State Forest; BSF), 11 % is owned by municipalities (municipal forests; MF) and 2 % is owned by the federal government [
6]. The BSF is the main forest owner in the two re-investigated study sites Falkenberg (FB) and Münchsgrün (MG), both under the responsibility of the forest division Waldsassen; smaller shares belong to other forest owners, e.g., private or municipal forests [
7].
Many natural (e.g., climate change, loss of biodiversity) and management influences (e.g., intensive use) affect forest resilience, i.e., the ability of forests to recover after disturbances [
8]. Sustainable forest management (SFM) aims to support the forest and the multifunctional forest sector to a) secure and improve the growth and stability of forest stands, b) develop sustainable, economically and ecologically balanced management practices, and c) ensure the long-term integrity of forest ecosystems [
3]. Sustainability and regional conservation concepts, such as the one developed for the Waldsassen Forestry Division in 2010 and updated in 2019, describe a large range of conservation measures as an integral part of near-natural forest management. The primary objectives are to preserve and promote biodiversity, and to create site-appropriate, near-natural, climate-resilient, productive and species-rich mixed forests. In addition, more than one third of the state forest area in Bavaria is currently classified as at least of the following protection categories: Natura 2000, nature reserves or natural forests [
6,
9,
10].
Dead wood (DW) is considered an indicator of sustainable forest management. Dead wood is not only recognized as an important habitat for a wide range of organisms, but also play key roles in carbon, nutrient, and hydrological cycles and influences ecosystem processes. In most managed forests, DW such as decaying standing or downed trees, roots, and branches are scarce due to their removal [
11].
Red wood ants (
Formica rufa-group; hereafter RWA), a key ecological group of forest species (e.g., [
12]), interact with their environment in the most diverse ways, e.g., contributing to habitat biodiversity (e.g., [
13]), regulating pests (e.g., [
14]), and are bio-indicators of undetected tectonic activity [
15,
16,
17,
18,
19] and geogenic gases (“GeoBio-Interactions”), which play a critical role in their settlement [
20,
21,
22,
23]. Declines in insect diversity, species, and biomass are attributed to, for example, habitat loss, invasive species or climate change (e.g., [
24,
25]). For RWA, declines (e.g., [
26,
27]) but also population increases are reported(e.g., [
28,
29]). RWA are also considered to be endangered in Germany, although statistically sound long-term surveys are lacking, as continuous, standardized and systematic monitoring of RWA has been discontinued since the 1980s [
30,
31]. Therefore, it is currently not possible to make any statements about an urgently needed protection status for this species.
For a comparative analysis of presence/absence data of RWA nests with respect to different forest factors, we applied our developed, area-wide, systematic, reproducible, and integrated approach [
7,
18,
32,
33] and re-inventoried two study sites, a) Münchsgrün (MG), and b) Falkenberg (FB) after seven (MG) and four years (FB), respectively [
7], in the tectonically active Oberpfalz (NE Bavaria, Germany).
In this approach, not only the number of RWA nests was counted, but also the entire ecosystem around a RWA nest, e.g., tree species, tree age, natural restocking, dead wood, herb layer, and woodpecker cavities, was monitored and documented in a photo database [
7,
18,
32]. We asked six interrelated questions: (1) Are presence/absence data of RWA nests comparable to the two previous inventories? (2) What influence do the aforementioned forest factors have on RWA nests? (3) Does dead wood (standing and/or downed trees, lying branches) in, on and around a nest influence RWA nest settlements?, (4) Do woodpeckers influence RWA nests?, (5) What influence do tectonic structures have on RWA nest distribution (GeoBio-Interactions), and (6) What are the best time intervals for re-inventories? These results will further improve the understanding of the “GeoBio-Interactions”, contribute to better protection of RWA in forest management in the Oberpfalz, and to the German government’s Insect Conservation Action Program [
34].
Figure 1.
Tectonic setting of both study sites with (
a), major tectonic units, faults (black lines), earthquake events (blue dots) taken from literature [
35,
36,
37,
38,
39,
40,
41,
42,
43]. Inlet shows location of both sites within Bavaria, Germany close to the Czech border; (
b) Münchsgrün (MG) and (
c)Falkenberg (FB) study areas detailing mapped sites for the 2016, 2019 and 2023 inventories [
7].
Figure 1.
Tectonic setting of both study sites with (
a), major tectonic units, faults (black lines), earthquake events (blue dots) taken from literature [
35,
36,
37,
38,
39,
40,
41,
42,
43]. Inlet shows location of both sites within Bavaria, Germany close to the Czech border; (
b) Münchsgrün (MG) and (
c)Falkenberg (FB) study areas detailing mapped sites for the 2016, 2019 and 2023 inventories [
7].
Figure 2.
Gain and loss of RWA nests (ntot) for (a) MG and (b) FB study area for the 2016/2019 inventories and the 2023 re-inventory.
Figure 2.
Gain and loss of RWA nests (ntot) for (a) MG and (b) FB study area for the 2016/2019 inventories and the 2023 re-inventory.
Figure 3.
Results of MANOVA for active nest parameters NH, ND, medium tree age (TSprime) and woodpecker cavities (WpC) showing a grouped plot matrix for the 2016, 2019 and 2023 inventories for (a) MGBSF and (b) FBBSF. Tree age with a -100 signature represent no trees or clearings around a nest.
Figure 3.
Results of MANOVA for active nest parameters NH, ND, medium tree age (TSprime) and woodpecker cavities (WpC) showing a grouped plot matrix for the 2016, 2019 and 2023 inventories for (a) MGBSF and (b) FBBSF. Tree age with a -100 signature represent no trees or clearings around a nest.
Figure 4.
Visual representation of the forest composition (log-normed) observed in the field at active RWA nests (nact): pure deciduous trees (DT), pure coniferous trees (CT), combination of deciduous trees and coniferous trees (DT & CT), and no trees (none) for the 2016, 2019 and 2023 inventories for a) MG and b) FB. Legend – CT included: pine (Pinus sylvestris), spruce (Picea abies), larch (Larix decidua), fir (Abies alba), and Douglas fir (Pseudotsuga menzisii). DT included: alder (Alnus glutinosa), beech (Fagus silvatica), birch (Betula pendula), bloody dogwood (Cornus sanguinea), common hazel (Corylus avellana), common hornbeam (Carpinus betulus), red elderberry (Sambucus racemosa), rowan (Sorbus aucuparia), oak (Quercus robur), and willow (Salix).
Figure 4.
Visual representation of the forest composition (log-normed) observed in the field at active RWA nests (nact): pure deciduous trees (DT), pure coniferous trees (CT), combination of deciduous trees and coniferous trees (DT & CT), and no trees (none) for the 2016, 2019 and 2023 inventories for a) MG and b) FB. Legend – CT included: pine (Pinus sylvestris), spruce (Picea abies), larch (Larix decidua), fir (Abies alba), and Douglas fir (Pseudotsuga menzisii). DT included: alder (Alnus glutinosa), beech (Fagus silvatica), birch (Betula pendula), bloody dogwood (Cornus sanguinea), common hazel (Corylus avellana), common hornbeam (Carpinus betulus), red elderberry (Sambucus racemosa), rowan (Sorbus aucuparia), oak (Quercus robur), and willow (Salix).
Figure 5.
Visual representation of the qualitative composition of main herbs (log-normed) around and on active RWA nests (nact) for the 2016, 2019 and 2023 inventories for a) MG and b) FB. Legend – Highly abundant herbs included Blb: European blueberry [Vaccinium myrtillus], CrB: cranberry [Vaccinium vitis-idaea], Fer: eagle fern [Pteridium aquilinum], FoG: foxglove [Digitalis purpurea], Gr: grass [Poaceae], Mo: moss [Bryophta], and RGr: reed grass [Calamagrostis]. Herbs categorized as difH: bellflowers (Campanula), blackberry (Rubus), broom (Genista), cattail (Typha), chickweed (Stellaria media), cleavers (Gallium aparine), coltsfoot (Tussilago farfara), cranesbills (Geranium pratense), dandelions (Taraxacum officinale), dead-nettles (Lamium), field pansy (Viola tricolor), field pennycress (Thlaspi arvense), flatweed (Hypochaeris radicata), ground elder (Aegopodium podagraria), groundsel (Senecio vulgaris), hairy bittercress (Cardamine hirsuta), horsetail (Equisetum), lady’s fern (Athyrium filix-femina), lupin (Lupinus), May lily (Maianthemum bifolium), melde (Chenopodium album), mullein (Verbascum), nettles (Urtica), red dead nettle (Lamium purpureum), red sorrel (Rumex acetosella), rushes (Juncus acutus), St John’s wort (Hypericum perforatum), thale cress (Arabidopsis thaliana), thistle (Cirsium vulgare), wild strawberry (Fragaria vesca), wood sorrel (Oxalis acetosella), and Yarrows (Achillea filipendulina) and (Achillea millefolium).
Figure 5.
Visual representation of the qualitative composition of main herbs (log-normed) around and on active RWA nests (nact) for the 2016, 2019 and 2023 inventories for a) MG and b) FB. Legend – Highly abundant herbs included Blb: European blueberry [Vaccinium myrtillus], CrB: cranberry [Vaccinium vitis-idaea], Fer: eagle fern [Pteridium aquilinum], FoG: foxglove [Digitalis purpurea], Gr: grass [Poaceae], Mo: moss [Bryophta], and RGr: reed grass [Calamagrostis]. Herbs categorized as difH: bellflowers (Campanula), blackberry (Rubus), broom (Genista), cattail (Typha), chickweed (Stellaria media), cleavers (Gallium aparine), coltsfoot (Tussilago farfara), cranesbills (Geranium pratense), dandelions (Taraxacum officinale), dead-nettles (Lamium), field pansy (Viola tricolor), field pennycress (Thlaspi arvense), flatweed (Hypochaeris radicata), ground elder (Aegopodium podagraria), groundsel (Senecio vulgaris), hairy bittercress (Cardamine hirsuta), horsetail (Equisetum), lady’s fern (Athyrium filix-femina), lupin (Lupinus), May lily (Maianthemum bifolium), melde (Chenopodium album), mullein (Verbascum), nettles (Urtica), red dead nettle (Lamium purpureum), red sorrel (Rumex acetosella), rushes (Juncus acutus), St John’s wort (Hypericum perforatum), thale cress (Arabidopsis thaliana), thistle (Cirsium vulgare), wild strawberry (Fragaria vesca), wood sorrel (Oxalis acetosella), and Yarrows (Achillea filipendulina) and (Achillea millefolium).
Figure 6.
Examples of dead wood classes 1-3 in (a) MG and (b) FB study area: DW-1: Biotope trees (weathered red signature “A” in MG is a marker for the RWA nest) and dying, standing tree, DW-2: branches of different size and diameter as downed DW, and DW-3: downed DW with additionally tree stems. Photo credit: M.B. Berberich.
Figure 6.
Examples of dead wood classes 1-3 in (a) MG and (b) FB study area: DW-1: Biotope trees (weathered red signature “A” in MG is a marker for the RWA nest) and dying, standing tree, DW-2: branches of different size and diameter as downed DW, and DW-3: downed DW with additionally tree stems. Photo credit: M.B. Berberich.
Figure 7.
Visual representation of nest height (NH) classes 1–5, woodpecker classes (WpC 1-5, WpC 6-10, WpC>10), tree species and medium tree age (mTA) versus numbers of active nests (nact), and for (a,b) MGBSF and (c,d) FBBSF for the 2016, 2019 and 2023 inventories.
Figure 7.
Visual representation of nest height (NH) classes 1–5, woodpecker classes (WpC 1-5, WpC 6-10, WpC>10), tree species and medium tree age (mTA) versus numbers of active nests (nact), and for (a,b) MGBSF and (c,d) FBBSF for the 2016, 2019 and 2023 inventories.
Figure 8.
Density plots of active RWA nests (n
act) in (
a) MG and (
b) FB study area and (
c), tectonic stress directions at both sites (yellow; © World Stress Map 2016 [
47]). Data for 2016 (MG) and 2019 (FB) were taken from [
7].
Figure 8.
Density plots of active RWA nests (n
act) in (
a) MG and (
b) FB study area and (
c), tectonic stress directions at both sites (yellow; © World Stress Map 2016 [
47]). Data for 2016 (MG) and 2019 (FB) were taken from [
7].
Table 1.
Descriptive statistics of total nest numbers (ntot), numbers of active nests (nact), numbers of active re-identified nests (ntotR), nest height (NH) and diameter classes (ND) for the MG (2016/2023) and FB (2019/2023) inventories. Increase in active nest numbers (Δ nact) and percentages are set in bold. – = not present.
Table 1.
Descriptive statistics of total nest numbers (ntot), numbers of active nests (nact), numbers of active re-identified nests (ntotR), nest height (NH) and diameter classes (ND) for the MG (2016/2023) and FB (2019/2023) inventories. Increase in active nest numbers (Δ nact) and percentages are set in bold. – = not present.
Year |
Study Area |
Numbers |
Difference for MG (2016/2023)and FB (2019/2023)
|
a) Nest Height (NH) Classes of Active Nests (nact) |
b) Nest Diameter (ND) Classes of Active Nests (nact) |
Start-Ups |
Short |
Medium |
Tall |
Very Tall |
Small |
Medium |
Large |
Very large |
Extra-large |
ntot
|
nact
|
nactR
|
Δ ntot
|
% |
Δ nact
|
% |
Δ nactR%
|
0.01–0.10 |
0.11–0.50 |
0.51–1.00 |
1.01–1.50 |
1.51–2.00 |
0.01–0.50 |
0.51–1.00 |
1.01–1.50 |
1.51–2.00 |
>2.01 |
2016 |
MG |
2326 |
2292 |
– |
– |
– |
– |
– |
– |
277 |
1208 |
632 |
153 |
22 |
947 |
696 |
393 |
173 |
83 |
2023 |
MG |
2555 |
2513 |
1336 |
229 |
9.8 |
221 |
9.6 |
58.3 |
117 |
1496 |
706 |
165 |
29 |
836 |
829 |
594 |
199 |
55 |
2019 |
FB |
2830 |
2607 |
– |
– |
– |
– |
– |
– |
406 |
1453 |
607 |
138 |
3 |
1175 |
738 |
493 |
129 |
72 |
2023 |
FB |
2838 |
2763 |
1712 |
8 |
0.3 |
156 |
6.0 |
65.7 |
353 |
1683 |
622 |
96 |
9 |
1131 |
822 |
591 |
168 |
51 |
Table 2.
Descriptive statistics of total mapped area (hatot;), forest owners (BSF, MF, PF) that hold share of the mapped area, numbers of active nests (nact) in BSF, MF, PF, and Fauna-Flora-Habitat (FFH) areas, numbers of active nests (nact) in natural restocking (NR) and on partly cleared and clearing plots (CP) for the 2016, 2019 and 2023 inventories for (a) MG and (b) FB study area; – = not present.
Table 2.
Descriptive statistics of total mapped area (hatot;), forest owners (BSF, MF, PF) that hold share of the mapped area, numbers of active nests (nact) in BSF, MF, PF, and Fauna-Flora-Habitat (FFH) areas, numbers of active nests (nact) in natural restocking (NR) and on partly cleared and clearing plots (CP) for the 2016, 2019 and 2023 inventories for (a) MG and (b) FB study area; – = not present.
Year |
Mapped area |
State forest (BSF) |
Municipal forest (MF) |
Private forest (PF) |
FFH area* |
Number of active nests (nact) |
Number of nact in BSF, MF & PF for |
|
hatot
|
ha |
% |
ha |
% |
ha |
% |
ha |
% |
∑nact
|
BSF |
MF |
PF |
FFH* |
NR |
CP |
a) MG |
2016 |
149 |
128 |
85.9 |
7.4 |
5.0 |
13.6 |
9.1 |
7.7 |
5.2 |
2292 |
2110 |
89 |
93 |
89 |
44 |
120 |
2023 |
149 |
131 |
87.9 |
6.8 |
4.6 |
11.2 |
7.5 |
7.7 |
5.2 |
2513 |
2213 |
121 |
179 |
105 |
280 |
357 |
b) FB |
2019 |
200 |
167 |
83.5 |
– |
– |
33 |
16.5 |
– |
– |
2607 |
2221 |
– |
386 |
– |
331 |
57 |
2023 |
203 |
167 |
82.3 |
– |
– |
36 |
17.7 |
– |
– |
2763 |
2332 |
– |
431 |
– |
294 |
182 |
Table 3.
Descriptive statistics of total mapped area (hatot;) for forest owners (BSF, MF, PF), numbers of active nests (nact) in BSF, MF, PF, numbers of re-identified active nests (nactR) in BSF, MF, PF, dead wood classes 1– 3 (DW), numbers of active nests with DW (nactDW), and ratio of active nests with dead wood per ha for the 2016, 2019 and 2023 inventories for (a) MG and (b) FB. – = not present.
Table 3.
Descriptive statistics of total mapped area (hatot;) for forest owners (BSF, MF, PF), numbers of active nests (nact) in BSF, MF, PF, numbers of re-identified active nests (nactR) in BSF, MF, PF, dead wood classes 1– 3 (DW), numbers of active nests with DW (nactDW), and ratio of active nests with dead wood per ha for the 2016, 2019 and 2023 inventories for (a) MG and (b) FB. – = not present.
Year |
BSF |
MF |
PF |
|
hatot
|
n |
n/DW classes |
Sum nactDW
|
nactDW/ha |
hatot
|
nactDW
|
n/DW classes |
Sum nactDW
|
nactDW/ha |
hatot
|
n |
n/DW classes |
Sum nactDW
|
nactDW /ha |
|
1 |
2 |
3 |
1 |
2 |
3 |
1 |
2 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a) MG (nact) |
2016 |
128 |
2110 |
57 |
495 |
410 |
962 |
8 |
7.4 |
89 |
0 |
9 |
14 |
23 |
3 |
13.6 |
93 |
1 |
29 |
15 |
45 |
3 |
2023 |
131 |
2213 |
63 |
436 |
454 |
953 |
7 |
6.8 |
121 |
2 |
12 |
17 |
31 |
5 |
11.2 |
179 |
4 |
64 |
27 |
95 |
8 |
b) MG (nactR) |
2023 |
131 |
1336 |
31 |
220 |
232 |
483 |
4 |
6.8 |
70 |
0 |
5 |
11 |
16 |
2 |
11.2 |
51 |
0 |
13 |
9 |
22 |
2 |
c) FB (nact) |
2019 |
167 |
2221 |
40 |
880 |
447 |
1367 |
8 |
– |
– |
– |
– |
– |
– |
– |
33 |
386 |
5 |
181 |
15 |
201 |
6 |
2023 |
167 |
2332 |
55 |
653 |
426 |
1134 |
7 |
– |
– |
– |
– |
– |
– |
– |
36 |
431 |
4 |
141 |
24 |
169 |
5 |
d) FB (nactR) |
2023 |
167 |
1446 |
34 |
412 |
268 |
714 |
4 |
– |
– |
– |
– |
– |
– |
– |
36 |
266 |
3 |
96 |
10 |
109 |
3 |
Table 4.
Descriptive statistics of number of active nests (nact) with WpC, number (n) of woodpecker cavities (WpC) in active nests for (a) MG and (b) FB for the 2016, 2019 and 2023 inventories. – = not present.
Table 4.
Descriptive statistics of number of active nests (nact) with WpC, number (n) of woodpecker cavities (WpC) in active nests for (a) MG and (b) FB for the 2016, 2019 and 2023 inventories. – = not present.
Year |
Study site |
Mapped nests (nact)
|
Nests (nact) with WpC (n) |
Sum nact
|
Numbers of WpC (n) in nact
|
Sum WpC (n) |
BSF |
MF |
PF |
BSF |
MF |
PF |
2016 |
MG |
2292 |
291 |
16 |
14 |
321 |
743 |
32 |
36 |
811 |
2023 |
MG |
2513 |
338 |
31 |
24 |
393 |
921 |
37 |
55 |
1049 |
2019 |
FB |
2607 |
272 |
– |
36 |
308 |
849 |
– |
120 |
969 |
2023 |
FB |
2763 |
331 |
– |
68 |
399 |
908 |
– |
246 |
1154 |