Abridged English version
Today, the theory of the chemostat has been widely studied in the literature (see, for instance, the textbook [14]). The classical model describing the competition of several species on a single resource s is the following one:
(1) |
In this work, we propose to replace in (1) the classical functional responses , that are only resource dependent, by growth functions that depend on both the available resource s and the quantity of consumers . Then, the ith equation of the model (1) becomes:
(2) |
We introduce the concept of steady-state characteristic for each species i, which maps for each fixed s the steady state of the single equation (2) for initial condition . Then, we show that all the essential information for determining the asymptotic behavior of the species in competition with other species on the resource s relies on the knowledge of the characteristics .
Our main result is the following: assume that the function is bounded and negative for s larger than a certain number . Let be a solution of the equation
This persistence result generalizes previous results where the 's are of ‘Contois’ type [6], or where is linear [9]. Finally, considering functional responses that depend on both the resource and the consumer is also a way of generalizing ratio-dependent models, such as those introduced by Arditi and Ginsburg in prey–predator interaction models [2].
1 Introduction
Le chémostat a été introduit il y a plus d'un demi-siècle dans l'étude de la croissance des micro-organismes, et le modèle mathématique ci-dessous :
(3) |
Un modèle analogue de compétition entre plusieurs organismes dans le même dispositif a été ensuite étudié, soit :
(4) |
D'autres travaux ont considéré le modèle plus général :
(5) |
(6) |
(7) |
En 1976, Armstrong et McGehee [3,12] ont fait remarquer que la coexistence ne se fait pas nécessairement à l'équilibre et que des oscillations forcées (modèle du chémostat (4) où ou peuvent être des fonctions périodiques du temps) ou endogènes (cas du système (5)) peuvent expliquer la coexistence. Ces travaux sont à l'origine de nombreux résultats sur les modèles (4) et (5), dont certains, mathématiquement sophistiqués, ont contribué à élucider une partie des mécanismes qui assurent la coexistence des espèces (toujours dans le contexte des modèles de compétition pour une ressource).
En 1989, Arditi et Ginzburg [2], dans le contexte de l'étude de la relation « proie–prédateur » ont proposé de remplacer la réponse fonctionnelle par une fonction « densité-dépendante », c'est-à-dire de la forme . Cette approche a été très fructueuse pour l'amélioration de la compréhension de la relation proie–prédateur et des chaînes trophiques, mais ne semble pas avoir été explorée dans le cadre de l'étude de la compétition. C'est dans cet esprit que nous proposons de modifier le modèle de croissance d'une population de façon très simple, en introduisant une dépendance par rapport à la densité de la population dans la réponse fonctionnelle des modèles (4) et (5). Nous démontrons ainsi un résultat simple assurant, sous des hypothèses assez larges, l'existence d'un équilibre localement stable pour lequel il y a coexistence des espèces.
2 Le modèle « un consommateur – une ressource »
La quantité de ressource disponible à l'instant t est désignée par et la quantité de consommateur par . Le taux de croissance de la population est une fonction continue et différentiable de s et de x, que nous désignons par . Nous supposons également que le taux de disparition (la disparition pouvant être la conséquence de la mortalité, de migrations, de prélèvements, etc.) du consommateur est constant et égal à d. Ainsi le modèle de croissance de la population du consommateur est donné par :
(8) |
Hypothèses H1
Sous l'
Hypothèse H1
, les propriétés suivantes sont vérifiées.
Proposition 1
Soit fixé. D'après nos hypothèses, la fonction est strictement décroissante, et strictement négative pour x suffisamment grand. Ainsi, soit s'annule pour une valeur positive ou nulle, que l'on note , soit et l'on pose . Le point (i) de la Proposition 1 est alors une application immédiate de la théorie élémentaire des équations différentielles ordinaires [4]. Soit défini par (et par +∞ si d est plus grand que pour tout s). D'après nos hypothèses de croissance sur , la fonction est identiquement nulle pour s plus petit que , et devient strictement positive et croissante après. Le point (ii) de la proposition est ainsi démontré. □Preuve
Nous appelons caractéristique à l'équilibre la fonction , définie par la Proposition 1 ci-dessus (cf. Fig. 1). Exemple de caractéristique à l'équilibre.Définition 1
La caractéristique à l'équilibre que nous venons de définir est une partie de ce qu'on appelle usuellement l' « isocline du prédateur » dans la relation « proie–prédateur » du modèle (9) ci-dessousRemarque 1
Nous avons préféré lui donner ici un autre nom car, d'une part il ne s'agit pas d'une isocline au sens strict du terme et, d'autre part, nous allons en faire un usage différent de celui classiquement fait dans l'étude des systèmes dynamiques en dimension deux.(9)
En théorie mathématique des systèmes, la quantité s s'appellerait « entrée » du système décrit par l'Éq. (8). Dans un article récent [1], l'importance des systèmes pour lesquels à chaque entrée on peut associer un unique équilibre a été démontrée. Cette correspondance a été appelée input-to-state characteristic. Bien que notre fonction ne soit qu'un cas particulier de cette notion, nous n'avons pas conservé cette appellation qui nous paraît obscure dans le contexte de la dynamique des populations. De même, nous aurions pu nous référer à certaines des méthodes générales mises au point par ces auteurs, mais nous avons préféré garder une approche mathématique élémentaire. Cela dit, nous devons reconnaître ici que la théorie mathématique des systèmes développée pour les besoins de l'automatique est notre principale source d'inspiration [11], et que les méthodes proposées dans [1] ont également des applications en écologie [7].Remarque 2
Exemple 1
On reconnaît dans la fonction le modèle classique de croissance de Monod. Cette fonction de croissance est multipliée par une fonction qui réduit la croissance d'autant plus que la taille de la population est grande. Un calcul élémentaire montre que la fonction caractéristique de croissance est, dans ce cas :(10)
où si , et sinon.(11)
Exemple 2
Dans ce cas, la fonction est identiquement nulle lorsque . Sinon,(12)
où . Le modèle « ratio dépendant » est le cas limite obtenu pour .(13)
Exemple 3
Dans ce cas, les Hypothèses H1 ne sont pas satisfaites (puisque la fonction est constante par rapport à la variable x et la fonction n'est pas définie). Revenons à l'Éq. (8) et constatons que trois cas sont possibles, lorsque .(14)
avec défini, comme dans l'Exemple 1, par l'équation . Ce cas peut être vu comme le cas limite de l'Exemple 1 lorsque la fonction possède une décroissance vers 0 de moins en moins forte, c'est-à-dire lorsque α tend vers 0.
Considérons le système(9), dont nous supposons que la réponse fonctionnellesatisfait lesHypothèses H1 (donc dont la caractéristique à l'équilibre est bien définie), et supposons de plus que la fonction f est bornée, strictement positive pour et strictement négative pour . Si est une solution de l'équationThéorème 1
telle que et , alors avec est un équilibre localement (exponentiellement) stable du système (9).(15)
Elle est extrêmement simple. Nous devons d'abord montrer que est une solution du système de deux équations :Preuve
(16)
L'Éq. (17) est satisfaite par définition de la fonction . Si nous remplaçons dans (16) par d et par , nous obtenons , ce qui est précisemment l'hypothèse faite dans le théorème. La stabilité se démontre en calculant la matrice jacobienne du système (9) au point (voir Annexe A) et en vérifiant que l'équation caractéristique de cette dernière possède bien deux racines à parties réelles strictement négatives, ce qui, dans le cas d'une matrice , est immédiat. Nous ne poursuivons pas, car le Théorème 1 est un cas particulier du Théorème 2 énoncé ci-dessous. □(17)
Le Théorème 1 ne présente aucune originalité et n'est là que pour rappeler ce que nous généralisons. Il est d'autre part connu que, sans hypothèse supplémentaire, cet équilibre n'a aucune raison d'être globalement stable.
3 Le modèle de compétition « n consommateurs – une ressource »
Nous considérons dans ce paragraphe le modèle :
(18) |
Considérons le système(18)dont nous supposons que tous lessatisfont lesHypothèses H1 (donc dont les caractéristiques à l'équilibre sont bien définies). Supposons de plus que la fonction est bornée, strictement positive pour et strictement négative pour . Si est une solution de l'équationThéorème 2
telle que et alors est un équilibre localement (exponentiellement) stable du système (18).(19)
Nous devons montrer quePreuve
est une solution du système de équations(20)
Par définition des , les n Éqs. (21) sont toutes satisfaites quand on remplace par . Remplaçons alors dans (20) les par . Il vient :(21)
qui est précisément l'équation dont l'hypothèse du théorème suppose qu'elle possède une solution . Il reste à montrer que cet équilibre est stable. Pour cela, il suffit de calculer la matrice jacobienne et de montrer que ses valeurs propres sont à parties réelles négatives. Ceci exige quelques calculs élémentaires sur les nombres complexes, que nous repoussons en Annexe A. □(22)
Le Théorème 1 est clairement le cas particulier du Théorème 2 dans le cas où .
Il est clair que pour tous les i tels que , les à l'équilibre correspondant sont strictement positifs, alors que pour tous les autres indices i, les correspondant sont nuls et les espèces correspondantes sont éliminées. On en déduit que, les caractéristiques des consommateurs étant données par les et les , on pourra toujours trouver un assez grand pour que toutes les espèces soient présentes à l'équilibre. Cette proposition peut s'interpréter inversement : pour un donné et un entier n arbitraire, on peut définir des et des tels que toutes les espèces soient présentes à l'équilibre.
4 Discussion
4.1 Réalisme biologique des Hypothèses H1
Pour de nombreuses d'espèces, il est réaliste d'imaginer que le taux de croissance d'une population, toutes choses égales par ailleurs, diminue quand la taille de la population augmente. C'est un phénomène de compétition « intraspécifique » pour l'accès à la ressource. Toutefois, il peut exister, pour certaines espèces, une forme de mutualisme qui contredirait cette hypothèse. Il existe un cas, celui de la croissance d'une espèce dans un milieu liquide, où la diffusion du substrat peut devenir un facteur limitant si les organismes ont tendance à constituer des colonies plus ou moins sphériques. Dans ce cas, seule la zone en périphérie a accès à la ressource et alors la considération du rapport volume/surface conduit à introduire une réduction du taux de croissance tel que dans l'Exemple 2. Le cas classique de la « loi d'action de masse » pure où la fonction ne dépend pas de x (Exemple 3) apparaît comme un cas limite dégénéré de notre modèle. C'est une dégénérescence forte, car la caractéristique à l'équilibre cesse d'exister en tant que fonction. C'est pourquoi les conclusions de notre Théorème 2 ne sont plus vraies dans le modèle classique « proie dépendant ».
Le cas de la « ratio dépendance », dès lors qu'on introduit un petit paramètre α comme dans l'Exemple 3, apparaît comme un cas particulier de notre modèle. L'introduction de ce petit paramètre ne change rien d'essentiel à la philosophie développée par les partisans de la ratio-dépendance. Une partie des arguments en faveur de la ratio dépendance peut donc être reprise en faveur de notre modèle.
Nous avons supposé un taux de disparition d constant. Poursuivant la même logique il serait judicieux de le rendre lui aussi dépendant de la densité et d'écrire
(23) |
Il est possible de le faire sans changer les conclusions de cette note, mais la présentation aurait été un petit peu plus compliquée. Toujours par souci de simplicité, nous avons supposé que la réponse fonctionnelle est égale au taux de croissance. On sait que, si ces derniers sont proportionnels, on peut, par simple changement d'unités, éliminer la constante de proportionnalité. Mais, même dans le cas où il n'y a pas proportionnalité, on peut sans peine étendre les résultats de la note.
Dans la Réf. [10] sont développés, dans le contexte de l'écologie microbienne, des arguments en faveur de réponses fonctionnelles du type de notre Exemple 3.
4.2 Critique du modèle de compétition proposé
Le modèle défini par les Éqs. (18) n'est pas la conséquence logique du modèle défini par les Éqs. (9), comme les Éqs. (4) le sont de (3). En effet, dans ce dernier cas, il était supposé que tous les individus d'une espèce ont à tout moment le même accès à la ressource, et il y a une certaine logique à supposer qu'il en est de même pour un mélange d'individus d'espèces différentes. Dans notre cas, dans la mesure où nous supposons une certaine forme de compétition entre les individus d'une même espèce pour accéder à la ressource, il serait peu logique de ne pas supposer qu'elle existe aussi entre les individus des espèces différentes. Cela se traduirait mathématiquement par l'introduction de fonctions dépendant de toutes les variables :
(24) |
Toutefois, nous pouvons justifier notre modèle avec l'argument suivant. Dans tous les cas où les individus de la même espèce ont tendance à rester ensemble (que ce soit la conséquence d'un comportement grégaire ou la conséquence de la division cellulaire dans le cas de micro-organismes), les effets de la compétition intraspécifique seront plus importants que ceux de la compétition interspécifique. Mathématiquement, grâce aux théorèmes généraux de stabilité structurelle [5], nous pouvons affirmer que, si les dérivées partielles de par rapport aux variables (j différent de i) sont assez petites en valeur absolue devant la valeur absolue de la dérivée partielle par rapport à , les conclusions du Théorème 2 restent identiques. En d'autres termes, si la compétition interspécifique est faible devant la compétition intraspécifique, la coexistence des espèces est possible. On retrouve le fait élémentaire bien connu dans le cas des équations de Lotka Volterra. Cependant, il reste souhaitable que des résultats mathématiques plus poussés viennent quantifier ce résultat.
4.3 Critique de l'étude à l'équilibre
Les remarques de Armstrong et McGhee restent valables, mais il faudrait continuer à examiner le modèle (18), non plus du seul point de vue des équilibres, mais aussi de celui de la persistance « tout court ». Dans la note [9], il est montré que, dans le cas du chémostat ( linéaire et égaux), seule la persistance à l'équilibre est possible. Mais les travaux récents de [7] montrent le même phénomène dans un cas ou la mortalité des espèces dépend de la densité de la population des espèces.
Dans [6], une condition suffisante de persistance est démontrée pour le modèle particulier où la fonction est une logistique et où les sont des fonctions de type « Contois ». Dans la note [13], un résultat du même type est proposé dans le cadre plus général que nous proposons ici.
Annexe A
Soit un système de n équations différentielles
(A.1) |
(A.2) |
La matrice jacobienne du système (A.1) à l'équilibre est la matrice définie par :Définition 2
(A.3)
Soit λ une valeur propre de la matrice J, c'est-à-dire un nombre réel ou complexe tel qu'il existe un vecteur V non nul vérifiant . Rappelons que si pour toute valeur propre λ de J, la partie réelle de λ est strictement négative, alors l'équilibre est exponentiellement stable (cf. [4]). Dans notre problème la matrice jacobienne (de taille ) prend la forme suivante :
(A.4) |
(A.5) |
(A.6) |
(A.7) |
- – soit , et comme un au moins des doit être non nul, on doit avoir pour cet indice i. La valeur propre correspondante à un tel vecteur propre est donc réelle négative ;
- – soit , et alors nous pouvons calculer en fonction de et reporter dans (A.6), ce qui donne, après avoir fait passer une partie du premier membre dans le second :
(A.8)
(A.9) |
(A.10) |
(A.11) |
(A.12) |