Research Article |
Corresponding author: Zdeněk Faltýnek Fric ( zdfric@gmail.com ) Academic editor: Martin Wiemers
© 2024 Terezie Peškařová, Alois Pavlíčko, Tomáš Kuras, Zdeněk Faltýnek Fric, Martin Konvička.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Peškařová T, Pavlíčko A, Kuras T, Fric ZF, Konvička M (2024) Population status of the highly endangered Lycaena helle (Papilionoidea, Lycaenidae) in the Šumava Mts. two decades after establishment. Nota Lepidopterologica 47: 171-186. https://doi.org/10.3897/nl.47.126025
|
The process of translocation of endangered species represents a useful conservation tool, but subsequent monitoring of translocated populations is often neglected, although it may supply critical information for conservation work. The EU-protected glacial relic Lycaena helle (Lycaenidae) went extinct in the Czech Republic, where a few lowland populations historically existed, in the 1950s. In 2002, a mountain-dwelling population was established in Šumava Mts., among grasslands and fenlands surrounding the abandoned village of Nové Údolí, from 2/33 ♂♂/♀♀ originating in the Türnitz Alps, Austria. The transferred population was surveyed a decade later and again 21 years after the transfer. During the second (2023) survey, we visited 212 grassland patches to record the species’ presence and relative abundance, and to ascertain its habitat preferences. A third of the surveyed patches were occupied by L. helle, in a total area of occupied grasslands was 66.6 ha scattered over 4.5 km2, with approximately 2,400 adult individuals, recorded up to 4.2 km from the original release point. L. helle typically occupies wet tussocky grasslands with closed edges, minimum management and high cover of its larval host plant, the mid-to late-successional forb Bistorta officinalis. It avoids drier managed grasslands rich in flowering plants. The expansion rate from the release point (0.1 km/year during the first decade, 0.2 km/year over two decades) was slower than in two earlier L. helle transfers, targeting Morvan and Forez Mts., France (≈ 0.5 km/year in both cases). We attribute this to the relief of the release locality, a shallow depression surrounded by higher-elevated forested ridges, but we expect faster expansion once the continuous woodlands are crossed.
Owing to their aesthetic appeal, popularity, good life history knowledge and armies of enthusiasts, butterflies have always ranked high among animals intentionally translocated to novel localities (
The perception of butterfly translocations by experts also varies greatly. On the one hand, deliberate unauthorised translocations are discouraged because they interfere with species’ natural ranges and evolutionary processes (
For the Violet Cooper, Lycaena helle (Denis & Schiffermüller, 1775) (Lepidoptera: Lycaenidae), a species protected by the EU Habitats Directive (
The Eurosiberian range of L. helle encompasses the Hercynian mountains in Western Europe (Westerwald, Eifel, Ardennes, Vosges, Black Forest, Jura, Madeleine), Massif Central, Pyrenees and the Alps (a few scattered populations) (
Left: Aerial map of the section of Šumava National Park showing Czechia/Germany state border (narrow white line), the sections of the park surveyed for presence of the transferred butterfly Lycaena helle (bounded by narrow red lines), the original release site (blue dot), the most distant observation in 2023 (yellow dot), and all the occupied grassland patches (red dots). Right, top: Distribution map of L. helle in the Czech Republic, with blue dots standing for pre-1950 records, green dots for pre-1980 records, and red dots the introduced Šumava NP population. Right, middle, shows a section of Central Europe with sites of origin and release of the L. helle transfer. Right, bottom shows approximate distribution of the donor population in Türnitz Alps, Austria. All maps are from mapy.cz, licenced by Creative Commons.
Larvae develop on the the bistorts Bistorta officinalis Delarbre, or B. vivipara (L.) Delarbre in the extreme North. When used as the larval host, B. officinalis is also the main nectar source (
In Czech Republic, only lowland populations of L. helle existed historically, inhabiting fens along major rivers (
Here, we report results of a survey, carried out in 2023, i.e., two decades after the transfer, targeting distribution extent and habitat selection of the transferred population. We aim to answer two specific questions: 1) What habitat is selected by the population, does it differ from other populations in Western and Central Europe, and what implications emerge for habitat management? 2) What is the size of currently occupied area and how rapidly had the population expanded, if compared with the two previous transfers?
Lycaena helle is univoltine in mountains, with adults on wings in May–June, but forms two to three generations in lowlands (
The source of the transferred butterflies was a population inhabiting a mountain valley in Türnitz Alps, part of Northern limestone Alps, Austria (along local road L101 Schmelz – Ulreichberg; 47.84N, 15.41E, elevation 900 m; “Alps: Mariazell” in
The release site (Nové Údolí, 48.82N, 13.80E, 860 m elevation; 160 km aerial distance from the source locality; Figs
Habitats of the newly established population of Lycaena helle in the Šumava Mts. Left: A grassland patch in most suitable stage, with high representation of blooming host plant Bistorta officinalis. Right, top: A grassland patch in more advanced stage of succession, with progressing encroachment by Salix spp. shrubs. Right, bottom: Patch at the banks of Studená Vltava River with Carex brisoides sedges suppressing B. officinalis.
The transferred butterflies were released at two grassland patches with abundant B. officinalis, separated by 200 m. The central point between them (48.8243 N, 13.8014 E) is hereafter considered as the “release point”. The presence of the butterfly was checked annually by A.P., who is a local resident. Until the third year after the transfer, the butterflies were only observed at the release sites. Subsequently, they expanded to adjacent grasslands. In 2011 and 2012, a decade after the transfer,
As we aimed on both ascertaining the current distribution extent and elucidating habitat selection of the transferred population, we selected a minimalistic approach, which allowed addressing both questions during the short flight period of the butterfly. Hence, each grassland in our area of interests was visited only once, by zigzagging its entire area (cf.
The survey covered all grasslands and fenlands in wider environs of the release site: around the Nové Údolí deserted village, downstream (eastwards) and upstream (westwards) along Studená Vltava River, along its right-side tributary Hraniční potok brook, and still further north along the state border. We also surveyed grasslands and fenlands around Stožec village (4 km NE aerial distance from the release site, separated by 1.2 km of mature spruce forests) and Černý Kříž homestead (6.4 km E, through 4.7 km of mature forests) (Fig.
The surveys were carried out June 1‒10, 2023, 9–16 hours CEST, if weather was suitable for butterfly activity. As a “grassland patch”, we understand a patch of grassland or fenland discernible in the field by clear borders (woodland edge, a strip of shrubs or trees, a brook, road, stone wall), or by a distinct management (a fenced pasture, unmanaged fen). Using detailed aerial maps, we surveyed 212 such patches. The centroid of each patch was marked into mapy.cz mobile application, and the patch was zigzagged to detect the butterfly and assess its abundance (ordinal scale, 1: single individual, 2: ≤5, 3: ≤10, 4: ≤20, 5: >20 individuals). Each visit to a patch was further characterised by: Julian Date; closest Hour; Clouds (ordinal 1–3 scale, from clear sky to overcast); Wind (1: none, 2: some wind, 3: light breeze); Woody cover (trees or shrubs >1.5 m, % of total area); Slope (1: flat, 2: mild, 3: steep); Impenetrability (subjective assessment considering tussocks, height of sward and litter accumulation, ordinal 1–3); Wetness (1–3, from dry to waterlogged); Bistorta cover (% projection of total patch area); Nectar abundance (1: none, 2: some nectar, 3: abundant nectar); Nectar diversity (the number of flowering forbs’ species, ordinal scale, 1: ≤5 spp., 2: ≤10 spp., 3: >10 spp.); and Management (1/0, site mown or grazed vs. unmanaged; for this we utilised local knowledge of A.P.). Following fieldwork, we measured the area of each grassland patch (in square metres, square roots were used in computations); and the lengths of its edges distinguishing open, partly closed and closed edges (expressed as % of the entire perimeter length). Finally, we measured the distance of each patch centroid from the original release site (in metres).
Two response variables analysed in parallel were L. helle presence during the survey, coded binomially (present/absent), and abundance, coded at ordinal 1–5 scale.
Prior to analyses, we transformed all numeric and ordinal predictors to zero means and unit variance. We then visualised the relationships among thus transformed predictors using the principal component analysis (PCA) in Canoco, v. 5 (
We then used the general linear models (glm) in R 4.2.2 to model responses of L. helle presence and abundance to properties of the grassland patches. We used binomial link for presence and Poisson link for abundance and followed the information theory approach (= AIC values, which detect optimal points between competing models’ precision and complexity) for comparing the models’ performance.
We first constructed single-term regressions y ~ x, plus polynomial terms for Serial date and Hour, with all the potential predictors. Next, we defined two covariate models, the first accounting for mutual positions of the sites in space (Positional model: built by backward elimination from latitude, longitude and altitude, their interactions, and polynomials), the second just considering aerial distances from the release point (Distance model). Then, we assessed the predictors with “significant” effects (ΔAIC between null and fitted model ≤2.0) in single-term regressions against these two covariate models. Finally, we constructed multiple regressions, using stepwise elimination (drop1 function in R) from models, containing either Position or Distance, plus all predictors which were significant in the single- term regressions, until we obtained models with the lowest AICs.
We surveyed 212 grassland patches with a summed area of 285.7 ha (average patch 1.4 ± 1.04 SD, range 0.1–5.9 ha) (Figs
The PCA ordination (Fig.
PCA ordination biplot showing the mutual relationships among characteristics of the 212 individual grassland patches surveyed for presence and abundance of Lycaena helle in Šumava National Park. The narrow black arrows (and triangles, for the factor Management) are predictors used in subsequent regression analysis. The red arrows, dependent variables in the regressions, are visualised as “supplementary variables”, not influencing the mutual positions of predictors.
The regression models corroborated this. Both presence and abundance of the butterfly (Table
Results of single-term regressions, relating presence and ordinally coded abundance of adult Lycaena helle butterflies to predictors describing the 212 visited grassland patches around the Nové údolí site of release in Šumava Mts. Models with ΔAIC ≥2.0 compared with the null model are shown in bold.
Model | Presence | Abundance | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Coefficients | D.f. | Dev. | D2% | AIC | Coefficients | D.f. | Dev. | D2% | AIC | |
Null (y ~+1) | – | 211 | 270.3 | 272.3 | – | 211 | 382.4 | 574.1 | ||
~Serial date | -1.47x | 210 | 198.2 | 26.7 | 202.2 | -1.15x | 210 | 245.6 | 35.8 | 439.2 |
~ Serial date2 | -23.52x -5.20x2 | 209 | 194.6 | 28.0 | 200.6 | -18.72x -2.60x2 | 209 | 241.7 | 36.8 | 437.4 |
~Hour | -0.02x | 210 | 270.3 | 0.00 | 274.3 | 0.06x | 210 | 381.9 | 0.1 | 575.5 |
~Hour2 | -0.32x +3.17x2 | 209 | 268.0 | 0.9 | 274.0 | 0.71x +1.50x2 | 209 | 380.0 | 0.6 | 575.6 |
~Cloudiness | -0.65x | 210 | 254.1 | 6.0 | 258.1 | -0.56x | 210 | 346.3 | 9.4 | 540.0 |
~Wind | -0.34x | 210 | 265.1 | 1.9 | 269.1 | -0.35x | 210 | 363.2 | 5.0 | 556.8 |
~Area | -0.92x | 210 | 248.9 | 7.9 | 252.9 | -0.82x | 210 | 334.2 | 12.6 | 527.8 |
~Bistorta cover | 0.89x | 210 | 236.0 | 12.7 | 240.0 | 0.62x | 210 | 316.2 | 17.3 | 509.9 |
~Nectar abundance | -0.82x | 210 | 246.5 | 8.8 | 250.5 | -0.55x | 210 | 346.9 | 9.3 | 540.5 |
~Impenetrability | 0.17x | 210 | 268.9 | 0.5 | 272.9 | 0.08x | 210 | 381.2 | 0.3 | 574.9 |
~Slope | -0.30x | 210 | 266.2 | 1.5 | 270.2 | -0.14x | 210 | 379.4 | 0.8 | 573.0 |
~Woody cover | 0.08x | 210 | 270.0 | 0.1 | 274.0 | 0.08x | 210 | 381.3 | 0.3 | 574.9 |
~Wetness | 0.41x | 210 | 262.6 | 2.9 | 266.6 | 0.32x | 210 | 365.9 | 4.3 | 559.5 |
~Management | -1.42x | 210 | 258.6 | 4.3 | 262.6 | -1.05x | 210 | 363.4 | 5.0 | 557.0 |
~Nectar diversity | -0.54x | 210 | 259.1 | 4.1 | 263.1 | -0.35x | 210 | 366.8 | 4.1 | 560.4 |
~Closed edges | 1.36x | 210 | 264.3 | 2.2 | 268.3 | 0.23x | 210 | 373.7 | 2.3 | 567.3 |
~Partly closed edges | -0.27x | 210 | 267.1 | 1.2 | 271.1 | -0.15x | 210 | 379.2 | 0.8 | 572.8 |
~Open edges | -0.07x | 210 | 270.1 | 0.1 | 274.1 | -0.09x | 210 | 381.3 | 0.3 | 574.9 |
Position and Distance models explained even higher amounts of variation (Table
Components of the multiple-regression model Lycaena helle presence ~ Distance + Serial day +Serial day2 +Bistorta cover +Wetness. Serial day and Distance from the transfer site were input a priori, before stepwise elimination from all predictors that achieved ΔAIC >2.0 in single-term regressions presented in Table
Results of multiple regression models, relating presence and ordinal abundance of Lycaena helle to characteristics of 212 grassland patches in the Šumava Mts. distribution area of the introduced populations. The table shows both terms of position and distance models, used as covariates in further analysis; results of adding single predictors onto the covariate models, and resulting multiple-regression models.
Model | Presence | Abundance | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Terms coefficients | D.f. | Dev. | D2% | AIC | Terms coefficients | D.f. | Dev. | D2% | AIC | |
Null | – | 211 | 270.3 | 272.3 | – | 211 | 382.4 | 574.1 | ||
Positional | -28.59 Latitude -14.29 Latitude2 +2.14 Longitude -47.29 Longitude2 | 207 | 169.3 | 37.4 | 179.3 | -19.93 Latitude -7.06 Latitude2 -0.35 Longitude -34.83 Longitude2 | 207 | 203.0 | 46.9 | 402.6 |
Serial date | -0.82 Serial date | 206 | 165.1 | 38.9 | 177.1 | – | ||||
Bistorta cover | 1.39 Bistorta cover | 206 | 132.0 | 51.2 | 144.0 | 0.49 Bistorta cover | 206 | 163.0 | 57.4 | 364.7 |
Wetness | 0.53 Wetness | 206 | 162.6 | 39.8 | 174.6 | 0.29 Wetness | 206 | 191.5 | 49.9 | 393.1 |
Multiple regression | -25.40 Latitude -16.83 Latitude2 -4.48 Longitude -58.66 Longitude2 -0.93 Serial date +1.43 Bistorta cover | 205 | 127.7 | 53.1 | 141.7 | -15.92 Latitude -4.62 Latitude2 +0.09 Longitude -32.19 Longitude2 +0.49 Bistorta cover | 206 | 163.0 | 57.4 | 364.7 |
Distance | -2.05 Distance | 210 | 183.0 | 32.3 | 187.0 | -1.37 Distance | 210 | 228.2 | 40.3 | 421.8 |
Date2 | -9.58 Serial date -11.04 Serial date2 | 208 | 171.7 | 36.5 | 179.7 | -7.81 Serial date -6.43 Serial date2 | 208 | 208.9 | 45.4 | 406.6 |
Bistorta cover | 1.27 Bistorta cover | 209 | 146.6 | 45.8 | 152.6 | 0.49 Bistorta cover | 209 | 185.7 | 51.4 | 381.4 |
Nectar abundance | – | -0.23 Nectar abundance | 209 | 223.6 | 41.5 | 419.2 | ||||
Penetrability | – | 0.19 Penetrability | 209 | 222.9 | 41.6 | 418.6 | ||||
Wetness | 0.60 Wetness | 209 | 173.2 | 35.9 | 179.2 | 0.33 Wetness | 209 | 211.0 | 44.8 | 406.6 |
Slope | -0.33 Slope | 209 | 179.9 | 1.8 | 185.9 | – | ||||
Management | -1.21 Management | 209 | 178.8 | 33.8 | 184.8 | -0.64 Management | 209 | 222.0 | 42.0 | 417.7 |
Multiple regression | -2.98 Distance -10.28 Serial date -13.54 Serial date2 -0.41 Slope +1.32 Bistorta cover | 206 | 134.1 | 50.4 | 146.1 | -1.32 Distance -6.64 Serial date -4.77 Serial date2 +0.42 Bistorta cover +0.18 Wetness | 206 | 168.1 | 56.0 | 369.7 |
Twenty years after its establishment, the transferred Šumava Mts. population of L. helle inhabits 66.6 ha of waterlogged grassland patches, scattered across 4.5 km2 of mountain grasslands. Identically to other occupied areas in Western and Central Europe, the Šumava L. helle prefers waterlogged grasslands with high representation of its host plant, Bistorta officinalis, sheltered from winds by contiguous woody boundaries, typically with flat relief and no active management. It avoids drier grasslands on steeper slopes managed by grazing or mowing and containing higher diversity of nectaring forbs. Surveying grasslands in wider environs revealed that the butterfly has not yet expanded through wider stretches of contiguous forests.
The association with waterlogged sheltered grasslands with high host plant supply corroborates the findings from other regions such as Westerwald and Eifel Mts., Germany (
Interestingly, nectar diversity correlated positively with active management and negatively with high Bistorta cover. This suggests that management interventions increase the number of flowering plants, while suppressing the L. helle host plant (cf.
In 2011 and 2012, i.e. one decade after establishment the Šumava population,
Perhaps more interesting is the expansion rate following the transfer. Compared to the 900 m reported by
The 0.2 km/year expansion rate calls for comparing with the two earlier L. helle transfers, both carried out in France, to Morvan Mts. (Bourgogne-Franche-Comté, highest point Haut-Folin, 901 m elevation) and to Forez Mts. (northern projection of the Massif Central, highest point Pierre-sur-Haute, 1,631 m elevation) (
The situation also deserves comparison with Bog Fritillary Boloria eunomia (Esper, 1800) (Nymphalidae), another glacial relic (
The newly established L. helle population in the Šumava Mts. thus has been expanding at a slower rate than the previously transferred populations of the same species and both transferred and natural populations of an unrelated butterfly with similar habitat requirements. The Šumava L. helle has not yet trespassed barriers of contiguous spruce forests (cf. Fig.
Still, such species-specific circumstances do not explain the slower expansion of Šumava L. helle, when compared to the translocated populations in Morvan and Forez. A possible reason relates to geomorphology of the release localities. The Nové Údolí area discussed here is a shallow pan with high representation of grasslands, surrounded at all sites by higher and mostly forest-covered ridges. Comparing this with Morvan and Forez would require a detailed analysis of local orographies. Likely it took L. helle 20 seasons to colonise all the suitable grasslands within this area. Given that in 2023, L. helle was recorded at the highest point of the grasslands between Nové Údolí and Strážný, the further expansion to the northwardly situated grasslands will likely accelerate.
As the introduced population inhabits sites directly adjoining the border with Germany, expansion into Germany is fully expectable, and was indeed confirmed by records of L. helle from the village Haidmühle in 2021. The landscape and spatial configuration of habitats at the German side of the border differ from the Czech side, however, as there were no post-war transfers of inhabitants, a majority of the area is farmed, and potential habitats are more fragmented. A detailed survey of L. helle performance on the German side of the border is desirable and can provide an interesting comparison to its performance on abandoned land.
The newly established population of Lycaena helle in Šumava Mts., Czech Republic, is apparently viable, having expanded from the original release point to a wide area along the course of Studená Vltava River, and has potential for further expansion. Its preference for minimally managed successionally advanced grasslands do not differ from that of other mountain populations in Central and Western Europe. The transfer, initially motivated by fears for the donor population in the Türnitz Alps, Austria, established the butterfly in a previously unoccupied area, thus enhancing its total numbers in mountains of Central Europe, was a butterfly conservation success. The target area does not suffer the common problems encountered in insect conservation, such as intensive land use and resulting habitat homogenisation and landscape impenetrability. Ultimately, addressing these aspects in broad landscapes will be more decisive for butterfly biodiversity than single-species transfers. Encouragingly, the donor population was also still extant in 2023, and measures to sustain its habitats, including partial removal of the young spruce cultures from the former valley grasslands, have been adopted in the interim period.
The administration of National Park issued us all necessary permits for fieldwork in the strictly protected area. Matthias Dolek, Jan Christian Habel and Thomas Schmitt provided much valued comments in their reviews. The study was supported by Technology Agency of the Czech Republic (SS01010526).