scholarly journals The Influence of Trophic Interactions in the Plankton Community on Its Spatiotemporal Dynamics

Author(s):  
E.E. Giricheva

The present work considers a model of spatiotemporal dynamics for plankton community in a vertical column of water. Usually mesozooplankton (e.g. copepods) is considered as the main grazer in models of plankton systems. However, numerous studies have appeared recently that microzooplankton, not copepods, are the primary grazers in oceanic ecosystems. This work presents the model with microzooplankton and copepods as phytoplankton predators. Functional response for copepod feeding includes its preference for phytoplankton and microzooplankton. The effects of half saturation constants of zooplankton and diet variants on the stability of spatial homogeneous equilibriums was studied. In the case of effective grazing by mesozooplankton the system demonstrates stable coexistence of both zooplankton groups if phytoplankton carrying capacity is sufficient for it. In the case if microzooplankton protists are ineffective grazers – they are displaced by copepods or survive in a small number if their portion in the copepods diet is insignificant. The phenomenon of multistability has been identified in the system for copepod diet with microzooplankton preference. It means that several equilibria are stable at the same time, and the dynamics of the system are determined by the initial conditions. Allowing for spatial dependence in the modeled plankton community causes spatial structures appearance induced by Turing instability. It means that the system with diffusion only is able to induce stationary vertical profiles of plankton biomass. Such system state is possible if both predators are effective grazers and the microzooplankton prevails in the copepods diet.

2013 ◽  
Vol 2013 ◽  
pp. 1-14 ◽  
Author(s):  
Hunki Baek ◽  
Dong Ick Jung ◽  
Zhi-wen Wang

We investigate spatiotemporal dynamics of a semi-ratio-dependent predator-prey system with reaction-diffusion and zero-flux boundary. We obtain the conditions for Hopf, Turing, and wave bifurcations of the system in a spatial domain by making use of the linear stability analysis and the bifurcation analysis. In addition, for an initial condition which is a small amplitude random perturbation around the steady state, we classify spatial pattern formations of the system by using numerical simulations. The results of numerical simulations unveil that there are various spatiotemporal patterns including typical Turing patterns such as spotted, spot-stripelike mixtures and stripelike patterns thanks to the Turing instability, that an oscillatory wave pattern can be emerged due to the Hopf and wave instability, and that cooperations of Turing and Hopf instabilities can cause occurrence of spiral patterns instead of typical Turing patterns. Finally, we discuss spatiotemporal dynamics of the system for several different asymmetric initial conditions via numerical simulations.


2021 ◽  
Vol 3 (3) ◽  
Author(s):  
Satish Kumar Tiwari ◽  
Ravikant Singh ◽  
Nilesh Kumar Thakur

AbstractWe propose a model for tropic interaction among the infochemical-producing phytoplankton and non-info chemical-producing phytoplankton and microzooplankton. Volatile information-conveying chemicals (infochemicals) released by phytoplankton play an important role in the food webs of marine ecosystems. Microzooplankton is an ecologically important grazer of phytoplankton for coexistence of a large number of phytoplankton species. Here, we discuss how information transferred by dimethyl sulfide shapes the interaction of phytoplankton. Phytoplankton deterrents may lead to propagation of IPP bloom. The interaction between IPP and microzooplankton follows the Beddington–DeAngelis-type functional response. Analytically, we discuss boundedness, stability and Turing instability of the model system. We perform numerical simulation for temporal (ODE model) as well as a spatial model system. Our numerical investigation shows that microzooplankton grazing refuse of IPP leads to oscillatory dynamics. Increasing diffusion coefficient of microzooplankton shows Turing instability. Time evolution also plays an important role in the stability of system dynamics. The results obtained in this paper are useful to understand the dominance of algal bloom in coastal and estuarine ecosystem.


1992 ◽  
Vol 02 (01) ◽  
pp. 1-9 ◽  
Author(s):  
YOHANNES KETEMA

This paper is concerned with analyzing Melnikov’s method in terms of the flow generated by a vector field in contrast to the approach based on the Poincare map and giving a physical interpretation of the method. It is shown that the direct implication of a transverse crossing between the stable and unstable manifolds to a saddle point of the Poincare map is the existence of two distinct preserved homoclinic orbits of the continuous time system. The stability of these orbits and their role in the phenomenon of sensitive dependence on initial conditions is discussed and a physical example is given.


1999 ◽  
Author(s):  
O. Burak Ozdoganlar ◽  
William J. Endres

Abstract This paper presents a mathematical perspective, to complement the intuitive or practice-oriented perspective, to classifying machining operations as parallel-process (simultaneous) or single-process in nature. Illustrative scenarios are provided to demonstrate how these two perspectives may lead in different situations to the same or different conclusions regarding process parallelism. A model representation of a general parallel-process machining system is presented, based on which the general parallel-process stability eigenvalue problem is formulated. For a special simplified case of the general system, analytical methods are employed to derive a fully analytical stability solution. Thorough study of this solution through eigenvector analysis sheds light on some fundamental phenomena of parallel-process machining stability, such as dependence of the stability solution on phasing of the initial conditions (disturbances). This establishes the importance, when employing numerical time-domain simulation for such analyses, of specifying initial conditions for the multiple processes to be arbitrarily phased so that correct results are achieved across all spindle speeds.


2021 ◽  
Vol 31 (08) ◽  
pp. 2150143
Author(s):  
Zunxian Li ◽  
Chengyi Xia

