The classic modularity index for community detection in complex networks was recently tailored to water distribution networks (WDNs) and extended in order to be cut-position sensitive. Next, the WDN-oriented modularity index was enhanced in order to overcome the resolution limit of the classic modularity. Nonetheless, the modularity-based metrics developed so far allow the networks to be segmented into modules/segments that are similar to each other according to specific pipe characteristics (e.g., pipe lengths, distributed demand, background leakages, etc.). The present work extends and proves the strategy to overcome the resolution limits focusing on an infrastructure index that drives WDN segmentation toward modules that are internally similar with respect to given attributes (e.g., pipe diameters, average pipe pressures, average pipe elevations, etc.), since this aim is suitable for several practical purposes. The introduction of the *attribute*-based infrastructure index permits a comprehensive discussion and comparison of the metrics for infrastructure network segmentation through simple examples. Finally, the practical implications of increasing the resolution of internally similar modules are demonstrated on a well-known benchmark WDN considering various pipe attributes.

## INTRODUCTION

Water distribution networks (WDNs) are essential for all human activities in urban areas. The complexity in analyzing, managing and planning works on such infrastructures stems from their large size (up to thousands of pipes), the underlying hydraulics, as well as the alteration of asset conditions from their original installation. A pragmatic approach to understand WDN real behavior and support effective decisions resorts to segmenting the system into smaller portions (named *districts, segments* or *modules*) suited for different technical purposes including monitoring (e.g., district metering areas), control (e.g., pressure control zones) or even WDN modeling (e.g., calibration). Nonetheless, segmenting a real WDN is not a trivial task due to its size, the looped topology, its hydraulic functioning and because cuts that separate modules are actually costly devices (e.g., valves, flow/pressure gauges) to be installed at technically feasible locations (i.e., in accessible vaults/manholes at pipes end nodes). In addition, a WDN segmentation designed to match a technical purpose might not be adequate for another different scope.

The interest in this topic is documented by many contributions where segmentation is analyzed for various purposes including reliability analysis (e.g., Jacobs & Goulter 1988; Yang *et al.* 1996), location of isolation valves (e.g., Walski 1983), and analysis of contaminant spread (e.g., Davidson *et al.* 2005). Other studies exploited concepts from graph theory to identify the main structure in WDN for monitoring and control goals (e.g., Deuerlein 2008; Perelman & Ostfeld 2011; Alvisi & Franchini 2014) like model calibration, metering water consumption, early contaminant detection, control of pressure/leakages, and network vulnerability analysis (Yazdani & Jeffrey 2012).

A recent approach (Scibetta *et al.* 2013; Diao *et al.* 2013) faced the problem of segmenting WDNs by applying the concepts of community detection (e.g., Fortunato 2010) from complex network theory (Albert & Barabasi 2002; Newman 2010). The initial contributions in this area analyzed the application of the classic modularity concept (Newman & Girvan 2004) to identify WDN modules. However, it was observed by some authors (e.g., Barthélemy 2011; Giustolisi & Ridolfi 2014a) that strong differences exist between the immaterial networks (e.g., food web, trade, World Wide Web), for which the classical modularity concept was conceived, and infrastructure networks (e.g., gas, electricity, water). In particular, WDNs have material links (pipes), spatial constraints (two dimensionality, urban layout, location of water sources/demands, etc.), and material devices installed along links (e.g., valves, pumps, meters, etc.).

Giustolisi & Ridolfi (2014a) introduced a WDN-oriented modularity index that accounted for the WDN infrastructural peculiarities. In short, WDN-oriented modularity was (i) sensitive to position of cuts (since devices in WDN are installed close to nodes instead of middle of links, as originally assumed for immaterial networks); (ii) its formulation was based on the number of pipes into the modules instead of nodal degree; and (iii) pipe features were introduced in terms of *weights* of network links to drive the segmentation process.

Consistently with classic modularity concept, the *weight*-based WDN-oriented modularity index allows identification of modules that are similar to each other. From a technical perspective, it was reported to be of direct relevance for length-based pipe characteristics (e.g., total water demand distributed along pipes, propensity to pipe background leakages, etc.). Indeed, maximizing such WDN-oriented modularity index is likely to return modules that are suited for technical purposes such as, for example, water consumption metering and leakage monitoring.

