Community completeness as a measure of restoration success: multiple-study comparisons across ecosystems and ecological groups

When restoring habitat for biodiversity, the most effective outcome will be achieved by restoration projects which target several organism groups or ecosystem types. Such integrated approaches require direct comparisons among different ecological communities while evaluating success of restoration. The Community Completeness Index (CCI) is a recently developed metric that allows such comparisons by accounting for both present and absent but otherwise suitable taxa. We empirically evaluated the applicability of CCI for assessing the outcome of ecological restoration. We analyzed how species richness and the completeness of ecological communities recover after restoration, for different ecological groups and ecosystem types, and how it develops over time after restoration. Analyses were performed on 18 datasets with per site presence-absence data from Northern Europe. Each dataset represented one of the three habitat types (mire, forest, grassland) and different ecological groups (plants, flying insects, epigeic invertebrates). Datasets contained pristine, degraded and restored sites. We calculated the dark diversity and subsequently CCI based on species co-occurrences. Our multiple-study analyses revealed that CCI of grassland plant communities increased faster after restoration than invertebrate communities or plant communities in forests and mires. In addition, flying insect communities demonstrated significantly highest CCI in pristine mires. Some results were significant only for richness but not for CCI indicating species pool effect. Finally, completeness and species richness of restored communities increased with time since restoration. As such, our study demonstrated that CCI is a useful tool in evaluating restoration success across different organism groups and ecosystem types.


Introduction
The Anthropocene is characterized by rapid human-caused ecosystem loss and degradation (Lewis and Maslin 2015;Ellis 2018). To meet international targets on halting biodiversity loss, the importance of ecological restoration is widely emphasized (Suding 2011). In order to reveal if restoration targets are met, restoration actions require evaluations of their success, which is overall a key subject for the development of restoration ecology (Suding 2011). Nevertheless, how to evaluate restoration success is highly debatable (see e.g., Hobbs et al. 2009;McAlpine et al. 2016). The traditional use of past and present ecosystems as targets in restoration projects can be inappropriate due to rapid and uncertain environmental change (Corlett 2016). Furthermore, evaluations of restoration actions should optimally use Before-After-Control-Impact (BACI) study design and long term post-restoration monitoring (Christie et al. 2019). Due to high costs, these are rarely done. Therefore, success evaluation methods reflecting the resilience of restored ecosystem to future changes should be considered (Suding 2011;Corlett 2016).
It has been suggested that more effective outcomes of restoration for biodiversity should be achieved via integrated plant-animal (i.e., targeting several taxa) restoration approaches (McAlpine et al. 2016), and when whole landscapes are being restored (e.g., rewilding, Corlett 2016). Both approaches require comparisons across separate ecological communities while evaluating restoration success, either within the same ecosystem (across taxa comparisons) or within the same landscape (across habitat comparison of some representative taxa), respectively. Across taxa comparisons within focal ecosystem would show if the entire ecosystem is restored, not just its separate elements (Pärtel et al. 2013). In order to compare restoration success across taxa, standardised metrics that can easily be compared between taxa are needed.
Usually, traditional diversity measures (e.g., species richness) are limited in possibilities to perform direct comparisons across communities in different locations, as there is a variation in the size of the species pool. The habitat-specific species pool is defined as the set of species which potentially can inhabit a site (Zobel 2016). Community completeness (Pärtel et al. 2013) can be used for direct comparisons, as it standardizes local richness to available species pool. The community completeness is describing the proportion of habitat-specific species pool which is present in the site, i.e. explores local diversity which is mathematically independent of species pool size. The rationale behind using community completeness for assessing restoration success is that the aim of any restoration effort is to achieve a biological community that contains as high proportion of the habitat-specific species pool as possible. The species pool for a given site includes two species diversity measures: observed diversity (habitat-specific species actually present in a site) and dark diversity (the absent portion of habitat-specific species pool (Pärtel et al. 2011)). While calculation of observed diversity is straightforward, dark diversity cannot be measured directly and can only be estimated, e.g. from species co-occurrence patterns (Lewis et al. 2016). These two measures can be used to calculate the relatively recently proposed Community Completeness Index (CCI), expressed as ln(observed diversity/dark diversity) (Pärtel et al. 2013;Lewis et al. 2016). CCI indicates whether the focal community is becoming more (higher CCI) or less (lower CCI) complete following restoration activities. As a relative measure, CCI is comparable among different communities. If restoration is successful, we might expect communities to become more complete with time, thus revealing higher ecosystem or landscape resilience. However, this remains to be evaluated empirically.
The properties of the CCI suggests that it could be used in formulating clear goals of most ecological restoration projects and refer restoration targets to the mechanisms of larger processes (Török and Helm 2017), but it was not tested empirically. Restoration projects should aim to achieve higher CCI by encouraging all habitat-specific species available from the regional species pool to disperse and remain in the communities after restoration actions have been implemented. Thus, changes in CCI reflect whether desirable and available species are colonizing restored sites. In contrast, traditional diversity measures such as richness do not consider species identity and thus do not necessarily reveal if habitat was really improved for the target habitat-specialist species. For example, higher richness could be achieved merely by influx of generalist and invasive species, as well as specialists. Post-restoration community completeness should depend on spatial and temporal constraints, i.e. site connectivity and species dispersal abilities, and the time since restoration. Indeed, plant and hoverfly species with poor dispersal abilities tend to be part of the dark diversity more often than good dispersers (Riibak et al. 2015;Moeslund et al. 2017;Miličić et al. 2020). Further, communities should become more complete with increasing time since restoration actions, as potentially more species should continuously arrive into the restored sites under successful restoration (e.g. ''Field of Dreams'' hypothesis, Palmer et al. 1997). Longer time should be required to reach higher postrestoration completeness values for the communities dominated by poor dispersing species than those formed mostly of good dispersers. The possibility to restore an ecosystem might also depend on to what extent species are adapted to disturbances. Recovery from major disturbance should be easier for ecosystems which experience continuous but low intensity disturbances than for those rarely but intensively or almost never disturbed (see e.g., Betts et al. 2019). Finally, communities of different trophic levels could recover based on their trophic position in the food web, thus reflecting the sequence of species arrival after restoration: carnivores are dependent on herbivores, while herbivores are dependent on autotrophs (Lindeman 1942).
The aims of our study were to (a) evaluate the applicability of the CCI for assessing the outcome of restoration actions, and (b) evaluate how the completeness of ecological communities recovers after restoration, for different ecological groups and ecosystem types and how CCI recovery differs from the recovery of species richness. We hypothesize that: (1) Species richness will show similar but non-identical patterns of responses to restoration in comparison to CCI (due to the lack of standardization to local species pool size). (2) Restoration increases community completeness, such that CCI values are greater for restored than degraded communities and pristine communities are expected to have the highest completeness, and that community completeness in restored sites increases with time since restoration. (3) Across ecosystem comparisons should show that communities dominated by better dispersers and experiencing frequent, but low intensity disturbance (semi-natural grasslands), have higher completeness after restoration than those dominated by poor dispersers and affected rarely by severe disturbance (forests) and the lowest completeness should be seen for those almost never disturbed communities (mires). (4) For across ecological group comparisons, we formulate two different hypotheses for the completeness of investigated groups, depending on whether mobility or trophic position is more important: (i) if mobility is more important, we would expect that CCI values from highest to lowest would follow the order: flying insects, epigeic insects and plants; (ii) if trophic level is more important, we would expect CCI to be ordered from highest to lowest as plants, herbivores and predators.

