Plant-plant interactions influence post-fire recovery depending on fire history and nurse growth form

Backgrounds Plant-plant interactions are among the most important factors affecting the natural recovery of vegetation. While the impacts of nurse plants on species composition and biodiversity are well documented, the effects of different nurse’s growth forms on all biodiversity components including taxonomic, functional, and phylogenetic diversity have been less studied and compared, especially for their effects on different times after fire disturbance. This research was focused on comparing the effects of a perennial grass ( Elymus hispidens ), a perennial herb ( Phlomis can-cellata ), and a high shrub species ( Lonicera nummulariifolia ) on species composition and the biodiversity components, and how these impacts change across five sites with short-term (1 and 4 years sites), long-term (10 and 20 years sites) times since last fire and a control site where no fire was known in recorded history in semi-arid shrublands of Fereizi Chenaran located in Northeast of Iran. The changes of species composition and taxonomic, functional, and phylogenetic diversity were calculated with respect to the presence/absence of nurse’s growth forms, fire history, and their interactions. Results Nurse shrubs affected species composition and all biodiversity components, whereas all indices were reduced when considering Elymus grass as nurse plant. On the other hand, the herb Phlomis enhanced species composition and taxonomic diversity, while it had a negative effect on functional and phylogenetic diversity. Such specific effects of nurse types were mostly observed under long timescales (i.e., 10-and 20-year sites). Interestingly, the relative importance of nurse types and time since the last fire largely explained the variation of species composition and biodiversity components, with larger effects of nurse types on all biodiversity components. However, we found a significant contribution of fire explaining variation of species composition and phylogenetic diversity. Conclusions These results indicated nurse plants can affect the post-fire recovery of vegetation by providing specific mechanisms controlling beneficiary relatedness depending on their growth forms and time scales since the last fire. Therefore, these findings suggest perennial plants in the form of nurse species as a useful factor to develop techniques of active restoration in burned ecosystems.


