3. La forme de Bussière

Bussière raisonne non pas sur les densités de population mais sur la population cumulée, P(x), dont la forme est la suivante :

message URL form42.gif

La méthode d’estimation diffère de celle retenue pour la fonction exponentielle négative. Compte tenu de sa forme, l’ajustement se fait par une méthode itérative d’optimisation qui permet de déterminer les deux paramètres D0 et γ.

La population des communes des communes est agrégée par pas de distance au centre de 500 mètres. Un pas inférieur, de l’ordre de 250 mètres, ne serait pas systématiquement renseigné faute de communes. Un pas supérieur de l’ordre de 1 000 mètres, les cumuls sont plus grossiers et conduisent à une perte de qualité d’ajustement.

Tableau 4.10 : Résultats de la forme fonctionnelle proposée par Bussière sur les aires urbaines de 1975-1990
NOM AIRE URBAINE D075 γ75 D082 γ82 D090 γ90
1 PARIS 34 308,0 -0,151 0,98 30 000,0 -0,140 0,98 25 000,0 -0,124 0,98
2 LYON 20 000,0 -0,317 0,90 17 000,0 -0,290 0,92 17 100,0 -0,284 0,91
3 MARSEILLE-AIX-EN-PROVENCE 280 000,0 -1,264 0,72 270 000,0 -1,235 0,65 252 849,9 -1,210 0,58
4 LILLE 15 000,0 -0,304 0,95 13 000,0 -0,279 0,97 10 571,1 -0,244 0,98
5 BORDEAUX 6 000,0 -0,233 0,95 5 500,0 -0,216 0,95 5 000,0 -0,196 0,96
6 TOULOUSE 10 000,0 -0,322 0,48 9 000,0 -0,296 0,63 8 600,0 -0,270 0,71
7 NANTES 9 900,0 -0,353 0,89 9 100,0 -0,329 0,91 7 500,0 -0,285 0,93
8 NICE 15 000,0 -0,449 0,39 13 000,0 -0,408 0,50 12 000,0 -0,377 0,57
9 STRASBOURG 10 000,0 -0,379 0,81 9 000,0 -0,352 0,84 8 000,0 -0,322 0,85
10 GRENOBLE 10 000,0 -0,385 0,92 9 000,0 -0,359 0,93 8 500,0 -0,341 0,94
11 ROUEN 11 714,8 -0,414 0,97 11 000,0 -0,402 0,97 10 500,0 -0,387 0,97
12 TOULON 5 748,2 -0,312 0,93 5 200,0 -0,282 0,94 4 900,0 -0,265 0,95
13 RENNES 5 000,0 -0,312 0,68 4 500,0 -0,281 0,75 4 100,0 -0,254 0,8
14 NANCY 13 000,0 -0,486 0,87 12 000,0 -0,462 0,88 11 100,0 -0,439 0,89
15 MONTPELLIER 18 000,0 -0,663 0,57 15 000,0 -0,566 0,62 14 000,0 -0,506 0,67
16 VALENCIENNES 2 143,5 -0,180 0,97 2 300,0 -0,193 0,97 2 127,1 -0,187 0,97
17 TOURS 7 000,0 -0,395 0,81 6 500,0 -0,366 0,86 6 000,0 -0,340 0,89
19 CLERMONT-FERRAND 18 605,7 -0,654 0,56 17 000,0 -0,617 0,62 16 500,0 -0,606 0,65
20 SAINT-ETIENNE 54 606,5 -1,057 0,67 48 000,0 -1,014 0,68 46 100,0 -0,997 0,69
21 LENS 3 232,2 -0,219 0,97 3 081,9 -0,215 0,97 2 869,1 -0,206 0,7
22 CAEN 7 000,0 -0,415 0,69 6 500,0 -0,402 0,65 6 000,0 -0,369 0,7
23 ORLEANS 9 000,0 -0,496 0,79 7 700,0 -0,439 0,84 7 300,0 -0,405 0,88
24 DIJON 8 700,0 -0,468 0,69 8 000,0 -0,440 0,76 7 100,0 -0,396 0,77
25 HAVRE 77 839,6 -1,345 0,86 68 000,0 -1,277 0,81 62 000,0 -1,216 0,79
26 ANGERS 9 000,0 -0,508 0,82 8 300,0 -0,468 0,84 7 800,0 -0,434 0,85
27 BETHUNE 821,8 -0,118 0,96 820,1 -0,119 0,96 845,7 -0,121 0,97
28 MULHOUSE 12 353,7 -0,563 0,87 10 000,0 -0,499 0,88 8 000,0 -0,438 0,89
29 BREST 130 000,0 -1,938 0,68 110 000,0 -1,767 0,62 100 000,0 -1,676 0,57
30 REIMS 56 100,6 -1,265 0,80 52 000,0 -1,194 0,77 48 000,0 -1,121 0,76
31 DUNKERQUE 3 600,0 -0,314 0,75 3 189,4 -0,284 0,78 3 200,0 -0,285 0,79
32 MANS 119 563,5 -1,884 0,65 114 000,0 -1,819 0,61 110 000,0 -1,772 0,58
33 METZ 6 662,9 -0,434 0,76 6 000,0 -0,400 0,78 5 500,0 -0,372 0,78
34 LIMOGES 70 000,0 -1,532 0,71 65 000,0 -1,446 0,68 59 000,0 -1,372 0,65
35 AMIENS 10 000,0 -0,584 0,55 8 000,0 -0,511 0,57 6 000,0 -0,433 0,58
36 DOUAI 4 280,0 -0,352 0,92 4 150,0 -0,355 0,92 4 125,0 -0,359 0,93
37 AVIGNON 25 400,2 -1,022 0,40 22 000,0 -0,919 0,46 20 000,0 -0,857 0,48
38 BESANCON 10 000,0 -0,620 0,69 8 300,0 -0,559 0,62 7 700,0 -0,523 0,63
39 PERPIGNAN 60 941,6 -1,661 0,79 51 000,0 -1,429 0,75 43 000,0 -1,290 0,71
40 NIMES 30 000,0 -1,149 0,84 25 000,0 -1,034 0,76 21 000,0 -0,911 0,69
41 BAYONNE 2 216,8 -0,299 0,86 2 273,3 -0,292 0,88 2 308,5 -0,282 0,91
42 PAU 4 764,9 -0,442 0,95 4 700,0 -0,420 0,95 4 224,9 -0,387 0,95
43 GENEVE(CH)-ANNEMASSE 628,4 -0,174 0,85 638,8 -0,163 0,85 660,0 -0,149 0,89
44 LORIENT 3 206,3 -0,358 0,83 2 800,0 -0,323 0,88 2 727,0 -0,308 0,91
45 MONTBELIARD 1 584,3 -0,220 0,98 1 583,7 -0,220 0,98 1 610,0 -0,230 0,98
46 THIONVILLE 1 500,0 -0,197 0,93 1 406,5 -0,199 0,94 1 410,7 -0,203 0,94
47 POITIERS 4 500,0 -0,461 0,76 4 000,0 -0,416 0,80 3 600,0 -0,381 0,82
NOM AIRE URBAINE D075 γ75 D082 γ82 D090 γ90
48 TROYES 9 174,9 -0,630 0,86 8 000,0 -0,581 0,88 6 000,0 -0,500 0,89
49 ANNECY 12 000,0 -0,800 0,80 11 000,0 -0,730 0,83 10 100,0 -0,658 0,85
50 SAINT-NAZAIRE 4 071,7 -0,442 0,87 3 800,0 -0,417 0,87 3 200,0 -0,378 0,88
51 VALENCE 16 015,9 -0,941 0,71 13 500,0 -0,843 0,69 12 000,0 -0,778 0,68
52 ANGOULEME 2 681,0 -0,367 0,91 2 509,7 -0,344 0,92 2 200,0 -0,317 0,93
53 ROCHELLE 25 747,9 -1,186 0,73 21 000,0 -1,046 0,71 16 000,0 -0,903 0,71
54 BOURGES 8 000,0 -0,693 0,82 7 100,0 -0,637 0,82 6 500,0 -0,599 0,82
55 BOULOGNE-SUR-MER 5 092,2 -0,515 0,93 4 800,0 -0,499 0,93 4 200,0 -0,466 0,94
56 MAUBEUGE 2 751,9 -0,374 0,92 2 747,2 -0,373 0,92 2 700,0 -0,372 0,92
57 CHAMBERY 4 500,0 -0,540 0,73 3 996,1 -0,488 0,80 4 000,0 -0,469 0,83
58 CHALON-SUR-SAONE 6 000,0 -0,624 0,69 5 000,0 -0,556 0,73 3 930,1 -0,482 0,79
60 CHERBOURG 7 000,0 -0,690 0,98 6 200,0 -0,628 0,97 5 700,0 -0,570 0,97
61 HAGONDANGE-BRIEY 1 521,7 -0,265 0,93 1 450,0 -0,275 0,94 1 400,0 -0,279 0,94
62 CHARTRES 7 350,8 -0,735 0,94 6 500,0 -0,657 0,95 5 900,0 -0,593 0,94
63 MELUN 5 447,4 -0,596 0,92 5 100,0 -0,546 0,93 5 150,0 -0,506 0,93
64 CALAIS 17 753,5 -1,048 0,91 15 000,0 -0,957 0,90 13 000,0 -0,884 0,9
65 NIORT 4 500,0 -0,560 0,80 3 673,9 -0,489 0,81 3 000,0 -0,430 0,83
66 BEZIERS 22 253,1 -1,189 0,89 18 000,0 -1,090 0,86 16 000,0 -1,036 0,81
67 SAINT-BRIEUC 31 982,2 -1,522 0,80 27 000,0 -1,373 0,76 25 000,0 -1,320 0,73
68 COLMAR 3 319,5 -0,471 0,83 2 950,0 -0,439 0,84 2 600,0 -0,404 0,85
69 ARRAS 5 669,7 -0,602 0,93 5 000,0 -0,556 0,95 4 800,0 -0,540 0,96
70 BLOIS 3 737,4 -0,533 0,84 2 524,8 -0,417 0,86 2 504,4 -0,398 0,88
71 TARBES 12 811,9 -0,936 0,79 10 500,0 -0,839 0,79 9 000,0 -0,783 0,79
72 ROANNE 9 200,8 -0,780 0,85 10 200,0 -0,840 0,81 10 999,9 -0,889 0,76
73 SAINT-QUENTIN 8 600,0 -0,745 0,83 8 200,0 -0,735 0,84 8 000,0 -0,741 0,84
74 CREIL 5 000,0 -0,600 0,96 5 200,0 -0,591 0,96 4 793,5 -0,560 0,97
75 VANNES 3 000,0 -0,530 0,88 2 665,1 -0,462 0,89 2 500,0 -0,417 0,9
76 BELFORT 14 001,9 -1,016 0,94 12 000,0 -0,933 0,94 10 100,0 -0,845 0,93
77 CHARLEVILLE-MEZIERES 14 000,0 1,041 0,73 13 500,0 -1,024 0,72 13 000,0 -1,001 0,68
78 SAINT-OMER 969,6 -0,271 0,98 918,2 -0,254 0,98 880,0 -0,243 0,98
79 NEVERS 3 000,0 -0,476 0,88 2 249,2 -0,403 0,89 2 060,9 -0,382 0,9
80 QUIMPER 6 663,0 -0,768 0,94 6 900,0 -0,753 0,93 7 600,0 -0,768 0,92
81 EVREUX 10 000,0 -0,967 0,98 9 100,0 -0,900 0,86 7 800,0 -0,781 0,85
82 BEAUVAIS 10 000,0 -0,936 0,83 9 700,0 -0,922 0,81 9 900,0 -0,899 0,79
83 LAVAL 9 200,0 -0,903 0,91 8 400,0 -0,835 0,88 7 800,0 -0,780 0,88
84 COMPIEGNE 3 738,7 -0,582 0,86 3 649,4 -0,548 0,86 3 700,0 -0,531 0,87
85 ROCHE-SUR-YON 4 256,7 -0,649 0,80 4 499,9 -0,635 0,78 4 500,0 -0,615 0,76
86 SAINT-CHAMOND 2 981,2 -0,484 0,92 2 950,0 -0,475 0,93 2 800,0 -0,465 0,92
87 CHATEAUROUX 4 689,9 -0,622 0,94 4 200,0 -0,580 0,93 3 800,0 -0,546 0,93
88 BOURG-EN-BRESSE 36 816,9 -1,966 0,42 34 500,0 -1,853 0,39 33 000,0 -1,770 0,38
89 PERIGUEUX 5 000,0 -0,662 0,89 4 350,0 -0,603 0,90 3 890,0 -0,560 0,92
90 ALES 5 000,0 -0,663 0,93 4 600,0 -0,620 0,94 4 380,0 -0,590 0,93
91 MACON 4 342,4 -0,629 0,62 4 050,0 -0,593 0,68 3 900,0 -0,577 0,71
92 AGEN 3 100,0 -0,534 0,84 2 433,4 -0,458 0,87 2 100,0 -0,408 0,88
93 BRIVE-LA-GAILLARDE 4 000,0 -0,599 0,92 3 955,0 -0,580 0,93 3 700,0 -0,550 0,85
94 ALBI 3 864,7 -0,595 0,92 3 500,0 -0,552 0,92 3 011,3 -0,497 0,92
95 EPINAL 5 500,0 -0,702 0,81 5 200,0 -0,681 0,81 4 920,0 -0,656 0,82
96 CHALONS-SUR-MARNE 9 000,0 -0,894 0,95 8 670,0 -0,866 0,92 8 200,0 -0,850 0,91
98 AUXERRE 9 319,5 -1,010 0,75 9 300,0 -0,966 0,73 9 300,0 -0,951 0,72
100 CARCASSONNE 7 000,0 -0,879 0,80 6 700,0 -0,847 0,77 7 000,0 -0,821 0,77
109 DIEPPE 4 000,0 -0,669 0,91 3 300,0 -0,604 0,90 3 290,0 -0,593 0,9
118 ALENCON 5 600,0 -0,840 0,68 5 000,0 -0,765 0,75 4 800,0 -0,745 0,78
121 PUY-EN-VELAY 7 433,7 -0,968 0,90 7 000,0 -0,920 0,89 6 650,0 -0,880 0,88
135 SENS 10 966,3 -1,334 0,88 10 100,0 -1,228 0,87 9 800,0 -1,173 0,85
NOM AIRE URBAINE D075 γ75 D082 γ82 D090 γ90
139 LAON 2 923,6 -0,662 0,79 2 500,0 -0,606 0,81 2 125,0 -0,551 0,87
151 LONS-LE-SAUNIER 3 514,1 -0,773 0,83 2 797,5 -0,678 0,84 2 400,0 -0,620 0,84
153 LISIEUX 7 416,3 -1,146 0,82 7 000,0 -1,105 0,80 6 600,0 -1,076 0,78
155 DOLE 9 757,4 -1,291 0,91 9 300,0 -1,290 0,89 8 854,0 -1,270 0,92
163 SAINT-LO 11 182,9 -1,513 0,72 11 104,0 -1,465 0,68 10 978,0 -1,420 0,75
184 VITRY-LE-FRANCOIS 4 000,0 -0,916 0,90 3 900,0 -0,920 0,90 3 800,0 -0,920 0,91
186 BAR-LE-DUC 5 401,3 -1,093 0,87 5 578,0 -1,126 0,85 5 300,0 -1,105 0,84
218 BAYEUX 2 418,7 -0,861 0,92 2 497,0 -0,824 0,92 2 490,0 -0,807 0,92
242 LOURDES 4 842,3 -1,197 0,90 4 625,0 -1,170 0,89 4 523,0 -1,180 0,87
254 SAINT-GAUDENS 1 890,0 -0,811 0,67 1 850,0 -0,803 0,68 1 813,0 -0,791 0,68
294 GRAY 890,0 -0,545 0,89 850,0 -0,550 0,94 838,9 -0,563 0,95
310 LANGRES 2 838,7 -1,101 0,90 2 798,7 -1,111 0,89 2 623,0 -1,090 0,87