Data
On November 2018 we searched for ecological restoration studies from the region of Northern Europe by entering different combinations of three keywords into Google Scholar search field: restor* AND habitat type [= grassland, meadow, forest, woodland, mire, bog] AND country [= Finland, Sweden, Norway, Denmark, Netherlands, Belgium, Estonia, Latvia, Lithuania]. We only searched for studies published in 2004 or later. The reason was both that we considered possible differences in ecological conditions due to global climate change among studies over a longer period, and that we expected it to be difficult to get access the full data for studies older than * 15 years. We restricted our search to the northern part of Europe (Nordic-Baltic countries and Benelux), since this region would share similar climatic conditions and biogeographical history, thus resulting in similar habitat-specific species pools. After the search we selected publications which included ecological restoration datasets with such minimal requirements: species presence-absence data of plants, insects or spiders, one of three habitat types (semi-natural grasslands, mires, forests), three (pristine, restored, degraded) or two (pristine and restored, or degraded and restored) treatments, at least 20 total replications in a study system (= sites or sites*treatment), and data on time since restoration, i.e. the time between restoration action were performed and the time of the evaluation, in years (if available). We contacted authors asking them to provide species presence-absence data per site, applied treatment for each site and indication of time since restoration for restored sites. This resulted in 18 datasets with per site data from four countries (Table 1). 16 datasets contained data of time since restoration in years for restored sites. We organized all datasets so that each site would represent one of the three treatments (space-for-time substitution): pristine, degraded and restored. Most of the datasets included all three treatments, but some had only either pristine and restored or degraded and restored treatments. One paper (Marozas et al. 2007) did not frame the study in a restoration context, but compared burned and unburned managed forest sites. As burning is a common practice in forest biodiversity restoration projects (see e.g. Kärvemo et al. 2017), we considered this a restoration treatment and included these data into our restoration analyses.
Datasets were grouped for comparison across ecosystems and across ecological organism groups. We grouped the data sets on plant communities into three ecosystem types which reflects different intensity and regularity of disturbances and thus the respective communities are dominated by good or poor dispersing species: (1) semi-natural grasslands (frequent low intensity disturbance, good dispersing plants), (2) forests (rare high intensity disturbance, poor dispersing plants), (3) mires (no regular disturbance, poor dispersing plants). The same grouping was also done for flying insects, but it covered only two ecosystem types: grasslands and mires. For ecological group comparisons we grouped our datasets within grasslands and mires separately based on organism movement and dispersal abilities: (1) flying insects (highly mobile, good dispersers; it included bees, hoverflies, craneflies, lepidopterans), (2) epigeic invertebrates (intermediate dispersal abilities; incl. ants, carabid beetles, spiders), (3) plants (sessile organisms, limited dispersal abilities; incl. vascular plants, bryophytes). We are aware that such grouping might be oversimplistic from the perspective of individual species (as e.g. the seeds of some plant species can be carried long distances by wind), but here we consider the dominant mode of *Numbers refer to respective analyses: 1-treatment effects on plant richness and CCI across grasslands, mires and forests (10 datasets), 2-treatment effects on flying insect richness and CCI across grasslands and mires (5 datasets), 3-treatment effects on richness and CCI across communities of plants, epigeic invertebrates and flying insects in mires (7 datasets), 4-treatment effects on richness and CCI across communities of plants and flying insects in grasslands (9 datasets), 5-overall time since restoration effect on richness and CCI (16 datasets) dispersal of most species forming specific assemblage and assume that active behaviordriven animal dispersal is more effective than passive plant dispersal from habitat restoration perspective. At least partly, the same grouping could be seen from trophic level perspective: autotrophs or producers (= plants), primary consumers (= flying insects; majority of species, but not necessarily all as, e.g., some hoverflies are predators in larval stage), secondary consumers (= epigeic invertebrates; again, appropriate for most of the species, but some are herbivores, e.g. some carabid beetle species).

