A general framework for predicting delayed responses of ecological communities to habitat loss

Although biodiversity crisis at different spatial scales has been well recognised, the phenomena of extinction debt and immigration credit at a crossing-scale context are, at best, unclear. Based on two community patterns, regional species abundance distribution (SAD) and spatial abundance distribution (SAAD), Kitzes and Harte (2015) presented a macroecological framework for predicting post-disturbance delayed extinction patterns in the entire ecological community. In this study, we further expand this basic framework to predict diverse time-lagged effects of habitat destruction on local communities. Specifically, our generalisation of KH’s model could address the questions that could not be answered previously: (1) How many species are subjected to delayed extinction in a local community when habitat is destructed in other areas? (2) How do rare or endemic species contribute to extinction debt or immigration credit of the local community? (3) How will species differ between two local areas? From the demonstrations using two SAD models (single-parameter lognormal and logseries), the predicted patterns of the debt, credit, and change in the fraction of unique species can vary, but with consistencies and depending on several factors. The general framework deepens the understanding of the theoretical effects of habitat loss on community dynamic patterns in local samples.

corresponding empirical values are rarely available for use. Using Diamond's biogeographic kinetic model 10 as an example, the relaxation parameter is key to predict time-dependent extinction of species; however, thus far, few empirical values have been reported for different taxa. Of course, no models are perfect; KH's model can only be applied to predict extinction debt in the whole remaining intact habitat when the original model is employed. Moreover, the assumptions made by HK's model are too restricted. Regional SAD may shift after habitat destruction, in particular when the destruction is very intense (for simplicity, this relaxation has not been investigated here. We will show that most post-disturbance patterns could be predicted either from fixed SAD or from shifted SAAD scenarios).
It is clear that KH's original method only measures the change of regional richness (or gamma diversity). This information is insufficient to provide a complete understanding of the consequences of ecological communities in response to habitat loss. In the local field, it is often observed that local diversity can either increase, decrease, or remain unchanged after habitat destruction and can depend on sampling locations as well 13 , which may either remain consistent or contradictory to the temporal pattern of regional diversity change 14,15 . How can one resolve this local-regional diversity change paradox? Moreover, after habitat loss, either biotic homogeneity (fraction of unique species declines) or heterogeneity (fraction of unique species increases) may emerge [16][17][18] . In fact, a growing number of previous studies 13,15,19 have empirically documented these diverse, but distinct, ecological change patterns. However, a parsimonious, tractable, and synthetic framework that can link and explain these temporal ecological phenomena is unavailable until date. Other than regional species richness, KH's original method is unable to quantify any other diversity patterns directly, in particular for community structure (local diversity change and the changes in the fraction of unique species). However, by using simple modifications and relaxing some assumptions (SAAD can vary before and after habitat loss), we can demonstrate that KH's model can be generalised to predict very diverse temporal community change patterns across different spatial scales. By doing so, the generalised statistical framework can perfectly resolve the local-regional diversity change paradox and explain the reason behind the increase or decrease in the local species richness or the fraction of unique species between two local communities after habitat loss.
In the present study, we address the following questions: (1) How many species are subjected to delayed extinction in an intact local community during habitat destruction in other areas? (2) How many rare or endemic species will go extinct or immigrate into the focused local community and how do they contribute to the overall magnitude of extinction debt or immigration credit of the whole local community? (3) How will community structure change after habitat destruction? By answering these questions theoretically using the proposed framework, we can conduct diverse predictions and obtain a relatively complete picture on the delayed responses of ecological communities caused by habitat loss. The general framework takes a big step forward to forecast the temporal biodiversity change and project community structure alteration. The framework can not only assess the effects of habitat destruction but is also identically applicable to evaluate the effects of climate change.