Background
Wildfires are major natural disturbances impacting plant communities and ecosystems (Bowd et al. 2018).A large body of evidence suggests that fire has pronounced effects on plant biodiversity and the functioning of terrestrial ecosystems (Ghasempour et al. 2022;Omidipour et al. 2021;Bowman et al. 2016;Hurteau et al. 2014;Pausas and Ribeiro 2013).For instance, fire can significantly change taxonomic composition and promote functional/phylogenetic divergence or convergence (Gómez-González et al. 2011;Hernandez-Serrano´ et al. 2013) through diversification of plant lineages (Bytebier et al. 2011;He et al. 2012;2019).However, such important impacts of fire on vegetation are closely related to different time scales (Keeley et al. 2012;Gosper et al. 2013), presence of dominant perennial plants (Padilla and Pugnaire 2006), and ecosystem type (McLauchlan et al. 2020;He et al. 2019).The complex relationships between these important factors have posed challenges for fire ecologists in formulating clear theories and hypotheses regarding post-fire vegetation recovery (Mucina et al. 2019).Therefore, the evidence suggesting key factors and mechanisms that promote natural post-fire vegetation recovery has been widely embraced in fire ecology, fostering ecological knowledge and facilitating the development of active restoration practices.
Natural vegetation recovery is a lengthy process that occurs without support and has significant impacts on biodiversity, ecosystem functions, and services in degraded ecosystems (Sandell Festin et al. 2019;Bashirzadeh et al. 2022b;Wortley et al. 2013;Wilgan et al. 2011).In this context, plant biodiversity and composition have been identified as crucial indicators of the natural recovery of vegetation (Wortley et al. 2013;Badano and Cavieres 2006).Understanding key factors and mechanisms of this process through plant biodiversity and species composition can potentially inform the development of active restoration practices in burned ecosystems (Padilla and Pugnaire 2006).However, plant ecologists face major challenges when considering plant biodiversity, which encompasses various components such as taxonomic, functional, and phylogenetic diversity (Swenson et al. 2006;Cadotte and Davies 2016).Making connections between the results of these components and their interpretations can be difficult.Therefore, numerous studies in fire ecology have focused only on the floristic composition and taxonomic component (Bahalkeh et al. 2021), while less attention has been given to other components such as functional diversity (Abedi et al. 2022) and phylogenetic diversity (Bashirzadeh et al. 2023;Hipp et al. 2015;Keeley and Fotheringham 2003).It has been shown that natural vegetation recovery is observable not only in variations of floristic composition, but also in functional traits and evolutionary relatedness of plant species within communities (Keeley et al. 2012;Dantas et al. 2013).Thus, the overall outcomes of post-fire vegetation recovery are significantly dependent on biodiversity components.
Natural post-fire vegetation recovery is influenced by several major factors such as the presence of perennial plants (Ren et al. 2008;Niknam et al. 2018), spatial scale considered (Gomez-Aparicio 2009), and time scales since the last fire occurrence (He et al. 2019).A wide range of studies have recognized the positive effects of perennial plants, acting as nurse plants, on vegetation structure and soil seed bank in severe ecosystems such as alpine environments and global drylands (Bashirzadeh et al. 2022a;Erfanzadeh et al. 2022;Pashirzad et al. 2019a;Pashirzad et al. 2019b;Soliveres et al. 2015;Soliveres et al. 2014;Soliveres and Maestre 2014;He et al. 2012;Gomez-Gonz ´ alez ´ et al. 2011) and degraded ones (Gomez-Aparicio 2009; Padilla and Pugnaire 2006;Bashirzadeh et al. 2022b).The findings of such studies have provided valuable insights into natural vegetation recovery (Bond et al. 2005;Calitz et al. 2015;Bashirzadeh et al. 2022b), and have contributed to the development of restoration practices (Padilla and Pugnaire 2006;Ren et al. 2008).These findings and theories have primarily defined the impacts of nurse plants on the recruitment and survival of less tolerant plant species by creating suitable microenvironments under their canopies.
As a result, they contribute to increasing the taxonomic diversity of the community (e.g., Soliveres and Maestre 2014;Le Bagousse-Pinguet et al. 2019;Madrigal-Gonzalez et al. 2020).However, the specific effects of nurse species are not solely limited to adding new species and increasing the species pool.Differences in functional traits of plant species have been observed between plant species that grow under nurse species and those that can thrive in their absence (Rafiee et al. 2023;Soliveres et al. 2014).Furthermore, the maintenance of such traits across phylogeny has been observed in plant communities in relation to the presence of nurse species (Keeley et al. 2016;2012;Abedi et al. 2022, He et al. 2019).However, our understanding of the mechanisms promoted by nurse species in the recovery of natural vegetation in burned ecosystems, particularly in relation to fires, still has major gaps.This is because, in addition to the general effects of fire as a major disturbance factor on vegetation, there are also specific effects of fire that mediate facilitative plant interactions in burned ecosystems (Vega-Alvarez et al. 2019;Keeley et al. 2012).These fire effects may activate adaptive mechanisms that promote natural recovery of vegetation by increasing the distribution of traits related to post-fire recruitment and growth strategies (Alam et al. 2020;Lavorel et al. 2007;Pausas et al. 2003).
Interestingly, the impacts of nurse species on vegetation recovery in disturbed ecosystems can vary depending on their growth forms, study scale considered, and time intervals after disturbance (Gomez-Aparicio 2009;Erfanzadeh et al. 2014;He et al. 2012).In the case of alpine environments, there is strong evidence that shrubs have the most important effects on taxonomic and phylogenetic diversity (Butterfield et al. 2013;Cavieres et al. 2013;Camac et al. 2013;Svriz et al. 2013).On the other hand, in arid ecosystems, nurse grasses can strongly affect functional diversity by providing suitable microhabitats for beneficiary species with distant relatedness (Gastauer et al. 2018;Navarro-Cano et al. 2018;Bashirzadeh et al. 2022b).However, such effects of nurse growth forms are highly scale-dependent and should be considered across the scale of vegetation patches to the composition of the entire community, where the most effective theories and empirical perspectives in plant ecology are structured (Soliveres et al. 2015).Additionally, besides considering nurse growth forms and scale, the time interval since the last likely fire influences vegetation recovery in burned ecosystems (Paniw et al. 2017).In this context, species diversity is expected to recover significantly over time following a fire, as the time since the last fire increases (Abedi et al. 2022;Bashirzadeh et al. 2023).Herath et al. (2009) found that nearly all plant species in fire-prone ecosystems that were present in the pre-fire communities, were evident for 12 months after the fire.This phenomenon is consistent with the "autosuccession" model (Hanes 1971) and the "initial floristic composition" model (Enright et al. 2007) of community succession (He et al. 2019).However, it is not clear which forms of nurse growth and at what intervals contribute to this phenomenon in burned ecosystems.Therefore, further research studies are needed to clarify the relative impacts of nurse types on the natural recovery of all biodiversity components considering times scales and the proper scale of burned ecosystems.
Here, we analyzed the effects of three dominant perennial plants with different growth forms (Elymus hispidus subsp.hispidus, Phlomis cancellata, and Lonicera nummularifolia with perennial grass, perennial forb, and sprawling high shrub growth forms, respectively) on species composition and taxonomic, functional and phylogenetic diversity, over different time periods since the last fire (1, 4, 10, and 20 years since the last fire), and compared these to a control site, where no fire was known in recorded history.This research study was conducted in semi-arid scrublands of the Fereizi region in northeastern Iran, which can serve as an empirical example of post-fire vegetation recovery in scrublands.In addition, previous studies only addressed either the effects of fire in terms of fire intensity (Abedi et al. 2022;Heydari et al. 2017) or the presence of nurse plants on taxonomic diversity and species composition at point basis scales (i.e., vegetation patch scales) (Abedi et al. 2022;Bahalkeh et al. 2021;Paniw et al. 2017).Therefore, some important issues such as how nurse species affect all biodiversity components (i.e., taxonomic, functional, and phylogenetic diversity), which mechanisms are used to facilitate their natural recovery, and relevant contributions of the nurse species across time scales in burned ecosystems have remained unclear.In this regard, it is of great importance to consider some key functional characteristics of nurse species, including growth forms, when assessing their impacts in order to fill the gaps in understanding of the natural vegetation recovery after fire, and to provide important insights for the development of active restoration practices.Therefore, we obtained information on nurse plant growth forms, functional trait values, and phylogenetic relationships among plant species, and addressed the following questions: (i) Does each biodiversity component exhibit significant variation in relation to nurse plant types across different times since the last fire?(ii) Are the specific effects of each nurse plant on biodiversity correlated to times since the last fire, the specific biodiversity facet considered, and functional groups of plant species?(iii) To what extent do nurse plant types, time since the last fire, and their interaction explain the variance in different biodiversity components and species composition?

