Temporal Variation in Abundance and Distribution of Butterflies in a Southern Nigerian National Park

This study was carried out in Okomu National Park, Edo State, Nigeria, one of the parks created and managed by the Nigerian government. Okomu is about 200km² and located in Ovia South West Local Government Area of Edo State, Southern Nigeria. It lies within Latitude 06⁰ 20′ N and Longitude 0.05⁰ 16′ E (Fig. 1). The Park is drained by the Osse River, which is located at the eastern boundary and the Okomu River, from which the name of the park was derived and located at the western boundary. The land area covered by the Park makes it the smallest of the seven (7) National Parks in Nigeria (NNPS, Ojianwuna CC).

The climate of Okomu National Park and environs is relatively stable. Typically, the southern region has an important feature of the humid tropical wet and dry climate governed primarily by rainfall, which occurs between February and October, with the highest occurring in June, July, August, and September. The driest period is December, January, and March. Annually, temperature in Nigeria range between 11^oC and 38^oC (mean of 27.2^oC) and precipitation range between 2mm and 240mm (mean of 1162.5mm) (NNPS).

Based on the visual level of disturbance and the degree of canopy closure that was adopted from Barau et al.,the ranges were categorised into:

1. Scattered/Open forest (SOF) which indicated features of selective logging in disturbed forests and agricultural activities (plantation).

2. Closed Forest (CF) which showed moderately disturbed forests, trees in closed canopies, and thick vegetation.

3. Swamp Forest (SF), which showed moist patches along road and riparian habitats

4. Shrub (S) which showed grasses and herbaceous plants (Shrubs).

The straight line transect method was used where eight (8) transects were selected in the range of features based on the works of Barau et al. and [Ojianwuna and Akpan ]thus;

