Urban Resilience and Climate Change: Developing a Multidimensional Index to Adapt against Climate Change in the Iranian Capital City of Tehran

: Urban resilience studies have increased during recent years due to the signiﬁcance of climate change as an alarming issue in centralized and highly populated cities where urban functionalities are disintegrated. Towards this aim, an Exploratory Factor Analysis (EFA) was enrolled to streamline the urban resilience to climate change over the 22 districts in Tehran after assessing the resilience objectivity. Based on the results, the city coverage was classiﬁed into best (41%), moderate (15%), low (14%), and least resilient (30%). In addition, the urban municipal districts were classiﬁed into ﬁve functional zones including Wellbeing-wealth (WWZ), Ecological Conservation (ECZ), Core (CZ), Downtown (DZ), and Neutral Zone (NZ) after evaluating the concept of urban functionality in the resilience framework. The results indicated that the socio-cultural component is considered as the fundamental necessity, along with eco-environmental and economic components in capacity building to urban climate resilience. In fact, more than half of the Tehran coverage is regarded as resilient. Thus, the rest should be prioritized, despite the need to inspire from top-ranked districts, especially D4. However, downtown and neutral zones, especially D9 and D21, which account for up to 12% of the least resilient areas, should be evaluated seriously. Finally, the robustness of the proposed methodology was compared to the studies addressing the same concept, and we offer some preparatory and adaptive measures for urban planners and policymakers.


