Search

In the current manuscript, we consider a predator-prey model where the predator is modeled as a generalist using a modified Leslie-Gower scheme, and the prey exhibits group defence via a generalised response. We show that the model could exhibit finite time blow-up, contrary to the current literature (Eur. Phys. J. Plus 137, 28). We also propose a new concept via which the predator population blows up in finite time while the prey population quenches in finite time. The blow-up and quenching times are proved to be one and the same. Our analysis is complemented by numerical findings. This includes a numerical description of the basin of attraction for large data blow-up solutions, as well as several rich bifurcations leading to multiple limit cycles, both in co-dimension one and two. Lastly, we posit a delayed version of the model with globally existing solutions for any initial data.

Population ecology theory is replete with density dependent processes. However trait-mediated or behavioral indirect interactions can both reinforce or oppose density-dependent effects. This paper presents the first two species competitive ODE and PDE systems where an Allee effect, which is a density dependent process and the fear effect, which is non-consumptive and behavioral are both present. The stability of the equilibria is discussed analytically using the qualitative theory of ordinary differential equations. It is found that the Allee effect and the fear effect change the extinction dynamics of the system and the number of positive equilibrium points, but they do not affect the stability of the positive equilibria. We also observe some special dynamics that induce bifurcations in the system by varying the Allee or fear parameter. Interestingly we find that the Allee effect working in conjunction with the fear effect, can bring about several qualitative changes to the dynamical behavior of the system with only the fear effect in place, in regimes of small fear. That is, for small amounts of the fear parameter, it can change a competitive exclusion type situation to a strong competition type situation. It can also change a weak competition type situation to a bi-stability type situation. However for large fear regimes the Allee effect reinforces the dynamics driven by the fear effect. The analysis of the corresponding spatially explicit model is also presented. To this end the comparison principle for parabolic PDE is used. The conclusions of this paper have strong implications for conservation biology, biological control as well as the preservation of biodiversity.

We consider the effect of finite time extinction mechanisms (FTEM) such as (1) semi-linear harvesting terms, and (2) quasi-linear fast diffusion terms on two species Lokta-Volterra competition models. We show that these mechanisms can alter classical dynamics of competitive exclusion, and weak and strong competition by acting only on a \emph{small} portion of the weaker competitors' population, analogous to small defector populations in game theory \cite{DC23}. In particular, a stronger competitors population, with a few individuals dispersing (``defecting") very quickly, could exhibit bi-stability, as well as competitive exclusion \emph{reversal}. The non-linear harvesting is applied to aphid-soybean crop systems, wherein novel dynamics are observed. Applications to bio-control of invasive pests such as the soybean aphid are discussed., Comment: 30 pages

The literature posits that an introduced predator population, is able to drive it's target pest population extinct, if supplemented with high quality additional food of quantity $ξ> ξ_{critical}$, \cite{SP11, SPV18, SPD17, SPM13}. We show this approach leads to infinite time blow-up of the predator population. We propose an alternate model in which the additional food induces predator competition. Analysis herein indicates that there are threshold values $c^{*}_{1} < c^{*}_{2} < c^{*}_{3}$ of the competition parameter $c$, s.t. when $c < c^{*}_{1}$, the pest free state is globally stable, when $c^{*}_{2} < c < c^{*}_{3}$, bi-stability is possible, and when $c^{*}_{3} < c$, up to three interior equilibriums could exist. As $c$ and $ξ$-$c$ are varied, standard co-dimension one and co-dimension two bifurcations are observed. The recent dynamical systems literature involving predator competition, report several non-standard bifurcations such as the saddle-node-transcritical bifurcation (SNTC) occurring in co-dimension two \cite{KSV10, BS07}, and cusp-transcritical bifurcation (CPTC) in co-dimension three, \cite{D20, BS07}. We show that in our model structural symmetries can be exploited to construct a SNTC in co-dimension two, and a CPTC also in co-dimension two. We further use these symmetries to construct a novel pitchfork-transcritical bifurcation (PTC) in co-dimension two, thus explicitly characterizing a new organizing center of the model. Dynamics such as homoclinic orbits, concurrently occurring limit cycles, and competition driven Turing patterns are also observed. Our findings indicate that increasing additional food in predator-pest models, can hinder bio-control, contrary to some of the literature. However, additional food that also induces predator competition, leads to novel bio-control scenarios, and complements the work in \cite{H21, B98, K04, D20, BS07, VH19}., 30 pages

We investigate the classical two species ODE and PDE Lotka-Volterra competition models, where one of the competitors could potentially go extinct in finite time. We show that in this setting, classical theories and intuitions do not hold, and various counter intuitive dynamics are possible. In particular, the weaker competitor could avoid competitive exclusion, and the slower diffuser may not win. Numerical simulations are performed to verify our analytical findings., 21 pages, 27 figures

The Trojan Y Chromosome strategy (TYC) is a genetic biocontrol strategy designed to alter the sex ratio of a target invasive population by reducing the number of females over time. Recently an alternative strategy is introduced, that mimics the TYC strategy by harvesting females whilst stocking males . We consider the FHMS strategy, with a weak Allee effect, and show that the extinction boundary need note be hyperbolic. To the best of our knowledge, this is the first example of a non-hyperbolic extinction boundary in mating models, structured by sex. Next, we consider the spatially explicit model and show that the weak Allee effect is both sufficient and necessary for Turing patterns to occur. We discuss the applicability of our results to large scale biocontrol, as well as compare and contrast our results to the case with a strong Allee effect.

