L'encyclopédie des Sciences
  Dynamique des Populations
 

Einstein, à la fin de sa vie, avait coutume d'évoquer les trois explosions qui allaient sous peu menacer l'humanité: l'explosion des bombes nucléaires, l'explosion de nos savoirs, l'explosion de l'effectif des hommes. Le rôle des scientifiques à ce niveau n'est pas seulement d'améliorer leur savoir, mais de la partager et surtout de diffuser la conscience qu'ils ont acquise des conséquences qu'ils savent discerner. Einstein avait à coeur de jouer ce rôle. Au-delà des péripéties immédiates, il savait voir les évolutions à long terme, Par ces évolutions, celle du nombre des hommes lui apparaissait à l'évidence comme grosse des plus grands dangers. Aujourd'hui, en début de ces années 2000, nous pouvons constater sa perspicacité. L'humanité est réellement prise à la gorge par l'accroissement de son effectif ou de son maintient actuel sans la forcer à la baisse.

Pour introduire les mathématiques sociales, il nous faut d'abord déterminer les caractéristiques qui décrivent la dynamique du nombre d'individus avant de formaliser les propriétés uniques qui les caractérisent.

TABLES DE MORTALITÉ ET NATALITÉ (FONCTION BIOMÉTRIQUES)

Nul ne sait ni le jour ni l'heure de sa mort. Cette évidence individuelle n'est plus pertinente si nous nous intéressons non à telle personne, mais à une collectivité assez nombreuse. Alors jouent les compensations entre ceux qui succombent à des accidents prématurées et ceux qui échappent quasi miraculeusement aux pires dangers. Nous pouvons alors décrire la façon dont est payé globalement le tribut à la mort en considérant un grand nombre d'enfants nés la même année et en précisant, grâce aux données de l'état civil, comment leur effectif diminue peu à peu pour un jour s'annuler.

Un tel ensemble de conscrits est appelé par les démographes une "cohorte". Considérons donc la cohorte des Français mâles nés en 1900. Chaque année nous pouvons, en regroupant les indications de l'état civil, calculer le nombre de ceux qui sont décédés dans la cohorte.

Représentons par  le quotient du nombre des décès entre les anniversaires a et a+1 par l'effectif initial de la cohorte.

Remarque: Nous allons dans ce qui suit définir les bases des probabilités nécessaires au calcul actuariel. En effet, les mathématiques actuarielles réunissent le calcul financier et le calcul des probabilités. Le paiement d'un capital n'est plus certain et dépend par exemple de la survie d'une personne.

La suite des ces nombres contient la totalité de l'information nécessaire pour étudier la mortalité de cette cohorte.

Nous pouvons en déduire la proportion d'entre eux survivant à l'âge a:

  (1)

Nous pouvons également caractériser l'intensité de la mortalité à chaque âge en divisant le nombre des décès entre les âges a et a+1 par le nombre de survivants. Ce nombre est alors appelée le "quotient de mortalité":

  (2)

La liste âge par âge de ces trois paramètres d, S, et q est la "table de mortalité" de la cohorte étudiée. Le tableau ci-dessous en donne le résumé pour les âges multiples de 5.

Âge
a

Survivants

Décès de l'âge a à a+5

Probabilité des décès

0

1

0.228

0.228

5

0.772

0.13

0.017

10

0.759

0.9

0.012

15

0.750

0.23

0.031

20

0.727

0.23

0.032

25

0.704

0.20

0.028

30

0.684

0.21

0.031

35

0.663

0.25

0.038

40

0.638

0.29

0.045

45

0.609

0.26

0.043

50

0.583

0.36

0.062

55

0.547

0.47

0.086

60

0.500

0.62

0.124

65

0.438

0.78

0.178

70

0.360

   
  (3)

Choisissons au hasard dans ce tableau le nom d'un enfant inconnu dans la liste des naissances de l'année 1900. La grande question lors de cette naissance était: combien d'années vivra-t-il? Aujourd'hui, nous sommes en mesure de répondre rétroactivement à cette question, tout au moins en évoquant des probabilités, car nous connaissons la table de mortalité de cette cohorte.

Si la seule information dont nous disposons aujourd'hui à propose de cet enfant est le fait qu'il est né en 1900, nous pouvons déclarer que la probabilité qu'il ait été encore vivant à l'âge de 5 ans est égale à 0.772, à l'âge de 50 ans de 0.583… Autrement dit que la probabilité pour qu'il soit mort avant 5 ans est égale à 0.228.

Si nous désignons aujourd'hui un individu inconnu sur la liste de ceux qui ont été incorporés au cours de l'année 1920 à l'âge de 20 ans, nous pouvons de même calculer les probabilités des diverses durées de sa vie, mais nous avons une information supplémentaire: il était encore vivant à 20 ans, il a évité les risques de mort avant cet âge. La probabilité qu'il atteigne alors l'âge de 50 ans est devenue:

  (4)

Donc la probabilité pour une personne d'âge a d'être en vie à l'âge a+1 est égal à:

  (5)

La probabilité pour une personne d'âge a de décéder entre l'âge a et a+n est logiquement donne par:

  (6)

Ainsi à chaque âge nous pouvons donner la loi de la variable "durée encore à vivre". Cette loi peut être résumée en indiquant son espérance (cf. chapitre de Statistiques). Un calcul immédiat permet d'en donner la valeur à chaque âge en fonction de la table de mortalité.

Ce résultat peut être passionnant pour un historien, mais il donner une réponse à une question posée il y a longtemps et depuis oubliée. Ce qui nous intéresse est le présent. Cet enfant qui vient de naître, quelle est son espérance de vie? Pour répondre, il faudrait connaître la table de mortalité de sa génération, or le don de prémonition n'existe pas. L'attitude probabiliste permet de contourner cette difficulté à condition de bien préciser les hypothèses sous-jacentes.

Ainsi, pour répondre à la question: quelle est l'espérance de vie des nouveaux nés de 1990, nous faisons alors l'hypothèse, tout à fait gratuite, qu'ils rencontreront à chaque âge, à l'avenir, les conditions qu'ont rencontrées en 1990 les individus de ces âges: en l'an 2000 ils subiront la même mortalité que celle subie en 1990 par ceux qui sont nés en 1980. Bien sûr, personne n'imagine que la réalité sera conforme à cette hypothèse, mais le calcul qu'elle permet fournit une image synthétique des conditions actuelles de la lutte contre la mort.

