A System Dynamics Model and Analytic Network Process: An Integrated Approach to Investigate Urban Resilience

: During the last decade, the concept of urban resilience has been increasingly implemented in urban planning, with the main aim to design urban development strategies. Urban resilience is a multi-dimensional and dynamic concept. When applied to urban planning, it consists of studying cities as complex socio-economic systems. Municipalities are currently working to undertake appropriate actions to enrich the resilience of cities. Moreover, several difficulties concern the evaluation of the impacts over time of the strategies designed to enhance urban resilience. The present paper proposes an integrated approach based on the System Dynamics Model (SDM) and the Analytic Network Process (ANP). The objective of this research is to describe the method and to illustrate its application to the area called Basse di Stura, located in the city of Turin, Italy. The method is applied to evaluate the possible impacts of two different urban scenarios in terms of the change of urban resilience performance over time. The final result is represented by an index that describes urban resilience performance.


Introduction
Cities across the world are extremely exposed to different hazards and risks, both natural and man-made [1][2][3][4][5]. These risks include all of the dimensions of cities, from the environmental to the socio-economic, which are mainly related to social vulnerability [6]. Thus, the concept of resilience is recognized as a critical paradigm in the actual urban agenda to make cities able to face and react to these uncertainties [1, [7][8][9][10]. When applied to the urban planning field, urban resilience does not correspond to the ability of the urban system to return to its previous state as defined by the general meaning of resilience [9]. Instead, it has to be interpreted as the capacity of the urban system to absorb, to adapt and to respond to changes [6,7,11,12]. However, only a few cities have undertaken actions and strategies to improve their resilience [8,[13][14][15][16][17]. This condition is determined by the difficulties related to the definition of urban resilience and its measurement [18,19]. This concept has been analyzed from the point of view of several disciplines, giving rise to several definitions in the literature [20,21], such as in climate change [22] or risk reduction [23] studies. However, a literature review clearly shows that urban resilience is a multi-dimensional and dynamic phenomenon that characterizes cities as complex systems [7,12,20,[24][25][26]. Starting from the assumption "if you can't measure it, you can't improve it" [11,27], urban resilience measurement has become a fundamental aspect concerning the definition of strategies to enhance urban resilience [28][29][30][31][32]. Several methods have been proposed in the literature, that are both qualitative and quantitative [8,18,19,33,34]. Moreover, existing methods are considered to be underdeveloped in terms of evaluating urban resilience within its multi-dimensionality, complexity and dynamic behavior over time [14,18,19]. The main problem recognized in the existing methods is the lack of mutual interdependences among the different dimensions of urban systems [14,18,19,33,35,36]. The challenge is to provide an assessment method that is able to foresee the resilience of urban system in the future, starting from the current state and according to development strategies [37][38][39].
In this paper, an integrated evaluation approach based on Analytical Network Processes (ANP) [40,41] and System Dynamics Model (SDM) [42][43][44] is illustrated, according to their methodological background and operative characteristics. The main objective of this paper is to explore the ability of this approach to analyze the relationships among factors that afflict urban resilience and its role in the evaluation of actions to enhance urban resilience [45,46].
For this purpose, this paper illustrates the application of the proposed approach to analyze and evaluate two different scenarios (Inertial and Strategic) developed for a downgraded area in the city of Turin (Northern Italy) [47][48][49], in terms of their possible impacts on urban resilience performance over time [50].
The paper is structured as follows: Section 2 illustrates the integrated approach; section 3 is focused on the application of the proposed integrated approach to the real case study; section 4 is related to the discussion of the results; and section 5 concerns the conclusions and the final remarks related to future research perspectives.

Method
The proposed method is an integrated evaluation approach based on the combination of the System Dynamics Model (SDM) and the Analytic Network Process (ANP). In this section, these two techniques are illustrated, highlighting their methodological background and operative characteristics.

Properties and State of the Art
The Analytic Network Process (ANP) is a special method of the wide set of multi-criteria analysis methods (MCA). ANP is often used to analyze complex issues related to the decision-making process, due to its power to deal with qualitative and quantitative attributes and to identify the relationships between decision elements [51][52][53]. It is recognized as a generalization of the Analytic Hierarchy Process (AHP) [54][55][56]. Generally, multi-criteria techniques are based on the principle of independency of the elements involved in the evaluation. In contrast, the peculiarity of the ANP is its capability to identify and to consider the mutual dependencies between the criteria of the evaluation process [57][58][59]. In detail, the interrelations that occur between the elements are defined as "feedbacks". The basic structure of the ANP is represented by a network of clusters and nodes, i.e., the elements, that are contained in the clusters (Figure 1). The clusters represent both the criteria and the alternatives [53].   The network structure of the ANP allows us to represent the decision problem within its complexity [54][55][56]. As a matter of fact, interdependence among different elements of the problem can occur, such as: alternatives, criteria, sub-criteria and goals. In detail, interconnections are represented by lines, connecting elements or nodes to each other. Figure 1 illustrates the general structure of an ANP network where: (1) feedback connect components to other elements and (2) loops connect a component to itself [40]. Furthermore, in the network it is possible to identify the source node, that is the origin of the path of influence, and to identify the sink node, that is the final destination of influence [54][55][56]60]. ANP is recognized as a suitable decision-support tool, able to face those decision problems characterized by high level of interactions and interdependences between the elements. It is employed to analyze those issues which cannot be structured hierarchically [51][52][53]57]. However, the literature on the implementation of an ANP in the assessment of urban transformation strategies is not vast [61]. Specifically, the literature about the integration of ANP with SDM is very limited [62]. Furthermore, the current literature about the application of ANP in resilience evaluation is mainly related to the risk and vulnerability assessment [63], focusing on risk management and risk reduction [64,65]. In this context, this paper illustrates one of the first applications of the integration of ANP with SDM to the urban resilience assessment field.

Methodology
The main steps of the ANP are summarized as follows:

Decision Problem Arrangement
The decision problem has to be included in the network structure ( Figure 1). The first operation of this phase is to identify the evaluation goal. Once the goal has been defined, the second step is to recognize criteria (clusters/nodes) and alternatives. Furthermore, the relationships which occur among the elements must be identified during the development of the network. Different methods can be applied to obtain the network structure, such as focus groups and brainstorming [53].

Pairwise Comparison
In ANP, the pairwise comparison is applied to compare the elements and to establish their preferences [55,56,59]. Indeed, the different priorities are identified using the pairwise comparison and judgment both for clusters and nodes. Specifically, clusters are weighted by considering both their single role and impact and their mutual dependence among the clusters. The relative elements' relevance is evaluated through the Saaty nine-point scale (Table 1). In this phase, the internal importance vector , that stands for the relative importance of elements or clusters, is calculated through (Equation (1)): where A is the pairwise comparison matrix of the criteria, w is an eigenvector and max is the largest eigenvalue. Moderately / more important One decision element is moderately more influential than the other.

5
Much more important One decision element has more influence than the other 7 Very much more important One decision element has significantly more influence over the other.

9
Extremely more important The difference between influences of the two decision elements is extremely significant.

2,4,6,8
The difference between influences of the two decision elements is extremely significant.
Judgment values between equally, moderately, much, very much and extremely influential.
3. Supermatrix formulation The preferences obtained from the pairwise comparison are illustrated through the supermatrix ( Figure 2). This matrix is defined as a block matrix in which the elements of each block are defined by the relationships between and within the levels, as represented in the network [53].  Figure 2 illustrates the general structure of the supermatrix, in which Ci represents the ith cluster; for i=1,…N, eij indicates the jth element within Ci, for j=1,… , Wij is the block matrix that involves the priority weight vectors (w) of the influence of the elements in the ith cluster with respect to the jth cluster. If the ith cluster has no influence on the jth cluster itself (a case of inner dependence), then the elements of Wij are zero.
The supermatrix can support both the identification and the understanding of the relationships among the network. Furthermore, it provides the information about the priority of the elements that characterize the decision problem [54][55][56]59].

