1. Motivation The growth and degradation of populations in the nature and the struggle of one species to dominate other species have been an interesting topic for a long time. The application of mathematical concepts to explain these phenomena have been documented for centuries. The founders of mathematical-based modeling are Malthus (1798), Verhulse (1838), Pearl and Reed (1903), especially Lotka and Volterra whose most important results were published in the 1920s and 1930s. Lotka and Volterra modeled, independently of each other, the competition between predator and prey. Their work has important meaning for the population biology eld. They are the rst to study the phenomenon of species interactions by introducing simpli ed conditions that lead to solvable problems that have meaning until today. The proposed model is given by where N i 8 > > < > > : dN (t) dt dN 1 (t) dt 2 = r = r 1 2 N N 1 2 (t) (t) 1 N 1 N (t) K 1 2 1 (t) K 2 a a 12 21 N (t) K N 2 1 1 (t) K 2 (t) is the density of species i; i 2 f1; 2g, at time t. Parameters r ; ; are the growth rate and the carrying capacity of the species i; i 2 f1; 2g, respectively. Parameter a is interspeci c competitive coecient representing the negative e ect of species j on the growth of species i, i = j; i; j 2 f1; 2g. ij The ecological meaning of this model is that two species coexist only if the e ects of their competition are small. When the competing e ects of two species are large, one of the two species will be extinct. This famous principle is called the competitive exclusion principle. Today, this model is still applicable to competitions between a number of biophysical species in practice and in empirical observations [22, 63]. However, there are many other competing bio-systems, which cannot be explained by using the classic competition model of Lotka-Volterra (or the competitive exclusion principle). We present here two examples. In the rst example, Atkinson and Shorrocks [12] studied the competition of two species for having phytoplankton (food) in multiple environments. Competition is noted when one of the two species is absent, resulting in an increase the remaining species. Although the measured competing e ect is signi cant, the two species coexist. This result is contrary to the exclusion principle of the classical competition model. In the second example, Lei i and K iand Hanski [63] studied two species of parasites on the same Melitaeacinxia butter y. The results showed that the more competitive and less hostile species (Cotesia melitaearum) were not founded in some host species, while the less competitive (Hyposotherapy horticola) were found in all hosts of Melitaeacinxia. This result is also contrary to the principle of competitive exclusion. The main reason for the limitation of Lotka-Volterra''s classic competition model is that there are too many assumptions in the model, such as the assumption that the environment is homogeneous and stable (expressed by the carrying capacities K for the specie i; i 2 f1; 2g), the behavior of the individual species is the same and the competition is expressed only by interspeci c competitive coecient a i . Meanwhile, these factors appear frequently and play a very important role. For example, the migration behavior of individual species is a very important factor for species survival [80, 104]. Individuals of the same species or of di erent species may have di erent behaviors. Aggressive behavior is also used by individuals of wild species to compete for accommodation, to ght against their partners, etc. In addition, individuals may also change their behaviors frequently according to the change of the environment as studied in [110, 111]. Therefore, the development of new models that take into account the complex environments and the behaviors of individuals has been interested by many mathematicians. Following are some recent approaches. The complex environment and individual migration behavior in competitive ecosystems. The competition process and the migration process have the same time scale or di erent time scales. Aggressive behavior of individuals in competitive system. Age structure (adult group and immature group) in the competitive system. ij
MINISTRY OF EDUCATION AND TRAINING HANOI UNIVERSITY OF SCIENCE AND TECHNOLOGY ——————————- NGUYEN PHUONG THUY COMPETITIVE ECOSYSTEMS: CONTINUOUS AND DISCRETE MODELS DOCTORAL DISSERTATION OF MATHEMATICS HANOI - 2018 Contents DECLARATION OF AUTHORSHIP iii iv LIST OF ABBREVIATIONS LIST OF FIGURES ACKNOWLEDGEMENTS LIST OF TABLES INTRODUCTION LITERATURE REVIEW 10 1.1 Competition in ecology systems 10 1.2 Continuous models 11 1.3 Discrete models 13 1.4 Lyapunov’s methods and LaSalle’s invariance principle 16 1.5 Aggregation method 18 CONTINUOUS MODELS FOR COMPETITIVE SYSTEMS WITH 21 STRATEGY 2.1 Introduction on competitive systems 21 2.2 The classical competition model without individuals’ strategy 24 2.3 A model with an avoiding strategy 25 2.4 A model with an aggressive strategy 2.5 Discussion and Conclusion 39 32 DISCRETE MODELS FOR PREDATOR-PREY SYSTEMS 46 3.1 Introduction 46 3.2 Individual-based predator-prey model 47 3.3 Generating graph of the individual-based predator-prey model 50 3.4 3.3.1 Graph model for complex systems 50 3.3.2 Graph model for predator-prey system 52 3.3.3 Analysis of the generating graph 53 Conclusion and Perspectives 54 i APPLICATION: MODELING OF SOME REFERENCE ECOSYS57 TEMS 4.1 4.2 Modeling of the thiof-octopus system 57 4.1.1 Introduction 57 4.1.2 Model presentation 59 4.1.3 Analysis and Discussion 69 Modeling the brown plant-hopper system 74 4.2.1 Introduction 75 4.2.2 Modeling 4.2.3 Analysis and Discussion 79 76 CONCLUSION 95 BIBLIOGRAPHY 97 LIST OF PUBLICATIONS 107 ii LIST OF ABBREVIATIONS EBM : Equation-Based Model IBM : Individual-Based Model GBM : Graph-Based Model LSE : Local Superior resource Exploiter LIE : Local Inferior resource Exploiter BPH : Brown Plant Hopper LIST OF TABLES Table 2.1 Equilibria of aggregated model (2.13) and local stability analysis 40 Table 3.1 The statistics for several complex systems Table 3.2 The statistics for several steps of the simulation of the predator- 51 prey competition system 54 Table 3.3 Statistics about the cliques of the graphs at step of the simulation of the predator-prey competition system 55 Table 3.4 Statistics about the cliques of the graphs at step 530 of the simulation of the predator-prey competition system 55 LIST OF FIGURES Figure 1.1 Principle of equation-based modeling N1 and N2 are variables (compartments) F is the mathematical function which represents general laws applied to all members of the compartments [83] 13 Figure 1.2 Principle of individual-based modeling [83] 14 Figure 1.3 Principle of disk graph-based modeling [83] 15 Figure 2.1 Comparison of solutions of system (2.3) with their approximations through the aggregated system (2.10) for the both biotic and abiotic resource cases This figure shows the evolutions in time of each of the four state variables of system (2.3) (R, C1 , C1C and C1R ) and their approximations obtained from the aggregated system (2.10) (R , C1 , kC2 /H(C1 ) and (αC1 + α0 )C2 /H(C1 ) , respectively), for the same parameter values (r = 3; K = 20; S = 20; a1 = 0.8; e1 = 0.1; a2 = 0.6; e2 = 0.2; d1 = 0.4; d2C = 0.8; d2N = 0.8; α = 1.5; α0 = and k = 1) and initial conditions R(0) = 30; C1 (0) = 20; C2C (0) = 15 and C2R (0) = 10 31 Figure 2.2 Comparison of solutions of system (2.11) with their approximations through the aggregated system (2.13) for the both biotic and abiotic resource cases This figure shows the evolutions in time of each of the four state variables of system (2.11) (R, C1C , C1R and C2 ) and their approximations obtained from the aggregated system (2.13) (R, mC1 /L(C2 ), (βC2 + β0 )C1 /L(C2 ) and C2 , respectively), for the same parameter values (r = 5; K = 7; S = 7; a1 = 0.9; e1 = 0.1; a2 = 0.7; e2 = 0.2; d2 = 0.5; d1C = 0.2; d1N = 0.2; β = 5; α0 = 1, l = 0.2 and m = 0.4) and initial conditions R(0) = 30; C2 (0) = 20; C1C (0) = 15 and C1N (0) = 10 Figure 2.3 34 The outcomes of model (2.11) with the biotic resource 41 Figure 2.4 The outcomes of model (2.11) with the abiotic resource In each corresponding simulation, parametershave the same values as in the case of biotic resource and the values of S and the values of K are exactly the same 42 Figure 2.5 The left panel is about domains of the space (l, d1N , β) for the different outcomes of model 2.13 of the abiotic resource case Domain (I): LIE wins, domain (II): extinction, domain (III): LSE wins and domain (IV): exclusion via priority effects 43 Figure 2.6 The left panel is about domains of the space (l, d1N , β0 ) for the different outcomes of model 2.13 of the biotic resource case Domain (I): LIE wins, domain (II): extinction, domain (III): LSE wins and domain (IV): exclusion via priority effects 44 Figure 3.1 Species individual behavior at each simulation step 48 Figure 3.2 Distribution of individuals in several simulation steps Red, blue and green grid cells represent respectively Predator, Prey and Grass individuals a) at step 10, b) at step 100, c) at step 200, d) at step 300 Figure 3.3 49 Evolution of the number of individuals of each species The red, blue and green curves represent respectively the evolution of Predator, Prey and Grass 50 Figure 3.4 Individual Based Model (on the left) and the corresponding Disk Graph Based Model (on the right) 54 Figure 3.5 Distribution of degree in several simulation steps: a) at step 1, b) at step 530, c) at step 1000, d) at step 2500 56 Figure 4.1 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.7; r2 = 1.3; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3; E = 0.9 61 Figure 4.2 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.9; r2 = 0.7; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3; E = 0.9; d1 = 0.2; d2 = 0.4; k = 5; k = 7; m = 6; m = 0.2 64 Figure 4.3 Example of the case where the inferior competitor wins globally in model 2: A comparison between the aggregated model (blue dots) and the complete model (red curve) The parameters are the same as in Figure 4.2 except for ε = 0.01 65 Figure 4.4 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.9; r2 = 1.3; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3; E = 0.9; d1 = 0.2; d2 = 0.4; k = 5; k = 7; m = 6; α = 1; α0 = Figure 4.5 68 Example of the case where the inferior competitor wins globally in model A comparison between the aggregated model (blue dots) and the complete model (red curve) The parameters are the same as in Figure 4.4 except for ε = 0.01 68 Figure 4.6 Two cases where there exists a strictly positive equilibrium: (a) the case where (n∗1 , n∗2 ) is stable, (b) the case where (n∗1 , n∗2 ) is saddle 71 Figure 4.7 A photo of BPH-the predator of rice 76 Figure 4.8 Rice and brown plant-hopper system ni is the densities of rice respectively in patch i, i ∈ {1, 2} piA , piJ are the densities of brown plant-hopper in mature stage and in egg stage respectively in patch i, i ∈ {1, 2} m, m are the dispersal rates of brown planthopper in mature stage from region to region and opposite 77 Figure 4.9 Compare the density of rice on patch between the original model and the reduced one The case: rice wins globally in competition Parameters values are chosen as follows: r1 = 0.7; r2 = 0.2; K = 40; a1 = 0.2; a2 = 0.2; e1 = 0.05; e2 = 0.05; α1 = 0.2; α2 = 0.3; m = 0.3; m = 0.7; d1A = 0.3; d2A = 0.5; d1J = 0.2; d2J = 0.3 79 Figure 4.10 Equilibria and local stability analysis of the reduced model 87 Figure 4.11 The case: rice wins globally in competition Parameters values are chosen as follows: r1 = 0.7; r2 = 0.2; K = 40; a1 = 0.2; a2 = 0.2; e1 = 0.05; e2 = 0.05; α1 = 0.2; α2 = 0.3; m = 0.3; m = 0.7; d1A = 0.3; d2A = 0.5; d1J = 0.2; d2J = 0.3 92 Figure 4.12 The case: rice disappears on patch Parameters values are chosen as follows: r1 = 0.7; r2 = 0.2; K = 40; a1 = 0.5; a2 = 0.7; e1 = 0.6; e2 = 0.3; α1 = 0.1; α2 = 0.2; m = 0.3; m = 0.7; d1A = 0.3; d2A = 0.5; d1J = 0.2; d2J = 0.3 93 Figure 4.13 The case: the existence of rice and BPH on both patches Parameters values are chosen as follows: r1 = 0.3; r2 = 0.9; K = 40; a1 = 0.7; a2 = 0.1; e1 = 0.9; e2 = 0.5; α1 = 0.1; α2 = 0.1; m = 0.8; m = 0.2; d1A = 0.3; d2A = 0.5; d1J = 0.1; d2J = 0.3 94 INTRODUCTION Motivation The growth and degradation of populations in the nature and the struggle of one species to dominate other species have been an interesting topic for a long time The application of mathematical concepts to explain these phenomena have been documented for centuries The founders of mathematical-based modeling are Malthus (1798), Verhulse (1838), Pearl and Reed (1903), especially Lotka and Volterra whose most important results were published in the 1920s and 1930s Lotka and Volterra modeled, independently of each other, the competition between predator and prey Their work has important meaning for the population biology field They are the first to study the phenomenon of species interactions by introducing simplified conditions that lead to solvable problems that have meaning until today The proposed model is given by dN1 (t) N1 (t) N2 (t) − a12 = r1 N1 (t) − , dt K1 K1 N2 (t) N1 (t) dN2 (t) = r2 N2 (t) − − a21 , dt K2 K2 where Ni (t) is the density of species i, i ∈ {1, 2}, at time t Parameters ri and Ki are the growth rate and the carrying capacity of the species i, i ∈ {1, 2}, respectively Parameter aij is interspecific competitive coefficient representing the negative effect of species j on the growth of species i, i = j, i, j ∈ {1, 2} The ecological meaning of this model is that two species coexist only if the effects of their competition are small When the competing effects of two species are large, one of the two species will be extinct This famous principle is called the competitive exclusion principle Today, this model is still applicable to competitions between a number of biophysical species in practice and in empirical observations [22, 63] However, there are many other competing bio-systems, which cannot be explained by using the classic competition model of Lotka-Volterra (or the competitive exclusion principle) We present here two examples In the first example, Atkinson and Shorrocks [12] studied the competition of two species for having phytoplankton (food) in multiple environments Competition is noted when one of the two species is absent, resulting in an increase the remaining species Although the measured competing effect is significant, the two species coexist This result is contrary to the exclusion principle of the classical competition model In the second example, Lei and Hanski [63] studied two species of parasites on the same Melitaeacinxia butterfly The results showed that the more competitive and less hostile species (Cotesia melitaearum) were not founded in some host species, while the less competitive (Hyposotherapy horticola) were found in all hosts of Melitaeacinxia This result is also contrary to the principle of competitive exclusion The main reason for the limitation of Lotka-Volterra’s classic competition model is that there are too many assumptions in the model, such as the assumption that the environment is homogeneous and stable (expressed by the carrying capacities Ki for the specie i, i ∈ {1, 2}), the behavior of the individual species is the same and the competition is expressed only by interspecific competitive coefficient aij Meanwhile, these factors appear frequently and play a very important role For example, the migration behavior of individual species is a very important factor for species survival [80, 104] Individuals of the same species or of different species may have different behaviors Aggressive behavior is also used by individuals of wild species to compete for accommodation, to fight against their partners, etc In addition, individuals may also change their behaviors frequently according to the change of the environment as studied in [110, 111] Therefore, the development of new models that take into account the complex environments and the behaviors of individuals has been interested by many mathematicians Following are some recent approaches • The complex environment and individual migration behavior in competitive ecosystems The competition process and the migration process have the same time scale or different time scales • Aggressive behavior of individuals in competitive system • Age structure (adult group and immature group) in the competitive system Objective The objective of this thesis is to develop models for analyzing the effects of the environment, the behaviors of individuals (aggressive behavior, hunting habits) and the age structure (adults and juveniles) on the two species of competitive ecosystems To reach this goal, we divide this thesis into main work packages: - Developing models analyzing the effects of complex environments and aggressive behavior of the two competing ecosystems - Developing models analyzing the effect of age structure (adult and juvenile) the studied competing ecosystems - Building disk-graph based models to study competing ecosystems Figure 4.12: The case: rice disappears on patch Parameters values are chosen as follows: r1 = 0.7; r2 = 0.2; K = 40; a1 = 0.5; a2 = 0.7; e1 = 0.6; e2 = 0.3; α1 = 0.1; α2 = 0.2; m = 0.3; m = 0.7; d1A = 0.3; d2A = 0.5; d1J = 0.2; d2J = 0.3 Firstly, we re-write condition (4.19) of theorem as follows E1 x + E2 (1 − x) ≤ where E1 = a1 e1 K d1Jα+α − d1A and E2 = a2 e2 K d2Jα+α − d2A E1 (reps E2 ) are the functions of birth and death rates as well as the rate describing being matures of juveniles Therefore, E1 (reps E2 ) can be considered as an evolution function of BPH on patch (reps patch 2) There are several following situations: If E1 , E2 < then condition (4.19) holds In fact, in this case BPH does not grow on both patches Consequently, BPH globally gets extinct If E1 , E2 > then condition (4.19) does not hold In this case, BPH grows on both patches It is therefore survival If E1 < < E2 , then condition (4.19) is equivalent to > µ1 > E2 −E1 +E2 It means that if BPH is more likely to stay on patch (where it cannot grow), it eventually get extinct If E2 < < E1 , then condition (4.19) is equivalent to > µ2 > E1 −E2 +E1 Similarly, BPH will get extinct when it distributes more enough on patch Effects of age-structure parameters In order to study effects of age-structure parameters,we re-write condition (4.19) in the following way: G1 µ1 α1 + G2 µ2 α2 − (d1J d1A µ1 + d2J d2A µ2 ) ≤ where G1 = a1 e1 K − d1A and G2 = a2 e2 K − d2A G1 (reps G2 ) are the functions of birth and death rates That is the reason why we call G1 (reps G2 ) growth function of BPH in patch (reps patch 2) If G1 , G2 < then (∗) is satisfied This condition means 93 Figure 4.13: The case: the existence of rice and BPH on both patches Parameters values are chosen as follows: r1 = 0.3; r2 = 0.9; K = 40; a1 = 0.7; a2 = 0.1; e1 = 0.9; e2 = 0.5; α1 = 0.1; α2 = 0.1; m = 0.8; m = 0.2; d1A = 0.3; d2A = 0.5; d1J = 0.1; d2J = 0.3 that when BPH’s density decreases in both patches then BPH will die as a consequence If G1 , G2 > and (d1J d1A µ1 + d2J d2A µ2 ) ≈ then (4.19) is not satisfied This condition means that the energy that BPH gets from rice is higher than the loosing one due to death process BPH is therefore survival in both patches The present model takes into account a simple system of density-independence dispersal It would be interesting to consider density dependent dispersals in the model This would lead to a more complicated model and will be the subject of our future work Conclusion We have shown, in this chapter, some models for two ecological phenomena For the thiof-octopus system at the coast of Senegal, three models corresponding to three case with increasing complexity have proposed: (1) the case with refuge, (2) the case with refuge and density-independent migration and (3) the case with refuge and density-dependent migration For the rice-BPH system, a model is given by the equation-based approach We studied the local stability analysis and the global properties of their reduced models to get the knowledge about the complete model of these ecological phenomena 94 CONCLUSION Summary of contributions Competitive ecosystems have been under investigation for a long time Many models have been built to get the knowledge and to explain about these ecological phenomena in reality In this thesis, we have developed some continuous and discrete models for studying the effects of the environment, the local behaviors of individuals and the age structure of population on the competitive ecosystems both in theoretical and practical point of views The concrete results are given as follows: In term of theoretical point of view, Chapter dealt with the model with two opposite behaviors (aggressive and avoiding strategies) based on migration of individuals in a patchy (biotic and abiotic) environment revealed that under certain conditions, aggressiveness is efficient for survival of local inferior resource exploiter and even provokes global extinction of the local superior resource exploiter A new methodology of graph generating from individual-based models (a case study in predation dynamics) was proposed in Chapter A comparison with common graphs as well as the integration in term of biology point of view were reported In term of practical point of view, some effective models for two concrete ecological phenomena have been built in Chapter The competition rice-brown plant hopper model with stage structure of population showed some emerge results which support for decision makers for their management The competition thiof-octopus model coupling with fishing pressure figured out the strong increase of the fishing pressure in some areas leads to the depletion of the thiof and the invasion of its competitor, the octopus Futures works There are numerous potential research directions that we could investigate for improving the results in this thesis Here are some on which it would be nice to investigate In the current models, the behaviors of species were simple introduced by taking into account of only single competitive and non-competitive patches It would be also interesting to consider several competitive patches connected by migrations That could lead to much more complicated model but more interesting to investigate We considered, by using the discrete model, the competition of only two predatorprey species in a homogeneous environment As a perspective, we would like to consider more complex case studies of more than two species, for example, a system 95 of one prey and two predator species We would like also to consider more complex behaviors of individuals such as the migration behavior as well as take into account of the dynamics of the environment in the predator-prey system The present competition model about rice and brown plant hopper ecological system just took into account of the density-independent migration It would be interesting to consider density dependent migrations in the model for future study 96 BIBLIOGRAPHY [1] Aleksei, V (2004) Disk Graph, A Short Survey Springerlink, Oxford University Press [2] Alm, E and Arkin, A (2003) Biological networks Current Opinion in Structural Biology, 13: pp 193–202 [3] Amarasekare, P (2000) Spatial dynamics in a host-multi parasitoid community J Anim Ecol 69: pp 201–213 [4] Amarasekare, P (2000) Coexistence of competing parasitoids on a patchily distributed host: local vs spatial mechanisms Eco 81: pp 1286–1296 [5] Amarasekare, P (2002) Interference competition and species coexistence Proc R Soc Lond 269: pp 2541–2550 [6] Amarasekare, P., Nisbet R.M (2001) Spatial heterogeneity, source-sink dynamics, and the local coexistence of competing species The American Naturalist, 158: pp 572–584 [7] Amouroux, E., Chu, T., Boucher, A., and Drogoul, A (2007) Gama: an environment for implementing and running spatially explicit multi-agent simulations In PRIMA [8] Amouroux, E., Desvaux, S., and Drogoul, A (2008) Towards virtual epidemiology: An agent-based approach to the modeling of h5n1 propagation and persistence in north-vietnam In PRIMA08 [9] Amouroux, E., Gaudou, B., Desvaux, S., and Drogoul, A (2010) O.d.d.: a promising but incomplete formalism for individual-based model specification In RIVF2010, pp 123–126 [10] Arditi, R., Lobry, C., Sari, T (2015) Is dispersal always beneficial to carrying capacity? New insights from the multi-patch logistic equation? Theoretical Population Biology, 106: pp 45–59 [11] Arnold, V.I (1998) Ordinary Differential Equations Translated and Edited by Richard A Silverman The M.I.T Press 97 [12] Atkinson, W and Shorrocks, B (1981) Competition on a divided and ephemeral resource: a simulation model J Ani Eco., 50: pp 461–471 [13] Auger, P., Kooi, B.W., Bravo de la Parra, R., Poggiale, J.C (2006) Bifurcation analysis of a predatorprey model with predators using hawk and dove tactics J.Theo Biology, pp 597–607 [14] Auger, P., Bravo de la Parra R, Poggicale J C, Sanchez E, Nguyen Huu T (2008) Aggregation of Variables and Applications to Population Dynamics In Magal P, Ruan S (eds.), Structured Population Models in Biology and Epidemiology, Lecture Notes in Mathematics, Mathematical Biosciences Subseries, Springer, Berlin, 1936: pp 209–263 [15] Auger, P., Bravo de la Parra R, Poggiale JC, Sanchez E, Sanz L (2008) Aggregation methods in dynamical systems variables and applications in population and community dynamics Physics of Life Reviews, 5: pp 79–105 [16] Auger, P., Poggiale, J.C (1996) Emergence of population growth models: fast migration and slow growth Journal of Theoretical Biology, 182: pp 99–108 [17] Auger, P., Roussarie, R (1994) Complex ecological models with simple dynamics: From individuals to populations Acta Biotheoretica, 42: pp 1–11 [18] Auger, P., Bravo de la Parra, R., Morand, S., Sanchez, E (2002) A predatorprey model with predators using hawk and dove tactics Math Biosci 177, 178: pp 185–200 [19] Bae, S.H., and M.D Pathak (1970) Life history of Nilaparvatalugens (Homoptera: Delphacidae) and susceptibility of rice varieties to its attacks Ann.Entomol Soc Am., 63(1): pp 149-155 [20] Balgueras, E., Quintero, M.E and Hernndez-Gonzlez, C.L (2000) The origin of the Saharan Bank cephalopod fishery ICES J Mar Sci., 57(1): pp 15–23 [21] Barabasi, L and Oltvai, Z (2004) Network biology: understanding the cells functional organization Nature Reviews-Genetics, 5: pp 101–113 [22] Begon, M., Harper, J., and Townsend, C (1990) Ecology: Individuals, Populations and Communities Blackwell Scientific Publications, second edition [23] Begon, M., Colin R Townsend, John L Harper (2006) Ecology: From Individuals to Ecosystems Wiley-Blackwell Scientific Publications 98 [24] Bravo de Laguna (1982) Distribution and abundance of demersal resources of the CINECA region Rapports et Procs-verbaux des Runions du Conseil International pour lExploration de la Mer, 180: pp 432–446 [25] Caddy, J., and Rodhouse, P G (1998) Comparison of recent trends in cephalopod and groundfish landings: an indicator of widespread ecological change in global fisheries Reviews in Fish Biology and Fisheries, 8: pp 431–444 [26] Castella, F., Madec, S and Lagadeuc, Y (2016) Global behavior of N competing species with strong diffusion: diffusion leads to exclusion Applicable Analysis, Taylor and Francis, 95(2): pp 341–372 [27] Chen, Q., Chang, H., Govindan, R., Jamin, S., Shenker, S and Willinger, W (2002) The origin of power laws in internet topologies revisited In INFOCOM [28] Clark, B., Colbournm, C., and Johnson, D (1990b) Unit disk graphs Discrete Mathematics, 86: pp 165–177 [29] Cressman, R (1992) The Stability Concept of Evolutionary Game Theory Springer, Berlin [30] Christophe Crespelle, Matthieu Latapy and Thi Ha Duong Phan (2015) On the Termination of Some Biclique Operators on Multipartite Graphs Discrete Applied Mathematics, 195: pp 59–73 [31] Dale, M.R.T and Fortin, M.-J.(2010) From graph to spatial graph Annual Review of Ecology, Evolution and Systematics, 41: pp 21-38 [32] Demarcq, H (2009) Trends in primary production, sea surface temperature and wind in upwelling systems (1998-2007) Progress in Oceanography, 83(1-4): pp 376-385 [33] Diallo, M., Jouffre, D., Caverivire, A., Thiam, M (2012) The demographic explosion of Octopus vulgaris in Senegal during the 1999 summer Bulletin of Marine Science, 71(2): pp 1063–1065 [34] Diatta, Y., Bou-Ain, A., Clotilde-Ba, F and Capap, C (2003) Diet of four serranid species from the Senegalese coast (eastern tropical Atlantic) Acta Adriatica 44(2): pp 175–182 [35] Domain, F., Jouffre, D., Caverivire, A (2000) Growth of Octopus vulgaris from tagging in Senegalese waters J Mar Biol Ass U.K., 80(4): pp 699–706 99 [36] Drogoul, A (1993) When ants play chess (or can strategies emerge from tactical behaviors?) In Actes de MAAMAW93 Lausanne [37] Drogoul, A., Ferger, J., and Jacopin, E (1991) Pengi: Applying eco-problemsolving for behavior modelling in an abstract ecosystem In Modelling and Simulation: Proceedings of ESM91, Simulation Councils, Copenhague, pp 337–342 [38] Efendiev M.A., Eberl H.J (2007) On positivity of solutions of semi-linear convection-diffusion-reaction systems, with applications in ecology and environmental engineering RIMS Kyoto Kokyuroko, 1542: pp 92–101 [39] Eubank, S., Guclu, H., Kumar, V S A., Marathe, M V., Srinivasan, A., Toroczkai, Z., and Wang, N (2004) Modelling disease outbreaks in realistic urban social networks Nature, 429: 180–184 [40] Fabian, K., Roger, W., and Aaron, Z (2008) Ad-hoc networks beyond unit disk graphs Wireless Networks, 14: pp.715–729 [41] Ferrer, R and Sole, R.V.(2001) The small-world of human language In Proceedings of the Royal Society of London, volume B268, pp 2261–2265 [42] Garca-Ramos, G., Kirkpatrick, M (1997) Genetic models of adaptation and gene flow in peripheral populations Evolution, 51: pp 21–28 [43] Gaudou, B., Ho, T., and Marilleau, N (2009) Introduce collaboration in methodologies of modeling and simulation of complex systems In International Conference on Intelligent Networking and Collaborative Systems, Barcelone, Spain, pp 1–8 [44] Jean-Loup Guillaume and Matthieu Latapy (2004) Bipartite Graphs as Models of Complex Networks Combinatorial and Algorithmic Aspects of Networking, pp 127–139 [45] Jeong, H., Oltvai, Z., and Barabasi, A (2003) Prediction of protein essentiality based on genomic data Journal of Theoretical Biology, 1: pp 19–28 [46] Griebeler, E and Seitz, A (2002) An individial-based model for the conservation of the endangered large blue butterfly, maculineaarion (lepidoptera: Lycaenidea) Ecological Modeling, 156: pp 43–60 [47] Grimm, V., Berger, U., Bastiansen, F., Eliassen, S., Ginot, V., Giske, J., GossCustard, J., Grand, T., S.K Heinz, G H., Huth, A., Jepsen, J., Mooij, W., Piou, C., Railsback, S., Robbins, A., Robbins, M., Rossmanith, E., N, N R., 100 Strand, E., Souissi, S., Stillman, R., Visser, U., and DeAngelis, D (2006) A standard protocol for describing individual-based and agent-based models Ecological Modelling 198(1-2): pp 115–126 [48] Grimm, V and Railsback, S (2005) Individual-based Modeling and Ecology Princeton University Press [49] Hanski, I (1999) Metapopulation Ecology Oxford University Press, Oxford [50] Ho, H.S and Liu, T.H (1969) Ecological investigation on brown planthopper in Taichung district in Chinese, English Summary Plant Prot.Bull (Taiwan), 11(1): pp 33–42 [51] Hofbauer, J., Sigmund K (1998) Evolutionary Games and Population Dynamics Cambridge University, Cambridge [52] Hoppe, A.H (1975) The brown plant-hopper Rat Fighter, 7: pp 38–39 [53] Huson, M.L and Sen, A (1995) Broadcast scheduling algorithms for radio networks Military Communications Conference IEEE MILCOM, pp 647–651 [54] Huston, M., DeAngelis, D., and Post, W (1988) New computer models unify ecological theory BioScience, 38: pp 682–691 [55] Jeong, H., Tombor, B., Albert, R., Oltvai, Z and Barabllasi, A.-L (2000) The large-scale organization of metabolic networks Nature, 407: pp 651–654 [56] Jouffre, D (1998) Octopus vulgaris as a component of the benthic fauna of the NW African coast: A note on an investigation of species community organisation using multifactorial analysis S Afr J mar Sci: 20: pp 93–100 [57] Khajanchi, S.(2014) Dynamic behavior of a Beddington-DeAngelis type stage structured predatorprey model Applied Mathematics and Computation, 244: pp 344–360 [58] Khalil, H.K.(2002) Nonlinear Systems, 3rd ed., Upper Saddle River, New Jersey: Prentice Hall [59] Kleinberg, J.M., Kumar, R., Raghavan, P., Rajagopalan, S and Tomkins, A.S.(1999) The Web as a graph: Measurements, models, and methods In T Asano, H Imai, D T Lee, S Nakano, and T Tokuyama, editors, Proc 5th Annual Int Conf Computing and Combinatorics, COCOON, Springer-Verlag, number 1627, pp 1–17 101 [60] Kuznetsov, Y.A (1998) Elements of applied bifurcation theory 2nd Edition, Springer [61] LaSalle, J.P (1976) The stability of dynamical systems Regional Conference Series in Applied Matheamtics, SIAM (Philadelphia) [62] Lam, Pham Van (2013) The small insects and spiders damage plants discovered in Vietnam Agricultural Publisher: pp 266–267 [63] Lei, G and Hanski, I (1998) Spatial dynamics of two competing specialist parasitoids in a host metapopulation J Ani Eco., 67: pp 422–433 [64] Levin, S (1992) The problem of pattern and scale in ecology Ecology, 73: pp 1943–1967 [65] Levin, S.A., Pimentel D (1981) Selection of intermediate rates of increase in parasitehost systems American Naturalist, 117: pp 308–315 [66] Lott, D.F (1991) Intraspecific Variation in the Social Systems of Wild Vertebrates Cambridge University, New York [67] Liu, X and Zhao, X.-Q (2011) A Periodic Epidemic Model with Age Structure in A Patchy, Environment SIAM J Appl Math., 71: pp 1896–1917 [68] Lotka, A.J (1932) The growth of mixed populations: two species competing for a common food supply Journal of the Washington Academy of Sciences, 22: pp 461–469 [69] Meissa, B., Gascuel, D and Rivot, E (2013) Assessing stocks in data-poor African fisheries: a case study on the white grouper Epinephelus aeneus of Mauritania African Journal of Marine Science, 35(2): pp 253–267 [70] MacArthur, R.H., Levins, R (1964) Competition, habitat selection, and character displacement in a patchy environment Proceedings of the National Academy of Sciences of the USA, 51: pp 1207–1210 [71] Marva, M., Moussaoui, A., Bravo de la Parra, R., Auger, P (2013) A density dependent model describing age-structured population dynamics using hawk-dove tactics Journal of difference equations and applications, 19(6): pp 1022–1034 [72] Misra, B.C., Israel, P (1970) The leaf and plant-hopper problems of high yielding varieties of rice oryza 7: pp 127–130 102 [73] Moon, J and Moser, L (1965) On cliques in graphs Israel Journal of Math, pp 23–28 [74] Moussaoui, A., Nguyen-Ngoc, D., Auger, P (2015) Analysis of a model describing stage-structured population dynamics using hawk-dove tactics ARIMA, 20: pp 127–143 [75] Moussaoui, A., Nguyen-Ngoc, D., Auger, P (2012) A stage structured population model with mature individuals using hawk and dove tactics 11th African Conference on Research in Computer Science and Applied Mathematics, CARI, 2012, pp 251–258 [76] Mouquet, N., Loreau, M (2002) Coexistence in meta communities: the regional similarity hypothesis Am Nat 159: pp 420–426 [77] Morozov A.Y., Arashkevich E.G (2010) Toward a correct description of zoo plankton feeding in models: taking into account food-mediated unsynchronized vertical migration Journal of Theoretical Biology, 262: pp 346–360 [78] Morozov A.Y., Arashkevich E.G., Nikishina A., Solovyev, K (2011) Nutrientrich plankton communities stabilized via predatorprey interactions: revisiting the role of vertical heterogeneity Mathematical Medicine and Biology, 28(2): pp 185–215 [79] Murrell, D and Law, R (2000) Beetles in fragmented wood lands: a formal framework for dynamics of movement in ecological landscapes Journal of Animal Ecology, 69(3): pp 471–483 [80] Nee, S., May, R.M (1992) Dynamics of metapopulations: habitat destruction and competitive coexistence J Anim Ecol 61: pp 37–40 [81] Newman, M.E.J (2001) Scientific collaboration networks: I Network construction and fundamental results Phys Rev E, 64 [82] Newman, M.E.J., Watts, D.J and Strogatz, S.H (2002) Random graph models of social networks Proc Natl Acad Sci USA, 99 (Suppl.1): pp 2566–2572 [83] Nguyen-Ngoc Doanh (2010) PhD thesis: Coupling Equation-Based and Individual-Based Models in the Study of Complex Systems: A Case Study in Theoretical Population Ecology [84] Nguyen-Ngoc Doanh, Phan Thi Ha Duong, Nguyen Thi Ngoc Anh, Alexis Drogoul, Jean-Daniel Zucker (2010) Disk Graph-Based Model: a graph theoretical 103 approach for linking agent-based models and dynamical systems International Conference on Computing and Communication Technologies Research, Innovation and Vision for the Future, IEEE, DOI: 10.1109/RIVF.2010.5633033 [85] Nguyen-Ngoc D., Nguyen-Huu T., Auger P (2012) Effects of refuges and density dependent dispersal on interspecific competition dynamics International Journal of Bifurcation and Chaos, 22(2): pp 125-129 [86] Nguyen-Ngoc D., Nguyen-Huu, T., Auger, P (2012) Effects of fast density dependent dispersal on pre-emptive competition dynamics Ecol Complex, 10: pp 26–33 [87] Nguyen-Ngoc D., Bravo de la Parra, R., Zavala, M., Auger, P (2010) Competition and species coexistence in a metapopulation model: Can fast asymmetric migration reverse the outcome of competition in a homogeneous environment Journal of Theoretical Biology, 266(2): pp 256–263 [88] Ngoan, N.D (1971) Recent progress in rice insect research in Viet Nam In proceedings of a symposium on rice insects Tokyo, Japan Trop Agric Res Ser.5, Tropical Agriculture Rereach Center, Tokyo, pp 133–141 [89] Ndiaye, W., Thiaw, M., Diouf, K., Ndiaye, P., Thiaw, O.T and Panfili, J (2013) Changes in population structure of the white grouper Epinephelus aeneus as a result of long-term overexploitation in Senegalese waters African Journal of Marine Science 35.4: pp 465–472 [90] Pachpatte, B.G (1998) Inequalities for differential and integral equations San Diego: Academic Press ISBN 9780080534640 [91] Pachpatte, B.G (1996) Comparison theorems related to a certain inequality used in the theory of differential equations Soochow Journal of Mathematics, 22(3): pp 383-394 [92] Palla, G., Dernyi, I., Farkas, I and Vicsek, T (2005) Uncovering the overlapping community structure of complex networks in nature and society Nature, 435: pp 814–818 [93] Polis, G.A., Myers, C.A., Holt, R.D (1989) The ecology and evolution of intraguild predation: potential competitors that eat each other Annual Review of Ecological Systems, 20: pp 297–330 [94] Railsback, S F (2001) Concepts from complex adaptative systems as a framework for individual-based modeling Ecological Modeling, 139: pp 47–62 104 [95] Rathjen, W F., et Voss, G L (1987) The cephalopod fisheries: a review In: Cephalopod Life Cycles, Comparative Reviews, P.R Boyle (Ed), Academic Press (London), Vol II: pp 253–275 [96] Salihu, A.(2012) New method to calculate determinants of n×n (n ≥ 3) matrix, by reducing determinants to 2nd order International Journal of Algebra, 6(9): pp 913–917 [97] Samanta, M and Liang, S (2003) Predicting protein functions from redundancies in large-scale protein interaction networks In Proceedings of the National Academy of Sciences, 100: pp 12579–12583 [98] Stephen K Brown, Peter J Auster, Liz Lauck, and Michael Coyne (1998) National Oceanic and Atmospheric Administration (NOAA) (on-line).”Ecological Effects of Fishing” NOAA’s State of the Coast Report, Silver Spring [99] Steven H.Strogatz (2007) Nonlinear Dynamics And Chaos Sarat Book House [100] Taillandier, P and Buard, E (2009) Designing agent behaviour in agent-based simulation through participatory method In The 12th International Conference on Principles and Practices in Multi-Agent Systems (PRIMA), pp 571–578 [101] Thiao, D., Chaboud, C., Samba, A., Lalo, F., Cury, P (2012) Economic dimension of the collapse of the false cod Epinephelus aeneus in a context of ineffective management of the small-scale fisheries in Senegal Afr J Mar Sci., 34: pp 305–311 [102] Thiaw, M., Gascuel, M., Thiao, D., Thiaw, O.T and Jouffre, D (2011) Analyzing environmental and fishing effects on a short-lived species stock: the dynamics of octopus population, Octopus vulgaris, in the Senegalese waters African Journal of Marine Science, N33-2, pp 209–222 [103] Tilman, D (1982) Resource competition and community structure, Princeton University Press [104] Tilman, D., May, R M., Lehman, C L., and Nowak., M A (1994) Habitat destruction and the extinction debt Nature, 371: pp 65–66 [105] Vance, R.R (1984) Interference competition and the coexistence of two competitors on a single limiting resource Ecology, 65: pp 1349–1357 [106] Verhulst, F (2007) Singular perturbation methods for slow-fast dynamics Nonlinear Dynamics, 50(4): pp 747–753 105 [107] Vlachos, C., Paton, R., Saunders, J R., and Wu, Q H (2006) A rule-based approach to the modeling of bacterial ecosystems Biosystems, 84: pp 49–72 [108] Watts, D.J and Strogatz, S.H (1998) Collective dynamics of small-world networks Nature, 393: pp 440–442 [109] Wang, B.G., Li, W.T and Zhang, L (2016) An almost periodic epidemic model with age structure in a patchy environment Discrete and continuous dynamical systems, 21(1): pp 291–311 [110] Wolf, L.L., Waltz, E (1993) Alternative mating tactics in male white- faced dragonflies: Experimental evidence for a behavioural assessment ess Anim Behav 46: pp 325–334 [111] Yamane, A., Doi, T., Ono, Y (1996) Mating behaviors, courtship rank and mating success of male feral cat (Felis catus) J Ethol 14: pp 35–44 106 LIST OF PUBLICATIONS Thuy Nguyen-Phuong, Doanh Nguyen-Ngoc (2016) Effects of Behavioural Strategy on the Exploitative Competition Dynamics Acta Biotheoretica, 64, pp 495-517 (SCI) Thuy Nguyen-Phuong, Doanh Nguyen-Ngoc, Pierre Auger, Sidy Ly, Didier Jouffre (2016) Can Fishing Pressure Invert the Outcome of Interspecific Competition? The Case of the Thiof and of the Octopus Along the Senegalese Coast Acta Biotheoretica, 64, pp 519-536 (SCI) Thuy Nguyen-Phuong, Oanh Tran-Thi-Kim, Doanh Nguyen-Ngoc, Effects of Fast Dispersal and Stage-Structured on Predator-Prey Dynamics: A Case Study of Brown Plant-Hopper Ecological System (accepted for publication in Vietnam Journal of Mathematical Application) Thuy Nguyen-Phuong, Doanh Nguyen-Ngoc, Duong Phan-Thi-Ha, On the Generating Graph of an Individual-Based Predator-Prey Model (submitted) 107 ... Figure 4.1 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.7; r2 = 1.3; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3;... Figure 4.2 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.9; r2 = 0.7; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3;... Figure 4.4 Example of the case where the inferior competitor wins globally in model Parameters are chosen as follows r1 = 0.9; r2 = 1.3; K1 = 100; K2 = 100; a12 = 0.8; a21 = 2.2; q1 = 0.6; q2 = 0.3;