Spatiotemporal pattern of wild boar rooting in a Central European dry oak forest

Wild boar (Sus scrofa) density has significantly increased worldwide. In Europe, oak acorns are basic diet items for the species. However, regeneration of temperate oak forests has become excessively low. Thus, better understanding of the patterns and dynamics of wild boar rooting and the spatiotemporal relationship between the acorn density and the rooting has special importance. In our study, the acorn density, the presence and intensity of the rooting were measured monthly in a 28-ha oak forest stand between 2016 October and 2019 April. Study site was divided to 400 m2 grid cells to investigate the spatiotemporal patterns of the rooting. Comparing the acorn densities, we stated that in the first 2 years, but not in the third one, the oaks were masting. During the acorn-rich periods (October–April), significant inter-annual differences were found in the proportions of the disturbed area between the non-mast and the other two mast years. Throughout our observations, 9.12% of the cells stayed undisturbed. On average, in 12.59% of the cells new rooted patches appeared, and 15.6% of them was repeatedly disturbed. Rooting were mainly litter disturbances in 71.75% of cases. Throughout acorn fall periods (September–November), the most intensively rooted sites were the poorest in acorns, while the unrooted sites were the richest. Our results demonstrate that wild boars have a great effect on acorn density and on the whole forest by rooting almost the entire ground surface at least once, but also the masting of the oaks has a crucial role in their rooting dynamics.


Introduction
Nowadays, the regeneration of the temperate oak (Quercus spp.) forests has become excessively weak on the Northern Hemisphere, often reflected by the missing or scarce regeneration layer (Brudvig 2008;Katona et al. 2013;Leonardsson et al. 2015). This fact raises a number of serious economic and conservation issues (Brudvig 2008;Schelhaas et al. 2012;de Rigo et al. 2016;Eaton et al. 2016).
The outstanding biological importance of the oak forests is unquestionable; they act as important biodiversity hot spots and provide habitat for a great number of species of insects, fungi, bats and other organisms including several threatened species (Bobiec et al. 2011;Leonardsson et al. 2015;de Rigo et al. 2016;Eaton et al. 2016).
During the whole regeneration cycle of the woody plant species, seedlings represent one of the most vulnerable phases. The mortality rate of the seedlings can be particularly high, and they are susceptible to herbivore browsing and other biotic and abiotic circumstances (Gómez and Hódar 2008;Perea and Gil 2014).
The appropriate quantity and quality of seedlings, however, largely depends on the processes in the previous stages of the regeneration. The acorn production and the germination of the seeds can have a critical influence on the oak regeneration (Gómez and Hódar 2008;Perea and Gil 2014;Massei et al. 2015;Kamler et al. 2016). Additionally, seed predators also have a significant role in decelerating the regeneration process. The most typical forest seed predators include different birds, rodents and ungulates, especially Communicated by Claus Bässler.
The survival of the oak acorns is driven by many factors after falling from the tree. Apart from predation, they can often be ruined by infections and other abiotic effects (van Ginkel et al. 2013). But the primer determinant of their survival chance is the natural density of the oak acorns, which is enhanced by periodically occurring mast years (Sork 1993;Bisi et al. 2018;Bogdziewicz et al. 2018). According to the studies by Cutini et al. (2013) and Nováková and Štípek (2011), there is a strong positive correlation between the reproductive success of wild boar and acorn production. The tree can defend against its seed predators by pulsing the availability of resources-the acorns-creating lagged responses in the successive trophic levels. This can allow a better survival for the acorns in a mast year along with a decreased seed predator population (Haas and Heske 2005;Bisi et al. 2018).
In the recent decades, the number of wild boars has significantly increased worldwide. Simultaneously, the number of conflicts related to the species has also escalated (Massei and Genov 2004;Barrios-Garcia and Ballari 2012;Kamler et al. 2016;Bongi et al. 2017). Wild boar is an opportunistic omnivore species, although 90% of its diet consists of plant materials (Massei and Genov 2004;Barrios-Garcia and Ballari 2012;Bongi et al. 2017;Katona and Heltai 2018).
The most important food components are the large fruits and seeds of high energy content. In the temperate zone of Europe, the acorn of oak and beech can constitute as such food items (Groot Bruinderink and Hazebroek 1996;Massei and Genov 2004;Ballari and Barrios-García 2014;Kamler et al. 2016;Bisi et al. 2018). The acorn means the main diet component from autumn to spring, and it is not uncommon that the stomach content of wild boars consists of acorns around 90% in mast years (Groot Bruinderink et al. 1994;Bongi et al. 2017). Other studies also showed that acorn is the preferred food source for wild boars even if other food sources were available, e.g. anthropogenic or agricultural ones (Katona and Heltai 2018).
Acorn and other below-ground materials can be found and ingested through the rooting of the litter and the soil surface; therefore, rooting is the obvious sign of the feeding behaviour of wild boar (Welander 2000). Wild boars can consume a significant amount of oak affecting the regeneration of oak species, but by their rooting they also drive the changes in the composition of the above-ground herbaceous layer (Horčičková et al. 2019). Eventually, they can impede the regeneration of the forest and alter succession process and change their environment (Groot Bruinderink and Hazebroek 1996;Massei and Genov 2004;Gómez and Hódar 2008;Barrios-Garcia and Ballari 2012;van Ginkel et al. 2013;Kamler et al. 2016;Bongi et al. 2017;Ramirez et al. 2018). Based on this, wild boar can be considered as an ecosystem engineer, being able to profoundly change directly or indirectly the characteristics of the habitat and the availability of resources (Barrios-Garcia and Ballari 2012; Mráz and Katona 2014).
Some previous studies have already revealed the patterns of the rooting and the environmental factors driving them (Groot Bruinderink et al. 1994;Welander 2000;Bisi et al. 2018). The temporal patterns are relatively better described (Dardaillon 1986;Groot Bruinderink and Hazebroek 1996;Ballari and Barrios-García 2014;Kamler et al. 2016). However, we have much less knowledge about the spatial distribution patterns of the rooting and its dynamics along the year, particularly in temperate oak forests in periods with different intensities of mast-fall. Therefore, our questions were: 1. What is the proportion of the rooted areas of an oak forest in different periods of the year of different acorn availability? 2. How does the distribution of the rooted areas change from month to month in mast and non-mast years? 3. Is there any spatial relationship between the rooting intensity and acorn density?

