Analyzing minimum viable populations in deterministic community models using viability space decomposition

Minimum viable populations (MVPs) are population levels large enough to surmount risk from demographic, environmental, and genetic stochasticity. MVPs are estimated by biologists to guide conservation practices. However, MVPs are generally estimated for a target population without regard for how they interact with intra- and inter-species population dynamics in the broader ecological community. Thus, how and why population dynamics interact with MVPs imposed by conservation biologists remain unclear. When MVPs are imposed on a continuous population model, traditional analyses fail to capture the range of possible outcomes those MVPs create. Here, we describe viability space decomposition (VSD) as a mathematical tool to systematically analyze the potential crossing of MVPs during population dynamics. We demonstrate that different extinction and survival outcomes can be recovered from a model with imposed MVPs using three VSD concepts in junction with a traditional phase portrait: mortality manifolds which separate conditions that lead to different existential outcomes, ordering manifolds which determine the order of extinction events for multiple populations, and collapse manifolds which determine the survival or extinction of one species given the loss of another. We employ these methods with a standard consumer-resource model, and the methods can be scaled to systems with more species. VSD is a useful tool for conservation biologists and community ecologists concerned with boundary crossing problems in any dynamical system.