Materials and Methods
Predicting delayed species loss and gain for a single local community. Suppose there is a region with area size A 0 , species number S 0 , and total individual number of all the species N 0 . Further, there is a local community with area size A 1 that is nested within the region ( ⊂ A A 1 0 ) (Fig. 1). Then, the number of species found in the local community before habitat loss is given by the product of SAD and SAAD as follows 9, 20, 21 : where φ θ | n ( ) A 0 represents regional SAD with parameter θ A 0 estimated for the whole region A 0 , and ω P n A A (0 , , , represents SAAD with parameter ω. Accordingly, the number of endemic species that are only found in the community is given by , SAAD denotes the probability of observing m individuals for a species in the local community A 1 , which has a total of n individuals in the whole region A 0 , which will be explored in detail later. In particular, φ θ | n ( ) A 0 here is not a standard probability density function, because . However, if one utilises a standard probability density function φ θ | n ( ) Suppose there is also another local community with area size A t 1 in the region ( ), which will be completely destructed at some time point. The observed species richness in the intact local community A 1 is thus given by Eq. (1) immediately after habitat loss. However, after a long run since the irreversible loss of A t 1 , the local community A 1 (in the remaining intact region denoted by A A / t 0 1 ) will reach a new equilibrium of species richness, which can be estimated as follows: and ω | P n A A A (0 , , / , ) denote the re-fitted SAD and SAAD when the community reaches a new equilibrium after habitat destruction. If SAD, SAAD, or both are shifted to new probability functions with new parameters (these scenarios are also investigated later in this work), it is more appropriate to write the equation as φ θ , respectively. Accordingly, the predicted extinction debt for the local community A 1 in the remaining habitat can be computed as follows: Scientific RepoRts | 7: 998 | DOI:10.1038/s41598-017-01070-2 , immigration credit occurs. When the area size of the local community is large to include the whole remaining habitat , the above equation actually becomes KH's original model.
Delayed loss and gain of endemic species. The loss or gain of the endemic species in the local community A 1 is the part of the magnitude of extinction debt or immigration credit contributed by endemic species that is only found in the community. Based on this definition, their debt or credit magnitude is calculated using Eq.
(2) as follows: Delayed loss and gain of species with small population sizes. After the loss of A t 1 , part of the extinction debt or immigration credit derived from species with small population sizes in the local community A 1 can be calculated. For a given small-population threshold n c , the delayed loss or gain of species with population size n c found in the local community A 1 is calculated as follows: Model configuration for predicting extinction debt in a single local community (A) or changes in the fraction of unique species for two non-overlapped local communities (B) based on two community patterns, namely regional species abundance distribution and species spatial distribution. Grey area A t 1 denotes the destructed area of the region. Species loss and diversity in the circled areas A 1 and A 2 are the intact local communities in the remaining region without habitat destruction. The whole region is not close but open to immigration and speciation. ) in the intact region ( ), we can compute the changes in the area-based fraction of unique species number for two non-overlapped local communities after habitat loss as follows: Two SAD models as demonstrations. Canonical lognormal distribution. Suppose the regional SAD for the whole region A 0 is represented by Preston's canonical lognormal distribution 9,22 , may be used to estimate the parameter ∆ A 0 . However, in practice, the latter equality overwhelmingly determines the numerical estimation of the parameter.
Furthermore, in our study, SAAD is supposed to follow finite negative binomial distribution (NBD) 23 as follows: where k represents the aggregation degree of species distribution. The above probability model is flexible, because it reaches random distribution as k → ∞ and maximal aggregation of species distribution as k → 0. We assess and compare the effects of spatial distribution patterns of species on the estimation of extinction debt and immigration credit. Therefore, we substitute k = 1 and k = 100 in the above equation to represent distributional aggregation and randomness of species, respectively. We also report the comparative results for two more extreme values k = 0.001 and k = 1000, representing extreme aggregation and randomness of spatial distribution, respectively.
For example, when species distribution is aggregate (k = 1), delayed species loss and gain can be estimated as follows: is numerically determined by the following equation, The last equality in Eq. (11) is conditioned on the fact that the total number of individuals species is proportional to the area size, which is a widely applied prerequisite in many macro-ecological studies 9,24,25 .
For lognormal SAD model, KH presented an asymptotic formula for estimating species number after the community reaches a new equilibrium since the point of habitat destruction as = . This is a classic power-law species-area relation with a widely used slope 0.25 for a log transformation. This formula may be used for large Δ, but for guaranteeing computational accuracy, we numerically estimate the parameter Δ in φ ∆ n ( ).
Scientific RepoRts | 7: 998 | DOI:10.1038/s41598-017-01070-2 Single-parameter logseries distribution. Fisher's logseries distribution was a widely used SAD model in ecology, because it can predict a high volume of rare species in the ecological community, which has a form as follows 26 : However, because it contains two unknown parameters α and x, it is not suitable for estimating the number of species after habitat loss based on one observation on the number of individuals. Whereas, it is possible if one parameter (e.g. α here) is assumed to be constant before and after habitat loss and the estimation of species number after habitat loss can be achieved, as in KH's study. Based on this assumption, the loss of species richness is proportional to the habitat loss.
The same issue is applied to other variant logseries models, such as the one derived from the maximum entropy theory 24,27,28 , which is given by When N 0 goes to infinity and N 0 /S 0 approaches to a constant, Eq. (12b) becomes very similar to Fisher's original version (Eq. 12a). Since N 0 /S 0 can be observed, Eq. (12b) contains only one unknown parameter λ A 0 , which can be estimated before habitat loss. However, after habitat loss, Eq. (12b) would contain two unknown parameters ), but we only have the number of individuals after habitat loss in hand. Thus, it is not directly applicable to the present model unless some assumptions have been made. For example, like the discussion above, the parameter λ A 0 can be assumed to be constant (so ) or linearly related to the percent of lost habitat.
To break through the limit, a promising way is to introduce another variant logseries model containing only a single unknown parameter. Thus, it is applicable to predict the number of species after habitat loss. It is possible to extend the logseries model developed by R.A. Fisher, in which he refers the original model to another widely used probability model in ecology, namely the ordinary NBD. The ordinary NBD actually was used in modelling species abundance and later introduced to describe species distributional patterns (aggregation, random, or regular) 23,26,29 . However, ordinary NBD-driven logseries model (Eq. 12a) was not a proper probability model and degenerate or undefined at the zero point. To overcome this problem, we introduced a truncated NBD model, and used that to induce a single-parameter logseries model as follows: where γ A 0 is the unknown parameter, and C 1 is a normalisation factor to make . Details of the derivation of this model could be found in the Supplementary information. However, we did not fit γ A 0 directly; instead, as shown in the Supplementary information, it is easy to show that it is area-dependent, i.e.
, in which ω is a hyper-parameter that is assumed to be a constant. The parameter γ A 0 (actually we estimate the hyper-parameter ω here) can be numerically estimated by solving the following equality . As a result, the new species number at equilibrium after habitat loss can be estimated from Eq. (13). Apparently, the parameters γ A 0 and γ A A are unequal (though the hyper-parameter ω is fixed), strictly following our relaxed assumption that parameters can vary before and after habitat loss.
In our study, we utilised the above two regional SADs (Eqs 8 and 12c) and the SAAD model, which follows the finite negative binomial distribution, to demonstrate the effects of habitat destruction on local intact communities in the remaining habitats. We specifically investigated the effects of habitat loss on the delayed species gain or loss of a single local community and tested whether endemic and rare species may contribute to the gain or loss of species. We also investigated the effects of habitat loss on the alteration of community structure (quantified as the changes in the fraction of unique species between two random local communities). Moreover, we tested the contribution of endemic or small-population species to the loss and gain of biodiversity. For these analyses, we compared two distinct forms of species distribution, namely randomness (k = 100 and 1000, respectively) and aggregation (k = 1 and 0.001, respectively). We assessed whether shifting species distribution pattern (difference in distributional aggregation or randomness degree of species before and after habitat destruction) result in new species and community dynamic patterns in local communities.