Study area
The study area is located in Hungary, on the north of the lake Balaton, in Veszprém county, in the region of Zánka village (46°52′32.8′′ N 17°38′57.4′′ E). The area is in the dryer parts of the humid continental climatic region, with an average annual temperature between 10.2 and 10.5 °C. The annual average rainfall is between 600 and 640 mm. The bedrock of the area is dominated by new red sandstone; meanwhile, the soil can mainly be characterized by red clay rendzina.
Our investigation was carried out in a private dry Turkeysessile oak forest (Fraxino orno-Quercetum cerridis) managed by Hercegerdő Forest Corporation. The study area was a 28-ha forest unit surrounded directly by similar oak forests. Agricultural fields, dominated by vineyards, can be found in less than 500 m afar in the north-northeast neighbourhood. The forest unit had a quite even Turkey-sessile oak abundance (cca. 48% and 52% among oak trees, respectively). The forest stand is fairly mixed, and cca. 30% of the tree population is subdominant tree species, such as flowering ash (Fraxinus ornus), field maple (Acer campestre), checker tree (Sorbus torminalis) and common hornbeam (Carpinus betulus). The forest is managed by single tree selection and group selection patch cuts since 1994. The main product of the area is firewood. Mitigation of the oak regeneration due to acorn predation by wild boar is considered as a serious problem by the local forest managers, reflected in the fenced small gaps, established by previous timber utilization, to enhance the development of acorns and defend the emerging saplings.

Field data collection
After we chose our study area, we recorded several border points with a GPS device (GARMIN eTrex 20). The data for the monitoring have been collected once in about every 35 days between 2016 October and 2019 April. We designated 12 transects parallel with each other and the southern border of the area. All the transects were about 2 m wide, 500 to 1400 m long and 20 m from each other. On these transects, one sampling point per every 100 m was appointed; altogether, we obtained 118 sampling points (except in 2016 October when there were 122 and in 2017 September when there were 71, due to bad weather conditions) (Fig. 1).
On all these sampling points, we measured the density of acorns by oak species in a 1 m × 1 m quadrant. During this examination, we moved over the litter at the sampling points to find and count all the acorns. Then we arranged the litter back to decrease our disturbing effect. Along the full length of the transects, we also recorded by the GPS device the locations of all patches rooted by wild boar. When our transects crossed the rooted areas, we recorded an entry and an exit point. We considered a patch as rooted by wild boar when the disturbed site heavily differed from the undisturbed surroundings. Wild boar rooting often means the disturbance of the top of the soil layer. Thereby litter and soil can mingle. The size of the rooted patch usually varies between 1 and 100 m 2 , and the disturbed soil layer can be up to 20 to 40 cm deep.
During our data collection, we categorized these disturbed areas by the intensity of the rooting into three main categories: • Litter disturbances (further on LD): Only the litter was slightly tossed and moved over. • Upper soil disturbances (further on USD): More intensively disturbed patches, where both the litter and soil were rooted and slightly agitated. Only the top 5 cm of the soil was disturbed. • Deeper soil disturbances (further on DSD): When the rooting was relatively intensive and disturbed the soil deeper than 5 cm.

