Population ecology and behaviour of two Afrotropical forest butterflies

Over the last decades, numerous natural habitats have been converted into settlement areas, agricultural land, and tree plantations on a large spatial scale. As a result, natural ecosystems have been destroyed. In consequence, many ecosystems exist today as small and geographically isolated remnants. To what extent the original species diversity can persist in such small habitat patches is questionable and strongly depends on the ecology of the species. A prominent example of severe habitat destruction are the species-rich tropical cloud forests of Taita Hills in southern Kenya, which have been deforested almost completely during past decades. However, there still exist typical forest species in the few remaining forest fragments. In this study, we investigate the population ecology and behaviour of two butterfly species present in the cloud forest remnants of Taita Hills, Protogoniomorpha parhassus and Precis tugela. Over a period of one month, we conducted Mark-Release-Recapture to study population sizes and demographic structures, lifespan, dispersal, and behaviour. We found that both species exhibited medium population sizes and are sedentary. However, some individuals performed dispersal throughout the forest. The behaviour of the two species differs: While P. tugela was mostly observed basking with open wings, P. parhassus was mostly sitting under leaves with closed wings. The life span was rather long for butterflies. This study documents the population ecology and behaviour of these two Afrotropical butterflies and underlines the relevance of the conservation of cloud forest remnants to preserve species, which mainly depends on these habitat remnants.