Results
Delayed biodiversity loss and biodiversity accrual in a single local community. Regardless of the regional SAD and SAAD models being used (as long as they are always fixed before and after habitat loss), an intact local community A 1 with a larger area size would have a higher magnitude of extinction debt for a given percent of habitat loss ( Fig. 2A,B). For both SADs, an increase in the habitat loss would lead to higher extinction debts in the intact local community ( Fig. 2A,B). Moreover, the debt becomes higher when species distribution is less aggregate (Fig. 2A,B).

Figure 2.
Prediction of extinction debt or immigration credit for a single local community and changes in the fraction of unique species for a pair of local communities in the remaining intact region. All scenarios here assume that both SAD and SAAD are fixed before and after habitat loss. Species distribution is assumed to be random (k = 100) or aggregate (k = 1). (A,B) Extinction debt or immigration credit as a function of the area size of a single local community; (C,D) extinction debt or immigration credit for endemic species in a single local community; (E,F) extinction debt or immigration credit contributed by small-population species (the area of the local community is fixed to have a fraction 0.4 over the region); and (G,H) changes in the fraction of unique species as a function of the area size of one local community in a pair (the area size of the other community is fixed to be 10% of the whole region). Before habitat destruction, the region was assumed to have 30 species and 200 individuals in the model. In all the subplots, dashed lines represent the outcomes when 50% of the region is destructed, while solid lines represent the outcomes when only 10% of the region is destructed. Interestingly, no extinction debts for endemic species are observed in both the SAD models (Fig. 2C,D). More endemic species are found in the local community (thus immigration credit) when species distribution becomes aggregate and habitat destruction degree is high (Fig. 2C,D). By contrast, for species with different population sizes (not necessarily endemic), when the population threshold is small and only rare species are included, the immigration credit is observed in both the SAD models (Fig. 2E,F). When more species with larger or intermediate population sizes are taken into account (by setting a higher population threshold), immigration credits disappear, while extinction debts emerge, regardless of species distribution pattern and habitat loss degree (Fig. 2E,F).

