Research Article  Open Access
A Century of Topological Coevolution of Complex Infrastructure Networks in an Alpine City
Abstract
In this paper, we used complex network analysis approaches to investigate topological coevolution over a century for three different urban infrastructure networks. We applied network analyses to a unique timestamped network data set of an Alpine case study, representing the historical development of the town and its infrastructure over the past 108 years. The analyzed infrastructure includes the water distribution network (WDN), the urban drainage network (UDN), and the road network (RN). We use the dual representation of the network by using the Hierarchical Intersection Continuity Negotiation (HICN) approach, with pipes or roads as nodes and their intersections as edges. The functional topologies of the networks are analyzed based on the dual graphs, providing insights beyond a conventional graph (primal mapping) analysis. We observe that the RN, WDN, and UDN all exhibit heavy tailed node degree distributions with high dispersion around the mean. In 50 percent of the investigated networks, can be approximated with truncated [Pareto] powerlaw functions, as they are known for scalefree networks. Structural differences between the three evolving network types resulting from different functionalities and system states are reflected in the and other complex network metrics. Smallworld tendencies are identified by comparing the networks with their random and regular lattice network equivalents. Furthermore, we show the remapping of the dual network characteristics to the spatial map and the identification of criticalities among different network types through colocation analysis and discuss possibilities for further applications.
1. Introduction
Many complex systems can be described as networks [1], and with recent increases in computing power it is now feasible to investigate the topologies of entire networks consisting of highresolution data [2]. Examples of these types of investigations range from molecular interaction networks (e.g., protein interactions of cells) and social networks (e.g., communication between humans) to global transportation systems and individual human mobility [3–6].
Despite the differences in various types and representations of these networks, important commonalities exist. The analysis of complex networks gives insight to structural morphologies, similarities, recurring patterns, and scaling laws [7, 8]. The applications are multifaceted: identification of central nodes; prediction of future developments and network growth; information transfer; identification of vulnerabilities to enhance security [9]; and improvement of network resilience [10, 11]. Complex network analyses of critical infrastructure, such as water distribution networks (WDNs) and urban drainage networks (UDNs), provide valuable insights beyond the traditional engineering approaches, to design and operate systems in a more reliable way and to help buildup structural resiliency [12, 13].
In the past, most structural features in complex networks were investigated based on a conventional graph representation (socalled “primal space”), where pipes or conduits are the edges and their intersections the vertices of a mathematical graph [14, 15]. Conversely, different approaches, based for example on common attribute classification (i.e., road name or pipe size) or intersection continuity (i.e., maximum angle of deflection), consider the network structure in its “dual space”, i.e., functional components (e.g., pipes with same diameter) which belong together, represent the vertices and their intersection the edges of the graph [16, 17]. Further explanations are provided in the next section. Unlike the conventional primal representation, dual mapping approaches may also consider the continuity of links (pipes or conduits) over a variety of edges and hierarchy (e.g., pipe diameter; isolation valves; maximum designed flow; speed limits; road class) for further graph analysis.
There exist different ways of creating the dual graph of a network, taking into account physical (e.g., geometric) and/or behavioral (e.g., symbolic) considerations. The street name (SN) approach, for example, uses the historical naming conventions to create the dual graph, but neglects the geometrical properties of the network. Hybrid approaches, like the Hierarchical Intersection Continuity Negotiation (HICN) [16], combine geometric (e.g., maximum angle of deflection of connected roads) and hierarchical (e.g., road class) attributes, to better capture the structural network topology resulting from topdown (centralized designs) and bottomup localplanning actions (selforganization).
Previous studies using the dual mapping approach were mainly performed on road networks (RNs) [16–18], but some also on water distribution and urban drainage networks [19–21]. In principle, an extension to each network type is possible. Masucci et al. [16] investigated the road network growth for the city of London and found stable statistical properties to describe the topological network dynamics. Krueger et al. [20] applied the HICN principle for the first time to the evolving sewer networks in a large Asian city with 4 million people. The authors found that sewer network types quickly evolve to become scalefree in space and time. In Jun and Loganathan [19] a dual mapping approach was used to describe the connectivity of isolation zones in water distribution networks.
Klinkhamer et al. [21] examined the colocation of existing road and sewer networks in a large Midwestern US city and homoscedasticity of subnets across the city but did not examine temporal evolution of these networks. In Mair et al. [22] the geospatial colocation of roads, pipes, and sewers was investigated using data set for three Alpine case studies, finding strong similarities between these networks. Studies on the coevolution of water infrastructure networks (water distribution and urban drainage) and road network are crucial when investigating functional interdependencies and cascading vulnerabilities across multiplex network layers. Examples are the floodinduced change in road traffic or the collapse of entire road segments causing flow disruptions in all networks to different extents.
In this paper, we present for the first time a topological analysis of three infrastructure networks coevolving over a century. The results of the dual mapping for a unique dataset of 11 timestamped water distribution and urban drainage network states and 8 timestamped road networks of the mediumsize Alpine case study city, as the town and its infrastructure, evolved during the past 108 years, and the population tripled from about 40,000 to about 130,000. First results of this case study are presented in Zischg [23]. We investigated network topological metrics using the HICN dual mapping approach [16]. We observe that some infrastructure networks show node degree distributions that behave like truncated powerlaws under the dual representation. However, this “scalefree” network characteristics depend on the network type and change over time. With the presented methodology, differences and similarities of patterns (e.g., vertex connectivity) and trends for the infrastructure development are obtained. This study includes an investigation of the sensitivity of the dual mapping approach, using different criteria to build the new graph. The reflected structural features, such as the backbone of the networks, were uncovered for each network type and remapped to the spatial map. A further analysis shows the pairwise colocation of high node degree components (“network hubs”) across different infrastructure network types, which builds the basis for analyzing disturbances and structural resilience.
2. Data Analyses
2.1. Network Connectivity
Node degree distribution is a significant topological property of complex networks. The degree of a node in an undirected network describes the number intersecting links and is calculated through the network’s adjacency matrix , where the degree of node is defined by the sum of the th row of . For example, the node degree in social networks represents the number of contacts. Scalefree networks show node degree distributions that follow a Pareto powerlaw distribution [20, 21], with for , whereas random networks have Poisson distributed node degrees. We use the method proposed by Clauset et al. [24] to test the powerlaw hypothesis and determine scaling parameters of the node degree distributions for the various network states. By calculating the p value, an indicator for the goodnessoffit is determined. In case the p value is greater than 0.1, the powerlaw is a plausible hypothesis for the data within given ranges. However, the definitive recipe to fit powerlaw distributions does not yet exist [25]. The mean node degree for undirected networks is defined as , where is the total number of edges and is the total number of vertices. In the limits a mean node degree of 2 indicates a treelike network structure, and grid patterns or cyclic structures have mean node degrees around 4 [26]. Higher statistical moments of are also important, including the variance that reflects the dispersion around the mean [27].
Along with the node degree distribution, the characteristic path length (or average path length) is an important and robust measure of network topology. It quantifies the level of integration/segregation throughout the network. In water infrastructure and power grid networks energy losses are dependent on the characteristic path length. It is calculated by the average shortest path distance between all couples of nodes as follows:where is the number of vertices and denotes the shortest path between vertex and [3]. The probability density function of the shortest path lengths (between all couples of nodes for RNs; between all terminal nodes and the source or sink node, for WDNs and UDNs, respectively), can, for example, be considered as the approximation of the traveltime distribution with a nearly consistent distribution of flow velocities.
The local clustering coefficient of node describes the connectivity (number of edges m) among its k neighbors. A perfect cluster/clique indicates a full connection of all nodes/individuals. If an isolation of one node in the cluster occurs, the other nodes remain connected. Conversely, a indicates that node holds together all its neighboring nodes. Regarding infrastructures, a higher clustering coefficient indicates the existence of local and alternative flow paths in the network. In undirected graphs, for an individual node is defined as follows: The overall (average) clustering coefficient can be determined by averaging the local clustering across all nodes. When the network reveals a high , which is typical for regular lattice networks (high local efficiency), and has a small , as found in random networks (high global efficiency), then the network can be characterized as smallworld network [30]. Telesford et al. [31] developed methodology to test the “small worldness” by comparing the network with its equivalent random and regular lattice networks, which have the same node degree distribution, as follows:where is the characteristic path length of the randomized network equivalent and is the clustering coefficient of the regularized network equivalent. A value of 0 indicates a network that is in perfect balance between normalized values of high clustering and low characteristic path length. Negative values indicate a graph with more regular characteristics, whereas positive values indicate more random graph characteristics [31]. Smallworld networks are significantly more clustered (segregated) than random networks and have the same characteristic path length as random (integrated) networks, making them locally and globally efficient, for example, for optimal information transfer [30].
2.2. HICN Principle for Dual Mapping
The HICN approach emphasizes the functional topology of the network by aggregating components (e.g., pipes, conduits, and roads) with identical attributes (e.g., pipe diameter, pipe segments, and road type), while also maintaining a certain level of straightness (e.g., road sections) [16]. After reducing the network complexity with this “generalization model,” the aggregated edges are converted into vertices and the intersections are converted into edges. The resulting graph is the socalled “dual” (mapped) representation of the “primal” graph (see Figure 1). In addition to the edge class, the angular threshold is a second criterion used for the generalization model. It defines the maximum exterior convex angle of connected edges being merged [17].
The HICN allows for reducing the network complexity of the primal map and considers the hierarchy of network elements (e.g., different level of detail of the pipe representation). Identical cohorts of edges are considered as a single component, the dual node. However, the dual mapped network circumvents this issue through generalization and still preserves the connectivity information of the original network. Another advantage is the detachment from the geographical embedding, allowing the network to be nonplanar. With this methodology, the underlying hierarchy (e.g., highly connected components) of the network can be uncovered.
3. Alpine Case Study
To investigate the coevolving topology of three urban infrastructure networks with the HICN dual mapping methodology, we utilize available, highresolution network data for a mediumsize Alpine city. The temporal evolution of the urban infrastructure networks is defined through timestamped system states at 10year intervals, starting with the year 1910 for water distribution and urban drainage networks. The road network data starts with the year 1940, since historical orthophotos to reconstruct the network were only available from that time. The city has grown from approximately 40,000 inhabitants in 1910 to 130,894 in 2016. The historical data set describes the expansion of the networks and includes pipe rehabilitation, changing source (e.g., reservoirs), and sink nodes (e.g., sewer outfalls), altering population densities and variations of water consumption patterns of the water distribution and urban drainage systems. The detailed description of the network reconstruction for this case study can be found in the works of Sitzenfrei et al. [32] and Glöckner [33].
Figure 2 shows the timestamped networks at selected stages in the primal map. For the UDN a greater thickness and a darker color of the edges indicate a larger conduit size, which connect to one large (biological) waste water treatment plant (WWTP) in the eastern part of the city after the 1970s. Before the wastewater was discharged to the river after mechanical treatment (sedimentation, rack), outside the urban areas at the eastern parts of the city. When actually building the WWTP, there were no major changes in the combined drainage network necessary (combined transportation of sewage and storm water) and the old outlets were transferred to combined sewer overflows (CSOs). For security reasons the pipe diameter for the WDN cannot be shown. The color shading at the road network in 2010 indicates the different road types, ranging from residential (light grey), tertiary, secondary, and primary roads to motorways (black). Furthermore, the water demand which directly influences the flows of the WDN and the UDN during dry weather periods is illustrated. Clearly noticeable is the general reduction of the average water demand after the 1970s mainly due to water saving measures.
Table 1 gives a short narrative introduction on the history of the infrastructure networks of the investigated case study during the last century.