(1) T1(Disturbed forest at Compartment 54 Arakhuan (DIS) was at Latitude 06⁰ 20′ 43″N, Longitude 0.005⁰ 21′ 17″ E and 62m above sea level

(2) T2 (Forest Edge/Plantation Compartment 29 Julius Creek (FE/P) was at Latitude 06⁰ 24′ 23″ N, Longitude 0.05⁰ 19′ 23″ E and 82m above sea level for SOF,

(3) T3 (Moderately Disturbed habitat at Compartment 64 Arakhuan (MOF) was at Latitude 06⁰ 24′ 54″ N, Longitude 0.005⁰ 19′ 56″ E o AND 64 m above sea level

Fig. (1). A Map of Okomu National Park showing sampling points [Ojianwuna CC, et al.].

(4) T4 (Closed Forest at Julius Creek Compartment 30 (CC) was at Latitude 06⁰ 23′ 28″N, Longitude 0.005⁰ 17′ 06″ E and 39m above sea level for CF.

(5) T5 (Moist Patches on Road at compartment 54, Arakhuan (FAR/MP) was at Latitude 06⁰ 21′ 02″ N, Longitude 0.005⁰ 20′ 27″ E and 53m above sea level

(6) T6 (Riparian Habitat at Compartment 31, Julius Creek (RIP) was at Latitude 06⁰ 24′ 54″N, 0.005⁰ 18′ 00″ E and 10m above sea level for SF.

(7) T7 (Shrub at Compartment 89, Babui Range (SHR) was at Latitude 06⁰ 17′ 15″N, Longitude 0.005⁰ 21′ 36″E and 40m above sea level

(8) T8 (Grassland at compartment 28, Igwonwan Range (GRA) was at Latitude 06⁰ 23′ 49″N, Longitude 0.005⁰ 21′ 36″E and 82m above sea level at S.

Two observers walked the transect for 30mins, counting and obtaining a butterfly specimens.

Butterfly sampling was done between 9.00 am in the morning and 3.00pm in the afternoon when the activities of butterflies were high due to suitable sunny weather conditions. Butterflies were collected differently from the four ranges including Scattered/Open forest (SOF), Closed Forest (CF), Swamp Forest (SF), and Shrub (S). A total of 4 days was devoted for sampling in each month for 24 months (July 2012 to June 2014) under similar sampling techniques. Collection of butterflies was done twice in a month for two years using the Butterfly Trap Technique (BTT) [Venters, N.] and Butterfly Net Technique in the Transect (BNT) [Walepole MJ]. The butterfly trap technique was a passive method for butterfly collection while the butterfly net technique served as an active method. The cylindrical-like trap was made of green coloured net which was closed top and bottom with cylindrical ply board. The trap measured 500 cm length by 200 cm diameter of ply wood. The bottom ply wood of the trap had an opening size, enough to occupy the green bowl containing a mixture of banana and yeast which served as butterfly bait, and a slight opening to allow butterflies’ access to baits. These traps were strategically positioned under trees within the transect. Butterfly nets were used for active collection within transect. This was done at the beginning and end of the month for each sampled months to determine the richness and abundance of species or individuals. Butterfly net was used to catch butterflies in flight and was added to the associating transect. These methods were adopted from Fitzherbert et al. and Aduse-Poku et al. .All data from different days of sampling in each transect were pooled to give data for a month and were expressed as the number of species or individuals per trap /month. Butterflies were taken to the laboratory of Ambrose Alli University for identification using standard taxonomic keys by (Larsen).

Physical factors such as temperature, humidity, and rainfall were measured using the methods of [Brooks & Kyker-Snowman], with slight modifications. Randomness of sampling for physical variables within each transect was achieved by sampling from any point within (either sides or central point) the transect. RH/Temperature Data Logger (Lascor-EL-USB-2) was programmed to record the temperature and relative humidity on each forest floor transect. Three data loggers were installed in each of the transects before noon on the first day and later removed after three days. By removing and installing the units on consecutive days, all transects were sampled thrice. Rainfall was measured using the rain gauge Monsun radio-controlled rain gauge. This was placed 30m from the forest edge in an open forest about 50m to 100m wide cleared at the Rangers quarters in Arakhuan.

Environmental permit to use Okomu National Park Forest was obtained from the National Park Service, Abuja.

Physical attributes (Abundance and Distribution) of the transects were compared using the mixed model analysis of variance (ANOVA). Descriptive statistics were used to present temporal variations and variations in the taxonomic families. Generalized linear model analysis was done using the butterflies sampled from the different transects for the 24 months to ascertain R² and differences set at α= 0.05. Alpha diversity analysis was done and the magnitude of change between the transects seasons was compared. Shannon_H index, eveness, and Simpson_1-D index was used to determine species richness and abundance. All analysis was done using Graph Pad Prism and XL STAT (2021 Version).

The overall individuals of butterflies obtained from July 2012 to June 2014 in Okomu National Park, Edo State, are summarized in Table 1. A total of six thousand three hundred and ten individuals (6310) and 143 species were captured in Okomu National Park. According to decreasing abundance of butterfly family, Nymphalidae (3103: 49.2%), Pieridae (2269: 36.0%), Papilionidae (554: 8.8%), Lycaenidae (376: 6.0%) and Hesperiidae (8: 0.1%) dominated the sampled transects.

The summary of the generalized linear models of butterflies in various microhabitats in Okomu National Park, Edo State, Nigeria, is shown in Table 2. It showed that butterflies in the various microhabitats followed Poisson distribution except in RIP. Significant difference were observed in the butterflies in the various sampled transects (p<0.0001). Furthermore, the abundance and diversity index of butterflies by season in Okomu National Park, Edo State, Nigeria, is shown in Table 3. Of the 143 species sampled from the transects, higher species abundance was recorded in the wet season than the dry season. Similar trend was observed in the evenness, Shannon diversity, and Simpson dominance. There was a significant difference in the abundance of butterflies in the wet months (April-October) compared to the dry months (November- March) (P<0.05).

Relative abundance of butterflies between the sampled transects in Okomu National Park, Edo State, Nigeria, were generally significant (F _(Transects) = 24.66, P<0.05). However, within the different transects, significant abundances were equally observed. In the disturbed open forest (DIS) (Fig. 2a), Lycaenidae were highest collected mainly in the month August during the break of the raining period. Pieridae and Nymphalidae occurred throughout the months of sampling. The mean abundance of Pieridae and Nymphalidae were higher in July and August respectively when compared to the other months where a decline in the mean abundance was observed. Papilionidae was scarcely recorded while Hesperiidae was completely absent in this habitat. There was a significant difference in the butterflies obtained within sampled months (R²= 0.936; p<0.05).

In the forest edge/ plantation (FE/P) (Fig. 2b), individuals of Nymphalidae were collected throughout the year with peaks in August and Lycaenidae followed closely. Pieridae were equally sampled throughout the year and peaked in November. There was a decline in abundance during the dry season. Hesperiidae occurred only in July and Lycaenidae in the moths of July to September. The difference between butterflies obtained within sampled months was significant (R²= 0.928; p<0.05). In the moderately open forest (MOF) (Fig. 2c), all families of butterflies were collected throughout the years of sampling except Hesperidae which was not encountered. Most butterflies were abundant in the wet season and a decline was observed in the dry season. The abundance of Numphalidae in August was the highest and closely followed by their abundance in May. The difference between butterflies obtained within sampled months was significant (R²= 0.432; p<0.05).

In the Forest along Road with Moist Patches (FAR/MP) (Fig. 2d), Papilionidae, Pieridae, and Nymphalidae occurred throughout the years of sampling. Nymphalidae was highest in August and equally in May and June. Pieridae was highest in November. Abundances of Papilionid and Lycaenid butterflies were high in August. Hesperiidae was absent in this Transect. The differences in abundance between butterfly families were significant (R²= 0.825; p<0.05). The trend in abundance was the reverse in Riparian Forest around Rivers (RIP) (Fig. 2e). Within the butterfly families sampled in this transect, Pieridae was the highest and in December. Papilionidae and Nymphalidae were obtained in the dry season, November precisely. Hesperiidae were collected during the wet season in a very small abundance. Lycaenidae was highly occuring in August. The differences were significant (R²= 0.388; p=0.001).

The Closed Canopy (CC) transect (Fig. 2f) showed high abundance of Nymphalid butterflies in May and August. Although, a decline was observed afterwards. Pieridae in July was high and a rise in abundance was observed in November. Papilioniodae, Pieridae, and Nymphalidae were encountered throughout the sampling period. Papilionidae and Lycaenidae were mainly obtained in the wet season but in low numbers. Hesperiidae was also absent in this transect. The difference was significant (R²= 0.808; p< 0.05). In shrub habitat (SHR) (Fig. 2g), Nymphalidae was highest in September. Pieridae and Nymphalidae were collected throughout the sampling. Apart from Lycaenidae recorded in July, Lycaenidae and Hesperidae were absent in this transect. The difference was significant (R²= 0.707; p<0.05).

In grassland habitat (GRA) (Fig. 2h), the abundance of Nymphalidae was highest in August and was followed closely with their abundance in May. Pieridae was high in November. Pieridae and Nymphalidae were collected throughout the sampled period. Papilionidae occurred in lower numbers which cut across the wet and dry seasons. Hesperiidae and Lycaenidae were completely absent in this transect. The difference was significant (R²= 0.812; p<0.05).

Fig. (2a). Mean abundance of Butterfly per Families in DIS of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.936; F= 6.00; p<0.0001).

Fig. (2b). Mean abundance of Butterfly per Families in FE/P of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.928; F= 8.15; p<0.0001).

Fig. (2c). Mean abundance of Butterfly per Families in MOF of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.432; F=25.01 p<0.0001).

Fig. (2d). Mean abundance of Butterfly per Families in FAR/MP of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.825; F= 5.68 p<0.0001).

Fig. (2e). Mean abundance of Butterfly per Families in RIP of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.388; F= 11.35; p=0.001).

Fig. (2f). Mean abundance of Butterfly per Families in CC of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.808; F= 15.92; p<0.0001).

Fig. (2g). Mean abundance of Butterfly per Families in SHR of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.707; F= 3.79; p<0.0001).

Fig. (2h). Mean abundance of Butterfly per Families in GRA of Okomu National Park, Edo State, July 2012 to June 2014. (R2= 0.812; F= 4.58; p<0.0001).

Fig. (3). Relationship between the Physical parameters (Rainfall, Relative Humidity and Temperature) and butterflies of Okomu National Park from July 2012 to June 2014.

The relationship between monthly variations in rainfall, temperature, relative humidity, and butterfly abundance was recorded in Okomu National Park, Edo State, between July 2012 and June 2014 (Fig. 3). Rainfall ranged between 45mm to 785mm. Relative high rainfalls were recorded in the months of June, July, and August. The amount of rainfall reduced greatly in March and low rainfalls were observed in December, January, and February. Rainfall was positively correlated with butterfly abundance throughout the months of sampling (0.924).

The pattern of variation in relative humidity is similar to that of rainfall. It ranged from 65% in March to 90% in July. Being a rainforest and protected area, high percentage saturation was recorded in the months of July to September. March recorded the lowest saturation. Other months also recorded slightly higher relative humidity. Relative humidity was equally positively correlation with butterfly abundance through the sampling months (0.722).

Temperature ranged between 22.0^oC and 32.0^oC but slightly varied within the months, and there was no significant correlation between mean temperature and relative butterfly abundance throughout the sampling period.

Environmental factors influence the abundance and distribution of butterflies in tropical and subtropical regions. Rainfall is the most important physical factor directly affecting butterflies in the tropics and several studies have proven this [Kunte K, et al, 2009]. Rainfall is the key driving parameter for butterfly abundance and survival. According to [Ludecke et al. 2021], the main drivers of rainfall in Nigeria are the North Atlantic Oscillation (NAO, March), the Indian Oceanic Dipole (IOD, April and November), Atlantic Multidecadal Oscillation (AMO, June and September) and the Pacific Decadal Oscillation (PDO, September). Other factors such as temperature and moisture gradient have been reported to have indirect effect on butterflies, especially variation in abundance of adults [Barau K 2010, and Wolda H, et al, 1989].Often fluctuations in temperature between wet and dry months are minimal compared to rainfall which shows remarkable gradients. The study revealed that throughout the months of sampling, rainfall positively influenced butterfly abundance. It was observed that as the rains increased during the wet seasons, butterfly abundance in DIS, FE/P, MOF, FAR/MP, CC, SHR and GRA increased. The observation of [Ludecke et al. 2021] on rainfall pattern in months coincides with months with high butterfly abundance. However, a reverse trend of abundance was observed in RIP with low abundance of butterflies in the wet season. During the dry seasons, butterfly abundance was low in DIS, FE/P, MOF, FAR/MP, CC, SHR and GRA while RIP showed high abundance of butterflies. The seasonal trend could be attributed to synchrony with the phenology of food /host plants and less predation during the wet season which ensured the survival of larvae and adult butterflies [Wolda H et al., 1989, Spitzer K et al.]. Moreover, a significant reduction in rainfall that occurred during the dry season may have led to a decline in food availability and increased plant defoliation. Studies have shown that once there is a decrease in food supplies, there is an end to the annual cycle of most butterfly species, some which may migrate to more humid areas, aestivate or enter into diapauses, [Janzen DH et al.,1983, Martins CF et al. 1994]. This may be the reason there was low abundance of butterflies in the dry seasons, especially in the months of March when rainfall was lowest. The overall abundance of butterflies in RIP was highest in the dry seasons of February and March of the sampling period compared to the wet seasons. The presence of the river even in the dry season may be a contributing factor to the abundance of butterflies as they were observed in cluster around wet sands taking water, nutrients, and minerals for survival during the dry periods. They also tend to avoid the environment with extreme heat and seek cooler areas during this period. Studies have shown that riparian forests support more butterflies during dry seasons because of the vegetation around water, mud, wet sand with nutrients, and leaf litter compared to the wet months [Janzer DH et al.,1983, Vu LV et al., 2011].

High rainfalls in the study were obtained in the months of June, July, and August, but the highest abundance of butterflies was in August. The small windows of sunshine (August break) in between showers in August may be a factor responsible for the high abundance of butterflies. [Larsen et al.] explains that butterflies are encountered mainly in suitable sunny weather with adequate rainfall. The high butterfly abundance recorded in the months of June and July, which were the peak of the wet season, were not as high as the butterflies obtained in August. Studies have shown that some adults and larvae of butterflies may be killed or removed and eggs may be dislodged by heavy rainfall causing overflooding of the area. This may depress reproductive activities such as mating and egg laying and cause water logging of host plants. These may have led to the decrease in their numbers in those months of very high rainfall [Hammer KC et al.,2000], [Barau K, et al., 2010 - Arya MK, et al, 2014].

The abundance pattern of families of butterflies with respect to months indicated that Papilionidae and Pieridae species were seen all the year around and was significantly higher in the wet season (April to October) (Figures 2a-d and 2f-h). However, their numbers reduced during the dry season (November to March) when the areas became hot and the individuals searched for cooler areas. There was also a reverse trend of butterfly family abundance in RIP as most families were obtained in the dry season and low in the wet season. There were significant differences in the butterflies sampled in the various transects (p< 0.05). Again, the presence of moisture in RIP may be the reason butterfly families, especially Pieridae, Papilionidae, and Nymphalidae, were obtained during dry season in high numbers. Similar findings from [Nobre et al. 2012] and [Arya et al, 2014] showed high butterfly species in the wet season and significantly low species in the dry season. Papilio species were also observed in all transects. Researches have shown that physical variables like light, humidity can affect the abundance of Papilla species because their flight range gives them the opportunity to visit flowers, open habitats, and forest canopies. This may also be the reason they were seen in all habitats. Butterfly species were higher in the wet season than in the dry as well as the number of individuals, Shannon, evenness, and Simpson dominance index.

Lycaenidae had large numbers of individuals occurring at a particular season (July to September), probably when the host plants were abundant for the caterpillars to grow. [Kunte et al. 2000], [Kunte K et al. 2000] and [Rajagopal et al. 2011] state that the abundance of Lycaenidae at a particular season could be the time when the plants were in a suitable phenophase for the growth of the caterpillar. Hesperiidae were obtained in RIP and FE/P during wet season (April to May) in very low numbers. It may be as a result of their habits and time of sampling, usually found very early in the morning. [Kunte et al. 2000] also reported that they are seen active in early day time under shades of jungle or out in the open during cloudy/wet weather. This may be the reason they were found around RIP and FE/P habitats. Another reason for the low species of Lycaenidae and Hesperiidae may be as a result of their small size which may have restrained dispersal [Lodh JMR et al. 2016].