Calculations of dark diversity, completeness and species richness
From the data provided by the authors, we calculated habitat-specific species pools and dark diversity for each site within each dataset by employing the co-occurrence method proposed by Lewis et al. (2016) and available online at https://shiny.botany.ut.ee/DarkDiv/ . The general idea of the method is that present species are used as indicators for all absent species (targets): those absent species which elsewhere likely co-occur with many present species are selected to the dark diversity. Calculations are based on Beals index which determines the proportion of indicator species frequency where it co-occurs together with the target species (ratio ranging between 0 and 1). For each target species these proportions are averaged over all indicator species present in the community. The Beals index depends mathematically on the frequency of target species (common species have higher values) and therefore species specific thresholds are used. For the thresholds we calculated Beals index for all species across all sites. The threshold for a species is determined from the Beals index values at sites where the species is actually present by selecting a typical low value needed to have the species present. We used a narrow definition of dark diversity which includes fewer species but with more confidence (i.e. the 10th percentile of the index found in the sites where the species was actually present; the broad definition uses the 1st percentile). In order to include as much co-occurrence information as possible, we used the sites of all treatments within each dataset for determining co-occurrences. As such, dark diversity was determined by a statistical technique, rather than relying on subjective decisions on whether species should belong to particular habitat-specific species pool or not.
We calculated community completeness for each site within each dataset as follows: Completeness = ln(observed diversity/dark diversity). Such log-ratio serves as a logistic transformation of percentages or proportions. Therefore, it is not bounded between 0 and 100% (or 0 and 1) and could be used in statistics that expect unlimited sample space (e.g., mean calculations, models). If site's dark diversity was equal to zero, we assigned 0.5 value instead to make completeness calculations possible.
As a measure of species richness, we used the observed number of species per site. Species richness was log-transformed before the analyses to keep it on the same scale as completeness.

