Enhancing semi-aquatic species representativeness beyond protected areas: dragonflies in networks of conservation corridors

Complementarity is crucial when prioritizing sites for biodiversity conservation. Networks of conservation corridors (CCs) can contribute to regional representativeness by complementing biodiversity features included in existing protected areas (PAs). We ask whether criteria important for CC management and design are effective at prioritizing complementary sites, and how the consideration of species represented in PAs influence criteria performance. We focused on species turnover of generalist and specialist dragonflies across 88 riverine sites. Criteria assessed included site-level estimates of dragonfly species richness, estimates of local habitat quality and corridor width. Measures of local habitat quality were based on either dragonfly indicator species or proportion of alien vegetation. Results showed that CCs complement dragonfly diversity in PAs by contributing unrepresented generalist species. Of the criteria, corridor width was the most efficient at prioritizing complementary sites, while prioritization based on dragonfly indicator species or species richness underperformed. When aiming to prioritize CC sites that also complement sites situated in established PAs, wide corridors with low levels of alien vegetation should be favoured.


Introduction
Off-reserve areas are important for ongoing biodiversity conservation (Dahal et al. 2014;van Schalkwyk et al. 2019;Kane et al. 2022). Areas beyond protected area (PA) boundaries, such as conservation corridors (CCs), can add to regional connectivity and reduce the negative impacts of habitat loss and fragmentation by supporting dispersal among otherwise isolated habitats (Gilbert-Norton et al. 2010). Networks of CCs (also referred to as ecological networks) (Jongman and Pungetti 2004;Samways 2007) are especially beneficial when regionally employed alongside other ecological components such as stepping-stone habitats, core areas, and buffer zones. Together, these components enable ecosystems to exist in human-dominated landscapes (Forman 1995;Hess and Fischer 2001).
Corridors can further help maintain regional biodiversity by complementing biodiversity features (e.g., taxonomic diversity) already represented in established reserves. Important for identifying complementary sites is the concept of species turnover: the substitution (i.e., replacement) of species by other species across sites (Baselga and Orme 2012). For example, significant differences in species assemblages between networks of CCs and PAs suggest potential complementary value of CC sites . Complementarity can aid in identifying areas for conservation investment, i.e. conservation prioritization (Wilson et al. 2009). Importantly, complementarity can improve representativeness in an efficient manner, e.g., enhance representativeness with minimum management cost (Kukkala and Moilanen 2013). Representativeness is related to the need for reserves to represent the full variety of biodiversity, and involves the extent to which biodiversity features occur within a set of sites (Margules and Sarkar 2007;Heink 2009). Improving representativeness can be invaluable in regions where PA establishment is biased (e.g., concentrated in economically unviable areas), and where existing PAs add to the conservation of certain features that are already well-represented at the expense of others.
Site-level estimates of diversity are routinely used to guide CC management and planning van Schalkwyk et al. 2022). For example, wide corridors representing good quality habitat are associated with high species richness (Pryke and Samways 2003;van Schalkwyk et al. 2022). However, complementarity requires information on biodiversity content (rather than summary statistics such as species richness) as well as the joint consideration of multiple sites (Wilson et al. 2009). Dragonflies and damselflies (Odonata) (from here on referred to as dragonflies) have been extensively studied as indicators of ecosystem change, both terrestrial and riparian (Clausnitzer 2003;French and McCauley 2019;Vorster et al. 2020), and have also been incorporated in CC assessments (Kietzka et al. 2015;Mossman et al. 2015;Pryke et al. 2015). To date, however, not many assessments incorporate dragonfly species turnover (but see Pryke et al. (2015), and it is unknown whether promoting site-scale species diversity in CCs will encourage complementarity. In addition, conservation priority of potential sites will change as other priority areas are established and expanded, and so the process of conservation prioritization should be informed by what is already conserved (Wilson et al. 2009). The evaluation of CC site complementarity should therefore consider the influence of species represented in existing PAs.
We evaluate the efficiency of different criteria for prioritizing complementary corridors aimed at enhancing taxonomic representativeness of South African dragonflies. We focused on adult male dragonflies in streams of varying quality (in both corridors dissecting forestry plantations and PAs) (Kietzka et al. 2021a) and assess criteria that are often used at the site-scale to guide assessments, management, and planning of CC sites: species richness, local habitat quality, and corridor width (Hansen et al. 2018;van Schalkwyk et al. 2020;Kietzka et al. 2021a). Corridors were ranked according to these criteria, and we compared how different criteria influenced the gain in species representation with increasing number of selected sits. The efficiency of criteria for prioritizing complementary CC sites was assessed across two planning contexts: species represented across CC sites, and species represented across both CC and PA sites together.
Dragonfly species differ in their sensitivity to habitat degradation and disturbance. Some species are generalists that are tolerant and even thrive in transformed habitats, others are highly sensitive specialists that occur mostly in wild areas (Samways and Simaika 2016). We expect CC complementarity and criteria efficiency to differ between specialist and generalist species. By contributing to the decline of specialist species, landscape-scale land use change may contribute to assemblage homogenization and loss of turnover (Ball-Damerow et al. 2014). Corridors are situated in more transformed landscapes than PAs, and we therefore predict that CCs will complement the dragonfly species present in existing PAs through better representation of generalist species. Overall and irrespective of species sensitivity, we expect corridor width to be the most efficient criterion at selecting complementary sites, as wider corridors likely experience lower extinction rates, higher immigration rates and greater habitat heterogeneity than narrower corridors (MacArthur and Wilson 1967;Lees and Peres 2008). Furthermore, we predict that criteria representing changes to local habitat quality will be more efficient at prioritizing complementary CC sites for specialist species than generalist species.

Study area and sites
The study region is situated in the KwaZulu-Natal Midlands, South Africa ( Fig. 1). Networks of CCs in this region consist of natural and semi-natural vegetation that intersect plantation forestry blocks of alien trees. The data used were collected for a previous study (Kietzka et al. 2021a). A methodological description is provided below, and additional details can be found in Kietzka et al. (2021a).
In total, 110 study sites were selected along rivers situated in CCs (55 sites) and PAs (55 sites). These sites varied in factors previously shown to be important for South African dragonfly diversity, including corridor width (in the case of CC sites) and invasive plant species cover (Samways and Sharratt 2010;Kietzka et al. 2015Kietzka et al. , 2017Deacon et al. 2019). Not all 110 sites were considered in the present study. We removed six PA sites that were situated far from other sites so that inter-site geographical distances were more similar among CC and PA sites, so as to not artificially inflate species turnover between CCs and PAs. We removed six sites where no dragonflies were observed (see "Dragonfly sampling and species sensitivity" section below for details on sampling methods) and 10 sites that were situated close to each other (< 500 m). This left 43 PA sites, and 45 CC sites (88 sites in total). Minimum distance between any two sites was ~ 574 m, the mean was ~ 46 km, and the maximum distance was 106 km.

Dragonfly sampling and species sensitivity
A site consisted of a 100 m transect parallel to a river and included the 5 m area to each side of the transect. Each site was visited twice between 09h00 and 15h00 on sunny, windless days during the summer months (between January 2018 and March 2018). During each visit, male dragonflies were surveyed by two observers for one hour using close-focus binoculars (total of four observation hours per site). Where specimens were difficult to identify in flight, they were captured using a hand-net and kept as voucher specimens, which are held at Stellenbosch University, South Africa.
Species were split into two groups that differ in sensitivity to changes in the natural environment. This classification was based on the Habitat Sensitivity Sub-Index (HSS) of the Dragonfly Biotic Index, which uses the presence of dragonfly species (mostly adults but also includes complementary data on other life stages where it is necessary to establish total local species composition) to assess freshwater conditions (Samways and Simaika 2016). The HSS relates to how sensitive a species is to changes in its natural environment and is based on historic occurrence records in natural versus disturbed sites. Tolerant generalist species include species that are more common in artificial waterbodies than natural environments, as well as species equally common across both natural and disturbed habitats (species that have a HSS score of either 0 or 1). Sensitive specialist species include species that are mostly observed in fully natural environments, as well as species highly sensitive to environmental change (species that have a HSS score of either 2 or 3).

Criteria for site prioritization
We compared four criteria for prioritizing CC sites: (1) dragonfly species richness, (2) local habitat quality based on dragonfly indicator taxa, (3) local habitat quality based on the percentage alien vegetation at a site, and (4) corridor width ( Table 1).
Species richness was based on the observed species richness for each site (summed over both site surveys). A local estimate of habitat quality was based on the HSS (Samways and Simaika 2016). For each site, the HSS score for all observed species were summed and then divided by the number of species observed at that site. More details regarding the sub-index scores of the Dragonfly Biotic Index are provided by Samways and Simaika (2016).
Invasive alien vegetation remains an important threat to insect diversity (McGeoch 2002;Cardoso et al. 2020), woody species being especially problematic for dragonflies in South Africa (Samways and Sharratt 2010). A visual estimate of alien vegetation cover for each site was based on the most common species in the area, which included (in order of prevalence across sites): American bramble (Rubus cuneifolius), black wattle (Acacia mearnsii), bugweed (Solanum mauritianum), and scotch thistle (Cirsium vulgare). Each 100 m transect was subdivided into 10 plots. At each plot, two observers independently estimated alien vegetation cover, according to the aforementioned categories. Estimates recorded by the two observers were averaged and used to calculate the percentage of alien vegetation cover for each 100 m site.
Corridor width is an important design variable in conservation planning (Noss 1987;Beier 2018). Wider corridors are associated with higher terrestrial and aquatic arthropods species diversity (Pryke and Samways 2012;Kietzka et al. 2021a). We measured corridor width from aerial photographs as the linear minimum distance between two opposing corridor edges that intersects the centre of a transect.
None of the criteria showed spatial autocorrelation ( Figure S1), nor were any correlated with each other ( Figure S2).

Analyses
All analyses were conducted in R version 3.6.0 (R Core Team 2020).

Sampling efficiency
To determine dragonfly sampling efficiency, we used sampling completeness, which was calculated using the final slope of the species accumulation curve. This slope indicates the expected gain in diversity when sampling a new individual, and was calculated with the following formula: Here, n s denotes the total number of individuals, and n s−1 represents the number of individuals at the point in the curve where the last single species was added (Cardoso et al. 2011). This final slope was calculated using the R package BAT (Cardoso et al. 2021).

Criteria efficiency for representation
The efficiency of criteria to prioritize complementary CC sites was assessed using a patternbased approach based on the Environmental Diversity strategy of Faith and Walker (1996a). Sites were projected in a two-dimensional space (or 'ordination space') where the inter-site distances represent species turnover (i.e., sites situated further apart have greater species turnover) (Fig. 2). These distances were used to infer the relative complementary value of a site (e.g., sites that have unique species compositions would separate out further from the other sites and so have higher complementary value) and measure representativeness with increasing number of selected sites. We consider two scenarios when prioritizing CC sites: (1) representativeness of only networks of CCs and (2) representativeness of both CC and PA sites together. In the first scenario, we assess whether criteria can prioritize complementary sites across CCs. Here the ordination space represents only CCs (45 sites). In the second, we assess whether criteria prioritizing CC sites improve efficiency of existing PAs. Here criteria efficiency is assessed in reference to representativeness of CCs and PAs combined, and the ordination space represents both CCs and PAs (88 sites).
As input to ordination, species turnover was measured as between-site dissimilarities due to species replacement (Simpson dissimilarity calculated from species presence/absence data) (Baselga 2010). To define the two-dimensional space, we used multidimensional scaling (with 1000 iterations and 20 minimum random starts). The borders of this ordination space (or polytope) was defined using a convex hull (Faith and Walker 1996b). Within this polytope space, we placed hypothetical continuously spaced locations (or 'demand points'). The distance between all these hypothetical demand points and their nearest site was used to infer site complementarity. We placed twice as many demand points as sites included in the ordination (90 when ordination considers only CC sites, 176 when ordination considers both CC and PA sites).
To maximize representation, site selection should aim to minimise the distance between all hypothetical demand points and their nearest site, which represents the amount of 'forgone biodiversity' (Faith 1995). Sites were prioritized for inclusion in the selection process according to their criteria ranking (species richness, HSS, alien vegetation cover, or corridor width). As each site was added to the selection, we calculated the sum of the distance between each demand point and its nearest site. This distance was plotted against the number of selected sites.
To assess the complementary value of CCs to existing PAs, we considered the ordination space representing all sites (CCs and PAs), and compared random selection of CCs sites, PA sites, as well as CC and PA sites combined. Here we do not consider random site selection a worst-case scenario, but rather better than what would be expected under real world conditions, where conservation sites are rarely selected at random (Ferrier 2002;Delso et al. 2021). We also use random site selection as a benchmark to compare the ability of the four different criteria to select complementary sites. If a criterion performs better than random site selection, we consider this an indication of high efficiency. Where representativeness only considers CCs, random site selection of CC sites was the reference. Where representativeness considers both CCs and PAs, random selection of CC and PA sites combined was the reference. To this last scenario, we also considered random selection of only PA sites to assess criteria efficiency in reference to existing PAs.
For each of the two scenarios, the procedure described above was repeated 10 times, each time reprojecting the between-site dissimilarities in ordination space. Where criteria had similar values across different sites, we randomized the order of these same-value sites in each repetition. These analyses were repeated for all species, as well as subsets representing specialist and generalist species. Between-site dissimilarities in species replacement were calculated using the R package betapart (Baselga et al. 2021) and multidimensional scaling conducted using the vegan R package (Oksanen et al. 2019).

Sampling efficiency
In total, 36 species and 1175 individuals of adult male dragonflies were observed across the 88 sites. Table S1 provides a species abundance list. The species accumulation curve reached an asymptote ( Figure S3). Based on the final slope, the expected gain in diversity was low (0.003). The slope suggests that an additional ~ 333 dragonfly specimens would need to be recorded to increase the number of observed species.
We discarded from the data species with less than 10 observed individuals. These sparse observations could be due to low species detectability and we therefore consider them unreliable. Their removal also discouraged convergence issues during multidimensional scaling. Removal of these species did not significantly alter the ranking of sites according to species richness or the HSS. This reduced the number of species used in analyses to 20. Of these, nine were classified as specialists and 11 as generalists (Table S1). Figure 3 shows the curves tracing successive increase in species representation as more CC sites are included in the selection and compares random site selection of CCs, PAs, as well as CCs and PAs combined. In all cases, the ordination space (representativeness) considers species sampled from all sites (CCs and PAs). The graphs show that CCs complement Fig. 3 Forgone biodiversity with increasing number of randomly selected sites. Representation considered both conservation corridor (CC) and protected area (PA) sites. All-all species. Specialist-Specialist species. Generalist-Generalist species. The curve represents the gain in representation that corresponds to a decrease in forgone biodiversity as more sites are selected dragonfly diversity in PAs by contributing to the representation of generalist species, indicated by the more rapid decline in forgone biodiversity when CCs and PAs are considered together, compared to PA sites alone. Figure 4 compares efficiency of site prioritization according to the four criteria (species richness, HSS, alien vegetation cover, or corridor width) across two contexts: species present across CCs, and species present across both CCs and PAs.

Criteria efficiency
For all species, the efficiency of criteria for prioritizing complementary sites was generally comparable to random site selection (Figures S4 and S5). However, prioritization based on the HSS did not perform as well as random selection at representing species across CCs, and prioritization based on species richness performed more poorly than random selection at representing species across CCs and PAs combined.
For subsets of the species assemblage representing specialist and generalist species, the efficiency of criteria for prioritizing complementary sites was generally comparable to random site selection ( Figures S6 and S7). For sensitive species across CCs, site prioritization based on species richness or the HSS did not perform as well as random site selection at the start of the selection procedure.
Of the four criteria, corridor width was the best performer in prioritizing complementary CC sites (for all species, as well as specialist and generalist species separately), and occasionally performed better than random site selection. This was especially the case for generalist species when considering representation across both CC and PA sites. Fig. 4 Values of forgone biodiversity with increasing number of selected conservation corridor (CC) sites, for all species (All), specialist species (Specialist), and generalist species (Generalist). Two contexts are considered: species present across CCs, and species present across both CCs and protected areas (PAs). The curve represents the gain in representation that corresponds to a decrease in forgone biodiversity as more sites are selected according to four criteria: species richness (Spp), habitat quality based on the Habitat Sensitivity Sub-Index (HSS), habitat quality based on alien vegetation cover (AlienVeg), and corridor width (Width)

Discussion
We found CCs complement dragonfly assemblages in existing PAs through contributing to the representation of generalist species (i.e., species that are less responsive to habitat degradation). Of the criteria assessed, we found corridor width to be the most effective at prioritizing complementary sites, across either CCs or CCs and PAs combined. Furthermore, we did not find habitat quality to be more efficient at selecting complementary CC sites for specialist species compared to generalist species.
Studies have reported a global decrease in specialist species, which has been related to their higher sensitivity to degraded and fragmented habitats (Korkeamäki and Suhonen 2002;Colles et al. 2009;Morelli et al. 2020). These specialist species therefore tend to be located in less fragmented and less disturbed landscapes (Devictor et al. 2008). At the same time, specialist species can be more competitive in these natural environments (Büchi and Vuilleumier 2014), which could explain why generalist species are better represented across networks of CCs. The results show that by contributing to the representation of generalist species, CCs are aiding in the coexistence of sensitive and tolerant species across the region, thereby maximizing total species richness.
Of the criteria assessed, prioritization based on corridor width was the most efficient at selecting complementary CC sites, especially for generalist species. Research suggests that generalist species can benefit from habitat destruction, but they will perform better in landscapes of intermediate quality than landscapes of poor quality (Ramiadantsoa et al. 2018). Compared to narrow corridors, wide corridors experience less edge effects (van Schalkwyk et al. 2020), are better at maintaining connectivity (Hilty et al. 2006), contain better quality habitats (Pryke and Samways 2001), and experience less shading from alien trees, which is especially important for dragonflies in these CCs (Remsburg et al. 2008;Samways and Sharratt 2010). The higher conservation value of these wide corridors aligns with previous assessments based on site-level estimates of insect diversity (Pryke and Samways 2012;van Schalkwyk et al. 2017;Schalkwyk et al. 2020Schalkwyk et al. , 2021. This is the first work, however, to show the complementary value of wide corridors, both in relation to diversity represented in CCs and in relation to diversity represented in existing PAs. In general, prioritization criteria selected complementary sites that were comparable to random site selection, indicating unbiased species representation. Conservation site selection in the real-world is often biased and influenced by socio-political, economical, or aesthetic values (Kukkala and Moilanen 2013). Random site selection can therefore be expected to provide a good representation of a region's biodiversity and we did not consider it a worst-case scenario (Ferrier 2002).
Prioritization based on the HSS was less effective than random site selection when representing species across CCs, and prioritization based on species richness was less effective when representing species across both CCs and PAs. We only focused on the sensitivity sub-index of the Dragonfly Biotic Index as our focus was on habitat quality, but our results reflect other studies that also incorporated information on species distribution and threat status when calculating site scores, which found high levels of redundancy (the same species were repeatedly represented) (Simaika and Samways 2009). Species richness is often used for prioritizing conservation actions (Simaika et al. 2013) but it is also known to underperform at improving representativeness (Veach et al. 2017). Alternative ways of calculating species richness that shift the focus to rare or specialist species have been proposed (Astudillo-Scalia and de Albuquerque 2019; Jalkanen et al. 2020b). However, as we have shown, the value of CCs for increasing regional species richness lies in contributing to the representation of generalist species.
Variation in species detectability can distort results. We tried to control for this by removing the very rarest of species (< 10 observed individuals). The fact remains that rare species make up a large number of species in a community and could be rare for reasons other than low detectability. Future assessments should consider the influence of these rare species, while controlling for detection bias. Although our assessment suggests relatively good sampling completeness, we must bear in mind that this is only a temporal snapshot of the assemblages present across the whole region. Some sites were excluded from the analyses, as there were no dragonflies sampled at these sites. This implies that these sites have no value, and this is likely not true should more sampling be done over a longer time (Kietzka et al. 2021b) or if we also incorporated larval stages into our sampling scheme (Khelifa 2019;Kietzka et al. 2021c). We also did not consider the influence of inter-site geographical distance on criteria performance, which could be important considering the positive association between species turnover and geographical distance (i.e., distance decay) (Nekola and White 1999). Furthermore, species-level representation can be only one among many different prioritization goals (Lehtomäki and Moilanen 2013). Successful conservation plans need to go beyond only representing biodiversity to also ensuring biodiversity persistence by considering important processes, such as dispersal, species interactions, migration, and distributional range shifts in response to climate change (Sarkar et al. 2006;Economo 2011;Chase et al. 2020).

Conclusion and management implementations
Protected areas are often the focus of conservation assessments (Cantú-Salazar et al. 2013), but areas outside of PAs can contribute significantly to conservation. The potential of CCs for adding unrepresented features are increasingly being recognised in conservation plans, e.g., Schmid and Alonso (2005) and Jalkanen et al. (2020). Our results show that networks of CCs enhance regional representativeness by complementing dragonfly diversity represented in existing PAs. Wide conservation corridors with little alien vegetation cover provide the most complementarity for maximizing regional representativeness.