The problem of designing WDN modules was formulated as a multi-objective problem where the WDN-oriented modularity index should be maximized with the minimum number of cuts (i.e., costly devices to be installed).

Unfortunately, both classic and WDN-oriented modularity indexes are known to suffer from a resolution limit (Fortunato & Barthélemy 2007), which prevents further maximizing the value of the metric by increasing the number of cuts, after a threshold number of modules in the network is reached. Giustolisi & Ridolfi (2014b) analyzed the resolution limit of the *weight*-based modularity and proposed a new infrastructure modularity index to overcome such limit.

In addition, the first work by Giustolisi & Ridolfi (2014a) reported an *attribute*-based variant of the WDN-oriented modularity index that was conceived to maximize the similarity of pipes within each module. From a technical perspective, the *attribute*-based index was conceived to exploit pipe features not strictly related to pipe length (e.g., pipe diameter, average elevation, average pipe pressure, etc.). Thus, modules identifiable by maximizing the *attribute*-based index are better suited for other practical purposes such as, for example, WDN model calibration, pressure control, or leakage control.

Although the referenced works provided relevant innovations on the modularity-based approach for WDN segmentation, the framework of the WDN-oriented modularity indices need to be completed and explicitly framed from a technical perspective.

The present work aims at filling this gap by introducing and discussing the *attribute-oriented infrastructure index* that extends the strategy for mitigating the resolution limit to the *attribute*-based index.

A comprehensive framework of the segmentation metrics (directly based on modularity index or simply recalling the structure of that index) is presented along with a discussion on practical implications from WDN management perspectives.

Simple examples clarify the key concepts and the differences among the infrastructure segmentation metrics and provide thoughtful criteria for practitioners to select the best one according to specific technical purposes. Finally, the TOWN-C (Ostfeld *et al.* 2012) water distribution network is used to discuss the practical implications of increasing the resolution of modules by using the *attribute-oriented infrastructure segmentation index*, considering diameters or average elevation as pipe attributes.

## BRIEF ON WDN-ORIENTED MODULARITY INDEX

*n*is the number of network links/pipes,

_{p}*A*are the elements of the adjacency matrix,

_{ij}*P*is the expected fraction of pipes between vertices/nodes

_{ij}*i*and

*j*in the null/random network (i.e., the expected number of pipes in the network if they were randomly distributed), M

*is the identifier of network modules,*

_{i}*δ*is the Kronecker's delta function to apply the summation to the elements of the same module (i.e.,

*δ*= 1 if M

*M*

_{j}=*and*

_{i}*δ*= 0 otherwise), and summation runs on all the possible node couples (

*i*,

*j*), with

*i*≠

*j*. In Equation (1), the expected fraction

*P*is computed using node degree

_{ij}*k*(

_{i}*k*) of the

_{j}*i-*th (

*j-*th) node, i.e., the number of pipes incident in the node. The metric behind the modularity index measures the strength of a network or graph division into communities/modules. Hereinafter, we will use the word ‘module’ as it is more usual for infrastructure networks.

*n*is the number of pipes linking modules of the infrastructure, namely the number of ‘cuts' in the network (i.e., the decision variables of the WDN segmentation problem) and

_{c}*n*is the number of network modules. The summation inside the square brackets is related to pipe weights stored in the vector

_{m}**w**

*, whose sum is*

_{p}*W*, and Kronecker's delta function

*δ*makes that the sum refer only to the weights of pipes belonging to the

*m*-th module (i.e.,

*δ*= 1 if M

*=M*

_{m}*and*

_{k}*δ*= 0 otherwise). The WDN-oriented modularity is cut-position sensitive and counts pipes belonging to modules instead of the nodal degrees of the classic modularity. In fact, in WDN the pipe cuts occur in the actual position of the devices segmenting the network instead of in the middle of links, differently from what is assumed in classic modularity index when immaterial networks are considered.

*n*(Giustolisi & Ridolfi 2014a) and the following constraint holdsEquation (3) allows us to explain that the WDN-oriented modularity index, as well as the classic modularity index, measures the similarity of the modules to each other. In fact, if we writeit is evident that maximizing the metric