Plusieurs commentaires sont envisageables sur la base de ces estimations.

Les valeurs obtenues tant pour la densité centrale estimée que pour le gradient de densité en valeur absolue diminuent dans le temps, correspondant bien aux deux composantes de l’évolution de ces configurations : diminution de la population résidant au centre et dispersion de celle-ci en périphérie. La baisse des densités centrales estimées obtenues avec la forme de Bussière ne se retrouve pas sur la base des estimations sur les densités, pour lesquelles nous avions une tendance à la hausse. Contrairement aux résultats obtenus sur les densités, l’étalement urbain est manifeste au cours de la période 1975-1982 avec des taux de variation de gradient négatifs. La variation du gradient est plus ample que celle enregistrée pour la période ultérieure.

Les gradients moyens comme les densités centrales pour l’ensemble des aires urbaines sont significativement supérieurs à ceux obtenus sur la base des estimations retenant les densités.

Tableau 4.11 : Densités centrales estimées et gradients moyens selon des exponentielles sur la population cumulée et sur les densités
1975 1982 1990 82/75 90/82
Exponentielle sur densités D 0 329,26 398,36 442,46 16,0% 7,6%
γ -0,1424 -0,1443 -0,1417 1,3% -1,8%
Exponentielle sur population cumulée D 0 15 407,38 13 961,48 12 873,51 -9,4% -7,8%
γ -0,697 -0,671 -0,636 -3,7% -5,2%