C'est cette hypothèse qui fait que l'espérance de vie rajeunit parfois les vieux (leur espérance de vie ayant tendance à augmenter au fur et à mesure qu'ils deviennent plus âgés grâce aux progrès de la science…).

Prenons pour calcul l'espérance de vie la table de mortalité pour les hommes en Suisse en 1983-1993:

Âge
n

Survivants

Espérance de vie
E(a)

Âge
n

Survivants

Espérance de vie
E
(a)

0

1

73.68823

55

0.90224

22.81093

1

0.99246

73.24806

56

0.89581

21.97466

2

0.99183

72.29459

57

0.88875

21.14922

3

0.99148

71.32011

58

0.88099

20.33551

4

0.99117

70.34241

59

0.87247

19.53409

5

0.9909

69.36158

60

0.86312

18.7457

6

0.99066

68.37838

61

0.85288

17.97077

7

0.99044

67.39357

62

0.84169

17.20969

8

0.99022

66.40855

63

0.82948

16.46301

9

0.99001

65.42263

64

0.81618

15.73128

10

0.98971

64.44246

65

0.80174

15.01462

11

0.9896

63.44963

66

0.78609

14.31354

12

0.98938

62.46373

67

0.76918

13.62821

13

0.98915

61.47826

68

0.75096

12.95887

14

0.98889

60.49442

69

0.73138

12.30579

15

0.9886

59.51217

70

0.7104

11.66921

16

0.98823

58.53445

71

0.68798

11.04949

17

0.9877

57.56586

72

0.66409

10.44699

18

0.98692

56.61136

73

0.63873

9.861773

19

0.98581

55.6751

74

0.6119

9.294182

20

0.98439

54.75541

75

0.58362

8.744543

21

0.98285

53.84121

76

0.55393

8.21324

22

0.98131

52.9257

77

0.52291

7.700465

23

0.97975

52.00997

78

0.49066

7.206599

24

0.97815

51.09505

79

0.45733

6.731813

25

0.97653

50.17981

80

0.42309

6.276608

26

0.9749

49.26371

81

0.38819

5.840903

27

0.97328

48.34571

82

0.35288

5.425357

28

0.97166

47.42631

83

0.31748

5.030301

29

0.97007

46.50405

84

0.28234

4.656372

30

0.9685

45.57943

85

0.24785

4.304337

31

0.96694

44.65297

86

0.21446

3.974494

32

0.96541

43.72373

87

0.18263

3.667196

33

0.96388

42.79314

88

0.15281

3.382828

34

0.96236

41.86073

89

0.12546

3.120277

35

0.96082

40.92782

90

0.10091

2.879397

36

0.95927

39.99395

91

0.07942

2.658524

37

0.95768

39.06035

92

0.06111

2.455081

38

0.95604

38.12736

93

0.04593

2.266492

39

0.95434

37.19528

94

0.03373

2.086273

40

0.95257

36.26439

95

0.02414

1.915079

41

0.95071

35.33534

96

0.0168

1.751786

42

0.94874

34.40871

97

0.01134

1.595238

43

0.94665

33.48468

98

0.00739

1.447903

44

0.94441

32.5641

99

0.00464

1.306034

45

0.94201

31.64706

100

0.00279

1.172043

46

0.93942

30.73431

101

0.0016

1.04375

47

0.93662

29.82619

102

0.00087

0.91954

48

0.93356

28.92396

103

0.00045

0.777778

49

0.93023

28.0275

104

0.00021

0.666667

50

0.92659

27.1376

105

0.00009

0.555556

51

0.92259

26.25526

106

0.00004

0.25

52

0.91821

25.3805

107

0.00001

0

53

0.91338

24.51471

108

0

 

54

0.90808

23.65779

     
  (7)

Donc la courbe représentative, appelée "l'ordre des vivants" est:


  
(8)

Voyons comment calculer l'espérance de vie. Pour cela considérons un homme en vie lors de son a-ème anniversaire. Le nombre d'années qui lui reste à vivre est une variable aléatoire dont nous pouvons calculer l'espérance mathématique (cf. chapitre de Probabilités). En négligeant les fractions d'années, cette espérance peut s'écrire:

  (9)

Si a est pris comme étant égal à zéro, les démographes parlent de EDVN (espérance de vie à la naissance).

Voici pour information l'augmentation de l'espérance de vie chez les hommes (source: Institut National Français d'Études Démographiques) depuis 1996 à 2006:

1996

74.1

1997

74.5

1998

74.8

1999

75.0

2000

75.3

2001

75.5

2002

75.7

2003

75.9

2004

76.7

2005

76.8

2006

77.2

  (10)

Nous pouvons donc observer que l'espérance de vie augment d'un peu plus d'un an tous les quatre ans depuis plus de 50 ans (et jusqu'à quand…)


Source: Wikipédia   
(11)

Pour faire un enfant, il faut être deux certes les démographes le savent. Mais cette double source de chacun des nouveaux nés leur pose de tels problèmes de description et d'analyse de la fécondité qu'ils préfèrent en général l'ignorer. Leur attitude est justifiée par le fait que seule la conception nécessite l'intervention de deux acteurs. Au moment de la naissance, la mère agit seule. Or la démographie ne s'intéresse pas aux conceptions, inaccessibles à l'observation, mais seulement aux naissances.

Nous avons vu comment nous pouvons suivre un ensemble d'hommes ou de femmes nés une année donnée et enregistrer les décès successifs, ce qui permet d'établir la table de mortalité de cette cohorte. De même, nous pouvons noter les effectifs des enfants auxquels ils donnent naissance année après année. Nous obtenons ainsi la "table de fécondité".

Il suffit de suivre une cohorte de femmes de 15 jusqu'à 50 ans pour avoir une description complète de son comportement procréateur. Les données fournies par l'état civil permettent de calculer chaque année le nombre des enfants auxquels ont donné naissance les femmes de cette cohorte, regroupées par âge ou par groupes d'âge. En divisant par le taux des femmes survivants à cet âges, nous obtenons le "taux de fécondité" .

Si nous additionnons l'ensemble de ces taux, nous obtenons le nombre d'enfant qu'auraient eu, en moyenne, les femmes de cette cohorte si leur mortalité avait été nulle. Tel n'a évidemment pas été le cas. Pour caractériser la façon dont elles ont assuré le renouvellement de leur génération, il faut additionner les nombres réels moyens de naissances, produits du taux  par le taux de survie .

L'ensemble de ces données est présentés dans une table de fécondité dont voici deux exemples où  est le nombre d'effectif de naissances chez les femmes survivants à cette tranche d'âge et où  est le nombre de naissances de 1'000 femmes de cette tranche d'âge.

Femmes françaises nées vers 1830:

Âge
n

Taux de survie

Effectifs naissances

Taux de fécondité

15

0.672

91

135

20

0.645

464

720

25

0.616

589

955

30

0.587

475

810

35

0.558

328

565

40

0.528

153

290

Total

1

2'100

3'475

  (12)

Par exemple, sur 1'000 petites filles nées en 1830, 645 on atteint l'âge de 20 ans et ont eu 464 enfants entre 20 et 25 ans. L'intensité de la fécondité est mesurée par le nombre de naissances qu'auraient eu 1000 femmes de cet âge:

  (13)

Et la moyenne des naissances:

Âge
n

Survivantes

Effectifs naissances

15

672

91

20

645

464

25

616

589

30

587

475

35

558

328

40

528

153

Total

 

2'100

  (14)

Nous pouvons observer qu'en 1830, de 20 à 30 ans les femmes donnaient naissance en moyenne à:

  (15)

ce chiffre étant assimilé par le grand public aux "taux de fécondité" (donc attention à ne pas confondre avec ).

Pour faire le calcul sur l'ensemble de tous les âges il suffit de prendre le rapport du nombre d'enfants sur le nombre initial de femmes. Cela donne mondialement des chiffres indiqués dans la carte ci-dessous:


Source: Wikipédia   
(16)

RENOUVELLEMENT DE LA POPULATION

La question essentielle pour un ensemble humain en renouvellement permanent en raison des flux d'entrées et de sorties que sont les naissances et les décès est: notre effectif est-il en décroissance ou en croissance?

La table de natalité féminine permet de réponde, grâce au rapport du nombre de garçons à la naissance sur le nombre de filles. Ainsi, en ce début du 21ème siècle: il naît en moyenne 105 garçons pour 100 filles dans le monde. Ainsi, la proportion des filles et donc de:

  (17)

Vers quarante ans (49 en France selon l'INED), la prépondérance s'inverse et le nombre de femmes l'emporte généralement sur le nombre d'hommes, malgré de notables disparités régionales.


Source: Wikipédia   
(18)

La table de fécondité de 1830 montre donc qu'en moyenne une femme de la cohorte de 1830 a produit:

fille. Soit une augmentation que nous noterons k (en analogie avec le modèle exponentiel que nous verrons plus loin) de 3%. A tort, certains politiciens avancent la valeur (scolaire) 2.1 comme étant le taux de fécondité… qu'il faut pour assurer le renouvellement des générations ce qui n'est donc pas tout à fait exact.

L'effectif féminin était donc en accroissement. Donc la population pouvait assurer son renouvellement (tant le rapport est supérieur à 1).

Le nombre ainsi obtenu est le "taux net de reproduction". Ce taux est donc normalement constitué par le rapport entre le nombre de filles mises au monde par cent femmes, rapport corrigé par la mortalité prévue entre la naissance de ces filles et l'âge moyen à la reproduction, car une partie des filles n'atteindra pas l'âge de la reproduction, étant donné les décès survenus parmi elles entre leur naissance et leur âge à la maternité. L'âge moyen à la reproduction est donné par l'âge moyen des mères à la naissance.

MODÈLES DE POPULATIONS

Il existe de nombreux modèles mathématiques permettant d'étudier la croissance d'une population. Le terme "population" est utilisé ici au sens le plus large - il peut s'agir d'une population d'humains, d'animaux, de plantes, de personnes infectées par un virus etc.

Pour construire un modèle mathématique, il est nécessaire de faire des hypothèses. Ces hypothèses jouent deux rôles : préserver certaines caractéristiques essentielles de la réalité et simplifier suffisamment cette réalité afin qu'elle puisse être étudiée par la mathématique.

MODÈLE EXPONENTIEL

Dans ce modèle de dynamique des populations (un des plus simples), l'hypothèse sera la suivante : le taux de variation de la population est proportionnel, en tout temps t, à la population P(t) présente au temps t.

Nous pouvons penser, à priori, que cette hypothèse est raisonnable pour une foule de situations. Par exemple plus la population humaine est grande et plus le taux de variation de cette population, exprimé en nombre de personnes qui s'ajoutent par unité de temps, sera grand. De même, plus il y a de personnes infectées par un virus et plus, dans les semaines qui viennent, il y aura de nouveaux cas de personnes infectées.

Mathématiquement, cette hypothèse peut se traduire à l'aide de l'équation différentielle (cf. chapitre de Calcul Différentiel Et Intégral) :

  (19)

Cette équation différentielle est un modèle mathématique représentant une situation où le taux de croissance de la population est proportionnel à la grandeur de la population en tout temps . Dans ce cas, k est une constante appelée "taux d'accroissement" et nous verrons plus loin comment nous peuvons la déterminer. Dans certaines situations, la valeur de k est négative indiquant le fait que la population diminue avec le temps au lieu de croître. Il est évident qu'une solution à cette équation différentielle (cf. chapitre de Calcul Différentiel Et Intégral) est :

  (20)

En premier lieu nous déterminerons la valeur de la constante , à partir de données démographiques pour l'année 1965. À cette époque, il y avait 3 milliards de personnes sur la planète Terre. De plus, à cette époque, la population augmentait de 54 millions par année.

Ainsi, en 1965 :

  (21)

Ce qui nous donne :

  (22)

Soit 1.8%. L'équation différentielle est alors :

  (23)

Pour déterminer la constant multiplicative, il suffit de poser et de la choisir en conséquence (puisqu'elle correspond à la condition initiale). Ainsi, en 1965 nous avions :

  (24)

Si ce modèle mathématique est conforme à la réalité, la solution trouvée nous permettra d'estimer la population d'humains sur la terre pour des temps ultérieurs à 1965. Voici le graphique de la fonction :


  
(25)

Si nous évaluons P(37), cela nous fournira la prédiction pour la population en 2002. Nous trouvons 5.84 milliards ce qui est assez proche de la réalité.

Le lecteur pourra également facilement vérifier le tableau suivant selon les valeurs de k:

Taux d'accroissement par an
k

Temps de doublement en années

0.5%

139

1%

69

1.5%

46

2%

35

2.5%

28

3%

23

3.5%

20

4%

17

  (26)

Nous voyons alors immédiatement, que selon le modèle exponentiel, que une croissance qui peut paraître lente, de l'ordre de 3% par an, est en fait véritablement explosive puisqu'elle entraîne un doublement tous les 23 ans. Soit une multiplication par plus de 17 en un siècle. Or, un siècle est vite passé!

Cette simple constatation montre à quel point dans ce modèle sont mensongers les arguments qui nous font chercher la solution d'un problème économique ou social dans la croissance, et même, comme il est souvent admis, dans une croissance durable. Il est clair qu'aucune croissance ne peut véritablement durer, elle n'est qu'un épisode transitoire, nécessairement suivi d'un palier ou même d'une décroissance. Résoudre une difficulté par la croissance ou le maintient de celle-ci c'est reporter le problème à plus tard, à la période où il faudra trouver les moyens à la fois d'arrêter cette croissance et de résoudre autrement cette difficulté restée un temps camouflée.

Le modèle de croissance que nous avons présenté ne peut par contre être valable sur de très longues périodes de temps. En effet, si on calculait, en utilisant l'équation précédente, la population dans 7 siècles le résultat serait que sur chaque mètre carré de la terre, excluant l'eau, il y aurait, en moyenne, dix humains! De même, une population de personnes infectées par un virus ne peut pas vraiment être décrite par un tel modèle.

Ces derniers résultats nous enseignent que si on veut élaborer une modélisation de la croissance d'une population, qui soit davantage conforme à la réalité, il va falloir modifier nos hypothèses initiales. Ce que nous allons de suite voir.

MODÈLE LOGISTIQUE déterministe

Nous allons maintenant nous intéresser à un autre type de modèle autre qu'exponentiel (où la population explose) qui à l'avant d'avoir d'avoir un comportement asymptotique plutôt que divergent.

Ce type de comportement est intéressant car les ressources sont normalement limitées et qu'il y a une compétition entre individus. Le modèle logistique, également appelé "modèle de Verlhust" permet de rendre compte de cela relativement bien.

Soit N(t) la population au temps t. Posons :

  (27)

r est le taux d'accroissement qui cette fois-ci ne sera pas constant sera défini comme valant :

  (28)

K est la capacité maximale du milieu. Nous voyons que si K est infini que nous retombons immédiatement sur le modèle exponentiel et que si N(t) égale K alors r est nul.

Finalement, nous avons :

  (29)

Soit :

  (30)

Soit sous forme mathématique :

  (31)

avec comme condition initiale que .

Il s'agit maintenant de résoudre cette équation différentielle (cf. chapitre de Calcul Différentiel Et Intégral).

Donc nous avons à résoudre l'équation différentielle :

  (32)

Nous posons . L'équation différentielle devient alors :

  (33)

ce qui se simplifie en :

  (34)

Cette équation différentielle se résout comme à d'habitude. Nous écrivons d'abord l'équation homogène :

  (35)

Une solution particulière immédiate est :

  (36)

A est une constante. En injectant cette solution dans l'équation différentielle initiale donne :

  (37)

Nous voyons donc immédiatement que pour l'égalité soit satisfaite notre solution particulière devient la solution générale si nous l'écrivons sous la forme :

    (38)

La solution de l’équation différentielle du début est donc,

  (39)

où nous avons remplacé AK par A qui dépend bien évidemment des conditions initiales.

Donc dans notre cas d'étude, la solution s'écrira :

  (40)

qui est donc la relation finale du modèle logistique déterministe. Nous voyons d'ailleurs aisément que si t tend vers l'infini alors l'asymptote horizontale est K.


  
(41)

moDÈLE LOGISTIQUE chaotique

Nous allons supposer que chaque génération est proportionnelle à la précédente. D'une période à la suivante, l'évolution de la population peut alors se traduire par une suite du type:

  (42)

 représente la population à la période t,  la population à la période  et k le taux de reproduction.

Mais on se rend très rapidement compte qu'un tel modèle est irréaliste: la terre aurait depuis longtemps été submergée par une marée humaine. En effet, le modèle exponentiel ignore des facteurs importants, comme celui des ressources qui ne sont pas illimitées

Pour rétablir un modèle plus satisfaisant, il faut tenir compte du fait que les ressources disponibles sont limitées et que, pour tout territoire, il existe une population maximale au-delà de laquelle la population décroît, quelle que soit l'espèce. Pour trouver une fonction qui traduirait de façon plus réaliste l'évolution d'une population, énumérons, en les simplifiant un peu, les propriétés que doit vérifier cette évolution:

1. C'est un phénomène itératif: si  représente l'effectif de la population de la période , il dépend de celui de la période précédente .

2. Une population ne peut croître indéfiniment sur un territoire délimité: il existe un maximum après lequel elle décroît. Il faut donc prévoir un facteur rétroactif limitant la hausse de population quand sa densité devient trop élevée. Pour simplifier les calculs, x ne représentera pas la population en nombre absolu, mais en pourcentage de ce maximum correspondant à un territoire donné; x ne peut donc fluctuer qu'entre 0 et 1.

3. k représente le taux de croissance effectif d'une période à la suivante.

L'effectif de la population de la période  sera la valeur , exprimée en pourcentage de la population maximale que peut accueillir le territoire donné, et sera obtenu de l'effectif de la population  de la période précédente t par l'équation:

  (43)

le facteur  représentant l'effet rétroactif. Cette suite est souvent appelée "suite logistique".

En effet, quand la densité de la population est élevée, proche de la saturation, alors  est proche de 1 et, par conséquent,  sera proche de 0. Donc, ce facteur rétroactif aura tendance à minimiser la population .

Nous allons voir que selon la valeur de son taux de croissance effectif k, une population animale, peut tendre vers un état d'équilibre, ou fluctuer entre deux ou quatre ou huit valeur, ou varier de façon totalement aléatoire.

Ainsi, pour certaines espèces animales, il est normal que les populations varient régulièrement, tandis que pour d'autres, il est normal de tendre vers une situation d'équilibre. Cette variété, ce "chaos" dans certains cas, est liées aux propriétés mathématiques même de l'équation:

  (44)

La complexité devient la règle et non l'exception, et ce qui semble "chaotique" découle des propriétés bien précises d'une fonction bien précise.

Ce sont ces résultats que l'on va retrouver. Pour cela, on va étudier l'évolution de populations animales dont les taux de croissance effectifs k varient de 1 à 4. On étudiera ces évaluations par calcul et par l'étude du graphe de la fonction logistique.

Voyons quelques cas particuliers:

- Si la valeur initiale  vaut 0.1 et  l'évolution de la population pour les quarante premières périodes est représentée par la courbe ci-dessous:


  
(45)

On constate que pour ce taux de croissance et cette condition initiale, la population décroît et va tendre vers zéro.

- Si la valeur initiale  vaut 0.1 et  l'évolution de la population pour les quarante premières périodes est représentée par la courbe ci-dessous:


  
(46)

On constate que pour ce taux de croissance, la population va se stabiliser autour d'un nombre correspondant à la moitié de la population que le territoire pourrait accueillir.

- Si la valeur initiale  vaut 0.1 et  l'évolution de la population pour les quarante premières périodes est représentée par la courbe ci-dessous:


  
(47)

Avec un taux de croissance effectif de 3, le nombre d'individus de cette population se met à osciller entre deux valeurs que l'on peut calculer à l'aide de l'ordinateur. On obtient alors les valeurs 0.64 et 0.68.

- Si la valeur initiale  vaut 0.1 et  l'évolution de la population pour les quarante premières périodes est représentée par la courbe ci-dessous:


  
(48)

Avec un taux de croissance effectif de 3.5 la population oscille entre quatre valeurs: 0.39, puis 0.83, puis 049, et enfin 0.87. L'évolution d'une population qui a un tel taux de croissance effectif est nettement cyclique.

- Si la valeur initiale  vaut 0.1 et  l'évolution de la population pour les quarante premières périodes est représentée par la courbe ci-dessous:


  
(49)

On constate que, lorsque le taux de croissance effectif k est égal à 4, le nombre d'individus des populations des périodes successives semble osciller irrégulièrement, de façon "chaotique" entre les deux extrêmes: la saturation quand x tend vers 1 et l'extinction quand x tend vers 0.

La question maintenant, c'est qu'elle est le comportement de cette fonction pour d'autres valeurs initiales de population.

Pour le vérifier, on reprend les calculs précédents avec une population initiale représentant 0.8 de l population maximale pour un territoire donné et comparons les résultats obtenus.

La croissance de la population est toujours donnée par l'équation:

  (50)

- Selon que la population initiale soit  ou , on obtient pour :


  
(51)

On constate que, quelle que soit la population initiale, elle tend vers zéro.

- Selon que la population initiale soit  ou , on obtient pour :


  
(52)

On constate que, quelle que soit la population initiale, elle tend vers zéro.

-  Selon que la population initiale soit  ou , on obtient pour :


  
(53)

On constate ici que l'évolution des populations est très différente. On a pris des populations initiales éloignées. Si on avait pris des populations initiales très proche, quelle aurait été l'évolution des deux populations ?

-  Selon que la population initiale soit  ou , on obtient pour :


  
(54)

On constate ici que l'évolution des populations est très différente même si elles étaient très proches initialement et même si l'évolution de ces population d'écoulent d'une fonction très simple. Si on avait pris comme population initiale  on aurait constaté également une évolution très différente de l'évolution des populations.

Conclusion: en dehors de la phase chaotique, la valeur initiale n'a aucune importance, mais dans cette phase chaotique, au contraire, la plus petite variation de valeur initiale change du tout au tout les valeurs suivantes. C'est ce qu'on appelle "l'effet papillon".

DIAGRAMME DE FEIGENBAUM

Pour comprendre l'évolution d'une population selon le modèle logistique, on a représenté l'évolution dans le temps, quarante périodes, d'une population correspondant à une valeur initiale précise et à une constante k déterminée. On a constaté, en prenant quelques valeurs particulières de k, que pour ces différentes valeurs, la population avait une évolution différente. On va donc étudier le comportement de la fonction logistique en prenant k comme variable.

En donnant à k des valeurs comprises entre 0 et 4 avec un pas de 0.02, nous allons, pour chacune des ces valeurs, calculer ce que deviendra la population lapopulation pour chaque période comprise entre la 30ème et la 130ème .

Le graphique ci-dessous représente donc pour chaque valeur du taux de croissance effectif k en abscisse, cent valeurs successives de  en ordonnées pour t variant de 30 à 130.


  
(55)

Ce graphique (où l'on peut observer des bifurcations), est appelé "diagramme de Feigenbaum", du nom du physicien Mitchell Feigenbaum qui l'a étudié en profondeur et a montré qu'on le retrouvait dans de nombreux phénomènes naturels.

On va maintenant examiner deux propriétés remarquables de ce diagramme: le doublement de période et sa dimension fractale.

Si on revient à notre exemple où la fonction logistique donne l'évolution d'une population animale selon son taux de croissance effectif k, le diagramme de Feigenbaum indique que, quand ce taux est inférieur à 3, le système tend vers un état final stable. Cela correspond en général à notre intuition influencée par notre désir inconscient d'ordre et de simplicité.

Mais à partir de 3, cela se complique: le nombre d'individus par génération se met à osciller entre 2 puis 4, puis 8 valeurs… pour enfin entrer dans une zone "chaotique" où toute valeur semble possible. Il y a déjà quelque chose de fascinant. Mais si on observe de plus près cet intervalle de doublement qui précède la zone chaotique, on va découvrir des résultats encore plus curieux.

On prend ainsi d'abord l'intervalle des trois premiers doublements avec 2, puis 4, puis 8 branches, pour un taux de croissance effectif de  et on obtient le graphique suivante:


  
(56)

Il est remarquable d'observer également sur le graphique ci-dessus qu'il y a deux points qui sont totalement stables !

Si on prend maintenant l'intervalle  on obtient le diagramme de Feigenbaum suivant:


  
(57)

On constate sur le graphique ci-dessus que les bifurcations se multiplient à partir de points de plus en plus rapprochés et sur des intervalles de plus en plus courts. Le physicien Feigenbaum a démontré deux résultats curieux:

- les bifurcations vont se multiplier à l'infini (!) sur un intervalle qui ne dépassera pas le point de d'abscisse 3.5699456 appelé le "point de Feigenbaum" ou "porte d'entrée sur le chaos", car, après ce point, le système devient chaotique. Il se met à fluctuer entre des valeurs "imprévisibles" et devient extrêmement sensible aux conditions initiales.

- la longueur des intervalles propres aux différentes classes de bifurcation (2, 4, 8, 16,…) diminue dans un rapport constant, 4.6692… appelé, bien sûr, "constante de Feigenbaum".

Peitgen indique que ces bifurcations, points et constante de Feigenbaum, ne se retrouvent pas seulement dans le cas de la fonction logistique étudiée par May, mais dans de nombreux phénomènes physiques comme l'hydrodynamique, l'électronique, les lasers ou l'acoustique.

Remarquons enfin un dernier phénomène bien curieux, celui des "fenêtres". Si on examine plus attentivement la zone chaotique comprise entre le point de Feigenbaum et , on constate qu'il existe des zones étroites où la fonction logistique se remet à osciller entre un nombre fini de bifurcations avant de replonger dans le chaos.


  
(58)

On constate que, dans cette fenêtre, il y des bifurcations qui ressembleraient au diagramme de Feigenbaum lui-même si on l'agrandissait.

On aurait pu reprendre les calculs entre  et on aurait alors constaté que on retrouve la même figure et que, si on l'agrandissait à l'infini, on retrouvait toujours (!) la même figure toujours et toujours:

             
  
(59)

On vient ainsi de découvrir une figure fractale: l'attracteur étrange de la fonction logistique.

Grâce à la théorie mathématique du chaos appliquée à la dynamique des populations, l'écologie reçoit une stimulations décisive. Jusqu'au début des années 1960, le débat sur la dynamique des populations opposait les tenants d'une théorie déterministe, voyant des évolutions régulières des populations ne subissant qu'exceptionnellement de brutales variations, à ceux pensant ces évolutions comme purement aléatoires. Mais de nombreux faits restaient mal expliqués. En particulier les explosions cycliques de certaines populations et leur périodicité étrange n'entraient dans aucune des deux explications. En montrant que des modèles déterministes peuvent donner naissance à un  comportement aléatoire, Robert May, réconcilie ces points de vue à partir d'une théorie plus profonde. Ce qui apparaît à niveau d'appréhension comme une instabilité généralisée peut se concevoir à un autre comme un chaos stable. Une situation mathématique chaotique peut se révéler stable au point de vue écologique.

Révolution aux conséquence profondes dans la théorie de l'écologie, et aux implications pratiques non moins importante, l'approche de Robert May est un trait de lumière dans la situation théorique passablement confuse qui règne encore dans la dynamique des populations.

Si, à cause de la montée en puissance des ordinateurs, la théorie du chaos a permis des avancées importantes dans différentes disciplines, elle soulève un débat majeur, celui du déterminisme. Dans quelle mesure la science permet-elle de prédire l'avenir ? Pour certains, les résultats obtenus jusqu'ici dans le cadre de la théorie du chaos prouvent l'importance aux conditions initiales. Pour ces personnes, les équations déterministes n'ont qu'une portée limitée et l'avenir restera imprévisible. Pour d'autres au contraire, les résultats obtenus mont que l'on peut trouver un ordre et des lois dans ce qui peut sembler chaotique. Ces lois sont tout simplement plus complexes.

C'est un débat ouvert et vif qui déborde largement les milieux scientifiques.

Pour les utilisateur de Maple voici le code qui permet de générer le diagramme de Feigenbaum:

with(plots): with(plottools)
feigenbaum:=proc(début,fin,pas) local k,itéré,a,b,s;
s:={}; a:=début;
while a<=fin do itéré:=0.1;
  for k to 50 do itéré:=a*itéré*(1-itéré) od;
    for k to 100 do itéré:=a*itéré*(1-itéré);
      s:=s union {[a,evalf(itéré,4)]};
    od;
    a:=a+pas
od;
plot([op(s)],'a'=début..fin,style=POINT,symbol=POINT)
end:
feigenbaum(1,4,0.01);

LOI DE MALTHUS

Nous avons vu précédemment quelques modèles déterministes et chaotiques utilisant des taux d'accroissements pour la simulation. Introduisons maintenant un autre type de modèle utilisant le numéraire des enfants, des femmes et de leur fécondité respective.

Pour ce modèle, nous énonçons les hypothèses suivantes :

H1. Au départ N individus dont N/2 sont des femelles.

H2. Le taux d'accroissement r est supposé constant à cause d'un taux de fécondité  f constant aussi.

H3. La population s'accroît

Nous avons donc comme données la population au temps t notée N(t), la population féminine notée  et  le nombre d'enfants.

Avec les relations suivantes :

  (60)

et donc :

  (61)

donc le taux d'accroissement est :

  (62)

Finalement nous obtenons la loi de Malthus :

  (63)

ce modèle est donc continu, prend en compte la fécondité mais diverge…

MODÈLE DE LESLIE

Le modèle de Leslie est un peu plus perfectionné que les autres modèles déterministes mais tout aussi empirique (il est possible comme dans tout domaine de la science de construire des modèles théoriques valides et toujours plus complet et complexes).

Outre le taux de fécondité et de mortalité, il permet de prendre en compte les tranches d'âges de la population et certaines de leurs propriétés relativement aux deux facteurs précités. L'inconvénient de ce modèle est cependant les trop nombreux paramètres à déterminer pour l'ensemble des classes d'âges…

Le système se base sur la découpe de la population en tranches d'âges tels que par exemple :

-  : est le nombre d'individus de 1 an

-  : est le nombre d'individus de 2 ans

-  : est le nombre d'individus de plus 10 ans

Ensuite, l'idée est que l'évolution d'une classe d'âge dépende des autres classes d'âges. Par exemple, les naissances sont données par le taux de reproduction r sommé sur toutes les classes d'âges (bien évidemment pour certaines d'entre elles le taux est nul…) tel que :

  (64)

Il est de tradition et raisonnable pour l'être humain d'admettre que normalement seulement les classes d'âges pour lesquelles  et . Soit :

  (65)

Ensuite, le vieillissement et la mortalité m seront prises en compte par les relations :

  (66)

Il est relativement aisé de voir que ces équations peuvent être mises sous forme matricielle (cf. chapitre d'Algèbre Linéaire) de la manière suivante :

  (67)

où la matrice contenant les coefficients de mortalité et de reproduction est appelée "matrice de Leslie".

De manière plus compacte cela s'écrit :

  (68)

ou en partant de la population initiale :

N(t)=AN(0)   (69)

Il est possible de faire des analyses très intéressantes sur ce modèle relativement à l'âge de faire des enfants et aux conséquences y relatives. Ce modèle est relativement beaucoup utilisé en biologie marine.

PROPAGATION DES ÉPIDÉMIES

La "théorie des épidémies" fournit de nombreux systèmes intuitifs d'équations différentielles ou aux dérivées partielles en temps discret dans lesquelles interviennent des phénomènes de contamination, de diffusion, etc.

Nous allons ici prendre l'exemple de la diffusion de la rage dans une population de renards, et en présenter un modèle assez simple où nous considérons que les individus (les renards) peuvent se trouver dans trois états différents :

1. Il y a des renards sains S

2. Des renards infectés I

3. Et des renards morts R

Cette hypothèse des classification en trois catégories fait que ce modèle de propagation des épidémies est souvent appelé "modèle SIR".

Pour les plus optimistes, nous pouvons aussi utiliser le modèle qui suppose que certains individus sont immunisés et ne peuvent dons plus être infectés. Finalement, à trois états différents, nous pouvons lier trois équations d'évolutions propres telles que :

  (70)

avec :

  (71)

où les S, I, R sont donnés en % de la population totale.

Dans ce modèle, la première équation correspond au phénomène de contamination : lorsque cohabitent des individus sains et des individus infectés, un certain nombre d'individus sains sont infectés. Il est par ailleurs naturel de considérer ce terme comme proportionnel au produit . En effet, la quantité de microbes dans le milieu défini (donc la facteur r pour un individu sain donné d'être infecté en présence d'un individu infecté) est proportionnelle à I selon le facteur. Il nous faut ensuite multiplier cette probabilité par le nombre d'individus sains, c'est-à-dire par S. Le signe négatif est quant à lui présent pour signifier la diminution de la population saine.

Pour ce qui est des variations de I, le premier terme correspond aux individus contaminés (qui augmentent I). Le deuxième terme, en -aI, correspond aux individus qui meurrent.

Le tracé ci-dessous, montre un calcul de l'évolution des populations S, I, R en fonction du temps (une seconde étant considérée comme une boucle de l'algorithme). Selon les valeurs des paramètres a, r les comportements sont complètement différents : si a est trop élevé par exemple, les individus infectés meurent presque tout de suite et non pas le temps d'en infecter beaucoup d'autres (…). Voici un tracé de l'évolution des population au cours du temps, effectué avec l'algorithme suivant sur Matlab avec


  
(72)

L'algorithme étant du type suivant :

Pour i=1 à (temps max)/(par du temps) faire
  dS=-r*S(i)*I(i);
  dI=r*S(i)*I(i)-a*I(i);
  dR=a*I(i);
  S(i+1)=S(i)+dt*dS;
  I(i+1)=I(i)+dt*dI;
  R(i+1)=R(i)+dt*dR;
Fin Pour

Nous pouvons bien sûr arriver à de nombreux modèles suivant les hypothèse que nous faisons sur la maladie étudiée. Il sera souvent intéressant de définir plus de catégories de populations. Ainsi, pour des modèles de M.S.T. (Maladies Sexuellement Transmissibles), nous différencierons les populations suivant leur sexe (il faudrait aussi prendre l'âge en compte), ce qui multiplie le nombre de possibilités déjà par deux. Pour définir le système différentiel sous-jacent, il est devient ainsi préférable de modéliser le processus par un schéma :

Ainsi, pour étudier une infection de type SIDA (Syndrome d'Immunodéficience Acquise), nous pouvons par exemple faire un modèle dans lequel nous considérons quatre groupes de population : les individus sains X, les individus infectés Y, les porteurs sains Z, et enfin les individus malades A. Dans le schéma qui suit, les flèches entre les différents états représentent les transitions possible pour un individu donné avec les probabilités correspondantes :


  
(73)

En effet, du point de vue de l’étude du SIDA (et uniquement du SIDA), un individu donné, quel que soit son état, peut mourir d’autre chose que le SIDA, c’est ce que nous rangeons dans "mort naturelle". Nous considérons qu’il y a un taux de mortalité fixe, noté n, ce pour toutes les populations considérées (sains, infectés, malades ou porteurs sains). D’autre part, un individu sain X peut aussi être contaminé, auquel cas il devient un individu infecté Y. Un individu infecté peut lui évoluer vers deux stades possibles : malade A ou porteur sain Z.

Le nombre s représente la vitesse dévolution de la maladie : plus ce taux est grand, plus les individus infectés évoluent vite vers l’un des stades suivants. Le nombre p, compris entre 0 et 1, représente la probabilité pour que cette évolution se fasse vers l’état de malade. Bien sûr, si p diminue, la quantité de porteurs sains, immunisés contre la maladie, augmente. La maladie est alors moins dangereuse globalement. Enfin, un malade peut mourir de sa maladie, ce avec un taux de mortalité d considéré comme fixe, alors qu’un porteur sain, immunisé, ne peut plus évoluer du point de vue de cette maladie. Toutes ces hypothèses conduisent à poser des équations du système de la forme :

  (74)

où :

  (75)

En effet, ce coefficient représente un taux de contamination, qu'il est naturel de relier au rapport , probabilité lors d'une rencontre pour que la personne rencontrée soit infectée, donc contaminante.

MODÈLE PROIES-PRÉDATEURS (DE LOTKA-VOLTERRA)

Ce modèle d'interaction proies-prédateurs, a été proposé par Volterra après la première guerre mondiale. Il s'agissait d'expliciter la dynamique des population de sardines et de requins en mer Adriatique (et hop! un peu d'océanographie sur le site…); expliquer notamment pourquoi les quantités de sardines pêchées après l'interruption due à la guerre n'étaient plus aussi importantes que précédemment (ce qui peut sembler contre intuitif) et pourquoi à la reprise de la pêche la proportion observée de requins avait augmenté.

Reprenant ce que nous avions vu tout au début de ce chapitre et en notant N(t) le nombre d'individus proies et P(t) le nombre d'individus prédateurs nous avons (croissance exponentielle en absence de pêche) :

  (76)

et en l'absence de proies (dans le cas d'une pêche), le nombre de prédateurs diminuant exponentiellement, nous avons respectivement :

  (77)

A ce point du discours, nous devons considérer deux espèces (proies et prédateurs) qui ne sont pas isolés mais en interaction. Pour quantifier la contribution de l'interaction entre espèces, nous considérerons seulement la prédation, en assumant que sa valeur ou intensité (de l'interaction) est fonction de la probabilité de rencontre des proies-prédateurs qui sera elle supposée proportionnelle au produit du pourcentage des deux populations.

Les effets de ces rencontres n'ont pas les mêmes effets sur les deux espèces. Premièrement, bien sûr, chaque proie mangée par un prédateur est un gain net pour la population de ce dernier est une perte nette pour le premier. Ainsi, si l'effet des interactions est accepté comme étant proportionnel à , les signes de l'influence d'interaction différeront pour les deux espèces selon :

  (78)

Avant d'aller plus loin, cherchons les valeurs pour lesquelles les dérivées s'annulent (ce qui nous donnera au fait le point d'équilibre du système) :

  (79)

d'où :

  (80)

Une solution triviale est la "solution d'extinction" donnée par . Sinon, nous avons :

  (81)

Maintenant, normalisons les équations de Lotka-Volterra en écrivant (ainsi elles sont sans dimensions) :

  (82)

Avec cette normalisation, le modèle de Lotka-Volterra s'écrit :

  (83)

En réarrangeant les coefficients, le système s'écrit :

  (84)

pour lequel les dérivées s'annulent aux point unitaires.

Le tracé discret de ce système d'équations (dans lequel nous reconnaissons les termes logistiques vu plus haut), donne avec et les conditions initiales :


  
(85)

En comparaison, voici un exemple pratique (réel) de mesure de proies-prédateur (Lynx – Lièvre) par la Hudson Bay-Company :


  
(86)

Les deux diagrammes précédents représentent donc les variables x, y en fonction du temps. Cependant, ce qui peut être intéressant pour un scientifique (ou un écologiste…) c'est la représentation de y en fonction de x et inversement. Ainsi, nous obtenons pour les mêmes conditions initiales précédentes et pour des conditions diverses (les prédateurs sont en ordonnées et les proies en abscisses) :


  
(87)

Nous voyons ainsi (dans cette représentation de l'espace des phases) que pour des conditions initiales fixes, le système est périodique et a un point d'équilibre en . Ce qui correspond aux points où :

  (88)

Finalement, nous avons deux couples de points d'équilibre (c'est trivial en regardant le système d'équation) :

et   (89)

La question qui se poser est le sens de rotation (représentation) du plan des phases. Ainsi, en représentant les directions à l'aide d'un champ de vecteurs, nous obtenons la représentation (plus intéressante pour comprendre les observations de la pêche après la guerre…).


  
(90)

Pour savoir dans quelle direction nous nous dirigeons dans l'espace des phases à un moment donné, il suffit donc de connaître la dérivée dy/dx (ou réciproquement dx/dy). Ainsi nous avons :

  (91)

Ceci dit, nous voyons bien sur le diagramme des phases dans sa forme de champ de vecteurs qu'il arrive un moment dans le cycle de ce modèle que le nombre de prédateurs soit très élevé pour un faible nombre de proies. Donc le modèle mathématique (théorique) explique bien ce qui peut être à priori contre intuitif pour l'être humain. Il devient dès lors évident de deviner en quels points du cycle se trouvaient l'écosystème à la reprise de la pêche après la guerre.

Cependant, nous pouvons (devons) nous poser la question de ce qu'il se passe après un petite perturbation autour des points d'équilibres (ce qui est de la plus haute importance en écologie).

Nous avons donc le système :

  (92)

En y mettant une perturbation infiniment petite, celui-ci s'écrit :

  (93)

En négligeant les termes quadratiques nous obtenons :

  (94)

Dès lors, proche du point d'extinction , ce système s'écrit :

  (95)

Ce qui nous montre que proche du point d'extinction, le nombre d'individus des prédateurs diminue exponentiellement alors que les proies augmentent elles exponentiellement. Ceci à un sens biologique : quand il y a peu d'individus prédateurs (respectivement proies), les proies se multiplient alors qu'au fur et à mesure que le nombre de proies augmentent, les prédateurs se multiplient et se concentrent de plus en plus sur leurs proies (ahhh la nature…).

Proche du point d'équilibre , nous aurons :

  (96)

Ce cas n'est plus très trivial car nous avons alors des équations différentielles couplées. Pour résoudre ce système, différentions l'équation des proies par rapport au temps tel que :

  (97)

et en y injectant dy/dt (prédateurs) :

  (98)

Nous obtenons donc une petit équation différentielle du deuxième ordre (cf. chapitre de Calcul Différentiel Et Intégral). Dons la solution type est :

  (99)

En injectant cette solution dans l'équation différentielle, nous obtenons après simplification des exponentielles une simple équation du deuxième degré (cf. chapitre de Calcul Algébrique) :

  (100)

Dont la solution est triviale :

  (101)

Ainsi, la solution générale de l'équation différentielle est la combinaisons linéaire des deux solutions tel que :

  (102)

Mais nous avons donc :

  (103)

Dès lors, connaissant x(t) nous obtenons facilement :

  (104)

Utilisons maintenant la relation d'Euler (voir le traitement des nombres complexes dans le chapitre des nombres de la section arithmétique) :

  (105)

Ainsi, nous avons :

  (106)

et comme (cf. chapitre de Trigonométrie) , nous avons alors :

  (107)

et de manière similaires, nous obtenons pour les prédateurs :

  (108)

Ainsi, autour du point d'équilibre , les perturbations suffisamment petites pour valider la linéarisation (annulation des termes quadratiques) oscilles comme des ellipses (ou cercles) dont les axes sont définis par les deux équations précédentes.

MODÈLE DE HARDY-WEINBERG

Le modèle de Hardy-Weinberg (développé en 1908) décrit la diversité (plutôt que l'évolution) du vivant dans le sens de la génétique des populations. Cependant, ce modèle permet quand même de comprendre l'impact d'un agent de mutation sur une population d'individus et ainsi de résoudre les problèmes liés au refus de modèle de la sélection (évolution) naturelle de Darwin. Effectivement, les travaux de Darwin portaient sur des caractères qui variaient de façon continue. Or, il ne pouvait expliquer comment les individus se transmettent ces variations (c'est effectivement gênant). La génétique de Mendel, quant à elle imposait que seuls les caractères discontinus étaient héréditaires.

Il y avait donc conflit entre le modèle continu (de Darwin) et discontinu (de Mendel) et ce fut la génétique qui au vingtième siècle permis d'établir qu'au fait c'est le mélange de gênes qui crée la diversité et in extenso la sélection naturelle.

Remarque: Nous démontrerons en méthodes numériques (section d'informatique théorique) dans le cadre de l'étude des algorithmes génétiques que le chromosome ayant la valeur de fonction de "fitness" la plus grande sera celui qui aura le plus de chance de se reproduire. Ceci fait partie du théorème nommé "The Schema Theorem"...

Conclusion : L'évolution ne se fait pas à l'échelle de l'individu mais de la population (c'est une bonne leçon à retenir ça…)

Définitions:

D1. Un gène ou "génotype" est composé "d'allèles" (deux habituellement, donc nous considérerons les individus comme "diploïdes")

D2. Un "pool génique" est un ensemble de gênes présents dans une population

D3. Une "micro-évolution" est un changement de la fréquence des allèles d'un pool génétique d'une population.

D4. Un allèle est dit "allèle fixe" lorsque tous les membres d'un population portent deux allèles identiques.

D5. Une "spéciation" est l'apparition à long terme d'une micro-évolution faisant apparaître une nouvelle espèce.

D6. Une "espèce" est un ensemble de populations dont tous les membres ont le potentiel de se reproduire entre eux dans un environnement naturel (nous supposons dans ce modèle que les individus se reproduisent selon un mode sexué).

D6. La "mutation" est un agent évolutif qui modifie le pool génique d'une population en produisant de nouveaux gènes (nous considérons donc qu'il n'y a pas d'agents évolutifs à l'œuvre dans la population dans ce modèle).

Remarque: La mutation est un événement rare et le plus souvent néfaste. Ces effets quantitatifs sont plus importants chez les organismes à temps de génération court (bactéries et virus), l'effet étant nettement moins marqué chez les organismes à temps de génération plus long (animaux et végétaux).

Considérons maintenant une population composée de 13 individus (1 gêne, chacun composé de deux allènes pris parmi 3) :

  (109)

La fréquence de l'allène X1 est :

  (110)

La fréquence de l'allène X2 est :

  (111)

Enfin, la fréquence de l'allène X3 est :

  (112)

Et bien évidemment :

  (113)

Pour résumer, si nous considérons maintenant une population d'individus possédant des gènes composés de seulement deux type d'allène différents, nous avons dès lors (cas plus concret par rapport à la biologie du vivant sur notre Terre) :

  (114)

Nous avons donc pour la fréquence du nombre d'allènes B:

  (115)

et pour la fréquenced'allènes b :

  (116)

Avec donc :

  (117)

Nous pouvons faire quelques observations suite à la manière de construire ce modèle. Effectivement, la fréquence des allèles (équilibre génétique) dans une population demeure approximativement constante génération après génération si :

1. La population est de très grande taille (une perturbation infime ne change que très peu la fréquence des types génotypes)

2. Pas d’émigration ou d’immigration (il n'y pas de nouveaux type d'allènes dans la population)

3. Pas de mutations modifiant les allèles (cela contient le point 2)

4. Accouplements au hasard (non influencés par le type d’allèle étudié)

5. Pas de sélection naturelle (cela contient le point 4)

Il peut être pertinent (pour la culture générale) de mettre en évidence le fait que nous distinguons deux types d’accouplements non-aléatoires :

1. L’endogamie : croisement entre individus du voisinage, donc ayant des liens de parenté ; chez les espèces se dispersant peu, particulièrement dans le règne végétal.

2. L’homogamie : choix des partenaires qui se ressemblent pour certains caractères ; répandu chez les animaux.

L'endogamie et l'homogamie ne sont pas une cause directe de microévolution puisqu’ils ne modifient pas le pool génique d’une population. Cependant, si certains des individus issus de ces croisements non aléatoires ont plus de chance de s'accoupler que d'autres (par sélection naturelle), il s'ensuivra une modification de la fréquence allélique dans la population descendante et par conséquent, de la microévolution. En conséquence, les accouplements non aléatoires sont un facteur potentiel de microévolution.

Finalement, pour en revenir à notre modèle, la probabilité qu’un individu possède l’un des trois génotypes possibles (relativement à un génotype) dans la population (relativement à notre exemple) est :

- Probabilité qu'un individu soit BB :

  (118)

- Probabilité qu'un individu soit bb :

  (119)

- Probabilité qu'un individu soit Bb ou bB:

  (120)

Finalement, nous obtenons "l'équation d'Hardy-Weinberg" :

  (121)

Que nous pouvons représenter par le graphique :


  
(122)

Après cette étude, nous pouvons enfin (re)définir la "loi de Darwin" de manière à ce qu'elle se soustraie à la grande question (problématique) de la modification continue ou discrète :

La sélection naturelle est un agent évolutif qui modifie le pool génique d'une population en augmentant la fréquence des allèles produisant les phénotypes (génotypes définissant le comportement d'un individu) les mieux adaptés au milieu. Ainsi, les individus les mieux adaptés à leur environnement se reproduisent plus que les autres et contribuent davantage au patrimoine génétique de la descendance. La fréquence des bons gènes augmente graduellement dans la population, de génération en génération. Ainsi, la sélection naturelle dirige l'adaptation d'une population à son milieu en accumulant les génotypes qui favorisent la survie dans le milieu.

Remarques:

R1. Le modèle d'Hardy-Weinberg n'existe pas dans la nature. Effectivement, il est impossible que toutes les conditions nécessaires soient présentes en même temps. Cependant, c'est un modèle théorique qui permet d'évaluer s'il y a de la microévolution dans une population. Ainsi, on mesure la fréquence des allèles dans la population mère et l’on applique l'équation. Si la structure génétique de la population fille s'écarte de celle prédite par la loi, alors on sait qu'au moins un agent de l'évolution est à l'œuvre.

R2. Ce modèle théorique permet de prédire les coûts sociaux des maladies génétiques. Effectivement, si on connaît le pourcentage d'individus qui présentent une maladie récessive dans une population, on peut appliquer l'équation et évaluer le nombre de transmetteurs sains puis prévoir la probabilité d'apparition de la maladie dans les générations futures. Si celle-ci est grande (ce raisonnement n'est valable que dans un cadre libéral et capitaliste), on peut décider d'investir dans la recherche ou prévoir des budgets pour soigner les malades.

Exemples:

E1. Soit une génotype de séquence possible et considérons que . Nous aurons alors :

  (123)

Ainsi, dans une population de 1'000'000 d'individus, nous devrions avoir :

  (124)

E2. Si au Québec (6 millions d'habitant), une personne sur cinq à les yeux bleus , combien de personnes sont de génotype BB et combien sont Bb?

  (125)

Donc :

  (126)

Ainsi nous avons :

  (127)

E3. Nous introduisons dans un étant 1'000 grenouilles tachetées, homozygotes pour ce caractère (TT) et 250 grenouilles sans taches également homozygotes pour ce caractère (tt). Si nous laissons les grenouilles se reproduire pendant quelques années, en supposant que la population demeures stable, quel nombre de grenouilles TT, Tt, tt devrait-on alors observer ?

Nous avons donc une population de 1'250 individus dans laquelle nous avons 1'000 TT (2'000 allèles T) et 250 tt (500 allèles t). Donc :

  (128)

Si la population est de 1'250 individus, nous aurons donc :

  (129)

TAUX DE CROISSANCE AVEC LA TEMPÉRATURE

Pour de nombreuses espèces (mammifères, poissons, micro-organismes), il est raisonnable de considérer qu'en première approximation, la relation du taux de croissance de la population avec le température de l'environnement est un polynôme du seconde degré:

  (130)

où les valeurs a, b, c vont dépendre de l'espèce considérée.

Nous savons en plus qu'il existe pour chaque individu un optimum de croissance à une température donnée , ainsi que des températures minimale  et maximale  en deçà et au-delà desquelles il n'y a plus de croissance.

Ainsi, les paramètres a, b, c doivent vérifier les relations suivantes:

  (131)

Ce qui implique pour les trois températures cardinales  de respecter les relations démontrées dans le cadre de l'étude des polynômes du deuxième degré:

  (132)

Donc si nous connaissons pour un objet vivant  et  il est facile d'en déduire théoriquement une approximation de .

 
 
  nombre de visiteurs venus 745135 visiteurs (2556056 hits) Ici!

Tracked by Histats.com
Recherche personnalisée
$value) { if ($param == 'client') { google_append_url($google_ad_url, $param, 'ca-mb-' . $GLOBALS['google'][$param]); } else if (strpos($param, 'color_') === 0) { google_append_color($google_ad_url, $param); } else if ((strpos($param, 'host') === 0) || (strpos($param, 'url') === 0)) { google_append_url($google_ad_url, $param, $google_scheme . $GLOBALS['google'][$param]); } else { google_append_globals($google_ad_url, $param); } } google_append_url($google_ad_url, 'dt', round(1000 * array_sum(explode(' ', microtime())))); return $google_ad_url; } $google_ad_handle = @fopen(google_get_ad_url(), 'r'); if ($google_ad_handle) { while (!feof($google_ad_handle)) { echo fread($google_ad_handle, 8192); } fclose($google_ad_handle); } ?>
 
 
Ce site web a été créé gratuitement avec Ma-page.fr. Tu veux aussi ton propre site web ?
S'inscrire gratuitement