_{m}*Q*means finding the best set of cuts (in terms of number and positions) that generates segments that maximize

*Q*

_{2}. Since

*Q*

_{2}is the negative summation of the squares of

*n*numbers whose sum is unitary (see Equation (3)), it is maximized if the modules are similar to each other as much as possible depending on the topological distribution of the pipe weights in the network (Giustolisi & Ridolfi 2014b).

_{m}The WDN-oriented modularity index in Equation (2) is known to suffer from the resolution limit that occurs because there is a bound of the metrics of Equations (1) and (2) to the identification of small size modules. In fact, the two components, *Q _{1}* and

*Q*are conflicting with respect to the number of modules

_{2}*n*and a mathematical dominance of

_{m}*Q*

_{1}with respect to

*Q*

_{2}(namely the value of

*Q*

_{1}is always larger than

*Q*

_{2}), always occurs. This fact generates a sort of barrier for the identification of small modules whose value depends on the size of the network (Fortunato & Barthélemy 2007; Giustolisi & Ridolfi 2014b).

*infrastructure modularity index*It is obtained starting from the WDN-oriented modularity index and adding the term (

*n*– 1)/

_{m}*n*that represents the minimum theoretical fraction of pipes to be cut to obtain

_{p}*n*modules. This biases the term

_{m}*Q*

_{1}in Equation (4) by a function of the number of modules (

*n*) and eliminates the mathematical dominance of

_{m}*Q*

_{1}vs.

*Q*

_{2}. The

*infrastructure modularity index*has further good features with respect to the WDN-oriented modularity index which makes it more effective for multi-objective segmentation design for WDNs, as discussed by Giustolisi & Ridolfi (2014b).

*attributes*which are similar inside each module. It yieldswhere

**a**

*is the vector of pipe attributes,*

_{p}*ā*(N) is the mean value of the pipe attributes of the network N, i.e., of

**a**

*, and*

_{p}*ā*(M

*) is the mean value of the pipe attributes in M*

_{m}*. Function*

_{m}*δ*limits the summation of the pipe attributes to the elements belonging to the same module.

Indeed, for some WDN management purposes, it is technically advisable to segment the network by searching for cuts (i.e., devices) which generate modules with similar internal pipe characteristics like pipe diameters, pipe average pressures/elevations, etc.

Thus, the *attribute*-based index in Equation (6) has only a similar mathematical expression as the modularity-based index of Equation (2), but different properties.

In summary, Equation (2) is the WDN-oriented modularity (*weight*-based) index *measuring the similarity of modules to each other* and Equation (5) is the infrastructure modularity index having the same feature but aimed at eliminating the resolution limit drawback. Equation (6) is a further WDN-oriented (*attribute*-based) index *measuring similarity within each module with respect to a specified attribute*. In this latter case, we use the word *attribute*, instead of *weight*, to indicate a specific pipe characteristic in order to stress the different aim of Equation (6) with respect to Equations (2) and (5).

It should be noted that the constraint to unit of Equation (3), that is the driver for similarity among modules in Equations (2) and (5) does not hold for the *attribute*-oriented index in Equation (6). However, also the *attribute*-oriented index could be affected by the resolution limit drawback. Therefore, the aim of the next section is to extend the infrastructure index to Equation (6) and demonstrate that, although the resolution limit does not strictly exist for Equation (6), the modification of adding the term (*n _{m}* – 1)/

*n*is also beneficial for that

_{p}*attribute*-oriented index.

This is of technical relevance since, according to the multi-objective strategy for WDN segment design, the increase of number of cuts (i.e., costly valves/devices) is justified by an increased value of the adopted metric.

## ATTRIBUTE-BASED INFRASTRUCTURE SEGMENTATION INDEX

The resolution limit concerns the actual possibility to increase the value of the WDN-oriented modularity considering the increase of segmentation by one module using one cut, i.e., the minimum possible number of cuts, starting from *n _{m}* modules. Actually, this means to assume a sequential search of optimal cuts provided that, for generality of discussion, the segmentation with

*n*modules is a global optimum. Thus, the question is if with one cut it is always possible to obtain