Study area
The study area consists of semi-arid shrublands of the Fereizi region located in northeastern Iran (between 36° 35′ and 36° 48′ N, and 58° 24′ and 58° 52′ E) (Bashirzadeh et al. 2023) (Fig. 1a).Elevation ranges from 1980 to 2323 m, generally increasing from south to north.Mean annual precipitation ranges from 285 to 312 mm with a cold semi-arid climate based on Emberger's method (Additional file 1: Appendix S1).The rain season is from Autumn (late October) to early Spring, with the main rain season occurring from February through April, with an average temperature of 7.5° C (Memariani et al. 2009).
Altitude,, and temperature values were obtained from the WorldClim dataset in 30″ resolution (Hijmans et al. 2005).Plant communities are typically dominated by herbaceous plants and shrubs (Bashirzadeh et al. 2023) (Bashirzadeh et al. 2023;Memariani et al. 2009).Our analyses focused on the three most dominant plant species as nurse plants, L. nummulariifolia, P. cancellata, and E. hispidens, which occurred at all sample sites.L. nummulariifolia is a sprawling high shrub species (plant height: 3-5 m) with a wide crown (Barzegaran et al. 2022).Especially under moderate disturbance conditions, it became clear that such broad-crowned shrubs have a positive effect on understory plant species (Jankju and Ejtehadi 2016;Soliveres and Maestre 2014;Soliveres et al. 2012).Sprawling shrubs generally create suitable microenvironments for understory species by increasing soil moisture and reducing radiation and air temperature compared to open areas.P. cancellata was considered as a perennial herbaceous plant based on Raunkiaer's system of classification, previous studies conducted in northeastern Iran (Erfanian et al. 2021;Atashgahi et al. 2018;Memariani et al. 2016), and also according to our assessments of this plant species in the study region although it was highlighted as a sub-scrub plant in a reference (https:// powo.scien ce.kew.org/ taxon/ urn: lsid: ipni.org: names: 453832-1).We included this plant species as a nurse plant in our further analyses because it promotes the growth, recruitment, and survival of other species (Soliveres et al. 2014;Jankju 2013;Jankju et al. 2008).E. hispidens is a perennial grass species that can grow up to 1.5 m tall, and its limited crown provides suitable habitat for some plant species with specific niches (Gomez-Aparicio 2009;2004).