Selection of the preferable alternative
The final step of the ANP development are the calculation of the weighted supermatrix and final supermatrix. This matrix is processed to get a stable set of weights, that characterize the different elements [58]. The supermatrix is structured to get the global priority vector (Equation (2)) lim → (2) where W k is the supermatrix. System Dynamics Models (SDMs) are both computer-aided and theoretical tools that are employed to analyze dynamic and complex systems, characterized by mutual feedbacks and interdependences [66]. SDMs are based on the approach of the System Dynamics (SD). SD was introduced by Forrester [42][43][44] at the end of the 1950s to improve the understanding of the strategic processes related to complex systems [67]. The theory of SD is a suitable tool to describe complex and dynamic systems through the identification of the mutual relationships that occur among variables [45,68]. Furthermore, this approach is employed to analyze the impacts of different factors on defined objectives, providing useful information for the decision-makers [69]. The theory of System Dynamics gives the fundamental guidelines to structure the model that describes the system through the employment of feedback loops, stocks and flows. The System Dynamics Model is grounded on feedback concepts to manage non-linearity, multi-loop and time-lag characteristics of complex dynamic systems. SDM can be considered as an efficient scenario simulation approach that is able to reveal the likely temporal behavior of complex systems [68]. Therefore, it can be also used to support the definition of policies for sustainable development [70,71]. As a matter of fact, an increasing amount of SDM applications have been recognized in the literature in the last few decades [72]. In detail, SDMs are nowadays applied for dynamic and complex systems that involve a high level of uncertainty. Important fields of application are: transportation [73], land use [74], environment management [75], waste management [76], project management [77] and sustainable urban development [46,72,78,79]. However, there are few existing studies that have applied SDMs to simulate and to evaluate urban resilience within its evolution over time [38,80]. Furthermore, despite the ability of SDM to deal with scenario assessment [78], its integration with MCA in urban development and urban resilience assessment fields is quite limited. For this reason, this paper aims to highlight the potentiality of SDM integrated with MCA [81] in the context of urban resilience measurement.