Introduction
Transformation of still intact ecosystems into agricultural land has accelerated greatly over the last decades and has caused the destruction of habitats, particularly in the tropics (Bradshaw et al. 2009). Thus, numerous ecosystems today only exist as small and geographically isolated habitat remnants. As a consequence of these changes in habitat configuration, edge effects negatively impact habitat quality (Thomas 2016). Although simply being small remnants, they frequently host a large proportion of the original biodiversity. However, it is questionable whether species are able to persist in such small habitat remnants in the long term .
Whether a species is able to persist under these conditions of habitat fragmentation and habitat deterioration strongly depends on its ecological demands and behaviour, in particular on its degree of ecological specialisation, adaptability, and dispersal (Thomas 2016). It has been shown that species with specific habitat requirements, comparatively narrow ecological variability, and little adaptability are suffering greatly under environmental changes, such as habitat fragmentation and reduction of habitat quality (Mangels et al. 2017). This means that dispersal ability and behaviour are pivotal for species' survival, especially when they formerly existed in interconnected ecosystems (Ehl et al. 2019). In this context, mobile species are assumed to cope much better with habitat fragmentation than sedentary and highly specialized taxa do (Thomas 2016). Additionally, the life span of individuals also has a strong effect on the respective species' persistence in a region or habitat because species with a longer life span of their individuals can bridge poor environmental conditions better than those with short-lived individuals, such as many insect species (Morris et al. 2008, Rodriguez-Caro et al. 2021. The cloud forests of Taita Hills in southern Kenya represent one of these cases of dramatic losses of natural tropical forests with extraordinarily high diversity (Pellikka et al. , 2013Teucher et al. 2020). The Taita Hills cloud forests are one of the world's biodiversity hotspots (Myers et al. 2000). Originally, cloud forests covered the entire higher elevations of this mountain ridge. However, almost all of these species-rich forests have been transformed into tree plantations composed of exotic species (such as Eucalyptus, Grevillea robusta, Cupressus lusitanica) or converted into arable land for subsistence agriculture ). Most of these transformations occurred within the last century ). As a result, only few remnants of the species-rich cloud forest remain today, mostly on mountain peaks and rugged slopes with limited accessibility. These small forest pockets still hold a remarkable flora and fauna, some of which are endemic to the Taita Hills (Brooks et al. 1998).
While numerous studies have intensively investigated the species composition of animal and plant species, also considering habitat quality and invasive species, dedicated studies on population ecology and behavioural biology of specific species are still underrepresented, especially for tropical ecosystems. In this study, we selected butterflies as study organisms to analyse their population structure and behaviour. Butterfly species are excellent model organisms as most of them respond to and are highly sensitive to environmental changes (Bonebrake et al. 2010, Rákosy andSchmitt 2011). This is because they often require specific larval food plants and (micro)habitat structures as imagoes, as well as particular microhabitats for the developing larvae. Consequently, butterflies are also often used as indicators to evaluate the status of an ecosystem and/or of entire landscapes (Fleishman and Murphy 2009). Here we studied the population ecology and behaviour of butterflies. For this we used Mark-Release-Recapture (MRR), a powerful tool to analyse population sizes, demographic structures (including sex and age), and the behaviour of the individuals including dispersal (e.g., Ehl et al. 2018). Such information is highly suitable and reliable for the evaluation of species' persistence in small and isolated habitats and for the development of efficient nature conservation measures.
For our MRR study we selected two nymphalid butterflies, the forest mother-of-pearl Protogoniomorpha parhassus (Drury 1782) and the African leaf butterfly Precis tugela (Trimen 1879). Both species occur in high densities in Taita Hills and can be observed mainly in and around the cloud forest patches and adjoining gardens (Kioko et al. 2021). Both species are rather large butterflies (with wing-spans of 75-100 mm for P. parhassus and 55-64 mm for P. tugela) and do not show sexual dimorphisms. These butterfly species are common, distributed across major parts of Subsahara Africa (Larsen 1991(Larsen , 2005, and show comparatively strong dispersal behaviour (Kioko et al. 2021). They are mainly found in forest and woody habitats (Larsen 1991(Larsen , 2005Kioko et al. 2021). Their larvae feed on a variety of plants of the understorey (Larsen 1991(Larsen , 2005. In addition, P. parhassus belongs to the species whose imagoes at least occasionally feed on rotting fruits (Molleman et al. 2006); their perching behaviour is reducing with increasing habitat degradation and perching sites are more quickly reoccupied in degraded than in non-degraded forest fragments (Bonte and van Dyck 2009).
Here, we address the following research questions: 1. What is the population size and demographic structure of the local occurrence of the two species? 2. What is the dispersal and behaviour of these two species, and do they differ in these aspects? 3. What conclusions can be drawn from these findings for practical conservation, especially in context of ongoing forest destruction and the deterioration of habitat quality?

Study site
Study sites were selected in the Chawia forest fragment of Taita Hills (3° 29′ 06′′ S, 38° 20′ 45′′ E), where both species still occur in high densities Kioko et al. 2021). This forest fragment covers an area of 86 ha of cloud forest and is located at the southern edge of Taita Hills, at an elevation of 1400-1600 m a.s.l. (highest point: 1607 m) ( Fig. 1). Chawia forest is not gazetted under the forest conservation and management act of 2016, and thus is a community forest managed by the community forest associations (Himberg et al. 2009;Teucher et al. 2020). In consequence, this forest remnant suffers from diverse human activities and disturbances such as a dense network of roads and paths, collection of dead wood, planting of exotic tree species, infrastructure programmes (damming of small streams, affecting the hydrology of the forest), and hunting of animals (e.g. monkeys) (Rülke et al. 2020;Nzau et al. 2022). Thus, the Chawia forest fragment represents a highly disturbed cloud forest patch (see Pellikka et al. 2009). Nevertheless, the forest is still largely characterised by original cloud forest vegetation. It is surrounded by fields for subsistence agriculture (mostly annual maize and interspersed perennial banana cultivation) and tree plantations of mainly exotic and fast-growing tree species such as Eucalyptus spp., Grevillea robusta, Cupressus lusitanica (Pellikka et al. 2013).

Mark-release-recapture
MRR was conducted along already existing paths running along the forest edge and throughout the forest interior (in total about 3.5 km) (Fig. 1). The observation trail was walked at slow pace once a day (between 9:00 h and 16:00 h) from 2 to 30 September 2021, under largely sunny weather conditions and temperatures of at least 20 °C. Later, from 1 to 30 October non-systematic MRR was continued in the area, but data were only included in the calculation of dispersal distances and the longest life-span of individuals observed. Encountered individuals of P. parhassus and P. tugela were observed and captured with a butterfly net. We determined sex (male/female) and age according to wing condition considering the following three categories: 1 (completely undamaged wings including all marginal cilia), 2 (undamaged wings but without cilia), 3 (wings with moderate to heavily damaged wings) (see Zimmermann et al. 2005;Wendt et al. 2021). Individuals were then marked with a consecutive number at the underside of the hind wing using a waterproof pen and then released. For each marked butterfly, we noted date and time of capture, as well as respective weather conditions. Its exact location was determined using a GPS device (GARMIN GPSMap 64). Recaptures were treated similarly but labelled with "R" (= recapture) in our field book. In addition, the observed behaviour of each individual prior to capture was noted. Here we used the following categories: Flying, basking with open wings, resting with closed wings, or nectaring.

Statistics
The demographic structure was analysed under an open population formulation of the Jolly-Seber model, allowing entries and exits of the individuals to and from the population (POPAN model; Schwarz and Arnason 1996). The model estimates the following parameters: probability of an individual to enter the population (p ent ), survival probability (φ), capture probability (p), and gross population size (N). The specimen's sex was treated as individual covariate. The parameters φ, p ent , and N were tested for their dependence on sex by fitting multiple models with all combinations of sex-dependent and -independent parameters. φ, p, and p ent were additionally tested for linear time dependence. Quadratic dependence on time and interactions of sex and time were tested for φ, and p ent (a full list of models is given in Online Appendix 1). In total, we ran 144 models for each of the two species (6 φ × 6 p ent × 2 p × 2 N). We chose the model that fitted our data best by comparing corrected AIC scores. All analyses were conducted in RMark v.2.2.7 (Laake 2013), which served as a frontend to software Mark (Linux v. 9.0, White and Burnham 1999). We tested for differences in mean dispersal distances between females and males per species with two-sided t-tests in R (v. 4.2.0; R Core Team 2022). Distance measures between one and a subsequent capture locality of the same specimen were therefore log-transformed. Furthermore, the travelled distances were divided into distance classes (20 m, 30 m, and 50 m intervals, respectively), separately for each species and sex. To check for any potential artefacts caused by the chosen interval sizes, we analysed and compared these three interval size classes. The inverse cumulative percentage of these classes was determined, which corresponds to the probability density function, i.e. the dispersal kernel.
Based on these classes, the probabilities of dispersal flights were assessed by distance extrapolation. We used two frequently applied regression analyses, the negative exponential function (NEF) and the inverse power function (IPF). NEF tends to underestimate rare long-distance movements, whereas the IPF may encounter problems with 'zero' movements. The data were linearly transformed with a semi-ln plot for the NEF analyses or with a double-ln plot for the IPF analyses. In both equations, 'P' stands for the proportional probability that an individual will travel at least as far as the distance D, and 'a' for the intercept of the regression. NEF operates with the dispersal constant K as slope, whereas IPF uses the variable n as slope, which represents the effect of distance on dispersal (Ehl et al. 2018).
We selected the best model and the most suitable interval size, based on calculated stability indices (i.e. R 2 ) of the calculated curves, which corresponds to the proportion of explained variance of the dependent variable by the independent variable. This allowed extrapolations of the population's proportion that should travel distances exceeding the extent of the study area. The calculations were performed separately for males and females.
In both species, the models that fitted the data best suggested the entrance probability p ent to decrease over time (linear estimate in P. parhassus, quadratic in P. tugela; see Online Appendix 1). While no sex-dependence of any parameter was estimated for P. tugela, the parameters population size N and survival rate φ differed among sexes in P. parhassus. In contrast to P. tugela (607 females ± 90 SE; 610 males ± 90 SE), a strong excess in males was estimated for P. parhassus (205 females ± 23 SE; 793 males ± 108 SE). Mean minimum ages (time difference between first and last capture) were slightly higher in males of both species. Maximum ages were 23 days and 26 days for P. tugela and P NEF = ae −kD or ln P = ln a − k D P IPF = aD −n or ln P = ln a − n(lnD) P. parhassus, respectively. Daily survival rates of P. parhassus varied over time (estimate quadratic; median females 1.00 ± 0.02; males 1.00 ± 0.06). Median daily survival rates of P. tugela were 0.91 ± 0.01, irrespective of sex (see Fig. 2 ). The age structure assessed by wing conditions did not change remarkably over time (Fig. 3).
Individuals largely remained in the same area throughout the study period. Individuals of P. tugela were generally recaptured farther away from their place of capture (female median = 33.9 m; male median 85.7 m) than individuals of P. parhassus (female median = 20.1 m; male median 27.1 m) (Fig. 4). Males tended to higher mobility in both species. However, no significant differences were found in any comparison. Few individuals dispersed exceptionally far (identified as outliers in box-whisker plots, i.e., deviating more than 1.5 times the interquartile range from the upper quartile). While both sexes showed such disperser individuals, they were more frequent and more extreme in males of both species (maximum distances P. tugela: 1,368 m; P. parhassus: 1,372 m) than in females (maximum distances P. tugela: 721 m; P. parhassus: 675 m).
We applied the NEF and IPF functions for extrapolating the potential for long-distance dispersal. IPF returned higher stability indices (R 2 ) in all pairs of functions. In all cases, explained variance was higher in males (range R 2 values based on IPF: 0.94-0.97) than in females (0.79-0.90); consequently, estimates for males should be more reliable than for females. R 2 values within both functions were quite similar for the three different distance classes within sexes and within species (Table 2). So, as no relevant differences existed between these classes, only the results for the 30 m intervals (as the intermediate one) are given below; all results are given in Online Appendix 2 and 3.
Extrapolations for 1 km were roughly similar for both functions within both species but estimates declined much stronger for NEF than for IPF with longer distances. Thus, following NEF, longer dispersal events would be little likely for both species, while they are expected under the IPF model. In general, following the better supported IPF model, the likelihood of reaching a defined distance was twice as high for P. tugela than for P. parhassus. Thus, about 10% of P. tugela males and females should reach distances of 1 km or more, and still 3 to 4% should travel 5 km or more. No difference in dispersal behaviour was obtained between sexes. For P. parhassuss, about 5.1% of males and 6.4% of females are assumed to disperse at least 1 km. With distances of 5 km or more, these values decreased to 1.3% (males) and 2.2% (females). Thus, the slightly female-biased dispersal in P. parhassus seems to increase with distance ( Table 3, Fig. 4). A Spearman rank correlation did not reveal a correlation between dispersal distance and the elapsed time between capture and first recaptures for females or males.
P. tugela and P. parhassus did not show markedly different behaviour in terms of flight, interaction, or feeding. A conspicuous difference, however, was observed for the resting behaviour. While P. parhassus was mainly observed resting with closed wings, P. tugela rested mainly with open wings (Fig. 5).

Population size and demography
Population sizes and respective densities estimated for both study species are largely consistent with results from previous MRR studies of tropical forest butterflies (e.g. Marín et al. 2009;Vlasanek et al. 2013;Vlasanek and Novotny 2015;Seixas et al. 2017). Only few of these so far studied species exhibited remarkably higher densities than in our study (e.g. Marín et al. 2009;Vlasanek et al. 2013;Vlasanek and Novotny 2015). Butterfly species, which occur predominantly in forests and woody landscapes, show similar values, such as Euphydryas maturna in Czech Republic (Konvicka et al. 2005). In tropical forests, however, butterflies can also accumulate at very high densities at certain locations, such as around watering holes, spots with high numbers of blossoms, locations with overripe fruits , as well as in forest gaps (Hill et al. 2001). This would also apply to P. parhassus, which can be found in high densities at watering places and feasting on overripe fruits (Molleman et al. 2006).
The medium sized population densities found for these two species contrasts with the sometimes very high densities of species, which occur restricted to very specific habitat types such as bogs and wetlands, or in various types of grassland ecosystems Zimmermann et al. 2011;Weyer and Schmitt 2013;Ehl et al. 2017Ehl et al. , 2018. Such high population densities are rarely reached by butterflies predominantly found in forests and wooded ecosystems. These differences in butterfly densities might be due to the spatial distribution of resources, which are rather homogenously distributed across forested areas (Mbuthia 2003). Consequently, the butterflies might be less abundant and are more widely distributed across the entire ecosystem and landscape. However, the comparatively low population densities found here might also reflect a methodological artefact: Since in forest observers are largely bound to walk along existing pathways (and thus only cover a limited proportion of the ecosystem, also excluding the highly relevant canopy ecosystem), observers may cover the entire area of a butterfly's habitat in open ecosystems and thus may detect a much larger fraction of the individuals which are present. Another factor shaping population densities are shifting environmental conditions. In temperate regions, life spans of butterfly individuals as imagoes are mostly restricted to 1 3 some days to 2-3 weeks (e.g. Zimmermann et al. 2011;Pennekamp et al. 2014;Ehl et al. 2017). Hence, they build up spatially and temporally constricted populations (e.g. Fric et al. 2010;Weyer and Schmitt 2013). Similarly, arthropods of the arid and semi-arid tropical ecosystems also follow annual population cycles, triggered by dry and rainy seasons. These seasonal shifts frequently result in strong fluctuations in population densities (see Schmitt et al. 2021). However, most butterfly species of wet tropical forest ecosystems occur year-round and form much less or even no distinct generations, in contrast to species of temperate regions (Molleman et al. 2007;Seixas et al. 2017). Furthermore, they are often much more long-lived (see Beck & Fiedler 2009). Studies on Afrotropical butterflies in western Uganda showed comparatively long life spans (mean 39 days) and a longevity   (Molleman et al. 2007). This trend corresponds with our own results for both, P. parhassus and P. tugela, showing rather long life-times (including a longevity record of 52 days for P. parhassus). Interestingly, our data analyses yield a higher population size for P. parhassus males than for females. We assume this difference in our analyses being an artefact resulting from the much higher number of marked males in combination with lower recapture rates in females (14.3%) than in males (17.6%). This apparently has led to lower female population size estimates and a generally worse coverage of females. This contradicts the results of breeding experiments, which generally showed a balanced sex-ratio for butterflies (Ashoff and Schmitt 2014;Seixas et al. 2017). The observed imbalance in sex ratios could be due to several reasons: 1. For other butterfly species, females show a rather cryptic behaviour while males actively patrol and thus are more visible (Seixas et al. 2017, with references therein); distinct patrolling behaviour has already been demonstrated for P. parhassus, especially in disturbed forest fragments of Taita Hills (see Bonte and Dyck 2009). 2. Females may be more abundant in the high canopies and thus less likely to encounter observers on the ground. 3. Population dynamics could be subject to protandry (at the beginning of the generation males are much more abundant than females, gaining dominance only from the second half of the generation period, e.g. Weyer and Schmitt 2013;Pennekamp et al. 2014;Ehl et al. 2019); however, such protandry is the exception for tropical butterflies, and our results also show neither a clear temporal population dynamic nor significant changes in sex ratios over time. Further studies are essential to clarify these contingencies.

Mobility and behaviour
Our data provide important insights into the mobility behaviour of the two studied species. Most marked and recaptured individuals were sedentary and remained close to their location of first capture. Most individuals did not move farther than 100 m, i.e. P. parhassus males: 72%, females: 87%; P. tugela males: 69%, females: 65%. This finding is consistent with the results from numerous other MRR studies on butterflies, from the tropics and the temperate regions (Konvicka et al. 2008, Fric et al. 2010, Vlasanek et al. 2013, Vlasnek and Nowotny 2015, Seixas et al. 2017, Ehl et al. 2018, Wendt et al. 2021. Most individuals of the studied non-migratory butterfly species remained close to their point of capture with fairly restricted movements except for some few individuals performing long distance dispersal (Lewis 2001, Marín et al. 2009, Vlasanek et al. 2013, Vlasnek and Nowotny 2015, Seixas et al. 2017. A few individuals in our study also dispersed throughout the forest, with maximum distances of 1,372 m (P. parhassus) and 1,368 m (P. tugela). However, extrapolation of data calls for considerably higher dispersal capacity of the more generalist P. tugela than in P. parhassus, being more strongly linked to forests. Compared to many other butterfly species (e.g. Konvicka et al. 2005Konvicka et al. , 2008Zimmermann et al. 2005Zimmermann et al. , 2011Fric et al. 2010;Weyer and Schmitt 2013;Ehl et al. 2017Ehl et al. , 2018Ehl et al. , 2019Wendt et al. 2021), and in contrast to the sedentary behaviour of the majority of the individuals assessed in our study, the proportion of P. parhassus individuals dispersing over distances of kilometres is high and even very high in P. tugela. Long-distance dispersal even across non-habitat areas therefore seems likely for the less specialized species P. tugela. This species therefore might interconnect its populations all over the range of Taita Hills by dispersal. In case of a local extinction, recolonization should be realised quickly for this species. In general, such long-distance movements (even if they rarely take place) are highly relevant, as they enable habitat (re)colonization and ensure the exchange of individuals among single forest fragments. These dispersal events create the pre-condition for long-term persistence of local populations. The lower dispersal power and stronger habitat link of P. parhassus makes this scenario less likely in this second case, making exchange between forest fragments and recolonization considerably less likely for this species.
Because we captured and marked butterflies in this study exclusively along a defined path within and at the immediate edge of the forest, it is possible that movement distances were clearly underestimated. This is a general problem that occurs in MRR studies, since trapping is usually only done where there is a high probability of capturing many individuals of a particular species. Thus, the situation arises that only within the habitat, and rarely outside in the landscape matrix, data are collected-and thus movement distances are generally significantly underestimated. We also found that the two studied species differed in their behaviour. While P. parhassus can be observed sitting under/on vegetation with closed wings inside the shady forest, (especially under large leaves of various shrub and tree species), P. tugela spends a lot of time basking in the sun with open wings. This butterfly species was much more prevalent along the sunny forest edges as well as in adjacent gardens and shrub structures in the agricultural landscape (Schmitt et al. 2020, own observations JCH, TS). Thus, the two species can be clearly distinguished from each other with respect of their ecological demands and behaviour. P. parhassus is a typical species of humid forests and thus occurs mostly in shady areas in the forest of our study area where it mostly shows a sedentary behaviour. On the contrary, P. tugela is more mobile and flies longer distances, often encountered beyond the forest ), making it a less typical forest species with a variety of alternative habitats (Larsen 1991;Kioko et al. 2021).
To get a more realistic picture of the population ecology, dispersal and behaviour of these two butterfly species, additional studies with further methods would be useful. For example, larger areas could be equipped with bait traps in order to record larger movement distances of the butterflies. In addition, this method could be used to attract more butterflies and thus significantly increase the sample size. Thus, data on density and movement distance would become more reliable. Here, attention should also be paid to longer data recording in order to be able to make better statements about the lifespan of the butterflies. However, the results assessed with bait traps could also provide various shortcomings which has to be taken into consideration while interpreting the outcomes. For example, butterflies would be actively attracted with bait traps and thus might bridge greater distances than they do on average .

Conservation consequences
Our study shows that most individuals of both butterfly species (with a few exceptions) are site-faithful. P. parhassus occurs in restricted areas and seems to depend on cloud forest (see also Schmitt et al. 2020). The estimated dispersal behaviour of this species suggests that an exchange of individuals between the existing forest fragments is unlikely due to the geographic distance between fragments, and even dispersal via plantations with cloud forest elements might be difficult to cross. This might be different for P. tugela because of its higher dispersal capacity and its lower habitat specialisation . Nevertheless, strict protection and reconnection of the cloud forest remnants and increasing the sizes of forest fragments would be an important conservation measure (see Aben et al. 2016), which may help to increase population and species persistence, especially of forest species. The removal of Eucalyptus trees (and  other alien plant species) would also be an important measure to restore this disturbed and degraded forest. This might increase habitat quality and support survival probability of all species living there.