Delayed changes in the fraction of unique species for two local communities.
There are some interesting patterns for the changes in the fraction of unique species between two local communities. Despite species distribution status and habitat loss degree, the fraction of unique species between two non-overlapped local communities would decrease for the lognormal model and increase for the logseries model after habitat loss. In other words, the changes in the fraction of unique species number between two non-overlapped local communities have positive and negative values as shown in Fig. 2G,H, respectively. For the lognormal model, the reduction in the fraction of unique species becomes larger when the area size of a local community in a pair increases (the size of the other community is fixed) or when the habitat destruction degree is higher (Fig. 2G). By contrast, for the logseries model, the magnitude of change in the fraction of unique species of two non-overlapped local communities would always increase when the area size of either local community increases (Fig. 2H). Moreover, when species distribution is random, there is almost no change on the fraction of unique species for the logseries model, no matter how the habitat loss degree is (Fig. 2H).
Effects of shifting SAAD. When species distribution pattern is allowed to vary before and after habitat loss, both consistent and new patterns are observed (Fig. 3). The consistent aspects include immigration credit patterns for endemic species of a single local community (Fig. 3C,D), which are similar to those results when SAAD is fixed (Fig. 2C,D). The contribution of species with different population sizes to the overall debt or credit magnitude of the whole community is similar (Figs 2E,F versus 3E,F).
Other than these, some new interesting patterns emerge. Firstly, immigration credit can occur in the local community when species distribution shifts from aggregation to randomness (Fig. 3A,B). Moreover, when habitat destruction degree is low, immigration credit always persists, no matter how the area size of the local community varies (Fig. 3A,B). Secondly, the fraction of unique species can decrease initially when the local community size is small but increase later when SAAD is shifted from randomness to aggregation (Fig. 3G,H). Thirdly, irrespective of the population threshold set and the SAD model used, immigration credit persists when spatial distribution is shifted from aggregation to randomness and habitat loss percent is 10% (Fig. 3E,F). Finally, immigration credits may occur for endemic species in both SAD models, which were not observed when SAAD was fixed (Figs 2C,D versus 3C,D).

Effects of extreme species distribution.
Most results for the extreme aggregation (k = 0.001) and randomness (k = 1000) scenarios ( Supplementary Fig. S1) are similar to those reported above when k = 1 and k = 100 are assumed (Fig. 2 versus Supplementary Fig. S1). Moreover, akin to the random distribution case in the logseries model (Fig. 2H), no (or nearly no) changes in the fraction of unique species have been found in the extreme aggregation case (k = 0.001), regardless of the SADs used ( Supplementary Fig. S1G,H).