Data analysis
The data have been processed and analysed in Microsoft Office Excel 365, StatSoft Statistica 12 and in QGIS 2.18 Las Palmas (in HD72/EOV (EPSG: 23700) coordinate reference system).
We draw the borders of our study area from the collected points, and then we converted it to a polygon. A new layer has been added which was divided to 20 m by 20-400 m 2 cells. We kept the cells which were on our polygon-therefore within the borderline of the study area-that gave us 702 cells of 400 m 2 , altogether 28.08 ha. The entry and exit points have been sorted by the intensity of the rooting and then converted to litter, upper soil or deeper soil disturbance lines by linking the points. We analysed how these newly drawn lines distributed over the study area by marking the cells if there was any rooting in them. Although, it is important to mention, that it is some kind of overestimation of the rooted surface as the recognition of rooting along our transect did not mean that the whole-400 m 2 -sample unit, the full cell was disturbed. We were also able to separately mark the cells by the different kind of rooting types that they contained at a given time. In this way, we investigated the intensity of area use by wild boar population (its impact on the habitat) and not the home range use of single individuals.
We divided the study to different temporal sample units in two different ways. On the one hand, we appointed the acorn fall period, which was between 2016 October to November, 2017 September to November and 2018 September to November, when we expected the oak acorns to fall from the trees. On the other hand, we also appointed an acornrich and an acorn-weak period. The acorn-rich period was between October and April in every study year, when we expected that there are acorns available on the ground to the seed predators before germinating and becoming a seedling or ruined and disappearing. The acorn-weak period meant the rest of the study interval. We analysed the changes of the acorn density and wild boar rooting at the study site in the above-mentioned periods. We observed the inter-annual effects and determined the mast and non-mast years.

Statistics
We compared the acorn density, the proportion of the rooted cells, the proportion of rooting categories based on the temporal within-cell changes of disturbances between consecutive assessments-if the rooting in the cell appeared, remained, disappeared or the cell was repeatedly undisturbed-and the proportion of the classes of different rooting intensities between the acorn fall versus the non-fall and the acorn-rich versus the acorn-weak periods with Mann-Whitney U test. Kruskal-Wallis tests with Dunn's post hoc tests were applied to measure the significant differences among the values of all these above-mentioned variables among years or their sub-intervals. Our acorn density datasets showed non-normal distribution by Shapiro-Wilk W test (W = 0.19 to 0.67 and P < 0.0001 for all the measured periods) and in other cases typically formed by nominal and categorical variables. That was the reason to rely on nonparametric statistical tests during the comparisons.
We investigated how the proportions of the different kinds of earlier mentioned rooting intensity types changed during these time intervals and how this intensity interfered with the acorn density. In cases when different rooting types were revealed in the same cell (see as combined rooting categories later, e.g. LD + USD + DSD), we classified it to the highest intensity type based on what the cell contained. We categorized the acorn densities from different assessments separately; first we identified the sampling points with outlying high acorn density by the ROUT method (Motulsky and Brown 2006) (Q = 0.1%) in the case of every assessment. This gave us the upper acorn density category. Then if it was possible, we split up the rest to two other groups (middle and lower categories) with Jenks optimization method (Jenks 1967). To analyse the distribution of the different rooting intensity types along with the acorn density categories, Pearson or maximum likelihood ratio Chi square was used.