_{m}*Q*(

_{a}*n*+ 1,

_{c}*n*+ 1) >

_{m}*Q*(

_{a}*n*,

_{c}*n*) assuming starting from a global optimal division in

_{m}*n*modules.

_{m}*attribute*-oriented metric, it is always possible that dividing one of

*n*modules of the global optimal segmentation generates modules with the same attributes: i.e.,

_{m}**a**

*=*

_{p}*ā*(M

*), thus null values of summation of distances |*

_{m}**a**

*=*

_{p}*ā*(M

*)|*

_{m}*in Equation (7). Therefore, Equation (7) does not strictly hold for*

_{k}*attribute*-oriented metric. In any case, it can be argued that Equation (7) is a ‘soft barrier’ to the identification of small modules also for the

*attribute*-oriented metric.

*attribute*-oriented

*infrastructure index*adding the term (

*n*− 1)/

_{m}*n*In fact, Equation (7) becomesThus, the modification of the

_{p}*attribute*-oriented index is beneficial for the identification of small modules and also has the same effective features for multi-objective optimal WDN segmentation discussed in Giustolisi & Ridolfi (2014b). (i)

*IQ*is conflicting with the objective functions that increase monotonically with the number of cuts (modules), thus making it appropriate for providing a solution to the cost–benefit optimization problems of WDN segment design. (ii) Maximizing

_{a}*IQ*results in the identification of modules that are actually nested in those achievable maximizing the

_{a}*attribute*-based index

*Q*(Equation (6)). (iii) The possibility of increasing the resolution starting from any number of existing modules makes

_{a}*IQ*not influenced by existing (even non-optimal) segmentations.

_{a}## SOME SIMPLE EXEMPLIFYING NETWORKS

*ɛ*is a small value (e.g., the value of the precision of the 32 bit computing environment 2.2204 × 10

^{−16}) which is useful in order to avoid a null denominator when the values in

**a**

*are equal (i.e.,*

_{p}**a**

*(N)).*

_{p}= āEquation (10(a)) and (10(b)) define metrics allowing the division of the network into modules which are similar to each other according to the internal sum of a vector of pipe *weights* which are user-defined in **w*** _{p}*. The metric of Equation (10(b)) overcomes the resolution limit barrier in identifying small size modules during optimal segmentation.

Differently, Equation (10(c)) and (10(d)) are metrics allowing the division of the network into modules whose internal *attributes* are similar to each other according to their distance from the mean value. The pipe *attributes* are user-defined in **a*** _{p}*. The metric of Equation (10(d)) generally overcomes the resolution limit barrier in identifying small size modules during optimal segmentation search.

It is worth recalling that we here distinguish pipe *attributes* from *weights*. For example, the pipe lengths can be seen as *weights* when we sum them in the case of the metrics Equation (10(a)) and (10(b)), while they become *attributes* when we use them in order to segment the network modules with the same internal characteristics using the statistical distance from the mean value.

We refer to the two simple networks shown in Figure 1. The system in Figure 1(a) is a linear network (i.e., fully branched) composed of 48 pipes (*n _{p}*) and 49 nodes (

*n*) delivering water from a reservoir. The system in Figure 1(b) is derived from that in Figure 1(a) closing the odd pipes around a loop. It is composed of 24 couples of loops connected by one pipe, 192 pipes (

_{n}*n*) and 145 nodes (

_{p}*n*).

_{n}### Case study I: *IQ* vs. *Q* using the linear network

The multi-objective segmentation is performed on the network in Figure 1(a) using the metrics defined in Equation (10(a)) and (10(b)) assuming vector **w*** _{p}* equal to the identity vector. Therefore, the segmentation is based on topology and the solution using Equation (10(c)) and (10(d)) provides the trivial segmentation, i.e., the network is already a module with the same internal, constant, attributes

**a**

*=*

_{p}**w**

*.*

_{p}In fact, the segmentation based on the maximization of *IQ* vs. minimization of the number of cuts *n _{c}* already provides the optimal values of

*Q*because

*IQ*is a metric shifted from

*Q*by means of the term (

*n*− 1)/

_{m}*n*depending on the number of modules

_{p}*n*(Giustolisi & Ridolfi 2014b).