Les spécificités des gradients et des densités centrales estimées selon la taille des aires urbaines mises à jour au niveau de la forme fonctionnelle sur les densités se retrouvent en grande partie pour les estimations privilégiant la population cumulée.

La relation croissante entre taille des aires urbaines et densité centrale estimée reste valable (Tableau 4.12). Plus les aires sont peuplées, plus la densité centrale moyenne est importante.

A contrario, les taux de variation de ces densités estimées sur les deux périodes ont un sens de variation opposé. Avec les estimations retenant la population cumulée, les variations sont décroissantes. Cette diminution est particulièrement significative sur les aires urbaines de taille intermédiaire. Comme précédemment mais à la baisse cette fois, on note une convergence des taux de variation entre les différentes aires urbaines sur les deux périodes.

Tableau 4.12 :Valeurs des densités centrales moyennes selon la taille des aires urbaines (estimations sur la population cumulée)
Taille des aires urbaines 1975 1982 1990 82/75 90/82
< ;50 000 5 233,97 4 992,32 4 780,38 -4,6% -4,2%
50 000< ;x< ;100 000 7 195,65 6 731,79 6 406,28 -6,4% -4,8%
100 000< ;x< ;150 000 9 329,57 8 069,13 7 265,92 -13,5% -10,0%
150 000< ;x< ;300 000 24 818,83 22 052,21 20 039,71 -11,1% -9,1%
300 000< ;x< ;500 000 12 696,22 11 367,84 10 768,94 -10,5% -5,3%
> ;500 000 44 467,56 41 733,33 38 513,44 -6,1% -7,7%
Ensemble 15 407,38 13 961,48 12 873,52 -9,4% -7,8%

