Complexity-multistability relationships: How does species diversity shape the number of alternative stable states?

Ecologists have long investigated how community complexity affects ecological stability, yet how community complexity influences multistability, defined as the presence of alternative stable states, remains poorly understood. We developed a novel framework integrating stochastic community assembly with stability landscape analysis to quantify multistability from species interaction matrices. Using this framework, we systematically explored how species interaction properties shape the relationship between species diversity (species pool size) and the number of alternative stable states. Mean interaction strength was the primary determinant: competitive interactions amplified the positive relationship between species diversity and the number of alternative stable states. In competitive communities, a greater number of alternative stable states was associated with lower community uncertainty, a measure of the long-term unpredictability of community assembly dynamics. These results highlight the importance of characterizing the entire stability landscape. Our framework provides a general approach for understanding and quantifying multistability in complex ecological communities.