_{m}Figure 2 shows that the maximum resolution of the infrastructure modularity *IQ* corresponds to 48 modules, each composed of one pipe as in Figure 3(a), obtained with the minimum number of cuts, i.e., 47 ( =48 − 1), because the network is fully branched. The value of *IQ* (=0.979) depends on the number of pipes (i.e., modules). It is interesting to note that the segmentation considering *Q* allows the division of the network into seven modules with six cuts (see the maximum value of *Q* in Figures 2 and 3(b)), because the resolution limit generates a mathematical barrier for the identification of smaller modules.

In summary, the exercise shows that the *IQ* does not have the resolution limit to allow the identification of any module generated by one further cut in the network, while *Q* has a strong resolution limit increasing with the number of pipes.

It is worth noting that the segmentation solution with seven modules and six cuts, corresponding to the maximum value of metric *Q* (Figure 3(b)), is the same as the solution achievable with six cuts using the metric *IQ*.

### Case study II: *IQ* vs. *Q* using the looped network

This test is similar to the previous one but applied to the network of Figure 1(b). Figure 4 shows that the maximum of the curve *Q* corresponds to 12 cuts while *IQ* corresponds to 24 cuts.

Figure 5(a) and 5(b) show the maximum division into 25 modules with 24 cuts considering *IQ* and into 13 modules with 12 cuts considering *Q*. *IQ* allows each couple of loops to be divided while *Q* does not. This fact is a different confirmation of the effectiveness of *IQ* in identifying modules overcoming the resolution limit also in the presence of loops.

### Case study III: *IQ*_{a} vs. *Q*_{a} using both the networks

_{a}

_{a}

Here, we perform two tests about the linear and the looped networks of Figure 1, but the metric *IQ _{a}* of Equation (10(d)) is applied instead of that in Equation (10(b)). To this purpose the selected attribute in

**a**

*are the pipe diameters. Dummy diameters in the range [1, 12] are assumed. They are constant for each of four contiguous pipes starting from one ending node of the linear network. The looped network has the diameters of the loops constant and equal to that of the basis linear network.*

_{p}Figure 6(a) and 6(b) show that the maximum number of modules of *IQ _{a}* is much higher than

*Q*demonstrating the effectiveness of adding the term (

_{a}*n*– 1)/

_{m}*n*also to

_{p}*Q*. In addition, the maximum number of cuts of

_{a}*IQ*is 11 which corresponds to 12 modules in both cases, as better reported in Figures 7(a) and 8(a). Therefore,

_{a}*IQ*is effective to identify each group of identical pipes both in linear and in looped networks while

_{a}*IQ*would identify the same 24 modules described in cases I and II since the topological condition prevails.

Figures 7(b) and 8(b) shows that *Q _{a}* is unable to identify the 12 groups of pipes because of the resolution limit occurring also for the attribute-based modularity of Equation (10(c)). It follows the effectiveness of the modularity of Equation (10(d)) to separate modules with the same diameters, overcoming the modularity of Equation (10(c)).

In conclusion, the three case studies demonstrate that both the metrics – namely, the pipe *weight*- and pipe *attribute*-based indices – are enhanced by adding the term (*n _{m}* – 1)/

*n*. Furthermore, case study III demonstrates the effectiveness of the

_{p}*attribute*-based metric to separate modules based on the assumed pipe characteristics.

## TOWN-C CASE STUDY

This section compares the results of segment design achievable by using the *attribute*-based index (*Q _{a}*) and the

*attribute*-based infrastructure index (

*IQ*) metrics on the TOWN-C water distribution network. This network is chosen because it is well-known in the technical literature (data are available as supplementary material of Ostfeld

_{a}*et al.*(2012)) and its layout allows segments and relevant cuts to be clearly visualized. In addition, the presence of hydraulic devices (i.e., pumping stations and tanks) already installed in the network results in a scenario that is closer to the real context of designing segments into an existing system.

The analysis considers two attributes focusing on different technical purposes of WDN segmentation: the pipe diameter and the average elevation.

As for previous examples, the segmentation solutions descend from the two-objective optimization where the chosen WDN-oriented modularity index (*IQ _{a}* or

*Q*) should be maximized with the minimum number of cuts (devices) (see Giustolisi & Ridolfi (2014a) for optimization problem formulation).

