arxiv: v1 [physics.soc-ph] 25 Mar 2008

Pulic trnsport networks: empiricl nlysis nd modeling rxiv:0803.3514v1 [physics.soc-ph] 25 Mr 2008 C. von Ferer, 1, 2, T. Holovtch, 3, 1, Yu. Holovtch, 4, 5, nd V. Plchykov 4, 1 Applied Mthemtics Reserch Centre, Coventry University, Coventry CV1 5FB, UK 2 Theoretische Polymerphysik, Universität Freiurg, D-79104 Freiurg, Germny 3 Lortoire de Physique des Mtériux, Université Henri Poincré, Nncy 1, 54506 Vndœuvre les Nncy Cedex, Frnce 4 Institute for Condensed Mtter Physics, Ntionl Acdemy of Sciences of Ukrine, UA 79011 Lviv, Ukrine 5 Institut für Theoretische Physik, Johnnes Kepler Universität Linz, A-4040, Linz, Austri (Dted: Octoer 22, 2018) We use complex network concepts to nlyze sttisticl properties of urn pulic trnsport networks (PTN). To this end, we present comprehensive survey of the sttisticl properties of PTNs sed on the dt of fourteen cities of so fr unexplored network size. Especilly helpful in our nlysis re different network representtions. Within comprehensive pproch we clculte PTN chrcteristics in ll of these representtions nd perform comprtive nlysis. The stndrd network chrcteristics otined in this wy often correspond to fetures tht re of prcticl importnce to pssenger using pulic trffic in given city. Specific fetures re ddressed tht re unique to PTNs nd networks with similr trnsport functions (such s networks of neurons, cles, pipes, vessels emedded in 2D or 3D spce). Bsed on the empiricl survey, we propose model tht leit eing simple enough is cple of reproducing mny of the identified PTN properties. A centrl ingredient of this model is growth dynmics in terms of routes represented y self-voiding wlks. PACS numers: 02.50.-r, 07.05.Rm, 89.75.Hc I. INTRODUCTION The generl interest in networks of mn-mde nd nturl systems hs led to creful nlysis of vrious network instnces using empiricl, simultionl, nd theoreticl tools. The emergence of this field is sometimes referred to s the irth of network science [1, 2, 3, 4, 5]. In this pper, we use complex network concepts to nlyze the sttisticl properties of pulic trnsport networks (PTN) of lrge cities. These constitute n exmple of trnsporttion networks [3] nd shre generl fetures of these systems: evolutionry dynmics, optimiztion, emedding in two dimensionl (2D) spce. Other exmples of trnsporttion networks re given y the irport [6, 7, 8, 9, 10, 11, 12, 13], rilwy [14], or power grid networks [6, 15, 16]. While the evolution of PTN of given city is closely relted to the city growth itself nd therefore is influenced y numerous fctors of geogrphicl, historicl, nd socil origin, there is mple evidence tht PTNs of different cities shre common sttisticl properties tht rise due to their functionl purposes [17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. Some of these properties hve een nlyzed in former studies, however, the ojective of the present study is to present comprehensive survey of chrcteristics of PTNs nd to provide comprtive nlysis. Bsed on this empiricl survey we re in the C.vonFerer@coventry.c.uk holtrs@lpm.u-nncy.fr hol@icmp.lviv.u plchykov@icmp.lviv.u FIG. 1: One of the networks we nlyze in this study. The Los Angeles PTN consists of R = 1881 routes nd N = 44629 sttions, some of them re shown in this mp (color online). position to propose growth model tht cptures mny of the min (sttisticl) fetures of PTNs. A further distinct feture of our study is tht the PTNs we will consider re networks of ll mens of pulic trnsport of city (uses, trms, suwy, etc.) regrdless of the specific mens of trnsport. A numer of studies hve nlyzed specific su-networks of PTNs [17, 18, 19, 20, 23, 24, 27]. Exmples re the Boston [17, 18, 19, 20] nd Vienn [20] suwy networks nd the us networks of severl cities in Chin [24, 27]. However, ech prticulr trffic system (e.g. the network of uses or trms, or the suwy network) is not closed system: it is sugrph of wider trnsporttion system of

2 city, or s we cll it here, of PTN. Therefore to understnd nd descrie the properties of pulic trnsport in city s whole, one should nlyze the complete network, without restriction to specific prts. Indeed, for the cse of Boston it hs een shown tht chnging from the suwy system to the network suwy + us the network properties chnge drsticlly [18, 19]. Urn pulic trnsport networks of generl type hve so fr een nlyzed minly in two previous studies [21, 22]. In the first one, Ref. [21], the PTNs of Berlin, Düsseldorf, nd Pris were exmined, wheres the suject of Ref. [22] were pulic trnsport systems of 22 Polish cities. Ref. [21] concentrted on the scle-free properties. For the cities considered, the node degree distriution ws shown to follow power lw. Moreover power lws were found for numer of other specific fetures descriing the trffic lod on the PTN. However, the sttistics in this study ws too smll for definite conclusions. In Ref. [22] it ws found tht the node degree distriution my follow power lw or e descried y n exponentil function, depending on the ssumed network representtion. Besides, numer of other network chrcteristics (clustering, etweenness, ssorttivity) were extensively nlyzed. In the present pper, we nlyze PTNs of numer of mjor cities of the world (see tle I) [30, 31]. Our choice for this dt se ws motivted y the requirement to collect network smples from cities of different geogrphicl, culturl, nd economicl ckground. Our current nlysis extends former studies [21, 22] y considering cities with lrger pulic trnsport systems (the typicl numer of stops in the systems considered in Ref. [22] ws severl hundreds) s well s y systemticlly nlyzing different representtions. The ide of different network representtions nturlly rises in the network science [1, 2, 3, 4, 5]. For the PTN the primry network topology is given y the set of routes ech servicing n ordered series of given sttions (see Fig. 1 s n exmple). For the trnsporttion networks studied so fr minly two different neighorhood reltions were used. In the first one, two sttions re defined s neighors only if one sttion is the successor of the other in the series serviced y this route [18, 19]. In the second one, two sttions re neighors whenever they re serviced y common route [14]. We will exploit oth representtions in our study. Moreover, we introduce further nturl representtions (descried in detil in Section II) which mke the description of the PTNs of tle I comprehensive. In prticulr, this includes iprtite grph representtion of trnsporttion network tht reflects its intrinsic fetures [20, 24, 26]. There is nother reson to seek scle-free properties of PTNs considering lrger dt se of more cities with lrger pulic trnsport communictions involved. A currently well ccepted mechnism to explin the undnt occurrence of power lws is tht of preferentil ttchment or rich gets richer [32, 33, 34]. As fr s PTNs oviously re evolving networks, their evolution my e ex- City A P N R S Type Berlin 892 3.7 2992 211 29.4 BSTU Dlls 887 1.2 5366 117 59.9 B Düsseldorf 217 0.6 1494 124 28.5 BST Hmurg 755 1.8 8084 708 25.5 BFSTU Hong Kong 1052 7.0 2024 321 39.6 B Istnul 1538 11.1 4043 414 31.7 BST London 1577 8.3 10937 922 34.2 BST Los Angeles 1214 3.8 44629 1881 52.9 B Moscow 1081 10.5 3569 679 22.2 BEST Pris 2732 10.0 3728 251 38.2 BS Rome 5352 4.0 3961 681 26.8 BT Sõ Polo 1523 10.9 7215 997 58.3 B Sydney 1687 3.6 1978 596 16.3 B Tipei 2457 6.8 5311 389 70.5 B TABLE I: Cities nlyzed in this study. A: urn re (km 2 ); P : popultion (million inhitnts); N: numer of PTN sttions; R: numer of PTN routes; S: men route length. Type of trnsport which is tken into ccount in the PTN dtse: Bus, Electric trolleyus, Ferry, Suwy, Trm, Urn trin. pected to follow similr underlying mechnisms. However, scle-free networks hve lso een shown to rise when minimizing oth the effort for communiction nd the cost for mintining connections [35, 36]. Moreover, this kind of n optimiztion ws shown to led to smll world properties [37] nd to explin the ppernce of power lws in generl context [38]. Therefore, sclefree ehvior of PTNs my lso e relted to ovious ojectives to optimize their opertion. This pper is orgnized s follows. In the next section (II) we define different representtions in which the PTN will e nlyzed, sections III-V explore the network properties in these representtions. We seprtely nlyze locl chrcteristics, such s node degrees nd clustering coefficients (section III), nd glol chrcteristics, such s pth length distriutions nd centrlities (section IV). Specil ttention is pid to chrcteristics tht re unique to PTNs nd networks with similr construction principles. An exmple is given y the nlysis of sequences of routes which go in prllel long given sequence of sttions, feture we cll hrness effect. A description of correltions etween the properties of neighoring nodes in terms of generlized ssorttivities is performed in section V. Our findings for the sttistics of rel-word PTNs re supported y simultions of n evolutionry model of PTNs s displyed in section VI. Conclusions nd n outlook re given in section VII. Some of our results hve een preliminry nnounced in Ref. [25].

3 A B C D E E B A c F C D 1 2 3 e F A B C D E E B A d F C D 1 2 FIG. 2: : piece of pulic trnsport mp. Sttions A-F re serviced y the trm lines No 1 (solid line), No 2 (dshed line), No 3 (dotted line). Tking ll the lines to e indistinguishle, we cll such representtion L -spce. : L-spce. c: P -spce. All sttions tht elong to the sme route re connected. d: P-spce. e: C -spce. Ech route is represented y node, ech link corresponds to common sttion shred y the route nodes it connects. f: C-spce. 3 f F in L-spce represent ll sttions tht re within rech of single sttion trip. For nlyzing PTNs, the L-spce representtion hs een used in Refs. [18, 21, 22, 23, 27]. Extending the notion of L-spce one my either introduce multiple links etween nodes depending on the numer of services etween them or ssocite corresponding weight to single link. We will refer to such representtion s L -spce (c.f. Fig. 2). A prticulrly useful concept for the description of connectivity in trnsport networks which we refer to s P- spce [22] ws introduced in ref. [14] nd used in PTN nlysis in Refs. [20, 22, 27]. In this representtion the network is grph where sttions re represented y nodes tht re linked if they re serviced y t lest one common route. In P-spce representtion the neighorhood of given node represents ll sttions tht cn e reched without chnging mens of trnsport. The P-spce concept my e extended to include multiple or weighted links. Such representtion we refer to s P - spce (c.f. Figs. 2c nd 2d, correspondingly). A somewht different concept is tht of iprtite spce which is useful in the nlysis of coopertion networks [3, 57]. In this representtion which we cll B- spce oth routes nd sttions re represented y nodes [24, 25, 26]. Ech route node is linked to ll sttion nodes tht it services. No direct links etween nodes of sme type occur (see Fig. 3). Oviously, in B-spce the neighors of given route node re ll sttions tht it services while the neighors of given sttion node re ll routes tht service it. 1 2 3 II. PT NETWORK TOPOLOGY Although everyone hs n intuitive ide out wht PTN is, it ppers tht there re numerous wys to define its topology. Let us descrie some of them, defining different spces in which pulic trnsport networks will e nlyzed. A strightforwrd representtion of PT mp in the form of grph represents every sttion y node while ny two nodes tht re successively serviced y t lest one route re linked y n edge s shown in Fig. 2. Let us note, tht the full informtion out the network of N sttions nd R routes is given y the set of ordered lists ech corresponding to one route or to one of the two directions of given route. These simply list ll sttions serviced y tht route in the order of service etween two terminl sttions or in the course of round trip. Note tht multiple entries of given sttion in such list re possile nd do occur. Let us first introduce simple grph to represent this sitution. In the following we will refer to this grph s L-spce [22]. This grph represents ech sttion y node, link etween nodes indictes tht there is t lest one route tht services the corresponding sttion consecutively. No multiple links re llowed (see Fig. 2). The neighors of given node A B C D E F FIG. 3: A iprtite grph of trm lines (filled circles) nd sttions (circles) which corresponds to the pulic trnsport mp of Fig. 2. For the ske of illustrtion, lines corresponding to different trm routes re shown in differing wy. However, neither line type nor the order of the sttions mtter in this grph. Note tht Figs. 2c - 2f re the one-mode projections of this iprtite grph. We note tht the one mode projections of the iprtite grph of B-spce to the set of sttion nodes results in P- spce or in P -spce spce if we retin multiple links. The complementry projection to route nodes leds to grph which we cll C-spce (C -spce if multiple links re retined). In this spce ll nodes represent routes nd the neighors of ny route node re those routes with which it shres common sttion, see Figs. 2e, 2f. Below, we will study different fetures of the PT networks s they pper when represented in the ove defined spces. It is worthwhile to mention here, tht stndrd network chrcteristics eing represented in different spces turn out to e nturl chrcteristics one is

4 City k L z L l mx L l L C L c L k P z P l mx P l P C P c P k C z C l mx C l C C C c C Berlin 2.58 1.96 68 18.5 2.6 10 4 52.8 56.61 11.47 5 2.9 2.9 10 3 41.9 27.56 4.43 5 2.2 1.2 10 2 4.75 Dlls 2.18 1.28 156 52.0 1.4 10 5 55.0 100.58 11.23 8 3.2 5.9 10 3 48.6 11.09 3.45 7 2.7 9.2 10 1 5.34 Düsseldorf 2.57 1.96 48 12.5 8.6 10 3 24.4 59.01 10.56 5 2.6 1.2 10 3 19.7 32.18 2.47 4 1.8 4.9 10 1 2.23 Hmurg 2.65 1.85 156 39.7 1.4 10 5 254.7 50.38 7.96 11 4.7 1.4 10 4 132.2 17.51 4.49 10 4.0 9.9 10 2 28.3 Hong Kong 3.59 3.24 60 11.0 1.0 10 4 60.3 125.67 10.20 4 2.2 1.3 10 3 11.7 98.98 2.12 3 1.7 1.2 10 2 2.14 Istnul 2.30 1.54 131 29.7 5.7 10 4 41.0 76.88 10.59 6 3.1 4.2 10 3 41.5 52.81 3.86 5 2.3 2.6 10 2 5.00 London 2.60 1.87 107 26.5 1.4 10 5 320.6 90.60 16.97 6 3.3 1.2 10 4 90.0 49.91 6.80 6 2.6 7.4 10 2 11.1 Los Angeles 2.37 1.59 210 37.1 7.9 10 5 645.3 97.99 17.21 11 4.4 7.4 10 4 399.6 40.11 8.42 10 3.6 2.3 10 3 22.1 Moscow 3.32 6.25 27 7.0 1.1 10 4 127.4 65.47 26.48 5 2.5 2.7 10 3 38.0 109.37 4.57 4 1.9 3.2 10 2 3.59 Pris 3.73 5.32 28 6.4 1.0 10 4 78.5 50.92 24.06 5 2.7 3.1 10 3 59.6 39.95 4.67 4 1.9 1.1 10 2 2.72 Rome 2.95 2.02 87 26.4 5.0 10 4 163.4 69.05 11.34 6 3.1 4.2 10 3 41.4 59.40 4.86 5 2.5 5.1 10 2 7.04 Sõ Polo 3.21 4.17 33 10.3 3.4 10 4 268.0 137.46 19.61 5 2.7 6.0 10 3 38.2 151.72 4.25 4 2.0 5.2 10 2 4.27 Sydney 3.33 2.54 34 12.3 7.3 10 3 82.9 42.88 7.79 7 3.0 1.3 10 3 33.6 65.02 2.92 6 2.4 3.5 10 2 6.30 Tipei 3.12 2.42 74 20.9 5.3 10 4 186.2 236.65 12.96 6 2.4 3.6 10 3 15.4 93.33 2.95 5 1.8 1.6 10 2 2.44 TABLE II: PTN chrcteristics in different spces (suscripts refer to L, P, nd C-spces, correspondingly). k: node degree (nerest neighors numer z (1) ); z = z (2) / z (1) (z (2) eing the next nerest neighors numer); l mx, l : mximl nd men shortest pth length (10); C : etweenness centrlity (23); c: reltion of the men clustering coefficient to tht of the clssicl rndom grph of equl size (8). Averging hs een performed with respect to corresponding network, only the men shortest pth l is clculted with respect to the lrgest connected component. interested in when judging out the pulic trnsport of given city. To give n exmple, the verge length of shortest pth l in L-spce, l L gives the numer of stops one hs to pss on verge to trvel etween ny two sttions. When represented in P spce, l P tells out how mny chnges one hs to do to trvel etween ny two sttions. And, finlly, l C rings out the numer of chnges one hs to do to pss etween ny two routes. Another exmple is given y the node degree k: k L tells to how mny directions pssenger cn trvel t given sttion; k L is the numer of stops in the direct neighorhood; k P is the numer of other sttions rechle without chnging line; wheres k C tells how mny routes re directly ccessile from the given one. Tle II lists some of the PTN chrcteristics we otined for the cities under considertion using pulicly ville dt from the we pges of locl trnsport orgniztions [30, 31]. A detiled nlysis nd discussion is given in the following sections III - V. III. LOCAL NETWORK CHARACTERISTICS Let us first exmine locl properties of the PTNs under discussion. Insted of looking for chrcteristics of individul nodes we will e interested in their men vlues nd sttisticl distriutions. This pproch llows us to derive conclusions tht re significnt for the glol ehvior of the given network. The simplest ut highly importnt properties re those concerning the node degrees of network nd in prticulr their distriution. Erly ttempts to model complex networks were performed y mthemticins using the concept of rndom networks [40, 41] in which correltions re sent. A welth of insight ws gined y elorting the theory on rigorous grounds developing mny concepts which remin mong the core of network nlysis. A rndom grph is given y set of N nodes nd M links. The nodes to which the two ends of ech link re connected re chosen with constnt proility 2M/N. In cse tht multiple links re excluded the verge numer of neighors z 1 is equl to the verge node degree k which is: z 1 = k = 2M/N. (1) For the node degree k nd its moments k m the verge (1) cn lso e considered s n verge with respect to the node degree distriution p(k): k m = k mx k=1 p(k)k m, (2) with the ovious nottion k mx for the mximl node degree. In (2),... stnds for n ensemle verge over different network configurtions. In the following nlyzing empiricl dt we will often use the sme nottion for n verge over lrge network instnce. For clssicl rndom grphs of finite size the node degree distriution p(k) is inomil, in the infinite cse it ecomes Poisson distriution. In Figs. 4, 5 we show node degree distriutions for PTNs of severl cities in L, P, nd C-spces. Note tht to get smoother curves we plot in the cse of P nd C- spces the cumultive distriutions defined s: P (k) = k mx q=k p(q). (3)

5 c FIG. 4: : Node degree distriutions of PTN of severl cities in L-spce. : Cumultive node degree distriution in P-spce. c: Cumultive node degree distriution in C-spce. Berlin (circles, ˆk L = 1.24, ˆk P = 39.7), Düsseldorf (squres, ˆk L = 1.43, ˆk P = 58.8), Hong Kong (strs, ˆk L = 2.50, ˆk P = 125.1). c FIG. 5: : Node degree distriutions of PTN of severl cities in L-spce. : Cumultive node degree distriutions in P-spce. c: Cumultive node degree distriution in C-spce. London (circles, γ L = 4.48, γ P = 4.39), Los Angeles (strs, γ L = 4.85, γ P = 3.92), Pris (squres, γ L = 2.62, γ P = 3.70). In Fig. 4 the dt is shown in log-liner plot together with fits (for L nd P-spces) to n exponentil decy: p(k) exp( k/ˆk), (4) where ˆk is of the order of the men node degree. Within the ccurcy of the dt oth L nd P-spce distriutions for the cities nlyzed in Fig. 4 re nicely fitted y n exponentil decy. As fr s the L-spce dt is concerned, we find evidence for n exponentil decy for out hlf of the cities nlyzed, while the other prt rther demonstrte power lw decy of the form: p(k) 1/k γ. (5) Figs. 5, 5 show the corresponding plots for three other cities on log-log scle. Numericl vlues of the fit prmeters ˆk nd γ (4), (5) for different cities re given in Tle III. There, rcketed vlues indicte less relile fit. Note tht for L-spce the fit ws done directly for the node degree distriution p(k), wheres due to n essentil scttering of dt in P-spce the cumultive distriution (3) ws fitted nd the corresponding vlues for the fit prmeters γ P, ˆk P were extrcted from those for the cumultive distriutions. While the node degree distriution of lmost hlf of the cities in the L-spce representtion disply power lw decy (5), this is not the cse for the P-spce. So fr, the nlysis of PTNs of smller cities never showed ny power-lw ehvior in P-spce [22, 27]. The dt for the three cities shown in Fig. 5 gives first evidence of power lw ehvior of P (k) in the P-spce representtion. Previous results concerning node-degree distriutions of PTNs in L nd P-spces seemed to indicte tht in generl the degree distriution is power-lw like in L- spce nd exponentil in P-spce. This ws interpreted [22] s indicting strong correltions in L-spce nd rndom connections etween the routes explining P-spce ehvior. Our present study, which includes much less homogeneous selection of cities (Ref. [22] ws sed on exclusively Polish cities) shows tht lmost ny comintion of different distriutions in L nd P-spces my occur. However, the three cities tht show power lw

6 City γ L ˆkL γ P ˆkP Berlin (4.30) 1.24 (5.85) 39.7 Dlls 5.49 (0.78) (4.67) 76.9 Du sseldorf (3.76) 1.43 (4.62) 58.8 Hmurg (4.74) 1.46 (4.38) 60.7 Hong Kong (2.99) 2.50 (4.40) 125.1 Istnul 4.04 (1.13) (2.70) 71.4 London 4.48 (1.44) 4.39 (143.3) Los Angeles 4.85 (1.52) 3.92 (201.0) Moscow (3.22) 2.15 (2.91) 50.0 Pris 2.62 (3.30) 3.70 (100.0) Rome 3.95 (1.71) (5.02) 54.8 Sõ Polo 2.72 (4.20) (4.06) 225.0 Sydney 4.03 (1.88) (5.66) 38.7 Tipei (3.74) 1.75 (5.16) 201.0 TABLE III: Prmeters of the PTN node degree distriutions fit to n exponentil (4) nd power lw (5) ehvior. Brcketed vlues indicte less relile fits. Suscripts refer to L nd P-spces [31]. distriution in P-spce lso exhiit power lw ehvior in L-spce, s one cn see compring Figs. 5 nd 5. In C-spce the decy of the node degree distriution is exponentil or fster, s one cn see from the plots in Fig. 4c nd 5c. From the cities presented there, only the PTNs of Berlin, London, nd Los Angeles re governed y n exponentil decy nd their node degree distriutions cn e pproximted y stright line in the figures. For most cities tht show power lw degree distriution in L-spce the corresponding exponent γ L is γ L 4. Note tht the exponents found for the PTNs of Polish cities of similr size N lso lie in this region: γ L = 3.77 for Krkow (with numer of sttions N = 940), γ L = 3.9 for Lodz (N = 1023), γ L = 3.44 for Wrsw (N = 1530) [22]. According to the generl clssifiction of scle-free networks [2] this indictes tht in mny respect these networks re expected to ehve similr to those with exponentil node degree distriution. Prominent exceptions to this rule re provided y the PTNs of Pris (γ L = 2.62) nd Sõ-Polo (γ L = 2.72). Furthermore, vlues of γ L in the rnge 2.5 3.0 were recently reported for the us networks of three cities in Chin: Beijing (N = 3938), Shnghi (N = 2063), nd Nnjing (N = 1150) [27]. The connectivity within the closest neighorhood of given node is descried y the clustering coefficient defined s C i = 2y i k i (k i 1), (6) where y i is the numer of links etween the k i nerest neighors of the node i. The clustering coefficient of node my lso e defined s the proility of ny two of its rndomly chosen neighors to e connected. For the men vlue of the clustering coefficient of rndom FIG. 6: Men clustering coefficient C P(k) of severl PTN in P-spce. Berlin (strs), London (tringles), Tipei (circles). grph one finds C R = k R N = 2M N 2. (7) In Tle II we give the vlues of the men clustering coefficient in L, P, nd C-spces. The highest solute vlues of the clustering coefficient re found in P-spce, where their rnge is given y C P = 0.7 0.9 (c.f. with C L = 0.02 0.1). This is due to the fct tht in this spce ech route gives rise to fully connected sugrph (complete grph). In order to mke numers comprle we normlize the vlue of C y the men clustering coefficient (7) of rndom grph of the sme size: c = N 2 C /(2M). (8) In L nd P-representtions we find the men clustering coefficient to e lrger y orders of mgnitude reltive to the rndom grph. This difference is less pronounced in C-spce indicting lower degree of orgniztion in these networks. Furthermore, we find these vlues to vry strongly within the smple of the 14 cities. This suggests tht the concepts ccording to which vrious PTNs re structured led to mesurle difference in their orgniztion. In P-spce the clustering coefficient of node is strongly correlted with the node degree. In Fig. 6 we show the men clustering coefficient of nodes of degree k, C P (k), s function of k for severl PTNs. Its ehvior cn e understood s follows. Recll tht the P- spce the degree of node (sttion) equls the numer of sttions tht cn e reched from given one. Ech route enters the network s complete grph, within which every node hs clustering coefficient of one. A smll numer k of neighors of given sttion indictes tht the sttion elongs to single route (i.e. C P (k)

7 c FIG. 7: Men shortest pth length distriution P (l) for the PTN of Berlin (strs), Hong Kong (circles), Pris (oxes) nd Rome (tringles). Solid line shows fit to the function (12). : L-spce; : P-spce; c: C-spce. is most proly equl to one). For nodes with higher degrees k it is more prole tht they elong to more thn one route. Consequently, C P (k) decreses with k. The chnge in the ehvior of C P (k) should occur t some vlue of k which is of the order of the men numer of stops of the routes. The prominent feture of the function C P (k) in P-spce is tht it decys following power lw C P (k) k β. (9) Within simple model of networks with str-like topology this exponent is found to e of vlue β = 1 [22]. In trnsport networks. This ehvior ws first oserved for the Indin rilwy network [14] nd then for the Polish PTNs [22]. In our cse, the vlues of the exponent β for the networks studied lie in the rnge from 0.65 (Sõ Polo) to 0.96 (Los Angeles) with men vlue of 0.82. IV. GLOBAL CHARACTERISTICS A. Pth length distriution Let l i,j e the length of shortest pth etween sites i nd j in given spce. The men shortest pth is defined s l = 2 N(N 1) N i>j=1 l ij. (10) Note tht l is well-defined only if nodes i nd j elong to the sme connected component of the network. In the following ny expression s given in Eq. (10) will e restricted to this cse. Furthermore, relted network chrcteristics will e clculted for the lrgest (or gint) connected component, GCC. Correspondingly, N denotes the numer of constituting nodes of this component. Denoting the pth length distriution s P (l), the verge (10) reds l = l mx l=1 P (l)l, (11) where l mx is mximl shortest pth length found on the connected component. In Fig. 7 we plot the men shortest pth length distriutions otined in different spces for severl selected cities. Together with the dt we plot fit to the symmetric unimodl distriution [22]: P (l) = Al exp ( Bl 2 + Cl), (12) where A, B, C re fit prmeters. As cn e seen from the figures, the dt is generlly nicely reproduced y this nstz. However, in certin networks dditionl fetures my led to devition from this ehvior s cn e seen from Fig. 8, which shows the men shortest pth length distriution in L-spce P L (l) for Los Angeles. One oserves second locl mximum on the right shoulder of the distriution. Qulittively this ehvior my e explined y ssuming tht the PTN consists of more thn one community. For the simple cse of one lrge community nd second smller one t some distnce this sitution will result in short intr-community pths which will give rise to glol mximum nd set of longer pths tht connect the lrger to the smller community resulting in dditionl locl mxim. Such sitution definitely ppers to e present in the cse of the Los Angeles PTN, see Fig. 1. Let us introduce chrcteristic tht informs how remote given node is from the other nodes of the networks. For the node i this my e chrcterized y the vlue: l i = 1 N 1 l ij. (13) j i Now, the men shortest pth (10) cn e defined in terms

8 FIG. 8: Men shortest pth length distriution in L-spce, P L(l), for the PTN of Los Angeles. of l i s: l = 1 l i. (14) N In order to look for correltions etween l i nd the node degree k i let us introduce the vlue: i l(k) = 1 N l i δ k,ki, (15) N k i=1 where N k is numer of nodes of degree k nd δ k,ki is the Kronecker delt. Consequently, l(k) is the men shortest pth length etween ny node of degree k nd other nodes of the network. For the mjority of the nlyzed cities the dependence of the men pth l L (k) (15) on the node degree k in L-spce cn e pproximted y power lw l L (k) k α L. (16) The vlue of the exponent vries in the rnge α L = 0.17 0.27. We show this dependence for severl cities in Fig. 9. A prticulr reltion etween pth lengths nd node degrees cn e shown to hold relting the men pth length etween two nodes to the product of their node degrees. To this end let us define l(k, q) = N l i,j δ kik j,kq. (17) i,j=1 As hs een shown in [42], this reltion cn e pproximted y l(k, q) = A B log(kq). (18) c FIG. 9: Men pth l L(k) (15) in the L-spce s function of the node degree k with fit to the power lw decy (16). : Berlin, α L = 0.23; : Hong Kong, α L = 0.25; c: Pris, α L = 0.15; d: Tipei, α L = 0.23. For rndom networks the coefficients A nd B cn e clculted exctly [43]. The vlidity of Eq. (18) ws checked on the se of PTNs of some Polish cities nd rther good greement for the mjority of the cities ws found in L-spce. In our nlysis which concerns PTNs of much lrger size, we do not oserve the sme good greement for ll cities. The suggested logrithmic dependence (18) ws found y us in L-spce lso for the lrger cities, however with much more pronounced sctter of dt for lrge vlues of the product kq. In Fig. 10 we plot the men pth l L (k, q) in the L-spce for the PTN of Berlin, Hong Kong, Rome, nd Tipei. Note, however, tht due to the sctter of dt logrithmic dependence frequently is indistinguishle from power lw with smll exponent. In P-spce, the shortest pth length l ij gives the miniml numer of routes required to e used in order to rech site j strting from the site i. In turn, l i, Eq. (13), defines the numer of routes one uses on verge trveling from the site i to ny node of the network. The higher the node degree, the esier it is to ccess other routes in the network. Therefore, lso in P-spce one expects decrese of l P (k) when k increses. This is shown for severl cities in Fig. 11. Besides the expected decrese of l P (k), one cn see tendency to power-lw decy l P (k) k α P. (19) The vlue of the exponent α P vries in the intervl α P = 0.09 (for Sydney) to α P = 0.17 (for Dlls) nd is centered round α P = 0.12 0.13 s shown for the cities in Fig. 11. The men pth l P (k, q) s function of d

9 FIG. 12: Men pth l P(k, q) in P-spce for PTN of Berlin () nd Pris () s function of kq. FIG. 10: Men pth l L(k, q) (17) in the L-spce s function of kq for the PTN of Berlin (strs), Hong Kong (circles), Rome (tringles), nd Tipei (squres). mesuring pth lengths to other nodes or on counting the numer of pths etween other nodes medited y this node. The closeness C c (i) [45] nd grph C g (i) [46] centrlities of node i re sed on the shortest pth lengths l ij to other nodes j: C c (i) = C g (i) = 1 j i l, i,j (20) 1. mx j i l i,j (21) c FIG. 11: Men pth l P(k) in P-spce s function of the node degree k nd its fit to the power lw decy (19). : Berlin, α P = 0.13; : Hong Kong, α P = 0.12; c: Pris, α P = 0.13; d: Tipei, α L = 0.12. kq for severl cities is given in P-spce in Fig. 12. The scttering of dt is much more pronounced thn in L- spce. However one distinguishes tendency of l P (k, q) to decrese with n increse of kq. The red lines in Figs. 12 re the guides to the eye chrcterizing the decy. B. Centrlities To mesure the importnce of given node with respect to different properties of grph numer of so-clled centrlity mesures hve een introduced [44, 45, 46, 47, 48]. Most of them re sed on either d Only nodes j tht elong to the sme connected component s i contriute to (20), (21). For given node these properties oviously depend on the size of the connected component to which the node elongs. The importnce of the node i with respect to the connectivity within the grph my e mesured in terms of the numer of shortest pths σ jk (i) etween nodes j nd k tht go vi node i. Denoting y σ jk the overll numer of shortest pths etween nodes j nd k one defines stress C s (i) [47] nd etweenness C (i) [48] centrlities y: C s (i) = σ jk (i), (22) C (i) = j i k j i k σ jk (i) σ jk. (23) Numericl vlues of the etweenness centrlity (23) re given in Tle I in L, P nd C-spces. Averging the two centrlities tht re sed on pth length (20), (21) one otins vlues tht re closely relted to the verge shortest pth length on the GCC. As fr s this reltion is independent of the representtion of the PTN, we find very similr correspondence etween l nd the men centrlities C c, C g in ll spces considered s shown in Fig. 13. The fct tht these centrlities re sed on the inverse pth length is reflected y the negtive slope of the curves shown in the figures. The etweenness centrlity (23) nd the relted stress centrlity (22) of given node mesure the shre of the men pths etween nodes tht re medited y tht node. It is ovious tht node with high degree hs

10 10 1 0.1 0.01 0.001 0.0001 100 10 1 0.1 0.01 0.001 100 10 1 0.1 1e-005 5 10 15 20 25 30 35 40 45 50 55 0.0001 2 2.5 3 3.5 4 4.5 5 0.01 1.5 2 2.5 3 3.5 4 4.5 c FIG. 13: Correspondence etween the men shortest pth l nd men centrlities C c (open circles), C g (filled circles) for ll fourteen PTN listed in Tle II in () L, () P, nd (c) C-spces. higher proility to e prt of ny pth connecting other nodes. This reltion etween C nd the node degree my e quntified y oserving their correltion. In Figs. 14 we plot the men etweenness centrlity C (k) of ll nodes tht hve given degree k. There, we present results for the PTN of Pris in L, C nd P, nd B-spces. Especilly well expressed is the etweenness-degree correltion in L-spce (Fig. 14) nd with somewht less precision in C-spce (Fig. 14). In oth cses there is cler tendency to power lw C (k) k η with n exponent η = 2 3. Let us note here, tht this power lw together with the scle free ehvior of the degree distriution implies tht lso the etweenness distriution should follow power lw with n exponent δ. This ehvior is clerly identified in Fig.15 for the L-spce etweenness distriution of the Pris PTN, for which we find n exponent δ 1.5. The resulting scling reltion [49] η = (γ 1)/(δ 1) (24) is fulfilled within the ccurcy for these exponents. In the plots for oth B nd P-spces we oserve the occurrence of two regimes which correspond to smll nd lrge degrees k. This seprtion however hs different origin in ech of these cses. In the B-spce representtion, the network consists of nodes of two types, route nodes nd sttion nodes. Typiclly, sttion nodes re connected only to low numer of routes while there is miniml numer of sttions per route. One my thus identify the low degree ehvior s descriing the etweenness of sttion nodes, while the high degree ehvior corresponds to tht of route nodes. In the overlp region of the two regimes one my oserve tht when hving the sme degree sttion nodes hve higher etweenness thn route nodes. The occurrence of two regimes in the P-spce representtion hs similr origin s the chnge of ehvior oserved for the men clustering coefficient C P (k), see Fig.6. Nmely, sttions with low degrees in generl elong only to single route nd thus re of low importnce for the connectivity within the network resulting in low etweenness centrlity. Compring our results with those of Ref. [22] we do not however find sturtion for the low k region, s oserved there. Similr etweenness c FIG. 14: Men etweenness centrlity C (k) - degree k correltions for the PTN of Pris in () L, () C, (c) P, nd (d) B-spces. FIG. 15: Betweenness centrlity C - distriution for PTN of Pris in L - spce. C (k) - degree reltions s oserved in Fig. 14 for the PTN of Pris we lso find for most of the other cities, however, with different qulity of expression. d

11 FIG. 16: Cumultive hrness distriutions for Istnul () nd for Moscow () PTN. FIG. 17: Cumultive hrness distriutions for Los Angeles. : log-log scle; : log-liner scle. C. Hrness Besides the locl nd glol properties of networks descried ove which cn e defined in ny type of network, there re some chrcteristics tht re unique for PTNs nd networks with similr construction principles. A prticulrly striking exmple is the fct tht s fr s the routes shre the sme grid of streets nd trcks often numer of routes will proceed in prllel long shorter or longer sequences of sttions. Similr phenomen re oserved in networks uilt with rel spce consuming links such s cles, pipes, neurons, etc. In the present cse this ehvior my e esily worked out on the sis of sequences of sttions serviced y ech route. To quntify this ehvior recently the notion of network hrness hs een introduced [25]. It is descried y the hrness distriution P (r, s): the numer of sequences of s consecutive sttions tht re serviced y r prllel routes. Similrly to the node-degree distriutions, we oserve tht the hrness distriution for some cities (Hong Kong, Istnul, Pris, Rome, Sõ Polo, Sydney) my e fitted y power lw: P (r, s) r γs, for fixed s, (25) wheres the PTNs of other cities (Berlin, Dlls, Düsseldorf, London, Moscow) re etter fitted to n exponentil decy: P (r, s) exp ( r/ˆr s ), for fixed s. (26) As exmples we show the hrness distriution for Istnul (Fig. 16) nd for Moscow (Fig. 16). Moreover, sometimes (we oserve this for Los Angeles nd Tipei), for lrger s the regime (25) chnges to (26). We show this for the PTN of Los Angeles in Fig. 17. There, one cn see tht for smll vlues of s the curves re etter fitted to power lw dependence (25). With incresing s tendency to n exponentil decy (26) ppers: Fig. 17. As one cn oserve from the Figs. 16, 17 the slope of the hrness distriution P (r, s) s function of the numer of routes r increses with n increse of the sequence length s. For PTNs for which the hrness distriution follows power lw (25) the corresponding exponents γ s re found in the rnge of γ s = 2 4. For those distriutions with n exponentil decy the scle ˆr s (26) vries in the rnge ˆr s = 1.5 4. The power lws oserved for the ehvior of P (r, s) indicte certin level of orgniztion nd plnning which my e driven y the need to minimize the costs of infrstructure nd secondly y the fct tht points of interest tend to e clustered in certin loctions of city. Note tht this effect my e seen s result of the strong interdependence of the evolutions of oth the city nd its PTN. We wnt to emphsize tht the hrness effect is feture of the network given in terms of its routes ut it is invisile in ny of the grph representtions presented so fr. In prticulr PTN representtion in terms of simple grph which do not contin multiple links (such s L, P, C nd B-spces) cn not e used to extrct hrness ehvior. Furthermore, the multi-grph representtions (such s L, P, nd C -spces) would need to e extended to ccount for the continuity of routes. As noted ove, the notion of hrness my e useful lso for the description of other networks with similr properties. On the one hnd, the hrness distriution is closely relted to distriutions of flow nd lod on the network. On the other hnd, in the sitution of spce-consuming links (such s trcks, cles, neurons, pipes, vessels) the informtion out the hrness ehvior my e importnt with respect to the sptil optimiztion of networks. A generliztion my e redily formulted to ccount for rel-world networks in which links (such s cles) re orgnized in prllel over certin sptil distnce. While for the PTN this distnce is simply mesured y the length of sequence of sttions, more generl mesure would e the length of the contour long which these links proceed in prllel. V. GENERALIZED ASSORTATIVITIES To descrie correltions etween the properties of neighoring nodes in network the notion of ssorttivity ws introduced mesuring the correltion etween the node degrees of neighoring nodes in terms of the men Person correltion coefficient [50, 51]. Here, we propose to generlize this concept to lso mesure corre-

12 ltions etween the vlues of other node chrcteristics (other oservles). For ny link i let X i nd Y i e the vlues of the oservle t the two nodes connected y this link. Then the correltion coefficient is given y: M 1 i r = X iy i [M 1 i 1 2 (X i + Y i )] 2 M 1 i 1 2 (X2 i + Y i 2) [M 1 i 1 2 (X i + Y i )] 2 (27) where summtion is performed with respect to the M links of the network. Tking X i nd Y i to e the node degrees Eq. (27) is equivlent to the usul formul for the ssorttivity of network [50]. Here we will cll this specil cse the degree ssorttivity r (1). In the following we will investigte correltions etween other network chrcteristics such s the oservles considered ove, z 2, C i (6), C c (20), C g (21), C s (22), C (23). Consequently, this results in generlized ssorttivities of next nerest neighors (r (2) ), clustering coefficients (r cl ), closeness (r c ), grph (r g ), stress (r s ), nd etweenness (r ) centrlities. The numericl vlues of the ove introduced ssorttivities r (1) nd r (2) for the PTN under discussion re listed in Tle IV in L, P nd C-spces. With respect to the vlues of the stndrd node degree ssorttivity r (1) L in L-spce, we find two groups of cities. The first is chrcterized y vlues r (1) L = 0.1 0.3. Although these vlues re still smll they signl finite preference for ssorttive mixing. Tht is, links tend to connect nodes of similr degree. In the second group of cities these vlues re very smll r (1) L = 0.02 0.08 showing no preference in linkge etween nodes with respect to node degrees. PTNs of oth lrge nd medium sizes re present in ech of the groups. This indictes the sence of correltions etween network size nd degree ssorttivity r (1) L in L-spce. Mesuring the sme quntity in the P nd C-spces, we oserve different ehvior. In P-spce lmost ll cities re chrcterized y very smll (positive or negtive) vlues of r (1) P with the exception of the PTNs of Istnul (r (1) P = 0.12) nd Los Angeles (r (1) P = 0.12). On the contrry, in C-spce PTNs demonstrte cler ssorttive mixing with r (1) C = 0.1 0.5. An exception is the PTN of Pris with r (1) C = 0.06. As we hve seen ove, the PTNs demonstrte ssorttive (r (1) > 0) or neutrl (r (1) 0) mixing with respect to the node degree (first nerest neighors numer) k. Clculting ssorttivity with respect to the second next nerest neighor numer r (2) we explore the correltion of wider environment of djcent nodes. Due to the fct tht in this cse the two connected nodes shre t lest prt of this environment (the first nerest neighors of node form prt of the second nerest neighors of the djcent node) one my expect the ssorttivity r (2) to e non-negtive. Results for r (2) shown in Tle IV pper to confirm this ssumption. In ll the spces considered, we find tht ll PTNs tht elong to the group of neutrl mixing with respect to k lso elong to the sme City r (1) L r (2) L r (1) P r (2) P r (1) C r (2) C Berlin 0.158 0.616 0.065 0.441 0.086 0.318 Dlls 0.150 0.712 0.154 0.728 0.290 0.550 Düsseldorf 0.083 0.650 0.041 0.494 0.244 0.180 Hmurg 0.297 0.697 0.087 0.551 0.246 0.605 Hong Kong 0.205 0.632-0.067 0.238 0.131 0.087 Istnul 0.176 0.726-0.124 0.378 0.282 0.505 London 0.221 0.589 0.090 0.470 0.395 0.620 Los Angeles 0.240 0.728 0.124 0.500 0.465 0.753 Moscow 0.002 0.312-0.041 0.296 0.208 0.011 Pris 0.064 0.344-0.010 0.258 0.060-0.008 Rome 0.237 0.719 0.044 0.525 0.384 0.619 Sõ Polo -0.018 0.437-0.047 0.266 0.211 0.418 Sydney 0.154 0.642 0.077 0.608 0.458 0.424 Tipei 0.270 0.721 0.009 0.328 0.100 0.041 TABLE IV: Nerest neighors nd next nerest neighors ssorttivities r (1) nd r (2) in different spces for the whole PTN. group with respect to the second nerest neighors. For those PTNs tht disply significnt nerest neighors ssorttivity r (1) we find tht the second nerest neighor ssorttivity r (2) is in generl even stronger in line with the ove resoning. Recll tht oth closeness nd grph centrlities C c nd C g re mesured in terms of pth lengths, Eqs. (20), (21). It is nturl to expect tht djcent nodes will hve very similr (or lmost identicl) centrlities C c nd C g. In turn this will led to strong ssorttive mixing with high ssorttivities r c nd r g. This ssumption holds only if the verge pth length in the network is sufficiently lrge. The ltter is certinly the cse for PTNs in L-spce ut it does not hold in P nd even less in C- spces. Indeed, in L-spce, where most PTNs disply men pth length l L > 10 (see Tle II) we find vlues of rl c in the rnge rc L = 0.904 0.998 (rg L = 0.914 0.999). Exceptions re the two PTNs of cities with the smllest men pths. These re Moscow (rl c = 0.865, rg L = 0.870) with l L = 7.0 nd Pris (rl c = 0.831, rg L = 0.800) with l L = 6.4. In P nd C-spces where the men pth lengths re much shorter (of the order of three in P nd of the order of two in C-spces) the one-step difference in pth length etween djcent nodes leds to much reduced ssorttive mixing. Numericlly this is reflected in much lower (however positive) vlues of corresponding ssorttivities for PTNs where l is especilly smll. Indeed, for ll PTNs tht disply in P-spce men pth length l P < 2.7 we find rp c < 0.5 (rg P < 0.4). At the sme time, PTNs with lrger l P my disply lrger ssorttivities even in P-spce. The extreme exmple is Los Angeles with l P = 4.3 nd rp c = 0.914, rg P = 0.844. In C- spce, where vertices re routes the men pth length is even smller nd further reduction of closeness nd grph

13 centrlity ssorttivities is oserved. For five PTNs we find in C-spce l C < 2 (see Tle II) nd for these rc c < 0.3, rg C < 0.3. Agin the lrgest vlues re ttined in the Los Angeles PTN with l C = 3.4 nd rc c = 0.828, r g C = 0.648. For the other generlized ssorttivities (stress nd etweenness centrlity ssorttivities r s nd r nd clustering coefficient ssorttivity r cl we in generl find no evidence for ny (positive or negtive) correltion in ny of the spces considered. The only exception re the stress nd etweenness centrlity ssorttivities in L-spce, rl s nd rl. There, smll ut significntly positive vlues of rl s nd r L re found. The ltter is explined y the reltively lrge men pth length in this spce in conjunction with reltively smll node degree vlues. Let us recll tht stress nd etweenness centrlities essentilly count the numer of shortest pths medited y given node. If selected node is prt of mny such long pths while hving low degree, there is high proility tht ny of its neighors will lso e prt of these pths. Consequently, positive vlue of r c (r g ) will rise. The nlogous conclusion cn e drwn for nodes with low etweenness (or stress) centrlities. For most PTNs the vlues of the ssorttivities under considertion chnge in the rnge rl s = 0.26 0.64, r L = 0.20 0.61. Exceptions re the PTNs which in L-spce hve men pth length l L < 10, nmely Moscow, Pris nd Sõ Polo. There we find rl s = 0.02 0.10, r L = 0.02 0.10. VI. MODELING PTNS A. Motivtion nd description of the model Hving t hnd the ove descried welth of empiricl dt nd nlysis with respect to typicl scenrios found in vriety of rel-world PTNs we feel in the position to propose model tht leit eing simple my cpture the chrcteristic fetures of these networks. Nonetheless it should e cple of discriminting etween some of the vrious scenrios oserved. If we were only to reproduce the degree distriution of the network, stndrd models such s rndom networks [4, 52] or preferentil ttchment type models [6, 34, 53, 54, 55, 56] would suffice. The evolution of such networks however is sed on the ttchment of nodes. For description of PTNs the concept of routes s finite sequences of sttions is essentil [5, 23, 25, 28] nd llows for the representtion with respect to the spces defined ove. Moreover, tking route s the essentil element of PTN growth llows to ccount for the essentil iprtite structure of this network [20, 24, 26, 57]. Therefore, the growth dynmics in terms of routes will e centrl ingredient of our model. Another ovious requirement is the emedding of this model in two-dimensionl spce. To simplify mtters we will restrict the model to twodimensionl grid, in prticulr to squre lttice. Both the oservtions of power lw degree distriutions s well s the occurrence of the corresponding hrness distriutions descried ove indicte preference of routes to service common sttions (i.e. n ttrction etween routes). Let us descrie our model in more detil. As noticed ove, route will e modeled s sequence of sttions tht re djcent nodes on two-dimensionl squre lttice. Noting tht in generl loops in PTN routes re lmost sent, most simple choice to model PTN route is self-voiding wlk (SAW). It my sound less ovious tht route prt from eing non self-intersecting proceeds rndomly. However, the nlysis of geogrphicl dt [25] hs shown tht the frctl dimension of PTN routes closely coincides with tht of two-dimensionl SAW, d f = 4/3 [58]. To incorporte ll the ove fetures the model is set up s follows. A model PTN consists of R routes of S sttions ech constructed on possily periodic X X squre lttice. The dynmics of the route genertion dheres to the following rules: 1. Construct the first route s SAW of S lttice sites. 2. Construct the R 1 susequent routes s SAWs with the following preferentil ttchment rules: ) choose terminl sttion t x 0 with proility p k x0 + /X 2 ; (28) ) choose ny susequent sttion x of the route with proility p k x +. (29) In (28), (29) k x is the numer of times the lttice site x hs een visited efore (the numer of routes tht pss through x). Note tht to ensure the SAW property ny route tht intersects itself is discrded nd its construction is restrted with step 2). B. Glol topology of model PTN Let us first investigte the glol topology of this model s function of its prmeters. We first fix oth the numer of routes R nd the numer of sttions S per route s well s the size of the lttice X. This leves us with essentilly two prmeters nd, Eqs. (28), (29). Dependencies on R, S, nd X will e studied elow. For the rel-world PTNs s studied in the previous sections, lmost ll sttions elong to single component, GCC, with the possile exception of very smll numer of routes. Within the network however we often oserve wht ove we clled the hrness effect of severl routes proceeding in prllel for sequence of sttions. Let us first investigte from glol point of view which prmeters nd reproduce relistic mps of PTNs. In Fig. 18 we show simulted PTNs on lttices 300 300 for R = 1024, S = 64 nd different vlues of the prmeters