Introduction
In general, cities can hardly be pushed to be denser, because half of the world's population lives in urban areas.Therefore, the space left is gradually nibbled regardless of environmental concerns, leading to grounds with more impervious surfaces and less manageable situations.In addition, the interdependence of functioning units in cities plays a significant role in their resilience against changing environments [1].A broad consensus is observed for making the cities more resilient and bundling sustainable development efforts with urban resilience, although no agreement has ever covered the entire concept related to resilience among scientists.Different resilient patterns are reported within or across cities, resulting in increasing the costs of implementing resilient strategies due to the inequalities which lead to unevenly distributed urban services [2].
Urban resilience has attracted a lot of attention within the studies focused on cities, climate change, and sustainable development.Cities should be able to bounce back from disturbances due to the risks related to natural or anthropogenic incidences.Policy and Urban Sci.2023, 7, 7 2 of 18 decision makers always consider climate-driven events such as flooding and droughts seriously due to the level of damage and casualties, as well as their shocks and stresses [3].Climate and meteorological events occur less often compared to the past, and the return period is regarded just as a figure, not necessarily indicating their occurrence.Thus, the individual strategies should be transformed into curriculum-oriented ones derived from a robust climate framework in order to move towards sustainability [4].
Based on the Asian Cities Climate Change Resilience Network (ACCCRN) definition, climate resiliency in urban areas should be operationalized by a sweeping framework considering theoretical and empirical knowledge.Such a framework integrates urban systems, people, and institutions in a governable method, without which individual efforts cannot contribute to adaptation significantly.Therefore, people become capable through an improved urban system and better accept the policies and measures made by institutions, resulting in building a resilient city under the exposure of climate change.Although the aforementioned framework is inspired in all manners, resilience against natural hazards generated by climate change cannot be created easily as a better representative for disasters [5].The adaptation and mitigation strategies designed in compliance with the aforementioned framework guarantee the best outcomes while considering the possible drawbacks on other sectors, especially in arid and semi-arid areas.Thus, the urban systems with more protection and preparedness (resistance), less eco-hydrologic disintegration (absorption), and sustaining economy (recovery) are regarded as intrinsically more resilient against climate extremes [2,6].
Tehran, which is located on the south side of the Alborz Mountains, has exhibited a decrease in annual average rainfall and an increase in temperatures on average during the past decades [7].The Iranian capital of Tehran was selected as a case study to examine the resilience capacity in its 22 municipalities to climate change, especially drought, flooding, and scorching temperatures, considering the continuous environmental pressures, demographic growth, and urban complexity.The Climate Resilience Index (CRI), as a multidimensional concept, remains a great challenge amid the studies conducted.Accordingly, contributing to creating management instruments to be useful in decision-making processes resulted in developing a composite CRI bearing several dimensions of the urban systems, including Socio-cultural Resilience (SCR), Economic Resilience (ER), Inst-infrastructural Resilience (IIR), and Eco-environmental Resilience (EER) dimensions/components.
Based on the Organization for Economic Co-operation and Development (OECD) [8], the types of obstacles, nature of data, and objectives are among the parameters to be included while constructing a multidimensional index.A large number of researchers proposed an index for multi-hazards [9][10][11][12], along with the current resilience metrics designed to deal with a particular hazard [6,[13][14][15].Filling the knowledge gap and better understanding urban resilience have been addressed in the literature through engaging different spatial scales and varied approaches.Batica [16] investigated urban resilience at different scales including block cells, districts, and a whole city using the FAST approach over the case studies of Homburg, Beijing, and Barcelona.Ordination and Spatial Multi-Criteria Decision-Making (SMCDM) methods were recognized to be popular in the literatures reviewed.For example, Kotzee and Reyers [17] analyzed PCA to weight underpinning factors in order to reach a flexible multidimensional index.The spatial distribution of flood resilience over the landscape was mapped corresponding to the most contributing factors including social, ecologic, and economic components.
In addition, Ghaedamini Harouni et al. [18] proposed a multidimensional health index for Tehran utilizing factor analysis followed by a sensitivity assessment.Based on the results, the most significant health factors in Tehran included mental, social, and physical health status, respectively.Furthermore, Movahhed et al. [19] studied poverty based on the concept of political economy in the capital city of Tehran.Combining MCDM (VIKOR) and orientation methods indicated that housing, economic, social, and educational poverty are among the most influential components affecting citizens in the city.Furthermore, Zheng et al. [20] reviewed urban resilience to climate change in Beijing, applying an EFA as an Urban Sci.2023, 7, 7 3 of 18 objective scoring method, and iterpreted the derived loading of urban areas (municipal districts) on each factor according to their origins.The researchers adopted varied MCDM approaches to deal with urban resilience, while the AHP was considered the common method to adjust the indicator weights.Disaster vulnerability [21], earthquake preparedness [22], and urban physical resilience [23] were analyzed using the AHP, FAHP, and ANP approaches.In another study, Zhang et al. [24] analyzed the conceptual frameworks as a perfect example related to engaging time-series in urban resilience and assessed the Chinese cities of Beijing, Chongqing, and Yiwu over a period of six years.Chen and Leandro [15] evaluated time-series indicators for event and recovery phases, respectively.It is worth noting that the share of driving factors was regarded as disparate among the aforementioned areas.
In addition, Satour et al. [6] examined flood resilience in three coastal areas placed in the Moroccan semi-arid region considering social, physical, economic, and natural dimensions.The geographical information system (GIS) was utilized uniquely since mapping was performed in finer-resolution cells than that of the urban districts.Furthermore, Brunetta and Salata [25] investigated the fine-scale resilience mapping of Turin, Italy by engaging a holistic approach to detect vulnerable areas and generated state and pressure variables including shocks and disturbances according to ecosystem service features.Furthermore, Zheng et al. [20] studied urban resilience in the Chinese city of Panzhihua as a representative of resource base cities.A nearly two-decade assessments revealed that the city's resilience including health and well-being, infrastructure and environment, economy and society, as well as government capacity increase and a move toward stability, despite some momentary fluctuations.In another study, Chen and Quan [26] reviewed urban resilience to the Covid-19 pandemic applying TOPSIS and AHP and reported that the spatial distribution of the resilience across the city faces challenges including agglomerating clusters of high-resilience areas versus dispersed moderate-to low-resilience ones, despite an ever-growing increase in resilience toward the pandemic.The aforementioned researchers analyzed the significance of governance during pandemics.Moghadas et al. [27] prioritized municipal areas in Tehran from the perspective of flash flood resilience using TOPSIS.The hybrid AHP-TOPSIS model ranked the resilience components based on their contribution in a descending order including critical infrastructures, environment, economy, social, institutional, and community resilience.
The dimensions of urban resilience to change climate have not been well understood, mostly due to uncertainty in data-scarce areas, although a large number of studies have been conducted in this field [28].As indicated, the urban resilience across the Tehran municipalities was assessed utilizing the EFA compared to the hybrid MCDM method integrating ANP and DEMATEL, as well as enrolling in the recent submission [29].This study therefore aims to develop a multidimensional index to examine the resilience of cities to climate change.Preparing a list of different indicators related to the subject and combining them in a logical way is considered in this way.The Tehran metropolis has been chosen to test this new index.

Materials and Methods
In the present study, the resilience to climate change in Tehran was evaluated basically considering the capabilities and limitations in its 22 municipal districts.The proposed composite index aims to bear multidimensionality under the four resilience macro criteria including Socio-cultural, Economic, Inst-infrastructural, and Eco-environmental dimensions, as well.Accordingly, a couple of recommended resilience metrics implemented by researchers were selected, and the final list was formed in compliance with data availability and robustness.Then, a correlation analysis was performed to examine whether there is any sizable and explicable dependence between manifest variables (resilience indicators) before prioritizing urban districts followed by extracting resilience factors (construct variables).Finally, the statistical contribution of sampling units (urban districts) in the final resilience score was interpreted based on their respective functional zone.The degree to which an urban district belongs to a given functional zone is regarded as an abstract measure considering human development criteria and marginal indices defined by official reports such as land use, commercial activities, taxes paid, and the like.Therefore, there are five different functional zones in Tehran including Wellbeing-welfare (WWZ), Ecological Conservation (ECZ), Core (CZ), Downtown (DTZ), and Neutral Zone (NZ), as well [30].

Conceptual Framework
The conceptual framework in the study is based on reducing the risk of climate change by making the cities more resilient [2,6].Towards this aim, the countertendency changed to interdependency, and the realistic framework was drawn (Figure 1).
whether there is any sizable and explicable dependence between manifest variables (resilience indicators) before prioritizing urban districts followed by extracting resilience factors (construct variables).Finally, the statistical contribution of sampling units (urban districts) in the final resilience score was interpreted based on their respective functional zone.The degree to which an urban district belongs to a given functional zone is regarded as an abstract measure considering human development criteria and marginal indices defined by official reports such as land use, commercial activities, taxes paid, and the like.Therefore, there are five different functional zones in Tehran including Wellbeing-welfare (WWZ), Ecological Conservation (ECZ), Core (CZ), Downtown (DTZ), and Neutral Zone (NZ), as well [30].

Conceptual Framework
The conceptual framework in the study is based on reducing the risk of climate change by making the cities more resilient [2,6].Towards this aim, the countertendency changed to interdependency, and the realistic framework was drawn (Figure 1).

Study Area
Tehran is the capital of Iran, and Tehran province, with a population of around 8.7 million, is host to 22 districts (Figure 2) and has been a destination for mass migrations from all over the country since the 20th century.In addition, Tehran is considered the most populous city in Iran and Western Asia, with the second largest metropolitan area in the Middle East after Cairo, and is ranked 24th in the world by the population of its metropolitan area [31].Furthermore, Tehran is regarded as a mountainside city spreading at the south domain of the Alborz Mountains with an altitude of 900-1700 meters above sea level and is exposed to earthquakes, floods, droughts, storms, and surrounding watery watersheds.The city has not experienced much flooding, despite being prone to such hazards.Floods usually occur in spring after heavy rainfalls mostly through water-logged streets, and earthquakes remain the most harmful threat.It is noteworthy that industrialtraffic-induced air pollution is considered the most recurrent hazard in the city.Meanwhile, informal settlements, the homeless, especially in winter, and drug abuse are left as common social concerns [32].

Study Area
Tehran is the capital of Iran, and Tehran province, with a population of around 8.7 million, is host to 22 districts (Figure 2) and has been a destination for mass migrations from all over the country since the 20th century.In addition, Tehran is considered the most populous city in Iran and Western Asia, with the second largest metropolitan area in the Middle East after Cairo, and is ranked 24th in the world by the population of its metropolitan area [31].Furthermore, Tehran is regarded as a mountainside city spreading at the south domain of the Alborz Mountains with an altitude of 900-1700 meters above sea level and is exposed to earthquakes, floods, droughts, storms, and surrounding watery watersheds.The city has not experienced much flooding, despite being prone to such hazards.Floods usually occur in spring after heavy rainfalls mostly through water-logged streets, and earthquakes remain the most harmful threat.It is noteworthy that industrialtraffic-induced air pollution is considered the most recurrent hazard in the city.Meanwhile, informal settlements, the homeless, especially in winter, and drug abuse are left as common social concerns [32].

Data in Use, Variable Selection
In the present study, the data were mostly collected from the Statistical Center (SCI; amar.org.ir),Municipality (MOF; tehran.ir),and Provincial Government Office of Tehran (PGOT; ostan-th.ir)through exploring annals and atlases related to the capital city of Tehran.However, a large amount of information was obtained by reviewing the literature related to Tehran from different perspectives including green space [33], life expectancy [34], extreme weather parameters [7], spatial inequality [35], earthquake preparedness [22], physical resilience [23], and poverty [19].Table 1 indicates the complete list of indicators to form the multidimensional index.In addition, a descriptive and analytical view of the datasets used in the study is reflected in the following sections (Figure 3 and Table 1).

Data in Use, Variable Selection
In the present study, the data were mostly collected from the Statistical Center (SCI; amar.org.ir),Municipality (MOF; tehran.ir),and Provincial Government Office of Tehran (PGOT; ostan-th.ir)through exploring annals and atlases related to the capital city of Tehran.However, a large amount of information was obtained by reviewing the literature related to Tehran from different perspectives including green space [33], life expectancy [34], extreme weather parameters [7], spatial inequality [35], earthquake preparedness [22], physical resilience [23], and poverty [19].Table 1 indicates the complete list of indicators to form the multidimensional index.In addition, a descriptive and analytical view of the datasets used in the study is reflected in the following sections (Figure 3 and Table 1).

Data in Use, Variable Selection
In the present study, the data were mostly collected from the Statistical Center (SCI; amar.org.ir),Municipality (MOF; tehran.ir),and Provincial Government Office of Tehran (PGOT; ostan-th.ir)through exploring annals and atlases related to the capital city of Tehran.However, a large amount of information was obtained by reviewing the literature related to Tehran from different perspectives including green space [33], life expectancy [34], extreme weather parameters [7], spatial inequality [35], earthquake preparedness [22], physical resilience [23], and poverty [19].Table 1 indicates the complete list of indicators to form the multidimensional index.In addition, a descriptive and analytical view of the datasets used in the study is reflected in the following sections (Figure 3 and Table 1).The multidimensional index of urban resilience to climate change was formed to provide urban managers and policymakers with some generic knowledge and solid understandings, resulting in creating practical adaptation strategies.A large number of expert meetings were held to screen the most appropriate and applicable criteria to the city of Tehran in addition to the aforementioned studies [29].The architecture of the resilience framework was based upon the factor analysis as a more tangible mechanism, unlike previous studies which focused on expert-based weighting.In fact, the disparity or similarity amongst the indicators was determined by extracting common factors which define the loading of each manifest indicator and the contribution of sampling units in creating the construct factors, which in turn avoids conducting subjective resilience weighting all the way through generating an objective scoring process.It is worth noting that CRI for each urban district was retrieved by summing the amount of variance explained multiplied by the given factor loading.Accordingly, the high-speed factor covering the most explained variance during the rotation process plays a significant role in forming the significance and direction of the resilient indicators.The aforementioned factor can be explicated by reflecting the share of the main resilience dimensions such as SCR, ER, IIR, and EER on the extracted factors, as well.

Multivariate Analysis
Factor analysis outperforms its multivariate rivals such as PCA, CA, and the like by limiting the data boundary issues such as normal distribution supposition, giving the researchers the ability to investigate whether the data are in compliance with their null structure (confirmatory) or to discover multidimensionality inside the data, as well as providing novel extraction techniques, such as a rotated component matrix.According to Moghadas et al. (2019), the factor analysis reduces data dimensionality and standard errors reflected on component loadings simultaneously.The statistical structure describing the factor analysis conducted during the study is as follows [20].
where, X i is n original indicators, f i (j = 1, 2, 3, . . ., k) represents the k latent factors, and e i is regarded as the difference of the ith indicator factor.In addition, ∝_ij indicates the loading coefficient of the ith indicator on the jth latent factor (f i ), reflecting the correlation degree between the original indicator and respective common factor.The lower ∝_ij decreases the weight of X i in the common factor (f i ).The routine analytical process is described as follows.

Data Pre-Processing
The statistics, including descriptive, correlation, and confirmatory factor analysis, as well as respective fitness tests, were applied using SPSS version 16 (released 2007).The resilience indicators, which are usually normalized based on Equations ( 2) and (3) for positive and negative contributions, respectively, were not modified while obtaining the CRI regarding leaving the factor and sample loads intact.
where i represents the ith sample (i = 1, 2, 3, . . ., n), j is considered as the jth indicator (j = 1, 2, 3, . . ., m), Max indicates the maximum value of a certain indicator, and Min is regarded as the minimum value of the same indicator.

Conducting the Factor Analysis
The factor analysis populates the construct factors based on the contribution of manifest indicators.In this regard, the variables exceeding the loading brink of 0.50 can be attributed to the construct factors.Accordingly, a Kaiser-Meyer-Olkin (KMO) test and Bartlett's test of sphericity were implemented to find whether the data in use are considered as adequate and co-dependent.A Principal Component Analysis (PCA) powered by the rotation method of Varimax and Kaiser Normalization was performed after the factor analysis to extract the main factors.In addition, the objective technique of weighting was based on the rotation sums of the squared loading.

Visualization
All of the spatial operations including visualizing resilient indicators and CRI classes, as well as ranking, were conducted enrolling ArcGIS (Ver.10.6).

Results
A multidimensional metric populated by the manifest indicators was designed following the conceptual framework proposed in the study to review urban resilience to climate change in the capital city of Tehran.

Factor Analysis
To conduct the intended exploratory factor analysis, the assumptions related to dependent observations and normal distribution fully complied, meaning that almost all of the data had a normal distribution save for I17 (Infrastructure vulnerability) and I21 (Slope).Then, the Kaiser-Meyer-Olkin (KMO) test and Bartlett's test of sphericity were implemented to find whether the data in use are regarded as adequate and co-dependent.Finally, the EFA analysis was performed after complying with the interdependence (rejecting the Bartlett null hypothesis) and adequacy results (0.75 < KMO = 0.89 < 1) (Table 3).The construct factors followed by the percent of the variance explained were retrieved after forming the EFA.In this respect, the most contributing manifest indicators can be observed loading the respective factor(s) in bold (Table 4).To illustrate a more vivid picture of the conducted EFA, a 3D component plot and the causal network were drawn, as well (Figures 4 and 5).Based on the results, 16 out of 22 resilience indicators, including I22-f1, I21-f1, I13-f1, I4-f1, I6-f1, I9-f1, I18-f1, I11-f1, I17-f2, I15-f2, I14-f3, I20-f3, I3-f4, I2-f4, I16-f5, and I7-f6 were directly loaded on their most corresponding construct factor, and the rest, including Air quality index (I19), Commercial land use (I8), Death rate (I5), Employment (I12), Public awareness (I1), and Urban worn-out areas (I10) could not heavily contribute to any factors.

Sample Loadings (Districts)
To obtain the CRI, urban district loadings including shares of sampling units were calculated and summarized (Table A2-Appendix A and Figure 6).Cross-investigating the contribution of urban districts based on their respective resilience class and functional zone is considered as a step toward providing an analytical view regarding the amount of resilience in Tehran.The resilience classes including best, moderate, low, and least were mainly populated based on the CRI range observed.In addition, the functional zones including WWZ, ECZ, CZ, DTZ, and NZ represented commerce, wealth, residential, and

Sample Loadings (Districts)
To obtain the CRI, urban district loadings including shares of sampling units were calculated and summarized (Table A2-Appendix A and Figure 6).Cross-investigating the contribution of urban districts based on their respective resilience class and functional zone is considered as a step toward providing an analytical view regarding the amount of resilience in Tehran.The resilience classes including best, moderate, low, and least were mainly populated based on the CRI range observed.In addition, the functional zones including WWZ, ECZ, CZ, DTZ, and NZ represented commerce, wealth, residential, and infrastructures across the city.

Sample Loadings (Districts)
To obtain the CRI, urban district loadings including shares of sampling units were calculated and summarized (Table A2-Appendix A and Figure 6).Cross-investigating the contribution of urban districts based on their respective resilience class and functional zone is considered as a step toward providing an analytical view regarding the amount of resilience in Tehran.The resilience classes including best, moderate, low, and least were mainly populated based on the CRI range observed.In addition, the functional zones including WWZ, ECZ, CZ, DTZ, and NZ represented commerce, wealth, residential, and infrastructures across the city.

Mapping CRI
All of the municipal districts in Tehran were ranked from top to bottom and classified into the intended resilience classes from most to least resilient (Figure 7) after visualizing resilience indicators based on their entity (Figure 8).

Mapping CRI
All of the municipal districts in Tehran were ranked from top to bottom and classified into the intended resilience classes from most to least resilient (Figure 7) after visualizing resilience indicators based on their entity (Figure 8).

Discussion
Urban resilience is considered as a complex and challenging issue, though measurable, which can be addressed subjectively, unlike most studies conducted in this field.The studies which internalize the resilience theory by their own language are increasingly assorted.In this regard, the paradigm of scale, dimensionality, and timeline is regarded as a controversial matter among the scientific community.The present study aims to analyze the current knowledge and understandings by assessing Tehran's resilience over its 22 municipal districts.Conducting a single-hazard study plays a significant role in focusing on urban resilience to climate change and measuring the levels of preparedness as long as the hazards turn up.Accordingly, an EFA was enrolled based on PCA and Varimax as factor extraction and rotation methods.
The CRI proved to be a robust and informative framework due to its flexibility, integration, and applicability, although the multidimensional index (CRI) presented is not considered as the perfect instrument in characterizing climate resilience in urban areas.Based on the aforementioned studies, a list of resilience indicators was shortlisted based on data availability and veracity.Towards this aim, the urban system was defined having four resilience components including SCR, ER, IIR, and EER, and the 22 manifest indicators were assigned.Then, pre-processing was conducted to examine whether the interdependence and the adequacy tests complied well (KMO = 0.89, Bartlett sig .00).
The datasets utilized to characterize the resilience of Tehran, as the capital and the most populous city in Iran, demonstrated a normal pattern related to distributing all of the resilience indicators except for critical infrastructure vulnerability (I17) and slope (I21).Few outliers were observed on SCR and EER components after evaluating the data descriptively.In fact, the ER was regarded as a nearly zero-content outlier component.Based on the coefficient of variation (CV), life expectancy (I6) was considered as the most centralized variable, while access to urban services (I16) and infrastructure vulnerability (I17) were regarded as the most decentralized ones.Some remarkable relationships were detected after investigating the resilience indicators in paired format.For example, the positive link between migration and life expectancy is considered as a direct matter related to wealth and welfare, as the privileged tend to reside from other areas than the WWZ zone such as D1, D2, D3, and D5.Similarly, the welfare-elevation relationship illustrates that the privileged are genuinely interested in uptown with higher elevations.High classes of slope improve the quality of the residents'

Discussion
Urban resilience is considered as a complex and challenging issue, though measurable, which can be addressed subjectively, unlike most studies conducted in this field.The studies which internalize the resilience theory by their own language are increasingly assorted.In this regard, the paradigm of scale, dimensionality, and timeline is regarded as a controversial matter among the scientific community.The present study aims to analyze the current knowledge and understandings by assessing Tehran's resilience over its 22 municipal districts.Conducting a single-hazard study plays a significant role in focusing on urban resilience to climate change and measuring the levels of preparedness as long as the hazards turn up.Accordingly, an EFA was enrolled based on PCA and Varimax as factor extraction and rotation methods.
The CRI proved to be a robust and informative framework due to its flexibility, integration, and applicability, although the multidimensional index (CRI) presented is not considered as the perfect instrument in characterizing climate resilience in urban areas.Based on the aforementioned studies, a list of resilience indicators was shortlisted based on data availability and veracity.Towards this aim, the urban system was defined having four resilience components including SCR, ER, IIR, and EER, and the 22 manifest indicators were assigned.Then, pre-processing was conducted to examine whether the inter-dependence and the adequacy tests complied well (KMO = 0.89, Bartlett sig 0.00).
The datasets utilized to characterize the resilience of Tehran, as the capital and the most populous city in Iran, demonstrated a normal pattern related to distributing all of the resilience indicators except for critical infrastructure vulnerability (I17) and slope (I21).Few outliers were observed on SCR and EER components after evaluating the data descriptively.In fact, the ER was regarded as a nearly zero-content outlier component.Based on the coefficient of variation (CV), life expectancy (I6) was considered as the most centralized variable, while access to urban services (I16) and infrastructure vulnerability (I17) were regarded as the most decentralized ones.
Urban Sci.2023, 7, 7 13 of 18 Some remarkable relationships were detected after investigating the resilience indicators in paired format.For example, the positive link between migration and life expectancy is considered as a direct matter related to wealth and welfare, as the privileged tend to reside from other areas than the WWZ zone such as D1, D2, D3, and D5.Similarly, the welfare-elevation relationship illustrates that the privileged are genuinely interested in uptown with higher elevations.High classes of slope improve the quality of the residents' water.Ecologically, dwellers might actually benefit from extremely better water and air quality in higher elevations in uptown.The mixture of slope/elevation addresses the issue likewise.
The direct positive correlation between accident insurance and life expectancy indicates the resilience perfectly once hazards are guaranteed, at least financially to those with better living standards.In addition, the positive correlation between urban worn-out areas and commercial land use represents low resilience, because the commercial zone is placed in hazard-prone areas.Population density is regarded as another significant issue making the situation worse during the crisis.The positive correlation between vulnerable infrastructure and commercial land use proves the necessity to implement adaptive measurements and retrofit such areas.Urban planners and policymakers should focus on the concurrence of the poverty line and urban worn-out areas, since expecting some climate-oriented hazards and being in a basket case lead to bad consequences.Further, the negative relationship between life expectancy and the poverty line can be apprehended intuitively.Therefore, the lack of social and economic advantages reduces the preparedness to hazards.
Reviewing the resilience indicators independently reveals interesting results, since such indicators which positively/negatively impact the construct factors hardly contradict the studies conducted so far.For instance, references [6,19,21,27] found that migration (I4) contributed to urban resilience negatively, although the multivariate analysis displayed a positive contribution of migrants based on available data.The migration of the rich was reported to be the case regarding the dataset applied, although there are possible sources of uncertainty which lead to systematic errors.The employment rate (I12) was regarded as another contradictory issue detected during investigations, which populated the construct factors negatively according to the statistical analysis, despite being classified under the positive contributors.To justify the phenomenon, low-to moderate-income households are supposed to embark on jobs even under the super-competitive economic and job market situation ongoing in the city.Accordingly, the reverse contribution of the indicator can only be justified by the rich relying on their deposits and less applying to the jobs based on official surveys conducted by the CSI.Other studies indicated an adverse contribution compared to the present one [19,21,24,27,36].Access to urban services (I16) was considered as another indicator which was classified as a positive contributor, despite loading the construct factors negatively.The indicator was no longer regarded, since its statistical significance could be neglected completely.The slope appears to be a negative resilience contributor [22], while data indicate a direct correlation between the aforementioned indicator and other strong contributors.Finally, the EFA classified the 22 indicators into the two groups including I2, I3, I5, I8, I9, I10, I12, I16, I17, I19 and I1, I4, I6, I7, I11, I13, I14, I15, I18, I20, I21, I22 with negative and positive directions, respectively, based on the direction of communication in urban resilience.
Based on the rotation in squared sums of loadings, the extracted factors accounted for up to 85% of the total variance.It is noteworthy that f1 and f2, as the first two factors, accounted for nearly 50% of variance combined.As Table 4 shows, the factors extracted appear in a decent order.However, the notable distinction between the Extraction and Rotation sums of squared loadings is observed in the values of the first factor (f1) where the Rotation falls behind the Extraction method by an approximately 6% margin.Furthermore, the EFA was analyzed by dividing the results into the following phases of variable and sample loadings.

Variable Loadings
Based on the Rotated Component Matrix, the strong communication of 16 out of 22 resilience indicators demonstrated that SCR played an active role in forming f1, f4, and f6.Inst-infrastructural resilience (f2, f5) and Eco-Environmental resilience components (f1, f3) were considered as the second active resilient components, while the Economic resilience (ER) only appeared in factor f1.In other words, f1 is influenced by SCR, ER, and EER, while f2 and f5 are affected by IIR.In addition, f3 is influenced by IIR and ER, while f4 and f6 are affected by the SCR.
The main outcomes of the resilience-concerned studies mostly designed to bear natural hazards such as earthquakes, flooding, and climate change are discussed as follows.Kotzee and Reyers [17] prioritized flood-resilience-contributing dimensions as social > ecologic > economic components, which is surprisingly consistent with climate resilience findings (f1: SCR > EER > ER) (Table 4).The observed compliance stems from the similar nature of the hazards, including climate change and floods.
However, Zhang et al. [24] reached the priority of Inst-infrastructural > eco-environmental > economy > social components in preparing for climate resilience, while Moghadas et al. [27] prioritized flash flood resilience components, the result of which is in line with that of the previous studies (critical infrastructures > environment > economy > social > institutional and community resilience).Zheng et al. [20] found that economic resilience is considered as far more significant compared to other resilience dimensions including infrastructural, institutional, managerial, and environmental.In addition, Movahhed et al. [19] reported that housing, economy, as well as social and cultural sectors were the most contributing dimensions in the resilience to poverty, respectively.Furthermore, references [6,26] used equal weights in combining resilience indicators, while others utilized at least one objective or subjective weighting method based on statistical analysis and panel decisions, respectively, to prioritize the resilience components.

Sample Loadings
The CRI ranking followed by resilience classes in Tehran (Figure 8) demonstrated that the most resilient districts are located in northern Tehran, mainly surrounded by WWZ and EC (D4, D1, D5, D2, and D3).D4 was recognized as the top-ranked district mostly due to its ecological status and resulting functionality across the city.District 22 failed to pass the first-level standards required, despite being situated in the ECZ.Being situated within the CZ and ECZ helped moderate-resilience-level districts (D15, D12, D14, D6, and D11) to partly resist, absorb, and adapt to the climate-oriented crises.
However, the low-resilience urban districts including D8, D7, D16, D20, and D13 are densely populated, while they benefit from accessing services due to their mixed situation in CZ and DTZ.The least resilient class, including D22, D18, D17, D19, D21, and D9, is regarded as a combination of NZ and DTZ, in which the districts are not regarded to be rich enough to bounce back fast and not privileged with infrastructures and services available throughout the city.It is worth noting that the corresponding districts have been less included in land use planning as far as districts D9 and D21 are isolated in part.Analyzing the situation showed that district 22 suffers from a lack of services and infrastructures, despite benefiting from less population density, well-structured highways, and enough green space.
Based on Table A2-Appendix A, more than half of the city coverage is placed in the best/moderate resilience classes, and urban planners and policymakers should focus on the rest of the city.The low-resilience areas surrounded by DTZ and NZ should benefit from adaptive measures first.Comparing the previous study with the present one indicates some compliance regarding the unquestioned role of environment and social dimensions.However, some inconsistencies are observed in the objectivity of the approach (ANP vs. EFA) and the data used.Based on the results, both studies show consistency on the most and least resilient areas, and inconsistency is observed in the middle classes.In addition, district 22 was considered as the most challenging area in assigning a resilience class.

Conclusions
Urban resilience to climate change increased decisively due to climate-oriented hazards in a large number of centralized and highly populated cities where urban functional units are not well-distributed across the cities' domain.The present study seeks to enhance the understanding of climate resilience in the metropolitan area of Tehran in the absence of enrolling a sweeping framework.Thus, a tangible approach equipped with GIS was adopted to outline and illustrate the resilience hotspots over the 22 municipal districts in Tehran.Interpreting the EFA results coupled with CRI classification and functional zoning gave more insight than applying a simple subjective MCDM method.The repetitiousness of the proposed multidimensional metric was regarded as the other concern following its applicability in varied scales.Further studies can be conducted by engaging other resilience components such as physical and ecosystem-based service dimensions due to the high flexibility observed in the framework of the study.However, a finer scale and time-series should be utilized where applicable.
Based on the results, socio-cultural dimension including eco-environmental and economic components with other degrees of significance should be considered during capacity building for urban climate resilience.More than half of the Tehran coverage is observed in resilient areas, among which D4 should be further studied, because the literature is unanimously consistent as to the foremost ranking among other urban districts in Tehran.However, nearly 14 and 30% of the total areas in Tehran were classified as low-and least-resilience classes, respectively, indicating below-zero resilience scores.Based on the assessments conducted on corresponding functional zones, districts 9 and 21 should be involved in immediate and futuristic development planning.In other words, the existing neutral functional zone covering almost 12% of the total urban areas should not be disregarded as usual.
The adaptive measures proposed for empowering the capabilities and enhancing resilience status include increasing public awareness towards building capacity among moderate-income families to the same extent as the resilient rich, planning to repair and retrofit critical urban infrastructures while granting incentives and tax exemptions to those in critical risk of climate hazards, especially low-income households, updating hazard headquarters' workforce and fleet, as well as land use and population engineering.Finally, the conformity and robustness of the proposed framework were compared with other studies, despite the affirmative compliance in most and least resilient areas.The middleresilience-ranged districts can only be focused on by finer-scale information.This study may have some potential limitations.The main limitation associated with these models is their reliance on a combination of criteria.Different criteria can lead to different results.Moreover, dealing with data uncertainties in such models is difficult due to large amounts of data.The most important limitation of this study was access to the necessary data.There is also a lack of real data to evaluate the accuracy of the model.

Figure 1 .
Figure 1.Conceptual framework applied in the study.

Figure 1 .
Figure 1.Conceptual framework applied in the study.

Figure 3 .
Figure 3.The standardized distribution of resilience indicators.* Plots in circle and star indicate outliers based on their divergence from mean.

Figure 3 .
Figure 3.The standardized distribution of resilience indicators.* Plots in circle and star indicate outliers based on their divergence from mean.

Figure 3 .
Figure 3.The standardized distribution of resilience indicators.* Plots in circle and star indicate outliers based on their divergence from mean.

Figure 4 .
Figure 4.The causal network of Tehran resilience to climate change derived from construct factors.

Figure 4 . 19 Figure 5 .
Figure 4.The causal network of Tehran resilience to climate change derived from construct factors.Urban Sci.2023, 7, x FOR PEER REVIEW 11 of 19

Figure 8 .
Figure 8.The visualized resilience indicators, in which red and green ramps account for negative and positive contributions to the resilience process, respectively.

Table 1 .
A list of dimensions, indicators, and their definition.) is defined as a function of solid wastes and wastewater per capita, which loads reversely on the supply chain during long-run crises.-Populationdensity(I 3 ); n per hectare Population density (I 3 ) is negatively connected with sustainability.Slope (I 21 ) and elevation (I 22 ) were selected as indicators related to stability, which contribute to urban resilience positively from the perspective of eco-environmental resilience.
+ Migration (I 4 ); n Migration (I 4 ) usually occurs among vulnerable and low-resilience populations.However, most migrants in Tehran are considered wealthy families willing to reside in the northern states of the city.-Deathrate(I 5 ); per 100,000Rising death rate (I 5 ) increases the intricacy of the life system in the city.+Lifeexpectancy(I 6 ); year Life expectancy (I 6 ) is explicitly related to the willingness to recover from crises.A hopeful community bears the aftermath of disasters better.+Healthoverall index (I 7 ); multidimensional A health overall index (I 7 ) reflects the total functionality of a city health system.Economic (ER)-Commercial land use (I 8 ); percent Commercial land use (I 8 ): The possibility of proximity between dangerous and safe uses is regarded as an appropriate indicator for urban resilience.Thus, commercial zones are considered as vulnerable due to the density and usual unsafe environments.[6,19,20,24,36,38] -Poverty line (I 9 ); multidimensional The poverty line (I 9 ) indicates the ability of a given community to incur and bounce back during backlashes.-Urban worn-out areas (I 10 ); n Urban worn-out texture (I 10 ) is regarded as another indicator contributing to resilience in line with other economic ones.Such areas pose the most potential threats from the perspective of energy consumption and residents' safety.+ Accident insurance (I 11 ); n The proportion of people covered by insurance (I 11 ) allows for recovering as fast as possible.+ Employment (I 12 ); % + Welfare (I 13 ); multidimensional Employment (I 12 ) and welfare (I 13 ) positively contribute to the readiness of the society for climate resilience.

Table 2 .
Descriptive statistics.As represented in Table A1-Appendix A, there are a large number of cases paired with a strong correlation greater than 0.7.However, the present study seeks to analyze those indicating solid grounds, in which the positive relationships between Migration

Table 3 .
Rotation in sums of squared loadings.

Table 4 .
Cont.Data in bold face indicate that their corresponding indicators have a common factor greater than 0.7.
Data in bold face indicate that their corresponding indicators have a common factor greater than 0.7.

Table A2 .
Municipal districts loadings in Tehran.
* Resilient classes include best, moderate, low, and least.