In this paper, we explore the dynamical behaviors of the 1D two-grid coupled cellular neural networks. Assuming the boundary conditions of zero-flux type, the stability of the zero equilibrium is discussed by analyzing the relevant eigenvalue problem with the aid of the decoupling method, and the conditions for the occurrence of Turing instability and Hopf bifurcation at the zero equilibrium are derived. Furthermore, the approximate expressions of the bifurcating periodic solutions are also obtained by using the Hopf bifurcation theorem. Finally, numerical simulations are provided to demonstrate the theoretical results.


2014 ◽  
Vol 755 ◽  
pp. 705-731 ◽  
Author(s):  
Sasan Sarmast ◽  
Reza Dadfar ◽  
Robert F. Mikkelsen ◽  
Philipp Schlatter ◽  
Stefan Ivanell ◽  
...  

AbstractTwo modal decomposition techniques are employed to analyse the stability of wind turbine wakes. A numerical study on a single wind turbine wake is carried out focusing on the instability onset of the trailing tip vortices shed from the turbine blades. The numerical model is based on large-eddy simulations (LES) of the Navier–Stokes equations using the actuator line (ACL) method to simulate the wake behind the Tjæreborg wind turbine. The wake is perturbed by low-amplitude excitation sources located in the neighbourhood of the tip spirals. The amplification of the waves travelling along the spiral triggers instabilities, leading to breakdown of the wake. Based on the grid configurations and the type of excitations, two basic flow cases, symmetric and asymmetric, are identified. In the symmetric setup, we impose a 120° symmetry condition in the dynamics of the flow and in the asymmetric setup we calculate the full 360° wake. Different cases are subsequently analysed using dynamic mode decomposition (DMD) and proper orthogonal decomposition (POD). The results reveal that the main instability mechanism is dispersive and that the modal growth in the symmetric setup arises only for some specific frequencies and spatial structures, e.g. two dominant groups of modes with positive growth (spatial structures) are identified, while breaking the symmetry reveals that almost all the modes have positive growth rate. In both setups, the most unstable modes have a non-dimensional spatial growth rate close to $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\pi /2$ and they are characterized by an out-of-phase displacement of successive helix turns leading to local vortex pairing. The present results indicate that the asymmetric case is crucial to study, as the stability characteristics of the flow change significantly compared to the symmetric configurations. Based on the constant non-dimensional growth rate of disturbances, we derive a new analytical relationship between the length of the wake up to the turbulent breakdown and the operating conditions of a wind turbine.


2006 ◽  
Vol 6 (6) ◽  
pp. 13307-13321
Author(s):  
S. Ceccherini ◽  
C. Belotti ◽  
B. Carli ◽  
P. Raspollini ◽  
M. Ridolfi

Abstract. The retrieval of concentration vertical profiles of atmospheric constituents from spectroscopic measurements is often an ill-conditioned problem and regularization methods are frequently used to improve its stability. Recently a new method, that provides a good compromise between precision and vertical resolution, was proposed to determine analytically the value of the regularization parameter. This method is applied for the first time to real measurements with its implementation in the operational retrieval code of the satellite limb-emission measurements of the MIPAS instrument and its performances are quantitatively analyzed. The adopted regularization improves the stability of the retrieval providing smooth profiles without major degradation of the vertical resolution. In the analyzed measurements the retrieval procedure provides a vertical resolution that, in the troposphere and low stratosphere, is smaller than the vertical field of view of the instrument.


2009 ◽  
Vol 23 (10) ◽  
pp. 2303-2322 ◽  
Author(s):  
XIAO-FENG PANG ◽  
MEI-JIE LIU

The dynamic features of soliton transporting the bio-energy in the α-helix protein molecules with three channels under influences of temperature of systems and chain–chain interaction among these channels have been numerically studied by using the dynamic equations in a new model and the fourth-order Runge–Kutta method. This result obtained shows that the chain–chain interaction depresses the stability of the soliton due to the dispersed effect, but the stability of the soliton in the case of simultaneous motivation of three channels by an initial conditions is better than that in another initial condition. We also find from this investigation that the new soliton can transport steadily over 1000 amino acid residues in the cases of motion of long time of 120 ps, and retain their shapes and energies to travel towards the protein molecules after mutual collision of the solitons at the biological temperatures of 300 K. Therefore the soliton is very robust against the thermal perturbation of the α-helix protein molecules at 300 K. From the investigation of changes of features of the soliton with increasing temperature, we find that the amplitudes and velocities of the solitons decrease with increasing temperature of proteins, but the soliton disperses in the cases of higher temperature of 325 K and larger structure disorders. Thus we find that the critical temperature of the soliton occurring in the α-helix protein molecules is about 320 K. Therefore we can conclude that the soliton in the new model can play an important role in the bio-energy transport in the α-helix protein molecules with three channels at biological temperature, and the new model is possibly a candidate for the mechanism of this transport.


2021 ◽  
Vol 31 (15) ◽  
Author(s):  
Penghe Ge ◽  
Hongjun Cao

The existence of chaos in the Rulkov neuron model is proved based on Marotto’s theorem. Firstly, the stability conditions of the model are briefly renewed through analyzing the eigenvalues of the model, which are very important preconditions for the existence of a snap-back repeller. Secondly, the Rulkov neuron model is decomposed to a one-dimensional fast subsystem and a one-dimensional slow subsystem by the fast–slow dynamics technique, in which the fast subsystem has sensitive dependence on the initial conditions and its snap-back repeller and chaos can be verified by numerical methods, such as waveforms, Lyapunov exponents, and bifurcation diagrams. Thirdly, for the two-dimensional Rulkov neuron model, it is proved that there exists a snap-back repeller under two iterations by illustrating the existence of an intersection of three surfaces, which pave a new way to identify the existence of a snap-back repeller.


Sign in / Sign up

Export Citation Format

Share Document