_{a}### TOWN-C case (a): diameter-based metrics

Identifying WDN modules composed of homogeneous pipe diameters has technical relevance from both asset management and hydraulic modeling purposes. In fact, pipe diameters generally reflect the hydraulic functioning ranging from larger trunks, mainly to transport water (e.g., from the water sources), to smaller pipes, mainly used to distribute water to users (e.g., in the peripheral WDN areas). Thus, looking for modules composed of similar pipe diameters is expected to return WDN portions with different preeminent hydraulic functioning. This information can be useful for many practical purposes including, for example, the selection of candidate location of gate valves to isolate WDN portions where water is distributed (e.g., in case of malfunctioning) without affecting the main transport lines.

In addition, internal hydraulic resistances of pipes strongly depends on pipe diameters (i.e., on power 5 based on Darcy-Weisbach head loss formulation), thus it is of direct relevance for WDN model calibration purposes. Grouping pipes that are expected to show similar values of hydraulic resistance per unit length is a pragmatic way to reduce the number of unknowns of the calibration problem by introducing technical/engineering insight. Identifying WDN modules with homogeneous pipe diameters allows preservation of information on the network topology, returning modules with similar and contiguous pipes; thus, going beyond strategies for clustering the pipe database only. From such perspectives, the cuts that separate modules represent the most effective points to allocate pressure/flow sampling devices in order to maximize the observability of pipe hydraulic resistances (Giustolisi & Berardi 2011).

Figure 9(a) and 9(b) report WDN segments obtained by minimizing the number of cuts (e.g., pressure/flow meters) and maximizing the index *Q _{a}* (a) and

*IQ*(b) using pipe diameters as attribute. The metric

_{a}*Q*is maximized for 18 modules and 18 cuts. It is evident (see Figure 10(a)) that links representing pumps (i.e., with the same diameter) are in the same modules, as expected. In addition, the largest trunks, mainly used to transport water from tanks/reservoirs, are identified as belonging to the same modules; while some branched WDN portions, mainly used for water distribution to customers, are identified as separate segments. Nonetheless, in some cases the identified modules are visibly composed of portions with different diameters that should be separated from the rest of the network; this is the case of some branches connected with the upstream looped network portions.

_{a}Figure 9(b) shows that maximizing the *attribute*-oriented infrastructure index *IQ _{a}* allows this resolution limit to be overcome, returning up to 83 modules. It is worth noting that, consistently with the expected WDN hydraulic behavior, three types of sub-modules are identified: (i) branched portions separated by one cut only; (ii) looped inner portions, separated by multiple cuts from the rest of the WDN; (iii) transport modules, linking the WDN to tanks/reservoir or to other distribution modules.

From a WDN model calibration perspective, the layout of modules achievable by maximizing *IQ _{a}* is consistent with the general criteria for collecting pressure/measurement at nodes where such measures are maximally informative for the calibration of variables. This is the case of many cuts separating branches, since the measurements collected at these points actually maximize the topological observability of the unit hydraulic resistance for homogeneous pipes in that module (Walski 1983; Giustolisi & Berardi 2011).

### TOWN-C case (b): elevation-based metrics

Elevation is an essential driver to model and run water distribution systems fed by gravity or by pumps, since differences in elevations affect pressure regimes and, in turn, the capacity of satisfying water demand as well as the leakage outflow from pipes. Accordingly, considering average pipe elevation as the attribute for identifying modules has two technical purposes. On the one hand, this is a way to design WDN modules that are expected to experience similar pressure regime. On the other hand, due to the multi-objective strategy used to design modules, resulting solutions entail the minimum number of cuts separating WDN zones with homogeneous average elevation (and expected pressure regime) inside the modules but different elevation between contiguous modules. Thus, cuts are suited to representing the optimal location of pressure control devices (e.g., PCV) that are usually required to be the fewest possible, due to their cost of installation and maintenance. Although minimizing such devices (i.e., cuts) directly reflects economic criteria of water utilities, increasing the resolution of elevation-based modules is expected to provide additional information on the most suitable candidate nodes to be used as pressure set points of PCV (i.e., controlled by remote pressure readings).