Statistical analyses
In order to evaluate restoration treatment and time since restoration effects on community completeness and species richness, we ran linear-mixed models (lme function from package nlme (Pinheiro et al. 2020)) in R version 3.6 (R Core Team 2019). Schematic overview of all performed analyses is given in Fig. 1. Two alternative response variables were community completeness and log-species richness. Sites nested within each dataset were included as a random effect.
The treatment analyses were performed by making comparisons across habitats (for separate ecological groups) and across ecological groups (within separate habitat types). Treatment effect models contained three treatments (pristine, restored, degraded) and treatment interaction with ecosystem type (in across ecosystem comparisons) or ecological group (in across ecological group comparisons). At least two datasets for specific category (ecosystem type or ecological group) were required to include it in the respective analyses.
Models testing the effect of time since restoration on community completeness and species richness included only restored site data from the respective datasets. In these models, the actual number of years between the restoration and biodiversity survey served as an explanatory variable (one year per site). These analyses were performed at two levels: (1) Overall time since restoration effect on completeness and species richness of communities. Here 16 datasets were included with completeness and species richness as alternative response variables and time since restoration as an explanatory variable. There was no data grouping by ecological groups or habitat types. (2) Time since restoration effect on community completeness and species richness of ecological groups within particular habitat type (models with time since restoration*ecological group interaction) and on CCI and richness of plant communities across habitat types (models with time since Fig. 1 Schematic overview of performed analyses. Ecosystem treatment and time since restoration (for restored sites) effects on community indices (CCI and Species Richness) were evaluated by employing linear-mixed models. Separate analyses are numbered (1-5) and corresponds to the numbers of analyses in Table 1. The number of datasets used for each separate analysis is written in square brackets restoration*habitat type interaction). These analyses aimed to reveal if time since restoration effect depends on ecological group or habitat type.

Treatment effects on species richness and community completeness of different habitat types
We had enough replicates (= datasets) per habitat type to perform analyses on treatment effects across habitat types on species richness and CCI of plants and flying insects (Table 2; Fig. 2).
Treatment effects on plant species richness and community completeness across grasslands, mires and forests (10 datasets, Table 1). Plant species richness was highest in pristine and restored grasslands in comparison to any other differently treated habitat type, demonstrated by two significant interactions: pristine*grassland (overall highest richness) and restored*grassland (Table 2; Fig. 2a). Similarly, the completeness was highest for plants in restored grasslands, with lower levels of completeness in restored forests and mires, and in degraded sites of all three habitat types (Fig. 2b), indicated by a positive restored*grassland interaction (Table 2).
Treatment effects on flying insect species richness and community completeness across grasslands and mires (5 datasets, Table 1). Flying insect species richness were significantly lower in mires than in grasslands (Table 2; Fig. 2c). Moreover, two interactions were positively significant: pristine*mire and restored*mire, indicating highest flying insect species richness in pristine sites and positive restoration effect on it. Only the interaction pristine*mire were significantly positive for community completeness of flying insects (Table 2; Fig. 2d).

Treatment effects on species richness and community completeness of different ecological groups
We had enough replicates (= datasets) per ecological group to perform analyses on treatment effects across ecological groups within mires and grasslands (Table 3; Fig. 3).
Treatment effects on species richness and community completeness across communities of plants, epigeic invertebrates and flying insects in mires (7 datasets, Table 1). Species richness across ecological groups in mires showed several significant effects (Table 3; Fig. 3a). Pristine mire communities were associated with significantly lower species richness. However, both plant and flying insect richness was significantly positively associated with pristine sites (significant Pristine*Flying insects and Pristine*Plants interactions). Moreover, flying insect richness in mires were positively affected by restoration (significant Restored*Flying insects interaction). In terms of community completeness across ecological groups in mires, only interaction Pristine*Flying insects were nearly significant and positive (Fig. 3b).
Treatment effects on species richness and community completeness across communities of plants and flying insects in grasslands (9 datasets, Table 1). Grassland plant species  richness demonstrated significant Pristine*Plants and Restored*Plants interactions (Table 3; Fig. 3c). However, overall species richness in restored grasslands were significantly lower. Grassland plant community completeness was highest for restored sites (significant interaction Restored*Plants), followed by pristine sites (significant interaction Pristine*Plants) ( Table 2; Fig. 3d).

Time since restoration effect on species richness and community completeness
The analyses on overall time since restoration effect (16 datasets, Table 1) revealed a positive effect of time since restoration on community completeness (est.=0.017, p = 0.031, Fig. 4a), as well as on species richness (est.=0.016, p = 0.001, Fig. 4b). Thus, generally, restored communities become more species rich and more complete (gain more species from available species pool) with increasing time since restoration. Unexpectedly, this does not depend on ecological group nor habitat type, as the analyses of time since restoration effect on ecological groups within habitats and on habitats within ecological groups did not provide significant interactions (time since restoration*ecological group and time since restoration*habitat type, respectively).