Pour les gradients moyens de densité, la relation décroissante entre niveau du gradient et taille de l’aire urbaine est également observée au niveau de ces estimations. Les gradients diminuent à mesure que l’aire urbaine est peuplée. Les grandes aires urbaines seraient davantage étalées que les petites. On note comme précédemment une exception à cette observation les aires dont la population est comprise entre 50 000 et 100 000 habitants.

La tendance à la dispersion des populations est particulièrement manifeste pour les aires urbaines de grande taille (>500 000 hab.) et de taille intermédiaire (150 000 – 300 000 hab.) pour lesquels les taux de variation des gradients à la baisse sont les plus élevés.

Tableau 4.13 : Valeurs des gradients moyens de densités selon la taille des aires urbaines (estimations sur la population cumulée)
Taille des aires urbaines (hab.) 1975 1982 1990 82/75 90/82
< ;50000 -1,02 -0,99 -0,97 -2,8% -2,4%
50 000< ;x< ;100 000 -0,69 -0,72 -0,69 4,4% -4,7%
100 000< ;x< ;150 000 -0,71 -0,66 -0,62 -7,9% -5,8%
150 000< ;x< ;300 000 -0,74 -0,68 -0,64 -7,7% -6,8%
300 000< ;x< ;500 000 -0,46 -0,43 -0,41 -6,5% -5,3%
> ;500 000 -0,42 -0,39 -0,37 -6,0% -6,6%
Ensemble -0,70 -0,67 -0,64 -3,7% -5,2%