Methodology
The main steps to structure a SDM can be summarized as follow [67]:  Identify the problem. This is the first step to implement the SDM. This phase concerns the identification of the problem, in order to set the boundaries of the system and define the specific objective;  Develop a dynamic hypothesis. This step is referred to a theory that allows us to consider the mode behavior over time. In detail, the dynamic hypothesis is a conceptual model that consists of a causal loop diagram, stock and flow diagram or their combination;  Create the basic structure of a causal graph. This phase concerns the processing of the causal loop diagram with the identification of the most important elements of the system and the relationships that occur among them;  Augment the causal graph with more information. In this step, secondary and tertiary elements of the system are identified, in order to describe it in accurately;  Convert the augmented causal graph to a system dynamics flow graph. This step refers to the translation of the causal loop diagram into the stock and flow diagram that is used to simulate the dynamic behavior of the system over time;  Translate a system dynamics flow graph into programs or equations. This phase is fundamental to get the representation of the dynamic behavior of the system over time. In detail, this section describes the phases of creating a basic structure of causal graph and a system dynamics flow graph. The basic structure of a causal graph is represented by the causal loop diagram ( Figure 3) [78]. This diagram is used both to recognize the mutual influences which occur among the variables of the systems and to develop the feedback loops [45,69,74]. The processing of the causal loop diagram is the first phase to describe the dynamic behavior over time of the complex system [38,82]. Moreover, the mutual influences among variables can be either positive or negative [67]. A positive relationship means that the two variables will change in the same way. As an example, if variable A increases, then variable B will also increase. Within a negative relationship, the variables will change inversely. For instance, if variable C increases, variable B will decrease. Figure  3 illustrates a general structure of the causal loop diagram, identifying its main elements. The evolution dynamics of the system is represented by stocks and flows diagrams in Figure 4 [9, 10,26]. Specifically, stocks stand for those variables that both describe and characterize the state of the system at any particular time [45,69,83,84]. The second elements of the SDM are flows. These elements illustrate how stocks can change over time [67]. Furthermore, flows can impact stocks via the inflow or outflow and also interconnect the stocks within a system [68]. The final element of SDM model are the conveyors that represent the external inputs to the model. In fact, they are also used to represent different actions for the policy simulation and analysis [46,67]. Figure 4 illustrates the basic elements of a stock and flow diagram.
Furthermore, the stock and flow diagram represents integral finite differential equations that consider the feedback loops. The dynamic behavior of the system is obtained through their integration [67]. In this sense, stock and flow diagram can support the examination of the possible effects over time of different strategies through the simulation model that is grounded on differential equations [38,72,78].

