Conservation initiatives like the Natura 2000, the largest and the oldest programme in the European Union, aim to preserve local biodiversity through extensive land management practices [Evans 2012]. Wet meadows and calcareous grasslands of northwestern Europe are considered an important habitat for many plants and animals. Their biological diversity is renowned, including a variety of species listed under the habitat directives [Van Swaay 2002; Pe’er et al. 2014].
In Central Europe, various protected insects, birds and plants rely on such management, including mowing and cattle grazing. By employing low-intensity management on grasslands, heathlands and peatlands, threatened species like the butterflies,
The serious threat for biodiversity are invasive alien plants. Giant hogweed
Another particularly troublesome invasive alien plants are the goldenrods (
The aim of this study was to assess the impact of invasive plants – alien goldenrods on the population size of
Population sizes with 95% confidence intervals of
M | 56 | 39 | 80 | 553 | ||
F | 32 | 15 | 67 | |||
M | 62 | 34 | 113 | 2193 | ||
F | 63 | 35 | 114 | |||
J3A | 0.062 | M | 24 | 12 | 45 | 871 |
F | 30 | 18 | 51 | |||
J4A | 0.109 | M | 31 | 19 | 51 | 532 |
F | 27 | 15 | 48 | |||
M | 208 | 119 | 363 | 918 | ||
F | 171 | 84 | 348 | |||
J6A | 0.072 | M | 41 | 33 | 50 | 1306 |
F | 53 | 46 | 63 | |||
M | 33 | 22 | 49 | 293 | ||
F | 16 | 8 | 35 | |||
J2B | 0.204 | M | x | x | 123 | |
F | x | x | ||||
M | 61 | 23 | 165 | 453 | ||
F | 63 | 24 | 166 | |||
J4B | 0.315 | M | x | x | 35 | |
F | x | x | ||||
M | 57 | 40 | 81 | 141 | ||
F | 58 | 41 | 82 |
In cases where it was not possible to estimate the population size, the numbers of captured individuals were italicized. “Density” denotes relative butterfly density (abundance) expressed as number of individuals per 1 ha of area.
Population sizes with 95% confidence intervals of
M | 213 | 133 | 343 | 2547 | ||
F | 192 | 113 | 326 | |||
M | 114 | 61 | 212 | 3930 | ||
F | 110 | 58 | 210 | |||
J3A | 0.062 | M | 37 | x | x | 1113 |
F | 32 | x | x | |||
J4A | 0.109 | M | 130 | 82 | 206 | 1716 |
F | 57 | 36 | 91 | |||
M | 433 | 367 | 510 | 2019 | ||
F | 401 | 336 | 479 | |||
J6A | 0.072 | M | 84 | 64 | 110 | 2375 |
F | 87 | 67 | 113 | |||
M | 134 | 87 | 206 | 1473 | ||
F | 112 | 67 | 188 | |||
J2B | 0.204 | M | 20 | 13 | 29 | 328 |
F | 47 | 33 | 66 | |||
M | 135 | 64 | 282 | 942 | ||
F | 123 | 56 | 267 | |||
J4B | 0.315 | M | x | x | 76 | |
F | x | x | ||||
M | 137 | 91 | 206 | 285 | ||
F | 95 | 63 | 143 |
In cases where it was not possible to estimate the population size, the numbers of captured individuals were italicized. “Density” denotes relative butterfly density (abundance) expressed as number of individuals per 1 ha of area.
Surveys took place at the time of peak of butterfly flight period and just afterwards between 21 July 2022 and 5 August 2022 between 09:00 and 17:00, on rain-free days, with temperatures no lower than 20 °C. In the JA complex, capturing of individuals was carried out for eight days, while in the JB complex trapping lasted for seven days, due to reduced butterfly numbers and the mowing of the J4B meadow (Fig. 1 c).
Butterflies were captured on each habitat patch for one hour. Individuals were captured with an entomological net, marked on the wing with a unique identifier and finally released. In addition, for each habitat patch, data were recorded on the relative cover of the host plant
Each habitat patch was mapped. The boundaries of the habitat patches were marked using a GPS receiver. A polygon layer was created to calculate the area of the patches and prepare overview maps.
Using the POPAN model in programme MARK 10.0 [White, Burnham 1999], abundances of local butterfly populations were estimated with 95% confidence intervals. Abundances for each species were calculated by sex. In some cases it was not possible to estimate population abundances or confidence intervals, due to very low probability of capture and very low numbers of captured butterfly individuals. In such cases, only numbers of individuals caught are presented (Table 1, Table 2, Figs. 2 a, 2 b).
For further analysis, the estimated butterfly abundances were converted to relative densities of individuals per 1 ha of area (Table 1, Table 2, Figs. 2 a, 2 b). Butterfly densities, separately for butterfly species, were correlated with goldenrod cover and the ratio of expansive plant cover (sum of goldenrod and other aforementioned plants) to host plant cover. In addition, the effect of goldenrod and host plant on the summed density of both butterfly species was tested using ANOVA in the Statistica 13 software. In this analysis, plant cover was included as a categorical factor: low cover (up to 10%) and high cover (above 10%).
The capturing of Heteroptera bugs was carried out in the area of the village of Krempna (Podkarpackie Voivodeship) on sites with
Insects were caught using a sweep net. A 50-m section consisted of 100 scoops. Each step was followed by one blow with the sweep net. A total of 100 scoop strokes were made on each section. The above methodology ensures that the results obtained are standardised, due to the equal sample size and identical effort on each section. Material from each transect was poured into insect bags and labelled with the transect ID and date of trapping. The collected material was placed in a freezer to kill the insects and store them until sorting.
Heteroptera bugs from a given trapping and transect were separated from other captured insects and plant fragments, then sorted into small containers and described. The insects were then labelled to species and the numbers of individuals per bug species were counted. Finally, the number of species and the number of individuals present on each transect section were summarised.
The effect of environment type – transects with giant hogweed vs. reference transects – on the number of species (Poisson model) and the number of individuals (negative binomial model) was tested. Next, permutation methods, based on the Bray-Curtis dissimilarity measure, were tested whether species assemblages differed between transects with the invasive plant and reference transects. The analyses were performed in R 4.1.3 [R Core Team 2022] using a vegan package [Oksanen et al. 2022].
The JA complex of patches (Fig. 1a, 1b) was characterised by higher cover of
Both butterfly species surveyed occurred in each of the habitat patches studied.
Populations of both species inhabiting the JA complex were more numerous than those living in the JB complex. The most abundant populations of both species occur sequentially in patches J5A, J1A and J2A in the JA complex and J3B, J1B and J5B in the JB complex. Patches J2A and J6A were characterised by high relative butterfly densities in the JA complex, and patches J1B and J3B in the JB complex (Table 1, Table 2).
There were trends of decreasing butterfly densities with increasing goldenrod cover and increasing the ratio of expansive plant to host plant cover (Figs. 2a, 2b). However, these trends did not reach statistical significance. For goldenrod and
A total of 60 species Heteroptera bugs were found. Forty-two species were found on sites with hogweed present, while 44 species were found on reference sites (Table 3). The mean number of individuals was more than twice as high at the reference sites compared to the mean number of individuals at the hogweed sites (Wald Stat. = 9.24, P = 0.002) (Fig. 4). A similar effect was observed for the mean number of species. However, this effect was at the limit of statistical significance (Wald Stat. = 3.39, P = 0.065) (Fig. 4). Species assemblages of Heteroptera bugs occurring in the reference sites differed significantly from those in the sites where invasive giant hogweed was growing (pseudo-F = 1.66, P = 0.046). Mean species abundances in both habitat types are shown in Table 3. Species are ranked from those that occurred most abundantly. Species with a cumulative weight of 0.8 are most responsible for the difference found. Species whose abundances were statistically different have been highlighted. Of the most abundant species, all those whose contribution was statistically significant reached higher abundances at the reference sites (Table 3). Other species that were statistically significantly different were also more abundant at the reference sites (Table 3). It should be noted that the group of reference sites included more species found only in this environment compared to the opposite situation – the number of species found only in sites with invasive giant hogweed (Table 3).
Mean difference in Heteroptera bug species assemblages along 50m transect sections between reference sites and sites with invasive hogweed. Species with significantly different mean abundances are highlighted.
22.1 | 6.8 | 0.22 | 0.153 | |
16.3 | 1.5 | 0.39 | 0.016 | |
5.0 | 2.7 | 0.48 | 0.85 | |
3.1 | 0.6 | 0.54 | 0.04 | |
1.0 | 2.5 | 0.59 | 0.995 | |
3.3 | 0.4 | 0.62 | 0.446 | |
0.8 | 0.9 | 0.65 | 0.931 | |
0.5 | 1.1 | 0.67 | 0.993 | |
0.9 | 1.0 | 0.69 | 0.861 | |
0.9 | 1.2 | 0.72 | 0.995 | |
2.8 | 0.1 | 0.74 | 0.045 | |
0.8 | 0.1 | 0.76 | 0.003 | |
0.9 | 0.2 | 0.77 | 0.024 | |
0.0 | 0.8 | 0.78 | 0.992 | |
0.9 | 0.0 | 0.80 | 0.002 | |
0.4 | 0.5 | 0.81 | 0.994 | |
0.1 | 0.4 | 0.82 | 0.965 | |
0.1 | 0.5 | 0.83 | 0.996 | |
0.0 | 0.5 | 0.84 | 0.955 | |
0.4 | 0.2 | 0.85 | 0.608 | |
0.3 | 0.5 | 0.86 | 0.994 | |
0.5 | 0.0 | 0.87 | 0.01 | |
0.3 | 0.5 | 0.88 | 1 | |
0.0 | 0.4 | 0.89 | 0.989 | |
0.3 | 0.0 | 0.90 | 0.004 | |
0.1 | 0.3 | 0.91 | 0.979 | |
0.1 | 0.2 | 0.91 | 0.945 | |
0.3 | 0.1 | 0.92 | 0.569 | |
0.0 | 0.3 | 0.93 | 0.986 | |
0.1 | 0.2 | 0.93 | 0.976 | |
0.1 | 0.1 | 0.94 | 0.839 | |
0.3 | 0.0 | 0.94 | 0.01 | |
0.1 | 0.0 | 0.95 | 0.006 | |
0.3 | 0.0 | 0.95 | 0.016 | |
0.1 | 0.0 | 0.95 | 0.006 | |
0.1 | 0.0 | 0.96 | 0.016 | |
0.1 | 0.0 | 0.96 | 0.016 | |
0.1 | 0.0 | 0.96 | 0.016 | |
0.0 | 0.2 | 0.97 | 0.988 | |
0.1 | 0.1 | 0.97 | 0.962 | |
0.1 | 0.0 | 0.97 | 0.01 | |
0.1 | 0.0 | 0.97 | 0.018 | |
0.1 | 0.0 | 0.98 | 0.017 | |
0.0 | 0.1 | 0.98 | 0.978 | |
0.0 | 0.1 | 0.98 | 0.978 | |
0.4 | 0.0 | 0.98 | 0.031 | |
0.0 | 0.1 | 0.99 | 0.98 | |
0.0 | 0.1 | 0.99 | 0.983 | |
0.0 | 0.1 | 0.99 | 0.988 | |
0.0 | 0.1 | 0.99 | 0.988 | |
0.0 | 0.1 | 0.99 | 0.987 | |
0.0 | 0.1 | 0.99 | 0.987 | |
0.0 | 0.1 | 0.99 | 0.987 | |
0.0 | 0.1 | 1.00 | 0.987 | |
0.0 | 0.1 | 1.00 | 0.991 | |
0.1 | 0.0 | 1.00 | 0.029 | |
0.1 | 0.0 | 1.00 | 0.031 | |
0.1 | 0.0 | 1.00 | 0.031 | |
0.1 | 0.0 | 1.00 | 0.031 |
The study showed negative impacts of invasive plants on insects: two umbrella species, the butterflies
The habitat of the surveyed butterflies is in the Molinion meadows [Sielezniew, Dziekańska 2010]. These butterflies require two key resources in their life cycle: the host plant bloodroot and the host ants of the genus
For the inventoried Natura 2000 site Łąki w Jaworznie (Meadows in Jaworzno), the most important threats are listed in order (according to the Order of the Regional Director for Environmental Protection): abandonment/non-mowing, alien invasive species, problematic native species, land filling, land reclamation and drainage, and change in vegetation composition (succession). The threats listed continue to exist, and the state of habitat conservation and its impact on current butterfly population numbers is a result of land use regime, insufficient or misfunctioning conservation measures or even a lack of such measures from the establishment of the Natura 2000 site to the present.
The higher abundance of
Part 4 of the site with the JA complex is the best preserved in terms of habitat, and the most abundant local butterfly populations were recorded there. These patches have not been mown for many years and are overgrown with goldenrod and shrubs. There is also a succession towards vegetation from the association Filipendulion [Matuszkiewicz 2008]. The largest, central fragment of part 4 is regularly mown, however, not at the recommended time, but at a time convenient for the plot owners. With the current management – centre of the site mowed regularly but at the wrong time, and overgrowing edges of the Natura 2000 site – butterfly populations may still persist for a long time, even in a 20-year perspective.
In order to maintain abundant butterfly populations in the long-term perspective, parts of the area that have not been mown so far should be mown in rotation, every other year [Sielezniew, Dziekańska 2010; Kalarus et al. 2013]. This frequency of mowing is fully sufficient, as goldenrod is present sparsely in part 4. In addition, the negative impact of goldenrod on butterflies at the initial stage of the invasion can be partly compensated by the high cover of the host plant. However, the time of mowing is very important. Much of the area in the centre of part 4 is mown in June or early July before the flight period of the butterflies, thus sterilising the habitat and making it inaccessible to the butterflies. Butterflies of the genus
In parts 1 to 3 of the Natura 200 site, the problem of goldenrod invasion and succession is much greater than in part 4, which has had a negative impact on the abundance of butterflies. In addition, the problem of land drainage is also visible – in parts 1 and 2. In the JB complex, less abundant local butterfly populations were found with critically low abundance in patch J2B, where goldenrod cover reached 70%. In parts 1–3, mowing was mostly abandoned, leading to a situation where goldenrod coverage in part 3 reaches up to 100% and no butterflies could be found there. In the area where goldenrod cover is higher than 60%, mowing should be carried out at least twice a year at times before goldenrod flowering and later when regrowth of the cut invasive plant occurs [Kalarus et al. 2013]. Where goldenrod cover is lower patches of habitat should be mown once a year at the correct autumn date.
For the smallest part of the Natura 2000 site, part 2, there is an additional threat not previously considered. Part 2 is threatened by succession, invasion of goldenrod and deterioration of habitat conditions. On the other hand, in its buffer zone there are regularly mown, large patches of habitat that attract butterflies – for example, patches J4B and J5B. When habitat condition deteriorates in part 2 of the Natura 2000 site, butterflies disperse in search of the resources they need to survive [Nowicki, Vrabec 2011; Nowicki et al. 2014] and end up in the habitat patches in the buffer zone. There, the butterflies lay their eggs on the flowering host plant. However, their reproductive efforts are wasted, and a new generation of butterflies will not emerge due to the mowing of these patches when the butterflies fly, leading to the death of the eggs laid and the young larvae already feeding on the host plant. Patches such as meadows J4B and J5B are an ecological trap for the investigated butterflies. In order to ensure the survival of the butterflies in the long term, the habitat within the Natura 2000 site should be maintained in the best possible condition. The mowing schedule must be unconditionally aligned with conservation goals, even if other objectives appear profitable [cf. Szyszko-Podgórska 2019; Wróbel et al. 2022]. One should also remember about proper drainage, which will prevent the drying out of the butterfly habitat [cf. Syed, Shuqi 2022].
It is worth noting that even small patches of poorer quality habitat are important for butterflies. High relative densities of the butterflies have been observed in such patches. Such habitat patches can play an important role in the butterfly migration process and act as stepping-stones, thereby improving landscape connectivity [Kalarus, Nowicki 2015].
Invasive hogweed is a major problem in the Magurski National Park. The hogweed appears on meadows, forest paths, forest edges, roadsides, along drainage ditches and even in the gardens of private properties and often reaches large coverages. The study was carried out in the first year that active hogweed removal treatments were introduced in the study area. Despite this fact, a significant negative impact of hogweed on the abundance and number of Heteroptera bug species was found. Hogweed, also affects the bugs at the assemblage level. The most abundant species in the assemblages occurred less abundantly on sites with hogweed. Similarly, among the less abundant species, there were more species (especially those with a statistically significant proportion) that occurred only on reference sites.
The flowers of invasive hogweeds were characterised as simple habitus and not specialised for any particular type of pollination. They relied on insects for pollination and were attractive to a diverse group of generalist pollinators. A wide range of insects, including Hymenoptera (bees, wasps, ants), Diptera (flies), Coleoptera (beetles), and Hemiptera (bugs), visited these flowers for nectar and pollen [Grace, Nelson 1981; Hansen et al. 2007; Nielsen et al. 2008]. Previous studies have shown that hogweeds are used mostly by polyphagous insect species or species specific for
Hogweed appears to exclude stenotopic species and those with specialist habitat requirements. The negative impact of hogweed can be considered a strong impact, as both abundant and less abundant species of Heteroptera bugs were found in lower abundance on sites with hogweed compared to reference sites. In addition, the negative impact of hogweed has emerged, despite the removal treatments carried out.
Removal treatments have only recently been introduced, and the reduction of hogweed and its negative effects on biodiversity may take several years to become apparent. Efforts should be made to increase the abundance and diversity of Heteroptera bugs, on sites where hogweed was growing.
A major constraint that will affect the success of hogweed removal efforts is the fact that not all landowners and municipalities will undertake simultaneous hogweed removal. The hogweed removed in one location may be replaced by hogweed that arrives/spreads from another nearby location where no removal action has been carried out. In order to reduce the negative effects of hogweed, consistent, long-term management to remove the invasive plant is essential.
Giant hogweed is known for its substantial size, reaching heights of up to 4–5 m, featuring extensive clusters of leaves measuring 2–3 m in length, and exhibiting stable flowering shoots that resemble woody plants [Panasenko 2017]. Additionally, it produces large inflorescences. Its invasion has resulted in a complex system that encompasses various aspects of ecosystems, including environment, agriculture, forestry, land use, hydrosphere, soil system, global change ecology, biodiversity, management, and conservation. The understanding of its invasion requires a multidisciplinary research system that attract specialists from diverse fields who have conducted extensive studies on these plants.
In the JA and JB grassland complexes both species occur, and
Part 4 of the Natura 2000 site has the best habitat condition and the most abundant butterfly populations.
Parts 1–3 of the site are seriously threatened by succession and invasion of goldenrods.
In the smallest part 2, its buffer zone attracts butterflies and provides an ecological trap for them, which negatively affects butterfly populations.
The negative effect of goldenrod on butterflies at the initial stage of invasion may be partially offset by the high cover of the butterfly host plant.
Meadows providing habitat for the butterflies must be mown in the autumn from mid-September to the end of November.
Mowing meadows during butterfly flight, or just before butterfly flight, is not allowed and is contrary to conservation objectives.
In areas where goldenrod cover reaches more than 60%, more frequent mowing should be used – the first before this invasive plant flowers.
Small, poorer quality patches of habitat are also important for butterfly conservation, as they support butterfly dispersion and improve habitat connectivity.
Invasive hogweed has a negative impact on the abundance and assemblage structure of the assemblage of Heteroptera bugs.
The negative impact of hogweed has been detected despite the application of the first eradication treatments.
The effects of hogweed control may only become apparent after several years.
Success in eradicating hogweed may be limited due to the spread of hogweed from locations where removal is not applied.