4. Results and Discussion
In this section, we present the results of the historical coevolution in the dual representation of the three infrastructure networks (WDN, UDN, and RN). In a sensitivity analysis of the HICN method, we determine the effects on based on the variations of the angular threshold, the edge class, and the network partition (entire network vs. largest connected component (LCC)), before we present the development of the network characteristics over time. Finally, we also show the remapping of dual network characteristics to the primal map, investigate correlations between dual node degree and edge class (e.g., road type), and identify pairwise spatial colocations of dual nodes among the three network types.
4.1. Dual Mapping
The application of the HICN dual mapping to the historical infrastructure networks for the first (year 1910 for WDN and UDN, and 1940 for RN) and last stages (year 2010) is illustrated in Figure 3. The dual graphs show the node degree (darker and larger nodes represent central network “hubs” with high node degree). For the application of the HICN method to the water infrastructure networks, we used identical pipe diameters (for WDN) and conduit sizes (for UDN) for the edge class criterion, in combination with an angular threshold of 180 degrees; i.e., we ignore the curvature of pipe/conduit segments. On the other hand, for the dual mapping of the road networks, we used the of 45 degrees [35] and the road type (residential road to motorway) as criteria for the generalization model. We discuss the application of different parameter values in the next section.
(a)
(b)
(c)
4.2. Sensitivity of the Dual Mapping
For the HICN approach different angular thresholds and edge attributes can be considered to construct the dual graph. In Figure 4 we show the results of the sensitivity analysis on the node degree distribution for three parameters used for the HICN approach. First, we identify the sensitivity of the angular threshold using the WDN of 2000 as reference. Second, we choose the UDN 1950 with 10 subnetworks and investigate the graph characteristics for the entire network and the largest connected component (LCC) only. Finally, we show the effect of neglecting the edge class criterion (road type) for the RN in 2010 as reference.
(a)
(b)
(c)
(d)
(e)
(f)
With increasing angular threshold (see Figure 1) and the associated less strict criterion for aggregation of edges, we observe that more pipes (of the same class) are merged together and therefore the size of the dual graph is reduced (see Figure 4(a)). For example, the reference network has 7,827 edges in the primal space, which are generalized to 3,050 and 1,764 dual nodes for angular thresholds of 15 and 180 degrees, respectively. Figure 4(d) shows the outcome of 7 variations of the angular threshold from 15 to 180 degrees. All resulting node degree distributions are heavy tailed and by visual inspection relatively similar. However, when applying the method to test the powerlaw hypothesis as proposed by Clauset et al. [24], only for degrees Pareto distributions , are plausible data fits (p value ≥ 0.1). At this point it should be mentioned that powerlaw fitting is still a controversial issue and a definitive recipe to fit powerlaw distributions and to distinguish between powerlaw and powerlawlike distributions does not yet exist [25]. A minimum node degree of 3 for curve fitting was determined to be the most likely one, for all three network types. For higher threshold angles slightly decreasing slopes are observed. The results using angular thresholds of 90 and 180 degrees are identical, meaning that no sharp inner angles between connected pipes of the same class are found in the graph. Unlike for road networks, where also is a criterion for visibility and navigation, we suggest that restricting in WDNs and UDNs is less important to find and aggregate pipe segments with unique identity because of their underground locations. We conclude that for investigating the node degree distribution of the historical networks, the angular threshold is of minor importance, but must be consistently applied between the types and the states of the networks.
Figure 4(b) shows the results when investigating the entire graph (all subnetworks) and the largest connected component (LCC), which contains 87% of the total nodes. For the reference network UDN 1950, no significant changes in the Pareto exponent are seen (see Figure 4(e)). Finally, the neglect of the road class as additional criterion for the HICN method slightly decreases the size of the network but has a low effect on the probability density function (see Figures 4(c) and 4(f)). This can be interpreted with the low incidence of changing road types across straight road segments.
In this study, we used the pipe diameter to identify functionally identical pipe segments. Other examples could be the age or the material of the pipe, but also classified (design) flows or (measured or modeled) head losses might be used to construct the dual graph. The advantage of using the pipe diameter is that it is a surrogate measure of the flow and is independent of a hydraulic simulation. The black dots in Figure 4 represent the network configuration as presented in the subsequent section.
4.3. Characteristics of the Dual Graph
Resulting node degree distributions, , for the WDN, UDN, and RN in the dual representation are presented in Figure 5, plotted for all network types on loglog axes for the timestamped states. A darker node color indicates a younger and more mature network state. We observe heavytailed node degree distributions for 30 network states with dispersion indices , indicating high variance around the mean. Detailed analyses showed that 50% of the timestamped networks follow truncated powerlaw [Pareto] distributions across all network types. This does not necessarily mean, however, that alternative (e.g., lognormal, exponential) distributions are more plausible. The scaling parameters γ fall in the range between 3.23 and 3.50 for WDN ( = 2 or 3), 3.25 and 3.76 for UDN ( = 3 or 5), and 2.67 and 2.78 for RN ( = 3). For the WDN (mean: 3.35 ± 0.09) and UDN (mean: 3.49 ± 0.12) larger slopes γ and a stronger decrease of the number of leaf nodes ( = 1) are found compared with the RN (mean: 2.70 ± 0.04). A possible explanation of the truncation could be the missing house (low degree) connections for both water infrastructure networks. For the detailed parameters and the statistical tests, we refer to Appendix A.
(a)
(b)
(c)
The truncated powerlaw [Pareto] distribution also indicates that the probability of finding nodes with many connecting links (“hubs”) is much lower than of nodes with few connections (terminal dual nodes). According to the literature this behavior is typical for scalefree networks, which are dominant in most natural networks. The extent to which these distributions fit a powerlaw can be a useful marker of network resilience [36].
The “scalefree” similar characteristics, within the observed range , are also indicated with the significant higher maximum degrees (representing a “network hub”) compared to the mean degree (see Appendix B). Although the mean degree of the UDN is larger than that for the WDN, the maximum degree of both water infrastructure networks is similar. One reason for that is fewer changes in the conduit diameters, resulting in the aggregation of more conduits and thus having higher connectivity. The highest connectivity and dispersion index is found for the road network, indicated by high and high (see Appendix B and Figure 6(c)).
(a)
(b)
(c)
(d)
(e)
(f)
The growth of the networks in terms of total number of dual nodes over time is illustrated in Figure 6(a). Highest growth rates of the networks are seen in the 1960s and 1970s, which can be partly related to the economy boom and the implementation of the waste water treatment plant [32]. Figure 6(b) presents the scaling parameter of the node degree distributions over time. A change in the lower bound is indicated in round brackets. Although the powerlaw hypothesis could only be clearly proven in 50% of the timestamped networks (see Appendix A), we show the scaling parameter for all network states for comparison purposes. The statistical analysis shows that the RN tends towards a clear scalefree behavior during evolution, whereas the opposite is observed for the UDN. No clear trends were identified for the WDN (see p values ≥ 0.1 in Appendix A). According to Achard et al. [36] the extent to which these distributions fit a powerlaw can be a useful indicator of network resilience. Furthermore, when comparing γ with previous studies in the literature, similar ranges between 2 and 4 are reported [17, 20, 21].
In this study, the lowest exponents are reported for the RN compared to those for the WDN and UDN. Peaks of γ for the WDN and UDN can be explained with the treelike expansion of the network to new parts of the city and without a strong network densification at those times (see description in Table 1). During the last part of the century the for the WDN tends to decrease, whereas the UDN and RN evolution is characterized by a network growth with a nearly constant powerlaw exponent . This could indicate that the networks are now topologically “mature”; i.e., a similar behavior is expected when the network grows further. The ratios of variance and the mean of the node degrees are illustrated in Figure 6(c). For WDN and UDN it is nearly constant during their evolution [1 ≤ ≤ 2]. An increasing bilinear trend of dispersion is observed for the RN, indicating the preferential attachment tendencies of new nodes to already well connected nodes. All networks show dispersion indices greater than 1 and thus exceed the value of expected random graphs following a Poisson distribution (D = 1).
As a measure of functional segregation (local efficiency), high average clustering coefficients describe the presents of cliques, where neighboring nodes are well connected among each other. During the last half of the century, C remained nearly constant for all infrastructure networks, with highest values for the RN, followed by UDN and WDN at year 2010.
The characteristic path length of the network states in dual representation is shown in Figure 6(e). In the dual representation, the path length defines the number of changing edge attributes (e.g., diameter changes) between two dual vertices. Small indicate a global efficiency meaning that every vertex is connected to every other through a short distance. In general, increases with expanding geographical boundaries of the network (Figure 2). During the 1980s there is a significant increase of for the UDN, possibly because of the treelike connections of peripheral zones and neighboring villages to the central wastewater treatment plant. In contrast, for the WDN decreases during the past few years, as a result of the WDN densification and the construction of alternative flow paths for redundancy purposes. Furthermore, it is remarkable that for the UDNs at the early stages of the century is lower compared to that for WDNs. One reason for that is that the UDN at the historical center of the Alpine city had fewer alterations of conduit diameters (see Figure 2, middle left), resulting from the coarse design concepts and material limitations at that time. As the road network grows over time, we observe the lowest values for are nearly constant trend over time, indicating a high network integrity.
Figure 6(f) provides a quantitative measure of the smallworld properties over a spectrum of network topologies. Telesford et al. [31] describe in proximity to zero as smallworld network; however, no sharp boundary of the smallworld region exists. Positive indicate a graph with more random characteristics, as we observe for all infrastructure during the first half of the century. However, over time, the networks tend to become less random and closer to small worldness. While is always larger than 0.40 for the WDN, UDN and RN have  values less than 0.16 after the first network state. Therefore, we claim the smallworld property of UDN and RN, while the former in recent decades can be considered as the closest to the smallworld optimum.
Figure 7 shows an example of the random and regular lattice equivalents for the WDN 1910, needed to determine the small worldness. By preserving the same node degree distribution, edges of the initial network (see Figure 7(c)) are rewired in a stochastic process [37, 38]. Differences between the networks in Figures 7(a) and 7(b), as well as in Figures 7(d) and 7(e), are the preservation of the network connectivity. The effective number of rewiring per edge for creating the network equivalents is around 10 (for the regularization approximately 100 times more iterations needed), determined though parameter convergence (see Appendix C).
Previous studies have shown that understanding the network structure gives insights to vulnerabilities and structural resilience of the systems (i.e., scalefree networks are found to be highly resilient against random failure but vulnerable to targeted attacks) [9, 39]. Implications to the infrastructure management could, for example, relate to system operation to emphasize special protection or increased maintenance of critical network components (“hubs”). While this is relatively obvious at key points (such as water or wastewater treatment plants), in most cases less attention has been paid to individual pipe sections due to their complexity. A first step towards more resilient networks is already proposed by Mair et al. [40] and Zischg et al. [41], to use a “lesscritical” subset of the colocated road network with strong similarities to generate possible WDNs. Second, when considering all infrastructure networks as an entire system (“multiplex network”), severe failure cascades through the network layers should be prevented through avoiding certain interlayer links (e.g., colocation of “hubs”) [42]. Resulting resilient infrastructure systems should be capable of minimizing the failure impacts and a fast recovery to a stable system state [11]. While the former factor certainly depends on the (multiplex) network topology and on the type of network disruption (see e.g., [43, 44]), for the latter one several influencing factors exist, which in most cases are more difficult to quantify (e.g., resources availability, societal needs, preparedness actions, economy adjustments, etc.) [6]. However, the assessment of infrastructure system resilience goes beyond the scope of this study.
4.4. Remapping from Dual to Primal Space
Remapping the dual graph characteristics to the primal graph representation allows for georeferenced visualization and further spatial analyses among the different network types. Figure 8(a) shows the example of remapping the dual (node) degree to the original physical embedded infrastructure networks at year 2010. The highlighted bold edges represent the dual network “hubs”, i.e., the elements with the most interconnections. In case those elements fail, the network gets easily fragmented, making for example the emergency supply very difficult through the central WDN. A comparison of the dual degree and the edge class (pipe diameter, conduit size, and road type) shows that they do not necessarily correlate. This implies that edges with the highest capacity (large diameter in WDN, large conduits sizes in UDN, and motorways (mot) in the RN) do not have the most intersections.
(a)
(b)
(c)
For the RN the “high degree hubs” are identified to be tertiary and primary roads. The reason for this finding for the WDN and UDN is that the source and sink nodes (water sources and WWTP) are located outside the city, requiring “high capacity transmission edges” to the inner parts of the network, which usually have fewer connections. However, low degree “terminal dual nodes” are mostly found in the edge classes with the lowest capacity (see bivariate histograms in Figure 8(b)).
4.5. Geospatial Colocation
We present the results of the geospatial colocation analysis to identify the spatial relationships of dual node degrees among the three network types. Mair et al. [22] found that approximately 90% of the WDN and UDN are located below the RN for this case study, however without considering their topological characteristics. Here we pairwise compare the normalized dual (node) degrees of the colocated infrastructure networks (see Figure 8(c)). Nodes in the upper right corner of the scatter plots indicate that high node degree components of both networks are colocated. For example, one conduit segment (UDN) with the highest is colocated with one pipe segment (WDN) with the second highest (indicated with the arrow and a red circular marker). This could be an indicator for an increased cascading vulnerability when failures occur across multiple networks.
These findings provide a first step towards the assessment of structural resilience and network interdependencies. Besides the identification of the “connectors” (high connectivity nodes), the “carriers” (high capacity node) should be further addressed. One future direction of this work could be the analysis of cascading failure across multiple network layers (e.g., a pipe break occurs which affects the water supply but also influences the colocated road and urban drainage network parts due to traffic rerouting and additional inflow to sewers). These analyses could provide helpful insights in the resilient (re)design of networks, providing an integrated view across the usually separated systems.
5. Conclusions
The historical data set of water distribution, urban drainage, and road networks of an Alpine city was investigated using a dual mapping approach (Hierarchical Intersection Continuity Negotiation method) and complex network analysis metrics were estimated. From a sensitivity analysis for the angular threshold using the HICN approach it is concluded that, unlike for road networks, where also is a criterion for visibility, we suggest that restricting in WDNs and UDNs is less important than finding and aggregating pipe segments with unique identity because of their underground locations.
The node degree distributions are heavy tailed with greater variance than expected for random networks. Across all three network types and for 50% of the timestamped networks truncated [Pareto] powerlaw functions , ( between 2.67 and 3.76; mostly ranging between 2 and 3) were proven and were in the range of previous studies on urban infrastructure networks. While the RN tends towards a clear scalefree behavior during evolution, the opposite is observed for the UDN. No clear trends were identified for the WDN. All networks tend to fall in the smallworld region characterized by high global and local efficiencies, during their spatiotemporal evolution. We conclude that similar to other “selforganized” networks, infrastructure networks in their dual representation can also exhibit scalefree properties and thus exhibit the failure dynamics typical of scalefree networks.
Previous studies showed that networks with highly variable node degree distributions, such as scalefree networks, are highly resilient to random failures but have high vulnerabilities to targeted attacks. Furthermore, the extent to which these distributions fit a powerlaw can be a useful indicator of network resilience. With the presented dual mapping methodology, patterns (e.g., vertex connectivity) and trends for the future infrastructure development were obtained. Furthermore, the reflected structural features, such as the “highly connected” components of the networks, were identified and remapped to the primal map. It was shown that the “highly connected” components do not necessarily correlate with “high capacity” components. Through the pairwise comparison of the dual degree of different network types, we identified the location where “high degree” components are colocated. These findings could be a measure of assessing the interdependencies and cascading vulnerabilities across multilayer networks.
Appendix
A.
See Table 2.
 