Discussion
The present study develops a very flexible framework to predict the delayed responses of species diversity and community structure due to habitat destruction in local communities with various area sizes. The results provide primary theoretical evidence to demonstrate how extinction debt or immigration credit emerges; how the fraction of unique species changes; and how they vary along with changing area size of the local community, changing levels of habitat destruction, shifting species distribution, and shifting regional species abundance. In particular, this study is the first to explicitly demonstrate that habitat destruction can result in more endemic or rare species. Correspondingly, endemic and rare species do not necessarily contribute to extinction debt. Instead, they can contribute to immigration credit. Our study proves that there will be no change in the fraction of unique species when species distribution is highly aggregate.
The present results provide some interesting patterns. For example, for the single-parameter logseries SAD model, extinction debts can occur in local intact communities when habitat loss happens elsewhere (Fig. 3B). For the lognormal model, extinction debt is expected when species distribution is random (Figs 2A and 3A) 9 , but immigration credit can also occur as long as SAAD can change (Fig. 3A and B). These findings suggest that there are diverse delayed responses of ecological communities, depending on a variety of ecological factors (Figs 2 and 3 and Supplementary Figs S1 and S2). More importantly, these diverse area-dependent change in local diversity and community structure predicted by the general framework perfectly fits the possible empirical observations, such as climate change or habitat loss that may increase or decrease alpha diversity, the fraction of unique species, or both (Table 1).
Contradictory to the argument made in KH's paper, the present study theoretically demonstrates that logseries SAD can predict either extinction debt or immigration credit in the local community (Figs 2B and 3B, and Supplementary Fig. S1B). In particular, for the most restricted case when both SAD and SAAD are fixed (which was investigated in KH's study), extinction debt occurs ( Fig. 2B and Supplementary Fig. S1B). So how can KH show that logseries model predicts zero extinction debt in their study? A key reason is that their study assumed the alpha parameter to be constant when applying Fisher's logseries model (Eq. 12a). Finally, it is worth mentioning that, in another study 30 , Fisher's logseries model actually has been employed to predict delayed extinction of Amazonian trees when equipped with a neutral model that they developed in their study. As such, the present Scientific RepoRts | 7: 998 | DOI:10.1038/s41598-017-01070-2 Figure 3. Prediction of extinction debt or immigration credit for a single local community and changes in the fraction of unique species for a pair of local communities in the remaining intact region. All scenarios here assume SAD to be fixed, while SAAD can shift (exchangeable between random and aggregated distributions) before and after habitat loss. Species distribution is assumed to be random (k = 100) or aggregate (k = 1). (A,B) Extinction debt or immigration credit as a function of the area size of a single local community; (C,D) extinction debt or immigration credit for endemic species in a single local community; (E,F) extinction debt or immigration credit contributed by small-population species (the area of the local community is fixed to have a fraction 0.4 over the region); and (G,H) changes in the fraction of unique species as a function of the area size of one local community in a pair (the area size of the other community is fixed to be 10% of the whole region). Before habitat destruction, the region is assumed to have 30 species and 200 individuals in the model. In all the subplots, dashed lines represent the outcomes when 50% of the region is destructed, while solid lines represent the outcomes when only 10% of the region is destructed.
findings and the previous study 30 jointly support the fact that logseries SAD model shall predict extinction debts, depending on factors like model formulation method.
An interesting finding of this study is that endemic species (subplots C and D in Figs 2, 3 and Supplementary  Fig. S1) or species with small population sizes (subplots E and F in Figs 2, 3 and Supplementary Fig. S1) are predicted to increase in local communities in many cases. Thus, immigration credit is further strengthened, particularly for endemic species, when habitat destruction degree is high (subplots C and D in Figs 2, 3 and Supplementary Fig. S1). One possible reason is that habitat destruction drives many common species to go extinct and become rare or endemic species (because some of their populations in other areas have gone extinct),

Possible empirical patterns
Possible explanations using the proposed framework Local diversity versus regional diversity change 1.
Local diversity increases; Regional diversity increases Immigration credit occurs across all the spatial scales. This can happen when SAAD is shifted from aggregation to randomness and habitat destruction is low (Fig. 3A,B).

2.
Local diversity increases; Regional diversity decreases Immigration credit occurs at local communities with small area sizes, while extinction debt occurs for the whole remaining region. This may happen when SAD is fixed to be lognormal, while SAAD is allowed to shift from aggregation to randomness and when habitat destruction is high (Fig. 3A).

3.
Local diversity decreases; Regional diversity decreases Extinction debt occurs at all spatial scales. This may happen when (1) both SAD and SAAD are fixed before and after habitat loss ( Fig. 2A,B) or (2) when SAD is fixed, while SAAD is shifted from randomness to aggregation (Fig. 3A,B).
Changes in the fraction of unique species 1.

Fraction of unique species increases after habitat loss
This may occur (1) when the sum of area sizes of a pair of local communities is large enough and when SAAD is shifted from randomness to aggregation, while SAD is fixed (Fig. 3G,H) or (2) both SAD and SAAD are fixed in logseries model (Fig. 2H).

