Timescale separation and models of symbiosis: state space reduction, multiple attractors and initialization

Dynamic Energy Budget models relate whole organism processes such as growth, reproduction and mortality to suborganismal metabolic processes. Much of their potential derives from extensions of the formalism to describe the exchange of metabolic products between organisms or organs within a single organism, for example the mutualism between corals and their symbionts. Without model simplification, such models are at risk of becoming parameter-rich and hence impractical. One natural simplification is to assume that some metabolic processes act on 'fast' timescales relative to others. A common strategy for formulating such models is to assume that 'fast' processes equilibrate immediately, while 'slow' processes are described by ordinary differential equations. This strategy can bring a subtlety with it. What if there are multiple, interdependent fast processes that have multiple equilibria, so that additional information is needed to unambiguously specify the model dynamics? This situation can easily arise in contexts where an organism or community can persist in a 'healthy' or an 'unhealthy' state with abrupt transitions between states possible. To approach this issue, we offer the following: (a) a method to unambiguously complete implicitly defined models by adding hypothetical 'fast' state variables; (b) an approach for minimizing the number of additional state variables in such models, which can simplify the numerical analysis and give insights into the model dynamics; and (c) some implications of the new approach that are of practical importance for model dynamics, e.g. on the bistability of flux dynamics and the effect of different initialization choices on model outcomes. To demonstrate those principles, we use a simplified model for root-shoot dynamics of plants and a related model for the interactions between corals and endosymbiotic algae that describes coral bleaching and recovery.