Basse Di Stura Area
The approach illustrated above was applied to analyze the area of Basse di Stura, located in Turin (Northern Italy) ( Figure 5). This area extends over 150 hectares and it is delimited by three main arterial roads and by the Stura river ( Figure 6). The neighboring areas of Basse di Stura count 126,666 inhabitants. It is considered a critical part of the city of Turin due to its social and environmental conditions [85], meaning any serious threat would be a major concern for a large part of the city's community [86,87]. 3. Roma camps. These represent the main social challenge for the Basse di Stura area ( Figure 6).
Several Roma camps are located in the Northern part of the Basse di Stura area. Most of these settlements are not regular and they create several health and hygienic issues [86]. The areas surrounding Basse di Stura are residential areas with a high population density [86]. In this context, the issues and the pressures that characterize this area can afflict a great portion of local community ( Figure 6).
For these reasons, the public administration is actually interested in intervening to improve the general condition of this area [85]. Thus, this paper proposes the application of the integrated approach to simulate the possible effects and impacts of an Inertial and a Strategic scenario, in terms of changing of the urban resilience performance over time [50]. The area of Basse di Stura has been examined through a multi-dimensional and multi-perspective approach [49]. In detail, this choice has been made considering the peculiarity of the examined area and the effects that its conditions could have on different aspects of the city of Turin, especially in environmental and social dimensions (e.g., the health and safety of inhabitants). The final objective of this application is identifying which strategy could improve on the initial conditions over time, enhancing urban resilience.

Strategic Scenario
The strategic scenario describes the possible actions that could be undertaken to improve the current conditions of the Basse di Stura area [88]. These interventions have been defined through an in-depth analysis of the conditions of the area. The strategic scenario is based on the following actions: 1. Retrieve the brownfield areas. The target is reclaiming some of the brownfield areas ( Figure 7); 2. Restoring some of the abandoned existing buildings, to insert some new activities to improve both social and economic conditions ( Figure 7); 3. Reducing the soil consumption. The main aim is trying to maintain minimum extensions of the permeable surface. This objective is mainly due to the proximity of this area to the Stura river, while also considering its environmental potentiality; 4. Enhancing better conditions for the socio-economic dimension, through the hypothesis of creating new businesses in the restored buildings. Figure 7 illustrates these actions through a schematic representation. The parameters through which the strategic scenario has been characterized in the SDM model for the scenario simulation will be explained in the following section (section 4).

Application
As mentioned in the previous section, the integrated approach has been applied to analyze and evaluate the possible effects of two different scenarios on the area of Basse di Stura, according to the changing of the urban resilience performance. The specific application is detailed and described in the following sections.

Structuring of the model
The proposed integrated approach has been performed through a multi-steps structure. Figure  8a shows the general framework and Figure 8b illustrates how the general framework has been applied and adapted for the specific case study.

Identification of Indicators
After the goal of the evaluation is defined (Figure 8), the subsequent step foresees the identification of the indicators related to urban resilience concept to structure the SDM.
In line with the literature of urban resilience indicators [6,89,90], in this paper we consider the indicators that refer to the following five dimensions: (1) economy, (2) society, (3) built environment, (4) environment (5) governance. To better identify specific and accurate indicators for the area of Basse di Stura, experts followed the questions proposed by Meerow [20] as conceptual guidelines. Table 2 lists the set of 15 multidimensional indicators specifically fitted for this case study of Basse di Stura.

