Species-specific traits plus stabilizing processes best explain coexistence in biodiverse fire-prone plant communities

Coexistence in fire-prone Mediterranean-type shrublands has been explored in the past using both neutral and niche-based models. However, distinct differences between plant functional types (PFTs), such as fire-killed vs resprouting responses to fire, and the relative similarity of species within a PFT, suggest that coexistence models might benefit from combining both neutral and niche-based (stabilizing) approaches. We developed a multispecies metacommunity model where species are grouped into two PFTs (fire-killed vs resprouting) to investigate the roles of neutral and stabilizing processes on species richness and rank-abundance distributions. Our results show that species richness can be maintained in two ways: i) strictly neutral species within each PFT, or ii) species within PFTs differing in key demographic properties, provided that additional stabilizing processes, such as negative density regulation, also operate. However, only simulations including stabilizing processes resulted in structurally realistic rank-abundance distributions over plausible time scales. This result underscores the importance of including both key species traits and stabilizing (niche) processes in explaining species coexistence and community structure.