Temporal pattern of the oak acorn density
The average acorn density throughout the study was 2.75 ± 7.81 pcs./m 2 . However, if broken down to different periods, we got a more diverse picture. Figure 2 shows that our data were pretty skewed, which means in a large number of cases the oak acorn density was fairly low or even there was no acorn at all at the sampling points.
But if we examine the inter-annual effects on the acorn density, the different acorn fall periods varied. The medians were 2.5 pcs./m 2 (IQR:7) in 2016, 5 pcs./m 2 (IQR:13) in 2017 and 0 pcs./m 2 (IQR:0) in 2018 autumn. We found that the acorn density significantly differed among the annual acorn fall periods (KW = 215.9, p < 0.001); in 2018 autumn, the acorn density was lower than in the two previous acorn fall periods (Dunn's post hoc test: p < 0.001). Through the study, the highest acorn density median was registered in 2017 November (6.00 pcs./m 2 , IQR: 11.75). Therefore, we can assume that in 2016 and 2017 the oaks were masting, but not in 2018.
During the three investigated acorn-rich periods, the average acorn density was 2.96 ± 7.27 pcs./m 2 (Mdn: 0.00, IQR:3) and was higher (U = 431,434.5, p < 0.001) by contrast in the acorn-weak periods, when it was 1.88 ± 9.71 pcs./ m 2 (Mdn: 0.00, IQR:0). If we observe the annual acorn-rich periods, the medians of the acorn density were 1.00 pcs./m 2 (IQR:4) in the first, 2.00 pcs./m 2 (IQR:5) in the second and 0.00 pcs./m 2 (IQR:0) in the third year. The patterns were quite similar to the acorn fall periods, and the density was significantly lower in the third acorn-rich period than in the previous two (KW = 471.71, p < 0.001; Dunn's post hoc test: p < 0.001).

Total wild boar rooting disturbance of the study area
Through the spell of the study, the average proportion of the disturbed cells was 28.11 ± 22.01%. The maximum proportion of disturbed area was 80.31% and was registered in 2018 February, while the minimum was 0% and was observed in 2018 December (Fig. 3).
But during the acorn-rich periods, the average proportions of the disturbed area were the following: chronologically: 34.5 ± 14.12% (Mdn: 31.53%), 52.57 ± 21.94% (Mdn: 54.21%) and 5.42 ± 4.11% (Mdn: 5.28%). This inter-annual difference was significant between the third and the other two acorn-rich periods (KW = 13.46, p = 0.001; Dunn post hoc test: p < 0.05). The inter-annual changes of rooting were similar to the yearly pattern of the acorn density.

The spatiotemporal dynamics of the wild boar rooting
Throughout our observations 9.12% of the cells stayed undisturbed all along while we monitored the study area. Meanwhile, the highest disturbance frequency of a cell was 79.17%. The average disturbance frequency was 30.15 ± 17.75%. During our study, on average in 12.59 ± 9.58% (Mdn: 10.54%) of the cells new rooted patches appeared from one sampling to the next one ("appeared"), 15.6 ± 17.21% (Mdn: 7.83%) was repeatedly "disturbed", from 13.22 ± 11.05% (Mdn: 9.97%) of the cells the previously rooted patches "disappeared" and 58.59 ± 25.13% (Mdn: 54.84%) remained "undisturbed" (Fig. 4). The proportion of the "undisturbed" cells was significantly higher than that of the others (KW = 44.62, p < 0.001; Dunn's post hoc test: p < 0.001).
During the acorn-rich period of the non-mast year (2018-2019), there were significantly fewer "appeared" Fig. 2 Measured acorn densities at sampling points during the study interval. Acorn symbols indicate the acorn fall periods, while snowflake symbols indicate when the study site was covered by snow during the sampling. The breaks on the horizontal axis refer to gaps in data collection 1 3 Fig. 3 Proportion of the disturbed and undisturbed cells by rooting during the study interval. Acorn symbols indicate the acorn fall periods, while snowflake symbols indicate when the study site was cov-ered by snow during the sampling. The breaks on the horizontal axis refer to gaps in data collection Fig. 4 Spatiotemporal dynamics of rooting disturbance on the study site. Acorn symbols indicate the acorn fall periods, while snowflake symbols indicate when the study site was covered by snow during the sampling. The breaks on the horizontal axis refer to gaps in data collection cells than in the acorn-rich periods of the two mast years (KW = 8.27, p < 0.05; Dunn's post hoc test p < 0.05). Also, in the 2018-2019 acorn-rich period there were significantly fewer "disturbed" cells than in 2017-2018 (KW = 13.04, p = 0.0015; Dunn's post hoc test: p = 0.002) and there were more "undisturbed" cells than in the other acorn-rich periods (KW = 12.62, p = 0.0018; Dunn's post hoc test: p < 0.05).
Similarly to this, if we compare the acorn-rich and acornweak periods of the first two masting years, the proportion of the "appeared" cells was higher in the acorn-rich periods (U = 16.5, p < 0.05), while the proportion of the "undisturbed" cells was higher in the acorn-weak periods (U = 14.5, p < 0.05).

Temporal pattern of disturbance intensity in the wild boar rooted cells
The average proportions of different kinds of the main disturbance types throughout the study are as follows: LD: 71.75 ± 33.04%, USD: 12.84 ± 14.68%, DSD: 10.72 ± 18.46% and between 0.2 and 2.08% for the cells that were affected jointly by different rooting intensities, respectively (Fig. 5). The proportion of LD was significantly higher than that of all the other intensity types (KW = 85.04, p < 0.001; Dunn's post hoc test: p < 0.001). Moreover, the proportion of USD was also statistically higher than the values of three combined categories (LD + DSD; USD + DSD; LD + USD + DSD) (Dunn's post hoc test p < 0.05).

Fig. 5
Proportion of rooting intensity types during the study interval. Acorn symbols indicate the acorn fall periods, while snowflake symbols indicate when the study site was covered by snow during the sampling. The breaks on the horizontal axis refer to gaps in data collection The proportion of different rooting intensity types did not show any significant differences in inter-annual comparisons of the acorn-rich periods (Kruskal-Wallis tests: p > 0.05).

The interplay between the acorn density and the rooting intensity
The result of the acorn density categorization is shown in Table 1. In the first two mast periods of the study, we were able to categorize the acorn density into three categories, while in the last non-mast period we were only able to categorize it into two classes (except 2018 October).
Both, in the acorn fall and non-fall period there was a significant relationship between the acorn density and the intensity of the rooting (fall period: M-L X 2 = 22.49, df = 6, p < 0.001; non-fall period: Pearson X 2 = 92.61, df = 6, p < 0.001), meaning that the acorn density varied along the differently rooted areas.
During the acorn fall period, when there was no rooting in a cell ("undisturbed"), the proportion of the upper acorn density category was the highest (15.93%). But towards the more intensive rooting categories this proportion decreased and even missed from most of deeper soil disturbed cells. The proportion of the middle acorn density was the highest when there was litter disturbance (17.42%).
But if we examine the acorn fall periods separately, significant relationship was found only in the last, non-mast acorn fall period between the acorn density and the intensity of the rooting (M-L X 2 = 14.47, df = 6, p = 0.025). Moreover, cells with deeper soil disturbances were also only found in this period, but all with lower acorn density. Higher acorn density categories only occurred in the undisturbed and litter disturbed cells (upper acorn density: 20.19% and 25%; middle acorn density: 3.11% and 0%, respectively). On the contrary, in the first two-masting-acorn fall periods no Table 1 Acorn density categories Italics indicate the acorn fall periods, while snowflake symbols indicate when the study site was covered by snow during the sampling a Determined by ROUT method b Determined by Jenks optimization method; GVF: "goodness of variance fit" in Jenks optimization, 0 means worst fit, 1 means perfect fit; ND: "no data"

GVF b
Lower acorn density class b In the non-fall period, despite the significant difference among the cells of different rooting intensities, rooting impact did not cause so obvious decrease in the lower acorn supply (Fig. 6a, b).

Discussion
This study indicates that wild boar is a powerful ecosystem engineer species which can strongly influence its habitat due to its most peculiar attribute the rooting behaviour (Genov 1981;Massei and Genov 2004;Barrios-Garcia and Ballari 2012;Sandom et al. 2013;Genov et al. 2017;Sondej and Kwiatkowska-Falińska 2017). Approximately 90% of the study area was rooted at least once-which is much more than what Welander (2000) found in another deciduous forest-and only about 10% of the forest was exempt from rooting. Our study showed that wild boar rooting affected on average about 30% of the forest at a given time consisting of 15% of newly (within 1 month) created surface disturbances. Previous studies have already shown that wild boar can affect the soil nutrient cycle (Singer et al. 1984;Wirthner et al. 2012;Don et al. 2019) or the cover and the diversity of the forest herbaceous layer (Brunet et al. 2016;Sondej and Kwiatkowska-Falińska 2017) and of course the regeneration of the oak species (Groot Bruinderink and Hazebroek 1996;Gómez and Hódar 2008;van Ginkel et al. 2013;Perea and Gil 2014;Kamler et al. 2016;Zeman et al. 2016;Bongi et al. 2017).
We expected that the wild boars would root most actively during the acorn fall, even so their activity showed some lagging and they only became really active towards the end of the fall period. Moreover, their most active phases were usually during the winter and spring months-quite similar to the findings by Sandom et al. (2013)-but only if there was an adequate mast before. Though, the only time when we did not find any rooting was in winter (2018 December), these could have been due to the fresh snow cover and the lack of mast. Our studies showed that the rooting behaviour of the wild boar is fundamentally driven by the acorn density. During the first two autumns in our study, the oaks were masting, unlike in the last one. This phenomenon essentially determined the characteristics of the study area for the next few months from the perspective of available food source and therefore wild boar rooting. It also presumably affects the population density of the wild boars in the area, because oak acorns are the main nutrition source for the species starting from autumn, affecting their home range size and population dynamics (Bieber and Ruf 2005;Cutini et al. 2013;Morelle et al. 2015;Zeman et al. 2016;Bisi et al. 2018). It was also represented well in this study by the lower range of the rooted ground surface in the non-mast year, reflecting the lower degree of use of this habitat.
The intensity of the rooting also showed similar patterns; wild boars tended to root mainly the litter layer, and then rooting became more intensive during the summer months when the acorns hardly occur. The main target of this more intensive rooting was therefore probably not the acorns but more likely earthworms, plant roots or bulbs Fig. 6 Proportion of the acorn density categories in cells of different rooting intensities, a during the acorn fall period, b during the non-fall period (Welander 2000;Baubet et al. 2003Baubet et al. , 2004Massei and Genov 2004;Ballari and Barrios-García 2014).
Acorn density quickly decreased after the fall period. Therefore, we supposed a parallel decrease in the presence of rooting in the area. Even so we could not reveal any significant difference in the proportion of rooted cells between fall and non-fall periods. But we found that acorn density was lower in cells where rooting of different intensities was registered than in undisturbed cells. Areas with the most intensive rooting only had sampling points with lower acorn density, while sampling points with upper acorn density occurred in the largest proportion mainly in the undisturbed areas. The less intense litter disturbances can slightly diminish the acorn supply, but more intensive rooting will abolish the acorns. From the other side of this relationship, only in the third-nonmast-year we found deep rooting of the highest intensity, causing the entire depletion of acorns from these patches. All these suggest that wild boars likely show a less extent, but patchier rooting behaviour when searching for food through the acorn fall, and the proportion of acorns consumed by wild boar will be higher in the case of lower acorn density.
Wild boars are substantial seed predators (Groot Bruinderink and Hazebroek 1996; Cutini et al. 2013;van Ginkel et al. 2013;Kamler et al. 2016), and therefore they have a significant effect on the acorn density. Occasionally, they can even consume almost all of the produced acorns and therefore can impede the regeneration of the oaks (Cutini et al. 2013;van Ginkel et al. 2013;Zeman et al. 2016). Based on this study, it is still undecidable if the impact of the wild boars in the studied forest was sufficiently large to impede the oak regeneration completely. But we can assume that seed predators and principally wild boars had a great effect on acorn density and moreover the whole forest, by rooting almost the whole studied area at least once.
However, the effects of the wild boar rooting are not inherently negative for the forest; on the contrary, it can be a quite important function in an ecosystem. Rooting at an adequate level can provide an intermediate disturbance leading to an increase in biodiversity (Shea et al. 2004). Although the tolerable level of disturbance depends on the purpose of forest management. Even moderate wild boar rooting might be undesirable in commercial forests and can decelerate regeneration significantly, while in forests where the natural processes are more prioritized, a slower, patchier regeneration can also be sufficient. We believe that this methodology can help us to get a greater understanding of the effects and the dynamics of the wild boar rooting in temperate oak forests and to identify the favourable patterns of wild boar disturbances.

Conclusions
Understanding the spatiotemporal dynamics and effects of the wild boar rooting in oak forests is essential to manage the occurring conflicts in the future that can emerge with the rising wild boar populations. Our results suggest that wild boar can have a significant effect on oak regeneration, and the masting of the oaks has a crucial role in their rooting dynamics. We suggest the development of a predictive system for wild boar impacts (by using data on mast trees, acorn production and spatiotemporal characteristics of rooting site selection of wild boar) which can help the managers to identify when and where to expect large proportions of rooting and problems in the regeneration. This system can help intervene in a more efficient way, e.g. through proactive local population regulation of wild boar or fencing of sensitive regeneration patches.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.