Discussion
The applicability of Community Completeness Index (CCI) in evaluating success of restoration actions Assessment of restoration success is a key subject in restoration ecology (Wortley et al. 2013;Prach et al. 2019), which requires clearly stated and justified goals (Gann et al. 2019;Prach et al. 2019). As such, restoration targets should include measurable indicators which contain mechanistic relationships to larger ecosystem structure and processes (Prach et al. 2019). Our multiple-study comparisons revealed that community completeness could be highly relevant in formulating goals of restoration projects, as well as it could serve as useful indicator while evaluating their success. We demonstrated that CCI could be used in making comparisons of restoration treatments across ecological groups (plants, epigeic invertebrates and flying insects) and across habitats (grasslands, mires and forests). When comparing analyses with CCI and traditional species richness, there were many similarities. This is expected since local richness is the product of species pool size and community completeness. However, some results were significant only for richness, but not for CCI. These differences likely originate from variation in the species pool. For example, species richness might increase, but completeness remains the same if there are also more species in the dark diversity of specific habitat. In our study, species pool Higher or lower CCI values indicate that, respectively, more or less species from habitat-specific species pool is realized within a site variation effect seem to be especially associated with flying insect communities in restored mires. Many significant results were found for species richness, but not for CCI in across habitat type analyses for flying insects and in across ecological group analyses in mires. Likely, conditions in restored mire sites have changed so much for flying insects that species pools are different from pristine, indicating high specialisation of pristine mire insects and preference for stable habitat conditions. In contrast, the results which were evident both for richness and CCI seem to be more associated with plants in restored and pristine grasslands. These indicate directly the effect of local site-level process (restoration), as species pool effect was filtered out. Grassland plant communities seem to be improved by the restoration, as they prefer continuous disturbance. Finally, none of the site-level effects were completely masked by the species pool variation, as no results were statistically significant for CCI only. CCI based restoration targets depend on the chosen principles of species pool and observed/dark diversity estimations. Here we used species presence-absence data from all sites (pristine, restored, degraded) within a study system to define estimates for CCI. In this case, any changes in CCI show if available species generally are being gained or lost due to restoration and, importantly, still allow across-taxa comparisons. However, this approach does not give information on whether restored habitats are approaching a specific target community, for example a historic or pristine state, as has traditionally been the main purpose of restoration actions (Jackson and Hobbs 2009). If the pristine state is a specific target, it would be more relevant to calculate CCI from species co-occurrence patterns in the pristine sites only or from other way defined habitat specific species pool of the region. However, it is increasingly recognized that in the era of rapid human-caused environmental changes, finding a suitable reference sites could be difficult or impossible because habitat degradation can be so widespread that none of the desired habitats exist (Murphy 2018;Prach et al. 2019) and ecosystems are changing continuously due to anthropogenic disturbances. As a result of these changes, restored ecosystems can turn into some alternative stable states with novel mixtures of species, thus resulting in hybrid or novel ecosystems ). While evaluating restoration success, comparisons of restoration interventions against the counterfactual, rather than comparing against some fixed state, should be preferred (Corlett 2016). However, if the dynamic nature of ecosystems is considered while evaluating restoration success, then it appears problematic to define the desired restoration target. For example, which species should be considered as part of novel or hybrid ecosystem species pool and thus used in its CCI calculations? Clearly, actual species preferences to co-exist within whole study system should be used in such cases.
Despite differences in gaining estimates for CCI, a clear restoration goal could be to encourage all habitat-specific species available from the regional species pool to disperse and remain in the communities changed by restoration actions (= achieve higher CCI due to restoration). As such, CCI creates a link with mechanisms of larger processes: it reflects if suitable (or desirable) and available (from the regional pool) species are re-colonizing restored sites. Finally, as demonstrated in the current study, CCI allows comparisons across ecological group and across habitat types in accordance with restoration actions. Such comparisons are highly useful in both large and small-scale restoration projects. When entire landscapes are restored, then across-habitat comparisons would reveal if separate elements of the landscape are equally successfully restored. Similarly, when single habitat is restored, then across ecological group comparisons reveal if different elements of the whole ecosystem are being restored (Pärtel et al. 2013). In case some landscape or ecosystem elements are less successfully restored, the respective adaptive management techniques addressing specific elements could be applied.
Even though dark diversity and community completeness serve as promising tools in restoration ecology, they are relatively novel concepts (Pärtel et al. 2011) and thus their methodology is still under development (Brown et al. 2019;Carmona et al. 2019). For CCI estimations all species present within site must be recorded, or at least the possible bias of non-detected species should be the same throughout the study system. Thus, such bias is likely avoided when CCI is used for across habitat comparisons. However, it could be more difficult to avoid sampling bias in across ecological group comparisons. Differences in detectability and mobility between taxa can result in higher estimated completeness e.g. in plants which are sessile, and highly mobile animal (e.g., flying insect) communities. Also, flying insects might use several habitats during their life cycle (e.g., for feeding and nesting, Potts et al. 2005;Ö ckinger et al. 2018), thus the observed community composition is more likely to be just a snapshot of particular moment. However, at least bias in species detectability could be reduced by using environmental DNA as species detection method (see e.g. Boussarie et al. 2018).