recalculated to maximize p value from/to: = (5)/(3); (9)/(3); (4)/(3); p values ≥ 0.1 proving the powerlaw hypothesis are shown in bold (50% of the networks). 
B.
See Table 3.

C.
See Figure 9.
(a)
(b)
(c)
(d)
(e)
Data Availability
The current and historical georeferenced data of the water distribution and urban drainage networks, used to support the findings of this study, have not been made available because they are protected for security reasons. The road network data which was used to support the findings of this study was taken from open street maps. The dual mapped representation of the water distribution, urban drainage, and road networks, used to support the findings of this study, is available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was initiated at the international synthesis workshop “Dynamics of Structure and Functions of Complex Networks”, held at Korea University, Seoul, South Korea, during summer 2015. The authors gratefully acknowledge Sybille Glöckner, who provided the historical water distribution network, reconstructed during her master thesis. This research was partially funded by the Austrian Research Promotion Agency (FFG) within the research project ORONET [Project number: 858557], the Austrian Science Fund (FWF) [Project number: P 31104N29], the Austrian Marshall Plan Foundation, and the NSF Collaborative Research RIPS Type 2: Resilient Simulation for Water, Power, and Road Networks [Project number: 1441188].
References
 P. Blanchard and D. Volchenkov, Mathematical Analysis of Urban Spatial Networks, Understanding Complex Systems, 1, Springer, Berlin, Germany, 2009. View at: Publisher Site
 S. H. Strogatz, “Exploring complex networks,” Nature, vol. 410, no. 6825, pp. 268–276, 2001. View at: Publisher Site  Google Scholar
 Y. Assenov, F. Ramírez, S.E. Schelhorn, T. Lengauer, and M. Albrecht, “Computing topological parameters of biological networks,” Bioinformatics, vol. 24, no. 2, pp. 282–284, 2008. View at: Publisher Site  Google Scholar
 M. C. González, C. A. Hidalgo, and A.L. Barabási, “Understanding individual human mobility patterns,” Nature, vol. 453, no. 7196, pp. 779–782, 2008. View at: Publisher Site  Google Scholar
 S. Uhlmann, H. Mannsperger, J. D. Zhang et al., “Global microRNA level regulation of EGFRdriven cellcycle protein network in breast cancer,” Molecular Systems Biology, vol. 8, article 570, 2012. View at: Publisher Site  Google Scholar
 D. A. Eisenberg, J. Park, and T. P. Seager, “Sociotechnical Network Analysis for Power Grid Resilience in South Korea,” Complexity, vol. 2017, 14 pages, 2017. View at: Google Scholar
 A. Barabasi and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 A. P. Riascos, “Universal scaling of the distribution of land in urban areas,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 96, no. 3, Article ID 032302, 2017. View at: Google Scholar
 K. A. Zweig and K. Zimmermann, “Wanderer between the worlds  Selforganized network stability in attack and random failure scenarios,” in Proceedings of the 2nd IEEE International Conference on SelfAdaptive and SelfOrganizing Systems, SASO 2008, pp. 309–318, Italy, October 2008. View at: Google Scholar
 J. P. G. Sterbenz, E. K. Çetinkaya, M. A. Hameed, A. Jabbar, S. Qian, and J. P. Rohrer, “Evaluation of network resilience, survivability, and disruption tolerance: analysis, topology generation, simulation, and experimentation: Invited paper,” Telecommunication Systems, vol. 52, no. 2, pp. 705–736, 2013. View at: Publisher Site  Google Scholar
 H. Klammler, P. S. C. Rao, and K. Hatfield, “Modeling dynamic resilience in coupled technologicalsocial systems subjected to stochastic disturbance regimes,” Environment Systems and Decisions, vol. 38, no. 1, pp. 140–159, 2018. View at: Publisher Site  Google Scholar
 A. Yazdani, R. A. Otoo, and P. Jeffrey, “Resilience enhancing expansion strategies for water distribution systems: a network theory approach,” Environmental Modelling & Software, vol. 26, no. 12, pp. 1574–1582, 2011. View at: Publisher Site  Google Scholar
 A. Yazdani and P. Jeffrey, “Applying network theory to quantify the redundancy and structural robustness of water distribution systems,” Journal of Water Resources Planning and Management, vol. 138, no. 2, pp. 153–161, 2012. View at: Publisher Site  Google Scholar
 F. Zeng, X. Li, and K. Li, “Modeling complexity in engineered infrastructure system: Water distribution network as an example,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 27, no. 2, Article ID 023105, 2017. View at: Publisher Site  Google Scholar
 S. Porta, P. Crucitti, and V. Latora, “The network analysis of urban streets: a primal approach,” Environment and Planning B: Planning and Design, vol. 33, no. 5, pp. 705–725, 2006. View at: Publisher Site  Google Scholar
 A. P. Masucci, K. Stanilov, and M. Batty, “Exploring the evolution of London's street network in the information space: A dual approach,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 89, no. 1, Article ID 012805, 2014. View at: Google Scholar
 S. Porta, P. Crucitti, and V. Latora, “The network analysis of urban streets: A dual approach,” Physica A: Statistical Mechanics and its Applications, vol. 369, no. 2, pp. 853–866, 2006. View at: Publisher Site  Google Scholar
 M. Hu, R. Jiang, R. Wang, and Q. Wu, “Urban traffic simulated from the dual representation: Flow, crisis and congestion,” Physics Letters A, vol. 373, no. 2324, pp. 2007–2011, 2009. View at: Publisher Site  Google Scholar
 H. Jun and G. V. Loganathan, “Valvecontrolled segments in water distribution systems,” Journal of Water Resources Planning and Management, vol. 133, no. 2, pp. 145–155, 2007. View at: Publisher Site  Google Scholar
 E. Krueger, C. Klinkhamer, C. Urich, X. Zhan, and P. S. C. Rao, “Generic patterns in the evolution of urban water networks: Evidence from a large Asian city,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 95, no. 3, Article ID 032312, 2017. View at: Google Scholar
 C. Klinkhamer, E. Krueger, X. Zhan, F. Blumensaat, S. V. Ukkusuri, and P. S. C. Rao, Functionally Fractal Urban Networks: Geospatial Colocation and Homogeneity of Infrastructure, 2017, https://export.arxiv.org/abs/1712.03883v1.
 M. Mair, J. Zischg, W. Rauch, and R. Sitzenfrei, “Where to find water pipes and sewers?On the correlation of infrastructure networks in the urban environment,” Water (Switzerland), vol. 9, no. 2, p. 146, 2017. View at: Google Scholar
 J. Zischg, Understanding patterns, dependencies and resilience in complex urban water infrastructure networks, 2018, https://static1.squarespace.com/static/559921a3e4b02c1d7480f8f4/t/5af00302758d4670efbdc46d/1525678950226/Zischg+Jonatan_782.pdf.
 A. Clauset, C. R. Shalizi, and M. E. Newman, “Powerlaw distributions in empirical data,” SIAM Review, vol. 51, no. 4, pp. 661–703, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Á. Corral and Á. González, Powerlaw distributions in geoscience revisited, https://arxiv.or/pdf/1810.07868.pdf, 2018.
 M. Barthélemy, “Spatial networks,” Physics Reports, vol. 499, no. 1–3, pp. 1–101, 2011. View at: Publisher Site  Google Scholar
 E. Weisstein, “Pareto distribution,” in mathworlda wolfram web resource, 2018, http://mathworld.wolfram.com/ParetoDistribution.html. View at: Publisher Site  Google Scholar
 J. Zischg, C. Klinkhamer, X. Zhan et al., “Evolution of Complex Network Topologies in Urban Water Infrastructure,” in Proceedings of the 17th World Environmental and Water Resources Congress 2017, pp. 648–659, USA, May 2017. View at: Google Scholar
 P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software Environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site  Google Scholar
 D. J. Watts and S. H. Strogatz, “Collective dynamics of “smallworld” networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998. View at: Publisher Site  Google Scholar
 Q. K. Telesford, K. E. Joyce, S. Hayasaka, J. H. Burdette, and P. J. Laurienti, “The Ubiquity of SmallWorld Networks,” Brain Connectivity, vol. 1, no. 5, pp. 367–375, 2011. View at: Publisher Site  Google Scholar
 R. Sitzenfrei, M. Mair, F. TscheiknerGratl, B. Hupfauf, and W. Rauch, “What can we learn from historical water network transition?” in Proceedings of the World Environmental and Water Resources Congress 2015: Floods, Droughts, and Ecosystems, pp. 907–916, USA, May 2015. View at: Google Scholar
 S. Gloeckner, Die Historie der Hochdruckwasserversorgungsanlage in Innsbruck  Hydraulische Modellierung und Ausblick (in German) [Master Thesis], University of Innsbruck, Innsbruck, Austria, 2017.
 M. Schulze, Historische entwicklung der kanalisation innsbrucks und die entwicklung der technischen regelwerke (in German), University of Innsbruck, Innsbruck, Austria, 2012.
 X. Zhan, S. V. Ukkusuri, and P. S. Rao, “Dynamics of functional failures and recovery in complex road networks,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 96, no. 5, 2017. View at: Publisher Site  Google Scholar
 S. Achard, R. Salvador, B. Whitcher, J. Suckling, and E. Bullmore, “A resilient, lowfrequency, smallworld human brain functional network with highly connected association cortical hubs,” The Journal of Neuroscience, vol. 26, no. 1, pp. 63–72, 2006. View at: Publisher Site  Google Scholar
 M. Rubinov and O. Sporns, “Complex network measures of brain connectivity: Uses and interpretations,” NeuroImage, vol. 52, no. 3, pp. 1059–1069, 2010. View at: Publisher Site  Google Scholar
 S. Maslov and K. Sneppen, “Specificity and stability in topology of protein networks,” Science, vol. 296, no. 5569, pp. 910–913, 2002. View at: Publisher Site  Google Scholar
 Z. J. Bao, Y. J. Cao, L. J. Ding, and G. Z. Wang, “Comparison of cascading failures in smallworld and scalefree networks subject to vertex and edge attacks,” Physica A: Statistical Mechanics and its Applications, vol. 388, no. 20, pp. 4491–4498, 2009. View at: Publisher Site  Google Scholar
 M. Mair, W. Rauch, and R. Sitzenfrei, “Spanning TreeBased Algorithm for Generating Water Distribution Network Sets by Using Street Network Data Sets,” in Proceedings of the World Environmental and Water Resources Congress 2014: Water Without Borders, pp. 465–474, USA, June 2014. View at: Google Scholar
 J. Zischg, W. Rauch, and R. Sitzenfrei, “Morphogenesis of Urban Water Distribution Networks: A Spatiotemporal Planning Approach for CostEfficient and Reliable Supply,” Entropy, vol. 20, no. 9, p. 708, 2018. View at: Publisher Site  Google Scholar
 S. V. Buldyrev, R. Parshani, G. Paul, H. E. Stanley, and S. Havlin, “Catastrophic cascade of failures in interdependent networks,” Nature, vol. 464, no. 7291, pp. 1025–1028, 2010. View at: Publisher Site  Google Scholar
 C. D. Brummitt, R. M. D'Souza, and E. A. Leicht, “Suppressing cascades of load in interdependent networks,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 109, no. 12, pp. E680–E689, 2012. View at: Publisher Site  Google Scholar
 A. Vespignani, “Complex networks: The fragility of interdependency,” Nature, vol. 464, no. 7291, pp. 984985, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Jonatan Zischg et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.