Figure 10(a) reports the 11 modules entailing the maximum *attribute*-based index *Q _{a}* (i.e., Equation (10(c)) achievable with the minimum number of cuts (equal to 10). It is evident that modules correspond to differences in elevation and in some cases (shadowed in Figure 10(a)) cuts indicate feasible candidate location of pressure reduction valves to control branched WDN portions. Nonetheless, the resolution limits of

*Q*result in some inconsistencies such as, for example, pipes joining tanks (i.e., at higher elevation) still belong to modules that have lower average elevation.

_{a}When the *attribute-based infrastructure index IQ _{a}*, i.e., Equation (10(d)) is used, the increased resolution identifies 88 modules. It is evident that a number of branches are identified as different modules. However, this happens only if branches serve areas with different elevation from the nearest upstream WDN module. In fact, in some cases, like the shadowed module in Figure 10(b), the differences in elevation are not statistically significant to justify the creation of a new module; thus this module is the same as that of Figure 10(a) (i.e., based on

*Q*). A similar behavior occurs in the area near to tank T1 that is located on a flat area, where the increased resolution allows identification of branches as separate modules while the main looped WDN portion still belongs to the same module.

_{a}It should be mentioned that tanks and reservoir are all located in modules apart from the rest of the network, consistently with their different elevation and WDN hydraulic behavior. In addition, the long pipeline connecting the main pumping station (near to the reservoir) to the rest of the network is divided into three modules since its elevation drops from 56 m to about 12 m above sea level.

## CONCLUSIONS

The modularity concept has been recently borrowed from complex network theory to infrastructure networks and has been tailored for WDN analysis and management. The work by Giustolisi & Ridolfi (2014a) introduced the base formulation for WDN-oriented modularity indexes that were aimed at matching the peculiarities of WDN infrastructure as well as the technical meaning of pipes and modules.

The *weight*-based modularity index was introduced to maximize the similarity of modules with each other. Starting from the *weight*-based modularity index, Giustolisi & Ridolfi (2014b) proposed an *infrastructure modularity index* in order to overcome the resolution limit that resulted from the original formulation of the classic modularity index upon which the WDN-oriented modularity was developed.

In addition, the same authors introduced an *attribute*-based index that was suited to identify modules with the maximum similarity of the attributes within each module. This means that returned modules do not simply entail groups of similar pipes through the network, but implicitly preserve the information on WDN topology (i.e., modules with similar and contiguous pipes), which is essential for an infrastructure management perspective.

The present contribution demonstrates that also the *attribute*-based index can suffer from resolution limit and extend the concepts of *infrastructure index* variant also to *attribute*-based index.

The infrastructure segmentation metrics are better suited to identify segments in a multi-objective optimization paradigm, where the *weight*-based and *attribute*-based indices should be maximized while minimizing the number of cuts (i.e., costly devices). Also, it was found that the modules identified using the *Q _{a}* index are identified when the

*IQ*is used; this means that the additional modules identified by maximizing

_{a}*IQ*are actually nested in the previous ones.

_{a}The proposed didactical examples demonstrate the advantages in identifying modules by using the *infrastructure* metrics *IQ* and *IQ _{a}* instead of the

*weight*- and

*attribute*-based indexes

*Q*and

*Q*, respectively. Finally, the well-known TOWN-C literature network is used to discuss the practical implications of the increased resolution of segmentation by using the

_{a}*attribute*-based infrastructure index

*IQ*. Pipe diameters and average elevations were selected as pipe attribute, respectively. In all cases, the resolution limits that are typical of

_{a}*Q*clearly prevent identifying modules that are more suited for possible final technical purposes of the segmentations (e.g., ranging from WDN model calibration to pressure control and leakage detection plans). Differently,

_{a}*IQ*overcomes such limits and gives a very detailed and technically sound network segmentation.

_{a}## ACKNOWLEDGEMENT

The research reported in this paper was founded by the Italian Scientific Research Program of National Interest ‘Tools and procedures for an advanced and sustainable management of water distribution systems’ – PRIN2012 (Prot. 20127PKJ4X) – Italian Ministry of University and Research (MIUR).