The Trojan Y Chromosome Strategy (TYC) is an extremely well investigated biological control method for controlling invasive populations with an XX-XY sex determinism. In \cite{GP12, WP14} various dynamical properties of the system are analyzed, including well posedness, boundedness of solutions, and conditions for extinction or recovery. These results are derived under the assumption of positive solutions. In the current manuscript, we show that if the introduction rate of trojan fish is zero, under certain large data assumptions, negative solutions are possible for the male population, which in turn can lead to finite time blow-up in the female and male populations. A comparable result is established for \emph{any} positive initial condition if the introduction rate of trojan fish is large enough. Similar finite time blow-up results are obtained in a spatial temporal TYC model that includes diffusion. Lastly, we investigate improvements to the TYC modeling construct that may dampen the mechanisms to the blow-up phenomenon or remove the negativity of solutions. The results draw into suspect the reliability of current TYC models under certain situations., 18 pages, 8 figures

In the current manuscript, an attempt has been made to understand the dynamics of a time-delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis type functional responses for large initial data. In \cite{RK15}, we have seen that the model does possess globally bounded solutions, for small initial conditions, under certain parametric restrictions. Here, we show that actually solutions to this model system can blow-up in finite time, for large initial condition, \emph{even} under the parametric restrictions derived in \cite{RK15}. We prove blow-up in the delayed model, as well as the non delayed model, providing sufficient conditions on the largeness of data, required for finite time blow-up. Numerical simulations show, that actually the initial data does not have to be very large, to induce blow-up. The spatially explicit system is seen to possess Turing instability. We have also studied Hopf-bifurcation direction in the spatial system, as well as stability of the spatial Hopf-bifurcation using the central manifold theorem and normal form theory.

In this work we consider a four compartment Brusselator system. The reaction terms of this system are of non constant sign, thus components of the solution are not bounded apriori, and functional means to derive apriori bounds will fail. We prove global existence of solutions, via construction of an appropriate lyapunov functional. Furthermore due to the sign changing nonlinearities, the asymptotic sign condition is also not satisfied, causing further difficulties in proving the existence of global attractor. These difficulties are circumvented via the use of the lyapunov functional constructed along with the use of the uniform gronwall lemma. We are able to prove the existence of an attractor for the system, improving previous results in the literature.The Hausdorff and fractal dimensions of the attractor are also shown to be finite.In particular we derive a lower bound on the Hausdorff dimension of the global attractor. We use numerical simulations,as well as numerical attractor reconstruction methods via nonlinear time series analysis, to validate our results.

Systems of reaction-diffusion equations are commonly used in biological models of food chains. The populations and their complicated interactions present numerous challenges in theory and in numerical approximation. In particular, self-diffusion is a nonlinear term that models overcrowding of a particular species. The nonlinearity complicates attempts to construct efficient and accurate numerical approximations of the underlying systems of equations. In this paper, a new nonlinear splitting algorithm is designed for a partial differential equation that incorporates self-diffusion. We present a general model that incorporates self-diffusion and develop a numerical approximation. The numerical analysis of the approximation provides criteria for stability and convergence. Numerical examples are used to illustrate the theoretical results.

While language competition models of diachronic language shift are increasingly sophisticated, drawing on sociolinguistic components like variable language prestige, distance from language centers and intermediate bilingual transitionary populations, in one significant way they fall short. They fail to consider contact-based outcomes resulting in mixed language practices, e.g. outcome scenarios such as creoles or unmarked code switching as an emergent communicative norm. On these lines something very interesting is uncovered in India, where traditionally there have been monolingual Hindi speakers and Hindi/English bilinguals, but virtually no monolingual English speakers. While the Indian census data reports a sharp increase in the proportion of Hindi/English bilinguals, we argue that the number of Hindi/English bilinguals in India is inaccurate, given a new class of urban individuals speaking a mixed lect of Hindi and English, popularly known as "Hinglish". Based on predator-prey, sociolinguistic theories, salient local ecological factors and the rural-urban divide in India, we propose a new mathematical model of interacting monolingual Hindi speakers, Hindi/English bilinguals and Hinglish speakers. The model yields globally asymptotic stable states of coexistence, as well as bilingual extinction. To validate our model, sociolinguistic data from different Indian classes are contrasted with census reports: We see that purported urban Hindi/English bilinguals are unable to maintain fluent Hindi speech and instead produce Hinglish, whereas rural speakers evidence monolingual Hindi. Thus we present evidence for the first time where an unrecognized mixed lect involving English but not "English", has possibly taken over a sizeable faction of a large global population., This paper has been withdrawan as the model has now been modified and the existing model has some errors

In this work we use functional methods to prove the boundedness and global existence of solutions for a class of strongly coupled parabolic systems. We apply the results to deduce the global existence of solutions for a classic Shigesada-Kawasaki-Teramoto (SKT) type model, for an extended range of the self-diffusion and cross-diffusion coefficients than those available in the current literature. We provide numerical simulations in 2D, via a spectral Galerkin method to verify our global existence results, as well as to visualize the dynamics of the system.

We consider nonlinear parabolic SPDEs of the form $\partial_t u=-(-\Delta)^{\alpha/2} u + b(u) +\sigma(u)\dot w$, where$\dot w$ denotes space-time white noise. The functions $b$ and $\sigma$ are both locally Lipschitz continuous. Under some suitable conditions on the parameters of these SPDEs, we show that the second moment of their solutions blow up in finite time. This complements recent works of Khoshnevisan and his coauthors; see for instance Foondun and Khoshnevisan (2009), Foondun and Khoshnevisan (to appear, 2012) and Conus, Joseph, and Khoshnevisan, as well as those of Chow (Chow, 2009, and Chow, 2011). Furthermore, upon comparing our stochastic equations with their deterministic counterparts, we find that our results indicates that the presence of noise might affect the occurrence of blow-up.