Causal Loop Diagram
Once the indicators have been identified, the subsequent step consists of the implementation of the SDM. The first phase is to build up the causal loop diagram (Figure 9). The causal loop diagram ( Figure 9) has been processed during an expert focus group, in which the basis of the cognitive mapping was recalled [82,91]. In fact, experts were asked to identify the mutual relationships among the considered indicators (Table 2), while also considering their mutual influence [46]. In detail, this diagram has been built both on the literature review on urban dynamics [78,92,93] and on the experts' knowledge. Figure 9 illustrates the causal loop referred to the Basse di Stura area, in which the relationships among social, economic, environmental and governance indicators have been identified. As an example, the variable "Soil Consumption" has a negative influence on "Permeable Surface".

Structuring the ANP model
In this application, ANP has been selected for its ability to consider the dependences between the evaluation criteria. Moreover, the considered indicators (Table 2) have been organized in cluster and nodes, as suggested by the literature. Figure 10 illustrates the ANP network developed for the case study of Basse di Stura. In detail, the network is not only a final result of experts' brainstorming, but it is also built up on a causal loop diagram (Figure 9). In fact, all the relationships recognized in the causal loop have been also inserted in the ANP network. Thus, it is possible to observe the influence between indicators and clusters, indicators and goals and also feedback on the same elements. As an example, the ANP network shows the feedback of the "Environment" cluster, that is determined by the influences of its nodes as determined in the causal loop diagram.

Determination of the weights of indicators
As mentioned in Section 2, ANP requires the pairwise comparison to establish the priorities between indicators involved in the evaluation. Experts were asked to assess the relevance of the considered indicators in accordance with the Saaty's fundamental scale (Table 1). The first step of the evaluation concerns the pairwise comparisons among the clusters. Thus, the matrixes with the assessment performed through the pairwise comparison have been developed. These matrixes show the numerical values that represent the influence identified for the elements of the network. Once we established the priorities of the clusters, the subsequent step was performing the pairwise comparison for the nodes. As for the examination of the clusters, the judgments have been made in accordance with the influences and interdependencies recognized in the network. Both the comparison of the clusters and nodes have been obtained through discussion with a multidisciplinary panel of urban planning and assessment experts. As an example, Table 3 and Table 4 illustrate the pairwise comparison for clusters and nodes. Table 3 shows the pairwise comparison between "Built Environment" and "Environment". In detail, the assigned value "3" refers to the fundamental Saaty's scale (Table 1). This means that experts considered the dimension of the environment as being much more important than the built environment.   Table 4 shows the pairwise comparison between the nodes belonging to the environmental cluster. As for the previous pairwise comparison, the value assigned are used with Saaty's fundamental scale to determine the relevance of the different indicators. In detail, the assigned values mean that (1) the soil consumption is more important than the permeable surface (value 3); (2) the brownfield is very much important than the permeable surface (value 7); and (3) the brownfield is much more important than soil consumption (value 5). A great importance has been given to the brownfield. This is because the Basse di Stura area is a brownfield area. In this context, it is possible to underline that the pairwise comparison has been made considering the specific condition of the examined decision problem.  Once all the pairwise comparison matrixes for both the clusters and nodes are obtained, the totality of these related priorities develops the unweighted supermatrix (Appendix B). The supermatrix represents the relationships among the nodes that compose the network. Moreover, the unweighted supermatrix has been multiplied for the final vector priorities (Table 5) to get the weighted supermatrix (Appendix C).  Table 5 illustrates the final vector priorities obtained through the evaluations. These priorities give a great importance to (1) the brownfield, (2) health and safety of inhabitants, (3) permeable surface and (4) regeneration, while little importance has been given to the economic aspect. These priorities can be considered as the result of the combination of two different reasoning, that are (1) causal-loop diagram development and (2) pairwise comparison. It is important to underline that through the processing of the causal loop diagram, it was possible both to understand the mutual interdependences between the considered criteria and their possible effects on the Basse di Stura Area. As an example, the second most important criteria is the "health and safety of inhabitants", which is also strictly related to the extension of the brownfield area.