Fraction of unique species decreases after habitat loss
This may occur (1) when both SAD and SAAD are fixed in lognormal model ( Fig. 2G) or (2) when SAAD is shifted from aggregation to randomness, while SAD is fixed (Fig. 3G,H).

3.
No fraction of unique species change after habitat loss This may occur when species distribution is highly aggregate, regardless of SAD models used ( Supplementary Figs. S1G,H).
Alpha diversity versus change in the fraction of unique species 1.

Local alpha diversity increases; fraction of unique species increases
This may occur when species distribution pattern is extremely aggregated in lognormal model (Fig, S1A versus S1G).

2.
Local alpha diversity increases; fraction of unique species decreases This may occur when SAD is fixed, while SAAD is shifted from aggregation to randomness (Fig. 3A Table 1. Explaining diverse empirical patterns on the local-regional diversity change and changes in the fraction of unique species caused by habitat destruction (or climate change) using the general framework proposed in the present study.
resulting in a higher degree of endemism and rarity in local communities in the remaining habitat after destruction. This evidence is supported by high extinction debts in non-endemic species in the focused local community for many scenarios (Supplementary Fig. S2). This can happen when the whole community has a high debt ( Fig. 2A,B), while endemic species have an immigration credit (Fig. 2C,D). The subtraction of both will lead to a great loss or small gain in non-endemic species after habitat destruction ( Supplementary Fig. S2). Greater loss of non-endemic species is found when the sampling area becomes larger (Supplementary Fig. S2). This implied that common species may be more vulnerable to habitat destruction and may undergo delayed extinction. This can be reasonable, if these common species prefer to inhabit the destructed habitats, or the destructed habitats become keystone ecosystems to support all species.
Habitat destruction degree, area size of the local community, SAD, and SAAD models jointly play deterministic roles affecting the delayed changes in local diversity and community structure. In particular, shifting SAAD and extreme species distribution can generate new patterns ( Fig. 3 and Supplementary Fig. S1). For example, under extreme aggregation condition (k = 0.001), there is no net change in the fraction of unique species (Supplementary Fig. S1G,H). The key fact is that the fraction of unique species computed for a pair of local communities under maximal aggregation condition is one, regardless of regional SADs and habitat destruction degree used (see the proof of Theorem 2 in the Supplementary Information). Thus, these results imply that under extreme aggregation condition, each local community would be predicted to have completely different species composition (and all of them are endemic) with respect to any other communities.
It is worth mentioning that zero change on the fraction of unique species does not merely occur in the extreme aggregation condition. Actually, when species distribution is very random, the fraction of unique species in the single-parameter logseries SAD model is also unchanged, regardless of habitat loss degree (Fig. 2H). The reason is greatly due to the extreme right-skewed curve shape of the logseries SAD model (thus predicting rare species everywhere) and the asymptotic binomial distribution pattern of species when k is large (see the proof of Theorem 3 in the Supplementary Information). Finally, in the proof of Theorem 3, it can be seen that the change on the fraction of unique species is not related to the area size of the lost habitat after some algebra. This well explains why the zero net change occurs as long as species distribution is random.
In conclusion, the present study has theoretically discovered very diverse and interesting post-destruction responses of species diversity and community structure, which can match and explain a variety of empirical observations (Table 1). We have shown that there are area-and distribution-dependent transitions between extinction debts and immigration credits. Endemic species increase in the local community after habitat loss in many situations (Fig. 2C,D and Supplementary Fig. S1C,D). The fraction of unique species may increase, decrease, or remain unchanged depending on multiple factors, including the SAD model used, the area size of the local community, habitat destruction, and species distribution pattern (Figs 2G,H and 3G,H and Supplementary  Fig. S1G,H). There is no change in the fraction of unique species number in extremely aggregate distribution situation (k → 0), regardless of regional SAD models used ( Supplementary Fig. S1G,H).
For future implications, the succeeding steps are to further verify these diverse theoretical patterns predicted by the proposed framework at local-scale field experiments. Moreover, because there have been some empirical studies reporting the evidence of extinction debts 6, 31 , data collected from previous studies (if available) might be re-used to examine the possible alteration of community structure and change in the fraction of unique species after habitat loss. In addition, there are many other intriguing areas open for future research. For example, if observational data from the fields are available, it will be interesting to verify whether the fitted SAD models before and after habitat loss would have the same parametric form. It is also interesting to experimentally test the alteration in species distribution after habitat loss and its influence on the local community structure.