14 = 0.1 = 0.2 = 0.5 = 15 = 20 = 500 FIG. 18: PTN mps of different simulted cities of size 300 300 with R = 1024 routes of S = 64 sttions ech (color online). First row: = 0, = 0.1 0.5. Second row: = 0.5, = 15 500. With n increse of routes cover more nd more re. Increse of leds to clusteristion of the network. R S k L z L l mx L l L C L k P z P l mx P l P C P c P k C z C l mx C l C C C c C 256 16 0.5 2.92 1.66 61 20.8 4.7 10 3 44.15 3.18 7 3.0 4.7 10 2 7.98 86.39 1.36 6 1.9 1.2 10 2 2.22 256 16 5.0 2.99 1.74 80 21.7 7.5 10 3 42.95 3.76 9 3.4 8.8 10 2 11.7 59.96 1.99 8 2.2 1.5 10 2 2.79 256 32 0.5 2.76 1.60 127 38.1 3.0 10 4 84.45 4.32 8 3.3 1.9 10 3 13.6 60.51 1.75 7 2.2 1.6 10 2 2.90 256 32 5.0 2.90 1.72 177 43.1 5.3 10 4 74.24 5.22 10 4.0 3.8 10 3 23.7 33.06 2.69 9 2.8 2.3 10 2 4.55 512 16 0.5 2.95 1.68 73 22.5 6.7 10 3 50.07 3.39 7 3.1 6.5 10 2 9.14 169.7 1.44 6 1.9 2.3 10 2 2.25 512 16 5.0 3.12 1.78 80 23.3 1.0 10 4 51.56 3.79 10 3.5 1.2 10 3 12.3 115.3 2.24 9 2.1 2.9 10 2 2.88 512 32 0.5 2.83 1.63 166 44.2 4.7 10 4 99.53 4.56 10 3.6 2.8 10 3 15.7 118.4 2.03 9 2.2 3.0 10 2 2.92 512 32 5.0 3.12 1.79 175 44.6 7.2 10 4 97.05 5.37 9 3.9 4.7 10 3 22.2 60.36 3.08 8 2.7 4.4 10 2 5.04 1024 64 0.5 2.86 1.66 325 80.7 3.3 10 5 242.2 6.32 9 3.7 1.1 10 4 23.4 213.3 2.42 8 2.2 6.1 10 2 3.10 1024 64 1.0 2.97 1.72 355 88.5 4.8 10 5 222.2 6.74 12 4.2 1.7 10 4 32.4 143.9 2.97 11 2.5 7.9 10 2 4.39 TABLE V: Chrcteristics of the simulted PTN with X = 300, = 0 for different prmeters R, S, nd. nottions s in Tle II. The rest of nd. Ech route is represented y continuous line trcing the pth long its sequence of sttions. For representtion purposes, prllel routes re shown slightly shifted. Thus, the line thickness nd intensity of colors indicte the density of the routes. Prmeter governs the evolution of ech single susequent route. If = 0 ech susequent route is restricted to follow the previous one. The chnge of simulted PTNs with for fixed = 0 is shown in the first row of Fig. 18. For smll vlues of = 0 0.1 the PTNs otined result in lmost ll routes following the sme pth with only few devitions. Incresing from = 0.1 to = 0.2 the re covered y the routes increses while the mjority of the routes re concentrted on smll numer of pths. Further incresing to = 0.5 nd eyond we find wider distriuted coverge with the centrl prt of the network remining the most densely covered re. This is due to the non-equilirium growth process descried y Eqs. (28), (29). The prmeter quntifies the possiility to strt new route outside the existing network. For vnishing = 0 the resulting network lwys consists of single connected component, while for finite vlues of few or mny disconnected components my occur. The results found for = 0 nd vrying prmeters re completely independent of the lttice size X provided X is