Gibbs, T; Grilli, J; Rogers, T; Allesina, S

Random matrix theory successfully connects the structure of interactions of large ecological communities to their ability to respond to perturbations. One of the most debated aspects of this approach is that so far studies have neglected the role of population abundances on stability. While species abundances are well studied and empirically accessible, studies on stability have so far failed to incorporate this information. Here we tackle this question by explicitly including population abundances in a random matrix framework. We derive an analytical formula that describes the spectrum of a large community matrix for arbitrary feasible species abundance distributions. The emerging picture is remarkably simple: while population abundances affect the rate to return to equilibrium after a perturbation, the stability of large ecosystems is uniquely determined by the interaction matrix. We confirm this result by showing that the likelihood of having a feasible and unstable solution in the Lotka-Volterra system of equations decreases exponentially with the number of species for stable interaction matrices.