Positive restoration effects on CCI: different ecological groups and ecosystem types, main restricting factors
Restoration actions had a clear positive overall effect on completeness of communities, as CCI increased with the time since restoration. Furthermore, somewhat surprisingly, grassland plant communities had higher CCI in restored than in pristine sites. This could be linked with preference for disturbance of plant species in such habitats. In contrast, flying insect greatest completeness was associated with pristine mires, thus indicating their preference for very stable conditions. Likely, high CCI in restored sites indicates the dynamic nature of ecosystems, as our calculations of community completeness do not rely purely on species pools from reference sites, but instead consider general species cooccurrence patterns. As such, it is possible that hybrid ecosystems with mixed species pools (characteristic to both degraded and pristine state) were created by restoration actions (at least for plants in grasslands). Under this scenario, even after long time CCI should remain high for restored communities. Alternatively, it could be that poor competitor species with good dispersal abilities (r-strategists) arrived first to the restored sites (e.g., Woodcock et al. 2012), which temporally results in communities with high CCI values. However, in this case, K-strategists should outcompete r-strategists with time and the restored communities should become more similar in terms of their completeness values with those of pristine sites.
Indeed, post-restoration community completeness should depend on both spatial and temporal constraints, e.g., site connectivity, species dispersal abilities and the time since restoration, as these are among the most important factors for re-establishment after restoration (Jones et al. 2018;Rotchés-Ribalta et al. 2018). Species with poor dispersal abilities are also known to be part of dark diversity more often (Riibak et al. 2015;Moeslund et al. 2017;Miličić et al. 2020). In our study, across-habitat comparisons of plant communities indicated that regularly but with low intensity disturbed, and thus dominated by better dispersing species, communities (grasslands) became more ''complete'' faster than communities rarely but severely (forests) or almost never (mires) disturbed. However, across ecological group comparison clearly indicated that plant communities respond faster to restoration than better dispersing invertebrates (significantly positive restored*plants interaction in grasslands). Although this result could be linked to dispersal abilities (as plants are sessile, they need to respond immediately to environmental changes in a site), but more likely this is an indication of invertebrate dependency on plants. First autotrophs (lower trophic level in a foodweb) should establish after restoration, while heterotrophs (= higher trophic level) can respond only after changes in lower trophic level (e.g., Henson et al. 2009). This is also supported by the fact that greater CCI of flying insect communities was associated with pristine mires (i.e., with very stable conditions). Furthermore, fast plant response could be a result of seed bank effect: many plants could be present in degraded sites in dormant seed bank or as small remnant populations (Eriksson 1996; Kalamees et al. 2012;Vandvik et al. 2016), thus allowing plants to establish immediately after restoration. In contrast, invertebrates need to recolonize restored sites from nearby source populations (Ö ckinger et al. 2018).

Conclusions
Our study revealed that CCI derived from species co-occurrences could be used in evaluating restoration success if its target is to increase the number of ecologically suitable species present in the habitat-specific species pool, in the restored sites. In such case, CCI would indicate if restoration overall increases diversity and would also allow direct comparisons across separate communities. The usefulness of CCI estimated from species co-occurrences is especially clear when reference sites are not available, for example, when restoration actions result in hybrid or novel ecosystems. However, this method is more limited when the purpose of restoration is to get back to some desired ''pristine'' state. Here estimates for CCI (habitat-specific species pool, dark/observed diversities) from co-occurrences must be derived by only using reference sites. In practice this means that higher number of pristine sites should be sampled in order to derive reliable estimates.