Current Urban Resilience Index
As underlined in the previous sections, the objective of this application is evaluating the inertial and the strategic scenario in order to assess their evolution over time, in terms of urban resilience performance. To pursue this objective, it is necessary to analyze the current urban resilience performance of the area. For this reason, the current urban resilience index has been calculated. This index has been obtained through the Equation (3): where:  Si stands for the final score;  eij represents the normalized value of each indicator;  wi represents the weight of each criterion The equation (3) allows to combine the current performance of the area (Appendix A) for each indicator with the corresponding weights, obtained through the ANP (Table 5).
To obtain the current urban resilience index through the Equation (3), the performance of each indicator has been firstly normalized through the Equations (4) and (5): Normalization allows us to make all the indicators comparable in an interval between 0 and 1. In detail, Equation (4) is applied to normalize those indicators that have positive impacts on the considered objective. In addition, Equation (5) is used to normalize those indicators with negative impacts on the evaluation objective. The index obtained for the current urban resilience performance for Basse di Stura is 0.07. This low performance is based on the bad performance of the Basse di Stura area in (1) Brownfield, (2) Health and Safety of Inhabitants and (3) Buildings exposed to hazards, which are all considered as important in this context (Table 5).

Stock and Flow Diagram
The causal loop diagram has been applied as the basis to develop the stock and flow diagram. As shown in Figure 11, some indicators have been characterized as stocks and others as flows [78]. Strategy inputs have been converted into the conveyors, because these elements are used to describe the external conditions the actions/policies in the model [67]. The variables described as stocks are: (1) Brownfield, (2) Soil consumption, (3) Health and Safety of the Inhabitants, (4) Buildings exposed to hazards and (5) Number of Jobs. In detail, the stocks are those variables that characterize the evolution of the system. In this context, these variables have been translated into stocks in order to assess how they can characterize the Basse di Stura area over time.
To obtain the SDM simulation (Figures 12-14), different balance equations have been developed. These formulas can be generally described as follow Equation (6): As an example, the equation of the variable of soil consumption can be expressed as: Soil Consumption(t) = soil consumption (t -dt) + (new construction)*dt.
In this sense, the changing over time of the soil consumption simulation is the result of the combination of the current soil consumption at the year 0 ( Figure 12) plus the new construction. Furthermore, other specific equations have been developed to implement the SDM to simulate both the possible impacts and influence in enhancing urban resilience over time of the two scenarios considered (Strategic and Inertial) [50].

Results
SDM scenario analysis is useful to provide a possible range of futures, based on policies and actions [78]. In SDM, the values of variables and parameters can be adjusted to describe different actions and scenarios. The results of the simulations will be different based on the different inputs [78]. As previously mentioned, this paper considers two different scenarios: the inertial scenario and the strategic scenario [50]. The inertial scenario stands for no actions policy, while the strategic scenario represents the possible action to undertake for the improvement of the condition of the Bassa di Stura area. To characterize these two scenarios in the SDM model [45,78], different values have been established for the parameters: (1) recovery and restoration; (2) risk reduction; and (3) multistakeholder decision planning. In detail, low values have been established for the inertial scenario, to represent the no intervention policy. Also, medium values have been established for these parameters to represents the actions (Section 3.2) of the strategic scenario in the SDM model. Scenarios simulations have been performed for 15 years through SDM (STELLA, ISEE Systems Inc.), starting from the stock and flow diagram ( Figure 11). Specifically, these simulations have been developed for the variables identified as stocks. The main objective is testing their probable effects over time and how they could affect the urban resilience performance of the Basse di Stura area over time [50]. Figure 12, Figure 13 and Figure 14 illustrate the simulations, through which it is possible to highlight the differences in the evolution over time of these considered scenarios.

