Global redistribution and local migration in semi-discrete host-parasitoid population dynamic models

Fundamental limits of parasitoid-driven host population suppression: Implications for biological control

Parasitoid wasps are increasingly being used to control insect pest populations, where the pest is the host species parasitized by the wasp. Here we use the discrete-time formalism of the Nicholson-Bailey model to investigate a fundamental question-are there limits to parasitoid-driven suppression of the host population density while still ensuring a stable coexistence of both species? Our model formulation imposes an intrinsic self-limitation in the host's growth resulting in a carrying capacity in the absence of the parasitoid. Different versions of the model are considered with parasitism occurring at a developmental stage that is before, during, or after the growth-limiting stage. For example, the host's growth limitation may occur at its larval stage due to intraspecific competition, while the wasps attack either the host egg, larval or pupal stage. For slow-growing hosts, models with parasitism occurring at different life stages are identical in terms of their host suppression dynamics but have contrasting differences for fast-growing hosts. In the latter case, our analysis reveals that wasp parasitism occurring after host growth limitation yields the lowest pest population density conditioned on stable host-parasitoid coexistence. For ecologically relevant parameter regimes we estimate this host suppression to be roughly 10-20% of the parasitoid-free carrying capacity. We further expand the models to consider a fraction of hosts protected from parasitism (i.e., a host refuge). Our results show that for a given host reproduction rate there exists a critical value of protected host fraction beyond which, the system dynamics are stable even for high levels of parasitism that drive the host to arbitrary low population densities. In summary, our systematic analysis sheds key insights into the combined effects of density-dependence in host growth and parasitism refuge in stabilizing the host-parasitoid population dynamics with important implications for biological control.

A comparative approach to stabilizing mechanisms between discrete- and continuous-time consumer-resource models

There is rich literature on using continuous-time and discrete-time models for studying population dynamics of consumer-resource interactions. A key focus of this contribution is to systematically compare between the two modeling formalisms the stabilizing/destabilizing impacts of diverse ecological processes that result in a density-dependent consumer attack rate. Inspired by the Nicholson-Bailey/Lotka-Volterra models in discrete-time/continuous-time, respectively, we consider host-parasitoid interactions with an arbitrary parasitoid attack rate that is a function of both the host/parasitoid population densities. Our analysis shows that a Type II functional response is stabilizing in both modeling frameworks only when combined with other mechanisms, such as mutual interference between parasitoids. A Type III functional response is by itself stabilizing, but the extent of attack-rate acceleration needed is much higher in the discrete-time framework, and its stability regime expands with increasing host reproduction. Finally, our results show that while mutual parasitoid interference can stabilize population dynamics, cooperation between parasitoids to handle hosts is destabilizing in both frameworks. In summary, our comparative analysis systematically characterizes diverse ecological processes driving stable population dynamics in discrete-time and continuous-time consumer-resource models.

Attack by a common parasitoid stabilizes population dynamics of multi-host communities

We model the population dynamics of two host species attacked by a common parasitoid using a discrete-time formalism that captures their population densities from year to year. It is well known starting from the seminal work of Nicholson and Bailey that a constant parasitoid attack rate leads to an unstable host-parasitoid interaction. However, a Type III functional response, where the parasitoid attack rate accelerates with increasing host density stabilizes the population dynamics. We first consider a scenario where both host species are attacked by a parasitoid with the same Type III functional response. Our results show that sufficient fast acceleration of the parasitoid attack rate stabilizes the population dynamics of all three species. For two symmetric host species, the extent of acceleration needed to stabilize the three-species equilibrium is exactly the same as that needed for a single host-parasitoid interaction. However, asymmetry can lead to scenarios where the removal of a host species from a stable interaction destabilizes the interaction between the remaining host species and the parasitoid. Next, we consider a situation where one of the host species is attacked at a constant rate (i.e., Type I functional response), and the other species is attacked via a Type III functional response. We identify parameter regimes where a Type III functional response to just one of the host species stabilizes the three species interaction. In summary, our results show that a generalist parasitoid with a Type III functional response to one or many host species can play a key role in stabilizing population dynamics of host-parasitoid communities in apparent competition.

Stochasticity in host-parasitoid models informs mechanisms regulating population dynamics

Population dynamics of host-parasitoid interactions have been traditionally studied using a discrete-time formalism starting from the classical work of Nicholson and Bailey. It is well known that differences in parasitism risk among individual hosts can stabilize the otherwise unstable equilibrium of the Nicholson-Bailey model. Here, we consider a stochastic formulation of these discrete-time models, where the host reproduction is a random variable that varies from year to year and drives fluctuations in population densities. Interestingly, our analysis reveals that there exists an optimal level of heterogeneity in parasitism risk that minimizes the extent of fluctuations in the host population density. Intuitively, low variation in parasitism risk drives large fluctuations in the host population density as the system is on the edge of stability. In contrast, high variation in parasitism risk makes the host equilibrium sensitive to the host reproduction rate, also leading to large fluctuations in the population density. Further results show that the correlation between the adult host and parasitoid densities is high for the same year, and gradually decays to zero as one considers cross-species correlations across different years. We next consider an alternative mechanism of stabilizing host-parasitoid population dynamics based on a Type III functional response, where the parasitoid attack rate accelerates with increasing host density. Intriguingly, this nonlinear functional response makes qualitatively different correlation signatures than those seen with heterogeneity in parasitism risk. In particular, a Type III functional response leads to uncorrelated adult and parasitoid densities in the same year, but high cross-species correlation across successive years. In summary, these results argue that the cross-correlation function between population densities contains signatures for uncovering mechanisms that stabilize consumer-resource population dynamics.

Stochastic dynamics of predator-prey interactions

The interaction between a consumer (such as, a predator or a parasitoid) and a resource (such as, a prey or a host) forms an integral motif in ecological food webs, and has been modeled since the early 20th century starting from the seminal work of Lotka and Volterra. While the Lotka-Volterra predator-prey model predicts a neutrally stable equilibrium with oscillating population densities, a density-dependent predator attack rate is known to stabilize the equilibrium. Here, we consider a stochastic formulation of the Lotka-Volterra model where the prey's reproduction rate is a random process, and the predator's attack rate depends on both the prey and predator population densities. Analysis shows that increasing the sensitivity of the attack rate to the prey density attenuates the magnitude of stochastic fluctuations in the population densities. In contrast, these fluctuations vary non-monotonically with the sensitivity of the attack rate to the predator density with an optimal level of sensitivity minimizing the magnitude of fluctuations. Interestingly, our systematic study of the predator-prey correlations reveals distinct signatures depending on the form of the density-dependent attack rate. In summary, stochastic dynamics of nonlinear Lotka-Volterra models can be harnessed to infer density-dependent mechanisms regulating predator-prey interactions. Moreover, these mechanisms can have contrasting consequences on population density fluctuations, with predator-dependent attack rates amplifying stochasticity, while prey-dependent attack rates countering to buffer fluctuations.