La qualité d’ajustement est relativement proche selon la taille des aires urbaines.

Les R² sont très bons et même meilleurs que ceux obtenus dans le cadre de l’exponentielle sur les densités. Il est possible de les améliorer en lissant la courbe de la population cumulée par la ventilation de la population de communes importantes sur plusieurs pas de distance, en particulier celle de la commune centre.

Cependant quelle que soit la méthode de répartition de la population, cette bonne performance est imputable à la construction de la variable explicative qui correspond au cumul de valeurs qui se suivent. La sommation des valeurs successives en fonction de la distance au centre équivaut à réduire la dispersion du nuage de points. A chaque pas de distance correspond une et une seule valeur.

Pour confronter les résultats entre ces deux méthodes d’estimation, l’une fondée sur les densités, l’autre sur la population cumulée, il est nécessaire en conséquence d’ajuster la fonction exponentielle sur la base de surface et de population agrégée par pas de distance pour déterminer des densités moyennes, comme cela est fait pour l’estimation de l’exponentielle sur la base de la population cumulée. Les résultats suivant cette méthode d’estimation sont présentés dans leTableau 4.14.

Tableau 4.14 : Résultats des estimations de la fonction exponentielle de densité sur la base de densités agrégées par pas de distance
Nom des aires urbaines 1975 1982 1990
Log D o γ Log D o γ Log D o γ
PARIS 4,65
(27,0)
-0,91
(-24,9)
0,88
4,66
(28,6)
-0,87
(-24,9
0,89
4,67
(29,8)
-0,83
(-4,2)
0,89
LYON 3,16
(13,2)
-0,11
(-13,9)
0,77 3,15
(12,9)
-0,99
(-10,3)
0,73 3,17
(13,2))
-0,92
(-9,7)
0,70
MARSEILLE-AIX-EN-PROVENCE 1,86
(4,7)
-0,83
(-4,6)
0,41 2,10
(5,9)
-0,81
(-5,0)
0,44 2,25
(6,6)
-0,76
(-5,6)
0,43
LILLE 3,96
(16,3)
-0,15
(-10,3)
0,79
3,98
(17,9)
-0,15
(-10,6)
0,81 3,97
(19,0)
-0,14
(-10,2)
0,82
BORDEAUX 2,78
(10,3)
-0,12
(-13,1)
0,79
2,95
(11,9)
-0,12
(-12,8)
0,82 3,04
(13,0)
-0,12
(-13,4)
0,83
T de Student entre parenthèses