Discussion
We have shown how the inertial scenario and the strategic scenario have significant differences in the performances of "soil consumption" and "brownfields" (Figure 12), "health and safety of inhabitants" (Figure 13) and "number of buildings exposed to hazards" (Figure 14).
In detail, the strategic scenario shows a more limited increase in soil consumption than the inertial scenario. This behavior is due to the different values assigned to the "regeneration and requalification rate" parameter. In fact, this parameter has a decreasing impact on "new construction" variable, as highlighted on the causal loop diagram ( Figure 9). Moreover, the parameter of regeneration also has a decreasing impact on the "brownfields" and the "number of buildings exposed to hazards". In fact, the strategic scenario also shows decreasing tendencies over time in these variables. Moreover, the improvement of these conditions also implies a positive evolution in the social variable of the "health and safety of the inhabitants" (Figure 13), as identified in the causal loop diagram (Figure 9). However, no considerable differences have been recognized for the economic variable "number of jobs". These specific results can be explained through two different aspects. Firstly, the actions of the strategic scenarios are mainly focused on intervening on the environment and built environment dimensions. Secondly, the general tendencies of the economic sector of the city have also been to determine the possible tendencies of the number of jobs at this scale of analysis. This analysis has been useful to highlight that the improvement of (1) Brownfield, (2) Buildings exposed to risks and (3) Soil consumption can decrease the pressures that the conditions of the area of Basse di Stura can afflict on the city as a whole, especially in social terms (the Health and Safety of the Inhabitants). Furthermore, it also underlines that for the improvement of the economic conditions of this area, more integrated actions and policies are required.
Moreover, to analyze in depth the possible impacts of the considered scenarios in terms of changing urban resilience performance over time, final indexes referring to time 15 (t15) have been calculated for both scenarios. These indexes are calculated using Equation (3) that combines the normalized performance at time 15 of these scenarios obtained by the relevant simulations (Appendix C) with the weights of each criterion, which is obtained through the ANP (Table 5). Figure 15 illustrates and compares the obtained indexes. The final urban resilience index for the strategic scenario is 0.85. At the same time, the index obtained by the inertial scenario is 0.10. These indexes underline that the strategic scenario could improve the resilience performance of the area over time, including through a reduction of the soil consumption, the decontamination of brownfields, an improvement of the level of the safety of inhabitants and the adaptation of many buildings exposed to risks.

Conclusions
This paper illustrates the application of an integrated approach based on the combination of SDM and ANP to simulate the effects on urban resilience performance over time of two different scenarios. The main challenge of this application was trying to provide an evaluation tool that is able to assess urban resilience within its complexity, multidimensionality and dynamic behavior over time [18,37] in response to the current lack of an urban resilience measurement tool [94]. The choice of combining these two techniques arises from their common aspects, which can be summarized as follow: (1) They both examine the nature of the decision problems under analysis. Both SDM and ANP are applied to describe (SDM) and analyze (ANP) complex and dynamic decision problems in many different fields, such as urban transformation.
(2) They both assess the typology of mutual relationships. In real-world decision making, some problems have to be investigated considering the non-linear and feedback relationships between the variables involved. Both SDM and ANP are used based on this principle and for this reason they have been integrated to create a unique framework to analyze complex decision processes; (3) They both assess data in similar ways. Both SDM and ANP consider both qualitative and quantitative data. This paper allows us to underline both the strengths and the weaknesses of the proposed integrated approach. The main strengths can be summarized as follow:  The causal loop diagram can be considered as a suitable supporting tool in finding the relationships between the criteria to include in the ANP network. Experts were able to determine how the system functioned during the causal loop development (Figure 9).
 The integrated approach based on the combination of SDM and ANP can be considered suitable to respond to the need of an evaluation tool that is able to consider urban resilience based on its complexity, multidimensionality and dynamical behavior over time. In fact, it was possible to analyze the urban system as a complex and adaptive system, in which interventions on a specific urban area can elicit impacts (either positive or negative) on a city in a wide perspective. However, the main limitation of the proposed approach is the lack of the spatial scale both in the simulation and evaluation of the impacts of the strategies. Thus, an interesting future implementation approach involves the integration of this combined approach with GIS, to also assess the effects of urban development strategies in the spatial scale [45,95]. The final aim is providing an integrated tool able to show both the evolution of the impacts of strategies in a unique map and to support decision-makers in defining strategies based on their effects both on a spatial and a temporal scale.
Author Contributions: All authors equally contributed to the article. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.