Site selection
We selected five sites based on different times since the last fire in the study area.They included a control site where no fire was known in recorded history with natural vegetation, a 1-year site in which the last fire occurred in the last year (i.e., in 2021), a 4-year site, a 10-year site and a 20-year site, in which fires occurred 4, 10, and 20 years ago, respectively (Fig. 1a).These sites encompass similar ecological and topographical conditions, ensuring minimal variations in terms of climatic and topographical variables.These variables include annual precipitation and altitude, which range from 285 mm and 1980 m in the 4-year site to 312 mm and 2323 m in the control site (Additional file 1: Appendix S1).Furthermore, all the sites are located on the same geological substrate composed of limestone (Memariani et al. 2009).From all these sites, we extracted comparable data on altitude and precipitation from the WorldClim dataset (http:// www.world clim.org; Hijmans et al. 2005).To examine fire occurrence, we extracted the time since the last fire for each site from aerial photos and from information on fire time recorded by the staff of the environment department in the Fereizi region.
We estimated fire intensity for each site based on classifications from Heydari et al. (2017).Within these classifications, fire intensity is calculated using some specific indicators.Such specific indicators are (i) apparent fire; (ii) herbaceous cover; (iii) the percentage of shrubs, herbs, and trees killed; and (iv) stem charring of tree height (Heydari et al. 2017).In the 1-year plot, we found ground fires and burning lower tree branches.Herbaceous plants and some shrubs were killed.The stems of some plants were burned up to 70% of the whole total tree height.Additionally, between 25 and 30% of shrubs and herbs had been dried because of the fire (Fig. 1b).Such observations were also made at other burned sites by the staff working in the environment department at the Fereizi region.Finally, in the control site, no evidence of vegetation destroyed by fire was found, and no charring on tree stems was found as well.Based on such findings, all burned sites indicated an intermediate or moderate level of fire (Fig. 1b).

Vegetation sampling
To assess nurse plant effects on vegetation recovery across different times since the last fire, two transects in each site with 100 m length and approximately 500 m apart were established (Fig. 1a) on north-facing slopes, where similar ecological and topographical conditions were mostly found across transects.Then, on each transect, 20 individuals of each dominant plant species (i.e., L. nummulariifolia, E. hispidens, and P. cancellata) were randomly selected as potential nurses, and 20 1 m × 1 m quadrats were surveyed under the canopies of these dominant plant species, and the same number of quadrats in open areas far from such potential nurses.Then, the vegetation within the quadrats was surveyed as a percent coverage of present taxa from early May 2021 to late June 2021.To identify different microclimatic enhancement mechanisms of nurse species, it is essential to survey the vegetation beneath nurse plants and their open spaces using the same number and size of quadrats.This is because the measured plant biodiversity varies greatly depending on the sampling effort (Erfanzadeh et al. 2022;Pashirzad et al. 2019a;2019b;Soliveres et al. 2012).Therefore, the number of transects and quadrats laid out was determined based on two main criteria: (i) the extent of study sites and (ii) species-area curves (Bashirzadeh et al. 2022b).We identified a total of 75 plant taxa in our study sites.To avoid the effects of rare species, we excluded 20 (out of the 75) with a frequency below five in the entire dataset (Bahalkeh et al. 2021) (see details in Additional file 2: Appendix S2 in Supporting Information).

Functional trait information
We assessed eight functional traits for plant species in our sites, including leaf dry matter content (LDMC) (g), specific leaf area (SLA) (m 2 kg − 1), leaf nitrogen (mg g − 1) and carbon content (mg g − 1), plant height (m), seed mass (mg), growth form, and life span.These functional traits were introduced as key functional traits for plant species in burned ecosystems, and are dependent on biotic interactions and disturbance conditions (Rafiee et al. 2022;Lienin and Kleyer 2012;Lavorel et al. 2007;Meng et al. 2015;Abedi et al. 2022).For example, SLA and LDMC reflect plant growth rates in fire-prone ecosystems.Fire generally increases SLA and decreases LDMC because fire increases the abundance of fastgrowing and ruderal species compared with perennial species (Davies et al. 2009;Rhodes et al. 2010).In addition, other traits related to plant regeneration including seed mass, plant height, and plant growth also represent plant response to fire (Cruz et al. 2022).Life span represents plant strategies and abilities including competitive or facilitative abilities, and leaf nitrogen/carbon content and life span can indicate post-fire resprouting in shrubs (Anacker et al. 2011;Sakschewski et al. 2015), and differentiate annual from perennials plant species in herbaceous communities (Keeley et al. 2012).
We obtained information on these plant functional traits by measuring some traits including plant height and life span in field and others from publicly available trait datasets (BIEN package in R (Maitner et al. 2018), TRY (Kattge et al. 2020), LEDA (Kleyer et al. 2008) and TR8.To measure plant height in the field, we selected 10 individuals of each plant species and measured plant height as the mean of these individuals, focusing on the vegetative part.This approach was chosen to minimize bias introduced by the reproductive parts, which can vary at different growth stages, particularly in grasses (Cavieres et al. 2013).For other functional traits, we averaged observations of plant species individuals to create a species mean trait value, when multiple measurements per species were available in these databases.For some plant species, we used genus-level means when either species-level data were not available or plants were identified only to genus level (Lamanna et al. 2014) (see Supporting Information Additional file 2: Appendix S2).

Phylogenetic information
We obtained a phylogeny of 55 plant species present in all study sites based on the most up-to-date megaphylogeny for seed plants.We standardized species names in our dataset according to The Plant List using the R package "Taxonstand" (Cayuela et al. 2012).Then, we used the R function V.PhyloMaker (Jin and Qian 2019) to link the species names in our dataset with those in the megaphylogeny, and the scenario 3 approach (Qian and Jin 2016), to add species to the phylogeny.This phylogenetic tree was then used as reference lists from which phylogenetic diversity could be calculated for our communities in the dataset.

Measures of taxonomic, functional, and phylogenetic diversity
To measure the effects of nurse's growth forms on biodiversity components in different times since the last fire, we created some communities in relation to the presence of nurse species such as L. nummulariifolia, P. cancellata, and E. hispidens by randomly substituting 50% of the quadrats sampled in open microsites per 50 m transect, with the same number of randomly selected quadrats at the same transect beneath Lonicera, Phlomis, and Elymus areas, respectively.While the latter does not consider the relative proportion of nurse/open areas at each site, using this standard proportion allowed us to create different communities on transects (two communities on each transect) and compare our results across sites with different times since the last fire (Bashirzadeh et al. 2022a).In addition, this procedure allows us to understand the effects of nurse plants in communities, which is a field of plant ecology where robust theories and findings are often made.We achieve this by comparing communities that include nurse quadrats to exclusively open communities (i.e., communities without nurse effects) (Bashirzadeh et al. 2022a).In this context, we considered communities with the presence of L. nummulariifolia, P. cancellata, and E. hispidens as Lonicera, Phlomis, and Elymus communities and their absences as Open communities for further analyses (i.e., four communities per 50 m transect).Then, we measured biodiversity components for these communities in each site.To measure taxonomic diversity, we used the first two Hill numbers to estimate species richness (q = 0) and species diversity as the exponential of Shannon's entropy (q = 1; referring to Shannon diversity) (Chao et al. 2014) for Lonicera, Phlomis, Elymus, and Open communities.The calculation was based on the percent coverage of plant taxa, which is less affected by differences in total sampling effort than other methods (Chao and Jost 2012).
To measure functional and phylogenetic diversity, we assessed abundance-weighted mean pairwise distance (MPD) (Tucker et al. 2017;Webb et al. 2002) as the most robust measure of the phylogenetic and functional relatedness for Lonicera, Phlomis, Elymus, and Open communities.This index in the abundance-weighted case is equivalent to Rao's Q and Hill numbers (Webb et al. 2002).We calculated the standardized effect sizes of abundance-weighted MPD for functional (FSES.mpd)and phylogenetic (PSES.mpd)diversity to produce a phylogenetic and functional index of diversity that is independent of species richness.We used the independent-swap algorithm to draw a null distribution based on 999 randomizations and created standard effect sizes of mpd by comparing the observed community diversity to the null distribution of randomly assembled communities.Positive (SES values > 1.96) and negative (SES values < − 1.96) values of SES.mpd indicate significant functional/phylogenetic over-dispersion and clustering patterns in the community, respectively.We calculated SES.mpd with the "ses.mpd"functions in R package picante (Kembel et al. 2010).

Measures of relative interaction index (RII) in communities considering different nurse types
We calculated the Relative Interaction Index (RII) for species cover to quantify the effect of each nurse species on other species in different combinations of the three treatments (time since the last fire, forms of nurse growth, and their interactions), following Armas et al. (2004): where X Lonicera , X Phlomis , X Elymus , and X open represent the total plant cover in the four different patches.This index is symmetric about zero (no significant interaction) and varies between 0 and 1 for facilitation and between 0 and − 1 for competition.It was calculated using the summed total plant cover of all species present in communities considering the presence/absence of nurse types.

Measures of species composition and functional groups
We performed a Correspondence Analysis (CA) on Lonicera, Phlomis, Elymus, and Open communities to measure species composition for these plant communities.Correspondence Analysis was performed using the "ca" function in R package ca (Greenacre et al. 2016).Lonicera, Phlomis, Elymus, and open communities were plotted at CA as the distribution of their quadrats (i.e., nurse and open quadrats) with respect to time since the last fire.Then, the amount of variation among the communities on CA axis 1 and CA axis 2 (i.e., percentage of Eigen values for CA axis 1 and CA axis 2) was plotted with respect to times since the last fire.
In addition to the CA ordination of communities in terms of the dispersion of their quadrats, variations in the species composition of communities in CA spaces were determined based on the dispersion of their plant species (Bahalkeh et al. 2021).We then performed a cluster analysis of the functional traits matrix, where the plant species present in the communities were presented in rows and the functional traits in columns to identify groups of species with similar functional identities.Cluster analysis was performed using the "hclust" function based on the Ward.D method in the R package vegan (Oksanen et al. 2016).Before performing the cluster analysis, we measured the functional dissimilarity between plant species using the Gower distance index.Finally, cluster groups were plotted on CA species space to assess the presence of plant species and cluster groups in the communities in relation to time since the last fire.

Statistical analyses
We analyzed variation in RIIs of nurse types and the taxonomic, functional, and phylogenetic diversity relative to the presence/absence of nurse types, different times since the last fire, and their interactions.We analyzed the response of RIIs and biodiversity facets (i.e., taxonomic, functional, and phylogenetic diversity) (using separate analyses for each) relative to the presence/absence of Lonicera nummulariifolia, Phlomis cancellata, and Elymus hispidus in communities by developing linear mixed effect models.The models included the presence/absence of nurse's growth forms, different times since the last fire and their interactions as fixed effects, and study area as a random factor to account for the lack of independency between communities including different nurse's growth forms and open communities.The models were fitted using the lmer () function in the R package lmerTEST (Kuznetsova et al. 2017), and their marginal R 2 values were obtained using the function "r.squaredGLMM" provided in the package "MuMIn" (Barton 2013).Marginal R 2 is the proportion of variance explained solely by our fixed factors (presence of nurse types and times since the last fire).Finally, we calculated p-value for individual treatments (i.e., fixed effects) from F statistics of type III sum of squares with the Satterthwaite approximation to estimate the denominator degrees of freedom using function "summary" provided in the package "lmerTEST" (Kuznetsova et al. 2017).
To examine the relative importance of times since the last fire and the presence of nurse's growth forms in explaining the variation in each diversity index, we performed variation partitioning based on partial linear regression using the "varpart" function (Oksanen et al. 2016).The total percentage of variation explained was divided into a unique and shared contribution for two sets of predictors: (i) times since last fire (i.e., fire factor), (ii) presence of nurse's growth forms (i.e.nurse factor), and (iii) their interactions (i.e., shared environment between these factors).
To test whether the presence/absence of nurse's growth forms, different times since the last fire and their interactions significantly influence species composition, a linear mixed effect model was performed on CA quadrat scores for Lonicera, Phlomis, and Elymus communities (i.e., 50% nurse + 50% open quadrats) and Open communities (i.e., only open quadrats), separately for each CA axis, with presence/absence of nurse's growth forms, different times since last fire and their interactions treated as fixed effects, and the communities nested within the sites as the random effects.Finally, to assess the effect of these treatments on the cover of plant species in each cluster group, a linear mixed effect model, similar to previous ones, was performed on CA species scores.All linear mixed-effects models were fitted with the function lmer in the R package lmerTest (Kuznetsova et al. 2017), and P-values were calculated from F-statistics of type III sum of squares with the Satterthwaite approximation to estimate the denominator degrees of freedom.Dependent variables were checked for normality and log-transformed before performing linear mixed effect models (Legendre and Legendre 2012).

Effects of nurse's growth forms on relative interaction index values across times since the last fire
Time scale since the last fire, the presence/absence of nurse's growth forms and their interactions had highly significant effects on RIIs values for the whole understory community which refers to the presence of all beneficiary species in a given area (Fig. 2).We found negative RII values for the entire understory species in Elymus communities across all times since the last fire, including both short-term and long-term times.However, RII was much more negative in earlier time periods (i.e., 1-year and 4-year sites) compared to long-term times (i.e., 10and 20-year sites) in the Elymus communities (Fig. 2; Elymus).In contrast, positive RII values were found for understory communities in Phlomis and Lonicera communities across all time scales (Fig. 2; Phlomis and Lonicera).Notably, more positive RII values were obtained for understory communities in Lonicera compared to Phlomis communities across long-term times (10-year and 20-year sites) (Fig. 2).

Effects of nurse's growth forms on species composition across times since the last fire
All Lonicera communities located in different times since the last fire were scattered on the left-hand side of CA axis 1, and showed significant correlations with CA axis 2 (Fig. 3).In contrast, Elymus and Phlomis communities were scattered across CA axis 1, with more positive CA scores than negative scores (only observed in Elymus and Phlomis across 20-year site).Additionally, they exhibited stronger correlations with CA axis 1 than CA axis 2 (Fig. 3).Overall, plant communities, with respect to the presence of nurse's growth forms, occupied a wider range of spaces with diverse CA scores depending on times since the last fire.In this regard, communities at longer times (i.e., 10 and 20 years) exhibited larger CA spaces with diverse scores compared to communities at short-term times (i.e., 1 and 4 years) and the control site, whereas all nurse communities (except Elymus communities) were more scattered along CA axis 2 resulting in narrower CA spaces at short-term communities and the control site (Fig. 3).We found significant variation in species composition in relation to presence of nurse plants with different growth forms and time scales since last fire (Fig. 3).In terms of species composition on CA axis 1, the presence of nurse's growth forms, times since last fire and their interaction significantly explained the observed variation.On the other hand, the time scale and its interaction with nurse plants were the most important factors in explaining variations in species composition on CA axis 2 (Fig. 3; Axis 1 and Axis 2).
Five groups were identified in the cluster analysis conducted on the functional trait matrix (see details in Additional file 3: Appendix S3).Group A included only annual herbaceous species with plant height less than 0.50 m (e.g., Callipeltis cucularis, Galium aparine, Asperula glomerata, and Ixilirion tataricum), located on the upper negative side of the CA axis 1, where Elymus and Phlomis communities with long-term times (i.e., 10-and 20-year communities) were scattered on the ordination space (Fig. 3, Additional file 2: Appendix S2).Group B included annual herbaceous and grass species ranging in height from 0.014 to 0.53 m, including Bromus tectorum, Scandix stellata, Filago arvensis, Stipa arabica, Achillea wilhelmsii, Taeniathrum caput-medisae, and Vicia ervilia, most of which were located on the relatively upper left side of CA -axis 1, where mainly Lonicera was scattered at long-times and Elymus and Phlomis at short-term times, as well as the control area in the ordination space (Fig. 3, Additional file 2: Appendix S2).Group C included perennial shrub and cushion species with plant height greater than 0.50 m (e.g., Lonicera nummularifolia, Rosa beggeriana, Acanthophyllum glandulosum, Ziziphora clinopodioides m, and Lactuca orientalis), scattered mainly on upper positive side of CA axis 1, where Lonicera and Phlomis communities occupied the ordination space in long-term times (Fig. 3, Group C).Groups D and E included perennial herbaceous species such as Phlomis cancellata, Eryngium bungei, Hymenocrater calycinus, Noaea mucrunata, and Thalictrum sultanabadense, with perennial herbs greater than 0.5 m tall in group D and those with less than 0.5 m were located in group E (Fig. 3, Additional file 2: Appendix S2).In addition, these plant species were located on the lower positive side of CA -axis 1, where communities were found in short-term periods (i.e., 1, 4-year communities) than in long-term times (Fig. 3).The linear mixed effect models indicated important effects of nurse plants and their interactions with time since the last fire on variation in functional groups on CA axis 1 (Fig. 3).In addition, time scale and its interaction with nurse plants had highly significant effects on CA axis 2 (Fig. 3).

Effects of nurse's growth forms on biodiversity facets across times since the last fire
Taxonomic indices responded significantly to the presence of nurse plants and time scales since the last fire (Fig. 4; q0 = mR 2 adj = 0.172, q1 = mR 2 adj = 0.276).In this context, an increase in taxonomic diversity based on both diversity indices (i.e., q0 and q1) was observed in Lonicera and Phlomis communities at longer times (10 years and 20 years, respectively).In contrast, Elymus communities indicated a decrease in taxonomic diversity in compared to open communities under long-term times since the last fire.The control site illustrated higher taxonomic diversity in open communities (Fig. 4; q0 and q1 panels).Interestingly, nurse plants and their interactions with time scale since the last fire were recognized as the most important predictors of variations in q0 (species richness) and q1 indices between sites (Fig. 4; q0 and q1 panels).
Significant and consistent effects of nurse's growth forms and their interaction were observed for functional and phylogenetic diversity.In addition, time since the last fire indicated the significant effects on these biodiversity indices especially on variations of phylogenetic diversity (Fig. 4; FSES.mpd = mR 2 adj = 0.433; PSES.mpd = mR 2 adj = 0.389; please see mR 2 adj values and significance levels for the treatments on biodiversity indices in Fig. 4).In this regard, higher functional and phylogenetic diversity (i.e., functional and phylogenetic overdispersion) were observed in Lonicera communities in all times since last fire except in the control site (Fig. 4; Lonicera).In contrast, Phlomis and Elymus communities indicated a significant decrease in functional and phylogenetic diversity (i.e., functional/phylogenetic clustering) under all times especially in short-term times (i.e., 1 and 4-year sites) for Phlomis communities (Fig. 3; Phlomis) and under long-term times for Elymus communities (Fig. 3; Elymus and Control).However, stronger and consistent decreases in functional and phylogenetic diversity were observed for Elymus communities than Phlomis communities, especially under long times (Fig. 4).

Relative importance of fire and nurse's growth forms on biodiversity facets
The amount of variance explained for all diversity facets increased when considering nurse plants (orange fraction in Fig. 5).This was particularly true for all biodiversity indices, especially for functional and phylogenetic diversity (please see orange fraction (nurse factor) in Fig. 5).In addition, time since last fire was a particularly strong predictor of variation in all biodiversity facets important for phylogenetic diversity (Fig. 5; yellow fraction (Fire factor)).Interaction between the presence of nurse plants and time since the last fire (see shared area between orange and yellow fractions) indicated contributions to explain variations in biodiversity components (Fig. 5; see shared area between orange and yellow fractions for all biodiversity facets).Although such contributions were not the strongest for biodiversity components, this suggests that the effects of nurse's growth forms on biodiversity facets are relatively sensitive to time since the last fire.

Discussion
The influence of nurse plants on biodiversity and species composition, particularly in severe ecosystems such as drylands and alpine and mountainous semi-arid ecosystems, has received some attention in the last two decades (Bashirzadeh et al. 2022a;Rafiee et al. 2022;Pashirzad et al. 2019a;2019b;Butterfield and Briggs 2011;Valiente-Banuet and Verdu 2007;Cavieres et al. 2005).However, we still do not have enough information regarding nurse plant effects on post-fire recovery of vegetation in relation to important factors such as growth form, time since fire, and at an applicable scale for burned ecosystems and plant successions (He et al. 2019;Keeley et al. 2012).Here, we analyzed the effects of select dominant perennial plants representing grass, herb and shrub growth forms on species composition and taxonomic, functional and phylogenetic diversity across different time-scales in a semi-arid shrublands of the Fereizi region located in Northeast of Iran.
Our results show the specific effects of perennial plants on the natural recovery of vegetation depending on their growth forms and time intervals.For example, Lonicera nummulariifolia shrubs significantly enhanced species composition and all components of biodiversity especially under long times since the last fire, whereas Elymus hispidus grasses showed a negative effect on all biodiversity components through recruiting the only grasses and herbaceous plants.Phlomis communities indicated an increase in taxonomic diversity, while they had a negative effect on functional and phylogenetic diversity through including herbaceous and shrub species.Interestingly, even though the specific effects of nurse's growth forms on biodiversity components depending on their morphological characteristics, all the effects of nurse's growth forms peaked under long times since the last fire (i.e., 10-and 20-year sites).In addition, fire solely and in the form of interaction with nurse types was recognized as one of the important factors in explaining the variation of species composition and biodiversity facets in our study area.Based on these results, these factors have served as strong explanations for the major importance of nurse plants with different growth forms as drivers of vegetation restoration in burned ecosystems (Paniw et al. 2017;He et al. 2019).In addition, our results strongly indicate that fire, as an agent of disturbance and significant

Nurse plants influence species composition and biodiversity components depending on their growth forms in burned communities
The effects of different nurse plants on species composition and some specific biodiversity components were associated with growth form, suggesting that the morphological features of nurse plants, as well as some processes controlling beneficiary relatedness beneath canopies of nurse plants, likely have impacts (Gomez-Aparicio 2009).This finding suggests the presence of different ecological conditions within these communities, which can be attributed to the varying effects of nurse plants in creating favorable microclimates based on their specific characteristics (He et al. 2012).There are hypotheses suggesting facilitative impacts of nurse shrubs driven by their specific morphological features, such as a low root-to-shoot ratio (i.e., a less competitive strategy) (Piston et al. 2016;Gomez-Aparicio 2009), to support a high abundance of beneficiaries with different functional identities.Following this hypothesis, our results indicated an increase in taxonomic diversity consistent with the presence of post-fire recruiting plant species, which are usually short-lived herbs, grasses, and subshrubs that were mostly observed in Lonicera communities.Such observations could introduce drought-resilient shrubs as ecosystem engineers, promoting biodiversity by facilitating the survival and growth of other long-lived shrubs, herbs, and grasses (reviewed in Gómez-Aparicio et al. 2004;Calvo et al. 1999).On the other hand, these positive effects of shrubs on the recruitment of other plant species can provide new insights into developing plant restoration practices (Gómez-Aparicio et al. 2004;Castro et al. 2006;Siles et al. 2008;Padilla and Pugnaire 2006) by understanding the mechanisms promoted by nurse shrubs in the recovery of plant biodiversity.In this context, we found that nurse shrubs allow for a specific trade-off; competition among closely related taxa and facilitation among distantly related taxa (Soliveres et al. 2012).The significant increase in functional and phylogenetic diversity among beneficiary species in communities including Lonicera nummulariifolia shrubs represents the presence of plant species with different functional identities (i.e., those with complete differentiation in realized niches) (Bannister et al. 2019;Pashirzad et al. 2018;Bahalkeh et al. 2017).
In contrast, the nurse grass (Elymus hispidus) and herb (Phlomis cancellata) strongly influenced species composition, leading to a decrease in all biodiversity components, especially over long periods of time.Additionally, the nurse herb had a less extensive impact on functional and phylogenetic diversity, particularly at short-term sites.These results highlight the differences in mechanisms promoted by nurse grass compared to nurse herb, depending on their morphological characteristics.Grasses generally have a larger negative effect than herbs, which supports previous studies showing grasses to be better competitors with higher establishment and growth rates than herbs (Goldberg et al. 1999;Pywell et al. 2003).Grasses have fibrous roots and large root-to-shoot ratios, which allow them to compete efficiently for numerous resources.Due to this strong competitive ability, nurse grasses tend to restrict the presence of beneficiary species, ultimately leading to the dominance of stress-tolerant species.These stress-tolerant species are capable of tolerating the dual environmental filtering imposed by nurse grasses, such as annual and perennial forbs and grasses (Harrington and Johns 1990;Puerta-Piñero et al. 2007).These specific mechanisms alter the overall species composition of plant communities by introducing plant species that would not otherwise be present (Butterfield et al. 2013;He et al. 2012;Soliveres et al. 2015).Thus, our findings have revealed innovative features of nurse grass in the recovery of vegetation.These benefits of using grasses as nurse plants have not been commonly identified in restoration practices, except in wetlands (Egerova et al. 2003;McKee et al. 2007) and semi-arid steppes (Maestre et al. 2001).

Nurse plants influence species composition and biodiversity facets depending on time since last fire
We found that nurse species with different growth forms had a significant impact on species composition and biodiversity, especially at intermediate and long-time intervals since the last fire (i.e., at 10 and 20-year sites).These results provide strong empirical support for the intermediate disturbance hypothesis (IDH) (Huston 2014), which suggests the positive effects of nurse plants on biodiversity, resulting in an increase in biodiversity, consistent with an increase in time since the last fire.Furthermore, the positive effects of nurse plants on biodiversity reach their peak under moderate levels (Dee and Menges 2014;Beckage and Stout 2000).It is worth noting that nurse plants may only effectively restore vegetation in burned ecosystems under moderate and long-time intervals since the last fire (Tepley et al. 2017).Therefore, the occurrence of fire and its impacts on vegetation play a crucial role in determining the post-fire recovery of vegetation by influencing the mechanisms promoted by nurse plants (Castro et al. 2006;Gomez-Aparicio 2009;2004).In this regard, our results demonstrate the positive effects of fire on increasing the competitive and facilitative abilities of nurse grasses and shrubs, respectively.This finding contradicts previous research, which mostly argues that fire decreases competition (Bahalkeh et al. 2021).Therefore, we conclude that nurse types and time since the last fire are important drivers of the variance in species composition and biodiversity facets.
Furthermore, it was expected that fire would explain a larger amount of variance in phylogenetic diversity compared to other biodiversity facets, despite the strong influence of nurse plants on phylogenetic diversity (e.g., He et al. 2019;Paniw et al. 2017;Khapugin et al. 2016).Previous studies have highlighted the presence of strong phylogenetic conservatism among plant species in postfire recovery, indicating that these lineages have evolved in response to fire disturbances (He et al. 2019).Therefore, fire acts as a significant driver of diversification by promoting mutation, selection, divergence, and speciation (Stevens-Rumann et al. 2018).This leaves a strong fingerprint on current phylogenetic diversity and may explain the stronger effect of fire on this biodiversity component compared to other components (see variance partitioning analyses in Fig. 3; PSES.mpd).Consequently, analyzing biodiversity responses to nurse plants collectively and considering the time since the last fire offers strong potential for linking fire to the evolutionary history of post-fire vegetation and the mechanisms of nurse plants (Keeley et al. 2012;Verdu and Pausas 2007).

Management implications
The findings of this research study, which primarily focus on the significance of nurse plant growth forms in the natural recovery of vegetation, are consistent with the recent resurgence of facilitation as the dominant interaction in intermediate disturbance levels (Callaway 2007;Brooker et al. 2008).These findings could provide strong support for the use of plant-based techniques in the restoration practices of burned ecosystems.Implementing such procedures would be considered cost-effective and environmentally friendly, aligning with recent conceptual models that hypothesize plant facilitation as an exploitable interaction in active restoration (Navarro-Cano et al. 2018;Cuevas et al. 2013;Ren et al. 2008).In this regard, nurse plants were identified as a key factor in managing the conservation of vegetation and preventing the loss of plant biodiversity, particularly in the short time since the last fire.Therefore, we recommend the planting of seedlings or seeds of dominant perennial plants such as shrubs and herbaceous species as an effective measure for restoring burned ecosystems.This recommendation is based on our comprehensive evaluation of the effects of nurse plants on the natural recovery of vegetation, taking into account the different growth forms of nurse species, the appropriate scale of vegetation (i.e., community scale, where the most influential theories were developed), and the various time scales since the last fire.

Conclusion
In this study, we observed that all growth forms of nurse plants influenced species composition and various biodiversity components, regardless of the time since the last fire.Lonicera communities exhibited enhanced biodiversity facets alongside a significant change in species composition.Conversely, we found a decrease in taxonomic, functional, and phylogenetic diversity in Elymus communities.Additionally, we found that nurse herb species had a positive effect on species composition and taxonomic diversity in Phlomis communities, while it had a negative effect on functional and phylogenetic diversity.Interestingly, these specific effects of nurse types were predominantly observed in intermediate and long time since the last fire sites (i.e., 10 and 20 years).Based on these findings, we conclude that nurse plants have stronger and more consistent effects on species composition and overall biodiversity than the time since the last fire.However, fire made a significant contribution to phylogenetic diversity.These findings demonstrate the substantial effects of fire on the evolutionary history of vegetation during postfire recovery in burned ecosystems.Therefore, this study suggests that nurse plants can restore species composition and plant biodiversity, particularly under long times since the last fire, aligning with the intermediate disturbance hypothesis.Consequently, nurse plants can be considered a practical factor for developing active restoration practices in burned ecosystems.

Fig. 1 a
Fig. 1 a Characteristics of the study sites across different times since last fire.Within each study site, two transects 200 m in length and approximately 500 m apart were established, and 20 plots each nurse were placed on each ones.We sampled vegetation beneath nurse individuals and nearby paired open areas.b Five study sites were surveyed in the Fereizi-Chenaran region including (a) unburned site (i.e.control site without fire occurrence), (b) 1-year sites (i.e.site with most recent fire occurring in the last year), (c) 4-year site (i.e.site with most recent fire occurring 4 years prior), (d) 10-year site and (e) 20-year site (i.e.sites with last fire occurring between 10 and 20 years prior)

Fig. 2
Fig. 2 Mean ± SE (n = 20) of different nurse types (i.e.Elymus, Lonicera and Phlomis) Relative Interaction Index (RII) cover of the whole understory community at sites with different times since last fire (i.e.1-year, 4-year, 10-year and 20-year sites) and unburned site (control site).Significant treatment effects in two-way linear mixed-effects models on the effects of the time since last fire (TF), presence/absence of nurse type (N) treatments and their interactions (TF x N) are shown in the graphs p < 0.05; ***, p < 0.0001

Fig. 4
Fig. 4 Bar plots showing taxonomic (i.e.species richness (q0) and Shannon's diversity) (q1)) (above panels), functional (FSES.mpd)and phylogenetic (PSES.mpd)(below panels) diversity at the studied sites with respect to presence of nurse plants with grass (Elymus hispidens), herb (Phlomis cancellata) and shrub (Lonicera nummulariifolia) growth forms and different times since last fire.Variation in taxonomic, functional and phylogenetic indices across times since last fire and presence of nurse types were measured via linear mixed effect models.The marginal adjusted R 2 shows the variation explained by the fixed factors in our statistical models (mR 2 adj.are shown in each case).Significant treatment effects shown from two-way linear mixed-effects models on the effects of the time since last fire (TF), presence/absence of nurse type (N) treatments and their interactions (TF x N) are shown in the graphs.*, p < 0.05; **, p < 0.01; ***, p < 0.0001

Fig. 5
Fig. 5 The relative contribution of nurse types (NURSE in orange) and times since last fire (FIRE in yellow) to taxonomic (q0 and q1 indices), functional (FSES.mpd)and phylogenetic (PSES.mpd)diversity.Values represent the adjusted R 2 -values