Les résultats de l’ajustement de la fonction exponentielle négative sur ces densités calculées et agrégées par pas de distance sont améliorés par rapport aux précédents résultats de la même fonction mais déterminée sur la base des différentes densités communales en fonction de la distance au centre. Cette amélioration porte principalement sur le R² plus que sur les T de Student. Les densités centrales estimées adoptent la même tendance d’évolution à la hausse que celles obtenues avec une estimation exponentielle négative sur densités. Le gradient en revanche est orienté à la baisse sur la période comme cela apparaît dans le cadre de l’estimation sur la population cumulée.

Dans l’estimation de la forme fonctionnelle proposée par Bussière proposée, il apparaît également que les résidus de cet ajustement, ordonnés suivant la distance, traduisent un phénomène d’autocorrélation. Si l’on considère que l’ensemble des points de l’espace sont ramenés à une succession de points ordonnés dans l’espace comme les dates le sont dans des séries temporelles, le test de Durbin-Watson utilisé pour ces séries chronologiques offre alors une première mesure de cette autocorrélation. Les valeurs de ce test sont proches de 0,2 et confirment la diagnostic réalisé sur la base des graphiques (Greaphique 4.6).

message URL graphique46.gif
Graphique 4.6 : Résidus en fonction de la distance au centre pour des estimations de la fonction exponentielle sur des populations cumulées

Enfin, la méthode d’estimation fondée sur une recherche itérative d’optimisation pour déterminer les deux paramètres D0 et γ est très sensible aux valeurs d’initialisation retenues.

Outre la forme exponentielle négative sur les densités et sur la population cumulée, la littérature, comme nous l’avons mentionné dans le chapitre sur l’état de l’art envisage d’autres formes fonctionnelles, qui offrent sur les cas retenus généralement de meilleurs résultats d’estimation. Il n’est alors pas inintéressant d’envisager un comparatif entre des aires urbaines françaises. Les formes fonctionnelles proposées sont nombreuses et diverses. Parmi celles-ci nous retiendrons la forme Log normal pour laquelle nous procéderons à des estimations.